The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition

This article describes an unexpected phenomenon encountered during MD simulations: velocity rescaling using standard protocols can systematically change the proportion of total kinetic energy (KE) found in motions associated with the various degrees of freedom. Under these conditions, the simulation violates the principle of equipartition of energy, which requires a mean kinetic energy of RT/2 in each degree of freedom. A particularly pathological form of this problem occurs if one does not periodically remove the net translation of (and rotation about) the center of mass. In this case, almost all of the kinetic energy is converted into these two kinds of motion, producing a system with almost no kinetic energy associated with the internal degrees of freedom. We call this phenomenon “the flying ice cube.” We present a mathematical analysis of a simple diatomic system with two degrees of freedom, to document the origin of the problem. We then present examples from three kinds of MD simulations, one being an in vacuo simulation on a diatomic system, one involving a low resolution model of DNA in vacuo, and the third using a traditional all‐atom DNA model with full solvation, periodic boundary conditions, and the particle mesh Ewald method for treating long‐range electrostatics. Finally, we discuss methods for avoiding the problem. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 726–740, 1998

[1]  Quench echoes in molecular dynamics--A new phonon spectroscopy , 1980 .

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

[3]  Bernard R. Brooks,et al.  New spherical‐cutoff methods for long‐range forces in macromolecular simulation , 1994, J. Comput. Chem..

[4]  S C Harvey,et al.  Action at a distance in supercoiled DNA: effects of sequence on slither, branching, and intramolecular concentration. , 1996, Biophysical journal.

[5]  G. Grest,et al.  Normal-Mode Analysis by Quench-Echo Techniques: Localization in an Amorphous Solid , 1981 .

[6]  T Schlick,et al.  Biomolecular dynamics at long timesteps: bridging the timescale gap between simulation and experimentation. , 1997, Annual review of biophysics and biomolecular structure.

[7]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[8]  S. Harvey,et al.  Molecular modeling of closed circular DNA thermodynamic ensembles. , 1998, Biopolymers.

[9]  K. Schulten,et al.  Difficulties with multiple time stepping and fast multipole algorithm in molecular dynamics , 1997 .

[10]  Klaus Schulten,et al.  Temperature quench echoes in proteins , 1995 .

[11]  Peter A. Kollman,et al.  Molecular dynamics simulations highlight the structural differences among DNA: DNA, RNA:RNA, and DNA:RNA hybrid duplexes , 1997 .

[12]  K. Schulten,et al.  Velocity reassignment echoes in proteins , 1995 .

[13]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[14]  Stephen C. Harvey,et al.  Yammp: Development of a molecular mechanics program using the modular programming method , 1993, J. Comput. Chem..

[15]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[16]  S C Harvey,et al.  Flexible DNA: Genetically Unstable CTG·CAG and CGG·CCG from Human Hereditary Neuromuscular Disease Genes* , 1997, The Journal of Biological Chemistry.

[17]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[18]  Robert D. Skeel,et al.  Dangers of multiple time step methods , 1993 .

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

[20]  P A Kollman,et al.  Observation of the A-DNA to B-DNA transition during unrestrained molecular dynamics in aqueous solution. , 1996, Journal of molecular biology.

[21]  S. Harvey,et al.  Molecular dynamics simulations of small DNA plasmids: effects of sequence and supercoiling on intramolecular motions. , 1998, Biopolymers.

[22]  S. Harvey,et al.  Molecular mechanics model of supercoiled DNA. , 1989, Journal of molecular biology.

[23]  Karplus,et al.  Temperature echoes in molecular dynamics simulations of proteins. , 1993, Physical review letters.