• 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.
  • 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_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.
  • 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 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.
  • 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).
  • 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 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.
  • 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.
  • 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.
  • 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 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.
  • 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.
  • 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.
  • 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.
  • Fengpeng An, Guangpeng An, Qi An, Vito Antonelli, Eric Baussan, John Beacom, Leonid Bezrukov, Simon Blyth, Riccardo Brugnera, Margherita Buizza Avanzini, Jose Busto, Anatael Cabrera, Hao Cai, Xiao Cai, Antonio Cammi, Guofu Cao, Jun Cao, Yun Chang, Shaomin Chen, Shenjian Chen, Yixue Chen, Davide Chiesa, Massimiliano Clemenza, Barbara Clerbaux, Janet Conrad, Davide D'Angelo, Herve De Kerret, Zhi Deng, Ziyan Deng, Yayun Ding, Zelimir Djurcic, Damien Dornic, Marcos Dracos, Olivier Drapier, Stefano Dusini, Stephen Dye, Timo Enqvist, Donghua Fan, Jian Fang, Laurent Favart, Richard Ford, Marianne Goger-Neff, Haonan Gan, Alberto Garfagnini, Marco Giammarchi, Maxim Gonchar, Guanghua Gong, Hui Gong, Michel Gonin, Marco Grassi, Christian Grewing, Mengyun Guan, Vic Guarino, Gang Guo, Wanlei Guo, Xin-Heng Guo, Caren Hagner, Ran Han, Miao He, Yuekun Heng, Yee Hsiung, Jun Hu, Shouyang Hu, Tao Hu, Hanxiong Huang, Xingtao Huang, Lei Huo, Ara Ioannisian, Manfred Jeitler, Xiangdong Ji, Xiaoshan Jiang, Cecile Jollet, Li Kang, Michael Karagounis, Narine Kazarian, Zinovy Krumshteyn, Andre Kruth, Pasi Kuusiniemi, Tobias Lachenmaier, Rupert Leitner, Chao Li, Jiaxing Li, Weidong Li, Weiguo Li, Xiaomei Li, Xiaonan Li, Yi Li, Yufeng Li, Zhi-Bing Li, Hao Liang, Guey-Lin Lin, Tao Lin, Yen-Hsun Lin, Jiajie Ling, Ivano Lippi, Dawei Liu, Hongbang Liu, Hu Liu, Jianglai Liu, Jianli Liu, Jinchang Liu, Qian Liu, Shubin Liu, Shulin Liu, Paolo Lombardi, Yongbing Long, Haoqi Lu, Jiashu Lu, Jingbin Lu, Junguang Lu, Bayarto Lubsandorzhiev, Livia Ludhova, Shu Luo, Vladimir Lyashuk, Randolph Mollenberg, Xubo Ma, Fabio Mantovani, Yajun Mao, Stefano M. Mari, William F. McDonough, Guang Meng, Anselmo Meregaglia, Emanuela Meroni, Mauro Mezzetto, Lino Miramonti, Thomas Mueller, Dmitry Naumov, Lothar Oberauer, Juan Pedro Ochoa-Ricoux, Alexander Olshevskiy, Fausto Ortica, Alessandro Paoloni, Haiping Peng, Jen-Chieh Peng, Ezio Previtali, Ming Qi, Sen Qian, Xin Qian, Yongzhong Qian, Zhonghua Qin, Georg Raffelt, Gioacchino Ranucci, Barbara Ricci, Markus Robens, Aldo Romani, Xiangdong Ruan, Xichao Ruan, Giuseppe Salamanna, Mike Shaevitz, Valery Sinev, Chiara Sirignano, Monica Sisti, Oleg Smirnov, Michael Soiron, Achim Stahl, Luca Stanco, Jochen Steinmann, Xilei Sun, Yongjie Sun, Dmitriy Taichenachev, Jian Tang, Igor Tkachev, Wladyslaw Trzaska, Stefan van Waasen, Cristina Volpe, Vit Vorobel, Lucia Votano, Chung-Hsiang Wang, Guoli Wang, Hao Wang, Meng Wang, Ruiguang Wang, Siguang Wang, Wei Wang, Yi Wang, Yi Wang, Yifang Wang, Zhe Wang, Zheng Wang, Zhigang Wang, Zhimin Wang, Wei Wei, Liangjian Wen, Christopher Wiebusch, Bjorn Wonsak, Qun Wu, Claudia-Elisabeth Wulz, Michael Wurm, Yufei Xi, Dongmei Xia, Yuguang Xie, Zhi-zhong Xing, Jilei Xu, Baojun Yan, Changgen Yang, Chaowen Yang, Guang Yang, Lei Yang, Yifan Yang, Yu Yao, Ugur Yegin, Frederic Yermia, Zhengyun You, Boxiang Yu, Chunxu Yu, Zeyuan Yu, Sandra Zavatarelli, Liang Zhan, Chao Zhang, Hong-Hao Zhang, Jiawen Zhang, Jingbo Zhang, Qingmin Zhang, Yu-Mei Zhang, Zhenyu Zhang, Zhenghua Zhao, Yangheng Zheng, Weili Zhong, Guorong Zhou, Jing Zhou, Li Zhou, Rong Zhou, Shun Zhou, Wenxiong Zhou, Xiang Zhou, Yeling Zhou, Yufeng Zhou, Jiaheng Zou
    Oct. 18, 2015 hep-ex, physics.ins-det
    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.
  • 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 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 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.
  • 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 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.
  • 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.
  • 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 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.
  • 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.