Li+ Solvation and Transport Properties in Ionic Liquid/Lithium Salt Mixtures: A Molecular Dynamics Simulation Study

Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a protocol to enhance the quality of group 1 monatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting preexisting individual parameters for ions, model compounds, and water, we employ atom-pair specific Lennard-Jones (LJ) (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methylacetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.

Anomalous diffusion of lithium-anion clusters in ionic liquids

Lithium-ion transport is significantly retarded in ionic liquids (ILs). In this work, we performed extensive molecular dynamics simulations to mimic the kinetics of lithium ions in ILs using [N-methyl-N-propylpyrrolidium (pyr[Formula: see text])][bis(trifluoromethanesulfonyl)imide (Ntf[Formula: see text])] with added LiNtf[Formula: see text] salt. And we analyzed their transport, developing a two-state model and comparing it to the machine learning-identified states. The transport of lithium ions involves local shell exchanges of the Ntf[Formula: see text] in the medium. We calculated train size distributions over various time scales. The train size distribution decays as a power law, representing non-Poissonian bursty shell exchanges. We analyzed the non-Poissonian processes of lithium ions transport as a two-state (soft and hard) model. We analytically calculated the transition probability of the two-state model, which fits well to the lifetime autocorrelation functions of LiNtf[Formula: see text] shells. To identify two states, we introduced the graph neutral network incorporating local molecular structure. The results reveal that the shell-soft state mainly contributes to the transport of the lithium ions, and their contribution is more important in low temperatures. Hence, it is the key for enhanced lithium ion transport to increase the fraction of the shell-soft state.

Polymer Architecture-Induced Trade-off between Conductivities and Transference Numbers in Salt-Doped Polymeric Ionic Liquids

Recent experiments have demonstrated that polymeric ionic liquids that share the same cation and anion but possess different architectures can exhibit markedly different conductivity and transference number characteristics when doped with lithium salt. In this study, we used atomistic molecular simulations on polymer chemistries inspired by the experiments to probe the mechanistic origins underlying the competition between conductivity and transference numbers. Our results indicate that the architecture of the polycationic ionic liquid plays a subtle but crucial role in modulating the anion-cation interactions, especially their dynamical coordination characteristics. Chemistries leading to longer-lived anion-cation coordinations relative to lithium-anion coordinations lead to lower conductivities and higher transference numbers. Our results suggest that higher conductivities are accompanied by lower transference numbers and vice versa, revealing that alternative approaches may need to be considered to break this trade-off in salt-doped polyILs.

Regulating the Solvation Shell Structure of Lithium Ions for Smooth Li Metal Deposition in Quasi-Solid-State Batteries

Gel polymer electrolytes (GPE) are promising next-generation electrolytes for high-energy batteries, combining the multiple advantages of liquid and all-solid-state electrolytes. Herein, we a synthesized GPE using poly(ethylene glycol)acrylate (PEGDA) in order to understand how the GPE efficiently inhibits lithium dendrite formation and growth. The effects of PEGDA on the solvation shell structure of the lithium ion are investigated using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations, which are also supported by Raman spectroscopy. The GPE electrolytes with optimal PEGDA concentration exhibit high transference numbers (t Li + ${{_{{\rm Li}{^{+}}}}}$ =0.72) and ionic conductivity (σ=3.24 mS cm^-1 ). A symmetric lithium ion battery using GPE can be stably cycled for 1200 h in comparison to 320 h in a liquid electrolyte (LE), possibly owing to the high content of LiF (17.9 %) in the solid-electrolyte interphase film of the GPE cell. The observed concentration/electric field gradient observed through the finite element method also accounts for the good cycling performance. In addition, a LiCoO₂ |GPE|Li cell demonstrates excellent capacity retention of 87.09 % for 200 cycles; this approach could present promising guidelines for the design of high-energy lithium batteries.

Regulating the solvation chemistry of non-flammable high voltage electrolyte through salt-solvent ratio modulation

Boosting the energy density and safety issue of lithium-ion batteries has become ever more important to satisfy the diverse applications such as energy storage and mobile electronic devices. Herein, we present a new high voltage polyether-based electrolyte (HVPEE) by solvation structure design that can endure high-voltage operations and also possess non-flammable features. Especially, HVPEEs show better compatibility and stability with electrode than conventional electrolyte. We find that the solvent separated ion pair (SSIP) and contact ion pair (CIP) dominate the ion-solvent structure of HVPEEs, rather than the free solvent and ions. In this way, the oxidative decomposition of HVPEE on the cathode interface can be suppressed significantly due to the reduced highest occupied molecular orbital of SSIP complex structure than that of free TFSI^-. As a result, the oxidation voltage can achieve as high as 5.35 V when the ether group/Li is optimized at 10/1 in the HVPEE, enabling the LiFePO₄//Li full cells deliver a capacity of 165 mA h g^-1 with a capacity retention of 98 % after 200 cycles. Moreover, when the cut-off voltage is 4.5 V, the discharge capacity of the LiNi_0.6Mn_0.2Co_0.2O₂//Li full cell can reach 170 mA h g^-1.

How NaFTA salt affects the structural landscape and transport properties of Pyrr1,3FTA ionic liquid

Recently, it has been demonstrated that ionic liquids (ILs) with an asymmetric anion render a wider operational temperature range and can be used as a solvent in sodium ion batteries. In the present study, we examine the microscopic structure and dynamics of pure 1-methyl-1-propylpyrrolidinium fluorosulfonyl(trifluoromethylsulfonyl)amide (Pyrr_1,3FTA) IL using atomistic molecular dynamics simulations. How the addition of the sodium salt (NaFTA) having the same anion changes the structural landscape and transport properties of the pure IL has also been explored. The simulated x-ray scattering structure functions reveal that the gradual addition of NaFTA salt (up to 1.2 molal) suppresses the charge alternating feature of the pure IL because of the replacement of the Pyrr⁺ cations with the Na⁺ ions. The Na⁺ ions are majorly found near the oxygen atoms of the anions, but the probability of finding the Na⁺ ions near these atoms slightly decreases with increasing salt concentration. As expected, the Na⁺ ions stay away from the Pyrr⁺ cations. However, the probability of finding the anions around anions increases with increasing salt concentration. The simulated self-diffusion coefficients of the ions in the pure IL reveal slightly faster diffusion of the Pyrr⁺ cations as compared to the FTA^- anions. Interestingly, in the salt solution, despite having smaller size, the diffusion of the Na⁺ ions is found to be lesser than the Pyrr⁺ cations and the FTA^- anions. The analysis of the ionic conductivity and transport numbers reveals that the fractional contribution of the FTA^- anion to the overall conductivity remains nearly constant with increasing salt concentration, but the contribution of Pyrr⁺ cation decreases and Na⁺ ion increases.

Modulation of Diffusion Mechanism and Its Correlation with Complexation in Aqueous Deep Eutectic Solvents

Aqueous mixtures of deep eutectic solvents (DESs) have gained traction recently as an effective template to tailor their physicochemical properties. But detailed microscopic insights into the effects of water on the molecular relaxation phenomenon in DESs are not entirely understood. DESs are strong network-forming liquids due to the extensive hydrogen bonding and complex formation between their species, and therefore, water can behave as a controlled disruptor altering the microscopic structure and dynamics in DESs. In this study, the role of water in the diffusion mechanism of acetamide in the aqueous mixtures of DESs synthesized using acetamide and lithium perchlorate is investigated using molecular dynamics (MD) simulation and quasielastic neutron scattering (QENS). The acetamide dynamics comprises localized diffusion within transient cages and a jump diffusion process across cages. The jump diffusion process is observed to be strongly enhanced by about a factor of 10 as the water content in the system is increased. Meanwhile, the geometry of the localized dynamics is unaltered by addition of water, but the localized diffusion becomes significantly faster and more heterogeneous with increasing water concentration. The accelerating effects of water on localized diffusion are also substantiated by QENS experiments. The water concentration in the DES is observed to control the solvation structure of lithium ions, with the ions becoming significantly hydrated at 20 wt % water. The formation of interwater and water-acetamide hydrogen bonds is observed. The increase in water concentration is found to increase the number of H-bonds; however, their lifetimes are found to decrease substantially. Similarly, the lifetimes of acetamide-lithium complexes are also found to be diminished by increasing water concentration. A power-law scaling relationship between lifetimes and diffusion constants is established, elucidating the extent of coupling between diffusive processes and hydrogen bonding and microscopic complexation. This study demonstrates the ability to use water as an agent to probe the role of structural relaxation and complex lifetimes of diffusive processes at different time and length scales.

Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions

Polarizable molecular dynamics simulations are a fast progressing field in the scientific research of ionic liquids. The fundamentals of polarizable simulations, as well as their application to ionic liquids, were summarized in a review [Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, Jr., A. D.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem. Rev. 2019, 119, 7940–7995] in 2019. Since then, new methods to treat intermolecular interaction of induced dipoles in these highly charged systems were developed. This concerns the damping of these interactions and additional charge transfer as well as the prediction of ionic materials with ultrahigh refractive indices. In addition to the progress of the polarizable force fields, also thermostats and barostats for polarizable simulations evolved recently.

Ion Transport in the EMITFSI/PVDF System at Different Temperatures: A Molecular Dynamics Simulation

We used all-atom molecular dynamics simulations to study the ion transport in the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide/poly(vinylidene fluoride) (EMITFSI/PVDF) system with 40.05 wt % EMITFSI at different temperatures. The glass-transition temperature (T _g = 204 K) of this system shows a good agreement with the experimental value (200 K). With the increase of temperature, the peaks of the pair correlation function show an increasing trend. Interestingly, the coordination numbers of ion pairs and the degree of independent ion motion are mainly affected by the binding energy between ion pairs as the temperature increases. In addition, the ion transport properties with increasing temperature can be studied by the ion-pair relaxation times, ion-pair lifetimes, and diffusion coefficients. The simulation results illustrate that the ion transport is intensified. Especially, the cations can always diffuse faster than the anions. The power law shows that mobilities of anions and cations are seen to exhibit a "superionic" behavior. With the increase of temperature, transference numbers of anions decrease first and then increase and transference numbers of cations show the opposite changes; ionic conductivity increases gradually; and viscosity decreases gradually, indicating that the diffusion resistance of ions decreases. In general, after adding PVDF into the EMITFSI system, the glass-transition temperature and viscosity increase, the ionic conductivity and degree of independent ion motion decrease, and diffusion coefficients of cations decrease faster than those of the anions.

The effect of explicit counterion binding on the transference number of polyelectrolyte solutions

Polyelectrolyte solutions have been proposed as a method to improve the efficiency of lithium-ion batteries by increasing the cation transference number because the polymer self-diffusion coefficient is much lower than that of the counterion. However, this is not necessarily true for the polymer mobility. In some cases, negative transference numbers have been reported, which implies that the lithium ions are transporting to the same electrode as the anion, behavior that is often attributed to a binding of counterions to the polyion. We use a simple model where we bind some counterions to the polymer via harmonic springs to investigate this phenomenon. We find that both the number of bound counterions and the strength of their binding alter the transference number, and, in some cases, the transference number is negative. We also investigate how the transference number depends on the Manning parameter, the ratio of the Bjerrum length to charge separation along the chain. By altering the Manning parameter, the transference number can almost be doubled, which suggests that charge spacing could be a way to increase the transference number of polyelectrolyte solutions.

Controlling Li+ transport in ionic liquid electrolytes through salt content and anion asymmetry: a mechanistic understanding gained from molecular dynamics simulations

In this work, we report the results from molecular dynamics simulations of lithium salt-ionic liquid electrolytes (ILEs) based either on the symmetric bis[(trifluoromethyl)sulfonyl]imide (TFSI^-) anion or its asymmetric analogue 2,2,2-(trifluoromethyl)sulfonyl-N-cyanoamide (TFSAM^-). Relating lithium's coordination environment to anion mean residence times and diffusion constants confirms the remarkable transport behaviour of the TFSAM^--based ILEs that has been observed in recent experiments: for increased salt doping, the lithium ions must compete for the more attractive cyano over oxygen coordination and a fragmented landscape of solvation geometries emerges, in which lithium appears to be less strongly bound. We present a novel, yet statistically straightforward methodology to quantify the extent to which lithium and its solvation shell are dynamically coupled. By means of a Lithium Coupling Factor (LCF) we demonstrate that the shell anions do not constitute a stable lithium vehicle, which suggests for this electrolyte material the commonly termed "vehicular" lithium transport mechanism could be more aptly pictured as a correlated, flow-like motion of lithium and its neighbourhood. Our analysis elucidates two separate causes why lithium and shell dynamics progressively decouple with higher salt content: on the one hand, an increased sharing of anions between lithium limits the achievable LCF of individual lithium-anion pairs. On the other hand, weaker binding configurations naturally entail a lower dynamic stability of the lithium-anion complex, which is particularly relevant for the TFSAM^--containing ILEs.

Salt-in-Ionic-Liquid Electrolytes: Ion Network Formation and Negative Effective Charges of Alkali Metal Cations

Salt-in-ionic liquid electrolytes have attracted significant attention as potential electrolytes for next generation batteries largely due to their safety enhancements over typical organic electrolytes. However, recent experimental and computational studies have shown that under certain conditions alkali cations can migrate in electric fields as if they carried a net negative effective charge. In particular, alkali cations were observed to have negative transference numbers at small mole fractions of alkali-metal salt that revert to the expected net positive transference numbers at large mole fractions. Simulations have provided some insights into these observations, where the formation of asymmetric ionic clusters, as well as a percolating ion network, could largely explain the anomalous transport of alkali cations. However, a thermodynamic theory that captures such phenomena has not been developed, as ionic associations were typically treated via the formation of ion pairs. The theory presented herein, based on the classical polymer theories, describes thermoreversible associations between alkali cations and anions, where the formation of large, asymmetric ionic clusters and a percolating ionic network are a natural result of the theory. Furthermore, we present several general methods to calculate the effective charge of alkali cations in ionic liquids. We note that the negative effective charge is a robust prediction with respect to the parameters of the theory and that the formation of a percolating ionic network leads to the restoration of net positive charges of the cations at large mole fractions of alkali metal salt. Overall, we find excellent qualitative agreement between our theory and molecular simulations in terms of ionic cluster statistics and the effective charges of the alkali cations.

Computational Study of the Structure and Transport in Pyrrolidinium-Li-TFSI-Silica Ionogels

Ionogels (IGs) are a unique class of composite materials with attributes that make them promising materials for applications in electrochemical energy storage. Due to the solid porous matrix that confines the ionic liquid (IL) in the IG, they can be used as self-supporting electrolytes. Furthermore, interactions of the IL with the porous matrix can have beneficial effects on transport, such as lowering the freezing/glass transition temperature of the conducting IL. In this work, we employ molecular dynamics simulations to investigate the influence of the porous morphology and solid volume fraction on ionic conductivity and Li⁺ diffusivity using a representative 0.5 M Li-bis(trifluoromethane)sulfonimide (TFSI)-pyrrolidinium (Pyr1.3) IL confined in a nanoporous silica matrix. The effect of the morphology of the confining matrix is compared using the pure IL as a baseline. We find that the tracer and collective Li⁺ diffusion and ionic conductivity of all the model IGs have significantly lower temperature dependence than the corresponding pure IL. In general, low-silica IGs with wide pores displayed the best transport properties at high temperatures, but the trends with the morphology for the nested set of transport coefficients we examined changed as the collective behavior of the Li⁺ ions and the molecular IL components were considered. Remarkably, some of the model IGs displayed better transport properties on a volume of fluid basis at low temperatures than the constituent IL. These trends were tied to structural changes revealed by the radial distribution functions of the IL components and the silica surface, including a decreasing Li⁺ adsorption peak of the surface silica indicating a change in the relative contributions of bulk-like and surface-like transport in the confined IL.

Can the microscopic and macroscopic transport phenomena in deep eutectic solvents be reconciled?

Deep eutectic solvents (DESs) have become ubiquitous in a variety of industrial and pharmaceutical applications since their discovery. However, the fundamental understanding of their physicochemical properties and their emergence from the microscopic features is still being explored fervently. Particularly, the knowledge of transport mechanisms in DESs is essential to tune their properties, which shall aid in expanding the territory of their applications. This perspective presents the current state of understanding of the bulk/macroscopic transport properties and microscopic relaxation processes in DESs. The dependence of these properties on the components and composition of the DES is explored, highlighting the role of hydrogen bonding (H-bonding) interactions. Modulation of these interactions by water and other additives, and their subsequent effect on the transport mechanisms, is also discussed. Various models (e.g. hole theory, free volume theory, etc.) have been proposed to explain the macroscopic transport phenomena from a microscopic origin. But the formation of H-bond networks and clusters in the DES reveals the insufficiency of these models, and establishes an antecedent for dynamic heterogeneity. Even significantly above the glass transition, the microscopic relaxation processes in DESs are rife with temporal and spatial heterogeneity, which causes a substantial decoupling between the viscosity and microscopic diffusion processes. However, we propose that a thorough understanding of the structural relaxation associated to the H-bond dynamics in DESs will provide the necessary framework to interpret the emergence of bulk transport properties from their microscopic counterparts.

Water in Protic Ionic Liquid Electrolytes: From Solvent Separated Ion Pairs to Water Clusters

The large electrochemical and cycling stability of "water-in-salt" systems have rendered promising prospective electrolytes for batteries. The impact of addition of water on the properties of ionic liquids has already been addressed in several publications. In this contribution, we focus on the changes in the state of water. Therefore, we investigated the protic ionic liquid N-butyl-pyrrolidinium bis(trifluoromethanesulfonyl)imide with varying water content at different temperatures with the aid of molecular dynamics simulations. It is revealed that at very low concentrations, the water is well dispersed and best characterized as shared solvent molecules. At higher concentrations, the water forms larger aggregates and is increasingly approaching a bulk-like state. While the librational and rotational dynamics of the water molecules become faster with increasing concentration, the translational dynamics are found to become slower. Further, all dynamics are found to be faster if the temperature increases. The trends of these findings are well in line with the experimental measured conductivities.

On the Coordination Chemistry of the lanthanum(III) Nitrate Salt in EAN/MeOH Mixtures

A thorough structural characterization of the La(NO₃)₃ salt dissolved into several mixtures of ethyl ammonium nitrate (EAN) and methanol (MeOH) with EAN molar fraction χ_EAN ranging from 0 to 1 has been carried out by combining molecular dynamics (MD) and X-ray absorption spectroscopy (XAS). The XAS and MD results show that changes take place in the La³⁺ first solvation shell when moving from pure MeOH to pure EAN. With increasing the ionic liquid content of the mixture, the La³⁺ first-shell complex progressively loses MeOH molecules to accommodate more and more nitrate anions. Except in pure EAN, the La³⁺ ion is always able to coordinate both MeOH and nitrate anions, with a ratio between the two ligands that changes continuously in the entire concentration range. When moving from pure MeOH to pure EAN, the La³⁺ first solvation shell passes from a 10-fold bicapped square antiprism geometry where all the nitrate anions act only as monodentate ligands to a 12-coordinated icosahedral structure in pure EAN where the nitrate anions bind the La³⁺ cation both in mono- and bidentate modes. The La³⁺ solvation structure formed in the MeOH/EAN mixtures shows a great adaptability to changes in the composition, allowing the system to reach the ideal compromise among all of the different interactions that take place into it.

The structural properties of the La(NO₃)₃ salt dissolved into EAN/methanol mixtures were characterized by molecular dynamics and X-ray absorption spectroscopy. The La³⁺ solvation shell undergoes significant changes with increasing the ionic liquid content of the mixture, progressively losing methanol molecules to accommodate more and more nitrate anions. The La³⁺ solvation structure shows great adaptability to composition changes, passing from a 10-fold bicapped square antiprism geometry in pure methanol to a 12-coordinated icosahedral complex in EAN.

Water accelerates the hydrogen-bond dynamics and abates heterogeneity in deep eutectic solvent based on acetamide and lithium perchlorate

Deep eutectic solvents (DESs) have become a prevalent and promising medium in various industrial applications. The addition of water to DESs has attracted a lot of attention as a scheme to modulate their functionalities and improve their physicochemical properties. In this work, we study the effects of water on an acetamide based DES by probing its microscopic structure and dynamics using classical molecular dynamics simulation. It is observed that, at low water content, acetamide still remains the dominant solvate in the first solvation shell of lithium ions, however, beyond 10 wt. %, it is replaced by water. The increase in the water content in the solvent accelerates the H-bond dynamics by drastically decreasing the lifetimes of acetamide-lithium H-bond complexes. Additionally, water-lithium H-bond complexes are also found to form, with systematically longer lifetimes in comparison to acetamide-lithium complexes. Consequently, the diffusivity and ionic conductivity of all the species in the DES are found to increase substantially. Non-Gaussianity parameters for translational motions of acetamide and water in the DES show a conspicuous decrease with addition of water in the system. The signature of jump-like reorientation of acetamide is observed in the DES by quantifying the deviation from rotational Brownian motion. However, a notable decrease in the deviation is observed with an increase in the water content in the DES. This study demonstrates the intricate connection between H-bond dynamics and various microscopic dynamical parameters in the DES, by investigating the modulation of the former with addition of water.

Multicomponent Phase Separation in Ternary Mixture Ionic Liquid Electrolytes

We investigate the phase behavior of ternary mixtures of ionic liquid, organic solvent, and lithium salt by molecular dynamics simulations. We find that at room temperature, the electrolyte separates into distinct phases with specific compositions; an ion-rich domain that contains a fraction of solvent molecules and a second domain of pure solvent. The phase separation is shown to be entropy-driven and is independent of lithium salt concentration. Phase separation is only observed at microsecond time scales and greatly affects the transport properties of the electrolyte.

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.

The Hydrotropic Effect of Ionic Liquids in Water-in-Salt Electrolytes*

Water-in-salt electrolytes have successfully expanded the electrochemical stability window of aqueous electrolytes beyond 2 V. Further improvements in stability can be achieved by partially substituting water with either classical organic solvents or ionic liquids. Here, we study ternary electrolytes composed of LiTFSI, water, and imidazolium ionic liquids. We find that the LiTFSI solubility strongly increases from 21 mol kg^-1 in water to up to 60 mol kg^-1 in the presence of ionic liquid. The solution structure is investigated with Raman and NMR spectroscopy and the enhanced LiTFSI solubility is found to originate from a hydrotropic effect of the ionic liquids. The increased reductive stability of the ternary electrolytes enables stable cycling of an aqueous lithium-ion battery with an energy density of 150 Wh kg^-1 on the active material level based on commercially relevant Li₄ Ti₅ O₁₂ and LiNi_0.8 Mn_0.1 Co_0.1 O₂ electrode materials.

Spatial-Decomposition Analysis of Electrical Conductivity in Mixtures of Ionic Liquid and Sodium Salt for Sodium-Ion Battery Electrolytes

Ionic liquid (IL)-based electrolytes are a promising material for the development of sodium-ion batteries, and their performance can be quantified by electrical conductivity. In this highly concentrated ionic system, the correlated motions of ion pairs are influential on the ionic transport properties. Herein, all-atom analyses are conducted through molecular dynamics simulations to bridge the macroscopically observable electrical conductivity with the molecular pictures of correlated motion of ion pairs. The analysis is applied to three mixtures of IL with sodium salt that are relevant to the electrolyte for a sodium-ion battery: [1-ethyl-3-methylimidazolium, Na][bis(fluorosulfonyl)amide] ([C₂C₁im, Na][FSA]), [N-methyl-N-propylpyrrolidinium, Na][FSA] ([C₃C₁pyrr, Na][FSA]), and [K, Na][FSA]. The computational results on electrical conductivities are in agreement with the experimental reports, and their dependency on temperature and sodium-ion composition is reproduced well. The overall contributions from cross-correlated motions are found to be negative in all the IL mixtures; thus, the total conductivities are less than their Nernst-Einstein estimates. The spatial view of cross-correlated motions is further obtained by decomposing the time correlation functions of velocities according to the distances between ion pairs. It is observed that ion pairs are moving in the same direction for ∼0.3 ps when they were initially within the first coordination shell, followed by motions toward opposite directions. The cross-correlation terms are also dissected into local and nonlocal components, and it is commonly seen for all the ion pairs that the local component is negative for cation-anion pairs and is positive for cation-cation and anion-anion pairs. The motions of ion pairs are accompanied by a "backflow" that manifests in the form of the nonlocal component whose sign is opposite to the corresponding local component. In fact, the contributions of the correlated motions of ions to the electrical conductivity are not localized to contact pairs and extend spatially beyond the first coordination shell of the cation-anion pairs.

Hopping in High Concentration Electrolytes - Long Time Bulk and Single-Particle Signatures, Free Energy Barriers, and Structural Insights

Although ion-hopping is believed to be a significant mode of transport for small ions in liquid high concentration electrolytes (HCE), its bulk signatures over sufficiently long time intervals are yet to be shown. We computationally establish the long and short time imprints of hopping in HCEs using LiBF₄-in-sulfolane mixtures as models. The high viscosity of this electrolyte leads to significant dynamic heterogeneity in Li-ion transport. Li-ions exhibit a preference to transit to previously occupied Li-ion-sites, bridged through anion or solvent molecules. Hopping in the liquid matrix was found to be an activated process, whose free energy barrier and transition state structure have been determined. Evidence for nanoscale compositional heterogeneity at high salt concentrations is also presented. The simulations shed light on the composition, stiffness, and lifetime of the solvation shell of Li ions. The understanding of HCEs gleaned from this study will spearhead the choice, engineering and applicability of this class of electrolytes.

Effect of salt concentration on properties of mixed carbonate-based electrolyte for Li-ion batteries: a molecular dynamics simulation study

In this work, a computational framework is proposed by utilizing molecular dynamics simulation to explore the existing relation between molecular structure and ionic conductivity of the electrolyte system [LiPF₆+(EC+DMC 1:1)] consisting of a mixture of cyclic ethylene carbonate (EC) and acyclic dimethyl carbonate (DMC) solvents and lithium hexafluorophosphate (LiPF₆) salt to propose as a novel mixed organic solvent-based electrolytes to promote the performance of lithium-ion batteries (LIBs). To acquire a clear understanding of the structural and transport properties of the designed electrolytes, quantum chemistry (QC) calculations and molecular dynamics (MD) simulation are used. In the first step, the accurate molecular structures of the studied electrolytes in addition to their corresponding atomic partial charges are evaluated. The MD simulations are performed at 330 K varying the LiPF₆ concentration (0.5 M to 2.2 M). Analysis of the obtained results indicated that ionic diffusivity and conductivity of the electrolytes are dependent on the structure of solvated ions and lithium salt (LiPF₆) concentration. It is found that the obtained MD simulation results are in reasonable agreement with experimental results. Graphical abstract A representation of dependence of transport properties of electrolyte system [LiPF6 +(EC+DMC 1:1)] as function of salt concentration to be used in Lithium-ion batteries (LIBs).

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Impact of Anion Asymmetry on Local Structure and Supercooling Behavior of Water-in-Salt Electrolytes

Salts with asymmetric (fluorosulfonyl)(trifluoromethanesulfonyl)imide (FTFSI) anions have recently been shown to suppress crystallization of water-in-salt electrolytes, enabling low-temperature operation of high-voltage aqueous rechargeable batteries. To clarify the underlying mechanism for the kinetic suppression of crystallization, we investigate the local solution structures and dynamic behaviors of water-in-salt electrolytes based on the asymmetric FTFSI anion and its symmetric anion analogues by Raman spectroscopy and molecular dynamics simulations. We find that monodentate coordination of FTFSI to cations leads to high rotational mobility of the uncoordinated SO₂CF₃ group. We conclude that the peculiar, coordination-dependent, local dynamics in the asymmetric FTFSI anion, manifested by enhanced intramolecular bond rotation, enables the strong supercooling behavior.

Multiscale Lithium-Battery Modeling from Materials to Cells

New experimental technology and theoretical approaches have advanced battery research across length scales ranging from the molecular to the macroscopic. Direct observations of nanoscale phenomena and atomistic simulations have enhanced the understanding of the fundamental electrochemical processes that occur in battery materials. This vast and ever-growing pool of microscopic data brings with it the challenge of isolating crucial performance-decisive physical parameters, an effort that often requires the consideration of intricate interactions across very different length scales and timescales. Effective physics-based battery modeling emphasizes the cross-scale perspective, with the aim of showing how nanoscale physicochemical phenomena affect device performance. This review surveys the methods researchers have used to bridge the gap between the nanoscale and the macroscale. We highlight the modeling of properties or phenomena that have direct and considerable impact on battery performance metrics, such as open-circuit voltage and charge/discharge overpotentials. Particular emphasis is given to thermodynamically rigorous multiphysics models that incorporate coupling between materials' mechanical and electrochemical states.

Quantitative Investigation of Ion Clusters in a Double Salt Ionic Liquid by Both Vibrational Spectroscopy and Molecular Dynamics Simulation

Vibrational spectroscopy and molecular dynamics simulations are powerful tools frequently used to elucidate interactions among ions in ionic liquid electrolyte solutions. We apply these techniques to characterize ionic interactions in mixtures of 1-butyl-1-methylpyrrolidinium trifluoromethansulfonate, [C₁C₄pyr][CF₃SO₃], and lithium trifluoromethanesulfonate, LiCF₃SO₃, namely, [Li]_0.091[C₁C₄pyr]_0.909[CF₃SO₃] and [Li]_0.167[C₁C₄pyr]_0.833[CF₃SO₃]. The computational and experimental data indicate that extensive, LiCF₃SO₃-rich regions exist within the solutions, and most of the anionic species that are composed of these domains are either [Li₂CF₃SO₃]⁺ or LiCF₃SO₃ moieties. The [Li]_0.167[C₁C₄pyr]_0.833[CF₃SO₃] system contains a larger number of [Li₂CF₃SO₃]⁺ and [Li₃CF₃SO₃]²⁺ species than [Li]_0.091[C₁C₄pyr]_0.909[CF₃SO₃], which may explain, in part, the reduction in ionic conductivity when LiCF₃SO₃ is added to [C₁C₄pyr][CF₃SO₃]. The charge-organized liquid structure inherent to [C₁C₄pyr][CF₃SO₃] supports the dynamic coupling of vibrationally induced dipole moments to form optical phonons. Consequently, intense, IR-active vibrational modes are split into transverse optical and longitudinal optical components. Band splitting is reduced when LiCF₃SO₃ is added to the ionic liquid, suggesting that ionically associated anions impede the ability of the ionic liquid to support optical phonons.

The Use of Succinonitrile as an Electrolyte Additive for Composite-Fiber Membranes in Lithium-Ion Batteries

In the present work, the effect of temperature and additives on the ionic conductivity of mixed organic/ionic liquid electrolytes (MOILEs) was investigated by conducting galvanostatic charge/discharge and ionic conductivity experiments. The mixed electrolyte is based on the ionic liquid (IL) (EMI/TFSI/LiTFSI) and organic solvents EC/DMC (1:1 v/v). The effect of electrolyte type on the electrochemical performance of a LiCoO₂ cathode and a SnO₂/C composite anode in lithium anode (or cathode) half-cells was also investigated. The results demonstrated that the addition of 5 wt.% succinonitrile (SN) resulted in enhanced ionic conductivity of a 60% EMI-TFSI 40% EC/DMC MOILE from ~14 mS·cm^-1 to ~26 mS·cm^-1 at room temperature. Additionally, at a temperature of 100 °C, an increase in ionic conductivity from ~38 to ~69 mS·cm^-1 was observed for the MOILE with 5 wt% SN. The improvement in the ionic conductivity is attributed to the high polarity of SN and its ability to dissolve various types of salts such as LiTFSI. The galvanostatic charge/discharge results showed that the LiCoO₂ cathode with the MOILE (without SN) exhibited a 39% specific capacity loss at the 50th cycle while the LiCoO₂ cathode in the MOILE with 5 wt.% SN showed a decrease in specific capacity of only 14%. The addition of 5 wt.% SN to the MOILE with a SnO₂/C composite-fiber anode resulted in improved cycling performance and rate capability of the SnO₂/C composite-membrane anode in lithium anode half-cells. Based on the results reported in this work, a new avenue and promising outcome for the future use of MOILEs with SN in lithium-ion batteries (LIBs) can be opened.

Stabilizing lithium metal anode by octaphenyl polyoxyethylene-lithium complexation

Lithium metal is an ideal anode for lithium batteries due to its low electrochemical potential and high theoretical capacity. However, safety issues arising from lithium dendrite growth have significantly reduced the practical applicability of lithium metal batteries. Here, we report the addition of octaphenyl polyoxyethylene as an electrolyte additive to enable a stable complex layer on the surface of the lithium anode. This surface layer not only promotes uniform lithium deposition, but also facilitates the formation of a robust solid-electrolyte interface film comprising cross-linked polymer. As a result, lithium|lithium symmetric cells constructed using the octaphenyl polyoxyethylene additive exhibit excellent cycling stability over 400 cycles at 1 mA cm⁻², and outstanding rate performance up to 4 mA cm⁻². Full cells assembled with a LiFePO₄ cathode exhibit high rate capability and impressive cyclability, with capacity decay of only 0.023% per cycle.

Despite the large theoretical promise of Li metal anode, the dendrite growth poses a serious safety challenge. Here the authors address this issue by adding octaphenyl polyoxyethylene as an electrolyte additive which facilitates the formation of a dual-functional layer and excellent performance.

Current Status of AMOEBA-IL: A Multipolar/Polarizable Force Field for Ionic Liquids

Computational simulations of ionic liquid solutions have become a useful tool to investigate various physical, chemical and catalytic properties of systems involving these solvents. Classical molecular dynamics and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of IL systems have provided significant insights at the atomic level. Here, we present a review of the development and application of the multipolar and polarizable force field AMOEBA for ionic liquid systems, termed AMOEBA–IL. The parametrization approach for AMOEBA–IL relies on the reproduction of total quantum mechanical (QM) intermolecular interaction energies and QM energy decomposition analysis. This approach has been used to develop parameters for imidazolium– and pyrrolidinium–based ILs coupled with various inorganic anions. AMOEBA–IL has been used to investigate and predict the properties of a variety of systems including neat ILs and IL mixtures, water exchange reactions on lanthanide ions in IL mixtures, IL–based liquid–liquid extraction, and effects of ILs on an aniline protection reaction.

Unraveling the solvation geometries of the lanthanum(iii) bistriflimide salt in ionic liquid/acetonitrile mixtures

La(Tf₂N)₃ in C₈(mim)₂(Tf₂N)₂/acetonitrile mixtures forms 10-fold coordination complexes composed of both acetonitrile molecules and Tf₂N⁻ anions.