pith. machine review for the scientific record. sign in

arxiv: 2602.20969 · v2 · submitted 2026-02-24 · 🌌 astro-ph.HE

Recognition: no theorem link

Turbulent AGN coronae as the origin of diffuse neutrinos up to PeV energies

Authors on Pith no claims yet

Pith reviewed 2026-05-15 19:50 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords AGN coronaediffuse neutrinosturbulence accelerationIceCubePeV neutrinosmagnetization parameterSeyfert galaxies
0
0 comments X

The pith

AGN coronae with high magnetization can explain the full IceCube diffuse neutrino flux up to PeV energies.

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

The paper shows that turbulence acceleration of protons inside the coronae of active galactic nuclei can produce neutrinos across the entire observed range if enough AGNs have coronae with magnetization parameter sigma of order 1 to 10. Earlier calculations found that typical low-magnetization coronae cap neutrino energies at tens of TeV, requiring extra source populations for the higher-energy events seen by IceCube. Allowing a broad distribution of magnetization values across the AGN population raises the maximum proton energy and hardens the spectrum enough to reach PeV neutrinos while also reproducing the spectral peak near 30 TeV. The result holds as long as the distribution extends to high sigma and is not too steeply falling at the upper end.

Core claim

A population of AGN coronae with a sufficiently wide range of magnetization parameters sigma approximately 1-10 produces the entire diffuse neutrino spectrum measured by IceCube up to PeV energies through turbulence acceleration, and the fit is insensitive to the precise shape of the sigma distribution provided the high-magnetization tail is populated.

What carries the argument

Turbulence acceleration of protons in magnetized AGN coronae, where the magnetization parameter sigma sets both the maximum proton energy and the hardness of the accelerated spectrum.

If this is right

  • The model naturally produces the observed peak in the diffuse neutrino spectrum near 30 TeV without extra tuning.
  • No separate high-energy source class is required to explain neutrinos above 100 TeV.
  • The accumulated flux is dominated by the subset of AGNs whose coronae reach sigma of several to ten.

Where Pith is reading between the lines

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

  • Individual detections such as the >100 TeV neutrino from NGC 7469 become representative rather than exceptional once the population includes high-sigma coronae.
  • Correlations between neutrino arrival directions and X-ray bright Seyfert galaxies could be used to measure the actual fraction of high-magnetization coronae.
  • The same coronae would accelerate cosmic-ray protons to energies that might also contribute to the extragalactic gamma-ray background through photopion or proton synchrotron processes.

Load-bearing premise

A substantial fraction of the overall AGN population has coronae with magnetization parameter sigma between roughly 1 and 10.

What would settle it

A sharp cutoff in the diffuse neutrino spectrum well below 1 PeV, or a failure to detect neutrinos above 100 TeV from a statistically large sample of Seyfert galaxies whose X-ray and radio properties indicate high corona magnetization.

Figures

Figures reproduced from arXiv: 2602.20969 by Qi-Rui Yang, Ruo-Yu Liu, Xiang-Yu Wang.

Figure 1
Figure 1. Figure 1: shows that highly magnetized coronae with lower X-ray luminosities can accelerate protons to higher en￾ergies, with the maximum proton energy constrained by the Hillas condition. The dark region with σ ≲ 0.04 on the very left side of [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Diffuse all-flavor neutrino fluxes from the AGN coronae assuming that the magnetization distribution follows the log-normal distribution with parameters µ = −1 and s = 1. The colored dashed curves denote the contributions from AGNs in the different ranges of the X-ray luminosity, while the black solid curve represents the sum of them. The blue and black data points represent IceCube observations, while the… view at source ↗
Figure 3
Figure 3. Figure 3: Diffuse all-flavor neutrino fluxes from the AGN population assuming power-law distributions for the magne￾tization with different power-law indices, as shown in col￾ored dot-dashed curves. For comparison, we also plot the neutrino fluxes from the AGN coronae assuming that the magnetization distribution follows the log-normal distribu￾tion, as shown by the solid black curve. The others keep the same with Fi… view at source ↗
Figure 4
Figure 4. Figure 4: Same as the Fig.2 but characterizing the con￾tribution of each magnetization range. The colored dashed curves denote the contributions from AGNs in the different ranges of the magnetization σ, while the black solid curve represents the sum of them. The population of relatively low X-ray luminosity AGNs with LX ∼ 1041 − 1042 erg s−1 constitutes the primary contribution to the diffuse neutrino flux above ∼ 1… view at source ↗
read the original abstract

It has been shown that the turbulence acceleration in AGN coronae can account for 1-10 TeV neutrinos from some AGNs, such as the Seyfert galaxy NGC 1068. Based on this, there are attempts to explain the diffuse neutrinos observed by IceCube with the accumulated contribution from a population of AGNs, but it is found that the maximum neutrino energy is less than tens of TeV, and as a result, additional source classes are needed to explain the high-energy component above this energy. Recently, motivated by the detection of $>100$ TeV neutrinos from the Seyfert galaxy NGC 7469, it was shown that the turbulence acceleration in the corona can explain $>$100 TeV neutrinos given a larger magnetization parameter ($\sigma\sim 1$) in the corona, which leads to a larger maximum proton energy and a hard proton spectrum. In this paper, we extend this assumption to the population of AGNs and study whether the population of AGNs with a wide range of magnetization can explain the entire diffuse neutrino flux. We find that AGN coronae could account for the diffuse neutrinos up to PeV energies if a significant fraction of AGNs have magnetizations as large as $\sigma\sim 1-10$. This conclusion is insensitive to the shape of the magnetization parameter distribution as long as the range of the magnetization parameter is sufficiently wide and the distribution is flat towards high magnetization. Interestingly, this model can also explain the peak of the diffuse neutrino spectrum at $\sim30$ TeV.

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.

Referee Report

3 major / 2 minor

Summary. The manuscript claims that turbulent particle acceleration in AGN coronae can explain the full IceCube diffuse neutrino spectrum up to PeV energies (including the observed peak near 30 TeV) provided a significant fraction of the AGN population has coronae with magnetization parameters σ ∼ 1–10; the result is stated to be insensitive to the precise shape of the σ distribution so long as it is wide and flat at the high-σ end.

Significance. If the central assumption can be independently justified, the result would be significant: it would unify the entire observed neutrino flux (TeV to PeV) under a single source class and mechanism, removing the need to invoke separate populations for the high-energy tail. The work correctly extends the turbulence-acceleration framework previously applied to NGC 7469 and NGC 1068 to the population level and demonstrates that the spectral shape can be reproduced under the stated conditions.

major comments (3)
  1. [§4] §4 (population synthesis): the magnetization distribution is constructed to be wide and flat at high σ precisely so that the integrated flux and high-energy cutoff match IceCube data; the PeV component is therefore a direct consequence of the chosen parameter range rather than an independent prediction.
  2. [§2.2–2.3] §2.2–2.3 (model assumptions): no derivation from MHD corona simulations, X-ray polarization measurements, or variability statistics is supplied to justify why a non-negligible fraction of AGNs should reach σ ∼ 1–10 rather than the lower values (σ ≪ 1) expected from equipartition; the prevalence of high-σ coronae remains an unconstrained population statistic.
  3. [§5] §5 (results): quantitative outputs for the neutrino spectrum are shown only for the tuned distribution; the claimed insensitivity to distribution shape is not demonstrated by explicit comparison of multiple functional forms with the same high-σ tail.
minor comments (2)
  1. [Abstract] The abstract and §1 could more explicitly state the minimum fraction of high-σ sources required to match the observed normalization.
  2. [Figures] Figure captions should clarify the exact functional forms and parameter ranges used for the magnetization distributions shown in the plots.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We have revised the manuscript to address the concerns raised, particularly by adding explicit comparisons and expanded discussion of model assumptions. Our point-by-point responses follow.

read point-by-point responses
  1. Referee: [§4] §4 (population synthesis): the magnetization distribution is constructed to be wide and flat at high σ precisely so that the integrated flux and high-energy cutoff match IceCube data; the PeV component is therefore a direct consequence of the chosen parameter range rather than an independent prediction.

    Authors: We acknowledge that the distribution parameters, including the width and flat high-σ tail, are chosen such that the integrated flux matches IceCube observations. However, the upper range σ ∼ 1–10 is not arbitrary: it is directly motivated by the turbulence-acceleration modeling of NGC 7469, where σ ∼ 1 is required to produce neutrinos above 100 TeV. The population-level calculation then tests whether a distribution that includes such objects (as required by at least one detected source) can unify the entire TeV–PeV spectrum. We have revised §4 to clarify this motivation from individual sources and to emphasize that the spectral extension to PeV energies follows from the acceleration physics at high magnetization rather than from fine-tuning alone. revision: partial

  2. Referee: [§2.2–2.3] §2.2–2.3 (model assumptions): no derivation from MHD corona simulations, X-ray polarization measurements, or variability statistics is supplied to justify why a non-negligible fraction of AGNs should reach σ ∼ 1–10 rather than the lower values (σ ≪ 1) expected from equipartition; the prevalence of high-σ coronae remains an unconstrained population statistic.

    Authors: We agree that the manuscript does not derive the required fraction of high-σ coronae from first-principles MHD simulations or observational statistics. The assumption is motivated by the need to explain the >100 TeV neutrinos detected from NGC 7469. In the revised manuscript we have expanded §2.3 with a new paragraph that discusses possible physical pathways to elevated magnetization (e.g., localized magnetic-flux accumulation in reconnection layers, as suggested by some corona MHD simulations) and notes that future X-ray polarization observations can provide independent constraints. We present the high-σ fraction as a testable prediction rather than a derived result. revision: partial

  3. Referee: [§5] §5 (results): quantitative outputs for the neutrino spectrum are shown only for the tuned distribution; the claimed insensitivity to distribution shape is not demonstrated by explicit comparison of multiple functional forms with the same high-σ tail.

    Authors: We thank the referee for this observation. The revised §5 now includes a new figure that explicitly compares the diffuse neutrino spectra obtained with three different σ distributions (uniform, power-law with index −1, and log-normal), each normalized to have the same high-σ tail extending to σ = 10. The resulting fluxes agree to within ∼20 % above 10 TeV, including the location of the ∼30 TeV peak and the extension to PeV energies, thereby demonstrating the claimed insensitivity to the precise functional form. revision: yes

Circularity Check

1 steps flagged

Magnetization distribution chosen to reproduce IceCube flux up to PeV

specific steps
  1. fitted input called prediction [Abstract]
    "We find that AGN coronae could account for the diffuse neutrinos up to PeV energies if a significant fraction of AGNs have magnetizations as large as σ∼1-10. This conclusion is insensitive to the shape of the magnetization parameter distribution as long as the range of the magnetization parameter is sufficiently wide and the distribution is flat towards high magnetization."

    The wide range and flat high-σ tail are selected so the integrated flux matches IceCube data; the PeV component and overall spectrum shape are therefore direct consequences of the fitted population distribution rather than an independent prediction from corona physics.

full rationale

The central claim extends a turbulence-acceleration model (previously applied to individual sources) to the full AGN population by positing a wide, flat distribution of σ reaching 1–10. This population assumption directly determines both the normalization and the high-energy cutoff of the cumulative neutrino spectrum, with the paper noting the result is insensitive to exact shape provided the high-σ tail is present. No independent constraint from MHD simulations or observations is supplied to justify the required fraction of high-σ sources; the match to data therefore follows from the input choice rather than an independent derivation.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model relies on a tunable distribution of the magnetization parameter σ across the AGN population; this distribution is the primary free parameter introduced to match the observed neutrino spectrum.

free parameters (1)
  • magnetization parameter distribution
    The range and shape of σ values (flat toward high σ) are chosen so the integrated neutrino output reproduces the IceCube flux up to PeV.
axioms (1)
  • domain assumption Turbulent acceleration in AGN coronae produces neutrinos via proton-photon or proton-gas interactions
    Inherited from earlier single-source studies of NGC 1068 and NGC 7469.

pith-pipeline@v0.9.0 · 5581 in / 1308 out tokens · 18831 ms · 2026-05-15T19:50:31.682136+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

65 extracted references · 65 canonical work pages · 4 internal anchors

  1. [1]

    2021, PhRvD, 104, 022002, doi: 10.1103/PhysRevD.104.022002

    Abbasi, R., Ackermann, M., Adams, J., et al. 2021, PhRvD, 104, 022002, doi: 10.1103/PhysRevD.104.022002

  2. [2]

    2022, ApJ, 928, 50, doi: 10.3847/1538-4357/ac4d29

    Abbasi, R., Ackermann, M., Adams, J., et al. 2022, ApJ, 928, 50, doi: 10.3847/1538-4357/ac4d29

  3. [3]

    2025a, arXiv e-prints, arXiv:2507.22233, doi: 10.48550/arXiv.2507.22233 10

    Abbasi, R., Ackermann, M., Adams, J., et al. 2025a, arXiv e-prints, arXiv:2507.22233, doi: 10.48550/arXiv.2507.22233 10

  4. [4]

    2025b, arXiv e-prints, arXiv:2510.13403

    Abbasi, R., Ackermann, M., Adams, J., et al. 2025b, arXiv e-prints, arXiv:2510.13403. https://arxiv.org/abs/2510.13403

  5. [5]

    2025c, ApJ, 988, 141, doi: 10.3847/1538-4357/addd05

    Abbasi, R., Ackermann, M., Adams, J., et al. 2025c, ApJ, 988, 141, doi: 10.3847/1538-4357/addd05

  6. [6]

    B., et al

    Ackermann, M., Ajello, M., Atwood, W. B., et al. 2012, ApJ, 750, 3, doi: 10.1088/0004-637X/750/1/3

  7. [7]

    2024, JCAP, 2024, 075, doi: 10.1088/1475-7516/2024/09/075

    Ambrosone, A. 2024, JCAP, 2024, 075, doi: 10.1088/1475-7516/2024/09/075

  8. [8]

    2025, arXiv e-prints, arXiv:2507.06002, doi: 10.48550/arXiv.2507.06002

    Basu, V., Balagopal V., A., & Karle, A. 2025, arXiv e-prints, arXiv:2507.06002, doi: 10.48550/arXiv.2507.06002

  9. [9]

    Beloborodov, A. M. 2017, ApJ, 850, 141, doi: 10.3847/1538-4357/aa8f4f

  10. [10]

    2022, PhRvD, 106, 023028, doi: 10.1103/PhysRevD.106.023028

    Bresci, V., Lemoine, M., Gremillet, L., et al. 2022, PhRvD, 106, 023028, doi: 10.1103/PhysRevD.106.023028

  11. [11]

    2019, ApJ, 886, 122, doi: 10.3847/1538-4357/ab4c33

    Comisso, L., & Sironi, L. 2019, ApJ, 886, 122, doi: 10.3847/1538-4357/ab4c33

  12. [12]

    Studying the mirror acceleration via kinetic simulations of relativistic plasma turbulence

    Das, S., Xu, S., & Nättilä, J. 2025, arXiv e-prints, arXiv:2506.04212, doi: 10.48550/arXiv.2506.04212 Di Matteo, T., Blackman, E. G., & Fabian, A. C. 1997, MNRAS, 291, L23, doi: 10.1093/mnras/291.1.L23

  13. [13]

    Eichmann, B., Oikonomou, F., Salvatore, S., Dettmar, R.-J., & Tjus, J. B. 2022, ApJ, 939, 43, doi: 10.3847/1538-4357/ac9588

  14. [14]

    C., Lohfink, A., Kara, E., et al

    Fabian, A. C., Lohfink, A., Kara, E., et al. 2015, MNRAS, 451, 4375, doi: 10.1093/mnras/stv1218

  15. [15]

    Fiorillo, D. F. G., Comisso, L., Peretti, E., Petropoulou, M., & Sironi, L. 2024, ApJ, 974, 75, doi: 10.3847/1538-4357/ad7021

  16. [16]

    Fiorillo, D. F. G., Comisso, L., Peretti, E., Petropoulou, M., & Sironi, L. 2025, ApJ, 989, 215, doi: 10.3847/1538-4357/adec9c Grošelj, D., Hakobyan, H., Beloborodov, A. M., Sironi, L., & Philippov, A. 2024, PhRvL, 132, 085202, doi: 10.1103/PhysRevLett.132.085202

  17. [17]

    1991, ApJL, 380, L51, doi: 10.1086/186171

    Haardt, F., & Maraschi, L. 1991, ApJL, 380, L51, doi: 10.1086/186171

  18. [18]

    2023, arXiv e-prints, arXiv:2305.07086, doi: 10.48550/arXiv.2305.07086

    Halzen, F. 2023, arXiv e-prints, arXiv:2305.07086, doi: 10.48550/arXiv.2305.07086

  19. [19]

    2022, arXiv e-prints, arXiv:2202.00694, doi: 10.48550/arXiv.2202.00694

    Halzen, F., & Kheirandish, A. 2022, arXiv e-prints, arXiv:2202.00694, doi: 10.48550/arXiv.2202.00694

  20. [20]

    F., Richards, G

    Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 731, doi: 10.1086/509629 IceCube Collaboration. 2013, Science, 342, 1242856, doi: 10.1126/science.1242856 IceCube Collaboration, Abbasi, R., Ackermann, M., et al. 2022, Science, 378, 538, doi: 10.1126/science.abg3395

  21. [21]

    2022, arXiv e-prints, arXiv:2207.02097, doi: 10.48550/arXiv.2207.02097

    Inoue, S., Cerruti, M., Murase, K., & Liu, R.-Y. 2022, arXiv e-prints, arXiv:2207.02097, doi: 10.48550/arXiv.2207.02097

  22. [22]

    2018, ApJ, 869, 114, doi: 10.3847/1538-4357/aaeb95

    Inoue, Y., & Doi, A. 2018, ApJ, 869, 114, doi: 10.3847/1538-4357/aaeb95

  23. [23]

    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

  24. [24]

    2019, ApJ, 880, 40, doi: 10.3847/1538-4357/ab2715

    Inoue, Y., Khangulyan, D., Inoue, S., & Doi, A. 2019, ApJ, 880, 40, doi: 10.3847/1538-4357/ab2715

  25. [25]

    M., & Davis, S

    Jiang, Y.-F., Stone, J. M., & Davis, S. W. 2014, ApJ, 784, 169, doi: 10.1088/0004-637X/784/2/169

  26. [26]

    Karavola, D., Petropoulou, M., Fiorillo, D. F. G., et al. 2026, arXiv e-prints, arXiv:2601.01533, doi: 10.48550/arXiv.2601.01533

  27. [27]

    2010, PASJ, 62, 621, doi: 10.1093/pasj/62.3.621

    Kawabata, R., & Mineshige, S. 2010, PASJ, 62, 621, doi: 10.1093/pasj/62.3.621

  28. [28]

    R., Aharonian, F

    Kelner, S. R., Aharonian, F. A., & Bugayov, V. V. 2006, PhRvD, 74, 034018, doi: 10.1103/PhysRevD.74.034018

  29. [29]

    Kheirandish, A., Murase, K., & Kimura, S. S. 2021, ApJ, 922, 45, doi: 10.3847/1538-4357/ac1c77

  30. [30]

    S., Murase, K., & Mészáros, P

    Kimura, S. S., Murase, K., & Mészáros, P. 2021, Nature Communications, 12, 5615, doi: 10.1038/s41467-021-25111-7

  31. [31]

    S., Murase, K., & Toma, K

    Kimura, S. S., Murase, K., & Toma, K. 2015, ApJ, 806, 159, doi: 10.1088/0004-637X/806/2/159

  32. [32]

    2022, Annual Review of Nuclear and Particle Science, 72, 365, doi: 10.1146/annurev-nucl-011122-061547

    Kurahashi, N., Murase, K., & Santander, M. 2022, Annual Review of Nuclear and Particle Science, 72, 365, doi: 10.1146/annurev-nucl-011122-061547

  33. [33]

    2021, PhRvD, 104, 063020, doi: 10.1103/PhysRevD.104.063020

    Lemoine, M. 2021, PhRvD, 104, 063020, doi: 10.1103/PhysRevD.104.063020

  34. [34]

    2022, PhRvL, 129, 215101, doi: 10.1103/PhysRevLett.129.215101

    Lemoine, M. 2022, PhRvL, 129, 215101, doi: 10.1103/PhysRevLett.129.215101

  35. [35]

    2024, PhRvD, 109, 063006, doi: 10.1103/PhysRevD.109.063006

    Lemoine, M., Murase, K., & Rieger, F. 2024, PhRvD, 109, 063006, doi: 10.1103/PhysRevD.109.063006

  36. [36]

    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

  37. [37]

    Liska, M. T. P., Musoke, G., Tchekhovskoy, A., Porth, O., & Beloborodov, A. M. 2022, ApJL, 935, L1, doi: 10.3847/2041-8213/ac84db

  38. [38]

    F., Mineshige, S., & Shibata, K

    Liu, B. F., Mineshige, S., & Shibata, K. 2002, ApJL, 572, L173, doi: 10.1086/341877

  39. [39]

    2016, MNRAS, 456, L94, doi: 10.1093/mnrasl/slv178

    Marinucci, A., Bianchi, S., Matt, G., et al. 2016, MNRAS, 456, L94, doi: 10.1093/mnrasl/slv178

  40. [40]
  41. [41]

    2023, ApJ, 944, 122, doi: 10.3847/1538-4357/acaefe

    Meringolo, C., Cruz-Osorio, A., Rezzolla, L., & Servidio, S. 2023, ApJ, 944, 122, doi: 10.3847/1538-4357/acaefe

  42. [42]

    Merloni, A., Di Matteo, T., & Fabian, A. C. 2000, MNRAS, 318, L15, doi: 10.1046/j.1365-8711.2000.03943.x 11

  43. [43]

    Merloni, A., & Fabian, A. C. 2001, MNRAS, 321, 549, doi: 10.1046/j.1365-8711.2001.04060.x

  44. [44]

    Evaporation: The change from accretion via a thin disk to a coronal flow

    Meyer, F., Liu, B. F., & Meyer-Hofmeister, E. 2000, A&A, 361, 175, doi: 10.48550/arXiv.astro-ph/0007091

  45. [45]

    2022, ApJL, 941, L17, doi: 10.3847/2041-8213/aca53c

    Murase, K. 2022, ApJL, 941, L17, doi: 10.3847/2041-8213/aca53c

  46. [46]

    S., & Mészáros, P

    Murase, K., Kimura, S. S., & Mészáros, P. 2020a, PhRvL, 125, 011101, doi: 10.1103/PhysRevLett.125.011101

  47. [47]

    2020b, ApJ, 902, 108, doi: 10.3847/1538-4357/abb3c0

    Petropoulou, M. 2020b, ApJ, 902, 108, doi: 10.3847/1538-4357/abb3c0

  48. [48]

    2024, A&A, 684, L21, doi: 10.1051/0004-6361/202450025

    Henningsen, F. 2024, A&A, 684, L21, doi: 10.1051/0004-6361/202450025

  49. [49]

    2015, MNRAS, 452, 1877, doi: 10.1093/mnras/stv1467

    Padovani, P., Petropoulou, M., Giommi, P., & Resconi, E. 2015, MNRAS, 452, 1877, doi: 10.1093/mnras/stv1467

  50. [50]

    Pringle, J. E. 1981, ARA&A, 19, 137, doi: 10.1146/annurev.aa.19.090181.001033 Różańska, A., & Czerny, B. 2000, A&A, 360, 1170, doi: 10.48550/arXiv.astro-ph/0004158

  51. [51]

    2026, PhRvD, 113, 023019, doi: 10.1103/f66p-k6z9

    Saurenhaus, L., Capel, F., Oikonomou, F., & Buchner, J. 2026, PhRvD, 113, 023019, doi: 10.1103/f66p-k6z9

  52. [52]

    Sironi, L., & Beloborodov, A. M. 2020, ApJ, 899, 52, doi: 10.3847/1538-4357/aba622

  53. [53]

    2025, ApJ, 981, 103, doi: 10.3847/1538-4357/adb031

    Sommani, G., Franckowiak, A., Lincetto, M., & Dettmar, R.-J. 2025, ApJ, 981, 103, doi: 10.3847/1538-4357/adb031

  54. [54]

    W., Done, C., Salamon, M

    Stecker, F. W., Done, C., Salamon, M. H., & Sommers, P. 1991, PhRvL, 66, 2697, doi: 10.1103/PhysRevLett.66.2697

  55. [55]

    Svensson, R., & Zdziarski, A. A. 1994, ApJ, 436, 599, doi: 10.1086/174934

  56. [56]

    Thorne, K. S. 1974, ApJ, 191, 507, doi: 10.1086/152991

  57. [57]

    J., et al

    Trakhtenbrot, B., Ricci, C., Koss, M. J., et al. 2017, MNRAS, 470, 800, doi: 10.1093/mnras/stx1117

  58. [58]

    Watson, M. G. 2014, ApJ, 786, 104, doi: 10.1088/0004-637X/786/2/104

  59. [59]

    A., Werner, G

    Wong, K., Zhdankin, V., Uzdensky, D. A., Werner, G. R., & Begelman, M. C. 2020, ApJL, 893, L7, doi: 10.3847/2041-8213/ab8122

  60. [60]

    W., Zhdankin, V., Uzdensky, D

    Wong, K. W., Zhdankin, V., Uzdensky, D. A., Werner, G. R., & Begelman, M. C. 2025, MNRAS, 543, 1842, doi: 10.1093/mnras/staf1589

  61. [61]

    2025, ApJ, 995, 166, doi: 10.3847/1538-4357/ae1ba6

    Lemoine, M. 2025, ApJ, 995, 166, doi: 10.3847/1538-4357/ae1ba6

  62. [62]

    Yuan, C., Fiorillo, D. F. G., Petropoulou, M., & Liu, Q. 2025, arXiv e-prints, arXiv:2508.08233, doi: 10.48550/arXiv.2508.08233

  63. [63]

    Begelman, M. C. 2018, ApJL, 867, L18, doi: 10.3847/2041-8213/aae88c

  64. [64]

    2025, arXiv e-prints, arXiv:2511.16869, doi: 10.48550/arXiv.2511.16869

    Zhou, S., Sun, M., Mou, G., et al. 2025, arXiv e-prints, arXiv:2511.16869, doi: 10.48550/arXiv.2511.16869

  65. [65]

    Miller, J. M. 2005, ApJ, 618, 832, doi: 10.1086/426071