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

Spent nuclear fuel (SNF) antineutrino flux is an important source of
uncertainties for a reactor neutrino flux prediction. However, if one want to
determine the contribution of spent fuel, many data are needed, such as the
amount of spent fuel in the pool, the time after discharged from the reactor
core, the burnup of each assembly, and the antineutrino spectrum of the
isotopes in the spend fuel. A method to calculate the contribution of SNF is
proposed in this study. In this method, reactor simulation code verified by
experiment have been used to simulate the fuel depletion by taking into account
more than 2000 isotopes and fission products, the quantity of SNF in each six
spend fuel pool, and the antineutrino spectrum of SNF varying with time after
SNF discharged from core. Results show that the contribution of SNF to the
total antineutrino flux is about 0.26%~0.34%, and the shutdown impact is about
20%. The SNF spectrum would distort the softer part of antineutrino spectra,
and the maximum contribution from SNF is about 3.0%, but there is 18\%
difference between line evaluate method and under evaluate method. In addition,
nonequilibrium effects are also discussed, and the results are compatible with
theirs considering the uncertainties.

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 fluxweighted
baselines) and one far (1640~m fluxweighted 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 46~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 modelindependent predictions.

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

Uncertainties of fission fraction is an important uncertainty source for the
antineutrino flux prediction in a reactor antineutrino experiment. A new
MCbased method of evaluating the covariance coefficients between isotopes was
proposed. It was found that the covariance coefficients will varying with
reactor burnup and which may change from positive to negative because of
fissioning balance effect, for example, the covariance coefficient between
$^{235}$U and $^{239}$Pu changes from 0.15 to 0.13. Using the equation between
fission fraction and atomic density, the consistent of uncertainty of fission
fraction and the covariance matrix were obtained. The antineutrino flux
uncertainty is 0.55\% which does not vary with reactor burnup, and the new
value is about 8.3\% smaller.

Reactor antineutrino experiment are used to study neutrino oscillation,
search for signatures of nonstandard neutrino interaction, and monitor reactor
operation for safeguard application. Reactor simulation is an important source
of uncertainties for a reactor neutrino experiment. Commercial code is used for
reactor simulation to evaluate fission fraction in Daya Bay neutrino
experiment, but the source code doesn't open to our researcher results from
commercial secret. In this study, The open source code DRAGON was improved to
calculate the fission rates of the four most important isotopes in fissions,
$^{235}$U,$^{238}$U,$^{239}$Pu and $^{241}$Pu, and then was validated for PWRs
using the Takahama3 benchmark. The fission fraction results are consistent
with those of MIT's results. Then, fission fraction of Daya Bay reactor core
was calculated by using improved DRAGON code, and the fission fraction
calculated by DRAGON agreed well with these calculated by SCIENCE. The average
deviation less than 5\% for all the four isotopes. The correlation coefficient
matrix between $^{235}$U,$^{238}$U,$^{239}$Pu and $^{241}$Pu were also studied
using DRAGON, and then the uncertainty of the antineutrino flux by the fission
fraction was calculated by using the correlation coefficient matrix. The
uncertainty of the antineutrino flux by the fission fraction simulation is
0.6\% per core for Daya Bay antineutrino experiment. The uncertainties source
of fission fraction calculation need further to be studied in the future.

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}$ masssquared 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
threeneutrino framework.
This value of $\Delta m^{2}_{ee}$ is consistent with $\Delta
m^{2}_{\mu\mu}$ measured by muon neutrino disappearance, supporting the
threeflavor oscillation model.