We discuss the statistical foundations of morphological star-galaxy separation. We show that many of the star-galaxy separation metrics in common use today (e.g. by SDSS or SExtractor) are closely related both to each other, and to the model odds ratio derived in a Bayesian framework by Sebok (1979). While the scaling of these algorithms with the noise properties of the sources varies, these differences do not strongly differentiate their performance. We construct a model of the performance of a star-galaxy separator in a realistic survey to understand the impact of observational signal-to-noise ratio (or equivalently, 5-sigma limiting depth) and seeing on classification performance. The model quantitatively demonstrates that, assuming realistic densities and angular sizes of stars and galaxies, 10% worse seeing can be compensated for by approximately 0.4 magnitudes deeper data to achieve the same star-galaxy classification performance. We discuss how to probabilistically combine multiple measurements, either of the same type (e.g., subsequent exposures), or differing types (e.g., multiple bandpasses), or differing methodologies (e.g., morphological and color-based classification). These methods are increasingly important for observations at faint magnitudes, where the rapidly rising number density of small galaxies makes star-galaxy classification a challenging problem. However, because of the significant role that the signal-to-noise ratio plays in resolving small galaxies, surveys with large-aperture telescopes, such as LSST, will continue to see improving star-galaxy separation as they push to these fainter magnitudes.

Recent work has shown that Milky Way-mass galaxies display an incredible range of stellar halo properties. Yet, the origin of this diversity is unclear. The nearby galaxy M81 $-$ currently interacting with M82 and NGC 3077 $-$ sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations in M81's halo. Using a unique three-filter strategy, and numerous HST calibration fields, we reveal M81's halo in never-before-seen detail. We resolve the halo to unprecedented $V$-band equivalent surface brightnesses of $>$34 mag arcsec$^{-2}$, and produce the first-ever global stellar mass density map for a Milky Way-mass stellar halo outside of the Local Group (LG). Using the minor axis, we confirm the previous assessment of M81 as one of the lowest mass and metal-poorest stellar halos known ($M_{\star} \sim 10^9 M_{\odot}$, [Fe/H] $\sim {-}1.2$) $-$ indicating a relatively quiet prior accretion history. Yet, in our global stellar mass census we find that tidally unbound material from M82 and NGC 3077 provides a substantial infusion of metal-rich material ($M_{\star} \simeq 6{\times}10^8$ $M_{\odot}$, [Fe/H] $\simeq {-}$0.9). We further show that, following the accretion of its massive satellite M82 (and the LMC-like NGC 3077), M81 will host one of the most massive and metal-rich stellar halos in the nearby universe. Thus, the saga of M81: following a relatively passive history, M81's merger with M82 will completely transform its halo from a low-mass, anemic halo rivaling the MW, to a metal-rich behemoth rivaled only by systems such as M31. This dramatic transformation indicates that the observed diversity in stellar halo properties is primarily driven by diversity in the largest mergers these galaxies have experienced.

We introduce AXS (Astronomy eXtensions for Spark), a scalable open-source astronomical data analysis framework built on Apache Spark, a widely used industry-standard engine for big data processing. Building on capabilities present in Spark, AXS aims to enable querying and analyzing almost arbitrarily large astronomical catalogs using familiar Python/AstroPy concepts, DataFrame APIs, and SQL statements. We achieve this by i) adding support to Spark for efficient on-line positional cross-matching and ii) supplying a Python library supporting commonly-used operations for astronomical data analysis. To support scalable cross-matching, we developed a variant of the ZONES algorithm (Gray et al. 2004) capable of operating in distributed, shared-nothing architecture. We couple this to a data partitioning scheme that enables fast catalog cross-matching and handles the data skew often present in deep all-sky data sets. The cross-match and other often-used functionalities are exposed to the end users through an easy-to-use Python API. We demonstrate AXS' technical and scientific performance on SDSS, ZTF, Gaia DR2, and AllWise catalogs. Using AXS we were able to perform on-the-fly cross-match of Gaia DR2 (1.8 billion rows) and AllWise (900 million rows) data sets in ~ 30 seconds. We discuss how cloud-ready distributed systems like AXS provide a natural way to enable comprehensive end-user analyses of large datasets such as LSST.

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion.

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber.

We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package.

The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithms running in these pipelines are individually sophisticated and interact in subtle ways. This paper provides an overview of those pipelines, focusing more on those interactions than the details of any individual algorithm.

The dwarf satellites of `giant' Milky Way (MW)-mass galaxies are our primary probes of low-mass dark matter halos. The number and velocities of the satellite galaxies of the MW and M31 initially puzzled galaxy formation theorists, but are now reproduced well by many models. Yet, are the MW's and M31's satellites representative? Were galaxy formation models `overfit'? These questions motivate deep searches for satellite galaxies outside the Local Group. We present a deep survey of the `classical' satellites ($M_{\star}$$\geqslant$4$\times$10$^5 M_{\odot}$) of the MW-mass galaxy M94 out to 150 kpc projected distance. We find $only\ two$ satellites, each with $M_{\star}{\sim}10^6 M_{\odot}$, compared with 6-12 such satellites in the four other MW-mass systems with comparable data (MW, M31, M81, M101). Using a `standard' prescription for occupying dark matter halos (taken from the fully hydrodynamical EAGLE simulation) with galaxies, we find that such a sparse satellite population occurs in $<$0.2% of MW-mass systems $-$ a $<$1% probability among a sample of five (known systems + M94). In order to produce an M94-like system more frequently we make satellite galaxy formation much more stochastic than is currently predicted by dramatically increasing the slope and scatter of the stellar mass-halo mass (SMHM) relation. Surprisingly, the SMHM relation must be altered even for halos masses up to 10$^{11}M_{\odot} -$ significantly above the mass scales predicted to have increased scatter from current hydrodynamical models. The sparse satellite population of this `lonely giant' thus advocates for an important modification to ideas of how the satellites around MW-mass galaxies form.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017).

Using the most recent prototypes, design, and as-built system information, we test and quantify the capability of the Large Synoptic Survey Telescope (LSST) to discover Potentially Hazardous Asteroids (PHAs) and Near-Earth Objects (NEOs). We empirically estimate an expected upper limit to the false detection rate in LSST image differencing, using measurements on DECam data and prototype LSST software and find it to be about $450$~deg$^{-2}$. We show that this rate is already tractable with current prototype of the LSST Moving Object Processing System (MOPS) by processing a 30-day simulation consistent with measured false detection rates. We proceed to evaluate the performance of the LSST baseline survey strategy for PHAs and NEOs using a high-fidelity simulated survey pointing history. We find that LSST alone, using its baseline survey strategy, will detect $66\%$ of the PHA and $61\%$ of the NEO population objects brighter than $H=22$, with the uncertainty in the estimate of $\pm5$ percentage points. By generating and examining variations on the baseline survey strategy, we show it is possible to further improve the discovery yields. In particular, we find that extending the LSST survey by two additional years and doubling the MOPS search window increases the completeness for PHAs to $86\%$ (including those discovered by contemporaneous surveys) without jeopardizing other LSST science goals ($77\%$ for NEOs). This equates to reducing the undiscovered population of PHAs by additional $26\%$ ($15\%$ for NEOs), relative to the baseline survey.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We present the discovery of d1005+68, a new faint dwarf galaxy in the M81 Group, using observations taken with the Subaru Hyper Suprime-Cam. d1005+68's color-magnitude diagram is consistent with a distance of $3.98_{-0.43}^{+0.39}$ Mpc, establishing group membership. We derive an absolute $V$-band magnitude, from stellar isochrone fitting, of $M_{V} = -7.94_{-0.50}^{+0.38}$, with a half-light radius of $r_{h} = 188_{-41}^{+39}$ pc. These place d1005+68 within the radius-luminosity locus of Local Group and M81 satellites and among the faintest confirmed satellites outside the Local Group. Assuming an age of 12 Gyr, d1005+68's red giant branch is best fit by an isochrone of [Fe/H] $= -1.90 \pm 0.24$. It has a projected separation from nearby M81 satellite BK5N of only 5 kpc. As this is well within BK5N's virial radius, we speculate that d1005+68 may be a satellite of BK5N. If confirmed, this would make d1005+68 one of the first detected satellites-of-a-satellite.

We present a deep imaging survey of the Virgo Cluster, designed to study the connection between cluster galaxies and Virgo's diffuse intracluster light (ICL). Our observations span roughly 16 square degrees and reach a 3-sigma depth of mu(B)=29.5 and mu(V)=28.5 mag/arcsec^2. At these depths, the limiting systematic uncertainties are astrophysical: scattered starlight from foreground Galactic dust, and variations in faint background sources. The dust-scattered starlight is well-traced by deep far-infrared imaging, making it distinguishable from true Virgo diffuse light. Our imaging maps the Virgo core around M87 and the adjacent M86/M84 region, in subcluster B around M49, and in the more distant W' cloud around NGC 4365. Most of the detected ICL is found in the Virgo core and within the W' cloud, with little evidence for extensive ICL in subcluster B. The large amount of diffuse light seen in the infalling W' cloud likely illustrates the importance of the group environment for generating ICL. The bulk of the detected ICL is fairly red (B-V=0.7-0.9), indicative of old stellar populations. We estimate a total Virgo ICL fraction of 7-15%, somewhat smaller than expected for massive, evolved clusters, suggesting that Virgo is still growing its ICL component. We trace M87's extremely boxy halo out to ~ 150 kpc, and show that the current stripping rate of low luminosity galaxies is insufficient to have built M87's outer halo over a Hubble time. Finally, we identify another large ultra-diffuse galaxy in Virgo, likely in the process of being shredded by the cluster tidal field.

The number, distribution, and properties of dwarf satellites are crucial probes of the physics of galaxy formation at low masses and the response of satellite galaxies to the tidal and gas dynamical effects of their more massive parent.To make progress, it is necessary to augment and solidify the census of dwarf satellites of galaxies outside the Local Group. M\"uller et al. (2015) presented 16 dwarf galaxy candidates near M83, but lacking reliable distances, it is unclear which candidates are M83 satellites. Using red giant branch stars from the HST/GHOSTS survey in conjunction with ground-based images from VLT/VIMOS, we confirm that one of the candidates, dw1335-29-- with a projected distance of 26 kpc from M83 and a distance modulus of $(m - M)_0 = 28.5^{+0.3}_{-0.1}$ -- is a satellite of M83. We estimate an absolute magnitude $M_V = -10.1 \pm{0.4}$, an ellipticity of $0.40^{+0.14}_{-0.22}$, a half light radius of $656^{+121}_{-170 }$ pc, and [Fe/H] = $-1.3^{+0.3}_{-0.4}$. Owing to dw1335-29's somewhat irregular shape and possible young stars, we classify this galaxy as a dwarf irregular or transition dwarf. This is curious, as with a projected distance of 26 kpc from M83, dw1335-29 is expected to lack recent star formation. Further study of M83's dwarf population will reveal if star formation in its satellites is commonplace (suggesting a lack of a hot gas envelope for M83 that would quench star formation) or rare (suggesting that dw1335-29 has a larger M83-centric distance, and is fortuitously projected to small radii).

We use extensive gravity-sensitive DDO 51 photometry over 5100 square degrees, combined with SDSS broadband photometry, to select a catalog of $\sim 4,000$ giant stars covering a large fraction of the high Galactic latitude sky and reaching out to $\sim 80$ kpc in the Galactic halo. This sample of bright and unbiased tracers enables us to measure the radial profile and 3D structure of the stellar halo to large distance which had previously only been measured with sparse tracers or small samples. Using population synthesis models to reproduce the observed giant star luminosity function, we find that the halo maintains a $r^{-3.5}$ profile from $30$ to $80$ kpc with no signs of a truncation or sharp break over this range. The radial profile measurement is largely insensitive to individual halo substructure components, but we find that attempting to measure the shape of the halo is overwhelmed by the Sagittarius stream such that no ellipsoidal shape is a satisfactory description in this region. These measurements allow us to begin placing the Milky Way in context with the growing sample of external galaxies where similar halo profile measurements are available, with the goal of further linking the properties of stellar halos to the accretion histories that formed them.

We study the vertical structure of a stellar disk obtained from a fully cosmological high-resolution hydrodynamical simulation of the formation of a Milky Way-like galaxy. At the present day, the disk's mean vertical height shows a well-defined and strong pattern, with amplitudes as large as 3 kpc in its outer regions. This pattern is the result of a satellite - host halo - disk interaction and reproduces, qualitatively, many of the observable properties of the Monoceros Ring. In particular we find disk material at the distance of Monoceros ($R \sim$ 12-16 kpc, galactocentric) extending far above the mid plane ( 30$^{\circ}$, $\langle Z \rangle \sim$ 1-2 kpc) in both hemispheres, as well as well-defined arcs of disk material at heliocentric distances $\gtrsim 5$ kpc. The pattern was first excited $\approx 3$ Gyr ago as an $m=1$ mode that later winds up into a leading spiral pattern. Interestingly, the main driver behind this perturbation is a low-mass low-velocity fly-by encounter. The satellite has total mass, pericentre distance and pericentric velocity of $\sim 5\%$ of the host, $\sim 80$ kpc, and 215 km/s, respectively. The satellite is not massive enough to directly perturb the galactic disk but we show that the density field of the host dark matter halo responds to this interaction resulting in a strong amplification of the perturbative effects. This subsequently causes the onset and development of the Monoceros-like feature.

Using the Pan-STARRS1 survey, we derive limiting magnitude, spatial completeness and density maps that we use to probe the three dimensional structure and estimate the stellar mass of the so-called Monoceros Ring. The Monoceros Ring is an enormous and complex stellar sub-structure in the outer Milky Way disk. It is most visible across the large Galactic Anticenter region, 120 < l < 240 degrees, -30 < b < +40 degrees. We estimate its stellar mass density profile along every line of sight in 2 X 2 degree pixels over the entire 30,000 square degree Pan-STARRS1 survey using the previously developed MATCH software. By parsing this distribution into a radially smooth component and the Monoceros Ring, we obtain its mass and distance from the Sun along each relevant line of sight. The Monoceros Ring is significantly closer to us in the South (6 kpc) than in the North (9 kpc). We also create 2D cross sections parallel to the Galactic plane that show 135 degrees of the Monoceros Ring in the South and 170 degrees of the Monoceros Ring in the North. We show that the Northern and Southern structures are also roughly concentric circles, suggesting that they may be a wave rippling from a common origin. Excluding the Galactic plane, we observe an excess stellar mass of 4 million solar masses across 120 < l < 240 degrees. If we interpolate across the Galactic plane, we estimate that this region contains 8 million solar masses. If we assume (somewhat boldly) that the Monoceros Ring is a set of two Galactocentric rings, its total stellar mass is 60 million solar masses. Finally, if we assume that it is a set of two circles centered at a point 4 kpc from the Galactic center in the anti-central direction, as our data suggests, we estimate its stellar mass to be 40 million solar masses.

The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline.

The Ophiuchus stellar stream presents a dynamical puzzle: its old stellar populations ($\sim 12$ Gyr) cannot be reconciled with (1) its orbit in a simple model for the Milky Way potential and (2) its short angular extent, both of which imply that the observed stream formed within the last $<1$ Gyr. Recent theoretical work has shown that streams on chaotic orbits may abruptly fan out near their apparent ends; stars in these fans are dispersed in both position and velocity and may be difficult to associate with the stream. Here we present the first evidence of such stream-fanning in the Ophiuchus stream, traced by four blue horizontal branch (BHB) stars beyond the apparent end of the stream. These stars stand out from the background by their high velocities ($v_{\rm los} > 230$ km s$^{-1}$) against $\sim 40$ other stars: their velocities are comparable to those of the stream, but would be exceptional if they were unrelated halo stars. Their positions and velocities are, however, inconsistent with simple extrapolation of the observed cold, high-density portion of the stream. These observations suggest that stream-fanning may be a real, observable effect and, therefore, that Ophiuchus may be on a chaotic orbit. They also show that the Ophiuchus stream is more extended and hence dynamically older than previously thought, easing the stellar population vs. dynamical age tension.

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3pi survey. The first satellite, Laevens 3, is located at a heliocentric distance of d=67+/-3 kpc. With a total magnitude of Mv=-4.4+/-0.3 and a half-light radius rh=7+/-2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4 (Dra II), is a closer and fainter satellite (d~20 kpc, Mv =-2.9+/-0.8), whose uncertain size (rh = 19 +8/-6 pc) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (rh = 37 +9/-8 pc and Mv=-5.2+/-0.4). At a heliocentric distance of 67+/-5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph.

We present the results of a deep study of the isolated dwarf galaxies Andromeda XXVIII and Andromeda XXIX with Gemini/GMOS and Keck/DEIMOS. Both galaxies are shown to host old, metal-poor stellar populations with no detectable recent star formation, conclusively identifying both of them as dwarf spheroidal galaxies (dSphs). And XXVIII exhibits a complex horizontal branch morphology, which is suggestive of metallicity enrichment and thus an extended period of star formation in the past. Decomposing the horizontal branch into blue (metal poor, assumed to be older) and red (relatively more metal rich, assumed to be younger) populations shows that the metal rich are also more spatially concentrated in the center of the galaxy. We use spectroscopic measurements of the Calcium triplet, combined with the improved precision of the Gemini photometry, to measure the metallicity of the galaxies, confirming the metallicity spread and showing that they both lie on the luminosity-metallicity relation for dwarf satellites. Taken together, the galaxies exhibit largely typical properties for dSphs despite their significant distances from M31. These dwarfs thus place particularly significant constraints on models of dSph formation involving environmental processes such as tidal or ram pressure stripping. Such models must be able to completely transform the two galaxies into dSphs in no more than two pericentric passages around M31, while maintaining a significant stellar populations gradient. Reproducing these features is a prime requirement for models of dSph formation to demonstrate not just the plausibility of environmental transformation but the capability of accurately recreating real dSphs.

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 pi imaging data and confirmed with follow-up wide-field photometry from the Large Binocular Cameras. The stellar system, with an absolute magnitude of M_V=-1.8 +/-0.5, a heliocentric distance of 30 +2/-2 kpc, and a half-mass radius of 34 +9/-8 pc, shows remarkable similarity to faint, nearby, small satellites such as Willman 1, Segue 1, Segue 2, and Bo\"otes II. The discovery of Laevens 2/Triangulum II further populates the region of parameter space for which the boundary between dwarf galaxies and globular clusters becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar structures.

We combine observations of the Local Group with data from the NASA-Sloan Atlas to show the variation in the quenched fraction of satellite galaxies from low mass dwarf spheroidals and dwarf irregulars to more massive dwarfs similar to the Magellanic clouds. While almost all of the low mass ($M_\star \lesssim 10^7$ $M_\odot$) dwarfs are quenched, at higher masses the quenched fraction decreases to approximately 40-50%. This change in the quenched fraction is large, and suggests a sudden change in the effectiveness of quenching that correlates with satellite mass. We combine this observation with models of satellite infall and ram pressure stripping to show that the low mass satellites must quench within 1-2 Gyr of pericenter passage to maintain a high quenched fraction, but that many more massive dwarfs must continue to form stars today even though they likely fell in to their host >5 Gyr ago. We also characterize how the susceptibility of dwarfs to ram pressure must vary as a function of mass if it is to account for the change in quenched fractions. Though neither model predicts the quenching effectiveness a priori, this modeling illustrates the physical requirements that the observed quenched fractions place on possible quenching mechanisms.

We present a panoptic view of the stellar structure in the Galactic disk's outer reaches commonly known as the Monoceros Ring, based on data from Pan-STARRS1. These observations clearly show the large extent of the stellar overdensities on both sides of the Galactic disk, extending between b = -25 and b = +35 degrees and covering over 130 degrees in Galactic longitude. The structure exhibits a complex morphology with both stream-like features and a sharp edge to the structure in both the north and the south. We compare this map to mock observations of two published simulations aimed at explaining such structures in the outer stellar disk, one postulating an origin as a tidal stream and the other demonstrating a scenario where the disk is strongly distorted by the accretion of a satellite. These morphological comparisons of simulations can link formation scenarios to observed structures, such as demonstrating that the distorted-disk model can produce thin density features resembling tidal streams. Although neither model produces perfect agreement with the observations--the tidal stream predicts material at larger distances which is not detected while in the distorted disk model the midplane is warped to an excessive degree--future tuning of the models to accommodate these latest data may yield better agreement.

We report the discovery of a thin stellar stream found in Pan-STARRS1 photometry near the Galactic bulge in the constellation of Ophiuchus. It appears as a coherent structure in the colour-selected stellar density maps produced to search for tidal debris around nearby globular clusters. The stream is exceptionally short and narrow; it is about 2.5{\deg} long and 6' wide in projection. The colour-magnitude diagram of this object, which harbours a blue horizontal-branch, is consistent with an old and relatively metal-poor population ([Fe/H]~-1.3) located 9.5 +/- 0.9 kpc away at (l,b) ~ (5{\deg},+32{\deg}), and 5.0 +/- 1.0 kpc from the Galactic centre. These properties argue for a globular cluster as progenitor. The finding of such a prominent, nearby stream suggests that many streams could await discovery in the more densely populated regions of our Galaxy.

We present the fiducial sequences of a sample of Galactic star clusters in the five bands of the Pan-STARRS1 (PS1) photometric system (g_P1, r_P1, i_P1, z_P1, y_P1). These empirical sequences -- which include the red giant and sub-giant branches, the main sequence, and the horizontal branch -- were defined from deep colour-magnitude diagrams reaching below the oldest main-sequence turn-offs of 13 globular and 3 old open clusters covering a wide range of metallicities (-2.4 < [Fe/H] < +0.4). We find excellent agreement for the nine clusters in common with previous studies in similar photometric systems when transformed to the PS1 system. Because the photometric and spectroscopic properties of these stellar populations are accurately known, the fiducials provide a solid basis for the interpretation of observations in the PS1 system, as well as valuable constraints to improve the empirical colour--$T_{eff}$ relations.

We present a new satellite in the outer halo of the Galaxy, the first Milky Way satellite found in the stacked photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1) Survey. From follow-up photometry obtained with WFI on the MPG/ESO 2.2m telescope, we argue that the object, located at a heliocentric distance of 145+/-17 kpc, is the most distant Milky Way globular cluster yet known. With a total magnitude of M_V=-4.3+/-0.2 and a half-light radius of 20+/-2 pc, it shares the properties of extended globular clusters found in the outer halo of our Galaxy and the Andromeda galaxy. The discovery of this distant cluster shows that the full spatial extent of the Milky Way globular cluster system has not yet been fully explored.

Using deep 21-cm HI data from the Green Bank Telescope we have detected an ~18.3 kpc-long gaseous extension associated with the starbursting dwarf galaxy IC 10. The newly-found feature stretches 1.3 deg to the northwest and has a large radial velocity gradient reaching to ~65 km/s lower than the IC 10 systemic velocity. A region of higher column density at the end of the extension that possesses a coherent velocity gradient (~10 km/s across ~26 arcmin) transverse to the extension suggests rotation and may be a satellite galaxy of IC 10. The HI mass of IC 10 is 9.5x10^7 (d/805 kpc)^2 Msun and the mass of the new extension is 7.1x10^5 (d/805 kpc)^2 Msun. An IC 10-M31 orbit using known radial velocity and proper motion values for IC 10 show that the HI extension is inconsistent with the trailing portion of the orbit so that an M31-tidal or ram pressure origin seems unlikely. We argue that the most plausible explanation for the new feature is that it is the result of a recent interaction (and possible late merger) with another dwarf galaxy. This interaction could not only have triggered the origin of the recent starburst in IC 10, but could also explain the existence of previously-found counter-rotating HI gas in the periphery of the IC 10 which was interpreted as originating from primordial gas infall.

We present the discovery of a new dwarf galaxy, Perseus I/Andromeda XXXIII, found in the vicinity of Andromeda (M31) in stacked imaging data from the Pan-STARRS1 3{\pi} survey. Located 27.9{\deg} away from M31, Perseus I has a heliocentric distance of 785 +/- 65 kpc, compatible with it being a satellite of M31 at 374 +14/-10 kpc from its host. The properties of Perseus I are typical for a reasonably bright dwarf galaxy (M_V = -10.3 +/- 0.7), with an exponential half-light radius of r_h = 1.7 +/- 0.4 arcminutes or r_h = 400 +105/-85 pc at this distance, and a moderate ellipticity (\epsilon = 0.43 +0.15/-0.17). The late discovery of Perseus I is due to its fairly low surface brightness (\mu_0=25.7 +1.0/-0.9 mag/arcsec^2), and to the previous lack of deep, high quality photometric data in this region. If confirmed to be a companion of M31, the location of Perseus I, far east from its host, could place interesting constraints on the bulk motion of the satellite system of M31.

A number of mechanisms have been proposed to connect star-forming dwarf irregular galaxies with the formation of non-star-forming dwarf spheroidal galaxies, but distinguishing between these mechanisms has been difficult. We use the Via Lactea dark matter only cosmological simulations to test two well-motivated simple hypotheses---transformation of irregulars into dwarf spheroidal galaxies by tidal stirring and ram pressure stripping following a close passage to the host galaxy, and transformation via mergers between dwarfs---and predict the radial distribution and inferred formation times of the resulting dwarf spheroidal galaxies. We compare this to the observed distribution in the Local Group and show that 1) the observed dSph distribution far from the Galaxy or M31 can be matched by the VL halos that have passed near the host galaxy at least once, though significant halo-to-halo scatter exists, 2) models that require two or more pericenter passages for dSph-formation cannot account for the dSphs beyond 500 kpc such as Cetus and Tucana, and 3) mergers predict a flat radial distribution of dSphs and cannot account for the high dSph fraction near the Galaxy, but are not ruled out at large distances. The models also suggest that for dSphs found today beyond 500 kpc, mergers tend to occur significantly earlier than dwarf--host encounters, thus leading to a potentially observable difference in stellar populations. We argue that tidal interactions are sufficient to reproduce the observed distribution of dSphs if and only if a single pericenter passage is sufficient to form a dSph.

We report the discovery of two new dwarf galaxies, Lacerta I/Andromeda XXXI (Lac I/And XXXI) and Cassiopeia III/Andromeda XXXII (Cas III/And XXXII), in stacked Pan-STARRS1 r_P1- and i_P1-band imaging data. Both are luminous systems (M_V ~ -12) located at projected distances of 20.3{\deg} and 10.5{\deg} from M31. Lac I and Cas III are likely satellites of the Andromeda galaxy with heliocentric distances of 756^{+44}_{-28} kpc and 772^{+61}_{-56} kpc, respectively, and corresponding M31-centric distances of 275+/-7 kpc and 144^{+6}_{-4} kpc . The brightest of recent Local Group member discoveries, these two new dwarf galaxies owe their late discovery to their large sizes (r_h = 4.2^{+0.4}_{-0.5} arcmin or 912^{+124}_{-93} pc for Lac I; r_h = 6.5^{+1.2}_{-1.0} arcmin or 1456+/-267 pc for Cas III), and consequently low surface brightness (\mu_0 ~ 26.0 mag/arcsec^2), as well as to the lack of a systematic survey of regions at large radii from M31, close to the Galactic plane. This latter limitation is now alleviated by the 3{\pi} Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies.

We use data from the Pan-STARRS1 survey to present a panoramic view of the Sagittarius tidal stream in the southern Galactic hemisphere. As a result of the extensive sky coverage of Pan-STARRS1, the southern stream is visible along more than 60 degrees of its orbit, nearly double the length seen by the SDSS. The recently discovered southern bifurcation of the stream is also apparent, with the fainter branch of the stream visible over at least 30 degrees. Using a combination of fitting both the main sequence turn-off and the red clump, we measure the distance to both arms of the stream in the south. We find that the distances to the bright arm of the stream agree very well with the N-body models of Law & Majewski (2010). We also find that the faint arm lies ~5 kpc closer to the Sun than the bright arm, similar to the behavior seen in the northern hemisphere.

We report the discovery of a new dwarf galaxy, Andromeda XXVIII, using data from the recently-released SDSS DR8. The galaxy is a likely satellite of Andromeda, and, at a separation of $365^{+17}_{-1}$ kpc, would be one of the most distant of Andromeda's satellites. Its heliocentric distance is $650^{+150}_{-80}$ kpc, and analysis of its structure and luminosity show that it has an absolute magnitude of $M_V = -8.5^{+0.4}_{-1.0}$ and half-light radius of $r_h = 210^{+60}_{-50}$ pc, similar to many other faint Local Group dwarfs. With presently-available imaging we are unable to determine if there is ongoing or recent star formation, which prevents us from classifying it as a dwarf spheroidal or dwarf irregular.

We report the discovery of a new dwarf galaxy, Andromeda XXIX, using data from the recently-released Sloan Digital Sky Survey DR8, and confirmed by Gemini North telescope Multi-Object Spectrograph imaging data. And XXIX appears to be a dwarf spheroidal galaxy, separated on the sky by a little more than 15 degrees from M31, with a distance inferred from the tip of the red giant branch of 730kpc+/-75kpc, corresponding to a three dimensional separation from M31 of between 205kpc and 227kpc (close to M31's virial radius). Its absolute magnitude, as determined by comparison to the red giant branch luminosity function of the Draco dwarf spheroidal, is M_V = -8.3+/-0.4. And XXIX's stellar populations appear very similar to Draco's; consequently, we estimate a metallicity for And XXIX of [Fe/H]\sim-1.8. The half-light radius of And XXIX is 360pc+/-60pc and its ellipticity is 0.35+/-0.06, typical of dwarf satellites of the Milky Way and M31 at this absolute magnitude range.

We present a study of the dust properties of 12 classical and superbubble HII regions in the Large Magellanic Cloud. We use infrared photometry from Spitzer (8, 24, 70, and 160 \mum bands), obtained as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) program, along with archival spectroscopic classifications of the ionizing stars to examine the role of stellar sources on dust heating and processing. Our infrared observations show surprisingly little correlation between the emission properties of the dust and the effective temperatures or bolometric magnitudes of stars in the HII regions, suggesting that the HII region evolutionary timescale is not on the order of the dust processing timescale. We find that the infrared emission of superbubbles and classical HII regions shows little differentiation between the two classes, despite the significant differences in age and morphology. We do detect a correlation of the 24 \mum emission from hot dust with the ratio of 70 to 160 \mum flux. This correlation can be modeled as a trend in the temperature of a minority hot dust component, while a majority of the dust remains significantly cooler.

We present a means of characterizing and removing internal reflections between the CCD and other optical surfaces in an astronomical camera. The stellar reflections appear as out-of-focus images and are not necessarily axisymmetric about the star. Using long exposures of very bright stars as calibration images we are able to measure the position, size, and intensity of reflections as a function of their position on the field. We also measure the extended stellar point-spread function out to one degree. Together this information can be used to create an empirical model of the excess light from bright stars and reduce systematic artifacts in deep surface photometry. We then reduce a set of deep observations of the Virgo cluster with our method to demonstrate its efficacy and to provide a comparison with other strategies for removing scattered light.