Constant temperature simulations using the Langevin equation with velocity Verlet integration

Abstract 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.

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

[2]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[3]  S Doniach,et al.  Multiple time step diffusive Langevin dynamics for proteins , 1998, Proteins.

[4]  B. Brooks,et al.  Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamide , 1992, Biopolymers.

[5]  D Thirumalai,et al.  The nature of folded states of globular proteins , 1992, Biopolymers.

[6]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[7]  D. Ferguson,et al.  Isothermal-isobaric molecular dynamics simulations with Monte Carlo volume sampling , 1995 .

[8]  D. Grant,et al.  Generalized Langevin equations for molecular dynamics in solution , 1991 .

[9]  Sebastian Doniach,et al.  Long‐time overdamped Langevin dynamics of molecular chains , 1994, J. Comput. Chem..

[10]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

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

[12]  M. Karplus,et al.  Inertial effects in butane stochastic dynamics , 1989 .

[13]  W. Setzer,et al.  Conformations of saturated six-membered ring phosphorus heterocycles. X-ray crystallographic and H NMR study of cis-2-oxo-2-(dimethylamino)-3,5-diphenyl-1,3,2-oxozaphosphorinane, a cyclophosphamide-like molecule in a twist conformation , 1984 .

[14]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[15]  B. Brooks,et al.  An analysis of the accuracy of Langevin and molecular dynamics algorithms , 1988 .

[16]  DAVID M. FERGUSON,et al.  Parameterization and evaluation of a flexible water model , 1995, J. Comput. Chem..