Density Fluctuation Effects on Collective Neutrino Oscillations in O-Ne-Mg Core-Collapse Supernovae

We investigate the effect of matter density fluctuations on supernova collective neutrino flavor oscillations. In particular, we use full multi-angle, 3-flavor, self-consistent simulations of the evolution of the neutrino flavor field in the envelope of an O-Ne-Mg core collapse supernova at shock break-out (neutrino neutronization burst) to study the effect of the matter density “bump” left by the He-burning shell. We find a seemingly counterintuitive increase in the overall electron neutrino survival probability created by this matter density feature. We discuss this behavior in terms of the interplay between the matter density profile and neutrino collective effects. While our results give new insights into this interplay, they also suggest an immediate consequence for supernova neutrino burst detection: it will be difficult to use a burst signal to extract information on fossil burning shells or other fluctuations of this scale in the matter density profile. Consistent with previous studies, our results also show that the interplay of neutrino self-coupling and matter fluctuation could cause a significant increase in the electron neutrino survival probability at very low energy

💡 Research Summary

This paper investigates how small-scale matter density fluctuations influence collective neutrino flavor oscillations in the early neutrino burst of an O‑Ne‑Mg core‑collapse supernova. The authors focus on the “bump” in the matter density profile produced by the helium‑burning shell, which lies at a radius of roughly 1080–1100 km and creates a localized increase in electron density. Using fully self‑consistent, multi‑angle, three‑flavor simulations (implemented in two independent codes, BULB and FLAT), they compare two scenarios: the original density profile containing the bump (“Bump”) and a synthetic profile where the bump has been removed (“No Bump”).

The neutrino emission model is a pure νₑ neutronization burst with an average energy of 11 MeV (degeneracy parameter η = 3), a total luminosity of 10⁵³ erg s⁻¹, and a neutrinosphere radius of 60 km. The mixing parameters are taken as Δm²₁₂ = 7.6 × 10⁻⁵ eV², Δm²₃₁ = 2.4 × 10⁻³ eV², θ₁₂ = 0.59, θ₂₃ = π/4, θ₁₃ = 0.1, and δ = 0, with the normal mass hierarchy assumed throughout.

The simulations reveal a counter‑intuitive result: the survival probability of νₑ in the heaviest mass eigenstate (state 3) is larger when the bump is present (P_H ≈ 0.852) than when it is removed (P_H ≈ 0.759). This runs opposite to the naive expectation that a smoother density profile should make the evolution more adiabatic and thus increase survival. The effect is most pronounced at low neutrino energies (≲ 5 MeV), where the bump induces multiple MSW resonances at the atmospheric Δm² scale. Consequently, the low‑energy νₑ spectrum shows an enhanced survival probability for the Bump case, while the swap (spectral split) energy shifts upward in the No Bump case.

To interpret these findings, the authors employ the neutrino flavor isospin (NFIS) formalism. They write the evolution equations for each NFIS vector s_ω, which precesses under the combined influence of the vacuum Hamiltonian H_v, the matter potential H_e, and the self‑interaction term µ(r). By integrating over the energy spectrum they define a collective polarization vector S. In the limit where the spectral distribution g(ω) is sharply peaked (appropriate for the narrow νₑ burst), the dynamics reduce to the precession of S around an effective field that is the sum of H_v weighted by the average oscillation frequency and H_e.

The key insight is that the matter scale height H = |d ln n_e/dr|⁻¹ at the resonance determines whether the evolution follows the “synchronous MSW” regime (where the ensemble behaves like a single neutrino experiencing only the matter potential) or the “regular precession” mode (where collective self‑coupling dominates and all neutrinos rotate together in flavor space). The bump creates rapid variations in H, causing the system to repeatedly cross between these regimes. Each crossing can temporarily “freeze” the collective precession, allowing a larger fraction of neutrinos to remain in the heavy mass state. When the bump is removed, the scale height varies more smoothly, the system spends more time in the regular precession mode, and a larger portion of νₑ transitions to lower mass states, reducing P_H.

The authors also examine the observational implications. Current and planned large‑volume detectors (e.g., Hyper‑Kamiokande, DUNE, JUNO) will be capable of measuring the νₑ burst, but extracting subtle features such as the 10 % level differences in survival probability caused by the bump is challenging. Statistical uncertainties, detector energy resolution, and background contamination will likely mask the modest shift in swap energy and the low‑energy survival enhancement. Therefore, the paper concludes that, despite the clear theoretical impact of small‑scale density fluctuations on collective oscillations, the neutrino burst alone is unlikely to serve as a reliable probe of fossil burning shells or similar fine‑grained structure in the supernova envelope.

In summary, the study provides a detailed numerical demonstration that a modest density bump can increase νₑ survival through a non‑adiabatic interplay of matter‑induced MSW resonances and collective self‑interaction effects. It clarifies the conditions under which the “synchronous MSW” and “regular precession” regimes dominate, and it highlights the difficulty of using supernova neutrino burst data to infer detailed density profiles. The work thus advances both the theoretical understanding of flavor evolution in dense astrophysical environments and the practical assessment of what information can realistically be extracted from future supernova neutrino observations.