Pith. sign in

REVIEW 2 major objections 5 minor 27 references

Eleven Year Search for Supernovae with the IceCube Neutrino Observatory

T0 review · 2 major / 5 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 An eleven-year IceCube search will set a neutrino-based limit that covers 99 percent of Galactic core-collapse supernovae at or above a conservative low-mass flux.

desk verdict A methods-readiness conference paper for a still-blinded 11-year IceCube supernova search; the projected 99% coverage is the one substantive claim, and it hinges entirely on an ASTERIA/GEANT-4 equivalence that is asserted, not shown. read the letter →

arxiv 1908.07249 v1 pith:RFDBLDKL submitted 2019-08-20 astro-ph.HE

Robert Cross , Alexander Fritz , Spencer Griswold (for the IceCube Collaboration) This is my paper · ORCID
classification astro-ph.HE
keywords core-collapsesupernovaesupernovaneutrinosIceCubeGalacticratefaileddust-obscuredneutrinoburstdetectionSNDAQ
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

IceCube was built for high-energy neutrinos, but its unusually quiet optical sensors also let it see the tens-of-MeV neutrino burst of a Galactic core-collapse supernova as a collective rise in the count rate of all its photomultipliers. This paper presents the method and simulation tools for an eleven-year search over 3,911 days of data (3,670 days after quality cuts) aimed at supernovae that are hidden by dust or that fail to explode, since neutrinos escape where light does not. Its central claim is that once the data are unblinded, and if no burst is seen, IceCube will set an upper limit valid for 99 percent of all Galactic core-collapse supernovae whose neutrino flux is at least as high as the conservative 8.8-solar-mass benchmark model. A cut on the search statistic at 9.3 sigma, with detector systematics included, is what achieves that 99-percent retention. If correct, this would be the strongest neutrino-based constraint yet on the rate of obscured and failed stellar collapses in the Milky Way.

What carries the argument

The argument runs on the moving-average test statistic $\xi=\Delta\mu/\sigma_{\Delta\mu}$, computed in overlapping 1.5-second windows stepped by 500 milliseconds, with the largest value in each 10-second interval kept. Here $\Delta\mu$ is the most likely collective deviation of all optical-module hit rates from their running average, and $\sigma_{\Delta\mu}$ is estimated from the data themselves, so non-Poissonian dark-noise fluctuations are folded into the significance. The expected signal shapes come from ASTERIA, a fast parameterized simulation of the detector response that the paper states, without showing verification, reproduces the full simulation; simulated bursts are added to recorded rates and then passed through the same trigger logic that runs online. The retention claim follows from the separation between the simulated $\xi$ distributions for the five benchmark models and the background-only distribution in the data.

What would settle it

Re-run the fast simulation and the full detector simulation on identical 8.8-solar-mass supernovae at 10 kiloparsecs and compare the per-module and total hit rates; if the disagreement is larger than the 14-percent systematic band used in the paper, the $\xi>9.3\sigma$ cut would not in fact retain 99 percent of supernovae. The same check can be done on the unblinded data: the measured background $\xi$ distribution should match the simulated background, with no excess tail above 9 $\sigma$ once muon-correlated hits are accounted for.

Watch

Extended reading notes

Core claim

The central result is a projected sensitivity statement rather than a detection. For the five benchmark core-collapse models considered — a low-mass 8.8-solar-mass electron-capture collapse, an 11.2-solar-mass star, a 27-solar-mass star, a forced explosion of a 30-solar-mass star, and a failed collapse that forms a black hole — the expected values of the test statistic $\xi$ separate cleanly from the background distribution measured in real IceCube data. Only the lowest-mass model shows any overlap. The black-hole model, which produces a high average neutrino energy and a sharp cutoff after about one second, yields $\xi$ values that lie entirely beyond the background. Because the simulated supernova positions are drawn from published radial distributions of Galactic structure and the signals are added to recorded detector rates, the projected $\xi>9.3\sigma$ cut retains 99 percent of the simulated population under the conservative flux assumption. The search deliberately requires no optical counterpart, so it is sensitive to dust-obscured and failed supernovae that optical surveys would miss.

Load-bearing premise

The 99-percent coverage and the $\xi>9.3\sigma$ cut assume the fast parameterized simulation produces the same detector response as the full simulation, and that adding simulated bursts to recorded rates faithfully mimics real backgrounds; the paper asserts both without showing verification.

Editorial extensions

If this is right

  • If no burst is found after unblinding, the result will set an upper limit on the rate of Galactic core-collapse supernovae that is valid for 99 percent of all collapses with neutrino fluxes at or above the conservative low-mass model, independent of any optical detection.
  • The same test statistic already runs online in real time, so a genuine Galactic supernova would be flagged immediately and the buffered full waveforms would allow follow-up timing measurements with other neutrino detectors, potentially triangulating the source.
  • The clean separation of the black-hole model's $\xi$ distribution from background means that a failed collapse in the Milky Way should be detected with high significance if it happens during the live period; a null result therefore constrains the rate of such collapses.
  • The search is independent of external information, so it covers dust-obscured supernovae and complements optical and gravitational-wave based rate estimates, including the 1.7-to-2.5 per century expectation from stellar evolution.

Reading between the lines

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

  • An implicit consequence the paper does not spell out: a null result would bound the rate of failed collapses specifically, because the black-hole model produces the strongest signal of all five benchmarks; separating that rate from the total would require a model-dependent fit.
  • The quoted 9.3-sigma cut depends on the assumptions of a normal neutrino mass hierarchy and a particular Galactic progenitor distribution; if either changes, the 99-percent coverage cut shifts, so comparisons with other experiments should be made model by model rather than on a single number.
  • If the equivalence between the fast simulation and the full simulation were publicly verified, the same tool could be extended to estimate IceCube's sensitivity to other short MeV-scale neutrino transients, such as pre-supernova neutrinos, where the collective-rate method would also apply.
  • A reader could test the pipeline on the Magellanic Clouds by applying the same simulation with muon-subtracted rates, which the paper identifies as future work; this would tell whether the 99-percent coverage extends beyond the Milky Way.
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.

Referee Report

2 major / 5 minor

Summary. This manuscript, an IceCube Collaboration proceedings paper submitted to ICRC 2019, describes the methodology for an eleven-year search for Galactic core-collapse supernovae using 3911 days of IceCube data from April 2008 to December 2018. The paper introduces the detector response simulation chain (GEANT-4, ASTERIA, SNOwGLoBES), defines the SNDAQ test statistic ξ, and presents a projected sensitivity study using five progenitor models with distances drawn from Galactic radial distributions. The central forward claim is that, if no signal is observed upon unblinding, the analysis will provide an upper limit valid for 99% of all Galactic core-collapse supernovae with neutrino fluxes at or above the conservative 8.8 solar-mass Hüdepohl model, corresponding to a significance cut at ξ > 9.3σ with systematics included (Section 3, Figures 5 and 6).

Significance. If the projected sensitivity is realized, the analysis would deliver the most stringent neutrino-based constraint on the rate of obscured and failed core-collapse supernovae in the Milky Way, improving on the Baksan and LVD limits. The paper has strengths: it uses a large live-time dataset, a clear definition of the test statistic on real background data, multiple progenitor models spanning an order of magnitude in predicted flux, and an explicit, falsifiable projection of 99% signal retention. The use of recorded DOM rates as background and the inclusion of detector systematic uncertainties in the projected cut are positive features. However, the projected 99% coverage and the ξ > 9.3σ cut rest entirely on an asserted equivalence between the fast parameterized simulation ASTERIA and the full GEANT-4 Monte Carlo; this equivalence is not demonstrated in the manuscript, and the paper itself notes order-of-magnitude uncertainties in supernova flux models. The central claim is therefore plausible but not yet fully substantiated.

major comments (2)
  1. [Section 3, Figures 5-6] The projected 99% retention and the ξ > 9.3σ cut depend on the sentence 'The parameterized simulation [10], verified to produce the same results as the GEANT-4 Monte Carlo.' This verification is load-bearing for every sensitivity number in the paper: the ξ distributions for all five progenitor models are generated with ASTERIA, and the 99% coverage threshold is read directly from these distributions. Yet the manuscript provides no comparison plots, residuals, validation statistics, or public code/repository to support the equivalence. ASTERIA is described as a faster and less sophisticated parameterized simulation that uses GEANT-4 effective volumes, so a small bias in effective volume or energy response at ~10 MeV could shift the 99% retention threshold and change the projected upper limit. Please provide a quantitative ASTERIA-vs-GEANT-4 comparison for the full analysis chain that produces ξ, including the 1.5 s maximization, the 10 s selection, and the 250 μs deadtime handling, and state the energy range over which the verification holds. Without this, the central claim is an unsupported calibration step.
  2. [Section 3, simulated signal injection] The paper states that 'The simulated hit rates were added to uniformly sampled DOM rates recorded by SNDAQ between 2008 and 2018.' This injection step is not self-evidently valid with respect to the artificial deadtime. The deadtime fraction in IceCube depends on the total PMT rate, so adding simulated signal hits to recorded background rates changes the total rate and hence the deadtime loss. The text does not state whether the ASTERIA hit rates are corrected for deadtime as a function of the combined signal-plus-background rate, or whether the recorded DOM rates already include the deadtime suppression in a way that makes simple addition correct. If the deadtime coupling is not modeled, the effective signal size is biased and the 99% retention cut derived from Figure 6 will be incorrect. Please clarify the deadtime treatment in the injection procedure and justify the additive model quantitatively, or modify the simulation to account for the rate-dependent deadtime.
minor comments (5)
  1. [Section 3, benchmark models] The text says a 30 M⊙ progenitor model 'yields 1.97×10^53 M⊙'; the unit should presumably be erg, not solar masses.
  2. [Section 2] There is a typo 'artificial deadime' where 'deadtime' is intended.
  3. [Figure 6, left panel] The black hole model is described as lying entirely beyond the right edge of the plot, but the caption does not explicitly state this; please add a note in the caption so that readers do not infer the model is absent from the analysis.
  4. [References] Reference [10] is given only as 'ASTERIA: A Supernova TEst Routine for IceCube Analysis, 2019' with no author list, publication venue, or DOI. Since the manuscript relies on ASTERIA for the central projection, a full citation or a publicly accessible code repository is needed.
  5. [Section 1, deadtime formula] The formula for the deadtime suppression factor '0.87/(1+Rdark(t)/NDOM·τ)' is stated without derivation or a citation; please clarify the origin of the 0.87 factor and define all quantities precisely, as this factor directly affects the simulated rates.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the projected sensitivity is derived from external supernova models, a data-driven test statistic, and simulated signal injections, with no fitted parameter recycled into the claim.

full rationale

The paper's central forward claim is a projected upper limit for Galactic core-collapse supernovae, obtained by simulating 18,000 supernovae from five external benchmark models (Hüdepohl, Mirizzi, Nakazato, Sumiyoshi) with distances drawn from published Galactic distributions, adding the simulated hit rates to recorded SNDAQ rates, and then computing the test statistic xi from the data. The 99% retention cut at xi > 9.3 sigma is read off the simulated cumulative distribution, not fitted to any observed signal. The only arguably load-bearing internal element is the statement that ASTERIA 'was verified to produce the same results as the GEANT-4 Monte Carlo,' but this is a simulation-accuracy claim, not a definitional equivalence between the paper's inputs and its output; even if the verification were incomplete, that would be a correctness or validation risk rather than circularity. The benchmark flux models, progenitor distributions, and comparison limits (Baksan, LVD) are all external to this paper. No parameter is fitted to a subset of data and then renamed a prediction, and no uniqueness theorem or prior self-citation is invoked to force the analysis choice. The paper is explicitly a methods and sensitivity presentation for a blinded search, so the absence of a measured limit is consistent with the stated plan rather than a hidden reuse of inputs.

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

The central sensitivity projection adds no new physics entities and no data-fitted constants, but it relies on a simulation-calibrated significance cut, a Gaussian background assumption, and an asserted equivalence between ASTERIA and GEANT-4. These are the loaded assumptions the reader pays for upstream.

free parameters (1)
  • SN retention significance cut = xi > 9.3 sigma (after systematics)
    Chosen so that 99% of simulated Galactic supernovae are retained in the mock-injection study; it is a simulation-calibrated threshold that defines the projected limit and is not an independent prediction.
assumptions (3)
  • domain assumption The test statistic xi is distributed as a zero-mean unit Gaussian when no correlations are present in the DOM rates.
    Section 3 states this and uses it to interpret xi in units of Gaussian sigma. The paper itself shows muon bundles broaden the distribution (Figure 6), so the assumption is only approximate and material for threshold setting.
  • domain assumption ASTERIA reproduces the GEANT-4 detector response for supernova neutrino signals.
    Section 3 asserts this verification without presenting plots, code, or a public repository; all projected sensitivities in Figures 5 and 6 depend on the equivalence.
  • domain assumption A normal neutrino mass hierarchy and the Ahlers et al. progenitor distribution give the conservative sensitivity estimate.
    Section 3 says 'To err on the side of caution, we will assume the neutrino mass hierarchy to be normal and apply the progenitor distribution of [16]'. The final upper limit's meaning depends on this choice.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Eleven Year Search for Supernovae with the IceCube Neutrino Observatory." pith.science (2026). https://pith.science/paper/RFDBLDKL

@misc{pith2026190807249,
  author       = {Pith},
  title        = {Pith review of: Eleven Year Search for Supernovae with the IceCube Neutrino Observatory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RFDBLDKL}},
  note         = {Machine review of arXiv:1908.07249}
}
abstract

The IceCube Neutrino Observatory, which instruments 1$\,$km$^3$ of clear ice at the geographic South Pole, was mainly designed to detect particles with energies in the multi-GeV to PeV range. Due to ice temperatures between $-20^\circ$C to $-43^\circ$C and the low radioactivity of the ice, the dark noise rates of the 5160 photomultiplier tubes forming the IceCube lattice are of order 500 Hz, which is particularly low for 10 inch photomultipliers. Therefore, IceCube can extend its searches to bursts of $\mathcal{O}$(10$\,$MeV) neutrinos lasting several seconds, which are expected to be produced by Galactic core collapse supernovae. By observing a uniform rise in all photomultiplier rates, IceCube can provide a particularly high statistical precision for the neutrino rate from supernovae in the inner part of our Galaxy ($<$ 20 kpc). In this paper, the tools and the method to study potential obscured or failed core collapse supernovae in our Galaxy are presented. The analysis will be based on 3911 days of IceCube data taken between April 17, 2008 and December 31, 2018.

Figures

Figures reproduced from arXiv: 1908.07249 by the authors.

Figure 1
Figure 1. Top and side view of ∼ 3.4×105 simulated supernova ν interaction vertices registered by IceCube DOMs. The dust layer between -1950m and -2050m and the denser DeepCore subarray are clearly visible. Construction of IceCube finished in 2011, and since 2015 the trigger-capable uptime of the detector has averaged 99.7% around the clock. Due to the non-Poissonian character of the dark noise in the IceCube DOMs [4], the da… view at source ↗
Figure 2
Figure 2. Simulated DOM hits in IceCube for a supernova from a 13 M star located 10 kpc from Earth [12], assuming several different neutrino oscillation scenarios [10]. instrumented IceCube and DeepCore volumes yields at least one registered Cherenkov photon. The corresponding positions of neutrino interactions are shown in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Comparison of effective total volumes of ideal and background free detectors for 8.8 M [13], 20 M [14], and 40 M [15] stellar mod￾els. The solid curves show the IceCube effective volume when the artificial deadtime of 250 µs is applied; the dashed curves show the results with￾out the deadtime. For comparison, the estimated fiducial volumes of Super-Kamiokande, JUNO and Hyper-Kamiokande are shown as dashed orange lin… view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Left: Changes in the rate of simulated DOM hits from a 13 M supernova [12] 10 kpc from Earth found by Bayesian Blocks [6]. Right: RMS error in the estimate of the start time t0 vs. supernova distance [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: Test statistic ξ (in units of Gaussian σ) vs. progenitor dis￾tance, simulated with ASTERIA and SNDAQ for the five models discussed in this paper: the O-Ne-Mg core from Hüdepohl et al. [13]; an 11.2 M star [23]; a 27 M star [23]; a forced explosion of a 30 M star [12]; …
Figure 6
Figure 6. Figure 6: Left: The distribution of the test statistic ξ is shown for data taken in 2014 with and without muon subtraction as well as for four supernova models. A slight overlap occurs only for the model with the lowest progenitor mass [13]. Right: Fraction of supernovae missed …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 25 canonical work pages

  1. [10]

    ASTERIA: A Supernova TEst Routine for IceCube Analysis , 2019

    work page 2019

  2. [1]

    Achterberg et al., Astropart

    IceCube Collaboration, A. Achterberg et al., Astropart. Phys. 26 (2006) 155–173

    work page 2006

  3. [2]

    Halzen, J

    F. Halzen, J. E. Jacobsen, and E. Zas, Phys. Rev. D53 (1996) 7359–7361

    work page 1996

  4. [3]

    A. W. Alsabti and P. Murdin, eds., Handbook of Supernovae. Springer, New York, 2017

    work page 2017

  5. [4]

    Abbasi et al., Astron

    IceCube Collaboration, R. Abbasi et al., Astron. Astrophys. 535 (2011) A109

    work page 2011

  6. [5]

    M. T. Keil. PhD thesis, Munich, Max Planck Inst., 2003. astro-ph/0308228

    work page Pith review arXiv 2003

  7. [6]

    IceCube Collaboration, PoS(ICRC2017)936 (2018)

    work page 2018

  8. [7]

    IceCube Collaboration, V . Baum, D. Heereman, and R. Bruijn,Braz. J. Phys. 44 (2014) 0444

    work page 2014

Show all 27 references
  1. [8]

    Antonioli et al., New J

    P. Antonioli et al., New J. Phys. 6 (2004) 114

    2004

  2. [9]

    IceCube Collaboration, PoS(ICRC2019)865 (these proceedings)

  3. [11]

    IceCube Collaboration, PoS(ICRC2019)975 (these proceedings)

  4. [12]

    Nakazato et al., Astrophys

    K. Nakazato et al., Astrophys. J. Suppl. 205 (2013) 2

    2013

  5. [13]

    Hüdepohl et al., Phys

    L. Hüdepohl et al., Phys. Rev. Lett. 104 (2010) 251101

    2010

  6. [14]

    Totani et al., Astrophys

    T. Totani et al., Astrophys. J. 496 (1998) 216–225

    1998

  7. [15]

    Sumiyoshi, S

    K. Sumiyoshi, S. Yamada, and H. Suzuki, Astrophys. J. 667 (2007) 382–394

    2007

  8. [16]

    Ahlers, P

    M. Ahlers, P. Mertsch, and S. Sarkar, Phys. Rev. D80 (2009) 123017

    2009

  9. [17]

    MÃijhlbeier, H

    T. MÃijhlbeier, H. Nunokawa, and R. Zukanovich Funchal, Phys. Rev. D88 (2013) 085010

    2013

  10. [18]

    Brdar, M

    V . Brdar, M. Lindner, and X.-J. Xu,JCAP 1804 (2018) 025

    2018

  11. [19]

    Giunti and C

    C. Giunti and C. W. Kim, Fundamentals of Neutrino Physics and Astrophysics . Oxford Univ. Press, Oxford, UK, 2007

    2007

  12. [20]

    Petkov et al.,PoS(ICRC2017)960 (2018)

    Baksan Collaboration, V . Petkov et al.,PoS(ICRC2017)960 (2018)

    2018

  13. [21]

    LVD Collaboration, C. F. Vigorito et al., PoS(ICRC2017)1017 (2018)

    2018

  14. [22]

    S. J. Smartt, Ann. Rev. Astron. Astrophys. 47 (2009) 63–106

    2009

  15. [23]

    Mirizzi et al., Riv

    A. Mirizzi et al., Riv. Nuovo Cim. 39 (2016) 1–112

    2016

  16. [24]

    Mattila et al., Astrophys

    S. Mattila et al., Astrophys. J. 756 (2012) 111

    2012

  17. [25]

    Mirizzi, G

    A. Mirizzi, G. G. Raffelt, and P. D. Serpico, JCAP 0605 (2006) 012

    2006

  18. [26]

    IceCube Collaboration, PoS(ICRC2019)855 (these proceedings)

  19. [27]

    IceCube Collaboration, PoS(ICRC2017)1052 (2018). 8

    2018

Pith tools

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