skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on July 28, 2026

Title: Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids
We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. Our model also shows that long-lived ion pairs result in higher viscosities.  more » « less
Award ID(s):
1944211
PAR ID:
10618432
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
AIP Publishing
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
163
Issue:
4
ISSN:
0021-9606
Subject(s) / Keyword(s):
ionic liquid, dielectric constant, glass transition, Stockmater fluid, molelar dynamics simulation
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Metal-fullerene compounds are characterized by significant electron transfer to the fullerene cage, giving rise to an electric dipole moment. We use the method of electrostatic beam deflection to verify whether such reactions take place within superfluid helium nanodroplets between an embedded C 60 molecule and either alkali (heliophobic) or rare-earth (heliophilic) atoms. The two cases lead to distinctly different outcomes: C 60 Na n ( n = 1–4) display no discernable dipole moment, while C 60 Yb is strongly polar. This suggests that the fullerene and small alkali clusters fail to form a charge-transfer bond in the helium matrix despite their strong van der Waals attraction. The C 60 Yb dipole moment, on the other hand, is in agreement with the value expected for an ionic complex. 
    more » « less
  2. Ionic liquids are currently being considered as potential electrolyte candidates for next-generation batteries and energy storage devices due to their high thermal and chemical stability. However, high viscosity and low conductivity at lower temperatures have severely hampered their commercial applications. To overcome these challenges, it is necessary to develop structure–property models for ionic liquid transport properties to guide the ionic liquid design. This work expands our previous effort in developing a machine learning model on imidazolium-based ionic liquids to now include ten different cation families, representing structural and chemical diversity. The model dataset contains 2869 ionic conductivity values over a temperature range of 238–472 K collected from the NIST ILThermo database and literature values for 397 unique ionic liquids. The database covers 214 unique cations and 68 unique anions. Three machine learning models, namely multiple linear regression, random forest, and extreme gradient boosting are applied to correlate the ionic liquid conductivity data with cation and anion features. Shapely additive analysis is performed to glean insights into cation and anion features with significant impact on ionic conductivity. Finally, the extreme gradient boosting model is used to predict the ionic conductivity of ionic liquids from all the possible combinations of unique cations and anions to identify ionic liquids crossing the ionic conductivity threshold of 2.0 S m −1 . 
    more » « less
  3. Magnetic ionic liquids (MILs) are a subclass of ionic liquids that possess a paramagnetic metal within their chemical structure, making them susceptible to external magnetic fields. A total of twenty-four (24) MILs were prepared and characterized to investigate the effect of the ligand, cation and anion on the physiochemical properties of acetylacetonate-based MILs. It was found that thermal stabilities as high as 260 °C could be achieved by incorporating aromatic moieties in the anion structure. Additionally, the magnetic moment could be modulated by simply changing the transition metal in the anion. Magnetic moment values of 2.8 μ B , 4.5 μ B and 5.6 μ B were obtained by using Ni( ii ), Co( ii ), and Mn( ii ) as the metal centers, respectively. Furthermore, the viscosity of the MILs could be tailored from a few hundred centipoise to several thousand centipoise, increasing their potential applications in numerous interdisciplinary fields. Moreover, the MILs synthesized in this study were found to be insoluble in water at a MIL-to-solvent ratio of 0.01% (w/v), making them potentially useful in targeted separations, where very hydrophobic solvents are highly desired. 
    more » « less
  4. Abstract The electron-induced secondary electron emission (SEE) yields of imidazolium-based ionic liquids are presented for primary electron beam energies between 30 and 1000 eV. These results are important for understanding plasma synthesis of nanoparticles in plasma discharges with an ionic liquid electrode. Due to their low vapor pressure and high conductivity, ionic liquids can produce metal nanoparticles in low-pressure plasmas through reduction of dissolved metal salts. In this work, the low vapor pressure of ionic liquids is exploited to directly measure SEE yields by bombarding the liquid with electrons and measuring the resulting currents. The ionic liquids studied are [BMIM][Ac], [EMIM][Ac], and [BMIM][BF4]. The SEE yields vary significantly over the energy range, with maximum yields of around 2 at 200 eV for [BMIM][Ac] and [EMIM][Ac], and 1.8 at 250 eV for [BMIM][BF4]. Molecular orbital calculations indicate that the acetate anion is the likely electron donor for [BMIM][Ac] and [EMIM][Ac], while in [BMIM][BF4], the electrons likely originate from the [BMIM]+cation. The differences in SEE yields are attributed to varying ionization potentials and molecular structures of the ionic liquids. These findings are essential for accurate modeling of plasma discharges and understanding SEE mechanisms in ionic liquids. 
    more » « less
  5. We use molecular dynamics simulations of the primitive model of electrolytes to study the ionic structure in aqueous monovalent electrolyte solutions confined by charged planar interfaces over a wide range of electrolyte concentrations, interfacial separations, surface charge densities, and ion sizes. The investigations are inspired by recent experiments that have directly measured the increase in the decay length for highly concentrated electrolytes with an increase in concentration. The behavior of ions in the nanoconfinement created by the interfaces is probed by evaluating the ionic density profiles, net charge densities, integrated charges, and decay lengths associated with the screening of the charged interface. The results show the presence of two distinct regimes of screening behavior as the concentration is changed from 0.1M to 2.5M for a wide range of electrolyte systems generated by tuning the interfacial separation, surface charge density, and ionic size. For low concentrations, the integrated charge exhibits a monotonic decay to 0 with a decay length that decreases sharply with increasing concentration. For high concentrations (≳1M), the integrated charge has a non-monotonic behavior signaling charge inversion and formation of structured layers of ions near the interfaces. The decay length under these conditions rises with increasing concentration. To complement the simulation results, a variational approach is developed that produces charge densities with characteristics consistent with those observed in simulations. The results demonstrate the relation between the rise in the strength of steric correlations and the changes in the screening behavior. 
    more » « less