Daniel J. Grant

Coordination of aminoborane, NH2BH2, dictates selectivity and extent of H2 release in metal-catalysed ammonia borane dehydrogenation

In situ(11)B NMR monitoring, computational modeling, and external trapping studies show that selectivity and extent of H(2) release in metal-catalysed dehydrogenation of ammonia borane, NH(3)BH(3), are determined by coordination of reactive aminoborane, NH(2)BH(2), to the metal center.

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.

Thermochemistry and Electronic Structure of Small Boron Clusters (Bn, 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.

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.

Thermochemistry and electronic structure of small boron and boron oxide clusters and their anions.

Thermochemical properties of a set of small boron (B(n)) and boron oxide (B(n)O(m)) clusters, with n = 1-4 and m = 0-3, their anions, and the B(4)(2-) dianion, were calculated by using coupled-cluster theory CCSD(T) calculations with the aug-cc-pVnZ (n = D, T, Q, 5) basis sets extrapolated to the complete basis set limit with additional corrections. Enthalpies of formation, bond dissociation energies, singlet-triplet or doublet-quartet separation gaps, adiabatic electron affinities (EA), and both vertical electron attachment and detachment energies were evaluated. The predicted heats of formation show agreement close to the error bars of the literature results for boron oxides with the largest error for OBO. Our calculated adiabatic EAs are in good agreement with recent experiments: B (calc, 0.26 eV; exptl, 0.28 eV), B(2) (1.95, 1.80), B(3) (2.88, 2.820 +/- 0.020), B(4) (1.68, 1.60 +/- 0.10), BO (2.50, 2.51), BO(2) (4.48, 4.51), BOB (0.07), B(2)O(2) (0.37), B(3)O (2.05), B(3)O(2) (2.94, 2.94), B(4)O (2.58), and B(4)O(2) (3.14, 3.160 +/- 0.015). The BO bond is strong, so this moiety is maintained in most of the clusters. Thermochemical parameters of clusters are not linearly additive with respect to the number of B atoms. The EA tends to be larger in the dioxides. The growth mechanism of small boron oxides should be determined by a number of factors: (i) formation of BO bonds, (ii) when possible, formation of a cyclic B(3) or B(4), and (iii) combination of a boron cycle and a BO bond. When these factors compete, the strength of the BO bonds tends to compensate the destabilization arising from a loss of binding in the cyclic boron clusters, in such a way that a linear boron oxide prevails. When the B(2) moiety is present in these linear clusters, the oxide derivatives prefer a high spin state.

Heats of Formation of the H1,2OmSn (m, n = 0−3) Molecules from Electronic Structure Calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted from high level ab initio electronic structure calculations using the coupled cluster CCSD(T) method with augmented correlation-consistent basis sets extrapolated to the complete basis set (CBS) limit for the H(1,2)O(m)S(n) (m, n = 0-3) compounds, as well as various radicals involved in different bond breaking processes. To achieve near chemical accuracy (+/-1.0 kcal/mol), additional corrections were added to the CBS binding energies based on the frozen core CCSD(T) energies including corrections for core-valence, scalar relativistic, and first-order atomic spin-orbit effects. Geometries were optimized up through the CCSD(T)/aV(T+d)Z level. Vibrational zero point energies were computed at the MP2/aV(T+d)Z level. The calculated heats of formation are in excellent agreement with the available experimental data and allow the prediction of adiabatic bond dissociation energies (BDEs) to within +/-1.0 kcal/mol. The decomposition mechanisms were largely determined by a preference to maintain a strong S=O bond in the dissociated products as opposed to O=O and S=S bonds, exactly matching the ordering of the BDEs in the diatomics. For the H(2)X(2) and H(2)X(3) systems, as well as the HX(3) radicals, the energetically favorable decomposition pathway leads to the formation of XH radicals and breaking the X-X bond as opposed to breaking the X-H bond. For the HX(2) radicals, however, the more thermodynamically favorable pathway leads to a breaking of the H-X bond and forming X(2) molecules.

Lewis Acidities and Hydride, Fluoride, and X- Affinities of the BH3-nXn 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.

Bond Dissociation Energies in Second-Row Compounds

Heats of formation at 0 and 298 K are predicted for PF3, PF5, PF3O, SF2, SF4, SF6, SF2O, SF2O2, and SF4O as well as a number of radicals derived from these stable compounds on the basis of coupled cluster theory [CCSD(T)] calculations extrapolated to the complete basis set limit. In order to achieve near chemical accuracy (+/-1 kcal/mol), additional 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 vibrational zero-point energies. The calculated values substantially reduce the error limits for these species. A detailed comparison of adiabatic and diabatic bond dissociation energies (BDEs) is made and used to explain trends in the BDEs. Because the adiabatic BDEs of polyatomic molecules represent not only the energy required for breaking a specific bond but also contain any reorganization energies of the bonds in the resulting products, these BDEs can be quite different for each step in the stepwise loss of ligands in binary compounds. For example, the adiabatic BDE for the removal of one fluorine ligand from the very stable closed-shell SF6 molecule to give the unstable SF5 radical is 2.8 times the BDE needed for the removal of one fluorine ligand from the unstable SF5 radical to give the stable closed-shell SF4 molecule. Similarly, the BDE for the removal of one fluorine ligand from the stable closed-shell PF3O molecule to give the unstable PF2O radical is higher than the BDE needed to remove the oxygen atom to give the stable closed-shell PF3 molecule. The same principles govern the BDEs of the phosphorus fluorides and the sulfur oxofluorides. In polyatomic molecules, care must be exercised not to equate BDEs with the bond strengths of given bonds. The measurement of the bond strength or stiffness of a given bond represented by its force constant involves only a small displacement of the atoms near equilibrium and, therefore, does not involve any reorganization energies, i.e., it may be more appropriate to correlate with the diabatic product states.

Thermochemistry for the Dehydrogenation of Methyl-Substituted Ammonia Borane Compounds

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for (CH3)H2N-BH3, (CH3)HN=BH2, (BH3)HN=CH2, (CH3)H2B-NH3, (CH3)HB=NH2, and (NH3)HB=CH2, as well as various molecules involved in the different bond-breaking processes, from coupled cluster theory (CCSD(T)) calculations. In order to achieve near-chemical accuracy (+/-1 kcal/mol), three corrections were added to the complete basis set binding energies based on frozen core CCSD(T) energies, corrections for core-valence, scalar relativistic, and first-order atomic spin-orbit effects. Scaled vibrational zero-point energies were computed with the MP2 method. The heats of formation were predicted for the respective dimethyl- and trimethyl-substituted ammonia boranes, their dehydrogenated derivatives, and the various molecules involved in the different bond breaking processes, based on isodesmic reaction schemes calculated at the G3(MP2) level. Thermodynamics for dehydrogenation pathways in the monomethyl-substituted molecules were predicted. Dehydrogenation across the B-N bond is more favorable as opposed to dehydrogenation across the B-C and N-C bonds. Methylation at N reduces the exothermocity of the dehydrogenation reaction and makes the reaction more thermoneutral, while methylation at B moves it away from thermoneutral. Various mixtures of CH3NH2BH3 and NH3BH3 were made, and their melting points were measured. The lowest melting mixture contained approximately 35% NH3BH3 by weight and melted at 35-37 degrees C.

Structure and heats of formation of iodine fluorides and the respective closed-shell ions from CCSD(T) electronic structure calculations and reliable prediction of the steric activity of the free-valence electron pair in ClF6-, BrF6-, and IF6-.

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for IF, IF2-, IF2+, IF3, IF4-, IF4+, IF5, IF6-, IF6+, IF7, IF8-, BrF6-, and ClF6- from coupled cluster theory [CCSD(T)] calculations with effective-core potential correlation-consistent basis sets for I. In order 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 coupled-cluster level of theory except for IF6-, IF7, and IF8-. The calculated heats of formation for the neutral and ionic IFn fluorides were used to predict fluoride affinities. It is shown that high-level calculations are required to predict correctly the steric activity of the free-valence electron pair on the central atoms in IF6- (C3v), BrF6- (Oh), and ClF6- (Oh ). The vibrational spectrum of IF8- was reanalyzed, and complete mode descriptions for square-antiprismatic XF8 species of D4d symmetry are given.