Modeling Lithium Transport and Electrodeposition in Ionic-Liquid Based Electrolytes

Purely ionic electrolytes—wherein ionic liquids replace neutral solvents—have been proposed to improve lithium-ion-battery performance, on the basis that the unique microscopic characteristics of polarized ionic-liquid/electrode interfaces may improve the selectivity and kinetics of interfacial lithium-exchange reactions. Here we model a “three-ion” ionic-liquid electrolyte, composed of a traditional ionic liquid and a lithium salt with a common anion. Newman's concentrated-solution theory is extended to account for space charging and chemomechanical coupling. We simulate electrolytes in equilibrium and under steady currents. We find that the local conductivity and lithium transference number in the diffuse double layers near interfaces differ considerably from their bulk values. The mechanical coupling causes ion size to play a crucial role in the interface's electrical response. Interfacial kinetics and surface charge on the electrodes both affect the apparent transport properties of purely ionic electrolytes near interfaces. Larger ionic-liquid cations and anions may facilitate interfacial lithium-exchange kinetics.


INTRODUCTION
Improved lithium transport in separators can enhance the energy efficiency and power density of batteries. The use of ionic-liquid-based electrolytes in lithium-ion batteries has been shown to improve both current tolerance and cycle life (Lu et al., 2014;Piper et al., 2015). Numerous electrolytes based on ionic liquids have been characterized experimentally. Ionic-liquid-based electrolytes pose a challenge for modeling because they contain more than two ions, whose pairwise interactions strongly impact the electrolyte's overall performance. Although ionic-liquid solvents generally support solutions with higher ionic conductivity than conventional solvents, they tend to induce higher Haven ratios (Frömling et al., 2008).
The characterization of mass and charge transport in ionic liquids is impeded by the lack of satisfactory theoretical models. Conventional Nernst-Planck theory (dilute-solution theory) is based on an assumption that ions sit in a background of a neutral solvent species that is present in great excess, and therefore cannot be applied to fully ionic electrolytic solutions. For typical liquid electrolytes, characterization generally involves a series of measurements of individual transport properties (e.g., Hittorf transference number, Fickian lithium-salt diffusivity) whose very definitions are based on the notion that a neutral solvent is present. Even if this ambiguity is ignored, transport properties cannot be measured independently in concentrated solutions, because the thermodynamic properties and transport characteristics of different species are coupled. Moreover, as shown experimentally by Popovic et al. (2015), electrolytes at moderate-to-high salt concentrations contain multiple mobile species that are present with long lifetimes-including free ions, ion pairs, and higher-level structures-all of which contribute to ionic conductivity. Accurate characterization requires a group of orthogonal experiments and data analysis using coupled transport equations in a consistent way (Ma et al., 1995;Valoen and Reimers, 2005;Nyman et al., 2008;Hou and Monroe, 2020). A multispecies concentrated-solution theory for ionicliquid-based lithium electrolytes would help to link experimental characterization with continuum modeling. The foundational work of Pollard and Newman provides a useful starting point: they developed a model of molten-salt electrolytes based in concentrated-solution theory, applicable to domains where local electroneutrality holds (Pollard and Newman, 1979). We aim to extend that approach here.
In addition to their favorable bulk properties, ionic liquids have been suggested to provide a more stable electrolyte/electrode interface with regard to dendrite formation. Atomic force microscopy has been used to measure the force field near the ionic-liquid/electrode interface, showing a layered structure with enhanced mechanical strength (Gebbie et al., 2017;Hoffmann et al., 2018). Understanding how ionic liquids behave in close proximity to charged solid surfaces is crucial for electrolyte design.
Molecular dynamics has been widely deployed to study electrolyte structure near interfaces, especially ion arrangements in Stern layers . Selectivity for ion adsorption at electrode surfaces is highly variable: for example, anions can be chosen that either adsorb or desorb at a positively charged electrode . In spacecharge regions near interfaces, extremely high or low cationto-anion ratios may occur. Steric effects, rather than solubility considerations, determine the maximum ion concentrations in space-charge domains (McEldrew et al., 2018). Small variations of ion size have been observed to disrupt double-layer structure substantially (Bresme et al., 2012).
Continuum models have also successfully elucidated equilibrium structures of space-charge domains outside the Stern layer (i.e., diffuse double layers) near ionic-liquid/solid interfaces. These approaches tend to incorporate ion/ion interactions by specifying appropriate functionalities for the ions' electrochemical potentials (Kornyshev, 2007;Bazant et al., 2011;Hoffmann et al., 2018;de Souza et al., 2020). Short-range interactions and ion-size effects have been shown to rationalize the formation of layered structures near solid surfaces. Charge screening in the diffuse double layer can be accompanied by a significant pressure field near the interface (Monroe, 2020), and it remains unclear for fully ionic electrolytes how lithium cations would redistribute in response to this compressive force.
It should be emphasized that the practical use of fully ionic electrolytes in lithium-anode batteries is often accompanied by the formation of resistive solid-electrolyte interphase (SEI) layers (Xu, 2014). SEIs have been studied in detail for several fully ionic lithium-electrolyte systems (Xiong et al., 2014;Karimi et al., 2021). Highly resistive SEI layers can strongly impact both apparent reaction kinetics and interfacial capacitance. It is of great interest how fully ionic electrolytes respond as the SEI evolves, since SEI formation and growth during cycling is a significant mode of performance degradation.
A dynamical model of ionic-liquid-based electrolytes that considers mechanical effects is needed to study lithium transport in the diffuse double layer. Goyal and Monroe (2017) extended Newman's concentrated-solution theory by taking into account mechanics and electrodynamics. This paper deploys that framework to model a three-ion electrolyte system. As a case study we explore the baseline system comprising a binary solution of lithium tetrafluoroborate salt (LiBF 4 ) in the ionicliquid solvent 3-butyl-1-methyl imidaziolium tetrafluoroborate (bmimBF 4 ), which combine together to form the three-ion lithium electrolyte Li x bmim 1−x BF 4 . In later parts of the study we also explore the effects of changing both the anion and the ionic-liquid cation.
First we simulate the equilibrated diffuse double layer, to illustrate how pressure and ion partial molar volumes can affect property distributions near the interface. Then, we investigate how conductivity and transference number vary locally in the electrolyte in response to a steady current. Finally, we analyze how interfacial kinetics can affect lithium transport near the electrode surface. We also discuss how ionic-liquid-based electrolytes could be tuned to facilitate plating and stripping of lithium metal. We will emphasize how bulk ion transport is affected by changes in interfacial properties as practical systems age.

MODEL
For the baseline case, we consider an isothermal ternary ionic liquid system composed of Li + cations, bmim + cations, and BF 4 − anions, respectively indexed by subscripts k ∈ {Li, +, −}. Species molarities c k sum up to the total concentration c T . The mass density ρ = M k c k relates to species molarities through their molar masses M k . The number-average velocity of species k, v k , relates to its molar flux through N k = c k v k . The massaverage velocity of the fluid, v, is defined through the relation ρ v = M k N k . Introduction of Faraday's constant F and the equivalent species charges z k allow definitions of the local excess charge density ρ e = F z k c k and current density i = F z k N k .
Excess charge density produces an external electric field E governed by Poisson's equation, in which ǫ is the local dielectric permittivity, assumed to be constant. Each species follows a continuity equation, so that for all k ∈ {Li, +, −} A momentum balance relates the mass-average velocity to the hydrodynamic external pressure p, viscous stress τ , and electrostatic body forces ρ e E, in which both gravity and the diffusion stress discussed by Goyal and Monroe (2017) have been neglected. Hirschfelder et al. (1964) derived thermodynamic driving forces d k , which can be used to underpin Onsager-Stefan-Maxwell (OSM) equations that account for the energy dissipation through drag coefficients that resist relative species motion. In an isothermal system the OSM equations can be written as (Newman and Thomas-Alyea, 2004).
where R is the gas constant, T is the absolute temperature, and D kj is the Stefan-Maxwell diffusivity of species k through species j. Onsager reciprocity guarantees that D kj = D jk . Because the Gibbs-Duhem relation requires that k d k = 0, only two independent OSM equations suffice to model a ternary diffusion system. Inspection of the dissipation function for a liquid electrolyte suggests that the viscous stress is given by in which effects due to bulk viscosity have been discarded. We further take the electrolyte to be a Newtonian fluid, with constant viscosity η. For species k, composition variation, pressure gradients, and the electric field can all contribute to the electrochemicalpotential gradient, through More complex electrochemical-potential functionalities can be used to account for additional details about ion behavior, such as solvation energetics. Several electrochemical-potential functions have been proposed to model short-range ion/ion interactions (Bazant et al., 2011;Hoffmann et al., 2018;de Souza et al., 2020).
Here we restrict our discussion to ideal mixtures, as specified by the ln(c k /c T ) term in Equation (6). A thermodynamic analysis of the equation of state for volume also shows that the set of molar species concentrations {c k } is constrained by the set of species partial molar volumes {V k }, through Assuming a constant composition-independent bulk modulus K, each partial molar volume takes the form where V θ k is the partial molar volume at reference pressure p θ . Experimental values of V θ k are expected to reflect information about not only central species k, but also secondary and tertiary structures involving k, such as its solvation shell or coordination environment.
It is useful to express composition through ion k's particle fraction y k = c k /c T , which is convenient because y − = 1 − y + − y Li . Each species molarity c k can be written as a function of {y + , y Li , p} through which derives from Equations (7) and (8). Thus, all thermodynamic properties of each species are uniquely determined by the pressure (p), external electric field ( E), and the particle fractions of lithium and bmim cations (y Li and y + ).
The system of Equations (1)-(9) is closed by stating constitutive laws for the electrical properties ǫ and z k , the thermodynamic properties M k , V k , and K, and the transport properties D kj and η. Following Monroe and Delcourt's process (Monroe and Delacourt, 2013), the OSM transport laws from Equation (4) can be transformed into flux-explicit forms that include a generalized form of Ohm's law, as well as flux-explicit transport laws for bmim cations, and lithium ions, (12) These transport equations have been written in forms that use v − as a reference velocity for convection. Three transport properties-conductivity κ − , lithium transference number relative to the hexafluorophosphate velocity t − + and cross diffusivity for cation exchange D − × -link the thermodynamic driving forces to species fluxes and the ionic current density. In all cases the superscript − is included to emphasize that the values of these properties depend on the choice of reference velocity.
The three macroscopic transport properties that appear in Equations (10)-(12) map bijectively into the three Stefan-Maxwell diffusivities D +− , D Li+ and D Li− : Note that the first-order pressure dependence of these properties is determined entirely by the total molarity c T = 1/ y k V k p . The ions' molar masses (M + = 139.2 g/mol, M − = 86.8 g/mol, and M Li = 6.9 g/mol) are well known. Partial molar volumes at ambient pressure have also been reported by Slattery et al. (2007) for bmim + and BF − 4 (V θ + = 118.034 mL/mol and V θ − = 43.962 mL/mol). For lithium cations, the partial molar volume is assumed to be the same as the molar volume of lithium in its pure metal state, V θ Li = 13.109 mL/mol. The bulk modulus of the electrolyte is taken to be that of pure liquid bmimBF 4 , K = 2.8831 GPa (Sunkara et al., 2015); and the relative permittivity ǫ is 11.7 at 303 K and 0.1 MPa (Wakai et al., 2005). At 303 K, the ionic conductivity of pure bmimBF 4 is κ 0 = 0.295 S/m (Yu et al., 2009), which can be used to estimate the Stefan-Maxwell diffusivity D +− = 5.5017 × 10 −12 m 2 /s according to the relation κ 0 = F 2 z 2 + c T D +− /(RT) obtained by assuming y Li = 0 in Equation (13). D Li+ and D Li− have not been reported systematically, so we assume D Li− = D +− = D Li+ for simplicity, unless stated otherwise.
The present analysis focuses on transport phenomena within diffuse-double-layer and bulk-electrolyte domains. Regions between the electrode surfaces and the outer Helmholtz planes-including the SEI, Stern layer, etc.-are treated phenomenologically, in terms of the interfacial characteristic frequency f int , a property whose value can be deduced from the interfacial impedance.

Equilibrium Double Layers Near Blocking Electrodes
We first investigate the double layer structure in a Li 0.5 bmim 0.5 BF 4 system between two planar blocking electrodes biased by a 20 mV voltage difference. A one-dimensional system of thickness L = 4 nm is simulated. Although this electrode spacing is unrealistically small, it helps to make the distributions of properties within the interfacial double layers more readily visible in plots. Generally the structures of double layers computed with this model do not change when electrodes are moved further apart.
As Figure 1A shows, ionic liquids between blocking electrodes form an electrolytic capacitor, in which two double layers form near the electrode surfaces and the electric field is screened from the bulk liquid. Figure 1B shows that the fractions of both types of cation rise (fall) near the negative (positive) electrode. The external field imposes forces on the excess charge in the double layers, and a positive gauge pressure around 3.5 MPa builds up near both electrodes, as shown in Figure 1C.
At equilibrium, Poisson's Equation (1) and the momentum balance, Equation (3), combine to give a relationship between the electric field and the pressure, derived using the observation that E bulk ≈ 0.
When electrodes are blocking, the field is equal at the negative and positive electrodes' surfaces. Equation 16 requires the peak pressure to be identical at both edges of the electrolyte, as seen in Figure 1C. It is worth noting that the pressure distribution is not symmetric-its curvature differs at the electrolyte's boundaries. Due to the disparity in ion sizes, the positive and negative outer Helmholtz planes achieve different total ion concentrations, as shown in Figure 1B, causing the screening lengths at the positive and negative boundaries to differ. This effect leads to the pressure-distribution asymmetry in Figure 1C, as well as the difference between interfacial potential drops in Figure 1A.
Since Li + has both a smaller size and a lower mass than bmim + , the substitution of Li + for bmim + alters the energetics associated with double-layer formation. The lithium-ion fraction is apparently enhanced in both double layers, a chemomechanical effect illustrated in Figure 1D through the distribution of the lithium excess factor ζ = y Li /ȳ Li y + /ȳ + − 1, whereȳ Li andȳ + are the global-average particle fractions of Li + and bmim + . A positive (negative) ζ means that the ratio of lithium ions to bmim cations is higher (lower) that its average value. Even at the relatively low potential bias of 20 mV, there are about 60% more Li + ions than bmim + ions near the positive electrode and about 30% more near the negative electrode.
In the absence of current, Equations (4) and (6) together give equilibrium conditions for the distributions of bmim + and Li + , Subtracting Equation (17) from (18) produces a relationship between the gradients of lithium excess factor and pressure, In all cases studied here, the local mass density ρ satisfies so double-layer regions will have a higher lithium fraction than the global average. It is straightforward to simulate the equilibrium states of different mixtures of lithium salts Li + Y − and ionic liquids X + Y − that differ from the baseline case shown in Figure 1. In these simulations we continue to fix the potential bias at 20 mV. The ionic-liquid cations are selected from C n mim + with n = 2, 3, 4, 5, 6, and 8 (note that bmim + corresponds to n = 4), and the anions are varied between BF 4 − and Tf 2 N − [i.e., bis(trifluoromethane)sulfonimide, a.k.a. TFSI], for which partial molar volumes are available from the extensive report by Slattery et al. (2007). The partial molar volume of Tf 2 N − is 139.71 mL mol −1 , around three times that of BF 4 − and similar to bmim + . The global-average particle fraction of Li + was varied from 0.05 to 0.25, corresponding to a molarity range of approximately 0.6-5 M for the lithium salt. Figure 2A shows that given a particular choice of X + and Y − , a higher lithium fraction produces a higher surface pressure. Smaller partial molar volumes for X + and Y − also produce higher surface pressure given a fixed lithium fraction. This behavior owes to the fact that smaller ion sizes and larger lithium fractions give higher ion concentrations according to Equation (9), thereby producing thinner double layers. At a given potential bias, a thinner double layer implies a larger electric field, and a consequent higher surface pressure according to Equation (16). Figures 2B,C show the lithium excess factor ζ at the surfaces of the positive and negative electrodes, respectively. The lithium excess is positive at both interfaces, ranging from 25 to 80% at the positive electrode and 15-35% at the negative. For a given X + , a higher lithium fraction produces a higher lithium excess factor ζ near both electrodes. A larger anion can enhance ζ near the positive electrode up to 20%, whereas anion size has negligible effect on ζ near the negative electrode due to the low anion fraction there. In systems with the same lithium fraction and same choice of anion, a larger ionic-liquid cation produces a higher lithium excess near the negative electrode but a lower ζ near the positive. Changing cation size (from C 2 mim + to C 8 mim + ) can change ζ by 5%, but the effect goes in opposed directions at the positive and negative electrode.

Lithium Transport Between Reactive Electrodes
Interfacial impedance correlates the electric field at the boundary (i.e., the outer Helmholtz plane) with the faradaic current density. This property gives a phenomenological description of all processes happening between the electrode and the outer Helmholtz plane. In an earlier paper (Li and Monroe, 2019), the authors modeled faradaic interfaces with the relationship in which i is the applied current density, and R int and C int are the effective kinetic resistance and interfacial capacitance, respectively. The combined effects of resistance and capacitance are summarized by the single characteristic frequency f int = 1/(2πR int C int ), which can be measured with electrochemical impedance spectroscopy and controls how interfaces respond to a potential bias. Interfaces with small f int values behave more like electrolytic capacitors; large f int interfaces behave like resistors. Salminen et al. (2006) performed an interfacial impedance measurement from 0.1 to 65535 Hz for Li|bmimBF 4 +1 or 2 mol/kg LiBF 4 |Li sandwich cells, providing a complete semicircle on the Nyquist plot. The interfacial frequency is given by the peak frequency of the interfacial semicircle (Li and Monroe, 2019), so a practical interfacial frequency f int is expected to lie between 1 Hz and 10 kHz in a three-ion lithium electrolyte. In this range, the interface will be somewhat capacitive due to the likely formation of a highly resistive SEI, leading to the distinguishable presence of a double layer in which the potential drop is steep. It is also worth noting that Stern-layer phenomena merely change the relationship between the electric field at the outer Helmholtz plane and the faradaic current, and are therefore encapsulated by f int . Figure 3 summarizes simulations of a 10-µm-thick Li 0.5 bmim 0.5 BF 4 slab sandwiched between two reactive electrodes with f int = 1 kHz (2πR int C int = 10 −3 s), subjected to a potential bias of 20 mV. Charge transfer at the electrode/electrolyte interface is assumed to owe only to Faradaic reactions involving Li + . The cation/cation diffusivity D Li+ controls the drag interactions between lithium ions and ionicliquid cations, and is studied across values of D +− /D Li+ = −0.5, 0, and 1, all of which are thermodynamically consistent (Monroe and Delacourt, 2013). When D +− /D Li+ = 0, cations place no drag on each other; Li + and bmim + flow independently. The flow of bmim + slows the flow of Li + when D Li+ > 0, and speeds lithium flux when D Li+ < 0.
As shown in Figure 3A, the double layer does not entirely screen the interfacial electric field near a reactive electrode. Some unscreened field is associated with Ohmic conduction across the cell. For the f int value chosen, the potential drop across the double layers is fairly large in comparison to the drop across the bulk. Figure 3B shows that the concentrations of both cations (Li + and bmim + ) are polarized with opposite slopes in the bulk. The negative concentration gradient of the reactive species, Li + , coincides with the direction of the ionic current, while that of bmim + is opposed; anions distribute relatively uniformly, maintaining bulk electroneutrality. As the ratio D +− /D Li+ gets more positive, the trends remain similar and the magnitude of the polarization effect grows larger. The local values of conductivity and transference number that result from ion polarization can be calculated from Equations (13) to (14). Figure 3C shows that the electrolyte's conductivity very near the cathode can be significantly larger than the bulk, while that near the anode is smaller. The local conductivity depends strongly on compositionespecially on Li + fraction-whereas the transference number is primarily determined by the lithium excess factor. Figure 3D shows that lithium transference decreases across the bulk due to the opposed polarizations of the Li + and bmim + particle-fraction distributions. In the double layer, the lithium transference number is enhanced due to the pressure effect discussed earlier.
Faradaic currents induce bulk cation polarization and reduce the lithium concentration near the negative electrode, which could encourage side reactions. Three phenomena specific to double layers may mitigate the adverse consequences of bulk polarization at the cathode: conductivity is enhanced, Li + concentration is increased by double-layer polarization, and lithium transference is enhanced by the pressure effect. At the anode surface-somewhat counterintuitively-the lithium concentration can be substantially depleted if double-layer polarization is too strong. In a purely ionic electrolyte, the pressure effect that accompanies charge screening enhances both lithium concentration and transference number at both electrodes, potentially improving interfacial kinetics for any reactions involving lithium. Figure 4 shows how an electrolyte slab with baseline properties responds to a potential bias of 20 mV when the character of the electrode/electrolyte interfaces changes. Three f int values are studied: 20 kHz, 1 kHz, and 0 Hz. A small interfacial frequency corresponds to a situation with sluggish interfacial kinetics and a very capacitive interface. In this case, most of the potential drop occurs in the double layers, as shown in Figure 4A. Smaller bulk potential drop is associated with less bulk polarization and smaller currents (cf. Figure 4B) and more significant surface pressure ( Figure 4C). Local conductivity and local lithium transference number are both substantially affected by the interfacial frequency, as can be seen from the plots on Figures 4D,E. In the case of constant potential bias, interfacial kinetics impacts the cell response in a coupled way, changing both bulk polarization and double-layer structure.

Figure 5
shows simulations of four 10-µm-thick systems (Li 0.5 C 4 mim 0.5 BF 4 , Li 0.5 C 8 mim 0.5 BF 4 , Li 0.5 C 4 mim 0.5 Tf 2 N, and Li 0.5 C 8 mim 0.5 Tf 2 N) with different interfacial frequencies in the 6.3 Hz to 1 kHz range. A practical current density of 0.1 mA/cm 2 was used for the simulations. These data illustrate how transport near the electrodes is expected to vary in response to a change in interfacial kinetics. Practically, battery degradation is usually associated with increasing interfacial resistance due to SEI growth, morphology change, passivation, etc., which would correspond in these simulations to a lower interfacial frequency f int . Thus, as practical systems age, their responses would be expected to move from high f int to low f int along the curves in Figure 5.
The surface pressure induced by applied current can be found by combining Equations (16) and (21). After assuming E bulk ≪ E surf , one finds that For the same applied current, a smaller f int produces a higher surface pressure, and should therefore raise the lithium fraction near the electrodes more. In the high interfacial frequency limit f int → ∞, the properties at the electrode surface match their properties in the bulk; double-layer effects become insignificant. Conversely, the effects at low interfacial frequencies arise primarily from double-layer polarization. It is worth noting that the concentrations of both cations (Li + and C n mim + ) can be dramatically reduced at the positive electrode as a consequence of space charging as f int → 0. Figure 5A compares the local conductivity near the interface as interfacial frequency varies. As interfaces get more capacitive (smaller f int ), higher double-layer polarization can enhance the conductivity adjacent to the cathode by nearly half an order of magnitude, accompanied by a raised lithium transference number in the range 0.6-0.9. The lithium depletion in the double layer near the anode reduces the conductivity by a factor of ∼ 5, but the pressure effect in the double layer still significantly enhances the lithium fraction and the lithium transference number there, as can be seen from Figure 5B. Thus, the use of ionic liquids could improve lithium transport-especially via the surface transference-number enhancement-at both anode and cathode surfaces in systems where interfacial kinetics slows down over time (as it might due to interfacial degradation). Comparing cells run at the same current density, the interfacial performance exhibited by ionic liquids with larger cations and anions is always expected to be better than ionic liquids with smaller ions, as shown in Figure 5.

CONCLUSION
Concentrated-solution theory was exploited to develop a transport model for a three-ion, ionic-liquid-solvated lithium electrolyte, i.e., a mixture of a lithium salt and ionic liquid that share an anion, such as the ternary Li x bmim 1−x BF 4 system. The theory included three pairwise ion/ion diffusional drag interactions, as well as allowing for electrical non-neutrality and accounting for mechanical effects due to space charge. Simulations explored the steady-state response of this ternary ionic system to a constant potential bias or a constant applied current. The double-layer structure near the electrode surfaces was investigated in detail. A significant surface pressure (∼3.5 MPa for a 20 mV potential bias) can build up in the double layer as a consequence of the Lorentz force that acts on excess charge there. This surface pressure was observed to change the relative amounts of Li + and bmim + near interfaces significantly due to the disparities in their molar masses and sizes, greatly enhancing the surface Li + concentration. In a Li 0.5 bmim 0.5 BF 4 system under a 20 mV potential bias, 60% more Li + than bmim + was present near the positive electrode, and 30% more was present near the negative electrode. A comparison of various electrolyte formulations showed that the use of large anions can enhance surface lithium concentrations even more. Large cations have more of an effect at the cathode than at the anode. The use of ionic-liquid solvents could greatly facilitate lithium plating and stripping processes via this pressure effect.
Local conductivity and lithium transference number were also studied in a Li 0.5 bmim 0.5 BF 4 system. The interaction between Li + and bmim + was found to affect bulk cation polarization significantly. The response of electrolytes to currents is determined by the interfacial impedance. If an interface becomes more capacitive as the cell ages, as might be expected if electrode passivation increases over time, currents will induce a stronger double-layer polarization and a more significant pressure effect. Consequently, the local conductivity and lithium transference number near the negative electrode will increase to assist plating processes, while the pressure effect can raise the lithium transference number on the stripping side to compensate for conductivity reduction there. In a comparison of different purely ionic lithium electrolytes at constant currents, a combination of larger ionic-liquid cations and anions was always found to produce higher local lithium transference numbers at both interfaces.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article; further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
GL and CWM designed the research. GL implemented simulations and drafted the manuscript. CWM revised the manuscript and supervised. Both authors contributed to the article and approved the submitted version.