Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modelling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass haloes. When a typical surviving satellite fell in ($3.4\!-\!9.7\, \rm {Gyr}$ ago), the host halo mass and radius were typically 26–86 per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $r\lesssim 50\, \rm {kpc}$, opposite to dark matter only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at z = 0 and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $1.6\!-\!4.7\, \rm {Gyr}$ ago. Most orbital properties are minimally biased, ≲10 per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are ≈20 per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are ≈40–80 per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at z = 0.
- NSF-PAR ID:
- 10258059
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 504
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 1379 to 1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
ABSTRACT In the currently favoured cosmological paradigm galaxies form hierarchically through the accretion of satellites. Since a satellite is less massive than the host, its stars occupy a smaller volume in action space. Actions are conserved when the potential of the host halo changes adiabatically, so stars from an accreted satellite would remain clustered in action space as the host evolves. In this paper, we identify recently disrupted accreted satellites in three Milky Way-like disc galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellites through simulation snapshots. We try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the z = 0 snapshot. Even with completely error-free mock data we find that only 35 per cent (14/39) satellites are well recovered while the rest (25/39) are poorly recovered (i.e. either contaminated or split up). Most (10/14 ∼70 per cent) of the well-recovered satellites have infall times <7.1 Gyr ago and total mass >4 × 108M⊙ (stellar mass more than 1.2 × 106 M⊙, although our upper mass limit is likely to be resolution dependent). Since cosmological simulations predict that stellar haloes include a population of in situ stars, we test our ability to recover satellites when the data include 10–50 per cent in situ contamination. We find that most previously well-recovered satellites stay well recovered even with 50 per cent contamination. With the wealth of 6D phase space data becoming available we expect that cluster analysis in action space will be useful in identifying the majority of recently accreted and moderately massive satellites in the Milky Way.more » « less
-
more » « less
ABSTRACT The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($M_{*}\lesssim 10^7\, \rm {M}_{\odot }$) are mostly quiescent and higher mass satellites ($M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$) are mostly star forming, with intermediate-mass satellites ($M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced ≥1–2 pericentre passages (≈2.5–5 Gyr) within a MW-mass halo. Most galaxies with $M_{*}\gtrsim 10^{6.5}\, \rm {M}_{\odot }$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the Satellites Around Galactic Analogs (SAGA) survey’s quiescent fraction at $M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$, they do not generally reproduce SAGA’s turnover at lower masses.
-
more » « less
ABSTRACT While many tensions between Local Group (LG) satellite galaxies and Λ cold dark matter cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with $M_*\gt 10^5$ M⊙ around eight isolated Milky Way (MW) mass host galaxies and four hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with $M_*\gtrsim 10^5$ M⊙. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances ≲100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic discs of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhaloes in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at $M_*\gtrsim 10^5$ M⊙: we predict 2–10 such satellites to be discovered around the MW and possibly 6–9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogues in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.
-
more » « less
ABSTRACT We extend the Ultra-Diffuse Galaxy (UDG) abundance relation, NUDG − M200, to lower halo mass hosts $(M_{200}\sim 10^{11.6-12.2}\, \mathrm{M}_{\odot })$. We select UDG satellites from published catalogues of dwarf satellite galaxies around Milky Way analogues, namely the Exploration of Local Volume Satellites (ELVES) survey, the Satellite Around Galactic Analogs (SAGA) survey, and a survey of Milky Way-like systems conducted using the Hyper-Suprime Cam. Of the 516 satellites around a total of 75 Milky Way-like hosts, we find that 41 satellites around 33 hosts satisfy the UDG criteria. The distributions of host halo masses peak around $M_{200}\sim 10^{12}\, \mathrm{M}_{\odot }$, independent of whether the host has a UDG satellite or not. We use literature UDG abundances and those derived here to trace the NUDG − M200 relation over three orders of magnitude down to $M_{200}=10^{11.6}\, \mathrm{M}_{\odot }$ and find the best-fitting linear relation of $N_{\mathrm{ UDG}} = (38\pm 5) (\frac{M_{200}}{10^{14}})^{0.89\,\,\pm ~ 0.04}$. This sub-linear slope is consistent with earlier studies of UDG abundances as well as abundance relations for brighter dwarf galaxies, excluding UDG-formation mechanisms that require high-density environments. However, we highlight the need for further homogeneous characterization of UDGs across a wide range of environments to properly understand the NUDG − M200 relation.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stab955
