An official website of the United States government
ABSTRACT The James Webb Space Telescope will have the power to characterize high-redshift quasars at z ≥ 6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e. with bolometric luminosity $$L_{\rm bol}\geqslant 10^{46}\, \rm erg/s$$) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with $$L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$$ opens a path to understand the build-up of more normal BHs at z ≥ 6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z ≥ 4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the MBH − M⋆ scaling relation at z = 0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z ≥ 4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with $$L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$$ observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z ∼ 6, with the same diversity of subgrid physics, in order to gain statistics on the most extreme objects at high redshift.
more »
« less
Supermassive black holes in cosmological simulations I: M BH − M ⋆ relation and black hole mass function
ABSTRACT The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the MBH − M⋆ relation and the BH mass function. All simulations predict tight MBH − M⋆ relations, and struggle to produce BHs of $$M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$$ in galaxies of $$M_{\star }\sim 10^{10.5}\!-\!10^{11.5}\, \rm M_{\odot }$$. While the time evolution of the mean MBH − M⋆ relation is mild ($$\rm \Delta M_{\rm BH}\leqslant 1\, dex$$ for 0 $$\leqslant z \leqslant$$ 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for $$M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$$. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with $$M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$$ in most of the simulations. The BH mass function is dominated by efficiently accreting BHs ($$\log _{10}\, f_{\rm Edd}\geqslant -2$$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.
more »
« less
- PAR ID:
- 10237288
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 503
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 1940 to 1975
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT While the first “seeds” of supermassive black holes (BH) can range from $$\sim 10^2-10^6 \rm ~{\rm M}_{\odot }$$, the lowest mass seeds ($$\lesssim 10^3~\rm {\rm M}_{\odot }$$) are inaccessible to most cosmological simulations due to resolution limitations. We present our new BRAHMA simulations that use a novel flexible seeding approach to predict the $$z\ge 7$$ BH populations for low-mass seeds. We ran two types of boxes that model $$\sim 10^3~\rm {\rm M}_{\odot }$$ seeds using two distinct but mutually consistent seeding prescriptions at different simulation resolutions. First, we have the highest resolution $$[9~\mathrm{Mpc}]^3$$ (BRAHMA-9-D3) boxes that directly resolve $$\sim 10^3~\rm {\rm M}_{\odot }$$ seeds and place them within haloes with dense, metal-poor gas. Second, we have lower resolution, larger volume $$[18~\mathrm{Mpc}]^3$$ (BRAHMA-18-E4), and $$\sim [36~\mathrm{Mpc}]^3$$ (BRAHMA-36-E5) boxes that seed their smallest resolvable $$\sim 10^4~\&~10^5~\mathrm{{\rm M}_{\odot }}$$ BH descendants using new stochastic seeding prescriptions calibrated using BRAHMA-9-D3. The three boxes together probe key BH observables between $$\sim 10^3\,\mathrm{ and}\,10^7~\rm {\rm M}_{\odot }$$. The active galactic nuclei (AGN) luminosity function variations are small (factors of $$\sim 2-3$$) at the anticipated detection limits of potential future X-ray facilities ($$\sim 10^{43}~ \mathrm{ergs~s^{-1}}$$ at $$z\sim 7$$). Our simulations predict BHs $$\sim 10-100$$ times heavier than the local $$M_*$$ versus $$M_{\mathrm{ bh}}$$ relations, consistent with several JWST-detected AGN. For different seed models, our simulations merge binaries at $$\sim 1-15~\mathrm{kpc}$$, with rates of $$\sim 200-2000$$ yr−1 for $$\gtrsim 10^3~\rm {\rm M}_{\odot }$$ BHs, $$\sim 6-60$$ yr−1 for $$\gtrsim 10^4~\rm {\rm M}_{\odot }$$ BHs, and up to $$\sim 10$$ yr−1 amongst $$\gtrsim 10^5~\rm {\rm M}_{\odot }$$ BHs. These results suggest that Laser Interferometer Space Antenna mission has promising prospects for constraining seed models.more » « less
-
ABSTRACT The origin of the ‘seeds’ of supermassive black holes (BHs) continues to be a puzzle, as it is currently unclear if the imprints of early seed formation could survive to today. We examine the signatures of seeding in the local Universe using five $$[18~\mathrm{Mpc}]^3$$BRAHMA simulation boxes run to $z=0$. They initialize $$1.5\times 10^5~\rm {M}_{\odot }$$ BHs using different seeding models. The first four boxes initialize BHs as heavy seeds using criteria that depend on dense and metal-poor gas, Lyman–Werner radiation, gas spin, and environmental richness. The fifth box initializes BHs as descendants of lower mass seeds ($$\sim 10^3~\rm {M}_{\odot }$$) using a new stochastic seed model built in our previous work. In our simulations, we find that the abundances and properties of $$\sim 10^5-10^6~\rm {M}_{\odot }$$ local BHs hosted in $$M_*\lesssim 10^{9}~\rm {M}_{\odot }$$ dwarf galaxies, are sensitive to the assumed seeding criteria. This is for two reasons: (1) there is a substantial population of local $$\sim 10^5~\rm {M}_{\odot }$$ BHs that are ungrown relics of early seeds from $$z\sim 5-10$$; (2) BH growth up to $$\sim 10^6~\rm {M}_{\odot }$$ is dominated by mergers in our simulations all the way down to $$z\sim 0$$. As the contribution from gas accretion increases, the signatures of seeding start to weaken in more massive $$\gtrsim 10^6~\rm {M}_{\odot }$$ BHs, and they are erased for $$\gtrsim 10^7~\rm {M}_{\odot }$$ BHs. The different seed models explored here predict abundances of local $$\sim 10^6~\rm {M}_{\odot }$$ BHs ranging from $$\sim 0.01-0.05~\mathrm{Mpc}^{-3}$$ with occupation fractions of $$\sim 20-100~{{\ \rm per\ cent}}$$ for $$M_*\sim 10^{9}~\rm {M}_{\odot }$$ galaxies. These results highlight the potential for placing constraints on seeding models using local $$\sim 10^5-10^6~\rm {M}_{\odot }$$ BHs hosted in dwarf galaxies. Since merger dynamics and accretion physics impact the persistence of seeding signatures, and both high and low mass seed models can produce similar local BH populations, disentangling their roles will require combining high and low redshift constraints.more » « less
-
ABSTRACT We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (∼1000 galaxies with $$M_{\rm star}\gt 10^8\, M_\odot$$ at z = 0) at a resolution ($$\Delta {}x\sim {}20\, {\rm pc}$$ , $$m_{\rm b}\sim {}6\times {}10^4\, M_\odot$$ ) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L = 22.1 Mpc) thanks to its exceptionally high dynamic range (≳106) and the inclusion of multichannel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that star formation rates, gas masses, and metallicities of simulated galaxies with $$M_{\rm star}\lt 10^{10.5-11}\, M_\odot$$ broadly agree with observations. These galaxy scaling relations extend to low masses ($$M_{\rm star}\sim {}10^7\, M_\odot$$ ) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z ∼ 0–5. At low z , FIREbox predicts a peak in the stellar-mass–halo-mass relation but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a ΛCDM Universe while also highlighting modelling challenges to be addressed in next-generation galaxy simulations.more » « less
-
ABSTRACT The early growth of black holes (BHs) in atomic-cooling haloes is likely influenced by feedback on the surrounding gas. While the effects of radiative feedback are well-documented, mechanical feedback, particularly from active galactic nucleus (AGN) jets, has been comparatively less explored. Building on our previous work that examined the growth of a 100 $${\rm M_\odot }$$ BH in a constant density environment regulated by AGN jets, we expand the initial BH mass range from 1 to $$10^4\, {\rm M_\odot }$$ and adopt a more realistic density profile for atomic-cooling haloes. We reaffirm the validity of our analytic models for jet cocoon propagation and feedback regulation. We identify several critical radii – namely, the terminal radius of jet cocoon propagation, the isotropization radius of the jet cocoon, and the core radius of the atomic-cooling halo – that are crucial in determining BH growth given specific gas properties and jet feedback parameters. In a significant portion of the parameter space, our findings show that jet feedback substantially disrupts the halo’s core during the initial feedback episode, preventing BH growth beyond $$10^4 \, {\rm M_\odot }$$. Conversely, conditions characterized by low jet velocities and high gas densities enable sustained BH growth over extended periods. We provide a prediction for the BH mass growth as a function of time and feedback parameters. We found that, to form a supermassive BH ($$\gt 10^6 \, {\rm M_\odot }$$) within 1 Gyr entirely by accreting gas from an atomic-cooling halo, the jet energy feedback efficiency must be $$\lesssim 10^{-4} \dot{M}_{\rm BH} c^2$$ even if the seed BH mass is $$10^4 \, {\rm M_\odot }$$.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stab496
