skip to main content


Title: A dark matter profile to model diverse feedback-induced core sizes of ΛCDM haloes
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
Award ID(s):
1715216 1757983 1715101 1910346 1752913 1715847 1910965
NSF-PAR ID:
10184311
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
497
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2393 to 2417
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We perform a consistent comparison of the mass and mass profiles of massive (M⋆ > 1011.4 M⊙) central galaxies at z ∼ 0.4 from deep Hyper Suprime-Cam (HSC) observations and from the Illustris, TNG100, and Ponos simulations. Weak lensing measurements from HSC enable measurements at fixed halo mass and provide constraints on the strength and impact of feedback at different halo mass scales. We compare the stellar mass function (SMF) and the Stellar-to-Halo Mass Relation (SHMR) at various radii and show that the radius at which the comparison is performed is important. In general, Illustris and TNG100 display steeper values of α where $M_{\star } \propto M_{\rm vir}^{\alpha }$. These differences are more pronounced for Illustris than for TNG100 and in the inner rather than outer regions of galaxies. Differences in the inner regions may suggest that TNG100 is too efficient at quenching in situ star formation at Mvir ≃ 1013 M⊙ but not efficient enough at Mvir ≃ 1014 M⊙. The outer stellar masses are in excellent agreement with our observations at Mvir ≃ 1013 M⊙, but both Illustris and TNG100 display excess outer mass as Mvir ≃ 1014 M⊙ (by ∼0.25 and ∼0.12 dex, respectively). We argue that reducing stellar growth at early times in $M_\star \sim 10^{9-10} \, \mathrm{M}_{\odot }$ galaxies would help to prevent excess ex-situ growth at this mass scale. The Ponos simulations do not implement AGN feedback and display an excess mass of ∼0.5 dex at r < 30 kpc compared to HSC which is indicative of overcooling and excess star formation in the central regions. The comparison of the inner profiles of Ponos and HSC suggests that the physical scale over which the central AGN limits star formation is r ≲ 20 kpc. Joint comparisons between weak lensing and galaxy stellar profiles are a direct test of whether simulations build and deposit galaxy mass in the correct dark matter haloes and thereby provide powerful constraints on the physics of feedback and galaxy growth. Our galaxy and weak lensing profiles are publicly available to facilitate comparisons with other simulations. 
    more » « less
  2. ABSTRACT

    We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$), run to z = 0 with $30\, \mathrm{M}_{\odot }$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core ${\gt}200\rm \, pc$, while lower mass UFDs exhibit cusps down to ${\lesssim}100\rm \, pc$, as expected from energetic arguments. Our dwarfs with $M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

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

    Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy–halo size relation (GHSR) for low-mass ($M_\star \sim 10^{7-9}\, {\rm M}_\odot$) central galaxies over the past 12.5 billion years with the help of cosmological volume simulations (FIREbox) from the Feedback in Realistic Environments (FIRE) project. We find a nearly linear relationship between the half-stellar mass galaxy size R1/2 and the parent dark matter halo virial radius Rvir. This relation evolves only weakly since redshift z = 5: $R_{1/2}\, [{\rm kpc}] = (0.053\pm 0.002)(R_{\rm vir}/35\, {\rm kpc})^{0.934\pm 0.054}$, with a nearly constant scatter $\langle \sigma \rangle = 0.084\, [{\rm dex}]$. While this ratio is similar to what is expected from models where galaxy disc sizes are set by halo angular momentum, the low-mass galaxies in our sample are not angular momentum supported, with stellar rotational to circular velocity ratios vrot/vcirc ∼ 0.15. Introducing redshift as another parameter to the GHSR does not decrease the scatter. Furthermore, this scatter does not correlate with any of the halo properties we investigate – including spin and concentration – suggesting that baryonic processes and feedback physics are instead critical in setting the scatter in the GHSR. Given the relatively small scatter and the weak dependence of the GHSR on redshift and halo properties for these low-mass central galaxies, we propose using galaxy sizes as an independent method from stellar masses to infer halo masses.

     
    more » « less
  5. null (Ed.)
    ABSTRACT We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction cross-section per unit mass, (σ/m), and the dimensionless degree of dissipation, 0 < fdiss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $M_{\rm halo} \sim 10^{10-11}{\, \rm M_\odot }$ and $M_{\ast } \sim 10^{5-8}{\, \rm M_\odot }$. Central density profiles (parametrized as ρ ∝ rα) of simulated dwarfs become cuspy when $(\sigma /m)_{\rm eff} \gtrsim 0.1\, {\rm cm^{2}\, g^{-1}}$ (and fdiss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ −1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter ‘cooling flow’. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $(\sigma /m)_{\rm eff} \ll 0.1\, {\rm cm^{2}\, g^{-1}}$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $(\sigma /m) \gtrsim 10\, {\rm cm^{2}\, g^{-1}}$ and the fiducial fdiss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin ‘dark discs’ often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ. 
    more » « less