Mixed Quantum Classical Reaction Rates based on the Phase Space Formulation of the Hierarchical Equations of Motion

In a previous work J. Chem. Phys. 140 , 174105 (2014)], we have shown that a mixed quantum classical (MQC) rate theory can be derived to investigate the quantum tunneling effects in the proton transfer reactions. However, the method is based on the high temperature approximation of the hierarchical equation of motion (HEOM) with the Debye-Drude spectral density, and results in a multi-state Zusman type of equation. We now extend this theory to include quantum effects of the bath degrees of freedom. By writing the full HEOM into a multidimensional partial differential equation in phase space, we can define a new reaction coordinate, and the previous method can be generalized to the full quantum regime. The validity of the new method is demonstrated by using numerical examples, including the spin-Boson model, and the double well model for proton transfer reaction. The new method is found to resolve some key problems of the previous theory based on high temperature approximation, including possible numerical instability in long time simulation and wrong rate constant at low temperatures.