An official website of the United States government
We measure the dark matter halo mass function, with backsplash halos removed, from a wide range of cosmological-box and zoom-in simulations. These include the MultiDark Planck boxes, along with a suite of zoom-in simulations of Group, Milky Way, and LMC-mass halos. The Milky Way simulations include both CDM and non-CDM initial conditions. Using these measurements, we calibrate the parameters of flexible fitting functions for the halo mass function and the window function, along with parameterized models for various systematics, including finite box size effects, halo isolation criteria, halo detection efficiency, and contamination by artificial halos (objects forming from particle noise in the initial conditions). We show that this model shows remarkable consistency with N-body simulations over a broad range of redshifts, and ten orders of magnitude in halo mass (106M⊙ to 1016M⊙). Our model typically maintains a high precision of 12% and captures complex behaviors, including small-scale cut-offs, oscillations, and enhancements. In specific mass intervals for certain power spectra, we see larger deviations of 40-50%. Furthermore, when integrated with a simple model for environmental dependence, this fitting function provides a robust description of how environmental density influences the halo mass function. This precision model captures a wide variety of dark matter paradigms (including thermal relics, axions, and models with dark-sector interactions), is accurate for halo masses down to 107M⊙, and is a critical ingredient for model-independent dark-matter inference from forthcoming data.
more »
« less
This content will become publicly available on April 21, 2027
The DREAMS Project: Disentangling the Impact of Halo-to-halo Variance and Baryonic Feedback on Milky Way Dark Matter Density Profiles
Abstract In this work, we utilize a new suite of Milky Way–mass halos from the DREAMS Project, simulated with cold dark matter (CDM), to quantify the influence of baryon feedback and intrinsic halo-to-halo variance on dark matter density profiles. Our suite of 1024 halos varies over supernova and black hole feedback parameters from the IllustrisTNG model, as well as variations in two cosmological parameters. We find that, for the DREAMS parameter variations, Milky Way–mass dark matter density profiles in the IllustrisTNG model are largely insensitive to astrophysics and cosmology variations, with the dominant source of scatter instead arising from halo-to-halo variance. However, most of the (comparatively minor) feedback-driven variations come from the changes to supernova prescriptions. By comparing to dark-matter-only simulations, we find that the strongest supernova wind energies are so effective at preventing galaxy formation that the halos are nearly entirely collisionless dark matter. Finally, regardless of physics variation, all of the DREAMS halos are roughly consistent with a halo contracting adiabatically from the presence of baryons, unlike models that have bursty stellar feedback. This work represents a step toward assessing the uncertainty in Milky Way dark matter profiles, with direct implications for dark matter searches where systematic uncertainty in the density profile remains a major challenge.
more »
« less
- PAR ID:
- 10705068
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- The Astrophysical Journal
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 1002
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 8
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Recent observations suggest that the extended stellar halos of low-redshift massive galaxies are tightly connected to the assembly of their dark matter halos. In this paper, we use theIllustris,IllustrisTNG100, andIllustrisTNG300 simulations to compare how different stellar aperture masses trace halo mass. For massive central galaxies (M⋆≥ 1011.2M⊙), we find that a 2D outskirt stellar mass measured between 50 and 100 kpc (M⋆,[50,100]) consistently outperforms other aperture-based stellar masses. We further show thatM⋆,[50,100]correlates better with halo mass than the total mass of accreted stars (the ex situ mass), which suggests that not all accreted stars connect to halo assembly equally. While the galaxy formation recipes differ betweenIllustrisandIllustrisTNG100, the two simulations yield consistent ex situ outskirt fractions for massive galaxies (∼70% inM⋆,[50,100]). These results demonstrate the potential of using the outskirt stellar mass to deepen our understanding of the galaxy–halo connection in massive dark matter halos and to trace dark matter halos better.more » « less
-
ABSTRACT We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at z = 0 with stellar masses $$10^{10.7}\, \mathrm{M}_{\odot } \leqslant M_{\ast } \leqslant 10^{11.9}\, \mathrm{M}_{\odot }$$, the total power-law slope has a mean of 〈γ′〉 = 2.011 ± 0.007 and a scatter of $$\sigma _{\gamma ^{\prime }} = 0.171$$ over the radial range 0.4–4 times the stellar half-mass radius. Several correlations between γ′ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction, and in situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, γ′ is almost constant with redshift below z = 2. Through analysing IllustrisTNG model variations, we show that black hole kinetic winds are crucial to lowering γ′ and matching observed galaxy correlations. The effects of stellar winds on γ′ are subdominant compared to active galactic nucleus (AGN) feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter haloes are well described by steeper than NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter-only (DMO) run. Their inner density slopes anticorrelate (remain constant) with the halo mass in the FP (DMO) run, and anticorrelate with the halo concentration parameter c200 in both the types of runs. The dark matter haloes of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.more » « less
-
Abstract We analyze circular velocity profiles of seven ultradiffuse galaxies (UDGs) that are isolated and gas-rich. Assuming that the dark matter halos of these UDGs have a Navarro–Frenk–White (NFW) density profile or a Read density profile (which allows for constant-density cores), the inferred halo concentrations are systematically lower than the cosmological median, even as low as −0.6 dex (about 5 σ away) in some cases. Alternatively, similar fits can be obtained with a density profile that scales roughly as 1/ r 2 for radii larger than a few kiloparsecs. Both solutions require the radius where the halo circular velocity peaks ( R max ) to be much larger than the median expectation. Surprisingly, we find an overabundance of such large- R max halos in the IllustrisTNG dark-matter-only simulations compared to the Gaussian expectation. These halos form late and have higher spins compared to median halos of similar masses. The inner densities of the most extreme among these late-forming halos are higher than their NFW counterparts, leading to a ∼1/ r 2 density profile. However, the two well-resolved UDGs in our sample strongly prefer lower dark matter densities in the center than the simulated ones. Comparing to IllustrisTNG hydrodynamical simulations, we also find a tension in getting both low enough circular velocities and high enough halo mass to accommodate the measurements. Our results indicate that the gas-rich UDGs present a significant challenge for galaxy formation models.more » « less
-
ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $$0.5{{\ \rm per\ cent}}$$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($$M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$$) to the largest spirals ($$M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $$M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $$\gtrsim 100\ \rm pc$$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $$M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation.more » « less
