Dihydrogen vs . hydrogen bonding in the solvation of ammonia borane by tetrahydrofuran and liquid ammonia †

The solvation structures of two systems rich in hydrogen and dihydrogen bonding interactions have been studied in detail experimentally through neutron diffraction with hydrogen/deuterium isotopic substitution. The results were analysed by an atomistic Monte Carlo simulation employing refinement to the experimental scattering data. The systems studied were the hydrogen storage material ammonia borane (NH3BH3, AB) dissolved in tetrahydrofuran (THF), and liquid ammonia (NH3), the latter in which AB shows unusually high solubility (260 g AB per 100 g NH3) and potential regeneration properties. The full orientational and positional manner in which AB–AB, AB–THF and AB–NH3 pairs interact with each other were successfully deciphered from the wide Q-range total neutron scattering data. This provided an unprecedented level of detail into such highly (di)hydrogen bonding solute–solvent interactions. In particular this allowed insight into the way in which H–B acts as a hydrogen bond acceptor. The (di)hydrogen bonding was naturally determined to dictate the intermolecular interactions, at times negating the otherwise expected tendency for polar molecules to align themselves with anti-parallel dipole moments. Several causes for the extreme solubility of AB in ammonia were determined, including the ability of ammonia to (di)hydrogen bond to both ends of the AB molecule and the small size of the ammonia molecule relative to AB and THF. The AB B–H to ammonia H dihydrogen bond was found to dominate the intermolecular interactions, occurring almost three times more often than any other hydrogen or dihydrogen bond in the system. The favourability of this interaction was seen on the bulk scale by a large decrease in AB clustering in ammonia compared to in the dihydrogen bond-less THF.


Introduction
A classic example of an important intermolecular interaction is the hydrogen bond. 1 This follows the pattern X-HÁ Á ÁY, where the donor atom X, being more electronegative than hydrogen, causes said hydrogen to possess a partial positive charge.The acceptor, Y, is an atom capable of stabilising this charge, usually taking the form of an atom with a free lone or nonbonding pair of electrons, commonly a nitrogen or oxygen.Hydrogen bonds are common, and influence the properties of molecules as small as water, up to those as large and complex as proteins and DNA.
A somewhat rarer type of intermolecular interaction is dihydrogen bonding.This is quite similar to regular hydrogen bonding in terms of the donor, but the difference lies in that the acceptor is no longer a single atom.A dihydrogen bond instead follows the pattern X-HÁ Á ÁH-Y, where the H bound to Y possesses a partial negative charge, caused normally by the relatively lower electronegativity of Y.This leads to a weak electrostatic attraction between the two hydrogens involved.Dihydrogen bonds have been observed within transition metal hydrides, playing a role in processes such as proton transfer and H 2 evolution. 2Dihydrogen bonding has also frequently been observed in main group hydrides, especially those containing boron, where they have been shown to enable stereochemical control of reactions like reduction. 3,4mine-boranes are a family of compounds known to be rich in dihydrogen bonding, owing to the hydridic and protonic natures of the hydrogens bound to boron and nitrogen respectively.6][7][8][9] In AB, the dihydrogen bonds were distinguished from other intermolecular interactions (such as van der Waals forces) and classified as a type of hydrogen bond using ''atoms in molecules'' theory. 10The H N Á Á ÁH B dihydrogen bond is responsible for many of the properties of AB, including the high melting point (110-114 1C, cf. the isoelectronic ethane À182.8 1C). 11B has received wide-spread attention due to its potential as an air-stable and safe hydrogen storage material.AB contains 19.6 wt% hydrogen, approximately 13-16 wt% (2-2.5 eq.H 2 per AB) of which can be extracted by heating to mild temperatures (o150 1C) (Fig. 1).The dihydrogen bonds have to be overcome in order for the material to melt, which occurs simultaneously with the onset of hydrogen release. 11In addition to these intermolecular interactions, homopolar H B Á Á ÁH B interactions have been theorised to exist during the dehydrogenation process, due to an isotopic substitution experiment (which observed, for example, D 2 on dehydrogenation of NH 3 BD 3 ). 12 factor holding back the use of AB is in the difficulty of recycling the waste material post-dehydrogenation.This waste consists of a variety of partially hydrogenated hexagonal boron nitride-like products.Many attempts have been made to efficiently regenerate AB waste.[13][14][15][16][17] A promising one-pot route involves treating the waste with hydrazine in liquid ammonia, which was shown to regenerate an AB waste surrogate successfully in 24 hours at 40 1C.14 A further curiosity in terms of treating AB waste surrogates with liquid ammonia was demonstrated by Davis et al., who showed that some Lewis bases (specifically, pyridine and ammonia) were able to scavenge some of the residual high energy B-H bonds in the waste and rearrange them to form BH 3 .If the Lewis base used was ammonia, this led to direct regeneration of some of the waste back into AB.18 AB also has an extremely high solubility in ammonia (260 g AB per 100 g NH 3 at 298 K, cf. 25 g AB per 100 g tetrahydrofuran).11 This behaviour raises the question of how AB interacts with ammonia in solution, and has implications for future attempts at the regeneration of AB waste.
While there have been studies into the dynamics of dihydrogen bonding using IR and NMR spectroscopy, there has not yet to our knowledge been a detailed empirical analysis of the full structural and orientational way in which dihydrogen bonding occurs in solution. 19,20][23][24][25][26][27] In addition to the AB-ammonia interactions, a more fundamental interest in the nature of solvation in such a (di)hydrogen bondrich environment sparks our study of these systems.Insight into the AB-AB intermolecular interactions is a further goal for this work.Due to the concurrent release of hydrogen when AB melts, the interactions within molten AB cannot be directly studied; a concentrated solution of AB in liquid ammonia, especially considering the structural similarities, is probably the best analogue for studying such a system.
In this paper the solvation structure of AB in both liquid ammonia and tetrahydrofuran (THF), another common solvent, has been investigated through the use of neutron diffraction.9][30] Isotopic substitution of both the solvent and solute hydrogens for deuterium was performed to obtain additional, complementary constraints for the data refinement process.The results were compared with those from solid state diffraction and ab initio calculations of dimers, providing important insight into the distribution of dihydrogen bonding structures in the liquid state.A description of the theory behind the neutron diffraction experiments and the EPSR procedure is provided in the ESI.†

Results and discussion
The neutron scattering experiments were performed on 7 samples for AB in THF and 3 samples for AB in NH 3 (details in Experimental section).Table 1 shows the seed parameters used for the EPSR simulations, with the atom nomenclature shown in Fig. 2. The Lennard-Jones parameters for AB were adapted from the CHARMM parameters for methylamine used in Chen et al., as parameters for AB do not exist to our knowledge. 31It is important to mention that the seed potentials are only a starting point, and differences in structural behaviour induced in the simulation would be compensated for by the refinement.The AB charges were sourced from a Mulliken analysis of AB itself. 32For ammonia, the nitrogen parameters were taken from Gao et al., with the hydrogen values adjusted to match the AB N-Hs. 33[36][37][38] Fig. 1 The thermal dehydrogenation of AB leads to a multitude of unsaturated cyclic products.The normalised structure factors for the AB in THF and AB in ammonia experiments are shown in Fig. 3 and 4 respectively, together with the simulated EPSR fits.Scattering for Q o 2 Å À1 is dominated by the intermolecular interactions, weighted by each atom's neutron scattering length and concentration.The opposite occurs at high Q, with the intramolecular interactions dominant.
The EPSR fits are of good quality, with the deviations at Q o 0.5 Å À1 caused most likely by errors in the inelasticity correction, which would explain the larger magnitude in hydrogen (rather than deuterium) rich samples. 39,40he Fourier transformation of the total structure factors produces the total radial distribution functions, shown in Fig. S1 and S2 (ESI †).These show that the intramolecular bond lengths are modelled accurately in our simulations.
More specific structural information can be obtained by looking at the partial radial distribution functions (g(r), RDF) for atom pairs of interest in our systems.

Ammonia borane-ammonia borane interactions
We have used the output of the EPSR simulations to study both hetero-and homopolar interactions between AB molecules.Given the likelihood for AB-AB interactions to be less rigid in solution than in the solid state (and arguably more similar to melted/dehydrogenating AB), this allows an understanding of the interactions suggested to exist during dehydrogenation. 12ig. 5 shows the partial RDF for heteropolar H N Á Á ÁH B and homopolar H B Á Á ÁH B interactions.The sharp peak at 1.85 Å in the heteropolar function suggests the presence of strong dihydrogen bonding between AB molecules.This distance is smaller than the sum of the van der Waals radii for two hydrogen atoms (2.4 Å) and characteristic of, though slightly shorter than, the dihydrogen bonding present in the solid state (2.02 Å). 11,20 This structural motif appears in both ammonia and THF.In ammonia however, the magnitude of the peak is much lower suggesting AB-AB interactions are less numerous (both graphs were normalised to the bulk concentration of AB).The second peaks at 3.4 Å are not due to a second coordination shell, but the non-bonding H B s on the same AB molecule responsible for the 1.85 Å peak.
The difference in the number of AB-AB interactions between the two solvents can further be seen in an analysis of the AB-AB clustering (a cluster defined as having an H N Á Á ÁH B distance of 2.7 Å or less) present in the solutions, shown in Fig. 6.Overall 86% of the AB molecules in THF were in a cluster of size two or    more, compared to 72% in liquid ammonia.This lower overall clustering in ammonia also manifested as a tendency towards smaller clusters, with the largest involving about 20 AB molecules.In THF however, larger clusters were often observed containing over 30 AB molecules.Each H N on average dihydrogen bonds to 0.701 H B s in THF, compared to only 0.438 in ammonia.This was determined from the average coordination numbers obtained by integrating the partial radial distribution functions.The results are shown in Table 2.More detailed probability distributions of the dihydrogen bonding in the first coordination shell (o2.7 Å) are shown in Fig. 7.These show that although the average coordination number is less than one, situations where multiple H N s are dihydrogen bonded to each H B are not uncommon, more so in THF than in liquid ammonia, reflecting the increased clustering.
The lower clustering in ammonia suggests that AB-ammonia interactions are relatively favourable compared to AB-THF.This may be linked to the ability of ammonia to partially break down the polymeric structure of the waste material, redistributing the residual B-H bonds to form AB as mentioned before. 18This is a process that would be promoted by stronger AB-solvent interactions, and is not observed in THF.
In Fig. 5, the peak signifying homopolar H B Á Á ÁH B interactions does not appear until 4.5 Å.This long range suggests that there are no direct correlations in these solutions.The density observed in these functions at this distance is likely instead due to the heteropolar interactions.This suggests that the presence of D 2 during the dehydrogenation of NH 3 BD 3 , as observed by Wolstenholme et al., was more likely to be a result of hydrogen scrambling than homopolar coordination. 12ig. 8 shows the spatial density functions (SDF) showing where neighbouring AB molecule density is located around a central AB in ammonia (similar SDFs are seen in THF).The density is broadly located side-on to the central AB rather than end-on, and adopts the three-fold symmetry expected for intermolecular interactions primarily based on the hydrogen bonding sites.Fig. 9 shows the angular RDF for the AB-AB interactions, with y corresponding to the angle between the AB dipole moments.The distance r is that between H N and an H B .The function shown is for AB dissolved in THF, the same function calculated in liquid ammonia (Fig. S3, ESI †) showed the same features.This implies that the only difference caused by the solvent is in the number of AB-AB interactions rather than the structure.A strong bias towards anti-parallel dipole moments is observed.The side-on configuration shown in Fig. 9 is the only configuration where (heteropolar) dihydrogen bonding and anti-parallel dipole moments can occur together.The second peak, at 3.5 Å, is assigned to the distance between H N and the other H B s on the dihydrogen bonding AB.This journal is © the Owner Societies 2018 In the angular RDF in Fig. S4 (ESI †) with y equal to the angle between the N-H bonds, a preference for both a parallel and anti-parallel orientation of the N-H bonds is seen.The resulting AB-AB pairs fulfilling all of the above criteria are shown in Fig. 9.The N-H N bonds point towards the B-H B bonds, rather than directly at the H B s.This structure, observed before in theoretical studies of dimers, has been suggested to be favourable by allowing the protonic hydrogen to get nearer to the relatively negative boron atom. 41As a result, it was proposed that the s-bond itself was the acceptor in the dihydrogen bond rather than just the hydridic hydrogen. 20This may also explain how multiple dihydrogen bonds can form with a particular H B , as H N s can bond from many angles surrounding the s bond rather than from only pointing directly at the H B (regular hydrogen bonds typically being linear, see Section 2.4).
In considering the H N density in Fig. 8 near to the nitrogen end of AB, we note that there is no close-range H N Á Á ÁH N homopolar correlation either, the peak in the partial RDF appearing at approximately 4.2 Å (Fig. S6, ESI †).

Ammonia borane-tetrahydrofuran interactions
The H to O partial RDFs for AB dissolved in THF are shown in Fig. 10.Only H N shows strong interactions to O1 with two contacts at 2.05 and 3.44 Å respectively.The first is attributed to a direct H N Á Á ÁO1 hydrogen bond while the second peak is due to a longer range interaction between the oxygen and another H N on the same nitrogen.As expected, given the hydridic nature of H B , the H B Á Á ÁO1 RDF shows only a much longer range interaction, indicated by a broad peak centred at 4.83 Å. Fig. 11 shows the SDFs plotting the regions of highest probability density around the AB and THF in which the other molecule is located.These show that the density is centred around the hydrogen bonding sites O1 and H N .Fig. 12 shows the angular RDF, where the angle y is that between the AB and THF dipoles, and the distance r is that between H N and O1.A roughly equal, high probability is seen for orientations with y between 90 and 1801.Together with the SDFs and other angular RDFs (Fig. S7-S9, ESI †), we can determine that THF occupies a range around the nitrogen end of AB, with O1 always pointing at H N .This can involve a disruption of the usual tendency for pairs to align with anti-parallel dipole moments, which speaks to the strength of the hydrogen bonding interactions.It could also be a result of the increased degree of disorder present in liquid samples.
On average, each H N coordinates to 0.361 O1s, suggesting it is about twice as likely for an AB (via H N ) to bond to another AB (0.701) than to a THF molecule.A more detailed distribution of the AB-THF coordination is shown in Fig. 13, which also shows that it is much less likely for each H N to be within hydrogen bonding distance of more than one O1, than with multiple H B s (Fig. 7).Looking at the distribution of H N s around O1s also reveals that the THF oxygen can be within hydrogen bonding range of two H N s at once, suggesting that both of the oxygen non-bonding electron pairs could be involved in hydrogen bonding at the same time.

Ammonia borane-ammonia interactions
Fig. 14 shows the AB H to ammonia N A and H A partial RDFs.H A predominantly interacts with H B , as can be seen by the defined peak at 1.98 Å.This is a slightly longer dihydrogen bond than the corresponding bond between two ABs (H N Á Á ÁH B , 1.9 Å), suggesting it may be weaker.H N , however, does not seem to dihydrogen bond directly with H A , the 'shoulder' peak at 2.7 Å instead likely due to AB bonded via the H B , as was observed in Ricci et al. 42 The opposite was determined for bonding to N A , with H N showing a clear peak at 2.19 Å. H B however does not show any bonding to N A before the first peak at 2.88 Å, likely again caused by bonding to the other end of AB.The extra dihydrogen bonding interactions observed between AB and ammonia, and not present in THF, are likely responsible for the large difference in solubilities.
The average first shell (o2.5 Å) coordination number of N A around H N , and of H A around H B is 0.798, and 1.905 respectively.The latter is remarkably high, possibly due to the large number of H A s in the system compared to N A s. Something of note though in comparing to THF, is the effective doubling of the average coordination number of O1 around H N (0.361) in regards to N A .While it is likely this is at least partly due to the higher concentration of AB in ammonia, the smaller size of the ammonia molecule may also make it less sterically unfavourable for multiple H N s on an AB to be coordinated.Fig. 15 shows the more detailed analysis of the coordination, which reveals just how frequently H B Á Á ÁH A dihydrogen bonds occur.As many as five H A s were observed to be within dihydrogen bonding distance of one H B , eclipsing the coordination observed in THF.The coordination is also weighted more towards one or two bonds rather than zero as was observed in the AB-AB coordination in THF.
The SDFs for AB-ammonia pairs are shown in Fig. 16.The density around AB largely mirrors that seen in the previous pairs, a three-fold symmetrical pattern of density around the H N s, and between the H B s.In this case, the lobes of density between the H B s also extend down to meet end-on to the boron atom, unlike in the AB-AB pair.This could reflect the extreme prevalence of H B Á Á ÁH A dihydrogen bonds, and the smaller size of the ammonia molecule relative to AB enabling it to bond from positions where the larger AB could not.This large volume    of probability density may also be the reason that a coordination number as high as 5 can be observed, albeit very rarely.The location of the probability density being between the H B s rather than pointing at them again highlights that it is the B-H B s bond that is acting as the dihydrogen bond acceptor.
The angular RDFs were split into those calculated based on H N Á Á ÁN A , and H B Á Á ÁH A interactions.These are shown in Fig. S10 and S11 (ESI †) respectively.Orientations that satisfy the ARDFs and the SDFs, are shown in Fig. 17.The additional angular RDFs needed to narrow down the most probable orientations are shown in Fig. S12 and S13 (ESI †).

Directionality of dihydrogen bonding in solution
The hydrogen and dihydrogen bonding angles X-HÁ Á ÁY observed throughout the simulations, where X is the hydrogen bond donor and Y the acceptor, were determined with the resulting probability distributions shown in Fig. 18.All the interactions are predominantly linear, matching what would be expected for hydrogen bonds. 1 The dihydrogen bonds however, show wider peaks than the regular hydrogen bonds, indicating more flexibility.This could reflect a slightly weaker bond, or be a product of the larger number of these interactions forcing some to adopt a less ideal orientation.The other angles XHÁ Á ÁY-A, where A is an atom covalently bound to Y, were also determined and shown in Fig. 19.The peaks are much wider for these angles, reflecting less rigidity.The H N Á Á ÁH B -B angle however predominantly appears at 901.This is due to the dihydrogen bond acceptor formally being the H B -B s bond rather than just the hydridic hydrogen, as was seen in the SDF (Fig. 8).For the H A Á Á ÁH B -B dihydrogen bond, a slightly larger angle of 1101 is seen.This may be due to H N having a slightly more positive charge than H A , and therefore experiencing a stronger attraction to the electron density in the s bond.
Analysing the angles of the (di)hydrogen bonds also enables us to distinguish these from other intermolecular interactions (such as van der Waals forces).Fig. S14 (ESI †) shows the angular distribution for the non-(di)hydrogen bonding B-H B Á Á ÁH B .The distribution is very flat, with the reduction in probability at o1201 likely due to steric reasons.This is very different to what we see in the (di)hydrogen bonds, which all seem to show directionality.

Comparison of ammonia borane-ammonia results to a theoretical approach
After observing that (di)hydrogen bonding to both ends of the AB molecule by ammonia was occurring, the possibility was raised that a single ammonia could bond to both ends of the    same AB.Indeed in theoretical studies run on systems of AB dissolved in ammonia by Kar and Scheiner, their fully relaxed AB-NH 3 dimer (Fig. 20) showed H N Á Á ÁN A hydrogen bonding as well as H A dihydrogen bonding to two H B s. 43 Looking for pairs where (di)hydrogen bonding occurs on both ends of an AB molecule to one NH 3 in our empirically refined simulation showed that this occurred 32% of the time (out of all bonding pairs).Dihydrogen bonding was determined by the atom pairs satisfying distance constraints of o2.7 Å for H N Á Á ÁN A and o2.5 Å for H B Á Á ÁH A .These constraints were chosen to include the full breadth of the corresponding peak in the partial RDFs.Angular constraints were not included as any orientation which involves bonding to both sides of the AB would be radically different to what we saw in the prior analysis.Searching specifically for the bonding pattern in the theoretically optimised dimer (one H A dihydrogen bonding to two H B s, and N A hydrogen bonding to one H N ), showed it as present only 4.2% of the time.The average geometric parameters were also quite different (H N Á Á ÁN A lit.  43 The closest match was the last angle, which the authors claimed to show energy minima at 100 and 1801; in this work we see broad corresponding probability peaks at about 90 and 1801 (Fig. 19).The broadness of the peaks and the relatively small differences between experimental and theoretical results reflects the disordered nature of the liquid samples.The theoretical model also only had to minimise the energy of two molecules, while in the liquid multiple interactions had to be globally minimised.The chelate effect would suggest that solvation structures where one ammonia bonds to both ends of AB would be favoured, considering the increase in entropy gained from displacing other ammonia molecules.Nevertheless the probabilities were quite low, this may be because such orientations diverge from the most favourable angles (Fig. S10 and S11, ESI †); the multiply bonded pairs have H A -N A running parallel to the AB dipole rather than anti-parallel.

Conclusions
The solvation structure of AB by THF and liquid ammonia has been determined via refinement of a Monte Carlo simulation to neutron diffraction data, constrained by complementary data sets of isotopically substituted samples.The large capacity of AB to hydrogen bond was found to dominate the solute-solvent interactions.In THF, hydrogen bonding via the oxygen was usually observed to happen only once per THF, with instances where both non-bonding electron pairs were involved occurring five times less frequently.The low number of hydrogen bonds formed with THF, along with the larger and more rigid size relative to another AB, caused AB-AB contacts to be favoured over AB-THF.This manifested as a noticeable increase in AB clustering.
The extreme solubility of AB in liquid ammonia could then be seen in the sheer number of (di)hydrogen bonds formed on average between AB and ammonia.The first shell coordination number of N A around H N was more than double that of the THF oxygen around H N .This could be due to the increased availability of the nitrogen lone pair compared to the oxygen non-bonding electron pairs.The H B Á Á ÁH A dihydrogen bond in particular was prevalent, with a coordination number roughly 2.5 times greater than even N A around H N .This also came with a relative decrease in AB clustering.
The structure of the dihydrogen bonding between AB and either ammonia or another AB in solution was determined to be very directional, with the N-H N/A Á Á ÁH B angle centred around 1801.This was still slightly less rigid than regular hydrogen bonds.The frequency of these interactions was attributed to the preference for a 901 H N/A Á Á ÁH B -B angle, which enabled the bonding to occur from a much larger volume of space surrounding each H B .This may also be the reason for the lower rigidity of N-H N/A Á Á ÁH B , as some distortion of the linearity to point towards the B-H B s bond would be expected.
The fact that ammonia can hydrogen bond to both ends of the AB molecule undoubtedly leads to an increased solubility of AB compared to in THF.In comparing to AB-AB contacts however, which can also occur on both ends of the molecule, AB-ammonia contacts were still preferred.We determined this could be due to the smaller size of the ammonia molecule opening up new volumes of space in which hydrogen bonding could occur from.This was mainly observed in an end-on position next to the AB boron.
It is the combination of these factors that cause the high solubility of AB in liquid ammonia.An interesting test would be to study in a similar way the solvation structure of AB in water, as water also has a small size and could bond to both ends of AB.The solubility in water (33.6 g AB per 100 g solvent) however is more similar to that in THF (25 g AB) than ammonia (260 g AB). 11his could perhaps place more weight on the availability of the lone pair as well as the number of hydrogens involved in terms of factors contributing to solubility.It could also be possible that the three-fold symmetry of the ammonia is complementary to the AB symmetry, and enables a larger degree of dihydrogen bonding between the two.
In terms of implications towards AB regeneration, the sheer favourability in ammonia of the H B Á Á ÁH A dihydrogen bonds may be the reason for the observed B-H bond redistribution in AB waste treated with ammonia. 18Typically, the reduction of Fig. 20 Theoretically optimised AB-ammonia dimer in which bonding to both sides of the AB molecule from one ammonia is observed. 43X 3 (X = O, S, or halide) to BH 3 is the highest energy step in regeneration.In ammonia however, the formation of BH 3 would be stabilised by the numerous dihydrogen bonds formed with the H A s, making the reduction more facile.Ammonia can also stabilise the formation of AB via interacting with the other end of the molecule through a regular hydrogen bond.This is the only interaction possible in other coordinating solvents like THF and other ethers, which lack the multitude of partially positive hydrogens ammonia has and therefore cannot stabilise BH 3 .An example of a different solvent that does is water, but that raises issues in terms of incompatibility with reducing agents like sodium borohydride, that ammonia does not have.This explains the relative success of regeneration processes in ammonia, and suggests that any energetically efficient regeneration would be wise to use it as a solvent.

Experimental
11 B enriched ammonia borane was synthesised following the procedure described in Ramachandran and Gagare. 44 11  enriched sodium borohydride and borodeuteride were sourced from Katchem, ammonium sulfate from Alfa Aesar, and the THF from Sigma-Aldrich.Deuteration of the nitrogen hydrogens was achieved by dissolving in D 2 O and removing the solvent under reduced pressure three times.All the chemicals were used as received without further purification.For the neutron scattering experiments, the THF and ammonia were both sourced from Sigma-Aldrich and again used without further purification.
The neutron scattering experiments were run on the Small Angle Neutron Diffractometer for Amorphous and Liquid Samples (SANDALS) beamline at the spallation neutron source ISIS, at the Rutherford Appleton Laboratory, Didcot, UK. 45 The neutron diffraction patterns of ND 3 BD 3 , NH 3 BH 3 and ND 3 BH 3 were measured dissolved in THF and THF-d 8 , as well as ND 3 BD 3 in a mixture of 50 : 50 THF : THF-d 8 .The ND 3 BD 3 in THF-d 8 sample was 12 wt% AB, with the other samples in the same mole ratio.For the ammonia experiments, less combinations were possible due to exchanging of the ammonia protons with the AB N-Hs.ND 3 BD 3 and ND 3 BH 3 were dissolved and measured in ND 3 along with NH 3 BD 3 in NH 3 .The ND 3 BD 3 in ND 3 sample was 20 wt% AB, with the other samples in the same mole ratio.
The samples in THF were prepared under an inert atmosphere before being sealed in a cell made of a null-scattering titanium/zirconium alloy.The cell had a flat-plate geometry facing the neutron beam to reduce multiple scattering and absorption. 39The diffraction pattern was measured at room temperature, monitored by a thermocouple, and the length of each experiment was approximately 8 hours.The samples in ammonia were measured in the same cell at 220 K, with the ammonia being condensed directly into the cell.The data correction method followed was that implemented by the Gudrun program. 46This involved taking a background measurement of the empty instrument, the empty sample cell, and an incoherently scattering vanadium slab.Gudrun was also used to correct for the inelastic self-scattering using 'top-hat' deconvolution and an iterative method before the analysis by EPSR. 40,47Analysis of the final simulation box after EPSR was performed using the dlputils software. 48

Fig. 2
Fig. 2 Atom labels used throughout the analysis of the EPSR simulations of AB dissolved in THF and liquid ammonia.

Fig. 3
Fig. 3 Diffraction pattern (black) and EPSR fitted structure factor (red) for various isotope combinations of ammonia borane dissolved in THF.

Fig. 4
Fig. 4 Diffraction pattern (black) and EPSR fitted structure factor (red) for various isotope combinations of ammonia borane dissolved in ammonia.

Fig. 5
Fig. 5 Partial radial distribution functions for AB-AB correlations via heteropolar H N Á Á ÁH B and homopolar H B Á Á ÁH B interactions.The red line corresponds to the EPSR simulation of AB dissolved in ammonia and the black line that of AB dissolved in THF.

Fig. 6
Fig.6Cluster size distribution for AB-AB clusters in THF (black) and liquid ammonia (red).Clusters were defined as having an H N Á Á ÁH B distance of less than or equal to 2.7 Å, showing that long-range AB-AB interactions are more common in THF than ammonia.The inset shows the much larger range in cluster size present in THF.

Fig. 7
Fig. 7 Probability distributions of the number of H N s in dihydrogen bonding distance (o2.7 Å) to an H B , in the models of AB dissolved in THF (left) and liquid ammonia (right).These again highlight the greater AB-AB interaction in THF.

Fig. 8
Fig. 8 Spatial density functions for AB-AB interactions via H N Á Á ÁH B dihydrogen bonds, for AB dissolved in liquid ammonia.The density shows the 20% most likely positions for the neighbouring molecule's (left) H B s and (right) H N s in the first coordination shell, defined as having a distance of 2.0-4.0Å between the centre of the three-fold symmetric hydrogens and the neighbouring molecule's heteropolar hydrogen.The blue atom is the AB nitrogen and the green the boron.

Fig. 9
Fig.9Angular radial distribution function for AB-AB pairs in THF, showing the anti-parallel nature of the AB dipoles.The angle y is defined above, along with AB-AB orientations that satisfies the graph as well as the spatial density functions.The distance is that between H N and an H B .The blue atom is the AB nitrogen and green the AB boron.

Fig. 10
Fig. 10 Partial radial distribution functions for H N (black) and H B (red) to O1, calculated from the EPSR fitted system, for AB dissolved in THF.

Fig. 11
Fig. 11 Spatial density functions showing the location of the 20% most likely positions for molecules in the first coordination shell (2.0-3.6 Å) around AB and THF.The red density is the position of the THF O1s and the blue density the AB H N s.The blue AB atom is the AB nitrogen.

Fig. 12
Fig.12Angular radial distribution function for AB-THF interactions, illustrating the hydrogen bonding between H N and O1.y is the angle between the AB and THF dipole moments.The distance is calculated between H N and O1.AB-THF orientations that fit these peaks and the spatial density functions are shown above.

Fig. 13
Fig. 13 Probability distribution of (left) the number of O1s in hydrogen bonding distance (o2.7 Å) to each H N and (right) the number of H N s around O1, for the model of AB dissolved in THF.

Fig. 14
Fig. 14 Partial radial distribution functions for H N (black) and H B (red) to H A as well as to N A , calculated from the EPSR fitted system, for AB dissolved in liquid ammonia.

Fig. 15
Fig. 15 Probability distributions of the number of (left) H A s within dihydrogen bonding range (o2.5 Å) of H B s and (right) the number of N A s in the same distance of H N , for AB dissolved in liquid ammonia.

Fig. 16
Fig. 16 Spatial density functions showing the location of the 20% most likely positions for ammonia or AB molecules in the first coordination shell around the other, for the model solution of AB (20 wt%) in ammonia.The red density shows the H N to N A pair in a range of 2.2-3.5 Å, while the orange density shows the H B to H A pair in a range of 2.2-4.2Å.

Fig. 17
Fig.17AB-ammonia orientations that satisfy the spatial density functions and angular radial distribution functions.

Fig. 18
Fig. 18 Normalised probability distributions for the angle X-HÁ Á ÁYA, where X is the hydrogen bond donor and Y the acceptor, in (top) AB dissolved in THF and (bottom) AB dissolved in liquid ammonia.

Fig. 19
Fig. 19 Normalized probability distributions for the angle A-YÁ Á ÁHX, where Y is the hydrogen bond acceptor and X the donor, in (top) AB dissolved in THF and (bottom) AB dissolved in liquid ammonia.

Table 1
EPSR seed parameters used in the simulations