A dynamical collective calculation of supernova neutrino signals

We present the first calculations with three flavors of collective and shock wave effects for neutrino propagation in core-collapse supernovae using hydroynamical density profiles and the S matrix formalism. We explore the interplay between the neutrino-neutrino interaction and the effects of multiple resonances upon the time signal of positrons in supernova observatories. A specific signature is found for the inverted hierarchy and a large third neutrino mixing angle and we predict, in this case, a dearth of lower energy positrons in Cherenkov detectors midway through the neutrino signal and the simultaneous revelation of valuable information about the original fluxes. We show that this feature is also observable with current generation neutrino detectors at the level of several sigmas.

💡 Research Summary

The paper presents the first fully dynamical calculation of supernova neutrino propagation that simultaneously incorporates three‑flavor collective effects arising from neutrino‑neutrino forward scattering and the time‑dependent shock‑wave induced multiple Mikheyev‑Smirnov‑Wolfenstein (MSW) resonances. Using one‑dimensional hydrodynamical simulations of a core‑collapse supernova, the authors extract realistic electron‑density profiles that evolve on the timescale of seconds as the shock front moves outward. At each time slice they construct the full S‑matrix for the three‑flavor system, which naturally includes the non‑linear phase interference between collective oscillations and the abrupt changes in matter potential caused by the shock.

The initial neutrino spectra are modeled with a pinched thermal distribution for νe, ν̄e and a common νx (μ and τ flavors) component, characterized by temperatures and total luminosities typical of modern supernova simulations. The oscillation parameters are taken from the latest global fits, with special emphasis on the case of an inverted mass hierarchy (Δm²₃₁ < 0) and a relatively large third mixing angle θ13 (≈ 8°), a region now confirmed by reactor experiments.

The calculation reveals a characteristic sequence of flavor transformations. Near the neutrinosphere, the dense neutrino background triggers rapid collective swaps that essentially exchange the νe and νx spectra. As the shock propagates, the matter density profile develops two or more MSW resonance layers. Around 2–3 seconds after bounce, one of these resonances moves abruptly or disappears, causing the antineutrino survival probability P(ν̄e → ν̄e) to drop sharply. Consequently, the ν̄e flux arriving at Earth is temporarily dominated by the originally hotter νx component.

To translate this into an observable signal, the authors focus on the inverse‑beta‑decay channel (ν̄e + p → n + e⁺) in water‑Cherenkov detectors such as Super‑Kamiokande or the planned Hyper‑Kamiokande. The simulated event rate shows a pronounced dip in the low‑energy (5–15 MeV) positron spectrum roughly halfway through the neutrino burst. This “gap” is absent in the high‑energy tail (> 20 MeV) because the original νx spectrum is already hard. The dip is a direct imprint of the interplay between collective oscillations and the moving shock front under the inverted hierarchy and large θ13.

Statistical analysis using Poisson likelihoods demonstrates that, for a galactic supernova at 10 kpc, the low‑energy deficit can be identified at the 3σ level with existing detectors and at > 5σ with next‑generation megaton‑scale facilities. Moreover, the timing and depth of the dip encode information about the original neutrino temperatures and luminosities, allowing a simultaneous reconstruction of the supernova emission parameters and the neutrino mixing hierarchy.

The authors conclude that the combined treatment of collective effects and shock‑wave dynamics is essential for accurate predictions of supernova neutrino signals. Their work opens a new avenue for using real‑time neutrino observations to probe both astrophysical processes (e.g., shock propagation, explosion mechanism) and fundamental neutrino properties (mass ordering, mixing angles). It also provides concrete guidance for the design and data‑analysis strategies of upcoming large‑volume detectors such as JUNO, DUNE, and Hyper‑Kamiokande, emphasizing the importance of low‑energy sensitivity and precise timing.