REVIEW 2 major objections 5 minor 37 references
A new calculation of Earth-skimming very- and ultra-high energy tau neutrinos
T0 review · 2 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper calculates end-to-end how often Earth-skimming tau neutrinos yield detectable upward air showers for a POEMMA-like satellite, finding an effective aperture near 1 km² sr at $10^9$ GeV with 30–40% model uncertainty.
desk verdict A workmanlike proceedings calculation of POEMMA's tau-neutrino sensitivity that is genuinely new as a curve set, but the absolute aperture hangs on shower-detection modeling deferred to companion papers and excluded from the quoted uncertainty band. 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 object is the observation-probability differential $dP_{\rm obs} = ds'\, P_{\rm exit}(E_{\nu_\tau},\theta_{\rm tr})\, p_{\rm decay}(s')\, P_{\rm det}(E_{\nu_\tau},\theta_v,\theta_{\rm tr},s')$, integrated over the visible surface patch and the trigger solid angle to give the effective aperture $\langle A\Omega\rangle(E_{\nu_\tau}) = \int dP_{\rm obs}\, \hat{r}\cdot\hat{n}\, dS\, d\Omega_{\rm tr}$. Here $P_{\rm exit}$ encodes the neutrino charged-current conversion and tau energy loss in the Earth, $p_{\rm decay}$ is the tau decay probability along the path, and $P_{\rm det}$ encodes the Cherenkov geometry and photoelectron count: the tau decay must fall within the effective Cherenkov cone, $\theta_{\rm Ch}^{\rm eff}\sim 1.5^\circ$, and the shower, whose energy is taken as half the tau energy, must produce at least 10 photoelectrons in POEMMA's 2.5 m² effective optical aperture. The uncertainty band in the result is generated by switching among the ALLM and BDHM parameterizations of tau electromagnetic energy loss and the associated high-energy neutrino cross-section extrapolations. This machinery is what turns an assumed neutrino flux into a predicted rate of upward air showers.
What would settle it
Compare the photon-density and effective-Cherenkov-angle model against calibrated measurements of Cherenkov light from air showers of known energy and geometry; on the observational side, a detected diffuse flux of Earth-skimming tau neutrinos whose inferred rate disagrees with the predicted effective aperture by more than the quoted uncertainty band would falsify the sensitivity calculation.
Extended reading notes
Core claim
The paper's central claim is a new calculation of tau exit probabilities $P_{\rm exit}$ and detection probabilities $P_{\rm det}$ for Earth-skimming tau neutrinos, folded into an effective aperture and effective area for the POEMMA mission concept. With the standard charged-current neutrino interaction, tau propagation through the Earth with electromagnetic energy loss, tau decay in the atmosphere, and the Cherenkov photon density of the resulting extensive air shower, the effective aperture reaches around 1 km² sr at $E_{\nu_\tau}=10^9$ GeV for full 360° azimuthal coverage. Comparing the ALLM and BDHM tau energy-loss models and standard versus ALLM/BDHM-extrapolated neutrino cross sections gives a factor-of-two spread at the highest energies and a 30–40% uncertainty at the best energy; a top layer of water helps slightly at the highest energies, while rock gives more target nucleons at lower energies. Under five years of observation at 20% duty cycle, the full-coverage configuration is competitive with, and at the highest energies exceeds, the current IceCube/Auger/ANITA limits and the projected sensitivities of ground- and ice-based instruments, while the 30°-coverage configuration is not competitive for diffuse fluxes but remains the relevant mode for transient sources.
Load-bearing premise
The calculation leans on numerical models deferred to two companion papers: the air-shower Cherenkov photon density as a function of altitude and angle, the effective Cherenkov angle's dependence on shower parameters, and the approximation that the air-shower energy is half the tau energy; if those are inaccurate, the absolute sensitivity shifts outside the stated 30–40% band.
Editorial extensions
If this is right
- With full 360° azimuthal coverage and five years of observation, POEMMA's diffuse all-flavor sensitivity reaches below the current IceCube, Auger, and ANITA limits around $10^9$ GeV, so a null detection would begin to disfavor the higher cosmogenic flux models.
- With only 30° coverage, POEMMA would not compete in diffuse mode; the mission's case for diffuse neutrinos depends on expanding to full azimuthal coverage.
- For transient sources, the effective-area calculation gives concrete reach: long blazar flares modeled by RFGBW out to roughly 43 Mpc, tidal disruption events out to about 100 Mpc, and the extended-emission short gamma-ray burst model of KMMK out to about 117 Mpc for a single event.
- At the highest energies, the choice between ALLM and BDHM tau energy-loss models changes the expected rate by up to a factor of two, so the calculation identifies tau photonuclear energy loss as the input that most needs pinning down.
Reading between the lines
- The same $P_{\rm exit}\otimes P_{\rm det}$ machinery is detector-agnostic: if the air-shower and Cherenkov models are validated, any future satellite or high-altitude balloon can reuse the aperture integral with only its altitude, orbit, collection area, and threshold changed.
- Because the exit probability is directly proportional to the charged-current cross section at energies where attenuation is small, an observed rate of Earth-skimming tau showers could be inverted to constrain neutrino-nucleon cross sections around $10^8$–$10^{10}$ GeV, beyond direct accelerator reach.
- The paper's 30–40% uncertainty band is dominated by energy-loss and cross-section extrapolations, not by detector parameters, which suggests that improving those particle-physics inputs is currently more valuable than increasing the telescope's collection area.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a calculation of the acceptance of POEMMA to very- and ultra-high-energy Earth-skimming tau neutrinos. The authors write the diffuse effective aperture in Eq. (2.1) as an integral over the Earth's surface and tau emergence angles of the observation probability, which factorizes in Eq. (2.2) into tau exit probability, tau decay probability, and detection probability. They show the resulting effective aperture for POEMMA at 525 km altitude (Fig. 1), tau exit fluxes for several elevation angles (Fig. 2, left), diffuse all-flavor sensitivity curves (Fig. 2, right), and transient-source fluence sensitivities (Fig. 3). The main quantitative results are an effective aperture of order 1 km^2 sr near 10^9 GeV and a stated +/-30-40% sensitivity variation from neutrino cross-section and tau energy-loss model choices.
Significance. If the calculation is correct, this is a useful public quantification of POEMMA's tau-neutrino sensitivity, with an explicit breakdown of the dominant model uncertainties from neutrino cross sections and tau electromagnetic energy loss. The authors should be credited for a transparent formulation, for stating their approximations (theta_v approximately equal to theta_tr, E_shower = E_tau/2, theta_Ch^eff approximately 1.5 degrees), and for comparing against Auger, IceCube, ANITA, ARIANNA, ARA-37 and GRAND10k sensitivity estimates. The paper's principal value is as an independent cross-check of the POEMMA concept-study sensitivity; however, the absolute normalization rests on details of the air-shower detection model that are deferred to companion papers and are not yet validated at the level required by the quoted uncertainty band.
major comments (2)
- [Section 2, Eq. (2.3)] Eq. (2.3) as written is not consistent with the definition of dPobs in Eq. (2.2). In Eq. (2.2), dPobs is a differential probability in the decay path s' (dimensionless after integration), while in Eq. (2.3) it is multiplied by the disk area pi (v-s)^2 (theta_Ch^eff)^2 to obtain an effective area. Since the text states that Pdet already includes whether the detector lies within the Cherenkov cone, multiplying by the disk area appears to double-count the Cherenkov acceptance; if dPobs is meant to denote a different (e.g., per-unit-area) quantity in Eq. (2.3), that needs to be stated explicitly. Figure 3 and all transient-source fluence sensitivities are based on Eq. (2.3), so the authors should provide a derivation of the point-source effective area from Eq. (2.1) or otherwise clarify the formula.
- [Section 2 (Pdet paragraph) and final paragraph] The absolute value of the effective aperture and of the transient sensitivities depends on the air-shower detection model, which is not developed in this paper. The paragraph after Eq. (2.2) states that Pdet depends on the photon density as a function of altitude, theta_tr and shower energy; that the effective Cherenkov angle is 'on the order of 1.5 degrees, though its exact value depends on the altitude of the decay, theta_tr, and the shower energy [19]'; and that 'we have assumed that the shower energy is half of the tau energy.' These inputs enter multiplicatively through Pdet in Eq. (2.2), so a modest error in the photon-density model or in the shower-energy normalization would shift the central sensitivity outside the +/-30-40% band quoted for E_nu = 10^9 GeV, a band that covers only neutrino cross-section and tau energy-loss variations. The final paragraph's statement that 'Modeling of air showers, the impact of cloud cover and other variables are being reviewed and improved [37]' explicitly flags this part of the calculation as unfinished. The authors should either fold these systematics into the quoted uncertainty or validate the air-shower model (e.g., against full EAS simulations) before the absolute sensitivity curves are presented as quantitative predictions.
minor comments (5)
- [Section 3, second paragraph] The word 'intstrument' should be 'instrument'.
- [Author affiliations] The affiliation 'CRNS' should be 'CNRS'.
- [Eq. (2.1) and Eq. (2.2)] The notation dPobs is used as a differential in s' in Eq. (2.2), but Eq. (2.1) integrates it over dS and dOmega_tr; for clarity, define the observation probability density explicitly or indicate the integration variables in the definition of dPobs.
- [Fig. 1 (right)] The solid and dot-dashed curves are described in the text, but the figure itself has no legend; adding a legend or directly labeling the curves would improve readability.
- [Fig. 3 caption] The caption uses 'left figure' and 'right figure'; 'left panel' and 'right panel' would be clearer.
Circularity Check
No significant circularity: the sensitivity curves are computed outputs, not fitted inputs.
full rationale
The paper's derivation chain is a convolution of independently sourced physics inputs. Eq. (2.2) defines dPobs from P_exit (neutrino cross section and tau energy loss, with ALLM/BDHM and SM/ALLM/BDHM cross-section variants imported from external references), p_decay (standard tau decay kinematics), and P_det (POEMMA optical response parameters from Ref. [19]). The effective aperture in Eq. (2.1) and effective area in Eq. (2.3) are integrals over these inputs; they are outputs, not quantities used to set the inputs. The paper is transparent about the assumptions entering P_det: 'The effective Cherenkov angle theta_Ch^eff is on the order of 1.5 deg, though its exact value depends on the altitude of the decay, theta_tr, and the shower energy [19]' and 'In our evaluations below, we have assumed that the shower energy is half of the tau energy.' These are model inputs, possibly uncertain, but not circular: no target quantity is defined in terms of itself, and no fitted parameter is renamed as a prediction. The closing caveat ('Modeling of air showers, the impact of cloud cover and other variables are being reviewed and improved [37]') is a limitation statement, not a circularity. The quoted +/-30-40% uncertainty refers to cross-section and energy-loss variations; it does not indicate that the sensitivity was derived from itself. Refs. [19] and [22] are self-citations for technical details, but they are not used as an external uniqueness theorem or to forbid alternative choices, and the central claim remains an independent convolution of these inputs. Therefore no circularity step can be exhibited.
Assumptions & free parameters
free parameters (2)
- Shower energy fraction =
0.5 (E_shower = 0.5 * E_tau)
- Effective Cherenkov angle (theta_eff_Ch) =
~1.5 degrees
assumptions (5)
- domain assumption Neutrino-nucleon charged-current cross sections (sigma_SM from PDFs, plus sigma_ALLM and sigma_BDHM extrapolations) are correct at PeV-EeV energies.
- domain assumption Tau electromagnetic energy loss is described by the ALLM or BDHM structure-function parameterizations.
- domain assumption The Earth density profile is modeled with a 3 km water outer layer over rock of mass density 2.65 g/cm^3.
- domain assumption The tau decay length follows the standard gamma c tau = 5 km * E_tau / 10^8 GeV.
- domain assumption POEMMA detector response parameters (2.5 m^2 optical area, 0.2 quantum efficiency, 10 photoelectron threshold) represent the concept study design.
Cite this review
Pith. "Pith review of A new calculation of Earth-skimming very- and ultra-high energy tau neutrinos." pith.science (2026). https://pith.science/paper/RLQDWRQ4
@misc{pith2026190803603,
author = {Pith},
title = {Pith review of: A new calculation of Earth-skimming very- and ultra-high energy tau neutrinos},
year = {2026},
howpublished = {\url{https://pith.science/paper/RLQDWRQ4}},
note = {Machine review of arXiv:1908.03603}
}
read the original abstract
Cosmic neutrinos above a PeV are produced either within astrophysical sources or when ultra-high energy cosmic rays interact in transit through the cosmic background radiation. Detection of these neutrinos will be essential for understanding cosmic ray acceleration, composition and source evolution. By using the Earth as a tau neutrino converter for upward-going extensive air showers from tau decays, balloon-borne and space-based instruments can take advantage of a large volume and mass of the terrestrial neutrino target. The theoretical inputs and uncertainties in determining the tau lepton exit probabilities and their translation to detection acceptance will be discussed in the context of a new calculation we have performed. We quantify the experimental detection capability based on our calculation, including using the Probe of Extreme Multi-Messenger Astrophysics (POEMMA) concept study response parameters for optical air Cherenkov detection. These case studies are used to illustrate the features and uncertainties in upward tau air shower detection.
Figures
Reference graph
Works this paper leans on
-
[19]
M. H. Reno, J. F. Krizmanic and T. M. Venters, arXiv:1902.11287 [astro-ph.HE]
arXiv 1902
-
[37]
J. F. Krizmanic et al., PoS ICRC 2019, 936 (2019). 7
work page 2019
-
[1]
K. Kotera and A. V . Olinto, Ann. Rev. Astron. Astrophys.49, 119 (2011) [arXiv:1101.4256 [astro-ph.HE]]
arXiv 2011
-
[2]
L. A. Anchordoqui, Phys. Rep. 801, 1 (2019) [arXiv:1807.09645 [astro-ph.HE]]
arXiv 2019
-
[3]
M. G. Aartsen et al. [IceCube Collaboration], Phys. Rev. D 98, 062003 (2018) [arXiv:1807.01820 [astro-ph.HE]]
arXiv 2018
-
[4]
Sanguineti [ANTARES and KM3NeT Collaborations], Universe 5, 65 (2019)
M. Sanguineti [ANTARES and KM3NeT Collaborations], Universe 5, 65 (2019)
work page 2019
-
[5]
P. Allison et al., Astropart. Phys. 35, 457 (2012) [arXiv:1105.2854 [astro-ph.IM]]
arXiv 2012
-
[6]
S. W. Barwick et al. [ARIANNA Collaboration], Astropart. Phys. 70, 12 (2015) [arXiv:1410.7352 [astro-ph.HE]]
work page Pith review arXiv 2015
Show all 37 references
- [7]
-
[8]
Zas [Pierre Auger Collaboration], PoS ICRC 2017, 972 (2018)
E. Zas [Pierre Auger Collaboration], PoS ICRC 2017, 972 (2018)
2018
-
[9]
Abreu et al
P. Abreu et al. [Pierre Auger Collaboration], Astrophys. J. 755, L4 (2012) [arXiv:1210.3143 [astro-ph.HE]]
2012 arXiv
-
[10]
P. W. Gorham et al. [ANITA Collaboration], Phys. Rev. D99, no. 12, 122001 (2019) [arXiv:1902.04005 [astro-ph.HE]]
2019 arXiv
-
[11]
M. G. Aartsen et al. [IceCube, Fermi-LAT, MAGIC, AGILE, ASAS-SN, HAWC, H.E.S.S., INTEGRAL, Kanata, Kiso, Kapteyn, Liverpool Telescope, Subaru, Swift, NuSTAR, VERITAS and VLA/17B-403 Collaborations], Science 361, eaat1378 (2018) [arXiv:1807.08816 [astro-ph.HE]]
2018 arXiv
-
[12]
A. N. Otte, Phys. Rev. D 99, 083012 (2019) [arXiv:1811.09287 [astro-ph.IM]]
2019 arXiv
-
[13]
Fang et al., PoS ICRC 2017, 996 (2018) [arXiv:1708.05128 [astro-ph.IM]]
K. Fang et al., PoS ICRC 2017, 996 (2018) [arXiv:1708.05128 [astro-ph.IM]]. 6 Earth-skimming tau neutrinos Mary Hall Reno
2018 arXiv
-
[14]
Alvarez-Muñiz et al
J. Alvarez-Muñiz et al. [GRAND Collaboration], arXiv:1810.09994 [astro-ph.HE]
- [15]
-
[16]
Bertou, P
X. Bertou, P. Billoir, O. Deligny, C. Lachaud and A. Letessier-Selvon, Astropart. Phys. 17, 183 (2002) [astro-ph/0104452]
2002 arXiv
-
[17]
J. L. Feng, P. Fisher, F. Wilczek and T. M. Yu, Phys. Rev. Lett.88, 161102 (2002) [hep-ph/0105067]
2002 arXiv
-
[18]
Halzen and D
F. Halzen and D. Saltzberg, Phys. Rev. Lett. 81, 4305 (1998) [hep-ph/9804354]
1998 arXiv
-
[20]
A. V . Olinto et al., PoS ICRC 2017, 542 (2018) [arXiv:1708.07599 [astro-ph.IM]]
2018 arXiv
-
[21]
A. V . Olinto et al., PoS ICRC 2019, 378 (2019)
2019
-
[22]
T. M. Venters, M. H. Reno, J. F. Krizmanic, L. A. Anchordoqui, C. Guépin and A. V . Olinto, arXiv:1906.07209 [astro-ph.HE]
1906 arXiv
-
[23]
Guépin, F
C. Guépin, F. Sarazin, J. Krizmanic, J. Loerincs, A. Olinto and A. Piccone, JCAP 1903, 021 (2019) [arXiv:1812.07596 [astro-ph.IM]]
2019 arXiv
-
[24]
Motloch, N
P. Motloch, N. Hollon and P. Privitera, Astropart. Phys. 54, 40 (2014) [arXiv:1309.0561 [astro-ph.IM]]
2014 arXiv
-
[25]
Alvarez-Muñiz, W
J. Alvarez-Muñiz, W. R. Carvalho, A. L. Cummings, K. Payet, A. Romero-Wolf, H. Schoorlemmer and E. Zas, Phys. Rev. D 97, 023021 (2018) Erratum: [Phys. Rev. D 99, 069902 (2019)] [arXiv:1707.00334 [astro-ph.HE], arXiv:1901.08498 [astro-ph.HE]]
2018 arXiv
-
[26]
Y . S. Jeong, M. V . Luu, M. H. Reno and I. Sarcevic, Phys. Rev. D96, 043003 (2017) [arXiv:1704.00050 [hep-ph]]
2017 arXiv
-
[27]
Abramowicz, E
H. Abramowicz, E. M. Levin, A. Levy and U. Maor, Phys. Lett. B 269, 465 (1991)
1991
- [28]
-
[29]
M. M. Block, L. Durand and P. Ha, Phys. Rev. D 89, 094027 (2014) [arXiv:1404.4530 [hep-ph]]
2014 arXiv
-
[30]
Palomares-Ruiz, A
S. Palomares-Ruiz, A. Irimia and T. J. Weiler, Phys. Rev. D 73, 083003 (2006) [astro-ph/0512231]
2006 arXiv
-
[31]
Kotera, D
K. Kotera, D. Allard and A. V . Olinto, JCAP 1010, 013 (2010) [arXiv:1009.1382 [astro-ph.HE]]
2010 arXiv
-
[32]
Fang and B
K. Fang and B. D. Metzger, Astrophys. J. 849, 153 (2017) [arXiv:1707.04263 [astro-ph.HE]]
2017 arXiv
-
[33]
Rodrigues, A
X. Rodrigues, A. Fedynitch, S. Gao, D. Boncioli and W. Winter, Astrophys. J. 854, 54 (2018) d[arXiv:1711.02091 [astro-ph.HE]]
2018 arXiv
-
[34]
S. S. Kimura, K. Murase, P. Mészáros and K. Kiuchi, Astrophys. J. 848, L4 (2017) [arXiv:1708.07075 [astro-ph.HE]]
2017 arXiv
-
[35]
Albert et al
A. Albert et al. [ANTARES and IceCube and Pierre Auger and LIGO Scientific and Virgo Collaborations], Astrophys. J. 850, L35 (2017) [arXiv:1710.05839 [astro-ph.HE]]
2017
-
[36]
Lunardini and W
C. Lunardini and W. Winter, Phys. Rev. D 95, 123001 (2017) [arXiv:1612.03160 [astro-ph.HE]]
2017 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.