skip to main content


Title: Why Do Semianalytic Models Predict Higher Scatter in the Stellar Mass–Halo Mass Relation Than Cosmological Hydrodynamic Simulations?
Abstract

Semianalytic models (SAMs) systematically predict higher-stellar mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAMDark Sage. We design two black hole formation models that approximate results from theIllustrisTNG 300-1hydrodynamical simulation. In the first model, we assign a fixed black hole mass of 106Mto every halo that reaches 1010.5M. In the second model, we disregard any black hole growth as implemented in the standardDark Sagemodel. Instead, we force all black hole masses to follow the medianz= 0 black hole mass–halo mass relation inIllustrisTNG 300-1with an imposed fixed scatter. We find that each model on its own does not significantly reduce the scatter in stellar mass. To explore the effects of active galactic nucleus (AGN) feedback in addition to black hole seeding, we replace the native Dark Sage AGN feedback model with a simple model where we turn off cooling for galaxies with black hole masses above 108M. With the additional modification in AGN feedback, we find that the supermassive black hole seeding and fixed conditional distribution models create a significant reduction in the scatter in stellar mass at halo masses between 1011–14M. These results suggest that AGN feedback in SAMs acts in a qualitatively different way than feedback implemented in cosmological simulations. Either or both may require substantial modification to match the empirically inferred scatter in the stellar mass–halo mass relation.

 
more » « less
PAR ID:
10556048
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
976
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 148
Size(s):
Article No. 148
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Quenching of star formation in the central galaxies of cosmological halos is thought to result from energy released as gas accretes onto a supermassive black hole. The same energy source also appears to lower the central density and raise the cooling time of baryonic atmospheres in massive halos, thereby limiting both star formation and black hole growth, by lifting the baryons in those halos to greater altitudes. One predicted signature of that feedback mechanism is a nearly linear relationship between the central black hole’s mass (MBH) and the original binding energy of the halo’s baryons. We present the increasingly strong observational evidence supporting a such a relationship, showing that it extends up to halos of massMhalo∼ 1014M. We then compare current observational constraints on theMBHMhalorelation with numerical simulations, finding that black hole masses in IllustrisTNG appear to exceed those constraints atMhalo< 1013Mand that black hole masses in EAGLE fall short of observations atMhalo∼ 1014M. A closer look at IllustrisTNG shows that quenching of star formation and suppression of black hole growth do indeed coincide with black hole energy input that lifts the halo’s baryons. However, IllustrisTNG does not reproduce the observedMBHMhalorelation because its black holes gain mass primarily through accretion that does not contribute to baryon lifting. We suggest adjustments to some of the parameters in the IllustrisTNG feedback algorithm that may allow the resulting black hole masses to reflect the inherent links between black hole growth, baryon lifting, and star formation among the massive galaxies in those simulations.

     
    more » « less
  2. 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
  3. Abstract

    The properties of warm-hot gas around ∼L*galaxies can be studied with absorption lines from highly ionized metals. We predict Neviiicolumn densities from cosmological zoom-in simulations of halos with masses in ∼1012and ∼1013Mfrom the Feedback in Realistic Environments (FIRE) project. Neviiitraces the volume-filling, virial-temperature gas in ∼1012Mhalos. In ∼1013Mhalos the Neviiigas is clumpier, and biased toward the cooler part of the warm-hot phase. We compare the simulations to observations from the COS Absorption Survey of Baryon Harbors (or CASBaH) and COS Ultraviolet Baryon Survey (or CUBS). We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass–halo-mass relation, especially atM≳ 1010.5M. Median Neviiicolumns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Neviiiprofiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) ∼0.1, indicating that observations are consistent with plausible circumgalactic medium temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or active galactic nuclei feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Neviiicolumns. The circumgalactic Neviiiobservations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties.

     
    more » « less
  4. Abstract

    The fraction of dwarf galaxies hosting central, intermediate-mass black holes (IMBHs) at low redshifts is an important observational probe of black hole seeding at high redshift. Detections of nuclear accretion signatures in dwarf galaxies provides strong evidence for the presence of these IMBHs. We develop a Bayesian model to infer the black hole occupation fraction assuming a broken power-law Eddington ratio distribution function. Our approach accounts for nondetections, incompleteness, and contamination from star-forming-related emission. We apply this model to galaxies with X-ray data from the Chandra Source Catalog at distances <50 Mpc, radio data from the Very Large Array Sky Survey at <50 Mpc, and optical variability data from the Palomar Transient Factory atz< 0.055. We find a black hole occupation fraction of at least 90% at stellar masses ofM= 108Mand at least 39% atM= 107M(95% confidence intervals). We show the resulting black hole mass function. These constraints on the IMBH population have implications for the Laser Interferometer Space Antenna mission and for cosmological models of black hole seeding and growth. We also constrain the extremely low-luminosity end (Lbol≲ 1040erg s−1) of the active galactic nucleus (AGN) luminosity functions atz= 0. Our AGN luminosity functions are broadly consistent with an extrapolation of the shallow slope of the AGN luminosity functions from previous work.

     
    more » « less
  5. Abstract

    We present radial profiles of luminosity-weighted age (ageL) and ΔΣSFRfor various populations of high- and low-mass central and satellite galaxies in the TNG100 cosmological simulation. Using these profiles, we investigate the impact of intrinsic and environmental factors on the radial distribution of star formation. For both central galaxies and satellites, we investigate the effects of black hole mass, cumulative active galactic nucleus (AGN) feedback energy, morphology, halo mass, and local galaxy overdensity on the profiles. In addition, we investigate the dependence of radial profiles of the satellite galaxies as a function of the redshifts at which they joined their hosts, as well as the net change in star-forming gas mass since the satellites joined their host. We find that high-mass (M*> 1010.5M) central and satellite galaxies show evidence of inside-out quenching driven by AGN feedback. Effects from environmental processes only become apparent in averaged profiles at extreme halo masses and local overdensities. We find that the dominant quenching process for low-mass galaxies (M*< 1010M) is environmental, generally occurring at low halo mass and high local galaxy overdensity for low-mass central galaxies and at high host halo masses for low-mass satellite galaxies. Overall, we find that environmental processes generally drive quenching from the outside-in.

     
    more » « less