Pith. sign in

REVIEW 2 major objections 2 minor 32 references

The temporal distribution of SN1987A neutrino events

T0 review · 2 major / 2 minor · reviewed 2026-05-08 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read By aligning SN1987A neutrino detector timestamps using relative time offsets, the analysis determines clock shifts and shows the first Kamiokande-II event likely originated in the accretion phase rather than the neutronization burst.

desk verdict The paper introduces a new relative-time-offset alignment that sharpens the SN1987A chronology to sub-second level, but the evidence favoring an accretion origin for the first Kamiokande event remains moderate because detector response differences are not fully addressed. read the letter →

arxiv 2605.06401 v1 submitted 2026-05-07 hep-ph astro-ph.HEhep-ex

Riccardo Maria Bozza , Vigilante di Risi , Veronica Oliviero , Giulia Ricciardi , Francesco Vissani This is my paper · ORCID
classification hep-phastro-ph.HEhep-ex
keywords SN1987AneutrinoeventstimingalignmentKamiokande-IIBaksanIMBinversebetadecayaccretionphase
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 aligns the absolute times of neutrino events from the Kamiokande-II, Baksan, and IMB detectors for supernova 1987A by using the relative time offsets between events. This reduces the timing uncertainty from minutes to sub-second precision. The resulting unified timeline indicates that Baksan's clocks were advanced by 30.4 seconds and Kamiokande-II's were delayed by 6.4 seconds. With this alignment, the first Kamiokande-II event is tested for origin and found more consistent with the accretion phase of electron antineutrinos than with the neutronization burst, by a factor of 3 to 6 in likelihood. This supports the conclusion that the detectors observed only inverse beta decay events.

What carries the argument

The relative time offset (RTO) method for aligning absolute timestamps across detectors to achieve sub-second precision in the unified timeline.

What would settle it

Re-analysis of the raw detector data with alternative methods or a new supernova event with known absolute timing that yields different clock shifts would falsify the claimed 30.4 s and 6.4 s corrections.

Watch

Extended reading notes

Core claim

The chi-square analysis of relative time offsets shows that Baksan's absolute timestamps require an advancement of 30.4 s, while those of Kamiokande-II require a delay of about 6.4 s. This provides a unified timeline that favors an accretion phase electron anti-neutrino origin for the first Kamiokande-II event over a neutronization-burst origin, with a likelihood ratio of 3-6 depending on the MSW oscillation scenario, corroborating that only inverse beta decay events were detected.

Load-bearing premise

That the relative time offsets between events recorded at different detectors can be used to determine absolute clock shifts without introducing unaccounted systematic uncertainties.

Editorial extensions

If this is right

  • The analysis provides the most stringent constraints to date on the SN1987A chronology.
  • It establishes a precision benchmark for future Galactic supernova observations.
  • The result supports the standard interpretation that only inverse beta decay events were detected in the SN1987A data.
  • The likelihood ratio of 3-6 favors the accretion phase origin depending on the MSW oscillation scenario.

Reading between the lines

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

  • This clock alignment technique could resolve similar timing uncertainties in data from future multi-detector supernova neutrino observations.
  • The findings increase the reliability of using SN1987A events to test models of core-collapse supernova neutrino emission phases.
  • Further studies could check if the angular distribution inconsistencies are fully resolved by these timing corrections.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 2 minor

Summary. The paper claims to perform a high-precision alignment of the SN1987A neutrino event timestamps from Kamiokande-II and Baksan to the IMB clock by minimizing a chi-square statistic based on relative time offsets (RTO) for the first time, yielding an advancement of 30.4 s for Baksan and a delay of 6.4 s for Kamiokande-II at sub-second precision. Using this unified timeline, it tests the origin of the first Kamiokande-II event and reports a likelihood ratio of 3-6 favoring an accretion-phase electron anti-neutrino (inverse beta decay) origin over a neutronization-burst origin, depending on the MSW scenario; this is presented as corroborating that only IBD events were detected and as providing the most stringent constraints on the SN1987A chronology.

Significance. If the alignment procedure is robust, the work provides a valuable quantitative framework for reconciling the SN1987A datasets, reducing absolute timing uncertainties by two orders of magnitude and supplying a reproducible basis to address the long-standing angular-distribution tension. The explicit chi-square minimization on RTO and the resulting concrete shifts constitute a clear methodological advance that could serve as a benchmark for future Galactic supernova neutrino observations.

major comments (2)
  1. [§3.2] §3.2 (chi-square minimization of RTO): The procedure extracts absolute shifts by treating the observed event times in each detector as samples from the same underlying temporal distribution up to a constant offset. However, the detectors have distinct energy thresholds and efficiencies (KII ~7.5 MeV, Baksan ~10 MeV, IMB ~20 MeV), and the neutrino spectrum hardens during the accretion phase; without explicit correction or marginalization over detector response functions, the fitted shifts of +30.4 s and -6.4 s can be biased at the sub-second level claimed.
  2. [§4.3] §4.3 (likelihood-ratio test): The reported LR of 3-6 for accretion-phase IBD versus neutronization-burst origin for the first KII event is computed after applying the RTO-derived timing correction. Any systematic offset bias identified in the chi-square step propagates directly into this ratio and into the claim that timing alone resolves the angular-distribution tension in favor of pure IBD events.
minor comments (2)
  1. [Abstract] The abstract states that the LR depends on the 'specific MSW oscillation scenario' but the main text does not list the exact scenarios or oscillation parameters used; this should be tabulated or referenced to an equation in §4.
  2. [Figure 1] Figure 1 or 2 (event time distributions): The error model and exclusion criteria for events entering the chi-square fit are not fully specified in the caption or accompanying text; adding a brief statement on how background and efficiency are handled would improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment in turn below, indicating the revisions we will make where appropriate.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (chi-square minimization of RTO): The procedure extracts absolute shifts by treating the observed event times in each detector as samples from the same underlying temporal distribution up to a constant offset. However, the detectors have distinct energy thresholds and efficiencies (KII ~7.5 MeV, Baksan ~10 MeV, IMB ~20 MeV), and the neutrino spectrum hardens during the accretion phase; without explicit correction or marginalization over detector response functions, the fitted shifts of +30.4 s and -6.4 s can be biased at the sub-second level claimed.

    Authors: We agree that the differing energy thresholds and efficiencies, together with spectral evolution, mean the detected event samples are not drawn from identical distributions. The chi-square procedure therefore relies on an approximation that the observed times share a common temporal profile up to a constant offset. This approximation is data-driven and does not explicitly fold in detector response functions. We will revise §3.2 to state this limitation explicitly and to provide a brief estimate of the resulting bias using a simple model of spectrum hardening; the estimate indicates the bias remains well below the reported sub-second statistical precision. The quoted shifts are therefore retained, but the text will now qualify their interpretation. revision: partial

  2. Referee: [§4.3] §4.3 (likelihood-ratio test): The reported LR of 3-6 for accretion-phase IBD versus neutronization-burst origin for the first KII event is computed after applying the RTO-derived timing correction. Any systematic offset bias identified in the chi-square step propagates directly into this ratio and into the claim that timing alone resolves the angular-distribution tension in favor of pure IBD events.

    Authors: The likelihood ratio is evaluated on the aligned timeline, so any bias in the RTO shifts would affect the numerical value. As noted in our response to the §3.2 comment, we will add an explicit bias estimate showing that the effect is sub-dominant to the statistical uncertainty already reflected in the LR range of 3–6. We will also revise §4.3 to state that the LR is conditional on the alignment procedure and to clarify that the conclusion favoring an accretion-phase IBD origin is robust within the quoted range. The claim that the aligned chronology helps address the angular-distribution tension will be qualified accordingly. revision: partial

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity in the derivation chain

full rationale

The paper's central procedure applies standard chi-square minimization to relative time offsets extracted from independent detector event lists to determine absolute clock shifts, then uses the resulting unified timeline for a separate likelihood-ratio comparison of two physical hypotheses for the first event. No equation reduces the claimed offsets or likelihood ratios to fitted parameters by construction, no self-citation supplies a load-bearing uniqueness theorem or ansatz, and the analysis remains self-contained against external observational data without renaming known results.

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

The central claim rests on two fitted clock offsets obtained from chi-square minimization and on standard domain assumptions about neutrino detection channels and oscillation physics.

free parameters (2)
  • Baksan absolute time shift = 30.4 s
    Obtained by chi-square minimization against the IMB reference clock.
  • Kamiokande-II absolute time shift = 6.4 s
    Obtained by chi-square minimization against the IMB reference clock.
assumptions (2)
  • domain assumption Inverse beta decay is the dominant detection process for the observed events
    Invoked when interpreting all events as electron anti-neutrinos and when stating that only inverse beta decay events were detected.
  • domain assumption MSW oscillation parameters govern the flavor evolution from the supernova to Earth
    Used to compute the likelihood ratio that depends on the specific oscillation scenario.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The temporal distribution of SN1987A neutrino events." pith.science (2026). https://pith.science/paper/2605.06401

@misc{pith2026260506401,
  author       = {Pith},
  title        = {Pith review of: The temporal distribution of SN1987A neutrino events},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2605.06401}},
  note         = {Machine review of arXiv:2605.06401}
}
read the original abstract

The interpretation of the SN1987A neutrino data continues to be hindered by significant absolute timing uncertainties and a well-known tension in the angular distributions. We perform a quantitative, high-precision alignment of the Kamiokande-II and Baksan data with the IMB clock, using for the first time the relative time offset (RTO) and reducing the temporal uncertainty by two orders of magnitude to the sub-second level. The chi-square analysis shows that Baksan's absolute timestamps require an advancement of 30.4 s, while those of Kamiokande-II require a delay of about 6.4 s. The knowledge of the unified timeline provides a necessary and rigorous basis for testing whether the first Kamiokande-II event arose from the neutronization burst - a possibility motivated by its angular properties. Our analysis favors an accretion phase electron anti-neutrino origin over a neutronization-burst origin, with a likelihood ratio of 3-6, depending on the specific MSW oscillation scenario. This result corroborates the standard interpretation - that only inverse beta decay events were detected. Our framework yields the most stringent constraints to date on the SN1987A chronology and establishes a precision benchmark for future Galactic supernova observations.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    Hirata et al.,Observation of a neutrino burst from the supernova SN1987A,Phys

    K. Hirata et al.,Observation of a neutrino burst from the supernova SN1987A,Phys. Rev. Lett. 58(1987) 1490

    work page 1987

  2. [2]

    Hirata et al.,Observation in the Kamiokande-II detector of the neutrino burst from supernova SN1987A,Phys

    K. Hirata et al.,Observation in the Kamiokande-II detector of the neutrino burst from supernova SN1987A,Phys. Rev. D38(1988) 448

    work page 1988

  3. [3]

    Bionta et al.,Observation of a neutrino burst in coincidence with supernova 1987A in the large magellanic cloud,Phys

    R.M. Bionta et al.,Observation of a neutrino burst in coincidence with supernova 1987A in the large magellanic cloud,Phys. Rev. Lett.58(1987) 1494

    work page 1987

  4. [4]

    Bratton et al.,Angular distribution of events from SN1987A,Phys

    C. Bratton et al.,Angular distribution of events from SN1987A,Phys. Rev. D37(1988) 3361

    work page 1988

  5. [5]

    Alexeyev et al.,Detection of the neutrino signal from SN 1987A in the LMC using the INR baksan underground scintillation telescope,Phys

    E. Alexeyev et al.,Detection of the neutrino signal from SN 1987A in the LMC using the INR baksan underground scintillation telescope,Phys. Lett. B205(1988) 209

    work page 1988

  6. [6]

    Fiorillo, M

    D. Fiorillo, M. Heinlein, H.-T. Janka, G. Raffelt, E. Vitagliano and R. Bollig,Supernova simulations confront SN 1987A neutrinos,Phys. Rev. D108(2023) 083040

    work page 2023

  7. [7]

    Krivoruchenko,A statistival analysis of angular distribution of neutrino events observed in Kamiokande-II and IMB detectors from supernova SN1987A,Z

    M.I. Krivoruchenko,A statistival analysis of angular distribution of neutrino events observed in Kamiokande-II and IMB detectors from supernova SN1987A,Z. Phys. C44(1989) 633

    work page 1989

  8. [8]

    Bahcall, T

    J.N. Bahcall, T. Piran, W.H. Press and D.N. Spergel,Neutrino Temperatures and Fluxes From the LMC Supernova,Nature327(1987) 682

    work page 1987

Show all 32 references
  1. [9]

    Krauss,Neutrino Spectroscopy of the Supernova Sn1987a,Nature329(1987) 689

    L.M. Krauss,Neutrino Spectroscopy of the Supernova Sn1987a,Nature329(1987) 689

    1987

  2. [10]

    Costantini, A

    M.L. Costantini, A. Ianni and F. Vissani,SN1987A and the properties of neutrino burst,Phys. Rev. D70(2004) 043006

    2004

  3. [11]

    LoSecco,The Case for Neutrinos From Sn1987a,Phys

    J.M. LoSecco,The Case for Neutrinos From Sn1987a,Phys. Rev. D39(1989) 1013

    1989

  4. [12]

    Vissani,Comparative analysis of SN1987A antineutrino fluence,J

    F. Vissani,Comparative analysis of SN1987A antineutrino fluence,J. Phys. G42(2015) 013001

    2015

  5. [13]

    Malgin,Analysis of integral and averaged characteristics of the IMB and Kamioka signals from SN1987A,Nuovo Cim

    A. Malgin,Analysis of integral and averaged characteristics of the IMB and Kamioka signals from SN1987A,Nuovo Cim. C21(1998) 317

    1998

  6. [14]

    Bozza, V

    R.M. Bozza, V. di Risi, G. Matteucci, V. Oliviero, G. Ricciardi and F. Vissani,The flux of electron antineutrinos from supernova SN1987A data,JCAP05(2025) 027 (2025)

    2025

  7. [15]

    Vissani and A

    F. Vissani and A. Gallo Rosso,On the time distribution of supernova antineutrino flux, Symmetry13(2021)

    2021

  8. [16]

    Loredo and D.Q

    T.J. Loredo and D.Q. Lamb,Bayesian analysis of neutrinos observed from supernova SN-1987A,Phys. Rev. D65(2002) 063002. – 12 –

    2002

  9. [17]

    Pagliaroli, F

    G. Pagliaroli, F. Vissani, M.L. Costantini and A. Ianni,Improved analysis of SN1987A antineutrino events,Astropart. Phys.31(2009) 163

    2009

  10. [18]

    Pagliaroli, F

    G. Pagliaroli, F. Vissani, E. Coccia and W. Fulgione,Neutrinos from Supernovae as a Trigger for Gravitational Wave Search,Phys. Rev. Lett.103(2009) 031102

    2009

  11. [19]

    Neto, M.V

    P.D. Neto, M.V. d. Santos, P.C. de Holanda and E. Kemp,SN1987A neutrino burst: limits on flavor conversion,Eur. Phys. J. C83(2023) 459

    2023

  12. [20]

    Sanduleak,A deep objective-prism survey for Large Magellanic Cloud members, Contribution from CTIO89(1970)

    N. Sanduleak,A deep objective-prism survey for Large Magellanic Cloud members, Contribution from CTIO89(1970)

    1970

  13. [21]

    Hosaka et al.,Solar neutrino measurements in Super-Kamiokande-I,Phys

    J. Hosaka et al.,Solar neutrino measurements in Super-Kamiokande-I,Phys. Rev. D73(2006) 112001

    2006

  14. [22]

    W. Gajewski et al.,A Proposal for a Long Baseline Oscillation Experiment using a High Intensity Neutrino Beam from the Fermilab Main Injector to the IMB water Cherenkov Detector,Fermilab-proposal-0805(1990)

    1990

  15. [23]

    Kuzminov,The Baksan Neutrino Observatory,Eur

    V. Kuzminov,The Baksan Neutrino Observatory,Eur. Phys. J. Plus127(2012) 113

    2012

  16. [24]

    Lunardini and A

    C. Lunardini and A. Smirnov,Supernova neutrinos: Earth matter effects and neutrino mass spectrum,Nucl. Phys. B616(2001) 307

    2001

  17. [25]

    Bouchet et al.,JWST MIRI Imager Observations of Supernova SN 1987A,Astrophys

    P. Bouchet et al.,JWST MIRI Imager Observations of Supernova SN 1987A,Astrophys. J. 965(2024) 51

    2024

  18. [26]

    Yaron and A

    O. Yaron and A. Gal-Yam,WISeREP—An Interactive Supernova Data Repository,PASP124 (2012) 668

    2012

  19. [27]

    Ianni, G

    A. Ianni, G. Pagliaroli, A. Strumia, F.R. Torres, F. Villante and F. Vissani,The Likelihood for supernova neutrino analyses,Phys. Rev. D80(2009) 043007

    2009

  20. [28]

    Keil, G.G

    M.T. Keil, G.G. Raffelt and H.-T. Janka,Monte carlo study of supernova neutrino spectra formation,Astrophysical Journal590(2003) 971

    2003

  21. [29]

    Mirizzi, I

    A. Mirizzi, I. Tamborra, H.-T. Janka, N. Saviano, K. Scholberg, R. Bollig et al.,Supernova Neutrinos: Production, Oscillations and Detection,Riv. Nuovo Cim.39(2016) 1

    2016

  22. [30]

    Dighe and A

    A.S. Dighe and A. Smirnov,Identifying the neutrino mass spectrum from the neutrino burst from a supernova,Phys. Rev. D62(2000) 033007

    2000

  23. [31]

    Volpe,Neutrinos from dense environments: Flavor mechanisms, theoretical approaches, observations, and new directions,Rev

    M.C. Volpe,Neutrinos from dense environments: Flavor mechanisms, theoretical approaches, observations, and new directions,Rev. Mod. Phys.96(2024) 025004 [2301.11814]

    2024

  24. [32]

    Data and Codes Repository

    R.M. Bozza, V. di Risi, V. Oliviero, G. Ricciardi and F. Vissani, “Data and Codes Repository.” https://github.com/VOliviero/SN1987A_Analysis, 2026. – 13 –

    2026

Pith tools

Reviewed May 8, 2026 · model on record in the stance chip above.