Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Can the ebselen derivatives catalyze the isomerization of peroxynitrite to nitrate?

The reaction of ebselen and its derivatives (1-7) with peroxynitrite anion (ONOO(-); PN) has been studied in gas phase and in aqueous, dichloromethane, benzene, and cyclohexane solutions using B3LYP/6-311+G(d,p)//B3LYP/6-311G(d,p) and PCM-B3LYP/6-311+G(d,p)//B3LYP/6-311G(d,p) approaches, respectively. It was shown that the reaction of 2 (R=H) with PN proceeds via 2 + PN --> 2-PN --> 2-TS1 (O-O activation) --> 2-O(NO(2)(-)()) --> 2-SeO + NO(2)(-) pathway with a rate-determining barrier of 25.3 (14.8) kcal/mol at the NO(2)(-) dissociation step (numbers presented without parentheses are enthalpies, and those in parentheses are Gibbs free energies). The NO(3)(-) formation process, starting from the complex 2-O(NO(2)(-)()), requires by (7.9) kcal/mol more energy than the NO(2)(-) dissociation process and is unlikely to compete with the latter. Thus, in the gas phase, the peroxynitrite --> nitrate isomerization catalyzed by complex 2 is unlikely to occur. It is shown that the NO(3)(-) formation process is slightly more favorably than the NO(2)(-) dissociation process for complex 4, with a strongest electron-withdrawing ligand R=CF(3). Therefore, complex 4 (as well as complex 6 with R=OH) is predicted to be a good catalyst for peroxynitrite <--> nitrite isomerization in the gas phase. Solvent effects (a) change the rate-determining step of the reaction 2 + PN from NO(2)(-) dissociation in the gas phase to O-O activation, which occurs with barriers of (13.9), (8.4), (8.4), and (8.2) kcal/mol in water, dichloromethane, benzene, and cyclohexane, respectively, and (b) significantly reduce the NO(2)(-) dissociation energy, while only slightly destabilizing the NO(3)(-) formation barrier, and make the peroxynitrite <--> nitrate isomerization process practically impossible, even for complex 4.     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Mechanism of HCS + O2 reaction: hydrogen- or oxygen-transfer?

In spite of the potential importance of the HCS radical in both combustion and interstellar processes, its chemical reactivity has not been tackled previously. In the present paper, the oxidation reaction of the HCS radical is theoretically investigated for the first time at the CCSD(T)/6-311++G(3df,2p)//BH&HLYP/6-311++G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels. It is shown that the most feasible pathway is the O2 addition to the HCS radical forming the intermediate SC(H)OO which can undergo a subsequent O-extrusion leading to SC(H)O + 3O. This features an indirect O-transfer mechanism with the overall barrier of 4.4 and 3.5 kcal mol(-1), respectively, at the two levels. However, formation of the H-transfer product CS + HO2 is kinetically much less feasible, i.e., the direct mechanism has barriers of 14.3 and 8.7 kcal mol(-1), whereas the indirect mechanism has barriers of 12.6 and 10.7 kcal mol(-1), respectively. This result is in sharp contrast to the analogous HCO + O2 reaction, where the direct (with a barrier of 2.98 kcal mol(-1)) and indirect (2.26 kcal mol(-1)) H-transfer processes are highly competitive over the indirect O-transfer process (the least endothermicity is 19.9 kcal mol(-1)). The possible explanations and implications of the present results are provided.     

Gaseous reaction mechanism of C2F radical with water

The kinetic properties of the carbon-fluorine radicals are little understood except those of CFn (n =1-3). In this article, a detailed mechanistic study was reported on the gas-phase reaction between the simplest pi-bonded C2F radical and water as the first attempt to understand the chemical reactivity of the C2F radical. Various reaction channels are considered. The most kinetically competitive channel is the quasi-direct hydrogen-abstraction route forming P5 HCCF + OH. At the CCSD(T)/6-311+G(2d,2p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels, the overall H-abstraction barriers (4.5, 4.7, and 4.2 kcal/mol) for the C2F + H2O reaction are comparable to the corresponding values (5.5, 3.7, and 5.7 kcal/mol) for the analogous C2H + H2O reaction. This suggests that C2F is a reactive radical like the extensively studied C2H, in contrast to the situation of the CF and CF2 radicals that have much lower reactivity than the corresponding hydrocarbon species. Thus, the C2F radical is expected to play an important role in the combustion processes of the carbon-fluorine chemistry. Furthermore, addition of a second H2O can catalyze the reaction with the H-abstraction barrier significantly reduced to a marginally zero value (0.5 kcal/mol). This is also indicative of the potential relevance of the title reactions in the low-temperature atmospheric chemistry.     

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.     

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.     

Gas phase SN2 reactions of halide ions with trifluoromethyl halides: front- and back-side attack vs. complex formation

Density functional theory computations and pulsed-ionization high-pressure mass spectrometry experiments have been used to explore the potential energy surfaces for gas-phase S(N)2 reactions between halide ions and trifluoromethyl halides, X(-) + CF(3)Y --> Y(-) + CF(3)X. Structures of neutrals, ion-molecule complexes, and transition states show the possibility of two mechanisms: back- and front-side attack. From pulsed-ionization high-pressure mass spectrometry, enthalpy and entropy changes for the equilibrium clustering reactions for the formation of Cl(-)(BrCF(3)) (-16.5 +/- 0.2 kcal mol(-1) and -24.5 +/- 1 cal mol(-1) K(-1)), Cl(-)(ICF(3)) (-23.6 +/- 0.2 kcal mol(-1)), and Br(-)(BrCF(3)) (-13.9 +/- 0.2 kcal mol(-1) and -22.2 +/- 1 cal mol(-1) K(-1)) have been determined. These are in good to excellent agreement with computations at the B3LYP/6-311+G(3df)//B3LYP/6-311+G(d) level of theory. It is shown that complex formation takes place by a front-side attack complex, while the lowest energy S(N)2 reaction proceeds through a back-side attack transition state. This latter mechanism involves a potential energy profile which closely resembles a condensed phase S(N)2 reaction energy profile. It is also shown that the Cl(-) + CF(3)Br --> Br(-) + CF(3)Cl S(N)2 reaction can be interpreted using Marcus theory, in which case the reaction is described as being initiated by electron transfer. A potential energy surface at the B3LYP/6-311+G(d) level of theory confirms that the F(-) + CF(3)Br --> Br(-) + CF(4) S(N)2 reaction proceeds through a Walden inversion transition state.     

Corannulene as a Lewis base: computational modeling of protonation and lithium cation binding.

A computational modeling of the protonation of corannulene at B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p) and of the binding of lithium cations to corannulene at B3LYP/6-311G(d,p)//B3LYP/6-31G(d,p) has been performed. A proton attaches preferentially to one carbon atom, forming a sigma-complex. The isomer protonated at the innermost (hub) carbon has the best total energy. Protonation at the outermost (rim) carbon and at the intermediate (bridgehead rim) carbon is less favorable by ca. 2 and 14 kcal mol(-)(1), respectively. Hydrogen-bridged isomers are transition states between the sigma-complexes; the corresponding activation energies vary from 10 to 26 kcal mol(-)(1). With an empirical correction obtained from calculations on benzene, naphthalene, and azulene, the best estimate for the proton affinity of corannulene is 203 kcal mol(-)(1). The lithium cation positions itself preferentially over a ring. There is a small energetic preference for the 6-ring over the 5-ring binding (up to 2 kcal mol(-)(1)) and of the convex face over the concave face (3-5 kcal mol(-)(1)). The Li-bridged complexes are transition states between the pi-face complexes. Movement of the Li(+) cation over either face is facile, and the activation energy does not exceed 6 kcal mol(-)(1) on the convex face and 2.2 kcal mol(-)(1) on the concave face. In contrast, the transition of Li(+) around the corannulene edge involves a high activation barrier (24 kcal mol(-)(1) with respect to the lowest energy pi-face complex). An easier concave/convex transformation and vice versa is the bowl-to-bowl inversion with an activation energy of 7-12 kcal mol(-)(1). The computed binding energy of Li(+) to corannulene is 44 kcal mol(-)(1). Calculations of the (7)Li NMR chemical shifts and nuclear independent chemical shifts (NICS) have been performed to analyze the aromaticity of the corannulene rings and its changes upon protonation.     

Mechanisms of antibiotic resistance: QM/MM modeling of the acylation reaction of a class A beta-lactamase with benzylpenicillin

Understanding the mechanisms by which beta-lactamases destroy beta-lactam antibiotics is potentially vital in developing effective therapies to overcome bacterial antibiotic resistance. Class A beta-lactamases are the most important and common type of these enzymes. A key process in the reaction mechanism of class A beta-lactamases is the acylation of the active site serine by the antibiotic. We have modeled the complete mechanism of acylation with benzylpenicillin, using a combined quantum mechanical and molecular mechanical (QM/MM) method (B3LYP/6-31G+(d)//AM1-CHARMM22). All active site residues directly involved in the reaction, and the substrate, were treated at the QM level, with reaction energies calculated at the hybrid density functional (B3LYP/6-31+Gd) level. Structures and interactions with the protein were modeled by the AM1-CHARMM22 QM/MM approach. Alternative reaction coordinates and mechanisms have been tested by calculating a number of potential energy surfaces for each step of the acylation mechanism. The results support a mechanism in which Glu166 acts as the general base. Glu166 deprotonates an intervening conserved water molecule, which in turn activates Ser70 for nucleophilic attack on the antibiotic. This formation of the tetrahedral intermediate is calculated to have the highest barrier of the chemical steps in acylation. Subsequently, the acylenzyme is formed with Ser130 as the proton donor to the antibiotic thiazolidine ring, and Lys73 as a proton shuttle residue. The presented mechanism is both structurally and energetically consistent with experimental data. The QM/MM energy barrier (B3LYP/ 6-31G+(d)//AM1-CHARMM22) for the enzymatic reaction of 9 kcal mol(-1) is consistent with the experimental activation energy of about 12 kcal mol(-1). The effects of essential catalytic residues have been investigated by decomposition analysis. The results demonstrate the importance of the "oxyanion hole" in stabilizing the transition state and the tetrahedral intermediate. In addition, Asn132 and a number of charged residues in the active site have been identified as being central to the stabilizing effect of the enzyme. These results will be potentially useful in the development of stable beta-lactam antibiotics and for the design of new inhibitors.     

Radical-molecule reactions HCO/HOC + C2H2: mechanistic study

A detailed computational study is performed on the unknown radical-molecule reactions between HCO/HOC and acetylene (C2H2) at the CCSD(T)/6-311G(2d,p)//B3LYP/6-311G(d,p)+ZPVE, Gaussian-3//B3LYP/6-31G(d), and Gaussian-3//MP2(full)/6-31G(d) levels. For the HCO + C2H2 reaction, the most favorable pathway is direct C-addition forming the intermediate HC=CHCH=O followed by a 1,3-H-shift leading to H2C=CHC=O, which finally dissociates to the product C2H3 + CO. The overall reaction barrier is 13.8, 10.5, and 11.3 kcal/mol, respectively, at the three levels. The quasi-direct H-donation process to produce C2H3 + CO with barriers of 14.0, 14.1, and 14.1 kcal/mol is less competitive. Thus only at higher temperatures could the HCO + C2H2 reaction play a role. In contrast, the HOC + C2H2 reaction can barrierlessly generate C2H3 + CO via the quasi-direct H-donation mechanism proceeding via a prereactive complex with OH...C2 hydrogen bonding. This is suggestive of the potential importance of the HOC + C2H2 reaction in both combustion and interstellar processes. However, the direct C-addition channel is much less competitive. For both reactions, the possible formation of the intriguing interstellar molecules propadiene and propynal is also discussed. The present theoretical study represents the first attempt to probe the reaction mechanism between HOC and pi-systems. Future laboratory investigations on both reactions (particularly HOC + C2H2) are recommended.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

Ab initio model study on acetylcholinesterase catalysis: potential energy surfaces of the proton transfer reactions

Ab initio molecular orbital (MO) and hybrid density functional theory (DFT) calculations have been applied to the initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE), which is the nucleophiric addition of Ser200 in catalytic triads to a neurotransmitter acetylcholine (ACh). We focus our attention mainly on the effects of oxyanion hole and Glu327 on the potential energy surfaces (PESs) for the proton transfer reactions in the catalytic triad Ser200-His440-Glu327. The activation barrier for the addition reaction of Ser200 to ACh was calculated to be 23.4 kcal/mol at the B3LYP/6-31G(d)//HF/3-21G(d) level of theory. The barrier height under the existence of oxyanion hole, namely, Ser200-His440-Glu327-ACh-(oxyanion hole) system, decreased significantly to 14.2 kcal/mol, which is in reasonable agreement with recent experimental value (12.0 kcal/mol). Removal of Glu327 from the catalytic triad caused destabilization of both energy of transition state for the reaction and tetrahedral intermediate (product). PESs calculated for the proton transfer reactions showed that the first proton transfer process is the most important in the stabilization of tetrahedral intermediate complex. The mechanism of addition reaction of ACh was discussed on the basis of theoretical results.     

Proton and sodium cation affinities of harpagide: a computational study

The aim of this work was to estimate the proton and sodium cation affinities of harpagide (Har), an iridoid glycoside responsible for the antiinflammatory properties of the medicinal plant Harpagophytum. Monte Carlo conformational searches were performed at the semiempirical AM1 level to determine the most stable conformers for harpagide and its protonated and Na+-cationized forms. The 10 oxygen atoms of the molecule were considered as possible protonation and cationization sites. Geometry optimizations were then refined at the DFT B3LYP/6-31G level from the geometries of the most stable conformers found. Final energetics were obtained at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G level. The proton and sodium ion affinities of harpagide have been estimated at 223.5 and 66.0 kcal/mol, respectively. Since harpagide mainly provides HarNa+ ions in electrospray experiments, the DeltarG298 associated with the reaction of proton/sodium exchange between Har and methanol, MeOHNa+ + HarH+ --> MeOH2+ + HarNa+ (1), has been calculated; it has been estimated to be 1.9 kcal/mol. Complexing a methanol molecule to each reagent and product of reaction 1 makes the reaction become exothermic by 1.7 kcal/mol. These values are in the limit of the accuracy of the method and do not allow us to conclude definitely whether the reaction is endo- or exothermic, but, according to these very small values, the cation exchange reaction is expected to proceed easily in the final stages of the ion desolvation process.     

Quantum chemistry study on the mechanism of the reaction between ozone and 2,3,7,8‐TCDD

Despite of its fundamental importance, the mechanism of the reaction between ozone and dioxins are still lack detailed investigation so far. It is well-known that quantum chemical calculation is a well-established method for investigating chemical reactions. In this article, quantum chemical calculation was employed to investigate the mechanism of the reaction between ozone and dioxins, as exemplified by 2,3,7,8-TCDD. The theoretical study showed that, 2,3,7,8-TCDD was gradually destructed by ozone via six cleavages of the CC bonds. All the six cleavages of the CC bonds were calculated and discussed in detail based on the theoretical calculations by the UB3LYP/6-31G(d) method. At the same time, the energies of stationary points along the reaction process were calculated by the UMP2/6-311g(d,p)//UB3LYP/6-31G(d) method and the activation energy was obtained. The obtained activation energy was 12.25 kcal/mol, which was lower than that of the reaction between benzene and O₃(16.64 kcal/mol). This indicated that, by comparison with benzene, 2,3,7,8-TCDD could be more efficiently destructed by O₃. The reason for this result was also discussed. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

An experimental and theoretical study of the basicity of tetra-tert-butyltetrahedrane

The gas-phase basicity (GB) of tetra-tert-butyltetrahedrane (tBu4THD) was determined by FT-ICR mass spectrometry and comparison with reference compounds of known basicity. Its GB, 1035+/-10 kJ x mol(-1), makes tetra-tert-butyltetrahedrane one of the strongest bases reported so far. Ab initio calculations [B3LYP/6-31G(d) and B3LYP/6-311 + G(d,p)//6-31G(d)] have been carried out in order to compare the high experimental basicity of tBu4THD with that estimated theoretically. Both B3LYP/6-31G(d) and QCISD(T) calculations were used to determine the reaction path which connects the initial tetrahedrane-ammonium complex with the final products, protonated cyclobutadiene (CBDH+) and ammonia.     

Enthalpy of the gas-phase CO2 + Mg reaction from ab initio total energies

Various highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3(MP2), G3//B3LYP, G3(MP2)//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction enthalpies of the ground-state reaction of CO2 with Mg. All model chemistries predict highly endothermic reactions, with DeltaH(298) = 63.6-69.7 kcal x mol(-1). The difference between the calculated reaction enthalpies and the experimental value, evaluated with recommended experimental standard enthalpies of formation for products and reactants, is more than 20 kcal x mol(-1) for all methods. This difference originates in the incorrect experimental enthalpy of formation of gaseous MgO given in thermochemical databases. When the theoretical formation enthalpy for MgO calculated by a particular method is used, the deviation is reduced to 1.3 kcal x mol(-1). The performance of the methodologies used to calculate the heat of this particular reaction and the enthalpy of formation of MgO are discussed.     

Role of base assisted proton transfer in N-heterocyclic carbene-catalyzed intermolecular Stetter reaction

The mechanism of the NHC-catalyzed intermolecular Stetter reaction between benzaldehyde and cyclopropene has been investigated using the PCM-M062X/6-311++G(3df,2p)//M062X/6-31+G(d,p) level of DFT. Compared to the direct reaction, a substantial reduction in the activation free energy by 10.6–14.4 kcal/mol is observed when the reaction is performed in the presence of water, 1,5,7-triazabicyclo[4.4.0]dec-5-ene (TBD), and 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU). The bases promote the proton transfer step of the reaction to yield the Breslow intermediate. An early concerted transition state has been located for the stereocontrolling C–C bond formation step (ΔG^# = 26.6 kcal/mol) which is used to explain the diastereomeric ratio observed in the experiment.     

Accurate prediction of heats of formation by a combined method of B3LYP and neural network correction

Recently, we proposed the X1 method which combines the B3LYP/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) method with a neural network correction for an accurate yet efficient prediction of heats of formation (Wu and Xu, J Chem Phys 2007, 127, 214105). In this contribution, we discuss in detail how to set up the X1 neural network. We give examples, showing how to apply the X1 method and how the applicability of X1 can be extended. The overall mean absolute deviation of the X1 method from experiment for the 488 heats of formation is 1.52 kcal/mol compared with 9.44 kcal/mol for the original B3LYP results. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.