Hydrogen transfer between sulfuric acid and hydroxyl radical in the gas phase: competition among hydrogen atom transfer, proton-coupled electron-transfer, and double proton transfer

In an attempt to assess the potential role of the hydroxyl radical in the atmospheric degradation of sulfuric acid, the hydrogen transfer between H2SO4 and HO* in the gas phase has been investigated by means of DFT and quantum-mechanical electronic-structure calculations, as well as classical transition state theory computations. The first step of the H2SO4 + HO* reaction is the barrierless formation of a prereactive hydrogen-bonded complex (Cr1) lying 8.1 kcal mol(-1) below the sum of the (298 K) enthalpies of the reactants. After forming Cr1, a single hydrogen transfer from H2SO4 to HO* and a degenerate double hydrogen-exchange between H2SO4 and HO* may occur. The single hydrogen transfer, yielding HSO4* and H2O, can take place through three different transition structures, the two lowest energy ones (TS1 and TS2) corresponding to a proton-coupled electron-transfer mechanism, whereas the higher energy one (TS3) is associated with a hydrogen atom transfer mechanism. The double hydrogen-exchange, affording products identical to reactants, takes place through a transition structure (TS4) involving a double proton-transfer mechanism and is predicted to be the dominant pathway. A rate constant of 1.50 x 10(-14) cm(3) molecule(-1) s(-1) at 298 K is obtained for the overall reaction H2SO4 + HO*. The single hydrogen transfer through TS1, TS2, and TS3 contributes to the overall rate constant at 298 K with a 43.4%. It is concluded that the single hydrogen transfer from H2SO4 to HO* yielding HSO4* and H2O might well be a significant sink for gaseous sulfuric acid in the atmosphere.     

Mechanism for the gas-phase reaction between formaldehyde and hydroperoxyl radical. A theoretical study

We present a high-level theoretical study on the gas-phase reaction between formaldehyde and hydroperoxyl radical carried out using the DFT-B3LYP, QCISD, and CCSD(T) theoretical approaches in connection with the 6-311+G(d,p), 6-311+G(2df,2p), and aug-cc-pVTZ basis sets. The most favorable reaction path begins with the formation of a pre-reactive complex and produces the peroxy radical CH(2)(OO)OH in a process that is computed to be exothermic by 16.8 kcal/mol. This reaction involves a process in which the oxygen terminal of the HO(2) moiety adds to the carbon of formaldehyde, and, simultaneously, the hydrogen of the hydroperoxyl group is transferred to the oxygen of the carbonyl in a proton-coupled electron-transfer mechanism. Our calculations show that this transition state lies below the sum of the energy of the reactants, and we computed a rate constant at 300 K of 9.29 x 10(-14) cm(3) molecule(-1) s(-1), which is in good agreement with the experimental results. Also of interest in combustion chemistry, we studied the hydrogen abstraction process by HO(2), the result of which is the formation of HCO + H(2)O(2). We found two reaction paths with activation enthalpies close to 12 kcal/mol. For this process, we computed a rate constant of 1.48 x 10(-16) cm(3) molecule(-1) s(-1) at 700 K, which also agrees quite well with experimental results.     

Reaction mechanism between carbonyl oxide and hydroxyl radical: a theoretical study

The reaction mechanism of carbonyl oxide with hydroxyl radical was investigated by using CASSCF, B3LYP, QCISD, CASPT2, and CCSD(T) theoretical approaches with the 6-311+G(d,p), 6-311+G(2df, 2p), and aug-cc-pVTZ basis sets. This reaction involves the formation of H2CO + HO2 radical in a process that is computed to be exothermic by 57 kcal/mol. However, the reaction mechanism is very complex and begins with the formation of a pre-reactive hydrogen-bonded complex and follows by the addition of HO radical to the carbon atom of H2COO, forming the intermediate peroxy-radical H2C(OO)OH before producing formaldehyde and hydroperoxy radical. Our calculations predict that both the pre-reactive hydrogen-bonded complex and the transition state of the addition process lie energetically below the enthalpy of the separate reactants (DeltaH(298K) = -6.1 and -2.5 kcal/mol, respectively) and the formation of the H2C(OO)OH adduct is exothermic by about 74 kcal/mol. Beyond this addition process, further reaction mechanisms have also been investigated, which involve the abstraction of a hydrogen of carbonyl oxide by HO radical, but the computed activation barriers suggest that they will not contribute to the gas-phase reaction of H2COO + HO.     

Theoretical characterization of the gas-phase O(3)HO hydrogen-bonded complex.

We report a theoretical study on two gas-phase hydrogen-bonded complexes formed between ozone and hydroxyl radical that have relevance to atmospheric chemistry. This study was carried out by using CASSCF, CASPT2, QCISD, and CCSD(T) theoretical approaches in conjunction with the 6-311+G(2df,2p) and aug-cc-pVTZ basis sets. Both complexes have a planar structure and differ from each other in the orientation of the electronic density of the unpaired electron associated with the HO radical moiety. Our calculations predict their stabilities to be 0.87 and 0.67 kcal mol(-1), respectively, at 0 K and show the importance of anharmonic effects in computing the red shift of the HO stretch originating from the hydrogen-bonding interaction. We also report two transition states involving the movement of the HO moiety on the potential energy surfaces of these hydrogen-bonded complexes.     

Theoretical studies on reactions of the stabilized H2COO with HO2 and the HO2···H2O complex

The reactions of H(2)COO with HO(2) and the HO(2)···H(2)O complex are studied by employing the high-level quantum chemical calculations with B3LYP and CCSD(T) theoretical methods, the conventional transition-state theory (CTST), and the Rice-Ramsperger-Kassel-Marcus (RRKM) with Eckart tunneling correction. The calculated results show that the proton transfer plus the addition reaction channel (TS1A) is preferable for the reaction of H(2)COO with HO(2) because the barriers are -10.8 and 1.6 kcal/mol relative to the free reactants and the prereactive complex, respectively, at the CCSD(T)/6-311++G(3df,2p)//B3LYP/6-311++G(d,p) level of theory. Furthermore, the rate constant via TS1A (2.23 × 10(-10) cm(3) molecule(-1) s(-1)) combined with the concentrations of the species in the atmosphere demonstrates that the HO(2) radical would be the dominant sink of H(2)COO in some areas, where the concentration of water is less than 10(17) molecules cm(-3). In addition, although the single water molecule would lower the activated barrier of TS1A from 1.0 to 0.1 kcal/mol with respect to the respective complexes, the rate constant is lower than that of the reaction of HO(2) with H(2)COO.     

Ab initio study of the mechanism of the atmospheric reaction: NO2 + O3 --> NO3 + O2

The atmospheric reaction NO2 + O3 --> NO3 + O2 (1) has been investigated theoretically by using the MP2, G2, G2Q, QCISD, QCISD(T), CCSD(T), CASSCF, and CASPT2 methods with various basis sets. The results show that the reaction pathway can be divided in two different parts at the MP2 level of theory. At this level, the mechanism proceeds along two transition states (TS1 and TS2) separated by an intermediate, designated as A. However, when the single-reference higher correlated QCISD methodology has been employed, the minimum A and the transition state TS2 are not found on the hypersurface of potential energy, which confirms a direct reaction mechanism. Single-reference high correlated and multiconfigurational methods consistently predict the barrier height of reaction (1) to be within the range 2.5-6.1 kcal mol(-1), in reasonable agreement with experimental data. The calculated reaction enthalpy is -24.6 kcal mol(-1) and the reaction rate calculated at the highest CASPT2 level, of k = 6.9 x 10(-18) cm(3) molecule(-1) s(-1). Both results can be regarded also as accurate predictions of the methodology employed in this article.     

The "somersault" mechanism for the p-450 hydroxylation of hydrocarbons. The intervention of transient inverted metastable hydroperoxides

A series of model theoretical calculations are described that suggest a new mechanism for the oxidation step in enzymatic cytochrome P450 hydroxylation of saturated hydrocarbons. A new class of metastable metal hydroperoxides is described that involves the rearrangement of the ground-state metal hydroperoxide to its inverted isomeric form with a hydroxyl radical hydrogen bonded to the metal oxide (MO-OH --> MO....HO). The activation energy for this somersault motion of the FeO-OH group is 20.3 kcal/mol for the P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] to produce the isomeric ferryl oxygen hydrogen bonded to an *OH radical [Por(SH)Fe(III)-O....HO(-)]. This isomeric metastable hydroperoxide, the proposed primary oxidant in the P450 hydroxylation reaction, is calculated to be 17.8 kcal/mol higher in energy than the ground-state iron(III) hydroperoxide Cpd 0. The first step of the proposed mechanism for isobutane oxidation is abstraction of a hydrogen atom from the C-H bond of isobutane by the hydrogen-bonded hydroxyl radical to produce a water molecule strongly hydrogen bonded to anionic Cpd II. The hydroxylation step involves a concerted but nonsynchronous transfer of a hydrogen atom from this newly formed, bound, water molecule to the ferryl oxygen with a concomitant rebound of the incipient *OH radical to the carbon radical of isobutane to produce the C-O bond of the final product, tert-butyl alcohol. The TS for the oxygen rebound step is 2 kcal/mol lower in energy than the hydrogen abstraction TS (DeltaE() = 19.5 kcal/mol). The overall proposed new mechanism is consistent with a lot of the ancillary experimental data for this enzymatic hydroxylation reaction.     

On the nature of the unusually long OO bond in HO(3) and HO(4) radicals

The HO(3) and HO(4) polyoxide radicals have attracted some attention due to their potential role in ozone chemistry. Experimentally, the geometrical structure of HO(3) is known whereas that of HO(4) is not. Moreover, the existence of the latter radical has been questioned. Theoretical calculations on the two species have been reported before, showing important structural differences depending on the computational level. Both radicals present an unusually long OO bond (around 1.7-1.8 A) that can be associated with an intricate interaction between HO, or HO(2), with O(2). The nature of such interaction is investigated in detail using large scale ab initio methods (CASSCF, CASPT2, MRCI, QCISD, CCSD(T)) and density functional techniques (B3LYP) in connection with extended basis sets. Stabilization enthalpies at 298 K with respect to HO (or HO(2)) and O(2) have been computed amounting to -3.21 kcal mol(-1) for HO(3) (trans conformation) and 11.33 kcal mol(-1) for HO(4) (cis conformation). The corresponding formation enthalpies are 6.12 and 11.83 kcal mol(-1). The trans conformation of HO(4) is less stable than the cis one by 6.17 kcal mol(-1). Transition states for HO(4) dissociation and for cis/trans conversion are also described.     

Thermal decomposition of ethanol. III. A computational study of the kinetics and mechanism for the CH3+C2H5OH reaction

The mechanism for the CH3+C2H5OH reaction has been investigated by the modified Gaussian-2 method based on the geometric parameters of the stationary points optimized at the B3LYP/6-311+G(d,p) level of theory. Five transition states have been identified for the production of CH4+CH3CHOH (TS1), CH4+CH3CH2O (TS2), CH4+CH2CH2OH (TS3), CH3OH+CH3CH2 (TS4), and CH3CH2OCH3+H (TS5) with the corresponding barriers 12.0, 13.2, 16.0, 44.7, and 49.9 kcal/mol, respectively. The predicted rate constants and branching ratios for the three lower-energy H-abstraction reactions were calculated using the conventional and variational transition state theory with quantum-mechanical tunneling corrections for the temperature range 300-3000 K. The predicted total rate constant, kt=8.36 x 10(-76) T(20.00) exp(5258/T) cm3 mol(-1) s(-1) (300-600 K) and 6.10 x 10(-25) T(4.10)exp(-4058/T) cm3 mol(-1) s(-1) (600-3000 K), agrees closely with existing experimental data in the temperature range 403-523 K. Similarly, the predicted rate constants for CH3+CH3CD2OH and CD3+C2H5OD are also in reasonable agreement with available low temperature kinetic data.     

Hydroxy peroxy nitrites and nitrates from OH initiated reactions of isoprene

The reaction of hydroxy peroxy radicals (RO(2)) with NO represents one of the most crucial tropospheric processes, leading to terrestrial ozone formation or NO(x)() removal and chain termination. We investigate the formation of hydroxy peroxy nitrites (ROONO) and nitrates (RONO(2)) from the OH-isoprene reactions using DFT and ab initio theories and variational RRKM/master equation (vRRKM/ME) formalism. The binding energies of ROONO from NO addition to RO(2) are determined to be in the range of 20-22 kcal mol(-)(1), and the bond dissociation energies of ROONO to form an alkoxy radical (RO) and NO(2) range from 6 to 9 kcal mol(-)(1). Isomerization of ROONO to RONO(2) is exothermic by 22-28 kcal mol(-)(1). The entrance and exit channels of the RO(2)-NO reaction are found to be barrierless, and the rate constants to form ROONO are calculated to be 3 x 10(-)(12) to 2 x 10(-)(11) cm(3) molecule(-)(1) s(-)(1) using the canonical variational transition state theory. The vRRKM/ME analysis reveals negligible stabilization of excited ROONO and provides an assessment of ROONO isomerization to RONO(2).     

Falloff curves for the recombination reaction Cl + FC(O)O + M --> FC(O)OCl + M

The pressure dependence of the recombination reaction Cl + FC(O)O + M --> FC(O)OCl + M has been investigated at 296 K. FC(O)O radicals and Cl atoms were generated by laser flash photodissociation of FC(O)OO(O)CF at 193 nm in mixtures with Cl2 and He or SF6 over the total pressure range 8-645 Torr. The measured FC(O)O radical and F atom yields in the photolysis are 0.33 +/- 0.06 and 0.67 +/- 0.06. The reaction lies in the falloff range approaching the high-pressure limit. The extrapolations toward the limiting low- and high-pressure ranges were carried out using a reduced falloff curves formalism, which includes a recent implementation for the strong-collision broadening factors. The resulting values for the low-pressure rate coefficients are (2.2 +/- 0.4) x 10(-28)[He], (4.9 +/- 0.9) x 10(-28)[SF6], (1.9 +/- 0.3) x 10(-28)[Cl2] and (5.9 +/- 1.1) x 10(-28)[FC(O)OO(O)CF] cm3 molecule(-1) s(-1). The derived high-pressure rate coefficient is (4.4 +/- 0.8) x 10(-11) cm3 molecule(-1) s(-1). For the reaction Cl + FC(O)OCl --> Cl2 + FC(O)O a rate coefficient of (1.6 +/- 0.3) x 10(-11) cm3 molecule(-1) s(-1) was determined. The high-pressure rate coefficient was theoretically interpreted using SACM/CT calculations on an ab initio electronic potential computed at the G3S level of theory. Standard heat of formation values of -99.9 and -102.5 kcal mol(-1) were computed at the G3//B3LYP/6-311++G(3df,3pd) level of theory for cis-FC(O)OCl and trans-FC(O)OCl, respectively. The computed electronic barrier for the conversion between the trans and cis conformers is 8.9 kcal mol(-1). On the basis of the present results, the above reactions are expected to have a negligible impact on stratospheric ozone levels.     

Determination of equilibrium constants for the reaction between acetone and HO2 using infrared kinetic spectroscopy

The reaction between the hydroperoxy radical, HO(2), and acetone may play an important role in acetone removal and the budget of HO(x) radicals in the upper troposphere. We measured the equilibrium constants of this reaction over the temperature range of 215-272 K at an overall pressure of 100 Torr using a flow tube apparatus and laser flash photolysis to produce HO(2). The HO(2) concentration was monitored as a function of time by near-IR diode laser wavelength modulation spectroscopy. The resulting [HO(2)] decay curves in the presence of acetone are characterized by an immediate decrease in initial [HO(2)] followed by subsequent decay. These curves are interpreted as a rapid (<100 μs) equilibrium reaction between acetone and the HO(2) radical that occurs on time scales faster than the time resolution of the apparatus, followed by subsequent reactions. This separation of time scales between the initial equilibrium and ensuing reactions enabled the determination of the equilibrium constant with values ranging from 4.0 × 10(-16) to 7.7 × 10(-18) cm(3) molecule(-1) for T = 215-272 K. Thermodynamic parameters for the reaction determined from a second-law fit of our van't Hoff plot were Δ(r)H°(245) = -35.4 ± 2.0 kJ mol(-1) and Δ(r)S°(245) = -88.2 ± 8.5 J mol(-1) K(-1). Recent ab initio calculations predict that the reaction proceeds through a prereactive hydrogen-bonded molecular complex (HO(2)-acetone) with subsequent isomerization to a hydroxy-peroxy radical, 2-hydroxyisopropylperoxy (2-HIPP). The calculations differ greatly in the energetics of the complex and the peroxy radical, as well as the transition state for isomerization, leading to significant differences in their predictions of the extent of this reaction at tropospheric temperatures. The current results are consistent with equilibrium formation of the hydrogen-bonded molecular complex on a short time scale (100 μs). Formation of the hydrogen-bonded complex will have a negligible impact on the atmosphere. However, the complex could subsequently isomerize to form the 2-HIPP radical on longer time scales. Further experimental studies are needed to assess the ultimate impact of the reaction of HO(2) and acetone on the atmosphere.     

Peroxy radical isomerization in the oxidation of isoprene

We report experimental evidence for the formation of C(5)-hydroperoxyaldehydes (HPALDs) from 1,6-H-shift isomerizations in peroxy radicals formed from the hydroxyl radical (OH) oxidation of 2-methyl-1,3-butadiene (isoprene). At 295 K, the isomerization rate of isoprene peroxy radicals (ISO2?) relative to the rate of reaction of ISO2? + HO2 is k(isom)(295)/(k(ISO2?+HO2)(295)) = (1.2 ± 0.6) x 10(8) mol cm(-3), or k(isom)(295) ? 0.002 s(-1). The temperature dependence of this rate was determined through experiments conducted at 295, 310 and 318 K and is well described by k(isom)(T)/(k(ISO2?+HO2)(T)) = 2.0 x 10(21) exp(-9000/T) mol cm(-3). The overall uncertainty in the isomerization rate (relative to k(ISO2?+HO2)) is estimated to be 50%. Peroxy radicals from the oxidation of the fully deuterated isoprene analog isomerize at a rate ～15 times slower than non-deuterated isoprene. The fraction of isoprene peroxy radicals reacting by 1,6-H-shift isomerization is estimated to be 8-11% globally, with values up to 20% in tropical regions.     

Role of the somersault rearrangement in the oxidation step for flavin monooxygenases (FMO). A comparison between FMO and conventional xenobiotic oxidation with hydroperoxides

Model quantum mechanical calculations presented for C-4a-flavin hydroperoxide (FlHOOH) at the B3LYP/6-311+G(d,p) level suggest a new mechanism for flavoprotein monooxygenase (FMO) oxidation involving a concerted homolytic O-O bond cleavage in concert with hydroxyl radical transfer from the flavin hydroperoxide rather than an S(N)2-like displacement by the substrate on the C-4a-hydroperoxide OOH group. Homolytic O-O bond cleavage in a somersault-like rearrangement of hydroperoxide C-4a-flavinhydroperoxide (1) (FLHO-OH → FLHO···HO) produces an internally hydrogen-bonded HO(?) radical intermediate with a classical activation barrier of 27.0 kcal/mol. Model hydroperoxide 1 is used to describe the transition state for the key oxidation step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH). A comparison of the electron distribution in the transition structures for the PHBH hydroxylation of p-hydroxybenzoic acid (ΔE(?) = 23.0 kcal/mol) with that of oxidation of trimethylamine (ΔE(?) = 22.3 kcal/mol) and dimethyl sulfide (ΔE? = 14.1 kcal/mol) also suggests a mechanism involving a somersault mechanism in concert with transfer of an HO(?) radical to the nucleophilic heteroatom center with a hydrogen transfer back to the FLH-O residue after the barrier is crossed to produce the final product, FLH-OH. In each case the hydroxylation barrier was less than that of the O-O rearrangement barrier in the absence of a substrate supporting an overall concerted process. All three transition structures bear a resemblance to the TS for the comparable hydroxylation of isobutane (ΔE(?) = 29.2 kcal/mol) and for simple Fenton oxidation by aqueous iron(III) hydroperoxides. To our surprise the oxidation of N- and S-nucleophiles with conventional oxidants such as alkyl hydroperoxides and peracids also proceeds by HO(?) radical transfer in a manner quite similar to that for tricyclic hydroperoxide 1. Stabilization of the developing oxyradical produced by somersault rearrangement for concerted enzymatic oxidation with tricyclic hydroperoxide 1 results in a reduced overall activation barrier.     

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.     

Thermodynamics of the hydroxyl radical addition to isoprene

Oxidation of isoprene by the hydroxyl radical leads to tropospheric ozone formation. Consequently, a more complete understanding of this reaction could lead to better models of regional air quality, a better understanding of aerosol formation, and a better understanding of reaction kinetics and dynamics. The most common first step in the oxidation of isoprene is the formation of an adduct, with the hydroxyl radical adding to one of four unsaturated carbon atoms in isoprene. In this paper, we discuss how the initial conformations of isoprene, s-trans and s-gauche, influences the pathways to adduct formation. We explore the formation of pre-reactive complexes at low and high temperatures, which are often invoked to explain the negative temperature dependence of this reaction's kinetics. We show that at higher temperatures the free energy surface indicates that a pre-reactive complex is unlikely, while at low temperatures the complex exists on two reaction pathways. The theoretical results show that at low temperatures all eight pathways possess negative reaction barriers, and reaction energies that range from -36.7 to -23.0 kcal x mol(-1). At temperatures in the lower atmosphere, all eight pathways possess positive reaction barriers that range from 3.8 to 6.0 kcal x mol(-1) and reaction energies that range from -28.8 to -14.4 kcal x mol(-1).     

Mechanism of 1,3-migration in allylperoxyl radicals: computational evidence for the formation of a loosely bound radical-dioxygen complex

The three pathways postulated for 1,3-migration of the peroxyl group in the allylperoxyl radical (1a), a key reaction involved in the spontaneous autoxidation of unsaturated lipids of biological importance, have been investigated by means of quantum mechanical electronic structure calculations. According to the barrier heights calculated from RCCSD(T)/6-311+G(3df,2p) energies with optimized molecular geometries and harmonic vibrational frequencies determined at the UMP2/6-311+G(3df,2p) level, the allylperoxyl rearrangement proceeds by fragmentation of 1a through a transition structure (TS1) with a calculated DeltaH++(298 K) of 21.7 kcal/mol to give an allyl radical-triplet dioxygen loosely bound complex (CX). In a subsequent step, the triplet dioxygen moiety of CX recombines at either end of the allyl radical moiety to convert the complex to the rearranged peroxyl radical (1a') or to revert to the starting peroxyl radical 1a. CX shows an electron charge transfer of 0.026 e in the direction allyl --> O(2). The dominant attractive interactions holding in association the allyl radical-triplet dioxygen pair in CX are due chiefly to dispersion forces. The DeltaH(298 K) for dissociation of CX in its isolated partners, allyl radical and triplet dioxygen, is predicted to be at least 1 kcal/mol. The formation of CX prevents the diffusion of its partners and maintains the stereocontrol along the fragmentation-recombination processes. The concerted 1,3-migration in allylperoxyl radical is predicted to take place through a five-membered ring peroxide transition structure (TS2) showing two long C-O bonds. The DeltaH++(298 K) calculated for this pathway is less favorable than the fragmentation-recombination pathway by 1.9 kcal/mol. The cyclization of 1a to give a dioxolanyl radical intermediate (2a) is found to proceed through a five-membered ring transition structure (TS3) with a calculated DeltaH++(298 K) of 33.9 kcal/mol. Thus, the sequence of ring closure 1a --> 2a and ring opening 2a --> 1a' is unlikely to play any significant role in allylperoxyl rearrangement 1a --> 1a'. In the three pathways investigated, the energy of the transition structure is predicted to be somewhat lower in either heptane or aqueous solution than in the gas phase. Although the energy lowering calculated for TS1 is smaller than the calculated for TS2 and TS3, it is very unlikely that the solvent effects may reverse the predicted preference of the fragmentation-recombination pathway over the concerted and stepwise ring closure-ring opening mechanisms.     

Transition states for the dimerization of 1,3-cyclohexadiene: a DFT, CASPT2, and CBS-QB3 quantum mechanical investigation

Quantum mechanical calculations using restricted and unrestricted B3LYP density functional theory, CASPT2, and CBS-QB3 methods for the dimerization of 1,3-cyclohexadiene (1) reveal several highly competitive concerted and stepwise reaction pathways leading to [4 + 2] and [2 + 2] cycloadducts, as well as a novel [6 + 4] ene product. The transition state for endo-[4 + 2] cycloaddition (endo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 28.7 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 19.0 kcal/mol) is not bis-pericyclic, leading to nondegenerate primary and secondary orbital interactions. However, the C(s) symmetric second-order saddle point on the B3LYP energy surface is only 0.3 kcal/mol above endo-2TS. The activation enthalpy for the concerted exo-[4 + 2] cycloaddition (exo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 30.1 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 21.1 kcal/mol) is 1.4 kcal/mol higher than that of the endo transition state. Stepwise pathways involving diallyl radicals are formed via two different C-C forming transition states (rac-5TS and meso-5TS) and are predicted to be competitive with the concerted cycloaddition. Transition states were located for cyclization from intermediate rac-5 leading to the endo-[4 + 2] (endo-2) and exo-[2 + 2] (anti-3) cycloadducts. Only the endo-[2 + 2] (syn-3) transition state was located for cyclization of intermediate meso-5. The novel [6 + 4] "concerted" ene transition state (threo-4TS, DeltaH(double dagger)(UB3LYP(0K)) = 28.3 kcal/mol) is found to be unstable with respect to an unrestricted calculation. This diradicaloid transition state closely resembles the cyclohexadiallyl radical rather than the linked cyclohexadienyl radical. Several [3,3] sigmatropic rearrangement transition states were also located and have activation enthalpies between 27 and 31 kcal/mol.     

Experimental and theoretical study of reaction of OH with 1,3-butadiene

The kinetics of the reaction of hydroxyl radical with 1,3-butadiene at 240-340 K and a total pressure of approximately 1 Torr has been studied using relative rate combined with the discharge flow and mass spectrometer technique. The reaction dynamics of the same reaction has also been investigated using ab initio molecular orbital theory. The rate constant for this reaction was found to be negatively dependent on temperature, with an Arrhenius expression of k1 = (1.58 +/- 0.07) x 10(-11) exp[(436 +/- 13)/T] cm3 molecule(-1) s(-1) (uncertainties taken as 2sigma), which was in good agreement with that reported by Atkinson et al. and Liu et al. at 299-424 K. Mass spectral evidences were found for the addition of OH to both the terminal and the internal carbons of 1,3-butadiene. Our computational results suggest that both addition of OH to 1,3-butadiene and the abstraction of hydrogen atom from 1,3-butadiene by the OH radical are exothermic processes and that the addition of OH to the terminal carbon of the 1,3-butadiene is predicted to have an activation energy of 0.7 kcal mol(-1), being the most energetically favored reaction pathway.     

Effects of a single water molecule on the OH + H2O2 reaction

The effect of a single water molecule on the reaction between H(2)O(2) and HO has been investigated by employing MP2 and CCSD(T) theoretical approaches in connection with the aug-cc-PVDZ, aug-cc-PVTZ, and aug-cc-PVQZ basis sets and extrapolation to an ∞ basis set. The reaction without water has two elementary reaction paths that differ from each other in the orientation of the hydrogen atom of the hydroxyl radical moiety. Our computed rate constant, at 298 K, is 1.56 × 10(-12) cm(3) molecule(-1) s(-1), in excellent agreement with the suggested value by the NASA/JPL evaluation. The influence of water vapor has been investigated by considering either that H(2)O(2) first forms a complex with water that reacts with hydroxyl radical or that H(2)O(2) reacts with a previously formed H(2)O·OH complex. With the addition of water, the reaction mechanism becomes much more complex, yielding four different reaction paths. Two pathways do not undergo the oxidation reaction but an exchange reaction where there is an interchange between H(2)O(2)·H(2)O and H(2)O·OH complexes. The other two pathways oxidize H(2)O(2), with a computed total rate constant of 4.09 × 10(-12) cm(3) molecule(-1) s(-1) at 298 K, 2.6 times the value of the rate constant of the unassisted reaction. However, the true effect of water vapor requires taking into account the concentration of the prereactive bimolecular complex, namely, H(2)O(2)·H(2)O. With this consideration, water can actually slow down the oxidation of H(2)O(2) by OH between 1840 and 20.5 times in the 240-425 K temperature range. This is an example that demonstrates how water could be a catalyst in an atmospheric reaction in the laboratory but is slow under atmospheric conditions.