We have measured the muon flux and production rate of muon-induced neutrons at a depth of 611 m water equivalent. Our apparatus comprises three layers of crossed plastic scintillator hodoscopes for tracking the incident cosmic-ray muons and 760 L of a gadolinium-doped liquid scintillator for producing and detecting neutrons. The vertical muon intensity was measured to be Iμ=(5.7±0.6)×10-6 cm-2 s-1 sr-1. The yield of muon-induced neutrons in the liquid scintillator was determined to be Yn=(1.19±0.08(stat)±0.21(syst))×10-4 neutrons/(μ·g·cm-2). A fit to the recently measured neutron yields at different depths gave a mean muon energy dependence of 0.76±0.03 for liquid-scintillator targets.

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.

We report a new measurement of electron antineutrino disappearance using the fully constructed Daya Bay Reactor Neutrino Experiment. The final two of eight antineutrino detectors were installed in the summer of 2012. Including the 404 days of data collected from October 2012 to November 2013 resulted in a total exposure of 6.9×10^{5} GW_{th}  ton days, a 3.6 times increase over our previous results. Improvements in energy calibration limited variations between detectors to 0.2%. Removal of six ^{241}Am-^{13}C radioactive calibration sources reduced the background by a factor of 2 for the detectors in the experimental hall furthest from the reactors. Direct prediction of the antineutrino signal in the far detectors based on the measurements in the near detectors explicitly minimized the dependence of the measurement on models of reactor antineutrino emission. The uncertainties in our estimates of sin^{2}2θ_{13} and |Δm_{ee}^{2}| were halved as a result of these improvements. An analysis of the relative antineutrino rates and energy spectra between detectors gave sin^{2}2θ_{13}=0.084±0.005 and |Δm_{ee}^{2}|=(2.42±0.11)×10^{-3}  eV^{2} in the three-neutrino framework.

A measurement of electron antineutrino oscillation by the Daya Bay Reactor Neutrino Experiment is described in detail. Six 2.9-GWth nuclear power reactors of the Daya Bay and Ling Ao nuclear power facilities served as intense sources of νe's. Comparison of the νe rate and energy spectrum measured by antineutrino detectors far from the nuclear reactors (∼1500-1950 m) relative to detectors near the reactors (∼350-600 m) allowed a precise measurement of νe disappearance. More than 2.5 million νe inverse beta-decay interactions were observed, based on the combination of 217 days of operation of six antineutrino detectors (December, 2011-July, 2012) with a subsequent 1013 days using the complete configuration of eight detectors (October, 2012-July, 2015). The νe rate observed at the far detectors relative to the near detectors showed a significant deficit, R=0.949±0.002(stat)±0.002(syst). The energy dependence of νe disappearance showed the distinct variation predicted by neutrino oscillation. Analysis using an approximation for the three-flavor oscillation probability yielded the flavor-mixing angle sin22θ13=0.0841±0.0027(stat)±0.0019(syst) and the effective neutrino mass-squared difference of |Δmee2|=(2.50±0.06(stat)±0.06(syst))×10-3 eV2. Analysis using the exact three-flavor probability found Δm322=(2.45±0.06(stat)±0.06(syst))×10-3 eV2 assuming the normal neutrino mass hierarchy and Δm322=(-2.56±0.06(stat)±0.06(syst))×10-3 eV2 for the inverted hierarchy.

This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 2.9 GWth nuclear reactors with six detectors deployed in two near (effective baselines 512 and 561 m) and one far (1579 m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296 721 and 41 589 inverse β decay (IBD) candidates were detected in the near and far halls, respectively. The measured IBD yield is (1.55±0.04) ×10(-18)  cm(2) GW(-1) day(-1) or (5.92±0.14) ×10(-43)  cm(2) fission(-1). This flux measurement is consistent with previous short-baseline reactor antineutrino experiments and is 0.946±0.022 (0.991±0.023) relative to the flux predicted with the Huber-Mueller (ILL-Vogel) fissile antineutrino model. The measured IBD positron energy spectrum deviates from both spectral predictions by more than 2σ over the full energy range with a local significance of up to ∼4σ between 4-6 MeV. A reactor antineutrino spectrum of IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions.

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.

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.

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.