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

Iterative path integral formulation of equilibrium correlation functions for quantum dissipative systems

Jiushu Shao, Nancy Makri

Published Jan 8, 2002• 79 citations
Publisher

Abstract

We present an iterative path integral algorithm for computing multitime correlation functions of a quantum system coupled to a dissipative bath of harmonic oscillators. By splitting the Boltzmann operator into two parts and reordering the propagators in the expression for canonical correlation functions, we are able to transform the evolution time contour into a symmetric one so that a forward propagation and a backward one are specified. Because the memory induced by the bath through the Feynman–Vernon influence functional decays rapidly in the complex time plane, long-time correlations are negligible. Taking advantage of this fact, we show that the correlation function can be obtained via an iterative procedure. The method is used to calculate three-time correlation functions of a dissipative two-level system.

DOI
10.1063/1.1423936

Related papers

Path integral renormalization for quantum dissipative dynamics with multiple timescales

Nancy Makri

The reduced density matrix of a system coupled to a harmonic dissipative bath can be obtained accurately using an iterative path integral methodology that involves a propagator tensor, whose rank is roughly equal to the number of time steps that span the bath-induced memory. The present work shows that even in situations where the bath is characterized by high-frequency modes (which necessitate a small time step) and slow oscillators (which give rise to very long memory), efficient propagation that takes into account the memory on all time scales is feasible by renormalizing the propagator matrix to include successively larger-length influence functional correlations. At its crudest level, the renormalized propagator is an ordinary matrix, thus the computational demands of the method are very modest, similar to those for propagating the bare system. Yet, test applications on dissipative two-level systems show that the method produces accurate results in demanding situations. By including memory on successively increasing time scales, the method is capable of capturing nonexponential kinetics induced by sluggish solvents.

Molecular Physics•May 10, 2012•30 citations

Iterative blip-summed path integral for quantum dynamics in strongly dissipative environments

Nancy Makri

The iterative decomposition of the blip-summed path integral [N. Makri, J. Chem. Phys. 141, 134117 (2014)] is described. The starting point is the expression of the reduced density matrix for a quantum system interacting with a harmonic dissipative bath in the form of a forward-backward path sum, where the effects of the bath enter through the Feynman-Vernon influence functional. The path sum is evaluated iteratively in time by propagating an array that stores blip configurations within the memory interval. Convergence with respect to the number of blips and the memory length yields numerically exact results which are free of statistical error. In situations of strongly dissipative, sluggish baths, the algorithm leads to a dramatic reduction of computational effort in comparison with iterative path integral methods that do not implement the blip decomposition. This gain in efficiency arises from (i) the rapid convergence of the blip series and (ii) circumventing the explicit enumeration of between-blip path segments, whose number grows exponentially with the memory length. Application to an asymmetric dissipative two-level system illustrates the rapid convergence of the algorithm even when the bath memory is extremely long.

The Journal of Chemical Physics•Apr 3, 2017•49 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

Iterative path integral calculations of non-linear spectra and entanglement dynamics

Mohammad Mehdi Sahrapour

We review an iterative path integral method that allows the computation of exact, long-time dynamics of small systems interacting with a dissipative bath. The method takes advantage of the finite memory of large environments at finite temperature to allow an iterative evaluation of the dynamics, thereby replacing an exponential scaling in simulation time with a linear scaling. This method is applied to calculate the dynamics of two model systems. The first consists of two-qubits interacting with a common bath. In this case we observe a variety of entanglement effects. If the qubits are initially separable, through an indirect coupling, the bath can create steady state entanglement between the qubits. This effect is due to the existence of a decoherence-free subspace as a result of the form of the qubit Hamiltonian and system-bath coupling. Entanglement created by the bath is shown to decrease with increasing temperature and system-bath coupling strength. However large system-bath coupling causes a faster increase in the entanglement. Initially entangled qubits lose their entanglement as a result of interactions with the bath, an effect that is heightened at higher temperatures. Direct coupling between the qubits is shown to slow the decay of entanglement and preserve some entanglement at long times; however at high temperatures this steady state entanglement becomes negligible. The second system we consider is vibrational degree of freedom coupled to a bath of harmonic oscillators or two-level systems. We compute four-time correlation functions which are used to calculate response functions relevant to third order infrared or seventh order Raman experiments for harmonic, Morse, and quadratic-quartic potentials. Our calculations reveal the role of potential features (anharmonicity and eigenvalue spectrum), both on short and long time scales, on the response function. Further, thermal excitation causes dramatic changes in the appearance of the response function, introducing symmetry with respect to the main diagonal. Finally, coupling to harmonic dissipative baths leads to decay of the response function (primarily along the t3 direction) and a broadening of the peaks in its Fourier transform. At high temperatures two-level-system baths are less efficient in destroying coherence than harmonic baths of similar parameters.

•May 24, 2013•0 citations

Influence functionals with semiclassical propagators in combined forward–backward time

Keiran C. Thompson, Nancy Makri

In a recent letter [Chem. Phys. Lett. 291, 101 (1998)] we presented a semiclassical methodology for calculating influence functionals arising from many-body anharmonic environments in the path integral formulation of quantum dynamics. Taking advantage of the trace operation associated with the unobservable medium, we express the influence functional in terms of a single propagator along a combined forward–backward system path. This propagator is evaluated according to time-dependent semiclassical theory in a coherent state initial value representation. Because the action associated with propagation in combined forward and backward time is governed by the net force experienced by the environment due to its interaction with the system, the resulting propagator is generally a smooth function of coordinates and thus amenable to Monte Carlo sampling; yet, the interference between forward and reverse propagators is fully accounted for. In the present paper we present a more elaborate version of the semiclassical influence functional formalism, along with algorithms for evaluating the coherent state transform of the Boltzmann operator that enters the influence functional. This factor is evaluated by performing an imaginary time path integral, and various approximations of the resulting expression as well as sampling schemes are discussed. The feasibility of the approach is demonstrated via numerous test calculations involving a two-level system coupled to (a) a dissipative harmonic bath and (b) a ten-dimensional bath of coupled anharmonic oscillators.

The Journal of Chemical Physics•Jan 15, 1999•170 citations

Iterative evaluation of the path integral for a system coupled to an anharmonic bath

Nancy Makri

An iterative algorithm is presented for evaluating the path integral expression for the reduced density matrix of a quantum system interacting with an anharmonic dissipative bath whose influence functional is obtained via numerical methods. The method allows calculation of the reduced density matrix over very long time periods.

The Journal of Chemical Physics•Oct 8, 1999•60 citations
We present an iterative path integral algorithm for computing multitime correlation functions of a quantum system coupled to a dissipative bath of harmonic oscillators. By splitting the Boltzmann operator into two parts and reordering the propagators in the expression for canonical correlation functions, we are able to transform the evolution time contour into a symmetric one so that a forward propagation and a backward one are specified. Because the memory induced by the bath through the Feynman–Vernon influence functional decays rapidly in the complex time plane, long-time correlations are negligible. Taking advantage of this fact, we show that the correlation function can be obtained via an iterative procedure. The method is used to calculate three-time correlation functions of a dissipative two-level system.
The reduced density matrix of a system coupled to a harmonic dissipative bath can be obtained accurately using an iterative path integral methodology that involves a propagator tensor, whose rank is roughly equal to the number of time steps that span the bath-induced memory. The present work shows that even in situations where the bath is characterized by high-frequency modes (which necessitate a small time step) and slow oscillators (which give rise to very long memory), efficient propagation that takes into account the memory on all time scales is feasible by renormalizing the propagator matrix to include successively larger-length influence functional correlations. At its crudest level, the renormalized propagator is an ordinary matrix, thus the computational demands of the method are very modest, similar to those for propagating the bare system. Yet, test applications on dissipative two-level systems show that the method produces accurate results in demanding situations. By including memory on successively increasing time scales, the method is capable of capturing nonexponential kinetics induced by sluggish solvents.
The iterative decomposition of the blip-summed path integral [N. Makri, J. Chem. Phys. 141, 134117 (2014)] is described. The starting point is the expression of the reduced density matrix for a quantum system interacting with a harmonic dissipative bath in the form of a forward-backward path sum, where the effects of the bath enter through the Feynman-Vernon influence functional. The path sum is evaluated iteratively in time by propagating an array that stores blip configurations within the memory interval. Convergence with respect to the number of blips and the memory length yields numerically exact results which are free of statistical error. In situations of strongly dissipative, sluggish baths, the algorithm leads to a dramatic reduction of computational effort in comparison with iterative path integral methods that do not implement the blip decomposition. This gain in efficiency arises from (i) the rapid convergence of the blip series and (ii) circumventing the explicit enumeration of between-blip path segments, whose number grows exponentially with the memory length. Application to an asymmetric dissipative two-level system illustrates the rapid convergence of the algorithm even when the bath memory is extremely long.
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.
We review an iterative path integral method that allows the computation of exact, long-time dynamics of small systems interacting with a dissipative bath. The method takes advantage of the finite memory of large environments at finite temperature to allow an iterative evaluation of the dynamics, thereby replacing an exponential scaling in simulation time with a linear scaling. This method is applied to calculate the dynamics of two model systems. The first consists of two-qubits interacting with a common bath. In this case we observe a variety of entanglement effects. If the qubits are initially separable, through an indirect coupling, the bath can create steady state entanglement between the qubits. This effect is due to the existence of a decoherence-free subspace as a result of the form of the qubit Hamiltonian and system-bath coupling. Entanglement created by the bath is shown to decrease with increasing temperature and system-bath coupling strength. However large system-bath coupling causes a faster increase in the entanglement. Initially entangled qubits lose their entanglement as a result of interactions with the bath, an effect that is heightened at higher temperatures. Direct coupling between the qubits is shown to slow the decay of entanglement and preserve some entanglement at long times; however at high temperatures this steady state entanglement becomes negligible. The second system we consider is vibrational degree of freedom coupled to a bath of harmonic oscillators or two-level systems. We compute four-time correlation functions which are used to calculate response functions relevant to third order infrared or seventh order Raman experiments for harmonic, Morse, and quadratic-quartic potentials. Our calculations reveal the role of potential features (anharmonicity and eigenvalue spectrum), both on short and long time scales, on the response function. Further, thermal excitation causes dramatic changes in the appearance of the response function, introducing symmetry with respect to the main diagonal. Finally, coupling to harmonic dissipative baths leads to decay of the response function (primarily along the t3 direction) and a broadening of the peaks in its Fourier transform. At high temperatures two-level-system baths are less efficient in destroying coherence than harmonic baths of similar parameters.
In a recent letter [Chem. Phys. Lett. 291, 101 (1998)] we presented a semiclassical methodology for calculating influence functionals arising from many-body anharmonic environments in the path integral formulation of quantum dynamics. Taking advantage of the trace operation associated with the unobservable medium, we express the influence functional in terms of a single propagator along a combined forward–backward system path. This propagator is evaluated according to time-dependent semiclassical theory in a coherent state initial value representation. Because the action associated with propagation in combined forward and backward time is governed by the net force experienced by the environment due to its interaction with the system, the resulting propagator is generally a smooth function of coordinates and thus amenable to Monte Carlo sampling; yet, the interference between forward and reverse propagators is fully accounted for. In the present paper we present a more elaborate version of the semiclassical influence functional formalism, along with algorithms for evaluating the coherent state transform of the Boltzmann operator that enters the influence functional. This factor is evaluated by performing an imaginary time path integral, and various approximations of the resulting expression as well as sampling schemes are discussed. The feasibility of the approach is demonstrated via numerous test calculations involving a two-level system coupled to (a) a dissipative harmonic bath and (b) a ten-dimensional bath of coupled anharmonic oscillators.