A comparison of hydrogen release kinetics from 5- and 6-membered 1,2-BN-cycloalkanes.

The reaction order and Arrhenius activation parameters for spontaneous hydrogen release from cyclic amine boranes, i.e., BN-cycloalkanes, were determined for 1,2-BN-cyclohexane (1) and 3-methyl-1,2-BN-cyclopentane (2) in tetraglyme. Computational analysis identified a mechanism involving catalytic substrate activation by a ring-opened form of 1 or 2 as being consistent with experimental observations. This journal is © The Royal Society of Chemistry.

Amine boranes have long been targeted as a promising class of potential hydrogen storage materials.[1] The simplest amine borane, ammonia borane (H3NBH3), has attracted significant attention by virtue of its particularly high gravimetric hydrogen density (19.6 wt% H2).[2] Thermal decomposition of ammonia borane to release H2, however, frequently generates mixtures of oligomeric and polymeric products, which in turn can complicate efforts to regenerate the fully saturated starting material for reuse.[1] Extensive studies of variously substituted amine boranes have thus been conducted in the interest of identifying a system that releases H2 more selectively.[3,4]

Our approach in this regard has been to incorporate the amine borane H2NBH2 unit into carbocyclic structures to form saturated carbon–boron–nitrogen (CBN) heterocycles.[5] We have previously discovered that two of these compounds, 1,2-BN-cyclohexane (1)[6] and its constitutional isomer 3-methyl-1,2-BN-cyclopentane (2)[7] undergo full, thermally-induced BN-dehydrogenation to afford only the trimeric products 3 and 4 through the intermediacy of monomeric BN-“cycloalkene” species 1–H2 (ref. 8) and 2–H2, respectively (Scheme 1).

Scheme 1

Selective thermal decomposition of 1 and 2 to trimers 3 and 4, respectively.

While 1 and 2 appeared to exhibit the same general dehydrogenation selectivity, in subsequent experiments using neat material, we observed an intriguing difference in the thermal stability of these compounds. Specifically, thermo-gravimetric analysis-mass spectrometry (TGA-MS) revealed significant loss of H2 from 2 initiated at ∼50 °C,[9] but an analogous measure of decomposition of 1 occurred only upon heating to ∼70 °C (Fig. S2; see also Fig. S3† for the temperature-programmed desorption-mass spectrometry (TPD-MS) of 1). Thus, of the two materials, only 1 would meet the minimum requirement for operational stability set by the US Department of Energy's Hydrogen and Fuel Cells Program for on-board vehicular applications.[10]

A mechanistic investigation of the origin of this dichotomy of dehydrogenation reactivity would aid in the design of better amine borane-based hydrogen storage materials and add to our fundamental understanding of the reactivity of amine boranes. Herein, we provide a solution-phase kinetic analysis of initial H2 release from 1 and 2 using ReactIR. The kinetic data establish a second-order decomposition pathway that is in agreement with a computationally-derived mechanistic model. We also provide evidence that the ring strain associated with the 5-membered heterocycle 2 is ultimately responsible for its faster decomposition rates.[11,12]

The reaction order for the first step of dehydrogenation of 1 and 2 was determined using the initial rates method. Both 1 and 2 exhibit a characteristic IR frequency at 1600 cm−1 (attributed to an N–H bending mode); the disappearance of starting material can thus be readily monitored in situ by ReactIR.[13] The initial concentrations ([ ]0) of either 1 or 2 in a tetraglyme solution at 140 °C were varied between 0.560 M and 1.283 M, and initial rates (ri) were estimated based on linear regression of the respective portions of the substrate disappearance trends representing up to 20% conversion. As shown in Fig. 1, a linear fit of ln(ri) vs. ln([ ]0) yielded a slope of 2 for both 1 and 2, indicating the rate of initial H2 loss from each follows a second-order concentration dependence.[14]

Fig. 1

Abbreviated substrate disappearance trends measured at 140 °C for various initial concentrations of 1 (A) and 2 (C), and reaction order determination by the initial rates method for 1 (B) and 2 (D).

A subsequent Arrhenius analysis was enabled by maintaining a constant initial concentration of 0.742 M, and varying the reaction temperature from 120–160 °C for 1 and from 100–160 °C for 2. As shown in Table 1, the activation energy (Ea) for 2 is lower than that of 1 by 5.0 kcal mol−1. However, the pre-exponential factor (A) of 2 is an order of magnitude smaller than that of 1. Of these two opposing trends, the difference in Ea ultimately dominates within the temperature range studied, and so the reaction rate constants (k) of 2 are more than an order of magnitude greater than those of 1 (Fig. 2). (Extrapolation to lower temperatures, e.g., 50 °C, reveals an even wider gap of approximately three orders of magnitude.) It thus appears that enthalpic factors, as represented by Ea values,[15] are primarily responsible for the faster kinetics of H2 release from 2 relative to those of 1.

Experimentally determined activation energies (Ea) and pre-exponential factors (A) for 1 and 2 based on Arrhenius analyses

Reaction	E a (kcal mol−1)	A (M−1 s−1)
	23.8	2.25 × 108
	18.8	2.02 × 107

Fig. 2

Comparison of rate constants (k) for 1 (red) and 2 (blue) from 120 °C to 160 °C. Solid lines represent predicted values based on the Arrhenius activation terms listed in Table 1.

To further refine our understanding of the mechanism of the initial H2 release from either 1 or 2, we turned to computational modeling. The experimental evidence (vide supra) for a bimolecular process involving two molecules of 1 or 2 in the rate-determining step led us to investigate four distinct possible scenarios (Scheme 2): (A) simultaneous loss of one H2 equivalent each directly from two substrate molecules interacting in a “head-to-tail” fashion;[16] (B) formation of an intermediate prior to H2 release analogous to the diammoniate of diborane (DADB) species observed in thermal ammonia borane decomposition;[14] or generation of linear isomers 1′ or 2′ through reversible, heterolytic B–N bond cleavage and subsequent catalysis of H2 formation by interaction of these isomers' free BH2 (C) or NH2 groups (D) with another still cyclized molecule of either 1 or 2.[17] Ultimately, the activation barriers associated with pathways (A), (B), and (D) (Fig. S6–S9†) were calculated to be significantly higher for both 1 and 2 than those predicted for pathway (C), which is shown in detail in Fig. 3 and described in terms of 1 below.

Scheme 2

Proposed bimolecular decomposition scenarios.

Fig. 3

Potential energy surfaces for pathway (C): bimolecular H2 release from 1 (top) and 2 (bottom) involving catalysis by ring-opened intermediates. Gas-phase enthalpy (black) and free energy [green] values (kcal mol−1) were calculated using DFT (298 K) at the M06/DZVP2 level of theory.

Starting from two ring-closed molecules, B–N bond cleavage first generates one unit of 1′. A bridging hydride interaction between the BH2 groups of 1 and 1′ results in formation of complex [1–1′], which is lower in energy than the separate mixed species, but still less stable than the ring-closed starting materials. The lowest energy transition state (TS-1) for H2 release from 1 involves intramolecular transfer of the bridging hydride from 1 completely to 1′, while a proton from the NH2 group of 1 simultaneously combines with another hydride of the 1′ BH2 unit to form free H2.

The above gas-phase model predicts a lower reaction barrier for H2 release from 2 than from 1 (ΔG(2‡) = 38.7 kcal mol−1vs. ΔG(1‡) = 43.1 kcal mol−1), consistent with the trend observed experimentally. (This relative trend was also observed in calculations with an implicit tetraglyme solvent model (Fig. S10†)). Interestingly, the calculated ΔG298 K value for B–N bond dissociation in 2 is also lower than for 1 (ΔΔG298 K = 2.8 kcal mol−1), such that the equilibrium constants for this preliminary step differ by two orders of magnitude at 25 °C (Keq(2) = 8.15 × 10−16 and Keq(1) = 7.08 × 10−18). While our model does not predict B–N bond cleavage itself to be rate-limiting, a higher equilibrium concentration of 2′ relative to 1′ would nonetheless influence the observed kinetics of thermal decomposition in favor of faster apparent rates for 2.[18]

We propose both the greater facileness of initial B–N bond dissociation and the overall more rapid dehydrogenation of 2 occurs primarily as a result of greater molecular strain energy in 2 as compared to that in 1. Using the homodesmotic reaction scheme devised by Gilbert,[19] we calculated the strain energy of 2 to be 3.1 kcal mol−1 greater than that of 1. Furthermore, experimental evidence in support of this trend arises from a comparison of the length of the B–N bond in five- and six-membered BN-cycloalkanes as determined by single crystal X-ray analysis (see the ESI† for acquisition parameters and detailed structural data). Relative elongation of this bond was presumed to serve as an indicator of greater ring strain, and indeed, the average B–N bond length of a series of BN-cyclopentanes (1.630 Å)[20] was found to be longer than that measured for 1 (1.614(1) Å).[21,22] This additional strain destabilizes 2 relative to 1, such that the enthalpic contributions required to reach TS2 are correspondingly diminished compared to those needed to reach TS1 starting from 1. These relationships are reflected in both experimental and computational results, which respectively yield lower Ea and ΔH‡ values for 2 compared to those for 1.

To probe the intermediacy of ring-opened species such as 1′, we prepared a close analogue of 1 that contains a stereochemical label: trans-1,2-dimethyl-1,2-BN-cyclohexane (trans-5).[23] Sequential hydride/proton addition to cyclic aminoborane 5–H2 (ref. 24) furnishes a mixture of trans- and cis-5 (Scheme 3). The major trans diastereomer could be isolated by recrystallization, and its structure was unambiguously confirmed by single crystal X-ray diffraction analysis.

Scheme 3

Synthesis and structure of trans-5.

Mild heating of pure trans-5 leads to a partial isomerization back to the cis isomer as evidenced by 11B and 1H NMR (Fig. S1†). Based on previous studies of B–N bond cleavage in related cyclic systems,[25] the observed formation of cis-5 is consistent with a mechanism that involves B–N bond dissociation, B–(C3) bond rotation (or nitrogen inversion and N–(C6) rotation), and finally B–N bond re-formation. The studies with 5 thus provide indirect evidence for ring-opened species 1′ and 2′ as viable intermediates in the thermal decomposition of 1 and 2.

In summary, we have experimentally measured the kinetics of the release of the first H2 equivalent from 1,2-BN-cyclohexane (1) and 3-methyl-1,2-BN-cyclopentane (2) using ReactIR. A second-order concentration dependence was determined for both 1 and 2. Arrhenius analysis revealed a lower reaction barrier for 2 due to a smaller requisite activation energy (Ea). These trends were replicated in a computational model of a bimolecular dehydrogenation mechanism involving substrate activation catalyzed by a ring-opened form of 1 or 2. This mechanistic study sheds light on the origin of the differing thermal stability exhibited by two isomeric cyclic amine boranes and suggests that ring-strain as well as the strength of the B–N bond in cyclic CBN compounds need to be considered in the next generation of materials to provide sufficient thermal stability.

Conflicts of interest

There are no conflicts to declare.