Accurate potential energy surface and quantum reaction rate calculations for the H+CH4-->H2+CH3 reaction

Calculations for the cumulative reaction probability N(E) (for J=0) and the thermal rate constant k(T) of the H+CH(4)-->H(2)+CH(3) reaction are presented. Accurate electronic structure calculations and a converged Shepard-interpolation approach are used to construct a potential energy surface which is specifically designed to allow the precise calculation of k(T) and N(E). Accurate quantum dynamics calculations employing flux correlation functions and multiconfigurational time-dependent Hartree wave packet propagation compute N(E) and k(T) based on this potential energy surface. The present work describes in detail the various convergence test performed to investigate the accuracy of the calculations at each step. These tests demonstrate the predictive power of the present calculations. In addition, approximate approaches for reaction rate calculations are discussed. A quite accurate approximation can be obtained from a potential energy surface which includes only interpolation points on the minimum energy path.     

Singlet and triplet potential surfaces for the O2+C2H4 reaction

Electronic structure calculations at the CASSCF and UB3LYP levels of theory with the aug-cc-pVDZ basis set were used to characterize structures, vibrational frequencies, and energies for stationary points on the ground state triplet and singlet O(2)+C(2)H(4) potential energy surfaces (PESs). Spin-orbit couplings between the PESs were calculated using state averaged CASSCF wave functions. More accurate energies were obtained for the CASSCF structures with the MRMP2/aug-cc-pVDZ method. An important and necessary aspect of the calculations was the need to use different CASSCF active spaces for the different reaction paths on the investigated PESs. The CASSCF calculations focused on O(2)+C(2)H(4) addition to form the C(2)H(4)O(2) biradical on the triplet and singlet surfaces, and isomerization reaction paths ensuing from this biradical. The triplet and singlet C(2)H(4)O(2) biradicals are very similar in structure, primarily differing in their C-C-O-O dihedral angles. The MRMP2 values for the O(2)+C(2)H(4)→C(2)H(4)O(2) barrier to form the biradical are 33.8 and 6.1 kcal/mol, respectively, for the triplet and singlet surfaces. On the singlet surface, C(2)H(4)O(2) isomerizes to dioxetane and ethane-peroxide with MRMP2 barriers of 7.8 and 21.3 kcal/mol. A more exhaustive search of reaction paths was made for the singlet surface using the UB3LYP/aug-cc-pVDZ theory. The triplet and singlet surfaces cross between the structures for the O(2)+C(2)H(4) addition transition states and the biradical intermediates. Trapping in the triplet biradical intermediate, following (3)O(2)+C(2)H(4) addition, is expected to enhance triplet→singlet intersystem crossing.     

An improved potential energy surface for the F+H2 reaction

A new ground state potential energy surface has been developed for the F+H₂ reaction. Using the UCCSD(T) method, ab initio calculations were performed for 786 geometries located mainly in the exit channel of the reaction. The new data was used to correct exit channel errors that have become apparent in the potential energy surface of Stark and Werner [J. Chem. Phys. 104 (1996) 6515]. While the entrance channel and saddlepoint properties of the Stark–Werner surface are unchanged on the new potential, the exit channel behavior is more satisfactory. The exothermicity on the new surface is much closer to the experimental value. The new surface also greatly diminishes the exit channel van der Waals well that was too pronounced on the Stark–Werner surface. Several preliminary dynamical scattering calculations were carried out using the new surface for total angular momentum equal to zero for F+H₂ and F+HD. It is found that gross features of the reaction dynamics are quite similar to those predicted by the Stark–Werner surface, in particular the reactive resonance for F+HD and F+H₂ survive. However, the most of the exit channel van der Waals resonances disappear on the new surface. It is predicted that the differential cross-sections at low collision energy for the F+H₂ reaction may be drastically modified from the predictions based on the Stark–Werner surface.     

Quasiclassical trajectory study of the reaction H+CH4(nu3 = 0,1)-->CH3+H2 using a new ab initio potential energy surface

Detailed quasiclassical trajectory calculations of the reaction H+CH4(nu3 = 0,1)-->CH3 + H2 using a slightly updated version of a recent ab initio-based CH5 potential energy surface [X. Zhang et al., J. Chem. Phys. 124, 021104 (2006)] are reported. The reaction cross sections are calculated at initial relative translational energies of 1.52, 1.85, and 2.20 eV in order to make direct comparison with experiment. The relative reaction cross section enhancement ratio due to the excitation of the C-H antisymmetric stretch varies from 2.2 to 3.0 over this energy range, in good agreement with the experimental result of 3.0 +/- 1.5 [J. P. Camden et al., J. Chem. Phys. 123, 134301 (2005)]. The laboratory-frame speed and center-of-mass angular distributions of CH3 are calculated as are the vibrational and rotational distributions of H2 and CH3. We confirm that this reaction occurs with a combination of stripping and rebound mechanisms by presenting the impact parameter dependence of these distributions and also by direct examination of trajectories.     

Thermochemistry and metastable energy release for the dissociation [C3H7]+ → [C3H5]+ + H2

The reverse activation energy, E_rev, for the dissociation [C₃H₇]⁺ → [C₃H₅]⁺ + H₂ has been determined as 0.24 ± 0.06 eV from measurements of the AP of [C₃H₅]+ produced by electron-impact from thermally generated sec-C₃H₇ radicals. The energy release observed in the corresponding metastable dissociation is 0.21 ± 0.01 eV, indicating that virtually all of E_rev is partitioned as translational Kinetic energy of the fragmentation products. The metastable ion peak shape is also discussed with respect to the evaluation of the energy release.     

Quasiclassical trajectory calculations of the reaction C+C2H2-->l-C3H, c-C3H+H, C3+H2 using full-dimensional triplet and singlet potential energy surfaces

Full-dimensional, density functional theory (B3LYP/6-311g(d,p))-based potential energy surfaces (PESs) are reported and used in quasi-classical calculations of the reaction of C with C(2)H(2). For the triplet case, the PES spans the region of the reactants, the complex region (with numerous minima and saddle points) and the products, linear(l)-C(3)H+H, cyclic(c)-C(3)H+H and c-(3)C(3)+H(2). For the singlet case, the PES describes the complex region and products l-C(3)H+H, c-C(3)H+H and l-(1)C(3)+H(2). The PESs are invariant under permutation of like nuclei and are fit to tens of thousands of electronic energies. Energies and harmonic frequencies of the PESs agree well the DFT ones for all stationary points and for the reactant and the products. Dynamics calculations on the triplet PES find both l-C(3)H and c-C(3)H products, with l-C(3)H being dominant at the energies considered. Limited unimolecular reaction dynamics on the singlet PES find both products in comparable amounts as well as the C(3)+H(2) product.     

The energy surface for isomeric [C2H3O]+ ions: Further experimental evidence

Mass spectra from collisionally activated dissociation (CAD) of [C₂H₃O]⁺ ions, including isotopically labeled analogs, provide further information on the isomers [CH₃C?O⁺] (a), [CH₂?C?O⁺H] (b), [⁺CH₂CH?O] (c) and (d). Our data generally support the recent conclusions from theory by Radom and coworkers and from experiment by Terlouw, Holmes and coworkers. Most acetyl-containing molecular ions form a ions in high purity only at low energies, consistent with isomerization of higher energy molecular ions to form the more stable enol which dissociates to b. Isomer d, prepared from (CICH₂)₂CHOH, undergoes facile hydrogen scrambling, presumably through a degenerate 1,2-hydrogen shift. Theory suggests that c undergoes spontaneous isomerization to a and d; although [C₂H₃O]⁺ ions from BrCH₂CHO appear to consist of a and ～15% d, the latter are formed without substantial hydrogen scrambling.     

Ab initio study of the potential energy surface for the OH+CO-->H+CO2 reaction

Potential energy surface for the reaction OH+CO-->H+CO2 has been calculated using the complete active space self-consistent-field and multireference configuration interaction methods with the correlation consistent triple-, quadruple-, and quintuple-zeta basis sets. A specific- reaction-parameters density functional theory has been suggested, in which the B3LYP functional is reoptimized to give the highly accurate potential energy surface with less computational efforts.     

Theoretical study of [Si,O,C,O] species: Prediction of new species on triplet potential energy surface

The intermediates [Si,O,C,O] of the Si + CO₂ reaction have been studied in detail using high level ab iniitio methods. Both singlet and triplet [Si,O,C,O] species are characterized structurally and energetically. On the singlet potential energy surface (PES), the vdw‐OSi–CO isomer and in the triplet PES, the bent‐SiOCO isomer is found to be thermodynamically as well as kinetically most stable species. All possible isomerization transition states (TS) are located on both singlet and triplet potential surfaces. On the triplet surface, the stability of the trans‐OSiCO isomer is comparable with that of the bent‐SiOCO isomer. A non‐planar cis‐SiOCO isomer is located on the triplet PES, which is predicted for the first time. Heats of formation at 0 K (Δ_fH°, 0 K) for all singlet and triplet species are computed using G3B3, G3MP2, and CBS‐Q theories. The discrepancy between G3B3 and the other two methods for the heat of formation value for triplet trans‐OSiCO is discussed. The PESs for singlet as well as triplet species with their dissociation asymptotes are explored at the CCSD(T)/6‐311G(d,p)//MP2/6‐311G(d,p) level of theory. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Deuterium isotope effects and mechanism of the gas phase reaction [C3H7]+→[C3H5]+ + H2

The deuterium kinetic isotope effect and the deuterium isotope effect upon kinetic energy release have been calculated for the loss of H₂ from the [C₃H₇]⁺ ion. The calculations are based on the transition state structure suggested recently from ab initio calculations on the reaction pathway. The results obtained are in good agreement with experimental data.     

Accurate double many-body expansion potential energy surface for triplet H3+. II. The upper adiabatic sheet (2(3)A')

We report on a global potential energy hypersurface for the upper sheet of the lowest triplet state of H3+. The analytic representation is based on the double many-body expansion theory. The ab initio data points, calculated with a large cc-pV5Z basis, are represented with a root mean square deviation of only 5.54 cm(-1) in the energy region below the H(+)+2H(2S) dissociation threshold. The quasi-bound vibronic states supported by this surface have also been calculated.     

An eight-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction on a new modified potential energy surface

An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES.     

Pi-donation and stabilizing effects of pnicogens in carbenium and silicenium ions: a theoretical study of [C(XH2)3]+ and [Si(XH2)3]+ (X = N, P, As, Sb, Bi)

Quantum chemical calculations at the MP2 and CCSD(T) levels of theory are reported for cations of the general type [A(XH2)3]+ with A = C, Si and X = N, P, As, Sb, Bi. Population analysis, methyl stabilization energies (MSEs), and structural criteria were used to predict the p(pi)-donor ability of and the pi-stabilization energy exerted by this series of pnicogens. All of the substituents XH2 considered in these studies invariably stabilize the triply substituted carbenium as well as the silicenium ions. The calculated data show that the intrinsic p(pi)-donation of the group 15 atoms follows the order N < P < As < Sb < Bi. However, the trend of the stabilization energies is fully reversed. The intrinsic stabilization energies of the planar carbenium ions decrease monotonically from 161.2 kcal mol(-1) for X = NH2 to 98.0 kcal mol(-1) for X = BiH2. The effective stabilization of the pnicogens in the equilibrium structures, which also includes the energy-demanding pyramidalization of the XH2 substituents, follows the same trend, although the absolute numbers are reduced to 145.6 kcalmol(-1) for X = NH2 and 53.2 kcalmol(-1) for X = BiH2. This seemingly contrasting behavior of increasing p(pi) charge donation and decreasing stabilization has already been found for other substituents. Previous studies have shown that carbenium ions substituted by chalcogens up to the fourth row also stabilize C+ less effectively with respect to heavier substituents. Of the ions investigated in this study, only the silicenium ions that are stabilized by pnicogens from the third to the sixth row of the periodic system yield increased stabilizing energies that follow the corresponding intrinsic p(pi)-donor abilities of the respective substituent.     

Exact quantum mechanical reaction probabilities for the collinear H + H2 reaction on a porter-karplus potential energy surface

Exact quantum mechanical calculation of the reaction probability for the collinear H + H₂ reaction on a Porter-Karplus potential energy surface are carried out by the finite-difference boundary value method at 6 energes in the threshold region and compared to close coupling, distorted wave, classical S matrix, transition state theory, and vibrational adiabatic calculations.     

Theoretical study of [Li(H2O)n]+ and [K(H2O)n]+ (n = 1−4) complexes

The geometries, successive binding energies, vibrational frequencies, and infrared intensities are calculated for the [Li(H₂O)_n]⁺ and [K(H₂O)_n]⁺ (n = 1?4) complexes. The basis sets used are 6-31G* and LANL 1DZ (Los Alamos ECP +DZ ) at the SCF and MP 2 levels. There is an agreement for calculated structures and frequencies between the MP 2/6-31G* and MP 2/LANL 1DZ basis sets, which indicates that the latter can be used for calculations of water complexes with heavier ions. Our results are in a reasonable agreement with available experimental data and facilitate experimental study of these complexes. © 1995 John Wiley & Sons, Inc.     

Kinetic study on the H + SiH4 abstraction reaction using an ab initio potential energy surface

Variational transition state theory calculations with the correction of multidimensional tunneling are performed on a 12-dimensional ab initio potential energy surface for the H + SiH(4) abstraction reaction. The surface is constructed using a dual-level strategy. For the temperature range 200-1600 K, thermal rate constants are calculated and kinetic isotope effects for various isotopic species of the title reaction are investigated. The results are in very good agreement with available experimental data.     

H+NH3反应势能面的量子化学研究

本文采用UCCSD(T)/aug-cc-pVTZ方法研究了H NH3反应势能面,获得了夺氢反应和交换反应过渡态的的几何结构和振动频率。夺氢反应的过渡态具有Cs对称性,其能垒为61.92 kJ/mol。交换反应的过渡态具有C3v对称性,其能垒为39.69 kJ/mol。H NH3发生形成Td对称性的反应中间体NH4里德堡自由基。与夺氢反应相比,交换反应具有更低的反应能垒,并且NH4自由基在反应中可形成长寿命的共振态,和夺氢反应形成竞争关系,因此在H NH3反应的量子动力学研究中必须同时考虑这两类反应。本文还采用更大的基组aug-cc-pVQZ和aug-cc-pV5Z研究了势能面对基组的收敛行为。     

Potential energy surface for the CCl4 + H --> CCl3 + ClH reaction: kinetics and dynamics study

An analytical potential energy surface for the gas-phase CCl4 + H --> CCl3 + ClH reaction was constructed with suitable functional forms to represent vibrational modes. This surface is completely symmetric with respect to the permutation of the four chlorine atoms and is calibrated with respect to experimental thermal rate constants available over the temperature range 297-904 K. On this surface, the thermal rate constants were calculated using variational transition-state theory with semiclassical transmission coefficients over a wider temperature range 300-2500 K, therefore obtaining kinetics information at higher temperatures than are experimentally available. This surface was also used to analyze dynamical features, such as tunneling and reaction-path curvature. In the first case, the influence of the tunneling factor is very small since a heavy chlorine atom has to pass through the barrier. In the second, it was found that vibrational excitation of the Cl-H stretching mode can be expected in the exit channel.     

Exact quantum scattering study of the H + HS reaction on a new ab initio potential energy surface H2S (3A")

We present an exact quantum dynamical study and quasi-classical trajectory (QCT) calculations for the exchange and abstraction processes for the H + HS reaction. These calculations were based on a newly constructed high-quality potential energy surface for the lowest triplet state of H(2)S ((3)A"). The ab initio single-point energies were computed using complete active space self-consistent field and multi-reference configuration interaction method with a basis set of aug-cc-pV5Z. The time-dependent wave packet (TDWP) method was used to calculate the total reaction probabilities and integral cross sections over the collision energy (E(col)) range of 0.0-2.0 eV for the reactant HS initially at the ground state and the first vibrationally excited state. It was found that the initial vibrational excitation of HS enhances both abstraction and exchange processes. In addition, a good agreement is found between QCT and TDWP reaction probabilities at the total momentum J = 0 as a function of collision energy for the H + HS (v = 0, j = 0) reaction.     

A computational investigation on the potential energy surface of thiosulfeno with O(3P) reaction

The reaction paths of thiosulfeno radical (HS₂) with O(³P) have been investigated at the UB3LYP/aug-cc-pV(T + d)Z and UCCSD(T)/aug-cc-pV(T + d)Z//B3LYP levels. Two stable collision intermediates, HSSO and SS(H)O, have been considered for the HS₂ + O(³P) reaction. Four products of S + HSO, H + SSO, HS + SO, and S₂ + OH are obtained by starting from HSSO and SS(H)O. The calculated results show that the most feasible paths for the formation of S + HSO, H + SSO, and HS + SO products include no transition states in reaction path, while that of S₂ + OH product includes relatively high energy barriers of 23.0 kcal/mol. Therefore, S + HSO, H + SSO, and HS + SO are main products (with the stability other of HS + SO > H + SSO > S + HSO) and S₂ + OH is the second product in HS₂ + O(³P) reaction. Because, all intermediates, transition states, and products involved in the reaction paths lie below the initial reactants, the HS₂ + O(³P) reaction is expected to be rapid even at low temperatures.