
Neutrons produced by cosmic ray muons are an important background for
underground experiments studying neutrino oscillations, neutrinoless double
beta decay, dark matter, and other rareevent 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
meterswaterequivalent. 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 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.

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 fueldependent 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 energydependent, 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 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.

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 Bugey3 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 fourneutrino flavor framework.
Stringent limits on $\sin^2 2\theta_{\mu e}$ are set over 6 orders of magnitude
in the sterile masssquared splitting $\Delta m^2_{41}$. The sterileneutrino
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.9GW$_{\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$15001950 m) relative to
detectors near the reactors ($\sim$350600 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. 2011Jul. 2012) with a
subsequent 1013 days using the complete configuration of eight detectors (Oct.
2012Jul. 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 threeflavor
oscillation probability yielded the flavormixing angle
$\sin^22\theta_{13}=0.0841 \pm 0.0027(\mathrm{stat.}) \pm
0.0019(\mathrm{syst.})$ and the effective neutrino masssquared 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
threeflavor 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 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.

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 dataset 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 217day data set with six detectors. The
dominant statistical uncertainty was reduced by 49%. Intensive studies of the
cosmogenic muoninduced $^9$Li and fast neutron backgrounds and the
neutroncapture 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
threeneutrinooscillation framework. The combination of this result with the
gadoliniumcapture result is also reported.

We analysis the radial equations for massive Dirac fields in Schwarzschild
black hole spacetimes. Different approximation formulae under the WKB scheme
are developed for the transmission probability ${\cal T}$ of the radial
wavefunction with $E^{2}\gg V_{m}$, $E^{2}\approx V_{m}$, and $E^{2}\ll V_{m}$,
where $E$ is the energy of the field and $V_{m}$ is the maximum value of the
effective potential. Explicit results of ${\cal T}$ in these approximations are
given for various values of $E$, the mass $m$, and the angular momentum
parameter $\kappa$ of the fields. We also discuss the dependence of ${\cal T}$
on these parameters.

The effect of quenched disorder on the lowenergy and lowtemperature
properties of various two and threedimensional Heisenberg models is studied
by a numerical strong disorder renormalization group method. For strong enough
disorder we have identified two relevant fixed points, in which the gap
exponent, omega, describing the lowenergy tail of the gap distribution,
P(Delta) ~ Delta^omega is independent of disorder, the strength of couplings
and the value of the spin. The dynamical behavior of nonfrustrated random
antiferromagnetic models is controlled by a singletlike fixed point, whereas
for frustrated models the fixed point corresponds to a large spin formation and
the gap exponent is given by omega ~ 0. Another type of universality classes is
observed at quantum critical points and in dimerized phases but no infinite
randomness behavior is found, in contrast to onedimensional models.

The effect of quenched disorder on the lowenergy properties of various
antiferromagnetic spin ladder models is studied by a numerical strong disorder
renormalization group method and by density matrix renormalization. For strong
enough disorder the originally gapped phases with finite topological or dimer
order become gapless. In these quantum Griffiths phases the scaling of the
energy, as well as the singularities in the dynamical quantities are
characterized by a finite dynamical exponent, z, which varies with the strength
of disorder. At the phase boundaries, separating topologically distinct
Griffiths phases the singular behavior of the disordered ladders is generally
controlled by an infinite randomness fixed point.

The quantum critical behavior and the GriffithsMcCoy singularities of random
quantum Ising ferromagnets are studied by applying a numerical implementation
of the MaDasguptaHu renormalization group scheme. We check the procedure for
the analytically tractable onedimensional case and apply our code to the
quasionedimensional double chain. For the latter we obtain identical critical
exponents as for the simple chain implying the same universality class. Then we
apply the method to the twodimensional case for which we get estimates for the
exponents that are compatible with a recent study in the same spirit.