🔭 astrophysics

AGN Disks as Supernova Mufflers I: 3D Local Hydrodynamic Models

This paper employs 3D hydrodynamic simulations and analytic criteria to demonstrate that supernova shocks in active galactic nucleus accretion disks are muffled at radii dependent on the local disk scale height and black hole mass, with specific muffling thresholds identified for different disk models.

Original authors: Harrison E. Cook, Wladimir Lyra, Mordecai-Mark Mac Low, K. E. Saavik Ford, Barry McKernan

Published 2026-04-01

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Harrison E. Cook, Wladimir Lyra, Mordecai-Mark Mac Low, K. E. Saavik Ford, Barry McKernan

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 a supermassive black hole at the center of a galaxy, surrounded by a swirling, dense disk of gas and dust. This is an Active Galactic Nucleus (AGN). Think of this disk not just as a flat pancake, but as a thick, turbulent atmosphere, like a stormy ocean made of gas.

Now, imagine a massive star living inside this gas ocean. Eventually, that star explodes as a Supernova—a cosmic bomb releasing energy equivalent to a billion trillion nuclear detonations.

The big question this paper asks is: What happens when that bomb goes off underwater?

Does the explosion punch a hole straight through the gas ocean and blast into space? Or does the thick gas act like a heavy blanket, smothering the explosion and trapping all its energy inside the disk?

The authors call this phenomenon "muffling."

The Two Main Characters: The "Puffy" Disk and the "Thin" Disk

To figure this out, the scientists looked at two different theoretical models of what these gas disks look like:

The "Sirko & Goodman" (SG) Model: Think of this as a puffy, fluffy cloud. The gas gets thicker and puffier the closer you get to the black hole, but it gets very thin and spread out as you move away.
The "Thompson, Quataert & Murray" (TQM) Model: Think of this as a dense, flat sheet that stays thick and uniform for a long time before thinning out.

The "Muffling Radius": The Safety Zone

The researchers discovered a specific boundary line called the "Muffling Radius."

Inside the Radius (The Inner Disk): The gas is thin enough (or the explosion is close enough to the center) that the supernova has enough "oomph" to punch through the gas layers. It breaks out, sending a shockwave and light into the universe.
Outside the Radius (The Outer Disk): The gas is too thick, or the disk is too "puffy." The explosion hits the gas, slows down, and gets trapped. The energy doesn't escape as a bright flash; instead, it heats up the gas, creating turbulence and pressure, but no one outside sees the explosion.

The Analogy: Imagine throwing a stone into a shallow puddle versus a deep, thick swamp.

In the shallow puddle (inner disk), the stone splashes water high into the air (breakout).
In the thick swamp (outer disk), the stone sinks, creating a small bubble and ripples, but the water absorbs the impact (muffling).

Key Findings in Plain English

1. The Size of the Black Hole Matters (For Puffy Disks)
In the "puffy" disk model, the size of the black hole changes the rules.

If the black hole is small (1 million times the mass of our Sun), the "muffling zone" is huge. Even far out in the disk, the gas is thick enough to smother explosions.
If the black hole is giant (1 billion times the mass of our Sun), the disk is stretched out and thinner. The "muffling zone" shrinks down to a tiny area right next to the black hole. Most of the disk is thin enough for explosions to escape.

2. The "Flat" Disk is Different
In the "flat" disk model, the muffling zone is always in the same spot, no matter how big the black hole is. It's like a universal "danger zone" where explosions always get smothered.

3. Height is a Game-Changer
If a star explodes right in the middle of the gas disk (the "midplane"), it has to push through the thickest part of the gas to get out. It's like trying to swim out of a pool while holding your breath at the bottom.
However, if the star explodes slightly above or below the center line, it has a shorter path to the surface. It's like jumping out of the water from the side of the pool. The explosion escapes much easier and faster.

4. The "One-Way" Street
When an explosion happens off-center, it behaves strangely. The shockwave shooting up and away from the disk flies out quickly. But the shockwave shooting down has to fight through the dense gas to get to the other side. By the time it reaches the other side, it might have lost all its speed and become a slow, gentle flow rather than a violent shock. This means an observer on one side of the galaxy might see a bright flash, while an observer on the other side sees nothing.

Why Should We Care?

1. Finding the "Ghost" Explosions
Astronomers are looking for flashes of light that signal black holes merging (which creates gravitational waves). Sometimes, they see a flash that doesn't look like a normal supernova. This paper suggests those might be "muffled" supernovas trapped in AGN disks. Knowing this helps astronomers avoid false alarms.

2. Heating the Galaxy
Even if the explosion is "muffled" and doesn't send light out, it's not useless. The energy gets trapped in the gas, heating it up and keeping the disk from collapsing under its own gravity. It's like a radiator keeping a room warm. These trapped explosions might be the secret engine keeping these galactic disks stable.

3. Metal Delivery
Supernovas create heavy elements (like gold and iron). If the explosion is trapped, those heavy elements stay in the disk, enriching the gas. If it escapes, they get scattered into the galaxy. This helps explain why the centers of galaxies are so rich in heavy metals.

The Bottom Line

This paper uses powerful computer simulations to show that where a star explodes in a galactic disk determines its fate.

Inner/Thin regions: The explosion escapes, lighting up the universe.
Outer/Thick regions: The explosion is smothered, heating the gas and staying hidden.

It turns out that the "noise" of these hidden explosions might be just as important as the "signal" of the visible ones in understanding how galaxies grow and evolve.

1. Problem Statement

Active Galactic Nuclei (AGN) accretion disks are dense, stratified environments where stars and stellar-mass black holes can become embedded. While binary black hole mergers in these disks are a primary target for gravitational wave detection (e.g., GW190521), other energetic events, such as core-collapse Supernovae (SNe), may occur within the disk.

The Core Question: Do SNe occurring within an AGN disk break free and produce observable electromagnetic counterparts, or are they "muffled" (trapped), depositing their kinetic energy into the disk?
The Challenge: Unlike the uniform Interstellar Medium (ISM), AGN disks have non-uniform density, variable thickness (scale height), and strong shear. These factors significantly alter shock morphology and breakout conditions compared to standard ISM SNe.
Motivation: Distinguishing between AGN flares caused by binary mergers versus muffled SNe is critical for interpreting transient surveys (like LSST) and constraining the dynamics of nuclear star clusters.

2. Methodology

The authors employ a two-pronged approach combining analytic scaling arguments with high-resolution 3D hydrodynamic simulations.

A. Analytic Breakout Criterion

The authors derive a dimensionless parameter $C$ to predict shock breakout:
$C \equiv \frac{\epsilon_{SN}}{P_D} = \frac{3 m_p}{4 \pi k} \frac{E_{SN}}{\rho T H^3}$
Where:

$E_{SN}$ is the SN kinetic energy (assumed $10^{51}$ erg).
$\rho, T$ are the local gas density and temperature.
$H$ is the local vertical scale height of the disk.
$C > 1$ : The SN pressure exceeds the disk pressure; Breakout occurs.
$C < 1$ : The SN is Muffled; energy is deposited as heat into the disk.

They apply this criterion to two distinct AGN disk models generated by the pAGN code:

Sirko & Goodman (SG) Model: A standard thin-disk model.
Thompson, Quataert & Murray (TQM) Model: A model often used for radiation-pressure dominated disks.

B. 3D Hydrodynamic Simulations

To verify the analytic predictions, the authors perform 3D local shearing box simulations using the Athena code.

Setup: The simulations model a stratified Gaussian disk in hydrostatic equilibrium.
Initial Conditions: A symmetric 3D Gaussian perturbation is injected to represent the SN, adding mass ( $10 M_\odot$ ) and energy ( $10^{51}$ erg) to the grid.
Scenarios:
- Midplane Detonations: $z_0 = 0$ .
- Offset Detonations: $z_0 = 1H$ (simulating inclined orbits).
Parameters: Simulations are run for a Supermassive Black Hole (SMBH) mass of $M_{BH} = 10^8 M_\odot$ at radii $R = 10^3 R_s$ (expected breakout) and $R = 10^5 R_s$ (expected muffling).
Boundaries: Shear-periodic in radial/azimuthal directions; outflow in vertical directions.

3. Key Results

A. Analytic Predictions (Breakout Parameter $C$ )

SG Disks: The muffling radius ( $R_m$ $R_{m}$ , where $C=1$ $C = 1$ ) is inversely proportional to the SMBH mass.
- For $M_{BH} = 10^6 M_\odot$ , $R_m \approx 5 \times 10^6 R_s$ .
- For $M_{BH} = 10^9 M_\odot$ , $R_m \approx 10^2 R_s$ .
- Reasoning: Higher mass SMBHs have "puffier" disks (larger $H$ ) at a given radius, increasing the volume the shock must traverse ( $H^3$ dependence), making muffling more likely.
TQM Disks: The muffling radius is independent of SMBH mass, fixed at $R_m \approx 10^6 R_s$ $R_{m} \approx 1 0^{6} R_{s}$ .
- The TQM model features a distinct temperature drop and scale height trough around $R \sim 10^4 R_s$ , creating a localized peak in breakout probability before muffling sets in at larger radii.

B. Hydrodynamic Simulation Findings

Breakout Confirmation: At $R = 10^3 R_s$ ( $C \approx 2600$ ), the shock expands vertically, breaks through the disk surface within $\sim 30$ days, and evolves into a hemispherical plume.
Muffling Confirmation: At $R = 10^5 R_s$ ( $C \approx 0.04$ ), the shock fails to break out. Over 1200 years, it remains confined to $z \sim 1H$ . The disk shear elongates the cavity, eventually closing it. The kinetic energy is thermalized into the disk.
Offset Detonations ( $z_0 = 1H$ ):
- Breakout occurs faster (half the time) for the upward-propagating shock due to lower ambient density.
- The downward-propagating shock must climb the density gradient; it often becomes subsonic before reaching the midplane, failing to punch through the opposite side.
- This creates an asymmetric breakout, potentially leading to different observational signatures for observers on opposite sides of the disk.

4. Key Contributions

Predictive Analytic Criterion: Established a scaling relation ( $C$ ) that identifies the "muffling radius" based on local disk density, temperature, and scale height.
3D Verification: Provided the first 3D hydrodynamic confirmation that SN shocks in AGN disks can be completely muffled, depositing energy rather than escaping, validating the analytic $C < 1$ condition.
Model Differentiation: Demonstrated that the SG and TQM disk models yield fundamentally different breakout landscapes (mass-dependent vs. mass-independent muffling radii).
Inclination Effects: Quantified how inclined orbits ( $z > 0$ ) shift the muffling boundary outward and create asymmetric shock propagation.

5. Significance and Implications

Observational Constraints: The study suggests that many AGN flares previously attributed to binary mergers could be muffled SNe, or conversely, that true merger counterparts might be obscured if they occur in the outer, dense regions of the disk.
Disk Physics: Muffled SNe act as a localized heating source, potentially providing the pressure support necessary to prevent self-gravitating disks from collapsing.
Chemical Enrichment: Muffled SNe inject metals directly into the disk midplane, while breakout SNe can be swept up by AGN winds, spreading nucleosynthetic products to the outer disk.
Rate Estimates: Assuming a disk lifetime of $10^6$ years and a standard SN rate, the authors estimate that only a tiny fraction ( $\sim 10^{-6}$ ) of SNe in SG disks will actually break free and be observable as distinct transients. This implies that non-detection of AGN-SNe could indicate large, dense, "puffed-up" disks rather than a lack of stellar activity.

Conclusion: The paper establishes that the "muffling" of supernovae is a dominant physical process in AGN disks, heavily dependent on the SMBH mass and the specific disk model. This has profound implications for interpreting transient surveys and understanding the feedback loops between stars and AGN accretion disks.