Comparison of the local binding motifs in the imidazolium-based ionic liquids [EMIM][BF4] and [EMMIM][BF4] through cryogenic ion vibrational predissociation spectroscopy: Unraveling the roles of anharmonicity and intermolecular interactions

What quantum chemical simulations tell us about the infrared spectra, structure and interionic interactions of a bulk ionic liquid

The recently developed efficient protocol to explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is applied to ionic liquid 1-ethyl-3-methyl-imidazolium tetrafluoroborate [Emim][BF₄], and its C2-deuterated analog [Emim-d][BF₄]. It is shown that the solvation strongly modifies the frequencies and IR intensities of both cationic and anionic components of the ionic liquids. The main features of the bulk spectra are reproduced by the simulations for cluster ([Emim][BF₄])₈, representing an ion pair solvated by the first solvation shell. The geometry of the cluster closely resembles the solid-state structure of the actual ionic liquid and is characterized by short contacts of all CH moieties of the imidazolium ring with [BF₄]^- anions. Both structural and spectroscopic analyses allow the contacts to be interpreted as hydrogen bonds of approximately equal strength. The enthalpies of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to 1.2-1.5 kcal mol^-1, while the analogous estimates obtained for the gas-phase charged species [Emim][BF₄]₂^- and [Emim]₂[BF₄]⁺ increase to 3.6-3.9 kcal mol^-1.

1-Alkyl-3-methylimidazolium cation binding preferences in hexafluorophosphate ionic liquid clusters determined using competitive TCID measurements and theoretical calculations

Ionic liquids (ILs) exhibit unique properties that have led to their development and widespread use for a variety of applications. Development efforts have generally focused on achieving desired macroscopic properties via tuning of the IL through variation of the cations and anions. Both the macroscopic and microscopic properties of an IL influence its tunability and thus feasibility of use for selected applications. Works geared toward a microscopic understanding of the nature and strength of the intrinsic cation-anion interactions of ILs have been limited to date. Specifically, the intrinsic strength of the cation-anion interactions in ILs is largely unknown. In previous work, we employed threshold collision-induced dissociation (TCID) approaches supported and enhanced by electronic structure calculations to determine the bond dissociation energies (BDEs) and characterize the nature of the cation-anion interactions in a series of four 2 : 1 clusters of 1-alkyl-3-methylimidazolium cations with the hexafluorophosphate anion, [2C_nmim:PF₆]⁺. To examine the effects of the 1-alkyl chain on the structure and energetics of binding, the cation was varied over the series: 1-ethyl-3-methylimidazolium, [C₂mim]⁺, 1-butyl-3-methylimidazolium, [C₄mim]⁺, 1-hexyl-3-methylimidazolium, [C₆mim]⁺, and 1-octyl-3-methylimidazolium, [C₈mim]⁺. The variation in the strength of binding among these [2C_nmim:PF₆]⁺ clusters was found to be similar in magnitude to the average experimental uncertainty in the measurements. To definitively establish an absolute order of binding among these [2C_nmim:PF₆]⁺ clusters, we extend this work again using TCID and electronic structure theory approaches to include competitive binding studies of three mixed 2 : 1 clusters of 1-alkyl-3-methylimidazolium cations and the hexafluorophosphate anion, [C_n-2mim:PF₆:C_nmim]⁺ for n = 4, 6, and 8. The absolute BDEs of these mixed [C_n-2mim:PF₆:C_nmim]⁺ clusters as well as the absolute difference in the strength of the intrinsic binding interactions as a function of the cation are determined with significantly improved precision. By combining the thermochemical results of the previous independent and present competitive measurements, the BDEs of the [2C_nmim:PF₆]⁺ clusters are both more accurately and more precisely determined. Comparisons are made to results for the analogous [2C_nmim:BF₄]⁺ and [C_n-2mim:BF₄:C_nmim]⁺ clusters previously examined to elucidate the effects of the [PF₆]^- and [BF₄]^- anions on the binding.

Nature and strength of intrinsic cation-anion interactions of 1-alkyl-3-methylimidazolium hexafluorophosphate clusters

Imidazolium-based cations and the hexafluorophosphate anion are among the most commonly used ionic liquids (ILs). Yet, the nature and strength of the intrinsic cation-anion interactions, and how they influence the macroscopic properties of these ILs are still not well understood. Threshold collision-induced dissociation is utilized to determine the bond dissociation energies (BDEs) of the 2 : 1 clusters of 1-alkyl-3-methylimidazolium cations and the hexafluorophosphate anion, [2Cnmim:PF6]+. The cation, [Cnmim]+, is varied across the series, 1-ethyl-3-methylimidazolium [C2mim]+, 1-butyl-3-methylimidazolium [C4mim]+, 1-hexyl-3-methylimidazolium [C6mim]+, 1-octyl-3-methylimidazolium [C8mim]+, to examine the structural and energetic effects of the size of the 1-alkyl substituent of the cation on the binding to [PF6]-. Complementary electronic structure methods are employed for the [Cnmim]+ cations, (Cnmim:PF6) ion pairs, and [2Cnmim:PF6]+ clusters to elucidate details of the cation-anion interactions and their impact on structure and energetics. Multiple levels of theory are benchmarked with the measured BDEs including B3LYP, B3LYP-GD3BJ, and M06-2X each with the 6-311+G(d,p) basis set for geometry optimizations and frequency analyses and the 6-311+G(2d,2p) basis set for energetic determinations. The modest structural variation among the [Cnmim]+ cations produces only minor structural changes and variation in the measured BDEs of the [2Cnmim:PF6]+ clusters. Present results are compared to those previously reported for the analogous 1-alkyl-3-methylimidazolium tetrafluoroborate IL clusters to compare the effects of these anions on the nature and strength of the intrinsic binding interactions.

Ionic Liquid Clusters Generated from Electrospray Thrusters: Cold Ion Spectroscopic Signatures of Size-Dependent Acid-Base Interactions

We determine the intramolecular distortions at play in the 2-hydroxyethylhydrazinium nitrate (HEHN) ionic liquid (IL) propellant, which presents the interesting case that the HEH⁺ cation has multiple sites (i.e., hydroxy, primary amine, and secondary ammonium groups) available for H-bonding with the nitrate anion. These interactions are quantified by analyzing the vibrational band patterns displayed by cold cationic clusters, (HEH⁺)_n(NO₃^-)_n-1, n = 2-6, which are obtained using IR photodissociation of the cryogenically cooled, mass-selected ions. The strong interaction involving partial proton transfer of the acidic N-H proton in HEH⁺ cation to the nitrate anion is strongly enhanced in the ternary n = 2 cluster but is suppressed with increasing cluster size. The cluster spectra recover the bands displayed by the bulk liquid by n = 5, thus establishing the minimum domain required to capture this aspect of macroscopic behavior.

Gas-Phase Binding Energies and Dissociation Dynamics of 1-Alkyl-3-Methylimidazolium Tetrafluoroborate Ionic Liquid Clusters

Ionic liquids (ILs) have become increasingly popular due to their useful and unique properties, yet there are still many unanswered questions regarding their fundamental interactions. In particular, details regarding the nature and strength of the intrinsic cation-anion interactions and how they influence the macroscopic properties of ILs are still largely unknown. Elucidating the molecular-level details of these interactions is essential to the development of better models for describing ILs and enabling the purposeful design of ILs with properties tailored for specific applications. Current uses of ILs are widespread and diverse and include applications for energy storage, electrochemistry, designer/green solvents, separations, and space propulsion. To advance the understanding of the energetics, conformations, and dynamics of gas-phase IL clustering relevant to space propulsion, threshold collision-induced dissociation approaches are used to measure the bond dissociation energies (BDEs) of the 2:1 clusters of 1-alkyl-3-methylimidazolium cations and tetrafluoroborate, [2C_nmim:BF₄]⁺. The cation, [C_nmim]⁺, is varied across the series, 1-ethyl-3-methylimidazolium [C₂mim]⁺, 1-butyl-3-methylimidazolium [C₄mim]⁺, 1-hexyl-3-methylimidazolium [C₆mim]⁺, and 1-octyl-3-methylimidazolium [C₈mim]⁺, to examine the structural and energetic effects of the size of the 1-alkyl substituent on binding. Complementary electronic structure calculations are performed to determine the structures and energetics of the [C_nmim]⁺ and [BF₄]^- ions and their binding preferences in the (C_nmim:BF₄) ion pairs and [2C_nmim:BF₄]⁺ clusters. Several levels of theory, B3LYP, B3LYP-GD3BJ, and M06-2X, using the 6-311+G(d,p) basis set for geometry optimizations and frequency analyses and the 6-311+G(2d,2p) basis set for energetics, are benchmarked to examine their abilities to properly describe the nature of the binding interactions and to reproduce the measured BDEs. The modest structural variation among these [C_nmim]⁺ cations produces only minor structural changes and variation in the measured BDEs of the [2C_nmim:BF₄]⁺ clusters. Present findings indicate that the dominant cation-anion interactions involve the 3-methylimidazolium moieties and that these clusters are sufficiently small that differences in packing effects associated with the variable length of the 1-alkyl substituents are not yet significant.

Photoinduced Intermolecular Electron Transfer in Gas Phase Ion Pairs of the 1-Ethyl-3-methylimidazolium Cation and the Bis(trifluoromethylsulfonyl)imide Anion

In this study, the UV photodissociation of gas phase ion pairs of the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [emim]⁺[tf2n]^-, is shown to proceed primarily through radical intermediates. [emim]⁺[tf2n]^- ion pairs have been shown previously to undergo two-photon-dependent dissociation, but the mechanisms of this have not been probed in detail. By employing a two-laser pump probe spectroscopy and time-dependent density functional theory (TD-DFT) calculations, we have illustrated that one of the major UV photodissociation pathways in [emim]⁺[tf2n]^- ion pairs is an intermolecular electron transfer wherein the anion transfers an electron to the cation resulting in two neutral open-shelled products. These products were observable for at least 1.6 μs post photodissociation, the experimental limit, via detection of the [emim]⁺ cation. This data demonstrates that the likely photoproducts of [emim]⁺[tf2n]^- UV photodissociation are two neutral species that separate spatially, demonstrated through lack of observed relaxation pathways such as electron recombination. TD-DFT and frontier molecular orbital analysis calculations at the MN15/6-311++G(d,p) level are employed to aid in identifying excited state characteristics and support the interpretations of the experimental data. The energetic onset of the intermolecular electron transfer is consistent with previously observed [emim]⁺[tf2n]^- absorption spectra in the bulk and gas phases. The similarities between bulk and gas phase UV spectra imply that this electron-transfer pathway may be a major photodissociation channel in both phases.

Ion Dissociation in Ionic Liquids and Ionic Liquid Solutions

The extent to which cations and anions in ionic liquids (ILs) and ionic liquid solutions are dissociated is of both fundamental scientific interest and practical importance because ion dissociation has been shown to impact viscosity, density, surface tension, volatility, solubility, chemical reactivity, and many other important chemical and physical properties. When mixed with solvents, ionic liquids provide the unique opportunity to investigate ion dissociation from infinite dilution in the solvent to a completely solvent-free state, even at ambient conditions. The most common way to estimate ion dissociation in ILs and IL solutions is by comparing the molar conductivity determined from ionic conductivity measurements such as electrochemical impedance spectroscopy (EIS) (which measure the movement of only the charged, i.e., dissociated, ions) with the molar conductivity calculated from ion diffusivities measured by pulse field gradient nuclear magnetic resonance spectroscopy (PFG-NMR, which gives movement of all of the ions). Because the NMR measurements are time-consuming, the number of ILs and IL solutions investigated by this method is relatively limited. We have shown that use of the Stokes-Einstein equation with estimates of the effective ion Stokes radii allows ion dissociation to be calculated from easily measured density, viscosity, and ionic conductivity data (ρ, η, λ), which is readily available in the literature for a much larger number of pure ILs and IL solutions. Therefore, in this review, we present values of ion dissociation for ILs and IL solutions (aqueous and nonaqueous) determined by both the traditional molar conductivity/PFG-NMR method and the ρ, η, λ method. We explore the effect of cation and anion alkyl chain length, structure, and interaction motifs of the cation and anion, temperature, and the strength of the solvent in IL solutions.

Structure and Stability of the Ionic Liquid Clusters [EMIM]n[BF4]n+1- (n = 1-9): Implications for Electrochemical Separations

Precise functionalization of electrodes with size-selected ionic liquid (IL) clusters may improve the application of ILs in electrochemical separations. Herein we report our combined experimental and theoretical investigation of the IL clusters 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIM]_n[BF₄]_n+1^- (n = 1-9) and demonstrate their selectivity and efficiency toward targeted adsorption of ions from solution. The structures and energies of the IL clusters, predicted with global optimization, agree with and help interpret the ion abundances and stabilities measured by high-mass-resolution electrospray ionization mass spectrometry and collision-induced dissociation experiments. The [EMIM][BF₄]₂^- cluster, which was identified as the most stable IL cluster, was selectively soft-landed onto a working electrode. Electrochemical impedance spectroscopy revealed a lower charge transfer resistance on the soft-landed electrode containing [EMIM][BF₄]₂^- compared with an electrode prepared by drop-casting of an IL solution containing the full range of IL clusters. Our findings indicate that specific IL clusters may be used to increase the efficiency of electrochemical separations by lowering the overpotentials involved.

Chain Length Dependence of Hydrogen Bond Linkages between Cationic Constituents in Hydroxy-Functionalized Ionic Liquids: Tracking Bulk Behavior to the Molecular Level with Cold Cluster Ion Spectroscopy

Hydroxy functionalization of cations in ionic liquids (ILs) can lead to formation of contacts between their OH groups [so-called (c-c) interactions]. One class of these linkages involves cooperatively enhanced hydrogen bonds to anionic partners that are sufficiently strong to overcome the repulsion between two positively charged centers. Herein, we clarify how the propensity for the formation of (c-c) contacts depends on the alkyl chain length between two cationic rings and their OH groups by analyzing the temperature-dependent IR spectra of bulk ILs as well as the vibrational predissociation spectra of ∼35 K complexes comprised of two cations and one anion. This study compares the behavior of two cationic derivatives with ethyl and propyl chains complexed with two different anions: bis(trifluoromethylsulfonyl)imide and tetrafluoroborate. Only the bulk ILs with the longer chain propyl derivative [HPMPip⁺ = 1-(3-hydroxypropyl)-1-methylpiperidinium] display (c-c) interactions. Molecular-level aspects of this docking arrangement are revealed by analyzing the OH stretching fundamentals displayed by the ternary complexes.

Cooperatively enhanced hydrogen bonds in ionic liquids: closing the loop with molecular mimics of hydroxy-functionalized cations

The combined experimental and theoretical approach for the gas and the liquid phases provides a quantitative understanding of the competition between differently H-bonded and charged constituents in liquids.

Like-likes-Like: Cooperative Hydrogen Bonding Overcomes Coulomb Repulsion in Cationic Clusters with Net Charges up to Q=+6e

Quantum chemical calculations have been employed to study kinetically stable cationic clusters, wherein the monovalent cations are trapped by hydrogen bonding despite strongly repulsive electrostatic forces. We calculated linear and cyclic clusters of the hydroxy-functionalized cation N-(3-hydroxypropyl) pyridinium, commonly used as cation in ionic liquids. The largest kinetically stable cluster was a cyclic hexamer that very much resembles the structural motifs of molecular clusters, as known for water and alcohols. Surprisingly, strong cooperative hydrogen bonds overcome electrostatic repulsion and result in cationic clusters with a high net charge up to Q=+6e. The structural, spectroscopic, and electronic signatures of the cationic and related molecular clusters of 3-phenyl-1-propanol could be correlated to NBO parameters, supporting the existence of "anti-electrostatic" hydrogen bonds (AEHB), as recently suggested by Weinhold. We also showed that dispersion forces enhance the cationic cluster formation and compensate the electrostatic repulsion of one additional positive charge.

Structural Motifs in Cold Ternary Ion Complexes of Hydroxyl-Functionalized Ionic Liquids: Isolating the Role of Cation-Cation Interactions

We address the competition between intermolecular forces underlying the recent observation that ionic liquids (ILs) with a hydroxyl-functionalized cation can form domains with attractive interactions between the nominally repulsive positively charged constituents. Here we show that this behavior is present even in the isolated ternary (HEMIm⁺)₂NTf₂^- complex (HEMIm⁺ = 1-(2-hydroxyethyl)-3-methylimidazolium) cooled to about 35 K in a photodissociation mass spectrometer. Of the three isomers isolated by double resonance techniques, one is identified to exhibit direct contact between the cations. This linkage involves a cooperative H-bond wherein the OH group on one cation binds to the OH group on the other, which then attaches to the basic N atom of the anion. Formation of this motif comes at the expense of the usually dominant interaction of the acidic C₍₂₎H group on the Im ring with molecular anions, as evidenced by isomer-dependent shifts in the C₍₂₎H vibrational fundamentals.

Theoretical and Experimental Insights into the Dissociation of 2-Hydroxyethylhydrazinium Nitrate Clusters Formed via Electrospray

Ionic liquids are used for myriad applications, including as catalysts, solvents, and propellants. Specifically, 2-hydroxyethylhydrazinium nitrate (HEHN) has been developed as a chemical propellant for space applications. The gas-phase behavior of HEHN ions and clusters is important in understanding its potential as an electrospray thruster propellant. Here, the unimolecular dissociation pathways of two clusters are experimentally observed, and theoretical modeling of hydrogen bonding and dissociation pathways is used to help rationalize those observations. The cation/deprotonated cation cluster [HEH₂ - H]⁺, which is observed from electrospray ionization, is calculated to be considerably more stable than the complementary cation/protonated anion adduct, [HEH + HNO₃]⁺, which is not observed experimentally. Upon collisional activation, a larger cluster [(HEHN)₂HEH]⁺ undergoes dissociation via loss of nitric acid at lower collision energies, as predicted theoretically. At higher collision energies, additional primary and secondary loss pathways open, including deprotonated cation loss, ion-pair loss, and double-nitric-acid loss. Taken together, these experimental and theoretical results contribute to a foundational understanding of the dissociation of protic ionic liquid clusters in the gas phase.

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.

Vibrational Spectroscopy of Ionic Liquids

Vibrational spectroscopy has continued use as a powerful tool to characterize ionic liquids since the literature on room temperature molten salts experienced the rapid increase in number of publications in the 1990's. In the past years, infrared (IR) and Raman spectroscopies have provided insights on ionic interactions and the resulting liquid structure in ionic liquids. A large body of information is now available concerning vibrational spectra of ionic liquids made of many different combinations of anions and cations, but reviews on this literature are scarce. This review is an attempt at filling this gap. Some basic care needed while recording IR or Raman spectra of ionic liquids is explained. We have reviewed the conceptual basis of theoretical frameworks which have been used to interpret vibrational spectra of ionic liquids, helping the reader to distinguish the scope of application of different methods of calculation. Vibrational frequencies observed in IR and Raman spectra of ionic liquids based on different anions and cations are discussed and eventual disagreements between different sources are critically reviewed. The aim is that the reader can use this information while assigning vibrational spectra of an ionic liquid containing another particular combination of anions and cations. Different applications of IR and Raman spectroscopies are given for both pure ionic liquids and solutions. Further issues addressed in this review are the intermolecular vibrations that are more directly probed by the low-frequency range of IR and Raman spectra and the applications of vibrational spectroscopy in studying phase transitions of ionic liquids.

Controlling the kinetic and thermodynamic stability of cationic clusters by the addition of molecules or counterions

We discuss the stability of cationic clusters when adding molecules or counterions, and predict their occurrence in gas phase experiments.

Observation of isolated ionic liquid cations and water molecules in an inert solvent

1-Octyl-3-methyl imidazolium halides ([OMIM]I and [OMIM]Cl) were loaded on top of CCl₄, and an in situ inclusion process was monitored from the CCl₄ phase as time elapses by infrared absorption spectroscopy. Absorption from the bands corresponding to the C(2)-H and C(4,5)-H stretch modes in the imidazolium cation was reduced significantly compared to the bulk IL spectra. This indicates that (1) the [OMIM] cation exists in CCl₄ as a monomer, dissociated from the anion and other cations, and (2) hydrogen bonding between the anion and the cation increases the dipole strength of the CH moieties in the imidazolium ring. [OMIM]I was found to transfer into the CCl₄ phase much faster than [OMIM]Cl, and this instigated us to compare the transfer of aqueous solutions of ionic liquids, 1-butyl-3-methyl imidazolium halides ([BMIM]I and [BMIM]Cl) into the CCl₄ matrix. Not only [BMIM]I but also water molecules transferred faster compared to those in the [BMIM]Cl aqueous solution. Water molecules in CCl₄ were shown to form clusters in [BMIM]I; presumably, I^- anions work as nucleation centers of water clusters.

Water bridges anchored by a C–H⋯O hydrogen bond: the role of weak interactions in molecular solvation

Hydrogen-bonded water bridges are re-directed from a polar NH bond to a weakly activated C(2)–H bond upon N-methylation. Infrared spectra, supported by ab initio calculations, provide direct evidence of the role of the C(2)–H donor in the solvation of the imidazole ring.

Interaction between ionic liquid cation and water: infrared predissociation study of [bmim]+·(H2O)n clusters

The infrared predissociation spectra of [bmim]⁺·(H₂O)_n, n = 1–8, in the 2800–3800 cm⁻¹ region are presented and analyzed with the help of electronic structure calculations.

Identification of multiple conformers of the ionic liquid [emim][tf2n] in the gas phase using IR/UV action spectroscopy

We have identified three families of conformers in gas phase ion pairs of [emim][tf2n] using IR/UV spectroscopy.

Understanding the ionic liquid [NC4111][NTf2] from individual building blocks: an IR-spectroscopic study

This study explores the interactions underlying the IR spectra of the ionic liquid [NC4111][NTf2] and its deuterated isotopomer [d9-NC4111][NTf2] by first isolating the spectra of charged ionic building blocks using mass-selective CIVP spectroscopy and then following the evolution of these bands upon sequential assembly of the ionic constituents. The spectra of the (1,1) and (2,2) neutral ion pairs are recorded using superfluid helium droplets as well as a solid neon matrix, while those of the larger charged aggregates are again obtained with CIVP. In general, the cluster spectra are similar to that of the bulk, with the (2,2) system displaying the closest resemblance. Analysis of the polarization-dependent band intensities of the neutral ion pairs in liquid droplets as a function of external electric field yields dipole moments of the neutral aggregates. This information allows a coarse assessment of the packing structure of the neutral pairs to be antiparallel at 0.37 K, in contrast to the parallel arrangement found for the assembly of small, high-dipole neutral molecules with large rotational constants (e.g., HCN). The role of an extra anion or cation attached to both the (1,1) and the (2,2) ion pairs to form the charged clusters is discussed in the context of an additional remote, more unfavorable binding site intrinsic to the nature of the charged IL clusters and as such not anticipated in the bulk phase. Whereas for the anion itself only the lowest energy trans conformer was observed, the higher clusters showed an additional population of the cis conformer. The interactions are found to be consistent with a minimal role of hydrogen bonding.

Conformational preferences of monohydrated clusters of imidazole derivatives revisited

IR-UV double resonance spectroscopy was used to identify the conformers of monohydrated benzimidazole andN-methylbenzimidazole in a supersonic jet. A new OH–N bound conformer relevant to histidine containing proteins was discovered. The long standing differences in the literature about the relative energies and abundance of the monohydrated imidazole derivatives have also been resolved.