ab initio molecular orbital and density functional analysis of acetylene + O2 reactions with CHEMKIN evaluation

A number of researchers have indicated that a direct reaction of acetylene with oxygen needs to be included in detailed reaction mechanisms in order to model observed flame speeds and induction times. Four pathways for the initiation of acetylene oxidation to chain propagation are considered and the rate constants are compared with values used in the mechanisms:

• 1 ³O₂ + HCCH to triplet adduct and reaction on the triplet surface
• 2 ³O₂ + HCCH to triplet adduct, conversion of triplet adduct to singlet adduct via collision in the reaction environment, with further reaction of the singlet adduct
• 3 ¹O₂ + HCCH to singlet adduct
• 4 Isomerization of HCCH to vinylidene and then vinylidene insertion reaction with ³O₂

Elementary reaction pathways for oxidation of acetylene by addition reaction of O₂(³Σ) on the triplet surface are analyzed. ab initio molecular orbital and density functional calculations are employed to estimate the thermodynamic properties of the reactants, transition states, and products in this system. Acetylene oxidation reaction over the triplet surface is initiated by addition of molecular oxygen, O₂(³Σ), to a carbon atom, forming a triplet peroxy‐ethylene biradical. The reaction path to major products, either two formyl radicals or glyoxal radical plus hydrogen atom, involves reaction through three transition states: O₂(³Σ) addition to acetylene (TS1), peroxy radical addition at the ipso‐carbon to form a dioxirane (TS2), and cleavage of O? O bond in a three‐member ring (TS3). Single‐point QCISD(T) and B3LYP calculations with large basis sets were performed to try to verify barrier heights on important transition states. A second pathway to product formation is through spin conversion of the triplet peroxy‐ethylene biradical to the singlet by collision with bath gas. Rapid ring closure of the singlet peroxy‐ethylene biradical to form a four‐member ring is followed by breaking of the peroxy bond to form glyoxal, which further dissociates to either two formyl radicals or a glyoxal radical plus hydrogen atom. The overall forward rate constant through this pathway is estimated to be k_f = 2.21 × 10⁷ T^1.46e^{?33.1(kcal/mol)/RT}. Two additional pathways from the literature, HCCH + O₂(¹Δ) and pressure‐dependent isomerization of acetylene to vinylidene and then vinylidene reaction with O₂(³Σ), are also evaluated for completeness. CHEMKIN modeling on each of the four proposed pathways is performed and concentration profiles from these reactions are evaluated at 0.013 atm and 1 atm over 35 milliseconds. Through reaction on the triplet surface is evaluated to be not important. Formation of the triplet adduct with conversion (via collision) to a singlet and the vinylidene paths show similar and lower rates than those used in mechanisms, respectively. Our implementation of the HCCH + O₂(¹Δ) pathway of Benson suggests the need to include: (i) reverse reaction, (ii) barriers to further reaction of the initial adduct plus (iii) further evaluation of the O₂(¹Δ) addition barrier. The pathways from triplet adduct with conversion to singlet and from vinylidene are both recommended for initiation of acetylene oxidation. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 623–641, 2000