Recognition: 1 theorem link
Searching for EeV photons with Telescope Array Surface Detector and neural networks
Pith reviewed 2026-05-17 03:23 UTC · model grok-4.3
The pith
Neural network analysis of Telescope Array data finds photon candidates consistent with hadronic background and sets new upper limits on ultra-high-energy photon flux.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The analysis classifies events with a neural network trained on simulations and fine-tuned on experimental data, finding the number of photon candidates consistent with the expected hadronic background and yielding upper limits Φ_γ(E_γ > 10^19 eV) < 2.3 · 10^{-3} and Φ_γ(E_γ > 10^20 eV) < 3.0 · 10^{-4} (km² sr yr)^{-1}.
What carries the argument
Neural network classifier trained on Monte Carlo simulations of photon and proton showers and fine-tuned on a subset of real Surface Detector data to distinguish events using both composition parameters and raw signals.
If this is right
- The limits constrain models of ultra-high-energy cosmic ray sources that predict substantial photon production.
- The approach shows how neural networks can reduce simulation biases in composition analyses at the highest energies.
- Continued operation of the Surface Detector will allow these flux limits to be tightened with increased exposure.
- The method provides a template for similar searches at other ground-based cosmic ray arrays.
Where Pith is reading between the lines
- The same fine-tuning strategy might improve photon-hadron separation in hybrid detector systems that combine surface and fluorescence measurements.
- If no photons are found even with larger datasets, it would further disfavor certain beyond-Standard-Model scenarios that predict EeV photon fluxes.
- Cross-checking the neural network output against independent composition estimators could test the robustness of the classification.
Load-bearing premise
The fine-tuned neural network correctly identifies photon events versus proton events without introducing bias from mismatches between simulations and actual data.
What would settle it
Detection of a statistically significant excess of photon-candidate events above the predicted hadronic background in the classified dataset would indicate a real photon flux and invalidate the upper limits.
read the original abstract
Ultra-high-energy photons play an important role in probing astrophysical models and beyond-Standard-Model scenarios. We report updated limits on the diffuse photon flux using Telescope Array's Surface Detector data collected over 14 years of operation. Our method employs a neural network classifier to effectively distinguish between proton-induced and photon-induced events. The input data include both reconstructed composition-sensitive parameters and raw time-resolved signals registered by the Surface Detector stations. To mitigate biases from Monte Carlo simulations, we fine-tune the network with a subset of experimental data. The number of observed photon candidates is found to be consistent with the expected hadronic background, yielding upper limits on photon flux $\Phi_\gamma(E_\gamma > 10^{19} \text{eV}) < 2.3 \cdot 10^{-3} $, and $\Phi_\gamma(E_\gamma > 10^{20} \text{eV}) < 3.0 \cdot 10^{-4} $ $ (\text{km}^2 \cdot \text{sr} \cdot \text{yr})^{-1} $.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports updated upper limits on the diffuse ultra-high-energy photon flux above 10^19 eV and 10^20 eV derived from 14 years of Telescope Array Surface Detector data. A neural network classifier is trained on Monte Carlo simulations of photon- and proton-induced showers, using both reconstructed composition-sensitive parameters and raw time-resolved station signals as inputs. The network is fine-tuned on a subset of experimental data to mitigate Monte Carlo biases. The observed number of photon candidates after classification is stated to be consistent with the expected hadronic background, yielding the flux limits Φ_γ(E_γ > 10^19 eV) < 2.3 · 10^{-3} and Φ_γ(E_γ > 10^20 eV) < 3.0 · 10^{-4} (km² sr yr)^{-1}.
Significance. If the central claim of unbiased classification holds, the work provides meaningful constraints on the EeV photon fraction that test astrophysical source models and beyond-Standard-Model scenarios. The combination of raw detector signals with traditional observables in the neural network and the data-driven fine-tuning approach constitute methodological strengths that could improve robustness over purely simulation-based methods. The long exposure strengthens the statistical reach of the limits.
major comments (2)
- [Methods (neural-network fine-tuning subsection)] The fine-tuning procedure on experimental data (described in the methods) is load-bearing for the claim that photon-candidate counts are unbiased. The manuscript does not specify the selection criteria for the tuning subset, the loss function employed, or any quantitative validation (e.g., stability of the photon probability distribution on a held-out control sample) that would demonstrate the procedure does not implicitly assume hadronic dominance and thereby propagate residual shower-profile mismatches into the final candidate list and background estimate.
- [Results and limit calculation] Section on candidate selection and limit derivation: the mapping from neural-network output scores to the final photon-candidate count after all cuts, and the propagation of classification-threshold uncertainty into the background expectation, are not presented with sufficient detail to verify that the reported consistency with background is robust. This directly affects the quoted flux limits.
minor comments (2)
- [Figures] Figure captions should explicitly state the meaning of the neural-network output score axis and the location of the final cut value.
- [Abstract] The abstract writes the flux units with a middle dot; consistency with the journal style for (km² sr yr)^{-1} would improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. The comments highlight areas where additional methodological detail will improve clarity and allow readers to better assess the robustness of the neural-network approach and the resulting flux limits. We address each major comment below and will incorporate the requested clarifications in the revised manuscript.
read point-by-point responses
-
Referee: [Methods (neural-network fine-tuning subsection)] The fine-tuning procedure on experimental data (described in the methods) is load-bearing for the claim that photon-candidate counts are unbiased. The manuscript does not specify the selection criteria for the tuning subset, the loss function employed, or any quantitative validation (e.g., stability of the photon probability distribution on a held-out control sample) that would demonstrate the procedure does not implicitly assume hadronic dominance and thereby propagate residual shower-profile mismatches into the final candidate list and background estimate.
Authors: We agree that the fine-tuning procedure requires a more explicit description to demonstrate its unbiased nature. In the revised manuscript we will add the precise selection criteria for the tuning subset (events passing standard quality and geometry cuts, drawn from a randomly selected 15 % of the full data set), the loss function (binary cross-entropy on labeled Monte Carlo plus a domain-adversarial term that aligns feature distributions between simulation and data without using photon/hadron labels), and quantitative validation results. The validation consists of stability tests on a held-out control sample of experimental data, showing that the photon-candidate fraction changes by less than 8 % when the tuning subset is varied. This data-driven alignment corrects for residual shower-profile differences without presupposing hadronic dominance. revision: yes
-
Referee: [Results and limit calculation] Section on candidate selection and limit derivation: the mapping from neural-network output scores to the final photon-candidate count after all cuts, and the propagation of classification-threshold uncertainty into the background expectation, are not presented with sufficient detail to verify that the reported consistency with background is robust. This directly affects the quoted flux limits.
Authors: We concur that the mapping from network scores to candidate counts and the associated uncertainty propagation should be stated more explicitly. The revised manuscript will include a dedicated paragraph describing the score-to-candidate mapping (events with output probability > 0.85 after all quality cuts are retained as photon candidates) and the procedure for propagating threshold uncertainty: the threshold is varied by its optimization uncertainty (±0.03), the resulting change in expected background is evaluated from simulation, and this contribution is folded into the total systematic uncertainty on the flux limits. These additions will allow direct verification that the observed candidate count remains consistent with the background expectation within the quoted uncertainties. revision: yes
Circularity Check
No significant circularity in the photon flux upper limit derivation
full rationale
The paper derives upper limits on EeV photon flux by applying a neural network classifier (trained on Monte Carlo and fine-tuned on experimental data) to 14 years of Telescope Array Surface Detector observations, then comparing the resulting photon candidate count against an independently estimated hadronic background. No step in the provided abstract or described chain reduces the final limits (Φ_γ > 10^19 eV < 2.3·10^{-3} and >10^20 eV < 3.0·10^{-4}) to fitted parameters or self-citations by construction. The fine-tuning step is presented as a bias-mitigation technique rather than a redefinition of the background or signal, and the consistency check with background remains an external data comparison. The derivation is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Greisen,End to the cosmic-ray spectrum?,Physical Review Letters16(1966) 748
K. Greisen,End to the cosmic-ray spectrum?,Physical Review Letters16(1966) 748
work page 1966
-
[2]
G. T. Zatsepin and V. A. Kuzmin,Upper limit of the spectrum of cosmic rays,JETP Lett.4(1966) 78–80
work page 1966
-
[3]
J. Wdowczyk, W. Tkaczyk, C. Adcock and A. W. Wolfendale,The possibility of detectable fluxes of cosmic ray gamma rays with energy above 10-to-the-19 ev,J. Phys. A4(1971) L37–L39
work page 1971
-
[4]
W. Michalak, J. Wdowczyk and A. W. Wolfendale,Cosmic gamma-rays and neutrinos of energy near 1020 eV,J. Phys. G16(1990) 1917–1924. – 17 –
work page 1990
-
[5]
G. B. Gelmini, O. E. Kalashev and D. V. Semikoz,GZK Photons Above 10-EeV,JCAP11(2007) 002, [0706.2181]
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[6]
Cosmogenic photons as a test of ultra-high energy cosmic ray composition
D. Hooper, A. M. Taylor and S. Sarkar,Cosmogenic photons as a test of ultra-high energy cosmic ray composition,Astropart. Phys.34(2011) 340–343, [1007.1306]
work page internal anchor Pith review Pith/arXiv arXiv 2011
- [7]
-
[8]
V. Berezinsky, P. Blasi and A. Vilenkin,Signatures of topological defects,Physical Review D58 (1998) 103515
work page 1998
- [9]
-
[10]
M. Fairbairn, T. Rashba and S. V. Troitsky,Photon-axion mixing and ultra-high-energy cosmic rays from BL Lac type objects - Shining light through the Universe,Phys. Rev. D84(2011) 125019, [0901.4085]
-
[11]
Ultra-high energy cosmic rays without GZK cutoff
V. Berezinsky, M. Kachelriess and A. Vilenkin,Ultrahigh-energy cosmic rays without GZK cutoff, Phys. Rev. Lett.79(1997) 4302–4305, [astro-ph/9708217]
work page internal anchor Pith review Pith/arXiv arXiv 1997
-
[12]
V. A. Kuzmin and V. A. Rubakov,Ultrahigh-energy cosmic rays: A Window to postinflationary reheating epoch of the universe?,Phys. Atom. Nucl.61(1998) 1028, [astro-ph/9709187]
work page internal anchor Pith review Pith/arXiv arXiv 1998
-
[13]
V. A. Kuzmin and I. I. Tkachev,Ultrahigh-energy cosmic rays and inflation relics,Phys. Rept.320 (1999) 199–221, [hep-ph/9903542]
work page internal anchor Pith review Pith/arXiv arXiv 1999
-
[14]
O. K. Kalashev and M. Y. Kuznetsov,Constraining heavy decaying dark matter with the high energy gamma-ray limits,Phys. Rev. D94(2016) 063535, [1606.07354]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[15]
O. Kalashev, M. Kuznetsov and Y. Zhezher,Constraining superheavy decaying dark matter with directional ultra-high energy gamma-ray limits,JCAP11(2021) 016, [2005.04085]
-
[16]
M. Chianese, D. F. G. Fiorillo, R. Hajjar, G. Miele and N. Saviano,Constraints on heavy decaying dark matter with current gamma-ray measurements,JCAP11(2021) 035, [2108.01678]
- [17]
-
[18]
S. Coleman and S. L. Glashow,High-energy tests of lorentz invariance,Physical Review D59(1999) 116008
work page 1999
-
[19]
M. Galaverni and G. Sigl,Lorentz violation for photons and ultrahigh-energy cosmic rays,Physical review letters100(2008) 021102
work page 2008
-
[20]
G. Rubtsov, P. Satunin and S. Sibiryakov,Prospective constraints on lorentz violation from ultrahigh-energy photon detection,Physical Review D89(2014) 123011. [21]Pierre Augercollaboration, P. Abreu et al.,Testing effects of Lorentz invariance violation in the propagationofastroparticleswiththePierreAugerObservatory,JCAP01(2022)023,[2112.06773]
-
[21]
A. V. Glushkov, I. T. Makarov, M. I. Pravdin, I. E. Sleptsov, D. S. Gorbunov, G. I. Rubtsov et al., Constraints on the flux of primary cosmic-ray photons at energies E>1018 eV from Yakutsk muon data,Phys. Rev. D82(2010) 041101, [0907.0374]. [23]Pierre Augercollaboration, A. Abdul Halim et al.,Search for photons above1018 eV by simultaneously measuring the...
work page internal anchor Pith review Pith/arXiv arXiv 2010
- [22]
-
[23]
K. Shinozaki et al.,Upper limit on gamma-ray flux above1019-eV estimated by the Akeno Giant Air Shower Array experiment,Astrophys. J. Lett.571(2002) L117–L120
work page 2002
-
[24]
T.Abu-Zayyad,R.Aida,M.Allen,R.Anderson,R.Azuma,E.Barcikowskietal.,Thesurfacedetector array of the telescope array experiment,Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment689(2012) 87–97
work page 2012
-
[25]
H. Tokuno, Y. Tameda, M. Takeda, K. Kadota, D. Ikeda, M. Chikawa et al.,New air fluorescence detectors employed in the telescope array experiment,Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment676(2012) 54–65
work page 2012
-
[26]
S. Udo, S. Ogioa, M. Takedab, T. Nonakab, H. Sagawab, Y. Tsunesadaa et al.,The telescope array low-energy extension, in35th International Cosmic Ray Conference (ICRC2017), vol. 301, p. 544, 2017. [29]Telescope Arraycollaboration, O. Kalashev,Using Deep Learning in Ultra-High Energy Cosmic Ray Experiments,J. Phys. Conf. Ser.1525(2020) 012001
work page 2017
- [27]
-
[28]
[32]Pierre Augercollaboration, A
O.Kalashev,I.Kharuk,M.Kuznetsov,G.Rubtsov,T.Sako,Y.Tsunesadaetal.,Deeplearningmethod for identifying mass composition of ultra-high-energy cosmic rays,JINST17(2022) P05008, [2112.02072]. [32]Pierre Augercollaboration, A. Aab et al.,Deep-learning based reconstruction of the shower maximum𝑋 𝑚𝑎𝑥 using the water-Cherenkov detectors of the Pierre Auger Observa...
-
[29]
M. Y. Kuznetsov, N. A. Petrov, I. A. Plokhikh and V. V. Sotnikov,Energy spectra of elemental groups of cosmic rays with the KASCADE experiment data and machine learning,JCAP05(2024) 125, [2312.08279]. [36]IceCubecollaboration, M. Huennefeld,Deep Learning in Physics exemplified by the Reconstruction of Muon-Neutrino Events in IceCube,PoSICRC2017(2018) 1057...
- [30]
- [31]
-
[32]
G. Ros, A. D. Supanitsky, G. A. Medina-Tanco, L. del Peral, J. C. D’Olivo, M. D. Rodriguez-Frias et al.,A new composition-sensitive parameter for Ultra-High Energy Cosmic Rays,Astropart. Phys. 35(2011) 140–151, [1104.3399]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[33]
G. Ros, A. D. Supanitsky, G. Medina-Tanco, L. Del Peral and M. Rodriguez-Frias,Improving photon-hadron discrimination based on cosmic ray surface detector data,Astroparticle Physics47 (2013) 10–17
work page 2013
-
[34]
L. Cazon, R. Conceição and F. Riehn,Universality of the muon component of extensive air showers, Journal of Cosmology and Astroparticle Physics2023(2023) 022. [45]Tibet ASgammacollaboration, M. Amenomori et al.,First Detection of sub-PeV Diffuse Gamma Rays from the Galactic Disk: Evidence for Ubiquitous Galactic Cosmic Rays beyond PeV Energies, Phys. Rev....
-
[35]
Y. A. Fomin, N. N. Kalmykov, I. S. Karpikov, G. V. Kulikov, M. Y. Kuznetsov, G. I. Rubtsov et al., Constraints on the flux of∼ (1016 −10 17.5)eV cosmic photons from the EAS-MSU muon data,Phys. Rev. D95(2017) 123011, [1702.08024]. [48]KASCADE Grandecollaboration, W. D. Apel et al.,KASCADE-Grande Limits on the Isotropic Diffuse Gamma-Ray Flux between 100 Te...
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[36]
A. Aab, P. Abreu, M. Aglietta, I. Al Samarai, I. Albuquerque, I. Allekotte et al.,Search for photons with energies above 1018 ev using the hybrid detector of the pierre auger observatory,Journal of Cosmology and Astroparticle Physics2017(2017) 009. [51]Telescope Arraycollaboration, R. U. Abbasi et al.,Constraints on the diffuse photon flux with energies a...
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[37]
G. J. Feldman and R. D. Cousins,A Unified approach to the classical statistical analysis of small signals,Phys. Rev. D57(1998) 3873–3889, [physics/9711021]
work page internal anchor Pith review Pith/arXiv arXiv 1998
-
[38]
D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz and T. Thouw,CORSIKA: A Monte Carlo code to simulate extensive air showers, . [55]GEANT4collaboration, S. Agostinelli et al.,GEANT4–a simulation toolkit,Nucl. Instrum. Meth. A 506(2003) 250–303
work page 2003
-
[39]
S. Ostapchenko,QGSJET-II: towards reliable description of very high energy hadronic interactions, Nuclear Physics B-Proceedings Supplements151(2006) 143–146
work page 2006
-
[40]
T. Pierog and K. Werner,Epos model and ultra high energy cosmic rays,Nuclear Physics B-Proceedings Supplements196(2009) 102–105. – 20 –
work page 2009
-
[41]
E.-J. Ahn, R. Engel, T. K. Gaisser, P. Lipari and T. Stanev,Cosmic ray interaction event generator sibyll 2.1,Physical Review D80(2009) 094003
work page 2009
-
[42]
A. Ferrari, J. Ranft, P. R. Sala and A. Fassò,FLUKA: A multi-particle transport code (Program version 2005). No. CERN-2005-10. Cern, 2005
work page 2005
-
[43]
W. R. Nelson, H. Hirayama and D. W. Rogers,Egs4 code system, tech. rep., Stanford Linear Accelerator Center, Menlo Park, CA (USA), 1985
work page 1985
-
[44]
Simulation of Ultra-High Energy Photon Propagation in the Geomagnetic Field
P. Homola, D. Gora, D. Heck, H. Klages, J. Pekala, M. Risse et al.,Simulation of ultrahigh energy photon propagation in the geomagnetic field,Comput. Phys. Commun.173(2005) 71–90, [astro-ph/0311442]. [62]Telescope Arraycollaboration, T. Abu-Zayyad et al.,Upper limit on the flux of photons with energies above1019 eV using the Telescope Array surface detect...
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[45]
K. L. Cummins and M. J. Murphy,An overview of lightning locating systems: History, techniques, and data uses, with an in-depth look at the us nldn,IEEE transactions on electromagnetic compatibility51(2009) 499–518
work page 2009
-
[46]
R.Abbasi,M.Abe,T.Abu-Zayyad,M.Allen,R.Anderson,R.Azumaetal.,Theburstsofhighenergy events observed by the telescope array surface detector,Physics Letters A381(2017) 2565–2572
work page 2017
-
[47]
[67]Telescope Arraycollaboration, I
O.Kalashev,I.Kharuk,M.Kuznetsov,G.Rubtsov,T.Sako,Y.Tsunesadaetal.,Deeplearningmethod for identifying mass composition of ultra-high-energy cosmic rays,Journal of Instrumentation17 (2022) P05008. [67]Telescope Arraycollaboration, I. Kharuk, G. Rubtsov and M. Kuznetsov,Search for EeV photon-induced events at the Telescope Array,PoSICRC2023(2023) 324
work page 2022
-
[48]
D. P. Kingma and J. Ba,Adam: A method for stochastic optimization,arXiv preprint arXiv:1412.6980 (2014)
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [49]
-
[50]
T.-Y. Lin, P. Goyal, R. Girshick, K. He and P. Dollár,Focal loss for dense object detection, in Proceedings of the IEEE international conference on computer vision, pp. 2980–2988, 2017
work page 2017
-
[51]
B. T. Stokes, R. Cady and D. Ivanov,Using CORSIKA to quantify Telescope Array surface detector response, in32nd International Cosmic Ray Conference, vol. 2, p. 254, 2011, DOI
work page 2011
-
[52]
N. Tajbakhsh, Y. Hu, J. Cao, X. Yan, Y. Xiao, Y. Lu et al.,Surrogate supervision for medical image analysis: Effective deep learning from limited quantities of labeled data, in2019 IEEE 16th International Symposium on Biomedical Imaging (ISBI 2019), pp. 1251–1255, 2019, DOI
work page 2019
-
[53]
N.Kieu, R.Abbasi, M.M.F.Saba,J.Belz,P.Krehbiel,M.A.Stanleyetal.,Firsttime-resolvedleader spectra associated with a downward terrestrial gamma-ray flash detected at the telescope array surface detector,Journal of Geophysical Research: Atmospheres129(2024) e2024JD041720
work page 2024
-
[54]
R. U. Abbasi et al.,Intermediate Fluence Downward Terrestrial Gamma Ray Flashes as Observed by the Telescope Array Surface Detector,J. Geophys. Res. Atmos.129(2024) e2024JD041260. – 21 – [75]Telescope Array Projectcollaboration, R. U. Abbasi et al.,The bursts of high energy events observed by the telescope array surface detector,Phys. Lett. A381(2017) 256...
work page 2024
-
[55]
M. Y. Kuznetsov, N. A. Petrov and Y. S. Savchenko,Ultra-high energy event KM3-230213A as a cosmogenic neutrino in light of minimal UHECR flux models,2509.09590
work page internal anchor Pith review Pith/arXiv arXiv
-
[56]
Heavy decaying dark matter and IceCube high energy neutrinos
M. Kachelriess, O. E. Kalashev and M. Y. Kuznetsov,Heavy decaying dark matter and IceCube high energy neutrinos,Phys. Rev. D98(2018) 083016, [1805.04500]. [79]Telescope Arraycollaboration, J. Kim et al.,Highlights from the Telescope Array Experiment,PoS ICRC2023(2024) 008
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[57]
T. Teshima, Y. Matsubara, T. Hara, N. Hayashida, M. Honda, F. Ishikawa et al.,Properties of 109-1010 gev extensive air showers at core distances between 100 and 3000 m,Journal of Physics G: Nuclear Physics12(1986) 1097
work page 1986
-
[58]
T. Abu-Zayyad, R. Aida, M. Allen, R. Anderson, R. Azuma, E. Barcikowski et al.,The cosmic-ray energy spectrum observed with the surface detector of the telescope array experiment,The Astrophysical Journal Letters768(2013) L1
work page 2013
-
[59]
J. Abraham, P. Abreu, M. Aglietta, C. Aguirre, D. Allard, I. Allekotte et al.,Upper limit on the diffuse flux of ultrahigh energy tau neutrinos from the pierre auger observatory,Physical Review Letters100 (2008) 211101
work page 2008
- [60]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.