Robust nonadiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of nonadiabatic electronic transitions and test it against exact integration of the time-dependent Schrodinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the nonadiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.

Quantum annealing of an ising spin-glass by Green's function Monte Carlo

We present an implementation of quantum annealing (QA) via lattice Green's function Monte Carlo (GFMC), focusing on its application to the Ising spin glass in transverse field. In particular, we study whether or not such a method is more effective than the path-integral Monte Carlo- (PIMC) based QA, as well as classical simulated annealing (CA), previously tested on the same optimization problem. We identify the issue of importance sampling, i.e., the necessity of possessing reasonably good (variational) trial wave functions, as the key point of the algorithm. We performed GFMC-QA runs using such a Boltzmann-type trial wave function, finding results for the residual energies that are qualitatively similar to those of CA (but at a much larger computational cost), and definitely worse than PIMC-QA. We conclude that, at present, without a serious effort in constructing reliable importance sampling variational wave functions for a quantum glass, GFMC-QA is not a true competitor of PIMC-QA.

Optimization through quantum annealing: theory and some applications

Quantum annealing is a promising tool for solving optimization problems, similar in some ways to the traditional ( classical) simulated annealing of Kirkpatrick et al. Simulated annealing takes advantage of thermal fluctuations in order to explore the optimization landscape of the problem at hand, whereas quantum annealing employs quantum fluctuations. Intriguingly, quantum annealing has been proved to be more effective than its classical counterpart in many applications. We illustrate the theory and the practical implementation of both classical and quantum annealing - highlighting the crucial differences between these two methods - by means of results recently obtained in experiments, in simple toy-models, and more challenging combinatorial optimization problems ( namely, Random Ising model and Travelling Salesman Problem). The techniques used to implement quantum and classical annealing are either deterministic evolutions, for the simplest models, or Monte Carlo approaches, for harder optimization tasks. We discuss the pro and cons of these approaches and their possible connections to the landscape of the problem addressed.

Monte Carlo studies of quantum and classical annealing on a double well

We present results for a variety of Monte Carlo annealing approaches, both classical and quantum, benchmarked against one another for the textbook optimization exercise of a simple one-dimensional double well. In classical (thermal) annealing, the dependence upon the move chosen in a Metropolis scheme is studied and correlated with the spectrum of the associated Markov transition matrix. In quantum annealing, the path integral Monte Carlo approach is found to yield nontrivial sampling difficulties associated with the tunneling between the two wells. The choice of fictitious quantum kinetic energy is also addressed. We find that a "relativistic" kinetic energy form, leading to a higher probability of long real-space jumps, can be considerably more effective than the standard nonrelativistic one.

Optimization by quantum annealing: Lessons from simple cases

We investigate the basic behavior and performance of simulated quantum annealing (QA) in comparison with classical annealing (CA). Three simple one-dimensional case study systems are considered: namely, a parabolic well, a double well, and a curved washboard. The time-dependent Schrodinger evolution in either real or imaginary time describing QA is contrasted with the Fokker-Planck evolution of CA. The asymptotic decrease of excess energy with annealing time is studied in each case, and the reasons for differences are examined and discussed. The Huse-Fisher classical power law of double-well CA is replaced with a different power law in QA. The multiwell washboard problem studied in CA by Shinomoto and Kabashima and leading classically to a logarithmic annealing even in the absence of disorder turns to a power-law behavior when annealed with QA. The crucial role of disorder and localization is briefly discussed.

Theory of phonon dissipation in the conduction of stressed Au nanowires

Recent experiments on Au break junctions [Phys. Rev. Lett. 88 (2002) 216803] have characterized the nonlinear conductance of stretched short Au nanowires. They reveal in the voltage range 10-20 meV the signatures of dissipation effects, likely due to phonons in the nanowire, reducing the conductance below the quantized value of 2e(2)/h. We present here a theory, based on a model tight-binding Hamiltonian and on non-equilibrium Green's function techniques, which accounts for the main features of the experiment. The theory helps in revealing details of the experiment which need to be addressed with a more realistic, less idealized, theoretical framework. (C) 2004 Elsevier B.V. All rights reserved.

On transition rates in surface hopping

Trajectory surface hopping (TSH) is one of the most widely used quantum-classical algorithms for nonadiabatic molecular dynamics. Despite its empirical effectiveness and popularity, a rigorous derivation of TSH as the classical limit of a combined quantum electron-nuclear dynamics is still missing. In this work, we aim to elucidate the theoretical basis for the widely used hopping rules. Naturally, we concentrate thereby on the formal aspects of the TSH. Using a Gaussian wave packet limit, we derive the transition rates governing the hopping process at a simple avoided level crossing. In this derivation, which gives insight into the physics underlying the hopping process, some essential features of the standard TSH algorithm are retrieved, namely (i) non-zero electronic transition rate ("hopping probability") at avoided crossings; (ii) rescaling of the nuclear velocities to conserve total energy; (iii) electronic transition rates linear in the nonadiabatic coupling vectors. The well-known Landau-Zener model is then used for illustration. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4770280]

Predictors and dynamics of postpartum relapses in women with multiple sclerosis

Several studies have shown that pregnancy reduces multiple sclerosis (MS) relapses, which increase in the early postpartum period. Postpartum relapse risk has been predicted by pre-pregnancy disease activity in some studies.

Structure and thermal properties of salicylate-based-protic ionic liquids as new heat storage media. COSMO-RS structure characterization and modeling of heat capacities

During this research, we present a study on the thermal properties, such as the melting, cold crystallization, and glass transition temperatures as well as heat capacities from 293.15 K to 323.15 K of nine in-house synthesized protic ionic liquids based on the 3-(alkoxymethyl)-1H-imidazol-3-ium salicylate ([H-Im-C₁OC_n][Sal]) with n = 3–11. The 3D structures, surface charge distributions and COSMO volumes of all investigated ions are obtained by combining DFT calculations and the COSMO-RS methodology. The heat capacity data sets as a function of temperature of the 3-(alkoxymethyl)-1H-imidazol-3-ium salicylate are then predicted using the methodology originally proposed in the case of ionic liquids by Ge et al. 3-(Alkoxymethyl)-1H-imidazol-3-ium salicylate based ionic liquids present specific heat capacities higher in many cases than other ionic liquids that make them suitable as heat storage media and in heat transfer processes. It was found experimentally that the heat capacity increases linearly with increasing alkyl chain length of the alkoxymethyl group of 3-(alkoxymethyl)-1H-imidazol-3-ium salicylate as was expected and predicted using the Ge et al. method with an overall relative absolute deviation close to 3.2% for temperatures up to 323.15 K.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated 1/42,000, 1/43,700 and 1/49,500 SNPs explained 1/421%, 1/424% and 1/429% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/I 2-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Tiotropium and olodaterol fixed-dose combination versus mono-components in COPD (GOLD 2-4)

Efficacy and safety of tiotropium+olodaterol fixed-dose combination (FDC) compared with the mono-components was evaluated in patients with moderate to very severe chronic obstructive pulmonary disease (COPD) in two replicate, randomised, double-blind, parallel-group, multicentre, phase III trials. Patients received tiotropium+olodaterol FDC 2.5/5 μg or 5/5 μg, tiotropium 2.5 μg or 5 μg, or olodaterol 5 μg delivered once-daily via Respimat inhaler over 52 weeks. Primary end points were forced expiratory volume in 1 s (FEV1) area under the curve from 0 to 3 h (AUC0-3) response, trough FEV1 response and St George's Respiratory Questionnaire (SGRQ) total score at 24 weeks. In total, 5162 patients (2624 in Study 1237.5 and 2538 in Study 1237.6) received treatment. Both FDCs significantly improved FEV1 AUC0-3 and trough FEV1 response versus the mono-components in both studies. Statistically significant improvements in SGRQ total score versus the mono-components were only seen for tiotropium+olodaterol FDC 5/5 μg. Incidence of adverse events was comparable between the FDCs and the mono-components. These studies demonstrated significant improvements in lung function and health-related quality of life with once-daily tiotropium+olodaterol FDC versus mono-components over 1 year in patients with moderate to very severe COPD.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ∼2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Can the scaling behavior of electric conductivity be used to probe the self-organizational changes in solution with respect to the ionic liquid structure? The case of [C8MIM][NTf2]

Electrical conductivity of the supercooled ionic liquid [C8MIM][NTf2], determined as a function of temperature and pressure, highlights strong differences in its ionic transport behavior between low and high temperature regions. To date, the crossover effect which is very well known for low molecular van der Waals liquids has been rarely described for classical ionic liquids. This finding highlights that the thermal fluctuations could be dominant mechanisms driving the dramatic slowing down of ion motions near Tg. An alternative way to analyze separately low and high temperature dc-conductivity data using a density scaling approach was then proposed. Based on which a common value of the scaling exponent [gamma] = 2.4 was obtained, indicating that the applied density scaling is insensitive to the crossover effect. By comparing the scaling exponent [gamma] reported herein along with literature data for other ionic liquids, it appears that [gamma] decreases by increasing the alkyl chain length on the 1-alkyl-3-methylimidazolium-based ionic liquids. This observation may be related to changes in the interaction between ions in solution driven by an increase in the van der Waals type interaction by increasing the alkyl chain length on the cation. This effect may be related to changes in the ionic liquid nanostructural organization with the alkyl chain length on the cation as previously reported in the literature based on molecular dynamic simulations. In other words, the calculated scaling exponent [gamma] may be then used as a key parameter to probe the interaction and/or self-organizational changes in solution with respect to the ionic liquid structure.