🔭 astrophysics

Can magnetic reconnection power neutrino emission from AGN coronae?

This paper proposes that magnetic reconnection in the transrelativistic, low- $\beta$ coronae of active galactic nuclei, specifically modeled for NGC 1068, can efficiently accelerate protons to PeV energies to produce the observed high-energy neutrino emission without requiring fine-tuned spectral parameters.

Original authors: Omar French, Gregory R. Werner, Mitchell C. Begelman

Published 2026-05-12

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Omar French, Gregory R. Werner, Mitchell C. Begelman

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, like the one in NGC 1068. Around this black hole is a super-hot, super-dense "corona" (like the sun's atmosphere, but much more extreme). This corona is a chaotic soup of magnetic fields and particles, constantly churning like a violent storm.

For years, scientists have been puzzled by a mystery: Where do high-energy neutrinos come from?

Neutrinos are tiny, ghost-like particles that can pass through almost anything. The IceCube observatory in Antarctica has detected a steady stream of these high-energy neutrinos coming from NGC 1068. But here's the catch: we don't see the high-energy gamma rays that usually accompany them. This suggests the neutrinos are being made in a very crowded, "opaque" place where gamma rays get trapped and eaten before they can escape. The authors of this paper propose that the black hole's corona is exactly that place.

The Big Question: How do you get a particle moving that fast?

To make the neutrinos IceCube sees, you need protons (the nuclei of hydrogen atoms) moving at nearly the speed of light, with energies in the "PeV" range (quadrillions of electron volts). It's like trying to accelerate a bowling ball to the speed of a bullet using only a gentle breeze.

The paper asks: Can magnetic reconnection do the job?

The Mechanism: The "Magnetic Rubber Band"

Think of the magnetic fields in the corona as tangled rubber bands.

The Tangle: The corona is turbulent, meaning these magnetic "rubber bands" are constantly twisting and snapping.
The Snap (Reconnection): Sometimes, these bands snap and reconnect. When they do, they release a massive burst of energy, creating a "current sheet" (a thin, intense layer of electricity).
The Acceleration: As protons drift through this chaotic environment, they occasionally get caught in these snapping magnetic sheets. It's like a surfer catching a wave. Every time a proton hits one of these "magnetic waves," it gets a little kick.
The Staircase: Because the corona is so turbulent, the proton doesn't just get one kick; it gets hit by thousands of these reconnecting sheets over and over again. This repeated battering accelerates the proton to incredible speeds.

The Energy Budget: Solving the Puzzle

The authors didn't just guess; they built a mathematical model to see if this idea holds water. They used three known facts about NGC 1068 as "clues":

The Neutrino Output: How many neutrinos IceCube sees.
The X-Ray Output: How much X-ray light the corona emits.
The Density: How thick the "fog" of particles is in the corona.

By plugging these clues into their equations, they found that the corona's magnetic field strength and particle density are locked together in a specific relationship. They call this a "one-parameter family," which basically means: If you know the magnetic field strength, you automatically know the density and size of the corona.

The Result: It Works!

When they ran the numbers, they found:

The Magnetization: The magnetic energy in the corona is strong enough to be "transrelativistic" (a fancy way of saying the magnetic fields are so strong they act like they have mass). This is the perfect regime for the magnetic "kicks" to work efficiently.
The Speed Limit: The protons get accelerated up to tens of PeV. However, they can't go faster because they start crashing into photons (light particles) and losing energy. This "speed limit" naturally matches the energy range of the neutrinos IceCube detects.
The Spectrum: The distribution of particle energies they predicted matches the shape of the neutrino signal IceCube sees. They didn't have to "tune" the model to make it fit; the physics of the magnetic reconnection naturally produced the right result.

Why No Gamma Rays?

You might wonder, "If protons are crashing and making neutrinos, why don't we see the gamma rays?"
The answer is the crowded room. The corona is so dense with matter and radiation that any gamma rays produced are immediately absorbed by other particles and reprocessed into lower-energy light. The neutrinos, being ghosts, slip right out. This explains why IceCube sees the neutrinos but telescopes like MAGIC see no high-energy gamma rays.

The Bottom Line

This paper argues that the chaotic, magnetic storm around the black hole in NGC 1068 acts as a natural particle accelerator. Through a process called magnetic reconnection, where magnetic field lines snap and release energy, protons are repeatedly kicked to extreme speeds. These super-fast protons then collide to create the neutrinos we detect on Earth. The model fits the observations perfectly without needing to invent new physics, suggesting that magnetic reconnection is likely the engine powering these cosmic neutrino factories.

Technical Summary: Can Magnetic Reconnection Power Neutrino Emission from AGN Coronae?

Problem Statement
The IceCube Neutrino Observatory has detected a persistent high-energy neutrino excess from the Seyfert galaxy NGC 1068, with energies ranging from $\sim1.5$ to $15$ TeV. The absence of accompanying TeV $\gamma$ -rays implies that the neutrino source is a compact, $\gamma$ -ray-opaque region, likely the supermassive black hole (SMBH) corona. While hadronic processes (proton-proton and proton-photon collisions) can produce these neutrinos, the specific mechanism capable of accelerating protons to the required energies ( $\sim30–300$ TeV) within the dense, magnetized corona remains unconstrained. Previous models often tuned the proton spectral slope to fit observations rather than deriving it from first principles. This paper investigates whether magnetic reconnection in small-scale current sheets within a transrelativistic, turbulent corona can supply the necessary nonthermal protons to power the observed neutrino emission.

Methodology
The authors model the AGN corona as a strongly turbulent, low- $\beta$ , collisionless hydrogen plasma characterized by four parameters: magnetic field strength ( $B$ ), proton density ( $n_p$ ), radiation energy density ( $u_{rad}$ ), and coronal size ( $r_{co}$ ).

Constraining the Environment: The model utilizes the observed IceCube-band neutrino luminosity ( $L_\nu$ ), the X-ray luminosity ( $L_X$ ), and an assumed Thomson optical depth ( $\tau_T$ ) as inputs. These constraints reduce the four-parameter coronal model to a one-parameter family defined by a scaling factor $\lambda$ (where $B = \lambda B_0$ ).
Acceleration Mechanism: The authors model proton acceleration via intermittent encounters with reconnecting current sheets embedded in a turbulent cascade. They adopt a tearing-mediated cascade picture where the magnetohydrodynamic (MHD) cascade transitions to kinetic scales near the proton skin depth. Protons are energized by the reconnection electric field within these sheets. The acceleration rate is calculated by averaging over a power-law distribution of current sheet lengths ( $dN/dL \sim L^{-7/2}$ ), accounting for the geometric probability of a proton encountering a sheet large enough to confine its gyroradius.
Cooling and Escape: The steady-state proton spectrum is determined by balancing injection, radiative cooling (synchrotron, inverse-Compton, $pp$, photomeson, and Bethe-Heitler processes), and advective escape into the black hole. The maximum proton energy is set by the competition between acceleration and photomeson cooling.
Spectral Synthesis: Motivated by Particle-in-Cell (PIC) simulations of strong turbulence at comparable magnetization, the authors adopt a nonthermal proton injection spectrum with a fixed index ( $s=3.5$ ) and a high-energy cutoff determined by the acceleration-cooling balance. They then compute the resulting neutrino spectral energy distribution (SED) without fitting the proton spectral slope to the data.

Key Results

Transrelativistic Regime: Across the viable parameter family for NGC 1068, the proton magnetization is found to be transrelativistic, with $\sigma_p \equiv B^2/(4\pi n_p m_p c^2) \sim 0.3$ (ranging from 0.18 to 0.70).
Acceleration Limits: In this regime, repeated encounters with reconnecting current sheets can energize suprathermal protons up to tens of PeV. However, photomeson cooling limits the maximum attainable energy, creating a cutoff consistent with the observed neutrino band.
Spectral Consistency: Using a proton injection index of $s=3.5$ (derived from PIC simulations of transrelativistic turbulence), the model predicts a TeV neutrino spectral shape that is broadly consistent with the IceCube observations of NGC 1068. This agreement is achieved without tuning the proton spectral slope to the data.
Parameter Constraints: The viable coronal parameters span a range where the magnetic field is $6–26$ kG, proton density is $1.1–5.1 \times 10^{10}$ cm $^{-3}$ , and the coronal size is $6–29$ gravitational radii ( $r_g$ ). The required magnetic cascade power is approximately $0.3–0.4$ times the X-ray luminosity.

Significance and Claims
The paper claims to demonstrate that magnetic reconnection in a turbulent, transrelativistic SMBH corona is a viable mechanism for powering the high-energy neutrino emission from NGC 1068. A primary contribution is the reduction of the coronal parameter space to a one-parameter family using observational constraints, thereby moving beyond models that arbitrarily select spectral indices. The study shows that the transrelativistic magnetization naturally selected by the energy budget places the system in a regime where PIC simulations predict steep proton spectra ( $s \gtrsim 3$ ), which are required to match the observed neutrino SED.

The authors note that their model assumes a single-zone corona with uniform parameters and an advective escape mechanism. They acknowledge that realistic radial gradients, black hole spin effects (which could expand the viable parameter range), and more complex coupling between reconnection injection and stochastic turbulence remain open areas for refinement. However, the work establishes that reconnection-powered models can plausibly produce neutrino fluxes consistent with IceCube observations without fine-tuning the particle acceleration spectrum.