Oscar Cabeza

Molecular Dynamics Simulation of the Structure and Dynamics of Water–1-Alkyl-3-methylimidazolium Ionic Liquid Mixtures

We have performed extensive molecular dynamic simulations to analyze the influence of cation and anion natures, and of water concentration, on the structure and dynamics of water-1-alkyl-3-methylimidazolium ionic liquid mixtures. The dependence on water concentration of the radial distribution functions, coordination numbers, and hydrogen bonding degree between the different species has been systematically analyzed for different lengths of the cation alkyl chain (alkyl = ethyl, butyl, hexyl, and octyl) and several counterions. These include two halogens of different sizes and positions in Hoffmeister series, Cl(-) and Br(-), and the highly hydrophobic inorganic anion PF(6)(-) throughout its whole solubility regime. The formation of water clusters in the mixture has been verified, and the influences of both anion hydrophobicity and cation chain length on the structure and size of these clusters have been analyzed. The water cluster size is shown to be relatively independent of the cation chain length, but strongly dependent on the hydrophobicity of the anion, which also determines critically the network formation of water and therefore the miscibility of the ionic liquid. The greater influence of the anion relative to the cation one is seen to be reflected in all the analyzed physical properties. Finally, single-particle dynamics in IL-water mixtures is considered, obtaining the self-diffusion coefficients and the velocity autocorrelation functions of water molecules in the mixture, and analyzing the effect of cation, anion, and water concentration on the duration of the ballistic regime and on the time of transition to the diffusive regime. Complex non-Markovian behavior was detected at intermediate times within an interval progressively shorter as water concentration increases.

Effect of temperature and cationic chain length on the physical properties of ammonium nitrate-based protic ionic liquids.

We report a systematic study of the effect of the cationic chain length and degree of hydrogen bonding on several equilibrium and transport properties of the first members of the alkylammonium nitrate protic ionic liquids (PILs) family (ethylammonium, propylammonium, and butylammonium nitrate) in the temperature range between 10 and 40 °C. These properties were observed by means of several experimental techniques, including density, surface tension, refractometry, viscosimetry, and conductimetry. The dilatation coefficients and compressibilities, as well as the Rao coefficients, were calculated, and an increase of these magnitudes with alkyl chain length was detected. Moreover, the surface entropies and enthalpies of the studied PILs were analyzed, and the temperature dependence of the surface tension was observed to be describable by means of a harmonic oscillator model with surface energies and critical temperatures that are increasing functions of the cationic chain length. Moreover, the refractive indexes were measured and the thermo-optic coefficient and Abbe numbers were calculated, and the contribution of the electrostrictive part seemed to dominate the temperature dependence of the electric polarization. The electric conductivity and the viscosity were measured and the influence of the degree of hydrogen bonding in the supercooled liquid region analyzed. Hysteresis loops were detected in freezing-melting cycles and the effect of the length of the alkyl chain of the cation on the size of the loop analyzed, showing that longer chains lead to a narrowing of the supercooled region. The temperature dependence of the conductivity was studied in the Vogel-Fulcher-Tamman (VFT) framework and the fragility indices, the effective activation energies, and the Vogel temperatures obtained. A high-temperature Arrhenius analysis was also performed, and the activation energies of conductivity and viscosity were calculated, showing that these transport processes are governed by two distinct mechanisms. The exponents of the fractional Walden rule for the different compounds were obtained. Finally, the ionicities and fragilities of the studied PILs were analyzed, proving that all the studied PILs are subionic and fragile liquids, with propylammonium nitrate showing the lowest fragility and the greater ionicity of all the studied compounds.

MD simulations of the formation of stable clusters in mixtures of alkaline salts and imidazolium-based ionic liquids.

Structural and dynamical properties of room-temperature ionic liquids containing the cation 1-butyl-3-methylimidazolium ([BMIM](+)) and three different anions (hexafluorophosphate, [PF6](-), tetrafluoroborate, [BF4](-), and bis(trifluoromethylsulfonyl)imide, [NTf2](-)) doped with several molar fractions of lithium salts with a common anion at 298.15 K and 1 atm were investigated by means of molecular dynamics simulations. The effect of the size of the salt cation was also analyzed by comparing these results with those for mixtures of [BMIM][PF6] with NaPF6. Lithium/sodium solvation and ionic mobilities were analyzed via the study of radial distribution functions, coordination numbers, cage autocorrelation functions, mean-square displacements (including the analysis of both ballistic and diffusive regimes), self-diffusion coefficients of all the ionic species, velocity and current autocorrelation functions, and ionic conductivity in all the ionic liquid/salt systems. We found that lithium and sodium cations are strongly coordinated in two different positions with the anion present in the mixture. Moreover, [Li](+) and [Na](+) cations were found to form bonded-like, long-lived aggregates with the anions in their first solvation shell, which act as very stable kinetic entities within which a marked rattling motion of salt ions takes place. With very long MD simulation runs, this phenomenon is proved to be on the basis of the decrease of self-diffusion coefficients and ionic conductivities previously reported in experimental and computational results.

Solvation of lithium salts in protic ionic liquids: a molecular dynamics study.

The structure of solutions of lithium nitrate in a protic ionic liquid with a common anion, ethylammonium nitrate, at room temperature is investigated by means of molecular dynamics simulations. Several structural properties, such as density, radial distribution functions, hydrogen bonds, spatial distribution functions, and coordination numbers, are analyzed in order to get a picture of the solvation of lithium cations in this hydrogen-bonded, amphiphilically nanostructured environment. The results reveal that the ionic liquid mainly retains its structure upon salt addition, the interaction between the ammonium group of the cation and the nitrate anion being only slightly perturbed by the addition of the salt. Lithium cations are solvated by embedding them in the polar nanodomains of the solution formed by the anions, where they coordinate with the latter in a solid-like fashion reminiscent of a pseudolattice structure. Furthermore, it is shown that the average coordination number of [Li](+) with the anions is 4, nitrate coordinating [Li](+) in both monodentate and bidentate ways, and that in the second coordination layer both ethylammonium cations and other lithiums are also found. Additionally, the rattling motion of lithium ions inside the cages formed by their neighboring anions, indicative of the so-called caging effect, is confirmed by the analysis of the [Li](+) velocity autocorrelation functions. The overall picture indicates that the solvation of [Li](+) cations in this amphiphilically nanostructured environment takes place by means of a sort of inhomogeneous nanostructural solvation, which we could refer to as nanostructured solvation, and which could be a universal solvation mechanism in ionic liquids.

On the Interactions between Amino Acids and Ionic Liquids in Aqueous Media

The understanding of the molecular-level interactions between biomolecules and ionic liquids (ILs) in aqueous media is crucial for the optimization of a number of relevant biotechnological processes. In this work, the influence of a series of amino acids on the liquid-liquid equilibria between 1-butyl-3-methylimidazolium tricyanomethane and water was studied to evaluate the preferential interactions between these three compounds. The solubility effects observed are dependent on the polarity, size, and charge distribution of the amino acid side chains and are explained in terms of a refined version of the model proposed earlier (Freire et al. J. Phys. Chem. B 2009, 113, 202; Tome et al. J. Phys. Chem. B 2009, 113, 2815) for ion specific effects on aqueous solutions of imidazolium-based ILs. Although acting through different mechanisms, salting-in and salting-out phenomena possess a common basis which is the competition between water-amino acid side chain, IL-amino acid side chain, and water-IL interactions. The delicate balance between these interactions is dependent on the relative affinities of the biomolecules to water molecules or to IL cation and anion and determines the trend and magnitude of the solubility effect observed.

Surface Tension of Binary Mixtures of 1-Alkyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide Ionic Liquids: Experimental Measurements and Soft-SAFT Modeling

Ionic liquids have attracted a large amount of interest in the past few years. One approach to better understand their peculiar nature and characteristics is through the analysis of their surface properties. Some research has provided novel information on the organization of pure ionic liquids at the vapor-liquid interface; yet, a systematic study on the surface properties of mixtures of ionic liquids and their organization at the surface has not previously been carried out in the literature. This work reports, for the first time, a comprehensive analysis of the surface organization of mixtures of ionic liquids constituted by 1-alkyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquids, [C(n)mim][NTf(2)]. The surface tension of mixtures composed of [C(4)mim][NTf(2)] + [C(n)mim][NTf(2)] (n = 1, 2, 5, 6, 8, and 10) was experimentally determined, at 298.2 K and atmospheric pressure, in the whole composition range. From the experimental data, the surface tension deviations and the relative Gibbs adsorption isotherms were estimated showing how the surface composition of an ionic liquid mixture differs from that of the liquid bulk and that the surface is enriched by the ionic liquid with the longest alkyl chain length. Finally, the soft-SAFT equation of state coupled with the density gradient theory (DGT) was used, for the first time, to successfully reproduce the surface tension experimental data of binary mixtures of ionic liquids using a molecular-based approach. In addition, the DGT was used to compute the density profiles of the two components across the interface, confirming the experimental results for the components distribution at the bulk and at the vapor-liquid interface.

Mixtures of protic ionic liquids and molecular cosolvents: a molecular dynamics simulation.

In this work, the effect of molecular cosolvents (water, ethanol, and methanol) on the structure of mixtures of these compounds with a protic ionic liquid (ethylammonium nitrate) is analyzed by means of classical molecular dynamics simulations. Included are as-yet-unreported measurements of the densities of these mixtures, used to test our parameterized potential. The evolution of the structure of the mixtures throughout the concentration range is reported by means of the calculation of coordination numbers and the fraction of hydrogen bonds in the system, together with radial and spatial distribution functions for the various molecular species and molecular ions in the mixture. The overall picture indicates a homogeneous mixing process of added cosolvent molecules, which progressively accommodate themselves in the network of hydrogen bonds of the protic ionic liquid, contrarily to what has been reported for their aprotic counterparts. Moreover, no water clustering similar to that in aprotic mixtures is detected in protic aqueous mixtures, but a somehow abrupt replacing of [NO3](-) anions in the first hydration shell of the polar heads of the ionic liquid cations is registered around 60% water molar concentration. The spatial distribution functions of water and alcohols differ in the coordination type, since water coordinates with [NO3](-) in a bidentate fashion in the equatorial plane of the anion, while alcohols do it in a monodentate fashion, competing for the oxygen atoms of the anion. Finally, the collision times of the different cosolvent molecules are also reported by calculating their velocity autocorrelation functions, and a caging effect is observed for water molecules but not in alcohol mixtures.

Molecular dynamics simulations of the structural and thermodynamic properties of imidazolium-based ionic liquid mixtures.

In this work, extensive molecular dynamics simulations of mixtures of alcohols of several chain lengths (methanol and ethanol) with the ionic liquids (ILs) composed of the cation 1-hexyl-3-methylimidazolium and several anions of different hydrophobicity degrees (Cl(-), BF(4)(-), PF(6)(-)) are reported. We analyze the influence of the nature of the anion, the length of the molecular chain of the alcohol, and the alcohol concentration on the thermodynamic and structural properties of the mixtures. Densities, excess molar volumes, total and partial radial distribution functions, coordination numbers, and hydrogen bond degrees are reported and analyzed for mixtures of the ILs with methanol and ethanol. The aggregation process is shown to be highly dependent on the nature of the anion and the size of the alcohol, since alcohol molecules tend to interact predominantly with the anionic part of the IL, especially in mixtures of the halogenated IL with methanol. Particularly, our results suggest that the formation of an apolar network similar to that previously reported in mixtures of ILs with water does not take place in mixtures with alcohol when the chloride anion is present, the alcohol molecules being instead homogeneously distributed in the polar network of IL. Moreover, the alcohol clusters formed in mixtures of [HMIM][PF(6)] with alcohol were found to have a smaller size than in mixtures with water. Additionally, we provide a semiquantitative analysis of the dependence of the hydrogen bonding degree of the mixtures on the alcohol concentration.

Surfactant self-assembly nanostructures in protic ionic liquids.

The existence and properties of mesoscopic self-assembly structures formed by surfactants in protic ionic liquid solutions are reported. Micellar aggregates of n-alkyltrimethylammonium (n = 10, 12, 14, 16) chlorides and bromides and of n-alkylpyridinium (n = 12, 16) chlorides in ethylammonium nitrate and propylammonium nitrate were observed by means of several experimental techniques, including surface tension, transmission electron micrography, dynamic light scattering, and potentiometry using surfactant-selective electrodes. The effect of the alkyl chain length of both solute and solvent molecules on the critical micelle concentration is discussed, and a Stauff-Klevens law is seen to apply to surfactant solutions in both protic ionic liquids. The counterion role is also a matter of study in the case of alkyltrimethylammonium-based surfactants, and the presently reported evidence suggests that the place of the surfactant counterion in the Hoffmeisters series could determine its effect on micellization in IL solution. The size distribution of the aggregates is also analyzed together with the Gibbs free energies of micellization and the minimum surface area per monomer in all of the studied cases. All of the hereby reported evidence suggests that the negative entropic contribution arising from the release of the solvent layer upon micellization is also the driving force of conventional surfactant self-association in protic ionic liquids.

Electrical conductivity of seven binary systems containing 1-ethyl-3-methyl imidazolium alkyl sulfate ionic liquids with water or ethanol at four temperatures.

We present experimental measurements of specific electrical (or ionic) conductivity of seven binary systems of 1-ethyl-3-methyl imidazolium alkyl sulfate (EMIM-C(n)S) with water or ethanol. Electrical conductivity was measured at 298.15 K in all ranges of concentrations and selected mixtures also at 288.15, 308.15, and 318.15 K. The alkyl chains of the anions used are ethyl (EMIM-ES), butyl (EMIM-BS), hexyl (EMIM-HS), and, only for mixtures with ethanol, octyl (EMIM-OS). Let us note that the four ionic liquids (ILs) measured are miscible in water and ethanol at those temperatures and atmospheric pressure in all ranges of concentrations, but EMIM-OS jellifies for a given range of concentration with water. We compare the measured data in terms of the alkyl chain length and solvent nature. Data are compared with previously scarce results for these same systems and also for other aqueous and ethanol mixtures with ILs. In addition, we verify that our data fit the universal theoretical expression with no fitting parameters given by the pseudolattice-based Bahe-Varela model, except for IL concentrated mixtures. To fit well all ranges of concentrations, we add to the original equation two phenomenological terms with one fitting parameter each. Finally, we calculate the molar conductivity and fit it successfully with an expression derived from Onsager theory.