Heats of formation of xenon fluorides and the fluxionality of XeF(6) from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(+), XeF(-), XeF(2), XeF(4), XeF(5)(-), and XeF(6) from coupled cluster theory (CCSD(T)) calculations with new correlation-consistent basis sets for Xe. To achieve near chemical accuracy (+/-1 kcal/mol), up to four corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, a correction for first-order atomic spin-orbit effects, and in some cases, a second-order spin-orbit correction. Vibrational zero-point energies were computed at the coupled cluster level of theory. The structure of XeF(6) is difficult to obtain with the C(3)(v)() and O(h)() structures having essentially the same energy. The O(h)() structure is only 0.19 kcal/mol below the C(3)(v)() one at the CCSD(T)/CBS level using an approximate geometry for the C(3)(v)() structure. With an optimized C(3)(v)() geometry, the C(3)(v)() structure would probably become slightly lower in energy than the O(h)() one. The calculated heats of formation for the neutral XeF(n)() fluorides are less negative than the experimental values from the equilibrium measurements by 2.0, 7.7, and 12.2 kcal/mol for n = 2, 4, and 6, respectively. For the experimental values, derived from the photoionization measurements, this discrepancy becomes even larger, suggesting a need for a redetermination of the experimental values. Evidence is presented for the fluxionality of XeF(6) caused by the presence of a sterically active, free valence electron pair on Xe.     

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.     

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.     

Accurate thermochemical properties for energetic materials applications. I. Heats of formation of nitrogen-containing heterocycles and energetic precursor molecules from electronic structure theory

The heats of formation of 1H-imidazole, 1H-1,2,4-trizazole, 1H-tetrazole, CH3NO2, CH3N3, CH3NH2, CH2CHNO2, HClO4, and phenol, as well as cations and anions derived from some of the molecules have been calculated using ab initio molecular orbital theory. These molecules are important as models for compounds used for energetic materials synthesis. The predicted heats of formation of the heterocycle-based compounds are in excellent agreement with available experimental values and those derived from proton affinities and deprotonation enthalpies to <1 kcal/mol. The predicted value for the tetrazolium cation differs substantially from the experimental value, likely due to uncertainty in the measurement. The heats of formation of the nitro and amino molecules, as well as phenol/phenolate, also are in good agreement with the experimental values (<1.5 kcal/mol). The heat of formation of CH3N3 is predicted to be 72.8 kcal/mol at 298 K with an estimated error bar of +/-1 kcal/mol on the basis of the agreement between the calculated and experimental values for DeltaH(f)(HN3). The heat of formation at 298 K of HClO4 is -0.4 kcal/mol, in very good agreement with the experimental value, as well as a W2 literature study. An extrapolation of the CCSD(T)/aug-cc-pV(Q,5) energies was required to obtain this agreement. This result suggests that very large basis sets (> or =aug-cc-pV5Z) may be needed to fully recover the valence correlation energy contribution in compounds containing elements with high formal oxidation states at the central atom. In addition tight d functions are needed for the geometry predictions. Douglas-Kroll-Hess (DKH) scalar relativistic corrections for HClO4 and ClO4- at the MP2 level with correlation-consistent DKH basis sets were predicted to be large, likely due to the high formal oxidation state at the Cl.     

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.     

The heats of formation of diazene, hydrazine, N2H3+, N2H5+, N2H, and N2H3 and the Methyl Derivatives CH3NNH, CH3NNCH3, and CH3HNNHCH3

The heats of formation of N(2)H, diazene (cis- and trans-N(2)H(2)), N(2)H(3), and hydrazine (N(2)H(4)), as well as their protonated species (diazenium, N(2)H(3)(+), and hydrazinium, N(2)H(5)(+)), have been calculated by using high level electronic structure theory. Energies were calculated by using coupled cluster theory with a perturbative treatment of the triple excitations (CCSD(T)) and employing augmented correlation consistent basis sets (aug-cc-pVnZ) up to quintuple-zeta, to perform a complete basis set extrapolation for the energy. Geometries were optimized at the CCSD(T) level with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. Core-valence and scalar relativistic corrections were included, as well as scaled zero point energies. We find the following heats of formation (kcal/mol) at 0 (298) K: DeltaH(f)(N(2)H) = 60.8 (60.1); DeltaH(f)(cis-N(2)H(2)) = 54.9 (53.2); DeltaH(f)(trans-N(2)H(2)) = 49.9 (48.1) versus >/=48.8 +/- 0.5 (exptl, 0 K); DeltaH(f)(N(2)H(4)) = 26.6 (23.1) versus 22.8 +/- 0.2 (exptl, 298 K); DeltaH(f)(N(2)H(3)) = 56.2 (53.6); DeltaH(f)(N(2)H(3)(+)) = 231.6 (228.9); and DeltaH(f)(N(2)H(5)(+)) = 187.1 (182.7). In addition, we calculated the heats of formation of CH(3)NH(2), CH(3)NNH, and CH(3)HNNHCH(3) by using isodesmic reactions and at the G3(MP2) level. The calculated results for the hydrogenation reaction RNNR + H(2) --> RHNNHR show that substitution of an organic substituent for H improved the energetics, suggesting that these types of compounds may be possible to use in a chemical hydrogen storage system.     

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.     

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.     

Enthalpies of formation of gas-phase N3, N3-, N5+, and N5- from Ab initio molecular orbital theory, stability predictions for N5(+)N3(-) and N5(+)N5(-), and experimental evidence for the instability of N5(+)N3(-)

Ab initio molecular orbital theory has been used to calculate accurate enthalpies of formation and adiabatic electron affinities or ionization potentials for N3, N3-, N5+, and N5- from total atomization energies. The calculated heats of formation of the gas-phase molecules/ions at 0 K are DeltaHf(N3(2Pi)) = 109.2, DeltaHf(N3-(1sigma+)) = 47.4, DeltaHf(N5-(1A1')) = 62.3, and DeltaHf(N5+(1A1)) = 353.3 kcal/mol with an estimated error bar of +/-1 kcal/mol. For comparison purposes, the error in the calculated bond energy for N2 is 0.72 kcal/mol. Born-Haber cycle calculations, using estimated lattice energies and the adiabatic ionization potentials of the anions and electron affinities of the cations, enable reliable stability predictions for the hypothetical N5(+)N3(-) and N5(+)N5(-) salts. The calculations show that neither salt can be stabilized and that both should decompose spontaneously into N3 radicals and N2. This conclusion was experimentally confirmed for the N5(+)N3(-) salt by low-temperature metathetical reactions between N5SbF6 and alkali metal azides in different solvents, resulting in violent reactions with spontaneous nitrogen evolution. It is emphasized that one needs to use adiabatic ionization potentials and electron affinities instead of vertical potentials and affinities for salt stability predictions when the formed radicals are not vibrationally stable. This is the case for the N5 radicals where the energy difference between vertical and adiabatic potentials amounts to about 100 kcal/mol per N5.     

Thermodynamic properties of molecular borane amines and the [BH4-][NH4+] salt for chemical hydrogen storage systems from ab initio electronic structure theory

The heats of formation for the borane amines BH3NH3, BH2NH2, and HBNH, tetrahedral BH4-, and the BN molecule have been calculated by using ab initio molecular orbital theory. Coupled cluster calculations with single and double excitations and perturbative triples (CCSD(T)) were employed for the total valence electronic energies. Correlation consistent basis sets were used, up through the augmented quadruple-zeta, to extrapolate to the complete basis set limit. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. Geometries were calculated at the CCSD(T) level up through at least aug-cc-pVTZ and frequencies were calculated at the CCSD(T)/aug-cc-pVDZ level. The heats of formation (in kcal/mol) at 0 K in the gas phase are Delta Hf(BH3NH3) = -9.1, Delta Hf(BH2NH2) = -15.9, Delta Hf(BHNH) = 13.6, Delta Hf(BN) = 146.4, and Delta Hf(BH4-) = -11.6. The reported experimental value for Delta Hf(BN) is clearly in error. The heat of formation of the salt [BH4-][NH4+](s) has been estimated by using an empirical expression for the lattice energy and the calculated heats of formation of the two component ions. The calculations show that both BH3NH3(g) and [BH4-][NH4+](s) can serve as good hydrogen storage systems which release H2 in a slightly exothermic process. The hydride affinity of BH3 is calculated to be 72.2 kcal/mol, in excellent agreement with the experimental value at 298 K of 74.2 +/- 2.8 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.     

Thermochemical properties of selenium fluorides, oxides, and oxofluorides

The bond dissociation energies (BDEs), fluoride and fluorocation affinities, and electron affinities of SeF(n) (n = 1-6), SeOF(n) (n = 0-4), and SeO(2)F(n) (n = 0-2) have been predicted with coupled cluster CCSD(T) theory extrapolated to the complete basis set limit. To achieve near chemical accuracy, additional corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies. These included corrections for core-valence effects, scalar relativistic effects, for first-order atomic spin-orbit effects, and vibrational zero point energies. The adiabatic BDEs contain contributions from product reorganization energies and, therefore, can be much smaller than the diabatic BDEs and can vary over a wide range. For thermochemical calculations, the adiabatic values must be used, whereas for bond strength and kinetic considerations, the diabatic values should be used when only small displacements of the atoms without change of the geometry of the molecule are involved. The adiabatic Se-F BDEs of SeF(n) (n = 1-6) are SeF(6) = 90, SeF(5) = 27, SeF(4) = 93, SeF(3) = 61, SeF(2) = 86, and SeF = 76 kcal/mol, and the corresponding diabatic values are SeF(6) = 90, SeF(5) = 88, SeF(4) = 93, SeF(3) = 74, SeF(2) = 86, and SeF = 76 kcal/mol. The adiabatic Se-O BDEs of SeO(n) (n = 1-3), SeOF(n) (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 23 to 107 kcal/mol, whereas the diabatic ones range from 62 to 154 kcal/mol. The adiabatic Se-F BDEs of SeOF(n) (n = 1-4) and SeO(2)F(n) (n = 1,2) range from 20 to 88 kcal/mol, whereas the diabatic ones range from 73 to 112 kcal/mol. The fluoride affinities of SeF(n), (n = 1-6), SeO(n), (n = 1-3), SeOF(n), (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 15 to 121 kcal/mol, demonstrating that the Lewis acidity of these species covers the spectrum from very weak (SeF(6)) to very strong (SeO(3)) acids. The electron affinities which are a measure of the oxidizing power of a species, span a wide range from 1.56 eV in SeF(4) to 5.16 eV in SeF(5) and for the free radicals are much higher than for the neutral molecules. Another interesting feature of these molecules and ions stems from the fact that many of them possess both a Se free valence electron pair and a free unpaired valence electron, raising the questions of their preferred location and their influence on the Se-F and Se═O bond strengths.     

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

Carboranes and baskets from reaction of B4H10 with allene

The reaction of the boron hydride B4H10 with allene was studied at the CCSD(T)/6-311+G(d)//MP2/6-31G(d) level. The mechanism is surprisingly complex with 44 transition states and several branching points located. The four carboranes and one basket that have been observed experimentally are all connected by pathways that have very similar free energies of activation. In addition, two new structures, a basket (2,4-(CH2CH2CH2)B4H8, 5a) and a "classical" structure (1,4-(Me2C)bisdiborane, 7), which might be obtained from the B4H10 + C3H4 reaction under the right conditions (hot/cold, quenched, etc.) have been identified. The first branch point in the reaction is the competition between H2 elimination from B4H10 (DeltaG(298 K) = 32.2 kcal/mol) and the hydroboration of allene by B4H10 (DeltaG(298 K) = 31.3 kcal/mol). The next branch point in the hydroboration mechanism controls the formation of 2,4-(MeCHCH2)B4H8 (1) (DeltaG(298 K) = 31.5 kcal/mol) and arachno-1,2/arachno-1,3-Me2-1-CB4H7 (8 and 8a) (DeltaG(298 K) = 34.3 kcal/mol). Another branch point in the H2-elimination mechanism controls the formation of 1-Me-2,5-micro-CH2-1-CB4H7 (29) (DeltaG(298 K) = 0.1 kcal/mol) and 2,5-micro-CHMe-1-CB4H7 (25/26) (DeltaG(298 K) = 7.3 kcal/mol). Formation of 2-Me-2,3-C2B4H7, a carborane observed in the reaction of methylacetylene with B4H10, is calculated to be blocked by a high barrier for H2 elimination. All free energies are relative to B4H10 + allene. An interesting reaction step discovered is the "reverse hydroboration step" in which a hydrogen atom is transferred from carbon back to boron, which allows a CH hydrogen to shuttle between the terminal and central carbon of allene.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

Magnitude of the CH/pi interaction in the gas phase: experimental and theoretical determination of the accurate interaction energy in benzene-methane

The accurate CH/pi interaction energy of the benzene-methane model system was experimentally and theoretically determined. In the experiment, mass analyzed threshold ionization spectroscopy was applied to the benzene-methane cluster in the gas phase, prepared in a supersonic molecular beam. The binding energy in the neutral ground state of the cluster, which is regarded as the CH/pi interaction energy for this model system, was evaluated from the dissociation threshold measurements of the cluster cation. The experimentally determined binding energy (D(0)) was 1.03-1.13 kcal/mol. The interaction energy of the model system was calculated by ab initio molecular orbital methods. The estimated CCSD(T) interaction energy at the basis set limit (D(e)) was -1.43 kcal/mol. The calculated binding energy (D(0)) after the vibrational zero-point energy correction (1.13 kcal/mol) agrees well with the experimental value. The effects of basis set and electron correlation correction procedure on the calculated CH/pi interaction energy were evaluated. Accuracy of the calculated interaction energies by DFT methods using BLYP, B3LYP, PW91 and PBE functionals was also discussed.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.     

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.     

Thermodynamic properties of the C5, C6, and C8 n-alkanes from ab initio electronic structure theory

The heats of formation for the n-alkanes C(n)H(n+2) for n = 5, 6, and 8 have been calculated using ab initio molecular orbital theory. Coupled-cluster calculations with perturbative triples (CCSD(T)) were employed for the total valence electronic energies. Correlation-consistent basis sets were used, up through the augmented quadruple zeta, to extrapolate to the complete basis set limit. Geometries were optimized at the B3LYP/TZVP and MP2/aug-cc-pVTZ levels. The MP2 geometries were used in the CCSD(T) calculations. Frequencies were determined at the density functional level (B3LYP/TZVP), and scaled zero point energies were calculated from the B3LYP frequencies. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. The core/valence corrections are not small, (approximately 1.1 kcal/mol per carbon unit) and cannot be neglected for chemical accuracy. The calculated deltaH(298)f values are -35.0, -40.2, and -50.2 kcal/mol for C5H12, C6H14, and C8H18, respectively, in excellent agreement with the respective experimental values of -35.11 +/- 0.19, -39.89 +/- 0.19, and -49.90 +/- 0.31 kcal/mol. Isodesmic reaction energies are presented for some simple reactions involving C8H18 and are shown not to be strongly method dependent.     

A photoelectron and TPEPICO investigation of the acetone radical cation

The valence shell photoelectron spectrum, threshold photoelectron spectrum, and threshold photoelectron photoion coincidence (TPEPICO) mass spectra of acetone have been measured using synchrotron radiation. New vibrational progressions have been observed and assigned in the X 2B2 state photoelectron bands of acetone-h6 and acetone-d6, and the influence of resonant autoionization on the threshold electron yield has been investigated. The dissociation thresholds for fragment ions up to 31 eV have been measured and compared to previous values. In addition, kinetic modeling of the threshold region for CH3* and CH4 loss leads to new values of 78 +/- 2 kJ mol(-1) and 75 +/- 2 kJ mol(-1), respectively, for the 0 K activation energies for these two processes. The result for the methyl loss channel is in reasonable agreement with, but slightly lower than, that of 83 +/- 1 kJ mol(-1) derived in a recent TPEPICO study by Fogleman et al. The modeling accounts for both low-energy dissociation channels at two different ion residence times in the mass spectrometer. Moreover, the effects of the ro-vibrational population distribution, the electron transmission efficiency, and the monochromator band-pass are included. The present activation energies yield a Delta(f)H298 for CH3CO+ of 655 +/- 3 kJ mol(-1), which is 4 kJ mol(-1) lower than that reported by Fogleman et al. The present Delta(f)H298 for CH3CO+ can be combined with the Delta(f)H298 for CH2CO (-47.5 +/- 1.6 kJ mol(-1)) and H+ (1530 kJ mol(-1)) to yield a 298 K proton affinity for ketene of 828 +/- 4 kJ mol(-1), in good agreement with the value (825 kJ mol(-1)) calculated at the G2 level of theory. The measured activation energy for CH4 loss leads to a Delta(f)H298 (CH2CO+*) of 873 +/- 3 kJ mol(-1).