We derive the evolution equations for a system of neutrinos interacting among themselves and with a matter background, based upon the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. This theoretical framework gives an (unclosed) set of first-order coupled integro-differential equations governing the evolution of the reduced density matrices. By employing the hierarchy, we first rederive the mean-field evolution equations for the neutrino one-body density matrix associated with a system of neutrinos and antineutrinos interacting with matter and with an anisotropic neutrino background. Then, we derive extended evolution equations to determine neutrino flavor conversion beyond the commonly used mean-field approximation. To this aim we include neutrino-antineutrino pairing correlations to the two-body density matrix. The inclusion of these new contributions leads to an extended evolution equation for the normal neutrino density and to an equation for the abnormal one involving the pairing mean-field. We discuss the possible impact of neutrino-antineutrino correlations on neutrino flavor conversion in the astrophysical and cosmological environments, and possibly upon the supernova dynamics. Our results can be easily generalized to an arbitrary number of neutrino families.

Non-Standard Interactions (NSI) of neutrinos with matter can significantly alter neutrino flavor evolution in supernovae with the potential to impact explosion dynamics, nucleosynthesis and the neutrinos signal. In this paper we explore, both numerically and analytically, the landscape of neutrino flavor transformation effects in supernovae due to NSI and find new, heretofore unseen transformation processes can occur. These new transformations can take place with NSI strengths well below current experimental limits. Within a broad swath of NSI parameter space we observe Symmetric and Standard Matter-Neutrino Resonances (MNRs) for supernovae neutrinos, a transformation effect previously only seen in compact object merger scenarios; in another region of the parameter space we find the NSI can induce neutrino collective effects in scenarios where none would appear with only the standard case of neutrino oscillation physics; and in a third region the NSI can lead to the disappearance of the high (H) density Mikheyev-Smirnov-Wolfenstein (MSW) resonance. Using a variety of analytical tools we are able to describe quantitatively the numerical results allowing us to partition the NSI parameter according to the transformation processes observed. Our results indicate non-standard interactions of supernova neutrinos provide a sensitive probe of Beyond the Standard Model physics complementary to present and future terrestrial experiments.

Core-collapse supernova neutrinos undergo a variety of phenomena when they travel from the high neutrino density region and large matter densities to the Earth. We perform analytical calculations of the supernova neutrino fluxes including collective effects due to the neutrino-neutrino interactions, the Mikheev-Smirnov-Wolfenstein (MSW) effect due to the neutrino interactions with the background matter and decoherence of the wave packets as they propagate in space. We predict the numbers of one- and two-neutron charged and neutral-current electron-neutrino scattering on lead events. We show that, due to the energy thresholds, the ratios of one- to two-neutron events are sensitive to the pinching parameters of neutrino fluxes at the neutrinosphere, almost independently of the presently unknown neutrino properties. Besides, such events have an interesting sensitivity to the spectral split features that depend upon the presence/absence of energy equipartition among neutrino flavors. Our calculations show that a lead-based observatory like the Helium And Lead Observatory (HALO) has the potential to pin down important characteristics of the neutrino fluxes at the neutrinosphere, and provide us with information on the neutrino transport in the supernova core.