Association rate constants for reactions between resonance-stabilized radicals: C3H3 + C3H3, C3H3 + C3H5, and C3H5 + C3H5

Reactions between resonance-stabilized radicals play an important role in combustion chemistry. The theoretical prediction of rate coefficients and product distributions for such reactions is complicated by the fact that the initial complex-formation steps and some dissociation steps are barrierless. In this paper direct variable reaction coordinate transition state theory (VRC-TST) is used to predict accurately the association rate constants for the self and cross reactions of propargyl and allyl radicals. For each reaction, a set of multifaceted dividing surfaces is used to account for the multiple possible addition channels. Because of their resonant nature the geometric relaxation of the radicals is important. Here, the effect of this relaxation is explicitly calculated with the UB3LYP/cc-pvdz method for each mutual orientation encountered in the configurational integrals over the transition state dividing surfaces. The final energies are obtained from CASPT2/cc-pvdz calculations with all pi-orbitals in the active space. Evaluations along the minimum energy path suggest that basis set corrections are negligible. The VRC-TST approach was also used to calculate the association rate constant and the corresponding number of states for the C(6)H(5) + H --> C(6)H(6) exit channel of the C(3)H(3) + C(3)H(3) reaction, which is also barrierless. For this reaction, the interaction energies were evaluated with the CASPT2(2e,2o)/cc-pvdz method and a 1-D correction is included on the basis of CAS+1+2+QC/aug-cc-pvtz calculations for the CH(3) + H reference system. For the C(3)H(3) + C(3)H(3) reaction, the VRC-TST results for the energy and angular momentum resolved numbers of states in the entrance channels and in the C(6)H(5) + H exit channel are incorporated in a master equation simulation to determine the temperature and pressure dependence of the phenomenological rate coefficients. The rate constants for the C(3)H(3) + C(3)H(3) and C(3)H(5) + C(3)H(5) self-reactions compare favorably with the available experimental data. To our knowledge there are no experimental rate data for the C(3)H(3) + C(3)H(5) reaction.     

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

A study of the hydrogen abstraction reactions of C2H radical with CH3CN, C2H5CN, and C3H7CN by dual-level generalized transition state theory

The hydrogen abstraction reactions C2H + CH3CN --> products (R1), C2H + CH3CH2CN --> products (R2), and C2H + CH3CH2CH2CN --> products (R3) have been investigated by dual-level generalized transition state theory. Optimized geometries and frequencies of all the stationary points and extra points along the minimum-energy path (MEP) are performed at the BH&H-LYP and MP2 methods with the 6-311G(d, p) basis set, and the energy profiles are further refined at the MC-QCISD level of theory. The rate constants are evaluated using canonical variational transition state theory (CVT) with a small-curvature tunneling correction (SCT) over a wide temperature range 104-2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values. Our calculations show that for reaction R2, the alpha-hydrogen abstraction channel and beta-hydrogen abstraction channel are competitive over the whole temperature range. For reaction R3, the gamma-hydrogen abstraction channel is preferred at lower temperatures, while the contribution of beta-hydrogen abstraction will become more significant with a temperature increase. The branching ratio to the alpha-hydrogen abstraction channel is found negligible over the whole temperature range.     

Theoretical study of the C6H3 potential energy surface and rate constants and product branching ratios of the C2H(2Sigma+) + C4H2(1Sigma(g)+) and C4H(2Sigma+) + C2H2(1Sigma(g)+) reactions

Ab initio CCSD(T)cc-pVTZ//B3LYP6-311G(**) and CCSD(T)/complete basis set (CBS) calculations of stationary points on the C(6)H(3) potential energy surface have been performed to investigate the reaction mechanism of C(2)H with diacetylene and C(4)H with acetylene. Totally, 25 different C(6)H(3) isomers and 40 transition states are located and all possible bimolecular decomposition products are also characterized. 1,2,3- and 1,2,4-tridehydrobenzene and H(2)CCCCCCH isomers are found to be the most stable thermodynamically residing 77.2, 75.1, and 75.7 kcal/mol lower in energy than C(2)H + C(4)H(2), respectively, at the CCSD(T)/CBS level of theory. The results show that the most favorable C(2)H + C(4)H(2) entrance channel is C(2)H addition to a terminal carbon of C(4)H(2) producing HCCCHCCCH, 70.2 kcal/mol below the reactants. This adduct loses a hydrogen atom from the nonterminal position to give the HCCCCCCH (triacetylene) product exothermic by 29.7 kcal/mol via an exit barrier of 5.3 kcal/mol. Based on Rice-Ramsperger-Kassel-Marcus calculations under single-collision conditions, triacetylene+H are concluded to be the only reaction products, with more than 98% of them formed directly from HCCCHCCCH. The C(2)H + C(4)H(2) reaction rate constants calculated by employing canonical variational transition state theory are found to be similar to those for the related C(2)H + C(2)H(2) reaction in the order of magnitude of 10(-10) cm(3) molecule(-1) s(-1) for T = 298-63 K, and to show a negative temperature dependence at low T. A general mechanism for the growth of polyyne chains involving C(2)H + H(C[triple bond]C)(n)H --> H(C[triple bond]C)(n+1)H + H reactions has been suggested based on a comparison of the reactions of ethynyl radical with acetylene and diacetylene. The C(4)H + C(2)H(2) reaction is also predicted to readily produce triacetylene + H via barrierless C(4)H addition to acetylene, followed by H elimination.     

Calculations on the competition between association and reaction for C3H(+) + H2

A potential energy surface has been calculated for the competing associative and reactive ion-molecule processes involving the reactants C3H(+) + H2. Our ab initio results show that the linear ion C3H+ and H2 can directly access the deep potential well of the propargyl ion H2CCCH+, which is calculated to lie 390 kJ mol-1 below the zero-point energy of the reactants. Isomerization between the propargyl ion and the lower energy, cyclic C3H3+ ion, calculated to lie 501 kJ mol-1 below the zero-point energy of reactants, can subsequently occur via two pathways. One of these pathways involves a transition state lying 22 kJ mol-1 below the energy of the reactants while the other, which occurs at much lower energies, involves two transition states and an intermediate. The dissociation of c-C3H3+ into c-C3H2(+) + H is calculated to occur directly, without any intermediate potential energy maximum, but the energy of the products lies 7.3 kJ mol-1 above the energy of the reactants. Using the minimum energy potential pathway and properties of the stationary point structures determined via ab initio methods, we have calculated both the association rate coefficient to produce C3H3+ as a function of density and the branching ratio between the propargyl and cyclic structures of the ion. Our results are in good agreement with some experimental results and in conflict with others. Specifically, we agree with the 1:1 branching ratio measured for the propargyl and cyclic isomers of C3H3+ at 80 and 300 K and we agree with the rate coefficient for radiative association measured at 80 K. We cannot reproduce reported measurements that the reactive channel (C3H2(+) + H) is the dominant channel at 80 K and at low gas densities, or that the association channel at high densities saturates at an effective rate coefficient well below the Langevin value -2x10(-11) cm3 s-1 at 300K and 1x10(-10) cm3 s-1 at 80K.     

New calculations on the ion-molecule processes C2H2(+) + H2 --> C2H3(+) + H and C2H2(+) + H2 --> C2H4+

New high-level quantum chemical calculations have been undertaken to understand the rates and mechanisms of the reactive and associative channels for the reactants C2H2(+) + H2. The reactive channel, which produces C2H3(+) + H, has been shown to be slightly endothermic, confirming earlier calculations at a somewhat lower level and in agreement with some recent experimental work. The associative channel, leading to C2H4+, has been shown to proceed via a transition state with negative energy relative to the reactants, so that association is predicted to be efficient. This result is in conflict with an earlier theoretical study but in agreement with low-temperature experimental measurements.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.     

Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H+NH3-->H2+NH2 reaction using a seven-dimensional model and an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated for the initial ground and seven excited states of NH3 with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants are calculated for the temperature range 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the total reaction probabilities are overall very small, (b) the symmetric and asymmetric NH stretch excitations enhance the reaction significantly and almost all of the excited energy deposited was used to reduce the reaction threshold, (c) the excitation of the umbrella and bending motion have a smaller contribution to the enhancement of reactivity, (d) the main contribution to the thermal rate constants is thought to come from the ground state at low temperatures and from the stretch excited states at high temperatures, and (e) the calculated thermal rate constants are three to ten times smaller than the experimental data and transition state theory results.     

Thermochemistry and accurate quantum reaction rate calculations for H2/HD/D2 + CH3

Accurate quantum-mechanical results for thermodynamic data, cumulative reaction probabilities (for J = 0), thermal rate constants, and kinetic isotope effects for the three isotopic reactions H2 + CH3 --> CH4 + H, HD + CH3 --> CH4 + D, and D2 + CH3 --> CH(3)D + D are presented. The calculations are performed using flux correlation functions and the multiconfigurational time-dependent Hartree (MCTDH) method to propagate wave packets employing a Shephard interpolated potential energy surface based on high-level ab initio calculations. The calculated exothermicity for the H2 + CH3 --> CH4 + H reaction agrees to within 0.2 kcal/mol with experimentally deduced values. For the H2 + CH3 --> CH4 + H and D2 + CH3 --> CH(3)D + D reactions, experimental rate constants from several groups are available. In comparing to these, we typically find agreement to within a factor of 2 or better. The kinetic isotope effect for the rate of the H2 + CH3 --> CH4 + H reaction compared to those for the HD + CH3 --> CH4 + D and D2 + CH3 --> CH(3)D + D reactions agree with experimental results to within 25% for all data points. Transition state theory is found to predict the kinetic isotope effect accurately when the mass of the transferred atom is unchanged. On the other hand, if the mass of the transferred atom differs between the isotopic reactions, transition state theory fails in the low-temperature regime (T < 400 K), due to the neglect of the tunneling effect.     

Chemically activated reactions on the C7H5 energy surface: propargyl + diacetylene, i-C5H3 + acetylene, and n-C5H3 + acetylene

This study uses computational chemistry and statistical reaction rate theory to investigate the chemically activated reaction of diacetylene (butadiyne, C(4)H(2)) with the propargyl radical (C˙H(2)CCH) and the reaction of acetylene (C(2)H(2)) with the i-C(5)H(3) (CH(2)CCCC˙H) and n-C(5)H(3) (CHCC˙HCCH) radicals. A detailed G3SX-level C(7)H(5) energy surface demonstrates that the C(3)H(3) + C(4)H(2) and C(5)H(3) + C(2)H(2) addition reactions proceed with moderate barriers, on the order of 10 to 15 kcal mol(-1), and form activated open-chain C(7)H(5) species that can isomerize to the fulvenallenyl radical with the highest barrier still significantly below the entrance channel energy. Higher-energy pathways are available leading to other C(7)H(5) isomers and to a number of C(7)H(4) species + H. Rate constants in the large multiple-well (15) multiple-channel (30) chemically activated system are obtained from a stochastic solution of the one-dimensional master equation, with RRKM theory for microcanonical rate constants. The dominant products of the C(4)H(2) + C(3)H(3) reaction at combustion-relevant temperatures and pressures are i-C(5)H(3) + C(2)H(2) and CH(2)CCHCCCCH + H, along with several quenched C(7)H(5) intermediate species below 1500 K. The major products in the n-C(5)H(3) + C(2)H(2) reaction are i-C(5)H(3) + C(2)H(2) and a number of C(7)H(4) species + H, with C(7)H(5) radical stabilization at lower temperatures. The i-C(5)H(3) + C(2)H(2) reaction predominantly leads to C(7)H(4) + H and to stabilized C(7)H(5) products. The title reactions may play an important role in polycyclic aromatic hydrocarbon (PAH) formation in combustion systems. The C(7)H(5) potential energy surface developed here also provides insight into several other important reacting gas-phase systems relevant to combustion and astrochemistry, including C(2)H + the C(3)H(4) isomers propyne and allene, benzyne + CH, benzene + C((3)P), and C(7)H(5) radical decomposition, for which some preliminary analysis is presented.     

Thermal decomposition of ethanol. III. A computational study of the kinetics and mechanism for the CH3+C2H5OH reaction

The mechanism for the CH3+C2H5OH reaction has been investigated by the modified Gaussian-2 method based on the geometric parameters of the stationary points optimized at the B3LYP/6-311+G(d,p) level of theory. Five transition states have been identified for the production of CH4+CH3CHOH (TS1), CH4+CH3CH2O (TS2), CH4+CH2CH2OH (TS3), CH3OH+CH3CH2 (TS4), and CH3CH2OCH3+H (TS5) with the corresponding barriers 12.0, 13.2, 16.0, 44.7, and 49.9 kcal/mol, respectively. The predicted rate constants and branching ratios for the three lower-energy H-abstraction reactions were calculated using the conventional and variational transition state theory with quantum-mechanical tunneling corrections for the temperature range 300-3000 K. The predicted total rate constant, kt=8.36 x 10(-76) T(20.00) exp(5258/T) cm3 mol(-1) s(-1) (300-600 K) and 6.10 x 10(-25) T(4.10)exp(-4058/T) cm3 mol(-1) s(-1) (600-3000 K), agrees closely with existing experimental data in the temperature range 403-523 K. Similarly, the predicted rate constants for CH3+CH3CD2OH and CD3+C2H5OD are also in reasonable agreement with available low temperature kinetic data.     

A systematic computational study of the reactions of HO2 with RO2: the HO2 + C2H5O2 reaction

The reaction of HO2 with C2H5O2 has been studied using the density functional theory (B3LYP) and the coupled-cluster theory [CCSD(T)]. The reaction proceeds on the triplet potential energy surface via hydrogen abstraction to form ethyl hydroperoxide and oxygen. On the singlet potential energy surface, the addition-elimination mechanism is revealed. Variational transition state theory is used to calculate the temperature-dependent rate constants in the range 200-1000 K. At low temperatures (e.g., below 300 K), the reaction takes place predominantly on the triplet surface. The calculated low-temperature rate constants are in good agreement with the experimental data. As the temperature increases, the singlet reaction mechanism plays more and more important role, with the formation of OH radical predominantly. The isotope effect of the reaction (DO2 + C2D5O2 vs HO2 + C2H5O2) is negligible. In addition, the triplet abstraction energetic routes for the reactions of HO2 with 11 alkylperoxy radicals (CnHmO2) are studied. It is shown that the room-temperature rate constants have good linear correlation with the activation energies for the hydrogen abstraction.     

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.     

Theoretical study on the rate constants for the C2H5 + HBr --> C2H6 + Br reaction

The reaction C(2)H(5) + HBr --> C(2)H(6) + Br has been theoretically studied over the temperature range from 200 to 1400 K. The electronic structure information is calculated at the BHLYP/6-311+G(d,p) and QCISD/6-31+G(d) levels. With the aid of intrinsic reaction coordinate theory, the minimum energy paths (MEPs) are obtained at the both levels, and the energies along the MEP are further refined by performing the single-point calculations at the PMP4(SDTQ)/6-311+G(3df,2p)//BHLYP and QCISD(T)/6-311++G(2df,2pd)//QCISD levels. The calculated ICVT/SCT rate constants are in good agreement with available experimental values, and the calculate results further indicate that the C(2)H(5) + HBr reaction has negative temperature dependence at T < 850 K, but clearly shows positive temperature dependence at T > 850 K. The current work predicts that the kinetic isotope effect for the title reaction is inverse in the temperature range from 200 to 482 K, i.e., k(HBr)/k(DBr) < 1.     

Computational study on the kinetics and mechanism for the unimolecular decomposition of C6H5NO2 and the related C6H5 + NO2 and C6H5O + NO reactions

The kinetics and mechanisms for the unimolecular dissociation of nitrobenzene and related association reactions C(6)H(5) + NO(2) and C(6)H(5)O + NO have been studied computationally at the G2M(RCC, MP2) level of theory in conjunction with rate constant prediction with multichannel RRKM calculations. Formation of C(6)H(5) + NO(2) was found to be dominant above 850 K with its branching ratio > 0.78, whereas the formation of C(6)H(5)O + NO via the C(6)H(5)ONO intermediate was found to be competitive at lower temperatures, with its branching ratio increasing from 0.22 at 850 K to 0.97 at 500 K. The third energetically accessible channel producing C(6)H(4) + HONO was found to be uncompetitive throughout the temperature range investigated, 500-2000 K. The predicted rate constants for C(6)H(5)NO(2) --> C(6)H(5) + NO(2) and C(6)H(5)O + NO --> C(6)H(5)ONO under varying experimental conditions were found to be in good agreement with all existing experimental data. For C(6)H(5) + NO(2), the combination processes producing C(6)H(5)ONO and C(6)H(5)NO(2) are dominant at low temperature and high pressure, while the disproportionation process giving C(6)H(5)O + NO via C(6)H(5)ONO becomes competitive at low pressure and dominant at temperatures above 1000 K.     

Full dimensional time-dependent quantum dynamics study of the H + NH3 --> H2 + NH2 reaction

A rigorous full dimensional time-dependent wave packet method has been developed for the reactive scattering between an atom and a tetra-atomic molecule. The method has been applied to the hydrogen abstraction reaction H+NH(3)-->H(2)+NH(2). Initial state-selected total reaction probabilities are investigated for the reactions from the ground vibrational state and from four excited vibrational states of ammonia. The total reaction probabilities from two lowest "tunneling doublets" due to the inversion barrier for the umbrella bending motion of NH(3) and from two pairs of doubly degenerate vibrational states of NH(3) are also inspected. Integral cross sections and rate constants are calculated for the reaction from the ground state with the centrifugal-sudden approximation. The calculated results are compared with those from the previous seven dimensional calculations [M. Yang and J. C. Corchado, J. Chem. Phys. 126, 214312 (2007)]. This work shows that the full dimensional rate constants are a factor of 3 larger than the corresponding seven dimensional calculated values at T=200 K and are overall smaller than those obtained from the variational transition state theory in the whole temperature region. The work also reveals that nonreactive NH bonds of NH(3) cannot be treated as spectators due to the fact that three NH bonds are coupled with each other during the reaction process.     

A Computational Study of the Kinetics and Mechanism for the C2H3 + CH3OH Reaction

The mechanism for the C₂H₃ + CH₃OH reaction has been investigated by the Gaussian‐4 (G4) method based on the geometric parameters of the stationary points optimized at the B3LYP/6–31G(2df, p) level of theory. Four transition states have been identified for the production of C₂H₄ + CH₃O (TSR/P1), C₂H₄ + CH₂OH (TSR/P2), C₂H₃OH + CH₃ (TSR/P3), and C₂H₃OCH₃ + H (TSR/P4) with the corresponding barriers 8.48, 9.25, 37.62, and 34.95 kcal/mol at the G4 level of theory, respectively. The rate constants and branching ratios for the two lower energy H‐abstraction reactions were calculated using canonical variational transition state theory with the Eckart tunneling correction at the temperature range 300–2500 K. The predicted rate constants have been compared with existing literature data, and the uncertainty has been discussed. The branching ratio calculation suggests that the channel producing CH₃O is dominant up to about 1070 K, above which the channel producing CH₂OH becomes very competitive.     

HNCS与CH2(X2Π)反应微观动力学的理论研究

用量子化学密度泛函理论的UB3LYP/6-311+G**方法和高级电子相关的UQCISD(T)/6-311+G**方法研究了异硫氰酸(HNCS)与乙炔基自由基(C2H(X2Π))反应的微观机理. 采用双水平直接动力学方法IVTST-M, 获取反应的势能面信息, 应用正则变分过渡态理论并考虑小曲率隧道效应, 计算了在250～2500 K温度范围内反应的速率常数. 研究结果表明, HNCS与C2H(X2Π)反应为多通道、多步骤的复杂反应, 共存在三个可能的反应通道, 主反应通道为通过分子间H原子迁移, 生成主要产物NCS+C2H2. 反应速率常数随温度升高而增大, 表现为正温度效应. 速率常数计算中变分效果很小. 在低温区隧道效应对反应速率的贡献较大, 反应为放热反应.     

Ab initio studies of alkyl radical reactions: Combination and disproportionation reactions of CH3 with C2H5, and the decomposition of chemically activated C3H8

This paper reports the first quantitative ab initio prediction of the disproportionation/combination ratio of alkyl+alkyl reactions using CH3+C2H5 as an example. The reaction has been investigated by the modified Gaussian-2 method with variational transition state or Rice-Ramsperger-Kassel-Marcus calculations for several channels producing (1) CH4+CH2CH2, (2) C3H8, (3) CH4CH3CH, (4) H2+CH3CHCH2, (5) H2+CH3CCH3, and (6) C2H6+CH2 by H-abstraction and association/decomposition mechanisms through singlet and triplet potential energy paths. Significantly, the disproportionation reaction (1) producing CH4+C2H4 was found to occur primarily by the lowest energy path via a loose hydrogen-bonding singlet molecular complex, H3CHC2H4, with a 3.5 kcal/mol binding energy and a small decomposition barrier (1.9 kcal/mol), instead of a direct H-abstraction process. Bimolecular reaction rate constants for the formation of the above products have been calculated in the temperature range 300-3000 K. At 1 atm, formation of C3H8 is dominant below 1200 K. Over 1200 K, the disproportionation reaction becomes competitive. The sum of products (3)-(6) accounts for less than 0.3% below 1500 K and it reaches around 1%-4% above 2000 K. The predicted rate constant for the disproportionation reaction with multiple reflections above the complex well, k1=5.04 x T(0.41) exp(429/T) at 200-600 K and k1=1.96 x 10(-20) T(2.45) exp(1470/T) cm3 molecule(-1) s(-1) at 600-3000 K, agrees closely with experimental values. Similarly, the predicted high-pressure rate constants for the combination reaction forming C3H8 and its reverse dissociation reaction in the temperature range 300-3000 K, k2(infinity)=2.41 x 10(-10) T(-0.34) exp(259/T) cm3 molecule(-1) s(-1) and k(-2)(infinity)=8.89 x 10(22) T(-1.67)exp(-46 037/T) s(-1), respectively, are also in good agreement with available experimental data.     

Ab initio study on the kinetics of hydrogen abstraction for the H+alkene-->H2+alkenyl reaction class

Kinetics of the hydrogen abstraction reaction class of the H+alkene has been studied using the reaction class transition state theory (RC-TST) combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approach. The rate constants for the reference reaction, H+C2H4, were obtained by the canonical variational transition state theory (CVT) with the small curvature tunneling (SCT) correction in the temperature range of 300-3000 K. Combined with these data, both the RC-TST/LER, where only reaction energy is needed, and RC-TST/BHG, where no other information is needed, are found to be promising methods for predicting rate constants for a large number of reactions in this reaction class. Our analysis indicates that less than 50% systematic errors on the average exist in the predicted rate constants using the RC-TST/LER or RC-TST/BHG method while in comparison to explicit rate calculations the differences are less than 100% or a factor of 2 on the average.