Real-time and imaginary-time quantum hierarchal Fokker-Planck equations

Abstract

We consider a quantum mechanical system represented in phase space (referred to hereafter as "Wigner space"), coupled to a harmonic oscillator bath. We derive quantum hierarchal Fokker-Planck (QHFP) equations not only in real time but also in imaginary time, which represents an inverse temperature. This is an extension of a previous work, in which we studied a spin-boson system, to a Brownian system. It is shown that the QHFP in real time obtained from a correlated thermal equilibrium state of the total system possesses the same form as those obtained from a factorized initial state. A modified terminator for the hierarchal equations of motion is introduced to treat the non-Markovian case more efficiently. Using the imaginary-time QHFP, numerous thermodynamic quantities, including the free energy, entropy, internal energy, heat capacity, and susceptibility, can be evaluated for any potential. These equations allow us to treat non-Markovian, non-perturbative system-bath interactions at finite temperature. Through numerical integration of the real-time QHFP for a harmonic system, we obtain the equilibrium distributions, the auto-correlation function, and the first-and second-order response functions. These results are compared with analytically exact results for the same quantities. This provides a critical test of the formalism for a non-factorized thermal state and elucidates the roles of fluctuation, dissipation, non-Markovian effects, and systembath coherence. Employing numerical solutions of the imaginary-time QHFP, we demonstrate the capability of this method to obtain thermodynamic quantities for any potential surface. It is shown that both types of QHFP equations can produce numerical results of any desired accuracy. The FORTRAN source codes that we developed, which allow for the treatment of Wigner space dynamics with any potential form (TanimuranFP15 and ImTanimuranFP15), are provided as the supplementary material. C 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License. [http://dx

Analysis of 2D THz-Raman spectroscopy using a non-Markovian Brownian oscillator model with nonlinear system-bath interactions

“…38 Although the contribution is minor, we plot theR RTT component in Appendix B. The relative intensities of each component evaluated in the harmonic and anharmonic SL BO cases are presented in Table I.…”

Generalized quantum Fokker-Planck equation for photoinduced nonequilibrium processes with positive definiteness condition

“…Thus, generalization of QFPE to include nonMarkovian and nonequilibrium effects remains an important and interesting theoretical issue to be addressed. Although more general hierarchical equations [18][19][20][21] may be used to this end, the benefit of having a single closed form equation, which can account for such non-Markovian and nonequilibrium effects, cannot be overestated.…”

Spins Dynamics in a Dissipative Environment: Hierarchal Equations of Motion Approach Using a Graphics Processing Unit (GPU)

“…We consider the case that the hierarchal elements are expressed as the M × M density matrix with the elements ρĵ 1 ,...,j K (n) (t) jk , where 1 ≤ j, k ≤ M. The HEOM are the simultaneous differential equations in time that expressed in terms of ADOs in the matrix forma as (12) where ρ(t) is the vector determined by…”