Heats of formation of boron hydride anions and dianions and their ammonium salts [BnHmy-][NH4+]y with y=1-2

Thermochemical parameters of the closo boron hydride BnHn2- dianions, with n=5-12, the B3H8- and B11H14- anions, and the B5H9 and B10H14 neutral species were predicted by high-level ab initio electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T)/complete basis set (CBS) extrapolated energies, plus additional corrections were used to predict the heats of formation of the simplest BnHmy- species in the gas phase in kcal/mol at 298 K: DeltaHf(B3H8-)=-23.1+/-1.0; DeltaHf(B5H52-)=119.4+/-1.5; DeltaHf(B6H62-)=64.1+/-1.5; and DeltaHf(B5H9)=24.1+/-1.5. The heats of formation of the larger species were evaluated by the G3 method from hydrogenation reactions (values at 298 K, in kcal/mol with estimated error bars of+/-3 kcal/mol): DeltaHf(B7H72-)=51.8; DeltaHf(B8H82-)=46.1; DeltaHf(B9H92-)=24.4; DeltaHf(B10H102-)=-12.5; DeltaHf(B11H112-)=-11.8; DeltaHf(B12H122-)=-86.3; DeltaHf(B11H14-)=-57.3; and DeltaHf(B10H14)=18.7. A linear correlation between atomization energies of the dianions and energies of the BH units was found. The heats of formation of the ammonium salts of the anions and dianions were predicted using lattice energies (UL) calculated from an empirical expression based on ionic volumes. The UL values (0 K) of the BnHn2- dianions range from 319 to 372 kcal/mol. The values of UL for the B3H8- and B11H14- anions are 113 and 135 kcal/mol, respectively. The calculated lattice energies and gas-phase heats of formation of the constituent ions were used to predict the heats of formation of the ammonium crystal salts [BnHmy-][NH4+]y. These results were used to evaluate the thermodynamics of the H2 release reactions from the ammonium hydro-borate salts.     

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.     

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.     

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.     

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.     

Thermochemistry of new molecular species: SBr and HSBr

Total energies, optimized geometries, and vibrational frequencies of SBr and HSBr have been evaluated at the coupled cluster level of theory with the correlation consistent basis sets. Extrapolated to the complete basis set limit and with corrections for core-valence, scalar relativistic, and spin-orbit effects, atomization energies were computed and then combined with the experimental heats of formation of the atomic species to generate very accurate heats of formation for the species SBr and HSBr. For SBr, we predict 37.45 and 36.07 kcal/mol for DeltaHf(0 K) and DeltaHf (298.15 K), respectively, in very good agreement with the inferred experimental values of 37.98 and 36.15 kcal/mol. For HSBr, the estimate turns out to be 10.38 and 8.29 kcal/mol for DeltaHf (0 K) and DeltaHf (298.15 K), respectively. Using the more recent HBrO experimental heat of formation at 298.15 K of Lock et al., [J. Phys. Chem. 100, 7972 (1996)] the inferred experimental value for HSBr is predicted to be 8.15 kcal/mol, compared with 8.65 kcal/mol derived from the data of Ruscic and Berkowitz [J. Chem. Phys. 101, 7795 (1994)]. Considering the better agreement of the result with that predicted using the experimental value of DeltaHf(298.15 K) of Lock et al., the author also supports the suggestion made by Denis [J. Phys. Chem. A. 110, 5887 (2006)] that the result of Lock et al. should be preferred over the one of Ruscic and Berkowitz. For DeltaHf(0 K), the author found 10.38 and 10.56 kcal/mol, respectively, for the theoretical and inferred experimental estimates.     

Heats of formation of diphosphene, phosphinophosphinidene, diphosphine, and their methyl derivatives, and mechanism of the borane-assisted hydrogen release

The heats of formation of diphosphene (cis- and trans-P2H2), phopshinophosphinidene (singlet and triplet H2PP) and diphosphine (P2H4), as well as those of the P2H and P2H3 radicals resulting from PH bond cleavages, have been calculated by using high-level ab initio electronic structure theory. Energies were calculated using coupled-cluster theory with a perturbative treatment for triple excitations (CCSD(T)) and employing augmented correlation consistent basis sets with additional tight d-functions on P (aug-cc-pV(n+d)Z) up to quadruple- or quintuple-zeta, to perform a complete basis set extrapolation for the energy. Geometries and vibrational frequencies were determined with the CCSD(T) method. 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 298 [0] K: DeltaH(degree)(f)(P2H) = 53.4 [54.4]; DeltaH(degree)(f)(cis-P2H2) = 32.0 [33.9]; DeltaH(degree)(f)(trans-P2H2) = 28.7 [30.6]; DeltaH(degree)(f)(H2PP) = 53.7 [55.6]; DeltaH(degree)(f)(3H2PP) = 56.5 [58.3]; DeltaH(degree)(f)(P2H3) = 32.3 [34.8]; DeltaH(degree)(f)(P2H4) = 5.7 [9.1] (expt, 5.0 +/- 1.0 at 298 K); and DeltaH(degree)(f)(CH3PH2) = -5.0 [-1.4]. We estimate these values to have an accuracy of +/-1.0 kcal/mol. In contrast to earlier results, we found a singlet ground state for phosphinophosphinidene (H2PP) with a singlet-triplet energy gap of 2.8 kcal/mol. We calculated the heats of formation of the methylated derivatives CH3PPH, CH3HPPH2, CH3PPCH3, CH3HPP, (CH3)2PP, (CH3)2PPH2, and CH3HPPHCH3 by using isodesmic reactions at the MP2/CBS level. The calculated results for the hydrogenation reactions RPPR + H2 --> RHPPHR and R2PP + H2 --> R2PPH2 show that substitution of an organic substituent for H improves the energetics, suggesting that secondary diphosphines and diphosphenes are potential candidates for use in a chemical hydrogen storage system. A comparison with the nitrogen analogues is given. The mechanism for H2-generation from diphosphine without and with BH3 as a catalyst was examined. Including tunneling corrections, the rate constant for the catalyzed reaction is 4.5 x 1015 times faster than the uncatalyzed result starting from separated catalyst and PH2PH2.     

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.     

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.     

The ozonolysis of acetylene--a quantum chemical investigation.

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

Trimethoxybenzene complexes of pentafluorophenylchlorocarbene

Pentafluorophenylchlorocarbene, generated by laser flash photolysis (LFP) of pentafluorophenylchlorodiazirine, formed π-type complexes with 1,3,5-trimethoxybenzene in pentane. The carbene and carbene complexes were in equilibrium with K = 3.21 × 10(5) M(-1) at 294 K. From the temperature dependence of K, ΔH° = -10.2 kcal/mol, ΔS° = -9.5 eu, and ΔG° = -7.4 kcal/mol at 298 K. The carbene complexes were characterized by UV-vis spectroscopy and computational analysis.     

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.     

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.     

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.     

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

Cheletropic decomposition of cyclic nitrosoamines revisited: the nature of the transition states and a critical role of the ring strain

The cheletropic decompositions of 1-nitrosoaziridine (1), 1-nitroso-Delta(3)-pyrroline (2), 7-nitroso-7-azabicyclo[2.2. 1]hepta-2,5-diene (3), and 6-nitroso-6-azabicyclo[2.1.1]hexa-4-ene (4) have been studied theoretically using high level ab initio computations. Activation parameters of the decomposition of nitrosoaziridine 1 were obtained experimentally in heptane (DeltaH()(298) = 18.6 kcal mol(-)(1), DeltaS()(298) = -7.6 cal mol(-)(1) K(-)(1)) and methanol (20.3 kcal mol(-)(1), 0.3 cal mol(-)(1) K(-)(1)). Among employed theoretical methods (B3LYP, MP2, CCD, CCSD(T)//CCD), the B3LYP method in conjunction with 6-31+G, 6-311+G, and 6-311++G(3df,2pd) basis sets gives the best agreement with experimental data. It was found that typical N-nitrosoheterocycles 2-4 which have high N-N bond rotation barriers (>16 kcal mol(-)(1)) extrude nitrous oxide via a highly asynchronous transition state with a planar ring nitrogen atom. Nitrosoaziridine 1, with a low rotation barrier (<9 kcal mol(-)(1)) represents a special case. This compound can eliminate N(2)O via a low energy linear synperiplanar transition state (DeltaH()(298) = 20.6 kcal mol(-)(1), DeltaS()(298) = 2.5 cal mol(-)(1) K(-)(1)). Two higher energy transition states are also available. The B3LYP activation barriers of the cheletropic fragmentation of nitrosoheterocycles 2-4 decrease in the series: 2 (58 kcal mol(-)(1)) > 3 (18 kcal mol(-)(1)) > 4 (12) kcal mol(-)(1). The relative strain energies increase in the same order: 2 (0 kcal mol(-)(1)) < 3 (39 kcal mol(-)(1)) < 4 (52 kcal mol(-)(1)). Comparison of the relative energies of 2-4 and their transition states on a common scale where the energy of nitrosopyrroline 2 is assumed as reference indicates that the thermal stability of the cyclic nitrosoamines toward cheletropic decomposition is almost entirely determined by the ring strain.     

Thermodynamic properties of molecular borane phosphines, alane amines, and phosphine alanes and the [BH(4)(-)][PH(4)(+)], [AlH(4)(-)][NH(4)(+)], and [AlH(4)(-)][PH(4)(+)] salts for chemical hydrogen storage systems from ab initio electronic structure theory

The heats of formation for the molecules BH(3)PH(3), BH(2)PH(2), HBPH, AlH(3)NH(3), AlH(2)NH(2), HAlNH, AlH(3)PH(3), AlH(2)PH(2), HAlPH, AlH(4)(-), PH(3), PH(4), and PH(4)(+), as well as the diatomics BP, AlN, and AlP, have been calculated by using ab initio molecular orbital theory. The coupled cluster with single and double excitations and perturbative triples method (CCSD(T)) was 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. Additional d core functions were used for Al and P. 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 of the salts [BH(4)(-)][PH(4)(+)](s), [AlH(4)(-)][NH(4)(+)](s), and [AlH(4)(-)][PH(4)(+)](s) have 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 AlH(3)NH(3)(g) and [AlH(4)(-)][NH(4)(+)](s) can serve as good hydrogen storage systems that release H(2) in a slightly exothermic process. In addition, AlH(3)PH(3) and the salts [AlH(4)(-)][PH(4)(+)] and [BH(4)(-)][PH(4)(+)] have the potential to serve as H(2) storage systems. The hydride affinity of AlH(3) is calculated to be -70.4 kcal/mol at 298 K. The proton affinity of PH(3) is calculated to be 187.8 kcal/mol at 298 K in excellent agreement with the experimental value of 188 kcal/mol. PH(4) is calculated to be barely stable with respect to loss of a hydrogen to form PH(3).     

Accurate and efficient method for predicting thermochemistry of polycyclic aromatic hydrocarbons - bond-centered group additivity

A self-consistent estimation method for the thermochemical properties of polycyclic aromatic hydrocarbons (PAH) is presented. This method is based on enthalpies of formation (DeltaHf(degrees), entropies (S(degrees)298, and heat capacities (C(degrees)p obtained from B3LYP/6-31G(d) calculations of the total energies and frequencies for 139 PAHs, including C(60) and C(70) fullerenes. The enthalpies of formation were calculated using an optimized set of homodesmic reactions given the available experimental DeltaHf(degrees) of PAHs. The theoretical entropies were compared with the existing experimental entropies, and some inconsistencies in the experimental data were identified. The estimation method presented here is a systematic extension of the widely employed atom-centered group additivity method originally proposed by Benson. This new method is based on bond-centered groups that define bonds linking two atom-centered groups and specify the size of the rings to which they belong. In addition, a term to describe the resonance energy is included. The thermochemical properties of PAHs up to C(70) fullerene are estimated with a mean average deviation of 2.8 kcal mol(-1) in DeltaHf(degrees), 0.7 cal K(-1) mol(-1) in S(degrees)298, and about 0.5 cal K(-1) mol(-1) in the C(degrees)p. This bond-centered group additivity method for the thermochemical properties of PAHs significantly expands both the range of systems that can be estimated and the accuracy of the estimations. The results of this work also allow us to assess the quality of available experimental data. For example, there are strong indications that the literature DeltaHf(degrees)of benzo[k]fluoranthene is about 10 kcal mol(-1) too low.     

The enthalpies of formation of o-, m-, and p-benzoquinone: gas-phase ion energetics, combustion calorimetry, and quantum chemical computations combined

Radical anions of o-, m-, and p-benzoquinone were produced in a Fourier transform mass spectrometer by low energy electron attachment or collision-induced dissociation and were differentiated. Classical derivatization experiments also were carried out to authenticate the ortho and meta anions. Gas-phase techniques were used to measure the proton affinities of all three radical anions and the electron affinities of o- and m-benzoquinone. By combining these results in thermodynamic cycles, we derived heats of hydrogenation of o-, m-, and p-benzoquinone (Delta(hyd)H degrees (1o, 1m, and 1p) = 42.8 +/- 4.1, 74.8 +/- 4.1, and 38.5 +/- 3.0 kcal mol(-)(1), respectively) and their heats of formation (Delta(f)H degrees (1o, 1m, and 1p) = -23.1 +/- 4.1, 6.8 +/- 4.1, and -27.7 +/- 3.0 kcal mol(-)(1), respectively). Good accord with the literature value for the para derivative was obtained. Combustion calorimetry and heats of sublimation also were measured for benzil and 3,5-di-tert-butyl-o-benzoquinone. The former heat of formation agreed with previous determinations, while the latter result (Delta(f)H degrees (g) = -73.09 +/- 0.87 kcal mol(-)(1)) was transformed to Delta(f)H degrees (1o) = -18.9 +/- 2.2 kcal mol(-)(1) by removing the effect of the tert-butyl groups via isodesmic reactions. This led to a final value of Delta(f)H degrees (1o) = -21.0 +/- 3.1 kcal mol(-)(1). Additivity was found to work well for m-benzoquinone, but BDE1 and BDE2 for 1,2- and 1,4-dihydroxybenzene differed by a remarkably small 14.1 +/- 4.2 and 23.5 +/- 3.7 kcal mol(-)(1), respectively, indicating that o- and p-benzoquinone should be excellent radical traps.     

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.