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

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

Jianshu Cao, Gregory A. Voth

Published Apr 1, 1994• 364 citations
Publisher

Abstract

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.

DOI
10.1063/1.467175

Related papers

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

Seogjoo Jang, Gregory A. Voth

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.

The Journal of Chemical Physics•Aug 8, 1999•216 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

On the Feynman path centroid density for Bose-Einstein and Fermi-Dirac statistics

Pierre–Nicholas Roy, Gregory A. Voth

Expressions are derived for the Feynman path centroid density of systems of identical particles, namely, Bosons and Fermions. The resulting expressions are applied to a system of two independent particles in a quartic oscillator and the corresponding centroid densities are computed, plotted and compared to the distinguishable particle case. In regions where the particle centroids are close to each other, the Boson centroid density displays an enhanced amplitude in comparison to the distinguishable case. This behavior is attributed to the attractive correlations in Boson systems. The Fermion centroid density, however, can have the peculiar property of being negative in those regions where the particle centroids are in close proximity. This feature is related to the Pauli exclusion principle. This property of the Fermion centroid density rules out its strict interpretation as a probability distribution. Equilibrium properties such as the canonical partition function and the average position were accurately computed using both densities. The Boson density was also used to compute the position autocorrelation function using the Centroid Molecular Dynamics method and the results are in excellent agreement with those of an exact quantum calculation.

The Journal of Chemical Physics•Feb 22, 1999•52 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

A reductionist perspective on quantum statistical mechanics: Coarse-graining of path integrals

Anton V. Sinitskiy, Gregory A. Voth

Computational modeling of the condensed phase based on classical statistical mechanics has been rapidly developing over the last few decades and has yielded important information on various systems containing up to millions of atoms. However, if a system of interest contains important quantum effects, well-developed classical techniques cannot be used. One way of treating finite temperature quantum systems at equilibrium has been based on Feynman's imaginary time path integral approach and the ensuing quantum-classical isomorphism. This isomorphism is exact only in the limit of infinitely many classical quasiparticles representing each physical quantum particle. In this work, we present a reductionist perspective on this problem based on the emerging methodology of coarse-graining. This perspective allows for the representations of one quantum particle with only two classical-like quasiparticles and their conjugate momenta. One of these coupled quasiparticles is the centroid particle of the quantum path integral quasiparticle distribution. Only this quasiparticle feels the potential energy function. The other quasiparticle directly provides the observable averages of quantum mechanical operators. The theory offers a simplified perspective on quantum statistical mechanics, revealing its most reductionist connection to classical statistical physics. By doing so, it can facilitate a simpler representation of certain quantum effects in complex molecular environments.

The Journal of Chemical Physics•Sep 2, 2015•16 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
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.
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.
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.
Expressions are derived for the Feynman path centroid density of systems of identical particles, namely, Bosons and Fermions. The resulting expressions are applied to a system of two independent particles in a quartic oscillator and the corresponding centroid densities are computed, plotted and compared to the distinguishable particle case. In regions where the particle centroids are close to each other, the Boson centroid density displays an enhanced amplitude in comparison to the distinguishable case. This behavior is attributed to the attractive correlations in Boson systems. The Fermion centroid density, however, can have the peculiar property of being negative in those regions where the particle centroids are in close proximity. This feature is related to the Pauli exclusion principle. This property of the Fermion centroid density rules out its strict interpretation as a probability distribution. Equilibrium properties such as the canonical partition function and the average position were accurately computed using both densities. The Boson density was also used to compute the position autocorrelation function using the Centroid Molecular Dynamics method and the results are in excellent agreement with those of an exact quantum calculation.
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.
Computational modeling of the condensed phase based on classical statistical mechanics has been rapidly developing over the last few decades and has yielded important information on various systems containing up to millions of atoms. However, if a system of interest contains important quantum effects, well-developed classical techniques cannot be used. One way of treating finite temperature quantum systems at equilibrium has been based on Feynman's imaginary time path integral approach and the ensuing quantum-classical isomorphism. This isomorphism is exact only in the limit of infinitely many classical quasiparticles representing each physical quantum particle. In this work, we present a reductionist perspective on this problem based on the emerging methodology of coarse-graining. This perspective allows for the representations of one quantum particle with only two classical-like quasiparticles and their conjugate momenta. One of these coupled quasiparticles is the centroid particle of the quantum path integral quasiparticle distribution. Only this quasiparticle feels the potential energy function. The other quasiparticle directly provides the observable averages of quantum mechanical operators. The theory offers a simplified perspective on quantum statistical mechanics, revealing its most reductionist connection to classical statistical physics. By doing so, it can facilitate a simpler representation of certain quantum effects in complex molecular environments.
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.