REVIEW 4 major objections 5 minor 50 references
Search for active-sterile neutrino transitions using Pierre Auger Observatory data
T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The Pierre Auger Observatory's non-observation of ultrahigh-energy neutrino candidates yields 90% confidence-level limits on the transition magnetic moment coupling active neutrinos to sterile neutrinos with masses between 1 TeV and 100 TeV
desk verdict First Auger exclusion of dipole-portal sterile neutrinos, undercut a bit by a flavor-branching oversimplification in the Earth-skimming channel. 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 mechanism is the dipole-portal upscattering process $\nu n \to N X$, driven by the transition magnetic moment $\mu_{\nu N}$ between an active neutrino and a heavy sterile state. In deep inelastic scattering on sea quarks, the cross section grows like $\mu_{\nu N}^2$ times a polynomial $F(\hat{s},t;m_N)$ dominated by small-$x$ parton distribution functions, and switches on above the kinematic threshold $E \gtrsim m_N^2/(2m_n)$. The dipole interaction also induces a $Z$-boson exchange term that adds about 15% to the photon-only result. The two detection channels then act as a differential diagnostic: down-going exposure scales with the enhanced cross section until the efficien
What would settle it
Measure the ultrahigh-energy neutrino rate in the down-going channel under an assumed flux and compare it with the Earth-skimming tau rate over the same exposure. If both rates rise together relative to Standard Model expectations, the excess is a flux effect; the dipole-portal interpretation specifically predicts a relative suppression of the Earth-skimming channel. A measured channel ratio consistent with the Standard Model while the total rate exceeds the prediction, or an event count in the down-going channel above the Poisson bound that the paper uses ($n_{\rm exp} > 2.44$), would conflic
Extended reading notes
Core claim
The paper's central claim is that the Pierre Auger Observatory, by detecting no ultrahigh-energy neutrino candidates under three benchmark astrophysical fluxes, excludes transition-magnetic-moment couplings $\mu_{\nu N}$ between active and heavy sterile neutrinos of mass $m_N \in [1,100]$ TeV at 90% confidence. The exclusion is driven by the dipole-portal upscattering process $\nu n \to N X$, whose cross section rises steeply above the kinematic threshold $E \gtrsim m_N^2/(2m_n)$ because it is dominated by deep-inelastic scattering on small-$x$ sea quarks; the produced $N$ then decays promptly to $\gamma\nu$ or $Z\nu$. The paper computes modified exposures for both down-going and Earth-skimm
Load-bearing premise
The limits rest on extrapolating the proton's small-$x$ parton distribution functions via standard DGLAP evolution down to Bjorken-$x \sim 10^{-12}$; if sea-quark densities saturate or deviate there, the upscattering cross sections and the derived bounds shift accordingly.
Editorial extensions
If this is right
- If the bounds are right, any dipole portal coupling in the 1–100 TeV sterile-mass window must lie below roughly $10^{-5}$–$10^{-1}\,\mu_B$ depending on the assumed flux, closing parameter space that collider mono-photon searches cannot reach.
- The predicted channel asymmetry—enhanced down-going and suppressed Earth-skimming rates—provides a practical way to tell an astrophysical flux excess from new physics in future data.
- Because the limits are flavor-independent, they extend to electron- and tau-neutrino channels and complement the muon-flavor-only reach of the existing event-based bound.
- With continued exposure, the null result either yields the first UHE neutrino detection or sharpens the exclusions as improved cosmic-ray composition measurements narrow the allowed flux range.
- The same non-observation can be re-derived under any future improved flux constraints, so the excluded region in the $(m_N,\mu_{\nu N})$ plane will become more sharply defined.
Reading between the lines
- A larger-exposure future detector could use the same down-going/Earth-skimming asymmetry as a generic diagnostic for any promptly decaying upscattering portal, since a flux excess raises both channels while a cross-section enhancement suppresses the Earth-skimming one.
- The small-$x$ parton distribution extrapolation is the main hidden lever: if sea-quark densities saturate instead of following the standard DGLAP extrapolation, the limits weaken most in the high-energy, low-$m_N$ corner; an EeV-scale measurement of the inclusive $\nu N$ cross section would resolve this.
- The flavor-independent coverage means this null result also constrains sources opaque to muon neutrinos, because electron- and tau-neutrino channels are bounded even when the muon channel is absorbed.
- Although the paper assumes prompt decay of $N$, the same data could be recast for long-lived or displaced-decay regimes, where the shower development and the down-going/Earth-skimming balance would differ.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper investigates magnetic-moment-induced transitions between active and heavy sterile neutrinos, focusing on the process ν n → N X at ultrahigh energies. The authors derive an effective dipole cross section including photon and Z-boson exchange (Appendix A), compute the Pierre Auger Observatory exposure in both down-going and Earth-skimming channels, and use the non-observation of UHE neutrino candidates together with three benchmark fluxes to derive 90% C.L. exclusion regions in the (m_N, μνN) plane for m_N between 1 TeV and 100 TeV. The central numerical claim is that, depending on the assumed flux, values of μνN above roughly 10^-5–10^-1 μ_B are excluded for sterile neutrino masses in that range.
Significance. If the result holds, this would be the first UHE-neutrino-based constraint on the dipole portal for sterile neutrino masses above 1 TeV, extending existing IceCube and LHC bounds. The paper has several concrete strengths: the cross-section derivation in Appendix A is explicit and reproduces the photon-only limit; the analysis involves no fitted parameters; the flux dependence is stated honestly; and the exposure estimates use the collaboration's established simulation framework. The main risks are the small-x PDF extrapolation, the flavor normalization in Eq. (11), and the ad hoc upper-bound criterion of Eq. (6). The Earth-skimming branching concern does not affect the central limit because, as written, Eq. (11) uses only the down-going exposure.
major comments (4)
- [Section II, paragraph after Fig. 1] The kinematic threshold is misquoted: the text states E ≳ m_N^2/(2m_n) ≈ 53 (m_N/10 TeV)^2 TeV. For m_N = 10 TeV this evaluates to 1e4 TeV^2/(2×0.938 GeV) ≈ 5.3e4 TeV, i.e. 53 PeV, not 53 TeV. The error is a factor of 1000 in the central energy-scale claim. Please correct the prefactor and confirm that the numerical cross sections and Fig. 4 used the exact kinematics of Eq. (1) rather than this approximate formula.
- [Section IV, Eq. (11)] The flux scenarios are introduced as '(single-flavor) flux' values, while the down-going exposure E_DG in Eq. (5) is described in Sec. III.A as accounting for the sum over neutrino flavors. If E_DG is an all-flavor exposure and Φν is per-flavor, the expected number of events in Eq. (11) is missing a factor of 3. If, instead, E_DG is already per-flavor, the wording in Sec. III.A is misleading. This is load-bearing: a factor of 3 in n_exp translates into a shift of roughly √3 in the μνN exclusion contours. Please clarify and correct the calculation if needed.
- [Section III.A, systematic-uncertainty paragraph] The BSM cross-section enhancement is driven by sea-quark densities at small Bjorken-x, but the PDFs are used down to x ∼ 10^-12 with only a DGLAP-based extrapolation. The quoted 5–9% systematic uncertainty from Ref. [41] does not include non-linear QCD/saturation effects or the absence of experimental constraints in this region. Since the central limits depend directly on this extrapolation, please provide an estimate of how the exclusion regions change under alternative small-x behavior, or explicitly state that the results are conditional on the DGLAP extrapolation. As written, the robustness claim in the introduction is stronger than the quoted uncertainty justifies.
- [Section III.A, Eq. (6)] The upper boundary of the excluded regions is determined by the ad hoc condition m_n^{-1} σνn X_ground(θ) ≥ 1, which sets the down-going efficiency to zero as a step function. This criterion is not derived from the detector simulation, yet it directly controls the upper edge of the contours in Fig. 4. Please justify this criterion with dedicated simulation evidence or perform a sensitivity test (e.g., varying the threshold by a factor of 2) to show how the upper boundary shifts. Without this, the upper part of the exclusion region is a modeling assumption rather than an observable result.
minor comments (5)
- [Section III.B, Eq. (8)–(10)] The Earth-skimming simulation treats the BSM interaction as NC-like assuming the produced heavy neutrino decays promptly with a ντ in the final state. This is inconsistent with the stated flavor-independent μνN, under which N would decay to νe, νμ, and ντ with roughly equal probability. The stress-test concern that this overestimates the Earth-skimming exposure does not affect the central limits, because Eq. (11) uses only the down-going exposure. Still, for consistency and for a quantitative interpretation of Fig. 3, the branching should either be included or the figure should be labeled as illustrative.
- [Section IV, Eq. (11)] The paper should state explicitly that the final 90% C.L. limits are based only on the down-going channel and that the Earth-skimming channel is not included in n_exp. The abstract and introduction emphasize both channels, so the reader should not have to infer this from Eq. (11).
- [Section III.A and Eq. (5)] Please provide the numerical inputs for the exposure calculation: the observation time ΔT, effective area, the particular LHAPDF PDF set used, and the integration ranges. Figures 2 and 3 give results in cm^2 sr s, but without these numbers the calculation is not reproducible from the text.
- [Figure 1 and caption] There is an inconsistency in the line-style description: the text says m_N = 10 TeV is solid and m_N = 30 TeV is dash-dotted, while the caption says dashed and plain lines. Please align them.
- [Eq. (1) and Appendix A] Please specify whether q_i(x,Q^2) denotes the sum of quark and antiquark PDFs and whether the sum in Eq. (1) runs over both. This affects the numerical value of the cross section and the interpretation of the small-x enhancement.
Circularity Check
No significant circularity: the exclusion limits follow from a null observation and externally specified flux inputs, with no fitted parameter being renamed as a prediction.
full rationale
The derivation chain is self-contained and non-circular. The paper assumes an effective dipole portal Lagrangian (Eq. A1) and computes the neutrino–nucleon upscattering cross section from it (Eq. 1, Appendix A). This cross section enters the down-going exposure (Eq. 5) and, together with three externally motivated flux scenarios (Section IV, citing Refs. [4], [8], [1,9,10,44]), yields the expected event count (Eq. 11). The 90% C.L. exclusion is then obtained by requiring n_exp > 2.44 against the observed zero neutrino candidates. No parameter is fitted to the neutrino data and then used to "predict" that same data; the μ_νN dependence is explicit (σ ∝ μ^2), so the limit is a direct inversion of the null count. Self-citations to prior Pierre Auger exposure and cosmic-ray papers are experimental inputs (detector response, observed spectra), not theoretical conclusions that presuppose the target bound. The only notable assumptions — DGLAP extrapolation to x ~ 10^-12 (Section III A) and the treatment of heavy-neutrino decay flavors in the Earth-skimming simulation (Section III B) — are modeling uncertainties, not circular reductions. Moreover, the final exclusion contours in Fig. 4 are based on Eq. (11), which uses only the down-going exposure E_DG, so the Earth-skimming flavor-branching issue does not enter the central limit. The paper is therefore assigned a circularity score of 0.
Assumptions & free parameters
free parameters (3)
- Phi_nu-low flux normalization =
E^2 Phi ~ 2e-10 GeV cm^-2 sr^-1 s^-1 at 100 PeV
- Phi_nu-high flux normalization =
E^2 Phi ~ 2e-9 GeV cm^-2 sr^-1 s^-1 at 100 PeV
- Phi_nu-with-UHE-p flux normalization =
E^2 Phi ~ 2e-9 GeV cm^-2 sr^-1 s^-1 constant from 100 PeV to 30 EeV
assumptions (6)
- domain assumption Effective dipole-portal operator of Eq. (A1) with only hypercharge coupling is a valid low-energy description of active-sterile transitions.
- domain assumption The sterile neutrino N decays promptly via N -> gamma nu and N -> Z nu with no significant mixing with active neutrinos.
- domain assumption PDFs from LHAPDF 6.5.3 with DGLAP evolution are valid down to Bjorken-x ~ 10^-12.
- domain assumption The three benchmark flux scenarios (Phi_nu-low, Phi_nu-high, Phi_nu-with-UHE-p) bracket the true astrophysical neutrino flux at ultrahigh energies.
- domain assumption The detection efficiencies from the collaboration's simulations (Refs. [28-30]) are reliable for both SM and BSM interactions.
- ad hoc to paper The condition m_n^-1 sigma_nuN X_ground(theta) >= 1 marks the loss of down-going detection efficiency.
Cite this review
Pith. "Pith review of Search for active-sterile neutrino transitions using Pierre Auger Observatory data." pith.science (2026). https://pith.science/paper/OFOSP2AL
@misc{pith2026260801496,
author = {Pith},
title = {Pith review of: Search for active-sterile neutrino transitions using Pierre Auger Observatory data},
year = {2026},
howpublished = {\url{https://pith.science/paper/OFOSP2AL}},
note = {Machine review of arXiv:2608.01496}
}
read the original abstract
We investigate the sensitivity of the Pierre Auger Observatory to physics beyond the Standard Model arising from magnetic-moment-induced transitions between active and heavy sterile neutrinos. Such dipole portal interactions can enhance neutrino-nucleon cross sections above a kinematic threshold set by the sterile neutrino mass, leading to observable modifications of neutrino detection rates at ultrahigh energies (UHE). We estimate the impact of these interactions on both down-going and Earth-skimming neutrino detection channels, the contrasting responses of which enable discrimination between an enhanced neutrino flux and a modified interaction cross section. Using the non-observation of UHE neutrino candidates, we derive neutrino-flux-dependent constraints with 90% confidence-level on the transition magnetic moment for sterile neutrino masses in the range 1 TeV-100 TeV. Under the assumed flux scenarios, the resulting flavor-independent limits extend existing bounds into previously unexplored parameter space.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[41]
Predictions for high energy neutrino cross-sections from the ZEUS global PDF fits
A. Cooper-Sarkar and S. Sarkar, JHEP01, 075 (2008), arXiv:0710.5303 [hep-ph]
work page Pith review arXiv 2008
-
[1]
Abdul Halimet al.(Pierre Auger), JCAP05, 024 (2023), arXiv:2211.02857 [astro-ph.HE]
A. Abdul Halimet al.(Pierre Auger), JCAP05, 024 (2023), arXiv:2211.02857 [astro-ph.HE]
arXiv 2023
-
[2]
R. Aloisio, V. Berezinsky, and P. Blasi, JCAP10, 020 (2014), arXiv:1312.7459 [astro-ph.HE]
arXiv 2014
-
[3]
D. Biehl, D. Boncioli, C. Lunardini, and W. Winter, Sci. Rep.8, 10828 (2018), arXiv:1711.03555 [astro-ph.HE]
work page Pith review arXiv 2018
-
[4]
R. Alves Batista, R. M. de Almeida, B. Lago, and K. Kotera, JCAP01, 002 (2019), arXiv:1806.10879 [astro-ph.HE]
arXiv 2019
-
[5]
D. Boncioli, D. Biehl, and W. Winter, Astrophys. J.872, 110 (2019), arXiv:1808.07481 [astro-ph.HE]
arXiv 2019
-
[6]
B. T. Zhang and K. Murase, Phys. Rev. D100, 103004 (2019), arXiv:1812.10289 [astro-ph.HE]
arXiv 2019
-
[7]
A. Condorelli, D. Boncioli, E. Peretti, and S. Petr- era, Phys. Rev. D107, 083009 (2023), arXiv:2209.08593 [astro-ph.HE]
work page Pith review arXiv 2023
Show all 50 references
-
[8]
M. S. Muzio and G. R. Farrar, Astrophys. J. Lett.942, L39 (2023), arXiv:2209.08068 [astro-ph.HE]
2023 arXiv
-
[9]
Rodrigues, J
X. Rodrigues, J. Heinze, A. Palladino, A. van Vliet, and W. Winter, Phys. Rev. Lett.126, 191101 (2021), arXiv:2003.08392 [astro-ph.HE]
2021 arXiv
-
[10]
M. S. Muzio, M. Unger, and S. Wissel, Phys. Rev. D107, 103030 (2023), arXiv:2303.04170 [astro-ph.HE]
2023 arXiv
-
[11]
B´ erat, A
C. B´ erat, A. Condorelli, O. Deligny, F. Mon- tanet, and Z. Torres, Astrophys. J.966, 186 (2024), arXiv:2402.04759 [astro-ph.HE]
2024 arXiv
-
[12]
Botneret al.(IceCube), Nuclear Physics B - Proceed- ings Supplements143, 367 (2005)
O. Botneret al.(IceCube), Nuclear Physics B - Proceed- ings Supplements143, 367 (2005)
2005
-
[13]
Aabet al.(Pierre Auger), Nucl
A. Aabet al.(Pierre Auger), Nucl. Instrum. Meth. A 798, 172 (2015), arXiv:1502.01323 [astro-ph.IM]
2015 arXiv
-
[14]
Aielloet al.(KM3NeT), Nature638, 376 (2025), [Er- ratum: Nature 640, E3 (2025)]
S. Aielloet al.(KM3NeT), Nature638, 376 (2025), [Er- ratum: Nature 640, E3 (2025)]
2025
-
[15]
Adrianiet al.(KM3NeT), Astrophys
O. Adrianiet al.(KM3NeT), Astrophys. J. Lett.984, L41 (2025), arXiv:2502.08508 [astro-ph.HE]
2025 arXiv
-
[16]
K. S. Babu, S. Jana, and M. Lindner, JHEP10, 040 (2020), arXiv:2007.04291 [hep-ph]
2020 arXiv
-
[17]
Kusenko and T
A. Kusenko and T. J. Weiler, Phys. Rev. Lett.88, 161101 (2002), arXiv:hep-ph/0106071
2002 arXiv
-
[18]
L. A. Anchordoqui, J. L. Feng, H. Goldberg, and A. D. Shapere, Phys. Rev. D65, 124027 (2002), arXiv:hep- ph/0112247
2002
-
[19]
L. A. Anchordoqui, H. Goldberg, D. Gora, T. Paul, M. Roth, S. Sarkar, and L. L. Winders, Phys. Rev. D 82, 043001 (2010), arXiv:1004.3190 [hep-ph]
2010 arXiv
-
[20]
K. S. Capelle, J. W. Cronin, G. Parente, and E. Zas, Astropart. Phys.8, 321 (1998), arXiv:astro-ph/9801313
1998 arXiv
-
[21]
Bertou, P
X. Bertou, P. Billoir, O. Deligny, C. Lachaud, and A. Letessier-Selvon, Astropart. Phys.17, 183 (2002), arXiv:astro-ph/0104452
2002 arXiv
-
[22]
Domokos and S
G. Domokos and S. Kovesi-Domokos, Phys. Rev. D55, 2526 (1997), arXiv:hep-ph/9603242
1997 arXiv
-
[23]
Buckley, J
A. Buckley, J. Ferrando, S. Lloyd, K. Nordstr¨ om, B. Page, M. R¨ ufenacht, M. Sch¨ onherr, and G. Watt, Eur. Phys. J. C75, 132 (2015), arXiv:1412.7420 [hep-ph]
2015 arXiv
-
[24]
Gandhi, C
R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic, As- tropart. Phys.5, 81 (1996), arXiv:hep-ph/9512364
1996 arXiv
-
[25]
Armesto, C
N. Armesto, C. Merino, G. Parente, and E. Zas, Phys. Rev. D77, 013001 (2008), arXiv:0709.4461 [hep-ph]
2008 arXiv
-
[26]
Connolly, R
A. Connolly, R. S. Thorne, and D. Waters, Phys. Rev. D 83, 113009 (2011), arXiv:1102.0691 [hep-ph]
2011 arXiv
-
[27]
Abbasiet al.(IceCube), Phys
R. Abbasiet al.(IceCube), Phys. Rev. D104, 022001 (2021), arXiv:2011.03560 [hep-ex]
2021
-
[28]
Aabet al.(Pierre Auger), JCAP10, 022 (2019), arXiv:1906.07422 [astro-ph.HE]
A. Aabet al.(Pierre Auger), JCAP10, 022 (2019), arXiv:1906.07422 [astro-ph.HE]
2019 arXiv
-
[29]
Abreuet al.(Pierre Auger), Phys
P. Abreuet al.(Pierre Auger), Phys. Rev. D84, 122005 (2011), [Erratum: Phys.Rev.D 84, 029902 (2011)], arXiv:1202.1493 [astro-ph.HE]
2011 arXiv
-
[30]
Aabet al.(Pierre Auger), Phys
A. Aabet al.(Pierre Auger), Phys. Rev. D91, 092008 (2015), arXiv:1504.05397 [astro-ph.HE]
2015 arXiv
-
[31]
Zhang and W
Y. Zhang and W. Liu, Phys. Rev. D107, 095031 (2023), arXiv:2301.06050 [hep-ph]
2023 arXiv
-
[32]
Brdar, Y.-Y
V. Brdar, Y.-Y. Li, S. R. Mir, and Y.-L. Wang, JHEP 10, 230 (2025), arXiv:2502.07024 [hep-ph]
2025
-
[33]
C. W. Bauer, N. L. Rodd, and B. R. Webber, JHEP06, 121 (2021), arXiv:2007.15001 [hep-ph]
2021 arXiv
-
[34]
Abrahamet al.(Pierre Auger), Phys
J. Abrahamet al.(Pierre Auger), Phys. Rev. D79, 102001 (2009), arXiv:0903.3385 [astro-ph.HE]
2009 arXiv
-
[35]
Corcella, I
G. Corcella, I. Knowles, G. Marchesini, S. Moretti, K. Odagiri, P. Richardson, M. Seymour, and B. Webber, JHEP01, 010, hep-ph/0011363
- [36]
-
[37]
S. J. Sciutto, (1999), arXiv:astro-ph/9911331
1999 arXiv
-
[38]
D. Heck, J. Knapp, J. Capdevielle, G. Schatz, and T. Thouw, CORSIKA: A Monte Carlo code to simulate extensive air showers (1998), forschungszentrum Karl- sruhe Report FZKA 6019
1998
-
[39]
Fletcher, T
R. Fletcher, T. Gaisser, P. Lipari, and T. Stanev, Phys. Rev. D50, 5710 (1994)
1994
-
[40]
Kalmykov, S
N. Kalmykov, S. Ostapchenko, and A. Pavlov, Nucl. Phys. B Proc. Suppl.52, 17 (1997)
1997
-
[42]
H. K. Gupta (ed.),Encyclopedia of solid Earth geo- physics, 2nd edition(Springer, 2021)
2021
-
[43]
S. I. Dutta, M. H. Reno, I. Sarcevic, and D. Seckel, Phys. Rev. D63, 094020 (2001), arXiv:hep-ph/0012350
2001 arXiv
-
[44]
Ehlert, A
D. Ehlert, A. van Vliet, F. Oikonomou, and W. Winter, JCAP02, 022, arXiv:2304.07321 [astro-ph.HE]
-
[45]
Aabet al.(Pierre Auger), (2016), arXiv:1604.03637 [astro-ph.IM]
A. Aabet al.(Pierre Auger), (2016), arXiv:1604.03637 [astro-ph.IM]
2016 arXiv
-
[46]
Huang, S
G.-Y. Huang, S. Jana, M. Lindner, and W. Rodejohann, Phys. Lett. B840, 137842 (2023), arXiv:2204.10347 [hep- ph]
2023 arXiv
-
[47]
Magill, R
G. Magill, R. Plestid, M. Pospelov, and Y.-D. Tsai, Phys. Rev. D98, 115015 (2018), arXiv:1803.03262 [hep-ph]
2018 arXiv
-
[48]
Mu˜ noz-Ovalle and S
A. Mu˜ noz-Ovalle and S. Vogl, Phys. Rev. D113, 055004 (2026), arXiv:2510.02450 [hep-ph]. 11
2026
-
[49]
Barducci and A
D. Barducci and A. Dondarini, JHEP10, 165 (2024), arXiv:2404.09609 [hep-ph]
2024 arXiv
-
[50]
O. G. Miranda, D. K. Papoulias, O. Sanders, M. T´ ortola, and J. W. F. Valle, JHEP12, 191 (2021), arXiv:2109.09545 [hep-ph]
2021 arXiv
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.