Langevin dynamics

Langevin dynamics is an approach to mechanics using simplified models and using stochastic differential equations to account for omitted degrees of freedom.A molecular system in the real world is unlikely to be present in vacuum. Jostling of solvent or air molecules causes friction, and the occasional high velocity collision will perturb the system. Langevin dynamics attempts to extend molecular dynamics to allow for these effects. Also, Langevin dynamics allow controlling the temperature like a thermostat, thus approximating the canonical ensemble.

Langevin dynamics mimic the viscous aspect of a solvent. In itself, it is not a complete implicit solvent, i.e. it does not account for either the electrostatic screening nor the hydrophobic effect. For a system of $N$ particles with masses $M$, coordinates $X$ and velocities $V$, the continuous form of the simplest Langevin equation is [cite book

first=Tamar | last=Schlick | year=2002 | title=Molecular Modeling and Simulation | publisher=Springer | id=ISBN 0-387-95404-X | pages = 435-438

]

$F(X)\; =\; -\; abla\; U(X)-\; gamma\; Mdot\{X\}+R(t)=Mdot\{V\}\; (t)$

Where $F(X)$ is the force felt by the particles, $U(X)$ is the potential energy (e.g., the force field), $gamma$ is the collision parameter or damping constant (reciprocal time units). The dots represent first derivatives with respect to time (e.g. $dot\{X\}$ is the velocity of the particles).

R(t) is a random force vector, which is a stationary Gaussian process with zero-mean:

:$leftlangle\; R(t)\; ight\; angle\; =0$:$leftlangle\; R(t)R(t\text{'})^T\; ight\; angle\; =2\; gamma\; k\_B\; T\; M\; delta(t-t\text{'})$

Where T is the target temperature, k_B is Boltzmann's constant, and $delta$ is the Dirac delta.

If the main objective is to control temperature, care should be exercised to use a small damping constant $gamma$. As $gamma$ grows, it spans the inertial all the way to the diffusive (Brownian) regime. The Langevin dynamics limit of non-inertia is commonly described as Brownian dynamics.

Langevin differential equations which govern a random variable X can bereformulated as "Fokker-Planck" differential equations (Fokker-Planck Eq'n), or "masterequations," which govern the probability distribution of X.

References

ee also

* Hamiltonian mechanics
* Molecular dynamics
* Brownian dynamics
* Statistical mechanics
* Implicit solvation
* stochastic differential equations

External links

* [http://cmm.cit.nih.gov/intro_simulation/node24.html Langevin Dynamics (LD) Simulation]