Hybrid density functional theory with specific reaction parameter: hydrogen abstraction reaction of fluoromethane by the hydroxyl radical

Three potential energy surfaces with specific reaction parameters are developed and tested for the OH + CH(3)F --> H(2)O + CH(2)F reaction. The goal of this work is to determine surfaces that provide calculated reaction rate constants that are comparable to the experimental data. The potential energy surfaces are constructed using hybrid and hybrid meta density functional theory methods, and the levels of electronic structure theory used in this study are mPW1PW91, B1B95, and mPW1B95 in conjunction with the 6-31+G(d,p) basis set. The reaction rate constants are calculated over the range 200-1,500 K using variational transition state theory with multidimensional tunneling contributions. New specific-reaction-parameter Hartree-Fock contributions are determined, and the hybrid density functional theory methods with these new contributions (35.5 +/- 1.2% for mPW1PW91, 36.6 +/- 1.2% for B1B95, and 40.7 +/- 1.2% for mPW1B95, respectively) reproduce accurate rate constants over an extended temperature range. On these potential energy surfaces, the classical barrier height for the hydrogen abstraction reaction is determined to be 3.4-3.8 kcal/mol, with a best estimate value of 3.6 kcal/mol.     

Canonical variational transition-state theory study of the CF3CH2CH3 + OH reaction

Variational transition-state theory rate constants with multidimensional tunneling contributions using the small curvature method have been calculated for the CF3CH2CH3 (HFC-263fb) + OH reaction over a temperature range from 200 to 373 K. The mPW1B95-41.0 hybrid functional, parametrized by Albu and Swaminathan to generate theoretical rate constants nearly identical to the experimental values for the CH3F + OH reaction, has been used in conjunction with the 6-31+G** basis set to explore the potential energy surface of the title reaction. The good agreement found between theoretical predictions and the experimental data available suggests that the present approach is an excellent option to obtain high-quality results at low computational cost for direct dynamics studies of hydrogen abstraction reactions from complex hydrofluorocarbons. The reliability of the structure activity relationship used to estimate rate constant values for OH reactions with hydrofluorocarbons is also discussed in detail.     

Hybrid density functional theory with a specific reaction parameter: hydrogen abstraction reaction of trifluoromethane by the hydroxyl radical

Potential energy surfaces are developed and tested for the OH + CHF₃ → H₂O + CF₃ reaction. The objective is to obtain surfaces that give calculated rate constants comparable to the experimental ones. The potential energy surfaces are constructed using hybrid and hybrid meta density functional theory methods (mPW1PW91, B1B95, and mPW1B95) with specific reaction parameters in conjunction with the 6–31+ G(d,p) basis set. The rate constants are calculated over the temperature range 200–1,500 K using variational transition state theory with multidimensional tunneling contributions. The hybrid density functional theory methods with specific-reaction-parameter Hartree-Fock exchange contributions (32.8–34.8% for mPW1PW91, 34.2–36.0% for B1B95, and 37.8–39.7% for mPW1B95, respectively) provide accurate rate constants over an extended temperature range. The classical barrier height for the hydrogen abstraction reaction is determined to be 6.5–6.9 kcal/mol on these potential energy surfaces, and the best estimate value is 6.77 kcal/mol. Electronic Supplementary Material  Supplementary material is available for this article at and is accessible for authorized users.     

Deuterium enrichment of interstellar methanol explained by atom tunneling

We calculate, down to low temperature and for different isotopes, the reaction rate constants for the hydrogen abstraction reaction H + H(3)COH → H(2) + CH(2)OH/CH(3)O. These explain the known abundances of deuterated forms of methanol in interstellar clouds, where CH(2)DOH can be almost as abundant as CH(3)OH. For abstraction from both the C- and the O-end of methanol, the barrier-crossing motion involves the movement of light hydrogen atoms. Consequently, tunneling plays a dominant role already at relatively high temperature. Our implementation of harmonic quantum transition state theory with on the fly calculation of forces and energies accounts for these tunneling effects. The results are in good agreement with previous semiclassical and quantum dynamics calculations (down to 200 K) and experimental studies (down to 295 K). Here we extend the rate calculations down to lower temperature: 30 K for abstraction from the C-end of methanol and 80 K for abstraction from the OH-group. At all temperatures, abstraction from the C-end is preferred over abstraction from the O-end, more strongly so at lower temperature. Furthermore, the tunneling behavior strongly affects the kinetic isotope effects (KIEs). D + H(3)COH → HD + CH(2)OH has a lower vibrationally adiabatic barrier than H + H(3)COH → H(2) + CH(2)OH, giving rise to an inverse KIE (k(H)/k(D) < 1) at high temperature, in accordance with previous experiments and calculations. However, since tunneling is more facile for the light H atom, abstraction by H is favored over abstraction by D below ~135 K, with a KIE k(H)/k(D) of 11.2 at 30 K. The H + D(3)COD → HD + CD(2)OD reaction is calculated to be much slower than the D + H(3)COH → HD + CH(2)OH, in agreement with low-temperature solid-state experiments, which suggests the preference for H (as opposed to D) abstraction from the C-end of methanol to be the mechanism by which interstellar methanol is deuterium-enriched.     

Enzyme dynamics and tunneling enhanced by compression in the hydrogen abstraction catalyzed by soybean lipoxygenase-1

A fully microscopical simulation of the rate-limiting hydrogen abstraction catalyzed by soybean lipoxygenase-1 (SLO-1) has been carried out. This enzyme exhibits the largest, and weakly temperature dependent, experimental H/D kinetic isotope effect (KIE) reported for a biological system. The theoretical model used here includes the complete enzyme with a solvation shell of water molecules, the Fe(III)-OH- cofactor, and the linoleic acid substrate. We have used a hybrid QM(PM3/d-SRP)/MM method to describe the potential energy surface of the whole system, and the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) to calculate the rate constant and the primary KIE. The computational results show that the compression of the wild-type active site enzyme results in the huge contribution of tunneling (99%) to the rate of the hydrogen abstraction. Importantly, the active site becomes more flexible in the Ile553Ala mutant reactant complex simulation (for which a markedly temperature dependent KIE has been experimentally determined), thus justifying the proposed key role of the gating promoting mode in the reaction catalyzed by SLO-1. Finally, the results indicate that the calculated KIE for the wild-type enzyme has an important dependence on the barrier width.     

Multidimensional tunneling in the [1,5] shift in (Z)-1,3-pentadiene: how useful are Swain-Schaad exponents at detecting tunneling?

Rates, kinetic isotope effects (KIE), and Swain-Schaad exponents (SSE) have been calculated for a variety of isotopologues for the [1,5] shift in (Z)-1,3-pentadiene using mPW1K/6-31+G(d,p). Quantum mechanical effects along the reaction coordinate were incorporated with the zero-curvature tunneling (ZCT) model and with the multidimensional small curvature tunneling (SCT) model, which allows for coupling of modes perpendicular to the reaction coordinate. The latter model gives the best agreement with experimental rates and primary KIEs. The small quasiclassical primary KIE (2.6) is rationalized in terms of a nonlinear transition state. For sp3 to sp2 rehybridization, the quasiclassical alpha-secondary KIE shows an unusual inverse effect due to compression of the nonbonding hydrogens in the suprafacial transition state. SCT transmission coefficients (kappa) increase the rates by as much as one order of magnitude. Tunneling allows the reactant to evade 1-2.5 kcal/mol of the barrier depending on the isotope. Inclusion of tunneling in the secondary KIE increases it beyond the equilibrium isotope effect and converts the inverse effect (0.95) into a normal KIE (1.12). Tunneling was found to deflate the primary y SSE but by an amount too small to distinguish it from the quasiclassical SSE. On the other hand, when a specific labeling pattern is used, the difference between the quasiclassical secondary SSE (4.1) and the tunneling secondary SSE (2.3) may be sufficiently large to detect tunneling. The mixed secondary SSE shows even larger differences.     

Small temperature dependence of the kinetic isotope effect for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase

The H/D primary kinetic isotope effect (KIE) for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase (ecDHFR) is calculated as a function of temperature employing ensemble-averaged variational transition-state theory with multidimensional tunneling. The calculated KIEs display only a small temperature dependence over the temperature range of 5 to 45 degrees C. We identify two key features that contribute to canceling most of the temperature dependence of the KIE that would be expected on the basis of simpler models. Related issues such as the isotope effects on Arrhenius preexponential factors, large differences between free energies of activation and Arrhenius activation energy, and fluctuations of effective barriers are also discussed.     

Theoretical prediction on the thermal stability of cyclic ozone and strong oxygen tunneling

Dual-level dynamics calculation with variational transition state theory including multidimensional tunneling has been performed on the isomerization reaction of cyclic ozone → normal (open) ozone, which was believed to be the stability-determining reaction of the elusive cyclic ozone molecule under thermal condition. The high-level potential energy surface data were obtained from the calculation using the MRCISD+Q theory with the aug-cc-pVQZ basis set, while the low-level reaction path information was obtained using the hybrid density functional theory B3LYP with the cc-pVTZ basis set. The calculated results showed very significant tunneling effects below 300 K (a factor of ～200 at 300 K and over 10(7) at 200 K). Because of the strong tunneling effects and the potential energy surface crossing of the 1A(1) and 1A(2) states, the isomerization reactions were found to be significantly faster than previously believed. The half-life of the cyclic ozone was estimated only ～10 s at 200 K and ～70 s below 100 K, which might partly explain the unsuccessful attempts for its experimental identification. The kinetic isotope effects (KIEs) for various (18)O substitution reactions were also calculated as a function of temperature and were as high as 10 at very low temperature. Because of the large KIEs, the experimental identification of the cyclic (18)O(3) seems more promising.     

Quantum dynamics of hydride transfer catalyzed by bimetallic electrophilic catalysis: synchronous motion of Mg(2+) and H(-) in xylose isomerase

Xylose isomerase exhibits a bridged-bimetallic active-site motif in which the substrate is bound to two metals connected by a glutamate bridge, and X-ray crystallographic studies suggest that metal movement is involved in the hydride transfer rate-controlling catalytic step. Here we report classical/quantal dynamical simulations of this step that provide new insight into the metal motion. The potential energy surface is calculated by treating xylose with semiempirical molecular orbital theory augmented by a simple valence bond potential and the rest of the system by molecular mechanics. The rate constant for the hydride-transfer step was calculated by ensemble-averaged dynamical simulations including both variational transition-state theory for determination of the statistically averaged dynamical bottleneck and optimized multidimensional tunneling calculations. The dynamics calculations include 25 317 atoms, with quantized vibrational free energy in 89 active-site degrees of freedom, and with 32 atoms moving through static secondary zone transition-state configurations in the quantum tunneling simulation. Our simulations show that the average Mg-Mg distance R increases monotonically as a function of the hydride-transfer progress variable z. The range of the average R along the reaction path is consistent with the X-ray structure, thus providing a dynamical demonstration of the postulated role of Mg in catalysis. We also predicted the primary deuterium kinetic isotope effect (KIE) for the chemical step. We calculated a KIE of 3.8 for xylose at 298 K, which is consistent with somewhat smaller experimentally observed KIEs for glucose substrate at higher temperatures. More than half of our KIE is due to tunneling; neglecting quantum effects on the reaction coordinate reduces the calculated KIE to 1.8.     

Kinetics of the hydrogen abstraction *CH3 + alkane --> CH4 + alkyl reaction class: an application of the reaction class transition state theory

Kinetics of the hydrogen abstraction reaction (*)CH(3) + CH(4) --> CH(4) + (*)CH(3) is studied by a direct dynamics method. Thermal rate constants in the temperature range of 300-2500 K are evaluated by the canonical variational transition state theory (CVT) incorporating corrections from tunneling using the multidimensional semiclassical small-curvature tunneling (SCT) method and from the hindered rotations. These results are used in conjunction with the Reaction Class Transition State Theory/Linear Energy Relationship (RC-TST/LER) to predict thermal rate constants of any reaction in the hydrogen abstraction class of (*)CH(3) + alkanes. Our analyses indicate that less than 40% systematic errors on the average exist in the predicted rate constants using the RC-TST/LER method while comparing to explicit rate calculations the differences are less than 100% or a factor of 2 on the average.     

Static and dynamic second hyperpolarizability calculated by time-dependent density functional cubic response theory with local contribution and natural bond orbital analysis

The static and dynamic second hyperpolarizability gamma has been investigated by time-dependent density functional cubic response theory. The third-order coupled perturbed Kohn-Sham equations were solved to obtain the third-order perturbed charge density. Calculations on a number of small molecules (N(2), CO(2), C(2)H(4), CO, HF, H(2)O, and CH(4)), paradisubstituted oligoacetylene chains, benzene, and eight paradisubstituted benzenes were performed to verify the implementation and to assess the accuracy of the nonhybrid and hybrid time-dependent density functional theory computations. Nitroaniline and a derivative were taken as examples to investigate the distribution of the "gamma density" and to demonstrate the feasibility of analyzing cubic response functions in terms of contributions from natural bond orbitals (NBOs) and natural localized molecular orbitals (NLMOs). The results highlight the contributions from atoms and bonds on different functional groups to the total value of gamma based on the NBO/NLMO analysis, which might be helpful for new nonlinear optical materials design.     

Temperature dependence of the kinetic isotope effects in thymidylate synthase. A theoretical study

In recent years, the temperature dependence of primary kinetic isotope effects (KIE) has been used as indicator for the physical nature of enzyme-catalyzed H-transfer reactions. An interactive study where experimental data and calculations examine the same chemical transformation is a critical means to interpret more properly temperature dependence of KIEs. Here, the rate-limiting step of the thymidylate synthase-catalyzed reaction has been studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) simulations in the theoretical framework of the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) combined with Grote-Hynes theory. The KIEs were calculated across the same temperature range examined experimentally, revealing a temperature independent behavior, in agreement with experimental findings. The calculations show that the H-transfer proceeds with ～91% by tunneling in the case of protium and ～80% when the transferred protium is replaced by tritium. Dynamic recrossing coefficients are almost invariant with temperature and in all cases far from unity, showing significant coupling between protein motions and the reaction coordinate. In particular, the relative movement of a conserved arginine (Arg166 in Escherichia coli ) promotes the departure of a conserved cysteine (Cys146 in E. coli ) from the dUMP by polarizing the thioether bond thus facilitating this bond breaking that takes place concomitantly with the hydride transfer. These promoting vibrations of the enzyme, which represent some of the dimensions of the real reaction coordinate, would limit the search through configurational space to efficiently find those decreasing both barrier height and width, thereby enhancing the probability of H-transfer by either tunneling (through barrier) or classical (over-the-barrier) mechanisms. In other words, the thermal fluctuations that are coupled to the reaction coordinate, together with transition-state geometries and tunneling, are the same in different bath temperatures (within the limited experimental range examined). All these terms contribute to the observed temperature independent KIEs in thymidylate synthase.     

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.     

Large tunneling effect on the hydrogen transfer in bis(μ-oxo)dicopper enzyme: a theoretical study

Type-III copper-containing enzymes have dicopper centers in their active sites and exhibit a novel capacity for activating aliphatic C-H bonds in various substrates by taking molecular oxygen. Dicopper enzyme models developed by Tolman and co-workers reveal exceptionally large kinetic isotope effects (KIEs) for the hydrogen transfer process, indicating a significant tunneling effect. In this work, we demonstrate that variational transition state theory allows accurate prediction of the KIEs and Arrhenius parameters for such model systems. This includes multidimensional tunneling based on state-of-the-art quantum-mechanical calculations of the minimum-energy path (MEP). The computational model of bis(μ-oxo)dicopper enzyme consists of 70 atoms, resulting in a 204-dimensional potential energy surface. The calculated values of E(a)(H) - E(a)(D), A(H)/A(D), and the KIE at 233 K are -1.86 kcal/mol, 0.51, and 28.1, respectively, for the isopropyl ligand system. These values agree very well with experimental values within the limits of experimental error. For the representative tunneling path (RTP) at 233 K, the pre- and post-tunneling configurations are 3.3 kcal/mol below the adiabatic energy maximum, where the hydrogen travels 0.54 ? by tunneling. We found that tunneling is very efficient for hydrogen transfer and that the RTP is very different from the MEP. It is mainly heavy atoms that move as the reaction proceeds from the reactant complex to the pretunneling configuration, and the hydrogen atom suddenly hops at that point.     

Theoretical Simulations of Kinetic Isotope Effects on Decarboxylation of 3-Carboxybenzisoxazole

Kinetic isotope effect values on the decarboxylation of 3-carboxybenzisoxazole have been computed using the second-order Kleinert's variational perturbation theory in the framework of Feynman's path integrals along with the potential energy surface obtained at the MP2/6-31+G(d) level. Good agreement with the experimental data was obtained, demonstrating that this novel computational approach for computing KIE values of organic reaction is a viable alternative to the traditional method employing the Bigeleisen equation and harmonic vibrational frequencies. Compared with the experimental measurements, consideration of anharmonicity and tunneling effects can significantly improve the calculated KIE values, reducing the root-mean-square deviation from 1.19 % for traditional method to 0.20 % for path-integral method.     

Influence of end group and surface structure on the current-voltage characteristics of alkanethiol monolayers on Au(111)

We present density functional theory calculations of the electronic structure and tunneling characteristics of alkanethiolate monolayers on Au(111). We systematically analyze radical3 x radical3 full coverage monolayers of SC6H12X molecules with different terminal groups, X=CH3, NH2, SH, OH, COOH, OCH3, on defect-free ("perfect") Au(111). We also study the influence of the surface-molecule bonding structure by comparing the properties of monolayers of SC6H12CH3 molecules on the perfect surface and on Au(111) surfaces with vacancies or adatoms. The tunneling currents (I) through the adsorbed monolayers with a single chemical contact have been calculated within the Tersoff-Hamann approach for voltages between -1 and +1 V. Computed currents are found to depend linearly on V at low voltage, with typical values of approximately 60 and 150 pA/molecule at 0.2 and 0.5 V, respectively, in good agreement with several experimental data. Computed tunneling currents show also a significant dependence on both the terminal group X and the surface structure. In particular, in order of decreasing intensities, currents for the different end groups are NH2 approximately SH>CH3>OH>OCH3>COOH. The relationships between the tunneling current, the work function of the surface+SAM, and the lineup of the HOMO with respect to the Fermi energy of the metal surface are examined.     

Adsorption Behavior of CH2 and CH3 on Metal Clusters Cun (n=1-6)

Using density functional theory with generalized gradient approximation and hybrid functional, we studied the properties of energy, charge population, and vibration of CH2 and CH3 adsorbed on Cun (n=1-6) clusters. The results show that the DFT calculation with the hybrid functional matches the experimental results better in both cases. The calculation results indicate that the adsorption of CH2 is stronger than that of CH3. During adsorption, the charges transfer from Cu to CH2 or CH3. The obtained vibrational frequencies for different modes of CH2 and CH3 adsorbed on Cun agree well with the experimental results for the adsorption on Cu(111) surface.     

Deep tunneling dominates the biologically important hydride transfer reaction from NADH to FMN in morphinone reductase

The temperature dependence of the primary kinetic isotope effect (KIE), combined temperature-pressure studies of the primary KIE, and studies of the alpha-secondary KIE previously led us to infer that hydride transfer from nicotinamide adenine dinucleotide to flavin mononucleotide in morphinone reductase proceeds via environmentally coupled hydride tunneling. We present here a computational analysis of this hydride transfer reaction using QM/MM molecular dynamics simulations and variational transition-state theory calculations. Our calculated primary and secondary KIEs are in good agreement with the corresponding experimental values. Although the experimentally observed KIE lies below the semiclassical limit, our calculations suggest that approximately 99% of the reaction proceeds via tunneling: this is the first "deep tunneling" reaction observed for hydride transfer. We also show that the dominant tunneling mechanism is controlled by the isotope at the primary rather than the secondary position: with protium in the primary position, large-curvature tunneling dominates, whereas with deuterium in this position, small-curvature tunneling dominates. Also, our study is consistent with tunneling being preceded by reorganization: in the reactant, the rings of the nicotinamide and isoalloxazine moieties are stacked roughly parallel to each other, and as the system moves toward a "tunneling-ready" configuration, the nicotinamide ring rotates to become almost perpendicular to the isoalloxazine ring.     

Dynamics of the concerted triple proton transfer in cyclic water trimer using the multiconfiguration molecular mechanics algorithm

Cyclic water clusters are important molecular species to understand the nature of hydrogen bonded networks. Theoretical studies for the dynamics of triple proton transfer in the cyclic water trimer were performed. The potential energy surface (PES) of triple proton transfer is generated by the multiconfiguration molecular mechanics (MCMM) algorithm. We have used the MP2/6-31G(d,p) level for high-level ab initio data (energies, gradients, and Hessians), which are used in the Shepard interpolation. Eight high-level reference points were added step by step, including two points for the critical configurations of the large curvature tunneling paths. The more high-level points are used, the better the potential energy surfaces become. The rate constant and kinetic isotope effect (KIE) for the triple proton transfer at 300 K, which have been calculated by the canonical variational transition-state theory with microcanonical optimized multidimensional semiclassical tunneling approximation, are 1.6 x 10(-3) s(-1) and 230, respectively. Tunneling is very important not only for the triple proton transfer but also for the triple deuterium transfer. The MCMM results show good agreement with those from the direct ab initio dynamics calculations.