A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations
نویسندگان
چکیده
A method for sampling efficiently the free energy landscape of a complex system with respect to some given collective variables is proposed. Inspired by metadynamics [A. Laio, M. Parrinello, Proc. Nat. Acad. Sci. USA 99 (2002) 12562], we introduce an extended system where the collective variables are treated as dynamical ones and show that this allows to sample the free energy landscape of these variables directly. The sampling is accelerated by using an artificially high temperature for the collective variables. The validity of the method is established using general results for systems with multiple time-scales, and its numerical efficiency is also discussed via error analysis. We also show how the method can be modified in order to sample the reactive pathways in free energy space, and thereby analyze the mechanism of a reaction. Finally, we discuss how the method can be generalized and used as an alternative to the Kirkwood generalized thermodynamic integration approach for the calculation of free energy differences. 2006 Elsevier B.V. All rights reserved.
منابع مشابه
Temperature accelerated Monte Carlo (TAMC): a method for sampling the free energy surface of non-analytical collective variables.
We introduce a new method to simulate the physics of rare events. The method, an extension of the Temperature Accelerated Molecular Dynamics, comes in use when the collective variables introduced to characterize the rare events are either non-analytical or so complex that computing their derivative is not practical. We illustrate the functioning of the method by studying the homogeneous crystal...
متن کاملMetastability, conformation dynamics, and transition pathways in complex systems
We present a systematic introduction to the basic concepts and techniques for determining transition pathways and transition rates in systems with multiple metastable states. After discussing the classical transition state theory and its limitations, we derive a new set of equations for the optimal dividing surfaces. We then discuss transition path sampling, which is the most general technique ...
متن کاملAdaptive-boost molecular dynamics simulation of carbon diffusion in iron
We have developed an accelerated molecular dynamics (MD) method to model atomic-scale rare events. In this method, a smooth histogram of collective variables is first estimated by canonical ensemble molecular dynamics calculations, and then a temperature-dependent boost potential is iteratively constructed to accelerate the MD simulation. This method not only allows us to observe the rare event...
متن کاملGaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation
A Gaussian accelerated molecular dynamics (GaMD) approach for simultaneous enhanced sampling and free energy calculation of biomolecules is presented. By constructing a boost potential that follows Gaussian distribution, accurate reweighting of the GaMD simulations is achieved using cumulant expansion to the second order. Here, GaMD is demonstrated on three biomolecular model systems: alanine d...
متن کاملCalculating free energies using scaled-force molecular dynamics algorithm
One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or proteindrug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules betwee...
متن کامل