A method for calculating the quantum canonical rate constant of chemical reactions in a many body system by means of a short-time flux autocorreltion function combined with a maximum entropy numerical analytic continuation scheme is presented. The rate constant is expressed as the time integral of the real-time flux autocorreltion function. The real-time flux autocorrelation function is evaluated for short times fully quantum mechanically by path integral Monte Carlo simulations. The maximum entropy approach is then used to extract the rate from the short real-time fux autocorreltion data. We present two numerical tests, one for proton transfer in harmoic dissipative environments in the photosynthetic reaction center. The results obtained using the flux autocorreltion data up to the time of no more than #beta#h are in excellent agreement with the eact quantum calculation over a wide range of parameters including even the tunneling regime.
展开▼
