
Linear mixed models (LMMs) are used extensively to model dependecies of
observations in linear regression and are used extensively in many application
areas. Parameter estimation for LMMs can be computationally prohibitive on big
data. Stateoftheart learning algorithms require computational complexity
which depends at least linearly on the dimension $p$ of the covariates, and
often use heuristics that do not offer theoretical guarantees. We present
scalable algorithms for learning highdimensional LMMs with sublinear
computational complexity dependence on $p$. Key to our approach are novel dual
estimators which use only kernel functions of the data, and fast computational
techniques based on the subsampled randomized Hadamard transform. We provide
theoretical guarantees for our learning algorithms, demonstrating the
robustness of parameter estimation. Finally, we complement the theory with
experiments on large synthetic and real data.

It has been shown that neutrino masses can be determined under the particle
ansatz. In this paper, we give the general formulas of neutrino masses related
to the neutrino oscillation parameters which show that there is a mass
hierarchy transition. Using the best fit values measured by electron and muon
neutrino oscillation experiments, the total neutrino mass is about 0.25 eV.
According to the standard neutrino cosmology, the neutrino energy density is
about 55 eV/cm$^3$ and $\Omega_m h^2>\Omega_\nu h^2\sim 0.0055$.

We propose an ansatz, named the "particle ansatz", that in the rest frame of
a flavor neutrino all its mass eigenstate momentums should simultaneously be
zero. Under the particle ansatz, the mass measurable phase differences (mmPhDs)
can be derived as $\Delta\phi_{ji}\simeq [2m_{\nu_\alpha}(m_jm_i)/4
E_{\nu_\alpha}]L$. We show that using mmPhDs all three neutrino masses can be
determined. In the foreseeable future all neutrino masses are expected to be
determined by the spectral analysis of neutrino experiments.

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.

Hawkes processes are a class of simple point processes that are selfexciting
and have clustering effect, with wide applications in finance, social networks
and many other fields. This paper considers a selfexciting Hawkes process
where the baseline intensity is timedependent, the exciting function is a
general function and the jump sizes of the intensity process are independent
and identically distributed nonnegative random variables. This Hawkes model is
nonMarkovian in general. We obtain closedform formulas for the Laplace
transform, moments and the distribution of the Hawkes process. To illustrate
the applications of our results, we use the Hawkes process to model the
clustered arrival of trades in a dark pool and analyze various performance
metrics including timetofirstfill, timetocompletefill and the expected
fill rate of a resting dark order.

We introduce varbvs, a suite of functions written in R and MATLAB for
regression analysis of largescale data sets using Bayesian variable selection
methods. We have developed numerical optimization algorithms based on
variational approximation methods that make it feasible to apply Bayesian
variable selection to very large data sets. With a focus on examples from
genomewide association studies, we demonstrate that varbvs scales well to data
sets with hundreds of thousands of variables and thousands of samples, and has
features that facilitate rapid data analyses. Moreover, varbvs allows for
extensive model customization, which can be used to incorporate external
information into the analysis. We expect that the combination of an easytouse
interface and robust, scalable algorithms for posterior computation will
encourage broader use of Bayesian variable selection in areas of applied
statistics and computational biology. The most recent R and MATLAB source code
is available for download at Github (https://github.com/pcarbo/varbvs), and the
R package can be installed from CRAN
(https://cran.rproject.org/package=varbvs).

Here we present a multiscale method to calculate the saddle point associated
with the effective dynamics arising from a stochastic system which couples slow
deterministic drift and fast stochastic dynamics. This problem is motivated by
the transition states on free energy surfaces in chemical physics. Our method
is based on the gentlest ascent dynamics which couples the position variable
and the direction variable and has the local convergence to saddle points. The
dynamics of the direction vector is derived in terms of the covariance function
with respective to the equilibrium distribution of the fast stochastic process.
We apply the multiscale numerical methods to efficiently solve the obtained
multiscale gentlest ascent dynamics, {and discuss the acceleration techniques
based on the adaptive idea.} The examples of stochastic ordinary and partial
differential equations are presented.

We consider the asymptotic solutions of an interface problem corresponding to
an elliptic partial differential equation with Dirich let boundary condition
and transmission condition, subject to the small geometric perturbation and the
high contrast ratio of the conductivity. We consider two types of
perturbations: the first corresponds to a thin layer coating a fixed bounded
domain and the second is the per perturbation of the interface. As the
perturbation size tends to zero and the ratio of the conductivities in two
subdomains tends to zero, the twoparameter asymptotic expansions on the fixed
reference domain are derived to any order after the single parameter expansions
are solved be forehand. Our main tool is the asymptotic analysis based on the
Taylor expansions for the properly extended solutions on fixed domains. The
Neumann boundary condition and Robin boundary condition arise in twoparameter
expansions, depending on the relation of the geometric perturbation size and
the contrast ratio.

Nonlinear kernel regression models are often used in statistics and machine
learning because they are more accurate than linear models. Variable selection
for kernel regression models is a challenge partly because, unlike the linear
regression setting, there is no clear concept of an effect size for regression
coefficients. In this paper, we propose a novel framework that provides an
effect size analog of each explanatory variable for Bayesian kernel regression
models when the kernel is shiftinvariant  for example, the Gaussian kernel.
We use function analytic properties of shiftinvariant reproducing kernel
Hilbert spaces (RKHS) to define a linear vector space that: (i) captures
nonlinear structure, and (ii) can be projected onto the original explanatory
variables. The projection onto the original explanatory variables serves as an
analog of effect sizes. The specific function analytic property we use is that
shiftinvariant kernel functions can be approximated via random Fourier bases.
Based on the random Fourier expansion we propose a computationally efficient
class of Bayesian approximate kernel regression (BAKR) models for both
nonlinear regression and binary classification for which one can compute an
analog of effect sizes. We illustrate the utility of BAKR by examining two
important problems in statistical genetics: genomic selection (i.e. phenotypic
prediction) and association mapping (i.e. inference of significant variants or
loci). Stateoftheart methods for genomic selection and association mapping
are based on kernel regression and linear models, respectively. BAKR is the
first method that is competitive in both settings.

The light unflavoured meson $\eta/\eta^{\prime}$ decays are valuable for
testing nonperturbation quantum chromodynamics and exploring new physics
beyond the Standard Model. The proper $\eta/\eta'$ generators will help to
obtain detection efficiencies and suppress backgrounds in BESIII physics. In
this paper, we develop a series of event generators on the BESIII detector to
study the decays of $\eta/\eta^{\prime}\to\gamma l^{+}l^{}$,
$\eta/\eta^{\prime}\to\gamma \pi^{+}\pi^{}$, $\eta^{\prime}\to\omega
e^{+}e^{}$, $\eta\to\pi^{+}\pi^{}\pi^{0}$,
$\eta/\eta^{\prime}\to\pi^{0}\pi^{0}\pi^{0}$, $\eta^{\prime}\to\eta\pi\pi$ and
$\eta'\to\pi^+\pi^\pi^+\pi^/\pi^+\pi^\pi^0\pi^0$. Monte Carlo simulation
distributions show that all event generators are performing very well. The
consistence between our simulation distributions and the BESIII results examine
the usability of the $\eta/\eta^{\prime}$ generators which can be also used in
other $\eta/\eta^{\prime}$ factories and improved by further experiment
results.

Partial differential equations with random inputs have become popular models
to characterize physical systems with uncertainty coming from, e.g., imprecise
measurement and intrinsic randomness. In this paper, we perform asymptotic rare
event analysis for such elliptic PDEs with random inputs. In particular, we
consider the asymptotic regime that the noise level converges to zero
suggesting that the system uncertainty is low, but does exists. We develop
sharp approximations of the probability of a large class of rare events.

In this article, we derive the asymptotic expansion, up to an arbitrary order
in theory, for the solution of a twodimensional elliptic equation with
strongly anisotropic diffusion coefficients along different directions, subject
to the Neumann boundary condition and the Dirichlet boundary condition on
specific parts of the domain boundary, respectively. The illposedness arising
from the Neumann boundary condition in the strongly anisotropic diffusion limit
is handled by the decomposition of the solution into a mean part and a
fluctuation part. The boundary layer analysis due to the Dirichlet boundary
condition is conducted for each order in the expansion for the fluctuation
part. Our results suggest that the leading order is the combination of the mean
part and the composite approximation of the fluctuation part for the general
Dirichlet boundary condition.

Among numerical methods for partial differential equations arising from
steepest descent dynamics of energy functionals (e.g., AllenCahn and
CahnHilliard equations), the convex splitting method is wellknown to maintain
unconditional energy stability for a large time step size. In this work, we
show how to use the convex splitting idea to solve the problem of finding
transition states, i.e., index1 saddle points of the same energy functionals.
Based on the previous work of iterative minimization formulation (IMF) for
saddle points (SIAM J. Numer. Anal., vol. 53, p1786, 2015), we introduce the
convex splitting idea to minimize the auxiliary functional at each cycle of the
IMF. We present a general principle of constructing convex splitting forms for
these auxiliary functionals and show how to avoid solving nonlinear equations.
The new numerical scheme based on the convex splitting method allows for a much
larger time step size than a traditional semiimplicit scheme we tested and it
turns out that a large time step has a huge advantage in computational
efficiency. The new methods are tested numerically for the one dimensional
GinzburgLandau energy functional in the search of the AllenCahn or
CahnHilliard types of transition states.

Jinping Neutrino Experiment (Jinping) is proposed to significantly improve
measurements on solar neutrinos and geoneutrinos in China Jinping Laboratory 
a lab with a number of unparalleled features, thickest overburden, lowest
reactor neutrino background, etc., which identify it as the worldbest
lowenergy neutrino laboratory. The proposed experiment will have target mass
of 4 kilotons of liquid scintillator or waterbased liquid scintillator, with a
fiducial mass of 2 kilotons for neutrinoelectron scattering events and 3
kilotons for inversebeta interaction events. A number of initial sensitivities
studies have been carried out, including on the transition phase for the solar
neutrinos oscillation from the vacuum to the matter effect, the discovery of
solar neutrinos from the carbonnitrogenoxygen (CNO) cycle, the resolution of
the high and low metallicity hypotheses, and the unambiguous separation on U
and Th cascade decays from the dominant crustal antielectron neutrinos in
China.

We consider the noiseinduced transitions in the randomly perturbed discrete
logistic map from a linearly stable periodic orbit consisting of T periodic
points. The traditional large deviation theory and asymptotic analysis for
small noise limit as well as the derived quasipotential can not distinguish
the quantitative difference in noiseinduced stochastic instabilities of these
T periodic points. We generalize the transition path theory to the
discretetime continuousspace stochastic process to attack this problem. As a
first criterion of quantifying the relative instability among T periodic
points, we compare the distribution of the last passage locations in the
transitions from the whole periodic orbit to a prescribed set far away. This
distribution is related to the contributions to the transition rate from each
periodic points. The second criterion is based on the competency of the
transition paths associated with each periodic point. Both criteria utilise the
reactive probability current in the transition path theory. Our numerical
results for the logistic map reveal the transition mechanism of escaping from
the stable periodic orbit and identify which peri odic point is more prone to
lose stability so as to make successful transitions under random perturbations.

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.

We discuss importance sampling schemes for the estimation of finite time exit
probabilities of small noise diffusions that involve escape from an
equilibrium. A factor that complicates the analysis is that rest points are
included in the domain of interest. We build importance sampling schemes with
provably good performance both preasymptotically, that is, for fixed size of
the noise, and asymptotically, that is, as the size of the noise goes to zero,
and that do not degrade as the time horizon gets large. Simulation studies
demonstrate the theoretical results.

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.

Liquid scintillator (LS) will be adopted as the detector material in JUNO
(Jiangmen Underground Neutrino Observatory). The energy resolution requirement
of JUNO is 3%, which has never previously been reached. To achieve this energy
resolution, the light yield of liquid scintillator is an important factor. PPO
(the fluor) and bisMSB (the wavelength shifter) are the two main materials
dissolved in LAB. To study the influence of these two materials on the
transmission of scintillation photons in LS, 25 and 12 cmlong quartz vessels
were used in a light yield experiment. LS samples with different concentration
of PPO and bisMSB were tested. At these lengths, the light yield growth is not
obvious when the concentration of PPO is higher than 4 g/L. The influence from
bisMSB becomes insignificant when its concentration is higher than 8 mg/L.
This result could provide some useful suggestions for the JUNO LS.

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.

Manipulation of antiferromagnetic (AFM) spins by electrical means is on great
demand to develop the AFM spintronics with low power consumption. In spite of
the electrical modulation of insulated AFMs through coupling between their
intrinsic ferroelectricity and antiferromagnetism, direct electrical control of
AFM metals remains challenging due to the screening effect by the surface
charge, and the manipulation is confined to a limited depth of atomic
dimensions, which is insufficient to form a stable AFM exchange spring. In the
present letter we primarily report a reversible electrical control of exchange
spring in AFM metals, using an ionic liquid to exert a substantial
electricfield effect. The exchange spring could transfer the force to the
ferromagnet/antiferromagnet interface, enabling a deeper modulation depth in
AFM metals. The manipulation of AFM moments by gate voltage is demonstrated in
[Co/Pt]/IrMn model system and a single IrMn layer with the IrMn thickness up to
5 nm. Besides the fundamental significance of modulating the spin structures in
metallic AFM via allelectrical fashion, the present finding would advance the
development of lowpowerconsumption AFM spintronics.

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.

We study the anomalous Halllike effect (AHLE) and the effective anisotropic
magnetoresistance (EAMR) in antiferromagnetic {\gamma} IrMn3/Y3Fe5O12(YIG) and
Pt/YIG heterostructures. For {\gamma} IrMn3/YIG, the EAMR and the AHLE
resistivity change sign with temperature due to the competition between the
spin Hall magnetoresistance (SMR) and the magnetic proximity effect (MPE)
induced by the interfacial antiferromagnetic uncompensated magnetic moment. In
contrast, for Pt/YIG the AHLE resistivity changes sign with temperature whereas
no sign change is observed in the EAMR. This is because the MPE and the SMR
play a dominant role in the AHLE and the EAMR, respectively. As new types of
galvanomagnetic property, the AHLE and the EAMR have proved vital in
disentangling the MPE and the SMR in metal/insulatingferromagnet
heterostructures.

We consider noiseinduced transition paths in randomly perturbed dynami cal
systems on a smooth manifold. The classical FreidlinWentzell large devia tion
theory in Euclidean spaces is generalized and new forms of action functionals
are derived in the spaces of functions and the space of curves to accommodate
the intrinsic constraints associated with the manifold. Numerical meth ods are
proposed to compute the minimum action paths for the systems with constraints.
The examples of conformational transition paths for a single and double rod
molecules arising in polymer science are numerically investigated.