Dihydrogen Bonding in Main Group Elements: A Case Study of Complexes of LiH, BH3, and AlH3 with Third-Row Hydrides

The borazine dimer: the case of a dihydrogen bond competing with a classical hydrogen bond

Dimers of borazine were studied using matrix isolation infrared spectroscopy and ab initio quantum chemical calculations. Computations were performed at the MP2 and M06-2X levels of theory using the 6-311++G(d,p) and aug-cc-pVDZ basis sets for the various homodimers. At both levels of theory, an aligned stacked structure was found to be the global minimum, which was nearly isoenergetic to a parallel displaced structure. A T-shaped structure, where the N-H of one borazine pointed towards the N of the second borazine, was found to be a local minimum. In addition to these structures, a dihydrogen bonded structure, where the hydrogen attached to the nitrogen of borazine interacted with the hydrogen attached to the boron atom of another borazine, was also indicated. Experimentally, we observed the T-shaped dimer and the dihydrogen bonded dimer. This is one of the rare examples of experimental evidence for a dihydrogen bond, in a system other than in a metal hydride. These results for the borazine dimer were clearly different from the benzene dimer where the parallel displaced structure was found to be the global minimum followed by the T-shaped structure at the MP2/aug-cc-pVDZ level of theory. AIM, EDA and NBO analyses were carried out for all the structures to explore the nature of interactions.

Hydrogen and Dihydrogen Bonds in the Reactions of Metal Hydrides

The dihydrogen bond-an interaction between a transition-metal or main-group hydride (M-H) and a protic hydrogen moiety (H-X)-is arguably the most intriguing type of hydrogen bond. It was discovered in the mid-1990s and has been intensively explored since then. Herein, we collate up-to-date experimental and computational studies of the structural, energetic, and spectroscopic parameters and natures of dihydrogen-bonded complexes of the form M-H···H-X, as such species are now known for a wide variety of hydrido compounds. Being a weak interaction, dihydrogen bonding entails the lengthening of the participating bonds as well as their polarization (repolarization) as a result of electron density redistribution. Thus, the formation of a dihydrogen bond allows for the activation of both the MH and XH bonds in one step, facilitating proton transfer and preparing these bonds for further transformations. The implications of dihydrogen bonding in different stoichiometric and catalytic reactions, such as hydrogen exchange, alcoholysis and aminolysis, hydrogen evolution, hydrogenation, and dehydrogenation, are discussed.

Time-dependent density functional theory study on excited-state dihydrogen bonding O-H···H-Ge of the dihydrogen-bonded phenol-triethylgermanium complex

Intermolecular dihydrogen bond O-H···H-Ge in the electronically excited state of the dihydrogen-bonded phenol-triethylgermanium (TEGH) complex was studied theoretically using time-dependent density functional theory. Analysis of the frontier molecular orbitals revealed a locally excited S(1) state in which only the phenol moiety is electronically excited. In the predicted infrared spectrum of the dihydrogen-bonded phenol-TEGH complex, the O-H stretching vibrational mode shifts to a lower frequency in the S(1) state in comparison with that in ground state. The Ge-H stretching vibrational mode demonstrates a relatively smaller redshift than the O-H stretching vibrational mode. Upon electronic excitation to the S(1) state, the O-H and Ge-H bonds involved in the dihydrogen bond both get lengthened, whereas the C-O bond is shortened. With an increased binding energy, the calculated H···H distance significantly decreases in the S(1) state. Thus, the intermolecular dihydrogen bond O-H···H-Ge of the dihydrogen-bonded phenol-TEGH complex becomes stronger in the electronically excited state than that in the ground state.

The B-H...H-P dihydrogen bonding in ion pair complexes [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)] (n = 0-3) and its implication in H(2) elimination and activation reactions

The B-H(delta-)...(delta+)H-P dihydrogen bonding (DHB) in ion pair complexes [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)] (n = 0-3) and its role in the combination of proton and hydride with the release of H(2) or, reversibly, the heterolytic activation of H(2) by Lewis pairs (CF(3))(3)BPH(3-n)(Me)(n) have been theoretically investigated at the MP2 and DFT levels. It is found that the B-H...H-P bonds behave similarly to those in neutral pairs and ion-molecule complexes in most respects, such as the linearity of the H...H-P moiety, the characteristics of the electron transfer and rearrangement, and the topological properties of the DHB critical point, except that in certain cases, a blue-shifting of the H-bond vibrational frequency is observed. In [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)], the proton shifting within the complexes leads to the formation of the dihydrogen complex B(CF(3))(3)(eta(2)-H(2)), which is followed by a subsequent H(2) release. The stability of B(CF(3))(3)(eta(2)-H(2)) (D(e)/D(0) = 10.8/6.0 kcal/mol) makes the proton-hydride combination proceed in a fashion similar to the protonation reactions in transition-metal hydrides rather than those in group 13 hydrides EH(4)(-) (E = B, Al, Ga). As for the H(2)-splitting reaction R(3)BPR'(3) + H(2) --> [R(3)BH(-)][HPR'(3)(+)], classical Lewis pair (CLP) (CF(3))(3)BPH(3) exhibits a high barrier and results in an unstable ion pair product [(CF(3))(3)BH(-)][HPH(3)(+)] compared with the "frustrated Lewis pair" (FLP) (C(6)F(5))(3)BP(tBu)(3). A detailed analysis of the mechanistic aspects of H(2) activation by (CF(3))(3)BPH(3) and (C(6)F(5))(3)BP(tBu)(3), supported by another CLP (CF(3))(3)BP(tBu)(3) which has a binding energy comparable to (CF(3))(3)BPH(3) but a reaction exothermicity comparable to (C(6)F(5))(3)BP(tBu)(3), allows us to suggest that the low stability of FLP (C(6)F(5))(3)BP(tBu)(3) is the determining factor for the low reaction barrier. The relative stability and other properties of the ion pair products [R(3)BH(-)][HPR'(3)(+)] have also been analyzed. Results strongly support the view proposed by Rokob et al. [ Rokob , T. A. ; Hamza , A. ; Stirling , A. ; Soos , T. ; Papai , I. Angew. Chem., Int. Ed. 2008 , 47 , 2435 ] that the frustration energy lowers the energy barrier and increases the exothermicity of the reaction.

Dihydrogen bonding: donor-acceptor bonding (AH...HX) versus the H2 molecule (A-H2-X)

Dihydrogen bonds (DHBs) play a role in, among others, crystal packing, organometallic reaction mechanisms, and potential hydrogen-storage materials. In this work we have analyzed the central H-H bond in linear H(4), LiH...HX, BH(4)(-)...HX, and AlH(4)(-)...HX complexes with various X by using the quantitative molecular orbital model contained in Kohn-Sham density functional theory at the BP86/TZ2P level of theory. First, we address the questions of if and how one can distinguish, in principle, between a H...H donor-acceptor DHB and the formation of an H(2) molecule by using the simple H(4) model system. The results of these analyses have been used to gain an understanding of the bonding in more realistic model systems (some of which have been studied experimentally), and how this differs from the bonding in H(4).

Infrared-optical double-resonance measurements on O-H...H-Ge dihydrogen-bonded phenol-triethylgermanium hydride complex in the gas phase

Spectroscopic investigation of a dihydrogen-bonded complex between phenol and triethylgermanium hydride is reported here. Laser-induced fluorescence excitation, fluorescence-detected infrared, and IR-UV hole-burning spectroscopic studies were carried out in supersonic jet to investigate the complex formation between phenol and triethylgermanium hydride. The lowering of the O-H stretching frequency of the phenol moiety in the complex with triethylgermanium hydride clearly establishes the role of phenol as hydrogen bond donor. The experimental results together with the ab-initio calculations unambiguously confirm formation of an O-H...H-Ge dihydrogen-bonded complex between phenol and triethylgermanium hydride.

Electron density topological properties are useful to assess the difference between hydrogen and dihydrogen complexes

B3LYP/6-31++G(d,p) and MP2/6-31++G(d,p) calculations for a series of hydrogen- and dihydrogen-bonded systems have been carried out in order to analyze the topology of the electron density and the energy densities at the respective energy-optimized bond critical points. Even though there are no significant differences when these properties are represented as a function of the dimerization energy, they can be separated into two well-defined sets if those properties are correlated with intermolecular distances. When analyzing the dependence of various properties with equilibrium bond lengths, the specific trends of dihydrogen bond systems consist of (a) lower electron density at the bond critical point, and (b) lower concentration/depletion of that density which can be translated in a different behavior for the Laplacian components. Furthermore, the sets of molecules form two different plots which allow for a valuable classification between hydrogen- and dihydrogen-bonded systems.

Hydrogen Bonds with π and σ Electrons as the Multicenter Proton Acceptors: High Level ab Initio Calculations

Pi-electrons of acetylene and sigma-electrons of molecular hydrogen were investigated as Lewis bases in different complexes. Hence high level ab initio calculations were performed up to the MP2/6-311++G(3df,3pd) level of approximation. It was found that species analyzed possess characteristics typical for H-bonded systems. The Bader theory was additionally applied; bond paths between proton and pi-electrons of acetylene or sigma-electrons of molecular hydrogen were detected with the corresponding bond critical points attributed to the proton-acceptor interactions. Numerous correlations between topological, geometrical and energetic parameters were also found. For example, the H...pi or H...sigma interaction is stronger for the shorter corresponding distance between the proton and the middle of C[triple bond]C or H-H bond. It is connected with the greater elongation of C[triple bond]C or H-H bonds and the greater transfer of electron charge from the Lewis base to the Lewis acid.

Proton Affinity Correlations between Hydrogen and Dihydrogen Bond Acceptors

Several series of hydrogen- and dihydrogen-bonded complexes with HCN, C2H2, HF, H2O, CH3CONH2, and CH3COOH as donors and H2O, MeOH, EtOH, MeOMe, NH3, NH2Me, NHMe2, NMe3, NEtMe2, and BH3-NMe3 as acceptors were investigated using the MP2/6-311++G(d,p) level of theory. The total lowering of the X-H stretching frequencies in the hydrogen-bonded complexes were linearly correlated with the proton affinities of the accepting bases. From comparison of hydrogen- and dihydrogen-bonded complexes, a scaling factor to estimate the exact proton affinity of a dihydrogen bond acceptor was developed. Further, the scaling factor involving linear donors (1.204) is marginally higher than that involving nonlinear donor molecules (1.162). Finally, it was found that, given identical conditions, a hydrogen bond will be about 16-20% stronger than a corresponding dihydrogen bond.

Reversible, metal-free hydrogen activation

Although reversible covalent activation of molecular hydrogen (H2) is a common reaction at transition metal centers, it has proven elusive in compounds of the lighter elements. We report that the compound (C6H2Me3)2PH(C6F4)BH(C6F5)2 (Me, methyl), which we derived through an unusual reaction involving dimesitylphosphine substitution at a para carbon of tris(pentafluorophenyl) borane, cleanly loses H2 at temperatures above 100 degrees C. Preliminary kinetic studies reveal this process to be first order. Remarkably, the dehydrogenated product (C6H2Me3)2P(C6F4)B(C6F5)2 is stable and reacts with 1 atmosphere of H2 at 25 degrees C to reform the starting complex. Deuteration studies were also carried out to probe the mechanism.

Complete infrared spectroscopic characterization of phenol-borane-trimethylamine dihydrogen-bonded complex in the gas phase

In the paper we report the first observation of the vibrational spectrum in the B-H stretching region in the gas phase for a dihydrogen bonded complex. The appearance of three transitions for the B-H stretching modes of a (di)hydrogen-bonded complex involving borane-trimethylamine indicates the lowering of the symmetry on the BH3 group upon interaction with phenol. Further, the shift in the O-H stretching frequency indicates that phenol is hydrogen bonded to borane-trimethylamine. The two sets of the present data establish, unequivocally, the formation of O-H...H-B dihydrogen-bonded complex between phenol and borane-trimethylamine.

Electron bridging dihydrogen bond in the imidazole-contained anion derivatives

The large contact distance of electron bridging dihydrogen bond (EBDB), which is over 2.4 A, is the most prominent characteristic for the imidazole-contained anion derivatives. The elongation of N-H bond and the shortening of H...H distance can be observed upon hydration and hydrogenation. Transformation from EBDB to dissociative H2 is convenient upon sequential hydrogenation. The H...H distance decreases with the enhancement of the electronegativity of the heavy atom which contacts directly with one of these two hydrogen atoms. NMR shielding of the bonding N varies significantly upon hydration and hydrogenation. The spin-spin coupling constants, 1J(H-H), is dominated predominantly by the paramagnetic spin-orbit and diamagnetic spin-orbit contributions instead of the Fermi-contact term. Enhancement of electronegativity of the heavy atom leads to the increase of 1J(H-H) coupling constants. The stabilization is enhanced upon hydration predominantly for the formation of O-H...N H bond, while it is reversed upon hydrogenation for the cleavage of big pi bond, Pi5(6). Enhancement of the stability is demonstrated by the increase of stabilization energy and vertical electron detachment energy with the electronegativity of the heavy atom. The dominant contributions for the formation of such electron bridging dihydrogen bond are the high polarity of each fragment, large electron density between two fragments, and strong bonding interaction of the bridging electron with H(N) atoms. The H...H interaction can be formed by X-Hdelta+ and Hdelta- -Y polar molecules in Hdelta+...Hdelta- and Hdelta+...e...Hdelta+ of two forms.

Thermodynamic Properties of Molecular Borane Phosphines, Alane Amines, and Phosphine Alanes and the [BH4-][PH4+], [AlH4-][NH4+], and [AlH4-][PH4+] 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).

Theoretical Study on the Small Clusters of LiH, NaH, BeH2, and MgH2

High-level ab initio molecular orbital theory is used to calculate the geometries, vibrational frequencies, atomic charges, and binding energies of the small clusters (LiH)(n), (NaH)(n), (BeH(2))(n), and (MgH(2))(n) (n = 1-4). For (LiH)(n) and (NaH)(n), there are planar cyclic structures when n = 2, 3. We have found the cubic structure T(d) in addition to the planar cyclic D(4)(h) when n = 4. The D(4)(h) is less stable than the T(d) geometry. For (BeH(2))(n) and (MgH(2))(n), when n = 3, there are three kinds of structures: chain C(2)(v), planar cyclic D(3)(h), and hat-like C(2)(v). The C(2)(v) geometry is more stable than the others. When n = 4, there are four kinds of structures: chain D(2)(h), cubic T(d), string-like C(2), and cubic transformation C(1). The most stable compounds in the families of (LiH)(n), (NaH)(n), (BeH(2))(n), and (MgH(2))(n) are cubic T(d), cubic T(d), chain D(2)(h), and string-like C(2) geometries, respectively, when n = 4. Calculated binding energies range from -24 to -37 kcal/mol for (LiH)(n) and --19 to -30 kcal/mol for (NaH)(n), (BeH(2))(n), and (MgH(2))(n). The hydrogen atoms in hydride clusters always have negative charges. The atomic charges of planar cyclic structures are weaker than those of cubic structures, and there is a tendency of reducing along with the increase of the cluster size. The vibrational frequencies of planar cyclic structures have consistent tendency, too. It indicates that the bond distance increases with the ionic character of the bond.

How Short Can the H···H Intermolecular Contact Be? New Findings that Reveal the Covalent Nature of Extremely Strong Interactions

Ab initio calculations at the MP2/6-311++G(d,p) and MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ levels have been performed for the following complexes: H2OH+...HBeH, H2OH+...HBeBeH, H2OH+...HBeF, HClOH+...HBeH, Cl2OH+...HBeH, and Cl2OH+...HBeF. For all dimers considered, extremely short H...H intermolecular contacts (1.0-1.3 A) were obtained. These are the shortest intermolecular distances which have ever been reported, with binding energies within the range of 13.7-24.3 kcal/mol (MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ level). The interaction energies of the complexes analyzed were also extrapolated to the complete basis set (CBS) limit. To explain the nature of such strong interactions, the Bader theory was applied, and the characteristics of the bond critical points (BCPs) were analyzed. It was pointed out that for the major part of the H...H contacts considered here the Laplacian of the electron density at H...H BCP is negative indicating the partly covalent nature of such a connection. The term "covalent character of the hydrogen bond" used sometimes in recent studies is discussed. An analysis of the interaction energy components for dihydrogen bonded systems considered indicates that in contrast to conventional hydrogen bonded systems the attractive electrostatic term is outweighed by the repulsive exchange energy term and that the higher order delocalization energy term is the most important attractive term.

Proton Affinities of Borane−Amines: Consequences on Dihydrogen Bonding

The calculated proton affinities of four borane-amines using Gaussian-2 theory have been found to be comparable to conventional bases such as water, methanol, and ammonia. On the other hand the structure of protonated borane-ammonia, [HBH(3)-NH(3)](+), is found to be drastically different from that of protonated ammonia, [HNH(3)](+), and can appropriately be described as a eta(2)-H(2) complex with [BH(2)-NH(3)](+) molecular cation. Further, the proton affinities of borane-amines are related to the ease of H(2) elimination.

C-H···H-C interactions in organoammonium tetraphenylborates: another look at dihydrogen bonds

The crystal structures of the tetraphenylborates of the dabcoH+, guanidinium (MeCN solvate), and biguanidinium cations are shown to contain a variety of C-H···H-C dihydrogen (DB) bonds of nominally zero polarity, as well as a variety of N-H···N, C-H···N, N-H···Ph, and C-H···Ph hydrogen (HB) bonds. These intermolecular bonds have been characterized topologically after multipole refinement of the structures. The coexistence of the DBs and HBs in each of the structures makes it possible to establish their relative strength hierarchy. It also illustrates the importance of the DBs in satisfying the tendency of these structures to maximize the total intermolecular bonding engagement. To compare the above DBs with other DBs, the results of an extensive set of MP2/6-31G(d,p) calculations (supplied by I. Alkorta) were analyzed for reference correlations between the bond-critical parameters. Thus, for an X-H···H-Y bond, the difference Δε(H)m between the Mulliken charges on the H atoms in the uncomplexed X-H and H-Y components correlates quite well with the X-H···H-Y parameters and can be used for predicting the topological strength of an X-H···H-Y bond. The use of the difference Δε(H)c in the bond does not appear to change the correlation significantly; closer correlations are observed when the amount of charge transferred on formation of the H···H bond is used instead of Δε(H)m or Δε(H)c. Bonding interactions are obtained even between like or symmetry-related H atoms as a consequence of induced-dipole interactions, which accounts for the existence of the above intermolecular C-H···H-C bonds with d(H···H) = 2.18–2.57 Å, electron density at the bond-critical point of ~0.05–0.08 e/Å3, and a rough estimate of the H···H binding energy of ~1-5 kcal/mol. Examination of the bond-critical parameters of X-H···H-Y bonds also suggests a criterion of stability of these bonds with respect to the transition from non-shared (closed-shell) X-H···H-Y interaction to covalent (shared-shell) X···H-H···Y interaction. This transition appears to be discontinuous.Key words: bond-critical parameters, bond topology, dihydrogen bonds, hydrogen bonds, organoammonium tetraphenylborates.