Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $\nu_e$ and $\bar\nu_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $\nu_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

We present a linear stability analysis of the fast-pairwise neutrino flavor conversion based on a result of our latest axisymmetric core-collapse supernova (CCSN) simulation with full Boltzmann neutrino transport. In the CCSN simulation, coherent asymmetric neutrino emissions of electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$), in which the asymmetry of $\nu_{\rm e}$ and $\bar{\nu}_{\rm e}$ is anti-correlated with each other, occur at almost the same time as the onset of aspherical shock expansion. We find that the asymmetric neutrino emissions play a crucial role on occurrences of fast flavor conversions. The linear analysis shows that unstable modes appear in both pre- and post-shock flows; for the latter they appear only in the hemisphere of higher $\bar{\nu}_{\rm e}$ emissions (the same hemisphere with stronger shock expansion). We analyze in depth the characteristics of electron-lepton-number (ELN) crossing by closely inspecting the angular distributions of neutrinos in momentum space. The ELN crossing happens in various ways, and the property depends on the radius: in the vicinity of neutron star, $\bar{\nu}_{\rm e}$ ($\nu_{\rm e}$) dominates over $\nu_{\rm e}$ ($\bar{\nu}_{\rm e}$) in the forward (backward) direction: at the larger radius the ELN crossing occurs in the opposite way. We also find that the non-radial ELN crossing occurs at the boundary between no ELN crossing and the radial one, which is an effect of genuine multi-D transport. Our findings indicate that the collective neutrino oscillation may occur more commonly in CCSNe and suggest that the CCSN community needs to accommodate these oscillations self-consistently in the modelling of CCSNe.

We make a strong case that the fast neutrino-flavor conversion, one of the collective flavor oscillation modes, commonly occurs in core-collapse supernovae (CCSNe). It is confirmed in the numerical data obtained in realistic simulations of CCSNe but the argument is much more generic and applicable universally: the coherent neutrino-nucleus scattering makes the electron lepton number (ELN) change signs at some inward direction and trigger the flavor conversion in the outward direction in the pre-shock region. Although the ELN crossing is tiny and that is why it has eluded recognition so far, it is still large enough to induce the flavor conversion. Our findings will have an important observational consequences for CCSNe neutrinos.

Neutrinos are densely populated deep inside the core of massive stars after their gravitational collapse to produce supernova explosions and form compact stars such as neutron stars (NS) and black holes (BH). It has been considered that they may change their flavor identities through so-called fast-pairwise conversions induced by mutual forward scatterings. If that is really the case, the dynamics of supernova explosion will be influenced, since the conversion may occur near the neutrino sphere, from which neutrinos are effectively emitted. In this paper, we conduct a pilot study of such possibilities based on the results of fully self-consistent, realistic simulations of a core-collapse supernova explosion in two spatial dimensions under axisymmetry. As we solved the Boltzmann equations for neutrino transfer in the simulation not as a post-process but in real time, the angular distributions of neutrinos in momentum space for all points in the core at all times are available, a distinct feature of our simulations. We employ some of these distributions extracted at a few selected points and times from the numerical data and apply linear analysis to assess the possibility of the conversion. We focus on the vicinity of the neutrino sphere, where different species of neutrinos move in different directions and have different angular distributions as a result. This is a pilot study for a more thorough survey that will follow soon. We find no positive sign of conversion unfortunately at least for the spatial points and times we studied in this particular model. We hence investigate rather in detail the condition for the conversion by modifying the neutrino distributions rather arbitrarily by hand.

Collective neutrino oscillations are induced by the presence of neutrinos themselves. As such they are intrinsically nonlinear phenomena and are much more complex than linear counterparts such as the vacuum or MSW oscillations. They obey integro-differential equations, numerical solutions of which are also very challenging. If one focuses on the onset of the collective oscillations, on the other hand, the equations can be linearized and the technique of linear analysis can be employed. Unfortunately, however, it is well known that such an analysis, when applied with discretizations of continuous angular distributions, suffers from the appearance of so-called spurious modes, unphysical eigenmodes of the discretized linear equations. In this paper, we analyze in detail the origin of these unphysical modes and present a simple solution to this annoying problem. We have found that the spurious modes originate from the artificial production of pole singularities instead of a branch cut in the Riemann surface by the discretizations. The branching point singularities in the Riemann surface for the original undiscretized equations can be recovered by approximating the angular distributions with polynomials and then performing the integrals analytically. We demonstrate for some examples that this simple prescription removes the spurious modes indeed. We also propose an even simpler method: a piecewise linear approximation to the angular distribution. It is shown that the same methodology is applicable to the multi-energy case as well as to the dispersion relation approach that was proposed very recently.