Neutrino flavor conversions may dramatically affect the inner working of compact astrophysical objects as well as the synthesis of the heavier elements. We present the first sophisticated numerical solution of the neutrino flavor conversion within a (2+1+1) dimensional setup: we include the advective term in the neutrino equations of motion and track the flavor evolution in two spatial dimensions, one angular variable, and time. Notably, the advective term hinders the development of neutrino pairwise conversions, if the conditions favoring such conversions (i.e., crossings between the angular distributions of $\nu_e$ and $\bar\nu_e$ or a non-negligible flux of neutrinos traveling backward with respect to the main propagation direction) exist for time scales shorter than the typical time scale of the advective term. As a consequence, fast pairwise conversions can only occur when the conditions favoring flavor conversions are self-sustained and global, such as the ones induced by the lepton emission self-sustained asymmetry (LESA) in core-collapse supernovae. Our work highlights the major impact of the dynamical evolution of the neutrino field on the growth of flavor instabilities and the strong interplay between classical and quantum effects. Critical limitations of the linear stability analysis, used to predict neutrino flavor instabilities, are also pointed out.

Neutrinos offer a window to physics beyond the Standard Model. In particular, high-energy astrophysical neutrinos, with TeV-PeV energies, may provide evidence of new, "secret" neutrino-neutrino interactions that are stronger than ordinary weak interactions. During their propagation over cosmological distances, high-energy neutrinos could interact with the cosmic neutrino background via secret interactions, developing characteristic energy-dependent features in their observed energy distribution. For the first time, we look for signatures of secret neutrino interactions in the diffuse flux of high-energy astrophysical neutrinos, using 6 years of publicly available IceCube High Energy Starting Events (HESE). We find no significant evidence for secret neutrino interactions, but place competitive upper limits on the coupling strength of the new mediator through which they occur, in the mediator mass range of 1-100 MeV.

Beyond-the-Standard-Model interactions of neutrinos among themselves -- {\it secret interactions} -- in the supernova core may prevent the shock revival, halting the supernova explosion. Besides, if supernova neutrinos en route to Earth undergo secret interactions with relic neutrinos, the neutrino burst reaching Earth may be down-scattered in energy, falling below the detection threshold. We probe secret neutrino interactions through supernova neutrinos and apply our findings to the supernova SN 1987A. We place the most stringent bounds on secret interactions occurring through a new mediator with mass between 10 MeV and 15 GeV.

Neutrino fast pairwise conversions have been postulated to occur in the dense core of a core-collapse supernova (SN), possibly having dramatic consequences on the SN mechanism and the observable neutrino signal. One crucial condition favoring pairwise conversions is the presence of crossings between the electron neutrino and antineutrino angular distributions (i.e., electron neutrino lepton number crossings, ELN crossings). A stationary and spherically symmetric SN toy-model is constructed to reproduce the development of the neutrino angular distributions in the dense SN core in the absence of perturbations induced by hydrodynamical instabilities. By iteratively solving the neutrino Boltzmann equations including the collisional term, our model predicts that ELN crossings can develop only in the proximity of the decoupling region and for a sharp radial evolution of the baryon density, when the electron neutrino and antineutrino number densities are comparable. Such conditions are likely to occur only in the late SN stages. Interestingly, flavor instabilities induced by spatial or temporal perturbations are unlikely to generate ELN crossings dynamically within our simplified setup.

We consider the flavor transformation of neutrinos through oscillatory matter profiles. We show that the neutrino oscillation Hamiltonian in this case describes a Rabi system with an infinite number of Rabi modes. We further show that, in a given physics problem, the majority of the Rabi modes have too small amplitudes to be relevant. We also go beyond the rotating wave approximation and derive the relative detuning of the Rabi resonance when multiple Rabi modes with small amplitudes are present. We provide an explicit criterion of whether an off-resonance Rabi mode can affect the parametric flavor transformation of the neutrino.

We perform a self-consistent calculation of collective neutrino oscillations including the effect of back scattered neutrinos (Halo effect) in the `single-angle' approximation, within a spherically symmetric supernova model. We find that due to the Halo effect the onset of flavor transformations is pushed to smaller radii, by a few kilometers. The celebrated phenomenon of the spectral split is found to be robust under the present inclusion of the Halo effect.

We perform a numerical study of the matter-neutrino resonance in a multi-angle calculation in the vicinity of an accretion disk. We assume thermally distributed neutrino and anti-neutrino fields emitted by two-dimensional disk that is homogeneous and isotropic; the electrons are assumed to be at constant density. We compare the the result of this computation to that obtained using single-angle approximation. We investigate the robustness of matter-neutrino resonance in environment surrounding accretion disks by progressively relaxing the single angle approximation. We find that the multi-angle results in the present simplified model do not support a robust resonance mechanism as suggested by the single angle treatment. We also discuss the context under which matter-neutrino resonance may be important in future studies.

We investigate the stability, with respect to spatial inhomogeneity, of a two-dimensional dense neutrino gas. The system exhibits growth of seed inhomogeneity due to nonlinear coherent neutrino self-interactions. In the absence of incoherent collisional effects, we observe a dependence of this instability growth rate on the neutrino mass spectrum: the normal neutrino mass hierarchy exhibits spatial instability over a larger range of neutrino number density compared to that of the inverted case. We further consider the effect of elastic incoherent collisions of the neutrinos with a static background of heavy, nucleon-like scatterers. At small scales, the growth of flavor instability can be suppressed by collisions. At large length scales we find, perhaps surprisingly, that for inverted neutrino mass hierarchy incoherent collisions fail to suppress flavor instabilities, independent of the coupling strength.

We study indirect detection signals from solar annihilation of dark matter (DM) particles into light right-handed (RH) neutrinos with a mass in a $1-5$ GeV range. These RH neutrinos can have a sufficiently long lifetime to allow them to decay outside the Sun and their delayed decays can result in a signal in gamma rays from the otherwise `dark' solar direction, and also a neutrino signal that is not suppressed by the interactions with solar medium. We find that the latest Fermi-LAT and IceCube results place limits on the gamma ray and neutrino signals, respectively. Combined photon and neutrino bounds can constrain the spin-independent DM-nucleon elastic scattering cross section better than direct detection experiments for DM masses from 200 GeV up to several TeV. The bounds on spin-dependent scattering are also much tighter than the strongest limits from direct detection experiments.

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.

A dense neutrino medium can experience collective flavor oscillations through nonlinear neutrino-neutrino refraction. To make this multi-dimensional flavor transport problem more tractable, all existing studies have assumed certain symmetries (e.g., the spatial homogeneity and directional isotropy in the early universe) to reduce the dimensionality of the problem. In this work we show that, if both the directional and spatial symmetries are not enforced in the neutrino line model, collective oscillations can develop in the physical regimes where the symmetry-preserving oscillation modes are stable. Our results suggest that collective neutrino oscillations in real astrophysical environments (such as core-collapse supernovae and black-hole accretion discs) can be qualitatively different from the predictions based on existing models in which spatial and directional symmetries are artificially imposed.

We demonstrate a multipole expansion method to calculate collective neutrino oscillations in supernovae using the neutrino bulb model. We show that it is much more efficient to solve multi-angle neutrino oscillations in multipole basis than in angle basis. The multipole expansion method also provides interesting insights into multi-angle calculations that were accomplished previously in angle basis.

We demonstrate the non-negligible effect of transition magnetic moments on three-flavor collective oscillations of Majorana neutrinos in the core of type-II supernovae, within the single-angle approximation. We argue that data from a galactic supernova in conjunction with terrestrial experiments can potentially give us clues about the non-zero nature of neutrino transition magnetic moments if these are Majorana fermions, even if their values are as small as those predicted by the Standard Model augmented by nonzero neutrino Majorana masses.

We study the effect of Majorana transition magnetic moments on the flavor evolution of neutrinos and antineutrinos inside the core of Type-II supernova explosions. We find non-trivial collective oscillation effects relating neutrinos and antineutrinos of different flavors, even if one restricts the discussion to Majorana transition electromagnetic moment values that are not much larger than those expected from standard model interactions and nonzero neutrino Majorana masses. This appears to be, to the best of our knowledge, the only potentially observable phenomenon sensitive to such small values of Majorana transition magnetic moments. We briefly comment on the effect of Dirac transition magnetic moments and on the consequences of our results for future observations of the flux of neutrinos of different flavors from a nearby supernova explosion.