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

Title: Next generation extended Lagrangian first principles molecular dynamics

ORCiD logo [1]
  1. Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
Publication Date:
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 147; Journal Issue: 5; Related Information: CHORUS Timestamp: 2018-02-14 12:09:51; Journal ID: ISSN 0021-9606
American Institute of Physics
Country of Publication:
United States

Citation Formats

Niklasson, Anders M. N. Next generation extended Lagrangian first principles molecular dynamics. United States: N. p., 2017. Web. doi:10.1063/1.4985893.
Niklasson, Anders M. N. Next generation extended Lagrangian first principles molecular dynamics. United States. doi:10.1063/1.4985893.
Niklasson, Anders M. N. 2017. "Next generation extended Lagrangian first principles molecular dynamics". United States. doi:10.1063/1.4985893.
title = {Next generation extended Lagrangian first principles molecular dynamics},
author = {Niklasson, Anders M. N.},
abstractNote = {},
doi = {10.1063/1.4985893},
journal = {Journal of Chemical Physics},
number = 5,
volume = 147,
place = {United States},
year = 2017,
month = 8

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on August 4, 2018
Publisher's Accepted Manuscript

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

Save / Share:
  • An extended first-principles molecular dynamics (FPMD) method based on Kohn-Sham scheme is proposed to elevate the temperature limit of the FPMD method in the calculation of dense plasmas. The extended method treats the wave functions of high energy electrons as plane waves analytically and thus expands the application of the FPMD method to the region of hot dense plasmas without suffering from the formidable computational costs. In addition, the extended method inherits the high accuracy of the Kohn-Sham scheme and keeps the information of electronic structures. This gives an edge to the extended method in the calculation of mixtures ofmore » plasmas composed of heterogeneous ions, high-Z dense plasmas, lowering of ionization potentials, X-ray absorption/emission spectra, and opacities, which are of particular interest to astrophysics, inertial confinement fusion engineering, and laboratory astrophysics.« less
  • A set of interatomic potentials for hydrocarbons that are based upon the self-consistent charge transfer tight-binding approximation to density functional theory have been developed and implemented into the quantum molecular dynamics code ''LATTE''. The interatomic potentials exhibit an outstanding level of transferability and have been applied in molecular dynamics simulations of tert-butylacetylene under thermodynamic conditions that correspond to its single-shock Hugoniot. We have achieved precise conservation of the total energy during microcanonical molecular dynamics trajectories under incomplete convergence via the extended Lagrangian Born-Oppenheimer molecular dynamics formalism. In good agreement with the results of a series of flyer-plate impact experiments, ourmore » SCC-TB molecular dynamics simulations show that tert-butylactylene molecules polymerize at shock pressures around 6.1 GPa.« less
  • We present a generalization to our previously developed quantum wavepacket ab initio molecular dynamics (QWAIMD) method by using multiple diabatic electronic reduced single particle density matrices, propagated within an extended Lagrangian paradigm. The Slater determinantal wavefunctions associated with the density matrices utilized may be orthogonal or nonorthogonal with respect to each other. This generalization directly results from an analysis of the variance in electronic structure with quantum nuclear degrees of freedom. The diabatic electronic states are treated here as classical parametric variables and propagated simultaneously along with the quantum wavepacket and classical nuclei. Each electronic density matrix is constrained tomore » be N-representable. Consequently two sets of new methods are derived: extended Lagrangian-QWAIMD (xLag-QWAIMD) and diabatic extended Lagrangian-QWAIMD (DxLag-QWAIMD). In both cases, the instantaneous potential energy surface for the quantum nuclear degrees of freedom is constructed from the diabatic states using an on-the-fly nonorthogonal multireference formalism. By introducing generalized grid-based electronic basis functions, we eliminate the basis set dependence on the quantum nucleus. Subsequent reuse of the two-electron integrals during the on-the-fly potential energy surface computation stage yields a substantial reduction in computational costs. Specifically, both xLag-QWAIMD and DxLag-QWAIMD turn out to be about two orders of magnitude faster than our previously developed time-dependent deterministic sampling implementation of QWAIMD. Energy conservation properties, accuracy of the associated potential surfaces, and vibrational properties are analyzed for a family of hydrogen bonded systems.« less
  • We present an efficient general approach to first principles molecular dynamics simulations based on extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The reduction of the optimization requirement reduces the computational cost to a minimum, but without causing any significant loss of accuracy or long-term energy drift. The optimization-free first principles molecular dynamics requires only one single diagonalization per time step, but is still able to provide trajectories at the same level of accuracy as “exact,” fully converged, Born-Oppenheimer molecular dynamics simulations. The optimization-free limitmore » of extended Lagrangian Born-Oppenheimer molecular dynamics therefore represents an ideal starting point for robust and efficient first principles quantum mechanical molecular dynamics simulations.« less