Neutrino-neutrino refraction causes self-induced flavor conversion in dense neutrino fluxes. For the first time, we include the azimuth angle of neutrino propagation as an explicit variable and find a new generic multi-azimuth-angle (MAA) instability which, for simple spectra, occurs in the normal neutrino mass hierarchy. Matter suppression of this instability in supernovae requires larger densities than the traditional bimodal case. The new instability shows explicitly that solutions of the equations for collective flavor oscillations need not inherit the symmetries of initial or boundary conditions. This change of paradigm requires reconsideration of numerous results in this field.

It has been recently pointed out that removing the axial symmetry in the "multi-angle effects" associated with the neutrino-neutrino interactions for supernova (SN) neutrinos, a new multi-azimuthal-angle (MAA) instability would arise. In particular, for a flux ordering $F_{\nu_e} > F_{\bar\nu_e} > F_{\nu_x}$, as expected during the SN accretion phase, this instability occurs in the normal neutrino mass hierarchy. However, during this phase the ordinary matter density can be larger than the neutrino one, suppressing the self-induced conversions. At this regard, we investigate the matter suppression of the MAA effects, performing a linearized stability analysis of the neutrino equations of motion, in the presence of realistic SN density profiles. We compare these results with the numerical solution of the SN neutrino non-linear evolution equations. Assuming axially symmetric distributions of neutrino momenta we find that the large matter term strongly inhibits the MAA effects. In particular, the hindrance becomes stronger including realistic forward-peaked neutrino angular distributions. As a result, in our model for a $10.8$ $M_{\odot}$ iron-core SNe, MAA instability does not trigger any flavor conversion during the accretion phase. Instead, for a $8.8$ $M_{\odot}$ O-Ne-Mg core SN model, with lower matter density profile and less forward-peaked angular distributions, flavor conversions are possible also at early times.

We construct a simple example for self-induced flavor conversion in dense neutrino gases showing new solutions that violate the symmetries of initial conditions. Our system consists of two opposite momentum modes 1 and 2, each initially occupied with equal densities of nu_e and anti-nu_e. Restricting solutions to symmetry under 1 <-> 2 allows for the usual bimodal instability ("flavor pendulum") in the inverted neutrino mass hierarchy (IH) and stability (no self-induced flavor conversion) in the normal hierarchy (NH). Lifting this symmetry restriction allows for a second pendulum-like solution that occurs in NH where the modes 1 and 2 swing in opposite directions in flavor space. Any small deviation from 1-2 symmetry in the initial condition triggers the new instability in NH. This effect corresponds to the recently identified multi-azimuth angle (MAA) instability of supernova neutrino fluxes. Both cases show explicitly that solutions of the equations of collective flavor oscillations need not inherit the symmetries of initial conditions, although this has been universally assumed.

The dense neutrino flux streaming from a core-collapse supernova can undergo self-induced flavor conversion caused by neutrino-neutrino refraction. Numerical studies of these nonlinear effects are challenging because representing the neutrino radiation field by discrete energy and angle bins can easily lead to unphysical solutions. In particular, if the number of angle bins N_a is too small, flavor conversion begins too deep and produces completely spurious results. At the same time, N_a=1 (single-angle approximation) can be a good proxy for the N_a -> infinity limit. Based on a linearized stability analysis, we explain some of the puzzling effects of discrete angle distributions.