Pith. sign in

REVIEW 2 major objections 8 minor 32 references

Search for gravitational waves associated with high-energy neutrinos detected by IceCube during the third observing run of LIGO-Virgo

T0 review · 2 major / 8 minor · reviewed 2026-07-31 · grok-4.5

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

Pith's one-line read No statistically significant gravitational-wave transient is found with IceCube high-energy neutrinos in LIGO-Virgo O3; the search sets 90% distance exclusion limits for several emission models.

desk verdict Clean O3×IceCat-1 null with usable model-dependent exclusion distances; incremental but correctly scoped and methodologically solid. read the letter →

arxiv 2607.28539 v1 pith:2WXV3ZJJ submitted 2026-07-30 astro-ph.HE gr-qc

Matthias Vereecken , Matteo Pracchia , Kara Merfeld , Iara Tosta e Melo , Wasim Javed , Patrick Sutton This is my paper · ORCID
classification astro-ph.HEgr-qc
keywords gravitationalwaveshigh-energyneutrinosIceCubeLIGO-VirgoO3multimessengerastronomyunmodeledsearchexclusiondistancesGRBprogenitors
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 asks whether any of the high-energy neutrino alerts IceCube issued during LIGO-Virgo’s third observing run were accompanied by a generic gravitational-wave transient. The authors run an unmodeled, direction- and time-targeted search that can pick up weaker signals than all-sky catalogs or real-time alerts, so it can catch coincidences ordinary neutrino follow-ups would miss. They find no candidate louder than background expectation across 23 Bronze and Gold alerts, and the distribution of p-values matches pure chance. From that null result they convert search sensitivity into lower bounds on how far away a source of each tested waveform family would have to lie. A sympathetic reader cares because a joint gravitational-wave plus high-energy-neutrino detection would be direct evidence of hadronic acceleration in GRB-like engines, including choked-jet systems that may be dim in gamma rays.

What carries the argument

X-Pipeline unmodeled targeted search: coherent excess-power analysis of multi-detector strain in a fixed [−500, +500] s on-source window and sky pixels compatible with each neutrino localization, ranked against time-slide background and tuned on injections of the tested waveform families.

What would settle it

A gravitational-wave candidate from one of these neutrino directions whose ranking statistic exceeds the loudest background trials at high significance, or an independently confirmed source lying well inside a quoted 90% exclusion distance that should have produced a recoverable injection-level signal.

Watch

Extended reading notes

Core claim

Across 15 Bronze and 8 Gold IceCat-1 track alerts that overlap usable O3 interferometer data, the loudest on-source gravitational-wave candidates are consistent with noise; the smallest p-value is 0.032 and the cumulative p-value distribution follows the null hypothesis. The paper therefore reports no evidence for associated gravitational-wave emission, either event-by-event or in the population, and quotes 90% exclusion distances for binary neutron-star, neutron-star–black-hole, accretion-disk-instability, and circular sine-Gaussian models that typically range from tens to a few hundred megaparsecs.

Load-bearing premise

The quoted distance limits assume any associated gravitational waves radiate the fixed energies and morphologies the authors inject, and that the emission falls inside a thousand-second window around the neutrino.

Editorial extensions

If this is right

  • If any of the analyzed neutrinos is astrophysical, a BNS-like gravitational-wave counterpart is excluded inside roughly 30–60 Mpc and an NSBH-like counterpart inside roughly 70–140 Mpc under the search assumptions.
  • Generic burst models with higher radiated energy or lower characteristic frequency are excluded to larger distances (up to a few hundred megaparsecs for the most favorable ADI and CSG cases).
  • Targeted unmodeled searches remain sensitive to weaker or morphologically unexpected signals that real-time and catalog-based neutrino follow-ups can miss.
  • The same pipeline can be reapplied to future higher-purity neutrino streams and improved localizations without changing the core method.

Reading between the lines

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

  • Cascade neutrino alerts, unavailable in O3 but active in later runs, would raise the average astrophysical purity of the sample and therefore tighten population-level statements even with a continued null.
  • Because the search is deliberately unmodeled, a future detection would immediately constrain engine physics (rotational instability, disk clumps, or compact-binary merger) rather than only confirming a pre-chosen template family.
  • Improved angular uncertainties in later IceCube catalogs would shrink the sky grids the pipeline must tile, reducing trial factors and modestly improving effective reach for the same detector network.
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

2 major / 8 minor

Summary. The manuscript reports an unmodeled, targeted search with X-Pipeline for generic gravitational-wave (GW) transients coincident with IceCube IceCat-1 track alerts during LIGO–Virgo O3. Data are analyzed in a fixed [−500, +500] s on-source window and over a sky grid derived from each neutrino localization; significance is assessed from off-source time-slide background trials of equal length, and 90% exclusion distances are obtained from injections of BNS, NSBH, accretion-disk instability (ADI), and circular sine-Gaussian (CSG) waveforms. No on-source candidate is significant (lowest p = 0.032 among 23 events; cumulative p-value distribution consistent with the null), and model-dependent 90% distance limits are tabulated (Tables 2–3, Fig. 2).

Significance. This is a clean, complementary multimessenger null result: it inverts the usual GW-triggered neutrino follow-up and can in principle recover weaker or morphologically unexpected GW emission associated with high-energy neutrinos. Methods follow established X-Pipeline targeted-search practice (coherent time–frequency maps, consistency tests, automatic cut tuning on the same families used for limits, injection-based 90% amplitude/distance conversion). The null association and the conditional Mpc limits for the stated emission models are useful reference results for O3 and for planning O4 cascade-alert and IceCat-2 analyses. Strengths include transparent conditioning of limits on explicit waveform families (Table 1, §3) and a population-level p-value check against the uniform null (Fig. 1).

major comments (2)
  1. [Abstract, §1, §5] Abstract and §1 claim the search is “sensitive to gravitational-wave signals weaker than those reported in real time or in the gravitational-wave transient catalog.” For generic bursts this is plausible (targeted coherent gain), but §5 explicitly states that the BNS/NSBH 90% exclusion distances (∼30–60 Mpc and ∼70–140 Mpc) are lower than O3 matched-filter single-detector BNS ranges. The abstract wording should be qualified so it does not imply a sensitivity gain over CBC catalog searches for inspiral-like signals; the correct selling point is the inverse, unmodeled, neutrino-triggered design that can catch non-CBC or sub-threshold burst emission missed by GW-first neutrino follow-ups.
  2. [§2, §4, Table 3] §4 circularizes IceCat-1 HEALPix 90% contours to a single Gaussian σ for the X-Pipeline sky grid, with only a brief consistency check against GCN angular errors. For several events σ ≳ 1–2° (Table 3), and IceCat maps can be non-Gaussian. Because coherent consistency tests are direction-dependent (§2), residual mismatch between the true localization and the circularized grid can both reduce efficiency and bias the loudest-event ranking. A short quantitative validation (e.g., injection recovery on a subset of events using the native HEALPix support, or a statement of the maximum efficiency loss) is needed to underwrite the tabulated distances for the larger-error alerts.
minor comments (8)
  1. [Abstract] Typo in abstract and elsewhere: “gravitional-wave” → “gravitational-wave”.
  2. [Title, Abstract] Title/affiliation block and running header say “LIGO-Virgo” while the abstract mentions KAGRA; O3 KAGRA data are not used. Align wording for consistency.
  3. [§3] Eq. (1)–(2): define ι and the inclination sampling cuts in the equation block or immediately adjacent text; the CSG inclination range (cos ι ∈ [0.996, 1]) appears only in prose.
  4. [Table 1] Table 1 caption says the accretion-disk mass is “always equal to 1.5 M⊙” but this is not a column; state it once in the table note and confirm it applies to all five ADI labels.
  5. [Fig. 2] Fig. 2 is hard to read in grayscale (many overlapping waveform classes). Consider separate panels per model or a colorblind-safe palette and explicit median markers matching Table 2.
  6. [§5] §5: “Seeing one such p-value in 23 analyses is well within the expected range” — a one-line binomial or look-elsewhere statement (e.g., expected number of p ≤ 0.032 under the null) would make the claim quantitative.
  7. [§1, References] References: GWTC-5 is cited as Abac et al. 2026 [arXiv:2605.27223]; ensure the bibliographic entry matches the version the journal will accept and that “version 5” wording remains accurate at publication.
  8. [References] IceCube Collaboration 2020 appears twice in the reference list with different contents (cascade-alert tech note vs GCN 27235); disambiguate the keys.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: standard targeted null search with external alerts, off-source background, and stated injection models.

full rationale

The paper’s central claims are an observational null (loudest on-source candidates consistent with time-slide off-source background; cumulative p-values match uniform) and conditional 90% exclusion distances obtained by injecting fixed waveform families (BNS/NSBH mass priors, CSG with EGW fixed to 10^{-2} M⊙c², five discrete van Putten ADI sets) into the same X-Pipeline pipeline. Significance is not fitted from the target quantity; model energies and the [-500,+500] s window are explicit inputs used only to convert amplitude limits into Mpc, not re-derived from the null. Citations to X-Pipeline, prior LVK GRB follow-ups, and ADI/CSG morphologies are methodological scaffolding, not load-bearing uniqueness theorems that force the result. The derivation chain is self-contained against external IceCube alerts and public GW strain; no step reduces by construction to its own inputs.

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

The null claim rests on standard GW data-analysis assumptions and external catalogs. The quantitative Mpc limits additionally rest on hand-chosen emission energies, discrete ADI parameter sets, a fixed on-source window, and circularized sky errors. No new physical entities are introduced.

free parameters (5)
  • EGW for CSG models = 10^-2 M⊙c²
    Fixed to 10^-2 M⊙c² as an ‘optimistic estimate’ to convert hrss limits into distance (§3, Eq. 2); directly sets CSG exclusion distances.
  • ADI model parameter sets (M, χ, ε) = ADI-A through ADI-E as in Table 1
    Five discrete combinations from van Putten models (Table 1) chosen to span morphologies; EGW and duration follow from those choices and drive ADI distance limits.
  • On-source time window = 1000 s total
    Fixed to [-500, +500] s around each neutrino (§2), motivated by GRB duration literature rather than fit to these data; defines what counts as a coincidence.
  • Search frequency band = 20–500 Hz
    20–500 Hz chosen by analogy to LVK GRB follow-up (§2); excludes higher-frequency burst content from the analysis.
  • CSG quality factor / inclination cuts = Q~9; cos ι in [0.996,1] (CSG), [0.866,1] (CBC)
    Q≈9 and tight face-on cos ι ranges for CSG and CBC injections (§3) are analysis choices that affect recovered sensitivity.
assumptions (5)
  • domain assumption Detector noise is uncorrelated across sites so time slides produce valid background trials for loudest-event p-values.
    Stated basis for off-source background enlargement in §2.
  • domain assumption If a neutrino alert is astrophysical and has an associated GW transient detectable by this search, emission occurs within the adopted on-source window and 20–500 Hz band from a direction inside the (circularized) neutrino error region.
    Required for interpreting both the null and the exclusion distances as constraints on joint sources (§2, §4).
  • domain assumption X-Pipeline coherent excess-power ranking plus consistency tests is a valid detection statistic for generic short GW transients.
    Method foundation citing Sutton et al. 2010; Was et al. 2012 (§2).
  • domain assumption Bronze/Gold astrophysical-signal probabilities and IceCat-1 reconstructions are reliable external inputs.
    Event selection and sky grids taken from IceCat-1 / GCN (§4).
  • ad hoc to paper Circularizing HEALPix 90% contours to a single Gaussian σ is adequate for the coherent sky grid.
    Explicit approximation in §4 because the pipeline defaults to Gaussian error regions; checked only against GCN angular errors.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Search for gravitational waves associated with high-energy neutrinos detected by IceCube during the third observing run of LIGO-Virgo." pith.science (2026). https://pith.science/paper/2WXV3ZJJ

@misc{pith2026260728539,
  author       = {Pith},
  title        = {Pith review of: Search for gravitational waves associated with high-energy neutrinos detected by IceCube during the third observing run of LIGO-Virgo},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2WXV3ZJJ}},
  note         = {Machine review of arXiv:2607.28539}
}
read the original abstract

We search for generic gravitational-wave transients associated with high-energy neutrinos detected by the IceCube Neutrino Observatory during the third Observing Run (O3) of Advanced LIGO, Advanced Virgo, and KAGRA. We perform an unmodeled, targeted search, which is sensitive to gravitational-wave signals weaker than those reported in real time or in the gravitional-wave transient catalog, and can thus uncover coincidences missed by typical neutrino follow-up searches. We find no statistically significant gravitational-wave signal and set lower bounds on the distance of possible gravitational-wave sources for different emission models.

Figures

Figures reproduced from arXiv: 2607.28539 by the authors.

Figure 1
Figure 1. Cumulative p-value distribution for the analyzed events (blue). The dashed black line shows the median expectation obtained from p￾values drawn from a uniform distribution, while the gray lines indicate the 2σ-band around this median [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Exclusion distances at 90% confidence level obtained for the different waveform classes for Bronze (top row) and Gold (bottom row) events. from IceCat-1 detected by IceCube, using data collected by LIGO and Virgo during Observing run O3. Typical sources for such coincident emission are short GRBs, expected mostly from the merger of a binary neutron star system, or long GRBs, ex￾pected mostly from the core collapse o… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 2 linked inside Pith

  1. [1]

    G., Ackermann, M., Adams, J., et al

    Aartsen, M. G., Ackermann, M., Adams, J., et al. 2017, JINST, 12, P03012, [Erratum: JINST 19, E05001 (2024)]

    2017

  2. [2]

    G., Abbott, R., Abouelfettouh, I., et al

    Abac, A. G., Abbott, R., Abouelfettouh, I., et al. 2024, Astrophys. J., 977, 255

    2024

  3. [3]

    G., Abe, A., Abouelfettouh, I., et al

    Abac, A. G., Abe, A., Abouelfettouh, I., et al. 2026 [arXiv:2605.27223]

    arXiv 2026

  4. [4]

    2023a, Astrophys

    Abbasi, R., Ackermann, M., Adams, J., et al. 2023a, Astrophys. J., 959, 96, [Er- ratum: Astrophys.J. 971, 192 (2024)]

    2024

  5. [5]

    2024, Astrophys

    Abbasi, R., Ackermann, M., Adams, J., et al. 2024, Astrophys. J., 964, 126

    2024

  6. [6]

    Abbasi, R. et al. 2025, PoS, ICRC2025, 1224

    2025

  7. [7]

    2016, Phys

    Abbott, B., Abbott, R., Abbott, T., et al. 2016, Phys. Rev. Lett., 116, 061102

    2016

  8. [8]

    D., Abraham, S., et al

    Abbott, R., Abbott, T. D., Abraham, S., et al. 2021, Astrophys. J., 915, 86

    2021

Show all 32 references
  1. [9]

    D., Acernese, F., et al

    Abbott, R., Abbott, T. D., Acernese, F., et al. 2022, Astrophys. J., 928, 186

    2022

  2. [10]

    2024, Astrophys

    Abbott, R., Abe, H., Acernese, F., et al. 2024, Astrophys. J., 966, 137

    2024

  3. [11]

    2021, Astrophys

    Abe, K., Bronner, C., Hayato, Y ., et al. 2021, Astrophys. J., 918, 78

    2021

  4. [12]

    2013, JCAP, 06, 008

    Adrian-Martinez, S., Al Samarai, I., Albert, A., et al. 2013, JCAP, 06, 008

    2013

  5. [13]

    2024, JCAP, 04, 026

    Aiello, S., Albert, A., Alves Garre, S., et al. 2024, JCAP, 04, 026

    2024

  6. [14]

    2023, JCAP, 04, 004

    Albert, A., Alves, S., André, M., et al. 2023, JCAP, 04, 004

    2023

  7. [15]

    2017, Astrophys

    Albert, A., André, M., Anghinolfi, M., et al. 2017, Astrophys. J. Lett., 850, L35

    2017

  8. [16]

    2013, Rev

    Ando, S., Baret, B., Bartos, I., et al. 2013, Rev. Mod. Phys., 85, 1401

    2013

  9. [17]

    2011, Astropart

    Baret, B., Bartos, I., Bouhou, B., et al. 2011, Astropart. Phys., 35, 1

    2011

  10. [18]

    Fynbo, J. P. U., Watson, D., Thöne, C. C., et al. 2006, Nature, 444, 1047

    2006

  11. [19]

    P., Barthelmy, S

    Gehrels, N., Norris, J. P., Barthelmy, S. D., et al. 2006, Nature, 444, 1044 Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, Astrophys. J., 622, 759

    2006

  12. [20]

    2020, Mon

    Guetta, D., Rahin, R., Bartos, I., & Della Valle, M. 2020, Mon. Not. Roy. Astron. Soc., 492, 843 IceCube Collaboration. 2020, High Energy Neutrino Cascade Alerts, Tech. rep., NASA, accessed: 2025-11-24 IceCube Collaboration. 2020, GRB Coordinates Network, 27235, 1

    2020

  13. [21]

    S., Murase, K., Bartos, I., et al

    Kimura, S. S., Murase, K., Bartos, I., et al. 2018, Phys. Rev. D, 98, 043020

    2018

  14. [22]

    P., et al

    Levan, A., Malesani, D., Gompertz, B. P., et al. 2023, Nature Astron., 7, 976

    2023

  15. [23]

    J., Gompertz, B

    Levan, A. J., Gompertz, B. P., Salafia, O. S., et al. 2024, Nature, 626, 737

    2024

  16. [24]

    & Bartos, I

    Murase, K. & Bartos, I. 2019, Ann. Rev. Nucl. Part. Sci., 69, 477

    2019

  17. [25]

    & Ohme, F

    Pannarale, F. & Ohme, F. 2014, Astrophys. J. Lett., 791, L7

    2014

  18. [26]

    C., Gompertz, B

    Rastinejad, J. C., Gompertz, B. P., Levan, A. J., et al. 2022, Nature, 612, 223

    2022

  19. [27]

    2016, Phys

    Senno, N., Murase, K., & Meszaros, P. 2016, Phys. Rev. D, 93, 083003

    2016

  20. [28]

    Sutton, P. J. 2013 [arXiv:1304.0210]

    2013 arXiv

  21. [29]

    J., Jones, G., Chatterji, S., et al

    Sutton, P. J., Jones, G., Chatterji, S., et al. 2010, New J. Phys., 12, 053034

    2010

  22. [30]

    2018, Astron

    Tanga, M., Krühler, T., Schady, P., et al. 2018, Astron. Astrophys., 615, A136 van Putten, M. H. P. M. 2001, Phys. Rev. Lett., 87, 091101 van Putten, M. H. P. M., Levinson, A., Lee, H. K., et al. 2004, Phys. Rev. D, 69, 044007 von Kienlin, A., Meegan, C. A., Paciesas, W. S., e...

    2018

  23. [31]

    J., Jones, G., & Leonor, I

    Was, M., Sutton, P. J., Jones, G., & Leonor, I. 2012, Phys. Rev. D, 86, 022003

    2012

  24. [32]

    R., Biwer, C., Fairhurst, S., et al

    Williamson, A. R., Biwer, C., Fairhurst, S., et al. 2014, Phys. Rev. D, 90, 122004 Article number, page 6 Vereecken et al.: Search for GWs associated with IceCube neutrinos during O3 Table 3.Neutrino events analyzed in this work, and their 90% exclusion distance for each sourc...

    2014

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.