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.

Measurement of the ν_{e}-Nucleus Charged-Current Double-Differential Cross Section at ⟨E_{ν}⟩=2.4 GeV Using NOvA

The inclusive electron neutrino charged-current cross section is measured in the NOvA near detector using 8.02×10^{20} protons-on-target in the NuMI beam. The sample of GeV electron neutrino interactions is the largest analyzed to date and is limited by ≃17% systematic rather than the ≃7.4% statistical uncertainties. The double-differential cross section in final-state electron energy and angle is presented for the first time, together with the single-differential dependence on Q^{2} (squared four-momentum transfer) and energy, in the range 1  GeV≤E_{ν}<6  GeV. Detailed comparisons are made to the predictions of the GENIE, GiBUU, NEUT, and NuWro neutrino event generators. The data do not strongly favor a model over the others consistently across all three cross sections measured, though some models have especially good or poor agreement in the single differential cross section vs Q^{2}.

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.

Scintillation light detection in the 6-m drift-length ProtoDUNE Dual Phase liquid argon TPC

DUNE is a dual-site experiment for long-baseline neutrino oscillation studies, neutrino astrophysics and nucleon decay searches. ProtoDUNE Dual Phase (DP) is a 6  ×  6  ×  6 m 3 liquid argon time-projection-chamber (LArTPC) that recorded cosmic-muon data at the CERN Neutrino Platform in 2019-2020 as a prototype of the DUNE Far Detector. Charged particles propagating through the LArTPC produce ionization and scintillation light. The scintillation light signal in these detectors can provide the trigger for non-beam events. In addition, it adds precise timing capabilities and improves the calorimetry measurements. In ProtoDUNE-DP, scintillation and electroluminescence light produced by cosmic muons in the LArTPC is collected by photomultiplier tubes placed up to 7 m away from the ionizing track. In this paper, the ProtoDUNE-DP photon detection system performance is evaluated with a particular focus on the different wavelength shifters, such as PEN and TPB, and the use of Xe-doped LAr, considering its future use in giant LArTPCs. The scintillation light production and propagation processes are analyzed and a comparison of simulation to data is performed, improving understanding of the liquid argon properties.

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.

Search for Active-Sterile Antineutrino Mixing Using Neutral-Current Interactions with the NOvA Experiment

This Letter reports results from the first long-baseline search for sterile antineutrinos mixing in an accelerator-based antineutrino-dominated beam. The rate of neutral-current interactions in the two NOvA detectors, at distances of 1 and 810 km from the beam source, is analyzed using an exposure of 12.51×10^{20} protons-on-target from the NuMI beam at Fermilab running in antineutrino mode. A total of 121 of neutral-current candidates are observed at the far detector, compared to a prediction of 122±11(stat.)±15(syst.) assuming mixing only between three active flavors. No evidence for ν[over ¯]_{μ}→ν[over ¯]_{s} oscillation is observed. Interpreting this result within a 3+1 model, constraints are placed on the mixing angles θ_{24}<25° and θ_{34}<32° at the 90% C.L. for 0.05  eV^{2}≤Δm_{41}^{2}≤0.5  eV^{2}, the range of mass splittings that produces no significant oscillations at the near detector. These are the first 3+1 confidence limits set using long-baseline accelerator antineutrinos.

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}.

First measurement of neutrino oscillation parameters using neutrinos and antineutrinos by NOvA

The NOvA experiment has seen a 4.4σ signal of ν[over ¯]_{e} appearance in a 2 GeV ν[over ¯]_{μ} beam at a distance of 810 km. Using 12.33×10^{20} protons on target delivered to the Fermilab NuMI neutrino beamline, the experiment recorded 27 ν[over ¯]_{μ}→ν[over ¯]_{e} candidates with a background of 10.3 and 102 ν[over ¯]_{μ}→ν[over ¯]_{μ} candidates. This new antineutrino data are combined with neutrino data to measure the parameters |Δm_{32}^{2}|=2.48_{-0.06}^{+0.11}×10^{-3}  eV^{2}/c^{4} and sin^{2}θ_{23} in the ranges from (0.53-0.60) and (0.45-0.48) in the normal neutrino mass hierarchy. The data exclude most values near δ_{CP}=π/2 for the inverted mass hierarchy by more than 3σ and favor the normal neutrino mass hierarchy by 1.9σ and θ_{23} values in the upper octant by 1.6σ.

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.

Measurement of the Neutrino Mixing Angle θ_{23} in NOvA

This Letter reports new results on muon neutrino disappearance from NOvA, using a 14 kton detector equivalent exposure of 6.05×10^{20} protons on target from the NuMI beam at the Fermi National Accelerator Laboratory. The measurement probes the muon-tau symmetry hypothesis that requires maximal θ_{23} mixing (θ_{23}=π/4). Assuming the normal mass hierarchy, we find Δm_{32}^{2}=(2.67±0.11)×10^{-3}  eV^{2} and sin^{2}θ_{23} at the two statistically degenerate values 0.404_{-0.022}^{+0.030} and 0.624_{-0.030}^{+0.022}, both at the 68% confidence level. Our data disfavor the maximal mixing scenario with 2.6σ significance.

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.

First Measurement of Electron Neutrino Appearance in NOvA

We report results from the first search for ν_{μ}→ν_{e} transitions by the NOvA experiment. In an exposure equivalent to 2.74×10^{20} protons on target in the upgraded NuMI beam at Fermilab, we observe 6 events in the Far Detector, compared to a background expectation of 0.99±0.11(syst) events based on the Near Detector measurement. A secondary analysis observes 11 events with a background of 1.07±0.14(syst). The 3.3σ excess of events observed in the primary analysis disfavors 0.1π<δ_{CP}<0.5π in the inverted mass hierarchy at the 90% C.L.