
Reactor neutrinos have been an important tool for both discovery and
precision measurement in the history of neutrino studies. Since the first
generation of reactor neutrino experiments in the 1950s, the detector
technology has been greatly advanced. New ideas, new knowledge, and modern
software also enhanced the power of the experiments. The current reactor
neutrino experiments, Daya Bay, Double Chooz, and RENO have led neutrino
physics into the precision era. In this article, we will review these
developments and accumulations, address the key issues in designing a
stateofart reactor neutrino experiment, and explain how the challenging
requirements of determining the neutrino mass hierarchy with the next
generation experiment JUNO could be realized in the near future.

Linear alkylbenzene (LAB) based liquid scintillator is adopted as the central
detector for the Jiangmen Underground Neutrino Observatory (JUNO) liquid
scintillator detectors. A quenching factor measurement instrument is designed
based on the Compton scattering process. Two different quenchers for the liquid
scintillator have been investigated and the result shows that the scintillation
light of the JUNO liquid scintillator can be quenched to a level. The emission
spectrum with the absence of the quencher is also showing a desired behavior.

For any $p\in(0,\,1]$, let $H^{\Phi_p}(\mathbb{R}^n)$ be the MusielakOrlicz
Hardy space associated with the MusielakOrlicz growth function $\Phi_p$,
defined by setting, for any $x\in\mathbb{R}^n$ and $t\in[0,\,\infty)$,
$$ \Phi_{p}(x,\,t):= \begin{cases} \frac{t}{\log(e+t)+[t(1+x)^n]^{1p}} &
\qquad \text{when } n(1/p1)\notin \mathbb{N} \cup \{0\}; \\
\frac{t}{\log(e+t)+[t(1+x)^n]^{1p}[\log(e+x)]^p} & \qquad \text{when }
n(1/p1)\in \mathbb{N}\cup\{0\},\\ \end{cases} $$ which is the sharp target
space of the bilinear decomposition of the product of the Hardy space
$H^p(\mathbb{R}^n)$ and its dual. Moreover, $H^{\Phi_1}(\mathbb{R}^n)$ is the
prototype appearing in the realvariable theory of general MusielakOrlicz
Hardy spaces. In this article, the authors find a new structure of the space
$H^{\Phi_p}(\mathbb{R}^n)$ by showing that, for any $p\in(0,\,1]$,
$H^{\Phi_p}(\mathbb{R}^n)=H^{\phi_0}(\mathbb{R}^n) +H_{W_p}^p(\mathbb{R}^n)$
and, for any $p\in(0,\,1)$, $H^{\Phi_p}(\mathbb{R}^n)=H^{1}(\mathbb{R}^n)
+H_{W_p}^p(\mathbb{R}^n)$, where $H^1(\mathbb{R}^n)$ denotes the classical real
Hardy space, $H^{\phi_0}(\mathbb{R}^n)$ the OrliczHardy space associated with
the Orlicz function $\phi_0(t):=t/\log(e+t)$ for any $t\in [0,\infty)$ and
$H_{W_p}^p(\mathbb{R}^n)$ the weighted Hardy space associated with certain
weight function $W_p(x)$ that is comparable to $\Phi_p(x,1)$ for any
$x\in\mathbb{R}^n$. As an application, the authors further establish an
interpolation theorem of quasilinear operators based on this new structure.

In the past few decades, numerous searches have been made for the
neutrinoless doublebeta decay (0$\nu\beta\beta$) process, aiming to establish
whether neutrinos are their own antiparticles (Majorana neutrinos), but no
0$\nu\beta\beta$ decay signal has yet been observed. A number of new
experiments are proposed but they ultimately suffer from a common problem: the
sensitivity may not increase indefinitely with the target mass. We have
performed a detailed analysis of the physics potential by using the Jiangmen
Underground Neutrino Observatory (JUNO) to improve the sensitivity to
0$\nu\beta\beta$ up to a few meV, a major step forward with respect to the
experiments currently being planned. JUNO is a 20 kton lowbackground liquid
scintillator (LS) detector with 3\%/$\sqrt{E \text{(MeV)}}$ energy resolution,
now under construction. It is feasible to build a balloon filled with enriched
xenon gas (with $^{136}$Xe up to 80\%) dissolved in LS, inserted into the
central region of the JUNO LS. The energy resolution is $\sim$1.9\% at the
$Q$value of $^{136}$Xe 0$\nu\beta\beta$ decay. Ultralow background is the key
for 0$\nu\beta\beta$ decay searches. Detailed studies of background rates from
intrinsic 2$\nu\beta\beta$ and $^{8}$B solar neutrinos, natural radioactivity,
and cosmogenic radionuclides (including light isotopes and $^{137}$Xe) were
performed and several muon veto schemes were developed. We find that JUNO has
the potential to reach a sensitivity (at 90\% C. L.) to
$T^{0\nu\beta\beta}_{1/2}$ of $1.8\times10^{28}$ yr ($5.6\times10^{27}$ yr)
with $\sim$50 tons (5 tons) of fiducial $^{136}$Xe and 5 years exposure, while
in the 50ton case the corresponding sensitivity to the effective neutrino
mass, $m_{\beta\beta}$, could reach (512) meV, covering completely the
allowed region of inverted neutrino mass ordering.

The discovery of the neutrino mixing angle $\theta_{13}$ opens new
opportunities for the discovery of the leptonic CP violation for high intensity
neutrino beams. MOMENT a future neutrino facility with a highpower proton beam
of 15 MW from a continuouswave linac is focused on that discovery. The high
power of the proton beam causes extreme radiation conditions for the facility
and especially for the target station where the pion capture system of five
superconducting solenoids is located. In this paper initial studies are
performed for the effects of the radiation on the solenoid structure and the
area surrounding it. A concept cooling system is also proposed.

One of the most widelyimplemented service standards provided by the Open
Geospatial Consortium (OGC) to the user community is the Web Map Service (WMS).
WMS is widely employed globally, but there is limited knowledge of the global
distribution, adoption status or the service quality of these online WMS
resources. To fill this void, we investigated global WMSs resources and
performed distributed performance monitoring of these services. This paper
explicates a distributed monitoring framework that was used to monitor 46,296
WMSs continuously for over one year and a crawling method to discover these
WMSs. We analyzed server locations, provider types, themes, the spatiotemporal
coverage of map layers and the service versions for 41,703 valid WMSs.
Furthermore, we appraised the stability and performance of basic operations for
1210 selected WMSs (i.e., GetCapabilities and GetMap). We discuss the major
reasons for request errors and performance issues, as well as the relationship
between service response times and the spatiotemporal distribution of client
monitoring sites. This paper will help service providers, end users and
developers of standards to grasp the status of global WMS resources, as well as
to understand the adoption status of OGC standards. The conclusions drawn in
this paper can benefit geospatial resource discovery, service performance
evaluation and guide service performance improvements.

The Daya Bay Reactor Neutrino Experiment discovered an unexpectedly large
neutrino oscillation related to the mixing angle $\theta_{13}$ in 2012. This
finding paved the way to the next generation of neutrino oscillation
experiments. In this article, we review the history, featured design, and
scientific results of Daya Bay. Prospects of the experiment are also described.

The aim of this article is to give a complete solution to the problem of the
bilinear decompositions of the products of some Hardy spaces
$H^p(\mathbb{R}^n)$ and their duals in the case when $p<1$ and near to $1$, via
wavelets, paraproducts and the theory of bilinear Calder\'onZygmund operators.
Precisely, the authors establish the bilinear decompositions of the product
spaces $H^p(\mathbb{R}^n)\times\dot\Lambda_{\alpha} (\mathbb{R}^n)$ and
$H^p(\mathbb{R}^n)\times\Lambda_{\alpha}(\mathbb{R}^n)$, where, for all
$p\in(\frac{n}{n+1},\,1)$ and $\alpha:=n(\frac{1}{p}1)$, $H^p(\mathbb{R}^n)$
denotes the classical real Hardy space, and $\dot\Lambda_{\alpha}$ and
$\Lambda_{\alpha}$ denote the homogeneous, respectively, the inhomogeneous
Lipschitz spaces. Sharpness of these two bilinear decompositions are also
proved. As an application, the authors establish some divcurl lemmas at the
endpoint case.

The flux of geoneutrinos at any point on the Earth is a function of the
abundance and distribution of radioactive elements within our planet. This flux
has been successfully detected by the 1kt KamLAND and 0.3kt Borexino
detectors with these measurements being limited by their low statistics. The
planned 20kt JUNO detector will provide an exciting opportunity to obtain a
high statistics measurement, which will provide data to address several
questions of geological importance. This paper presents the JUNO detector
design concept, the expected geoneutrino signal and corresponding backgrounds.
The precision level of geoneutrino measurements at JUNO is obtained with the
standard leastsquares method. The potential of the Th/U ratio and mantle
measurements is also discussed.

The Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton multipurpose
underground liquid scintillator detector, was proposed with the determination
of the neutrino mass hierarchy as a primary physics goal. It is also capable of
observing neutrinos from terrestrial and extraterrestrial sources, including
supernova burst neutrinos, diffuse supernova neutrino background, geoneutrinos,
atmospheric neutrinos, solar neutrinos, as well as exotic searches such as
nucleon decays, dark matter, sterile neutrinos, etc. We present the physics
motivations and the anticipated performance of the JUNO detector for various
proposed measurements. By detecting reactor antineutrinos from two power plants
at 53km distance, JUNO will determine the neutrino mass hierarchy at a 34
sigma significance with six years of running. The measurement of antineutrino
spectrum will also lead to the precise determination of three out of the six
oscillation parameters to an accuracy of better than 1\%. Neutrino burst from a
typical corecollapse supernova at 10 kpc would lead to ~5000
inversebetadecay events and ~2000 allflavor neutrinoproton elastic
scattering events in JUNO. Detection of DSNB would provide valuable information
on the cosmic starformation rate and the average corecollapsed neutrino
energy spectrum. Geoneutrinos can be detected in JUNO with a rate of ~400
events per year, significantly improving the statistics of existing geoneutrino
samples. The JUNO detector is sensitive to several exotic searches, e.g. proton
decay via the $p\to K^++\bar\nu$ decay channel. The JUNO detector will provide
a unique facility to address many outstanding crucial questions in particle and
astrophysics. It holds the great potential for further advancing our quest to
understanding the fundamental properties of neutrinos, one of the building
blocks of our Universe.

Based on the standard Gaisser's formula, a modified parametrization for the
sealevel cosmicray muon flux is introduced. The modification is verified
against experimental results. The average vertical cosmicray muon intensity as
a function of depth of standard rock is simulated using the modified formula as
input to the MUSIC code. The calculated muon intensities is consistent with the
experimental measurements.

Rayleigh scattering poses an intrinsic limit for the transparency of organic
liquid scintillators. This work focuses on the Rayleigh scattering length of
linear alkylbenzene (LAB), which will be used as the solvent of the liquid
scintillator in the central detector of the Jiangmen Underground Neutrino
Observatory. We investigate the anisotropy of the Rayleigh scattering in LAB,
showing that the resulting Rayleigh scattering length will be significantly
shorter than reported before. Given the same overall light attenuation, this
will result in a more efficient transmission of photons through the
scintillator, increasing the amount of light collected by the photosensors and
thereby the energy resolution of the detector.

A highstatistics measurement of the neutrinos from a galactic corecollapse
supernova is extremely important for understanding the explosion mechanism, and
studying the intrinsic properties of neutrinos themselves. In this paper, we
explore the possibility to constrain the absolute scale of neutrino masses
$m^{}_\nu$ via the detection of galactic supernova neutrinos at the Jiangmen
Underground Neutrino Observatory (JUNO) with a 20 kiloton liquidscintillator
detector. In assumption of a nearlydegenerate neutrino mass spectrum and a
normal mass ordering, the upper bound on the absolute neutrino mass is found to
be $m^{}_\nu < (0.83 \pm 0.24)~{\rm eV}$ at the 95% confidence level for a
typical galactic supernova at a distance of 10 kpc, where the mean value and
standard deviation are shown to account for statistical fluctuations. For
comparison, we find that the bound in the SuperKamiokande experiment is
$m^{}_\nu < (0.94 \pm 0.28)~{\rm eV}$ at the same confidence level. However,
the upper bound will be relaxed when the model parameters characterizing the
time structure of supernova neutrino fluxes are not exactly known, and when the
neutrino mass ordering is inverted.

Let $L$ be a homogeneous divergence form higher order elliptic operator with
complex bounded measurable coefficients and $(p_(L),\, p_+(L))$ be the maximal
interval of exponents $q\in[1,\,\infty]$ such that the semigroup
$\{e^{tL}\}_{t>0}$ is bounded on $L^q(\mathbb{R}^n)$. In this article, the
authors establish the nontangential maximal function characterizations of the
associated Hardy spaces $H_L^p(\mathbb{R}^n)$ for all $p\in(0,\,p_+(L))$,
which, when $p=1$, answers a question asked by Deng et al. in [J. Funct. Anal.
263 (2012), 604674]. Moreover, the authors characterize $H_L^p(\mathbb{R}^n)$
via various versions of square functions and Lusinarea functions associated to
the operator $L$.

We has set up a light scattering spectrometer to study the depolarization of
light scattering in linear alkylbenzene. From the scattering spectra it can be
unambiguously shown that the depolarized part of light scattering belongs to
Rayleigh scattering. The additional depolarized Rayleigh scattering can make
the effective transparency of linear alkylbenzene much better than it was
expected. Therefore sufficient scintillation photons can transmit through the
large liquid scintillator detector of JUNO. Our study is crucial to achieving
the unprecedented energy resolution 3\%/$\sqrt{E\mathrm{(MeV)}}$ for JUNO
experiment to determine the neutrino mass hierarchy. The spectroscopic method
can also be used to judge the attribution of the depolarization of other
organic solvents used in neutrino experiments.

Linear alkylbenzene (LAB) is adopted to be the organic solvent for the
Jiangmen Underground Neutrino Observatory (JUNO) liquid scintillator detectors
due to the ultratransparency. However the current Rayleigh scattering length
calculation disagrees with the measurement. The present paper for the first
time reports the Rayleigh scattering of LAB being anisotropic and the
depolarization ratio being 0.31+0.01(stat.)+0.01(sys.). We proposed an
indirectly method for Rayleigh scattering measurement with
EinsteinSmoluchowskiCabannes formula, and the Rayleigh scattering length of
LAB is determined to be 28.2+1.0 m at 430 nm.

We report the measurements of the densities of linear alkylbenzene at three
temperatures over 4 to 23 Celsius degree with pressures up to 10 MPa. The
measurements have been analysed to yield the isobaric thermal expansion
coefficients and, so far for the first time, isothermal compressibilities of
linear alkylbenzene. Relevance of results for current generation (i.e. Daya
Bay) and next generation (i.e. JUNO) large liquid scintillator neutrino
detectors are discussed.

Let $\varphi$ be a MusielakOrlicz function satisfying that, for any
$(x,\,t)\in\mathbb{R}^n\times(0,\,\infty)$, $\varphi(\cdot,\,t)$ belongs to the
Muckenhoupt weight class $A_\infty (\mathbb{R}^n)$ with the critical weight
exponent $q(\varphi)\in[1,\,\infty)$ and $\varphi(x,\,\cdot)$ is an Orlicz
function with $0<i(\varphi)\le I(\varphi)\le 1$ which are, respectively, its
critical lower type and upper type. In this article, the authors establish the
Riesz transform characterizations of the MusielakOrliczHardy spaces
$H_\varphi (\mathbb{R}^n)$ which are generalizations of weighted Hardy spaces
and OrliczHardy spaces. Precisely, the authors characterize $H_\varphi
(\mathbb{R}^n)$ via all the first order Riesz transforms when
$\frac{i(\varphi)}{q(\varphi)}>\frac{n1}{n}$, and via all the Riesz transforms
with the order not more than $m\in\mathbb{N}$ when
$\frac{i(\varphi)}{q(\varphi)}>\frac{n1}{n+m1}$. Moreover, the authors also
establish the Riesz transform characterizations of $H_\varphi(\mathbb{R}^n)$,
respectively, by means of the higher order Riesz transforms defined via the
homogenous harmonic polynomials or the odd order Riesz transforms. Even if when
$\varphi(x,t):=tw(x)$ for all $x\in{\mathbb R}^n$ and $t\in [0,\infty)$, these
results also widen the range of weights in the known Riesz characterization of
the classical weighted Hardy space $H^1_w({\mathbb R}^n)$ obtained by R. L.
Wheeden from $w\in A_1({\mathbb R}^n)$ into $w\in A_\infty({\mathbb R}^n)$ with
the sharp range $q(w)\in [1,\frac n{n1})$, where $q(w)$ denotes the critical
index of the weight $w$.

Let $L$ be a onetoone operator of type $\omega$ having a bounded $H_\infty$
functional calculus and satisfying the $k$DaviesGaffney estimates with
$k\in{\mathbb N}$. In this paper, the authors introduce the weak Hardy space
$WH_L^p(\mathbb{R}^n)$ associated to $L$ for $p\in (0,\,1]$ via the
nontangential square function $S_L$ and establish a weak molecular
characterization of $WH_L^p(\mathbb{R}^n)$. Typical examples of such operators
include the $2k$order divergence form homogeneous elliptic operator
$L_1:=(1)^k\sum_{\alpha=k=\beta}\partial^\beta(a_{\alpha,\beta}\partial^\alpha)$,
where $\{a_{\alpha,\beta}\}_{\alpha=k=\beta}$ are complex bounded
measurable functions, and the $2k$order Schr\"odinger type operator $L_2:=
(\Delta)^k+V^k$, where $\Delta$ is the Laplacian operator and $0\le V\in
L^k_{\mathop\mathrm{loc}}(\mathbb{R}^n)$. As applications, for $i\in\{1,\,2\}$
and $p\in(\frac{n}{n+k},\,1]$, the authors prove that the associated Riesz
transform $\nabla^k (L_i^{1/2})$ is bounded from $WH^p_{L_i}(\mathbb{R}^n)$ to
the classical weak Hardy space $WH^p(\mathbb{R}^n)$ and, for all $0<p<r\le1$
and $\alpha=n(\frac{1}{p}\frac{1}{r})$, the fractional power
$L_i^{\frac{\alpha}{2k}}$ is bounded from $WH_{L_i}^p(\mathbb{R}^n)$ to
$WH_{L_i}^r(\mathbb{R}^n)$. Furthermore, the authors find the dual space of
$WH_L^p(\mathbb{R}^n)$ for $p\in(0,\,1]$, which can be defined via mean
oscillations based on some subtle coverings of bounded open sets and, even when
$L:=\Delta$, are also previously unknown. In particular, if $L$ is a
nonnegative selfadjoint operator in $L^2({\mathbb R}^n)$ satisfying the
DaviesGaffney estimates, the authors further establish the weak atomic
characterization of $WH_L^p(\mathbb{R}^n)$.

We design compact waveguide couplers via impedancetunable transformation
optics. By tuning impedance coefficients in the original space, twodimensional
metallic and dielectric waveguide couplers are designed with a high efficiency.
Through tuning refractive index simultaneously, we find that the transformation
medium inside a designed metallic waveguide coupler can be a reducedparameter
material for coupling waves between waveguides with arbitrary different cross
sections and embedded media. In the design of dielectric waveguide couplers, we
apply two different schemes: one is that both core region and its cladding
region are contained in a transformed space, the other is only core region
contained in the transformed space. The former has a very high efficiency near
100%; the latter is less efficient with a very small decline which can be a
simplified candidate in the design of nearperfect dielectric waveguide
couplers. The transformation medium for dielectric waveguide couplers can also
be reducedparameter material by selecting appropriate refractive index
coefficients. Twodimensional numerical simulations confirm our design with
good performances.

A system employing a desktop FADC has been developed to investigate the
features of 8 inches Hamamatsu PMT. The system stands out for its highspeed
and informative results as a consequence of adopting fast waveform sampling
technology. Recording full waveforms allows us to perform digital signal
processing, pulse shape analysis, and precision timing extraction. High
precision after pulse time and charge distribution characteristics are
presented in this manuscript. Other photomultipliers characteristics, such as
dark rate and transit time spread, can also be obtained by exploiting waveform
analysis using this system.

Nonlinearity of the liquid scintillator energy response is a key to measuring
the neutrino energy spectrum in reactor neutrino experiments such as Daya Bay
and JUNO. We measured in laboratory the nonlinearity of the Linear Alkyl
Benzene based liquid scintillator, which is used in Daya Bay and will be used
in JUNO, via Compton scattering process. By tagging the scattered gamma from
the liquid scintillator sample simultaneously at seven angles, the instability
of the system was largely cancelled. The accurately measured nonlinearity will
improve the precision of the $\theta_{13}$, $\Delta m^2$, and reactor neutrino
spectrum measurements at Daya Bay.

Neutrino beam with about 300 MeV in energy, highflux and medium baseline is
considered a rational choice for measuring CP violation before the more
powerful Neutrino Factory will be built. Following this concept, a unique
neutrino beam facility based on muondecayed neutrinos is proposed. The
facility adopts a continuouswave proton linac of 1.5 GeV and 10 mA as the
proton driver, which can deliver an extremely high beam power of 15 MW. Instead
of piondecayed neutrinos, unprecedentedly intense muondecayed neutrinos are
used for better background discrimination. The schematic design for the
facility is presented here, including the proton driver, the assembly of a
mercuryjet target and capture superconducting solenoids, a pion/muon beam
transport line, a long muon decay channel of about 600 m and the detector
concept. The physics prospects and the technical challenges are also discussed.

We report on the production and characterization of liquid scintillators for
the detection of electron antineutrinos by the Daya Bay Reactor Neutrino
Experiment. One hundred eightyfive tons of gadoliniumloaded (0.1% by mass)
liquid scintillator (GdLS) and two hundred tons of unloaded liquid
scintillator (LS) were successfully produced from a linearalkylbenzene (LAB)
solvent in six months. The scintillator properties, the production and
purification systems, and the quality assurance and control (QA/QC) procedures
are described.

We report a new strategy to remove the reflections resulted from the finite
embedded transformationoptical design by proposing an impedancetunable
coordinate transformation, on which the functions of impedance coefficients can
be derived in the original space without changing the refractive index. Based
on the new approach, twodimensional (2D) reflectionless beam
compressors/expanders, bends, shifters and splitters are designed using the
modified anisotropic medium. It is found that the reflections can be removed in
magnetic response materials for TE polarization or dielectric response
materials for TM polarization. The numerical simulations confirm that various
reflectionless optical elements can be realized in the pure transformation
optics. The impedancetunable coordinate transformation can be generalized to
threedimensional (3D) cases and be applied to other transformationoptical
designs.