Time-reversible Born-Oppenheimer molecular dynamics

In Born-Oppenheimer molecular dynamics (BOMD) simulations based on density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which in practice never is fully converged. The calculated energies and the forces are therefore only approximate, which may lead to an unphysical energy drift and instabilities. Here we discuss a...

In principle, we should not need the time-dependent extension of density-functional theory (TDDFT) for excitations, and in particular not for Molecular Dynamics (MD) studies: the theorem by Hohenberg and Kohn teaches us that for any observable that we wish to look at (including dynamical properties or observables dependent on excited states) there is a corresponding functional of the ground-state density. Yet the unavailability of such magic functionals in many cases (the the...

We have derived equations for nonadiabatic Ehrenfest molecular dynamics which conserve the total energy in the case of time-dependent discretization for electrons. A discretization is time-dependent in all cases where it or part of it depends on the positions of the nuclei, for example, in atomic orbital basis sets, and in the projector augmented-wave (PAW) method, where the augmentation functions depend on the nuclear positions. We have derived, implemented, and analyzed the...

In this paper we present the Uppsala Quantum Chemistry package (UQUANTCHEM), a new and versatile computational platform with capabilities ranging from simple Hartree-Fock calculations to state of the art First principles Extended Lagrangian Born Oppenheimer Molecular Dynamics (XL- BOMD) and diffusion quantum Monte Carlo (DMC). The UQUANTCHEM package is distributed under the general public license and can be directly downloaded from the code web-site. Together with a presentat...

We prove that for a combined system of classical and quantum particles, it is possible to write a dynamics for the classical particles that incorporates in a natural way the Boltzmann equilibrium population for the quantum subsystem. In addition, these molecular dynamics do not need to assume that the electrons immediately follow the nuclear motion (in contrast to any adiabatic approach), and do not present problems in the presence of crossing points between different potenti...

Mixed-quantum-classical molecular dynamics simulation implies an effective measurement on the electronic states owing to continuously tracking the atomic forces.Based on this insight, we propose a quantum trajectory mean-field approach for nonadiabatic molecular dynamics simulations. The new protocol provides a natural interface between the separate quantum and classical treatments, without invoking artificial surface hopping algorithm. Moreover, it also bridges two widely ad...

It is shown that the adiabatic Born-Oppenheimer expansion does not satisfy the necessary condition for the applicability of perturbation theory. A simple example of an exact solution of a problem that can not be obtained from the Born-Oppenheimer expansion is given. A new version of perturbation theory for molecular systems is proposed.

This article has been published as a chapter in "Chemical Reactivity Theory: A Density Functional View", ed. P. K. Chattaraj (CRC Press, New York, 2009), ch. 8, p. 105. In it, an overview of the relationship between time-dependent DFT and quantum hydrodynamics is presented, showing the role that Bohmian mechanics can play within the ab-initio methodology as both a numerical and an interpretative tool.

Severe methodological and numerical problems of the traditional quantum mechanical approach to the description of molecular systems are outlined. To overcome these, a simple alternative to the Born-Oppenheimer approximation is presented on the basis of taking the nuclei as classical particles.

In our previous paper [J. Chem.Phys. {\bf 137}, 22A544 (2012)] we argued that the Born-Oppenheimer approximation could not be based on an exact transformation of the molecular Schr\"{o}dinger equation. In this Comment we suggest that the fundamental reason for the approximate nature of the Born-Oppenheimer model is the lack of a complete set of functions for the electronic space, and the need to describe the continuous spectrum using spectral projection.