REVIEW 3 major objections 5 minor 40 references
Pre-supernova neutrino timing can warn of core collapse 14 hours ahead, roughly doubling the lead time of conventional rate-only alarms.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 08:47 UTC pith:EM2C5ZSY
load-bearing objection Useful extension of Sheshukov's pre-SN alarm with realistic detector setups, but the equal-FAR claim is undercut by the uncalibrated multi-template minimum. the 3 major comments →
Development of an early warning method incorporating pre-supernova neutrino light curves
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a log-likelihood-ratio test which incorporates the time evolution of the pre-SN neutrino event rate—rather than only the total count—can detect an imminent core collapse earlier without sacrificing false-alarm control. The test computes the likelihood that the observed event times come from a background-plus-signal process using each theoretical light curve as a template, and treats the core-collapse time t* as a nuisance parameter, reporting the maximum likelihood ratio over t*. Calibrating the global false-alarm rate with toy Monte Carlo under background-only, the method reaches a false-alarm rate of one per century at 14.5 hours before collapse for KamLAND alone
What carries the argument
The central object is a log-likelihood-ratio (LLR) test statistic for a sliding time window, built from an inhomogeneous Poisson process whose time-varying rate is the sum of background and a predicted signal light curve R_s(t - t*). The unknown core-collapse time t* is profiled out by maximizing the LLR; the query compares the resulting statistic against thresholds calibrated by Monte Carlo to yield a global false-alarm rate. Multiple theoretical light-curve models (Odrzywolek, Yoshida, Kato, Patton) serve as templates, with the reported sensitivity taken as the minimum false-alarm rate across the model set, and the mass-ordering dependence enters through the MSW survival probability.
Load-bearing premise
The improvement relies on the real star's neutrino light curve being one of the library of templates used in the likelihood; if the true pre-supernova emission has a different time shape, the reported 14-hour warning is not guaranteed.
What would settle it
Inject a simulated pre-supernova signal from a model not in the reference library (for example, a 12-solar-mass star, or a model with different nuclear reaction rates) into the KamLAND/SK background streams, and measure the alarm lead time at a 1-per-century false-alarm rate; if it falls well below 14 hours, the library-robustness claim fails. Similarly, run the alarm for the equivalent of a century of background-only data and count alarms—more than one would falsify the false-alarm calibration.
If this is right
- The rate+time alarm issues alerts 14.0–14.7 h before core collapse for a 15 M_sun Betelgeuse-like star at 150 pc, versus 8.2–12.3 h for rate-only, at the same false-alarm rate.
- Each detector individually with rate+time matches the conventional combined rate-only alarm, so the system remains effective if one detector is down.
- The method extends sensitivity to more distant or weaker pre-SN signals because it uses timing information beyond total counts.
- It maintains a global false-alarm rate of one per century, calibrated via toy Monte Carlo including trials factors.
- The combined rate+time configuration yields the earliest alerts among all tested configurations.
Where Pith is reading between the lines
- The profiled-likelihood template approach could transfer to other detectors (such as JUNO or Hyper-Kamiokande) and to other transient signals with predictable time profiles, such as gravitational-wave chirps or gamma-ray bursts.
- Because the likelihood uses the characteristic oxygen- and silicon-burning features, a triggered alert also encodes the burning stage, giving astrophysical information about the star's final hours—not just a yes/no alarm.
- The method's lead-time gain depends on the library covering the true stellar models; an observed event could be used to validate and update the library, turning the alarm system into a measurement tool.
- The robustness of the quoted false-alarm rate to the multi-template trials factor is a testable extension: a long background-only run should produce no more than one alarm per century.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a pre-supernova neutrino early-warning method that uses the time profile of the expected event rate, not just the total rate. The test statistic is a profile log-likelihood ratio over the unknown core-collapse time, applied to sliding windows in KamLAND, SK-Gd, and their combination. The authors simulate signals from existing pre-SN models, calibrate the false-alarm rate with background-only toy Monte Carlo, and report warning lead times at a 1-per-century FAR threshold. The headline result is that the rate+time analysis outperforms the conventional rate-only analysis at the same false-alarm rate, with combined KamLAND+SK giving the earliest alerts (Table 2).
Significance. If the equal-FAR comparison is properly established, this is a valuable and practical improvement for real-time supernova early warning, directly relevant to SNEWS 2.0. The paper has clear strengths: the inhomogeneous-Poisson likelihood in Eqs. (2)-(5) is standard and correctly formulated; the global FAR is addressed with toy MC rather than naive single-window p-values; the detector backgrounds are realistic; and the performance is summarized in a concrete table with multiple models and mass orderings. The main risk is not the likelihood construction but the handling of the multi-template model library, which the authors themselves flag by reporting the minimum FAR over the ensemble. A second, independent issue is an inconsistency in the combined-analysis claim. These are fixable with a clear calibration statement and a corrected comparison, but they are load-bearing for the paper's central claim.
major comments (3)
- [Section 6, final paragraph; Table 2] The paper states that for the rate+time analysis it evaluates an ensemble of reference light-curve models and 'report[s] the minimum FAR across the model set.' If the alarm is triggered whenever any of the library templates crosses the threshold, the operational false-alarm rate is not the minimum per-template FAR but the union over all templates, which is larger. The text mentions a trials factor for the sliding window but does not state that the toy-MC calibration is performed on an ensemble statistic (e.g., the maximum of the per-template LLR or the minimum p-value). Without such a calibration, the claimed 'same false alarm rate' comparison between rate-only and rate+time in Table 2 is not established. Please either specify that the global FAR is calibrated on the ensemble-max statistic, or recompute the lead times with the correct trials-factor correction.
- [Section 7 and Table 2] There is an internal inconsistency in the combined-analysis claim. For the Odrzywolek 15 M_sun NO case, Table 2 lists KamLAND rate+time as 14.5 h and Combined rate+time as 14.0 h. However, the text in Section 7 states that the combined analysis 'delivers earlier alerts than either detector alone,' and the Discussion says the combined rate+time 'yields the earliest alerts among all configurations.' A combined analysis that includes KamLAND data cannot be worse than KamLAND alone if the test statistic is properly constructed; the discrepancy suggests either a typo in Table 2 or a difference in the treatment of the two analyses. This must be corrected and explained, since the 'combined is best' result is one of the paper's headline conclusions.
- [Section 7, first paragraph; Section 9] The paper promises robustness to model mismatch: 'To assess robustness to model mismatch, we evaluate sensitivity when the injected true pre-SN light curve model differs from the reference model.' No such results are shown in Section 7 or anywhere else; the only injected signals are from Odrzywolek and Patton, both of which are also included in the reference template library (Table 1). Thus the reported lead times are matched-filter results, and the claimed robustness from using 'multiple time profiles' is unsupported. Either add the mismatch-study results or explicitly limit the claim to the matched-template case.
minor comments (5)
- [Eq. (2)] The likelihood factorization is correct, but it would help to specify the domain of t* (e.g., t* > t0) and to note that R_B(t) is treated as constant over the window in this implementation.
- [Abstract] The abstract says 'we propose an alarm method,' but Section 1 attributes the rate+time likelihood approach to Sheshukov et al. (2021). The novel contribution here is the treatment of the core-collapse time as a profiled nuisance parameter and the realistic multi-detector evaluation; wording such as 'we develop and evaluate' would be more accurate.
- [Table 2] The KamLAND rate-only IO entry is listed as 'N/A' while the text does not explain whether no alarm is issued before collapse or whether the number is too small to report. A brief note would avoid confusion.
- [Section 3.1] The selection criteria are given in detail, but the text says 'sub-MeV energy threshold' while the prompt energy cut is 0.9 MeV; this is consistent but could be clarified as the effective analysis threshold after selection.
- [Section 4] The background rates are quoted in the text (0.19 day^-1 for KamLAND and 12.4 day^-1 for SK-Gd) but not listed in a table; a small table of background components would improve reproducibility.
Circularity Check
No circular derivation: the rate+time LLR is built from an inhomogeneous Poisson model and calibrated under H0; the main caveat (multi-template minimum-FAR selection) is a calibration/soundness issue, not a circular one.
full rationale
The central test statistic is constructed in Eqs. (2)-(5) from first principles: the background-only likelihood in Eq. (3) and signal+background likelihood in Eq. (2) are products of Poisson count terms and time PDFs, and the test statistic is the profile likelihood ratio over the core-collapse time t*. Nothing in this construction is defined in terms of the warning times in Table 2. The false-alarm rate is then obtained by toy Monte Carlo under H0, with a trials factor for the sliding window. Injecting Odrzywolek or Patton models and testing with the same library is a matched-filter sensitivity study: the injected signal is known and the templates are fixed external stellar-evolution models; this does not fit any parameter to the quantity being reported. The one caveat is that the paper states "In this study, we evaluate an ensemble of models and report the minimum F AR across the model set" without explicitly stating that the multi-template selection is included in the H0 calibration. If not included, the effective global FAR would exceed 1/century, making the quoted lead times optimistic. But that is a statistical-calibration concern, not a reduction of the output to the input by construction. Self-citations to previous KamLAND/SK papers (e.g., S. Abe et al. 2024) provide detector configurations, background rates, and the existence of a combined alarm; these are empirical inputs, not an unverified uniqueness theorem, and they are not load-bearing for the LLR derivation.
Axiom & Free-Parameter Ledger
free parameters (3)
- Sliding analysis window T =
KamLAND: 200 h; SK: 48 h
- Reference signal template =
15 Msun, 150 pc, normal mass ordering (default); 15-25 Msun mass range
- Template-selection rule over the model ensemble =
minimum FAR across the reference library
axioms (6)
- domain assumption Pre-SN neutrino light-curve models (Odrzywolek, Yoshida, Kato, Patton) accurately represent the true neutrino emission of a collapsing 15-25 Msun star.
- domain assumption MSW flavor conversion with survival probabilities p = 0.680 (NO) and p = 0.023 (IO).
- domain assumption Background rates are known and constant at 0.19 day^-1 (KamLAND) and 12.4 day^-1 (SK) per the 2024 medium reactor-activity scenario.
- standard math Events are independent draws from an inhomogeneous Poisson process with known expected rate.
- standard math KamLAND and SK data streams are statistically independent.
- domain assumption Detector response, cuts, and efficiencies follow the published KamLAND/SK-Gd analyses (IBD selection, BDT, mini-balloon cut).
read the original abstract
Massive stars ($M>8\mathrm{M_\odot}$) emit neutrinos known as pre-supernova (pre-SN) neutrinos through thermal and nuclear interactions for cooling the stellar core during the final stage of stellar evolution. Real-time monitoring of their pre-SN neutrino interaction rate offers a crucial opportunity to issue an early warning to a core-collapse supernova. Some neutrino detectors, including KamLAND and Super-Kamiokande already operate pre-SN alarm systems based on a statistically significant excess of the observed event rate over the expected background. To improve alarm sensitivity, an alarm method which incorporates the time evolution of the observed pre-SN neutrino event rate was proposed in A. Sheshukov et al. (2021). We evaluate the performance of the light-curve likelihood approach under realistic KamLAND and SK operating conditions, including realistic background rates, global false-alarm-rate calibration, and the combined alarm. The results demonstrate a significant improvement in the alarm time and distance compared to the conventional rate-only method, while maintaining the same false alarm rate.
Figures
Reference graph
Works this paper leans on
-
[1]
2018, arXiv e-prints, arXiv:1805.04163, doi: 10.48550/arXiv.1805.04163
Abe, K., Abe, K., Aihara, H., et al. 2018, arXiv e-prints, arXiv:1805.04163, doi: 10.48550/arXiv.1805.04163
-
[2]
Abe, K., et al. 2022, Nucl. Instrum. Meth. A, 1027, 166248, doi: 10.1016/j.nima.2021.166248
arXiv 2022
-
[3]
Abe, K., Bronner, C., Hayato, Y., et al. 2024, Nuclear Instruments and Methods in Physics Research A, 1065, 169480, doi: 10.1016/j.nima.2024.169480
arXiv 2024
-
[4]
Abe, S., et al. 2008, Phys. Rev. Lett., 100, 221803, doi: 10.1103/PhysRevLett.100.221803
-
[5]
Abe, S., et al. 2024, Astrophys. J., 973, 140, doi: 10.3847/1538-4357/ad5fee
-
[6]
Abe, S., Araki, T., Chiba, K., et al. 2025, Phys. Rev. Lett., 135, 262501, doi: 10.1103/jkf6-48j8
-
[7]
Abi, B., et al. 2021, Eur. Phys. J. C, 81, 423, doi: 10.1140/epjc/s10052-021-09166-w
-
[8]
2024, JCAP, 2024, 057, doi: 10.1088/1475-7516/2024/01/057 Al Kharusi, S., BenZvi, S
Abusleme, A., Adam, T., Ahmad, S., et al. 2024, JCAP, 2024, 057, doi: 10.1088/1475-7516/2024/01/057 Al Kharusi, S., BenZvi, S. Y., Bobowski, J. S., et al. 2021, New Journal of Physics, 23, 031201, doi: 10.1088/1367-2630/abde33
-
[9]
Albanese, V., Alves, R., Anderson, M. R., et al. 2021, arXiv e-prints, arXiv:2104.11687, doi: 10.48550/arXiv.2104.11687
-
[10]
Krivosheina, I. V. 1987, JETP Lett., 45, 589
1987
-
[11]
Antonioli, P., et al. 2004, New J. Phys., 6, 114, doi: 10.1088/1367-2630/6/1/114
-
[12]
Asakura, K., et al. 2016, Astrophys. J., 818, 91, doi: 10.3847/0004-637X/818/1/91
-
[13]
Bionta, R. M., et al. 1987, Phys. Rev. Lett., 58, 1494, doi: 10.1103/PhysRevLett.58.1494
-
[14]
Dolan, M. M., Mathews, G. J., Lam, D. D., et al. 2016, ApJ, 819, 7, doi: 10.3847/0004-637X/819/1/7
-
[15]
2007, Earth and Planetary Science Letters, 258, 147, doi: https://doi.org/10.1016/j.epsl.2007.03.038
Enomoto, S., Ohtani, E., Inoue, K., & Suzuki, A. 2007, Earth and Planetary Science Letters, 258, 147, doi: https://doi.org/10.1016/j.epsl.2007.03.038
-
[16]
Fukuda, Y., et al. 2003, Nucl. Instrum. Meth. A, 501, 418, doi: 10.1016/S0168-9002(03)00425-X
-
[17]
Gando, A., et al. 2012, Phys. Rev. C, 85, 045504, doi: 10.1103/PhysRevC.85.045504
-
[18]
Gando, A., et al. 2013, Phys. Rev. D, 88, 033001, doi: 10.1103/PhysRevD.88.033001
-
[19]
Hirata, K., et al. 1987, Phys. Rev. Lett., 58, 1490, doi: 10.1103/PhysRevLett.58.1490
-
[20]
Huber, P. 2011, Phys. Rev. C, 84, 024617, doi: 10.1103/PhysRevC.84.024617
-
[21]
2020, The Astrophysical Journal, 902, 63, doi: 10.3847/1538-4357/abb8db
Joyce, M., Leung, S.-C., Moln´ ar, L., et al. 2020, The Astrophysical Journal, 902, 63, doi: 10.3847/1538-4357/abb8db
-
[22]
2015, ApJ, 808, 168, doi: 10.1088/0004-637X/808/2/168
Kato, C., Delfan Azari, M., Yamada, S., et al. 2015, ApJ, 808, 168, doi: 10.1088/0004-637X/808/2/168
-
[23]
Kato, C., Ishidoshiro, K., & Yoshida, T. 2020, Annual Review of Nuclear and Particle Science, 70, 121, doi: 10.1146/annurev-nucl-040620-021320
-
[24]
Kato, C., Nagakura, H., Furusawa, S., et al. 2017, Astrophys. J., 848, 48, doi: 10.3847/1538-4357/aa8b72
-
[25]
Machado, L. N., et al. 2022, Astrophys. J., 935, 40, doi: 10.3847/1538-4357/ac7f9c
-
[26]
A., Lhuillier, D., Fallot, M., et al
Mueller, T. A., Lhuillier, D., Fallot, M., et al. 2011, Phys. Rev. C, 83, 054615, doi: 10.1103/PhysRevC.83.054615
-
[27]
2010, Acta Phys
Odrzywolek, A., & Heger, A. 2010, Acta Phys. Polon. B, 41, 1611
2010
-
[28]
Odrzywolek, A., Misiaszek, M., & Kutschera, M. 2004, Astropart. Phys., 21, 303, doi: 10.1016/j.astropartphys.2004.02.002
-
[29]
Pagliaroli, G., Vissani, F., Costantini, M. L., & Ianni, A. 2009, Astroparticle Physics, 31, 163, doi: 10.1016/j.astropartphys.2008.12.010
-
[30]
Patton, K. M., Lunardini, C., Farmer, R. J., & Timmes, F. X. 2017, Astrophys. J., 851, 6, doi: 10.3847/1538-4357/aa95c4
-
[31]
Paxton, B., Bildsten, L., Dotter, A., et al. 2011, Astrophys. J. Suppl., 192, 3, doi: 10.1088/0067-0049/192/1/3 Pietrzy´ nski, G., Graczyk, D., Gallenne, A., et al. 2019, Nature, 567, 200, doi: 10.1038/s41586-019-0999-4
-
[32]
Raj, N., Takhistov, V., & Witte, S. J. 2020, Phys. Rev. D, 101, 043008, doi: 10.1103/PhysRevD.101.043008
-
[33]
2012, Annual Review of Nuclear and Particle Science, 62, 81, doi: 10.1146/annurev-nucl-102711-095006
Scholberg, K. 2012, Annual Review of Nuclear and Particle Science, 62, 81, doi: 10.1146/annurev-nucl-102711-095006
-
[34]
2021, JCAP, 2021, 053, doi: 10.1088/1475-7516/2021/12/053
Sheshukov, A., Vishneva, A., & Habig, A. 2021, JCAP, 2021, 053, doi: 10.1088/1475-7516/2021/12/053
-
[35]
Suzuki, A. 2014, Eur. Phys. J. C, 74, 3094, doi: 10.1140/epjc/s10052-014-3094-x
-
[36]
2013, ApJ, 771, 28, doi: 10.1088/0004-637X/771/1/28
Takahashi, K., Yoshida, T., & Umeda, H. 2013, ApJ, 771, 28, doi: 10.1088/0004-637X/771/1/28
-
[37]
2016, MNRAS, 456, 1320, doi: 10.1093/mnras/stv2649
Yamada, S. 2016, MNRAS, 456, 1320, doi: 10.1093/mnras/stv2649
-
[38]
Vogel, P., Schenter, G. K., Mann, F. M., & Schenter, R. E. 1981, Phys. Rev. C, 24, 1543, doi: 10.1103/PhysRevC.24.1543
-
[39]
Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015, doi: 10.1103/RevModPhys.74.1015
-
[40]
2016, PhRvD, 93, 123012, doi: 10.1103/PhysRevD.93.123012
Yoshida, T., Takahashi, K., Umeda, H., & Ishidoshiro, K. 2016, PhRvD, 93, 123012, doi: 10.1103/PhysRevD.93.123012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.