Recognition: 3 theorem links
· Lean TheoremSensitivity of the As-Built Askaryan Radio Array to Ultra-High Energy Neutrinos
Pith reviewed 2026-05-08 17:58 UTC · model grok-4.3
The pith
The Askaryan Radio Array reaches world-leading sensitivity to ultra-high energy neutrinos above 10^19 eV after ten years of operation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the as-built ARA, with its 2013-2023 exposure and an enhanced simulation pipeline that adds secondary particle production and realistic detector response, achieves world-leading sensitivity to ultra-high energy neutrinos above about 10^19 eV. Under the most optimistic flux models, up to 13 neutrinos are predicted to have been recorded at trigger level. Secondary interactions account for as much as 30 percent of the acceptance starting at 10^19 eV, and the work maps the expected rates and signatures of both multi-pulse events (from direct, refracted, and secondary pulses) and multi-station coincidences.
What carries the argument
The enhanced simulation pipeline that adds data-driven detector simulations and fully incorporates secondary particle production to compute the array-wide trigger-level acceptance.
If this is right
- ARA provides the tightest existing limits on the ultra-high energy neutrino flux above 10^19 eV.
- Up to 13 neutrinos may already have been recorded at trigger level, depending on the final event selection.
- Secondary particle production must be included in acceptance calculations because it contributes up to 30 percent of the total at the highest energies.
- Multi-pulse and multi-station signatures become observable and can be used to confirm neutrino candidates.
- The results directly guide the station spacing, trigger design, and analysis strategy for next-generation radio neutrino detectors such as IceCube-Gen2 Radio.
Where Pith is reading between the lines
- If no events are ultimately identified in the data, the flux models themselves would need revision rather than the detector performance.
- Other ice-based radio arrays could improve their own sensitivity forecasts by adopting similar secondary-inclusive simulations.
- Multi-station coincidences may enable crude source localization even with sparse arrays, opening a path to cross-correlation with other cosmic-ray or gamma-ray observatories.
- The 30 percent contribution from secondaries suggests that future designs might deliberately optimize for detecting these lower-energy secondaries to increase overall event rates.
Load-bearing premise
The simulation pipeline accurately models radio emission from both primary and secondary particles, the properties of the Antarctic ice, and the response of the as-built detector, while the external neutrino flux models used for rate predictions are realistic.
What would settle it
A search of the actual 2013-2023 ARA data that finds a number of neutrino candidates statistically inconsistent with the predicted range (zero to 13 events) under the optimistic flux models would directly test the sensitivity calculation.
Figures
read the original abstract
The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino observatory designed to detect the impulsive radio waves produced by relativistic particle cascades in the Antarctic glacial ice. Using a significantly enhanced simulation pipeline, which adds data-driven detector simulations and fully incorporates secondary particle production, we calculate the trigger-level acceptance of the entire array. We compare the resulting trigger-level sensitivity to constraints on the UHE neutrino flux from other detectors. Given its exposure from 2013 to 2023, we find that ARA achieves a world-leading sensitivity above about $10^{19}$ eV, depending on the details of the event selection used in a search. Moreover, we find that up to 13 neutrinos are predicted to have been observed in this period at trigger-level, assuming the most optimistic neutrino flux models. We show that observations of secondary particles account for up to 30\% of the total acceptance starting at $10^{19}$ eV, and we explore the potential signatures and implications of both multi-pulse (from direct and refracted pulses and/or from secondary particle interactions) and multi-station events. Finally, we comment on the implications of this study for the design of next-generation UHE neutrino experiments, in particular IceCube-Gen2 Radio.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a significantly enhanced Monte Carlo simulation pipeline for the Askaryan Radio Array (ARA) that incorporates data-driven detector simulations and full secondary particle production (muons, taus, and cascades). It computes the trigger-level acceptance for the as-built array over 2013–2023, compares the resulting sensitivity to existing UHE neutrino flux constraints from other experiments, and concludes that ARA achieves world-leading sensitivity above ~10^{19} eV. The work predicts up to 13 trigger-level events under the most optimistic flux models, attributes up to 30% of acceptance to secondaries at the highest energies, and discusses multi-pulse and multi-station signatures along with implications for IceCube-Gen2 Radio.
Significance. If the simulation accurately captures radio emission from primaries and secondaries, ice properties, and as-built detector response, the results meaningfully update ARA's exposure and highlight the non-negligible role of secondaries at E > 10^{19} eV. The quantitative acceptance curves, event-rate forecasts, and design implications for next-generation radio arrays constitute a useful contribution to the UHE neutrino field.
major comments (2)
- [Abstract and simulation pipeline] Abstract and simulation pipeline description: the headline claims of world-leading sensitivity above 10^{19} eV and up to 13 predicted events rest on the reported 30% secondary contribution to trigger-level acceptance. The manuscript states that the pipeline is 'significantly enhanced' and 'data-driven' but provides no cross-validation of the radio-emission modeling for secondary particles (muons, taus, cascades) against independent codes or in-situ calibration data; without such checks the acceptance (and therefore both sensitivity and event-rate) numbers could be systematically biased.
- [Results on event predictions] Results section on event predictions: the upper bound of 13 events is quoted only for the most optimistic neutrino flux models. The manuscript should quantify how the predicted event numbers and the 'world-leading' status change under more conservative flux assumptions (e.g., those already excluded by other experiments) to demonstrate robustness of the central claims.
minor comments (2)
- [Abstract] The abstract's qualifier 'depending on the details of the event selection used in a search' is imprecise; the main text should explicitly list the trigger and analysis selections for which the quoted sensitivity curves and event rates are computed.
- [Methods] Notation for acceptance, effective volume, and sensitivity should be defined consistently in the methods section to avoid potential confusion between trigger-level and analysis-level quantities.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive review. We address each major comment below and have incorporated revisions to improve clarity and robustness.
read point-by-point responses
-
Referee: [Abstract and simulation pipeline] Abstract and simulation pipeline description: the headline claims of world-leading sensitivity above 10^{19} eV and up to 13 predicted events rest on the reported 30% secondary contribution to trigger-level acceptance. The manuscript states that the pipeline is 'significantly enhanced' and 'data-driven' but provides no cross-validation of the radio-emission modeling for secondary particles (muons, taus, cascades) against independent codes or in-situ calibration data; without such checks the acceptance (and therefore both sensitivity and event-rate) numbers could be systematically biased.
Authors: We agree that explicit discussion of validation strengthens the claims. The radio-emission modeling for secondaries employs the same established formalisms (ZHS and Alvarez-Muñiz) and Monte Carlo frameworks used for primary cascades in prior ARA analyses and cross-checked against NuRadioMC and other independent codes. In the revised manuscript we add a dedicated paragraph in the simulation pipeline section that cites these validations, references in-situ ARA calibration data for ice properties, and notes the assumptions specific to secondary-induced cascades. This directly addresses the concern of potential systematic bias. revision: yes
-
Referee: [Results on event predictions] Results section on event predictions: the upper bound of 13 events is quoted only for the most optimistic neutrino flux models. The manuscript should quantify how the predicted event numbers and the 'world-leading' status change under more conservative flux assumptions (e.g., those already excluded by other experiments) to demonstrate robustness of the central claims.
Authors: We concur that robustness under varied flux assumptions is important. The revised manuscript includes an expanded table and accompanying text that reports predicted trigger-level event numbers for a representative set of flux models, ranging from the most optimistic (yielding up to 13 events) to more conservative models consistent with existing IceCube and ANITA limits (yielding 0–3 events). The world-leading sensitivity claim above ~10^{19} eV remains unchanged because it is driven by the acceptance curve itself rather than any single flux normalization; the added quantification shows the result is robust across the allowed parameter space. revision: yes
Circularity Check
No significant circularity: forward simulation from external models
full rationale
The paper derives trigger-level acceptance and event-rate predictions via forward Monte Carlo simulation of radio emission, ice properties, secondary particles, and as-built detector response, then folds in external neutrino flux models. No parameters are fitted to the target sensitivity or counts; no self-citation chain or uniqueness theorem reduces the central claims to inputs by construction. The derivation remains self-contained and externally benchmarkable.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Radio emission from electromagnetic cascades in ice follows the standard Askaryan mechanism with known coherence and angular distribution
Reference graph
Works this paper leans on
-
[1]
First, AraSim contains a native Monte Carlo generator of neutrino interactions
Event Generation AraSim supports three ways of generating events. First, AraSim contains a native Monte Carlo generator of neutrino interactions. Second, users may provide a list of particle interactions and locations to be simulated. Third, users may directly provide the electric field at each 9 https://github.com/ara-software/AraSim 10 https://github.co...
-
[2]
time domain
Askaryan Emission Model AraSim supports two Askaryan emission models. They are described in detail elsewhere [49], so we only briefly recount them here. The first is a “time domain” model, where the electric field is found by time-differentiating the vector potential for the particle shower. The vec- tor potential A(θ, t) under a far-field assumption has ...
-
[3]
refracted
Modeling the Environment After signals are generated, they must be propagated from the interaction vertex position to each of the sta- tion’s channels. The properties of the Antarctic glacial ice have important effects on the propagation of radio waves, namely refractive, attenuation, and birefringence effects, all of which AraSim can simulate. Below, we ...
-
[4]
traditional
Modeling the Detector After an electric field arrives at the station, AraSim must model the response of the detector. a. Array and Station Layout AraSim can be used to simulate up to 37 identical ARA stations in a hexagonal layout with user-defined inter-station spacing. This capability supported the orig- inal design studies for the array. Simulation of ...
2018
-
[5]
trees,” where each tree contains “branches
Event Readout AraSim saves all information using the ROOT file for- mat. ROOT organizes files into “trees,” where each tree contains “branches” that store individual variables or arrays. Each entry in a tree corresponds to a sin- gle event, allowing efficient, column-oriented access to large datasets without loading the entire file into mem- ory. This str...
-
[6]
Aiello et al
S. Aiello et al. (KM3NeT), Nature 638, 376 (2025), [Er- ratum: Nature 640, E3 (2025)]
2025
-
[7]
Greisen, Phys
K. Greisen, Phys. Rev. Lett. 16, 748 (1966)
1966
-
[8]
G. T. Zatsepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966)
1966
-
[9]
The Pierre Auger Cosmic Ray Observatory
A. Aab et al. (Pierre Auger), Nucl. Instrum. Meth. A 798, 172 (2015), arXiv:1502.01323 [astro-ph.IM]
work page Pith review arXiv 2015
-
[10]
R. Abbasi et al. ((IceCube Collaboration) §, IceCube), Phys. Rev. Lett. 135, 031001 (2025), arXiv:2502.01963 [astro-ph.HE]
- [11]
-
[12]
Barwick, D
S. Barwick, D. Besson, P. Gorham, and D. Saltzberg, J. Glaciol. 51, 231 (2005)
2005
-
[13]
G. A. Askar’yan, Zh. Eksp. Teor. Fiz. 41, 616 (1961)
1961
-
[14]
G. A. Askar’yan, Soviet Journal of Experimental and Theoretical Physics 21, 658 (1965)
1965
-
[15]
M. Muzio et al. (ARA), PoS ICRC2025, 1126 (2025), arXiv:2509.16414 [astro-ph.HE]
- [16]
-
[17]
Observation of In-ice Askaryan Radiation from High-Energy Cosmic Rays
N. Alden et al. (ARA), Phys. Rev. Lett. 136, 151003 (2026), arXiv:2510.21104 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[18]
L. D. Landau and I. Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 92, 535 (1953)
1953
-
[19]
A. B. Migdal, Phys. Rev. 103, 1811 (1956)
1956
-
[20]
J. Alvarez-Mu˜ niz, P. M. Hansen, A. Romero-Wolf, and E. Zas, Phys. Rev. D 101, 083005 (2020), arXiv:2003.09705 [astro-ph.HE]
-
[21]
A. Cummings, R. Krebs, S. Wissel, J. Alvarez-Mu˜ niz, W. R. Carvalho, Jr., A. Romero-Wolf, H. Schoorlem- mer, and E. Zas, Phys. Rev. D 111, 023012 (2025), arXiv:2311.03646 [hep-ph]
-
[23]
A. Coleman, O. Ericsson, C. Glaser, and M. Bustamante, Phys. Rev. D 110, 023044 (2024), arXiv:2402.02432 [astro-ph.HE]
-
[24]
Cummings, A
A. Cummings, A. Bishop, R. Krebs, and W. Luszczak, PoS ARENA2022, 014 (2023)
2023
-
[25]
A. Bishop et al. (ARA), PoS ICRC2023, 1169 (2023), arXiv:2308.07401 [astro-ph.HE]
- [26]
-
[27]
P. Allison et al. , Astropart. Phys. 35, 457 (2012), arXiv:1105.2854 [astro-ph.IM]
-
[28]
P. Allison et al. (ARA), Astropart. Phys. 70, 62 (2015), arXiv:1404.5285 [astro-ph.HE]
-
[29]
P. Allison et al., Nucl. Instrum. Meth. A 930, 112 (2019), arXiv:1809.04573 [astro-ph.IM]
-
[30]
A. M. Dziewonski and D. L. Anderson, Phys. Earth Planet. Interiors 25, 297 (1981)
1981
-
[31]
A. Connolly, R. S. Thorne, and D. Waters, Phys. Rev. D 83, 113009 (2011), arXiv:1102.0691 [hep-ph]
-
[32]
Kuhlmann, Nucl
S. Kuhlmann, Nucl. Phys. B Proc. Suppl. 79, 108 (1999)
1999
-
[33]
P. Allison et al. (ARA), Phys. Rev. D102, 043021 (2020), arXiv:1912.00987 [astro-ph.HE]
-
[34]
P. Allison et al. (ARA), Phys. Rev. D105, 122006 (2022), arXiv:2202.07080 [astro-ph.HE]
-
[35]
G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998), arXiv:physics/9711021
work page Pith review arXiv 1998
- [36]
-
[37]
M. Ahlers and F. Halzen, Phys. Rev. D 86, 083010 (2012), arXiv:1208.4181 [astro-ph.HE]
- [38]
-
[39]
R. Abbasi et al. (IceCube), Phys. Rev. D 113, 062002 (2026), arXiv:2507.22234 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[40]
Abdul Halim et al
A. Abdul Halim et al. (Pierre Auger), PoS ICRC2023, 1488 (2023)
2023
- [41]
-
[42]
C. Glaser et al. , Eur. Phys. J. C 80, 77 (2020), arXiv:1906.01670 [astro-ph.IM]
-
[43]
D. Garc´ ıa-Fern´ andez, A. Nelles, and C. Glaser, Phys. Rev. D 102, 083011 (2020), arXiv:2003.13442 [astro- ph.HE]
-
[44]
S. Agarwal et al. (RNO-G), JINST 20 (04), P04015, arXiv:2411.12922 [astro-ph.IM]
- [45]
-
[46]
Hokanson-Fasig (RNO), PoS ICRC2019, 913 (2020)
B. Hokanson-Fasig (RNO), PoS ICRC2019, 913 (2020)
2020
-
[47]
R. Abbasi et al. (IceCube-Gen2), PoS ICRC2021, 1183 (2021), arXiv:2107.08910 [astro-ph.HE]
- [48]
-
[49]
Brun and F
R. Brun and F. Rademakers, Nucl. Instrum. Meth. A 389, 81 (1997)
1997
-
[50]
A. Connolly, Phys. Rev. D 105, 123012 (2022), arXiv:2110.09015 [astro-ph.IM]
-
[51]
S. De Kockere, D. Van den Broeck, U. A. Latif, K. D. de Vries, N. van Eijndhoven, T. Huege, and S. Buitink, Phys. Rev. D 110, 023010 (2024), arXiv:2403.15358 [astro-ph.HE]
- [52]
- [53]
-
[54]
E. S. Hong, Searching for Ultra-high Energy Neutrinos with Data from a Prototype Station of the Askaryan Radio Array , Ph.D. thesis, The Ohio State University (2014). 22
2014
-
[55]
J. Alvarez-Muniz, W. R. Carvalho, Jr., and E. Zas, As- tropart. Phys. 35, 325 (2012), arXiv:1107.1189 [astro- ph.HE]
-
[56]
J. Alvarez-Muniz, A. Romero-Wolf, and E. Zas, Phys. Rev. D 84, 103003 (2011), arXiv:1106.6283 [astro- ph.HE]
-
[57]
J. Alvarez-Muniz, R. A. Vazquez, and E. Zas, Phys. Rev. D 62, 063001 (2000), arXiv:astro-ph/0003315
-
[58]
J. Alvarez-Muniz and E. Zas, Phys. Lett. B 434, 396 (1998), arXiv:astro-ph/9806098
-
[59]
J. Alvarez-Muniz and E. Zas, Phys. Lett. B 411, 218 (1997), arXiv:astro-ph/9706064
-
[60]
Salcedo Gomez, J
A. Salcedo Gomez, J. Flaherty, A. Connolly, et al. (ARA), PoS ARENA2024, 009 (2024)
2024
-
[61]
P. B. Price, O. V. Nagornov, R. Bay, D. Chirkin, Y. He, P. Miocinovic, A. Richards, K. Woschnagg, B. Koci, and V. Zagorodnov, Proceedings of the National Academy of Sciences 99, 7844 (2002)
2002
-
[62]
G. D. Q. Robin, S. Evans, and J. T. Bailey, Philosophical Transactions of the Royal Society of London, Series A: Mathematical and Physical Sciences 265, 437 (1969)
1969
- [63]
-
[64]
Kravchenko, D
I. Kravchenko, D. Besson, and J. Meyers, Journal of Glaciology 50, 522–532 (2004)
2004
-
[65]
V. V. Bogorodsky, C. R. Bentley, and P. Gudmand- sen, Radioglaciology (Springer Science & Business Media, 2012)
2012
-
[66]
D. E. Voigt, c-axis fabric of the south pole ice core, (spc14), U.S. Antarctic Program (USAP) Data Center (2017)
2017
- [67]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.