Improved Measurement of the Evolution of the Reactor Antineutrino Flux and Spectrum at Daya Bay

Reactor neutrino experiments play a crucial role in advancing our knowledge of neutrinos. In this Letter, the evolution of the flux and spectrum as a function of the reactor isotopic content is reported in terms of the inverse-beta-decay yield at Daya Bay with 1958 days of data and improved systematic uncertainties. These measurements are compared with two signature model predictions: the Huber-Mueller model based on the conversion method and the SM2018 model based on the summation method. The measured average flux and spectrum, as well as the flux evolution with the ^{239}Pu isotopic fraction, are inconsistent with the predictions of the Huber-Mueller model. In contrast, the SM2018 model is shown to agree with the average flux and its evolution but fails to describe the energy spectrum. Altering the predicted inverse-beta-decay spectrum from ^{239}Pu fission does not improve the agreement with the measurement for either model. The models can be brought into better agreement with the measurements if either the predicted spectrum due to ^{235}U fission is changed or the predicted ^{235}U, ^{238}U, ^{239}Pu, and ^{241}Pu spectra are changed in equal measure.

Precision Measurement of Reactor Antineutrino Oscillation at Kilometer-Scale Baselines by Daya Bay

We present a new determination of the smallest neutrino mixing angle θ_{13} and the mass-squared difference Δm_{32}^{2} using a final sample of 5.55×10^{6} inverse beta-decay (IBD) candidates with the final-state neutron captured on gadolinium. This sample is selected from the complete dataset obtained by the Daya Bay reactor neutrino experiment in 3158 days of operation. Compared to the previous Daya Bay results, selection of IBD candidates has been optimized, energy calibration refined, and treatment of backgrounds further improved. The resulting oscillation parameters are sin^{2}2θ_{13}=0.0851±0.0024, Δm_{32}^{2}=(2.466±0.060)×10^{-3}  eV^{2} for the normal mass ordering or Δm_{32}^{2}=-(2.571±0.060)×10^{-3}  eV^{2} for the inverted mass ordering.

First Measurement of High-Energy Reactor Antineutrinos at Daya Bay

This Letter reports the first measurement of high-energy reactor antineutrinos at Daya Bay, with nearly 9000 inverse beta decay candidates in the prompt energy region of 8-12 MeV observed over 1958 days of data collection. A multivariate analysis is used to separate 2500 signal events from background statistically. The hypothesis of no reactor antineutrinos with neutrino energy above 10 MeV is rejected with a significance of 6.2 standard deviations. A 29% antineutrino flux deficit in the prompt energy region of 8-11 MeV is observed compared to a recent model prediction. We provide the unfolded antineutrino spectrum above 7 MeV as a data-based reference for other experiments. This result provides the first direct observation of the production of antineutrinos from several high-Q_{β} isotopes in commercial reactors.

Joint Determination of Reactor Antineutrino Spectra from ^{235}U and ^{239}Pu Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of ^{235}U and ^{239}Pu has been carried out by the Daya Bay and PROSPECT Collaborations. This Letter reports the level of consistency of ^{235}U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant ^{235}U and ^{239}Pu isotopes and improves the uncertainty of the ^{235}U spectral shape to about 3%. The ^{235}U and ^{239}Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the ^{235}U and ^{239}Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

Improved Constraints on Sterile Neutrino Mixing from Disappearance Searches in the MINOS, MINOS+, Daya Bay, and Bugey-3 Experiments

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the θ_{μe} mixing angle are derived that constitute the most constraining limits to date over five orders of magnitude in the mass-squared splitting Δm_{41}^{2}, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL_{s} for Δm_{41}^{2}<13  eV^{2}. Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL_{s} for Δm_{41}^{2}<1.6 eV^{2}.

Extraction of the ^{235}U and ^{239}Pu Antineutrino Spectra at Daya Bay

This Letter reports the first extraction of individual antineutrino spectra from ^{235}U and ^{239}Pu fission and an improved measurement of the prompt energy spectrum of reactor antineutrinos at Daya Bay. The analysis uses 3.5×10^{6} inverse beta-decay candidates in four near antineutrino detectors in 1958 days. The individual antineutrino spectra of the two dominant isotopes, ^{235}U and ^{239}Pu, are extracted using the evolution of the prompt spectrum as a function of the isotope fission fractions. In the energy window of 4-6 MeV, a 7% (9%) excess of events is observed for the ^{235}U (^{239}Pu) spectrum compared with the normalized Huber-Mueller model prediction. The significance of discrepancy is 4.0σ for ^{235}U spectral shape compared with the Huber-Mueller model prediction. The shape of the measured inverse beta-decay prompt energy spectrum disagrees with the prediction of the Huber-Mueller model at 5.3σ. In the energy range of 4-6 MeV, a maximal local discrepancy of 6.3σ is observed.

Measurement of the Electron Antineutrino Oscillation with 1958 Days of Operation at Daya Bay

We report a measurement of electron antineutrino oscillation from the Daya Bay Reactor Neutrino Experiment with nearly 4 million reactor ν[over ¯]_{e} inverse β decay candidates observed over 1958 days of data collection. The installation of a flash analog-to-digital converter readout system and a special calibration campaign using different source enclosures reduce uncertainties in the absolute energy calibration to less than 0.5% for visible energies larger than 2 MeV. The uncertainty in the cosmogenic ^{9}Li and ^{8}He background is reduced from 45% to 30% in the near detectors. A detailed investigation of the spent nuclear fuel history improves its uncertainty from 100% to 30%. Analysis of the relative ν[over ¯]_{e} rates and energy spectra among detectors yields sin^{2}2θ_{13}=0.0856±0.0029 and Δm_{32}^{2}=(2.471_{-0.070}^{+0.068})×10^{-3}  eV^{2} assuming the normal hierarchy, and Δm_{32}^{2}=-(2.575_{-0.070}^{+0.068})×10^{-3}  eV^{2} assuming the inverted hierarchy.

Evolution of the Reactor Antineutrino Flux and Spectrum at Daya Bay

The Daya Bay experiment has observed correlations between reactor core fuel evolution and changes in the reactor antineutrino flux and energy spectrum. Four antineutrino detectors in two experimental halls were used to identify 2.2 million inverse beta decays (IBDs) over 1230 days spanning multiple fuel cycles for each of six 2.9 GW_{th} reactor cores at the Daya Bay and Ling Ao nuclear power plants. Using detector data spanning effective ^{239}Pu fission fractions F_{239} from 0.25 to 0.35, Daya Bay measures an average IBD yield σ[over ¯]_{f} of (5.90±0.13)×10^{-43}  cm^{2}/fission and a fuel-dependent variation in the IBD yield, dσ_{f}/dF_{239}, of (-1.86±0.18)×10^{-43}  cm^{2}/fission. This observation rejects the hypothesis of a constant antineutrino flux as a function of the ^{239}Pu fission fraction at 10 standard deviations. The variation in IBD yield is found to be energy dependent, rejecting the hypothesis of a constant antineutrino energy spectrum at 5.1 standard deviations. While measurements of the evolution in the IBD spectrum show general agreement with predictions from recent reactor models, the measured evolution in total IBD yield disagrees with recent predictions at 3.1σ. This discrepancy indicates that an overall deficit in the measured flux with respect to predictions does not result from equal fractional deficits from the primary fission isotopes ^{235}U, ^{239}Pu, ^{238}U, and ^{241}Pu. Based on measured IBD yield variations, yields of (6.17±0.17) and (4.27±0.26)×10^{-43}  cm^{2}/fission have been determined for the two dominant fission parent isotopes ^{235}U and ^{239}Pu. A 7.8% discrepancy between the observed and predicted ^{235}U yields suggests that this isotope may be the primary contributor to the reactor antineutrino anomaly.

Limits on Active to Sterile Neutrino Oscillations from Disappearance Searches in the MINOS, Daya Bay, and Bugey-3 Experiments

Searches for a light sterile neutrino have been performed independently by the MINOS and the Daya Bay experiments using the muon (anti)neutrino and electron antineutrino disappearance channels, respectively. In this Letter, results from both experiments are combined with those from the Bugey-3 reactor neutrino experiment to constrain oscillations into light sterile neutrinos. The three experiments are sensitive to complementary regions of parameter space, enabling the combined analysis to probe regions allowed by the Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE experiments in a minimally extended four-neutrino flavor framework. Stringent limits on sin^{2}2θ_{μe} are set over 6 orders of magnitude in the sterile mass-squared splitting Δm_{41}^{2}. The sterile-neutrino mixing phase space allowed by the LSND and MiniBooNE experiments is excluded for Δm_{41}^{2}<0.8  eV^{2} at 95%  CL_{s}.

Improved Search for a Light Sterile Neutrino with the Full Configuration of the Daya Bay Experiment

This Letter reports an improved search for light sterile neutrino mixing in the electron antineutrino disappearance channel with the full configuration of the Daya Bay Reactor Neutrino Experiment. With an additional 404 days of data collected in eight antineutrino detectors, this search benefits from 3.6 times the statistics available to the previous publication, as well as from improvements in energy calibration and background reduction. A relative comparison of the rate and energy spectrum of reactor antineutrinos in the three experimental halls yields no evidence of sterile neutrino mixing in the 2×10^{-4}≲|Δm_{41}^{2}|≲0.3  eV^{2} mass range. The resulting limits on sin^{2}2θ_{14} are improved by approx imately a factor of 2 over previous results and constitute the most stringent constraints to date in the |Δm_{41}^{2}|≲0.2  eV^{2} region.