skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Cosmic Evolution of the Supermassive Black Hole Population: A Hybrid Observed Accretion and Simulated Mergers Approach
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
Award ID(s):
2407089 2106990
PAR ID:
10593035
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
976
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
6
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. ABSTRACT The concurrent growth of supermassive black holes (SMBHs) and their host galaxies remains to be fully explored, especially at high redshift. While often understood as a consequence of self-regulation via AGN feedback, it can also be explained by alternative SMBH accretion models. Here, we expand on previous work by studying the growth of SMBHs with the help of a large suite of cosmological zoom-in simulations (MassiveFIRE) that are part of the Feedback in Realistic Environments (FIRE) project. The growth of SMBHs is modelled in post-processing with different black hole accretion models, placements, and merger treatments, and validated by comparing to on-the-fly calculations. Scaling relations predicted by the gravitational torque-driven accretion (GTDA) model agree with observations at low redshift without the need for AGN feedback, in contrast to models in which the accretion rate depends strongly on SMBH mass. At high redshift, we find deviations from the local scaling relations in line with previous theoretical results. In particular, SMBHs are undermassive, presumably due to stellar feedback, but start to grow efficiently once their host galaxies reach M* ∼ 1010M⊙. We analyse and explain these findings in the context of a simple analytic model. Finally, we show that the predicted scaling relations depend sensitively on the SMBH location and the efficiency of SMBH merging, particularly in low-mass systems. These findings highlight the relevance of understanding the evolution of SMBH-galaxy scaling relations to predict the rate of gravitational wave signals from SMBH mergers across cosmic history. 
    more » « less
  3. Abstract We analyze the dynamics of low-mass black hole (BH) seeds in the high-redshift (z ≳ 5) Universe using a suite of [4.5 Mpc]3and [9 Mpc]3BRAHMAcosmological hydrodynamic simulations. The simulations form seeds with massMseed = 2.2 × 103Min halos that exceed critical thresholds of dense and metal-poor gas mass (5–150Mseed) and the halo mass (1000–10,000Mseed). While the initialBRAHMAboxes pinned the BHs to the halo centers, here we implement a subgrid dynamical friction (DF) model. We also compare simulations where the BH is allowed to wander without the added DF. We investigate the spatial and velocity offsets of BHs in their host subhalos, as well as BH merger rates. We find that subgrid DF is crucial to ensure that a significant fraction of BHs effectively sink to halo centers byz ∼ 5, thereby enabling them to get gravitationally bound and merge with other BHs at separations close to the spatial resolution (∼0.2–0.4 kpc) of the simulation. For the BHs that merge, the associated merger timescales lag between ∼100 and 1000 Myr after their host halos merge. Compared to predictions using BH repositioning, the overallz ≳ 5 BH merger rates under subgrid DF decrease by a factor of ∼4–10. Under subgrid DF, the different seed models predict merger rates between ∼100 and 1000 events per year atz ≳ 5. These mergers dominate early BH growth, assembling BHs up to ∼104–105Mbyz ∼ 5, wherein ≲2% of their mass is assembled via gas accretion. Our results highlight the promise for constraining seeding mechanisms using gravitational waves from future facilities such as the Laser Interferometer Space Antenna. 
    more » « less
  4. ABSTRACT Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies, however, remain open questions. Recent galaxy formation simulations predict a delayed, or ‘two-phase,’ growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs on to the local scaling relations. The predicted two-phase growth implies a steep increase, or ‘kink,’ in BH-galaxy scaling relations at a stellar mass $$\rm {M}_{*}\sim 5\times 10^{10}$$ M⊙. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at z ∼ 0.5−4. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and active galactic nuclei activity at the low-mass end, which could more directly test two-phase SMBH growth. 
    more » « less
  5. Abstract From the luminous quasars atz∼ 6 to the recentz∼ 9–11 active galactic nuclei (AGN) revealed by JWST, observations of the earliest black hole (BH) populations can provide unique constraints on BH evolution. We use theBRAHMAsimulations with constrained initial conditions to investigate BH assembly in extreme overdense regions. The simulations implement heavy ∼104–105Mseeds forming in dense, metal-poor gas exposed to sufficient Lyman–Werner flux. With gas accretion modeled via the Bondi–Hoyle formalism and BH dynamics with a subgrid dynamical friction scheme, we isolate the impact of seeding, dynamics, accretion, and feedback on BH evolution. With fiducial stellar and AGN feedback inherited fromIllustrisTNG, accretion is suppressed atz≳ 9, leaving mergers as the dominant growth channel. Gas accretion dominates atz≲ 9, where permissive models (super-Eddington or low radiative efficiency) build ∼109MBHs powering quasars byz∼ 6, while stricterIllustrisTNG-based prescriptions yield much smaller BHs (∼106–108M). Our seed models strongly affect mergers atz≳ 9: only the most lenient models (with ∼105Mseeds) produce enough BH mergers to reach ≳106Mbyz∼ 10, consistent with current estimates for GN-z11. Our dynamical friction model gives low merger efficiencies. Therefore, even in such extreme regions, we are unable to produce ≳107MBHs byz∼ 9–10, as currently inferred for GHZ9, UHZ1, and CAPERS-LRD-z9. If the BH-to-stellar mass ratios of these sources are indeed so extreme, they would require either very short BH merger timescales or reduced AGN thermal feedback. Weaker stellar feedback boosts both star formation and BH accretion and cannot raise these ratios. 
    more » « less