skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory

Abstract

Iterative energy minimization with the aim of achieving self-consistency is a common feature of Born-Oppenheimer molecular dynamics (BOMD) and classical molecular dynamics with polarizable force fields. In the former, the electronic degrees of freedom are optimized, while the latter often involves an iterative determination of induced point dipoles. The computational effort of the self-consistency procedure can be reduced by re-using converged solutions from previous time steps. However, this must be done carefully, as not to break time-reversal symmetry, which negatively impacts energy conservation. Self-consistent schemes based on the extended Lagrangian formalism, where the initial guesses for the optimized quantities are treated as auxiliary degrees of freedom, constitute one elegant solution. We report on the performance of two integration schemes with the same underlying extended Lagrangian structure, which we both employ in two radically distinct regimes—in classical molecular dynamics simulations with the AMOEBA polarizable force field and in BOMD simulations with the Onetep linear-scaling density functional theory (LS-DFT) approach. Furthermore, both integration schemes are found to offer significant improvements over the standard (unpropagated) molecular dynamics formulation in both the classical and LS-DFT regimes.

Authors:
 [1];  [2];  [3]; ORCiD logo [4];  [5];  [1]
  1. Univ. of Southampton, Southampton (United Kingdom)
  2. Univ. of Southampton, Southampton (United Kingdom); Gdansk Univ. of Technology, Gdansk (Poland)
  3. Univ. of California, Berkeley, CA (United States)
  4. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  5. Univ. of California, Berkeley, CA (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC). Basic Energy Sciences (BES) (SC-22); USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
1414137
Alternate Identifier(s):
OSTI ID: 1476472
Report Number(s):
LA-UR-17-20731
Journal ID: ISSN 0021-9606; TRN: US1800661
Grant/Contract Number:  
AC52-06NA25396; AC02-05CH11231
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 146; Journal Issue: 12; Journal ID: ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; Computer Science; Mathematics; Material Science; molecular dynamics; polarizable force fields; density functional theory

Citation Formats

Vitale, Valerio, Dziedzic, Jacek, Albaugh, Alex, Niklasson, Anders M. N., Head-Gordon, Teresa, and Skylaris, Chris -Kriton. Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory. United States: N. p., 2017. Web. doi:10.1063/1.4978684.
Vitale, Valerio, Dziedzic, Jacek, Albaugh, Alex, Niklasson, Anders M. N., Head-Gordon, Teresa, & Skylaris, Chris -Kriton. Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory. United States. doi:10.1063/1.4978684.
Vitale, Valerio, Dziedzic, Jacek, Albaugh, Alex, Niklasson, Anders M. N., Head-Gordon, Teresa, and Skylaris, Chris -Kriton. Tue . "Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory". United States. doi:10.1063/1.4978684. https://www.osti.gov/servlets/purl/1414137.
@article{osti_1414137,
title = {Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory},
author = {Vitale, Valerio and Dziedzic, Jacek and Albaugh, Alex and Niklasson, Anders M. N. and Head-Gordon, Teresa and Skylaris, Chris -Kriton},
abstractNote = {Iterative energy minimization with the aim of achieving self-consistency is a common feature of Born-Oppenheimer molecular dynamics (BOMD) and classical molecular dynamics with polarizable force fields. In the former, the electronic degrees of freedom are optimized, while the latter often involves an iterative determination of induced point dipoles. The computational effort of the self-consistency procedure can be reduced by re-using converged solutions from previous time steps. However, this must be done carefully, as not to break time-reversal symmetry, which negatively impacts energy conservation. Self-consistent schemes based on the extended Lagrangian formalism, where the initial guesses for the optimized quantities are treated as auxiliary degrees of freedom, constitute one elegant solution. We report on the performance of two integration schemes with the same underlying extended Lagrangian structure, which we both employ in two radically distinct regimes—in classical molecular dynamics simulations with the AMOEBA polarizable force field and in BOMD simulations with the Onetep linear-scaling density functional theory (LS-DFT) approach. Furthermore, both integration schemes are found to offer significant improvements over the standard (unpropagated) molecular dynamics formulation in both the classical and LS-DFT regimes.},
doi = {10.1063/1.4978684},
journal = {Journal of Chemical Physics},
number = 12,
volume = 146,
place = {United States},
year = {Tue Mar 28 00:00:00 EDT 2017},
month = {Tue Mar 28 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 5 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Semiempirical GGA-type density functional constructed with a long-range dispersion correction
journal, January 2006

  • Grimme, Stefan
  • Journal of Computational Chemistry, Vol. 27, Issue 15, p. 1787-1799
  • DOI: 10.1002/jcc.20495