We explore the energy and entropy transport as well as the lepton number variation induced from the mixing between electron and sterile neutrinos with keV mass in the supernova core. We develop a radial- and time-dependent treatment of the sterile-electron neutrino mixing, by including ordinary matter effects, reconversions between sterile and electron antineutrinos, as well as the collisional production of sterile particles for the first time. The dynamical feedback due to the production of sterile particles on the composition and thermodynamic properties of the core only leads to major implications for the supernova physics for large mixing angles ($\sin^2 2 \theta > 10^{-10}$). Our findings suggest that a self-consistent appraisal of the electron-sterile conversion physics in the supernova core would relax the bounds on the sterile neutrino mixing parameters reported in the literature for mixing angles smaller than $10^{-6}$, leaving the parameter space of the mass and mixing angles of sterile neutrinos relevant to dark matter searches unconstrained by supernovae.

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.

We propose the RES-NOVA project which will hunt neutrinos from core-collapse supernovae (SN) via coherent elastic neutrino-nucleus scattering (CE$\nu$NS) using an array of archaeological lead (Pb) based cryogenic detectors. The high CE$\nu$NS cross-section on Pb and the ultra-high radiopurity of archaeological Pb enable the operation of a high statistics experiment equally sensitive to all neutrino flavors with reduced detector dimensions in comparison with existing Neutrino Observatories, and easy scalability to larger detector volumes. RES-NOVA is planned to operate according to three phases with increasing detector volumes: (60 cm)$^3$, (140 cm)$^3$, and ultimately 15$\times$(140 cm)$^3$. It will be sensitive to SN bursts up to Andromeda with 5$\sigma$ sensitivity with already existing technology and will have excellent energy resolution with $1$ keV threshold. Within our Galaxy, it will be possible to discriminate core-collapse SNe from black hole forming collapses with no ambiguity even in the first phase of RES-NOVA. The average neutrino energy of all flavors, the SN neutrino light curve, and the total energy emitted in neutrinos can potentially be constrained with a precision of few $\%$ in the final detector phase. RES-NOVA will be sensitive to flavor-blind neutrinos from the diffuse SN neutrino background with an exposure of $620$ ton $\cdot$ y. The proposed RES-NOVA project has the potential to lay down the foundations for a new generation of neutrino telescopes, while relying on a very simple technological setup.

Neutrinos are unique probes of core-collapse supernova dynamics, especially in the case of black hole (BH) forming stellar collapses, where the electromagnetic emission may be faint or absent. By investigating two 3D hydrodynamical simulations of BH-forming stellar collapses of mass 40 and 75 M_sun, we identify the physical processes preceding BH formation through neutrinos, and forecast the neutrino signal expected in the existing IceCube and Super-Kamiokande detectors, as well as in the future generation DUNE facility. Prior to the abrupt termination of the neutrino signal corresponding to BH formation, both models develop episodes of strong and long-lasting activity by the spiral standing accretion shock instability (SASI). We find that the spiral SASI peak in the Fourier power spectrum of the neutrino event rate will be distinguishable at 3 sigma above the detector noise for distances up to O(30) kpc in the most optimistic scenario, with IceCube having the highest sensitivity. Interestingly, given the long duration of the spiral SASI episodes, the spectrograms of the expected neutrino event rate carry clear signs of the evolution of the blue spiral SASI frequency as a function of time, as the shock radius and post-shock fluid velocity evolve. Due to the high accretion luminosity and its large-amplitude SASI-induced modulations, any contribution from asymmetric (dipolar or quadrupolar) neutrino emission associated with the lepton emission self-sustained asymmetry (LESA) is far subdominant in the neutrino signal.

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.

The mixing of active neutrinos with their sterile counterparts with keV mass is known to have a potentially major impact on the energy loss from the supernova core. By relying on a set of three static hydrodynamical backgrounds mimicking the early accretion phase and the Kelvin-Helmoltz cooling phase of a supernova, we develop the first self-consistent, radial- and time-dependent treatment of sterile and tau neutrinos mixing in the dense stellar core. We follow the flavor evolution by including ordinary matter effects, collisional production of sterile neutrinos, as well as reconversions of sterile states into active ones. The dynamical feedback of the sterile neutrino production on the matter background leads to the development of a tau lepton asymmetry that grows in time until it reaches a value larger than 0.15. Our results hint towards significant implications for the supernova physics, and call for a self-consistent modeling of the sterile neutrino transport in the supernova core to constrain the mixing parameters of sterile neutrinos.

We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity.

We briefly review the dominant neutrino fluxes at Earth from different sources and present the {\em Grand Unified Neutrino Spectrum\/} ranging from meV to PeV energies. For each energy band and source, we discuss both theoretical expectations and experimental data. This compact review should be useful as a brief reference to those interested in neutrino astronomy, fundamental particle physics, dark-matter detection, high-energy astrophysics, geophysics, and other related topics.

Rotation of core-collapse supernovae (SNe) affects the neutrino emission characteristics. By comparing the neutrino properties of three three-dimensional SN simulations of a 15 M_sun progenitor (one non-rotating model and two models rotating at different velocities), we investigate how the neutrino emission varies with the flow dynamics in the SN core depending on the different degrees of rotation. The large-amplitude sinusoidal modulations due to the standing accretion-shock instability (SASI) are weaker in both the rotating models than in the non-rotating case. The SN progenitor rotation reduces the radial velocities and radial component of the kinetic energy associated with convection interior to the proto-neutron star. This effect seems to disfavor the growth of the hemispheric neutrino-emission asymmetries associated with the lepton-emission self-sustained asymmetry (LESA). An investigation of the multipole expansion of the neutrino luminosity and the electron neutrino lepton number flux shows a dominant quadrupolar mode in rotating SN models. Our findings highlight the power of using neutrinos as probes of SN hydrodynamics.

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.

The sources of cosmic rays with energies above 55 EeV are still mysterious. A guaranteed associated flux of ultra high energy neutrinos known as the cosmogenic neutrino flux will be measured by next generation radio facilities, such as the proposed Giant Radio Array for Neutrino Detection (GRAND). By using the orthogonal information provided by the cosmogenic neutrino flux, we here determine the prospects of GRAND to constrain the source redshift evolution and the chemical composition of the cosmic ray sources. If the redshift evolution is known, independently on GRAND's energy resolution, GRAND with 200,000 antennas will constrain the proton/iron fraction to the $\sim5-10\%$ level after one year of data taking; on the other hand, if hints on the average source composition are given, GRAND will measure the redshift evolution of the sources to a $\sim 10\%$ uncertainty. However, the foreseen configuration of GRAND alone will not be able to break the degeneracy between redshift evolution of the sources and their composition. Our findings underline the discriminating potential of next generation radio array detectors and motivate further efforts in this direction.

Measuring the rotation of core-collapse supernovae (SN) and of their progenitor stars is extremely challenging. Here it is demonstrated that neutrinos may potentially be employed as stellar gyroscopes, if phases of activity by the standing accretion-shock instability (SASI) affect the neutrino emission prior to the onset of the SN explosion. This is shown by comparing the neutrino emission properties of self-consistent, three-dimensional (3D) SN simulations of a 15 M_sun progenitor without rotation as well as slow and fast rotation compatible with observational constraints. The explosion of the fast rotating model gives rise to long-lasting, massive polar accretion downflows with stochastic time-variability, detectable e.g. by the IceCube Neutrino Observatory for any observer direction. While spectrograms of the neutrino event rate of non-rotating SNe feature a well-known sharp peak due to SASI for observers located in the proximity of the SASI plane, the corresponding spectrograms of rotating models show activity over a wide range of frequencies, most notably above 200 Hz for rapid rotation. In addition, the Fourier power spectra of the event rate for rotating models exhibit a SASI peak with lower power than in non-rotating models. The spectra for the rotating models also show secondary peaks at higher frequencies with greater relative heights compared to the main SASI peak than for non-rotating cases. These rotational imprints will be detectable for SNe at 10 kpc or closer.

The IceCube Neutrino Observatory detects high energy astrophysical neutrinos in two event topologies: tracks and cascades. Since the flavor composition of each event topology differs, tracks and cascades can be used to test the neutrino properties and the mechanisms behind the neutrino production in astrophysical sources. Assuming a conventional model for the neutrino production, the IceCube data sets related to the two channels are in $>3\sigma$ tension with each other. Invisible neutrino decay with lifetime $\tau/m=10^2$ s/eV solves this tension. Noticeably, it leads to an improvement over the standard non-decay scenario of more than $3\sigma$ while remaining consistent with all other multi-messenger observations. In addition, our invisible neutrino decay model predicts a reduction of $59\%$ in the number of observed $\nu_\tau$ events which is consistent with the current observational deficit.

The detection of the diffuse supernova neutrino background (DSNB) will preciously contribute to gauge the properties of the core-collapse supernova population. We estimate the DSNB event rate in the next-generation neutrino detectors, Hyper-Kamiokande enriched with Gadolinium, JUNO, and DUNE. The determination of the supernova unknowns through the DSNB will be heavily driven by Hyper-Kamiokande, given its higher expected event rate, and complemented by DUNE that will help in reducing the parameters uncertainties. Meanwhile, JUNO will be sensitive to the DSNB signal over the largest energy range. A joint statistical analysis of the expected rates in 20 years of data taking from the above detectors suggests that we will be sensitive to the local supernova rate at most at a 20-33% level. A non-zero fraction of supernovae forming black holes will be confirmed at a 90% CL, if the true value of that fraction is larger than 20%. On the other hand, the DSNB events show extremely poor statistical sensitivity to the nuclear equation of state and mass accretion rate of the progenitors forming black holes.

The increasing statistics of the high-energy neutrino flux observed by the IceCube Observatory points towards an excess of events above the atmospheric neutrino background in the 30--400 TeV energy range. Such an excess is compatible with the findings of the ANTARES Telescope and it would naturally imply the possibility that more than one source class contributes to the observed flux. Electromagnetically hidden sources have been invoked to interpret this excess of events at low energies. By adopting a unified model for the electromagnetically bright and choked gamma-ray bursts and taking into account particle acceleration at the internal and collimation shock radii, we discuss whether bright and choked bursts are viable candidates. Our findings suggest that, although producing a copious neutrino flux, choked and bright astrophysical jets cannot be the dominant sources of the excess of neutrino events. A fine tuning of the model parameters or distinct scenarios for choked jets should be invoked in order to explain the low-energy neutrino data of IceCube and ANTARES.

Long duration gamma-ray bursts (GRBs) have often been considered as the natural evolution of some core-collapse supernova (CCSN) progenitors. However, the fraction of CCSNe linked to astrophysical jets and their properties are still poorly constrained. While any successful astrophysical jet harbored in a CCSN should produce high energy neutrinos, photons may be able to successfully escape the stellar envelope only for a fraction of progenitors, possibly leading to the existence of high-luminosity, low-luminosity and not-electromagnetically bright ("choked") GRBs. By postulating a CCSN-GRB connection, we accurately model the jet physics within the internal-shock GRB model and assume scaling relations for the GRB parameters that depend on the Lorentz boost factor $\Gamma$. The IceCube high energy neutrino flux is then employed as an upper limit of the neutrino background from electromagnetically bright and choked GRBs to constrain the jet and the progenitor properties. The current IceCube data set is compatible with up to $1\%$ of all CCSNe harboring astrophysical jets. Interestingly, those jets are predominantly choked. Our findings suggest that neutrinos can be powerful probes of the burst physics and can provide major insights on the CCSN-GRB connection.

The remnant of neutron star mergers is dense in neutrinos. By employing inputs from one hydrodynamical simulation of a binary neutron star merger remnant with a black hole of $3\ M_\odot$ in the center, dimensionless spin parameter $0.8$ and an accretion torus of $0.3\ M_\odot$, the neutrino emission properties are investigated as the merger remnant evolves. Initially, the local number density of $\bar{\nu}_e$ is larger than that of $\nu_e$ everywhere above the remnant. Then, as the torus approaches self-regulated equilibrium, the local abundance of neutrinos overcomes that of antineutrinos in a funnel around the polar region. The region where the fast pairwise flavor conversions can occur shrinks accordingly as time evolves. Still, we find that fast flavor conversions do affect most of the neutrino-driven ejecta. Assuming that fast flavor conversions lead to flavor equilibration, a significant enhancement of nuclei with mass numbers $A>130$ is found as well as a change of the lanthanide mass fraction by more than a factor of a thousand. Our findings hint towards a potentially relevant role of neutrino flavor oscillations for the prediction of the kilonova (macronova) lightcurves and motivate further work in this direction.

The massive neutron star (NS) or black hole (BH) accretion disk resulting from NS--NS or NS--BH mergers is dense in neutrinos. We present the first study on the role of angular distributions in the neutrino flavor conversion above the remnant disk. In particular, we focus on "fast" pairwise conversions whose rate depends on the local angular intensity of the electron lepton number carried by neutrinos. Because of the emission geometry and the flux density of $\bar{\nu}_e$ being larger than that of $\nu_e$, fast conversions prove to be a generic phenomenon in NS--NS and NS--BH mergers for physically motivated disturbances in the mean field of flavor coherence. Our findings suggest that, differently from the core-collapse supernova case, fast flavor conversions seem to be unavoidable in compact mergers and could have major consequences for the jet dynamics and the synthesis of elements above the remnant disk.

We perform the first one-point fluctuation analysis of the high-energy neutrino sky. This method reveals itself to be especially suited to contemporary neutrino data, as it allows to study the properties of the astrophysical components of the high-energy flux detected by the IceCube telescope, even with low statistics and in the absence of point source detection. Besides the veto-passing atmospheric foregrounds, we adopt a simple model of the high-energy neutrino background by assuming two main extra-galactic components: star-forming galaxies and blazars. By leveraging multi-wavelength data from Herschel and Fermi, we predict the spectral and anisotropic probability distributions for their expected neutrino counts in IceCube. We find that star-forming galaxies are likely to remain a diffuse background due to the poor angular resolution of IceCube, and we determine an upper limit on the number of shower events that can reasonably be associated to blazars. We also find that upper limits on the contribution of blazars to the measured flux are unfavourably affected by the skewness of the blazar flux distribution. One-point event clustering and likelihood analyses of the IceCube HESE data suggest that this method has the potential to dramatically improve over more conventional model-based analyses, especially for the next generation of neutrino telescopes.

According to recent studies, the collective flavor evolution of neutrinos in core-collapse supernovae depends strongly on the flavor-dependent angular distribution of the local neutrino radiation field, notably on the angular intensity of the electron-lepton number carried by neutrinos. To facilitate further investigations of this subject, we study the energy and angle distributions of the neutrino radiation field computed with the Vertex neutrino-transport code for several spherically symmetric (1D) supernova simulations (of progenitor masses 11.2, 15 and 25 M_sun) and explain how to extract this information from additional models of the Garching group. Beginning in the decoupling region ("neutrino sphere"), the distributions are more and more forward peaked in the radial direction with an angular spread that is largest for $\nu_e$, smaller for $\bar\nu_e$, and smallest for $\nu_x$, where $x=\mu$ or $\tau$. While the energy-integrated $\nu_e$ minus $\bar\nu_e$ angle distribution has a dip in the forward direction, it does not turn negative in any of our investigated cases.

Constraints on the number and luminosity of the sources of the cosmic neutrinos detected by IceCube have been set by targeted searches for point sources. We set complementary constraints by using the 2MASS Redshift Survey (2MRS) catalogue, which maps the matter distribution of the local Universe. Assuming that the distribution of the neutrino sources follows that of matter we look for correlations between `warm' spots on the IceCube skymap and the 2MRS matter distribution. Through Monte Carlo simulations of the expected number of neutrino multiplets and careful modelling of the detector performance (including that of IceCube-Gen2) we demonstrate that sources with local density exceeding $10^{-6} \, \text{Mpc}^{-3}$ and neutrino luminosity $L_{\nu} \lesssim 10^{42} \, \text{erg} \, \text{s}^{-1}$ ($10^{41} \, \text{erg} \, \text{s}^{-1}$) will be efficiently revealed by our method using IceCube (IceCube-Gen2). At low luminosities such as will be probed by IceCube-Gen2, the sensitivity of this analysis is superior to requiring statistically significant direct observation of a point source.

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.

Dark matter detectors that utilize liquid xenon have now achieved tonne-scale targets, giving them sensitivity to all flavours of supernova neutrinos via coherent elastic neutrino-nucleus scattering. Considering for the first time a realistic detector model, we simulate the expected supernova neutrino signal for different progenitor masses and nuclear equations of state in existing and upcoming dual-phase liquid xenon experiments. We show that the proportional scintillation signal (S2) of a dual-phase detector allows for a clear observation of the neutrino signal and guarantees a particularly low energy threshold, while the backgrounds are rendered negligible during the supernova burst. XENON1T (XENONnT and LZ; DARWIN) experiments will be sensitive to a supernova burst up to 25 (35; 65) kpc from Earth at a significance of more than 5 sigma, observing approximately 35 (123; 704) events from a 27 Msun supernova progenitor at 10 kpc. Moreover, it will be possible to measure the average neutrino energy of all flavours, to constrain the total explosion energy, and to reconstruct the supernova neutrino light curve. Our results suggest that a large xenon detector such as DARWIN will be competitive with dedicated neutrino telescopes, while providing complementary information that is not otherwise accessible.

Neutrinos play a key role in core-collapse supernova explosions. Carrying information from deep inside the stellar core, neutrinos are direct probes of the supernova mechanism. Intriguing recent developments on the role of neutrinos in supernovae are reviewed, as well as our current understanding of the flavor conversions in the stellar envelope, and the detection perspectives of the next burst.

Long-duration gamma-ray bursts (GRBs) have been often considered as the natural evolution of some core-collapse supernovae (SNe). While GRBs with relativistic jets emit an electromagnetic signal, GRBs with mildly relativistic jets are opaque to photons and, therefore, could be detectable through neutrinos only. We discuss the possibility that successful GRBs and mildly relativistic jets belong to the same class of astrophysical transients with different Lorentz factor Gamma_b and study the production of high-energy neutrinos as a function of Gamma_b, by including both proton-photon and proton-proton interactions. By assuming a SN-GRB connection, we find that the diffuse neutrino emission from optically thick jets with Lorentz factors lower than the ones of successful GRBs can be one of the main components of the observed IceCube high-energy neutrino flux. Moreover, under the assumption that all these jets belong to the same class of astrophysical transients, we show that the IceCube high-energy neutrino data provide indirect constraints on the rate of non-successful jets, favoring a local rate lower than tens of percent of the local SN rate. These limits are currently comparable to the ones obtained in dedicated searches on choked sources and are expected to become tighter with accumulation of more high-energy neutrino data.

Neutrinos play a crucial role in the collapse and explosion of massive stars, governing the infall dynamics of the stellar core, triggering and fueling the explosion and driving the cooling and deleptonization of the newly formed neutron star. Due to their role neutrinos carry information from the heart of the explosion and, due to their weakly interacting nature, offer the only direct probe of the dynamics and thermodynamics at the center of a supernova. In this paper, we review the present status of modelling the neutrino physics and signal formation in collapsing and exploding stars. We assess the capability of current and planned large underground neutrino detectors to yield faithful information of the time and flavor dependent neutrino signal from a future Galactic supernova. We show how the observable neutrino burst would provide a benchmark for fundamental supernova physics with unprecedented richness of detail. Exploiting the treasure of the measured neutrino events requires a careful discrimination of source-generated properties from signal features that originate on the way to the detector. As for the latter, we discuss self-induced flavor conversions associated with neutrino-neutrino interactions that occur in the deepest stellar regions; matter effects that modify the pattern of flavor conversions in the dynamical stellar envelope; neutrino-oscillation signatures that result from structural features associated with the shock-wave propagation as well as turbulent mass motions in post-shock layers. Finally, we highlight our current understanding of the formation of the diffuse supernova neutrino background and we analyse the perspectives for a detection of this relic signal that integrates the contributions from all past core-collapse supernovae in the Universe.

Mounting evidence suggests that the TeV-PeV neutrino flux detected by the IceCube telescope has mainly an extragalactic origin. If such neutrinos are primarily produced by a single class of astrophysical sources via hadronuclear ($pp$) interactions, a similar flux of gamma-ray photons is expected. For the first time, we employ tomographic constraints to pinpoint the origin of the IceCube neutrino events by analyzing recent measurements of the cross correlation between the distribution of GeV gamma rays, detected by the Fermi satellite, and several galaxy catalogs in different redshift ranges. We find that the corresponding bounds on the neutrino luminosity density are up to one order of magnitude tighter than those obtained by using only the spectrum of the gamma-ray background, especially for sources with mild redshift evolution. In particular, our method excludes any hadronuclear source with a spectrum softer than $E^{-2.1}$ as a main component of the neutrino background, if its evolution is slower than $(1+z)^3$. Starburst galaxies, if able to accelerate and confine cosmic rays efficiently, satisfy both spectral and tomographic constraints.

Cosmic-ray (CR) protons can accumulate for cosmological times in clusters of galaxies. Their hadronic interactions with protons of the intra-cluster medium (ICM) generate secondary electrons, gamma-rays and high-energy neutrinos. In light of the high-energy neutrino events recently discovered by the IceCube observatory, we estimate the contribution from galaxy clusters to the diffuse gamma-ray and neutrino backgrounds. For the first time, we consistently take into account the synchrotron emission generated by secondary electrons and require the clusters radio counts to be respected. For a choice of parameters respecting current constraints from radio to gamma-rays, and assuming a proton spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the IceCube flux, and much less to the total extragalactic gamma-ray background observed by Fermi. They account for less than 1% for spectral indexes <-2. The high-energy neutrino flux observed by IceCube can be reproduced without violating radio constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the high-energy IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with the particle energy. We stress that our results are valid for all kind of sources injecting CR protons into the ICM, and that, while IceCube can test the most optimistic scenarios for spectral indexes >=-2.2 by stacking few nearby massive objects, clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds in any realistic scenario.

Gamma-ray bursts (GRBs) have been suggested as possible sources of the high-energy neutrino flux recently detected by the IceCube telescope. We revisit the fireball emission model and elaborate an analytical prescription to estimate the high-energy neutrino prompt emission from pion and kaon decays, assuming that the leading mechanism for the neutrino production is lepto-hadronic. To this purpose, we include hadronic, radiative and adiabatic cooling effects and discuss their relevance for long- (including high- and low-luminosity) and short-duration GRBs. The expected diffuse neutrino background is derived, by requiring that the GRB high-energy neutrino counterparts follow up-to-date gamma-ray luminosity functions and redshift evolutions of the long and short GRBs. Although dedicated stacking searches have been unsuccessful up to now, we find that GRBs could contribute up to a few % to the observed IceCube high-energy neutrino flux for sub-PeV energies, assuming that the latter has a diffuse origin. Gamma-ray bursts, especially low-luminosity ones, could however be the main sources of the IceCube high-energy neutrino flux in the PeV range. While high-luminosity and low-luminosity GRBs have comparable intensities, the contribution from the short-duration component is significantly smaller. Our findings confirm the most-recent IceCube results on the GRB searches and suggest that larger exposure is mandatory to detect high-energy neutrinos from high-luminosity GRBs in the near future.

Neutrino oscillations, especially to light sterile states, can affect the nucleosynthesis yields because of their possible feedback effect on the electron fraction (Ye). For the first time, we perform nucleosynthesis calculations for neutrino-driven wind trajectories from the neutrino-cooling phase of an 8.8 Msun electron-capture supernova, whose hydrodynamic evolution was computed in spherical symmetry with sophisticated neutrino transport and whose Ye evolution was post-processed by including neutrino oscillations both between active and active-sterile flavors. We also take into account the alpha-effect as well as weak magnetism and recoil corrections in the neutrino absorption and emission processes. We observe effects on the Ye evolution which depend in a subtle way on the relative radial positions of the sterile MSW resonances, of collective flavor transformations, and on the formation of alpha particles. For the adopted supernova progenitor, we find that neutrino oscillations, also to a sterile state with eV-mass, do not significantly affect the element formation and in particular cannot make the post-explosion wind outflow neutron rich enough to activate a strong r-process. Our conclusions become even more robust when, in order to mimic equation-of-state dependent corrections due to nucleon potential effects in the dense-medium neutrino opacities, six cases with reduced Ye in the wind are considered. In these cases, despite the conversion of active neutrinos to sterile neutrinos, Ye increases or is not significantly lowered compared to the values obtained without oscillations and active flavor transformations. This is a consequence of a complicated interplay between sterile-neutrino production, neutrino-neutrino interactions, and alpha-effect.

Cosmic-ray protons accumulate for cosmological times in clusters of galaxies as their typical radiative and diffusive escape times are longer than the Hubble time. Their hadronic interactions with protons of the intra-cluster medium generate secondary electrons, gamma-rays and neutrinos. We here estimate the contribution from clusters to the diffuse gamma-ray and neutrino backgrounds. We model the cluster population by means of their mass function, using a phenomenological luminosity-mass relation applied to all clusters, as well as a detailed semi-analytical model. Additionally, we consider observationally-motivated values for the cluster magnetic field. This is a crucial parameter since the observed radio counts due to synchrotron emission by secondary electrons need to be respected. For a choice of parameters respecting all current constraints, and assuming a spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the extragalactic gamma-ray background observed by Fermi and to the IceCube flux. They account for less than 1% for spectral indices <=-2. The neutrino flux observed by IceCube can be reproduced without violating constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with energy. In the case of proton-photon interactions, we find that very likely protons do not reach sufficiently high energies to produce neutrinos in clusters. We argue that our results are optimistic due to our assumptions, and that clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds. Finally, we find that the cluster contribution to the angular fluctuations in the gamma-ray background is subdominant, less than 10%. [abridged]

Star-forming galaxies have been predicted to contribute considerably to the diffuse gamma-ray background as they are guaranteed reservoirs of cosmic rays. Assuming that the hadronic interactions responsible for high-energy gamma rays also produce high-energy neutrinos and that O(100) PeV cosmic rays can be produced and confined in starburst galaxies, we here discuss the possibility that star-forming galaxies are also the main sources of the high-energy neutrinos observed by the IceCube experiment. First, we compute the diffuse gamma-ray background from star-forming galaxies, adopting the latest Herschel PEP/HerMES luminosity function and relying on the correlation between the gamma-ray and infrared luminosities reported by Fermi observations. Then we derive the expected intensity of the diffuse high-energy neutrinos from star-forming galaxies including normal and starburst galaxies. Our results indicate that starbursts, including those with active galactic nuclei and galaxy mergers, could be the main sources of the high-energy neutrinos observed by the IceCube experiment. We find that assuming a cosmic-ray spectral index of 2.1-2.2 for all starburst-like galaxies, our predictions can be consistent with both the Fermi and IceCube data, but larger indices readily fail to explain the observed diffuse neutrino flux. Taking the starburst high-energy spectral index as free parameter, and extrapolating from GeV to PeV energies, we find that the spectra harder than E^(-2.15) are likely to be excluded by the IceCube data, which can be more constraining than the Fermi data for this population.

The neutrino emission characteristics of the first full-scale three-dimensional supernova simulations with sophisticated three-flavor neutrino transport for three models with masses 11.2, 20 and 27 M_sun are evaluated in detail. All the studied progenitors show the expected hydrodynamical instabilities in the form of large-scale convective overturn. In addition, the recently identified LESA phenomenon (lepton-number emission self-sustained asymmetry) is generic for all our cases. Pronounced SASI (standing accretion-shock instability) activity appears in the 20 and 27 M_sun cases, partly in the form of a spiral mode, inducing large but direction and flavor-dependent modulations of neutrino emission. These modulations can be clearly identified in the existing IceCube and future Hyper-Kamiokande detectors, depending on distance and detector location relative to the main SASI sloshing direction.

During the stalled-shock phase of our 3D hydrodynamical core-collapse simulations with energy-dependent, 3-flavor neutrino transport, the lepton-number flux (nue minus antinue) emerges predominantly in one hemisphere. This novel, spherical-symmetry breaking neutrino-hydrodynamical instability is termed LESA for "Lepton-number Emission Self-sustained Asymmetry." While the individual nue and antinue fluxes show a pronounced dipole pattern, the heavy-flavor neutrino fluxes and the overall luminosity are almost spherically symmetric. Initially, LESA seems to develop stochastically from convective fluctuations, it exists for hundreds of milliseconds or more, and it persists during violent shock sloshing associated with the standing accretion shock instability. The nue minus antinue flux asymmetry originates mainly below the neutrinosphere in a region of pronounced proto-neutron star (PNS) convection, which is stronger in the hemisphere of enhanced lepton-number flux. On this side of the PNS, the mass-accretion rate of lepton-rich matter is larger, amplifying the lepton-emission asymmetry, because the spherical stellar infall deflects on a dipolar deformation of the stalled shock. The increased shock radius in the hemisphere of less mass accretion and minimal lepton-number flux (antinue flux maximum) is sustained by stronger convection on this side, which is boosted by stronger neutrino heating because the average antinue energy is higher than the average nue energy. Asymmetric heating thus supports the global deformation despite extremely nonstationary convective overturn behind the shock. While these different elements of LESA form a consistent picture, a full understanding remains elusive at present. There may be important implications for neutrino-flavor oscillations, the neutron-to-proton ratio in the neutrino-heated supernova ejecta, and neutron-star kicks, which remain to be explored.

The first full-scale three-dimensional (3D) core-collapse supernova (SN) simulations with sophisticated neutrino transport show pronounced effects of the standing accretion shock instability (SASI) for two high-mass progenitors (20 and 27 M_sun). In a low-mass progenitor (11.2 M_sun), large-scale convection is the dominant nonradial hydrodynamic instability in the postshock accretion layer. The SASI-associated modulation of the neutrino signal (80 Hz in our two examples) will be clearly detectable in IceCube or the future Hyper-Kamiokande detector, depending on progenitor properties, distance, and observer location relative to the main SASI sloshing direction. The neutrino signal from the next galactic SN can therefore diagnose the nature of the hydrodynamic instability.

To study the capabilities of supernova neutrino detectors, the instantaneous spectra are often represented by a quasi-thermal distribution of the form f(E) = E^alpha e^{-(alpha+1)E/E_{av}} where E_{av} is the average energy and alpha a numerical parameter. Based on a spherically symmetric supernova model with full Boltzmann neutrino transport we have, at a few representative post-bounce times, re-converged the models with vastly increased energy resolution to test the fit quality. For our examples, the spectra are well represented by such a fit in the sense that the counting rates for a broad range of target nuclei, sensitive to different parts of the spectrum, are reproduced very well. Therefore, the mean energy and root-mean-square energy of numerical spectra hold enough information to provide the correct alpha and to forecast the response of multi-channel supernova neutrino detection.

Neutrino oscillations in the Earth matter may introduce peculiar modulations in the supernova (SN) neutrino spectra. The detection of this effect has been proposed as diagnostic tool for the neutrino mass hierarchy at "large" 1-3 leptonic mixing angle theta13. We perform an updated study on the observability of this effect at large next-generation underground detectors (i.e., 0.4 Mton water Cherenkov, 50 kton scintillation and 100 kton liquid Argon detectors) based on neutrino fluxes from state-of-the-art SN simulations and accounting for statistical fluctuations via Montecarlo simulations. Since the average energies predicted by recent simulations are lower than previously expected and a tendency towards the equalization of the neutrino fluxes appears during the SN cooling phase, the detection of the Earth matter effect will be more challenging than expected from previous studies. We find that none of the proposed detectors shall be able to detect the Earth modulation for the neutrino signal of a typical galactic SN at 10 kpc. It should be observable in a 100 kton liquid Argon detector for a SN at few kpc and all three detectors would clearly see the Earth signature for very close-by stars only (d ~ 0.2 kpc). Finally, we show that adopting IceCube as co-detector together with a Mton water Cherenkov detector is not a viable option either.

Recent cosmological data favour additional relativistic degrees of freedom beyond the three active neutrinos and photons, often referred to as 'dark' radiation. Light sterile neutrinos is one of the prime candidates for such additional radiation. However, constraints on sterile neutrinos based on the current cosmological data have been derived using simplified assumptions about thermalisation of the sterile neutrino at the Big Bang Nucleosynthesis (BBN) epoch. These assumptions are not necessarily justified and here we solve the full quantum kinetic equations in the (1 active + 1 sterile) scenario and derive the number of thermalised species just before BBN begins (T~1MeV) for null (L=0) and large (L=0.01) initial lepton asymmetry and for a range of possible mass-mixing parameters. We find that the full thermalisation assumption during the BBN epoch is justified for initial small lepton asymmetry only. Partial or null thermalisation occurs when the initial lepton asymmetry is large.

We estimate the diffuse supernova neutrino background (DSNB) using the recent progenitor-dependent, long-term supernova simulations from the Basel group and including neutrino oscillations at several post-bounce times. Assuming multi-angle matter suppression of collective effects during the accretion phase, we find that oscillation effects are dominated by the matter-driven MSW resonances, while neutrino-neutrino collective effects contribute at the 5-10% level. The impact of the neutrino mass hierarchy, of the time-dependent neutrino spectra and of the diverse progenitor star population is 10% or less, small compared to the uncertainty of at least 25% of the normalization of the supernova rate. Therefore, assuming that the sign of the neutrino mass hierarchy will be determined within the next decade, the future detection of the DSNB will deliver approximate information on the MSW-oscillated neutrino spectra. With a reliable model for neutrino emission, its detection will be a powerful instrument to provide complementary information on the star formation rate and for learning about stellar physics.

Neutrinos streaming from a supernova (SN) core occasionally scatter in the envelope, producing a small "neutrino halo" with a much broader angle distribution than the primary flux originating directly from the core. Cherry et al. (2012) have recently pointed out that, during the accretion phase, the halo actually dominates neutrino-neutrino refraction at distances exceeding some 100 km. However, the multiangle matter effect (which increases if the angle distribution is broader) still appears to suppress self-induced flavor conversion during the accretion phase.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

Motivated by recent hints for sterile neutrinos from the reactor anomaly, we study active-sterile conversions in a three-flavor scenario (2 active + 1 sterile families) for three different representative times during the neutrino-cooling evolution of the proto-neutron star born in an electron-capture supernova. In our "early model" (0.5 s post bounce), the nu_e-nu_s MSW effect driven by Delta m^2=2.35 eV^2 is dominated by ordinary matter and leads to a complete nu_e-nu_s swap with little or no trace of collective flavor oscillations. In our "intermediate" (2.9 s p.b.) and "late models" (6.5 s p.b.), neutrinos themselves significantly modify the nu_e-nu_s matter effect, and, in particular in the late model, nu-nu refraction strongly reduces the matter effect, largely suppressing the overall nu_e-nu_s MSW conversion. This phenomenon has not been reported in previous studies of active-sterile supernova neutrino oscillations. We always include the feedback effect on the electron fraction Y_e due to neutrino oscillations. In all examples, Y_e is reduced and therefore the presence of sterile neutrinos can affect the conditions for heavy-element formation in the supernova ejecta, even if probably not enabling the r-process in the investigated outflows of an electron-capture supernova. The impact of neutrino-neutrino refraction is strong but complicated, leaving open the possibility that with a more complete treatment, or for other supernova models, active-sterile neutrino oscillations could generate conditions suitable for the r-process.

We study collective oscillations of a two-flavor neutrino system with arbitrary but fixed density. In the vacuum limit, modes with different energies quickly de-phase (kinematical decoherence), whereas in the limit of infinite density they lock to each other (synchronization). For intermediate densities, we find different classes of solutions. There is always a phase transition in the sense of partial synchronization occurring only above a density threshold. For small mixing angles, partial or complete decoherence can be induced by a parametric resonance, introducing a new time scale to the problem, the final outcome depending on the spectrum and mixing angle. We derive an analytic relation that allows us to calculate the late-time degree of coherence based on the spectrum alone.

Precision cosmology and big-bang nucleosynthesis mildly favor extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to derive credible regions for the assumed common mass scale m_s and effective number N_s of thermally excited sterile neutrino states. Our results are compatible with the existence of one or perhaps two sterile neutrinos, as suggested by LSND and MiniBooNE, if m_s is in the sub-eV range.

It was recently realized that three-flavor effects could peculiarly modify the development of spectral splits induced by collective oscillations, for supernova neutrinos emitted during the cooling phase of a protoneutron star. We systematically explore this case, explaining how the impact of these three-flavor effects depends on the ordering of the neutrino masses. In inverted mass hierarchy, the solar mass splitting gives rise to instabilities in regions of the (anti)neutrino energy spectra that were otherwise stable under the leading two-flavor evolution governed by the atmospheric mass splitting and by the 1-3 mixing angle. As a consequence, the high-energy spectral splits found in the electron (anti)neutrino spectra disappear, and are transferred to other flavors. Imperfect adiabaticity leads to smearing of spectral swap features. In normal mass hierarchy, the three-flavor and the two-flavor instabilities act in the same region of the neutrino energy spectrum, leading to only minor departures from the two-flavor treatment.

Collective flavor transformations in supernovae, caused by neutrino-neutrino interactions, are essentially a two-flavor phenomenon driven by the atmospheric mass difference and the small mixing angle theta_13. In the two-flavor approximation, the initial evolution depends logarithmically on theta_13 and the system remains trapped in an unstable fixed point for theta_13 = 0. However, any effect breaking exact nu_mu-nu_tau equivalence triggers the conversion. Such three-flavor perturbations include radiative corrections to weak interactions, small differences between the nu_mu and nu_tau fluxes, or non-standard interactions. Therefore, extremely small values of theta_13 are in practice equivalent, the fate of the system depending only on the neutrino spectra and their mass ordering.

In core-collapse supernovae, the nu_e and anti-nu_e species may experience collective flavor swaps to non-electron species nu_x, within energy intervals limited by relatively sharp boundaries ("splits"). These phenomena appear to depend sensitively upon the initial energy spectra and luminosities. We investigate the effect of generic variations of the fractional luminosities (l_e, l_{anti-e}, l_x) with respect to the usual "energy equipartition" case (1/6, 1/6, 1/6), within an early-time supernova scenario with fixed thermal spectra and total luminosity. We represent the constraint l_e+l_{anti-e}+4l_x=1 in a ternary diagram, which is explored via numerical experiments (in single-angle approximation) over an evenly-spaced grid of points. In inverted hierarchy, single splits arise in most cases, but an abrupt transition to double splits is observed for a few points surrounding the equipartition one. In normal hierarchy, collective effects turn out to be unobservable at all grid points but one, where single splits occur. Admissible deviations from equipartition may thus induce dramatic changes in the shape of supernova (anti)neutrino spectra. The observed patterns are interpreted in terms of initial flavor polarization vectors (defining boundaries for the single/double split transitions), lepton number conservation, and minimization of potential energy.

Neutrino and antineutrino fluxes from a core-collapse galactic supernova are studied, within a representative three-flavor scenario with inverted mass hierarchy and tiny 1-3 mixing. The initial flavor evolution is dominated by collective self-interaction effects, which are computed in a full three-family framework along an averaged radial trajectory. During the whole time span considered (t=1-20 s), neutrino and antineutrino spectral splits emerge as dominant features in the energy domain for the final, observable fluxes. Some minor or unobservable three-family features (e.g., related to the muonic-tauonic flavor sector) are also discussed for completeness. The main results can be useful for SN event rate simulations in specific detectors.

In the dense supernova core, self-interactions may align the flavor polarization vectors of neutrinos and antineutrinos, and induce collective flavor transformations. Different alignment ansatzes are known to describe approximately the phenomena of synchronized or bipolar oscillations, and the split of neutrino energy spectra. We discuss another phenomenon observed in some numerical experiments in inverted hierarchy, showing features akin to a low-energy split of antineutrino spectra. The phenomenon appears to be approximately described by another alignment ansatz which, in the considered scenario, reduces the (nonadiabatic) dynamics of all energy modes to only two neutrino plus two antineutrino modes. The associated spectral features, however, appear to be fragile when passing from single- to multi-angle simulations.