The Journal of Chemical Physics

Influence functionals with semiclassical propagators in combined forward–backward time

Keiran C. Thompson, Nancy Makri

Published Jan 15, 1999• 170 citations

Abstract

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.

DOI: 10.1063/1.478011

Related papers

Rigorous forward-backward semiclassical formulation of many-body dynamics

Keiran C. Thompson, Nancy Makri

A rigorous, yet practical semiclassical formulation of time correlation functions or expectation values is presented. The main idea is to combine the forward and backward propagation steps into a single semiclassical propagator for those degrees of freedom that are not probed in a calculation, while retaining an explicit two-propagator description of the observable low-dimensional system. The combined forward-backward treatment of the environment naturally leads to extensive cancellation and results in an action that is small on the scale of Planck's constant. As a consequence, the integrand is smooth and thus amenable to Monte Carlo sampling. At the same time, important nonclassical effects terms, arising from interference among multiple bounce solutions of the system component, are fully accounted for.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topicsMay 1, 1999112 citations

Real time path integral methods for a system coupled to an anharmonic bath

Gregory Ilk, Nancy Makri

We propose an efficient path integral scheme for calculating the quantum dynamics of an arbitrary one-dimensional system coupled nonlinearly to many anharmonic noninteracting ‘‘bath’’ degrees of freedom. The starting point is an improved discretization of the path integral in terms of numerically constructed propagators [Chem. Phys. Lett. 193, 435 (1992)]. The resulting influence functional is comprised of one-dimensional correlation functions with step-structured time-dependent potentials and therefore is similar in structure to that employed in the spin-boson calculations of Coalson [J. Chem. Phys. 86, 995 (1987)]. In the present case, though, the influence functional is nonlinear and is computed using numerical iterative wave function propagation methods. Numerical tests on a system coupled to ten anharmonic oscillators demonstrate the efficiency of the proposed scheme, which requires numerical effort that scales only linearly with the number of anharmonic bath degrees of freedom.

The Journal of Chemical PhysicsOct 15, 199450 citations

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

Jiushu Shao, Nancy Makri

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 Journal of Chemical PhysicsJan 8, 200279 citations

Forward–backward semiclassical dynamics in the interaction representation

Jiushu Shao, Nancy Makri

The forward–backward semiclassical dynamics methodology [J. Phys. Chem. 103, 7753, 9479 (1999)] is reformulated in the interaction representation. The new version of the method allows for a fully quantum mechanical description of a low-dimensional subsystem of interest, along with a semiclassical forward–backward treatment of the solvent coordinates and their coupling to the reference subsystem. Application to the long-time tunneling dynamics in a symmetric double-well system coupled to a harmonic bath shows that the interaction FBSD is capable of capturing quantitatively the tunneling and decoherence effects induced by weakly dissipative environments.

The Journal of Chemical PhysicsSep 1, 200036 citations

Time correlation functions for mixed quantum-classical systems

J. Liam McWhirter

We consider the time correlation function of observables pertaining to a (quantum subsystem +bath), where the bath is coupled to a reservoir with many degrees of freedom. Integrating over the coordinates of this reservoir and assuming no initial correlations between the (quantum subsystem+bath) and the reservoir, we obtain an expression for the time correlation function that contains an influence functional. We then take the semiclassical and Fokker–Planck limits while modeling the reservoir with an Ohmic continuum of harmonic oscillators coupled bilinearily to the coordinates of the bath. The semiclassical limit is taken using a variant of Pechukas’ stationary phase analysis of the reduced propagator that yields a time correlation function written in terms of connected “classical” paths. These paths are got by solving the concatenation of several short-time interval Pechukas equations; as a result, the determination of these paths is more feasible than the determination of the “classical” path associated with a single long-time interval Pechukas equation. This concatenation includes the dissipative and stochastic forces associated with a classical Brownian particle. We then use decoherence arguments derived from an inspection of the influence functional to eliminate the phase interference structure of the bath. This elimination yields a mixed quantum-classical time correlation function that can be evaluated using nonadiabatic mixed quantum-classical dynamics schemes similar to those proposed recently by Webster and Tully.

The Journal of Chemical PhysicsApr 8, 19988 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, 20050 citations