Path-integral calculations of heavy atom kinetic isotope effects in condensed phase reactions using higher-order trotter factorizations

A convenient approach to compute kinetic isotope effects (KIEs) in condensed phase chemical reactions is via path integrals (PIs). Usually, the primitive approximation is used in PI simulations, although such quantum simulations are computationally demanding. The efficiency of PI simulations may be greatly improved, if higher-order Trotter factorizations of the density matrix operator are used. In this study, we use a higher-order PI method, in conjunction with mass-perturbation, to compute heavy-atom KIE in the decarboxylation of orotic acid in explicit sulfolane solvent. The results are in good agreement with experiment and show that the mass-perturbation higher-order Trotter factorization provides a practical approach for computing condensed phase heavy-atom KIE.