Geometry of the transition state of radical abstraction reactions involving Si-H, Ge-H, and Sn-H Bonds

Transition-state interatomic distances in the reactions C˙H₃+SiH₄, Si˙H₃+SiH₄, C˙H₃+GeH₄, and C˙H₃+SnH₄ are calculated by the B3LYP density functional and intersecting parabolas methods. A semiempirical algorithm is developed for the calculation of the Y...H and C...H distances in the transition state of the radical abstraction reactions R˙+YH involving silanes, germanes, and stannanes and the reverse reactions of silyl, germanyl, and stannyl radicals with hydrocarbons. This algorithm is used to calculate interatomic distances in these reactions. An analysis of the calculated data shows that the Y...H and C...H distances in these reactions depend on the following factors: the enthalpy of reaction, the radius of the Y atom (Y = C, Si, Ge, Sn), and four-electron repulsion during the attack of a radical on the C-H bond adjacent to the double bond. Empirical equations relating the interatomic distances to the enthalpy of reaction and to the Y-R bond length are set up.     

Kinetic parameters and geometry of the transition state of acyl radical decarbonylation

Experimental data on acyl radical decomposition reactions (RC^·O → R^· + CO, where R = alkyl or aryl) are analyzed in terms of the intersecting parabolas method. Kinetic parameters characterizing these reactions are calculated. The transition state of methyl radical addition to CO at the C atoms is calculated using the DFT method. A semiempirical algorithm is constructed for calculating the transition state geometry for the decomposition of acyl radicals and for the reverse reactions of R^· addition to CO. Kinetic parameters (activation energy and rate constant) and geometry (interatomic distances in the transition state) are calculated for 18 decomposition reactions of structurally different acyl radicals. A linear correlation between the interatomic distance r ^#(C…C) (or r ^#(C…O)) in the transition state the enthalpy of the reaction (δH _e) is established for acyl decomposition reactions (at br _e = const). A comparative analysis of the enthalpies, activation energies, and interatomic distances in the transition state is carried out for the decomposition and formation of acyl, carboxyl, and formyl radicals.     

Geometric parameters of transition states of radical abstraction reactions with C...H...C reaction center

An algorithm of calculations of interatomic distances in the transition states (TS) of reactions of hydrogen abstraction by alkyl, allyl, and benzyl radicals from C—H bonds of organic molecules using the enthalpies of the corresponding reactions is proposed. The geometric parameters of the TS of the reactions involving carbon-centered radicals with the C...H...C reaction center, calculated using experimental data, are compared with other characteristics of the reactions and reactants. The r(C...H...C) distance in the TS of the reactions of alkyl radicals with alkanes remains unchanged as the enthalpies of reactions vary, being a characteristic parameter of a given class of reactions. -Bonds adjacent to the reaction center are responsible for an increase in the parameter r(C...H...C) in the TS.     

Radical addition reactions: factors determining the transition-state geometry

Interatomic distances in the reaction centers of the addition reactions of (i) H^· to the C=C, C=O, N≡C, and C≡C bonds, (ii) ^·CH₃ radical to the C=C, C=O, and C≡C bonds, and (iii) alkyl, aminyl, and alkoxyl radicals to olefin C=C bonds were determined using a new semiempirical method for calculating transition-state geometries of radical reactions. For all reactions of the type X^· + Y=Z → X— Y—Z^· the r ^# _X...Y distance in the transition state is a linear function of the enthalpy of reaction. Parameters of this dependence were determined for seventeen classes of radical addition reactions. The bond elongation, Δr ^# _X...Y, in the transition state decreases as the triplet repulsion, electronegativity difference between the atoms X and Y in the reaction center, and the force constant of the attacked multiple bond increase. __________ Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 894–902, April, 2005.     

The energy and geometric characteristics of the transition state in reactions of RO<Stack><Subscript>2</Subscript><Superscript>•</Superscript></Stack> with carbonyl compound C-H bonds

The energy and geometry of the transition state in reactions of the ethyl peroxyl radical with ethane, ethanol (its α and β C-H bonds), acetone, butanone-2, and acetaldehyde were calculated by the density functional theory method. In all these reactions (except EtO_2/? + ethanol α C-H bond), the C…H…O reaction center has an almost linear configuration (φ = 176° ± 2°); polar interaction only influences the r ^≠ (C…O) interatomic bond. In the reaction of EtO_2/? with the ethanol α C-H bond, it is the O-H…O H-bond formed in the transition state that determines the configuration of the reaction center with the angle φ(C…H…O) = 160°. The results were used to estimate the r ^≠ (C…H) and r ^≠ (O…H) interatomic bonds in the transition state by the method of intersecting parabolas and the contribution of polar interaction to the activation energy of reactions between peroxyl radicals and aldehydes and ketones.     

Variational analysis of the phenyl + O2 and phenoxy + O reactions

Variational transition state analysis was performed on the barrierless phenyl + O2 and phenoxy + O association reactions. In addition, we also calculated rate constants for the related vinyl radical (C2H3) + O2 and vinoxy radical (C2H3O) + O reactions and provided rate constant estimates for analogous reactions in substituted aromatic systems. Potential energy scans along the dissociating C-OO and CO-O bonds (with consideration of C-OO internal rotation) were obtained at the O3LYP/6-31G(d) density functional theory level. The CO-O and C-OO bond scission reactions were observed to be barrierless, in both phenyl and vinyl systems. Potential energy wells were scaled by G3B3 reaction enthalpies to obtain accurate activation enthalpies. Frequency calculations were performed for all reactants and products and at points along the potential energy surfaces, allowing us to evaluate thermochemical properties as a function of temperature according to the principles of statistical mechanics and the rigid rotor harmonic oscillator (RRHO) approximation. The low-frequency vibrational modes corresponding to R-OO internal rotation were omitted from the RRHO analysis and replaced with a hindered internal rotor analysis using O3LYP/6-31G(d) rotor potentials. Rate constants were calculated as a function of temperature (300-2000 K) and position from activation entropies and enthalpies, according to canonical transition state theory; these rate constants were minimized with respect to position to obtain variational rate constants as a function of temperature. For the phenyl + O2 reaction, we identified the transition state to be located at a C-OO bond length of between 2.56 and 2.16 A (300-2000 K), while for the phenoxy + O reaction, the transition state was located at a CO-O bond length of 2.00-1.90 A. Variational rate constants were fit to a three-parameter form of the Arrhenius equation, and for the phenyl + O2 association reaction, we found k(T) = 1.860 x 1013T-0.217 exp(0.358/T) (with k in cm3 mol-1 s-1 and T in K); this rate equation provides good agreement with low-temperature experimental measurements of the phenyl + O2 rate constant. Preliminary results were presented for a correlation between activation energy (or reaction enthalpy) and pre-exponential factor for heterolytic O-O bond scission reactions.     

Model of the radical addition reaction as the superposition of three potential curves

A new semiempirical model of the reaction of radical addition to molecules with multiple bonds has been developed. In the framework of this model, the transition state (TS) of the reaction X^. + Y=Z → XYZ^. is considered as the result of the intersection of the potential curve of the stretching vibration of the forming bond X-Y with the curve that is the difference between the amplitudes of stretching vibrations of the Y-Z and Y=Z bonds and the stretching vibrations are considered harmonic. The kinetic parameters describing the activation energy as a function of the enthalpy of the reaction were calculated for 34 classes of addition reactions using the new model. The factors determining the activation energy of the addition reactions are analyzed: triplet repulsion in the TS, the π electrons in the α position to the reaction center, the electronegativity of atoms of the reaction center of the TS, the steric factor, the interaction of polar groups in the TS, and the force constants of the reacting bonds. The increments characterizing the contribution of these factors to the activation energy are calculated. The model is also used to describe the energy of 12 classes of cyclization reactions and 16 classes of radical decomposition reactions. The parameters that make it possible to estimate the activation energy of the reaction from its enthalpy are calculated for these classes of reactions.     

Reactivity of C-H and O-H bonds in reactions with aminyl and nitroxyl radicals and cyclic mechanisms of chain termination in oxidizable alcohols and olefins

The parabolic model of transition state has been used to analyze the problem of why aromatic amines and nitroxyl radicals cause the cyclic mechanism of chain termination in oxidizable alcohols and olefins (where HO₂· and >C(OH)0₂· radicals participate in chain propagation) and not in oxidizable hydrocarbons. The differences are caused by the existence of a weak triplet repulsion in transition states with the N...H...0 and 0...H...0 reaction centers, while the triplet repulsion is strong in transition states with reaction centers of the C...H...0 and C...H...N type in oxidizable hydrocarbons.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1972–1976, August, 1996.     

Reactions of Peroxy Radicals with Hydrocarbons and Alcohols: Transition State Geometry and the Polar Effect

The reaction EtOO + EtH EtOOH + Et was studied by the intersecting parabolas method and calculated using density functional theory. The interatomic C–H, O–H, and C–O distances of the transition state for this reaction were calculated using these methods. The formulas for calculating these distances from experimental data were obtained. Similar calculations and comparisons were carried out for the reaction EtOO + MeCH₂OH EtOOH + MeCHOH. The polar effect of the hydroxy groups on the transition state manifested itself in a decrease in the activation energy and in the formation of a nonlinear structure of the transition state. An empirical formula for estimating the C–H–O angle in the transition state from the enthalpy and activation energy was derived.     

Reaction of the hydroxyl radical with phenol in water up to supercritical conditions

The rate constants for the reactions of phenol with the hydroxyl radical (OH*) in water have been measured from room temperature to 380 degrees C using electron pulse radiolysis and transient absorption spectroscopy. The reaction scheme designed to fit the data shows the importance of an equilibrium, giving back reactants (OH* radical and phenol) from the dihydroxycyclohexadienyl radical formed by their reaction, and the non-negligible contribution of the hydroxycyclohexadienyl radical absorption from H* atom addition. The accuracy of the reaction scheme and the reaction rate constants determined from it have been determined by the analysis of two different experiments, one under pure N2O atmosphere and the second under a mixture a N2O and O2. We report reaction rates for the H* and OH* radical addition to phenol, the formation of phenoxyl, the second-order recombination, the reaction of dihydroxycyclohexadienyl with O2, and the decay of the peroxyl adduct. Nearly all of the reaction rates deviate strongly from Arrhenius behavior.     

CH2ClO与NO反应机理的理论研究

采用B3LYP, MP2方法在6-31＋(d,p)和6-311＋＋G(d,p)水平研究了CH₂ClO自由基与NO反应的微观机理, 找到了三个可能的反应通道. 并得到了各反应通道的反应物、中间体、过渡态和产物的优化构型、谐振频率. 成功地解释了Wu等的实验结论. 从电子密度拓扑分析的角度, 讨论了化学反应过程中化学键的变化规律, 为实验研究大气化学反应提供理论依据. 找到了该反应的结构过渡态(结构过渡区)和能量过渡态, 发现了反应热与结构过渡区之间的关系.     

Reaction dynamics simulations of the identity S(N)2 reaction H(2)O + HOOH(2)(+)--> H(2)OOH(+)+ H(2)O. Requirements for reaction and competition with proton transfer

Despite the fact that the transition structure of the gas phase S(N)2 reaction H(2)O + HOOH(2)(+)--> HOOH(2)(+)+ H(2)O is well below the reactants in potential energy, the reaction has not yet been observed by experiment. Variational transition state RRKM theory reveals a strong preference for the competing proton transfer reaction H(2)O + HOOH(2)(+)--> H(3)O(+)+ HOOH due to entropy factors. Born-Oppenheimer reaction dynamics simulations confirm these results. However, by increasing the collision energy to around 7.5 eV the probability for nucleophilic substitution increases relative to proton transfer. These observations are explained by the presence of the key common intermediate HOO(H)[dot dot dot]H-OH(2)(+) which leads to effective proton transfer, but can be avoided with increasing collision energy. However, the S(N)2 probability remains below 0.2 since successful passage through the TS requires optimum initial orientation of the reactants, excitation of the relative translational motion and good phase correlation between the O-O vibration and the motion of the incoming water.     

Reaction rates for O + HD → OH + D and O + HD → OD + H

We have calculated reaction rates for the reactions O + HD → OH + D and O + DH → OD + H using improved canonical variational transition state theory and least-action ground-state transmission coefficients with an ab initio potential energy surface. The kinetic isotope effects are in good agreement with experiment. The optimized tunneling paths and properties of the variational transition states and the rate enhancement for vibrationally excited reactants are also presented and compared with those for the isotopically unsubstituted reaction O + H₂ → OH + H. The thermal reactions at low and room temperature are predicted to occur by tunneling at extended configurations, i.e., to initiate early on the reaction path and to avoid the saddle point regions. Tunneling also dominates the low and room temperature reactions for excited vibrational states, but in these cases the results are not as sensitive to the nature of the tunneling path. Overbarrier mechanisms dominate for both thermal and excited-vibrational state reactions for T > 600 K. For the excited-state reaction (with initial vibrational quantum number n > 0) a transition state switch occurs for T > 1000 K for the O + HD(n = 1) → OD + H case and for T > 1500 K for the O + DH(n = 1) → OD + H reaction, and this may be a general phenomenon for excited-state reactions at higher temperature. In the present case the switch occurs from an early variational transition state where the vibrationally adiabatic approximation is expected to be valid to a tighter variational transition state where nonadiabatic effects are probably important and should be included.     

CH2NH与O(3P)反应的量子化学及电子密度拓扑研究

采用B3LYP、MP2(full)和 QCISD 三种方法在6-311G(d, p)和aug-cc-pVDZ基组水平上对三线态O(3P)原子与CH2NH(s)的反应进行了详细的理论研究. 采用B3LYP和MP2(full)方法对反应势能面上的各驻点进行了几何构型优化, 通过振动频率分析证实了过渡态的真实性, 内禀反应坐标(IRC)跟踪验证了过渡态与反应物和产物的连接关系, 用上述三种方法计算得到了各反应通道的反应势垒. 对反应过程中的一些重要点进行了电子密度拓扑分析研究. 研究结果表明, O(3P)原子进攻CH2NH(s)中的N2原子和C1原子生成CH2NHO(t)和OCH2NH(t), CH2NHO(t)中N2上的H5可迁移到C1上异构化为CH3NO(t).     

Dye-sensitized photooxygenation of the C=N bond. 5. substituent effects on the cleavage of the C=N bond of C-aryl-N-aryl-N-methylhydrazones

The title compounds are cleaved cleanly at the C=N bond by singlet oxygen ((1)O(2), (1)Delta(g)) yielding arylaldehydes and N-aryl-N-methylnitrosamines. These reactions take place more rapidly at -78 degrees C than at room temperature. The effects of substituent variation at both the C-aryl and N-aryl groups were studied using a competitive method. Good correlations of the resulting rate ratios with substituent constants (sigma(-) or sigma(+)) were obtained yielding small to very small rho values indicative of small to very small changes in charge distribution between the reactant and the rate determining transition state. Electron withdrawing groups on the C-aryl moiety retard reaction somewhat by preferential stabilization of the hydrazone. Electron donors on the other hand slightly stabilize the rate determining transition state. Substituents on the N-aryl group have almost no effect. Inhibition by 3,5-di-tert-butylphenol was not observed showing that free (uncaged) radical intermediates are not involved in the mechanism. We postulate a mechanism in which the initial event is exothermic electron transfer from the hydrazone to (1)O(2) leading to an ion-radical caged pair. Subsequent covalent bond formation between the hydrazone carbon and an oxygen atom is rate controlling. The transition state for this step also has a lower enthalpy than the starting reactants, but the free energy of activation is dominated by a large negative TDeltaS++term leading to the negative temperature dependence. Direct formation of a C-O bond in the initial step is not unambiguously ruled out. Subsequent steps lead to C-N cleavage.     

Polar effect and geometry of the transition state in the reactions of ozone with C-H bonds of polar molecules

A large body of experimental data on the reactions of ozone with C-H bonds of polar molecules in the liquid and gas phases is analyzed in the framework of the intersecting parabolas model. The reactions are considered as the abstraction reaction O₃ + RH → HOOO^. + R^.. The contribution from the polar effect to the activation energy of such reactions is calculated. This contribution is ?6.8 kJ/mol for the reactions of ozone with aliphatic alcohols, and is ?8.1, ?11.7, ?6.8, and ?2.2 kJ/mol for the reactions of ozone with ketones, ethers, 1,3-dioxolanes, and 1,3-dioxanes, respectively. The contribution is insignificant in the reactions of ozone with aldehydes. The interatomic distances in the transition state of these reactions r ^#(C…H) and r ^#(O…H) and the angle between the C…H and O…H bonds are calculated. For the reactions in polar solvents, the contribution from solvation to the activation energy is calculated. In most of the systems considered, this contribution is insignificant (from ?1 to ?3 kJ/mol). The reactions involving ozone are compared to the reactions of peroxy radicals with the same classes of compounds.     

煤炭燃烧过程中N2O消除反应机理的密度泛函理论研究

用密度泛函理论B3LYP方法对煤炭燃烧过程中N₂O的消除反应进行研究。选用6-311++G**和aug-cc-pVTZ基组,优化了反应通道上反应物、过渡态和产物的几何构型。预测了它们的热力学性质(总能量、焓、熵和吉布斯自由能)及其随温度的变化。预测N₂O+CO反应的活化能为200 kJ·mol^-1,与实验值193±8 kJ·mol^-1较一致。计算了500~1 800 K 温度范围的反应速率常数。在N₂O的分解中,N₂O与H和CN自由基的反应为动力学优先进行的反应,其活化能为50~55 kJ·mol^-1。在B3LYP/aug-cc-pVTZ level水平下,N₂O+CN反应是热力学最有利的自发反应,其吉布斯自由能变化为-407 kJ·mol^-1。     

The reactivity of phenoxyl radicals of bioantioxidants in the abstraction reactions

The enthalpies, activation energies, and rate constants for the reactions of 15 phenoxyl radicals derived from natural bioantioxidants with hydroperoxides, C-H bonds of linoleic acid, SH-groups of l-cysteine, and O-H bonds of α-tocopherol (60 reactions) were calculated. The activation energies were calculated using the model of intersecting parabolas. The interatomic distances in the reaction sites of the transition states of the studied reactions were calculated. The factors affecting the reactivity of these radicals are discussed. The activation energy of the reaction of oxygen with the O-H bond of the 1,2-dihydroxybenezene semiquinone radical was estimated.     

Thermodynamic origin of the increased rate of hydrolysis of phosphate and phosphorothioate esters in DMSO/water mixtures

The hydrolysis rates of the dianions of phosphate and phosphorothioate esters are substantially accelerated by the addition of polar aprotic solvents such as DMSO and acetonitrile. The activation barrier DeltaG is smaller due to a lower enthalpy of activation. The enthalpy of transfer of p-nitrophenyl phosphate (pNPP) and p-nitrophenyl phosphorothioate (pNPPT), from water to 0.6 (mol) aq DMSO (60 mol % water in DMSO) were measured calorimetrically. The enthalpies of activation for the hydrolysis reactions in the two solvents permitted the calculation of the enthalpy of transfer of the transition states. This transfer is thermodynamically favorable for both the reactants and the transition states but is more favorable for the transition states. In the case of pNPP, the enthalpy of transfer of the reactant is -23.9 kcal/mol, compared to -28.3 for the transition state. The difference is greater for pNPPT, where the enthalpy of transfer of the reactant is -23.2 kcal/mol and that for the transition state is -35.3. The results show that the reduced enthalpies of activation in both hydrolysis reactions arise not from a destabilization of the reactants in the mixed solvent, but from the fact that the enthalpy of transfer of the transition states to the mixed solvent is significantly more negative than the enthalpy of transfer of the reactants.     

Thermochemical and kinetic analyses on oxidation of isobutenyl radical and 2-hydroperoxymethyl-2-propenyl radical

In recognition of the importance of the isobutene oxidation reaction in the preignition chemistry associated with engine knock, the thermochemistry, chemical reaction pathways, and reaction kinetics of the isobutenyl radical oxidation at low to intermediate temperature range were computationally studied, focusing on both the first and the second O2 addition to the isobutenyl radical. The geometries of reactants, important intermediates, transition states, and products in the isobutenyl radical oxidation system were optimized at the B3LYP/6-311G(d,p) and MP2(full)/6-31G(d) levels, and the thermochemical properties were determined on the basis of ab initio, density functional theory, and statistical mechanics. Enthalpies of formation for several important intermediates were calculated using isodesmic reactions at the DFT and the CBS-QB3 levels. The kinetic analysis of the first O2 addition to the isobutenyl radical was performed using enthalpies at the CBS-QB3 and G3(MP2) levels. The reaction forms a chemically activated isobutenyl peroxy adduct which can be stabilized, dissociate back to reactants, cyclize to cyclic peroxide-alkyl radicals, and isomerize to the 2-hydroperoxymethyl-2-propenyl radical that further undergoes another O2 addition. The reaction channels for isomerization and cyclization and further dissociation on this second O2 addition were analyzed using enthalpies at the DFT level with energy corrections based on similar reaction channels for the first O2 addition. The high-pressure limit rate constants for each reaction channel were determined as functions of temperature by the canonical transition state theory for further kinetic model development.