Filtered propagator functional for iterative dynamics of quantum dissipative systems

Abstract We describe a Fortran program which calculates the reduced density matrix of a one-dimensional quantum mechanical continuous or discrete system coupled to a harmonic dissipative environment. The algorithm is based on Feynman's path integral formulation of time-dependent quantum mechanics. An adiabatic reference is employed to obtain accurate propagators and the harmonic bath is replaced by an influence functional which is discretized by optimal discrete variable representations. A propagator functional of statistically significant path segments is constructed which allows iterative evaluation of the path integral over long time periods. High efficiency is achieved with the aid of sorting and filtering criteria. The appended program is executable in either serial or parallel mode.

[1]  R. Zwanzig.  Nonlinear generalized Langevin equations , 1973 .

[2]  J. Doll,et al.  Equilibrium and dynamical Fourier path integral methods , 2007 .

[3]  R. Feynman,et al.  Quantum Mechanics and Path Integrals , 1965 .

[4]  John C. Light,et al.  Theoretical Methods for Rovibrational States of Floppy Molecules , 1989 .

[5]  William H. Press,et al.  Numerical recipes , 1990 .

[6]  Robert E. Wyatt,et al.  Optical potential for laser induced dissociation , 1983 .

[7]  Nancy Makri,et al.  System-specific discrete variable representations for path integral calculations with quasi-adiabatic propagators , 1993 .

[8]  N. Makri.  Improved Feynman propagators on a grid and non-adiabatic corrections within the path integral framework , 1992 .

[9]  R. Feynman,et al.  Space-Time Approach to Non-Relativistic Quantum Mechanics , 1948 .

[10]  TENSOR PROPAGATOR WITH WEIGHT-SELECTED PATHS FOR QUANTUM DISSIPATIVE DYNAMICS WITH LONG-MEMORY KERNELS , 1996 .

[11]  W. Domcke,et al.  Recursive evaluation of the real-time path integral for dissipative systems. The spin-boson model , 1995 .

[12]  A. Leggett,et al.  Dynamics of the dissipative two-state system , 1987 .

[13]  N. Makri.  Feynman path integration in quantum dynamics , 1991 .

[14]  Ronnie Kosloff,et al.  Solution of the time-dependent Liouville-von Neumann equation: dissipative evolution , 1992 .

[15]  N. Makri,et al.  Tensor propagator for iterative quantum time evolution of reduced density matrices. II. Numerical methodology , 1995 .

[16]  B. Shore,et al.  Modelling laser ionisation , 1978 .

[17]  Nancy Makri,et al.  Path integrals for dissipative systems by tensor multiplication. Condensed phase quantum dynamics for arbitrarily long time , 1994 .

[18]  N. Makri,et al.  TENSOR PROPAGATOR FOR ITERATIVE QUANTUM TIME EVOLUTION OF REDUCED DENSITY MATRICES. I: THEORY , 1995 .

[19]  N. Makri,et al.  Time-dependent discrete variable representations for quantum wave packet propagation , 1995 .