Dynamics of the isotope exchange reaction of D with H3 +, H2D+, and D2H. J Chem Phys 2021;154:084307. [PMID: 33639774 DOI: 10.1063/5.0038434]

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

Motivated by a recent semiclassical analysis of chemical reaction thresholds [Bonnet et al., J. Chem. Phys., 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 et al., Phys. Chem. Chem. Phys., 2022, 24, 5489]. The rate constants computed with the new, general algorithm closely match recent Ring Polymer Molecular Dynamics (RPMD) [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] and experimentally derived [Bowen et al., J. Chem. Phys., 2021, 154, 084307] ones spanning ∼4 orders of magnitude from 70 to 1500 K.

Low-temperature reaction dynamics of paramagnetic species in the gas phase

Radicals are abundant in a range of important gas-phase environments. They are prevalent in the atmosphere, in interstellar space, and in combustion processes. As such, understanding how radicals react is essential for the development of accurate models of the complex chemistry occurring in these gas-phase environments. By controlling the properties of the colliding reactants, we can also gain insights into how radical reactions occur on a fundamental level. Recent years have seen remarkable advances in the breadth of experimental methods successfully applied to the study of reaction dynamics involving paramagnetic species-from improvements to the well-known crossed molecular beams approach to newer techniques involving magnetically guided and decelerated beams. Coupled with ever-improving theoretical methods, quantum features are being observed and interesting insights into reaction dynamics are being uncovered in an increasingly diverse range of systems. In this highlight article, we explore some of the exciting recent developments in the study of chemical dynamics involving paramagnetic species. We focus on low-energy reactive collisions involving neutral radical species, where the reaction parameters are controlled. We conclude by identifying some of the limitations of current methods and exploring possible new directions for the field.

Capturing quantum effects with quasi-classical trajectories in the D + H+3 → H2D+ + H reaction

We present quasi-classical trajectory (QCT) cross sections, rate constants, and product state distributions for the D + H+3 → H₂D⁺ + H reaction. Using the same H+4 potential surface, the rate constants obtained from several QCT-based methods correcting for zero-point effects by Gaussian binning the product H₂D⁺ are compared to ring polymer molecular dynamics (RPMD) rate constants [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] which include quantum effects and to recent experimentally derived rate constants [Bowen et al., J. Chem. Phys., 2021, 154, 084307]. QCT with standard binning predicts rate constants that increase slowly as the temperature decreases from 1500 to 100 K. In contrast, the RPMD rate constants decrease rapidly with decreasing temperature. By 100 K, the QCT standard binning rate constant is more than 3 orders of magnitude larger than the RPMD rate constant. We show that QCT with Gaussian binning and proper normalization captures the zero-point effects and reproduces the RPMD rate constants over a large temperature range. Furthermore, the simple technique of counting only reactive trajectories with vibrational energy above the product zero-point energy matches the RPMD results well down to ∼300 K. The present Gaussian binned rate constants are in fair agreement with new experimentally derived rate constants from 100 to 1500 K. However, because the Gaussian binned rate constants do not include tunneling, important at lower temperatures, and the RPMD and experimentally derived rate constants have significant differences, the roles of the competing effects of zero-point energy, internal excitation of the H+3, and quantum tunneling are not simple and require further study for a consistent picture of the dynamics. Since rate constants for complex forming reactions, such as the title reaction, are difficult to converge with RPMD, alternative QCT-based methods, which include quantum effects and in addition provide product state distributions as described here, are highly desirable.