Formation of a Criegee intermediate in the low-temperature oxidation of dimethyl sulfoxide

Dimethyl sulfoxide (DMSO) is the major sulfur-containing constituent of the Marine Boundary Layer. It is a significant source of H2SO4 aerosol/particles and methane sulfonic acid via atmospheric oxidation processes, where the mechanism is not established. In this study, several new, low-temperature pathways are revealed in the oxidation of DMSO using CBS-QB3 and G3MP2 multilevel and B3LYP hybrid density functional quantum chemical methods. Unlike analogous hydrocarbon peroxy radicals the chemically activated DMSO peroxy radical, [CH3S(=O)CH2OO*]*, predominantly undergoes simple dissociation to a methylsulfinyl radical CH3S*(=O) and a Criegee intermediate, CH2OO, with the barrier to dissociation 11.3 kcal mol(-1) below the energy of the CH3S(=O)CH2* + O2 reactants. The well depth for addition of O2 to the CH3S(=O)CH2 precursor radical is 29.6 kcal mol(-1) at the CBS-QB3 level of theory. We believe that this reaction may serve an important role in atmospheric photochemical and irradiated biological (oxygen-rich) media where formation of initial radicals is facilitated even at lower temperatures. The Criegee intermediate (carbonyl oxide, peroxymethylene) and sulfinyl radical can further decompose, resulting in additional chain branching. A second reaction channel important for oxidation processes includes formation (via intramolecular H atom transfer) and further decomposition of hydroperoxide methylsulfoxide radical, *CH2S(=O)CH2OOH over a low barrier of activation. The initial H-transfer reaction is similar and common in analogous hydrocarbon radical + O2 reactions; but the subsequent very low (3-6 kcal mol(-1)) barrier (14 kcal mol(-1) below the initial reagents) to beta-scission products is not common in HC systems. The low energy reaction of the hydroperoxide radical is a beta-scission elimination of *CH2S(=O)CH2OOH into the CH2=S=O + CH2O + *OH product set. This beta-scission barrier is low, because of the delocalization of the *CH2 radical center through the -S(=O) group, to the -CH2OOH fragment in the transition state structure. The hydroperoxide methylsulfoxide radical can also decompose via a second reaction channel of intramolecular OH migration, yielding formaldehyde and a sulfur-centered hydroxymethylsulfinyl radical HOCH2S*(=O). The barrier of activation relative to initial reagents is 4.2 kcal mol(-1). Heats of formation for DMSO, DMSO carbon-centered radical and Criegee intermediate are evaluated at 298 K as -35.97 +/- 0.05, 13.0 +/- 0.2 and 25.3 +/- 0.7 kcal mol(-1) respectively using isodesmic reaction analysis. The [CH3S*(=O) + CH2OO] product set is shown to form a van der Waals complex that results in O-atom transfer reaction and the formation of new products CH3SO2* radical and CH2O. Proper orientation of the Criegee intermediate and methylsulfinyl radical, as a pre-stabilized pre-reaction complex, assist the process. The DMSO radical reaction is also compared to that of acetonyl radical.     

Thermochemical properties, DeltafH degrees (298), S degrees (298), and Cp degrees (T), for n-butyl and n-pentyl hydroperoxides and the alkyl and peroxy radicals, transition states, and kinetics for intramolecular hydrogen shift reactions of the peroxy radicals

Alkyl radicals in atmospheric and combustion environments undergo a rapid association with molecular oxygen (3O2) to form an alkyl peroxy radical (ROO*). One important reaction of these peroxy radicals is the intramolecular H-shift (intramolecular abstraction) to form a hydroperoxide alkyl radical (R'*COOH), where the hydroperoxide alkyl radical may undergo chemical activation reaction with O2 and result in chain branching at moderate to low temperatures. The thermochemistry and trends in kinetic parameters for the hydrogen shift reactions from each carbon (4-8-member-ring TST's) in n-butyl and n-pentyl peroxy radicals (CCCCOO* and CCCCCOO*) are analyzed using density functional and ab initio calculation methods. Thermochemical properties, DeltafH degrees (298 K), C-H bond energies, S degrees (298 K), and Cp degrees (T) of saturated linear C4 and C5 aliphatic peroxides (ROOH), as well as the corresponding hydroperoxide alkyl radicals (R'*COOH), are determined. DeltafH degrees (298 K) are obtained from isodesmic reactions and the total energies of the CBS-QB3 and B3LYP computational methods. Contributions to the entropy and the heat capacity from translation, vibration, and external rotation are calculated using the rigid-rotor-harmonic-oscillator approximation based on the CBS-QB3 frequencies and structures. The results indicate that pre-exponential factors, A(T), decrease with the increase of the ring size (4-8-member-ring TS, H-atom included). The DeltaH for 4-, 5-, 6-, and 7-member rings in n-butyl (and n-pentyl) peroxy are 40.8 (40.8), 31.4 (31.5), 20.5 (20.0), 22.6-p (19.4) kcal mol(-1), respectively. The DeltaH for the 8-member ring in n-pentylperoxy is 23.8-p kcal mol(-1), All abstractions are from secondary (-CH2-) groups except those marked (-p), which are from primary sites. Enthalpy and barrier values from the B3LYP/6-311++G(2d,p) and BHandHLYP/6-311G(d,p) methods are compared with CBS-QB3 results. The B3LYP results show good agreement with the higher level CBS-QB3 calculation method; the BHandH barriers for the intramolecular peroxy H-shifts are not acceptable.     

Unusually fast 1,6-h shifts of enolic hydrogens in peroxy radicals: formation of the first-generation C2 and C3 carbonyls in the oxidation of isoprene

In a theoretical investigation using the CBS-QB3//UB3LYP/6-31+G** method supported by higher-level computations such as CBS-QB3//UQCISD/6-31+G**, the 1,6-H shifts of the enolic hydrogen in peroxy radicals of the type Z-HO-CH═CH-CH(2)-OO(?) were found to face exceptionally low energy barriers of only about 11 kcal mol(-1)--i.e., 6-9 kcal mol(-1) lower than the barriers for similar shifts of alkane hydrogens--such that they can proceed at unequaled rates of order 10(5) to 10(6) s(-1) at ambient temperatures. The unusually low barriers for enolic 1,6-H shifts in peroxy radicals, characterized here for the first time to our knowledge, are rationalized. As cases in point, the secondary peroxy radicals Z-HO-CH═C(CH(3))-CH(OO(?))-CH(2)OH (case A) and Z-HO-CH═CH-C(CH(3))(OO(?))-CH(2)OH (case B) derived from the primary Z-δ-hydroxy-peroxy radicals in the oxidation of isoprene, are predicted to undergo 1,6-H shifts of their enolic hydrogens at TST-calculated rates in the range 270-320 K of k(T)(A) = 5.4 × 10(-4) × T(5.04) × exp(-1990/T) s(-1) and k(T)(B) = 109 × T(3.13) × exp(-3420/T) s(-1), respectively, i.e., 2.0 × 10(6) and 6.2 × 10(4) s(-1), respectively, at 298 K, far outrunning in all relevant atmospheric and laboratory conditions their reactions with NO proposed earlier as their dominant pathways (Dibble J. Phys. Chem. A 2004, 108, 2199). These fast enolic-H shifts are shown to provide the explanation for the first-generation formation of methylglyoxal + glycolaldehyde, and glyoxal + hydroxyacetone in the oxidation of isoprene under high-NO conditions, recently determined by several groups. However, under moderate- and low-NO atmospheric conditions, the fast interconversion and equilibration of the various thermally labile, initial peroxy conformers/isomers from isoprene and the isomerization of the initial Z-δ-hydroxy-peroxy radicals, both recently proposed by us (Peeters et al. Phys. Chem. Chem. Phys. 2009, 11, 5935), are expected to substantially reduce the yields of the small carbonyls at issue.     

Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.     

12]Annulene radical anions revisited: evaluation of structure assignments based on computed energetic and electron spin resonance data

We report density functional and coupled cluster calculations on numerous monocyclic and bicyclic (CH)12(*-) isomers. At the RCCSD(T)/cc-pVDZ//UB3LYP/6-31+G* level, a nearly planar, bond-equalized radical anion of 1,7-di-trans-[12]annulene (4a(*-)) is lowest in energy; several other isomers and conformations lie within 3 kcal/mol of 4a(*-). RCCSD(T)/AUG-cc-pVDZ//UB3LYP/6-31+G* results place the all-cis isomer 3(*-) slightly below 4a(*-) in energy. Validation studies on the heptalene radical anion, [16]annulene radical anion, and tri-trans-[12]annulene radical anion indicate that electron spin resonance (ESR) hyperfine coupling constants (aH values) computed at the BLYP/EPR-III level on DFT geometries give much better agreement with experimental values than those computed using B3LYP/6-31G*. We were unable to locate any C12H12(*-) isomer that could account for the ESR spectrum previously attributed to a highly twisted structure for the 1,7-di-trans-[12]annulene radical anion. Our computed energetic and ESR data for [12]annulene radical anions and their valence isomers suggest that 4a(*-) may have been made, yet its ESR spectrum was incorrectly assigned to the bicyclic isomer 6b(*-). Finally, the computed (1)H NMR shift values of the dianion of 4 reveal a distinct diatropic ring current that should aid in its characterization.     

Reactions of hydrated electron with various radicals: spin factor in diffusion-controlled reactions

The reactions of hydrated electron (eaq-) with various radicals have been studied in pulse radiolysis experiments. These radicals are hydroxyl radical (*OH), sulfite radical anion (*SO3-), carbonate radical anion (CO3*-), carbon dioxide radical anion (*CO2-), azidyl radical (*N3), dibromine radical anion (Br2*-), diiodine radical anion (I2*-), 2-hydroxy-2-propyl radical (*C(CH3)2OH), 2-hydroxy-2-methyl-1-propyl radical ((*CH2)(CH3)2COH), hydroxycyclohexadienyl radical (*C6H6OH), phenoxyl radical (C6H5O*), p-methylphenoxyl radical (p-(H3C)C6H4O*), p-benzosemiquinone radical anion (p-OC6H4O*-), and phenylthiyl radical (C6H5S*). The kinetics of eaq- was followed in the presence of the counter radicals in transient optical absorption measurements. The rate constants of the eaq- reactions with radicals have been determined over a temperature range of 5-75 degrees C from the kinetic analysis of systems of multiple second-order reactions. The observed high rate constants for all the eaq- + radical reactions have been analyzed with the Smoluchowski equation. This analysis suggests that many of the eaq- + radical reactions are diffusion-controlled with a spin factor of 1/4, while other reactions with *OH, *N3, Br2*-, I2*-, and C6H5S* have spin factors significantly larger than 1/4. Spin dynamics for the eaq-/radical pairs is discussed to explain the different spin factors. The reactions with *OH, *N3, Br2*-, and I2*- have also been found to have apparent activation energies less than that for diffusion control, and it is suggested that the spin factors for these reactions decrease with increasing temperature. Such a decrease in spin factor may reflect a changing competition between spin relaxation/conversion and diffusive escape from the radical pairs.     

A temperature-dependent relative-rate study of the OH initiated oxidation of n-butane: the kinetics of the reactions of the 1- and 2-butoxy radicals

The kinetics of the reactions of 1-and 2-butoxy radicals have been studied using a slow-flow photochemical reactor with GC-FID detection of reactants and products. Branching ratios between decomposition, CH3CH(O*)CH2CH3 --> CH3CHO + C2H5, reaction (7), and reaction with oxygen, CH3CH(O*)CH2CH3+ O2 --> CH3C(O)C2H5+ HO2, reaction (6), for the 2-butoxy radical and between isomerization, CH3CH2CH2CH2O* --> CH2CH2CH2CH2OH, reaction (9), and reaction with oxygen, CH3CH2CH2CH2O* + O2 --> C3H7CHO + HO2, reaction (8), for the 1-butoxy radical were measured as a function of oxygen concentration at atmospheric pressure over the temperature range 250-318 K. Evidence for the formation of a small fraction of chemically activated alkoxy radicals generated from the photolysis of alkyl nitrite precursors and from the exothermic reaction of 2-butyl peroxy radicals with NO was observed. The temperature dependence of the rate constant ratios for a thermalized system is given by k7/k6= 5.4 x 10(26) exp[(-47.4 +/- 2.8 kJ mol(-1))/RT] molecule cm(-3) and k9/k8= 1.98 x 10(23) exp[(-22.6 +/- 3.9 kJ mol(-1))/RT] molecule cm(-3). The results agree well with the available experimental literature data at ambient temperature but the temperature dependence of the rate constant ratios is weaker than in current recommendations.     

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.     

Structures, internal rotor potentials, and thermochemical properties for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals

Structures, enthalpy (Δ(f)H°(298)), entropy (S°(T)), and heat capacity (C(p)(T)) are determined for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals using the density functional B3LYP and composite CBS-QB3 calculations. Enthalpies of formation (Δ(f)H°(298)) are determined at the B3LYP/6-31G(d,p), B3LYP/6-31+G(2d,2p), and composite CBS-QB3 levels using several work reactions for each species. Entropy (S) and heat capacity (C(p)(T)) values from vibration, translational, and external rotational contributions are calculated using the rigid-rotor-harmonic-oscillator approximation based on the vibration frequencies and structures obtained from the density functional studies. Contribution to Δ(f)H(T), S, and C(p)(T) from the analysis on the internal rotors is included. Recommended values for enthalpies of formation of the most stable conformers of nitroacetone cc(═o)cno2, acetonitrite cc(═o)ono, nitroacetate cc(═o)no2, and acetyl nitrite cc(═o)ono are -51.6 kcal mol(-1), -51.3 kcal mol(-1), -45.4 kcal mol(-1), and -58.2 kcal mol(-1), respectively. The calculated Δ(f)H°(298) for nitroethylene c═cno2 is 7.6 kcal mol(-1) and for vinyl nitrite c═cono is 7.2 kcal mol(-1). We also found an unusual phenomena: an intramolecular transfer reaction (isomerization) with a low barrier (3.6 kcal mol(-1)) in the acetyl nitrite. The NO of the nitrite (R-ONO) in CH(3)C(═O')ONO moves to the C═O' oxygen in a motion of a stretching frequency and then a shift to the carbonyl oxygen (marked as O' for illustration purposes).     

Preparation and reactivity of [D3d]-octahedrane: the most stable (CH)12 hydrocarbon

The synthesis of the (CH)12 hydrocarbon [D(3d)]-octahedrane (heptacyclo[6.4.0.0(2,4).0(3,7).0(5,12).0(6,10).0(9,11)]dodecane) 1 and its selective functionalization retaining the hydrocarbon cage is described. The B3LYP/6-311+G* strain energy of 1 is 83.7 kcal mol(-1) (4.7 kcal mol(-1) per C-C bond) which is significantly higher than that of the structurally related (CH)16 [D(4d)]-decahedrane 2 (75.4 kcal mol(-1); 3.1 kcal mol(-1) per C-C bond) and (CH)20 [I(h)]-dodecahedrane 3 (51.5 kcal mol(-1); 1.7 kcal mol(-1) per C-C bond); the heats of formation for 1-3 computed according to homodesmotic equations are 52, 35, and 4 kcal mol(-1). Catalytic hydrogenation of 1 leads to consecutive opening of the two cyclopropane rings to give C2-bisseco-octahedrane (pentacyclo[6.4.0.0(2,6).0(3,11).0(4,9)]dodecane) 16 as the major product. Although 1 is highly strained, its carbon skeleton is kinetically quite stable: Upon heating, 1 does not decompose until above 180 degrees C. The B3LYP/6-31G* barriers for the S(R)2 attack of the tBuO. and Br3C. radicals on a carbon atom of one of the cyclopropane fragments (Delta(298) = 27-28 kcal mol(-1)) are higher than those for hydrogen atom abstraction. The latter barriers are virtually identical for the abstraction from the C1-H and C2-H positions with the tBuO. radical (DeltaG(298) = 17.4 and 17.9 kcal mol(-1), respectively), but significantly different for the reaction at these positions with the Br3C. radical (DeltaG(298) = 18.8 and 21.0 kcal mol(-1)). These computational results agree well with experiments, in which the chlorination of 1 with tert-butyl hypochlorite gave a mixture of 1- and 2-chlorooctahedranes (ratio 3:2). The bromination with carbon tetrabromide under phase-transfer catalytic (PTC) conditions (nBu4NBr/NaOH) selectively gave 1-bromooctahedrane in 43 % isolated yield. For comparison, the PTC bromination was also applied to 2,4-dehydroadamantane yielding 54 % 7-bromo-2,4-dehydroadamantane.     

Ab Initio Kinetics for Thermal Decomposition of CH(3)N(•)NH(2), cis-CH(3)NHN(•)H, trans-CH(3)NHN(•)H, and C(•)H(2)NNH(2) Radicals

The thermal decomposition of the CH(3)N(?)NH(2), cis-CH(3)NHN(?)H, trans-CH(3)NHN(?)H, and C(?)H(2)NNH(2) radicals, which are the four radical products from the H-abstraction reactions of monomethylhydrazine, were theoretically studied by using ab initio Rice-Ramsperger-Kassel-Marcus (RRKM) transition-state theory and master equation analysis. Various decomposition pathways were identified by using either the QCISD(T)/cc-pV∞Z//CASPT2/aug-cc-pVTZ or the QCISD(T)/cc-pV∞Z//B3LYP/6-311++G(d,p) quantum chemistry methods. The results reveal that the β-scission of NH(2) to form methyleneimine is the predominant channel for the decomposition of the C(?)H(2)NNH(2) radical due to its small energy barrier of 13.8 kcal mol(-1). The high pressure limit rate coefficient for the reaction is fitted by 3.88 × 10(19)T(-1.672) exp(-9665.13/T) s(-1). In addition, the pressure dependent rate coefficients exhibit slight temperature dependence at temperatures of 1000-2500 K. The cis-CH(3)NHN(?)H and trans-CH(3)NHN(?)H radicals are the two distinct spatial isomers with an energy barrier of 26 kcal mol(-1) for their isomerization. The β-scission of CH(3) from the cis-CH(3)NHN(?)H radical to form trans-diazene has an energy barrier of 35.2 kcal mol(-1), and the β-scission of CH(3) from the trans-CH(3)NHN(?)H radical to form cis-diazene has an energy barrier of 39.8 kcal mol(-1). The CH(3)N(?)NH(2) radical undergoes the β-scission of methyl hydrogen and amine hydrogen to form CH(2)═NNH(2), trans-CH(3)N═NH, and cis-CH(3)N═NH products, with the energy barriers of 42.8, 46.0, and 50.2 kcal mol(-1), respectively. The dissociation and isomerization rate coefficients for the reactions were calculated via the E/J resolved RRKM theory and multiple-well master equation analysis at temperatures of 300-2500 K and pressures of 0.01-100 atm. The calculated rate coefficients associated with updated thermochemical property data are essential components in the development of kinetic mechanisms for the pyrolysis and oxidation of MMH and its derivatives.     

Rate and equilibrium constant of the reaction of 1-methylvinoxy radicals with O2: CH3COCH2 + O2<--> CH3COCH2O2

The reaction of 1-methylvinoxy radicals, CH3COCH2, with molecular oxygen has been investigated by experimental and theoretical methods as a function of temperature (291-520 K) and pressure (0.042-10 bar He). Experiments have been performed by laser photolysis coupled to a detection of 1-methylvinoxy radicals by laser-induced fluorescence LIF. The potential energy surface calculations were performed using ab inito molecular orbital theory at the G3MP2B3 and CBSQB3 level of theory based on the density function theory optimized geometries. Derived molecular properties of the characteristic points of the potential energy surface were used to describe the mechanism and kinetics of the reaction under investigation. At 295 K, no pressure dependence of the rate constant for the association reaction has been observed: k(1,298K) = (1.18 +/- 0.04) x 10(-12) cm3 s(-1). Biexponential decays have been observed in the temperature range 459-520 K and have been interpreted as an equilibrium reaction. The temperature-dependent equilibrium constants have been extracted from these decays and a standard reaction enthalpy of deltaH(r,298K) = -105.0 +/- 2.0 kJ mol(-1) and entropy of deltaS(r,298K) = -143.0 +/- 4.0 J mol(-1) K(-1) were derived, in excellent agreement with the theoretical results. Consistent heats of formation for the vinoxy and the 1-methylvinoxy radical as well as their O2 adducts are recommended based on our complementary experimental and theoretical study deltaH(f,298K) = 13.0 +/- 2.0, -32. 9+/- 2.0, -85.9 +/- 4.0, and -142.1 +/- 4.0 kJ mol(-1) for CH2CHO, CH3COCH2 radicals, and their adducts, respectively.     

Cis-trans isomerization of chemically activated 1-methylallyl radical and fate of the resulting 2-buten-1-peroxy radical

The cis-trans isomerization of chemically activated 1-methylallyl is investigated using RRKM/Master Equation methods for a range of pressures and temperatures. This system is a prototype for a large range of allylic radicals formed from highly exothermic (～35 kcal/mol) OH + alkene reactions. Energies, vibrational frequencies, anharmonic constants, and the torsional potential of the methyl group are computed with density functional theory for both isomers and the transition state connecting them. Chemically activated radicals are found to undergo rapid cis-trans isomerization leading to stabilization of significant amounts of both isomers. In addition, the thermal rate constant for trans → cis isomerization of 1-methylallyl is computed to be high enough to dominate reaction with O(2) in 10 atm of air at 700 K, so models of the chemistry of the (more abundant and more commonly studied) trans-alkenes may need to be modified to include the cis isomers of the corresponding allylic radicals. Addition of molecular oxygen to 1-methylallyl radical can form 2-butene-1-peroxy radical (CH(3)CH═CHCH(2)OO(?)), and quantum chemistry is used to thoroughly explore the possible unimolecular reactions of the cis and trans isomers of this radical. The cis isomer of the 2-butene-1-peroxy radical has the lowest barrier (via 1,6 H-shift) to further reaction, but this barrier appears to be too high to compete with loss of O(2).     

Exchange of organic radicals with organo-cobalt complexes formed in the living radical polymerization of vinyl acetate

Exchange of organic radicals between solution and organo-cobalt complexes is experimentally observed and the reaction pathway is probed through DFT calculations. Cyanoisopropyl radicals from AIBN (2,2'-azobisisobutyronitrile) enter solutions of cobalt(II) tetramesityl porphyrin ((TMP)Co(II)*, 1) and vinyl acetate (VAc) in benzene and react to produce transient hydride (TMP)Co-H and radicals (*CH(OAc)CH2C(CH3)2CN (R1*)) that proceed on to form organo-cobalt complexes (TMP)Co-CH(OAc)CH3 (4, Co-R2) and (TMP)Co-CH(OAc)CH2C(CH3)2CN (3, Co-R1), respectively. Rate constants for cyanoisopropyl radical addition with vinyl acetate and hydrogen atom transfer to (TMP)Co(II)* are reported through kinetic studies for the formation and transformation of organo-cobalt species in this system. Rate constants for near-degenerate exchanges of radicals in solution with organo-cobalt complexes are deduced from (1)H NMR studies and kinetic modeling. DFT computations revealed formation of an unsymmetrical adduct of (TMP)Co-CH(OAc)CH3 (4) with *CH(OAc)CH3 (R2*) and support an associative pathway for radical interchange through a three-centered three-electron transition state [R...Co...R]. Associative radical interchange of the latent radical groups in organo-cobalt porphyrin complexes with freely diffusing radicals in solution that is observed in this system provides a pathway for mediation of living radical polymerization of vinyl acetate.     

The cysteine radical cation: structures and fragmentation pathways

A theoretical study on the structures, relative energies, isomerization reactions and fragmentation pathways of the cysteine radical cation, [NH(2)CH(CH(2)SH)COOH].+, is reported. Hybrid density functional theory (B3LYP) has been used in conjunction with the 6-311++G(d,p) basis set. The isomer at the global minimum, Captodative-1, has the structure NH(2)C.(CH(2)SH)C(OH)(2)+; the stability of this ion is attributed to the captodative effect in which the NH(2) functions as a powerful pi-electron donor and C(OH)(2)+ as a powerful pi-electron acceptor. Ion Distonic-S-1, H(3)N(+)CH(CH(2)S.)COOH, in which the radical is formally situated on the S atom, is higher in enthalpy (DeltaH degrees (0)) than Captodative-1 by 6.1 kcal mol(-1), but is lower in enthalpy than another isomer Distonic-C-1, H(3)N(+)C.(CH(2)SH)COOH, by 8.2 kcal mol(-1). Isomerization of the canonical radical cation of cysteine, [H(2)NCH(CH(2)SH)COOH].+, (Canonical-1), to Captodative-1 has an enthalpy of activation of 25.8 kcal mol(-1), while the barrier against isomerization of Canonical-1 to Distonic-S-1 is only 9.6 kcal mol(-1). Two additional transient tautomers, one with the radical located at C(alpha) and the charge on SH(2), and the other a carboxy radical with the charge on NH(3), are reported. Plausible fragmentation pathways (losses of small molecules, CO(2), CH(2)S, H(2)S and NH(3), and neutral radicals COOH. , HSCH(2). and NH(2).) from Canonical-1 are examined.     

Unimolecular reactions in the CF3CH2Cl ↔ CF2ClCH2F system: isomerization by interchange of Cl and F atoms

The recombination of CF(2)Cl and CH(2)F radicals was used to prepare CF(2)ClCH(2)F* molecules with 93 ± 2 kcal mol(-1) of vibrational energy in a room temperature bath gas. The observed unimolecular reactions in order of relative importance were: (1) 1,2-ClH elimination to give CF(2)═CHF, (2) isomerization to CF(3)CH(2)Cl by the interchange of F and Cl atoms and (3) 1,2-FH elimination to give E- and Z-CFCl═CHF. Since the isomerization reaction is 12 kcal mol(-1) exothermic, the CF(3)CH(2)Cl* molecules have 105 kcal mol(-1) of internal energy and they can eliminate HF to give CF(2)═CHCl, decompose by rupture of the C-Cl bond, or isomerize back to CF(2)ClCH(2)F. These data, which provide experimental rate constants, are combined with previously published results for chemically activated CF(3)CH(2)Cl* formed by the recombination of CF(3) and CH(2)Cl radicals to provide a comprehensive view of the CF(3)CH(2)Cl* ? CF(2)ClCH(2)F* unimolecular reaction system. The experimental rate constants are matched to calculated statistical rate constants to assign threshold energies for the observed reactions. The models for the molecules and transition states needed for the rate constant calculations were obtained from electronic structures calculated from density functional theory. The previously proposed explanation for the formation of CF(2)═CHF in thermal and infrared multiphoton excitation studies of CF(3)CH(2)Cl, which was 2,2-HCl elimination from CF(3)CH(2)Cl followed by migration of the F atom in CF(3)CH, should be replaced by the Cl/F interchange reaction followed by a conventional 1,2-ClH elimination from CF(2)ClCH(2)F. The unimolecular reactions are augmented by free-radical chemistry initiated by reactions of Cl and F atoms in the thermal decomposition of CF(3)CH(2)Cl and CF(2)ClCH(2)F.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

Branching ratios for the reaction of selected carbonyl-containing peroxy radicals with hydroperoxy radicals

An important chemical sink for organic peroxy radicals (RO(2)) in the troposphere is reaction with hydroperoxy radicals (HO(2)). Although this reaction is typically assumed to form hydroperoxides as the major products (R1a), acetyl peroxy radicals and acetonyl peroxy radicals have been shown to undergo other reactions (R1b) and (R1c) with substantial branching ratios: RO(2) + HO(2) → ROOH + O(2) (R1a), RO(2) + HO(2) → ROH + O(3) (R1b), RO(2) + HO(2) → RO + OH + O(2) (R1c). Theoretical work suggests that reactions (R1b) and (R1c) may be a general feature of acyl peroxy and α-carbonyl peroxy radicals. In this work, branching ratios for R1a-R1c were derived for six carbonyl-containing peroxy radicals: C(2)H(5)C(O)O(2), C(3)H(7)C(O)O(2), CH(3)C(O)CH(2)O(2), CH(3)C(O)CH(O(2))CH(3), CH(2)ClCH(O(2))C(O)CH(3), and CH(2)ClC(CH(3))(O(2))CHO. Branching ratios for reactions of Cl-atoms with butanal, butanone, methacrolein, and methyl vinyl ketone were also measured as a part of this work. Product yields were determined using a combination of long path Fourier transform infrared spectroscopy, high performance liquid chromatography with fluorescence detection, gas chromatography with flame ionization detection, and gas chromatography-mass spectrometry. The following branching ratios were determined: C(2)H(5)C(O)O(2), Y(R1a) = 0.35 ± 0.1, Y(R1b) = 0.25 ± 0.1, and Y(R1c) = 0.4 ± 0.1; C(3)H(7)C(O)O(2), Y(R1a) = 0.24 ± 0.15, Y(R1b) = 0.29 ± 0.1, and Y(R1c) = 0.47 ± 0.15; CH(3)C(O)CH(2)O(2), Y(R1a) = 0.75 ± 0.13, Y(R1b) = 0, and Y(R1c) = 0.25 ± 0.13; CH(3)C(O)CH(O(2))CH(3), Y(R1a) = 0.42 ± 0.1, Y(R1b) = 0, and Y(R1c) = 0.58 ± 0.1; CH(2)ClC(CH(3))(O(2))CHO, Y(R1a) = 0.2 ± 0.2, Y(R1b) = 0, and Y(R1c) = 0.8 ± 0.2; and CH(2)ClCH(O(2))C(O)CH(3), Y(R1a) = 0.2 ± 0.1, Y(R1b) = 0, and Y(R1c) = 0.8 ± 0.2. The results give insights into possible mechanisms for cycling of OH radicals in the atmosphere.     

Dissociation of the ground state vinoxy radical and its photolytic precursor chloroacetaldehyde: electronic nonadiabaticity and the suppression of the H+ketene channel

This work is a study of the competition between the two unimolecular reaction channels available to the vinoxy radical (CH(2)CHO), C-H fission to form H+ketene, and isomerization to the acetyl radical (CH(3)CO) followed by C-C fission to form CH(3) + CO. Chloroacetaldehyde (CH(2)ClCHO) was used as a photolytic precursor to the vinoxy radical in its ground state; photodissociation of chloroacetaldehyde at 193 nm produces vinoxy radicals with internal energies spanning the G3//B3LYP calculated barriers to the two available unimolecular reaction channels. The onset of the CH(3) + CO channel, via isomerization to the acetyl radical, was found to occur at an internal energy of 41 +/- 2 kcal/mol, agreeing well with our calculated isomerization barrier of 40.8 kcal/mol. Branching to the H+ketene channel was too small to be detected; we conclude that the branching to the H+ketene channel must be at least a factor of 200 lower than what is predicted by a RRKM analysis based on our electronic structure calculations. This dramatic result may be explained in part by the presence of a conical intersection at planar geometries along the reaction coordinate leading to H+ketene, which results in electronically nonadiabatic recrossing of the transition state.     

Thermochemistry and kinetics of the 2-butanone-4-yl CH3C(=O)CH2CH2• + O2 reaction system

Thermochemistry and kinetic pathways on the 2-butanone-4-yl (CH₃C(=O)CH₂CH₂•) + O₂ reaction system are determined. Standard enthalpies, entropies, and heat capacities are evaluated using the G3MP2B3, G3, G3MP3, CBS-QB3 ab initio methods, and the B3LYP/6-311g(d,p) density functional calculation method. The CH₃C(=O)CH₂CH₂• radical + O₂ association reaction forms a chemically activated peroxy radical with 35 kcal mol⁻¹ excess of energy. The chemically activated adduct can undergo RO−O bond dissociation, rearrangement via intramolecular hydrogen transfer reactions to form hydroperoxide-alkyl radicals, or eliminate HO₂ and OH. The hydroperoxide-alkyl radical intermediates can undergo further reactions forming ketones, cyclic ethers, OH radicals, ketene, formaldehyde, or oxiranes. A relatively new path showing a low barrier and resulting in reactive product sets involves peroxy radical attack on a carbonyl carbon atom in a cyclic transition state structure. It is shown to be important in ketones when the cyclic transition state has five or more central atoms.