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.

Optimized Structural Data at the Complete Basis Set Limit via Successive Quadratic Minimizations.

Two variants of a successive quadratic minimization method (SQM and c-SQM) are suggested to calculate the structural properties of molecular systems at the complete basis set (CBS) limit. When applied to H3+, H2O, CH2O, SH2, and SO2, they revealed CBS/(x1, x2) structural parameters that significantly surpass the raw ones calculated at the x2 basis set level. Such a performance has also been verified for the intricate case of the water dimer. Because the c-SQM method is system specific, thus showing somewhat enhanced results relative to the general SQM protocol, it can be of higher cost depending on the level of calibration used. Yet, it hardly surpasses the general quality of the results obtained with the cost-effective SQM method. Since the number of cycles required to reach convergence is relatively small, both schemes are simple to use and easily adaptable to any of the existing extrapolation schemes for the Hartree-Fock and correlation energies.

Accurate DMBE potential-energy surface for CNO(2A″) and rate coefficients in C(3P)+NO collisions.

A realistic double many-body expansion potential energy surface (PES) is developed for the 2A″ state of the carbon-nitrogen-oxygen (CNO) system based on MRCI-F12/cc-pVQZ-F12 ab initio energies. The new PES reproduces the fitted points with chemical accuracy (root mean square deviation up to 0.043 eV) and explicitly incorporates long range energy terms that can accurately describe the electrostatic and dispersion interactions. Thermal rate coefficients were computed for the C(3P) + NO(2Π) reaction for temperatures ranging from 15 K to 10 000 K, and the values are compared to previously reported results. The differences are rationalized, and the major importance of long range forces in predicting the rate coefficients for barrierless reactions is emphasized.

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 HN2 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 HN2 predicted to have a bent global minimum 36 kcal mol-1 below the N(4S) + NH(X̃3Σ-) 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.

Quasiclassical Study of the C(3P) + NO(X2Π) and O(3P) + CN(X2Σ+) Collisional Processes on an Accurate DMBE Potential Energy Surface.

The predicted rate constants for C + NO and O + CN collisions in three potential energy surfaces (PESs) for the 2A' state of the CNO molecule are compared using quasiclassical trajectories. Different temperature dependencies are obtained for the C + NO reaction, which are explained in terms of the long-range properties of the PESs. Recommended values and mechanistic details are also reported. For O + CN collisions, a better agreement between the theoretical results is found, except for temperatures below 100 K.

The O + NO( v) Vibrational Relaxation Processes Revisited.

We have carried out a quasiclassical trajectory study of the O + NO( v) energy transfer process using DMBE potential energy surfaces for the ground-states of the 2A' and 2A″ manifolds. State-to-state vibrational relaxation rate constants have been computed over the temperature range 298 and 3000 K and initial vibrational states between v = 1 and 9. The momentum-Gaussian binning approach has been employed to calculate the probability of the vibrational transitions. A comparison of the calculated state-to-state rate coefficients with the results from experimental studies and previous theoretical calculations shows the relevance of the 1 2A″ potential energy surface to the title vibrational relaxation process.

Accurate Explicit-Correlation-MRCI-Based DMBE Potential-Energy Surface for Ground-State CNO.

We report a new global double many-body expansion potential energy surface for the ground state of the CNO(2A') manifold, calculated by the explicit correlation multireference configuration interaction method. The functional form was accurately fitted to 3701 ab initio points with a root mean squared deviation of 0.99 kcal mol-1. All stationary points reported in previous forms are systematically described and improved, in addition to three new ones and a characterization of an isomerization transition state between the CNO and NCO minima. The novel proposed form gives a realistic description of both short-range and long-range interactions and hence is commended for dynamics studies.

Modeling cusps in adiabatic potential energy surfaces.

A method for modeling cusps on adiabatic potential energy surfaces without the need for any adiabatic-to-diabatic transformation is presented and shown to be successfully applied to the (2)A″ state of NO2. The more complicated case of a system with permutationally equivalent crossing seams is also examined and illustrated by considering the two first (2)A' states of the nitrogen trimer.

Single-sheeted double many-body expansion potential energy surface for ground-state ClO2.

A global single-sheeted double many-body expansion potential energy surface is reported for the ground electronic state of ClO2. The potential energy surface is obtained by fitting 3200 energy points that map all atom-diatom dissociation channels as well as all relevant stationary points, including the well-known OClO and ClOO structures. The ab initio calculations are obtained at the multireference configuration interaction level of theory, employing the cc-pVXZ (X = D, T) Dunning basis sets, and then extrapolated to the complete basis set limit with the generalized uniform singlet- and triplet-pair protocol. The topographical features of the novel global potential energy surface are examined in detail.

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.

Ab initio-based global double many-body expansion potential energy surface for the first 2A" electronic state of NO2.

An ab initio-based global double many-body expansion potential energy surface is reported for the first electronic state of the NO(2)((2)A") manifold. Up to 1700 ab initio energies have been employed to map the full configuration space of the title molecule, including stationary points and asymptotic channels. The calculated grid energies have been scaled to account for the incompleteness of the basis set and truncation of the MRCI expansion and fitted analytically with a total root-mean-squared-deviation smaller than 1.0 kcal mol(-1). The lowest point of the potential energy surface corresponds to the (2)B(1) linear minimum, which is separated from the C(s) distorted minimum by a C(2v) barrier of ≈9.7 kcal mol(-1). As usual, the proposed form includes a realistic representation of long-range interactions. Preliminary work indicates its reliability for reaction dynamics.

A theoretical study of rate coefficients for the O + NO vibrational relaxation.

Quasi-classical trajectories have been integrated to study the vibrational relaxation of the O + NO(v) process as a function of the initial vibrational quantum number for T = 298 K, 1500 K, and 3000 K. Two reliable potential energy surfaces have been employed for the A' and A'' doublet states of NO2. The calculated vibrational relaxation rate constants show a nearly v-independent behavior at room temperature and a moderate increase with v for higher temperatures. Although deviating significantly from the recommended values, good agreement with recent experimental results has been obtained. The importance of multi-quantum transitions is also analyzed.