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\pm 0.86)\times 10^{-5}$, $(10.22\pm 0.87)\times 10^{-5}$, and $(17.03\pm 1.22)\times 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 \pm 0.03$ for the dependence of the neutron yield on muon energy.

The Daya Bay Experiment consists of eight identically designed detectors located in three underground experimental halls named as EH1, EH2, EH3, with 250, 265 and 860 meters of water equivalent vertical overburden, respectively. Cosmic muon events have been recorded over a two-year period. The underground muon rate is observed to be positively correlated with the effective atmospheric temperature and to follow a seasonal modulation pattern. The correlation coefficient $\alpha$, describing how a variation in the muon rate relates to a variation in the effective atmospheric temperature, is found to be $\alpha_{\text{EH1}} = 0.362\pm0.031$, $\alpha_{\text{EH2}} = 0.433\pm0.038$ and $\alpha_{\text{EH3}} = 0.641\pm0.057$ for each experimental hall.

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$_{\textrm{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, $\bar{\sigma}_f$, of $(5.90 \pm 0.13) \times 10^{-43}$ cm$^2$/fission and a fuel-dependent variation in the IBD yield, $d\sigma_f/dF_{239}$, of $(-1.86 \pm 0.18) \times 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 was 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$\sigma$. This discrepancy indicates that an overall deficit in 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 \pm 0.17)$ and $(4.27 \pm 0.26) \times 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 yield suggests that this isotope may be the primary contributor to the reactor antineutrino anomaly.

A new measurement of the reactor antineutrino flux and energy spectrum by the Daya Bay reactor neutrino experiment is reported. The antineutrinos were generated by six 2.9~GW$_{\mathrm{th}}$ nuclear reactors and detected by eight antineutrino detectors deployed in two near (560~m and 600~m flux-weighted baselines) and one far (1640~m flux-weighted baseline) underground experimental halls. With 621 days of data, more than 1.2 million inverse beta decay (IBD) candidates were detected. The IBD yield in the eight detectors was measured, and the ratio of measured to predicted flux was found to be $0.946\pm0.020$ ($0.992\pm0.021$) for the Huber+Mueller (ILL+Vogel) model. A 2.9~$\sigma$ deviation was found in the measured IBD positron energy spectrum compared to the predictions. In particular, an excess of events in the region of 4-6~MeV was found in the measured spectrum, with a local significance of 4.4~$\sigma$. A reactor antineutrino spectrum weighted by the IBD cross section is extracted for model-independent predictions.

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 LSND and MiniBooNE experiments in a minimally extended four-neutrino flavor framework. Stringent limits on $\sin^2 2\theta_{\mu e}$ are set over 6 orders of magnitude in the sterile mass-squared splitting $\Delta m^2_{41}$. The sterile-neutrino mixing phase space allowed by the LSND and MiniBooNE experiments is excluded for $\Delta m^2_{41} < 0.8$ eV$^2$ at 95% CL$_s$.

A measurement of electron antineutrino oscillation by the Daya Bay Reactor Neutrino Experiment is described in detail. Six 2.9-GW$_{\rm th}$ nuclear power reactors of the Daya Bay and Ling Ao nuclear power facilities served as intense sources of $\overline{\nu}_{e}$'s. Comparison of the $\overline{\nu}_{e}$ rate and energy spectrum measured by antineutrino detectors far from the nuclear reactors ($\sim$1500-1950 m) relative to detectors near the reactors ($\sim$350-600 m) allowed a precise measurement of $\overline{\nu}_{e}$ disappearance. More than 2.5 million $\overline{\nu}_{e}$ inverse beta decay interactions were observed, based on the combination of 217 days of operation of six antineutrino detectors (Dec. 2011--Jul. 2012) with a subsequent 1013 days using the complete configuration of eight detectors (Oct. 2012--Jul. 2015). The $\overline{\nu}_{e}$ rate observed at the far detectors relative to the near detectors showed a significant deficit, $R=0.949 \pm 0.002(\mathrm{stat.}) \pm 0.002(\mathrm{syst.})$. The energy dependence of $\overline{\nu}_{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 $\sin^22\theta_{13}=0.0841 \pm 0.0027(\mathrm{stat.}) \pm 0.0019(\mathrm{syst.})$ and the effective neutrino mass-squared difference of $\left|{\Delta}m^2_{\mathrm{ee}}\right|=(2.50 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$. Analysis using the exact three-flavor probability found ${\Delta}m^2_{32}=(2.45 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$ assuming the normal neutrino mass hierarchy and ${\Delta}m^2_{32}=(-2.56 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$ for the inverted hierarchy.

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\times10^{-4} \lesssim |\Delta m^{2}_{41}| \lesssim 0.3$ eV$^{2}$ mass range. The resulting limits on $\sin^{2}2\theta_{14}$ are improved by approximately a factor of 2 over previous results and constitute the most stringent constraints to date in the $|\Delta m^{2}_{41}| \lesssim 0.2$ eV$^{2}$ region.

The disappearance of reactor $\bar{\nu}_e$ observed by the Daya Bay experiment is examined in the framework of a model in which the neutrino is described by a wave packet with a relative intrinsic momentum dispersion $\sigma_\text{rel}$. Three pairs of nuclear reactors and eight antineutrino detectors, each with good energy resolution, distributed among three experimental halls, supply a high-statistics sample of $\bar{\nu}_e$ acquired at nine different baselines. This provides a unique platform to test the effects which arise from the wave packet treatment of neutrino oscillation. The modified survival probability formula was used to fit Daya Bay data, providing the first experimental limits: $2.38 \cdot 10^{-17} < \sigma_{\rm rel} < 0.23$. Treating the dimensions of the reactor cores and detectors as constraints, the limits are improved: $10^{-14} \lesssim \sigma_{\rm rel} < 0.23$, and an upper limit of $\sigma_{\rm rel} <0.20$ is obtained. All limits correspond to a 95\% C.L. Furthermore, the effect due to the wave packet nature of neutrino oscillation is found to be insignificant for reactor antineutrinos detected by the Daya Bay experiment thus ensuring an unbiased measurement of the oscillation parameters $\sin^22\theta_{13}$ and $\Delta m^2_{32}$ within the plane wave model.

This article reports an improved independent measurement of neutrino mixing angle $\theta_{13}$ at the Daya Bay Reactor Neutrino Experiment. Electron antineutrinos were identified by inverse $\beta$-decays with the emitted neutron captured by hydrogen, yielding a data-set with principally distinct uncertainties from that with neutrons captured by gadolinium. With the final two of eight antineutrino detectors installed, this study used 621 days of data including the previously reported 217-day data set with six detectors. The dominant statistical uncertainty was reduced by 49%. Intensive studies of the cosmogenic muon-induced $^9$Li and fast neutron backgrounds and the neutron-capture energy selection efficiency, resulted in a reduction of the systematic uncertainty by 26%. The deficit in the detected number of antineutrinos at the far detectors relative to the expected number based on the near detectors yielded $\sin^22\theta_{13} = 0.071 \pm 0.011$ in the three-neutrino-oscillation framework. The combination of this result with the gadolinium-capture result is also reported.

The Daya Bay experiment was the first to report simultaneous measurements of reactor antineutrinos at multiple baselines leading to the discovery of $\bar{\nu}_e$ oscillations over km-baselines. Subsequent data has provided the world's most precise measurement of $\rm{sin}^22\theta_{13}$ and the effective mass splitting $\Delta m_{ee}^2$. The experiment is located in Daya Bay, China where the cluster of six nuclear reactors is among the world's most prolific sources of electron antineutrinos. Multiple antineutrino detectors are deployed in three underground water pools at different distances from the reactor cores to search for deviations in the antineutrino rate and energy spectrum due to neutrino mixing. Instrumented with photomultiplier tubes (PMTs), the water pools serve as shielding against natural radioactivity from the surrounding rock and provide efficient muon tagging. Arrays of resistive plate chambers over the top of each pool provide additional muon detection. The antineutrino detectors were specifically designed for measurements of the antineutrino flux with minimal systematic uncertainty. Relative detector efficiencies between the near and far detectors are known to better than 0.2%. With the unblinding of the final two detectors' baselines and target masses, a complete description and comparison of the eight antineutrino detectors can now be presented. This paper describes the Daya Bay detector systems, consisting of eight antineutrino detectors in three instrumented water pools in three underground halls, and their operation through the first year of eight detector data-taking.

The Jiangmen Underground Neutrino Observatory (JUNO) is proposed to determine the neutrino mass hierarchy using an underground liquid scintillator detector. It is located 53 km away from both Yangjiang and Taishan Nuclear Power Plants in Guangdong, China. The experimental hall, spanning more than 50 meters, is under a granite mountain of over 700 m overburden. Within six years of running, the detection of reactor antineutrinos can resolve the neutrino mass hierarchy at a confidence level of 3-4$\sigma$, and determine neutrino oscillation parameters $\sin^2\theta_{12}$, $\Delta m^2_{21}$, and $|\Delta m^2_{ee}|$ to an accuracy of better than 1%. The JUNO detector can be also used to study terrestrial and extra-terrestrial neutrinos and new physics beyond the Standard Model. The central detector contains 20,000 tons liquid scintillator with an acrylic sphere of 35 m in diameter. $\sim$17,000 508-mm diameter PMTs with high quantum efficiency provide $\sim$75% optical coverage. The current choice of the liquid scintillator is: linear alkyl benzene (LAB) as the solvent, plus PPO as the scintillation fluor and a wavelength-shifter (Bis-MSB). The number of detected photoelectrons per MeV is larger than 1,100 and the energy resolution is expected to be 3% at 1 MeV. The calibration system is designed to deploy multiple sources to cover the entire energy range of reactor antineutrinos, and to achieve a full-volume position coverage inside the detector. The veto system is used for muon detection, muon induced background study and reduction. It consists of a Water Cherenkov detector and a Top Tracker system. The readout system, the detector control system and the offline system insure efficient and stable data acquisition and processing.

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$\times$10$^5$ GW$_{\rm 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 two 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\theta_{13}$ and $|\Delta m^2_{ee}|$ were halved as a result of these improvements. Analysis of the relative antineutrino rates and energy spectra between detectors gave $\sin^{2}2\theta_{13} = 0.084\pm0.005$ and $|\Delta m^{2}_{ee}|= (2.42\pm0.11) \times 10^{-3}$ eV$^2$ in the three-neutrino framework.

This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 2.9~GW$_{th}$ nuclear reactors with six detectors deployed in two near (effective baselines 512~m and 561~m) and one far (1,579~m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296,721 and 41,589 inverse beta decay (IBD) candidates were detected in the near and far halls, respectively. The measured IBD yield is (1.55 $\pm$ 0.04) $\times$ 10$^{-18}$~cm$^2$/GW/day or (5.92 $\pm$ 0.14) $\times$ 10$^{-43}$~cm$^2$/fission. This flux measurement is consistent with previous short-baseline reactor antineutrino experiments and is $0.946\pm0.022$ ($0.991\pm0.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$\sigma$ over the full energy range with a local significance of up to $\sim$4$\sigma$ between 4-6 MeV. A reactor antineutrino spectrum of IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions.

A search for light sterile neutrino mixing was performed with the first 217 days of data from the Daya Bay Reactor Antineutrino Experiment. The experiment's unique configuration of multiple baselines from six 2.9~GW$_{\rm th}$ nuclear reactors to six antineutrino detectors deployed in two near (effective baselines 512~m and 561~m) and one far (1579~m) underground experimental halls makes it possible to test for oscillations to a fourth (sterile) neutrino in the $10^{\rm -3}~{\rm eV}^{2} < |\Delta m_{41}^{2}| < 0.3~{\rm eV}^{2}$ range. The relative spectral distortion due to electron antineutrino disappearance was found to be consistent with that of the three-flavor oscillation model. The derived limits on $\sin^22\theta_{14}$ cover the $10^{-3}~{\rm eV}^{2} \lesssim |\Delta m^{2}_{41}| \lesssim 0.1~{\rm eV}^{2}$ region, which was largely unexplored.

A new measurement of the $\theta_{13}$ mixing angle has been obtained at the Daya Bay Reactor Neutrino Experiment via the detection of inverse beta decays tagged by neutron capture on hydrogen. The antineutrino events for hydrogen capture are distinct from those for gadolinium capture with largely different systematic uncertainties, allowing a determination independent of the gadolinium-capture result and an improvement on the precision of $\theta_{13}$ measurement. With a 217-day antineutrino data set obtained with six antineutrino detectors and from six 2.9 GW$_{th}$ reactors, the rate deficit observed at the far hall is interpreted as $\sin^22\theta_{13}=0.083\pm0.018$ in the three-flavor oscillation model. When combined with the gadolinium-capture result from Daya Bay, we obtain $\sin^22\theta_{13}=0.089\pm0.008$ as the final result for the six-antineutrino-detector configuration of the Daya Bay experiment.

A measurement of the energy dependence of antineutrino disappearance at the Daya Bay Reactor Neutrino Experiment is reported. Electron antineutrinos ($\overline{\nu}_{e}$) from six $2.9$ GW$_{\rm th}$ reactors were detected with six detectors deployed in two near (effective baselines 512 m and 561 m) and one far (1579 m) underground experimental halls. Using 217 days of data, 41589 (203809 and 92912) antineutrino candidates were detected in the far hall (near halls). An improved measurement of the oscillation amplitude $\sin^{2}2\theta_{13} = 0.090^{+0.008}_{-0.009} $ and the first direct measurement of the $\overline{\nu}_{e}$ mass-squared difference $|\Delta m^{2}_{ee}|= (2.59_{-0.20}^{+0.19}) \times 10^{-3}\ {\rm eV}^2 $ is obtained using the observed $\overline{\nu}_{e}$ rates and energy spectra in a three-neutrino framework. This value of $|\Delta m^{2}_{ee}|$ is consistent with $|\Delta m^{2}_{\mu\mu}|$ measured by muon neutrino disappearance, supporting the three-flavor oscillation model.