Computational study of the reaction mechanism of the methylperoxy self-reaction

To provide insight on the reaction mechanism of the methyperoxy (CH(3)O(2)•) self-reaction, stationary points on both the spin-singlet and the spin-triplet potential energy surfaces of 2(CH(3)O(2)•) have been searched at the B3LYP/6-311++G(2df,2p) level. The relative energies, enthalpies, and free energies of these stationary points are calculated using CCSD(T)/cc-pVTZ. Our theoretical results indicate that reactions on a spin-triplet potential energy surface are kinetically unfavorable due to high free energy barriers, while they are more complicated on the spin-singlet surface. CH(3)OOCH(3) + O(2)(1) can be produced directly from 2(CH(3)O(2)•), while in other channels, three spin-singlet chain-structure intermediates are first formed and subsequently dissociated to produce different products. Besides the dominant channels producing 2CH(3)O• + O(2) and CH(3)OH + CH(2)O + O(2) as determined before, the channels leading to CH(3)OOOH + CH(2)O and CH(3)O• + CH(2)O + HO(2)• are also energetically favorable in the self-reaction of CH(3)O(2)• especially at low temperature according to our results.

Near-infrared kinetic spectroscopy of the HO2 and C2H5O2 self-reactions and cross reactions

The self-reactions and cross reactions of the peroxy radicals C2H5O2 and HO2 were monitored using simultaneous independent spectroscopic probes to observe each radical species. Wavelength modulation (WM) near-infrared (NIR) spectroscopy was used to detect HO2, and UV absorption monitored C2H5O2. The temperature dependences of these reactions were investigated over a range of interest to tropospheric chemistry, 221-296 K. The Arrhenius expression determined for the cross reaction, k2(T) = (6.01(-1.47)(+1.95)) x 10(-13) exp((638 +/- 73)/T) cm3 molecules(-1) s(-1) is in agreement with other work from the literature. The measurements of the HO2 self-reaction agreed with previous work from this lab and were not further refined. The C2H5O2 self-reaction is complicated by secondary production of HO2. This experiment performed the first direct measurement of the self-reaction rate constant, as well as the branching fraction to the radical channel, in part by measurement of the secondary HO2. The Arrhenius expression for the self-reaction rate constant is k3(T) = (1.29(-0.27)(+0.34)) x 10(-13)exp((-23 +/- 61)/T) cm3 molecules(-1) s(-1), and the branching fraction value is alpha = 0.28 +/- 0.06, independent of temperature. These values are in disagreement with previous measurements based on end product studies of the branching fraction. The results suggest that better characterization of the products from RO2 self-reactions are required.

High-resolution cavity ringdown spectroscopy of the jet-cooled ethyl peroxy radical C2H5O2

We have recorded high resolution, partially rotationally resolved, jet-cooled cavity ringdown spectra of the origin band of the A-X electronic transition of both the G and T conformers of the perproteo and perdeutero isotopologues of the ethyl peroxy radical, C(2)H(5)O(2). This transition, located in the near infrared, was studied using a narrow band laser source (< or approximately 250 MHz) and a supersonic slit-jet expansion coupled with an electric discharge allowing us to obtain rotational temperatures of about 15 K. All four spectra have been successfully simulated using an evolutionary algorithm approach with a Hamiltonian including rotational and spin-rotational terms. Excellent agreement with the experimental spectra was obtained by fitting seven molecular parameters in each ground and the first excited electronic states as well as the band origin of the electronic transition. This analysis unambiguously confirms the assignment of the lower frequency origin band to the G conformer and the higher frequency one to the T conformer.

Effect of halogenation on the mechanism of the atmospheric reactions between methylperoxy radicals and NO. A computational study

The mechanism of the reactions between the halogenated methylperoxy radicals, CHX(2)O(2) (X = F, Cl), and NO is investigated by using ab initio and density functional quantum mechanical methods. Comparison is made with the mechanism of the CH(3)O(2) + NO reaction. The most important energy minima in the potential energy surface are found to be the two conformers of the halogenated methyl peroxynitrite association adducts, CHX(2)OONOcp and CHX(2)OONOtp, and the halogenated methyl nitrates, CHX(2)ONO(2). The latter are suggested to be formed through the one-step isomerization of the peroxynitrite adduct and may lead upon decomposition to carbonylated species, CX(2)O + HONO and CHXO + XNO(2). The ambiguous issue of the unimolecular peroxynitrite to nitrate isomerization is reconsidered, and the possibility of a triplet transition state involvement in the ROONOtp <--> RONO(2) rearrangement is examined. The overall calculations and the detailed correlation with the methyl system show the significant effect of the halogenation on the lowering of the entrance potential energy well which corresponds to the formation of the peroxynitrites. The increased attractive character of the potential energy surface found upon halogenation combined with the increased exothermicity of the CHX(2)O(2) + NO --> CHX(2)O + NO(2) reaction are suggested to be the important factors contributing to the enhanced reactivity of the halogenated reactions relative to CH(3)O(2) + NO. The calculated heat of formation values indicate the large stabilization of the fluorinated derivatives.

Formation of volatile organic compounds in the heterogeneous oxidation of condensed-phase organic films by gas-phase OH

The yield of volatile organic compounds (VOCs) from the heterogeneous oxidation of condensed-phase organic and hydrocarbon soot films by gas-phase OH has been studied in a coated-wall flow tube at room temperature. Simultaneously, OH concentrations are measured using a chemical-ionization mass spectrometer (CIMS) operated in negative ion mode and VOCs are measured using a commercial proton-transfer-reaction mass spectrometer (PTR-MS). It is observed that a variety of aldehydes/carbonyls and carboxylic acids are formed. Specifically, detailed experiments were conducted with stearic acid, where products are observed that contain as many as 13 carbon atoms with the average carbon number of the products between 3 and 5. The yield of VOCs, relative to the loss of OH radicals, is strongly dependent on the partial pressure of O2 in the carrier gas, ranging from 0.08 +/- 0.03 in a nominally pure He carrier gas to 0.34 +/- 0.14 in 6 Torr of pure O2. Yields from other organics are somewhat lower than those from stearic acid, ranging in conditions of pure O2 from 0.10 +/- 0.04 for BES (bis(ethylhexyl)sebacate), to 0.03 +/- 0.01 for n-hexane soot, to 0.01 +/- 0.005 for pyrene. Under atmospheric conditions, OH oxidation of select organics may be an efficient source of small VOCs. In particular, formic acid is formed in significant yield from all the surfaces.

Atmospheric chemistry of 2-ethoxy-3,3,4,4,5-pentafluorotetrahydro-2,5-bis[1,2,2,2-tetrafluoro-1-(trifluoromethyl)ethyl]-furan: kinetics, mechanisms, and products of Cl atom and OH radical initiated oxidation

Smog chamber/FTIR techniques were used to study the atmospheric chemistry of the title compound which we refer to as RfOC2H5. Rate constants of k(Cl + RfOC2H5) = (2.70 +/- 0.36) x 10(-12), k(OH + RfOC2H5) = (5.93 +/- 0.85) x 10(-14), and k(Cl + RfOCHO) = (1.34 +/- 0.20) x 10(-14) cm3 molecule(-1') s(-1) were measured in 700 Torr of N2, or air, diluent at 294 +/- 1 K. From the value of k(OH + RfOC2H5) the atmospheric lifetime of RfOC2H5 was estimated to be 1 year. Two competing loss mechanisms for RfOCH(O*)CH3 radicals were identified in 700 Torr of N2/O2 diluent at 294 +/- 1 K; decomposition via C-C bond scission giving a formate (RfOCHO), or reaction with 02 giving an acetate (RfOC(O)CH3). In 700 Torr of N2/O2 diluent at 294 +/- 1 K the rate constant ratio k(O2)/k(diss) = (1.26 +/- 0.74) x 10(-19) cm3 molecule(-1). The OH radical initiated atmospheric oxidation of RfOC2H5 gives Rf0CHO and RfOC(O)CH3 as major products. RfOC2H5 has a global warming potential of approximately 55 for a 100 year horizon. The results are discussed with respect to the atmospheric chemistry and environmental impact of RfOC2H5.

Rovibronic bands of the Ã←X̃ transition of CH3OO and CD3OO detected with cavity ringdown absorption near 1.2–1.4μm

We have recorded several rovibronic bands of CH3OO and CD3OO in their A<--X transitions in the range of 1.18-1.40 microm with the cavity ringdown technique. While the electronic origins for these species have been reported previously, many newly observed rovibronic bands are described here. The experimental vibrational frequencies (given as nu in the unit cm(-1) in this paper) for the COO bending (nu8) and COO symmetric stretching (nu7) modes in the A state are 378 and 887 cm(-1) for CH3OO, and 348 and 824 cm(-1) for CD3OO, respectively. In addition, two other vibrational frequencies were observed for the A state of CD3OO, namely, nu5 (954 cm(-1)) and nu6 (971 cm(-1)). These experimental vibrational frequencies for the A state of both CH3OO and CD3OO are in good agreement with predictions from quantum-chemical calculations at the UB3LYP/aug-cc-pVTZ level. The enhanced activity of the nu5 vibrational mode in CD3OO is rationalized by mode mixing with the nu7 mode, as supported by calculations of multidimensional Franck-Condon factors. In addition, many hot bands involving the methyl torsional mode (nu12) are observed for both normal and deuterated methyl peroxy. These bands include the "typical" sequence transitions and some "atypical" ones due to the nature of the eigenvalues and eigenfunctions which are a consequence of the low, but very different, torsional barriers in the X and A states. In addition, the 12(2)2 band in CH3OO and the 12(3)3 band in CD3OO show quite different structures than the origin bands, an effect which results from tunneling splittings comparable to the rotational contour.

Quantum chemical and master equation studies of the methyl vinyl carbonyl oxides formed in isoprene ozonolysis

Methyl vinyl carbonyl oxide is an important intermediate in the reaction of isoprene and ozone and may be responsible for most of the (*)OH formed in isoprene ozonolysis. We use CBS-QB3 calculations and RRKM/master equation simulations to characterize all the pathways leading to the formation of this species, all the interconversions among its four possible conformers, and all of its irreversible isomerizations. Our calculations, like previous studies, predict (*)OH yields consistent with experiment if thermalized syn-methyl carbonyl oxides form (*)OH quantitatively. Natural bond order analysis reveals that the vinyl group weakens the C=O bond of the carbonyl oxide, making rotation about this bond accessible to this chemically activated intermediate. The vinyl group also allows one conformer of the carbonyl oxide to undergo electrocyclization to form a dioxole, a species not previously considered in the literature. Dioxole formation, which has a CBS-QB3 reaction barrier of 13.9 kcal/mol, is predicted to be favored over vinyl hydroperoxide formation, dioxirane formation, and collisional stabilization. Our calculations also predict that two dioxole derivatives, 1,2-epoxy-3-butanone and 3-oxobutanal, should be major products of isoprene ozonolysis.

A systematic computational study on the reactions of HO2 with RO2: The HO2 + CH3O2(CD3O2) and HO2 + CH2FO2 reactions

A systematic theoretical study of the reactions of HO2 with RO2 has been carried out. The major concern of the present work is to gain insight into the reaction mechanism and then to explain experimental observations and to predict new product channels for this class of reactions of importance in the atmosphere. In this paper, the reaction mechanisms for two reactions, namely, HO2 + CH3O2 and HO2 + CH2FO2, are reported. Both singlet and triplet potential energy surfaces are investigated. The complexity of the present system makes it impossible to use a single ab initio method to map out all the reaction paths. Various ab initio methods including MP2, CISD, QCISD(T), CCSD(T), CASSCF, and density function theory (B3LYP) have been employed with the basis sets ranging from 6-31G(d) to an extrapolated complete basis set (CBS) limit. It has been established that the CCSD(T)/cc-pVDZ//B3LYP/6-311G(d,p) scheme represents the most feasible method for our systematic study. For the HO2 + CH3O2 reaction, the production of CH3OOH is determined to be the dominant channel. For the HO2 + CH2FO2 reaction, both CH2FOOH and CHFO are major products, whereas the formation of CHFO is dominant in the overall reaction. The computational findings give a fair explanation for the experimental observation of the products.

Measurements and Modeling of DO2 Formation in the Reactions of C2D5 and C3D7 Radicals with O2

Time-resolved production of HO2 and DO2 from the reactions of nondeuterated and deuterated ethyl and propyl radicals with O2 are measured as a function of temperature and pressure in the "transition region" between 623 and 748 K using the technique of laser photolysis/long path frequency modulation spectroscopy. Experimental measurements, using both pulsed-photolytic Cl-atom-initiated oxidation of ethane and propane and direct photolysis of ethyl, n-propyl, and isopropyl iodides, are compared to kinetic models based on the results of time-dependent master equation calculations with ab initio characterization of stationary points. The formation of DO2 and HO2 from the subsequent reaction of the alkyl radicals with O2 is followed by infrared frequency modulation spectroscopy. The concentration of I atoms is simultaneously monitored by direct absorption of a second laser probe on the spin-orbit transition. The kinetic models accurately describe the time scale and amplitude of the DO2 and HO2 formation resulting from C2D5 + O2, n-C3D7 + O2, i-C3D7 + O2, and i-C3H7 + O2. Overall, a very good level of agreement is found between theory and experiments over a wide range of temperatures, pressures, and O2 concentrations. Good agreement is also found between previous literature studies and the theory presented in this work except in the case of the high-temperature rate coefficients for the reaction of i-C3H7 + O2 to form propene. A reinvestigation of the high-temperature kinetics of the i-C3H7 + O2 reaction appears warranted. The results from the present work suggest that the theory for formation of HO2 from the reactions of ethyl and both isomeric forms of propyl radicals with O2 are very well established at this time. It is hoped that these reactions can now form the groundwork for the study and interpretation of larger and more complex R + O2 systems.

Investigation of the radical product channel of the CH3C(O)O2 + HO2 reaction in the gas phase

The reaction of CH(3)C(O)O(2) with HO(2) has been investigated at 296 K and 700 Torr using long path FTIR spectroscopy, during photolysis of Cl(2)/CH(3)CHO/CH(3)OH/air mixtures. The branching ratio for the reaction channel forming CH(3)C(O)O, OH and O(2) (reaction ) has been determined from experiments in which OH radicals were scavenged by addition of benzene to the system, with subsequent formation of phenol used as the primary diagnostic for OH radical formation. The dependence of the phenol yield on benzene concentration was found to be consistent with its formation from the OH-initiated oxidation of benzene, thereby confirming the presence of OH radicals in the system. The dependence of the phenol yield on the initial peroxy radical precursor reagent concentration ratio, [CH(3)OH](0)/[CH(3)CHO](0), is consistent with OH formation resulting mainly from the reaction of CH(3)C(O)O(2) with HO(2) in the early stages of the experiments, such that the limiting yield of phenol at high benzene concentrations is well-correlated with that of CH(3)C(O)OOH, a well-established product of the CH(3)C(O)O(2) + HO(2) reaction (via channel (3a)). However, a delayed source of phenol was also identified, which is attributed mainly to an analogous OH-forming channel of the reaction of HO(2) with HOCH(2)O(2) (reaction ), formed from the reaction of HO(2) with product HCHO. This was investigated in additional series of experiments in which Cl(2)/CH(3)OH/benzene/air and Cl(2)/HCHO/benzene/air mixtures were photolysed. The various reaction systems were fully characterised by simulations using a detailed chemical mechanism. This allowed the following branching ratios to be determined: CH(3)C(O)O(2) + HO(2)--> CH(3)C(O)OOH + O(2), k(3a)/k(3) = 0.38 +/- 0.13; --> CH(3)C(O)OH + O(3), k(3b)/k(3) = 0.12 +/- 0.04; --> CH(3)C(O)O + OH + O(2), k(3c)/k(3) = 0.43 +/- 0.10: HOCH(2)O(2) + HO(2)--> HCOOH + H(2)O + O(2), k(17b)/k(17) = 0.30 +/- 0.06; --> HOCH(2)O + OH + O(2), k(17c)/k(17) = 0.20 +/- 0.05. The results therefore provide strong evidence for significant participation of the radical-forming channels of these reactions, with the branching ratio for the title reaction being in good agreement with the value reported in one previous study. As part of this work, the kinetics of the reaction of Cl atoms with phenol (reaction (14)) have also been investigated. The rate coefficient was determined relative to the rate coefficient for the reaction of Cl with CH(3)OH, during the photolysis of mixtures of Cl(2), phenol and CH(3)OH, in either N(2) or air at 296 K and 760 Torr. A value of k(14) = (1.92 +/- 0.17) x 10(-10) cm(3) molecule(-1) s(-1) was determined from the experiments in N(2), in agreement with the literature. In air, the apparent rate coefficient was about a factor of two lower, which is interpreted in terms of regeneration of phenol from the product phenoxy radical, C(6)H(5)O, possibly via its reaction with HO(2).

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.

Kinetics and branching ratio studies of the reaction of C2H5O2 + HO2 using chemical ionisation mass spectrometry

The overall rate coefficient for the reaction of C(2)H(5)O(2) with HO(2) was determined using a turbulent flow chemical ionization mass spectrometer (TF-CIMS) system over the pressure range of 75 to 200 Torr and temperatures between 195 and 298 K. The temperature dependence of the overall rate coefficient for the reaction between C(2)H(5)O(2) and HO(2) was fitted using the following Arrhenius expression: k(T) = (2.08) x 10(-13) exp [(864 +/- 79)/T] cm(-3) molecule(-1) s(-1). The upper limits for the branching ratios for reactive channels leading to O(3) and OH production were quantified for the first time. A tropospheric model has been used to assess the impact of the experimental error of the rate coefficients determined in this study on predicted concentrations of a number of key species, including O(3), OH, HO(2), NO and NO(2). In all cases it is found that the propagated error is very small and will not in itself be a major cause of uncertainty in modelled concentrations. However, at low temperatures, where there is a wide discrepancy between existing kinetic studies, modelling using the range of kinetic data in the literature shows a small but significant variation for [C(2)H(5)O(2)], [C(2)H(5)OOH], [NO(x)] and the HO(2) : OH ratio. Furthermore, a structure-activity relationship (SAR) was developed to rationalise the reactivity of the reaction between RO(2) and HO(2).

Photodissociation spectroscopy of the Mg+-acetic acid complex

We have studied the structure and photodissociation of Mg(+)-acetic acid clusters. Ab initio calculations suggest four relatively strongly bound ground state isomers for the [MgC(2)H(4)O(2)](+) complex. These isomers include the cis and trans forms of the Mg(+)-acetic acid association complex with Mg(+) bonded to the carbonyl O atom of acetic acid, the Mg(+)-acetic acid association complex with Mg(+) bonded to the hydroxyl O atom of acetic acid, or to a Mg(+)-ethenediol association complex. Photodissociation through the Mg(+)-based 3p<--3s absorption bands in the near UV leads to direct (nonreactive) and reactive dissociation products: Mg(+), MgOH(+), Mg(H(2)O)(+), CH(3)CO(+), and MgCH(3) (+). At low energies the dominant reactive quenching pathway is through dehydration to Mg(H(2)O)(+), but additional reaction channels involving C-H and C-C bond activation are also open at higher energies.