A high precision calibration of the nonlinearity in the energy response of the Daya Bay Reactor Neutrino Experiment's antineutrino detectors is presented in detail. The energy nonlinearity originates from the particle-dependent light yield of the scintillator and charge-dependent electronics response. The nonlinearity model is constrained by γ calibration points from deployed and naturally occurring radioactive sources, the β spectrum from 12B decays, and a direct measurement of the electronics nonlinearity with a new flash analog-to-digital converter readout system. Less than 0.5% uncertainty in the energy nonlinearity calibration is achieved for positrons of kinetic energies greater than 1 MeV.

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.

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.

A search for a time-varying νe signal was performed with 621 days of data acquired by the Daya Bay Reactor Neutrino Experiment over 704 calendar days. The time spectrum of the measured νe flux normalized to its prediction was analyzed with a Lomb-Scargle periodogram, which yielded no significant signal for periods ranging from 2 hours to nearly 2 years. The normalized time spectrum was also fit for a sidereal modulation under the Standard Model Extension (SME) framework to search for Lorentz and CPT violation (LV-CPTV). Limits were obtained for all six flavor pairs ēμ,ēτ, μτ, ēē,μμ and ττ by fitting them one at a time, constituting the first experimental constraints on the latter three. Daya Bay's high statistics and unique layout of multiple directions from three pairs of reactors to three experimental halls allowed the simultaneous constraint of individual SME LV-CPTV coefficients without assuming others contribute negligibly, a first for a neutrino experiment.

The prediction of reactor antineutrino spectra will play a crucial role as reactor experiments enter the precision era. The positron energy spectrum of 3.5 million antineutrino inverse beta decay reactions observed by the Daya Bay experiment, in combination with the fission rates of fissile isotopes in the reactor, is used to extract the positron energy spectra resulting from the fission of specific isotopes. This information can be used to produce a precise, data-based prediction of the antineutrino energy spectrum in other reactor antineutrino experiments with different fission fractions than Daya Bay. The positron energy spectra are unfolded to obtain the antineutrino energy spectra by removing the contribution from detector response with the Wiener-SVD unfolding method. Consistent results are obtained with other unfolding methods. A technique to construct a data-based prediction of the reactor antineutrino energy spectrum is proposed and investigated. Given the reactor fission fractions, the technique can predict the energy spectrum to a 2% precision. In addition, we illustrate how to perform a rigorous comparison between the unfolded antineutrino spectrum and a theoretical model prediction that avoids the input model bias of the unfolding method.

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.

The Daya Bay Collaboration responds to comments posted by S. Parke and R. Zukanovich Funchal regarding our use of Delta(m^2_ee).

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.

The prediction of reactor antineutrino spectra will play a crucial role as reactor experiments enter the precision era. The positron energy spectrum of 3.5 million antineutrino inverse beta decay reactions observed by the Daya Bay experiment, in combination with the fission rates of fissile isotopes in the reactor, is used to extract the positron energy spectra resulting from the fission of specific isotopes. This information can be used to produce a precise, data-based prediction of the antineutrino energy spectrum in other reactor antineutrino experiments with different fission fractions than Daya Bay. The positron energy spectra are unfolded to obtain the antineutrino energy spectra by removing the contribution from detector response with the Wiener-SVD unfolding method. Consistent results are obtained with other unfolding methods. A technique to construct a data-based prediction of the reactor antineutrino energy spectrum is proposed and investigated. Given the reactor fission fractions, the technique can predict the energy spectrum to a 2% precision. In addition, we illustrate how to perform a rigorous comparison between the unfolded antineutrino spectrum and a theoretical model prediction that avoids the input model bias of the unfolding method.

The establishment of a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using the data collected from December 2011 to August 2017, a search was performed for electron-Antineutrino signals that coincided with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of 10, 500, and 1000 s relative to the occurrence of the GW events and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-Antineutrino candidates were consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) of the electron-Antineutrino fluence of (1.13-2.44)×1011 cm-2 at 5 MeV to 8.0×107 cm-2 at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be (5.4-7.0)×109 cm-2 for the three time windows.