We present the first spectroscopic measurements of the ATLAS and Aliqa Uma streams from the Southern Stellar Stream Spectroscopic Survey ($S^5$), in combination with the photometric data from the Dark Energy Survey and astrometric data from $Gaia$. From the coherence of spectroscopic members in radial velocity and proper motion, we find out that these two systems are physically one stream with discontinuity in morphology and density on the sky (the "kink" feature). We refer to this entire stream as ATLAS-Aliqa Uma stream, or AAU stream. We perform a comprehensive exploration of the effect of baryonic substructures and find that only an encounter with the Sagittarius dwarf $\sim 0.5$ Gyr ago can create a feature similar to the observed "kink". In addition, we also identify two gaps in the ATLAS component associated with the broadening in the stream width (the "broadening" feature). These gaps have likely been created by small mass perturbers, such as dark matter halos, as AAU stream is the most distant cold stream known with severe variations in both the stream surface density and the stream track on the sky. With the stream track, stream distance and kinematic information, we determine the orbit of the AAU stream and find that it has been affected by the Large Magellanic Cloud, resulting in a misalignment between the proper motion and stream track. Together with the Orphan-Chenab Stream, AAU is the second stream pair that has been found to be a single stream separated into two segments by external perturbation.

The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.

We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

Stars slingshotted by the supermassive black hole at the Galactic centre will escape the Milky Way so quickly that their trajectories will be almost straight lines. Previous works have shown how these `hypervelocity stars' are subsequently deflected by the gravitational field of the Milky Way and the Large Magellanic Cloud (LMC), but have neglected to account for the reflex motion of the Milky Way in response to the fly by of the LMC. A consequence of this motion is that the hypervelocity stars we see on the outskirts of the Milky Way today were ejected from where the Milky Way centre was hundreds of millions of years ago. This change in perspective causes large apparent deflections in the trajectories of the hypervelocity stars, which are of the same order as the deflections caused by the gravitational force of the Milky Way and LMC. We quantify these deflections by simulating the production of hypervelocity stars in an isolated Milky Way (with a spherical or flattened dark matter halo), in a fixed-in-place Milky Way with a passing LMC, and in a Milky Way which responds to the passage of the LMC. The proper motion precision necessary to measure these deflections will be possible with the combination of Gaia with the proposed GaiaNIR successor mission, and these measurements will unlock the hypervelocity stars as probes of the shape of the Milky Way, the mass of the LMC, and of the dance of these two galaxies.

Due to their close proximity, the Large and Small Magellanic Clouds (SMC/LMC) provide natural laboratories for understanding how galaxies form and evolve. With the goal of determining the structure and dynamical state of the SMC, we present new spectroscopic data for $\sim$ 3000 SMC red giant branch stars observed using the AAOmega spectrograph at the Anglo-Australian Telescope. We complement our data with further spectroscopic measurements from previous studies that used the same instrumental configuration and proper motions from the \textit{Gaia} Data Release 2 catalogue. Analysing the photometric and stellar kinematic data, we find that the SMC centre of mass presents a conspicuous offset from the velocity centre of its associated $\mbox{H\,{\sc i}}$ gas, suggesting that the SMC gas is likely to be far from dynamical equilibrium. Furthermore, we find evidence that the SMC is currently undergoing tidal disruption by the LMC within 2\,kpc of the centre of the SMC, and possibly all the way in to the very core. This is evidenced by a net outward motion of stars from the SMC centre along the direction towards the LMC and apparent tangential anisotropy at all radii. The latter is expected if the SMC is undergoing significiant tidal stripping, as we demonstrate using a suite of $N$-body simulations of the SMC/LMC system disrupting around the Milky Way. These results suggest that dynamical models for the SMC that assume a steady state will need to be revisited.

We present Magellan/IMACS spectroscopy of three recently discovered ultra-faint Milky Way satellites, Grus II, Tucana IV, and Tucana V. We measure systemic velocities of V_hel = -110.0 +/- 0.5 km/s, V_hel = 15.9 +/- 1.8 km/s, and V_hel = -36.2 +/-2.5 km/s for the three objects, respectively. Their large relative velocity differences demonstrate that the satellites are unrelated despite their close physical proximity. We determine a velocity dispersion for Tuc IV of sigma = 4.3^+1.7_-1.0 km/s, but we cannot resolve the velocity dispersions of the other two systems. For Gru II we place an upper limit (90% confidence) on the dispersion of sigma < 1.9 km/s, and for Tuc V we do not obtain any useful limits. All three satellites have metallicities below [Fe/H] = -2.1, but none has a detectable metallicity spread. We determine proper motions for each satellite based on Gaia astrometry and compute their orbits around the Milky Way. Gru II is on a tightly bound orbit with a pericenter of 25 kpc and orbital eccentricity of 0.45. Tuc V likely has an apocenter beyond 100 kpc, and could be approaching the Milky Way for the first time. The current orbit of Tuc IV is similar to that of Gru II, with a pericenter of 25 kpc and an eccentricity of 0.36. However, a backward integration of the position of Tuc IV demonstrates that it collided with the Large Magellanic Cloud at an impact parameter of 4 kpc ~120 Myr ago, deflecting its trajectory and possibly altering its internal kinematics. Based on their sizes, masses, and metallicities, we classify Gru II and Tuc IV as likely dwarf galaxies, but the nature of Tuc V remains uncertain.

Recent measurements suggest that the Large Magellanic Cloud (LMC) may weigh as much as 25\% of the Milky Way. In this work we explore how such a large satellite affects mass estimates of the Milky Way based on equilibrium modelling of the stellar halo or other tracers. In particular, we show that if the LMC is ignored, the Milky Way mass is overestimated by as much as 50\%. This bias is due to the bulk motion in the outskirts of the Galaxy's halo and can be, at least in part, accounted for with a simple modification to the equilibrium modelling. Finally, we show that the LMC has a substantial effect on the orbit Leo I which acts to increase its present day speed relative to the Milky Way. We estimate that accounting for a $1.5\times10^{11} M_\odot$ LMC would lower the inferred Milky Way mass to $\sim10^{12} M_\odot$.

We present new chemo--kinematics of the Hercules dwarf galaxy based on Keck II-- DEIMOS spectroscopy. Our 21 confirmed members have a systemic velocity of $v_{\mathrm{Herc}}=46.4\pm1.3$ kms$^{-1}$ and a velocity dispersion $\sigma_{v,\mathrm{Herc}}=4.4^{+1.4}_{-1.2}$ kms$^{-1}$. From the strength of the Ca II triplet, we obtain a metallicity of [Fe/H]= $-2.48\pm0.19$ dex and dispersion of $\sigma_{\rm{[Fe/H]}}= 0.63^{+0.18}_{-0.13}$ dex. This makes Hercules a particularly metal--poor galaxy, placing it slightly below the standard mass--metallicity relation. Previous photometric and spectroscopic evidence suggests that Hercules is tidally disrupting and may be on a highly radial orbit. From our identified members, we measure no significant velocity gradient. By cross--matching with the second \textit{Gaia} data release, we determine an uncertainty--weighted mean proper motion of $\mu_{\alpha}^*=\mu_{\alpha}\cos(\delta)=-0.153\pm{0.074}$ mas yr$^{-1}$, $\mu_{\delta}=-0.397\pm0.063$ mas yr$^{-1}$. This proper motion is slightly misaligned with the elongation of Hercules, in contrast to models which suggest that any tidal debris should be well aligned with the orbital path. Future observations may resolve this tension.

The GD-1 stream is one of the longest and coldest stellar streams discovered to date, and one of the best objects for constraining the dark matter properties of the Milky Way. Using data from {\it Gaia} DR2 we study the proper motions, distance, morphology and density of the stream to uncover small scale perturbations. The proper motion cleaned data shows a clear distance gradient across the stream, ranging from 7 to 12 kpc. However, unlike earlier studies that found a continuous gradient, we uncover a distance minimum at $\varphi_{1}\approx$-50 deg, after which the distance increases again. We can reliably trace the stream between -85$<\varphi_{1}<$15 deg, showing an even further extent to GD-1 beyond the earlier extension of \citet{Price-Whelan18a}. We constrain the stream track and density using a Boolean matched filter approach and find three large under densities and find significant residuals in the stream track lining up with these gaps. In particular, a gap is visible at $\varphi_{1}$=-3 deg, surrounded by a clear sinusoidal wiggle. We argue that this wiggle is due to a perturbation since it has the wrong orientation to come from a progenitor. We compute a total initial stellar mass of the stream segment of 1.58$\pm$0.07$\times$10$^{4}$ M$_{\odot}$. With the extended view of the spur in this work, we argue that the spur may be unrelated to the adjacent gap in the stream. Finally, we show that an interaction with the Sagittarius dwarf can create features similar to the spur.

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.

We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the star in the Galactic frame is $1755\pm50$ km/s. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A* with a velocity of $\sim 1800$ km/s and travelled for $4.8$ Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion $V_{y,\odot}= 246.1\pm 5.3$ km/s or position $R_0=8.12\pm 0.23$ kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disk of young stars at the Galactic centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time.

We analyze the distribution of stars along the GD-1 stream with a combination of data from the ${\it Gaia}$ satellite and the Pan-STARRS survey, and we show that the population of subhalos predicted by the cold dark matter paradigm are necessary and sufficient to explain the perturbations observed in the linear density of stars. This allows us to set novel constraints on alternative dark matter scenarios that predict a suppression of the subhalo mass function on scales smaller than the mass of dwarf galaxies. A combined analysis of the density perturbations in the GD-1 and Pal 5 streams leads to a $95\%$ lower limit on the mass of warm dark matter thermal relics $m_{\rm WDM}>4.6$ keV; adding dwarf satellite counts strengthens this to $m_{\rm WDM}>6.3$ keV.

New data from the $\textit{Gaia}$ satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyze the stellar density perturbations in the GD-1 stream and show that they cannot be due to known baryonic structures like giant molecular clouds, globular clusters, or the Milky Way's bar or spiral arms. A joint analysis of the GD-1 and Pal 5 streams instead requires a population of dark substructures with masses $\approx 10^{7}$ to $10^9 \ M_{\rm{\odot}}$. We infer a total abundance of dark subhalos normalised to standard cold dark matter $n_{\rm sub}/n_{\rm sub, CDM} = 0.4 ^{+0.3}_{-0.2}$ ($68 \%$), which corresponds to a mass fraction contained in the subhalos $f_{\rm{sub}} = 0.14 ^{+0.11}_{-0.07} \%$, compatible with the predictions of hydrodynamical simulation of cold dark matter with baryons.

We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ravi, Wambelong, Willka Yaku, and Turbio. We find evidence for a misalignment between stream tracks and the systemic proper motion of streams that may suggest a systematic gravitational influence from the Large Magellanic Cloud. These proper motions, when combined with radial velocity measurements, will allow for detailed orbit modeling which can be used to constrain properties of the LMC and its affect on nearby streams, as well as global properties of the Milky Way's gravitational potential.

We study the orbits of ultra-faint dwarf galaxies in the combined presence of the Milky Way and LMC and we find 6 dwarfs which were likely accreted with the LMC (Car 2, Car 3, Hor 1, Hyi 1, Phe 2, Ret 2), in addition to the SMC, representing strong evidence of dwarf galaxy group infall. This procedure depends on the gravitational pull of the LMC, thus allowing us to place a lower bound on the Cloud's mass of $M_{\rm LMC} > 1.24\times10^{11} M_\odot$. This mass estimate is validated by applying the technique to a cosmological zoom-in simulation of a Milky Way-like galaxy with an LMC analogue where we find that while this lower bound may be overestimated, it will improve in the future with smaller observational errors. We apply this technique to dwarf galaxies lacking radial velocities and find that Eri 3 has a broad range of radial velocities for which it has a significant chance ($> 0.4$) of having being bound to the Cloud. We study the non-Magellanic classical satellites and find that Fornax has an appreciable probability of being an LMC satellite if the LMC is sufficiently massive. In addition, we explore how the orbits of the Milky Way satellites change in the presence of the LMC and find a significant change for several objects. Finally, we find that the LMC satellites are slightly smaller than the Milky Way satellites at a fixed luminosity, possibly due to the different tidal environments they have experienced.

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 introduce a novel abundance matching technique that produces a more accurate estimate of the pre-infall halo mass, $M_{200}$, for satellite galaxies. To achieve this, we abundance match with the mean star formation rate, averaged over the time when a galaxy was forming stars, $\langle {\rm SFR}\rangle$, instead of the stellar mass, $M_*$. Using data from the Sloan Digital Sky Survey, the GAMA survey and the Bolshoi simulation, we obtain a statistical $\langle {\rm SFR}\rangle-{\rm M}_{200}$ relation in $\Lambda{\rm CDM}$. We then compare the pre-infall halo mass, $M^{\rm abund}_{200}$, derived from this relation with the pre-infall dynamical mass, $M^{\rm dyn}_{200}$, for 21 nearby dSph and dIrr galaxies, finding a good agreement between the two. As a first application, we use our new $\langle {\rm SFR}\rangle-{\rm M}_{200}$ relation to empirically measure the cumulative mass function of a volume-complete sample of bright Milky Way satellites within 280 kpc of the Galactic centre. Comparing this with a suite of cosmological 'zoom' simulations of Milky Way-mass halos that account for subhalo depletion by the Milky Way disc, we find no missing satellites problem above $M_{200} \sim 10^9\,{\rm M}_\odot$ in the Milky Way. We discuss how this empirical method can be applied to a larger sample of nearby spiral galaxies.

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/ .

Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the in-fall of the Magellanic Clouds.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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.

We present a spectroscopic study of the tidal tails and core of the Milky Way satellite Tucana III, collectively referred to as the Tucana III stream, using the 2dF+AAOmega spectrograph on the Anglo-Australian Telescope and the IMACS spectrograph on the Magellan/Baade Telescope. In addition to recovering the brightest 9 previously known member stars in the Tucana III core, we identify 22 members in the tidal tails. We observe strong evidence for a velocity gradient of 8.0 km/s/deg (or 18.3 km/s/kpc over at least 3$^{\circ}$ (or 1.3 kpc) on the sky. Based on the continuity in velocity we confirm that the Tucana III tails are real tidal extensions of Tucana III. The large velocity gradient of the stream implies that Tucana III is likely on a radial orbit. We successfully obtain metallicities for 4 members in the core and 12 members in the tails. We find that members close to the ends of the stream tend to be more metal-poor than members in the core, indicating a possible metallicity gradient between the center of the progenitor halo and its edge. The spread in metallicity suggests that the progenitor of the Tucana III stream is a dwarf galaxy rather than a star cluster. Furthermore, we find that with the precise photometry of the Dark Energy Survey data, there is a discernible color offset between metal-rich disk stars and metal-poor stream members. This metallicity-dependent color offers a more efficient method to recognize metal-poor targets and will increase the selection efficiency of stream members for future spectroscopic follow-up programs on stellar streams.

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.

We carried out a spectroscopic follow-up program of the four new stellar stream candidates detected by Belokurov & Koposov (2016) in the outskirts of the Large Magellanic Cloud (LMC) using FORS2 (VLT). The medium-resolution spectra were used to measure the line-of-sight velocities, estimate stellar metallicities and to classify stars into Blue Horizontal Branch (BHB) and Blue Straggler (BS) stars. Using the 4-D phase-space information, we attribute approximately one half of our sample to the Magellanic Clouds, while the rest is part of the Galactic foreground. Only two of the four stream candidates are confirmed kinematically. While it is impossible to estimate the exact levels of MW contamination, the phase-space distribution of the entire sample of our Magellanic stars matches the expected velocity gradient for the LMC halo and extends as far as 33 deg (angular separation) or 29 kpc from the LMC center. Our detections reinforce the idea that the halo of the LMC seems to be larger than previously expected, and its debris can be spread in the sky out to very large separations from the LMC center. Finally, we provide some kinematic evidence that many of the stars analysed here have likely come from the Small Magellanic Cloud.

Iorio et al. (2018) mapped out the Milky Way halo using a sample of RR Lyrae stars drawn from a cross-match of Gaia with 2MASS. We investigate the significant residual in their model which we constrain to lie at Galactocentric radii $12<R<27\;\mathrm{kpc}$ and extend over $2600\;\mathrm{deg}^2$ of the sky. A counterpart of this structure exists in both the Catalina Real Time Survey and the sample of RR Lyrae variables identified in Pan-STARRS by Hernitschek et al. (2016), demonstrating that this structure is not caused by the spatial inhomogeneity of Gaia. The structure is likely the Virgo Stellar Stream and/or Virgo Over-Density. We show the structure is aligned with the Magellanic Stream and suggest that it is either debris from a disrupted dwarf galaxy that was a member of the Vast Polar Structure or that it is SMC debris from a tidal interaction of the SMC and LMC $3\;\mathrm{Gyr}$ ago. If the latter then the sub-structure in Virgo may have a Magellanic origin.

We report the discovery of a thin stellar stream - which we name the Jet stream - crossing the constellations of Hydra and Pyxis. The discovery was made in data from the SLAMS survey, which comprises deep $g$ and $r$ imaging for a $650$ square degree region above the Galactic disc performed by the CTIO Blanco + DECam. SLAMS photometric catalogues will be made publicly available. The stream is approximately 0.18 degrees wide and 10 degrees long, though it is truncated by the survey footprint. Its colour-magnitude diagram is consistent with an old, metal-poor stellar population at a heliocentric distance of approximately 29 kpc. We corroborate this measurement by identifying a spatially coincident overdensity of likely blue horizontal branch stars at the same distance. There is no obvious candidate for a surviving stream progenitor.

We use astrometry and broad-band photometry from Data Release 2 of the ESA's Gaia mission to map out low surface-brightness features in the stellar density distribution around the Large and Small Magellanic Clouds. The LMC appears to have grown two thin and long stellar streams in its Northern and Southern regions, highly reminiscent of spiral arms. We use computer simulations of the Magellanic Clouds' in-fall to demonstrate that these arms were likely pulled out of the LMC's disc due to the combined influence of the SMC's most recent fly-by and the tidal field of the Milky Way.

Using proper motion measurements from Gaia DR2, we probe the origin of 26 previously known hypervelocity stars (HVSs) around the Milky Way. We find that a significant fraction of these stars have a high probability of originating close to the Milky Way centre, but there is one obvious outlier. HVS3 is highly likely to be coming almost from the centre of the Large Magellanic Cloud (LMC). During its closest approach, $21.1^{+6.1}_{-4.6}$ Myr ago, it had a relative velocity of $870^{+69}_{-66}$ kms$^{-1}$ with respect to the LMC. This large kick velocity is only consistent with the Hills mechanism, requiring a massive black hole at the centre of the LMC. This provides strong direct evidence that the LMC itself harbours a massive black hole of at least $4\times 10^3 -10^4 M_\odot$.

We present results of the first dynamical stream fits to the recently discovered Tucana III stream. These fits assume a fixed Milky Way potential and give proper motion predictions, which can be tested with the upcoming Gaia Data Release 2. These fits reveal that Tucana III is on an eccentric orbit around the Milky Way and, more interestingly, that Tucana III passed within 15 kpc of the Large Magellanic Cloud (LMC) approximately 75 Myr ago. Given this close passage, we fit the Tucana III stream in the combined presence of the Milky Way and the LMC. We find that the predicted proper motions depend on the assumed mass of the LMC and that the LMC can induce a substantial proper motion perpendicular to the stream track. A detection of this misalignment will directly probe the extent of the LMC's influence on our Galaxy, and has implications for nearly all methods which attempt to constraint the Milky Way potential. Such a measurement will be possible with the upcoming Gaia DR2, allowing for a measurement of the LMC's mass.

We report the discovery of a nearby dwarf galaxy in the constellation of Hydrus, between the Large and the Small Magellanic Clouds. Hydrus 1 is a mildy elliptical ultra-faint system with luminosity $M_V\sim$ -4.7 and size $\sim$ 50 pc, located 28 kpc from the Sun and 24 kpc from the LMC. From spectroscopy of $\sim$ 30 member stars, we measure a velocity dispersion of 2.7 km/s and find tentative evidence for a radial velocity gradient consistent with 3 km/s rotation. Hydrus 1's velocity dispersion indicates that the system is dark matter dominated, but its dynamical mass-to-light ratio M/L $\sim$ 66 is significantly smaller than typical for ultra-faint dwarfs at similar luminosity. The kinematics and spatial position of Hydrus~1 make it a very plausible member of the family of satellites brought into the Milky Way by the Magellanic Clouds. While Hydrus 1's proximity and well-measured kinematics make it a promising target for dark matter annihilation searches, we find no evidence for significant gamma-ray emission from Hydrus 1. The new dwarf is a metal-poor galaxy with a mean metallicity [Fe/H]=-2.5 and [Fe/H] spread of 0.4 dex, similar to other systems of similar luminosity. Alpha-abundances of Hyi 1 members indicate that star-formation was extended, lasting between 0.1 and 1 Gyr, with self-enrichment dominated by SN Ia. The dwarf also hosts a highly carbon-enhanced extremely metal-poor star with [Fe/H] $\sim$ -3.2 and [C/Fe] $\sim$ +3.0.

We use a new panoramic imaging survey, conducted with the Dark Energy Camera, to map the stellar fringes of the Large and Small Magellanic Clouds to extremely low surface brightness V $\gtrsim$ 32 mag arcsec$^{-2}$. Our results starkly illustrate the closely interacting nature of the LMC-SMC pair. We show that the outer LMC disk is strongly distorted, exhibiting an irregular shape, evidence for warping, and significant truncation on the side facing the SMC. Large diffuse stellar substructures are present both to the north and south of the LMC, and in the inter-Cloud region. At least one of these features appears co-spatial with the bridge of RR Lyrae stars that connects the Clouds. The SMC is highly disturbed -- we confirm the presence of tidal tails, as well as a large line-of-sight depth on the side closest to the LMC. Young, intermediate-age, and ancient stellar populations in the SMC exhibit strikingly different spatial distributions. In particular, those with ages $\sim$ 1.5-4 Gyr exhibit a spheroidal distribution with a centroid offset from that of the oldest stars by several degrees towards the LMC. We speculate that the gravitational influence of the LMC may already have been perturbing the gaseous component of the SMC several Gyr ago. With careful modeling, the variety of substructures and tidal distortions evident in the Magellanic periphery should tightly constrain the interaction history of the Clouds.

Using a large sample of Main Sequence stars with 7-D measurements supplied by Gaia and SDSS, we study the kinematic properties of the local (within ~10 kpc from the Sun) stellar halo. We demonstrate that the halo's velocity ellipsoid evolves strongly with metallicity. At the low [Fe/H] end, the orbital anisotropy (the amount of motion in the radial direction compared to the tangential one) is mildly radial with 0.2<beta<0.4. However, for stars with [Fe/H]>-1.7 we measure extreme values of beta~0.9. Across the metallicity range considered, i.e. -3<[Fe/H]<-1, the stellar halo's spin is minimal, at the level of 20<v_theta (km/s) <30. Using a suite of cosmological zoom-in simulations of halo formation, we deduce that the observed acute anisotropy is inconsistent with the continuous accretion of dwarf satellites. Instead, we argue, the stellar debris in the inner halo were deposited in a major accretion event by a satellite with Mvir>10^10 Msun around the epoch of the Galactic disc formation, i.e. between 8 and 11 Gyr ago. The radical halo anisotropy is the result of the dramatic radialisation of the massive progenitor's orbit, amplified by the action of the growing disc.

Efficient expansions of the gravitational field of (dark) haloes have two main uses in the modelling of galaxies: first, they provide a compact representation of numerically-constructed (or real) cosmological haloes, incorporating the effects of triaxiality, lopsidedness or other distortion. Secondly, they provide the basis functions for self-consistent field expansion algorithms used in the evolution of $N$-body systems. We present a new family of biorthogonal potential-density pairs constructed using the Hankel transform of the Laguerre polynomials. The lowest-order density basis functions are double-power-law profiles cusped like $\rho \sim r^{-2 + 1/\alpha}$ at small radii with asymptotic density fall-off like $\rho \sim r^{-3 -1/(2\alpha)}$. Here, $\alpha$ is a parameter satisfying $\alpha \ge 1/2$. The family therefore spans the range of inner density cusps found in numerical simulations, but has much shallower -- and hence more realistic -- outer slopes than the corresponding members of the only previously-known family deduced by Zhao (1996) and exemplified by Hernquist & Ostriker (1992). When $\alpha =1$, the lowest-order density profile has an inner density cusp of $\rho \sim r^{-1}$ and an outer density slope of $\rho \sim r^{-3.5}$, similar to the famous Navarro, Frenk & White (1997) model. For this reason, we demonstrate that our new expansion provides a more accurate representation of flattened NFW haloes than the competing Hernquist-Ostriker expansion. We utilize our new expansion by analysing a suite of numerically-constructed haloes and providing the distributions of the expansion coefficients.

Using deep photometric data from CFHT/Megacam, we study the morphology and density of the GD-1 stream, one of the longest and coldest stellar streams in the Milky Way. Our deep data recovers the lower main sequence of the stream with unprecedented quality, clearly separating it from Milky Way foreground and background stars. An analysis of the distance to different parts of the stream shows that GD-1 lies at a heliocentric distance between 8 and 10 kpc, with only a shallow gradient across 45 deg on the sky. Matched filter maps of the stream density show clear density variations, such as deviations from a single orbital track and tentative evidence for stream fanning. We also detect a clear under-density in the middle of the stream track at $\varphi_{1}$=-45 deg surrounded by overdense stream segments on either side. This location is a promising candidate for the elusive missing progenitor of the GD-1 stream. We conclude that the GD-1 stream has clearly been disturbed by interactions with the Milky Way disk or other sub-halos.

We present a new method for embedding a stellar disc in a cosmological dark matter halo and provide a worked example from a {\Lambda}CDM zoom-in simulation. The disc is inserted into the halo at a redshift z = 3 as a zero-mass rigid body. Its mass and size are then increased adiabatically while its position, velocity, and orientation are determined from rigid-body dynamics. At z = 1, the rigid disc is replaced by an N-body disc whose particles sample a three-integral distribution function (DF). The simulation then proceeds to z = 0 with live disc and halo particles. By comparison, other methods assume one or more of the following: the centre of the rigid disc during the growth phase is pinned to the minimum of the halo potential, the orientation of the rigid disc is fixed, or the live N-body disc is constructed from a two rather than three-integral DF. In general, the presence of a disc makes the halo rounder, more centrally concentrated, and smoother, especially in the innermost regions. We find that methods in which the disc is pinned to the minimum of the halo potential tend to overestimate the amount of adiabatic contraction. Additionally, the effect of the disc on the subhalo distribution appears to be rather insensitive to the disc insertion method. The live disc in our simulation develops a bar that is consistent with the bars seen in late-type spiral galaxies. In addition, particles from the disc are launched or "kicked up" to high galactic latitudes.

We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of $g > 23.5$ with a relative photometric calibration uncertainty of $< 1 \%$. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.

We explore the connection between galaxies and dark matter halos in the Milky Way (MW) and quantify the implications on properties of the dark matter particle and the phenomenology of low-mass galaxy formation. This is done through a probabilistic comparison of the luminosity function of MW dwarf satellite galaxies to models based on two suites of zoom-in simulations. One suite is dark-matter-only while the other includes a disk component, therefore we can quantify the effect of the MW's baryonic disk on our results. We apply numerous Stellar-Mass-Halo-Mass (SMHM) relations allowing for multiple complexities: scatter, a characteristic break scale, and subhalos which host no galaxy. In contrast to previous works we push the model/data comparison to the faintest dwarfs by modeling observational incompleteness, allowing us to draw three new conclusions. Firstly, we constrain the SMHM relation for $10^2<M_*/M_\odot<10^8$ galaxies, allowing us to bound the peak halo mass of the faintest MW satellite to $M_\mathrm{vir}>2.4\times10^8M_\odot$ ($1\sigma$). Secondly, by translating to a Warm Dark Matter (WDM) cosmology, we bound the thermal relic mass $m_\mathrm{WDM}>2.9$ keV at 95\% confidence, on a par with recent constraints from the Lyman-$\alpha$ forest. Lastly, we find that the observed number of ultra-faint MW dwarfs is in tension with the theoretical prediction that reionisation prevents galaxy formation in almost all $10^8M_\odot$ halos. This can be tested with the next generation of deep imaging surveys. To this end, we predict the likely number of detectable satellite galaxies in the Subaru/HSC survey and the LSST. Confronting these predictions with future observations will be amongst our strongest tests of WDM and the effect reionisation on low-mass systems.

We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.

We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.

We use a new contiguous imaging survey conducted using the Dark Energy Camera to investigate the distribution and properties of young stellar populations in the Magellanic inter-Cloud region. These young stars are strongly spatially clustered, forming a narrow chain of low-mass associations that trace the densest HI gas in the Magellanic Bridge and extend, in projection, from the SMC to the outer disk of the LMC. The associations in our survey footprint have ages $\lesssim 30$ Myr, masses in the range $\sim 100-1200\,{\rm M}_\odot$, and very diffuse structures with half-light radii of up to $\sim 100$ pc. The two most populous are strongly elliptical, and aligned to $\approx 10^{{\rm o}}$ with the axis joining the centres of the LMC and SMC. These observations strongly suggest that the young inter-Cloud populations formed in situ, likely due to the compression of gas stripped during the most recent close LMC-SMC encounter. The associations lie at distances intermediate between the two Clouds, and we find no evidence for a substantial distance gradient across the imaged area. Finally, we identify a vast shell of young stars surrounding a central association, that is spatially coincident with a low column density bubble in the HI distribution. The properties of this structure are consistent with a scenario where stellar winds and supernova explosions from massive stars in the central cluster swept up the ambient gas into a shell, triggering a new burst of star formation. This is a prime location for studying stellar feedback in a relatively isolated environment.

Only in the Milky Way is it possible to conduct an experiment which uses stellar streams to detect low-mass dark matter subhaloes. In smooth and static host potentials, tidal tails of disrupting satellites appear highly symmetric. However, perturbations from dark subhaloes, as well as from GMCs and the Milky Way bar, can induce density fluctuations that destroy this symmetry. Motivated by the recent release of unprecedentedly deep and wide imaging data around the Pal~5 stellar stream, we develop a new probabilistic, adaptive and non-parametric technique which allows us to bring the cluster's tidal tails into clear focus. Strikingly, we uncover a stream whose density exhibits visible changes on a variety of angular scales. We detect significant bumps and dips, both narrow and broad: two peaks on either side of the progenitor, each only a fraction of a degree across, and two gaps, $\sim2^{\circ}$ and $\sim9^{\circ}$ wide, the latter accompanied by a gargantuan lump of debris. This largest density feature results in a pronounced inter-tail asymmetry which cannot be made consistent with an unperturbed stream according to a suite of simulations we have produced. We conjecture that the sharp peaks around Pal 5 are epicyclic overdensities, while the two dips are consistent with impacts by subhaloes. Assuming an age of 3.4 Gyr for Pal 5, these two gaps would correspond to the characteristic size of gaps created by subhaloes in the mass range of $10^6-10^7 M_\odot$ and $10^7-10^8 M_\odot$ respectively. In addition to dark substructure, we find that the bar of the Milky Way can plausibly produce the asymmetric density seen in Pal 5 and that GMCs could cause the smaller gap.

We explore the possibility that the observed population of Galactic hypervelocity stars (HVSs) originate as runaway stars from the Large Magellanic Cloud (LMC). Pairing a binary evolution code with an N-body simulation of the interaction of the LMC with the Milky Way, we predict the spatial distribution and kinematics of an LMC runaway population. We find that runaway stars from the LMC can contribute Galactic HVSs at a rate of $3 \times 10^{-6}\;\mathrm{yr}^{-1}$. This is composed of stars at different points of stellar evolution, ranging from the main-sequence to those at the tip of the asymptotic giant branch. We find that the known B-type HVSs have kinematics which are consistent with an LMC origin. There is an additional population of hypervelocity white dwarfs whose progenitors were massive runaway stars. Runaways which are even more massive will themselves go supernova, producing a remnant whose velocity will be modulated by a supernova kick. This latter scenario has some exotic consequences, such as pulsars and supernovae far from star-forming regions, and a small rate of microlensing from compact sources around the halo of the LMC.

Tidal streams in the Milky Way are sensitive probes of the population of dark-matter subhalos predicted in cold-dark-matter (CDM) simulations. We present a new calculus for computing the effect of subhalo fly-bys on cold tidal streams based on the action-angle representation of streams. The heart of this calculus is a line-of-parallel-angle approach that calculates the perturbed distribution function of a given stream segment by undoing the effect of all impacts. This approach allows one to compute the perturbed stream density and track in any coordinate system in minutes for realizations of the subhalo distribution down to 10^5 Msun, accounting for the stream's internal dispersion and overlapping impacts. We study the properties of density and track fluctuations with suites of simulations. The one-dimensional density and track power spectra along the stream trace the subhalo mass function, with higher-mass subhalos producing power only on large scales, while lower mass subhalos cause structure on smaller scales. The time-dependence of impacts and of the evolution of the stream after an impact gives rise to bispectra. We further find that tidal streams are essentially corrugated sheets in the presence of subhalo perturbations: different projections of the track all reflect the same pattern of perturbations, facilitating their observational measurement. We apply this formalism to density data for the Pal 5 stream and make a first rigorous determination of 10^{+11}_{-6} dark-matter subhalos with masses between 3x10^6 and 10^9 Msun within 20 kpc from the Galactic center (corresponding to 1.4^{+1.6}_{-0.9} times the number predicted by CDM-only simulations or to f_{sub}(r<20 kpc) ~ 0.2%). Improved data will allow measurements of the subhalo mass function down to 10^5 Msun, thus definitively testing whether dark matter clumps on the smallest scales relevant for galaxy formation.

We exploit the first data release from the Gaia mission to identify candidate Mira variables in the outskirts of the Magellanic Clouds. The repeated observations of sources during the initial phase of the Gaia mission is used to identify stars that show signs of variability. This variability information, combined with infrared photometry from 2MASS and WISE, allows us to select a clean sample of giants in the periphery of the LMC. We find evidence for Miras surrounding the LMC out to ~20 deg in all directions, apart from the North-West quadrant. Our sample does not generally follow the gas distribution of the Magellanic system; Miras are notably absent in the gaseous bridge between the LMC and SMC, but they are likely related to the stellar RR Lyrae bridge reported by Belokurov et al. (2016). The stellar stream discovered by Mackey et al. (2016) to the North of the LMC is almost perfectly delineated by our Mira variables, and likely extends further East toward the Galactic plane. The presence of an intermediate-age population in this stream advocates an LMC disc origin. We also find a significant excess of Miras to the East of the LMC; these more diffusely distributed stars are likely stripped SMC stars due to interactions with the LMC. Miras are also identified in regions of the sky away from the Clouds; we locate stars likely associated with known massive substructures, and also find potential associations with stripped SMC debris above the Galactic plane.

We present the discovery of stellar tidal tails around the Large and the Small Magellanic Clouds in the Gaia DR1 data. In between the Clouds, their tidal arms are stretched towards each other to form an almost continuous stellar bridge. Our analysis relies on the exquisite quality of the Gaia's photometric catalogue to build detailed star-count maps of the Clouds. We demonstrate that the Gaia DR1 data can be used to detect variable stars across the whole sky, and in particular, RR Lyrae stars in and around the LMC and the SMC. Additionally, we use a combination of Gaia and Gale to follow the distribution of Young Main Sequence stars in the Magellanic System. Viewed by Gaia, the Clouds show unmistakable signs of interaction. Around the LMC, clumps of RR Lyrae are observable as far as ~20 degrees, in agreement with the most recent map of Mira-like stars reported in Deason et al (2016). The SMC's outer stellar density contours show a characteristic S-shape, symptomatic of the on-set of tidal stripping. Beyond several degrees from the center of the dwarf, the Gaia RR Lyrae stars trace the Cloud's trailing arm, extending towards the LMC. This stellar tidal tail mapped with RR Lyrae is not aligned with the gaseous Magellanic Bridge, and is shifted by some ~5 degrees from the Young Main Sequence bridge. We use the offset between the bridges to put constraints on the density of the hot gaseous corona of the Milky Way.

Ample observational capabilities exist today to detect the small density perturbations that low-mass dark matter subhaloes impart on stellar streams from disrupting Galactic satellites. In anticipation of these observations, we investigate the expected number and size of gaps by combining an analytic prescription for gap evolution on circular orbits with the flux of subhaloes near the stream. We explore the distribution of gap sizes and depths for a typical cold stream around the Milky Way and find that for a given stream age and gap depth, each subhalo mass produces a characteristic gap size. For a stream with an age of a few Gyr, orbiting at a distance of 10-20 kpc from the Galactic center, even modest subhaloes with a mass of $10^6-10^7 M_\odot$ produce gaps with sizes that are on the order of several degrees. We consider the number and distribution of gap sizes created by subhaloes with masses $10^5-10^9 M_\odot$, accounting for the expected depletion of subhaloes by the Milky Way disk, and present predictions for six cold streams around the Milky Way. For Pal 5, we forecast 0.7 gaps with a density depletion of at least 25\% and a typical gap size of $8^\circ$. Thus, there appears to be {\it no tension} between the recent non-detection of density depletions in the Pal 5 tidal tails and $\Lambda$CDM expectations. These predictions can be used to guide the scale of future gap searches.

In aspherical potentials orbital planes continuously evolve. The gravitational torques impel the angular momentum vector to precess, that is to slowly stray around the symmetry axis, and nutate, i.e. swing up and down periodically in the perpendicular direction. This familiar orbital pole motion - if detected and measured - can reveal the shape of the underlying gravitational potential, the quantity only crudely gauged in the Galaxy so far. Here we demonstrate that the debris poles of stellar tidal streams show a very similar straying and swinging behavior, and give analytic expressions to link the amplitude and the frequency of the pole evolution to the flattening of the dark matter distribution. While these results are derived for near-circular orbits, we show they are also valid for eccentric orbits. Most importantly, we explain how the differential orbital plane precession leads to the broadening of the stream and show that streams on polar orbits ought to scatter faster. We provide expressions for the stream width evolution as a function of the axisymmetric potential flattening and the angle from the symmetry plane and prove that our models are in good agreement with streams produced in N-body simulations. Interestingly, the same intuition applies to streams whose progenitors are on short or long-axis loops in a triaxial potential. Finally, we present a compilation of the Galactic cold stream data, and discuss how the simple picture developed here, along with stream modelling, can be used to constrain the symmetry axes and flattening of the Milky Way.

We establish the connection between the Magellanic Clouds (MCs) and the dwarf galaxy candidates discovered in the Dark Energy Survey (DES) by building a dynamical model of the MC satellite populations, based on an extensive suite of tailor-made numerical simulations. Our model takes into account the response of the Galaxy to the MCs infall, the dynamical friction experienced by the MCs and the disruption of the MC satellites by their hosts. The simulation suite samples over the uncertainties in the MC's proper motions, the masses of the MW and the Clouds themselves and allows for flexibility in the intrinsic volume density distribution of the MC satellites. As a result, we can accurately reproduce the DES satellites' observed positions and kinematics. Assuming that Milky Way (MW) dwarfs follow the distribution of subhaloes in $\Lambda$CDM, we further demonstrate that, of 14 observed satellites, the MW halo contributes fewer than 4 (8) of these with 68% (95%) confidence and that 7 (12) DES dwarfs have probabilities greater than 0.7 (0.5) of belonging to the LMC. Marginalising over the entire suite, we constrain the total number of the Magellanic satellites at ~70, the mass of the LMC around $10^{11}M_\odot$ and show that the Clouds have likely endured only one Galactic pericentric passage so far. Finally, we give predictions for the line-of-sight velocities and the proper motions of the satellites discovered in the vicinity of the LMC.

We report the discovery of a substantial stellar overdensity in the periphery of the Large Magellanic Cloud (LMC), found using public imaging from the first year of the Dark Energy Survey. The structure appears to emanate from the edge of the outer LMC disk at a radius $\approx 13.5$ degrees due north of its centre, and stretches more than $10$ kpc towards the east. It is roughly $1.5$ kpc wide and has an integrated $V$-band luminosity of at least $M_V = -7.4$. The stellar populations in the feature are indistinguishable from those in the outer LMC disk. We attempt to quantify the geometry of the outer disk using simple planar models, and find that only a disk with mild intrinsic ellipticity can simultaneously explain the observed stellar density on the sky and the azimuthal line-of-sight distance profile. We also see possible non-planar behaviour in the outer disk that may reflect a warp and/or flare, as well as deviations that resemble a ring-like structure between $\sim9-12$ degrees from the LMC centre. Based on all these observations, we conclude that our remote, stream-like feature is likely comprised of material that has been stripped from the outskirts of the LMC disk, although we cannot rule out that it represents a transient overdensity in the disk itself. We conduct a simple $N$-body simulation to show that either type of structure could plausibly arise due to the tidal force of the Milky Way; however we also recognize that a recent close interaction between the LMC and the SMC may be the source of the stripping or perturbation. Finally, we observe evidence for extremely diffuse LMC populations extending to radii of $\sim 18.5$ kpc in the disk plane ($\approx 20$ degrees on the sky), corroborating previous spectroscopic detections at comparable distances.

We develop a formalism for modelling the impact of dark matter subhaloes on cold thin streams. Our formalism models the formation of a gap in a stream in angle-frequency space and is able to handle general stream and impact geometry. We analyse an $N$-body simulation of a cold stream formed from a progenitor on an eccentric orbit in an axisymmetric potential, which is perturbed by a direct impact from a $10^8 M_\odot$ subhalo, and produce a complete generative model of the perturbed stream that matches the simulation well at a range of times. We show how the results in angle-frequency space can be related to physical properties of the gaps and that previous results for more constrained simulations are recovered. We demonstrate how our results are dependent upon the mass of the subhalo and the location of the impact along the stream. We find that gaps formed far downstream grow more rapidly than those closer to the progenitor due to the more ordered nature of the stream members far from the progenitor. Additionally, we show that the minimum gap density plateaus in time at a value that decreases with increasing subhalo mass.

The shape of the Galactic dark halo can, in principle, be inferred through modelling of stellar tidal streams in the Milky Way halo. The brightest and the longest of these, the Sagittarius stream, reaches out to large Galactocentric distances and hence can deliver the tightest constraints on the Galaxy's potential. In this contribution, we revisit the idea that the Sagittarius Stream was formed from a rotating progenitor. We demonstrate that the angle between the disk's angular momentum and the progenitor's orbital angular momentum does not remain constant throughout the disruption. Instead, it undergoes a dramatic evolution caused, in part, by the changes in the progenitor's moment of inertia tensor. We show that, even in a spherical potential, the streams produced as a result of a disky dwarf disruption appear to be "precessing". Yet, this debris plane evolution is illusory as it is solely caused by the swaying and wobbling of the progenitor's disk. Stream plane precession is therefore not an unambiguous indicator of asphericity of the dark halo.

In this work, I investigate the properties of Lyman limit systems (LLSs) using state-of-the-art zoom-in cosmological galaxy formation simulations with on the fly radiative transfer, which includes both the cosmic UV background (UVB) and local stellar sources. I compare the simulation results to observations of the incidence frequency of LLSs and the HI column density distribution function over the redshift range $z=2-5$ and find good agreement. I explore the connection between LLSs and their host halos and find that LLSs reside in halos with a wide range of halo masses with a nearly constant covering fraction within a virial radius. Over the range $z=2-5$, I find that more than half of the LLSs reside in halos with $M < 10^{10}h^{-1}M_\odot$, indicating that absorption line studies of LLSs can probe these low-mass galaxies which H$_2$-based star formation models predict to have very little star formation. I study the physical state of individual LLSs and test a simple model (Schaye 2001) which encapsulates many of their properties. I confirm that LLSs have a characteristic absorption length given by the Jeans length and that they are in photoionization equilibrium at low column densities. Finally, I investigate the self-shielding of LLSs to the UVB and explore how the non-sphericity of LLSs affects the photoionization rate at a given $N_{\rm HI}$. I find that at $z\approx 3$, LLSs have an optical depth of unity at a column density of $\sim 10^{18} {\rm cm}^{-2}$ and that this is the column density which characterizes the onset of self-shielding.

The dynamics of stellar streams in rotating barred potentials is explained for the first time. Naturally, neighbouring stream stars reach pericentre at slightly different times. In the presence of a rotating bar, these neighbouring stream stars experience different bar orientations during pericentric passage and hence each star receives a different torque from the bar. These differing torques reshape the angular momentum and energy distribution of stars in the stream, which in turn changes the growth rate of the stream. For a progenitor orbiting in the same sense as the bar's rotation and satisfying a resonance condition, the resultant stream can be substantially shorter or longer than expected, depending on whether the pericentric passages of the progenitor occur along the bar's minor or major axis respectively. We present a full discussion of this phenomenon focusing mainly on streams confined to the Galactic plane. In stark contrast with the evolution in static potentials, which give rise to streams that grow steadily in time, rotating barred potentials can produce dynamically old, short streams. This challenges the traditional viewpoint that the inner halo consists of well phase-mixed material whilst the tidally-disrupted structures in the outer halo are more spatially coherent. We argue that this mechanism plays an important role in explaining the mysteriously short Ophiuchus stream that was recently discovered near the bulge region of the Milky Way.

This paper explores the effect of the LMC on the mass estimates obtained from the timing argument. We show that accounting for the presence of the LMC systematically lowers the Local Group mass ($M_{\rm LG}$) derived from the relative motion of the Milky Way--Andromeda pair. Motivated by this result we apply a Bayesian technique devised by Pe\~narrubia et al. (2014) to simultaneously fit (i) distances and velocities of galaxies within 3~Mpc and (ii) the relative motion between the Milky Way and Andromeda derived from HST observations, with the LMC mass ($M_{\rm LMC}$) as a free parameter. Our analysis returns a Local Group mass $M_{\rm LG}=2.64^{+0.42}_{-0.38}\times 10^{12}M_\odot$ at a 68\% confidence level. The masses of the Milky Way, $M_{\rm MW}=1.04_{-0.23}^{+0.26}\times 10^{12}M_\odot$, and Andromeda, $M_{\rm M31}=1.33_{-0.33}^{+0.39}\times 10^{12}M_\odot$, are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find a (total) LMC mass $M_{\rm LMC}=0.25_{-0.08}^{+0.09}\times 10^{12}M_\odot$, which is indicative of an extended dark matter halo and supports the scenario where this galaxy is just past its first pericentric approach. Consequently, these results suggest that the LMC may induce significant perturbations on the Galactic potential.

Cold or Warm, the Dark Matter substructure spectrum must extend to objects with masses as low as $10^7 M_\odot$, according to the most recent Lyman-$\alpha$ measurements. Around a Milky Way-like galaxy, more than a thousand of these subhaloes will not be able to form stars but are dense enough to survive even deep down in the potential well of their host. There, within the stellar halo, these dark pellets will bombard tidal streams as they travel around the Galaxy, causing small but recognizable damage to the stream density distribution. The detection and characterization of these stream ruptures will allow us to constrain the details of the subhalo-stream interaction. In this work, for the first time, we will demonstrate how the properties of a subhalo, most importantly its mass and size, can be reliably inferred from the gap it produces in a tidal stream. For a range of realistic observational setups, mimicking e.g. SDSS, DES, Gaia and LSST data, we find that it is possible to measure the {\it complete set} of properties (including the phase-space coordinates during the flyby) of dark perturbers with $M>10^7 M_\odot$, up to a 1d degeneracy between the mass and velocity.

There is hope to discover dark matter subhalos free of stars (predicted by the current theory of structure formation) by observing gaps they produce in tidal streams. In fact, this is the most promising technique for dark substructure detection and characterization as such gaps grow with time, magnifying small perturbations into clear signatures observable by ongoing and planned Galaxy surveys. To facilitate such future inference, we develop a comprehensive framework for studies of the growth of the stream density perturbations. Starting with simple assumptions and restricting to streams on circular orbits, we derive analytic formulae that describe the evolution of all gap properties (size, density contrast etc) at all times. We uncover complex, previously unnoticed behavior, with the stream initially forming a density enhancement near the subhalo impact point. Shortly after, a gap forms due to the relative change in period induced by the subhalo's passage. There is an intermediate regime where the gap grows linearly in time. At late times, the particles in the stream overtake each other, forming caustics, and the gap grows like $\sqrt{t}$. In addition to the secular growth, we find that the gap oscillates as it grows due to epicyclic motion. We compare this analytic model to N-body simulations and find an impressive level of agreement. Importantly, when analyzing the observation of a single gap we find a large degeneracy between the subhalo mass, the impact geometry and kinematics, the host potential and the time since flyby.

We study the high column density regime of the HI column density distribution function and argue that there are two distinct features: a turnover at NHI ~ 10^21 cm^-2 which is present at both z=0 and z ~ 3, and a lack of systems above NHI ~ 10^22 cm^-2 at z=0. Using observations of the column density distribution, we argue that the HI-H2 transition does not cause the turnover at NHI ~ 10^21 cm^-2, but can plausibly explain the turnover at NHI > 10^22 cm^-2. We compute the HI column density distribution of individual galaxies in the THINGS sample and show that the turnover column density depends only weakly on metallicity. Furthermore, we show that the column density distribution of galaxies, corrected for inclination, is insensitive to the resolution of the HI map or to averaging in radial shells. Our results indicate that the similarity of HI column density distributions at z=3 and z=0 is due to the similarity of the maximum HI surface densities of high-z and low-z disks, set presumably by universal processes that shape properties of the gaseous disks of galaxies. Using fully cosmological simulations, we explore other candidate physical mechanisms that could produce a turnover in the column density distribution. We show that while turbulence within GMCs cannot affect the DLA column density distribution, stellar feedback can affect it significantly if the feedback is sufficiently effective in removing gas from the central 2-3 kpc of high-redshift galaxies. Finally, we argue that it is meaningful to compare column densities averaged over ~ kpc scales with those estimated from quasar spectra which probe sub-pc scales due to the steep power spectrum of HI column density fluctuations observed in nearby galaxies.