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

The Journal of Chemical Physics

Thermal Gaussian molecular dynamics for quantum dynamics simulations of many-body systems: Application to liquid para-hydrogen

Ionuţ Georgescu, Jason Deckman, Laura Fredrickson, Vladimir A. Mandelshtam

Published May 5, 2011• 25 citations
Publisher

Abstract

A new method, here called thermal Gaussian molecular dynamics (TGMD), for simulating the dynamics of quantum many-body systems has recently been introduced [I. Georgescu and V. A. Mandelshtam, Phys. Rev. B 82, 094305 (2010)]. As in the centroid molecular dynamics (CMD), in TGMD the N-body quantum system is mapped to an N-body classical system. The associated both effective Hamiltonian and effective force are computed within the variational Gaussian wave-packet approximation. The TGMD is exact for the high-temperature limit, accurate for short times, and preserves the quantum canonical distribution. For a harmonic potential and any form of operator \documentclass[12pt]{minimal}\begin{document}$\hat{A}$\end{document}Â, it provides exact time correlation functions CAB(t) at least for the case of \documentclass[12pt]{minimal}\begin{document}$\hat{B}$\end{document}B̂, a linear combination of the position, \documentclass[12pt]{minimal}\begin{document}$\hat{x}$\end{document}x̂, and momentum, \documentclass[12pt]{minimal}\begin{document}$\hat{p}$\end{document}p̂, operators. While conceptually similar to CMD and other quantum molecular dynamics approaches, the great advantage of TGMD is its computational efficiency. We introduce the many-body implementation and demonstrate it on the benchmark problem of calculating the velocity time auto-correlation function for liquid para-hydrogen, using a system of up to N = 2592 particles.

DOI
10.1063/1.3585648

Related papers

Molecular dynamics with quantum fluctuations

Ionuţ Georgescu, Vladimir A. Mandelshtam

A quantum dynamics approach, called Gaussian molecular dynamics, is introduced. As in the centroid molecular dynamics, the $N$-body quantum system is mapped to an $N$-body classical system with an effective Hamiltonian arising within the variational Gaussian wave-packet approximation. The approach is exact for the harmonic oscillator and for the high-temperature limit, accurate in the short-time limit and is computationally very efficient.

Physical Review B•Sep 9, 2010•12 citations

Quantum Molecular Dynamics

Riccardo Fantoni

We formulate a new quantum many body simulation method for a general quantum fluid at any given temperature. Unlike path integral Monte Carlo our method evolves in imaginary time the density matrix from its initial delta function condition to its final thermal form at time time equal to the inverse temperature. It does this with a molecular dynamics scheme on a classical Hamiltonian that has the same functional form as the one of the quantum mechanical Hamiltonian according to the properties of the continuous representation of John R. Klauder. We then end up with the thermal density matrix that can be used to extract thermal averages of observables using the Monte Carlo method equally well in any statistics.

arXiv:2606.01389•Quantum Reports•Feb 6, 2026•0 citations

Representing the thermal state in time-dependent density functional theory

Normand A. Modine, Ryan Hatcher

Classical molecular dynamics (MD) provides a powerful and widely used approach to determining thermodynamic properties by integrating the classical equations of motion of a system of atoms. Time-Dependent Density Functional Theory (TDDFT) provides a powerful and increasingly useful approach to integrating the quantum equations of motion for a system of electrons. TDDFT efficiently captures the unitary evolution of a many-electron state by mapping the system into a fictitious non-interacting system. In analogy to MD, one could imagine obtaining the thermodynamic properties of an electronic system from a TDDFT simulation in which the electrons are excited from their ground state by a time-dependent potential and then allowed to evolve freely in time while statistical data are captured from periodic snapshots of the system. For a variety of systems (e.g., many metals), the electrons reach an effective state of internal equilibrium due to electron-electron interactions on a time scale that is short compared to electron-phonon equilibration. During the initial time-evolution of such systems following electronic excitation, electron-phonon interactions should be negligible, and therefore, TDDFT should successfully capture the internal thermalization of the electrons. However, it is unclear how TDDFT represents the resulting thermal state. In particular, the thermal state is usually represented in quantum statistical mechanics as a mixed state, while the occupations of the TDDFT wavefunctions are fixed by the initial state in TDDFT. We work to address this puzzle by (A) reformulating quantum statistical mechanics so that thermodynamic expectations can be obtained as an unweighted average over a set of many-body pure states and (B) constructing a family of non-interacting (single determinant) TDDFT states that approximate the required many-body states for the canonical ensemble.

The Journal of Chemical Physics•May 28, 2015•11 citations

The formulation of quantum statistical mechanics based on the Feynman path centroid density. IV. Algorithms for centroid molecular dynamics

Jianshu Cao, Gregory A. Voth

Numerical algorithms are developed for the centroid molecular dynamics (centroid MD) method to calculate dynamical time correlation functions for general many-body quantum systems. Approaches based on the normal mode path integral molecular dynamics and staging path integral Monte Carlo methods are described to carry out a direct calculation of the force on the centroid variables in the centroid MD algorithm. A more efficient, but approximate, scheme to compute the centroid force is devised which is based on the locally optimized harmonic approximation for the centroid potential. The centroid MD equations in the latter method can be solved with the help of an iterative procedure or through extended Lagrangian dynamics. A third algorithm introduces an effective centroid pseudopotential to approximate the full many-body centroid mean force potential by effective pairwise centroid interactions. Numerical simulations for both prototype models and more realistic many-body systems are performed to explore the feasibility and limitations of each algorithm.

The Journal of Chemical Physics•Oct 1, 1994•379 citations

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

Scott Habershon

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

The Journal of Chemical Physics•Sep 12, 2013•8 citations

"Almost classical" many-body systems: The quantum-mechanical corrections to the moments of a general spectrum

F. Barocchi, Massimo Moraldi, Marco Zoppi

The Wigner method in quantum statistical mechanics, which allows the derivation of the quantum-mechanical behavior of the properties of an "almost classical" system, has been applied to the determination of series expansions with respect to $h$ for a general correlation function and the moments related to its spectrum. Explicit expressions of the terms up to the order ${h}^{4}$ are given for the zeroth, second, and fourth moment in the case of a many-body system with a Hamiltonian $\mathit{H}=\frac{{\mathit{P}}^{2}}{2m}+\ensuremath{\Phi}(\mathit{R})$ and for variables which are functions only either of ${\mathit{P}}$ or ${\mathit{R}}$ coordinates. From these expressions the corrections to the classical behavior can be calculated via the classical molecular-dynamics simulation technique.

Physical Review A•Oct 1, 1982•33 citations
A new method, here called thermal Gaussian molecular dynamics (TGMD), for simulating the dynamics of quantum many-body systems has recently been introduced [I. Georgescu and V. A. Mandelshtam, Phys. Rev. B 82, 094305 (2010)]. As in the centroid molecular dynamics (CMD), in TGMD the N-body quantum system is mapped to an N-body classical system. The associated both effective Hamiltonian and effective force are computed within the variational Gaussian wave-packet approximation. The TGMD is exact for the high-temperature limit, accurate for short times, and preserves the quantum canonical distribution. For a harmonic potential and any form of operator \documentclass[12pt]{minimal}\begin{document}A^\hat{A}A^\end{document}Â, it provides exact time correlation functions CAB(t) at least for the case of \documentclass[12pt]{minimal}\begin{document}B^\hat{B}B^\end{document}B̂, a linear combination of the position, \documentclass[12pt]{minimal}\begin{document}x^\hat{x}x^\end{document}x̂, and momentum, \documentclass[12pt]{minimal}\begin{document}p^\hat{p}p^​\end{document}p̂, operators. While conceptually similar to CMD and other quantum molecular dynamics approaches, the great advantage of TGMD is its computational efficiency. We introduce the many-body implementation and demonstrate it on the benchmark problem of calculating the velocity time auto-correlation function for liquid para-hydrogen, using a system of up to N = 2592 particles.
A quantum dynamics approach, called Gaussian molecular dynamics, is introduced. As in the centroid molecular dynamics, the NNN-body quantum system is mapped to an NNN-body classical system with an effective Hamiltonian arising within the variational Gaussian wave-packet approximation. The approach is exact for the harmonic oscillator and for the high-temperature limit, accurate in the short-time limit and is computationally very efficient.
We formulate a new quantum many body simulation method for a general quantum fluid at any given temperature. Unlike path integral Monte Carlo our method evolves in imaginary time the density matrix from its initial delta function condition to its final thermal form at time time equal to the inverse temperature. It does this with a molecular dynamics scheme on a classical Hamiltonian that has the same functional form as the one of the quantum mechanical Hamiltonian according to the properties of the continuous representation of John R. Klauder. We then end up with the thermal density matrix that can be used to extract thermal averages of observables using the Monte Carlo method equally well in any statistics.
Classical molecular dynamics (MD) provides a powerful and widely used approach to determining thermodynamic properties by integrating the classical equations of motion of a system of atoms. Time-Dependent Density Functional Theory (TDDFT) provides a powerful and increasingly useful approach to integrating the quantum equations of motion for a system of electrons. TDDFT efficiently captures the unitary evolution of a many-electron state by mapping the system into a fictitious non-interacting system. In analogy to MD, one could imagine obtaining the thermodynamic properties of an electronic system from a TDDFT simulation in which the electrons are excited from their ground state by a time-dependent potential and then allowed to evolve freely in time while statistical data are captured from periodic snapshots of the system. For a variety of systems (e.g., many metals), the electrons reach an effective state of internal equilibrium due to electron-electron interactions on a time scale that is short compared to electron-phonon equilibration. During the initial time-evolution of such systems following electronic excitation, electron-phonon interactions should be negligible, and therefore, TDDFT should successfully capture the internal thermalization of the electrons. However, it is unclear how TDDFT represents the resulting thermal state. In particular, the thermal state is usually represented in quantum statistical mechanics as a mixed state, while the occupations of the TDDFT wavefunctions are fixed by the initial state in TDDFT. We work to address this puzzle by (A) reformulating quantum statistical mechanics so that thermodynamic expectations can be obtained as an unweighted average over a set of many-body pure states and (B) constructing a family of non-interacting (single determinant) TDDFT states that approximate the required many-body states for the canonical ensemble.
Numerical algorithms are developed for the centroid molecular dynamics (centroid MD) method to calculate dynamical time correlation functions for general many-body quantum systems. Approaches based on the normal mode path integral molecular dynamics and staging path integral Monte Carlo methods are described to carry out a direct calculation of the force on the centroid variables in the centroid MD algorithm. A more efficient, but approximate, scheme to compute the centroid force is devised which is based on the locally optimized harmonic approximation for the centroid potential. The centroid MD equations in the latter method can be solved with the help of an iterative procedure or through extended Lagrangian dynamics. A third algorithm introduces an effective centroid pseudopotential to approximate the full many-body centroid mean force potential by effective pairwise centroid interactions. Numerical simulations for both prototype models and more realistic many-body systems are performed to explore the feasibility and limitations of each algorithm.
We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.
The Wigner method in quantum statistical mechanics, which allows the derivation of the quantum-mechanical behavior of the properties of an "almost classical" system, has been applied to the determination of series expansions with respect to hhh for a general correlation function and the moments related to its spectrum. Explicit expressions of the terms up to the order h4{h}^{4}h4 are given for the zeroth, second, and fourth moment in the case of a many-body system with a Hamiltonian $\mathit{H}=\frac{{\mathit{P}}^{2}}{2m}+\ensuremath{\Phi}(\mathit{R})$ and for variables which are functions only either of P{\mathit{P}}P or R{\mathit{R}}R coordinates. From these expressions the corrections to the classical behavior can be calculated via the classical molecular-dynamics simulation technique.