Heterogeneous Li coordination in solvent-in-salt electrolytes enables high Li transference numbers

The transport properties and the underlying coordination structure of a ternary electrolyte consisting of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), 1,2-dimethoxyethane (DME), and 1,3-dioxolane (DOL) is studied over a wide concentration range, up to that of a Solvent-in-Salt (SiS) electrolyte. Among other advantages for next-generation battery applications, SiS electrolytes offer a high lithium transference number (tLi) of 0.73. We analyze the transport mechanism by electrophoretic NMR (eNMR), providing the mobilities (μi) of all species. Intriguingly, in the SiS region, the mobility of the neutral species DME exceeds the cation mobility (μDME > μLi), suggesting a heterogeneous transport mechanism, where the Li+ mobility is averaged over different species. Based on Raman spectroscopy, NMR spectroscopy and MD simulations, we derive a model for a concentration-dependent Li+ coordination environment with a heterogeneous Li+ coordination in the SiS region, where the 1st coordination shell either consists of TFSI- and DOL only, or of DME, TFSI-, and DOL. Lithium ions partially coordinated by DME migrate faster in an electric field, in contrast to lithium ions solely coordinated by anions and DOL molecules, explaining the peculiarity of the rapidly migrating neutral DME molecules. Further, DME is identified as an exclusively bidentate ligand, while TFSI- and DOL act as bridging ligands coordinating different Li+ ions. Thus, Li+ coordination heterogeneity is the basis for Li+ transport heterogeneity and for achieving very high Li+ transference numbers. In addition, an effective dynamic decoupling of Li+ and anions occurs with an Onsager coefficient σ+- ≈ 0. These results provide a deeper understanding of the very efficient lithium-ion transport in SiS electrolytes, with the potential to bring further improvements for battery applications.

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li+) and bis(trifluoromethanesulfonyl)imide (TFSI-) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force-torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion-ion and ion-solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1-4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Overall, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

A Computational and Spectroscopic Analysis of Solvate Ionic Liquids Containing Anions with Long and Short Perfluorinated Alkyl Chains

Anion-driven, nanoscale polar-apolar structural organization is investigated in a solvate ionic liquid (SIL) setting by comparing sulfonate-based anions with long and short perfluorinated alkyl chains. Representative SILs are created from 1,2-bis(2-methoxyethoxy)ethane ("triglyme" or "G3"), lithium nonafluoro-1-butanesulfonate, and lithium trifluoromethanesulfonate. Molecular dynamics simulations, density functional theory computations, and vibrational spectroscopy provide insight into the overall liquid structure, cation-solvent interactions, and cation-anion association. Significant competition between G3 and anions for cation-binding sites characterizes the G3-LiC4F9SO3 mixtures. Only 50% of coordinating G3 molecules form tetradentate complexes with Li+ in [(G3)1Li][C4F9SO3]. Moreover, the SIL is characterized by extensive amounts of ion pairing. Based on these observations, [(G3)1Li][C4F9SO3] is classified as a "poor" SIL, similar to the analogous [(G3)1Li][CF3SO3] system. Even though the comparable basicity of the CF3SO3- and C4F9SO3- anions leads to similar SIL classifications, the hydrophobic fluorobutyl groups support extensive apolar domain formation. These apolar moieties permeate throughout [(G3)1Li][C4F9SO3] and persist even at relatively low dilution ratios of [(G3)10Li][C4F9SO3]. By way of comparison, the CF3 group is far too short to sustain polar-apolar segregation. This demonstrates how chemically modifying the anions to include hydrophobic groups can impart unique nanoscale organization to a SIL. Moreover, tuning these nano-segregated fluorinated domains could, in principle, control the presence of dimensionally ordered states in these mixtures without changing the coordination of the lithium ions.

Accelerated lithium-ion diffusion via a ligand 'hopping' mechanism in lithium enriched solvate ionic liquids

Solvate ionic liquids (SILs), equimolar amounts of lithium salts and polyether glymes, are well studied highly customisable "designer solvents". Herein the physical, thermal and ion mobility properties of SILs with increased LiTFSI (LiTFSA) concentration, with ligand 1 : >1 LiTFSI stoichiometric ratios, are presented. It was found that between 60-80 °C, the lithium cation diffuses up to 4 times faster than the corresponding anion or ligand (glyme). These systems varied from viscous liquids to self-supporting gels, though were found to thin exponentially when heated to mild temperatures (50-60 °C). They were also found to be thermally stable, up to 200 °C, well in excess of normal operating temperatures. Ion mobility, assessed under an electric potential via ionic conductivity, showed the benefit of SIL optimisation for attaining greater concentrations of Li+ cations to store charge during supercapacitor charging and discharging. Molecular dynamics simulations interrogate the mechanism of enhanced diffusion at high temperatures, revealing a lithium hopping mechanism that implicates the glyme in bridging two lithiums through changes in the denticity.

Infrared Spectroscopy of Li+ Solvation in Diglyme: Ab Initio Molecular Dynamics and Experiment

Infrared (IR) spectra of solutions of the lithium salt LiBF4 in diglyme, CH3O(CH2CH2O)2CH3, are studied via IR spectroscopy and ab initio molecular dynamics (AIMD) simulations. Experiments show that the major effects of LiBF4, compared to neat diglyme, are the appearance of a new broad band in the 250-500 cm-1 frequency region and a broadening and intensity enhancement of the diglyme band in the 900-1150 cm-1 region accompanied by a red-shift. Computational analysis indicates that hindered translational motions of Li+ in its solvation cage are mainly responsible for the new far-IR band, while the changes in the mid-IR are due to Li+-coordination-dependent B-F stretching vibrations of BF4- anions coupled with diglyme vibrations. Molecular motions in these and lower frequency regions are generally correlated, revealing the collective nature of the vibrational dynamics, which involve multiple ions/molecules. Herein, a detailed analysis of these features via AIMD simulations of the spectrum and its components, combined with analysis of the generalized normal modes of the solution components, is presented. Other minor spectral changes as well as diglyme conformational changes induced by the lithium salt are also discussed.

Accurate Calculation of Solvation Properties of Lithium Ions in Nonaqueous Solutions

We perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

Dynamic Heterogeneity of Solvent Motion and Ion Transport in Concentrated Electrolytes

Molecular-level understanding of the cation transference number t₊⁰, an important property that characterizes the transport of working cations, is critical to the bottom-up design of battery electrolytes. We quantify t₊⁰ in a model tetraglyme-based electrolyte using molecular dynamics simulation and the Onsager approach. t₊⁰ exhibits a concentration dependence in three distinct regimes: dilute, intermediate, and concentrated. The cluster approximation uncovers dominant correlations and dynamic heterogeneity in each regime. In the dilute regime, t₊⁰ decreases sharply as increasing numbers of solvent molecules become coordinated with Li⁺. The crossover to the intermediate regime, t₊⁰ ≈ 0, occurs when all solvent molecules become coordinated, and a plateau is obtained because anions enter the Li⁺ solvation shell, resulting in ion pairs that do not contribute to t₊⁰. Transference in concentrated electrolytes is dominated by the presence of cations in a variety of negatively charged and solvent-excluded clusters, resulting in t₊⁰ < 0.

Polypropylene carbonate-based electrolytes as model for a different approach towards improved ion transport properties for novel electrolytes

Linear poly(alkylene carbonates) such as polyethylene carbonate (PEC) and polypropylene carbonate (PPC) have gained increasing interest due to their remarkable ion transport properties such as high Li⁺ transference numbers. The cause of these properties is not yet fully understood which makes it challenging to replicate them in other polymer electrolytes. Therefore, it is critical to understand the underlying mechanisms in polycarbonate electrolytes such as PPC. In this work we present insights from impedance spectroscopy, transference number measurements, PFG-NMR, IR and Raman spectroscopy as well as molecular dynamics simulations to address this issue. We find that in addition to plasticization, the lithium ion coordination by the carbonate groups of the polymer is weakened upon gelation, leading to a rapid exhange of the lithium ion solvation shell and consequently a strong increase of the conductivity. Moreover, we study the impact of the anions by employing different conducting salts. Interestingly, while the total conductivity decreases with increasing anion size, the reverse trend can be observed for the lithium ion transference numbers. Via our holistic approach, we demonstrate that this behavior can be attributed to differences in the collective ion dynamics.

Electric-Field-Induced Spatially Dynamic Heterogeneity of Solvent Motion and Cation Transference in Electrolytes

While electric fields primarily result in migration of charged species in electrolytic solutions, the solutions are dynamically heterogeneous. Solvent molecules within the solvation shells of the cation will be dragged by the field while free solvent molecules will not. We combine electrophoretic NMR measurements of ion and solvent velocities under applied electric fields with molecular dynamics simulations to interrogate different solvation motifs in a model liquid electrolyte. Measured values of the cation transference number (t_{+}^{0}) agree quantitatively with simulation-based predictions over a range of electrolyte concentrations. Solvent-cation interactions strongly influence the concentration-dependent behavior of t_{+}^{0}. We identify a critical concentration at which most of the solvent molecules lie within solvation shells of the cations. The dynamic heterogeneity of solvent molecules is minimized at this concentration where t_{+}^{0} is approximately equal to zero.

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.

Ionic Liquids and Water: Hydrophobicity vs. Hydrophilicity

Many chemical processes rely extensively on organic solvents posing safety and environmental concerns. For a successful transfer of some of those chemical processes and reactions to aqueous media, agents acting as solubilizers, or phase-modifiers, are of central importance. In the present work, the structure of aqueous solutions of several ionic liquid systems capable of forming multiple solubilizing environments were modeled by molecular dynamics simulations. The effect of small aliphatic chains on solutions of hydrophobic 1-alkyl-3-methylimidazolium bis(trifluoromethyl)sulfonylimide ionic liquids (with alkyl = propyl [C₃C₁im][NTf₂], butyl [C₄C₁im][NTf₂] and isobutyl [iC₄C₁im][NTf₂]) are covered first. Next, we focus on the interactions of sulphonate- and carboxylate-based anions with different hydrogenated and perfluorinated alkyl side chains in solutions of [C₂C₁im][C_nF_2n+1SO₃], [C₂C₁im][C_nH_2n+1SO₃], [C₂C₁im][CF₃CO₂] and [C₂C₁im][CH₃CO₂] (n = 1, 4, 8). The last system considered is an ionic liquid completely miscible with water that combines the cation N-methyl-N,N,N-tris(2-hydroxyethyl)ammonium [N_{1 2OH 2OH 2OH}]⁺, with high hydrogen-bonding capability, and the hydrophobic anion [NTf₂]^-. The interplay between short- and long-range interactions, clustering of alkyl and perfluoroalkyl tails, and hydrogen bonding enables a wealth of possibilities in tailoring an ionic liquid solution according to the needs.

NaFSI and NaTFSI Solutions in Ether Solvents from Monoglyme to Poly(ethylene oxide)-A Molecular Dynamics Study

Classical molecular dynamics simulations have been performed for a series of electrolytes based on sodium bis(fluorosulfonyl)imide or sodium bis(trifluoromethylsulfonyl)imide salts and monoglyme, tetraglyme, and poly(ethylene oxide) as solvents. Structural properties have been assessed through the analysis of coordination numbers and binding patterns. Residence times for Na-O interactions have been used to investigate the stability of solvation shells. Diffusion coefficients of ions and electrical conductivity of the electrolytes have been estimated from molecular dynamics trajectories. Contributions to the total conductivity have been analyzed in order to investigate the role of ion-ion correlations. It has been found that the anion-cation interactions are more probable in the systems with NaTFSI salts. Accordingly, the degree of correlations between ion motions is larger in NaTFSI-based electrolytes.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

MeTFSI (Me = Li, Na) Solvation in Ethylene Carbonate and Fluorinated Ethylene Carbonate: A Molecular Dynamics Study

Classical and ab initio molecular dynamics (MD) simulations have been performed for electrolytes based on LiTFSI and NaTFSI solutions in ethylene carbonate and its mono- and difluoro derivatives. Differences between electrolytes with Li⁺ or Na⁺ ions and the effect of fluorination on the structure and transport properties have been analyzed. The observed differences are related to the strength of Me⁺-carbonate binding, which is weaker for the Na⁺ cation and/or fluorinated solvents. Infrared spectra have been computed from ab initio MD and density functional tight binding (DFTB) MD trajectories. The changes of vibrational frequencies have been related to the local structure of the electrolyte and to interactions between salt cations and solvent molecules. The frequency shifts obtained from the AIMD simulations agree with experimental data, whereas DFTB underestimates Na⁺-carbonate interactions.

Quantification of cation-cation, anion-anion and cation-anion correlations in Li salt/glyme mixtures by combining very-low-frequency impedance spectroscopy with diffusion and electrophoretic NMR

Directional correlations between the movements of cations and anions exert a strong influence on the charge and mass transport properties of concentrated battery electrolytes. Here, we combine, for the first time, very-low-frequency impedance spectroscopy on symmetrical Li|electrolyte|Li cells with diffusion and electrophoretic NMR in order to quantify cation-cation, anion-anion and cation-anion correlations in Li salt/tetraglyme (G4) mixtures with Li salt to G4 ratios between 1 : 1 and 1 : 2. We find that all correlations are negative, with like-ion anticorrelations (cation-cation and anion-anion) being generally stronger than cation-anion anticorrelations. In addition, we observe that like-ion anticorrelations are stronger for the heavier type of ion and that all anticorrelations become weaker with decreasing Li salt to G4 ratio. These findings are in contrast to theories considering exclusively anion-cation correlations in form of ion pairs, as the latter imply positive cation-anion correlations. We analyze in detail the influence of anticorrelations on Li+ transference numbers and on the Haven ratio. In order to rationalize our results, we derive linear response theory expressions for all ion correlations. These expressions show that the Li+ ion transport under anion-blocking conditions in a battery is governed by equilibrium center-of-mass fluctuations in the electrolytes. This suggests that in future electrolyte theories and computer simulations, more attention should be paid to equilibrium center-of-mass fluctuations.

Active learning and neural network potentials accelerate molecular screening of ether-based solvate ionic liquids

Solvate ionic liquids (SIL) have promising applications as electrolyte materials. Despite the broad design space of oligoether ligands, most reported SILs are based on simple tri- and tetraglyme. Here, we describe a computational search for complex ethers that can better stabilize SILs. Through active learning, a neural network interatomic potential is trained from density functional theory data. The learned potential fulfills two key requirements: transferability across composition space, and high speed and accuracy to find low-energy ligand-ion poses across configurational space. Candidate ether ligands for Li⁺, Mg²⁺ and Na⁺ SILs with higher binding affinity and electrochemical stability than the reference compounds are identified. Lastly, their properties are related to the geometry of the coordination sphere.