skip to main content


Title: Unravelling the mass spectrum of destroyed dwarf galaxies with the metallicity distribution function
ABSTRACT

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
Author(s) / Creator(s):
; ; ;
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
  1. ABSTRACT

    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.

     
    more » « 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 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
  3. 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 oftemplateCARDs from disrupted dwarfs, with different stellar massesMand quenching timest100. 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 oft100, 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.

     
    more » « less
  4. 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.

     
    more » « less
  5. null (Ed.)
    ABSTRACT Recent progress in constraining the massive accretions (>1:10) experienced by the Milky Way (MW) and the Andromeda galaxy (M31) offers an opportunity to understand the dwarf galaxy population of the Local Group. Using zoom-in dark matter-only simulations of MW-mass haloes and concentrating on subhaloes that are thought to be capable of hosting dwarf galaxies, we demonstrate that the infall of a massive progenitor is accompanied with the accretion and destruction of a large number of subhaloes. Massive accreted progenitors do not increase the total number of infalling subhaloes on to a MW-mass host, but instead focus surrounding subhaloes on to the host causing a clustering in the infall time of subhaloes. This leads to a temporary elevation in the number of subhaloes as well as changes in their cumulative radial profile within the virial radius of the host. Surviving subhaloes associated with a massive progenitor have a large diversity in their orbits. We find that the star formation quenching times of Local Group dwarf spheroidal galaxies ($10^{5} \mathrm{\, M_{\odot }} \lesssim \mathrm{\mathit{ M}}_{*} \lesssim 10^{7} \mathrm{\, M_{\odot }}$) are clustered around the times of the most massive accretions suffered by the MW and M31. Our results imply that (a) the quenching time of dwarf spheroidals is a good proxy of their infall time and b) the absence of recently quenched satellites around M31 suggests that M33 is not on its first infall and was accreted much earlier. 
    more » « less