An efficient algorithm for capturing quantum effects in classical reactive scattering: application to D + H+3 → H 2 D + + H.

Motivated by a recent semiclassical analysis of chemical reaction thresholds [Bonnet , , 2022, 157, 094114], we present an efficient algorithm for including zero-point energy (ZPE) effects in classical reactive scattering. The algorithm is an extension of the quasi-classical trajectory (QCT) Gaussian binning method. We apply it to the astrophysically important D + H+3 reaction, where there are significant quantum effects and where application of other methods is problematic [Braunstein , , 2022, , 5489]. The rate constants computed with the new, general algorithm closely match recent Ring Polymer Molecular Dynamics (RPMD) [Bulut , , 2019, , 8766] and experimentally derived [Bowen , , 2021, , 084307] ones spanning ∼4 orders of magnitude from 70 to 1500 K.