An official website of the United States government
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
Building a Semi-analytic Black Hole Seeding Model using IllustrisTNG Host Galaxies
ABSTRACT A major open question in astrophysics is the mechanisms by which massive black holes (BHs) form in the early Universe, which pose constraints on seeding models. We study BH formation and evolution in a flexible model combining the cosmological IllustrisTNG (TNG) simulations with semi-analytic modelling in post-processing. We identify our TNG model hosts based on various criteria including a minimum gas mass of $10^7$$–$$10^9$${\rm M}_{\odot }$$, total host mass of $$10^{8.5}$$–$$10^{10.5}$${\rm M}_{\odot }$$, and a maximum gas metallicity of 0.01–0.1 $$\mathrm{Z}_{\odot }$$. Each potential host is assigned a BH seed with a probability of 0.01–1. The populations follow the TNG galaxy merger tree. This approach improves upon the predictive power of the simple TNG BH seeding prescription, narrowing down plausible seeding parameter spaces, and it is readily adaptable to other cosmological simulations. Several model realizations predict $$z\lesssim 4$$ BH mass densities that are consistent with empirical data as well as the TNG BHs. However, high-redshift BH number densities can differ by factors of $$\sim$$ 10 to $$\gtrsim$$ 100 between seeding parameters. In most model realizations, $$\lesssim 10^5$${\rm M}_{\odot }$$ BHs substantially outnumber heavier BHs at high redshifts. Mergers between such BHs are prime targets for gravitational-wave detection with Laser Interferometer Space Antenna. The $z=0$ BH mass densities in most realizations of the model agree well with observations, but our strictest seeding criteria fail at high redshift. Our findings strongly motivate the need for better empirical constraints on high-z BHs, and they underscore the significance of recent active galactic nucleus discoveries with JWST.
more »
« less
- Award ID(s):
- 1909933
- PAR ID:
- 10562676
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 536
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2783-2804
- Size(s):
- p. 2783-2804
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT JWST has revealed a large population of accreting black holes (BHs) in the early Universe. Recent work has shown that even after accounting for possible systematic biases, the high-z$$M_*{\!-\!}M_{\rm \rm bh}$$ relation can be above the local scaling relation by $$\gt 3\sigma$$. To understand the implications of these overmassive high-z BHs, we study the BH growth at $$z\sim 4{\!-\!}7$$ using the $$[18~\mathrm{Mpc}]^3$$BRAHMA cosmological simulations with systematic variations of heavy seed models that emulate direct collapse black hole (DCBH) formation. In our least restrictive seed model, we place $$\sim 10^5~{\rm M}_{\odot }$$ seeds in haloes with sufficient dense and metal-poor gas. To model conditions for direct collapse, we impose additional criteria based on a minimum Lyman Werner flux (LW flux $$=10~J_{21}$$), maximum gas spin, and an environmental richness criterion. The high-z BH growth in our simulations is merger dominated, with a relatively small contribution from gas accretion. The simulation that includes all the above seeding criteria fails to reproduce an overmassive high-z$$M_*{\!-\!}M_{\rm bh}$$ relation consistent with observations (by factor of $$\sim 10$$ at $$z\sim 4$$). However, more optimistic models that exclude the spin and environment based criteria are able to reproduce the observed relations if we assume $$\lesssim 750~\mathrm{Myr}$$ delay times between host galaxy mergers and subsequent BH mergers. Overall, our results suggest that current JWST observations may be explained with heavy seeding channels if their formation is more efficient than currently assumed DCBH conditions. Alternatively, we may need higher initial seed masses, additional contributions from lighter seeds to BH mergers, and / or more efficient modes for BH accretion.more » « less
-
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
-
ABSTRACT We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $$z$$ ≳ 6 quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in the IllustrisTNG simulations. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $$\gtrsim 10^{12}\, {\rm M}_{\odot }\,h^{-1}$$ haloes by $$z$$ = 7) using a (9 Mpc h−1)3 simulated volume. BH growth is maximal within haloes that are compact and have a low tidal field. For these haloes, we consider an array of gas-based seeding prescriptions wherein $$M_{\mathrm{seed}}=10^4\!-\!10^6\, {\rm M}_{\odot }\,h^{-1}$$ seeds are inserted in haloes above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $$\tilde{M}_{\mathrm{h}}$$ and $$\tilde{M}_{\mathrm{sf,mp}}$$, respectively, in units of Mseed). We find that a seed model with $$\tilde{M}_{\mathrm{sf,mp}}=5$$ and $$\tilde{M}_{\mathrm{h}}=3000$$ successfully produces a $$z$$ ∼ 6 quasar with $$\sim 10^9\, {\rm M}_{\odot }$$ mass and ∼1047 erg s−1 luminosity. BH mergers play a crucial role at $$z$$ ≳ 9, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. With more stringent seeding conditions (e.g. $$\tilde{M}_{\mathrm{sf,mp}}=1000$$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $$z$$ ≳ 6 quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors ≳10) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results demonstrate that progenitors of $$z$$ ∼ 6 quasars have distinct BH merger histories for different seeding models, which will be distinguishable with Laser Interferometer Space Antenna observations.more » « less
-
null (Ed.)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
