BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michele Parrinello

Michele Parrinello appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.comp-phcond-mat.stat-mechphysics.chem-phcond-mat.mtrl-scicond-mat.softBiological Physicsquant-phBiomoleculesComputational Engineering, Finance, and Sciencecond-mat.quant-gasphysics.data-anphysics.geo-ph

Catalog footprint

What is connected

42works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Exploration vs Convergence Speed in Adaptive-bias Enhanced Sampling

In adaptive-bias enhanced sampling methods, a bias potential is added to the system to drive transitions between metastable states. The bias potential is a function of a few collective variables and is gradually modified according to the underlying free energy surface. We show that when the collective variables are suboptimal, there is an exploration-convergence tradeoff, and one must choose between a quickly converging bias that will lead to fewer transitions, or a slower to converge bias that can explore the phase space more efficiently but might require a much longer time to produce an accurate free energy estimate. The recently proposed On-the-fly Probability Enhanced Sampling (OPES) method focuses on fast convergence, but there are cases where fast exploration is preferred instead. For this reason, we introduce a new variant of the OPES method that focuses on quickly escaping metastable states, at the expense of convergence speed. We illustrate the benefits of this approach on prototypical systems and show that it outperforms the popular metadynamics method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kinetics of drug release from clay using enhanced sampling methods

A key step in the development of a new drug is the design of drug-excipient complexes that lead to an optimal drug release kinetics. Computational chemistry, specifically enhanced sampling molecular dynamics methods, can play a key role in this context by minimizing the need for expensive experiments and reducing cost and time. Here, we show that recent advances in enhanced sampling methodologies can be brought to fruition in this area. We demonstrate the potential of these methodologies by simulating the release kinetics of the complex praziquantel-montmorillonite in water. Praziquantel finds promising applications in the treatment of schistosomiasis, but its biopharmaceutical profile needs to be improved and a cheap material like the montmorillonite clay would be a very convenient excipient. We simulate the drug release both from surface and interlayer space and find that the time it takes to the drug to diffuse from the interlayer space to the solution controls the release time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Solubility of organic salts in solvent-antisolvent mixtures: A combined experimental and molecular dynamics simulations approach

We combine molecular dynamics simulations with experiments to estimate solubilities of organic salts in complex growth environments. We predict the solubility by simulations of the growth and dissolution of ions at the crystal surface kink sites at different solution concentrations. Thereby, the solubility is identified as the solution's salt concentration, where the energy of the ion pair dissolved in solution equals the energy of the ion pair crystallized at the kink sites. The simulation methodology is demonstrated for the case of anhydrous sodium acetate crystallized from various solvent-antisolvent mixtures. To validate the predicted solubilities, we have measured the solubilities of sodium acetate in-house, using an experimental setup and measurement protocol that guarantees moisture-free conditions, which is key for a hygroscopic compound like sodium acetate. We observe excellent agreement between the experimental and the computationally evaluated solubilities for sodium acetate in different solvent-antisolvent mixtures. Given the agreement and the rich data the simulations produce, we can use them to complement experimental tasks which in turn will reduce time and capital in the design of complicated industrial crystallization processes of organic salts.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Deep learning the slow modes for rare events sampling

The development of enhanced sampling methods has greatly extended the scope of atomistic simulations, allowing long-time phenomena to be studied with accessible computational resources. Many such methods rely on the identification of an appropriate set of collective variables. These are meant to describe the system's modes that most slowly approach equilibrium. Once identified, the equilibration of these modes is accelerated by the enhanced sampling method of choice. An attractive way of determining the collective variables is to relate them to the eigenfunctions and eigenvalues of the transfer operator. Unfortunately, this requires knowing the long-term dynamics of the system beforehand, which is generally not available. However, we have recently shown that it is indeed possible to determine efficient collective variables starting from biased simulations. In this paper, we bring the power of machine learning and the efficiency of the recently developed on-the-fly probability enhanced sampling method to bear on this approach. The result is a powerful and robust algorithm that, given an initial enhanced sampling simulation performed with trial collective variables or generalized ensembles, extracts transfer operator eigenfunctions using a neural network ansatz and then accelerates them to promote sampling of rare events. To illustrate the generality of this approach we apply it to several systems, ranging from the conformational transition of a small molecule to the folding of a mini-protein and the study of materials crystallization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Prediction of a supersolid phase in high-pressure deuterium

Supersolid is a mysterious and puzzling state of matter whose possible existence has stirred a vigorous debate among physicists for over 60 years. Its elusive nature stems from the coexistence of two seemingly contradicting properties, long-range order and superfluidity. We report computational evidence of a supersolid phase of deuterium under high pressure ($p >800$ GPa) and low temperature (T $<$ 1.0 K). In our simulations, that are based on bosonic path integral molecular dynamics, we observe a highly concerted exchange of atoms while the system preserves its crystalline order. The exchange processes are favoured by the soft core interactions between deuterium atoms that form a densely packed metallic solid. At the zero temperature limit, Bose-Einstein condensation is observed as the permutation probability of $N$ deuterium atoms approaches $1/N$ with a finite superfluid fraction. Our study provides concrete evidence for the existence of a supersolid phase in high-pressure deuterium and could provide insights on the future investigation of supersolid phases in real materials.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The role of water in host-guest interaction

One of the main applications of atomistic computer simulations is the calculation of ligand binding energies. The accuracy of these calculations depends on the force field quality and on the thoroughness of configuration sampling. Sampling is an obstacle in modern simulations due to the frequent appearance of kinetic bottlenecks in the free energy landscape. Very often this difficulty is circumvented by enhanced sampling techniques. Typically, these techniques depend on the introduction of appropriate collective variables that are meant to capture the system's degrees of freedom. In ligand binding, water has long been known to play a key role, but its complex behaviour has proven difficult to fully capture. In this paper we combine machine learning with physical intuition to build a non-local and highly efficient water-describing collective variable. We use it to study a set of of host-guest systems from the SAMPL5 challenge. We obtain highly accurate binding energies and good agreement with experiments. The role of water during the binding process is then analysed in some detail.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-Driven Collective Variables for Enhanced Sampling

Designing an appropriate set of collective variables is crucial to the success of several enhanced sampling methods. Here we focus on how to obtain such variables from information limited to the metastable states. We characterize these states by a large set of descriptors and employ neural networks to compress this information in a lower-dimensional space, using Fisher's linear discriminant as an objective function to maximize the discriminative power of the network. We test this method on alanine dipeptide, using the non-linearly separable dataset composed by atomic distances. We then study an intermolecular aldol reaction characterized by a concerted mechanism. The resulting variables are able to promote sampling by drawing non-linear paths in the physical space connecting the fluctuations between metastable basins. Lastly, we interpret the behavior of the neural network by studying its relation to the physical variables. Through the identification of its most relevant features, we are able to gain chemical insight into the process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How collective phenomena impact CO2 reactivity and speciation in different media

CO2 has attracted considerable attention in the recent years due to its role in the greenhouse effect and environmental management. While its reaction with water has been studied extensively, the same cannot be said for reactivity in supercritical CO2 phase, where the conjugate acid/base equilibria proceed through different mechanisms and activation barriers. In spite of the apparent simplicity of the CO2 + H2O reaction, the collective effect of different environments has drastic influence on the free energy profile. Enhanced sampling techniques and well-tailored collective variables provide a detailed picture of the enthalpic and entropic drivers underscoring the differences in the formation mechanism of carbonic acid in the gas, aqueous and supercritical CO2 phases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Metadynamics of paths

We present a method to sample reactive pathways via biased molecular dynamics simulations in trajectory space. We show that the use of enhanced sampling techniques enables unconstrained exploration of multiple reaction routes. Time correlation functions are conveniently computed via reweighted averages along a single trajectory and kinetic rates are accessed at no additional cost. These abilities are illustrated analyzing a model potential and the umbrella inversion of NH$_3$ in water. The algorithm allows a parallel implementation and promises to be a powerful tool for the study of rare events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations

Accurate prediction of a gas solubility in a liquid is crucial in many areas of chemistry, and a detailed understanding of the molecular mechanism of the gas solvation continues to be an active area of research. Here, we extend the idea of constant chemical potential molecular dynamics (CμMD) approach to the calculation of the gas solubility in the liquid under constant gas chemical potential conditions. As a representative example, we utilize this method to calculate the isothermal solubility of carbon dioxide in water. Additionally, we provide microscopic insight into the mechanism of solvation that preferentially occurs in areas of the surface where the hydrogen network is broken.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality

We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to 7 electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rethinking Metadynamics: from bias potentials to probability distributions

Metadynamics is an enhanced sampling method of great popularity, based on the on-the-fly construction of a bias potential that is function of a selected number of collective variables. We propose here a change in perspective that shifts the focus from the bias to the probability distribution reconstruction, while keeping some of the key characteristics of metadynamics, such as the flexible on-the-fly adjustments to the free energy estimate. The result is an enhanced sampling method that presents a drastic improvement in convergence speed, especially when dealing with suboptimal and/or multidimensional sets of collective variables. The method is especially robust and easy to use, in fact it requires only few simple parameters to be set, and it has a straightforward reweighting scheme to recover the statistics of the unbiased ensemble. Furthermore it gives more control on the desired exploration of the phase space, since the deposited bias is not allowed to grow indefinitely and it does not push the simulation to uninteresting high free energy regions. We demonstrate the performance of the method in a number of representative examples.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A microscopic description of acid-base equilibrium

Acid-base reactions are ubiquitous in nature. Understanding their mechanisms is crucial in many fields, from biochemistry to industrial catalysis. Unfortunately, experiments only give limited information without much insight into the molecular behaviour. Atomistic simulations could complement experiments and shed precious light on microscopic mechanisms. The large free energy barriers connected to proton dissociation however make the use of enhanced sampling methods mandatory. Here we perform an ab initio molecular dynamics (MD) simulation and enhance sampling with the help of methadynamics. This has been made possible by the introduction of novel descriptors or collective variables (CVs) that are based on a conceptually new outlook on acid-base equilibria. We test successfully our approach on three different aqueous solutions of acetic acid, ammonia, and bicarbonate. These are representative of acid, basic, and amphoteric behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Calculation of phase diagrams in the multithermal-multibaric ensemble

From the Ising model and the Lennard-Jones fluid, to water and the iron-carbon system, phase diagrams are an indispensable tool to understand phase equilibria. In spite of the effort of the simulation community the calculation of a large portion of a phase diagram using computer simulation is still today a significant challenge. Here we propose a method to calculate phase diagrams involving liquid and solid phases by the reversible transformation of the liquid and the solid. To this end we introduce an order parameter that breaks the rotational symmetry and we leverage our recently introduced method to sample the multithermal-multibaric ensemble. In this way in a single molecular dynamics simulation we are able to compute the liquid-solid coexistence line for entire regions of the temperature and pressure phase diagram. We apply our approach to the bcc-liquid phase diagram of sodium and the fcc-bcc-liquid phase diagram of aluminum.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Gaussian Mixture Based Enhanced Sampling

Many processes in chemistry and physics take place on timescales that cannot be explored using standard molecular dynamics simulations. This renders the use of enhanced sampling mandatory. Here we introduce an enhanced sampling method that is based on constructing a model probability density from which a bias potential is derived. The model relies on the fact that in a physical system most of the configurations visited can be grouped into isolated metastable islands. To each island we associate a distribution that is fitted to a Gaussian mixture. The different distributions are linearly combined together with coefficients that are computed self consistently. Remarkably, from this biased dynamics, rates of transition between different metastable states can be straightforwardly computed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Path Integral Molecular Dynamics for Bosons

Trapped Bosons exhibit fundamental physical phenomena and are potentially useful for quantum technologies. We present a method for simulating Bosons using path integral molecular dynamics. A main challenge for simulations is including all permutations due to exchange symmetry. We show that evaluation of the potential can be done recursively, avoiding explicit enumeration of permutations, and scales cubically with system size. The method is applied to Bosons in a 2D trap and agrees with essentially exact results. An analysis of the role of exchange with decreasing temperature is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Predicting polymorphism in molecular crystals using orientational entropy

We introduce a computational method to discover polymorphs in molecular crystals at finite temperature. The method is based on reproducing the crystallization process starting from the liquid and letting the system discover the relevant polymorphs. This idea, however, conflicts with the fact that crystallization has a time scale much longer than that of molecular simulations. In order to bring the process within affordable simulation time, we enhance the fluctuations of a collective variable by constructing a bias potential with well tempered metadynamics. We use as collective variable an entropy surrogate based on an extended pair correlation function that includes the correlation between the orientation of pairs of molecules. We also propose a similarity metric between configurations based on the extended pair correlation function and a generalized Kullback-Leibler divergence. In this way, we automatically classify the configurations as belonging to a given polymorph using our metric and a hierarchical clustering algorithm. We find all relevant polymorphs for both substances and we predict new polymorphs. One of them is stabilized at finite temperature by entropic effects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chemical Potential Calculations In Dense Liquids Using Metadynamics

The calculation of chemical potential has traditionally been a challenge in atomistic simulations. One of the most used approaches is Widom's insertion method in which the chemical potential is calculated by periodically attempting to insert an extra particle in the system. In dense systems this method fails since the insertion probability is very low. In this paper we show that in a homogeneous fluid the insertion probability can be increased using metadynamics. We test our method on a supercooled high density binary Lennard-Jones fluid. We find that we can obtain efficiently converged results even when Widom's method fails.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dimer Metadynamics

Sampling complex potential energies is one of the most pressing challenges of contemporary computational science. Inspired by recent efforts that use quantum effects and discretized Feynman's path integrals to overcome large barriers we propose a replica exchange method. In each replica two copies of the same system with halved potential strengths interact via inelastic springs. The strength of the spring is varied in the different replicas so as to bridge the gap between the infinitely strong spring, that corresponds to the Boltzmann replica and the less tight ones. We enhance the spring length fluctuations using Metadynamics. We test the method on simple yet challenging problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Overcoming timescale and finite-size limitations to compute nucleation rates from small scale Well Tempered Metadynamics simulations

Condensation of a liquid droplet from a supersaturated vapour phase is initiated by a prototypical nucleation event. As such it is challenging to compute its rate from atomistic molecular dynamics simulations. In fact at realistic supersaturation conditions condensation occurs on time scales that far exceed what can be reached with conventional molecular dynamics methods. Another known problem in this context is the distortion of the free energy profile associated to nucleation due to the small, finite size of typical simulation boxes. In this work the problem of time scale is addressed with a recently developed enhanced sampling method while contextually correcting for finite size effects. We demonstrate our approach by studying the condensation of argon, and showing that characteristic nucleation times of the order of magnitude of hours can be reliably calculated, approaching realistic supersaturation conditions, thus bridging the gap between what standard molecular dynamics simulations can do and real physical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A perturbative solution to metadynamics ordinary differential equation

Metadynamics is a popular enhanced sampling scheme wherein by periodic application of a repulsive bias, one can surmount high free energy barriers and explore complex landscapes. Recently metadynamics was shown to be mathematically well founded, in the sense that the biasing procedure is guaranteed to converge to the true free energy surface in the long time limit irrespective of the precise choice of biasing parameters. A differential equation governing the post-transient convergence behavior of metadynamics was also derived. In this short communication, we revisit this differential equation, expressing it in a convenient and elegant Riccati-like form. A perturbative solution scheme is then developed for solving this differential equation, which is valid for any generic biasing kernel. The solution clearly demonstrates the robustness of metadynamics to choice of biasing parameters and gives further confidence in the widely used method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular Dynamics Simulations of Solutions at Constant Chemical Potential

Molecular Dynamics studies of chemical processes in solution are of great value in a wide spectrum of applications, which range from nano-technology to pharmaceutical chemistry. However, these calculations are affected by severe finite-size effects, such as the solution being depleted as the chemical process proceeds, which influence the outcome of the simulations. To overcome these limitations, one must allow the system to exchange molecules with a macroscopic reservoir, thus sampling a Grand-Canonical ensemble. Despite the fact that different remedies have been proposed, this still represents a key challenge in molecular simulations. In the present work we propose the Constant Chemical Potential Molecular Dynamics (C$μ$MD) method, which introduces an external force that controls the environment of the chemical process of interest. This external force, drawing molecules from a finite reservoir, maintains the chemical potential constant in the region where the process takes place. We have applied the C$μ$MD method to the paradigmatic case of urea crystallization in aqueous solution. As a result, we have been able to study crystal growth dynamics under constant supersaturation conditions, and to extract growth rates and free-energy barriers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variationally Optimized Free Energy Flooding for Rate Calculation

We propose a new method to obtain kinetic properties of infrequent events from molecular dynamics simulation. The procedure employs a recently introduced variational approach [Valsson and Parrinello, Phys. Rev. Lett. 113, 090601 (2014)] to construct a bias potential as a function of several collective variables that is designed to flood only the associated free energy surface up to a predefined level. The resulting bias potential effectively accelerates transitions between metastable free energy minima while ensuring bias-free transition states, thus allowing accurate kinetic rates to be obtained. We test the method on a few illustrative systems for which we obtain an order of magnitude improvement in efficiency relative to previous approaches, and several orders of magnitude relative to unbiased molecular dynamics. We expect an even larger improvement in more complex systems. This and the ability of the variational approach to deal efficiently with a large number of collective variables will greatly enhance the scope of these calculations. This work is a vindication of the potential that the variational principle has if applied in innovative ways

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Variational Approach to Enhanced Sampling and Free Energy Calculations

The ability of widely used sampling methods, such as molecular dynamics or Monte Carlo, to explore complex free energy landscapes is severely hampered by the presence of kinetic bottlenecks. A large number of solutions have been proposed to alleviate this problem. Many are based on the introduction of a bias potential which is a function of a small number of collective variable. However constructing such a bias is not simple. Here we introduce a functional of the bias potential and an associated variational principle. The bias that minimizes the functional relates in a simple way to the free energy surface. This variational principle can be turned into a practical, efficient and flexible sampling method. A number of numerical examples are presented which include the determination of a three dimensional free energy surface. We argue that, beside being numerically advantageous, our variational approach provides a convenient standpoint for looking with novel eyes at the sampling problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evaluating functions of positive-definite matrices using colored noise thermostats

Many applications in computational science require computing the elements of a function of a large matrix. A commonly used approach is based on the the evaluation of the eigenvalue decomposition, a task that, in general, involves a computing time that scales with the cube of the size of the matrix. We present here a method that can be used to evaluate the elements of a function of a positive-definite matrix with a scaling that is linear for sparse matrices and quadratic in the general case. This methodology is based on the properties of the dynamics of a multidimensional harmonic potential coupled with colored noise generalized Langevin equation (GLE) thermostats. This "$f-$thermostat" (FTH) approach allows us to calculate directly elements of functions of a positive-definite matrix by carefully tailoring the properties of the stochastic dynamics. We demonstrate the scaling and the accuracy of this approach for both dense and sparse problems and compare the results with other established methodologies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From Metadynamics to Dynamics

Metadynamics is a commonly used and successful enhanced sampling method. By the introduction of a history dependent bias which depends on a restricted number of collective variables(CVs) it can explore complex free energy surfaces characterized by several metastable states separated by large free energy barriers. Here we extend its scope by introducing a simple yet powerful method for calculating the rates of transition between different metastable states. The method does not rely on a previous knowledge of the transition states or reaction co-ordinates, as long as CVs are known that can distinguish between the various stable minima in free energy space. We demonstrate that our method recovers the correct escape rates out of these stable states and also preserves the correct sequence of state-to-state transitions, with minimal extra computational effort needed over ordinary metadynamics. We apply the formalism to three different problems and in each case find excellent agreement with the results of long unbiased molecular dynamics runs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accelerating the convergence of path integral dynamics with a generalized Langevin equation

The quantum nature of nuclei plays an important role in the accurate modelling of light atoms such as hydrogen, but it is often neglected in simulations due to the high computational overhead involved. It has recently been shown that zero-point energy effects can be included comparatively cheaply in simulations of harmonic and quasi-harmonic systems by augmenting classical molecular dynamics with a generalized Langevin equation (GLE). Here we describe how a similar approach can be used to accelerate the convergence of path integral (PI) molecular dynamics to the exact quantum mechanical result in more strongly anharmonic systems exhibiting both zero point energy and tunnelling effects. The resulting PI-GLE method is illustrated with applications to a double-well tunnelling problem and to liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Colored-noise thermostats à 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Metadynamics with adaptive Gaussians

Metadynamics is an established sampling method aimed at reconstructing the free-energy surface relative to a set of appropriately chosen collective variables. In standard metadynamics the free-energy surface is filled by the addition of Gaussian potentials of pre-assigned and typically diagonal covariance. Asymptotically the free-energy surface is proportional to the bias deposited. Here we consider the possibility of using Gaussians whose variance is adjusted on the fly to the local properties of the free-energy surface. We suggest two different prescriptions: one is based on the local diffusivity and the other on the local geometrical properties. We further examine the problem of extracting the free-energy surface when using adaptive Gaussians. We show that the standard relation between the bias and the free energy does not hold. In the limit of narrow Gaussians an explicit correction can be evaluated. In the general case we propose to use instead a relation between bias and free energy borrowed from umbrella sampling. This relation holds for all kinds of incrementally deposited bias. We illustrate on the case of alanine dipeptide the advantage of using adaptive Gaussians in conjunction with the new free-energy estimator both in terms of accuracy and speed of convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Microscopic origins of the anomalous melting behaviour of high-pressure sodium

Recent experiments have shown that sodium, a prototype simple metal at ambient conditions, exhibits unexpected complexity under high pressure. One of the most puzzling phenomena in the behaviour of dense sodium is the pressure-induced drop in its melting temperature, which extends from 1000 K at ~30GPa to as low as room temperature at ~120GPa. Despite significant theoretical effort to understand the anomalous melting its origins have remained unclear. In this work, we reconstruct the sodium phase diagram using an ab-initio-quality neural-network potential. We demonstrate that the reentrant behaviour results from the screening of interionic interactions by conduction electrons, which at high pressure induces a softening in the short-range repulsion. It is expected that such an effect plays an important role in governing the behaviour of a wide range of metals and alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Momentum distribution, vibrational dynamics and the potential of mean force in ice

By analyzing the momentum distribution obtained from path integral and phonon calculations we find that the protons in hexagonal ice experience an anisotropic quasi-harmonic effective potential with three distinct principal frequencies that reflect molecular orientation. Due to the importance of anisotropy, anharmonic features of the environment cannot be extracted from existing experimental distributions that involve the spherical average. The full directional distribution is required, and we give a theoretical prediction for this quantity that could be verified in future experiments. Within the quasi-harmonic context, anharmonicity in the ground state dynamics of the proton is substantial and has quantal origin, a finding that impacts the interpretation of several spectroscopies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nucleation mechanism for the direct graphite-to-diamond phase transition

Graphite and diamond have comparable free energies, yet forming diamond from graphite is far from easy. In the absence of a catalyst, pressures that are significantly higher than the equilibrium coexistence pressures are required to induce the graphite-to-diamond transition. Furthermore, the formation of the metastable hexagonal polymorph of diamond instead of the more stable cubic diamond is favored at lower temperatures. The concerted mechanism suggested in previous theoretical studies cannot explain these phenomena. Using an ab initio quality neural-network potential we performed a large-scale study of the graphite-to-diamond transition assuming that it occurs via nucleation. The nucleation mechanism accounts for the observed phenomenology and reveals its microscopic origins. We demonstrated that the large lattice distortions that accompany the formation of the diamond nuclei inhibit the phase transition at low pressure and direct it towards the hexagonal diamond phase at higher pressure. The nucleation mechanism proposed in this work is an important step towards a better understanding of structural transformations in a wide range of complex systems such as amorphous carbon and carbon nanomaterials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A self-learning algorithm for biased molecular dynamics

A new self-learning algorithm for accelerated dynamics, reconnaissance metadynamics, is proposed that is able to work with a very large number of collective coordinates. Acceleration of the dynamics is achieved by constructing a bias potential in terms of a patchwork of one-dimensional, locally valid collective coordinates. These collective coordinates are obtained from trajectory analyses so that they adapt to any new features encountered during the simulation. We show how this methodology can be used to enhance sampling in real chemical systems citing examples both from the physics of clusters and from the biological sciences.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio quality neural-network potential for sodium

An interatomic potential for high-pressure high-temperature (HPHT) crystalline and liquid phases of sodium is created using a neural-network (NN) representation of the ab initio potential energy surface. It is demonstrated that the NN potential provides an ab initio quality description of multiple properties of liquid sodium and bcc, fcc, cI16 crystal phases in the P-T region up to 120 GPa and 1200 K. The unique combination of computational efficiency of the NN potential and its ability to reproduce quantitatively experimental properties of sodium in the wide P-T range enables molecular dynamics simulations of physicochemical processes in HPHT sodium of unprecedented quality.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio quality study of the graphite-diamond phase coexistence

An interatomic potential for the diamond and graphite phases of carbon has been created using a neural-network (NN) representation of the ab initio potential energy surface. The NN potential combines the accuracy of a first-principle description of both phases with the efficiency of empirical force fields and allows one to perform, for the first time, a molecular dynamics study, of ab initio quality, of the thermodynamics of graphite-diamond coexistence. Good agreement between the experimental and calculated coexistence curves is achieved if nuclear quantum effects are included in the simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Displaced path integral formulation for the momentum distribution of quantum particles

The proton momentum distribution, accessible by deep inelastic neutron scattering, is a very sensitive probe of the potential of mean force experienced by the protons in hydrogen-bonded systems. In this work we introduce a novel estimator for the end to end distribution of the Feynman paths, i.e. the Fourier transform of the momentum distribution. In this formulation, free particle and environmental contributions factorize. Moreover, the environmental contribution has a natural analogy to a free energy surface in statistical mechanics, facilitating the interpretation of experiments. The new formulation is not only conceptually but also computationally advantageous. We illustrate the method with applications to an empirical water model, ab-initio ice, and one dimensional model systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient stochastic thermostatting of path integral molecular dynamics

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high-frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently-developed colored-noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nosé-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Langevin equation with colored noise for constant-temperature molecular dynamics simulations

We discuss the use of a Langevin equation with a colored (correlated) noise to perform constant-temperature molecular dynamics simulations. Since the equations of motion are linear in nature, it is easy to predict the response of a Hamiltonian system to such a thermostat and to tune at will the relaxation time of modes of different frequency. This allows one to optimize the time needed to thermalize the system and generate independent configurations. We show how this frequency-dependent response can be exploited to control the temperature of Car-Parrinello-like dynamics, keeping at low temperature the electronic degrees of freedom, without affecting the adiabatic separation from the vibrations of the ions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nuclear quantum effects in ab initio dynamics: theory and experiments for lithium imide

Owing to their small mass, hydrogen atoms exhibit strong quantum behavior even at room temperature. Including these effects in first principles calculations is challenging, because of the huge computational effort required by conventional techniques. Here we present the first ab-initio application of a recently-developed stochastic scheme, which allows to approximate nuclear quantum effects inexpensively. The proton momentum distribution of lithium imide, a material of interest for hydrogen storage, was experimentally measured by inelastic neutron scattering experiments and compared with the outcome of quantum thermostatted ab initio dynamics. We obtain favorable agreement between theory and experiments for this purely quantum mechanical property, thereby demonstrating that it is possible to improve the modelling of complex hydrogen-containing materials without additional computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Static disorder and structural correlations in the low temperature phase of lithium imide

Based on ab-initio molecular dynamics simulations, we investigate the low temperature crystal structure of Li2NH which in spite of its great interest as H-storage material is still matter of debate. The dynamical simulations reveal a precise correlation in the fractional occupation of Li sites which leads average atomic positions in excellent agreement with diffraction data and solves inconsistencies of previous proposals.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A hybrid approach to Fermi operator expansion

In a recent paper we have suggested that the finite temperature density matrix can be computed efficiently by a combination of polynomial expansion and iterative inversion techniques. We present here significant improvements over this scheme. The original complex-valued formalism is turned into a purely real one. In addition, we use Chebyshev polynomials expansion and fast summation techniques. This drastically reduces the scaling of the algorithm with the width of the Hamiltonian spectrum, which is now of the order of the cubic root of such parameter. This makes our method very competitive for applications to ab-initio simulations, when high energy resolution is required.

0 citations0 reviewsNo code linkedOpen access