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

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

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

Published Dec 22, 2000• 17 citations
Publisher

Abstract

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.

DOI
10.1063/1.1327293

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

Invariant operators for quadratic Hamiltonians

Mark Andrews

The motion of wave packets can be easily determined for any Hamiltonian that is quadratic in position and momentum, even if the coefficients of the terms in the Hamiltonian vary with time. The method is based on the existence of an invariant operator, linear in position and momentum, the coefficients of these operators being solutions of the corresponding classical system. This immediately yields a set of very simple wave packets whose evolution is easily determined, and whose magnitude has the same form as the energy eigenfunctions of the harmonic oscillator (i.e., Gaussian or Hermite–Gaussian). This set provides a complete basis for finding the evolution of any state and the exact propagator is readily determined. For the harmonic oscillator, this set includes coherent states, squeezed states, displaced number states, and squeezed number states as special cases. The following important properties hold for every quadratic Hamiltonian, no matter what time dependence is present in the Hamiltonian: (1) The motion of the centroid of any wave packet separates from that of any moments relative to the centroid. (2) Every detail of the evolution of the quantum system can be calculated from the solutions of the corresponding classical system. (3) Wave functions of Gaussian (or Hermite–Gaussian) form will retain that form as they evolve.

American Journal of Physics•Apr 1, 1999•23 citations

A semiclassical approach to the dynamics of many-body Bose/Fermi systems by the path integral centroid molecular dynamics

Kenichi Kinugawa, Hidemi Nagao, Koji Ohta

We present a formalism of the path integral centroid molecular dynamics (CMD) extended to Bose and Fermi statistics as a semiclassical approach to explore the dynamics of quantum many-body systems. The validity of the method is examined in relation to the time correlation functions. The presently proposed scheme, refined from our previous derivation [Chem. Phys. Lett. 307, 187 (1999)], is aimed at the calculations of not the exact quantum-mechanical dynamics but the semiclassical dynamics under certain approximations. The formalism is based on the projection operator with which the Bose/Fermi system is mapped onto a particular type of pseudo-Boltzmann system. In the pseudo-Boltzmann system the correlation due to the Bose/Fermi statistics is introduced via an extra pseudopotential called the permutation potential and its relevant operator. Using the present semiclassical formalism, the time correlation function of centroid position, which is evaluated from the CMD trajectories in the pseudo-Boltzmann system, is an approximation to the Kubo canonical correlation function of position operator of the exact quantum-statistical system composed of bosons or fermions. There is no such apparent relation between the momentum operator and the corresponding momentum centroid.

The Journal of Chemical Physics•Jan 22, 2001•22 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

Path integral formulation of dissipative quantum dynamics

Alexey Novikov

In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.

•Jan 20, 2005•0 citations

Elementary derivation of the quantum propagator for the harmonic oscillator

Jiushu Shao

Operator algebra techniques are employed to derive the quantum evolution operator for the harmonic oscillator. The derivation begins with the construction of the annihilation and creation operators and the determination of the wave function for the coherent state as well as its time-dependent evolution, and ends with the transformation of the propagator in a mixed position-coherent-state representation to the desired one in configuration space. Throughout the entire procedure, besides elementary operator manipulations, it is only necessary to solve linear differential equations and to calculate Gaussian integrals.

American Journal of Physics•Sep 21, 2016•9 citations
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.
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 motion of wave packets can be easily determined for any Hamiltonian that is quadratic in position and momentum, even if the coefficients of the terms in the Hamiltonian vary with time. The method is based on the existence of an invariant operator, linear in position and momentum, the coefficients of these operators being solutions of the corresponding classical system. This immediately yields a set of very simple wave packets whose evolution is easily determined, and whose magnitude has the same form as the energy eigenfunctions of the harmonic oscillator (i.e., Gaussian or Hermite–Gaussian). This set provides a complete basis for finding the evolution of any state and the exact propagator is readily determined. For the harmonic oscillator, this set includes coherent states, squeezed states, displaced number states, and squeezed number states as special cases. The following important properties hold for every quadratic Hamiltonian, no matter what time dependence is present in the Hamiltonian: (1) The motion of the centroid of any wave packet separates from that of any moments relative to the centroid. (2) Every detail of the evolution of the quantum system can be calculated from the solutions of the corresponding classical system. (3) Wave functions of Gaussian (or Hermite–Gaussian) form will retain that form as they evolve.
We present a formalism of the path integral centroid molecular dynamics (CMD) extended to Bose and Fermi statistics as a semiclassical approach to explore the dynamics of quantum many-body systems. The validity of the method is examined in relation to the time correlation functions. The presently proposed scheme, refined from our previous derivation [Chem. Phys. Lett. 307, 187 (1999)], is aimed at the calculations of not the exact quantum-mechanical dynamics but the semiclassical dynamics under certain approximations. The formalism is based on the projection operator with which the Bose/Fermi system is mapped onto a particular type of pseudo-Boltzmann system. In the pseudo-Boltzmann system the correlation due to the Bose/Fermi statistics is introduced via an extra pseudopotential called the permutation potential and its relevant operator. Using the present semiclassical formalism, the time correlation function of centroid position, which is evaluated from the CMD trajectories in the pseudo-Boltzmann system, is an approximation to the Kubo canonical correlation function of position operator of the exact quantum-statistical system composed of bosons or fermions. There is no such apparent relation between the momentum operator and the corresponding momentum centroid.
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.
In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.
Operator algebra techniques are employed to derive the quantum evolution operator for the harmonic oscillator. The derivation begins with the construction of the annihilation and creation operators and the determination of the wave function for the coherent state as well as its time-dependent evolution, and ends with the transformation of the propagator in a mixed position-coherent-state representation to the desired one in configuration space. Throughout the entire procedure, besides elementary operator manipulations, it is only necessary to solve linear differential equations and to calculate Gaussian integrals.