Energy diffusion controlled reaction rate in dissipative Hamiltonian systems

In this paper the energy diffusion controlled reaction rate in dissipative Hamiltonian systems is investigated by using the stochastic averaging method for quasi Hamiltonian systems. The boundary value problem of mean first- passage time （MFPT） of averaged system is formulated and the energy diffusion controlled reaction rate is obtained as the inverse of MFPT. The energy diffusion controlled reaction rate in the classical Kramers bistable potential and in a two-dimensional bistable potential with a heat bath are obtained by using the proposed approach respectively. The obtained results are then compared with those from Monte Carlo simulation of original systems and from the classical Kraraers theory. It is shown that the reaction rate obtained by using the proposed approach agrees well with that from Monte Carlo simulation and is more accurate than the classical Kramers rate.

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 multistate 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.