Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2609.03419

In defence of the Ehrenfest mean-field molecular dynamics

Vladimir U. Nazarov

quant-phSubmitted Sep 3, 2026• Announced Sep 4, 2026
v1 · Sep 3, 2026v2 · Sep 5, 2026
arXivHTMLPDF

Abstract

Re-visiting the problem of the coupled electronic--nuclear motion, we prove that the so called `mean-field' Ehrenfest picture is an exact theory, rather than an approximation, in the classical limit for nuclei. By establishing its full equivalence to the much more involved Exact Factorization approach to the problem (Abedi {\it et al.}, Europhys. Lett. {\bf 106}, 33001 (2014)), our conclusion rehabilitates the mean-field Ehrenfest molecular dynamics as a simple, effective, while formally exact scheme, amenable to the use within the framework of the time-dependent density functional theory.

Comments
10 pages

Related papers

Time-dependent mean-field approximation for nuclear dynamical problems

Shimon Levit

A general method for a microscopic description of nuclear dynamical problems is discussed. The method is based on a functional integral representation for the many-body time evolution operator $U({t}_{f},{t}_{i})$. In the stationary phase limit, a time-dependent mean-field approximation for the matrix elements of $U$ is obtained. Using standard procedures this allows for extraction of quantum mechanical information about bound states, tunneling, or scattering phenomena in a many-nucleon system. The approximation represents a natural generalization of the time-dependent Hartree-Fock method.NUCLEAR STRUCTURE, NUCLEAR REACTIONS derived mean-field approximation for matrix elements of $\mathrm{exp}(\ensuremath{-}iHt)$ from exact functional integral; discussed relation to static HF, TDHF, and scattering calculations.

Physical Review C•Apr 1, 1980•64 citations

Stochastic and deterministic molecular dynamics derived from the time-independent Schrödinger equation

Anders Szepessy

Ehrenfest, Born-Oppenheimer, Langevin and Smoluchowski dynamics are shown to be accurate approximations of time-independent Schrödinger 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ödinger equation and stability of the corresponding Hamilton-Jacobi equation, bypasses the usual separation of nuclei and electron wave functions, includes crossing electron eigenvalues, and gives a different perspective on the Born-Oppenheimer approximation, Schrödinger Hamiltonian systems, stochastic electron equilibrium states and numerical simulation in molecular dynamics modeling.

arXiv:0812.4338•Dec 23, 2008•3 citations

Quantum dynamics of molecular motions based on a functional method

A. K. Rajagopal, F. A. Buot

A time-dependent functional formalism of molecular positions is developed in order to study non- equilibrium problems. It is based on the Liouvillean quantum field dynamics and is different from the density-functional framework. It provides a direct self-consistent approach to molecular dynamics of a system which is initially in thermodynamic equilibrium and is subjected to time-dependent external force.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Aug 1, 1994•3 citations

Time-Dependent Density Functional Theory from a Bohmian Perspective

Pratim Kumar Chattaraj

Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg, Kohn, and Sham [4,5] in the mid-1960s. According to this theory, the fundamental physical information about a many-body system is provided by single-particle densities in a three-dimensional space, which are obtained variationally within a time-independent framework. When compared with other previous formalisms, DFT presents two clear advantages: (1) it is able to treat many-body problems in a sufficiently accurate way and (2) it is computationally simple. This explains why it is one of the most widely used theories to deal with electronic structure-the electronic ground-state energy as a function of the position of the atomic nuclei determines molecular structures and solids, providing at the same time the forces acting on the atomic nuclei when they are not at their equilibrium positions. At present, DFT is used routinely to solve many problems in gas phase and condensed matter. Furthermore, it has made possible the development of accurate molecular dynamics schemes in which the forces are evaluated quantum mechanically ‘‘on the fly.’’ Nonetheless, DFT is a fundamental tool provided the systems studied are relatively large; for small systems, standard methods based on the use of the wave function render quite accurate results [6]. Moreover, it is also worth stressing that all practical applications of DFT rely on essentially uncontrolled approximations [7] (e.g., the local density approximation [4,5], the local spin-density approximation, or generalized gradient approximations [8]), and therefore the validity of DFT is conditioned to its capability to provide fairly good values of the experimental data.

CRC Press eBooks•Feb 23, 2009•3 citations

Density functional approaches to collective phenomena in nuclei: Time-dependent density functional theory for perturbative and non-perturbative nuclear dynamics

Takashi Nakatsukasa

We present the basic concepts and our recent developments in density functional approaches with Skyrme functionals for describing nuclear dynamics at low energy.Time-dependent density functional theory (TDDFT) is utilized for the exact linear response with an external perturbation.For a description of collective dynamics beyond the perturbative regime, we present a theory of a decoupled collective submanifold to describe a slow motion based on TDDFT.Selected applications are shown to demonstrate the quality of their performance and feasibility.The advantages and disadvantages of the numerical aspects are also discussed.. . . . . .

arXiv:1210.0378•Progress of Theoretical and Experimental Physics•Jan 1, 2012•47 citations

Thawed Gaussian Ehrenfest dynamics

Jiří J. L. Vaníček

Ehrenfest dynamics is a widely used mixed quantum--classical approach for nonadiabatic molecular dynamics, whereas thawed Gaussian wavepacket dynamics provides an efficient semiclassical description of adiabatic nuclear quantum dynamics. Here we describe thawed Gaussian Ehrenfest dynamics (TGED), which unifies and generalizes these two methods to capture both electronic nonadiabaticity and nuclear quantum effects within a single framework. The fully variational formulation of TGED is derived by applying the time-dependent variational principle to a Hartree product of electronic and Gaussian nuclear wavepackets. Replacing the effective locally quadratic molecular potential obtained from this variational treatment by alternative effective locally quadratic potentials yields an infinite family of TGED methods, of which we present several members. We analyze the limiting cases of the general formalism and show, in particular, that it reduces to conventional Ehrenfest dynamics in the classical limit for the nuclei and to thawed Gaussian wavepacket dynamics in the absence of electronic coupling. Finally, we present explicit geometric integrators for the entire family of methods and identify the conditions under which the different approximations become exact.

physics.chem-phCross-listedarXiv:2607.10847•Jul 14, 2026
Re-visiting the problem of the coupled electronic--nuclear motion, we prove that the so called `mean-field' Ehrenfest picture is an exact theory, rather than an approximation, in the classical limit for nuclei. By establishing its full equivalence to the much more involved Exact Factorization approach to the problem (Abedi {\it et al.}, Europhys. Lett. {\bf 106}, 33001 (2014)), our conclusion rehabilitates the mean-field Ehrenfest molecular dynamics as a simple, effective, while formally exact scheme, amenable to the use within the framework of the time-dependent density functional theory.
A general method for a microscopic description of nuclear dynamical problems is discussed. The method is based on a functional integral representation for the many-body time evolution operator U(tf,ti)U({t}_{f},{t}_{i})U(tf​,ti​). In the stationary phase limit, a time-dependent mean-field approximation for the matrix elements of UUU is obtained. Using standard procedures this allows for extraction of quantum mechanical information about bound states, tunneling, or scattering phenomena in a many-nucleon system. The approximation represents a natural generalization of the time-dependent Hartree-Fock method.NUCLEAR STRUCTURE, NUCLEAR REACTIONS derived mean-field approximation for matrix elements of $\mathrm{exp}(\ensuremath{-}iHt)$ from exact functional integral; discussed relation to static HF, TDHF, and scattering calculations.
Ehrenfest, Born-Oppenheimer, Langevin and Smoluchowski dynamics are shown to be accurate approximations of time-independent Schrödinger 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ödinger equation and stability of the corresponding Hamilton-Jacobi equation, bypasses the usual separation of nuclei and electron wave functions, includes crossing electron eigenvalues, and gives a different perspective on the Born-Oppenheimer approximation, Schrödinger Hamiltonian systems, stochastic electron equilibrium states and numerical simulation in molecular dynamics modeling.
Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg, Kohn, and Sham [4,5] in the mid-1960s. According to this theory, the fundamental physical information about a many-body system is provided by single-particle densities in a three-dimensional space, which are obtained variationally within a time-independent framework. When compared with other previous formalisms, DFT presents two clear advantages: (1) it is able to treat many-body problems in a sufficiently accurate way and (2) it is computationally simple. This explains why it is one of the most widely used theories to deal with electronic structure-the electronic ground-state energy as a function of the position of the atomic nuclei determines molecular structures and solids, providing at the same time the forces acting on the atomic nuclei when they are not at their equilibrium positions. At present, DFT is used routinely to solve many problems in gas phase and condensed matter. Furthermore, it has made possible the development of accurate molecular dynamics schemes in which the forces are evaluated quantum mechanically ‘‘on the fly.’’ Nonetheless, DFT is a fundamental tool provided the systems studied are relatively large; for small systems, standard methods based on the use of the wave function render quite accurate results [6]. Moreover, it is also worth stressing that all practical applications of DFT rely on essentially uncontrolled approximations [7] (e.g., the local density approximation [4,5], the local spin-density approximation, or generalized gradient approximations [8]), and therefore the validity of DFT is conditioned to its capability to provide fairly good values of the experimental data.
We present the basic concepts and our recent developments in density functional approaches with Skyrme functionals for describing nuclear dynamics at low energy.Time-dependent density functional theory (TDDFT) is utilized for the exact linear response with an external perturbation.For a description of collective dynamics beyond the perturbative regime, we present a theory of a decoupled collective submanifold to describe a slow motion based on TDDFT.Selected applications are shown to demonstrate the quality of their performance and feasibility.The advantages and disadvantages of the numerical aspects are also discussed.. . . . . .
Ehrenfest dynamics is a widely used mixed quantum--classical approach for nonadiabatic molecular dynamics, whereas thawed Gaussian wavepacket dynamics provides an efficient semiclassical description of adiabatic nuclear quantum dynamics. Here we describe thawed Gaussian Ehrenfest dynamics (TGED), which unifies and generalizes these two methods to capture both electronic nonadiabaticity and nuclear quantum effects within a single framework. The fully variational formulation of TGED is derived by applying the time-dependent variational principle to a Hartree product of electronic and Gaussian nuclear wavepackets. Replacing the effective locally quadratic molecular potential obtained from this variational treatment by alternative effective locally quadratic potentials yields an infinite family of TGED methods, of which we present several members. We analyze the limiting cases of the general formalism and show, in particular, that it reduces to conventional Ehrenfest dynamics in the classical limit for the nuclei and to thawed Gaussian wavepacket dynamics in the absence of electronic coupling. Finally, we present explicit geometric integrators for the entire family of methods and identify the conditions under which the different approximations become exact.