Direct collapse black holes (BHs) are promising candidates for producing massive z ≳ 6 quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: (1) low gas angular momentum (spin), and (2) a minimum incident Lyman–Werner (LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4\!-\!10^6\, \mathrm{M}_{\odot }\, h^{-1})$ in haloes with a total mass >3000 × Mseed and a dense, metal-poor gas mass >5 × Mseed. Within this framework, we find that the seed-forming haloes have a prior history of star formation and metal enrichment, but they also contain pockets of dense, metal-poor gas. When seeding is further restricted to haloes with low gas spins, the number of seeds formed is suppressed by factors of ∼6 compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal-poor gas is required to have a critical LW flux (Jcrit). Even for Jcrit values as low as 50J21, no $8\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$ seeds are formed. While lower mass ($1.25\times 10^{4},1\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$) seeds do form, they are strongly suppressed (by factors of ∼10–100) compared to the baseline model at gas mass resolutions of $\sim 10^4~\mathrm{M}_{\odot }\, h^{-1}$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, none of the seeds are able to grow to the supermassive regime ($\gtrsim 10^6~\mathrm{M}_{\odot }\, h^{-1}$) by z = 7. Our results hint that producing the bulk of the z ≳ 6 supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities.
Formation of supermassive black holes (BHs) remains a theoretical challenge. In many models, especially beginning from stellar relic ‘seeds,’ this requires sustained super-Eddington accretion. While studies have shown BHs can violate the Eddington limit on accretion disc scales given sufficient ‘fuelling’ from larger scales, what remains unclear is whether or not BHs can actually capture sufficient gas from their surrounding interstellar medium (ISM). We explore this in a suite of multiphysics high-resolution simulations of BH growth in magnetized, star-forming dense gas complexes including dynamical stellar feedback from radiation, stellar mass-loss, and supernovae, exploring populations of seeds with masses $\sim 1\!-\!10^{4}\, \mathrm{M}_{\odot }$. In this initial study, we neglect feedback from the BHs: so this sets a strong upper limit to the accretion rates seeds can sustain. We show that stellar feedback plays a key role. Complexes with gravitational pressure/surface density below $\sim 10^{3}\, \mathrm{M}_{\odot }\, {\rm pc^{-2}}$ are disrupted with low star formation efficiencies so provide poor environments for BH growth. But in denser cloud complexes, early stellar feedback does not rapidly destroy the clouds but does generate strong shocks and dense clumps, allowing $\sim 1{{\ \rm per\ cent}}$ of randomly initialized seeds to encounter a dense clump with low relative velocity and produce runaway, hyper-Eddington accretion (growing by orders of magnitude). Remarkably, mass growth under these conditions is almost independent of initial BH mass, allowing rapid intermediate-mass black hole (IMBH) formation even for stellar-mass seeds. This defines a necessary (but perhaps not sufficient) set of criteria for runaway BH growth: we provide analytic estimates for the probability of runaway growth under different ISM conditions.
more » « less- NSF-PAR ID:
- 10383259
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 518
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3606-3621
- Size(s):
- p. 3606-3621
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
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
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 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.
-
ABSTRACT Possible formation scenarios of supermassive black holes (BHs) in the early universe include rapid growth from less massive seed BHs via super-Eddington accretion or runaway mergers, yet both of these scenarios would require seed BHs to efficiently sink to and be trapped in the Galactic Centre via dynamical friction. This may not be true for their complicated dynamics in clumpy high-z galaxies. In this work, we study this ‘sinking problem’ with state-of-the-art high-resolution cosmological simulations, combined with both direct N-body integration of seed BH trajectories and post-processing of randomly generated test particles with a newly developed dynamical friction estimator. We find that seed BHs less massive than $10^8\, \mathrm{M}_\odot$ (i.e. all but the already-supermassive seeds) cannot efficiently sink in typical high-z galaxies. We also discuss two possible solutions: dramatically increasing the number of seeds such that one seed can end up trapped in the Galactic Centre by chance, or seed BHs being embedded in dense structures (e.g. star clusters) with effective masses above the mass threshold. We discuss the limitations of both solutions.more » « less
