Search for cascade events with Baikal-GVD

A.A. Doroshenko; A.G. Solovjev; A.I. Panfilov; A.N. Dyachok; A.P. Koshechkin; A.R. Gafarov; A.V. Avrorin; A.V. Korobchenko; A.V. Kozhin; A.V. Skurikhin

REVIEW 3 major objections 6 minor 8 references

Search for cascade events with Baikal-GVD

T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read A partially built Baikal-GVD already sets a neutrino flux upper limit only about three times above IceCube's.

desk verdict A genuinely new but still preliminary GVD cascade search: three candidates above 100 TeV and a limit within a factor of three of IceCube, with the paper itself admitting the atmospheric background estimate is not done yet. read the letter →

arxiv 1908.05430 v1 pith:KED6LQVF submitted 2019-08-15 astro-ph.HE astro-ph.IM

Baikal-GVD Collaboration: A.D. Avrorin , A.V. Avrorin , V.M. Aynutdinov , R. Bannash , I.A Belolaptikov , V.B. Brudanin , N.M. Budnev , G.V. Domogatsky

show 44 more authors

A.A. Doroshenko R. Dvornicky A.N. Dyachok Zh.-A.M. Dzhilkibaev L. Fajth S.V Fialkovsky A.R. Gafarov K.V. Golubkov N.S. Gorshkov T.I. Gress R. Ivanov K.G. Kebkal O.G. Kebkal E.V. Khramov M.M. Kolbin K.V. Konischev A.V. Korobchenko A.P. Koshechkin A.V. Kozhin M.V. Kruglov M.K. Kryukov V.F. Kulepov M.B. Milenin R.A. Mirgazov V. Nazari A.I. Panfilov D.P. Petukhov E.N. Pliskovsky M.I. Rozanov E.V. Rjabov V.D. Rushay G.B. Safronov B.A. Shaybonov M.D. Shelepov F. Simkovic A.V. Skurikhin A.G. Solovjev M.N. Sorokovikov I. Stekl O.V. Suvorova E.O. Sushenok V.A. Tabolenko B.A. Tarashansky S.A. Yakovlev

This is my paper · ORCID

classification astro-ph.HEastro-ph.IM

keywords neutrinoastronomyBaikal-GVDcascadeeventsastrophysicalneutrinosdiffusefluxupperlimitCherenkovdetectorhigh-energy

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

The paper's goal is to show that Baikal-GVD, a kilometer-scale neutrino telescope still under construction in Lake Baikal, can already search for astrophysical neutrinos through cascade events. Using 872 cluster-days of data recorded by one cluster in 2016 and three clusters in 2018, the analysis selected 417 cascade-like events, of which 18 had more than 20 hit optical modules and three had reconstructed energies above 100 TeV. Taking those three events as neutrino candidates, the collaboration derives a 90% confidence single-flavor upper limit on an $E^{-2.46}$ diffuse flux that is about three times higher than the flux measured by IceCube. The authors conclude that the detector is already in the same sensitivity ballpark as IceCube and, with more data and refined analysis, should identify the first astrophysical neutrinos. The result is explicitly preliminary: the probability that atmospheric muons or neutrinos produce the three events is still being calculated.

What carries the argument

The load-bearing mechanism is the cascade-event selection and two-step reconstruction chain. A cascade event is a neutrino interaction that produces a shower of particles and a roughly spherical Cherenkov light front, and the analysis isolates such events by demanding a high hit multiplicity, $N_{\mathrm{hit}}>20$, and reconstructed energy above 100 TeV, where atmospheric backgrounds are strongly suppressed. In the first reconstruction step the shower vertex is fit by a chi-squared minimum using photon arrival times under the assumption that the shower is a point-like light source; in the second step the direction and energy are fit by maximum likelihood using amplitude probabilities calibrated against Monte Carlo simulations of Cherenkov light propagation in water. This machinery matters because it lets a sparse, partially deployed array retain energy resolution of about 30%, a median direction resolution of about 4 degrees, and a vertex resolution of about 2 meters, making early-construction data astrophysically competitive.

What would settle it

Compute, via Monte Carlo simulation, the expected number of atmospheric muon and neutrino background events in 872 cluster-days passing the same $N_{\mathrm{hit}}>20$ and $E>100~\mathrm{TeV}$ selection; if that expectation is of order one or more, the three candidates are consistent with atmospheric background and the astrophysical interpretation fails.

Watch

Extended reading notes

Core claim

The central claim is that cascade-mode neutrino searches with a partially built Baikal-GVD array reach the astrophysical flux scale measured by IceCube. From 872 cluster-days collected in 2016 and 2018, the collaboration reconstructed 417 cascade-like events, restricted to 18 with $N_{\mathrm{hit}}>20$, and found three events with reconstructed energy above 100 TeV. Under the assumption of an $E^{-2.46}$ spectrum and single-flavor normalization, the three events are used to set a 90% confidence level upper limit on the diffuse astrophysical neutrino flux; the limit is only a factor of about three above the IceCube flux. The paper presents the three high-energy events as candidates for astrophysical neutrinos and states that more data and a refined analysis will allow the first astrophysical neutrinos to be identified.

Load-bearing premise

The result collapses if the three recorded high-energy events are actually atmospheric muons or atmospheric neutrinos rather than astrophysical neutrinos; the paper explicitly notes that the calculation of the probabilities for obtaining such high-multiplicity events from atmospheric backgrounds is still in progress.

Editorial extensions

If this is right

If the limit is correct, a partially built Baikal-GVD can independently probe the same diffuse astrophysical neutrino flux that IceCube discovered, from the Northern hemisphere.
With the planned first-stage array of nine clusters and more live time, the three candidates should grow into a statistically significant excess, allowing the first claimed astrophysical neutrinos at Baikal-GVD.
The 90% CL limit, being only a factor of three above the IceCube flux, provides a useful cross-check on the IceCube measurement in an energy range around 100 TeV.
More data will allow the atmospheric background probability to be quantified, turning the current candidates into either a confirmed signal or a background-dominated sample.

Reading between the lines

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

If the three candidates are genuine astrophysical cascades, their zenith distribution offers an early glimpse at the Northern-sky neutrino sky, complementing IceCube's largely Southern-sky sensitivity.
The sensitivity could be improved ahead of full construction by stronger rejection of atmospheric backgrounds, for example using the outer strings as a veto, rather than by waiting for more deployed volume.
A dedicated search combining 2019's five-cluster data with the 2016 and 2018 samples could test whether the three events cluster in a particular sky region or are isotropically distributed.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

What you should know: this is a real measurement, not a simulation. Using 872 cluster-days of Baikal-GVD data from 2016 and 2018, the collaboration finds three events with reconstructed energy above 100 TeV and Nhit > 20, and derives a 90% CL single-flavor E^-2.46 upper limit about three times the IceCube flux. The paper explicitly says the atmospheric muon/neutrino background probabilities for these three events are still being calculated. That is the key caveat, and the authors deserve credit for putting it in the text rather than burying it.

What is good: the data are new, the three candidate events are new, and this is the first cascade-channel search from a partially built GVD. The analysis pipeline is standard — vertex reconstruction by time minimization, energy/direction by maximum likelihood, then quality cuts — but it is applied to a new detector and gives an independent cross-check in the Northern Hemisphere. The spectral index is fixed to the IceCube value, no free parameter is fitted to force agreement, and the limit is a straightforward Feldman-Cousins count. There is no circularity here. The authors also openly label the result as preliminary, which is the right tone for an ICRC proceedings.

Where it is soft: the load-bearing gap is the missing atmospheric muon bundle background. The atmospheric neutrino background is small (about 0.08 events/yr above 100 TeV), but atmospheric muons dominate below 100 TeV and are not quantified at all. All three candidates sit close to the energy threshold (107, 153, 155 TeV), and only one is described as fully contained. A misreconstructed or partially contained muon bundle could easily populate that region. The quoted 30% energy resolution is also not propagated into the limit, though for an upper limit based on counts this likely would not change the conclusion dramatically. Minor points: the optical module count is inconsistent (1440 in the introduction vs 1140 in the conclusion), and references [4] and [5] are duplicated. None of these are fatal, but the background estimate is essential before the events can be called astrophysical candidates.

Who should read it: people working on neutrino telescopes, especially those tracking Baikal-GVD progress and early science. It is a status report, not a discovery claim, and it should be read as such. I would send it to peer review: the experimental result is worth refereeing, with the clear instruction that the background probabilities need to be added or the candidate interpretation softened. The paper is honest, the work is reproducible in principle, and the gap is stated by the authors themselves — that is exactly the kind of paper a good referee can help improve.

Referee Report

3 major / 6 minor

Summary. The paper reports a search for high-energy cascade neutrino events with the partially constructed Baikal-GVD detector, using 872 cluster-days of live time from one cluster in 2016 and three clusters in 2018. After a standard two-step reconstruction (vertex fit, then maximum-likelihood energy/direction fit) and quality cuts, the analysis selects 417 events with hit multiplicity Nhit > 13, of which 18 have Nhit > 20 and 3 have reconstructed energy above 100 TeV. The paper derives a 90% CL upper limit on a single-flavor E^-2.46 diffuse astrophysical neutrino flux using the Feldman-Cousins prescription [9], taking the three events as signal, and reports that the limit is only about a factor of three above the IceCube flux. The central interpretive claim is that Baikal-GVD, at partial construction, is already probing the astrophysical neutrino flux and that these three events are candidates for astrophysical neutrinos.

Significance. If the three events are indeed astrophysical neutrino cascades, this is an important early result: it would show that a modular, under-construction km^3-scale detector can approach IceCube sensitivity before completion. The analysis uses standard and reproducible tools: Monte Carlo simulation of signal and background, a well-defined reconstruction chain, and a conventional upper-limit prescription with an externally fixed spectral index. The paper is honest about its preliminary nature and explicitly states that the critical background probability calculation is still in progress. However, because that background calculation is the load-bearing element for the candidate interpretation, the significance of the result currently rests on an unquantified assumption rather than on demonstrated background rejection. The paper's value is therefore as an early sensitivity demonstration and a progress report, with the astrophysical interpretation not yet established.

major comments (3)

[Section 3, Table 2, Conclusion] The paper explicitly states that 'The calculations of the probabilities to obtain such high multiplicity events from atmospheric muons and neutrinos are in progress,' and it quantifies only the atmospheric-neutrino background (0.08 events/yr above 100 TeV), not atmospheric muon bundles, which dominate the event population below 100 TeV according to Fig. 2. Only one of the three selected events is described as contained, and the reconstructed energies are 107, 153, and 155 TeV, i.e., all three are just above the 100 TeV threshold. Without a quantitative estimate of the atmospheric muon-bundle background in this final sample, the identification of these three events as astrophysical neutrino candidates is not established. This identification is the load-bearing step for the paper's claim that Baikal-GVD is 'in the same ballpark' as IceCube and that these are 'candidates for events from astrophysical neutrinos.' The upper limit is conservative if the events are treated as signal, but the candidate interpretation requires the background to be negligible, which is precisely the unquantified point.
[Section 2.2, Table 2, Section 3] The paper quotes an energy resolution of about 30% for cascades but does not propagate this resolution into the selection or the limit. All three candidates are within a factor of 1.5 of the 100 TeV threshold (107, 153, and 155 TeV), so an atmospheric event with true energy below 100 TeV could migrate above the cut, while a genuine astrophysical neutrino at the threshold could migrate below. This migration matters because the expected atmospheric background rises steeply with decreasing energy, as shown in Fig. 2, and because the upper limit is based on the number of observed events above threshold. Please provide an estimate of the migration effect, either from unfolded MC distributions or from a smeared signal-plus-background model, or justify quantitatively that the effect is negligible for the reported limit.
[Section 3, Fig. 7 (right)] The 90% CL upper limit is derived from three observed events with no background subtraction and no systematic uncertainties assigned to the effective area, the energy scale, or the Monte Carlo simulation. The paper should state explicitly that the limit is a conservative upper limit under the assumption that all three events are signal, and that it is not an astrophysical flux constraint if a substantial fraction of the events is atmospheric. This clarification is necessary because the text moves from the limit to the claim that the detector is 'in the same ballpark' as IceCube without distinguishing a flux measurement from an upper limit with unquantified background.

minor comments (6)

[Section 3, 'Our limit is only a three higher...'] This sentence should read 'a factor of three higher than the IceCube flux'; the current phrasing is unclear and should be corrected.
[References] References [4] and [5] are identical (Astropart.Phys. 25, 140 (2006)); one of them is likely intended to be a different paper and should be corrected.
[Intro and Conclusion] The introduction says that four additional clusters commissioned in 2017-2019 bring the total to 1440 optical modules, while the conclusion says that five clusters with 1140 OMs are taking data since April 2019. These numbers are inconsistent and should be reconciled.
[Section 2.2, 'contained event'] The paper states that one of the three candidates is a contained event but does not define the containment criterion. Please specify how containment is determined (e.g., distance from cluster boundary or vertex position relative to instrumented volume).
[Figures 5 and 6] The right panels show gamma-ray source sky maps with 2-degree circles around the reconstructed event directions, but the figures lack clear axis labels and a legend for the color scale. The reader cannot tell whether the circles are centered on the sources or on the reconstructed directions, or what correlation is being claimed.
[Section 2.2, Eq. (1)] Equation (1) has a nonstandard layout with the normalization written as 1/(N_hit - 4) inside the sum; while the meaning is understandable, the notation should be cleaned up for publication.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the flux limit is derived from observed event counts with a standard statistical prescription and an externally fixed spectral index; IceCube is used only as a comparison benchmark.

full rationale

The paper's central result, the 90% CL single-flavor E^-2.46 flux limit, is obtained by applying the Feldman-Cousins prescription [9] to the three selected events with reconstructed energy above 100 TeV and Nhit > 20, as stated in Section 3: 'Taking into account the three recorded events with energies above 100 TeV and assuming an E^-2.46 spectrum single-flavor limit to the flux has been derived according to [9].' The input to this derivation is the observed event count (3 events in 872 cluster-days) plus an assumed spectral index; no parameter is fitted to reproduce the IceCube flux. The IceCube flux is invoked only as an external benchmark for comparison ('Our limit is only a three higher than the IceCube flux, so we are in the same ballpark'), not as an input that determines the limit. The expected signal and atmospheric-neutrino background rates in Section 2.2 use the IceCube normalization [6,8], but these expectations are not used to define or tune the three-event selection in a way that forces the resulting limit. The paper's own admission that 'The calculations of the probabilities to obtain such high multiplicity events from atmospheric muons and neutrinos are in progress' is a completeness or background-estimation limitation, not a circularity: it weakens the astrophysical-candidate interpretation but does not make any prediction equivalent to an input by construction. No self-citation is load-bearing for the limit derivation: references [1,4,5,7] supply detector description and prior reconstruction/selection methods, while the statistical procedure and IceCube flux are external. The finding is therefore that the paper is self-contained with respect to its flux-limit claim, with the caveat about unquantified atmospheric background belonging to correctness risk rather than circularity.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The analysis depends on the detector simulation, the reconstruction efficiency, and the Feldman-Cousins method. The only hand-chosen inputs are selection thresholds. No new particles, forces, or entities are introduced, and the IceCube spectral index is an external input used as a reference.

free parameters (3)

Hit multiplicity selection threshold = Nhit > 20
Chosen to suppress atmospheric muon bundles. The number of candidates and the resulting limit depend directly on this cut; it is an analysis choice rather than a fitted model parameter.
Reconstructed energy selection threshold = E > 100 TeV
Used to define astrophysical neutrino candidates because the expected atmospheric background is strongly suppressed above this energy. The threshold is stated before the data results, but its value affects the candidate count and limit.
Optical module hit charge threshold = > 1.5 ph.el.
Used for cascade vertex and energy reconstruction. It is a hand-chosen reconstruction cut that affects which events enter the analysis.

assumptions (3)

domain assumption Detector simulation accurately models light propagation in Baikal water and optical module response.
The effective areas and expected event rates in Fig. 1 rely on Monte Carlo simulation from this paper and prior work; no direct data-to-MC validation is shown in this preprint.
domain assumption Atmospheric backgrounds above 100 TeV are negligible or do not dominate the three candidates.
The paper asserts that background from atmospheric neutrinos is strongly suppressed above 100 TeV, but states that the candidate background probabilities are still in progress. The flux limit implicitly assumes signal origin for the three events.
standard math The Feldman-Cousins method is correctly applied to derive the 90% CL upper limit.
The limit is derived by citing reference [9]; no derivation is shown in the paper, and standard confidence-interval construction is assumed.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Search for cascade events with Baikal-GVD." pith.science (2026). https://pith.science/paper/KED6LQVF

@misc{pith2026190805430,
  author       = {Pith},
  title        = {Pith review of: Search for cascade events with Baikal-GVD},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KED6LQVF}},
  note         = {Machine review of arXiv:1908.05430}
}

read the original abstract

Baikal-GVD is a next generation, kilometer-scale neutrino telescope currently under construction in Lake Baikal. GVD is formed by multi-megaton sub-arrays (clusters) and is designed for the detection of astrophysical neutrino fluxes at energies from a few TeV up to 100 PeV. The design of the Baikal-GVD allows one to search for astrophysical neutrinos with flux values measured by IceCube already at early phases of the array construction. We present here preliminary results of the search for high-energy neutrinos via the cascade mode with the Baikal-GVD neutrino telescope.

Figures

Figures reproduced from arXiv: 1908.05430 by the authors.

**Figure 6.** Left: The events recorded in 21.08.2018 (see text). Right: The sky map of gamma-ray sources with E > 1 GeV as well as a 2 circle around the reconstructed positions of this event. 6 PoS(ICRC2019)873 [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

8 extracted references · 7 canonical work pages

[5]

Aynutdinov et al., Astropart.Phys

V . Aynutdinov et al., Astropart.Phys. 25, 140 (2006)

work page 2006
[9]

Feldman and R.D

G.J. Feldman and R.D. Cousins, Phys. Rev. D57, 3873 (1998). 8 PoS(ICRC2019)873

work page 1998
[1]

Aynutdinov et al., NIM A742 82-88 (2014)

V . Aynutdinov et al., NIM A742 82-88 (2014)

work page 2014
[2]

M. G. Aartsen et al., IceCube Coll., Science 342, 1242856 (2013), [arXiv:1311.5238]

arXiv 2013
[3]

I. Taboada, Talk at XXVIII International Conference on Neutrino Physics and Astrophysics, 4-9 June 2018, Heidelberg, Germany, DOI: 10.5281/zenodo.1286918, URL: https://doi.org/10.5281/zenodo.1286918

work page doi:10.5281/zenodo.1286918 2018
[6]

M.G.Aartsen et al. , Phys. Rev. Lett. 113 101101 (2014)

work page 2014
[7]

Avrorin et al., PoS (ICRC2017)962, (2017)

A.D. Avrorin et al., PoS (ICRC2017)962, (2017)

work page 2017
[8]

M. G. Aartsen et al., IceCube Coll., Phys. Rev. D 91, 022001 (2015)

work page 2015

Pith tools

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

[1] [5]

Aynutdinov et al., Astropart.Phys

V . Aynutdinov et al., Astropart.Phys. 25, 140 (2006)

work page 2006

[2] [9]

Feldman and R.D

G.J. Feldman and R.D. Cousins, Phys. Rev. D57, 3873 (1998). 8 PoS(ICRC2019)873

work page 1998

[3] [1]

Aynutdinov et al., NIM A742 82-88 (2014)

V . Aynutdinov et al., NIM A742 82-88 (2014)

work page 2014

[4] [2]

M. G. Aartsen et al., IceCube Coll., Science 342, 1242856 (2013), [arXiv:1311.5238]

arXiv 2013

[5] [3]

I. Taboada, Talk at XXVIII International Conference on Neutrino Physics and Astrophysics, 4-9 June 2018, Heidelberg, Germany, DOI: 10.5281/zenodo.1286918, URL: https://doi.org/10.5281/zenodo.1286918

work page doi:10.5281/zenodo.1286918 2018

[6] [6]

M.G.Aartsen et al. , Phys. Rev. Lett. 113 101101 (2014)

work page 2014

[7] [7]

Avrorin et al., PoS (ICRC2017)962, (2017)

A.D. Avrorin et al., PoS (ICRC2017)962, (2017)

work page 2017

[8] [8]

M. G. Aartsen et al., IceCube Coll., Phys. Rev. D 91, 022001 (2015)

work page 2015