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

Journal of Physics A Mathematical and General

Coherent-state path integral approach to the damped harmonic oscillator

Alexey Novikov, Ulrich Kleinekathöfer, Michael Schreiber

Published Feb 11, 2004• 4 citations
Publisher

Abstract

The time evolution of a system of bilinearly coupled bosonic modes is investigated using the real-time path integral technique in the coherent-state representation. 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 present treatment is applied to a dissipative harmonic oscillator within the Caldeira–Leggett model. 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 contrast to many perturbative approaches, the bath degrees of freedom are treated explicitly. 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 can the classical trajectories for the system coordinate be recovered.

DOI
10.1088/0305-4470/37/8/013

Related papers

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

Exact time evolution and master equations for the damped harmonic oscillator

Robert Karrlein, Hermann Grabert

Using the exact path integral solution for the damped harmonic oscillator it is shown that in general there does not exist an exact dissipative Liouville operator describing the dynamics of the oscillator for arbitrary initial bath preparations. Exact nonstationary Liouville operators can be found only for particular preparations. Three physically meaningful examples are examined. An exact master equation is derived for thermal initial conditions. Second, the Liouville operator governing the time evolution of equilibrium correlations is obtained. Third, factorizing initial conditions are studied. Additionally, one can show that there are approximate Liouville operators independent of the initial preparation describing the long-time dynamics under appropriate conditions. The general form of these approximate master equations is derived and the coefficients are determined for special cases of the bath spectral density including the Ohmic, Drude, and weak coupling cases. The connection with earlier work is discussed.

arXiv:physics/9610001•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1997•196 citations

Exact diagonalization of two quantum models for the damped harmonic oscillator

M. Rosenau da Costa, A. O. Caldeira, S. M. Dutra +1 more

The damped harmonic oscillator is a workhorse for the study of dissipation in quantum mechanics. However, despite its simplicity, this system has given rise to some approximations whose validity and relation to more refined descriptions deserve a thorough investigation. In this work, we apply a method that allows us to diagonalize exactly the dissipative Hamiltonians that are frequently adopted in the literature. Using this method, we derive the conditions of validity of the rotating-wave approximation (RWA) and show how this approximate description relates to more general ones. We also show that the existence of dissipative coherent states is intimately related to the RWA. Finally, through the evaluation of the dynamics of the damped oscillator, we notice an important property of the dissipative model that has not been properly accounted for in previous works, namely the necessity of new constraints to the application of the factorizable initial conditions.

arXiv:quant-ph/9903022•Physical Review A•Jan 10, 2000•72 citations

Time Evolution of Bath Properties in Spin-Boson Dynamics

Sohang Kundu, Nancy Makri

The dynamical behaviors of a two-level system (TLS) coupled to a harmonic dissipative bath has been studied extensively using a variety of analytical and numerical methods. The focus of the vast majority of these studies has been on the properties of the TLS, averaged with respect to the bath degrees of freedom. In this work, we use real-time path integral methods to probe the behavior of select bath degrees of freedom during the dynamics of a symmetric two-level system (TLS) coupled to a dissipative bath by calculating system-bath densities (SBD) and coordinate expectation values. Overall, the SBD motion on each diabatic state is simpler than the motion of the total density. In the weak coupling regime, which characterizes the parameters of oscillators that comprise such a bath, the SBD on each TLS state remains primarily compact and Gaussian-like, such that its peak is well characterized by the mode expectation value. In the absence of a dissipative environment, nonadiabatic density depletion leads to spikes in coordinate expectation values. The evolution of the SBD peak trajectory for two discrete modes exhibits Lissajous patterns with frequency-dependent shapes that strongly resemble classical trajectory motion on a torus. These patterns become more complex when the coupling of the mode to the TLS is increased outside of this regime, leading to persistent small amplitude oscillations in the TLS populations characterized by a very slow decay and SBD trajectories that exhibit behaviors reminiscent of chaotic classical systems. Indirect coupling to a dissipative bath has a stabilizing effect on the dynamics, eliminating spikes, synchronizing the SBD motion on the two diabatic states and regularizing the SBD trajectory to simple rectangular Lissajous-like shapes with a slowly shrinking boundary, regardless of the mode frequencies.

The Journal of Physical Chemistry B•Jul 20, 2021•11 citations

Dissipative dynamics of a biased qubit coupled to a harmonic oscillator: analytical results beyond the rotating wave approximation

Johannes Hausinger, Milena Grifoni

We study the dissipative dynamics of a biased two-level system (TLS) coupled to a harmonic oscillator (HO), the latter interacting with an Ohmic environment. Using the Van Vleck perturbation theory and going to second order in the coupling between the TLS and the HO, we show how the Hamiltonian of the TLS–HO system can be diagonalized analytically. Our model represents an improvement on the usually used Jaynes–Cummings Hamiltonian as an initial rotating wave approximation is avoided. By assuming a weak coupling to the thermal bath, analytical expressions for the time evolution of the populations of the TLS are found: the population is characterized by a multiplicity of damped oscillations together with a complex relaxation dynamics towards thermal equilibrium. The long time evolution is characterized by a single relaxation rate, which is largest at resonance and whose expression can be given in a closed analytic form.

arXiv:0806.3387•New Journal of Physics•Nov 20, 2008•46 citations

Dissipative dynamics of a biased qubit coupled to a harmonic oscillator: Analytical results beyond the rotating wave approximation

Johannes Hausinger, Milena Grifoni

We study the dissipative dynamics of a biased two-level system (TLS) coupled to a harmonic oscillator (HO), the latter interacting with an Ohmic environment. Using Van-Vleck perturbation theory and going to second order in the coupling between TLS and HO, we show how the Hamiltonian of the TLS-HO system can be diagonalized analytically. Our model represents an improvement to the usually used Jaynes-Cummings Hamiltonian as an initial rotating wave approximation is avoided. By assuming a weak coupling to the thermal bath, analytical expressions for the time evolution of the populations of the TLS are\nfound: the population is characterized by a multiplicity of damped oscillations together with a complex relaxation dynamics towards thermal equilibrium. The long time evolution is characterized by a single relaxation rate, which is largest at resonance and whose expression can be given in closed analytic form.

•Jan 1, 2008•43 citations
The time evolution of a system of bilinearly coupled bosonic modes is investigated using the real-time path integral technique in the coherent-state representation. 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 present treatment is applied to a dissipative harmonic oscillator within the Caldeira–Leggett model. 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 contrast to many perturbative approaches, the bath degrees of freedom are treated explicitly. 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 can the classical trajectories for the system coordinate be recovered.
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.
Using the exact path integral solution for the damped harmonic oscillator it is shown that in general there does not exist an exact dissipative Liouville operator describing the dynamics of the oscillator for arbitrary initial bath preparations. Exact nonstationary Liouville operators can be found only for particular preparations. Three physically meaningful examples are examined. An exact master equation is derived for thermal initial conditions. Second, the Liouville operator governing the time evolution of equilibrium correlations is obtained. Third, factorizing initial conditions are studied. Additionally, one can show that there are approximate Liouville operators independent of the initial preparation describing the long-time dynamics under appropriate conditions. The general form of these approximate master equations is derived and the coefficients are determined for special cases of the bath spectral density including the Ohmic, Drude, and weak coupling cases. The connection with earlier work is discussed.
The damped harmonic oscillator is a workhorse for the study of dissipation in quantum mechanics. However, despite its simplicity, this system has given rise to some approximations whose validity and relation to more refined descriptions deserve a thorough investigation. In this work, we apply a method that allows us to diagonalize exactly the dissipative Hamiltonians that are frequently adopted in the literature. Using this method, we derive the conditions of validity of the rotating-wave approximation (RWA) and show how this approximate description relates to more general ones. We also show that the existence of dissipative coherent states is intimately related to the RWA. Finally, through the evaluation of the dynamics of the damped oscillator, we notice an important property of the dissipative model that has not been properly accounted for in previous works, namely the necessity of new constraints to the application of the factorizable initial conditions.
The dynamical behaviors of a two-level system (TLS) coupled to a harmonic dissipative bath has been studied extensively using a variety of analytical and numerical methods. The focus of the vast majority of these studies has been on the properties of the TLS, averaged with respect to the bath degrees of freedom. In this work, we use real-time path integral methods to probe the behavior of select bath degrees of freedom during the dynamics of a symmetric two-level system (TLS) coupled to a dissipative bath by calculating system-bath densities (SBD) and coordinate expectation values. Overall, the SBD motion on each diabatic state is simpler than the motion of the total density. In the weak coupling regime, which characterizes the parameters of oscillators that comprise such a bath, the SBD on each TLS state remains primarily compact and Gaussian-like, such that its peak is well characterized by the mode expectation value. In the absence of a dissipative environment, nonadiabatic density depletion leads to spikes in coordinate expectation values. The evolution of the SBD peak trajectory for two discrete modes exhibits Lissajous patterns with frequency-dependent shapes that strongly resemble classical trajectory motion on a torus. These patterns become more complex when the coupling of the mode to the TLS is increased outside of this regime, leading to persistent small amplitude oscillations in the TLS populations characterized by a very slow decay and SBD trajectories that exhibit behaviors reminiscent of chaotic classical systems. Indirect coupling to a dissipative bath has a stabilizing effect on the dynamics, eliminating spikes, synchronizing the SBD motion on the two diabatic states and regularizing the SBD trajectory to simple rectangular Lissajous-like shapes with a slowly shrinking boundary, regardless of the mode frequencies.
We study the dissipative dynamics of a biased two-level system (TLS) coupled to a harmonic oscillator (HO), the latter interacting with an Ohmic environment. Using the Van Vleck perturbation theory and going to second order in the coupling between the TLS and the HO, we show how the Hamiltonian of the TLS–HO system can be diagonalized analytically. Our model represents an improvement on the usually used Jaynes–Cummings Hamiltonian as an initial rotating wave approximation is avoided. By assuming a weak coupling to the thermal bath, analytical expressions for the time evolution of the populations of the TLS are found: the population is characterized by a multiplicity of damped oscillations together with a complex relaxation dynamics towards thermal equilibrium. The long time evolution is characterized by a single relaxation rate, which is largest at resonance and whose expression can be given in a closed analytic form.
We study the dissipative dynamics of a biased two-level system (TLS) coupled to a harmonic oscillator (HO), the latter interacting with an Ohmic environment. Using Van-Vleck perturbation theory and going to second order in the coupling between TLS and HO, we show how the Hamiltonian of the TLS-HO system can be diagonalized analytically. Our model represents an improvement to the usually used Jaynes-Cummings Hamiltonian as an initial rotating wave approximation is avoided. By assuming a weak coupling to the thermal bath, analytical expressions for the time evolution of the populations of the TLS are\nfound: the population is characterized by a multiplicity of damped oscillations together with a complex relaxation dynamics towards thermal equilibrium. The long time evolution is characterized by a single relaxation rate, which is largest at resonance and whose expression can be given in closed analytic form.