Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Signal Model and Energy Reconstruction for the Radio Detection of Inclined Air Showers in the 50-200 MHz Frequency Band

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

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

Pith's one-line read Adapting a radio signal model to the 50–200 MHz band, this paper reconstructs electromagnetic air-shower energies with intrinsic resolution below 5% at both the Auger and GRANDProto300 sites, and below 10% on realistic sparse antenna…

desk verdict Useful frequency-band extension of the Schlüter–Huege radio signal model, but the <5% resolution is an in-sample fit on the tuning libraries; the sparse-array <10% results are the stronger, out-of-sample evidence. read the letter →

arxiv 2507.06698 v1 pith:EYWG5IN6 submitted 2025-07-09 astro-ph.IM astro-ph.HEhep-ex

Lukas Gülzow , Ralph Engel , Tim Huege , Markus Roth , Felix Schlüter This is my paper · ORCID
classification astro-ph.IMastro-ph.HEhep-ex
keywords cosmic-rayairshowersradiodetectionenergyreconstructionlateraldistributionfunctiongeomagneticradiationchargeexcessCoREASGRANDProto300
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

This paper extends an existing radio-based energy reconstruction method for inclined cosmic-ray air showers from the 30–80 MHz band up to the 50–200 MHz band used by the GRAND detector family. The method isolates the geomagnetic component of the radio emission, fits a lateral distribution function to it, and applies geometry and air-density corrections to link the measured radiation energy to the electromagnetic shower energy. On CoREAS simulations, the reconstruction achieves an intrinsic energy resolution below 5% at the Pierre Auger site in Argentina and at the GRANDProto300 site in China, despite the latter having a magnetic field almost three times stronger. When tested on realistic sparse antenna layouts with added noise, the resolution stays below 10% with negligible bias, which is what matters for building practical radio arrays.

What carries the argument

The central object is the lateral distribution function (LDF) of the geomagnetic energy fluence, Eq. (3), a sum of a Gaussian peak and a sigmoid term whose seven shape parameters are parametrised as functions of the distance to shower maximum, $d_{\max}$, so that the fit reduces to four degrees of freedom: the radiation energy $E_{\mathrm{geo}}$ and the two radio-core coordinates. The LDF is fitted after an early-late correction removes geometric asymmetries and a per-site parametrisation of the charge-excess fraction isolates the geomagnetic component. The second load-bearing element is the modified density correction, Eq. (7), which compensates coherence loss before the corrected radiation energy is converted to electromagnetic energy with a power law, Eq. (8).

What would settle it

Generate independent CoREAS libraries with a different hadronic interaction model (e.g., EPOS-LHC vs. Sibyll2.3) or with an independent simulation code such as ZHAireS, re-run the fixed tuning of this paper, and check whether the intrinsic resolution stays below 5% and the sparse-array resolution below 10%. If the reconstruction degrades significantly, the quoted resolutions are tuned to one simulation library rather than to the physics.

Watch

Extended reading notes

Core claim

The paper's central claim is that the geomagnetic radiation energy of an inclined air shower, reconstructed with a parametrised lateral distribution function fitted to the charge-excess-subtracted radio fluence, determines the electromagnetic shower energy to better than 5% intrinsic resolution at two benchmark sites with very different magnetic fields. The same reconstruction applied to simulated sparse arrays, including instrumental noise, smearing and timing jitter, yields a resolution better than 10% with negligible bias. To reach this, the authors re-fit the charge-excess fraction parametrisation for the higher-frequency band, express all lateral-distribution shape parameters as functions of the distance to shower maximum, and introduce a modified density correction that absorbs the coherence loss caused by shorter wavelengths and, especially in China, the strong magnetic field.

Load-bearing premise

The central calibration and the headline resolutions both come from CoREAS simulations, so the whole argument collapses if CoREAS does not correctly predict 50–200 MHz radio emission, particularly the coherence loss in strong magnetic fields.

Editorial extensions

If this is right

  • If the simulation results carry over to data, GRANDProto300 can achieve sub-10% energy resolution for inclined showers in the $10^{17}$–$10^{20}$ eV range on its sparse antenna grid.
  • The same reconstruction works for a very large array of 10,000 km$^2$ with 1 km spacing, again with $<10\%$ resolution, which is relevant for the full GRAND design.
  • The fit requiring only five antennas with signal implies that even partial or prototype arrays can provide per-shower energy reconstruction.
  • At low cosmic-ray energies, where signal approaches the noise level, the resolution degrades, setting the practical energy threshold of the method.

Reading between the lines

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

  • The $<5\%$ figure is computed on the same star-shaped simulation libraries used to tune the LDF parameters; an independent test on simulations generated with a different interaction model or code is needed before treating it as a guaranteed error budget.
  • The realistic benchmark adds Gaussian noise, amplitude smearing and timing jitter, but does not simulate real radio-frequency interference, antenna gain errors or atmospheric uncertainties, so real-data resolution is likely to be somewhat worse.
  • Since the model's $S_{19}$ parameter, the reference radiation energy at 10 EeV, differs between the two sites as expected from the magnetic field scaling, a cross-check against an absolute energy scale from fluorescence or surface detectors could validate the model's normalisation.
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 / 5 minor

Summary. The paper adapts the radio-detection signal model of Schlüter and Huege from the 30-80 MHz band to the 50-200 MHz band, for the Pierre Auger site in Argentina and the GRANDProto300 site in China. It introduces a refitted charge-excess parametrization, a modified lateral distribution function with shape parameters parametrized by the distance to shower maximum, and a density correction to account for coherence loss in the stronger Chinese magnetic field. Using CoREAS simulations, the authors report intrinsic energy resolutions below 5% for ideal star-shaped antenna layouts at both sites, and below 10% for sparse, noisy arrays representing GRANDProto300 and a 10,000 km^2 array.

Significance. If supported, an energy resolution below 5% for inclined air showers in the 50-200 MHz band, with a separate demonstration below 10% on realistic sparse arrays, would be a valuable extension of the established 30-80 MHz method and directly relevant to GRAND and other high-frequency radio detectors. The paper's concrete two-site parameter tables, the explicit LDF functional form, and the independent simulation sets with noise, amplitude smearing, and time jitter are useful contributions; the GP300 sparse-array benchmark is a genuine out-of-sample test of the reconstruction pipeline. However, the headline <5% claim is currently evaluated on the same simulation libraries used to tune the model parameters, so it is best interpreted as an in-sample consistency measure rather than a demonstrated predictive resolution.

major comments (3)
  1. [Section 4, first paragraph and Fig. 4] The <5% resolutions in Fig. 4 are obtained on 'the simulation libraries with star-shaped antenna layout with which we tuned the parameters.' All model ingredients, including the charge-excess coefficients in Eqs. (1)-(2), the LDF shape coefficients in Table 1, and the density and energy calibration parameters in Eqs. (7)-(8), are fitted to these same libraries. With roughly 4,000 events and tens of fitted parameters, the quoted scatter is not an unbiased estimate of predictive performance. Please add an explicit train/test split or K-fold cross-validation, for example by fitting on disjoint subsets in energy and zenith angle and reporting the held-out resolution; alternatively, if the in-sample 'intrinsic' resolution is the intended quantity, the claim should be reworded and supported by a stability analysis such as bootstrap uncertainties on the fitted parameters.
  2. [Section 4, Fig. 5] The sparse-array simulations are independent of the tuning libraries, which is a strength, but they are performed only for the China site and they reuse the star-shaped China calibration values shown in Fig. 4, including S19=14.16 GeV and gamma=1.9970. The sparse-array results therefore validate the LDF fit and the full reconstruction pipeline for one magnetic-field configuration only; they do not by themselves support the conclusion that the method is readily adaptable to any magnetic field configuration at sparse arrays. Please add a comparable out-of-sample sparse-array test for the Argentina site, or clearly qualify the conclusion to state that sparse-array performance has been demonstrated only for China.
  3. [Section 4, Figs. 4-5] No statistical uncertainties or per-bin event counts are reported for the resolution and bias values. The resolution is derived from bin-wise distributions, and the '<5%' and '<10%' statements may depend on the binning and on finite Monte Carlo statistics. Without uncertainties or event counts, it is difficult to judge whether differences between the two sites or between the ideal and sparse configurations are significant. Please add estimates of the statistical uncertainty on each resolution and bias point, or at least provide the number of events per bin.
minor comments (5)
  1. [Section 1] The sentence 'the atmosphere models are provided by the radiotools package [11])' contains an extra closing parenthesis after the citation; please remove it.
  2. [Section 2, Eq. (4)] The expression for p(r) in the r >= r0 branch is ambiguous: it should likely read p(r) = 2*(r0/r)^(b/1000) rather than the printed '2*(r0/r)^b/1000'. Please add explicit parentheses so that the mild decrease of the exponent with distance is unambiguous.
  3. [Section 2, Table 1] The table caption states that the first four parameters use the cpar(dmax) term of Eq. (6), but it would help to state the units of dmax explicitly in the caption or in the text around Eq. (6), since the coefficients carry km and km^2 units.
  4. [Section 3, Eq. (7)] The density-correction formula mixes parameters with different dimensions inside one expression; a short sentence defining the units of p2 and p3, and noting that the denominator is dimensionless, would improve readability.
  5. [General] The paper would benefit from a data-availability or reproducibility statement describing how the CoREAS libraries and analysis scripts can be accessed, since the numerical parameter tables alone do not allow the fits to be reproduced.

Circularity Check

1 steps flagged · score 6.0 of 10

The headline <5% intrinsic resolution is a training-set residual: Fig. 4 is evaluated on the same CoREAS libraries used to fit the model parameters, so it measures self-consistency, not predictive accuracy.

  1. fitted input called prediction [Section 4, first paragraph and Fig. 4; parameters fitted in Section 3 via Eqs. (7)-(8)]
    "The first benchmark for our method is the intrinsic reconstruction performance on the simulation libraries with star-shaped antenna layout with which we tuned the parameters, with no added noise. ... We determine the parameter values of p0, p1, p2, p3, c_alpha, S19 and gamma in a joint fit of Eqs. (7) and (8)."

    The <5% resolution in Fig. 4 is computed on the same ~4,000-event CoREAS libraries on which the charge-excess parametrisations (Eqs. 1-2), the LDF shape parameters (Table 1), the density correction, and the energy-scale parameters S19 and gamma in Eq. (8) were all fitted. Because S19 and gamma are determined by a joint fit of Eqs. (7)-(8) to exactly these events, the reconstructed-to-true energy ratio shown as 'resolution' is the training residual of that fit, not an out-of-sample prediction. The paper reports no train/test split and no parameter uncertainties for the star-shaped libraries, so the <5% number is a self-consistency statement.

full rationale

The paper is transparent that the first benchmark uses the tuning libraries, which is precisely the problem: every load-bearing parameter of the reconstruction, including the power-law energy scale S19 and exponent gamma of Eq. (8), is fitted to those libraries, and the quoted <5% resolution is the scatter of those same events around the fitted relation. That makes the headline intrinsic-resolution claim an in-sample goodness-of-fit rather than a predictive result. The sparse-array benchmarks with added noise are genuine out-of-sample simulations for the site in China, which is why the circularity is partial (score 6) rather than complete; however they do not isolate the <5% star-shaped claim and use parameters fixed on the star-shaped libraries. Concerns that CoREAS may not describe real 50-200 MHz emission, or that parameter uncertainties are absent, are correctness/robustness issues rather than circularity and are not counted toward the score beyond the in-sample benchmark.

Assumptions & free parameters 13 free parameters · 7 assumptions · 0 invented entities

The central claim rests on a large number of fitted constants: charge-excess coefficients, six LDF shape parameterizations per site, a fixed sigmoid slope, density-correction parameters, and energy-scale parameters. These are not derived from first principles; they are tuned to CoREAS output. The modeling assumptions (LDF form, density-correction form, power-law Eem-Sgeo relation, early-late correction) come from prior work or are introduced phenomenologically. No new physical entities are invented.

free parameters (13)
  • Charge-excess parametrization coefficients, Argentina (Eq. 1) = 0.302, 729 km, 682 m, 2.98, 0.178
    Fitted to CoREAS star-shaped libraries to describe the charge-excess fraction in the 50-200 MHz band; used to isolate the geomagnetic fluence.
  • Charge-excess parametrization coefficients, China (Eq. 2) = 0.229, 1106 km, 614 m, 1.43, 0.166
    Fitted to CoREAS star-shaped libraries for the Dunhuang site; used to isolate the geomagnetic fluence.
  • LDF shape-parameter coefficients for r02 (Table 1) = ARG: 0.666, 771.9 km, 98.65 km^2; CHN: 0.586, 1176 km, 166.3 km^2
    Fitted iteratively to LDF shape behavior versus d_max; the sigmoid length scale in Eq. (3).
  • LDF shape-parameter coefficients for p_inner (Table 1) = ARG: 1.464, 18982 km, 32.94 km^2; CHN: 1.541, 3810 km, 64.99 km^2
    Fitted Gaussian exponent inside the Cherenkov radius, Eq. (4).
  • LDF shape-parameter coefficients for a_rel (Table 1) = ARG: 0.233, 4848 km, 3.79 km^2; CHN: 0.281, 2173 km, 34.78 km^2
    Relative sigmoid amplitude in Eq. (3).
  • LDF shape-parameter coefficients for b (Table 1) = ARG: 282.2, 2.73 km, 6457 km^2; CHN: 249.7, 4.21 km, 14561 km^2
    Fitted falloff exponent for r >= r0 in Eq. (4).
  • LDF shape-parameter coefficients for sigma (Table 1) = ARG: 0.027, 0.805, 61.97; CHN: 0.035, 0.770, 62.94
    Gaussian width of the LDF, Eq. (5).
  • LDF shape-parameter coefficients for r0 (Table 1) = ARG: 0.941, 4536 km, 15.96 km^2; CHN: 0.818, 1776 km, 36.82 km^2
    Cherenkov-ring radius correction, Eq. (6).
  • Sigmoid slope s in Eq. (3) = Fixed to a constant, value not stated in the text
    Chosen to make the sigmoid significant only inside the Cherenkov radius; the exact value is not reported, hampering exact reproduction.
  • Density correction parameters p0, p1, p2, p3, c_alpha, Argentina (Eq. 7, Fig. 3) = p0=0.96, p1=-1.30/(kg m^-3), p2=0.00 kg m^-3, p3=0.00, c_alpha=2.00
    Determined in a joint fit of Eqs. (7) and (8) to CoREAS libraries.
  • Density correction parameters p0, p1, p2, p3, c_alpha, China (Eq. 7, Fig. 3) = p0=47.90, p1=-0.01/(kg m^-3), p2=0.03 kg m^-3, p3=0.19, c_alpha=1.06
    Determined in a joint fit of Eqs. (7) and (8) to CoREAS libraries; captures coherence loss in the strong Chinese magnetic field.
  • Energy calibration S19 and gamma, Argentina (Eq. 8, Fig. 4) = S19=5.82 GeV, gamma=1.9933
    Fitted jointly with density-correction parameters to the Sgeo-Eem correlation.
  • Energy calibration S19 and gamma, China (Eq. 8, Fig. 4) = S19=14.16 GeV, gamma=1.9970
    Fitted jointly with density-correction parameters to the Sgeo-Eem correlation.
assumptions (7)
  • domain assumption CoREAS accurately simulates radio emission of inclined air showers at 50-200 MHz, including geomagnetic and charge-excess components and coherence loss.
    All parameterizations and performance numbers are derived from CoREAS libraries; no real-data validation is reported (Sections 1 and 4).
  • domain assumption The early-late correction of Ref. [6] removes geometry-dependent asymmetries in the 50-200 MHz band.
    Applied as the first step of the signal model (Section 2, Fig. 1).
  • ad hoc to paper The functional form of the charge-excess fraction from Ref. [1] (Eqs. 4.10 and 4.12) remains valid; only its coefficients need refitting.
    The model reuses the analytic form and refits coefficients (Eqs. 1-2); no derivation of the form is given.
  • ad hoc to paper The lateral distribution function of Eq. (3), a Gaussian plus sigmoid with power-law exponent, is an adequate model of geomagnetic fluence versus axis distance.
    The LDF is introduced phenomenologically and its shape parameters are fitted (Section 2).
  • ad hoc to paper The density correction of Eq. (7) with the chosen p0-p3 and c_alpha functional form captures the coherence-loss dependence on density and geomagnetic angle.
    Explicitly modified from Eq. (5.2) of Ref. [1]; parameters determined in a joint fit (Section 3).
  • domain assumption The relation between corrected radiation energy and electromagnetic energy is a single power law (Eq. 8).
    Standard assumption from previous work; slope and normalization are fitted.
  • domain assumption The simulated noise model (64 uV/m Gaussian noise, 7.5% amplitude smearing, 5 ns time jitter) represents realistic detector conditions for the sparse arrays.
    Used in the sparse-array benchmarks (Section 4); no comparison with measured noise is shown.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Signal Model and Energy Reconstruction for the Radio Detection of Inclined Air Showers in the 50-200 MHz Frequency Band." pith.science (2026). https://pith.science/paper/EYWG5IN6

@misc{pith2026250706698,
  author       = {Pith},
  title        = {Pith review of: Signal Model and Energy Reconstruction for the Radio Detection of Inclined Air Showers in the 50-200 MHz Frequency Band},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EYWG5IN6}},
  note         = {Machine review of arXiv:2507.06698}
}
read the original abstract

The radio emission of cosmic-ray air-showers changes significantly depending on parameters like signal frequency, magnetic field configuration and observing altitude. We use CoREAS simulations to adapt an existing signal model for the radio emission of inclined showers in the 30-80 MHz frequency band to the wide 50-200 MHz band. Our model uses a parametrisation of the charge excess fraction to isolate the geomagnetic emission component. We reconstruct the geomagnetic radiation energy by fitting a lateral distribution function, provided by the model, to the geomagnetic energy fluence distribution of the shower. After we correct for the shower geometry and air density, we correlate the radiation energy with the electromagnetic energy of the shower. We show that the method intrinsic energy resolutions < 5% for the sites of the Pierre Auger Observatory and GRANDProto300. For GRANDProto300, we test the reconstruction with simulations of a realistic, sparse antenna grid and with added noise, and find an energy resolution of < 10% with negligible bias. We do a similar study for a much larger array of 10, 000 km2 with 1 km antenna spacing. We find an intrinsic energy resolution of < 10%.

Figures

Figures reproduced from arXiv: 2507.06698 by the authors.

Figure 1
Figure 1. The first two steps of our signal model which eliminate geometric asymmetries and to isolate the geomagnetic energy fluence, respectively. Shown for an example simulation of the GRAND@Auger site with a star-shaped antenna pattern (white dots). The regions with 90% or more of the maximum energy fluence are outlined in grey. Axis distances 𝑟 are normalised with respect to the Cherenkov radius 𝑟0. Left: Early-late corr… view at source ↗
Figure 2
Figure 2. Fit of 𝑓LDF(𝑟) (solid black line) to the geomagnetic fluence 𝑓 par geo (blue points) of an example air-shower simulation with a star-shaped antenna pattern using the properties of the GP300 site in China. We display the Gaussian and sigmoid components of the LDF as dot-dashed and dashed lines, respectively. The coloured area under the curve represents the geomagnetic radiation energy 𝐸geo of the air-shower. and the … view at source ↗
Figure 3
Figure 3. Modelled density correction term of Eq. (7) for the calculation of 𝑆geo plotted against 𝜌max. The colour map shows the dependency on the geomagnetic angle 𝛼. The black line is determined in a joint fit of Eqs. (7) and (8), which also provides values for 𝑆19 and 𝛾 (see [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Reconstruction performances for the simulation sets without noise and star-shaped antenna layouts. The left column of each panel shows the correlation of the reconstructed 𝑆geo against the true 𝐸 MC em . The dashed line represents the power law from Eq. (8). The right …
Figure 5
Figure 5. Figure 5: Reconstruction performances for the simulation sets with artificial noise and sparse antenna layouts. Both panels have the same properties as the panels in [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Detector-level simulation closure of radio air-shower reconstruction with native RF-chain inversion and out-of-fold endpoint interpolation

    astro-ph.IM 2026-07 conditional novelty 6.0 of 10

    Detector-level RF-chain-triggered ZHAireS packages reconstruct energy to ~14.3–14.5% scatter and direction to 0.05° median residual using fixed equal-prior iron–proton geomagnetic templates, with strong residual zenit...

Reference graph

Works this paper leans on

13 extracted references · 7 canonical work pages · cited by 1 Pith paper

  1. [1]

    Schlüter and T

    F. Schlüter and T. Huege, JCAP 01(2023) 008 [2203.04364]

    arXiv 2023

  2. [2]

    Hörandel et al., PoS(ICRC2025)294,(2025)

    J. Hörandel et al., PoS(ICRC2025)294,(2025)

    work page 2025

  3. [3]

    The Radio Detector of the Pierre Auger Observatory -- status and expected performance

    Pierre Augercollaboration, EPJ Web Conf. 283(2023) 06002 [2305.10104]

    work page Pith review arXiv 2023

  4. [4]

    China Phys

    GRANDcollaboration, Sci. China Phys. Mech. Astron. 63 (2020) 219501 [1810.09994]

    arXiv 2020

  5. [5]

    Huege, Phys

    T. Huege, Phys. Rept. 620(2016) 1 [1601.07426]

    arXiv 2016

  6. [6]

    A Rotationally Symmetric Lateral Distribution Function for Radio Emission from Inclined Air Showers

    T. Huege, L. Brenk and F. Schlüter, EPJ Web Conf. 216(2019) 03009 [1808.00729]

    work page Pith review arXiv 2019

  7. [7]

    Huege, M

    T. Huege, M. Ludwig and C.W. James, AIP Conf. Proc. 1535(2013) 128 [1301.2132]

    arXiv 2013

  8. [8]

    de Mello Neto et al., PoS(ICRC2025)1024,(2025)

    J.T. de Mello Neto et al., PoS(ICRC2025)1024,(2025)

    work page 2025

Show all 13 references
  1. [9]

    Verpoest et al., PoS(ICRC2025)428,(2025)

    S. Verpoest et al., PoS(ICRC2025)428,(2025)

    2025

  2. [10]

    Ma et al., PoS(ICRC2025)453,(2025)

    P. Ma et al., PoS(ICRC2025)453,(2025)

    2025

  3. [11]

    Glaser, A

    C. Glaser, A. Nelles, I. Plaisier et al.The European Physical Journal C 79 (2019)

    2019

  4. [12]

    Chiche, C

    S. Chiche, C. Zhang, F. Schlüter et al., Phys. Rev. Lett. 132(2024) 231001 [2404.14541]

    2024 arXiv

  5. [13]

    Schlüter, M

    F. Schlüter, M. Gottowik, T. Huege et al., Eur. Phys. J. C 80 (2020) 643 [2005.06775]. 8 Radio Signal Model and Energy Reconstruction for Inclined Air Showers in 50-200 MHz Tim Huege Acknowledgments This work is part of the NUTRIG project, supported by the Agence Nationale de ...

    2020 arXiv

Pith tools

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