Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.