We use the GRUMPY galaxy formation model based on a suite of zoom-in, high-resolution, dissipationless Λ Cold Dark Matter (ΛCDM) simulations of the Milky Way (MW) sized haloes to examine total matter density within the half-mass radius of stellar distribution, ρtot(< r1/2), of satellite dwarf galaxies around the MW hosts and their mass assembly histories. We compare model results to ρtot(< r1/2) estimates for observed dwarf satellites of the Milky Way spanning their entire luminosity range. We show that observed MW dwarf satellites exhibit a trend of decreasing total matter density within a half-mass radius, ρtot(< r1/2), with increasing stellar mass. This trend is in general agreement with the trend predicted by the model. None of the observed satellites are overly dense compared to the results of our ΛCDM-based model. We also show that although the halo mass of many satellite galaxies is comparable to the halo mass of the MW progenitor at z ≳ 10, at these early epochs halos that survive as satellites to z = 0 are located many virial radii away from the MW progenitors and thus do not have a chance to merge with it. Our results show that neither the densities estimated in observed Milky Way satellites nor their mass assembly histories pose a challenge to the ΛCDM model. In fact, the broad agreement between density trends with the stellar mass of the observed and model galaxies can be considered as yet another success of the model.
Accreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the ‘stellar haloes’ of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass–metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L > Lhost/100. We find that the MW stellar halo has N ∼ 1−3 massive progenitors (L ≳ 108L⊙) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of 2. Future data sets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass scale for both the MW and its satellites.
more » « less- NSF-PAR ID:
- 10399520
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 520
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 6091-6103
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
ABSTRACT -
Abstract In the era of large-scale spectroscopic surveys in the Local Group, we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the chemical abundance ratio distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies’ assembly histories by modeling the CARDs of the stellar halos of the Latte galaxies as a linear combination of
template CARDs from disrupted dwarfs, with different stellar massesM ⋆and quenching timest 100. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of <10%. For the fraction of mass accreted as a function oft 100, we find the residuals of 20%–30% for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories that dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Although the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos. -
ABSTRACT We use GRUMPY, a simple regulator-type model for dwarf galaxy formation and evolution, to forward model the dwarf galaxy satellite population of the Milky Way (MW) using the Caterpillar zoom-in simulation suite. We show that luminosity and distance distributions of the model satellites are consistent with the distributions measured in the DES, PS1, and SDSS surveys, even without including a model for the orphan galaxies. We also show that our model for dwarf galaxy sizes can simultaneously reproduce the observed distribution of stellar half-mass radii, r1/2, of the MW satellites and the overall r1/2–M⋆ relation exhibited by observed dwarf galaxies. The model predicts that some of the observed faint stellar systems with r1/2 < 10 pc are ultra-faint dwarf galaxies. Scaling of the stellar mass M⋆ and peak halo mass Mpeak for the model satellites is not described by a power law, but has a clear flattening of M⋆–Mpeak scaling at $M_{\rm peak}\lt 10^8\, \, M_{\odot }$ imprinted by reionization. As a result, the fraction of low mass haloes ($M_{\rm peak}\lt 10^8 \, M_{\odot }$) hosting galaxies with MV < 0 is predicted to be 50 per cent at $M_{\rm peak}\sim 3.6 \times 10^7\, \, M_{\odot }$. We find that such high fraction at that halo mass helps to reproduce the number of dwarf galaxies discovered recently in the HSC-SSP survey. Using the model we forecast that there should be the total of $440^{+201}_{-147}$ (68 per cent confidence interval) MW satellites with MV < 0 and r1/2 > 10 pc within 300 kpc and make specific predictions for the HSC-SSP, DELVE-WIDE, and LSST surveys.
-
null (Ed.)ABSTRACT We investigate the internal structure of elliptical galaxies at z ∼ 0.2 from a joint lensing–dynamics analysis. We model Hubble Space Telescope images of a sample of 23 galaxy–galaxy lenses selected from the Sloan Lens ACS (SLACS) survey. Whereas the original SLACS analysis estimated the logarithmic slopes by combining the kinematics with the imaging data, we estimate the logarithmic slopes only from the imaging data. We find that the distribution of the lensing-only logarithmic slopes has a median 2.08c ± 0.03 and intrinsic scatter 0.13 ± 0.02, consistent with the original SLACS analysis. We combine the lensing constraints with the stellar kinematics and weak lensing measurements, and constrain the amount of adiabatic contraction in the dark matter (DM) haloes. We find that the DM haloes are well described by a standard Navarro–Frenk–White halo with no contraction on average for both of a constant stellar mass-to-light ratio (M/L) model and a stellar M/L gradient model. For the M/L gradient model, we find that most galaxies are consistent with no M/L gradient. Comparison of our inferred stellar masses with those obtained from the stellar population synthesis method supports a heavy initial mass function (IMF) such as the Salpeter IMF. We discuss our results in the context of previous observations and simulations, and argue that our result is consistent with a scenario in which active galactic nucleus feedback counteracts the baryonic-cooling-driven contraction in the DM haloes.more » « less
-
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