Abstract The accretion and feedback processes governing supermassive black hole (SMBH) growth span an enormous range of spatial scales, from the Event Horizon to the circumgalactic medium. Recent general relativistic magnetohydrodynamic (GRMHD) simulations demonstrate that strong magnetic fields can substantially suppress gas accretion onto black holes. These simulations show that magnetic fields create magnetically arrested disk states, reducing inflow rates by up to 2 orders of magnitude relative to classical predictions. We incorporate this magnetic suppression prescription from recent GRMHD studies into DarkSage, a semianalytic model that tracks SMBH and galaxy coevolution over cosmic time. Implementing the suppression across different accretion rate regimes, we explore its impact on the distribution of black hole masses, stellar masses in galaxies, and active galactic nucleus (AGN) luminosities. We find that restricting suppression to sub-Eddington accretors (fEdd < 3 × 10−3) and rescaling AGN feedback efficiencies gives simultaneous agreement with the observed local distributions of both galaxy and black hole masses. At early cosmic times (z > 6), super-Eddington growth episodes dominate in our model, reproducing the high number densities of luminous AGN recently discovered by the James Webb Space Telescope. Our results highlight the critical sensitivity of galaxy assembly to the coupling between small-scale accretion physics and large-scale feedback regulation. Magnetic suppression of hot gas accretion can reconcile low-redshift constraints while preserving the rapid black hole growth required at early cosmic epochs, thereby providing a physically motivated bridge between horizon-scale GRMHD simulations and cosmological galaxy-formation models.
more »
« less
Bridging Scales: How Much Do Supermassive Black Holes Grow in the Suppressed Bondi Regime?
Abstract The coevolution of supermassive black holes (SMBHs) and their host galaxies remains one of the central open questions in cosmology, rooted in the coupling between accretion, feedback, and the multiscale physics that links the event horizon to the circumgalactic medium. Here we bridge these scales by embedding a first-principles, GRMHD-informed prescription for black hole accretion and feedback—derived from multizone simulations that self-consistently connect inflows and outflows from the horizon to the Bondi radius—within cosmological magnetohydrodynamic zoom-in simulations of ∼1014M⊙halos. These GRMHD results predict a “suppressed Bondi” regime in which magnetic stresses and relativistic winds strongly reduce effective accretion rates in a spin-dependent manner. We find that black holes cannot grow efficiently by accretion until they exceed ∼107M⊙, regardless of the feedback strength. Beyond this threshold, systems bifurcate: low-spin (η ∼ 0.02) black holes continue to accrete without quenching star formation, while high-spin (η ≳ 0.3) black holes quench effectively but become starved of further growth. Early, massive seeding partially alleviates this tension through merger-driven assembly, yet an additional cold or super-Eddington accretion mode appears essential to reproduce the observed SMBH population and the empirical black hole–galaxy scaling relations. Our results demonstrate that GRMHD-informed feedback models can account for the maintenance-mode behavior of low-luminosity active galactic nuclei like M87*, but cannot by themselves explain the full buildup of SMBH mass across cosmic time. A unified, multiregime framework is required to capture the evolving interplay between spin-dependent feedback, cold inflows, and mergers in driving coevolution.
more »
« less
- Award ID(s):
- 2502826
- PAR ID:
- 10678484
- Publisher / Repository:
- American Astronomical Society
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 998
- Issue:
- 1
- ISSN:
- 2041-8205
- Page Range / eLocation ID:
- L18
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Tracking the Assembly of Supermassive Black Holes: A Comparison of Diverse Models across Cosmic TimeAbstract Galaxies grow alongside central supermassive black holes (SMBHs) through fueling and feedback. However, the origins of this coevolution remain unclear and vary across modeling frameworks. Using semianalytic models (SAMs), we trace SMBH mass assembly acrossMBH ∼ 106−10M⊙. We find significant discrepancies between observations and physics-based models of the local black hole mass function (BHMF), likely from differences in the stellar mass function and scaling relations used to infer the BHMF. Most physics-based models agree atz ∼ 1–4 and broadly match the JWST broad-line active galactic nucleus (AGN) BHMFs atz = 4–5. These models also reproduce the observed bolometric AGN luminosity evolution, except the SAMDark Sage, which predicts an excess. Interestingly, this pronounced “knee” in the bolometric AGN luminosity function predicted byDark SagearoundLbol ∼ 1046erg s−1is consistent with the inferred abundance and luminosity of “little red dots” atz = 5–6, under the assumption that they are powered entirely by AGN activity. In contrast to other models,Dark Sagedeploys multiple growth channels for SMBHs that include mergers, hot-mode accretion, merger-driven cold-accretion, and secular-instability-driven accretion. We analyze the black hole mass buildup and accretion histories inDark Sage, which, unlike other models, also allows for super-Eddington accretion, and we find that, on average, SMBHs primarily grow through secular disk instabilities and merger-driven cold gas accretion modes. We also find that black hole mergers contribute the majority of the growth of ∼60% of the total mass budget only for the most massive SMBHs byz= 0.more » « less
-
Abstract Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofRB≈ 2 × 105GM•/c2, whereM•is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r−1rather thanr−3/2and the accretion rate is consequently suppressed by over 2 orders of magnitude below the Bondi rate . We find continuous energy feedback from the accretion flow to the external medium at a level of . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.more » « less
-
Abstract General relativistic magnetohydrodynamics (GRMHD) simulations are an indispensable tool in studying accretion onto compact objects. The Event Horizon Telescope (EHT) frequently uses libraries of ideal GRMHD simulations to interpret polarimetric, event-horizon-scale observations of supermassive black holes at the centers of galaxies. In this work, we present a library of 10 nonradiative, ideal GRMHD simulations that were utilized by the EHT Collaboration in their analysis of Sagittarius A*. The parameter survey explores both low (SANE) and high (MAD) magnetization states across five black hole spinsa* = −15/16, −1/2, 0, +1/2, +15/16 where each simulation was run out to 30,000GM/c−3. We find the angular momentum and energy flux in SANE simulations closely matches the thin-disk value, with minor deviations in prograde models due to fluid forces. This leads to spin equilibrium arounda* ∼ 0.94, consistent with previous studies. We study the flow of conserved quantities in our simulations and find mass, angular momentum, and energy transport in SANE accretion flows to be primarily inward and fluid dominated. MAD models produce powerful jets with outflow efficiency >1 fora* = + 0.94, leading to black hole spin-down in prograde cases. We observe outward directed energy and angular momentum fluxes on the horizon, as expected for the Blandford–Znajek mechanism. MAD accretion flows are sub-Keplerian and exhibit greater variability than their SANE counterpart. They are also hotter than SANE disks withinr ≲ 10GM/c−2. This study is accompanied by a public release of simulation data athttp://thz.astro.illinois.edu/.more » « less
-
Abstract Supermassive black holes (SMBHs) can grow through both accretion and mergers. It is still unclear how SMBHs evolve under these two channels from high redshifts to the SMBH population we observe in the local Universe. Observations can directly constrain the accretion channel but cannot effectively constrain mergers yet, while cosmological simulations provide galaxy merger information but can hardly return accretion properties consistent with observations. In this work, we combine the observed accretion channel and the simulated merger channel, taking advantage of observations and cosmological simulations, to depict a realistic evolution pattern of the SMBH population. With this methodology, we can derive the scaling relation between the black hole mass (MBH) and host-galaxy stellar mass (M⋆), and the local black hole mass function (BHMF). Our scaling relation is lower than those based on dynamically measuredMBH, supporting the claim that dynamically measured SMBH samples may be biased. We show that the scaling relation has little redshift evolution. The BHMF steadily increases fromz= 4 toz= 1 and remains largely unchanged fromz= 1 toz= 0. The overall SMBH growth is generally dominated by the accretion channel, with possible exceptions at high mass (MBH≳ 108M⊙orM⋆≳ 1011M⊙) and low redshift (z≲ 1). We also predict that around 25% of the total SMBH mass budget in the local Universe may be locked within long-lived, wandering SMBHs, and the wandering mass fraction and wandering SMBH counts increase withM⋆.more » « less
An official website of the United States government

