Kinetics of the unimolecular decomposition of the 2-chloroallyl radical

The thermal decomposition of the 2-chloroallyl radical, CH(2)CClCH(2) --> CH(2)CCH(2) + Cl (1), was studied using the laser photolysis/photoionization mass spectrometry technique. Rate constants were determined in time-resolved experiments as a function of temperature (720-840 K) and bath gas density ([He] = (3-12) x 10(16), [N(2)] = 6 x 10(16) molecule cm(-3)). C(3)H(4) was observed as a primary product of reaction 1. The rate constants of reaction 1 are in the falloff, close to the low-pressure limit, under the conditions of the experiments. The potential energy surface (PES) of reaction 1 was studied using a variety of quantum chemical methods. The results of the study indicate that the minimum energy path of the CH(2)CClCH(2) dissociation proceeds through a PES plateau corresponding to a weakly bound Cl-C(3)H(4) complex; a PES saddle point exists between the equilibrium CH(2)CClCH(2) structure and the Cl-C(3)H(4) complex. The results of quantum chemical calculations, the rate constant values obtained in the experimental study, and literature data on the reverse reaction of addition of Cl to allene were used to create a model of reactions 1 and -1. The experimental dependences of the rate constants on temperature and pressure were reproduced in RRKM/master equation calculations. The reaction model provides expressions for the temperature dependences of the high-pressure-limit and the low-pressure-limit rate constants and the falloff broadening factors (at T = 300-1600 K): k(infinity)(1) = 1.45 x 10(20)T(-1.75) exp(-19609 K/T) s(-1), k(infinity)(-)(1) = 8.94 x 10(-10)T(-0.40) exp(481 K/T) cm(3) molecule(-1) s(-1), k(1)(0)(He) = 5.01 x 10(-32)T(-12.02) exp(-22788 K/T) cm(3) molecule(-1) s(-1), k(1)(0)(N(2)) = 2.50 x 10(-32)T(-11.92) exp(-22756 K/T) cm(3) molecule(-1) s(-1), F(cent)(He) = 0.46 exp(-T/1001 K) + 0.54 exp(-T/996 K) + exp(-4008 K/T), and F(cent)(N(2)) = 0.37 exp(-T/2017 K) + 0.63 exp(-T/142 K) + exp(-4812 K/T). The experimental data are not sufficient to specify all the parameters of the model; consequently, some of the model parameters were obtained from quantum chemical calculations and from analogy with other reactions of radical decomposition. Thus, the parametrization is most reliable under conditions close to those used in the experiments.     

Low-energy paths for the unimolecular decomposition of CH3OH: a G2M/statistical theory study

The potential energy surface (PES) of the CH3OH system has been characterized by ab initio molecular orbital theory calculations at the G2M level of theory. The mechanisms for the decomposition of CH3OH and the related bimolecular reactions, CH3 + OH and 1CH2 + H2O, have been elucidated. The rate constants for these processes have been calculated using variational RRKM theory and compared with available experimental data. The total decomposition rate constants of CH3OH at the high- and low-pressure limits can be represented by k infinity = 1.56 x 10(16) exp(-44,310/T) s-1 and kAr0 = 1.60 x 10(36) T-12.2 exp(-48,140/T) cm3 molecule-1 s-1, respectively, covering the temperature range 1000-3000 K, in reasonable agreement with the experimental values. Our results indicate that the product branching ratios are strongly pressure dependent, with the production of CH3 + OH and 1CH2 + H2O dominant under high (P > 10(3) Torr) and low (P < 1 atm) pressures, respectively. For the bimolecular reaction of CH3 and OH, the total rate constant and the yields of 1CH2 + H2O and H2 + HCOH at lower pressures (P < 5 Torr) could be reasonably accounted for by the theory. For the reaction of 1CH2 with H2O, both the yield of CH3 + OH and the total rate constant could also be satisfactorily predicted theoretically. The production of 3CH2 + H2O by the singlet to triplet surface crossing, predicted to occur at 4.3 kcal mol-1 above the H2C...OH2 van der Waals complex (which lies 82.7 kcal mol-1 above CH3OH), was neglected in our calculations.     

Experimental and computational studies of the phenyl radical reaction with propyne.

The kinetics for the gas-phase reaction of phenyl radical with propyne has been measured by cavity ring-down spectrometry (CRDS), and the mechanism and initial product branching have been elucidated with the help of quantum chemical calculations. Absolute rate constants measured by the CRDS technique can be expressed by the following Arrhenius equation: (k/cm(3) mol(-1) s(-1)): k(propyne)(T=301-428 K)=(3.68+/-0.92) x 10(11)exp[-(1685+/-80)/T]. The experiment is unable to distinguish between the possible reactive channels, but theory indicates that phenyl radicals preferably add to the unsaturated terminal carbon atom in propyne under our experimental conditions. Theoretical kinetic calculations, employing high-level G2M(RCC, RMP2) and G3 energetic and IRCMax(RCCSD(T)//B3LYP-DFT) molecular parameters, reproduce the total experimental rate constants within a factor of three. Calculated total and branching rate constants are provided for high-T kinetic modeling. Addition reactions of phenyl to C3H4 are estimated to be less important molecular-growth pathways in high-T conditions (T>1000 K) in comparison to the C6H5 + C2H2 reaction.     

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.     

Variational analysis of the phenyl + O2 and phenoxy + O reactions

Variational transition state analysis was performed on the barrierless phenyl + O2 and phenoxy + O association reactions. In addition, we also calculated rate constants for the related vinyl radical (C2H3) + O2 and vinoxy radical (C2H3O) + O reactions and provided rate constant estimates for analogous reactions in substituted aromatic systems. Potential energy scans along the dissociating C-OO and CO-O bonds (with consideration of C-OO internal rotation) were obtained at the O3LYP/6-31G(d) density functional theory level. The CO-O and C-OO bond scission reactions were observed to be barrierless, in both phenyl and vinyl systems. Potential energy wells were scaled by G3B3 reaction enthalpies to obtain accurate activation enthalpies. Frequency calculations were performed for all reactants and products and at points along the potential energy surfaces, allowing us to evaluate thermochemical properties as a function of temperature according to the principles of statistical mechanics and the rigid rotor harmonic oscillator (RRHO) approximation. The low-frequency vibrational modes corresponding to R-OO internal rotation were omitted from the RRHO analysis and replaced with a hindered internal rotor analysis using O3LYP/6-31G(d) rotor potentials. Rate constants were calculated as a function of temperature (300-2000 K) and position from activation entropies and enthalpies, according to canonical transition state theory; these rate constants were minimized with respect to position to obtain variational rate constants as a function of temperature. For the phenyl + O2 reaction, we identified the transition state to be located at a C-OO bond length of between 2.56 and 2.16 A (300-2000 K), while for the phenoxy + O reaction, the transition state was located at a CO-O bond length of 2.00-1.90 A. Variational rate constants were fit to a three-parameter form of the Arrhenius equation, and for the phenyl + O2 association reaction, we found k(T) = 1.860 x 1013T-0.217 exp(0.358/T) (with k in cm3 mol-1 s-1 and T in K); this rate equation provides good agreement with low-temperature experimental measurements of the phenyl + O2 rate constant. Preliminary results were presented for a correlation between activation energy (or reaction enthalpy) and pre-exponential factor for heterolytic O-O bond scission reactions.     

Experimental and modeling studies of the pressure and temperature dependences of the kinetics and the OH yields in the acetyl + O2 reaction

The acetyl + O(2) reaction has been studied by observing the time dependence of OH by laser-induced fluorescence (LIF) and by electronic structure/master equation analysis. The experimental OH time profiles were analyzed to obtain the kinetics of the acetyl + O(2) reaction and the relative OH yields over the temperature range of 213-500 K in helium at pressures in the range of 5-600 Torr. More limited measurements were made in N(2) and for CD(3)CO + O(2). The relative OH yields were converted into absolute yields by assuming that the OH yield at zero pressure is unity. Electronic structure calculations of the stationary points of the potential energy surface were used with a master equation analysis to fit the experimental data in He using the high-pressure limiting rate coefficient for the reaction, k(∞)(T), and the energy transfer parameter, (ΔE(d)), as variable parameters. The best-fit parameters obtained are k(∞) = 6.2 × 10(-12) cm(-3) molecule(-1) s(-1), independent of temperature over the experimental range, and (ΔE(d))(He) = 160(T/298?K) cm(-1). The fits in N(2), using the same k(∞)(T), gave (ΔE(d))(N(2)) = 270(T/298?K) cm(-1). The rate coefficients for formation of OH and CH(3)C(O)O(2) are provided in parametrized form, based on modified Troe expressions, from the best-fit master equation calculations, over the pressure and temperature ranges of 1 ≤ p/Torr ≤ 1.5 × 10(5) and 200 ≤ T/K ≤ 1000 for He and N(2) as the bath gas. The minor channels, leading to HO(2) + CH(2)CO and CH(2)C(O)OOH, generally have yields <1% over this range.     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Kinetics of the CH2Cl + CH3 and CHCl2 + CH3 radical-radical reactions

The CH2Cl + CH3 (1) and CHCl2 + CH3 (2) cross-radical reactions were studied by laser photolysis/photoionization mass spectroscopy. Overall rate constants were obtained in direct real-time experiments in the temperature region 301-800 K and bath gas (helium) density (6-12) x 10(16) atom cm(-3). The observed rate constant of reaction 1 can be represented by an Arrhenius expression k1 = 3.93 x 10(-11) exp(91 K/T) cm3 molecule(-1) s(-1) (+/-25%) or as an average temperature-independent value of k1= (4.8 +/- 0.7) x 10(-11) cm3 molecule(-1) s(-1). The rate constant of reaction 2 can be expressed as k2= 1.66 x 10(-11) exp(359 K/T) cm3 molecule(-1) s(-1) (+/-25%). C2H4 and C2H3Cl were detected as the primary products of reactions 1 and 2, respectively. The experimental values of the rate constant are in reasonable agreement with the prediction based on the "geometric mean rule." A separate experimental attempt to determine the rate constants of the high-temperature CH2Cl + O2 (10) and CHCl2 + O2 (11) reaction resulted in an upper limit of 1.2 x 10(-16) cm(3) molecule(-1) s(-1) for k10 and k11 at 800 K.     

A computational study on the kinetics and mechanism for the unimolecular decomposition of o-nitrotoluene

The kinetics and mechanism for the unimolecular decomposition of o-nitrotoluene (o-CH(3)C(6)H(4)NO(2)) have been studied computationally at the G2M(RCC, MP2)//B3LYP/6-311G(d, p) level of theory in conjunction with rate constant predictions with RRKM and TST calculations. The results of the calculations reveal 10 decomposition channels for o-nitrotoluene and its six isomeric intermediates, among them four channels give major products: CH(3)C(6)H(4) + NO(2), C(6)H(4)C(H)ON (anthranil) + H(2)O, CH(3)C(6)H(4)O (o-methyl phenoxy) + NO, and C(6)H(4)C(H(2))NO + OH. The predicted rate constants in the 500-2000 K temperature range indicate that anthranil production, taking place initially by intramolecular H-abstraction from the CH(3) group by NO(2) followed by five-membered ring formation and dehydration, dominates at temperatures below 1000 K, whereas NO(2) elimination becomes predominant above 1100 K and CH(3)C(6)H(4)O formation by the nitro-nitrite isomerization/decomposition process accounts for only 5-11% of the total product yield in the middle temperature range 800-1300 K. The branching ratio for CH(2)C(6)H(4)NO formation by the decomposition process of CH(2)C(6)H(4)N(O)OH is negligible. The predicted high-pressure-limit rate constants with the rate expression of 4.10 x 10(17) exp[-37000/T] s(-1) for the NO(2) elimination channel and 9.09 x 10(12) exp[-25800/T] s(-1) for the H(2)O elimination channel generally agree reasonably with available experimental data. The predicted high-pressure-limit rate constants for the NO and OH elimination channels are represented as 1.49 x 10(14) exp[-30000/T] and 1.31 x 10(15) exp[-38000/T] s(-1), respectively.     

Theoretical study of the reaction of ethane with oxygen molecules in the ground triplet and singlet delta States

Quantum chemical calculations are carried out to study the reaction of ethane with molecular oxygen in the ground triplet and singlet delta states. Transition states, intermediates, and possible products of the reaction on the triplet and singlet potential energy surfaces are identified on the basis of the coupled-cluster method. The basis set dependence of coupled-cluster energy values is estimated by the second-order perturbation theory. The values of energy barriers are also refined by using the compound CBS-Q and G3 techniques. It was found that the C(2)H(6) + O(2)(X(3)Σ(g)(-)) reaction leads to the formation of C(2)H(5) and HO(2) products, whereas the C(2)H(6) + O(2)(a(1)Δ(g)) process produces C(2)H(4) and H(2)O(2) molecules. The appropriate rate constants of these reaction paths are estimated on the basis of variational and nonvariational transition-state theories assuming tunneling and possible nonadiabatic transitions in the temperature range 500-4000 K. The calculations showed that the rate constant of the C(2)H(6) + O(2)(a(1)Δ(g)) reaction path is much greater than that of the C(2)H(6) + O(2)(X(3)Σ(g)(-)) one. At the same time, the singlet and triplet potential surface intersection is detected that leads to the appearance of the nonadiabatic quenching channel O(2)(a(1)Δ(g)) + C(2)H(6) → O(2)(X (3)Σ(g)(-)) + C(2)H(6). The rate constant of this process is estimated with the use of the Landau-Zener model. It is demonstrated that, in the case of the existence of thermal equilibrium in the distribution of molecules over the electronic states, at low temperatures (T < 1200 K) the main products of the reaction of C(2)H(6) with O(2) are C(2)H(4) and H(2)O(2), rather than C(2)H(5) and HO(2). At higher temperature (T > 1200 K) the situation is inverted.     

Experimental and theoretical studies of the kinetics of the reactions of OH and OD with 2-methyl-3-buten-2-ol between 300 and 415 K at low pressure

The rate constants for the reactions of OH and OD with 2-methyl-3-buten-2-ol (MBO) have been measured at 2, 3, and 5 Torr total pressure over the temperature range 300-415 K using a discharge-flow system coupled with laser induced fluorescence detection of OH. The measured rate constants at room temperature and 5 Torr for the OH + MBO reaction in the presence of O2 and the OD + MBO reaction are (6.32 +/- 0.27) and (6.61 +/- 0.66) x 10(-11) cm3 molecule(-1) s(-1), respectively, in agreement with previous measurements at higher pressures. However, the rate constants begin to show a pressure dependence at temperatures above 335 K. An Arrhenius expression of k0 = (2.5 +/- 7.4) x 10(-32) exp[(4150 +/- 1150)/T] cm6 molecule(-2) s(-1) was obtained for the low-pressure-limiting rate constant for the OH + MBO reaction in the presence of oxygen. Theoretical calculations of the energetics of the OH + MBO reaction suggest that the stability of the different HO-MBO adducts are similar, with predicted stabilization energies between 27.0 and 33.4 kcal mol(-1) relative to the reactants, with OH addition to the internal carbon predicted to be 1-4 kcal mol(-1) more stable than addition to the terminal carbon. These stabilization energies result in estimated termolecular rate constants for the OH + MBO reaction using simplified calculations based on RRKM theory that are in reasonable agreement with the experimental values.     

The association reaction between C2H and 1-butyne: a computational chemical kinetics study

The potential energy surfaces (PES) for the reaction of the C(2)H radical with 1-butyne (C(4)H(6)) have been studied using the CBS-QB3 method. Density functional B3LYP/cc-pVTZ and M06-2X/6-311++G(d,p) calculations have also been performed to analyze the reaction energetics. For detailed theoretical calculation on the total reaction mechanism, the initial association reactions on more and less substituted C atoms of 1-butyne are treated separately followed by a variational transition state theory (VTST) calculation to obtain reaction rates. The successive unimolecular reactions from the association reaction complexes are subjected to Rice-Ramsperger-Kassel-Marcus (RRKM) calculations for reaction rate constants and product branching ratios. The calculated rate constants in the temperature range 70-295 K for both the association reactions are found to be highly temperature dependent at low temperatures, which is contrary to the experimental findings of temperature independent association rates. We have explained this observation with the help of variational nature of the transition states, and we found a "loose" transition state at low temperatures. The calculated product branching ratios for the unimolecular reactions generally agree with the available experimental data, although some channels show a significant method dependency and therefore the correlation with experiment is lost to some extent. Our detailed reaction energetics calculations confirm that the C(2)H + C(4)H(6) reaction proceeds without an entrance barrier and leads to the important products ethynylallene + CH(3), 1,3-hexadiyne + H, 3,4-hexadiene-1-yne + H, 2-ethynyl-1,3-butadiene + H, 3,4-dimethylenecyclobut-1-ene + H and fulvene + H exothermic by 25-75 kcal mol(-1), with strong dependence of the product distribution on the association mode of C(2)H with C(4)H(6), making these reactions fast under low temperature conditions of Titan's atmosphere. Therefore this study can provide a detailed picture of the complex hydrocarbon formation mechanism in the upper atmosphere.     

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.     

CH(A2Delta) Formation in hydrocarbon combustion: the temperature dependence of the rate constant of the reaction C2H + O2--> CH(A2Delta) + CO2

The temperature dependence of the rate constant of the chemiluminescence reaction C2H + O2 --> CH(A) + CO2, k1e, has been experimentally determined over the temperature range 316-837 K using pulsed laser photolysis techniques. The rate constant was found to have a pronounced positive temperature dependence given by k1e(T) = AT(4.4) exp(1150 +/- 150/T), where A = 1 x 10(-27) cm(3) s(-1). The preexponential factor for k1e, A, which is known only to within an order of magnitude, is based on a revised expression for the rate constant for the C2H + O(3P) --> CH(A) + CO reaction, k2b, of (1.0 +/- 0.5) x 10(-11) exp(-230 K/T) cm3 s(-1) [Devriendt, K.; Van Look, H.; Ceursters, B.; Peeters, J. Chem. Phys. Lett. 1996, 261, 450] and a k2b/k1e determination of this work of 1200 +/- 500 at 295 K. Using the temperature dependence of the rate constant k1e(T)/k1e(300 K), which is much more accurately and precisely determined than is A, we predict an increase in k(1e) of a factor 60 +/- 16 between 300 and 1500 K. The ratio of rate constants k2b/k1e is predicted to change from 1200 +/- 500 at 295 K to 40 +/- 25 at 1500 K. These results suggest that the reaction C2H + O2 --> CH(A) + CO2 contributes significantly to CH(A-->X) chemiluminescence in hot flames and especially under fuel-lean conditions where it probably dominates the reaction C2H + O(3P) --> CH(A) + CO.     

Kinetics of the O + ICN reaction

The kinetics of the O + ICN reaction was studied using a relative rate method, with O + C(2)H(2) as the competing reaction. Carbon monoxide products formed in the competing reaction and subsequent secondary chemistry were detected as a function of reagent ICN pressure to obtain total rate constants for the O + ICN reaction. Analysis of the experimental data yields rate constants of k(1) = (3.7 ± 1.0 to 26.2 ± 4.0) × 10(-14) cm(3) molecule(-1) s(-1) over the total pressure range 1.5-9.5 Torr. Product channel NCO + I, the only bimolecular exothermic channel of the reaction, was investigated by detection of N(2)O in the presence of NO and found to be insignificant. An ab initio calculation of the potential energy surface (PES) of the reaction at the CCSD(T)/CEP-31G//DFT-B3LYP/CEP-31G level of theory was also performed. The pathways leading to bimolecular product channels are kinetically unfavorable. Formation and subsequent stabilization of an ICNO adduct species appears to dominate the reaction, in agreement with the experimentally observed pressure dependent rate constants.     

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.     

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.     

Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

Theoretical investigation of the hydrogen abstraction reaction of the OH radical with CH2FCH2F (HFC-152): a dual-level direct dynamics study

The hydrogen abstraction reaction of the OH radical with CH(2)FCH(2)F (HFC-152) is studied theoretically over the 150-3000 K temperature range. In this study, the two most recently developed hybrid density functional theories, namely, BB1K and MPWB1K, are applied, and their efficiency in reaction dynamics calculation is discussed. The BB1K/6-31+G(d,p) method gives the best result for the potential energy surface (PES) calculations, including barrier heights, reaction path information (the first and second derivatives of PES), geometry of transition state structures, and even weak hydrogen bond orientations. The rate constants were obtained by the dual-level direct dynamics with the interpolated single-point energy method (VTST-ISPE) using the BB1K/MG3S//BB1K/6-31+G(d,p) quantum model. The canonical variational transition state theory (CVT) with the small-curvature tunneling correction methods are used to calculate the rate constants in comparison to the experimental data. The total rate constant and its temperature dependency in the form of a fitted three-parameter Arrhenius expression is k(T) = 5.4 x 10(-13)(T/298)3.13 exp{-322/T} cm3 molecule(-1) s(-1). A significant variational effect, which is not common generally for hydrogen-transfer reactions, is reported and analyzed.