Hydrogen abstraction from copolymers of fluorinated olefins and allyl or vinyl ethers by hydroxyl radicals: a computational quantum semi-empirical study

The theoretical study of hydrogen abstraction by hydroxyl radicals on two substrates (copolymers of fluorinated olefins and allyl or vinyl ethers) was carried using the MNDO, AM1 and PM3 quantum semi-empirical methods. This study was performed as a function of the site of hydrogen abstraction and of the computational method. The results of the calculations clearly show that the transition state is early along the reaction coordinate and pinpoint that the reactions are not under enthalpic control. The results provide evidence of the importance of the polar effects due to the fluorine atoms.     

Hydrogen abstraction from poly(propylene) and poly(propylene oxide) by hydroxyl radicals: a computational quantum semi-empirical study

The reaction of hydrogen abstraction by hydroxyl radicals on two substrates (poly(propylene) and poly(propylene oxide)) was studied using three quantum semi-empirical methods (MNDO, AM1, PM3). The calculations were performed as a function of the site of abstraction (hydrogen atom on a secondary or tertiary carbon atom), and of the calculation method. In each case, we localised the transition state and showed that this transition state occurs early along the reaction coordinate. The results concerning the activation energies depend on the sites and the calculation methods. The calculated results were compared to experimental ones.     

interaction of nitro compounds with radicals: a quantum chemical study

Photoreduction of nitro compounds is accompanied by formation of various radical products that can react with the starting nitro compound, thus causing deviation of the decomposition kinetics from the first-order kinetics with respect to the nitro compound. The results of quantum chemical modeling of the reactions of nitro compounds with radicals and the pathways of further transformations of radical adducts formed in the reactions are presented. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 2, pp. 202–206, February, 2006.     

Reactivities of hydroxyl and perhydroxyl radicals toward cytosine and thymine: A comparative study

As the hydroxyl (OH) and perhydroxyl (OOH) radicals are known to play important roles in biological systems, their reactions with cytosine and thymine were studied. Addition reactions of these radicals at different sites of cytosine and thymine, and hydrogen abstraction reactions by each of the two radicals from the different sites of the two molecules were studied at the B3LYP/6‐31G(d,p), B3LYP/AUG‐cc‐pVDZ and BHandHLYP/AUG‐cc‐pVDZ levels of density functional theory. Effect of solvation in aqueous media on the reactions was studied at all these levels of theory using single point energy calculations using the polarizable continuum model. The present study shows that whereas the OH radical would abstract H atoms from the various sites of cytosine and thymine efficiently, the OOH radical would have poor reactivity in this regard. The OH radical is also predicted to be much more reactive than the OOH radical with regard to addition reactions at the C5 and C6 sites of both thymine and cytosine, though the OOH radical is also predicted to have significant reactivity in this respect. © 2012 Wiley Periodicals, Inc.     

Reactivities of superoxide and hydroperoxyl radicals with disubstituted cyclic nitrones: a DFT study

The unique ability of nitrone spin traps to detect and characterize transient free radicals by electron paramagnetic resonance (EPR) spectroscopy has fueled the development of new spin traps with improved properties. Among a variety of free radicals in chemical and biological systems, superoxide radical anion (O(2)(?-)) plays a critical role as a precursor to other more oxidizing species such as hydroxyl radical (HO(?)), peroxynitrite (ONOO(-)), and hypochlorous acid (HOCl), and therefore the direct detection of O(2)(?-) is important. To overcome the limitations of conventional cyclic nitrones, that is, poor reactivity with O(2)(?-), instability of the O(2)(?-) adduct, and poor cellular target specificity, synthesis of disubstituted nitrones has become attractive. Disubstituted nitrones offer advantages over the monosubstituted ones because they allow bifunctionalization of spin traps, therefore accommodating all the desired spin trap properties in one molecular design. However, because of the high number of possible disubstituted analogues as candidate, a systematic computational study is needed to find leads for the optimal spin trap design for biconjugation. In this paper, calculation of the energetics of O(2)(?-) and HO(2)(?) adduct formation from various disubstituted nitrones at PCM/B3LYP/6-31+G(d,p)//B3LYP/6-31G(d) level of theory was performed to determine the most favorable disubstituted nitrones for this reaction. In addition, our results provided general trends of radical reactivity that is dependent upon but not exclusive to the charge densities of nitronyl-C, the position of substituents including stereoselectivities, and the presence of intramolecular H-bonding interaction. Unusually high exoergic ΔG(298K,aq)'s for O(2)(?-) and HO(2)(?) adduct formation were predicted for (3S,5S)-5-methyl-3,5-bis(methylcarbamoyl)-1-pyrroline N-oxide (11-cis) and (4S,5S)-5-dimethoxyphosphoryl-5-methyl-4-ethoxycarbonyl-1-pyrroline N-oxide (29-trans) with ΔG(298K,aq) = -3.3 and -9.4 kcal/mol, respectively, which are the most exoergic ΔG(298K,aq) observed thus far for any nitrone at the level of theory employed in this study.     

Radical ions and radicals in electron-irradiated acrylates and N-isopropylacrylamide: a quantum chemical study

The semiempirical quantum chemical methods MNDO, AM1 and PM3 were used to investigate the performance of the single excited configuration interaction (SCI) approximation for calculating low energy excitation energies of open-shell systems. Systematic calculations were done for eight radicals formed by reactions of H√, OH√ and e⁻_aq with various acrylates and N-isopropylacrylamide. The calculated electronic spectra show a reasonable correlation with experimental data for both neutral radicals and radical ions. The AM1 as well as the PM3 formalism can be successfully applied to calculate the low energy excited states of these types of open shell systems. The best correlation between experimental and calculated excitation energies was obtained using the PM3 method (correlation coefficient 0.96, overall average error 0.16 eV).     

Interaction between methyl and hydroxyl radicals: a theoretical study

In this paper, we conduct a computational quantum mechanic study of two molecules interaction between methyl (·CH₃) and hydroxyl (·OH) radicals. The molecular exploration has been focused on the possibility of finding non-bonding interactions (formation of complexes by weak-bond interaction among C–O–H atoms) and the potential energy reaction between those two molecules. It shows that the formed complexes presents in slightly lower potential energy than that of the reactants and/or products. The existence of these complexes could proceed to the further interaction, e.g. reaction of the molecules, as those complexes formed a particular configuration.     

Reaction of ethylene with hydroxyl radicals: a theoretical study

Ab initio calculations of portions of the C2H5O potential energy surface critical to the title reaction are presented. These calculations are based on QCISD geometries and frequencies and RQCISD(T) energies extrapolated to the complete-basis-set limit. Rate coefficients for the reaction of C2H4 with OH are calculated using this surface and the two transition-state model of Greenwald and co-workers [J. Phys. Chem. A 2005, 109, 6031] for the association of OH with C2H4. The present calculations reproduce most of the experimental data, including the temperature and pressure dependence of the rate coefficients, with only a small (0.4 kcal/mol) adjustment to the energy barrier for direct hydrogen abstraction. We confirm the importance of this channel above 800 K and find that a significant fraction of the total rate coefficient (approximately 10%) is due to the formation of vinyl alcohol above this temperature. Calculations of the vinyl alcohol channel are consistent with the recent observation of this molecule in low-pressure flames [Taatjes, C. A.; Hansen, N.; McIlroy, A.; Miller, J. A.; Senosiain, J. P.; Klippenstein, S. J.; Qi, F.; Sheng, L.; Zhang, Y.; Cool, T. A.; Wang, J.; Westmoreland, P. R.; Law, M. E.; Kasper, T.; Kohse-H?inghaus, K. Science 2005, 308, 1887] and suggest that this reaction should be included in hydrocarbon oxidation mechanisms.     

Theoretical study on the mechanism and kinetics of addition of hydroxyl radicals to fluorobenzene

Geometries, frequencies, reaction barriers, and reaction rates were calculated for the addition of OH radical to fluorobenzene using Möller–Plesset second‐order perturbation (MP2) and G3 methods. Four stationary points were found along each reaction path: reactants, prereaction complex, transition state, and product. A potential for association of OH radical and fluorobenzene into prereaction complex was calculated, and the associated transition state was determined for the first time. G3 calculations give higher reaction barriers than MP2, but also a significantly deeper prereaction complex minimum. The rate constants, calculated with Rice–Ramsperger–Kassel–Marcus (RRKM) theory using G3 energies, are much faster and in much better agreement with the experiment than those calculated with MP2 method, as the deeper well favors the formation of prereaction complex and also increases the final relative populations of adducts. The discrepancies between the experimental and calculated rate constants are attributed to the errors in calculated frequencies as well as to the overestimated G3 reaction barriers and underestimated prereaction complex well depth. It was possible to rectify those errors and to reproduce the experimental reaction rates in the temperature range 230–310 K by treating the relative translation of OH radical and fluorobenzene as a two‐dimensional particle‐in‐the‐box approximation and by downshifting the prereaction complex well and reaction barriers by 0.7 kcal mol^?1. The isomeric distribution of fluorohydroxycyclohexadienyl radicals is calculated from the reaction rates to be 30.9% ortho, 22.6% meta, 38.4% para, and 8.3% ipso. These results are in agreement with experiment that also shows dominance of ortho and para channels. © 2012 Wiley Periodicals, Inc.     

Growing graphene sheets from reactions with methyl radicals: a quantum chemical study.

Hydrogen abstraction reactions by methyl radicals on the zigzag and armchair edges of perylene are studied by density functional theory (DFT) to explore various growth pathways that seem to be in line with experimental observations. The DFT approach is validated by comparing the results obtained from calculations with six different functionals with those obtained from correlated ab initio methods, thereby emphasizing the calculation of reaction barriers. A useful compromise between accuracy and computational cost is provided by DFT, and possible pathways are studied in detail at this level of calculation. Our computational study is carried out by ordering, as a first step, all of the isomers that arise from the abstraction of one or two H atoms from 1,12-dimethyl-1,12-dihydroperylene and 3,4-dimethyl-3,4-dihydroperylene with respect to their energies. Subsequently, only those pathways that connect low-energy isomers are investigated. The calculations reveal that the selected pathways are favored thermodynamically, and also that the reaction barriers are somewhat higher than the energy locally available for the respective reaction. Notably, in the case of 3,4-dimethyl-3,4-dihydroperylene, the first two reaction steps have no or only a very low reaction barrier. The final conclusion of our study is that a cascade of reactions is possible that leads to the growth of a graphene sheet on a graphite surface.     

A computational study of the reaction kinetics of methyl radicals with trifluorohalomethanes

Ab initio calculations have been used to characterize the transition states for halogen abstraction by CH₃ in reactions with CF₄, CF₃Cl, CF₃Br, and CF₃I (1–4). Geometries and frequencies were obtained at the HF/6-31G(d) and MP2=full/6-31G(d) levels of theory. Energy barriers were computed via the Gaussian-2 methodology, and the results were employed in transition state theory analyses to obtain the rate constants over 298–2500 K. There is good accord with literature measurements in the approximate temperature range 360–500 K for reactions (2–4), and the computed activation energies are accurate to within ±6 kJ mol⁻¹. Recommended rate constant expressions for use in combustion modeling are k;₁=1.6×10⁻¹⁹ (T/K)^2.41 exp(−13150 K/T), k₂=8.4×10⁻²⁰(T/K)^2.34 exp(−5000 K/T), k₃=4.6×10⁻¹⁹ (T/K)^2.05 exp(−3990 K/T), and k₄=8.3×10⁻¹⁹ (T/K)^2.18 exp(−1870 K/T) cm³ molecule⁻¹ s⁻¹. The results are discussed in the context of flame suppression chemistry. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet 30: 179–184, 1998.     

On the reactivity of diethyl hydroxyl amine toward free radicals

Diethyl hydroxyl amine is an efficient trap for alkyl, alkoxy, and peroxy radicals. The specific rate constant for the reaction of ethyl radicals (gas phase, 25°C), tert-butoxy radicals (benzene solution, 115°C), and poly (peroxystyryl) peroxy radicals (styrene solution, 50°C) were evaluated as 7.2 × 10⁵, 7.7 × 10⁷, and 2.9 × 10⁵ M^?1·sec^?1, respectively. Several possible secondary reactions of the nitroxide radicals are discussed.     

The association reaction between C2H and 1-butyne: a computational chemical kinetics study

The potential energy surfaces (PES) for the reaction of the C(2)H radical with 1-butyne (C(4)H(6)) have been studied using the CBS-QB3 method. Density functional B3LYP/cc-pVTZ and M06-2X/6-311++G(d,p) calculations have also been performed to analyze the reaction energetics. For detailed theoretical calculation on the total reaction mechanism, the initial association reactions on more and less substituted C atoms of 1-butyne are treated separately followed by a variational transition state theory (VTST) calculation to obtain reaction rates. The successive unimolecular reactions from the association reaction complexes are subjected to Rice-Ramsperger-Kassel-Marcus (RRKM) calculations for reaction rate constants and product branching ratios. The calculated rate constants in the temperature range 70-295 K for both the association reactions are found to be highly temperature dependent at low temperatures, which is contrary to the experimental findings of temperature independent association rates. We have explained this observation with the help of variational nature of the transition states, and we found a "loose" transition state at low temperatures. The calculated product branching ratios for the unimolecular reactions generally agree with the available experimental data, although some channels show a significant method dependency and therefore the correlation with experiment is lost to some extent. Our detailed reaction energetics calculations confirm that the C(2)H + C(4)H(6) reaction proceeds without an entrance barrier and leads to the important products ethynylallene + CH(3), 1,3-hexadiyne + H, 3,4-hexadiene-1-yne + H, 2-ethynyl-1,3-butadiene + H, 3,4-dimethylenecyclobut-1-ene + H and fulvene + H exothermic by 25-75 kcal mol(-1), with strong dependence of the product distribution on the association mode of C(2)H with C(4)H(6), making these reactions fast under low temperature conditions of Titan's atmosphere. Therefore this study can provide a detailed picture of the complex hydrocarbon formation mechanism in the upper atmosphere.     

Reaction of hydroxyl radicals with azacytosines: a pulse radiolysis and theoretical study

Pulse radiolysis and density functional theory (DFT) calculations at B3LYP/6-31+G(d,p) level have been carried out to probe the reaction of the water-derived hydroxyl radicals (*OH) with 5-azacytosine (5Ac) and 5-azacytidine (5Acyd) at near neutral and basic pH. A low percentage of nitrogen-centered oxidizing radicals, and a high percentage of non-oxidizing carbon-centered radicals were identified based on the reaction of transient intermediates with 2,2'-azinobis(3-ethylbenzothiazoline-6-sulfonate), ABTS2-. Theoretical calculations suggests that the N3 atom in 5Ac is the most reactive center as it is the main contributor of HOMO, whereas C5 atom is the prime donor for the HOMO of cytosine (Cyt) where the major addition site is C5. The order of stability of the adduct species were found to be C6-OH_5Ac*>C4-OH_5Ac*>N3-OH_5Ac*>N5-OH_5Ac* both in the gaseous and solution phase (using the PCM model) respectively due to the additions of *OH at C6, C4, N3, and N5 atoms. These additions occur in direct manner, without the intervention of any precursor complex formation. The possibility of a 1,2-hydrogen shift from the C6 to N5 in the nitrogen-centered C6-OH_5Ac* radical is considered in order to account for the experimental observation of the high yield of non-oxidizing radicals, and found that such a conversion requires activation energy of about 32 kcal/mol, and hence this possibility is ruled out. The hydrogen abstraction reactions were assumed to occur from precursor complexes (hydrogen bonded complexes represented as S1, S2, S3, and S4) resulted from the electrostatic interactions of the lone pairs on the N3, N5, and O8 atoms with the incoming *OH radical. It was found that the conversion of these precursor complexes to their respective transition states has ample barrier heights, and it persists even when the effect of solvent is considered. It was also found that the formation of precursor complexes itself is highly endergonic in solution phase. Hence, the abstraction reactions will not occur in the present case. Finally, the time dependent density functional theory (TDDFT) calculations predicted an absorption maximum of 292 nm for the N3-OH_5Ac* adduct, which is close to the experimentally observed spectral maxima at 290 nm. Hence, it is assumed that the addition to the most reactive center N3, which results the N3-OH_5Ac* radical, occurs via a kinetically driven process.     

Antioxidant activity of thiamine and its structural analogs in reactions with electrochemically generated hydroxyl radicals and hydrogen peroxide

We have used differential pulse voltammetry on platinum and copper cathodes to study the antioxidant activity of thiamine (vitamin B₁) and its structural analogs. The results indicate that the activity of thiamine relative to the hydroxyl radical is lower than the activity of 3-benzyl-4-methyl-5-(2-hydroxyethyl)thiazolium chloride. At the same time, thiamine can more efficiently react with hydrogen peroxide and bind divalent iron ions.     

Thermodynamics and kinetics of glyoxal dimer formation: a computational study

Density functional theory (B3LYP//6-311+G) calculations including Poisson-Boltzmann implicit solvent were used to study the hydration of glyoxal and subsequent formation of dimeric species in solution. Our calculations show that the dioxolane ring dimer is the thermodynamic sink among all monomers and dimers with varying degrees of hydration. Although fully hydrated species are thermodynamically favored over their less hydrated counterparts, we find that a preliminary dehydration step precedes dimerization and ring closure. Ring closure of the open dimer monohydrate to the dioxolane ring dimer is kinetically favored over both hydration to the open dimer dihydrate and ring closure to form the dioxane ring dimer. The kinetic barriers for different geometric approaches for dimerization suggest an explanation why oligomerization stops after the formation of a dioxolane ring trimer as observed experimentally.     

On the electrophilicity of hydroxyl radical: a laser flash photolysis and computational study

The rate coefficients for reactions of hydroxyl radical with aromatic hydrocarbons were measured in acetonitrile using a novel laser flash photolysis method. Comparison of kinetic data obtained in acetonitrile with those obtained in aqueous solution demonstrates an unexpected solvent effect on the reactivity of hydroxyl radical. In particular, reactions of hydroxyl radical with benzene were faster in water than in acetonitrile, and by a significant factor of 65. Computational studies, at the B3LYP and CBS-QB3 levels, have confirmed the rate enhancement of hydroxyl radical addition to benzene via calculation of the transition states in the presence of explicit solvent molecules as well as a continuum dielectric field. The origin of the rate enhancement lies entirely in the structures of the transition states and not in the pre-reactive complexes. The calculations reveal that the hydroxyl radical moiety becomes more anionic in the transition state and, therefore, looks more like hydroxide anion. In the transition states, solvation of the incipient hydroxide anion is more effective with water than with acetonitrile and provides the strong energetic advantage for a polar solvent capable of hydrogen bonding. At the same time, the aromatic unit looks more like the radical cation in the transition state. The commonly held view that hydroxyl radical is electrophilic in its reactions with DNA bases is, therefore, strongly dependent on the ability of the organic substrate to stabilize the resulting radical cation.     

Reaction of phenyl radicals with acetylene: quantum chemical investigation of the mechanism and master equation analysis of the kinetics

The mechanism of the C(6)H(5) + C(2)H(2) reaction has been investigated by various quantum chemical methods. Electrophilic addition to the CC triple bond is found to be the only important mode of phenyl radical attack on acetylene. The initially formed chemically activated C(6)H(5)C(2)H(2) adducts may follow several isomerization pathways in competition with collisional stabilization and H-elimination. Thermochemistry of various decomposition and isomerization channels is evaluated by the G2M method. For key intermediates, the following standard enthalpies of formation have been deduced from isodesmic reactions: 94.2 +/- 2.0 kcal/mol (C(6)H(5)CHCH), 86.4 +/- 2.0 kcal/mol (C(6)H(5)CCH(2)), and 95.5 +/- 1.8 kcal/ mol (o-C(6)H(4)C(2)H(3)). The accuracy of theoretical predictions was examined through extensive comparisons with available experimental and theoretical data. The kinetics and product branching of the C(6)H(5) + C(2)H(2) reaction have been evaluated by weak collision master equation/Rice-Ramsperger-Kassel-Marcus (RRKM) analysis of the truncated kinetic model including only kinetically important transformations of the isomeric C(8)H(7) radicals. Available experimental kinetic data can be quantitatively reproduced by calculation with a minor adjustment of the C(6)H(5) addition barrier from 3.7 to 4.1 kcal/mol. Our predicted total rate constant, k(R1) = (1.29 x 10(10))T(0.834) exp(-2320/T) cm(3) mol(-)(1) s(-)(1), is weakly dependent on P and corresponds to the phenylation process under combustion conditions (T > 1000 K).     

Fast kinetics of the reactions of hydroxyl radicals with nitrone spin traps

The technique of spin trapping with nitrone spin traps nas gained wide acceptance as a method for estimating·OH yields in ESR studies. In our study, fast optical kinetic techniques applied to a series of these traps (PBN, 2-PyBN, 3-PyBN, 4-PyBN, 3-PyOBN and 4-PyOBN) reveal relaxation spectra that indicate two absorption maxima with different time constants, with all except 4-PyOBN showing second order behavior. These two spectral regions show different kinetics. Thus, two reaction sites are indicated, only one of which is necessarily a measure of initial · OH when ESR methods are used. One other trap (DMPO) after · OH reaction decays in one mode suggesting that its final product might be useful as a measure of initial · OH. Also, our ESR evidence shows that OH detection can be improved significantly by spin trapping -hydroxyalkyl radicals formed by · OH attack on alcohols.     

Sensitivity analysis in chemical kinetics: Recent developments and computational comparisons

The advantages and disadvantages of various methods of parametric sensitivity analysis in chemical kinetic modeling are discussed. Particular attention is given to estimates of computational labor for realistic problems, and quantitative comparisons are made utilizing a 52-reaction, 11-species CO oxidation mechanism. The authors′ CHEMSEN/AIM program compares favorably to other techniques in many circumstances, and provides the additional convenience of accepting input information in familiar chemical notation. This paper also reviews recent developments in theory of sensitivity analysis, relevant to chemical kinetic modeling.