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:2606.01389 · Quantum Reports

Quantum Molecular Dynamics

Riccardo Fantoni

Published Feb 6, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.3390/quantum8010012

Related papers

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

Quantum Monte Carlo Simulation

Arthur McGurn

An introduction is given to the methodology of the path-integral formulation of quantum Monte Carlo simulation. This allows us to reformulate a d −dimensional partition function of a quantum mechanical system into the form of a d +1−dimensional path-integral written in classical mechanical variables. In addition, the quantum mechanical partition function is shown to be related by a process of analytical continuation to the time-development operator of the quantum mechanical system. The evaluation of the path-integral representation is discussed for the Maxwell-Boltzmann ideal gas, a chain of atoms, and liquid 4 He . The properties of liquid Helium I and Helium II are discussed in terms of the two-fluid model, the Ginzburg theory, and recent quantum Monte Carlo Simulation studies. The properties of liquid 3 He are discussed in the context of recent analytical p-wave pairing theories which are compared to the theory of s-wave pairing in the theory of superconducting electron.

CRC Press eBooks•Jul 22, 2026•0 citations

Algorithms for quantum molecular dynamics: from matrix product states to path integrals

Dmitri Iouchtchenko

This thesis describes several novel approaches in quantum molecular dynamics for obtaining properties of molecular systems in different regimes. We investigate ground state properties of chains of linear rotors with dipole–dipole interactions via the density matrix renormalization group (DMRG), by deriving the appropriate form of the interaction operator and implementing it in ITensor. This provides us with further evidence of a quantum phase transition in this system. We also improve the sampling of Gaussian mixture distributions for finite temperature path integral Monte Carlo (PIMC) of vibronic Hamiltonians. To do this, we replace random sampling by quasi-random sampling, and improve sampling distributions by optimizing their parameters. Finally, we introduce estimators and integrators for constrained free energy simulations in path integral molecular dynamics (PIMD). This method is applied to the study of a water dimer, for which we obtain a quantum potential of mean force.

•Jan 13, 2021•0 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•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

Fixed Point Quantum Monte Carlo

Romain Chessex, Massimo Borrelli, Hans Christian Öttinger

We present a new approach to the study of equilibrium properties in many-body quantum physics. Our method takes inspiration from Density Matrix Quantum Monte Carlo and incorporates new crucial features. First of all, the dynamics is transferred to the Laplace representation where an exact equation can be derived and solved using a simulation-step that, unlike most Monte Carlo methods, is not a priori physically bounded. Moreover, the spawning events are formulated in terms of two-process stochastic unravellings of quantum master equations, a formalism that is particularly useful when working with density matrices. And last, this is equivalent to an interaction picture, where the free part is integrated exactly and the convergence rate can be greatly increased if the interaction parameter is small. We benchmark our method by applying it to two case-studies in condensed matter physics, show its accuracy and further discuss its efficiency.

arXiv:2201.01383•Jan 4, 2022•0 citations
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.
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.
An introduction is given to the methodology of the path-integral formulation of quantum Monte Carlo simulation. This allows us to reformulate a d −dimensional partition function of a quantum mechanical system into the form of a d +1−dimensional path-integral written in classical mechanical variables. In addition, the quantum mechanical partition function is shown to be related by a process of analytical continuation to the time-development operator of the quantum mechanical system. The evaluation of the path-integral representation is discussed for the Maxwell-Boltzmann ideal gas, a chain of atoms, and liquid 4 He . The properties of liquid Helium I and Helium II are discussed in terms of the two-fluid model, the Ginzburg theory, and recent quantum Monte Carlo Simulation studies. The properties of liquid 3 He are discussed in the context of recent analytical p-wave pairing theories which are compared to the theory of s-wave pairing in the theory of superconducting electron.
This thesis describes several novel approaches in quantum molecular dynamics for obtaining properties of molecular systems in different regimes. We investigate ground state properties of chains of linear rotors with dipole–dipole interactions via the density matrix renormalization group (DMRG), by deriving the appropriate form of the interaction operator and implementing it in ITensor. This provides us with further evidence of a quantum phase transition in this system. We also improve the sampling of Gaussian mixture distributions for finite temperature path integral Monte Carlo (PIMC) of vibronic Hamiltonians. To do this, we replace random sampling by quasi-random sampling, and improve sampling distributions by optimizing their parameters. Finally, we introduce estimators and integrators for constrained free energy simulations in path integral molecular dynamics (PIMD). This method is applied to the study of a water dimer, for which we obtain a quantum potential of mean force.
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
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.
We present a new approach to the study of equilibrium properties in many-body quantum physics. Our method takes inspiration from Density Matrix Quantum Monte Carlo and incorporates new crucial features. First of all, the dynamics is transferred to the Laplace representation where an exact equation can be derived and solved using a simulation-step that, unlike most Monte Carlo methods, is not a priori physically bounded. Moreover, the spawning events are formulated in terms of two-process stochastic unravellings of quantum master equations, a formalism that is particularly useful when working with density matrices. And last, this is equivalent to an interaction picture, where the free part is integrated exactly and the convergence rate can be greatly increased if the interaction parameter is small. We benchmark our method by applying it to two case-studies in condensed matter physics, show its accuracy and further discuss its efficiency.