Ab initio and direct quasiclassical-trajectory study of the F+CH4-->HF+CH3 reaction

We present an electronic structure and dynamics study of the F+CH4-->HF+CH3 reaction. CCSD(T)/aug-cc-pVDZ geometry optimizations, harmonic-frequency, and energy calculations indicate that the potential-energy surface is remarkably isotropic near the transition state. In addition, while the saddle-point F-H-C angle is 180 degrees using MP2 methods, CCSD(T) geometry optimizations predict a bent transition state, with a 153 degrees F-H-C angle. We use these high-quality ab initio data to reparametrize the parameter-model 3 (PM3) semiempirical Hamiltonian so that calculations with the improved Hamiltonian and employing restricted open-shell wave functions agree with the higher accuracy data. Using this specific-reaction-parameter PM3 semiempirical Hamiltonian (SRP-PM3), we investigate the reaction dynamics by propagating quasiclassical trajectories. The results of our calculations using the SRP-PM3 Hamiltonian are compared with experiments and with the estimates of two recently reported potential-energy surfaces. The trajectory calculations using the SRP-PM3 Hamiltonian reproduce quantitatively the measured HF vibrational distributions. The calculations also agree with the experimental HF rotational distributions and capture the essential features of the excitation function. The results of the SRP semiempirical Hamiltonian developed here clearly improve over those using the two prior potential-energy surfaces and suggest that reparametrization of semiempirical Hamiltonians is a promising strategy to develop accurate potential-energy surfaces for reaction dynamics studies of polyatomic systems.     

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.     

Theoretical study of the dynamics of the H+CH4 and H+C2H6 reactions using a specific-reaction-parameter semiempirical Hamiltonian

We present a theoretical study of the reactions of hydrogen atoms with methane and ethane molecules and isotopomers. High-accuracy electronic-structure calculations have been carried out to characterize representative regions of the potential-energy surface (PES) of various reaction pathways, including H abstraction and H exchange. These ab initio calculations have been subsequently employed to derive an improved set of parameters for the modified symmetrically-orthogonalized intermediate neglect of differential overlap (MSINDO) semiempirical Hamiltonian, which are specific to the H+alkane family of reactions. The specific-reaction-parameter (SRP) Hamiltonian has then been used to perform a quasiclassical-trajectory study of both the H+CH4 and H+C2H6 reactions. The calculated values of dynamics properties of the H+CH4-->H2+CH3 reaction and isotopologues, including alkyl product speed distributions, diatomic product internal-state distributions, and cross sections, are generally in good agreement with experiment and with the results provided by the ZBB3 PES [Z. Xie et al., J. Chem. Phys. 125, 133120 (2006)]. The results of trajectories propagated with the SRP Hamiltonian for the H+C2H6-->H2+C2H5 reaction also agree with experiment. The level of agreement between the results calculated with the SRP Hamiltonian and experiment in both the H+methane and H+ethane reactions indicates that semiempirical Hamiltonians can be improved for not only a specific reaction but also a family of reactions.     

Classical trajectory study of the dynamics of the reaction of Cl atoms with ethane

We present an electronic-structure and dynamics study of the Cl + C2H6 --> HCl + C2H5 reaction. The stationary points of the ground-state potential energy surface have been characterized using various electronic-structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis limit, using geometries and harmonic frequencies obtained at the MP2/aug-cc-pVTZ level, are in agreement with the experimental reaction energy. Ab initio information has been used to reparameterize a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in key regions of the potential energy surface. The improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories. Computed kinetic energy release and scattering angle distributions at a collision energy of approximately 5.5 kcal mol(-1) are in reasonable agreement with experiments, but no evidence was found for the low translational energy HCl products scattered in the backward hemisphere reported in recent experiments.     

Potential energy surface, kinetics, and dynamics study of the Cl+CH4-->HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.     

Reduced dimensionality quantum dynamics of Cl + CH4 --> HCl + CH3 on an ab initio potential

We study the reaction Cl + CH(4)--> HCl + CH(3) using a 2-D potential energy surface obtained by fitting a double Morse analytical function to high level (CCSD(T)/cc-pVTZ//MP2/cc-pVTZ)ab initio data. Dynamics simulations are performed in hyperspherical coordinates with the close-coupled equations being solved using R-matrix propagation. Quantum contributions from spectator modes are included via a harmonic zero-point correction to the ab initio data prior to fitting the potential. This is the first time this method has been applied to a heavy-light-heavy reaction and the first time it has been used to study differential cross sections. We find thermal rate constants and state-to-state differential cross sections which are in good agreement with experimental data. We discuss the applicability of our method to the study of kinetic isotope effects (KIEs), which we derive for the CH(4)/CD(4) substitution. The calculated KIE compares favourably with experiment. Finally, we discuss the sensitivity of the results of dynamics simulations on the accuracy of the fitted potential.     

An ab initio based global potential energy surface describing CH5+ --> CH3+ + H2

A full-dimensional, ab initio based potential energy surface (PES) for CH(5)(+), which can describe dissociation is reported. The PES is a precise fit to 36173 coupled-cluster [CCSD(T)] calculations of electronic energies done using an aug-cc-pVTZ basis. The fit uses a polynomial basis that is invariant with respect to permutation of the five H atoms, and thus describes all 120 equivalent minima. The rms fitting error is 78.1 cm(-1) for the entire data set of energies up to 30,000 cm(-1) and a normal-mode analysis of CH(5)(+) also verifies the accuracy of the fit. Two saddle points have been located on the surface as well and compared with previous theoretical work. The PES dissociates correctly to the fragments CH(3)(+) + H(2) and the equilibrium geometry and normal-mode analyses of these fragments are also presented. Diffusion Monte Carlo calculations are done for the zero-point energies of CH(5)(+) (and some isotopologs) as well as for the separated fragments of CH(5)(+), CH(3)(+) + H(2) and those of CH(4)D(+), CH(3)(+) + HD and CH(2)D(+) + H(2). Values of D(0) are reported for these dissociations. A molecular dynamics calculation of CH(4)D(+) dissociation at one total energy is also performed to both validate the applicability of the PES for dynamics studies as well as to test a simple classical statistical prediction of the branching ratio of the dissociation products.     

Crossed beams and theoretical studies of the O(3P)+CH4-->H+OCH3 reaction excitation function

The excitation function for the reaction, O(3P)+CH4-->H+OCH3, has been measured in a crossed molecular beams experiment and determined with direct dynamics calculations that use the quasiclassical trajectory method in conjunction with a recently developed semiempirical Hamiltonian. Good agreement is found between experiment and theory, enabling us to address two fundamental issues for the O(3P)+CH4 reaction that arise for all O(3P)+saturated hydrocarbon reactions: (1) the importance of triplet excited states that correlate adiabatically to ground-state reactants and products and (2) the importance of intersystem crossing processes involving the lowest singlet surface [corresponding to reaction with O(1D)]. Our results indicate that the first excited triplet surface contributes substantially to the cross section when the collision energy exceeds the reaction barrier (approximately 2 eV) by more than 0.5 eV. Although triplet-singlet crossings may occur at all energies, we have found that their effect on the excitation function is negligible for the collision energies studied-up to 1.5 eV above threshold.     

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.     

An ab initio potential surface describing abstraction and exchange for H+CH4

We report an ab initio-based global potential energy surface for H+CH4 that describes the abstraction and exchange reactions. The PES, which is invariant with respect to any permutation of five H atoms, is a fit to 20,728 electronic energies calculated using the partially spin-restricted coupled-cluster method (RCCSD(T)) with a moderately large basis (aug-cc-pVTZ). A first set of quasiclassical trajectory calculations using this PES are reported for the H+CD4-->HD+CD3 reaction at collision energies of 0.65 and 1.52 eV and are compared to experiment and recent direct dynamics calculations done with density functional theory.     

Kinetics of CH + N2 revisited with multireference methods

The potential energy surface for the CH + N2 reaction was reexamined with multireference ab initio electronic structure methods employing basis sets up to aug-cc-pvqz. Comparisons with related CCSD(T) calculations were also made. The multireference ab initio calculations indicate significant shortcomings in single reference based methods for two key rate-limiting transition states. Transition state theory calculations incorporating the revised best estimates for the transition state properties provide order of magnitude changes in the predicted rate coefficient in the temperature range of importance to the mechanism for prompt NO formation. At higher temperatures, two distinct pathways make a significant contribution to the kinetics. A key part of the transition state analysis involves a variable reaction coordinate transition state theory treatment for the formation of H + NCN from HNCN. The present predictions for the rate coefficients resolve the discrepancy between prior theory and very recent experimental measurements.     

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.     

Combined valence bond-molecular mechanics potential-energy surface and direct dynamics study of rate constants and kinetic isotope effects for the H + C2H6 reaction

This article presents a multifaceted study of the reaction H+C(2)H(6)-->H(2)+C(2)H(5) and three of its deuterium-substituted isotopologs. First we present high-level electronic structure calculations by the W1, G3SX, MCG3-MPWB, CBS-APNO, and MC-QCISD/3 methods that lead to a best estimate of the barrier height of 11.8+/-0.5 kcal/mol. Then we obtain a specific reaction parameter for the MPW density functional in order that it reproduces the best estimate of the barrier height; this yields the MPW54 functional. The MPW54 functional, as well as the MPW60 functional that was previously parametrized for the H+CH(4) reaction, is used with canonical variational theory with small-curvature tunneling to calculate the rate constants for all four ethane reactions from 200 to 2000 K. The final MPW54 calculations are based on curvilinear-coordinate generalized-normal-mode analysis along the reaction path, and they include scaled frequencies and an anharmonic C-C bond torsion. They agree with experiment within 31% for 467-826 K except for a 38% deviation at 748 K; the results for the isotopologs are predictions since these rate constants have never been measured. The kinetic isotope effects (KIEs) are analyzed to reveal the contributions from subsets of vibrational partition functions and from tunneling, which conspire to yield a nonmonotonic temperature dependence for one of the KIEs. The stationary points and reaction-path potential of the MPW54 potential-energy surface are then used to parametrize a new kind of analytical potential-energy surface that combines a semiempirical valence bond formalism for the reactive part of the molecule with a standard molecular mechanics force field for the rest; this may be considered to be either an extension of molecular mechanics to treat a reactive potential-energy surface or a new kind of combined quantum-mechanical/molecular mechanical (QM/MM) method in which the QM part is semiempirical valence bond theory; that is, the new potential-energy surface is a combined valence bond molecular mechanics (CVBMM) surface. Rate constants calculated with the CVBMM surface agree with the MPW54 rate constants within 12% for 534-2000 K and within 23% for 200-491 K. The full CVBMM potential-energy surface is now available for use in variety of dynamics calculations, and it provides a prototype for developing CVBMM potential-energy surfaces for other reactions.     

Rate constants and hydrogen isotope substitution effects in the CH3 + HCl and CH3 + Cl2 reactions

The kinetics of the CH3 + Cl2 (k2a) and CD3 + Cl2 (k2b) reactions were studied over the temperature range 188-500 K using laser photolysis-photoionization mass spectrometry. The rate constants of these reactions are independent of the bath gas pressure within the experimental range, 0.6-5.1 Torr (He). The rate constants were fitted by the modified Arrhenius expression, k2a = 1.7 x 10(-13)(T/300 K)(2.52)exp(5520 J mol(-1)/RT) and k2b = 2.9 x 10(-13)(T/300 K)(1.84)exp(4770 J mol(-1)/RT) cm(3) molecule(-1) s(-1). The results for reaction 2a are in good agreement with the previous determinations performed at and above ambient temperature. Rate constants of the CH3 + Cl2 and CD3 + Cl2 reactions obtained in this work exhibit minima at about 270-300 K. The rate constants have positive temperature dependences above the minima, and negative below. Deuterium substitution increases the rate constant, in particular at low temperatures, where the effect reaches ca. 45% at 188 K. These observations are quantitatively rationalized in terms of stationary points on a potential energy surface based on QCISD/6-311G(d,p) geometries and frequencies, combined with CCSD(T) energies extrapolated to the complete basis set limit. 1D tunneling as well as the possibility of the negative energies of the transition state are incorporated into a transition state theory analysis, an approach which also accounts for prior experiments on the CH3 + HCl system and its various deuterated isotopic substitutions [Eskola, A. J.; Seetula, J. A.; Timonen, R. S. Chem. Phys. 2006, 331, 26].     

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.     

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.     

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.     

Theoretical study of isomerization and decomposition of propenal

We investigated the dynamics of isomerization and multi-channel dissociation of propenal (CH(2)CHCHO), methyl ketene (CH(3)CHCO), hydroxyl propadiene (CH(2)CH(2)CHOH), and hydroxyl cyclopropene (cyclic-C(3)H(3)-OH) in the ground potential-energy surface using quantum-chemical calculations. Optimized structures and vibrational frequencies of molecular species were computed with method B3LYP∕6-311G(d,p). Total energies of molecules at optimized structures were computed at the CCSD(T)∕6-311+G(3df,2p) level of theory. We established the potential-energy surface for decomposition to CH(2)CHCO + H, CH(2)CH + HCO, CH(2)CH(2)∕CH(3)CH + CO, CHCH∕CH(2)C + H(2)CO, CHCCHO∕CH(2)CCO + H(2), CHCH + CO + H(2), CH(3) + HCCO, CH(2)CCH + OH, and CH(2)CC∕cyclic-C(3)H(2) + H(2)O. Microcanonical rate coefficients of various reactions of trans-propenal with internal energies 148 and 182 kcal mol(-1) were calculated using Rice-Ramsperger-Kassel-Marcus and Variational transition state theories. Product branching ratios were derivable using numerical integration of kinetic master equations and the steady-state approximation. The concerted three-body dissociation of trans-propenal to fragments C(2)H(2) + CO + H(2) is the prevailing channel in present calculations. In contrast, C(3)H(3)O + H, C(2)H(3) + HCO and C(2)H(4) + CO were identified as major channels in the photolysis of trans-propenal. The discrepancy between calculations and experiments in product branching ratios indicates that the three major photodissociation channels occur mainly on an excited potential-energy surface whereas the other channels occur mainly on the ground potential-energy surface. This work provides profound insight in the mechanisms of isomerization and multichannel dissociation of the system C(3)H(4)O.     

Quasi-classical trajectory calculations analyzing the role of bending mode excitations of methane in the Cl + CH4 reaction

The effects of the methane torsional (nu(2)), umbrella (nu(4)), and the combination nu(2)+nu(4) bending mode excitations on the reactivity and dynamics of the gas-phase Cl + CH(4) --> HCl + CH(3) reaction were analyzed. Quasi-classical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were used on an analytical potential energy surface previously developed by our group. With respect to the reactivity, we found that excitation of either bending mode independently gave similar increases in the reactivity, while the increase observed upon excitation of both modes was larger than the sum of the effect of exciting them independently. Both results agree with recent experimental measures. With respect to the dynamics (rotovibrational and angular distributions of the products), the two bending modes and their combination gave very similar pictures, reproducing the experimental behavior. The satisfactory agreement obtained with a great variety of experimental data (always qualitatively acceptable and sometimes even quantitatively) of the present QCT study lends confidence to the potential energy surface constructed by our group.