The Diffuse Supernova Neutrino Background

About one core-collapse supernova occurs every second somewhere in the observable universe. Each releases approximately $3\times 10^{53}$ erg of energy as neutrinos, with the burst lasting about 10 seconds. Integrated over cosmological time and across all the supernovae that have ever exploded, the resulting accumulated neutrino flux at Earth is approximately 10–20 antineutrinos per square centimetre per second above 10 MeV — the Diffuse Supernova Neutrino Background, or DSNB.

The DSNB is not yet detected. At standard predictions it is just below the sensitivity of the current generation of detectors. But by the late 2020s, with the Super-Kamiokande Gadolinium upgrade and JUNO’s commissioning, the first definitive detection is widely expected. The DSNB measurement would become the first observation of the integrated neutrino history of the universe — a cosmological probe of the supernova rate over redshifts up to z ~ 5, completely independent of optical surveys.

This post describes the physics of the DSNB, how detectors are searching for it, and what its eventual detection will tell us.

Where the DSNB comes from

A core-collapse supernova produces neutrinos in three temporal phases:

• Neutronization burst (first 10 ms): A short, intense ν_e pulse from electron capture on protons during core collapse
• Accretion phase (10 ms – 0.5 s): Continued neutrino emission as material falls onto the proto-neutron star
• Cooling phase (0.5 s – 10 s): The proto-neutron star radiates its thermal energy as a near-equal mix of all six neutrino flavors

The total energy carried away by neutrinos is about $3\times 10^{53}$ erg per supernova — 99% of the gravitational binding energy released by the collapse. The energy is roughly equipartitioned among the six species (ν_e, ν̄_e, ν_μ, ν̄_μ, ν_τ, ν̄_τ), with ν̄_e carrying about 17% of the total, with mean energy ~15 MeV and a quasi-thermal spectrum extending to roughly 30 MeV.

Now consider this happening at every redshift in the observable universe. Locally, the volumetric core-collapse-supernova rate is about $10^{-4}$ per Mpc³ per year. Over the past 10 billion years, accounting for the cosmic star-formation history that peaks at $z\sim 2$, the integrated accumulated supernova count is approximately $10^{19}$ — a vast population.

The neutrino flux from any single distant supernova at Earth is tiny (the SN1987A burst, at 50 kpc, gave 24 events worldwide). But the cumulative isotropic flux from all past supernovae is the DSNB.

The expected spectrum and rate

The DSNB spectrum at Earth is a redshifted superposition of all the past supernova emission spectra: $\frac{d{\Phi }_{\nu }}{dE}=\frac{c}{H_0}{\int }_0^{z_{max}}\frac{R_{SN}(z)}{1+z}\cdot \frac{d{N}_{\nu }}{d{E}^{\prime }}\cdot \frac{d{E}^{\prime }}{dE}\cdot \frac{dz}{E(z)}$ where $R_{SN}(z)$ is the volumetric supernova rate at redshift $z$, $d{N}_{\nu }/d{E}^{\prime }$ is the per-supernova spectrum at the source, $E^{\prime }=E(1+z)$ is the rest-frame energy, $E(z)$ is the cosmological expansion function, and the integration runs over all redshifts where supernovae could have occurred.

Numerical evaluation of this integral, using a typical supernova spectrum and the standard cosmic star-formation history, gives:

• ${\Phi }_{\nu }\approx 10$–$20$ cm⁻² s⁻¹ for ν̄_e above 10 MeV
• Spectrum peaks around 8–10 MeV (after redshifting from the source-frame peak of ~15 MeV)
• Tail extends to ~30 MeV

The most relevant detection energy window is 10–30 MeV. Below 10 MeV, atmospheric neutrino backgrounds and reactor antineutrinos dominate. Above 30 MeV, atmospheric neutrinos dominate. The 10–30 MeV window is the “sweet spot” for DSNB detection.

Detection: inverse beta decay on hydrogen

The standard detection reaction is the same as for reactor antineutrinos: ${\bar{\nu }}_e+p\to e^{+}+n$ For DSNB, the antineutrino energies are higher (10–30 MeV vs. reactor 1–10 MeV), but the kinematic structure is similar. Each event produces a delayed coincidence of (positron annihilation) + (neutron capture).

Two backgrounds dominate the DSNB-energy region:

Atmospheric neutrinos: Cosmic-ray-produced atmospheric neutrinos contribute several events per year to the same energy window. They can be partially discriminated by their characteristic event topology.

Spallation backgrounds: Cosmic-ray muons interacting with detector materials produce isotopes whose decays mimic the inverse-beta-decay coincidence. These backgrounds dominate the lowest-energy region (below 16 MeV) and limit the practical energy threshold.

Reactor antineutrinos: Tail of the reactor spectrum extends slightly into the DSNB region, but is sub-dominant above 10 MeV.

The combined background rate, integrated over the 10–30 MeV window, is comparable to or exceeds the predicted DSNB rate at most current detectors. Substantial background suppression is needed.

Super-Kamiokande Gd era

A breakthrough came in 2020 when Super-Kamiokande began doping its 50 kt of pure water with gadolinium sulfate. The Gd captures thermal neutrons with a much higher probability than hydrogen, releasing a characteristic 8 MeV gamma-ray cascade that is easily identified. With Gd doping, neutrons from inverse beta decay become taggable with high efficiency.

This dramatically improves background rejection: most of the spallation backgrounds do not produce associated neutrons, so the Gd-tagging requirement removes them. Atmospheric backgrounds with neutrons are reduced by topology cuts.

Super-Kamiokande’s first dedicated DSNB analysis with Gd doping was published in 2024. It set an upper limit on the DSNB ν̄_e flux at: ${\Phi }_{DSNB}<2.0{\text{ cm}}^{-2}{\text{s}}^{-1}\text{ above 16 MeV}(90\%\text{ C.L.})$ This is approximately a factor of 2 above the standard prediction of about 0.8–1.2 cm⁻²s⁻¹ in the same energy window. With additional running and improved analysis, Super-Kamiokande is expected to either claim a first detection or set a stringent constraint within a few years.

JUNO’s role

The Jiangmen Underground Neutrino Observatory, scheduled for first physics in 2026, will be a 20-kiloton liquid-scintillator detector with excellent energy resolution and intrinsic neutron tagging through gadolinium loading or boron-10. JUNO’s dedicated DSNB analysis is expected to detect the standard signal at ~10 events per year against ~3 background events per year — a definitive detection within 1–2 years of operation.

JUNO’s combination of sensitivity and energy resolution will allow not just DSNB detection but spectral measurements that constrain:

• The supernova rate at high redshift
• The mean ν̄_e energy from a typical supernova (currently estimated from 1987A and theory)
• Failed-supernova fraction (cores that collapse to black holes without producing visible explosions)

Hyper-Kamiokande and DUNE

Hyper-Kamiokande, scheduled for 2027, will scale up the Super-K Gd approach by a factor of 11.5 in target mass. Its DSNB sensitivity reaches predicted levels by a wide margin. Within 5 years of operation, Hyper-K could measure the DSNB flux to ~10% precision.

DUNE, with liquid argon, has a different detection channel: ${\nu }_e+{}^{40}\text{Ar}\to e^{-}+{}^{40}{\text{K}}^{\ast }$. This gives DUNE specific sensitivity to the ν_e component of the DSNB — complementing the inverse-beta-decay-only ν̄_e measurement at water and scintillator detectors. Combined ν_e + ν̄_e measurements will constrain the supernova-population spectral parameters.

What the DSNB tells us

A detected DSNB will provide unique scientific value:

First measurement of the supernova rate at high z. Optical supernova surveys are sensitive only to z < 2 (redshift around 10 billion years ago). The DSNB integrates over all redshifts and provides an independent constraint on $R_{SN}$ at z = 2–5, where the cosmic star-formation rate peaks but optical observations are difficult.

Constraints on the failed-supernova fraction. Some massive stars collapse to black holes without producing a visible supernova (“failed supernovae”). These events still produce a neutrino burst — actually a brighter ν_e burst since the proto-neutron star is shorter-lived. The DSNB integrates these contributions; comparing DSNB to optical-survey supernova rates constrains the failed-supernova fraction.

Precision measurement of supernova spectra. The integrated DSNB shape encodes the average neutrino emission spectrum from a typical supernova. Fitting it constrains parameters like the mean neutrino energy, the energy hierarchy between flavors, and the neutrino-driven wind temperature.

Tests of beyond-Standard-Model physics. The DSNB provides a long-baseline neutrino-propagation experiment over cosmological distances. Effects like neutrino decay, non-standard interactions in matter, or sterile-neutrino oscillations modify the DSNB spectrum in calculable ways. A precise DSNB measurement is a unique test of these scenarios.

Outlook

The DSNB is the next major neutrino discovery within reach of current technology. By 2027–2030, multiple experiments with complementary detection channels will have either confirmed the standard prediction or revealed unexpected physics.

The detection itself will be a quiet milestone — no single dramatic event, just a slow accumulation of statistical excess above background — but it will represent the first measurement of an integrated neutrino flux of cosmic origin. Solar neutrinos come from one star. Supernova-1987A neutrinos came from one supernova. Reactor neutrinos come from a few human-built reactors. The DSNB comes from the entire universe of supernovae that have ever happened.

That makes it, in the sense of cosmological context, the most expansive neutrino measurement that human technology has ever attempted.