DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

Abstract

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly relatedmore » to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.« less

Authors:
ORCiD logo [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Uppsala Univ. (Sweden)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES); USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1630868
Alternate Identifier(s):
OSTI ID: 1604109
Report Number(s):
LA-UR-19-31455
Journal ID: ISSN 0021-9606; TRN: US2200706
Grant/Contract Number:  
89233218CNA000001; FWP LANLE8AN
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 152; Journal Issue: 10; Conference: 2019 MRS Fall Meeting & Exhibit, Boston, MA (United States), 1-6 Dec 2019; Journal ID: ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; Mathematics; Molecular dynamics; Krylov subspace approximation; Integration methods; Quantum based molecular dynamics; Non-linear systems of equations; Nonlinear systems; Numerical linear algebra; Optimization problems; Newton-Krylov method; Self consistent field methods; Born-Oppenheimer molecular dynamics; Density-matrix

Citation Formats

Niklasson, Anders Mauritz N. Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation. United States: N. p., 2020. Web. doi:10.1063/1.5143270.
Niklasson, Anders Mauritz N. Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation. United States. https://doi.org/10.1063/1.5143270
Niklasson, Anders Mauritz N. Wed . "Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation". United States. https://doi.org/10.1063/1.5143270. https://www.osti.gov/servlets/purl/1630868.
@article{osti_1630868,
title = {Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation},
author = {Niklasson, Anders Mauritz N.},
abstractNote = {It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.},
doi = {10.1063/1.5143270},
journal = {Journal of Chemical Physics},
number = 10,
volume = 152,
place = {United States},
year = {Wed Mar 11 00:00:00 EDT 2020},
month = {Wed Mar 11 00:00:00 EDT 2020}
}

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

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

Save / Share:

Works referenced in this record:

Time-reversible ab initio molecular dynamics
journal, April 2007

  • Niklasson, Anders M. N.; Tymczak, C. J.; Challacombe, Matt
  • The Journal of Chemical Physics, Vol. 126, Issue 14
  • DOI: 10.1063/1.2715556

Energy conserving, linear scaling Born-Oppenheimer molecular dynamics
journal, October 2012

  • Cawkwell, M. J.; Niklasson, Anders M. N.
  • The Journal of Chemical Physics, Vol. 137, Issue 13
  • DOI: 10.1063/1.4755991

Fractional occupations and density-functional energies and forces
journal, June 1992


Hybrid Krylov Methods for Nonlinear Systems of Equations
journal, May 1990

  • Brown, Peter N.; Saad, Youcef
  • SIAM Journal on Scientific and Statistical Computing, Vol. 11, Issue 3
  • DOI: 10.1137/0911026

Approximate Self‐Consistent Molecular‐Orbital Theory. V. Intermediate Neglect of Differential Overlap
journal, September 1967

  • Pople, J. A.; Beveridge, D. L.; Dobosh, P. A.
  • The Journal of Chemical Physics, Vol. 47, Issue 6
  • DOI: 10.1063/1.1712233

Graph-based linear scaling electronic structure theory
journal, June 2016

  • Niklasson, Anders M. N.; Mniszewski, Susan M.; Negre, Christian F. A.
  • The Journal of Chemical Physics, Vol. 144, Issue 23
  • DOI: 10.1063/1.4952650

Efficient and Accurate Born–Oppenheimer Molecular Dynamics for Large Molecular Systems
journal, October 2017

  • Peters, Laurens D. M.; Kussmann, Jörg; Ochsenfeld, Christian
  • Journal of Chemical Theory and Computation, Vol. 13, Issue 11
  • DOI: 10.1021/acs.jctc.7b00937

Next generation extended Lagrangian first principles molecular dynamics
journal, August 2017

  • Niklasson, Anders M. N.
  • The Journal of Chemical Physics, Vol. 147, Issue 5
  • DOI: 10.1063/1.4985893

Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004


The Role and Perspective of Ab Initio Molecular Dynamics in the Study of Biological Systems
journal, June 2002

  • Carloni, Paolo; Rothlisberger, Ursula; Parrinello, Michele
  • Accounts of Chemical Research, Vol. 35, Issue 6
  • DOI: 10.1021/ar010018u

Extended Lagrangian Excited State Molecular Dynamics
journal, January 2018

  • Bjorgaard, J. A.; Sheppard, D.; Tretiak, S.
  • Journal of Chemical Theory and Computation, Vol. 14, Issue 2
  • DOI: 10.1021/acs.jctc.7b00857

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals
journal, June 2001

  • Schlegel, H. Bernhard; Millam, John M.; Iyengar, Srinivasan S.
  • The Journal of Chemical Physics, Vol. 114, Issue 22
  • DOI: 10.1063/1.1372182

Hybrid Extended Lagrangian, Post-Hartree–Fock Born–Oppenheimer ab Initio Molecular Dynamics Using Fragment-Based Electronic Structure
journal, May 2016

  • Li, Junjie; Haycraft, Cody; Iyengar, Srinivasan S.
  • Journal of Chemical Theory and Computation, Vol. 12, Issue 6
  • DOI: 10.1021/acs.jctc.6b00001

Analytic hyperpolarizability and polarizability derivative with fractional occupation numbers for large extended systems
journal, February 2017

  • Nishimoto, Yoshio
  • The Journal of Chemical Physics, Vol. 146, Issue 8
  • DOI: 10.1063/1.4976551

Fock matrix dynamics
journal, March 2004


Dyanmics of organic reactions
journal, November 1973

  • Wang, I. S. Y.; Karplus, M.
  • Journal of the American Chemical Society, Vol. 95, Issue 24
  • DOI: 10.1021/ja00805a033

Generalized extended Lagrangian Born-Oppenheimer molecular dynamics
journal, October 2014

  • Niklasson, Anders M. N.; Cawkwell, Marc J.
  • The Journal of Chemical Physics, Vol. 141, Issue 16
  • DOI: 10.1063/1.4898803

Geometric integration in Born-Oppenheimer molecular dynamics
journal, December 2011

  • Odell, Anders; Delin, Anna; Johansson, Börje
  • The Journal of Chemical Physics, Vol. 135, Issue 22
  • DOI: 10.1063/1.3660689

Stable and Efficient Linear Scaling First-Principles Molecular Dynamics for 10000+ Atoms
journal, November 2014

  • Arita, Michiaki; Bowler, David R.; Miyazaki, Tsuyoshi
  • Journal of Chemical Theory and Computation, Vol. 10, Issue 12
  • DOI: 10.1021/ct500847y

Car-Parrinello molecular dynamics: Car-Parrinello molecular dynamics
journal, September 2011

  • Hutter, Jürg
  • Wiley Interdisciplinary Reviews: Computational Molecular Science, Vol. 2, Issue 4
  • DOI: 10.1002/wcms.90

GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

QM/MM Methods for Biomolecular Systems
journal, January 2009

  • Senn, Hans Martin; Thiel, Walter
  • Angewandte Chemie International Edition, Vol. 48, Issue 7
  • DOI: 10.1002/anie.200802019

Classical trajectories using the full a b i n i t i o potential energy surface H +CH 4 →CH 4 +H
journal, May 1978

  • Leforestier, Claude
  • The Journal of Chemical Physics, Vol. 68, Issue 10
  • DOI: 10.1063/1.435520

Accurate Classical Polarization Solution with No Self-Consistent Field Iterations
journal, April 2017

  • Albaugh, Alex; Niklasson, Anders M. N.; Head-Gordon, Teresa
  • The Journal of Physical Chemistry Letters, Vol. 8, Issue 8
  • DOI: 10.1021/acs.jpclett.7b00450

A note on the Pulay force at finite electronic temperatures
journal, December 2008

  • Niklasson, Anders M. N.
  • The Journal of Chemical Physics, Vol. 129, Issue 24
  • DOI: 10.1063/1.3036203

Thermostating extended Lagrangian Born-Oppenheimer molecular dynamics
journal, April 2015

  • Martínez, Enrique; Cawkwell, Marc J.; Voter, Arthur F.
  • The Journal of Chemical Physics, Vol. 142, Issue 15
  • DOI: 10.1063/1.4917546

Energy versus free-energy conservation in first-principles molecular dynamics
journal, May 1992

  • Wentzcovitch, Renata M.; Martins, José Luís; Allen, Philip B.
  • Physical Review B, Vol. 45, Issue 19
  • DOI: 10.1103/physrevb.45.11372

On the theory underlying the Car-Parrinello method and the role of the fictitious mass parameter
journal, January 2006

  • Tangney, Paul
  • The Journal of Chemical Physics, Vol. 124, Issue 4
  • DOI: 10.1063/1.2162893

An efficient and stable hybrid extended Lagrangian/self-consistent field scheme for solving classical mutual induction
journal, November 2015

  • Albaugh, Alex; Demerdash, Omar; Head-Gordon, Teresa
  • The Journal of Chemical Physics, Vol. 143, Issue 17
  • DOI: 10.1063/1.4933375

Theory of ab initio molecular-dynamics calculations
journal, November 1991


Using Neighbor Diversity to Detect Fraudsters in Online Auctions
journal, May 2014


Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory
journal, March 2017

  • Vitale, Valerio; Dziedzic, Jacek; Albaugh, Alex
  • The Journal of Chemical Physics, Vol. 146, Issue 12
  • DOI: 10.1063/1.4978684

Spin eigenstate-dependent Hartree—Fock molecular dynamics
journal, February 1992


Accelerated, energy-conserving Born–Oppenheimer molecular dynamics via Fock matrix extrapolation
journal, January 2005

  • Herbert, John M.; Head-Gordon, Martin
  • Physical Chemistry Chemical Physics, Vol. 7, Issue 18
  • DOI: 10.1039/b509494a

Ab initio molecular dynamics: basic concepts, current trends and novel applications
journal, December 2002


Higher-Order Extended Lagrangian Born–Oppenheimer Molecular Dynamics for Classical Polarizable Models
journal, January 2018

  • Albaugh, Alex; Head-Gordon, Teresa; Niklasson, Anders M. N.
  • Journal of Chemical Theory and Computation, Vol. 14, Issue 2
  • DOI: 10.1021/acs.jctc.7b01041

Convergence acceleration of iterative sequences. the case of scf iteration
journal, July 1980


GFN2-xTB—An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions
journal, January 2019

  • Bannwarth, Christoph; Ehlert, Sebastian; Grimme, Stefan
  • Journal of Chemical Theory and Computation, Vol. 15, Issue 3
  • DOI: 10.1021/acs.jctc.8b01176

Linear-scaling first-principles molecular dynamics of complex biological systems with the Conquest code
journal, September 2016

  • Otsuka, Takao; Taiji, Makoto; Bowler, David R.
  • Japanese Journal of Applied Physics, Vol. 55, Issue 11
  • DOI: 10.7567/jjap.55.1102b1

Molecular dynamics without effective potentials via the Car-Parrinello approach
journal, August 1990


Extended Lagrangian Born–Oppenheimer molecular dynamics with dissipation
journal, June 2009

  • Niklasson, Anders M. N.; Steneteg, Peter; Odell, Anders
  • The Journal of Chemical Physics, Vol. 130, Issue 21
  • DOI: 10.1063/1.3148075

Canonical density matrix perturbation theory
journal, December 2015

  • Niklasson, Anders M. N.; Cawkwell, M. J.; Rubensson, Emanuel H.
  • Physical Review E, Vol. 92, Issue 6
  • DOI: 10.1103/physreve.92.063301

A New Method for Treating Drude Polarization in Classical Molecular Simulation
journal, October 2017

  • Albaugh, Alex; Head-Gordon, Teresa
  • Journal of Chemical Theory and Computation, Vol. 13, Issue 11
  • DOI: 10.1021/acs.jctc.7b00838

A mathematical investigation of the Car-Parrinello method
journal, January 1998

  • Bornemann, Folkmar A.; Schütte, Christof
  • Numerische Mathematik, Vol. 78, Issue 3
  • DOI: 10.1007/s002110050316

How well do Car–Parrinello simulations reproduce the Born–Oppenheimer surface? Theory and examples
journal, January 2002

  • Tangney, P.; Scandolo, S.
  • The Journal of Chemical Physics, Vol. 116, Issue 1
  • DOI: 10.1063/1.1423331

Canonical-ensemble extended Lagrangian Born–Oppenheimer molecular dynamics for the linear scaling density functional theory
journal, September 2017

  • Hirakawa, Teruo; Suzuki, Teppei; Bowler, David R.
  • Journal of Physics: Condensed Matter, Vol. 29, Issue 40
  • DOI: 10.1088/1361-648x/aa810d

Broyden's method for self-consistent field convergence acceleration
journal, April 1984


Lagrangian formulation with dissipation of Born-Oppenheimer molecular dynamics using the density-functional tight-binding method
journal, July 2011

  • Zheng, Guishan; Niklasson, Anders M. N.; Karplus, Martin
  • The Journal of Chemical Physics, Vol. 135, Issue 4
  • DOI: 10.1063/1.3605303

DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB)
journal, March 2011

  • Gaus, Michael; Cui, Qiang; Elstner, Marcus
  • Journal of Chemical Theory and Computation, Vol. 7, Issue 4
  • DOI: 10.1021/ct100684s

Higher-order symplectic integration in Born–Oppenheimer molecular dynamics
journal, December 2009

  • Odell, Anders; Delin, Anna; Johansson, Börje
  • The Journal of Chemical Physics, Vol. 131, Issue 24
  • DOI: 10.1063/1.3268338

Extended Lagrangian Density Functional Tight-Binding Molecular Dynamics for Molecules and Solids
journal, June 2015

  • Aradi, Bálint; Niklasson, Anders M. N.; Frauenheim, Thomas
  • Journal of Chemical Theory and Computation, Vol. 11, Issue 7
  • DOI: 10.1021/acs.jctc.5b00324

First principles molecular dynamics without self-consistent field optimization
journal, January 2014

  • Souvatzis, Petros; Niklasson, Anders M. N.
  • The Journal of Chemical Physics, Vol. 140, Issue 4
  • DOI: 10.1063/1.4862907

Recursive Factorization of the Inverse Overlap Matrix in Linear-Scaling Quantum Molecular Dynamics Simulations
journal, June 2016

  • Negre, Christian F. A.; Mniszewski, Susan M.; Cawkwell, Marc J.
  • Journal of Chemical Theory and Computation, Vol. 12, Issue 7
  • DOI: 10.1021/acs.jctc.6b00154

Real-World Predictions from Ab Initio Molecular Dynamics Simulations
book, January 2011

  • Kirchner, Barbara; di Dio, Philipp J.; Hutter, Jürg
  • Multiscale Molecular Methods in Applied Chemistry
  • DOI: 10.1007/128_2011_195

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. II. Generalizations based on mass-weighting, idempotency, energy conservation and choice of initial conditions
journal, January 2001

  • Iyengar, Srinivasan S.; Schlegel, H. Bernhard; Millam, John M.
  • The Journal of Chemical Physics, Vol. 115, Issue 22
  • DOI: 10.1063/1.1416876

An extended-Lagrangian scheme for charge equilibration in reactive molecular dynamics simulations
journal, July 2015

  • Nomura, Ken-ichi; Small, Patrick E.; Kalia, Rajiv K.
  • Computer Physics Communications, Vol. 192
  • DOI: 10.1016/j.cpc.2015.02.023

Approximate Self‐Consistent Molecular Orbital Theory. I. Invariant Procedures
journal, November 1965

  • Pople, J. A.; Santry, D. P.; Segal, G. A.
  • The Journal of Chemical Physics, Vol. 43, Issue 10
  • DOI: 10.1063/1.1701475

Time-reversible always stable predictor-corrector method for molecular dynamics of polarizable molecules
journal, January 2003

  • Kolafa, Ji?�
  • Journal of Computational Chemistry, Vol. 25, Issue 3
  • DOI: 10.1002/jcc.10385

Structural and dynamical properties of hot dense matter by a Thomas-Fermi-Dirac molecular dynamics
journal, September 2006


Iterative Procedures for Nonlinear Integral Equations
journal, October 1965