We present extremely deep imaging of the well-known edge-on galaxy NGC 4565 (the "Needle Galaxy"), observed as part of the Dragonfly Edge-on Galaxies Survey. Our imaging traces starlight over scales of ~100 kpc to surface brightnesses of ~32 mag/arcsec^2 in azimuthally averaged radial profiles. In narrow slice profiles along the major axis (with bin sizes ranging from 1.7 x 0.5 kpc to 1.7 x 7.8 kpc), we achieve a depth of ~29-30 mag/arcsec^2. We confirm the previously observed asymmetric disc truncation. In addition, the sharp northwest truncation turns over to a shallower component that coincides with a fan-like feature seen to wrap around the northwest disc limb. We propose that the fan may be a tidal ribbon, and qualitatively replicate the fan with simple simulations. Alternative explanations of the fan and the disc's asymmetry are considered as well. We conclude that we are most likely seeing evidence for accretion-based outer disk growth in NGC 4565.

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.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

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.

Deep observations of faint surface brightness stellar tidal streams in external galaxies with LSST are addressed in this White Paper contribution. We propose using the Wide--Fast--Deep survey that contains several nearby galaxies (at distances where the stars themselves are not resolved, i.e., beyond 20 Mpc). In the context of hierarchical galaxy formation, it is necessary to understand the prevalence and properties of tidal substructure around external galaxies based on integrated (i.e., unresolved) diffuse light. This requires collecting observations on much larger samples of galaxies than the Milky Way and M31. We will compare the observed structures to the predictions of cosmological models of galactic halo formation that inform us about the number and properties of streams around Milky Way-like galaxies. The insight gained from these comparisons will allow us to infer the properties of stream progenitors (masses, dynamics, metallicities, stellar populations). The changes in the host galaxies caused by the interactions with the dissolving companion galaxies will be another focus of our studies. We conclude by discussing synergies with WFIRST and Euclid, and also provide concrete suggestions for how the effects of scattered light could be minimized in LSST images to optimize the search for low surface brightness features, such as faint unresolved stellar tidal streams.

Tidal debris structures formed from disrupted satellites contain important clues about the assembly histories of galaxies. To date, studies of these structures have been hampered by reliance on by-eye identification and morphological classification which leaves their interpretation significantly uncertain. In this work we present a new machine-vision technique based on the Subspace-Constrained Mean Shift (SCMS) algorithm which can perform these tasks automatically. SCMS finds the location of the high-density `ridges' that define substructure morphology. After identification, the coefficients of an orthogonal series density estimator are used to classify points on the ridges as part of a continuum between shell-like or stream-like debris, from which a global morphological classification can be determined. We dub this procedure Subspace Constrained Unsupervised Detection of Structure (SCUDS). By applying this tool to controlled N--body simulations of minor mergers we demonstrate that the extracted classifications correspond to the well--understood underlying physics of phase mixing. The application of SCUDS to resolved stellar population data from near-future surveys will inform our understanding of the buildup of galaxies stellar halos.

We present a new distance estimation for the Milky Way dwarf spheroidal satellite Sculptor obtained from multi-epoch mid-infrared observations of RR Lyrae stars. The 3.6 {\mu}m observations have been acquired with the Infrared Array Camera on board the Spitzer Space Telescope as part of the SMHASH Program. Mid-infrared light curves for 42 RRL were obtained, from which we measured Sculptor's distance modulus to be {\mu} = 19.60 $\pm$ 0.02 (statistical) $\pm$ 0.04 (photometric) mag (with $\sigma_{sys}=$ = 0.09 mag), using the 3.6 {\mu}m empirical period-luminosity relations derived from the Galactic globular cluster M4, or {\mu} = 19.57 $\pm$ 0.02 (statistical) $\pm$ 0.04 (photometric) mag (with $\sigma_{sys}=$ = 0.11 mag) using empirical relations in the same passband recently derived from the Large Magellanic Cloud globular cluster Reticulum. Both these measurements are in good agreement with values presented in previous works with Sculptor RR Lyrae stars in optical bands, and are also consistent with recent near-infrared RR Lyrae results. Best agreement with the literature is found for the latter modulus which is equivalent to a distance of d = 82 $\pm$ 1 (statistical) $\pm$ 2 (photometric) kpc (with $\sigma_{sys}=$ = 4 kpc). Finally, using a subsample of RR Lyrae stars with spectroscopic metallicities, we demonstrate that these distance estimates are not affected by metallicity effects.

We present 1,201 galaxies at $0.05<z<0.45$ that host tidal features, detected from the first $\sim\! 200$ deg$^2$ of imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). All galaxies in the present sample have spectroscopic observations from the Sloan Digital Sky Survey (SDSS) spectroscopic campaigns, generating a sample of 21208 galaxies. Of these galaxies, we identify 214 shell systems and 987 stream systems. For 575 of these systems, we are additionally able to measure the $(g-i)$ colors of the tidal features. We find evidence for star formation in a subset of the streams, with the exception of streams around massive ellipticals, and find that stream host galaxies span the full range of stellar masses in our sample. Galaxies which host shells are predominantly red and massive: we find that observable shells form more frequently around ellipticals than around disc galaxies of the same stellar mass. Although the majority of the shells in our sample are consistent with being formed by minor mergers, $15\% \pm 4.4\%$ of shell host galaxies have $(g-i)$ colors as red as their host galaxy, consistent with being formed by major mergers. These "red shells" are additionally preferentially aligned with the major axis of the host galaxy, as previously predicted from simulations. We suggest that although the bulk of the observable shell population originates from fairly minor mergers, which preferentially form shells that are not aligned with the major axis of the galaxy, major mergers produce a significant number of observable shells.

We study tidal features (TFs) around galaxies in the REsolved Spectroscopy of a Local VolumE (RESOLVE) survey. Our sample consists of 1048 RESOLVE galaxies that overlap with the DECam Legacy Survey, which reaches an r-band 3$\sigma$ depth of $\sim$27.9 mag arcsec$^{-2}$ for a 100 arcsec$^{2}$ feature. Images were masked, smoothed, and inspected for TFs like streams, shells, or tails/arms. We find TFs in 17$^{\pm 2} \%$ of our galaxies, setting a lower limit on the true frequency. The frequency of TFs in the gas-poor (gas-to-stellar mass ratio $<$ 0.1) subsample is lower than in the gas-rich subsample (13$^{\pm 3} \%$ vs. 19$^{\pm 2} \%$). Within the gas-poor subsample, galaxies with TFs have higher stellar and halo masses, $\sim 3\times$ closer distances to nearest neighbors (in the same group), and possibly fewer group members at fixed halo mass than galaxies without TFs, but similar specific star formation rates. These results suggest TFs in gas-poor galaxies are typically streams/shells from dry mergers or satellite disruption. In contrast, the presence of TFs around gas-rich galaxies does not correlate with stellar or halo mass, suggesting these TFs are often tails/arms from resonant interactions. Similar to TFs in gas-poor galaxies, TFs in gas-rich galaxies imply 1.7x closer nearest neighbors in the same group; however, TFs in gas-rich galaxies are associated with diskier morphologies, higher star formation rates, and higher gas content. In addition to interactions with known neighbors, we suggest that TFs in gas-rich galaxies may arise from accretion of cosmic gas and/or gas-rich satellites below the survey limit.

Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 \mu m}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals.

We study the relationships between galaxy environments and galaxy properties related to disk (re)growth, considering two highly complete samples that are approximately baryonic mass limited into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. We quantify galaxy environments using both group identification and smoothed galaxy density field methods. We use by-eye and quantitative morphological classifications plus atomic gas content measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galaxies, and UV-bright disk host galaxies all become distinctly more common below group halo mass ~10^11.5 Msun, implying that this low group halo mass regime may be a preferred regime for significant disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue early types can regrow late-type disks. In fact, we find that the only significant difference in the typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different colors, where at fixed baryonic mass red early and late types have higher typical group halo masses than blue early and late types. More generally, we argue that the traditional morphology-environment relation (i.e., that denser environments tend to have more early types) can be largely attributed to the morphology-galaxy mass relation for centrals and the color-environment relation for satellites.

How do galaxies move relative to one another? While we can examine the motion of dark matter subhalos around their hosts in simulations of structure formation, determining the orbits of satellites around their parent galaxies from observations is impossible except for a small number of nearby cases. In this work we outline a novel approach to probing the orbital distributions of infalling satellite galaxies using the morphology of tidal debris structures. It has long been understood that the destruction of satellites on near-radial orbits tends to lead to the formation of shells of debris, while those on less eccentric orbits produce tidal streams. We combine an understanding of the scaling relations governing the orbital properties of debris with a simple model of how these orbits phase-mix over time to produce a `morphology metric' that more rigorously quantifies the conditions required for shells to be apparent in debris structures as a function of the satellite's mass and orbit and the interaction time. Using this metric we demonstrate how differences in orbit distributions can alter the relative frequency of shells and stream structures observed around galaxies. These experiments suggest that more detailed modeling and careful comparisons with current and future surveys of low surface brightness features around nearby galaxies should be capable of actually constraining orbital distributions and provide new insights into our understanding of structure formation.

The dark matter halo of the Milky Way is expected to be triaxial and filled with substructure. It is hoped that streams or shells of stars produced by tidal disruption of stellar systems will provide precise measures of the gravitational potential to test these predictions. We develop a method for inferring the Galactic potential with tidal streams based on the idea that the stream stars were once close in phase space. Our method can flexibly adapt to any form for the Galactic potential: it works in phase-space rather than action-space and hence relies neither on our ability to derive actions nor on the integrability of the potential. Our model is probabilistic, with a likelihood function and priors on the parameters. The method can properly account for finite observational uncertainties and missing data dimensions. We test our method on synthetic datasets generated from N-body simulations of satellite disruption in a static, multi-component Milky Way including a triaxial dark matter halo with observational uncertainties chosen to mimic current and near-future surveys of various stars. We find that with just four well-measured stream stars, we can infer properties of a triaxial potential with precisions of order 5-7 percent. Without proper motions we obtain 15 percent constraints on potential parameters and precisions around 25 percent for recovering missing phase-space coordinates. These results are encouraging for the eventual goal of using flexible, time-dependent potential models combined with larger data sets to unravel the detailed shape of the dark matter distribution around the Milky Way.