Modeling the voltage loss mechanisms in lithium-sulfur cells : the importance of electrolyte resistance and precipitation kinetics

Lithium-sulfur (Li-S) cells have the potential to deliver practical energy densities exceeding 400 Wh/kg at cell level in the next few years. However, a lack of understanding of the complex electrochemical, transport, and phase-change phenomena in Li-S cells is arguably holding back development. Acquiring this knowledge requires experimental characterization in tandem with mechanistic modeling. However, existing Li-S models currently contradict some key features of experimental findings, particularly the evolution of cell resistance during discharge. In this study, we demonstrate that, by introducing a concentrationdependent electrolyte conductivity, the correct trends in voltage drop due to electrolyte resistance and activation overpotentials are retrieved. In addition, we reveal the existence of an often overlooked potential drop mechanism in the low voltage-plateau which originates from the limited rate of Li2S precipitation. Understanding these voltage drop mechanisms is important for correctly interpreting experimental results and for parametrizing advanced mechanistic models. While these effects are exemplified here via a lumped mechanistic model, any existing Li-S models should benefit from including them.


Introduction
The lithium-sulfur (Li-S) cell could provide the next step-change in battery technology with a promising practical energy density of 500-600 Wh/kg.However, broader uptake of Li-S technology is currently hampered by relatively fast capacity fade, incomplete sulfur utilization, and low Columbic efficiency due to the polysulfide shuttle, among other issues [1,2,3,4,5].Improvements in battery design and materials research require a good understanding of the complicated electrochemical processes in the Li-S system, which involve multistep reactions, multiscale transport phenomena, and morphology variations during battery operation.Modelling in tandem with experimental testing has been shown to accelerate research and development in the field.[6,7].
Mikhaylik et al [8] developed the first Li-S model for studying polysulfide shuttle.This lumped model computes reduction potentials in the high and low voltage-plateaus based on the Nernst equation, but it neglects activation overpotentials, electrolyte resistance, and dissolution/precipitation reactions, all which influence the operating voltage.A more detailed one-dimensional (1D) mechanistic model was later described by Kumaresan et al [9] based on the Nernst-Plank equations for dilute solutions.Kumaresan's model considers multicomponent transport in the porous cathode and separator, charge-transfer kinetics, dissolution/precipitation reactions, and changes in porosity and electrochemically active surface area as a result of dissolution/precipitation. While the model qualitatively reproduces some essential features of a typical Li-S discharge profile -such as the two characteristic voltage plateaus and the voltage dip in-between, it also requires the input of a large number of physical and phenomenological parameters that are not easily obtained experimentally.In view of the complexity of the Kumaresan model, Ghaznavi and Chen [10,11,12] performed a sensitivity analysis of this model.They concluded that, in order for this model to predict charging, further development in modelling precipitation/dissolution reactions is required.Bessler et al [13,14] developed a 1D Li-S model based on a similar framework, which additionally considers the electrochemical double-layer, polysulfide shuttle, and irreversible precipitation at the anode as a capacity fading mechanism.In addition to charge/discharge curves, electrochemical impedance spectra at different depths of discharge (DoD) were predicted.
However, while the mechanistic models by Kumaresan et al [9] and Bessler et al [13,14] can reproduce the general discharge profile of a Li-S cell, they fail to capture more intricate cell behavior such as the changing electrolyte resistance (R s ) during discharge.The variation of R s with DoD is a characteristic feature of Li-S cells that has been observed in various electrolyte chemistries by electrochemical impedance spectroscopy (EIS) [15,16,17,18].It is generally agreed that the high-frequency intercept of the EIS curve is dominated by the electrolyte resistance of the measured Li-S cell.As shown in Fig. 1a, during discharge R s increases in the high plateau and reaches a maximum at the transition between the two voltage plateaus; it then decreases throughout the low-plateau.Furthermore, EIS measurements suggest that R s accounts for most of the voltage-drop in the high plateau, as the low-frequency resistanceswhich are frequently associated with charge-transfer -only become significant in the low-plateau [15,18,19].In the mechanistic models, however, the predicted voltage drop due to electrolyte resistance is two to three orders of magnitude smaller than the predicted voltage-drop due to activation overpotentials.Such model prediction, directly contradicting EIS data, is visible in the simulated EIS curves by Bessler et al [13].The disagreement is also seen for predictions produced from Kumaresan's model [9], shown in Fig. 1b, where the simulated R s exhibits a very different evolution with DoD compared to measurements.
The R s profile during discharge is typically explained with the variation of electrolyte concentration (or viscosity) caused by the dissolution and subsequent precipitation of polysulfide species.Indeed, it was experimentally observed that the electrolyte conductivity of Li-S cells strongly depends on the concentration of both lithium salt [20] and lithium polysulfides [21].At high salt or polysulfide concentrations (typically > 1 mol/L), increased ionic interactions reduce the electrolyte conductivity.Existing Li-S models, however, rely on dilute solution theory in which ionic conductivities are independent of ionic concentrations.
We propose that introducing a concentration dependence of the electrolyte conductivity is necessary to retrieve the experimentally documented trends in the voltage-drop in Li-S cells during discharge.This feature is included in a lumped mechanistic model that describes electrochemical and precipitation reactions, electrode charge-transfer kinetics, as well as morphology variations due to precipitation in a Li-S cell.Compared to the more detailed Kumaresan model, the lumped model does not consider mass-transport and charge-localization effects, and therefore it cannot predict transport limitations.We note, however, that the sensitivity analysis by Ghaznavi and Chen [12] indicates that mass transport does not have a significant impact on predictions of the Kumaresan model unless ionic diffusion coefficients are reduced by more than an order of magnitude.Consequently, the discharge curves produced by the lumped model closely resemble those obtained with the Kumaresan model.The advantage of the lumped approach is that the model requires fewer fitting parameters and reduced computational resources.

Model formulation
The model considers six electrochemical reactions and one precipitation reaction during discharge: 0.5S 8 + e -→ 0.5S 2-

8
(2) The model considers six electrochemical reactions and one precipitation reaction during discharge: Eq. 1-6 describes the typical reactions in a catholyte-type cell in which sulfur is initially dissolved in the electrolyte prior to discharge.However, the model can be easily modified to include sulfur dissolution as well as additional precipitation reactions.The cell voltage can be written as the contribution of three terms: where E j and η j are the reduction potential and the activation overpotential for a cathodic reaction j respectively, E 1 and η 1 are the anodic reduction potential and overpotential, and IR s is the potential drop due to electrolyte resistance.We assume that the Li + dissolution kinetics are sufficiently fast such that the anode overpotential η 1 is negligible [13].All three components are dependent on the species concentrations, C i , which vary with time due to electrochemical reactions: Here, ε is the spatially-averaged cell porosity, i j is the current density due to electrochemical reaction j, s i,j is the stoichiometric coefficient of species i in reaction j, n j is the number of electrons transferred in reaction j, F is the Faraday constant, and a v is the specific surface area for electrochemical reactions.
As this is a 0D model, there is no spatial variation of concentration due to mass transport.The concentration variation of the two species participating in the precipitation reaction in Eq. 7 is given by: where r p represents the precipitation rate of Li 2 S, C Li + ,0 is the initial Li + concentration, and the Li + concentration derives from the charge conservation.We employ the expression proposed by Kumaresan et al [9] to describe the rate of precipitation: Here, v Li2S is the volume fraction of Li 2 S with respect to cell volume, k p is the precipitation rate constant, K sp is the solubility product, and V Li2S is the molar volume of Li 2 S. According to Eq. 11, precipitation occurs when the concentration product, C 2 Li + C S 2− , exceeds the solubility product, K sp .Furthermore, the rate of precipitation is taken to be proportional to the amount of precipitated Li 2 S that provides the necessary surfaces for solid-phase nucleation and growth.
The current densities are related to overpotentials through the Butler-Volmer equation.If the anodic and cathodic transfer coefficients are assumed to both equal 0.5, the current-overpotential relation can be written as: in which i 0 j is the exchange current density for reaction j, R is the ideal gas constant and T is the temperature.Furthermore, the volumetric current densities due to electrochemical reactions are constrained by the applied current I via: where A is the apparent geometric area of the cell and l is the cell thickness.
The reduction potential for a reaction j is given by the Nernst equation where E 0 j stands for the standard reduction potential for reaction j at the reference concentration of 1 mol L −1 .
The total electrolyte resistance across the cell can be written as where σ is the electrolyte conductivity.During discharge, electrochemical and precipitation reactions significantly alter the polysulfide concentrations, which in turn influence the electrolyte conductivity.While the exact relation between ionic concentration and electrolyte conductivity has not been established for Li-S cells, it is reasonable to assume that upon increasing Li + concentration, the ionic conductivity of the electrolyte first increases and then decreases.This effect is observed in the electrolytes in both lithium-ion and Li-S cells [20,21,22].We further assume cells are usually built with the electrolyte salt concentration close to the optimal value for the maximum electrolyte conductivity (private communication with Oxis Energy LTD).With these assumptions, we propose a linear phenomenological function for σ at high Li + concentrations: In this expression C Li + is taken to represent the total ionic concentration of the electrolyte as it equals the sum of anion concentrations.In writing Eq. 15 we have assumed that the electrolyte conductivity is only a function of the total anion concentration (represented by C Li + ) instead of the concentrations of each individual ionic species.This limitation is due to the unknown transport properties of the various polysulfide dianions that may exist during cell discharge.
Measuring the properties of individual polysulfide species is challenging due to the complex chemical equilibria among polysulfide species [23].As shown in Fig. 2a, the electrolyte conductivity reaches the maximum value of σ 0 at Li + concentration C Li + ,0 , then decreases with slope b.Finally, we employ the phenomenological expression used by Kumaresan et al [9] to describe the change in specific electrochemical surface area with cathode porosity: where ε 0 is the initial porosity and ξ is a fitting parameter.The model contains algebraic and ordinary differential equations that can be solved numerically.The model parameters and their assumed values are listed in Table 1.Due to the lack of reported data on standard reduction potentials, exchange current densities, and precipitation kinetics, the parameters in the lumped model -as with the parameters in other Li-S models in the literature -are obtained from calibrating the model with measured discharge curves.Compared to the Kumaresan model [9], the lumped model does not contain the numerous parameters associated with ionic transport, but requires two additional parameters to describe the concentration-conductivity relation.

Results and discussion
The simulated discharge curves and electrolyte resistances at 0.15C (0.34A) and 0.03C are illustrated in Fig. 2c-d.Similar to the Kumaresan model [9], the lumped mechanistic model is able to qualitatively reproduce the essential features of the discharge profile of Li-S batteries, e.g. a sloping high-plateau voltage, a flat low-plateau voltage, and a voltage dip in-between.However, the lumped model also reproduces the correct trend and magnitude of the electrolyte resistance during discharge due to the introduction of concentration-dependent electrolyte conductivity.The shape of the resistance profile follows closely the evolution of Li + concentration during discharge as shown in Fig. 2b.Initially, as an increasing amount of Li + dissolves into the solution to form Li 2 S x , the electrolyte conductivity reduces according to the proposed conductivityconcentration relation Eq. 16.In the low plateau, as Li + precipitates out of the solution as Li 2 S, the electrolyte conductivity increases.The peak of R s therefore corresponds to the onset of Li 2 S precipitation.It follows that the voltage dip between the voltage plateaus is not only due to the super-saturation of S 2-as described by Kumaresan et al [9], but also a consequence of electrolyte resistance peaking at the transition between the two voltage plateaus.
According to Fig. 2d, the peak in electrolyte resistance increases with discharge current.At higher currents, the electrochemical production rates of Li + and S 2− are faster, whereas their chemical precipitation rate remains the same at the onset of Li 2 S precipitation.Consequently, the concentrations of Li + and S 2− are higher, which gives rise to larger electrolyte resistance for a higher discharge current.It is clear from our analysis that the exact shape of the R s profile is strongly influenced by the precipitation rate of Li 2 S. In the present model, however, the precipitation rate follows Kumaresan's phenomenological expressions Eq. 11, which are known to not reproduce re-dissolution of Li 2 S upon charging [11,12].A more accurate precipitation/dissolution model would require detailed solid phase nucleation and growth mechanisms as well as experimentally determined solubility products and precipitation rates, values currently not well established in Li-S literature.
The simulated activation overpotentials (defined in Eq. 12) for the highplateau process reaction Eq. 2 and the low-plateau process reaction Eq. 6 are shown in Fig. 3.It is clear that the smaller high-plateau overpotential η 2 remains relatively constant in the high plateau whereas the larger low-plateau overpotential η 6 increases with increasing DoD.It is found that the increase in low-plateau overpotential is primarily attributed to the reduction in available electrochemical surface area due to Li 2 S precipitation, which is also shown in Fig. 3.As the insulating Li 2 S gradually covers a larger portion of the conductive cathode surface, larger activation overpotentials are required to drive the increasing electrochemical current densities.We note that the simulated trend in activation overpotentials resembles that of the low-frequency resistances measured by EIS, which also remain small in the high plateau but rise quickly in the low plateau [15,18,19].Since the activation overpotentials qualitatively reflect the resistance due to charge-transfer, the model agrees with the hypothesis that the low-frequency resistances in EIS measurements arise from charge-transfer processes [15,18,19].
In addition to the potential drop associated with the electrolyte resistance and activation overpotentials, a potential shift also occurs due to Li 2 S precipitation that alters the reduction potentials in the Nernst equation (Eq.14).In the absence of precipitation, the S 2− concentration increases continuously during discharge in the low plateau, which causes the reduction potentials to drop gradually as dictated by the Nernst equation.This scenario is depicted in Fig. 4 for the low-plateau reduction potential E 6 without precipitation, which is similar to the low-plateau reduction potential calculated by Mikhaylik et al [8].In the presence of precipitation, the electrochemical production rate of S 2− reaches a dynamic equilibrium with its removal rate due to precipitation, thereby holding the S 2− concentration and the reduction potential E 6 at relatively a constant values.As illustrated in Fig. 4, the precipitation effectively increases the reduction potential in the low-plateau and gives rise to a flat low-plateau potential.Furthermore, since the precipitation rate is finite, the equilibrium concentration of S 2− is higher at higher currents and the reduction potential is correspondingly lower.The limited rate of precipitation effectively manifests as a 'precipitation resistance' that leads to a drop in E 6 at higher currents, as demonstrated in Fig. 4. Due to the quasi steady-state conditions of the testing procedure, EIS measurements do not explicitly reflect reduction potential drop.However, the shift in equilibrium potential has been reported experimentally by the galvanostatic intermittent titration technique (GITT), for which a long voltage relaxation (∼ 40 mV over 20 h) was observed in the low plateau after a current pulse had been removed [24].This slow potential equilibration was attributed to the slow precipitation/dissolution kinetics in Li-S cells.Since the magnitude of the reduction potential shift is sensitive to the precipitation rate, the GITT technique might be useful for estimating the precipitation rate constants and solubility products needed for more detailed Li-S models.
In the graphical abstract, the three voltage-drop mechanisms considered in the model are shown together for the representative low-plateau reaction Eq. 6.
We note that all three mechanisms are related to the precipitation of Li 2 S: the precipitation affects the electrolyte resistance and reduction potentials through the change in polysulfide concentrations, and influences activation overpoten-tials via the change in available electrochemical surface area.For practical, high energy-density Li-S cells, the electrolyte resistance is particularly important since these cells generally contain less electrolyte and therefore higher ionic concentrations [25].It is therefore important for high energy-density Li-S cells to employ solvents with high ionic conductivities as well as a thin separator to minimize the Ohmic voltage loss.The voltage-drop due to activation overpotentials and the apparent 'precipitation resistance' are more difficult to quantify since they are dependent on the poorly understood reaction mechanisms and precipitation kinetics in Li-S cells.While solvents with low Li 2 S solubility could facilitate precipitation thereby potentially reducing both the electrolyte resistance and the 'precipitation resistance', the increased amount of insulating precipitates induce larger activation overpotentials.Large activation overpotentials could be mitigated through the use of nano-structured cathodes with high conductive surface area.
In the present work only a qualitative comparison between the model prediction and experimental data of cell discharge and electrolyte resistance curves can be made.This limitation is due to: (i) values for many physical parameters, especially those related to Li 2 S precipitation and electrolyte conductivity, have not been obtained experimentally, nor are they established in the literature; (ii) the present model neglects the polysulfide shuttle as well as transport limitations that could occur at high currents, therefore it does not sufficiently capture the variation of discharge capacity with discharge current.However, improvements to existing modelling approaches have been presented which can qualitatively reproduce more features of a typical Li-S cell behavior than could be reproduced with previous models.These improvements are easy to add to any existing Li-S model to ensure better agreement with the observed cell performance.

Conclusion
In summary, we have demonstrated a lumped mechanistic model for Li-S cells considering electrochemical and precipitation reactions, activation overpotentials, electrolyte resistance, as well as morphology variations during discharge.The model is the first to the authors' knowledge to qualitatively reproduce the evolution of electrolyte resistance reported in the literature.The change in electrolyte resistance during discharge can be explained as the result of the concentration dependence of electrolyte conductivity in conjunction with ionic concentration variations due to precipitation.The increase in activation overpotential in the low plateau can be explained by the reduced electrochemical surface area, also associated with Li 2 S precipitation.In addition, the limited rate of precipitation is shown to cause a drop in reduction potential, which manifests itself as a 'precipitation resistance'.In view of the central role that Li 2 S precipitation plays in determining the operating voltage of Li-S batteries, future efforts should seek to develop an accurate understanding of precipitation mechanisms, as well as to measure the polysulfide solubility, precipitation rates, and the concentration-conductivity relation.Fig. 3: Simulated activation overpotentials η 2 (for reaction Eq. 2) and η 6 (for Eq. 6) at 0.2C, and the specific electrochemical surface area (a v ) during discharge.
Fig. 4: The simulated reduction potential E 6 at 0.15C (solid line) and 0.03C (dashed line), and same potential calculated without considering Li 2 S precipitation (symbols).∆E 6 indicates the shift in E 6 caused by the difference in discharge current.

Fig. 1 :
Fig.1:(a) Electrolyte resistance (R s ) measured from a pouch Li-S cell manufactured by Oxis Energy Ltd. Galvanostatic EIS measurements were performed during a 0.34 A cell discharge, and the series resistance was extracted as the high-frequency x-axis intercept of the Nyquist plots at different depths of discharge.(b) Simulated R s from the Kumaresan model[9] for a Li-S cell with similar energy capacity; model parameters are taken from Ghaznavi and Chen[12].Fig. 2: (a) Electrolyte conductivity as a function of Li + concentration according to Eq. 16.(b) Simulated concentration of Li + during a 0.15C (0.34 A) discharge.(c) Simulated cell voltages and (d) simulated electrolyte resistances during discharge at 0.15 C and 0.03 C.Fig.3: Simulated activation overpotentials η 2 (for reaction Eq. 2) and η 6 (for Eq. 6) at 0.2C, and the specific electrochemical surface area (a v ) during discharge.Fig.4:The simulated reduction potential E 6 at 0.15C (solid line) and 0.03C (dashed line), and same potential calculated without considering Li 2 S precipitation (symbols).∆E 6 indicates the shift in E 6 caused by the difference in discharge current.

Figure 1 :Figure 2 : 1 Figure 3 : 6 Figure 4 :
Figure1: (a) Electrolyte resistance (Rs) measured from a pouch Li-S cell manufactured by Oxis Energy Ltd. Galvanostatic EIS measurements were performed during a 0.34 A cell discharge, and the series resistance was extracted as the high-frequency x-axis intercept of the Nyquist plots at different depths of discharge.(b) Simulated Rs from the Kumaresan model[9] for a Li-S cell with similar energy capacity; model parameters were taken from Ghaznavi and Chen[12].