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.
In this study, we modify the semi-analytic model galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass–metallicity relation for dwarf galaxies requires us to include H2 cooling, an updated ultraviolet background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size–mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy–halo connection. We predict a total of $300 ^{+75} _{-99}$ satellites with an absolute V-band magnitude (MV) less than 0 within 300 kpc from our Milky Way analogues. The fraction of subhaloes that host a galaxy at least this bright drops to 50 per cent by a halo peak mass of ∼8.9 × 107 M⊙, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population.
more » « less- NSF-PAR ID:
- 10496657
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 529
- Issue:
- 4
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3387-3407
- Size(s):
- ["p. 3387-3407"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT While dwarf galaxies observed in the field are overwhelmingly star forming, dwarf galaxies in environments as dense or denser than the Milky Way are overwhelmingly quenched. In this paper, we explore quenching in the lower density environment of the Small-Magellanic-Cloud-mass galaxy NGC 3109 (M$_* \sim 10^8 \, \text{M}_\odot$), which hosts two known dwarf satellite galaxies (Antlia and Antlia B), both of which are ${\rm H}\, \rm{\small I}$ deficient compared to similar galaxies in the field and have recently stopped forming stars. Using a new semi-analytic model in concert with the measured star formation histories and gas masses of the two dwarf satellite galaxies, we show that they could not have been quenched solely by direct ram pressure stripping of their interstellar media, as is common in denser environments. Instead, we find that separation of the satellites from pristine gas inflows, coupled with stellar-feedback-driven outflows from the satellites (jointly referred to as the starvation quenching model), can quench the satellites on time-scales consistent with their likely infall times into NGC 3109’s halo. It is currently believed that starvation is caused by ‘weak’ ram pressure that prevents low-density, weakly bound gas from being accreted on to the dwarf satellite, but cannot directly remove the denser interstellar medium. This suggests that star-formation-driven outflows serve two purposes in quenching satellites in low-mass environments: outflows from the host form a low-density circumgalactic medium that cannot directly strip the interstellar media from its satellites, but is sufficient to remove loosely bound gaseous outflows from the dwarf satellites driven by their own star formation.
-
more » « less
ABSTRACT The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low-mass haloes are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass (M*)–halo mass (M200) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the hydrogen-cooling limit. In the second model, M* scales as a steep power law of M200 with no explicit cut-off, as suggested by recent semi-analytical work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultrafaint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the M*–M200 relation at the faint end and to elucidate the mechanisms that determine which low-mass haloes ‘light up’ or remain dark in the LCDM scenario.
-
more » « less
ABSTRACT We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding time-scales of τin = 1.01 ± 0.13 Gyr for GSE and $\tau _\text{in} = 3.08^{+3.19}_{-1.16}$ Gyr for Wukong/LMS-1. GSE formed stars for $\tau _\text{tot} = 5.40^{+0.32}_{-0.31}$ Gyr, sustaining star formation for ∼1.5–2 Gyr after its first infall into the Milky Way ∼10 Gyr ago. Our fit suggests that star formation lasted for $\tau _\text{tot} = 3.36^{+0.55}_{-0.47}$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass M⋆ predicted by simulations and semi-analytic models of galaxy formation. Our inferred values of the outflow mass-loading factor reasonably match $\eta \propto M_\star ^{-1/3}$ as predicted by galactic wind models. Our fitting method is based only on Poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes (20 ≤ N ≤ 2000) and measurement uncertainties (0.01 ≤ σ[α/Fe], σ[Fe/H] ≤ 0.5; $0.02 \le \sigma _{\log _{10}(\text{age})} \le 1$). Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space.
-
more » « less
Abstract Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy’s star-forming efficiency. Mass, energy, and metal outflows (
η M ,η E , andη Z , here normalized by the star formation rate, the SNe energy, and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas, which enters the interstellar medium (ISM), although whether such high mass loadings agree with observations is still unclear. To better understand how the relative importance of ejective (i.e., high mass loading) versus preventative (i.e., high energy loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from 1010to 1012M ⊙, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low-mass-loaded (η M ∼ 0.1–10) but high-energy-loaded (η E ∼ 0.1–1) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass loading of outflows but are quite sensitive to our choice of the energy loading, preferringη E ∼ 1 for the lowest-mass halos and ∼0.1 for Milky Way–like halos.
