A nine-dimensional ab initio global potential energy surface for the H₂O⁺ + H₂ → H₃O⁺ + H reaction

An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81,000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm(-1). The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H4O(+) well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H2O(+) rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H2O(+) reactant is strongly coupled with the reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction.