A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and are available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectroscopically observed in the literature. We derive systemic PMs for 39 galaxies and galaxy candidates out to 420 kpc, and generally find good consistency for the subset with measurements available from other studies. We derive the implied Galactocentric velocities, and calculate orbits in canonical MW halo potentials of low (0.8 × 10 12 M ) and high mass (1.6 × 10 12 M ). Comparison of the distributions of orbital apocenters and 3D velocities to the halo virial radius and escape velocity, respectively, suggests that the satellite kinematics are best explained in the high-mass halo. Tuc III, Crater II, and additional candidates have orbital pericenters small enough to imply significant tidal influences. Relevant to the missing satellite problem, the fact that fewer galaxies are observed to be near apocenter than near pericenter implies that there must be a population of distant dwarf galaxies yet to be discovered. Of the 39 dwarfs: 12 have orbital poles that do not align with the MW plane of satellites (given reasonable assumptions about its intrinsic thickness); 10 have insufficient PM accuracy to establish whether they align; and 17 satellites align, of which 11 are co-orbiting and (somewhat surprisingly, in view of prior knowledge) 6 are counter-orbiting. Group infall might have contributed to this, but no definitive association is found for the members of the Crater-Leo group.

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 present a catalog of Fe, Mg, Si, Ca, and Ti abundances for 2961 red giant stars that are likely members of eight dwarf satellite galaxies of the Milky Way (MW): Sculptor, Fornax, Leo I, Sextans, Leo II, Canes Venatici I, Ursa Minor, and Draco. For the purposes of validating our measurements, we also observed 445 red giants in MW globular clusters and 21 field red giants in the MW halo. The measurements are based on Keck/DEIMOS medium-resolution spectroscopy combined with spectral synthesis. We estimate uncertainties in [Fe/H] by quantifying the dispersion of [Fe/H] measurements in a sample of stars in monometallic globular clusters. We estimate uncertainties in Mg, Si, Ca, and Ti abundances by comparing our medium-resolution spectroscopic measurements to high-resolution spectroscopic abundances of the same stars. For this purpose, our DEIMOS sample included 132 red giants with published high-resolution spectroscopy in globular clusters, the MW halo field, and dwarf galaxies. The standard deviations of the differences in [Fe/H] and [α/Fe] (the average of [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe]) between the two samples is 0.15 and 0.16, respectively. This catalog represents the largest sample of multi-element abundances in dwarf galaxies to date. The next papers in this series draw conclusions on the chemical evolution, gas dynamics, and star formation histories from the catalog presented here. The wide range of dwarf galaxy luminosity reveals the dependence of dwarf galaxy chemical evolution on galaxy stellar mass.

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M vir,MW ). Despite Leo I's large Galacto-centric space velocity (200 km s −1 ) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M vir,MW > 10 12 M at 95% confidence for a variety of Bayesian priors on M vir,MW . In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M vir,MW would increase by 30%. Imposing a mass weighted ΛCDM prior, we find a median Milky Way virial mass of M vir,MW = 1.6 × 10 12 M , with a 90% confidence interval of [1.0 − 2.4] × 10 12 M . We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

According to LCDM theory, hierarchical evolution occurs on all mass scales, implying that satellites of the Milky Way should also have companions. The recent discovery of ultra-faint dwarf galaxy candidates in close proximity to the Magellanic Clouds provides an opportunity to test this theory. We present proper motion (PM) measurements for 13 of the 32 new dwarf galaxy candidates using Gaia data release 2. All 13 also have radial velocity measurements. We compare the measured 3D velocities of these dwarfs to those expected at the corresponding distance and location for the debris of an LMC analog in a cosmological numerical simulation. We conclude that 4 of these galaxies (Hor1, Car2, Car3 and Hyi1) have come in with the Magellanic Clouds, constituting the first confirmation of the type of satellite infall predicted by LCDM. Ret2, Tuc2 and Gru1 have velocity components that are not consistent within 3 sigma of our predictions and are therefore less favorable. Hya2 and Dra2 could be associated with the LMC and merit further attention. We rule out Tuc3, Cra2, Tri2 and Aqu2 as potential members. Of the dwarfs without measured PMs, 5 of them are deemed unlikely on the basis of their positions and distances alone as being too far from the orbital plane expected for LMC debris (Eri2, Ind2, Cet2, Cet3 and Vir1). For the remaining sample, we use the simulation to predict PMs and radial velocities, finding that Phx2 has an overdensity of stars in DR2 consistent with this PM prediction.