Proper Motions, Orbits, and Tidal Influences of Milky Way Dwarf Spheroidal Galaxies

We combine Gaia EDR3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we Monte Carlo sample over the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare the dSph radial orbital phase to subhalos in MW-like N-body simulations and infer that there is not an excess of satellites near their pericenter. With projections of future Gaia data releases, we find dSph orbital precision will be limited by uncertainties in the distance and/or MW potential rather than proper motion precision. Finally, we provide our membership catalogs to enable community follow-up.Comment: 26 pages, 13 figures + appendix with extra figures. ApJ accepted. Catalogs with membership, additional figures, and a machine readable compilation of tables 1-4 are available at https://zenodo.org/record/653329

Under Pressure: Quenching Star Formation in Low-Mass Satellite Galaxies via Stripping

Recent studies of galaxies in the local Universe, including those in the Local Group, find that the efficiency of environmental (or satellite) quenching increases dramatically at satellite stellar masses below ~ $10^8\ {\rm M}_{\odot}$. This suggests a physical scale where quenching transitions from a slow "starvation" mode to a rapid "stripping" mode at low masses. We investigate the plausibility of this scenario using observed HI surface density profiles for a sample of 66 nearby galaxies as inputs to analytic calculations of ram-pressure and viscous stripping. Across a broad range of host properties, we find that stripping becomes increasingly effective at $M_{*} < 10^{8-9}\ {\rm M}_{\odot}$, reproducing the critical mass scale observed. However, for canonical values of the circumgalactic medium density ($n_{\rm halo} < 10^{-3.5}$${\rm cm}^{-3}$), we find that stripping is not fully effective; infalling satellites are, on average, stripped of < 40 - 70% of their cold gas reservoir, which is insufficient to match observations. By including a host halo gas distribution that is clumpy and therefore contains regions of higher density, we are able to reproduce the observed HI gas fractions (and thus the high quenched fraction and short quenching timescale) of Local Group satellites, suggesting that a host halo with clumpy gas may be crucial for quenching low-mass systems in Local Group-like (and more massive) host halos.Comment: updated version after review, now accepted to MNRAS; Accepted 2016 August 22. Received 2016 August 18; in original form 2016 June 2

Multiple Chemodynamic Stellar Populations of the Ursa Minor Dwarf Spheroidal Galaxy

We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $\sigma_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($\sigma_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($\epsilon=0.75\pm0.03$) and the metal-poor population is closer to spherical ($\epsilon=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.Comment: 20 pages, 11 figures, data included. Comments welcome. Accepted to MNRA

Robust velocity dispersion and binary population modeling of the ultra-faint dwarf galaxy Reticulum II

We apply a Bayesian method to model multi-epoch radial velocity measurements in the ultra-faint dwarf galaxy Reticulum II, fully accounting for the effects of binary orbital motion and systematic offsets between different spectroscopic datasets. We find that the binary fraction of Ret II is higher than 0.5 at the 90% confidence level, if the mean orbital period is assumed to be 30 years or longer. Despite this high binary fraction, we infer a best-fit intrinsic dispersion of 2.8$_{-1.2}^{+0.7}$ km/s, which is smaller than previous estimates, but still indicates Ret II is a dark-matter dominated galaxy. We likewise infer a $\lesssim$ 1% probability that Ret II's dispersion is due to binaries rather than dark matter, corresponding to the regime $M_{\odot}/L_{\odot} \lesssim$ 2. Our inference of a high close binary fraction in Ret II echoes previous results for the Segue 1 ultra-faint dwarf and is consistent with studies of Milky Way halo stars that indicate a high close binary fraction tends to exist in metal-poor environments.Comment: 9 pages, 5 figures, to be submitted to MNRA