Translational Diffusion in a Set of Imidazolium-Based Ionic Liquids [bmim]+A- and Their Mixtures with Water

Do Electrostatics Control the Diffusive Dynamics of Solitary Water? NMR and MD Studies of Water Translation and Rotation in Dipolar and Ionic Solvents

NMR-based measurements of the diffusion coefficients and rotation times of solitary water and benzene at 300 K are reported in a diverse collection of 13 conventional organic solvents and 10 imidazolium ionic liquids. Proton chemical shifts of water are found to be correlated to water OH-stretching frequencies, confirming the importance of electrostatic interactions in these shifts. However, the influence of magnetic interactions in aromatic solvents renders chemical shifts a less reliable indicator of electrostatics. Diffusion coefficients (DB) and rotational correlation times (τB) of benzene in the solvents examined are accurately described as functions of viscosity (η) by DB ∝ η-0.81 and τB ∝ η0.64. Literature values of DB and τB in alkane and normal alcohols, which were not included among the solvents studied here, are systematically faster than predicted by these correlations, indicating that factors beyond solvent viscosity play a role in determining the friction on benzene. In contrast to benzene, water diffusion and rotation are poorly described in terms of viscosity alone, even in the dipolar and ionic solvents measured here. The present data and the substantial literature data already available on dilute water diffusion show a systematic dependence of DW on solvent polarity among isoviscous solvents. The aspect of solvent polarity most relevant to water dynamics is the ability of a solvent to accept hydrogen bonds from water, as conveniently quantified by the frequency of water's OH stretching band, ΔνOH. The friction on translation, ζtr = kBT/DW, and rotation, ζrot = kBTτW, are both well correlated by functions of the form ζ(η, ΔνOH) = a1ηa2 exp (a3ΔνOH), where the ai are adjustable parameters. Molecular dynamics simulations reveal a strong coupling between electrostatic and nonelectrostatic water-solvent interactions, which makes it impossible to dissect the friction on water into additive dielectric and hydrodynamic components. Simulations also provide a tentative explanation for the unusual form of the correlating function ζ(η, ΔνOH), at least in the case of ζrot.

The State of Water in "Ionic Liquid [bmim]Cl/AlCl3/H2O" Systems

The unique physical and chemical properties of ionic liquids (ILs) determine their numerous applications in "green" chemistry and material science. Recently, systems based on ILs have been considered to be promising for use in a new generation of electrochemical devices. The results of a nuclear magnetic resonance (NMR) study of the microstructure of 1-butyl-3-methylimidazolium chloride (IL)/water mixtures in the presence of Al3+ cations are presented. For the first time, the splitting of spectral lines of water in such systems has been recorded. Comparing the 1H and 27Al NMR data, we have detected the existence of different solvate complexes of Al3+ with Cl- and estimated the characteristic times of exchange processes. For the system under study, a model of the Al3+ cation environment and its evolution with temperature and water content has been described. Quantum-chemical calculations have been performed to substantiate the model.

Anomalous Dependence of Translational Diffusion on the Water Mole Fraction for Solute Molecules Dissolved in a 1-Butyl-3-methylimidazolium Tetrafluoroborate/Water Mixture

Translational diffusion coefficients of carbon monoxide (CO), diphenylacetylene (DPA), and diphenylcyclopropenone (DPCP) were determined in mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate ([C₄mim]BF₄) and water using transient grating spectroscopy at different mole fractions of water (x_w). While DPA exhibited a larger diffusion coefficient than DPCP at low water mole fractions (x_w < 0.7), as observed for conventional liquids and ionic liquids (ILs), it was smaller at high mole fractions (x_w > 0.9). The apparent molecular radius of DPA determined using the Stokes-Einstein equation at x_w > 0.9 is close to the radius of an IL cluster in a water pool as determined from small-angle neutron scattering experiments (J. Bowers et al., Langmuir, 2004, 20, 2192-2198), suggesting that the DPA molecules are trapped in IL clusters in the water pool and move together. The solvation state of DPCP in the mixture was studied using Raman spectroscopy. Dramatically strong water/DPCP hydrogen bonding was observed at higher water mole fractions, suggesting that DPCP is located near the cluster interfaces. The large diffusion coefficient of DPCP suggests that hopping of DPCP between IL clusters occurs through hydrogen bonding with water.

Solvation structure and dynamics of coumarin 153 in an imidazolium-based ionic liquid with chloroform, benzene, and propylene carbonate

In order to determine the self-diffusion coefficients D of all the species in the solutions at 298.2 K, ¹H and ¹⁹F NMR diffusion ordered spectroscopy (DOSY) has been conducted on coumarin 153 (C153) in binary mixed solvents of an imidazolium-based ionic liquid (IL), 1-dodecyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide (C₁₂mimTFSA), with three molecular liquids (MLs) of chloroform (CL), benzene (BZ), and propylene carbonate (PC) as a function of ML mole fraction x_ML. Below x_ML ≈ 0.8, the D values of each species do not significantly depend on the MLs. However, above this mole fraction, the diffusion of C153 becomes smoother in the order of BZ ≈ CL > PC systems. The interactions among C153, C₁₂mim⁺, TFSA^-, and ML molecules have been investigated using infrared (IR) and ¹H and ¹³C NMR spectroscopic techniques. The relations of the diffusion of the species with the interactions among them have been discussed on the molecular scale. In the IL solution, the C153 carbonyl oxygen atom is hydrogen-bonded with the imidazolium ring C²-H atom of C₁₂mim⁺. C₁₂mim⁺ also forms an ion pair with TFSA^-. Thus, C153, C₁₂mim⁺, and TFSA^- cooperatively move in the CL and BZ solutions at a lower ML content, x_ML < ∼0.8. On the other hand, at a higher ML content, x_ML > ∼0.8, the C153 molecule diffuses with CL and BZ molecules because of the hydrogen bonding between the C153 carbonyl O atom and the CL H atom and the π-π interaction between the C153 and BZ ring planes, respectively. For the PC system, the change in the relative self-diffusion coefficients of each species with increasing x_ML differs from those for the CL and BZ systems because of both hydrogen bonding donor H and acceptor O atoms of PC for C153, the IL cation and anion, and PC themselves.

Recent advances in NMR spectroscopy of ionic liquids

This review presents recent developments in the application of NMR spectroscopic techniques in the study of ionic liquids. NMR has been the primary tool not only for the structural characterization of ionic liquids, but also for the study of dynamics. The presence of a host of NMR active nuclei in ionic liquids permits widespread use of multinuclear NMR experiments. Chemical shifts and multinuclear coupling constants are used routinely for the structure elucidation of ionic liquids and of products formed by their covalent interactions with other materials. Also, the availability of a multitude of NMR techniques has facilitated the study of dynamical processes in them. These include the use of NOESY to study inter-ionic interactions, pulsed-field gradient techniques for probing transport properties, and relaxation measurements to elucidate rotational dynamics. This review will focus on the application of each of these techniques to investigate ionic liquids.

Structural Features of the [C4mim][Cl] Ionic Liquid and Its Mixtures with Water: Insight from a 1H NMR Experimental and QM/MD Study

The 1H NMR chemical shift of water exhibits non-monotonic dependence on the composition of an aqueous mixture of 1-butyl-3-methylimidazolium chloride, [C4mim][Cl], ionic liquid (IL). A clear minimum is observed for the 1H NMR chemical shift at a molar fraction of the IL of 0.34. To scrutinize the molecular mechanism behind this phenomenon, extensive classical molecular dynamics simulations of [C4mim][Cl] IL and its mixtures with water were carried out. A combined quantum mechanics/molecular mechanics approach based on the density functional theory was applied to predict the NMR chemical shifts. The proliferation of strongly hydrogen-bonded complexes between chloride anions and water molecules is found to be the reason behind the increasing 1H NMR chemical shift of water when its molar fraction in the mixture is low and decreasing. The model shows that the chemical shift of water molecules that are trapped in the IL matrix without direct hydrogen bonding to the anions is considerably smaller than the 1H NMR chemical shift predicted for the neat water. The structural features of neat IL and its mixtures with water have also been analyzed in relation to their NMR properties. The 1H NMR spectrum of neat [C4mim][Cl] was predicted and found to be in very reasonable agreement with the experimental data. Finally, the experimentally observed strong dependence of the chemical shift of the proton at position 2 in the imidazolium ring on the composition of the mixture was rationalized.

Unravelling the Interactions of Magnetic Ionic Liquids by Energy Decomposition Schemes: Towards a Transferable Polarizable Force Field

This work aims at unravelling the interactions in magnetic ionic liquids (MILs) by applying Symmetry-Adapted Perturbation Theory (SAPT) calculations, as well as based on those to set-up a polarisable force field model for these liquids. The targeted MILs comprise two different cations, namely: 1-butyl-3-methylimidazolium ([Bmim]⁺) and 1-ethyl-3-methylimidazolium ([Emim]⁺), along with several metal halides anions such as [FeCl₄]⁻, [FeBr₄]⁻, [ZnCl₃]⁻ and [SnCl₄]²⁻ To begin with, DFT geometry optimisations of such MILs were performed, which in turn revealed that the metallic anions prefer to stay close to the region of the carbon atom between the nitrogen atoms in the imidazolium fragment. Then, a SAPT study was carried out to find the optimal separation of the monomers and the different contributions for their interaction energy. It was found that the main contribution to the interaction energy is the electrostatic interaction component, followed by the dispersion one in most of the cases. The SAPT results were compared with those obtained by employing the local energy decomposition scheme based on the DLPNO-CCSD(T) method, the latter showing slightly lower values for the interaction energy as well as an increase of the distance between the minima centres of mass. Finally, the calculated SAPT interaction energies were found to correlate well with the melting points experimentally measured for these MILs.

Rapid relaxation NMR measurements to predict rate coefficients in ionic liquid mixtures. An examination of reaction outcome changes in a homologous series of ionic liquids

A series of ionic liquids based on the 1-alkyl-3-methylimidazolium cations were examined as components of the solvent mixture for a bimolecular substitution process. The effects on both the rate coefficient of the process and the NMR spin-spin relaxation of the solvent components of changing either the alkyl chain length or the amount of ionic liquid in the reaction mixture were determined. At a constant mole fraction, a shorter alkyl chain length resulted in a greater rate coefficient enhancement and a longer relaxation time, with the opposite effects for a longer alkyl chain length. For a given ionic liquid, increasing the proportion of salt in the reaction mixture resulted in a greater rate coefficient and a shorter relaxation time. The microscopic origins of the rate coefficient enhancement were determined and a step change found in the activation parameters on increasing the alkyl chain length from hexyl to octyl, suggesting notable structuring in solution. Across a range of ionic liquids and solvent compositions, the relaxation time from NMR measurements was shown to relate to the reaction rate coefficient. The approach of using fast and simple NMR relaxation measurements to predict reaction outcomes was exemplified using a morpholinium-based ionic liquid.

Structure, Molecular Interactions, and Dynamics of Aqueous [BMIM][BF4] Mixtures: A Molecular Dynamics Study

Molecular dynamics simulations with many-body polarizable force fields were carried out to investigate the thermodynamic, structural, and dynamic properties of aqueous solutions of 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF₄]). The radial distribution functions exhibit well-defined features, revealing favored structural correlations between [bmim]⁺, [BF₄]^-, and H₂O. The addition of water is shown to alter ionic liquid structural organizations by replacing counterions in the coordination shells and disrupt the cation-anion network. At low water concentration, the majority of water molecules are isolated from each other and have lower average dipole moment than that in pure water. With increasing hydration level, while [bmim][BF₄] ionic network breaks up and becomes isolated ion pairs or free ions in the dilute limit, water begins to form clusters of increasing sizes and eventually forms a percolating network. As a result, the average water dipole moment increases and approaches its bulk value. Water is also observed to have a substantial influence on the dynamics of ionic liquids. At low water content, the cation and anion have similar diffusion coefficients due to the correlated ionic motion of long-lived ion pairs. As the water concentration increases, both ions exhibit greater mobility and faster rotations from the breakup of ionic network. Consequently, the ionic conductivity of [bmim][BF₄] aqueous solutions rises with increasing water composition.

Voronoi Polyhedra as a Tool for the Characterization of Inhomogeneous Distribution in 1-Butyl-3-methylimidazolium Cation-Based Ionic Liquids

The inhomogeneity distribution in four imidazolium-based ionic liquids (ILs) containing the 1-butyl-3-methylimidazolium (C₄mim) cation, coupled with tetrafluoroborate (BF₄), hexafluorophosphate (PF₆), bis(trifluoromethanesulfonyl)amide (TFSA), and trifluoromethanesulfonate (TfO) anions, was characterized using Voronoi polyhedra. For this purpose, molecular dynamic simulations have been performed on the isothermal-isobaric (NpT) ensemble. We checked the ability of the potential models to reproduce the experimental density, heat of vaporization, and transport properties (diffusion and viscosity) of these ionic liquids. The inhomogeneity distribution of ions around the ring, methyl, and butyl chain terminal hydrogen atoms of the C₄mim cation was investigated by means of Voronoi polyhedra analysis. For this purpose, the position of the C₄mim cation was described successively by the ring, methyl, and butyl chain terminal hydrogen atoms, while that of the anions was described by their F or O atom. We calculated the Voronoi polyhedra distributions of the volume, the density, and the asphericity parameters as well as that of the radius of the spherical intermolecular voids. We carried out the analysis in two steps. In the first step, both ions were taken into account. The calculated distributions gave information on the neighboring ions around a reference one. In the second step, to distinguish between like and oppositely charged ions and then to get information on the inhomogeneity distribution of the like ions, we repeated the same calculations on the same sample configurations and removed one of the ions and considered only the other one. Detailed analysis of these distributions has revealed that the ring hydrogen atoms are mainly solvated by the anions, while the methyl and butyl terminal H atoms are surrounded by like atoms. The extent of this inhomogeneity was assessed by calculating the cluster size distribution that shows that the dimers are the most abundant ones.

Revisiting Ionic Liquid Structure-Property Relationship: A Critical Analysis

In the last few years, ionic liquids (ILs) have been the focus of extensive studies concerning the relationship between structure and properties and how this impacts their application. Despite a large number of studies, several topics remain controversial or not fully answered, such as: the existence of ion pairs, the concept of free volume and the effect of water and its implications in the modulation of ILs physicochemical properties. In this paper, we present a critical review of state-of-the-art literature regarding structure–property relationship of ILs, we re-examine analytical theories on the structure–property correlations and present new perspectives based on the existing data. The interrelation between transport properties (viscosity, diffusion, conductivity) of IL structure and free volume are analysed and discussed at a molecular level. In addition, we demonstrate how the analysis of microscopic features (particularly using NMR-derived data) can be used to explain and predict macroscopic properties, reaching new perspectives on the properties and application of ILs.

Molecular Dynamics Insights into the Nanoscale Structural Organization and Local Interaction of Aqueous Solutions of Ionic Liquid 1-Butyl-3-methylimidazolium Nitrate

Considering the growing number of applications of the aqueous ionic liquids (ILs), atomistic molecular dynamics (MD) simulations were used to probe the effect of water molar fraction, x_w, ranging from 0.00 to 0.90, on the nanoscale local structure of 1-butyl-3-methylimidazolium nitrate, [bmim][NO₃], IL. The results prove that, with water addition, the cation-anion, cation-cation, and anion-anion structural correlations are weakened, while strong anion-water and unconventional cation-water hydrogen bonds are formed in the solutions. Water molecules were detected as bridges between nitrate anions, and the water cluster size distribution at different x_w's was investigated. Simulation shows a similar pattern of probability densities for water and anion around the acidic hydrogen atoms of the reference cation ring, while both species move away from the cation butyl chain. Increasing the water concentration to x_w = 0.90 causes decreasing of the local arrangement of the nearest-neighboring cations, because of the weakening of cation-cation π-π stacking. In addition, this dilution reduces the probability of the in-plane cation-anion conformation, disrupts both the polar ionic network and nonpolar domains, and diminishes the nanoaggregation of the cation butyl chains compared to those of the neat IL. These results can rationalize the origins of the fluidity enhancements and transport property trends upon adding water to the imidazolium-based ILs. The current study proposes a deep atomistic-level insight into the complex coupling between water concentration, microscopic structure, and local interactions of aqueous imidazolium-based ILs with hydrophilic anions.