Pith. sign in

REVIEW 3 major objections 4 minor 14 references

Probing the high energy spectrum of neutral pions in ultra-high energy proton-Air interactions

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

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

Pith's one-line read The slope of the low muon-number tail at the ground is a direct probe of the high-energy neutral-pion spectrum of the first proton-air collision.

desk verdict Competent proceedings paper with a useful calibration result, but the abstract's 'direct link' to the neutral-pion spectrum overstates what the simulation actually shows. read the letter →

arxiv 1908.09668 v1 pith:XNJRRR7L submitted 2019-08-26 hep-ph hep-ex

Lorenzo Cazon , Ruben Conceição , Miguel Martins , Felix Riehn This is my paper · ORCID
classification hep-phhep-ex
keywords ultra-highenergycosmicraysextensiveairshowersmuonnumberfluctuationsneutralpionproductionfirsthadronicinteractionmodelsproblemforwardLHCphysics
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

Ultra-high-energy cosmic-ray air showers hide their first collision: by the time a proton hits an air nucleus at center-of-mass energy around 100 TeV, the details of that interaction are buried under tens of generations of secondary particles. This paper argues that one fragment of the ground signal escapes that burial. Using simulations with three post-LHC hadronic models, it shows that the exponential low-value tail of the shower-to-shower distribution of muon number, $N_\mu$, is set by the energy flowing into the hadronic sector in the very first interaction, not by the later cascade. The paper then connects that first-interaction energy split to the inclusive production of high-energy neutral pions, and shows that measuring the tail slope $\Lambda_\mu$ at a ground observatory would constrain the neutral-pion spectrum at $\sqrt{s}\sim100$ TeV, beyond accelerator reach. The practical payoff is a ground-based window into multiparticle production at energies no collider currently accesses.

What carries the argument

The paper's central object is the exponential low tail of the ground muon-number distribution for proton-induced showers, quantified by its slope $\Lambda_\mu$. That slope is interpreted through a weighted first-interaction energy variable $\alpha_1 = \sum_i (E_{{\rm had},i}/E_0)^\beta$, with $\beta = \log(m)/\log(m_{\rm tot})$ from the Heitler-Matthews cascade model; this variable encodes how the first interaction splits energy between the hadronic and electromagnetic sectors. A calibration curve between $\Lambda_\mu$ and $\Lambda_{\rm had}$, obtained by reweighting simulated showers, is what lets a ground measurement stand in for the first-interaction energy flow. The neutral-pion spectrum enters through $E_{\rm had} = 1 - E_{\rm em}$, since $\pi^0 \to \gamma\gamma$ feeds the electromagnetic component and the highest-energy pions carry the information about fast leading particles.

What would settle it

If, on real data with good statistics, the low-$N_\mu$ tail slope changed with atmospheric depth or with the amount of later shower development while the first-interaction energy flow were held fixed, the claimed direct link would fail; a direct test would be to compare the $\Lambda_{\rm had}$ inferred from $\Lambda_\mu$ at approximately $10^{17}$ eV with the same quantity measured by forward LHC experiments at 13 TeV.

Watch

Extended reading notes

Core claim

The central claim is that the shape of the shower-to-shower muon-number distribution, specifically its exponential tail at low $N_\mu$, is controlled by the first hadronic interaction of the primary proton. The paper defines $\Lambda_\mu$ as the slope of this tail and shows, by reweighting simulated showers from a large ensemble, that $\Lambda_\mu$ is tied to $\Lambda_{\rm had}$, the slope of the distribution of the hadronic energy fraction $E_{\rm had}/E_0$ deposited by that first interaction. Because the electromagnetic sector is fed almost entirely by neutral pions, fluctuations of $E_{\rm had}/E_0$ are the same as fluctuations of the electromagnetic energy fraction, so $\Lambda_{\rm had}$ in turn reflects the inclusive high-energy tail of the neutral-pion spectrum. Simulation with QGSJET-II.04, EPOS-LHC, and SIBYLL 2.3c shows that modifying the neutral-pion inclusive cross-section at large $x_L$ moves $\Lambda_\mu$, and that a mixed composition of 25% protons, 50% helium, and 25% nitrogen still permits a clean measurement of the proton tail when the muon number is smeared by 20%.

Load-bearing premise

The link between the measured low-muon tail slope and the first-interaction hadronic-energy slope is universal across hadronic models and remains valid in real mixed-composition data; the paper demonstrates it only by reweighting simulated showers, not by analytic derivation.

Editorial extensions

If this is right

  • A ground array that collects on the order of 3000 proton-like showers could distinguish among the tested hadronic interaction models using the low-$N_\mu$ tail slope alone.
  • The same measurement at cosmic-ray energies around $10^{17}$ eV could be confronted with LHC forward measurements at $\sqrt{s}=13$ TeV, offering a direct accelerator-to-cosmic-ray cross-check.
  • If the calibration holds, $\Lambda_\mu$ becomes a measurement of the fluctuation of the hadronic energy fraction in the first interaction at $\sqrt{s}\sim100$ TeV, where no collider data exist.
  • A suppression of neutral-pion production at large $x_L$ shifts the muon tail, so the tail slope provides a test of high-rapidity multiparticle production and of possible violations of longitudinal scaling at ultra-high energies.
  • The measurement remains feasible even under a pessimistic mixed composition with abundant helium, provided enough showers are recorded.

Reading between the lines

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

  • A natural extension the authors do not spell out: the same tail analysis could be applied to the electromagnetic component, such as shower-maximum depth or ground electrons, to separate first-interaction energy flow from subsequent shower physics and cross-check the muon-based calibration.
  • If the model-dependent step between $\Lambda_{\rm had}$ and the neutral-pion spectrum cannot be sharpened, the practical reach of the method may be limited to detecting qualitative deviations from standard-model expectations, such as an anomalous leading-particle energy fraction, rather than a precise spectrum.
  • The technique effectively converts the atmosphere into a forward detector for $\sqrt{s}\sim100$ TeV proton-air collisions; a dedicated measurement of $\Lambda_\mu$ as a function of primary energy could map the onset of any new high-energy behavior.
  • An obvious testable extension is to apply the tail-slope method to existing high-statistics data sets and compare the derived $\Lambda_{\rm had}$ at overlapping energies with forward LHC measurements; any mismatch would pinpoint either the calibration or new physics.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. This proceedings paper claims that the slope of the low-number tail of the muon-number distribution, Λ_mu, measured at the ground in ultra-high-energy cosmic-ray air showers is a direct probe of properties of the first hadronic interaction, specifically the fraction of primary energy transferred to the hadronic component and, more strongly, the high-energy tail of the neutral-pion energy spectrum. The study uses CONEX simulations of proton-induced showers at 10^19 eV and 67° zenith angle with a 1 GeV muon threshold. Section 3.1 establishes a simulation-based calibration between Λ_mu and the slope Λ_had of the first-interaction hadronic-energy fraction distribution, and shows that Λ_mu can be extracted in a mixed-composition scenario with 20% muon-number smearing. Section 3.2 reports that suppressing high-x_L neutral-pion production in SIBYLL 2.3c changes the N_mu tail, and the authors argue that this connects Λ_mu to the pion spectrum. The paper also discusses experimental precision and connections to LHC forward measurements.

Significance. If the claimed connection holds, this observable would provide a new way to access first-interaction physics at center-of-mass energies around 100 TeV, beyond current accelerators, using existing or planned cosmic-ray observatories such as the Pierre Auger Observatory. The paper has notable strengths: the simulation setup is clearly stated, statistical bands are shown, the mixed-composition and detector-smearing analysis is a useful feasibility check, and the comparison across three post-LHC hadronic models is informative. The central limitation, however, is that the evidence supports sensitivity to the first-interaction hadronic/electromagnetic energy split, not yet a model-independent sensitivity to the neutral-pion spectral shape. The abstract's phrase 'direct link' overstates what the presented simulations establish.

major comments (3)
  1. [Sec. 3.2, Fig. 5] The simulation that suppresses high-x_L neutral-pion production changes not only the spectral shape of the π0 distribution but also the total electromagnetic energy fraction E_em = 1 - E_had, and hence α1 in Eq. (2.1). The observed change in the N_mu tail is therefore fully compatible with Λ_mu tracking only the hadronic/electromagnetic energy split, with no demonstrated sensitivity to the shape of the π0 spectrum beyond its contribution to that split. To support the central claim, the authors should show a control test where the π0 spectral shape is varied while E_had/E0 (or α1) is held fixed, or provide a quantitative decomposition showing that Λ_mu carries information about the shape independently of the first moment.
  2. [Sec. 3.1, Fig. 3 (left)] The calibration between Λ_mu and Λ_had is obtained by selecting simulated showers from a large ensemble and fitting the response, rather than from a physical derivation or an out-of-sample test. The claim that the relation is independent of hadronic models is supported only by three post-LHC models that share common assumptions about the leading-particle and energy-flow behavior. The figure shows model-dependent lines, so the 'independently of the hadronic interaction models' assertion requires a stronger demonstration, for example a model with a deliberately different first-interaction energy-flow distribution or a cross-validation on independent simulation sets.
  3. [Abstract vs. Sec. 3.2, final paragraph] The abstract claims that the slope of the low-N_mu tail is 'a direct link to the high energy spectrum of neutral pions,' but the text explicitly states that the derived function relating Λ_mu to the π0 spectrum 'has some dependence on the details of the hadronic interaction models.' This is an internal inconsistency in the level of claim. Either the abstract should be softened to reflect the model-dependent calibration, or the paper should provide an estimate of the systematic uncertainty and demonstrate that the model dependence does not affect the qualitative conclusion.
minor comments (4)
  1. [Sec. 4] There is an encoding artifact in the text: 'protons.¢aMoreover' should read 'protons. Moreover'. Also, 'the slope of theNµ' is missing a space.
  2. [Fig. 3 (left), caption] The caption reads 'Conversion between Λµ and Λα,' but the text and axes use Λ_had; please make the notation consistent.
  3. [Eq. (2.2)] The symbols m and m_tot are not defined precisely enough; it is unclear whether they refer to charged multiplicity, total hadronic multiplicity, or something else. Please clarify in the text.
  4. [Introduction and Conclusions] The paper uses the words 'prove' and 'proved' for statements obtained from a finite set of simulations; 'demonstrate' or 'show' would be more accurate and appropriate for a proceedings contribution.

Circularity Check

1 steps flagged · score 4.0 of 10

The central premise that the N_mu tail shape is governed by first-interaction hadronic energy flow is load-bearing and imported from the authors' own prior paper [8]; the new calibration and pion-spectrum simulations add independent content but do not break the self-citation chain.

  1. self citation load bearing [Section 2, paragraph after Eq. (2.2): 'It is shown in [8] that α1 distribution...']
    "It is shown in [8] that α1 distribution of the first interaction is enough to describe the main features (the width1 and exponential tail) while the rest of the shower contributes only to the overall value of Nµ."

    The paper's derivation chain starts from the assertion that the shape of the N_mu distribution, especially its exponential tail, is determined by the first-interaction hadronic-energy flow. That foundational assertion is not re-derived here; it is lifted verbatim from [8], a prior paper by the same authors (Cazon, Conceição, Riehn). The subsequent argument uses this premise as established: Sec. 3.1 identifies the measured slope Lambda_mu with Lambda_had, and Sec. 3.2 extends this to the neutral-pion spectrum. Because the load-bearing premise is a self-citation rather than an independently validated, externally benchmarked result, the first link of the claimed 'direct link' reduces to the authors' own earlier work. The new simulation calibrations in Sec.

full rationale

No equation-level tautology was found: alpha1 is defined from first-interaction hadronic energies, and Lambda_had is defined as the slope of the Ehad/E0 distribution; the mapping from Lambda_mu to Lambda_had is a simulation-based calibration, not a definitional identity. Similarly, the pion-spectrum sensitivity shown in Sec. 3.2 is a Monte Carlo check, not a fitted parameter renamed as a prediction. The main circularity-adjacent issue is the load-bearing self-citation to [8]: the paper's starting point that N_mu fluctuations and the exponential tail are governed by the first interaction is imported from the same group's earlier work. This is not a fully circular derivation because the calibration curves and the pion-spectrum simulation are new and falsifiable in principle, and the paper even admits that the Lambda_mu-to-pion-spectrum function retains model dependence. That admission weakens the abstract's 'direct link' claim, but it is an overstatement/correctness concern rather than a definitional circularity. Overall, the central claim retains independent simulation content, but the foundational premise is self-cited, so a moderate score of 4 is appropriate.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The quantitative framework rests on simulation inputs and on model-dependent choices: beta in alpha1 comes from Heitler-Matthews with model multiplicities; the calibration between Lambda_mu and Lambda_had is empirical; the 20% smearing is a hand-chosen detector assumption. The paper introduces no new physical entities, so invented_entities is empty.

free parameters (3)
  • Beta (power-law index in alpha1) = Not stated numerically; model-dependent via first-interaction multiplicities
    Defines alpha1 = sum_i (Ehad_i/E0)^beta, used for correlation with N_mu. Beta = log(m)/log(mtot) from Eq. (2.2) relies on multiplicities from hadronic models, not on an independent measurement.
  • N_mu tail fit range = Not specified; indicated by vertical bars in Fig. 3 (right)
    The measured slope Lambda_mu depends on the chosen fitting interval, which is not given in the text; this is a hand-chosen analysis choice.
  • Detector N_mu smearing = 20% (Gaussian smearing assumed)
    Chosen conservatively to mimic the main experimental uncertainty; not derived from a specific detector simulation, and the exact smearing procedure is not described.
assumptions (5)
  • domain assumption The shape of the N_mu distribution, especially its low tail, is essentially determined by the first interaction, while later interactions only affect the average.
    Invoked in Sec. 2 to justify reading first-interaction properties from the tail; supported only by simulation correlation, not by an analytic proof.
  • domain assumption Ehad/E0 is strongly correlated with N_mu in the low tail, and the tail of the N_mu distribution is dominated by Ehad/E0.
    Sec. 3.1 ansatz; no closed-form relation is provided.
  • ad hoc to paper A calibration curve between Lambda_mu and Lambda_had derived from simulations is independent of the hadronic interaction model.
    The model independence is an empirical claim from Fig. 3 (left) for three models, not a derived result; the paper itself notes pion-level calibration has model dependence.
  • domain assumption The Heitler-Matthews relation beta = log(m)/log(mtot) describes the energy/multiplicity scaling of muon production.
    Used in Eq. (2.2) to build alpha1; this is a simplified shower model, not a measured law.
  • domain assumption CONEX and the three hadronic interaction models adequately describe air-shower development and muon transport for this purpose.
    All results are simulation-based; no comparison to data is shown.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Probing the high energy spectrum of neutral pions in ultra-high energy proton-Air interactions." pith.science (2026). https://pith.science/paper/XNJRRR7L

@misc{pith2026190809668,
  author       = {Pith},
  title        = {Pith review of: Probing the high energy spectrum of neutral pions in ultra-high energy proton-Air interactions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XNJRRR7L}},
  note         = {Machine review of arXiv:1908.09668}
}
read the original abstract

The interaction of ultra-high energy cosmic rays with the atmosphere nuclei has long been seen as a unique opportunity to study hadronic interactions above energies attainable by accelerators. However, so far the multiparticle production properties of the first interaction have been difficult to assess as they are masked by the many interactions that outline the shower development. In this work, we demonstrate, that relevant properties of the ultra-high energy first interaction can be accessed through the analysis of the shower-to-shower distribution of muons arriving at the ground. In particular, it is shown that the slope of the low-tail of the number of muon distribution measured at the ground is a direct link to the high energy spectrum of neutral pions produced in the interaction of the primary protons. In this presentation, it will also address the experimental feasibility of such measurements and their connection with physical quantities being currently measured at the Large Hadron Collider.

Figures

Figures reproduced from arXiv: 1908.09668 by the authors.

Figure 1
Figure 1. (left) Distribution of the number of muons at the ground in simulated EAS induced by protons with E = 1019 eV and θ = 67◦ . The curve were generated with CONEX and so the muon energy threshold is 1GeV. (right) Scheme of the energy flow in the first interaction into the electromagnetic and hadronic component of the shower. α1 = m ∑ i=1  E had i E0 β (2.1) where E had i /E0 is the fraction of energy carried by each … view at source ↗
Figure 2
Figure 2. (left) Distribution of the hadronic energy, Ehad/E0, and Nµ and (right) distribution of α1 versus Nµ . Both plots were generated with showers induced by protons having an energy of E = 1019 eV and a zenith angle of θ = 67◦ . Let us now consider that in the first interaction, only a small amount of energy was passed into the hadronic component. In this case, the size of the hadronic shower would be smaller, regardles… view at source ↗
Figure 3
Figure 3. (left) Conversion between Λµ and Λα. The filled circles indicate the predictions by the different interaction models. The lines show how Λµ changes for each model if Λhad is changed. (right) Number of showers with a given number of muons at ground for the nuclear primaries: p, He and N. The vertical bars indicate the fitting range.The total number of events is ≈ 200 k. The measurement of Λhad has two major requireme… view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: (left) Relative fluctuation of the reconstructed Λµ as a function of the number of events in the experiments for the mass composition scenario 1:2:1:0. The x-bottom scale shows the number of events that fall in the fit region while the axis on top shows the total numbe…
Figure 5
Figure 5. Figure 5: (left) Energy spectrum of neutral pions as a function of the lab. energy fraction, xL. The nominal cross-section in Sibyll 2.3c is shown in gray. The yellow curve represents a modified cross-section where production of neutral pions at large xL is suppressed. (right) D…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 5 canonical work pages

  1. [1]

    A model for the transport of muons in extensive air showers

    L. Cazon, R. Conceição, M. Pimenta and E. Santos, A model for the transport of muons in extensive air showers, Astropart. Phys. 36 (2012) 211–223, [1201.5294]

    work page Pith review arXiv 2012

  2. [2]

    Andringa, L

    S. Andringa, L. Cazon, R. Conceição and M. Pimenta, The Muonic longitudinal shower profiles at production, Astropart. Phys. 35 (2012) 821–827, [1111.1424]

    arXiv 2012

  3. [3]

    Aab et al., Muons in air showers at the Pierre Auger Observatory: Mean number in highly inclined events, Phys

    P IERRE AUGER collaboration, A. Aab et al., Muons in air showers at the Pierre Auger Observatory: Mean number in highly inclined events, Phys. Rev. D91 (2015) 032003, [1408.1421]

    arXiv 2015

  4. [4]

    T ELESCOPE ARRAY collaboration, R. U. Abbasi et al., Study of muons from ultrahigh energy cosmic ray air showers measured with the Telescope Array experiment, Phys. Rev. D98 (2018) 022002, [1804.03877]

    work page Pith review arXiv 2018

  5. [5]

    Ostapchenko, Monte Carlo treatment of hadronic interactions in enhanced Pomeron scheme: I

    S. Ostapchenko, Monte Carlo treatment of hadronic interactions in enhanced Pomeron scheme: I. QGSJET-II model, Phys. Rev. D83 (2011) 014018, [1010.1869]

    arXiv 2011

  6. [6]

    Pierog, I

    T. Pierog, I. Karpenko, J. M. Katzy, E. Yatsenko and K. Werner, EPOS LHC: Test of collective hadronization with data measured at the CERN Large Hadron Collider, Phys. Rev. C92 (2015) 034906, [1306.0121]

    arXiv 2015

  7. [7]

    Riehn, H

    F. Riehn, H. P. Dembinski, R. Engel, A. Fedynitch, T. K. Gaisser and T. Stanev,The hadronic interaction model SIBYLL 2.3c and Feynman scaling, PoS ICRC2017 (2017) 301, [1709.07227]

    arXiv 2017

  8. [8]

    Probing the energy spectrum of hadrons in proton air interactions at ultrahigh energies through the fluctuations of the muon content of extensive air showers

    L. Cazon, R. Conceição and F. Riehn, Probing the energy spectrum of hadrons in proton air interactions at ultrahigh energies through the fluctuations of the muon content of extensive air showers, Phys. Lett. B784 (2018) 68–76, [1803.05699]

    work page Pith review arXiv 2018

Show all 14 references
  1. [9]

    Matthews, A Heitler model of extensive air showers, Astropart

    J. Matthews, A Heitler model of extensive air showers, Astropart. Phys. 22 (2005) 387–397

    2005

  2. [10]

    Abreu et al., Measurement of the proton-air cross-section at√s = 57 TeV with the Pierre Auger Observatory, Phys

    P IERRE AUGER collaboration, P. Abreu et al., Measurement of the proton-air cross-section at√s = 57 TeV with the Pierre Auger Observatory, Phys. Rev. Lett. 109 (2012) 062002, [1208.1520]

    2012 arXiv

  3. [11]

    Aab et al., Depth of maximum of air-shower profiles at the Pierre Auger Observatory

    P IERRE AUGER collaboration, A. Aab et al., Depth of maximum of air-shower profiles at the Pierre Auger Observatory. II. Composition implications, Phys. Rev. D90 (2014) 122006, [1409.5083]

    2014 arXiv

  4. [12]

    Chatrchyan et al., Measurement of pseudorapidity distributions of charged particles in proton-proton collisions at√s = 8 TeV by the CMS and TOTEM experiments, Eur

    CMS, TOTEM collaboration, S. Chatrchyan et al., Measurement of pseudorapidity distributions of charged particles in proton-proton collisions at√s = 8 TeV by the CMS and TOTEM experiments, Eur. Phys. J.C74 (2014) 3053, [1405.0722]

    2014 arXiv

  5. [13]

    Adriani et al., Measurements of longitudinal and transverse momentum distributions for neutral pions in the forward-rapidity region with the LHCf detector, Phys

    LHC F collaboration, O. Adriani et al., Measurements of longitudinal and transverse momentum distributions for neutral pions in the forward-rapidity region with the LHCf detector, Phys. Rev. D94 (2016) 032007, [1507.08764]

    2016 arXiv

  6. [14]

    Aab et al., Prototype muon detectors for the AMIGA component of the Pierre Auger Observatory, JINST 11 (2016) P02012, [1605.01625]

    P IERRE AUGER collaboration, A. Aab et al., Prototype muon detectors for the AMIGA component of the Pierre Auger Observatory, JINST 11 (2016) P02012, [1605.01625]. 7

    2016 arXiv

Pith tools

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