Theoretical study of oxidation of cyclohexane diol to adipic anhydride by [Ru(IV)(O)(tpa)(H2O)]2+ complex (tpa ═ tris(2-pyridylmethyl)amine)

The catalytic conversion of 1,2-cyclohexanediol to adipic anhydride by Ru(IV)O(tpa) (tpa ═ tris(2-pyridylmethyl)amine) is discussed using density functional theory calculations. The whole reaction is divided into three steps: (1) formation of α-hydroxy cyclohexanone by dehydrogenation of cyclohexanediol, (2) formation of 1,2-cyclohexanedione by dehydrogenation of α-hydroxy cyclohexanone, and (3) formation of adipic anhydride by oxygenation of cyclohexanedione. In each step the two-electron oxidation is performed by Ru(IV)O(tpa) active species, which is reduced to bis-aqua Ru(II)(tpa) complex. The Ru(II) complex is reactivated using Ce(IV) and water as an oxygen source. There are two different pathways of the first two steps of the conversion depending on whether the direct H-atom abstraction occurs on a C-H bond or on its adjacent oxygen O-H. In the first step, the C-H (O-H) bond dissociation occurs in TS1 (TS2-1) with an activation barrier of 21.4 (21.6) kcal/mol, which is followed by abstraction of another hydrogen with the spin transition in both pathways. The second process also bifurcates into two reaction pathways. TS3 (TS4-1) is leading to dissociation of the C-H (O-H) bond, and the activation barrier of TS3 (TS4-1) is 20.2 (20.7) kcal/mol. In the third step, oxo ligand attack on the carbonyl carbon and hydrogen migration from the water ligand occur via TS5 with an activation barrier of 17.4 kcal/mol leading to a stable tetrahedral intermediate in a triplet state. However, the slightly higher energy singlet state of this tetrahedral intermediate is unstable; therefore, a spin crossover spontaneously transforms the tetrahedral intermediate into a dione complex by a hydrogen rebound and a C-C bond cleavage. Kinetic isotope effects (k(H)/k(D)) for the electronic processes of the C-H bond dissociations calculated to be 4.9-7.4 at 300 K are in good agreement with experiment values of 2.8-9.0.     

Catalytic mechanism of class A beta-lactamase. I. The role of Glu166 and Serl30 in the deacylation reaction

The tetrahedral intermediate formation process, which is the first step in the deacylation reaction by class A beta-lactamase, was investigated by the ab initio molecular orbital method. In this study, benzyl penicillin was used as the substrate. From the results of our molecular dynamics study of the structure of beta-lactam antibiotics-beta-lactamase complex, the substrate, Ser70, Lys73, Ser130, Glu166 and a water molecule for the deacylation reaction were considered for construction of a model for calculation. The calculation results indicated that Glu166 plays a role in holding a water molecule, which is necessary for the deacylation reaction, and that the hydrogen bond network among Lys73Nzeta, Ser130Ogamma, and the carboxyl group of the beta-lactam antibiotics was formed by the uptake of beta-lactam antibiotics by beta-lactamase. The activation energy for this reaction was 33.3 kcal/mol, and it is very likely that the reaction occurred at body temperature. Subsequent calculation results obtained by using the model excluding Ser130 and the carboxyl group of the substrate indicated that the activation energy for this reaction was 40.8 kcal/mol, which is 7.5 kcal/mol higher than that of the previous reaction. It was found that the hydrogen bond network plays an important role in decreasing the activation energy for the tetrahedral intermediate formation reaction. Lys73Nzeta, which is located at the edge of the hydrogen bond network, played a role in forming a hydrogen bond with Glu166Oepsilon in order to help the deacylation reaction. The role of amino acid residues around the active site of class A beta-lactamase was also discussed.     

Hydrogen and deuterium exchange on Pt(111) and its poisoning by carbon monoxide studied by surface sensitive high-pressure techniques

Catalytic hydrogen/deuterium exchange on a platinum (111) single crystal and its poisoning with carbon monoxide was studied using scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), sum frequency generation vibrational spectroscopy (SFG), and mass spectrometry under reaction conditions at pressures in the mTorr to atmospheric range. At 298 K and in the presence of 200 mTorr of hydrogen and 20 mTorr of deuterium the surface is catalytically active, producing HD with activation energy of approximately 5.3 kcal/mol. Addition of 5 mTorr of CO stops the reaction completely. In situ STM images reveal an ordered surface structure of chemisorbed CO. At 353 K the addition of 5 mTorr of CO slows the reaction by 3 orders of magnitude, but HD production continues with an activation energy of 17.4 kcal/mol. Changes in coverage and adsorption site of CO were followed by XPS and SFG up to a temperature of 480 K. From these data, a CO dominated, mobile and catalytically active surface model is proposed.     

Chemistry of the t-butoxyl radical: evidence that most hydrogen abstractions from carbon are entropy-controlled

Absolute rate constants and Arrhenius parameters for hydrogen abstractions (from carbon) by the t-butoxyl radical ((t) BuO.) are reported for several hydrocarbons and tertiary amines in solution. Combined with data already in the literature, an analysis of all the available data reveals that most hydrogen abstractions (from carbon) by (t) BuO. are entropy controlled (i.e., TdeltaS > deltaH, in solution at room temperature). For substrates with C-H bond dissociation energies (BDEs) > 92 kcal/mol, the activation energy for hydrogen abstraction decreases with decreasing BDE in accord with the Evans-Polanyi equation, with alpha approximately 0.3. For substrates with C-H BDEs in the range from 79 to 92 kcal/mol, the activation energy does not vary significantly with C-H BDE. The implications of these results in the context of the use of (t) BuO. as a chemical model for reactive oxygen-centered radicals is discussed.     

The cocondensation reaction of lithium atoms and anisole: an experimental and theoretical study of the reaction pathway

The cocondensation reaction of lithium atoms and pure anisole leads to an ortho CH activation and the formation of lithium hydride. This simple two-component system allows the investigation of the reaction mechanism with included donor molecules. Therefore two anisole and one dilithium molecule, which was identified in an earlier spectroscopic study, were considered for the reaction pathway calculations.Firstly, two intermediates can be found along the reaction pathway, which show the reaction before and after the critical CH activation step. Secondly, a low-lying transition state can be identified, which allows the carbon hydrogen bond to be broken with an activation energy of less than 20 kcal/mol instead of more than 100 kcal/mol, if a free radical mechanism is employed. All calculations were performed at the B3LYP/6-31G** level of theory.     

The "somersault" mechanism for the p-450 hydroxylation of hydrocarbons. The intervention of transient inverted metastable hydroperoxides

A series of model theoretical calculations are described that suggest a new mechanism for the oxidation step in enzymatic cytochrome P450 hydroxylation of saturated hydrocarbons. A new class of metastable metal hydroperoxides is described that involves the rearrangement of the ground-state metal hydroperoxide to its inverted isomeric form with a hydroxyl radical hydrogen bonded to the metal oxide (MO-OH --> MO....HO). The activation energy for this somersault motion of the FeO-OH group is 20.3 kcal/mol for the P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] to produce the isomeric ferryl oxygen hydrogen bonded to an *OH radical [Por(SH)Fe(III)-O....HO(-)]. This isomeric metastable hydroperoxide, the proposed primary oxidant in the P450 hydroxylation reaction, is calculated to be 17.8 kcal/mol higher in energy than the ground-state iron(III) hydroperoxide Cpd 0. The first step of the proposed mechanism for isobutane oxidation is abstraction of a hydrogen atom from the C-H bond of isobutane by the hydrogen-bonded hydroxyl radical to produce a water molecule strongly hydrogen bonded to anionic Cpd II. The hydroxylation step involves a concerted but nonsynchronous transfer of a hydrogen atom from this newly formed, bound, water molecule to the ferryl oxygen with a concomitant rebound of the incipient *OH radical to the carbon radical of isobutane to produce the C-O bond of the final product, tert-butyl alcohol. The TS for the oxygen rebound step is 2 kcal/mol lower in energy than the hydrogen abstraction TS (DeltaE() = 19.5 kcal/mol). The overall proposed new mechanism is consistent with a lot of the ancillary experimental data for this enzymatic hydroxylation reaction.     

Theoretical prediction of S-H bond rupture in methanethiol upon interaction with gold

Organic thiols are known to react with gold surface to form self-assembled monolayers (SAMs), which can be used to produce materials with highly attractive properties. Although the structure of various SAMs is widely investigated, some aspects of their formation still represent a matter of debate. One of these aspects is the mechanism of S-H bond dissociation in thiols upon interaction with gold. This work presents a new suggestion for this mechanism on the basis of DFT study of methanethiol interaction with a single gold atom and a Au(20) cluster. The reaction path of dissociation is found to be qualitatively independent of the model employed. However, the highest activation barrier of S-H bond dissociation on the single gold atom (12.9 kcal/mol) is considerably lower than that on the Au(20) cluster (28.9 kcal/mol), which can be attributed to the higher extent of gold unsaturation. The energy barrier of S-H cleavage decreases by 4.6 kcal/mol in the presence of the second methanethiol molecule at the same adsorption site on the model gold atom. In the case of the Au(20) cluster we have observed the phenomenon of hydrogen transfer from one methanethiol molecule to another, which allows reducing the energy barrier of dissociation by 9.1 kcal/mol. This indicates the possibility of the "relay" hydrogen transfer to be the key step of the thiol adsorption observed for the SAMs systems.     

Phosphorylation reaction in cAPK protein kinase-free energy quantum mechanical/molecular mechanics simulations

We present results of a theoretical analysis of the phosphorylation reaction in cAMP-dependent protein kinase using a combined quantum mechanical and molecular mechanics (QM/MM) approach. Detailed analysis of the reaction pathway is provided using a novel QM/MM implementation of the nudged elastic band method, finite temperature fluctuations of the protein environment are taken into account using free energy calculations, and an analysis of hydrogen bond interactions is performed on the basis of calculated frequency shifts. The late transfer of the substrate proton to the conserved aspartate (D166), the activation free energy of 15 kcal/mol, and the slight exothermic (-3 kcal/mol) character of the reaction are all consistent with the experimental data. The near attack conformation of D166 in the reactant state is maintained by interactions with threonine-201, asparagine-177, and most notably by a conserved water molecule serving as a strong structural link between the primary metal ion and the D166. The secondary Mg ion acts as a Lewis acid, attacking the beta-gamma bridging oxygen of ATP. This interaction, along with a strong hydrogen bond between the D166 and the substrate, contributes to the stabilization of the transition state. Lys-168 maintains a hydrogen bond to a transferring phosphoryl group throughout a reaction process. This interaction increases in the product state and contributes to its stabilization.     

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9

To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reaction in aqueous solution, respectively. Our calculated activation free energy barrier for the enzyme catalyzed reaction is 22.5 kcal/mol, which agrees very well with the experimental value of 20.9 kcal/mol. The difference in potential of mean force between a corresponding prereaction state and the transition state for the solution reaction is computed to be 30.9 kcal/mol. Thus, our simulations indicate that the enzyme SET7/9 plays an essential catalytic role in significantly lowering the barrier for the methyl-transfer reaction step. For the reaction in solution, it is found that the hydrogen bond network near the reaction center undergoes a significant change, and there is a strong shift in electrostatic field from the prereaction state to the transition state, whereas for the enzyme reaction, such an effect is much smaller and the enzyme SET7/9 is found to provide a preorganized electrostatic environment to facilitate the methyl-transfer reaction. Meanwhile, we find that the transition state in the enzyme reaction is a little more dissociative than that in solution.     

Determination of the C4-H bond dissociation energies of NADH models and their radical cations in acetonitrile

Heterolytic and homolytic bond dissociation energies of the C4-H bonds in ten NADH models (seven 1,4-dihydronicotinamide derivatives, two Hantzsch 1,4-dihydropyridine derivatives, and 9,10-dihydroacridine) and their radical cations in acetonitrile were evaluated by titration calorimetry and electrochemistry, according to the four thermodynamic cycles constructed from the reactions of the NADH models with N,N,N',N'-tetramethyl-p-phenylenediamine radical cation perchlorate in acetonitrile (note: C9-H bond rather than C4-H bond for 9,10-dihydroacridine; however, unless specified, the C9-H bond will be described as a C4-H bond for convenience). The results show that the energetic scales of the heterolytic and homolytic bond dissociation energies of the C4-H bonds cover ranges of 64.2-81.1 and 67.9-73.7 kcal mol(-1) for the neutral NADH models, respectively, and the energetic scales of the heterolytic and homolytic bond dissociation energies of the (C4-H)(.+) bonds cover ranges of 4.1-9.7 and 31.4-43.5 kcal mol(-1) for the radical cations of the NADH models, respectively. Detailed comparison of the two sets of C4-H bond dissociation energies in 1-benzyl-1,4-dihydronicotinamide (BNAH), Hantzsch 1,4-dihydropyridine (HEH), and 9,10-dihydroacridine (AcrH(2)) (as the three most typical NADH models) shows that for BNAH and AcrH(2), the heterolytic C4-H bond dissociation energies are smaller (by 3.62 kcal mol(-1)) and larger (by 7.4 kcal mol(-1)), respectively, than the corresponding homolytic C4-H bond dissociation energy. However, for HEH, the heterolytic C4-H bond dissociation energy (69.3 kcal mol(-1)) is very close to the corresponding homolytic C4-H bond dissociation energy (69.4 kcal mol(-1)). These results suggests that the hydride is released more easily than the corresponding hydrogen atom from BNAH and vice versa for AcrH(2), and that there are two almost equal possibilities for the hydride and the hydrogen atom transfers from HEH. Examination of the two sets of the (C4-H)(.+) bond dissociation energies shows that the homolytic (C4-H)(.+) bond dissociation energies are much larger than the corresponding heterolytic (C4-H)(.+) bond dissociation energies for the ten NADH models by 23.3-34.4 kcal mol(-1); this suggests that if the hydride transfer from the NADH models is initiated by a one-electron transfer, the proton transfer should be more likely to take place than the corresponding hydrogen atom transfer in the second step. In addition, some elusive structural information about the reaction intermediates of the NADH models was obtained by using Hammett-type linear free-energy analysis.     

Is the peptide bond formation activated by Cu(2+) interactions? Insights from density functional calculations

The catalytic role that Cu(2+) cations play in the peptide bond formation has been addressed by means of density functional calculations. First, the Cu(2+)-(glycine)2 --> Cu(2+)-(glycylglycine) + H2O reaction was investigated since mass spectrometry low collision activated dissociation (CAD) spectra of Cu(2+)-(glycine)2 led to the elimination of a water molecule, which suggested that an intracomplex peptide bond formation might have occurred. Results show that this intracomplex condensation is associated to a very high free energy barrier (97 kcal mol(-1)) and reaction free energy (66 kcal mol(-1)) because of the loss of metal coordination during the reaction. Second, on the basis of the salt-induced peptide formation theory, the condensation reaction between two glycines was studied in aqueous solution using discrete water molecules and the conductor polarized continuum model (CPCM) continuous method. It is found that the synergy between the interaction of glycines with Cu(2+) and the presence of water molecules acting as proton-transfer helpers significantly lower the activation barrier (from 55 kcal/mol for the uncatalyzed system to 20 kcal/mol for the Cu(2+) solvated system) which largely favors the formation of the peptide bond.     

Mechanistic study on the gold-catalyzed C-S bond formation of α-thioallenes to form 2,5-dihydrothiophenes

The reaction mechanism of the AuCl-catalyzed reaction of the α-thioallenes to give 2,5-dihydrothiophenes has been computationally studied using DFT (B3LYP/6-31G*, SDD for Au). Calculations indicate the complexation of α-thioallene with AuCl occurs preferentially at the distal double bond, followed by the C-S bond formation, the proton transfer from the sulfur to the carbon "b", and the [1,2]-hydride shift to give the 2,5-dihydrothiophene gold complex. The proton transfer is the rate-limiting step with very high activation energy in the gas phase. In the presence of one water molecule, the activation free energy of the proton transfer was lowered by as much as 19.9 kcal/mol. Furthermore, one dichloromethane molecule stabilized all of the transition structures by its hydrogen bonds.     

Theoretical study of ethynylbenzene adsorption on Au(111) and implications for a new class of self-assembled monolayer

Density functional calculations of the adsorption of ethynylbenzene on the Au(111) surface show that, after cleavage of the C-H bond, the terminal carbon makes a strong covalent bond to the surface. The bond energy is shown to be about 70 kcal.mol(-1) with the fcc hollow site being most stable and the molecule oriented perpendicular to the surface. Adsorption without elimination of hydrogen is also possible via a hydrogen 1,2 shift to form a vinylidene surface-bound species, or opening of the C-C triple bond and adsorption through the two carbon atoms in a flat conformation. The reaction energy for formation of the surface-bound vinylidene is estimated to be 5 kcal.mol(-1) exothermic relative to the isolated ethynylbenzene and gold substrate.     

Identification of decomposition reactions for HMDSO organosilicon using quantum chemical calculations

Hexamethyldisiloxane [HMDSO, (CH₃)₃-SiOSi-(CH₃)₃] is an important precursor for SiO₂ formation during flame-based silica material synthesis. As a result, HMDSO reactions in flame have been widely investigated experimentally, and many results have indicated that HMDSO decomposition reactions occur very early in this process. In this paper, quantum chemical calculations are performed to identify the initial decomposition of HMDSO and its subsequent reactions using the density functional theory at the level of B3LYP/6-311+G (d, p). Four reaction pathways—(a) Si O bond dissociation of HMDSO, (b) Si C bond dissociation of HMDSO, (c) dissociation and recombination of Si O and Si C bonds, and (d) elimination of a methane molecule from HMDSO—have been examined and identified. From the results, it is found that the barrier of 84.38 kcal/mol and Si O bond dissociation energy of 21.55 kcal/mol are required for the initial decomposition reaction of HMDSO in the first pathway, but the highest free energy barrier (100.69 kcal/mol) is found in the third reaction pathway. By comparing the free energy barriers and reaction rate constants, it is concluded that the most possible initial decomposition reaction of HMDSO is to eliminate the CH₃ radical by Si C bond dissociation.     

Isolation, characterization of an intermediate in an oxygen atom-transfer reaction, and the determination of the bond dissociation energy

The synthesis and structure of a phosphine oxide-bound intermediate molecule originating from a dioxo-molybdenum(VI) complex is described. The loss of phosphine oxide has been followed by surface-induced dissociation mass spectrometry that gave the bond dissociation energy of 29.5 (+/- 3.5) kcal/mol. Considering the bond dissociation energy for a Mo=O bond to be 100 kcal/mol, this value suggests that the primary driving force for the reaction is the formation of the intermediate complex.     

Why urea eliminates ammonia rather than hydrolyzes in aqueous solution

A joint QM/MM and ab initio study on the decomposition of urea in the gas phase and in aqueous solution is reported. Numerous possible mechanisms of intramolecular decomposition and hydrolysis have been explored; intramolecular NH3 elimination assisted by a water molecule is found to have the lowest activation energy. The solvent effects were elucidated using the TIP4P explicit water model with free energy perturbation calculations in conjunction with QM/MM Monte Carlo simulations. The explicit representation of the solvent was found to be essential for detailed resolution of the mechanism, identification of the rate-determining step, and evaluation of the barrier. The assisting water molecule acts as a hydrogen shuttle for the first step of the elimination reaction. The forming zwitterionic intermediate, H3NCONH, participates in 8-9 hydrogen bonds with water molecules. Its decomposition is found to be the rate-limiting step, and the overall free energy of activation for the decomposition of urea in water is computed to be approximately 37 kcal/mol; the barrier for hydrolysis by an addition/elimination mechanism is found to be approximately 40 kcal/mol. The differences in the electronic structure of the transition states of the NH3 elimination and hydrolysis were examined via natural bond order analysis. Destruction of urea's resonance stabilization during hydrolysis via an addition/elimination mechanism and its preservation in the rearrangement to the H3NCONH intermediate were identified as important factors in determining the preferred reaction route.     

Reactivity of alkanes on zeolites: a computational study of propane conversion reactions

In this work, quantum chemical methods were used to study propane conversion reactions on zeolites; these reactions included protolytic cracking, primary hydrogen exchange, secondary hydrogen exchange, and dehydrogenation reactions. The reactants, products, and transition-state structures were optimized at the B3LYP/6-31G level and the energies were calculated with CBS-QB3, a complete basis set composite energy method. The computed activation barriers were 62.1 and 62.6 kcal/mol for protolytic cracking through two different transition states, 30.4 kcal/mol for primary hydrogen exchange, 29.8 kcal/mol for secondary hydrogen exchange, and 76.7 kcal/mol for dehydrogenation reactions. The effects of basis set for the geometry optimization and zeolite acidity on the reaction barriers were also investigated. Adding extra polarization and diffuse functions for the geometry optimization did not affect the activation barriers obtained with the composite energy method. The largest difference in calculated activation barriers is within 1 kcal/mol. Reaction activation barriers do change as zeolite acidity changes, however. Linear relationships were found between activation barriers and zeolite deprotonation energies. Analytical expressions for each reaction were proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts, as long as the deprotonation energies are first acquired.     

On the mechanisms of oxidation of organic sulfides by H2O2 in aqueous solutions

The mechanism of oxidation of organic sulfides in aqueous solutions by hydrogen peroxide was investigated via ab initio calculations. Specifically, two reactions, hydrogen transfer of hydrogen peroxide to form water oxide and the oxidation of dimethyl sulfide (DMS) by hydrogen peroxide to form dimethyl sulfoxide, were studied as models of these processes in general. Solvent effects are included both via including explicitly water molecules and via the polarizable continuum model. The former was found to have a much more significant effect than the latter. When explicit water molecules are included, a mechanism different from those proposed in the literature was found. Specific interactions including hydrogen bonding with 2-3 water molecules can provide enough stabilization for the charge separation of the activation complex. The energy barrier of the oxidation of DMS by hydrogen peroxide was estimated to be 12.7 kcal/mol, within the experimental range of the oxidation of analogous compounds (10-20 kcal/mol). The major reaction coordinates of the reaction are the breaking of the O-O bond of H2O2 and the formation of the S-O bond, the transfer of hydrogen to the distal oxygen of hydrogen peroxide occurring after the system has passed the transition state. Reaction barriers of the hydrogen transfer of H2O2 are an average of 10 kcal/mol or higher than the reaction barriers of the oxidation of DMS. Therefore, a two-step oxidation mechanism in which, first, the transfer of a hydrogen atom occurs to form water oxide and, second, the transfer of oxygen to the substrate occurs is unlikely to be correct. Our proposed oxidation mechanism does not suggest a pH dependence of oxidation rate within a moderate range around neutral pH (i.e., under conditions in which hydronium and hydroxide ions do not participate directly in the reaction), and it agrees with experimental observations over moderate pH values. Also, without including a protonated solvent molecule, it has activation energies that correspond to measured activation energies.     

Reaction mechanism and tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution: a combined Monte Carlo and quantum mechanics study

A combined Monte Carlo and quantum mechanical study was carried out to analyze the tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution. Second- and fourth-order M?ller-Plesset perturbation theory calculations indicate that in the gas phase thiol (Pym-SH) is more stable than the thione (Pym-NH) by ca. 8 kcal/mol. In aqueous solution, thermodynamic perturbation theory implemented on a Monte Carlo NpT simulation indicates that both the differential enthalpy and Gibbs free energy favor the thione form. The calculated differential enthalpy is DeltaH(SH)(-->)(NH)(solv) = -1.7 kcal/mol and the differential Gibbs free energy is DeltaG(SH)(-->)(NH)(solv) = -1.9 kcal/mol. Analysis is made of the contribution of the solute-solvent hydrogen bonds and it is noted that the SH group in the thiol and NH group in the thione tautomers act exclusively as a hydrogen bond donor in aqueous solution. The proton transfer reaction between the tautomeric forms was also investigated in the gas phase and in aqueous solution. Two distinct mechanisms were considered: a direct intramolecular transfer and a water-assisted mechanism. In the gas phase, the intramolecular transfer leads to a large energy barrier of 34.4 kcal/mol, passing through a three-center transition state. The proton transfer with the assistance of one water molecule decreases the energy barrier to 17.2 kcal/mol. In solution, these calculated activation barriers are, respectively, 32.0 and 14.8 kcal/mol. The solvent effect is found to be sizable but it is considerably more important as a participant in the water-assisted mechanism than the solvent field of the solute-solvent interaction. Finally, the calculated total Gibbs free energy is used to estimate the equilibrium constant.     

Reactivity of isobutane on zeolites: a first principles study

In this work, ab initio and density functional theory methods are used to study isobutane protolytic cracking, primary hydrogen exchange, tertiary hydrogen exchange, and dehydrogenation reactions catalyzed by zeolites. The reactants, products, and transition-state structures are optimized at the B3LYP/6-31G* level, and the final energies are calculated using the CBS-QB3 composite energy method. The computed activation barriers are 52.3 kcal/mol for cracking, 29.4 kcal/mol for primary hydrogen exchange, 29.9 kcal/mol for tertiary hydrogen exchange, and 59.4 kcal/mol for dehydrogenation. The zeolite acidity effects on the reaction barriers are also investigated by changing the cluster terminal Si-H bond lengths. The analytical expressions between activation barriers and zeolite deprotonation energies for each reaction are proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts.