Time-reversible Born-Oppenheimer molecular dynamics

We present a new method which combines Car-Parrinello and Born-Oppenheimer molecular dynamics in order to accelerate density functional theory based ab-initio simulations. Depending on the system a gain in efficiency of one to two orders of magnitude has been observed, which allows ab-initio molecular dynamics of much larger time and length scales than previously thought feasible. It will be demonstrated that the dynamics is correctly reproduced and that high accuracy can be ...

Extended Lagrangian Born-Oppenheimer molecular dynamics [{\em Phys.\ Rev.\ Lett.\ } {\bf 2008}, {\em 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 evalua...

With the continuous growth of processing power for scientific computing, first principles Born-Oppenheimer molecular dynamics (MD) simulations are becoming increasingly popular for the study of a wide range of problems in materials science, chemistry and biology. Nevertheless, the computational cost still remains prohibitively large in many cases, particularly in comparison to classical MD simulations using empirical force fields. Here we show how to circumvent the major comp...

This report presents a new approach for treating the coupling of electrons and nuclei in quantum mechanical calculations for molecules and condensed matter. It includes the standard "Born-Oppenheimer approximation" as a special case but treats both adiabatic and non-adiabatic corrections using perturbation theory. The adiabatic corrections include all terms that do not explicitly involve the nuclear wavefunctions, so that the nuclei move on a single electronic potential surfa...

Computer simulation methods, such as Monte Carlo or Molecular Dynamics, are very powerful computational techniques that provide detailed and essentially exact information on classical many-body problems. With the advent of ab-initio molecular dynamics, where the forces are computed on-the-fly by accurate electronic structure calculations, the scope of either method has been greatly extended. This new approach, which unifies Newton's and Schr\"odinger's equations, allows for c...

Extended Lagrangian Born-Oppenheimer molecular dynamics [Phys. Rev. Lett., ${\bf 100}$, 123004 (2008)] is formulated for general Hohenberg-Kohn density functional theory and compared to the extended Lagrangian framework of first principles molecular dynamics by Car and Parrinello [Phys. Rev. Lett. ${\bf 55}$, 2471 (1985)]. It is shown how extended Lagrangian Born-Oppenheimer molecular dynamics overcomes several shortcomings of regular, direct Born-Oppenheimer molecular dynami...

Ehrenfest, Born-Oppenheimer, Langevin and Smoluchowski dynamics are shown to be accurate approximations of time-independent Schr\"odinger observables for a molecular system avoiding caustics, in the limit of large ratio of nuclei and electron masses, without assuming that the nuclei are localized to vanishing domains. The derivation, based on a Hamiltonian system interpretation of the Schr\"odinger equation and stability of the corresponding Hamilton-Jacobi equation, bypasses...

In this work, we describe various improved implementations of the mapping approach to surface hopping (MASH) for simulating nonadiabatic dynamics. These include time-reversible and piecewise-continuous integrators, which is only formally possible because of the deterministic nature of the underlying MASH equations of motion. The new algorithms allow for the use of either wave-function overlaps or nonadiabatic coupling vectors to propagate the spin, which encodes the electroni...

Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electro...

We present in detail the recently derived ab-initio molecular dynamics (AIMD) formalism [Phys. Rev. Lett. 101 096403 (2008)], which due to its numerical properties, is ideal for simulating the dynamics of systems containing thousands of atoms. A major drawback of traditional AIMD methods is the necessity to enforce the orthogonalization of the wave-functions, which can become the bottleneck for very large systems. Alternatively, one can handle the electron-ion dynamics within...