A patch that imparts unconditional stability to explicit integrators for Langevin-like equations

Nawaf Bou-Rabee, Eric Vanden Eijnden

Research output: Contribution to journalArticle

Abstract

This paper extends the results in [8] to stochastic differential equations (SDEs) arising in molecular dynamics. It implements a patch to explicit integrators that consists of a Metropolis-Hastings step. The 'patched integrator' preserves the SDE's equilibrium distribution and is accurate on finite time intervals. As a corollary this paper proves the integrator's accuracy in estimating finite-time dynamics along an infinitely long solution - a first in molecular dynamics. The paper also covers multiple time-steps, holonomic constraints and scalability. Finally, the paper provides numerical tests supporting the theory.

Original languageEnglish (US)
Pages (from-to)2565-2580
Number of pages16
JournalJournal of Computational Physics
Volume231
Issue number6
DOIs
StatePublished - Mar 20 2012

Fingerprint

integrators
Molecular dynamics
molecular dynamics
Scalability
Differential equations
estimating
differential equations
intervals

Keywords

  • Metropolis-Hastings
  • Molecular dynamics
  • RATTLE
  • RESPA
  • Verlet

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

A patch that imparts unconditional stability to explicit integrators for Langevin-like equations. / Bou-Rabee, Nawaf; Vanden Eijnden, Eric.

In: Journal of Computational Physics, Vol. 231, No. 6, 20.03.2012, p. 2565-2580.

Research output: Contribution to journalArticle

@article{9f68698d82fd4fd2839a4d2bb05f843e,
title = "A patch that imparts unconditional stability to explicit integrators for Langevin-like equations",
abstract = "This paper extends the results in [8] to stochastic differential equations (SDEs) arising in molecular dynamics. It implements a patch to explicit integrators that consists of a Metropolis-Hastings step. The 'patched integrator' preserves the SDE's equilibrium distribution and is accurate on finite time intervals. As a corollary this paper proves the integrator's accuracy in estimating finite-time dynamics along an infinitely long solution - a first in molecular dynamics. The paper also covers multiple time-steps, holonomic constraints and scalability. Finally, the paper provides numerical tests supporting the theory.",
keywords = "Metropolis-Hastings, Molecular dynamics, RATTLE, RESPA, Verlet",
author = "Nawaf Bou-Rabee and {Vanden Eijnden}, Eric",
year = "2012",
month = "3",
day = "20",
doi = "10.1016/j.jcp.2011.12.007",
language = "English (US)",
volume = "231",
pages = "2565--2580",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "6",

}

TY - JOUR

T1 - A patch that imparts unconditional stability to explicit integrators for Langevin-like equations

AU - Bou-Rabee, Nawaf

AU - Vanden Eijnden, Eric

PY - 2012/3/20

Y1 - 2012/3/20

N2 - This paper extends the results in [8] to stochastic differential equations (SDEs) arising in molecular dynamics. It implements a patch to explicit integrators that consists of a Metropolis-Hastings step. The 'patched integrator' preserves the SDE's equilibrium distribution and is accurate on finite time intervals. As a corollary this paper proves the integrator's accuracy in estimating finite-time dynamics along an infinitely long solution - a first in molecular dynamics. The paper also covers multiple time-steps, holonomic constraints and scalability. Finally, the paper provides numerical tests supporting the theory.

AB - This paper extends the results in [8] to stochastic differential equations (SDEs) arising in molecular dynamics. It implements a patch to explicit integrators that consists of a Metropolis-Hastings step. The 'patched integrator' preserves the SDE's equilibrium distribution and is accurate on finite time intervals. As a corollary this paper proves the integrator's accuracy in estimating finite-time dynamics along an infinitely long solution - a first in molecular dynamics. The paper also covers multiple time-steps, holonomic constraints and scalability. Finally, the paper provides numerical tests supporting the theory.

KW - Metropolis-Hastings

KW - Molecular dynamics

KW - RATTLE

KW - RESPA

KW - Verlet

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

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

U2 - 10.1016/j.jcp.2011.12.007

DO - 10.1016/j.jcp.2011.12.007

M3 - Article

VL - 231

SP - 2565

EP - 2580

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 6

ER -