Identifying Anode and Cathode Contributions in Li-Ion Full-Cell Impedance Spectra

: Measured impedance spectra of Li-ion battery cells are often reproduced with equivalent circuits or physical models to determine losses due to charge transfer processes at the electrodes. The identiﬁed model parameters can usually not readily or unambiguously be assigned to the anode and the cathode. A new measurement method is presented that enables the assignment of features of impedance spectra of full cells to single electrodes. To this end, temperature gradients are imprinted perpendicular to the electrode layers of a single-layered Li-ion battery cell while impedance spectra are measured. The method exploits different dependences of the charge transfer processes at the electrodes on temperature. An equivalent circuit model of RC-elements and the effect of temperature on the related electrode properties is discussed to demonstrate the feasibility of the method. A reliable assignment of the change of impedance spectra to the electrode processes is shown to be possible. The assignment can be used to identify if changes in an impedance spectrum originate from the anode or the cathode. than the full cell impedance. These studies show that the time constants of the main charge transfer reactions at anode and cathode are most often


Introduction
Measured electrochemical impedance spectra are usually evaluated with physicochemical models [1,2] or, more frequently, with equivalent circuits [3,4] to analyse internal processes of Li-ion battery (Lib) cells. The simulations use physicochemical parameters, such as reaction constants, or circuit elements, including the charge transfer resistance, to reproduce the measured impedance spectra of full cells. It is very difficult to distinguish between effects of the anode and the cathode and to assign them to specific features of a spectrum. Usually half-cell measurements are used for this purpose but the results found on half cells can only be applied to full cells to some extent.
Many spectra of full cells include a distorted arc that is usually assigned to a superposition of two semicircles which represent charge transfer impedances at the anode and the cathode [1,5]. Independent of the shape of the spectra, unambiguous allocation of the semicircles to their respective electrodes requires prior knowledge, e.g., via half-cell measurements or property variation of single electrodes. Correct allocation, however, is crucial if anode and cathode parameters of equivalent circuit models have to be identified unambiguously. Numerous studies have analysed half cells [6][7][8][9][10][11], symmetrical cells [12][13][14][15][16][17] or cells including a reference electrode [18][19][20][21][22] in order to investigate the impedances of an individual electrode rather than the full cell impedance. These studies show that the time constants of the main charge transfer reactions at anode and cathode are most often in the same order of magnitude, but not identical. This leads to the above-mentioned two overlapping arcs in an impedance spectrum which are somewhat separated from each other. To the best of our knowledge, no study has yet shown how to identify anode and cathode impedance contributions from full-cell spectra without using reference electrodes or modifying one electrode. Such information, however, will be extremely helpful to understand processes in the battery: changes of impedance spectra, e.g., due to a change of the State of Charge (SoC) or the State of Health (SoH), can be assigned to changes of the properties of the anode or the cathode [23]. Figure 1 shows measured impedance spectra of an anode half cell, a cathode half cell, and the corresponding full cell of a Li-ion battery cell. The spectrum is reduced by fitting and subtracting inductive and diffusive parts, so that only the charge transfer impedances of anode and cathode, represented by the arcs, are visible. The sum of the half-cell impedances equals the full cell impedance. In order to illustrate the issue, the full cell spectrum is simulated by an equivalent circuit model with two constant phase elements (CPE) each in parallel to a charge transfer resistance R, which we denote by CPE||R (solid blue line). Simulations on this simple model are able to reproduce the measured spectrum, but without half-cell measurements, it cannot be decided which of the two CPE||R elements must be assigned to the anode and which to the cathode. In this study, we present a new measurement method to address this problem by using Peltier elements to imprint a temperature gradient across the width of a singlelayered Li-ion cell. The basic idea to distinguish anode and cathode in the modelling of a full-cell impedance spectrum is to measure impedance spectra Z( f ) at different anode and cathode temperatures. Cooling down one side of the cell while heating up the other side decelerates processes in the cold electrode while processes in the warm electrode are accelerated. Both effects will be visible in the resulting full-cell impedance spectrum. A simple equivalent circuit model, consisting of two RC-elements in series has been used to analyse the corresponding full-cell impedance spectra after they have been reduced to the impedance of charge transfer processes. It will be demonstrated how to assign the semicircles of the individual charge transfer reactions to the anode and the cathode. With charge transfer reaction, we always mean the electrochemical reaction with subsequent intercalation or deintercalation. Troxler et al. [24] used a similar setup. However, they put their focus on multi-layered cells which did not allow them to heat or cool the two electrodes separately. Therefore, an assignment of electrode contributions to anode and cathode has not been done yet. We verify the robustness of our method with measurements of two batteries of the same type. The method is further supported by additional half-cell measurements. This investigation demonstrates a method to unambiguously identify the parameters of charge transfer processes at the electrodes in impedance-based models. We will demonstrate that the new method can be used to assign parameters of simulations of full-cell spectra to the anode and the cathode. The study will demonstrate the feasibility of the method by assigning the temperature-dependent relaxation times of two RC elements to the respective temperature-dependent semi-circles of simulated impedance spectra.

Materials and Methods
The experimental realisation has been done with the measurement setup depicted in Figure 2. A single-layered battery cell (graphite/NiMnCo-Oxide, NMC) has been placed between two Peltier elements. In Figure 2, the anode is heated up and the cathode is cooled down. Temperatures have been measured with two PT 100 temperature sensors, which have been positioned between the Peltier elements and the anode and cathode side of the battery cell, respectively. Cooling rips and fans have been attached to the Peltier element at the cold side to increase the achievable temperature difference applied across the width of the cell. All measurements have been performed in a temperature chamber at 23 • C to ensure stable temperature. The temperature gradient across the cell is adjusted by the set temperature of the chamber and the current applied to the Peltier elements.  Two 60-mAh high-energy density Lib cells produced by Custom Cells, denoted as cell I and cell II, with roughly 17 cm² NMC and graphite electrodes respectively have been analysed with a Gamry Reference 3000 potentiostat. Both cells had a Lithium titanate (LTO) reference electrode integrated for the half-cell measurements. Measurements for the full cell were conducted as well on the cells with reference electrode. All measurements have been performed at 50% SoC referred to the nominal capacity. Impedance spectra have been measured from 0.1 Hz-10 5 Hz with an amplitude of 6 mA. The cells have been stored at room temperature and 50% SoC between the measurements.
In a first measurement, the cathode was heated up while the anode was cooled down. An impedance spectrum has been measured after the measured temperatures had stabilised. This measurement is called Meas 1. It took about 2 min to stabilise the temperature measured by the Peltier elements. For the results shown, we have waited 30 min to ensure a steady temperature profile across the cell. Afterwards, the warm and cold sides have been switched by turning the cell by 180°and a second impedance spectrum has been measured, called Meas 2.
The actual temperature profile across the cell could not be measured but has been approximated from the temperature measurements at the battery surfaces. The Peltier elements were supplied with 6.6 A and 12 V, resulting in 74 • C at the hot side and 39 • C at the cold side.

Temperature Distribution
Using the temperatures at the battery surface boundary conditions, the temperature profile within the cell can roughly be estimated. The results depend on the thermal conductivities and the thicknesses of the pouch cell bag and of the electrolyte-soaked electrodes and the separator. Thermal conductivities for similar cell materials have been used from the literature. Based on the parameters presented in Table 1, the temperature distribution shown in Figure 3 has been calculated to roughly estimate the temperature distribution. Even though the temperature difference between anode and cathode is lower compared to the measured temperatures at the Peltier elements, it is still sufficient for our purpose because the impedance of a Li-ion cell is very sensitive to changes in temperature.  Figure 3 shows that the estimated temperature of the anode (59 • C) is 5 • C higher than the temperature of the cathode (54 • C). The thermal conductivities of anode and cathode are similar, and therefore, the temperatures of the anode and cathode will approximately be switched after the pouch cell is turned between the Peltier elements.

Data Evaluation Method
The procedure starts with measurement of Z Meas1 and Z Meas2 (Step1). This is followed by a modelling step (Step 2) and a model reduction step (Step 3), before the assignment of electrodes is done (Step 4). The impedance measurement results from Meas 1 and Meas 2, i.e., Z Meas1 and Z Meas2 , cannot be evaluated directly. In the upcoming sections, necessary simplifications and all steps of the data evaluation are discussed. The data evaluation strategy is displayed schematically in Figure 4.

Simplifications
In order to be able to evaluate the experimental data, assumptions about the underlying physical processes have to be made to simplify the analysis. We assume that charge transfer at the anode and cathode depends differently on temperature. As a consequence, we expect that the change of the impedance at the anode and the change of the impedance at the cathode are likewise different if the temperatures of both electrodes are switched. Therefore, we calculated the differences of the impedance spectra of the full cells Meas 1 and Meas 2 in terms of the moduli |Z( f )| and the phases ϕ( f ) at each frequency f : and Note that henceforth we omit writing the frequency dependence of the impedances throughout this paper to improve readability. Following the definition of Meas 1 and Meas 2, the differences represent the effect of a warm cathode that is cooled down, superimposed by a cold anode that is heated up.
Data evaluation is focused on the arc of the impedance spectrum that represents charge transfer across the electrolyte/electrode interfaces at the anode and the cathode. The shape of the measured arc is affected by several physicochemical processes and properties of the cell. Basically, it is formed by the charge transfer resistances and the double layer capacitances, which would result in two superimposed semi-circles. However, solid state diffusion affects the low-frequency part of the semi-circle and inductive and capacitive effects, induced by the wires or stray capacitances, can affect the high-frequency part of the semicircle. Additionally, geometrical and material inhomogeneities lead to a distribution of time constants of the RC elements which flattens the arc [29], and the solid electrolyte interface at the electrode surface may further distort the arc. As a consequence, a comprehensive equivalent circuit to model the processes at the electrode is rather complex. A typical equivalent circuit is shown in Figure 5a and the original spectra are given in Figure 5b. For the purpose of this investigation, we have to simplify the analysis of the spectra.
Reduction of measured spectra to charge transfer processes: Here, we start by modelling the spectra Z Meas1 and Z Meas2 with the equivalent circuit given in Figure 5a to obtain Z Sim1 and Z Sim2 (step 2). After simulation, we reduce Z Sim to extract the impedances of the charge transfer processes at the electrodes Z corM1 and Z corM2 . For this, we subtract impedance contributions caused by the internal resistance, diffusion, and induction. Hereby, the frequency range is unchanged. These contributions were calculated with a parameter identification algorithm in Matlab which minimises the difference between the measured absolute impedances and the sum of the assumed impedance contributions. The charge transfer impedances are represented by a resistor in parallel to a constant phase element for this purpose [4]. The equivalent circuit for this extraction process is illustrated in Figure 5a. We are aware of the fact that the surface layer on the electrodes contributes to the cell impedance as well [30]. However, the procedure works in both cases, with and without surface layer effects in impedance spectra. Figure 5b shows the measurement results for Meas 1 and Meas 2 and the same results after subtracting the impedance contributions of the internal resistance, diffusion, and induction (Z CorM1 and Z CorM2 ). The respective steps are depicted in the flow diagram in Figure 4a.
We have further simplified the investigation hereafter by replacing the CPE, used to perform the reduction of the spectra, with two RC elements representing the two electrodes. Each of them consists of a charge transfer resistance R in parallel to a capacitance C. The latter is usually attributed to the double layer at the electrolyte/electrode interface. This simplification is done as the effect of a temperature change on the charge transfer processes at the electrodes can be understood more easily in terms of the charge transfer resistance and double layer capacitance compared to the rather abstract coefficients of CPE elements which cannot easily be assigned to physical quantities. Note that as the equivalent circuit used to extract the electrode impedance (Figure 5a) is only an approximation of the actual impedances of the processes in the cell, the reduced spectra that are supposed to reflect just the electrode impedances, are subject to uncertainty. Likewise, the subsequent use of the simple RC elements introduces further uncertainty. These uncertainties are acceptable as the investigation only aims to identify the arcs of an impedance spectrum with respect to the anode and the cathode rather than providing an exact model of the cell impedances. We consider the reduced spectra and the simple RC model sufficient for the proof of concept in this manuscript, in particular, because we have substantiated the results by further validation steps in Section 3.1.

Assumptions for Temperature Influence
In the following, we discuss the expected effects of a temperature change on electrodes that are modelled by RC elements.
Nonlinear temperature effects: Firstly, charge transfer reactions at an electrode and diffusion in the vicinity of the electrode/electrolyte interface accelerate with increasing temperature which decreases the charge transfer resistance of the electrode. We have measured the resistances of the charge transfer reactions of the electrodes of two cells with half-cell measurements at various temperatures and found a nonlinear temperature dependence for the anode and the cathode. The results for temperatures between 5 • C and 26 • C are depicted in Figure 6. It can be seen that the change in resistance due to a change in temperature is higher for high resistances than for low resistances. It has to be noted that the resistance of the two cells differs. We attribute this to the fact that they have been made exclusively for the purpose of this investigation. Cells produced in smaller numbers usually have higher cell-to-cell variations. Since the two cells are analysed independently of each other, the difference can be accepted for this study. The assignment of the electrode impedance contributions will be done for both cells.
Temperature dependence of parameters: Secondly, we assume that the temperature change mainly affects the charge transfer resistance rather than the double layer capacitance. The effect of temperature on the double layer capacitance is difficult to assess. The double layer forms to compensate differences in the chemical potentials of the electrode and the electrolyte at the electrode/electrolyte interface. These potentials also depend on temperature. However, the relative dependence of their difference on temperature is assumed significantly smaller than that of the charge transfer resistance. This assumption is supported by the findings of Ahmed and Bade Shresthra [31]. We did comparable calculations as Ahmed and Bade Shresthra for the half-cell data of the cells used within this research and found a change of the double layer capacitance of only ≈5% for each electrode after changing the temperature from 5 • C to 17 • C. Temperature impact on time constants: Finally, we assume that the time constant related to the charge transfer process at the electrolyte/anode interface and that related to the electrolyte/cathode interface are different but they are in the same order of magnitude. In other words, the relevant frequency ranges of the involved impedances which form the arc overlap significantly, but they are shifted against each other to some extent. This is a common property of full-cell impedance spectra which has been shown with half-cell measurements [8].

Identification of the Electrodes in Full-Cell Spectra
Based on the assumptions and the approach introduced in the previous sections, we now describe the general procedure of how to identify the cathode and the anode in Z Meas1 and Z Meas2 , respectively. The measurements with the different temperature configurations Z Meas1 and Z Meas2 each include the impedance contributions of both, the anode and the cathode. As seen in Equation (1), this leads to a total of four impedances affecting ∆|Z| and ∆ϕ. This is exemplary shown in Figure 7 for cell I. Z Meas1 is modelled with two RC-elements Z RC1a and Z RC1b in series, representing the two electrodes. Likewise, Meas 2 is modelled with Z RC2a and Z RC2b . A nonlinear least square fitting procedure as proposed in [28] has been applied to minimise the difference between Z Meas1 and the sum of Z RC1a and Z RC1b as well as the difference between Z Meas2 and the sum of Z RC2a and Z RC2b . Please note that only the full-cell experiments were used, not the half-cell measurements. The simple model reproduces the experiments with sufficient accuracy for the analysis. Up to this point, Z RC1a , Z RC1b , Z RC2a , and Z RC2b cannot readily be assigned to the anode and cathode. This could easily be accomplished directly if the frequency ranges of the impedances contributing to the arcs do not overlap significantly. In this case, the investigation could be limited to one of the two frequency ranges. The RC-element with the impedance decreasing from Meas 1 to Meas 2 would be assigned to the anode and vice versa, because the anode has been heated up and the cathode been cooled down. However, usually, the frequency ranges overlap significantly and the identification is more complex so that it cannot be decided unambiguously how the semi-circles have changed. Therefore, further steps must be considered in order to assign Z RC1a , Z RC1b , Z RC2a , and Z RC2b to the electrodes and their corresponding temperature conditions. At low frequencies, the impedances of all four RC elements are dominated by the real part of the complex impedance. Hence, the behaviour of the real part of the impedances in this frequency range can aid in the identification process. The largest resistance must be assigned to a cold electrode, since its value must become smaller after heating up the electrode. Hence, if the largest resistance of Z RC1a , Z RC1b , Z RC2a , and Z RC2b corresponds to measurement 1, where the cathode is heated, the cold anode must be assigned to the respective impedance. This is Z RC1a in the example illustrated in Figure 7. Consequently, Z RC1b must obviously be assigned to the warm cathode. In turn, if the largest resistance corresponds to measurement 2, the cold cathode will be assigned accordingly. Next, the two remaining RC elements (Z RC2a and Z RC2b in our example) will be assigned which, at this point, cannot be uniquely attributed to cathode or anode since both fulfil the criterion that they are smaller than Z RC1a and larger than Z RC1b . To this end, the relaxation times τ = RC corresponding to each of the four impedances will be calculated from the resistances (R RC1a , R RC1b , R RC2a , R RC2b ) and capacitances (C RC1a , C RC1b , C RC2a , C RC2b ) that have been obtained from the fit. Since we assume that the double layer capacitance does not change significantly with temperature, we assign those impedances to the same electrode which has a similar capacitance. In our example, the values shown in Table 2 have been obtained. Thus, Z RC1a and Z RC2a are assigned to the anode and Z RC1b and Z RC2b are assigned to the cathode. The described procedure is schematically shown in the block diagram in Figure 4b.

Consistency Checks
Because of the uncertainties of the procedure that have been mentioned above, it is reasonable to conduct a few consistency checks to substantiate the results. First of all, it should be checked if the assigned impedance values comply with the expected temperature behaviour that we have described above. For frequencies f ≈ 0.1 Hz, Equation (3) can be approximated and rearranged such that ∆|Z| = (|Z cold anode + Z warm cathode |) − (|Z warm anode + Z cold cathode |) The first term on the right-hand side of Equation (5) expresses the change of the resistance at the anode and the second term that of the cathode. Consequently, the sign of ∆|Z| is positive if the change of the resistance of the anode is greater than the change of the resistance of the cathode and vice versa. Thus, the RC element showing the larger change of the resistance must be assigned to the anode if ∆|Z| is positive. If ∆|Z| is negative, it must be assigned to the cathode. Additionally, the cold electrode showing the larger resistance must show the larger change because of the nonlinear temperature dependency described above. Looking at Table 2, it can easily be checked that the assignments that we have made in our example comply with these requirements. ∆|Z| is positive at the low frequency end (calculated from Equation (3)); hence, R RC1a − R RC1b must be larger than R RC2a − R RC2b . This is schematically shown in Figure 4b.
Finally, it should be verified if the measured frequency behaviours of ∆|Z| and ∆ϕ are consistent with the modelled frequency behaviour of the impedances of the RC elements and their assignments to the electrodes. This should further substantiate the results since the investigations have been based on the Nyquist plots so far, which give, besides the relaxations times, no direct information on the frequency behaviour of the involved impedances. To this end, ∆|Z| and ∆ϕ will first be calculated from Z RC1a , Z RC1b , Z RC2a , and Z RC2b , and the assignments of these impedances to the warm and cold electrodes are performed, using Equation (3). The main features are depicted in Figure 8, which shows the moduli of Z 1 , Z 2 , and ∆|Z| on the left-hand side and their phases and ∆ϕ on the right-hand side. The features of the individual RC elements (not shown) are well known. The moduli of the impedances start from the associated resistance values at low frequencies and level off to zero at the high-frequency end, while the corresponding phases change from zero to π/2. The moduli and the phases of Z 1 and Z 2 show a similar behaviour. However, ∆|Z| and ∆ϕ can show quite different features, depending on the values of the involved resistances and capacitances. Figure 8 shows the features under the conditions that we have assumed, i.e.,

•
The resistance values are in the same order of magnitude; • The largest resistance is assigned to a cold anode; • The change of the resistance at the anode and the cathode is smaller than the difference of the anode and cathode resistance; • The double layer capacitance of the anode and the cathode are in the same order of magnitude; • The double layer capacitance at an electrode does not change significantly with temperature.
The main features that can be observed are: • ∆|Z| first decreases with increasing frequency. It is positive at the low-frequency end because RC1a has been assigned to the cold anode, i.e., the highest resistance. With increasing frequency, the moduli of Z 1 and Z 2 cross each other. As a consequence, their difference becomes negative, and then, it shows a minimum. At high frequencies, ∆|Z| gradually levels off to zero. If the highest resistance is assigned to the cold cathode, ∆|Z| will be negative at the low frequencies and it will show a maximum. • ∆ϕ approaches zero at the low-and high-frequency end. It has a minimum at lower frequencies and a maximum at larger frequencies.
At the high-frequency end, one could expect that ∆Z and ∆ϕ are dominated by the processes showing the smallest relaxation times, i.e., the warm electrodes since the slower processes do contribute less to the total impedance in this region. However, we have found that the significant overlap of the frequency ranges of the four impedances involved, including eight parameters, make it difficult to assign specific features seen in the highfrequency range to the warm electrodes. In the next section, we apply the described data evaluation procedure to our measured data.

Method Application to Experimental Data
We have performed the procedure described in Sections 2 and 2.2 with two cells of the same type in order to verify the reproducibility. The spectra already shown in Figures 5b  and 7 depict the results of cell I. The spectra of cell II (not shown) have a comparable behaviour. Tables 2 and 3 list the identified values of the four RC impedances Z RC1a , Z RC1b , Z RC2a , and Z RC2b for cell I and cell II, respectively. Following the reasoning of the simulations made in Section 2.2.3, Z RC1a and Z RC2a of the measurement have to be assigned to the anode and Z RC1b and Z RC2b have to be assigned to the cathode. The absolute resistance and capacity values of the two cells differ to some extent due to manufacturing variability and possible additional errors from the simplifications described in Section 2.2.2. Figure 9a,b shows the frequency behaviour of ∆|Z| and ∆ϕ of cell I and cell II, respectively. Qualitatively, their frequency behaviours fit well to the simulation results illustrated in Figure 8. ∆|Z| is positive at the low-frequency end, decreases with increasing frequency, changes its sign in the mid-frequency range and shows a minimum before its starts to level off. The phase difference shows a minimum and rises above zero to another maximum, which seems to be beyond the investigated frequency range for cell II. Therefore, the analysis of the frequency dependency of the measurement data supports or, at least, does not falsify the conclusion that the anodes in the measured cells have a higher charge transfer resistance and a higher time constant than the cathode. The frequency behaviour of ∆|Z| and ∆ϕ of both cells is qualitatively the same which substantiates the robustness of the procedure qualitatively.

Validation of the Assignments with Half-Cell Measurements
Half-cell measurements have been additionally performed to verify the attribution of electrodes and thus our methodology. To this end, the anode and cathode impedances have been simultaneously measured with the impedances of the full cell to determine the relaxation times of the related charge transfer impedance independently. It should be noted that the sum of the imaginary parts of the half-cell measurements equals those of the corresponding full-cell measurements, which also holds for the sum of the real parts. This proves that the three-electrode arrangement cell can be used to split the measured full-cell impedances into the contributions of the individual electrodes. The untreated results are shown in Figure 10a,b which depicts the imaginary parts of the impedances over the frequency. The reciprocal of the characteristic frequency at the minimum of the imaginary part corresponds to the relaxation time of the charge transfer impedance of an electrode, assuming that it is determined by the behaviour of an RC element as described in Section 2.2.1. Table 4 lists the relaxations times of the electrodes under warm and cold conditions and compares them to those derived by our method. The relaxation times derived from the half-cell measurements confirm the assignment of relaxations times to anode and cathode from the full-cell measurements. The differences in absolute values may be attributed to our simplifications as mentioned above. Nevertheless, the half-cell measurements support the simulated results, in particular, the result that the anode has a higher charge transfer resistance and a higher relaxation time.

Conclusions
In this work, we have presented a new method which enables to assign anode and cathode impedance contributions solely with measurements on full-cell LiBs. A temperature gradient is applied on a single-layered Li-ion battery cell and exploits different dependences of the anode and cathode charge transfer processes on temperature. Feasibility has been demonstrated by reproducing the measurement results with a simple model of RC elements showing that the method is suitable to decide if the identified parameters of an electrode process used for the simulation should be assigned to the anode or the cathode. The assignment has been substantiated by further consistency tests and been proven with half-cell measurements. Knowledge about the anode and cathode contributions to the overall full-cell behaviour improves the possibility to evaluate various kinds of characterisation measurements. Impact of changes of the measurement conditions, i.e., temperature, SoC or SoH can be assigned to a specific electrode more easily if the introduced method is applied in advance. The method might be optimised by increasing the temperature gradient across the cell, by decreasing the effects disturbing the electrode impedances, or by measuring the impedances at various temperatures. Future studies which comprise electrode manufacturers may also target to proof all cases shown in Figure 4, while this study is just be a proof of the principle. However, once the impedance contribution of a specific process is assigned to its related electrode, this information can be used for more sophisticated models to simulate the impedance spectra of other cells of the same cell chemistry, cell balancing, and electrode particle sizes. The introduced method might not be limited to Lithium-ion cells but could also be suitable for, e.g., Li-S cells [32,33].

Conflicts of Interest:
The authors declare no conflicts of interest.