An analysis of the accuracy of Langevin and molecular dynamics algorithms

Analytic expressions for mean squared positions and velocities of a harmonic oscillator are derived for Langevin dynamics algorithms valid in the high and low friction limits, and for the Verlet algorithm. For typical values of the parameters, errors in the positions are small. However, if the velocity is defined by the usual Verlet form, kinetic energies (and therefore calculated temperatures) can be in error by several per cent for the Langevin algorithms. If the Bunger-Brooks-Karplus algorithm is used to calculate positions, a simple redefinition of the velocity results greatly in improved kinetic energies. In addition, due to cancellation of errors in the velocities and the positions, the correct virial is obtained. The effect of including the force derivative in diffusive algorithms is examined. Positional and velocity averages are calculated for the Verlet algorithm for arbitrary initial conditions, and errors in the total energy and virial are analysed. Connection is made with the Langevin algorith...

[1]  R. Pastor,et al.  Frictional models for stochastic simulations of proteins , 1988, Biopolymers.

[2]  M. Karplus,et al.  Parametrization of the friction constant for stochastic simulations of polymers , 1988 .

[3]  William G. Hoover,et al.  Simple exact for well known molecular dynamics algorithms , 1987 .

[4]  R. W. Hockney,et al.  Time integration in particle models , 1987 .

[5]  J A McCammon,et al.  Large‐amplitude bending motions in phenylalanine transfer RNA , 1984, Biopolymers.

[6]  M. Karplus,et al.  Stochastic boundary conditions for molecular dynamics simulations of ST2 water , 1984 .

[7]  Herman J. C. Berendsen,et al.  ALGORITHMS FOR BROWNIAN DYNAMICS , 1982 .

[8]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[9]  Michael P. Allen,et al.  Brownian dynamics simulation of a chemical reaction in solution , 1980 .

[10]  D. Ermak,et al.  Numerical integration of the Langevin equation: Monte Carlo simulation , 1980 .

[11]  S. Orszag,et al.  Advanced Mathematical Methods For Scientists And Engineers , 1979 .

[12]  T. Schneider,et al.  Molecular-dynamics study of a three-dimensional one-component model for distortive phase transitions , 1978 .

[13]  David Beeman,et al.  Some Multistep Methods for Use in Molecular Dynamics Calculations , 1976 .

[14]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[15]  S. Chandrasekhar.  Stochastic problems in Physics and Astronomy , 1943 .