Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration

Alex Albaugh, Mark Tuckerman, Teresa Head-Gordon

Research output: Contribution to journalArticle

Abstract

In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ∼6-8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. The combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.

Original languageEnglish (US)
Pages (from-to)2195-2205
Number of pages11
JournalJournal of Chemical Theory and Computation
Volume15
Issue number4
DOIs
StatePublished - Apr 9 2019

Fingerprint

iteration
self consistent fields
Molecular dynamics
Polarization
Potential energy surfaces
polarization
Peptides
Salts
molecular dynamics
Proteins
Water
Computer simulation
field theory (physics)
peptides
simulation
potential energy
aqueous solutions
salts
proteins
water

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration. / Albaugh, Alex; Tuckerman, Mark; Head-Gordon, Teresa.

In: Journal of Chemical Theory and Computation, Vol. 15, No. 4, 09.04.2019, p. 2195-2205.

Research output: Contribution to journalArticle

@article{1691cbc0c51e47f49bdde45aff2515cb,
title = "Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration",
abstract = "In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ∼6-8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. The combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.",
author = "Alex Albaugh and Mark Tuckerman and Teresa Head-Gordon",
year = "2019",
month = "4",
day = "9",
doi = "10.1021/acs.jctc.9b00072",
language = "English (US)",
volume = "15",
pages = "2195--2205",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "4",

}

TY - JOUR

T1 - Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration

AU - Albaugh, Alex

AU - Tuckerman, Mark

AU - Head-Gordon, Teresa

PY - 2019/4/9

Y1 - 2019/4/9

N2 - In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ∼6-8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. The combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.

AB - In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ∼6-8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. The combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.

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

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

U2 - 10.1021/acs.jctc.9b00072

DO - 10.1021/acs.jctc.9b00072

M3 - Article

VL - 15

SP - 2195

EP - 2205

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 4

ER -