Time integrators for molecular dynamics

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

This paper invites the reader to learn more about time integrators for Molecular Dynamics simulation through a simple MATLAB implementation. An overview of methods is provided from an algorithmic viewpoint that emphasizes long-time stability and finite-time dynamic accuracy. The given software simulates Langevin dynamics using an explicit, second-order (weakly) accurate integrator that exactly reproduces the Boltzmann-Gibbs density. This latter feature comes from adding a Metropolis acceptance-rejection step to the integrator. The paper discusses in detail the properties of the integrator. Since these properties do not rely on a specific form of a heat or pressure bath model, the given algorithm can be used to simulate other bath models including, e.g., the widely used v-rescale thermostat.

Original languageEnglish (US)
Pages (from-to)138-162
Number of pages25
JournalEntropy
Volume16
Issue number1
DOIs
StatePublished - Jan 2014

Fingerprint

integrators
molecular dynamics
baths
thermostats
readers
acceptability
rejection
computer programs
heat
simulation

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)

Keywords

  • Ergodicity
  • Explicit integrators
  • Metropolis algorithm
  • Weak accuracy

Cite this

Bou-Rabee, Nawaf. / Time integrators for molecular dynamics. In: Entropy. 2014 ; Vol. 16, No. 1. pp. 138-162.
@article{912cb7c6a4254f019d5a3c3b88a640b0,
title = "Time integrators for molecular dynamics",
abstract = "This paper invites the reader to learn more about time integrators for Molecular Dynamics simulation through a simple MATLAB implementation. An overview of methods is provided from an algorithmic viewpoint that emphasizes long-time stability and finite-time dynamic accuracy. The given software simulates Langevin dynamics using an explicit, second-order (weakly) accurate integrator that exactly reproduces the Boltzmann-Gibbs density. This latter feature comes from adding a Metropolis acceptance-rejection step to the integrator. The paper discusses in detail the properties of the integrator. Since these properties do not rely on a specific form of a heat or pressure bath model, the given algorithm can be used to simulate other bath models including, e.g., the widely used v-rescale thermostat.",
keywords = "Ergodicity, Explicit integrators, Metropolis algorithm, Weak accuracy",
author = "Nawaf Bou-Rabee",
year = "2014",
month = "1",
doi = "10.3390/e16010138",
language = "English (US)",
volume = "16",
pages = "138--162",
journal = "Entropy",
issn = "1099-4300",
publisher = "Multidisciplinary Digital Publishing Institute (MDPI)",
number = "1",

}

Time integrators for molecular dynamics. / Bou-Rabee, Nawaf.

In: Entropy, Vol. 16, No. 1, 01.2014, p. 138-162.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Time integrators for molecular dynamics

AU - Bou-Rabee, Nawaf

PY - 2014/1

Y1 - 2014/1

N2 - This paper invites the reader to learn more about time integrators for Molecular Dynamics simulation through a simple MATLAB implementation. An overview of methods is provided from an algorithmic viewpoint that emphasizes long-time stability and finite-time dynamic accuracy. The given software simulates Langevin dynamics using an explicit, second-order (weakly) accurate integrator that exactly reproduces the Boltzmann-Gibbs density. This latter feature comes from adding a Metropolis acceptance-rejection step to the integrator. The paper discusses in detail the properties of the integrator. Since these properties do not rely on a specific form of a heat or pressure bath model, the given algorithm can be used to simulate other bath models including, e.g., the widely used v-rescale thermostat.

AB - This paper invites the reader to learn more about time integrators for Molecular Dynamics simulation through a simple MATLAB implementation. An overview of methods is provided from an algorithmic viewpoint that emphasizes long-time stability and finite-time dynamic accuracy. The given software simulates Langevin dynamics using an explicit, second-order (weakly) accurate integrator that exactly reproduces the Boltzmann-Gibbs density. This latter feature comes from adding a Metropolis acceptance-rejection step to the integrator. The paper discusses in detail the properties of the integrator. Since these properties do not rely on a specific form of a heat or pressure bath model, the given algorithm can be used to simulate other bath models including, e.g., the widely used v-rescale thermostat.

KW - Ergodicity

KW - Explicit integrators

KW - Metropolis algorithm

KW - Weak accuracy

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

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

U2 - 10.3390/e16010138

DO - 10.3390/e16010138

M3 - Article

AN - SCOPUS:84893205677

VL - 16

SP - 138

EP - 162

JO - Entropy

JF - Entropy

SN - 1099-4300

IS - 1

ER -