rate constant via flux-flux correlation function for linear potential energy surface

How does one calculate fully quantum mechanical rate ($\kappa$) in the golden-rule approximation for two linear potential energy surfaces?


Attempt:


Miller (83) proposes $\kappa=\int{Tr[\exp{(-\beta\hat{H})}\hat{F}\exp{(-i\hat{H}t/\hbar)}\hat{F}\exp{(i\hat{H}t/\hbar)}]}dt$


Where integrand is simply the flux-flux correlation function: $C\_{ff}(t)$. Which can be calculated under Feynman's path integral formalism. My attempt (which is in vain) at calculating $C\_{ff}(t)$ is as follows:


$C\_{ff}(t)=Tr[\exp{(-\beta\hat{H})}\hat{F}\exp{(-i\hat{H}t/\hbar)}\hat{F}\exp{(i\hat{H}t/\hbar)}]$


=$Tr[\exp{(-\beta\frac{\hat{H}}{2})}\hat{F}\exp{(-\beta\frac{\hat{H}}{2})}\exp{(-i\hat{H}t/\hbar)}\hat{F}\exp{(i\hat{H}t/\hbar)}]$


By cyclicly permuting the operators we reach at:


=$Tr[\exp{(i\hat{H}t/\hbar)}\exp{(-\beta\frac{\hat{H}}{2})}\hat{F}\exp{(-\beta\frac{\hat{H}}{2})}\exp{(-i\hat{H}t/\hbar)}\hat{F}]$


The boltzmann operator and quantum mechanical propagator can be combined as follows:


=$Tr[\exp{\hat{H}(\frac{it}{\hbar}-\frac{\beta}{2})}\hat{F}\exp{\hat{H}(\frac{-it}{\hbar}-\frac{\beta}{2})}\hat{F}]$


In the golden-rule (non-adiabatic) case, we have two electronic states 0 and 1. So F is simply a projection operator. Hence one can obtain:


=$Tr[\exp{\hat{H\_0}(\frac{it}{\hbar}-\frac{\beta}{2})}\exp{\hat{H\_1}(\frac{-it}{\hbar}-\frac{\beta}{2})}]$


This basically is kernel corresponding to two potential energy surfaces $V\_0$ and $V\_1$. For trajectory starting at $x\_a$ and ending at $x\_b$, we have


$C\_{ff}(t)=\int{\int{K\_0(x\_a,x\_b,\frac{it}{\hbar}-\frac{\beta}{2})K\_1(x\_b,x\_a,\frac{-it}{\hbar}-\frac{\beta}{2})}}dx\_adx\_b$ 


For a linear potential energy surfaces (PES), where my PES looks as follows:


$V\_0=k\_0 x$


$V\_1=k\_1 x$


My kernels are:


$K\_0=\sqrt{\frac{m}{2\pi t\_0}}\exp{(-S\_0)}$ 


$K\_1=\sqrt{\frac{m}{2\pi t\_1}}\exp{(-S\_1)}$


$S's$ correspond to action which is:


$S\_n(x\_a,x\_b,t\_n)=\frac{m(x\_a-x\_b)^2}{2 t\_n}-\frac{(x\_a+x\_b)k\_nt\_n}{2}-\frac{k\_n^2t\_n^3}{24m}$


The problem is the integral for flux flux correlation function doesn't seem to be converging with the imaginary argument for $t$'s. I am trying to integrate w.r.t $x\_a$, $x\_b$ and $t$ from -Inf to +Inf. My final answer for rate should look something like this:


$\exp{\frac{k\_0^2k\_1^2\hbar^2\beta^3}{24m(k\_0-k\_1)^2}}$


Is it a gaussian integral with respect to $x\_a$ and $x\_b$? One has to be careful because there is also an imaginary parts in the exponent. How does one reach the final answer for rate with those integrals? Really confused! Any help is appreciated.


