Molecular Dynamics Simulations of Lithium-Doped Ionic-Liquid Electrolytes

Voltage-induced modulation of interfacial ionic liquids measured using surface plasmon resonant grating nanostructures

We have used surface plasmon resonant metal gratings to induce and probe the dielectric response (i.e., electro-optic modulation) of ionic liquids (ILs) at electrode interfaces. Here, the cross-plane electric field at the electrode surface modulates the refractive index of the IL due to the Pockels effect. This is observed as a shift in the resonant angle of the grating (i.e., Δϕ), which can be related to the change in the local index of refraction of the electrolyte (i.e., Δnlocal). The reflection modulation of the IL is compared against a polar (D2O) and a non-polar solvent (benzene) to confirm the electro-optic origin of resonance shift. The electrostatic accumulation of ions from the IL induces local index changes to the gratings over the extent of electrical double layer (EDL) thickness. Finite difference time domain simulations are used to relate the observed shifts in the plasmon resonance and change in reflection to the change in the local index of refraction of the electrolyte and the thickness of the EDL. Simultaneously using the wavelength and intensity shift of the resonance enables us to determine both the effective thickness and Δn of the double layer. We believe that this technique can be used more broadly, allowing the dynamics associated with the potential-induced ordering and rearrangement of ionic species in electrode-solution interfaces.

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.

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.

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.

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

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

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.

Ether-Functionalized Pyrrolidinium-Based Room Temperature Ionic Liquids: Physicochemical Properties, Molecular Dynamics, and the Lithium Ion Coordination Environment

The physicochemical properties of room temperature ionic liquids (RTILs) consisting of bis(trifluoromethanesulfonyl)amide (TFSA^- ) combined with 1-hexyl-1-methylpyrrolidinium (Pyr_1,6 ⁺ ), 1-(butoxymethyl)-1-methylpyrrolidinium (Pyr_1,1O4 ⁺ ), 1-(4-methoxybutyl)-1-methyl pyrrolidinium (Pyr_1,4O1 ⁺ ), and 1-((2-methoxyethoxy)methyl)-1-methylpyrrolidinium (Pyr_1,1O2O1 ⁺ ) were investigated using both experimental and computational approaches. Pyr_1,1O2O1 TFSA, which contains two ether oxygen atoms, showed the lowest viscosity, and the relationship between its physicochemical properties and the position and number of the ether oxygen atoms was discussed by a careful comparison with Pyr_1,1O4 TFSA and Pyr_1,4O1 TFSA. Ab initio calculations revealed the conformational flexibility of the side chain containing the ether oxygen atoms. In addition, molecular dynamics (MD) calculations suggested that the ion distributions have a significant impact on the transport properties. Furthermore, the coordination environments of the Li ions in the RTILs were evaluated using Raman spectroscopy, which was supported by MD calculations using 1000 ion pairs. The presented results will be valuable for the design of functionalized RTILs for various applications.

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.

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.

Understanding Structural and Transport Properties of Dissolved Li2 S8 in Ionic Liquid Electrolytes through Molecular Dynamics Simulations

Lithium-sulfur batteries with high energy density are considered as one of the most promising future energy storage devices. However, the parasitic lithium polysulfides shuttle phenomenon severely hinders the commercialization of such batteries. Ionic liquids have been found to suppress the lithium polysulfides solubility, diminishing the shuttle effect effectively. Herein, we performed classical molecular dynamics simulations to explore the microscopic mechanism and transport behaviors of typical Li₂ S₈ species in ionic liquids and ionic liquid-based electrolyte systems. We found that the trifluoromethanesulfonate anions ([OTf]^- ) exhibit higher coordination strength with lithium ions compared with bis(trifluoromethanesulfonyl)imide anions ([TFSI]^- ) in static microstructures. However, the dynamical characteristics indicate that the presence of the [OTf]^- anions in ionic liquid electrolytes bring faster Li⁺ exchange rate and easier dissociation of Li⁺ solvation structures. Our simulation models offer a significant guidance to future studies on designing ionic liquid electrolytes for lithium-sulfur batteries.

Pseudohalogen Chemistry in Ionic Liquids with Non-innocent Cations and Anions

Within the second funding period of the SPP 1708 "Material Synthesis near Room Temperature",which started in 2017, we were able to synthesize novel anionic species utilizing Ionic Liquids (ILs) both, as reaction media and reactant. ILs, bearing the decomposable and non-innocent methyl carbonate anion [CO3 Me]- , served as starting material and enabled facile access to pseudohalide salts by reaction with Me3 Si-X (X=CN, N3 , OCN, SCN). Starting with the synthesized Room temperature Ionic Liquid (RT-IL) [nBu3 MeN][B(OMe)3 (CN)], we were able to crystallize the double salt [nBu3 MeN]2 [B(OMe)3 (CN)](CN). Furthermore, we studied the reaction of [WCC]SCN and [WCC]CN (WCC=weakly coordinating cation) with their corresponding protic acids HX (X=SCN, CN), which resulted in formation of [H(NCS)2 ]- and the temperature labile solvate anions [CN(HCN)n ]- (n=2, 3). In addition, the highly labile anionic HCN solvates were obtained from [PPN]X ([PPN]=μ-nitridobis(triphenylphosphonium), X=N3 , OCN, SCN and OCP) and HCN. Crystals of [PPN][X(HCN)3 ] (X=N3 , OCN) and [PPN][SCN(HCN)2 ] were obtained when the crystallization was carried out at low temperatures. Interestingly, reaction of [PPN]OCP with HCN was noticed, which led to the formation of [P(CN)2 ]- , crystallizing as HCN disolvate [PPN][P(CN⋅HCN)2 ]. Furthermore, we were able to isolate the novel cyanido(halido) silicate dianions of the type [SiCl0.78 (CN)5.22 ]2- and [SiF(CN)5 ]2- and the hexa-substituted [Si(CN)6 ]2- by temperature controlled halide/cyanide exchange reactions. By facile neutralization reactions with the non-innocent cation of [Et3 HN]2 [Si(CN)6 ] with MOH (M=Li, K), Li2 [Si(CN)6 ] ⋅ 2 H2 O and K2 [Si(CN)6 ] were obtained, which form three dimensional coordination polymers. From salt metathesis processes of M2 [Si(CN)6 ] with different imidazolium bromides, we were able to isolate new imidazolium salts and the ionic liquid [BMIm]2 [Si(CN)6 ]. When reacting [Mes(nBu)Im]2 [Si(CN)6 ] with an excess of the strong Lewis acid B(C6 F5 )3 , the voluminous adduct anion {Si[CN⋅B(C6 F5 )3 ]6 }2- was obtained.

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.

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.

Insights into the solvation and dynamic behaviors of a lithium salt in organic- and ionic liquid-based electrolytes

Fundamental molecular insights were provided to understand the advantages of IL solvent electrolytes with high conductivity over organic solvent electrolytes.