Masking Resonance Artifacts in Force-Splitting Methods for Biomolecular Simulations by Extrapolative Langevin Dynamics

Adrian Sandu, Tamar Schlick

Research output: Contribution to journalArticle

Abstract

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,Tmin), visible inaccuracies or instabilities can occur. Impulse-MTS schemes are vulnerable to these resonance errors since large energy pulses are introduced to the governing dynamics equations when the slow forces are evaluated. We recently showed that such resonance artifacts can be masked significantly by applying extrapolative splitting to stochastic dynamics. Theoretical and numerical analyses of force-splitting integrators based on the Verlet discretization are reported here for linear models to explain these observations and to suggest how to construct effective integrators for biomolecular dynamics that balance stability with accuracy. Analyses for Newtonian dynamics demonstrate the severe resonance patterns of the Impulse splitting, with this severity worsening with the outer timestep, Δt; Constant Extrapolation is generally unstable, but the disturbances do not grow with Δt. Thus, the stochastic extrapolative combination can counteract generic instabilities and largely alleviate resonances with a sufficiently strong Langevin heat-bath coupling (γ), estimates for which are derived here based on the fastest and slowest motion periods. These resonance results generally hold for nonlinear test systems: a water tetramer and solvated protein. Proposed related approaches such as Extrapolation/Correction and Midpoint Extrapolation work better than Constant Extrapolation only for timesteps less thanTmin/2. An effective extrapolative stochastic approach for biomolecules that balances long-timestep stability with good accuracy for the fast subsystem is then applied to a biomolecule using a three-class partitioning: the medium forces are treated byMidpoint Extrapolationvia position Verlet, and the slow forces are incorporated byConstant Extrapolation. The resulting algorithm (LN) performs well on a solvated protein system in terms of thermodynamic properties and yields an order of magnitude speedup with respect to single-timestep Langevin trajectories. Computed spectral density functions also show how the Newtonian modes can be approximated by using a small γ in the range of 5-20 ps-1.

Original languageEnglish (US)
Pages (from-to)74-113
Number of pages40
JournalJournal of Computational Physics
Volume151
Issue number1
DOIs
StatePublished - May 1 1999

Fingerprint

masking
artifacts
Extrapolation
extrapolation
simulation
integrators
Biomolecules
impulses
proteins
Proteins
Spectral density
Probability density function
baths
disturbances
Thermodynamic properties
thermodynamic properties
Trajectories
trajectories
heat
Computer simulation

Keywords

  • Force splitting
  • Molecular dynamics
  • Multiple timesteps
  • Numerical resonance
  • Symplecticness

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

Masking Resonance Artifacts in Force-Splitting Methods for Biomolecular Simulations by Extrapolative Langevin Dynamics. / Sandu, Adrian; Schlick, Tamar.

In: Journal of Computational Physics, Vol. 151, No. 1, 01.05.1999, p. 74-113.

Research output: Contribution to journalArticle

@article{c3e3c2daab3c4f96a9f798305266dae8,
title = "Masking Resonance Artifacts in Force-Splitting Methods for Biomolecular Simulations by Extrapolative Langevin Dynamics",
abstract = "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,Tmin), visible inaccuracies or instabilities can occur. Impulse-MTS schemes are vulnerable to these resonance errors since large energy pulses are introduced to the governing dynamics equations when the slow forces are evaluated. We recently showed that such resonance artifacts can be masked significantly by applying extrapolative splitting to stochastic dynamics. Theoretical and numerical analyses of force-splitting integrators based on the Verlet discretization are reported here for linear models to explain these observations and to suggest how to construct effective integrators for biomolecular dynamics that balance stability with accuracy. Analyses for Newtonian dynamics demonstrate the severe resonance patterns of the Impulse splitting, with this severity worsening with the outer timestep, Δt; Constant Extrapolation is generally unstable, but the disturbances do not grow with Δt. Thus, the stochastic extrapolative combination can counteract generic instabilities and largely alleviate resonances with a sufficiently strong Langevin heat-bath coupling (γ), estimates for which are derived here based on the fastest and slowest motion periods. These resonance results generally hold for nonlinear test systems: a water tetramer and solvated protein. Proposed related approaches such as Extrapolation/Correction and Midpoint Extrapolation work better than Constant Extrapolation only for timesteps less thanTmin/2. An effective extrapolative stochastic approach for biomolecules that balances long-timestep stability with good accuracy for the fast subsystem is then applied to a biomolecule using a three-class partitioning: the medium forces are treated byMidpoint Extrapolationvia position Verlet, and the slow forces are incorporated byConstant Extrapolation. The resulting algorithm (LN) performs well on a solvated protein system in terms of thermodynamic properties and yields an order of magnitude speedup with respect to single-timestep Langevin trajectories. Computed spectral density functions also show how the Newtonian modes can be approximated by using a small γ in the range of 5-20 ps-1.",
keywords = "Force splitting, Molecular dynamics, Multiple timesteps, Numerical resonance, Symplecticness",
author = "Adrian Sandu and Tamar Schlick",
year = "1999",
month = "5",
day = "1",
doi = "10.1006/jcph.1999.6202",
language = "English (US)",
volume = "151",
pages = "74--113",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "1",

}

TY - JOUR

T1 - Masking Resonance Artifacts in Force-Splitting Methods for Biomolecular Simulations by Extrapolative Langevin Dynamics

AU - Sandu, Adrian

AU - Schlick, Tamar

PY - 1999/5/1

Y1 - 1999/5/1

N2 - 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,Tmin), visible inaccuracies or instabilities can occur. Impulse-MTS schemes are vulnerable to these resonance errors since large energy pulses are introduced to the governing dynamics equations when the slow forces are evaluated. We recently showed that such resonance artifacts can be masked significantly by applying extrapolative splitting to stochastic dynamics. Theoretical and numerical analyses of force-splitting integrators based on the Verlet discretization are reported here for linear models to explain these observations and to suggest how to construct effective integrators for biomolecular dynamics that balance stability with accuracy. Analyses for Newtonian dynamics demonstrate the severe resonance patterns of the Impulse splitting, with this severity worsening with the outer timestep, Δt; Constant Extrapolation is generally unstable, but the disturbances do not grow with Δt. Thus, the stochastic extrapolative combination can counteract generic instabilities and largely alleviate resonances with a sufficiently strong Langevin heat-bath coupling (γ), estimates for which are derived here based on the fastest and slowest motion periods. These resonance results generally hold for nonlinear test systems: a water tetramer and solvated protein. Proposed related approaches such as Extrapolation/Correction and Midpoint Extrapolation work better than Constant Extrapolation only for timesteps less thanTmin/2. An effective extrapolative stochastic approach for biomolecules that balances long-timestep stability with good accuracy for the fast subsystem is then applied to a biomolecule using a three-class partitioning: the medium forces are treated byMidpoint Extrapolationvia position Verlet, and the slow forces are incorporated byConstant Extrapolation. The resulting algorithm (LN) performs well on a solvated protein system in terms of thermodynamic properties and yields an order of magnitude speedup with respect to single-timestep Langevin trajectories. Computed spectral density functions also show how the Newtonian modes can be approximated by using a small γ in the range of 5-20 ps-1.

AB - 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,Tmin), visible inaccuracies or instabilities can occur. Impulse-MTS schemes are vulnerable to these resonance errors since large energy pulses are introduced to the governing dynamics equations when the slow forces are evaluated. We recently showed that such resonance artifacts can be masked significantly by applying extrapolative splitting to stochastic dynamics. Theoretical and numerical analyses of force-splitting integrators based on the Verlet discretization are reported here for linear models to explain these observations and to suggest how to construct effective integrators for biomolecular dynamics that balance stability with accuracy. Analyses for Newtonian dynamics demonstrate the severe resonance patterns of the Impulse splitting, with this severity worsening with the outer timestep, Δt; Constant Extrapolation is generally unstable, but the disturbances do not grow with Δt. Thus, the stochastic extrapolative combination can counteract generic instabilities and largely alleviate resonances with a sufficiently strong Langevin heat-bath coupling (γ), estimates for which are derived here based on the fastest and slowest motion periods. These resonance results generally hold for nonlinear test systems: a water tetramer and solvated protein. Proposed related approaches such as Extrapolation/Correction and Midpoint Extrapolation work better than Constant Extrapolation only for timesteps less thanTmin/2. An effective extrapolative stochastic approach for biomolecules that balances long-timestep stability with good accuracy for the fast subsystem is then applied to a biomolecule using a three-class partitioning: the medium forces are treated byMidpoint Extrapolationvia position Verlet, and the slow forces are incorporated byConstant Extrapolation. The resulting algorithm (LN) performs well on a solvated protein system in terms of thermodynamic properties and yields an order of magnitude speedup with respect to single-timestep Langevin trajectories. Computed spectral density functions also show how the Newtonian modes can be approximated by using a small γ in the range of 5-20 ps-1.

KW - Force splitting

KW - Molecular dynamics

KW - Multiple timesteps

KW - Numerical resonance

KW - Symplecticness

UR - http://www.scopus.com/inward/record.url?scp=0002403441&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0002403441&partnerID=8YFLogxK

U2 - 10.1006/jcph.1999.6202

DO - 10.1006/jcph.1999.6202

M3 - Article

VL - 151

SP - 74

EP - 113

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 1

ER -