This content will become publicly available on February 3, 2027

Title: Tracking the Assembly of Supermassive Black Holes: A Comparison of Diverse Models across Cosmic Time
Abstract 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
Award ID(s):
2502826
PAR ID:
10678483
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
American Astronomical Society
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
998
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
48
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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≳ 108MorM≳ 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
  2. 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 ∼1014Mhalos. 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
  3. Abstract To facilitate new studies of galaxy-merger-driven fueling of active galactic nuclei (AGNs), we present a catalog of 387 AGNs that we have identified in the final population of over 10,000z< 0.15 galaxies observed by the Sloan Digital Sky Survey-IV (SDSS-IV) integral field spectroscopy survey Mapping Nearby Galaxies at Apache Point Observatory (MaNGA). We selected the AGNs via mid-infrared Wide-field Infrared Survey Explorer colors, Swift/Burst Alert Telescope ultra-hard X-ray detections, NRAO Very Large Array Sky Survey and Faint Images of the Radio Sky at Twenty centimeters radio observations, and broad emission lines in SDSS spectra. By combining the MaNGA AGN catalog with a new SDSS catalog of galaxy mergers that were identified based on a suite of hydrodynamical simulations of merging galaxies, we study the link between galaxy mergers and nuclear activity for AGNs above a limiting bolometric luminosity of 1044.4erg s−1. We find an excess of AGNs in mergers, relative to nonmergers, for galaxies with stellar mass ∼1011M, where the AGN excess is somewhat stronger in major mergers than in minor mergers. Further, when we combine minor and major mergers and sort by merger stage, we find that the highest AGN excess occurs in post-coalescence mergers in the highest-mass galaxies. However, we find no evidence of a correlation between galaxy mergers and AGN luminosity or accretion rate. In summary, while galaxy mergers overall do appear to trigger or enhance AGN activity more than nonmergers, they do not seem to induce higher levels of accretion or higher luminosities. We provide the MaNGA AGN Catalog and the MaNGA Galaxy Merger Catalog for the community here. 
    more » « less
  4. Abstract The James Webb Space Telescope has unveiled an abundant population of potential active galactic nuclei (AGN) at high redshift (z ≳ 4) known as little red dots (LRDs), which are likely hosted in relatively low-mass galaxies. However, previous theoretical models have highlighted the difficulty in continuously feeding massive black holes in the central regions of bursty, high-redshift galaxies because of repeated gas evacuation by stellar feedback. We analyze galaxies in high-redshift FIRE-2 simulations to understand whether they are capable of hosting the observed abundant population of high-redshift AGN. We use a gravitational torque-driven accretion (GTDA) model and a simple freefall accretion model to derive black hole accretion rates and construct predicted AGN bolometric luminosity functions forz = 5–7. The GTDA model and the freefall model with black holes accreting ≲1% of their central gas supply (<100 pc) per freefall time predict AGN abundances that are more than sufficient to explain the most recent LRD observations. The fiducial models, in fact, overpredict the number of low-luminosity AGN as compared with observations. We explore possible resolutions of this tension. A plausible, though likely not unique, scenario for alleviating the AGN overpredictions, which also provides a good match to the host-galaxy UV luminosity distribution, suggests that LRDs are super-Eddington-accreting, Eddington-luminosity-limited,MBH ≳ 2 × 105Mblack holes residing inM ≳ 2 × 107Mgalaxies. We show that, under simple assumptions, mock observations of such sources can reproduce key observed LRD characteristics. 
    more » « less
  5. Abstract Active galactic nuclei (AGNs) are a promising source of the binary black hole mergers observed in gravitational waves with LIGO-Virgo-Kagra (LVK). Constraining the AGN channel allows us to limit AGN parameter space (disk density, size, average lifetime) and nuclear star cluster (NSC) parameter space. Constraints on AGNs and NSCs have implications for Λ cold dark matter models of AGN feedback and models of AGN-driven supermassive black hole merger and growth. Here we present several qualitative studies of the AGN channel using new public, open-source, fast, reproducible codeMcFACTS(https://github.com/mcfacts/mcfacts): Monte Carlo for AGN channel Testing and Simulation. We demonstrate several important features for testing the AGN channel, including: (i) growth to large-mass intermediate-mass black hole is helped by the presence of migration traps or swamps, (ii) flat BH initial mass functions highlight hierarchical merger features in the mass spectrum, (iii) the (q, χeff) anticorrelation is a strong test of the bias to prograde mergers in the AGN channel, (iv) spheroid encounters can drive a fraction of mergers with high in-plane spin components (χp), (v) a high rate of extreme mass ratio inspirals are driven by an initial population of embedded retrograde BH, and (vi) both LVK and LISA are powerful probes of models of AGN disks and their embedded populations. 
    more » « less