Increasing the time step in molecular dynamics

Anna M. Nyberg, Tamar Schlick

Research output: Contribution to journalArticle

Abstract

We examine how the two properties of the Langevin/implicit-Euler scheme for molecular dynamics (MD) - stability over large time steps and damping of high-frequency vibrations greater than some cutoff frequency ωc - affect the dynamics of liquid butane model. For time steps Δt≤ 20 fs with a "natural" cutoff ωc=kBT/⇌h (kB and ⇌h are Boltzmann's and Planck's constants, T = temperature), the associated kinetic energy is greater than 85% of its value of Δt=1 fs, where the connection to a heat bath is weak and the generated trajectories by our scheme and conventional MD are similar. At larger Δt the connection to the heat bath is much stronger, and for Δt > 20 fs the intrinsic numerical damping of the scheme sets in. The comparison of bond-length and dihedral-angle distributions at three different time steps reveals a small, broadening trend at larger Δt. The differences in a dynamic property, the velocity autocorrelation function, are however much larger. There is a drastic difference for 2 and 20 fs, and for Δt > 40 fs most of the motion between time steps is damped and more random. Thus, while the Langevin equation per se gives a Boltzmann distribution, the expected configurational sampling can be obtained as long as the numerical damping does not disturb the balance between the random and damping terms. Consequently, for studying certain dynamic functions time steps in the same range as in conventional MD are needed, but for static properties larger time steps can be used.

Original languageEnglish (US)
Pages (from-to)538-546
Number of pages9
JournalChemical Physics Letters
Volume198
Issue number6
DOIs
StatePublished - Oct 23 1992

Fingerprint

Molecular dynamics
Damping
molecular dynamics
damping
baths
Cutoff frequency
Bond length
Dihedral angle
cut-off
Theophylline
Autocorrelation
Kinetic energy
Vibrations (mechanical)
dynamic stability
heat
Boltzmann distribution
time functions
Trajectories
butanes
Sampling

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Spectroscopy
  • Condensed Matter Physics
  • Atomic and Molecular Physics, and Optics
  • Surfaces and Interfaces

Cite this

Increasing the time step in molecular dynamics. / Nyberg, Anna M.; Schlick, Tamar.

In: Chemical Physics Letters, Vol. 198, No. 6, 23.10.1992, p. 538-546.

Research output: Contribution to journalArticle

Nyberg, Anna M. ; Schlick, Tamar. / Increasing the time step in molecular dynamics. In: Chemical Physics Letters. 1992 ; Vol. 198, No. 6. pp. 538-546.
@article{db6d7c075deb474c8c32be5cef229bd7,
title = "Increasing the time step in molecular dynamics",
abstract = "We examine how the two properties of the Langevin/implicit-Euler scheme for molecular dynamics (MD) - stability over large time steps and damping of high-frequency vibrations greater than some cutoff frequency ωc - affect the dynamics of liquid butane model. For time steps Δt≤ 20 fs with a {"}natural{"} cutoff ωc=kBT/⇌h (kB and ⇌h are Boltzmann's and Planck's constants, T = temperature), the associated kinetic energy is greater than 85{\%} of its value of Δt=1 fs, where the connection to a heat bath is weak and the generated trajectories by our scheme and conventional MD are similar. At larger Δt the connection to the heat bath is much stronger, and for Δt > 20 fs the intrinsic numerical damping of the scheme sets in. The comparison of bond-length and dihedral-angle distributions at three different time steps reveals a small, broadening trend at larger Δt. The differences in a dynamic property, the velocity autocorrelation function, are however much larger. There is a drastic difference for 2 and 20 fs, and for Δt > 40 fs most of the motion between time steps is damped and more random. Thus, while the Langevin equation per se gives a Boltzmann distribution, the expected configurational sampling can be obtained as long as the numerical damping does not disturb the balance between the random and damping terms. Consequently, for studying certain dynamic functions time steps in the same range as in conventional MD are needed, but for static properties larger time steps can be used.",
author = "Nyberg, {Anna M.} and Tamar Schlick",
year = "1992",
month = "10",
day = "23",
doi = "10.1016/0009-2614(92)85028-9",
language = "English (US)",
volume = "198",
pages = "538--546",
journal = "Chemical Physics Letters",
issn = "0009-2614",
publisher = "Elsevier",
number = "6",

}

TY - JOUR

T1 - Increasing the time step in molecular dynamics

AU - Nyberg, Anna M.

AU - Schlick, Tamar

PY - 1992/10/23

Y1 - 1992/10/23

N2 - We examine how the two properties of the Langevin/implicit-Euler scheme for molecular dynamics (MD) - stability over large time steps and damping of high-frequency vibrations greater than some cutoff frequency ωc - affect the dynamics of liquid butane model. For time steps Δt≤ 20 fs with a "natural" cutoff ωc=kBT/⇌h (kB and ⇌h are Boltzmann's and Planck's constants, T = temperature), the associated kinetic energy is greater than 85% of its value of Δt=1 fs, where the connection to a heat bath is weak and the generated trajectories by our scheme and conventional MD are similar. At larger Δt the connection to the heat bath is much stronger, and for Δt > 20 fs the intrinsic numerical damping of the scheme sets in. The comparison of bond-length and dihedral-angle distributions at three different time steps reveals a small, broadening trend at larger Δt. The differences in a dynamic property, the velocity autocorrelation function, are however much larger. There is a drastic difference for 2 and 20 fs, and for Δt > 40 fs most of the motion between time steps is damped and more random. Thus, while the Langevin equation per se gives a Boltzmann distribution, the expected configurational sampling can be obtained as long as the numerical damping does not disturb the balance between the random and damping terms. Consequently, for studying certain dynamic functions time steps in the same range as in conventional MD are needed, but for static properties larger time steps can be used.

AB - We examine how the two properties of the Langevin/implicit-Euler scheme for molecular dynamics (MD) - stability over large time steps and damping of high-frequency vibrations greater than some cutoff frequency ωc - affect the dynamics of liquid butane model. For time steps Δt≤ 20 fs with a "natural" cutoff ωc=kBT/⇌h (kB and ⇌h are Boltzmann's and Planck's constants, T = temperature), the associated kinetic energy is greater than 85% of its value of Δt=1 fs, where the connection to a heat bath is weak and the generated trajectories by our scheme and conventional MD are similar. At larger Δt the connection to the heat bath is much stronger, and for Δt > 20 fs the intrinsic numerical damping of the scheme sets in. The comparison of bond-length and dihedral-angle distributions at three different time steps reveals a small, broadening trend at larger Δt. The differences in a dynamic property, the velocity autocorrelation function, are however much larger. There is a drastic difference for 2 and 20 fs, and for Δt > 40 fs most of the motion between time steps is damped and more random. Thus, while the Langevin equation per se gives a Boltzmann distribution, the expected configurational sampling can be obtained as long as the numerical damping does not disturb the balance between the random and damping terms. Consequently, for studying certain dynamic functions time steps in the same range as in conventional MD are needed, but for static properties larger time steps can be used.

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

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

U2 - 10.1016/0009-2614(92)85028-9

DO - 10.1016/0009-2614(92)85028-9

M3 - Article

AN - SCOPUS:0142080331

VL - 198

SP - 538

EP - 546

JO - Chemical Physics Letters

JF - Chemical Physics Letters

SN - 0009-2614

IS - 6

ER -