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

Title: Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Abstract

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ KTK, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

Authors:
ORCiD logo [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1843167
Report Number(s):
LA-UR-20-22452
Journal ID: ISSN 1549-9618
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Chemical Theory and Computation
Additional Journal Information:
Journal Volume: 16; Journal Issue: 6; Journal ID: ISSN 1549-9618
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; Mathematics; Born-Oppenheimer molecular dynamics; First principles electronic structure theory; Hartree-Fock; Density matrix; Extended Lagrangian; Molecular dynamics; Algorithms; Optimization; Approximation; Electrical properties

Citation Formats

Niklasson, Anders Mauritz N. Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics. United States: N. p., 2020. Web. doi:10.1021/acs.jctc.0c00264.
Niklasson, Anders Mauritz N. Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics. United States. https://doi.org/10.1021/acs.jctc.0c00264
Niklasson, Anders Mauritz N. Mon . "Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics". United States. https://doi.org/10.1021/acs.jctc.0c00264. https://www.osti.gov/servlets/purl/1843167.
@article{osti_1843167,
title = {Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics},
author = {Niklasson, Anders Mauritz N.},
abstractNote = {Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ KTK, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.},
doi = {10.1021/acs.jctc.0c00264},
journal = {Journal of Chemical Theory and Computation},
number = 6,
volume = 16,
place = {United States},
year = {Mon May 04 00:00:00 EDT 2020},
month = {Mon May 04 00:00:00 EDT 2020}
}

Works referenced in this record:

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

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

Quantum Theory of the Dielectric Constant in Real Solids
journal, April 1962


On the existence of the optimal order for wavefunction extrapolation in Born-Oppenheimer molecular dynamics
journal, June 2016

  • Fang, Jun; Gao, Xingyu; Song, Haifeng
  • The Journal of Chemical Physics, Vol. 144, Issue 24
  • DOI: 10.1063/1.4954234

Linear-scaling density-functional theory with tens of thousands of atoms: Expanding the scope and scale of calculations with ONETEP
journal, July 2009

  • Hine, N. D. M.; Haynes, P. D.; Mostofi, A. A.
  • Computer Physics Communications, Vol. 180, Issue 7
  • DOI: 10.1016/j.cpc.2008.12.023

DFTB+, a Sparse Matrix-Based Implementation of the DFTB Method
journal, July 2007

  • Aradi, B.; Hourahine, B.; Frauenheim, Th.
  • The Journal of Physical Chemistry A, Vol. 111, Issue 26
  • DOI: 10.1021/jp070186p

Kernel Polynomial Approximations for Densities of States and Spectral Functions
journal, March 1996

  • Silver, R. N.; Roeder, H.; Voter, A. F.
  • Journal of Computational Physics, Vol. 124, Issue 1
  • DOI: 10.1006/jcph.1996.0048

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


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

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


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

Extended Lagrangian Born-Oppenheimer molecular dynamics in the limit of vanishing self-consistent field optimization
journal, December 2013

  • Souvatzis, Petros; Niklasson, Anders M. N.
  • The Journal of Chemical Physics, Vol. 139, Issue 21
  • DOI: 10.1063/1.4834015

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

Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties
journal, September 1998

  • Elstner, M.; Porezag, D.; Jungnickel, G.
  • Physical Review B, Vol. 58, Issue 11, p. 7260-7268
  • DOI: 10.1103/PhysRevB.58.7260

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

\mathcal{O}(N) methods in electronic structure calculations
journal, February 2012


Application of localized molecular orbitals to the solution of semiempirical self-consistent field equations
journal, January 1996


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

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

Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach
journal, April 2005

  • VandeVondele, Joost; Krack, Matthias; Mohamed, Fawzi
  • Computer Physics Communications, Vol. 167, Issue 2
  • DOI: 10.1016/j.cpc.2004.12.014

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

Ab initiomolecular dynamics for liquid metals
journal, January 1993


Analysis of Time Reversible Born-Oppenheimer Molecular Dynamics
journal, December 2013

  • Lin, Lin; Lu, Jianfeng; Shao, Sihong
  • Entropy, Vol. 16, Issue 1
  • DOI: 10.3390/e16010110

Ab initio molecular-dynamics techniques extended to large-length-scale systems
journal, January 1992


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

Unified Approach for Molecular Dynamics and Density-Functional Theory
journal, November 1985


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

Numerical Optimization of Density Functional Tight Binding Models: Application to Molecules Containing Carbon, Hydrogen, Nitrogen, and Oxygen
journal, November 2017

  • Krishnapriyan, A.; Yang, P.; Niklasson, A. M. N.
  • Journal of Chemical Theory and Computation, Vol. 13, Issue 12
  • DOI: 10.1021/acs.jctc.7b00762

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

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

Kohn−Sham Density Functional Theory Electronic Structure Calculations with Linearly Scaling Computational Time and Memory Usage
journal, December 2010

  • Rudberg, Elias; Rubensson, Emanuel H.; Sałek, Paweł
  • Journal of Chemical Theory and Computation, Vol. 7, Issue 2
  • DOI: 10.1021/ct100611z

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


Dielectric Constant with Local Field Effects Included
journal, January 1963


Efficient Linear Scaling Algorithm for Tight-Binding Molecular Dynamics
journal, July 1994


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

General atomic and molecular electronic structure system
journal, November 1993

  • Schmidt, Michael W.; Baldridge, Kim K.; Boatz, Jerry A.
  • Journal of Computational Chemistry, Vol. 14, Issue 11, p. 1347-1363
  • DOI: 10.1002/jcc.540141112

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

Extended Born-Oppenheimer Molecular Dynamics
journal, March 2008


The SIESTA method for ab initio order- N materials simulation
journal, March 2002

  • Soler, José M.; Artacho, Emilio; Gale, Julian D.
  • Journal of Physics: Condensed Matter, Vol. 14, Issue 11
  • DOI: 10.1088/0953-8984/14/11/302

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

Zur Quantentheorie der Molekeln
journal, January 1927


Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation
journal, March 2020

  • Niklasson, Anders M. N.
  • The Journal of Chemical Physics, Vol. 152, Issue 10
  • DOI: 10.1063/1.5143270

Time-Reversible Born-Oppenheimer Molecular Dynamics
journal, September 2006


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

Recent progress with large-scaleab initio calculations: the CONQUEST code
journal, April 2006

  • Bowler, D. R.; Choudhury, R.; Gillan, M. J.
  • physica status solidi (b), Vol. 243, Issue 5
  • DOI: 10.1002/pssb.200541386

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations
journal, March 2020

  • Hourahine, B.; Aradi, B.; Blum, V.
  • The Journal of Chemical Physics, Vol. 152, Issue 12
  • DOI: 10.1063/1.5143190

Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics
journal, February 2007


Communication: Constrained molecular dynamics for polarizable models
journal, November 2018

  • Coretti, Alessandro; Bonella, Sara; Ciccotti, Giovanni
  • The Journal of Chemical Physics, Vol. 149, Issue 19
  • DOI: 10.1063/1.5055704

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

New Developments in Molecular Orbital Theory
journal, April 1951


Extended Lagrangian free energy molecular dynamics
journal, October 2011

  • Niklasson, Anders M. N.; Steneteg, Peter; Bock, Nicolas
  • The Journal of Chemical Physics, Vol. 135, Issue 16
  • DOI: 10.1063/1.3656977

Molecular orbital propagation to accelerate self-consistent-field convergence in an ab initio molecular dynamics simulation
journal, March 2008

  • Atsumi, Teruo; Nakai, Hiromi
  • The Journal of Chemical Physics, Vol. 128, Issue 9
  • DOI: 10.1063/1.2839857

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

Wechselwirkung neutraler Atome und hom�opolare Bindung nach der Quantenmechanik
journal, June 1927

  • Heitler, W.; London, F.
  • Zeitschrift f�r Physik, Vol. 44, Issue 6-7
  • DOI: 10.1007/BF01397394

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

Some Recent Advances in Density Matrix Theory
journal, April 1960


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

Stability of the thermal Hartree-Fock approximation
journal, January 1963


Iterative Procedures for Nonlinear Integral Equations
journal, October 1965


Hybrid functionals for large periodic systems in an all-electron, numeric atom-centered basis framework
journal, July 2015

  • Levchenko, Sergey V.; Ren, Xinguo; Wieferink, Jürgen
  • Computer Physics Communications, Vol. 192
  • DOI: 10.1016/j.cpc.2015.02.021