pith. sign in

arxiv: 2605.06401 · v1 · submitted 2026-05-07 · ✦ hep-ph · astro-ph.HE· hep-ex

The temporal distribution of SN1987A neutrino events

Riccardo Maria Bozza , Vigilante di Risi , Veronica Oliviero , Giulia Ricciardi , Francesco Vissani This is my paper

Pith reviewed 2026-05-08 08:21 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.HEhep-ex
keywords SN1987Aneutrino eventstiming alignmentKamiokande-IIBaksanIMBinverse beta decayaccretion phase
0
0 comments X

The pith

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.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

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.

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.

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.

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.

Where Pith is reading between the lines

These are 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.

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.

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.

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

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

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.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 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.
axioms (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.

pith-pipeline@v0.9.0 · 5527 in / 1447 out tokens · 31371 ms · 2026-05-08T08:21:09.156854+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

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

  9. [9]

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

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

    work page 1987

  10. [10]

    Costantini, A

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

    work page 2004

  11. [11]

    LoSecco,The Case for Neutrinos From Sn1987a,Phys

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

    work page 1989

  12. [12]

    Vissani,Comparative analysis of SN1987A antineutrino fluence,J

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

    work page 2015

  13. [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

    work page 1998

  14. [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)

    work page 2025

  15. [15]

    Vissani and A

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

    work page 2021

  16. [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 –

    work page 2002

  17. [17]

    Pagliaroli, F

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

    work page 2009

  18. [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

    work page 2009

  19. [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

    work page 2023

  20. [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)

    work page 1970

  21. [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

    work page 2006

  22. [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)

    work page 1990

  23. [23]

    Kuzminov,The Baksan Neutrino Observatory,Eur

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

    work page 2012

  24. [24]

    Lunardini and A

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

    work page 2001

  25. [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

    work page 2024

  26. [26]

    Yaron and A

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

    work page 2012

  27. [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

    work page 2009

  28. [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

    work page 2003

  29. [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

    work page 2016

  30. [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

    work page 2000

  31. [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]

    work page arXiv 2024

  32. [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 –

    work page 2026