Effect of ion structure on conductivity in lithium-doped ionic liquid electrolytes: A molecular dynamics study

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.

Insights into the intermolecular interactions and temperature-concentration dependence of transport in ionic liquid-based EMI–TFSI/LiTFSI electrolytes

DFT calculations and MD simulations show that the EMI–TFSI/LiTFSI system is stabilized by strong nonbonded attractions, and the lithium-ion conductivity depends on the LiTFSI concentration and system temperature.

Spectroscopic analysis focusing on ionic liquid/metal electrode and organic semiconductor interfaces in an electrochemical environment

The solid-liquid interface forms an electric double layer that enables the function of electronic devices and, thus, represents an important area of electrochemical research. Because ionic liquids (ILs) are becoming prominent candidates for new high-performing electrolytes, their interface with solid substrates (e.g., metal electrodes or organic semiconductors) attracts substantial attention. An example of improvement achieved using ILs as electrolytes is a decrease in the operating voltage of transistors from >10 V in traditional SiO₂-gated transistors to <1 V in IL-gated electronic double-layer organic field-effect devices. This perspective discusses the investigation of poorly accessible IL/substrate interfaces using both attenuated total reflectance ultraviolet (ATR-UV) spectroscopy and a newly developed electrochemical setup combined with ATR-UV (EC-ATR-UV), which allows analysis of the interfacial area under the application of varying electric potential. The recent EC-ATR-UV applications in interfacial analytical chemistry are overviewed and compared to other spectroscopic methods described in the recent literature. Lastly, the supplementation of experimental data with theoretical calculations (e.g., quantum chemical calculations and molecular dynamics simulations) is also addressed.

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.

Interplay of Structure and Dynamics in Lithium/Ionic Liquid Electrolytes: Experiment and Molecular Simulation

Despite their promising use in electrochemical and electrokinetic devices, ionic-liquid-based electrolytes often exhibit complex behavior arising from a subtle interplay of their structure and dynamics. Here, we report a joint experimental and molecular simulation study of such electrolytes obtained by mixing 1-butyl 3-methylimidazolium tetrafluoroborate with lithium tetrafluoroborate. More in detail, experiments consisting of X-ray scattering, pulsed field gradient NMR, and complex impedance spectroscopy are analyzed in the light of molecular dynamics simulations to probe the structural, dynamical, and electrochemical properties of this ionic-liquid-based electrolyte. Lithium addition promotes the nanostructuration of the liquid as evidenced from the appearance of a scattering prepeak that becomes more pronounced. Microscopically, using the partial structure factors determined from molecular dynamics, this prepeak is shown to correspond to the formation of well-ordered positive/negative charge series and also large aggregates (Li_n(BF₄)_4-m)^(4-m+n)-, which develop upon lithium addition. Such nanoscale ordering entails a drastic decrease in both the molecular mobility and ionic conductivity. In particular, the marked association of Li⁺ cations with four BF₄^- anions and long ion pairing times, which are promoted upon lithium addition, are found to severely hinder the Li⁺ transport properties.

Diffusivity and Structure of Room Temperature Ionic Liquid in Various Organic Solvents

Room-temperature ionic liquids (RTILs) hold promise for applications in electric double layer capacitors (EDLCs), owing to a much wider potential window, lower vapor pressure, and better thermal and chemical stabilities compared to conventional aqueous and organic electrolytes. However, because the low diffusivity of ions in neat RTILs negates the EDLCs' advantage of high power density, the ionic liquids are often used in mixture with organic solvents. In this study, we measured the diffusivity of cations and anions in RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) ([BMIM⁺][TFSI^-]), mixed with 10 organic solvents, by using the pulsed-field gradient NMR method. The ion diffusivity was found to follow that of neat solvents and in most studied solvents showed an excellent agreement with the predicted values reported in the recent molecular dynamics (MD) study [Thompson, M. W.; J. Phys. Chem. B 2019, 123, 1340-1347]. In two solvents consisting of long-chain molecules, however, the MD simulations predictions slightly underestimated the ionic diffusivities. The degree of ion dissociation was also estimated for each solvent by comparing the ionic conductivity with the molar conductivity derived from the diffusion measurements. The degree of ion dissociation and the hydrodynamic radius of ions suggest that the ions are coordinated by ∼1 solvent molecule. The scarcity of solvent-ion interactions explains the fact that the diffusivity of ions in the mixture significantly depends on the viscosity of the solvent.

Interfacial Structures in Ionic Liquid-Based Ternary Electrolytes for Lithium-Metal Batteries: A Molecular Dynamics Study

Lithium-metal batteries are promising candidates to fulfill the future performance requirements for energy storage applications. However, the tendency to form metallic dendrites and the undesirable side reactions between the electrolyte and the Li electrode lead to poor performance and safety issues in these batteries. Therefore, understanding the interfacial properties and the Li-metal surface/electrolyte interactions is crucial to resolve the remaining obstacles and make these devices feasible. Here, we report a computational study on the interface effects in ternary polymer electrolytes composed by poly(ethylene oxide) (PEO), lithium salts, and different ionic liquids (ILs) confined between two Li-metal slabs. Atomistic simulations are used to characterize the local environment of the Li⁺ ions and the transport properties in the bulk and at the interface regions. Aggregation of ions at the metal surface is seen in all investigated systems; the structure and composition are directly correlated to the IL components. The strong interactions between the electrolyte species and the Li-metal atoms result in the structuring of the electrolyte at the interface region, in which comparatively small and flat ions result in a well-defined region with extensive Li⁺ populations and high self-diffusion coefficients. In contrast, large ions such as [P222mom]⁺ increase the PEO density in the bulk due to large steric effects at the interface. Therefore, the choice of specific ILs in ternary polymer electrolytes can tune the structure-dynamic properties at the Li-metal surface/electrolyte interface, controlling the SEI formation at the electrode surface, and thereby improve battery performance.

Self-Learning Molecular Design for High Lithium-Ion Conductive Ionic Liquids using Maze Game

A self-learning artificial intelligence system for an autonomous molecular search was recently utilized in place of laborious material development processes by humans. In this approach, because the evaluation of unsuitable or unrealistic candidates considerably decreases the search efficiency, prior knowledge of the chemistry and engineering requirements should be embedded into the molecular-generative algorithm. However, when using naive rule-based restrictions, one must implement the complex rule logic into the code each time, depending on the materials and potential applications. Herein, we propose a molecular-generative method using a maze game to control the allowable constituent fragments of molecules, which improves the flexibility and consistency to implement the rules. We performed an autonomous search for optimized cation structures of high Li-ion conductive ionic liquids evaluated by molecular dynamics simulations, in its practically reasonable scope defined by the maze game. From the search, we discover that acyl ammonium cations are favorable for high Li-ion conductivity because of the high association between the cations and Li ions. These results broaden our existing insight owing to the ability to explore beyond our practical experiences.

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.

Ultrafast vibrational dynamics of a trigonal planar anionic probe in ionic liquids (ILs): A two-dimensional infrared (2DIR) spectroscopic investigation

A major impediment limiting the widespread application of ionic liquids (ILs) is their high shear viscosity. Incorporation of a tricyanomethanide (TCM^-) anion in ILs leads to low shear viscosity and improvement of several characteristics suitable for large scale applications. However, properties including interactions of TCM^- with the local environment and dynamics of TCM^- have not been thoroughly investigated. Herein, we have studied the ultrafast dynamics of TCM^- in several imidazolium ILs using linear IR and two-dimensional infrared spectroscopy techniques. The spectral diffusion dynamics of the CN stretching modes of TCM^- in all ILs exhibit a nonexponential behavior with a short time component of ∼2 ps and a long time component spanning ∼9 ps to 14 ps. The TCM^- vibrational probe reports a significantly faster relaxation of ILs compared to those observed previously using linear vibrational probes, such as thiocyanate and selenocyanate. Our results indicate a rapid relaxation of the local ion-cage structure embedding the vibrational probe in the ILs. The faster relaxation suggests that the lifetime of the local ion-cage structure decreases in the presence of TCM^- in the ILs. Linear IR spectroscopic results show that the hydrogen-bonding interaction between TCM^- and imidazolium cations in ILs is much weaker. Shorter ion-cage lifetimes together with weaker hydrogen-bonding interactions account for the low shear viscosity of TCM^- based ILs compared to commonly used ILs. In addition, this study demonstrates that TCM^- can be used as a potential vibrational reporter to study the structure and dynamics of ILs and other molecular systems.

Restricted Ion Transport by Plasticizing Side Chains in Polycarbonate-Based Solid Electrolytes

Increasing the ionic conductivity has for decades been an overriding goal in the development of solid polymer electrolytes. According to fundamental theories on ion transport mechanisms in polymers, the ionic conductivity is strongly correlated to free volume and segmental mobility of the polymer for the conventional transport processes. Therefore, incorporating plasticizing side chains onto the main chain of the polymer host often appears as a clear-cut strategy to improve the ionic conductivity of the system through lowering of the glass transition temperature (T g). This intended correlation between T g and ionic conductivity is, however, not consistently observed in practice. The aim of this study is therefore to elucidate this interplay between segmental mobility and polymer structure in polymer electrolyte systems comprising plasticizing side chains. To this end, we utilize the synthetic versatility of the ion-conductive poly(trimethylene carbonate) (PTMC) platform. Two types of host polymers with side chains added to a PTMC backbone are employed, and the resulting electrolytes are investigated together with the side chain-free analogue both by experiment and with molecular dynamics (MD) simulations. The results show that while added side chains do indeed lead to a lower T g, the total ionic conductivity is highest in the host matrix without side chains. It was seen in the MD simulations that while side chains promote ionic mobility associated with the polymer chain, the more efficient interchain hopping transport mechanism occurs with a higher probability in the system without side chains. This is connected to a significantly higher solvation site diversity for the Li+ ions in the side-chain-free system, providing better conduction paths. These results strongly indicate that the side chains in fact restrict the mobility of the Li+ ions in the polymer hosts.

The Effect of Concentration of Lithium Salt on the Structural and Transport Properties of Ionic Liquid-Based Electrolytes

Ionic liquids (ILs) are used as electrolytes in high-performance lithium-ion batteries, which can effectively improve battery safety and energy storage capacity. All atom molecular dynamics simulation and experiment were combined to investigate the effect of the concentration of lithium salt on the performance of electrolytes of four IL solvents ([C_nmim][TFSI] and [C_nmim][FSI], n = 2, 4). The IL electrolytes exhibit higher density and viscosity; meanwhile, larger lithium ion transfer numbers as the concentration of LiTFSI increases. Furthermore, in order to explore the effect of the concentration of lithium salt on the ionic associations of Li⁺ and anion of IL, the microstructures of the lithium salt in various IL electrolytes at different concentrations were investigated. The structural analysis indicated that strong bidentate and monodentate coordination was found between Li⁺ and anion of all IL electrolytes. Both cis and trans isomerism of [FSI]⁻ were observed in [FSI]⁻-type IL electrolyte systems. Furthermore, the existence of the ion cluster [Li[anion]_x]^(x−1)− in IL electrolytes and the cluster became more closed and compact as the concentration of LiTFSI increases.

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.

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.

Attenuated total reflectance far-ultraviolet and deep-ultraviolet spectroscopy analysis of the electronic structure of a dicyanamide-based ionic liquid with Li+

Elucidating the unique electronic structure of ionic liquid molecules around Li⁺ using electronic absorption spectroscopy, theoretical calculations, and chemometric analyses.

Cation influence on heterocyclic ammonium ionic liquids: a molecular dynamics study

Four different ionic liquids (ILs) consisting of the bis(trifluoromethanesulfonyl)imide ([NTf₂]⁻) anion, with structurally similar systematically varying cations, are investigated herein through classical molecular dynamics.

Fluctuation enhancement of ion diffusivity in liquids

The diffusivity of ions in liquid solutions is known either to decrease with an increase in the ion size or to have a single maximum depending on the ion size. This article presents evidence for the appearance of multiple maxima and thus multiple ion sizes with enhanced diffusivity.

Dissolved chloride markedly changes the nanostructure of the protic ionic liquids propylammonium and ethanolammonium nitrate

The bulk nanostructure of 15 mol% propylammonium chloride (PACl) dissolved in propylammonium nitrate (PAN) and 15 mol% ethanolammonium chloride (EtACl) in ethanolammonium nitrate (EtAN) has been determined using neutron diffraction with empirical potential structure refinement fits.

Ab Initio Simulations and Electronic Structure of Lithium-Doped Ionic Liquids: Structure, Transport, and Electrochemical Stability

Density functional theory (DFT), density functional theory molecular dynamics (DFT-MD), and classical molecular dynamics using polarizable force fields (PFF-MD) are employed to evaluate the influence of Li(+) on the structure, transport, and electrochemical stability of three potential ionic liquid electrolytes: N-methyl-N-butylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([pyr14][TFSI]), N-methyl-N-propylpyrrolidinium bis(fluorosulfonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]). We characterize the Li(+) solvation shell through DFT computations of [Li(Anion)n]((n-1)-) clusters, DFT-MD simulations of isolated Li(+) in small ionic liquid systems, and PFF-MD simulations with high Li-doping levels in large ionic liquid systems. At low levels of Li-salt doping, highly stable solvation shells having two to three anions are seen in both [pyr14][TFSI] and [pyr13][FSI], whereas solvation shells with four anions dominate in [EMIM][BF4]. At higher levels of doping, we find the formation of complex Li-network structures that increase the frequency of four anion-coordinated solvation shells. A comparison of computational and experimental Raman spectra for a wide range of [Li(Anion)n]((n-1)-) clusters shows that our proposed structures are consistent with experiment. We then compute the ion diffusion coefficients and find measures from small-cell DFT-MD simulations to be the correct order of magnitude, but influenced by small system size and short simulation length. Correcting for these errors with complementary PFF-MD simulations, we find DFT-MD measures to be in close agreement with experiment. Finally, we compute electrochemical windows from DFT computations on isolated ions, interacting cation/anion pairs, and liquid-phase systems with Li-doping. For the molecular-level computations, we generally find the difference between ionization energy and electron affinity from isolated ions and interacting cation/anion pairs to provide upper and lower bounds, respectively, to experiment. In the liquid phase, we find the difference between the lowest unoccupied and highest occupied electronic levels in pure and hybrid functionals to provide lower and upper bounds, respectively, to experiment. Li-doping in the liquid-phase systems results in electrochemical windows little changed from the neat systems.

Direct Correlation between Ionic Liquid Transport Properties and Ion Pair Lifetimes: A Molecular Dynamics Study

Self-diffusivities as a function of temperature were computed for 29 different ionic liquids (ILs) covering a wide variety of cation and anion classes. Ideal ionic conductivities (σNE) were estimated from the self-diffusivities via the Nernst-Einstein relation. The ion pair (IP) lifetimes (τIP) and ion cage (IC) lifetimes (τIC) of each IL were also computed. A linear relationship between the calculated self-diffusivities and the inverse of IP or IC lifetimes was observed. A similar inverse linear relationship was also observed for ideal ionic conductivity. These relationships were found to be independent of temperature and the nature of the IL. These observations connect macroscopic dynamic properties with local atomic-level motions and strongly suggest that the dynamics of ILs are governed by a universal IP or IC forming and breaking mechanism. Thus, in order to design an ionic liquid with enhanced dynamics, one should consider how to minimize IP or IC lifetimes.

Structural and aggregate analyses of (Li salt + glyme) mixtures: the complex nature of solvate ionic liquids

The structure and interactions of different (lithium salt plus glyme) equimolar mixtures are probed by Molecular Dynamics simulations.

Refined method for predicting electrochemical windows of ionic liquids and experimental validation studies

A combined classical molecular dynamics (MD) and ab initio MD (AIMD) method was developed for the calculation of electrochemical windows (ECWs) of ionic liquids. In the method, the liquid phase of ionic liquid is explicitly sampled using classical MD. The electrochemical window, estimated by the energy difference between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), is calculated at the density functional theory (DFT) level based on snapshots obtained from classical MD trajectories. The snapshots were relaxed using AIMD and quenched to their local energy minima, which assures that the HOMO/LUMO calculations are based on stable configurations on the same potential energy surface. The new procedure was applied to a group of ionic liquids for which the ECWs were also experimentally measured in a self-consistent manner. It was found that the predicted ECWs not only agree with the experimental trend very well but also the values are quantitatively accurate. The proposed method provides an efficient way to compare ECWs of ionic liquids in the same context, which has been difficult in experiments or simulation due to the fact that ECW values sensitively depend on experimental setup and conditions.