The full data set of the Daya Bay reactor neutrino experiment is used to probe the effect of the charged current non-standard interactions (CC-NSI) on neutrino oscillation experiments. Two different approaches are applied and constraints on the corresponding CC-NSI parameters are obtained with the neutrino flux taken from the Huber-Mueller model with a $5\%$ uncertainty. For the quantum mechanics-based approach (QM-NSI), the constraints on the CC-NSI parameters $\epsilon_{e\alpha}$ and $\epsilon_{e\alpha}^{s}$ are extracted with and without the assumption that the effects of the new physics are the same in the production and detection processes, respectively. The approach based on the weak effective field theory (WEFT-NSI) deals with four types of CC-NSI represented by the parameters $[\varepsilon_{X}]_{e\alpha}$. For both approaches, the results for the CC-NSI parameters are shown for cases with various fixed values of the CC-NSI and the Dirac CP-violating phases, and when they are allowed to vary freely. We find that constraints on the QM-NSI parameters $\epsilon_{e\alpha}$ and $\epsilon_{e\alpha}^{s}$ from the Daya Bay experiment alone can reach the order $\mathcal{O}(0.01)$ for the former and $\mathcal{O}(0.1)$ for the latter, while for WEFT-NSI parameters $[\varepsilon_{X}]_{e\alpha}$, we obtain $\mathcal{O}(0.1)$ for both cases.

Daya Bay presents the first measurement of cosmogenic He-8 isotope production in liquid scintillator, using an innovative method for identifying cascade decays of He-8 and its child isotope, Li-8. We also measure the production yield of Li-9 isotopes using well-established methodology. The results, in units of 10(-8)mu-1g(-1) cm(2), are 0.307 +/- 0.042, 0.341 +/- 0.040, and 0.546 +/- 0.076 for He-8, and 6.73 +/- 0.73, 6.75 +/- 0.70, and 13.74 +/- 0.82 for 9Li at average muon energies of 63.9 GeV, 64.7 GeV, and 143.0 GeV, respectively. The measured production rate of He-8 isotopes is more than an order of magnitude lower than any other measurement of cosmogenic isotope production. It replaces the results of previous attempts to determine the ratio of He-8 to Li-9 production that yielded a wide range of limits from 0% to 30%. The results provide future liquid-scintillator-based experiments with improved ability to predict cosmogenic backgrounds.

This Letter reports the first measurement of the oscillation amplitude and frequency of reactor antineutrinos at Daya Bay via neutron capture on hydrogen using 1958 days of data. With over 3.6 million signal candidates, an optimized candidate selection, improved treatment of backgrounds and efficiencies, refined energy calibration, and an energy response model for the capture-on-hydrogen sensitive region, the relative ν[over ¯]_{e} rates and energy spectra variation among the near and far detectors gives sin^{2}2θ_{13}=0.0759_{-0.0049}^{+0.0050} and Δm_{32}^{2}=(2.72_{-0.15}^{+0.14})×10^{-3}  eV^{2} assuming the normal neutrino mass ordering, and Δm_{32}^{2}=(-2.83_{-0.14}^{+0.15})×10^{-3}  eV^{2} for the inverted neutrino mass ordering. This estimate of sin^{2}2θ_{13} is consistent with and essentially independent from the one obtained using the capture-on-gadolinium sample at Daya Bay. The combination of these two results yields sin^{2}2θ_{13}=0.0833±0.0022, which represents an 8% relative improvement in precision regarding the Daya Bay full 3158-day capture-on-gadolinium result.

This Letter presents results of a search for the mixing of a sub-eV sterile neutrino with three active neutrinos based on the full data sample of the Daya Bay Reactor Neutrino Experiment, collected during 3158 days of detector operation, which contains 5.55×10^{6} reactor ν[over ¯]_{e} candidates identified as inverse beta-decay interactions followed by neutron capture on gadolinium. The analysis benefits from a doubling of the statistics of our previous result and from improvements of several important systematic uncertainties. No significant oscillation due to mixing of a sub-eV sterile neutrino with active neutrinos was found. Exclusion limits are set by both Feldman-Cousins and CLs methods. Light sterile neutrino mixing with sin^{2}2θ_{14}≳0.01 can be excluded at 95% confidence level in the region of 0.01  eV^{2}≲|Δm_{41}^{2}|≲0.1  eV^{2}. This result represents the world-leading constraints in the region of 2×10^{-4}  eV^{2}≲|Δm_{41}^{2}|≲0.2  eV^{2}.

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.

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.

A joint determination of the reactor antineutrino spectra resulting from the fission of ^235U and ^239Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of ^235U 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 ^235U and ^239Pu isotopes and improves the uncertainty of the ^235U spectral shape to about 3%. The ^235U and ^239Pu 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 ^235U and ^239Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

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 Day a 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 Day a 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)×10 11 cm −2 at 5 MeV to 8.0×10 7 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)×10 9 cm −2 for the three time windows.

To maximize the light yield of the liquid scintillator (LS) for the Jiangmen Underground Neutrino Observatory (JUNO), a 20 t LS sample was produced in a pilot plant at Daya Bay. The optical properties of the new LS in various compositions were studied by replacing the gadolinium-loaded LS in one antineutrino detector. The concentrations of the fluor, PPO, and the wavelength shifter, bis-MSB, were increased in 12 steps from 0.5 g/L and <0.01 mg/L to 4 g/L and 13 mg/L, respectively. The numbers of total detected photoelectrons suggest that, with the optically purified solvent, the bis-MSB concentration does not need to be more than 4 mg/L. To bridge the one order of magnitude in the detector size difference between Daya Bay and JUNO, the Daya Bay data were used to tune the parameters of a newly developed optical model. Then, the model and tuned parameters were used in the JUNO simulation. This enabled to determine the optimal composition for the JUNO LS: purified solvent LAB with 2.5 g/L PPO, and 1 to 4 mg/L bis-MSB.

The DUNE IDR describes the proposed physics program and technical designs of the DUNE Far Detector modules in preparation for the full TDR to be published in 2019. It is intended as an intermediate milestone on the path to a full TDR, justifying the technical choices that flow down from the high-level physics goals through requirements at all levels of the Project. These design choices will enable the DUNE experiment to make the ground-breaking discoveries that will help to answer fundamental physics questions. Volume 1 contains an executive summary that describes the general aims of this document. The remainder of this first volume provides a more detailed description of the DUNE physics program that drives the choice of detector technologies. It also includes concise outlines of two overarching systems that have not yet evolved to consortium structures: computing and calibration. Volumes 2 and 3 of this IDR describe, for the single-phase and dual-phase technologies, respectively, each detector module's subsystems, the technical coordination required for its design, construction, installation, and integration, and its organizational structure.

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 $\theta_{\mu e}$ mixing angle are derived that constitute the most stringent limits to date over five orders of magnitude in the sterile mass-squared splitting $\Delta m^2_{41}$, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL$_s$ for $\Delta m^2_{41}<5\,$eV$^2$.Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL$_s$ for $\Delta m^2_{41}$ $<$ 1.2 eV$^2$.

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 $\gamma$ calibration points from deployed and naturally occurring radioactive sources, the $\beta$ spectrum from $^{12}$B 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.

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

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.

A search for a time-varying (nu) over bar (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 (nu) over bar (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 (e) over bar(mu) over bar, (e) over bar(tau) over bar, (mu) over bar(tau) over bar, (e) over bar(e) over bar, (mu) over bar(mu) over bar and (tau) over bar(tau) over bar 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.

Neutrons produced by cosmic ray muons are an important background for underground experiments studying neutrino oscillations, neutrinoless double beta decay, dark matter, and other rare-event signals. A measurement of the neutron yield in the three different experimental halls of the Daya Bay Reactor Neutrino Experiment at varying depth is reported. The neutron yield in Daya Bay's liquid scintillator is measured to be Y-n = (10.26 +/- 0.86) x 10(-5), (10.22 +/- 0.87) x 10(-5), and (17.03 +/- 1.22) x 10(-5) mu(-1)g(-1)cm(2) at depths of 250, 265, and 860 meters-water-equivalent. These results are compared to other measurements and the simulated neutron yield in FLUKA and GEANT4. A global fit including the Daya Bay measurements yields a power law coefficient of 0.77 +/- 0.03 for the dependence of the neutron yield on muon energy.

The DUNE IDR describes the proposed physics program and technical designs of the DUNE far detector modules in preparation for the full TDR to be published in 2019. It is intended as an intermediate milestone on the path to a full TDR, justifying the technical choices that flow down from the high-level physics goals through requirements at all levels of the Project. These design choices will enable the DUNE experiment to make the ground-breaking discoveries that will help to answer fundamental physics questions. Volume 3 describes the dual-phase module's subsystems, the technical coordination required for its design, construction, installation, and integration, and its organizational structure.

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.