Collective pair conversion $\nu_e\bar\nu_e\leftrightarrow \nu_{x}\bar\nu_{x}$ by forward scattering, where $x=\mu$ or $\tau$, may be generic for supernova neutrino transport. Depending on the local angular intensity of the electron lepton number carried by neutrinos, the conversion rate is "fast," i.e., of the order of $\sqrt{2}G_{\rm{F}}(n_{\nu_e}{-}\,n_{\bar\nu_e})\gg\Delta m^2_{\rm atm}/2E$. We present a novel approach to understand these phenomena: A dispersion relation for the frequency and wave number $(\Omega,\bf{K})$ of disturbances in the mean field of $\nu_e\nu_x$ flavor coherence. Run-away solutions occur in "dispersion gaps," i.e., in "forbidden" intervals of $\Omega$ and/or $\bf{K}$ where propagating plane waves do not exist. We stress that the actual solutions also depend on the initial and/or boundary conditions which need to be further investigated.

Cosmic ray accelerators like supernova and hypernova remnants in star forming galaxies are one of the most plausible sources of the IceCube observed diffuse astrophysical neutrinos. The neutrino producing hadronic processes will also produce a diffuse gamma ray flux, constrained by the Fermi-LAT bounds. The fact that point sources like blazars also contribute to the diffuse gamma ray flux implies large gamma opacity of the neutrino sources. Indeed, for these high redshift star forming galaxies the gamma absorption during the intergalactic propagation can be significant. In addition, large gamma attenuation inside these extreme source galaxies can reduce the cascade component of the diffuse flux. Under the current astrophysical uncertainties affecting these absorptions processes, we find the associated diffuse gamma ray flux can remain compatible with the current Fermi-LAT bounds.

Neutrino-neutrino refraction in dense media can cause self-induced flavor conversion triggered by collective run-away modes of the interacting flavor oscillators. The growth rates were usually found to be of order a typical vacuum oscillation frequency $\Delta m^2/2E$. However, even in the simple case of a $\nu_e$ beam interacting with an opposite-moving $\bar\nu_e$ beam, and allowing for spatial inhomogeneities, the growth rate of the fastest-growing Fourier mode is of order $\mu=\sqrt{2} G_{\rm F} n_{\nu}$, a typical $\nu$--$\nu$ interaction energy. This growth rate is much larger than the vacuum oscillation frequency and gives rise to flavor conversion on a much shorter time scale. This phenomenon of "fast flavor conversion" occurs even for vanishing $\Delta m^2/2E$ and thus does not depend on energy, but only on the angle distributions. Moreover, it does not require neutrinos to mix or to have masses, except perhaps for providing seed disturbances. We also construct a simple homogeneous example consisting of intersecting beams and study a schematic supernova model proposed by Ray Sawyer, where $\nu_e$ and $\bar\nu_e$ emerge with different zenith-angle distributions, the key ingredient for fast flavor conversion. What happens in realistic astrophysical scenarios remains to be understood.

Neutrino flavor evolution in core-collapse supernovae, neutron-star mergers, or the early universe is dominated by neutrino-neutrino refraction, often spawning "self-induced flavor conversion", i.e., shuffling of flavor among momentum modes. This effect is driven by collective run-away modes of the coupled "flavor oscillators" and can spontaneously break the initial symmetries such as axial symmetry, homogeneity, isotropy, and even stationarity. Moreover, the growth rates of unstable modes can be of the order of the neutrino-neutrino interaction energy instead of the much smaller vacuum oscillation frequency: self-induced flavor conversion does not always require neutrino masses. We illustrate these newly found phenomena in terms of simple toy models. What happens in realistic astrophysical settings is up to speculation at present.

Self-induced flavor conversion of supernova (SN) neutrinos is a generic feature of neutrino-neutrino dispersion. The corresponding run-away modes in flavor space can spontaneously break the original symmetries of the neutrino flux and in particular can spontaneously produce small-scale features as shown in recent schematic studies. However, the unavoidable "multi-angle matter effect" shifts these small-scale instabilities into regions of matter and neutrino density which are not encountered on the way out from a SN. The traditional modes which are uniform on the largest scales are most prone for instabilities and thus provide the most sensitive test for the appearance of self-induced flavor conversion. As a by-product we clarify the relation between the time evolution of an expanding neutrino gas and the radial evolution of a stationary SN neutrino flux. Our results depend on several simplifying assumptions, notably stationarity of the solution, the absence of a "backward" neutrino flux caused by residual scattering, and global spherical symmetry of emission.

IceCube has measured a diffuse astrophysical flux of TeV-PeV neutrinos. The most plausible sources are unique high energy cosmic ray accelerators like hypernova remnants (HNRs) and remnants from gamma ray bursts in star-burst galaxies, which can produce primary cosmic rays with the required energies and abundance. In this case, however, ordinary supernova remnants (SNRs), which are far more abundant than HNRs, produce a comparable or larger neutrino flux in the ranges up to 100-150 TeV energies, implying a spectral break in the IceCube signal around these energies. The SNRs contribution in the diffuse flux up to these hundred TeV energies provides a natural baseline and then constrains the expected PeV flux.

Mixing lepton doublets of the Standard Model can lead to lepton flavour asymmetries in the Early Universe. We present a diagrammatic representation of this recently identified source of $CP$ violation and elaborate in detail on the correlations between the lepton flavours at different temperatures. For a model where two sterile right-handed neutrinos generate the light neutrino masses through the see-saw mechanism, the lower bound on reheat temperatures in accordance with the observed baryon asymmetry turns out to be $\gsim 1.2\times 10^9\,{\rm GeV}$. With three right-handed neutrinos, substantially smaller values are viable. This requires however a tuning of the Yukawa couplings, such that there are cancellations between the individual contributions to the masses of the light neutrinos.