Neutrinos propagating in dense neutrino media such as those in core-collapse supernovae can experience fast flavor conversions on scales much shorter than those expected in vacuum. It is believed that a necessary condition for the occurrence of fast modes is that the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. However, most of the state-of-the-art supernova simulations do not provide such detailed angular information and instead, consider only a few moments of neutrino angular distributions. We here propose an efficient method to use these available few moments to search for fast modes in supernova simulations. Our method, which is based on searching for crossings in the angular distributions, can work with any number of moments provided by the simulation though a larger number of crossings can be captured when higher moments are available.

We study collective oscillations of Majorana neutrinos in some of the most extreme astrophysical sites such as neutron star merger remnants and magneto-rotational core-collapse supernovae which include dense neutrino media in the presence of strong magnetic fields. We show that neutrinos can reach flavor equilibrium if neutrino transition magnetic moment $\mu_\nu$ is strong enough, namely when $\mu_\nu/\mu_{\rm{B}} \gtrsim 10^{-14}-10^{-15}$ with $\mu_{\rm{B}}$ being the Bohr magneton. This sort of flavor equilibrium, which is not necessarily flavor equipartition, can occur on (short) scales determined by the strength of the magnetic term. Our findings can have interesting implications for the physics of such violent astrophysical environments.

A dense neutrino gas, such as the one anticipated in the supernova environment, can experience fast neutrino flavor conversions on scales much shorter than those expected in vacuum probably provided that the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. We perform a detailed investigation of the neutrino angular distributions obtained by solving the Boltzmann equations for fixed matter profiles of some representative snapshots during the post-bounce phase of core-collapse supernovae in multidimensional calculations of an $11.2\mathrm{M}_{\odot}$ and a $27\mathrm{M}_{\odot}$ progenitor models. Although the $11.2\mathrm{M}_{\odot}$ model features $\nu_e - \bar\nu_e$ angular crossings and the associated fast modes at different time snapshots, the $27\mathrm{M}_{\odot}$ model does not show any crossings within the decoupling region. We show that this can be understood by studying the multipole components of the neutrino distributions. In fact, there is a higher chance for the occurrence of $\nu_e - \bar\nu_e$ angular crossings for the zones where the multipole components of the neutrino distributions are strong enough. We also show that there can exist more than one crossings between the angular distributions of $\nu_e$ and $\bar{\nu}_e$. In addition, apart from the crossings within the neutrino decoupling region, there is a class of $\nu_e-\bar\nu_e$ angular crossings which appears very deep inside the proto-neutron star.

The dense neutrino medium in a core-collapse supernova or a neutron-star merger event can experience fast flavor conversions on time/distance scales that are much smaller than those of vacuum oscillations. It is believed that fast neutrino flavor transformation occurs in the region where the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. We present the first study of this crossing phenomenon and the fast neutrino flavor conversions in multidimensional (multi-D) supernova models. We examine the neutrino distributions obtained by solving the Boltzmann transport equation for several fixed profiles which are representative snapshots taken from separate 2D and 3D supernova simulations with an $11.2 M_\odot$ progenitor model. Our research shows that the spherically asymmetric patterns of the $\nu_e$ and $\bar\nu_e$ fluxes in multi-D models can assist the appearance of the crossing between the $\nu_e$ and $\bar\nu_e$ angular distributions. In the models that we have studied, there exist unstable neutrino oscillation modes in and beyond the neutrino decoupling region which have amplitude growth rates as large as an $e$-fold per nanosecond. This finding can have important consequences for the explosion mechanism, nucleosynthesis, and neutrino signals of core-collapse supernovae.

We present a numerical survey of the nonlinear flavor development of dense neutrino gases. This study is based on the stationary, two-dimensional ($x$ and $z$), two-beam, monochromatic neutrino line model with a periodic boundary condition along the $x$ direction. Similar to a previous work, we find that small-scale flavor structures can develop in a neutrino gas even if the physical conditions are nearly homogeneous along the $x$ axis initially. The power diffusion from the large-scale to small-scale structures increases with the neutrino density and helps to establish a semi-exponential dependence of the magnitudes of the Fourier moments on the corresponding wave numbers. The overall flavor conversion probabilities in the neutrino gases with small initial sinusoidal perturbations reach certain equilibrium values at large distances which are mainly determined by the neutrino-antineutrino asymmetry. Similar phenomena also exist in a neutrino gas with a localized initial perturbation, albeit only inside an expanding flavor conversion region. Our work suggests that a statistical treatment may be possible for the collective flavor oscillations of a dense neutrino gas in a multi-dimensional environment.

A dense neutrino medium could experience self-induced flavor conversions on relatively small scales in the presence of the so-called fast flavor conversion modes. Owing to the fact that fast conversion scales could be much smaller than the ones of the traditional collective neutrino oscillations, it has been speculated that fast modes could lead to some sort of flavor decoherence/equilibrium. We study the evolution of fast modes in the nonlinear regime and we show that not only fast modes are not guaranteed to lead to flavor equilibrium, but also they could lead to some sort of collective neutrino oscillations but on scales determined by neutrino number density. For the $\nu_e$ dominated case, we observe large amplitude collective oscillations, whereas a sort of flavor stabilization is reached for the $\bar{\nu}_e$ dominated case.

The flavor conversion of a neutrino usually occurs at densities $\lesssim G_F^{-1} \omega$, whether in the ordinary matter or the neutrino medium, and on time/distance scales of order $\omega^{-1}$, where $G_F$ is the Fermi weak coupling constant and $\omega$ is the vacuum oscillation frequency of the neutrino. In 2005 Sawyer and more recently both he and other groups have shown that neutrino flavor conversions can occur on scales much shorter than $\omega^{-1}$ in a very dense, anisotropic neutrino gas such as that in a core-collapse supernova or a binary neutron star merger. It has also been suggested that the fast neutrino flavor conversion require a crossing in the electron lepton number (ELN) angular distribution or spectrum. We study a few simple analytical models with which we elucidate the roles of the dense matter and the ELN spectrum crossing in neutrino fast flavor conversions. We show that a large matter density can induce fast neutrino flavor conversions in certain astrophysical scenarios such as at the early epoch of a core-collapse supernova.

A dense neutrino medium such as that inside a core-collapse supernova can experience collective flavor conversion or oscillations because of the neutral-current weak interaction among the neutrinos. This phenomenon has been studied in a restricted, stationary supernova model which possesses the (spatial) spherical symmetry about the center of the supernova and the (directional) axial symmetry around the radial direction. Recently it has been shown that these spatial and directional symmetries can be broken spontaneously by collective neutrino oscillations. In this paper we analyze the neutrino flavor instabilities in a time-dependent supernova model. Our results show that collective neutrino oscillations start at approximately the same radius in both the stationary and time-dependent supernova models unless there exist very rapid variations in local physical conditions on timescales of a few microseconds or shorter. Our results also suggest that collective neutrino oscillations can vary rapidly with time in the regimes where they do occur which need to be studied in time-dependent supernova models.

We use the quantum Monte Carlo (QMC) techniques to calculate the static structure function $S(q)$ of a one-component ion lattice and use it to calculate the thermal conductivity $\kappa$ of high-density solid matter expected in the neutron star crust. By making detailed comparisons with the results for the thermal conductivity obtained using standard techniques based on the one-phonon approximation (OPA) valid at low temperature, and the multi-phonon harmonic approximation expected to be valid over a wide range of temperatures, we asses the temperature regime where $S(q)$ from QMC can be used directly to calculate $\kappa$. We also compare the QMC results to those obtained using classical Monte Carlo to quantitatively asses the magnitude of the quantum corrections. We find that quantum effects became relevant for the calculation of $\kappa$ at temperature $T \lesssim 0.3 ~\Omega_\mathrm{P}$, where $\Omega_\mathrm{P}$ is the ion plasma frequency. At $T \simeq 0.1 ~\Omega_\mathrm{P}$ the quantum effects suppress $\kappa$ by about $30\%$. The comparison with the results of the OPA indicates that dynamical information beyond the static structure is needed when $T \lesssim 0.1~ \Omega_\mathrm{P}$. These quantitative comparisons help to establish QMC as a viable technique to calculate $\kappa$ at moderate temperatures in the range $T=0.1-1~\Omega_\mathrm{P}$ of relevance to the study of accreting neutron stars. This finding is especially important because QMC is the only viable technique so far for calculating $\kappa$ in multi-component systems at low-temperatures.

Neutrino flavor oscillations in the presence of ambient neutrinos is nonlinear in nature which leads to interesting phenomenology that has not been well understood. It was recently shown that, in the two-dimensional, two-beam neutrino Line model, the inhomogeneous neutrino oscillation modes on small distance scales can become unstable at larger neutrino densities than the homogeneous mode does. We develop a numerical code to solve neutrino oscillations in the multi-angle/beam Line model with a continuous neutrino angular distribution. We show that the inhomogeneous oscillation modes can occur at even higher neutrino densities in the multi-angle model than in the two-beam model. We also find that the inhomogeneous modes on sufficiently small scales can be unstable at smaller neutrino densities with ambient matter than without, although a larger matter density does shift the instability region of the homogeneous mode to higher neutrino densities in the Line model as it does in the one-dimensional supernova Bulb model. Our results suggest that the inhomogeneous neutrino oscillation modes can be difficult to treat numerically because the problem of spurious oscillations becomes more severe for oscillations on smaller scales.