CO(2) fixation by Rubisco: computational dissection of the key steps of carboxylation, hydration, and C-C bond cleavage.

Despite intensive experimental and computational studies, some important features of the mechanism of the photosynthetic CO(2)-fixing enzyme, Rubisco, are still not understood. To complement our previous investigation of the first catalytic step, the enolization of D-ribulose-1,5-bisphosphate (King et al., Biochemistry 1998, 44, 15414-15422), we present the first complete computational dissection of subsequent steps of the carboxylation reaction that includes the roles of the central magnesium ion and modeled residues of the active site. We investigated carboxylation, hydration, and C-C bond cleavage using the density functional method and the B3LYP/6-31G(d) level to perform geometry optimizations. The energies were determined by B3LYP/6-311+G(2d,p) single-point calculations. We modeled a fragment of the active site and substrate, taking into account experimental findings that the residues coordinated to the Mg ion, especially the carbamylated Lys-201, play critical roles in this reaction sequence. The carbamate appears to act as a general base, not only for enolization but also for hydration of the beta ketoacid formed by addition of CO(2) and, as well, cleavage of the C2-C3 bond of the hydrate. We show that CO(2) is added directly, without assistance of a Michaelis complex, and that hydration of the resultant beta ketoacid occurs in a separate subsequent step with a discrete transition state. We suggest that two conformations of the hydrate (gem-diol), with different metal coordination, are possible. The step with the highest activation energy during the carboxylation cycle is the C-C bond cleavage. Depending on the conformations of the gem-diol, different pathways are possible for this step. In either case, special arrangements of the metal coordination result in bond breaking occurring at remarkably low activation energies (between 28 and 37 kcal mol(-1)) which might be reduced further in the enzyme environment.     

The final catalytic step of cytochrome p450 aromatase: a density functional theory study

B3LYP density functional theory calculations are used to unravel the mysterious third step of aromatase catalysis. The feasibility of mechanisms in which the reduced ferrous dioxygen intermediate mediates androgen aromatization is explored and determined to be unlikely. However, proton-assisted homolysis of the peroxo hemiacetal intermediate to produce P450 compound I and the C19 gem-diol likely proceeds with a low energetic barrier. Mechanisms for the aromatization and deformylation sequence which are initiated by 1beta-hydrogen atom abstraction by P450 compound I are considered. 1beta-Hydrogen atom abstraction from substrates in the presence of the 2,3-enol encounters strikingly low barriers (5.3-7.8 kcal/mol), whereas barriers for this same process rise to 17.0-27.1 kcal/mol in the keto tautomer. Transition states for 1beta-hydrogen atom abstraction from enolized substrates in the presence of the 19-gem-diol decayed directly to the experimentally observed products. If the C19 aldehyde remains unhydrated, aromatization occurs with concomitant decarbonylation and therefore does not support dehydration of the C19 aldehyde prior to the final catalytic step. On the doublet surface, the transition state connects to a potentially labile 1(10) dehydrogenated product, which may undergo rapid aromatization, as well as formic acid. Ab initio molecular dynamics confirmed that the 1beta-hydrogen atom abstraction and deformylation or decarbonylation occur in a nonsynchronous, coordinated manner. These calculations support a dehydrogenase behavior of aromatase in the final catalytic step, which can be summarized by 1beta-hydrogen atom abstraction followed by gem-diol deprotonation.     

Reactions of 1,3-cyclohexadiene with singlet oxygen. A theoretical study.

A thorough study of the reaction of singlet oxygen with 1,3-cyclohexadiene has been made at the B3LYP/6-31G(d) and CASPT2(12e,10o) levels. The initial addition reaction follows a stepwise diradical pathway to form cyclohexadiene endoperoxide with an activation barrier of 6.5 kcal/mol (standard level = CASPT2(12e,10o)/6-31G(d); geometries and zero-point corrections at B3LYP/6-31G(d)), which is consistent with an experimental value of 5.5 kcal/mol. However, as the enthalpy of the transition structure for the second step is lower than the diradical intermediate, the reaction might also be viewed as a nonsynchronous concerted reaction. In fact, the concertedness of the reaction is temperature dependent since entropy differences create a free energy barrier for the second step of 1.8 kcal/mol at 298 K. There are two ene reactions; one is a concerted mechanism (DeltaH(double dagger) = 8.8 kcal/mol) to 1-hydroperoxy-2,5-cyclohexadiene (5), while the other, which forms 1-hydroperoxy-2,4-cyclohexadiene (18), passes through the same diradical intermediate (9) as found on the pathway to endoperoxide. The major pathway from the endoperoxide is O-O bond cleavage (22.0 kcal/mol barrier) to form a 1,4-diradical (25), which is 13.9 kcal/mol less stable than the endoperoxide. From the diradical, two low-energy pathways exist, one to epoxyketone (29) and the other to the diepoxide (27), where both products are known to be formed experimentally with a product ratio sensitive to the nature of substitutents. A significantly higher activation barrier leads to C-C bond cleavage and direct formation of maleic aldehyde plus ethylene.     

Heat of hydrogenation of 1,5-dehydroquadricyclane. A computational and experimental study of a highly pyramidalized alkene

The radical anion of the highly pyramidalized alkene 1,5-dehydroquadricyclane (1) was generated in the gas phase from the Squires reaction of 1,5-bis(trimethylsilyl)quadricyclane with F-/F2. The electron binding energy and proton affinity of 1*- were determined by bracketing experiments to be 0.6 +/- 0.1 eV and 386 +/- 5 kcal/mol, respectively. These values are in good agreement with values predicted by density functional theory (B3LYP/6-31+G*) and ab initio (CASPT2/6-31+G*) calculations. The experimental heat of hydrogenation of 1, obtained from a thermochemical cycle, was found to be 91 +/- 9 kcal/mol. This value of deltaH(H2) leads to values of 67 +/- 9 kcal/mol for the olefin strain energy (OSE) of 1, 172 +/- 9 kcal/mol for its heat of formation, and 23 +/- 9 kcal/mol for its pi bond dissociation enthalpy. Since the retro-Diels-Alder reaction of neutral 1 is computed to be highly exothermic, the finding that 1*- apparently does not undergo a retro-Diels-Alder reaction is of particular interest. The B3LYP/6-31+G* optimized geometry of 1 suggests that the bonding in this alkene is partially delocalized, presumably because the highly pyramidalized double bond in 1 interacts with the distal cyclopropane bonds in a manner that eventually leads to a retro-Diels-Alder reaction. The good agreement of the B3LYP and (2/2)CASPT2 values for the heat of hydrogenation and OSE of 1 with the experimentally derived values provides indirect evidence for the correctness of the B3LYP prediction that the equilibrium geometry of 1 lies part way along the reaction coordinate to the transition structure for the retro-Diels-Alder reaction.     

Mechanism of formation of the internal aldimine in pyridoxal 5'-phosphate-dependent enzymes

In this paper we studied the mechanism of formation of the internal aldimine, a common intermediate to most pyridoxal 5'-phosphate (PLP)-dependent enzymes. A large model based on the crystal structure from the human ornithine decarboxylase (ODC) enzyme was constructed and in total accounts for 504 atoms. The reaction mechanism was investigated using the ONIOM methodology (B3LYP/6-31G(d)//AM1), and the final energies were calculated with the M06/6-311++G(2d,2p)//B3LYP/6-31G(d) level of theory. It was demonstrated that the reaction is accomplished in three sequential steps: (i) the nucleophilic attack of Lysine69 to PLP, (ii) the carbinolamine formation, and (iii) a final dehydration step. For the carbinolamine formation, several mechanistic hypotheses were explored, and the preferred pathway assigns a key role for the conserved active site Cys360. The overall reaction is exergonic in -9.1 kcal/mol, and the rate-limiting step is the dehydration step (E(a) = 13.5 kcal/mol). For the first time, we provide an atomistic portrait of this mechanism in an enzymatic environment. Moreover, we were able to assign a novel role to Cys360 in the ODC reaction mechanism that was never proposed.     

Theoretical study of nucleophilic substitution at the disulfide bridge of cyclo-l-cystine

The reaction of cyclo-l-cystine with thiolate is examined at the B3LYP/6-31+G level. The two isomers of cyclo-l-cystine differ in their dihedral angle about the disulfide bond; the M isomer (with dihedral angle of -90.1 degrees) is found to be slightly lower in energy. The nucleophilic substitution reaction at sulfur follows the addition-elimination mechanism, exemplified by the hypercoordinate sulfur intermediate on the reaction surface. The reaction is exergonic (DeltaG = -6.16 kcal mol(-1)), and both the entrance and exit transition state lie below the reactant energies.     

Computational investigation of hydrogen abstraction from 2-aminoethanol by the 1,5-dideoxyribose-5-yl radical: a model study of a reaction occurring in the active site of ethanolamine ammonia lyase

Hydrogen abstraction from 2-aminoethanol by the 5'-deoxyadenosyl radical, which is formed upon Co--C bond homolysis in coenzyme B(12), was investigated by theoretical means with employment of the DFT (B3LYP) and ab initio (MP2) approaches. As a model system for the 5'-deoxyadenosyl moiety the computationally less demanding 1,5-dideoxyribose was employed; two conformers, which differ in ring conformation (C2- and C3-endo), were considered. If hydrogen is abstracted from "free" substrate by the C2-endo conformer of the 1,5-dideoxyribose-5-yl radical, the activation enthalpy is 16.7 kcal mol(-1); with the C3-endo counterpart, the value is 17.3 kcal mol(-1). These energetic requirements are slightly above the activation enthalpy limit (15 kcal mol(-1)) determined experimentally for the rate-determining step of the sequence, that is, hydrogen delivery from 5'-deoxyadenosine to the product radical. The activation enthalpy is lower when the substrate interacts with at least one amino acid from the active site. According to the computations, when a His model system partially protonates the substrate the activation enthalpy is 4.5 kcal mol(-1) for the C3-endo conformer and 5.8 kcal mol(-1) for the C2-endo counterpart. As hydrogen abstraction from the fully as well as the partially protonated substrate is preceded by the formation of quite stable encounter complexes, the actual activation barriers are around 13-15 kcal mol(-1). A synergistic interaction of 2-aminoethanol with two amino acids where His partially protonates the NH(2) group and Asp partially deprotonates the OH group of the substrate results in an activation enthalpy of 12.4 kcal mol(-1) for the C3-endo conformer and 13.2 kcal mol(-1) for the C2-endo counterpart. However, if encounter complexes exist in the active site, the actual activation barriers are much higher (>25 kcal mol(-1)) than that reported for the rate-determining step. These findings together with previous computations suggest that the energetics of the initial hydrogen abstraction decrease with an interaction of the substrate with only a protonating auxiliary, but for the rearrangement of the radical the synergistic effects of two auxiliaries are essential to pull the barrier below the limit of 15 kcal mol(-1).     

A DFT computational study of the bis-silylation reaction of acetylene catalyzed by palladium complexes

In this paper we have investigated at the DFT(B3LYP) level the catalytic cycle for the bis-silylation reaction of alkynes promoted by palladium complexes. A model-system formed by an acetylene molecule, a disilane molecule, and the Pd(PH(3))(2) complex has been used. The most relevant features of this catalytic cycle can be summarized as follows: (i) The first step of the cycle is an oxidative addition involving H(3)Si-SiH(3) and Pd(PH(3))(2). It occurs easily and leads to the cis (SiH(3))(2)Pd(PH(3))(2) complex that is 5.39 kcal mol(-1) lower in energy than reactants. (ii) The transfer of the two silyl groups to the C-C triple bond does not occur in a concerted way, but involves many steps. (iii) The cis (SiH(3))(2)Pd(PH(3))(2) complex, obtained from the oxidative addition, is involved in the formation of the first C-Si bond (activation barrier of 18.34 kcal mol(-1)). The two intermediates that form in this step cannot lead directly to the formation of the final bis(silyl)ethene product. However, they can isomerize rather easily (the two possible isomerizations have a barrier of 16.79 and 7.17 kcal mol(-1)) to new more stable species. In both these new intermediates the second silyl group is adjacent to the acetylene moiety and the formation of the second C-Si bond can occur rapidly leading to the (Z)-bis(silyl)ethene, as experimentally observed. (iv) The whole catalytic process is exothermic by 41.54 kcal mol(-1), in quite good agreement with the experimental estimate of this quantity (about 40 kcal mol(-1)).     

Molecule-induced alkane homolysis with dioxiranes.

The mechanisms of C-H and C-C bond activations with dimethyldioxirane (DMD) were studied experimentally and computationally at the B3LYP/6-311+G**//B3LYP/6-31G* density functional theory level for the propellanes 3,6-dehydrohomoadamantane (2) and 1,3-dehydroadamantane (3). The sigma(C-C) activation of 3 with DMD (Delta G(*) = 23.9 kcal mol(-1) and Delta G(r) = -5.4 kcal mol(-1)) is the first example of a molecule-induced homolytic C-C bond cleavage. The C-H bond hydroxylation observed for 2 is highly exergonic (Delta G(r) = -74.4 kcal mol(-1)) and follows a concerted pathway (Delta G(*) = 34.8 kcal mol(-1)), in contrast to its endergonic molecule-induced homolysis (Delta G(*) = 28.8 kcal mol(-1) and Delta G(r) = +9.2 kcal mol(-1)). The reactivities of 2 and 3 with CrO(2)Cl(2), which follow a molecule-induced homolytic activation mechanism, parallel the DMD results only for highly reactive 3, but differ considerably for more stable propellanes such as 4-phenyl-3,6-dehydrohomoadamantane (1) and 2.     

A theoretical study of reactivity and regioselectivity in the hydroxylation of adamantane by ferrate(VI)

The conversion of adamantane to adamantanols mediated by ferrate (FeO(4)(2)(-)), monoprotonated ferrate (HFeO(4)(-)), and diprotonated ferrate (H(2)FeO(4)) is discussed with the hybrid B3LYP density functional theory (DFT) method. Diprotonated ferrate is the best mediator for the activation of the C-H bonds of adamantane via two reaction pathways, in which 1-adamantanol is formed by the abstraction of a tertiary hydrogen atom (3 degrees ) and 2-adamantanol by the abstraction of a secondary hydrogen atom (2 degrees ). Each reaction pathway is initiated by a C-H bond cleavage via an H-atom abstraction that leads to a radical intermediate, followed by a C-O bond formation via an oxygen rebound step to lead to an adamantanol complex. The activation energies for the C-H cleavage step are 6.9 kcal/mol in the 1-adamantanol pathway and 8.4 kcal/mol in the 2-adamantanol pathway, respectively, at the B3LYP/6-311++G level of theory, whereas those of the second reaction step corresponding to the rebound step are relatively small. Thus, the rate-determining step in the two pathways is the C-H bond dissociation step, which is relevant to the regioselectivity for adamantane hydroxylation. The relative rate constant (3 degrees )/(2 degrees ) for the competing H-atom abstraction reactions is calculated to be 9.30 at 75 degrees C, which is fully consistent with an experimental value of 10.1.     

Theoretical Study on the Dehydrogenation Reaction of H2S by VS+ (3∑-)

The dehydrogenation reaction of H2S by the 3∑- ground state of VS+: VS+ + H2S → VS2+ + H2 has been studied by using Density Functional Theory (DFT) at the B3LYP/DZVP level. It is found that the reaction proceeds along two possible pathways (A and B) yielding two isomer dehydrogenation products VS2+-1 (3B2) and VS2+-2 (3A1), respectively. For both pathways,the reaction has a two-step-reaction mechanism that involves the migration of two hydrogen atoms from S2 to V+, respectively. The migration of the second hydrogen via TS3 and that of the first via TS4 are the rate-determining steps for pathways A and B, respectively. The activation energy is 17.4 kcal/mol for pathway A and 22.8 kcal/mol for pathway B relative to the reactants. The calculated reaction heat of 9.9 kcal/mol indicates the endothermicity of pathway A and that of -11.9 kcal/mol suggests the exothermicity of pathway B.     

Thermochemistry of Lewis adducts of BH3 and nucleophilic substitution of triethylamine on NH3BH3 in tetrahydrofuran

The thermochemistry of the formation of Lewis base adducts of BH(3) in tetrahydrofuran (THF) solution and the gas phase and the kinetics of substitution on ammonia borane by triethylamine are reported. The dative bond energy of Lewis adducts were predicted using density functional theory at the B3LYP/DZVP2 and B3LYP/6-311+G** levels and correlated ab initio molecular orbital theories, including MP2, G3(MP2), and G3(MP2)B3LYP, and compared with available experimental data and accurate CCSD(T)/CBS theory results. The analysis showed that the G3 methods using either the MP2 or the B3LYP geometries reproduce the benchmark results usually to within ~1 kcal/mol. Energies calculated at the MP2/aug-cc-pVTZ level for geometries optimized at the B3LYP/DZVP2 or B3LYP/6-311+G** levels give dative bond energies 2-4 kcal/mol larger than benchmark values. The enthalpies for forming adducts in THF were determined by calorimetry and compared with the calculated energies for the gas phase reaction: THFBH(3) + L → LBH(3) + THF. The formation of NH(3)BH(3) in THF was observed to yield significantly more heat than gas phase dative bond energies predict, consistent with strong solvation of NH(3)BH(3). Substitution of NEt(3) on NH(3)BH(3) is an equilibrium process in THF solution (K ≈ 0.2 at 25 °C). The reaction obeys a reversible bimolecular kinetic rate law with the Arrhenius parameters: log A = 14.7 ± 1.1 and E(a) = 28.1 ± 1.5 kcal/mol. Simulation of the mechanism using the SM8 continuum solvation model shows the reaction most likely proceeds primarily by a classical S(N)2 mechanism.     

An efficient methyltrioxorhenium(VII)-catalyzed transformation of hydrotrioxides (ROOOH) into dihydrogen trioxide (HOOOH)

Dihydrogen trioxide (HOOOH) is formed nearly quantitatively in the low-temperature (-70 degrees C) methyltrioxorhenium(VII) (MTO)-catalyzed transformation of silyl hydrotrioxides (R3SiOOOH), and some acetal hydrotrioxides, in various solvents, as confirmed by 1H, and 17O NMR spectroscopy. The calculated energetics (B3LYP) for the catalytic cycle, using H3SiOOOH as a model system, is consistent with the experimentally observed activation energy (9.5 +/- 2.0 kcal/mol) and a small kinetic solvent isotope effect (kH2O/kD2O = 1.1 +/- 0.1), indicating an initial concerted reaction between the silyl hydrotrioxide and MTO in the rate-determining step. With the addition of water in the next step, the intermediate undergoes a sigma-bond metathesis reaction to break the Re-OOOH bond and form HOOOH, together with the second dihydroxy intermediate. The final step in the catalytic cycle involves a second, catalytic water that lowers the barrier to form H3SiOH and MTO.     

Conformational selectivity in the diels-alder cycloaddition: predictions for reactions of s-trans-1,3-butadiene

Diels-Alder cycloaddition of s-trans-1,3-butadiene (1) should yield trans-cyclohexene (7), just as reaction of the s-cis conformer gives cis-cyclohexene (9). Investigation of this long-overlooked process with Hartree-Fock, Moller-Plesset, CASSCF, and DFT methods yielded in every case a C(2)-symmetric concerted transition state. At the B3LYP/6-31G (+ZPVE) level, this structure is predicted to be 42.6 kcal/mol above reactants, while the overall reaction is endothermic by 16.7 kcal/mol. A stepwise diradical process has been studied by UBLYP/6-31G theory and found to have barriers of 35.5 and 17.7 kcal/mol for the two steps. Spin correction lowers these values to 30.1 and 13.0 kcal/mol. The barrier to pi-bond rotation in cis-cyclohexene (9) is predicted (B3LYP theory) to be 62.4 kcal/mol, with trans-cyclohexene (7) lying 53.3 kcal/mol above cis isomer 9. Results suggest that pi-bond isomerization and concerted reaction may provide competitive routes for Diels-Alder cycloreversion. It is concluded that full understanding of the Diels-Alder reaction requires consideration of both conformers of 1,3-butadiene.     

Does Silica Surface Catalyse Peptide Bond Formation? New Insights from First‐Principles Calculations

The role that silica surface could have played in prebiotic chemistry as a catalyst for peptide bond formation has been addressed at the B3LYP/6-31+G(d,p) level for a model reaction involving glycine and ammonia on a silica cluster mimicking an isolated terminal silanol group present at the silica surface. Hydrogen-bond complexation between glycine and the silanol is followed by the formation of the mixed surface anhydride Si(surf)-O-C(=O)-R, which has been suggested in the literature to activate the C=O bond towards nucleophilic attack by a second glycine molecule, here simulated by the simpler NH3 molecule. However, B3LYP/6-31+G(d,p) calculations show that formation of the surface mixed anhydride Si(surf)-O-C(=O)-R is disfavoured (delta(r)G298 approximately 6 kcal mol(-1)), and that the surface bond only moderately lowers the free-energy barrier of the nucleophilic attack responsible for peptide bond formation (deltaG298(double dagger) approximately 48 kcal mol(-1)) in comparison with the uncatalysed reaction (deltaG298(double dagger) approximately 52 kcal mol(-1)). A further decrease of the free-energy barrier of peptide bond formation (deltaG298(double dagger) approximately 41 kcal mol(-1)) is achieved by a single water molecule close to the reaction centre acting as a proton-transfer helper in the activated complex. A possible role of strained silica surface defects on the formation of the surface mixed anhydride Si(surf)-O-C(=O)-R has also been addressed.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Quantum mechanical study of the ring-closing reaction of the hexatriene radical cation

The ring-closing reaction of hexatriene radical cation 1(*)(+) to 1,3-cyclohexadiene radical cation 2(*)(+) was studied computationally at the B3LYP/6-31G* and QCISD(T)/6-311G*//QCISD/6-31G* levels of theory. Both, concerted and stepwise mechanisms were initially considered for this reaction. Upon evaluation at the B3LYP level of theory, three of the possible pathways-a concerted C(2)-symmetric via transition structure 3(*)(+) and stepwise C(1)-symmetric pathways involving three-membered ring intermediate 5(*)(+) and four-membered ring intermediate 6(*)(+)-were rejected due to high-energy stationary points along the reaction pathway. The two remaining pathways were found to be of competing energy. The first proceeds through the asymmetric, concerted transition structure 4(*)(+) with an activation barrier E(a) = 16.2 kcal/mol and an overall exothermicity of -23.8 kcal/mol. The second pathway, beginning from the cis,cis,trans rotamer of 1(*)(+), proceeds by a stepwise pathway to the cyclohexadiene product with an overall exothermicity of -18.6 kcal/mol. The activation energy for the rate-determining step in this process, the formation of the intermediate bicyclo[3.1.0]hex-2-ene via transition structure 9(*)(+), was found to be 20.4 kcal/mol. More rigorous calculations of a smaller subsection of the potential energy hypersurface at the QCISD(T)//QCISD level confirmed these findings and emphasized the importance of conformational control of the reactant.     

DFT study on the mechanism of water-assisted dihydrogen elimination in group 6 octahedral metal hydride complexes

The reaction of water with octahedral bis-, tris- and tetrakis-(phosphine)tungsten, (phosphine)molybdenum and (phosphine)chromium complexes has been studied using B3LYP/def2-TZVPP level of DFT to elucidate dissociative, associative and hydride migratory insertion mechanisms for hydrogen elimination. In the dissociative mechanism, phosphine dissociation requires 19.3-28.5 kcal mol(-1) of energy. The phosphine-water ligand exchange is endergonic due to poor coordination ability of water to group 6 metals (binding energy 8.8-15.5 kcal mol(-1)). The ligand exchange leads to intermolecular M-HH(2)O dihydrogen interaction and facilitates dihydrogen elimination (G(act) = 6.8-15.5 kcal mol(-1)). In the associative mechanism, a water molecule in the first solvation shell interacts with the M-H bond through a dihydrogen bond (interaction energy 2.7-4.0 kcal mol(-1)) and leads to the elimination of H(2) by forming a hydroxide complex. Compared to the dissociative mechanism, G(act) of associative mechanisms are ～22 kcal mol(-1) higher. In the hydride migratory insertion mechanism, the hydride ligand shifts to the CO ligand (G(act) = 25.4-30.4 kcal mol(-1)) to afford a formyl complex and subsequently the H-H bond coupling occurs between formyl and water ligand (G(act) = 2.8-4.4 kcal mol(-1)). In many cases, the migratory insertion mechanism can simultaneously operate with the dissociative mechanism as a minor pathway, whereas owing to high G(act) value, the associative mechanism can be described as inactive in the reaction. The general argument that dihydrogen elimination is preceded by the formation of a dihydrogen intermediate is not applicable for the systems studied herein as the most favoured dissociative mechanism does not pass through such an intermediate. On the other hand, irrespective of the mechanisms, dihydrogen elimination invariably occurs with the formation of a dihydrogen bonded transition state. Our results also suggest that group 6 octahedral metal hydride complexes are attractive targets for the design of water splitting reactions.     

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.     

Experimental determination of the alpha and beta C--H bond dissociation energies in naphthalene

The acidities of the two different sites in naphthalene (1alpha and 1beta) and the electron affinities of the alpha- and beta-naphthyl radicals were measured using a Fourier transform mass spectrometer. Both carbon-hydrogen bond dissociation energies for naphthalene also were obtained, in this case via the application of a thermodynamic cycle. The final results are DeltaH(o)acid (1alpha) = 394.2+/-1.2 kcal mol(-1), DeltaH(o)acid (1beta) = 395.5+/-1.3 kcal mol(-1), EA(alpha) = 31.6+/-0.5 kcal mol(-1), EA(beta) = 31.6+/-0.5 kcal mol(-1), BDE(1alpha) = 112.2+/-1.3 kcal mol(-1) and BDE(1alpha) = 111.9+/-1.4 kcal mol(-1), and they are compared to benzene and phenyl radical as well as ab initio and density functional theory (B3LYP) calculations.