Reactor Antineutrino Oscillation Findings

• Reactor antineutrino oscillation is the measurement of electron antineutrino disappearance and spectral distortion over varying baselines using the three-flavor framework.
• It employs segmented liquid-scintillator detectors, inverse beta decay detection, and rigorous statistical analyses to extract parameters like θ13 and Δm²ee with sub-4% precision.
• The results benchmark the standard oscillation model and guide future experiments on neutrino CP-violation and mass-ordering in global neutrino research.

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 $\theta_{13}$0 (with $\theta_{13}$1 in eV$\theta_{13}$2, $\theta_{13}$3 in meters, $\theta_{13}$4 in MeV), and $\theta_{13}$5 is the electron-flavor-projected effective mass splitting (Bak et al., 2018).

For kilometer-scale baselines, the dominant oscillation is at the $\theta_{13}$6 "atmospheric" scale, governed by $\theta_{13}$7. At baselines of order 100–1000 km, the energy spectrum is modulated by the solar-scale mass splitting $\theta_{13}$8 and mixing angle $\theta_{13}$9. The standard detection channel is inverse beta decay (IBD): $\theta_{12}$0, with the positron producing a prompt signal and the neutron captured on Gd or H producing a delayed $\theta_{12}$1.

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 $\theta_{12}$2–$\theta_{12}$3 m and $\theta_{12}$4–$\theta_{12}$5 km (e.g., RENO, Daya Bay)—enables cancellation of correlated reactor and detector systematics (Bak et al., 2018, collaboration et al., 2015). High-statistics samples are accumulated over thousands of live days; for example, RENO recorded $\theta_{12}$6 (far) and $\theta_{12}$7 (near) candidates after all cuts with a $\theta_{12}$8 ($\theta_{12}$9) background fraction (Bak et al., 2018). Daya Bay's $\Delta m^2_{ee}$0-AD configuration collected $\Delta m^2_{ee}$1 IBDs over $\Delta m^2_{ee}$2 days (collaboration et al., 2018).

Critical detector parameters include:

• Target volume (e.g., $\Delta m^2_{ee}$3–$\Delta m^2_{ee}$4 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_{ee}$5, very large scintillator detectors (KamLAND, SNO+, JUNO) at several hundred km baselines, and water Cherenkov detectors (Super-Kamiokande with Gd doping), are used (Abreu et al., 14 Nov 2025, Gouvêa et al., 2020).

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: $\Delta m^2_{ee}$6 where $\Delta m^2_{ee}$7 is built by scaling the near spectrum to the far baseline without oscillation, thus suppressing reactor-model and detector-correlated systematic uncertainties (Bak et al., 2018).

Parameter extraction proceeds via binned or unbinned maximum-likelihood or $\Delta m^2_{ee}$8 methods, incorporating full statistical and systematic covariance matrices. Pull terms accommodate uncertainties on uncorrelated reactor flux ($\Delta m^2_{ee}$9 typical), detection efficiency ($\Delta m^2_{32}$0 relative), energy-scale ($\Delta m^2_{32}$1), and background rates (Bak et al., 2018).

Typical dominant systematic errors:

• For $\Delta m^2_{32}$2: detection efficiency and reactor flux,
• For $\Delta m^2_{32}$3: absolute and relative energy scale.

Error correlations are characterized: in RENO, the $\Delta m^2_{32}$4 vs. $\Delta m^2_{32}$5 contour is mildly anti-correlated ($\Delta m^2_{32}$6) (Bak et al., 2018).

Backgrounds are accurately measured from sideband data or reactor-off running, and include accidentals, cosmogenic $\Delta m^2_{32}$7Li/$\Delta m^2_{32}$8He, fast neutrons, and $\Delta m^2_{32}$9Cf contamination. In SNO+, cosmogenic and radiogenic backgrounds are further constrained by likelihood discriminators and pulse-shape analysis (Abreu et al., 14 Nov 2025, Collaboration et al., 7 May 2025).

4. Key Oscillation Parameter Measurements

The spectral distortion as a function of $\Delta m^2_{21}$0 directly yields the oscillation amplitude and frequency. Representative results for leading experiments:

Experiment	$\Delta m^2_{21}$1	$\Delta m^2_{21}$2 ($\Delta m^2_{21}$3 eV$\Delta m^2_{21}$4)	Statistics	Ref.
RENO ($\Delta m^2_{21}$5 d, Gd)	$\Delta m^2_{21}$6	$\Delta m^2_{21}$7	$\Delta m^2_{21}$8k IBDs	(Bak et al., 2018)
Daya Bay (nGd, 3158 d)	$\Delta m^2_{21}$9	$\overline{\nu}_e$0 (on $\overline{\nu}_e$1)	$\overline{\nu}_e$2M IBDs	(Li, 2024)
Daya Bay (nH, 1958 d)	$\overline{\nu}_e$3	$\overline{\nu}_e$4	$\overline{\nu}_e$5M IBDs	(collaboration et al., 2024)
Combined Daya Bay (nGd+nH)	$\overline{\nu}_e$6	—	—	(collaboration et al., 2024)
Double Chooz (101 d)	$\overline{\nu}_e$7	Fixed	4121 IBDs	(Abe et al., 2011)
SNO+ ($\overline{\nu}_e$8 kt·yr, 2022–2025)	—	$\overline{\nu}_e$9 ($L$0)	$L$1 events	(Abreu et al., 14 Nov 2025)

The observed prompt-energy spectral ratio $L$2 exhibits an oscillatory deficit dipping to $L$3 at $L$4–$L$5 MeV in RENO (Bak et al., 2018), and the $L$6 distribution reveals a survival probability minimum at $L$7 m/MeV.

The parameter $L$8 is now measured with sub-4% relative precision (Daya Bay nGd sample), and $L$9 with better than 3% relative error (Li, 2024, collaboration et al., 2018). The SNO+ measurement of $E$0, $E$1 eV$E$2, approaches KamLAND's precision and is consistent with the global fit (Abreu et al., 14 Nov 2025).

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): $E$3, $E$4 eV$E$5,
• RENO (2018): $E$6, $E$7,
• Double Chooz: $E$8,
• Accelerator (T2K, MINOS, NO$E$9A): $$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},$$0 eV$$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},$$1,
• SNO+ (2025): $$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},$$2 eV$$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},$$3 (Abreu et al., 14 Nov 2025).

All results are mutually consistent within uncertainties and now dominate the precision on $$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},$$4 and $$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},$$5; SNO+ and KamLAND set the standard for $$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},$$6 and, in global combinations, for $$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},$$7 (Abreu et al., 14 Nov 2025).

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

A persistent $$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},$$8 deficit relative to the reevaluated reactor flux model, termed the reactor antineutrino anomaly, motivated intensive scrutiny (Mention et al., 2011). Short-baseline experiments (STEREO, PROSPECT, SoLid) and radioactive-source deployments at KamLAND (CeLAND) exclude the parameter space ($$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},$$9, $\theta_{13}$00 eV$\theta_{13}$01) favored by the anomaly, with no significant evidence for eV-scale sterile neutrinos (Bernard, 2019, Gando et al., 2013). Daya Bay and RENO have set additional constraints by searching for anomalous spectral distortions at short $\theta_{13}$02 (Li, 2024).

A prominent, reactor-power-correlated excess near $\theta_{13}$03 MeV, observed in RENO, Daya Bay, and Double Chooz, constitutes a significant spectral deviation from predicted models, highlighting deficiencies in $\theta_{13}$04 production calculations for specific actinides ($\theta_{13}$05U, $\theta_{13}$06Pu) or beta-decay branch modeling (Seo, 2014). 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}$07 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 (Bak et al., 2018, collaboration et al., 2015). The simultaneous measurement of $\theta_{13}$08 and $\theta_{13}$09 sharpens global oscillation fits, reduces parameter correlations (notably with $\theta_{13}$10), and improves predictions for $\theta_{13}$11 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 $\theta_{13}$12 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 (Abreu et al., 14 Nov 2025). 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.