Numerical resonance artifacts have become recognized recently as a limiting factor to increasing the timestep in multiple-timestep (MTS) biomolecular dynamics simulations. At certain timesteps correlated to internal motions (e.g., 5 fs, around half the period of the fastest bond stretch, T min ), vi
Nonlinear Resonance Artifacts in Molecular Dynamics Simulations
โ Scribed by Tamar Schlick; Margaret Mandziuk; Robert D. Skeel; K. Srinivas
- Publisher
- Elsevier Science
- Year
- 1998
- Tongue
- English
- Weight
- 398 KB
- Volume
- 140
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
โฆ Synopsis
The intriguing phenomenon of resonance, a pronounced integrator-induced corruption of a system's dynamics, is examined for simple molecular systems subject to the classical equations of motion. This source of timestep limitation is not well appreciated in general, and certainly analyses of resonance patterns have been few in connection to biomolecular dynamics. Yet resonances are present in the commonly used Verlet integrator, in symplectic implicit schemes, and also limit the scope of current multiple-timestep methods that are formulated as symplectic and reversible. The only general remedy to date has been to reduce the timestep. For this purpose, we derive method-dependent timestep thresholds (e.g., Tables 1 and2) that serve as useful guidelines in practice for biomolecular simulations. We also devise closely related symplectic implicit schemes for which the limitation on the discretization stepsize is much less severe. Specifically, we design methods to remove third-order, or both the third-and fourth-order, resonances. These severe low-order resonances can lead to instability or very large energies. Our tests on two simple molecular problems (Morse and Lennard-Jones potentials), as well as a 22-atom molecule, N-acetylalanyl-Nmethylamide, confirm this prediction; our methods can delay resonances so that they occur only at larger timesteps (EW method) or are essentially removed (LIM2 method). Although stable for large timesteps by this approach, trajectories show large energy fluctuations, perhaps due to the coupling with other factors that induce instability in complex nonlinear systems. Thus, the methods developed here may be more useful for conformational sampling of biomolecular structures. The analysis presented here for the blocked alanine model emphasizes that one-dimensional analysis of resonances can be applied to a more complex, multimode system to analyze resonance behavior, but that resonance due to frequency coupling is more complex to pinpoint. More generally, instability, apparently due to numerically induced
๐ SIMILAR VOLUMES
The effect of scaling the molecular velocities to a fixed total energy in molecular dynamics simulations within the (N,V,E) ensemble has been investigated. The effect of using different time steps is also discussed. It is found that, even for small time steps, velocity scaling has a substantial infl
Molecular dynamics simulations were performed for an actin monomer without \(\mathrm{Ca}^{2+}\) and nucleotide for a period of \(100 \mathrm{psec}\) in order to study the stable structure and fluctuations of actin monomers in solution. For the initial structure for the molecular dynamics simulations
This article considers the treatment of long-range interactions in molecular dynamics simulations. We investigate the effects of using different cutoff distances, constant versus distancedependent dielectric, and different smoothing methods. In contrast to findings of earlier studies, we find that i