ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of $${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $$\gtrsim {\dot{M}}_{\rm crit}$$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare $${\dot{M}}_{\rm crit}$$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that $${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$.
more »
« less
The galaxy–halo size relation of low-mass galaxies in FIRE
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
- PAR ID:
- 10361414
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 510
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 3967-3985
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 apply the barred Schwarzschild method developed by Tahmasebzadeh et al. (2022) to a barred S0 galaxy, NGC 4371, observed by IFU instruments from the TIMER and ATLAS3D projects. We construct the gravitational potential by combining a fixed black hole mass, a spherical dark matter halo, and stellar mass distribution deprojected from 3.6 μm S$^4$G image considering an axisymmetric disc and a triaxial bar. We independently modelled kinematic data from TIMER and ATLAS3D. Both models fit the data remarkably well. We find a consistent bar pattern speed from the two sets of models with $$\Omega _{\rm p} = 23.6 \pm 2.8 \, \mathrm{km \, s^{-1} \, kpc^{-1} }$$ and $$\Omega _{\rm p} = 22.4 \pm 3.5 \, \mathrm{km \, s^{-1} \, kpc^{-1} }$$, respectively. The dimensionless bar rotation parameter is determined to be $$\mathcal {R} \equiv R_{\rm cor}/R_{\rm bar}=1.88 \pm 0.37$$, indicating a likely slow bar in NGC 4371. Additionally, our model predicts a high amount of dark matter within the bar region ($$M_{\rm DM}/ M_{\rm total}$$\sim 0.51 \pm 0.06$$), which, aligned with the predictions of cosmological simulations, indicates that fast bars are generally found in baryon-dominated discs. Based on the best-fitting model, we further decompose the galaxy into multiple 3D orbital structures, including a BP/X bar, a classical bulge, a nuclear disc, and a main disc. The BP/X bar is not perfectly included in the input 3D density model, but BP/X-supporting orbits are picked through the fitting to the kinematic data. This is the first time a real barred galaxy has been modelled utilizing the Schwarzschild method including a 3D bar.more » « less
-
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
-
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