We introduce a novel methodology for establishing the presence of Standing Accretion Shock Instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda $\it {et\, al.}$ (2017) as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to $\sim 6$ kpc for IceCube and and up to $\sim 3$ kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

We introduce a phenomenological concordance scenario with a relativistic jet for the Tidal Disruption Event (TDE) AT2019dsg, which has been proposed as source of the astrophysical neutrino event IceCube-191001A. Noting that AT2019dsg is one of the brightest (and few) TDEs observed in X-rays, we connect the neutrino production with the X-rays: an expanding cocoon causes the progressive obscuration of the X-rays emitted by the accretion disk, while at the same time it provides a sufficiently intense external target of back-scattered X-rays for photo-pion production off protons. We also describe the late-term emission of the neutrino (about 150 days after the peak), by scaling the production radius with the black body radius. Our energetics and assumptions for the jet and the cocoon are compatible with expectations from numerical simulations of TDEs. We predict 0.26 neutrino events in the right energy range in IceCube.

We explore the potential of current and future liquid scintillator neutrino detectors of O(10) kt mass to localize a presupernova neutrino signal in the sky. In the hours preceding the core collapse of a nearby star (at distance D < 1 kpc), tens to hundreds of inverse beta decay events will be recorded, and their reconstructed topology in the detector can be used to estimate the direction to the star. Although the directionality of inverse beta decay is weak (~8% forward-backward asymmetry for currently available liquid scintillators), we find that for a fiducial signal of 200 events (which is realistic for Betelgeuse), a positional error of ~60 degrees can be achieved, resulting in the possibility to narrow the list of potential stellar candidates to less than ten, typically. For a configuration with improved forward-backward asymmetry (~40%, as expected for a lithium-loaded liquid scintillator), the angular sensitivity improves to ~15 degrees, and - when a distance upper limit is obtained from the overall event rate - it is in principle possible to uniquely identify the progenitor star. Any localization information accompanying an early supernova alert will be useful to multi-messenger observations and to particle physics tests using collapsing stars.

We discuss the potential of detecting thermal neutrinos from matter-rich binary mergers, via a decades-long multi-messenger campaign involving a Mt-scale water Cherenkov neutrino detector and one or more next generation gravitational wave detectors, capable of observing mergers up to redshift $z\sim 2$. The search of neutrinos in time-coincidence with gravitational wave detections will allow to identify single neutrinos from individual mergers above the background, and to study their distributions in energy, redshift and type (double neutron-star or neutron-star-black hole merger) of the candidate sources. We find that, for merger rates consistent with current LIGO-Virgo constraints, and for a $100~{\rm Mt\cdot yr}$ exposure, between ${\mathcal O(10^{-1})}$ and ${\mathcal O(10)}$ neutrino events are expected. For extreme cases of mergers with more than $10^{52}$ ergs emitted in $\mathrel{{\bar \nu}_e}$, the number of events can be as large as $\sim 100$, with sensitivity to mergers up to redshift $z\sim 0.5$ or so. Such scenarios can already be tested with a $10~{\rm Mt\cdot yr}$ exposure, resulting in constraints on the post-merger evolution of the systems being considered.

A recent analysis of cosmic ray air showers observed at the Pierre Auger Observatory indicates that nearby starburst galaxies (SBGs) might be the cause of ~10% of the Ultra-High-Energy Cosmic Ray flux at energies E > 39 EeV. Since high energy neutrinos are a direct product of cosmic ray interactions, we investigate SBGs as a possible source of some of the ~ 10^-2 - 1 PeV neutrinos observed at IceCube. A statistical analysis is performed to establish the degree of positional correlation between the observed neutrinos and a set of 45 nearby radio- and infrared-bright SBGs. Our results are consistent with no causal correlation. However, a scenario where ~10% of the High Energy Starting Events (HESE) in the detector are coming from the candidate SBGs is not excluded. The same conclusion is reached for different data subsets, as well as two different subsets of SBGs motivated by the Pierre Auger Observatory analysis.

We study Primordial Black Holes (PBHs) as sources of massive neutrinos via Hawking radiation. Emphasizing that black holes emit neutrino mass eigenstates, we describe quantitatively how the PBH evolution and lifetime is affected by the mass and fermionic -- Dirac or Majorana -- nature of neutrinos. In the case of Dirac neutrinos, PBHs radiate right-handed and left-handed neutrinos in equal amounts, thus possibly increasing the effective number of neutrino species, $N_{\rm eff}$. Considering the full system of Friedmann equations, and the current constraint on $N_{\rm eff}$, we derive a novel bound on the initial PBH fraction $\beta^\prime$ for black hole masses $4\times 10^7\ {\rm g}\lesssim M_i \lesssim 10^9$ g. Future measurements of $N_{\rm eff}$ may be able to constraint the initial fraction for black hole masses as low as 1 g. If an excess in $N_{\rm eff}$ is found, PBHs with Dirac neutrinos could provide a minimal explanation of it. For example, for $10^7\ {\rm g} \lesssim M_i\lesssim 10^9$ g and $\beta^\prime \gtrsim 10^{-13}$, an excess radiation at the level of $0.2\lesssim \Delta N_{\rm eff}\lesssim 0.37$ is produced, which can alleviate the tension of the Hubble parameter measurements. Finally, we obtain the diffuse flux of right-helical neutrinos from PBHs at the Earth, and show that their detection in a PTOLEMY-like detector (using neutrino capture on Tritium) would be difficult.

Tidal Disruption Events (TDEs) are processes where stars are torn apart by the strong gravitational force near to a massive or supermassive black hole. If a jet is launched in such a process, particle acceleration may take place in internal shocks. We demonstrate that jetted TDEs can simultaneously describe the observed neutrino and cosmic ray fluxes at the highest energies if stars with heavier compositions, such as carbon-oxygen white dwarfs, are tidally disrupted and these events are sufficiently abundant. We simulate the photo-hadronic interactions both in the TDE jet and in the propagation through the extragalactic space and we show that the simultaneous description of Ultra-High Energy Cosmic Ray (UHECR) and PeV neutrino data implies that a nuclear cascade in the jet develops by photo-hadronic interactions.

We estimate the Diffuse Supernova Neutrino Background (DSNB) using the simulation results for neutron star-forming and black hole-forming stellar collapses from the Garching group. Scenarios with different distributions of black-hole forming collapses with the progenitor mass are discussed, and the uncertainty on the cosmological rate of collapses is included. The electron antineutrino component of the DSNB above 11 MeV of energy is found to be ~(1.4 - 3.7) cm^-2 s^-1; the contribution of black hole-forming collapses could dominate the flux above ~25 MeV. We calculate the potential of detecting the DSNB at SuperK-Gd and JUNO, in about a decade-long period of operation. We find that, in our model, it is likely that a significant excess above the background will be obtained at JUNO, while detection will be more difficult at SuperK-Gd. The potential when the two experimental results are examined jointly is discussed as well. We also consider an example of a future O(10) kt slow liquid scintillator detector, and show that there the odds of detection are very good. Our results motivate experimental efforts in reducing the backgrounds due to neutral current scattering of atmospheric neutrinos in SuperK-Gd.

We present a new calculation of the neutrino flux received at Earth from a massive star in the $\sim 24$ hours of evolution prior to its explosion as a supernova (presupernova). Using the stellar evolution code MESA, the neutrino emissivity in each flavor is calculated at many radial zones and time steps. In addition to thermal processes, neutrino production via beta processes is modeled in detail, using a network of 204 isotopes. We find that the total produced $\nu_{e}$ flux has a high energy spectrum tail, at $E \gtrsim 3 - 4$ MeV, which is mostly due to decay and electron capture on isotopes with $A = 50 - 60$. In a tentative window of observability of $E \gtrsim 0.5$ MeV and $t < 2$ hours pre-collapse, the contribution of beta processes to the $\nu_{e}$ flux is at the level of $\sim90\%$ . For a star at $D=1$ kpc distance, a 17 kt liquid scintillator detector would typically observe several tens of events from a presupernova, of which up to $\sim 30\%$ due to beta processes. These processes dominate the signal at a liquid argon detector, thus greatly enhancing its sensitivity to a presupernova.

We study the production of high energy neutrinos in jets from the tidal disruption of stars by supermassive black holes. The diffuse neutrino flux expected from these tidal disruption events (TDEs) is calculated both analytically and numerically, taking account the dependence of the rate of TDEs on the redshift and black hole mass. We find that ~10% of the observed diffuse flux at IceCube at an energy of about 1 PeV can come from TDEs if the characteristics of known jetted tidal disruption events are assumed to apply to the whole population of these sources. If, however, plausible scalings of the jet Lorentz factor or variability timescale with the black hole mass are taken into account, the contribution of the lowest mass black holes to the neutrino flux is enhanced. In this case, TDEs can account for most of the neutrino flux detected at IceCube, describing both the neutrino flux normalization and spectral shape with moderate baryonic loadings. While the uncertainties on our assumptions are large, a possible signature of TDEs as the origin of the IceCube signal is the transition of the flux flavor composition from a pion beam to a muon damped source at the highest energies, which will also result in a suppression of Glashow resonance events.

We present a new calculation of neutrino emissivities and energy spectra from a massive star going through the advanced stages of nuclear burning (presupernova) in the months before becoming a supernova. The contributions from beta decay and electron capture, pair annihilation, plasmon decay, and the photoneutrino process are modeled in detail, using updated tabulated nuclear rates. We also use realistic conditions of temperature, density, electron fraction and nuclear isotopic composition of the star from the state of the art stellar evolution code MESA. Results are presented for a set of progenitor stars with mass between 15 $M_\odot$ and 30 $M_\odot$. It is found that beta processes contribute substantially to the neutrino emissivity above realistic detection thresholds of few MeV, at selected positions and times in the evolution of the star.

I review the physics of the Diffuse Supernova Neutrino flux (or Background, DSNB), in the context of future searches at the next generation of neutrino observatories. The theory of the DSNB is discussed in its fundamental elements, namely the cosmological rate of supernovae, neutrino production inside a core collapse supernova, redshift, and flavor oscillation effects. The current upper limits are also reviewed, and results are shown for the rates and energy distributions of the events expected at future liquid argon and liquid scintillator detectors of O(10) kt mass, and water Cherenkov detectors up to a 0.5 Mt mass. Perspectives are given on the significance of future observations of the DSNB, both at the discovery and precision phases, for the investigation of the physics of supernovae and of the properties of the neutrino.

The IceCube Neutrino Observatory has provided the first map of the high energy (~ 0.01 -- 1 PeV) sky in neutrinos. Since neutrinos propagate undeflected, their arrival direction is an important identifier for sources of high energy particle acceleration. Reconstructed arrival directions are consistent with an extragalactic origin, with possibly a galactic component, of the neutrino flux. We present a statistical analysis of positional coincidences of the IceCube neutrinos with known astrophysical objects from several catalogs. When considering starburst galaxies with the highest flux in gamma-rays and infrared radiation, up to n=8 coincidences are found, representing an excess over the ~4 predicted for the randomized, or "null" distribution. The probability that this excess is realized in the null case, the p-value, is p=0.042. This value falls to p=0.003 for a partial subset of gamma-ray-detected starburst galaxies and superbubble regions in the galactic neighborhood. Therefore, it is possible that starburst galaxies, and the typically hundreds of superbubble regions within them, might account for a portion of IceCube neutrinos. The physical plausibility of such correlation is discussed briefly.

The Fermi Bubbles have been imaged in sub-TeV gamma rays at Fermi-LAT, and, if their origin is hadronic, they might have been seen with low statistics in $\sim 0.1- 1$ PeV neutrinos at IceCube. We discuss the detectability of these objects at the new High Altitude Water Cherenkov (HAWC) gamma ray detector. HAWC will view the North Bubble for $\sim 2-3$ hours a day, and will map its spectrum at 0.1-100 TeV. For the hard primary proton spectrum required to explain five events at IceCube, a high significance detection at HAWC will be achieved in less than 30 days. The combination of results at HAWC and IceCube will substantiate the hadronic model, or constrain its spectral parameters.

The IceCube Neutrino Observatory announced thirty-seven candidate events observed with deposited energies above $\sim$ 30 TeV with three-year dataset, more than expected from atmospheric backgrounds. We discuss the detectability of the Fermi Bubbles (FB) by IceCube and show that up to 6 -- 7 of the 37 events could originate from the FB. If the observed gamma rays from the FB are created due to baryonic mechanism, high-energy ($>$ GeV) neutrinos should be emitted as a counterpart. These neutrinos should be detectable as shower- or track-like events at a Km$^3$ neutrino detector. For a hard primary cosmic-ray proton spectrum, $E^{-2.1}$, and cutoff energy at or above 10 PeV, the FB flux substantially exceeds the atmospheric backgrounds. For a steeper spectrum, $E^{-2.3}$, and/or lower cutoff energy, detection with high significance will require a longer running time.

We study the physics potential of the detection of the Cosmic Neutrino Background via neutrino capture on tritium, taking the proposed PTOLEMY experiment as a case study. With the projected energy resolution of $\Delta \sim$ 0.15 eV, the experiment will be sensitive to neutrino masses with degenerate spectrum, $m_1 \simeq m_2 \simeq m_3 = m_\nu \gtrsim 0.1$ eV. These neutrinos are non-relativistic today; detecting them would be a unique opportunity to probe this unexplored kinematical regime. The signature of neutrino capture is a peak in the electron spectrum that is displaced by $2 m_{\nu}$ above the beta decay endpoint. The signal would exceed the background from beta decay if the energy resolution is $\Delta \lesssim 0.7 m_\nu $. Interestingly, the total capture rate depends on the origin of the neutrino mass, being $\Gamma^{\rm D} \simeq 4$ and $\Gamma^{\rm M} \simeq 8$ events per year (for a 100 g tritium target) for unclustered Dirac and Majorana neutrinos, respectively. An enhancement of the rate of up to $\mathcal{O}(1)$ is expected due to gravitational clustering, with the unique potential to probe the local overdensity of neutrinos. Turning to more exotic neutrino physics, PTOLEMY could be sensitive to a lepton asymmetry, and reveal the eV-scale sterile neutrino that is favored by short baseline oscillation searches. The experiment would also be sensitive to a neutrino lifetime on the order of the age of the universe and break the degeneracy between neutrino mass and lifetime which affects existing bounds.

We discuss the possibility that the IceCube neutrino telescope might be observing the Fermi Bubbles. If the bubbles discovered in gamma rays originate from accelerated protons, they should be strong emitters of high energy (> GeV) neutrinos. These neutrinos are detectable as shower- or track-like events at a Km^3 neutrino observatory. For a primary cosmic ray cutoff energy at or above 10 PeV, the Fermi Bubble flux substantially exceeds the atmospheric background, and could account for up to ~ 4-5 of the 28 events detected above ~ 30 TeV at IceCube. Running the detector for ~ 5-7 more years should be sufficient to discover this flux at high significance. For a primary cosmic ray flux with steeper spectrum, and/or lower cutoff energy, longer running times will be required to overcome the background.

We study absorption of ultra high energy neutrinos by the cosmic neutrino background, with full inclusion of the effect of the thermal distribution of the background on the resonant annihilation channel. For a hierarchical neutrino mass spectrum (with at least one neutrino with mass below $\sim 10^{-2}$ eV), thermal effects are important for ultra high energy neutrino sources at $z \gtrsim 16$. The neutrino transmission probability shows no more than two separate suppression dips since the two lightest mass eigenstates contribute as a single species when thermal effects are included. Results are applied to a number of models of ultra high energy neutrino emission. Suppression effects are strong for sources that extend beyond $z \sim 10$, which can be realized for certain top down scenarios, such as superheavy dark matter decays, cosmic strings and cosmic necklaces. For these, a broad suppression valley should affect the neutrino spectrum at least in the energy interval $10^{12} - 10^{13}$ GeV -- which therefore is disfavored for ultra high energy neutrino searches -- with only a mild dependence on the neutrino mass spectrum and hierarchy. The observation of absorption effects would indicate a population of sources beyond $z \sim 10$, and favor top-down mechanisms; it would also be an interesting probe of the physics of the relic neutrino background in the unexplored redshift interval $z \sim 10 -100$.

Tantalizing cosmological and terrestrial evidence suggests the number of light neutrinos may be greater than 3, motivating a careful re-examination of cosmological bounds on extra light species. Big Bang Nucleosynthesis constrains the number of relativistic neutrino species present during nucleosynthesis, $N_{eff}^{BBN}$, while measurements of the CMB angular power spectrum constrain the effective energy density in relativistic neutrinos at the time of matter-radiation equality, $N_{eff}^{CMB}$. There are a number of scenarios where new sterile neutrino species may have different contributions to $\Delta N_{eff}^{BBN}$ and $\Delta N_{eff}^{CMB}$, for masses that may be relevant to reconciling cosmological constraints with various terrestrial claims of neutrino oscillations. We consider a scenario with two sterile neutrinos and explore whether partial thermalization of the sterile states can ease the tension between cosmological constraints on $N_{eff}^{BBN}$ and terrestrial data. We then investigate the effect of a non-zero neutrino mass on their contribution to the radiation abundance, finding reductions in $\Delta N_{eff}^{CMB}$ of more than 5% for neutrinos with masses above 0.5 eV. While the effects we investigate here could play a role, we nevertheless find that two additional light sterile neutrinos species cannot fit all the data at the 95% confidence level.

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.

We study massive particle radiation from cosmic string kinks, and its observability in extremely high energy neutrinos. In particular, we consider the emission of moduli --- weakly coupled scalar particles predicted in supersymmetric theories --- from the kinks of cosmic string loops. Since kinks move at the speed of light on strings, moduli are emitted with large Lorentz factors, and eventually decay into many pions and neutrinos via hadronic cascades. The produced neutrino flux has energy $E \gtrsim 10^{11} \rm{GeV}$, and is affected by oscillations and absorption (resonant and non-resonant). It is observable at upcoming neutrino telescopes such as JEM-EUSO, and the radio telescopes LOFAR and SKA, for a range of values of the string tension, and of the mass and coupling constant of the moduli.

Recently the Fermi-LAT data have revealed two gamma-ray emitting bubble-shaped structures at the Galactic center. If the observed gamma rays have hadronic origin (collisions of accelerated protons), the bubbles must emit high energy neutrinos as well. This new, Galactic, neutrino flux should trace the gamma ray emission in spectrum and spatial extent. Its highest energy part, above 20-50 TeV, is observable at a kilometer scale detector in the northern hemisphere, such as the planned KM3NeT, while interesting constraints on it could be obtained by the IceCube Neutrino Observatory at the South pole. The detection or exclusion of neutrinos from the Fermi bubbles will discriminate between hadronic and leptonic models, thus bringing unique information on the still mysterious origin of these objects and on the time scale of their formation.

We discuss the diffuse flux of electron neutrinos and antineutrinos from cosmological failed supernovae, stars that collapse directly into a black hole, with no explosion. This flux has a hotter energy spectrum compared to regular, neutron-star forming collapses, and therefore it dominates the total diffuse flux from core collapses above 20-45 MeV of neutrino energy. Reflecting the features of the originally emitted neutrinos, the flux of nu_e and anti-nu_e at Earth is larger for larger survival probability of these species, and for stiffer equations of state of nuclear matter. In the energy window 19-29 MeV, the flux from failed supernovae is susbtantial, ranging from 7% to a dominant fraction of the total flux from all core collapses. It can be as large as phi = 0.38 s^{-1} cm^{-2} for anti-nu_e (phi = 0.28 s^{-1} cm^{-2} for nue), normalized to a local rate of core collapses of R_{cc}(0)=10^{-4} yr^{-1} Mpc^{-3}. In 5 years, a 0.45 Mt water Cherenkov detector should see 5-65 events from failed supernovae, while up to 160 events are expected for the same mass with Gadolinium addition. A 0.1 Mt liquid argon experiment should record 1-11 events. Signatures of neutrinos from failed supernovae are the enhancement of the total rates of events from core collapses (up to a factor of 2) and the appearance of high energy tails in the event spectra.

We study the detectability of neutrino bursts from nearby direct black hole-forming collapses (failed supernovae) at Megaton detectors. Due to their high energetics, these bursts could be identified - by the time coincidence of N >= 2 or N >= 3 events within a ~ 1 s time window - from as far as ~ 4-5 Mpc away. This distance encloses several supernova-rich galaxies, so that failed supernova bursts could be detected at a rate of up to one per decade, comparable to the expected rate of the more common, but less energetic, neutron star-forming collapses. Thus, the detection of a failed supernova within the lifetime of a Mt detector is realistic. It might give the first evidence of direct black hole formation, with important implications on the physics of this phenomenon.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

New data from WMAP have appeared, related to both the fractional energy density in relativistic species at decoupling and also the primordial helium abundance, at the same time as other independent observational estimates suggest a higher value of the latter than previously estimated. All the data are consistent with the possibility that the effective number of relativistic species in the radiation gas at the time of Big Bang Nucleosynthesis may exceed the value of 3, as expected from a CP-symmetric population of the known neutrino species. Here we explore the possibility that new neutrino physics accounts for such an excess. We explore different realizations, including neutrino asymmetry and new neutrino species, as well as their combination, and describe how existing constraints on neutrino physics would need to be relaxed as a result of the new data, as well as possible experimental tests of these possibilities.

In the delayed explosion scenario of core-collapse supernovae (SNe), the accretion phase shows pronounced convective overturns and a low-multipole hydrodynamic instability, the standing accretion shock instability (SASI). These effects imprint detectable fast time variations on the emerging neutrino flux. Among existing detectors, IceCube is best suited to this task, providing an event rate of ~1000 events per ms during the accretion phase for a fiducial SN distance of 10 kpc, comparable to what could be achieved with a megaton water Cherenkov detector. If the SASI activity lasts for several hundred ms, a Fourier component with an amplitude of 1% of the average signal clearly sticks out from the shot noise. We analyze in detail the output of axially symmetric hydrodynamical simulations that predict much larger amplitudes up to frequencies of a few hundred Hz. If these models are roughly representative for realistic SNe, fast time variations of the neutrino signal are easily detectable in IceCube or future megaton-class instruments. We also discuss the information that could be deduced from such a measurement about the physics in the SN core and the explosion mechanism of the SN.

I study the diffuse flux of electron antineutrinos from stellar collapses with direct black hole formation (failed supernovae). This flux is more energetic than that from successful supernovae, and therefore it might contribute substantially to the total diffuse flux above realistic detection thresholds. The total flux might be considerably higher than previously thought, and approach the sensitivity of SuperKamiokande. For more conservative values of the parameters, the flux from failed supernovae dominates for antineutrino energies above 30-45 MeV, with potential to give an observable spectral distortion at Megaton detectors.

We explore the possibility to observe the effects of electron neutrinos from past galactic supernovae, through a geochemical measurement of the amount of Technetium 97 produced by neutrino-induced reactions in a Molybdenum ore. The calculations we present take into account the recent advances in our knowledge of neutrino interactions, of neutrino oscillations inside a supernova, of the solar neutrino flux at Earth and of possible failed supernovae. The predicted Technetium 97 abundance is of the order of 10^7 atoms per 10 kilotons of ore, which is close to the current geochemical experimental sensitivity. Of this, 10-20% is from supernovae. Considering the comparable size of uncertainties, more precision in the modeling of neutrino fluxes as well as of neutrino cross sections is required for a meaningful measurement.

We analyze the 1496 days of SuperKamiokande data to put limits on the nue, anti-nue, numu + nutau and anti-numu + anti-nutau components of the diffuse flux of supernova neutrinos, in different energy intervals and for different neutrino energy spectra. By considering the presence of only one component at a time, we find the following bounds at 90% C.L. and for neutrino energy E>19.3$ MeV: Phi_{nue}<73.3-154 cm^-2 s^-1, Phi_{anti-nue} <1.4-1.9 cm^-2 s^-1, Phi_{numu+nutau} <(1.0-1.4) 10^3 cm^-2 s^-1, and Phi_{anti-numu+ anti-nutau} <(1.3-1.8) 10^3 cm^-2 s^-1, where the intervals account for varying the neutrino spectrum. In the interval E = 22.9 - 36.9 MeV, we find Phi_{nue}<39-54 cm^-2 s^-1, which improves on the existing limit from SNO in the same energy window. Our results for numu + nutau and anti-numu + anti-nutau improve by about four orders of magnitude over the previous best constraints from LSD.

We discuss the flavor conversion of neutrinos from core collapse supernovae that have oxygen-neon-magnesium (ONeMg) cores. Using the numerically calculated evolution of the star up to 650 ms post bounce, we find that, for the normal mass hierarchy, the electron neutrino flux in a detector shows signatures of two typical features of an ONeMg-core supernova: a sharp step in the density profile at the base of the He shell and a faster shock wave propagation compared to iron core supernovae. Before the shock hits the density step (t ~ 150 ms), the survival probability of electron neutrinos is about 0.68, in contrast to values of 0.32 or less for an iron core supernova. The passage of the shock through the step and its subsequent propagation cause a decrease of the survival probability and a decrease of the amplitude of oscillations in the Earth, reflecting the transition to a more adiabatic propagation inside the star. These changes affect the lower energy neutrinos first; they are faster and more sizable for larger theta_13. They are unique of ONeMg-core supernovae, and give the possibility to test the speed of the shock wave. The time modulation of the Earth effect and its negative sign at the neutronization peak are the most robust signatures in a detector.

I address the question of what can be learned from the observation of the diffuse supernova neutrino flux in the precision phase, at next generation detectors of Megaton scale. An analytical study of the spectrum of the diffuse flux shows that, above realistic detection thresholds of 10 MeV or higher, the spectrum essentially reflects the exponential-times-polynomial structure of the original neutrino spectrum at the emission point. There is only a weak (tens of per cent) dependence on the power \beta describing the growth of the supernova rate with the redshift. Different original neutrino spectra correspond to large differences in the observed spectrum of events at a water Cerenkov detector: for typical supernova rates, the ratio of the numbers of events in the first and second energy bins (of 5 MeV width) varies in the interval 1.5 - 4.3 for pure water (energy threshold 18 MeV) and in the range 1 - 2.5 for water with Gadolinium (10 MeV threshold). In the first case discrimination would be difficult due to the large errors associated with background. With Gadolinium, instead, the reduction of the total error down to 10-20 % level would allow spectral sensitivity, with a dramatic improvement of precision with respect to the SN1987A data. Even in this latter case, for typical neutrino luminosity the dependence on \beta is below sensitivity, so that it can be safely neglected in data analysis.

I review the status and perspectives of the research on the diffuse flux of (core collapse) supernova neutrinos (DF). In absence of a positive signal, several upper bounds exist in different detection channels. Of these, the strongest is the limit from SuperKamiokande (SK) of 1.2 electron antineutrinos cm^-2 s^-1 at 90% confidence level above 19.3 MeV of neutrino energy. The predictions of the DF depend on the cosmological rate of supernovae and on the neutrino emission in a individual supernova (spectrum, luminosity,..). Above the SK threshold, they range between 0.05 electron antineutrinos cm^-2 s^-1 up to touching the SK limit. The SK upper bound constrains part of the parameter space of the supernova rate - and indirectly of the star formation rate - only in models with relatively hard neutrino spectra, while predictions with softer spectra would need bounds stronger by about a factor of 4 to be tested. Experimentally, a feasible and very important goal for the future is the improvement of background discrimination and the resulting lowering of the detection threshold. Theory instead will benefit from reducing the uncertainties on the supernova neutrino emission (either with more precise numerical modeling or with data from a galactic supernova) and on the supernova rate. The latter will be provided precisely by next generation supernova surveys up to a normalization factor. Therefore, the detection of the DF is likely to be precious chiefly to constrain such normalization and to study the physics of neutrino emission in supernovae.

I calculate the diffuse flux of electron antineutrinos from all supernovae using the information on the neutrino spectrum from SN987A and the information on the rate of supernovae from direct supernova observations. The interval of flux allowed at 99 % confidence level is 0.05 - 0.35 cm^-2 s ^-1 above the SuperKamiokande (SK) energy cut of 19.3 MeV. This result is at least a factor of 4 smaller than the current SK upper limit of 1.2 cm^-2 s ^-1, thus motivating the experimental efforts to lower the detection energy threshold or to upgrade to higher volumes. A Megaton water Cherenkov detector with 90 % efficiency would record 2 - 44 inverse beta decay events a year depending on the energy cut.

I derive an upper bound on the electron neutrino component of the diffuse supernova neutrino flux from the constraint on the antineutrino component at SuperKamiokande. The connection between antineutrino and neutrino channels is due to the similarity of the muon and tau neutrino and antineutrino fluxes produced in a supernova, and to the conversion of these species into electron neutrinos and antineutrinos inside the star. The limit on the electron neutrino flux is 5.5 cm^-2 s^-1 above 19.3 MeV of neutrino energy, and is stronger than the direct limit from Mont Blanc by three orders of magnitude. It represents the minimal sensitivity required at future direct searches, and is intriguingly close to the reach of the Sudbury Neutrino Observatory (SNO) and of the ICARUS experiment. The electron neutrino flux will have a lower bound if the electron antineutrino flux is measured. Indicatively, the first can be smaller than the second at most by a factor of 2-3 depending on the details of the neutrino spectra at production.

After recent results from solar neutrino experiments and KamLAND we know that neutrinos from SN1987A underwent flavor conversion. This must be taken into account in the analysis of the data. Assuming the normal mass hierarchy of neutrinos we calculate the permutation factors p for the Kamiokande-2, IMB and Baksan detectors. The conversion inside the star leads to p = 0.28 - 0.32; regeneration effects inside the Earth reduce this factor down to 0.15 - 0.20 (at E = 40 MeV). We study in details the influence of conversion on the observed signal depending on the parameters of the original neutrino spectra. For a given set of these parameters, the conversion could increase the average energy of the observed events by up to 50 % and the number of events by a factor of 2 at Kamiokande-2 and of 3 - 5 at IMB. Inversely, we find that neglecting the conversion effects can lead up to 50 % error in the determination of the average energy of the original antinu_e flux and about 50 % error on its luminosity. By comparison with experimental data we conclude that the Kamiokande-2 data alone do not favor strong conversion effect. In contrast, the combined analysis of the Kamiokande and IMB results slightly favors strong conversion effects (that is, large difference of the original spectra). In comparison with the no oscillation case, the latter require lower average energy and higher luminosity of the original antinu_e flux.

We consider in details the effects of the 13-mixing (sin^2 theta_{13}) and of the type of mass hierarchy/ordering (sign[ Delta m^2_{13}]) on neutrino signals from the gravitational collapses of stars. The observables (characteristics of the energy spectra of nu_e and antinu_e events) sensitive to sin^2 theta_{13} and sign[Delta m^2_{13}] have been calculated. They include the ratio of average energies of the spectra, r_E = <E>/<\bar E >, the ratio of widths of the energy distributions, r_Gamma, the ratios of total numbers of nu_e and antinu_e events at low energies, S, and in the high energy tails, R_{tail}. We construct and analyze scatter plots which show the predictions for the observables for different intervals of sin^2 theta_{13} and signs of Delta m^2_{13}, taking into account uncertainties in the original neutrino spectra, the star density profile, etc.. Regions in the space of observables r_E, r_Gamma, S, R_{tail} exist in which certain mass hierarchy and intervals of sin^2 theta_{13} can be identified or discriminated. We elaborate on the method of the high energy tails in the spectra of events. The conditions are formulated for which sin^2 theta_{13} can be (i) measured, (ii) restricted from below, (iii) restricted from above. We comment on the possibility to determine sin^2 theta_{13} using the time dependence of the signals due to the propagation of the shock wave through the resonance layers of the star. We show that the appearance of the delayed Earth matter effect in one of the channels (nu_e or antinu_e) in combination with the undelayed effect in the other channel will allow to identify the shock wave appeareance and determine the mass hierarchy.

The formalism of flavor conversion of supernova neutrinos is generalized to include possible differences in the fluxes of the muon and tau neutrinos produced in the star. In this case the radiatively induced difference of the nu_mu and nu_tau potentials in matter becomes important. The nu_mu and nu_tau flux differences can manifest themselves in the effects of the Earth matter on the observed nu_e (antinu_e) signal if: (i) the neutrino mass hierarchy is normal (inverted); (ii) the solution of the solar neutrino problem is in the LMA region; (iii) the mixing U_{e3} is relatively large: |U_{e3}|>10^{-3}. We find that for differences in the nu_mu - nu_tau (antinu_mu - antinu_tau) average energies and/or integrated luminosities < 20 %, the relative deviation of the observed nu_e (antinu_e) energy spectrum at E > 50 MeV from that in the case of the equal fluxes can reach 20 - 30 % (10 - 15 %) for neutrinos crossing the Earth. It could be detected in future if large detectors sensitive to the nu_e (antinu_e) energy spectrum become available. The study of this effect would allow one to test the predictions of the nu_mu, nu_tau, antinu_mu, antinu_tau fluxes from supernova models and therefore give an important insight into the properties of matter at extreme conditions. It should be taken into account in the reconstruction of the neutrino mass spectrum and mixing matrix from the supernova neutrino observations. We show that even for unequal nu_mu and nu_tau fluxes, effects of leptonic CP violation can not be studied in the supernova neutrino experiments.

We perform a detailed study of the Earth matter effects on supernova neutrinos. The dependences of these effects on the properties of the original neutrino fluxes, on the trajectory of the neutrinos inside the Earth and on the oscillation parameters are described. We show that, for a large fraction (60 %) of the possible arrival times of the signal, the neutrino flux crosses a substantial amount of the matter of the Earth at least for one of the existing detectors. For oscillation parameters from the LMA solution of the solar neutrino problem the Earth matter effect consists in an oscillatory modulation of the $\barnue$ and/or $\nu_e$ energy spectra. The relative deviation with respect to the undistorted spectra can be as large as 20-30 % for E> 20 MeV and 70-100 % for E> 40 MeV. For parameters from the SMA and LOW solutions the effect is localized at low energies (E<10 MeV) and is not larger than 10 %. The Earth matter effects can be revealed (i) by the observation of oscillatory distortions of the energy spectra in a single experiment and (ii) by the comparison between the spectra at different detectors. For a supernova at distance D=10 Kpc, comparing the results of SuperKamiokande (SK), SNO and LVD experiments one can establish the effect at (2-3) \sigma level, whereas larger statistical significance ((4-5) \sigma) is obtained if two experiments of SK-size or larger are available. Studies of the Earth matter effect will select or confirm the solution of the solar neutrino problem, probe the mixing U_{e 3} and identify the hierarchy of the neutrino mass spectrum.

We study matter effects on oscillations of high-energy neutrinos in the Universe. Substantial effect can be produced by scattering of the neutrinos from cosmological sources ($z\gta 1$) on the relic neutrino background, provided that the latter has large CP-asymmetry: $\eta\equiv (n_\nu-n_{\bar{\nu}})/n_\gamma\gta 1$, where $n_\nu$, $n_{\bar{\nu}}$ and $n_\gamma$ are the concentrations of neutrinos, antineutrinos and photons. We consider in details the dynamics of conversion in the expanding neutrino background. Applications are given to the diffuse fluxes of neutrinos from GRBs, AGN, and the decay of super-heavy relics. We find that the vacuum oscillation probability can be modified by $\sim (10-20)%$ and in extreme cases allowed by present bounds on $\eta$ the effect can reach $\sim 100%$. Signatures of matter effects would consist (i) for both active-active and active-sterile conversion, in a deviation of the numbers of events produced in a detector by neutrinos of different flavours, $N_{\alpha}~(\alpha=e,\mu,\tau)$, and of their ratios from the values given by vacuum oscillations; such deviations can reach $\sim 5-15%$, (ii) for active-sterile conversion, in a characteristic energy dependence of the ratios $N_{e}/N_{\mu},N_{e}/N_{\tau},N_{\mu}/N_{\tau}$. Searches for these matter effects will probe large CP and lepton asymmetries in the universe.

We study properties of the oscillation effects in the matter of the Earth on antineutrino fluxes from supernovae. We show that these effects can provide explanation of the difference in the energy spectra of the events detected by Kamiokande-2 and IMB detectors from SN1987A as well as the absence of high-energy events with $E\gta 40$ MeV. This explanation requires the neutrino oscillation parameters $\Delta m^2$ and $\sin^2 2\theta $ to be in the region of the LMA solution of the solar neutrino problem and the normal mass hierarchy if $|U_{e 3}|^2\gta 10^{-3}$. The hierarchy can be inverted if $|U_{e 3}|^2\ll 10^{-3}$. The solution of the solar neutrino problem based on $\nu_e$-conversion to a pure sterile state is disfavoured by SN1987A data.