🔭 astrophysics

McFACTS IV: Electromagnetic Counterparts to AGN Disk Embedded Binary Black Hole Mergers

This paper introduces an enhanced version of the \texttt{McFACTS} simulation code to predict the bolometric electromagnetic luminosities of jets and shocks from binary black hole mergers in AGN disks, demonstrating that such events in dense disks with specific mass and spin characteristics are highly likely to produce observable electromagnetic counterparts.

Original authors: Emily McPike, Rosalba Perna, K. E. Saavik Ford, Barry McKernan, Vera Delfavero, Miranda McCarthy, Kaila Nathaniel, Jake Postiglione, Nicolas Posner, Varun Pritmani, Shawn Ray, Richard O'Shaughnessy

Published 2026-02-05

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Emily McPike, Rosalba Perna, K. E. Saavik Ford, Barry McKernan, Vera Delfavero, Miranda McCarthy, Kaila Nathaniel, Jake Postiglione, Nicolas Posner, Varun Pritmani, Shawn Ray, 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

The Big Picture: Black Holes in a Cosmic Soup

Imagine the center of a galaxy as a giant, swirling whirlpool of gas and dust, known as an Active Galactic Nucleus (AGN). In the middle of this whirlpool sits a supermassive black hole, like a massive anchor.

Floating in this "cosmic soup" are smaller black holes. Usually, these small black holes just drift around. But in this specific environment, the gas acts like a cosmic traffic cop. It slows the black holes down, forcing them to migrate toward specific "parking spots" (called migration traps) where they get crowded together.

When two of these black holes get close enough, they crash into each other and merge. This is a Binary Black Hole (BBH) merger. We already know these crashes happen because our gravitational wave detectors (LIGO/Virgo) hear the "thud" of the collision.

The Big Question: Can we see these crashes, too?

Usually, black hole mergers are invisible. They only send out ripples in space (gravitational waves). But because these mergers happen inside a thick, gas-filled disk, the crash might kick up a massive splash of light (an electromagnetic counterpart) that telescopes could see. This paper tries to predict when and how bright that splash will be.

The Tool: A Cosmic Simulator

The authors used a computer code called McFACTS (Monte Carlo For AGN Channel Testing and Simulation). Think of this code as a massive, high-speed video game simulator.

It creates 100 different galaxies.
It populates them with black holes of different sizes and spins.
It runs the simulation forward in time to see which black holes crash, how many times they merge (some merge, then merge again with a new partner), and what happens to the gas around them.

The Two Types of "Splashes"

When a black hole merger happens in this gas disk, the paper looks at two ways it might create light:

The Shockwave (The "Bow Wave"):
When the two black holes merge, the new, bigger black hole gets a sudden kick (like a cannonball firing). It shoots through the gas disk. As it plows through the gas, it creates a shockwave, similar to the bow wave of a speedboat cutting through water.
- The Paper's Finding: The authors calculated how bright this shockwave would be. They found that for the types of disks they modeled, this shockwave is usually too dim to be seen. It's like a speedboat making a small ripple in a massive ocean—it gets lost in the background noise.
The Jet (The "Firehose"):
If the new black hole is spinning very fast and is surrounded by enough gas, it can launch a powerful jet of energy straight up out of the disk, like a firehose spraying water.
- The Paper's Finding: This is the promising one. If the black hole is heavy enough and spinning fast enough, this jet can be incredibly bright—brighter than the entire galaxy's center!

The "Breakout" Challenge

Here is the tricky part: The galaxy's gas disk is thick and opaque (like a very dense fog). Even if a powerful jet is launched, it might get stuck inside the fog.

The Analogy: Imagine trying to blow a bubble through a thick blanket. If you don't blow hard enough or fast enough, the bubble pops inside the blanket, and you see nothing.
The Paper's Rule: For us to see the light, the jet has to punch a hole through the disk and escape before the engine (the black hole) runs out of fuel. The authors calculated that only the heaviest, fastest-spinning black holes (usually those that have merged multiple times) have enough power to break out of the "blanket" and be seen by our telescopes.

What Makes a Visible Flash?

The authors tested many different scenarios to see what conditions create the brightest flashes. They found three main ingredients are needed:

Heavy Black Holes: You need "top-heavy" black holes (ones that started very massive). It's easier to build a massive, powerful engine if you start with big parts.
Long-Lived Disks: The gas disk needs to stick around long enough (about a million years) for the black holes to have time to merge, merge again, and build up that massive size. If the disk disappears too fast, the black holes don't get a chance to grow.
Migration Traps: The gas needs to have specific "traffic jams" where black holes get stuck together. This helps them merge efficiently.

The "Who to Watch" List

Based on their simulations, the authors give a specific recommendation for astronomers:

Don't waste time looking for flashes from small, first-time mergers. They are likely too dim or get stuck in the gas.
Focus your telescopes on the massive mergers. Specifically, if a gravitational wave detector hears a crash involving black holes with a combined "chirp mass" (a measure of their weight) of 44 times the mass of our Sun or more, there is a very high chance (up to 100% in their best-case models) that a bright flash will follow.

The Catch: Geometry

Even if a bright flash happens, we might not see it.

The Analogy: Imagine a lighthouse. If the beam is pointing at you, you see it. If it's pointing away, or if a thick wall (the galaxy disk) is in the way, you see nothing.
The paper notes that because of the angle of the galaxy and the direction the black hole is kicked, we might only see about 1 out of every 4 of these potential flashes.

Summary

This paper is a "recipe book" for finding light from black hole crashes. It tells us that while small crashes in gas disks are likely invisible, the giant, multi-generational crashes in dense, long-lasting gas disks are the best candidates for producing a visible flash.

If we want to find these flashes, we should point our telescopes at the locations of the heaviest black hole mergers and watch them for a few months after the crash, hoping to catch the "firehose" of light breaking out of the cosmic fog.

Technical Summary: McFACTS IV: Electromagnetic Counterparts to AGN Disk Embedded Binary Black Hole Mergers

Problem Statement
Active Galactic Nuclei (AGN) accretion disks are a leading theoretical channel for the formation and merger of binary black holes (BBHs), a population increasingly detected by the LIGO-Virgo-KAGRA (LVK) network. While BBH mergers in AGN disks are unique among formation channels in their potential to produce electromagnetic (EM) counterparts via interactions with surrounding gas, the lack of predictive statistical models currently hinders the identification of these counterparts in archival data and real-time time-domain surveys (e.g., ZTF, LSST). Furthermore, the fraction of LVK events originating from the AGN channel ( $f_{BBH,AGN}$ ) remains uncertain, and distinguishing true EM counterparts from intrinsic AGN variability is challenging due to poor GW localization and the bright, variable nature of AGN.

Methodology
This work utilizes an updated version of the Monte Carlo For AGN Channel Testing and Simulation (McFACTS v.0.4.0) code, a 1D rapid population synthesis tool for dynamic compact binary formation. The study introduces a new capability to self-consistently generate distributions of bolometric EM luminosities (jets and shocks) corresponding to predicted LVK-detectable GW events.

The methodology involves:

Population Synthesis: Simulating BBH mergers within a Sirko-Goodman AGN disk model surrounding a $10^8 M_\odot$ supermassive black hole (SMBH). The code tracks hierarchical mergers, spin evolution, and migration.
EM Emission Modeling:
- Shock Luminosity ( $L_{shock}$ ): Calculated based on the energy dissipated when the GW recoil kick displaces the remnant's Hill sphere gas into the surrounding disk (McKernan et al. 2019).
- Jet Luminosity ( $L_{jet}$ ): Modeled assuming the remnant enters a Magnetically Arrested Disk (MAD) state, launching a jet via the Blandford-Znajek mechanism. The power is parameterized by the Bondi-Hoyle-Lyttleton accretion rate, black hole spin ( $\chi$ ), and jet efficiency.
Observability Criteria: A critical component is the calculation of jet breakout time ( $t_{bo}$ ). A jet is considered potentially observable only if it breaks out of the optically thick disk before the engine turns off ( $t_{bo} < T_{jet}$ ). If $t_{bo} > T_{jet}$ , the emission is "muffled" and emerges on a photon diffusion timescale, rendering it indistinguishable from AGN variability.
Parameter Space Exploration: The study varies key physical parameters to assess their impact on EM counterpart rates and luminosities, including:
- AGN disk lifetime ( $\tau_{AGN}$ ).
- Disk radial extent ( $R_{outer}$ ).
- Disk turbulence strength (stochastic migration).
- Black Hole Initial Mass Function (IMF) (varying power-law indices and peak masses).
- Initial spin distribution width.

Key Results

Dominance of Jets over Shocks: In the fiducial Sirko-Goodman disk model, shock luminosities are generally sub-dominant to the host AGN luminosity ( $L_{AGN} \sim 10^{45}$ erg s $^{-1}$ ) and are unlikely to be distinguishable. Conversely, accretion-powered jets can significantly outshine the AGN.
Role of Hierarchical Mergers: Dense regions in the disk, specifically migration traps and swamps (e.g., at $R \sim 700 R_g$ ), efficiently drive hierarchical mergers. These processes produce high-mass, high-spin remnants (generations 2 and higher) capable of powering luminous jets.
Chirp Mass Dependence: The likelihood of producing an observable EM counterpart is strongly correlated with the chirp mass ( $\mathcal{M}$ ). For the fiducial model, the probability of an observable flare rises with mass, reaching $\sim 100\%$ for mergers with $\mathcal{M} \in [44, 88] M_\odot$ (before accounting for geometric obscuration).
Parameter Sensitivity:
- Disk Lifetime: Longer-lived disks ( $\tau_{AGN} \ge 1$ Myr) allow for more extensive hierarchical merging, increasing the population of high-spin, high-mass remnants and thus the rate of bright jets.
- IMF: Top-heavy IMFs (e.g., mass segregation or peak-only distributions around $34 M_\odot$ ) significantly enhance the efficiency of hierarchical growth, leading to more massive remnants and brighter jets compared to standard power-law IMFs.
- Turbulence: Increased disk turbulence smears out migration traps, slightly reducing the efficiency of hierarchical mergers and the resulting jet luminosity variance.
Geometric Constraints: Even when a jet breaks out, geometric factors limit observability. Assuming equal numbers of Type I and Type II AGN, only $\sim 1/4$ of breakout events are geometrically detectable (face-on Type I with breakout toward the observer).

Significance and Claims
The paper claims to present the first population synthesis study of EM counterparts specifically for the AGN channel. Its primary contribution is the development of a framework that links GW observables (chirp mass, spin) to predicted EM properties (luminosity, breakout time) within a physically motivated disk environment.

The authors assert that:

Target Selection: EM follow-up campaigns should prioritize GW events with chirp masses $\mathcal{M} \gtrsim 44 M_\odot$ , as these have the highest statistical probability of producing observable EM counterparts.
Timescales: Observable emission may emerge as late as $\sim 4$ months post-merger, suggesting that monitoring campaigns should extend over timescales of weeks to months.
Model Constraints: The detection or non-detection of counterparts in future surveys (e.g., LSST, Einstein Telescope) will provide critical constraints on AGN disk properties (density, scale height, lifetime) and the nature of the BBH population (IMF, spin alignment).
Multi-messenger Potential: Confident associations between GW events and EM flares in AGN disks offer a pathway to independently constrain the Hubble constant ( $H_0$ ).

The study remains modest regarding specific detection guarantees, acknowledging that uncertainties in AGN variability, jet launching physics, and geometric viewing angles currently limit the ability to definitively associate specific GW events with EM counterparts. The work serves as a guide for future observational strategies and theoretical refinements.