Ionic collision processes in simple gas mixtures containing hydrogen

Ring-Polymer Molecular Dynamics Calculations of Thermal Rate Coefficients and Branching Ratios for the Interstellar H3+ + CO → H2 + HCO+/HOC+ Reaction and Its Deuterated Analogue

The reaction between H₃⁺ and CO is important in understanding the H₃⁺ destruction mechanism in the interstellar medium. In this work, thermal rate coefficients for the H₃⁺ + CO and D₃⁺ + CO reactions are calculated using ring-polymer molecular dynamics (RPMD) on a high-level machine-learning potential energy surface. The RPMD results agree well with the classical molecular dynamics results, where nuclear quantum effects are completely ignored, whereas the agreement between the RPMD results and the previous quasi-classical trajectory is good only at low temperatures. The calculated [HCO⁺]/[HOC⁺] product branching ratios decrease as the temperature increases, and the product branching is exclusively determined by the initial collisional orientation, which governs the formation of an ion-dipole complex, H₃⁺···CO or H₃⁺···OC, that dissociates into H₂ + HCO⁺ or H₂ + HOC⁺, respectively, via a direct mechanism. However, the contribution of the indirect mechanism via the rearrangement between H₃⁺···CO and H₃⁺···OC increases as the temperature increases, although its absolute fraction is small.