arxiv: 2601.01533 · v2 · submitted 2026-01-04 · 🌌 astro-ph.HE

Recognition: no theorem link

Diffuse neutrino flux from relativistic reconnection in AGN coronae

D. Karavola , M. Petropoulou , D. F. G. Fiorillo , A. Georgakakis , L. Comisso , L. Sironi

Authors on Pith no claims yet

Pith reviewed 2026-05-16 17:58 UTC · model grok-4.3

classification 🌌 astro-ph.HE

keywords diffuse neutrino fluxAGN coronaemagnetic reconnectionIceCubeproton accelerationpion coolingnon-jetted AGN

0 comments

The pith

Magnetic reconnection in AGN coronae explains IceCube's diffuse neutrino flux up to 1 PeV if 10 percent of coronae reach high magnetization.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates whether neutrinos produced by protons accelerated through magnetic reconnection in the coronae of non-jetted active galactic nuclei can account for the diffuse flux detected by IceCube. The authors generate spectral templates for neutrino emission across ranges of proton magnetization, X-ray luminosity, and black hole mass, including the effects of strong synchrotron cooling on pions and muons. They combine these templates with a mock catalog of AGN that matches observed X-ray and infrared samples from redshift 0 to 4. A reader would care because the work identifies a concrete parameter window in which coronal reconnection alone reproduces the measured flux below PeV energies without invoking jets.

Core claim

The central claim is that coronal emission satisfactorily explains the most recent IceCube measurements of the diffuse neutrino flux up to energies of approximately 1 PeV, provided that about 10 percent of the AGN coronae have proton plasma magnetization sigma_p approximately 10^5 while the rest follow a distribution of lower magnetizations. Synchrotron cooling of pions and muons in the strong coronal magnetic fields suppresses emission at higher energies, so an additional population such as jetted AGN is needed above that range.

What carries the argument

A library of neutrino spectral templates generated over a grid of proton plasma magnetization sigma_p, X-ray coronal luminosity, and black hole mass, which is then integrated against a mock AGN catalog consistent with observed samples to yield the diffuse flux.

If this is right

Coronal reconnection matches the observed diffuse flux up to 1 PeV under the stated 10 percent high-magnetization condition.
Pion and muon synchrotron cooling naturally cuts off the coronal contribution above 1 PeV.
An additional neutrino population, with jetted AGN as strong candidates, is required to explain any flux at higher energies.
The model relies on efficient proton acceleration via reconnection in the dense X-ray photon field of the corona.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the required 10 percent high-magnetization fraction is real, it constrains the typical magnetic conditions inside AGN coronae across the population.
Multi-messenger searches could test the scenario by looking for spatial correlations between high-energy neutrinos and AGN with bright, compact X-ray coronae.
The cooling-induced cutoff supplies a spectral feature that future detectors could use to separate coronal neutrinos from contributions of other source classes.

Load-bearing premise

That a fixed 10 percent fraction of AGN coronae can be assigned sigma_p around 10^5 while the remainder follow a lower-magnetization distribution, without independent observational evidence for that split or for the reconnection efficiency.

What would settle it

A direct measurement of the fraction of AGN coronae with magnetization near 10^5 that differs substantially from 10 percent, or an observed neutrino spectrum above 1 PeV that lacks the sharp suppression predicted by pion and muon cooling.

read the original abstract

IceCube observations point to Active Galactic Nuclei (AGN) as promising contributors to the observed astrophysical neutrino flux. Close to the central black hole, protons can be accelerated through magnetic reconnection to very high energies and subsequently interact with abundant X-ray photons in the source, leading to neutrino production. We investigate whether the diffuse neutrino flux observed by IceCube can originate, via proton acceleration, in reconnection-powered coronae of non-jetted AGN. We create a library of neutrino spectral templates, over a large grid of values for the three key model parameters: the proton plasma magnetization of the corona $\sigma_{\rm p}$, the X-ray coronal luminosity, and the black hole mass. Synchrotron cooling of pions and muons plays a significant role due to the large coronal magnetic fields. We couple the single-source model with a mock AGN catalog, consistent with the observed X-ray and mid-infrared AGN samples at redshifts $z=0-4$, to infer the diffuse neutrino flux. Coronal emission satisfactorily explains the most recent IceCube measurements of the diffuse neutrino flux up to energies of $\sim 1$~PeV, provided that $\sim$10\% of the AGN coronae have $\sigma_{\rm p} \sim 10^5$, while the rest are distributed over a range of lower magnetizations. Coronal emission is suppressed at higher energies by pion and muon cooling, so that another population is required, with jetted AGN being strong candidates.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper builds useful new neutrino spectral templates for AGN coronae that include pion and muon cooling, but the match to IceCube data up to 1 PeV rests on tuning 10% of sources to extreme magnetization without independent justification.

read the letter

This paper shows that relativistic reconnection in non-jetted AGN coronae can reproduce the IceCube diffuse neutrino flux up to about 1 PeV, but only after assigning roughly 10% of the population a proton magnetization of order 10^5 while the rest sit at lower values. That split is chosen to set the overall normalization rather than derived from separate data or simulations. The result is therefore conditional on that choice. The work is worth reading for the templates themselves, which are new. They compute neutrino spectra across a grid in sigma_p, coronal X-ray luminosity, and black hole mass, and they include synchrotron cooling of pions and muons. Those losses matter because the coronal magnetic fields are strong, and they naturally suppress the flux above a few hundred TeV. The authors then fold the templates into a mock catalog built to match observed X-ray and mid-IR AGN samples from z=0 to 4. That step is concrete and improves on earlier estimates that stayed at the single-source level. The high-energy cutoff from cooling is a clear takeaway and correctly points to jetted sources for anything above PeV. The main limitation is the 10% high-magnetization subpopulation. Neutrino luminosity rises steeply with sigma_p, so the fraction directly controls whether the total flux lands on the IceCube points. The paper supplies no separate constraint from coronal variability, spin distributions, or PIC runs to motivate exactly that 10% level. The mock catalog is also constructed to be consistent with the same X-ray and mid-IR samples that set the normalization, which introduces moderate circularity. These issues are real but not fatal; they are the sort of thing referees can press on. The templates and the catalog approach are solid enough that the paper belongs in the neutrino astrophysics literature. I would bring it to a reading group for the modeling details and would cite the spectral library in future work on coronae. It deserves peer review because the calculation is explicit and the cooling treatment is an advance, even if the parameter tuning needs more scrutiny.

Referee Report

2 major / 2 minor

Summary. The paper models the diffuse neutrino flux from proton acceleration via relativistic magnetic reconnection in the coronae of non-jetted AGN. It generates a library of neutrino spectral templates over a grid in proton magnetization σ_p, X-ray coronal luminosity, and black-hole mass, incorporating synchrotron cooling of pions and muons. These templates are convolved with a mock AGN catalog consistent with X-ray and mid-IR samples at z=0–4. The central claim is that the IceCube diffuse flux up to ~1 PeV is satisfactorily reproduced provided ~10% of coronae have σ_p ~10^5 while the remainder follow a lower-magnetization distribution; higher energies are suppressed by cooling and require a separate (e.g., jetted) population.

Significance. If the result survives scrutiny of the magnetization distribution, the work supplies a concrete, reconnection-based mechanism that can account for a substantial fraction of the observed astrophysical neutrino background from non-jetted sources. The inclusion of detailed pion/muon cooling and the population synthesis via a mock catalog constitute genuine strengths that allow falsifiable spectral predictions and multi-messenger tests.

major comments (2)

[Abstract and §4] Abstract and §4 (results): the requirement that ~10% of AGN coronae have σ_p ~10^5 is introduced to normalize the total neutrino flux to IceCube data. Because neutrino luminosity scales steeply with σ_p (via E_p,max and pγ optical depth), this fraction functions as a free normalization parameter whose value is selected post-hoc rather than derived from independent constraints such as coronal X-ray variability, black-hole spin statistics, or PIC reconnection simulations.
[§3] §3 (mock catalog): the catalog is constructed to be consistent with the same X-ray and mid-IR AGN samples that set the luminosity function and thus the overall normalization. This introduces circularity between the input population statistics and the output diffuse flux, weakening the claim that the model independently explains the IceCube measurements.

minor comments (2)

[Template construction] Notation: σ_p is defined in the abstract but its precise relation to the magnetic field and proton density should be restated explicitly in the template-construction section for clarity.
[Figures] Figures: the comparison plots of model flux versus IceCube data would be strengthened by showing the envelope arising from variations in the low-magnetization subpopulation rather than a single central curve.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the model's strengths. We address each major comment below, clarifying the role of model parameters and the use of observational inputs.

read point-by-point responses

Referee: [Abstract and §4] Abstract and §4 (results): the requirement that ~10% of AGN coronae have σ_p ~10^5 is introduced to normalize the total neutrino flux to IceCube data. Because neutrino luminosity scales steeply with σ_p (via E_p,max and pγ optical depth), this fraction functions as a free normalization parameter whose value is selected post-hoc rather than derived from independent constraints such as coronal X-ray variability, black-hole spin statistics, or PIC reconnection simulations.

Authors: We agree that the ~10% fraction is selected to reproduce the observed flux normalization, as the neutrino output depends sensitively on σ_p. This is not hidden in the manuscript; the abstract and §4 explicitly present it as a requirement for the model to match IceCube data up to ~1 PeV. The spectral shape and high-energy cutoff, however, arise directly from the inclusion of pion/muon synchrotron cooling and are independent of the precise fraction. We will revise the text to discuss possible physical motivations for a high-σ_p tail, referencing recent PIC simulations of relativistic reconnection and X-ray variability constraints, while acknowledging that a first-principles derivation of the exact fraction lies beyond the present scope. This parameter can be further tested with future multi-messenger observations. revision: partial
Referee: [§3] §3 (mock catalog): the catalog is constructed to be consistent with the same X-ray and mid-IR AGN samples that set the luminosity function and thus the overall normalization. This introduces circularity between the input population statistics and the output diffuse flux, weakening the claim that the model independently explains the IceCube measurements.

Authors: We disagree that circularity is present. The X-ray and mid-IR luminosity functions are fixed observational inputs that describe the demographics of the AGN population. The mock catalog samples sources from these observed distributions and then applies our reconnection-based neutrino emission model to compute the diffuse flux. This is the standard population-synthesis approach and allows a prediction of the neutrino background from the known AGN population; the overall normalization is set by the model's microphysical parameters (σ_p distribution, etc.), not by re-fitting the catalog. We will add a clarifying paragraph in §3 to emphasize this distinction between input demographics and model-derived emission. revision: partial

Circularity Check

1 steps flagged

10% high-σ_p subpopulation tuned to normalize the predicted neutrino flux

specific steps

fitted input called prediction [Abstract]
"Coronal emission satisfactorily explains the most recent IceCube measurements of the diffuse neutrino flux up to energies of ∼1 PeV, provided that ∼10% of the AGN coronae have σ_p ∼10^5, while the rest are distributed over a range of lower magnetizations."

The 10% fraction is introduced as a free parameter whose value is selected so that the integrated flux from the mock catalog matches the observed IceCube normalization. Because neutrino output scales strongly with σ_p (via higher E_p,max and pγ optical depth), this choice directly forces the result rather than deriving it from independent data or first-principles constraints on the magnetization distribution.

full rationale

The paper builds neutrino templates over a grid of σ_p, L_X and M_BH, then couples them to a mock AGN catalog normalized to the same X-ray and mid-IR samples that set the overall population. The central claim that coronal reconnection 'satisfactorily explains' the IceCube flux up to ~1 PeV is made only after assigning ~10% of sources σ_p~10^5; this fraction is chosen to set the absolute normalization because neutrino luminosity rises steeply with σ_p. No separate observational or simulation constraint on the high-magnetization tail is supplied, so the reported agreement is achieved by construction via the adjustable subpopulation fraction.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on three fitted or chosen quantities: the 10% high-magnetization fraction, the specific value σ_p = 10^5 for that subpopulation, and the assumed reconnection acceleration efficiency. Standard astrophysical assumptions (reconnection as the acceleration mechanism, X-ray photon field as target, standard AGN demographics) are invoked without new justification.

free parameters (2)

high-magnetization fraction
Chosen at ~10% to reproduce the IceCube flux normalization
σ_p for high-magnetization subpopulation
Set to 10^5 to reach the required neutrino energies before cooling cuts in

axioms (2)

domain assumption Magnetic reconnection efficiently accelerates protons to relativistic energies in AGN coronae
Invoked to justify the proton population that produces neutrinos
domain assumption The mock AGN catalog accurately represents the true distribution of coronal X-ray luminosities and black-hole masses at z=0-4
Used to integrate single-source templates into the diffuse flux

pith-pipeline@v0.9.0 · 5587 in / 1586 out tokens · 42481 ms · 2026-05-16T17:58:26.581241+00:00 · methodology

discussion (0)

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Particle Acceleration, Coronal Neutrino Production, and the Diffuse Extragalactic Neutrino Background from Supermassive Black Holes
astro-ph.HE 2026-05 unverdicted novelty 6.0

The cosmologically integrated neutrino emission from supermassive black hole coronae in Seyfert galaxies can account for the sub-PeV diffuse extragalactic neutrino flux observed by IceCube.
Turbulent AGN coronae as the origin of diffuse neutrinos up to PeV energies
astro-ph.HE 2026-02 unverdicted novelty 6.0

A population of AGN coronae with magnetization parameters spanning up to σ ~ 10 can reproduce the entire observed diffuse neutrino flux from TeV to PeV energies.