🔭 astrophysics

AGN-driven BBH mergers: Black hole populations and hierarchical growth across the AGN parameter space

This paper presents an updated semi-analytical framework demonstrating that active galactic nucleus (AGN) disks can efficiently drive hierarchical black hole mergers to form high-mass binaries consistent with events like GW190521 and GW231123, with merger rates and population characteristics strongly dependent on AGN parameters such as viscosity, Eddington ratio, and supermassive black hole mass.

Original authors: Maria Paola Vaccaro, Michela Mapelli, Alessandro Alberto Trani, Boyuan Liu

Published 2026-06-10

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Maria Paola Vaccaro, Michela Mapelli, Alessandro Alberto Trani, Boyuan Liu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the center of a galaxy as a giant, swirling cosmic dance floor. In the middle sits a supermassive black hole, the "bouncer" of the party, surrounded by a thick disk of gas and dust. This is an Active Galactic Nucleus (AGN).

This paper is like a detailed simulation of what happens when smaller black holes (the "dancers") get stuck in this gas disk. The authors, a team of astrophysicists, built a computer model to see how these black holes find partners, dance together, and eventually crash into one another to create even bigger black holes.

Here is the story of their findings, broken down into simple concepts:

1. The Dance Floor Rules (The Gas Disk)

Think of the gas disk as a thick, sticky syrup. When a black hole moves through it, the gas pushes on it, causing it to migrate (move) toward the center.

The Viscosity Factor: The authors tested two types of "syrup."
- Thin Syrup (Low Viscosity): The black holes move easily and quickly. They find partners fast, merge, and the resulting new black hole stays on the dance floor to find another partner. This leads to a "hierarchical" growth, like a snowball rolling down a hill, getting bigger and bigger.
- Thick Syrup (High Viscosity): The gas is too sticky and chaotic. It disrupts the black holes, making it hard for them to stay together long enough to merge repeatedly. The "snowball" effect stops early, and the black holes remain smaller.

2. The "Snowball" Effect (Hierarchical Growth)

In the low-viscosity (thin syrup) scenario, the black holes can merge over and over again.

First Generation: Two normal-sized black holes merge.
Second Generation: The new, bigger black hole merges with another one.
Result: This creates a "high-mass tail." The paper finds that this process can build black holes so massive they shouldn't exist according to standard star evolution rules (crossing the "pair-instability mass gap"). It's like building a tower of blocks that keeps getting taller than physics usually allows, but the gas disk provides the extra support needed.

3. The Dance Partners (Mass and Spin)

The paper looks at who dances with whom:

Uneven Couples: As the black holes get bigger through repeated mergers, they tend to pair up with much smaller partners. It's like a giant trying to dance with a child; the mass ratio becomes very unequal.
Spinning in Sync: Because the gas disk forces the black holes to spin in the same direction as the disk itself, their "spins" (how fast they rotate) tend to line up. If they spin the same way, it creates a specific signature that ground-based detectors (like LIGO) can spot.
The Eccentricity: Most of these pairs become perfectly round circles before they merge (low eccentricity), but a few remain in weird, stretched-out orbits (high eccentricity) because they merged so quickly they didn't have time to "smooth out."

4. The "Big Hits" (GW190521 and GW231123)

The authors compared their simulation to real-life events detected by gravitational wave observatories, specifically two very massive collisions known as GW190521 and GW231123.

The Verdict: The AGN disk model can produce black holes heavy enough to match these events. However, it's rare.
The Catch: While the model can create the right total mass, it struggles to perfectly match the specific combination of the two individual masses and their spins for these specific events. It suggests that while AGN disks are a likely candidate for these massive crashes, they might not be the only explanation, or our understanding of the "syrup" (the gas physics) needs fine-tuning.

5. The "Time Limit" (AGN Lifetime)

The dance floor doesn't last forever. The central black hole eventually stops feeding, and the gas disk dissipates.

The paper finds that the longer the party lasts (the AGN lifetime), the more time black holes have to merge repeatedly. If the party is short, the "snowball" never gets very big. If the party lasts a long time, you get the massive, intermediate-mass black holes.

Summary of the "Recipe"

The paper concludes that for AGN disks to be efficient factories for massive black hole mergers, three things are crucial:

Low Viscosity: The gas needs to be "thin" enough to let black holes migrate and merge efficiently.
Gas Hardening: The gas must actively push the binary pairs closer together. If you turn off the gas effect and rely only on random bumps into other stars, the massive black holes don't form.
Time: The AGN needs to stay active long enough for the chain of mergers to happen.

In short, the paper argues that active galactic nuclei are powerful cosmic workshops where gas acts as a glue and a conveyor belt, allowing black holes to build up massive sizes through repeated collisions, creating a unique population of mergers that we can now start to detect.

Technical Summary: AGN-driven BBH mergers: Black hole populations and hierarchical growth across the AGN parameter space

Problem and Motivation
Active galactic nuclei (AGN) disks have emerged as a promising dynamical channel for the formation and merger of binary black holes (BBHs). While stellar-mass black holes (BHs) embedded in these disks can migrate, pair up, and merge via gas torques and gravitational wave (GW) emission, the predicted properties of this population remain highly uncertain. Previous studies have highlighted a strong dependence on physical ingredients such as disk structure, migration torques, gas hardening efficiency, and AGN lifetimes, often adopting different assumptions that hinder direct comparison. Recent observational constraints suggest the AGN channel cannot account for the entirety of the BBH merger population detected by GW observatories, yet its specific contribution and distinctive signatures require a systematic exploration of the AGN parameter space to identify robust observational features and sources of uncertainty.

Methodology
The authors present an updated semi-analytical framework implemented within the fastcluster software environment to model the formation, evolution, and hierarchical growth of BBHs in AGN disks. The approach tracks the lifecycle of stellar-origin (first-generation, 1g) black holes through capture, radial migration, binary pairing, gas-driven hardening, dynamical perturbations, and eventual merger.

Key methodological components include:

Disk Modeling: The accretion disk is modeled using the steady-state analytic Sirko & Goodman (2003) (SG) model. The authors self-consistently solve one-dimensional disk equations for varying supermassive black hole (SMBH) masses ( $M_\bullet$ ), Eddington ratios ( $\lambda_\bullet$ ), and disk viscosity parameters ( $\alpha$ ). Two viscosity values are tested: $\alpha = 0.01$ (low viscosity) and $\alpha = 0.1$ (high viscosity).
Initial Conditions: Simulations initialize BH populations based on nuclear star cluster (NSC) scaling relations. 1g BHs are drawn from a population synthesis catalog (sevn) assuming solar metallicity. Spins are assigned based on a Maxwellian distribution, with orientations aligned to the disk angular momentum.
Evolutionary Physics:
- Migration: BHs migrate via Type I and Type II regimes, with prescriptions from Grishin et al. (2024) and Kanagawa et al. (2018).
- Pairing: Binary formation follows a gas-assisted criterion, accounting for the fraction of prograde orbits ( $f_{progr}$ ).
- Hardening and Encounters: Post-formation, binaries evolve under combined gas torques (circumbinary disk interaction) and three-body encounters with single BHs. The authors integrate orbital evolution equations combining gas-driven hardening (using Ishibashi & Gröbner 2024) and GW emission (Peters 1964). Encounters are modeled using a pre-computed grid of post-Newtonian scattering experiments (tsunami).
- Hierarchical Growth: Merger remnants are calculated using aligned-spin prescriptions and recoil kick velocities. Remnants that remain bound to the disk participate in subsequent merger generations ( $N$ -th generation), subject to AGN lifetime limits and ejection criteria.
Population Synthesis: To construct an intrinsic BBH population representative of the local Universe ( $z \lesssim 0.3$ ), the authors weight simulation outputs according to observed AGN distributions of $M_\bullet$ and $\lambda_\bullet$ (Greene & Ho 2007).

Key Results
The study systematically explores the dependence of BBH merger populations on $M_\bullet$ , $\lambda_\bullet$ , and $\alpha$ , yielding the following findings:

Hierarchical Growth and Mass Distribution:
- AGN disks can produce repeated mergers, building a high-mass tail extending beyond the pair-instability supernova (PISN) mass gap ( $50 \lesssim m_1/M_\odot \lesssim 130$ ) into the intermediate-mass black hole (IMBH) range.
- Viscosity Dependence: Hierarchical growth is significantly more efficient in low-viscosity disks ( $\alpha = 0.01$ ), where higher gas surface densities shorten dynamical timescales, allowing for higher merger generations ( $N_{gen} > 100$ ). High-viscosity disks ( $\alpha = 0.1$ ) suppress the formation of massive remnants, typically saturating at lower generations ( $N_{gen} \approx 20$ ).
- Mass Ratios: The population evolves toward increasingly unequal mass ratios at high primary mass. While primary masses grow through successive mergers, secondary masses saturate at lower values (typically below the PISN gap), as high-generation primaries merge with lower-generation secondaries.
Merger Efficiency and Parameter Space:
- Merger efficiency ( $\eta$ ) generally increases with the Eddington ratio ( $\lambda_\bullet$ ).
- The dependence on $M_\bullet$ is non-trivial; efficiency and maximum remnant mass tend to peak at intermediate-to-large SMBH masses, but no single characteristic behavior spans the entire parameter space. Some high-mass configurations show no mergers, indicating that large $M_\bullet$ alone does not guarantee efficient assembly.
Spin and Eccentricity Signatures:
- Effective Spin ( $\chi_{eff}$ ): The distribution of $\chi_{eff}$ depends strongly on the initial prograde fraction ( $f_{progr}$ ). Prograde configurations yield positive $\chi_{eff}$ , while retrograde yields negative values. Intermediate fractions produce a bimodal distribution. A correlation exists between primary mass and $|\chi_{eff}|$ , as repeated mergers reinforce coherent spin orientations.
- Eccentricity: Most mergers enter the ground-based detector band (at $f_{peak} = 10$ Hz) with low eccentricities ( $e \sim 10^{-3}–10^{-2}$ ). However, a subpopulation retains high eccentricity ( $e \gtrsim 0.8$ ), having merged before significant circularization occurred.
Comparison to GW Events:
- The AGN channel can reproduce systems consistent with the massive events GW190521 and GW231123, though these occupy the tails of the predicted distributions.
- Systems matching the mass parameters of GW190521 occur at a rate of $\sim 2 \times 10^{-3}$ in low-viscosity models, while GW231123-like systems are rarer ( $\sim 3 \times 10^{-4}$ ).
- The model tends to underpredict the secondary mass ( $m_2$ ) for these specific events compared to inferred parameters, suggesting that while the AGN channel can produce the primary mass and spin characteristics, matching the full component mass distribution is challenging.
Model Robustness:
- The general properties of the population (high-mass tail, spin correlations) are robust across variations in migration torque formalisms (e.g., thermal corrections) and gas hardening prescriptions (IG20 vs. IG24 vs. C24).
- Critical Dependency: The formation of massive remnants and the high-mass tail are entirely suppressed when gas hardening is switched off, confirming that gas-driven binary evolution is a key ingredient for hierarchical growth in AGN disks.

Significance and Claims
The paper claims to provide a comprehensive and up-to-date framework that systematically maps the AGN parameter space to assess the robustness of BBH merger predictions. The authors assert that the AGN channel does not produce a single characteristic population but rather a broad diversity of outcomes. The study identifies distinctive population-level signatures—specifically the correlation between primary mass and $|\chi_{eff}|$ , the trend toward unequal mass ratios at high masses, and the retention of memory regarding initial orbital configurations in the spin distribution.

The work emphasizes that while the AGN channel can explain massive BBH events like GW190521 and GW231123, these events likely represent rare occurrences within the channel. The authors conclude that gas-driven binary evolution is essential for efficient hierarchical growth, and that the specific physical assumptions regarding disk structure and hardening significantly influence the efficiency of massive remnant formation, though the qualitative existence of a high-mass tail remains a robust feature of the channel.