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

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

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 kmbaselines. 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 datataking.

We report a new measurement of electron antineutrino disappearance using the
fullyconstructed 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}$tondays, 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 threeneutrino
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 shortbaseline 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 46 MeV. A reactor antineutrino spectrum of IBD reactions is extracted
from the measured positron energy spectrum for modelindependent 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
threeflavor 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
gadoliniumcapture result and an improvement on the precision of $\theta_{13}$
measurement. With a 217day 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 threeflavor oscillation model. When combined with the
gadoliniumcapture result from Daya Bay, we obtain
$\sin^22\theta_{13}=0.089\pm0.008$ as the final result for the
sixantineutrinodetector configuration of the Daya Bay experiment.

We report an improved measurement of the neutrino mixing angle $\theta_{13}$
from the Daya Bay Reactor Neutrino Experiment. We exclude a zero value for
$\sin^22\theta_{13}$ with a significance of 7.7 standard deviations. Electron
antineutrinos from six reactors of 2.9 GW$_{\rm th}$ were detected in six
antineutrino detectors deployed in two near (fluxweighted baselines of 470 m
and 576 m) and one far (1648 m) underground experimental halls. Using 139 days
of data, 28909 (205308) electron antineutrino candidates were detected at the
far hall (near halls). The ratio of the observed to the expected number of
antineutrinos assuming no oscillations at the far hall is $0.944\pm 0.007({\rm
stat.}) \pm 0.003({\rm syst.})$. An analysis of the relative rates in six
detectors finds $\sin^22\theta_{13}=0.089\pm 0.010({\rm stat.})\pm0.005({\rm
syst.})$ in a threeneutrino framework.

The Daya Bay Reactor Neutrino Experiment has measured a nonzero value for
the neutrino mixing angle $\theta_{13}$ with a significance of 5.2 standard
deviations. Antineutrinos from six 2.9 GW$_{\rm th}$ reactors were detected in
six antineutrino detectors deployed in two near (fluxweighted baseline 470 m
and 576 m) and one far (1648 m) underground experimental halls. With a 43,000
tonGW_{\rm th}day livetime exposure in 55 days, 10416 (80376) electron
antineutrino candidates were detected at the far hall (near halls). The ratio
of the observed to expected number of antineutrinos at the far hall is
$R=0.940\pm 0.011({\rm stat}) \pm 0.004({\rm syst})$. A rateonly analysis
finds $\sin^22\theta_{13}=0.092\pm 0.016({\rm stat})\pm0.005({\rm syst})$ in a
threeneutrino framework.

The Daya Bay Reactor Neutrino Experiment is designed to determine precisely
the neutrino mixing angle $\theta_{13}$ with a sensitivity better than 0.01 in
the parameter sin$^22\theta_{13}$ at the 90% confidence level. To achieve this
goal, the collaboration will build eight functionally identical antineutrino
detectors. The first two detectors have been constructed, installed and
commissioned in Experimental Hall 1, with steady datataking beginning
September 23, 2011. A comparison of the data collected over the subsequent
three months indicates that the detectors are functionally identical, and that
detectorrelated systematic uncertainties exceed requirements.

$Om$ diagnostic can differentiate between different models of dark energy
without the accurate current value of matter density. We apply this geometric
diagnostic to dilaton dark energy(DDE) model and differentiate DDE model from
LCDM. We also investigate the influence of coupled parameter $\alpha$ on the
evolutive behavior of $Om$ with respect to redshift $z$. According to the
numerical result of $Om$, we get the current value of equation of state
$\omega_{\sigma0}$=0.952 which fits the WMAP5+BAO+SN very well.

Applying the parametrization of dark energy density, we can construct
directly independentmodel potentials. In BornInfeld type phantom dark energy
model, we consider four special parametrization equation of state parameter.
The evolutive behavior of dark energy density with respect to redshift $z$,
potentials with respect to $\phi$ and $z$ are shown mathematically. Moreover,
we investigate the effect of parameter $\eta$ upon the evolution of the
constructed potential with respect to $z$. These results show that the
evolutive behavior of constructed BornInfeld type dark energy model is quite
different from those of the other models.

Statefinder diagnostic is a useful method which can differ one dark energy
model from the others. The Statefinder pair $\{r, s\}$ is algebraically related
to the equation of state of dark energy and its first time derivative. We apply
in this paper this method to the dilaton dark energy model based on WeylScaled
induced gravitational theory. We investigate the effect of the coupling between
matter and dilaton when the potential of dilaton field is taken as the Mexican
hat form. We find that the evolving trajectory of our model in the $rs$
diagram is quite different from those of other dark energy models.

Using a new methodstatefinder diagnostic which can differ one dark energy
model from the others, we investigate in this letter the dynamics of
BornInfeld(BI) type dark energy model. The evolutive trajectory of BI type
dark energy with Mexican hat potential model with respect to $efolding$ time
$N$ is shown in the $r(s)$ diagram. When the parameter of noncanonical kinetic
energy term $\eta\to0$ or kinetic energy $\dot{\phi}^2\to0$, BI type dark
energy(Kessence) model reduces to Quintessence model or $\Lambda$CDM model
corresponding to the statefinder pair $\{r, s\}$=$\{1, 0\}$ respectively. As a
result, the the evolutive trajectory of our model in the $r(s)$ diagram in
Mexican hat potential is quite different from those of other dark energy
models.

In this paper, we regard dilaton in Weylscaled induced gravitational theory
as coupled Quintessence, which is called DCQ model by us. Parametrization of
the dark energy model is a good method by which we can construct the scalar
potential directly from the effective equation of state function
$\omega_\sigma(z)$ describing the properties of the dark energy. Applying this
method to the DCQ model, we consider four parametrizations of $\omega(z)$ and
investigate the features of the constructed DCQ potentials, which possess two
different evolutive behaviors called "O" mode and "E" mode. Lastly, we comprise
the results of the constructed DCQ model with those of quintessence model
numerically.

We generally investigate the scalar field model with the lagrangian
$L=F(X)V(\phi)$, which we call it {\it General NonCanonical Scalar Field
Model}. We find that it is a special square potential(with a negative minimum)
that drives the linear field solution($\phi=\phi_0t$) while in Kessence
model(with the lagrangian $L=V(\phi)F(X)$) the potential should be taken as an
inverse square form. Hence their cosmological evolution are totally different.
We further find that this linear field solutions are highly degenerate, and
their cosmological evolutions are actually equivalent to the divergent model
where its sound speed diverges. We also study the stability of the linear field
solution. With a simple form of $F(X)=1\sqrt{12X}$ we indicate that our model
may be considered as a unified model of dark matter and dark energy. Finally we
study the case when the baryotropic index $\gamma$ is constant. It shows that,
unlike the Kessence, the detailed form of F(X) depends on the potential
$V(\phi)$. We analyze the stability of this constant $\gamma_0$ solution and
find that they are stable for $\gamma_0\leq1$. Finally we simply consider the
constant c_s^2 case and get an exact solution for F(X)

In this paper, we regard dilaton in Weylscaled induced gravitational theory
as a coupled quintessence. Based on this consideration, we investigate the
dilaton coupled quintessence(DCQ) model in $\omega\omega'$ plane, which is
defined by the equation of state parameter for the dark energy and its
derivative with respect to $N$(the logarithm of the scale factor $a$). We find
the scalar field equation of motion in $\omega\omega'$ plane, and show
mathematically the property of attractor solutions which correspond to
$\omega_\sigma\sim1$, $\Omega_\sigma=1$. Finally, we find that our model is a
tracking one which belongs to "freezing" type model classified in
$\omega\omega'$ plane.

Based on dilatonic dark energy model, we consider two cases: dilaton field
with positive kinetic energy(coupled quintessence) and with negative kinetic
energy(phantom). In the two cases, we investigate the existence of attractor
solutions which correspond to an equation of state parameter $\omega=1$ and a
cosmic density parameter $\Omega_\sigma=1$. We find that the coupled term
between matter and dilaton can't affect the existence of attractor solutions.
In the Mexican hat potential, the attractor behaviors, the evolution of state
parameter $\omega$ and cosmic density parameter $\Omega$, are shown
mathematically. Finally, we show the effect of coupling term on the evolution
of $X(\frac{\sigma}{\sigma_0})$ and $Y(\frac{\dot{\sigma}}{\sigma^2_0})$ with
respect to $N(lna)$ numerically.

Recent many physicists suggest that the dark energy in the universe might
result from the BornInfeld(BI) type scalar field of string theory. The
universe of BI type scalar field with potential can undergo a phase of
accelerating expansion. The corresponding equation of state parameter lies in
the range of $\displaystyle 1<\omega<{1/3}$. The equation of state parameter
of BI type scalar field without potential lies in the range of
$0\leq\omega\leq1$. We find that weak energy condition and strong energy
condition are violated for phantom BI type scalar field. The equation of state
parameter lies in the range of $\omega<1$.

On this paper we consider the classical wormhole solution of the BornInfeld
scalar field. The corresponding classical wormhole solution can be obtained
analytically for both very small and large $\dot{\phi}$. At the extreme limits
of small $\dot{\phi}$ the wormhole solution has the same format as one obtained
by Giddings and Strominger[10]. At the extreme limits of large $\dot{\phi}$ the
wormhole solution is a new one. The wormhole wavefunctions can also be obtained
for both very small and large $\dot{\phi}$. These wormhole wavefunctions are
regarded as solutions of quantummechanical WheelerDewitt equation with
certain boundary conditions.