Reaction Dynamics of the H + HeH⁺ → He + H₂ ⁺ System

The reaction dynamics for the H + HeH⁺ → He + H₂ ⁺ reaction in its electronic ground state is investigated using two different representations of the potential energy surface (PES). The first uses a combined kernel and neural network representation of UCCSD-(T) reference data whereas the second is a corrected PES (cR-PES) that eliminates an artificial barrier in the entrance channel appearing in its initial expansion based on full configuration interaction reference data. Despite the differences between the two PESs, both yield k _v=0,j=0 ≈ 2 × 10^-9 cm³/molecule/s at T = 10 K which is consistent with a T-independent Langevin rate k _L = 2.1 × 10^-9 cm³/molecule/s but considerably larger than the only experimentally reported value k _ICR = (9.1 ± 2.5) × 10^-10 cm³/molecule/s from ion cyclotron resonance experiments. Similarly, branching ratios for the reaction outcomes are comparable for the two PESs. However, when analyzing less averaged properties such as initial state-selected T-dependent rate coefficients and final vibrational states of the H₂ ⁺ product for low temperatures, the differences in the two PESs manifest themselves in the observables. Thus, depending on the property analyzed, accurate and globally valid representations of the PES are required, whereas more approximate and empirical construction schemes can be followed for state-averaged observables.