Dihydrogen intermolecular contacts in group 13 compounds : H ⋯ H or E ⋯ H ( E = B , Al , Ga ) interactions ? †

A systematic theoretical analysis of homopolar dihydrogen interactions in group 13 compounds is presented here. Ab initio calculations and structural analysis allow us to demonstrate that interactions involving B–H⋯H–B contacts are comparable in strength to the previously studied C–H⋯H–C ones, yet attractive and important for the stabilization of dimers of large molecules. We have also shown that a polyhedral skeleton enhances the B–H⋯H–B interaction strength with respect to non-polyhedral compounds, and it has also been proven that Al–H⋯H–Al and Ga–H⋯H–Ga interactions can be attractive in some cases. If H⋯E (B, Al and Ga) short contacts are present, the interaction is significantly strengthened, especially for Al and Ga. In general, H⋯H interactions combined with associated H⋯E (B, Al and Ga) short contacts are responsible for the stability of a large number of dimers of group 13 compounds and may play an important role in the packing of their crystal structures.


Introduction
2][3][4][5] Dihydrogen contacts are a type of intermolecular interaction that involves two hydrogen atoms with the same (X-H⋯H-X) or different polarities (X-H⋯H-Y). 6sually, the term dihydrogen bond is reserved for the second case whereas the first one is referred to as dihydrogen interaction. 7While X-H δ+ ⋯ δ− H-Y dihydrogen bonds have an electrostatic contribution due to the opposite charge of the two hydrogen atoms, the homopolar X-H δ+ ⋯ δ+ H-X dihydrogen interactions are generally classified as weak and governed by dispersion.The dihydrogen bonding of heteropolar units such as C-H⋯H-B 8 or N-H⋯H-B 9,10 has been widely studied and many experimental 11,12 and theoretical [13][14][15][16] works can be found in the literature.However, homopolar contacts have attracted less attention and only in recent years has there been an increase in the number of publications on the subject.For example, the surprising strength of some C-H⋯H-C contacts was for the first time evaluated via a theoretical analysis, 17,18 and confirmed experimentally by the stabilization of alkanes with very long C-C bonds. 19On the other hand, in group 13 compounds, B-H⋯H-B interactions have been observed in boron hydrides and complex metal hydrides 20 and lyotropic lamellar phases of carborane-cage amphiphiles, 21 and they seem to play an important role in the design of materials for hydrogen storage. 22Descending down the group 13, Al-H⋯H-Al intramolecular contacts have been found in several polymorphs of the crystal structure of alane, 23 whereas short Ga-H⋯H-Ga contacts can be found only in a handful of structures and, to our knowledge, there is no theoretical work in the literature devoted to them.However, short heteropolar N-H⋯H-Ga contacts have been detected, for example, in the crystal structure of methylamine-gallane. 24 Homopolar dihydrogen contacts have not been observed in compounds of the two heaviest elements of group 13, In and Tl.
From a theoretical point of view, dihydrogen bonding has been approached in several ways.The quantum theory of atoms in molecules (AIM) has been used by McGrady and Wolstenholme to prove the existence of a bond path typical of attractive interactions between two hydrogen atoms of the same polarity. 20The electron localization function (ELF) can be also employed to find and characterize dihydrogen bonds. 25odern Valence Bond (VB) 18 and, especially, Molecular Orbital (MO)-based methods have been widely used to investigate the nature of C-H⋯H-C interactions (see section 2.1 in ref. 1 for a comprehensive account of theoretical methods that take into account dispersion forces).
In this work, we undertake a combined structural and systematic ab initio analysis of intermolecular E-H⋯H-E contacts in group 13 compounds, focusing on those with a homopolar character.E⋯H (B, Al and Ga) short contacts will also be evaluated since they are often associated with dihydrogen contacts.Since Hartree-Fock theory completely neglects dispersion, we have employed second-order Møller-Plesset (MP2) perturbation theory combined with a large basis set to capture the non-covalent nature of the interaction.This theoretical approach has been extensively tested and compared to high level calculations such as CCSD(T)/CBS, giving very good results for the study of H⋯H contacts in alkanes 17 and other families of compounds. 26We have studied small boranes, as well as larger arachno, nido and closo polyhedral boranes, and also Al and Ga hydrides.The heaviest elements of group 13, In and Tl, have been excluded in the present study, because these atoms are affected by relativistic effects 27 and, thus, their hydrides behave differently from those of the lighter elements of the group (B, Al and Ga).

Boron hydrides
We have searched the Cambridge Structural Database (CSD) 28 for E-H⋯H-E (E = B, Al, Ga, In and Tl) contacts shorter than twice the van der Waals radius of hydrogen 29 plus 25%, i.e. 3 Å, in order to have a global picture of the occurrence of these interactions.B-H⋯H-B contacts were the most numerous ones with 3542 hits (Fig. 1a).We have also analyzed the B-H⋯H angles associated with those contacts.Since angles smaller than 80°are unlikely due to steric hindrance, values for B-H⋯H angles range from 80°to 180°and show a clear maximum at around 115°(Fig. 1b).This is similar to what was found for R 3 C-H⋯H-CR 3 contacts, where the C-H⋯H angle presents a peak at 125°. 18s for the quality of the analysed data, it is important to assess how the positions of the H atoms have been determined in the corresponding crystal structures.We have only found three examples of neutron diffraction structures, containing 37 short B-H⋯H-B contacts (CSD refcodes FUYYIH02, 30 KOMZAN, 31 SOHZIZ 32 ), which exclude the possibility of a meaningful statistical analysis.However, X-ray structures can also be useful, especially when the H positions have been refined.On the other hand, structures where the H atoms have been artificially placed at fixed distances (1.1 and 1.12 Å in the case of B-H) can introduce relatively large errors in the averaged H⋯H distances.In Fig. 1a, we compare the H⋯H distances obtained with and without the inclusion of this set of structures with fixed hydrogens, obtaining similar results.In both cases, the van der Waals peak can be adjusted to a Gaussian function centred at 2.91 Å ± 0.02.This can be due to the fact that the geometrical parameters associated with the interaction topologies, with B-H⋯H angles around 115°, compensate in part the error introduced by the fixed bond lengths.Only in the case of linear topologies (B-H⋯H angles close to 180°) would the error in the H⋯H distance be added.
An 84% of all the B-H⋯H-B short contacts was found to correspond to the crystal structures of polyhedral boranes.This is evidenced by the data analysis shown in Fig. 1a, where the B-H⋯H-B contacts shorter than 5 Å are presented.The blue bars represent the frequency of all B-H⋯H-B contacts, while the grey bars represent the subset of contacts in which both B atoms involved in the interaction are hexacoordinated, which is the most common coordination number for the B atom in polyhedral boranes.Within the remaining 16%, borane adducts with Lewis bases like amines and phosphines are abundant.In the crystal structures of these compounds, B-H⋯H-B interactions coexist with heteropolar dihydrogen interactions such as N-H⋯H-B, which are usually shorter and contribute to a considerable extent to the stabilization of the crystal structure.For example, a methylamine-borane crystal shows N-H⋯H-B contacts as short as 2.13 Å, whereas the B-H⋯H-B ones are at a distance of 2.97 Å. 33 Another interesting case is the crystal structure of 1-aza-closo-dodecaborane, in which two short N-H⋯H-B (1.99-2.57Å) and two B-H⋯H-B (2.70-2.97Å) intermolecular contacts can be found. 34he situation changes for polyhedral boranes, where B-H⋯H-B contacts usually predominate over other possible heteropolar dihydrogen interactions.The clusters of closo type, especially the derivatives of the neutral species dicarbadodecaborane (C 2 H 12 B 10 ), are able to establish mainly B-H⋯H-B interactions only between terminal hydrogen atoms (H t -H t ).For example, a carborane dimer found in the crystal structure of a mixed ytterbium(II) complexdicarbadodecaborane (CSD refcode QEQFUN) 35 presents three B-H⋯H-B short contacts (2.60-2.65 Å) with a 1 : 3 interaction topology (Fig. 2).However, this is not strictly a homopolar dihydrogen bond since the B-H units involved in the interactions in each monomer have different polarities due to the relative positions of the two carbon atoms.Arachno and nido clusters have, in addition, bridging hydrogen atoms that can establish short contacts with terminal hydrogens (H t ⋯H b ) or other bridging ones (H b ⋯H b ).
It must be noted that the three H⋯H contacts shown in Fig. 2 correspond to distances in the range 2.60-2.65 Å, slightly longer than twice the van der Waals radius of hydrogen (2.40 Å), 36 but still within the van der Waals peak shown in Fig. 1a.On the other hand, the corresponding H⋯B distances (3.00-3.04Å) are slightly shorter than the sum of van der Waals radii (3.11 Å).This means that the dispersion forces holding the two molecules together involve not only H⋯H, but also B⋯H interactions, that may be as strong or even more important, as was already commented upon in our previous study of C-H⋯H-C contacts. 17In what follows, thus, it must be understood that reference to E-H⋯H-E interactions will refer to the combined effect of H⋯H and H⋯E interactions.
The angular distribution of the short B-H⋯H-B contacts in the crystal structures of closo boranes (i.e., those contained within the van der Waals peak of Fig. 1a, with H⋯H ≤ 3.5 Å) is shown in Fig. 1b, where one can appreciate a narrow peak centered at around 110°-120°, and a broader one centered at around 150°.In our example structure, (Fig. 2) the two B-H⋯H angles associated with the 1 : 3 interaction topology, α and β, present values in the ranges 129-167°and 100-103°, respectively, and a structural database search for those types of contacts involving six-coordinated boron atoms (1) shows the presence of two distinct peaks corresponding to the two sets of B-H⋯H angles, α and β, whose distributions are centered at around 130 and 99°, respectively (Fig. 1b).Other contact topologies, such as 2 and 3 present a distribution that is more similar to the overall distribution, with all topologically independent sets of B-H⋯H angles (γ, δ, τ and ω in 2-4) presenting sharp maxima at around 110°followed by a shoulder that extends all the way to 180°(see Fig. S1 in the ESI †).The interaction topology 1 : 1 that involves one BH group from each molecule ( 4) is less common and only a few examples can be found among the family of polyboranes.
In the nido structures, the situation is somewhat more convoluted because of the lower symmetry.Hence, these boranes can establish contacts involving terminal and/or bridging hydrogen atoms (H t ⋯H t , H t ⋯H b and H b ⋯H b ).Moreover, there are more relative orientations of the two interacting molecules, which we broadly identify as involving B-H groups from the bottom (b), the sides (s) or the face (f ) of the nest (5).We will therefore classify somewhat crudely, the interaction topologies in dimers of nido boranes as b-f, f-f, b-b, s-s… It is worth pointing out that we have observed a B-H⋯H-B intermolecular distance of 2.28 Å between dodecaborate anions in the lithium salt, Li 2 (B 12 H 12 ), 37 which is shorter than those found between neutral molecules.This fact could be in part attributed to the Coulombic attraction caused on two neighboring anions by an intervening cation, since the contacts become longer as the size of the intervening cation increases, with a practically linear dependence on the ionic radius of the alkaline cation (2.42, 2.51, 2.64 and 2.73 Å for the Na, 38 K, 39 Rb 39 and Cs 40 salts, respectively).

Aluminium hydrides
As for the Al hydrides, alane is found experimentally as a polymer with the formula (AlH 3 ) n that presents several polymorphic crystal structures based on vertex-and/or edgesharing AlH 6 octahedra, 41 while molecular Al 2 H 6 has only been isolated in solid hydrogen due to its instability. 42Only 21 structures with Al-H⋯H-Al intermolecular short contacts (H⋯H < 3 Å) were found.We have observed first that, in general, aluminium tends to form very short Al-H⋯Al contacts Fig. 2 Interaction topology of the 1,2-dicarbadodecaborane dimer as found in the crystal structure of bis(η 5 -pentamethyl-cyclopentadienyl)ytterbium(II) 1,2-dicarbadodecaborane (QEQFUN). 35The B-H⋯H-B short contacts at 2.60-2.70Å are shown as multi-band cylinders.
(1.77-2.34Å) that indicate incipient formation of Al-H-Al bridges, consistent with the fact that these contacts appear only for four-coordinated Al atoms.Other contacts appear only at distances above 3.10, with a van der Waals peak centered at 3.45 Å, in agreement with a van der Waals radii sum of 3.45 Å (see Fig. S2 in the ESI †).For instance, in the crystal structure of dimethylamine-alane, 43 the two Al⋯H-Al distances (d 1 in Fig. 3) are 2.07 Å, resulting in a trigonal bipyramidal environment around Al.There are, however, a few examples where actual Al-H⋯H-Al contacts can be found.[48][49][50][51]

Gallium hydrides
Digallane, H 2 Ga(μ-H) 2 GaH 2 , was prepared and characterized for the first time in 1989 after several years of search, 52,53 but, to date, no crystal structure for this compound has been reported.Its structure, however, has been determined in the gas phase by electron diffraction, 53 and also the structure of the related compound, Me 2 Ga(μ-H) 2 GaMe 2 , is known in the gas phase. 54Our CSD survey yielded 13 structures with Ga-H⋯H-Ga intermolecular contacts shorter than 3 Å.In 10 of the structures, C-H⋯H-Ga 55 or N-H⋯H-Ga 43 heteropolar dihydrogen contacts are found along with the Ga-H⋯H-Ga ones.A larger number of structures (35) were found with Ga⋯H-Ga contacts at distances around the sum of the van der Waals radii (3.52 Å).In our CSD survey, only one structure was found with a short In-H⋯H-In contact (CSD refcode SITNAK10), 56 which corresponds also to the only short In⋯H contact, 2.96 Å, to be compared with a van der Waals sum of 3.63 Å.Finally, no Tl-H⋯H-Tl short intermolecular contacts were found in our CSD survey.
It is worth remarking that we are aware of the fact that among the thousands of structures taken into account here for statistical purposes, not all the E-H⋯H-E short contacts correspond to dihydrogen interactions.Some H⋯H short distances are the consequence of other types of intermolecular interactions, e.g.hydrogen bonds, B-H⋯π, Al⋯donor, etc.In the next section, we focus on those particular H⋯H contacts that are not accompanied by any other interaction and, thus, can be considered as actual dihydrogen interactions.

Dimers of EH 3 molecules
In this section we present the results of our computational analysis of several existing and hypothetical dimers of the group 13 hydrides, comparing the results with experimental references when possible.We start by looking at the dimers of trivalent hydrides of the type EH 3 , in order to evaluate the energetics of the simplest possible homopolar B-H⋯H-B interactions, even if such molecules are unlikely to be found in crystal structures since the formation of diborane type dimers ( 6) is strongly favored (Table 1).Since InH 3 and TlH 3 are unstable and their corresponding dihydrides, In 2 H 6 and Tl 2 H 6 , have been the subject of solely theoretical studies to date, 57 we have decided not to include them in the present analysis.The results of calculations on the dimers of EH 3 (E = B, Al, Ga) hydrides with several topologies (6, 7a-7d) are presented in Table 1.We note that dimers with topologies 7a, 7b and 7c are not likely and the very low dissociation energies indicate that they are held together by nothing else than very weak dispersion forces. 61However, these very weak interactions could play a role in the stabilization of supramolecular assemblies when working cooperatively, as previously observed in alkanes.It must be stressed that the topology 7a, which allows for only one H⋯H interaction, corresponds to very weak attractive or even repulsive interaction and yields optimized H⋯H and E⋯H distances significantly longer than the corresponding sum of van der Waals radii. 29Interestingly, for topology 7d, dissociation energies are significantly larger than for other topologies and follow the trend Al > Ga > B (Fig. 4).For this topology, the Al⋯H and Ga⋯H distances are much closer to the van der Waals sum than the B⋯H one.In fact, it is clearly seen that the dissociation energy increases as the E⋯H contact distance approaches the van der Waals distance, again pointing to a possible role of the E⋯H interactions in the intermolecular dissociation energy (Fig. 4).On the other hand, the calculated H⋯H distances are longer than the van der Waals distance (2.40 Å) but similar to those found in alkane dimers with considerably strong interactions (dissociation energy for methane dimer is 1.75 kJ mol −1 with a H⋯H distance of 3.165 Å). 17 Also hypothetical E⋯E interactions may favor stronger attractive intermolecular interactions.Such contacts appear at distances clearly shorter than the sum of the van der Waals radii for both Al and Ga (E⋯E = 3.99 and 4.17 Å, radii sum = 4.50 and 4.64 Å, respectively) with topology 7d (Fig. 4).
In order to gain further insight into the intermolecular interactions, we have carried out an analysis of the electron density of the EH 3 dimers with topology 7d by means of the AIM method.For the BH 3 dimer, only a bond critical point (BCP) was found in a bond path between the two B atoms, suggesting the absence of an attractive H⋯H interaction.However, for AlH 3 and GaH 3 , four BCPs are found in each dimer associated with four H⋯H bond paths corresponding to the four H⋯H contacts.The electron densities at these BCPs (0.0042 for Al and 0.0032 for Ga) are similar to those reported for CH⋯HC contacts 17 and the positive values of the Laplacian are characteristic of closed-shell interactions.
Thus, the relatively large dissociation energies calculated for AlH 3 and GaH 3 dimers with topology 7d can be due to the combination of H⋯H dihydrogen and E⋯H (E = Al, Ga) interactions.It is worth remarking that our goal here is to evaluate the strength of noncovalent interactions (7a-d) rather than to compare them with the covalent bonded dimers (6).It is clear that the formation of a covalent 3c-2e bond will be strongly preferred when possible (see the large values of D in Table 1).Only when there is no possibility of covalent bonding, noncovalent interactions, several orders of magnitude weaker, could participate in the stabilization of crystal structures.

Dimers of diborane and digallane
In our calculations, diborane and digallane also appear to be capable of forming supramolecular dimers with the formula E 4 H 12 (E = B and Ga).Both terminal and bridging hydrogen atoms are able to establish dihydrogen contacts where two different intermolecular interactions can be distinguished (8a and 8b).The first one (8a) is observed in the crystal structure of diborane 62 and involves two terminal hydrogen atoms at 2.74 Å, while in the same crystal structure, the contacts between a terminal and a bridging hydrogen atom (8b) are still shorter (2.58 Å).In an earlier crystal structure of diborane, 63 a short contact occurs only between a terminal and a bridging hydrogen atom at 2.73 Å (8b).The 8b topology implies a heteropolar dihydrogen interaction due to the different character of the two hydrogen atoms, which makes us expect a stronger attraction.The existence of the terminal-bridging interaction has been recently detected in an AIM analysis of several boron hydrides. 64ometry optimizations were carried out for dimers of diborane, dialane and digallane in topologies, 8a and 8b.While diborane and digallane resulted in stable dimers (Fig. 5 and Table 2), no supramolecular dimers were found for dialane (see below).Key intermolecular distances and dissociation energies for diborane and digallane shown in Table 2 indicate that the E-H⋯H-E terminal-terminal interactions are similarly weak for B and Ga.The terminal-bridging interaction topology 8b led to different results for diborane and digallane, as can be seen in Fig. 5.The optimized diborane dimer presents three terminal-bridging H⋯H interactions (2.711-2.875Å) and two short B-H⋯B contacts involving both terminal (3.204Å) and bridging (3.057 Å) hydrogen atoms.In turn, the interaction topology of digallane after geometry optimization shows two H⋯H contacts (H t ⋯H t at 2.900 Å, H t ⋯H b at 2.937 Å) and two shorter Ga⋯H-Ga contacts (2.880 Å).The latter interaction can be responsible for the quite larger dissociation energy of the digallane dimer 8b compared to the diborane analogue (Table 2).This is consistent with the natural atomic charges calculated by means of a NPA analysis of the H t ⋯H b digallane dimer (Fig. 5).In such a structure, Ga and terminal H atoms involved in the interaction present atomic charges of 0.957 and −0.319, respectively, whereas analogous atoms in diborane show atomic charges of 0.021 and −0.038.
It is also interesting to observe how the natural atomic charges are modified when forming dimers with respect to the monomers.In diborane, terminal H atoms have charges of −0.037 while the same atom's charge is −0.043 when involved in an H⋯H contact in the dimer.However bridging H and B atoms show the same charge both in the monomer and dimer.In digallane, charges for Ga, terminal H and bridging H are 0.942, −0.295 and −0.351, respectively.These charges are significantly modified in the atoms that establish interactions in the dimer, as can be seen in Fig. 5. Therefore, it seems plausible that an electrostatic effect, together with a small amount of charge transfer is contributing to the large dissociation energy observed in the digallane dimer.The different nature of the bridging hydrogen atoms in diborane and digallane is also worth noting.While their charges are slightly positive in diborane (0.043-0.045), they present a more hydridic character in digallane (natural charge of −0.351), in good agreement with the lower electronegativity of Ga relative to B and H.
Taking the four data sets in Table 2, a clear correlation appears between the shortest E⋯H distance in each dimer and its dissociation energy (Fig. 6), while no correlation at all can be found between the H⋯H distance and the dissociation energy.This latter fact does not mean that H⋯H interactions are irrelevant.In fact, an AIM analysis of the four dimers in Table 2 shows that BCPs are found between hydrogen atoms and no bond path appears that indicates an E⋯H (E = B, Ga) interaction.A dimer with topology 8a presents a BCP between two terminal hydrogen atoms (ρ = 0.0026), both for B and Ga.For topology 8b, the (B 2 H 6 ) 2 dimer shows two BCPs (ρ = 0.0037, 0.0041) between the terminal and bridging hydrogen atoms while the (Ga 2 H 6 ) 2 dimer shows only one BCP (ρ = 0.0092) between two terminal hydrogens (the bond path corresponds to the stripped cylinder in Fig. 5).
In the case of dialane, optimizations starting from either interaction topology led to dimerization of the molecule, resulting in a tetranuclear structure Al 4 H 12 (see Fig. S4 in the ESI †).The new formed tetramer showed all three Al-Al distances very close to 2.6 Å, and H⋯H contacts between the  a An extra H t ⋯H t contact at 2.900 Å appears after the geometry optimization of H t ⋯H b interaction topology (white and blue multi-band cylinders in Fig. 5).original molecules, of both terminal-terminal (H t ⋯H t ) and terminal-bridge (H t ⋯H b ) nature, with distances of 2.29 Å and 2.39 Å, respectively.The dissociation energy was calculated to be as high as 87.5 kJ mol −1 .The tendency of aluminium to form short Al-H⋯Al bonds in alane derivatives has been discussed above and previously stated to be greater than that of gallium.For example, dimethylamine-alane presents an Al-H⋯Al distance of 2.07 Å between neighbouring molecules in its crystal structure while in dimethylamine-gallane, the shortest Ga-H⋯Ga distance is 2.89 Å, both with the same topology. 65

Dimers of arachno and nido polyboranes
Let us now try to extend our analysis to some representative members of the large and important family of polyhedral boron compounds, which includes carboranes and metallacarboranes. 668][69][70] The geometries of these three compounds were directly retrieved from their crystal structures, finding two different conformations between molecular units that are associated to two different H⋯H interaction patterns for each of them.Then, the two dimeric topologies were fully optimized for each case (Fig. 7).All calculated dissociation energies were positive (see results in Table 3) and, thus, we have identified the interactions between monomers as attractive.The dissociation energies of the dimers of B 10 H 14 are remarkably high, with H⋯H contacts between bridging hydrogen atoms as short as 2.23 Å.We have observed that, in general, H⋯H contacts involving one bridging and one terminal hydrogen atom (i.e.non-homopolar interactions) are shorter than those involving two terminal or two bridging hydrogen atoms.We note that, in general, the molecules of nido and arachno boranes in their crystal structures tend to arrange in such a way that the number of terminal-bridging interactions are maximized.In fact, these interactions are expected to be more energetic than the terminalterminal ones due to the different character, acidic and hydridic, of the two H atoms involved.Our calculated H⋯H distances are in very good agreement with the experimental values.In the combined X-ray and neutron diffraction structure of dodecaborane (CSD refcode FUYYIH02), 30 we find   72 ).In all three cases the back to face (b-f ) arrangement of a dimer that is found in the crystal packing appears to be stable toward dissociation and the dissociation energy increases with the size of the borane.The distribution of B-H⋯H angles in these optimized structures follows the same general pattern seen in Fig. 1b.Optimization of dimers with the topology found for sideways interactions in the crystal structures, however, lead to other topologies that imply more hydrogen atoms into the intermolecular interaction.For instance, B 5 H 9 gives a bottom to side topology for a dimer (Fig. 7) that can be described as resulting from the tilting one of the molecules from the b-f geometry, in such a way that two more H⋯H short contacts are introduced without disrupting the 4 : 1 contacts, therefore increasing the dissociation energy from 11.7 to 16.5 kJ mol −1 .Obviously, we cannot pretend the interaction topology to be the same as that found in the solid state, because in the latter case somewhat weaker interactions can be compensated by the possibility of extending them in one or two dimensions.It is nevertheless rewarding that in the three cases studied here the b-f interactions that correspond to the minima for the dimers are consistent with those found in one direction of the crystal structure.In this way, our results clearly indicate that (i) those interactions are attractive and (ii) they increase with the size of the borane.The latter result is consistent with our previous finding that the C-H⋯H-C interactions between polyhedranes are enhanced by (a) the degree of substitution of the supporting C atom and (b) the loss of pyramidality of the C-H group as the size of the polyhedron increases. 17

Dimers of closo dicarbadodecaborane
Unfortunately, most of the members of the closo polyboranes family are charged species, which makes them unsuitable for a computational study of noncovalent intermolecular interactions that are undoubtedly much weaker than ionic forces.A neutral species of this family is the icosahedral 1,2-closo-dicarbadodecaborane, in which two B-H units are substituted by C-H ones.We present in Table 4 the geometrical parameters of four selected dimers of dicarbadodecaborane retrieved from four different CSD crystal structures.The dissociation energies (D in Table 4) are calculated on the geometry of only two interacting molecules (omitting all others) with the desired interaction topology as found in the crystal structures without further optimization.We are aware that more molecules may be present in the unit cell of each crystal structure that contribute to the stabilization of the ensemble and, thus, the calculated dimers may not be minima of the potential energy surface.However, in this way, by avoiding geometrical rearrangements like those observed in the optimization of nido polyborane dimers, we can have a picture of the interaction strength between two molecules as they are found in a crystal structure.
The calculated dissociation energies for dicarbadodecaborane dimers with topologies 1 : 3 (1) and 2 : 2 (2) are close to those found for dimers of arachno B 4 H 10 with a similar number of H⋯H interactions (Tables 3 and 4).Large nido dimers like B 10 H 14 present larger dissociation energies than the closo dicarbadodecaborane because the former can establish much more H⋯H and B⋯H contacts (for the latter, see Table 3).The dicarbadodecaborane dimer with 1 : 2 interaction topology (BEZTUW) shows in its crystal structure two heteropolar C-H⋯H-B contacts at 2.97 and 2.99 Å, in addition to the two B-H⋯H-B ones, which might also contribute to the total dissociation energy of 10.34 kJ mol −1 .On the other hand, the dissociation energy for TUQTEE, 4.26 kJ mol −1 , with only one  H⋯H interaction (2.90 Å) is much larger than those obtained for the 1 : 1 dimers of BH 3 and B 2 H 6 , 0.30 and 0.64 kJ mol −1 , respectively.This magnification of the dissociation energy by a polyhedral skeleton was already observed for alkanes, among which the dimer of dodecahedrane C 20 H 20 presented a dissociation energy of 4.14 kJ mol −1 , in contrast with the minute 0.52 kJ mol −1 found for methane with the same 1 : 1 interaction topology. 17An explanation to this remarkable behaviour is that E-H⋯H-E interactions in large systems (e.g.dimers of polyhedral alkanes or boranes) are stabilized by a different mechanism than those of small ones (e.g.CH 4 or BH 3 dimers).
While the latter are governed by interactions between oscillating dipoles E δ+ -H δ− ⋯H δ+ -E δ− and E δ− -H δ+ ⋯H δ− -E δ+ and, thus, the stabilization can be attributed to "classic" dispersion, the great stability of larger dimers can be explained in terms of a reorganization of the bonding electrons of two E-H bonds to create alternative H⋯H and E⋯E "bonds" and also to the charge transfer between the two interacting units to establish long-range E⋯H interactions. 18

Conclusions
We have carried out a systematic theoretical study of homopolar dihydrogen interactions in group 13 compounds.Several families of compounds containing E-H (E = B, Al and Ga) interacting units have been analysed, from the smallest trivalent hydrides EH 3 to the larger polyhedral boranes.B-H⋯H-B contacts are by far the most abundant ones and involve small dissociation energies, comparable to those associated to previously studied C-H⋯H-C contacts. 17In some cases, B-H⋯H-B interactions, combined with B-H⋯B ones, can be associated to surprisingly large dissociation energies, as for example in the calculated dimer of the nido polyborane B 10 H 14 (D = 27.14 kJ mol −1 ).We have identified two factors that strengthen the interaction: (a) the presence of a polyhedral skeleton associated to the E-H interacting unit (for example, a 1 : 1 H⋯H interaction is one order of magnitude larger in a dicarbadodecaborane dimer than in a diborane one) and, (b) a high number of H⋯H and E⋯H contacts between two monomers working cooperatively.This cooperative character of the intermolecular contacts was previously found in alkanes, 17 and now also observed in boron hydrides.We note that the 2 : 2 staggered (7d) and the 3 : 3 staggered topologies are the most stable ones for EH 3 and CH 4 noncovalent dimers, respectively, but for E = B the dissociation energy is almost negligible.In general, the calculated dissociation energies are larger for interaction topologies that involve multiple H⋯H and E⋯H short contacts (7d and 8b) for B, Al and Ga.This is in good agreement with the experimental data, since in many crystal structures, the real scenario is a combination of several short contacts of H⋯H, H⋯E and even E⋯E types, all contributing to the global stability.Based on our theoretical analysis, we have observed that the strength of the interaction within group 13 follows the trend Al > Ga > B, with Al showing a remarkable tendency to form very short Al⋯H bonds.

Fig. 1
Fig. 1 (a) Distribution of B-H⋯H-B distances shorter than 5 Å as found in the CSD.The blue bars represent all B-H⋯H-B contacts, the transparent bars with grey outline represent only the contacts involving two hexacoordinate B atoms, and the solid grey bars correspond to all the crystal structures excluding artificially fixed B-H bond lengths and R < 5%.(b) Distribution of the B-H⋯H angles for the H⋯H contacts (a) within the van der Waals peak (i.e., shorter than 3.5 Å) in blue, and the same distribution for the two distinct angles α and β (Fig. 2) in the 3 : 1 contacts between B 3 H 3 and H-B groups.See the ESI † for the angular distribution of other contact topologies and ref. 29 for a more detailed definition of the van der Waals peak concept.

Fig. 3
Fig.3Representation of the two short Al-H⋯Al contacts (dashed lines) as found in dimers of some AlH 3 adducts with Lewis bases, for example dimethylamine-alane.43

Fig. 4
Fig. 4 Dissociation energy as a function of the difference between the E⋯H (E = B, Al and Ga) intermolecular distance and the sum of the van der Waals radii in EH 3 dimers with interaction topology 2 : 2 (7d).

Fig. 5
Fig. 5 Optimized structures of dimers of diborane (above) and digallane (below) with terminal-bridging interaction topologies.Short H⋯H and E⋯H contacts are shown as multi-band cylinders and dashed lines, respectively.Natural atomic charges are given next to each atom.

Fig. 6 5 .
Fig. 6 Dissociation energy as a function of the difference between the H⋯H (red squares) and E⋯H (blue circles) intermolecular distances and the sum of the corresponding van der Waals radii in (E 2 H 6 ) 2 (E = B and Ga) dimers, with the optimized structures shown in Fig. 5.

Fig. 7
Fig. 7 Optimized geometries and shortest H⋯H contacts in dimers of an arachno and two nido polyboranes.The short intermolecular contacts are shown as striped cylinders, blue and white for H t ⋯H t , grey and white for H t ⋯H b , and red and white for H b ⋯H b contacts.Topological code: b-f = bottom-face, b-s = bottom-side, b-b = bottom-bottom.

Table 1
Intermolecular distances and dissociation energies (D) for group 13 trivalent hydride dimers.BSSE-corrected dissociation energies (D) are given in kJ mol −1 .The dissociation energies of the diborane-type E 2 H 6 dimers into two EH 3 monomers are given for comparison.Experimental dissociation energies for diborane and dialane are given in parenthesis, and the sum of van der Waals radii is given in square brackets 7d 3.656 [3.52] 3.244 [2.40] 3.77

Table 2
Intermolecular distances (in Å) and dissociation energies for group 13 dimers of E 2 H 6 molecules with terminal-terminal and terminal-bridge interaction topologies 8a and 8b.BSSE-corrected dissociation energies (D) are given in kJ mol−1 H t ⋯H t and H t ⋯H b contacts at 2.55 and 2.41 Å, respectively.Other dodecaborane derivatives also present H t ⋯H b contacts at 2.41-2.46Å (CSD refcode CCOBOR 71 ) and H t ⋯H t ones at 2.43-2.50Å (CSD refcode PIJTUX

Table 3
Number of calculated short B⋯H (<3.890Å) and H⋯H (<3.000Å) contacts between two terminal H t ⋯H t , two bridging H b ⋯H b or one terminal and one bridging H t ⋯H b hydrogen atoms in nido and arachno borane dimers.Interaction topologies involve the bottom (b), face (f ) or side (s) of the monomers.The minimum and maximum contact distances are given in parenthesis (Å) for all cases; BSSE-corrected dissociation energies (D) are given in kJ mol−1We consider contacts not exceeding the sum of the van der Waals radii by 25%.The van der Waals radii proposed by Alvarez 29 are 1.20 Å for hydrogen and 1.91 Å for boron.

Table 4
Experimental geometrical parameters and calculated dissociation energies (D in kJ mol −1 ) for the dimers of closo polyhedral Open Access Article.Published on 09 February 2017.Downloaded on 8/13/2019 9:06:03 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.