⚛️ phenomenology

Coupled Time-Dependent Proton Acceleration and Leptonic-Hadronic Radiation in Turbulent Supermassive Black Hole Coronae

This paper presents a time-dependent numerical framework that self-consistently couples proton acceleration with leptonic-hadronic radiation to successfully model multi-messenger signals from both steady sources like NGC 1068 and transient events such as tidal disruption events, revealing how cascade feedback can significantly delay electromagnetic and neutrino emissions.

Original authors: Chengchao Yuan, Damiano F. G. Fiorillo, Maria Petropoulou, Qinrui Liu

Published 2026-02-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Chengchao Yuan, Damiano F. G. Fiorillo, Maria Petropoulou, Qinrui 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 a supermassive black hole at the center of a galaxy not just as a cosmic vacuum cleaner, but as a chaotic, super-heated kitchen. In this kitchen, the "corona" is a swirling, turbulent cloud of hot gas and magnetic fields hovering just above the black hole. This paper presents a new, highly detailed recipe for simulating what happens in this kitchen, specifically focusing on how tiny particles (protons) get accelerated to incredible speeds and how they create a storm of light and invisible particles (neutrinos).

Here is a breakdown of the paper's main ideas using simple analogies:

1. The Problem: A Race Against Time

In this cosmic kitchen, three things happen at roughly the same speed:

Acceleration: Magnetic turbulence acts like a giant, chaotic pinball machine, bumping protons and speeding them up.
Cooling: As these protons speed up, they crash into photons (light particles), losing energy and creating new particles.
Cascades: These new particles crash into more things, creating a chain reaction (a cascade) that creates even more light and particles.

Previously, scientists struggled to model this because these processes happen so fast and influence each other so deeply. It's like trying to predict the weather while the wind, rain, and temperature are all changing every second based on each other.

2. The Solution: A New "Time-Traveling" Simulator

The authors built a new computer code (a numerical framework) that acts like a high-speed, time-traveling simulator. Instead of just guessing the final result, it watches the story unfold second by second.

The Engine: It uses a mathematical equation (the Fokker-Planck equation) to track how the protons move and speed up.
The Feedback Loop: Crucially, this simulator talks to another program (called AM3) that calculates the radiation. If the protons create a burst of light, that light immediately goes back and slows the protons down. The simulator updates the protons, then the light, then the protons again, over and over, in real-time.

3. Test Case A: The Steady Kitchen (NGC 1068)

The team first tested their simulator on a "steady" black hole called NGC 1068. This is a galaxy that has been spewing out high-energy neutrinos (ghostly particles that rarely interact with matter) for a long time.

The Result: The simulator successfully recreated the exact pattern of neutrinos detected by the IceCube telescope in Antarctica.
The Check: It also made sure the model didn't produce too much gamma-ray light, which would have contradicted what other telescopes see.
The Takeaway: The model proves that a turbulent "kitchen" near a black hole is a very likely place where these neutrinos are born.

4. Test Case B: The Explosive Kitchen (TDEs)

Next, they looked at a "transient" event called a Tidal Disruption Event (TDE). Imagine a star wandering too close to a black hole and getting ripped apart. This creates a temporary, violent flare-up. They used a specific event, AT 2019dsg, as their test subject.

The Surprise: In these weaker, temporary coronas, the "cascades" (the chain reactions of light) become very important. The light created by the protons doesn't just fly away; it bounces back and hits the protons, slowing them down significantly.
The Delay: Because of this feedback, the model predicts a strange delay. The black hole might start eating the star, but the resulting burst of light (in ultraviolet and X-rays) and neutrinos might not peak until 100 days later. It's like lighting a fuse on a firework, but the explosion happens long after the fuse is lit because the heat is building up slowly.

5. Why This Matters

The authors created a flexible tool that can be used for many different cosmic events, not just black holes.

Versatility: Whether the particles are being accelerated by turbulence, magnetic reconnection (like snapping rubber bands), or shockwaves, this tool can handle it.
Multi-Messenger Astronomy: It helps scientists connect the dots between different types of signals: light (optical, X-ray, gamma), particles (neutrinos), and gravity.
Future Proof: As new telescopes come online (like IceCube-Gen2), this tool will help astronomers interpret what they see, bridging the gap between the tiny physics of particles and the massive physics of black holes.

In summary: The paper introduces a powerful new way to simulate the chaotic dance between protons and light near black holes. It successfully explains the steady neutrino emissions from one galaxy and predicts a delayed, multi-wavelength explosion for another type of event, showing that the "echo" of light can significantly change how particles behave in space.

Technical Summary: Coupled Time-Dependent Proton Acceleration and Leptonic-Hadronic Radiation in Turbulent Supermassive Black Hole Coronae

Problem Statement
Turbulent coronae surrounding supermassive black holes (SMBHs) are potential sites for accelerating non-thermal particles to high energies, generating observable multi-messenger radiation (neutrinos and electromagnetic waves). However, modeling these processes is challenging due to the comparable timescales of particle acceleration, radiative cooling, and the development of electromagnetic (EM) cascades. A self-consistent treatment is required to incorporate the feedback of EM cascades on proton cooling and acceleration, particularly in compact coronal regions where photon densities are high. While steady-state models exist, there is a lack of systematic frameworks that couple time-dependent proton acceleration with radiative feedback across multiple spatial and temporal scales, especially for transient events like tidal disruption events (TDEs).

Methodology
The authors present a numerical framework that self-consistently couples time-dependent proton acceleration with leptonic-hadronic radiation modeling.

Proton Acceleration: The evolution of the proton distribution function $f(p, t)$ is governed by the time-dependent Fokker-Planck (FP) equation. This equation incorporates momentum diffusion (acceleration), escape, and cooling terms. The authors employ the Chang-Cooper method for spatial discretization and the Crank-Nicolson scheme for time discretization to ensure stability and positivity, even in cooling-dominant regimes.
Radiation Modeling: The FP solver is coupled with AM3, an open-source software for time-dependent lepto-hadronic radiation modeling. AM3 calculates the evolution of secondary particles (photons, neutrinos, $e^\pm$ , $\mu^\pm$ , $\pi^\pm$ , $\pi^0$ ) resulting from $p\gamma$ , $pp$, $\gamma\gamma$ , proton synchrotron, and Bethe-Heitler interactions.
Coupling Mechanism: To avoid double-counting cooling, the FP equation exclusively determines the evolution of the proton spectrum. At each time step, AM3 computes the cooling rates ( $t_{cool}^{-1} = t_{p\gamma}^{-1} + t_{pp}^{-1} + t_{BH}^{-1} + t_{p-sy}^{-1}$ ) based on the current proton distribution and photon fields. These rates are fed back into the FP solver to update the proton evolution.
Validation: The code is validated against steady-state solutions (reduced ordinary differential equations) and tested with specific injection functions and acceleration/cooling profiles to ensure convergence and independence from injection details in the high-energy regime.

Key Applications and Results
The framework is applied to two distinct astrophysical scenarios:

Steady Corona in NGC 1068:
- Setup: The model simulates a magnetized turbulent corona around the SMBH in the Seyfert galaxy NGC 1068, a known high-energy neutrino source.
- Results: The time-dependent simulation reproduces the IceCube neutrino spectrum (peaking at 1–100 TeV) while satisfying gamma-ray constraints imposed by Fermi-LAT. The model demonstrates that the system converges to a steady state within $\sim 300$ free-fall timescales ( $t_{fs}$ ).
- Findings: In this high-luminosity X-ray environment, proton cooling is dominated by interactions with coronal X-ray photons. The feedback from EM cascades on the proton distribution is found to be negligible, and the time-dependent solution aligns closely with previous steady-state results.
Transient Corona in TDE AT 2019dsg:
- Setup: A phenomenological model is developed for a transient corona powered by a TDE, characterized by a mass fallback rate $\propto t^{-5/3}$ . The model links the time-dependent accretion rate to corona properties (density, opacity).
- Results: The simulation predicts delayed multi-wavelength and neutrino emissions. Specifically, early-stage EM cascade feedback significantly impacts proton cooling in weaker coronae (where X-ray luminosity is low).
- Findings: The buildup of EM cascades, combined with proton accumulation during the super-Eddington phase, drives optical/UV, X-ray, and neutrino light curves to peak at a delay of approximately $O(100)$ days (specifically $\sim 140$ days in the model). This offers a potential explanation for time-delayed neutrino emission observed in TDEs.

Significance and Claims
The paper claims to provide the first systematic description coupling time-dependent proton acceleration (via Fokker-Planck) with radiative feedback in a fully time-dependent manner.

Versatility: The framework is presented as a flexible tool capable of modeling both steady sources (e.g., NGC 1068, other neutrino-emitting Seyferts) and transient phenomena (e.g., TDEs, gamma-ray bursts, AGN flares).
Mechanism Agnosticism: While the primary application focuses on turbulent acceleration, the authors note the framework can accommodate other mechanisms such as magnetic reconnection, shearing flows, and shock acceleration by modifying the diffusion and drift terms.
Multi-Messenger Insight: The work bridges plasma dynamics (acceleration microphysics) and radiation modeling, offering a tool to interpret high-energy astrophysical phenomena in the context of upcoming observatories like IceCube-Gen2 and KM3NeT.
Modest Limitations: The authors acknowledge that the current proton acceleration model relies on existing parameterizations (e.g., diffusion coefficients derived from PIC simulations) and does not yet self-consistently model the feedback of accelerated protons on the turbulence cascade itself, nor the detailed origin of coronal X-rays. They suggest these as directions for future refinement.