Absolute rate calculations for atom abstractions by radicals: energetic,structural and electronic factors

We calculate transition-state energies of atom-transfer reactions from reaction energies, electrophilicity indices, bond lengths, and vibration frequencies of the reactive bonds. Our calculations do not involve adjustable parameters and uncover new patterns of reactivity. The generality of our model is demonstrated comparing the vibrationally adiabatic barriers obtained for 100 hydrogen-atom transfers with the corresponding experimental activation energies, after correction for the heat capacities of reactants and transition state. The rates of half of these reactions are calculated using the Transition-State Theory with the vibrationally adiabatic path of the Intersecting-State Model and the semiclassical correction for tunneling (ISM/scTST). The calculated rates are within an order of magnitude of the experimental ones at room temperature. The temperature dependencies and kinetic isotope effects of selected systems are also in good agreement with the available experimental data. Our model elucidates the roles of the reaction energy, electrophilicity, structural parameters, and tunneling in the reactivity of these systems and can be applied to make quantitative predictions for new systems.     

Absolute rate calculations: atom and proton transfers in hydrogen-bonded systems.

We calculate energy barriers of atom- and proton-transfer reactions in hydrogen-bonded complexes in the gas phase. Our calculations do not involve adjustable parameters and are based on bond-dissociation energies, ionization potentials, electron affinities, bond lengths, and vibration frequencies of the reactive bonds. The calculated barriers are in agreement with experimental data and high-level ab initio calculations. We relate the height of the barrier with the molecular properties of the reactants and complexes. The structure of complexes with strong hydrogen bonds approaches that of the transition state, and substantially reduces the barrier height. We calculate the hydrogen-abstraction rates in H-bonded systems using the transition-state theory with the semiclassical correction for tunneling, and show that they are in excellent agreement with the experimental data. H-bonding leads to an increase in tunneling corrections at room temperature.     

Insights Gained by Modeling the Kinetics of Archetypal Hydrogen Atom Transfer Reactions: NH3 + ·H ⇆ H2N· + H2 and CH4 + ·H ⇆ H3C· + H2

Experimental measurements of the kinetics of the title reactions extend to temperature ranges of 1360 K for the ammonia‐hydrogen reaction and of 1602 K for the methane‐hydrogen reaction. Curved plots of ln(k) versus 1/T are obtained. Many theoretical calculations modeling these reactions routinely use tunneling corrections to match experiment. The steepness and curvatures of the plots are modeled successfully in this work and are shown to be caused solely by changes in the bond dissociation energies of the bonds involved in the reactions without invoking tunneling or any other adjustable parameters. The conclusion that tunneling does not contribute significantly to the rates in the temperature range of the measurements is in stark contrast with those theoretical calculations invoking large tunneling factors in the experimental temperature range. Support for the conclusion is provided by theoretical calculations of harmonic quantum transition state theory implementing instanton theory. There is direct experimental evidence that significant tunneling occurs in some H atom transfers, as with isotopomers of H₂ + ·H and other H transfers at very low temperatures. However, there is no direct experimental evidence of significant tunneling contributions to the rates of the title reactions in the temperature range of the measurements. Insights are gained into what specific forces must be overcome by the enthalpy of activation for reaction to occur.     

Proton mobility in the solid inorganic hydrates of acids and acid salts

The method for the calculation of the proton-transfer frequency (n_t) and its activation energies (E_t) was suggested. Results of the calculations were presented. The experimental data on the activation energy of proton-containing group rotation and protonconductivity values for some hydrates of inorganic acids and acid salts were compared with the calculated ones.     

Direct-dynamics VTST study of the [1,7] hydrogen shift in 7-methylocta-1,3(Z),5(Z)-triene. A model system for the hydrogen transfer reaction in previtamin D3

Direct-dynamics canonical variational transition-state theory calculations with microcanonically optimized multidimensional transmission coefficient (CVT/muOMT) for tunneling were carried out at the MPWB1K/6-31+G(d,p) level to study the [1,7] sigmatropic hydrogen rearrangement in 7-methylocta-1,3(Z),5(Z)-triene. This compound has seven conformers, of which only one leads to products, although all of them have to be included in the theoretical treatment. The calculated CVT/muOMT rate constants are in good agreement with the available experimental data. To try to understand the role of tunneling in the hydrogen shift reaction, we have also calculated the thermal rate constants for the monodeuterated compound in the interval T = 333.2-388.2 K. This allowed us to evaluate primary kinetic isotope effects (KIEs) and make a direct comparison with the experiment. Our calculations show that both the large measured KIE and the large measured difference in the activation energies between the deuterated and root compounds are due to the quantum tunneling. The tunneling contribution to the KIE becomes noticeable only when the coupling between the reaction coordinate and the transverse modes is taken into account. Our results confirm previous experimental and theoretical works, which guessed that the obtained kinetic parameters pointed to a reaction with an important contribution due to tunneling. The above conclusion would be essentially valid for the case of the [1,7] hydrogen shift in previtamin D3 because of the similarity to the studied model system.     

Theory of electron transfer reactions in blue-copper proteins

The rates of electron transfer reactions in azurin and plastocyanin are calculated with the Intersecting-State Model and compared with experimental data. The calculated distance, free-energy and temperature dependencies of the intraprotein rates in Ru-modified azurins are in good agreement with the experiment. These calculations do not require the fitting of any parameters. Significant tunneling contributions to the room temperature rate are found in some systems. In some cases the symmetry or the orientation of the donor and acceptor orbitals are not favorable and the ET rates are reduced by factors exceeding 4 orders of magnitude.     

Reaction path dynamics and theoretical rate constants for the SiH3Cl+H→SiH2Cl+H2 reaction by ab initio direct dynamics method

Ab initio direct dynamics method has been used to study the title reaction. Electronic structure information including geometries, gradients and force constants (Hessians) are calculated at the UQCISD/6-311+G^** level. Energies along the minimum energy path are improved by a series of single-point G2//QCISD calculations. The changes of the geometries, vibratioanal frequencies, potential energies and total curvature along the reaction path are discussed. The rate constants in the temperature range 200–3000 K are calculated by canonical variational transition state theory with small-curvature tunneling correction (CVT/SCT) method. The results show that the variational effect is small and in the lower temperature range, the small curvature tunneling effect is important for the reaction.     

Dual-level direct dynamics studies for the reactions of dimethyl ether with hydrogen atom and methyl radical

The dual-level direct dynamics approach is employed to study the dynamics of the CH(3)OCH(3) + H (R1) and CH(3)OCH(3) + CH(3) (R2) reactions. Low-level calculations of the potential energy surface are carried out at the MP2/6-311+G(d,p) level of theory. High-level energetic information is obtained at the QCISD(T) level of theory with the 6-311+G(3df,3pd) basis set. The dynamics calculations are performed using variational transition state theory (VTST) with the interpolated single-point energies (ISPE) method, and small-curvature tunneling (SCT) is included. It is shown that the reaction of CH(3)OCH(3) with H (R1) may proceed much easier and with a lower barrier height than the reaction with CH(3) radical (R2). The calculated rate constants and activation energies are in good agreement with the experimental values. The calculated rate constants are fitted to k(R1) = 1.16 x 10(-19) T(3) exp(-1922/T) and k(R2) = 1.66 x 10(-28) T(5) exp(-3086/T) cm(3) mol(-1) s(-1) over a temperature range 207-2100 K. Furthermore, a small variational effect and large tunneling effect in the lower temperature range are found for the two reactions.     

Ab initio direct dynamics studies on the reaction of H atom with CH3CH2Cl

The multiple channel reaction H + CH(3)CH(2)Cl --> products has been studied by the ab initio direct dynamics method. The potential energy surface information is calculated at the MP2/6-311G(d,p) level of theory. The energies along the minimum energy path are further improved by single-point energy calculations at the PMP4(SDTQ)/6-311+G(3df,2p) level of theory. For the reaction, four reaction channels (one chlorine abstraction, one alpha-hydrogen abstraction, and two beta-hydrogen abstractions) have been identified. The rate constants for each reaction channel are calculated by using canonical variational transition state theory incorporating the small-curvature tunneling correction in the temperature range 298-5000 K. The total rate constants, which are calculated from the sum of the individual rate constants, are in good agreement with the experimental data. The calculated temperature dependence of the branching fractions indicates that for the title reaction, H-abstraction reaction is the major reaction channel in the whole temperature range 298-5000 K.     

Density functional theory study of beta-hydride elimination of ethyl on flat and stepped Cu surfaces

Plane wave density functional theory calculations have been used to characterize the transition states for beta-hydride elimination of ethyl on Cu(100), Cu(110), Cu(111), and Cu(221). The reaction rates predicted by these calculations have been compared to experiments by including tunneling corrections within harmonic transition state theory. Tunneling corrections are found to be important in describing the peak temperatures observed using temperature programmed desorption experiments on Cu(110), Cu(111), and Cu(221). Once these corrections are included, the effective activation energies obtained from our calculations are in good agreement with previous experimental studies of this reaction on these four Cu surfaces. The transition states determined in our calculations are used to examine two general hypotheses that have been suggested to describe structure sensitivity in metal-catalyzed surface reactions.     

“Spectroscopic” Calculations of CH Bond Dissociation Energies for Ethane, Propane, Butane, Isobutane, Pentane, Hexane, and Neopentane Using Fundamental Vibration Frequencies

The fundamental spectrum and the parameters of the potential function of a number of saturated hydrocarbon molecules are calculated in an anharmonic approximation. The calculation is performed by the variational technique using a minimal Morse-harmonic basis. The potential function is taken as the sum of the Morse function for CH bonds and the harmonic function for the skeletal and deformation vibrations. The initial approximation for the potential function is found by ab initio calculations in a 6-31G basis and refined by solving the inverse problem. The calculated CH bond dissociation energies depend significantly on the molecular structure and on the position of CH bonds in the molecule. These energies correlate well with the experimental cleavage energies of these bonds. The changes in the dipole moment of the molecule induced by vibrations were found by ab initio calculations in a 6-31G basis. The calculated IR transmission curves are in good agreement with the experimental curves.     

Mechanisms of amine-catalyzed organosilicate hydrolysis at circum-neutral pH

Mono- and polyamines can catalyze the hydrolysis and condensation of organosilicate starting materials in biomimetic silica synthesis pathways at circum-neutral pHs and room temperature. Our study is focused on understanding the mechanistic role of amines in catalyzing the hydrolysis process that precedes condensation. We have conducted (29)Si NMR experimental studies over a range of temperature and pHs for the hydrolysis rates of trimethylethoxysilane (TMES), a model compound with only one hydrolyzable bond, combined with quantum mechanical hybrid density functional theory calculations of putative intermediate and transition-state structures for TMES and tetramethyl orthosilicate (TMOS). Comparison of calculated energies with experimentally determined activation energies indicates that amine catalysis of TMES is primarily a consequence of the amine's acidity at neutral pH. The proton released by the amine is transferred to the organosilicate, producing a protonated ethoxy leaving group that can be displaced by water in an S(N)2 reaction. For TMOS, the activation energy of proton-transfer coupled with S(N)2 substitution is comparable to that for Corriu's nucleophile-activated nucleophilic displacement, such that the mechanism of amine-catalyzed hydrolysis is dependent mostly on the ambient pH conditions as well as the type of amine. The relevance of our results to biological silica precipitation is discussed.     

星际媒介H2NCH2CN与H2O反应中的H2O分子偶合质子转移机理和氢隧道效应

H2NCH2CN+H2O→H2NCH2C(OH)NH是一个重要的反应, 涉及到星际媒介中甘氨酸的形成, 与早期地球上的氨基酸起源有关. 如果没有考虑氢隧道效应, 在MP2/6-311+G(d,p)级别上计算反应能垒是254.7 kJ·mol-1, 在星际媒介中该气相反应很难进行. 在星际媒介冰颗粒表面上, 水分子催化反应增强了该化学反应的活性. H2NCH2CN与(H2O)3反应中的两个水分子作为催化剂降低活化能77.5 kJ·mol-1和活化自由能70.9 kJ·mol-1, 并且通过氢键桥协同传递质子. 量子氢隧道对于该反应进行至关紧要,采用小弯曲隧道(SCT)近似和正则变分过渡态理论(CVT)方法研究. 温度50 K时, 速率常数kSCT/CVT为1.86×10-23 cm3·molecule-1·s-1, 表明在星际媒介中通过质子隧道机理该反应容易进行. 研究结果与地球上的氨基酸起源于地球本身物质的观点相一致.     

Tunneling in the 1,5-hydrogen shift reactions of 1,3-cyclopentadiene and 5-methyl-1,3-cyclopentadiene

MPW1K/6-31+G(d,p) calculations which include the effects of small curvature tunneling find that, around room temperature, thermally activated tunneling dominates the 1,5-hydrogen shift reactions of 1,3-cyclopentadiene (2a) and 5-methyl-1,3-cyclopentadiene (2c). The calculated temperature dependence of the H/D kinetic isotope effect (KIE) for the latter rearrangement agrees well with experimental measurements that were published nearly 40 years ago. It is argued that the experimental KIEs provide prima facie evidence for tunneling in this reaction. The calculations also predict that it should be possible, at least in principle, to confirm this conclusion by observing curvature in the Arrhenius plot for the rearrangement of 2c.     

Large-amplitude dynamics in vinyl radical: the role of quantum tunneling as an isomerization mechanism

We report tunneling splittings associated with the large amplitude 1,2 H-atom migration to the global minima in the vinyl radical. These are obtained using a recent full-dimensional ab initio potential energy surface (PES) [A. R. Sharma, B. J. Braams, S. Carter, B. C. Shepler, and J. M. Bowman, J. Chem. Phys. 130(17), 174301 (2009)] and independently, directly calculated "reaction paths." The PES is a multidimensional fit to coupled cluster single and double and perturbative treatment of triple excitations coupled-cluster single double triple (CCSD(T)) with the augmented correlation consistent triple zeta basis set (aug-cc-pVTZ). The reaction path potentials are obtained from a series of CCSD(T)/aug-cc-pVnTZ calculations extrapolated to the complete basis set limit. Approximate 1D calculations of the tunneling splitting for these 1,2-H atom migrations are obtained using each of these potentials as well as quite different 1D Hamiltonians. The splittings are calculated over a large energy ranges, with results from the two sets of calculations in excellent agreement. Though negligibly slow (>1 s) for the vibrational ground state, this work predicts tunneling-promoted 1,2 hydride shift dynamics in vinyl to exhibit exponential growth with internal vibrational excitation, specifically achieving rates on the sub-μs time scale at energies above E ≈ 7500 cm(-1). Most importantly, these results begin to elucidate the possible role of quantum isomerization through barriers without dissociation, in competition with the more conventional picture of classical roaming permitted over a much narrower window of energies immediately below the bond dissociation limit. Furthermore, when integrated over a Boltzmann distribution of thermal energies, these microcanonical tunneling rates are consistent with sub-μs time scales for 1,2 hydride shift dynamics at T > 1400 K. These results have potential relevance for combustion modeling of low-pressure flames, as well as recent observations of nuclear spin statistical mixing from high-resolution IR/microwave spectroscopy on vinyl radical.     

The electronic structure of peracids. Functional models for cytochrome p-450

An extensive set of ground state ab initio and semiempirical molecular orbital calculations has been performed on both peroxytrifluoroacetic and peroxyacetic acids. The equilibrium geometry of peroxyacetic acid was calculated with the STO-3G and MINDO/3 methods, and peroxide rotational barriers for both peracids were obtained with STO-3G and PCILO. Extended basis set calculations with the 6-31G^** basis were performed for both peracids to compare the electronic structures of these two compounds. Electrostatic potential maps in the region of the peroxide bonds of both peracids were also calculated using INDO wavefunetions. These results are discussed with respect to the enhanced reactivity of peroxytrifluoracetic acid relative to peroxyacetic acid and the nature of the oxygen electrophilicity in these compounds and by analogy in cytochrome P-450, for which peroxytrifluoroacetic acid is considered to be an effective chemical model.     

Tunneling and classical paths for proton transfer in an enzyme reaction dominated by tunneling: oxidation of tryptamine by aromatic amine dehydrogenase

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.     

Mechanism and kinetics properties for the reaction: Chloroethane with atomic O (P)

In this paper, we present direct dynamics calculations for the multiple-channel reaction of CH₃CH₂Cl with atomic O (³P) in a wide temperature range (200–3000 K), based on canonical variational transition state theory including small curvature corrections. Four distinct saddle points, one for α-abstraction and three for β-abstraction, have been located for this reaction. The potential energy surface information has been calculated at the MP2/6-311G(d,p) level. The energies along the minimum energy path have been further improved by single-point energy calculations at the G3MP2 level. In the β-abstraction channel, Jahn–Teller effect has been found. Changes of geometries, generalized normal-mode vibrational frequencies, and potential energies along the reaction paths for all channels have been discussed and compared. The calculated total rate constants match the available experimental values reasonable well over the measured temperature range. The results show the variational effect can be negligible and the small curvature tunneling contribution plays an important role for the calculation of the rate constant. At low temperature α-abstraction may be the major reaction channel, while β-abstraction will have more contribution to the whole reaction rate as the temperature increase.     

Direct dynamics studies on the hydrogen abstraction reactions CF3O+CH4(CD4)→CF3OH(CF3OD)+CH3(CD3)

The dynamics properties of the hydrogen abstraction reaction CF₃O+CH₄→CF₃OH+CH₃ are studied by dual-level direct dynamics method. Optimization calculations are preformed by B3LYP and MP2 with the 6-311G(d,p) basis set, and the single-point calculations are done at the multi-coefficient correction method based on quadratic configuration interaction with single and double excitations (MC-QCISD) method. The rate constants are evaluated by canonical variational transition-state theory with a small-curvature tunneling correction over a wide range of temperature 200–2000 K. The agreement between theoretical and experimental rate constants is good in the measured temperature range. The calculated results show that the variational effect is small and almost neglected over the whole temperature range, whereas, the tunneling correction plays a role in the lower temperature range. The kinetic isotope effect for the reaction is ‘normal’. The value of k_H/k_D is 2.38 at room temperature and it decreases with the temperature increasing.