Statistical modelling of NH⁺/ND⁺ + H₂/HD/D₂ branching ratios

In this work we study hydrogen isotope fractionation along the proposed gas-phase ammonia formation pathway at temperatures relevant to interstellar modelling. Specifically we study the isotopically substituted variants of the NH⁺ + H₂ reaction, for which the primary product is NH₂⁺ + H. Adiabatic capture theory calculations have been performed for the association reaction. A new potential-energy surface has been determined for the NH₂⁺ product. An extensive set of rovibrational energy levels has been calculated for the NH₂⁺ isotopologues described by this potential-energy surface. These rovibrational energy levels have been used to determine energy- and angular-momentum-resolved numbers of accessible product states from the NH₃⁺ isotopologue intermediates, which in turn have been averaged to give statistical branching fractions for all isotopically-substituted NH⁺ + H₂ reactions. It is determined that in all cases where both hydrogen and deuterium are present, the NHD⁺ product is preferred.