Collective motion artifacts arising in long-duration molecular dynamics simulations

Abstract

We tested a variety of molecular dynamics simulation strategies in long‐duration (up to several nanoseconds) constant‐temperature simulations of liquid water under periodic boundary conditions. Such long durations are necessary to achieve adequate conformational sampling in simulations of membrane assemblies and other large biomolecular systems. Under a variety of circumstances, serious artifacts arise in the form of spurious collective behavior that becomes obvious only after the simulation has gone at least several hundred picoseconds. The potential energy of the system drops and the system changes from a liquid to an icy or glassy state. The underlying cause is accumulated center‐of‐mass motion of the system, coupled with velocity rescaling associated with constant‐temperature control. The velocity rescaling in the constant‐temperature algorithm reduces the thermal velocity as the net center‐of‐mass velocity grows, effectively causing the kinetic energy of the system to drain from thermal motions into coordinated motions. We found that the incidence and magnitude of the underlying artifactual motion leading to the spurious transition is mediated by: choice of method for computing electrostatic interactions; choice of ensemble; size of the simulation cell; SHAKE tolerance; frequency of nonbonded pairlist updating; and closeness of coupling to the temperature bath. The appearance of the spurious transition can be avoided by periodically subtracting net center‐of‐mass motion during the dynamics, or by improving the accuracy of the simulation by means of tightening SHAKE tolerance and updating nonbonded pairlists every timestep. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 121–131, 2000

Thermostat Algorithms for Molecular Dynamics Simulations

“…For simulations started with vanishing overall momenta, one generally observes a very slow initial rise (often taking nanoseconds) followed by a very sudden burst of translational and rotational kinetic energy. The accumulation of kinetic energy in these degrees of freedom will effectively cool down the internal ones, giving rise to the so-called "flying ice cube effect" [138,139]. The most obvious remedy to this problem is to remove the overall center of mass motion from the atomic velocities at regular interval during the simulation.…”

Molecular dynamics simulation of nucleic acids: Successes, limitations, and promise

“…Other issues that the community should be aware of include the potential for inhibited sampling in minimal solvent models, the need for proper energy conservation and removal of net translational motion in periodic systems, 125,126 avoiding periodicity-induced artifacts in Ewald simulations, and biases of the available nucleic acid force fields. Further theoretical study of nucleic acids is necessary to better understand and overcome the limitations.…”

Establishing conditions for simulating hydrophobic solutes in electric fields by molecular dynamics

“…One can introduce this to molecular dynamics by the application of appropriate thermostats, which modify the Newtonian MD scheme such that a statistical ensemble is generated at a constant temperature. However, it is well documented that inappropriate choices of thermostats may significantly affect the fluctuations in the system and, in some cases, induce energy drifts caused by the accumulation of numerical errors [39,40].…”