Consequences of a Heavy-Metal Scenario of Ultra-High-Energy Cosmic Rays

Alena Bakalov\'a; Ana Laura M\"uller; Jakub V\'icha; Maximilian Stadelmaier; Olena Tkachenko

REVIEW 4 major objections 4 minor 24 references

Consequences of a Heavy-Metal Scenario of Ultra-High-Energy Cosmic Rays

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read Assuming the cosmic-ray beam above 40 EeV is pure iron, the paper shows that a single constant shift of the predicted shower-depth scale makes the ultra-high-energy data consistent.

desk verdict A clearly framed heavy-iron scenario with fitted Xmax shifts and a four-species composition model; the main caveat is that the constant shift is fitted to high-energy data and then applied everywhere else. read the letter →

arxiv 2502.08747 v2 pith:ETDXOTLE submitted 2025-02-12 astro-ph.HE

Jakub Vícha , Alena Bakalová , Olena Tkachenko , Ana Laura Müller , Maximilian Stadelmaier This is my paper · ORCID

classification astro-ph.HE

keywords ultra-high-energycosmicraysheavy-metalscenarioironcompositiondepthofshowermaximumXmaxmasshadronicinteractionmodelsmuonproblemGalacticmagneticfieldbacktracking

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 proposes an extreme scenario for ultra-high-energy cosmic rays: above 40 EeV the arriving beam is pure iron nuclei, and the depth of shower maximum $X_{\rm max}$ predicted by hadronic-interaction models is allowed a constant rescaling shift. Under these assumptions the measured $X_{\rm max}$ moments from the Pierre Auger Observatory become consistent for both QGSJet II-04 and Sibyll 2.3d, with fitted shifts of $\Delta X_{\rm max}=52\pm1^{+11}_{-8}\,\mathrm{g\,cm^{-2}}$ and $29\pm1^{+12}_{-7}\,\mathrm{g\,cm^{-2}}$, respectively. The same scenario yields a four-species mass-composition model, explains the spectral instep by the fading of nitrogen nuclei, reduces the muon deficit for QGSJet II-04, and, after magnetic-field backtracking, restricts possible extragalactic source directions. A sympathetic reader should care because it shows that one systematic shortcoming in the predicted $X_{\rm max}$ scale can simultaneously resolve several longstanding tensions in cosmic-ray data.

What carries the argument

The load-bearing device is a constant shift $X_{\rm max}\to X_{\rm max}+\Delta X_{\rm max}$ applied to the hadronic-model predictions. $X_{\rm max}$ is the atmospheric depth (in g cm$^{-2}$) at which an air shower reaches its maximum number of particles, a mass-sensitive observable. The shift is calibrated by fitting the mean $X_{\rm max}$ above 40 EeV to pure-iron expectations for each interaction model, while the measured $X_{\rm max}$ fluctuations certify that pure iron is the right assumption. Once fixed, the shift is applied to template distributions of four primary species, which are fitted to the public $X_{\rm max}$ distributions; the resulting fraction parameterizations are then used to decompose the energy spectrum, reinterpret muon measurements, and backtrack arrival directions.

What would settle it

Measure $X_{\rm max}$ of individual showers with an independent technique that does not rely on the same hadronic models (for example, direct fluorescence observations) and compare the energy dependence of the DNN-based $X_{\rm max}$ with the shifted model predictions; a scale offset that changes with energy or with the inferred primary mass below 40 EeV would rule out the constant-shift scenario.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central claim is that current hadronic-interaction models get the absolute position of $X_{\rm max}$ wrong by a single constant, while correctly predicting the elongation rate and fluctuations. Fixing this gauge freedom by demanding pure iron above $10^{19.6}$ eV, the paper derives the required shift for QGSJet II-04 and Sibyll 2.3d, then applies it at all lower energies down to about $10^{18.1}$ eV to fit the measured $X_{\rm max}$ distributions and obtain primary fractions of protons, helium, nitrogen, and iron. The resulting picture makes the measured mean and variance of $\ln A$ interpretable, connects the 15 EeV spectral instep to nitrogen fading, reduces the muon problem from about 50% to about 20%, and, after backtracking the highest-energy events in the Galactic magnetic field, points to sources mostly around the Galactic anti-center while the dipole anisotropy can be reproduced at the $2\sigma$ level.

Load-bearing premise

The load-bearing premise is that a single, constant shift of the predicted $X_{\rm max}$ scale applies at all energies and for all primary species, even though it is fitted only above 40 EeV and no physical mechanism for such a shift is provided.

Editorial extensions

If this is right

A single energy- and mass-independent deficit in the predicted $X_{\rm max}$ scale would explain why composition moments, $X_{\rm max}$ distributions, and the muon signal cannot be simultaneously described by unshifted models.
The energy spectrum's instep near 15 EeV would be a composition landmark: the rigidity at which nitrogen nuclei fade from the beam matches the iron suppression, suggesting a common origin tied to magnetic rigidity $E/Z$.
The muon deficit of QGSJet II-04 would shrink from roughly 50% to 20%, indicating that part of the muon problem is actually a misplaced $X_{\rm max}$ scale rather than missing muon production.
If the scenario is right, high-energy cosmic-ray sources are mostly visible in the Galactic anti-center direction, with low-longitude directions shadowed by the Galactic magnetic field.

Reading between the lines

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

A direct test would be to compare the DNN-inferred $X_{\rm max}$ scale with fluorescence-telescope measurements: if the offset between the two techniques varies with energy or primary mass, the constant-shift assumption would fail.
The scenario predicts that a future hadronic-interaction model that deepens $X_{\rm max}$ by the same constant while preserving fluctuations would automatically resolve the current tensions; conversely, a model that cannot be shifted without breaking elongation-rate agreement would disprove the gauge freedom.
If the composition fractions derived here are correct, the neutrino and photon fluxes expected from ultra-high-energy cosmic-ray sources would change because the per-nucleon energy budget shifts with the heavier primary mix, a prediction that next-generation observatories could check.

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.

Desk Editor's Note

read the letter

The paper does something useful: it takes the long-standing tension between Xmax moments and composition seriously and asks what would happen if the absolute Xmax scale from hadronic models is simply wrong by a constant. The package is new as a whole – fitted ΔXmax values for QGSJet II-04 and Sibyll 2.3d (52 and 29 g/cm2), a four-species composition model derived from public Auger data, and consequences for the energy spectrum, muon production, and GMF backtracking. The assumptions are stated up front and the paper does not oversell the scenario; the figures are marked PRELIMINARY and the text is careful about what is assumed versus derived.

The fits themselves look internally consistent and use public data with standard methods. The muon-problem alleviation from ~50% to ~20% for QGSJet II-04 is a concrete, checkable consequence. The backtracking exercise is honest about the strong GMF shadowing. The paper also gives credit where due: the scale-shift idea already appears in refs [9] and [13], and the elongation-rate and fluctuations are taken from the models, not re-fitted.

The load-bearing weakness is exactly what the stress-test note says: the constant ΔXmax is fitted only above 40 EeV, then applied down to 10^18.1 eV and to all four primaries. The paper provides no physical mechanism or independent calibration for this constancy, and if the true offset depends on energy or species, the composition fractions in Fig. 2 are not determined by the data. I also agree that the agreement between the shifted iron curve and the measured mean Xmax is by construction – the shift is chosen to make that happen. The composition fractions are fitted to the same Xmax distributions and then used to 'predict' the spectrum components, so the consistency claim is weaker than it appears. These are not fatal flaws for a scenario paper, but they should be stated in the abstract or summary as limits, not buried in the middle.

I would send this to a serious referee. The constancy assumption and the absence of systematic uncertainties on the composition fractions are exactly what a referee should probe. It is a well-written, honest scenario study that deserves engagement, not a desk rejection. I would cite it if I worked on UHECR phenomenology, and I would probably bring it to a reading group for the discussion of what counts as a consistent interpretation in this field.

Referee Report

4 major / 4 minor

Summary. The paper proposes an 'extreme' heavy-metal scenario in which the cosmic-ray beam above about 40 EeV is assumed to consist of pure iron nuclei and in which all hadronic interaction models have a constant but incorrect absolute Xmax scale. The authors fit the required constant shift DeltaXmax for QGSJet II-04 (52 g/cm2) and Sibyll 2.3d (29 g/cm2) to the Auger DNN-measured mean Xmax above 10^19.6 eV, then apply this same shift to Xmax distributions down to 10^18.1 eV in order to fit four-species composition fractions. From those fractions they derive per-species energy spectra, discuss the muon problem, and backtrack the arrival directions of events above 78 EeV in the Galactic magnetic field. The paper is framed as a consistency demonstration of an extreme scenario, not as a best-fit standard model.

Significance. If the scenario were established, it would offer a single coherent interpretation of the Xmax moments, composition, spectral features, and muon excess, and it would restrict possible source directions. The paper is transparent about its assumptions and uses public Auger data, which is a strength, and the GMF backtracking section adds a concrete phenomenological consequence. However, the central consistency is partly by construction: DeltaXmax is fitted to the very mean Xmax it is then used to reproduce, and its constancy in energy and primary species is assumed rather than tested. The derived composition and spectral decomposition are therefore conditional on an untested assumption. The paper is a useful scenario-level consistency study, but it does not yet establish the heavy-metal interpretation as a quantitatively supported mass-composition model.

major comments (4)

[Section 3] The values DeltaXmax = 52 g/cm2 and 29 g/cm2 are obtained by fitting the pure-iron prediction to the measured mean Xmax above 10^19.6 eV. The subsequent agreement of the shifted iron curve with the data is therefore enforced by the fit, not an independent success of the scenario. The non-tautological content in Fig. 1 is limited to the already-known fact that the Xmax fluctuations are compatible with pure iron and that sigma^2(ln A) becomes consistent with the model-independent estimator of Ref. [13]. The paper should explicitly separate these fitted elements from the predicted ones so that the 'consistent interpretation' claim is not overstated.
[Section 4, top panels of Fig. 2] The shift DeltaXmax is derived using only iron nuclei at E > 10^19.6 eV, yet it is applied to all four primary species and to energies down to 10^18.1 eV without any check of energy or species independence. If the true offset varies with energy or with primary mass, the template fit in Section 4 can absorb the wrong shift into the composition fractions, so the fractions shown in Fig. 2 are not determined by the data alone. The authors should test constancy, for example by fitting DeltaXmax separately in lower-energy bins, or by comparing the shifted QGSJet II-04 and Sibyll predictions with the unshifted measured Xmax moments below 10^19.6 eV, and they should show how the inferred fractions change under a conservative range of DeltaXmax.
[Sections 4 and 5] The composition fractions in Fig. 2 are presented without statistical or systematic uncertainties, and the functional form of the parametrization (Gaussian multiplied by exponential, with the sum normalized to 1) is imposed rather than derived from a physical model. As a consequence, the claim that the instep feature near 15 EeV is 'related to the fading of nitrogen nuclei' is a property of the chosen parametrization rather than an independent prediction. The paper should report goodness-of-fit values for the Xmax-distribution fits, parameter covariances, and a stability check against alternative parametrizations or against relaxing the normalization constraint.
[Section 6 and Fig. 3 (top-left panel)] The statement that the muon problem is alleviated from about 50% to about 20% is not sufficiently specified. A constant shift of the Xmax scale is not a modification of hadronic interactions, so it is unclear from the text and figure whether the comparison is a genuine prediction of the muon signal for fixed primary energy and mass, or whether it is a consistency loop in which the primary mass is first inferred from the shifted Xmax and then used to evaluate muon expectations. The exact quantity plotted and the role of DeltaXmax in the muon prediction need to be defined precisely.

minor comments (4)

[Figures 1-3] Most figures are marked PRELIMINARY. If this manuscript is intended as a journal publication rather than a proceedings contribution, the preliminary labels should be removed or the figures should be replaced by final versions with the same analysis.
[Reference [15]] The reference 'Eur. Phys. J. C210 (2020) 751' appears to contain a typo in the volume number; it should be checked against the published version.
[Section 2, bullet list] The statement that the measured Xmax fluctuations above 10^19.6 eV are consistent with pure iron with p(chi2 > 0.5) would benefit from a definition of the chi2 statistic and the number of degrees of freedom, since the DNN moments in Ref. [11] are not independent of the model assumptions used in their reconstruction.
[Section 3, last paragraph] The comparison with the model-independent sigma^2(ln A) from Ref. [13] is made only for the energy range 3-10 EeV, while the DeltaXmax shift is derived above 40 EeV; the paper should state whether this comparison is meant to constrain the shift at lower energies and what the implied uncertainty is.

Circularity Check

1 steps flagged · score 2.0 of 10

Only the 'consistent interpretation' of the mean Xmax after applying the fitted ΔXmax is tautological; external cross-checks and independent consequences keep the paper essentially non-circular.

fitted input called prediction [Section 3 (Adjustment of Xmax scale), paragraph discussing Fig. 1, p. 4]
"We then fit the energy evolution of ⟨Xmax⟩ above 10^19.6 eV using the prediction of these two models for iron nuclei allowing a shift of the Xmax scale: Xmax→Xmax+ΔXmax. We obtained ΔXmax= 52± 1+11−8 g/cm2 and ΔXmax = 29± 1+12−7 g/cm2 for QGSJet II-04 and Sibyll2.3d, respectively. ... After applying a shift of ΔXmax, the Xmax moments can be interpreted in a consistent manner, which is also in agreement with the results of the model-independent estimation of σ2(lnA)≈0.7−2.5 from the correlation between Xmax and ground signal in the energy range 3−10 EeV [13]."

The shift ΔXmax is the free parameter obtained by solving ⟨Xmax⟩_data ≈ ⟨Xmax⟩_Fe,model + ΔXmax above 40 EeV. Therefore, the statement that after applying ΔXmax the mean Xmax (and hence ⟨ln A⟩) 'can be interpreted in a consistent manner' as iron is the same equation re-read with the fitted value; it is forced by construction rather than an independent test. The independent content in this section is the pre-checked σ(Xmax) and the external agreements with the 3–10 EeV 2D fits [9] and σ2(ln A) from [13]; those keep the circularity minor.

full rationale

The paper transparently labels ΔXmax as a fitted quantity and the composition fractions as fits, so most of the derivation is honest parameter estimation rather than hidden prediction. Section 3's 'After applying a shift ... consistent manner' is the one place where a fitted parameter is presented as producing consistency: the agreement of ⟨Xmax⟩ with the shifted iron curve is guaranteed by the fit. However, the paper immediately cites independent support: the ΔXmax values are consistent with the 2D Xmax–ground-signal fits in 3–10 EeV [9], and σ2(ln A) agrees with model-independent estimates [13]; the σ(Xmax) consistency with pure iron was checked before the shift. Section 4 fits the Xmax distributions with shifted templates; the good tail description is a model-check, not a prediction. The spectral decomposition, muon-signal alleviation, and GMF backtracking are conditional consequences of the assumed scenario and are not claimed as first-principles predictions. The only author self-citation, [24], is an adopted method for dipole-direction restriction, not a load-bearing uniqueness theorem. Thus no central claim reduces to its inputs by construction beyond the one presentational tautology, and the overall circularity is minor.

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

The scenario rests on three explicit assumptions (Section 2), plus standard hadronic interaction models, GMF models, and detector corrections from the literature. The Xmax shift and the composition fractions are free parameters fitted to the data, and the 40 EeV threshold is a post hoc selection. No new particles, forces, or dimensions are introduced.

free parameters (4)

ΔXmax QGSJet II-04 = 52 +11/-8 g/cm2
Fitted to match the pure-iron prediction to the Auger mean Xmax above 40 EeV (Section 3).
ΔXmax Sibyll 2.3d = 29 +12/-7 g/cm2
Same fit as above for Sibyll 2.3d (Section 3).
Primary-fraction parametrizations = not fully specified
Gaussian times exponential functions fitted to the Auger Xmax distributions for p, He, N, and Fe (Section 4). The specific parameters are not listed.
Energy threshold 40 EeV = 10^19.6 eV
Chosen as the energy above which Xmax fluctuations match pure iron; this is a post hoc selection that is then used to derive the shift.

assumptions (6)

domain assumption Xmax fluctuations and the elongation rate from QGSJet II-04 and Sibyll 2.3d are correctly predicted.
Stated as the first of the three assumptions in Section 2; the paper provides no independent check of this for the shifted scenario.
ad hoc to paper All hadronic interaction models have a constant, energy-independent error in the absolute Xmax scale.
Assumption 2 in Section 2; the constant shift is fitted at high energy and applied at all lower energies without a physical mechanism.
ad hoc to paper The cosmic-ray beam above 40 EeV consists of pure iron nuclei.
Assumption 3 in Section 2; the main premise of the scenario, motivated by the Xmax fluctuation data.
domain assumption Template Xmax distributions for p, He, N, and Fe from simulations, corrected for detector effects according to ref [5], are valid for the fits.
Used in Section 4 for the composition fits and in the example distribution of Section 6.
domain assumption The total energy spectrum measured by Auger [14] can be decomposed into individual species using the fitted fractions.
Section 5 derives the species fluxes by multiplying the total spectrum by the fitted fractions; this assumes the fractions are unbiased.
domain assumption The Galactic magnetic field models UF23 and JF12 with Planck-tuned parameters are adequate for backtracking.
Section 7 uses these models for the source-direction and dipole studies; GMF uncertainties are not propagated.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Consequences of a Heavy-Metal Scenario of Ultra-High-Energy Cosmic Rays." pith.science (2026). https://pith.science/paper/ETDXOTLE

@misc{pith2026250208747,
  author       = {Pith},
  title        = {Pith review of: Consequences of a Heavy-Metal Scenario of Ultra-High-Energy Cosmic Rays},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ETDXOTLE}},
  note         = {Machine review of arXiv:2502.08747}
}

abstract

We assume an extreme scenario, in which the arriving cosmic rays are composed of only iron nuclei at energies above $10^{19.6}\,\text{eV}\simeq40\,\text{EeV}$, while allowing a freedom in the scale of the depth of shower maximum ($X_{\rm{max}}$) and preserving the elongation rate and fluctuations of $X_{\rm{max}}$ predicted by models of hadronic interactions. We derive the shift of the $X_{\rm{max}}$ scale for QGSJet II-04 and Sibyll 2.3d models using the public data from the Pierre Auger Observatory. We then propose a new mass-composition model for the energy evolution of four primary species at the ultra-high energies by fitting the publicly-available $X_{\rm{max}}$ distributions. We discuss the consequences of this Heavy-metal scenario on the energy spectrum of individual primary species, hadronic interaction studies, and the effect of the Galactic magnetic field on the arrival directions.

Figures

Figures reproduced from arXiv: 2502.08747 by the authors.

**Figure 2.** Top panels: The energy evolutions of four primary fractions fitted to the 𝑋max distributions from [5] together with curves corresponding to their parametrizations for QGSJet II-04 (left) and Sibyll 2.3d (right). Bottom panels: The differential flux measured at the Pierre Auger Observatory [14] and contributions of individual primaries derived from the parametrizations shown in the respective top panels. is consisten… view at source ↗

**Figure 3.** Top-left panel: The estimations of the muon signal measured directly at the Pierre Auger Observatory for three different energies [15, 16]. We show also the expectations for QGSJet II-04 between protons and iron nuclei (dashed lines), and after applying Δ𝑋max (full lines). Top-right panel: An example of the primary-fraction fit to the 𝑋max distribution measured by the Pierre Auger Observatory [5] using templates pre… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 14 canonical work pages

[9]

Pierre AugerCollaboration, Phys. Rev. D109 (2024) 102001 [2401.10740]

arXiv 2024
[13]

Pierre AugerCollaboration, Phys. Lett. B762 (2016) 288 [1609.08567]

arXiv 2016
[1]

Pierre AugerCollaboration, Nucl. Instrum. Meth. A798 (2015) 172 [1502.01323]

arXiv 2015
[2]

Pierog, I

T. Pierog, I. Karpenko, J.M. Katzy et al.,Phys. Rev. C92(2015) 034906

work page 2015
[3]

Ostapchenko,Phys

S. Ostapchenko,Phys. Rev. D83(2011) 014018

work page 2011
[4]

Riehn, R

F. Riehn, R. Engel, A. Fedynitch et al.,Phys. Rev. D102 (2020) 063002

work page 2020
[5]

Pierre AugerCollaboration, Phys. Rev. D90(2014) 122005

work page 2014
[6]

Pierre AugerCollaboration, PoS ICRC2023(2023) 319

work page 2023

Show all 24 references

[7]

Pierog and K

T. Pierog and K. Werner,PoS ICRC2023(2023) 230

2023
[8]

Ebr et al.,PoS ICRC2023(2023) 245

J. Ebr et al.,PoS ICRC2023(2023) 245

2023
[10]

Pierre AugerCollaboration, Phys. Rev. D96(2017) 122003

2017
[11]

Pierre AugerCollaboration, Phys. Rev. D111 (2025) 022003 [2406.06319]

2025 arXiv
[12]

Pierre AugerCollaboration, JCAP 02 (2013) 026

2013
[14]

Pierre AugerCollaboration, Phys. Rev. D102 (2020) 062005

2020
[15]

Pierre AugerCollaboration, Eur. Phys. J. C210 (2020) 751

2020
[16]

Pierre AugerCollaboration, Phys. Rev. Lett.126 (2021) 152002

2021
[17]

Pierre AugerCollaboration, Science 357 (2017) 1266 [1709.07321]

2017 arXiv
[18]

Pierre AugerCollaboration, Astrophys. J. Suppl.264 (2023) 50 [2211.16020]

2023 arXiv
[19]

Alves Batista, J

R. Alves Batista, J. Becker Tjus, J. Dörner et al.,JCAP 2022 (2022) 035 [2208.00107]

2022 arXiv
[20]

Unger and G.R

M. Unger and G.R. Farrar,The Astrophysical Journal970 (2024) 95 [2311.12120]

2024 arXiv
[21]

Jansson and G.R

R. Jansson and G.R. Farrar,The Astrophysical Journal757 (2012) 14 [1204.3662]

2012 arXiv
[22]

PlanckCollaboration, Astronomy and Astrophysics596 (2016) A103 [1601.00546]

2016 arXiv
[23]

Sommers,Astroparticle Physics14(2001) 271

P. Sommers,Astroparticle Physics14(2001) 271

2001
[24]

Bakalová, J

A. Bakalová, J. Vícha and P. Trávníček,JCAP 12(2023) 016 [2303.08766]. 8

2023 arXiv

Pith tools

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

[1] [9]

Pierre AugerCollaboration, Phys. Rev. D109 (2024) 102001 [2401.10740]

arXiv 2024

[2] [13]

Pierre AugerCollaboration, Phys. Lett. B762 (2016) 288 [1609.08567]

arXiv 2016

[3] [1]

Pierre AugerCollaboration, Nucl. Instrum. Meth. A798 (2015) 172 [1502.01323]

arXiv 2015

[4] [2]

Pierog, I

T. Pierog, I. Karpenko, J.M. Katzy et al.,Phys. Rev. C92(2015) 034906

work page 2015

[5] [3]

Ostapchenko,Phys

S. Ostapchenko,Phys. Rev. D83(2011) 014018

work page 2011

[6] [4]

Riehn, R

F. Riehn, R. Engel, A. Fedynitch et al.,Phys. Rev. D102 (2020) 063002

work page 2020

[7] [5]

Pierre AugerCollaboration, Phys. Rev. D90(2014) 122005

work page 2014

[8] [6]

Pierre AugerCollaboration, PoS ICRC2023(2023) 319

work page 2023