Quantum mechanical canonical rate theory: A new approach based on the reactive flux and numerical analytic continuation methods

We present the reactive flux analytic continuation (RFAC) method, based on the quantum reactive flux formalism combined with a numerical analytic cont inuation approach to calculate quantum canonical rates in condensed phase s ystems. We express the imaginary time reactive-flux correlation function in terms of a frequency dependent rate constant, and use path integral formal ism to derive a working expression suitable for Monte Carlo simulation tech niques. The imaginary time data obtained by simulation is analytically cont inued to the real time using the maximum entropy method to obtain the react ion rate. Motivated by the success of the method to predict the rates for a simple one dimensional parabolic barrier model, we assess its accuracy for a condensed phase reaction modeled by a double-well coupled to a harmonic bath. We note that the method is applicable to a more general Hamiltonian a s long as the reaction coordinate can be identified. The reaction rates com puted in this fashion are in very good agreement with analytic and numerica lly exact results. We demonstrate the applicability of the method for a wid e range of model parameters and temperatures. (C) 2000 American Institute o f Physics. [S0021-9606(00)50606-0].