pith. the verified trust layer for science. sign in

arxiv: 2509.09590 · v2 · submitted 2025-09-11 · 🌌 astro-ph.HE

Ultra-high energy event KM3-230213A as a cosmogenic neutrino in light of minimal UHECR flux models

Pith reviewed 2026-05-18 17:39 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords cosmogenic neutrinosultra-high energy cosmic raysKM3NeTTelescope Arrayneutrino fluxgamma-ray backgroundUHECR propagation
0
0 comments X p. Extension

The pith

The 220 PeV neutrino seen by KM3NeT can arise as a cosmogenic particle in minimal ultra-high energy cosmic ray models.

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

The paper tests whether the unusually energetic neutrino event KM3-230213A detected by KM3NeT originates from interactions of ultra-high energy cosmic rays with background photons during intergalactic travel. It adopts flux models fitted to Telescope Array data that assume light nuclear composition and a single source population evolving in the standard way with redshift. These models produce a cosmogenic neutrino flux that lies within roughly two standard deviations of both the KM3NeT-only measurement and the combined neutrino observatory data. The associated gamma-ray flux also remains below current limits from Fermi-LAT and ultra-high energy gamma-ray observations. A sympathetic reader would care because this offers a conventional astrophysical account that avoids invoking new physics or rare source classes to explain the event.

Core claim

We show that the predictions of the cosmogenic neutrino flux in these models are consistent with the measurements of the KM3NeT-only and with that of the 'global neutrino observatory' at approximately 2σ level. Notably, this result is achieved in a minimal version of the UHECR flux models that assume one source population with a standard cosmological evolution. We also estimate the corresponding cosmogenic gamma-ray flux and show that it is consistent with Fermi-LAT IGRB measurements and UHE gamma-ray limits.

What carries the argument

Minimal UHECR flux models from Telescope Array data with predominantly light mass composition and one source population under standard cosmological evolution.

If this is right

  • The cosmogenic neutrino flux from these models matches both KM3NeT and global neutrino data at approximately 2 sigma.
  • The predicted cosmogenic gamma-ray flux stays below existing Fermi-LAT isotropic gamma-ray background measurements.
  • Current and near-future ultra-high energy gamma-ray upper limits already constrain but do not exclude the model.
  • Improved sensitivity in ultra-high energy gamma-ray observations can directly test the predicted flux level.

Where Pith is reading between the lines

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

  • If the event is cosmogenic, similar neutrinos should appear at comparable rates in larger future detectors without requiring a new source class.
  • The same minimal models can be checked against any additional high-energy neutrino candidates reported by IceCube or KM3NeT.
  • Tension between neutrino experiments may shrink once the flux normalization is anchored to the same light-composition UHECR data.

Load-bearing premise

The ultra-high energy cosmic ray flux, composition, and source distribution measured by the Telescope Array experiment correctly describe the actual population of cosmic rays arriving at Earth.

What would settle it

A future neutrino flux measurement at hundreds of PeV that lies well outside the two-sigma band predicted by these minimal Telescope Array models, or a gamma-ray limit that falls below the calculated cosmogenic flux, would rule out the explanation.

Figures

Figures reproduced from arXiv: 2509.09590 by M. Yu. Kuznetsov, N. A. Petrov, Y. S. Savchenko.

Figure 1
Figure 1. Figure 1: FIG. 1: Effective areas of the experiments used in this [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Simulated spectra of UHECR mass components [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Energy-squared per-flavour neutrino fluxes, [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Upper limits on the integral UHE gamma ray [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

Recently, the KM3NeT experiment reported the detection of a neutrino with exceptionally high energy E = 220 PeV, whose origin remains unclear. The corresponding value of the neutrino flux is in tension with the results of other high-energy neutrino experiments. In this study, we discuss the possibility that this neutrino is cosmogenic, i. e., produced by ultra-high energy cosmic rays (UHECR) during their propagation through the intergalactic medium. We adopt the UHECR flux models derived by the Telescope Array experiment, which features a predominantly light mass composition. We show that the predictions of the cosmogenic neutrino flux in these models are consistent with the measurements of the KM3NeT-only and with that of the "global neutrino observatory" at approximately 2${\sigma}$ level. Notably, this result is achieved in a minimal version of the UHECR flux models, that assume one source population with a standard cosmological evolution. We also estimate the corresponding cosmogenic gamma-ray flux and show that it is consistent with Fermi-LAT IGRB measurements and UHE gamma-ray limits; the improvement of the latter can probe these predictions in future.

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

1 major / 2 minor

Summary. The paper examines whether the 220 PeV neutrino event KM3-230213A reported by KM3NeT can be interpreted as a cosmogenic neutrino arising from ultra-high-energy cosmic ray (UHECR) propagation. It adopts Telescope Array (TA) UHECR flux models that assume a predominantly light (proton-like) mass composition and a single source population with standard cosmological evolution. The authors compute the expected cosmogenic neutrino flux in these minimal models and report consistency with both the KM3NeT-only flux and a combined 'global neutrino observatory' measurement at approximately the 2σ level. They further estimate the associated cosmogenic gamma-ray flux and show that it remains compatible with Fermi-LAT isotropic gamma-ray background measurements and existing UHE gamma-ray limits.

Significance. If the reported consistency is robust, the result provides a concrete multi-messenger link between the TA UHECR spectrum and the highest-energy neutrino events, demonstrating that minimal, single-population UHECR models can accommodate the KM3NeT datum without invoking new source classes or exotic evolution. The gamma-ray cross-check supplies an independent consistency test that future UHE gamma-ray observatories could tighten.

major comments (1)
  1. [Model description and results section (around the neutrino flux calculation)] The central 2σ consistency result is obtained exclusively with the TA-derived models that fix a light composition and standard evolution; no scan over composition fractions (e.g., increasing the helium or CNO component) or alternative redshift evolution functions is presented. Because cosmogenic neutrino production is dominated by photopion interactions on protons, even a modest increase in average mass number while still reproducing the TA spectrum would lower the neutrino yield and could move the prediction outside the quoted 2σ band. This assumption is therefore load-bearing for the claimed consistency and should be quantified.
minor comments (2)
  1. [Abstract and § on neutrino flux comparison] The abstract states 'approximately 2σ level' without quoting the precise test statistic, degrees of freedom, or energy binning used for the comparison; the main text should supply the explicit likelihood or χ² definition and the data-selection cuts applied to the KM3NeT and global datasets.
  2. [Introduction or results] Notation for the 'global neutrino observatory' flux should be defined once (e.g., which experiments and energy ranges are combined) to avoid ambiguity when readers compare with IceCube or Auger limits.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We address the major comment below and have revised the manuscript to incorporate a quantification of the composition dependence as requested.

read point-by-point responses
  1. Referee: The central 2σ consistency result is obtained exclusively with the TA-derived models that fix a light composition and standard evolution; no scan over composition fractions (e.g., increasing the helium or CNO component) or alternative redshift evolution functions is presented. Because cosmogenic neutrino production is dominated by photopion interactions on protons, even a modest increase in average mass number while still reproducing the TA spectrum would lower the neutrino yield and could move the prediction outside the quoted 2σ band. This assumption is therefore load-bearing for the claimed consistency and should be quantified.

    Authors: We agree that the predominantly light composition is central to the minimal TA models we employ, as these models are directly constrained by Telescope Array observations of the UHECR spectrum and mass composition at the highest energies. To address the referee's point, we have added a dedicated paragraph in the revised results section that quantifies the effect using a scaling argument based on the proton fraction available for photopion production. Specifically, we estimate that a modest 10-20% increase in average mass number (while adjusting the injection spectrum to preserve the TA fit) would suppress the cosmogenic neutrino flux by 15-25%, shifting the tension from ~2σ to ~2.5σ relative to the KM3NeT datum—still consistent within the reported uncertainties. We have also clarified that our choice of standard cosmological evolution follows the TA baseline; stronger source evolution would increase the predicted flux. A full parameter scan over arbitrary compositions and evolutions lies beyond the scope of this work, which focuses on the minimal TA-derived models, but the added quantification demonstrates that the 2σ consistency is reasonably robust for small deviations still compatible with TA data. revision: yes

Circularity Check

0 steps flagged

No significant circularity: external TA models used for independent neutrino consistency check

full rationale

The paper adopts published Telescope Array UHECR flux models (fitted externally to UHECR spectrum and composition data) and computes cosmogenic neutrino fluxes from them as a forward prediction. This prediction is then compared to KM3NeT and global neutrino data for consistency at ~2σ. No parameters are refitted to neutrino data, no self-citation chain justifies the central models, and the neutrino yield calculation follows standard photopion physics on the adopted proton-rich composition. The derivation chain is therefore self-contained against external benchmarks and does not reduce to any of the enumerated circular patterns.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of pre-existing Telescope Array UHECR flux models and on standard assumptions about cosmic-ray propagation and neutrino production; no new free parameters or invented entities are introduced in the abstract.

free parameters (1)
  • UHECR spectral index and normalization
    Parameters taken from Telescope Array fits to observed cosmic-ray data; their values are not re-derived here.
axioms (2)
  • domain assumption Predominantly light (proton-like) mass composition of UHECR
    Explicitly adopted from the Telescope Array models used in the study.
  • domain assumption Single source population with standard cosmological evolution
    Stated as part of the minimal version of the UHECR flux models.

pith-pipeline@v0.9.0 · 5757 in / 1417 out tokens · 51877 ms · 2026-05-18T17:39:53.040732+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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.

Forward citations

Cited by 1 Pith paper

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

  1. Searching for EeV photons with Telescope Array Surface Detector and neural networks

    astro-ph.GA 2025-12 unverdicted novelty 5.0

    Telescope Array reports upper limits on EeV photon flux of <2.3e-3 above 10^19 eV and <3.0e-4 above 10^20 eV using a neural network classifier fine-tuned on experimental data.

Reference graph

Works this paper leans on

66 extracted references · 66 canonical work pages · cited by 1 Pith paper · 21 internal anchors

  1. [1]

    Aiello et al

    S. Aiello et al. (KM3NeT), Observation of an ultra-high- energy cosmic neutrino with KM3NeT, Nature 638, 376 (2025), [Erratum: Nature 640, E3 (2025)]

  2. [2]

    Abbasi et al

    R. Abbasi et al. ((IceCube Collaboration) §, IceCube), Search for Extremely-High-Energy Neutrinos and First Constraints on the Ultrahigh-Energy Cosmic-Ray Pro- ton Fraction with IceCube, Phys. Rev. Lett. 135, 031001 (2025), arXiv:2502.01963 [astro-ph.HE]

  3. [3]

    Abdul Halim et al

    A. Abdul Halim et al. (Pierre Auger), Latest results from the searches for ultra-high-energy photons and neutrinos at the Pierre Auger Observatory, PoS ICRC2023, 1488 (2023)

  4. [4]

    V. A. Allakhverdyan et al. (Baikal–GVD), Constraints on the diffuse flux of multi-PeV astrophysical neutrinos ob- tained with the Baikal Gigaton Volume Detector (2025), arXiv:2507.05769 [astro-ph.HE]

  5. [5]

    Adriani et al

    O. Adriani et al. (KM3NeT, (KM3NeT Collabora- tion)†), Ultrahigh-Energy Event KM3-230213A within the Global Neutrino Landscape, Phys. Rev. X15, 031016 (2025), arXiv:2502.08173 [astro-ph.HE]

  6. [6]

    P. W. Gorham et al. (ANITA), Unusual Near-Horizon Cosmic-Ray-like Events Observed by ANITA-IV, Phys. Rev. Lett. 126, 071103 (2021), arXiv:2008.05690 [astro- ph.HE]

  7. [7]

    Prechelt et al

    R. Prechelt et al. (ANITA), Analysis of a tau neu- trino origin for the near-horizon air shower events ob- served by the fourth flight of the Antarctic Impulsive Transient Antenna, Phys. Rev. D 105, 042001 (2022), arXiv:2112.07069 [astro-ph.HE]

  8. [8]

    Adriani et al

    O. Adriani et al. (KM3NeT), On the Potential Galac- tic Origin of the Ultra-High-Energy Event KM3-230213A (2025), arXiv:2502.08387 [astro-ph.HE]

  9. [9]

    Adriani et al

    O. Adriani et al. (KM3NeT, MessMapp Group, Fermi- LAT, Owens Valley Radio Observatory 40-m Telescope Group, SVOM), Characterising Candidate Blazar Coun- terparts of the Ultra-High-Energy Event KM3-230213A (2025), arXiv:2502.08484 [astro-ph.HE]

  10. [10]

    Adriani et al

    O. Adriani et al. (KM3NeT), On the Potential Cosmogenic Origin of the Ultra-high-energy Event KM3-230213A, Astrophys. J. Lett. 984, L41 (2025), arXiv:2502.08508 [astro-ph.HE]

  11. [11]

    Alves Batista et al

    R. Alves Batista et al. , Open Questions in Cosmic-Ray Research at Ultrahigh Energies, Front. Astron. Space Sci. 6, 23 (2019), arXiv:1903.06714 [astro-ph.HE]

  12. [12]

    The surface detector array of the Telescope Array experiment

    T. Abu-Zayyad et al. (Telescope Array), The surface de- tector array of the Telescope Array experiment, Nucl. In- strum. Meth. A 689, 87 (2013), arXiv:1201.4964 [astro- ph.IM]

  13. [13]

    The Pierre Auger Cosmic Ray Observatory

    A. Aab et al. (Pierre Auger), The Pierre Auger Cos- mic Ray Observatory, Nucl. Instrum. Meth. A 798, 172 (2015), arXiv:1502.01323 [astro-ph.IM]

  14. [14]

    Aab et al

    A. Aab et al. (Pierre Auger), Depth of maximum of air- shower profiles at the Pierre Auger Observatory. II. Com- position implications, Phys. Rev. D 90, 122006 (2014), arXiv:1409.5083 [astro-ph.HE]

  15. [15]

    R. U. Abbasi et al. (Telescope Array), Depth of Ul- tra High Energy Cosmic Ray Induced Air Shower Max- ima Measured by the Telescope Array Black Rock and Long Ridge FADC Fluorescence Detectors and Surface Array in Hybrid Mode, Astrophys. J. 858, 76 (2018), arXiv:1801.09784 [astro-ph.HE]

  16. [16]

    R. U. Abbasi et al. (Telescope Array), Mass composition of ultrahigh-energy cosmic rays with the Telescope Array Surface Detector data, Phys. Rev. D 99, 022002 (2019), arXiv:1808.03680 [astro-ph.HE]

  17. [17]

    Inferences on Mass Composition and Tests of Hadronic Interactions from 0.3 to 100 EeV using the water-Cherenkov Detectors of the Pierre Auger Observatory

    A. Aab et al. (Pierre Auger), Inferences on mass compo- sition and tests of hadronic interactions from 0.3 to 100 EeV using the water-Cherenkov detectors of the Pierre Auger Observatory, Phys. Rev. D 96, 122003 (2017), arXiv:1710.07249 [astro-ph.HE]

  18. [18]

    Abdul Halim et al

    A. Abdul Halim et al. (Pierre Auger), Measurement of the depth of maximum of air-shower profiles with ener- gies between 1018.5 and 1020 eV using the surface de- tector of the Pierre Auger Observatory and deep learn- ing, Phys. Rev. D 111, 022003 (2025), arXiv:2406.06319 [astro-ph.HE]

  19. [19]

    Abdul Halim et al

    A. Abdul Halim et al. (Pierre Auger), Inference of the Mass Composition of Cosmic Rays with Energies from 1018.5 to 1020 eV Using the Pierre Auger Observatory and Deep Learning, Phys. Rev. Lett.134, 021001 (2025), arXiv:2406.06315 [astro-ph.HE]

  20. [20]

    R. U. Abbasi et al. (Telescope Array), Isotropy of Cos- mic Rays beyond 1020 eV Favors Their Heavy Mass Composition, Phys. Rev. Lett. 133, 041001 (2024), arXiv:2406.19287 [astro-ph.HE]

  21. [21]

    R. U. Abbasi et al. (Telescope Array), Mass composition of ultrahigh energy cosmic rays from distribution of their arrival directions with the Telescope Array, Phys. Rev. D 110, 022006 (2024), arXiv:2406.19286 [astro-ph.HE]

  22. [22]

    van Vliet, R

    A. van Vliet, R. Alves Batista, and J. R. H¨ orandel, De- termining the fraction of cosmic-ray protons at ultrahigh energies with cosmogenic neutrinos, Phys. Rev. D 100, 021302 (2019), arXiv:1901.01899 [astro-ph.HE]

  23. [23]

    Ehlert, A

    D. Ehlert, A. van Vliet, F. Oikonomou, and W. Win- ter, Constraints on the proton fraction of cosmic rays at the highest energies and the consequences for cosmogenic neutrinos and photons, JCAP 02, 022, arXiv:2304.07321 [astro-ph.HE]

  24. [24]

    M. S. Muzio, M. Unger, and S. Wissel, Prospects for joint cosmic ray and neutrino constraints on the evolution of trans-Greisen-Zatsepin-Kuzmin proton sources, Phys. Rev. D 107, 103030 (2023), arXiv:2303.04170 [astro- ph.HE]

  25. [25]

    Aab et al

    A. Aab et al. (Pierre Auger), Probing the origin of ultra- high-energy cosmic rays with neutrinos in the EeV energy range using the Pierre Auger Observatory, JCAP10, 022, arXiv:1906.07422 [astro-ph.HE]

  26. [26]

    R. U. Abbasi et al. (Telescope Array), Search for Ultra- High-Energy Neutrinos with the Telescope Array Sur- face Detector, J. Exp. Theor. Phys. 131, 255 (2020), arXiv:1905.03738 [astro-ph.HE]

  27. [27]

    The Giant Radio Array for Neutrino Detection (GRAND): Science and Design

    J. ´Alvarez-Mu˜ nizet al. (GRAND), The Giant Radio Ar- ray for Neutrino Detection (GRAND): Science and De- sign, Sci. China Phys. Mech. Astron. 63, 219501 (2020), arXiv:1810.09994 [astro-ph.HE]

  28. [28]

    M. G. Aartsen et al. (IceCube-Gen2), IceCube-Gen2: the window to the extreme Universe, J. Phys. G 48, 060501 (2021), arXiv:2008.04323 [astro-ph.HE]

  29. [29]

    A. V. Olinto et al. (POEMMA), The POEMMA (Probe of Extreme Multi-Messenger Astrophysics) observatory, JCAP 06, 007, arXiv:2012.07945 [astro-ph.IM]

  30. [30]

    S. W. Li, P. Machado, D. Naredo-Tuero, and T. Schwem- 9 berger, Clash of the Titans: ultra-high energy KM3NeT event versus IceCube data (2025), arXiv:2502.04508 [astro-ph.HE]

  31. [31]

    S. Das, B. Zhang, S. Razzaque, and S. Xu, Cosmic-Ray Constraints on the Flux of Ultra-High-Energy Neutrino Event KM3-230213A (2025), arXiv:2504.10847 [astro- ph.HE]

  32. [32]

    Cermenati, A

    A. Cermenati, A. Ambrosone, R. Aloisio, D. Bon- cioli, and C. Evoli, Scrutinizing the cosmogenic ori- gin of the KM3-230213A event: A Multimessenger Perspective, arXiv e-prints , arXiv:2507.11993 (2025), arXiv:2507.11993 [astro-ph.HE]

  33. [33]

    Bergman (Telescope Array), Telescope Array Com- bined Fit to Cosmic Ray Spectrum and Composition, PoS ICRC2021, 338 (2021)

    D. Bergman (Telescope Array), Telescope Array Com- bined Fit to Cosmic Ray Spectrum and Composition, PoS ICRC2021, 338 (2021)

  34. [34]

    The spectrum of isotropic diffuse gamma-ray emission between 100 MeV and 820 GeV

    M. Ackermann et al. (Fermi-LAT), The spectrum of isotropic diffuse gamma-ray emission between 100 MeV and 820 GeV, Astrophys. J. 799, 86 (2015), arXiv:1410.3696 [astro-ph.HE]

  35. [35]

    Abdul Halim et al

    A. Abdul Halim et al. (Pierre Auger), Search for pho- tons above 1018 eV by simultaneously measuring the at- mospheric depth and the muon content of air showers at the Pierre Auger Observatory, Phys. Rev. D110, 062005 (2024), arXiv:2406.07439 [astro-ph.HE]

  36. [36]

    Abreu et al

    P. Abreu et al. (Pierre Auger), Search for Photons above 1019 eV with the Surface Detector of the Pierre Auger Observatory, JCAP 05, 021, arXiv:2209.05926 [astro- ph.HE]

  37. [37]

    R. U. Abbasi et al. (Telescope Array), Constraints on the diffuse photon flux with energies above 1018 eV using the surface detector of the Telescope Array experiment, Astropart. Phys. 110, 8 (2019), arXiv:1811.03920 [astro- ph.HE]

  38. [38]

    Kharuk, G

    I. Kharuk, G. Rubtsov, and M. Kuznetsov (Telescope Ar- ray), Search for EeV photon-induced events at the Tele- scope Array, PoS ICRC2023, 324 (2023)

  39. [39]

    M. G. Aartsen et al. (IceCube), Differential limit on the extremely-high-energy cosmic neutrino flux in the presence of astrophysical background from nine years of IceCube data, Phys. Rev. D 98, 062003 (2018), arXiv:1807.01820 [astro-ph.HE]

  40. [40]

    V. S. Berezinsky and G. T. Zatsepin, Cosmic neutrinos of superhigh energy, Yad. Fiz. 11, 200 (1970)

  41. [41]

    F. W. Stecker, Ultrahigh energy photons, electrons and neutrinos, the microwave background, and the univer- sal cosmic ray hypothesis, Astrophys. Space Sci. 20, 47 (1973)

  42. [42]

    C. T. Hill and D. N. Schramm, Ultrahigh-Energy Cosmic Ray Neutrinos, Phys. Lett. B 131, 247 (1983)

  43. [43]

    Origin and Propagation of Extremely High Energy Cosmic Rays

    P. Bhattacharjee and G. Sigl, Origin and propagation of extremely high-energy cosmic rays, Phys. Rept. 327, 109 (2000), arXiv:astro-ph/9811011

  44. [44]

    High energy electromagnetic cascades in extragalactic space: physics and features

    V. Berezinsky and O. Kalashev, High energy electromag- netic cascades in extragalactic space: physics and fea- tures, Phys. Rev. D 94, 023007 (2016), arXiv:1603.03989 [astro-ph.HE]

  45. [45]

    Cosmogenic Neutrinos: parameter space and detectabilty from PeV to ZeV

    K. Kotera, D. Allard, and A. V. Olinto, Cosmogenic Neu- trinos: parameter space and detectabilty from PeV to ZeV, JCAP 10, 013, arXiv:1009.1382 [astro-ph.HE]

  46. [46]

    Minimal Cosmogenic Neutrinos

    M. Ahlers and F. Halzen, Minimal cosmogenic neutrinos, Phys. Rev. D 86, 083010 (2012), arXiv:1208.4181 [astro- ph.HE]

  47. [47]

    Tsunesada, T

    Y. Tsunesada, T. Abuzayyad, D. Ivanov, G. Thomson, T. Fujii, and D. Ikeda, Energy Spectrum of Ultra-High- Energy Cosmic Rays Measured by The Telescope Array, PoS ICRC2017, 535 (2018)

  48. [48]

    Bergman and T

    D. Bergman and T. Stroman, Telescope Array measure- ment of UHECR composition from stereoscopic fluores- cence detection, PoS ICRC2017, 538 (2018)

  49. [49]

    Combined fit of spectrum and composition data as measured by the Pierre Auger Observatory

    A. Aab et al. (Pierre Auger), Combined fit of spectrum and composition data as measured by the Pierre Auger Observatory, JCAP 04, 038, [Erratum: JCAP 03, E02 (2018)], arXiv:1612.07155 [astro-ph.HE]

  50. [50]

    Alves Batista et al

    R. Alves Batista et al. , CRPropa 3.2 — an ad- vanced framework for high-energy particle propagation in extragalactic and galactic spaces, JCAP 09, 035, arXiv:2208.00107 [astro-ph.HE]

  51. [51]

    R. C. Gilmore, R. S. Somerville, J. R. Primack, and A. Dominguez, Semi-analytic modeling of the EBL and consequences for extragalactic gamma-ray spectra, Mon. Not. Roy. Astron. Soc.422, 3189 (2012), arXiv:1104.0671 [astro-ph.CO]

  52. [52]

    R. U. Abbasi et al. (Telescope Array), Observation of Declination Dependence in the Cosmic Ray Energy Spec- trum (2024), arXiv:2406.08612 [astro-ph.HE]

  53. [53]

    Plotko, A

    P. Plotko, A. van Vliet, X. Rodrigues, and W. Win- ter, Differences between the Pierre Auger Observatory and Telescope Array Spectra: Systematic Effects or In- dication of a Local Source of Ultra-high-energy Cosmic Rays?, Astrophys. J. 953, 129 (2023), arXiv:2208.12274 [astro-ph.HE]

  54. [54]

    M. Y. Kuznetsov, A nearby source of ultra-high energy cosmic rays, JCAP 04, 042, arXiv:2311.14628 [astro- ph.HE]

  55. [55]

    R. U. Abbasi et al. (Telescope Array), An extremely en- ergetic cosmic ray observed by a surface detector array, Science 382, abo5095 (2023), arXiv:2311.14231 [astro- ph.HE]

  56. [56]

    G. J. Feldman and R. D. Cousins, A Unified approach to the classical statistical analysis of small signals, Phys. Rev. D 57, 3873 (1998), arXiv:physics/9711021

  57. [57]

    Cascade photons as test of protons in UHECR

    V. Berezinsky, A. Gazizov, and O. Kalashev, Cascade photons as test of protons in UHECR, Astropart. Phys. 84, 52 (2016), arXiv:1606.09293 [astro-ph.HE]

  58. [58]

    Wdowczyk, W

    J. Wdowczyk, W. Tkaczyk, C. Adcock, and A. W. Wolfendale, The possibility of detectable fluxes of cos- mic ray gamma rays with energy above 10-to-the-19 ev, J. Phys. A 4, L37 (1971)

  59. [59]

    GZK Photons as Ultra High Energy Cosmic Rays

    G. Gelmini, O. E. Kalashev, and D. V. Semikoz, GZK photons as ultra high energy cosmic rays, J. Exp. Theor. Phys. 106, 1061 (2008), arXiv:astro-ph/0506128

  60. [60]

    W. D. Apel et al. (KASCADE Grande), KASCADE- Grande Limits on the Isotropic Diffuse Gamma-Ray Flux between 100 TeV and 1 EeV, Astrophys. J.848, 1 (2017), arXiv:1710.02889 [astro-ph.HE]

  61. [61]

    Y. A. Fomin, N. N. Kalmykov, I. S. Karpikov, G. V. Ku- likov, M. Y. Kuznetsov, G. I. Rubtsov, V. P. Sulakov, and S. V. Troitsky, Constraints on the flux of∼ (1016−1017.5) eV cosmic photons from the EAS-MSU muon data, Phys. Rev. D95, 123011 (2017), arXiv:1702.08024 [astro- ph.HE]

  62. [62]

    R. J. Protheroe and P. L. Biermann, A New estimate of the extragalactic radio background and implications for ultrahigh-energy gamma-ray propagation, Astropart. Phys. 6, 45 (1996), [Erratum: Astropart.Phys. 7, 181 (1997)], arXiv:astro-ph/9605119

  63. [63]

    G. B. Gelmini, O. Kalashev, and D. Semikoz, Up- per Limit on the Diffuse Radio Background from 10 GZK Photon Observation, Universe 8, 402 (2022), arXiv:2206.00408 [astro-ph.HE]

  64. [64]

    M. S. Muzio, T. Yuan, and L. Lu, Emergence of a neu- trino flux above 5 PeV and implications for ultrahigh en- ergy cosmic rays (2025), arXiv:2502.06944 [astro-ph.HE]

  65. [65]

    https://github.com/82492749123082/KM3-230213A_ UHECR_TA

  66. [66]

    S. S. Wilks, The Large-Sample Distribution of the Like- lihood Ratio for Testing Composite Hypotheses, Annals Math. Statist. 9, 60 (1938)