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

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

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.     

Thinking out of the black box: accurate barrier heights of 1,3-dipolar cycloadditions of ozone with acetylene and ethylene

Accurate barriers for the 1,3-dipolar cycloadditions of ozone with acetylene and ethylene have been determined via the systematic extrapolation of ab initio energies within the focal point approach of Allen and co-workers. Electron correlation has been accounted for primarily via coupled cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [CCSD, CCSD(T), CCSDT, and CCSDT(Q)]. For the concerted [4 + 2] cycloadditions, the final recommended barriers are DeltaH(0K) = 9.4 +/- 0.2 and 5.3 +/- 0.2 kcal mol(-1) for ozone adding to acetylene and ethylene, respectively. These agree with recent results of Cremer et al. and Anglada et al., respectively. The reaction energy for O3 + C2H2 exhibits a protracted convergence with respect to inclusion of electron correlation, with the CCSDT/cc-pVDZ and CCSDT(Q)/cc-pVDZ values differing by 2.3 kcal mol-1. Recommended enthalpies of formation (298 K) for cycloadducts 1,2,3-trioxole and 1,2,3-trioxolane are +32.8 and -1.6 kcal mol(-1), respectively. Popular composite ab initio approaches [CBS-QB3, CBS-APNO, G3, G3B3, G3(MP2)B3, G4, G4(MP3), and G4(MP2)] predict a range of barrier heights for these systems. The CBS-QB3 computed barrier for ozone and acetylene, DeltaH(0K) = 4.4 kcal mol(-1), deviates by 5 kcal mol(-1) from the focal point value. CBS-QB3 similarly underestimates the barrier for the reaction of ozone and ethylene, yielding a prediction of only 0.7 kcal mol(-1). The errors in the CBS-QB3 results are significantly larger than mean errors observed in application to the G2 test set. The problem is traced to the nontransferability of MP2 basis set effects in the case of these reaction barriers. The recently published G4 and G4(MP2) approaches perform substantially better for O3 + C2H2, predicting enthalpy barriers of 9.0 and 8.4 kcal mol(-1), respectively. For the prediction of these reaction barriers, the additive corrections applied in the majority of the composite approaches considered lead to worse agreement with the reference focal point values than would be obtained relying only on single point energies evaluated at the highest level of theory utilized within each composite method.     

Conformational studies in the cyclohexane series. 2. Phenylcyclohexane and 1-methyl-1-phenylcyclohexane

The structures and relative energies of the conformers of phenylcyclohexane, and 1-methyl-1-phenylcyclohexane have been calculated at theoretical levels including HF/6-31G, B3LYP/6-311G, MP2/6-311G, MP2/6-311(2df,p), QCISD/6-311G, and QCISD/6-311G(2df,p). The latter gives conformational enthalpy (DeltaH degrees ), entropy (DeltaS degrees ), and free energy (DeltaG degrees ) values for phenylcyclohexane that are in excellent agreement with the experimental data. The calculations for 1-methyl-1-phenylcyclohexane find a free energy difference of 1.0 kcal/mol at -100 degrees C, favoring the conformation having an axial phenyl group, that is in only modest agreement with the experimental value of 0.32 +/- 0.04 kcal/mol. The origin of the phenyl rotational profiles for the conformers of phenylcyclohexane and 1-methyl-1-phenylcyclohexane is discussed.     

Theoretical predictions of the initial decomposition steps of dimethylnitramine

The structures and energies of the reactants, products, and transition states of the initial steps in the gas-phase decomposition of dimethylnitramine (DMNA) have been determined by quantum chemical calculations at the B3LYP density-functional theory, MP2, and G2 levels. The pathways considered are NO2 elimination, HONO elimination, and nitro-nitrite rearrangement. The NO2 elimination is predicted to be the main channel of the gas-phase decomposition of DMNA in accord with experiment. The values of the Arrhenius parameters, log A=16.6+/-0.5 and Ea=40.0+/-0.6 kcal/mol, for the N-NO2 bond-fission reaction were obtained using a canonical variational theory with B3LYP energies and frequencies. The HONO-elimination channel has the next lowest activation energy of 44.7+/-0.5 kcal/mol (log A=13.6+/-0.5) and is characterized by a five-member transition-state configuration in which a hydrogen atom from one of the methyl groups is transferred to an oxygen atom of NO2. Tunneling contributions to the rate of this reaction have been estimated. The nitro-nitrite rearrangement reaction occurs via a transition state in which both oxygen atoms of NO2 are loosely bound to the central nitrogen atom, for which Rice-Ramsperger-Kassel-Marcus theory predicts log A=14.4+/-0.6 and Ea=54.1+/-0.8 kcal/mol.     

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.     

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.     

The ozonolysis of acetylene--a quantum chemical investigation.

The ozonolysis of acetylene was investigated using CCSD(T), CASPT2, and B3LYP-DFT in connection with a 6-311+G(2d,2p) basis set. The reaction is initiated by the formation of a van der Waals complex followed by a [4pi + 2pi] cycloaddition between ozone and acetylene (activation enthalpy DeltaH(a)(298) = 9.6 kcal/mol; experiment, 10.2 kcal/mol), yielding 1,2,3-trioxolene, which rapidly opens to alpha-ketocarbonyl oxide 5. Alternatively, an O atom can be transferred from ozone to acetylene (DeltaH(a)(298) = 15.6 kcal/mol), thus leading to formyl carbene, which can rearrange to oxirene or ketene. The key compound in the ozonolysis of acetylene is 5 because it is the starting point for the isomerization to the corresponding dioxirane 19 (DeltaH(a)(298) = 16.9 kcal/mol), for the cyclization to trioxabicyclo[2.1.0]pentane 10 (DeltaH(a)(298) = 19.5 kcal/mol), for the formation of hydroperoxy ketene 15 (DeltaH(a)(298) = 20.6 kcal/mol), and for the rearrangement to dioxetanone 9 (DeltaH(a)(298) = 23.6 kcal/mol). Compounds 19, 10, 15, and 9 rearrange or decompose with barriers between 13 and 16 kcal/mol to yield as major products formanhydride, glyoxal, formaldehyde, formic acid, and (to a minor extent) glyoxylic acid. Hence, the ozonolysis of acetylene possesses a very complicated reaction mechanism that deserves intensive experimental studies.     

Reductive elimination/oxidative addition of carbon-hydrogen bonds at Pt(IV)/Pt(II) centers: mechanistic studies of the solution thermolyses of Tp(Me2)Pt(CH3)2H

Reductive elimination of methane occurs upon solution thermolysis of kappa(3)-Tp(Me)2Pt(IV)(CH(3))(2)H (1, Tp(Me)2 = hydridotris(3,5-dimethylpyrazolyl)borate). The platinum product of this reaction is determined by the solvent. C-D bond activation occurs after methane elimination in benzene-d(6), to yield kappa(3)-Tp(Me)2Pt(IV)(CH(3))(C(6)D(5))D (2-d(6)), which undergoes a second reductive elimination/oxidative addition reaction to yield isotopically labeled methane and kappa(3)-Tp(Me)2Pt(IV)(C(6)D(5))(2)D (3-d(11)). In contrast, kappa(2)-Tp(Me)2Pt(II)(CH(3))(NCCD(3)) (4) was obtained in the presence of acetonitrile-d(3), after elimination of methane from 1. Reductive elimination of methane from these Pt(IV) complexes follows first-order kinetics, and the observed reaction rates are nearly independent of solvent. Virtually identical activation parameters (DeltaH(++)(obs) = 35.0 +/- 1.1 kcal/mol, DeltaS(++)(obs) = 13 +/- 3 eu) were measured for the reductive elimination of methane from 1 in both benzene-d(6) and toluene-d(8). A lower energy process (DeltaH(++)(scr) = 26 +/- 1 kcal/mol, DeltaS(++)(scr) = 1 +/- 4 eu) scrambles hydrogen atoms of 1 between the methyl and hydride positions, as confirmed by monitoring the equilibration of kappa(3)-Tp(Me)()2Pt(IV)(CH(3))(2)D (1-d(1)()) with its scrambled isotopomer, kappa(3)-Tp(Me)2Pt(IV)(CH(3))(CH(2)D)H (1-d(1'). The sigma-methane complex kappa(2)-Tp(Me)2Pt(II)(CH(3))(CH(4)) is proposed as a common intermediate in both the scrambling and reductive elimination processes. Kinetic results are consistent with rate-determining dissociative loss of methane from this intermediate to produce the coordinatively unsaturated intermediate [Tp(Me)2Pt(II)(CH(3))], which reacts rapidly with solvent. The difference in activation enthalpies for the H/D scrambling and C-H reductive elimination provides a lower limit for the binding enthalpy of methane to [Tp(Me)2Pt(II)(CH(3))] of 9 +/- 2 kcal/mol.     

Hydrogen atom transfer from iron(II)-tris[2,2'-bi(tetrahydropyrimidine)] to TEMPO: a negative enthalpy of activation predicted by the Marcus equation

The transfer of a hydrogen atom from iron(II)-tris[2,2'-bi(tetrahydropyrimidine)], [FeII(H2bip)3]2+, to the stable nitroxide, TEMPO, was studied by stopped-flow UV-vis spectrophotometry. The products are the deprotonated iron(III) complex [FeIII(H2bip)2(Hbip)]2+ and the hydroxylamine, TEMPO-H. This reaction can also be referred to as proton-coupled electron transfer (PCET). The equilibrium constant for the reaction is close to 1; thus, the reaction can be driven in either direction. The rate constants for the forward and reverse reactions at 298 K are k1 = 260 +/- 30 M-1 s-1 and k-1 = 150 +/- 20 M-1 s-1. Interestingly, the rate constant for the forward reaction decreases as reaction temperature is increased, implying a negative activation enthalpy: DeltaH1 = -2.7 +/- 0.4 kcal mol-1, DeltaS1 = -57 +/- 8 cal mol-1 K-1. Marcus theory predicts this unusual temperature dependence on the basis of independently measured self-exchange rate constants and equilibrium constants: DeltaHcalcd = -3.5 +/- 0.5 kcal mol-1, DeltaScalcd = -42 +/- 10 cal mol-1 K-1. This result illustrates the value of the Marcus approach for these types of reactions. The dominant contributor to the negative activation enthalpy is the favorable enthalpy of reaction, DeltaH1 degrees = -9.4 +/- 0.6 kcal mol-1, rather than the small negative activation enthalpy for the H-atom self-exchange between the iron complexes.