Radical-molecule reactions HCO/HOC + C2H2: mechanistic study

A detailed computational study is performed on the unknown radical-molecule reactions between HCO/HOC and acetylene (C2H2) at the CCSD(T)/6-311G(2d,p)//B3LYP/6-311G(d,p)+ZPVE, Gaussian-3//B3LYP/6-31G(d), and Gaussian-3//MP2(full)/6-31G(d) levels. For the HCO + C2H2 reaction, the most favorable pathway is direct C-addition forming the intermediate HC=CHCH=O followed by a 1,3-H-shift leading to H2C=CHC=O, which finally dissociates to the product C2H3 + CO. The overall reaction barrier is 13.8, 10.5, and 11.3 kcal/mol, respectively, at the three levels. The quasi-direct H-donation process to produce C2H3 + CO with barriers of 14.0, 14.1, and 14.1 kcal/mol is less competitive. Thus only at higher temperatures could the HCO + C2H2 reaction play a role. In contrast, the HOC + C2H2 reaction can barrierlessly generate C2H3 + CO via the quasi-direct H-donation mechanism proceeding via a prereactive complex with OH...C2 hydrogen bonding. This is suggestive of the potential importance of the HOC + C2H2 reaction in both combustion and interstellar processes. However, the direct C-addition channel is much less competitive. For both reactions, the possible formation of the intriguing interstellar molecules propadiene and propynal is also discussed. The present theoretical study represents the first attempt to probe the reaction mechanism between HOC and pi-systems. Future laboratory investigations on both reactions (particularly HOC + C2H2) are recommended.     

Reaction of the i-C4H5 (CH2CCHCH2) radical with O2

The resonantly stabilized radical i-C(4)H(5) (CH(2)CCHCH(2)) is an important intermediate in the combustion of unsaturated hydrocarbons and is thought to be involved in the formation of polycyclic aromatic hydrocarbons through its reaction with acetylene (C(2)H(2)) to form benzene + H. This study uses quantum chemistry and statistical reaction rate theory to investigate the mechanism and kinetics of the i-C(4)H(5) + O(2) reaction as a function of temperature and pressure, and unlike most resonantly stabilized radicals we show that i-C(4)H(5) is consumed relatively rapidly by its reaction with molecular oxygen. O(2) addition occurs at the vinylic and allenic radical sites in i-C(4)H(5), with respective barriers of 0.9 and 4.9 kcal mol(-1). Addition to the allenic radical form produces an allenemethylperoxy radical adduct with only around 20 kcal mol(-1) excess vibrational energy. This adduct can isomerize to the ca. 14 kcal mol(-1) more stable 1,3-divinyl-2-peroxy radical via concerted and stepwise processes, both steps with barriers around 10 kcal mol(-1) below the entrance channel energy. Addition of O(2) to the vinylic radical site in i-C(4)H(5) directly forms the 1,3-divinyl-2-peroxy radical with a small barrier and around 36.8 kcal mol(-1) of excess energy. The 1,3-divinyl-2-peroxy radical isomerizes via ipso addition of the O(2) moiety followed by O atom insertion into the adjacent C-C bond. This process forms an unstable intermediate that ultimately dissociates to give the vinyl radical, formaldehyde, and CO. At higher temperatures formation of vinylacetylene + HO(2), the vinoxyl radical + ketene, and the 1,3-divinyl-2-oxyl radical + O paths have some importance. Because of the adiabatic transition states for O(2) addition, and significant reverse dissociation channels in the peroxy radical adducts, the i-C(4)H(5) + O(2) reaction proceeds to new products with rate constant of around 10(11) cm(3) mol(-1) s(-1) at typical combustion temperatures (1000-2000 K). For fuel-rich flames we show that the reaction of i-C(4)H(5) with O(2) is likely to be faster than that with C(2)H(2), bringing into question the importance of the i-C(4)H(5) + C(2)H(2) reaction in initiating ring formation in sooting flames.     

First principle and ReaxFF molecular dynamics investigations of formaldehyde dissociation on Fe(100) surface

Detailed formaldehyde adsorption and dissociation reactions on Fe(100) surface were studied using first principle calculations and molecular dynamics (MD) simulations, and results were compared with available experimental data. The study includes formaldehyde, formyl radical (HCO), and CO adsorption and dissociation energy calculations on the surface, adsorbate vibrational frequency calculations, density of states analysis of clean and adsorbed surfaces, complete potential energy diagram construction from formaldehyde to atomic carbon (C), hydrogen (H), and oxygen (O), simulation of formaldehyde adsorption and dissociation reaction on the surface using reactive force field, ReaxFF MD, and reaction rate calculations of adsorbates using transition state theory (TST). Formaldehyde and HCO were adsorbed most strongly at the hollow (fourfold) site. Adsorption energies ranged from ?22.9 to ?33.9 kcal/mol for formaldehyde, and from ?44.3 to ?66.3 kcal/mol for HCO, depending on adsorption sites and molecular direction. The dissociation energies were investigated for the dissociation paths: formaldehyde → HCO + H, HCO → H + CO, and CO → C + O, and the calculated energies were 11.0, 4.1, and 26.3 kcal/mol, respectively. ReaxFF MD simulation results were compared with experimental surface analysis using high resolution electron energy loss spectrometry (HREELS) and TST based reaction rates. ReaxFF simulation showed less reactivity than HREELS observation at 310 and 523 K. ReaxFF simulation showed more reactivity than the TST based rate for formaldehyde dissociation and less reactivity than TST based rate for HCO dissociation at 523 K. TST‐based rates are consistent with HREELS observation. © 2013 Wiley Periodicals, Inc.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Mechanism of HCS + O2 reaction: hydrogen- or oxygen-transfer?

In spite of the potential importance of the HCS radical in both combustion and interstellar processes, its chemical reactivity has not been tackled previously. In the present paper, the oxidation reaction of the HCS radical is theoretically investigated for the first time at the CCSD(T)/6-311++G(3df,2p)//BH&HLYP/6-311++G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels. It is shown that the most feasible pathway is the O2 addition to the HCS radical forming the intermediate SC(H)OO which can undergo a subsequent O-extrusion leading to SC(H)O + 3O. This features an indirect O-transfer mechanism with the overall barrier of 4.4 and 3.5 kcal mol(-1), respectively, at the two levels. However, formation of the H-transfer product CS + HO2 is kinetically much less feasible, i.e., the direct mechanism has barriers of 14.3 and 8.7 kcal mol(-1), whereas the indirect mechanism has barriers of 12.6 and 10.7 kcal mol(-1), respectively. This result is in sharp contrast to the analogous HCO + O2 reaction, where the direct (with a barrier of 2.98 kcal mol(-1)) and indirect (2.26 kcal mol(-1)) H-transfer processes are highly competitive over the indirect O-transfer process (the least endothermicity is 19.9 kcal mol(-1)). The possible explanations and implications of the present results are provided.     

Direct versus Hydrogen‐Assisted CO Dissociation on the Fe (100) Surface: a DFT Study

CO dissociation: Three most probable pathways to CO dissociation on the Fe?(100) surface exist: a) direct, CO→C+O (-) and H-assisted b) H+CO?HCO→CH+O (-) or c) CO+H?COH→C+OH (-). Under high hydrogen pressure conditions and highly occupied surfaces the formation of HCO and subsequent dissociation to CH+O may at best compete with direct dissociation.     

Reaction mechanism of naphthyl radicals with molecular oxygen. 1. Theoretical study of the potential energy surface

Potential energy surfaces (PESs) of the reactions of 1- and 2-naphthyl radicals with molecular oxygen have been investigated at the G3(MP2,CC)//B3LYP/6-311G** level of theory. Both reactions are shown to be initiated by barrierless addition of O(2) to the respective radical sites of C(10)H(7). The end-on O(2) addition leading to 1- and 2-naphthylperoxy radicals exothermic by 45-46 kcal/mol is found to be more preferable thermodynamically than the side-on addition. At the subsequent reaction step, the chemically activated 1- and 2-C(10)H(7)OO adducts can eliminate an oxygen atom leading to the formation of 1- and 2-naphthoxy radical products, respectively, which in turn can undergo unimolecular decomposition producing indenyl radical + CO via the barriers of 57.8 and 48.3 kcal/mol and with total reaction endothermicities of 14.5 and 10.2 kcal/mol, respectively. Alternatively, the initial reaction adducts can feature an oxygen atom insertion into the attacked C(6) ring leading to bicyclic intermediates a10 and a10' (from 1-naphthyl + O(2)) or b10 and b10' (from 2-naphthyl + O(2)) composed from two fused six-member C(6) and seven-member C(6)O rings. Next, a10 and a10' are predicted to decompose to C(9)H(7) (indenyl) + CO(2), 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H, and 1-C(9)H(7)O (1-benzopyranyl) + CO, whereas b10 and b10' would dissociate to C(9)H(7) (indenyl) + CO(2), 2-C(9)H(7)O (2-benzopyranyl) + CO, and 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H. On the basis of this, the 1-naphthyl + O(2) reaction is concluded to form the following products (with the overall reaction energies given in parentheses): 1-naphthoxy + O (-15.5 kcal/mol), indenyl + CO(2) (-123.9 kcal/mol), 1-benzopyranyl + CO (-97.2 kcal/mol), and 1,2-naphthoquinone + H (-63.5 kcal/mol). The 2-naphthyl + O(2) reaction is predicted to produce 2-naphthoxy + O (-10.9 kcal/mol), indenyl + CO(2) (-123.7 kcal/mol), 2-benzopyranyl + CO (-90.7 kcal/mol), and 1,2-naphthoquinone + H (-63.2 kcal/mol). Simplified kinetic calculations using transition-state theory computed rate constants at the high-pressure limit indicate that the C(10)H(7)O + O product channels are favored at high temperatures, while the irreversible oxygen atom insertion first leading to the a10 and a10' or b10 and b10' intermediates and then to their various decomposition products is preferable at lower temperatures. Among the decomposition products, indenyl + CO(2) are always most favorable at lower temperatures, but the others, 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H (from a10 and b10'), 1-C(9)H(7)O (1-benzopyranyl) + CO (from a10'), and 2-C(10)H(7)O (2-benzopyranyl) + O (from b10 and minor from b10'), may notably contribute or even become major products at higher temperatures.     

A direct dynamics trajectory study of F- + CH(3)OOH reactive collisions reveals a major non-IRC reaction path

A direct dynamics simulation at the B3LYP/6-311+G(d,p) level of theory was used to study the F- + CH3OOH reaction dynamics. The simulations are in excellent agreement with a previous experimental study (J. Am. Chem. Soc. 2002, 124, 3196). Two product channels, HF + CH2O + OH- and HF + CH3OO-, are observed. The former dominates and occurs via an ECO2 mechanism in which F- attacks the CH3- group, abstracting a proton. Concertedly, a carbon-oxygen double bond is formed and OH- is eliminated. Somewhat surprisingly this is not the reaction path, predicted by the intrinsic reaction coordinate (IRC), which leads to a deep potential energy minimum for the CH2(OH)2...F- complex followed by dissociation to HF + CH2(OH)O-. None of the direct dynamics trajectories followed this path, which has an energy release of -63 kcal/mol and is considerably more exothermic than the ECO2 path whose energy release is -27 kcal/mol. Other product channels not observed, and which have a lower energy than that for the ECO2 path, are F- + CO + H2 + H2O (-43 kcal/mol), F- + CH2O + H2O (-51 kcal/mol), and F- + CH2(OH)2 (-60 kcal/mol). Formation of the CH3OOH...F- complex, with randomization of its internal energy, is important, and this complex dissociates via the ECO2 mechanism. Trajectories which form HF + CH3OO- are nonstatistical events and, for the 4 ps direct dynamics simulation, are not mediated by the CH3OOH...F- complex. Dissociation of this complex to form HF + CH3OO- may occur on longer time scales.     

Vibrational distributions of the CO(v) products of the C2H2 + O(3P) and HCCO + O(3P) reactions studied by FTIR emission

The C2H2 + O(3P) and HCCO + O(3P) reactions are investigated using Fourier transform infrared (FTIR) emission spectroscopy. The O(3P) radicals are produced by 193 nm photolysis of an SO2 precursor or microwave discharge in O2. The HCCO radical is either formed in the first step of the C2H2 + O(3P) reaction or by 193 nm photodissociation of ethyl ethynyl ether. Vibrationally excited CO and CO2 products are observed. The microwave discharge experiment [C2H2 + O(3P)] shows a bimodal distribution of the CO(v) product, which is due to the sequential C2H2 + O(3P) and HCCO + O(3P) reactions. The vibrational distribution of CO(v) from the HCCO + O(3P) reaction also shows its own bimodal shape. The vibrational distribution of CO(v) from C2H2 + O(3P) can be characterized by a Boltzmann plot with a vibrational temperature of approximately 2400 +/- 100 K, in agreement with previous results. The CO distribution from the HCCO + O(3P) reaction, when studied under conditions to minimize other processes, shows very little contamination from other reactions, and the distribution can be characterized by a linear combination of Boltzmann plots with two vibrational temperatures: 2320 +/- 40 and 10 300 +/- 600 K. From the experimental results and previous theoretical work, the bimodal CO(v) distribution for the HCCO + O(3P) reaction suggests a sequential dissociation process of the HC(O)CO++ --> CO + HCO; HCO --> H + CO.     

Theoretical study of isomerization and decomposition of propenal

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

C2H+H2CO: a new route for formaldehyde removal

The title unknown reaction is theoretically studied at various levels to probe the interaction mechanism between the ethynyl radical (HC triple bond C) and formaldehyde (H(2)C double bond O). The most feasible pathway is a barrier-free direct H-abstraction process leading to acetylene and formyl radical (C(2)H(2)+HCO) via a weakly bound complex, and then the product can take secondary dissociation to the final product C(2)H(2)+CO+H. The C-addition channel leading to propynal plus H-atom (HCCCHO+H) has the barrier of only 3.6, 2.9, and 2.1 kcal/mol at the CCSD(T)/6-311+G(3df,2p)MP2//6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE, and G3//MP2 levels, respectively [CCSD(T)--coupled cluster with single, double, and triple excitations; ZPVE--zero-point vibrational energy; QCISD--quadratic configuration interaction with single and double excitations; G3//MP2-Gaussian-3 based on Moller-Plesset geometry]. The O addition also leading to propynal plus H atom needs to overcome a higher barrier of 5.3, 8.7, and 3.0 kcalmol at the three corresponding levels. The title no-barrier reaction presents a new efficient route to remove the pollutant H(2)CO, and should be included in the combustion models of hydrocarbons. It may also represent the fastest radical-H(2)CO reaction among the available theoretical data. Moreover, it could play an important role in the interstellar chemistry where the zero- or minute-barrier reactions are generally favored. Discussions are also made on the possible formation of the intriguing propynal in space via the title reaction on ice surface.     

Study of the C(3P) + OH(X2Pi) --> CO(X1Sigma(g)+) + H(2S) reaction: a fully global ab initio potential energy surface of the X2A' state

The C((3)P) + OH(X (2)Pi) --> CO(X (1)Sigma(g)(+)) + H((2)S) reaction has been investigated by ab initio electronic structure calculations of the X(2)A' state based on the multireference (MR) internally contracted single and double configuration interaction (SDCI) method plus Davidson correction (+Q) using Dunning aug-cc-pVQZ basis sets. In particular, the multireference space is taken to be a complete active space (CAS). Improvement over previously proposed potential energy surfaces for HCO/COH is obtained in the sense that present surface describes also the potential part where the CO interatomic distance is large. A large number of geometries (around 2000) have been calculated and analytically fitted using the reproducing kernel Hilbert space (RKHS) method of Ho and Rabitz both for the two-body and three-body terms following the many-body decomposition of the total electronic energies. Results show that the global reaction is highly exothermic ( approximately 6.4 eV) and barrierless (relative to the reactant channel), while five potential barriers are located on this surface. The three minima and five saddle points observed are characterized and found to be in good agreement with previous work. The three minima correspond to the formation of HCO and COH complexes and to the CO + H products, with the COH complex being a metastable minimum relative to the product channel. The five saddle points correspond to potential barriers for both the dissociation/formation of HCO and COH into/from CO + H, to barriers for the isomerization of HCO into COH and to barriers for the inversion of HCO and COH through their respective linear configuration.     

Theoretical studies on cycloaddition reactions between the 2-aza-1,3-butadiene cation and olefins

Density functional (B3LYP) calculations, using the 6-31G basis set, have been employed to study the title reactions. For the model reaction (H(2)C=C-NH(+)=CH(2) + H(2)C=CH(2)), a complex has been formed with 6.2 kcal/mol of stabilization energy and the transition state is 4.0 kcal/mol above this complex, but 2.1 kcal/mol below the reactants. However, the substituent effects are quite remarkable. When ethene is substituted by electron-withdrawing group CN, the reaction could also yield six-membered-ring products, but the energy barriers are all more than 7 kcal/mol, which shows that CN group unfavors the reaction. The other substituents, such as CH(3)O and CH(3) groups, have also been considered in the present work, and the results show that they are favorable for the formation of six-membered-ring adducts. The calculated results have been rationalized with frontier orbital interaction and topological analysis.     

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.     

HCO+NO2反应势能面的理论研究

在B3LYP/6-311G(d,p)和CCSD(T)/6-311G(d,p)水平上给出了HCO+NO₂反应详细的势能面信息.计算结果表明,该反应采用两种无垒进攻方式,分别得到两种加合物H(O)CNO₂和H(O)CONO.找到7种能量低于反应物且合理的产物及相应的反应路径.通过对热力学和动力学的分析,产物HONO+CO(P2,P3),HNO+CO₂(P1)和H+CO₂+NO(P6)的形成更为有利.计算结果同实验相符,且有助于深入了解HCO自由基的化学行为.     

Computational study on kinetics and mechanisms of unimolecular decomposition of succinic acid and its anhydride

The mechanisms and kinetics of unimolecular decomposition of succinic acid and its anhydride have been studied at the G2M(CC2) and microcanonical RRKM levels of theory. It was shown that the ZsgsZ conformer of succinic acid, with the Z-acid form and the gauche conformation around the central C-C bond, is its most stable conformer, whereas the lowest energy conformer with the E-acid form, ECGsZ, is only 3.1 kcal/mol higher in energy than the ZsgsZ. Three primary decomposition channels of succinic acid producing H2O + succinic anhydride with a barrier of 51.0 kcal/mol, H2O + OCC2H3COOH with a barrier of 75.7 kcal/mol and CO2 + C2H5COOH with a barrier of 71.9 kcal/mol were predicted. The dehydration process starting from the ECGCZ-conformer is found to be dominant, whereas the decarboxylation reaction starting from the ZsgsZ-conformer is only slightly less favorable. It was shown that the decomposition of succinic anhydride occurs via a concerted fragmentation mechanism (with a 69.6 kcal/mol barrier), leading to formation of CO + CO2 + C2H4 products. On the basis of the calculated potential energy surfaces of these reactions, the rate constants for unimolecular decomposition of succinic acid and its anhydride were predicted. In addition, the predicted rate constants for the unimolecular decomposition of C2H5COOH by decarboxylation (giving C2H6 + CO2) and dehydration (giving H3CCHCO + H2O) are in good agreement with available experimental data.     

Reaction mechanism between carbonyl oxide and hydroxyl radical: a theoretical study

The reaction mechanism of carbonyl oxide with hydroxyl radical was investigated by using CASSCF, B3LYP, QCISD, CASPT2, and CCSD(T) theoretical approaches with the 6-311+G(d,p), 6-311+G(2df, 2p), and aug-cc-pVTZ basis sets. This reaction involves the formation of H2CO + HO2 radical in a process that is computed to be exothermic by 57 kcal/mol. However, the reaction mechanism is very complex and begins with the formation of a pre-reactive hydrogen-bonded complex and follows by the addition of HO radical to the carbon atom of H2COO, forming the intermediate peroxy-radical H2C(OO)OH before producing formaldehyde and hydroperoxy radical. Our calculations predict that both the pre-reactive hydrogen-bonded complex and the transition state of the addition process lie energetically below the enthalpy of the separate reactants (DeltaH(298K) = -6.1 and -2.5 kcal/mol, respectively) and the formation of the H2C(OO)OH adduct is exothermic by about 74 kcal/mol. Beyond this addition process, further reaction mechanisms have also been investigated, which involve the abstraction of a hydrogen of carbonyl oxide by HO radical, but the computed activation barriers suggest that they will not contribute to the gas-phase reaction of H2COO + HO.     

Probing the barrier for CH2CHCO --> CH2CH + CO by the velocity map imaging method

This work determines the dissociation barrier height for CH2CHCO --> CH2CH + CO using two-dimensional product velocity map imaging. The CH2CHCO radical is prepared under collision-free conditions from C-Cl bond fission in the photodissociation of acryloyl chloride at 235 nm. The nascent CH2CHCO radicals that do not dissociate to CH2CH + CO, about 73% of all the radicals produced, are detected using 157-nm photoionization. The Cl(2P(3/2)) and Cl(2P(1/2)) atomic fragments, momentum matched to both the stable and unstable radicals, are detected state selectively by resonance-enhanced multiphoton ionization at 235 nm. By comparing the total translational energy release distribution P(E(T)) derived from the measured recoil velocities of the Cl atoms with that derived from the momentum-matched radical cophotofragments which do not dissociate, the energy threshold at which the CH2CHCO radicals begin to dissociate is determined. Based on this energy threshold and conservation of energy, and using calculated C-Cl bond energies for the precursor to produce CH2CHC*O or C*H2CHCO, respectively, we have determined the forward dissociation barriers for the radical to dissociate to vinyl + CO. The experimentally determined barrier for CH2CHC*O --> CH2CH + CO is 21+/-2 kcal mol(-1), and the computed energy difference between the CH2CHC*O and the C*H2CHCO forms of the radical gives the corresponding barrier for C*H2CHCO --> CH2CH + CO to be 23+/-2 kcal mol(-1). This experimental determination is compared with predictions from electronic structure methods, including coupled-cluster, density-functional, and composite Gaussian-3-based methods. The comparison shows that density-functional theory predicts too low an energy for the C*H2CHCO radical, and thus too high a barrier energy, whereas both the Gaussian-3 and the coupled-cluster methods yield predictions in good agreement with experiment. The experiment also shows that acryloyl chloride can be used as a photolytic precursor at 235 nm of thermodynamically stable CH2CHC*O radicals, most with an internal energy distribution ranging from approximately 3 to approximately 21 kcal mol(-1). We discuss the results with respect to the prior work on the O(3P) + propargyl reaction and the analogous O(3P) + allyl system.     

The reaction of tricarbon with acetylene: an ab initio/RRKM study of the potential energy surface and product branching ratios

Ab initio calculations of the potential energy surface for the C3(1Sigmag+)+C2H2(1Sigmag+) reaction have been performed at the RCCSD(T)/cc-pVQZ//B3LYP/6-311G(d,p) + ZPE[B3LYP/6-311G(d,p)] level with extrapolation to the complete basis set limit for key intermediates and products. These calculations have been followed by statistical calculations of reaction rate constants and product branching ratios. The results show the reaction to begin with the formation of the 3-(didehydrovinylidene)cyclopropene intermediate i1 or five-member ring isomer i7 with the entrance barriers of 7.6 and 13.8 kcal/mol, respectively. i1 rearranges to the other C5H2 isomers, including ethynylpropadienylidene i2, singlet pentadiynylidene i3, pentatetraenylidene i4, ethynylcyclopropenylidene i5, and four- and five-member ring structures i6, i7, and i8 by ring-closure and ring-opening processes and hydrogen migrations. i2, i3, and i4 lose a hydrogen atom to produce the most stable linear isomer of C5H with the overall reaction endothermicity of approximately 24 kcal/mol. H elimination from i5 leads to the formation of the cyclic C5H isomer, HC2C3, +H, 27 kcal/ mol above C3+C2H2. 1,1-H2 loss from i4 results in the linear pentacarbon C5+H2 products endothermic by 4 kcal/mol. The H elimination pathways occur without exit barriers, whereas the H2 loss from i4 proceeds via a tight transition state 26.4 kcal/mol above the reactants. The characteristic energy threshold for the reaction under single collision conditions is predicted be in the range of approximately 24 kcal/mol. Product branching ratios obtained by solving kinetic equations with individual rate constants calculated using RRKM and VTST theories for collision energies between 25 and 35 kcal/mol show that l-C5H+H are the dominant reaction products, whereas HC2C3+H and l-C5+H2 are minor products with branching ratios not exceeding 2.5% and 0.7%, respectively. The ethynylcyclopropenylidene isomer i5 is calculated to be the most stable C5H2 species, more favorable than triplet pentadiynylidene i3t by approximately 2 kcal/mol.     

Unimolecular dissociation of the CH3OCO radical: an intermediate in the CH3O + CO reaction

This work investigates the unimolecular dissociation of the methoxycarbonyl, CH(3)OCO, radical. Photolysis of methyl chloroformate at 193 nm produces nascent CH(3)OCO radicals with a distribution of internal energies, determined by the velocities of the momentum-matched Cl atoms, that spans the theoretically predicted barriers to the CH(3)O + CO and CH(3) + CO(2) product channels. Both electronic ground- and excited-state radicals undergo competitive dissociation to both product channels. The experimental product branching to CH(3) + CO(2) from the ground-state radical, about 70%, is orders of magnitude larger than Rice-Ramsperger-Kassel-Marcus (RRKM)-predicted branching, suggesting that previously calculated barriers to the CH(3)OCO --> CH(3) + CO(2) reaction are dramatically in error. Our electronic structure calculations reveal that the cis conformer of the transition state leading to the CH(3) + CO(2) product channel has a much lower barrier than the trans transition state. RRKM calculations using this cis transition state give product branching in agreement with the experimental branching. The data also suggest that our experiments produce a low-lying excited state of the CH(3)OCO radical and give an upper limit to its adiabatic excitation energy of 55 kcal/mol.