---
title: Reactor Antineutrino Oscillation Findings
url: https://www.emergentmind.com/topics/reactor-antineutrino-oscillation-result
type: topic
---

# Reactor Antineutrino Oscillation Findings

A reactor antineutrino oscillation result is a measurement quantifying the disappearance (and at longer baselines, spectral distortion) of electron antineutrinos ($\overline{\nu}_e$) emitted by nuclear reactors, as they propagate over distances ranging from tens of meters to several hundred kilometers. These results probe the parameters of the lepton mixing matrix (PMNS), primarily the mixing angles $\theta_{13}$, $\theta_{12}$ and the mass-squared differences $\Delta m^2_{ee}$ (or $\Delta m^2_{32}$) and $\Delta m^2_{21}$, by measuring the survival probability of $\overline{\nu}_e$ as a function of energy and baseline. Reactor experiments, through optimized detection channels, background suppression, and relative or absolute spectral analysis, provide some of the most precise constraints on the neutrino sector, serving as cornerstones for three-flavor oscillation phenomenology and for the exploration of anomalies beyond the Standard Model.

## 1. Physical Framework and Survival Probability

In the three-neutrino framework, the reactor antineutrino survival probability at baseline $L$ and energy $E$ is given by
\[
P(\overline{\nu}_e \rightarrow \overline{\nu}_e) = 1
- \sin^2 2\theta_{13} \, \sin^2 \Delta_{ee}
- \cos^4 \theta_{13} \sin^2 2\theta_{12} \sin^2 \Delta_{21},
\]
where $\Delta_{ij} \equiv 1.267 \, \Delta m^2_{ij} \, L / E$ (with $\Delta m^2_{ij}$ in eV$^2$, $L$ in meters, $E$ in MeV), and $\Delta m^2_{ee} \equiv \cos^2\theta_{12} \, \Delta m^2_{31} + \sin^2\theta_{12} \, \Delta m^2_{32}$ is the electron-flavor-projected effective mass splitting [1806.00248].

For kilometer-scale baselines, the dominant oscillation is at the $|\Delta m^2_{ee}| \sim 2.5 \times 10^{-3}\,\mathrm{eV}^2$ "atmospheric" scale, governed by $\sin^2 2\theta_{13}$. At baselines of order 100–1000 km, the energy spectrum is modulated by the solar-scale mass splitting $\Delta m^2_{21} \sim 7.5 \times 10^{-5}\,\mathrm{eV}^2$ and mixing angle $\theta_{12}$. The standard detection channel is inverse beta decay (IBD): $\overline{\nu}_e + p \to e^+ + n$, with the positron producing a prompt signal and the neutron captured on Gd or H producing a delayed $\gamma$.

## 2. Experimental Apparatus and Data Acquisition

Most reactor oscillation measurements employ segmented liquid-scintillator detectors equipped with PMTs and optimized for delayed-coincidence background rejection. Near-far arrangement—simultaneously deploying identical detectors at $L_\text{near} \sim 300$–$500$ m and $L_\text{far} \sim 1.4$–$1.6$ km (e.g., RENO, Daya Bay)—enables cancellation of correlated reactor and detector systematics [1806.00248, 1505.03456]. High-statistics samples are accumulated over thousands of live days; for example, RENO recorded $103\,212$ (far) and $850\,666$ (near) candidates after all cuts with a $4.7\%$ ($2.0\%$) background fraction [1806.00248]. Daya Bay's $8$-AD configuration collected $>3.9 \times 10^6$ IBDs over $1958$ days [1809.02261].

Critical detector parameters include:
- Target volume (e.g., $16$–$20$ t LAB-based LS with 0.1% Gd loading),
- Overburden (minimizing cosmogenic backgrounds),
- muon veto systems,
- fine-grained energy calibration (calibration sources, cosmogenic spectra, spatial corrections).

For longer-baseline measurements of $\Delta m^2_{21}$, very large scintillator detectors (KamLAND, SNO+, JUNO) at several hundred km baselines, and water Cherenkov detectors (Super-Kamiokande with Gd doping), are used [2511.11856, 2006.01155].

## 3. Statistical Methodologies and Systematic Uncertainties

The canonical analysis method is a far-to-near prompt-energy spectral ratio, which cancels flux and detection-model uncertainties to leading order:
\[
R(E_p) = \frac{N_\text{far}(E_p) - B_\text{far}(E_p)}{N_\text{pred,far}(E_p)},
\]
where $N_\text{pred,far}(E_p)$ is built by scaling the near spectrum to the far baseline without oscillation, thus suppressing reactor-model and detector-correlated systematic uncertainties [1806.00248].

Parameter extraction proceeds via binned or unbinned maximum-likelihood or $\chi^2$ methods, incorporating full statistical and systematic covariance matrices. Pull terms accommodate uncertainties on uncorrelated reactor flux ($0.9\%$ typical), detection efficiency ($0.21\%$ relative), energy-scale ($0.15\%$), and background rates [1806.00248].

Typical dominant systematic errors:
- For $\sin^2 2\theta_{13}$: detection efficiency and reactor flux,
- For $|\Delta m^2_{ee}|$: absolute and relative energy scale.

Error correlations are characterized: in RENO, the $\sin^2 2\theta_{13}$ vs. $|\Delta m^2_{ee}|$ contour is mildly anti-correlated ($\rho \approx -0.2$) [1806.00248].

Backgrounds are accurately measured from sideband data or reactor-off running, and include accidentals, cosmogenic $^9$Li/$^8$He, fast neutrons, and $^{252}$Cf contamination. In SNO+, cosmogenic and radiogenic backgrounds are further constrained by likelihood discriminators and pulse-shape analysis [2511.11856, 2505.04469].

## 4. Key Oscillation Parameter Measurements

The spectral distortion as a function of $L/E$ directly yields the oscillation amplitude and frequency. Representative results for leading experiments:

| Experiment        | $\sin^2 2\theta_{13}$        | $|\Delta m^2_{ee}|$ ($10^{-3}$ eV$^2$) | Statistics         | Ref.         |
|-------------------|-----------------------------|-----------------------------------------|--------------------|--------------|
| RENO ($\sim2200$ d, Gd)       | $0.0896 \pm 0.0048_{stat} \pm 0.0047_{syst}$ | $2.68 \pm 0.12_{stat} \pm 0.07_{syst}$ | $>950$k IBDs       | [1806.00248] |
| Daya Bay (nGd, 3158 d)         | $0.0851 \pm 0.0024$                          | $\pm 0.060$ (on $|\Delta m^2_{32}|$)   | $5.55$M IBDs       | [2409.14523] |
| Daya Bay (nH, 1958 d)          | $0.0759^{+0.0050}_{-0.0049}$                 | $2.72^{+0.14}_{-0.15}$                 | $3.6$M IBDs        | [2406.01007] |
| Combined Daya Bay (nGd+nH)     | $0.0833 \pm 0.0022$                          | —                                     | —                  | [2406.01007] |
| Double Chooz (101 d)           | $0.086 \pm 0.041_{stat} \pm 0.030_{syst}$    | Fixed                                  | 4121 IBDs          | [1112.6353]  |
| SNO+ ($1.46$ kt·yr, 2022–2025) | —                                           | $0.00793^{+0.00021}_{-0.00024}$ ($\Delta m^2_{21}$) | $>100$ events | [2511.11856] |

The observed prompt-energy spectral ratio $R(E_p)$ exhibits an oscillatory deficit dipping to $R \sim 0.85$ at $E_p\approx 4$–$5$ MeV in RENO [1806.00248], and the $L/E$ distribution reveals a survival probability minimum at $L/E \approx 500$ m/MeV.

The parameter $\sin^2 2\theta_{13}$ is now measured with sub-4% relative precision (Daya Bay nGd sample), and $|\Delta m^2_{ee}|$ with better than 3% relative error [2409.14523, 1809.02261]. The SNO+ measurement of $\Delta m^2_{21}$, $7.93^{+0.21}_{-0.24} \times 10^{-5}$ eV$^2$, approaches KamLAND's precision and is consistent with the global fit [2511.11856].

## 5. Reactor Oscillation Results in Global and Comparative Context

The precision reactor measurements are consistent with the three-flavor oscillation paradigm and serve as reference points for global fits. Comparative values:

- Daya Bay (2017): $\sin^2 2\theta_{13} = 0.0841\pm0.0027$, $|\Delta m^2_{ee}|=(2.50\pm0.08)\times 10^{-3}$ eV$^2$,
- RENO (2018): $\sin^2 2\theta_{13} = 0.0896 \pm 0.0068$, $|\Delta m^2_{ee}| = 2.68 \pm 0.14$,
- Double Chooz: $\sin^2 2\theta_{13} = 0.090 \pm 0.035$,
- Accelerator (T2K, MINOS, NO$\nu$A): $|\Delta m_{32}^2| \approx (2.52\pm 0.07)\times 10^{-3}$ eV$^2$,
- SNO+ (2025): $\Delta m^2_{21} = (7.93^{+0.21}_{-0.24})\times 10^{-5}$ eV$^2$ [2511.11856].

All results are mutually consistent within uncertainties and now dominate the precision on $\theta_{13}$ and $\Delta m^2_{ee}$; SNO+ and KamLAND set the standard for $\Delta m^2_{21}$ and, in global combinations, for $\sin^2 \theta_{12}$ [2511.11856].

## 6. Reactor Antineutrino Anomaly, Sterile Oscillation Searches, and Spectral Features

A persistent $\sim 6\%$ deficit relative to the reevaluated reactor flux model, termed the reactor antineutrino anomaly, motivated intensive scrutiny [1101.2755]. Short-baseline experiments (STEREO, PROSPECT, SoLid) and radioactive-source deployments at KamLAND (CeLAND) exclude the parameter space ($\sin^2 2\theta_{new} \sim 0.1$, $\Delta m^2_{new} \sim 1$ eV$^2$) favored by the anomaly, with no significant evidence for eV-scale sterile neutrinos [1905.11896, 1309.6805]. Daya Bay and RENO have set additional constraints by searching for anomalous spectral distortions at short $L/E$ [2409.14523].

A prominent, reactor-power-correlated excess near $E_p\sim 5$ MeV, observed in RENO, Daya Bay, and Double Chooz, constitutes a significant spectral deviation from predicted models, highlighting deficiencies in $\overline{\nu}_e$ production calculations for specific actinides ($^{235}$U, $^{239}$Pu) or beta-decay branch modeling [1410.7987]. This excess is spectrally isolated and cancels in near-far differences, preserving the oscillation measurements but demanding further theoretical and experimental attention.

## 7. Impact, Prospects, and Future Experiments

Reactor antineutrino oscillation results determine the smallest PMNS angle $\theta_{13}$ and the atmospheric mass-squared splitting with the highest available precision, fundamentally enabling the CP-violation and mass-ordering programs of long-baseline accelerator (DUNE, T2HK) and medium-baseline reactor (JUNO, RENO-50) experiments [1806.00248, 1505.03456]. The simultaneous measurement of $\sin^2 2\theta_{13}$ and $\Delta m^2_{ee}$ sharpens global oscillation fits, reduces parameter correlations (notably with $\delta_{CP}$), and improves predictions for $\overline{\nu}_e$ appearance and disappearance at all energies and baselines.

The continued improvement in event statistics, energy calibration, and background modeling, together with next-generation detectors (full SNO+ with higher light yield, JUNO with 20 kt LS and $<1\%$ resolution), is expected to further reduce uncertainties on all mixing parameters, test the unitarity of the PMNS matrix, and probe new physics scenarios at subpercent levels [2511.11856]. The precision achieved sets a benchmark for the field and ensures robust cross-comparison between disappearance and appearance channels in the global neutrino oscillation program.

Source: https://www.emergentmind.com/topics/reactor-antineutrino-oscillation-result