Temperature Dependence of Volumetric and Dynamic Properties of Imidazolium-Based Ionic Liquids

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro-Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density-temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Unraveling the Morphology of [CnC1Im]Cl Ionic Liquids Combining Cluster and Aggregation Analyses

A characteristic feature of ionic liquids is their nanosegregation, resulting in the formation of polar and nonpolar domains. The influence of increasing the alkyl side chain on the morphology of ionic liquids has been the subject of many studies. Typically, the polar network (charged part of the cation and anion) constitutes a continuous subphase that partially breaks to allow the formation of a nonpolar domain with the increase of the alkyl chain. As the nonpolar network expands, the number of tails per aggregate increases until the ionic liquid percolates. In this work, we demonstrate how the complementary software packages TRAVIS and AGGREGATES can be employed in conjunction to gain insights into the size and morphology of the [CnC1Im]Cl family, with n ∈ {2, 4, 6, 8, 10, 12}. The combination of the two approaches rounds off the picture of the intricate arrangement and structural features of the alkyl chains.

Molecular Dynamics Simulation of the Influence of External Electric Fields on the Glass Transition Temperature of the Ionic Liquid 1-Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide

We present the results of molecular dynamics simulations of the ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [C₂C₁im][NTf₂] in the presence of external electric fields (EEFs) of varying strengths to understand the effects of EEFs on the glass transition temperature T_g. We compute T_g with an automated and objective method and observe a depression in T_g when cooling the IL within an EEF above a critical strength. The effect is reversible, and glasses prepared with EEFs recover their original zero-field T_g when heated. By examining the dynamics and structure of the liquid phase, we find that the EEF lowers the activation energy for diffusion, reducing the energetic barrier for movement and consequently T_g. We show that the effect can be leveraged to drive an electrified nonvapor compression refrigeration cycle.

Molecular Modelling of Ionic Liquids: Situations When Charge Scaling Seems Insufficient

Charge scaling as an effective solution to the experiment-computation disagreement in molecular modelling of ionic liquids (ILs) could bring the computational results close to the experimental reference for various thermodynamic properties. According to the large-scale benchmark calculations of mass density, solvation, and water-ILs transfer-free energies in our series of papers, the charge-scaling factor of 0.8 serves as a near-optimal option generally applicable to most ILs, although a system-dependent parameter adjustment could be attempted for further improved performance. However, there are situations in which such a charge-scaling treatment would fail. Namely, charge scaling cannot really affect the simulation outcome, or minimally perturbs the results that are still far from the experimental value. In such situations, the vdW radius as an additional adjustable parameter is commonly tuned to minimize the experiment-calculation deviation. In the current work, considering two ILs from the quinuclidinium family, we investigate the impacts of this vdW-scaling treatment on the mass density and the solvation/partition thermodynamics in a fashion similar to our previous charge-scaling works, i.e., scanning the vdW-scaling factor and computing physical properties under these parameter sets. It is observed that the mass density exhibits a linear response to the vdW-scaling factor with slopes close to -1.8 g/mL. By further investigating a set of physiochemically relevant temperatures between 288 K and 348 K, we confirm the robustness of the vdW-scaling treatment in the estimation of bulk properties. The best vdW-scaling parameter for mass density would worsen the computation of solvation/partition thermodynamics, and a marginal decrease in the vdW-scaling factor is considered as an intermediate option balancing the reproductions of bulk properties and solvation thermodynamics. These observations could be understood in a way similar to the charge-scaling situation. i.e., overfitting some properties (e.g., mass density) would degrade the accuracy of the other properties (e.g., solvation free energies). Following this principle, the general guideline for applying this vdW-tuning protocol is by using values between the density-derived choice and the solvation/partition-derived solution. The charge and current vdW scaling treatments cover commonly encountered ILs, completing the protocol for accurate modelling of ILs with fixed-charge force fields.

Tracing the Effect of Replacing [Gly]- with [Ala]- and Hydroxylation of [emim]+ on the Fine-Tuning of the Transport Properties of the Corresponding Amino Acid-Based Ionic Liquids Using MD Simulation

Natural amino acid-based ionic liquids (AAILs) composed of deprotonated amino acids, [AA]^-, as anions and hydroxylated imidazolium cations provide an eco-friendly nontoxic IL family with the growing number of chemical and biochemical revolutionary applications. In this paper, the transport properties of four AAILs composed of 1-(2-hydroxyethyl)-3-methylimidazolium ([HOemim]⁺) and 1-ethyl-3-methylimidazolium ([emim]⁺) cations with alaninate and glycinate anions were studied by molecular dynamics (MD) simulations. A nonpolarizable all-atom force field with the scaled charge (±0.8e) on each of the ions was applied and compared with the unit charge model in some cases. The tunable effects of the presence of the hydroxyl group in the side chain of the imidazolium cation, the type of amino acid anion, and the varied temperature on the dynamical behavior of AAILs were investigated in detail. The experimentally compatible trends of the simulated ionic self-diffusion coefficients, ionic conductivity, and ionicity were found to be inverse to the viscosity and ionic association of these ILs as [emim][Gly] > [emim][Ala] > [HOemim][Gly] > [HOemim][Ala]. The main reason behind these trends is the higher ability of the hydroxylated cation for the hydrogen-bond formation with [AA]^-. The mean square displacement (MSD), self-diffusion, and transference number of imidazolium cations are larger than those of [AA]^- anions, except in the case of [HOemim][Gly]. It was found that the activation energy for diffusion of [AA]^- is lower than that of [HOemim]⁺ but higher than that of [emim]⁺ in [HOemim][AA] and [emim][AA] ILs, respectively. The computed velocity autocorrelation function (VACF) showed that [Gly]^-, as the lightest ion, has the shortest mean collision time and velocity randomization time among the ions, especially in the [HOemim][Gly] IL. Replacing [emim]⁺ with [HOemim]⁺, similar to the effect of decreasing temperature, causes significant decreasing of the ionic self-diffusion and increasing of the well depth of the first minimum of the ionic VACFs. Current findings show that introducing suitable functional groups in the side chain of imidazolium cations can be a viable approach for efficient engineering design and fine-tuning of the transport properties of these AAILs.

How Temperature, Pressure, and Salt Concentration Affect Correlations in LiTFSI/EMIM-TFSI Electrolytes: A Molecular Dynamics Study

Classical polarizable molecular dynamics simulations have been performed for LiTFSI solutions in the EMIM-TFSI ionic liquid. Different temperature or pressure values and salt concentrations have been examined. The structure and dynamics of the solvation shell of Li+ cations, diffusion coefficients of ions, conductivities of the electrolytes, and correlations between motions of ions have been analyzed. The results indicated that regardless of the conditions, significant correlations are present in all systems. The degree of correlations depends mainly on the salt fraction in the electrolyte and is much less affected by temperature and pressure changes. A positive correlation between motions of Li+ cations and TFSI anions, leading to the occurrence of negative Li+ transference numbers, exists for all conditions, although temperature and pressure changes affect the speed of anion exchange in Li+ solvation shells.

Extending the timescale of molecular simulations by using time-temperature superposition: rheology of ionic liquids

Molecular dynamics simulations are used to determine the temperature dependence of the dynamic and rheological properties of a model imidazolium-based ionic liquid (IL). The simulation results for the volumetric properties of the IL are in good agreement with the experimental results. The temperature dependence of the diffusion coefficient of anions and cations follows the Vogel-Fulcher-Tammann equation over the range of the temperatures studied. The shear viscosity of the IL shows a Newtonian plateau at low shear rates and shear-thinning behavior at high shear rates. The dynamic modulus values indicate that the IL behaves like a viscous liquid at high temperatures and low frequencies, while its viscoelastic response becomes similar to that of an elastic solid at low temperatures and high frequencies. Using the time-temperature superposition (TTS) principle, the dynamic moduli, shear viscosity, and mean squared displacement of cations and anions in the diffusive regime can be collapsed onto master curves by applying a single set of shift factors. Due to the large mismatch in the timescale investigated by the atomistically detailed simulations and experiments, the glass transition temperature predicted in simulations shifts to higher values. When this timescale mismatch is accounted for by using appropriate shift factors, the master curves of the dynamic moduli obtained in simulations closely match those obtained in experiments. This result demonstrates the exciting ability of TTS to overcome the large timescale disparity between simulations and experiments which will enable the use of molecular simulations for quantitatively predicting the rheological property values at frequencies of practical interest.

The Dependence of Ionic Liquid Solvent Effects on the Nucleophilic Heteroatom in SN Ar Reactions. Highlighting the Potential for Control of Selectivity

Nucleophilic aromatic substitution (S_N Ar) reactions of 1-fluoro-4-nitrobenzene using similar nitrogen and sulfur nucleophiles were studied through extensive kinetic analysis in mixtures containing ionic liquids. The interactions of the ionic liquid components with the starting materials and transition state for each process were investigated in an attempt to construct a broad predictive framework for how ionic liquids affect reaction outcome. It was found that, based on the activation parameters, the microscopic interactions and thus the ionic liquid solvent effect were different for each of the nucleophiles considered. The results from this study suggest that it may be possible to rationally select a given ionic liquid mixture to selectively control reaction outcome of an S_N Ar reaction where multiple nucleophiles are present.

Protic Ionic Liquids Based on the Alkyl-Imidazolium Cation: Effect of the Alkyl Chain Length on Structure and Dynamics

Protic ionic liquids are known to form extended hydrogen-bonded networks that can lead to properties different from those encountered in the aprotic analogous liquids, in particular with respect to the structure and transport behavior. In this context, the present paper focuses on a wide series of 1-alkyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquids, [HC _nIm][TFSI], with the alkyl chain length ( n) on the imidazolium cation varying from ethyl ( n = 2) to dodecyl ( n = 12). A combination of several methods, such as vibrational spectroscopy, wide-angle X-ray scattering (WAXS), broadband dielectric spectroscopy, and ¹H NMR spectroscopy, is used to understand the correlation between local cation-anion coordination, nature of nanosegregation, and transport properties. The results indicate the propensity of the -NH site on the cation to form stronger H-bonds with the anion as the alkyl chain length increases. In addition, the position and width of the scattering peak q₁ (or the pre-peak), resolved by WAXS and due to the nanosegregation of the polar from the nonpolar domains, are clearly dependent on the alkyl chain length. However, we find no evidence from pulsed-field gradient NMR of a proton motion decoupled from molecular diffusion, hypothesized to be facilitated by the longer N-H bonds localized in the segregated ionic domains. Finally, for all protic ionic liquids investigated, the ionic conductivity displays a Vogel-Fulcher-Tammann dependence on inverse temperature, with an activation energy E_a that also depends on the alkyl chain length, although not strictly linearly.

Probing translational and rotational dynamics in hydrophilic/hydrophobic anion based imidazolium ionic liquid-water mixtures

In this investigation, we examine the effect of water concentration and temperature on the dynamical properties of [Hmim][Cl] and [Hmim][NTf2] ionic liquids (ILs). The dynamical properties such as translational diffusion coefficients, ion-pair lifetimes, and rotational correlation times are calculated using molecular dynamics simulations. The simulations predict that water concentration also significantly impacts the magnitude of dynamical properties. At low, intermediate and high water concentrations, the following trend in diffusion coefficients is seen: Cl- > Hmim+; Cl- > NTf2-; Hmim+ ([Hmim][Cl]) > Hmim+ ([Hmim] [NTf2]). At ultra-low water concentrations of [Hmim][Cl] IL, several bridge like configurations form between water molecules and Cl- anions, which are supported by a complex distribution of water clusters. The effect of an increase in the water concentration leads to a decrease in ion-pair lifetimes between the Hmim+ cations and Cl-/NTf2- anions, which strongly correlates with the trends observed from the diffusion coefficients. A biexponential function was found to be the best fit for the RACF at neat/ultra-low water concentrations of [Hmim][Cl] and [Hmim][NTf2] ILs, whereas a single exponential function was sufficient to fit the RACF at low, intermediate and high water concentrations. The rotational relaxation time of the Hmim+ cations is larger in neat [Hmim][Cl] compared to that in neat [Hmim][NTf2] with an opposite trend seen with hydration. The rotational correlation time of water molecules is larger in [Hmim][Cl] compared to that in [Hmim][NTf2] at low and intermediate water concentrations, with similar correlation times observed at high water concentrations.