
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.

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.

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.

The glueballtovacuum matrix elements of local gluonic operators in scalar,
tensor, and pseudoscalar channels are investigated numerically on several
anisotropic lattices with the spatial lattice spacing ranging from 0.1fm 
0.2fm. These matrix elements are needed to predict the glueball branching
ratios in $J/\psi$ radiative decays which will help identify the glueball
states in experiments. Two types of improved local gluonic operators are
constructed for a selfconsistent check and the finite volume effects are
studied. We find that lattice spacing dependence of our results is very weak
and the continuum limits are reliably extrapolated, as a result of improvement
of the lattice gauge action and local operators. We also give updated glueball
masses with various quantum numbers.

The glueballtovacuum matrix elements of local gluonic operators in scalar,
tensor, and pseudoscalar channels are investigated numerically on several
anisotropic lattices with the spatial lattice spacing in the range 0.1fm 
0.2fm. These matrix elements are needed to predict the glueball branching
ratios in $J/\psi$ radiative decays which will help to identify the glueball
states in experiments. Two types of improved local gluonic operators are
constructed for a selfconsistent check, and the finite volume effects are also
studied. The lattice spacing dependence of our results is very small and the
continuum limits are reliably extrapolated.

We describe a lattice calculation of the matrix elements relevant for
glueball production in $J / \psi$ radiative decays. The techniques for such a
calculation on anisotropic lattices with an improved action are outlined. We
present preliminary results showing the efficacy of the computational method.