pith. sign in

arxiv: 2606.06409 · v1 · pith:55MP5MHPnew · submitted 2026-06-04 · 🌌 astro-ph.HE

A high-energy neutrino flare associated with nearby bright interacting supernova SN 2021foa

Ming-Xuan Lu , Yun-Feng Liang , Xue-Rui Ouyang , Da-Bin Lin , Xiang-Gao Wang , Yi-Zhong Fan This is my paper

Pith reviewed 2026-06-28 00:21 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords neutrino astronomysupernovaeIceCubehigh-energy neutrinosinteracting supernovaecentral engineSN 2021foa
0
0 comments X

The pith

IceCube data show a 4-sigma neutrino excess from the interacting supernova SN 2021foa coinciding with its optical peak.

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

The paper presents evidence that high-energy neutrinos were emitted by the nearby supernova SN 2021foa. A time-dependent search in IceCube muon-track data found a cluster of events 16 to 22 days after discovery, matching the timing of the optical maximum. Monte Carlo tests put the probability of a random coincidence at about 6.7 times 10 to the minus 5, or roughly 4 sigma. SN 2021foa is notable for its flip-flop spectral changes between Type IIn and Ibn, indicating complex mass ejection. The neutrino luminosity greatly exceeds the supernova's optical and kinetic output, pointing to a central engine powering a jet that stalls in surrounding gas.

Core claim

Using IceCube's second data release, a time-dependent analysis identifies a neutrino clustering with maximum test statistic of approximately 28.2 that occurs 16-22 days after the discovery of SN 2021foa. Simulations show this excess has a chance probability of 6.7×10^{-5}, corresponding to 4.0 sigma significance. The spatial and temporal match indicates the neutrinos come from the supernova. Given that the inferred neutrino energy far surpasses the optical radiative energy and ejecta kinetic energy, the emission is attributed to a delayed central engine driving a choked jet within the dense circumstellar medium. SN 2021foa exhibits repeated transitions between hydrogen-rich and helium-rich s

What carries the argument

The time-dependent maximum-likelihood analysis for neutrino clustering in direction and arrival time, combined with Monte Carlo background simulations to assess significance.

If this is right

  • The neutrino output being orders of magnitude larger than optical and kinetic energies implies an additional energy source beyond standard supernova processes.
  • A choked jet from the central engine can explain the neutrino production in the dense circumstellar medium.
  • The flip-flop behavior of SN 2021foa highlights the role of episodic mass loss in enabling such neutrino emission.
  • Similar associations may exist for other nearby interacting supernovae.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the association holds, it opens the possibility that choked jets in supernovae are common neutrino sources at high energies.
  • Monitoring other bright interacting supernovae with neutrino telescopes could test the prevalence of such flares.
  • The complex mass-loss history inferred from spectra may be key to understanding why only certain supernovae produce detectable neutrinos.
  • This event could guide searches for multi-messenger signals in future supernova observations.

Load-bearing premise

The specific 16-22 day time window and the neutrino clustering criteria were fixed before looking at the data rather than chosen after seeing a signal.

What would settle it

Repeating the analysis with a larger set of background trials that include all possible time windows around the supernova would show whether the significance drops below 3 sigma when accounting for the look-elsewhere effect.

read the original abstract

While core-collapse supernovae have been widely discussed as potential neutrino sources, definitive observational evidence has remained elusive. In this work, we report evidence of high-energy neutrino emission in the direction of supernova SN 2021foa, which is one of the closest and brightest interacting supernovae observed to date. Using the second data release of muon track data from the IceCube Neutrino Observatory, we conducted a time-dependent analysis and identified a neutrino clustering that temporally coincides with the optical peak of SN 2021foa, occurring approximately 16 to 22 days after the discovery date, with a maximum test statistic of $\sim 28.2$. Monte Carlo simulations indicate that the probability of observing such a neutrino excess by chance is $\sim6.7 \times 10^{-5}$, corresponding to a significance of $\sim4.0\,\sigma$. The spatial and temporal correlation strongly suggests that the neutrinos originate from the supernova. SN 2021foa is a unique "flip-flop" supernova; its spectra repeatedly transitioned between hydrogen-rich (Type IIn) and helium-rich (Type Ibn) phases within 50 days post-peak, reflecting a violent and complex mass-loss history of its progenitor. The inferred neutrino energy exceeds the optical radiative energy and ejecta kinetic energy of the supernova by orders of magnitude, suggesting that the neutrino emission is likely powered by a delayed central engine driving a jet that is choked within the dense circumstellar medium.

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 / 0 minor

Summary. The manuscript reports evidence for high-energy neutrino emission from the nearby interacting supernova SN 2021foa using IceCube muon-track data. A time-dependent analysis identifies a neutrino clustering temporally coincident with the optical peak (16–22 days post-discovery) with maximum test statistic ~28.2; Monte Carlo simulations give a chance probability of ~6.7×10^{-5} (~4σ), which the authors interpret as indicating a supernova origin. The paper further notes the supernova’s unusual “flip-flop” spectral evolution and argues that the inferred neutrino energy implies a delayed central engine driving a choked jet.

Significance. If the reported association survives a complete trials correction and background-model validation, the result would constitute the first statistically compelling link between high-energy neutrinos and a core-collapse supernova, directly constraining neutrino-production mechanisms in dense circumstellar environments and supporting the existence of central engines in interacting supernovae.

major comments (2)
  1. [Abstract] Abstract (Monte Carlo significance paragraph): The claimed ~4.0σ significance rests on a Monte Carlo p-value of ~6.7×10^{-5} for a neutrino excess in the specific 6-day window 16–22 days post-discovery. The manuscript must explicitly state whether this window, the clustering metric, and the test-statistic definition were fixed a priori; if any optimization occurred after inspecting the data, the Monte Carlo ensemble must repeat the identical search procedure on background-only realizations to furnish the correct frequentist probability. Without this information the reported significance cannot be verified.
  2. [Time-dependent analysis] Time-dependent analysis section: The abstract supplies no description of the background model (atmospheric neutrinos, detector systematics, or data-selection cuts) or the number of trials implicit in the time-dependent search. These details are load-bearing for the central statistical claim and must be provided with sufficient quantitative information for an independent assessment of the p-value.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and for identifying key points that require clarification to strengthen the statistical claims in our manuscript. We address each major comment below and will revise the paper accordingly to ensure the significance and analysis details are fully transparent and verifiable.

read point-by-point responses

  1. Referee: [Abstract] Abstract (Monte Carlo significance paragraph): The claimed ~4.0σ significance rests on a Monte Carlo p-value of ~6.7×10^{-5} for a neutrino excess in the specific 6-day window 16–22 days post-discovery. The manuscript must explicitly state whether this window, the clustering metric, and the test-statistic definition were fixed a priori; if any optimization occurred after inspecting the data, the Monte Carlo ensemble must repeat the identical search procedure on background-only realizations to furnish the correct frequentist probability. Without this information the reported significance cannot be verified.

    Authors: We agree that explicit confirmation of the a priori status of the search parameters is essential for verifying the frequentist significance. The 6-day window (16–22 days post-discovery) was selected to align with the independently observed optical peak of SN 2021foa, providing a physically motivated choice rather than a post-hoc optimization on the neutrino data alone. The clustering metric and test-statistic definition were fixed before performing the analysis. In the revised manuscript, we will add a clear statement in the abstract and methods confirming these parameters were predetermined. The Monte Carlo ensemble already implements the identical full search procedure (including time-window scanning) on background-only realizations, properly incorporating any trials factors; we will emphasize this explicitly to allow independent verification of the ~4σ claim. revision: yes

  2. Referee: [Time-dependent analysis] Time-dependent analysis section: The abstract supplies no description of the background model (atmospheric neutrinos, detector systematics, or data-selection cuts) or the number of trials implicit in the time-dependent search. These details are load-bearing for the central statistical claim and must be provided with sufficient quantitative information for an independent assessment of the p-value.

    Authors: We acknowledge that the current description of the background model and trials factors in the time-dependent analysis section is insufficient for full independent assessment. The analysis uses the publicly released IceCube muon-track dataset with standard atmospheric neutrino background modeling, data-selection cuts, and detector response as documented in the second data release. To address this, the revised manuscript will expand the time-dependent analysis section with quantitative details, including the expected atmospheric neutrino rate in the region of interest, considerations of relevant systematics, and the effective number of trials arising from the time-dependent search procedure. This will enable direct evaluation of the reported p-value. revision: yes

Circularity Check

0 steps flagged

No circularity: observational excess from external IceCube data with MC p-value

full rationale

The paper reports a statistical neutrino excess (TS ~28.2, p ~6.7e-5 via MC) temporally aligned with SN 2021foa optical peak. This is an observational claim using external public IceCube muon-track data and independent optical light-curve information. No equation or result reduces by construction to a parameter fitted from the same neutrino sample; the MC ensemble is presented as the direct frequentist measure under the null. No self-citation chain, ansatz smuggling, or self-definitional loop is present in the derivation. The analysis is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard statistical methods for time-dependent point-source searches in neutrino data and the assumption that the chosen time window and event selection are unbiased. No free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Standard IceCube background estimation and Monte Carlo simulation procedures are valid and correctly implemented.
    The significance calculation depends on these procedures being free of unaccounted systematics.

pith-pipeline@v0.9.1-grok · 5817 in / 1263 out tokens · 23794 ms · 2026-06-28T00:21:12.469513+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

40 extracted references · 4 canonical work pages · 1 internal anchor

  1. [1]

    G.et al.Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector.Science342, 1242856 (2013)

    Aartsen, M. G.et al.Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector.Science342, 1242856 (2013)

    2013

  2. [2]

    G.et al.Observation of High-Energy Astrophysical Neutrinos in Three Years of IceCube Data.Phys

    Aartsen, M. G.et al.Observation of High-Energy Astrophysical Neutrinos in Three Years of IceCube Data.Phys. Rev. Lett.113, 101101 (2014)

    2014

  3. [3]

    G.et al.Multimessenger observations of a flaring blazar coincident with high- energy neutrino IceCube-170922A.Science361, eaat1378 (2018)

    Aartsen, M. G.et al.Multimessenger observations of a flaring blazar coincident with high- energy neutrino IceCube-170922A.Science361, eaat1378 (2018)

    2018

  4. [4]

    G.et al.Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert.Science361, 147–151 (2018)

    Aartsen, M. G.et al.Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert.Science361, 147–151 (2018)

    2018

  5. [5]

    Science378, 538–543 (2022)

    Abbasi, R.et al.Evidence for neutrino emission from the nearby active galaxy NGC 1068. Science378, 538–543 (2022)

    2022

  6. [6]

    Abbasi, R.et al.Observation of high-energy neutrinos from the Galactic plane.Science 380, adc9818 (2023)

    2023

  7. [7]

    Stein, R.et al.A tidal disruption event coincident with a high-energy neutrino.Nature Astronomy5, 510-518 (2021)

    2021

  8. [8]

    Reusch, S.et al.Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino.Phys. Rev. Lett.128, 221101 (2022)

    2022

  9. [9]

    van Velzen, S.et al.Establishing accretion flares from supermassive black holes as a source of high-energy neutrinos.Mon. Not. R. Astron. Soc.529, 2559-2576 (2024)

    2024

  10. [10]

    R. L. Li, C. Yuan, H. N. He, Y . Wang, B. Y . Zhu, Y . F. Liang, N. Jiang and D. M. Wei, Mon. Not. Roy. Astron. Soc.546, no.4, stag188 (2026)

    2026

  11. [11]

    & Ahlers, M

    Murase, K., Guetta, D. & Ahlers, M. Hidden Cosmic-Ray Accelerators as an Origin of TeV- PeV Cosmic Neutrinos.Phys. Rev. Lett.116, 071101 (2016)

    2016

  12. [12]

    & Waxman, E

    Murase, K. & Waxman, E. Constraining High-Energy Cosmic Neutrino Sources: Implications and Prospects.Phys. Rev. D94, 103006 (2016)

    2016

  13. [13]

    Sridhar, N., Metzger, B. B. & Fang, K. High-energy Neutrinos from Gamma-Ray-faint Accretion-powered Hypernebulae.Astrophys. J.960, 74 (2024)

    2024

  14. [14]

    A., Lacki, B

    Murase, K., Thompson, T. A., Lacki, B. C. & Beacom, J. F. New Class of High-Energy Transients from Crashes of Supernova Ejecta with Massive Circumstellar Material Shells. Phys. Rev. D84, 043003 (2011). 15

    2011

  15. [15]

    New Prospects for Detecting High-Energy Neutrinos from Nearby Supernovae

    Murase, K. New Prospects for Detecting High-Energy Neutrinos from Nearby Supernovae. Phys. Rev. D97, 081301 (2018)

    2018

  16. [16]

    & Meszaros, P

    Senno, N., Murase, K. & Meszaros, P. Choked Jets and Low-Luminosity Gamma-Ray Bursts as Hidden Neutrino Sources.Phys. Rev. D93, 083003 (2016)

    2016

  17. [17]

    Zhu, J.-P.et al.High-energy Neutrinos from Choked Gamma-Ray Bursts in Active Galactic Nucleus Accretion Disks.Astrophys. J. Lett.911, L19 (2021)

    2021

  18. [18]

    & Waxman, E

    Meszaros, P. & Waxman, E. TeV neutrinos from successful and choked gamma-ray bursts. Phys. Rev. Lett.87, 171102 (2001)

    2001

  19. [19]

    & Waxman, E

    Razzaque, S., Meszaros, P. & Waxman, E. TeV neutrinos from core collapse supernovae and hypernovae.Phys. Rev. Lett.93, 181101 (2004)

    2004

  20. [20]

    & Beacom, J

    Ando, S.’i. & Beacom, J. F. Revealing the supernova-gamma-ray burst connection with TeV neutrinos.Phys. Rev. Lett.95, 061103 (2005)

    2005

  21. [21]

    & Ioka, K

    Murase, K. & Ioka, K. TeV–PeV Neutrinos from Low-Power Gamma-Ray Burst Jets inside Stars.Phys. Rev. Lett.111, 121102 (2013)

    2013

  22. [22]

    J.961, 84 (2024)

    Abbasi, R.et al.Search for Galactic Core-collapse Supernovae in a Decade of Data Taken with the IceCube Neutrino Observatory.Astrophys. J.961, 84 (2024)

    2024

  23. [23]

    Abbasi, R.et al.Constraining High-energy Neutrino Emission from Supernovae with Ice- Cube.Astrophys. J. Lett.949, L12 (2023)

    2023

  24. [24]

    & Kamionkowski, M

    Chang, P.-W., Zhou, B., Murase, K. & Kamionkowski, M. High-energy neutrinos from choked-jet supernovae: Searches and implications.Phys. Rev. D109, 103041 (2024)

    2024

  25. [25]

    & Zheng, D

    Ji, S., Wang, Z., Zhu, L. & Zheng, D. High-energy neutrino emission from the Type IIn supernova SN 2017hcd.arXiv e-printsarXiv:2603.03679 (2026)

    work page arXiv 2026

  26. [26]

    G.et al.The Detection of a SN IIn in Optical Follow-up Observations of IceCube Neutrino Events.Astrophys

    Aartsen, M. G.et al.The Detection of a SN IIn in Optical Follow-up Observations of IceCube Neutrino Events.Astrophys. J.811, 52 (2015)

    2015

  27. [27]

    & Zhang, H.-Q

    Lu, M.-X., Liang, Y .-F., Wang, X.-G. & Zhang, H.-Q. SN 2023syz and SN 2025cbj: Two Type IIn Supernovae Associated with IceCube High-energy Neutrinos.Astrophys. J.997, 145 (2026)

    2026

  28. [28]

    Stein, R.et al.SN 2023uqf: An Interacting Supernova Coincident with a High-Energy Neutrino.Preprint at https://arxiv.org/abs/2508.08355(2025). 16

    work page arXiv 2025

  29. [29]

    Schlegel, E. M. A new subclass of type II supernovae ?.Mon. Not. R. Astron. Soc.244, 269-271 (1990)

    1990

  30. [30]

    Filippenko, A. V . Optical Spectra of Supernovae.Annu. Rev. Astron. Astrophys.35, 309-355 (1997)

    1997

  31. [31]

    Farias, D.et al.SN 2021foa: The ”Flip-Flop” Type IIn / Ibn supernova.arXiv e-prints arXiv:2409.01359 (2024)

    work page arXiv 2024

  32. [32]

    Gangopadhyay, A.et al.SN 2021foa: the bridge between SN IIn and Ibn.Mon. Not. Roy. Astron. Soc.537, 2898–2917 (2025)

    2025

  33. [33]

    Astrophys.662, L10 (2022)

    Reguitti, A.et al.SN 2021foa, a transitional event between a Type IIn (SN 2009ip-like) and a Type Ibn supernova.Astron. Astrophys.662, L10 (2022)

    2022

  34. [34]

    Abbasi, R.et al.IceCube Second Track Data Release IceTracks-DR2: Data from 2008-2022 for Neutrino Source Searches.arXiv e-printsarXiv:2605.19040 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  35. [35]

    Phys.29, 299–305 (2008)

    Braun, J.et al.Methods for point source analysis in high energy neutrino telescopes.As- tropart. Phys.29, 299–305 (2008)

    2008

  36. [36]

    Braun, J.et al.Time-dependent point source search methods in high energy neutrino astron- omy.Astroparticle Physics33, 175-181 (2010)

    2010

  37. [37]

    Hidden Hearts of Neutrino Active Galaxies.Astrophys

    Murase, K. Hidden Hearts of Neutrino Active Galaxies.Astrophys. J. Lett.941, L17 (2022)

    2022

  38. [38]

    Murase, K., Kimura, S. S. & Meszaros, P. Hidden Cores of Active Galactic Nuclei as the Origin of Medium-Energy Neutrinos: Critical Tests with the MeV Gamma-Ray Connection. Phys. Rev. Lett.125, 011101 (2020)

    2020

  39. [39]

    A.et al.The Zwicky Transient Facility Bright Transient Survey

    Perley, D. A.et al.The Zwicky Transient Facility Bright Transient Survey. II. A Public Statistical Sample for Exploring Supernova Demographics.Astrophys. J.904, 35 (2020)

    2020

  40. [40]

    G.et al.The IceCube Neutrino Observatory: Instrumentation and Online Sys- tems.JINST12, P03012 (2017)

    Aartsen, M. G.et al.The IceCube Neutrino Observatory: Instrumentation and Online Sys- tems.JINST12, P03012 (2017). 17 AcknowledgmentsWe acknowledge the IceCube Collaboration for the release of 14-year muon-track data and for providing theSkyLLHsoftware package for neutrino point-source analysis. This work is supported by the National Key Research and Deve...

    2017