Search

The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history including the recent pericentric passage of the LMC, and recent work has found that the MW disk is moving with a lower bound "travel velocity" of $\sim 32~{\rm km}~{\rm s}^{-1}$ with respect to the outer stellar halo. Previous Timing Argument measurements attempt to account for this non-equilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10--12\% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume., Comment: 14 pages, 3 figures; Published in ApJ; updated to match published version

We present analysis of the proper-motion (PM) field of the red clump stars in the Large Magellanic Cloud (LMC) disk using the Gaia Early Data Release 3 catalog. Using a kinematic model based on old stars with 3D velocity measurements, we construct the residual PM field by subtracting the center-of-mass motion and internal rotation motion components. The residual PM field reveals asymmetric patterns, including larger residual PMs in the southern disk. Comparisons between the observed residual PM field with those of five numerical simulations of an LMC analog that is subject to the tidal fields of the Milky Way and the Small Magellanic Cloud (SMC) show that the present-day LMC is not in dynamical equilibrium. We find that both the observed level of disk heating (PM residual root-mean-square of 0.057$\pm$0.002 mas yr$^{-1}$) and kinematic asymmetry are not reproduced by Milky Way tides or if the SMC impact parameter is larger than the size of the LMC disk. This measured level of disk heating provides a novel and important method to validate numerical simulations of the LMC-SMC interaction history. Our results alone put constraints on an impact parameter $\lesssim$10 kpc and impact timing $, Comment: 17 pages, 10 figures, resubmitted to AAS Journal

We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system., 32 pages, 15 figures, accepted for publication on ApJ

Gravitational interactions between the Large Magellanic Cloud (LMC) and the stellar and dark matter halo of the Milky Way are expected to give rise to disequilibrium phenomena in the outer Milky Way. A local wake is predicted to trail the orbit of the LMC, while a large-scale over-density is predicted to exist across a large area of the northern Galactic hemisphere. Here we present the detection of both the local wake and Northern over-density (hereafter the "collective response") in an all-sky star map of the Galaxy based on 1301 stars at 60, Comment: To appear in the 22 April 2021 issue of Nature

Recent advancements in astrometry and in cosmological models of dark matter halo growth have significantly changed our understanding of the dynamics of the Local Group. The most dramatic changes owe to a new picture of the structure and dynamics of the Milky Way’s most massive satellite galaxy, the Large Magellanic Cloud (LMC), which is most likely on its first passage about the Milky Way and ten times larger in mass than previously assumed. The LMC’s orbit through the Milky Way’s dark matter and stellar halo will leave characteristic signatures in both density and kinematics. Furthermore, the gravitational perturbations produced by both direct tidal forcing from the LMC and the response of the halo to its passage will together cause significant perturbations to the orbits of tracers of the Milky Way’s dark matter distribution. We advocate for the use of basis field expansion methods to fully capture and quantify these effects.

We present deep Hubble Space Telescope (HST) photometry of the ultra-faint dwarf (UFD) galaxies Pegasus III (Peg III) and Pisces II (Psc II), two of the most distant satellites in the halo of the Milky Way (MW). We measure the structure of both galaxies, derive mass-to-light ratios with newly determined absolute magnitudes, and compare our findings to expectations from UFD-mass simulations. For Peg III, we find an elliptical half-light radius of $a_h{{=}}1.88^{+0.42}_{-0.33}$ arcminutes ($118^{+31}_{-30}$ pc) and $M_V{=}{-4.17}^{+0.19}_{-0.22}$; for Psc II, we measure $a_h{=}1.31^{+0.10}_{-0.09}$ arcminutes ($69\pm8$ pc) and $M_V{=}{-4.28}^{+0.19}_{-0.16}$. We do not find any morphological features that indicate a significant interaction between the two has occurred, despite their close separation of only $\sim$40 kpc. Using proper motions (PMs) from Gaia early Data Release 3, we investigate the possibility of any past association by integrating orbits for the two UFDs in a MW-only and a combined MW and Large Magellanic Cloud (LMC) potential. We find that including the gravitational influence of the LMC is crucial, even for these outer-halo satellites, and that a possible orbital history exists where Peg III and Psc II experienced a close ($\sim$10-20 kpc) passage about each other just over $\sim$1 Gyr ago, followed by a collective passage around the LMC ($\sim$30-60 kpc) just under $\sim$1 Gyr ago. Considering the large uncertainties on the PMs and the restrictive priors imposed to derive them, improved PM measurements for Peg III and Psc II will be necessary to clarify their relationship. This would add to the rare findings of confirmed pairs of satellites within the Local Group., Comment: 18 pages, 11 figures; submitted to ApJ

Full list of authors: Nidever, David L.; Olsen, Knut; Choi, Yumi; Ruiz-Lara, Tomas; Miller, Amy E.; Johnson, L. Clifton; Bell, Cameron P. M.; Blum, Robert D.; Cioni, Maria-Rosa L.; Gallart, Carme; Majewski, Steven R.; Martin, Nicolas F.; Massana, Pol; Monachesi, Antonela; Noël, Noelia E. D.; Sakowska, Joanna D.; van der Marel, Roeland P.; Walker, Alistair R.; Zaritsky, Dennis; Bell, Eric F. Conn, Blair C.; de Boer, Thomas J. L.; Gruendl, Robert A.; Monelli, Matteo; Muñoz, Ricardo R.; Saha, Abhijit; Vivas, A. Katherina; Bernard, Edouard; Besla, Gurtina; Carballo-Bello, Julio A.; Dorta, Antonio; Martinez-Delgado, David; Goater, Alex; Rusakov, Vadim; Stringfellow, Guy S., We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances. © 2021. The American Astronomical Society. All rights reserved, Y.C., E.F.B., and A.M. acknowledge support from NSF grant AST 1655677. A.D., C.G., T.R.L., and M.M. acknowledge support by the Spanish Ministry of Economy and Competitiveness (MINECO) under the grants AYA2014-56795-P and AYA2017-89076-P as well as AYA2016-77237-C3-1-P. T.R.L. has support from a Spinoza grant (NWO) awarded to A. Helmi and acknowledges support by an MCIU Juan de la CiervaFormacion grant (FJCI-2016-30342). C.P.M.B. and M.-R.L.C. acknowledge support from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement No. 682115). A.M. acknowledges support from FONDECYT Regular 1181797. R. R.M. acknowledges partial support from project BASAL AFB170002 as well as FONDECYT project No. 1170364. D.M.D. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709). Image processing: Travis Rector (University of Alaska Anchorage), Mahdi Zamani, and Davide de Martin. Based on observations at Cerro Tololo Inter-American Observatory, NSF's National Optical-Infrared Astronomy Research Laboratory (NOAO Prop. ID: 2013A-0411 and 2013B-0440; PI: D. L. Nidever), which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey (DES) collaboration. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundacao Carlos Chagas Filho de Amparo, Financiadora de Estudos e Projetos, Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and the Ministerio da Ciencia, Tecnologia e Inovacao, the Deutsche Forschungsgemeinschaft, and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciencies de l'Espai (IEEC/CSIC), the Institut de Fisica d'Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universitat Munchen and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, the Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University. This work has made use of data from the European Space Agency (ESA) mission Gai a (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC,.https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

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., 25 pages, 13 figures. Accepted for publication in ApJ

CEMP-no stars, a subset of carbon enhanced metal poor (CEMP) stars ($\rm [C/Fe]\geq0.7$ and $\rm [Fe/H]\lesssim-1$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $M_{\rm vir} \sim 10^8$ Msun and $M_{\ast}\sim10^3-10^4$ Msun at $z=0$, as well as in the halo of the Milky Way (MW). These CEMP-no stars are local fossils that may reflect the properties of the first (Pop~III) and second (Pop~II) generation of stars. However, cosmological simulations have struggled to reproduce the observed level of carbon enhancement of the known CEMP-no stars. Here we present new cosmological hydrodynamic zoom-in simulations of isolated UFDs that achieve a gas mass resolution of $m_{\rm gas}\sim60$ Msun. We include enrichment from Pop~III faint supernovae (SNe), with $ E_{\rm SN}=0.6\times10^{51}$ erg, to understand the origin of CEMP-no stars. We confirm that Pop~III and Pop~II stars are mainly responsible for the formation of CEMP and C-normal stars respectively. New to this study, we find that a majority of CEMP-no stars in the observed UFDs and the MW halo can be explained by Pop~III SNe with normal explosion energy ($ E_{\rm SN}=1.2\times10^{51}$~erg) and Pop~II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $\rm [C/Fe]\gtrsim2$, corresponding to the absolute carbon abundance of $\rm A(C)\gtrsim6.0$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $\rm [C/Fe]\approx3-4$, by adopting faint SNe, it is still challenging to reproduce CEMP-no stars with extreme level of carbon abundance of $\rm A(C)\approx 7.0-7.5$, observed both in the MW halo and UFDs., Comment: 14 pages, 8 figures, submitted to MNRAS

Indications of disequilibrium throughout the Milky Way (MW) highlight the need for compact,flexible, non-parametric descriptions of phase--space distributions of galaxies. We present a new representation of the current Dark Matter (DM) distribution and potential derived from N-body simulations of the Milky Way and Large Magellanic Cloud (LMC) system using Basis Function Expansions (BFEs). We incorporate methods to maximize the physical signal in the representation. As a result, the simulations of $10^8$ DM particles representing the MW--LMC system can be described by 354 coefficients. We find that the LMC induces asymmetric perturbations (odd l, m) to the MW's halo, which are not well-described by oblate, prolate, or triaxial halos. Furthermore, the energy in high-order even modes (l,m $\geq$ 2) is similar to average triaxial halos found in cosmological simulations. As such, the response of the MW's halo to the LMC must be accounted for in order to recover the imprints of its assembly history. The LMC causes the outer halo ($\geq$ 30 kpc) to shift from the disk center of mass (COM) by $\sim$15-25 kpc at present day, manifesting as a dipole in the BFE and in the radial velocities of halo stars. The shift depends on the LMC's infall mass, the distortion of the LMC's halo and the MW halo response. Within 30 kpc, halo tracers are expected to orbit the COM of the MW's disk, regardless of LMC infall mass. The LMC's halo is also distorted by MW tides, we discuss the implications for its mass loss and the subsequent effects on current Magellanic satellites., Comment: 35 pages, 25 figures, Accepted for publication in ApJ. 3d movie can be found here: https://vimeo.com/546207117

With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: i.) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); ii.) Magellanic satellites that were recently captured by the MCs $, Comment: 32 pages, 11 figures. Accepted to ApJ

We present resolved HI synthesis maps from the Very Large Array (VLA) of three interacting dwarf systems: the NGC 3664 dwarf pair, the NGC 3264 dwarf pair, and the UGC 4638 dwarf triplet. All three dwarf systems are captured at various stages of interaction and span a range of environments. We detect clear hallmarks of tidal interactions through the presence of HI bridges, and diffuse HI extensions that surround the dwarfs. We overlay the HI data on Pan-STARRS r-band images and find further evidence of tidal interactions through coincident distorted HI and tidal stellar features in NGC 3264 and UGC 4638, and an unwound spiral arm pointing towards its smaller companion in NGC 3264. In UGC 4638, both the gas and diffuse stars are extended to similar radii east of the primary, which could indicate that the smaller dwarf in the system has already completed one pass through the primary. We additionally find that our three systems, and those from the Local Volume TiNy Titans survey, are not HI deficient and thus the interaction has not resulted in a loss of gas from the systems. A comparison with non-interacting dwarf galaxies shows that the interactions have a significant impact on the kinematics of the systems. Our new resolved HI kinematics, combined with detailed stellar and HI morphologies, provide crucial constraints for future dynamical modelling of hierarchical mergers and the baryon cycle at the low-mass scale., Comment: 17 pages, 6 figures, accepted to AJ

We present N-body simulations of a Sagittarius-like dwarf spheroidal galaxy (Sgr) that follow its orbit about the Milky Way (MW) since its first crossing of the Galaxy's virial radius to the present day. As Sgr orbits around the MW, it excites vertical oscillations, corrugating and flaring the Galactic stellar disc. These responses can be understood by a two-phase picture in which the interaction is first dominated by torques from the wake excited by Sgr in the MW dark halo before transitioning to tides from Sgr's direct impact on the disc at late times. We show for the first time that a massive Sgr model simultaneously reproduces the locations and motions of arc-like over densities, such as the Monoceros Ring and the Triangulum Andromeda stellar clouds, that have been observed at the extremities of the disc, while also satisfying the solar neighbourhood constraints on the vertical structure and streaming motions of the disc. In additional simulations, we include the Large Magellanic Cloud (LMC) self consistently with Sgr. The LMC introduces coupling through constructive and destructive interference, but no new corrugations. In our models, the excitation of the current structure of the outer disk can be traced to interactions as far back as 6-7 Gyr ago (corresponding to $z\leq1$). Given the apparently quiescent accretion history of the MW over this timescale, this places Sgr as the main culprit behind the vertical oscillations of the disc and the last major accretion event for the Galaxy with the capacity to modulate its chemodynamical structure.

The globular cluster (GC) system of the Milky Way (MW) provides important information on the MW’s present structure and past evolution. Full 3d motions, accessed through proper motions (PMs), are required to calculate accurate orbits of GCs in the MW halo. We present our HST program to create a PM database for 20 halo GCs. We demonstrate how the observed PMs of individual GCs can be used to study their origins, and we also describe how the PM measurements of our entire targets can be used to constrain the anisotropy profile. Finally, we describe how our PM results can be used for Gaia as an external check, and discuss prospects of PM measurements with HST and Gaia in the coming years.

Space observatories like the Hubble Space Telescope and Gaia are providing unprecedented 6D phase space information of satellite galaxies. Such measurements can shed light on the structure and assembly history of the Local Group, but improved statistical methods are needed to use them efficiently. Here we illustrate such a method using analogs of the Local Group's two most massive satellite galaxies, the Large Magellanic Cloud (LMC) and Triangulum (M33), from the Illustris dark-matter-only cosmological simulation. We use a Bayesian inference scheme combining measurements of positions, velocities, and specific orbital angular momenta ( j ) of the LMC/M33 with importance sampling of their simulated analogs to compute posterior estimates of the Milky Way (MW) and Andromeda's (M31) halo masses. We conclude the resulting host halo mass is more susceptible to bias when using measurements of the current position and velocity of satellites, especially when satellites are at short-lived phases of their orbits (i.e. at pericentre). Instead, the j value of a satellite is well-conserved over time and provides a more reliable constraint on host mass. The inferred virial mass of the MW (M31) using j of the LMC (M33) is $\rm M_{vir, MW}=1.02^{+0.77}_{-0.55}\times10^{12}\; M_{\odot}$ ($\rm M_{vir, M31}=1.37^{+1.39}_{-0.75}\times10^{12}\; M_{\odot}$). Choosing simulated analogs whose j values are consistent with the conventional picture of a previous (< 3 Gyr ago), close encounter (< 100 kpc) of M33 about M31 results in a very low virial mass for M31 ($\rm\sim\!10^{12}\; M_{\odot}$). This supports the new scenario put forth in Patel et al. (2017), wherein M33 is on its first passage about M31 or on a long period orbit. We conclude that this Bayesian inference scheme, utilising satellite j, is a promising method to reduce the current factor of two spread in the mass range of the MW and M31 moving forward., 24 pages, 9 figures. Submitted to MNRAS

High precision proper motion (PM) measurements are available for approximately 20% of all known dwarf satellite galaxies of the Milky Way (MW). Here we extend the Bayesian framework of Patel et al. (2017b) to include all MW satellites with measured 6D phase space information and apply it with the Illustris-Dark simulation to constrain the MW's mass. Using the properties of each MW satellite individually, we find that the scatter among mass estimates is reduced when the magnitude of specific orbital angular momentum (j) is adopted rather than their combined instantaneous positions and velocities. We also find that high j satellites (i.e. Leo II) constrain the upper limits for the MW's mass and low j satellites rather than the highest speed satellites (i.e. Leo I and LMC), set the lower mass limits. When j of all classical satellites is used to simultaneously estimate the MW's mass, we conclude the halo mass is $0.85^{+0.23}_{-0.26}\times 10^{12} \rm \, M_{\odot}$ (including Sagittarius dSph) and $0.96^{+0.29}_{-0.28} \times 10^{12} \rm \, M_{\odot}$ (excluding Sagittarius dSph), cautioning that low j satellites on decaying orbits like Sagittarius dSph may bias the distribution. These estimates markedly reduce the current factor of two spread in the mass range of the MW. We also find a well-defined relationship between host halo mass and satellite j distribution, which yields the prediction that upcoming PMs for ultra-faint dwarfs should reveal j within $ 5\times10^3 - 10^4$ kpc km s$^{-1}$. This is a promising method to significantly constrain the cosmologically expected mass range for the MW and eventually M31 as more satellite PMs become available., 20 pages, 7 figures, submitted to ApJ

Using N-body simulations of the Large Magellanic Cloud (LMC's) passage through the Milky Way (MW), tailored to reproduce observed kinematic properties of both galaxies, we show that the high-speed tail of the Solar Neighborhood dark matter distribution is overwhelmingly of LMC origin. Two populations contribute at high speeds: 1) Particles that were once bound to the LMC, and 2) MW halo particles that have been accelerated owing to the response of the halo to the recent passage of the LMC. These particles reach speeds of 700-900 km/s with respect to the Earth, above the local escape speed of the MW. The high-speed particles follow trajectories similar to the Solar reflex motion, with peak velocities reached in June. For low-mass dark matter, these high-speed particles can dominate the signal in direct-detection experiments, extending the reach of the experiments to lower mass and elastic scattering cross sections even with existing data sets. Our study shows that even non-disrupted MW satellite galaxies can leave a significant dark-matter footprint in the Solar Neighborhood., Comment: Submitted to JCAP, comments are welcome

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material., Comment: 9 pages, 5 figues, 1 table, submitted to ApJ

We report the discovery of 15 stars in the H3 survey that lie, in projection, near the tip of the trailing gaseous Magellanic Stream (MS). The stars have Galactocentric velocities $< -155$ km s$^{-1}$, Galactocentric distances of $\approx 40$ to 80 kpc (increasing along the MS), and [Fe/H] consistent with that of stars in the Small Magellanic Cloud. These 15 stars comprise 94% (15 of 16) of the H3 observed stars to date that have $R_{GAL} > 37.5$ kpc, $-$350 km s$^{-1} < V_{GSR} < -155$ km s$^{-1}$, and are not associated with the Sagittarius Stream. They represent a unique portion of the Milky Way's outer halo phase space distribution function and confirm that unrelaxed structure is detectable even at radii where H3 includes only a few hundred stars. Due to their statistical excess, their close association with the MS and H I compact clouds in the same region, both in position and velocity space, and their plausible correspondence with tidal debris in a published simulation, we identify these stars as debris of past Magellanic Cloud encounters. These stars are evidence for a stellar component of the tidal debris field far from the Clouds themselves and provide unique constraints on the interaction., Accepted for publication in ApJL

We explore the stellar structure of the Large Magellanic Cloud (LMC) disk using data from the Survey of the MAgellanic Stellar History (SMASH) and the Dark Energy Survey. We detect a ring-like stellar overdensity in the red clump star count map at a radius of ~6 degrees (~5.2 kpc at the LMC distance) that is continuous over ~270 degrees in position angle and is only limited by the current data coverage. The overdensity shows an amplitude up to 2.5 times higher than that of the underlying smooth disk. This structure might be related to the multiple arms found by de Vaucouleurs. We find that the overdensity shows spatial correlation with intermediate-age star clusters, but not with young (< 1 Gyr) main-sequence stars, indicating the stellar populations associated with the overdensity are intermediate in age or older. Our findings on the LMC overdensity can be explained by either of two distinct formation mechanisms of a ring-like overdensity: (1) the overdensity formed out of an asymmetric one-armed spiral wrapping around the LMC main body, which is induced by repeated encounters with the Small Magellanic Cloud (SMC) over the last Gyr, or (2) the overdensity formed very recently as a tidal response to a direct collision with the SMC. Although the measured properties of the overdensity alone cannot distinguish between the two candidate scenarios, the consistency with both scenarios suggests that the ring-like overdensity is likely a product of tidal interaction with the SMC, but not with the Milky Way halo., 15 pages, 4 figures, 1 table, Published in ApJ

Discoveries of low mass galaxy pairs and groups are increasing. Studies indicate that dwarf galaxy pairs are gas rich in the field and exhibit elevated star formation rates, suggestive of interactions. Lacking are dynamical models of observed dwarf galaxy pairs to disentangle the physical processes regulating their baryon cycles. We present new optical data and the first detailed theoretical model of an observed tidal encounter between two isolated low mass galaxies, NGC 4490 & NGC 4485. This system is an isolated analog of the Magellanic Clouds and is surrounded by a ~50 kpc extended HI envelope. We use hybrid $N$-body and test-particle simulations along with a visualization interface $Identikit$ to simultaneously reproduce the observed present-day morphology and kinematics. Our results demonstrate how repeated encounters between two dwarf galaxies can "park" baryons at very large distances, without the aid of environmental effects. Our best match to the data is an 8:1 mass ratio encounter where a one-armed spiral is induced in the NGC 4490-analog, which we postulate explains the nature of diffuse starlight presented in the new optical data. We predict that the pair will fully merge in ~370 Myr, but that the extended tidal features will continue to evolve and return to the merged remnant over ~5 Gyr. This pre-processing of baryons will affect the efficiency of gas stripping if such dwarf pairs are accreted by a massive host. In contrast, in isolated environments this study demonstrates how dwarf-dwarf interactions can create a long-lived supply of gas to the merger remnant., 30 pages, 11 figures, 2 tables, minor edits to reflect published version

The demographics of dwarf galaxy populations have long been in tension with predictions from the Λ cold dark matter (ΛCDM) paradigm1–4. If primordial density fluctuations were scale-free as predicted, dwarf galaxies should themselves host dark-matter subhaloes5, the most massive of which may have undergone star formation resulting in dwarf galaxy groups. Ensembles of dwarf galaxies are observed as sate­llites of more massive galaxies6–9, and there is observational10 and theoretical11 evidence to suggest that these satellites at redshift z = 0 were captured by the massive host halo as a group. However, the evolution of dwarf galaxies is highly susceptible to environment12–14, making these satellite groups imperfect probes of ΛCDM in the low-mass regime. Here we report one of the clearest examples yet of hierarchical structure formation at low masses: using deep multi-wavelength data, we identify seven isolated, spectroscopically confirmed groups of only dwarf galaxies. Each group hosts three to five known members, has a baryonic mass of ~4.4 × 109 to 2 × 1010 solar masses (M⊙), and requires a mass-to-light ratio of

Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution ($m_{dm} = 4 \times 10^4 M_{\odot}$ ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself., 41 pages, 27 figures. Accepted to ApJ. Some terminology was changed. High-resolution images and figures can be found at https://bit.ly/2S25YzC

We present the first detailed kinematic analysis of the proper motions (PMs) of stars in the Magellanic Bridge, from both the \textit{Gaia} Data Release 2 catalog and from \textit{Hubble Space Telescope} Advanced Camera for Surveys data. For the \textit{Gaia} data, we identify and select two populations of stars in the Bridge region, young main sequence (MS) and red giant stars. The spatial locations of the stars are compared against the known H {\small I} gas structure, finding a correlation between the MS stars and the H {\small I} gas. In the \textit{Hubble Space Telescope} fields our signal comes mainly from an older MS and turn-off population, and the proper motion baselines range between $\sim 4$ and 13 years. The PMs of these different populations are found to be consistent with each other, as well as across the two telescopes. When the absolute motion of the Small Magellanic Cloud is subtracted out, the residual Bridge motions display a general pattern of pointing away from the Small Magellanic Cloud towards the Large Magellanic Cloud. We compare in detail the kinematics of the stellar samples against numerical simulations of the interactions between the Small and Large Magellanic Clouds, and find general agreement between the kinematics of the observed populations and a simulation in which the Clouds have undergone a recent direct collision., 13 pages, 10 figures, 2 tables, submitted to ApJ, accepted February 8th, 2019

The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is confidently detected, consistent with the known line-of-sight rotation curves: $V_{\rm rot} = -206\pm86$ km s$^{-1}$ (counter-clockwise) for M31, and $V_{\rm rot} = 80\pm52$ km s$^{-1}$ (clockwise) for M33. We measure the center-of-mass PM of each galaxy relative to surrounding background quasars in DR2. This yields that $({\mu}_{\alpha*},{\mu}_{\delta})$ equals $(65 \pm 18 , -57 \pm 15)$ $\mu$as yr$^{-1}$ for M31, and $(31 \pm 19 , -29 \pm 16)$ $\mu$as yr$^{-1}$ for M33. In addition to the listed random errors, each component has an additional residual systematic error of 16 $\mu$as yr$^{-1}$. These results are consistent at 0.8$\sigma$ and 1.0$\sigma$ with the (2 and 3 times higher-accuracy) measurements already available from Hubble Space Telescope (HST) optical imaging and VLBA water maser observations, respectively. This lends confidence that all these measurements are robust. The new results imply that the M31 orbit towards the Milky Way is somewhat less radial than previously inferred, $V_{\rm tan, DR2+HST} = 57^{+35}_{-31}$ km s$^{-1}$, and strengthen arguments that M33 may be on its first infall into M31. The results highlight the future potential of Gaia for PM studies beyond the Milky Way satellite system., Comment: 15 pages, 7 figures, ApJ, in press

We present the discovery of a very faint stellar system, SMASH 1, that is potentially a satellite of the Large Magellanic Cloud. Found within the Survey of the MAgellanic Stellar History (SMASH), SMASH 1 is a compact ($r_h = 9.1^{+5.9}_{-3.4}$ pc) and very low luminosity (M_V = -1.0 +/- 0.9, $L_V=10^{2.3 +/- 0.4}$ Lsun) stellar system that is revealed by its sparsely populated main sequence and a handful of red-giant-branch candidate member stars. The photometric properties of these stars are compatible with a metal-poor ([Fe/H]=-2.2) and old (13 Gyr) isochrone located at a distance modulus of ~18.8, i.e. a distance of ~57 kpc. Situated at 11.3$^\circ$ from the LMC in projection, its 3-dimensional distance from the Cloud is ~13 kpc, consistent with a connection to the LMC, whose tidal radius is at least 16 kpc. Although the nature of SMASH 1 remains uncertain, its compactness favors it being a stellar cluster and hence dark-matter free. If this is the case, its dynamical tidal radius is only, Comment: 7 pages, 4 figures, 1 table. Accepted for publication in ApJL. v2: corrected figure 4 in which rh for globular clusters were miscalculated. The updated figure reflect the published version of the paper

In this paper, we introduce the Local Volume TiNy Titans sample (LV-TNT), which is a part of a larger body of work on interacting dwarf galaxies: TNT (Stierwalt et al. 2015). This LV-TNT sample consists of 10 dwarf galaxy pairs in the Local Universe (< 30 Mpc from Milky Way), which span mass ratios of M_(*,1)/M_(*,2) < 20, projected separations < 100 kpc, and pair member masses of log(M_*/M_Sun) < 9.9. All 10 LV-TNT pairs have resolved synthesis maps of their neutral hydrogen, are located in a range of environments and captured at various interaction stages. This enables us to do a comparative study of the diffuse gas in dwarf-dwarf interactions and disentangle the gas lost due to interactions with halos of massive galaxies, from the gas lost due to mutual interaction between the dwarfs. We find that the neutral gas is extended in the interacting pairs when compared to non-paired analogs, indicating that gas is tidally pre-processed. Additionally, we find that the environment can shape the HI distributions in the form of trailing tails and that the gas is not unbound and lost to the surroundings unless the dwarf pair is residing near a massive galaxy. We conclude that a nearby, massive host galaxy is what ultimately prevents the gas from being reaccreted. Dwarf-dwarf interactions thus represent an important part of the baryon cycle of low mass galaxies, enabling the "parking" of gas at large distances to serve as a continual gas supply channel until accretion by a more massive host., 21 pages, 10 figures, accepted for publication in MNRAS

We present N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position and velocity with respect to the Galactic Center. We find that the LMC can produce strong disturbances - warping of the stellar and gaseous discs - in the Galaxy, without violating constraints from the phase-space distribution of stars in the Solar Neighbourhood. All models correctly reproduce the phases of the warp and its anti-symmetrical shape about the disc's mid-plane. If the warp is due to the LMC alone, then the largest mass model is favoured ($2.5\times10^{11}\,\rm{M_{\odot}}$). Still, some quantitative discrepancies remain, including deficits in height of $\Delta Z=0.7 \, \rm{kpc}$ at $R=22 \,\rm{kpc}$ and $\Delta Z=0.7\, \rm{kpc}$ at $R=16 \,\rm{kpc}$. This suggests that even higher infall masses for the LMC's halo are allowed by the data. A comparison with the vertical perturbations induced by a heavy Sagittarius dSph model ($10^{11}\,\rm{M_{\odot}}$) suggest that positive interference with the LMC warp is expected at $R=16 \, \rm{kpc}$. We conclude that the vertical structure of the Galactic disc beyond the Solar Neighbourhood may jointly be shaped by its most massive satellites. As such, the current structure of the MW suggests we are seeing the process of disc heating by satellite interactions in action., Comment: 13 pages, 13 figures, submitted to MNRAS, comments welcome

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields along the Orphan Stream. We determine absolute PMs of several individual stars per target field using established techniques that utilize distant background galaxies to define a stationary reference frame. Five Orphan Stream stars are identified in one of the four fields based on combined color-magnitude and PM information. The average PM is consistent with the existing model of the Orphan stream by Newberg et al. In addition to the Orphan stream stars, we detect stars that likely belong to other stellar streams. To identify which stellar streams these stars belong to, we examine the 2-d bulk motion of each group of stars on the sky by subtracting the PM contribution of the solar motion (which is a function of position on the sky and distance) from the observed PMs, and comparing the vector of net motion with the spatial extent of known stellar streams. By doing this, we identify candidate stars in the Sagittarius and Lethe streams, and a newly-found stellar stream at a distance of ~17 kpc, which we tentatively name the "Parallel stream". Together with our Sagittarius stream study (Sohn et al., 2015, ApJ, 803, 56), this work demonstrates that even in the Gaia era, HST will continue to be advantageous in measuring PMs of old stellar populations on a star-by-star basis, especially for distances beyond ~10 kpc., Comment: 10 pages, 7 figures, 3 tables, accepted for publication in ApJ

We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high-precision proper motions from the Hubble Space Telescope. Given the updated velocities, cosmological simulations of hierarchical structure formation favour a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy. We illustrate here that the observed irregular morphology and internal kinematics of the Magellanic System (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. These conclusions provide further support that the MCs are completing their first infall to our system. In particular, we demonstrate that the Magellanic Stream, a band of H i gas trailing behind the Clouds 150° across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-centre, warped stellar bar of the LMC, and its one-armed spiral can be naturally explained by a recent direct collision with its lower mass companion, the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars, the majority of which are not associated with massive spiral galaxies. We infer that dwarf–dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are not only important for the evolution of dwarf galaxies but also have direct consequences for the build-up of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type haloes.

As new work on the proper motions (PMs) of the Large Magellanic Cloud (LMC) has come out, our view of the history of the Magellanic Clouds has evolved. We now believe they are on their first infall into the Milky Way (MW), having been tidally bound at the start of infall (though not necessarily now). Combining these observations with initial PMs of the Small Magellanic Cloud (SMC) suggests a new formation mechanism of the Magellanic Stream through the stripping of material from the SMC. However, large uncertainties remain in the exact mass of the LMC. We present a measurement of the systemic proper motions of the SMC from astrometry with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST), covering a ~3 year baseline of 30 fields with background QSOs. We find these motions to be μW = −0.82 ± 0.06 mas/yr and μN = −1.23 ± 0.07 mas/yr. Combining these measurements with previous efforts in studying the Clouds will help constrain their interactions with each other and the MW, including the mass of the LMC and the MW, as well as provide new insight into the internal dynamics of the SMC.

Proper motions (PMs) are required to calculate accurate orbits of globular clusters (GCs) in the Milky Way (MW) halo. We present our HST program to create a PM database for 20 GCs at distances of RGC = 10–100 kpc. Targets are discussed along with PM measurement methods. We also describe how our PM results can be used for Gaia as an external check, and discuss the synergy between HST and Gaia as astrometric instruments in the coming years.

We examine Milky Way–Magellanic Cloud systems selected from the Millennium-II Simulation in order to place the orbits of the Magellanic Clouds in a cosmological context. Our analysis shows that satellites massive enough to be LMC analogues are typically accreted at late times. Moreover, those that are accreted at early times and survive to the present have orbital properties that are discrepant with those observed for the LMC. The high velocity of the LMC, coupled with the dearth of unbound orbits seen in the simulation, argues that the mass of the MW’s halo is unlikely to be less than 2 × 10 12 M� . This conclusion is further supported by statistics of haloes hosting satellites with masses, velocities and separations comparable to those of the LMC. We further show that: (1) LMC and SMC-mass objects are not particularly uncommon in MW-mass haloes; (2) the apparently high angular momentum of the LMC is not cosmologically unusual; and (3) it is rare for a MW halo to host a LMC–SMC binary system at z = 0, but high speed binary pairs accreted at late times are possible. Based on these results, we conclude that the LMC was accreted within the past four Gyr and is currently making its first pericentric passage about the MW.

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 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 \times 10^{12} M_\odot$) and "high" mass ($1.6 \times 10^{12} M_\odot$). 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., Comment: 18 pages, 7 figures, 4 tables, accepted for publication in A&A

We use GAIA DR2 proper motions of the RIOTS4 field OB stars in the Small Magellanic Cloud (SMC) to study the kinematics of runaway stars. The data reveal that the SMC Wing has a systemic peculiar motion relative to the SMC Bar of (v_RA, v_Dec) = (62 +/-7, -18+/-5) km/s and relative radial velocity +4.5 +/- 5.0 km/s. This unambiguously demonstrates that these two regions are kinematically distinct: the Wing is moving away from the Bar, and towards the Large Magellanic Cloud with a 3-D velocity of 64 +/- 10 km/s. This is consistent with models for a recent, direct collision between the Clouds. We present transverse velocity distributions for our field OB stars, confirming that unbound runaways comprise on the order of half our sample, possibly more. Using eclipsing binaries and double-lined spectroscopic binaries as tracers of dynamically ejected runaways, and high-mass X-ray binaries (HMXBs) as tracers of runaways accelerated by supernova kicks, we find significant contributions from both populations. The data suggest that HMXBs have lower velocity dispersion relative to dynamically ejected binaries, consistent with the former corresponding to less energetic supernova kicks that failed to unbind the components. Evidence suggests that our fast runaways are dominated by dynamical, rather than supernova, ejections., Comment: Accepted to ApJ Letters. 10 pages, 4 figures

We present a new measurement of the systemic proper motion of the Small Magellanic Cloud (SMC), based on an expanded set of 30 fields containing background quasars and spanning a $\sim$3 year baseline, using the \textit{Hubble Space Telescope} (\textit{HST}) Wide Field Camera 3. Combining this data with our previous 5 \textit{HST} fields, and an additional 8 measurements from the \textit{Gaia}-Tycho Astrometric Solution Catalog, brings us to a total of 43 SMC fields. We measure a systemic motion of $\mu_{W}$ = $-0.82$ $\pm$ 0.02 (random) $\pm$ 0.10 (systematic) mas yr$^{-1}$ and $\mu_{N}$ = $-1.21$ $\pm$ 0.01 (random) $\pm$ 0.03 (systematic) mas yr$^{-1}$. After subtraction of the systemic motion, we find little evidence for rotation, but find an ordered mean motion radially away from the SMC in the outer regions of the galaxy, indicating that the SMC is in the process of tidal disruption. We model the past interactions of the Clouds with each other based on the measured present-day relative velocity between them of $103 \pm 26$ km s$^{-1}$. We find that in 97\% of our considered cases, the Clouds experienced a direct collision $147 \pm 33$ Myr ago, with a mean impact parameter of $7.5 \pm 2.5$ kpc., Comment: 13 pages, 12 figures, 3 tables, accepted to ApJ

We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{\rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{\rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{\rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $\beta = 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (

The outer regions of galactic disks have received increased attention since ultraviolet observations with GALEX demonstrated that nearly 30% of galaxies have UV emission beyond their optical extents, indicating star formation activity. These galaxies have been termed extended UV (XUV) disks. Here, we address whether these observations contradict the gas surface density threshold for star formation inferred from Halpha radial profiles of galaxies. We run smoothed particle hydrodynamics simulations of isolated disk galaxies with fiducial star formation prescriptions and show that over-densities owing to the presence of spiral structure can induce star formation in extended gas disks. For direct comparison with observations, we use the 3-D radiative transfer code Sunrise to create simulated FUV and K_s band images. We find that galaxies classified as Type I XUV disks are a natural consequence of spiral patterns, but we are unable to reproduce Type II XUV disks. We also compare our results to studies of the Kennicutt-Schmidt relation in outer disks., Comment: Published in ApJ

Previous models of galactic disk heating in interactions invoke restrictive assumptions not necessarily valid in modern LCDM contexts: that satellites and orbits are rigid and circular, with slow decay over many orbital times from dynamical friction. This leads to a linear scaling of disk heating with satellite mass: disk heights and velocity dispersions scale ~M_sat/M_disk. In turn, observed disk thicknesses present strong constraints on merger histories: the implication for the Milky Way is that, 14 pages, 3 figures, accepted to ApJ (minor revisions to match accepted version)

The HST proper motion (PM) measurements of the Clouds have severe implications for their interaction history with the Milky Way (MW) and with each other. The Clouds are likely on their first passage about the MW and the SMC's orbit about the LMC is better described as quasi-periodic rather than circular. Binary L/SMC orbits that satisfy observational constraints on their mutual interaction history (e.g., the formation of the Magellanic Bridge during a collision between the Clouds ~300 Myr ago) can be located within 1σ of the mean PMs. However, these binary orbits are not co-located with the Magellanic Stream (MS) when projected on the plane of the sky and the line-of-sight velocity gradient along the LMC's orbit is significantly steeper than that along the MS. These combined results ultimately rule out a purely tidal origin for the MS: tides are ineffective without multiple pericentric passages and can neither decrease the velocity gradient nor explain the offset stream in a polar orbit configuration. Alternatively, ram pressure stripping of an extended gaseous disk may naturally explain the deviation. The offset also suggests that observations of the little-explored region between RA 21h and 23h are crucial for characterizing the full extent of the MS.

Recent proper motion measurements of the Large and Small Magellanic Clouds (LMC and SMC, respectively) by Kallivayalil et al (2006a,b) suggest that the 3D velocities of the Clouds are substantially higher (~100 km/s) than previously estimated and now approach the escape velocity of the Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately modeled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds using the new velocities and a LCDM-motivated MW model with virial mass Mvir = 1e12 Msun (e.g. Klypin et al 2002). We conclude that the Clouds are either currently on their first passage about the MW or, if the MW can be accurately modeled by an isothermal sphere to distances >200 kpc (i.e., Mvir > 2e12 Msun), that their orbital period and apogalacticon distance must be a factor of two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss the implications of this orbital analysis for our understanding of the star formation history, the nature of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing the LMC and SMC that extends roughly 100 degrees across the sky. Specifically, as a consequence of the new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and ram pressure stripping models.

We describe several aspects of the variability of bright pulsating red giant (PRG) stars, with emphasis on the determination of pulsation modes and pulsation masses by using observed multiperiodicity and pulsation constants (Q- values). Our studies of bright PRGs provide a model for the analysis and interpretation of the thousands of PRGs being discovered in surveys of Local Group galaxies.

We report multiperiodicity in five small‐amplitude pulsating red giants: RZ Ari, V523 Mon, BC CMi, UX Lyn, and FS Com. For each of these stars, two or three periods recur in each season of our 5000 day database of V observations. The periods and their ratios are consistent with low‐order radial pulsation modes. The amplitudes of the modes change significantly on timescales of 1–5 yr; specifically, the amplitudes of the dominant periods vary on timescales of 2000–3500 days. Most often, the amplitudes of the modes in a given star rise and fall in unison.

Interactions between pairs of isolated dwarf galaxies provide a critical window into low-mass hierarchical, gas-dominated galaxy assembly and the buildup of stellar mass in low-metallicity systems. We present the first VLT/MUSE optical IFU observations of the interacting dwarf pair dm1647+21, selected from the TiNy Titans survey. The H$\alpha$ emission is widespread and corresponds to a total unobscured star formation rate (SFR) of 0.44 M$_{\odot}$ yr$^{-1}$, 2.7 times higher than the SFR inferred from SDSS data. The implied specific SFR (sSFR) for the system is elevated by more than an order of magnitude above non-interacting dwarfs in the same mass range. This increase is dominated by the lower-mass galaxy, which has a sSFR enhancement of $>$ 50. Examining the spatially-resolved maps of classic optical line diagnostics, we find the ISM excitation can be fully explained by star formation. The velocity field of the ionized gas is not consistent with simple rotation. Dynamical simulations indicate that the irregular velocity field and the stellar structure is consistent with the identification of this system as an ongoing interaction between two dwarf galaxies. The widespread, clumpy enhancements in star formation in this system point to important differences in the effect of mergers on dwarf galaxies, compared to massive galaxies: rather than the funneling of gas to the nucleus and giving rise to a nuclear starburst, starbursts in low-mass galaxy mergers may be triggered by large-scale ISM compression, and thus be more distributed., Comment: Accepted for publication in ApJ. 11 pages, 5 figures, 1 table. Figures slightly degraded to meet arXiv size restrictions. For more information about TiNy Titans see https://lavinia.as.arizona.edu/~tinytitans/

This paper reviews our current understanding of the orbital history and mass of the Large and Small Magellanic Clouds. Specifically, I will argue that the Clouds are on their first infall about our Milky Way and that their total masses are necessarily ∼ 10 times larger than traditionally estimated. This conclusion is based on the recently revised HST proper motions of the Clouds and arguments concerning the binary status of the LMC-SMC pair and their baryon fractions.

We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations., Comment: 21 pages, 14 figures, 4 tables; Accepted for publication in the Astrophysical Journal

We present a scenario for the production of dynamically cold rings around the Milky Way via a high-eccentricity, flyby encounter. These initial conditions are more cosmologically motivated than those considered in previous works. We find that the encounters we examine generically produce a series of nearly dynamically cold ring-like features on low-eccentricity orbits that persist over timescales of ~2-4 Gyr via the tidal response of the primary galaxy to the close passage of the satellite. Moreover, they are both qualitatively and quantitatively similar to the distribution, kinematics, and stellar population of the Monoceros ring. Therefore, we find that a high eccentricity flyby by a satellite galaxy represents a cosmologically appealing scenario for forming kinematically distinct ring-like features around the Milky Way., Comment: accepted to ApJL; 4 Figures

We introduce TiNy Titans (TNT), the first systematic study of star formation and the subsequent processing of the interstellar medium in interacting dwarf galaxies. Here we present the first results from a multiwavelength observational program based on a sample of 104 dwarf galaxy pairs selected from a range of environments within the spectroscopic portion of the Sloan Digital Sky Survey and caught in various stages of interaction. The TNT dwarf pairs span mass ratios of M/M 100 Å, occur in 20% of the TNT dwarf pairs, regardless of environment, compared to only 6%–8% of the matched unpaired dwarfs. Starbursts can be triggered throughout the merger (i.e., out to large pair separations) and not just approaching coalescence. Despite their enhanced star formation and triggered starbursts, most TNT dwarf pairs have similar gas fractions relative to unpaired dwarfs of the same stellar mass. Thus, there may be significant reservoirs of diffuse, non-star-forming neutral gas surrounding the dwarf pairs, or the gas consumption timescales may be long in the starburst phase. The only TNT dwarf pairs with low gas fractions (f) and the only dwarfs, either paired or unpaired, with Hα EQW < 2 Å are found near massive galaxy hosts. We conclude that dwarf–dwarf interactions are significant drivers of galaxy evolution at the low-mass end, but ultimately environment is responsible for the quenching of star formation. This preliminary study is a precursor to an ongoing high-resolution H i and optical imaging program to constrain the spatial distribution of star formation and gas throughout the course of the dwarf–dwarf merger sequence.

We present proper motions for the Large and Small Magellanic Clouds (LMC and SMC) based on three epochs of Hubble Space Telescope data, spanning a ~7 yr baseline, and centered on fields with background QSOs. The first two epochs, the subject of past analyses, were obtained with ACS/HRC, and have been reanalyzed here. The new third epoch with WFC3/UVIS increases the time baseline and provides better control of systematics. The three-epoch data yield proper-motion random errors of only 1%-2% per field. For the LMC this is sufficient to constrain the internal proper-motion dynamics, as will be discussed in a separate paper. Here we focus on the implied center-of-mass proper motions: μ W, LMC = –1.910 ± 0.020 mas yr–1, μ N, LMC = 0.229 ± 0.047 mas yr–1, and μ W, SMC = –0.772 ± 0.063 mas yr–1, μ N, SMC = –1.117 ± 0.061 mas yr–1. We combine the results with a revised understanding of the solar motion in the Milky Way to derive Galactocentric velocities: v tot, LMC = 321 ± 24 km s–1 and v tot, SMC = 217 ± 26 km s–1. Our proper-motion uncertainties are now dominated by limitations in our understanding of the internal kinematics and geometry of the Clouds, and our velocity uncertainties are dominated by distance errors. Orbit calculations for the Clouds around the Milky Way allow a range of orbital periods, depending on the uncertain masses of the Milky Way and LMC. Periods 4 Gyr are ruled out, which poses a challenge for traditional Magellanic Stream models. First-infall orbits are preferred (as supported by other arguments as well) if one imposes the requirement that the LMC and SMC must have been a bound pair for at least several Gyr.