Multiple Time-Step Dual-Hamiltonian Hybrid Molecular Dynamics - Monte Carlo Canonical Propagation Algorithm

Yunjie Chen, Seyit Kale, Jonathan Weare, Aaron R. Dinner, Benoît Roux

Research output: Contribution to journalArticle

Abstract

A multiple time-step integrator based on a dual Hamiltonian and a hybrid method combining molecular dynamics (MD) and Monte Carlo (MC) is proposed to sample systems in the canonical ensemble. The Dual Hamiltonian Multiple Time-Step (DHMTS) algorithm is based on two similar Hamiltonians: a computationally expensive one that serves as a reference and a computationally inexpensive one to which the workload is shifted. The central assumption is that the difference between the two Hamiltonians is slowly varying. Earlier work has shown that such dual Hamiltonian multiple time-step schemes effectively precondition nonlinear differential equations for dynamics by reformulating them into a recursive root finding problem that can be solved by propagating a correction term through an internal loop, analogous to RESPA. Of special interest in the present context, a hybrid MD-MC version of the DHMTS algorithm is introduced to enforce detailed balance via a Metropolis acceptance criterion and ensure consistency with the Boltzmann distribution. The Metropolis criterion suppresses the discretization errors normally associated with the propagation according to the computationally inexpensive Hamiltonian, treating the discretization error as an external work. Illustrative tests are carried out to demonstrate the effectiveness of the method.

Original languageEnglish (US)
Pages (from-to)1449-1458
Number of pages10
JournalJournal of Chemical Theory and Computation
Volume12
Issue number4
DOIs
StatePublished - Apr 12 2016

Fingerprint

Hamiltonians
Molecular dynamics
molecular dynamics
propagation
Boltzmann distribution
integrators
acceptability
differential equations
Differential equations

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Multiple Time-Step Dual-Hamiltonian Hybrid Molecular Dynamics - Monte Carlo Canonical Propagation Algorithm. / Chen, Yunjie; Kale, Seyit; Weare, Jonathan; Dinner, Aaron R.; Roux, Benoît.

In: Journal of Chemical Theory and Computation, Vol. 12, No. 4, 12.04.2016, p. 1449-1458.

Research output: Contribution to journalArticle

Chen, Yunjie ; Kale, Seyit ; Weare, Jonathan ; Dinner, Aaron R. ; Roux, Benoît. / Multiple Time-Step Dual-Hamiltonian Hybrid Molecular Dynamics - Monte Carlo Canonical Propagation Algorithm. In: Journal of Chemical Theory and Computation. 2016 ; Vol. 12, No. 4. pp. 1449-1458.
@article{49f6daea1733464faf27fdcc2483acf2,
title = "Multiple Time-Step Dual-Hamiltonian Hybrid Molecular Dynamics - Monte Carlo Canonical Propagation Algorithm",
abstract = "A multiple time-step integrator based on a dual Hamiltonian and a hybrid method combining molecular dynamics (MD) and Monte Carlo (MC) is proposed to sample systems in the canonical ensemble. The Dual Hamiltonian Multiple Time-Step (DHMTS) algorithm is based on two similar Hamiltonians: a computationally expensive one that serves as a reference and a computationally inexpensive one to which the workload is shifted. The central assumption is that the difference between the two Hamiltonians is slowly varying. Earlier work has shown that such dual Hamiltonian multiple time-step schemes effectively precondition nonlinear differential equations for dynamics by reformulating them into a recursive root finding problem that can be solved by propagating a correction term through an internal loop, analogous to RESPA. Of special interest in the present context, a hybrid MD-MC version of the DHMTS algorithm is introduced to enforce detailed balance via a Metropolis acceptance criterion and ensure consistency with the Boltzmann distribution. The Metropolis criterion suppresses the discretization errors normally associated with the propagation according to the computationally inexpensive Hamiltonian, treating the discretization error as an external work. Illustrative tests are carried out to demonstrate the effectiveness of the method.",
author = "Yunjie Chen and Seyit Kale and Jonathan Weare and Dinner, {Aaron R.} and Beno{\^i}t Roux",
year = "2016",
month = "4",
day = "12",
doi = "10.1021/acs.jctc.5b00706",
language = "English (US)",
volume = "12",
pages = "1449--1458",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "4",

}

TY - JOUR

T1 - Multiple Time-Step Dual-Hamiltonian Hybrid Molecular Dynamics - Monte Carlo Canonical Propagation Algorithm

AU - Chen, Yunjie

AU - Kale, Seyit

AU - Weare, Jonathan

AU - Dinner, Aaron R.

AU - Roux, Benoît

PY - 2016/4/12

Y1 - 2016/4/12

N2 - A multiple time-step integrator based on a dual Hamiltonian and a hybrid method combining molecular dynamics (MD) and Monte Carlo (MC) is proposed to sample systems in the canonical ensemble. The Dual Hamiltonian Multiple Time-Step (DHMTS) algorithm is based on two similar Hamiltonians: a computationally expensive one that serves as a reference and a computationally inexpensive one to which the workload is shifted. The central assumption is that the difference between the two Hamiltonians is slowly varying. Earlier work has shown that such dual Hamiltonian multiple time-step schemes effectively precondition nonlinear differential equations for dynamics by reformulating them into a recursive root finding problem that can be solved by propagating a correction term through an internal loop, analogous to RESPA. Of special interest in the present context, a hybrid MD-MC version of the DHMTS algorithm is introduced to enforce detailed balance via a Metropolis acceptance criterion and ensure consistency with the Boltzmann distribution. The Metropolis criterion suppresses the discretization errors normally associated with the propagation according to the computationally inexpensive Hamiltonian, treating the discretization error as an external work. Illustrative tests are carried out to demonstrate the effectiveness of the method.

AB - A multiple time-step integrator based on a dual Hamiltonian and a hybrid method combining molecular dynamics (MD) and Monte Carlo (MC) is proposed to sample systems in the canonical ensemble. The Dual Hamiltonian Multiple Time-Step (DHMTS) algorithm is based on two similar Hamiltonians: a computationally expensive one that serves as a reference and a computationally inexpensive one to which the workload is shifted. The central assumption is that the difference between the two Hamiltonians is slowly varying. Earlier work has shown that such dual Hamiltonian multiple time-step schemes effectively precondition nonlinear differential equations for dynamics by reformulating them into a recursive root finding problem that can be solved by propagating a correction term through an internal loop, analogous to RESPA. Of special interest in the present context, a hybrid MD-MC version of the DHMTS algorithm is introduced to enforce detailed balance via a Metropolis acceptance criterion and ensure consistency with the Boltzmann distribution. The Metropolis criterion suppresses the discretization errors normally associated with the propagation according to the computationally inexpensive Hamiltonian, treating the discretization error as an external work. Illustrative tests are carried out to demonstrate the effectiveness of the method.

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

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

U2 - 10.1021/acs.jctc.5b00706

DO - 10.1021/acs.jctc.5b00706

M3 - Article

VL - 12

SP - 1449

EP - 1458

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 4

ER -