skip to main content

Title: Planes of satellites around Milky Way/M31-mass galaxies in the FIRE simulations and comparisons with the Local Group
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 more » 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. « less
; ; ; ; ; ; ;
Award ID(s):
1752913 1910346 1715216 2108318
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
1379 to 1397
Sponsoring Org:
National Science Foundation
More Like this

    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 LGmore »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.

    « less
  2. 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 datamore »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.« less

    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 ormore »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.

    « less

    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 themore »NUDG − M200 relation.

    « less
  5. ABSTRACT Within lambda cold dark matter ($\Lambda$CDM), dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhaloes, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$$\gt 10^5$ M$_{\odot }$; Carina and Fornax) as well as several ultrafaint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$$\sim$6$\times 10^ {3}$ M$_{\odot }$) substructure population of isolated LMC-mass hosts ($M_{\text{200m}}$ = 1–3$\times 10^ {11}$ M$_{\odot }$) and place the Gaia  + DES results in a cosmological context. By comparing number counts of subhaloes in simulations with and without baryons, we find that, within 0.2 $r_{\text{200m}}$, LMC-mass hosts deplete $\sim$30 per cent of their substructure, significantly lower than the $\sim$70 per cent of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$  = 880 M$_{\odot }$), $\sim 5\!-\!10$ subhaloes form galaxies with $M_\ast$$\ge 10^{4}$ M$_{\odot }$ , in agreement with the seven observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faintmore »end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultrafaint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association.« less