Theoretical study of the addition and abstraction reactions of hydroxyl radical with uracil

The addition as well as abstraction reactions of hydroxyl radical (OH) with the nucleic acid base, uracil (U), in the gas phase has been explored at the B3LYP/6-31+G(d,p) level of density functional theory (DFT). The energy barrier of the OH addition to both the C5 and C6 positions of the uracil is less than 1 kcal/mol while the hydrogen abstractions (H-abstractions) from either the N1 or the N3 positions are ∼9.5 kcal/mol. Further the energetics of these reactions are assessed by applying the effect of aqueous medium through the polarizable continuum model (PCM). Both the gas and the solution phase data established that the thermodynamic and kinetic factors are more favorable for the OH addition to either C5 or C6 positions of the uracil than the H-abstraction reactions. Moreover, calculations at the MPW1K/6-31+G(d,p), CCSD(T)/6-31+G(d,p)//B3LYP/6-31+G(d,p) and CCSD(T)/6-31+G(d,p)//MPW1K/6-31+G(d,p) levels of theoretical methods qualitatively supported the B3LYP/6-31+G(d,p) results.     

Gaseous reaction mechanism of C2F radical with water

The kinetic properties of the carbon-fluorine radicals are little understood except those of CFn (n =1-3). In this article, a detailed mechanistic study was reported on the gas-phase reaction between the simplest pi-bonded C2F radical and water as the first attempt to understand the chemical reactivity of the C2F radical. Various reaction channels are considered. The most kinetically competitive channel is the quasi-direct hydrogen-abstraction route forming P5 HCCF + OH. At the CCSD(T)/6-311+G(2d,2p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels, the overall H-abstraction barriers (4.5, 4.7, and 4.2 kcal/mol) for the C2F + H2O reaction are comparable to the corresponding values (5.5, 3.7, and 5.7 kcal/mol) for the analogous C2H + H2O reaction. This suggests that C2F is a reactive radical like the extensively studied C2H, in contrast to the situation of the CF and CF2 radicals that have much lower reactivity than the corresponding hydrocarbon species. Thus, the C2F radical is expected to play an important role in the combustion processes of the carbon-fluorine chemistry. Furthermore, addition of a second H2O can catalyze the reaction with the H-abstraction barrier significantly reduced to a marginally zero value (0.5 kcal/mol). This is also indicative of the potential relevance of the title reactions in the low-temperature atmospheric chemistry.     

Increased reactivity of the *Cr(CO)3(C5Me5) radical with thiones versus thiols: a theoretical and experimental investigation

2-pyridinethione (2-mercaptopyridine, H-2mp) undergoes rapid oxidative addition with 2 mol of the 17-electron organometallic radical *Cr(CO)3Cp (where Cp*=C5Me5), yielding hydride H-Cr(CO)3Cp* and thiolate (eta1-2mp)Cr(CO)3Cp*. In a slower secondary reaction, (eta1-2mp)Cr(CO)3Cp* loses CO generating the N,S-chelate complex (eta2-2mp)Cr(CO)2Cp* for which the crystal structure is reported. The rate of 2-pyridine thione oxidative addition with *Cr(CO)3Cp* (abbreviated *Cr) in toluene best fits rate=kobs[H-2mp][*Cr]; kobs(288 K)=22 +/- 4 M(-1) s(-1); DeltaH++=4 +/- 1 kcal/mol; DeltaS++=- 40 +/- 5 cal/mol K. The rate of reaction is the same under CO or Ar, and the reaction of deuterated 2-pyridine thione (D-2mp) shows a negligible (inverse) kinetic isotope effect (kD/kH=1.06 +/- 0.10). The rate of decarbonylation of (eta1-2mp)Cr(CO)3Cp* forming (eta2-2mp)Cr(CO)2Cp* obeys simple first-order kinetics with kobs (288 K)=3.1x10(-4) s(-1), DeltaH++=23 +/- 1 kcal/mol, and DeltaS++=+ 5.0 +/- 2 cal/mol K. Reaction of 4-pyridine thione (4-mercaptopyridine, H-4mp) with *Cr(CO)3Cp* in THF and CH2Cl2 also follows second-order kinetics and is approximately 2-5 times faster than H-2mp in the same solvents. The relatively rapid nature of the thione versus thiol reactions is attributed to differences in the proposed 19-electron intermediate complexes, [*(S=C5H4N-H)Cr(CO)3Cp*] versus [*(H-S-C6H5)Cr(CO)3Cp*]. In comparison, reactions of pyridyl disulfides occur by a mechanism similar to that followed by aryl disulfides involving direct attack of the sulfur-sulfur bond by the metal radical. Calorimetric data indicate Cr-SR bond strengths for aryl and pyridyl derivatives are similar. The experimental conclusions are supported by B3LYP/6-311+G(3df,2p) calculations, which also provide additional insight into the reaction pathways open to the thione/thiol tautomers. For example, the reaction between H* radical and the 2-pyridine thione S atom yielding a thionyl radical is exothermic by approximately 30 kcal/mol. In contrast, the thiuranyl radical formed from the addition of H* to the 2-pyridine thiol S atom is predicted to be unstable, eliminating either H* or HS* without barrier.     

Effects of OH radical addition on proton transfer in the guanine-cytosine base pair

Double proton transfer (PT) reactions in guanine-cytosine OH radical adducts are studied by the hybrid density functional B3LYP approach. Concerted and stepwise proton-transfer processes are explored between N1(H) on guanine (G) and N3 on cytosine (C), and between N4(H) on C and O6 on G. All systems except GC6OH display a concerted mechanism. 8OHGC has the highest dissociation energy and is 1.2 kcal/mol more stable than the nonradical GC base pair. The origin of the interactions are investigated through the estimation of intrinsic acid-basic properties of the *OH-X monomer (X = G or C). Solvent effects play a significant role in reducing the dissociation energy. The reactions including *OH-C adducts have significantly lower PT barriers than both the nonradical GC pair and the *OH-G adducts. All reactions are endothermic, with the GC6OH --> GC6OHPT reaction has the lowest reaction energy (4.6 kcal/mol). In accordance with earlier results, the estimated NBO charges show that the G moiety carries a slight negative charge (and C a corresponding positive one) in each adduct. The formation of a partial ion pair may be a potential factor leading to the PT reactions being thermodynamically unfavored.     

The role of nucleobase carboradical and carbanion on DNA lesions: a theoretical study

DNA base release induced by H and OH radical addition to thymine and their corresponding electron adducts is studied at the DFT B3LYP/6-31+G(d,p) level in gas phase and in solution. H atom transfer after radical formation from C2' on the sugar to the C6 site on the base is shown to be prohibited for the radical species. Their corresponding electron adducts, albeit minor events in cellular systems, show excellent capabilities to proton transfer from C2' on the sugar to the C6 site on the base. The barriers for subsequent N-glycosidic bond dissociation range from 0.1 to 1.6 kcal mol(-1) at the B3LYP level and around 5 kcal mol(-1) using the BB1K functional, implying that these reactions can serve as a source to abasic sites. Analysis of bond dissociation energies show that all the reactions are exothermic, which is consistent with the changes in N-glycosidic bond lengths during the proton-transfer reactions. Bulk solvation plays a reverse influence on proton transfer and the bond rupture reactions. Molecular orbitals, NPA charges, and electron affinities are calculated to shed further light on the properties leading up to the intramolecular reactions.     

Theoretical study of the antioxidant properties of pyridoxine

Molecules acting as antioxidants capable of scavenging reactive oxygen species (ROS) are of utmost importance in the living cell. The antioxidative properties of pyridoxine (vitamin B6) have recently been discovered. In this study, we have analyzed the reactivity of pyridoxine toward the ROS (.-)OH, (.-)OOH, and (.-)O(2)- at the density functional theory level (functionals B3LYP and MPW1B95). Two reaction types have been studied as follows: addition to the aromatic ring atoms and hydrogen/proton abstraction. Our results show that (.-)OH is the most reactive species, while (.-)OOH displays low reactivity and (.-)O2(-) does not react at all with pyridoxine. The most exergonic reactions are those where (.-)H is removed from the CH(2)OH groups or the ring-bound OH group and range from -33 to -39 kcal/mol. The most exergonic addition reactions occur by attacking the carbon atoms bonded to nitrogen but with an energy gain of only 6 kcal/mol.     

Thermodynamic and kinetic study of ibuprofen with hydroxyl radical: A density functional theory approach

Ibuprofen, a frequently detected pharmaceutical in natural and engineered waters, was studied in both neutral and anionic forms using density functional theory at the B3LYP/6‐311++G**//B3LYP/6‐31G* level of theory in its reaction with hydroxyl radical ( ? OH). The reaction pathways included ? OH addition to aromatic ring, abstraction of a H‐atom, and nucleophilic attack on the carbonyl group. The results showed that H‐atom abstraction pathways are the most favorable. The free energy change for H‐atom abstraction reaction ranges from ?37.8 to ?15.9 kcal/mol; for ? OH addition ranges from ?3.85 to ?1.23 kcal/mol; and for nucleophilic attack on the carbonyl group is 13.9 kcal/mol. The calculated rate constant between neutral ibuprofen and ? OH, 6.72 × 10⁹ M^?1s^?1, is consistent with the experimental value, 6.5 ± 0.2 × 10⁹ M^?1s^?1. Our results provide direct evidence for byproduct formation and identification on the molecular level. © 2013 Wiley Periodicals, Inc.     

Gas-phase kinetics study of reaction of OH radical with CH3NHNH2 by second-order multireference perturbation theory

The gas-phase kinetics of H-abstraction reactions of monomethylhydrazine (MMH) by OH radical was investigated by second-order multireference perturbation theory and two-transition-state kinetic model. It was found that the abstractions of the central and terminal amine H atoms by the OH radical proceed through the formation of two hydrogen bonded preactivated complexes with energies of 6.16 and 5.90 kcal mol(-1) lower than that of the reactants, whereas the abstraction of methyl H atom is direct. Due to the multireference characters of the transition states, the geometries and ro-vibrational frequencies of the reactant, transition states, reactant complexes, and product complexes were optimized by the multireference CASPT2/aug-cc-pVTZ method, and the energies of the stationary points of the potential energy surface were refined at the QCISD(T)/CBS level via extrapolation of the QCISD(T)/cc-pVTZ and QCISD(T)/cc-pVQZ energies. It was found that the abstraction reactions of the central and two terminal amine H atoms of MMH have the submerged energy barriers with energies of 2.95, 2.12, and 1.24 kcal mol(-1) lower than that that of the reactants respectively, and the abstraction of methyl H atom has a real energy barrier of 3.09 kcal mol(-1). Furthermore, four MMH radical-H(2)O complexes were found to connect with product channels and the corresponding transition states. Consequently, the rate coefficients of MMH + OH for the H-abstraction of the amine H atoms were determined on the basis of a two-transition-state model, with the total energy E and angular momentum J conserved between the two transition-state regions. In units of cm(3) molecule(-1) s(-1), the rate coefficient was found to be k(1) = 3.37 × 10(-16)T(1.295) exp(1126.17/T) for the abstraction of the central amine H to form the CH(3)N(?)NH(2) radical, k(2) = 2.34 × 10(-17)T(1.907) exp(1052.26/T) for the abstraction of the terminal amine H to form the trans-CH(3)NHN(?)H radical, k(3) = 7.41 × 10(-20)T(2.428) exp(1343.20/T) for the abstraction of the terminal amine H to form the cis-CH(3)NHN(?)H radical, and k(4) = 9.13 × 10(-21)T(2.964) exp(-114.09/T) for the abstraction of the methyl H atom to form the C(?)H(2)NHNH(2) radical, respectively. Assuming that the rate coefficients are additive, the total rate coefficient of these theoretical predictions quantitatively agrees with the measured rate constant at temperatures of 200-650 K, with no adjustable parameters.     

Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1,3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions

Activation barriers and reaction energetics for the three main classes of 1,3-dipolar cycloadditions, including nine different reactions, were evaluated with the MPW1K and B3LYP density functional methods, MP2, and the multicomponent CBS-QB3 method. The CBS-QB3 values were used as standards for 1,3-dipolar cycloaddition activation barriers and reaction energetics, and the density functional theory (DFT) and MP2 methods were benchmarked against these values. The MPW1K/6-31G* method and basis set performs best for activation barriers, with a mean absolute deviation (MAD) value of 1.1 kcal/mol. The B3LYP/6-31G* method and basis set performs best for reaction enthalpies, with a MAD value of 2.4 kcal/mol, while the MPW1K method shows large errors for reaction energetics. The MP2 method gives the expected systematic underestimation of barriers. Concerted and nearly synchronous transition structures are predicted by all DFT and MP2 methods. Also reported are revised estimated 0 K experimental activation enthalpies for a standard set of hydrocarbon pericyclic reactions and updated comparisons to experiment for DFT, ab initio, and multicomponent methods. B3LYP and MPW1K methods with MAD values of 1.5 and 2.1 kcal/mol, respectively, fortuitously outperform the multicomponent CBS-QB3 method, which has a MAD value of 2.3. The MAD value of the O3LYP functional improves to 2.4 kcal/mol from the previously reported 3.0 kcal/mol.     

Mechanisms of formation of 8-oxoguanine due to reactions of one and two OH* radicals and the H2O2 molecule with guanine: A quantum computational study

Mechanisms of formation of the mutagenic product 8-oxoguanine (8OG) due to reactions of guanine with two separate OH* radicals and with H2O2 were investigated at the B3LYP/6-31G, B3LYP/6-311++G, and B3LYP/AUG-cc-pVDZ levels of theory. Single point energy calculations were carried out with the MP2/AUG-cc-pVDZ method employing the optimized geometries at the B3LYP/AUG-cc-pVDZ level. Solvent effect was treated using the PCM and IEF-PCM models. Reactions of two separate OH* radicals and H2O2 with the C2 position of 5-methylimidazole (5MI) were investigated taking 5MI as a model to study reactions at the C8 position of guanine. The addition reaction of an OH* radical at the C8 position of guanine is found to be nearly barrierless while the corresponding adduct is quite stable. The reaction of a second OH* radical at the C8 position of guanine leading to the formation of 8OG complexed with a water molecule can take place according to two different mechanisms, involving two steps each. According to one mechanism, at the first step, 8-hydroxyguanine (8OHG) complexed with a water molecule is formed ,while at the second step, 8OHG is tautomerized to 8OG. In the other mechanism, at the first step, an intermediate complexed (IC) with a water molecule is formed, the five-membered ring of which is open, while at the second step, the five-membered ring is closed and a hydrogen bonded complex of 8OG with a water molecule is formed. The reaction of H2O2 with guanine leading to the formation of 8OG complexed with a water molecule can also take place in accordance with two different mechanisms having two steps each. At the first step of one mechanism, H2O2 is dissociated into two OH* groups that react with guanine to form the same IC as that formed in the reaction with two separate OH* radicals, and the subsequent step of this mechanism is also the same as that of the reaction of guanine with two separate OH* radicals. At the first step of the other mechanism of the reaction of guanine with H2O2, the latter molecule is dissociated into a hydrogen atom and an OOH* group which become bonded to the N7 and C8 atoms of guanine, respectively. At the second step of this mechanism, the OOH* group is dissociated into an oxygen atom and an OH* group, the former becomes bonded to the C8 atom of guanine while the latter abstracts the H8 atom bonded to C8, thus producing 8OG complexed with a water molecule. Solvent effects of the aqueous medium on certain reaction barriers and released energies are appreciable. 5MI works as a satisfactory model for a qualitative study of the reactions of two separate OH* radicals or H2O2 occurring at the C8 position of guanine.     

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.     

Significant effects of phosphorylation on relative stabilities of DNA and RNA sugar radicals: remarkably high susceptibility of h-2' abstraction in RNA

The roles of nucleic acid radicals in DNA and RNA damage cannot be properly understood in the absence of knowledge of the C-H bond strengths depicting the energy cost to generate each of these radicals. However, previous theoretical studies on the relative energies of different nucleic acid radicals are not fully convincing mainly because of the use of oversimplified model compounds. In the present study we chose nucleoside 3',5'-bisphosphates as model compounds for DNA and RNA, in which the effects of both the nucleobase and phosphorylation were taken into consideration. Using the newly developed ONIOM-G3B3 methods, we calculated the gas-phase bond dissociation enthalpies and solution-phase bond dissociation free energies of all the carbohydrate C-H bonds in the model compounds. It was found that the monoanionic phosphate group (OPO3H-) was a better radical stabilization group than the OH group by 1.3 kcal/mol, whereas the neutral phosphate group (OPO3H2) was a significantly worse radical stabilization group than OH by 4.4 kcal/mol. Due to these reasons, the relative thermodynamic susceptibility of H-abstraction from deoxyribonucleotides and ribonucleotides varied considerably depending on the phosphorylation state and the charge carried by the phosphate groups. Strikingly, the bond dissociation free energy of C2'-H in ribonucleotides was dramatically lower than that of all the other C-H bonds by 5-6 kcal/mol regardless of the phosphorylation state and the charge carried by the phosphate group. This explained the previous experimental finding that radiation damage of RNA occurs mainly via H-abstraction at H-2'. A model study suggested that the strength of the hydrogen bonding interaction between the 2'-OH and 3-phosphate groups should dramatically increase from ribonucleoside 3',5'-bisphosphate to its C2' radical. The strengthened hydrogen bonding stabilized the C2' radical, rendering the C2'-H bond of RNA extraordinarily vulnerable to H-abstraction.     

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.     

Solution calorimetric and stopped-flow kinetic studies of the reaction of *Cr(CO)3C5Me5 with PhSe-SePh and PhTe-TePh. Experimental and theoretical estimates of the Se-Se, Te-Te, H-Se, and H-Te bond strengths

The kinetics of the oxidative addition of PhSeSePh and PhTeTePh to the stable 17-electron complex *Cr(CO)3C5Me5 have been studied utilizing stopped-flow techniques. The rates of reaction are first-order in each reactant, and the enthalpy of activation decreases in going from Se (deltaH(double dagger) = 7.0 +/- 0.5 kcal/mol, deltaS(double dagger) = -22 +/- 3 eu) to Te (deltaH(double dagger) = 4.0 +/- 0.5 kcal/mol, deltaS(double dagger) = -26 +/- 3 eu). The kinetics of the oxidative addition of PhSeH and *Cr(CO)3C5Me5 show a change in mechanism in going from low (overall third-order) to high (overall second-order) temperatures. The enthalpies of the oxidative addition of PhE-EPh to *Cr(CO)3C5Me5 in toluene solution have been measured and found to be -29.6, -30.8, and -28.9 kcal/mol for S, Se, and Te, respectively. These data are combined with enthalpies of activation from kinetic studies to yield estimates for the solution-phase PhE-EPh bond strengths of 46, 41, and 33 kcal/mol for E = S, Se, and Te, respectively. The corresponding Cr-EPh bond strengths are 38, 36, and 31 kcal/mol. Two methods have been used to determine the enthalpy of hydrogenation of PhSeSePh in toluene on the basis of reactions of HSPh and HSePh with either *Cr(CO)3C5Me5 or 2-pyridine thione. These data lead to a thermochemical estimate of 72 kcal/mol for the PhSe-H bond strength in toluene solution, which is in good agreement with kinetic studies of H atom transfer from HSePh at higher temperatures. The reaction of H-Cr(CO)3C5Me5 with PhSe-SePh is accelerated by the addition of a Cr radical and occurs via a rapid radical chain reaction. In contrast, the reaction of PhTe-TePh and H-Cr(CO)3C5Me5 does not occur at any appreciable rate at room temperature, even in the presence of added Cr radicals. This is in keeping with a low PhTe-H bond strength blocking the chain and implies that H-TePh < or = 63 kcal/mol. Structural data are reported for PhSe-Cr(CO)3C5Me5 and PhS-Cr(CO)3C5Me5. The two isostructural complexes do not show signs of an increase in steric strain in terms of metal-ligand bonds or angles as the Cr-EPh bond is shortened in going from Se to S. Bond strength estimates of the PhE-H and PhE-EPh derived from density functional theory calculations are in reasonable agreement with experimental data for E = Se but not for E = Te. The nature of the singly occupied molecular orbital of the *EPh radicals is calculated to show increasing localization on the chalcogenide atom in going from S to Se to Te.     

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.     

Kinetics of hydrogen abstraction reactions of butene isomers by OH radical

The rate coefficients of H-abstraction reactions of butene isomers by the OH radical were determined by both canonical variational transition-state theory and transition-state theory, with potential energy surfaces calculated at the CCSD(T)/6-311++G(d,p)//BH&HLYP/6-311G(d,p) level and CCSD(T)/6-311++G(d,p)//BH&HLYP/cc-pVTZ level and quantum mechanical tunneling effect corrected by either the small-curvature tunneling method or the Eckart method. While 1-butene contains allylic, vinylic, and alkyl hydrogens that can be abstracted to form different butene radicals, results reveal that s-allylic H-abstraction channels have low and broad energy barriers, and they are the most dominant channels which can occur via direct and indirect H-abstraction channels. For the indirect H-abstraction s-allylic channel, the reaction can proceed via forming two van der Waals prereactive complexes with energies that are 2.7-2.8 kcal mol(-1) lower than that of the entrance channel at 0 K. Assuming that neither mixing nor crossover occurs between different reaction pathways, the overall rate coefficient was calculated by summing the rate coefficients of the s-allyic, methyl, and vinyl H-abstraction paths and found to agree well with the experimentally measured OH disappearance rate. Furthermore, the rate coefficients of p-allylic H abstraction of cis-2-butene, trans-2-butene, and isobutene by the OH radical were also determined at 300-1500 K, with results analyzed and compared with available experimental data.     

Kinetics of alpha-hydroxy-alkylperoxyl radicals in oxidation processes. HO2*-initiated oxidation of ketones/aldehydes near the tropopause

A comparative theoretical study is presented on the formation and decomposition of alpha-hydroxy-alkylperoxyl radicals, Q(OH)OO* (Q = RR'C:), important intermediates in the oxidation of several classes of oxygenated organic compounds in atmospheric chemistry, combustion, and liquid-phase autoxidation of hydrocarbons. Detailed potential energy surfaces (PESs) were computed for the HOCH2O2* <==>HO2* + CH2O reaction and its analogues for the alkyl-substituted RCH(OH)OO* and R2C(OH)OO* and the cyclic cyclo-C6H10(OH)OO*. The state-of-the-art ab initio methods G3 and CBS-QB3 and a nearly converged G2M//B3LYP-DFT variant were found to give quasi-identical results. On the basis of the G2M//B3LYP-DFT PES, the kinetics of the approximately equal to 15 kcal/mol endothermal alpha-hydroxy-alkylperoxyl decompositions and of the reverse HO2*+ ketone/aldehyde reactions were evaluated using multiconformer transition state theory. The excellent agreement with the available experimental (kinetic) data validates our methodologies. Contrary to current views, HO2* is found to react as fast with ketones as with aldehydes. The high forward and reverse rates are shown to lead to a fast Q(OH)OO* <==>HO2* + carbonyl quasi-equilibrium. The sizable [Q(OH)OO*]/[carbonyl] ratios predicted for formaldehyde, acetone, and cyclo-hexanone at the low temperatures (below 220 K) of the earth's tropopause are shown to result in efficient removal of these carbonyls through fast subsequent Q(OH)OO* reactions with NO and HO2*. IMAGES model calculations indicate that at the tropical tropopause the HO2*-initiated oxidation of formaldehyde and acetone may account for 30% of the total removal of these major atmospheric carbonyls, thereby also substantially affecting the hydroxyl and hydroperoxyl radical budgets and contributing to the production of formic and acetic acids in the upper troposphere and lower stratosphere. On the other hand, an RRKM-master equation analysis shows that hot alpha-hydroxy-alkylperoxyls formed by the addition of O(2) to C(1)-, C(2)-, and C(3)-alpha-hydroxy-alkyl radicals will quasi-uniquely fragment to HO2* plus the carbonyl under all atmospheric conditions.     

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.     

DFT studies of the phosphinidene derivative HPNaF and its insertion reaction with R–H(R=F,OH, NH<Subscript>2</Subscript>, CH<Subscript>3</Subscript>)

The formations of the phosphinidene derivative HPNaF and its insertion reactions with R–H (R=F, OH, NH₂, CH₃) have been systematically investigated employing the density functional theory (DFT), such as the B3LYP and MPW1PW91 methods. A comparison with the results of MP2 calculations shows that MPW1PW91 underestimates the barrier heights for the four reactions considered. Similarly, the same is also true for the B3LYP method depending on the selected reactions, but by much less than MPW1PW91, where the barrier heights of the four reactions are 25.2, 85.7, 119.0, and 142.4 kJ/mol at the B3LYP/6-311+G* level of theory, respectively. All the mechanisms of the four reactions are identical to each other, i.e., an intermediate has been located during the insertion reaction. Then, the intermediate could dissociate to substituted phosphinidane(H₂RP) and NaF with a barrier corresponding to their respective dissociation energies. Correspondingly, the reaction energies for the four reactions are −92.2, −68.1, −57.2, and −44.3 kJ/mol at the B3LYP/6-311+G* level of theory, respectively, where both the B3LYP and MPW1PW91 methods underestimate the reaction energies compared with the MP2 results. The linear correlations between the calculated barrier heights and the reaction energies have also been observed. As a result, the relative reactivity among the four insertion reactions should be as follows: H–F > H–OH > H–NH₂ > H–CH₃.     

Hydrogen atom transfer reactions of iron-porphyrin-imidazole complexes as models for histidine-ligated heme reactivity

Hydrogen atom transfer (HAT) reactions of the bis(histidine) cytochrome active site models (TPP)FeII(ImH)2 (FeIIImH) and (TPP)Fe(Im)(ImH) (FeIIIIm) have been examined in acetonitrile solvent (TPP = tetraphenylporphyrin, ImH = 4-methylimidazole). The ascorbate derivative 5,6-isopropylidine ascorbate, hydroquinone, and the hydroxylamine TEMPOH all rapidly add H* to FeIIIIm to give FeIIImH. Similarly, the phenoxyl radical 2,4,6-tBu3C6H2O* and excess TEMPO* each oxidize FeIIImH to give FeIIIIm. On the basis of redox potential, pKa, and equilibrium measurements, the N-H bond in FeIIImH was found to have a bond dissociation free energy (BDFE) of 70 +/- 2 kcal mol(-1). A hydrogen atom transfer mechanism (concerted transfer of e- and H+) is indicated based on data for the ascorbate and TEMPO* reactions.