Molecular mechanics studies (MM4) of sulfides and mercaptans

The MM4 force field has been extended to the title class of compounds. The vibrational spectra, structures, conformational equilibria, and heats of formation have been studied for 47 conformers of 29 compounds. In general, the properties may be calculated with accuracy that is competitive with that for hydrocarbons. The structures are better fit than previously because of the inclusion of a torsion–bend interaction term, which has its origin in the lone pair (Bohlmann) effect. Available experimental data do not suffice to yield detailed torsional potentials, or geometries as a function of torsion angle, and these quantities were determined by ab initio calculations at the MP2/6-31G* level. The rms error in the calculated frequencies of seven representative structures (with a total of 64 experimental and 96 ab initio frequencies) is 25 cm⁻¹. The heats of formation for 23 compounds have a weighted rms error of 0.36 kcal/mol. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1827–1847, 1997     

Molecular mechanics (MM4) study of amines

The MM4 force field has been extended to include aliphatic amines. About 20 amines have been examined to obtain a set of useful molecular mechanics parameters for this class. The vibrational spectra of seven amines (172 frequencies) calculated by MM4 have an overall rms error of 27 cm(-1), compared with corresponding MM4 value of 24 cm(-1) for alkanes. The rms and signed average errors of the moments of inertia of nine simple amines compared with the experimental data were 0.18% and -0.004%, respectively. The heats of formation of 30 amines were also studied. The MM4 weighted standard deviation is 0.41 kcal/mol, compared with experiment. Electronegativity effects occur in the hydrocarbon portion of an amine from the nitrogen, and are accounted for by including electronegativity induced changes in bond lengths and angles, and induced dipole-dipole interactions in the molecule. Negative hyperconjugation results from the presence of the lone pair of electrons on nitrogen, and leads to the Bohlmann bands in the infrared, and also to strong and unusual geometric changes in the molecules (Bohlmann effect), all of which are fairly well accounted for. The conformational energies in amines appear to be less straightforward than those for most other classes of molecules, apparently because of the Bohlmann effect, and these are probably not yet completely understood. In general, the agreement between the MM4 calculated results and the available data is reasonably good.     

MM2' calculations on methylenecyclohexane,methylenecyclopentane, and cyclopentane. Pitfalls in the two-bond drive technique: How large should the ring be?

Molecular mechanics calculations (MM2') are reported on methylenecyclohexane, 1 , methylenecyclopentane, 2 , and cyclopentane itself, 3. The calculated torsional energy barrier for the chair/chair interconversion of 1 if 8.7 kcal/mol (experimental ΔH ^≠ = 8.4 ±0.1 kcal/mol); compounds 2 and 3 have virtually free pseudorotational pathways (calculated ΔH ^≠ = 2.33 and 0.008 kcal/mol, respectively).     

Molecular mechanics (MM4) calculations on carbonyl compounds part I: aldehydes

Aliphatic aldehydes have been studied with the aid of the MM4 force field. The structures, moments of inertia, vibrational spectra, conformational energies, barriers to internal rotation, and dipole moments have been examined for six compounds (nine conformations). MM4 parameters have been developed to fit the indicated quantities to the wide variety of experimental data. Ab initio (MP2) and density functional theory (B3LYP) calculations have been used to augment and/or replace experimental data, as appropriate. Because more, and to some extent, better, data have become available since MM3 was developed, it was anticipated that the overall accuracy of the information calculated with MM4 would be better than with MM3. The best single measure of the overall accuracy of a force field is the accuracy to which the moments of inertia of a set of compounds (from microwave spectroscopy) can be reproduced. For all of the 20 moments (seven conformations) experimentally known for the aldehyde compounds, the MM4 rms error is 0.30%, while with MM3, the most accurate force field presently available, the rms error over the same set is 1.01%. The calculation of the vibrational spectra was also improved overall. For the four aldehydes that were fully analyzed (over a total of 78 frequencies), the rms errors with MM4 and MM3 are 18 and 38 cm^?1, respectively. These improvements came from several sources, but the major ones were separate parameters involving the carbonyl carbon for formaldehyde, the alkyl aldehydes and the ketones, and new crossterms featured in the MM4 force field that are not present in the MM3 version. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1396–1425, 2001     

Alcohols,ethers, carbohydrates,and related compounds. I. The MM4 force field for simple compounds

Simple alcohols and ethers have been studied with the MM4 force field. The structures of 13 molecules have been well fit using the MM4 force field. Moments of inertia have been fit with rms percentage errors as indicated: 18 moments for ethers, 0.28%; 21 moments for alcohols, 0.22%. Rotational barriers and conformational equilibria have also been examined, and the experimental and ab initio results are reproduced substantially better with MM4 than they were with MM3. Much of the improvement comes from the use of additional interaction terms in the force constant matrix, of which the torsion-bend and torsion-torsion are particularly important. Induced dipoles are included in the calculation, and dipole moments are reasonably well fit. It has been possible for the first time to fit conformational energetic data for both open chain and cyclic alcohols (e.g., propanol and cyclohexanol) with the same parameter set. For vibrational spectra, over a total of 82 frequencies, the rms error is 27 cm(-1), as opposed to 38 cm(-1) with MM3. Both the alpha and beta bond shortening resulting from the presence of the electronegative oxygen atom in the molecule are well reproduced. The electronegativity of the oxygen is sufficient that one must also include not only the alpha and beta electronegativity effects on bond lengths, but also on angle distortions, if structures are to be well reproduced. The heats of formation of 32 alcohols and ethers were fit overall to within experimental error (weighted standard deviation error 0.26 kcal/mol).     

Conformational properties of permethylcyclohexane as compared to cyclohexane: A force field study

The conformational properties of the recently synthesized highly strained permethylcyclohexane molecule 2 have been studied by empirical force field calculations using three different potentials (CFF, MM2, MM2′) and second-derivative optimization methods. A comparison of the results with the conformational behavior of parent cyclohexane 1 leads to the following conclusions: The best conformation of 2 is a chair minimum whose six-membered ring is flatter than that of 1 , due to the strong H…H repulsions introduced by the methyl groups. The twist minimum of 2 is energetically less favorable than the chair by an amount similar to 1 . A potential energy barrier Δ V^# for the chair inversion of 2 of 15.32 kcal/mol results with the CFF, only about three kcal/mol higher than for 1 . The free energy of activation ΔG^# for this process obtained with the CFF is 16.96 kcal/mol (at 333 K) and agrees well with the experimental value of 16.7(2) kcal/mol.¹ MM2 and MM2′ give substantially lower and higher potential energy inversion barriers Δ V^# of 9.03 and 20.29 kcal/mol, respectively, which is attributed to inappropriate torsional energy terms in these force fields. The characteristic difference in the conformational behavior of 2 and 1 concerns the boat forms which are substantially less favorable in the per-methyl compound than in 1 . Expectedly, strong H…H repulsions between the 1,4 diaxial flagpole–bowsprit methyl groups in 2 are responsible for this difference. The particularly high strain of the boat forms of 2 leads to flexibility differences as compared to 1 which in turn affect the relative entropies of the various statiomers (stationary point conformations); e.g., the chair ring inversion activation entropies of 2 and 1 are predicted by the CFF calculations to have opposite signs (?4.82 and 3.41 cal/mol K, respectively, at 298 K). The twist and half-twist statiomers of 2 are much more rigid than those of 1 , which is a consequence of the substantially larger boat barriers along their pseudorotational interconversion paths. The boat transition state separating two enantiomeric twist minima represents a barrier calculated to be more than tenfold higher for 2 than for 1 (CFF Δ V^# values 11.14 and 0.92 kcal/mol, respectively); likewise the half-boat chair inversion barrier of 2 is calculated 5.07 kcal/mol less favorable than the respective half-twist barrier. These statiomers are practically equienergetic in the case of 1 . Except for the axial flagpole–bowsprit CH₃ substituents of the boat forms, the methyl groups of all the relevant calculated statiomers of 2 are more or less staggered. The rotational barrier of the equatorial methyl groups of the chair minimum of 2 is computationally predicted to be 5.78 kcal/mol (ΔG^#), i.e., unusually high. Interesting vibrational effects are brought about by the strong H…H repulsions in 2 ; thus the chair minimum has a largest C? H stretching frequency estimated to be 3050 cm^?1 and involves several particularly low frequencies which have a substantial influence on its entropy. CFF calculations for the lower homologue permethylcyclopentane 5 indicate that its pseudorotational properties are similar to those of cyclopentane 4 , in contradistinction to the pair 2/1 .     

Alcohols,ethers, carbohydrates,and related compounds. IV. Carbohydrates

Ab initio calculations [B3LYP/6-311++G(2d,2p)] have been carried out on 84 conformations of 12 different sugars (hexoses), in both pyranose and furanose forms, with the idea of generating a data base for carbohydrate structural energies that may be used for developing the predictive value of molecular mechanics calculations for carbohydrates. The average value for the apparent gas phase anomeric effect for a series of 31 pairs of pyranose conformations was found to be 1.83 kcal/mol (vs. 2.67 kcal/mol with a smaller basis set used in earlier calculations). In developing MM4 to reproduce these data, it was necessary first to have good energies for simple alcohols and ethers, together with an adequate treatment of hydrogen bonding, and then to include the anomeric effect, and the ethylene glycol type system, as was previously recognized. It was also found that the so-called delta-2 effect, long recognized in carbohydrates, must be explicitly included, in order to obtain acceptable results. When a force field that included all of these items as developed from the small molecules based on the MM4 hydrocarbon force field was applied without any parameter adjustment to the set of hexopyranose and furanose conformations mentioned earlier, the E(beta) - E(alpha) was found to have an average value of 1.88 kcal/mol, versus 1.74 for the quantum calculations. The signed average and RMS deviations of the MM4 from the QM results were +0.15 and 0.87 kcal/mol.     

Passive and active oxidation of Si(100) by atomic oxygen: a theoretical study of possible reaction mechanisms

Reaction mechanisms for oxidation of the Si(100) surface by atomic oxygen were studied with high-level quantum mechanical methods in combination with a hybrid QM/MM (Quantum mechanics/Molecular Mechanics) method. Consistent with previous experimental and theoretical results, three structures, "back-bond", "on-dimer", and "dimer-bridge", are found to be the most stable initial surface products for O adsorption (and in the formation of SiO(2) films, i.e., passive oxidation). All of these structures have significant diradical character. In particular, the "dimer-bridge" is a singlet diradical. Although the ground state of the separated reactants, O+Si(100), is a triplet, once the O atom makes a chemical bond with the surface, the singlet potential energy surface is the ground state. With mild activation energy, these three surface products can be interconverted, illustrating the possibility of the thermal redistribution among the initial surface products. Two channels for SiO desorption (leading to etching, i.e., active oxidation) have been found, both of which start from the back-bond structure. These are referred to as the silicon-first (SF) and oxygen-first (OF) mechanisms. Both mechanisms require an 89.8 kcal/mol desorption barrier, in good agreement with the experimental estimates of 80-90 kcal/mol. "Secondary etching" channels occurring after initial etching may account for other lower experimental desorption barriers. The calculated 52.2 kcal/mol desorption barrier for one such secondary etching channel suggests that the great variation in reported experimental barriers for active oxidation may be due to these different active oxidation channels.     

Restricted geometry optimization: a different way to estimate stabilization energies for aromatic molecules of various types

At RHF, MPn, and DFT levels, a procedure of geometry optimization under the restrictions of pi-orbital interactions (GOR) was developed, thus providing a conjugated molecule with the following two types of localized reference geometries: a "GL" geometry where all double bonds are localized, and n different "GE-n" geometries, in each of which only two double bonds were permitted to conjugate. Interestingly, the molecular energy differences between the corresponding pairs of GE-n and GL geometries were found to be additive in each of the acyclic polyenes, and these were not additive for benzene. As a result, an extra stabilization energy (ESE) value of -39.0 kcal/mol was found in benzene. Afterward, GOR was applied to benzene- and furan-like species, strained aromatic molecules, and substituted benzenes, and the calculated ESEs for these molecules were found to be in reasonable ranges. The GOR can isolate a specific group from other groups, and it has several special functions. First, with regard to the substituent effect, the ESE difference between substituted benzene and benzene can be partitioned into conjugative and inductive parts. Second, the behavior of strained aromatic molecules can be ascertained from the roles of their resonance interactions, strained-induced bond localization (SIBL), and inductive effects, indicating that it is resonance interactions, rather than SIBL, which are responsible for localizing double bonds. Emphatically, it is the GL and GE-n geometries of aromatic molecules, rather than nonaromatic compounds, which can be used as the reference structures for calculating ESE. Particularly, these localized geometries are no longer arbitrary.     

Modelling of carbohydrate–aromatic interactions: <Emphasis Type="Italic">ab initio </Emphasis>energetics and force field performance

Summary Aromatic amino acid residues are often present in carbohydrate-binding sites of proteins. These binding sites are characterized by a placement of a carbohydrate moiety in a stacking orientation to an aromatic ring. This arrangement is an example of CH/π interactions. Ab initio interaction energies for 20 carbohydrate–aromatic complexes taken from 6 selected ultra-high resolution X-ray structures of glycosidases and carbohydrate-binding proteins were calculated. All interaction energies of a pyranose moiety with a side chain of an aromatic residue were calculated as attractive with interaction energy ranging from −2.8 to −12.3 kcal/mol as calculated at the MP2/6-311+G(d) level. Strong attractive interactions were observed for a wide range of orientations of carbohydrate and aromatic ring as present in selected X-ray structures. The most attractive interaction was associated with apparent combination of CH/π interactions and classical H-bonds. The failure of Hartree–Fock method (interaction energies from +1.0 to −6.9 kcal/mol) can be explained by a dispersion nature of a majority of the studied complexes. We also present a comparison of interaction energies calculated at the MP2 level with those calculated using molecular mechanics force fields (OPLS, GROMOS, CSFF/CHARMM, CHEAT/CHARMM, Glycam/AMBER, MM2 and MM3). For a majority of force fields there was a strong correlation with MP2 values. RMSD between MP2 and force field values were 1.0 for CSFF/CHARMM, 1.2 for Glycam/AMBER, 1.2 for GROMOS, 1.3 for MM3, 1.4 for MM2, 1.5 for OPLS and to 2.3 for CHEAT/CHARMM (in kcal/mol). These results show that molecular mechanics approximates interaction energies very well and support an application of molecular mechanics methods in the area of glycochemistry and glycobiology.     

Derivation of class II force fields: V. Quantum force field for amides,peptides, and related compounds

As the field of biomolecular structure advances, there is an ever-growing need for accurate modeling of molecular energy surfaces to simulate and predict the properties of these important systems. To address this need, a second generation amide force field for use in simulations of small organics as well as proteins and peptides has been derived. The critical question of what accuracy can be expected from calculations in general, and with this class II force field in particular, is addressed for structural, dynamic, and energetic properties. The force field is derived from a recent methodology we have developed that involves the systematic use of quantum mechanical observables. Systematic ab initio calculations were carried out for numerous configurations of 17 amide and related compounds. Relative energies and first and second derivatives of the energy of 638 structures of these compounds resulted in 140,970 ab initio quantum mechanical observables. The class II peptide quantum mechanical force field (QMFF), containing 732 force constants and reference values, was parameterized against these observables. A major objective of this work is to help establish the role of anharmonicity and coupling in improving the accuracy of molecular force fields, as these terms have not yet become an agreed upon standard in the ever more extensive simulations being used to probe biomolecular properties. This has been addressed by deriving a class I harmonic diagonal force field (HDFF), which was fit to the same energy surface as the QMFF, thus providing an opportunity to quantify the effects of these coupling and anharmonic contributions. Both force field representations are assessed in terms of their ability to fit the observables. They have also been tested by calculating the properties of 11 stationary states of these amide molecules. Optimized structures, vibrational frequencies, and conformational energies obtained from the quantum calculations and from both the QMFF and the HDFF are compared. Several strained and derivatized compounds including urea, formylformamide, and butyrolactam are included in these tests to assess the range of applicability (transferability) of the force fields. It was found that the class II coupled anharmonic force field reproduced the structures, energies, and vibrational frequencies significantly more faithfully than the class I harmonic diagonal force field. An important measure, rms energy deviation, was found to be 1.06 kcal/mol with the class II force field, and 2.30 kcal/mol with the harmonic diagonal force field. These deviations represent the error in relative configurational energy differences for strained and distorted structures calculated with the force fields compared with quantum mechanics. This provides a measure of the accuracy that might be expected in applications where strain may be important such as calculating the energy of a system as it approaches a (rotational) barrier, in ligand binding to a protein, or effects of introducing substituents into a molecule that may induce strain. Similar results were found for structural properties. Protein dynamics is becoming of ever-increasing interest, and, to simulate dynamic properties accurately, the dynamic behavior of model compounds needs to be well accounted for. To this end, the ability of the class I and class II force fields to reproduce the vibrational frequencies obtained from the quantum energy surface was assessed. An rms deviation of 43 cm⁻¹ was achieved with the coupled anharmonic force field, as compared to 105 cm⁻¹ with the harmonic diagonal force field. Thus, the analysis presented here of the class II force field for the amide functional group demonstrates that the incorporation of anharmonicity and coupling terms in the force field significantly improves the accuracy and transferability with regard to the simulation of structural, energetic, and dynamic properties of amides. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 430–458, 1998     

Regioselective addition of Ti enolates to conjugated enones: New insights into the mechanism based on experimental and DFT investigation

The regioselective addition mechanism of the Ti(IV) enolates derived from α-diazo-β-keto carbonyl compounds and α-diazo-β-keto phosphonates to conjugated enones has been studied on the basis of a hypothetical bridging chloride-controlled theory, by density functional theory (DFT), and experimentally. The DFT results indicate that, for the Ti(IV) enolate 3 derived from α-diazo-β-keto carbonyl compounds, the free energy of the bridging chloride-controlled 1,2-addition transition state is 2.4 kcal/mol higher than that of 1,4-addition, and the calculated enthalpies of 1,2-addition is 4.36 kcal/mol more than that of 1,4-addition. For the Ti(IV) enolate 4 derived from α-diazo-β-keto phosphonates, in contrary, the free energy of the bridging chloride-controlled 1,2-addition transition state is 1.1 kcal/mol lower than that of 1,4-addition, and the calculated enthalpy of 1,2-addition is 3.46 kcal/mol less than that of 1,4-addition. Our findings demonstrate that the nucleophilic addition of these Ti(IV) enolates to conjugated enones was carried out not only kinetically but also irreversibly for the first time.     

Heats of formation of organic molecules calculated from ab initio theory and a group equivalent scheme: Alkenes

A bond and group equivalent scheme that allows the calculation of heats of formation of alkenes from ab initio 6-31G* energies has been developed. For a group of 26 compounds, the root mean square (rms) error for the calculated heat of formation was 0.78 kcal/mol. Heats of formation have been predicted for an additional nine compounds for which the experimental values are either unknown or suspect. The heats of hydrogenation of barrelene and related compounds are discussed. © 1994 by John Wiley & Sons, Inc.     

Alcohols,ethers, carbohydrates,and related compounds. III. The 1,2-dimethoxyethane system

Ethylene glycol, its dimethyl ether, and some related compounds have been studied using the MM4 molecular mechanics force field. The MM4 calculated structural and energetic results have been brought into satisfactory agreement with a considerable number of experimental data and MP2/6-311++G(2d,2p) ab initio calculations. The heats of formation of these compounds are also well calculated. The MM4 ethylene glycol conformations in particular are in good agreement, both geometrically and in terms of energy, with those from the ab initio calculations. The corresponding dimethyl ether is of special interest, because it has been suggested that the trans-gauche conformation is unusually stable due to the hydrogen bonding of a hydrogen on a methyl group with the more distant oxygen. It is shown in the present work that while this conformation is more stable than might have been expected, the energy is adequately calculated by MM4 without using any hydrogen bonding between the Cbond;H bond and the oxygen. If such hydrogen bonding occurs, it amounts to no more than about 0.5 kcal/mol in energy, and is too small to detect with certainty. Additionally, energetic relationships in trans-1,2-dimethoxycyclohexane, 1,3,5,7-tetraoxadecalin, and 3-methoxytetrahydropyran have been studied, and the calculated results are compared with experimental information, which is adequately reproduced.     

Improved prediction of heats of formation of energetic materials using quantum mechanical calculations

We present simple atom and group-equivalent methods that will convert quantum mechanical energies of molecules to gas phase heats of formation of CHNO systems. In addition, we predict heats of sublimation and vaporization derived from information obtained from the quantum-mechanically calculated electrostatic potential of each isolated molecule. The heats of sublimation and vaporization are combined with the aforementioned gas phase heats of formation to produce completely predicted condensed phase heats of formation. These semiempirical computational methods, calibrated using experimental information, were applied to a series of CHNO molecules for which no experimental information was used in the development of the methods. These methods improve upon an earlier effort of Rice et al. [Rice, B. M.; Pai, S. V.; Hare, J. Combust. Flame 1999, 118, 445] through the use of a larger basis set and the application of group equivalents. The root-mean-square deviation (rms) from experiment for the predicted group-equivalent gas phase heats of formation is 3.2 kcal/mol with a maximum deviation of 6.5 kcal/mol. The rms and maximum deviation of the predicted liquid heats of formation are 3.2 and 7.4 kcal/mol, respectively. Finally, the rms and maximum deviation of predicted solid heats of formation are 5.6 and 12.2 kcal/mol, respectively, an improvement in the rms of approximately 40% compared to the earlier Rice et al. predictions using atom equivalents and a smaller basis set (B3LYP/6-31G*).     

Heats of formation of organic molecules by ab initio methods: Thiaalkanes

The bond energy scheme for calculating heats of formation of organic molecules from ab initio data (6–31G*) has been extended to include 24 compounds containing sulfur in the sulfide oxidation state. The rms deviation from the experimental values for these compounds is 0.54 kcal/mol, which is approximately experimental error.     

ansa-Zirconocenium catalysis of syndiospecific polymerization of propylene: Theory and experiment

The syndiospecific propylene polymerizations catalyzed by isopropylidene(cyclopentadienyl)(fluorenyl)- and (2,2-dimethylpropylidene)(cyclopentadienyl)(fluorenyl)-zirconocenium ( 1 ⁺ and 2 ⁺) have been investigated theoretically and compared with experimental observations. With the ab initio calculated structures for the transition state (TS) of 1 ⁺(M)P and 2 ⁺(M)P (M = propylene, P = 2-methylpentyl), their steric energies (E°) have been computed using MM2 force-field. The difference between steric energies E°(m) and E°(r) for the meso and racemic enchainment of propylene, respectively, is defined as the stereocontrol energy [δE°(m ? r)] for syndiotactic propagation. The δE°(m ? r) for the TS of 1 ⁺ (M)P is about 2.1 kcal/mol, the value is 1 kcal/mol greater for 2 ⁺(M)P. The observed steric pentad distributions of the syndiotactic poly(propylene) obtained by these catalysts are consistent with smaller effective stereocontrol energy, which is about two-third as large as δE°(m ? r) values calculated for the MM2 optimized structure. Syndiotactic enchainment is favored over isotactic enchainment for all combinations of site configurations in the catalyst. α-Agostic interaction seems to enhance syndioselectivity, whereas γ-agostic interaction changes the stereoselectivity to meso enchainment. The mirror plane symmetry of the syndiotactic propagating species renders the stereoselectivity of the polymerization insensitive to reaction conditions. These catalysts are also highly regiospecific. © 1995 John Wiley & Sons, Inc.     

A multilayered representation,quantum mechanical and molecular mechanics study of the CH3F + OH− reaction in water

The bimolecular nucleophilic substitution (S_N2) reaction of CH₃F + OH^? in aqueous solution was investigated using a combined quantum mechanical and molecular mechanics approach. Reactant complex, transition state, and product complex along the reaction pathway were analyzed in water. The potentials of mean force were calculated using a multilayered representation with the DFT and CCSD(T) level of theory for the reactive region. The obtained free energy activation barrier for this reaction at the CCSD(T)/MM representation is 18.3 kcal/mol which agrees well with the experimental value at ～21.6 kcal/mol. Both the solvation effect and solute polarization effect play key roles on raising the activation barrier height in aqueous solution, with the former raising the barrier height by 3.1 kcal/mol, the latter 1.5 kcal/mol. © 2013 Wiley Periodicals, Inc.     

Heats of formation of organic molecules by Ab Initio calculations: Carboxylic acids and esters

A bond and group equivalent scheme that allows the calculation of heats of formation for carboxylic acids and esters from ab initio 6-31G* energies has been developed. For a group of 16 compounds, the rms error for the calculated heats of formation was 0.64 kcal/mol. Heats of formation have been predicted for an additional seven compounds for which the experimental values are either unknown or suspect. © 1992 by John Wiley & Sons, Inc.     

Accurate heats of formation of the "Arduengo-type" carbene and various adducts including H2 from ab initio molecular orbital theory

The heats of formation of saturated and unsaturated diaminocarbenes (imadazol(in)-2-ylidenes) have been calculated by using high levels of ab initio electronic structure theory. The calculations were done at the coupled cluster level through noniterative triple excitations with augmented correlation consistent basis sets up through quadruple. In addition, four other corrections were applied to the frozen core atomization energies: (1) a zero point vibrational correction; (2) a core/valence correlation correction; (3) a scalar relativistic correction; (4) a first-order atomic spin-orbit correction. The value of DeltaHf( 298) for the unsaturated carbene 1 is calculated to be 56.4 kcal/mol. The value of DeltaHf( 298) for the unsaturated triplet carbene (3)1 is calculated to be 142.8 kcal/mol, giving a singlet-triplet splitting of 86.4 kcal/mol. Addition of a proton to 1 forms 3 with DeltaHf( 298)(3) = 171.6 kcal/mol with a proton affinity for 1 of 250.5 kcal/mol at 298 K. Addition of a hydrogen atom to 1 forms 4 with DeltaHf( 298)(4) = 72.7 kcal/mol and a C-H bond energy of 35.8 kcal/mol at 298 K. Addition of H- to 1 gives 5 with DeltaHf( 298)(5) = 81.2 kcal/mol and 5 is not stable with respect to loss of an electron to form 4. Addition of H2 to the carbene center forms 6 with DeltaHf( 298)(6) = 41.5 kcal/mol and a heat of hydrogenation at 298 K of -14.9 kcal/mol. The value of DeltaHf( 298) for the saturated carbene 7 (obtained by adding H2 to the C=C bond of 1) is 47.4 kcal/mol. Hydrogenation of 7 to form the fully saturated imidazolidine, 8, gives DeltaHf( 298)(8) = 14.8 kcal/mol and a heat of hydrogenation at 298 K of -32.6 kcal/mol. The estimated error bars for the calculated heats of formation are +/-1.0 kcal/mol.