Pith. sign in

REVIEW 1 major objections 6 minor 11 references

Search for Ultra-High-Energy Neutrinos at the Pierre Auger Observatory: New Triggers, Methods, and Constraints

T0 review · 1 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read A blind seven-year search with two new electromagnetic triggers finds no ultra-high-energy neutrino candidates in Auger's 60–75 degree zenith band, improving the diffuse flux limit by 25% and point-source limits by up to 1.5-fold.

desk verdict Solid, incremental Auger proceedings with new DGL limits from EM triggers, but the blind-search claim underdocuments the training-sample exclusion. read the letter →

arxiv 2507.10214 v2 pith:QPKKNBWF submitted 2025-07-14 astro-ph.HE astro-ph.IM

Srijan Sehgal (for the Pierre Auger Collaboration) This is my paper
classification astro-ph.HEastro-ph.IM
keywords ultra-high-energyneutrinosPierreAugerObservatoryneutrino-inducedairshowerselectromagnetictriggersToTdtriggerMoPSdiffusefluxupperlimitspoint-sourceneutrino
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 reports a blind search for ultra-high-energy neutrinos in Pierre Auger Observatory data from 1 January 2014 to 31 December 2021, in the downward-going low-zenith band $60^\circ<\theta<75^\circ$. Two additional electromagnetic triggers, ToTd and MoPS, are folded into the event selection, together with a retrained Fisher discriminant built on the summed Area-over-Peak of the four earliest stations. No neutrino candidates survive in any of the five zenith sub-ranges. From this null result the paper derives a single-flavor 90% confidence upper limit on an $E_\nu^{-2}$ diffuse flux of $k_{90}^{\mathrm{DGL}}<1.2\times10^{-7}\,\mathrm{GeV}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ over $1.3\times10^{18}$ to $2.5\times10^{19.5}$ eV, a 25% improvement over the previous selection, plus declination-dependent point-source limits with up to a 1.5-fold gain at the most sensitive declinations. Because neutrinos are the only messengers that escape dense sources and travel undeflected, tightening these limits directly constrains models of cosmogenic and astrophysical neutrino production.

What carries the argument

The central objects are two hardware-level triggers, Time-over-Threshold-deconvolved (ToTd) and Multiplicity-of-Positive-Steps (MoPS), which catch the broad, low-amplitude, long-duration signals of electromagnetic cascades that the muon-tuned threshold and ToT triggers miss; they had been shown to boost neutrino-shower detection by a factor of 5–10 below $10^{19}$ eV. The neutrino selection adds a strict condition that 75% of the stations closest to the core carry a ToT, MoPS, or ToTd trigger, and uses a Fisher discriminant trained on five zenith sub-ranges, with the summed Area-over-Peak (ratio of integrated signal to peak, averaged over photomultipliers) of the four earliest stations as a key input. The exposure is computed from a hexagonal effective area of $A_{\mathrm{hex}}=1.95\,\mathrm{km}^2$ and the minute-by-minute count of active hexagons, folded with the neutrino-nucleon cross-section. Limits use the Feldman–Cousins prescription extended for systematic uncertainties, with $N_{90}=2.39$ and exposure systematics $[-13\%,+19.5\%]$.

What would settle it

Compare the event lists: retrieve the events in the 20% training subsample and check whether any of them pass the final Fisher cut and appear in the 2014–2021 search sample; if any do, the blind-search background estimate and the $N_{90}$-based limits would need to be recomputed.

Watch

Extended reading notes

Core claim

The central claim is that adding the ToTd and MoPS triggers to the surface-detector analysis, and retraining the neutrino identifier on the summed Area-over-Peak of the four earliest triggered stations, makes the Pierre Auger Observatory more sensitive to electromagnetic-rich, “young” air showers induced by neutrinos in the $60^\circ$–$75^\circ$ zenith range, especially below $10^{19}$ eV. Applying this selection as a blind search to seven years of data yields zero neutrino candidates. Assuming a 1:1:1 flavor ratio and an $E_\nu^{-2}$ spectrum, the same exposure implies $k_{90}^{\mathrm{DGL}}<1.2\times10^{-7}\,\mathrm{GeV}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ in the DGL channel alone, with the ToTd/MoPS triggers giving a 25% tighter limit and up to a factor-5 gain in exposure at lower energies, mainly through the $\nu_e$ charged-current channel. For point-like sources the new selection yields declination-dependent 90% limits up to 1.5 times stronger at the best declinations and extends coverage to declinations $\delta\lesssim -68^\circ$, which are accessible only through this channel in the Auger search program.

Load-bearing premise

The quoted background expectation of fewer than one event in twenty years—and hence the limit—relies on the 20% data sample used to train the Fisher discriminant being kept out of the candidate search; the paper does not state that this exclusion happened.

Editorial extensions

If this is right

  • The DGL diffuse limit tightens by about 25% relative to the earlier ToT+TH-only selection, to $k_{90}^{\mathrm{DGL}}<1.2\times10^{-7}\,\mathrm{GeV}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$.
  • In the energy range $1.3\times10^{18}$ to $2.5\times10^{19.5}$ eV, where about 90% of the expected event rate sits, the new EM triggers raise the exposure by up to a factor of 5 at low energies.
  • Point-source limits improve by up to 1.5-fold at the most sensitive declinations, and declinations $\delta\lesssim-68^\circ$ are probed only by this DGL channel.
  • The $\nu_e$ charged-current channel contributes about 85% of the total exposure, so the assumed 1:1:1 flavor ratio is a direct driver of the quoted limit.
  • The selection was tuned so the expected background is fewer than one event in twenty years, and zero candidates were observed, making the search effectively background-free at the quoted confidence.

Reading between the lines

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

  • A decisive re-analysis would be to train the Fisher discriminant on an explicitly disjoint subsample and publish the overlap check; that would settle whether the blind-search premise is fully satisfied.
  • Because ToTd and MoPS target electromagnetic-dominated showers, the same selection should extend naturally to ultra-high-energy photon searches, where the signal is similarly EM-rich; this is a cross-check the paper does not run.
  • The DGL band's exclusive access to southern declinations complements the northern-sky sensitivity of IceCube, so a joint declination-dependent limit combining both experiments should be tighter than either alone.
  • The factor-5 low-energy exposure gain suggests the same triggers could lower the energy threshold for real-time neutrino alerts for multi-messenger follow-up, a use the paper does not discuss.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 6 minor

Summary. This paper, from the Pierre Auger Collaboration, describes a search for ultra-high-energy neutrino-induced air showers in the downward-going low-zenith (DGL) range 60° < θ < 75° using the Surface Detector with the new Time-over-Threshold-deconvolved (ToTd) and Multiplicity-of-Positive-Steps (MoPS) triggers. The authors construct a Fisher-discriminant-based selection trained on 20% of the data, evaluate the detector exposure with a hexagonal-cell model, and apply the selection to data from 1 January 2014 to 31 December 2021. No neutrino candidates are found, yielding a single-flavor 90% C.L. integrated diffuse limit k_DGL90 < 1.2×10^-7 GeV cm^-2 s^-1 sr^-1 over an energy range stated as E_nu in [1.3×10^18, 2.5×10^19.5] eV, together with declination-dependent point-source limits that improve by up to a factor of 1.5 over the previous ToT+TH-only analysis at the most sensitive declinations.

Significance. The result is a defensible incremental contribution: it demonstrates that the ToTd and MoPS triggers improve the selection efficiency for electromagnetic-dominated low-energy neutrino showers, quantifies the gain in exposure (up to a factor of 5 at lower energies), and updates the Auger DGL constraints with a longer data set. The exposure formula and the explicit comparison against the ToT+TH-only configuration are strengths, and the non-detection is a falsifiable empirical statement. The paper does not provide machine-checked code or proofs, but the analysis is described in sufficient outline that the main limit calculation can be inspected. The overall significance to the field is moderate rather than transformative, and the central claim depends on the blind-search protocol being correctly implemented.

major comments (1)
  1. [Section 3.2] The blind-search claim is not adequately documented. The text states that the Fisher discriminant is trained on 20% of the data selected at random from the analysis period, and that the selection is then applied to Observatory data from 1 January 2014 to 31 December 2021, but it never states that the 20% training subsample is excluded from, or made statistically independent of, the final candidate search. If the same events are used both to set the Fisher threshold and to count candidates, the quoted background expectation of fewer than one event in 20 years is a property of the training set, not an independent prediction, and the use of N90 = 2.39 in Eq. (2) (a value appropriate only for zero expected background) would be invalid. This issue is load-bearing for the central non-detection claim and for both the diffuse limit in Section 4 and the point-source limits in Section 5. Please state explicitly whether the training sample was removed from the search sample (and, if applicable, that the search sample shown in Figure 2 is the complement), or provide a corrected background estimate and N90 if it was not removed.
minor comments (6)
  1. [Section 4] The energy range written as '[1.3×10^18−2.5×10^19.5] eV' uses a non-integer exponent in the upper bound; please clarify whether this is intentional and give the numerical value in standard scientific notation. The same notation appears in the abstract.
  2. [Equation (1)] The symbol m_N is not defined; please state that it is the nucleon mass and specify the units of A_hex and the resulting exposure units.
  3. [Figure 3 caption] The caption says 'The dashed lines are the exposures for this analysis, and the solid lines for the analysis performed with only ToT and TH triggers,' but the legend in panel (a) appears to use solid lines for 'All Triggers' and dashed for 'ToT+TH Only'; please reconcile the caption with the figure.
  4. [Reference [3]] The reference 'Nucl. Instrum. Meth. A 613 (2022) 29–39' appears to have an incorrect year; this is likely 2010. Please verify.
  5. [Section 3.1] The abstract quotes an improvement factor of 5–10 from the new triggers, while Section 4 reports 'up to a 5× increase' in exposure; please reconcile these numbers and specify whether they refer to reconstructed events or exposure.
  6. [Section 3.2] The phrase '75% of stations closest to the core must have ToT, MoPS, or ToTd triggers' is ambiguous for events with exactly four stations; clarify how the 75% requirement is applied when the number of stations is not divisible by four.

Circularity Check

1 steps flagged · score 6.0 of 10

The Fisher discriminant threshold is fitted to a 20% subset of the same search data, and the resulting 'expected background <1 event in 20 years' is used as a zero-background input to N90 in Eq. (2), without stating that the training subset is excluded from the final candidate search.

  1. fitted input called prediction [Section 3.2 (ν Selection) and Section 4, Eq. (2)]
    "The FDA is trained on 20% of the data selected at random from the analysis period to evaluate a detection threshold such that the expected background is fewer than one event in 20 years for the entire zenith angle range. ... Applying this selection, a search for neutrino-induced EASs was performed in the Observatory data between 1st January 2014 up to 31st December 2021. No neutrino candidates were found ... In this work, the Feldman–Cousins approach [6] extended to include systematic uncertainties [7, 8], is adopted yielding N90 = 2.39."

    As written, the Fisher cut is fitted to a random 20% subset of the same 2014–2021 data that is subsequently searched; the paper does not state that this training subset is removed or made statistically independent of the candidate search. The 'fewer than one event in 20 years' background is therefore imposed on the training sample by the choice of cut, not independently predicted for the search sample. Equation (2) then adopts N90 = 2.39, the zero-background Feldman–Cousins value, so the diffuse and point-source limits rest on a zero-background assumption that is partly guaranteed by construction rather than measured. If the training events remain in the search, the quoted background and limits are optimistically stated; if they are excluded, that exclusion is not documented.

full rationale

The paper's central non-detection and upper limits are not a pure tautology: the exposure in Eq. (1) is obtained from Monte Carlo simulations and detector-geometry counting, Eq. (2) is a standard limit inversion, and the null observation of candidates is a data result. The self-citations to prior Auger analyses [2,4] provide the baseline method, but the new trigger and exposure calculations are independent of those results and are not used to force the final limit. The one load-bearing circular step is the background construction in Section 3.2: the discriminant threshold is chosen on a 20% random subset of the analysis-period data so that the background is below one event in 20 years, and the same period is then searched with that threshold. Because the paper never states that the training subsample is excluded from or statistically independent of the search sample, the 'expected background < 1 event in 20 years' is a fitted property of the training data, not an a priori prediction. This fitted zero-background assumption is then inserted into Eq. (2) via N90 = 2.39, which is the Feldman–Cousins value for zero expected background, thereby making the quoted limits partly circular. The actual observed candidate count and the MC-based exposure retain independent content, so the circularity is partial rather than complete, but the blind-search justification for the limits is not documented. Score 6 reflects one 'prediction' (the background expectation) that reduces by construction, while the central limits still depend on independent exposure and observation inputs.

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

The central result is an experimental upper limit; converting observed event counts into flux limits relies on simulated detector efficiencies, assumed neutrino cross sections, an assumed spectral shape, and a statistical prescription. No new particles, forces, or entities are postulated.

free parameters (2)
  • Fisher discriminant cut value = 4.18 in the theta in (67.5, 70.5] degree sub-range; analogous thresholds in five sub-ranges
    Chosen so the expected background is fewer than one event in 20 years, using 20% of the analysis data. This is a data-tuned threshold, not a fixed physics constant.
  • Electromagnetic station fraction requirement = 75%
    Hand-chosen selection criterion: 75% of stations closest to the core must have ToT, MoPS, or ToTd triggers. This affects both signal efficiency and background rejection.
assumptions (7)
  • domain assumption Monte Carlo simulations of neutrino-induced showers and detector trigger and reconstruction response accurately model the efficiency epsilon_i,c(E_nu, X, theta).
    The exposure formula in Eq. (1) uses epsilon from simulations; any trigger-efficiency mis-modeling propagates directly into the reported limits.
  • domain assumption The effective detector geometry can be represented as hexagonal cells with A_hex = 1.95 km^2, with time-dependent count n_hex(t) capturing array changes.
    Used in Eq. (1) to convert live time into exposure; a coarse approximation of the triangular SD grid.
  • domain assumption Neutrino-nucleon cross sections sigma_i,c(E_nu) at ultra-high energies are known within the quoted theoretical uncertainty.
    The exposure integral folds in sigma; the quoted systematic uncertainty of [-13%, +19.5%] includes cross-section uncertainty.
  • domain assumption The flavor ratio at Earth is 1:1:1 for the total diffuse exposure.
    Summing channels assumes equal nu_e:nu_mu:nu_tau; oscillation averages make this standard, but it is not guaranteed for all source classes.
  • domain assumption The diffuse flux follows a power-law E^-2 spectrum for the integrated limit.
    The limit normalization k is defined for an E^-2 flux; different spectral indices would rescale the limit.
  • standard math Feldman-Cousins with systematic uncertainties yields N90 = 2.39 for zero background.
    Statistical method cited from [6,7,8]; provides the 90% C.L. conversion factor for zero observed candidates.
  • domain assumption Exposure is uniform in right ascension within plus or minus 0.6% for the point-source search.
    Assumed for the declination-dependent exposure calculation, citing [9].

how reviews work

0 comments
Cite this review

Pith. "Pith review of Search for Ultra-High-Energy Neutrinos at the Pierre Auger Observatory: New Triggers, Methods, and Constraints." pith.science (2026). https://pith.science/paper/QPKKNBWF

@misc{pith2026250710214,
  author       = {Pith},
  title        = {Pith review of: Search for Ultra-High-Energy Neutrinos at the Pierre Auger Observatory: New Triggers, Methods, and Constraints},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QPKKNBWF}},
  note         = {Machine review of arXiv:2507.10214}
}
abstract

The Pierre Auger Observatory has the capability to identify neutrino-induced extensive air showers above $10^{17}$ eV by using its large Surface Detector (SD) array. Data from the Observatory have been used to set some of the most stringent upper limits to the neutrino flux in the ultra-high energy (UHE) range. The data have also been used for follow-up detection of transient events in the context of multi-messenger astrophysics. In mid-2013, two additional SD triggers (Time-over-Threshold-deconvolved (ToTd) and Multiplicity-of-Positive Steps (MoPS)) were shown to increase the detection capability for the neutrino-induced air showers in the energy regime below $10^{19}$ eV by a factor of 5-10. This contribution will give an overview of the ongoing work regarding the searches for UHE neutrinos at the Pierre Auger Observatory. The impact of the ToTd and MoPS triggers for neutrino search in the zenith angle range of $60^{\circ} < \theta < 75^{\circ}$ is discussed. A novel neutrino identification method, which integrates these triggers, is applied to observational data to look for neutrino-like events using a $\textit{blind}$ search strategy. New constraints to point-like sources of UHE neutrinos will be presented for the angular range explored.

Figures

Figures reproduced from arXiv: 2507.10214 by the authors.

Figure 1
Figure 1. Reconstructed number of simulated neutrino events for all energies and channels (CC and NC) as a function of simulated zenith angle 𝜃MC for the sample with All (ToTd, MoPS + ToT, TH) triggers (purple bars) and only ToT+TH triggers (dark green bars). The bottom panel shows the ratio of the two samples. The increase in reconstructed events is most pronounced at energies below ∼ 1 EeV, largely due to the enhanced sensi… view at source ↗
Figure 2
Figure 2. Distribution of the Fisher variable after the DGL event selection for events with reconstructed zenith angle 𝜃rec ∈ [67.5 ◦ , 70.5 ◦ ]. The open histograms show the background training sample (purple) and the search sample (blue) and the filled histogram (orange) depicts the simulated DGL 𝜈 events. Events above the Fisher value indicated by the vertical red dashed line would be regarded as neutrino candidates popula… view at source ↗
Figure 3
Figure 3. Comparisons of exposure to UHE𝜈 in the DGL angular range for the time period 1 Jan 2014- 31 Dec 2021. The dashed lines are the exposures for this analysis, and the solid lines for the analysis performed with only ToT and TH triggers for the same time period. The single flavour 90% C.L. integrated limit gives 𝑘 DGL 90 < 1.2 × 10−7 GeV cm−2 s −1 sr−1 in the DGL channel only. It applies to an energy range of 𝐸𝜈 ∈ [1.3 … view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Comparison of the limits (1 Jan 2004-31 Dec 2013 with ToT+TH triggers and 1 Jan 2014-31 Dec 2021 with All triggers) to the current upper limits on the diffuse flux of UHE neutrinos. IceCube limits from [10] are scaled for a 𝐸 −2 𝜈 flux assumption. The predicted fluxes …
Figure 5
Figure 5. Figure 5: The upper limits (01.01.2004–31.12.2021) at 90% C.L. for different neutrino searches performed at the Pierre Auger Observatory of a single flavor point-like flux of UHE. The limit obtained in this analysis (purple) for the DGL channel is compared to the limits obtained…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

11 extracted references · 5 canonical work pages

  1. [1]

    Coleman et al.,Astropart

    A. Coleman et al.,Astropart. Phys.149 (2023) 102819 [arXiv:2205.05845]

    arXiv 2023

  2. [2]

    Aab et al

    A. Aab et al. [Pierre Auger Coll.],JCAP10(2019) 022 [arXiv:1906.07422]

    arXiv 2019

  3. [3]

    J.Abrahametal.[PierreAugerColl.], Nucl.Instrum.Meth.A 613(2022)29–39[ arXiv:1111.6764]

    work page Pith review arXiv 2022

  4. [4]

    J. L. Navarro Quirante,Thesis Univ. Granada(2013)

    work page 2013

  5. [5]

    R. A. Fisher,Annals Eugen.7 (1936) 179–188

    work page 1936

  6. [6]

    G. J. Feldman and R. D. Cousins,Phys. Rev. D57(1998) 3873–3889[arXiv:physics/9711021]

    arXiv 1998

  7. [7]

    Aab et al

    A. Aab et al. [Pierre Auger Coll.],Phys. Rev. D91 (9)(2015) 092008 [arXiv:1504.05397]

    arXiv 2015

  8. [8]

    Including Systematic Uncertainties in Confidence Interval Construction for Poisson Statistics

    J. Conrad, et al.,Phys. Rev. D67(2003) 092008 [arXiv:hep-ex/0202013]

    work page Pith review arXiv 2003

Show all 11 references
  1. [9]

    Aab et al

    A. Aab et al. [Pierre Auger Coll.],Science357 (6537)(2017) 1266–1270 [arXiv:1709.07321]

    2017 arXiv

  2. [10]

    Abbasi et al

    R. Abbasi et al. [IceCube Coll.], [arXiv:2502.01963]

  3. [11]

    Tor Vergata

    A. Aab et al. [Pierre Auger Coll.],JCAP11(2019) 004 [arXiv:1906.07419]. 8 Search for Ultra-High-Energy Neutrinos at the Pierre Auger Observatory with EM triggersSrijan Sehgal The Pierre Auger Collaboration A. Abdul Halim13, P. Abreu70, M. Aglietta53,51, I. Allekotte1, K. Almeid...

    2019 arXiv

Pith tools

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