Pith. sign in

REVIEW 3 major objections 4 minor 12 references

This paper predicts that cosmogenic photon fluxes from realistic cosmic-ray source compositions are more than 1.5 orders of magnitude below current experimental limits, with only proton-rich injection bringing them within reach.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 01:25 UTC pith:PH3JVHV7

load-bearing objection Solid, well-scoped simulation study of cosmogenic photons; the central ordering is robust, but the single-point normalization makes the abstract's '>1.5 orders' margin over-optimistic. the 3 major comments →

arxiv 2607.14670 v1 pith:PH3JVHV7 submitted 2026-07-16 astro-ph.HE

Cosmogenic photon fluxes at ultra-high energies

Marcus Niechciol , Chiara Jane Papior , Markus Risse This is my paper
classification astro-ph.HE
keywords cosmogenic photonsGZK effectultra-high-energy cosmic raysPierre Auger Observatoryphoton flux limitsCRPropacompositionsource scenarios
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper asks whether the ultra-high-energy photons produced when cosmic rays travel through intergalactic space are, in principle, observable. Simulating propagation with CRPropa and normalizing to the Pierre Auger cosmic-ray spectrum, the authors find that for the mixed heavy-nucleus compositions that best describe Auger data, the predicted photon flux is more than 1.5 orders of magnitude below current experimental upper limits. The only scenarios within reach of present and near-future observatories are those with a substantial proton component at the highest energies, and some pure-proton parameter choices are already excluded. They also provide 'proton max' and 'iron min' benchmark fluxes to bracket the uncertainty in the predictions.

Core claim

The paper's central claim is a quantitative prediction: the cosmogenic photon flux at Earth depends sharply on the mass composition and cutoff energy of the injected cosmic rays, and for the mixed compositions that reproduce the Pierre Auger spectrum and composition data, the flux lies more than 1.5 orders of magnitude below present limits. The highest realistic flux comes from a pure-proton scenario with flat spectrum and high cutoff; it is not yet excluded but would be testable with an order-of-magnitude sensitivity improvement. The lowest benchmark, pure iron with a low cutoff, still yields a nonzero EeV photon floor but far below any foreseeable sensitivity. Because the photon yield scal

What carries the argument

The central object is the photon yield, defined as the number of photons above a threshold (typically 10^18 eV) reaching Earth per injected primary, computed with the CRPropa 3.2 propagation code. The paper derives its key scaling from it: for nuclei of mass number A, the yield is approximately N_p(E)/A^2.3, following the superposition model. The absolute flux level is set by normalizing each simulated cosmic-ray scenario to the Pierre Auger-measured flux at 10^19 eV. Six source scenarios (pure proton maximum, proton II, mixed Ia/Ib/II, iron minimum) bracket the range of realistic and benchmark fluxes.

Load-bearing premise

The absolute flux scale in every scenario is set by normalizing the simulated cosmic-ray flux to the measured flux at exactly 10^19 eV; a factor-of-two error in that calibration would shift all predicted photon fluxes by the same factor and could erode the 'more than 1.5 orders of magnitude below limits' margin for mixed compositions.

What would settle it

If the Pierre Auger Observatory with AugerPrime achieves a factor-of-ten improvement in photon sensitivity above 10^18.7 eV and still sees no diffuse flux, then all pure-proton scenarios with flat spectra and high cutoffs (the 'proton max' region of Fig. 10) are excluded. Conversely, a detected diffuse photon flux at the level predicted for the 'mixed Ib' or 'proton II' scenarios would falsify the claim that proton-poor compositions keep the flux more than 1.5 orders below present limits.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • If the paper is right, detection of cosmogenic photons above 10^18 eV with current-generation observatories would require a proton-rich source component; the Auger-fitted proton-poor compositions are effectively invisible.
  • The cutoff-energy degeneracy in the mixed-II scenario is broken by photons: changing the LE component cutoff from 10^19.4 eV to infinity raises the >10^19 eV photon flux by more than 1.5 orders of magnitude, while leaving the cosmic-ray spectrum unchanged.
  • Present photon limits already exclude pure-proton injections with flat spectra (Gamma ≲ 1.8) and high cutoffs (E_max ≳ 10^20.5 eV), as shown in Fig. 10.
  • A floor of cosmogenic photons exists even in the pessimistic pure-iron scenario, though it is far below experimental reach.
  • The photon flux above ~2×10^19 eV is dominated by sources within 50 Mpc, so a single nearby proton-emitting source could dominate the observable signal.

Where Pith is reading between the lines

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

  • If a proton-emitting source sits near the 15 Mpc distance of maximum yield (the paper notes Virgo is close), its individual contribution could exceed the diffuse flux from the entire proton-poor population; a directional search or anisotropy test would be a cheaper way to find it than a diffuse-flux search.
  • The paper's single-point flux normalization at 10^19 eV could be replaced by a full-spectrum fit; doing so might shift the absolute photon flux by factors of order unity and should be tested before concluding the margin is robust.
  • Because the photon yield scales as A^-2.3, the photon flux acts as a magnifying glass on the highest-energy nucleons: measuring the flux just above 10^19 eV would directly probe the tail of the injection spectrum, which is otherwise poorly constrained.
  • The connection between photons and neutrinos from the same GZK interactions means that joint non-observations with future detectors can push the allowed composition toward heavier injection and lower cutoffs more strongly than either messenger alone.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This paper presents CRPropa 3.2 simulations of cosmogenic photon production during UHECR propagation. Sections 3–4 examine how the photon yield depends on primary energy, mass, source distance, spectral index, cutoff, and background models. Section 5 defines six source scenarios: three mixed-composition scenarios derived from Auger spectrum/composition fits (mixed Ia, mixed Ib, mixed II), one pure-proton scenario (proton II), and two extreme benchmarks (proton max, iron min) where the cosmic-ray flux is matched to the edge of the measured spectrum. The main predictions are that mixed scenarios without initial protons lead to photon fluxes more than 1.5 orders of magnitude below current limits, that proton-rich scenarios may be within reach of AugerPrime, and that some pure-proton parameter combinations are already constrained. The paper also shows that cosmogenic photons break a cutoff degeneracy in the mixed II scenario and provides a parameter map for pure-proton exclusion.

Significance. If the results hold, the paper provides updated, well-documented predictions for cosmogenic photon fluxes and a useful assessment of their detectability. Strengths include the use of a publicly available simulation code, the explicit scenario table, the comparison with previous calculations, and falsifiable statements about current and near-future experimental sensitivity. The main value is in quantifying that realistic Auger-consistent mixed compositions imply fluxes far below current limits, while a proton admixture near the GZK threshold could bring the flux within reach. However, the absolute scale is set by a single-point normalization, and the closest mixed scenario sits only about 1.5 orders below the limits; the paper's own quoted factor-of-two systematics therefore make the precise '>1.5 orders' claim fragile.

major comments (3)
  1. [Sec. 2, final paragraph; Sec. 6, first paragraph] The absolute flux scale is fixed by normalizing all scenarios to the Auger cosmic-ray flux at exactly 10^19 eV [9]. The paper itself reports factor ~2 differences with previous calculations [4,5], a factor ~1.5 change when the radio background is varied, and a factor ~2 reduction for a 1 nG magnetic field. For the closest mixed scenario ('mixed Ib'), the margin to the photon limits is only about 1.5 orders (factor ~30); a factor-2 upward shift would reduce this to ~1.2 orders. Thus the abstract's specific claim of 'more than 1.5 orders of magnitude below present limits' is not robust to the documented scale uncertainty. Please either propagate a normalization uncertainty into the photon fluxes or soften the quantitative claim.
  2. [Sec. 5, Table 1 and following paragraph] The sentence 'Other scenarios from [11,12] lead to photon fluxes bracketed by these three scenarios' is asserted without showing the comparison. Since the photon flux is highly sensitive to the proton fraction at the highest energies, this bracketing is load-bearing for the 'realistic range' conclusion. Please show the additional scenarios (or at least their proton fractions and resulting photon fluxes) or provide a bounding argument from the composition fits.
  3. [Sec. 6, Fig. 10] The pure-proton exclusion map scans only E_max and Gamma at the fiducial normalization and source evolution m. The text notes that changing m by ±6 changes the 10^18 eV photon flux by a factor ~1.5, which would shift the exclusion boundary in Fig. 10. The figure should either include this systematic band or explicitly state that the boundary is for the fiducial choice.
minor comments (4)
  1. [Sec. 5, text near Eq. (5.1)] The flux units in the sentence 'amounts here to ~1.3·10^-3 km−2 sr−2 yr−2' should be 'km−2 sr−1 yr−1'.
  2. [Secs. 4–5, Eqs. (4.1), (5.1), (5.2)] Three different cutoff functions are used. Please add a sentence in Sec. 5 explaining why the mixed scenarios use a plateau-shaped cutoff (e^{1-E/Emax} for E>Emax) while the proton scenarios use a pure exponential; otherwise the table is hard to interpret.
  3. [Sec. 5, 'proton II' paragraph] The acknowledged tensions with neutrino limits and with the EeV cosmic-ray flux should be repeated when the scenario is used to draw conclusions, so the reader does not mistake it for a fully realistic benchmark.
  4. [Fig. 7 caption] The experimental limits are labeled by experiment names but no line/color mapping is given in the caption; please clarify which curve corresponds to which limit.

Circularity Check

0 steps flagged

No significant circularity: photon fluxes are new outputs from external input scenarios and benchmark against genuine experimental limits.

full rationale

The derivation chain is self-contained. The mixed-composition scenarios (mixed Ia/Ib from [11], mixed II from [12]) and proton II (from [10]) are external inputs taken from Auger fits; they are not defined in terms of the photon flux this paper predicts. The simulation absolute scale is fixed by normalizing to the measured all-particle cosmic-ray flux at 10^19 eV (Sec. 2), an external calibration that is not derived from the cosmogenic-photon quantities being reported. The predicted integral photon fluxes (Fig. 7) are then compared with measured upper limits [13], an independent benchmark. The 'proton max' and 'iron min' scenarios are explicitly constructed benchmark extremes ('by construction'), so their extreme fluxes make no hidden predictive claim. The only self-citations ([15] for comparison details and [16] for a previously identified technical issue) are not load-bearing: no central conclusion rests on these citations. The sensitivity of the absolute photon scale to the 10^19 eV normalization and to radio-background choices is a systematic uncertainty, not circularity.

Axiom & Free-Parameter Ledger

7 free parameters · 6 axioms · 0 invented entities

The central results depend on (a) the accuracy of the CRPropa propagation code for UHE photon cascades, (b) the adopted background photon fields, (c) the external Auger-fit source scenarios, (d) the assumed homogeneous source distribution and rigidity-dependent cutoff, and (e) the single-point flux normalization at 10^19 eV. These are inputs from prior measurement/literature rather than results of this paper, so they are listed as axioms and parameters rather than claims.

free parameters (7)
  • proton max: Gamma, E_p,max = Gamma=2.3, lg(E_p,max/eV)=21.0
    Hand-chosen to maximize the photon flux while keeping the simulated CR spectrum within the Auger-measured band (Sec. 5).
  • iron min: Gamma, E_p,max = Gamma=2.5, lg(E_p,max/eV)=19.0
    Hand-chosen to minimize the photon flux while keeping the simulated CR spectrum within the Auger-measured band (Sec. 5).
  • benchmark normalization rescale (proton max / iron min) = 20-30% upward / downward
    Rescaled so the simulated cosmic-ray flux scratches the upper (lower) band of the measured spectrum at highest energies; this directly sets the absolute level of the extreme benchmark photon fluxes (Sec. 5).
  • proton II: Gamma, E_p,max, m = Gamma=2.25, lg(E_p,max/eV)=19.75, m=5, cutoff Eq. (5.2)
    Adopted from the Auger-fit scenario in [10]; load-bearing for the 'proton II' flux level but not fitted in this paper.
  • mixed Ia: Gamma, E_p,max, fractions = Gamma=0.96, lg(E_p,max/eV)=18.68, f_He=0.673, f_N=0.281, f_Si=0.046
    From the Auger combined fit [11]; no protons; load-bearing for the central '>1.5 orders below limits' claim.
  • mixed Ib: Gamma, E_p,max, fractions = Gamma=2.04, lg(E_p,max/eV)=19.88, f_N=0.798, f_Si=0.202
    From the Auger combined fit [11]; no protons or helium.
  • mixed II: LE/HE components = LE: Gamma=3.52, lg(E_p,max/eV)>19.4, f_H=0.487, f_He=0.073, f_N=0.44; HE: Gamma=-1.99, lg(E_p,max/eV)=18.15, f_He=0.236,
    From [12]; the LE proton fraction is responsible for the high-energy photon tail and the E_p,max,LE degeneracy band in Fig. 7.
axioms (6)
  • domain assumption CRPropa 3.2 correctly implements UHE photon production and EM-cascade physics (photo-pion, pair production, inverse Compton, synchrotron)
    The paper relies on the CRPropa description of photon interactions via Refs. [3,6]; no independent verification of the code's absolute cross-section implementation is provided.
  • domain assumption Extragalactic background light [7] and radio background [8] models used for propagation are accurate
    The authors test the radio background alternative [14], finding a factor ~1.5 change, so this assumption carries ~50% uncertainty; used in Sec. 2 and Sec. 6.
  • domain assumption Auger-fitted mixed-composition scenarios [11,12] genuinely describe UHECR injection (zero or low proton fractions)
    These fits are external results; the paper's central 'realistic' flux predictions inherit them. If the true injection has a larger proton component, fluxes could be higher (acknowledged in Sec. 6).
  • domain assumption Sources are homogeneously distributed with min distance 1 Mpc, max distance 1000 Mpc, and no strong local proton source
    Sec. 2 states the distance range; Sec. 6 notes a local proton source at ~15 Mpc could boost the flux; the contribution beyond 1000 Mpc is <3% for the studied thresholds.
  • domain assumption E_max of a nucleus scales with charge Z (E_max = Z E_p,max)
    Rigidity-dependent cutoff, standard in UHECR source modeling, stated in the Table 1 caption; affects the nucleon energies available for GZK photon production.
  • domain assumption Superposition model approximates nucleus photon yields as N^A_gamma(E) ≈ A N^p_gamma(E/A) ≈ N^p_gamma(E)/A^2.3
    Used in Sec. 3 to interpret A-scaling; not load-bearing for the absolute flux results, which use full simulations.

pith-pipeline@v1.3.0-alltime-deepseek · 10110 in / 23104 out tokens · 194963 ms · 2026-08-02T01:25:15.125882+00:00 · methodology

0 comments
Cite this review

Pith. "Pith review of Cosmogenic photon fluxes at ultra-high energies." pith.science (2026). https://pith.science/paper/PH3JVHV7

@misc{pith2026260714670,
  author       = {Pith},
  title        = {Pith review of: Cosmogenic photon fluxes at ultra-high energies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PH3JVHV7}},
  note         = {Machine review of arXiv:2607.14670}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

During their propagation, ultra-high energy cosmic rays produce cosmogenic photons. The expected flux level of these photons may vary by orders of magnitude depending on parameters such as the spectrum and composition of cosmic rays at injection or the source distance. We investigate the photon yields for various assumptions on injection parameters. The photon yield is largest for proton-emitting sources at about 15 Mpc distance for cosmogenic photons above $10^{18}$ eV. While the photon yield from protons always exceeds the one from heavier nuclei of the same total energy, the differences are reduced for larger source distances. Then, we quantify the cosmogenic photon fluxes for different source scenarios. We regard mixed-composition scenarios that were found to provide a reasonable description of the data from the Pierre Auger Observatory. In addition, benchmark scenarios leading to comparably high ("maximum") or low ("minimum") photon fluxes are determined assuming pure proton or pure iron primaries. For the mixed-composition scenarios that do not contain initial protons, the predicted cosmogenic photon fluxes are below present experimental limits by more than 1.5 orders of magnitude. In case of a substantial fraction of protons at the highest energies, the related photon flux might be in reach. Certain parameter combinations of pure proton scenarios are constrained already by present photon limits.

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

12 extracted references · 9 linked inside Pith

  1. [1]

    Risse and P

    M. Risse and P. Homola,Search for ultrahigh energy photons using air showers,Mod. Phys. Lett. A22(2007) 749 [astro-ph/0702632]. [2]Pierre AugerCollaboration,Searches for Ultra-High-Energy Photons at the Pierre Auger Observatory,Universe8(2022) 579 [2210.12959]

    Pith/arXiv arXiv 2007

  2. [3]

    Heiter, D

    C. Heiter, D. Kuempel, D. Walz and M. Erdmann,Production and propagation of ultra-high energy photons using CRPropa 3,Astropart. Phys.102(2018) 39 [1710.11406]

    Pith/arXiv arXiv 2018

  3. [4]

    Gelmini, O

    G.B. Gelmini, O. Kalashev and D. Semikoz,Upper Limit on the Diffuse Radio Background from GZK Photon Observation,Universe8(2022) 402 [2206.00408]

    Pith/arXiv arXiv 2022

  4. [5]

    Sarkar, K.-H

    B. Sarkar, K.-H. Kampert and J. Kulbartz,Ultra-High Energy Photon and Neutrino Fluxes in Realistic Astrophysical Scenarios, inProc. 32nd International Cosmic Ray Conference, vol. 2, p. 198, 2011, DOI

    2011

  5. [6]

    Alves Batista et al.,CRPropa 3.2 — an advanced framework for high-energy particle propagation in extragalactic and galactic spaces,JCAP09(2022) 035 [2208.00107]

    R. Alves Batista et al.,CRPropa 3.2 — an advanced framework for high-energy particle propagation in extragalactic and galactic spaces,JCAP09(2022) 035 [2208.00107]

    Pith/arXiv arXiv 2022

  6. [7]

    Gilmore, R.S

    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 (2012) 3189 [1104.0671]

    Pith/arXiv arXiv 2012

  7. [8]

    Protheroe and P.L

    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(1996) 45 [astro-ph/9605119]. [9]Pierre AugerCollaboration,The energy spectrum of cosmic rays beyond the turn-down around10 17 eV as measured with the surface detector of the Pierre Auger Observ...

    Pith/arXiv arXiv 1996

  8. [10]

    Boncioli for thePierre AugerCollaboration,A new view of UHECRs with the Pierre Auger Observatory,PoSUHECR2024(2025) 027 [2509.15862]

    D. Boncioli for thePierre AugerCollaboration,A new view of UHECRs with the Pierre Auger Observatory,PoSUHECR2024(2025) 027 [2509.15862]. [11]Pierre AugerCollaboration,Combined fit of spectrum and composition data as measured by the Pierre Auger Observatory,JCAP04(2017) 038 [1612.07155]. [12]Pierre AugerCollaboration,Constraining the sources of ultra-high-...

    arXiv 2025

  9. [14]

    Nit ¸u, H.T.J

    I.C. Nit ¸u, H.T.J. Bevins, J.D. Bray and A.M.M. Scaife,An updated estimate of the cosmic radio background and implications for ultra-high-energy photon propagation,Astropart. Phys. 126(2021) 102532 [2004.13596]

    Pith/arXiv arXiv 2021

  10. [15]

    Papior, M

    C.J. Papior, M. Niechciol and M. Risse,Estimating the GZK Photon Flux from Extragalactic Cosmic Rays,PoSICRC2025(2025) 353

    2025

  11. [16]

    Bobrikova, M

    A. Bobrikova, M. Niechciol, M. Risse and P. Ruehl,Predicting the UHE photon flux from GZK-interactions of hadronic cosmic rays using CRPropa 3,PoSICRC2021(2021) 449. – 14 – [17]Pierre AugerCollaboration,Search for photons above10 18 eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory, Ph...

    Pith/arXiv arXiv 2021

  12. [18]

    Schmidt for thePierre AugerCollaboration,AugerPrime: Status and first results,PoS ICRC2025(2025) 385 [2508.08056]

    D. Schmidt for thePierre AugerCollaboration,AugerPrime: Status and first results,PoS ICRC2025(2025) 385 [2508.08056]. – 15 –

    Pith/arXiv arXiv 2025