Theoretical study of the pyrolysis of methyltrichlorosilane in the gas phase. 2. Reaction paths and transition states

The kinetics for the previously proposed 114-reaction mechanism for the chemical vapor deposition (CVD) process that leads from methyltrichlorosilane (MTS) to silicon carbide (SiC) are examined. Among the 114 reactions, 41 are predicted to proceed with no intervening barrier. For the remaining 73 reactions, transition states and their corresponding barrier heights have been explored using second-order perturbation theory (MP2) with the aug-cc-pVDZ basis set. Final energies for the reaction barriers were obtained using both MP2 with the aug-cc-pVTZ basis set and coupled cluster theory (CCSD(T)) with the aug-cc-pVDZ basis set. CCSD(T)/aug-cc-pVTZ energies were estimated by assuming additivity of basis set and correlation effects. Partition functions for the computation of thermodynamic properties of the transition states were calculated with MP2/aug-cc-pVDZ. Forward and reverse Gibbs free energy barriers were obtained at 11 temperatures ranging from 0 to 2000 K. Important reaction pathways are illustrated at 0 and 1400 K.     

Basis set convergence of the coupled-cluster correction, δ(MP2)(CCSD(T)): best practices for benchmarking non-covalent interactions and the attendant revision of the S22, NBC10, HBC6, and HSG databases

In benchmark-quality studies of non-covalent interactions, it is common to estimate interaction energies at the complete basis set (CBS) coupled-cluster through perturbative triples [CCSD(T)] level of theory by adding to CBS second-order perturbation theory (MP2) a "coupled-cluster correction," δ(MP2)(CCSD(T)), evaluated in a modest basis set. This work illustrates that commonly used basis sets such as 6-31G*(0.25) can yield large, even wrongly signed, errors for δ(MP2)(CCSD(T)) that vary significantly by binding motif. Double-ζ basis sets show more reliable results when used with explicitly correlated methods to form a δ(MP2-F12)(CCSD(T(*))-F12) correction, yielding a mean absolute deviation of 0.11 kcal mol(-1) for the S22 test set. Examining the coupled-cluster correction for basis sets up to sextuple-ζ in quality reveals that δ(MP2)(CCSD(T)) converges monotonically only beyond a turning point at triple-ζ or quadruple-ζ quality. In consequence, CBS extrapolation of δ(MP2)(CCSD(T)) corrections before the turning point, generally CBS (aug-cc-pVDZ,aug-cc-pVTZ), are found to be unreliable and often inferior to aug-cc-pVTZ alone, especially for hydrogen-bonding systems. Using the findings of this paper, we revise some recent benchmarks for non-covalent interactions, namely the S22, NBC10, HBC6, and HSG test sets. The maximum differences in the revised benchmarks are 0.080, 0.060, 0.257, and 0.102 kcal mol(-1), respectively.     

Substituent effects in parallel-displaced pi-pi interactions

High-quality quantum-mechanical methods are used to examine how substituents tune pi-pi interactions between monosubstituted benzene dimers in parallel-displaced geometries. The present study focuses on the effect of the substituent across entire potential energy curves. Substituent effects are examined in terms of the fundamental components of the interaction (electrostatics, exchange-repulsion, dispersion and induction) through the use of symmetry-adapted perturbation theory. Both second-order M?ller-Plesset perturbation theory (MP2) with a truncated aug-cc-pVDZ' basis and spin-component-scaled MP2 (SCS-MP2) with the aug-cc-pVTZ basis are found to mimic closely estimates of coupled-cluster with perturbative triples [CCSD(T)] in an aug-cc-pVTZ basis. Substituents can have a significant effect on the electronic structure of the pi cloud of an aromatic ring, leading to marked changes in the pi-pi interaction. Moreover, there can also be significant direct interactions between a substituent on one ring and the pi-cloud of the other ring.     

An ab initio benchmark study of hydrogen bonded formamide dimers

The five singly and doubly hydrogen bonded dimers of formamide are calculated at the correlated level by using resolution of identity M?ller-Plesset second-order perturbation theory (RIMP2) and the coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] method. All structures are optimized with the Dunning aug-cc-pVTZ and aug-cc-pVQZ basis sets. The binding energies are extrapolated to the complete basis set (CBS) limit by using the aug-cc-pVXZ (X = D, T, Q) basis set series. The effect of extending the basis set to aug-cc-pV5Z on the geometries and binding energies is studied for the centrosymmetric doubly N-H...O bonded dimer FA1 and the doubly C-H...O bonded dimer FA5. The MP2 CBS limits range from -5.19 kcal/mol for FA5 to -14.80 kcal/mol for the FA1 dimer. The DeltaCCSD(T) corrections to the MP2 CBS limit binding energies calculated with the 6-31+G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets are mutually consistent to within < or =0.03 kcal/mol. The DeltaCCSD(T) correction increases the binding energy of the C-H...O bonded FA5 dimer by 0.4 kcal/mol or approximately 9% over the distance range +/-0.5 Angstrom relative to the potential minimum. This implies that the ubiquitous long-range C-H...O interactions in proteins are stronger than hitherto calculated.     

Stabilisation energy of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels

Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.     

Approximations to complete basis set-extrapolated, highly correlated non-covalent interaction energies

The first-principles calculation of non-covalent (particularly dispersion) interactions between molecules is a considerable challenge. In this work we studied the binding energies for ten small non-covalently bonded dimers with several combinations of correlation methods (MP2, coupled-cluster single double, coupled-cluster single double (triple) (CCSD(T))), correlation-consistent basis sets (aug-cc-pVXZ, X = D, T, Q), two-point complete basis set energy extrapolations, and counterpoise corrections. For this work, complete basis set results were estimated from averaged counterpoise and non-counterpoise-corrected CCSD(T) binding energies obtained from extrapolations with aug-cc-pVQZ and aug-cc-pVTZ basis sets. It is demonstrated that, in almost all cases, binding energies converge more rapidly to the basis set limit by averaging the counterpoise and non-counterpoise corrected values than by using either counterpoise or non-counterpoise methods alone. Examination of the effect of basis set size and electron correlation shows that the triples contribution to the CCSD(T) binding energies is fairly constant with the basis set size, with a slight underestimation with CCSD(T)∕aug-cc-pVDZ compared to the value at the (estimated) complete basis set limit, and that contributions to the binding energies obtained by MP2 generally overestimate the analogous CCSD(T) contributions. Taking these factors together, we conclude that the binding energies for non-covalently bonded systems can be accurately determined using a composite method that combines CCSD(T)∕aug-cc-pVDZ with energy corrections obtained using basis set extrapolated MP2 (utilizing aug-cc-pVQZ and aug-cc-pVTZ basis sets), if all of the components are obtained by averaging the counterpoise and non-counterpoise energies. With such an approach, binding energies for the set of ten dimers are predicted with a mean absolute deviation of 0.02 kcal/mol, a maximum absolute deviation of 0.05 kcal/mol, and a mean percent absolute deviation of only 1.7%, relative to the (estimated) complete basis set CCSD(T) results. Use of this composite approach to an additional set of eight dimers gave binding energies to within 1% of previously published high-level data. It is also shown that binding within parallel and parallel-crossed conformations of naphthalene dimer is predicted by the composite approach to be 9% greater than that previously reported in the literature. The ability of some recently developed dispersion-corrected density-functional theory methods to predict the binding energies of the set of ten small dimers was also examined.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

Geometries and stabilities of various configurations of benzene dimer: details of novel V-shaped structure revealed

Benzene dimer configurations namely T-shaped, parallel-displaced, sandwich, and V-shaped that were proposed by experimental studies are investigated using second- and fourth-order Møller–Plesset perturbation theory. The MP2 method with aug-cc-pVDZ and aug-cc-pVTZ basis sets unequivocally shows that the parallel-displaced configuration is considerably more stable than T-shaped structure. On the other hand, the MP4(SDTQ)/aug-cc-pVDZ level predicts that the T-shaped and parallel-displaced configurations are nearly isoenergetic, which is parallel to the previous results of estimated CCSD(T)/CBS level reported recently. The lowest energy T-shaped configuration is stabilized by 0.17 kcal/mol over the parallel-displaced configuration at the MP4(SDTQ)/aug-cc-pVDZ level. Although the structures of all the four different types of configurations are found to be stable at both MP2 and full MP4 methods, the V-shaped configuration is the least stable among them. The calculated interaction energy of ?2.3 kcal/mol for the lowest energy T-shaped structure at the MP4(SDTQ)/aug-cc-pVDZ level is in good agreement with the experimental value of ?2.4 ± 0.4 kcal/mol. We conclude that the MP4(SDTQ) with a reasonably good basis set can be used for systems involving π–π interactions to obtain qualitative and quantitative results.     

Ab initio MO study on the global minimum structure of C3H3 anion: propynl-1-yl versus allenyl anions

The geometrical structures of the C₃H⁻₃ anion are surveyed at the coupled-cluster doubles (CCD) level of theory with the aug-cc-pVDZ basis set. To clarify the CCD geometries, the stable two isomers -- propynl-l-yl 1 and allenyl 2 anions -- are further optimized at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory both with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. The final energies are calculated at the CCSD(T) and the complete active space self-consistent field (CASSCF) multi-reference internally contracted CI (MRCI) levels of theory with the aug-cc-pVTZ basis set. At the MRCI level of theory including both the corrections due to the cluster energies (MRCI+Q) and the zero-point vibrational energies, the allenyl anion 2 is about 1.3 kcal mol⁻¹ lower in energy than the propynl-l-yl anion 1. These results contrast with the previous theoretical estimates, where the propynl-l-yl anion 1 is 2-3 kcal mol⁻¹ lower in energy than the allenyl anion 2. The activation energies of the intramolecular hydrogen transfer in the 1 → 2 conversion reactions are 63.5 kcal mol⁻¹ at the MRCI+Q level of theory with the aug-cc-pVTZ basis set including the zero-point energy corrections. The adiabatic electron affinity of the planer propargyl (H₂CCCH) radical, which is the global minimum of the C₃H₃ radical, is calculated to be 0.976 eV (after correction for the zero-point energy changes) at the CCSD(T) level of theory with the aug-cc-pVTZ basis set. The present electron affinity is in fairly good agreement with the experimental one (0.893 eV) observed by Oakes and Ellison.     

Estimates of the ab initio limit for pi-pi interactions: the benzene dimer

State-of-the-art electronic structure methods have been applied to the simplest prototype of aromatic pi-pi interactions, the benzene dimer. By comparison to results with a large aug-cc-pVTZ basis set, we demonstrate that more modest basis sets such as aug-cc-pVDZ are sufficient for geometry optimizations of intermolecular parameters at the second-order M?ller-Plesset perturbation theory (MP2) level. However, basis sets even larger than aug-cc-pVTZ are important for accurate binding energies. The complete basis set MP2 binding energies, estimated by explicitly correlated MP2-R12/A techniques, are significantly larger in magnitude than previous estimates. When corrected for higher-order correlation effects via coupled cluster with singles, doubles, and perturbative triples [CCSD(T)], the binding energies D(e) (D(0)) for the sandwich, T-shaped, and parallel-displaced configurations are found to be 1.8 (2.0), 2.7 (2.4), and 2.8 (2.7) kcal mol(-1), respectively.     

Ab initio study of substituent effects in the interactions of dimethyl ether with aromatic rings

Ab initio calculations have been used to investigate the interaction energies of dimers of dimethyl ether with benzene, hexafluorobenzene, and several monosubstituted benzenes. The potential energy curves were explored at the MP2/aug-cc-pVDZ level for two basic configurations of the dimers, one in which the oxygen atom of the dimethyl ether was pointed towards the aromatic ring and the other in which the oxygen atom was pointed away from the aromatic ring. Once the optimum intermolecular distances between the dimethyl and the aromatic ring had been determined for each of the dimers in both configurations at the MP2/aug-cc-pVDZ level, single point energy calculations were performed at the MP2/aug-cc-pVTZ level. A CCSD(T) correction term to the energy was determined and this was combined with the MP2/aug-cc-pVTZ energies to estimate the CCSD(T)/aug-cc-pVTZ interaction energies of the dimers. The estimated CCSD(T)/aug-cc-pVTZ interaction energies are predicted to be attractive for all of the dimers in both configurations and dispersion interactions are found to be a large component of the stabilization of the dimers. For the dimers with the dimethyl ether oxygen pointing towards the aromatic ring, the strengths of interaction energies are found to increase as the aromatic ring becomes more electron deficient, while for the dimers with the dimethyl ether oxygen pointing away from the aromatic ring, they increase as the aromatic ring becomes more electron rich. In both cases, the trends can be explained in terms of the electrostatic potentials of the dimethyl ether and the aromatic rings.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

Symmetry breaking and electron correlation in O2−, O2, and O2+: A comparison of coupled cluster and quadratic configuration interaction approaches

Potential energy curves are calculated for O₂⁻, O₂, and O₂⁺ at the CCSD, QCISD, CCSD(T), and QCISD(T) levels of theory using aug-cc-pVDZ and aug-cc-pVTZ basis sets with electron correlation built onto inversion symmetry constrained and relaxed UHF wave functions. The spectroscopic constant r_e, w_e, w_e, x_e, D_j, and α_e, are determined from the potential curves using standard second-order perturbation theory expressions and are compared with experimental values to assess the relative accuracy of the theoretical approaches. Comparison of corresponding symmetry-constrained and symmetry-relaxed calculations indicates that the CCSD method is generally superior to CCSD(T), QCISD, and QCISD(T) in recovering from a symmetry-broken reference function. © 1996 John Wiley & Sons, Inc.     

State of the art theoretical study and comparison to experiment for the phenol...argon complex

The phenol...argon complex was studied by means of various high level ab initio quantum mechanics methods and high resolution threshold ionization spectroscopy. The structure and stabilization energy of different conformers were determined. Stabilization energy of van der Waals bonded and H-bonded PhOH...Ar complex determined at CCSD(T) complete basis set (CBS) level for CP-RI-MP2/cc-pVTZ/Ar aug-cc-pVTZ geometries amount to 434 and 285 cm(-1). The CCSD(T)/CBS were constructed either as a sum of MP2/CBS interaction energy and CCSD(T) correction term [difference between CCSD(T) and MP2 correlation energies determined with medium basis set] or directly from CCSD(T)/aug-cc-pVDZ and aug-cc-pVTZ energies. Both schemes provide very similar values. Harmonic vibrational analysis revealed that the H-bonded structure does not represent energy minimum but first order transition structure. The respective imaginary vibrational mode (16 cm(-1)) connects two possible argon locations -- above and below the phenol aromatic ring. Including the DeltaZPVE, we obtained stabilization enthalpy at 0 K of 389 cm(-1). This value is marginally higher (25-35 cm(-1), 0.07-0.10 kcal/mol) than the experimental value. The determination of DeltaZPVE constitutes the most significant error and possible improvements should come from more accurate evaluation of the (nonharmonic) vibrational frequencies.     

Intermolecular interactions in nitrogen-containing aromatic systems

π–π and CH···N interactions are vital in biological systems. In this study, stacking and hydrogen-bonded interactions in pyrazine and triazine dimers were investigated by density functional theory combined with symmetry-adapted perturbation theory (DFT-SAPT) and counterpoise (CP)-corrected supermolecular MP2, SCS-MP2, B3LYP-D and CCSD(T) calculations. All interaction energies were computed using the optimized structures at the CP-corrected SCS/aug-cc-pVDZ level, which gave 1–2 kJ/mol lower interaction energies than the ones computed at the MP2 level. For both dimers, doubly hydrogen-bonded and cross-(displaced) stacked orientations were found to be the lowest energy ones. The reference CCSD(T) calculations favored the former structure in both dimer systems, whereas MP2 and SCS-MP2 located the latter as the lowest energy isomer. In particular, the former was found to be lower in energy than the latter by 2.28 and 1.01 kJ/mol at the CCSD(T)/aug-cc-pVDZ level for pyrazine and triazine, respectively. B3LYP-D produced interaction energies in agreement with the CCSD(T) at the equilibrium geometries, but it overestimates them at the short range and underestimates at the long intermonomer separations. Furthermore, it tends to give smaller equilibrium distances compared to the CCSD(T). DFT-SAPT method was in a good agreement with the reference CCSD(T) calculations. This suggests that DFT-SAPT can be employed to compute the full potential energy surface of these dimers. Moreover, DFT-SAPT calculations showed that the electrostatic and dispersion contributions are the most important energy components stabilizing these dimers. The present study aims to show which theoretical method is the most promising one for the investigation of intermolecular interactions dominated by π–π and CH···N. Therefore, the findings obtained in this study can be used to unravel the structures of nucleic acid bases and other systems stabilized by π–π and CH···N interactions.     

Do N-heterocyclic aromatic rings prefer π-stacking?

The IR-UV double resonance spectroscopy of phenylacetylene complexes with triazine, pyrazine and pyridine in the acetylene C-H group of phenylacetylene was investigated. These spectra indicate that in the complexes of triazine, pyrazine and pyridine the acetylenic group is minimally perturbed and the predominant interaction is with the π electron density of the phenyl ring of phenylacetylene. Geometries of the complexes optimized at M06-2X/aug-cc-pVDZ and MP2/aug-cc-pVDZ levels, combined with highly accurate energy calculations at the complete basis set (CBS) limit of CCSD(T), indicate the formation of π-stacked complexes in all the three cases. Additionally, a C-H...N hydrogen-bonded complex between pyridine and phenylacetylene was also observed. The present results indicate that N-heterocyclic aromatic rings favor formation of π-stacked complexes.     

An ab initio theoretical prediction: an antiaromatic ring pi-dihydrogen bond accompanied by two secondary interactions in a "wheel with a pair of pedals" shaped complex FH . . . C4H4 . . . HF

By the counterpoise-correlated potential energy surface method (interaction energy optimization), the structure of the pi H-bond complex FH cdots, three dots, centered FH . . . C4H4 . . . HF has been obtained at the second-order M?ller-Plesset perturbation theory (MP2/aug-cc-pVDZ) level. Intermolecular interaction energy of the complex is calculated to be -7.8 kcal/mol at the coupled-cluster theory with single, double substitutions and perturbatively linked triple excitations CCSD (T)/aug-cc-pVDZ level. The optimized structure is a "wheel with a pair of pedals" shaped (1mid R:1) structure in which both HF molecules almost lie on either vertical line passing through the middle-point of the C[Double Bond]C bond on either side of the horizontal plane of the C4 ring for cyclobutadiene. In the structure, an antiaromatic ring pi-dihydrogen bond is found, in which the proton acceptor is antiaromatic 4 electron and 4 center pi bond and the donors are both acidic H atoms of HF molecules. In accompanying with the pi-dihydrogen bond, two secondary interactions are exposed. The first is a repulsive interaction between an H atom of HF and a near pair of H atoms of C4H4 ring. The second is the double pi-type H bond between two lone pairs on a F atom and a far pair of H atoms.     

Ab initio and analytic intermolecular potentials for Ar-CF4

Ab initio calculations at the CCSD(T) level of theory were performed to characterize the Ar + CF4 intermolecular potential. Potential energy curves were calculated with the aug-cc-pVTZ basis set, and with and without a correction for basis set superposition error (BSSE). Additional calculations were performed with other correlation consistent basis sets to extrapolate the Ar-CF4 potential energy minimum to the complete basis set (CBS) limit. Both the size of the basis set and BSSE have substantial effects on the Ar + CF4 potential. Calculations with the aug-cc-pVTZ basis set, and with a BSSE correction, appear to give a good representation of the BSSE corrected potential at the CBS limit. In addition, MP2 theory is found to give potential energies in very good agreement with those determined by the much higher level CCSD(T) theory. Two model analytic potential energy functions were determined for Ar + CF4. One is a fit to the aug-cc-pVTZ calculations with a BSSE correction. The second was derived by fitting an average BSSE corrected potential, which is an average of the CCSD(T)/aug-cc-pVTZ potentials with and without a BSSE correction. These analytic functions are written as a sum of two-body potentials and excellent fits to the ab initio potentials are obtained by representing each two-body interaction as a Buckingham potential.     

Kinetic isotope effects for Cl + CH4 ⇌ HCl + CH3 calculated using ab initio semiclassical transition state theory

Calculations were carried out for 25 isotopologues of the title reaction for various combinations of (35)Cl, (37)Cl, (12)C, (13)C, (14)C, H, and D. The computed rate constants are based on harmonic vibrational frequencies calculated at the CCSD(T)/aug-cc-pVTZ level of theory and X(ij) vibrational anharmonicity coefficients calculated at the CCSD(T) /aug-cc-pVDZ level of theory. For some reactions, anharmonicity coefficients were also computed at the CCSD(T)/aug-cc-pVTZ level of theory. The classical reaction barrier was taken from Eskola et al. [J. Phys. Chem. A 2008, 112, 7391-7401], who extrapolated CCSD(T) calculations to the complete basis set limit. Rate constants were calculated for temperatures from ～100 to ～2000 K. The computed ab initio rate constant for the normal isotopologue is in good agreement with experiments over the entire temperature range (～10% lower than the recommended experimental value at 298 K). The ab initio H/D kinetic isotope effects (KIEs) for CH(3)D, CH(2)D(2), CHD(3), and CD(4) are in very good agreement with literature experimental data. The ab initio (12)C/(13)C KIE is in error by ～2% at 298 K for calculations using X(ij) coefficients computed with the aug-cc-pVDZ basis set, but the error is reduced to ～1% when X(ij) coefficients computed with the larger aug-cc-pVTZ basis set are used. Systematic improvements appear to be possible. The present SCTST results are found to be more accurate than those from other theoretical calculations. Overall, this is a very promising method for computing ab initio kinetic isotope effects.