Potential of mean force

General description

The PMF is a type of simulation which can be done in Monte Carlo or Molecular Dynamics simulations which examines how a system's energy changes as a function of some specific reaction coordinate parameter. For example, it may examine how the system's energy changes as a function of the distance between two residues, or as a protein is pulled through a lipid bilayer. It can be a geometrical coordinate or a more general energetic (solvent) coordinate. Often PMF simulations are used in conjunction with umbrella sampling because typically the PMF simulation will fail to adequately sample the system space as it proceeds.[1]

Mathematical description

The Potential of Mean Force [2] of a system with N particles is strictly the potential that gives the average force over all the configurations of all the n+1...N particles acting on a particle j at any fixed configuration keeping fixed a set of particles 1...n

-\nabla _{j}w^{(n)}\,=\,{\frac {\int e^{-\beta V}(-\nabla _{j}V)dq_{n+1}...dq_{N}}{\int e^{-\beta V}dq_{n+1}....dq_{N}}},~j=1,2,....,n

Above, -\nabla _{j}w^{(n)} is the averaged force, i.e. "mean force" on particle j. And w^{(n)} is the so-called potential of mean force. For n=2, w^{(2)}(r) is the average work needed to bring the two particles from infinite separation to a distance r. It is also related to the radial distribution function of the system, g(r), by:[3]

g(r)=e^{-\beta w^{(2)}(r)}

Application

The potential of mean force w^{(2)} is usually applied in the Boltzmann inversion method as a first guess for the effective pair interaction potential that ought to reproduce the correct radial distribution function in a mesoscopic simulation.[4]

See also

References

  1. A. R. Leach, Molecular Modelling: Principles and Applications, 2001, ISBN 0-582-38210-6
  2. Kirkwood, J. G. Statistical Mechanics of fluid Mixtures. J. Chem. Phys. 1935, 3, 300; Statistical Mechanics of Liquid Solutions. Chemical Reviews 1936, 19, 275.
  3. See Chandler, section 7.3
  4. Reith, Dirk, Mathias Pütz, and Florian Müller‐Plathe. "Deriving effective mesoscale potentials from atomistic simulations." Journal of computational chemistry 24.13 (2003): 1624-1636.

Further reading

External links


This article is issued from Wikipedia - version of the Tuesday, December 15, 2015. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.