A transition-rate investigation by molecular dynamics with the Langevin/implicit-Euler scheme

Anna M. Nyberg, Tamar Schlick

Research output: Contribution to journalArticle

Abstract

We report results from molecular dynamics simulations for a bistable piecewise-harmonic potential. A new method for molecular dynamics-the Langevin/implicit-Euler scheme-is investigated here and compared to the common Verlet integration algorithm. The implicit scheme introduces new computational and physical features since it (1) does not restrict integration time step to a very small value, and (2) effectively damps vibrational modes ω≫ωc, where ωc is a chosen cutoff frequency. The main issue we explore in this study is how different choices of time steps and cutoff frequencies affect computed transition rates. The one-dimensional, double-well model offers a simple visual and computational opportunity for observing the two different damping forces introduced by the scheme-frictional and intrinsic-and for characterizing the dominating force at a given parameter combination. Another question we examine here is the choice of time step below which the Langevin/implicit-Euler scheme produces "correct" transition rates for a model potential whose energy distribution is "well-described" classically.

Original languageEnglish (US)
Pages (from-to)4986-4996
Number of pages11
JournalThe Journal of chemical physics
Volume95
Issue number7
StatePublished - 1991

Fingerprint

Cutoff frequency
Molecular dynamics
molecular dynamics
cut-off
Potential energy
Damping
vibration mode
Computer simulation
energy distribution
damping
potential energy
harmonics
simulation

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

A transition-rate investigation by molecular dynamics with the Langevin/implicit-Euler scheme. / Nyberg, Anna M.; Schlick, Tamar.

In: The Journal of chemical physics, Vol. 95, No. 7, 1991, p. 4986-4996.

Research output: Contribution to journalArticle

@article{b55e1c060a0042b28984608e7bb27e27,
title = "A transition-rate investigation by molecular dynamics with the Langevin/implicit-Euler scheme",
abstract = "We report results from molecular dynamics simulations for a bistable piecewise-harmonic potential. A new method for molecular dynamics-the Langevin/implicit-Euler scheme-is investigated here and compared to the common Verlet integration algorithm. The implicit scheme introduces new computational and physical features since it (1) does not restrict integration time step to a very small value, and (2) effectively damps vibrational modes ω≫ωc, where ωc is a chosen cutoff frequency. The main issue we explore in this study is how different choices of time steps and cutoff frequencies affect computed transition rates. The one-dimensional, double-well model offers a simple visual and computational opportunity for observing the two different damping forces introduced by the scheme-frictional and intrinsic-and for characterizing the dominating force at a given parameter combination. Another question we examine here is the choice of time step below which the Langevin/implicit-Euler scheme produces {"}correct{"} transition rates for a model potential whose energy distribution is {"}well-described{"} classically.",
author = "Nyberg, {Anna M.} and Tamar Schlick",
year = "1991",
language = "English (US)",
volume = "95",
pages = "4986--4996",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "7",

}

TY - JOUR

T1 - A transition-rate investigation by molecular dynamics with the Langevin/implicit-Euler scheme

AU - Nyberg, Anna M.

AU - Schlick, Tamar

PY - 1991

Y1 - 1991

N2 - We report results from molecular dynamics simulations for a bistable piecewise-harmonic potential. A new method for molecular dynamics-the Langevin/implicit-Euler scheme-is investigated here and compared to the common Verlet integration algorithm. The implicit scheme introduces new computational and physical features since it (1) does not restrict integration time step to a very small value, and (2) effectively damps vibrational modes ω≫ωc, where ωc is a chosen cutoff frequency. The main issue we explore in this study is how different choices of time steps and cutoff frequencies affect computed transition rates. The one-dimensional, double-well model offers a simple visual and computational opportunity for observing the two different damping forces introduced by the scheme-frictional and intrinsic-and for characterizing the dominating force at a given parameter combination. Another question we examine here is the choice of time step below which the Langevin/implicit-Euler scheme produces "correct" transition rates for a model potential whose energy distribution is "well-described" classically.

AB - We report results from molecular dynamics simulations for a bistable piecewise-harmonic potential. A new method for molecular dynamics-the Langevin/implicit-Euler scheme-is investigated here and compared to the common Verlet integration algorithm. The implicit scheme introduces new computational and physical features since it (1) does not restrict integration time step to a very small value, and (2) effectively damps vibrational modes ω≫ωc, where ωc is a chosen cutoff frequency. The main issue we explore in this study is how different choices of time steps and cutoff frequencies affect computed transition rates. The one-dimensional, double-well model offers a simple visual and computational opportunity for observing the two different damping forces introduced by the scheme-frictional and intrinsic-and for characterizing the dominating force at a given parameter combination. Another question we examine here is the choice of time step below which the Langevin/implicit-Euler scheme produces "correct" transition rates for a model potential whose energy distribution is "well-described" classically.

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

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

M3 - Article

VL - 95

SP - 4986

EP - 4996

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 7

ER -