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 we have explored the structural and energetic landscape of potential π(+)-π(+) stacked motifs, hydrogen-bonding arrangements and anion-π(+) interactions for gas-phase ion pair (IP) conformers and IP-dimers of 1,3-dimethylimidazolium chloride, [C1C1im]Cl. We classify cation-cation ring stacking as an electron deficient π(+)-π(+) interaction, and a competitive anion on-top IP motif as an anion-donor π(+)-acceptor interaction. 21 stable IP-dimers have been obtained within an energy range of 0-126 kJ mol(-1). The structures have been found to exhibit a complex interplay of structural features. We have found that low energy IP-dimers are not necessarily formed from the lowest energy IP conformers. The sampled range of IP-dimers exhibits new structural forms that cannot be recovered by examining the ion-pairs alone, moreover the IP-dimers are recovering additional key features of the local liquid structure. Including dispersion is shown to impact both the relative energy ordering and the geometry of the IPs and IP-dimers, however the impact is found to be subtle and dependent on the underlying functional.

A detailed investigation of hydrogen bonding in the pure ionic liquids [C4C1im]Cl and [C2C1im]Cl has been carried out using primarily molecular dynamics techniques. Analyses of the individual atom-atom pair radial distribution functions, and in particular those for C···Cl(-), have revealed that hydrogen bonding to the first methylene or methyl units of the substituent groups is important. Multiple geometric criteria for defining a hydrogen bond have been applied, and in particular the choice of the cutoff angle has been carefully examined. The interpretation of hydrogen bonding within these ionic liquids is highly angle dependent, and justification is provided for why it may be appropriate to employ a wider angle criteria than the 30° used for water or alcohol systems. The different types of hydrogen bond formed are characterized, and "top" conformations where the Cl anion resides above (or below) the imidazolium ring are investigated. The number of hydrogen bonds undertaken by each hydrogen atom (and the chloride anion) is quantified, and the propensity to form zero, one, or two hydrogen bonds is established. The effects of an increase in temperature on the static hydrogen bonding are also briefly examined.

A systematic electronic structure analysis of hydrogen bonding (H-bonding), anion-π(+) and π(+)-π(+) interactions present in [C1C1im]Cl ion-pairs (IPs) and selected [C1C1im]2Cl2 IP-dimers has been carried out. Interactions have been characterised using a combination of QTAIM, NCIPLOT, NBO and qualitative MO theory. IP-dimers form non-directional charge quadrupolar arrangements due to Coulombic interactions. These are found to associate either as clusters or as loosely associated IP-IP structures. Large conformational changes are found to occur for very little cost in energy, indicating that charge screening is essentially independent of the cation ring orientation. H-bond formation is accompanied by charge transfer and polarisation of the entire [C1C1im](+) ring. Charge transfer does not follow the same trend for the CHelpG, QTAIM and NBO methods. Weak "stacked" π(+)-π(+) interactions are stabilised in the presence of anions, which locate between and at the periphery of the rings, novel strongly bent H-bonds are also present. Primary (ring; C-H···Cl(-)) H-bonds and anion-π(+) (C(2)···Cl(-)) interactions are found to decay more rapidly with distance than secondary (aliphatic; C(M)-H···Cl(-)) H-bonds. This leads to an increase in the relative importance of secondary H-bond interactions in the IP-dimers. Moreover, rotation of the methyl groups within the "stacked" π(+)-π(+) IP-dimers facilitates the formation of (stronger) linear secondary H-bonds. Thus, compared to isolated IPs, secondary H-bonds may play an increased role within the condensed phase. Overall we find that structural fluidity is facilitated by fluctuating hydrogen bond, π(+)-π(+) and anion-π(+) interactions.

The structures of mixtures of ionic liquids (ILs) featuring a common 1-butyl-3-methylimidazolium ([C4C1im](+)) cation but different anions have been investigated both experimentally and computationally. (1)H and (13)C NMR of the ILs and their mixtures has been performed both on the undiluted liquids and those diluted by CD2Cl2. These experiments have been complemented by quantum chemical density functional theory calculations and molecular dynamics simulations. These techniques have identified the formation of preferential interactions between H(2) of the imidazolium cation and the most strongly hydrogen bond (H-bond) accepting anion. In addition, a preference for the more weakly H-bond accepting anion to interact above the imidazolium ring through anion-π(+) interactions has been identified. The modelling of these data has identified that the magnitude of these preferences are small, of the order of only a few kJ mol(-1), for all IL mixtures. No clustering of the anions around a specific cation could be observed, indicating that these interactions arise from the reorientation of the cation within a randomly assigned network of anions. π(+)-π(+) stacking of the imidazolium cations was also studied and found to be promoted by ILs with a strong H-bond accepting anion. Stacking interactions are easily disrupted by the introduction of small proportions (<50 mol%) of a weakly coordinating anion due to their propensity to form anion-π(+) interactions. These results suggest that the formation of IL mixtures with different anions leads to subtle structural changes of much lower energy than the Coulombic ordering of ions, accounting for why most IL mixtures exhibit ideal, or nearly ideal, behaviour.

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.

SAXS, 129Xe NMR and PALS were used to interrogate the relationship between the structure, free volume and physicochemical properties of ionic liquid mixtures.