Author: Pastor Richard Brooks Bernard Szabo Attila
Publisher: Taylor & Francis Ltd
ISSN: 1362-3028
Source: Molecular Physics, Vol.65, Iss.6, 1988-12, pp. : 1409-1419
Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.
Abstract
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 algorithms, and it is shown for harmonic oscillators that different definitions of the velocity are required to optimally calculate the temperature, pressure, and total energy, respectively.
Related content
Langevin and molecular dynamics simulations of Lennard-Jones liquids
By Giró A. Guardia E. Padró J.A.
Molecular Physics, Vol. 55, Iss. 5, 1985-08 ,pp. :
By Mehlig B. Heermann D.W. Forrest B.M.
Molecular Physics, Vol. 76, Iss. 6, 1992-08 ,pp. :
Algorithms for canonical molecular dynamics simulations
Molecular Physics, Vol. 72, Iss. 1, 1991-01 ,pp. :
An impulse integrator for Langevin dynamics
By Skeel Robert D. Izaguirre Jesús A.
Molecular Physics, Vol. 100, Iss. 24, 2002-12 ,pp. :