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

Path integral centroid variables and the formulation of their exact real time dynamics

Seogjoo Jang, Gregory A. Voth

Published Aug 8, 1999• 216 citations
Publisher

Abstract

A formalism is presented in this paper which, for the first time, establishes the theoretical basis for the quantum time evolution of path integral centroid variables and also provides clear motivation for using these variables to study condensed phase quantum dynamics. The equilibrium centroid distribution is first shown to be a well-defined distribution function which is specific to the canonical density operator. A quantum mechanical quasi-density operator (QDO) is associated with each value of the distribution so that, upon application of the standard quantum mechanical formalism, the QDO can be used to provide a rigorous definition of both static and dynamical centroid variables. Various properties of the dynamical centroid variables are derived, including the perspective that the centroid constraint on the imaginary time paths introduces a nonstationarity in the equilibrium ensemble which, in turn, can be shown to yield information on the correlations of spontaneous fluctuations. The analytic solution for the harmonic oscillator and a numerical solution for a double well system are provided which illustrate the various aspects of the theory. The theory contained herein provides the basis for a derivation of Centroid Molecular Dynamics, as well as the systematic improvements of that theory.

DOI
10.1063/1.479514

Related papers

The formulation of quantum statistical mechanics based on the Feynman path centroid density. I. Equilibrium properties

Jianshu Cao, Gregory A. Voth

A formulation of quantum statistical mechanics is discussed in which the Feynman path centroid density in Feynman path integration is recast as the central statistical distribution used to average equilibrium and dynamical quantities. In this formulation, the path integral centroid density occupies the same role as the Boltzmann density in classical statistical mechanics. Therefore, the statistical ensemble of imaginary time path centroid configurations provides the distribution which is used to average the appropriately formulated effective operators and imaginary time correlation functions. An accurate renormalized diagrammatic perturbation theory for the centroid density and centroid-constrained imaginary time propagator will also be described with particular emphasis given to the mathematical advantages arising from the centroid-based formulation. The present paper is concerned with the calculation of equilibrium properties from the centroid perspective, while the companion paper describes a centroid-based formalism for calculating dynamical time correlation functions.

The Journal of Chemical Physics•Apr 1, 1994•364 citations

On the Feynman path centroid density as a phase space distribution in quantum statistical mechanics

Rigoberto Hernandez, Jianshu Cao, Gregory A. Voth

The phase space formulation of quantum statistical mechanics using the Feynman path centroid density offers an alternative perspective to the standard Wigner prescription for the classical-like evaluation of equilibrium and/or dynamical quantities of statistical systems. The use of this formulation has been implicit in recent work on quantum rate theories, for example, in which the centroid density distribution replaces the classical Boltzmann distribution. In order to further understand the approximations involved in this and similar transcriptions, the present work elaborates and clarifies the issue of operator ordering in a rigorous centroid-based formulation. In particular, through the use of the Weyl correspondence, a precise definition of the centroid symbol of operators and their products is presented. Though we fall short of finding the algebraic structure tantamount to that found in the Weyl symbols—of which the Wigner distribution is an example— the resulting expressions have internal consistency and are amenable to approximate evaluation through cumulant expansions.

The Journal of Chemical Physics•Sep 22, 1995•75 citations

Spectral decomposition and Bloch equation of the operators represented by fixed-centroid path integrals

Telesforo López-Ciudad, Rafael Ramı́rez

Interesting approaches to study statical and dynamic properties of quantum systems, e.g., the quantum transition state theory and the centroid molecular dynamics, have been previously derived using fixed centroid path integrals. We show that these constrained propagators can be alternatively defined using an operator formalism. An interesting result is the finding of the differential equations that determine the temperature dependence of these propagators. One equation applies to path integrals with fixed-centroid position (i.e., those used in quantum transition state theory), and the other one to path integrals with fixed-centroid position and momentum (i.e., those used in centroid molecular dynamics). Both equations are solved for a harmonic oscillator, so that the spectral decomposition of the operators represented by fixed-centroid path integrals is derived for this particular case. Their eigenvalues build an alternating geometric series, showing explicitly the impossibility of considering such operators as true density operators, i.e., some eigenfunctions are associated to “negative probabilities.” The eigenfunctions are shown to be a generalization of the coherent and squeezed states of the harmonic oscillator. The physical meaning of centroid molecular dynamics, an approximation to study the time evolution of these mixed states, is clarified by considering the time evolution of the corresponding eigenfunctions. The mixed states constructed with “negative probabilities” display vanishing small position and momentum dispersion in the high temperature limit.

The Journal of Chemical Physics•Dec 22, 2000•17 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

Operator formulation of centroid dynamics for Bose–Einstein and Fermi–Dirac statistics

Nicholas V. Blinov, Pierre–Nicholas Roy

This paper is devoted to the development of an operator formulation of the recent extension of the centroid molecular dynamics method [J. Chem. Phys. 110, 3647 (1999); 111, 5303 (1999)] to boson and fermion systems. An operator calculus is used to rederive the basic equations of centroid dynamics. The following generalization to the case of systems of many indistinguishable particles is based on the use of a projection operator. Two different definitions of the quasi-density operator for bosonic and fermionic systems are suggested. The first definition allows an exact evaluation of equilibrium properties for systems with exchange effects using classical-like molecular dynamics calculations. The second one provides a formal justification of Bose–Einstein/Fermi–Dirac centroid dynamics with the same set of approximations as for Boltzmann statistics, and can be used to extract quantum dynamical information. In this case, the corresponding centroid correlation function can be related to a double Kubo transformed quantum mechanical one.

The Journal of Chemical Physics•Nov 1, 2001•28 citations

Open quantum systems beyond equilibrium: Lindblad equation and path integral molecular dynamics

Benedikt M. Reible, Somayeh Ahmadkhani, Luigi Delle Site

The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.

quant-pharXiv:2603.10839•Physical Review A•Mar 12, 2026•0 citations
A formalism is presented in this paper which, for the first time, establishes the theoretical basis for the quantum time evolution of path integral centroid variables and also provides clear motivation for using these variables to study condensed phase quantum dynamics. The equilibrium centroid distribution is first shown to be a well-defined distribution function which is specific to the canonical density operator. A quantum mechanical quasi-density operator (QDO) is associated with each value of the distribution so that, upon application of the standard quantum mechanical formalism, the QDO can be used to provide a rigorous definition of both static and dynamical centroid variables. Various properties of the dynamical centroid variables are derived, including the perspective that the centroid constraint on the imaginary time paths introduces a nonstationarity in the equilibrium ensemble which, in turn, can be shown to yield information on the correlations of spontaneous fluctuations. The analytic solution for the harmonic oscillator and a numerical solution for a double well system are provided which illustrate the various aspects of the theory. The theory contained herein provides the basis for a derivation of Centroid Molecular Dynamics, as well as the systematic improvements of that theory.
A formulation of quantum statistical mechanics is discussed in which the Feynman path centroid density in Feynman path integration is recast as the central statistical distribution used to average equilibrium and dynamical quantities. In this formulation, the path integral centroid density occupies the same role as the Boltzmann density in classical statistical mechanics. Therefore, the statistical ensemble of imaginary time path centroid configurations provides the distribution which is used to average the appropriately formulated effective operators and imaginary time correlation functions. An accurate renormalized diagrammatic perturbation theory for the centroid density and centroid-constrained imaginary time propagator will also be described with particular emphasis given to the mathematical advantages arising from the centroid-based formulation. The present paper is concerned with the calculation of equilibrium properties from the centroid perspective, while the companion paper describes a centroid-based formalism for calculating dynamical time correlation functions.
The phase space formulation of quantum statistical mechanics using the Feynman path centroid density offers an alternative perspective to the standard Wigner prescription for the classical-like evaluation of equilibrium and/or dynamical quantities of statistical systems. The use of this formulation has been implicit in recent work on quantum rate theories, for example, in which the centroid density distribution replaces the classical Boltzmann distribution. In order to further understand the approximations involved in this and similar transcriptions, the present work elaborates and clarifies the issue of operator ordering in a rigorous centroid-based formulation. In particular, through the use of the Weyl correspondence, a precise definition of the centroid symbol of operators and their products is presented. Though we fall short of finding the algebraic structure tantamount to that found in the Weyl symbols—of which the Wigner distribution is an example— the resulting expressions have internal consistency and are amenable to approximate evaluation through cumulant expansions.
Interesting approaches to study statical and dynamic properties of quantum systems, e.g., the quantum transition state theory and the centroid molecular dynamics, have been previously derived using fixed centroid path integrals. We show that these constrained propagators can be alternatively defined using an operator formalism. An interesting result is the finding of the differential equations that determine the temperature dependence of these propagators. One equation applies to path integrals with fixed-centroid position (i.e., those used in quantum transition state theory), and the other one to path integrals with fixed-centroid position and momentum (i.e., those used in centroid molecular dynamics). Both equations are solved for a harmonic oscillator, so that the spectral decomposition of the operators represented by fixed-centroid path integrals is derived for this particular case. Their eigenvalues build an alternating geometric series, showing explicitly the impossibility of considering such operators as true density operators, i.e., some eigenfunctions are associated to “negative probabilities.” The eigenfunctions are shown to be a generalization of the coherent and squeezed states of the harmonic oscillator. The physical meaning of centroid molecular dynamics, an approximation to study the time evolution of these mixed states, is clarified by considering the time evolution of the corresponding eigenfunctions. The mixed states constructed with “negative probabilities” display vanishing small position and momentum dispersion in the high temperature limit.
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.
This paper is devoted to the development of an operator formulation of the recent extension of the centroid molecular dynamics method [J. Chem. Phys. 110, 3647 (1999); 111, 5303 (1999)] to boson and fermion systems. An operator calculus is used to rederive the basic equations of centroid dynamics. The following generalization to the case of systems of many indistinguishable particles is based on the use of a projection operator. Two different definitions of the quasi-density operator for bosonic and fermionic systems are suggested. The first definition allows an exact evaluation of equilibrium properties for systems with exchange effects using classical-like molecular dynamics calculations. The second one provides a formal justification of Bose–Einstein/Fermi–Dirac centroid dynamics with the same set of approximations as for Boltzmann statistics, and can be used to extract quantum dynamical information. In this case, the corresponding centroid correlation function can be related to a double Kubo transformed quantum mechanical one.
The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.