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.     

Structure and electronic properties of amino acid ionic liquids

The interactions between eight amino acid based anions and four imidazolium-based cations have been investigated by density functional theory. The electronic and structural properties of the resulting amino acid ionic liquids (AAILs) have been unveiled by means of the atoms in molecules framework. The calculated interaction energy was found to increase in magnitude with decreasing alkyl chain length at imidazolium ring. Moreover, AAILs composed of an amino acid with some functional group such as aromatic ring had decreased interaction energy. Finally, several correlative relationships between glass transition temperature and interaction energy as well as density at bond critical point have been checked for 1-ethyl-3-methylimidazolium based ILs. Although the obtained correlations do not show excellent fits, a preliminary estimation of the glass transition temperature of different AAILs can be achieved by use of their electronic properties.     

Molecular simulation study of some thermophysical and transport properties of triazolium-based ionic liquids

Results of a molecular dynamics study of several triazolium-based ionic liquids are reported. Triazolium cations include 1,2,4-triazolium, 1,2,3-triazolium, 4-amino-1,2,4-triazolium, and 1-methyl-4-amino-1,2,4-triazolium. Each cation was paired with a nitrate or perchlorate anion. These materials are part of a class of ionic compounds that have been synthesized recently but for which little physical property data are available. Properties of the more common ionic liquid, 1-n-butyl-3-methylimidazolium nitrate, are also computed and compared with the properties of the triazolium-based compounds. A molecular mechanics force field was developed for these materials using a mix of ab initio calculations and parameter fitting using the molecular compound 1H-1,2,4-triazole as a basis for the triazolium cations. Liquid-phase properties that were computed include heat capacities, cohesive energy densities, gravimetric densities/molar volumes as a function of temperature and pressure, self-diffusivities, rotational time constants, and various pair correlation functions. In the solid phase, heat capacities and lattice parameters were computed. Of all of these properties, only lattice parameters have been measured experimentally (and only for four of the triazolium compounds). The agreement with the experimental crystal structures was good. When compared with that of the imidazolium-based ionic liquid, the triazolium-based materials have much smaller molar volumes, higher cohesive energy densities, and larger specific heat capacities. They also tend to be less compressible, have a higher gravimetric density, and have faster rotational dynamics but similar translational dynamics.     

Molecular modeling and experimental studies of the thermodynamic and transport properties of pyridinium-based ionic liquids

A combined experimental and molecular dynamics study has been performed on the following pyridinium-based ionic liquids: 1-n-hexyl-3-methylpyridinium bis(trifluoromethanesulfonyl)imide ([hmpy][Tf(2)N]), 1-n-octyl-3-methylpyridinium bis(trifluoromethanesulfonyl)imide ([ompy][Tf(2)N]), and 1-n-hexyl-3,5-dimethylpyridinium bis(trifluoromethanesulfonyl)imide ([hdmpy][Tf(2)N]). Pulsed field gradient nuclear magnetic resonance spectroscopy was used to determine the self-diffusivities of the individual cations and anions as a function of temperature. Experimental self-diffusivities range from 10(-11) to 10(-10) m(2)/s. Activation energies for diffusion are 44-49 kJ/mol. A classical force field was developed for these compounds, and molecular dynamics simulations were performed to compute dynamic as well as thermodynamic properties. Evidence of glassy dynamics was found, preventing accurate determination of self-diffusivities over molecular dynamics time scales. Volumetric properties such as density, isothermal compressibility, and volumetric expansivity agree well with experiment. Simulated heat capacities are within 2% of experimental values.     

Molecular dynamics simulation studies of the influence of imidazolium structure on the properties of imidazolium/azide ionic liquids

Atomistic molecular dynamics simulations were performed on 1-butyl-3-methyl-imidazolium azide [bmim][N(3)], 1-butyl-2,3-dimethylimidazolium azide [bmmim][N(3)], and 1-butynyl-3-methyl-imidazolium azide [bumim][N(3)] ionic liquids. The many-body polarizable APPLE&P force field was augmented with parameters for the azide anion and the bumim cation. Good agreement between the experimentally determined and simulated crystal structure of [bumim][N(3)] as well as the liquid-state density and ionic conductivity of [bmmim][N(3)] were found. Methylation of bmim (yielding bmmim) resulted in dramatic changes in ion structuring in the liquid and slowing of ion motion. Conversely, replacing the butyl group of bmim with the smaller 2-butynyl group resulted in an increase of ion dynamics.     

Liquid-phase structure of dialkylimidazolium ionic liquids from computer simulations

We present a detailed computational study of the structure of ionic liquids based on the imidazolium cation. Both imidazolium-ring stacking and hydrogen bonding behavior are investigated from radial and spatial orientational distribution functions, as well as orientational correlation functions. The alkyl chain size and anion effect on the liquid structure are provided and discussed. Our results support models for liquid organization comparable to those formulated on the basis of experimental observations.     

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.     

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.     

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.     

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.     

Molecular dynamics simulations of polymer transport in nanocomposites

Molecular dynamics simulations on the Kremer-Grest bead-spring model of polymer melts are used to study the effect of spherical nanoparticles on chain diffusion. We find that chain diffusivity is enhanced relative to its bulk value when polymer-particle interactions are repulsive and is reduced when polymer-particle interactions are strongly attractive. In both cases chain diffusivity assumes its bulk value when the chain center of mass is about one radius of gyration R(g) away from the particle surface. This behavior echoes the behavior of polymer melts confined between two flat surfaces, except in the limit of severe confinement where the surface influence on polymer mobility is more pronounced for flat surfaces. A particularly interesting fact is that, even though chain motion is strongly speeded up in the presence of repulsive boundaries, this effect can be reversed by pinning one isolated monomer onto the surface. This result strongly stresses the importance of properly specifying boundary conditions when the near surface dynamics of chains are studied.     

Amphotropic liquid-crystalline behaviour of glycolipids in amino acid ionic liquids

We examined lyotropic liquid-crystalline behaviour of glycolipids (GLs) with a normal alkyl chain or a diacetylene-functionalised alkyl chain in several amino acid ionic liquids (AAILs). It was found that the mixtures of GL and AAIL form various nanosegregated liquid-crystalline phases, such as smectic, bicontinuous cubic and hexagonal columnar phases, depending on the two-component ratio and AAIL species. The observed liquid-crystalline behaviours were summarised as phase diagrams. It is noteworthy that the employment of amino acid anions with superior hydrogen-bonding ability, such as aspartic and glutamic acid anions, gives a phase diagram with a wide liquid-crystalline region. Comparing with a phase diagram obtained for the GL/water mixtures, we gained insights on the similarity/dissimilarity between water and AAILs as self-organisation media of amphiphiles. For the diacetylene-functionalised molecule, UV irradiation was carried out to progress polymerisation. It is of interest that the polymerisation reaction progressed when the glycolipid formed a smectic phase in an AAIL while a reaction progress was not found when it formed a bicontinuous cubic phase in another AAIL. We believe that AAILs have a great potential to be a liquid media not only for amphiphiles but for various functional materials, such as polymers and colloids, to form novel assemblies.     

Quantification of amino acid ionic liquids using liquid chromatography-mass spectrometry

Amino acid ionic liquids (AAILs) containing imidazolium cations and amino acid (AA) anions, were synthesized and applied as task-specific ionic liquids. A sensitive and fast liquid chromatography-mass spectrometry (LC-MS) method was established for the quantitative analysis of 20 AAILs. Using ion pairing-reversed phase liquid chromatography technique, heptafluorobutyric acid was used as ion-pairing reagent to increase the retention of AAILs. Based on the zwitterionity of amino acid, this method was proposed to determine both the cation and the anion of AAILs simultaneously. The limit of detection of this method is down to 1-15ng/mL and the analysis time is less than 15min. According to the analytical data of seven selected AAILs, we found that the content of amino acid anion is always lower than that of butyl methyl imidazolium cation in AAILs. Moreover, the molar ratio of imidazolium cation to amino acid anion is dependent on the chemical property of the amino acid. These results supplied useful information on the interaction of imidazolium cation with acidic, basic, neutral and non-polar amino acids in AAILs.     

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.     

Effect of anionic structure on the phase formation and hydrophobicity of amino acid ionic liquids aqueous two-phase systems

Compared with the conventional ionic liquids, amino acid ionic liquids are more biodegradable and biocompatible, and can enhance stability of biomaterials. In this work, amino acid ionic liquids 1-butyl-3-methylimidazolium L-serine ([C(4)mim][Ser]), 1-butyl-3-methylimidazolium glycine ([C(4)mim][Gly]), 1-butyl-3-methylimidazolium L-alanine ([C(4)mim][Ala]) and 1-butyl-3-methylimidazolium L-leucine ([C(4)mim][Leu]) have been synthesized. These ionic liquids are found to form aqueous two-phase systems (ATPSs) by the salted-out of K(3)PO(4) in aqueous solutions. Phase diagram of the ATPSs and the Gibbs energies of transfer of methylene group from the bottom salt-rich phase to the top ionic liquid-rich phase have been determined at 298.15K and pH 14, and the effect of anionic structure of the ionic liquids on phase formation of the ATPSs and the relative hydrophobicity between the top and the bottom phases are then examined. In order to understand the effect of relative hydrophobicity of the phases in equilibrium in the ATPSs on the extraction/separation capability of biomolecules, the partition coefficients of cytochrome-c (as a model biomolecule) in the ATPSs are measured by spectrophotometry. It is suggested that hydrophobic interactions are mainly responsible for the higher partition coefficients of cytochrome-c in aqueous two-phase systems at pH 14, and the extraction and separation capacity of biomolecules can be improved by the modulation of the relative hydrophobicity of the phases and/or the pH of the system.     

Decoupling charge transport from the structural dynamics in room temperature ionic liquids

Light scattering and dielectric spectroscopy measurements were performed on the room temperature ionic liquid (RTIL) [C4mim][NTf2] in a broad temperature and frequency range. Ionic conductivity was used to estimate self-diffusion of ions, while light scattering was used to study structural relaxation. We demonstrate that the ionic diffusion decouples from the structural relaxation process as the temperature of the sample decreases toward T(g). The strength of the decoupling appears to be significantly lower than that expected for a supercooled liquid of similar fragility. The structural relaxation process in the RTIL follows well the high-temperature mode coupling theory (MCT) scenario. Using the MCT analysis we estimated the dynamic crossover temperature in [C4mim][NTf2] to be T(c) ~ 225 ± 5 K. However, our analysis reveals no sign of the dynamic crossover in the ionic diffusion process.     

Molecular dynamics simulations of ionic liquid-vapour interfaces: effect of cation symmetry on structure at the interface

The influence of alkyl chain symmetry of the imidazolium cation on the structure and properties of the ionic liquid-vapour interface has been addressed through molecular dynamics simulations. The anion chosen is bis(trifluoromethylsulfonyl)imide (NTf(2)). Profiles of number densities, orientation of cations, charge density, electrostatic potential, and surface tension have been obtained. At the interface, both cations and anions were present, and the alkyl chains of the former preferred to orient out into the vapour phase. A large fraction of cations preferred to be oriented with their ring-normal parallel to the surface and alkyl chains perpendicular to it. These orientational preferences are reduced in ionic liquids with symmetric cations. Although the charge densities at the interface were largely negative, an additional small positive charge density has been observed for systems with longer alkyl chains. The electrostatic potential difference developed between the liquid and the vapour phases were positive and decreased with increasing length of the alkyl group. The calculated surface tension of the liquids also decreased with increasing alkyl chain length, in agreement with experiment. The surface tension of an ionic liquid with symmetric cation was marginally higher than that of one with an asymmetric, isomeric cation.     

Molecular dynamics simulations of structure and dynamics of organic molecular crystals

A set of model compounds covering a range of polarity and flexibility have been simulated using GAFF, CHARMM22, OPLS and MM3 force fields to examine how well classical molecular dynamics simulations can reproduce structural and dynamic aspects of organic molecular crystals. Molecular structure, crystal structure and thermal motion, including molecular reorientations and internal rotations, found from the simulations have been compared between force fields and with experimental data. The MM3 force field does not perform well in condensed phase simulations, while GAFF, CHARMM and OPLS perform very similarly. Generally molecular and crystal structure are reproduced well, with a few exceptions. The atomic displacement parameters (ADPs) are mostly underestimated in the simulations with a relative error of up to 70%. Examples of molecular reorientation and internal rotation, observed in the simulations, include in-plane reorientations of benzene, methyl rotations in alanine, decane, isopropylcyclohexane, pyramidal inversion of nitrogen in amino group and rotation of the whole group around the C-N bond. Frequencies of such dynamic processes were calculated, as well as thermodynamic properties for reorientations in benzene and alanine. We conclude that MD simulations can be used for qualitative analysis, while quantitative results should be taken with caution. It is important to compare the outcomes from simulations with as many experimental quantities as available before using them to study or quantify crystal properties not available from experiment.     

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.     

The surface structure of ionic liquids: comparing simulations with x-ray measurements

The surface-normal electron density profile of an ionic liquid, [bmim][PF6], derived from x-ray reflectivity measurements, is compared with two independent molecular-dynamics simulations. It is shown that a meaningful comparison requires a detailed accounting for both thermal and nonthermal surface roughening effects. The former is due to thermally excited capillary waves, and the latter is due to the molecular zero-point motion and form. These quantities influence very significantly, but differently, the simulated and measured density profiles. Stripping off these effects from both types of profiles yields the intrinsic structure factor of the surface. The simulated intrinsic structure factors are found to deviate considerably from the measured one. The introduction of additional ad hoc surface roughness to the simulated profiles greatly reduces the deviation, however, no physical origin for this effect can be identified. The method employed in this study should prove useful for simulation-experiment comparisons of other liquid surfaces, provided they obey capillary-wave theory, as do almost all liquid surfaces studied to date by x-ray reflectivity.