High-level direct-dynamics variational transition state theory calculations including multidimensional tunneling of the thermal rate constants, branching ratios, and kinetic isotope effects of the hydrogen abstraction reactions from methanol by atomic hydrogen

We report a detailed theoretical study of the hydrogen abstraction reaction from methanol by atomic hydrogen. The study includes the analysis of thermal rate constants, branching ratios, and kinetic isotope effects. Specifically, we have performed high-level computations at the MC3BB level together with direct dynamics calculations by canonical variational transition state theory (CVT) with the microcanonically optimized multidimensional tunneling (μOMT) transmission coefficient (CVT/μOMT) to study both the CH(3)OH+H→CH(2)OH+H(2) (R1) reaction and the CH(3)OH+H→CH(3)O+H(2) (R2) reaction. The CVT/μOMT calculations show that reaction R1 dominates in the whole range 298≤T (K)≤2500 and that anharmonic effects on the torsional mode about the C-O bond are important, mainly at high temperatures. The activation energy for the total reaction sum of R1 and R2 reactions changes substantially with temperature and, therefore, the use of straight-line Arrhenius plots is not valid. We recommend the use of new expressions for the total R1 + R2 reaction and for the R1 and R2 individual reactions.     

Branching ratios in the hydroxyl reaction with propene

The branching ratios for the reactions of attachment of hydroxyl radical to propene and hydrogen-atom abstraction were measured at 298 K over the buffer gas pressure range 60-400 Torr (N(2)) using a subatmospheric pressure turbulent flow reactor coupled with a chemical ionization quadrupole mass spectrometer. Isotopically enriched water H(2)(18)O was used to produce (18)O-labeled hydroxyl radicals in reaction with fluorine atoms. The β-hydroxypropyl radicals formed in the attachment reactions 1a and 1b , OH + C(3)H(6) → CH(2)(OH)C(?)HCH(3) (eq 1a ) and OH + C(3)H(6) → C(?)H(2)CH(OH)CH(3) (eq 1b ), were converted to formaldehyde and acetaldehyde in a sequence of secondary reactions in O(2)- and NO-containing environment. The (18)O-labeling propagates to the final products, allowing determination of the branching ratio for the attachment channels of reaction 1. The measured branching ratio for attachment is β(1b) = k(1b)/(k(1a) + k(1b)) = 0.51 ± 0.03, independent of pressure over the 60-400 Torr pressure range. An upper limit on the hydrogen-abstraction channel, OH + C(3)H(6) → H(2)O + C(3)H(5) (eq 1c ), was determined by measuring the water yield in reactions of OH and OD radicals (produced via H(D) + NO(2) → OH(OD) + NO reactions) with C(3)H(6) as k(1c)/(k(1a) + k(1b) + k(1c)) < 0.05 (at 298 K, 200 Torr N(2)).     

Hydrogen atom reactivity toward aqueous tert-butyl alcohol

Through a combination of pulse radiolysis, purification, and analysis techniques, the rate constant for the H + (CH(3))(3)COH → H(2) + (?)CH(2)C(CH(3))(2)OH reaction in aqueous solution is definitively determined to be (1.0 ± 0.15) × 10(5) M(-1) s(-1), which is about half of the tabulated number and 10 times lower than the more recently suggested revision. Our value fits on the Polanyi-type, rate-enthalpy linear correlation ln(k/n) = (0.80 ± 0.05)ΔH + (3.2 ± 0.8) that is found for the analogous reactions of other aqueous aliphatic alcohols with n equivalent abstractable H atoms. The existence of such a correlation and its large slope are interpreted as an indication of the mechanistic similarity of the H atom abstraction from α- and β-carbon atoms in alcohols occurring through the late, product-like transition state. tert-Butyl alcohol is commonly contaminated by much more reactive secondary and primary alcohols (2-propanol, 2-butanol, ethanol, and methanol), whose content can be sufficient for nearly quantitative scavenging of the H atoms, skewing the H atom reactivity pattern, and explaining the disparity of the literature data on the H + (CH(3))(3)COH rate constant. The ubiquitous use of tert-butyl alcohol in pulse radiolysis for investigating H atom reactivity and the results of this work suggest that many other previously reported rate constants for the H atom, particularly the smaller ones, may be in jeopardy.     

Theoretical study on the OH + CH3NHCOOCH3 reaction

The multiple-channel reactions OH + CH3NHC(O)OCH3 --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-311+G(d,p) level, and energetic information is further refined by the BMC-CCSD (single-point) method. The rate constants for every reaction channels, R1, R2, R3, and R4, are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-1000 K. The total rate constants are in good agreement with the available experimental data and the two-parameter expression k(T) = 3.95 x 10(-12) exp(15.41/T) cm3 molecule(-1) s(-1) over the temperature range 200-1000 K is given. Our calculations indicate that hydrogen abstraction channels R1 and R2 are the major channels due to the smaller barrier height among four channels considered, and the other two channels to yield CH3NC(O)OCH3 + H2O and CH3NHC(O)(OH)OCH3 + H2O are minor channels over the whole temperature range.     

13C, 18O, and D fractionation effects in the reactions of CH3OH isotopologues with Cl and OH radicals

A relative rate experiment is carried out for six isotopologues of methanol and their reactions with OH and Cl radicals. The reaction rates of CH2DOH, CHD2OH, CD3OH, (13)CH3OH, and CH3(18)OH with Cl and OH radicals are measured by long-path FTIR spectroscopy relative to CH3OH at 298 +/- 2 K and 1013 +/- 10 mbar. The OH source in the reaction chamber is photolysis of ozone to produce O((1)D) in the presence of a large excess of molecular hydrogen: O((1)D) + H2 --> OH + H. Cl is produced by the photolysis of Cl2. The FTIR spectra are fitted using a nonlinear least-squares spectral fitting method with measured high-resolution infrared spectra as references. The relative reaction rates defined as alpha = k(light)/k(heavy) are determined to be: k(OH + CH3OH)/k(OH + (13)CH3OH) = 1.031 +/- 0.020, k(OH + CH3OH)/k(OH + CH3(18)OH) = 1.017 +/- 0.012, k(OH + CH3OH)/k(OH + CH2DOH) = 1.119 +/- 0.045, k(OH + CH3OH)/k(OH + CHD2OH) = 1.326 +/- 0.021 and k(OH + CH3OH)/k(OH + CD3OH) = 2.566 +/- 0.042, k(Cl + CH3OH)/k(Cl + (13)CH3OH) = 1.055 +/- 0.016, k(Cl + CH3OH)/k(Cl + CH3(18)OH) = 1.025 +/- 0.022, k(Cl + CH3OH)/k(Cl + CH2DOH) = 1.162 +/- 0.022 and k(Cl + CH3OH)/k(Cl + CHD2OH) = 1.536 +/- 0.060, and k(Cl + CH3OH)/k(Cl + CD3OH) = 3.011 +/- 0.059. The errors represent 2sigma from the statistical analyses and do not include possible systematic errors. Ground-state potential energy hypersurfaces of the reactions were investigated in quantum chemistry calculations at the CCSD(T) level of theory with an extrapolated basis set. The (2)H, (13)C, and (18)O kinetic isotope effects of the OH and Cl reactions with CH3OH were further investigated using canonical variational transition state theory with small curvature tunneling and compared to experimental measurements as well as to those observed in CH4 and several other substituted methane species.     

Shock tube/laser absorption measurements of the reaction rates of OH with ethylene and propene

Reaction rates of hydroxyl (OH) radicals with ethylene (C?H?) and propene (C?H?) were studied behind reflected shock waves. OH + ethylene → products (rxn 1) rate measurements were conducted in the temperature range 973-1438 K, for pressures from 2 to 10 atm, and for initial concentrations of ethylene of 500, 751, and 1000 ppm. OH + propene → products (rxn 2) rate measurements spanned temperatures of 890-1366 K, pressures near 2.3 atm, and initial propene concentrations near 300 ppm. OH radicals were produced by shock-heating tert-butyl hydroperoxide, (CH?)?-CO-OH, and monitored by laser absorption near 306.7 nm. Rate constants for the reactions of OH with ethylene and propene were extracted by matching modeled and measured OH concentration time-histories in the reflected shock region. Current data are in excellent agreement with previous studies and extend the temperature range of OH + propene data. Transition state theory calculations using recent ab initio results give excellent agreement with our measurements and other data outside our temperature range. Fits (in units of cm3/mol/s) to the abstraction channels of OH + ethylene and OH + propene are k? = 2.23 × 10? (T)(2.745) exp(-1115 K/T) for 600-2000 K and k? = 1.94 × 10? (T)(2.229) exp(-540 K/T) for 700-1500 K, respectively. A rate constant determination for the reaction TBHP → products (rxn 3) was also obtained in the range 745-1014 K using OH data from behind both incident and reflected shock waves. These high-temperature measurements were fit with previous low-temperature data, and the following rate expression (0.6-2.6 atm), applicable over the temperature range 400-1050 K, was obtained: k? (1/s) = 8.13 × 10?12 (T)(7.83) exp(-14598 K/T).     

H atom attack on propene

The reaction of propene (CH(3)CH═CH(2)) with hydrogen atoms has been investigated in a heated single-pulsed shock tube at temperatures between 902 and 1200 K and pressures of 1.5-3.4 bar. Stable products from H atom addition and H abstraction have been identified and quantified by gas chromatography/flame ionization/mass spectrometry. The reaction for the H addition channel involving methyl displacement from propene has been determined relative to methyl displacement from 1,3,5-trimethylbenzene (135TMB), leading to a reaction rate, k(H + propene) → H(2)C═CH(2) + CH(3)) = 4.8 × 10(13) exp(-2081/T) cm(3)/(mol s). The rate constant for the abstraction of the allylic hydrogen atom is determined to be k(H + propene → CH(2)CH═CH(2) + H(2)) = 6.4 × 10(13) exp(-4168/T) cm(3)/(mol s). The reaction of H + propene has also been directly studied relative to the reaction of H + propyne, and the relationship is found to be log[k(H + propyne → acetylene + CH(3))/k(H + propene → ethylene + CH(3))] = (-0.461 ± 0.041)(1000/T) + (0.44 ± 0.04). The results showed that the rate constant for the methyl displacement reaction with propene is a factor of 1.05 ± 0.1 larger than that for propyne near 1000 K. The present results are compared with relevant earlier data on related compounds.     

Implementation of a microcanonical variational transition state theory for direct dynamics calculations of rate constants

Calculation of microcanonical rate constants has been an important field in chemical dy-namic studies for many years because it can be used not only to give good prediction of rate con-stants in microcanonical assembly, but also to calculate rate constants with certain conserved quantum numbers such as the total angular momentum, and in turn, can be easily converted into thermal rate constants[1—3]. The widely used method for calculating microcanonical rate constants of unimolecular reac-tions…     

Absolute rate coefficient of the OH + CH(3)C(O)OH reaction at T = 287-802 K. The two faces of pre-reactive H-bonding

The rate constants for the reaction OH + CH3C(O)OH --> products (1) were determined over the temperature range 287-802 K at 50 and 100 Torr of Ar or N2 bath gas using pulsed laser photolysis generation of OH by CH3C(O)OH photolysis at 193 nm coupled with OH detection by pulsed laser-induced fluorescence. The rate coefficient displays a complex temperature dependence with a sharp minimum at 530 K, indicating the competition between a reaction proceeding through a pre-reactive H-bonded complex to form CH3C(O)O + H2O, expected to prevail at low temperatures, and a direct methyl-H abstraction channel leading to CH2C(O)OH + H2O, which should dominate at high temperatures. The temperature dependence of the rate constant can be described adequately by k1(287-802 K) = 2.9 x 10(-9) exp{-6030 K/T} + 1.50 x 10(-13) exp{515 K/T} cm3 molecule(-1)(s-1), with a value of (8.5 +/- 0.9) x 10-13 cm3 molecule(-1)(s-1) at 298 K. The steep increase in rate constant in the range 550-800 K, which is reported for the first time, implies that direct abstraction of a methyl-H becomes the dominant pathway at temperatures greater than 550 K. However, the data indicates that up to about 800 K direct methyl-H abstraction remains adversely affected by the long-range H-bonding attraction between the approaching OH radical and the carboxyl -C(O)OH functionality.     

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.     

High-temperature rate constants for CH3OH + Kr --> products, OH + CH3OH --> products, OH + (CH3)(2)CO --> CH2COCH3 + H2O, and OH + CH3 --> CH) + H2O

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the dissociation of methanol between 1591 and 2865 K. Rate constants for two product channels [CH3OH + Kr --> CH3 + OH + Kr (1) and CH3OH + Kr --> 1CH2 + H2O + Kr (2)] were determined. During the course of the study, it was necessary to determine several other rate constants that contributed to the profile fits. These include OH + CH3OH --> products, OH + (CH3)2CO --> CH2COCH3 + H2O, and OH + CH3 --> 1,3CH2 + H2O. The derived expressions, in units of cm(3) molecule(-1) s(-1), are k(1) = 9.33 x 10(-9) exp(-30857 K/T) for 1591-2287 K, k(2) = 3.27 x 10(-10) exp(-25946 K/T) for 1734-2287 K, kOH+CH3OH = 2.96 x 10-16T1.4434 exp(-57 K/T) for 210-1710 K, k(OH+(CH3)(2)CO) = (7.3 +/- 0.7) x 10(-12) for 1178-1299 K and k(OH+CH3) = (1.3 +/- 0.2) x 10(-11) for 1000-1200 K. With these values along with other well-established rate constants, a mechanism was used to obtain profile fits that agreed with experiment to within <+/-10%. The values obtained for reactions 1 and 2 are compared with earlier determinations and also with new theoretical calculations that are presented in the preceding article in this issue. These new calculations are in good agreement with the present data for both (1) and (2) and also for OH + CH3 --> products.     

Temperature dependence of carbon-13 kinetic isotope effects of importance to global climate change

We report here a theoretical study of the 13C kinetic isotope effect (KIE) and its temperature dependence for the reaction OH + CH4 --> H2O + CH3, the major sink of atmospheric methane in the troposphere. The KIE values at various atmospherically significant temperatures were determined by direct dynamics using variational transition state theory with multidimensional tunneling contributions (VTST/MT). The potential energy surfaces (PESs) were generated by hybrid density functional theory as well as by recently developed doubly hybrid density functional theory methods. Comparisons of our calculated KIEs with experimental data and theoretical values in the literature reveal the critical contributions due to multidimensional tunneling and torsion anharmonicity as well as the critical issue of the choice of internal rotational axis.     

Kinetics of the reaction of methyl radical with hydroxyl radical and methanol decomposition

The CH3 + OH bimolecular reaction and the dissociation of methanol are studied theoretically at conditions relevant to combustion chemistry. Kinetics for the CH3 + OH barrierless association reaction and for the H + CH2OH and H + CH3O product channels are determined in the high-pressure limit using variable reaction coordinate transition state theory and multireference electronic structure calculations to evaluate the fragment interaction energies. The CH3 + OH --> 3CH2 + H2O abstraction reaction and the H2 + HCOH and H2 + H2CO product channels feature localized dynamical bottlenecks and are treated using variational transition state theory and QCISD(T) energies extrapolated to the complete basis set limit. The 1CH2 + H2O product channel has two dynamical regimes, featuring both an inner saddle point and an outer barrierless region, and it is shown that a microcanonical two-state model is necessary to properly describe the association rate for this reaction over a broad temperature range. Experimental channel energies for the methanol system are reevaluated using the Active Thermochemical Tables (ATcT) approach. Pressure dependent, phenomenological rate coefficients for the CH3 + OH bimolecular reaction and for methanol decomposition are determined via master equation simulations. The predicted results agree well with experimental results, including those from a companion high-temperature shock tube determination for the decomposition of methanol.     

Theoretical investigation of the hydrogen abstraction reaction of the OH radical with CH3CHF2 (HFC152-a): a dual level direct density functional theory dynamics study

The hydrogen abstraction reaction of the OH radical with CH(3)CHF(2) (HFC152-a) has been studied theoretically over a wide temperature range, 200-3000 K. Two different reactive sites of the molecule, CH(3) and CHF(2) groups have been investigated precisely, and results confirm that CHF(2) position of the molecule is a highly reactive site. In this study, three recently developed hybrid density functional theories, namely, MPWB1K, MPW1B95, and MPW1K, are used. The MPWB1K/6-31+G(d,p) method gives the best result for kinetic calculations, including barrier heights, reaction path information and geometry of transition state structures and other stationary points. To refine the barrier height of each channel, a single point energy calculation was performed in MPWB1K/MG3S method. The obtained rate constants by dual level direct dynamics with the interpolated single point energy method (VTST-ISPE) using DFT quantum computational methods, are consistent with available experimental data. The canonical variational transition state theory (CVT) with the zero-curvature and also the small-curvature tunneling correction methods is used to calculate the rate constants. Over the temperature range 200-3000 K, the variation effect, tunneling contribution, branching ratio of each channel are calculated. The rate constants and their temperature dependency in the form of a fitted three-parameter Arrhenius expression are k(1)(T) = 2.00 x 10(-19)(T)(2.24) exp(-1273/T), k(2)(T) = 1.95 x 10(-19)(T)(2.46) exp(-2374/T), and k(T) = 3.13 x 10(-19)(T)(2.47) exp(- 1694/T) cm(3) molecule(-1) s(-1). For the H abstraction from the CHF(2) group, a nonclassical reflection effect is detected as a dominant quantum effect.     

Site-specific rate coefficients for reaction of OH with ethanol from 298 to 900 K

The rate coefficients for reactions of OH with ethanol and partially deuterated ethanols have been measured by laser flash photolysis/laser-induced fluorescence over the temperature range 298-523 K and 5-100 Torr of helium bath gas. The rate coefficient, k(1.1), for reaction of OH with C(2)H(5)OH is given by the expression k(1.1) = 1.06 × 10(-22)T(3.58)?exp(1126/T) cm(3) molecule(-1) s(-1), and the values are in good agreement with previous literature. Site-specific rate coefficients were determined from the measured kinetic isotope effects. Over the temperature region 298-523 K abstraction from the hydroxyl site is a minor channel. The reaction is dominated by abstraction of the α hydrogens (92 ± 8)% at 298 K decreasing to (76 ± 9)% with the balance being abstraction at the β position where the errors are 2σ. At higher temperatures decomposition of the CH(2)CH(2)OH product from β abstraction complicates the kinetics. From 575 to 650 K, biexponential decays were observed, allowing estimates to be made for k(1.1) and the fractional production of CH(2)CH(2)OH. Above 650 K, decomposition of the CH(2)CH(2)OH product was fast on the time scale of the measured kinetics and removal of OH corresponds to reaction at the α and OH sites. The kinetics agree (within ±20%) with previous measurements. Evidence suggests that reaction at the OH site is significant at our higher temperatures: 47-53% at 865 K.     

Temperature-dependent branching ratios of deuterated methoxy radicals (CH2DO•) reacting with O2

The methoxy radical is an intermediate in the atmospheric oxidation of methane, and the branching ratio (k(1a)/k(1b)) (CH(2)DO? + O(2) → CHDO + HO(2) (1a) and CH(2)DO? + O(2) → CH(2)O + DO(2) (1b)) strongly influences the HD/H(2) ratio in the atmosphere, which is widely used to investigate the global cycling of molecular hydrogen. By using the FT-IR smog chamber technique, we measured the yields of CH(2)O and CHDO from the reaction at 250-333 K. Kinetic modeling was used to confirm the suppression of secondary chemistry. The resulting branching ratios are well fit by an Arrhenius expression: ln(k(1a)/k(1b)) = (416 ± 152)/T + (0.52 ± 0.53), which agrees with the room-temperature results reported in the only previous study. The present results will be used to test our theoretical understanding of the role of tunneling in the methoxy + O(2) reaction, which is the prototype for the entire class of alkoxy + O(2) reactions.     

The first experimental test of the hypothesis that enzymes have evolved to enhance hydrogen tunneling

The literature hypothesis that "the optimization of enzyme catalysis may entail the evolutionary implementation of chemical strategies that increase the probability of quantum-mechanical tunneling" is experimentally tested herein for the first time. The system employed is the key to being able to provide this first experimental test of the "enhanced hydrogen tunneling" hypothesis, one that requires a comparison of the three criteria diagnostic of tunneling (vide infra) for the same, or nearly the same, reaction with and without the enzyme. Specifically, studied herein are the adenosylcobalamin (AdoCbl, also known as coenzyme B(12))-dependent diol dehydratase model reactions of (i). H(D)(*) atom abstraction from ethylene glycol-d(0) and ethylene glycol-d(4) solvent by 5'-deoxyadenosyl radical (Ado(*)) and (ii.) the same H(*) abstraction reactions by the 8-methoxy-5'-deoxyadenosyl radical (8-MeOAdo(*)). The Ado(*) and 8-MeOAdo(*) radicals are generated by Co-C thermolysis of their respective precursors, AdoCbl and 8-MeOAdoCbl. Deuterium kinetic isotope effects (KIEs) of the H(*)(D(*)) abstraction reactions from ethylene glycol have been measured over a temperature range of 80-120 degrees C: KIE = 12.4 +/- 1.1 at 80 degrees C for Ado(*) and KIE = 12.5 +/- 0.9 at 80 degrees C for 8-MeOAdo(*) (values ca. 2-fold that of the predicted maximum primary times secondary ground-state zero-point energy (GS-ZPE) KIE of 6.4 at 80 degrees C). From the temperature dependence of the KIEs, zero-point activation energy differences ([E(D) - E(H)]) of 3.0 +/- 0.3 kcal mol(-)(1) for Ado(*) and 2.1 +/- 0.6 kcal mol(-)(1) for 8-MeOAdo(*) have been obtained, both of which are significantly larger than the nontunneling, zero-point energy only maximum of 1.2 kcal mol(-)(1). Pre-exponential factor ratios (A(H)/A(D)) of 0.16 +/- 0.07 for Ado(*) and 0.5 +/- 0.4 for 8-MeOAdo(*) are observed, both of which are significantly less than the 0.7 minimum for nontunneling behavior. The data provide strong evidence for the expected quantum mechanical tunneling in the Ado(*) and 8-MeOAdo(*)-mediated H(*) abstraction reactions from ethylene glycol. More importantly, a comparison of these enzyme-free tunneling data to the same KIE, (E(D) - E(H)) and A(H)/A(D) data for a closely related, Ado(*)-mediated H(*) abstraction reaction from a primary CH(3)- group in AdoCbl-dependent methylmalonyl-CoA mutase shows the enzymic and enzyme-free data sets are identical within experimental error. The Occam's Razor conclusion is that at least this adenosylcobalamin-dependent enzyme has not evolved to enhance quantum mechanical tunneling, at least within the present error bars. Instead, this B(12)-dependent enzyme simply exploits the identical level of quantum mechanical tunneling that is available in the enzyme-free, solution-based H(*) abstraction reaction. The results also require a similar, if not identical, barrier width and height within experimental error for the H(*) abstraction both within, and outside of, the enzyme.     

Hydrogen atom transfer reactions of a ruthenium imidazole complex: hydrogen tunneling and the applicability of the Marcus cross relation

The reaction of Ru(II)(acac)2(py-imH) (Ru(II)imH) with TEMPO(*) (2,2,6,6-tetramethylpiperidine-1-oxyl radical) in MeCN quantitatively gives Ru(III)(acac)2(py-im) (Ru(III)im) and the hydroxylamine TEMPO-H by transfer of H(*) (H(+) + e(-)) (acac = 2,4-pentanedionato, py-imH = 2-(2'-pyridyl)imidazole). Kinetic measurements of this reaction by UV-vis stopped-flow techniques indicate a bimolecular rate constant k(3H) = 1400 +/- 100 M(-1) s(-1) at 298 K. The reaction proceeds via a concerted hydrogen atom transfer (HAT) mechanism, as shown by ruling out the stepwise pathways of initial proton or electron transfer due to their very unfavorable thermochemistry (Delta G(o)). Deuterium transfer from Ru(II)(acac)2(py-imD) (Ru(II)imD) to TEMPO(*) is surprisingly much slower at k(3D) = 60 +/- 7 M(-1) s(-1), with k(3H)/k(3D) = 23 +/- 3 at 298 K. Temperature-dependent measurements of this deuterium kinetic isotope effect (KIE) show a large difference between the apparent activation energies, E(a3D) - E(a3H) = 1.9 +/- 0.8 kcal mol(-1). The large k(3H)/k(3D) and DeltaE(a) values appear to be greater than the semiclassical limits and thus suggest a tunneling mechanism. The self-exchange HAT reaction between Ru(II)imH and Ru(III)im, measured by (1)H NMR line broadening, occurs with k(4H) = (3.2 +/- 0.3) x 10(5) M(-1) s(-1) at 298 K and k(4H)/k(4D) = 1.5 +/- 0.2. Despite the small KIE, tunneling is suggested by the ratio of Arrhenius pre-exponential factors, log(A(4H)/A(4D)) = -0.5 +/- 0.3. These data provide a test of the applicability of the Marcus cross relation for H and D transfers, over a range of temperatures, for a reaction that involves substantial tunneling. The cross relation calculates rate constants for Ru(II)imH(D) + TEMPO(*) that are greater than those observed: k(3H,calc)/k(3H) = 31 +/- 4 and k(3D,calc)/k(3D) = 140 +/- 20 at 298 K. In these rate constants and in the activation parameters, there is a better agreement with the Marcus cross relation for H than for D transfer, despite the greater prevalence of tunneling for H. The cross relation does not explicitly include tunneling, so close agreement should not be expected. In light of these results, the strengths and weaknesses of applying the cross relation to HAT reactions are discussed.     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

Computational study of the reaction of CH2(X3B1) with CH3OH

The reaction of triplet methylene with methanol is a key process in alcohol combustion but surprisingly this reaction has never been studied. The reaction mechanism is investigated by using various high-level ab initio methods, including the complete basis set extrapolation (CBS-QB3 and CBS-APNO), the latest Gaussian-n composite method (G4), and the Weizmann-1 method (W1U). A total of five product channels and six transition states are found. The dominant mechanism is direct hydrogen abstraction, and the major product channel is CH(3) + CH(3)O, involving a weak prereactive complex and a 7.4 kcal/mol barrier. The other hydrogen abstraction channel, CH(3) + CH(2)OH, is less important even though it is more exothermic and involves a similar barrier height. The rate coefficients are predicted in the temperature range 200-3000 K. The tunneling effect and the hindered internal rotational freedoms play a key role in the reaction. Moreover, the reaction shows significant kinetic isotope effect.