The Milky Way (MW) and M31 both harbour massive satellite galaxies, the Large Magellanic Cloud (LMC) and M33, which may comprise up to 10 per cent of their host's total mass. Massive satellites can ...change the orbital barycentre of the host-satellite system by tens of kiloparsec and are cosmologically expected to harbour dwarf satellite galaxies of their own. Assessing the impact of these effects crucially depends on the orbital histories of the LMC and M33. Here, we revisit the dynamics of the MW-LMC system and present the first detailed analysis of the M31-M33 system utilizing high-precision proper motions and statistics from the dark-matter-only Illustris cosmological simulation. With the latest Hubble Space Telescope proper motion measurements of M31, we reliably constrain M33's interaction history with its host. In particular, like the LMC, M33 is either on its first passage (tinf < 2 Gyr ago) or if M31 is massive ( greater than or equal to 2 x 10 super( 12) M...), it is on a long-period orbit of about 6 Gyr. Cosmological analogues of the LMC and M33 identified in Illustris support this picture and provide further insight about their host masses. We conclude that, cosmologically, massive satellites such as the LMC and M33 are likely completing their first orbits about their hosts. We also find that the orbital energies of such analogues prefer an MW halo mass ~1.5 x 10 super( 12) M... and an M31 halo mass greater than or equal to 1.5 x 10 super( 12) M... Despite conventional wisdom, we conclude it is highly improbable that M33 made a close (<100 kpc) approach to M31 recently (tperi < 3 Gyr ago). Such orbits are rare (<1 per cent) within the 4s error space allowed by observations. This conclusion cannot be explained by perturbative effects through four-body encounters amongst the MW, M31, M33, and the LMC. This surprising result implies that we must search for a new explanation for M33's strongly warped gas and stellar discs. (ProQuest: ... denotes formulae/symbols omitted.)
ABSTRACT
Highly r-process enhanced metal-poor stars (MP r-II; $\rm Eu/Fe\gt 1$ and $\rm Fe/H\lesssim -1.5$) have been observed in ultra-faint dwarf (UFD) galaxy, specifically in Reticulum II ...(Ret II). The fact that only a few UFDs contain such stars implies that the r-process site may reflect very rare but individually prolific events, such as neutron star mergers (NSMs). Considering the relatively short star formation history of UFDs, it is puzzling how they could experience such a rare phenomenon. In this work, we show the results of cosmological hydrodynamic zoom-in simulations of isolated UFDs ($M_{\rm vir}\approx 10^7\!-\!10^8{\, {\rm M}_\odot }$ and $M_{\ast }\approx 10^3\!-\!10^4{\, {\rm M}_\odot }$ at z = 0) to explain the formation of MP r-II stars in UFDs. We employ a simple toy model for NSM events, adopting parameters consistent with observations, such as the NSM rate (1 per $M_{\ast }\approx 10^5{\, {\rm M}_\odot }$) and europium (Eu) mass ($M_{\rm Eu}\approx 10^{-5}{\, {\rm M}_\odot }$). We identify only one simulated galaxy ($\rm M_{vir}\approx 4.6\times 10^7{\, {\rm M}_\odot }$, $M_{\ast }\approx 3.4\times 10^3{\, {\rm M}_\odot }$ at z = 0) with abundances similar to Ret II in a simulation volume that hosts ∼30 UFD analogues, indicating that such abundances are possible but rare. By exploring a range of key parameters, we demonstrate that the most important factor in determining the formation of MP r-II stars in UFDs is how quickly subsequent stars can be formed out of r-process enriched gas. We find that it takes between 10 and 100 Myr to form the first and second bursts of MP r-II stars. Over this period, Eu-polluted gas maintains the required high abundance ratios of $\rm Eu/Fe\gt 1$.
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. ...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 Large Magellanic Cloud), 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 × 1012 (including Sagittarius dSph) and 0.96+0.29−0.28 × 1012 (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 × 103-104 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.
ABSTRACT
Carbon enhanced metal poor (CEMP)-no stars, a subset of CEMP stars ($\rm C/Fe\ge 0.7$ and $\rm Fe/H\lesssim -1$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $M_{\rm ...vir}\approx 10^8{\, \mathrm{ M}_\odot }$ and $M_{\ast }\approx 10^3-10^4{\, \mathrm{ M}_\odot }$ 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}\approx 60{\, \mathrm{ M}_\odot }$. We include enrichment from Pop III faint supernovae (SNe), with ESN = 0.6 × 1051 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 (ESN = 1.2 × 1051 erg) and Pop II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $\rm C/Fe\gtrsim 2$, corresponding to the absolute carbon abundance of $\rm A(C)\gtrsim 6.0$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $\rm C/Fe\approx 3-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.
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 (mdm = 4 × 104M ) 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 a pronounced wake, which we decompose in Transient and Collective responses, in both the DM and stellar halos. The wake leads to overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Collective response will result in redshifted radial velocities of stars in the north and blueshifts in the south, at distances >45 kpc. The Transient response 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, its orbital trajectory, and the mass of the MW, and it may inform us about the nature of the DM particle itself.
With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest-mass, ultrafaint 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 six-dimensional phase-space information, we calculate the orbital histories of 13 ultrafaint 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 < 1 Gyr ago (Reticulum 2, Phoenix 2); (iii) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and (iv) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC's gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic 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 a Large Magellanic Cloud (LMC) analog in a cosmological numerical simulation. We conclude that four 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 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, five of them are deemed unlikely on the basis of their positions and distances alone 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.
Abstract
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 MW 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 distorted MW(MW+LMC) system can be described by ∼236(2067) coefficients. We find that the LMC induces asymmetric perturbations (odd
l, m
) to the MW’s halo, which are inconsistent with oblate, prolate, or triaxial halos. Furthermore, the energy in high order even modes (
l
,
m
> 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 (>30 kpc) to shift from the disk center of mass (COM) by ∼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.
This paper explores the effect of the Large Magellanic Cloud (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. 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_{\rm LG}=2.64^{+0.42}_{-0.38}\times 10^{12}\,\mathrm{M}_{\odot }$
at a 68 per cent confidence level. The masses of the Milky Way,
$M_{\rm MW}=1.04_{-0.23}^{+0.26}\times 10^{12}\,\mathrm{M}_{\odot }$
, and Andromeda,
$M_{{\rm M}31}=1.33_{-0.33}^{+0.39}\times 10^{12}\,\mathrm{M}_{\odot }$
, are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find a (total) LMC mass
$M_{\rm LMC}=0.25_{-0.08}^{+0.09}\times 10^{12}\,\mathrm{M}_{\odot }$
, which is indicative of an extended dark matter halo and supports the scenario where this galaxy is just past its first pericentric approach. Consequently, these results suggest that the LMC may induce significant perturbations on the Galactic potential.
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 GC 10 - 100 kpc, with a median ...per-coordinate PM uncertainty of 0.06 mas yr − 1 . Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combining existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with Sgr. These clusters and NGC 6101 have tangential velocity v tan > 290 km s−1, whereas all other clusters have v 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 GC = 10 to 40 kpc. The velocity dispersion anisotropy parameter β = 0.609 − 0.229 + 0.130 . The enclosed mass M ( < 39.5 kpc ) = 0.61 − 0.12 + 0.18 × 10 12 M , and the virial mass M vir = 2.05 − 0.79 + 0.97 × 10 12 M . These are consistent with, but on the high side among, recent mass estimates in the literature.