Combined Spectroscopic, Thermodynamic, and Theoretical Approach for Detecting and Quantifying Hydrogen Bonding and Dispersion Interaction in Ionic Liquids

ConspectusIonic liquids (ILs) are attracting increasing interest in science and engineering due to their unique properties that can be tailored for specific applications. Clearly, a better understanding of their behavior on the microscopic scale will help to elucidate macroscopic fluid phenomena and thereby promote potential applications. The advantageous properties of these innovative fluids arise from the delicate balance of Coulomb interactions, hydrogen bonding, and dispersion forces. The development of these properties requires a fundamental understanding of the strength, location, and direction of the different types of interactions and their contribution to the overall phase behavior. Contrary to expectations, hydrogen bonding and dispersion interactions have a significant influence on the structure, dynamics, and phase behavior of ILs.The synergy between experimental and theoretical methods has now advanced to a stage where hydrogen bonds and dispersion effects as well as the competition between the two can be studied in detail. In this account, we demonstrate that a suitable combination of spectroscopic, thermodynamic, and theoretical methods enables the detection, dissection, and quantification of noncovalent interactions, even in complex systems such as ionic liquids. This approach encompasses far-infrared vibrational spectroscopy (FIR), various thermodynamic methods for determining enthalpies of vaporization, and quantum chemical techniques that allow us to switch dispersion contributions on or off when calculating the energies and spectroscopic properties of clusters.We briefly discuss these experimental and theoretical methods, before providing various examples illustrating how the mélange of Coulomb interaction, hydrogen bonds, and dispersion forces can be analyzed, and their individual contributions quantified. First, we demonstrated that both hydrogen bonding and dispersion interactions are manifested in the FIR spectra and can be quantified by observed shifts of characteristic spectral signatures. Through the selection of suitable protic ionic liquids (PILs) featuring anions with varying interaction strengths and alkyl chain lengths, we were able to demonstrate that dispersion interactions can compete with hydrogen bonding. The resultant transition enthalpy serves as a measure of the dispersion interaction. Contrary to expectations, PILs possess lower enthalpies of vaporization compared with aprotic ILs (AILs). The reason for this is simple: In protic ILs, ion pairs carry both the hydrogen bond and attractive dispersion between the cation and anion into the gas phase. By utilizing a well-curated set of protic ILs and molecular analogues, we successfully disentangled Coulomb interaction, hydrogen bonding, and dispersion interaction through purely thermodynamic methods.

Long-Lived Weak Ion Pairs in Ionic Liquids: An Insight from All-Atom Molecular Dynamics Simulations

The microstructure and local dynamics of ions in room-temperature ionic liquids (RTILs) have drawn a lot of attention due to their extensive potential applications in numerous fields. It is well-known that the widely used definitions of ion pairs (IPs) cannot reflect the full picture of RTILs. In this study, we find a universal residence time (τ_MR), which is regardless of the number of counterions in the first solvation shell in RTILs. Inspired by this, we propose a weak IP (WIP) model from a spatiotemporal perspective and demonstrate that the WIPs are long-lived and that their lifetimes obey a log-normal distribution, which is different from the literature. In addition, the electrostatic interactions are the main factors in the formation of WIPs, and the reorientations of ions are vital to the ruptures of WIPs. This research provides a new perspective for understanding the microstructural and dynamical properties of RTILs.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

Can a gas phase contact ion pair containing a hydrocarbon carbocation be formed in the ground state?

So far, no conclusive evidence of a ground-state contact ion-pair containing a hydrocarbon carbocation has been given in the gas phase.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Heterogeneity in the microstructure and dynamics of tetraalkylammonium hydroxide ionic liquids: insight from classical molecular dynamics simulations and Voronoi tessellation analysis

Microscopic structural and dynamic heterogeneities were investigated for three ionic liquids (ILs), tetraethylammonium hydroxide, tetrapropylammonium hydroxide, and tetrabutylammonium hydroxide employing classical molecular dynamics (MD) simulations.

Fine probing the effect of replacing [PF6]- with [PF3(C2F5)3]- on the local structure and nanoscale organization of [bmim]+-based ionic liquids using MD simulation

Comparative all-atom molecular dynamics simulations are used to study the microscopic local structure and interionic interactions of two ionic liquids (ILs) composed of the 1-butyl-3-methylimidazolium cation, [bmim]+, coupled with the hexafluorophosphate, [PF6]-, or tris(pentafluoroethyl)trifluorophosphate, [FAP]-, anions. Respective distribution functions clearly reveal that the structural correlations between the cation and anion decrease when (i) replacing [PF6]- with [FAP]-, (ii) scaling the partial atomic charges, and (iii) considering the anion's structural flexibility versus rigidity. Replacement of [PF6]- with [FAP]- expands the nonpolar domains totally and causes the decreasing of the three-dimensional polar networks as well as the diminishing of the nano-aggregation of cation side chains. Current simulations show that with increasing the anion size and its charge delocalization, the probability of the in-plane cation-anion conformation, its related hydrogen bond acceptor ability, and the cation-cation π-π interaction decreases in accordance with the fluidity enhancements of the corresponding imidazolium-based IL. Hence, structural findings can explain and justify rationally the origins of the observed trends in the simulated dynamical properties of these ILs in our previous report. A complete understanding of the microscopic structure of ILs is necessary to control the outstanding properties of ILs as designer solvents that will support experimentalists for the best engineering design and a breakthrough efficiency of IL-related processes.

Prediction of 1H NMR chemical shifts for clusters of imidazolium-based ionic liquids

Nuclear magnetic resonance (NMR) has been widely used to elucidate the bulk structure of ionic liquids. In this work, we calculated ¹H NMR chemical shifts of 1-ethyl-3-methylimidazolium (C₂mim⁺) ionic liquids combined with various anions such as chloride (Cl), tetrafluoroborate (BF₄), hexafluorophosphate (PF₆), acetate (OAc), trifluoroacetate (TFA), and dicyanamide (DCA). The previously established level of theory, HF/6-311G+(3df,2p), was used for the accurate prediction of NMR chemical shifts both in gas phase and in solvents with varying dielectric constant such as CHCl₃ and ethanol. The following factors affecting the predicted proton chemical shifts were considered. Firstly, ionic clusters consisting of 2, 8 and 16 ion pairs were optimized to model interionic interactions present in the bulk of ionic liquids. In larger clusters the distribution of the calculated chemical shifts of individual protons in the C₂mim⁺ cation was examined with respect to the position of the cation in the cluster. We further confirmed that electronic properties of ionic liquids such as magnetic shielding had local nature, thus allowing us to accurately predict proton NMR chemical shifts of ionic liquids from relatively small-sized clusters. Secondly, solvent effects in single ion pairs as well as larger ionic clusters were accounted through a Conductor-like Polarisable Continuum Model (CPCM). Solvent effects generated through a dielectric constant of either chloroform or ethanol were found to be important in single ion pairs due to improved description of interionic distances. With increasing cluster size the difference between gas-phase and CPCM optimized structures became minimal, thus resulting in similar values for calculated ¹H NMR chemical shifts. We also established that the model size that produced the best results for imidazolium ionic liquids strongly depended on the anion type. Strongly coordinating anions such as chloride and acetate require calculations of clusters consisting of at least 8 ion pairs, whereas weakly coordinating anions produce excellent accuracy for single ion pairs optimized in the presence of solvent. The polarity of the solvent was found to play a minor role.

Mapping the Free Energy of Lithium Solvation in the Protic Ionic Liquid Ethylammonuim Nitrate: A Metadynamics Study

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 eV. Other less probable conformations having different solvation pattern are also investigated.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Trends in Two- and Three-Body Effects in Multiscale Clusters of Ionic Liquids

Applications of higher correlated levels of ab initio theory to condensed systems require a significant amount of computational resources. The recent development of the fragment molecular orbital (FMO) approach alleviates this issue by splitting the system into individual fragments and achieves the accuracy of the method by accounting for all possible two-body and three-body interactions. In this work a comprehensive application of the FMO approach in combination with a second order of Møller-Plesset perturbation theory method, MP2, is presented for multiscale clusters of ionic liquids such as [C₁mim]X, [C₁mpyr]X, [C₂py]X, and [NMe₄]X, where X = chloride and tetrafluoroborates, BF₄^-, with the clusters varying in size from 4, 8, 16, to 32 ion pairs. Reliable cutoff criteria for the inclusion of two-body and three-body interactions are identified for both HF energy and MP2 correlation energy to achieve the desired accuracy of 1 kJ mol^-1. The importance of two-body and three-body interactions in ionic liquids is also discussed.

Accurate prediction of energetic properties of ionic liquid clusters using a fragment-based quantum mechanical method

Accurate prediction of physicochemical properties of ionic liquids (ILs) is of great significance to understand and design novel ILs with unique properties.

Ion pairing in ionic liquids

In the present article we briefly review the extensive discussion in literature about the presence or absence of ion pair-like aggregates in ionic liquids. While some experimental studies point towards the presence of neutral subunits in ionic liquids, many other experiments cannot confirm or even contradict their existence. Ion pairs can be detected directly in the gas phase, but no direct method is available to observe such association behavior in the liquid, and the corresponding indirect experimental proofs are based on such assumptions as unity charges at the ions. However, we have shown by calculating ionic liquid clusters of different sizes that assuming unity charges for ILs is erroneous, because a substantial charge transfer is taking place between the ionic liquid ions that reduce their total charge. Considering these effects might establish a bridge between the contradicting experimental results on this matter. Beside these results, according to molecular dynamics simulations the lifetimes of ion-ion contacts and their joint motions are far too short to verify the existence of neutral units in these materials.

Infrared and Raman Spectroscopy from Ab Initio Molecular Dynamics and Static Normal Mode Analysis: The C-H Region of DMSO as a Case Study

Carbon-hydrogen (C-H) vibration modes serve as key probes in the chemical identification of hydrocarbons and in vibrational sum-frequency generation spectroscopy of hydrocarbons at the liquid/gas interface. Their assignments pose a challenge from a theoretical viewpoint. In this work, we present a detailed study of the C-H stretching region of dimethyl sulfoxide using a new ab initio molecular dynamics (AIMD) module that we have implemented in NWChem. Through a combination of AIMD simulations and static normal mode analysis, we interpret experimental infrared and Raman spectra and explore the role of anharmonic effects in this system. Comprehensive anharmonic normal mode analysis of the C-H stretching region casts doubt upon previous experimental assignments of the shoulder on the symmetric C-H stretching peak. In addition, our AIMD simulations also show significant broadening of the in-phase symmetric C-H stretching resonance, which suggests that the experimentally observed shoulder is due to thermal broadening of the symmetric stretching resonance.

Hydrogen bonding in ionic liquids

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.

The effect of dispersion forces on the interaction energies and far infrared spectra of protic ionic liquids

We could show by means of dispersion-corrected DFT calculations that the interaction energy in protic ionic liquids can be dissected into Coulomb interaction, hydrogen bonding and dispersion interaction. The H-bond energy as well as the dispersion energy can be quantified to be 50 kJ mol(-1) each representing ten percent of the overall interaction energy. The dispersion interaction could be dissected into two portions. One third could be related to the dispersion interaction within an ion-pair enhancing the H-bond strength, two thirds stem from dispersion interaction between the ion-pairs. This distribution of dispersion interaction is reflected in the far infrared (FIR) spectra. The H-bond band is shifted weaker than the low frequency band where the latter indicates diffuse cation-anion interaction and H-bond bending motions. Finally, we can dissect the different types of interaction energies indicating their characteristic influence on vibrational modes in the FIR.

Simulations of room temperature ionic liquids: from polarizable to coarse-grained force fields

This perspective article summarizes the recent advances in the classical molecular modelling of room temperature ionic liquids.

The effect of dispersion interactions on the properties of LiF in condensed phases

Classical molecular dynamics simulations are performed on LiF in the framework of the polarizable ion model. The overlap repulsion and polarization terms of the interaction potential are derived on a purely non-empirical, first-principles basis. For the dispersion, three cases are considered: a first one in which the dispersion parameters are set to zero and two others in which they are included, with different parametrizations. Various thermodynamic, structural and dynamic properties are calculated for the solid and liquid phases. The melting temperature is also obtained from direct coexistence simulations of the liquid and solid phases. Dispersion interactions appear to have an important effect on the densities of both phases and on the melting point, although the liquid properties are not affected when simulations are performed in the NVT ensemble at the experimental density.

Understanding the evaporation of ionic liquids using the example of 1-ethyl-3-methylimidazolium ethylsulfate

In this work we present a comprehensive temperature-dependence analysis of both the structural and the dynamic properties of a vaporized ionic liquid (1-ethyl-3-methylimidazolium ethylsulfate). This particular ionic liquid is known to be distillable from experimental studies and thus enables us to deepen the understanding of the evaporation mechanism of ionic liquids. We have used ab initio molecular dynamics of one ion pair at three different temperatures to accurately describe the interactions present in this model ionic liquid. By means of radial and spatial distribution functions a large impact on the coordination pattern at 400 K is shown which could explain the transfer of one ion pair from the bulk to the gas phase. Comparison of the free energy surfaces at 300 K and 600 K supports the idea of bulk phase-like and gas phase-like ion pairs. The different coordination patterns caused by the temperature, describing a loosening of the anion side chains, are also well reflected in the power spectra. The lifetime analysis of typical conformations for ionic liquids shows a characteristic behavior at 400 K (temperature close to the experimental evaporation temperature), indicating that conformational changes occur when the ionic liquid is evaporated.

Carbene formation in ionic liquids: spontaneous, induced, or prohibited?

We present a theoretical study of carbene formation from the 1-ethyl-3-methylimidazolium acetate ionic liquid in the absence and presence of CO2 in gas and liquid phase. Although CO2 physisorption constitutes a precursory step of chemisorption (the CO2's reaction with carbenes, which forms from cations via proton abstraction by anions), it also enables a very stable CO2-anion associate. However, this counteracts the chemical absorption by reducing the basicity of the anion and the electrophilicity of the CO2, which is reflected by charge transfer. Accordingly, the observable carbene formation in the gas phase is hindered in the presence of CO2. In the neat liquid, the carbene formation is also suppressed by the charge screening compared to the case of the gas phase; nevertheless, indications for carbene incidents appear. Interestingly, in the CO2-containing liquid we detect more carbene-like incidents than in the neat one, which is caused by the way CO2 is solvated. Despite the weakness of the CO2-cation interaction, the CO2-anion associate is distorted by cations, which can be seen in longer associate distances and reduced "binding" energies. While the single solvating anion is shifted away from CO2, many more solvating cations approach it compared to the case of the gas phase. This leads to the conclusion that while the ionic liquid effect stabilizes charged species, introducing neutral species such as CO2 provides an opposite trend, leading to an inverse ionic liquid effect with the facilitation of carbene formation and thus of chemical absorption.