D-optimal designs via a cocktail algorithm

A fast new algorithm is proposed for numerical computation of (approximate) D-optimal designs. This "cocktail algorithm" extends the well-known vertex direction method (VDM; Fedorov 1972) and the multiplicative algorithm (Silvey, Titterington and Torsney, 1978), and shares their simplicity and monotonic convergence properties. Numerical examples show that the cocktail algorithm can lead to dramatically improved speed, sometimes by orders of magnitude, relative to either the multiplicative algorithm or the vertex exchange method (a variant of VDM). Key to the improved speed is a new nearest neighbor exchange strategy, which acts locally and complements the global effect of the multiplicative algorithm. Possible extensions to related problems such as nonparametric maximum likelihood estimation are mentioned.Comment: A number of changes after accounting for the referees' comments including new examples in Section 4 and more detailed explanations throughou

First-principles study of the ferroelastic phase transition in CaCl_2

First-principles density-functional calculations within the local-density approximation and the pseudopotential approach are used to study and characterize the ferroelastic phase transition in calcium chloride (CaCl_2). In accord with experiment, the energy map of CaCl_2 has the typical features of a pseudoproper ferroelastic with an optical instability as ultimate origin of the phase transition. This unstable optic mode is close to a pure rigid unit mode of the framework of chlorine atoms and has a negative Gruneisen parameter. The ab-initio ground state agrees fairly well with the experimental low temperature structure extrapolated at 0K. The calculated energy map around the ground state is interpreted as an extrapolated Landau free-energy and is successfully used to explain some of the observed thermal properties. Higher-order anharmonic couplings between the strain and the unstable optic mode, proposed in previous literature as important terms to explain the soft-phonon temperature behavior, are shown to be irrelevant for this purpose. The LAPW method is shown to reproduce the plane-wave results in CaCl_2 within the precision of the calculations, and is used to analyze the relative stability of different phases in CaCl_2 and the chemically similar compound SrCl_2.Comment: 9 pages, 6 figures, uses RevTeX

Measurements of long-range near-side angular correlations in sNN=5\sqrt{s_{\text{NN}}}=5TeV proton-lead collisions in the forward region

Two-particle angular correlations are studied in proton-lead collisions at a nucleon-nucleon centre-of-mass energy of $\sqrt{s_{\text{NN}}}=5$TeV, collected with the LHCb detector at the LHC. The analysis is based on data recorded in two beam configurations, in which either the direction of the proton or that of the lead ion is analysed. The correlations are measured in the laboratory system as a function of relative pseudorapidity, $\Delta\eta$, and relative azimuthal angle, $\Delta\phi$, for events in different classes of event activity and for different bins of particle transverse momentum. In high-activity events a long-range correlation on the near side, $\Delta\phi \approx 0$, is observed in the pseudorapidity range $2.0<\eta<4.9$. This measurement of long-range correlations on the near side in proton-lead collisions extends previous observations into the forward region up to $\eta=4.9$. The correlation increases with growing event activity and is found to be more pronounced in the direction of the lead beam. However, the correlation in the direction of the lead and proton beams are found to be compatible when comparing events with similar absolute activity in the direction analysed.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-040.htm

Study of the production of Λb0\Lambda_b^0 and B‾0\overline{B}^0 hadrons in pppp collisions and first measurement of the Λb0→J/ψpK−\Lambda_b^0\rightarrow J/\psi pK^- branching fraction

The product of the $\Lambda_b^0$ ($\overline{B}^0$) differential production cross-section and the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ ($\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$) is measured as a function of the beauty hadron transverse momentum, $p_{\rm T}$, and rapidity, $y$. The kinematic region of the measurements is $p_{\rm T}<20~{\rm GeV}/c$ and $2.0<y<4.5$. The measurements use a data sample corresponding to an integrated luminosity of $3~{\rm fb}^{-1}$ collected by the LHCb detector in $pp$ collisions at centre-of-mass energies $\sqrt{s}=7~{\rm TeV}$ in 2011 and $\sqrt{s}=8~{\rm TeV}$ in 2012. Based on previous LHCb results of the fragmentation fraction ratio, $f_{\Lambda_B^0}/f_d$, the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ is measured to be \begin{equation*} \mathcal{B}(\Lambda_b^0\rightarrow J/\psi pK^-)= (3.17\pm0.04\pm0.07\pm0.34^{+0.45}_{-0.28})\times10^{-4}, \end{equation*} where the first uncertainty is statistical, the second is systematic, the third is due to the uncertainty on the branching fraction of the decay $\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$, and the fourth is due to the knowledge of $f_{\Lambda_b^0}/f_d$. The sum of the asymmetries in the production and decay between $\Lambda_b^0$ and $\overline{\Lambda}_b^0$ is also measured as a function of $p_{\rm T}$ and $y$. The previously published branching fraction of $\Lambda_b^0\rightarrow J/\psi p\pi^-$, relative to that of $\Lambda_b^0\rightarrow J/\psi pK^-$, is updated. The branching fractions of $\Lambda_b^0\rightarrow P_c^+(\rightarrow J/\psi p)K^-$ are determined.Comment: 29 pages, 19figures. All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-032.htm

BB flavour tagging using charm decays at the LHCb experiment

An algorithm is described for tagging the flavour content at production of neutral $B$ mesons in the LHCb experiment. The algorithm exploits the correlation of the flavour of a $B$ meson with the charge of a reconstructed secondary charm hadron from the decay of the other $b$ hadron produced in the proton-proton collision. Charm hadron candidates are identified in a number of fully or partially reconstructed Cabibbo-favoured decay modes. The algorithm is calibrated on the self-tagged decay modes $B^+ \to J/\psi \, K^+$ and $B^0 \to J/\psi \, K^{*0}$ using $3.0\mathrm{\,fb}^{-1}$ of data collected by the LHCb experiment at $pp$ centre-of-mass energies of $7\mathrm{\,TeV}$ and $8\mathrm{\,TeV}$. Its tagging power on these samples of $B \to J/\psi \, X$ decays is $(0.30 \pm 0.01 \pm 0.01) \%$.Comment: All figures and tables, along with any supplementary material and additional information, are available at http://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-027.htm

Evidence for the strangeness-changing weak decay Ξb−→Λb0π−\Xi_b^-\to\Lambda_b^0\pi^-

Using a $pp$ collision data sample corresponding to an integrated luminosity of 3.0~fb$^{-1}$, collected by the LHCb detector, we present the first search for the strangeness-changing weak decay $\Xi_b^-\to\Lambda_b^0\pi^-$. No $b$ hadron decay of this type has been seen before. A signal for this decay, corresponding to a significance of 3.2 standard deviations, is reported. The relative rate is measured to be ${{f_{\Xi_b^-}}\over{f_{\Lambda_b^0}}}{\cal{B}}(\Xi_b^-\to\Lambda_b^0\pi^-) = (5.7\pm1.8^{+0.8}_{-0.9})\times10^{-4}$, where $f_{\Xi_b^-}$ and $f_{\Lambda_b^0}$ are the $b\to\Xi_b^-$ and $b\to\Lambda_b^0$ fragmentation fractions, and ${\cal{B}}(\Xi_b^-\to\Lambda_b^0\pi^-)$ is the branching fraction. Assuming $f_{\Xi_b^-}/f_{\Lambda_b^0}$ is bounded between 0.1 and 0.3, the branching fraction ${\cal{B}}(\Xi_b^-\to\Lambda_b^0\pi^-)$ would lie in the range from $(0.57\pm0.21)\%$ to $(0.19\pm0.07)\%$.Comment: 7 pages, 2 figures, All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-047.htm

Lichenometric dating (lichenometry) and the biology of the lichen genus rhizocarpon:challenges and future directions

Lichenometric dating (lichenometry) involves the use of lichen measurements to estimate the age of exposure of various substrata. Because of low radial growth rates and considerable longevity, species of the crustose lichen genus Rhizocarpon have been the most useful in lichenometry. The primary assumption of lichenometry is that colonization, growth and mortality of Rhizocarpon are similar on surfaces of known and unknown age so that the largest thalli present on the respective faces are of comparable age. This review describes the current state of knowledge regarding the biology of Rhizocarpon and considers two main questions: (1) to what extent does existing knowledge support this assumption; and (2) what further biological observations would be useful both to test its validity and to improve the accuracy of lichenometric dates? A review of the Rhizocarpon literature identified gaps in knowledge regarding early development, the growth rate/size curve, mortality, regeneration, competitive effects, colonization, and succession on rock surfaces. The data suggest that these processes may not be comparable on different rock surfaces, especially in regions where growth rates and thallus turnover are high. In addition, several variables could differ between rock surfaces and influence maximum thallus size, including rate and timing of colonization, radial growth rates, environmental differences, thallus fusion, allelopathy, thallus mortality, colonization and competition. Comparative measurements of these variables on surfaces of known and unknown age may help to determine whether the basic assumptions of lichenometry are valid. Ultimately, it may be possible to take these differences into account when interpreting estimated dates

Thermodynamic model of hardness: Particular case of boron-rich solids

A number of successful theoretical models of hardness have been developed recently. A thermodynamic model of hardness, which supposes the intrinsic character of correlation between hardness and thermodynamic properties of solids, allows one to predict hardness of known or even hypothetical solids from the data on Gibbs energy of atomization of the elements, which implicitly determine the energy density per chemical bonding. The only structural data needed is the coordination number of the atoms in a lattice. Using this approach, the hardness of known and hypothetical polymorphs of pure boron and a number of boron-rich solids has been calculated. The thermodynamic interpretation of the bonding energy allows one to predict the hardness as a function of thermodynamic parameters. In particular, the excellent agreement between experimental and calculated values has been observed not only for the room- temperature values of the Vickers hardness of stoichiometric compounds, but also for its temperature and concentration dependencies