Metadynamics: Difference between revisions
Jump to navigation
Jump to search
no edit summary
No edit summary |
No edit summary |
||
| Line 28: | Line 28: | ||
Since the biasing force is not known prior to the simulation, it is calculated during the simulation and adaptively applied. To accelerate sampling, the biasing force is applied even if an inadequate number of samples is collected in a bin. In this case, the biasing force is scaled in the early stages to avoid artifacts from applications of overestimated biasing forces. The biasing force can also be smoothed to decrease noise in collected data. For further details, see feimode in [[ABF:Controls]]. | Since the biasing force is not known prior to the simulation, it is calculated during the simulation and adaptively applied. To accelerate sampling, the biasing force is applied even if an inadequate number of samples is collected in a bin. In this case, the biasing force is scaled in the early stages to avoid artifacts from applications of overestimated biasing forces. The biasing force can also be smoothed to decrease noise in collected data. For further details, see feimode in [[ABF:Controls]]. | ||
===Sampling Space Discretization=== | |||
Due to numerical reasons, mean forces are collected on a regular grid. The averaging of instantaneous collective force is then done in small intervals centered at discrete CV values: | |||
<center><math>\frac{\partial G(\boldsymbol \xi^{*})}{\partial \boldsymbol \xi} = - {\langle \bold F_{ICF} \rangle}_{\boldsymbol \xi=\boldsymbol \xi^{*} \pm \bold h/2} </math> ... (5)</center> | |||
with the standard error: | |||
<center><math>\sigma_e \left( \frac{\partial G(\boldsymbol \xi^{*})}{\partial \boldsymbol \xi} \right) = \sqrt{ \frac {N_{corr}} {N_b} } \sigma(\bold F_{ICF})_ {\boldsymbol \xi=\boldsymbol \xi^{*} \pm \bold h/2} </math> ... (6)</center> | |||
where the standard deviation is given by: | |||
<center><math>\sigma(\bold F_{ICF})_ {\boldsymbol \xi=\boldsymbol \xi^{*} \pm \bold h/2} = \sqrt { \frac{1}{N_b} \sum_{i=1}^{N_b} \left( F_{ICF}(\boldsymbol \xi_i) - \overline {F_{ICF}} \right)^2 \vert _{\boldsymbol \xi=\boldsymbol \xi^{*} \pm \bold h/2} }</math> ... (7)</center> | |||
where <math>\bold h</math> is the interval size also called a bin, <math>N_b</math> is the number of samples collected in a bin centered at <math>\boldsymbol \xi^{*}</math>, and <math>N_{corr}</math> is a statistical inefficiency due to correlation in time series. | |||
Therefore, each CV involved in ABF simulations must be discretized by specifying an interval in which the sampling is performed and the number of intervals (bins) for discretization, for further details, see [[ABF:Collective variables]]. The increasing number of bins improves the accuracy of Equation 5 and subsequently the quality of the integrated free energy (Equation 1) but it also increases the noise because of a smaller number of samples collected in a bin. A reasonable compromise is the number of bins, which leads to 0.1 Å or 1-2° bin widths. | |||