REVIEW 4 major objections 6 minor 26 references
Towards a Cosmic-Ray Energy Scale with the Auger Engineering Radio Array
T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The radio energy scale can cross-check the fluorescence-based cosmic-ray energy scale event by event.
desk verdict AERA's energy-scale cross-check plan is coherent, but the 'independent' radio scale rests entirely on uncalibrated CoREAS absolute fluence; this is a useful status report, not yet a result. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing device is the matched-simulation comparison: each measured shower is reconstructed twice, once from the surface detector (yielding $E_{\rm SD}$ on the fluorescence scale) and once from radio signals (yielding $E^{\rm data}_{\rm rad}$ on the radio scale). The surface-detector energy fixes the input energy of a CoREAS simulation, after subtracting the invisible energy, and the simulated event is passed through the same Offline reconstruction with measured noise added per station, giving $E^{\rm sim}_{\rm rad}$. Two calibration chains carry the absolute meaning of the comparison: on the measurement side, lab-measured electronics responses, drone-calibrated antenna patterns, and a per-station Galactic calibration with about 6% systematic uncertainty; on the simulation side, the ~3% agreement between CoREAS and ZHAireS on the energy scale. The radiation energy itself is obtained by fitting the GeoCeLDF lateral distribution to the per-station energy fluences.
What would settle it
Run the planned comparison on the 912 hybrid events and bin the ratio $E^{\rm data}_{\rm rad}/E^{\rm sim}_{\rm rad}$ by zenith angle, energy, and season; if the ratio leaves unity by more than the combined 6% calibration and 3% simulation uncertainties, the two scales disagree, and if the residual tracks a reconstruction ingredient such as the temperature correction or the lateral-distribution fit, the absolute radio-emission prediction of CoREAS would be the part called into question.
Extended reading notes
Core claim
The central claim is that the radio energy scale can provide an independent cross-check of the fluorescence-based cosmic-ray energy scale. For each measured hybrid event, the surface detector gives an energy $E_{\rm SD}$ on the fluorescence scale; this is converted to the calorimetric energy, fed into a CORSIKA/CoREAS simulation of the same shower, and the simulation is reconstructed with the same pipeline and the same measured noise as the real event. The result is a predicted radiation energy $E^{\rm sim}_{\rm rad}$ to compare with the measured $E^{\rm data}_{\rm rad}$. Because CoREAS derives the radio emission from first-principles classical electrodynamics without free parameters, the ratio of these two quantities directly probes whether the fluorescence and radio energy scales agree, without needing simultaneous fluorescence and radio observation of the same shower.
Load-bearing premise
The cross-check assumes that the radio-simulation code CoREAS predicts the absolute amount of radio emission from a real air shower to within a few percent, yet that prediction is backed only by agreement with another simulation code, not by an experimental calibration of absolute radio-emission efficiency that is independent of the fluorescence energy scale.
Editorial extensions
If this is right
- If $E^{\rm data}_{\rm rad}/E^{\rm sim}_{\rm rad}$ is consistent with unity, the radio and fluorescence cosmic-ray energy scales agree to within the combined systematic uncertainties of about 6% (radio calibration) and 3% (simulation agreement).
- A significant offset would localize the discrepancy between the calorimetric fluorescence scale and the first-principles radio scale, motivating a re-examination of either the fluorescence photometric calibration or the absolute radio-emission prediction.
- The 912-event hybrid dataset is large enough to make the cross-check competitive with the 14% systematic uncertainty currently assigned to the fluorescence energy scale.
- Using the surface-detector energy instead of the fluorescence energy as simulation input is what makes the comparison feasible, since simultaneous fluorescence and radio detection is too rare to provide sufficient statistics.
- The same matched-simulation methodology can be applied with the AugerPrime radio detector, extending the cross-check to higher energies and more inclined showers.
Reading between the lines
- A practical use of the data not spelled out in the paper: binning the ratio in zenith angle, energy, and season can separate a true energy-scale offset from a flawed reconstruction ingredient; a flat ratio strengthens the fluorescence scale, while a trend points to a specific step.
- Because the simulation input energy comes from the surface detector on the fluorescence scale, the analysis measures consistency between two scales but cannot, by itself, say which scale is wrong when they disagree.
- The same first-principles simulation chain could in principle be used to compare a radio-based scale with the fluorescence scale of the Telescope Array, which could clarify whether the known spectrum differences between the two observatories originate in their energy scales.
- If the data-simulation ratio differs between proton and iron simulations, the method would reveal sensitivity to the assumed primary mass, an effect not addressed in the comparison as presented.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings contribution from the Pierre Auger Collaboration describes the methodology and current status of an envisaged cross-check of the cosmic-ray energy scale using the Auger Engineering Radio Array (AERA). The plan is to compare, event by event, the radio radiation energy reconstructed from hybrid SD-AERA measurements with the radiation energy predicted by CoREAS simulations whose input energy is set by the SD energy estimator, which is calibrated to the fluorescence-detector (FD) energy scale. The paper reports a hybrid dataset of 912 events, summarizes the radio reconstruction chain (electronics calibration, antenna patterns, temperature correction, Galactic calibration, noise subtraction, and the GeoCeLDF lateral-distribution fit), and describes the simulation setup with GDAS atmospheres and measured noise injection. It does not present reconstructed E_rad values, a data-simulation comparison, or a final systematic uncertainty budget; the central claim is that this method will allow an independent cross-check of the FD energy scale.
Significance. If the proposed cross-check works, it would provide a valuable test of the FD-based cosmic-ray energy scale using a radio technique that does not rely on atmospheric fluorescence propagation and that can operate around the clock. The paper has genuine strengths: a large and carefully selected hybrid dataset (912 events), a realistic simulation setup with GDAS atmospheres and event-by-event measured noise injection, identical reconstruction for data and simulations, and an explicit status list for each ingredient. These are appropriate building blocks for a future result. The significance is currently prospective rather than demonstrated: no quantitative comparison is shown, so the paper's contribution at this stage is the validated methodology and dataset, not the energy-scale cross-check itself.
major comments (4)
- [§1 and §6] The claim that the radio energy scale is 'independent' rests on CoREAS's absolute fluence predictions, but the only validation cited is code-to-code agreement with ZHAireS at about 3% on the energy scale (§6). This validates internal consistency between two simulations, not the absolute radio-emission efficiency against an FD-independent reference. In the proposed event-by-event ratio E_rad^data/E_rad^sim, a constant multiplicative offset in CoREAS fluence is exactly degenerate with an apparent FD-radio energy-scale discrepancy. The analysis needs either an explicit closure test that exercises the absolute normalization, or an explicit statement of this degeneracy and how it will be broken.
- [§4 and §6] No reconstructed E_rad values, no data-versus-simulation comparison, and no combined systematic uncertainty are presented, so the headline 'will allow us to make an independent cross-check' (§6) is a statement of intent. A simulated closure test (for example, reconstructing simulated events through the full chain and showing that E_rad^sim reproduces the input radio energy within the quoted uncertainties on a subset of events) would make the readiness claim concrete and would also quantify the impact of the planned Rice-based signal estimator.
- [§4] Using E_SD to set the simulation input energy puts E_rad^sim on the FD scale by construction. The paper should specify how the systematic uncertainty of E_SD, including the FD absolute scale and the invisible-energy corrections described in [22] versus the Sibyll 2.3d parameterization, propagates through the simulation into E_rad^sim. Without this propagation, the ratio E_rad^data/E_rad^sim cannot be assigned an uncertainty, and the comparison may be dominated by the very FD-scale systematics it is intended to test.
- [§3] The signal fluence estimator is explicitly preliminary: the text states that it is 'currently estimated using a simple noise subtraction approach' and that a Rice-distribution-based method is planned. Since the per-station fluence is the input to the GeoCeLDF fit and to E_rad, the quoted 4% radiation-energy resolution and the final ratio are not yet fixed. Please state the expected impact of the planned estimator change on both the resolution and the systematic uncertainty.
minor comments (6)
- [§2] The sentence about the AugerPrime Radio Detector ends with a placeholder citation '[? ]'; this reference needs to be completed.
- [Fig. 2] The red/green color coding described in the caption is not visible in grayscale, and the diagram's flow is hard to follow; consider a labeled schematic with arrows that are explained in the text.
- [§3] The phrase 'validated successfully with a resolution of the radiation energy of 4%' should specify whether this is the per-event resolution of the integrated E_rad, over which energy and zenith range, and whether it refers to the GeoCeLDF fit alone or to the full reconstruction chain.
- [§5] It would be useful to report the number of events removed by each quality cut (signal-station requirement, GeoCeLDF fit success, thunderstorm flag) so that the 912-event sample is reproducible from the raw hybrid set.
- [§6] The bullet stating 'a systematic uncertainty on the CR energy scale of 6%' is ambiguous; it should say explicitly that this is the propagation of the Galactic-calibration uncertainty into E_rad and then into the primary energy scale.
- [References] Reference [12] is formatted as 'arXiv. [2407.18654]' with an odd period between 'arXiv' and the identifier; the formatting should be harmonized with the journal style.
Circularity Check
No circular step: AERA-measured E_rad is compared against CoREAS-predicted E_rad driven by FD-calibrated E_SD, so the test is a genuine cross-check; the only weakness is the lack of an FD-independent absolute normalization for CoREAS, which is a validation caveat, not a definitional circle.
full rationale
The paper's derivation chain compares an absolutely calibrated AERA measurement E_rad^data with a CoREAS simulation prediction E_rad^sim whose input energy is the FD-calibrated SD energy E_SD (Section 4). These two quantities are not defined in terms of one another: E_rad^data comes from the detector calibration chain (lab-measured electronics, drone-rescaled antenna patterns, Galactic calibration, Section 3), while E_rad^sim comes from a first-principles simulation of radio emission from an air shower with prescribed energy E_SD. The ratio tests whether the FD-scale input energy and the radio-emission prediction are mutually consistent. That is an independent cross-check, not a fitted prediction. The paper's own caveat in Section 6 — that the code-to-code agreement between CoREAS and ZHAireS at the 3% level is important because the simulation 'is at the foundations of the radio energy scale' — identifies a genuine limitation: the absolute fluence of CoREAS has not been experimentally calibrated against an FD-independent reference, so a constant offset in simulated fluence would be absorbed as an apparent FD-radio scale mismatch. But this is a scientific validation concern, not circularity by construction. The many Auger self-citations are for detector calibration, software, and the simulation codes; none of them simply asserts the target result, and the ZHAireS comparison provides an independent code benchmark. Therefore no circular step is present; the slight score elevation reflects the mild independence caveat rather than a definitional loop.
Assumptions & free parameters
free parameters (1)
- Per-event GeoCeLDF parameters =
not quoted (one fit per event)
assumptions (6)
- domain assumption CoREAS predicts absolute radio emission at the AERA band with sufficient accuracy for a few-percent-level energy-scale cross-check.
- domain assumption The SD energy estimator E_SD is an unbiased proxy for the FD-scale primary energy for the selected 912 events.
- domain assumption The invisible-energy parameterization for Sibyll 2.3d correctly converts total primary energy to calorimetric energy.
- domain assumption GDAS atmospheres with 3 h time resolution are adequate for the radio simulations.
- domain assumption The GeoCeLDF integral gives an unbiased estimate of radiation energy with 4 percent resolution.
- domain assumption Signal cleaning and processing biases affect data and simulation reconstructions in the same way.
Cite this review
Pith. "Pith review of Towards a Cosmic-Ray Energy Scale with the Auger Engineering Radio Array." pith.science (2026). https://pith.science/paper/G76KEVBE
@misc{pith2026241113387,
author = {Pith},
title = {Pith review of: Towards a Cosmic-Ray Energy Scale with the Auger Engineering Radio Array},
year = {2026},
howpublished = {\url{https://pith.science/paper/G76KEVBE}},
note = {Machine review of arXiv:2411.13387}
}
read the original abstract
Radio detection of cosmic-ray (CR) induced extensive air showers with digital antenna arrays is a matured technique by now. At the Pierre Auger Observatory, the Auger Engineering Radio Array (AERA) has been measuring air-shower signals in conjunction with the particle detectors of the surface detector (SD) for over ten years. For an absolute determination of the CR energy with the Auger baseline detectors, the shower size estimator from the SD is calibrated with the energy scale of the fluorescence detector (FD). However, AERA has an independent access to the energy scale through the reconstructed radio signals. The hybrid detectors at the Pierre Auger Observatory offer the unique opportunity to compare the two independent energy scales. In this contribution, we present our envisaged methodology for cross-checking the agreement between the energy scales of the FD and AERA using hybrid SD-AERA shower data and simulations. We show individual steps of our radio signal reconstruction and highlight the key ingredients for calibrated energy measurements.
Figures
Reference graph
Works this paper leans on
-
[22]
Pierre Augercollaboration, Phys. Rev. D100 (2019) 082003
work page 2019
- [1]
- [2]
- [3]
- [4]
-
[5]
Pierre Augercollaboration, Sci. Rev. - from the end of the world1 (2020) 8
work page 2020
- [6]
-
[7]
Pierre Augercollaboration, PoS ECRS (2023) 093
work page 2023
Show all 26 references
-
[8]
Abreu et al.,, Nucl
P. Abreu et al.,, Nucl. Instrum. Meth.A 635 (2011) 92
2011
-
[9]
Aab et al.,JINST 12 (2017) T10005
A. Aab et al.,JINST 12 (2017) T10005
2017
-
[10]
Reuzki et al.,PoS ARENA2024029
A. Reuzki et al.,PoS ARENA2024029
-
[11]
Correia dos Santos et al.,PoS ARENA2024030
D. Correia dos Santos et al.,PoS ARENA2024030
- [12]
-
[13]
Glaser et al.,Astrop
C. Glaser et al.,Astrop. Phys.104 (2019) 64
2019
-
[14]
Pierre Augercollaboration, Phys. Rev. Lett.116 (2016) 241101
2016
-
[15]
Pierre Augercollaboration, Phys. Rev. D93(2016) 122005
2016
-
[16]
Glaser et al.,JCAP 2016 (2016) 024
C. Glaser et al.,JCAP 2016 (2016) 024
2016
-
[17]
Pont, Ph.D
B. Pont, Ph.D. thesis, Radboud University, 2021.https://hdl.handle.net/2066/234170
2021
-
[18]
Gottowik, Ph.D
M. Gottowik, Ph.D. thesis, Wuppertal U., 2021.https://doi.org/10.25926/21hw-wq22
2021 doi
-
[19]
Abu-Zayyad et al.,Astrop
T. Abu-Zayyad et al.,Astrop. Phys.61(2015) 93
2015
-
[20]
Heck et al.,Tech
D. Heck et al.,Tech. Umw. Wis. B(1998) 6019
1998
-
[21]
Huege et al.,AIP Conf
T. Huege et al.,AIP Conf. Proc.1535 (2013) 128
2013
-
[23]
Riehn et al.,Phys
F. Riehn et al.,Phys. Rev. D102 (2020) 063002
2020
-
[24]
Kobal et al.,Astroparticle Physics15(2001) 259
M. Kobal et al.,Astroparticle Physics15(2001) 259
2001
- [25]
-
[26]
Mitra et al.,Astrop
P. Mitra et al.,Astrop. Phys.123(2020) 102470. 8
2020
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.