Relations between the Fractional Stokes−Einstein and Nernst−Einstein Equations and Velocity Correlation Coefficients in Ionic Liquids and Molten Salts

Development of Machine Learning Models for Ion-Selective Electrode Cation Sensor Design

Polyvinyl chloride (PVC) membrane-based ion-selective electrode (ISE) sensors are common tools for water assessments, but their development relies on time-consuming and costly experimental investigations. To address this challenge, this study combines machine learning (ML), Morgan fingerprint, and Bayesian optimization technologies with experimental results to develop high-performance PVC-based ISE cation sensors. By using 1745 data sets collected from 20 years of literature, appropriate ML models are trained to enable accurate prediction and a deep understanding of the relationship between ISE components and sensor performance (R 2 = 0.75). Rapid ionophore screening is achieved using the Morgan fingerprint based on atomic groups derived from ML model interpretation. Bayesian optimization is then applied to identify optimal combinations of ISE materials with the potential to deliver desirable ISE sensor performance. Na+, Mg2+, and Al3+ sensors fabricated from Bayesian optimization results exhibit excellent Nernst slopes with less than 8.2% deviation from the ideal value and superb detection limits at 10-7 M level based on experimental validation results. This approach can potentially transform sensor development into a more time-efficient, cost-effective, and rational design process, guided by ML-based techniques.

Effect of K+ Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion-ion and ion-solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K+ ion. Diffusion-based Nernst-Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

Ion Transport in Polymer Electrolytes: Building New Bridges between Experiment and Molecular Simulation

ConspectusPolymer electrolytes constitute a promising type of material for solid-state batteries. However, one of the bottlenecks for their practical implementation lies in the transport properties, often including restricted Li+ self-diffusion and conductivity and low cationic transference numbers. This calls for a molecular understanding of ion transport in polymer electrolytes in which molecular dynamics (MD) simulation can provide both new physical insights and quantitative predictions. Although efforts have been made in this area and qualitative pictures have emerged, direct and quantitative comparisons between experiment and simulation remain challenging because of the lack of a unified theoretical framework to connect them.In our work, we show that by computing the glass transition temperature (Tg) of the model system and using the normalized inverse temperature 1000/(T - Tg + 50), the Li+ self-diffusion coefficient can be compared quantitatively between MD simulations and experiments. This allows us to disentangle the effects of Tg and the polymer dielectric environment on ion conduction in polymer electrolytes, giving rise to the identification of an optimal solvating environment for fast ion conduction.Unlike Li+ self-diffusion coefficients and ionic conductivity, the transference number, which describes the fraction of current carried by Li+ ions, depends on the boundary conditions or the reference frame (RF). This creates a non-negligible gap when comparing experiment and simulation because the fluxes in the experimental measurements and in the linear response theory used in MD simulation are defined in different RFs. We show that by employing the Onsager theory of ion transport and applying a proper RF transformation, a much better agreement between experiment and simulation can be achieved for the PEO-LiTFSI system. This further allows us to derive the theoretical expression for the Bruce-Vincent transference number in terms of the Onsager coefficients and make a direct comparison to experiments. Since the Bruce-Vincent method is widely used to extract transference numbers from experimental data, this opens the door to calibrating MD simulations via reproducing the Bruce-Vincent transference number and using MD simulations to predict the true transference number.In addition, we also address several open questions here such as the time-scale effects on the ion-pairing phenomenon, the consistency check between different types of experiments, the need for more accurate force fields used in MD simulations, and the extension to multicomponent systems. Overall, this Account focuses on building new bridges between experiment and simulation for quantitative comparison, warnings of pitfalls when comparing apples and oranges, and clarifying misconceptions. From a physical chemistry point of view, it connects to concentrated solution theory and provides a unified theoretical framework that can maximize the power of MD simulations. Therefore, this Account will be useful for the electrochemical energy storage community at large and set examples of how to approach experiments from theory and simulation (and vice versa).

Ionic Conductivity of a Lithium-Doped Deep Eutectic Solvent: Glass Formation and Rotation-Translation Coupling

Deep eutectic solvents with admixed lithium salts are considered as electrolytes in electrochemical devices, such as batteries or supercapacitors. Compared to eutectic mixtures of hydrogen-bond donors and lithium salts, their raw-material costs are significantly lower. Not much is known about glassy freezing and rotational-translation coupling of such systems. Here, we investigate these phenomena by applying dielectric spectroscopy to the widely studied deep eutectic solvent glyceline, to which 1 and 5 mol % LiCl were added. Our study covers a wide temperature range, including a deeply supercooled state. The temperature dependences of the detected dipolar reorientation dynamics and ionic direct current (dc) conductivity reveal the signatures of glassy freezing. In comparison to pure glyceline, the lithium admixture leads to a reduction of ionic conductivity, which is accompanied by a reduction of the rotational dipolar mobility. However, this reduction is much smaller than that for deep eutectic solvents (DESs), where one main component is lithium salt, which we trace back to the lower glass-transition temperatures of lithium-doped DESs. In contrast to pure glyceline, the ionic and dipolar dynamics become increasingly decoupled at low temperatures and obey a fractional Debye-Stokes-Einstein relation, as previously found in other glass-forming liquids. The obtained results demonstrate the relevance of decoupling effects and glass transition to the enhancement of the technically relevant ionic conductivity in such lithium-doped DESs.

Effect of Ion Mass on Dynamic Correlations in Ionic Liquids

Ionic liquids (ILs) are a class of liquid salts with distinct properties such as high ionic conductivity, low volatility, and a broad electrochemical window, making them appealing for use in energy storage applications. The ion-ion correlations are some of the key factors that play a critical role in the ionic conductivity of ILs. In this work, we present the study of the impact of ion mass on ion-ion correlations in ILs, applying a combination of broadband dielectric spectroscopy measurements and molecular dynamics simulations. We examined three ILs with the same cation but different anions to consider three different cases of cation-anion masses: M+ > M-, M+ ≈ M-, and M+ < M-. We applied the momentum conservation approach to estimate the contribution of distinct ion-ion correlations from experimental data and obtained good agreement with direct calculations of distinct ion-ion correlations from molecular dynamics simulations. Our findings reveal that relative ion mass has a strong effect on the distinct ion-ion correlations, leading to swapping of the relative amplitude of distinct cation-cation and anion-anion correlations.

Insights into the structure and Ion transport of pectin-[BMIM][PF6] electrolytes

We investigate the effect of pectin on the structure and ion transport properties of the room-temperature ionic liquid electrolyte 1-n-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]) using molecular dynamics simulations. We find that pectin induces intriguing structural changes in the electrolyte that disrupt large ionic aggregates and promote the formation of smaller ionic clusters, which is a promising finding for ionic conductivity. Due to pectin in [BMIM][PF6] electrolytes, the diffusion coefficient of cations and anions is observed to decrease by a factor of four for a loading of 25 wt. % of pectin in [BMIM][PF6] electrolyte. A strong correlation between the ionic diffusivities (D) and ion-pair relaxation timescales (τc) is observed such that D ∼ τc-0.75 for cations and D ∼ τc-0.82 for anions. The relaxation timescale exponents indicate that the ion transport mechanisms in pectin-[BMIM][PF6] electrolytes are slightly distinct from those found in neat [BMIM][PF6] electrolytes (D∼τc-1). Since pectin marginally affects ionic diffusivities at the gain of smaller ionic aggregates and viscosity, our results suggest that pectin-ionic liquid electrolytes offer improved properties for battery applications, including ionic conductivity, mechanical stability, and biodegradability.

Search for a Grotthuss mechanism through the observation of proton transfer

The transport of protons is critical in a variety of bio- and electro-chemical processes and technologies. The Grotthuss mechanism is considered to be the most efficient proton transport mechanism, generally implying a transfer of protons between 'chains' of host molecules via elementary reactions within the hydrogen bonds. Although Grotthuss proposed this concept more than 200 years ago, only indirect experimental evidence of the mechanism has been observed. Here we report the first experimental observation of proton transfer between the molecules in pure and 85% aqueous phosphoric acid. Employing dielectric spectroscopy, quasielastic neutron, and light scattering, and ab initio molecular dynamic simulations we determined that protons move by surprisingly short jumps of only ~0.5-0.7 Å, much smaller than the typical ion jump length in ionic liquids. Our analysis confirms the existence of correlations in these proton jumps. However, these correlations actually reduce the conductivity, in contrast to a desirable enhancement, as is usually assumed by a Grotthuss mechanism. Furthermore, our analysis suggests that the expected Grotthuss-like enhancement of conductivity cannot be realized in bulk liquids where ionic correlations always decrease conductivity.

Anion and ether group influence in protic guanidinium ionic liquids

Ionic liquids are attractive liquid materials for many advanced applications. For targeted design, in-depth knowledge about their structure-property-relations is urgently needed. We prepared a set of novel protic ionic liquids (PILs) with a guanidinium cation with either an ether or alkyl side chain and different anions. While being a promising cation class, the available data is insufficient to guide design. We measured thermal and transport properties, nuclear magnetic resonance (NMR) spectra as well as liquid and crystalline structures supported by ab initio computations and were able to obtain a detailed insight into the influence of the anion and the ether substitution on the physical and spectroscopic properties. For the PILs, hydrogen bonding is the main interaction between cation and anion and the H-bond strength is inversely related to the proton affinity of the constituting acid and correlated to the increase of ¹H and ¹⁵N chemical shifts. Using anions from acids with lower proton affinity leads to proton localization on the cation as evident from NMR spectra and self-diffusion coefficients. In contrast, proton exchange was evident in ionic liquids with triflate and trifluoroacetate anions. Using imide-type anions and ether side groups decreases glass transitions as well as fragility, and accelerated dynamics significantly. In case of the ether guanidinium ionic liquids, the conformation of the side chain adopts a curled structure as the result of dispersion interactions, while the alkyl chains prefer a linear arrangement.

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.

Conduction Mechanism in Graphene Oxide Membranes with Varied Water Content: From Proton Hopping Dominant to Ion Diffusion Dominant

Proton conductors, particularly hydrated solid membranes, have various applications in sensors, fuel cells, and cellular biological systems. Unraveling the intrinsic proton transfer mechanism is critical for establishing the foundation of proton conduction. Two scenarios on electrical conduction, the Grotthuss and the vehicle mechanisms, have been reported by experiments and simulations. But separating and quantifying the contributions of these two components from experiments is difficult. Here, we present the conductive behavior of a two-dimensional layered proton conductor, graphene oxide membrane (GOM), and find that proton hopping is dominant at low water content, while ion diffusion prevails with increasing water content. This change in the conduction mechanism is attributable to the layers of water molecules in GOM nanosheets. The overall conductivity is greatly improved by forming one layer of water molecules. It reaches the maximum with two layers of water molecules, resulting from creating a complete hydrogen-bond network within GOM. When more than two layers of water molecules enter the GOM nanosheets, inducing the breakage of the ordered lamellar structure, protons spread in both in-plane and out-of-plane directions inside the GOM. Our results validate the existence of two conduction mechanisms and show their distinct contributions to the overall conductivity. Furthermore, these findings provide an optimization strategy for the design of realizing the fast proton transfer in materials with water participation.

Implications of Anion Structure on Physicochemical Properties of DBU-Based Protic Ionic Liquids

Protic ionic liquids (PILs) are potential candidates as electrolyte components in energy storage devices. When replacing flammable and volatile organic solvents, PILs are expected to improve the safety and performance of electrochemical devices. Considering their technical application, a challenging task is the understanding of the key factors governing their intermolecular interactions and physicochemical properties. The present work intends to investigate the effects of the structural features on the properties of a promising PIL based on the 1,8-diazabicyclo[5.4.0]undec-7-ene (DBUH⁺) cation and the (trifluoromethanesulfonyl)(nonafluorobutanesulfonyl)imide (IM14^–) anion, the latter being a remarkably large anion with an uneven distribution of the C–F pool between the two sides of the sulfonylimide moieties. For comparison purposes, the experimental investigations were extended to PILs composed of the same DBU-based cation and the trifluoromethanesulfonate (TFO^–) or bis(trifluoromethanesulfonyl)imide (TFSI^–) anion. The combined use of multiple NMR methods, thermal analyses, density, viscosity, and conductivity measurements provides a deep characterization of the PILs, unveiling peculiar behaviors in DBUH-IM14, which cannot be predicted solely on the basis of differences between aqueous pK_a values of the protonated base and the acid (ΔpK_a). Interestingly, the thermal and electrochemical properties of DBUH-IM14 turn out to be markedly governed by the size and asymmetric nature of the anion. This observation highlights that the structural features of the precursors are an important tool to tailor the PIL’s properties according to the specific application.

Bridging Database and Experimental Analysis to Reveal Super-hydrodynamic Conductivity Scaling Regimes in Ionic Liquids

Ion transport through electrolytes critically impacts the performance of batteries and other devices. Many frameworks used to model ion transport assume hydrodynamic mechanisms and focus on maximizing conductivity by minimizing viscosity. However, solid-state electrolytes illustrate that non-hydrodynamic ion transport can define device performance. Increasingly, selective transport mechanisms, such as hopping, are proposed for concentrated electrolytes. However, viscosity-conductivity scaling relationships in ionic liquids are often analyzed with hydrodynamic models. We report data-centric analyses of hydrodynamic transport models of viscosity-conductivity scaling in ionic liquids by merging three databases to bridge physical properties and computational descriptors. With this expansive database, we constrained scaling analyses using ion sizes defined from simulated volumes, as opposed to estimating sizes from activity coefficients. Remarkably, we find that many ionic liquids exhibit positive deviations from the Nernst-Einstein model, implying ions move faster than hydrodynamics should allow. We verify these findings using microrheology and conductivity experiments. We further show that machine learning tools can improve predictions of conductivity from molecular properties, including predictions from solely computational features. Our findings reveal that many ionic liquids exhibit super-hydrodynamic viscosity-conductivity scaling, suggesting mechanisms of correlated ion motion, which could be harnessed to enhance electrochemical device performance.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Controlling the Ion Transport Number in Solvent-in-Salt Solutions

Solvent-in-salt (SIS) systems present promising materials for the next generation of energy storage applications. The ion dynamics is significantly different in these systems from that of ionic liquids and diluted salt solutions. In this study, we analyze the ion dynamics of two salts, Li-TFSI and Li-FSI, in highly concentrated aqueous and acetonitrile solutions. We performed high-frequency dielectric measurements covering the range of up to 50 GHz and molecular dynamics simulations. The analysis of the conductivity spectra provides the characteristic crossover time between individual charge rearrangements and the normal charge diffusion regime resulting in DC conductivity. Analysis revealed that the onset of normal charge diffusion occurs at the scale of ∼1.5-3.5 Å, comparable to the average distance between the ions. Based on the idea of momentum conservation, distinct ion correlations were estimated experimentally and computationally. The analysis revealed that cation-anion correlations can be suppressed by changing the solvent concentration in SIS systems, leading to an increase of the light ion (Li⁺ in our case) transport number. This discovery suggests a way for improving the light cation transport number in SIS systems by tuning the solvent concentration.

Synthesis and Characterization of Fluorinated Phosphonium Ionic Liquids to Use as New Engineering Solvents

In this work, a set of novel fluorinated ionic liquids (FILs), based on different tetra-alkyl-phosphonium cations with perfluorobutanesulfonate and perfluoropentanoate anions, were synthesized and characterized in order to check their suitability to apply as engineering solvents. Thermophysical and thermal properties were both determined between 293.15 and 353.15 K, and the molecular volumes and free volumes and the coefficients of isobaric thermal expansion were determined from experimental values of refractive index and density. Lastly, the Walden plot was used to evaluate the ionicity of the novel FILs. The cytotoxicity of these compounds was also determined using the human hepatocellular carcinoma cells (HepG2) and the human colon carcinoma cells (Caco-2). Finally, the results were all discussed with the aim of understanding the behaviour of these compounds, considering the influence of the anion and the hydrogenated alkyl chain length. In summary, the new FILs synthesized in this work present adequate properties for their application in different industrial processes. Most of these compounds are liquid at room temperature with high decomposition temperatures. All phosphonium-based FILs have lower densities than conventional ionic liquids and common organic solvents, and the viscosity depends directly on the selected anion. Furthermore, the ionicity of FILs based on the sulfonate anion is higher than those based on the carboxylate anion. Finally, the phosphonium-based FILs have no significant effect on cellular viability at lower concentrations.

Fast evaluation technique for the shear viscosity and ionic conductivity of electrolyte solutions

With the growing need to obtain ideal materials for various applications, there is an increasing interest in computational methods to rapidly and accurately search for materials. Molecular dynamics simulation is one of the successful methods used to investigate liquid electrolytes with high transport properties applied in lithium-ion batteries. However, further reduction in computational cost is required to find a novel material with the desired properties from a large number of combinations. In this study, we demonstrate an effective fast evaluation technique for shear viscosity and ionic conductivity by molecular dynamics simulation for an exhaustive search of electrolyte materials with high transport properties. The proposed model was combined with a short-time correlation function of the stress tensor and empirical relationships to address the issues of inefficient and uncertain evaluation by conventional molecular dynamics methods. Because we focus on liquid electrolytes consisting of organic solvents and lithium salts, our model requires dissociation ratio and effective diffusion size of lithium salts. Our method is applied to search for the compositional combinations of electrolytes with superior transport properties even at low temperatures. These results correlate well with experimental results.

Highly-fluorinated Triaminocyclopropenium Ionic Liquids

A series of highly‐fluorinated triaminocyclopropenium salts, with up to six fluorous groups, were prepared and their properties as ionic liquids investigated. Reaction of pentachlorocyclopropane or tetrachlorocyclopropene with bis(2,2,2‐trifluoroethyl)amine, HN(CH₂CF₃)₂, occurs in the presence of a trialkylamine, NR₃, to give cations with two fluorinated amino groups, [C₃(N(CH₂CF₃)₂)₂(NR₂)]⁺ (R=Et, Pr, Bu, Hex), with traces of [C₃(N(CH₂CF₃)₂)₃]⁺. Use of appropriate reagent ratios and reaction times and subsequent addition of a dialkylamine, HNR'R”, gives cations with one fluorinated amino group, [C₃(N(CH₂CF₃)₂)(NR₂)(NR'R”)]⁺ ((NR₂)(NR'R”)=(NBu₂)₂, (NEt₂)(NPr₂), (NBu₂)(NBuMe)). These cations were isolated as chloride salts and some of these were converted to bistriflamide, dicyanamide and triflate salts to provide ionic liquids. These salts were characterised by thermal (DSC and TGA) measurements and miscibility/solubility properties (determined in a range of solvents). Ionic liquids (ILs) were also characterised by density, viscosity and conductivity measurements where possible. X‐ray diffraction studies of chloride salts showed the formation of fluorous regions and more hydrophilic ionic regions in the solid state.

Evidence of supercoolable nanoscale water clusters in an amorphous ionic liquid matrix

Nanoscale water clusters in an ionic liquid matrix, also called "water pockets," were previously found in some mixtures of water with ionic liquids containing hydrophilic anions. However, in these systems, at least partial crystallization occurs upon supercooling. In this work, we show for mixtures of 1-butyl-3-methylimidazolium dicyanamide with water that none of the components crystallizes up to a water content of 72 mol. %. The dynamics of the ionic liquid matrix is monitored from above room temperature down to the glass transition by combining depolarized dynamic light scattering with broadband dielectric and nuclear magnetic resonance spectroscopy, revealing that the matrix behaves like a common glass former and stays amorphous in the whole temperature range. Moreover, we demonstrate by a combination of Raman spectroscopy, small angle neutron scattering, and molecular dynamics simulation that, indeed, nanoscale water clusters exist in this mixture.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

Importance of the Ion-Pair Lifetime in Polymer Electrolytes

Ion pairing is commonly considered as a culprit for the reduced ionic conductivity in polymer electrolyte systems. However, this simple thermodynamic picture should not be taken literally, as ion pairing is a dynamical phenomenon. Here we construct model poly(ethylene oxide)-bis(trifluoromethane)sulfonimide lithium salt systems with different degrees of ion pairing by tuning the solvent polarity and examine the relation between the cation-anion distinct conductivity σ+-d and the lifetime of ion pairs τ+- using molecular dynamics simulations. It is found that there exist two distinct regimes where σ+-d scales with 1/τ+- and τ+-, respectively, and the latter is a signature of longer-lived ion pairs that contribute negatively to the total ionic conductivity. This suggests that ion pairs are kinetically different depending on the solvent polarity, which renders the ion-pair lifetime highly important when discussing its effect on ion transport in polymer electrolyte systems.

Structure-Property Relation of Trimethyl Ammonium Ionic Liquids for Battery Applications

Ionic liquids are attractive and safe electrolytes for diverse electrochemical applications such as advanced rechargeable batteries with high energy densities. Their properties that are beneficial for energy storage and conversion include negligible vapor-pressure, intrinsic conductivity as well as high stability. To explore the suitability of a series of ionic liquids with small ammonium cations for potential battery applications, we investigated their thermal and transport properties. We studied the influence of the symmetrical imide-type anions bis(trifluoromethanesulfonyl)imide ([TFSI]−) and bis(fluorosulfonyl)imide ([FSI]−), side chain length and functionalization, as well as lithium salt content on the properties of the electrolytes. Many of the samples are liquid at ambient temperature, but their solidification temperatures show disparate behavior. The transport properties showed clear trends: the dynamics are accelerated for samples with the [FSI]− anion, shorter side chains, ether functionalization and lower amounts of lithium salts. Detailed insight was obtained from the diffusion coefficients of the different ions in the electrolytes, which revealed the formation of aggregates of lithium cations coordinated by anions. The ionic liquid electrolytes exhibit sufficient stability in NMC/Li half-cells at elevated temperatures with small current rates without the need of additional liquid electrolytes, although Li-plating was observed. Electrolytes containing [TFSI]− anions showed superior stability compared to those with [FSI]− anions in battery tests.

Anion effects on Li ion transference number and dynamic ion correlations in glyme-Li salt equimolar mixtures

To achieve single-ion conducting liquid electrolytes for the rapid charge and discharge of Li secondary batteries, improvement in the Li+ transference number of the electrolytes is integral. Few studies have established a feasible design for achieving Li+ transference numbers approaching unity in liquid electrolytes consisting of low-molecular-weight salts and solvents. Previously, we studied the effects of Li+-solvent interactions on the Li+ transference number in glyme- and sulfolane-based molten Li salt solvates and clarified the relationship between this transference number and correlated ion motions. In this study, to deepen our insight into the design principles of single-ion conducting liquid electrolytes, we focused on the effects of Li+-anion interactions on Li ion transport in glyme-Li salt equimolar mixtures with different counter anions. Interestingly, the equimolar triglyme (G3)-lithium trifluoroacetate (Li[TFA]) mixture ([Li(G3)][TFA]) demonstrated a high Li+ transference number, estimated via the potentiostatic polarization method (tPPLi = 0.90). Dynamic ion correlation studies suggested that the high tPPLi could be mainly ascribed to the strongly coupled Li+-anion motions in the electrolytes. Furthermore, high-energy X-ray total scattering measurements combined with all-atom molecular dynamics simulations showed that Li+ ions and [TFA] anions aggregated into ionic clusters with a relatively long-range ion-ordered structure. Therefore, the collective motions of the Li ions and anions in the form of highly aggregated ion clusters, which likely diminish rather than enhance ionic conductivity, play a significant role in achieving high tPPLi in liquid electrolytes. Based on the dynamic ion correlations, a potential design approach is discussed to accomplish single-ion conducting liquid electrolytes with high ionic conductivity.

Microscopic origins of conductivity in molten salts unraveled by computer simulations

Molten salts are crucial materials in energy applications, such as batteries, thermal energy storage systems or concentrated solar power plants. Still, the determination and interpretation of basic physico-chemical properties like ionic conductivity, mobilities and transference numbers cause debate. Here, we explore a method for determination of ionic electrical mobilities based on non-equilibrium computer simulations. Partial conductivities are then determined as a function of system composition and temperature from simulations of molten LiF_αCl_βI_γ (with α + β + γ = 1). High conductivity does not necessarily coincide with high Li⁺ mobility for molten LiF_αCl_βI_γ systems at a given temperature. In salt mixtures, the lighter anions on average drift along with Li⁺ towards the negative electrode when applying an electric field and only the heavier anions move towards the positive electrode. In conclusion, the microscopic origin of conductivity in molten salts is unraveled here based on accurate ionic electrical mobilities and an analysis of the local structure and kinetics of the materials.

Ion Dissociation in Ionic Liquids and Ionic Liquid Solutions

The extent to which cations and anions in ionic liquids (ILs) and ionic liquid solutions are dissociated is of both fundamental scientific interest and practical importance because ion dissociation has been shown to impact viscosity, density, surface tension, volatility, solubility, chemical reactivity, and many other important chemical and physical properties. When mixed with solvents, ionic liquids provide the unique opportunity to investigate ion dissociation from infinite dilution in the solvent to a completely solvent-free state, even at ambient conditions. The most common way to estimate ion dissociation in ILs and IL solutions is by comparing the molar conductivity determined from ionic conductivity measurements such as electrochemical impedance spectroscopy (EIS) (which measure the movement of only the charged, i.e., dissociated, ions) with the molar conductivity calculated from ion diffusivities measured by pulse field gradient nuclear magnetic resonance spectroscopy (PFG-NMR, which gives movement of all of the ions). Because the NMR measurements are time-consuming, the number of ILs and IL solutions investigated by this method is relatively limited. We have shown that use of the Stokes-Einstein equation with estimates of the effective ion Stokes radii allows ion dissociation to be calculated from easily measured density, viscosity, and ionic conductivity data (ρ, η, λ), which is readily available in the literature for a much larger number of pure ILs and IL solutions. Therefore, in this review, we present values of ion dissociation for ILs and IL solutions (aqueous and nonaqueous) determined by both the traditional molar conductivity/PFG-NMR method and the ρ, η, λ method. We explore the effect of cation and anion alkyl chain length, structure, and interaction motifs of the cation and anion, temperature, and the strength of the solvent in IL solutions.

Role of Viscosity in Deviations from the Nernst-Einstein Relation

Deviations from the Nernst-Einstein relation are commonly attributed to ion-ion correlation and ion pairing. Despite the fact that these deviations can be quantified by either experimental measurements or molecular dynamics simulations, there is no rule of thumb to tell the extent of deviations. Here, we show that deviations from the Nernst-Einstein relation are proportional to the inverse viscosity by exploring the finite-size effect on transport properties under periodic boundary conditions. This conclusion is in accord with the established experimental results of ionic liquids.

Ion Mobilities, Transference Numbers, and Inverse Haven Ratios of Polymeric Ionic Liquids

We probe the ion mobilities, transference numbers, and inverse Haven ratio of ionic liquids and polymerized ionic liquids as a function of their molecular weight using a combination of atomistic equilibrium and nonequilibrium molecular dynamics simulations. In contrast to expectations, we demonstrate that the inverse Haven ratio increases with increasing degree of polymerization (N) and then decreases at larger N. For a fixed center of mass reference frame, we demonstrate that such results arise as a consequence of the strong cation-cation correlated motions, which exceed (in magnitude) the self-diffusivity of cations. Together, our findings challenge the premise underlying the pursuit of pure polymeric ionic liquids as high transference number, single-ion conducting electrolytes.

Rheological and Ionic Transport Properties of Nanocomposite Electrolytes Based on Protic Ionic Liquids and Silica Nanoparticles

In this study, the effect of hydrophilic silica nanoparticle (AEROSIL 200) addition on the rheological and transport properties of several protic ionic liquids (PILs) consisting of protonated 1,8-diazabicyclo[5.4.0]undec-7-ene cation (DBU) was studied. Interactions between the surface silanol groups of the silica nanoparticles and the ions of these PILs affected the nature of particle aggregation and the hydrogen bonding environment, which was reflected in the nonlinear rheological behaviors and transport properties of their colloidal suspensions. In contrast to shear-thinning gels formed by colloidal suspensions of the silica nanoparticles in [DBU][TFSA] ([TFSA] = [N(SO₂CF₃)₂]), [DBU][TfO] ([TfO] = [CF₃SO₃]), and [DBU][TFA] ([TFA] = [CF₃CO₂]), a shear-thickening stable suspension was formed in the [DBU][MSA] ([MSA] = [CH₃SO₃]) system. A relatively strong interaction between the silanol groups and the ions of [DBU][MSA] and the ability of this PIL to form a thicker solvation layer through hydrogen bonding were assumed to be responsible for this unique behavior. Moreover, the [DBU][MSA]-silica system showed a large enhancement in the conductivity at a certain silica concentration. This enhancement was not observed in the other PIL-silica composites that exhibited shear-thinning behavior. Even though diffusion of ions was found to be restricted in the presence of silica, a preferentially stronger interaction between [MSA] anions and the silica surface resulted in an increase in the number of charge carriers.

Impact of anion shape on Li+ solvation and on transport properties for lithium–air batteries: a molecular dynamics study

Here we report the influence of the anion shape over the solvation structure and transport properties over commonly employed Li–O₂ electrolytes and discuss their implications for the device.

Self-Consistent Scheme Combining MD and Order-N DFT Methods: An Improved Set of Nonpolarizable Force Fields for Ionic Liquids

The nonpolarizable force field of ionic liquids is tuned by using the self-consistent scheme of molecular dynamics (MD) simulation and first-principles calculation based on the order-N density functional theory (DFT). The atomic charges are determined by using the whole MD cell for DFT calculation and accounts effectively for the many-body effects of charge transfer and intramolecular polarization. The charges represent effective interactions in the condensed phase within the framework of the nonpolarizable force field and can be an alternative for an explicitly many-body model incorporating, for example, polarizability. Here we demonstrate the performance of nonpolarizable force field determined with the MD-DFT self-consistent scheme in imidazolium-, pyrrolidinium-, and ammonium-based ionic liquids. The variation ranges of molecular charges are much larger with the compositions of the ionic liquid than with the thermodynamic conditions, and the charge-ordering structures become systematically weaker with the effective charges. For energetic properties, while the calculated heat of vaporization depends on the atomic and molecular charges, the corresponding heat capacity is not strongly affected by the DFT-based variation. For transport properties, the self-diffusion coefficient, electrical conductivity, and viscosity vary much more in the self-consistent scheme. The effective DFT charge is observed to enhance the fluidity of ionic liquids and improve the accuracy of electrical conductivity and viscosity. This is due to the weakened interactions among the ions, and the too slow motions observed with a full-charge model are well corrected through the iteration of MD and DFT. We therefore conclude that the set of nonpolarizable force fields obtained with the MD-DFT self-consistent scheme leads to better description of transport properties of ionic liquids.

Structural features of selected protic ionic liquids based on a super-strong base

Protic ionic liquids (PIL) were prepared from a super-strong base 1,7-diazabicyclo[5.4.0]undec-7-ene (DBU) and super-strong acids, trifluoromethane sulfonic acid (TfOH), and (trifluoromethanesulfonyl)-(nonafluorobutylsulfonyl)imide, (IM14H), ([DBUH][TfO] and [DBUH][IM14], respectively; the latter for the first time) and their chemical and physical properties and structural features have been explored using a synergy of experimental and computational tools. The short range order in neat DBU, as well as the long range structural correlations induced by charge correlation and hydrogen bonding interactions in the ionic liquids, have been explored under ambient conditions, where these compounds are proposed for a variety of applications. Similar to other [DBUH]-based PILs, the probed ones behave as good ionic liquids. Molecular dynamics-rationalised X-ray diffraction patterns show the major role played by hydrogen bonding in affecting morphology in these systems. Additionally, we find further evidence for the existence of fluorous domains in [DBUH][IM14], thus potentially extending the range of applications for these PILs.