Quantum-classical Liouville dynamics of proton and deuteron transfer rates in a solvated hydrogen-bonded complex

Proton and deuteron transfer reactions in a hydrogen-bonded complex dissolved in a polar solution are studied using quantum-classical Liouville dynamics. Reactive-flux correlation functions that involve quantum-classical Liouville dynamics for species operators and quantum equilibrium sampling are used to calculate the rate constants. Adiabatic and nonadiabatic reaction rates are computed, compared, and analyzed. Large variations of the kinetic isotope effect (KIE) for this reaction have been observed in the literature, which depend on the nature of the approximate calculation used to estimate the proton and deuteron transfer rates. Our estimate of the KIE lies at the low end of the range of previously observed values, suggesting a rather small KIE for this reaction.