A theoretical and experimental chemist's joint view on hydrogen bonding in ionic liquids and their binary mixtures.

A combined experimental and theoretical approach including quantum chemistry tools and computational simulation techniques can provide a holistic description of the nature of the interactions present in ionic liquid media. The nature of hydrogen bonding in ionic liquids is an especially intriguing aspect, and it is affected by all types of interactions occurring in this media. Overall, these interactions represent a delicate balance of forces that influence the structure and dynamics, and hence the properties of ionic liquids. An understanding of the fundamental principles can be achieved only by a combination of computations and experimental work. In this contribution we show recent results shedding light on the nature of hydrogen bonding, for certain cases the formation of a three-dimensional network of hydrogen bonding, and its dynamics by comparing 1-ethyl-3-methylimidazolium based acetate, chloride and thiocyanate ionic liquids.A particularly interesting case to study hydrogen bonding and other interactions is the investigation of binary mixtures of ionic liquids of the type [cation1][anion1]/[cation1][anion2]. In these mixtures, competing interactions are to be expected. We present both a thorough property meta-analysis of the literature and new data covering a wide range of anions, i.e., mixtures of 1-ethyl-3-methylimidazolium acetate with either trifluoroacetate, tetrafluoroborate, methanesulfonate, or bis(trifluoromethanesulfonyl)imide. In most cases, ideal mixing behavior is found, a surprising result considering the multitude of interactions present. However, ideal mixing behavior allows for the prediction of properties such as density, refractive index, surface tension, and, in most cases, viscosity as function of molar composition. Furthermore, we show that the prediction of properties such as the density of binary ionic liquid mixtures is possible by making use of group contribution methods which were originally developed for less complex non-ionic molecules. Notwithstanding this ideal mixing behavior, several exciting applications are discussed where preferential solvation via hydrogen bonding gives rise to non-additive effects leading to performance improvements. The assessment of the excess properties and (1)H NMR spectroscopic studies provide information on these structural changes and preferential interactions occurring in binary mixtures of ionic liquid, that clearly support the conclusions drawn from the computational studies.

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.

Side chain fluorination and anion effect on the structure of 1-butyl-3-methylimidazolium ionic liquids.

We present a comprehensive molecular dynamics simulation study on 1-butyl-3-methylimidazolium ionic liquids and their fluorinated analogs. The work focused on the effect of fluorination at varying anions. The main findings are that the fluorination of the cations side chain increases overall structuring, especially the aggregation of cation side chain. Furthermore, large and weakly coordinating anions tend to occupy on-top positions of the cation and decrease the aggregation of cation side chains, most likely due to enhanced alkyl-anion interaction.

Effect of dispersion on the structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium thiocyanate.

We present a comprehensive density functional study, using the Perdew-Burke-Ernzerhof (PBE) functional, to elucidate the effect of including or neglecting the dispersion correction on the structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium thiocyanate. We have investigated the structure of the liquid phase and observed that specific interactions between the anions and cations of the ionic liquid were not accurately represented if the dispersion was neglected. The dynamics of the system is more accurately described if the dispersion correction is taken into account and its omission also leads to an incorrect representation of the hydrogen-bonding dynamics. Finally, the power spectrum is predicted and in good agreement with experimental results. Thus, we conclude that it is possible to represent the structure and dynamics of systems containing ionic liquids accurately using ab initio molecular dynamics and a correction for dispersion.

The bulk and the gas phase of 1-ethyl-3-methylimidazolium ethylsulfate: dispersion interaction makes the difference.

Ab initio molecular dynamics simulations were carried out on systems representing the gas and the bulk phase of 1-ethyl-3-methylimidazolium ethylsulfate [C(2)C(1)im][C(2)SO(4)]. The power spectra and cation-anion spatial distribution revealed different interactions of the anion and cation in the bulk phase versus the gas phase. In the bulk phase, all oxygen atoms of the anions are involved and interaction via the rear hydrogen atoms is possible, forming a closer packed system. The alkyl groups of cations and anions governed by dispersion interaction stick together in the bulk but are relatively free in the gas phase.

Proton transfer and polarity changes in ionic liquid-water mixtures: a perspective on hydrogen bonds from ab initio molecular dynamics at the example of 1-ethyl-3-methylimidazolium acetate-water mixtures--part 1.

The ionic liquid 1-ethyl-3-methylimidazolium acetate [C(2)C(1)Im][OAc] shows a great potential to dissolve strongly hydrogen bonded materials, related with the presence of a strong hydrogen bond network in the pure liquid. A first step towards understanding the solvation process is characterising the hydrogen bonding ability of the ionic liquid. The description of hydrogen bonds in ionic liquids is a question under debate, given the complex nature of this media. The purpose of the present article is to rationalise not only the existence of hydrogen bonds in ionic liquids, but also to analyse their influence on the structure of the pure liquid and how the presence of water, an impurity inherent to ionic liquids, affects this type of interaction. We perform an extensive study using ab initio molecular dynamics on the structure of mixtures of the ionic liquid 1-ethyl-3-methylimidazolium acetate with water, at different water contents. Hydrogen bonds are present in the pure liquid, and the presence of water modifies and largely disturbs the hydrogen bond network of the ionic liquid, and also affects the formation of other impurities (carbenes) and the dipole moment of the ions. The use of ab initio molecular dynamics is the recommended tool to explore hydrogen bonding in ionic liquids, as an explicit electronic structure calculation is combined with the study of the condensed phase.

On the ideality of binary mixtures of ionic liquids.

In this work, structural and dynamical properties of the binary mixture of 1-ethyl-3-methyl-imidazolium chloride and 1-ethyl-3-methyl-imidazolium thiocyanate are investigated from ab initio molecular dynamics simulations and compared to the pure ionic liquids. Furthermore, the binary mixture is simulated with two different densities to gain insight into how the selected density affects the different properties. In addition, a simple NMR experiment is carried out to investigate the changes of the chemical shifts of the hydrogen atoms due to the composition of the mixture.

On the density scaling of pVT data and transport properties for molecular and ionic liquids.

In this work, a general equation of state (EOS) recently derived by Grzybowski et al. [Phys. Rev. E 83, 041505 (2011)] is applied to 51 molecular and ionic liquids in order to perform density scaling of pVT data employing the scaling exponent γ(EOS). It is found that the scaling is excellent in most cases examined. γ(EOS) values range from 6.1 for ammonia to 13.3 for the ionic liquid [C(4)C(1)im][BF(4)]. These γ(EOS) values are compared with results recently reported by us [E. R. López, A. S. Pensado, M. J. P. Comuñas, A. A. H. Pádua, J. Fernández, and K. R. Harris, J. Chem. Phys. 134, 144507 (2011)] for the scaling exponent γ obtained for several different transport properties, namely, the viscosity, self-diffusion coefficient, and electrical conductivity. For the majority of the compounds examined, γ(EOS) > γ, but for hexane, heptane, octane, cyclopentane, cyclohexane, CCl(4), dimethyl carbonate, m-xylene, and decalin, γ(EOS) < γ. In addition, we find that the γ(EOS) values are very much higher than those of γ for alcohols, pentaerythritol esters, and ionic liquids. For viscosities and the self-diffusion coefficient-temperature ratio, we have tested the relation linking EOS and dynamic scaling parameters, proposed by Paluch et al. [J. Phys. Chem. Lett. 1, 987-992 (2010)] and Grzybowski et al. [J. Chem. Phys. 133, 161101 (2010); Phys. Rev. E 82, 013501 (2010)], that is, γ = (γ(EOS)/φ) + γ(G), where φ is the stretching parameter of the modified Avramov relation for the density scaling of a transport property, and γ(G) is the Grüneisen constant. This relationship is based on data for structural relaxation times near the glass transition temperature for seven molecular liquids, including glass formers, and a single ionic liquid. For all the compounds examined in our much larger database the ratio (γ(EOS)/φ) is actually higher than γ, with the only exceptions of propylene carbonate and 1-methylnaphthalene. Therefore, it seems the relation proposed by Paluch et al. applies only in certain cases, and is really not generally applicable to liquid transport properties such as viscosities, self-diffusion coefficients or electrical conductivities when examined over broad ranges of temperature and pressure.

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.

Density scaling of the transport properties of molecular and ionic liquids.

Casalini and Roland [Phys. Rev. E 69, 062501 (2004); J. Non-Cryst. Solids 353, 3936 (2007)] and other authors have found that both the dielectric relaxation times and the viscosity, η, of liquids can be expressed solely as functions of the group (TV (γ)), where T is the temperature, V is the molar volume, and γ a state-independent scaling exponent. Here we report scaling exponents γ, for the viscosities of 46 compounds, including 11 ionic liquids. A generalization of this thermodynamic scaling to other transport properties, namely, the self-diffusion coefficients for ionic and molecular liquids and the electrical conductivity for ionic liquids is examined. Scaling exponents, γ, for the electrical conductivities of six ionic liquids for which viscosity data are available, are found to be quite close to those obtained from viscosities. Using the scaling exponents obtained from viscosities it was possible to correlate molar conductivity over broad ranges of temperature and pressure. However, application of the same procedures to the self-diffusion coefficients, D, of six ionic and 13 molecular liquids leads to superpositioning of poorer quality, as the scaling yields different exponents from those obtained with viscosities and, in the case of the ionic liquids, slightly different values for the anion and the cation. This situation can be improved by using the ratio (D∕T), consistent with the Stokes-Einstein relation, yielding γ values closer to those of viscosity.

Interaction between the pi-system of toluene and the imidazolium ring of ionic liquids: a combined NMR and molecular simulation study.

The solute-solvent interactions and the site-site distances between toluene and ionic liquids (ILs) 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BMMIm][NTf2] and 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf2] at various molar ratios were determined by NMR experiments (1D NMR, rotating-frame Overhauser effect spectroscopy (ROESY)) and by molecular simulation using an atomistic force field. The difference in behavior of toluene in these ILs has been related to the presence of H-bonding between the C2-H and the anion in [BMIm][NTf2] generating a stronger association (>20 kJ.mol-1) than in the case of [BMMIm][NTf2]. Consequently, toluene cannot cleave this H-bond in [BMIm][NTf2] which remains in large aggregates of ionic pairs. However, toluene penetrates the less strongly bonded network of [BMMIm][NTf2] and interacts with [BMMIm] cations.

Relationship between viscosity coefficients and volumetric properties using a scaling concept for molecular and ionic liquids.

In this work, a scaling concept based on relaxation theories of the liquid state was combined with a relation previously proposed by the authors to provide a general framework describing the dependency of viscosity on pressure and temperature. Namely, the viscosity-pressure coefficient (partial differentialeta/partial differentialp)T was expressed in terms of a state-independent scaling exponent, gamma. This scaling factor was determined empirically from viscosity versus Tvgamma curves. New equations for the pressure- and temperature-viscosity coefficients were derived, which are of considerable technological interest when searching for appropriate lubricants for elastohydrodynamic lubrication. These relations can be applied over a broad range of thermodynamic conditions. The fluids considered in the present study are linear alkanes, pentaerythritol ester lubricants, polar liquids, associated fluids, and several ionic liquids, compounds selected to represent molecules of different sizes and with diverse intermolecular interactions. The values of the gamma exponent determined for the fluids analyzed in this work range from 1.45 for ethanol to 13 for n-hexane. In general, the pressure-viscosity derivative is well-reproduced with the values obtained for the scaling coefficient. Furthermore, the effects of volume and temperature on viscosity can be quantified from the ratio of the isochoric activation energy to the isobaric activation energy, Ev/Ep. The values of gamma and of the ratio Ev/Ep allow a classification of the compounds according to the effects of density and temperature on the behavior of the viscosity.

Bulk and liquid-vapor interface of pyrrolidinium-based ionic liquids: a molecular simulation study.

Using molecular dynamics simulations, we have studied the structure of three 1-butyl-1-methylpyrrolidinium ionic liquids whose anions are triflate, bis(trifluoromethanesulfonyl)imide, and tris(pentafluoroethyl)trifluorophosphate. The structure of the bulk phase of the three ionic liquids has been interpreted using radial and spatial distribution functions and structure factors that allows us to characterize the morphology of the polar and nonpolar domains present in this family of liquids. The size of the polar regions depends on the anion size, whereas the morphology of the nonpolar domains is anion-independent. Furthermore, the surface ordering properties of the ionic liquids and charge and density profiles were also studied using molecular simulations. The surface tension of the liquid-vapor interfaces of these ionic liquids was also predicted from our molecular simulations. In addition, microscopic structural analysis of orientational ordering at the interface and density profiles along the direction normal to the interface suggest that the alkyl chains of the cation tend to protrude toward the vacuum, and the presence of the interface leads to a strong organization of the liquid phase in the region close to the interface. In the interfacial area, the polar regions of the ionic liquids are more structured than those in the bulk phase, whereas the opposite behavior is observed for the nonpolar regions.

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.

2D or not 2D: structural and charge ordering at the solid-liquid interface of the 1-(2-hydroxyethyl)-3-methylimidazolium tetrafluoroborate ionic liquid.

Molecular dynamics simulations of a 5 nm-thick layer of the ionic liquid 1-(2-hydroxyethyl)-3-methylimidazolium tetrafluoroborate, [(OH)C2C1im][BF4], over silica, alumina and boro-silicate glass substrates have been performed. The structure of the ionic liquid at the solid-liquid interface has been interpreted taking into account the corresponding normal density profiles, lateral interfacial structure, orientational ordering and planar density contours. Comparisons with experimental data suggest that the adsorption and stratification process of ionic liquids over solid substrates can be correctly modeled using a realistic rendition of a non-uniform amorphous substrate such as a glass material.

Using molecular simulation to understand the structure of [C2C1im]+-alkylsulfate ionic liquids: bulk and liquid-vapor interfaces.

Using molecular dynamics simulations we have studied the structure of alkylsulfate-based ionic liquids: 1-ethyl-3-methylimidazolium n-alkylsulfate [C(2)C(1)im][C(n)SO(4)] (n = 2, 4, 6 and 8). The structure of the different ionic liquids have been interpreted taking into account radial and spatial distribution functions, and structure factors, that allowed us to characterize the morphology of the polar and nonpolar domains present in this family of liquids. The size of the nonpolar regions depends linearly on the anion alkyl chain length. Furthermore, properties of the surface of ionic liquids, such as surface tension, ordering, and charge and density profiles, were studied using molecular simulation. We were able to reproduce the experimental values of the surface tension with a maximum deviation of 10%, and it was possible to relate the values of the surface tension with the structure of the liquid-vacuum interfaces. Microscopic structural analysis of orientational ordering at the interface and density profiles along the direction normal to the interface suggest that the alkyl chains of the anions tend to protrude toward the vacuum, and the presence of the interface leads to a strong organization of the liquid phase in the region close to the interface, stronger when the side-chain length of the anions increases.

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.

Molecular dynamics simulations of the liquid surface of the ionic liquid 1-hexyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide: structure and surface tension.

Molecular dynamics simulations of the liquid-vacuum interface of the ionic liquid 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide were performed with an all-atom force field. Structural properties of the interface, such as orientational ordering and density profiles, were calculated. The hexyl side chain of the cation is likely to protrude outward from the surface layer. There is a region with enhanced density from that of the bulk where the cation preferably slants with the imidazolium ring tending to be perpendicular to the interface. The surface tensions are calculated using mechanical and thermodynamic definitions via profiles along the direction normal to the interface. We also discuss the different contributions to the surface tension due to the repulsion-dispersion and electrostatic interactions. The use of local pressure profiles provides an explanation to the systematic problems encountered by several researchers to obtain accurate values of the surface tension at low temperature. Even when macroscopically the system looks in equilibrium, locally this is not accomplished.