Metadynamics: Difference between revisions

From VASP Wiki
(Redirected page to Category:Metadynamics)
Tag: New redirect
(Removed redirect to Category:Metadynamics)
Tag: Removed redirect
 
Line 1: Line 1:
#REDIRECT [[Category:Metadynamics]]
In metadynamics,{{cite|laio:pnas:02}}{{cite|iannuzzi:prl:03}} the bias potential
that acts on a selected number of geometric parameters (collective variables) &xi;={&xi;<sub>1</sub>, &xi;<sub>2</sub>, ...,&xi;<sub>''m''</sub>} is constructed on-the-fly during the simulation. The Hamiltonian for the metadynamics <math>\tilde{H}(q,p)</math> can be written as:
:<math>
\tilde{H}(q,p,t) = H(q,p) + \tilde{V}(t,\xi),
</math>
where <math>H(q,p)</math> is the Hamiltonian for the original (unbiased) system, and <math>\tilde{V}(t,\xi)</math> is the time-dependent bias potential. The latter term is usually defined as a sum of Gaussian hills with height ''h'' and width ''w'':
:<math>
\tilde{V}(t,\xi) = h \sum_{i=1}^{\lfloor t/t_G \rfloor} \exp{\left\{ -\frac{|\xi^{(t)}-\xi^{(i \cdot t_G)}|^2}{2
w^2} \right\}}.
</math>
In practice, <math>\tilde{V}(t,\xi)</math> is updated by adding a new Gaussian with a time increment ''t''<sub>G</sub>, which is typically one or two orders of magnitude greater than the time step used in the MD simulation.
 
In the limit of infinite simulation time, the bias potential is related to the free energy via:
:<math>
A(\xi) = - \lim_{t \to \infty} \tilde{V}(t,\xi) + const.
</math>
Practical hints as how to adjust the parameters used in metadynamics (''h'', ''w'', ''t''<sub>G</sub>) are given in Refs.{{cite|ensing:jpc:05}} and {{cite|laio:jpc:05}}.
 
The error estimation in free-energy calculations with metadynamics is discussed in Ref.{{cite|laio:jpc:05}}.
 
== References ==
<references/>
 
<noinclude>
----
 
[[Category:Advanced molecular-dynamics sampling]][[Category:Theory]]

Latest revision as of 11:53, 16 October 2024

In metadynamics,[1][2] the bias potential that acts on a selected number of geometric parameters (collective variables) ξ={ξ1, ξ2, ...,ξm} is constructed on-the-fly during the simulation. The Hamiltonian for the metadynamics can be written as:

where is the Hamiltonian for the original (unbiased) system, and is the time-dependent bias potential. The latter term is usually defined as a sum of Gaussian hills with height h and width w:

In practice, is updated by adding a new Gaussian with a time increment tG, which is typically one or two orders of magnitude greater than the time step used in the MD simulation.

In the limit of infinite simulation time, the bias potential is related to the free energy via:

Practical hints as how to adjust the parameters used in metadynamics (h, w, tG) are given in Refs.[3] and [4].

The error estimation in free-energy calculations with metadynamics is discussed in Ref.[4].

References