X-ray Crystal Structure of Thorium Tetrahydroborate, Th(BH4)4, and Computational Studies of An(BH4)4 (An = Th, U).

The crystal structure of Th(BH4)4 is described. Two of the four BH4- ions are terminal and tridentate (κ3), whereas the other two bridge between neighboring ThIV centers in a κ2,κ2 (i.e., bis-bidentate) fashion. Thus, each thorium center is bound to six BH4- groups by 14 Th-H bonds. The six boron atoms describe a distorted octahedron in which the κ3-BH4- ions are mutually cis; the 14 ligating hydrogen atoms define a highly distorted bicapped hexagonal antiprism. The thorium centers are linked into a polymer consisting of interconnected helical chains wound about 4-fold screw axes. The structures of An(BH4)4 (An = Th, U) were also investigated by DFT. The geometries of [An(BH4)6]2-, [An3(BH4)16]4-, and [An5(BH4)26]6- fragments of the polymeric structures were optimized at the B3LYP and/or PBE levels. Most calculated geometries are 14-coordinate and agree with the experimental structures, but isolated [Th(BH4)6]2- units are predicted to feature 16-coordinate ThIV centers.

Infrared Signature of the Cation-π Interaction between Calcite and Aromatic Hydrocarbons.

The cation-π interaction is proposed as an important mechanism for the adsorption of aromatic hydrocarbons having non-zero quadrupole moments by mineral surfaces. Direct evidence supporting such a mechanism is, however, limited. Using the model mineral calcite, we probe the cation-π interaction with adsorbed benzene, toluene, and ethylbenzene (BTE) molecules using attenuated total reflectance Fourier transform infrared spectroscopy. We show that the presence of calcite increases the energy required to excite the synchronized bending of aromatic C-H bonds of BTE molecules. The unique conformation of this vibrational mode indicates that the planar aromatic rings of BTE molecules are constrained in a tilted face-down position by the cation-π interaction, as further confirmed by density functional theory calculations. Our results suggest that the shift of the excitation energy of the aromatic C-H bending may be used as an infrared signature for the cation-π interaction occurring on mineral surfaces.

A trigonal-prismatic hexanuclear technetium(II) bromide cluster: solid-state synthesis and crystallographic and electronic structure.

The compound Na{[Tc6Br12]2Br} has been obtained from the decomposition of TcBr4 under vacuum in a Pyrex ampule at 450 °C. The stoichiometry of the compound has been confirmed by energy-dispersive X-ray spectroscopy and its structure determined by single-crystal X-ray diffraction. The compound contains a trigonal-prismatic hexanuclear [Tc6Br12] cluster. The cluster is composed of two triangular Tc3Br6 units linked by multiple Tc-Tc bonds. In the Tc3Br6 unit, the average Tc-Tc distance [2.6845(5) Å] is characteristic of Tc-Tc single bonds, while the average Tc-Tc distance between the two triangular units [2.1735(5) Å] is characteristic of Tc≡Tc triple bonds. The electronic structure of the [Tc6Br12] cluster was studied by first-principles calculations, which confirm the presence of single and triple Tc-Tc bonds in the cluster.

Differentiating between trivalent lanthanides and actinides.

The reactions of LnCl(3) with molten boric acid result in the formation of Ln[B(4)O(6)(OH)(2)Cl] (Ln = La-Nd), Ln(4)[B(18)O(25)(OH)(13)Cl(3)] (Ln = Sm, Eu), or Ln[B(6)O(9)(OH)(3)] (Ln = Y, Eu-Lu). The reactions of AnCl(3) (An = Pu, Am, Cm) with molten boric acid under the same conditions yield Pu[B(4)O(6)(OH)(2)Cl] and Pu(2)[B(13)O(19)(OH)(5)Cl(2)(H(2)O)(3)], Am[B(9)O(13)(OH)(4)]·H(2)O, or Cm(2)[B(14)O(20)(OH)(7)(H(2)O)(2)Cl]. These compounds possess three-dimensional network structures where rare earth borate layers are joined together by BO(3) and/or BO(4) groups. There is a shift from 10-coordinate Ln(3+) and An(3+) cations with capped triangular cupola geometries for the early members of both series to 9-coordinate hula-hoop geometries for the later elements. Cm(3+) is anomalous in that it contains both 9- and 10-coordinate metal ions. Despite these materials being synthesized under identical conditions, the two series do not parallel one another. Electronic structure calculations with multireference, CASSCF, and density functional theory (DFT) methods reveal the An 5f orbitals to be localized and predominately uninvolved in bonding. For the Pu(III) borates, a Pu 6p orbital is observed with delocalized electron density on basal oxygen atoms contrasting the Am(III) and Cm(III) borates, where a basal O 2p orbital delocalizes to the An 6d orbital. The electronic structure of the Ce(III) borate is similar to the Pu(III) complexes in that the Ce 4f orbital is localized and noninteracting, but the Ce 5p orbital shows no interaction with the coordinating ligands. Natural bond orbital and natural population analyses at the DFT level illustrate distinctive larger Pu 5f atomic occupancy relative to Am and Cm 5f, as well as unique involvement and occupancy of the An 6d orbitals.

Periodic trends in lanthanide and actinide phosphonates: discontinuity between plutonium and americium.

The hydrothermal reactions of trivalent lanthanide and actinide chlorides with 1,2-methylenediphosphonic acid (C1P2) in the presence of NaOH or NaNO(3) result in the crystallization of three structure types: RE[CH(2)(PO(3)H(0.5))(2)] (RE = La, Ce, Pr, Nd, Sm; Pu) (A type), NaRE(H(2)O)[CH(2)(PO(3))(2)] (RE = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy; Am) (B type), or NaLn[CH(2)(PO(3)H(0.5))(2)]·(H(2)O) (Ln = Yb and Lu) (C type). These crystals were analyzed using single crystal X-ray diffraction, and the structures were used directly for detailed bonding calculations. These phases form three-dimensional frameworks. In both A and B, the metal centers are found in REO(8) polyhedra as parts of edge-sharing chains or edge-sharing dimers, respectively. Polyhedron shape calculations reveal that A favors a D(2d) dodecahedron while B adopts a C(2v) geometry. In C, Yb and Lu only form isolated MO(6) octahedra. Such differences in terms of structure topology and coordination geometry are discussed in detail to reveal periodic deviations between the lanthanide and actinide series. Absorption spectra for the Pu(III) and Am(III) compounds are also reported. Electronic structure calculations with multireference methods, CASSCF, and density functional theory, DFT, reveal localization of the An 5f orbitals, but natural bond orbital and natural population analyses at the DFT level illustrate unique occupancy of the An 6d orbitals, as well as larger occupancy of the Pu 5f orbitals compared to the Am 5f orbitals.

Synthesis of a uranyl persulfide complex and quantum chemical studies of formation and topologies of hypothetical uranyl persulfide cage clusters.

The compound Na(4)[(UO(2))(S(2))(3)](CH(3)OH)(8) was synthesized at room temperature in an oxygen-free environment. It contains a rare example of the [(UO(2))(S(2))(3)](4-) complex in which a uranyl ion is coordinated by three bidentate persulfide groups. We examined the possible linkage of these units to form nanoscale cage clusters analogous to those formed from uranyl peroxide polyhedra. Quantum chemical calculations at the density functional and multiconfigurational wave function levels show that the uranyl-persulfide-uranyl, U-(S(2))-U, dihedral angles of model clusters are bent due to partial covalent interactions. We propose that this bent interaction will favor assembly of uranyl ions through persulfide bridges into curved structures, potentially similar to the family of nanoscale cage clusters built from uranyl peroxide polyhedra. However, the U-(S(2))-U dihedral angles predicted for several model structures may be too tight for them to self-assemble into cage clusters with fullerene topologies in the absence of other uranyl-ion bridges that adopt a flatter configuration. Assembly of species such as [(UO(2))(S(2))(SH)(4)](4-) or [(UO(2))(S(2))(C(2)O(4))(4)](4-) into fullerene topologies with ~60 vertices may be favored by use of large counterions.

Elucidation of tetraboric acid with a new borate fundamental building block in a chiral uranyl fluoroborate.

A new neutral borate species, H(2)B(4)O(7) (also known as tetraboric acid), with a one-dimensional chain structure, is found in the interlayer spacing in Rb(2)[(UO(2))(2)B(8)O(12)F(6)]·H(2)B(4)O(7) (RbUBOF-2) derived from boric acid flux reaction of uranyl(VI) nitrate with RbBF(4). This new form of tetraboric acid possesses a novel borate fundamental building block with the symbol 4Δ:<3Δ>Δ.

Uranium and thorium hydride complexes as multielectron reductants: a combined neutron diffraction and quantum chemical study.

The unusual uranium reaction system in which uranium(4+) and uranium(3+) hydrides interconvert by formal bimetallic reductive elimination and oxidative addition reactions, [(C(5)Me(5))(2)UH(2)](2) (1) ⇌ [(C(5)Me(5))(2)UH](2) (2) + H(2), was studied by employing multiconfigurational quantum chemical and density functional theory methods. 1 can act as a formal four-electron reductant, releasing H(2) gas as the byproduct of four H(2)/H(-) redox couples. The calculated structures for both reactants and products are in good agreement with the X-ray diffraction data on 2 and 1 and the neutron diffraction data on 1 obtained under H(2) pressure as part of this study. The interconversion of the uranium(4+) and uranium(3+) hydride species was calculated to be near thermoneutral (~-2 kcal/mol). Comparison with the unknown thorium analogue, [(C(5)Me(5))(2)ThH](2), shows that the thorium(4+) to thorium(3+) hydride interconversion reaction is endothermic by 26 kcal/mol.

Heats of formation of MHxCly (M = Si, P, As, Sb) compounds and main group fluorides from high level electronic structure calculations.

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for the MH(x)Cl(y) compounds (M = Si, P, As, and Sb) and for a number of trivalent, tetravalent, and pentavalent fluorides (SbF(3), BiF(3), GeF(4), SnF(4), PbF(4), AsF(5), SbF(5)) from coupled cluster theory (CCSD(T)) calculations using correlation consistent basis sets and extrapolation to the complete basis set limit. Small-core, relativistic effective core potentials were used for the heavier elements (Ge, As, Sn, Sb, Pb, and Bi), including correlation of the outer core electrons. Additional scalar relativistic (for the lighter elements) and atomic spin-orbit corrections are included in order to achieve near chemical accuracy of ±1.5 kcal/mol. Vibrational zero point energies were computed from scaled harmonic frequencies at the second order Møller-Plesset perturbation theory (MP2) level where possible. Agreement between theory and the available experimental data is excellent. We present a revised heat of formation of the antimony atom in the gas phase. The calculated values will be of use in predicting the behavior of chemical vapor deposition systems.

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

Hydrogen storage by boron-nitrogen heterocycles: a simple route for spent fuel regeneration.

We describe a new hydrogen storage platform based on well-defined BN heterocyle materials. Specifically, we demonstrate that regeneration of the spent fuel back to the charged fuel can be accomplished using molecular H(2) and H(-)/H(+) sources. Crystallographic characterization of intermediates along the regeneration pathway confirms our structural assignments and reveals unique bonding changes associated with increasing hydrogen content on boron and nitrogen. Synthetic access to the fully charged BN cyclohexane fuels will now enable investigations of these materials in hydrogen desorption studies.

Resonance stabilization energy of 1,2-azaborines: a quantitative experimental study by reaction calorimetry.

Aromatic and single-olefin six-membered BN heterocycles were synthesized, and the heats of hydrogenation were measured calorimetrically. A comparison of the hydrogenation enthalpies of these compounds revealed that 1,2-azaborines have a resonance stabilization energy of 16.6 ± 1.3 kcal/mol, in good agreement with calculated values.

Heats of formation of XeF(3)(+), XeF(3)(-), XeF(5)(+), XeF(7)(+), XeF(7)(-), and XeF(8) from high level electronic structure calculations.

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(3)(+), XeF(3)(-), XeF(5)(+), XeF(7)(+), XeF(7)(-), and XeF(8) from coupled cluster theory (CCSD(T)) calculations with effective core potential correlation-consistent basis sets for Xe and including correlation of the nearest core electrons. Additional corrections are included to achieve near chemical accuracy of +/-1 kcal/mol. Vibrational zero point energies were computed at the MP2 level of theory. Unlike the other neutral xenon fluorides, XeF(8) is predicted to be thermodynamically unstable with respect to loss of F(2) with the reaction calculated to be exothermic by 22.3 kcal/mol at 0 K. XeF(7)(+) is also predicted to be thermodynamically unstable with respect to the loss of F(2) by 24.1 kcal/mol at 0 K. For XeF(3)(+), XeF(5)(+), XeF(3)(-), XeF(5)(-), and XeF(7)(-), the reactions for loss of F(2) are endothermic by 14.8, 37.8, 38.2, 59.6, and 31.9 kcal/mol at 0 K, respectively. The F(+) affinities of Xe, XeF(2), XeF(4), and XeF(6) are predicted to be 165.1, 155.3, 172.7, and 132.5 kcal/mol, and the corresponding F(-) affinities are 6.3, 19.9, 59.1, and 75.0 kcal/mol at 0 K, respectively.

Dinitrogen difluoride chemistry. improved syntheses of cis- and trans-N2F2, synthesis and characterization of N2F(+)Sn2F9(-), ordered crystal structure of N2F(+)Sb2F11(-), high-level electronic structure calculations of cis-N2F2, trans-N2F2, F2N=N, and N2F(+), and Mechanism of the trans-cis isomerization of N2F2.

N(2)F(+) salts are important precursors in the synthesis of N(5)(+) compounds, and better methods are reported for their larger scale production. A new, marginally stable N(2)F(+) salt, N(2)F(+)Sn(2)F(9)(-), was prepared and characterized. An ordered crystal structure was obtained for N(2)F(+)Sb(2)F(11)(-), resulting in the first observation of individual N[triple bond]N and N-F bond distances for N(2)F(+) in the solid phase. The observed N[triple bond]N and N-F bond distances of 1.089(9) and 1.257(8) A, respectively, are among the shortest experimentally observed N-N and N-F bonds. High-level electronic structure calculations at the CCSD(T) level with correlation-consistent basis sets extrapolated to the complete basis limit show that cis-N(2)F(2) is more stable than trans-N(2)F(2) by 1.4 kcal/mol at 298 K. The calculations also demonstrate that the lowest uncatalyzed pathway for the trans-cis isomerization of N(2)F(2) has a barrier of 60 kcal/mol and involves rotation about the N=N double bond. This barrier is substantially higher than the energy required for the dissociation of N(2)F(2) to N(2) and 2 F. Therefore, some of the N(2)F(2) dissociates before undergoing an uncatalyzed isomerization, with some of the dissociation products probably catalyzing the isomerization. Furthermore, it is shown that the trans-cis isomerization of N(2)F(2) is catalyzed by strong Lewis acids, involves a planar transition state of symmetry C(s), and yields a 9:1 equilibrium mixture of cis-N(2)F(2) and trans-N(2)F(2). Explanations are given for the increased reactivity of cis-N(2)F(2) with Lewis acids and the exclusive formation of cis-N(2)F(2) in the reaction of N(2)F(+) with F(-). The geometry and vibrational frequencies of the F(2)N=N isomer have also been calculated and imply strong contributions from ionic N(2)F(+) F(-) resonance structures, similar to those in F(3)NO and FNO.

Thermochemistry and electronic structure of small boron clusters (B(n), n = 5-13) and their anions.

Thermochemical parameters of a set of small-sized neutral (B(n)) and anionic (B(n)(-)) boron clusters, with n = 5-13, were determined using coupled-cluster theory CCSD(T) calculations with the aug-cc-pVnZ (n = D, T, and Q) basis sets extrapolated to the complete basis set limit (CBS) plus addition corrections and/or G3B3 calculations. Enthalpies of formation, adiabatic electron affinities (EA), vertical (VDE), and adiabatic (ADE) detachment energies were evaluated. Our calculated EAs are in good agreement with recent experiments (values in eV): B(5) (CBS, 2.29; G3B3, 2.48; exptl., 2.33 +/- 0.02), B(6) (CBS, 2.59; G3B3, 3.23; exptl., 3.01 +/- 0.04), B(7) (CBS, 2.62; G3B3, 2.67; exptl., 2.55 +/- 0.05), B(8) (CBS, 3.02; G3B3, 3.11; exptl., 3.02 +/- 0.02), B(9) (G3B3, 3.03; exptl., 3.39 +/- 0.06), B(10) (G3B3, 2.85; exptl., 2.88 +/- 0.09), B(11) (G3B4, 3.48;, exptl., 3.43 +/- 0.01), B(12) (G3B3, 2.33; exptl., 2.21 +/- 0.04), and B(13) (G3B3, 3.62; exptl., 3.78 +/- 0.02). The difference between the calculated adiabatic electron affinity and the adiabatic detachment energy for B(6) is due to the fact that the geometry of the anion is not that of the ground-state neutral. The calculated adiabatic detachment energies to the (3)A(u), C(2h) and (1)A(g), D(2h) excited states of B(6), which have geometries similar to the (1)A(g), D(2h) state of B(6)(-), are 2.93 and 3.06 eV, in excellent agreement with experiment. The VDEs were also well reproduced by the calculations. Partitioning of the electron localization functions into pi and sigma components allows probing of the partial and local delocalization in global nonaromatic systems. The larger clusters appear to exhibit multiple aromaticity. The binding energies per atom vary in a parallel manner for both neutral and anionic series and approach the experimental value for the heat of atomization of B. The resonance energies and the normalized resonance energies are convenient indices to quantify the stabilization of a cluster of elements.

Thermodynamic properties of the XO(2), X(2)O, XYO, X(2)O(2), and XYO(2) (X, Y = Cl, Br, and I) isomers.

High level ab initio electronic structure calculations at the coupled cluster level with a correction for triples extrapolated to the complete basis set limit have been made for the thermodynamics of the BrBrO(2), IIO(2), ClBrO(2), ClIO(2), and BrIO(2) isomers, as well as various molecules involved in the bond dissociation processes. Of the BrBrO(2) isomers, BrOOBr is predicted to be the most stable by 8.5 and 9.3 kcal/mol compared to BrBrO(2) and BrOBrO at 298 K, respectively. The weakest bond in BrOOBr is the O-Br bond with a bond dissociation energy (BDE) of 15.9 kcal/mol, and in BrBrO(2), it is the Br-Br bond of 19.1 kcal/mol. The smallest BDE in BrOBrO is for the central O-Br bond with a BDE of 12.6 kcal/mol. Of the IIO(2) isomers, IIO(2) is predicted to be the most stable by 3.3, 9.4, and 28.9 kcal/mol compared to IOIO, IOOI, and OIIO at 298 K, respectively. The weakest bond in IIO(2) is the I-I bond with a BDE of 22.2 kcal/mol. The smallest BDEs in IOIO and IOOI are the terminal O-I bonds with values of 19.0 and 5.2 kcal/mol, respectively.

Lewis acidities and hydride, fluoride, and X- affinities of the BH(3-n)Xn compounds for (X = F, Cl, Br, I, NH2, OH, and SH) from coupled cluster theory.

Atomization energies at 0 K and enthalpies of formation at 0 and 298 K are predicted for the BH(4-n)X(n)(-) and the BH(3-n)X(n)F(-) compounds for (X = F, Cl, Br, I, NH(2), OH, and SH) from coupled cluster theory (CCSD(T)) calculations with correlation-consistent basis sets and with an effective core potential on I. To achieve near chemical accuracy (+/-1.0 kcal/mol), additional corrections were added to the complete basis set binding energies. The hydride, fluoride, and X(-) affinities of the BH(3-n)X(n) compounds were predicted. Although the hydride and fluoride affinities differ somewhat in their magnitudes, they show very similar trends and are both suitable for judging the Lewis acidities of compounds. The only significant differences in their acidity strength orders are found for the boranes substituted with the strongly electron withdrawing and back-donating fluorine and hydroxyl ligands. The highest H(-) and F(-) affinities are found for BI(3) and the lowest ones for B(NH(2))(3). Within the boron trihalide series, the Lewis acidity increases monotonically with increasing atomic weight of the halogen, that is, BI(3) is a considerably stronger Lewis acid than BF(3). For the X(-) affinities in the BX(3), HBX(2), and H(2)BX series, the fluorides show the highest values, whereas the amino and mercapto compounds show the lowest ones. Hydride and fluoride affinities of the BH(3-n)X(n) compounds exhibit linear correlations with the proton affinity of X(-) for most X ligands. Reasons for the correlation are discussed. A detailed analysis of the individual contributions to the Lewis acidities of these substituted boranes shows that the dominant effect in the magnitude of the acidity is the strength of the BX(3)(-)-F bond. The main contributor to the relative differences in the Lewis acidities of BX(3) for X, a halogen, is the electron affinity of BX(3) with a secondary contribution from the distortion energy from planar to pyramidal BX(3). The B-F bond dissociation energy of X(3)B-F(-) and the distortion energy from pyramidal to tetrahedral BX(3)(-) are of less importance in determining the relative acidities. Because the electron affinity of BX(3) is strongly influenced by the charge density in the empty p(z) lowest unoccupied molecular orbital of boron, the amount of pi-back-donation from the halogen to boron is crucial and explains why the Lewis acidity of BF(3) is significantly lower than those of BX(3) with X = Cl, Br, and I.

Heats of formation and bond energies of the H(3-n)BX(n) compounds for (X = F, Cl, Br, I, NH2, OH, and SH).

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for the borane compounds H(3-n)BX(n) for (X = F, Cl, Br, I, NH2, OH, and SH) and various radicals from coupled cluster theory (CCSD(T)) calculations with an effective core potential correlation-consistent basis set for I. In order to achieve near chemical accuracy (+/-1.5 kcal/mol), three corrections were added to the complete basis set binding energies calculated from frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, and a correction for first-order atomic spin-orbit effects. Vibrational zero point energies were computed at the MP2 level. The calculated heats of formation are in excellent agreement with the available experimental data for the closed shell molecules, but show larger differences with the reported "experimental" values for the BX2 radicals. The heats of formation of the BX2 radicals were also calculated at the G3(MP2) level of theory, and the values were in excellent agreement with the more accurate CCSD(T) values. On the basis of extensive comparisons with experiment for a wide range of compounds, our calculated values for these radicals should be good to +/-1.5 kcal/mol and thus are to be preferred over the experimental values. The accurately calculated heats of formation allow us to predict the B-X and B-H adiabatic bond dissociation energies (BDEs) to within +/-1.5 kcal/mol. The B-F BDEs in the H(3-n)BF(n) compounds and in BF (1Sigma+) are the largest BDEs in comparison to the other substituents that were investigated. The second and third largest B-X BDEs in the H(3-n)BX(n) and BX compounds are predicted for X = OH and NH2, respectively. The substituents have a minimal effect on the B-H BDEs in HBX2 and H2BX compared to the first B-H BDE of borane. The differences in adiabatic and diabatic BDEs, which are related to the reorganization energy in the product, can be estimated from singlet-triplet splittings in these molecules, and can account for the large fluctuations in adiabatic BDEs observed, specifically for the BX2 and HBX radicals, during the stepwise loss of the respective substituents.

Heats of formation and bond dissociation energies of the halosilanes, methylhalosilanes, and halomethylsilanes.

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for SiH(3)X, SiH(2)XCH(3), and SiH(3)CH(2)X with X = F, Cl, Br, and I from coupled cluster theory (CCSD(T)) calculations with effective core potential correlation-consistent basis sets for Br and I. To achieve near chemical accuracy (+/-1 kcal/mol), three 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, and a correction for first order atomic spin-orbit effects. Vibrational zero point energies were computed at the CCSD(T) level of theory and the C-H and Si-H stretches scaled to experiment. The C-H, Si-H, Si-C, C-X, and Si-X (X = F, Cl, Br, and I) bond dissociation energies (BDEs) in the halosilanes, halomethysilanes, and methylhalosilanes were predicted. Except for methyliodosilane, methyl substitution leads to an increase in Si-X BDE when compared to the Si-X BDE in the halosilanes. Except for methyliodosilane, halide substitution leads to an increase in the Si-C BDE in comparison to the Si-C BDE in methylsilane of 86.9 kcal/mol at 0 K. Unlike the methylhalosilanes, the halomethylsilanes all show a decrease in the Si-C BDE when compared to the Si-C BDE in methylsilane. The trends correlate with the electronegativity of the substituent.