Reactions of C2H with 1- and 2-butynes: an ab initio/RRKM study of the reaction mechanism and product branching ratios

Ab initio CCSD(T)/cc-pVTZ(CBS)//B3LYP/6-311G** calculations of the C(6)H(7) potential energy surface are combined with RRKM calculations of reaction rate constants and product branching ratios to investigate the mechanism and product distribution in the C(2)H + 1-butyne/2-butyne reactions. 2-Ethynyl-1,3-butadiene (C(6)H(6)) + H and ethynylallene (C(5)H(4)) + CH(3) are predicted to be the major products of the C(2)H + 1-butyne reaction. The reaction is initiated by barrierless ethynyl additions to the acetylenic C atoms in 1-butyne and the product branching ratios depend on collision energy and the direction of the initial C(2)H attack. The 2-ethynyl-1,3-butadiene + H products are favored by the central C(2)H addition to 1-butyne, whereas ethynylallene + CH(3) are preferred for the terminal C(2)H addition. A relatively minor product favored at higher collision energies is diacetylene + C(2)H(5). Three other acyclic C(6)H(6) isomers, including 1,3-hexadiene-5-yne, 3,4-hexadiene-1-yne, and 1,3-hexadiyne, can be formed as less important products, but the production of the cyclic C(6)H(6) species, fulvene, and dimethylenecyclobut-1-ene (DMCB), is predicted to be negligible. The qualitative disagreement with the recently measured experimental product distribution of C(6)H(6) isomers is attributed to a possible role of the secondary 2-ethynyl-1,3-butadiene + H reaction, which may generate fulvene as a significant product. Also, the photoionization energy curve assigned to DMCB in experiment may originate from vibrationally excited 2-ethynyl-1,3-butadiene molecules. For the C(2)H + 2-butyne reaction, the calculations predict the C(5)H(4) isomer methyldiacetylene + CH(3) to be the dominant product, whereas very minor products include the C(6)H(6) isomers 1,1-ethynylmethylallene and 2-ethynyl-1,3-butadiene.     

An experimental and theoretical study of the product distribution of the reaction CH2 (X 3B1) + NO

Measurements of the product branching ratios of the reaction CH2 (X 3B1) + NO (1) are presented together with calculations of the thermal rate constant and branching ratios using unimolecular rate theory. The reaction was investigated experimentally at room temperature using FTIR spectroscopy. The yields of the main products HCNO and HCN were found to be gamma HCNO = 0.89 +/- 0.06, gamma HCN = 0.11 +/- 0.06. Other minor products could be rationalized by numerical simulations of the reaction system taking into account possible consecutive reactions. The potential energy surface for the reaction was characterized by quantum chemical calculations using ab initio and density functional methods. The proposed reaction pathways connecting reactants to products were explored by multi-channel unimolecular rate theory calculations to determine the CH2 (X) + NO capture rate constant and the rate constants for the different product channels as a function of temperature. The calculated capture rate constant of k = 2.3 x 10(13) cm3 mol-1 s-1 is in good agreement with experimental values at room temperature. Collisional stabilization of the initial H2CNO recombination complex was predicted to be negligible up to pressures of > 1 bar. For ambient pressures and temperatures up to 2000 K, HCNO + H were calculated as the dominating products, with gamma HCNO approximately 0.94 in agreement with the experiments. The channel to HCN + OH was calculated with 0.015 < or = gamma HCN < or = 0.05, only slightly below the experimental value.     

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.     

Dissociation channels of the 1-buten-2-yl radical and its photolytic precursor 2-bromo-1-butene

The work presented here is the first in a series of studies that use a molecular beam scattering technique to investigate the unimolecular reaction dynamics of C(4)H(7) radical isomers. Photodissociation of the halogenated precursor 2-bromo-1-butene at 193 nm under collisionless conditions produced 1-buten-2-yl radicals with a range of internal energies spanning the predicted barriers to the unimolecular reaction channels of the radical. Resolving the velocities of the stable C(4)H(7) radicals, as well as those of the products, allows for the identification of the energetic onset of each dissociation channel. The data show that radicals with at least 30.7 +/- 2 kcal/mol of internal energy underwent C-C fission to form allene + methyl, and radicals with at least 36.7 +/- 4 kcal/mol of internal energy underwent C-H fission to form H + 1-butyne and H + 1,2-butadiene; both of these observed barriers agree well with the G3//B3LYP calculations of Miller. HBr elimination from the parent molecule was observed, producing vibrationally excited 1-butyne and 1,2-butadiene. In the subsequent dissociation of these C(4)H(6) isomers, the major channel was C-C fission to form propargyl + methyl, and there is also evidence of at least one of the possible H + C(4)H(5) channels. A minor C-Br fission channel produces 1-buten-2-yl radicals in an excited electronic state and with low kinetic energy; these radicals exhibit markedly different dissociation dynamics than do the radicals produced in their ground electronic state.     

High-level direct-dynamics variational transition state theory calculations including multidimensional tunneling of the thermal rate constants, branching ratios, and kinetic isotope effects of the hydrogen abstraction reactions from methanol by atomic hydrogen

We report a detailed theoretical study of the hydrogen abstraction reaction from methanol by atomic hydrogen. The study includes the analysis of thermal rate constants, branching ratios, and kinetic isotope effects. Specifically, we have performed high-level computations at the MC3BB level together with direct dynamics calculations by canonical variational transition state theory (CVT) with the microcanonically optimized multidimensional tunneling (μOMT) transmission coefficient (CVT/μOMT) to study both the CH(3)OH+H→CH(2)OH+H(2) (R1) reaction and the CH(3)OH+H→CH(3)O+H(2) (R2) reaction. The CVT/μOMT calculations show that reaction R1 dominates in the whole range 298≤T (K)≤2500 and that anharmonic effects on the torsional mode about the C-O bond are important, mainly at high temperatures. The activation energy for the total reaction sum of R1 and R2 reactions changes substantially with temperature and, therefore, the use of straight-line Arrhenius plots is not valid. We recommend the use of new expressions for the total R1 + R2 reaction and for the R1 and R2 individual reactions.     

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.     

A study of the unimolecular dissociation of the 2-buten-2-yl radical via the 193 nm photodissociation of 2-chloro-2-butene

This work investigates the unimolecular dissociation of the 2-buten-2-yl radical. This radical has three potentially competing reaction pathways: C-C fission to form CH3 + propyne, C-H fission to form H + 1,2-butadiene, and C-H fission to produce H + 2-butyne. The experiments were designed to probe the branching to the three unimolecular dissociation pathways of the radical and to test theoretical predictions of the relevant dissociation barriers. Our crossed laser-molecular beam studies show that 193 nm photolysis of 2-chloro-2-butene produces 2-buten-2-yl in the initial photolytic step. A minor C-Cl bond fission channel forms electronically excited 2-buten-2-yl radicals and the dominant C-Cl bond fission channel produces ground-state 2-buten-2-yl radicals with a range of internal energies that spans the barriers to dissociation of the radical. Detection of the stable 2-buten-2-yl radicals allows a determination of the translational, and therefore internal, energy that marks the onset of dissociation of the radical. The experimental determination of the lowest-energy dissociation barrier gave 31 +/- 2 kcal/mol, in agreement with the 32.8 +/- 2 kcal/mol barrier to C-C fission at the G3//B3LYP level of theory. Our experiments detected products of all three dissociation channels of unstable 2-buten-2-yl as well as a competing HCl elimination channel in the photolysis of 2-chloro-2-butene. The results allow us to benchmark electronic structure calculations on the unimolecular dissociation reactions of the 2-buten-2-yl radical as well as the CH3 + propyne and H + 1,2-butadiene bimolecular reactions. They also allow us to critique prior experimental work on the H + 1,2-butadiene reaction.     

Calculational and experimental investigations of the pressure effects on radical-radical cross combination reactions: C2H5+C2H3

Pressure-dependent product yields have been experimentally determined for the cross-radical reaction C2H5 + C2H3. These results have been extended by calculations. It is shown that the chemically activated combination adduct, 1-C4H8*, is either stabilized by bimolecular collisions or subject to a variety of unimolecular reactions including cyclizations and decompositions. Therefore the "apparent" combination/disproportionation ratio exhibits a complex pressure dependence. The experimental studies were performed at 298 K and at selected pressures between about 4 Torr (0.5 kPa) and 760 Torr (101 kPa). Ethyl and vinyl radicals were simultaneously produced by 193 nm excimer laser photolysis of C2H5COC2H3 or photolysis of C2H3Br and C2H5COC2H5. Gas chromatograph/mass spectrometry/flame ionization detection (GC/MS/FID) were used to identify and quantify the final reaction products. The major combination reactions at pressures between 500 (66.5 kPa) and 760 Torr are (1c) C2H5+C2H3-->1-butene, (2c) C2H5 + C2H5-->n-butane, and (3c) C2H3+C2H3-->1,3-butadiene. The major products of the disproportionation reactions are ethane, ethylene, and acetylene. At moderate and lower pressures, secondary products, including propene, propane, isobutene, 2-butene (cis and trans), 1-pentene, 1,4-pentadiene, and 1,5-hexadiene are also observed. Two isomers of C4H6, cyclobutene and/or 1,2-butadiene, were also among the likely products. The pressure-dependent yield of the cross-combination product, 1-butene, was compared to the yield of n-butane, the combination product of reaction (2c), which was found to be independent of pressure over the range of this study. The [1-C4H8]/[C4H10] ratio was reduced from approximately 1.2 at 760 Torr (101 kPa) to approximately 0.5 at 100 Torr (13.3 kPa) and approximately 0.1 at pressures lower than about 5 Torr (approximately 0.7 kPa). Electronic structure and RRKM calculations were used to simulate both unimolecular and bimolecular processes. The relative importance of C-C and C-H bond ruptures, cyclization, decyclization, and complex decompositions are discussed in terms of energetics and structural properties. The pressure dependence of the product yields were computed and dominant reaction paths in this chemically activated system were determined. Both modeling and experiment suggest that the observed pressure dependence of [1-C4H8]/[C4H10] is due to decomposition of the chemically activated combination adduct 1-C4H8* in which the weaker allylic C-C bond is broken: H2C=CHCH2CH3-->C3H5+CH3. This reaction occurs even at moderate pressures of approximately 200 Torr (26 kPa) and becomes more significant at lower pressures. The additional products detected at lower pressures are formed from secondary radical-radical reactions involving allyl, methyl, ethyl, and vinyl radicals. The modeling studies have extended the predictions of product distributions to different temperatures (200-700 K) and a wider range of pressures (10(-3)-10(5) Torr). These calculations indicate that the high-pressure [1-C4H8]/[C4H10] yield ratio is 1.3+/-0.1.     

Formation and decomposition of chemically activated cyclopentoxy radicals from the c-C5H9 + O reaction

The formation and the decomposition of chemically activated cyclopentoxy radicals from the c-C5H9 + O reaction have been studied in the gas phase at room temperature. Two different experimental arrangements have been used. Arrangement A consisted of a laser-flash photolysis set up combined with quantitative Fourier transform infrared spectroscopy and allowed the determination of the stable products at 4 mbar. The c-C5H9 radicals were produced via the reaction c-C5H10 + Cl with chlorine atoms from the photolysis of CFCl3; the O atoms were generated by photolysis of SO2. Arrangement B, a conventional discharge flow-reactor with molecular beam sampling, was used to determine the rate coefficient. Here, the hydrocarbon radicals (c-C5H9, C2H5, CH2OCH3) were produced via the reaction of atomic fluorine with c-C5H10, C2H6, and CH3OCH3, respectively, and detected by mass spectrometry after laser photoionization. For the c-C5H9 + O reaction, the relative contributions of intermediate formation (c-C5H9O) and direct abstraction (c-C5H8 + OH) were found to be 68 +/- 5 and 32 +/- 4%, respectively. The decomposition products of the chemically activated intermediate could be identified, and the following relative branching fractions were obtained: c-C5H8O + H (31 +/- 2%), CH2CH(CH2)2CHO + H (40 +/- 5%), 2 C2H4 + H + CO (17 +/- 5%), and C3H4O + C2H4 + H (12 +/- 5%). Additionally, the product formation of the c-C5H8 + O reaction was studied, and the following relative yields were obtained (mol %): C2H4, 24%; C3H4O, 18%; c-C5H8O, 30%; c-C5H8O, 23%; 4-pentenal, 5%. The rate coefficient of the c-C5H9 + O reaction was determined relative to the reactions C2H5 + O and CH3OCH2 + O leading to k = (1.73 +/- 0.05) x 10(14) cm3 mol(-1) s(-1). The experimental branching fractions are analyzed in terms of statistical rate theory with molecular and transition-state data from quantum chemical calculations, and high-pressure limiting Arrhenius parameters for the unimolecular decomposition reactions of C5H9O species are derived.     

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.     

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.     

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.     

Ab initio energies and product branching ratios for the O+C3H6 reaction

Intermediate and transition-state energies have been calculated for the O+C3H6 (propene) reaction using the compound ab initio CBS-QB3 and G3 methods in combination with density functional theory. The lowest-lying triplet and singlet potential energy surfaces of the O-C3H6 system were investigated. RRKM statistical theory was used to predict product branching fractions over the 300-3000 K temperature and 0.001-760 Torr pressure ranges. The oxygen atom adds to the C3H6 terminal olefinic carbon in the primary step to generate a nascent triplet biradical, CH3CHCH2O. On the triplet surface, unimolecular dissociation of CH3CHCH2O to yield H+CH3CHCHO is favored over the entire temperature range, although the competing H2CO+CH3CH product channel becomes significant at high temperature. Rearrangement of triplet CH3CHCH2O to CH3CH2CHO (propanal) via a 1,2 H-atom shift has a barrier of 122.3 kJ mol(-1), largely blocking this reaction channel and any subsequent dissociation products. Intersystem crossing of triplet CH3CHCH2O to the singlet surface, however, leads to facile rearrangement to singlet CH3CH2CHO, which dissociates via numerous product channels. Pressure was found to have little influence over the branching ratios under most conditions, suggesting that the vibrational self-relaxation rates for p相似文献   

Associative charge transfer reactions. Temperature effects and mechanism of the gas-phase polymerization of propene initiated by a benzene radical cation

In associative charge transfer (ACT) reactions, a core ion activates ligand molecules by partial charge transfer. The activated ligand polymerizes, and the product oligomer takes up the full charge from the core ion. In the present system, benzene(+*) (Bz(+*)) reacts with two propene (Pr) molecules to form a covalently bonded ion, C(6)H(6)(+*) + 2 C(3)H(6) --> C(6)H(12)(+*) + C(6)H(6). The ACT reaction is activated by a partial charge transfer from Bz(+*) to Pr in the complex, and driven to completion by the formation of a covalent bond in the polymerized product. An alternative channel forms a stable association product (Bz.Pr)(+*), with an ACT/association product ratio of 60:40% that is independent of pressure and temperature. In contrast to the Bz(+*)/propene system, ACT polymerization is not observed in the Bz(+*)/ethylene (Et) system since charge transfer in the Bz(+*)(Et) complex is inefficient to activate the reaction. The roles of charge transfer in these complexes are verified by ab initio calculations. The overall reaction of Bz(+*) with Pr follows second-order kinetics with a rate constant of k (304 K) = 2.1 x 10(-12) cm(3) s(-1) and a negative temperature coefficient of k = aT(-5.9) (or an activation energy of -3 kcal/mol). The kinetic behavior is similar to sterically hindered reactions and suggests a [Bz(+*) (Pr)]* activated complex that proceeds to products through a low-entropy transition state. The temperature dependence shows that ACT reactions can reach a unit collision efficiency below 100 K, suggesting that ACT can initiate polymerization in cold astrochemical environments.     

Potential energy surface and product branching ratios for the reaction of dicarbon, C(2)(X(1)Sigma(g)(+)), with methylacetylene, CH(3)CCH(X(1)A(1)): an ab initio/RRKM study

Ab initio calculations of the potential energy surface for the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH(X(1)A(1)) reaction have been carried at the G2M level of theory. The calculations show that the dicarbon molecule in the ground singlet electronic state can add to methylacetylene without a barrier producing a three-member or a four-member ring intermediate, which can rapidly rearrange to the most stable H(3)CCCCCH isomer on the C(5)H(4) singlet surface. This isomer can then lose a hydrogen atom (H) or molecular hydrogen (H(2)) from the CH(3) group with the formation of H(2)CCCCCH and HCCCCCH, respectively. Alternatively, H atom migrations and three-member-ring closure/opening rearrangements followed by H and H(2) losses can lead to other isomers of the C(5)H(3) and C(5)H(2) species. According to the calculated energetics, the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH reaction is likely to be a major source of the C(5)H(3) radicals (in particular, the most stable H(2)CCCCCH and HCCCHCCH isomers, which are relevant to the formation of benzene through the reactions with CH(3)). Among heavy-fragment product channels, only C(3)H(3) + C(2)H and c-C(3)H(2) + C(2)H(2) might compete with C(5)H(3) + H and C(5)H(2) + H(2). RRKM calculations of reaction rate constants and product branching ratios depending on the reactive collision energy showed that the major reaction products are expected to be H(2)CCCCCH + H (64-66%) and HCCCHCCH + H (34-30%), with minor contributions from HCCCCCH + H(2) (1-2%), HCCCHCC + H(2) (up to 1%), C(3)H(3) + C(2)H (up to 1%), and c-C(3)H(2) + C(2)H(2) (up to 0.1%) if the energy randomization is complete. The calculations also indicate that the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH(X(1)A(1)) reaction can proceed by direct H-abstraction of a methyl hydrogen to form C(3)H(3) + C(2)H almost without a barrier.     

H2FCS单分子分解反应的直接动力学研究

采用量子化学QCISD(T)/6-311++G(d,p)//B3LYP/6-311+G(d,p)方法研究了H2FCS单分子分解反应的微观动力学性质, 构建了反应势能剖面. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT), 分别计算了在200～3000 K温度范围内的速率常数kTST、kCVT和kCVT/SCT. 计算结果表明, H2FCS可经过不同的反应通道生成10种小分子产物, 脱H反应和HF消去反应为标题反应的主反应通道, 其中HF消去反应产物HCS可由两条反应通道生成. 在200~3000 K温度区间内得到三条反应通道的表观反应速率常数三参数表达式分别为 , 和 . 速率常数计算结果显示, 量子力学隧道效应在低温区间对反应速率常数的影响显著, 而变分效应在计算温度范围内可以忽略.     

Statistical theory for the kinetics and dynamics of roaming reactions

We present a statistical theory for the effect of roaming pathways on product branching fractions in both unimolecular and bimolecular reactions. The analysis employs a separation into three distinct steps: (i) the formation of weakly interacting fragments in the long-range/van der Waals region of the potential via either partial decomposition (for unimolecular reactants) or partial association (for bimolecular reactants), (ii) the roaming step, which involves the reorientation of the fragments from one region of the long-range potential to another, and (iii) the abstraction, addition, and/or decomposition from the long-range region to yield final products. The branching between the roaming induced channel(s) and other channels is obtained from a steady-state kinetic analysis for the two (or more) intermediates in the long-range region of the potential. This statistical theory for the roaming-induced product branching is illustrated through explicit comparisons with reduced dimension trajectory simulations for the decompositions of H(2)CO, CH(3)CHO, CH(3)OOH, and CH(3)CCH. These calculations employ high-accuracy analytic potentials obtained from fits to wide-ranging CASPT2 ab initio electronic structure calculations. The transition-state fluxes for the statistical theory calculations are obtained from generalizations of the variable reaction coordinate transition state theory approach. In each instance, at low energy the statistical analysis accurately reproduces the branching obtained from the trajectory simulations. At higher energies, e.g., above 1 kcal/mol, increasingly large discrepancies arise, apparently due to a dynamical biasing toward continued decomposition of the incipient molecular fragments (for unimolecular reactions). Overall, the statistical theory based kinetic analysis is found to provide a useful framework for interpreting the factors that determine the significance of roaming pathways in varying chemical environments.     

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.     

Ab initio chemical kinetics for the ClOO + NO reaction: effects of temperature and pressure on product branching formation

The kinetics and mechanism for the reaction of ClOO with NO have been investigated by ab initio molecular orbital theory calculations based on the CCSD(T)/6-311+G(3df)//PW91PW91∕6-311+G(3df) method, employed to evaluate the energetics for the construction of potential energy surfaces and prediction of reaction rate constants. The results show that the reaction can produce two key low energy products ClNO + (3)O(2) via the direct triplet abstraction path and ClO + NO(2) via the association and decomposition mechanism through long-lived singlet pc-ClOONO and ClONO(2) intermediates. The yield of ClNO + O(2) ((1)△) from any of the singlet intermediates was found to be negligible because of their high barriers and tight transition states. As both key reactions initially occur barrierlessly, their rate constants were evaluated with a canonical variational approach in our transition state theory and Rice-Ramspergen-Kassel-Marcus/master equation calculations. The rate constants for ClNO + (3)O(2) and ClO + NO(2) production from ClOO + NO can be given by 2.66 × 10(-16) T(1.91) exp(341/T) (200-700 K) and 1.48 × 10(-24) T(3.99) exp(1711/T) (200-600 K), respectively, independent of pressure below atmospheric pressure. The predicted total rate constant and the yields of ClNO and NO(2) in the temperature range of 200-700 K at 10-760 Torr pressure are in close agreement with available experimental results.     

Thermochemistry and reaction paths in the oxidation reaction of benzoyl radical: C6H5C•(═O)

Alkyl substituted aromatics are present in fuels and in the environment because they are major intermediates in the oxidation or combustion of gasoline, jet, and other engine fuels. The major reaction pathways for oxidation of this class of molecules is through loss of a benzyl hydrogen atom on the alkyl group via abstraction reactions. One of the major intermediates in the combustion and atmospheric oxidation of the benzyl radicals is benzaldehyde, which rapidly loses the weakly bound aldehydic hydrogen to form a resonance stabilized benzoyl radical (C6H5C(?)═O). A detailed study of the thermochemistry of intermediates and the oxidation reaction paths of the benzoyl radical with dioxygen is presented in this study. Structures and enthalpies of formation for important stable species, intermediate radicals, and transition state structures resulting from the benzoyl radical +O2 association reaction are reported along with reaction paths and barriers. Enthalpies, ΔfH298(0), are calculated using ab initio (G3MP2B3) and density functional (DFT at B3LYP/6-311G(d,p)) calculations, group additivity (GA), and literature data. Bond energies on the benzoyl and benzoyl-peroxy systems are also reported and compared to hydrocarbon systems. The reaction of benzoyl with O2 has a number of low energy reaction channels that are not currently considered in either atmospheric chemistry or combustion models. The reaction paths include exothermic, chain branching reactions to a number of unsaturated oxygenated hydrocarbon intermediates along with formation of CO2. The initial reaction of the C6H5C(?)═O radical with O2 forms a chemically activated benzoyl peroxy radical with 37 kcal mol(-1) internal energy; this is significantly more energy than the 21 kcal mol(-1) involved in the benzyl or allyl + O2 systems. This deeper well results in a number of chemical activation reaction paths, leading to highly exothermic reactions to phenoxy radical + CO2 products.