Phosphine adsorption on the In-rich InP(001) surface: evidence of surface dative bonds at room temperature

Adsorption of phosphine on indium phosphide compound semiconductor surfaces is a key process during the chemical vapor deposition of this material. Recent experimental infrared studies of the In-rich InP surfaces exposed to phosphine show a complex vibrational pattern in the P-H stretch region, presumably due to overlapping contributions from several structural species. We have performed density functional calculations using finite-sized cluster models to investigate the dissociative adsorption of PH3 on the In-rich InP surface. We find that initially PH3 forms a dative bond with one of the surface In atoms with a binding energy of approximately 11 kcal mol-1 at 298 K. The In-PH3 bond length is 2.9 A, 0.3 A greater than the In-P covalent bond length computed for In-PH2 species produced by hydrogen migration to a neighboring atom. However, the dissociation process, though exothermic, involves a significant activation barrier of approximately 23 kcal mol-1, suggesting the possibility of metastable trapping of the dative bonded PH3 molecules. Indeed, a careful vibrational analysis of different P-H stretching modes of the surface-bound PH3 and PH2 units gives excellent agreement with the observed infrared frequencies and their relative intensities. Moreover, at higher temperatures the frequency modes associated with PH3 disappear either due to desorption or dissociation of this molecule, an observation also well supported from the computed thermochemical parameters at different temperatures. The computed energy parameters and infrared analysis provide direct evidence that PH3 is present as a dative bonded complex on the InP surface at room temperature.     

Quantum chemical study of adsorption and dissociation of H2S on the gallium-rich GaAs (001)-4 x 2 surface

H(2)S adsorption and dissociation on the gallium-rich GaAs(001)-4 x 2 surface is investigated using hybrid density functional theory. Starting from chemisorbed H(2)S on the GaAs(001)-4 x 2 surface, two possible reaction routes have been proposed. We find that H(2)S adsorbs molecularly onto GaAs(001)-4 x 2 via the formation of a dative bond, and this process is exothermic with adsorption energy of 6.6 kcal/mol. For the first reaction route, one of the H atoms from the chemisorbed H(2)S is transferred to a second-layer As atom and the dissociated SH is inserted into the Ga-As bond with an activation barrier of 8.2 kcal/mol, which is found to be 29.3 kcal/mol more stable than the reactants. For the second case, the dissociated species may insert themselves into the Ga-Ga dimer resulting in the Ga-H-Ga and Ga-HS-Ga bridge-bonded states, which are found to be 29.8 and 22.2 kcal/mol more stable than the reactants, respectively. However, the calculations also show that the activation barrier (16.1 kcal/mol) for chemisorbed H(2)S dissociation through the second route is higher than the transfer of one H atom into a second-layer As atom. As a result, we conclude that sulfur insertion into the Ga-As bond is more kinetically favorable.     

Automerization reaction of cyclobutadiene and its barrier height: an ab initio benchmark multireference average-quadratic coupled cluster study

The problem of the double bond flipping interconversion of the two equivalent ground state structures of cyclobutadiene (CBD) is addressed at the multireference average-quadratic coupled cluster level of theory, which is capable of optimizing the structural parameters of the ground, transition, and excited states on an equal footing. The barrier height involving both the electronic and zero-point vibrational energy contributions is 6.3 kcal mol(-1), which is higher than the best earlier theoretical estimate of 4.0 kcal mol(-1). This result is confirmed by including into the reference space the orbitals of the CC sigma bonds beyond the standard pi orbital space. It places the present value into the middle of the range of the measured data (1.6-10 kcal mol(-1)). An adiabatic singlet-triplet energy gap of 7.4 kcal mol(-1) between the transition state (1)B(tg) and the first triplet (3)A(2g) state is obtained. A low barrier height for the CBD automerization and a small DeltaE((3)A(2g),(1)B(1g)) gap bear some relevance on the highly pronounced reactivity of CBD, which is briefly discussed.     

A DFT study of SiH(4) activation by Cp(2)LnH

A theoretical study of SiH(4) activation by Cp(2)LnH complexes for the entire series of lanthanides has been carried out at the DFT-B3PW91 level of theory. The reaction paths corresponding to H/H exchange and silylation, formation of Cp(2)Ln(SiH(3)), have been computed. They both occur via a single-step sigma-bond metathesis mechanism. For the athermal H/H exchange reaction, the calculated activation barrier averages 1.8 kcal.mol(-)(1) relative to the precursor adduct Cp(2)LnH(eta(2)-SiH(4)) for all lanthanide elements. The silylation path is slightly exogenic (DeltaE approximately -6.5 kcal.mol(-1)) with an activation barrier averaging 5.2 kcal.mol(-1) relative to the precursor adduct where SiH(4) is bonded by two Si-H bonds. Both pathways are therefore thermally accessible. The H/H exchange path is calculated to be kinetically more favorable whereas the silylation reaction is thermodynamically preferred. The reactivity of this familly of lanthanide complexes with SiH(4) contrasts strongly with that obtained previously with CH(4). The considerably lower activation barrier for silylation relative to methylation is attributed to the ability of Si to become hypervalent.     

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.     

Theoretical study of the mechanism and rate constant of the B + CO2 reaction

The different stationary points on the potential energy surface relative to the title reaction have been reinvestigated at the B3LYP/aug-cc-pVDZ level with relative energies computed at the CCSD(T)/aug-cc-pVTZ level with B3LYP/aug-cc-pVDZ optimized geometries and by using the G3B3 composite method. Two entrance channels have been identified. The first one corresponds to boron addition at one of the oxygen atoms of the CO 2 molecule leading to trans-BOCO, which is found to be about 27 kcal/mol exothermic with a potential energy barrier of 16.4 kcal/mol (G3B3). The second channel, which has not been identified in previous theoretical works, corresponds to a direct insertion of the boron atom into a CO bond and leads to OBCO. The B + CO 2 --> OBCO step is found to be about 84 kcal/mol exothermic and needs to overcome a potential energy barrier of only 3.6 kcal/mol (G3B3). The rate constant at 300 K of the insertion step, calculated by using TST theory with G3B3 calculated activation energy value, is 5.4 10 (-14) cm (3) molecule (-1) s (-1), in very good agreement with the experimental data ((7.0 +/- 2.8) 10 (-14) cm (3) molecule (-1) s (-1), DiGiuseppe, T. G.; Davidovits, P. J. Chem. Phys. 1981, 74, 3287). The one corresponding to the addition process is found to be several orders of magnitude smaller because of a much higher potential energy barrier. The addition channel would not contribute to the title reaction even at high temperature. A modified Arrhenius equation has been fitted in the 300-1000 K temperature range, which might be useful for chemical models.     

Energetics of cyclohexene adsorption and reaction on Pt(111) by low-temperature microcalorimetry

The heat of adsorption and sticking probability of cyclohexene on Pt(111) were measured as a function of coverage using single-crystal adsorption calorimetry in the temperature range from 100 to 300 K. At 100 K, cyclohexene adsorbs as intact di-sigma bonded cyclohexene on Pt(111), and the heat of adsorption is well described by a second-order polynomial (130 - 47 theta - 1250 theta(2)) kJ/mol, yielding a standard enthalpy of formation of di-sigma bonded cyclohexene on Pt(111) at low coverages of -135 kJ/mol and a C-Pt sigma bond strength of 205 kJ/mol. At 281 K, cyclohexene dehydrogenates upon adsorption, forming adsorbed 2-cyclohexenyl (c-C6H(9,a)) and adsorbed hydrogen, and the heat of adsorption is well described by another second-order polynomial (174 - 700 theta + 761 theta(2)) kJ/mol. This yields a standard enthalpy of formation of adsorbed 2-cyclohexenyl on Pt(111) at a low coverage of -143 kJ/mol. At coverages below 0.10 ML, the sticking probability of cyclohexene on Pt(111) is close to unity (>0.95), independent of temperature.     

New insight into the gas-phase bimolecular self-reaction of the HOO radical

The singlet and triplet potential energy surfaces (PESs) for the gas-phase bimolecular self-reaction of HOO*, a key reaction in atmospheric environments, have been investigated by means of quantum-mechanical electronic structure methods (CASSCF and CASPT2). All the reaction pathways on both PESs consist of a first step involving the barrierless formation of a prereactive doubly hydrogen-bonded complex, which is a diradical species lying about 8 kcal/mol below the energy of the reactants at 0 K. The lowest energy reaction pathway on both PESs is the degenerate double hydrogen exchange between the HOO* moieties of the prereactive complex via a double proton transfer mechanism involving an energy barrier of only 1.1 kcal/mol for the singlet and 3.3 kcal/mol for the triplet at 0 K. The single H-atom transfer between the two HOO* moieties of the prereactive complex (yielding HOOH + O2) through a pathway keeping a planar arrangement of the six atoms involves a conical intersection between either two singlet or two triplet states of A' and A" symmetries. Thus, the lowest energy reaction pathway occurs via a nonplanar cisoid transition structure with an energy barrier of 5.8 kcal/mol for the triplet and 17.5 kcal/mol for the singlet at 0 K. The simple addition between the terminal oxygen atoms of the two HOO* moieties of the prereactive complex, leading to the straight chain H2O4 intermediate on the singlet PES, involves an energy barrier of 7.3 kcal/mol at 0 K. Because the decomposition of such an intermediate into HOOH + O2 entails an energy barrier of 45.2 kcal/mol at 0 K, it is concluded that the single H-atom transfer on the triplet PES is the dominant pathway leading to HOOH + O2. Finally, the strong negative temperature dependence of the rate constant observed for this reaction is attributed to the reversible formation of the prereactive complex in the entrance channel rather than to a short-lived tetraoxide intermediate.     

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.     

Study of the dative bond in 2-aminoethoxydiphenyl borate at various levels of theory: another poor performance of the B3LYP method for B-N dative bonds

Aminoethoxydiphenyl borate (2-APB), 1, is a potent inhibitor of store-operated calcium entry channels (SOCCs). Other SOCC inhibitors are being investigated as promising pharmacological agents for a variety of conditions. Though toxic, 2-APB could be useful in the development of additional inhibitors, but its preferred binding structure must first be determined. Thus, we performed ab initio calculations to study the conformers and the strength of the dative bond of 2-APB. As a first step, we performed a series of computations at various levels of theory. We obtained vastly different dissociation energies for the dative bond depending on whether we used MP2 or B3LYP (7-10 kcal/mol different). This discrepancy has previously been observed for other B-N dative bonds by Gilbert, who found that the MP2 values were in much better agreement with experimental values (Gilbert, T. M. J. Phys. Chem. A 2004, 108, 2550-2554). Since we lacked experimental data for comparison, we performed CCSD(T) calculations and found them to have similar results to those from MP2. Thus, we conclude that MP2 is more accurate for 2-APB. The dissociation free energy at the MP2 level is 7 kcal/mol and indicates that the dative bond conformer will be the predominant structure in the gas phase. The dissociation energy is comparatively low due to the electron donation from the oxygen atom to the boron atom and due to the ring strain in the dative bond conformer.     

Ab initio study of the pathways and barriers of tricyclo[4.1.0.0(2,7)]heptene isomerization

The thermal isomerization of tricyclo[4.1.0.0(2,7)]heptene has been studied using computational chemistry with structures determined at the MCSCF level and energies at the MRMP2 level. Both the allowed conrotatory and forbidden disrotatory pathways have been elucidated resulting in cycloheptatriene isomers. Four reaction channels are available for the conrotatory pathway depending on which bond breaks first in the bicyclobutane moiety leading to enantiomeric pairs of (E,Z,Z)-1,3,5-cycloheptatriene and (Z,E,Z)-1,3,5-cycloheptatriene intermediates. The activation barrier is calculated to be 31.3 kcal·mol?1 for two channels and 37.5 kcal·mol?1 for the other two. The lower activation barrier leading to the (E,Z,Z)-1,3,5-cycloheptatriene enantiomeric pair is proposed to be due to resonance within the transition state. The same behavior was observed for the disrotatory pathway with activation barriers of 42.0 kcal·mol?1 and 55.1 kcal·mol?1 for the two channels, again with one transition state resonance stabilized. The barriers for trans double bond rotation of the intermediate cycloheptatrienes are determined to be 17.1 and 17.4 kcal·mol?1, about 5 kcal·mol?1 more than that for the seven carbon diene (E,Z)-1,3-cycloheptadiene. The electrocyclic ring closure of the trans cycloheptatrienes have been modeled and barriers determined to be 11.1 and 11.9 kcal·mol?1 for the formation of bicyclo[3.2.0]hepta-2,6-diene. This structure was previously reported as the end product for thermolysis of the parent tricyclo[4.1.0.0(2,7)]heptene. The thermodynamically more stable cycloheptatriene can be formed from bicyclo[3.2.0]hepta-2,6-diene through a two step process with a calculated pseudo first-order barrier of 36.4 kcal·mol?1. The trans-cycloheptatrienes reported herein are the first characterization of a small seven-membered ring triene with a trans double bond.     

Bonding along a linear B...N...B triad

The synthesis of tris(2,6-dihydroxyphenyl)amine diborate, 4, is reported. This compound contains a linear B...N...B array for which a symmetrical three-center two-electron (3c-2e) bond is possible. The X-ray crystal structure of 4 shows that 3c-2e bonding is, in fact, absent. Rather, the B-N-B array of 4 is unsymmetrical, having a 2c-2e B-N dative bond with the remaining boron pyramidalized outward and bonded to the oxygen of THF, i.e., 4 x THF. In THF solution, 4 displays temperature-dependent 13C NMR spectra from which a DeltaG++ of 11.6 kcal/mol at 262 K may be calculated. The dynamic process observed in solution corresponds to a bond-switching equilibrium in which the B-N bond oscillates between the two borons ("bell clapper"). Ab initio calculations indicate that the most likely pathway for the bond switch does not involve a 3c-2e B...N...B bond, but rather occurs by nucleophilic attack of THF on the datively bonded boron to generate 4 x (THF)2, lacking any B-N interactions, followed by loss of one THF. The B-N-B system of 4 sans the perturbing effect of solvent was also investigated computationally. The form of 4 containing a 3c-2e bond is found to be a transition state in the solvent-free bond-switch reaction of 4, lying 2.66 kcal/mol above 4. The stability of three-center bonds to in-line distortion (viz., X...Y...X --> X-Y.........X) is discussed from the point of view of the second-order Jahn-Teller effect.     

Understanding selectivity in the oxidative addition of the carbon-sulfur bonds of 2-cyanothiophene to Pt(0)

The reaction of 2-cyanothiophene with a zerovalent platinum bisalkylphosphine fragment yields two thiaplatinacycles derived from the cleavage of the substituted and unsubstituted C-S bonds. While cleavage away from the cyano group is preferred kinetically, cleavage adjacent to the cyano group is preferred thermodynamically. Density functional theory using B3LYP level of theory on a model of this system is consistent with experimental results in that the transition state energy leading to the formation of the kinetically favored C-S bond cleavage product is lower by 5.3 kcal mol(-1) than the barrier leading to the thermodynamically favored product. There is a 6.7 kcal mol(-1) difference between these two products. The cyano substituent at the 2- position of thiophene did not substantially change the mechanism involved in the C-S bond cleavage of thiophene previously reported.     

Effect of geminal substitution at silicon on 1-sila- and 1,3-disilacyclobutanes' strain energies, their 2+2 cycloreversion enthalpies, and Si=C pi-bond energies in silenes.

To evaluate the effect of geminal substitution at silicon on 1-sila- and 1,3-disilacyclobutanes' strain energies, their 2+2 cycloreversion enthalpies, and Si=C pi-bond energies in silenes, an ab initio MO study of silenes, R2Si=CH2 (1), 1-silacyclobutanes, cyclo-R2Si(CH2)3 (2), and 1,3-disilacyclobutanes, cyclo-(R2SiCH2)2 (3), was performed using the level of theory denoted MP4/TZ(d)//MP2/6-31G(d) (TZ means the 6-311G(d) basis set for elements of the second period and hydrogen, and the McLean-Chandler (12s,9p)/[6s,5p](d) basis set for the third period elements). In the series R = H, CH3, SiH3, CH3O, NH2, Cl, F, the growth of the reaction enthalpies and strain energies is proportional to the substituents' electronegativities. 2+2 cycloreversion of 2 is endothermic by 40.6-63.1 kcal/mol, whereas that of 3 is endothermic by 72.7-114.2 kcal/mol. On going from a silicon to a fluorine substituent at the sp2-hybridized silicon atom, the pi-bond energy in 1 weakens by 11.3 kcal/mol, and the Si=C bond length shortens by 0.053 A. The effect of substituents' electronegativities at the double-bonded silicon atom in silenes is formulated as follows: the higher electronegativity, the shorter and the weaker the Si=C pi-bond. The latter is rationalized in terms of more strained geometry resulting from the energetic cost for planarizing the R2SiC moiety. The enthalpies of the ring-opening reaction are 68.0-80.1 kcal/mol (a cleavage of the Si-C bond in 3), 65.0-76.4 kcal/mol (a cleavage of the Si-C bond in 2), and 58.0-64.9 kcal/mol (a cleavage of the C-C bond in 2). The pronounced difference in the enthalpies of 2+2 cycloreversion of 1-sila- and 1,3-disilacyclobutanes is mainly due to the difference in the enthalpies of diradicals' decomposition. The decomposition of diradicals resulting from a cleavage of C-C and Si-C bonds in 2 is exothermic by 24.3-3.3 kcal/mol (apart from the difluoro derivative which is endothermic by 5.1 kcal/mol) and 27.0-13.3 kcal/mol, respectively. The decomposition of a 1,4-diradical resulting from ring opening of 3, apart from the disilyl derivative, is the endothermic process for which the enthalpy varies from 10.6 to 40.4 kcal/mol.     

DFT study of Co-C bond cleavage in the neutral and one-electron-reduced alkyl-cobalt(III) phthalocyanines

Density functional theory (DFT) has been applied to the analysis of the structural and electronic properties of the alkyl-cobalt(III) phthalocyanine complexes, [CoIIIPc]-R (Pc = phthalocyanine, R = Me or Et), and their pyridine adducts. The BP86/6-31G(d) level of theory shows good reliability for the optimized axial bond lengths and bond dissociation energies (BDEs). The mechanism of the reductive cleavage was probed for the [CoIIIPc]-Me complex which is known as a highly effective methyl group donor. In the present analysis, which follows a recent study on the reductive Co-C bond cleavage in methylcobalamin (J. Phys. Chem. B 2007, 111, 7638-7645), it is demonstrated that addition of an electron and formation of the pi-anion radical [CoIII(Pc*)]-Me- significantly lowers the energetic barrier required for homolytic Co-C bond dissociation. Such BDE lowering in [CoIII(Pc*)]-Me- arises from the involvement of two electronic states: upon electron addition, a quasi-degenerate pi*Pc state is initially formed, but when the cobalt-carbon bond is stretched, the unpaired electron moves to a sigma*Co-C state and the final cleavage involves the three-electron (sigma)2(sigma*)1 bond. As in corrin complexes, the pi*Pc-sigma*Co-C states crossing does not take place at the equilibrium geometry of [CoIII(Pc*)]-Me- but only when the Co-C bond is stretched to approximately 2.3 A. The DFT computed Co-C BDE of 23.3 kcal/mol in the one-electron-reduced phthalocyanine species, [CoIII(Pc*)]-Me-, is lowered by approximately 37% compared to the neutral Py-[CoIIIPc]-Me complex where BDE = 36.8 kcal/mol. A similar comparison for the corrin-containing complexes shows that a DFT computed BDE of 20.4 kcal/mol for [CoIII(corrin*)]-Me leads to approximately 45% bond strength reduction, in comparison to 37.0 kcal/mol for Im-[CoIII(corrin)]-Me+. These results suggest some preference by the alkylcorrinoids for the reductive cleavage mechanism.     

The Quadruple Bonding in C2 Reproduces the Properties of the Molecule

Ever since Lewis depicted the triple bond for acetylene, triple bonding has been considered as the highest limit of multiple bonding for main elements. Here we show that C₂ is bonded by a quadruple bond that can be distinctly characterized by valence‐bond (VB) calculations. We demonstrate that the quadruply‐bonded structure determines the key observables of the molecule, and accounts by itself for about 90 % of the molecule's bond dissociation energy, and for its bond lengths and its force constant. The quadruply‐bonded structure is made of two strong π bonds, one strong σ bond and a weaker fourth σ‐type bond, the bond strength of which is estimated as 17–21 kcal mol^?1. Alternative VB structures with double bonds; either two π bonds or one π bond and one σ bond lie at 129.5 and 106.1 kcal mol^?1, respectively, above the quadruply‐bonded structure, and they collapse to the latter structure given freedom to improve their double bonding by dative σ bonding. The usefulness of the quadruply‐bonded model is underscored by “predicting” the properties of the ³ state. C₂’s very high reactivity is rooted in its fourth weak bond. Thus, carbon and first‐row main elements are open to quadruple bonding!     

Adsorption mechanisms of isoxazole and oxazole on Si(100)-2 × 1 surface: Si-N dative bond addition vs. [4+2] cycloaddition

The surface reaction pathways of isoxazole and oxazole on Si(100)-2 × 1 surface were theoretically investigated. They both form a weakly bound Si-N dative bond adduct on Si(100)-2 × 1 surface. In the case of isoxazole, the barrierlessly formed Si-N adduct is the most important surface product, that cannot be easily converted into other species. On the other hand, a facile concerted [4+2](CC) cycloaddition without involving the initial Si-N dative bond adduct was also found in the case of oxazole adsorption. The existence of Diels-Alder reactions is attributed to the particular arrangement of the two heteroatoms of oxazole in such a way that the two Si-C σ-bonds can be formed in a [4+2] fashion. In short, the unique geometric arrangements and electronegativity of these similar heteroatomic molecules yielded distinctively different surface reaction characteristics.     

Theoretical study of the Cp2Zr-catalyzed hydrosilylation of ethylene. Reaction mechanism including new sigma-bond activation

The Cp(2)Zr-catalyzed hydrosilylation of ethylene was theoretically investigated with DFT and MP2-MP4(SDQ) methods, to clarify the reaction mechanism and the characteristic features of this reaction. Although ethylene insertion into the Zr-SiH(3) bond of Cp(2)Zr(H)(SiH(3)) needs a very large activation barrier of 41.0 (42.3) kcal/mol, ethylene is easily inserted into the Zr-H bond with a very small activation barrier of 2.1 (2.8) kcal/mol, where the activation barrier and the energy of reaction calculated with the DFT(B3LYP) method are given and in parentheses are those values which have been corrected for the zero-point energy, hereafter. Not only this ethylene insertion reaction but also the coupling reaction between Cp(2)Zr(C(2)H(4)) and SiH(4) easily takes place to afford Cp(2)Zr(H)(CH(2)CH(2)SiH(3)) and Cp(2)Zr(CH(2)CH(3))(SiH(3)) with activation barriers of 0.3 (0.7) and 5.0 (5.4) kcal/mol, respectively. This coupling reaction involves a new type of Si-H sigma-bond activation which is similar to metathesis. The important interaction in the coupling reaction is the bonding overlap between the d(pi)-pi bonding orbital of Cp(2)Zr(C(2)H(4)) and the Si-H sigma orbital. The final step is neither direct C-H nor Si-C reductive elimination, because both reductive eliminations occur with a very large activation barrier and significantly large endothermicity. This is because the d orbital of Cp(2)Zr is at a high energy. On the other hand, ethylene-assisted C-H reductive elimination easily occurs with a small activation barrier, 5.0 (7.5) kcal/mol, and considerably large exothermicity, -10.6 (-7.1) kcal/mol. Also, ethylene-assisted Si-C reductive elimination and metatheses of Cp(2)Zr(H)(CH(2)CH(2)SiH(3)) and Cp(2)Zr(CH(2)CH(3))(SiH(3)) with SiH(4) take place with moderate activation barriers, 26.5 (30.7), 18.4 (20.5), and 28.3 (31.5) kcal/mol, respectively. From these results, it is clearly concluded that the most favorable catalytic cycle of the Cp(2)Zr-catalyzed hydrosilylation of ethylene consists of the coupling reaction of Cp(2)Zr(C(2)H(4)) with SiH(4) followed by the ethylene-assisted C-H reductive elimination.     

Mechanisms of initial propane activation on molybdenum oxides: a density functional theory study

We report the first detailed density functional theory study on the mechanisms of initial propane activation on molybdenum oxides. We consider 6 possible mechanisms of the C-H bond activation on metal oxides, leading to 17 transition states. We predict that hydrogen abstraction by terminal Mo=O is the most feasible reaction pathway. The calculated activation enthalpy and entropy are 32.3 kcal/mol and -28.6 cal/(mol/K), respectively, in reasonably good agreement with the corresponding experimental values (28.0 kcal/mol and -29.1 cal/(mol/K)). We find that activating the methylene C-H bond is 4.7 kcal/mol more favorable than activating the methyl C-H bond. This regioselectivity is correlated with the difference in strength between a methylene C-H bond and a methyl C-H bond. Our calculations suggest that a combined effect from both the methylene and the methyl C-H bond cleavages leads to the experimentally observed overall kinetic isotopic effects from propane to propylene on the MoO(x)/ZrO(2) catalysts.