Prigogine-Defay Ratio for an Ionic Glass-Former: Molecular Dynamics Simulations

The pressure dependence of the glass-transition temperature, T(g)(P), of the ionic glass-former 2Ca(NO(3))(2) center dot 3KNO(3), CKN, has been obtained by molecular dynamics (MD) simulations The liquid-glass difference of thermal expansivity, Delta alpha, heat capacity, Delta C(p), and isothermal compressibility, Delta kappa, have been calculated as a function of pressure. It has been found that the Ehrenfest relation dT(g)/dP = TV Delta alpha/Delta C(p) predicts the pressure dependence of T, but the other Ehrenfest relation, dT(g)/dP = Delta kappa/Delta alpha, does not. Consequently, the Prigogine-Defay ratio, Pi = Delta C(p)Delta kappa/TV Delta alpha(2), is Pi similar to 1.2 at low pressures, but increases 1 order of magnitude at high pressures. The pressure dependence of the Prigogine-Defay ratio is interpreted in light of recent explanations for the finding Pi > 1.

Alkoxy Chain Effect on the Viscosity of a Quaternary Ammonium Ionic Liquid: Molecular Dynamics Simulations

The viscosity of ionic liquids based on quaternary ammonium cations is reduced when one of the alkyl chains is replaced by an alkoxy chain (Zhou et al. Chem. Eur. J. 2005, 11, 752.). A microscopic picture of the role played by the ether function in decreasing the viscosity of quaternary ammonium ionic liquids is provided here by molecular dynamics (MD) simulations. A model for the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM(2)E TFSI, is compared to the tetraalky-lammonium counterpart. The alkoxy derivative has lower viscosity, higher ionic diffusion coefficients, and higher conductivity than the tetraalkyl system at the same density and temperature. A clear signature of the ether function on the liquid structure is observed in cation-cation correlations, but not in anion-anion or anion-cation correlations. In both the alkyl and the alkoxy ionic liquids, there is aggregation of long chains of neighboring cations within micelle-like structures. The MD simulations indicate that the less effective assembly between the more flexible alkoxy chains, in comparison to alkyl chains, is the structural reason for higher ionic mobility in MOENM(2)E TFSI.

Molecular dynamics simulations of highly concentrated salt solutions : structural and transport effects in polymer electrolytes

Structural, thermodynamic and transport properties have been calculated in concentrated non-aqueous NaI solutions using molecular dynamics simulations. Although the solvent has been represented by a simplistic Stockmayer fluid (spherical particles with point dipoles), the general trends observed are still a useful indication of the behavior of real non-aqueous electrolyte systems. Results indicate that in low dielectric media, significant ion pairing and clustering occurs. Contact ion pairs become more prominent at higher temperatures, independent of the dielectric strength of the solvent. Thermodynamic analysis shows that this temperature behavior is predominantly entropically driven. Calculation of ionic diffusivities and conductivities in the NaI/ether system confirms the clustered nature of the salt, with the conductivities significantly lower than those predicted from the Nernst-Einstein relation. In systems where the solvent-ion interactions increase relative to ion-ion interactions (lower charge or higher solvent dipole moment), less clustering is observed and the transport properties indicate independent motion of the ions, with higher calculated conductivities. The solvent in this system is the most mobile species, in comparison with the polymer electrolytes where the solvent is practically immobile.

Highly concentrated salt solutions : molecular dynamics simulations of structure and transport

Molecular dynamics (MD) simulations in NaI solutions, where the solvent has been represented by the Stockmayer fluid, were performed as a function of temperature, salt concentration, and solvent dipole strength. At higher temperatures contact ion pairs become more prevalent, regardless of solvent strength. An examination of the temperature dependence of the potential of mean force demonstrates the entropic nature of this effect. The transport properties calculated in the simulations are dependent on the balance between solvent dielectric constant and ion charge. In systems with a large solvent dipole moment, the ions appear to be independently mobile, and deviations from Nernst–Einstein behavior are small. In systems of smaller solvent dipole moment or greater ion charge, the ions form clusters, and large deviations from Nernst–Einstein behavior are observed.

Ion clustering in molecular dynamics simulations of sodium iodide solutions

Model systems of sodium iodide dissolved in dimethyl ether or 1,2-dimethoxyethane (glyme) were studied in order to investigate the structural and dynamic properties of ionic solutions in small and polymeric ethers. Full molecular dynamics simulations were performed at a range of different salt concentrations. An algorithm was designed which assigns ions to clusters and then calculates all the terms which contribute to ionic conductivity. In dilute solutions, free ions are the most common ionic species, followed by ion pairs. As the concentration increases, pairs become the most common species, with significant concentrations of clusters with 3 through 6 ions. Changing the solvent from dimethyl ether to glyme significantly decreases the ion clustering due to the chelate effect in which the two oxygens on a solvent stabilize an associated cation. The conductivity in stable systems is shown to be primarily the result of the movement of free ions and the relative movement of ions within neutral pairs. The Nernst-Einstein relation, commonly used in the discussion of polymer electrolytes, is shown to be inadequate to quantitatively describe conductivity in the model systems.