Pith. sign in

REVIEW 2 major objections 5 minor 19 references

Constraints on the spread of nuclear masses in ultra-high-energy cosmic rays based on the Phase I hybrid data from the Pierre Auger Observatory

T0 review · 2 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read The correlation of shower maximum depth with ground signal in the Phase I hybrid data constrains the cosmic-ray mass spread to $\sigma(\ln A)\in[1.0,1.7]$ below the ankle, ruling out two-neighbouring-mass-group beams.

desk verdict Solid incremental Auger update on UHECR mass spread, credible but overclaims model-independence without showing all model conversions. read the letter →

arxiv 2507.08997 v1 pith:EU2MDHBH submitted 2025-07-11 astro-ph.HE hep-ph

Alexey Yushkov (for the Pierre Auger Collaboration) This is my paper · ORCID
classification astro-ph.HEhep-ph
keywords ultra-high-energycosmicraysmasscompositionair-showerobservableshadronicinteractionmodelsrankcorrelationanklePierreAugerObservatory
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to establish that the spread of nuclear masses in ultra-high-energy cosmic rays below the ankle is large, with $\sigma(\ln A)$ between 1.0 and 1.7, and that this conclusion is nearly independent of which hadronic interaction model is used to simulate air showers. It does so with a rank correlation between two scaled shower observables, the depth of shower maximum $X^*_{\max}$ and the water-Cherenkov signal $S^*_{38}$, measured in hybrid events that trigger both the fluorescence and surface detectors. If the claim holds, it rules out all pure compositions and all mixtures made of only two neighbouring mass groups near the ankle, giving a nearly model-independent handle on the cosmic-ray beam. The same correlation, studied as a function of energy and zenith angle, is used to test the newest post-LHC hadronic models: the mass mixes inferred from $X_{\max}$ fits alone are incompatible with the data for the Sibyll family of models, while EPOS LHC-R, with its deeper $X_{\max}$ scale, produces a larger mass spread that matches the observation.

What carries the argument

The central object is $r_G$, a rank correlation coefficient resistant to outliers, computed between $X^*_{\max}$ and $S^*_{38}$, the depth of shower maximum and the water-Cherenkov signal at 1000 m from the core, rescaled to a reference energy of 10 EeV and a zenith angle of 38°. The argument rests on two nearly model-independent facts: mass groups separate cleanly in the $(X^*_{\max},S^*_{38})$ plane, and the fluctuations of both observables are nearly identical across hadronic interaction models. This makes $r_G$ a monotonic proxy for the spread of primary masses $\sigma(\ln A)$: positive for pure beams, increasingly negative as the mix broadens, and most negative for a 50/50 proton-iron beam. The conversion from $r_G$ to $\sigma(\ln A)$ is done with simulated mixtures of proton, helium, oxygen and iron nuclei in 0.05 steps, using the newest post-LHC models EPOS LHC-R, QGSJet-III-01, Sibyll⋆ and Sibyll 2.3e.

What would settle it

A direct check is to measure $r_G$ in the same energy bin with an independent systematic treatment: if a reanalysis of the same Phase I events returned $r_G \approx 0$ rather than $-0.094 \pm 0.024$, the inferred $\sigma(\ln A)\in[1.0,1.7]$ would not follow. Alternatively, if a future hadronic interaction model that reproduces the observed $X_{\max}$ distributions predicts a pure-proton or pure-iron $r_G$ equal to the measured value, the exclusion of the neighbouring two-group mixes would be falsified.

Watch

Extended reading notes

Core claim

The central claim is that the rank correlation coefficient $r_G$ between $X^*_{\max}$ and $S^*_{38}$ is negative below the ankle and translates into a well-defined degree of mass mixing. In the energy bin $\lg(E/\mathrm{eV})\in[18.5,18.6]$ the data give $r_G = -0.094 \pm 0.024$, compatible with $\sigma(\ln A)\in[1.09,1.74]$ under Sibyll 2.3e and $\sigma(\ln A)\in[0.97,1.61]$ under EPOS LHC-R; across $\lg(E/\mathrm{eV})\in[18.3,18.7]$ the spread remains within $\sigma(\ln A)\in[1.0,1.7]$. Because the mass-group separation in the $(X^*_{\max},S^*_{38})$ plane and the size of the observable fluctuations are nearly model-independent, all pure compositions and all two-neighbouring-mass-group mixes are excluded below the ankle for every hadronic model considered. Above the ankle the mixing degree falls to $\sigma(\ln A)\in[0.0,1.3]$. The energy and zenith-angle dependence of $r_G$ further shows that $X_{\max}$-based composition fits are inconsistent with the observed correlation for the Sibyll models, while EPOS LHC-R, whose deeper $X_{\max}$ scale yields $\sigma(\ln A)\in[1.2,1.4]$, agrees with the data.

Load-bearing premise

The analysis assumes that the correlation between the two scaled shower observables, and how much they fluctuate, is almost the same in every hadronic interaction model; if a future model changes that correlation, the inferred $\sigma(\ln A)$ range and the exclusion of two-group mixes would become model-dependent.

Editorial extensions

If this is right

  • Below the ankle, the cosmic-ray beam must contain at least three distinct mass groups, so source models that produce only proton-helium, only helium-CNO, or only CNO-iron mixtures are excluded.
  • Because the $\sigma(\ln A)$ constraints are nearly model-independent, any mismatch between a hadronic model's predicted $r_G$ and the observed correlation tests the model itself rather than the composition inference.
  • Composition fits based on $X_{\max}$ alone are not sufficient below the ankle: they must also reproduce the observed $r_G$, and among current models only those with a deeper $X_{\max}$ scale, such as EPOS LHC-R, do so.
  • Above the ankle the inferred spread narrows to $\sigma(\ln A)\in[0.0,1.3]$, indicating that the beam becomes more homogeneous at the highest energies.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable consequence of the model-independence claim is that a future hadronic model which reproduces the observed $X_{\max}$ distributions but not the observed $r_G$ would reveal its deficiency in the muon or ground-signal sector, pointing to adjustments in the $X_{\max}$ scale or fluctuations rather than in composition.
  • Because $r_G$ is a rank statistic, the central constraint should be insensitive to global energy-scale offsets; a natural cross-check would be to repeat the analysis with a different reference energy and verify that the inferred $\sigma(\ln A)$ interval does not shift.
  • If the narrowing of $\sigma(\ln A)$ above the ankle persists in future data, the result could sharpen astrophysical inferences about the transition between galactic and extragalactic cosmic-ray source populations.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The paper presents a preliminary update of the Pierre Auger correlation analysis between the depth of shower maximum Xmax* and the surface-detector signal S38* (a proxy for the muon number), using Phase I hybrid data (2004–2021, 9661 events) and the latest post-LHC hadronic interaction models (EPOS LHC-R, QGSJet-III-01, Sibyll*, Sibyll 2.3e). The rank correlation coefficient rG is measured as a function of energy and zenith angle. Using a forward grid of simulated mixtures of (p, He, O, Fe), the authors convert rG into constraints on sigma(ln A). Below the ankle they quote sigma(ln A) in [1.0, 1.7], which excludes mixes consisting only of neighboring mass groups (p-He, He-CNO, CNO-Fe); above the ankle they quote sigma(ln A) in [0.0, 1.3]. They compare the observed rG with predictions derived from Xmax-fraction fits ('FD Xmax mix') and find that Sibyll 2.3e and Sibyll* are incompatible with the data, while EPOS LHC-R is compatible because of its deeper Xmax scale.

Significance. If the model-independence claim holds, this correlation analysis is a powerful composition diagnostic that largely bypasses the hadronic-model systematic uncertainty that plagues absolute Xmax measurements. The paper's strengths include a clear description of the rank-correlation method, a forward simulation grid with no parameters fitted to the observed rG, the use of a Phase I data set roughly 2.8 times larger than the previous one, and explicit estimates of statistical and systematic uncertainties on rG. The resulting sigma(ln A) constraints and the apparent tension with the Sibyll-family predictions are timely and relevant for current UHECR source and composition studies. However, the central generalization that the constraints hold for any hadronic interaction model is not supported by the evidence shown for the full model set, and this weakens the claim of model-independent exclusion of two-neighboring-group mixes.

major comments (2)
  1. [Section 2, Fig. 3, and Summary (Section 4)] The claim that the sigma(ln A) constraints hold for 'any hadronic interaction model' is not quantitatively established. The rG-to-sigma(ln A) conversion is shown only for Sibyll 2.3e and EPOS LHC-R, yielding allowed intervals [1.09, 1.74] and [0.97, 1.61] at lg(E/eV) = 18.5–18.6, i.e. a lower-bound difference of 0.12. The maximum sigma(ln A) attainable by a two-neighboring-group mix is about 0.69–0.70, so the margin relative to the quoted lower bound is only about 0.3. If QGSJet-III-01 or Sibyll* shifts the conversion in the opposite direction by an amount comparable to the observed Sibyll–EPOS difference, the exclusion of p-He, He-CNO, and CNO-Fe mixes would become model-dependent and could fail. The paper should either show the conversion for all four models or explicitly quantify the model-to-model spread in the conversion and propagate that uncertainty into the quoted sigma(ln A) interval before asserting model independence.
  2. [Section 2, Fig. 4] The energy evolution of the sigma(ln A) constraints in Fig. 4 is presented without any uncertainty band, so it is unclear whether the quoted [1.0, 1.7] interval includes the +0.01/−0.02 systematic uncertainty on rG and the finite grid spacing of 0.05 in the mixture fractions. The text states that experimental systematics have a minor impact, but the propagation to sigma(ln A) is not shown. A quantitative propagation of both the rG systematic uncertainty and the grid granularity would make the central interval claim more robust.
minor comments (5)
  1. [Section 1] The term 'Phase I hybrid data' is used in the abstract and introduction but is not defined; please specify the data-taking period and the meaning of 'Phase I' (here it appears to be 12/2004–12/2021).
  2. [Section 2, Fig. 2] The text states that the observed energy dependence 'hold[s] for all hadronic interaction models used in this study', but Fig. 2 shows only Sibyll 2.3e and EPOS LHC-R. Please either show the corresponding curves for QGSJet-III-01 and Sibyll* or soften the wording.
  3. [Section 2] The sentence describing the behavior of rG with mass spread should clarify that it refers to simulated rG values; the data are a single measurement and are negative in the relevant energy range.
  4. [Section 2, Fig. 3] Please specify how the allowed sigma(ln A) interval is derived from the grid: which simulated points are considered compatible with the observed rG, and whether the interval corresponds to 1-sigma statistical compatibility only.
  5. [Summary (Section 4)] The statement that the conclusions 'hold for any hadronic interaction model, including their pre-LHC versions' is stronger than what the current figures demonstrate; given the preliminary nature of the analysis, a more cautious wording would be appropriate.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: measured rG is mapped through a forward-simulated mass-mix grid to sigma(ln A), and the model-independence claim rests on explicit model comparisons plus a prior, independent published Auger analysis.

full rationale

The derivation chain is self-contained. rG is computed directly from the Phase I hybrid data (Section 2, Fig. 1) and is not fitted to sigma(ln A). The sigma(ln A) intervals are obtained by forward-simulating (p, He, O, Fe) mixtures in 0.05 steps with fixed hadronic models and reading off the simulated rG values compatible with the observed rG (Fig. 3); no parameter is adjusted to the observed correlation. The exclusion of p-He, He-CNO and CNO-Fe mixtures follows arithmetically from the fact that such two-neighboring-group mixtures have maximum sigma(ln A) of about 0.7, below the quoted lower bound of about 1.0, so it is not an input assumption. The model-independence claim is supported in this paper by explicit rG-to-sigma conversions for Sibyll 2.3e ([1.09, 1.74]) and EPOS LHC-R ([0.97, 1.61]) at lg(E/eV) = 18.5-18.6, and by the energy and zenith-angle comparisons in Figs. 2, 4 and 5. The extension to pre-LHC model versions cites the previous Auger publication [3], which is an independent, earlier-dataset peer-reviewed result rather than an unverified self-citation; no uniqueness theorem is imported and no ansatz is smuggled in via citation. The residual model spread in the conversion is a model-dependence/correctness concern, not circularity, and the paper's statement that uncertainties of the FD Xmax mix will be propagated in a forthcoming publication does not indicate a circular step.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The analysis introduces no free parameters fitted to the measured rG: the scaling to reference energy and zenith angle uses standard attenuation and elongation inputs, and the sigma(ln A) conversion is a forward grid. All axioms are domain assumptions about the validity of hadronic models and the model independence of the correlation, plus one standard mathematical property of rank correlations. No new entities are postulated.

assumptions (6)
  • domain assumption Post-LHC hadronic interaction models (EPOS LHC-R, QGSJet-III-01, Sibyll*, Sibyll 2.3e) provide valid predictions for Xmax and muon content at ultra-high energies.
    Section 1: the simulation library uses these models; all interpretations of data in terms of mass composition rely on their predictions.
  • domain assumption The correlation between Xmax* and S38* is nearly model-independent for mass-group separation and fluctuations.
    Section 2: 'nearly model-independent separation between mass groups in the (Xmax*, S38*) plane...' This is the load-bearing premise for the claimed model independence of sigma(ln A) constraints.
  • domain assumption S(1000) is a valid proxy for the muon number at ground.
    Introduction and Section 2: 'Using the signal in water-Cherenkov detectors at 1000 meters from the shower core as a proxy for the muon number'. The correlation's mass sensitivity relies on this.
  • domain assumption The selection cuts (in particular theta >= 35 degrees) do not bias the measured correlation.
    Section 1 and 3: 'In the correlation analysis, such bias is found to be minimal...' The justification for the 35-degree cut is based on simulations, not data.
  • standard math Rank correlation coefficients are invariant under monotonic transformations, so systematic scale offsets in Xmax or S(1000) do not affect rG.
    Section 2: 'Rank correlation coefficients are invariant to the absolute values of S(1000) and Xmax...'
  • domain assumption The grid of simulated mixtures with (p, He, O, Fe) in 0.05 steps adequately spans the possible mass compositions.
    Section 2 and Fig. 3: the conversion from rG to sigma(ln A) uses this grid; intermediate mass groups are not simulated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Constraints on the spread of nuclear masses in ultra-high-energy cosmic rays based on the Phase I hybrid data from the Pierre Auger Observatory." pith.science (2026). https://pith.science/paper/EU2MDHBH

@misc{pith2026250708997,
  author       = {Pith},
  title        = {Pith review of: Constraints on the spread of nuclear masses in ultra-high-energy cosmic rays based on the Phase I hybrid data from the Pierre Auger Observatory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EU2MDHBH}},
  note         = {Machine review of arXiv:2507.08997}
}
abstract

We present an analysis of the correlation between the depth of the maximum of air-shower profiles and the signal in water-Cherenkov stations in events registered simultaneously by the fluorescence and surface detectors of the Pierre Auger Observatory. The analysis enables us to place constraints on the spread of nuclear masses in ultra-high-energy cosmic rays with a minor impact from the experimental systematic uncertainties and uncertainties in air-shower simulations. Due to this unique feature, the correlation analysis has previously allowed us to exclude all pure and proton-helium compositions near the ankle in the cosmic-ray energy spectrum at 5$\sigma$ confidence level. The same property makes the correlation analysis an effective tool for testing the consistency of predictions of the post-LHC hadronic interaction models, including their latest versions such as EPOS LHC-R, QGSJet-III-01, Sibyll${}^\bigstar$ and Sibyll 2.3e. In this work, the correlation analysis using the Phase I hybrid data from the Pierre Auger Observatory is presented. The analysis uses the newest generation of hadronic interaction models and covers an extended energy range around the ankle in the cosmic-ray spectrum.

Figures

Figures reproduced from arXiv: 2507.08997 by the authors.

Figure 1
Figure 1. Correlation between 𝑋 ∗ max and 𝑆 ∗ 38. Left panel: proton and iron showers (samples of 2000 events each) simulated with Sibyll 2.3e; the legend contains 𝑟G values for pure beams and proton-iron equal mix (maximum mixing degree). Right panel: correlation in the data. Energy range: lg(𝐸/eV) ∈ [18.5, 18.6]; zenith angle range: 𝜃/deg ∈ [35, 60]. 18.4 18.6 18.8 19.0 19.2 19.4 lg(E/ eV) −0.5 −0.4 −0.3 −0.2 −0.1 0.0 0.1 0… view at source ↗
Figure 2
Figure 2. Energy dependence of the correlation in the data compared with correlations in simulations with Sibyll 2.3e and EPOS LHC-R for protons, iron nuclei and proton-iron equal mix. Additionally, the correlation expected for the mass composition inferred from fits of the FD 𝑋max distributions [17, 18] is shown. 𝜎(ln 𝐴) ∈ [0.0, 1.3] (also for lg(𝐸/eV) ∈ [19.0, 19.5], not shown here). These conclusions are valid for previous… view at source ↗
Figure 3
Figure 3. Conversion from 𝑟G to 𝜎(ln 𝐴) using simulations with Sibyll 2.3e (left) and EPOS LHC-R (right). The shaded area represents the observed value with its statistical errors. The 𝜎(ln 𝐴) ranges of the simulated mixes compatible with the correlation in the data are shown in each panel. See text for more details. 18.3 18.4 18.5 18.6 18.7 18.8 18.9 19.0 lg(E/ eV) 0.0 0.5 1.0 1.5 2.0 ) A (ln σ Sibyll2.3e θ/deg∈[35,60] Auger… view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Constraints on 𝜎(ln 𝐴) from the correlation analysis as a function of energy using Sibyll 2.3e and EPOS LHC-R. Dashed lines indicate the values for pure compositions and the proton-iron equal mix (maximum mixing degree). 10 20 30 40 50 60 θ/ deg −0.6 −0.5 −0.4 −0.3 −0.…
Figure 5
Figure 5. Figure 5: Zenith-angle dependence of 𝑟G. The 𝑟G values in the data are compared with Sibyll⋆ (left) and EPOS LHC-R (right) simulations for protons, iron nuclei, proton-iron equal mix, and the FD 𝑋max mix. Energy range: lg(𝐸/eV) ∈ [18.5, 18.7]. 5 [PITH_FULL_IMAGE:figures/full_fi…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 13 canonical work pages

  1. [1]

    Sensitivity of the correlation between the depth of shower maximum and the muon shower size to the cosmic ray composition

    P. Younk and M. Risse,Sensitivity of the correlation between the depth of shower maximum and the muon shower size to the cosmic ray composition,Astropart.Phys. 35(2012) 807 [1203.3732]

    work page Pith review arXiv 2012

  2. [2]

    Pierre Auger Collaboration,The Pierre Auger Cosmic Ray Observatory, Nucl. Instrum. Meth. A798(2015) 172 [1502.01323]

    arXiv 2015

  3. [3]

    Lett.B762(2016) 288 [1609.08567]

    Pierre Auger Collaboration,Evidence for a mixed mass composition at the ‘ankle’ in the cosmic-ray spectrum,Phys. Lett.B762(2016) 288 [1609.08567]

    arXiv 2016

  4. [4]

    T. Fitoussi,Depth of maximum of air-shower profiles above1017.8 eV measured with the fluorescence detector of the Pierre Auger Observatory and mass-composition implications, PoS ICRC2023(2023) 319

    work page 2023

  5. [5]

    Pierre Auger Collaboration,Reconstruction of events recorded with the surface detector of the Pierre Auger Observatory,JINST 15(2020) P10021 [2007.09035]

    arXiv 2020

  6. [6]

    Santos,Update on the Offline Framework for AugerPrime and production of reference simulation libraries using the VO Auger grid resources, PoS ICRC2023(2023) 248

    E. Santos,Update on the Offline Framework for AugerPrime and production of reference simulation libraries using the VO Auger grid resources, PoS ICRC2023(2023) 248

    work page 2023

  7. [7]

    Heck et al.,CORSIKA: A Monte Carlo code to simulate extensive air showers, Report FZKA-6019 (1998)

    D. Heck et al.,CORSIKA: A Monte Carlo code to simulate extensive air showers, Report FZKA-6019 (1998)

    work page 1998

  8. [8]

    Pierre Auger Collaboration,The Offline Software Framework of the Pierre Auger Observatory, Nucl. Instrum. Methods Phys. Res.580 (2007) 1485

    work page 2007

Show all 19 references
  1. [9]

    Pierog and K

    T. Pierog and K. Werner,EPOS LHC-R : up-to-date hadronic model for EAS simulations, PoS ICRC2023(2023) 230

    2023

  2. [10]

    Ostapchenko,QGSJET-III model of high energy hadronic interactions: The formalism, Phys

    S. Ostapchenko,QGSJET-III model of high energy hadronic interactions: The formalism, Phys. Rev. D109 (2024) 034002

    2024

  3. [11]

    Ostapchenko,QGSJET-III model of high energy hadronic interactions

    S. Ostapchenko,QGSJET-III model of high energy hadronic interactions. II. Particle production and extensive air shower characteristics, Phys. Rev. D109 (2024) 094019. 7 Spread of nuclear masses in UHECRs from the Pierre Auger Observatory data Alexey Yushkov

    2024

  4. [12]

    Riehn, R

    F. Riehn, R. Engel and A. Fedynitch,Sibyll ⋆: ad-hoc modifications for an improved description of muon data in extensive air showers, PoS ICRC2023(2023) 429 [2309.05390]

    2023 arXiv

  5. [13]

    Riehn et al.,Hadronic interaction model Sibyll 2.3d and extensive air showers,Phys

    F. Riehn et al.,Hadronic interaction model Sibyll 2.3d and extensive air showers,Phys. Rev. D 102 (2020) 063002 [1912.03300]

    2020 arXiv

  6. [14]

    Pierog,EPOS LHC-R: A Global Approach to Solve the Muon Puzzle, Presented at FZU – Institute of Physics of the Czech Academy of Sciences (2025), https://indico.fzu.cz/event/303/

    T. Pierog,EPOS LHC-R: A Global Approach to Solve the Muon Puzzle, Presented at FZU – Institute of Physics of the Czech Academy of Sciences (2025), https://indico.fzu.cz/event/303/

    2025

  7. [15]

    Pierre Auger Collaboration,Measurement of the cosmic-ray energy spectrum above 2.5×1018 eV using the Pierre Auger Observatory,Phys. Rev. D102(2020) 062005 [2008.06486]

    2020 arXiv

  8. [16]

    Gideon and R

    R. Gideon and R. Hollister,A rank correlation coefficient resistant to outliers, Journal of the American Statistical Association82 (1987) 656

    1987

  9. [17]

    O. Tkachenko,Studies of the mass composition of cosmic rays and proton-proton interaction cross-sections at ultra-high energies with the Pierre Auger Observatory, PoS ICRC2023 (2023) 438

    2023

  10. [18]

    Tkachenko,Measurement of the Inelastic Proton-Proton Cross Section at√𝑠≥ 40TeV using the Hybrid Data of the Pierre Auger Observatory, these proceedings

    O. Tkachenko,Measurement of the Inelastic Proton-Proton Cross Section at√𝑠≥ 40TeV using the Hybrid Data of the Pierre Auger Observatory, these proceedings

  11. [19]

    Tor Vergata

    Pierre Auger Collaboration,Testing hadronic-model predictions of depth of maximum of air-shower profiles and ground-particle signals using hybrid data of the Pierre Auger Observatory, Phys. Rev. D109 (2024) 102001 [2401.10740]. 8 Spread of nuclear masses in UHECRs from the Pie...

    2024 arXiv

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.