We investigate how sterile neutrinos with a range of masses influence cosmology in MeV-scale reheating temperature scenarios. By computing the production of sterile neutrinos through the combination of mixing and scattering in the early Universe, we find that light sterile neutrinos, with masses and mixings as inferred from short-baseline neutrino oscillation experiments, are consistent with big-bang nucleosynthesis (BBN) and cosmic microwave background (CMB) radiation for the reheating temperature of ${\cal O}(1)$ MeV if the parent particle responsible for reheating decays into electromagnetic components (radiative decay). In contrast, if the parent particle mainly decays into hadrons (hadronic decay), the bound from BBN becomes more stringent. In this case, the existence of the light sterile neutrinos can be cosmologically excluded, depending on the mass and the hadronic branching ratio of the parent particle.

From a theoretical point of view, there is a strong motivation to consider an MeV-scale reheating temperature induced by long-lived massive particles with masses around the weak scale, decaying only through gravitational interaction. In this study, we investigate lower limits on the reheating temperature imposed by big-bang nucleosynthesis assuming both radiative and hadronic decays of such massive particles. For the first time, effects of neutrino self-interactions and oscillations are taken into account in the neutrino thermalization calculations. By requiring consistency between theoretical and observational values of light element abundances, we find that the reheating temperature should conservatively be $T_{\rm RH} \gtrsim 1.8$ MeV in the case of the 100% radiative decay, and $T_{\rm RH} \gtrsim$ 4-5 MeV in the case of the 100% hadronic decays for particle masses in the range of 10 GeV to 100 TeV.

In supernovae neutrinos are emitted from a region with a width $r_{\rm eff}$ of a few kilometers (rather than from a surface of infinitesimal width). We study the effect of integration (averaging) over such an extended emission region on collective oscillations. The averaging leads to additional suppression of the correlation (off-diagonal element of the density matrix) by a factor $ \sim 1/r_{\rm eff} V_e \sim 10^{-10}$, where $V_e$ is the matter potential. This factor enters the initial condition for further collective oscillations and, consequently, leads to a delay of the strong flavour transitions. We justify and quantify this picture using a simple example of collective effects in two intersecting fluxes. We have derived the evolution equation for the density matrix elements integrated over the emission region and solved it both numerically and analytically. For the analytic solution we have used linearized equations. We show that the delay of the development of the instability and the collective oscillations depends on the suppression factor due to the averaging (integration) logarithmically. If the instability develops inside the production region, the integration leads not only to a delay but also to a modification of the exponential grow.

Neutrino-neutrino interactions can lead to collective flavor conversion in the dense parts of a core collapse supernova. Growing instabilities that lead to collective conversions have been studied intensely in the limit of two-neutrino species and occur for inverted mass ordering in the case of a perfectly spherical supernova. We examine two simple models of colliding and intersecting neutrino beams and show, that for three neutrino species instabilities exist also for normal mass ordering even in the case of a fully symmetric system. Whereas the instability for inverted mass ordering is associated with $\Delta m_{31}^2$, the new instability we find for normal mass ordering is associated with $\Delta m_{21}^2$. As a consequence, the growth rate of these new instabilities for normal ordering is smaller by about an order of magnitude compared to the rates of the well studied case of inverted ordering.

We study several methods for timing the neutrino signal of a galactic supernova for different detectors via Monte Carlo simulations. We find that, for the methods we studied, both Hyper-Kamiokande and IceCube can reach precisions of $\sim1\,$ms for the neutrino burst, while a potential IceCube Gen2 upgrade will reach sub-ms precision. In the case of a failed SN, we find that detectors like SK and JUNO can reach precisions of $\sim0.1\,$ms while HK could potentially reach a resolution of $\sim 0.01\,$ms so that the impact of the BH formation process itself becomes relevant. Two possible applications for this are the triangulation of a (failed) SN as well as the possibility to constrain neutrino masses via a time-of-flight measurement using a potential gravitational wave signal as reference.

Motivated by the discovery of the first high-energy astrophysical neutrino source, the blazar TXS 0506+056, we revisit the IceCube flavor ratio analysis. Assuming large statistics from identified blazars, collected in the forthcoming years by the IceCube detector and its successor IceCube-Gen2, we demonstrate that the constraints on several new physics scenarios in which the baseline dependent terms in neutrino oscillation probabilities are not averaged, can be improved. As a representative case, we consider pseudo-Dirac neutrinos while neutrino decay is also discussed.

Collective oscillations of supernova neutrinos above the neutrino sphere can be completely described by the propagation of individual neutrinos in external potentials and are in this sense a linear phenomenon. An effective theory of collective oscillations can be developed based on certain assumptions about time dependence of these potentials. General conditions for strong flavor transformations are formulated and these transformations can be interpreted as parametric resonance effects induced by periodic modulations of the potentials. We study a simplified and solvable example, where a probe neutrino is propagating in a flux of collinear neutrinos, such that $\nu \nu-$ interactions in the flux are absent. Still, this example retains the main feature - the coherent flavor exchange. Properties of the parametric resonance are studied, and it is shown that integrations over energies and emission points of the flux neutrinos suppress modulations of the potentials and therefore strong transformations. The transformations are also suppressed by changes in densities of background neutrinos and electrons.

Sterile neutrinos produced through oscillations are a well motivated dark matter candidate, but recent constraints from observations have ruled out most of the parameter space. We analyze the impact of new interactions on the evolution of keV sterile neutrino dark matter in the early Universe. Based on general considerations we find a mechanism which thermalizes the sterile neutrinos after an initial production by oscillations. The thermalization of sterile neutrinos is accompanied by dark entropy production which increases the yield of dark matter and leads to a lower characteristic momentum. This resolves the growing tensions with structure formation and X-ray observations and even revives simple non-resonant production as a viable way to produce sterile neutrino dark matter. We investigate the parameters required for the realization of the thermalization mechanism in a representative model and find that a simple estimate based on energy- and entropy conservation describes the mechanism well.

We consider several aspects related to the form, derivation and applications of the Liouville equation (LE) for flavour evolution of neutrinos. To take into account the quantum nature of neutrinos we derive the evolution equation for the matrix of densities using wave packets instead of Wigner functions. The obtained equation differs from the standard LE by an additional term which is proportional to the difference of group velocities. We show that this term describes loss of the propagation coherence in the system. In absence of inelastic collisions, the LE can be reduced to a single derivative equation over a trajectory coordinate. Additional time and spacial dependence may steam from initial (production) conditions. The transition from single neutrino evolution to the evolution of a neutrino gas is considered.

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

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.

Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave Background and large scale structure. This paper discusses how charging sterile neutrinos under a new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available data which is way better than the standard $\Lambda$CDM model with one additional fully thermalized sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the locally measured value.

Sterile neutrinos are thermalised in the early Universe via oscillations with the active neutrinos for certain mixing parameters. The most detailed calculation of this thermalisation process involves the solution of the momentum-dependent quantum kinetic equations, which track the evolution of the neutrino phase space distributions. Until now the collision terms in the quantum kinetic equations have always been approximated using equilibrium distributions, but this approximation has never been checked numerically. In this work we revisit the sterile neutrino thermalisation calculation using the full collision term, and compare the results with various existing approximations in the literature. We find a better agreement than would naively be expected, but also identify some issues with these approximations that have not been appreciated previously. These include an unphysical production of neutrinos via scattering and the importance of redistributing momentum through scattering, as well as details of Pauli blocking. Finally, we devise a new approximation scheme, which improves upon some of the shortcomings of previous schemes.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. Such sterile neutrinos are incompatible with cosmology because they suppress structure formation unless they can be prevented from thermalising in the early Universe or removed by subsequent decay or annihilation. Here we present a novel scenario in which both sterile neutrinos and dark matter are coupled to a new, light pseudoscalar. This can prevent thermalisation of sterile neutrinos and make dark matter sufficiently self-interacting to have an impact on galactic dynamics and possibly resolve some of the known problems with the standard cold dark matter scenario. Even more importantly it leads to a strongly self-interacting plasma of sterile neutrinos and pseudoscalars at late times and provides an excellent fit to CMB data. The usual cosmological neutrino mass problem is avoided by sterile neutrino annihilation to pseudoscalars. The preferred value of $H_0$ is substantially higher than in standard $\Lambda$CDM and in much better agreement with local measurements.

Neutrino-neutrino refraction can lead to non-periodic flavor oscillations in dense neutrino gases, and it has been hypothesized that some solutions are chaotic in nature. This is of particular interest in the case of neutrino emission from core-collapse supernovae where the measurement of the spectral shape for different flavors can provide crucial information about both neutrino physics and the physical conditions close to the proto-neutron star. Whether a system is chaotic or not can be assessed by the Lyapunov exponents which quantify the rate of divergence of nearby trajectories in the system. We have done a numerical case study for a simple toy model of two neutrino flavors with two momentum states traveling against each other which is known to exhibit flavor transition instabilities. We find the leading Lyapunov exponent to be positive in all cases, confirming the chaoticity of the system for both the normal and the inverted neutrino mass hierarchy. However, more Lyapunov exponents were approximately zero in the inverted hierarchy compared to the normal which has implications for the stability of the system. To investigate this, we have calculated a generalized set of normal modes, the so-called covariant Lyapunov vectors. The covariant Lyapunov vectors associated with vanishing Lyapunov exponents showed the existence of marginally stable directions in phase space for some cases. While our analysis was done for a toy model example, it should work equally well for more realistic cases of neutrinos streaming from a proto-neutron star and provide valuable insight into the nature of the flavor instability. We finally stress that our approach captures many more properties of the physical system than the linear stability analyses which have previously been performed.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. However, such neutrinos seem incompatible with cosmology because they have too large an impact on cosmic structure formation. Here we show that new interactions in the sterile neutrino sector can prevent their production in the early Universe and reconcile short baseline oscillation experiments with cosmology.

While the cosmic baryon asymmetry has been measured at high accuracy to be 6.1*10^-10, a corresponding lepton asymmetry could be as large as 10^-2 if it hides in the neutrino sector. It has been known for some time that a large asymmetry can be generated from a small initial asymmetry given the existence of a sterile neutrino with a mass less than the mass of the active neutrino. While the magnitude of the final lepton asymmetry is deterministic, its sign has been conjectured to be chaotic in nature. This has been proven in the single momentum approximation, also known as the quantum rate equations, but has up to now not been established using the full momentum dependent quantum kinetic equations. Here we investigate this problem by solving the quantum kinetic equations for a system of 1 active and 1 sterile neutrino on an adaptive grid. We show that by increasing the resolution oscillations in the lepton asymmetry are eliminated and that the sign of the lepton asymmetry is in fact deterministic. This paper also serves as a launch paper for the adaptive solver LASAGNA which is available at http://users-phys.au.dk/steen.