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 < β < 0.4. However, for stars with [Fe/H]> −1.7 we measure extreme values of β ∼ 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 10 10 M ⊙ 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.

In a companion paper by Koposov et al., RR Lyrae from 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 motiontrack 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 × 10 11 M 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 × 10 11 M . 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 perform a search for stellar streams around the Milky Way using the first 3 yr of multiband optical imaging data from the Dark Energy Survey (DES). We use DES data covering ∼5000 deg 2 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 ∼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 11 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 extratidal 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, and the large-and small-scale distribution of dark matter around the Milky Way.

Tidal streams in the Milky Way are sensitive probes of the population of low-mass 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 streams based on the action-angle representation of streams. The heart of this calculus is a lineof-parallel-angle approach that calculates the perturbed distribution function of a stream segment by undoing the effect of all relevant 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 M , accounting for the stream's internal dispersion and overlapping impacts. We study the statistical properties of density and track fluctuations with large suites of simulations of the effect of subhalo fly-bys. 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. We also find significant density and track bispectra that are observationally accessible. We further demonstrate that different projections of the track all reflect the same pattern of perturbations, facilitating their observational measurement. We apply this formalism to data for the Pal 5 stream and make a first rigorous determination of 10 +11 −6 dark-matter subhalos with masses between 10 6.5 M and 10 9 M within 20 kpc from the Galactic center (corresponding to 1.4 +1.6 −0.9 times the number predicted by CDMonly simulations or to f sub (r < 20 kpc) ≈ 0.2 %) assuming that the Pal 5 stream is 5 Gyr old. Improved data will allow measurements of the subhalo mass function down to 10 5 M , thus definitively testing whether dark matter is clumpy on the smallest scales relevant for galaxy formation.

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

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 LG ) derived from the relative motion of the Milky Way-Andromeda pair. Motivated by this result we apply a Bayesian technique devised by Peñarrubia 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 LMC ) as a free parameter. Our analysis returns a Local Group mass M LG = 2.64 +0.42 −0.38 × 10 12 M at a 68% confidence level. The masses of the Milky Way, M MW = 1.04 +0.26 −0.23 × 10 12 M , and Andromeda, M M31 = 1.33 +0.39 −0.33 × 10 12 M , are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find a (total) LMC mass M LMC = 0.25 +0.09 −0.08 × 10 12 M , 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.

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 < 10 8 galaxies, allowing us to bound the peak halo mass of the faintest MW satellite to M vir > 2.4 × 10 8 M (1σ). Secondly, by translating to a Warm Dark Matter (WDM) cosmology, we bound the thermal relic mass m WDM > 2.9 keV at 95% confidence, on a par with recent constraints from the Lyman-α 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 8 M 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.

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 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 , 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 • . Thus, there appears to be no tension between the recent non-detection of density depletions in the Pal 5 tidal tails and ΛCDM expectations. These predictions can be used to guide the scale of future gap searches.