Heats of formation of triplet ethylene, ethylidene, and acetylene

Heats of formation of the lowest triplet state of ethylene and the ground triplet state of ethylidene have been predicted by high level electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T) energies extrapolated to the complete basis set limit using correlation consistent basis sets (CBS), plus additional corrections predict the following heats of formation in kcal/mol: DeltaH0r(C2H4,3A1) = 80.1 at 0 K and 78.5 at 298 K, and DeltaH0t(CH3CH,3A' ') = 86.8 at 0 K and 85.1 at 298 K, with an error of less than +/-1.0 kcal/mol. The vertical and adiabatic singlet-triplet separation energies of ethylene were calculated as DeltaES-T,vert = 104.1 and DeltaES-T,adia = 65.8 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo (DMC) values of 103.5 +/- 0.3 and 66.4 +/- 0.3 kcal/mol. Both sets of computational values differ from the experimental estimate of 58 +/- 3 kcal/mol for the adiabatic splitting. The computed singlet-triplet gap at 0 K for acetylene is DeltaES-T,adia(C2H2) = 90.5 kcal/mol, which is in notable disagreement with the experimental value of 82.6 kcal/mol. The heat of formation of the triplet is DeltaH0tC2H2,3B2) = 145.3 kcal/mol. There is a systematic underestimation of the singlet-triplet gaps in recent photodecomposition experiments by approximately 7 to 8 kcal/mol. For vinylidene, we predict DeltaH0t(H2CC,1A1) = 98.8 kcal/mol at 298 K (exptl. 100.3 +/- 4.0), DeltaH0t(H2CC,3B2) = 146.2 at 298 K, and an energy gap DeltaES-T-adia(H2CC) = 47.7 kcal/mol.     

Heats of formation and singlet-triplet separations of hydroxymethylene and 1-hydroxyethylidene

Thermochemical parameters of hydroxymethylene (HC:OH) and 1-hydroxyethylidene (CH3C:OH) were evaluated by using coupled-cluster, CCSD(T), theory, in conjunction with the augmented correlation consistent, aug-cc-pVnZ, basis sets, with n = D, T, Q, and 5, extrapolated to the complete basis set limit. The predicted value at 298 K for Delta Hf(CH2O) is -26.0 +/- 1 kcal/mol, as compared to an experimental value of -25.98 +/- 0.01 kcal/mol, and for Delta Hf(CH:OH) it is 26.1 +/- 1 kcal/mol. The hydroxymethylene-formaldehyde energy gap is 52.1 +/- 0.5 kcal/mol, the singlet-triplet separation of hydroxymethylene is Delta E(ST)(HC:OH) = 25.3 +/- 0.5 kcal/mol, the proton affinity is PA(HC:OH) = 222.5 +/- 0.5 kcal/mol, and the ionization energy is IEa(HC:OH) = 8.91 +/- 0.04 eV. The predicted value at 298 K for Delta Hf(CH3CHO) is -39.1 +/- 1 kcal/mol as compared to an experimental value of -40.80 +/- 0.35 kcal/mol, and for Delta Hf(CH3C:OH) it is 11.2 +/- 1 kcal/mol. The hydroxyethylidene-acetaldehyde energy gap is 50.6 +/- 0.5 kcal/mol, the singlet-triplet separation of 1-hydroxyethylidene is Delta E(ST)(CH3C:OH) = 30.5 +/- 0.5 kcal/mol, the proton affinity is PA(CH3C:OH) = 234.7 +/- 0.5 kcal/mol, and the ionization energy is IEa(CH3C:OH) = 8.18 +/- 0.04 eV. The calculated energy differences between the carbene and aldehyde isomers, and, thus, the heats of formation of the carbenes, differ from the experimental values by 2.5 kcal/mol.     

Ab initio prediction of the gas- and solution-phase acidities of strong Brønsted acids: the calculation of pKa values less than -10

The intrinsic gas-phase acidities of a series of 21 Br?nsted acids have been predicted with G3(MP2) theory. The G3(MP2) results agree with high level CCSD(T)/CBS acidities for H(2)SO(4), FSO(3)H, CH(3)SO(3)H, and CF(3)SO(3)H to within 1 kcal/mol. The G3(MP2) results are in excellent agreement with experimental gas-phase acidities in the range 342-302 kcal/mol to within <1 kcal/mol for 14 out of 15 acids. Five of the six acids in the range of 302-289 kcal/mol had an average deviation of 5.5 kcal/mol and the strongest acid, (CF(3)SO(2))(3)CH, deviated by 15.0 kcal/mol. These high-level calculations strongly suggest that the experimental acidities in this very acidic part of the scale need to be remeasured. The CCSD(T)/CBS (mixed exponential Gaussian) additive approach for CH(3)CO(2)H, HNO(3), H(2)SO(4), CH(3)SO(3)H, FSO(3)H, and CF(3)SO(3)H gives excellent agreement (+/-1 kcal/mol) with experiment for the DeltaH(f)(0)'s of non-sulfur containing species, and supports the low end of the experimental values for H(2)SO(4) and FSO(3)H. Use of a larger basis set (aug-cc-pV5Z) in the CBS extrapolation improves the agreement with experiment for both H(2)SO(4) and FSO(3)H. The G3(MP2) heats of formation for RSO(3)H molecules tend to be underestimated as compared to the CCSD(T)/CBS approach by 2.5-7.0 kcal/mol. COSMO solvation calculations were used to predict solution free energies and pK(a) values with pK(a)'s up to -17.4. Including the solvation of the proton gives good agreement with experimental pK(a) values in the very acidic regime, whereas it is less reliable for weaker acids. The use of CH(3)CO(2)H and HNO(3) as reference acids in the less acidic and more acidic regions of the scale, respectively, provided improved results to within +/-2 pK(a) units in nearly all cases (+/-3 kcal/mol accuracy).     

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*).     

Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.     

Computational study on the mechanisms and energetics of trimethylindium reactions with H2O and H2S

The reactions of trimethylindium (TMIn) with H2O and H2S are relevant to the chemical vapor deposition of indium oxide and indium sulfide thin films. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/[6-31G(d,p)+Lanl2dz]//B3LYP/[6-31G(d,p)+Lanl2dz] and CCSD(T)/[6-31G(d,p)+Lanl2dz] //MP2/[6-31G(d,p)+Lanl2dz] methods. The results of both methods are in good agreement for the optimized geometries and relative energies. When TMIn reacts with H2O and H2S, initial molecular complexes [(CH3)3In:OH2 (R1)] and [(CH3)3In:SH2 (R2)] are formed with 12.6 and 3.9 kcal/mol binding energies. Elimination of a CH4 molecule from each complex occurs with a similar energy barrier at TS1 (19.9 kcal/mol) and at TS3 (22.1 kcal/mol), respectively, giving stable intermediates (CH3)2InOH and (CH3)2InSH. The elimination of the second CH4 molecule from these intermediate products, however, has to overcome very high and much different barriers of 66.1 and 53.2 kcal/mol, respectively. In the case of DMIn with H2O and H2S reactions, formation of both InO and InS is exothermic by 3.1 and 30.8 kcal/mol respectively. On the basis of the predicted heats of formation of R1 and R2 at 0 K and -20.1 and 43.6 kcal/mol, the heats of formation of (CH3)2InOH, (CH3)2InSH, CH3InO, CH3InS, InO, and InS are estimated to be -20.6, 31.8, and 29.0 and 48.4, 35.5, and 58.5 kcal/mol, respectively. The values for InO and InS are in good agreement with available experimental data. A similar study on the reactions of (CH3)2In with H2O and H2S has been carried out; in these reactions CH3InOH and CH3InSH were found to be the key intermediate products.     

Methyl cation affinities of rare gases and nitrogen and the heat of formation of diazomethane

The methyl cation affinities of the rare gases have been calculated at 0 and 298 K by using coupled cluster theory including noniterative, quasiperturbative triple excitations with the new correlation-consistent basis sets for Xe up through aug-cc-pV5Z in some cases. To achieve near chemical accuracy (+/-1 kcal/mol) in the thermodynamic properties, we add to the estimated complete basis set valence binding energies, based on frozen core coupled cluster theory energies, two corrections: (1) a core/valence correction and (2) a scalar relativistic correction. Vibrational zero-point energies were computed at the coupled cluster level of theory at the CCSD(T)/aug-cc-pVDZ level. The calculated rare gas methyl cation affinities (MCA in kcal/mol) at 298 K are the following: MCA(He) = 1.7, MCA(Ne) = 2.5, MCA(Ar) = 16.9, MCA(Kr) = 25.5, and MCA(Xe) = 36.6. Because of the importance of the MCA(N(2)) in the experimental measurements of the MCA scale, we calculated a number of quantities associated with CH(3)N(2)(+) and CH(2)N(2). The calculated values for diazomethane at 298 K are: DeltaH(f)(CH(2)N(2)) = 65.3 kcal/mol, PA(CH(2)N(2)) = 211.9 kcal/mol, and MCA(N(2)) = 43.2 kcal/mol.     

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).     

Unimolecular dissociation and thermochemistry of CH3OOH

The unimolecular dissociation of CH3OOH is investigated by exciting the molecule in the region of its 5nu(OH) band and probing the resulting OH fragments using laser-induced fluorescence. The measured OH fragment rotational and translational energies are used to determine the CH3O-OH bond dissociation energy, which we estimate to be approximately 42.6+/-1 kcal/mol. Combining this value with the known heats of formation of the fragments also gives an estimate for the heat of formation of CH3OOH which at 0 K we determine to be deltaH(f)0=-27+/-1 kcal/mol. This experimental value is in good agreement with the results of ab initio calculations carried out at the CCSD(T)/complete basis set limit which finds the heat of formation of CH3OOH at 0 K to be deltaH(f)0=-27.3 kcal/mol.     

Accurate thermochemical properties for energetic materials applications. II. Heats of formation of imidazolium-, 1,2,4-triazolium-, and tetrazolium-based energetic salts from isodesmic and lattice energy calculations

A computational approach to the prediction of the heats of formation (DeltaH(f)degrees' s of solid-state energetic salts from electronic structure and volume-based thermodynamics (VBT) calculations is described. The method uses as its starting point reliable DeltaH(f)degrees' s for energetic precursor molecules and ions. The DeltaH(f)degrees' s of more complex energetics species such as substituted imidazole, 1,2,4-triazole, and tetrazole molecules and ions containing amino, azido, and nitro (including methyl) substituents are calculated using an isodesmic approach at the MP2/complete basis set level. On the basis of comparisons to experimental data for neutral analogues, this isodesmic approach is accurate to <3 kcal/mol for the predicted cation and anion DeltaH(f)degrees' s. The DeltaH(f)degrees' s of the energetic salts in the solid state are derived from lattice energy (U(L)) calculations using a VBT approach. Improved values for the alpha and beta parameters of 19.9 (kcal nm)/mol and 37.6 kcal/mol for the U(L) equation were obtained on the basis of comparisons to experimental U(L)' s for a series of 23 salts containing ammonium, alkylammonium, and hydrazinium cations. The total volumes are adjusted to account for differences between predicted and experimental total volumes due to different shapes of the ions (flat vs spherical). The predicted DeltaH(f)degrees' s of the energetic salts are estimated to have error bars of 6-7 kcal/mol, on the basis of comparisons to established experimental DeltaH(f)degrees' s of a subset of the salts studied. Energetic salts with the highest positive DeltaH(f)degrees' s are predicted for azido-containing cations, coupled with heterocyclic anions containing nitro substituents. The substitution of functional groups on carbon versus nitrogen atoms of the heterocyclic cations has interesting stabilization and destabilization effects, respectively.     

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.     

Thermodynamic properties of arsenic compounds and the heat of formation of the As atom from high level electronic structure calculations

Structures, vibrational frequencies, atomization energies at 0 K, and heats of formation at 0 and 298 K are predicted for the compounds As(2), AsH, AsH(2), AsH(3), AsF, AsF(2), and AsF(3) from frozen core coupled cluster theory calculations performed with large correlation consistent basis sets, up through augmented sextuple zeta quality. The coupled cluster calculations involved up through quadruple excitations. For As(2) and the hydrides, it was also possible to examine the impact of full configuration interaction on some of the properties. In addition, adjustments were incorporated to account for extrapolation to the frozen core complete basis set limit, core/valence correlation, scalar relativistic effects, the diagonal Born-Oppenheimer correction, and atomic spin orbit corrections. Based on our best theoretical D(0)(As(2)) and the experimental heat of formation of As(2), we propose a revised 0 K arsenic atomic heat of formation of 68.86 ± 0.8 kcal/mol. While generally good agreement was found between theory and experiment, the heat of formation of AsF(3) was an exception. Our best estimate is more than 7 kcal/mol more negative than the single available experimental value, which argues for a re-examination of that measurement.     

Reliable predictions of the thermochemistry of boron-nitrogen hydrogen storage compounds: BxNxHy, x = 2, 3

Thermochemical data calculated using ab initio molecular orbital theory are reported for 16 BxNxHy compounds with x = 2, 3 and y > or = 2x. Accurate gas-phase heats of formation were obtained using coupled cluster with single and double excitations and perturbative triples (CCSD(T)) valence electron calculations extrapolated to the complete basis set (CBS) limit with additional corrections including core/valence, scalar relativistic, and spin-orbit corrections to predict the atomization energies and scaled harmonic frequencies to correct for zero point and thermal energies and estimate entropies. Computationally cheaper calculations were also performed using the G3MP2 and G3B3 variants of the Gaussian 03 method, as well as density functional theory (DFT) using the B3LYP functional. The G3MP2 heats of formation are too positive by up to approximately 6 kcal/mol as compared with CCSD(T)/CBS values. The more expensive G3B3 method predicts heats of formation that are too negative as compared with the CCSD(T)/CBS values by up to 3-4 kcal/mol. DFT using the B3LYP functional and 6-311+G** basis set predict isodesmic reaction energies to within a few kcal/mol compared with the CCSD(T)/CBS method so isodesmic reactions involving BN compounds and the analogous hydrocarbons can be used to estimate heats of formation. Heats of formation of c-B3N3H12 and c-B3N3H6 are -95.5 and -115.5 kcal/mol at 298 K, respectively, using our best calculated CCSD(T)/CBS approach. The experimental value for c-B3N3H6 appears to be approximately 7 kcal/mol too negative. Enthalpies, entropies, and free energies are calculated for many dehydrocoupling and dehydrogenation reactions that convert BNH6 to alicyclic and cyclic oligomers and H2(g). Generally, the reactions are highly exothermic and exergonic as well because of the release of 1 or more equivalents of H2(g). For c-B3N3H12 and c-B3N3H6, available experimental data for sublimation and vaporization lead to estimates of their condensed phase 298 K heats of formation: DeltaHf degrees [c-B3N3H12(s)] = -124 kcal/mol and DeltaHf degrees [c-B3N3H6(l)] = -123 kcal/mol. The reaction thermochemistries for the dehydrocoupling of BNH6(s) to c-B3N3H12(s) and the dehydrogenation of c-B3N3H12(s) to c-B3N3H6(l) are much less exothermic compared with the gas-phase reactions due to intermolecular forces which decrease in the order BNH6 > cyclo-B3N3H12 > cyclo-B3N3H6. The condensed phase reaction free energies are less negative compared with the gas-phase reactions but are still too favorable for BNH6 to be regenerated from either c-B3N3H12 or c-B3N3H6 by just an overpressure of H2.     

Dissociative photoionization and thermochemistry of dihalomethane compounds studied by threshold photoelectron photoion coincidence spectroscopy

The dissociative photoionization studies have been performed for a set of dihalomethane CH(2)XY (X,Y = Cl, Br, and I) molecules employing the threshold photoelectron photoion coincidence (TPEPICO) technique. Accurate dissociation onsets for the first and second dissociation limits have been recorded in the 10-13 eV energy range, and ionization potentials have been measured for these compounds. By using our experimental dissociation onsets and the known heat of formation of CH(2)Cl(2) molecule, it has been possible to derive the 0 and 298 K heats of formation of all six neutral dihalomethanes as well as their ionic fragments, CH(2)Cl(+), CH(2)Br(+), and CH(2)I(+), to a precision better than 3 kJ/mol. These new measurements serve to fill the lack of reliable experimental thermochemical information on these molecules, correct the old literature values by up to 19 kJ/mol, and reduce their uncertainties. From our thermochemical results it has also been possible to derive a consistent set of bond dissociation energies for the dihalomethanes.     

Optimization of parameters for semiempirical methods II. Applications

MNDO/AM1-type parameters for twelve elements have been optimized using a newly developed method for optimizing parameters for semiempirical methods. With the new method, MNDO-PM3, the average difference between the predicted heats of formation and experimental values for 657 compounds is 7.8 kcal/mol, and for 106 hypervalent compounds, 13.6 kcal/mol. For MNDO the equivalent differences are 13.9 and 75.8 kcal/mol, while those for AM1, in which MNDO parameters are used for aluminum, phosphorus, and sulfur, are 12.7 and 83.1 kcal/mol, respectively. Average errors for ionization potentials, bond angles, and dipole moments are intermediate between those for MNDO and AM1, while errors in bond lengths are slightly reduced.     

Extension of the PDDG/PM3 and PDDG/MNDO semiempirical molecular orbital methods to the halogens

The new semiempirical methods, PDDG/PM3 and PDDG/MNDO, have been parameterized for halogens. For comparison, the original MNDO and PM3 were also reoptimized for the halogens using the same training set; these modified methods are referred to as MNDO' and PM3'. For 442 halogen-containing molecules, the smallest mean absolute error (MAE) in heats of formation is obtained with PDDG/PM3 (5.6 kcal/mol), followed by PM3' (6.1 kcal/mol), PDDG/MNDO (6.6 kcal/mol), PM3 (8.1 kcal/mol), MNDO' (8.5 kcal/mol), AM1 (11.1 kcal/mol), and MNDO (14.0 kcal/mol). For normal-valent halogen-containing molecules, the PDDG methods also provide improved heats of formation over MNDO/d. Hypervalent compounds were not included in the training set and improvements over the standard NDDO methods with sp basis sets were not obtained. For small haloalkanes, the PDDG methods yield more accurate heats of formation than are obtained from density functional theory (DFT) with the B3LYP and B3PW91 functionals using large basis sets. PDDG/PM3 and PM3' also give improved binding energies over the standard NDDO methods for complexes involving halide anions, and they are competitive with B3LYP/6-311++G(d,p) results including thermal corrections. Among the semiempirical methods studied, PDDG/PM3 also generates the best agreement with high-level ab initio G2 and CCSD(T) intrinsic activation energies for S(N)2 reactions involving methyl halides and halide anions. Finally, the MAEs in ionization potentials, dipole moments, and molecular geometries show that the parameter sets for the PDDG and reoptimized NDDO methods reduce the MAEs in heats of formation without compromising the other important QM observables.     

F+ and F⁻ affinities of simple N(x)F(y) and O(x)F(y) compounds

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for the neutral and ionic N(x)F(y) and O(x)F(y) systems using coupled cluster theory with single and double excitations and including a perturbative triples correction (CCSD(T)) method with correlation consistent basis sets extrapolated to the complete basis set (CBS) limit. To achieve near chemical accuracy (±1 kcal/mol), three corrections to the electronic energy were added to the frozen core CCSD(T)/CBS binding energies: corrections for core-valence, scalar relativistic, and first order atomic spin-orbit effects. Vibrational zero point energies were computed at the CCSD(T) level of theory where possible. The calculated heats of formation are in good agreement with the available experimental values, except for FOOF because of the neglect of higher order correlation corrections. The F(+) affinity in the N(x)F(y) series increases from N(2) to N(2)F(4) by 63 kcal/mol, while that in the O(2)F(y) series decreases by 18 kcal/mol from O(2) to O(2)F(2). Neither N(2) nor N(2)F(4) is predicted to bind F(-), and N(2)F(2) is a very weak Lewis acid with an F(-) affinity of about 10 kcal/mol for either the cis or trans isomer. The low F(-) affinities of the nitrogen fluorides explain why, in spite of the fact that many stable nitrogen fluoride cations are known, no nitrogen fluoride anions have been isolated so far. For example, the F(-) affinity of NF is predicted to be only 12.5 kcal/mol which explains the numerous experimental failures to prepare NF(2)(-) salts from the well-known strong acid HNF(2). The F(-) affinity of O(2) is predicted to have a small positive value and increases for O(2)F(2) by 23 kcal/mol, indicating that the O(2)F(3)(-) anion might be marginally stable at subambient temperatures. The calculated adiabatic ionization potentials and electron affinities are in good agreement with experiment considering that many of the experimental values are for vertical processes.     

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).     

Experimental and quantum-chemical studies on photoionization and dissociative photoionization of CH2Br2

The dissociative photoionization of CH2Br2 in a region approximately 10-24 eV was investigated with photoionization mass spectroscopy using a synchrotron radiation source. An adiabatic ionization energy of 10.25 eV determined for CH2Br2 agrees satisfactorily with predictions of 10.26 and 10.25 eV with G2 and G3 methods, respectively. Observed major fragment ions CH2Br+, CHBr+, and CBr+ show appearance energies at 11.22, 12.59, and 15.42 eV, respectively; minor fragment ions CHBr2+, Br+, and CH2+ appear at 12.64, 15.31, and 16.80 eV, respectively. Energies for formation of observed fragment ions and their neutral counterparts upon ionization of CH2Br2 are computed with G2 and G3 methods. Dissociative photoionization channels associated with six observed fragment ions are proposed based on comparison of determined appearance energies and predicted energies. An upper limit of DeltaH0f,298(CHBr+) < or = 300.7 +/- 1.5 kcal mol(-1) is derived experimentally; the adiabatic ionization energy of CHBr is thus derived to be < or = 9.17 +/- 0.23 eV. Literature values for DeltaH0f,298(CBr+) = 362.5 kcal mol(-1) and ionization energy of 10.43 eV for CBr are revised to be less than 332 kcal mol(-1) and 9.11 eV, respectively. Also based on a new experimental ionization energy, DeltaH0f,298(CH2Br2+) is revised to be 236.4 +/- 1.5 kcal mol(-1).     

Preparation and reactivity of [D3d]-octahedrane: the most stable (CH)12 hydrocarbon

The synthesis of the (CH)12 hydrocarbon [D(3d)]-octahedrane (heptacyclo[6.4.0.0(2,4).0(3,7).0(5,12).0(6,10).0(9,11)]dodecane) 1 and its selective functionalization retaining the hydrocarbon cage is described. The B3LYP/6-311+G* strain energy of 1 is 83.7 kcal mol(-1) (4.7 kcal mol(-1) per C-C bond) which is significantly higher than that of the structurally related (CH)16 [D(4d)]-decahedrane 2 (75.4 kcal mol(-1); 3.1 kcal mol(-1) per C-C bond) and (CH)20 [I(h)]-dodecahedrane 3 (51.5 kcal mol(-1); 1.7 kcal mol(-1) per C-C bond); the heats of formation for 1-3 computed according to homodesmotic equations are 52, 35, and 4 kcal mol(-1). Catalytic hydrogenation of 1 leads to consecutive opening of the two cyclopropane rings to give C2-bisseco-octahedrane (pentacyclo[6.4.0.0(2,6).0(3,11).0(4,9)]dodecane) 16 as the major product. Although 1 is highly strained, its carbon skeleton is kinetically quite stable: Upon heating, 1 does not decompose until above 180 degrees C. The B3LYP/6-31G* barriers for the S(R)2 attack of the tBuO. and Br3C. radicals on a carbon atom of one of the cyclopropane fragments (Delta(298) = 27-28 kcal mol(-1)) are higher than those for hydrogen atom abstraction. The latter barriers are virtually identical for the abstraction from the C1-H and C2-H positions with the tBuO. radical (DeltaG(298) = 17.4 and 17.9 kcal mol(-1), respectively), but significantly different for the reaction at these positions with the Br3C. radical (DeltaG(298) = 18.8 and 21.0 kcal mol(-1)). These computational results agree well with experiments, in which the chlorination of 1 with tert-butyl hypochlorite gave a mixture of 1- and 2-chlorooctahedranes (ratio 3:2). The bromination with carbon tetrabromide under phase-transfer catalytic (PTC) conditions (nBu4NBr/NaOH) selectively gave 1-bromooctahedrane in 43 % isolated yield. For comparison, the PTC bromination was also applied to 2,4-dehydroadamantane yielding 54 % 7-bromo-2,4-dehydroadamantane.