Extension of the CL&Pol Polarizable Force Field to Electrolytes, Protic Ionic Liquids, and Deep Eutectic Solvents.

The polarizable CL&Pol force field presented in our previous study, Transferable, Polarizable Force Field for Ionic Liquids (J. Chem. Theory Comput. 2019, 15, 5858, DOI: http://doi.org/10.1021/acs.jctc.9b0068910.1021/acs.jctc.9b00689), is extended to electrolytes, protic ionic liquids (PIL), deep eutectic solvents (DES), and glycols. These systems are problematic in polarizable simulations because they contain either small, highly charged ions or strong hydrogen bonds, which cause trajectory instabilities due to the pull exerted on the induced dipoles. We use a Tang-Toennies (TT) function to dampen, or smear, the interactions between charges and induced dipole at a short range involving small, highly charged atoms (such as hydrogen or lithium), thus preventing the "polarization catastrophe". The new force field gives stable trajectories and is validated through comparison with experimental data on density, viscosity, and ion diffusion coefficients of liquid systems of the above-mentioned classes. The results also shed light on the hydrogen-bonding pattern in ethylammonium nitrate, a PIL, for which the literature contains conflicting views. We describe the implementation of the TT damping function, of the temperature-grouped Nosé-Hoover thermostat for polarizable molecular dynamics (MD) and of the periodic perturbation method for viscosity evaluation from non-equilibrium trajectories in the LAMMPS MD code. The main result of this work is the wider applicability of the CL&Pol polarizable force field to new, important classes of fluids, achieving robust trajectories and a good description of equilibrium and transport properties in challenging systems. The fragment-based approach of CL&Pol will allow ready extension to a wide variety of PILs, DES, and electrolytes.

Sodium diffusion in ionic liquid-based electrolytes for Na-ion batteries: the effect of polarizable force fields.

Understanding the transport of sodium ions in ionic liquids is key to designing novel electrolyte materials for sodium-ion batteries. In this work, we combine molecular dynamics simulation and experiments to study how molecular interactions and local ordering affect relevant physico-chemical properties. Ionic transport and local solvation environments are investigated in electrolytes composed of sodium bis(fluorosulfonyl)imide, (Na[FSI]), in N,N-methylpropylpyrrolidinium bis(fluorosulfonyl)imide, [C3C1pyr][FSI], at different salt concentrations. The electrolyte systems are modelled by means of molecular dynamic simulations using a polarizable force field. We show that including polarization effects explicitly in the molecular simulations is required in order to attain a reliable description of the transport properties of sodium in the [C3C1pyr][FSI] electrolyte. The validation of the computational results upon comparison with experimental data allows us to assess the suitability of polarizable force fields in describing and interpreting the structure and dynamics of the sodium salt-ionic liquid system, which is essential to enable the application of IL-based electrolytes in novel energy-storage technologies.

Self-assembled nanostructures in ionic liquids facilitate charge storage at electrified interfaces.

Driven by the potential applications of ionic liquids (ILs) in many emerging electrochemical technologies, recent research efforts have been directed at understanding the complex ion ordering in these systems, to uncover novel energy storage mechanisms at IL-electrode interfaces. Here, we discover that surface-active ILs (SAILs), which contain amphiphilic structures inducing self-assembly, exhibit enhanced charge storage performance at electrified surfaces. Unlike conventional non-amphiphilic ILs, for which ion distribution is dominated by Coulombic interactions, SAILs exhibit significant and competing van der Waals interactions owing to the non-polar surfactant tails, leading to unusual interfacial ion distributions. We reveal that, at an intermediate degree of electrode polarization, SAILs display optimum performance, because the low-charge-density alkyl tails are effectively excluded from the electrode surfaces, whereas the formation of non-polar domains along the surface suppresses undesired overscreening effects. This work represents a crucial step towards understanding the unique interfacial behaviour and electrochemical properties of amphiphilic liquid systems showing long-range ordering, and offers insights into the design principles for high-energy-density electrolytes based on spontaneous self-assembly behaviour.

Molecular dynamics simulations of polyethers and a quaternary ammonium ionic liquid as CO2 absorbers.

The properties of mixtures of butyltrimethylammonium bis(trifluoromethylsulfonyl)imide, [N4111][NTf2], with poly(ethyleneglycol) dimethyl ether, PEO, were described as a function of PEO chain size by molecular dynamics simulations. Both PEO chain size and mixture composition revealed to play a significant role in determining the structure and the dynamics of the fluids. The remarkably higher viscosity observed for mixtures composed by 0.25 mole fraction of PEO was attributed to the increase in the gauche population of OCCO dihedral of the polyether of longer chains. The negative solvation enthalpy (ΔsolH < 0) and entropy (ΔsolS < 0) revealed a favorable CO2 absorption by the neat and mixture systems. The CO2 absorption was higher in neat PEO, particularly considering longer chains. The gas solubility in the mixtures presented intermediate values in comparison to the neat PEO and neat ionic liquid. The CO2 solutions had their structures discussed in the light of the calculated radial and spatial distribution functions.

Interactions between water and 1-butyl-1-methylpyrrolidinium ionic liquids.

We report experimental results on the diffusivity of water in two ionic liquids obtained using the pulsed-gradient spin-echo NMR method. Both ionic liquids have the same cation, 1-butyl-1-methylpyrrolidinium, but different trifluoromethyl-containing anions. One has a strongly hydrophobic anion, bis(trifluoromethylsulfonyl)amide, while the second has a hydrophilic anion, trifluoromethylsulfonate. Transport of water in these ionic liquids is much faster than would be predicted from hydrodynamic laws, indicating that the neutral water molecules experience a very different friction than the anions and cations at the molecular level. Temperature-dependent viscosities, conductivities, and densities are reported as a function of water concentration to further analyze the properties of the ionic liquid-water mixtures. These results on the properties of water in ionic liquids should be of interest to researchers in diverse areas ranging from separations, solubilizing biomass and energy technologies.

Preliminary study on suitability of ionic liquids as potential passive-sampling media of polyaromatic-hydrocarbon (PAH) analyses in water.

Recently, ionic liquids (ILs) have been regarded as an attractive water-immiscible phase in liquid-liquid extraction. Because ILs have a wide range of polarity irrespective of their miscibility with water, the possibility of using them as an effective extraction phase for a broad range of contaminants means they are starting to be of particular interest. In this study we investigated a wide variety of ionic liquids, which are known to be hydrolytically stable and of a hydrophobic character, for their potential suitability as passive-sampling media for monitoring selected polyaromatic hydrocarbons. Preliminary research in this field has indicated very promising results using these novel extraction media. Because there is an enormous number of possible cation-anion combinations offering tuneable properties of ionic liquids with the potential for effective passive extraction, we hope this paper will encourage the scientific community to undertake further studies verifying the undoubted usefulness of these alternative solvents as passive samplers for many other groups of analytes. Additionally, because of the unusual solubility properties that have already been proved for ILs, it is very probable that it would soon be possible to deliver a very effective system able to extract analytes differing widely in polarity.

Solvation of a Cellulose Microfibril in Imidazolium Acetate Ionic Liquids: Effect of a Cosolvent.

The solvation and the onset of dissolution of a cellulose I(ß) microcrystal in ionic liquid media are studied by molecular simulation. Ionic liquids can dissolve large amounts of cellulose, which can later be regenerated from solution, but their high viscosity is an inconvenience. Hydrogen bonding between the anion of the ionic liquid and cellulose is the main aspect determining dissolution. Here we try to elucidate the role of a molecular cosolvent, dimethyl sulfoxide (DMSO), which is an aprotic polar compound, in the system composed of cellulose and the ionic liquid 1-butyl-3-methylimidazolium acetate. We calculated quantities related to specific interactions (mainly hydrogen bonds), conformations, and the structure of local solvation environments, both for a solvated oligomer chain of cellulose and for a model microfibril composed of 36 chains in the I(ß) crystal structure. We compare two solvent systems: the pure ionic liquid and a mixed solvent with an equimolar composition in ionic liquid and DMSO. All entities are represented by detailed all-atom, fully flexible force fields. The main conclusions are that DMSO behaves as an "innocent" cosolvent, lowering the viscosity and accelerating mass transport in the system, but without interacting specifically with cellulose or disrupting the interactions between cellulose with the anions of the ionic liquid. An understanding of solvation in mixed solvents composed of ionic liquids and molecular compounds can enable the design of high-performance media for the use of biomass materials.

Effect of unsaturation on the absorption of ethane and ethylene in imidazolium-based ionic liquids.

The influence of the presence of imidazolium side chain unsaturation on the solubility of ethane and ethylene was studied in three ionic liquids: 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide-saturated alkyl side-chain in the cation; 1-methyl-3-(buten-3-yl)imidazolium bis(trifluorosulfonyl)imide-double bond in the side-chain of the cation; and 1-methyl-3-benzylimidazolium bis(trifluorosulfonyl)imide-benzyl group in the side-chain of the cation. The solubility of both gases decreases when the side-chain of the cations is functionalized with an unsaturated group. This can be explained by a less favorable enthalpy of solvation. The difference of solubility between ethane and ethylene can be explained from a balance of enthalpic and entropic factors: for the ionic liquid with the saturated alkyl side-chain and the benzyl-substituted side-chain, it is the favorable entropy of solvation that explains the larger ethylene solubility, whereas in the case of the saturated side-chain, it is the more favorable enthalpy of solvation. Molecular simulation allowed the identification of the mechanisms of solvation and the preferential solvation sites for each gas in the different ionic liquids. Simulations have shown that the entropy of solvation is more favorable when the presence of the gas weakens the cation-anion interactions or when the gas can be solvated near different sites of the ionic liquid.

Effect of water on the carbon dioxide absorption by 1-alkyl-3-methylimidazolium acetate ionic liquids.

The absorption of carbon dioxide by the pure ionic liquids 1-ethyl-3-methylimidazolium acetate ([C(1)C(2)Im][OAc]) and 1-butyl-3-methylimidazolium acetate ([C(1)C(4)Im][OAc]) was studied experimentally from 303 to 343 K. As expected, the mole fraction of absorbed carbon dioxide is high (0.16 at 303 K and 5.5 kPa and 0.19 at 303 and 9.6 KPa for [C(1)C(2)Im][OAc] and [C(1)C(4)Im][OAc], respectively), does not obey Henry's law, and is compatible with the chemisorption of the gas by the liquid. Evidence of a chemical reaction between the gas and the liquid was found both by NMR and by molecular simulation. In the presence of water, the properties of the liquid absorber significantly change, especially the viscosity that decreases by as much as 25% (to 78 mPa s) and 30% (to 262 mPa s) in the presence of 0.2 mol fraction of water for [C(1)C(2)Im][OAc] and [C(1)C(2)Im][OAc] at 303 K, respectively. The absorption of carbon dioxide decreases when the water concentration increases: a decrease of 83% in CO(2) absorption is found for [C(1)C(4)Im][OAc] with 0.6 mol fraction of water at 303 K. It is proved in this work, by combining experimental data with molecular simulation, that the presence of water not only renders the chemical reaction between the gas and the ionic liquid less favorable but also lowers the (physical) solubility of the gas as it competes by the same solvation sites of the ionic liquid. The lowering of the viscosity of the liquid absorbent largely compensates these apparent drawbacks and the mixtures of [C(1)C(2)Im][OAc] and [C(1)C(2)Im][OAc] with water seem promising to be used for carbon dioxide capture.

Ligand effect on the catalytic activity of ruthenium nanoparticles in ionic liquids.

Suspensions of small sized (1-2.5 nm) ruthenium nanoparticles (RuNPs) have been obtained by decomposition, under H(2), of (η(4)-1,5-cyclooctadiene)(η(6)-1,3,5-cyclooctatriene)ruthenium(0), [Ru(COD)(COT)], in the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(4)Im][NTf(2)], and in the presence of different compounds acting as ligands: C(8)H(17)NH(2), PPhH(2), PPh(2)H and H(2)O. Previous and new liquid NMR experiments showed that the ligands are coordinated or in the proximity to the surface of the RuNPs. Herein is reported how the ligand affects the catalytic performance (activity and selectivity) compared to a ligand-free system of RuNPs, when RuNPs in [C(1)C(4)Im][NTf(2)] are used as catalysts for the hydrogenation of various unsaturated compounds (1,3-cyclohexadiene, limonene and styrene). It has been observed that σ-donor ligands increase the activity of the nanoparticles, contrarily to π-acceptor ones.

Relevant parameters for assessing the environmental impact of some pyridinium, ammonium and pyrrolidinium based ionic liquids.

Several physico-chemical properties relevant to determine the environmental impact of ionic liquids - aqueous solubility, octanol/water partition coefficient, chromatographically derived lipophilicity and infinite dilution diffusion coefficients in water - were measured in ionic liquids based on pyridinium, ammonium and pyrrolidinium cations with bis(trifluoromethylsulfonyl)imide anions. The influence of the presence of hydroxyl or ester groups in the physico-chemical properties of these liquids was checked. It appeared that the presence of functional oxygenated moieties reduces the lipophilicity of ionic liquids and so decreases the risk of bioaccumulation in environment.

Glycine in 1-butyl-3-methylimidazolium acetate and trifluoroacetate ionic liquids: effect of fluorination and hydrogen bonding.

The solvation of glycine in two ionic liquids (ILs), namely, 1-butyl-3-methylimidazolium acetate, [C(1)C(4)Im][OAc], and 1-butyl-3-methylimidazolium trifluoroacetate, [C(1)C(4)Im][TFA], was studied by a combination of experimental and theoretical methods. The solubility of glycine in both ILs was determined at 333.15 K to be (8.1±0.5) and (1.0±0.5) wt % in [C(1)C(4)Im][OAc] and [C(1)C(4)Im][TFA], respectively. By IR spectroscopy it was found that, when dissolved in the ILs, glycine was mainly present in its zwitterionic form. Structural and energetic aspects of the solvation of glycine in the ILs and in mixtures of ILs and water were investigated by ab initio calculations and molecular dynamic simulations. It was observed that the firstly solvation shell around glycine consisted predominantly of acetate or trifluoroacetate anions, which formed hydrogen bonds either with the carboxylic group of neutral glycine or with the protonated ammonium group of the zwitterionic form. When water is present in the solutions, hydrogen bonds between water and the anion prevail. The overall energy of the system was decomposed into its components between pairs of species. It was established that the dominant contribution to the interaction energy between glycine and the IL was due to hydrogen bonds with the anions and the statistics of hydrogen bonds were analysed.

Phase behaviour, interactions, and structural studies of (amines+ionic liquids) binary mixtures.

We present a study on the phase equilibrium behaviour of binary mixtures containing two 1-alkyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}imide-based ionic liquids, [C(n)mim] [NTf(2)] (n=2 and 4), mixed with diethylamine or triethylamine as a function of temperature and composition using different experimental techniques. Based on this work, two systems showing an LCST and one system with a possible hourglass shape are measured. Their phase behaviours are then correlated and predicted by using Flory-Huggins equations and the UNIQUAC method implemented in Aspen. The potential of the COSMO-RS methodology to predict the phase equilibria was also tested for the binary systems studied. However, this methodology is unable to predict the trends obtained experimentally, limiting its use for systems involving amines in ionic liquids. The liquid-state structure of the binary mixture ([C(2)mim] [NTf(2)]+diethylamine) is also investigated by molecular dynamics simulation and neutron diffraction. Finally, the absorption of gaseous ethane by the ([C(2)mim][NTf(2)]+diethylamine) binary mixture is determined and compared with that observed in the pure solvents.

Using ethane and butane as probes to the molecular structure of 1-alkyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl] imide ionic liquids.

In this work, we have studied the solubility and the thermodynamic properties of solvation, between 298 and 343 K and at pressures close to atmospheric, of ethane and n-butane in several ionic liquids based on the bis[(trifluoromethyl) sulfonyl]imide anion and on 1-alkyl-3-methylimidazolium cations, [CnC1Im] [NTf2], with alkyl side-chains varying from two to ten carbon atoms. The solubility of butane is circa one order of magnitude larger than that of ethane with mole fractions as high as 0.15 in [C10C1Im][NTf2] at 300 K. The solubilities of both n-butane and ethane gases are higher for ionic liquids with longer alkyl chains. The behaviour encountered is explained by the preferential solvation of the gases in the non-polar domains of the solvents, the larger solubility of n-butane being attributed to the dispersive contributions to the interaction energy. The rise in solubility with increasing size of the alkyl-side chain is explained by a more favourable entropy of solvation in the ionic liquids with larger cations. These conclusions are corroborated by molecular dynamics simulation studies.

Influence of ionic association, transport properties, and solvation on the catalytic hydrogenation of 1,3-cyclohexadiene in ionic liquids.

The influence of the nature of two different ionic liquids, namely 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(4)Im][NTf(2)], and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(1)C(4)Im][NTf(2)], on the catalytic hydrogenation of 1,3-cyclohexadiene with [Rh(COD)(PPh(3))(2)][NTf(2)] (COD = 1,5-cyclooctadiene) was studied. Initially, the effect of different concentrations of 1,3-cyclohexadiene on the molecular interactions and on the structure in two ionic liquids was investigated by NMR and by molecular dynamic simulations. It was found that in both ionic liquids 1,3-cyclohexadiene is solvated preferentially in the lipophilic regions. Furthermore, the higher solubility of 1,3-cyclohexadiene in [C(1)C(4)Im][NTf(2)] and the smaller positive values of the excess molar enthalpy of mixing for the 1,3-cyclohexadiene + [C(1)C(4)Im][NTf(2)] system in comparison with 1,3-cyclohexadiene + [C(1)C(1)C(4)Im][NTf(2)] indicate more favorable interactions between 1,3-cyclohexadiene and the C(1)C(4)Im(+) cation than with the C(1)C(1)C(4)Im(+) cation. Subsequently, diffusivity and conductivity measurements of the 1,3-cyclohexadiene + ionic liquid mixtures at different compositions allowed a characterization of mass and charge transport in the media and access to the ionicity of ionic liquids in the mixture. From the dependence of the ratio between molar conductivity and the conductivity inferred from NMR diffusion measurements, Λ(imp)/Λ(NMR), on concentration of 1,3-cyclohexadiene in the ionic liquid mixture, it was found that increasing the amount of 1,3-cyclohexadiene leads to a decrease in the ionicity of the medium. Finally, the reactivity of the catalytic hydrogenation of 1,3-cyclohexadiene using [Rh(COD)(PPh(3))(2)][NTf(2)] performed in [C(1)C(4)Im][NTf(2)] at different compositions of 1,3-cyclohexadiene and in [C(1)C(1)C(4)Im][NTf(2)] at one composition was related linearly to the viscosity, hence the reaction rate is determined by the mass transport properties of the media.

Effect of alkyl chain length and hydroxyl group functionalization on the surface properties of imidazolium ionic liquids.

Properties of the surface of ionic liquids, such as surface tension, ordering, and charge and density profiles, were studied using molecular simulation. Two types of modification in the molecular structure of imidazolium cations were studied: the length of the alkyl side chain and the presence of a polar hydroxyl group at the end of the side chain. Four ionic liquids were considered: 1-ethyl-3-methylimidazolium tetrafluoroborate, [C(2)C(1)im][BF(4)]; 1-(2-hydroxyethyl)-3-methylimidazolium tetrafluoroborate, [C(2)OHC(1)im][BF(4)]; 1-octyl-3-methylimidazolium tetrafluoroborate, [C(8)C(1)im][BF(4)] and 1-(8-hydroxyoctyl)-3-methylimidazolium tetrafluoroborate, [C(8)OHC(1)im][BF(4)]. The surface tension was calculated using both mechanical and thermodynamic definitions, with consistent treatment of the long-range corrections. The simulations reproduce the available experimental values of surface tension with a maximum deviation of ±10%. This energetic characterization of the interface is completed by microscopic structural analysis of orientational ordering at the interface and density profiles along the direction normal to the interface. The presence of the hydroxyl group modifies the local structure at the interface, leading to a less organized liquid phase. The results allow us to relate the surface tension to the structural ordering at the liquid-vacuum interface.

Ruthenium nanoparticles in ionic liquids: structural and stability effects of polar solutes.

Ionic liquids are a stabilizing medium for the in situ synthesis of ruthenium nanoparticles. Herein we show that the addition of molecular polar solutes to the ionic liquid, even in low concentrations, eliminates the role of the ionic liquid 3D structure in controlling the size of ruthenium nanoparticles, and can induce their aggregation. We have performed the synthesis of ruthenium nanoparticles by decomposition of [Ru(COD)(COT)] in 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(4)Im][NTf(2)], under H(2) in the presence of varying amounts of water or 1-octylamine. For water added during the synthesis of metallic nanoparticles, a decrease of the solubility in the ionic liquid was observed, showed by nanoparticles located at the interface between aqueous and ionic phases. When 1-octylamine is present during the synthesis, stable nanoparticles of a constant size are obtained. When 1-octylamine is added after the synthesis, aggregation of the ruthenium nanoparticles is observed. In order to explain these phenomena, we have explored the molecular interactions between the different species using (13)C-NMR and DOSY (Diffusional Order Spectroscopy) experiments, mixing calorimetry, surface tension measurements and molecular simulations. We conclude that the behaviour of the ruthenium nanoparticles in [C(1)C(4)Im][NTf(2)] in the presence of 1-octylamine depends on the interaction between the ligand and the nanoparticles in terms of the energetics but also of the structural arrangement of the amine at the nanoparticle's surface.

Influence of ester functional groups on the liquid-phase structure and solvation properties of imidazolium-based ionic liquids.

The incorporation of ester functions in the side chains in 1-alkyl-3-methylimidazolium cations seems to increase the biodegradability of these ionic liquids. We study here how the presence of ester functional groups affects the liquid-state structure (namely, the microphase segregation between polar and nonpolar domains in these ionic liquids) and the way in which the solvation of gases can be understood in these solvents. We use molecular simulation to study the structure of the ionic liquids 3-methyl-1-(pentoxycarbonylmethyl)imidazolium octylsulfate, [C(1)COOC(5)C(1)im][C(8)SO(4)]; and 3-methyl-1-(pentoxycarbonylmethyl)imidazolium bis(trifluoromethylsulfonyl)imide, [C(1)COOC(5)C(1)im][NTf(2)] in the liquid phase and to assess the molecular mechanisms of solvation of carbon dioxide and ethane. The presence of ester functions influences the relative size of the polar and nonpolar domains in the ionic liquids, but does not significantly affect the solvation of gases.

How do physical-chemical parameters influence the catalytic hydrogenation of 1,3-cyclohexadiene in ionic liquids?

The catalytic hydrogenation of 1,3-cyclohexadiene using [Rh(COD)(PPh(3))(2)]NTf(2) (COD = 1,5-cyclooctadiene) was performed in two ionic liquids: 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(4)Im][NTf(2)], and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(1)C(4)Im][NTf(2)]. It is observed that the reaction is twice as fast in [C(1)C(4)Im][NTf(2)] than in [C(1)C(1)C(4)Im][NTf(2)]. To explain the difference in reactivity, molecular interactions and the microscopic structure of ionic liquid +1,3-cyclohexadiene mixtures were studied by NMR and titration calorimetry experiments, and by molecular simulation in the liquid phase. Diffusivity and viscosity measurements allowed the characterization of mass transport in the reaction media. We could conclude that the diffusivity of 1,3-cyclohexadiene is 1.9 times higher in [C(1)C(4)Im][NTf(2)] than in [C(1)C(1)C(4)Im][NTf(2)] and that this difference could explain the lower reactivity observed in [C(1)C(1)C(4)Im][NTf(2)].

Assessing the dispersive and electrostatic components of the cohesive energy of ionic liquids using molecular dynamics simulations and molar refraction data.

Molecular dynamics simulations were used to calculate the density and the cohesive molar internal energy of seventeen different ionic liquids in the liquid phase. The results were correlated with previously reported experimental density and molar refraction data. The link between the dispersive component of the total cohesive energy of the fluid and the corresponding molar refraction was established in an unequivocal way. The results have shown that the two components of the total cohesive energy (dispersive and electrostatic) exhibit strikingly different trends and ratios along different families of ionic liquids, a notion that may help explain their diverse behavior toward different molecular solutes and solvents.