Ab initio rate constants from hyperspherical quantum scattering: application to H + CH4 --> H2 + CH3

A general and practical procedure is described for calculating rate constants for chemical reactions using a minimal number of ab initio calculations and quantum-dynamical computations. The method exploits a smooth interpolating functional developed in the hyperspherical representation. This functional is built from two Morse functions and depends on a relatively small number of parameters with respect to conventional functionals developed to date. Thus only a small number of ab initio points needs to be computed. The method is applied to the H + CH4 --> H2 + CH3 reaction. The quantum scattering calculations are performed treating explicitly the bonds being broken and formed. All the degrees of freedom except the breaking and forming bonds are optimized ab initio and harmonic vibrational frequencies and zero-point energies for them are calculated at the MP2(full) level with a cc-pVTZ basis set. Single point energies are calculated at a higher level of theory with the same basis set, namely CCSD(T, full). We report state-to-state cross sections and thermal rate constants for the title reaction and make comparisons with previous results. The calculated rate constants are in good agreement with experiments.     

Quasiclassical trajectory study of the F + CH4 reaction dynamics on a dual-level interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the F + CH4 reactive system has been constructed using the interpolation method of Collins and co-workers. The ab initio calculations have been performed using second-order M?ller-Plesset (MP2) perturbation theory to build the initial PES. Scaling all correlation (SAC) methodology has been employed to improve the ab initio calculations and to construct a dual-level PES. Using this PES, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations and internal energy distributions has been carried out for the F + CH4 and F + CD4 reactions and the theoretical results have been compared with the available experimental data.     

Ab initio potential energy surface, variational transition state theory, and quasiclassical trajectory studies of the F+CH4-->HF+CH3 reaction

In this work we present a study of the F+CH(4)-->HF+CH(3) reaction (DeltaHdegrees(298 K)=-32.0 kcal mol(-1)) using different methods of the chemical reaction theory. The ground potential energy surface (PES) is characterized using several ab initio methods. Full-dimensional rate constants have been calculated employing the variational transition state theory and using directly ab initio data. A triatomic analytical representation of the ground PES was derived from ab initio points calculated at the second- and fourth-order M?ller-Plesset levels with the 6-311+G(2df,2pd) basis set, assuming the CH(3) fragment to be a 15 a.m.u. pseudoatom in the fitting process. This is suggested from experiments that indicate that the methyl group is uncoupled to the reaction coordinate. A dynamics study by means of the quasiclassical trajectory (QCT) method and employing this analytical surface was also carried out. The experimental data available on the HF internal states distributions are reproduced by the QCT results. Very recent experimental information about the reaction stereodynamics is also borne out by our QCT calculations. Comparisons with the benchmark F+H(2) and analogous Cl+CH(4) reactions are established throughout.     

Reduced dimensionality spin-orbit dynamics of CH3 + HCl ⇌ CH4 + Cl on ab initio surfaces

A reduced dimensionality quantum scattering method is extended to the study of spin-orbit nonadiabatic transitions in the CH(3) + HCl ? CH(4) + Cl((2)P(J)) reaction. Three two-dimensional potential energy surfaces are developed by fitting a 29 parameter double-Morse function to CCSD(T)/IB//MP2/cc-pV(T+d)Z-dk ab initio data; interaction between surfaces is described by geometry-dependent spin-orbit coupling functions fit to MCSCF/cc-pV(T+d)Z-dk ab initio data. Spectator modes are treated adiabatically via inclusion of curvilinear projected frequencies. The total scattering wave function is expanded in a vibronic basis set and close-coupled equations are solved via R-matrix propagation. Ground state thermal rate constants for forward and reverse reactions agree well with experiment. Multi-surface reaction probabilities, integral cross sections, and initial-state selected branching ratios all highlight the importance of vibrational energy in mediating nonadiabatic transition. Electronically excited state dynamics are seen to play a small but significant role as consistent with experimental conclusions.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

Theoretical study of H-abstraction reactions from CH3Cl and CH3Br molecules by ClO and BrO radicals

The rate constants of the H-abstraction reactions from CH(3)Cl and CH(3)Br molecules by ClO and BrO radicals have been estimated over the temperature range of 300-2500 K using four different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using B3LYP and MP2 methods combined with the cc-pVTZ basis set. Single-point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using the cc-pVTZ and cc-pVQZ basis sets. Canonical transition-state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature. In order to choose the appropriate levels of theory with chlorine- and bromine-containing species, the reference reaction Cl ((2)P(3/2)) + CH(3)Cl → HCl + CH(2)Cl (R(ref)) was first theoretically studied because its kinetic parameters are well-established from numerous experiments, evaluation data, and theoretical studies. The kinetic parameters of the reaction R(ref) have been determined accurately using the CCSD(T)/cc-pVQZ//MP2/cc-pVTZ level of theory. This level of theory has been used for the rate constant estimation of the reactions ClO + CH(3)Cl (R(1)), ClO + CH(3)Br (R(2)), BrO + CH(3)Cl (R(3)), and BrO + CH(3)Br (R(4)). Six-parameter Arrhenius expressions have been obtained by fitting to the computed rate constants of these four reactions (including cis and trans pathways) over the temperature range of 300-2500 K.     

Jahn-Teller effect in tetrahedral symmetry: large-amplitude tunneling motion and rovibronic structure of CH4+ and CD4+

The energy level structures of the ground vibronic states of 12CH4+, 13CH4+, and 12CD4+ have been measured by pulsed-field-ionization zero-kinetic-energy photoelectron spectroscopy. The nuclear spin symmetries of the tunneling-rotational levels have been determined in double-resonance experiments via selected rotational levels of the v3=1 and v3=2 vibrational levels of the X 1A1 ground state of CH4. The energy level structures of 12CH4+, 13CH4+, and 12CD4+ have been analyzed with an effective tunneling-rotational Hamiltonian. The analysis together with a group theoretical treatment of the Tx(e+t2) Jahn-Teller effect in the Td(M) group prove that the equilibrium geometry of 12CH4+, 13CH4+, and 12CD4+ has C2v symmetry and characterize the pseudorotational dynamics in these fluxional cations. The tunneling behavior is discussed in terms of the relevant properties of the potential energy surface, some of which have been recalculated at the CCSD(T)/cc-pVTZ level of ab initio theory.     

Quasiclassical trajectory study of the C(1D) + H2 reaction and isotopomeric variants: kinetic isotope effect and CD/CH branching ratio

The recently proposed ab initio single-sheeted double many-body expansion potential energy for the methylene molecule has been used to perform quasiclassical trajectory (QCT) calculations for the title reaction. Thermal and initial state-specific (v = 0, j = 0) rate constants for the C((1)D) + H(2)/HD/D(2) reactions have been obtained over a wide range of temperatures. Cross sections for the reaction C((1)D) + H(2) and its deuterated isotopes have also been calculated, as well as the CD/CH branching ratios for the C((1)D) + HD reaction. It is found that the CD + H product channel in the C((1)D) + HD reaction is preferred relative to the CH + D channel. The estimated rate constants are predicted to be in the order k(H2) > k(HD) > k(D2) and the calculated cross sections and rate constants compared with available theoretical and experimental data.     

Searching for resonances in the reaction Cl+CH4-->HCl+CH3: quantum versus quasiclassical dynamics and comparison with experiments

A quantum-mechanical (QM) and quasiclassical trajectory (QCT) study was performed on the title reaction, using a pseudotriatomic ab initio based surface. Probabilities and integral cross sections present some clear peaks versus the collision energy E(col), which we assign to Feshbach resonances of the transition state, where the light H atom oscillates between the heavy Cl and CH(3) groups. For ground-state reactants, reactivity is essentially of quantum origin (QCT observables and oscillations are smaller, or much smaller, than QM ones), and the calculated integral cross section and product distributions are in reasonable agreement with the experiment. The reaction occurs through an abstraction mechanism, following both a direct and an indirect mechanism. The quasiclassical trajectory calculations show the participation of a short-lived collision complex in the microscopic reaction mechanism. Finally, QCT differential cross sections of Cl+CH(4)-->HCl (nu(')=0 and 1)+CH(3) oscillate versus E(col), whereas experimentally this only occurs for HCl (nu(')=1). This theoretical result and other oscillating properties found here could, however, be related to the existence of a Feshbach resonance for the production of HCl (nu(')=1), as suggested by experimentalists.     

An ab initio spin-orbit-corrected potential energy surface and dynamics for the F + CH4 and F + CHD3 reactions

We report an analytical ab initio three degrees of freedom (3D) spin-orbit-correction surface for the entrance channel of the F + methane reaction obtained by fitting the differences between the spin-orbit (SO) and non-relativistic electronic ground state energies computed at the MRCI+Q/aug-cc-pVTZ level of theory. The 3D model surface is given in terms of the distance, R(C-F), and relative orientation, Euler angles ? and θ, of the reactants treating CH(4) as a rigid rotor. The full-dimensional (12D) "hybrid" SO-corrected potential energy surface (PES) is obtained from the 3D SO-correction surface and a 12D non-SO PES. The SO interaction has a significant effect in the entrance-channel van der Waals region, whereas the effect on the energy at the early saddle point is only ～5% of that at the reactant asymptote; thus, the SO correction increases the barrier height by ～122 cm(-1). The 12D quasiclassical trajectory calculations for the F + CH(4) and F + CHD(3) reactions show that the SO effects decrease the cross sections by a factor of 2-4 at low collision energies and the effects are less significant as the collision energy increases. The inclusion of the SO correction in the PES does not change the product state distributions.     

Theoretical Studies of the Reaction Mechanisms of CH3S＋NO2

The potential energy surface for the CH3S NO2 reaction has been studied using the ab initio G3(MP2) method. A variety of possible complexes and saddle points along the minimum energy reaction paths have been characterized at UMP2 (full)/6-31G(d) level. The calculations reveal dominating reaction mechanisms of the title reaction: CH3S NO2 firstly produce intermediate CH3SONO,then break up into CH3SO NO. The results are valuable to understand the atmospheric sulfur compounds oxidation mechanism.     

Effect of the Ar-Ni(s) potential on the cross section for Ar+CH4/Ni{111} collision-induced desorption and the need for a more accurate CH4/Ni{111} potential

In a previous paper [L. Sun, P. de Sainte Claire, O. Meroueh, and W. L Hase, J. Chem. Phys. 114, 535 (2001)], a classical trajectory simulation was reported of CH(4) desorption from Ni{111} by Ar-atom collisions. At an incident angle theta(i) of 60 degrees (with respect to the surface normal), the calculated collision-induced desorption (CID) cross sections are in excellent agreement with experiment. However, for smaller incident angles the calculated cross sections are larger than the experimental values and for normal collisions, theta(i)=0 degrees , the calculated cross sections are approximately a factor of 2 larger. This trajectory study used an analytic function for the Ar+Ni(s) intermolecular potential which gives an Ar-Ni{111} potential energy minimum which is an order of magnitude too deep. In the work reported here, the previous trajectory study is repeated with an Ar+Ni(s) analytic intermolecular potential which gives an accurate Ar-Ni{111} potential energy minimum and also has a different surface corrugation than the previous potential. Though there are significant differences between the two Ar+Ni(s) analytic potentials, they have no important effects on the CID dynamics and the cross sections reported here are nearly identical to the previous values. Zero-point energy motions of the surface and the CH(4)-Ni(s) intermolecular modes are considered in the simulation and they are found to have a negligible effect on the CID cross sections. Calculations of the intermolecular potential between CH(4) and a Ni atom, at various levels of theory, suggest that there are substantial approximations in the ab initio calculation used to develop the CH(4)+Ni{111} potential. The implication is that the differences between the trajectory and experimental CID cross sections may arise from an inaccurate CH(4)+Ni{111} potential used in the trajectory simulation.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

Ab initio and direct quasiclassical-trajectory study of the Cl + CH4-->HCl + CH3 reaction

We present an electronic structure and dynamics study of the Cl + CH(4)--> HCl + CH(3) reaction. We have characterized the stationary points of the ground-state potential-energy surface using various electronic structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis-set limit based on geometries and harmonic frequencies obtained at the CCSD(T)/aug-cc-pvtz level, are in agreement with the experimental reaction energy and indirect measurements of the barrier height. Using ab initio information, we have reparametrized a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in various regions of the potential-energy surface. This improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories and characterize the reaction dynamics. The good agreement of the calculated HCl rotational and angular distributions with the experiment indicates that reparametrizing semiempirical Hamiltonians is a promising approach to derive accurate potential-energy surfaces for polyatomic reactions. However, excessive energy leakage from the initial vibrational energy of the CH(4) molecule to the reaction coordinate in the trajectory calculations calls into question the suitability of the standard quasiclassical-trajectory method to describe energy partitioning in polyatomic reactions.     

Quasiclassical trajectory study of energy transfer and collision-induced dissociation in hyperthermal Ar + CH4 and Ar + CF4 collisions

We present a study of energy transfer in collisions of Ar with methane and perfluoromethane at hyperthermal energies (E(coll) = 4-10 eV). Quasiclassical trajectory calculations of Ar + CX(4) (X = H, F) collisions indicate that energy transfer from reagents' translation to internal modes of the alkane molecule is greatly enhanced by fluorination. The reasons for the enhancement of energy transfer upon fluorination are shown to emerge from a decrease in the hydrocarbon vibrational frequencies of the CX(4) molecule with increasing the mass of the X atom, and to an increase of the steepness of the Ar-CX(4) intermolecular potential. At high collision energies, we find that the cross section of Ar + CF(4) collisions in which the amount of energy transfer is larger than needed to break a C-F bond is at least 1 order of magnitude larger than the cross sections of Ar + CH(4) collisions producing CH(4) with energy above the dissociation limit. In addition, collision-induced dissociation is detected in short time scales in the case of the fluorinated species at E(coll) = 10 eV. These results suggest that the cross section for degradation of fluorinated hydrocarbon polymers under the action of nonreactive hyperthermal gas-phase species might be significantly larger than that of hydrogenated hydrocarbon polymers. We also illustrate a practical way to derive intramolecular potential energy surfaces for bond-breaking collisions by improving semiempirical Hamiltonians based on grids of high-quality ab initio calculations.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single point calculations for all the stationary points were carried out at the QCISD(T)/ 6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

Direct dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2

We present a direct ab initio dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2, which is predicted to have four possible reaction channels caused by different attacking orientations of HO2 radical to CH2O. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of the four reaction channels are calculated at the B3LYP/cc-pVTZ level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of some single-point multilevel energy calculations (HL). The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that, in the whole temperature range, the more favorable reaction channels are Channels 1 and 3. The total ICVT/SCT rate constants of the four channels at the HL//B3LYP/cc-pVTZ level of theory are in good agreement with the available experiment data over the measured temperature ranges, and the corresponding three-parameter expression is k(ICVT/SCT) = 3.13 x 10(-20) T(2.70) exp(-11.52/RT) cm3 mole(-1) s(-1) in the temperature range of 250-3000 K. Additionally, the flexibility of the dihedral angle of H2O2 is also discussed to explain the different experimental values.     

One- and three-parameter density functional theory and high level ab initio theory study of the CH+CH disproportionation reaction

Complete basis set (CBS) ab initio computational studies were performed with the target being to explore the CH+CH potential energy surface. Several closed and open shell intermediates were located on the potential energy surface. Computed enthalpies for the branching reactions, as well as heats of formation are in excellent agreement. Although CBS computed energies are of high quality, this computational study is not capable of predicting the branching product ratio due to fact that neither the MP2 nor the 6-311G(2d,2p) basis set are sufficient to locate the reactant complexes and the transition state structures for the hydrogen and carbon transfer reactions in the reaction complexes. To properly explore the CH+CH potential energy surface a much higher ab initio theory level is required.     

Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

A systematic computational study on the reactions of HO2 with RO2: The HO2 + CH3O2(CD3O2) and HO2 + CH2FO2 reactions

A systematic theoretical study of the reactions of HO2 with RO2 has been carried out. The major concern of the present work is to gain insight into the reaction mechanism and then to explain experimental observations and to predict new product channels for this class of reactions of importance in the atmosphere. In this paper, the reaction mechanisms for two reactions, namely, HO2 + CH3O2 and HO2 + CH2FO2, are reported. Both singlet and triplet potential energy surfaces are investigated. The complexity of the present system makes it impossible to use a single ab initio method to map out all the reaction paths. Various ab initio methods including MP2, CISD, QCISD(T), CCSD(T), CASSCF, and density function theory (B3LYP) have been employed with the basis sets ranging from 6-31G(d) to an extrapolated complete basis set (CBS) limit. It has been established that the CCSD(T)/cc-pVDZ//B3LYP/6-311G(d,p) scheme represents the most feasible method for our systematic study. For the HO2 + CH3O2 reaction, the production of CH3OOH is determined to be the dominant channel. For the HO2 + CH2FO2 reaction, both CH2FOOH and CHFO are major products, whereas the formation of CHFO is dominant in the overall reaction. The computational findings give a fair explanation for the experimental observation of the products.