Molecular Simulation/Umbrella Sampling
Umbrella sampling is a sampling method used in computational physics and chemistry. This sampling can sample the rare states which normal molecular dynamic sampling ignored. Therefore, umbrella sampling can improve free energy calculation when a system is undergoing a systematic change.
Biased molecular dynamics simulations
[edit | edit source]Normal MD simulations samples system in equilibrium. In an MD simulation of the time series of the C-C-C-C dihedral angle of n-butane(aq), only gauche states and trans states are sampled. Because this simulation is only performed in 2 ns, states with high free energy (e.g. cis state) are less likely to happen. These configurations are ignored and it is impossible to calculate the free energy of these states from this simulation. An artificial bias potential is needed in this case to help the molecule cross the energy barrier. With bias potential, rare states can be effectively sampled.
In this case, a harmonic bias potential is needed to counteract the dihedral barrier.
High free energy states were captured by biased simulation. In order to calculate the free energy profile of these states, biased probability distribution has to be converted to an unbiased probability distribution.
Acquire free energy profile from biased simulations
[edit | edit source]The potential energy includes the bias potential at the reaction coordinate is
The probability distribution of this potential is
The probability distribution of unbiased potential is
From this equation, we can derive,
Free energy profile can be calculated from probability distribution by,
Using this relation, the PMF of the biased simulation can be converted to unbiased PMF by:
term is denoted as . It is generally a constant and in some cases does not affect the relative energy and no needed to calculate. It can be calculated by [1]:
This method provides free energy profile of all possible states. In umbrella sampling of n-butane(aq), the chosen bias potential covered all reaction coordinates. General cases are more complex, which leads to a more complex determination of bias potential.
Choice of Bias Potential
[edit | edit source]The previous section discussed the biased molecular dynamic simulation of n-butane(aq). The reaction coordinate is one-dimensional and periodic, and the bias potential was chosen to be the negative the dihedral potential of n-butane[2]. The optimum bias potential is the opposite of the free energy [1]. However, is unknown for most cases. For general cases, the bias potential needs to be adjusted along the reaction coordinate. Thus, a harmonic bias potential restrained on a reference point with respect to a window on the reaction coordinate is introduced[2]:
Therefore, a full umbrella sampling can be obtained by a series of biased MD simulation on different reference points on the reaction coordinate.
Calculation of the Potential of Mean Force from Umbrella Sampling Data
[edit | edit source]The Weighted Histogram Analysis Method (WHAM)[3] transferred a series of biased distribution functions to a single distribution function. The value of needs to be estimated to give the correct value of :
The true distribution P(s) is the weight average of each step[1]:
And , where is the total number of steps sampled for window [3].
Combined with , both and can be obtained.
The other way to analyze umbrella sampling is Umbrella Integration, see[1].
See also
[edit | edit source]For more information about umbrella sampling, see[4]
References
[edit | edit source]- ↑ a b c d Kästner, Johannes (2011). "Umbrella sampling". Wiley Interdisciplinary Reviews: Computational Molecular Science. 1 (6): 932–942.
- ↑ a b Invalid
<ref>
tag; no text was provided for refs namedr1
- ↑ a b Kumar, S; Rosenberg, JM; Bouzida, D; Swendsen, RH (1992). "The weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method". Journal of computational chemistry. 13 (8): 1011–1021.
{{cite journal}}
: Cite has empty unknown parameter:|1=
(help) - ↑ Torrie, GM; Valleau, JP (February 1977). "Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling". Journal of Computational Physics. 23 (2): 187–199. doi:10.1016/0021-9991(77)90121-8.