pith. sign in

arxiv: 2605.15762 · v1 · pith:YO3HJ5ZZnew · submitted 2026-05-15 · 🌌 astro-ph.HE · hep-ex· hep-ph

Average universal shower profile reconstruction using radio interferometry

J. Alvarez-Mu\~niz , W. R. Carvalho Jr. , R. Concei\c{c}\~ao , D. Dias This is my paper

Pith reviewed 2026-05-20 17:01 UTC · model grok-4.3

classification 🌌 astro-ph.HE hep-exhep-ph
keywords radio detectionextensive air showersuniversal shower profileGaisser-Hillas parameterscosmic ray compositionhadronic interaction modelsinterferometryultra-high energy
0
0 comments X

The pith

Radio interferometry reconstructs the average universal shower profile of air showers, yielding Gaisser-Hillas parameters that separate primary masses and hadronic models more effectively than fluorescence profiles.

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

The paper shows how radio data collected through interferometry can be used to build an average longitudinal profile for extensive air showers and then fit the Gaisser-Hillas shape parameters R and L to that profile. A reader would care because these parameters carry information on whether the primary particle was a proton or a heavier nucleus and on the details of high-energy hadronic collisions. The work demonstrates that the radio-derived average profiles give clearer separation between different primaries and different interaction models than the equivalent profiles obtained from fluorescence measurements.

Core claim

Using Monte Carlo simulations, the authors reconstruct the average longitudinal profile directly from radio data and extract the Gaisser-Hillas shape parameters (R, L). They find that the radio-derived average profiles provide enhanced separation between primary masses and hadronic interaction models compared to that obtained from fluorescence-equivalent longitudinal profiles.

What carries the argument

Radio interferometric reconstruction of the average Universal Shower Profile, from which the Gaisser-Hillas parameters R and L are extracted to quantify longitudinal development.

If this is right

  • Radio interferometry accesses higher-order information on longitudinal shower development beyond the depth of shower maximum.
  • Average USP parameters improve sensitivity to primary mass composition at ultra-high energies.
  • The same parameters improve sensitivity to the choice of hadronic interaction model.
  • Radio methods can supply continuous, weather-independent measurements that complement optical techniques.

Where Pith is reading between the lines

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

  • Future radio arrays could apply this averaging technique to increase statistics for rare ultra-high-energy events without waiting for clear nights.
  • The approach might allow tighter constraints on hadronic physics once real data from existing radio arrays are analyzed with the same method.
  • If the separation gain holds, radio-derived profiles could reduce systematic uncertainties in composition analyses that currently rely on fluorescence alone.

Load-bearing premise

The Monte Carlo simulations accurately model both the radio emission from real air showers and the longitudinal development that the Gaisser-Hillas parameters describe.

What would settle it

A side-by-side comparison, on the same observed showers, of radio-reconstructed average (R, L) values against those measured independently by fluorescence detectors would show whether the reported gain in separation persists outside simulations.

Figures

Figures reproduced from arXiv: 2605.15762 by D. Dias, J. Alvarez-Mu\~niz, R. Concei\c{c}\~ao, W. R. Carvalho Jr..

Figure 1
Figure 1. Figure 1: FIG. 1. Coherence fluence map in the 2-dimensional plane [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Average USP for [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Normalized average profiles for all datasets, extracted [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Maximum difference between the RIT-derived aver [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: shows that the RIT-derived average USP pro￾vides a clear separation between proton and iron pri￾maries, with a discrimination power comparable to that [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. L (top) and R (bottom) dependence on [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. RIT-derived average profile, applied fit to Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. RIT-derived (blue dots) and [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. RIT-derived average profile, applied fit to Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Correlation between [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. All main datasets averaged map of the maximum [PITH_FULL_IMAGE:figures/full_fig_p010_13.png] view at source ↗
read the original abstract

Radio detection of extensive air showers enables near-continuous observation and precise measurements of the shower geometry and the depth of shower maximum, $X_{\rm max}$. Beyond $X_{\rm max}$, the longitudinal shower development follows a Universal Shower Profile (USP), whose shape parameters contain information on primary mass composition and hadronic interaction models. While previous studies have focused on event-by-event reconstruction of profile parameters, in this work we investigate the reconstruction of the average USP using radio interferometric techniques. Using Monte Carlo simulations, we reconstruct the average longitudinal profile directly from radio data and extract the Gaisser-Hillas shape parameters $(R,L)$. We find that the radio-derived average profiles provide enhanced separation between primary masses and hadronic interaction models compared to that obtained from fluorescence-equivalent longitudinal profiles. These results demonstrate that radio interferometry can access higher-order information on the longitudinal shower development and that the use of the average USP significantly improves the sensitivity to composition and hadronic interaction studies at ultra-high energies.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper claims that radio interferometric reconstruction of the average Universal Shower Profile (USP) from Monte Carlo simulations of extensive air showers yields Gaisser-Hillas parameters (R, L) that provide enhanced separation between primary masses and hadronic interaction models, outperforming fluorescence-equivalent longitudinal profiles. The work focuses on accessing higher-order longitudinal development information beyond X_max using radio data.

Significance. If the central result holds after addressing the simulation details, the approach could meaningfully advance ultra-high-energy cosmic-ray composition and hadronic-interaction studies by leveraging radio arrays' near-continuous coverage and the statistical power of average profiles. The paper correctly identifies that averaging reduces event-by-event fluctuations and that radio interferometry can in principle probe profile shape parameters.

major comments (3)
  1. [Methods and Results sections] The reconstruction pipeline, error propagation, and quantitative separation metrics (e.g., Mahalanobis distance or overlap integrals in (R, L) space) are not described with sufficient technical detail to allow independent verification of the claimed improvement. This is load-bearing for the strongest claim.
  2. [Simulation setup and comparison to fluorescence profiles] All comparisons are performed inside the same hadronic-interaction models that simultaneously generate the input longitudinal profiles and the simulated radio signals. Any model-specific bias in radio yield, coherence, or the assumed universality of the USP therefore propagates directly into the reported separation gain, without an external anchor such as real data or an independent radio-emission calculation.
  3. [Monte Carlo description] The manuscript does not report the number of simulated showers, the energy range, zenith-angle distribution, or the precise definition of the fluorescence-equivalent profiles used as baseline, all of which are required to assess whether the separation improvement is robust or an artifact of the chosen Monte Carlo ensemble.
minor comments (2)
  1. [Abstract] The abstract states that radio-derived profiles 'provide enhanced separation' but supplies no numerical measure of the improvement; a single sentence quantifying the gain would improve clarity.
  2. [Introduction and profile reconstruction] The exact functional form and fitting procedure for the Gaisser-Hillas parameters R and L are introduced without an equation or reference, leaving the reader to infer the parametrization.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. The comments highlight important areas for improving clarity, reproducibility, and discussion of limitations. We address each major comment point by point below and indicate the planned revisions.

read point-by-point responses

  1. Referee: [Methods and Results sections] The reconstruction pipeline, error propagation, and quantitative separation metrics (e.g., Mahalanobis distance or overlap integrals in (R, L) space) are not described with sufficient technical detail to allow independent verification of the claimed improvement. This is load-bearing for the strongest claim.

    Authors: We agree that the manuscript requires additional technical detail for independent verification. In the revised version, we will expand the Methods section with a dedicated subsection describing the reconstruction pipeline in full, including the interferometric visibility simulation, the profile reconstruction algorithm from radio data, the fitting procedure to extract Gaisser-Hillas parameters (R, L), the error propagation approach (via ensemble resampling of simulated events), and the quantitative separation metrics (Mahalanobis distance between distributions in (R, L) space together with overlap integrals of the corresponding probability densities). These additions will directly support the claimed improvement. revision: yes

  2. Referee: [Simulation setup and comparison to fluorescence profiles] All comparisons are performed inside the same hadronic-interaction models that simultaneously generate the input longitudinal profiles and the simulated radio signals. Any model-specific bias in radio yield, coherence, or the assumed universality of the USP therefore propagates directly into the reported separation gain, without an external anchor such as real data or an independent radio-emission calculation.

    Authors: We acknowledge that the study is conducted entirely within a consistent Monte Carlo framework using the same hadronic models for both fluorescence profiles and radio signals. This design enables a controlled, apples-to-apples comparison of the information content accessible by each technique under identical assumptions about shower development and radio emission. We will add an explicit paragraph in the Discussion section noting this as an inherent limitation of simulation-based work and discussing how model-specific biases could influence the quantitative separation gains, while emphasizing that the relative improvement remains informative within current theoretical models. revision: partial

  3. Referee: [Monte Carlo description] The manuscript does not report the number of simulated showers, the energy range, zenith-angle distribution, or the precise definition of the fluorescence-equivalent profiles used as baseline, all of which are required to assess whether the separation improvement is robust or an artifact of the chosen Monte Carlo ensemble.

    Authors: We thank the referee for identifying these omissions. In the revised manuscript, the Simulation Setup section will be expanded to report the total number of simulated showers (broken down by primary type and hadronic model), the primary energy range, the zenith-angle distribution, and a precise definition of the fluorescence-equivalent profiles (constructed by integrating the simulated energy deposit along the shower axis with standard fluorescence yield models). These details will allow readers to evaluate the robustness of the reported separation improvement. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results are simulation-based comparisons

full rationale

The paper reports a Monte Carlo study in which average universal shower profiles are reconstructed from simulated radio interferometry data, Gaisser-Hillas parameters (R, L) are extracted, and separation power between primaries and hadronic models is compared against fluorescence-equivalent profiles. This outcome is a numerical result of applying the reconstruction pipeline to an ensemble of simulated showers; no equation or parameter is shown to be defined in terms of the claimed separation gain, no fitted input is relabeled as a prediction, and no load-bearing step reduces to a self-citation or ansatz that imports the target result. The derivation chain therefore remains self-contained against the external benchmark of the chosen Monte Carlo generators.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The reconstruction depends on standard domain assumptions about air-shower radio emission and Monte Carlo fidelity; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Monte Carlo simulations accurately represent real extensive air showers and their radio emission.
    The paper uses these simulations to demonstrate the reconstruction and the claimed separation improvement.

pith-pipeline@v0.9.0 · 5718 in / 1128 out tokens · 69603 ms · 2026-05-20T17:01:24.081011+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · 7 internal anchors

  1. [1]

    Engel, D

    R. Engel, D. Heck, and T. Pierog, Extensive air show- ers and hadronic interactions at high energy, Ann. Rev. Nucl. Part. Sci.61, 467 (2011)

    work page 2011

  2. [2]

    Aabet al.(Pierre Auger), Depth of maximum of air- shower profiles at the Pierre Auger Observatory

    A. Aabet al.(Pierre Auger), Depth of maximum of air- shower profiles at the Pierre Auger Observatory. II. Com- position implications, Phys. Rev. D90, 122006 (2014), arXiv:1409.5083 [astro-ph.HE]

    work page arXiv 2014

  3. [3]

    Measurements of the Cosmic Ray Composition with Air Shower Experiments

    K.-H. Kampert and M. Unger, Measurements of the Cos- mic Ray Composition with Air Shower Experiments, As- tropart. Phys.35, 660 (2012), arXiv:1201.0018 [astro- ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  4. [4]

    Aabet al.(Pierre Auger), Measurement of the Fluctu- ations in the Number of Muons in Extensive Air Showers with the Pierre Auger Observatory, Phys

    A. Aabet al.(Pierre Auger), Measurement of the Fluctu- ations in the Number of Muons in Extensive Air Showers with the Pierre Auger Observatory, Phys. Rev. Lett.126, 152002 (2021), arXiv:2102.07797 [hep-ex]

    work page arXiv 2021

  5. [5]

    A. Abdul Halimet al.(Pierre Auger), 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, 102001 (2024), arXiv:2401.10740 [astro-ph.HE]

    work page arXiv 2024

  6. [6]

    Albrechtet al., Global tuning of hadronic interaction models with accelerator-based and astroparticle data, Nature Rev

    J. Albrechtet al., Global tuning of hadronic interaction models with accelerator-based and astroparticle data, Nature Rev. Phys.8, 98 (2026), arXiv:2508.21796 [astro- ph.HE]

    work page arXiv 2026

  7. [7]

    Andringa, R

    S. Andringa, R. Concei¸ c˜ ao, and M. Pimenta, Mass com- position and cross-section from the shape of cosmic ray shower longitudinal profiles, Astropart. Phys.34, 360 (2011)

    work page 2011

  8. [8]

    The muonic longitudinal shower profiles at production

    S. Andringa, L. Cazon, R. Concei¸ c˜ ao, and M. Pimenta, The Muonic longitudinal shower profiles at production, Astropart. Phys.35, 821 (2012), arXiv:1111.1424 [hep- ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  9. [9]

    Concei¸ c˜ ao, S

    R. Concei¸ c˜ ao, S. Andringa, F. Diogo, and M. Pimenta, The average longitudinal air shower profile: exploring the shape information, J. Phys. Conf. Ser.632, 012087 (2015)

    work page 2015

  10. [10]

    A. Aabet al.(Pierre Auger), Measurement of the aver- age shape of longitudinal profiles of cosmic-ray air show- ers at the Pierre Auger Observatory, JCAP03, 018, arXiv:1811.04660 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv

  11. [11]

    Huege, Radio detection of cosmic ray air showers in the digital era, Physics Reports620, 1 (2016)

    T. Huege, Radio detection of cosmic ray air showers in the digital era, Physics Reports620, 1 (2016)

    work page 2016

  12. [12]

    F. D. Kahn and I. Lerche, Radiation from cosmic ray air showers, Proc. Roy. Soc. A289, 206 (1966)

    work page 1966

  13. [13]

    G. A. Askaryan, Excess negative charge of an electron- photon shower and its coherent radio emission, Sov. Phys. JETP14, 441 (1962)

    work page 1962

  14. [14]

    A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

    S. Buitinket al., A large light-mass component of cos- mic rays at 10ˆ17 - 10ˆ17.5 eV from radio observations, Nature531, 70 (2016), arXiv:1603.01594 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  15. [15]

    P. A. Bezyazeekovet al., Reconstruction of cosmic ray air showers with Tunka-Rex data using template fit- ting of radio pulses, Phys. Rev. D97, 122004 (2018), arXiv:1803.06862 [astro-ph.IM]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  16. [16]

    A. Abdul Halimet al.(Pierre Auger), Demonstrating Agreement between Radio and Fluorescence Measure- ments of the Depth of Maximum of Extensive Air Show- ers at the Pierre Auger Observatory, Phys. Rev. Lett. 132, 021001 (2024), arXiv:2310.19963 [astro-ph.HE]

    work page arXiv 2024

  17. [17]

    De Henau, S

    V. De Henau, S. Bouma, J. Bray, S. Buitink, A. Corstanje, M. Desmet, E. Dickinson, L. van Don- gen, B. Hare, and T. Huege, Investigating double bump air showers with the SKA, PoSICRC2025, 236 (2025), arXiv:2510.13788 [astro-ph.HE]

    work page arXiv 2025

  18. [18]

    de Errico and C

    B. de Errico and C. Timmermans, Air shower develop- ment through the time dependence of its induced electric field, (2026), arXiv:2603.05424 [hep-ex]

    work page arXiv 2026

  19. [19]

    Schoorlemmer and W

    H. Schoorlemmer and W. R. Carvalho, Radio interfer- ometry applied to the observation of cosmic-ray induced extensive air showers, The European Physical Journal C 81, 10.1140/epjc/s10052-021-09925-9 (2021)

    work page doi:10.1140/epjc/s10052-021-09925-9 2021

  20. [20]

    Schl¨ uter and T

    F. Schl¨ uter and T. Huege, Expected performance of air-shower measurements with the radio-interferometric technique, Journal of Instrumentation16(07), P07048

    work page

  21. [21]

    Corstanjeet al., LOFAR-style reconstruction of cosmic-ray air showers with SKA-Low, Phys

    A. Corstanjeet al., LOFAR-style reconstruction of cosmic-ray air showers with SKA-Low, Phys. Rev. D112, 023017 (2025), arXiv:2504.16873 [astro-ph.HE]

    work page arXiv 2025

  22. [22]

    S. J. Sciutto, Aires: A system for air shower simulations, arXiv preprint (1999), astro-ph/9911331

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  23. [23]

    E. Zas, F. Halzen, and T. Stanev, Electromagnetic pulses from high-energy showers: Implications for neutrino de- tection, Phys. Rev. D45, 362 (1992)

    work page 1992

  24. [24]

    Alvarez-Mu˜ niz, A

    J. Alvarez-Mu˜ niz, A. Romero-Wolf, and E. Zas, Practical and accurate calculations of askaryan radiation, Phys. Rev. D81, 123009 (2010)

    work page 2010

  25. [25]

    Alvarez-Mu˜ niz, W

    J. Alvarez-Mu˜ niz, W. R. Carvalho, and E. Zas, Monte carlo simulations of radio pulses in atmospheric showers using zhaires, Astroparticle Physics35, 325 (2012)

    work page 2012

  26. [26]

    Riehn, R

    F. Riehn, R. Engel, A. Fedynitch, T. K. Gaisser, and T. Stanev, Hadronic interaction model SIBYLL 2.3d and extensive air showers, Phys. Rev. D102, 063002 (2020)

    work page 2020

  27. [27]

    Monte Carlo treatment of hadronic interactions in enhanced Pomeron scheme: I. QGSJET-II model

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

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  28. [28]

    Abdul Halimet al.(Pierre Auger), Measuring the muon content of inclined air showers using AERA and the water-Cherenkov detectors of the Pierre Auger Observa- tory, Phys

    A. Abdul Halimet al.(Pierre Auger), Measuring the muon content of inclined air showers using AERA and the water-Cherenkov detectors of the Pierre Auger Observa- tory, Phys. Rev. D112, 123042 (2025), arXiv:2507.02558 [astro-ph.HE]

    work page arXiv 2025

  29. [29]

    An additional dataset of 1000 proton showers, generated with SIBYLL 2.3d (ten times larger thanSIB Proton), was used to assess statistical stability. Consistent with thejackkniferesults, the deviation in the (R, L) space be- tween datasets with different statistics is 0.60%, well be- low the 2.18% separation betweenSIB ProtonandQGS Proton. 9 APPENDIX In F...

    work page