TY - JOUR

T1 - Constant temperature simulations using the Langevin equation with velocity Verlet integration

AU - Paterlini, M. Germana

AU - Ferguson, David M.

N1 - Copyright:
Copyright 2018 Elsevier B.V., All rights reserved.

PY - 1998/9/15

Y1 - 1998/9/15

N2 - An algorithm, which reduces to velocity Verlet in the limit of zero friction, is obtained for the generalized Langevin equation. The formulation presented is unique in that the velocities are based on a direct second order Taylor expansion of the inertial forces. The resulting finite difference equation is compared with a previous formulation of Verlet-based Langevin dynamics. The equations are implemented to study the physical properties of dense neon and liquid water at constant temperatures as a function of the friction rate γ. The results show canonical ensembles can be obtained at moderate friction rates without recurring to switching functions to control energy losses, or overdamping. The LD methodology is further applied to a protein in solution and compared to MD simulations which use a more conventional scaling function for temperature control. Analyses of the MD trajectories show the protein remains at either lower or higher temperature compared to the solvent, depending upon the treatment of the potential function at the system's boundary. On the contrary, LD simulations show a rapid and even equilibration between protein and solvent in all cases.

AB - An algorithm, which reduces to velocity Verlet in the limit of zero friction, is obtained for the generalized Langevin equation. The formulation presented is unique in that the velocities are based on a direct second order Taylor expansion of the inertial forces. The resulting finite difference equation is compared with a previous formulation of Verlet-based Langevin dynamics. The equations are implemented to study the physical properties of dense neon and liquid water at constant temperatures as a function of the friction rate γ. The results show canonical ensembles can be obtained at moderate friction rates without recurring to switching functions to control energy losses, or overdamping. The LD methodology is further applied to a protein in solution and compared to MD simulations which use a more conventional scaling function for temperature control. Analyses of the MD trajectories show the protein remains at either lower or higher temperature compared to the solvent, depending upon the treatment of the potential function at the system's boundary. On the contrary, LD simulations show a rapid and even equilibration between protein and solvent in all cases.

UR - http://www.scopus.com/inward/record.url?scp=0032363602&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0032363602&partnerID=8YFLogxK

U2 - 10.1016/S0301-0104(98)00214-6

DO - 10.1016/S0301-0104(98)00214-6

M3 - Article

AN - SCOPUS:0032363602

VL - 236

SP - 243

EP - 252

JO - Chemical Physics

JF - Chemical Physics

SN - 0301-0104

IS - 1-3

ER -