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.
The orbits of satellite galaxies encode rich information about their histories. We investigate the orbital dynamics and histories of satellite galaxies around Milky Way (MW)-mass host galaxies using the FIRE-2 cosmological simulations, which, as previous works have shown, produce satellite mass functions and spatial distributions that broadly agree with observations. We first examine trends in orbital dynamics at z = 0, including total velocity, specific angular momentum, and specific total energy: the time of infall into the MW-mass halo primarily determines these orbital properties. We then examine orbital histories, focusing on the lookback time of first infall into a host halo and pericentre distances, times, and counts. Roughly 37 per cent of galaxies with $M_{\rm star}\lesssim 10^7\, {\rm M}_{\odot }$ were ‘pre-processed’ as a satellite in a lower-mass group, typically $\approx 2.7\, {\rm Gyr}$ before falling into the MW-mass halo. Half of all satellites at z = 0 experienced multiple pericentres about their MW-mass host. Remarkably, for most (67 per cent) of these satellites, their most recent pericentre was not their minimum pericentre: the minimum typically was ∼40 per cent smaller and occurred $\sim 6\, {\rm Gyr}$ earlier. These satellites with growing pericentres appear to have multiple origins: for about half, their specific angular momentum gradually increased over time, while for the other half, most rapidly increased near their first apocentre, suggesting that a combination of a time-dependent MW-mass halo potential and dynamical perturbations in the outer halo caused these satellites’ pericentres to grow. Our results highlight the limitations of idealized, static orbit modelling, especially for pericentre histories.
more » « less- Award ID(s):
- 2007232
- NSF-PAR ID:
- 10381343
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 518
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 1427-1447
- Size(s):
- ["p. 1427-1447"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
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.
-
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
-
null (Ed.)ABSTRACT We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the Feedback In Realistic Environments-2 simulations. We select the 14 most massive satellites by stellar mass within $d_\mathrm{host}\le 300\mathrm{\, kpc}$ of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1–2 per cent of snapshots during z = 0−0.2, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for <500 Myr. However, if we select hosts with a Large Magellanic Cloud-like satellite near first pericentre, the fraction of snapshots with MW-like planes increases dramatically to 7–16 per cent, with lifetimes of 0.7–1 Gyr, likely because of group accretion of satellites. We find that M31’s satellite distribution is much more common: M31’s satellites lie within ∼1σ of the simulation median for every plane metric we consider. We find no significant difference in average satellite planarity for isolated hosts versus hosts in LG-like pairs. Baryonic and dark matter-only simulations exhibit similar levels of planarity, even though baryonic subhaloes are less centrally concentrated within their host haloes. We conclude that planes of satellites are not a strong challenge to ΛCDM cosmology.more » « less
-
more » « less
ABSTRACT Using the eagle (Evolution and Assembly of GaLaxies and their Environments) suite of simulations, we demonstrate that both cold gas stripping and starvation of gas inflow play an important role in quenching satellite galaxies across a range of stellar and halo masses, M⋆ and M200. Quantifying the balance between gas inflows, outflows, and star formation rates, we show that even at z = 2, only $\approx 30{{\ \rm per\ cent}}$ of satellite galaxies are able to maintain equilibrium or grow their reservoir of cool gas – compared to $\approx 50{{\ \rm per\ cent}}$ of central galaxies at this redshift. We find that the number of orbits completed by a satellite on first-infall to a group environment is a very good predictor of its quenching, even more so than the time since infall. On average, we show that intermediate-mass satellites with M⋆ between will be quenched at $10^{9}\, {\rm M}_{\odot }\, {\rm and}\, 10^{10}\, {\rm M}_{\odot }$ first pericenter in massive group environments, $M_{200}\gt 10^{13.5}\, {\rm M}_{\odot }$; and will be quenched at second pericenter in less massive group environments, $M_{200}\lt 10^{13.5}\, {\rm M}_{\odot }$. On average, more massive satellites ($M_{\star }\gt 10^{10}\, {\rm M}_{\odot }$) experience longer depletion time-scales, being quenched between first and second pericenters in massive groups, while in smaller group environments, just $\approx 30{{\ \rm per\ cent}}$ will be quenched even after two orbits. Our results suggest that while starvation alone may be enough to slowly quench satellite galaxies, direct gas stripping, particularly at pericenters, is required to produce the short quenching time-scales exhibited in the simulation.
