• ### Scalable Algorithms for Learning High-Dimensional Linear Mixed Models(1803.04431)

March 12, 2018 stat.ML
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. State-of-the-art 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 high-dimensional 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.
• ### Estimation on the neutrino masses by using neutrino oscillation parameters(1802.10042)

Feb. 27, 2018 hep-ph, astro-ph.CO
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$.
• ### Mass measurable phase differences in neutrino oscillations(1802.09270)

Feb. 26, 2018 hep-ph, astro-ph.CO
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_j-m_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.
• ### Quenching of fluorescence for linear alkylbenzene(1801.04432)

Jan. 13, 2018 hep-ex, physics.ins-det
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.
• ### Transform Analysis for Hawkes Processes with Applications in Dark Pool Trading(1710.01452)

Oct. 4, 2017 math.PR, q-fin.TR
Hawkes processes are a class of simple point processes that are self-exciting and have clustering effect, with wide applications in finance, social networks and many other fields. This paper considers a self-exciting Hawkes process where the baseline intensity is time-dependent, the exciting function is a general function and the jump sizes of the intensity process are independent and identically distributed non-negative random variables. This Hawkes model is non-Markovian in general. We obtain closed-form 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 time-to-first-fill, time-to-complete-fill and the expected fill rate of a resting dark order.
• ### varbvs: Fast Variable Selection for Large-scale Regression(1709.06597)

Sept. 19, 2017 stat.CO, q-bio.QM
We introduce varbvs, a suite of functions written in R and MATLAB for regression analysis of large-scale 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 genome-wide 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 easy-to-use 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.r-project.org/package=varbvs).
• ### Multiscale Gentlest Ascent Dynamics for Saddle Point in Effective Dynamics of Slow-Fast System(1602.00953)

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.
• ### Two-parameter asymptotic expansions for elliptic equations with small geometric perturbation and high contrast ratio(1708.04385)

Aug. 15, 2017 math.AP
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 two-parameter 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 two-parameter expansions, depending on the relation of the geometric perturbation size and the contrast ratio.
• ### Bayesian Approximate Kernel Regression with Variable Selection(1508.01217)

June 10, 2017 q-bio.QM, stat.ME, stat.AP, stat.ML
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 shift-invariant --- for example, the Gaussian kernel. We use function analytic properties of shift-invariant 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 shift-invariant 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). State-of-the-art 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.
• ### Event generators for $\eta/\eta^{\prime}$ decays on BESIII(1703.09920)

March 30, 2017 hep-ph, hep-ex
The light unflavoured meson $\eta/\eta^{\prime}$ decays are valuable for testing non-perturbation 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.
• ### Moderate Deviation for Random Elliptic PDEs with Small Noise(1703.05285)

March 16, 2017 math.PR
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.
• ### Asymptotic Expansion with Boundary Layer Analysis for Strongly Anisotropic Elliptic Equations(1701.03353)

Jan. 9, 2017 math.AP
In this article, we derive the asymptotic expansion, up to an arbitrary order in theory, for the solution of a two-dimensional 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 ill-posedness 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.
• ### Convex Splitting Method for the Calculation of Transition States of Energy Functional(1610.07153)

Oct. 23, 2016 math.NA
Among numerical methods for partial differential equations arising from steepest descent dynamics of energy functionals (e.g., Allen-Cahn and Cahn-Hilliard equations), the convex splitting method is well-known 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., index-1 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 semi-implicit 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 Ginzburg-Landau energy functional in the search of the Allen-Cahn or Cahn-Hilliard types of transition states.
• ### Letter of Intent: Jinping Neutrino Experiment(1602.01733)

June 21, 2016 hep-ex, physics.ins-det
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 world-best low-energy neutrino laboratory. The proposed experiment will have target mass of 4 kilotons of liquid scintillator or water-based liquid scintillator, with a fiducial mass of 2 kilotons for neutrino-electron scattering events and 3 kilotons for inverse-beta 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 carbon-nitrogen-oxygen (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 anti-electron neutrinos in China.
• ### Explore Stochastic Instabilities of Periodic Points by Transition Path Theory(1601.01424)

Jan. 7, 2016 nlin.CD
We consider the noise-induced 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 quasi-potential can not distinguish the quantitative difference in noise-induced stochastic instabilities of these T periodic points. We generalize the transition path theory to the discrete-time continuous-space 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 multi-purpose 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 extra-terrestrial 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 53-km distance, JUNO will determine the neutrino mass hierarchy at a 3-4 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 core-collapse supernova at 10 kpc would lead to ~5000 inverse-beta-decay events and ~2000 all-flavor neutrino-proton elastic scattering events in JUNO. Detection of DSNB would provide valuable information on the cosmic star-formation rate and the average core-collapsed neutrino energy spectrum. Geo-neutrinos 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.
• ### Escaping from an attractor: Importance sampling and rest points I(1303.0450)

Sept. 9, 2015 math.PR, math.OC
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 pre-asymptotically, 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 of linear alkylbenzene in large liquid scintillator detectors(1504.00987)

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.
• ### Preliminary study of light yield dependence on LAB liquid scintillator composition(1506.00237)

May 31, 2015 hep-ex, physics.ins-det
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 bis-MSB (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 cm-long quartz vessels were used in a light yield experiment. LS samples with different concentration of PPO and bis-MSB 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 bis-MSB becomes insignificant when its concentration is higher than 8 mg/L. This result could provide some useful suggestions for the JUNO LS.
• ### Spectroscopic study of light scattering in linear alkylbenzene for liquid scintillator neutrino detectors(1504.00986)

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.
• ### Rayleigh scattering and depolarization ratio in linear alkylbenzene(1504.01001)

April 7, 2015 hep-ex, physics.ins-det
Linear alkylbenzene (LAB) is adopted to be the organic solvent for the Jiangmen Underground Neutrino Observatory (JUNO) liquid scintillator detectors due to the ultra-transparency. 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 Einstein-Smoluchowski-Cabannes formula, and the Rayleigh scattering length of LAB is determined to be 28.2+-1.0 m at 430 nm.
• ### Electrical control of exchange spring in antiferromagnetic metals(1504.01193)

April 6, 2015 cond-mat.mtrl-sci
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 electric-field 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 all-electrical fashion, the present finding would advance the development of low-power-consumption AFM spintronics.
• ### Densities, isobaric thermal expansion coefficients and isothermal compressibilities of linear alkylbenzene(1408.0877)

April 4, 2015 hep-ex, physics.ins-det
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.
• ### Nature of magnetotransport in metal/insulating-ferromagnet heterostructures: Spin Hall magnetoresistance or magnetic proximity effect(1503.07388)

We study the anomalous Hall-like 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/insulating-ferromagnet heterostructures.
• ### Finding Transition Pathways on Manifolds(1408.3602)

Aug. 13, 2014 math-ph, math.MP
We consider noise-induced transition paths in randomly perturbed dynami- cal systems on a smooth manifold. The classical Freidlin-Wentzell 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.