Quasi-Classical Trajectory Calculations of the Hydrogen Abstraction Reaction H + NH3

New Method To Extract Final-State Information of Polyatomic Reactions Based on Normal Mode Analysis

State-to-state reaction dynamics provides a comprehensive insight into reaction mechanisms of chemical reactions at the atomic level. A new scheme to extract final-state information based on normal mode analysis is proposed in this work. Different from the traditional scheme extracting the coordinates and momenta from the last step of each trajectory, they are taken in the new scheme from a specific step of each reactive trajectory within the last vibrational period of the product molecule by demanding the corresponding geometry of the step to have the minimum potential energy. Test calculations on the collisions between the atom H and the molecules H₂O, H₂S, and NH₃ show that the new scheme works much better than the traditional one. In addition, the new scheme is applied to calculate the vibrational state distribution of the product NH₂ in the reaction H + NH₃ → H₂ + NH₂.

A nine-dimensional global potential energy surface for NH4(X2A1) and kinetics studies on the H + NH3 ↔ H2 + NH2 reaction

A nine-dimensional global potential energy surface (PES) for the NH₄ system is developed from ∼10⁵ high-level ab initio points and the hydrogen abstraction kinetics on the PES agree with experiment.

Constructing Potential Energy Surfaces for Polyatomic Systems: Recent Progress and New Problems

Different methods of constructing potential energy surfaces in polyatomic systems are reviewed, with the emphasis put on fitting, interpolation, and analytical (defined by functional forms) approaches, based on quantum chemistry electronic structure calculations. The different approaches are reviewed first, followed by a comparison using the benchmark H + CH4 and the H + NH3 gas-phase hydrogen abstraction reactions. Different kinetics and dynamics properties are analyzed for these reactions and compared with the available experimental data, which permits one to estimate the advantages and disadvantages of each method. Finally, we analyze different problems with increasing difficulty in the potential energy construction: spin-orbit coupling, molecular size, and more complicated reactions with several maxima and minima, which test the soundness and general applicability of each method. We conclude that, although the field of small systems, typically atom-diatom, is mature, there still remains much work to be done in the field of polyatomic systems.

Kinetics and dynamics of the NH3 + H → NH2 + H2 reaction using transition state methods, quasi-classical trajectories, and quantum-mechanical scattering

On a recent analytical potential energy surface developed by two of the authors, an exhaustive kinetics study, using variational transition state theory with multidimensional tunneling effect, and dynamics study, using both quasi-classical trajectory and full-dimensional quantum scattering methods, was carried out to understand the reactivity of the NH(3) + H → NH(2) + H(2) gas-phase reaction. Initial state-selected time-dependent wave packet calculations using a full-dimensional model were performed, where the total reaction probabilities were calculated for the initial ground vibrational state and for four excited vibrational states of ammonia. Thermal rate constants were calculated for the temperature range 200-2000 K using the three methods and compared with available experimental data. We found that (a) the total reaction probabilities are very small, (b) the symmetric and asymmetric N-H stretch excitations enhance the reactivity, (c) the quantum-mechanical calculated thermal rate constants are about one order of magnitude smaller than the transition state theory results, which reproduce the experimental evidence, and (d) quasi-classical trajectory calculations, which were performed with the main goal of analyzing the influence of the zero-point energy problem on the final dynamics results, reproduce the quantum scattering calculations on the same surface.