The properties of ionic liquids depend on the chemical structure of the constituent ions. An important difference between molten inorganic salts and room temperature ionic liquids (RTILs) is that in RTILs the charge is frequently not located at the center of mass. This paper describes a molecular dynamics investigation of the influence of charge location on the structure and transport properties of ionic liquids. The model considered consists of univalent spherical ions with the cation charge moved away from its center of mass. It is shown that the charge location has an important influence on the liquid properties. As the charge is moved off center, the electrical conductivity initially increases, and the shear viscosity decreases. However, when the charge exceeds a certain displacement, this behavior is reversed. With further charge displacement, the conductivity decreases sharply and the viscosity increases rapidly. This behavior reversal can be traced to the formation of directional ion pairs that are present in sufficient numbers, and have lifetimes sufficiently long to strongly influence the liquid properties. We suggest that the influence of directional ion pairing can explain what appear to be anomalously low conductivities and high viscosities observed for some RTILs. The rotational and reorientational motions of the cations are examined, and shown to be strongly influenced by ion-pair formation when the charge is far off center. The temperature dependence of the transport properties is considered for selected systems, and deviations from Arrhenius behavior are found to be most important for the conductivity. Based on our results, this possibly indicates that directional ion pairs create an additional "barrier" to charge transport in some ionic liquids.

The influence of ion size disparity on structural and dynamical properties of ionic liquids is systematically investigated employing molecular dynamics simulations. Ion size ratios are varied over a realistic range (from 1:1 to 5:1) while holding other important molecular and system parameters fixed. In this way we isolate and identify effects that stem from size disparity alone. In strongly size disparate systems the larger species (cations in our model) tend to dominate the structure; the anion-anion distribution is largely determined by anion-cation correlations. The diffusion coefficients of both species increase, and the shear viscosity decreases with increasing size disparity. The influence of size disparity is strongest up to a size ratio of 3:1, then decreases, and by 5:1 both the diffusion coefficients and viscosity appear to be approaching limiting values. The conventional Stokes-Einstein expression for diffusion coefficients holds reasonably well for the cations but fails for the smaller anions as size disparity increases likely due to the neglect of strong anion-cation correlations. The electrical conductivity is not a simple monotonic function of size disparity; it first increases up to size ratios of 2:1, remains nearly constant until 3:1, then decreases such that the conductivities of the 1:1 and 5:1 systems are similar. This behavior is traced to the competing influences of ion diffusion (enhancing) and ion densities (reducing) on conductivities at constant packing fraction. The temperature dependence of the transport properties is examined for the 1:1 and 3:1 systems. In accord with experiment, the temperature dependence of all transport properties is well represented by the Vogel-Fulcher-Tammann equation. The dependence of the diffusion coefficients on the temperature/viscosity ratio is well described by the fractional Stokes-Einstein relation D proportional to (T/eta)(beta) with beta approximately = 0.8, consistent with the exponent observed for many molten inorganic salts.

Molecular dynamics simulations are used to investigate the influence of water on model ionic liquids. Several models, where the ions vary in size, and in the location of the charge with respect to the center of mass, are considered. Particular attention is focused on the variation in transport properties (diffusion coefficients, shear viscosity, and electrical conductivity) with water concentration. An effort is made to identify the underlying physical reasons for water's influence. The results for our model ionic liquids fall loosely into two categories, depending on the molecular characteristics of the constituent ions. If the ion size disparity is not too large (cation:anion diameter ratio < or approximately 2:1), and if the ion charge location is such that directional ion pair bonds are relatively weak, then we find that the ionic diffusion coefficients and the electrical conductivity increase, and the viscosity decreases with increasing water concentration. This agrees with what is commonly observed experimentally for room temperature ionic liquids (RTILs). For these systems, we do not find changes in the equilibrium structure that can account for the strong influence of water on the transport properties. Rather, by varying the molecular mass of water in our simulations, we demonstrate that the dominant effect of water can be dynamical in origin. In RTIL-water mixtures, the molecular mass of water is generally much less than that of the ions it replaces. These lighter water molecules tend to displace much heavier counterions from the ion coordination shells. This reduces caging and increases the diffusivity, which leads to higher conductivities and lower viscosities. For models with a larger ion size disparity (3:1), or in charge-off-center systems, where strong directional ion pairs are important in the pure ionic liquid, the behavior can be quite different. In these systems, the diffusion coefficients and electrical conductivity can still display conventional behavior and increase when water is added even though the reasons for this can be more complex than in the simpler cases noted above. However, in these systems the viscosity can increase, sometimes quite steeply, with increasing water concentration. We trace this unusual behavior to the formation of associated structures, extended anion-water chains that can weave among the cations in the size disparate case, and strongly bound cation-water-anion clusters in the charge-off-center systems.

Room temperature ionic liquids differ from molten salts in many ways, our work concentrates on two distinguishing features. These are large cation-anion size disparities and at least one ionic species where the center of mass and the center of charge do not coincide. In earlier work, we examined the influences of these features in isolation on simple spherical models. This paper extends this work to ionic liquid models where both features are present, and where the characteristic distance sigma(+-) (') determining the strength of the Coulombic attractions is unconstrained. We consider the interplay among these molecular features and elucidate their relative importance to the behavior of ionic liquids. Particular attention is focused on the transport properties. We find that size disparity, charge location, and sigma(+-) (') can all have large (often competing) effects. In our models, size disparity and small charge displacements lead to weakly bound, directional ion pairs, and the resulting asymmetric ion-counterion distribution gives rise to increased diffusion coefficients, consequently lower viscosity, and increased conductivity. These observations are analogous to effects reported in the literature, and we see similarities between the directional ion pairs in our models and directional cation-anion pairing through weak hydrogen bonding in room temperature ionic liquids. In our models, large charge displacements lead to strongly bound, long-lived, directional ion pairs, and in this regime the trends noted above are reversed, increased viscosities, and decreased conductivities are observed. Recently, creating more strongly hydrogen bonded, directional ion pairs has been put forward as possible means of achieving larger viscosity reductions. The trend reversal that we observe suggests that this might not work in practice.

A simple functional form for a general equation of state based on an effective near-neighbor pair interaction of an extended Lennard-Jones (12,6,3) type is given and tested against experimental data for a wide variety of fluids and solids. Computer simulation results for ionic liquids are used for further evaluation. For fluids, there appears to be no upper density limitation on the equation of state. The lower density limit for isotherms near the critical temperature is the critical density. The equation of state gives a good description of all types of fluids, nonpolar (including long-chain hydrocarbons), polar, hydrogen-bonded, and metallic, at temperatures ranging from the triple point to the highest temperature for which there is experimental data. For solids, the equation of state is very accurate for all types considered, including covalent, molecular, metallic, and ionic systems. The experimental pvT data available for solids does not reveal any pressure or temperature limitations. An analysis of the importance and possible underlying physical significance of the terms in the equation of state is given.