Recognition: 2 theorem links
· Lean TheoremCan magnetic reconnection power neutrino emission from AGN coronae?
Pith reviewed 2026-05-12 04:41 UTC · model grok-4.3
The pith
Reconnection in transrelativistic black hole coronae accelerates protons to tens of PeV, producing observed neutrinos from NGC 1068.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Across the one-parameter family of coronal conditions, the proton magnetization is transrelativistic with σ_p ∼ 0.3. In this regime, repeated encounters with intermittent reconnecting current sheets can energize suprathermal protons up to tens of PeV before photomeson cooling limits further acceleration. These injected particles may then be further processed by stochastic interactions with the turbulent cascade. Motivated by PIC simulations of strong turbulence at comparable magnetization, the adopted nonthermal proton spectrum produces a predicted TeV spectral shape that is broadly consistent with NGC 1068 without fitting the proton spectral slope.
What carries the argument
The one-parameter family of coronal conditions obtained by matching neutrino and X-ray luminosities plus Thomson depth, which fixes the proton magnetization (magnetic energy density over proton rest-mass energy density) near 0.3 and thereby permits repeated acceleration in reconnecting current sheets.
If this is right
- Suprathermal protons reach energies of tens of PeV before photomeson cooling becomes dominant.
- The injected particles undergo additional stochastic acceleration in the turbulent cascade.
- The resulting neutrino emission in the TeV band matches the observed spectrum of NGC 1068.
- Reconnection supplies the nonthermal protons needed for high-energy neutrino production in AGN coronae.
- The mechanism operates across the entire one-parameter family of allowed coronal parameters.
Where Pith is reading between the lines
- The same modeling framework could be applied to other neutrino-bright AGN to test whether their coronae fall into the same magnetization range.
- If the acceleration limit is general, reconnection in low-beta coronae may contribute to the diffuse astrophysical neutrino background.
- Multi-messenger observations that independently constrain coronal magnetic field or density would provide a direct test of the derived one-parameter family.
- Targeted turbulence simulations at exactly σ_p ∼ 0.3 and the inferred plasma beta would tighten the spectral prediction.
Load-bearing premise
The nonthermal proton spectrum shape is taken directly from particle-in-cell simulations of strong turbulence at comparable magnetization without further verification against the specific coronal conditions.
What would settle it
A measured neutrino spectrum from NGC 1068 whose TeV shape deviates significantly from the predicted distribution, or a direct constraint showing that coronal protons cannot reach tens of PeV because photomeson cooling is stronger than assumed.
Figures
read the original abstract
We investigate whether reconnection of small-scale current sheets in transrelativistic supermassive black hole (SMBH) coronae can supply the nonthermal protons needed for high-energy neutrino emission, using NGC 1068 as a test case. We model the corona as a strongly turbulent, low-$\beta$, collisionless hydrogen plasma with characteristic size $r_{\rm co}$, magnetic field strength $B$, proton density $n_p$, and radiation energy density $u_{\rm rad}$. Combining the observed IceCube-band neutrino luminosity with the X-ray luminosity and Thomson optical depth reduces these coronal quantities to a one-parameter family. Across this family, the proton magnetization $\sigma_p \equiv B^2/(4\pi n_p m_p c^2)$ is transrelativistic with $\sigma_p \sim 0.3$. In this regime, we show that repeated encounters with intermittent reconnecting current sheets can energize suprathermal protons up to tens of PeV before photomeson cooling limits further acceleration. These injected particles may then be further processed by stochastic interactions with the turbulent cascade. Motivated by PIC simulations of strong turbulence at comparable magnetization, we adopt a nonthermal proton spectrum with an independently specified index and find that the predicted TeV spectral shape is broadly consistent with NGC~1068 without fitting the proton spectral slope.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper investigates whether magnetic reconnection in transrelativistic, low-β AGN coronae can accelerate protons to produce the observed high-energy neutrinos, taking NGC 1068 as a test case. The corona is modeled with parameters r_co, B, n_p, and u_rad; combining the IceCube neutrino luminosity, X-ray luminosity, and Thomson optical depth reduces these to a one-parameter family in which σ_p ≡ B²/(4π n_p m_p c²) ≈ 0.3. In this regime, repeated encounters with reconnecting current sheets are shown to energize suprathermal protons to tens of PeV before photomeson cooling dominates; a nonthermal proton spectrum is then adopted from external PIC simulations of strong turbulence at comparable magnetization, yielding a predicted TeV neutrino spectral shape that is broadly consistent with NGC 1068 observations without fitting the spectral index.
Significance. If the central claims are robust, the work supplies a concrete, observationally constrained mechanism linking reconnection-driven acceleration in SMBH coronae to IceCube neutrinos. The reduction of four coronal quantities to a one-parameter family via multiple independent observables is a methodological strength that makes the model more falsifiable than fully free-parameter approaches. The result, if confirmed, would strengthen the case for AGN coronae as neutrino sources and motivate targeted PIC or hybrid simulations of transrelativistic reconnection under coronal radiation fields.
major comments (3)
- [Derivation of the one-parameter family (abstract and §3)] The one-parameter family is constructed using the observed IceCube-band neutrino luminosity as a constraint. The subsequent claim that the predicted TeV spectral shape is 'broadly consistent without fitting' therefore depends on the same luminosity datum that fixed the normalization; an explicit demonstration that the shape agreement survives when the neutrino luminosity is treated as a prediction (rather than an input) would remove the appearance of circularity.
- [Adoption of proton spectrum and comparison to NGC 1068 (abstract and §4)] The nonthermal proton spectrum (index and form) is imported directly from external PIC simulations of strong turbulence at σ_p ∼ 0.3 and described as 'independently specified.' The coronal conditions derived from the family include a specific u_rad and low-β plasma that may alter injection, reconnection rate, or stochastic processing relative to the cited PIC setups; without a dedicated comparison or additional simulation at the exact (r_co, B, n_p, u_rad) values of the family, this step remains load-bearing for the no-fitting claim.
- [Acceleration and cooling limits (abstract and §3.2)] The assertion that protons reach tens of PeV before photomeson cooling is central to the viability of the mechanism. The abstract and main text provide no error propagation on the remaining free parameter of the family, nor an explicit scan showing that the maximum energy stays above the photomeson threshold across the full allowed range; such a check is required to establish robustness.
minor comments (2)
- The abstract would be clearer if it stated the numerical range spanned by the remaining free parameter after the luminosity and optical-depth constraints are applied.
- Provide the exact magnetization values and turbulence regimes of the referenced PIC runs so readers can judge the degree of overlap with the derived coronal family.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments, which have prompted us to clarify key aspects of our analysis. We respond point by point to the major comments below.
read point-by-point responses
-
Referee: The one-parameter family is constructed using the observed IceCube-band neutrino luminosity as a constraint. The subsequent claim that the predicted TeV spectral shape is 'broadly consistent without fitting' therefore depends on the same luminosity datum that fixed the normalization; an explicit demonstration that the shape agreement survives when the neutrino luminosity is treated as a prediction (rather than an input) would remove the appearance of circularity.
Authors: The observed neutrino luminosity is used solely to set the overall normalization (amplitude) of the predicted spectrum. The TeV spectral shape, however, is fixed by the independently adopted nonthermal proton distribution (whose index and form are taken from PIC simulations at matching magnetization) together with the energy-dependent photomeson cooling. These shape-determining ingredients are independent of the luminosity constraint. Across the one-parameter family, σ_p remains ≈0.3 while the remaining parameter (e.g., coronal size) varies; the resulting neutrino spectral shape stays consistent with NGC 1068 data. We will revise the text to explicitly separate normalization from shape and to state that the proton spectral index was not adjusted to fit the neutrino observations. revision: partial
-
Referee: The nonthermal proton spectrum (index and form) is imported directly from external PIC simulations of strong turbulence at σ_p ∼ 0.3 and described as 'independently specified.' The coronal conditions derived from the family include a specific u_rad and low-β plasma that may alter injection, reconnection rate, or stochastic processing relative to the cited PIC setups; without a dedicated comparison or additional simulation at the exact (r_co, B, n_p, u_rad) values of the family, this step remains load-bearing for the no-fitting claim.
Authors: The magnetization σ_p ≈ 0.3 is the dominant parameter controlling the nonthermal spectrum in the referenced turbulence simulations, and this value is reproduced throughout our one-parameter family. The low-β, transrelativistic regime is also consistent with those setups. While a dedicated simulation at the precise (r_co, B, n_p, u_rad) values would be desirable, it lies outside the scope of the present study. We will expand the discussion in the revised manuscript to compare the physical conditions more explicitly with the cited PIC works and to justify the applicability of the adopted spectrum. revision: partial
-
Referee: The assertion that protons reach tens of PeV before photomeson cooling is central to the viability of the mechanism. The abstract and main text provide no error propagation on the remaining free parameter of the family, nor an explicit scan showing that the maximum energy stays above the photomeson threshold across the full allowed range; such a check is required to establish robustness.
Authors: We agree that a quantitative robustness check is needed. In the revised manuscript we will add an explicit scan over the allowed range of the remaining free parameter, together with error propagation, demonstrating that the maximum proton energy remains above the photomeson threshold for the entire family. revision: yes
Circularity Check
No significant circularity; spectrum index externally adopted and consistency check is independent of fitted normalization
full rationale
The derivation reduces four coronal quantities to a one-parameter family via observed neutrino luminosity, X-ray luminosity, and Thomson depth, then computes σ_p ≈ 0.3 across the family. The nonthermal proton spectrum index is explicitly adopted from external PIC simulations of strong turbulence at comparable magnetization and is described as independently specified without fitting to the NGC 1068 data. The TeV spectral shape is then checked for broad consistency. Because the index and functional form originate outside the present observational constraints and are not adjusted to match the target spectrum, the consistency statement does not reduce to the input luminosities by construction. No self-citations, self-definitional steps, or renamings of known results appear in the provided text.
Axiom & Free-Parameter Ledger
free parameters (2)
- remaining coronal parameter after luminosity constraints
- proton spectral index
axioms (2)
- domain assumption The corona is a strongly turbulent, low-β, collisionless hydrogen plasma.
- domain assumption Photomeson cooling limits proton energy before further acceleration or escape.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Combining the observed IceCube-band neutrino luminosity with the X-ray luminosity and Thomson optical depth reduces these coronal quantities to a one-parameter family. Across this family, the proton magnetization sigma_p ~0.3.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Motivated by PIC simulations of strong turbulence at comparable magnetization, we adopt a nonthermal proton spectrum with an independently specified index
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
2020a, Physical Review Letters, 124, 10.1103/physrevlett.124.051103
Aartsen, M. G., Ackermann, M., Adams, J., et al. 2020, PhRvL, 124, 051103, doi: 10.1103/PhysRevLett.124.051103
- [2]
-
[3]
Bauer, F. E., Ar´ evalo, P., Walton, D. J., et al. 2015, The Astrophysical Journal, 812, 116, doi: 10.1088/0004-637X/812/2/116
-
[4]
Begelman, M. C., Rudak, B., & Sikora, M. 1990, ApJ, 362, 38, doi: 10.1086/169241
-
[5]
1934, Proceedings of the Royal Society of London Series A, 146, 83, doi: 10.1098/rspa.1934.0140
Bethe, H., & Heitler, W. 1934, Proceedings of the Royal Society of London Series A, 146, 83, doi: 10.1098/rspa.1934.0140
-
[6]
Bethe, H. A., & Maximon, L. C. 1954, Phys. Rev., 93, 768, doi: 10.1103/PhysRev.93.768
-
[7]
Blanco, C., Hooper, D., Linden, T., & Pinetti, E. 2025, Phys. Rev. D, 112, 123016, doi: 10.1103/wnjh-7nwp
-
[8]
2006, PhRvL, 96, 115002, doi: 10.1103/PhysRevLett.96.115002
Boldyrev, S. 2006, PhRvL, 96, 115002, doi: 10.1103/PhysRevLett.96.115002
-
[9]
2019, ApJ, 886, 122, doi: 10.3847/1538-4357/ab4c33 —
Comisso, L., & Sironi, L. 2019, ApJ, 886, 122, doi: 10.3847/1538-4357/ab4c33 —. 2022, ApJL, 936, L27, doi: 10.3847/2041-8213/ac8422
-
[10]
Dahlin, J. T., Drake, J. F., & Swisdak, M. 2017, Physics of Plasmas, 24, 092110, doi: 10.1063/1.4986211
-
[11]
2026, The Astrophysical Journal, 999, 261, doi: 10.3847/1538-4357/ae3d94
Davis, Z., Comisso, L., Haggerty, C., & N¨ attil¨ a, J. 2026, The Astrophysical Journal, 999, 261, doi: 10.3847/1538-4357/ae3d94
-
[12]
2020, Physical Review D, 102, doi: 10.1103/physrevd.102.023003
Demidem, C., Lemoine, M., & Casse, F. 2020, Physical Review D, 102, doi: 10.1103/physrevd.102.023003
-
[13]
1979, ApJ, 232, 106, doi: 10.1086/157269
Eichler, D. 1979, ApJ, 232, 106, doi: 10.1086/157269
-
[14]
Fiorillo, D. F. G., Comisso, L., Peretti, E., Petropoulou, M., & Sironi, L. 2024a, The Astrophysical Journal, 974, 75, doi: 10.3847/1538-4357/ad7021
-
[15]
Fiorillo, D. F. G., Petropoulou, M., Comisso, L., Peretti, E., & Sironi, L. 2024b, The Astrophysical Journal Letters, 961, L14, doi: 10.3847/2041-8213/ad192b
-
[16]
French, O., Guo, F., Zhang, Q., & Uzdensky, D. A. 2023, ApJ, 948, 19, doi: 10.3847/1538-4357/acb7dd
-
[17]
French, O., Werner, G. R., & Uzdensky, D. A. 2026, Journal of Plasma Physics, 92, E10, doi: 10.1017/S0022377825101189
-
[18]
Gaisser, T. K., Engel, R., & Resconi, E. 2016, Cosmic Rays and Particle Physics (Cambridge University Press)
work page 2016
-
[19]
2022, Physical Review Letters, 129, 265101, doi: 10.1103/physrevlett.129.265101
Goodbred, M., & Liu, Y.-H. 2022, Physical Review Letters, 129, 265101, doi: 10.1103/physrevlett.129.265101
-
[20]
Greenhill, L. J., Gwinn, C. R., Antonucci, R., & Barvainis, R. 1996, ApJL, 472, L21, doi: 10.1086/310346 Groˇ selj, D., Philippov, A., Beloborodov, A. M., &
-
[21]
2026, ApJ, 1001, 64, doi: 10.3847/1538-4357/ae50fc
Mushotzky, R. 2026, ApJ, 1001, 64, doi: 10.3847/1538-4357/ae50fc
-
[22]
2014, Physical Review Letters, 113, doi: 10.1103/physrevlett.113.155005
Guo, F., Li, H., Daughton, W., & Liu, Y.-H. 2014, Physical Review Letters, 113, doi: 10.1103/physrevlett.113.155005
-
[23]
2020, Physics of Plasmas, 27, 080501, doi: 10.1063/5.0012094
Guo, F., Liu, Y.-H., Li, X., et al. 2020, Physics of Plasmas, 27, 080501, doi: 10.1063/5.0012094
-
[24]
Heitler, W. 1954, Quantum theory of radiation (Oxford University (Clarendon Press)) IceCube Collaboration, Abbasi, R., et al. 2026a, The Astrophysical Journal Letters, 1000, L26, doi: 10.3847/2041-8213/ae4aad —. 2026b, The Astrophysical Journal Letters, 1000, L37, doi: 10.3847/2041-8213/ae4aac IceCube Collaboration, Aartsen, M. G., Ackermann, M., et al. 2...
-
[25]
2020, ApJL, 891, L33, doi: 10.3847/2041-8213/ab7661
Inoue, Y., Khangulyan, D., & Doi, A. 2020, ApJL, 891, L33, doi: 10.3847/2041-8213/ab7661
-
[26]
2025, Journal of Cosmology and Astroparticle Physics, 2025, 075, doi: 10.1088/1475-7516/2025/04/075
Sironi, L. 2025, Journal of Cosmology and Astroparticle Physics, 2025, 075, doi: 10.1088/1475-7516/2025/04/075
-
[27]
Kelner, S. R., & Aharonian, F. A. 2008, Physical Review D, 78, doi: 10.1103/physrevd.78.034013
-
[28]
2021, Physical Review D, 104, doi: 10.1103/physrevd.104.063020
Lemoine, M. 2021, Physical Review D, 104, doi: 10.1103/physrevd.104.063020
-
[29]
Lemoine, M., & Malkov, M. A. 2020, MNRAS, 499, 4972, doi: 10.1093/mnras/staa3131
-
[30]
Lemoine, M., Murase, K., & Rieger, F. 2024, Phys. Rev. D, 109, 063006, doi: 10.1103/PhysRevD.109.063006
-
[31]
2025, A&A, 697, A124, doi: 10.1051/0004-6361/202453296
Lemoine, M., & Rieger, F. 2025, A&A, 697, A124, doi: 10.1051/0004-6361/202453296
-
[32]
2025, Space Science Reviews, 221, 16, doi: 10.1007/s11214-025-01142-0
Liu, Y.-H., Hesse, M., Genestreti, K., et al. 2025, Space Science Reviews, 221, 16, doi: 10.1007/s11214-025-01142-0
-
[33]
Loureiro, N. F., & Boldyrev, S. 2017a, PhRvL, 118, 245101, doi: 10.1103/PhysRevLett.118.245101 —. 2017b, ApJ, 850, 182, doi: 10.3847/1538-4357/aa9754
-
[34]
N., Inoue, Y., Sentoku, Y., & Sano, T
Ly, M. N., Inoue, Y., Sentoku, Y., & Sano, T. 2026, Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos. https://arxiv.org/abs/2601.01999 MAGIC Collaboration, Acciari, V. A., Ansoldi, S., et al. 2019, The Astrophysical Journal, 883, 135, doi: 10.3847/1538-4357/ab3a51
-
[35]
Mallet, A., Schekochihin, A. A., & Chandran, B. D. G. 2017, MNRAS, 468, 4862, doi: 10.1093/mnras/stx670
-
[36]
Mandelstam, S. 1958, Phys. Rev., 112, 1344, doi: 10.1103/PhysRev.112.1344
-
[37]
Marinucci, A., Bianchi, S., Matt, G., et al. 2016, Monthly Notices of the Royal Astronomical Society: Letters, 456, L94, doi: 10.1093/mnrasl/slv178
-
[38]
Mbarek, R., Philippov, A., Chernoglazov, A., Levinson, A., & Mushotzky, R. 2024, Phys. Rev. D, 109, L101306, doi: 10.1103/PhysRevD.109.L101306
-
[39]
Favre, J. M. 2014, Astronomy & Astrophysics, 570, A111, doi: 10.1051/0004-6361/201424083
-
[40]
Murase, K., Kimura, S. S., & M´ esz´ aros, P. 2020, Physical Review Letters, 125, 011101, doi: 10.1103/PhysRevLett.125.011101
-
[41]
Murase, K., & Stecker, F. W. 2023, High-Energy Neutrinos from Active Galactic Nuclei (WORLD SCIENTIFIC), 483–540, doi: 10.1142/9789811282645 0010
-
[42]
1995, ApJ, 452, 710, doi: 10.1086/176343 OpenAI
Narayan, R., & Yi, I. 1995, ApJ, 452, 710, doi: 10.1086/176343
-
[43]
Ricci, C., Ho, L. C., Fabian, A. C., et al. 2018, Monthly Notices of the Royal Astronomical Society, 480, 1819, doi: 10.1093/mnras/sty1879
-
[44]
, archivePrefix = "arXiv", eprint =
Rowan, M. E., Sironi, L., & Narayan, R. 2017, The Astrophysical Journal, 850, 29, doi: 10.3847/1538-4357/aa9380
-
[45]
F., N¨ attil¨ a, J., & Zhdankin, V
Serrano, R. F., N¨ attil¨ a, J., & Zhdankin, V. 2025, Journal of Plasma Physics, 91, E50, doi: 10.1017/S002237782500011X
- [46]
-
[47]
Sikora, M., Kirk, J. G., Begelman, M. C., & Schneider, P. 1987, The Astrophysical Journal, 320, L81, doi: 10.1086/184980
-
[48]
2014, ApJL, 783, L21, doi: 10.1088/2041-8205/783/1/L21
Sironi, L., & Spitkovsky, A. 2014, ApJL, 783, L21, doi: 10.1088/2041-8205/783/1/L21
-
[49]
Stecker, F. W. 1968, Phys. Rev. Lett., 21, 1016, doi: 10.1103/PhysRevLett.21.1016
- [50]
-
[51]
Werner, G. R., Uzdensky, D. A., Cerutti, B., Nalewajko, K., & Begelman, M. C. 2016, The Astrophysical Journal Letters, 816, L8, doi: 10.3847/2041-8205/816/1/L8
-
[52]
Wong, K., Zhdankin, V., Uzdensky, D. A., Werner, G. R., & Begelman, M. C. 2020, The Astrophysical Journal, 893, L7, doi: 10.3847/2041-8213/ab8122
-
[53]
Wong, K. W., Zhdankin, V., Uzdensky, D. A., Werner, G. R., & Begelman, M. C. 2025, Monthly Notices of the Royal Astronomical Society, 543, 1842, doi: 10.1093/mnras/staf1589
-
[54]
2001, The Astrophysical Journal, 562, L63–L66, doi: 10.1086/337972
Zenitani, S., & Hoshino, M. 2001, The Astrophysical Journal, 562, L63–L66, doi: 10.1086/337972
-
[55]
Zhdankin, V., Boldyrev, S., & Uzdensky, D. A. 2016, Physics of Plasmas, 23, 055705, doi: 10.1063/1.4944820
-
[56]
Begelman, M. C. 2018, ApJL, 867, L18, doi: 10.3847/2041-8213/aae88c
-
[57]
Begelman, M. C. 2019, Physical Review Letters, 122, doi: 10.1103/physrevlett.122.055101
-
[58]
Zou, M., Hakobyan, H., Mbarek, R., et al. 2024, A New Particle Pusher with Hadronic Interactions for Modeling Multimessenger Emission from Compact Objects. https://arxiv.org/abs/2410.22781
-
[59]
Zyla, P. A., et al. 2020, PTEP, 2020, 083C01, doi: 10.1093/ptep/ptaa104
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.