Pith. sign in

REVIEW 2 major objections 6 minor 15 references

A New Method for Measuring the Pion-Air Cross Section at Multi-TeV Energies Using Muon Bundle Properties in Deep Underground Detectors

T0 review · 2 major / 6 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 This paper argues that underground muon-bundle counts can measure the inelastic pion-air cross section at TeV energies, with sensitivity projected at 10-15% around 1-2 TeV.

desk verdict A clean, genuinely new sensitivity forecast for pion-air cross sections from underground muon bundles, but the projected 10-15% is for a mean-multiplicity observable that ignores shower-to-shower fluctuations. read the letter →

arxiv 2507.22463 v3 pith:ULTDYH3B submitted 2025-07-30 astro-ph.HE

Karolin Hymon , Anatoli Fedynitch This is my paper
classification astro-ph.HE
keywords pion-aircrosssectionmuonbundlesundergrounddetectorsIceCubeKM3NeTpuzzlehadronicinteractionmodelscosmic-rayairshowers
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 that the muon bundles recorded by deep underground detectors such as IceCube carry a measurable imprint of the inelastic pion-air cross section. Because accelerator beams cannot reach multi-TeV pion-air collisions, the authors develop a sensitivity forecast: after replacing the SIBYLL 2.3c cross section with a spline and perturbing it, they compute how the muon multiplicity distribution responds. With a realistic assumption on multiplicity resolution (an uncertainty of log10(ΔM/M)=0.1 per bin), they claim the pion-air inelastic cross section could be constrained to roughly 10-15% between 1 and 2 TeV, and to better than 50% between 200 GeV and 6 TeV at 1.5 km w.e. depth. If real IceCube or KM3NeT data are analyzed this way, it would provide the first multi-TeV constraint on this quantity and reduce the hadronic-model uncertainty behind the muon puzzle.

What carries the argument

The central object is the underground muon multiplicity distribution M(N_mu), computed by MCEq and MUTE as a function of cosmic-ray energy and depth. The argument is carried by a cubic-spline parameterization of the inelastic cross section, denoted \hat{\$\sigma$}_{c_i}, with 10 knots spaced logarithmically from 10 GeV to $10^{10}$ GeV; the 12 spline coefficients c_i are the tuning parameters. Finite-difference gradients \partial M/\partial \hat{c}_i (Eq. 4) form a linearized response model (Eq. 5), which is fitted to pseudo-data via a penalized chi-square (Eq. 6) with Gaussian priors; the covariance of the fitted coefficients is propagated through the cross-section Jacobian (Eq. 7). The physics mechanism is the competition between interaction and decay lengths near the critical energy (~100-200 GeV), which makes an increased cross section raise the multiplicity through low-energy muon pile-up.

What would settle it

Generate mock underground muon multiplicity data from a simulation in which the pion-air cross section is artificially increased by, say, 20%, then run the paper's spline-fitting pipeline on it: if the 20% increase is not recovered within the quoted 10-15% uncertainty, the method's projected sensitivity is overstated.

Watch

Extended reading notes

Core claim

The paper establishes, through a numerical sensitivity study, that the inclusive underground muon multiplicity distribution changes in a specific, calculable way when the inelastic pion-air cross section is increased: at energies above the ~100-200 GeV critical energy, a larger cross section suppresses higher-energy muons and boosts lower-energy pile-up, raising the multiplicity. Using MCEq with SIBYLL 2.3c at the South Pole location, the authors replace the nominal cross section with a cubic spline having 10 logarithmic knots, compute the gradient of the flux with respect to each knot (Eq. 4), and build a linearized model (Eq. 5). Pseudo-data generated at the nominal model with log10(ΔM/M)=0.1 per bin, a penalized chi-square fit with Gaussian priors, and covariance propagation yield the projected uncertainties: 10-15% for sigma_{\pi+air} around 1-2 TeV and below 50% from 200 GeV to 6 TeV at 1.5 km w.e., with a second sensitivity window near 200 TeV; sigma_{p+air} is constrained to ~10% around 80-200 TeV. The energy and depth dependence cleanly separates the pion and proton contributions.

Load-bearing premise

The method assumes that SIBYLL 2.3c's treatment of everything except the pion-air cross section is accurate, so perturbing just that cross section in the simulation reproduces what a real change would do to underground muon bundles.

Editorial extensions

If this is right

  • Real IceCube or KM3NeT bundle counts could yield the first constraints on the multi-TeV inelastic pion-air cross section.
  • If the quoted sensitivity holds, hadronic interaction models would gain a new anchor point at energies no accelerator can currently reach, directly feeding into air-shower simulations.
  • The method also constrains the proton-air cross section around 100 TeV, complementing LHC-based measurements.
  • Extending the analysis to multiple zenith angles could broaden the measurable energy window from a few TeV to several hundred TeV.
  • A measured cross section that disagrees with SIBYLL 2.3c would give concrete input for resolving the muon puzzle rather than only parameterizing it.

Reading between the lines

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

  • If real data return constraints near the projected precision, the same spline-perturbation technique could be applied to kaon-air cross sections or to other hadronic observables such as charm production, extending the method beyond pions.
  • The depth dependence the authors find (peak sensitivity shifting ~1 TeV per km w.e.) suggests that a single detector measuring bundles at multiple depths, e.g., via different overburdens or zenith angles, could map the cross section as a function of energy with one dataset.
  • The projection assumes SIBYLL 2.3c is otherwise correct; a comparison against predictions with other hadronic models would reveal how much of the claimed reach is model-dependent rather than physics.
  • Cross-section constraints from this method could also serve as a prior for atmospheric lepton flux calculations, reducing the systematic error in neutrino-telescope analyses.
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 / 6 minor

Summary. This paper proposes a new indirect method to constrain the inelastic pion-air cross section and proton-air cross section at multi-TeV energies using the muon multiplicity distribution of bundles detected in deep underground Cherenkov detectors such as IceCube and KM3NeT. The authors use MCEq and MUTE with the SIBYLL 2.3c hadronic interaction model to compute mean muon multiplicities as a function of primary cosmic-ray energy, replace the nominal cross section with a cubic-spline parametrization, and estimate the sensitivity of the bundle multiplicity distribution to each spline coefficient via a penalized chi-square Fisher forecast. They find that, assuming a per-bin multiplicity uncertainty of log10(Delta M/M) = 0.1, the pion-air cross section can be constrained to roughly 10-15% between about 0.7 and 2 TeV at a depth of 1.5 km w.e., with a second region of 25-30% sensitivity around 200 TeV, while the proton-air cross section can be constrained to about 10% near 100 TeV. The paper concludes that this method could provide the first multi-TeV pion-air cross-section constraints and help address the muon puzzle.

Significance. If the forecast is validated, the method would fill an important experimental gap: accelerator measurements of pion-air interactions currently stop at a few hundred GeV, so a multi-TeV constraint from cosmic-ray data would be genuinely novel and useful for air-shower simulations. The study is computationally concrete: it uses established tools (MCEq, MUTE, SIBYLL 2.3c), parameterizes the cross section with a spline, and propagates gradients and covariances in a standard Fisher framework. The depth-dependent sensitivity map is falsifiable and could guide future analyses with real IceCube or KM3NeT data. However, the numerical claims rest on two unverified statistical assumptions: the muon multiplicity is treated as a deterministic function of primary energy, and the per-bin uncertainty is specified without defining the binning or exposure. These issues make the current 10-15% figure an optimistic estimate rather than a robust projection.

major comments (2)
  1. [Sec. 2, Eq. (3)] Eq. (3) defines the bundle flux M(N_mu) as a deterministic change of variables from E_CR to N_mu, using the mean multiplicity <N_mu(E_CR)> computed by MCEq/MUTE. Real underground detectors observe event-by-event muon multiplicities, which fluctuate around <N_mu> because of shower-to-shower variations in the first interaction, inelasticity, and particle production. The observable is therefore a convolution of Phi_CR(E_CR) with the conditional distribution P(N_mu | E_CR), not the right-hand side of Eq. (3). Because both the pseudo-data and the model in Secs. 2.2 and 4 are generated from the same deterministic map, the Fisher forecast in Eqs. (6)-(7) is for a different observable unless the fluctuations are negligible. Please quantify the width of P(N_mu | E_CR) at the energies of interest, using e.g. full air-shower simulations, and either show that it is small relative to the assumed log10(Delta M/M)=0.1 per-bin uncertainty or propagate it through Eq. (3). This is a concrete and testable issue that directly affects the headline 10-15% uncertainty.
  2. [Sec. 4, Eqs. (6)-(7)] The chi-square in Eq. (6) assumes a per-bin uncertainty log10(Delta M/M)=0.1, but the number of multiplicity bins, their widths, and whether they are treated as independent are not specified. For a Fisher forecast, the parameter covariance scales with the number of independent bins and with the bin width, so the absolute uncertainties shown in Fig. 4 and quoted in the text (10-15% for sigma_pi+air) are not well-defined without this information. Please specify the binning scheme used to evaluate Eq. (6), and if the result is intended to be binning-independent, demonstrate that explicitly. If the 0.1 dex uncertainty includes statistical contributions, the assumed exposure and observation time should also be stated. Without these details, the numerical precision claims cannot be reproduced or independently assessed.
minor comments (6)
  1. [Sec. 2, Eq. (2)] The summation symbol in Eq. (2) is not defined; specify whether the sum runs over primary species, nucleons, or perturbed runs, and clarify the resulting units of the 'muon multiplicity flux'.
  2. [Sec. 2.2, Eq. (4)] The notation M((1+delta)c_i) is ambiguous because c_i are spline coefficients rather than direct cross-section values; clarify how a relative perturbation of a coefficient translates into the modified cross-section table that is inserted into MCEq.
  3. [Sec. 3, Fig. 3] The left panels of Fig. 3 plot the ratio of log10(modified) to log10(nominal), which compresses the visual magnitude of the effect; consider showing the relative difference (modified/nominal - 1) in linear scale for easier interpretation.
  4. [Sec. 4, Fig. 4] The two panels use different x-axis labels (E and E_p), and the 'ratio' subpanels are not defined in the caption; state what quantity is plotted in each subpanel and what the SIBYLL 2.3c curve represents.
  5. [Sec. 2.2] After mentioning 10 knots, the text says the spline is 'fully defined by 12 coefficients'; a one-sentence explanation of the two edge-padding coefficients would remove confusion.
  6. [References] Reference [15] is marked as 'Submitted to Phys. Rev. D (2025)'; if it is not yet published, this should be stated explicitly in the citation or updated once the published version is available.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the paper is a self-consistent Fisher sensitivity forecast that uses its own model as pseudo-data and claims only projected uncertainties, not an empirical cross-section measurement.

full rationale

The derivation chain is transparent: MCEq with SIBYLL 2.3c computes muon yields; Eq. (3) transforms them into a multiplicity distribution via a Jacobian; the pion-air and proton-air cross sections are reparameterized as a spline with coefficients c_i; central finite differences give the gradients dM/dc_i; pseudo-data are set to the nominal distribution M_nominal; a chi-square with log10(Delta M/M)=0.1 yields a covariance; Eq. (7) propagates that covariance to sigma(E). No step uses an external measurement of sigma_pi+air as an input, and no output is presented as an empirical measurement. The pseudo-data are generated from the same model used to compute the gradients, but this is the standard construction of a sensitivity or Fisher forecast: the paper explicitly states that 'the inferred uncertainties on the cross sections represent the anticipated sensitivity achievable given the assumed experimental precision in measuring muon multiplicity distributions.' The self-citations to MCEq, MUTE, and SIBYLL provide the simulation machinery; these are code-based, externally used tools whose predictions are not defined in terms of the target cross-section measurement, so they are independent support rather than a circular premise. The main caveats, such as SIBYLL fidelity and the deterministic treatment of the N_mu-E_CR relation in Eq. (3), are model-validity and statistical-modeling risks, not circularity. Therefore no circular step can be exhibited.

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

The central claim rests on three types of input: the spline coefficients and prior widths that drive the sensitivity forecast, the SIBYLL 2.3c model that generates both pseudo-data and gradients, and the assumed statistical uncertainty on measured muon multiplicities. No new particles, forces, or physical entities are introduced. The main unverified inputs are the accuracy of the simulation model and the feasibility of the assumed 0.1 log10 multiplicity uncertainty in the presence of systematic detector effects.

free parameters (3)
  • 12 spline coefficients c_i for the pion-air inelastic cross section = nominal 0, prior width 5
    These coefficients are the tuning parameters in Eq. 5 and the fitted parameters in the chi-square analysis of Eq. 6; the projected cross-section uncertainties are propagated from their covariance in Eq. 7.
  • Per-bin multiplicity uncertainty log10(Delta M/M) = 0.1
    Assumed experimental precision motivated by IceCube multiplicity reconstruction studies [15]; this number directly sets the scale of the projected uncertainties in Eq. 6.
  • Prior width delta c_i = 5
    Chosen by hand to regularize the spline fit and prevent unphysical boundary distortions; it affects the inverse Hessian and therefore the projected uncertainties.
assumptions (5)
  • domain assumption SIBYLL 2.3c provides an adequate nominal description of all hadronic interactions except the modified cross sections.
    The pseudo-data M_nominal and all gradients are computed with SIBYLL 2.3c; if the model misrepresents particle production, the projected sensitivity is biased. Invoked throughout Sec. 2.1 and Sec. 2.2.
  • domain assumption The superposition approximation is valid for cosmic-ray nuclei.
    Stated in Sec. 2.1 to justify simulating only protons and using a nucleon flux; this neglects composition fluctuations and any composition uncertainty that could degenerate with cross-section changes.
  • domain assumption The Global Spline Fit (GSF) cosmic-ray flux and composition model is accurate.
    Used in Eq. 2 and Fig. 2 to weight primary energies; CR flux and composition uncertainties are not propagated in the sensitivity calculation.
  • domain assumption The muon multiplicity is a monotonic function of primary energy.
    Needed for the transformation in Eq. 3, where dE_CR/dN_mu is computed; monotonicity is asserted without proof in Sec. 2.
  • domain assumption Kaon-air cross section variations contribute negligibly.
    Kaon effects are shown to be smaller in Fig. 3 and are excluded from the detailed sensitivity analysis; if kaon contributions grow at other depths or angles, they could bias the pion-air estimate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A New Method for Measuring the Pion-Air Cross Section at Multi-TeV Energies Using Muon Bundle Properties in Deep Underground Detectors." pith.science (2026). https://pith.science/paper/ULTDYH3B

@misc{pith2026250722463,
  author       = {Pith},
  title        = {Pith review of: A New Method for Measuring the Pion-Air Cross Section at Multi-TeV Energies Using Muon Bundle Properties in Deep Underground Detectors},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ULTDYH3B}},
  note         = {Machine review of arXiv:2507.22463}
}
read the original abstract

The interaction cross section of charged pions with air nuclei is a critical parameter for accurately simulating extensive air showers. Improving the modeling of high-energy pion interactions is essential for addressing the muon puzzle-the observed deficit of muons in simulations compared to indirect experimental estimates. As collider experiments cannot directly probe these interactions, we propose a novel measurement approach using muon bundles detected in deep-underground water Cherenkov detectors, such as IceCube and KM3NeT. This method aims to constrain the pion-air inelastic cross section, thereby reducing uncertainties in air shower simulations and advancing our understanding of cosmic ray interactions.

Figures

Figures reproduced from arXiv: 2507.22463 by the authors.

Figure 1
Figure 1. Top panel: Vertical muon flux at the surface for a primary proton with an energy of 1 PeV, computed at the South Pole using the SIBYLL 2.3c model. Relative contributions from different parent particle species are indicated. Bottom panel: Ratio of the modified muon flux when increasing the pion￾air cross section 𝜎𝜋+air by 50% in the ∼1 TeV region relative to the nominal cross section. 0.0 0.5 1.0 1.5 2.0 2.5 3.0 log1… view at source ↗
Figure 3
Figure 3. Illustration of the impact on the muon multiplicity distribution resulting from variations of the inelastic cross sections within specific energy ranges defined by spline coefficients, relative to the nominal SIBYLL 2.3c predictions. Left panels: Effects on the multiplicity distribution at a depth of 𝑑 = 1.5 km w.e. induced by a 20% increase (𝑐𝑖 = +0.2) in the pion-, kaon-, and proton-air cross sections. Right panel… view at source ↗
Figure 4
Figure 4. Projected sensitivities to 𝜎𝜋+air and 𝜎p+air at various underground depths. The shift in the energy range of maximum sensitivity with increasing depth is clearly evident. 6 [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

15 extracted references · 2 canonical work pages

  1. [1]

    Fedynitch, W

    A. Fedynitch, W. Woodley and M.-C. Piro, On the Accuracy of Underground Muon Intensity Calculations , https://doi.org/10.3847/1538-4357/ac5027 Astrophys. J. 928 (2022) 27 [ https://arxiv.org/abs/2109.11559 2109.11559 ]

    arXiv 2022

  2. [2]

    Fedynitch and M

    A. Fedynitch and M. Huber, Data-driven hadronic interaction model for atmospheric lepton flux calculations , https://doi.org/10.1103/PhysRevD.106.083018 Phys. Rev. D 106 (2022) 083018 [ https://arxiv.org/abs/2205.14766 2205.14766 ]

    arXiv 2022

  3. [3]

    Ya \ n ez and A

    J.P. Ya \ n ez and A. Fedynitch, Data-driven muon-calibrated neutrino flux , https://doi.org/10.1103/PhysRevD.107.123037 Phys. Rev. D 107 (2023) 123037 [ https://arxiv.org/abs/2303.00022 2303.00022 ]

    arXiv 2023

  4. [4]

    Woodley, A

    W. Woodley, A. Fedynitch and M.-C. Piro, Cosmic ray muons in laboratories deep underground , https://doi.org/10.1103/PhysRevD.110.063006 Phys. Rev. D 110 (2024) 063006 [ https://arxiv.org/abs/2406.10339 2406.10339 ]

    arXiv 2024

  5. [5]

    Riehn, A

    F. Riehn, A. Fedynitch and R. Engel, Sibyll ^ , https://doi.org/10.1016/j.astropartphys.2024.102964 Astropart. Phys. 160 (2024) 102964 [ https://arxiv.org/abs/2404.02636 2404.02636 ]

    arXiv 2024

  6. [6]

    Albrecht et al., The Muon Puzzle in cosmic-ray induced air showers and its connection to the Large Hadron Collider , https://doi.org/10.1007/s10509-022-04054-5 Astrophys

    J. Albrecht et al., The Muon Puzzle in cosmic-ray induced air showers and its connection to the Large Hadron Collider , https://doi.org/10.1007/s10509-022-04054-5 Astrophys. Space Sci. 367 (2022) 27 [ https://arxiv.org/abs/2105.06148 2105.06148 ]

    arXiv 2022

  7. [7]

    NA61/SHINE collaboration, Measurement of hadron production in --C interactions at 158 and 350 \, \, GeV/c with NA61/SHINE at the CERN SPS , https://doi.org/10.1103/PhysRevD.107.062004 Phys. Rev. D 107 (2023) 062004 [ https://arxiv.org/abs/2209.10561 2209.10561 ]

    work page Pith review arXiv 2023

  8. [8]

    Pierre Auger collaboration, Measurement of the proton-air cross-section at s =57 TeV with the Pierre Auger Observatory , https://doi.org/10.1103/PhysRevLett.109.062002 Phys. Rev. Lett. 109 (2012) 062002 [ https://arxiv.org/abs/1208.1520 1208.1520 ]

    arXiv 2012

Show all 15 references
  1. [9]

    R.U. Abbasi et al., Measurement of the proton-air cross section with Telescope Array s Black Rock Mesa and Long Ridge fluorescence detectors, and surface array in hybrid mode , https://doi.org/10.1103/PhysRevD.102.062004 Phys. Rev. D 102 (2020) 062004 [ https://arxiv.org/abs/2...

    2020

  2. [10]

    Fedynitch, R

    A. Fedynitch, R. Engel, T.K. Gaisser, F. Riehn and T. Stanev, Calculation of conventional and prompt lepton fluxes at very high energy , https://doi.org/10.1051/epjconf/20159908001 EPJ Web Conf. 99 (2015) 08001 [ https://arxiv.org/abs/1503.00544 1503.00544 ]

    2015

  3. [11]

    Fedynitch, F

    A. Fedynitch, F. Riehn, R. Engel, T.K. Gaisser and T. Stanev, Hadronic interaction model sibyll 2.3c and inclusive lepton fluxes , https://doi.org/10.1103/PhysRevD.100.103018 Phys. Rev. D 100 (2019) 103018 [ https://arxiv.org/abs/1806.04140 1806.04140 ]

    2019 arXiv

  4. [12]

    Dembinski, R

    H.P. Dembinski, R. Engel, A. Fedynitch, T. Gaisser, F. Riehn and T. Stanev, Data-driven model of the cosmic-ray flux and mass composition from 10 GeV to 10^ 11 GeV , https://doi.org/10.22323/1.301.0533 PoS ICRC2017 (2018) 533 [ https://arxiv.org/abs/1711.11432 1711.11432 ]

    2018 arXiv

  5. [13]

    IceCube collaboration, The IceCube Neutrino Observatory: Instrumentation and Online Systems , https://doi.org/10.1088/1748-0221/12/03/P03012 JINST 12 (2017) P03012 [ https://arxiv.org/abs/1612.05093 1612.05093 ]

    2017 arXiv

  6. [14]

    Riehn, R

    F. Riehn, R. Engel, A. Fedynitch, T.K. Gaisser and T. Stanev, Hadronic interaction model Sibyll 2.3d and extensive air showers , https://doi.org/10.1103/PhysRevD.102.063002 Phys. Rev. D 102 (2020) 063002 [ https://arxiv.org/abs/1912.03300 1912.03300 ]

    2020 arXiv

  7. [15]

    IceCube collaboration, Measurement of the mean number of muons with energies above 500 GeV in air showers detected with the IceCube Neutrino Observatory , Submitted to Phys. Rev. D (2025) [ https://arxiv.org/abs/2506.19241 2506.19241 ]

    2025

Pith tools

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