Ionic liquids (IL) and hydrogen bonding (H-bonding) are two diverse fields for which there is a developing recognition of significant overlap. Doubly ionic H-bonds occur when a H-bond forms between a cation and anion, and are a key feature of ILs. Doubly ionic H-bonds represent a wide area of H-bonding which has yet to be fully recognised, characterised or explored. H-bonds in ILs (both protic and aprotic) are bifurcated and chelating, and unlike many molecular liquids a significant variety of distinct H-bonds are formed between different types and numbers of donor and acceptor sites within a given IL. Traditional more neutral H-bonds can also be formed in functionalised ILs, adding a further level of complexity. Ab initio computed parameters; association energies, partial charges, density descriptors as encompassed by the QTAIM methodology (ρBCP), qualitative molecular orbital theory and NBO analysis provide established and robust mechanisms for understanding and interpreting traditional neutral and ionic H-bonds. In this review the applicability and extension of these parameters to describe and quantify the doubly ionic H-bond has been explored. Estimating the H-bonding energy is difficult because at a fundamental level the H-bond and ionic interaction are coupled. The NBO and QTAIM methodologies, unlike the total energy, are local descriptors and therefore can be used to directly compare neutral, ionic and doubly ionic H-bonds. The charged nature of the ions influences the ionic characteristics of the H-bond and vice versa, in addition the close association of the ions leads to enhanced orbital overlap and covalent contributions. The charge on the ions raises the energy of the Ylp and lowers the energy of the X-H σ* NBOs resulting in greater charge transfer, strengthening the H-bond. Using this range of parameters and comparing doubly ionic H-bonds to more traditional neutral and ionic H-bonds it is clear that doubly ionic H-bonds cover the full range of weak through to very strong H-bonds.

5Deep eutectic solvents (DESs) are exemplars of systems with the ability to form neutral, ionic and doubly ionic H-bonds. Herein, the pairwise interactions of the constituent components of the choline chlorideurea DES are examined. Evidence is found for a tripodal CH⋅⋅⋅Cl doubly ionic H-bond motif. Moreover it is found that the covalency of doubly ionic H-bonds can be greater than, or comparable with, neutral and ionic examples. In contrast to many traditional solvents, an "alphabet soup" of many different types 10 of H-bond (OH⋅⋅⋅O=C, NH⋅⋅⋅O=C, OH⋅⋅⋅Cl, NH⋅⋅⋅Cl, OH⋅⋅⋅NH, CH⋅⋅⋅Cl, CH⋅⋅⋅O=C, NH⋅⋅⋅OH and NH⋅⋅⋅NH) can form. These H-bonds exhibit substantial flexibility in terms of number and strength. It is anticipated that H-bonding will have a significant impact on the entropy of the system and thus could play an important role in the formation of the eutectic.

In this paper the structural and energetic landscapes of ion-pair dimer conformers of 1,3-dimethylimidazolium based ionic liquids have been explored ([C1C1im][A])2, A = Cl(-), [NO3](-), [MeSO4](-), [OTf](-) and [BF4](-)). A common low-energy conformer has been selected for full electronic structure analysis. We have compared and contrasted each cluster based on the relative hydrogen bonding ability (β-value) of the anion, which varies experimentally as Cl(-) > [NO3](-) ≈ [MeSO4](-) > [OTf](-) ≈ [BF4](-). Correlations between experimental β-values, computed binding energies, charge transfer and various hydrogen bonding data have been made and outliers have been explained in terms of environmental effects present in the liquid phase. This is most evident in the structurally similar [MeSO4](-) and [OTf](-) anions that have very similar hydrogen bonding motifs, but significantly different β-values. Moreover, detailed analysis of the cluster molecular orbitals, for each anion, reveals a subtle interplay between two modes of interaction, an in-plane traditional H-bonding and inter-planar anion-π interaction. Inter-planar anion-π interactions are particularly prominent for the [NO3](-) cluster. We have rationalized how the full range of interactions could impact on the structuring of ILs at surfaces and the effect these may have on viscosity.

A range of methods for the computational prediction of experimentally derived α and β Kamlet-Taft parameters, indicators of hydrogen bond (H-bond) acidity and basicity for ionic liquids (ILs) have been explored. Most usefully, a good correlation has been established between several simple and easily computed quantities which allow for a "quick bench-top" evaluation. More accurate, but also more sophisticated methods employing TD-DFT calculations involving the Kamlet-Taft dyes have been examined and evaluated. Importantly, these techniques open up the opportunity for pre-screening and a priori prediction of properties for ILs not yet synthesised. A key fundamental insight into IL H-bonds has been the determination of an estimate for the energy associated with replacing both neutral molecules in a H-bond with ionic molecules, thus forming the "doubly ionic" H-bond found in ILs.

Experimental near edge X-ray absorption fine structure (NEXAFS) spectra, X-ray photoelectron (XP) spectra and Auger electron spectra are reported for sulfur in ionic liquids (ILs) with a range of chemical structures. These values provide experimental measures of the atomic charge in each IL and enable the evaluation of the suitability of NEXAFS spectroscopy and XPS for probing the relative atomic charge of sulfur. In addition, we use Auger electron spectroscopy to show that when XPS binding energies differ by less than 0.5 eV, conclusions on atomic charge should be treated with caution. Our experimental data provides a benchmark for calculations of the atomic charge of sulfur obtained using different methods. Atomic charges were computed for lone ions and ion pairs, both in the gas phase (GP) and in a solvation model (SMD), with a wide range of ion pair conformers considered. Three methods were used to compute the atomic charges: charges from the electrostatic potential using a grid based method (ChelpG), natural bond orbital (NBO) population analysis and Bader's atoms in molecules (AIM) approach. By comparing the experimental and calculated measures of the atomic charge of sulfur, we provide an order for the sulfur atoms, ranging from the most negative to the most positive atomic charge. Furthermore, we show that both ChelpG and NBO are reasonable methods for calculating the atomic charge of sulfur in ILs, based on the agreement with both the XPS and NEXAFS spectroscopy results. However, the atomic charges of sulfur derived from ChelpG are found to display significant, non-physical conformational dependence. Only small differences in individual atomic charge of sulfur were observed between lone ion (GP) and ion pair IL(SMD) model systems, indicating that ion-ion interactions do not strongly influence individual atomic charges.