🔭 astrophysics

McFACTS II: Mass Ratio--Effective Spin Relationship of Black Hole Mergers in the AGN Channel

Using the McFACTS simulation code, this study demonstrates that dense, moderately short-lived AGN disks with a steep black hole mass function and a high fraction of prograde binaries are required to reproduce the observed anti-correlation between mass ratio and effective spin in LIGO-Virgo-KAGRA gravitational wave events.

Original authors: Harrison E. Cook, Barry McKernan, K. E. Saavik Ford, Vera Delfavero, Kaila Nathaniel, Jake Postiglione, Shawn Ray, Emily J. McPike, Richard O'Shaughnessy

Published 2026-04-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Harrison E. Cook, Barry McKernan, K. E. Saavik Ford, Vera Delfavero, Kaila Nathaniel, Jake Postiglione, Shawn Ray, Emily J. McPike, Richard O'Shaughnessy

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 massive, swirling cosmic dance floor. In the middle sits a supermassive black hole, the "DJ," spinning a giant disk of gas around it. This is an Active Galactic Nucleus (AGN). Around this dance floor, smaller black holes are floating like dancers waiting for a partner.

This paper is a simulation study using a computer program called McFACTS (which stands for Monte Carlo For AGN Channel Testing and Simulation). The researchers wanted to figure out: If black holes merge inside these galactic dance floors, what do the resulting "couples" look like?

Specifically, they are looking at two things:

Mass Ratio ( $q$ ): How similar in size are the two black holes? (Are they a "power couple" of similar size, or a giant with a tiny partner?)
Effective Spin ( $\chi_{eff}$ ): How fast are they spinning, and are they spinning in the same direction as the dance floor?

Here is the breakdown of their findings using everyday analogies:

1. The "Anti-Correlation" Mystery

The researchers noticed a pattern in real data from gravitational wave detectors (LIGO): There seems to be a rule where smaller partners (low mass ratio) tend to spin faster and in alignment (high spin), while similar-sized partners spin slower or randomly.

Think of it like a figure skater. If a heavy skater grabs a tiny partner, the heavy skater's momentum dominates, and they spin together in a very specific, fast way. If two skaters are the same size, they might struggle to agree on a spin direction, resulting in a messy, slower spin.

The paper asks: Can the AGN dance floor explain this pattern?

2. The Dance Floor Matters (Disk Models)

The team tested two different types of "dance floors" (gas disk models):

The "Sirko & Goodman" (SG) Floor: This is a dense, crowded, and short-lived party. It's chaotic but efficient.
The "Thompson et al." (TQM) Floor: This is a sparse, quiet, and long-lived gathering.

The Result: The dense, short-lived SG floor is the winner. It successfully creates the "anti-correlation" pattern seen in real data. The sparse TQM floor is too quiet; the black holes drift apart too much and don't merge often enough to create the pattern.

Analogy: You need a crowded mosh pit (SG) to get people bumping into each other and forming couples quickly. A sparse park (TQM) just isn't busy enough.

3. The "Generational" Effect

The simulation tracks "generations" of black holes:

1st Generation (1g): Black holes born from dying stars. They enter the dance floor with random spins.
2nd/3rd Generation (2g/3g): Black holes that were already the result of a previous merger.

The Result:

1st Gen couples tend to have random spins (average spin is zero).
2nd/3rd Gen couples are the "super-spinners." Because they were formed in a previous merger where the spin was aligned with the disk, they carry that momentum forward. They are also heavier.
The Pattern: The heavier, multi-merger black holes (2g/3g) tend to pair up with smaller partners, creating that "small partner = high spin" pattern.

4. The "Guest List" (Initial Mass)

The researchers changed the "guest list" of black holes entering the disk.

Steep List (More small black holes): This creates more couples of similar size.
Flat List (More big black holes): This creates more "Giant + Tiny" couples.
Finding: A "steep" list (more small black holes) fits the observed data better. It suggests the universe prefers forming black holes that are closer in size to each other.

5. The "Direction" Problem (Prograde vs. Retrograde)

In physics, "prograde" means spinning with the flow of the disk, and "retrograde" means spinning against it.

Finding: For the simulation to match reality, almost all (over 90%) of the black hole couples must be spinning with the flow (prograde).
Analogy: If too many dancers are spinning the wrong way (retrograde), the pattern gets messy and the "anti-correlation" disappears. The data suggests the dance floor forces everyone to spin in the same direction.

6. The "Warm-Up" (Eccentricity)

Black holes don't always start in perfect circles; some have wobbly, elliptical orbits.

Finding: If the black holes start with very wobbly orbits, it takes a long time for the gas to smooth them out into circles so they can merge.
The Sweet Spot: The best results come from a "dynamically cooled" population—black holes that have already settled down a bit before entering the main dance. If they are too wobbly, they don't merge in time.

The Bottom Line

The paper concludes that the dense, short-lived, and crowded AGN disk is the most likely place to find the specific type of black hole mergers we are seeing in the universe.

If the disk is too thin or too long-lived: You don't get enough mergers, and the pattern disappears.
If the disk is dense and short-lived: You get a perfect match for the "small partner = high spin" rule, driven by heavy, multi-merger black holes that have been "trained" to spin with the flow of the galaxy.

Essentially, the universe's "black hole dance floor" needs to be a high-energy, crowded, short-duration event to produce the couples we are currently detecting.

1. Problem Statement

Gravitational wave (GW) observations by LIGO-Virgo-KAGRA (LVK) have revealed a population of binary black hole (BBH) mergers with individual masses in the "upper mass gap" ( $\sim 50-120 M_\odot$ ), which are difficult to explain via standard stellar evolution. While hierarchical mergers in dense environments (like globular clusters or nuclear star clusters) are a proposed solution, the Active Galactic Nucleus (AGN) channel is particularly compelling due to its efficiency in forming BBHs and retaining merger remnants.

A key observational feature is an apparent anti-correlation between the binary mass ratio ( $q = M_2/M_1$ ) and the effective spin ( $\chi_{\text{eff}}$ ): systems with lower mass ratios tend to have higher effective spins. Standard dynamical channels struggle to reproduce this specific correlation. This paper investigates whether the AGN channel can naturally produce this $(q, \chi_{\text{eff}})$ anti-correlation and identifies which physical parameters of AGN disks and Nuclear Star Clusters (NSCs) are required to match LVK observations.

2. Methodology

The authors utilize McFACTS (Monte Carlo For AGN Channel Testing and Simulation), a new open-source code designed to simulate BBH evolution in AGN environments. The study employs a Monte Carlo approach, running simulations for 100 distinct galaxies to build statistically significant samples.

Key Physical Processes Modeled:

BH Population: Black holes (BHs) are populated from the NSC into the AGN disk.
Gas Dynamics: Gas damping circularizes eccentric orbits, followed by radial migration driven by migration torques.
Binary Formation & Hardening: BBHs form when orbital separations reach the Hill radius. Gas and dynamical encounters cause binaries to harden (shrink) or soften.
Merger & Recoil: Merged BHs receive natal spins ( $\chi \sim 0.7$ ) and recoil kicks, potentially ejecting them from the disk or placing them on new eccentric orbits for further hierarchical merging.
Replenishment: Gas torques replenish the embedded BH population from the NSC.

Parameters Varied:
The study systematically varies the following parameters to test their impact on the $(q, \chi_{\text{eff}})$ distribution:

Disk Models: Sirko & Goodman (2003, SG; dense) vs. Thompson et al. (2005, TQM; less dense).
Initial Mass Function (IMF): Power-law index $\gamma$ (steep $\gamma=2$ vs. shallow/top-heavy $\gamma=1$ ).
Disk Lifetime ( $\tau_{\text{AGN}}$ ): Ranging from 0.25 Myr to 5 Myr.
Disk Size ( $R_{\text{out}}$ ): Varying the outer radius to alter the interaction cross-section with the NSC.
Initial BH Eccentricity ( $e_{\text{max}}$ ): From circular ( $e=0$ ) to thermal distributions ( $e \sim 0.7$ ).
Retrograde Fraction ( $f_{\text{ret}}$ ): The fraction of BBHs with angular momentum anti-aligned with the disk.
Initial Spin Magnitude ( $\sigma_\chi$ ): Width of the initial spin distribution.

3. Key Contributions

Development of McFACTS: Introduction of a flexible, reproducible code that improves upon previous models by incorporating detailed disk-NSC interactions, migration traps, and radiative feedback.
Parameter Space Mapping: A comprehensive exploration of how AGN disk properties (density, lifetime, size) and NSC conditions (eccentricity, IMF) dictate the observable GW signatures.
Quantification of Hierarchical Signatures: Detailed analysis of how "generation" (1g, 2g, $\ge$ 3g) affects the mass-spin distribution, specifically isolating the contribution of hierarchical mergers to the observed anti-correlation.

4. Key Results

A. Reproducing the Anti-Correlation

The study finds that dense, moderately short-lived AGN disks (SG model) with a steep IMF ( $\gamma=2$ ) and dynamically cooled NSCs (low initial eccentricity) best reproduce the observed $(q, \chi_{\text{eff}})$ anti-correlation.

Default SG Model: Produces a slope of $d\chi/dq \approx -0.28$ , matching LVK O3 data (including GW190814).
Mechanism: 1st generation (1g) mergers cluster around $\chi_{\text{eff}} \approx 0$ . Higher-generation (2g, $\ge$ 3g) mergers inherit large positive spins from previous mergers and tend to have lower mass ratios ( $q$ ), creating the anti-correlation.

B. Impact of Disk Density and Lifetime

SG vs. TQM: Dense SG disks facilitate rapid migration and high merger rates, producing significant hierarchical populations. Low-density TQM disks produce far fewer mergers and struggle to generate hierarchical events unless the initial BH population is extremely circular ( $e \approx 0$ ) and the disk lifetime is extended.
Lifetime: Longer disk lifetimes in SG models allow for more hierarchical mergers, steepening the anti-correlation slope (more high- $\chi_{\text{eff}}$ , low- $q$ events).

C. Influence of Initial Conditions

IMF: A shallower IMF ( $\gamma=1$ ) increases the fraction of high-mass BHs, leading to more hierarchical mergers and a steeper anti-correlation slope.
Eccentricity: High initial eccentricity suppresses mergers in low-density disks due to long circularization timescales. In dense disks, high eccentricity can lead to spin-down (negative $\chi_{\text{eff}}$ ) if the circularization period is long.
Retrograde Fraction ( $f_{\text{ret}}$ ): This is a critical parameter.
- $f_{\text{ret}} = 0$ (all prograde): Produces the strong anti-correlation.
- $f_{\text{ret}} = 0.5$ (50% retrograde): Symmetrizes the $\chi_{\text{eff}}$ distribution around zero, flattening the slope to $d\chi/dq \approx 0$ .
- Constraint: To match observations, the fraction of retrograde binaries must be $< 10\%$ (preferably $>90\%$ prograde).

D. Specific Event Interpretations

The model successfully places specific LVK events in the context of AGN formation:

GW190521: Consistent with a hierarchical merger (2g-2g or 3g) in an AGN disk.
GW231123: Overlaps with the high- $\chi_{\text{eff}}$ , high- $q$ region populated by $\ge$ 3g events, suggesting a complex hierarchical origin.
GW190814: The low mass ratio ( $q \sim 0.1$ ) is difficult for other channels but achievable in AGN disks with specific IMF choices.

5. Significance and Conclusions

Validation of AGN Channel: The paper demonstrates that the AGN channel is a viable and robust explanation for the observed $(q, \chi_{\text{eff}})$ anti-correlation, provided the environment is a dense AGN disk with dynamically cooled black hole populations.
Constraints on Physics: Future LVK data (O4/O5) can tightly constrain:
- The retrograde fraction of BBHs in AGN (likely $<10\%$ ).
- The initial spin distribution of BHs ( $\sigma_\chi \le 0.2$ ).
- The density and lifetime of AGN disks (favoring dense, short-lived episodes).
Hierarchical Mergers: The study confirms that hierarchical mergers are the primary driver of high-spin, low-mass-ratio events in the AGN channel, distinguishing them from isolated binary evolution or gas-free dynamical channels.

In summary, McFACTS II provides a rigorous framework for interpreting GW data, suggesting that the observed population of BBH mergers strongly favors formation in dense, short-lived AGN disks where black holes are dynamically cooled and predominantly prograde.