Molecular dynamics algorithm for condensed systems with multiple time scales

Mark Tuckerman, Glenn J. Martyna, Bruce J. Berne

Research output: Contribution to journalArticle

Abstract

A frequently encountered problem in molecular dynamics simulations is the long runs required to study condensed systems consisting of both high frequency and low frequency degrees of freedom. Standard integrators require the choice of time step sufficiently small to guarantee stable solution of the highest frequency motion with the consequence that simulations require a very large number of central processing unit (CPU) cycles. In this note we present a new integrator that allows one to use a time step appropriate for the low frequency degrees of freedom without making any approximations related to the separation of time scales. This method is based on a choice of an analytically solvable reference system for the high frequency motion. We show how the analytical solution can be incorporated into a numerical integrator. The method is applied to two cases which are paradigms for this problem and it is shown that this approach and suitable generalizations should be very useful for future simulations of quantum and classical condensed matter systems.

Original languageEnglish (US)
Pages (from-to)1287-1291
Number of pages5
JournalThe Journal of chemical physics
Volume93
Issue number2
StatePublished - 1990

Fingerprint

integrators
Program processors
Molecular dynamics
molecular dynamics
Computer simulation
degrees of freedom
low frequencies
simulation
reference systems
central processing units
cycles
approximation

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Molecular dynamics algorithm for condensed systems with multiple time scales. / Tuckerman, Mark; Martyna, Glenn J.; Berne, Bruce J.

In: The Journal of chemical physics, Vol. 93, No. 2, 1990, p. 1287-1291.

Research output: Contribution to journalArticle

Tuckerman, Mark ; Martyna, Glenn J. ; Berne, Bruce J. / Molecular dynamics algorithm for condensed systems with multiple time scales. In: The Journal of chemical physics. 1990 ; Vol. 93, No. 2. pp. 1287-1291.
@article{b09789ae195449808982b3cc69067a23,
title = "Molecular dynamics algorithm for condensed systems with multiple time scales",
abstract = "A frequently encountered problem in molecular dynamics simulations is the long runs required to study condensed systems consisting of both high frequency and low frequency degrees of freedom. Standard integrators require the choice of time step sufficiently small to guarantee stable solution of the highest frequency motion with the consequence that simulations require a very large number of central processing unit (CPU) cycles. In this note we present a new integrator that allows one to use a time step appropriate for the low frequency degrees of freedom without making any approximations related to the separation of time scales. This method is based on a choice of an analytically solvable reference system for the high frequency motion. We show how the analytical solution can be incorporated into a numerical integrator. The method is applied to two cases which are paradigms for this problem and it is shown that this approach and suitable generalizations should be very useful for future simulations of quantum and classical condensed matter systems.",
author = "Mark Tuckerman and Martyna, {Glenn J.} and Berne, {Bruce J.}",
year = "1990",
language = "English (US)",
volume = "93",
pages = "1287--1291",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "2",

}

TY - JOUR

T1 - Molecular dynamics algorithm for condensed systems with multiple time scales

AU - Tuckerman, Mark

AU - Martyna, Glenn J.

AU - Berne, Bruce J.

PY - 1990

Y1 - 1990

N2 - A frequently encountered problem in molecular dynamics simulations is the long runs required to study condensed systems consisting of both high frequency and low frequency degrees of freedom. Standard integrators require the choice of time step sufficiently small to guarantee stable solution of the highest frequency motion with the consequence that simulations require a very large number of central processing unit (CPU) cycles. In this note we present a new integrator that allows one to use a time step appropriate for the low frequency degrees of freedom without making any approximations related to the separation of time scales. This method is based on a choice of an analytically solvable reference system for the high frequency motion. We show how the analytical solution can be incorporated into a numerical integrator. The method is applied to two cases which are paradigms for this problem and it is shown that this approach and suitable generalizations should be very useful for future simulations of quantum and classical condensed matter systems.

AB - A frequently encountered problem in molecular dynamics simulations is the long runs required to study condensed systems consisting of both high frequency and low frequency degrees of freedom. Standard integrators require the choice of time step sufficiently small to guarantee stable solution of the highest frequency motion with the consequence that simulations require a very large number of central processing unit (CPU) cycles. In this note we present a new integrator that allows one to use a time step appropriate for the low frequency degrees of freedom without making any approximations related to the separation of time scales. This method is based on a choice of an analytically solvable reference system for the high frequency motion. We show how the analytical solution can be incorporated into a numerical integrator. The method is applied to two cases which are paradigms for this problem and it is shown that this approach and suitable generalizations should be very useful for future simulations of quantum and classical condensed matter systems.

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

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

M3 - Article

AN - SCOPUS:0000173007

VL - 93

SP - 1287

EP - 1291

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 2

ER -