Determination of the Distribution of Relaxation Times by Means of Pulse Evaluation for Ofﬂine and Online Diagnosis of Lithium-Ion Batteries

: The distribution of relaxation times (DRT) analysis of impedance spectra is a proven method to determine the number of occurring polarization processes in lithium-ion batteries (LIBs), their polarization contributions and characteristic time constants. Direct measurement of a spectrum by means of electrochemical impedance spectroscopy (EIS), however, suffers from a high expenditure of time for low-frequency impedances and a lack of general availability in most online applications. In this study, a method is presented to derive the DRT by evaluating the relaxation voltage after a current pulse. The method was experimentally validated using both EIS and the proposed pulse evaluation to determine the DRT of automotive pouch-cells and an aging study was carried out. The DRT derived from time domain data provided improved resolution of processes with large time constants and therefore enabled changes in low-frequency impedance and the correlated degradation mechanisms to be identiﬁed. One of the polarization contributions identiﬁed could be determined as an indicator for the potential risk of plating. The novel, general approach for batteries was tested with a sampling rate of 10Hz and only requires relaxation periods. Therefore, the method is applicable in battery management systems and contributes to improving the reliability and safety of LIBs.


Introduction
The use of energy storage systems is essential for the transition to renewable energies. Due to the unsteady energy supply from renewable energies such as wind power and photovoltaics, energy storage devices are required to compensate for fluctuations in output and to reliably provide energy at all times. For the electrification of the transport sector, energy storage systems are required in order to move vehicles independently of an external energy supply and without burning fossil fuels. In battery electric vehicles in particular, LIBs are by far the most frequently used battery technology due to their high energy and power density.
Large batteries are required in electric vehicles (EVs) and stationary energy storage devices, which are therefore the main cost factor in these applications. Accordingly, the requirements for the service life of the LIBs are high and a high degree of reliability is a prerequisite.
Therefore, non-invasive investigation methods are used to achieve a comprehensive understanding of the underlying degradation mechanisms of LIBs to enable system optimizations with regard to the service life. In addition, the non-invasive methods are required for the online diagnosis of LIBs and to identify potential risks during operation in order to ensure safe operation.
In current studies, the continous determination of the changes in kinetic parameters of LIBs made it possible to draw conclusions about progressive aging and degradation mechanisms [1][2][3][4][5]. Compared to recording of the open-circuit-voltage (OCV) or to directly determining the capacity, measuring the dynamic behavior of LIBs is usually less timeconsuming [1].
EIS as a non-invasive technique is widely used for aging studies in order to measure the impedance directly and thus to determine the kinetic parameters of the LIBs [6][7][8][9]. Aging processes such as lithium plating [3] and the formation of the solid electrolyte interface (SEI) [4] could be investigated using EIS.
The derivation of the DRT on the basis of the EIS measurement data as well as the further DRT analysis is a common method for the investigation of LIBs and other batteries. In [10][11][12][13][14][15] methods are investigated and described to obtain the DRT from impedance data. In several publications, the number of dominant electrochemical processes of LIB as well as their time constants and polarization contributions could be quantified by analyzing the DRT of the measured spectra [14][15][16][17]. Changes in the time constants and polarization contributions during the aging of the LIBs under investigation can also be determined using the DRT, as recently shown by Sabet et al. [5].
The main disadvantage of EIS is the long measurement time, especially for low frequencies. The measurement time ranges from several hours for frequencies in the millihertz range to several days for frequencies in the microhertz range. Therefore the bandwidth and resolution of the impedance measurement data available for the DRT is limited [15]. Additionally, only small signal amplitudes are permitted to ensure that the steady state condition is not violated during the EIS. The long measuring times and low signal amplitudes cannot be implemented under operating conditions and with commercially available measuring electronics [18], which is why DRT analyzes have so far been limited to laboratory tests. But even under laboratory conditions, when measuring lower frequencies, maintaining the steady state and current drift are becoming increasingly critical [19].
Due to the general availability in online applications, the characterization of cell impedance based on pulse data is of great interest and the subject of some current studies (e.g., [20,21]). The time domain data of current and voltage curves can be evaluated in the event of pulse excitations in order to determine the impedance of the LIB and other battery types. The calculation on the basis of time domain data is generally faster than the direct measurement of the impedance, since several frequencies are excited simultaneously with excitations in the time domain [19]. Thus, on the basis of time domain data, a higher frequency resolution and wider bandwidth of the impedance measurements can be achieved, which is of particular interest with regard to the DRT.
An overview of the options for calculating the impedance using time domain data is given in [22]. In [19,23], the time domain data resulting from pulse measurements were brought into the frequency domain using a fast Fourier transform (FFT) in order to calculate the low frequency impedance. However, the excitation of a pulse is nonperiodic. A window function must therefore be used for the transformations, which in turn generates an offset of the signal [18]. To avoid additional sources of error, as an alternative to the FFT or Laplace transformation, a pre-selected equivalent circuit model (ECM) can be parameterized using the time domain data [19]. The impedance is obtained by transforming the ECM into the frequency domain. The ECM is chosen in a manner that the impedance behavior of the battery can be reproduced. Therefore, the model must be selected individually for each battery type, cell chemistry and the electrochemical processes to be expected [24]. In addition, the impedance of the batteries can change significantly with aging, state of charge (SOC) as well as temperature and thus also the dynamic behavior in the time domain. For this reason, prior knowledge of the electrochemical processes involved and of the aging behavior and further dependencies is required in order to select a suitable model.
The DRT can be derived from the calculated impedance. If the selected ECM consists only of RC elements (parallel connection of resistance and capacitor) and an internal resistance (Thévenin model), the DRT is given directly by the model in the time domain and a transformation to frequency domain is not necessary. Analytical methods have already been proposed for the online parameterization of a Thévenin model [20,25] but the number of RC elements is limited and too small to be able to derive a DRT with adequate resolution. The sufficient number of RC elements can be characterized by iterative numerical optimization methods such as genetic algorithms and method of least squares. Compared to the analytical methods, however, the computationally intensive iterative calculations are disadvantageous for online applications [25] and the number of RC elements is limited again due to the computing power available [26]. The initial values to be determined also require prior knowledge of the electrochemical processes and their approximate time constants and polarization contributions.
In summary, it can be stated that based on the analysis of the DRT, aging mechanisms have already been examined in numerous investigations. The DRT is determined according to the state of the art on the basis of the impedance. The direct measurement of the impedance by means of EIS cannot be implemented in current online applications without adaptations and at low frequencies, even under laboratory conditions, it is time-consuming and subject to inaccuracies. The indirect determination of the impedance based on the evaluation of time domain data leads to additional inaccuracies or requires an individual definition of a model for different cell chemistries and types. A direct derivation of the DRT based on time domain data, which can also be used online, would therefore be desirable.
In this study a method is proposed and described to derive the DRT directly from the measured voltage course during cell relaxation. The method is validated experimentally by comparing the results with a DRT obtained from EIS measurements. In an aging study, the introduced method and a subsequent DRT analysis are used to identify the time constants and polarization contributions of the electrochemical processes of LIBs. The method is shown to be sufficiently sensitive to quantify the changes in these parameters during aging without the need for EIS equipment or other devices not available in online applications. The study demonstrates that the DRT derived from time domain data offers significantly more insights into processes with large time constants respectively low-frequency impedance with less expenditure of time than a DRT derived from EIS measurements. In addition, a correlation between a polarization contribution of one of the processes with a large time constant and a degradation mechanism could be established. Thus, the proposed method offers a time-effective estimate of the kinetic parameters of a LIB, provides insights into battery degradation and is offline and online applicable.

Methods
Multiple electrochemical processes with a wide range of time constants occur in LIBs. In addition to the the charge transfer processes of the individual electrodes, diffusion processes such as solid state diffusion and diffusion of lithium ions through the surface films of the electrodes also influence the dynamic characteristics of LIBs [9]. Some processes have similar time constants and are therefore not easy to distiguish from each other only considering the spectra obtained by EIS. In Figure 1a, a spectrum of the cells examined in this research is shown.
The diffusion branch at low frequencies can be clearly assigned to the solid state diffusion as the time constant of these processes is usually a few decades larger compared to other dynamic processes [9]. The semicircle in the spectra is caused by charge transfer processes and the diffusion trough the SEI [27]. Since the time constants are of a comparable order of magnitude, these processes show overlapping effects in the impedance spectra and only a single semicircle is visible. Changes of the polarization contribution and time constants of the various processes during aging can therefore not be reliably separated and assigned to the specific processes. An approach to identify the single processes and their parameters is an analysis of the DRT [14].

DRT of Frequency Domain Data
Since the DRT analysis of frequency domain data is well-established and the focus of this study is the DRT of time domain data the method is only briefly described. We refer to [15] for a more detailed explanation of the procedure used. In order to calculate the DRT using frequency domain data, the resistive-capacitive part of the given spectra is reconstructed with a series of m RC-elements, a parallel circuit of an ohmic resistor and a capacitor.
where h(τ k ) are the RC-elements unknown polarization contributions. By pre-defining the time constants τ = R · C, the equation of the RC elements and thus the optimization problem become linear. In view of the nature of electrochemical processes, an equidistant distribution on a logarithmic scale is assumed for the time constants. The range of time constants is expanded by two decades and the number of RC-elements per decade is increased by a factor of three compared to the measured frequencies in order to enhance the accuracy and to obtain an adequate resolution of the DRT (see e.g., [15]). Since the number of RC-elements per decade is higher than the resolution of the impedance data, the optimization problem becomes ill-posed. Therefore, a regularization technique which was introduced by Tikhonov [28] is used to enable an analytical solution. The Tikhonov regularization for determining the DRT has already been shown in some publications [11,[13][14][15]. The regularization parameter was elected according to [15]. The value is kept constant for the evaluation of all spectra, since the results of the DRT analysis are sensitive to the parameter. Thus, the DRTs of the different spectra are comparable and changes in the dispersion due to the aging process which can be tracked.
The parts of the measured spectra with positive imaginary values cannot be modelled by a series of RC-elements and have to be discarded. Additionally, the internal resistance R i is subtracted from impedance values of the measured spectra prior to the fitting procedure. It is assumed that R i is equal to the zero crossing of the imaginary axis [9]. After preprocessing of the measurement data, the spectra mainly consist of resistive-capacitive contributions, which can be modelled by RC-elements.
The optimization problem is solved with a non-negative least squares fit because only positive resistances are physically meaningful. Figure 1b depicts the DRT of the spectrum shown in Figure 1a. The position of the peaks indicate the time constants of the processes involved and the area under the peaks correlates to the polarization contribution of the process [14][15][16]. In contrast to the Nyquist diagram, which shows only a single semicircle at moderate frequencies, various processes with different polarization contributions and time constants are visible in the same frequency range.  For the proposed method of an DRT analysis of time domain data, the voltage relaxation is evaluated. Similar to the DRT analysis of frequency domain data, an attempt is made to reconstruct the voltage behaviour with a series of RC-elements. The voltage course of a single RC-element during a relaxation period can be described as follows: where τ RC = R · C, the characteristic time constant of a RC-element. The voltage of the RC-element at the beginning of the relaxation period U 0 at t = 0 s can be calculated if the duration t cp and current I cp of the previous pulse are known: Equation (4) is only valid if the cell has been sufficiently relaxed before the current pulse and any overvoltages can be neglected. Assuming a finite number m of RC-elements, the overall voltage response can be computed numerically: By defining fixed time constants τ k , the optimization problem becomes linear. The values of the resistances R k must be set in such a way that the error between the voltage prediction by the series of RC-elements and the measured voltage data (or in this case the simulated voltage data) is minimized. Using the sum of squared errors as a cost function leads to the following optimization problem min{J}: where the vector R vec corresponds to the unknown polarization contributions of the RC-elements and U vec is the vector of the measured (or simulated) voltage course: The number n corresponds to the number of measurement points. The measurement data should be preprocessed in advance. In Section 2.2.3, a detailed description of the process is given. The dimension of the matrix A is n · m and the matrix is calculated for the different predefined constants and the time after the end of the pulse t which corresponds to the respective voltage measurement in U vec .

Predefinition of Time Constants τ
A logarithmically uniform distribution is defined for the time constants [τ 1 · · · τ m ]. The minimum and the maximum time constants τ min and τ max are determined based on the maximum sampling rate f s,max = 1 ∆t min and the duration of the relaxation process t max . The Shannon theorem states that frequencies that are higher than half of the sampling rate cannot be evaluated without further processing of the measured signal.
Taking into account the transient frequency f t = 1 2·π·τ of a RC-element, the following relationship result for the smallest observable time constant τ eval,min [29]: The largest evaluable time constant τ eval,max is limited by the duration of the relaxation phase according to [19].
The minimum and maximum time constants of the distribution τ min and τ max are chosen to be two decades smaller, respectively, andlarger than the evaluable time constants in order to increase the accuracy at both ends of the frequency dispersion: This procedure was similarly proposed by [15] for the DRT analysis of frequency domain data. In order to obtain a smooth DRT curve within the evaluated frequency range, the number of RC elements or time constants m must be selected to be high enough. A number of one hundred τ per decade was used for the analysis of the DRT in this study.

Pre-Processing of Measurement Data
Before calculating the DRT, the existing measurement data should be checked for usability and preprocessed. According to the Equation (12), the sampling rate must be high enough (or the timestep between two successive voltage measurements ∆t min small enough) to enable a meaningful evaluation of the DRT down to the time constant τ eval,min . Therefore, when selecting the the sampling rate beforehand, the minimum time constant relevant for the investigation should be taken into account. If the maximum sampling rate of the measurement device is lower, the DRT can only be evaluated for time constants that satisfy Equation (12). This consideration is particularly important for BMS, as they often have sample rates of 10 Hz or less.
The processes with low time constants only have measurable contributions at the beginning of the relaxation phase and quickly subside. Therefore, the sampling rate does not necessarily have to be constant over the entire measurement period and can be reduced in order to avoid large amounts of measurement data. With very high sampling frequencies and the evaluation of very small time constants, the speed of the current control needs to be considered. As long as the current has not dropped to zero, Equation (5) is not valid. These voltage measurements should therefore be discarded and only larger time constants evaluated. Apart from that, only the voltage values should be evaluated that were recorded at least one time step after the load drop. Since pure ohmic resistances cannot be modelled with Equation (5), the influence of internal resistance on the voltage response should not be taken into account.
The OCV U OCV must be subtracted from the voltage measurements u meas (t) before determining the DRT. Assuming a completely relaxed cell at the end of the relaxation phase, this can be achieved by subtracting the last measured voltage value u meas (t max ) from all measurements. A considerably long relaxation phase is therefore imperative.
At the beginning of the relaxation process, all processes contribute to changes in cell voltage. As relaxation progresses, only processes with larger time constants are relevant. Following the measurement data should be interpolated on a logarithmic scale in order to weigh the single measurement points accordingly. Due to the low impedance of automotive cells, the signal-to-noise ratio is usually rather low. As the cell becomes increasingly relaxed, the signal-to-noise ratio deteriorates further, as can be seen in Figure 2. Forming the moving average before interpolation smooths the measured values and can therefore improve the ratio.

Calculation of the DRT
Due to the long measurement time for low-frequency impedances, the number of time constants m usually exceeds the number of frequencies measured. In contrast to frequency domain data, the resolution in the time domain is comparatively high. The number of interpolation points can also be selected to be higher than m in order to use the additional information provided by the high-resolution pulse measurement. To solve the over-determined system, the Tikhonov regularization is applied again. Equation (6) is extented by the regularization term: where λ is the regularization parameter and I is the m · m indentity matrix. The value of the regularization parameter is optimized with respect to the sum of the square errors of the reconstructed voltage signal. A non-negative least square algorithm is proposed to solve the optimization problem. Due to the restriction to positive results, only physically possible polarizations can be calculated. Using the determined polarizations R vec and the predefined time constants τ, the spectrum can be reconstructed for different angular frequencies ω: If a measured or modelled spectrum is available, the sum of the squares of errors of the capacity-resistive part can be calculated and used as an additional selection criterion for the regularization parameter. This method for the evaluation of the parameter was already proposed by [15] for the DRT calculation of frequency domain data.
In Figure 2b, the DRT is given; this was derived from the voltage course during the relaxation period shown in Figure 2a. The DRT reveals three processes vizualized by the peaks. The area under the peaks corresponds to the polarization contribution of the process and the position of the peak to the time constant. The time constants as well as the polarizations of the identified processes are consistent with the parameters of the RCelements used in the battery model. Table 1 shows the calculated values and the original parameter set of the battery model as well as the relative deviation. The relative error of the parameters with larger time constants decreases due to the diminishing influence of the superimposed noise signal. The processes with short time constants only contribute to voltage changes for a short time and their voltage contribution quickly drops to zero. The voltage drop is also most pronounced in processes with the long time constants at the beginning of the relaxation, but it lasts longer. Therefore, the effect of these processes is visible at more measuring points, which enables an increasing averaging of the noise. In Figure 3a, the reconstructed voltage signal is plottet together with the simulated voltage course in order to prove the accuracy of the DRT. Figure 3b shows the reconstructed impedance, derived by using the calculated distribution function and Equation (19), alongside the impedance spetrum of the battery model. Regardless of whether the real part or the imaginary part is considered, the relative deviation between the two spectra of Figure 3b is less than 1% for the frequency range [ f eval,min · · · f eval,max ]. In general, the accuracy is higher for lower frequencies within the specified frequency range. Because of the added noise in the simulation, the relative error of the reconstructed voltage signal (Figure 3a) cannot be used as an indicator for the quality of the distribution function. However, the reconstructed voltage follows the mean value of the simulated signal, even with dynamic behavior in the first seconds. Despite the good reconstruction of the dynamics, the impact of the noise is substantially minimized, which proves that no overfitting has taken place.
Based on the sampling rate of commercially available BMS, the resolution of the simulated values was only 100 ms. Since the impedance of automotive sized cells is usually very low and hence the voltage response, a strong noise signal has been added compared to the excitation. Nevertheless, processes with time constants that are only slightly larger than the resolution can be identified and their contributions determined. The results therefore show that the method introduced is theoretically suitable for determining the DRT and the cell impedance and examining the dynamics and processes parameters of an LIB.

Experimental
In order to further investigate the applicability and accuracy of the introduced method under real conditions, an experimental study was carried out. To validate the introduced method, automotive LIBs were characterized with both EIS and pulse tests and the resulting DRT and derived process parameters were compared.
The automotive cells were then cyclically aged under various conditions. The changes in processes and their parameters were determined using both methods. Thus, the results obtained through time and frequency domain data can be directly compared and evaluated with regard to the traceability of degradation mechanisms.
For experimental investigations, large format pouch-type LIBs with a nominal capacity of 50 Ah were used. The anode of the examined cells are made of graphite and the cathode consists mainly of nickel-manganese-cobalt-oxide (NMC). All cells were taken from the same batch to minimize the impact of manufacturing tolerances.

Experimental Validation
For the experimental validation of the method introduced one of the cell was fully charged and subsequently discharged to determine the available cell capacity. In order to minimize the influence of the cell impedance, a constant-current-constant-voltage (CCCV) protocol was used for both charging and discharging. The CV phases were held until the current dropped below 50 mA. In the CC phase, 25 A were applied during charging and 50 A during discharging. According to the determined capacity the cell was then charged with a constant current of 25 A to certain SOC values (20/40/60/80%). After reaching the desired SOC level, the current was set to zero. The relaxation of the voltage was recorded in high resolution with a sampling rate of 2 MHz in the first six seconds after the load jump. After six seconds, the sampling rate was reduced down to 1 Hz and the voltage relaxation recorded for 4 h.
The pulse tests were followed by EIS measurements. Galvanostatic excitation with an amplitude of 3 A was applied. The frequency was varied in the range from 50 kHz to 0.5 mHz. For frequencies above 100 Hz, 16 measurements per decade were carried out, otherwise only eight were performed. According to Barai et al. [30] the relaxation phase of the previous pulse test is long enough to meet stationary conditions.

Aging Study
As part of the aging study, the LIBs were cycled under certain operating conditions that are within the manufacturer's cell specifications. A total of three aging scenarios were carried out, with the average SOC value being varied between the individual scenarios. With the exception of the charging current and the state of charge, the other operating conditions like the depth of discharge (DOD) (30%) and the discharge current (1.6 C) were set to be the same for all scenarios. Table 2 gives an overview about the deviating operating conditions for each scenario. The SOC range of each scenario was set in such a way that, according to [31], different degrees of capacity depletion and deterioration of the individual electrodes could be expected. The aging study thus enables an assessment of whether the sensitivity of the introduced method is sufficient to determine aging mechanisms qualitatively or even quantitatively.
Before starting the cycling, a checkup was carried out for each cell in order to determine the initial cell parameters. During the cycling of the cells, the long-term tests were occasionally interrupted and further checkups were carried out. The checkup procedure contains a CCCV capacity measurement. In addition, at 20/40/60 and 80% a pulse test and an EIS measurement are performed, following the procedure described in Section 3.1. However, the maximum sampling rate during the pulse test is limited to 10 Hz, which is more realistic for online applications such as battery management systems.

Experimental Validation
In this section, the results of the pulse tests and the EIS measurements of the experimental validation are shown and evaluated. Figure 4 depicts the recorded cell voltage during the relaxation at 20%. In addition to the originally measured signal, the moving average of the signal is given in order to be able to assess the accuracy of the reconstructed signal. It can be stated that the deviation between the reconstructed and the averaged signal over the entire range under consideration is less than 0.5 mV. The high measurement resolution reveals the inductive behavior of the cell voltage. Due to the inductance, the voltage over the cell connectors and current collectors drop to negative values immediately after the load jump. For less than 0.1 ms, the voltage rises again as the voltage across the inductive components drops quickly. Therefore, the voltage measurements for t < 0.1 ms have been omitted, as electrochemical processes (time constants of usually larger than 1 ms) are of interest for this study. However, it is pointed out that the inductive behavior can also be simulated by expanding the matrix A (Equation (9)) by the mathematical description of several parallel circuits of ohmic resistance and inductance (RL-elements).
Equation (19) is used to calculate the impedance of the cell for the pulse test at 20% percent. For better comparability with the spectrum measured directly after the pulse test, Figure 5a,b show the imaginary and real parts of both spectra over frequency.  Figure 5a shows that the imaginary part coincides with frequencies below 50 Hz. Due to the limited measuring range, only impedances at frequencies above 0.5 mHz can be specified for the measured spectrum.
The plot of the real part of the impedance in Figure 5b provides similar insights. The course of the reconstructed spectrum deviates at high frequencies, but follows the measured spectrum at frequencies of ≤5 Hz.
The increasing deviations at >5 Hz can be explained by the Butler-Volmer kinetics. The current excitation at relaxation period was set to zero, while the excitation during the direct measurement of the impedance was set between −3 A and 3 A. According to the Butler-Volmer equation, the kinetics of charge transfer processes are not linear. Thus, due to increasing charge transfer resistances, the potential response does not decrease linearly with the current excitation [32]. In addition, it has been experimentally proven in some publications that the relaxation time, which differs between the pulse test and the EIS, influences the charge transfer resistances [30,32]. It is known that charge transfer processes at the anode [4,15] and at the cathode [5,33] occur at moderate frequencies of 1-100 Hz. The deviations at frequencies of 5-100 Hz between pulse test and EIS were therefore to be expected.
The ranges of the spectrum at frequencies above 100 Hz could not be adequately reconstructed because the optimization function is limited to resistive-capacitive elements.
The impedance values measured at over 100 Hz were therefore discarded and not plotted in Figure 5. Figure 6 shows the DRT obtained from time domain data and the DRT determined from frequency domain data. The resolution of the DRT from time domain data is twice as high. Therefore, the heights of peaks in one DRT cannot be compared directly with those in the other DRT. In order to compare the polarization contributions of the identified processes, the sum of all polarization contributions within a peak must be determined instead. When calculating the polarization contributions for the first (peak at smallest time constant) to fifth peak, the same orders of magnitude result for the respective peak for both DRTs. In addition, both DRTs show that the contribution of the first and fifth peaks are greatest, followed by the fourth and second peaks. The third peak has a comparatively small contribution in both DRTs. Both distributions reveal four relevant processes with time constants in the range of 2 · 10 −3 s < τ < 2 · 10 0 s. According to f T = 1 2·πτ , these time constants approximately correspond to frequencies between 0.1 Hz and 100 Hz. Processes in this frequency range were assigned to charge transfer [4,5,15,33] and surface films of the electrodes [4,34]. Due to the Butler-Volmer kinetics, time constants and polarization contributions differ between the two distributions, which explains the discrepancy that occurs in the first four peaks.
For time constants >10 s, the DRT obtained from time domain data shows a comparatively better resolution of the processes. The DRT based on frequency domain data only indicates a single process for time constants between 10 and 1000 s. In contrast, the DRT obtained through time domain data reveals the existence of three different processes. In addition, time constants that are two orders of magnitude higher can be resolved. If the spectrum is recorded via EIS, the measurement takes more than a day to resolve large time constants. The evaluation of pulse data with the method introduced can therefore be advantageous for examining processes with large time constants such as solid state diffusion.

EIS Measurement Evaluation
In Figure 7 shows the measured spectra and the corresponding DRT against frequency of an examined cell as well as the development during aging. The measurements from one cell from aging scenario Sc 1 were selected for the exemplary representation. The upper diagrams show the values determined at 20% and the lower part at 80% SOC. In the Nyquist diagram, two major changes can be seen as the cell ages: • An increase of the internal cell resistance R i during aging (Shifting of the zero crossing of the imaginary axis to higher real parts) • Disappearance of the second semicircle at 20% SOC during the first 1050 cycles of the long-term test According to the second bullet point, the DRT at 20% reveals distinct changes in the polarization contribution and the time constant for a process at moderate frequencies. This peak can be assigned to the second semicircle in the spectrum because it occurs particularly when the SOC is low and the cell is new. As already mentioned, this peak can be assigned to a charge transfer process due to the range of the time constant.
In contrast to the observation of the spectra, the analysis of the DRT shows a slight increase in the contribution of this process during the last cycles. Furthermore, slight increases in the polarization of the charge transfer processes are visible at moderate frequencies and 80%. In addition, small changes in the parameters of the process with the smallest time constant are visible, but they do not follow a clear trend. Due to the frequency range, this process can be assigned to the transfer through the SEI of the anode [4,34].
For frequencies <0.1 Hz, a single process is visible. At low SOC, the polarization of the process remains almost constant during the first 4950 cycles. After the first 4950 cycles, a continuous increase can be observed. At 80% SOC, this continuous increase begins with the start of the long-term test.
The mentioned changes in the process parameters were recognizable for all cells, except for the increase in polarization contribution at low frequencies. This behavior was not pronounced in the cells of scenario Sc 2 . Changes of the process parameters were made at a similar rate within a scenario. However, there were deviations in the rate of change between the scenarios, especially with regard to the increase in polarization contribution at frequencies below 0.1 Hz. Figure 8 shows the cell voltage course during the relaxation period and the DRT derived from the time domain data. Again, the values at 20% (upper graphs) and 80% SOC (lower graphs) are given and the same cell as in Section 4.2.1 is evaluated. The voltage courses show a slowdown of the relaxation process as aging progresses. The DRT derived by the time domain data reveals changing parameters of processes with large time constants as the reason. A strong increase of polarization contributions can be observed for low frequencies. Compared to the DRT obtained from frequency domain data, not only a more pronounced increase can be observed. In addition, more processes with large time constants can be identified, since the DRT can also be specified for frequencies that are lower than 1 mHz.

Pulse Test Evaluation
According to Equation (10), the DRT only provides representative values for the frequencies ≤5 Hz. The DRT is cut off accordingly for frequencies above 5 Hz. The process with the smallest time constant shown in Figure 7 (assigned to the SEI) is therefore not visible.
The disappearance of a polarization contribution of a process at moderate frequencies and 20% SOC can be confirmed. The number of processes in the range of moderate frequencies between 0.1 and 5 Hz corresponds to the number determined on the basis of frequency domain data. However, the greater spread of the time constants and the polarization with no discernible trend shows that due to the lower sampling frequency the accuracy is too low to reliably quantize the parameters of the identified processes at moderate frequencies.

Comparison of DRT by Time and Frequency Domain Data
The results of the experimental studies can be summarized as follows: • Processes with characteristic time constants that do not meet Equation (12) cannot be identified with the set sampling rate using the time domain data alone.

•
Already with a maximum sampling rate of 10 Hz, which is realistic for online applications, the DRT by time domain data can be used to identify charge transfer processes and to provide a qualitative description. However, the change in process parameters during aging can only be traced to a limited extent.

•
Quantitative statements for charge transfer processes, even if higher sampling rates are used, differ between the two methods due to the Butler-Volmer kinetics. In fact, due to the strongly non-linear excitations in real applications, the better transferability of the results obtained from frequency domain data is questionable.

•
The DRT based on time domain data is more sensitive for processes with large characteristic time constants such as solid state diffusion. • When using frequency domain data, either longer measuring periods are required or the measuring range is limited to higher frequencies, since different frequencies have to be excited successively during EIS. Figure 9a shows the capacity course during the aging of a cell and the corresponding polarization contribution of the process with the largest time constant, which is determined by time domain data at 60% SOC. The polarization contribution is referred to below as R tau max . The results belong to the cell that has already been evaluated in the previous For a better overview of the correlation between the two parameters, the development of R tau max is plotted over the capacity. A strong negative correlation can be found for ≥2250 cycles. During the first cycles, the increase in polarization in relation to the capacity fade is less pronounced. Regardless of which scenario, a significantly weaker increase in R tau max could be confirmed for all cells at the beginning of the aging tests. Thus, the average correlation coefficients C coe f f for each scenario and for the cycles ≥2250 are given in Table 3. Taking into account the measurement data only for cycles ≥2250, a Pearson correlation coefficient of ≥98% can be specified for scenarios Sc1 and Sc3, which indicates a high negative correlation between the two parameters. For scenario Sc2, the correlation is weaker, but still noticeable. Figure 10 shows the capacity courses over the cycles for all cells examined. For scenario Sc2, a recovery of cell capacity can be seen. After 2250 cycles, the capacity remains relatively constant at around 48 A h. Because of the flat curve, the accurate correlation is more difficult to determine because of the corresponding small changes in polarization R tau max are difficult to detect. This could be one reason why the correlation coefficient is lower in scenario Sc2.

Correlation of the Capacity and the Identified Process Parameters
To further investigate the reasons for the lower correlation coefficient, the correlation between the capacity of the anode and the polarization R tau max was determined. The determination of the loss of active material at the anode (LAM An ) on the basis of the non-invasive measurements during the checkups was carried out in [27]. The correlation between the LAM An and R tau max was ≥ 96% for each cell of scenario Sc1 from cycle 0 to 9050. Taking into account the uncertainties in the estimation of the LAM An , 96% is a high value and it can be stated that the increase in R tau max is associated with the LAM An . For all cells of scenario Sc1 the decrease in capacity at the anode was more pronounced compared to the cells of Sc2. Gantenbein et al. has already experimentally confirmed a higher LAM An for cells that were cycled at SOC levels above 65%, while the LAM An for cells that were cycled below 45% was neglible [31]. The lower LAM An provides a further reason for the lower correlation in Sc2 and the lower capacity fade during aging. For the cells of Sc3, the LAM An after 3850 cycles could not be determined due to a change in the electrode balance [27]. According to the results of [31] a comparatively high LAM An is to be expected. This statement is supported by the correlation factor determined for scenario Sc3 (see Table 3).
In contrast to the correlation to the cell capacity (see Figure 9), the linear relationship between LAM An and increase in R tau max is visible from the first few cycles of the longterm test. It is known from the literature that the cell capacity can be limited by different degradation mechanisms in different aging stages of a LIB (e.g., [35,36]). In a previous publication, the loss of lithium inventory (LLI) was identified as the main aging mechanism during the first few cycles of the aging study described [27]. The LAM An becomes the limiting factor for the capacity of the cells of scenarios Sc1 and Sc3 after 2250 cycles. Therefore, a correlation between the cell capacity and the R tau max only becomes visible after 2250 cycles. It is thus assumed that the correlation between the cell capacity and R tau max originates from the correlation of LAM An and R tau max .

Discussion
In the context of the aging study, the method introduced has proven to be advantageous for investigating processes with time constants in the range of minutes to hours. In addition to the more precise resolution of the low-frequency processes compared to the DRT, which is derived from frequency domain data, the time required for the measurements is significantly lower. Besides higher resolution and lower time efforts, the steady state criterion with regard to the SOC is not violated during the relaxation period. In order to avoid a change in the SOC during the EIS at low frequencies, the excitation current must be limited, which leads to a comparatively poor signal-to-noise ratio.
However, high sampling rates are required in order to investigate frequency-dependent processes with characteristic time constants in the range of a few seconds or even milliseconds such as charge transfer and SEI. According to the Shannon theorem, time constants in the range of milliseconds should theoretically be able to be evaluated with a sampling rate of 10 Hz, which is realistic for an online application. In contrast to this, the aging study has shown that only the parameters of processes with time constants of almost two decades higher (τ ≥ 1.6 s) can be quantified with sufficient accuracy using the aforementioned sampling rate. In many applications, and especially in most online applications, the method presented is therefore limited to the aging monitoring of processes with time constants in the range of seconds to several hours, such as solid state diffusion.
Therefore, for applications in which both pulse tests and EIS are possible, a combined use of both techniques is proposed. By merging the impedance data obtained through time and frequency domain data as suggested by [19], the impedance spectrum and DRT can be determined for a wide frequency range.
Using the DRT derived from the time domain data, a correlation between the polarization contribution and LAM An could be determined at low frequencies around 35 µHz. Due to the limited frequency range of the EIS, this connection could not be established using the frequency domain data. A correlation between a degradation mechanisms and the impedance of LIBs could already be established in earlier studies. In a comprehensive experimental study, the correlations between kinetic parameters of an LIB and degradation mechanisms could be derived [1]. Recently, a correlation between low frequency impedance and the LLI was found in [24] and Schindler et al. [37] confirmed a correlation between the low frequency impedance and LAM. Using cell kinetics to gain insights into the degradation mechanisms is beneficial because measuring impedances at high and moderate frequencies is fast compared to deriving the mechanisms by recording open circuit voltage curves [1]. However, EIS measurements suffer from long measurement durations for low-frequency impedances and a lack of availability in most state-of-the-art online applications. In contrast, relaxation periods occur in every online application. As a result, the introduced method enables a time-effective investigation that can be used online in order to obtain the relevant parameters of the processes with large time constants and thus enable an online estimate of the degradation mechanisms.
It is known that LIBs can suffer from accelerated capacity loss towards the end of life due to lithium plating originally caused by a lack of anodic material [35,36]. For online applications, an early detection of an increasing capacity reduction and lithium plating is of great importance for the reliability and safety of LIBs. Many publications like [38] have therefore dealt with the prediction of the sudden reduction in capacity. As a further outcome of this study, the correlation found could be used to identify an advanced LAM An online and thus the risk of lithium plating and the resulting rapid loss of capacity at an early stage.
In summary, this work showed that the evaluation of the voltage relaxation after a current pulse is suitable for investigating the kinetic parameters of a LIB even without prior knowledge of the polarization processes involved. The DRT obtained from the time domain data offers a time-effective determination of processes with large characteristic time constants and enables the parameters to be estimated online. The correlations found provide insights into the degradation mechanisms both offline and online and thus support the future battery design and enable optimization of the operating strategies. The findings about the early detection of a potential risk of lithium plating are available for subsequent studies and contribute to improving the reliability and safety of LIBs.
In a future step, Post-Mortem analyses will be carried out to directly quantify the LAM An and validate the assumptions made. In addition, the transferability of the results to other cell chemistries and formats must be checked.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: