๐”– Bobbio Scriptorium
โœฆ   LIBER   โœฆ

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


Masking Resonance Artifacts in Force-Spl
โœ Adrian Sandu; Tamar Schlick ๐Ÿ“‚ Article ๐Ÿ“… 1999 ๐Ÿ› Elsevier Science ๐ŸŒ English โš– 420 KB

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

On velocity scaling in molecular dynamic
โœ Kersti Hermansson; George C. Lie; Enrico Clementi ๐Ÿ“‚ Article ๐Ÿ“… 1988 ๐Ÿ› John Wiley and Sons ๐ŸŒ English โš– 345 KB

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 for Actin
โœ Hitoshi Suda; Minoru Saito ๐Ÿ“‚ Article ๐Ÿ“… 1994 ๐Ÿ› Elsevier Science ๐ŸŒ English โš– 273 KB

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

Study of the electrostatics treatment in
โœ Robert Garemyr; Arne Elofsson ๐Ÿ“‚ Article ๐Ÿ“… 1999 ๐Ÿ› John Wiley and Sons ๐ŸŒ English โš– 355 KB ๐Ÿ‘ 3 views

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