skip to main content

Title: Black Hole Growth, Baryon Lifting, Star Formation, and IllustrisTNG

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
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 28
["Article No. 28"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Supermassive black holes (SMBHs) that reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star formation in massive galaxies. Here, we study the $10^{9-12.5}\, \mathrm{M_\odot }$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at z = 0, and show how the three components – SMBH, galaxy, and circumgalactic medium (CGM) – are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion black hole (BH) feedback mode, realized in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star formation quenching at stellar masses $\gtrsim 10^{10.5}\, \mathrm{M_\odot }$ but also to a change in thermodynamic properties of the (non-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\!-\!100\, \mathrm{Myr}$ to $1\!-\!10\, \mathrm{Gyr}$, effectively ceasing ongoing star formation and inhibiting radiative cooling and future gas accretion. In practice, the same active galactic nucleus (AGN) feedback channel is simultaneously ‘ejective’ and ‘preventative’ and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times. 
    more » « less
  2. Abstract

    We present measurements of black hole masses and Eddington ratios (λEdd) for a sample of 38 bright (M1450< −24.4 mag) quasars at 5.8 ≲z≲ 7.5, derived from Very Large Telescope/X–shooter near–IR spectroscopy of their broad Civand Mgiiemission lines. The black hole masses (on average,MBH∼ 4.6 × 109M) and accretion rates (0.1 ≲λEdd≲ 1.0) are broadly consistent with that of similarly luminous 0.3 ≲z≲ 2.3 quasars, but there is evidence for a mild increase in the Eddington ratio abovez≳ 6. Combined with deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 158μm line from the host galaxies and VLT/MUSE investigations of the extended Lyαhalos, this study provides fundamental clues to models of the formation and growth of the first massive galaxies and black holes. Compared to local scaling relations,z≳ 5.7 black holes appear to be over-massive relative to their hosts, with accretion properties that do not change with host galaxy morphologies. Assuming that the kinematics of theT∼ 104K gas, traced by the extended Lyαhalos, are dominated by the gravitational potential of the dark matter halo, we observe a similar relation between black hole mass and circular velocity as reported forz∼ 0 galaxies. These results paint a picture where the first supermassive black holes reside in massive halos atz≳ 6 and lead the first stages of galaxy formation by rapidly growing in mass with a duty cycle of order unity. The duty cycle needs to drastically drop toward lower redshifts, while the host galaxies continue forming stars at a rate of hundreds of solar masses per year, sustained by the large reservoirs of cool gas surrounding them.

    more » « less
  3. Abstract

    We study how supersonic streaming velocities of baryons relative to dark matter—a large-scale effect imprinted at recombination and coherent over ∼3 Mpc scales—affect the formation of dwarf galaxies atz≳ 5. We perform cosmological hydrodynamic simulations, including and excluding streaming velocities, in regions centered on halos withMvir(z= 0) ≈ 1010M; the simulations are part of the Feedback In Realistic Environments (FIRE) project and run with FIRE-3 physics. Our simulations comprise many thousands of systems with halo masses betweenMvir= 2 × 105Mand 2 × 109Min the redshift rangez= 20–5. A few hundred of these galaxies form stars and have stellar masses ranging from 100 to 107M. While star formation is globally delayed by approximately 50 Myr in the streaming relative to nonstreaming simulations and the number of luminous galaxies is correspondingly suppressed at high redshift in the streaming runs, these effects decay with time. Byz= 5, the properties of the simulated galaxies are nearly identical in the streaming versus nonstreaming runs, indicating that any effects of streaming velocities on the properties of galaxies at the mass scale of classical dwarfs and larger do not persist toz= 0.

    more » « less

    Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are undermassive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multichannel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fuelling: constant thresholds in these properties are typically crossed within ∼0.1 Hubble time of accelerated BH fuelling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold $\Sigma^\star _{1\,\rm kpc}\approx 10^{9.5} \, {\rm M_{\odot }}\,{\rm kpc}^{-2}$, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived thin gas discs, as well as with virialization of the inner circumgalactic medium. The halo mass Mhalo ∼ 1012 M⊙ and stellar mass M* ∼ 1010.5 M⊙ at which BH growth accelerates correspond to ∼L⋆ galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above L⋆.

    more » « less
  5. ABSTRACT We compare the star-forming main sequence (SFMS) of galaxies – both integrated and resolved on 1 kpc scales – between the high-resolution TNG50 simulation of IllustrisTNG and observations from the 3D-HST slitless spectroscopic survey at z ∼ 1. Contrasting integrated star formation rates (SFRs), we find that the slope and normalization of the star-forming main sequence in TNG50 are quantitatively consistent with values derived by fitting observations from 3D-HST with the Prospector Bayesian inference framework. The previous offsets of 0.2–1 dex between observed and simulated main-sequence normalizations are resolved when using the updated masses and SFRs from Prospector. The scatter is generically smaller in TNG50 than in 3D-HST for more massive galaxies with M*> 1010 M⊙, by ∼10–40 per cent, after accounting for observational uncertainties. When comparing resolved star formation, we also find good agreement between TNG50 and 3D-HST: average specific star formation rate (sSFR) radial profiles of galaxies at all masses and radii below, on, and above the SFMS are similar in both normalization and shape. Most noteworthy, massive galaxies with M*> 1010.5 M⊙, which have fallen below the SFMS due to ongoing quenching, exhibit a clear central SFR suppression, in both TNG50 and 3D-HST. In contrast, the original Illustris simulation and a variant TNG run without black hole kinetic wind feedback, do not reproduce the central SFR profile suppression seen in data. In TNG, inside-out quenching is due to the supermassive black hole (SMBH) feedback model operating at low accretion rates. 
    more » « less