Aqueous interfaces with hydrophobic room-temperature ionic liquids: a molecular dynamics study

We report a molecular dynamics study of the interface between water and (macroscopically) water-immiscible room-temperature ionic liquids "ILs", composed of PF6(-) anions and butyl- versus octyl-substituted methylimidazolium+ cations (noted BMI+ and OMI+). Because the parameters used to simulate the pure ILs were found to exaggerate the water/IL mixing, they have been modified by scaling down the atomic charges, leading to better agreement with the experiment. The comparison of [OMI][PF6] versus [BMI][PF6] ILs demonstrates the importance of the N-alkyl substituent on the extent of solvent mixing and on the nature of the interface. With the most hydrophobic [OMI][PF6] liquid, the "bulk" IL phase is dryer than with the [BMI][PF6] liquid. At the interface, the OMI+ cations retain direct contacts with the bulk IL, whereas the more hydrophilic PF6(-) anions gradually dilute in the local water micro-environment and are thus isolated from the "bulk" IL. The interfacial OMI+ cations are ordered with their imidazolium moiety pointing toward the aqueous side and their octyl chains toward the IL side of the interface. With the [BMI][PF6] liquid, the system gradually evolves from an IL-rich to a water-rich medium, leading to an ill-defined interfacial domain with high intersolvent mixing. As a result, the BMI+ cations are isotropically oriented "at the interface". Because the imidazolium cations are more hydrophobic than the PF6(-) anions, the charge distribution at the interface is heterogeneous, leading to a positive electrostatic potential at the interface with the two studied ILs. Mixing-demixing simulations on [BMI][PF6]/water mixtures are also reported, comparing Ewald versus reaction field treatments of electrostatics. Phase separation is very slow (at least 30 ns), in marked contrast with mixtures involving classical organic liquids, which separate in less than 0.5 ns at the microscopic level. The results allow us to better understand the specificity of the aqueous interfaces with hydrophobic ionic liquids, compared with classical organic solvents, which has important implications as far as the mechanism of liquid-liquid ion extraction is concerned.     

Molecular dynamics simulation of dual amino-functionalized imidazolium-based ionic liquids

In this work, the all-atom (AA) force fields were set up for three kinds of dual amino-functionalized imidazolium-based ionic liquids (ILs), composed by cations with different alkyl chain length and amino acid anion [Gly]⁻. The force field was based on our previous work and the default parameters were developed in this study. Molecular dynamics simulations were performed. Validation was carried out by comparing simulation densities with experimental data, and good agreement was obtained. Molar volume and heat capacity at constant pressure were predicted. Mean square displacements for these ILs were computed and these ILs were proved to move very slowly. It may be caused by hydrogen-bonded network between ions and the terminal azyl. To depict the microscopic structures of the ILs, many types of radial distribution functions were investigated. It is interesting to find that not only the cation and anion, but also the anions themselves will form hydrogen bonds.     

A comparison of ether- and alkyl-derivatized imidazolium-based room-temperature ionic liquids: a molecular dynamics simulation study

Molecular dynamics simulations of ether-derivatized imidazolium-based room-temperature ionic liquids (EDI-RTILs), [C(5)O(2)mim][TFSI] and [C(5)O(2)mim][BF(4)], have been performed and compared with simulations of alkyl-derivatized analogues (ADI-RTILs). Simulations yield RTIL densities, self-diffusion coefficients and viscosity in excellent agreement with experimental data. Simulations reveal that structure in the EDI-RTILs, quantified by the extent of nanoscale segregation of tails as well as cation-ion and cation-cation correlations, is reduced compared to that observed in the ADI-RTILs. Significant correlation between ether tail oxygen atoms and imidazolium ring hydrogen atoms was observed in the EDI-RTILs. This correlation is primarily intramolecular in origin but has a significant intermolecular component. Competition of ether oxygen atoms with oxygen atoms of TFSI(-) or fluorine atoms of BF(4)(-) for coordination of the ring hydrogen atoms was found to reduce the extent of cation-anion correlation in the EDI-RTILs compared to the ADI-RTILs. The reduction in intermolecular correlation, particularly tail-tail segregation, as well as weakening of cation-anion specific interactions due to the ether tail, may account for the faster dynamics observed in the EDI-RTILs compared to ADI-RTILs.     

Ultrafast solvation dynamics of coumarin 153 in imidazolium-based ionic liquids

The dynamic Stokes shift of coumarin 153 has been measured in two room-temperature ionic liquids, 1-(3-cyanopropyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide and 1-propyl-3-methylimidazolium tetrafluoroborate, using the fluorescence up-conversion technique with a 230 fs instrumental response function. A component of about 10-15% of the total solvation shift is found to take place on an ultrafast time scale < 10 ps. The amplitude of this component is substantially less than assumed previously by other authors. The origin of the difference in findings could be partly due to chromophore-internal conformational changes on the ultrafast time scale, superimposed to solvation-relaxation, or due to conformational changes of the chromophore ground state in polar and apolar environments. First three-pulse photon-echo peak-shift experiments on indocyanine green in room-temperature ionic liquids and in ethanol indicate a difference in the inertial component of the early solvent relaxation of <100 fs.     

Ab initio molecular dynamics simulations of a binary system of ionic liquids

This work presents first insights into the structural properties of a binary mixture of ionic liquids from the perspective of ab initio molecular dynamics simulations. Simulations were carried out for a one-to-one mixture of 1-ethyl-3-methyl-imidazolium thiocyanate and 1-ethyl-3-methyl-imidazolium chloride and compared to pure 1-ethyl-3-methyl-imidazolium thiocyanate.     

Nanoclusters in ionic liquids: evidence for N-heterocyclic carbene formation from imidazolium-based ionic liquids detected by (2)H NMR

The mystery of how 1,3-substituted imidazolium-based ionic liquids (ILs) can provide high stabilization for transition-metal(0) nanoclusters, that is, in the absence of the usual strongly coordinating anions, has been probed. 2H NMR product and kinetic studies of 1,3-substituted imidazolium ILs under D2 reveal that nanocluster-catalyzed H/D exchange occurs at the 2- (as well as at the 4-, 5-, and 8-) C-H positions of the imidazolium cation. The results (i) provide compelling evidence that N-heterocyclic carbene formation and ligation of nanoclusters is occurring in ILs; and (ii) argue that N-heterocyclic carbenes merit further investigation as heretofore unappreciated stabilizers of transition-metal nanoclusters.     

Thermodynamic studies of ionic interactions in aqueous solutions of imidazolium-based ionic liquids [Emim][Br] and [Bmim][Cl

Experimental measurements of density at different temperatures ranging from 293.15 to 313.15 K, the speed of sound and osmotic coefficients at 298.15 K for aqueous solution of 1-ethyl-3-methylimidazolium bromide ([Emim][Br]), and osmotic coefficients at 298.15 K for aqueous solutions of 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) in the dilute concentration region are taken. The data are used to obtain compressibilities, expansivity, apparent and limiting molar properties, internal pressure, activity, and activity coefficients for [Emim][Br] in aqueous solutions. Experimental activity coefficient data are compared with that obtained from Debye-Hückel and Pitzer models. The activity data are further used to obtain the hydration number and the osmotic second virial coefficients of ionic liquids. Partial molar entropies of [Bmim][Cl] are also obtained using the free-energy and enthalpy data. The distance of the closest approach of ions is estimated using the activity data for ILs in aqueous solutions and is compared with that of X-ray data analysis in the solid phase. The measured data show that the concentration dependence for aqueous solutions of [Emim][Br] can be accounted for in terms of the hydrophobic hydration of ions and that this IL exhibits Coulombic interactions as well as hydrophobic hydration for both the cations and anions. The small hydration numbers for the studied ILs indicate that the low charge density of cations and their hydrophobic nature is responsible for the formation of the water-structure-enforced ion pairs.     

Air-liquid interfaces of aqueous solutions containing ammonium and sulfate: spectroscopic and molecular dynamics studies

Investigations of the air-liquid interface of aqueous salt solutions containing ammonium (NH(4)(+)) and sulfate (SO(4)(2-)) ions were carried out using molecular dynamics simulations and vibrational sum frequency generation spectroscopy. The molecular dynamics simulations show that the predominant effect of SO(4)(2-) ions, which are strongly repelled from the surface, is to increase the thickness of the interfacial region. The vibrational spectra reported are in the O-H stretching region of liquid water. Isotropic Raman and ATR-FTIR (attenuated total reflection Fourier transform infrared) spectroscopies were used to study the effect of ammonium and sulfate ions on the bulk structure of water, whereas surface sum frequency generation spectroscopy was used to study the effect of these ions on the interfacial structure of water. Analysis of the interfacial and bulk vibrational spectra reveal that aqueous solutions containing SO(4)(2-) perturb the interfacial water structure differently than the bulk and, consistent with the molecular dynamics simulations, reveal an increase in the thickness of the interfacial region.     

On simulations of complex interfaces: molecular dynamics simulations of stationary phases

Methodological considerations for molecular dynamics simulations of complex interfaces are presented in this article. A slab geometry is examined in the context of stationary phases where selectivity occurs predominantly in pores within silica beads. Specifically, we examine the Whelk-O1 interface with n-hexane/2-propanol, the TMA-(Pro)(2)-N(CH(3))-tether interface with n-hexane/2-propanol, and the C(18)H(37)Si interface with water/methanol. The following methodological issues are considered in detail: The assessment of solvent density within the confined region and excluded volume of the interface; the structural equilibration of surface-bound moieties; solvent equilibration for binary mixtures; surface size effects, and periodic boundary conditions; the treatment of electrostatic interactions; and the impact of pore size.     

Transport properties of room-temperature ionic liquids from classical molecular dynamics

Room-temperature ionic liquids (RTILs) have attracted much attention in the scientific community in the past decade due their novel and highly customizable properties. Nonetheless, their high viscosities pose serious limitations to the use of RTILs in practical applications. To elucidate some of the physical aspects behind transport properties of RTILs, extensive classical molecular dynamics calculations are reported. Here, in particular, bulk viscosities and ionic conductivities of butyl-methyl-imidazole based RTILs are presented over a wide range of temperatures. The dependence of the properties of the liquids on simulation parameters, e.g., system-size effects or the choice of the interaction potential, is analyzed in detail.     

Picosecond time-resolved fluorescence study on solute-solvent interaction of 2-aminoquinoline in room-temperature ionic liquids: aromaticity of imidazolium-based ionic liquids

Time-resolved fluorescence spectra and fluorescence anisotropy decay of 2-aminoquinoline (2AQ) have been measured in eight room-temperature ionic liquids, including five imidazolium-based aromatic ionic liquids and three nonaromatic ionic liquids. The same experiments have also been carried out in several ordinary molecular liquids for comparison. The observed time-resolved fluorescence spectra indicate the formation of pi-pi aromatic complexes of 2AQ in some of the aromatic ionic liquids but not in the nonaromatic ionic liquids. The fluorescence anisotropy decay data show unusually slow rotational diffusion of 2AQ in the aromatic ionic liquids, suggesting the formation of solute-solvent complexes. The probe 2AQ molecule is likely to be incorporated in the possible local structure of ionic liquids, and hence the anisotropy decays only through the rotation of the whole local structure, making the apparent rotational diffusion of 2AQ slow. The rotational diffusion time decreases rapidly by adding a small amount of acetonitrile to the solution. This observation is interpreted in terms of the local structure formation in the aromatic ionic liquids and its destruction by acetonitrile. No unusual behavior upon addition of acetonitrile has been found for the nonaromatic ionic liquids. It is argued that the aromaticity of the imidazolium cation plays a key role in the local structure formation in imidazolium-based ionic liquids.     

Thermochemistry of imidazolium-based ionic liquids: experiment and first-principles calculations

In this work the molar enthalpy of formation of the ionic liquid 1-ethyl-3-methylimidazolium dicyanoamide in the gaseous phase [C(2)MIM][N(CN)(2)] was measured by means of combustion calorimetry and enthalpy of vaporization using transpiration. Available, but scarce, primary experimental results on enthalpies of formation of imidazolium based ionic liquids with the cation [C(n)MIM] (where n = 2 and 4) and anions [N(CN)(2)], [NO(3)] and [NTf(2)] were collected and checked for consistency using a group additivity procedure. First-principles calculations of the enthalpies of formation in the gaseous phase for the ionic liquids with the common cation [C(n)MIM] (where n = 2 and 4) and with the anions [N(CN)(2)], [NO(3)], [NTf(2)], [Cl], [BF(4)] and [PF(6)] have been performed using the G3MP2 theory. It has been established that the gaseous phase enthalpies of formation of these ionic liquids obey the group additivity rules.     

Aqueous solutions of ionic liquids: study of the solution/vapor interface using molecular dynamics simulations

We performed a detailed molecular dynamics study of the interfacial structure of aqueous solutions of 1-butyl-3-methylimidazolium tetrafluoroborate in order to explain the anomalous dependence of the surface tension on concentration. At low concentrations the surface tension decreases with concentration. At higher concentrations the surface becomes saturated; a plateau is observed in simulations with a non-polarizable force field while a possible increase is detected in simulations with a polarizable force field. The structure is characterized by a surplus of cations at the surface (with hydrophobic butyl chains pointing toward vacuum) which at low concentrations is only partly compensated by the anions because of asymmetric solvation. A more hydrophobic 1-butyl-3-methylimidazolium hexafluorophosphate is also simulated for comparison.     

Ab initio molecular dynamics simulation of ionic liquids

Ab initio Car-Parinnello molecular dynamics is used to simulate the structure and the dynamics of 1-butyl-3-methylimidazolium iodide ([bmim]I) ionic liquid at 300 K. Site-site pair correlation functions reveal that the anion has a strong interaction with any three C-H's of the imidazolium ring. The ring bends over and wraps around the anion such that the two nitrogen atoms take a distance to the anion. Electron donating butyl group contributes the electronic polarization in addition to geometrical (out-of-plane) polarization of the ring due to the liquid environment. This facilitates bending of the ring along the axis passing through nitrogen atoms. The average bending angle depends largely on the alkyl chain length and slightly on the anion type. Redistribution of electron density over the ring caused by the electron donating alkyl group provides additional independent evidence to the instability of lattice structure, hence the low melting point of the ionic liquid. Simulated viscosity and diffusion coefficients of [bmim]I are in quite agreement with the experiments.     

Molecular dynamics simulations of the surface tension of ionic liquids

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.     

Molecular dynamics simulations of nanoparticle self-assembly at ionic liquid-water and ionic liquid-oil interfaces

We have studied the self-assembly of hydrophobic nanoparticles at ionic liquid (IL)-water and IL-oil (hexane) interfaces using molecular dynamics (MD) simulations. For the 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF(6)])/water system, the nanoparticles rapidly approached the IL-water interface and equilibrated more into the IL phase although they were initially in the water phase. In contrast, when the nanoparticles were dispersed in the hexane phase, they slowly approached the IL-hexane interface but remained primarily in the hexane phase. Consequently, the IL-hexane interface was rather undisturbed by the nanoparticles whereas the IL-water interface changed significantly in width and morphology to accommodate the presence of the nanoparticles. The equilibrium positions of the nanoparticles were also supported and explained by potential of mean force (PMF) calculations. Interesting ordering and charge distributions were observed at the IL-liquid interfaces. At the IL-hexane interface, the [BMIM] cations preferentially oriented themselves so that they were immersed more in the hexane phase and packed efficiently to reduce steric hindrance. The ordering likely contributed to a heightened IL density and a slightly positive charge at the IL-hexane interface. In contrast, the cations at the IL-water interface were oriented isotropically unless in the presence of nanoparticles, where the cations aligned across the nanoparticle surfaces.     

Electrocatalytic reduction of CO2 in neat and water-containing imidazolium-based ionic liquids

Energetically efficient electrochemical reduction of CO₂ would offer the possibility of storing electricity from renewables in the form of fuels and other valuable chemicals. It may also help mitigate the increase of atmospheric CO₂ associated with global warming. However, the process suffers from a low energy efficiency because of the large overpotentials required. In aqueous electrolytes, the competing hydrogen evolution reaction also decreases the faradaic efficiency (which contributes to the low energy efficiency of the process). Recent claims of high faradaic efficiency and low overpotentials for the reduction of CO₂ in room-temperature ionic liquids (RTILs) and RTIL–water mixtures have spurred considerable research. Here, we offer a critical review of those claims and of recent work aimed at understanding the details of this important reaction in these nonconventional electrolytes.     

Effects of methylation at position 2 of cation ring on rotational dynamics of imidazolium-based ionic liquids investigated by NMR spectroscopy: [C4mim]Br vs [C4C1mim]Br

We investigated the rotational dynamics of two imidazolium-based ionic liquids, 1-butyl-3-methylimidazolium bromide ([C(4)mim]Br) and 1-butyl-2,3-dimethylimidazolium bromide ([C(4)C(1)mim]Br), to reveal the effects of methylation at position 2 of the imidazolium ring (C(2) methylation). The rotational correlation time (τ(local)) for each carbon in the cations is derived from the spin-lattice relaxation time of (13)C nuclear magnetic resonance. The τ(local) results obtained here provide three principle insights into the rotational dynamics of ionic liquids. First, all τ(local) values for [C(4)C(1)mim]Br are greater than those for [C(4)mim]Br owing to a viscosity increase due to C(2) methylation. Second, the rate of change in τ(local) on C(2) methylation differs among the carbons in the cation, which indicates that each carbon has a different microviscosity. Third, the τ(local) increase in the (13)C at the root of the butyl group on C(2) methylation is very small compared to both intuitive prediction and the results from quantum chemical calculations. This indicates that the motion of the butyl group root in [C(4)C(1)mim]Br is not significantly inhibited by the methyl group at the position 2 of the imidazolium ring. The finding provides conclusive information on the origin of the increases in the melting point on C(2) methylation. Hunt previously found through calculation that decreases in entropy are caused by two factors, namely, reductions in the rotational mobility of the butyl group and in the number of stable anion interaction sites with C(2) methylation, resulting in an increase in melting point and viscosity. Our finding experimentally illustrates that the origin of the increases in melting point is not the inhibition of butyl group motion and that the reduction in stable anion interaction sites plays a major role in the increases. Additionally, it is suggested that the viscosity increase on C(2) methylation can be interpreted in the same manner.     

Solubility of CO2 in imidazolium-based tetrafluoroborate ionic liquids

The solubility of CO₂ in imidazolium ionic liquids (ILs), 1-butyl-3-methyl imidazolium tetrafluoroborate ([bmim][BF₄]), 1-hexyl-3-methyl imidazolium tetrafluoroborate ([hmim][BF₄]) and 1-octyl-3-methyl imidazolium tetrtafluoroborate ([omim][BF₄]) was determined at 305-25 K and pressures from 1 to 9 MPa. The influence of chain length of alkyl substituents on the imidazolium cation on the solubility of CO₂ was investigated. The differences in solubility with chain length are in the sequence [omim][BF₄] > [hmim][BF₄] > [bmim][BF₄]. The solubility data were correlated by the extended Henry's law, and enthalpy, Gibbs free energy and entropy changes were obtained.     

Molecular dynamics simulations of charged nanoparticle self-assembly at ionic liquid-water and ionic liquid-oil interfaces

Nanoparticle self-assembly at liquid-liquid interfaces can be significantly affected by the individual nanoparticle charges. This is particularly true at ionic liquid (IL) based interfaces, where Coulombic forces play a major role. Employing 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF(6)]) as a model IL, we have studied the self-assembly of hydrophobic nanoparticles with different surface charges at the IL/water and IL/oil (hexane) interfaces using molecular dynamics simulations. In the IL/water system, the nanoparticles were initially dispersed in the water phase but quickly equilibrated at the interface, somewhat in favor of the IL phase. This preference was lessened with increased nanoparticle charge. In the IL/hexane system, all charged nanoparticles interacted with the IL to some extent, whereas the uncharged nanoparticles remained primarily in the hexane phase. Potential of mean force calculations supported the observations from the equilibrium studies and provided new insights into the interactions of the nanoparticles and ionic liquid based interfaces.