Hetero‐π‐systems from 2 + 2 cycloreversion,part 2.1Ab initio thermochemical study of heterocyclobutanes 2 + 2 cycloreversion to form heteroethenes H2C=X (X=NH,O, SiH2, PH,S)

Ab initio and DFT thermochemical study of diradical mechanism of 2 + 2 cycloreversion of parent heterocyclobutanes and 1,3‐diheterocyclobutanes, cyclo‐(CH₂CH₂CH₂X), and cyclo‐(CH₂XCH₂X), where X = NH, O, SiH₂, PH, S, was undertaken by calculating closed‐shell singlet molecules at three levels of theory: MP4/6‐311G(d)//MP2/6‐31G(d)+ZPE, MP4/6‐311G(d,p)//MP2/6‐31G (d,p)+ZPE, and B3LYP/6‐311+G(d,p)+ZPE. The enthalpies of 2 + 2 cycloreversion decrease on going from group 14 to group 16 elements, being substantially higher for the second row elements. Normally endothermic 2 + 2 cycloreversion is predicted to be exothermic for 1,3‐diazetidine and 1,3‐dioxtane. Strain energies of the four‐membered rings were calculated via the appropriate homodesmic reactions. The enthalpies of ring opening via the every possible one‐bond homolysis that results in the formation of the corresponding 1,4‐diradical were found by subtracting the strain energies from the central bond dissociation energies of the heterobutanes CH₃CH₂—CH₂XH, CH₃CH₂—XCH₃, and HXCH₂—XCH₃. The latter energies were determined via the enthalpies of the appropriate dehydrocondensation reactions, using C—H and X—H bond energies in CH₃XH calculated at G2 level of theory. Except 1,3‐disiletane, in which ring‐opening enthalpy attains 69.7 kcal/mol, the enthalpies of the most economical ring openings do not exceed 60.7 kcal/mol. The 1,4‐diradical decomposition enthalpies found as differences between 2 + 2 cycloreversion and ring‐opening enthalpies were negative, the least exothermicity was calculated for ⋅ CH₂SiH₂CH₂CH₂. The only exception was 1,3‐disiletane, which being diradical, CH₂SiH₂CH₂SiH₂, decomposed endothermically. Since decomposition of the diradical containing two silicon atoms required extra energy, raising the enthalpy of the overall reaction to 78.9 kcal/mol, 1,3‐disiletane was predicted to be highly resisting to 2 + 2 cycloreversion. © 2007 Wiley Periodicals, Inc. Heteroatom Chem 18:704–720, 2007; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20377     

Effect of the Electronic Structure of Substituents at the Silicon Atom on the Structure and Bond Energies of Methylhydrosilanes and Silenes

Ab initiocalculations with full geometry optimization were performed for methylhydrosilanes R₂HSiCH₃, dimethylsilanes C₂Si(CH₃)₂, and silenes R₂Si = CH₂ (R = H, CH₃, SiH₃, CH₃O, NH₂, Cl, F). The enthalpies of dehydrogenation methylhydrosilanes into silenes and of dehydrocondesation of methylhydrosilanes into dimethylsilanes were calculated. The enthalpies of dehydrogenation and dehydrocondensation increase with the electronegativity of substituent R. The Si-C and Si = C bond energies were calculated. As the electronegativity of the substituent increases, the Si-C bond shortens and strengthens, while the Si = C bond shortens and weakens.     

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.     

Selective functionalization of the Si(100) surface by a bifunctional alkynylamine molecule: density functional study of the switching adsorption linkage. 2

The reaction of the bifunctional organic molecule 1-(dimethylamino)-2-propyne (DMAP) on the Si(100) surface has been investigated by density functional calculations employing a two-dimer cluster model. We found that, once in the physisorbed dative bonded well (-20.0 kcal mol(-1)), DMAP can proceed via a number of pathways, involving the formation of Si-C sigma bonds, which lead to thermodynamically more stable configurations. We first considered the cycloaddition of the CC triple bond, leading to a Si-C di-sigma bonded product (-58.7 kcal mol(-1)), for which we computed an energy barrier of only 12.5 kcal mol(-1), consistently with the observed switching of DMAP adsorption linkage at 300 K. We also explored the dissociative pathway involving the methylene C-H bond cleavage on the dative bonded DMAP, leading to three adsorption products with one (-57.3 kcal mol(-1)) and three Si-C sigma bonds (-58.7 and -60.6 kcal mol(-1)). The energy barrier for this pathway is computed 24.7 kcal mol(-1) and may therefore compete at temperature above 300 K with the reaction pathway involving the addition of the alkyne unit.     

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH(3)OO and CH(3)CH(2)OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH(3)OO(-), CD(3)OO(-), and CH(3)CH(2)OO(-)) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH(3)OO, X(2)A' '] = 1.161 +/- 0.005 eV, EA[CD(3)OO, X(2)A' '] = 1.154 +/- 0.004 eV, and EA[CH(3)CH(2)OO, X(2)A' '] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: Delta E(X(2)A' '-A (2)A')[CH(3)OO] = 0.914 +/- 0.005 eV, Delta E(X(2)A' '-A (2)A')[CD(3)OO] = 0.913 +/- 0.004 eV, and Delta E(X(2)A' '-A (2)A')[CH(3)CH(2)OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube (FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta(acid)G(298)(CH(3)OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta(acid)G(298)(CD(3)OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta(acid)G(298)(CH(3)CH(2)OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta(acid)H(298)(CH(3)OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta(acid)H(298)(CD(3)OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta(acid)H(298)(CH(3)CH(2)OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH(298)(CH(3)OO-H) = 87.8 +/- 1.0 kcal mol(-1), DH(298)(CD(3)OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH(298)(CH(3)CH(2)OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH(3)OO and CH(3)CH(2)OO. Using experimental bond enthalpies, DH(298)(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta(f)H(298)[CH(3)OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta(f)H(298)[CH(3)CH(2)OO] = -6.8 +/- 2.3 kcal mol(-1).     

Structural changes, P-P bond energies, and homolytic dissociation enthalpies of substituted diphosphines from quantum mechanical calculations

The molecular structures of the diphosphines P(2)[CH(SiH(3))(2)](4), P(2)[C(SiH(3))(3)](4), P(2)[SiH(CH(3))(2)](4), and P(2)[Si(CH(3))(3)](4) and the corresponding radicals P[CH(SiH(3))(2)](2), P[C(SiH(3))(3)](2), P[SiH(CH(3))(2)](2), and P[Si(CH(3))(3)](2) were predicted by theoretical quantum chemical calculations at the HF/3-21G*, B3LYP/3-21G*, and MP2/6-31+G* levels. The conformational analyses of all structures found the gauche conformers of the diphosphines with C(2) symmetry to be the most stable. The most stable conformers of the phosphido radicals were also found to possess C(2) symmetry. The structural changes upon dissociation allow the release of some of the energy stored in the substituents and therefore contribute to the decrease of the P-P bond dissociation energy. The P-P bond dissociation enthalpies at 298 K in the compounds studied were calculated to vary from -11.4 kJ mol(-1) (P(2)[C(SiH(3))(3)](4)) to 179.0 kJ mol(-1) (P(2)[SiH(CH(3))(2)](4)) at the B3LYP/3-21G* level. The MP2/6-31+G* calculations predict them to be in the range of 52.8-207.9 kJ mol(-1). All the values are corrected for basis set superposition error. The P-P bond energy defined by applying a mechanical analogy of the flexible substituents connected by a spring shows less variation, between 191.3 and 222.6 kJ mol(-1) at the B3LYP/3-21G level and between 225.6 and 290.4 kJ mol(-1) at the MP2/6-31+G* level. Its average value can be used to estimate bond dissociation energies from the energetics of structural relaxation.     

Thermochemistry, bond energies, and internal rotor potentials of dimethyl tetraoxide

Thermochemical properties of dimethyl tetraoxide (CH(3)OOOOCH(3)), the dimer of the methylperoxy radical, are studied using ab initio and density functional theory methods. Methylperoxy radicals are known to be important intermediates in the tropospheric ozone cycle, and the self-reaction of methylperoxy radicals, which is thought to proceed via dimethyl tetraoxide, leads to significant chain radical termination in this process. Dimethyl tetraoxide has five internal rotors, three of them unique; the potential energy profiles are calculated for these rotors, as well as for those in the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals. The dimethyl tetraoxide internal rotor profiles show barriers to rotation of 2-8 kcal mol(-1). Using B3LYP/6-31(d) geometries, frequencies, internal rotor potentials, and moments of inertia, we determine entropy and heat capacity values for dimethyl tetraoxide and its radicals. Isodesmic work reactions with the G3B3 and CBS-APNO methods are used; we calculate this enthalpy as -9.8 kcal mol(-1). Bond dissociation energies (BDEs) are calculated for all C-O and O-O bonds in dimethyl tetraoxide, again with the G3B3 and CBS-APNO theoretical methods, and we suggest the following BDEs: 46.0 kcal mol(-1) for CH(3)-OOOOCH(3), 20.0 kcal mol(-1) for CH(3)O-OOOCH(3), and 13.9 kcal mol(-1) for CH(3)OO-OOCH(3). From the BDE calculations and the isodesmic enthalpy of formation for dimethyl tetraoxide, we suggest enthalpies of 2.1, 5.8, and 1.4 kcal mol(-1) for the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals, respectively. We evaluate the suitability of 10 different density functional theory (DFT) methods for calculating thermochemical properties of dimethyl tetraoxide and its radicals with the 6-31G(d) and 6-311++G(3df,3pd) basis sets, using a variety of work reaction schemes. Overall, the best-performed DFT methods of those tested were TPSSh, BMK, and B1B95. Significant improvements in accuracy were made by moving from atomization to isodesmic work reactions, with most DFT methods yielding errors of less than 2 kcal mol(-1) with the 6-311++G(3df,3pd) basis set for isodesmic calculations on the dimethyl tetraoxide enthalpy. These isodesmic calculations were basis set consistent, with a considerable reduction in error found by using the 6-311++G(3df,3pd) basis set over the 6-31G(d) basis set. This was not the case, however, for atomization and bond dissociation work reactions, where the two basis sets returned similar results. Improved group additivity terms for the O-O-O moiety (O/O2 central atom group) are also determined.     

Methyl and silyl mesolytic dissociations in the radical cations and radical anions of but-1-ene,allylsilane, hexa-1,3-diene,and penta-2,4-dienylsilane. CAS-MCSCF and coupled cluster theoretical study

Methyl or silyl dissociation in the CH(2)=CHCH(2)-XH(3) (a-XH(3)(*)(+)) and CH(2)=CHCH=CHCH(2)-XH(3) (p-XH(3)(*) (+)) radical cations (X = C, Si) yields a(+) or p(+) and XH(3)(*). Similarly, the radical anions a-CH(3)(*) (-) and p-CH(3)(*) (-) give the pi-delocalized anion and CH(3)(*) preferentially. In contrast, a-SiH(3)(*) (-) and p-SiH(3)(*-) prefer to dissociate into the pi-delocalized radical and silide. All reactions are endoergic: by 43-50 kcal mol(-)(1) in the radical cations, and easier to some extent in the radical anions, that require 29-33 (X = C) and 13-14 kcal mol(-)(1) (X = Si). The fragmentation energy profiles do not present significant barriers for the backward process in the case of the radical cations. All radical anions exhibit an energy maximum along the dissociation pathway, but the barrier is lower than the dissociation limit. Fragmentation is "activated" more in the anions than in the cations with respect to homolysis in the corresponding neutrals (that requires 72-81 kcal mol(-)(1)). Wave function analysis indicates that the C-X bond cleavage in the hydrocarbon radical ions, although formally comparable to a homolytic process, is at variance with this model, due to the spin recoupling of one of the two C-X bond electrons with the originally unpaired electron. This is basically true also for the silyl-substituted radical anions, in which the initial more delocalized charge distribution might suggest some heterolytic character of the bond cleavage.     

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.     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

Theoretical study on structures and stability of Si2CP isomers

The structures, energetics, spectroscopies, and isomerization of various doublet Si2CP species are explored theoretically. In contrast to the previously studied SiC2N and SiC2P radicals that have linear SiCCN and SiCCP ground states, the title Si2CP radical has a four-membered-ring form cSiSiPC 1 (0.0 kcal/mol) with Si-C cross-bonding as the ground-state isomer at the CCSD(T)/6-311G(2df)//B3LYP/6-311G(d)+ZPVE level, similar to the Si2CN radical. The second low-lying isomer 2 at 11.6 kcal/mol has a SiCSiP four-membered ring with C-P cross-bonding, yet it is kinetically quite unstable toward conversion to 1 with a barrier of 3.5 kcal/mol. In addition, three cyclic species with divalent carbene character, i.e., cSiSiCP 7, 7' with C-P cross-bonding and cSiCSiP 8 with Si-Si cross-bonding, are found to possess considerable kinetic stability, although they are energetically high lying at 44.4, 46.5, and 41.4 kcal/mol, respectively. Moreover, a linear isomer SiCSiP 5 at 44.3 kcal/mol also has considerable kinetic stability and predominantly features the interesting cumulenic /Si=C=Si=P/* form with a slight contribution from the silicon-phosphorus triply bonded form /Si=C*-Si[triple bond]P/. The silicon-carbon triply bonded form *Si[triple bond]C-Si[triple bond]P/ has negligible contribution. All five isomers are expected to be observable in low-temperature environments. Their bonding nature and possible formation strategies are discussed. For relevant species, the QCISD/6-311G(d) and CCSD(T)/6-311+G(2df) (single-point) calculations are performed to provide more reliable results. The calculated results are compared to those of the analogous C3N, C3P, SiC2N, and Si2CN radicals with 17 valence electrons. Implications in interstellar space and P-doped SiC vaporization processes are also discussed.     

Strain energy of small ring hydrocarbons. Influence of C-h bond dissociation energies

Ab initio calculations at the G2, G3, and CBS-Q levels of theory have been applied to the question of the origin of ring strain in a series of unsaturated hydrocarbons. In addition to the angular ring strain germane to all three-membered ring hydrocarbons, a general trend is in evidence that suggests that the increased ring strain (SE) of unsaturated small ring alkenes may be attributed in part to their relatively weak allylic C-H bonds. The high strain energy of cyclopropene (54.1 kcal/ mol) is attributed largely to angular strain. The anomalously low SE of cyclobutene relative to cyclobutane (DeltaSE = 4 kcal/mol) is a consequence of normal C-H bond dissociation energies for cyclobutane (100.6 kcal/mol) and very strong vinyl C-H bonds (111.9 kcal/mol) and a relatively strong pi-bond energy (63.5 kcal/mol) for cyclobutene. The greater SE of methylenecyclopropane (39.5 kcal/ mol), relative to methylcyclopropane (29.8 kcal/mol), can be attributed to the strong ring C-H bonds of methylcyclopropane (110.5 kcal/mol) and relatively weak allylic C-H bonds (99.3 kcal/mol) of methylenecyclopropane. The increased SE of 1-methylcyclopropene relative to isomeric methylenecyclopropane is ascribed to its weak ring C-H bonds and to angular strain. The relative thermodynamic stability of a series of small ring alkenes is determined by a measure of their hydrogenation enthalpies. Independent confirmation of the SEs of a series of substituted cyclopropenes is provided by their dimerization/combination with cyclopropane to form a six-membered ring reference compound.     

Energetics of the O-H bond and of intramolecular hydrogen bonding in HOC6H4C(O)Y (Y = H, CH3, CH2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds

The energetics of the phenolic O-H bond in a series of 2- and 4-HOC 6H 4C(O)Y (Y = H, CH3, CH 2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH 3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds and of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y, was investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of 2-hydroxybenzaldehyde (2HBA), 4-hydroxybenzaldehyde (4HBA), 2'-hydroxyacetophenone (2HAP), 2-hydroxybenzamide (2HBM), and 4-hydroxybenzamide (4HBM), at 298.15 K, were determined by micro- or macrocombustion calorimetry. The corresponding enthalpies of vaporization or sublimation were also measured by Calvet drop-calorimetry and Knudsen effusion measurements. The combination of the obtained experimental data led to Delta f H m (o)(2HBA, g) = -238.3 +/- 2.5 kJ.mol (-1), DeltafHm(o)(4HBA, g) = -220.3 +/- 2.0 kJ.mol(-1), Delta f H m (o)(2HAP, g) = -291.8 +/- 2.1 kJ.mol(-1), DeltafHm(o)(2HBM, g) = -304.8 +/- 1.5 kJ.mol (-1), and DeltafHm(o) (4HBM, g) = -278.4 +/- 2.4 kJ.mol (-1). These values, were used to assess the predictions of the B3LYP/6-31G(d,p), B3LYP/6-311+G(d,p), B3LYP/aug-cc-pVDZ, B3P86/6-31G(d,p), B3P86/6-311+G(d,p), B3P86/aug-cc-pVDZ, and CBS-QB3 methods, for the enthalpies of a series of isodesmic gas phase reactions. In general, the CBS-QB3 method was able to reproduce the experimental enthalpies of reaction within their uncertainties. The B3LYP/6-311+G(d,p) method, with a slightly poorer accuracy than the CBS-QB3 approach, achieved the best performance of the tested DFT models. It was further used to analyze the trends of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y evaluated by the ortho-para method and to compare the energetics of the phenolic O-H bond in 2- and 4-HOC 6H 4C(O)Y compounds. It was concluded that the O-H bond "strength" is systematically larger for 2-hydroxybenzoyl than for the corresponding 4-hydroxybenzoyl isomers mainly due to the presence of the intramolecular O...H hydrogen bond in the 2-isomers. The observed differences are, however, significantly dependent on the nature of the substituent Y, in particular, when an intramolecular H-bond can be present in the radical obtained upon cleavage of the O-H bond.     

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.     

New selective haloform-type reaction yielding 3-hydroxy-2,2-difluoroacids: theoretical study of the mechanism

Experimental results of an unprecedented haloform-type reaction in which 4-alkyl-4-hydroxy-3,3-difluoromethyl trifluoromethyl ketones undergo base-promoted selective cleavage of the CO-CF(3) bond, yielding 3-hydroxy-2,2-difluoroacids and fluoroform, are rationalized using DFT (B3LYP) calculations. The gas-phase addition of hydroxide ion to 1,1,1,3,3-pentafluoro-4-hydroxypentan-2-one (R) is found to be a barrierless process, yielding a tetrahedral intermediate (INT), involving a DeltaG(r)(298 K) of -61.4 kcal/mol. The CO-CF(3) bond cleavage in INT leads to a hydrogen-bonded [CH(3)CHOHCF(2)CO(2)H...CF(3)](-) complex by passage through a transition structure (TS1) with a DeltaG()(298 K) of 20.8 kcal/mol and a DeltaG(r)(298 K) of 9.8 kcal/mol. This complex undergoes a proton transfer between its components, yielding a hydrogen-bonded [CH(3)CHOHCF(2)CO(2)...CHF(3)](-) complex. This process has associated with it a DeltaG()(298 K) of only 3.1 kcal/mol and a DeltaG(r)(298 K) of -43.3 kcal/mol. The CO-CF(2) bond cleavage in INT leads to a hydrogen-bonded [CH(3)CHOHCF(2)...CF(3)CO(2)H](-) complex by passage through a transition structure (TS3) with a DeltaG()(298 K) of 29.2 kcal/mol and a DeltaG(r)(298 K) of 25.1 kcal/mol. The lower energy barrier found for CO-CF(3) bond cleavage in INT is ascribed to the larger number of fluorine atoms stabilizing the negative charge accumulated on the CF(3) moiety of TS1, as compared to the number of fluorine atoms stabilizing the negative charge on the CH(3)CHOHCF(2) moiety of TS3. The solvent-induced effects on the two pathways, introduced within the SCRF formalism through PCM calculations, do not reverse the predicted preference of the CO-CF(3) over the CO-CF(2) bond cleavage of R in the gas phase.     

Reactivity of alkyl versus silyl peroxides. The consequences of 1, 2-silicon bridging on the epoxidation of alkenes with silyl hydroperoxides and bis(trialkylsilyl)peroxides

The bond dissociation energies for a series of silyl peroxides have been calculated at the G2 and CBS-Q levels of theory. A comparison is made with the O-O BDE of the corresponding dialkyl peroxides, and the effect of the O-O bond strength on the activation barrier for oxygen atom transfer is discussed. The O-O bond dissociation enthalpies (DeltaH(298)) for bis (trimethylsilyl) peroxide (1) and trimethylsilyl hydroperoxide (2) are 54.8 and 53.1 kcal/mol, respectively at the G2 (MP2) and CBS-Q levels of theory. The O-O bond dissociation energies computed at G2 and G2(MP2) levels for bis(tert-butyl) peroxide and tert-butyl hydroperoxide are 45.2 and 48.3 kcal/mol, respectively. The barrier height for 1,2-methyl migration from silicon to oxygen in trimethylsilyl hydroperoxide is 47.9 kcal/mol (MP4//MP2/6-31G). The activation energy for the oxidation of trimethylamine to its N-oxide by bis(trimethylsilyl) peroxide is 28.2 kcal/mol (B3LYP/6-311+G(3df,2p)// B3LYP/6-31G(d)). 1,2-Silicon bridging in the transition state for oxygen atom transfer to a nucleophilic amine results in a significant reduction in the barrier height. The barrier for the epoxidation of E-2-butene with bis(dimethyl(trifluoromethyl))silyl peroxide is 25.8 kcal/mol; a reduction of 7.5 kcal/mol relative to epoxidation with 1. The activation energy calculated for the epoxidation of E-2-butene with F(3)SiOOSiF(3) is reduced to only 2.2 kcal/mol reflecting the inductive effect of the electronegative fluorine atoms.     

Theoretical study of (CH...C)- hydrogen bonds in CH(4-n)X(n) (X = F, Cl; n = 0, 1, 2) systems complexed with their homoconjugate and heteroconjugate carbanions

This work deals with a theoretical study of the (CH...C)- hydrogen bonds in CH4, CH3X, and CH2X2 (X = F, Cl) complexed with their homoconjugate and heteroconjugate carbanions. The properties of the complexes are calculated with the B3LYP method using the 6-311++G(d,p) or 6-311++G(2df,2p) basis sets. The deprotonation enthalpies (DPE) of the CH bond or the proton affinities of the carbanions (PA(C-) are calculated as well. All the systems with the exception of the CH4...CHCl2(-) one are characterized by a double minimum potential. In some of the complexes, the (CH(b)...C)- hydrogen bond is linear. In other systems, such as CH3F...CH2F- and CH3F...CHF2(-), there is a large departure from linearity, the systems being stabilized by electrostatic interactions between the nonbonded H of the neutral molecule and the F atom of the carbanion. In the transition state, the (CH(b)...C)- bond is linear, and there is a large contraction of the intermolecular C...C distance. The binding energies vary within a large range, from -1.4 to -11.1 kcal mol(-1) for the stable complexes and -8.6 to -44.1 kcal mol(-1) for the metastable complexes. The energy barriers to proton transfer are between 5 and 20 kcal mol(-1) for the heteroconjugate systems and between 3.8 and 8.3 kcal mol(-1) for the homoconjugate systems. The binding energies of the linear complexes depend exponentially on 1.5DPE - PA(C-), showing that the proton donor is more important than the proton acceptor in determining hydrogen bond strength. The NBO analysis indicates an important electronic reorganization in the two partners. The elongations of the CH bond resulting from the interaction with the carbanion depend on the occupation of the sigma*(CH(b)) antibonding orbitals and on the hybridization of the C bonded to H(b). The frequency shifts of the nu(CH)(A1) stretching vibration range between 15 and 1150 cm(-1). They are linearly correlated to the elongation of the CH(b) bond.     

1,3 Geminal interactions as the possible trend setting factors for C-H and C-C bond energies in alkanes. Support from a density functional theory based bond energy decomposition study

A bond energy decomposition analysis has been carried out to rationalize the well-established experimental fact that C-C and C-H bond energies decrease with increasing substitution on the carbon. It is shown that this trend is set by steric 1,3 repulsive interactions (geminal repulsion) that increase in the order 1,3 hydrogen-hydrogen < 1,3 hydrogen-carbon < 1,3 carbon-carbon. On the other hand, the radical stabilization energy has little influence on the observed trend for the C-H bond energy in H-CR(1)R(2)R(3) or the C-C bond energy in H(3)C-CR(1)R(2)R(3). Thus, it varies in H-CR(1)R(2)R(3) from -7.2 kcal/mol (H-CH3) to -6.5 kcal/mol (H-C(CH3)3) and in H(3)C-CR(1)R(2)R(3) from -19.0 kcal/mol (H(3)C-CH(3)) to -16.9 kcal/mol (H(3)C-C(CH(3))(3)). It was further found that the average intrinsic C-H bond energy in H-CR(1)R(2)R(3) of 129.2 kcal/mol is smaller than the average intrinsic C-C bond energy in H(3)C-CR(1)R(2)R(3) of 143.4 kcal/mol. However, after the inclusion of steric effects, the overall C-H bond becomes stronger than the C-C bond. The role of steric 1,3 repulsive interactions as the trend setting factor has most recently been suggested by Gronert (J. Org. Chem. 2006, 71, 1209) based on an empirical fit of alkane atomization energies.     

Why does perfluorination render bicyclo[2.2.0]hex-1(4)-ene stable toward dimerization? Calculations provide the answers

B3LYP calculations with two different basis sets have been performed to understand why bicyclo[2.2.0]hex-1(4)-ene (1a) undergoes dimerization with DeltaH(++) = 11.5 kcal/mol, but dimerization of perfluorobicyclo[2.2.0]hex-1(4)-ene (1b) has never been observed. The former reaction is computed to be exothermic by 37.2 kcal/mol, whereas the latter is calculated to be endothermic by 7.4 kcal/mol. The 44.6 kcal/mol difference between the enthalpies of these two reactions can be dissected into contributions of 24.5 kcal/mol for the difference between the enthalpies for forming diradical intermediates 2a and 2b and 20.1 kcal/mol for cyclization of 2a and 2b to, respectively, 3a and 3b. The latter enthalpy difference is largely attributable to repulsions between the endo-fluorines in the dimer, although the exo-fluorines also are found to contribute. The former enthalpy difference is attributable to the difference between the dissociation enthalpies of the pi bonds in 1a and 1b, which is shown to amount to 16 +/- 1 kcal/mol. About 25% of the stronger pi bond in fluoroalkene 1b is found to be due to hyperconjugation of the eight C-F bonds in 1b with the filled pi orbital. However, the major contributor to the stronger pi bond in 1b is shown to be the unfavorable interaction that results when a pyramidalized radical center is syn to a C-F bond. Both of these effects, which contribute to the greater strength of the pi bond in 1b, relative to that in 1a, are analyzed and discussed.     

Enthalpies of formation, bond dissociation energies, and molecular structures of the n-aldehydes (acetaldehyde, propanal, butanal, pentanal, hexanal, and heptanal) and their radicals

Aldehydes are important intermediates and products in a variety of combustion and gas-phase oxidation processes, such as in low-temperature combustion, in the atmosphere, and in interstellar media. Despite their importance, the enthalpies of formation and bond dissociation energies (BDEs) for the aldehydes are not accurately known. We have determined enthalpies of formation for acetaldehyde, propanal, and butanal from thermodynamic cycles, using experimentally measured reaction and formation enthalpies. All enthalpy values used for reference molecules and reactions were first verified to be accurate to within around 1 kcal mol-1 using high-level ab initio calculations. Enthalpies of formation were found to be -39.72 +/- 0.16 kcal mol-1 for acetaldehyde, -45.18 +/- 1.1 kcal mol-1 for propanal, and -49.27 +/- 0.16 kcal mol-1 for butanal. Enthalpies of formation for these three aldehydes, as well as for pentanal, hexanal, and heptanal, were calculated using the G3, G3B3, and CBS-APNO theoretical methods, in conjunction with bond-isodesmic work reactions. On the basis of the results of our thermodynamic cycles, theoretical calculations using isodesmic work reactions, and existing experimental measurements, we suggest that the best available formation enthalpies for the aldehydes acetaldehyde, propanal, butanal, pentanal, hexanal, and heptanal are -39.72, -45.18, -50.0, -54.61, -59.37, and -64.2 kcal mol-1, respectively. Our calculations also identify that the literature enthalpy of formation of crotonaldehyde is in error by as much as 1 kcal mol-1, and we suggest a value of -25.1 kcal mol-1, which we calculate using isodesmic work reactions. Bond energies for each of the bonds in the aldehydes up to pentanal were calculated at the CBS-APNO level. Analysis of the BDEs reveals the R-CH(2)CH=O to be the weakest bond in all aldehydes larger than acetaldehyde, due to formation of the resonantly stabilized vinoxy radical (vinyloxy radical/formyl methyl radical). It is proposed that the vinoxy radical as well as the more commonly considered formyl and acetyl radicals are important products of aldehyde combustion and oxidation, and the reaction pathways of the vinoxy, formyl, and acetyl radicals are discussed. Group additivity values for the carbon-oxygen-hydrogen groups common to the aldehydes are also determined. Internal rotor profiles and electrostatic potential surfaces are used to study the dipole induced dipole-dipole interaction in the synperiplanar conformation of propanal. It is proposed that the loss of this dipole-dipole interaction in RC(.-)HCH(2)CH=O radicals causes a ca. 1-2 kcal mol-1 decrease in the aldehyde C-H and C-C bond energies corresponding to RC(.-)HCH(2)CH=O radical formation.