Reaction of H + ketene to formyl methyl and acetyl radicals and reverse dissociations

Theoretical investigations on mechanisms and kinetics of methylketene with O(3P) reaction in the atmosphere

The O(³P)-initiated conversion mechanism and dynamics of CH₃CHCO were researched in atmosphere by executing density functional theory (DFT) computations. Optimizations of all the species and single-point energy computations were implemented at the B3LYP/6-311++G(d,p) and CCSD(T)/cc-pVTZ level, respectively. The explicit oxidation mechanism was introduced and discussed. The results state clearly that the O(³P) association was more energetically beneficial than the abstraction of H. The rate coefficients over the probable temperature range of 200-3000 K were forecasted by implementing Rice-Ramsperger-Kassel-Marcus (RRKM) theory. Specifically, the total rate coefficient of O(³P) association reactions is 1.19 × 10^-11 cm³ molecule^-1 s^-1 at 298 K, which is consistent with the experimental results (1.16 × 10^-11 cm³ molecule^-1 s^-1). The rate coefficients for the O(³P) with CH₂CO, CH₃CHCO, and (CH₃)₂CCO suggest that rate coefficient of ketene derivatives increase with the increase of methylation degree. Graphical abstract.

Products and mechanism of thermal decomposition of chlorpyrifos under inert and oxidative conditions

Chlorpyrifos (CPF) is a widely used pesticide; however, limited experimental work has been completed on its thermal decomposition. CPF is known to decompose into 3,5,6-trichloro-2-pyridinol (TCpyol) together with ethylene and HOPOS. Under oxidative conditions TCpyol can decompose into the dioxin-like 2,3,7,8-tetrachloro-[1,4]-dioxinodipyridine (TCDDPy). With CPF on the cusp of being banned in several jurisdictions worldwide, the question might arise as to how to safely eliminate large stockpiles of this pesticide. Thermal methods such as incineration or thermal desorption of pesticide-contaminated soils are often employed. To assess the safety of thermal methods, information about the toxicants arising from thermal treatment is essential. The present flow reactor study reports the products detected under inert and oxidative conditions from the decomposition of CPF representative of thermal treatments and of wildfires in CPF-contaminated vegetation. Ethylene and TCpyol are the initial products formed at temperatures between 550 and 650 °C, although the detection of HOPOS as a reaction product has proven to be elusive. During pyrolysis of CPF in an inert gas, the dominant sulfur-containing product detected from CPF is carbon disulfide. Quantum chemical analysis reveals that ethylene and HOPOS undergo a facile reaction to form thiirane (c-C₂H₄S) which subsequently undergoes ring opening reactions to form precursors of CS₂. At elevated temperatures (>650 °C), TCpyol undergoes both decarbonylation and dehydroxylation reactions together with decomposition of its primary product, TCpyol. A substantial number of toxicants is observed, including HCN and several nitriles, including cyanogen. No CS₂ is observed under oxidative conditions - sulfur dioxide is the fate of S in oxidation of CPF, and quantum chemical studies show that SO₂ formation is initiated by the reaction between HOPOS and O₂. The range of toxicants produced in thermal decomposition of CPF is summarised.

Temperature Dependence of the Reaction of Chlorine Atoms with CH3OH and CH3CHO

Rate constants of the reactions Cl + CH₃OH → CH₂OH + HCl ( k₁) and Cl + CH₃CHO → CH₃C(O) + HCl ( k₃) were measured at 100 Torr over the temperature range 230.3-297.1 K. Radical chemistry was initiated by pulsed laser photolysis of Cl₂ in mixtures of CH₃OH and CH₃CHO in a flow reactor. Heterodyne near-IR wavelength modulation spectroscopy was used to directly detect HO₂ produced from the subsequent reaction of CH₂OH with O₂ in real time to determine the rate of reaction of Cl with CH₃OH. The rate of Cl + CH₃CHO was measured relative to that of the Cl + CH₃OH reaction. Secondary chemistry, including that of the adducts HO₂·CH₃OH and HO₂·CH₃CHO, was taken into account. The Arrhenius expressions were found to be k₁( T) = 5.02_-1.5^+1.8 × 10^-11 exp[(20 ± 88)/ T] cm³ molecule^-1 s^-1 and k₃( T) = 6.38_-2.0^+2.4 × 10^-11 exp[(56 ± 90)/ T] cm³ molecule^-1 s^-1 (2σ uncertainties). The average values of the rate constants over this temperature range were k₁ = (5.45 ± 0.37) × 10^-11 cm³ molecule^-1 s^-1 and k₃ = (8.00 ± 1.27) × 10^-11 cm³ molecule^-1 s^-1 (2σ uncertainties), consistent with current literature values.

Mechanism and rate constants of the CH2 + CH2 CO reactions in triplet and singlet states: A theoretical study

Ab initio and density functional CCSD(T)-F12/cc-pVQZ-f12//B2PLYPD3/6-311G** calculations have been performed to unravel the reaction mechanism of triplet and singlet methylene CH₂ with ketene CH₂ CO. The computed potential energy diagrams and molecular properties have been then utilized in Rice-Ramsperger-Kassel-Marcus-Master Equation (RRKM-ME) calculations of the reaction rate constants and product branching ratios combined with the use of nonadiabatic transition state theory for spin-forbidden triplet-singlet isomerization. The results indicate that the most important channels of the reaction of ketene with triplet methylene lead to the formation of the HCCO + CH₃ and C₂ H₄ + CO products, where the former channel is preferable at higher temperatures from 1000 K and above. In the C₂ H₄ + CO product pair, the ethylene molecule can be formed either adiabatically in the triplet electronic state or via triplet-singlet intersystem crossing in the singlet electronic state occurring in the vicinity of the CH₂ COCH₂ intermediate or along the pathway of CO elimination from the initial CH₂ CH₂ CO complex. The predominant products of the reaction of ketene with singlet methylene have been shown to be C₂ H₄ + CO. The formation of these products mostly proceeds via a well-skipping mechanism but at high pressures may to some extent involve collisional stabilization of the CH₃ CHCO and cyclic CH₂ COCH₂ intermediates followed by their thermal unimolecular decomposition. The calculated rate constants at different pressures from 0.01 to 100 atm have been fitted by the modified Arrhenius expressions in the temperature range of 300-3000 K, which are proposed for kinetic modeling of ketene reactions in combustion. © 2018 Wiley Periodicals, Inc.

Comparing electronic structure predictions for the ground state dissociation of vinoxy radicals

This paper reports a series of electronic structure calculations performed on the dissociation pathways of the vinoxy radical (CH(2)CHO). We use coupled-cluster with single, double, and perturbative triple excitations (CCSD(T)), complete active space self-consistent field (CASSCF), multireference configuration interaction (MRCI), and MRCI with the Davidson correction (MRCI+Q) to calculate the barrier heights of the two unimolecular dissociation pathways of this radical. The effect of state averaging on the barrier heights is investigated at the CASSCF, MRCI, and MRCI+Q levels. The change in mixing angle along the reaction path is calculated as a measure of derivative coupling and found to be insufficient to suggest nonadiabatic recrossing. We also present a new analysis of previous experimental data on the unimolecular dissociation of ground state vinoxy. In particular, an error in the internal energy distribution of vinoxy radicals reported in a previous paper is corrected and a new analysis of the experimental sensitivity to the onset energy (barrier height) for the isomerization reaction is given. Combining these studies, a final "worst case" analysis of the product branching ratio is given and a statistical model using each of the calculated transition states is found to be unable to correctly reproduce the experimental data.

Toluene Combustion: Reaction Paths, Thermochemical Properties, and Kinetic Analysis for the Methylphenyl Radical + O2 Reaction

Aromatic compounds such as toluene and xylene are major components of many fuels. Accurate kinetic mechanisms for the combustion of toluene are, however, incomplete, as they do not accurately model experimental results such as strain rates and ignition times and consistently underpredict conversion. Current kinetic mechanisms for toluene combustion neglect the reactions of the methylphenyl radicals, and we believe that this is responsible, in part, for the shortcomings of these models. We also demonstrate how methylphenyl radical formation is important in the combustion and pyrolysis of other alkyl-substituted aromatic compounds such as xylene and trimethylbenzene. We have studied the oxidation reactions of the methylphenyl radicals with O2 using computational ab initio and density functional theory methods. A detailed reaction submechanism is presented for the 2-methylphenyl radical + O2 system, with 16 intermediates and products. For each species, enthalpies of formation are calculated using the computational methods G3 and G3B3, with isodesmic work reactions used to minimize computational errors. Transition states are calculated at the G3B3 level, yielding high-pressure limit elementary rate constants as a function of temperature. For the barrierless methylphenyl + O2 and methylphenoxy + O association reactions, rate constants are determined from variational transition state theory. Multichannel, multifrequency quantum Rice-Ramsperger-Kassel (qRRK) theory, with master equation analysis for falloff, provides rate constants as a function of temperature and pressure from 800 to 2400 K and 1 x 10(-4) to 1 x 10(3) atm. Analysis of our results shows that the dominant pathways for reaction of the three isomeric methylphenyl radicals is formation of methyloxepinoxy radicals and subsequent ring opening to methyl-dioxo-hexadienyl radicals. The next most important reaction pathway involves formation of methylphenoxy radicals + O in a chain branching process. At lower temperatures, the formation of stabilized methylphenylperoxy radicals becomes significant. A further important reaction channel is available only to the 2-methylphenyl isomer, where 6-methylene-2,4-cyclohexadiene-1-one (ortho-quinone methide, o-QM) is produced via an intramolecular hydrogen transfer from the methyl group to the peroxy radical in 2-methylphenylperoxy, with subsequent loss of OH. The decomposition of o-QM to benzene + CO reveals a potentially important new pathway for the conversion of toluene to benzene during combustion. A number of the important products of toluene combustion proposed in this study are known to be precursors of polyaromatic hydrocarbons that are involved in soot formation. Reactions leading to the important unsaturated oxygenated intermediates identified in this study, and the further reactions of these intermediates, are not included in current aromatic oxidation mechanisms.

Enthalpies of Formation, Bond Dissociation Energies, and Molecular Structures of the n-Aldehydes (Acetaldehyde, Propanal, Butanal, Pentanal, Hexanal, and Heptanal) and Their Radicals

Aldehydes are important intermediates and products in a variety of combustion and gas-phase oxidation processes, such as in low-temperature combustion, in the atmosphere, and in interstellar media. Despite their importance, the enthalpies of formation and bond dissociation energies (BDEs) for the aldehydes are not accurately known. We have determined enthalpies of formation for acetaldehyde, propanal, and butanal from thermodynamic cycles, using experimentally measured reaction and formation enthalpies. All enthalpy values used for reference molecules and reactions were first verified to be accurate to within around 1 kcal mol-1 using high-level ab initio calculations. Enthalpies of formation were found to be -39.72 +/- 0.16 kcal mol-1 for acetaldehyde, -45.18 +/- 1.1 kcal mol-1 for propanal, and -49.27 +/- 0.16 kcal mol-1 for butanal. Enthalpies of formation for these three aldehydes, as well as for pentanal, hexanal, and heptanal, were calculated using the G3, G3B3, and CBS-APNO theoretical methods, in conjunction with bond-isodesmic work reactions. On the basis of the results of our thermodynamic cycles, theoretical calculations using isodesmic work reactions, and existing experimental measurements, we suggest that the best available formation enthalpies for the aldehydes acetaldehyde, propanal, butanal, pentanal, hexanal, and heptanal are -39.72, -45.18, -50.0, -54.61, -59.37, and -64.2 kcal mol-1, respectively. Our calculations also identify that the literature enthalpy of formation of crotonaldehyde is in error by as much as 1 kcal mol-1, and we suggest a value of -25.1 kcal mol-1, which we calculate using isodesmic work reactions. Bond energies for each of the bonds in the aldehydes up to pentanal were calculated at the CBS-APNO level. Analysis of the BDEs reveals the R-CH(2)CH=O to be the weakest bond in all aldehydes larger than acetaldehyde, due to formation of the resonantly stabilized vinoxy radical (vinyloxy radical/formyl methyl radical). It is proposed that the vinoxy radical as well as the more commonly considered formyl and acetyl radicals are important products of aldehyde combustion and oxidation, and the reaction pathways of the vinoxy, formyl, and acetyl radicals are discussed. Group additivity values for the carbon-oxygen-hydrogen groups common to the aldehydes are also determined. Internal rotor profiles and electrostatic potential surfaces are used to study the dipole induced dipole-dipole interaction in the synperiplanar conformation of propanal. It is proposed that the loss of this dipole-dipole interaction in RC(.-)HCH(2)CH=O radicals causes a ca. 1-2 kcal mol-1 decrease in the aldehyde C-H and C-C bond energies corresponding to RC(.-)HCH(2)CH=O radical formation.

Thermodynamic Properties (Enthalpy, Bond Energy, Entropy, and Heat Capacity) and Internal Rotor Potentials of Vinyl Alcohol, Methyl Vinyl Ether, and Their Corresponding Radicals

Vinyl alcohols (enols) have been discovered as important intermediates and products in the oxidation and combustion of hydrocarbons, while methyl vinyl ethers are also thought to occur as important combustion intermediates. Vinyl alcohol has been detected in interstellar media, while poly(vinyl alcohol) and poly(methyl vinyl ether) are common polymers. The thermochemical property data on these vinyl alcohols and methyl vinyl ethers is important for understanding their stability, reaction paths, and kinetics in atmospheric and thermal hydrocarbon-oxygen systems. Enthalpies , entropies , and heat capacities (C(p)()(T)) are determined for CH(2)=CHOH, C(*)H=CHOH, CH(2)=C(*)OH, CH(2)=CHOCH(3), C(*)H=CHOCH(3), CH(2)=C(*)OCH(3), and CH(2)=CHOC(*)H(2). Molecular structures, vibrational frequencies, , and C(p)(T) are calculated at the B3LYP/6-31G(d,p) density functional calculation level. Enthalpies are also determined using the composite CBS-Q, CBS-APNO, and G3 methods using isodesmic work reactions to minimize calculation errors. Potential barriers for internal rotors are calculated at the B3LYP/6-31G(d,p) level and used to determine the hindered internal rotational contributions to entropy and heat capacity. The recommended ideal gas phase values calculated in this study are the following (in kcal mol(-1)): -30.0, -28.9 (syn, anti) for CH(2)=CHOH; -25.6, -23.9 for CH(2)=CHOCH(3); 31.3, 33.5 for C(*)H=CHOH; 27.1 for anti-CH(2)=C(*)OH; 35.6, 39.3 for C(*)H=CHOCH(3); 33.5, 32.2 for CH(2)=C(*)OCH(3); 21.3, 22.0 for CH(2)=CHOC(*)H(2). Bond dissociation energies (BDEs) and group additivity contributions are also determined. The BDEs reveal that the O-H, O-CH(3), C-OH, and C-OCH(3) bonds in vinyl alcohol and methyl vinyl ether are similar in energy to those in the aromatic molecules phenol and methyl phenyl ether, being on average around 3 kcal mol(-1) weaker in the vinyl systems. The keto-enol tautomerization enthalpy for the interconversion of vinyl alcohol to acetaldehyde is determined to be -9.7 kcal mol(-1), while the activation energy for this reaction is calculated as 55.9 kcal mol(-1); this is the simplest keto-enol tautomerization and is thought to be important in the reactions of vinyl alcohol. Formation of the formyl methyl radical (vinoxy radical/vinyloxy radical) from both vinyl alcohol and methyl vinyl ether is also shown to be important, and its reactions are discussed briefly.

A novel conical intersection topography and its consequences: The 1, 2A2 conical intersection seam of the vinoxy radical

A region of the 1, 2 2A seam of accidental conical intersections in the vinoxy radical exhibits a novel topography which has important consequences for both upper-state to lower-state and lower-state to upper-state nonadiabatic transitions. The impact of this topography on these nonadiabatic transitions is described. We also considered the possibility that this conical intersection seam provides a dynamical bottleneck to the photodissociation of vinoxy to H+ketene by facilitating nonadiabatic recrossing. Our analysis of the conical topographies and the proximity of the conical intersections to the transition state for dissociation to H+ketene does not support nonadiabatic recrossing as an effective dynamical bottleneck blocking the H+ketene channel.

Thermal Decomposition of Ethylene Oxide: Potential Energy Surface, Master Equation Analysis, and Detailed Kinetic Modeling

The unimolecular decomposition of ethylene oxide (oxirane) and the oxiranyl radial is examined by molecular orbital calculations, Rice-Ramsperger-Kassel-Marcus (RRKM)/Master Equation analysis, and detailed kinetic modeling of ethylene oxide pyrolysis in a single-pulse shock tube. It was found that the largest energy barrier to the decomposition of ethylene oxide lies in its initial isomerization to form acetaldehyde, and in agreement with previous studies, the isomerization was found to proceed through the *CH2CH2O* biradical. Because of the biradical nature of the transition states and intermediate, the energy barriers for the initial C-O rupture in ethylene oxide and the subsequent 1,2-H shift remain highly uncertain. An overall isomerization energy barrier of 59 +/- 2 kcal/mol was found to satisfactorily explain the available single pulse shock tube data. This barrier height is in line with the estimates made from an approximate spin-corrected procedure at the MP4/6-31+G(d) and QCISD(T)/6-31G(d) levels of theory. The dominant channel for the unimolecular decomposition of ethylene oxide was found to form CH3 + HCO at around the ambient pressure. It accounts for >90% of the total rate constant for T > 800 K. The high-pressure limit rate constant for the unimolecular decomposition of ethylene oxide was calculated as k(1,infinity)(s(-1)) = (3.74 x 10(10))T(1.298)e(-29990/T) for 600 < T < 2000 K.