A new formulation of three-body dynamical correlation energy for explicit potential functions

Quasiclassical Trajectory Study of the Si + SH Reaction on an Accurate Double Many-Body Expansion Potential Energy Surface

An accurate potential energy surface (PES) for the HSiS system based on MRCI+Q calculations extrapolated to the complete basis set limit is presented. Modeled with the double many-body expansion (DMBE) method, the PES provides an accurate description of the long-range interactions, including electrostatic and dispersion terms decaying as R^-4, R^-5, R^-6, R^-8, R^-10 that are predicted from dipole moments, quadrupole moments, and dipolar polarizabilities, which are also calculated at the MRCI+Q level. The novel PES is then used in quasiclassical trajectory calculations to predict the rate coefficients of the Si + SH → SiS + H reaction, which has been shown to be a major source of the SiS in certain regions of the interstellar medium. An account of the zero-point energy leakage based on various nonactive models is also given. It is shown that the reaction is dominated by long-range forces, with the mechanism Si + SH → SiSH → SSiH → SiS + H being the most important one for all temperatures studied. Although SSiH corresponds to the global minimum of the PES, the contribution from the direct reaction Si + SH → SSiH → SiS + H is less than 0.5% for temperatures higher than 500 K. The rovibrational distributions of the products are calculated using the momentum Gaussian binning method and show that as the temperature is increased the average vibrational quantum number decreases while the rotational distribution spreads up to larger values.

Accurate Potential Energy Surface for Quartet State HN2 and Interplay of N(4S) + NH(X̃3Σ-) versus H + N2(A3Σu+) Reactions

A global potential energy surface for the lowest quartet state of HN₂ is reported for the first time from accurate multireference ab initio calculations extrapolated to the complete basis set limit using the double many-body expansion method. All its stationary points are characterized, with the lowest quartet of HN₂ predicted to have a bent global minimum 36 kcal mol^-1 below the N(⁴S) + NH(X̃³Σ^-) asymptote, from which it is barrierlessly achievable. The entire set of calculated ab initio points has been fitted for energies up to 1000 kcal mol^-1 above the global minimum with an RMSD of 0.89 kcal mol^-1, a gap comprising all identified stationary points. Special care is taken in modeling the involved long-range forces and cusps caused by crossing seams. The novel PES prompts for the calculation of rate constants for several unexplored reactions that are relevant for combustion, plasma, and atmospheric chemistry.

Dynamics of the O + ClO reaction: reactive and vibrational relaxation processes

Classical trajectories have been integrated to study the O + ClO reaction, both reactive and vibrational energy transfer processes, for the range of temperatures 100 ≤ T/K ≤ 500 using momentum Gaussian binning. The employed potential energy surface is the recently proposed single-sheeted double many-body expansion potential energy surface for the (2)A" ground-state of ClO2 based on multireference ab initio data. A capture-type regime with a room-temperature rate constant of (17.8 ± 0.5) × 10(-12) cm(3) s(-1) and temperature dependence of k(T/K)/cm(3) s(-1) = 22.4 × 10(-12) × T(-0.81) exp(-39.2/T) has been found. Although the value reported here is half of the experimental and recommended one, tentative explanations are given. Other dynamical attributes are also examined for the title reaction, with state-to-all and state-to-state vibrational relaxation and excitation rate constants reported for temperatures of relevance in stratospheric chemistry.

Quantum wave packet and quasiclassical trajectory studies of the reaction H(2S) + CH(X2 Π; v = 0, j = 1) → C(1D) + H2 (X1 Σg+): Coriolis coupling effects and stereodynamics

The quantum mechanics (QM) and quasiclassical trajectory (QCT) calculations have been carried out for the title reaction with the ground minimal allowed rotational state of CH (j = 1) on the 1 (1)A' potential energy surface. For the reaction probability at total angular momentum J = 0, a similar trend of the QM and QCT calculations is observed, and the QM results are larger than the latter almost in the whole considered energy range (0.1-1.5 eV). The QCT integral cross sections are larger than the QM results with centrifugal sudden approximation, while smaller than those from QM method including Coriolis coupling for collision energies bigger than 0.25 eV. The quantum wave-packet computations show that the Coriolis coupling effects get more and more pronounced with increasing of J. In addition to the scalar properties, the stereodynamical properties, such as the average rotational alignment factor <P2 (j'•k)>, the angular distributions P(θr ), P(ϕr ), P(θr ,ϕr ), and the polarization-dependent generalized differential cross sections have been explored in detail by QCT approach.

QUASICLASSICAL TRAJECTORY STUDY OF PRODUCT POLARIZATIONS IN THE H(2S) + NH → N(4S) + H2 REACTION

Based on the DMBE potential surface of the 4 A ″ ground-state, the product rotational polarizations in the title reaction are studied by using quasiclassical trajectory (QCT) calculation method. Three angular distributions of P(θr), P(Φr), P(θr, Φr) and the four polarization-dependent differential cross sections (PDDCSs) were calculated for the collision energy range of 1–20 kcal/mol. The results revealed that the product is backward-scattering and the product rotational angular momentum j′ is aligned and oriented. With the increment of collision energy, the degree of the product alignment and orientation is enhanced, showing the collision energy-dependent behaviors of the product polarizations.

New Double Many-Body Expansion Potential Energy Surface for Ground-State HCN from a Multiproperty Fit to Accurate ab Initio Energies and Rovibrational Calculations

An accurate single-sheeted double many-body expansion potential energy surface has been obtained for the ground electronic state of the hydrogen cyanide molecule via a multiproperty fit to ab initio energies and rovibrational data. This includes 106 rovibrational levels and 2313 discrete points, which are fit with a rmsd of 4 cm(-1) and 2.42 kcal mol(-1), respectively, and seven zero first-derivatives that are reproduced at three stationary points. Since the potential also describes accurately the appropriate asymptotic limits at the various dissociation channels, it is commended both for the simulation of rovibrational spectra and reaction dynamics.

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.