IceCube search for high-energy neutrinos produced in the precursor stages of gamma-ray bursts

Nick van Eijndhoven (for the IceCube Collaboration); Paul Coppin

REVIEW 3 major objections 6 minor 17 references

IceCube search for high-energy neutrinos produced in the precursor stages of gamma-ray bursts

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 GRB precursor neutrinos may be the real signal, and a new IceCube search targets them

desk verdict Useful public GRB precursor catalog and a plausible sensitivity projection, but the blanket 'more sensitive for all Δt' claim needs an overlap check before it can be trusted. read the letter →

arxiv 1908.06653 v1 pith:67J6WM5M submitted 2019-08-19 astro-ph.HE

Paul Coppin , Nick van Eijndhoven (for the IceCube Collaboration) This is my paper · ORCID

classification astro-ph.HE

keywords gamma-rayburstsneutrinoastronomyGRBprecursorsIceCubeFermi-GBMunbinnedlikelihoodchokedjettimewindow

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 argues that the precursor stage of a gamma-ray burst, and not the prompt flash, may be the dominant time for high-energy neutrino production, and that IceCube's earlier GRB searches looked in the wrong window. To enable a search in the right window, the authors build an up-to-date public GRB catalog, identify precursor flashes in Fermi-GBM data from 2011 through 2018, and design an unbinned maximum-likelihood IceCube analysis over those precursor time windows. They find 172 precursors in 145 of 1843 bursts, about 8 percent, with roughly half of the overlapping events matching an earlier Fermi-GBM precursor search. The central quantitative claim is that their proposed stacked search, using 1527 GRBs, is more sensitive to the average neutrino fluence from GRB precursors than the previous IC40/59 IceCube analysis for every considered time window size.

What carries the argument

The load-bearing machinery is a two-stage pipeline. Stage one defines the target: Fermi-GBM light curves are background-fitted with a linear model, rebinned into Bayesian blocks (variable-width histograms that optimize bin edges to separate signal from background), and any bin whose rate exceeds the background by 30 Hz in at least two NaI detectors, is separated from the prompt episode by at least 2 s of quiescence, and carries less than one third of the prompt fluence is tagged as a precursor. Stage two quantifies the neutrino sensitivity: an unbinned Poisson maximum-likelihood ratio with spatial, energy, and temporal probability densities and a flat signal time profile inside the window. The key simplification is that for a single GRB with $\Delta t < 1000$ s the expected number of background events in the window is $\ll 1$, so the search is background-free and the required signal is a fixed roughly 2.3 events, independent of $\Delta t$.

What would settle it

Re-run the same precursor-identification algorithm on randomized Fermi-GBM background intervals, shifting the light curves so no real burst is present, and count how often the 30 Hz, two-detector, 2 s quiescence rule flags a precursor; if the false-positive rate is comparable to the reported 8% occurrence rate, the precursor sample is dominated by background fluctuations and the neutrino windows are not tied to real emission. A second decisive check would compute the stacked IceCube sensitivity with the same 1527 bursts but replacing the precursor windows with randomly shifted windows of equal length; a similar sensitivity curve would show the time information carries no power.

Watch

Extended reading notes

Core claim

The paper claims that high-energy neutrinos from GRB precursors have not yet been searched for, and that a search is now feasible. It presents a new public GRB catalog, a Fermi-GBM precursor sample of 172 episodes in 145 bursts, and an unbinned likelihood analysis for IceCube that combines spatial, energy, and temporal information. For a single GRB with a time window shorter than 1000 s, the expected number of background events is far below one, so the search is essentially background-free and the fluence sensitivity does not depend on the window size. In a stacked analysis over 1527 GRBs, the authors show their projected sensitivity to the average neutrino fluence is better than that of the earlier IC40/59 analysis for all time windows considered, primarily because of the full detector geometry, improved event selection, and roughly five times more bursts.

Load-bearing premise

The neutrino search is only as reliable as the precursor definition: a precursor is whatever passes the 30 Hz, two-detector, at-least-2-s-quiescence, less-than-one-third-fluence cuts, and if that tag actually picks up background fluctuations or the beginning of the prompt emission, the neutrino time windows are placed at the wrong moments.

Editorial extensions

If this is right

For time windows shorter than 1000 s, the single-GRB search is essentially background-free, so a coincident neutrino would be significant and the fluence limit does not depend on how the window size is chosen.
The stacked search over 1527 GRBs would place average per-GRB fluence limits below the IC40/59 limits for every window size plotted, making it the strongest available constraint on GRB precursor neutrino emission.
The public catalog and precursor table let other instruments or follow-up programs reuse the same windows for independent searches.
Choked-jet and other precursor-emission models that predict neutrino fluence above the projected sensitivity can be directly confirmed or excluded by this analysis.

Reading between the lines

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

An implication left implicit is that the same window-selection logic can be applied to subthreshold GRB-like bursts already stored in the database; if precursors occur without a prompt trigger, those neutrinos would otherwise be missed.
Because the precursor selection flags only gamma-ray-bright flashes, any choked-jet precursor that produces neutrinos without a detectable gamma-ray excess is invisible to this search; a neutrino follow-up of all GBM triggers, not just flagged precursors, would test that population.
The flat time profile in the likelihood ignores the internal light-curve shape of each precursor; folding in the observed time structure could increase sensitivity for long windows while also testing the choked-jet picture, but this is a modification the paper does not make.
The same unbinned likelihood with per-event time windows could be applied to other transient classes with known pre-main-emission activity, such as repeating fast radio bursts, if precursor-like gamma-ray counterparts are identified.

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

The real deliverable here is the public catalog: 172 precursor episodes across 145 GRBs from 2011–2018, extracted from Fermi-GBM with a reproducible Bayesian-block procedure, plus the GRBweb database and parsing tools. That is a concrete, citable contribution that others in the field can use. The sensitivity projection for a stacked IceCube search is also a reasonable next step, and the paper is honest that no neutrino search has been run yet.

The methods are standard and clearly described. The precursor definition follows the previous GBM search, which is a sensible anchor. The background characterization and Bayesian-block segmentation are standard and do not look problematic. The single-GRB sensitivity argument for Δt < 1000 s is sound: with background essentially zero, the required number of signal events is independent of window size.

Now the soft spots. The precursor selection depends on thresholds—30 Hz excess, two-detector coincidence, 2 s quiescence, one-third fluence ratio—that are arbitrary and not varied. The paper reports 8% with no error bar and admits only half the precursors overlap with the prior search; that alone suggests the selection is not yet well validated. The sensitivity projection also relies on an unpublished IceCube event selection, so the absolute numbers cannot be checked.

The bigger issue is the stacked sensitivity at large Δt. The likelihood in Eq. (4.1) is written for a single GRB, and the text says data 'can be combined in a stacking analysis' without explaining how overlapping pre-prompt windows are handled. With 1527 GRBs over 8 years, at Δt = 10^5 s the windows frequently overlap; a single neutrino event could fall in multiple GRB time windows and be counted multiple times if the stacking is a simple product of per-GRB likelihood terms. That would artificially inflate the expected signal for a given fluence. The previous IC40/59 analysis used window sizes of seconds, so this problem was negligible there. The paper gives no overlap check or union-of-windows treatment, so the claim that the new analysis is more sensitive for all considered Δt is not yet supported.

For small Δt, the comparison is probably fine. For large Δt, it is genuinely uncertain. This is fixable: show how the stacked likelihood treats overlapping windows, run pseudo-experiments, and validate the precursor catalog against independent data.

Who should read this? GRB and neutrino folks who want the precursor catalog or plan to do a precursor search. The catalog itself is worth a citation even if the sensitivity claim needs work. I would send this to peer review, but only after the overlap issue is addressed. It is a conference proceedings, so the bar is lower, but the advertised sensitivity claim should not stand as is.

Referee Report

3 major / 6 minor

Summary. This ICRC 2019 proceedings paper describes work toward an IceCube search for neutrinos from GRB precursors. The authors construct a new, publicly available GRB catalog, analyze Fermi-GBM data to identify precursor emission episodes (reporting 172 precursors in 145 of 1843 bursts, about 8%), and develop an unbinned-likelihood analysis for coincident IceCube neutrinos. They present a sensitivity comparison with the earlier IC40/59 GRB-neutrino analysis, claiming that their proposed analysis is more sensitive for all time-window sizes Δt up to 10^5 s. The paper is a methods and projections contribution rather than a completed measurement.

Significance. If the sensitivity projection is correct, the proposed analysis would substantially improve on previous IceCube constraints on GRB precursor neutrino fluence, providing a direct test of choked-jet precursor models. The public GRB catalog and the automated precursor identification are potentially useful resources for the community. The paper is preliminary, however; the central sensitivity claim rests on an unspecified stacking procedure and an unpublished event selection, and the precursor selection is not quantitatively validated. These issues need to be addressed before the stated conclusions can be accepted.

major comments (3)

[Section 4, Eq. (4.1), Fig. 3] The stacked sensitivity projection is not sufficiently defined to support the claim that the analysis is 'more sensitive for all considered time window sizes Δt'. Equation (4.1) gives the likelihood for a single GRB, and the text only says the data 'can be combined in a stacking analysis' without specifying how overlapping time windows [tp−Δt, tp] are treated. For 1527 GRBs at a mean rate of about 0.5 per day, a window of Δt = 10^5 s (about 1.16 days) covers a substantial fraction of the total observing time, so a single neutrino event can fall in the windows of several GRBs. If the stacked likelihood is a product of per-GRB terms, that event is counted multiple times, artificially altering both the signal and background expectations. The authors should present the joint likelihood, explain how overlapping windows are handled (e.g., by using the union of windows or a time-dependent background), and recompute the sensitivity, or alternatively justify that overlap effects are negligible over the plotted range. Without this, the central sensitivity claim is unsupported.
[Section 3] The precursor identification relies on hard thresholds (30 Hz rate excess, at least two NaI detectors, 2 s quiescence, and a precursor-to-prompt fluence ratio below 1/3) that are adopted from a previous search, but no systematic variation of these thresholds is presented, and no uncertainty is attached to the reported 8% fraction (145/1843 GRBs). The only validation is the statement that 'roughly half' of the identified precursors were also seen in the previous Fermi-GBM search, which is qualitative and gives no overlap fraction with errors. Since the precursor-tagged analysis in Section 4 uses these episodes to define the neutrino search windows, a misclassification rate would directly affect the expected background and the final sensitivity. The authors should provide quantitative validation, such as injection tests or a detailed comparison with the catalog of [13], or explicitly state that the precursor sample is preliminary and not yet used for the sensitivity claims.
[Section 4, Fig. 2] The sensitivity calculation is not reproducible from the information given. The event selection is described only as 'an improved event selection' and the figures are labeled 'IceCube preliminary', with no reference to the underlying event sample, effective area, or background rate as a function of declination. In addition, the statement that λ≈2.3 corresponds to detecting a signal in 90% of background-free trials is a statistical claim that should be justified by a pseudo-experiment or an analytic calculation. Without these details, a reader cannot verify the claimed improvement over IC40/59. The authors should either provide a reference to a public description of the event selection or give the necessary detector response information in the paper.

minor comments (6)

[Section 3] The fraction 145/1843 is 7.9%; the text rounds to 'about 8%'. Please quote the fraction with a Poisson uncertainty (e.g., 7.9% ± 0.6%) or state explicitly that the value is preliminary and not yet part of a full measurement.
[Figure 1] The caption does not explain the color coding of the Bayesian block histogram or the relationship between the top and bottom panels. Please add a sentence describing the yellow precursor bins and red prompt bins, and note that both panels show the same GRB.
[Section 2] The database is said to contain 'over six thousand GRBs', but Section 3 uses a sample of 1843 Fermi-GBM bursts. Please clarify the difference between the full catalog and the sample used for the precursor search, and explain why the IceCube analysis uses the latter.
[Eq. (4.1)] The Poisson prefactor (n_s+n_b)^N / N! is written in a compact form that may confuse readers; consider writing it as e^{-(n_s+n_b)} (n_s+n_b)^N / N! to make the normalization explicit.
[Reference [9]] The author list 'K. Murase and other' is incomplete; please use the full list or 'et al.'.
[Abstract / Section 5] The abstract says the paper will 'discuss the implications in case a significant signal is found', but the text contains no such discussion. Either add a short outlook or adjust the abstract.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the sensitivity projection is self-contained and uses no fitted target quantity as input.

full rationale

The paper does not use the target quantity, precursor neutrino fluence, as an input to the analysis. The unbinned likelihood in Eq. (4.1) contains a fitted signal count, but the paper presents only a projected sensitivity; no neutrino data are fit and no observed flux is derived. The precursor-selection thresholds (30 Hz excess, 2 s quiescence, fluence ratio below one third) are adopted from the earlier external Fermi-GBM precursor search [13] and cross-checked against that search for overlapping years; this is calibration of an input catalog, not a definition of the neutrino result. IceCube detector response, event selection, and the IC40/59 comparison line come from prior IceCube publications [5, 17], which are instrument characterizations and prior search results rather than assumptions of the precursor-neutrino signal. The GRB catalog is assembled from external observatory catalogs and GCN circulars. The central claim, that the stacked search is more sensitive than IC40/59, rests on detector simulation, effective area, and a larger GRB sample, not on a circular fit or on renaming a known result. An unaddressed statistical-overlap issue for very large time windows would be a correctness or validity concern, not circularity.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The analysis relies on hand-chosen precursor cuts and on standard IceCube detector and likelihood inputs from prior publications; no new physical entities are introduced. The robustness of the 8% precursor fraction and of the sensitivity projection depends on these assumptions.

free parameters (5)

Rate excess threshold for signal bins = 30 Hz
Section 3: bins are tagged as signal if the average rate exceeds the background parameterization by 30 Hz; this cut is chosen without a reported optimization or systematic scan and directly determines the precursor sample.
NaI detector coincidence requirement = at least 2 detectors
Section 3: an excess must be seen in at least two NaI detectors to count as signal; this suppresses background fluctuations but also removes faint precursors in one detector.
Quiescent period between episodes = 2 s
Section 3: emission episodes separated by at least 2 seconds are treated as separate; this defines whether a precursor is a separate episode rather than part of a continuous rise.
Precursor-to-prompt fluence ratio cutoff = 1/3
Section 3: preceding episodes are tagged as precursors only if their fluence is less than one third of the prompt fluence; the paper notes previous searches used a similar cut, but the value remains an arbitrary threshold.
Background fit interval = [t-35 s, t-10 s]
Section 3: a linear fit on this interval predicts the background rate at time t; the choice of interval affects all excess determinations and therefore the entire precursor list.

assumptions (6)

domain assumption GBM background is approximately linear on time scales shorter than about 50 s.
Section 3: the linear fit in [t-35 s, t-10 s] and interpolation assume this; if the background is strongly curved, for example during spacecraft maneuvers or radiation belts, the excess bins are mis-estimated.
domain assumption The prompt emission is the episode with the largest fluence, and any earlier weaker episode is a precursor.
Section 3: this ranking defines the sample; a separate early flare from the same central engine or a later rebrightening would not fit this classification.
domain assumption Precursor phases produce enhanced hadronic interactions and neutrino emission relative to the prompt phase.
Section 1 and Section 5: motivates the search via choked jet models [8,9]; the paper presents no data supporting this model choice, and the null prompt-phase results do not constrain it.
domain assumption IceCube detector response, event selection, and atmospheric background model from prior IceCube publications apply to the proposed search.
Section 4: the sensitivity estimate relies on the effective area and filter described in [5,15,16], but no simulation is run for the precursor search in this paper.
standard math The unbinned maximum likelihood (Eq. 4.1) correctly accounts for signal and background with Poisson statistics and factorized PDFs.
Section 4: the likelihood is standard for IceCube point-source searches; it assumes known background rate nb and separable spatial, energy, and temporal terms.
domain assumption GCN circulars can be reliably parsed with regular expressions to build a uniform GRB catalog.
Section 2: the new catalog is assembled by regex parsing of hand-written GCN circulars plus machine-readable catalogs; no validation of parsing accuracy is reported.

how reviews work

0 comments

Cite this review

Pith. "Pith review of IceCube search for high-energy neutrinos produced in the precursor stages of gamma-ray bursts." pith.science (2026). https://pith.science/paper/67J6WM5M

@misc{pith2026190806653,
  author       = {Pith},
  title        = {Pith review of: IceCube search for high-energy neutrinos produced in the precursor stages of gamma-ray bursts},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/67J6WM5M}},
  note         = {Machine review of arXiv:1908.06653}
}

read the original abstract

While gamma-ray bursts (GRBs) were once believed to be a dominant source of astrophysical neutrinos, the lack of any significant correlation between high-energy photons and neutrinos has resulted in stringent limits on their neutrino fluxes. Previous IceCube searches for neutrinos from GRBs have generally focused on the prompt phase of GRBs. However, the higher density during the GRB precursor stage could lead to more hadronic interactions, and therefore be the predominant time of neutrino emission. We present results from an analysis of Fermi-GBM data to identify and characterize individual precursor flashes. Together with an up-to-date catalog of GRBs, we have made our results available via an online tool. An IceCube analysis method to search for neutrinos arriving coincident with GRB precursors has been developed. In this presentation, we compare the sensitivity of this analysis to previous IceCube analyses and discuss the implications in case a significant signal is found.

Figures

Figures reproduced from arXiv: 1908.06653 by the authors.

**Figure 1.** Top: Characterization of the GBM background rate for GRB trigger bn120308588 and detector n8. The running mean of the rate (black) is compared to the prediction from a fit to prior data points (orange). Regions in which the two distributions match are used to fit the background rate (red). Bottom: Bayesian block histogram of the same GRB (bn120308588), using the combined data from detector n4 and n8. Bins belonging … view at source ↗

**Figure 2.** Sensitivity of our analysis to the neutrino fluence of a single GRB as a function of its declination and valid for any time window ∆t < 1000 s. 10 1 10 0 10 1 10 2 10 3 10 4 10 5 t (s) 10 5 10 4 10 3 E 2 d N / d E ( G e V c m 2 ) IceCube preliminary Sensitivity to the average fluence This analysis, 1527 GRBs IC40/59 analysis, 299 GRBs [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

[13]

S. Zhu. PhD dissertation, University of Maryland, 2015. https://drum.lib.umd.edu/handle/1903/17258

work page 2015
[1]

Zhang, The Physics of Gamma-Ray Bursts

B. Zhang, The Physics of Gamma-Ray Bursts . Cambridge University Press, 2018

work page 2018
[2]

Barnes et al., ApJ 860 (2018) 38

J. Barnes et al., ApJ 860 (2018) 38

work page 2018
[3]

B. P. Abbott et al., ApJL 848 (2017) L13

work page 2017
[4]

A. M. Hillas, ARA&A 22 (1984) 425

work page 1984
[5]

IceCube Collaboration, ApJ 843 (2017) 112

work page 2017
[6]

Burlon et al., ApJL 685 (Sept., 2008) L19

D. Burlon et al., ApJL 685 (Sept., 2008) L19

work page 2008
[7]

Lazzati, MNRAS 357 (Feb, 2005) 722

D. Lazzati, MNRAS 357 (Feb, 2005) 722

work page 2005

Show all 17 references

[8]

Wang and P

X.-Y . Wang and P. Mészáros,ApJ 670 (Dec, 2007) 1247

2007
[9]

Murase and other, Phys

K. Murase and other, Phys. Rev. D 93 (2016) 083003

2016
[10]

Meegan et al., ApJ 702 (2009) 791

C. Meegan et al., ApJ 702 (2009) 791

2009
[11]

S. D. Barthelmy et al., SSR 120 (2005) 143

2005
[12]

Aguilar, ICRC (2011) 235

J.-A. Aguilar, ICRC (2011) 235

2011
[14]

J. D. Scargle et al., ApJ 764 (2013) 167

2013
[15]

IceCube Collaboration, ASR 62 (2018) 2902

2018
[16]

IceCube Collaboration, JINST 12 (2017) P03012

2017
[17]

IceCube Collaboration, Nature 484 (2012) 351. 8

2012

Pith tools

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

[1] [13]

S. Zhu. PhD dissertation, University of Maryland, 2015. https://drum.lib.umd.edu/handle/1903/17258

work page 2015

[2] [1]

Zhang, The Physics of Gamma-Ray Bursts

B. Zhang, The Physics of Gamma-Ray Bursts . Cambridge University Press, 2018

work page 2018

[3] [2]

Barnes et al., ApJ 860 (2018) 38

J. Barnes et al., ApJ 860 (2018) 38

work page 2018

[4] [3]

B. P. Abbott et al., ApJL 848 (2017) L13

work page 2017

[5] [4]

A. M. Hillas, ARA&A 22 (1984) 425

work page 1984

[6] [5]

IceCube Collaboration, ApJ 843 (2017) 112

work page 2017

[7] [6]

Burlon et al., ApJL 685 (Sept., 2008) L19

D. Burlon et al., ApJL 685 (Sept., 2008) L19

work page 2008

[8] [7]

Lazzati, MNRAS 357 (Feb, 2005) 722

D. Lazzati, MNRAS 357 (Feb, 2005) 722

work page 2005