Mechanistic and kinetic study of the CH3CO+O2 reaction

The fundamental vibrational frequencies and spectroscopic constants of the C2O2H2 isomers: molecules known in simulated interstellar ice analogues

While trans-glyoxal may not be easily observable in astronomical sources through either IR or radioastronomy due to its C2h symmetry, its cis conformer along with the cyc-H2COCO epoxide isomer should be ready targets for astrochemical detection. The present quantum chemical study shows that not only are both molecular isomers strongly polar, they also have notable IR features and low isomerisation energies of 4.1 kcal mol-1 and 10.7 kcal mol-1, respectively. These three isomers along with two other C2O2H2 isomers have had their full set of fundamental vibrational frequencies and spectroscopic constants characterised herein. These isomers have previously been shown to occur in simulated astrophysical ices making them worthy targets of astronomical search. Furthermore, the hybrid quartic force field (QFF) approach utilized herein to produce the needed spectral data has a mean absolute percent error compared to the experimentally-available, gas phase fundamental vibrational frequencies of 0.6% and rotational constants to better than 0.1%. The hybrid QFF is defined from explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples level [CCSD(T)-F12b] including core electron correlation and a canonical CCSD(T) relativity correction for the harmonic (quadratic) terms in the QFF and simple CCSD(T)-F12b/cc-pVDZ energies for the cubic and quartic terms, the so-called "F12-TcCR+DZ QFF." This method is producing spectroscopically-accurate predictions for both fundamental vibrational frequencies and principal spectroscopic constants. Hence, the values computed in this work should be notably accurate and, hence, exceptionally useful to the spectroscopy and astrochemistry communities.

Comprehensive Theoretical Study on Four Typical Intramolecular Hydrogen Shift Reactions of Peroxy Radicals: Multireference Character, Recommended Model Chemistry, and Kinetics

Intramolecular hydrogen shift reactions in peroxy radicals (RO₂^• → ^•QOOH) play key roles in the low-temperature combustion and in the atmospheric chemistry. In the present study, we found that a mild-to-moderate multireference character of a potential energy surface (PES) is widely present in four typical hydrogen shift reactions of peroxy radicals (RO₂^•, R = ethyl, vinyl, formyl methyl, and acetyl) by a systematic assessment based on the T₁ diagnostic, %TAE diagnostic, M diagnostic, and contribution of the dominant configuration of the reference CASSCF wavefunction (C₀²). To assess the effects of these inherent multireference characters on electronic structure calculations, we compared the PESs of the four reactions calculated by the multireference method CASPT2 in the complete basis set (CBS) limit, single-reference method CCSD(T)-F12, and single-reference-based composite method WMS. The results showed that ignoring the multireference character will introduce a mean unsigned deviation (MUD) of 0.46-1.72 kcal/mol from CASPT2/CBS results by using the CCSD(T)-F12 method or a MUD of 0.49-1.37 kcal/mol by WMS for three RO₂^• reactions (R = vinyl, formyl methyl, and acetyl) with a stronger multireference character. Further tests by single-reference Kohn-Sham (KS) density functional theory methods showed even larger deviations. Therefore, we specifically developed a new hybrid meta-generalized gradient approximation (GGA) functional M06-HS for the four typical H-shift reactions of peroxy radicals based on the WMS results for the ethyl peroxy radical reaction and on the CASPT2/CBS results for the others. The M06-HS method has an averaged MUD of 0.34 kcal/mol over five tested basis sets against the benchmark PESs, performing best in the tested 38 KS functionals. Last, in a temperature range of 200-3000 K, with the new functional, we calculated the high-pressure-limit rate coefficients of these H-shift reactions by the multi-structural variational transition-state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and the thermochemical properties of all of the involved key radicals by the multi-structural torsional (MS-T) anharmonicity approximation method.

Unimolecular decomposition of acetyl peroxy radical: a potential source of tropospheric ketene

The unimolecular decomposition of acetyl peroxy radicals followed by subsequent nitration is known to lead to the formation of peroxy acetyl nitrate (PAN) in the troposphere. Using high level quantum chemical calculations, we show that the acetyl peroxy radical is a precursor in the formation of tropospheric ketene. The results show that the presence of a single or double water molecule(s) as a catalyst does not influence the decomposition reaction directly to form ketene and hydroperoxy radicals. The electronic excitation of the reactive and product complexes occurs in the wavelength range of ∼1400 nm, suggesting that the complexes undergo photoexcitation in the near IR region. The results ascertain that the dissociation of acetyl peroxy radicals into ketene and hydroperoxy radicals occurs more likely through the excitation route and the corresponding excitation wavelength reveals that the reactions are red-light driven. Three different product complexes, ketene·HO2, ketene·H2O·HO2 and ketene·(H2O)2·HO2, are formed from the reaction. The direct dynamics simulations show that the product complexes are more stable and possess a long lifetime. The calculated temperature dependent equilibrium constant of the product complexes reveals that their atmospheric abundances decrease with increasing altitudes.

Experimental Evidence of Dioxole Unimolecular Decay Pathway for Isoprene-Derived Criegee Intermediates

Ozonolysis of isoprene, one of the most abundant volatile organic compounds emitted into the Earth's atmosphere, generates two four-carbon unsaturated Criegee intermediates, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide). The extended conjugation between the vinyl substituent and carbonyl oxide groups of these Criegee intermediates facilitates rapid electrocyclic ring closures that form five-membered cyclic peroxides, known as dioxoles. This study reports the first experimental evidence of this novel decay pathway, which is predicted to be the dominant atmospheric sink for specific conformational forms of MVK-oxide (anti) and MACR-oxide (syn) with the vinyl substituent adjacent to the terminal O atom. The resulting dioxoles are predicted to undergo rapid unimolecular decay to oxygenated hydrocarbon radical products, including acetyl, vinoxy, formyl, and 2-methylvinoxy radicals. In the presence of O₂, these radicals rapidly react to form peroxy radicals (ROO), which quickly decay via carbon-centered radical intermediates (QOOH) to stable carbonyl products that were identified in this work. The carbonyl products were detected under thermal conditions (298 K, 10 Torr He) using multiplexed photoionization mass spectrometry (MPIMS). The main products (and associated relative abundances) originating from unimolecular decay of anti-MVK-oxide and subsequent reaction with O₂ are formaldehyde (88 ± 5%), ketene (9 ± 1%), and glyoxal (3 ± 1%). Those identified from the unimolecular decay of syn-MACR-oxide and subsequent reaction with O₂ are acetaldehyde (37 ± 7%), vinyl alcohol (9 ± 1%), methylketene (2 ± 1%), and acrolein (52 ± 5%). In addition to the stable carbonyl products, the secondary peroxy chemistry also generates OH or HO₂ radical coproducts.

THEORETICAL STUDIES ON PHOTOIONIZATION OF GUANINE TAUTOMERS AND INTERCONVERSION OF CATION RADICALS

The G3MP2B3 and P3 methods have been used to calculate the adiabatic and vertical ionization potentials (IPs) of the eight most stable tautomers of guanine. The calculated energy discrepancy between adiabatic and vertical IPs are in good agreement with the changes in geometry from neutral ground state to stable cation radicals. The geometries of imino-oxo form tautomers have no obvious change in the ionization process, which results in less energy discrepancy between vertical and adiabatic IPs. In the ionization process, the geometries of the amino-oxo and amino-hydroxy form tautomers change from nonplanar to planar structures. Hence the amino-oxo and amino-hydroxy form tautomers have larger energy discrepancy between vertical and adiabatic IPs. Further studies on the interconversion of the cation radicals shed further light on the transition process between the cation radicals and the main pathways are the hydrogen migrations and internal rotations of hydroxy (OH) and imino (NH) groups. The barriers of hydrogen rotations are lower than those of hydrogen migrations. Furthermore, the barriers for the hydrogen migrations between two rings are higher, which are about 3.0 eV.

Developments in Laboratory Studies of Gas-Phase Reactions for Atmospheric Chemistry with Applications to Isoprene Oxidation and Carbonyl Chemistry

Laboratory studies of gas-phase chemical processes are a key tool in understanding the chemistry of our atmosphere and hence tackling issues such as climate change and air quality. Laboratory techniques have improved considerably with greater emphasis on product detection, allowing the measurement of site-specific rate coefficients. Radical chemistry lies at the heart of atmospheric chemistry. In this review we consider issues around radical generation and recycling from the oxidation of isoprene and from the chemical reactions and photolysis of carbonyl species. Isoprene is the most globally significant hydrocarbon, but uncertainties exist about its oxidation in unpolluted environments. Recent experiments and calculations that cast light on radical generation are reviewed. Carbonyl compounds are the dominant first-generation products from hydrocarbon oxidation. Chemical oxidation can recycle radicals, or photolysis can be a net radical source. Studies have demonstrated that high-resolution and temperature-dependent studies are important for some significant species.

Mechanistic and kinetic study of CH2O+O3 reaction

Both singlet and triplet potential energy surfaces for the reaction of ground-state formaldehyde (CH(2)O) and ozone (O(3)) are theoretically investigated at the BMC-CCSD//BHandHLYP/6-311+G(d,p) level. Various possible isomerization and dissociation pathways are probed. Hydrogen abstraction, oxygen abstraction, and C-addition/elimination are found on both the singlet and the triplet surfaces. The major products for the total reaction are HCO and HOOO, which are generated via hydrogen abstraction. The transition state theory (TST) and multichannel RRKM calculations have been carried out for the total and individual rate constants for determinant channels over a wide range of temperatures and pressures.

Transient infrared absorption of t-CH3C(O)OO, c-CH3C(O)OO, and alpha-lactone recorded in gaseous reactions of CH3CO and O2

A step-scan Fourier-transform infrared spectrometer coupled with a multipass absorption cell was utilized to monitor the transient species produced in gaseous reactions of CH(3)CO and O(2); IR absorption spectra of CH(3)C(O)OO and alpha-lactone were observed. Absorption bands with origins at 1851+/-1, 1372+/-2, 1169+/-6, and 1102+/-3 cm(-1) are attributed to t-CH(3)C(O)OO, and those at 1862+/-3, 1142+/-4, and 1078+/-6 cm(-1) are assigned to c-CH(3)C(O)OO. A weak band near 1960 cm(-1) is assigned to alpha-lactone, cyc-CH(2)C(=O)O, a coproduct of OH. These observed rotational contours agree satisfactorily with simulated bands based on predicted rotational parameters and dipole derivatives, and observed vibrational wavenumbers agree with harmonic vibrational wavenumbers predicted with B3LYP/aug-cc-pVDZ density-functional theory. The observed relative intensities indicate that t-CH(3)C(O)OO is more stable than c-CH(3)C(O)OO by 3+/-2 kJ mol(-1). Based on these observations, the branching ratio for the OH+alpha-lactone channel of the CH(3)CO+O(2) reaction is estimated to be 0.04+/-0.01 under 100 Torr of O(2) at 298 K. A simple kinetic model is employed to account for the decay of CH(3)C(O)OO.

Ab initio study of the reaction of propionyl (C2H5CO) radical with oxygen (O2)

The reaction of propionyl radical with oxygen has been studied using the full coupled cluster theory with the complete basis set. This is the first time to gain a conclusive insight into the reaction mechanism and kinetics for this important reaction in detail. The reaction takes place via a chemical activation mechanism. The barrierless association of propionyl with oxygen produces the propionylperoxy radical, which decomposes to form the hydroxyl radical and the three-center alpha-lactone predominantly or the four-center beta-propiolactone. The oxidation of propionyl radical to carbon monoxide or carbon dioxide is not straightforward rather via the secondary decomposition of alpha-lactone and beta-propiolactone. Kinetically, the overall rate constant is almost pressure independent and it approaches the high-pressure limit around tens of torr of helium. At temperatures below 600 K, the rate constant shows negative temperature dependence. The experimental yields of the hydroxyl radical can be well reproduced, with the average energy transferred per collision -DeltaE=20-25 cm(-1) at 213 and 295 K (helium bath gas). At low pressures, together with the hydroxy radical, alpha-lactone is the major product, while beta-propiolactone only accounts for about one-fifth of alpha-lactone. At the high-pressure limit, the production of the propionylperoxy radical is dominant together with a fraction of the isomers. The infrared spectroscopy or the mass spectroscopy techniques are suggested to be employed in the future experimental study of the C2H5CO+O2 reaction.

Kinetics and mechanism of the reactions of CH3CO and CH3C(O)CH2 radicals with O2. Low-pressure discharge flow experiments and quantum chemical computations

The reactions CH(3)CO + O(2)--> products (1), CH(3)CO + O(2)--> OH +other products (1b) and CH(3)C(O)CH(2) + O(2)--> products (2) have been studied in isothermal discharge flow reactors with laser induced fluorescence monitoring of OH and CH(3)C(O)CH(2) radicals. The experiments have been performed at overall pressures between 1.33 and 10.91 mbar of helium and 298 +/- 1 K reaction temperature. OH formation has been found to be the dominant reaction channel for CH(3)CO + O(2): the branching ratio, Gamma(1b) = k(1b)/k(1), is close to unity at around 1 mbar, but decreases rapidly with increasing pressure. The rate constant of the overall reaction, k(2), has been found to be pressure dependent: the fall-off behaviour has been analysed in comparison with reported data. Electronic structure calculations have confirmed that at room temperature the reaction of CH(3)C(O)CH(2) with O(2) is essentially a recombination-type process. At high temperatures, the further reactions of the acetonyl-peroxyl adduct may yield OH radicals, but the most probable channel seems to be the O(2)-catalysed keto-enol transformation of acetonyl. Implications of the results for atmospheric modelling studies have been discussed.

Master equation simulations of competing unimolecular and bimolecular reactions: application to OH production in the reaction of acetyl radical with O2

Master equation calculations were carried out to simulate the production of hydroxyl free radicals initiated by the reaction of acetyl free radicals (CH3(C=O).) with molecular oxygen. In particular, the competition between the unimolecular reactions and bimolecular reactions of vibrationally excited intermediates was modeled by using a single master equation. The vibrationally excited intermediates (isomers of acetylperoxyl radicals) result from the initial reaction of acetyl free radical with O2. The bimolecular reactions were modeled using a novel pseudo-first-order microcanonical rate constant approach. Stationary points on the multi-well, multi-channel potential energy surface (PES) were calculated at the DFT(B3LYP)/6-311G(2df,p) level of theory. Some additional calculations were carried out at the CASPT2(7,5)/6-31G(d) level of theory to investigate barrierless reactions and other features of the PES. The master equation simulations are in excellent agreement with the experimental OH yields measured in N2 or He buffer gas near 300 K, but they do not explain a recent report that the OH yields are independent of pressure in nearly pure O2 buffer gas.