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.     

A novel salt bridge mechanism highlights the need for nonmobile proton conditions to promote disulfide bond cleavage in protonated peptides under low-energy collisional activation

The gas-phase fragmentation mechanisms of small models for peptides containing intermolecular disulfide links have been studied using a combination of tandem mass spectrometry experiments, isotopic labeling, structural labeling, accurate mass measurements of product ions, and theoretical calculations (at the MP2/6-311 + G(2d,p)//B3LYP/3-21G(d) level of theory). Cystine and its C-terminal derivatives were observed to fragment via a range of pathways, including loss of neutral molecules, amide bond cleavage, and S-S and C-S bond cleavages. Various mechanisms were considered to rationalize S-S and C-S bond cleavage processes, including charge directed neighboring group processes and nonmobile proton salt bridge mechanism. Three low-energy fragmentation pathways were identified from theoretical calculations on cystine N-methyl amide: (1) S-S bond cleavage dominated by a neighboring group process involving the C-terminal amide N to form either a protonated cysteine derivative or protonated sulfenyl amide product ion (44.3 kcal mol(-1)); (2) C-S bond cleavage via a salt bridge mechanism, involving abstraction of the alpha-hydrogen by the N-terminal amino group to form a protonated thiocysteine derivative (35.0 kcal mol(-1)); and (3) C-S bond cleavage via a Grob-like fragmentation process in which the nucleophilic N-terminal amino group forms a protonated dithiazolidine (57.9 kcal mol(-1)). Interestingly, C-S bond cleavage by neighboring group processes have high activation barriers (63.1 kcal mol(-1)) and are thus not expected to be accessible during low-energy CID experiments. In comparison to the energetics of simple amide bond cleavage, these S-S and C-S bond cleavage reactions are higher in energy, which helps rationalize why bond cleavage processes involving the disulfide bond are rarely observed for low-energy CID of peptides with mobile proton(s) containing intermolecular disulfide bonds. On the other hand, the absence of a mobile proton appears to "switch on" disulfide bond cleavage reactions, which can be rationalized by the salt bridge mechanism. This potentially has important ramifications in explaining the prevalence of disulfide bond cleavage in singly protonated peptides under MALDI conditions.     

Reactions of C5H5M (M = Fe, Ni) with substituted thiophenes

The ion molecule reactions between C₅H₅M⁺ (M = Fe, Ni) with some substituted thiophenes have been studied in an ion trap mass spectrometer. The reactions of halogen substituted thiophenes lead to the formation of a new C-C bond between the cyclopentadiene ring and the thiophene with the loss of a neutral HX. The reaction mechanism has been investigated by means of DFT calculations and it was found that the insertion of the metal atom in the C-X bond is the key step in the process.     

Three-fold barrier and normal coordinate analyses of (trihalomethyl)sulfenyl halides CX3-SX (X = F and Cl)

The structural stability of (trihalomethyl)sulfenyl halides CX3-SX (X is F and Cl) was investigated by DFT-B3LYP and ab initio MP2 calculations using 6-311 + G** basis set. Full energy optimizations were carried out from which the three-fold barrier about C-S bond was calculated to be about 3 kcal mol(-1) in (trifluoromethyl)sulfenyl fluoride and (trifluoromethyl)sulfenyl chloride and about 6 kcal mol(-1) in (trichloromethyl)sulfenyl fluoride and (trichloromethyl)sulfenyl chloride. The vibrational frequencies of the four molecules were computed at the DFT-B3LYP level and the vibrational assignments for the normal modes of the compounds in their ground state structure were made on the basis of normal coordinate calculations and reported experimental data.     

A density functional theory study on the mechanism of the permanganate oxidation of substituted alkenes.

Oxidation of substituted alkenes by permanganate follows a (3 + 2) cycloaddition mechanism. DFT calculations (Becke3LYP/6-31G(d)) strongly favor the (3 + 2) pathway against a (2 + 2) pathway that proceeds through a metallaoxetane. The difference in free activation energy between the two pathways is around 40 to 45 kcal/mol for the nine compounds covered. The results for the (3 + 2) cycloaddition mechanism are in agreement with experimentally observed kinetic data reported earlier. Symmetric transition states are calculated for alkenes with at least one CH2 group between the double bond and the acid group, while all others are unsymmetrical due to the repulsion between the permanganate oxygens and the carboxylic oxygens. Steric bulk, introduced by the number and position of substituents at the double bond, has no significant effect on the activation energies. A higher level of theory (Becke3LYP/6-311G(d,p)) leads to a reduction of the rmsd between experimental and calculated values, but has little influence on the calculated geometries.     

Addition of alkyl radicals to transition-metal-coordinated CO: calculation of the reaction of [Ru(CO)5] and related complexes and relevance to alkane carbonylation

Electronic structure methods have been used to study the transition state and products of the reaction between alkyl radicals and CO coordinated in transition-metal complexes. At the B3LYP DFT level, methyl addition to a carbonyl of [Ru(CO)5] or [Ru(CO)3(dmpe)] is calculated to be about 6 kcal/mol more exothermic than addition to free CO. In contrast, methyl addition to [Mo(CO)6] is 12 kcal/mol less exothermic than addition to CO, while the reaction enthalpy of methyl addition to [Pd(CO)4] is comparable to that of free CO. Related results are obtained at the CCSD-T level and for the reactions of the cyclohexyl radical. The transition state for these reactions is characterized by significant distortion of the geometry of the reactant complex, which include lengthening and bending of the M-CO bond, but with negligible C-C bond formation. Accordingly, the activation energy for addition to coordinated carbonyls is 2-10 kcal/mol greater than that of addition to free CO. Additional calculations were also carried out on representative unsaturated metal carbonyls. The calculated results afford an understanding of the mechanism of previously reported photochemical alkane carbonylation systems utilizing d(8)-ML5 metal carbonyls as cocatalysts. In particular, it is strongly indicated that such systems operate via direct attack by an alkyl radical at a CO ligand, a reaction that has not previously been proposed.     

On the nature of carbon-hydrogen bond activation at rhodium and related reactions

Over the past 20 years, substantial progress has been made in the understanding of the activation of C-H and other strong bonds by reactive metal complexes in low oxidation states. This paper will present an overview of the use of pentamethylcyclopentadienyl and trispyrazolylborate rhodium complexes for the activation of arene and alkane C-H bonds. Insights into bond strengths, kinetic and thermodynamic selectivities, and the nature of the intermediates involved will be reviewed. The role of eta-2 arene complexes will be shown to be critical to the C-H activation reactions. Some information about the fleeting alkane sigma-complexes will also be presented. In addition, use of these complexes with thiophenes has shown the ability to cleave C-S bonds. Mechanistic information has been obtained indicating coordination through sulfur prior to cleavage. Relevant examples of nickel-based C-S cleavage will also be given.     

Solution calorimetric and stopped-flow kinetic studies of the reaction of *Cr(CO)3C5Me5 with PhSe-SePh and PhTe-TePh. Experimental and theoretical estimates of the Se-Se, Te-Te, H-Se, and H-Te bond strengths

The kinetics of the oxidative addition of PhSeSePh and PhTeTePh to the stable 17-electron complex *Cr(CO)3C5Me5 have been studied utilizing stopped-flow techniques. The rates of reaction are first-order in each reactant, and the enthalpy of activation decreases in going from Se (deltaH(double dagger) = 7.0 +/- 0.5 kcal/mol, deltaS(double dagger) = -22 +/- 3 eu) to Te (deltaH(double dagger) = 4.0 +/- 0.5 kcal/mol, deltaS(double dagger) = -26 +/- 3 eu). The kinetics of the oxidative addition of PhSeH and *Cr(CO)3C5Me5 show a change in mechanism in going from low (overall third-order) to high (overall second-order) temperatures. The enthalpies of the oxidative addition of PhE-EPh to *Cr(CO)3C5Me5 in toluene solution have been measured and found to be -29.6, -30.8, and -28.9 kcal/mol for S, Se, and Te, respectively. These data are combined with enthalpies of activation from kinetic studies to yield estimates for the solution-phase PhE-EPh bond strengths of 46, 41, and 33 kcal/mol for E = S, Se, and Te, respectively. The corresponding Cr-EPh bond strengths are 38, 36, and 31 kcal/mol. Two methods have been used to determine the enthalpy of hydrogenation of PhSeSePh in toluene on the basis of reactions of HSPh and HSePh with either *Cr(CO)3C5Me5 or 2-pyridine thione. These data lead to a thermochemical estimate of 72 kcal/mol for the PhSe-H bond strength in toluene solution, which is in good agreement with kinetic studies of H atom transfer from HSePh at higher temperatures. The reaction of H-Cr(CO)3C5Me5 with PhSe-SePh is accelerated by the addition of a Cr radical and occurs via a rapid radical chain reaction. In contrast, the reaction of PhTe-TePh and H-Cr(CO)3C5Me5 does not occur at any appreciable rate at room temperature, even in the presence of added Cr radicals. This is in keeping with a low PhTe-H bond strength blocking the chain and implies that H-TePh < or = 63 kcal/mol. Structural data are reported for PhSe-Cr(CO)3C5Me5 and PhS-Cr(CO)3C5Me5. The two isostructural complexes do not show signs of an increase in steric strain in terms of metal-ligand bonds or angles as the Cr-EPh bond is shortened in going from Se to S. Bond strength estimates of the PhE-H and PhE-EPh derived from density functional theory calculations are in reasonable agreement with experimental data for E = Se but not for E = Te. The nature of the singly occupied molecular orbital of the *EPh radicals is calculated to show increasing localization on the chalcogenide atom in going from S to Se to Te.     

Activation of C-H / H-H bonds by rhodium(II) porphyrin bimetalloradicals

Reactivity, kinetic, and thermodynamic studies are reported for reactions of a rhodium(II) bimetalloradical with H(2), and with the methyl C-H bonds for a series of substrates CH(3)R (R = H, CH(3), OH, C(6)H(5)) using a m-xylyl diether tethered diporphyrin ligand. Bimolecular substrate reactions involving the intramolecular use of two metalloradical centers and preorganization of the four-centered transition state (M*...X...Y*...M) result in large rate enhancements as compared to termolecular reactions of monometalloradicals. Activation parameters and deuterium kinetic isotope effects for the substrate reactions are reported. The C-H bond reactions become less thermodynamically favorable as the substrate steric requirements increase, and the activation free energy (DeltaG++) decreases regularly as DeltaG degrees becomes more favorable. An absolute Rh-H bond dissociation enthalpy of 61.1 +/- 0.4 kcal mol(-1) is directly determined, and the derived Rh-CH(2)R BDE values increase regularly with the increase in the C-H BDE.     

DFT studies on the copolymerization of alpha-olefins with polar monomers: ethylene-methyl acrylate copolymerization catalyzed by a Pd-based diimine catalyst.

Gradient-corrected density functional theory has been used to study the elementary reactions for the copolymerization of ethylene with methyl acrylate catalyzed by Pd-based diimine catalysts, modeled by the generic complex N(wedge)N-Pd(n-C(3)H(7))(+), with N(wedge)N = -NHCHCHNH-. The steric effects in the real systems are discussed on the basis of the calculations for the catalyst with N(wedge)N = -NArCRCRNAr-, R = CH(3), and Ar = C(6)H(3)(i-Pr(2)) and the previous calculations on ethylene/propylene polymerization. Considerations have been given to the different possible acrylate complexes, as well as the transition states and the products (agostic complexes and the alternative chelates) for two acrylate insertion paths (1,2 and 2,1). The chelate-opening reactions have also been studied. The results revealed a strong electronic preference for the 2,1-insertion paths, with a barrier that is 4.5 kcal/mol lower than any other studied insertion pathway. In the real systems the 2,1-insertion of acrylate is preferred by 0.5 kcal/mol. The 2,1-insertion barrier calculated for the real system of 12.4 kcal/mol is in very good agreement with the experimental value of 12.1 kcal/mol. The six-member chelate is the most stable insertion product with an energy that is 21 kcal/mol lower than the kinetic insertion product. The reactions of the chelate opening by ethylene that start from the lowest energy complexes have the lowest barrier for the four-member ring (23 kcal/mol) and the highest for the six-member structure (30.4 kcal/mol). The high barrier for the opening of the six-member chelate suggests the possibility of a two-step chelate-opening mechanism. The internal barriers for the chelate-opening reactions starting from the higher energy complexes are lower then the one-step reaction that starts from the preferred complex and comparable to those of the ethylene insertion into the Pd-alkyl bond. While the chelate opening by a subsequent acrylate insertion seems to be facile for the generic catalyst, steric effects in the real catalyst are likely to decrease the acrylate pi-complexation energies and increase the insertion barriers to the extent where such a reaction becomes unfeasible.     

Oxidative addition of the fluoromethane C-F bond to Pd. An ab initio benchmark and DFT validation study

We have computed a state-of-the-art benchmark potential energy surface (PES) for two reaction pathways (oxidative insertion, OxIn, and S(N)2) for oxidative addition of the fluoromethane C-F bond to the palladium atom and have used this to evaluate the performance of 26 popular density functionals, covering LDA, GGA, meta-GGA, and hybrid density functionals, for describing these reactions. The ab initio benchmark is obtained by exploring the PES using a hierarchical series of ab initio methods (HF, MP2, CCSD, CCSD(T)) in combination with a hierarchical series of seven Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account through a full four-component all-electron approach. Our best estimate of kinetic and thermodynamic parameters is -5.3 (-6.1) kcal/mol for the formation of the reactant complex, 27.8 (25.4) kcal/mol for the activation energy for oxidative insertion (OxIn) relative to the separate reactants, 37.5 (31.8) kcal/mol for the activation energy for the alternative S(N)2 pathway, and -6.4 (-7.8) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). Our work highlights the importance of sufficient higher angular momentum polarization functions for correctly describing metal-d-electron correlation. Best overall agreement with our ab initio benchmark is obtained by functionals from all three categories, GGA, meta-GGA, and hybrid DFT, with mean absolute errors of 1.4-2.7 kcal/mol and errors in activation energies ranging from 0.3 to 2.8 kcal/mol. The B3LYP functional compares very well with a slight underestimation of the overall barrier for OxIn by -0.9 kcal/mol. For comparison, the well-known BLYP functional underestimates the overall barrier by -10.1 kcal/mol. The relative performance of these two functionals is inverted with respect to previous findings for the insertion of Pd into the C-H and C-C bonds. However, all major functionals yield correct trends and qualitative features of the PES, in particular, a clear preference for the OxIn over the alternative S(N)2 pathway.     

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.     

Mechanism of homogeneous Ir(III) catalyzed regioselective arylation of olefins

The mechanism of hydroarylation of olefins by a homogeneous Ph-Ir(acac)(2)(L) catalyst is elucidated by first principles quantum mechanical methods (DFT), with particular emphasis on activation of the catalyst, catalytic cycle, and interpretation of experimental observations. On the basis of this mechanism, we suggest new catalysts expected to have improved activity. Initiation of the catalyst from the inert trans-form into the active cis-form occurs through a dissociative pathway with a calculated DeltaH(0 K)() = 35.1 kcal/mol and DeltaG(298 K)() = 26.1 kcal/mol. The catalytic cycle features two key steps, 1,2-olefin insertion and C-H activation via a novel mechanism, oxidative hydrogen migration. The olefin insertion is found to be rate determining, with a calculated DeltaH(0 K)() = 27.0 kcal/mol and DeltaG(298 K)() = 29.3 kcal/mol. The activation energy increases with increased electron density on the coordinating olefin, as well as increased electron-donating character in the ligand system. The regioselectivity is shown to depend on the electronic and steric characteristics of the olefin, with steric bulk and electron withdrawing character favoring linear product formation. Activation of the C-H bond occurs in a concerted fashion through a novel transition structure best described as an oxidative hydrogen migration. The character of the transition structure is seven coordinate Ir(V), with a full bond formed between the migrating hydrogen and iridium. Several experimental observations are investigated and explained: (a) The nature of L influences the rate of the reaction through a ground-state effect. (b) The lack of beta-hydride products is due to kinetic factors, although beta-hydride elimination is calculated to be facile, all further reactions are kinetically inaccessible. (c) Inhibition by excess olefin is caused by competitive binding of olefin and aryl starting materials during the catalytic cycle in a statistical fashion. On the basis of this insertion-oxidative hydrogen transfer mechanism we suggest that electron-withdrawing substituents on the acac ligands, such as trifluoromethyl groups, are good modifications for catalysts with higher activity.     

Initial nitridation of the Ge(100)-2 x 1 surface by ammonia

The reaction of ammonia (NH(3)) on the Ge(100)-2 x 1 surface is investigated using density functional theory (DFT). We find that NH(3) adsorbs molecularly onto Ge(100)-2 x 1 via the formation of a dative bond. The calculations also show that, unlike Si(100)-2 x 1, the activation barrier for subsequent dissociation of NH(3) adsorbed on Ge(100)-2 x 1 is higher than that of reversible desorption, which indicates that NH(3) has a low reactive sticking probability on the Ge(100)-2 x 1 surface. We also predict that nitrogen insertion into the Ge-Ge dimer requires NH(3) overexposure because the activation barrier for NH(2) insertion into the Ge-Ge dimer is significantly above the entrance channel. The nitridation reaction pathway results in the N-H bridge-bonded state, which is found to be 17.4 kcal/mol more stable than the reactants. We find that the reactions of NH(3) on the Ge(100)-2 x 1 surface generally involve higher activation barriers and less stable intermediates than the analogous reactions on the Si(100)-2 x 1 surface.     

Substituent Effects on Site Selectivity (C=C vs C=N) in Heterocumulene-Heterodiene [4 + 2] Cycloadditions: Density Functional and Semiempirical AM1 Molecular Orbital Calculations

The effect of substituents on the site selectivity (C=C vs C=N) in the [4 + 2] cycloaddition between heterocumulenes (ketene imines) 2a-g with heterodienes (acroleines 9a-n and 4-acylfuran-2,3-diones 1a-d) is treated by semiempirical AM1 molecular orbital and density functional calculations using Becke's three-parameter hybrid method (B3LYP/6-31G). For some reactions calculations were also done at the B3LYP/6-31+G level of theory. For reaction of the oxa 1,3-dienes with ketene imines unsubstituted at the terminal carbon invariably addition across the C=C heterocumulene double bond has a lower activation energy than addition across the C=N double bond. Substitution of methyl or especially phenyl groups at the ketene imine C-terminus leads to a reversal of the respective activation energies. Incorporation of the oxa 1,3-diene system into the heterocyclic dione 1 substantially enhances the reactivity ( approximately 10 kcal mol(-1) lower activation energies) as compared to similarly substituted acroleins. At the DFT level of theory all reactions are found to proceed via a concerted asynchronous mechanism.     

Structure and C-S bond cleavage in aryl 1-methyl-1-arylethyl sulfide radical cations

Steady state and laser flash photolysis (LFP) of a series of p-X-cumyl phenyl sulfides (4-X-C(6)H(4)C(CH(3))(2)SC(6)H(5): 1, X = Br; 2, X = H; 3, X = CH(3); 4, X = OCH(3)) and p-X-cumyl p-methoxyphenyl sulfides (4-X-C(6)H(4)C(CH(3))(2)SC(6)H(4)OCH(3): 5, X = H; 6, X = CH(3); 7, X = OCH(3)) has been carried out in the presence of N-methoxy phenanthridinium hexafluorophosphate (MeOP(+)PF(6)(-)) under nitrogen in MeCN. Steady state photolysis showed the formation of products deriving from the C-S bond cleavage in the radical cations 1(+?)-7(+?) (2-aryl-2-propanols and diaryl disulfides). Formation of 1(+?)-7(+?) was also demonstrated by LFP experiments evidencing the absorption bands of the radical cations 1(+?)-3(+?) (λ(max) = 530 nm) and 5(+?)-7(+?) (λ(max) = 570 nm) mainly localized in the arylsulfenyl group and radical cation 4(+?) (λ(max) = 410, 700 nm) probably mainly localized in the cumyl ring. The radical cations decayed by first-order kinetics with a process attributable to the C-S bond cleavage. On the basis of DFT calculations it has been suggested that the conformations most suitable for C-S bond cleavage in 1(+?)-4(+?) and 7(+?) are characterized by having the C-S bond almost collinear with the π system of the cumyl ring and by a significant charge and spin delocalization from the ArS ring to the cumyl ring. Such a delocalization is probably at the origin of the observation that the rates of C-S bond cleavage result in very little sensitivity to changes in the C-S bond dissociation free energy (BDFE). A quite large reorganization energy value (λ = 43.7 kcal mol(-1)) has been calculated for the C-S bond scission reaction in the radical cation. This value is much larger than that (λ = 12 kcal mol(-1)) found for the C-C bond cleavage in bicumyl radical cations, a reaction that also leads to cumyl carbocations.     

Transition state characterization for the reversible binding of dihydrogen to bis(2,2'-bipyridine)rhodium(I) from temperature- and pressure-dependent experimental and theoretical studies

Thermodynamic and kinetic parameters for the oxidative addition of H2 to [Rh(I)(bpy)2]+ (bpy = 2,2'-bipyridine) to form [Rh(III)(H)2(bpy)2]+ were determined from either the UV-vis spectrum of equilibrium mixtures of [Rh(I)(bpy)2]+ and [Rh(III)(H)2(bpy)2]+ or from the observed rates of dihydride formation following visible-light irradiation of solutions containing [Rh(III)(H)2(bpy)2]+ as a function of H2 concentration, temperature, and pressure in acetone and methanol. The activation enthalpy and entropy in methanol are 10.0 kcal mol(-1) and -18 cal mol(-1) K(-1), respectively. The reaction enthalpy and entropy are -10.3 kcal mol(-1) and -19 cal mol(-1) K(-1), respectively. Similar values were obtained in acetone. Surprisingly, the volumes of activation for dihydride formation (-15 and -16 cm(3) mol(-1) in methanol and acetone, respectively) are very close to the overall reaction volumes (-15 cm(3) mol(-1) in both solvents). Thus, the volumes of activation for the reverse reaction, elimination of dihydrogen from the dihydrido complex, are approximately zero. B3LYP hybrid DFT calculations of the transition-state complex in methanol and similar MP2 calculations in the gas phase suggest that the dihydrogen has a short H-H bond (0.823 and 0.810 Angstroms, respectively) and forms only a weak Rh-H bond (1.866 and 1.915 Angstroms, respectively). Equal partial molar volumes of the dihydrogenrhodium(I) transition state and dihydridorhodium(III) can account for the experimental volume profile found for the overall process.     

A theoretical study on CO2 insertion into an M(bond)H bond (M = Rh and Cu)

Insertion of CO₂ into the transition metal-hydride bond of [Rh^IIIH₂(PH₃)₃]⁺, Cu^IH(PH₃)₂, and Rh^IH(PH₃)₃ was theoretically investigated with ab initio MO/MP 4, SD-CI , and CCD methods. The geometries of reactants, transition states (TS ), and products were optimized at the Hartree-Fock level, and then MP 4, SD-CI , and CCD calculations were performed on those optimized structures. The TS of the CO₂ insertion into the Cu^I(bond)H bond is the most reactantlike, while the TS of the CO₂ insertion into the Rh^III(bond)H bond is the most productlike. The activation energy (E_a) and the reaction energy (ΔE) were calculated to be 6.5 and −33.5 kcal/mol for the CO₂ insertion into the Cu¹(bond)H bond, 21.2 and −7.0 kcal/mol for the CO₂ insertion into the Rh¹(bond)H bond, and 51.3 and −1.1 kcal/mol for the Rh¹¹¹(bond)H bond at the SD-CI level, where negative ΔE represents exothermicity. These results are discussed in terms of the M(bond)H bond energy and the trans-influence of the hydride ligand. © 1996 John Wiley & Sons, Inc.     

Chemo- and periselectivity in the addition of [OsO2(CH2)2] to ethylene: a theoretical study

Quantum chemical calculations by using density functional theory at the B3LYP level have been carried out to elucidate the reaction course for the addition of ethylene to [OsO2(CH2)2] (1). The calculations predict that the kinetically most favorable reaction proceeds with an activation barrier of 8.1 kcal mol(-1) via [3+2] addition across the O=Os=CH2 moiety. This reaction is -42.4 kcal mol(-1) exothermic. Alternatively, the [3+2] addition to the H2C=Os=CH2 fragment of 1 leads to the most stable addition product 4 (-72.7 kcal mol(-1)), yet this process has a higher activation barrier (13.0 kcal mol(-1)). The [3+2] addition to the O=Os=O fragment yielding 2 is kinetically (27.5 kcal mol(-1)) and thermodynamically (-7.0 kcal mol(-1)) the least favorable [3+2] reaction. The formal [2+2] addition to the Os=O and Os=CH2 double bonds proceeds by initial rearrangement of 1 to the metallaoxirane 1 a. The rearrangement 1-->1 a and the following [2+2] additions have significantly higher activation barriers (>30 kcal mol(-1)) than the [3+2] reactions. Another isomer of 1 is the dioxoosmacyclopropane 1 b, which is 56.2 kcal mol(-1) lower in energy than 1. The activation barrier for the 1-->1 b isomerization is 15.7 kcal mol(-1). The calculations predict that there are no energetically favorable addition reactions of ethylene with 1 b. The isomeric form 1 c containing a peroxo group is too high in energy to be relevant for the reaction course. The accuracy of the B3LYP results is corroborated by high level post-HF CCSD(T) calculations for a subset of species.     

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.