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

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.

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts - often in collaboration with experimental groups - to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO2 capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.

Simulation of polymeric mixed ionic and electronic conductors with a combined classical and quantum mechanical model

In organic polymeric materials with mixed ionic and electronic conduction (OMIEC), the excess charge in doped polymers is very mobile and the dynamics of the polymer chain cannot be accurately described with a model including only fixed point charges. Ions and polymer are comparatively slower and a methodology to capture the correlated motions of excess charge and ions is currently unavailable. Considering a prototypical interface encountered in this type of materials, we constructed a scheme based on the combination of MD and QM/MM to evaluate the classical dynamics of polymer, water and ions, while allowing the excess charge of the polymer chains to rearrange following the external electrostatic potential. We find that the location of the excess charge varies substantially between chains. The excess charge changes across multiple timescales as a result of fast structural fluctuations and slow rearrangement of the polymeric chains. Our results indicate that such effects are likely important to describe the phenomenology of OMIEC, but additional features should be added to the model to enable the study of processes such as electrochemical doping.

Some Considerations about the Anodic Limit of Ionic Liquids Obtained by Means of DFT Calculations

Ionic liquids are good candidates as the main component of safe electrolytes for high-energy lithium-ion batteries. The identification of a reliable algorithm to estimate the electrochemical stability of ionic liquids can greatly speed up the discovery of suitable anions able to sustain high potentials. In this work, we critically assess the linear dependence of the anodic limit from the HOMO level of 27 anions, whose performances have been experimentally investigated in the previous literature. A limited r Pearson's value of ≈0.7 is found even with the most computationally demanding DFT functionals. A different model considering vertical transitions in a vacuum between the charged state and the neutral molecule is also exploited. In this case, the best-performing functional (M08-HX) provides a Mean Squared Error (MSE) of 1.61 V² on the 27 anions here considered. The ions which give the largest deviations are those with a large value of the solvation energy, and therefore, an empirical model that linearly combines the anodic limit calculated by vertical transitions in a vacuum and in a medium with a weight dependent on the solvation energy is proposed for the first time. This empirical method can decrease the MSE to 1.29 V² but still provides an r Pearson's value of ≈0.72.

Screening of the Role of the Chemical Structure in the Electrochemical Stability Window of Ionic Liquids: DFT Calculations Combined with Data Mining

Ionic liquids have attracted the attention of researchers as possible electrolytes for electrochemical energy storage devices. However, their properties, such as the electrochemical stability window (ESW), ionic conductivity, and diffusivity, are influenced both by the chemical structures of cations and anions and by their combinations. Most studies in the literature focus on the understanding of common ionic liquids, and little effort has been made to find ways to improve our atomistic understanding of those systems. The goal of this paper is to explore the structural characteristics of cations and anions that form ionic liquids that can expand the HOMO/LUMO gap, a property directly linked to the ESW of the electrolyte. For that, we design a framework for randomly generating new ions by combining their fragments. Within this framework, we generate about 10⁴ cations and 10⁴ anions and fully optimize their structures using density functional theory. Our calculations show that aromatic cations are less stable ionic liquids than aliphatic ones, an expected result if chemical rationale is used. More importantly, we can improve the gap by adding electron-donating and electron-withdrawing functional groups to the cations and anions, respectively. The increase can be about 2 V, depending on the case. This improvement is reflected in a wider ESW.

Giant and bidirectionally tunable thermopower in nonaqueous ionogels enabled by selective ion doping

Ionic thermoelectrics show great potential in thermal sensing owing to their ultrahigh thermopower, low cost, and ease in production. However, the lack of effective n-type ionic thermoelectric materials seriously hinders their applications. Here, we report giant and bidirectionally tunable thermopowers within an ultrawide range from −15 to +17 mV K⁻¹ in solid ionic liquid–based ionogels. Particularly, a record high negative thermopower of −15 mV K⁻¹ is achieved in the ternary ionogel, rendering it among the best n-type ionic thermoelectric materials under the same condition. A thermopower regulation strategy through ion doping to selectively induce ion aggregates to enhance ion-ion interactions is proposed. These selective ion interactions are found to be decisive in modulating the sign and magnitude of the thermopower in the ionogels. A prototype wearable device integrated with 12 p-n pairs is demonstrated with a total thermopower of 0.358 V K⁻¹, showing promise for ultrasensitive thermal detection.

Experimental measurement and prediction of ionic liquid ionisation energies

Ionic liquid (IL) valence electronic structure provides key descriptors for understanding and predicting IL properties. The ionisation energies of 60 ILs are measured and the most readily ionised valence state of each IL (the highest occupied molecular orbital, HOMO) is identified using a combination of X-ray photoelectron spectroscopy (XPS) and synchrotron resonant XPS. A structurally diverse range of cations and anions were studied. The cation gave rise to the HOMO for nine of the 60 ILs presented here, meaning it is energetically more favourable to remove an electron from the cation than the anion. The influence of the cation on the anion electronic structure (and vice versa) were established; the electrostatic effects are well understood and demonstrated to be consistently predictable. We used this knowledge to make predictions of both ionisation energy and HOMO identity for a further 516 ILs, providing a very valuable dataset for benchmarking electronic structure calculations and enabling the development of models linking experimental valence electronic structure descriptors to other IL properties, e.g. electrochemical stability. Furthermore, we provide design rules for the prediction of the electronic structure of ILs.

Comparison of the Performances of Different Computational Methods to Calculate the Electrochemical Stability of Selected Ionic Liquids

The electrochemical stability windows (ESW) of selected ionic liquids have been calculated by comparing different computational approaches previously suggested in the literature. The molecular systems under study are based on di-alkyl imidazolium and tetra-alkyl ammonium cations coupled with two different imide anions (namely, bis-fluorosulfonyl imide and bis-trifluoromethyl sulfonyl imide), for which an experimental investigation of the ESW is available. Thermodynamic oxidation and reduction potentials have here been estimated by different models based on calculations either on single ions or on ionic couples. Various Density Functional Theory (DFT) functionals (MP2, B3LYP, B3LYP including a polarizable medium and empirical dispersion forces) were exploited. Both vertical and adiabatic transitions between the starting states and the oxidized or reduced states were considered. The approach based on calculations on ionic couples is not able to reproduce the experimental data, whatever the used DFT functional. The best quantitative agreement is obtained by calculations on single ions when the MP2 functional in vacuum is considered and the transitions between differently charged states are vertical (purely electronic without the relaxation of the structure). The B3LYP functional underestimates the ESW. The inclusion of a polar medium excessively widens the ESW, while a large shrinkage of the ESW is obtained by adopting an adiabatic transition scheme instead of a vertical transition one.

Ionic Liquid Vapors in Vacuum: Possibility to Derive Anodic Stabilities from DFT and UPS

Ultraviolet photoelectron spectroscopy (UPS) investigations of several gas-phase ionic liquid (IL) ion pairs have been conducted. [EMIM][OTF], [PYR14][OTF], [EMIM][DCA], [PYR14][DCA], [PYR14][TCM], [PYR14][FSI], [PYR14][PF6], [S222][TFSI], [P4441][TFSI], and [EMMIM][TFSI] vapor UPS spectra are presented for the first time. The experimental low-binding-energy cutoff value (highest occupied molecular orbital, HOMO energy) of the ionic liquid ion pairs, which is of great interest, has been measured. Many studies use calculated gas-phase electronic properties to estimate the liquid-phase electrochemical stability. Hybrid density functional theory (DFT) calculations have been used to interpret the experimental data. The gas-phase photoelectron spectra in conjunction with the theoretical calculations are able to verify most HOMO energies and assign them to the cation or anion. The hybrid M06 functional is shown to offer a very good description of the ionic liquid electronic structure. In some cases, the excellent agreement between the UPS spectra and the M06 calculation validates the conformer found and constitutes as a first indirect experimental determination of ionic liquid ion-pair structure. Comparisons with recent theoretical studies are made, and implications for electrochemical applications are discussed. The new data provide a much-needed reference for future ab initio calculations and support the argument that modeling of IL cations and anions separately is incorrect.

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.

Engineering high-energy-density sodium battery anodes for improved cycling with superconcentrated ionic-liquid electrolytes

Non-uniform metal deposition and dendrite formation in high-density energy storage devices reduces the efficiency, safety and life of batteries with metal anodes. Superconcentrated ionic-liquid electrolytes (for example 1:1 ionic liquid:alkali ion) coupled with anode preconditioning at more negative potentials can completely mitigate these issues, and therefore revolutionize high-density energy storage devices. However, the mechanisms by which very high salt concentration and preconditioning potential enable uniform metal deposition and prevent dendrite formation at the metal anode during cycling are poorly understood, and therefore not optimized. Here, we use atomic force microscopy and molecular dynamics simulations to unravel the influence of these factors on the interface chemistry in a sodium electrolyte, demonstrating how a molten-salt-like structure at the electrode surface results in dendrite-free metal cycling at higher rates. Such a structure will support the formation of a more favourable solid electrolyte interphase, accepted as being a critical factor in stable battery cycling. This new understanding will enable engineering of efficient anode electrodes by tuning the interfacial nanostructure via salt concentration and high-voltage preconditioning.

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.

Boosting Rechargeable Batteries R&D by Multiscale Modeling: Myth or Reality?

This review addresses concepts, approaches, tools, and outcomes of multiscale modeling used to design and optimize the current and next generation rechargeable battery cells. Different kinds of multiscale models are discussed and demystified with a particular emphasis on methodological aspects. The outcome is compared both to results of other modeling strategies as well as to the vast pool of experimental data available. Finally, the main challenges remaining and future developments are discussed.

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.

Ion Dynamics in Ionic-Liquid-Based Li-Ion Electrolytes Investigated by Neutron Scattering and Dielectric Spectroscopy

A detailed understanding of the diffusion mechanisms of ions in pure and doped ionic liquids remains an important aspect in the design of new ionic-liquid electrolytes for energy storage. To gain more insight into the widely used imidazolium-based ionic liquids, the relationship between viscosity, ionic conductivity, diffusion coefficients, and reorientational dynamics in the ionic liquid 3-methyl-1-methylimidazolium bis(trifluoromethanesulfonyl)imide (DMIM-TFSI) with and without lithium bis(trifluoromethanesulfonyl)imide (Li-TFSI) was examined. The diffusion coefficients for the DMIM⁺ cation and the role of ion aggregates were investigated by using the quasielastic neutron scattering (QENS) and neutron spin echo techniques. Two diffusion mechanisms are observed for the DMIM⁺ cation with and without Li-TFSI, that is, translational and local. The data additionally suggest that Li⁺ ion transport along with ion aggregates, known as the vehicle mechanism, may play a significant role in the ion diffusion process. These dielectric-spectroscopy investigations in a broad temperature and frequency range reveal a typical α-β-relaxation scenario. The α relaxation mirrors the glassy freezing of the dipolar ions, and the β relaxation exhibits the signatures of a Johari-Goldstein relaxation. In contrast to the translational mode detected by neutron scattering, arising from the decoupled faster motion of the DMIM⁺ ions, the α relaxation is well coupled to the dc charge transport, that is, the average translational motion of all three ion species in the material. The local diffusion process detected by QENS is only weakly dependent on temperature and viscosity and can be ascribed to the typical fast dynamics of glass-forming liquids.

Valence electronic structure of [EMIM][BF4] ionic liquid: photoemission and DFT+D study

The ultraviolet photoelectron spectrum (UPS) of the [EMIM][BF₄] ionic liquid was recorded and analyzed. Together with the gas-phase UPS spectrum of the [EMIM][BF₄] vapor and ab initio calculation methods, detailed insight into the electronic structure of this simple ionic liquid is possible. The low binding energy tail in the UPS spectrum is about 7.4 eV, in agreement with previous estimations of the HOMO–LUMO gap of the [EMIM][BF₄] ion-pair. The bulk ab initio calculations are able to explain most of the features in the spectrum. However, DFT consistently lacks accuracy in the description of the top of the valence band. The dispersion corrected PBE calculation (PBE-D3) did offer very good agreement with the experimental structure, but the recently-developed vdW-DF functionals C09, optPBE, optB88 and CX were found to offer the best agreement in terms of the electronic structure.

The ultraviolet photoelectron spectrum (UPS) of the [EMIM][BF₄] ionic liquid was recorded and compared to previously measured vapor phase UPS spectrum.

Negative effective Li transference numbers in Li salt/ionic liquid mixtures: does Li drift in the “Wrong” direction?

Due to association with anions and an inverted drift direction in an electric field, Li⁺ cations have negative effective transference numbers.

The importance of transport property studies for battery electrolytes: revisiting the transport properties of lithium–N-methyl-N-propylpyrrolidinium bis(fluorosulfonyl)imide mixtures

All three ion–ion interactions contribute to transport properties in {Li[FSI]–[Pyr₁₃][FSI]} mixtures. Tracer diffusion coefficients of LI⁺ in [Pyr₁₃][FSI] are predicted.

A comparative study of room temperature ionic liquids and their organic solvent mixtures near charged electrodes

The structural properties of electrolytes consisting of solutions of ionic liquids in a polar solvent at charged electrode surfaces are investigated using classical atomistic simulations. The studied electrolytes consisted of tetraethylammonium tetrafluoroborate (NEt4-BF4), 1-ethyl-3-methylimidazolium tetrafluoroborate (c2mim-BF4) and 1-octyl-3-methylimidazolium tetrafluoroborate (c8mim-BF4) salts dissolved in acetonitrile solvent. We discuss the influence of electrolyte concentration, chemical structure of the ionic salt, temperature, conducting versus semiconducting nature of the electrode, electrode geometry and surface roughness on the electric double layer structure and capacitance and compare these properties with those obtained for pure room temperature ionic liquids. We show that electrolytes consisting of solutions of ions can behave quite differently from pure ionic liquid electrolytes.

Computational and Experimental Study of Li-Doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Ψ). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Ψ. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Ψ-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Ψ and suggest a source of impedance for Li+ to reach the electrode. Finally, we make predictions of the specific energy at ideal graphite utilizing the computed capacitance and previously derived electrochemical windows of the liquids.