Reactive Force Field for Proton Diffusion in BaZrO3 using an empirical valence bond approach

A new reactive force field to describe proton diffusion within the solid-oxide fuel cell material BaZrO3 has been derived. Using a quantum mechanical potential energy surface, the parameters of an interatomic potential model to describe hydroxyl groups within both pure and yttrium-doped BaZrO3 have been determined. Reactivity is then incorporated through the use of the empirical valence bond model. Molecular dynamics simulations (EVB-MD) have been performed to explore the diffusion of hydrogen using a stochastic thermostat and barostat whose equations are extended to the isostress-isothermal ensemble. In the low concentration limit, the presence of yttrium is found not to significantly influence the diffusivity of hydrogen, despite the proton having a longer residence time at oxygen adjacent to the dopant. This lack of influence is due to the fact that trapping occurs infrequently, even when the proton diffuses through octahedra adjacent to the dopant. The activation energy for diffusion is found to be 0.42 eV, in good agreement with experimental values, though the prefactor is slightly underestimated.Comment: Corrected titl

Canonical sampling through velocity-rescaling

We present a new molecular dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains constant during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. We illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liquid phases. Its performance is excellent and largely independent on the thermostat parameter also with regard to the dynamic properties

Molecular Simulations Matching Denaturation Experiments for N-6-Methyladenosine

Post-transcriptional modifications are crucial for RNA function and can affect its structure and dynamics. Force-field-based classical molecular dynamics simulations are a fundamental tool to characterize biomolecular dynamics, and their application to RNA is flourishing. Here, we show that the set of force-field parameters for N-6-methyladenosine (m(6)A) developed for the commonly used AMBER force field does not reproduce duplex denaturation experiments and, specifically, cannot be used to describe both paired and unpaired states. Then, we use reweighting techniques to derive new parameters matching available experimental data. The resulting force field can be used to properly describe paired and unpaired m(6)A in both syn and anti conformation, which thus opens the way to the use of molecular simulations to investigate the effects of N-6 methylations on RNA structural dynamics

Does Ménière's Disease in the Elderly Present Some Peculiar Features?

Object. Aim of our study was to establish some peculiar features of Ménière's Disease (MD) in a group of elderly MD patients, in which the first vertigo spell happened when over 65 years old. Material and Methods. We analyzed a group of 73 younger than 65-years-old and a group of 30 elderly MD patients. All patients underwent a neurotological evaluation, an anamnestic evaluation including a lifetime history of migraine, and blood withdrawal for autoantibody screening. Results. Some differences were found between elderly and younger MD patients. Elderly MD patients presented a higher prevalence of Tumarkin attacks and a lower prevalence of lifetime history of migraine; moreover, they presented a faster develop of hearing loss and vertigo spells than a subgroup of 32 younger patients matched for the duration of illness. Conclusions. Some clinical features of MD in elderly have been pointed out. Particularly, the lower rate of migrainous history and positivity for autoantibodies often associated with MD, in our opinion, support the hypothesis of a vascular disorder acting as a predisposing factor for MD in elderly

Molecular Dynamics of Solids at Constant Pressure and Stress Using Anisotropic Stochastic Cell Rescaling

Molecular dynamics simulations of solids are often performed using anisotropic barostats that allow the shape and volume of the periodic cell to change during the simulation. Most existing schemes are based on a second-order differential equation that might lead to undesired oscillatory behaviors and should not be used in the equilibration phase. We recently introduced stochastic cell rescaling, a first-order stochastic barostat that can be used for both the equilibration and production phases. Only the isotropic and semi-isotropic variants have been formulated and implemented so far. In this paper, we develop and implement the equations of motion of the fully anisotropic variant and test them on Lennard-Jones solids, ice, gypsum, and gold. The algorithm has a single parameter that controls the relaxation time of the volume, results in the exponential decay of correlation functions, and can be effectively applied to a wide range of systems

Accurate sampling using Langevin dynamics

We show how to derive a simple integrator for the Langevin equation and illustrate how it is possible to check the accuracy of the obtained distribution on the fly, using the concept of effective energy introduced in a recent paper [J. Chem. Phys. 126, 014101 (2007)]. Our integrator leads to correct sampling also in the difficult high-friction limit. We also show how these ideas can be applied in practical simulations, using a Lennard-Jones crystal as a paradigmatic case

Colored-noise thermostats \`a la carte

Recently, we have shown how a colored-noise Langevin equation can be used in the context of molecular dynamics as a tool to obtain dynamical trajectories whose properties are tailored to display desired sampling features. In the present paper, after having reviewed some analytical results for the stochastic differential equations forming the basis of our approach, we describe in detail the implementation of the generalized Langevin equation thermostat and the fitting procedure used to obtain optimal parameters. We discuss in detail the simulation of nuclear quantum effects, and demonstrate that, by carefully choosing parameters, one can successfully model strongly anharmonic solids such as neon. For the reader's convenience, a library of thermostat parameters and some demonstrative code can be downloaded from an on-line repository