Production of NH(ND) radicals in the reactions of N(22D) with H2(D2): Nascent vibrational distributions of NH(X 3Σ−) and ND(X 3Σ−)

An accurate NH2(X2A′′) CHIPR potential energy surface via extrapolation to the complete basis set limit and dynamics of the N(2D) + H2(X1+g) reaction

An accurate CHIPR potential energy surface for NH2(X2A′′) is structured to study the N(2D) + H2(X1Σ+g) reaction using the time-dependent wave packet and quasi-classical trajectory method.

Accurate Potential Energy Surfaces for the Three Lowest Electronic States of N(2D) + H2(X1∑g +) Scattering Reaction

The three lowest full three-dimensional adiabatic and three diabatic global potential energy surfaces are reported for the title system. The accurate ab initio method (MCSCF/MRCI) with larger basis sets (aug-cc-pVQZ) is used to reduce the adiabatic potential energies, and the global adiabatic potential energy surfaces are deduced by a three-dimensional B-spline fitting method. The conical intersections and the mixing angles between the lowest three adiabatic potential energy surfaces are precisely studied. The most possible nonadiabatic reaction pathways are predicted, i.e., N(2D) + H2(X1∑g +) → NH2(22A') → CI (12A'-22A') → NH2(12A') → CI (12A″-12A') → NH2(12A″) → NH(X3∑-) + H(2S). The products of the first excited state (NH(a1Δ) + H(2S)) and the second excited state (NH(b1∑g +) + H(2S)) can be generated in these nonadiabatic reaction pathways too.

The effect of intersystem crossings in N((2)D) + H2 collisions

The transitions between quartet and doublet states of the NH2 molecule are studied for the first time, allowing the evaluation of the N((4)S) + H2 reactive channel. High level ab initio calculations of the spin-orbit coupling are performed over the whole configurational space of the NH2 molecule and fitted to a proposed analytic form. Quasiclassical trajectories coupled with the surface hopping method are employed to calculate reaction cross section and rate constants. The reaction is largely affected by the initial rovibrational states of H2, while the formation of long-lived complexes enhances the reaction probability.

Accurate ab initio-based adiabatic global potential energy surface for the 2(2)A" state of NH2 by extrapolation to the complete basis set limit

A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N((2)D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N((2)D)+H2(X(1)Σg (+))(ν=0,j=0)→NH(a(1)Δ)+H((2)S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction.

Accurate double many-body expansion potential energy surface by extrapolation to the complete basis set limit and dynamics calculations for ground state of NH2

An accurate single-sheeted double many-body expansion potential energy surface is reported for the title system. A switching function formalism has been used to warrant the correct behavior at the H2(X1Σg+)+N(2D) and NH (X3Σ-)+H(2S) dissociation channels involving nitrogen in the ground N(4S) and first excited N(2D) states. The topographical features of the novel global potential energy surface are examined in detail, and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. The novel surface can be using to treat well the Renner-Teller degeneracy of the 12A″ and 12A' states of NH 2. Such a work can both be recommended for dynamics studies of the N(2D)+H2 reaction and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen-containing systems. In turn, a test theoretical study of the reaction N(2D)+H2(X1Σg+)(ν=0,j=0)→NH (X3Σ-)+H(2S) has been carried out with the method of quantum wave packet on the new potential energy surface. Reaction probabilities, integral cross sections, and differential cross sections have been calculated. Threshold exists because of the energy barrier (68.5 meV) along the minimum energy path. On the curve of reaction probability for total angular momentum J = 0, there are two sharp peaks just above threshold. The value of integral cross section increases quickly from zero to maximum with the increase of collision energy, and then stays stable with small oscillations. The differential cross section result shows that the reaction is a typical forward and backward scatter in agreement with experimental measurement result.

VECTOR CORRELATIONS AND PRODUCT POLARIZATIONS IN THE N(2D) + D2 → ND + D REACTIVE SYSTEM

The vector correlations between products and reagents of the N (2D) + D 2 reaction are investigated by employing quasi-classical trajectory (QCT) calculation on the accurate DMBE potential energy surface (PES) of the 2A″ state. Stereo-dynamic quantities, including the four generalized polarization-dependent differential cross-sections (PDDCSs), the angular distribution P(θr), the dihedral-angle distribution P(φr), as well as the product rotational angular distribution in the polar form of P(θr, φr), are calculated in the center-of-mass (CM) frame. The results indicate that the product rotational angular momentum j′ not only aligns along the y-axis, but also orients to the negative direction of the y-axis. The isotope effect in the context of chemical stereo-dynamics and influences of different versions of ground-state PESs on vector correlations are shown and discussed.

QUASI-CLASSICAL TRAJECTORY STUDY OF THE REACTIONS N(2D) WITH H2, D2, AND HD

Quasi-classical trajectory (QCT) calculations are carried out for the title reactions on the potential energy surface (PES) of Ho et al.1 Our calculated integral cross-section values have been compared with the recent two quantum mechanics (QM) ones: they are close to those of one QM calculation in the high collision energy range, but they approach to another one in the low collision energy range. The product rotational alignments 〈P2 (J' ⋅ K)〉 have also been calculated.

Differential and Integral Cross Sections of the N(2D) + H2 → NH + H Reaction from Exact Quantum and Quasi-Classical Trajectory Calculations

Exact quantum mechanical state-to-state differential and integral cross sections and their energy dependence have been calculated on an accurate NH2 potential energy surface (PES), using a newly proposed Chebyshev wave packet method. The NH product is found to have a monotonically decaying vibrational distribution and an inverted rotational distribution. The product angular distributions peak in both forward and backward directions, but with a backward bias. This backward bias is more pronounced than that observed previously on a less accurate PES. Both the differential and integral cross sections oscillate mildly with collision energy, indicating the dominance of short-lived resonances. The quantum mechanical results are compared with those obtained from quasi-classical trajectories. The agreement is generally reasonable and the discrepancies can be attributed to the neglect of quantum effects such as tunneling. Detailed analysis of the trajectories revealed that the backward bias in the differential cross section stems overwhelmingly from the fast insertion component of the reaction, augmented with some flux from the abstraction channel, particularly at high collision energies.

Dynamics of insertion reactions of H2 molecules with excited atoms

Recent progress in the study of insertion reactions of hydrogen molecules with excited atoms is reviewed in this article. In particular, the dynamics of the reaction of O(1D), N(2D), C(1D), and S(1D) with H2 and its isotopomers, which have received a great deal of attention over the past decade, are examined in detail. All of these systems have in common the existence of several potential energy surfaces (PES) correlating with the reagents' states, and consequently, they can give rise to reaction following different adiabatic and nonadiabatic pathways. The main contribution, however, arises from their ground singlet PESs which feature the existence of deep wells with small or null barriers for insertion. Accordingly, these reactions proceed mainly via formation of relatively long-lived collision complexes and display an overall nearly statistical behavior. In spite of their similarities, the various reactions have peculiar characteristics caused by important differences of their respective PESs. The contribution of excited PES to the global reactivity, which has also become an important issue and a challenge both for theory and experiment, is also examined. The different theoretical approaches are discussed in the text, along with the experimental results obtained by a variety of techniques. The recent exact quantum treatments of these reactive systems together with the development of a rigorous statistical model have contributed to a very accurate description which in many cases matches very well the detailed measurements. The quasi-classical trajectory (QCT) method has also provided a fairly accurate description of the reaction dynamics for these systems. In particular, the analysis in terms of collision times has yielded interesting clues about the reaction mechanisms.

Exact quantum dynamics of N(D2)+H2→NH+H reaction: Cross-sections, rate constants, and dependence on reactant rotation

Using an exact Chebyshev wave packet method, initial state-specified (upsilon(i)=0, j(i)=0,2) integral cross-sections and rate constants are obtained for the title reaction on the latest ab initio potential energy surface. Reaction probabilities up to J=29 are dependent on the reactant rotation and show mild oscillations superimposed on a broad background. Due to a barrier in the entrance channel, the cross sections increase with energy with clear thresholds and the rate constants vary with temperature in the Arrhenius form. The calculated canonical rate constant is in good agreement with the experimental measurements. Our results also indicate that the quasiclassical trajectory method underestimates the rate due to the neglect of tunneling, while the quantum statistical approach overestimates because of the short lifetime of the reaction intermediate.

Influence of rotation and isotope effects on the dynamics of the N(D2)+H2 reactive system and of its deuterated variants

Integral cross sections and thermal rate constants have been calculated for the N((2)D)+H(2) reaction and its isotopic variants N((2)D)+D(2) and the two-channel N((2)D)+HD by means of quasiclassical trajectory and statistical quantum-mechanical model methods on the latest ab initio potential-energy surface [T.-S. Ho et al., J. Chem. Phys. 119, 3063 (2003)]. The effect of rotational excitation of the diatom on the dynamics of these reactions has been investigated and interesting discrepancies between the classical and statistical model calculations have been found. Whereas a net effect of reagent rotation on reactivity is always observed in the classical calculations, only a very slight effect is observed in the case of the asymmetric N((2)D)+HD reaction for the statistical quantum-mechanical method. The thermal rate constants calculated on this Potential-Energy Surface using quasiclassical trajectory and statistical model methods are in good agreement with the experimental determinations, although the latter are somewhat larger. A reevaluation of the collinear barrier of the potential surface used in the present study seems timely. Further theoretical and experimental studies are needed for a full understanding of the dynamics of the title reaction.

A Quantum Wave Packet Dynamics Study of the N(2D) + H2 Reaction

We report a dynamics study of the reaction N((2)D) + H(2) (v=0, j=0-5) --> NH + H using the time-dependent quantum wave packet method and a recently reported single-sheeted double many-body expansion potential energy surface for NH(2)(1(2)A' ') which has been modeled from accurate ab initio multireference configuration-interaction calculations. The calculated probabilities for (v=0, j=0-5) are shown to display resonance structures, a feature also visible to some extent in the calculated total cross sections for (v=0, j=0). A comparison between the calculated centrifugal-sudden and coupled-channel reaction probabilities validate the former approximation for the title system. Rate constants calculated using a uniform J-shifting scheme and averaged over a Boltzmann distribution of rotational states are shown to be in good agreement with the available experimental values. Comparisons with other theoretical results are also made.

Computational study of the reaction of N(D2) atoms with CH2F radicals: An example of a barrier-free reaction involving very high internal energies

The singlet potential-energy surface for the N(2D)+CH2F(2A') reaction has been studied employing both second-order Møller-Plesset and density-functional theories. The energies of the involved species have been refined using the Gaussian-2, complete basis set, and coupled-cluster singles and doubles (triples) methods. The reaction proceeds through the formation of an initial intermediate, which does not involve any activation barrier. Based on the energy profile for the singlet potential-energy surface, the preferred product should be the most exothermic one, namely, HCN+HF, followed by HNC+HF and FCN+H2. This result seems in contradiction with a computational study of the kinetics of the title reaction in terms of the statistical theories, which leads to the prediction that the production of HNC+HF should be the dominant channel. Consequently, a limited molecular-dynamics study has been carried out, concluding that in fact the system behaves in a nonstatistical way. According to the molecular-dynamics study, the most exothermic channel, HCN+HF, should be the dominant one. An analysis of the possible role of the singlet surface in the reaction of N(4S) with CH2F(2A') has also been carried out. The computational study shows that the microcanonical coefficients for the nonadiabatic channels are much smaller than the competing adiabatic ones. Therefore, the reaction of N(4S) with CH2F(2A') should proceed on the triplet surface without spin change.

Repulsive double many-body expansion potential energy surface for the reactions N(4S)+ H2⇌ NH(X3Σ–)+ H from accurate ab initio calculations

A single-sheeted DMBE potential energy surface is reported for the reactions N(4S)+H2<-->NH(X3Sigma-)+H based on a fit to accurate multireference configuration interaction energies. These have been calculated using the aug-cc-pVQZ basis set of Dunning and the full valence complete active space wave function as reference, being semi-empirically corrected by scaling the two-body and three-body dynamical correlation energies. The topographical features of the novel global potential energy surface are examined in detail, including a conical intersection involving the two first 4A'' potential energy surfaces which has been transformed into an avoided crossing in the present single-sheeted representation.