A Method for Monitoring State-of-Charge of Lithium-Ion Cells Using Multi-Sine Signal Excitation

: In this paper, a method for monitoring SoC of a lithium-ion battery cell through continuous impedance measurement during cell operation is introduced. A multi-sine signal is applied to the cell operating current, and the cell SoH and SoC can be simultaneously monitored via impedance at each frequency. Unlike existing studies in which cell impedance measurement is performed ex situ through EIS equipment, cell state estimation is performed in situ. The measured impedance takes into account cell temperature and cell SoH, enabling accurate SoC estimation. The measurement system conﬁgured for the experiment and considerations for the selection of measurement parameters are described, and the accuracy of cell SoC estimation is presented. cannot be detected at higher frequencies, but the effects of fast kinetics can be revealed.


Introduction
A Li-ion battery is a type of secondary battery in which lithium ions move from a negative electrode to a positive electrode through an electrolyte during a discharge process. Lithium is the lightest solid element and has the lowest standard reduction potential. Therefore, when applied as an electrode material, the battery cell can obtain an electromotive force higher than 3 V, and a high energy density per weight and a high energy density per volume can be obtained. Li-ion batteries have less memory effect and less self-discharge as well. In addition to these advantages, due to their reliable performance, long life cycle, and advantages such as reduced pollution, Li-ion batteries are one of the most promising power sources for portable electronics, electric vehicles, renewable energy storage devices, et cetera [1,2]. However, high energy density batteries with low thermal stability electrode materials may have low safety performance [3][4][5]. In particular, battery failure caused by extreme conditions such as excessive external force, high temperature, low temperature, overcharge, and over-discharge becomes a serious problem [6][7][8].
Overcharging the Li-ion battery can be one of the most important safety issues. Overcharging occurs when a charging current is forced in after the battery reaches its upper cutoff voltage or state of charge (SoC) limit. It is usually caused by a malfunction of the battery charger or by an inaccurate estimate of the condition of the battery in the battery management system (BMS). When the Li-ion battery is overcharged, in addition to the increase in internal resistance, decomposition of the binder and electrolyte, formation of insoluble products, blocked electrode pores, and gas generation may occur sequentially [9]. Especially, severe expansion of the battery may occur due to gas accumulation [10], and when the internal pressure exceeds the limit, structural deformation, rupture, and an internal short circuit may occur [11,12]. Lithium metal and moisture in the air may react after the battery ruptures, and even flammable gases may ignite [8,13]. Heat generated by side reactions and internal short circuits can accelerate battery failure mechanisms through natural positive feedback, leading to thermal runaway, eventually causing battery ignition or explosion [8,14,15]. In the case of over-discharge, a significant loss of active lithium and positive electrode material occurs [16], and a decrease in battery capacity occurs [17,18].

•
Battery cells can be used longer by preventing overcharging and over-discharging, which can cause permanent damage.

•
It enables more aggressive cell operation. If the reliability of the SoC estimation is low, it must be operated conservatively in order to avoid overcharging and overdischarging the cell and to make the operation as safe as possible.

•
As it shows reliable estimation results for any usage profile of the cell, it improves the reliability of the use of the battery cell in the application system. • As the battery pack design does not have to be overengineered, it enables the production of smaller and lighter battery packs. It eventually lowers the price of the battery pack. In addition, the reliable battery system reduces battery maintenance costs.
Nonetheless, there is still no way to directly measure the SoC of a battery cell. Therefore, it must be estimated or inferred indirectly from the measured current, voltage, and temperature.

SoC estimation method based on open circuit voltage (OCV)
In order to represent the thermodynamic state of the energy of the cell electrode, the OCV according to the chemical composition of the electrode is used. OCV is the voltage of the electrode when it has been stable for a sufficient time without flowing current. Cell voltage is related to temperature and electrode particle surface concentration, whereas cell SoC is related to particle average concentration. In other words, OCV is measured when the electrodes of an electrochemical cell reach equilibrium and there is no voltage deviation depending on the position inside the electrode, which reflects the Gibbs free energy at thermodynamic equilibrium. The OCV represents a strong dependence on SoC in most batteries. Nevertheless, it is impractical for real-time or continuous state estimation since a relaxation time of several hours is generally required to reach electrochemical equilibrium. This is especially problematic for battery applications where resting time can never exist. Furthermore, Li-ion batteries have a flatter OCV compared to lead-acid batteries, making SoC estimation difficult. In addition, cell temperature and cell state of health (SoH) can also lead to SoC estimation errors based on OCV.
• SoC estimation method based on ampere counting In the SoC estimation method based on ampere counting, the accumulated amount of charging and discharging current relative to the discharging capacity is defined as ∆SoC and calculated with the set initial SoC. Since only cumulative current information is used for SoC estimation, it has the advantage of requiring relatively low performance for the hardware and software of the BMS. However, this method also has weaknesses in estimating cell SoC. It becomes a problem if the SoC of the battery cell is entirely dependent on the initial SoC. Since only ∆SoC is calculated, it is unavoidable to misestimate the SoC if the initial SOC setting is incorrect. In addition, it is impossible to know exactly the total capacity and coulombic efficiency of a cell, which should be approximated, and both of these approximations contribute to the error in cell SoC estimation. Self-discharge currents and leakage currents from electronic circuits for measuring cell performance increase errors as well. Furthermore, errors in voltage, current, and temperature measurements contribute to increasing the estimation error. These errors are more integrated and intensified as battery cell operating time increases. As a result, the uncertainty in the SOC estimate is exacerbated by the accumulated measurement errors. Therefore, this method can show reliable estimation results for a short period of time only if the initial conditions are well known, unless there is a feedback mechanism for error correction.

•
SoC estimation method based on heuristic data The heuristic based SoC estimation method is a method based on experimental data. Statistical rules or patterns found from data obtained through various cell charge and discharge experiments are used. These methods include fuzzy logic, neural networks, and support vector machines (SVM). Reliable estimation results can be obtained when the learning technique is implemented with a large amount of experimental data under different conditions. However, it takes a lot of time for the necessary experimental data to be properly secured.

•
SoC estimation method based on adaptive control An adaptive control based SoC estimation method such as a Kalman filter (KF) [23] or a sliding mode observer [24,25] compares the actual SoC measurement result with the estimated value and gradually reduces the difference according to the feedback principle. Although this method has a high estimation performance, it is relatively complicated to implement and thus has a disadvantage of high cost.

•
SoC estimation method based on equivalent circuit model The electrochemical impedance spectroscopy (EIS) method is a well-established technique for determining the dynamic behavior of electrochemical systems [26,27]. It is used to characterize battery impedance behavior over a wide frequency range [28,29]. Using EIS, the measured spectrum in an electrochemical system can be interpreted as an impedance spectrum of a lumped element model consisting of resistors, inductors, and capacitors. It can be used to implement dynamic simulation models [30]. The SoC [21,22,31] and SoH [28,32,33] of the battery cell can be estimated through the equivalent circuit model obtained through EIS. Moreover, the commonly used definition for battery end of life (EoL) is the predefined battery impedance increase at nominal conditions [34][35][36]. Nevertheless, EIS equipment is generally used in laboratories for general propose impedance measurements [3], which are not suitable for battery monitoring purposes. For EIS measurements, battery cells need to be detached from the operating load, and it is usually time consuming for impedance measurements in a wide frequency range. In addition, EIS equipment can be an excessive investment for battery monitoring systems, especially heavy and bulky for portable devices. Moreover, without algorithms for cell state estimation, EIS measurements alone cannot estimate cell SoH and SoC.
Unlike the method using EIS, the proposed method can be used to estimate the SoC of the cell during operation. There is a growing interest in the use of cell impedance to monitor the condition of batteries. The papers by Qahouq [37] and Waag et al. [38] deal with the impedance measurement of the operating cell but do not cover cell SoC and SoH estimation using the measured impedance. In the paper by Do et al. [39], the measured impedance is not adjusted according to the temperature change, and the computational complexity increases by using the extended KF. In addition, this paper does not deal with the method of estimating the cell SoH and SoC from the measured impedance. In the papers by Huang et al. [40] and Howey et al. [41], the cell SoC is estimated by measuring the cell impedance, but the effect of temperature on the cell impedance is not considered. The paper by Fleischer et al. [33] shows a good result of estimating the state of a cell in operation. However, a physical model and an equivalent circuit model are used, and so-called mutation algorithms are used. Since nonlinear differential equations and matrix operations are required, the complexity of the operation is increased. On the other hand, the proposed method is relatively simple but accurate. Using the cell temperature and impedance at two frequencies, SoH and SoC of the cells in operation can be estimated.

Proposed Cell State Estimation Method
Battery cell SoC can be estimated by continuously measured impedance during discharge. To estimate the cell state, a multi-sine signal with a small amplitude is applied to the cell operating current, and the cell impedance is measured through the amplitude of the voltage response. The multi-sine signal is the sum of two different frequencies, and each frequency is used to estimate the cell SoH and SoC.
If the current i expressed by Equation (1) flows through the battery cell, the cell voltage e can be expressed by Equation (2).
where I dc is direct current (DC) bias, ∆I f is the amplitude of the excited test frequency f, E dc is the offset voltage, ∆E f is the amplitude of the output voltage, and φ f is the phase difference. Dividing the voltage by the current as Equation (3) produces a complex impedance Z f .
The electrochemical impedance of batteries depends on frequency and characterized by its modulus |Z f | and phase angle e jφ . Another expression is given as the real and imaginary parts of the complex impedance.
EIS measurements usually use a single-sine signal in which individual frequencies are measured sequentially, which is also known as stepped sine or frequency sweep. Therefore, single-sine EIS has the disadvantage that it takes a long time to acquire impedance in a wide frequency range. This disadvantage can be overcome by measuring several frequencies simultaneously. The method of measuring multiple frequencies at the same time is called multi-sine EIS. Multi-sine signals have already been used for impedance spectroscopy and transfer function measurements in biomedical applications [42,43], material characterization [44], and other fields such as battery measurements [26,28]. Nonetheless, multi-sine EIS is not a common method for estimating the in situ state of a battery cell. In general, the multi-sine EIS, like the single-sine EIS, requires the cell to be separated from the application circuit for impedance measurements. In the proposed battery SoC monitoring method, the sum signal of two test frequencies is excited to the cell operating current and its response voltage is measured. A Fourier transform is used on the sampled cell voltage to obtain the amplitude at each test frequency. The impedance at each test frequency is obtained by substituting the amplitude of each response voltage into Equation (3).

Measurement System
The measurement system is configured to measure cell impedance by applying a multisine signal to the operating current. Table 1 shows the specifications of the Li-ion battery cell used and Figure 1 shows a block diagram and a picture of the measurement system.  The test battery cell is placed in the temperature chamber of the Binder GmbH set to 25 °C. NXP's silicon temperature sensor KTY 81-110 is attached to the cell to measure the actual cell surface temperature. The voltage signal output port of USB-6212, which is a data acquisition (DAQ) module of the Natural Instruments, is connected to an electronic load to apply test frequencies to the cell operating current. This DAQ module also acquires cell voltage, current, and temperature data. Experiments are controlled by adjusting measurement parameters through a graphical user interface (GUI) created using LabVIEW from National Instruments. Measured and calculated data is displayed on the monitor and can also be written to the hard disk. The analog inputs of the DAQ module have 16bit analogue to digital converter (ADC) resolution, a maximum sample rate of 400 kS/s, and an input range of ±10 V. The analog output has 16-bit digital to analog converter (DAC) resolution, an output range of ±10 V, and a maximum update rate of 250 kS/s. The operation of the electronic load is described in the author's previous paper [45].

Definition of Key Terms
There is no standard definition that defines SoC and SoH for battery cells and may be defined differently depending on the battery application system. Therefore, it is necessary to clarify the terms defined in this paper before discussing cell state estimation. •

Definition of C-rate
In general, C-rate is used to indicate the rate of charge and discharge of a battery cell.
The relationship shown in Equation (4) is given between the cell nominal capacity C and the operating current .
where h is the time (hour) taken to completely discharge the battery cell, is the current at which the cell operates, and C is the nominal capacity of the battery (Ah). Here, the reciprocal of the h value is defined as the C-rate. •

Definition of SoC
In this paper, the cell SoC is set through the following three steps: The test battery cell is placed in the temperature chamber of the Binder GmbH set to 25 • C. NXP's silicon temperature sensor KTY 81-110 is attached to the cell to measure the actual cell surface temperature. The voltage signal output port of USB-6212, which is a data acquisition (DAQ) module of the Natural Instruments, is connected to an electronic load to apply test frequencies to the cell operating current. This DAQ module also acquires cell voltage, current, and temperature data. Experiments are controlled by adjusting measurement parameters through a graphical user interface (GUI) created using LabVIEW from National Instruments. Measured and calculated data is displayed on the monitor and can also be written to the hard disk. The analog inputs of the DAQ module have 16-bit analogue to digital converter (ADC) resolution, a maximum sample rate of 400 kS/s, and an input range of ±10 V. The analog output has 16-bit digital to analog converter (DAC) resolution, an output range of ±10 V, and a maximum update rate of 250 kS/s. The operation of the electronic load is described in the author's previous paper [45].

Definition of Key Terms
There is no standard definition that defines SoC and SoH for battery cells and may be defined differently depending on the battery application system. Therefore, it is necessary to clarify the terms defined in this paper before discussing cell state estimation.

•
Definition of C-rate In general, C-rate is used to indicate the rate of charge and discharge of a battery cell. The relationship shown in Equation (4) is given between the cell nominal capacity C N and the operating current i.
where h is the time (hour) taken to completely discharge the battery cell, i is the current at which the cell operates, and C N is the nominal capacity of the battery (Ah). Here, the reciprocal of the h value is defined as the C-rate.

• Definition of SoC
In this paper, the cell SoC is set through the following three steps: (1) Full charge of the battery cell. The constant-current (CC) charging procedure at 1C and the constant-voltage (CV) charging procedure at 4.2 V are followed in sequence, where 4.2 V is the upper cut-off voltage specified by the cell manufacturer. This first step is complete when the charging current drops below C/10. Subsequently, each cell is given 90 min of relaxation time.
Batteries 2021, 7, 76 6 of 18 (2) Discharge the cell at 1C until it reaches a lower cut-off voltage. From this step, the actual cell capacity C cell is obtained.
(3) Set the cell to the target SoC. To fully charge the cell, Stage 1 is performed again, and then the charge of the cell is consumed until the target SoC is reached. Since different cell relaxation times can introduce errors, the same relaxation times must be observed [46,47]. Each cell is given 90 min of relaxation time after the target SoC is set and before it is used in the experiment.
The SoC is then defined as Equation (5).
The SoC of a cell is the ratio of the residual capacity (C residual ) to the total available capacity when the battery cell is fully charged (C cell ). The C cell can be obtained from the second step of the procedure for setting up the cell SoC mentioned above. Depth of discharge (DoD) has the exact opposite definition of SoC, i.e., a cell of 80% SoC has the same meaning as a cell of 20% DoD.

Definition of SoH
Battery cell SoH is defined as the ratio of the cell actual total available capacity C cell to the cell initial total available capacity C init and expressed by Equation (6).
In this paper, cell impedances from 100% to 80% SoH of cells are measured and compared. The EoL of a battery is reached when the energy content or power capacity is no longer sufficient for the application. This can vary depending on battery application; thus, there is no universal definition of how many cell SoHs have reached EoL. Nonetheless, publications assume that EoL is reached when the battery cell SoH is less than 80% [48-50].

Selection of Test Amplitude
Depending on the battery cell application, the amplitude of the test frequency should be selected taking into account the trade-off between measurement accuracy and investment costs for battery management. In order for the appropriate test amplitude to be selected for impedance measurements, it must be selected to be small enough not to violate the linear criteria of the electrochemical battery system but large enough to obtain a suitable signal-to-noise ratio for robust measurements. The amplitude of the test signal must be small enough to satisfy the pseudo-linearity of the cell response voltage. This is crucial because the current-voltage curve of a Li-ion cell shows a nonlinear relationship that follows the Butler-Volmer kinetics. An excitation amplitude of 5 to 10 mV is generally recommended [51], e.g., if it is used for CC measurement (galvanostatic), the excitation amplitude of the current must be set so that the maximum absolute value of the voltage amplitude does not exceed 10 mV. Especially when a multi-signal signal is used, the maximum amplitude of the applied signal should be considered. Because the amplitudes at each frequency increase the total amplitude of the test signal, pseudo-linearity can be violated. For more test frequencies to be used, the amplitude of each frequency must be lowered, i.e., if a multi-sine signal of N frequencies is used, the output voltage amplitude at each frequency must be reduced to a maximum 10/N·mV. Amplitudes that are too small can be a problem as well. This is due to the finite resolution of the hardware that acquires and generates the test signal. Moreover, impedance measurements are not possible when the output voltage amplitude at individual frequencies becomes less than the system noise.

Selection of Test Frequencies
In a Li-ion battery cell, charged particles containing lithium ions are always involved in the electrode process. It is generally believed that the time constants of the movement or reaction of charged particles in these processes are different. Thus, impedance at different frequencies is related to different processes [52]. The influence of mass transfer, whose time constant is slower, is dominant at impedance at lower frequencies. Mass transfer (including diffusion) becomes faster as the temperature is higher and as the concentration of Li-ions increases (as DoD is lower). Faster transfer speed can be interpreted as lower impedance. On the other hand, the lower the cell temperature or the higher the DoD, the slower the mass transport [53], which appears to be a higher impedance in the lower frequency. This means that impedance at lower frequency can be more advantageous as it is used for cell SoC estimation. Studies on battery cell impedance consistently show that impedance at lower frequencies is more dependent on changes in cell SoC. However, since the frequency is the reciprocal of the period, the lower the frequency, the longer the period required to measure each continuous impedance. For example, if 1 mHz is used for the measurement, it takes about 17 min to collect only one period of signal. Therefore, the test frequency should be selected considering the time interval required for state estimation depending on the battery application. There is one more problem: the longer the measurement of one period, the greater the change in the internal state of the cell such as SoC and cell internal temperature. This also increases the error in impedance measurements. In this paper, 1 Hz where one period of signal can be collected per second is used as an example for SoC estimation.
There are also some considerations for higher frequencies to be used. In this case, the achievable sampling rate in the hardware can be a constraint. In theory, satisfying the Nyquist-Shannon criteria is sufficient for the signal to be restored, but, generally, oversampling is required due to noise and non-ideal properties in analog filters. Another problem is that the cell impedance is generally lower at higher frequencies, which can lead to the same problems when the test amplitude is too low. Most cell degradation results in an increase in internal resistance [54][55][56][57], and the 1 kHz impedance of the cell used in this paper represents the ohmic resistance of the cell in which the imaginary part of the impedance becomes 0. In paper [46], the accuracy of SoH estimation with impedance at different frequencies is compared, which shows that impedance at higher frequencies is more advantageous for use in SoH estimation, which is less affected by cell SoC and temperature. In this paper, a frequency of 1 kHz is used for SoH estimation as an example.

Measurement Results
The measurement system and parameters described in Section 2 are used for cell impedance measurements. The overall process of cell state estimation in this paper is shown in Figure 2 as a flow chart. Table 2 shows the experimental conditions. Cell impedance is measured every second while discharging from 0% to 100% DoD for each cell with 100% to 80% SoH. A DC bias of 1C and an excitation amplitude of 130 mA are applied, i.e., an amplitude of 65 mA is given at each frequency of 1 Hz and 1 kHz. This amplitude is chosen for the cell output voltage amplitude to be less than 10 mV. The multi-sine test signal, the sum signal of the two frequencies, applies to the cell DC bias, and the impedance at each frequency is measured simultaneously.  Table 2 shows the experimental conditions. Cell impedance is measured every second while discharging from 0% to 100% DoD for each cell with 100% to 80% SoH. A DC bias of 1C and an excitation amplitude of 130 mA are applied, i.e., an amplitude of 65 mA is given at each frequency of 1 Hz and 1 kHz. This amplitude is chosen for the cell output voltage amplitude to be less than 10 mV. The multi-sine test signal, the sum signal of the two frequencies, applies to the cell DC bias, and the impedance at each frequency is measured simultaneously.     Figure 3 shows continuous 1 Hz and 1 kHz impedance during discharging of cells with different SoHs.  Table 2 shows the experimental conditions. Cell impedance is measured every second while discharging from 0% to 100% DoD for each cell with 100% to 80% SoH. A DC bias of 1C and an excitation amplitude of 130 mA are applied, i.e., an amplitude of 65 mA is given at each frequency of 1 Hz and 1 kHz. This amplitude is chosen for the cell output voltage amplitude to be less than 10 mV. The multi-sine test signal, the sum signal of the two frequencies, applies to the cell DC bias, and the impedance at each frequency is measured simultaneously.   As shown in Figure 3, the lower the cell SoH, the higher the cell impedance, and while 1 kHz impedance is relatively constant for DoD changes, 1 Hz impedance is more affected by cell DoD. While the cell is completely discharged, the increased rate of the highest value to the lowest value of 1 kHz impedance is 3.74%. At the same time, for 1 Hz impedance, the increased rate of the highest value to the lowest value between DoD 5% and 95% is 10.16%. The 1 kHz impedance shows a smaller increased rate, and it indicates less susceptibility to changes in cell DoD. Therefore, the 1 kHz impedance is more suitable to be used for SoH estimation of cells for which the cell SoC is unknown. Contrastively, 1 Hz impedance shows larger increase rate, which means that it is more suitable for cell SoC estimation.
For the calculation of the increase rate, 1 Hz impedance between DoD 5% and 95% in Figure 3 are used (if the full range of DoD 0% to 100% is considered, the difference in impedance deviation at 1 Hz and 1 kHz will be more prominent.). The reason is that at the beginning and end of cell discharge, the 1 Hz impedance is not measured correctly and cannot be used for cell SoC estimation. The voltage response of an electrochemical cell is not only affected by the cell impedance. The cell discharge process causes an overpotential, which is the voltage loss described by polarization. Polarization refers to a phenomenon in which the electrode potential becomes excessive or insufficient in an equilibrium state. During the reaction process in a battery cell, the rate of charge transfer process occurring in each cell component is not the same. If this rate is relatively slow for a particular process, it becomes the rate-limiting process for the entire reaction of the cell. As the cell discharges, current flows between both terminals, causing the voltage to be measured below its equilibrium potential. In this case, the difference between the voltage at both terminals and the equilibrium voltage is called overpotential, indicating the degree of polarization. Figure 4 shows the three polarization regions appearing in the typical discharge curve of a Li-ion battery cell. In the cell discharge process, polarization can be classified into three categories: activation polarization, concentration polarization, and ohmic polarization [58]. Figure 3, the lower the cell SoH, the higher the cell impedance, and while 1 kHz impedance is relatively constant for DoD changes, 1 Hz impedance is more affected by cell DoD. While the cell is completely discharged, the increased rate of the highest value to the lowest value of 1 kHz impedance is 3.74%. At the same time, for 1 Hz impedance, the increased rate of the highest value to the lowest value between DoD 5% and 95% is 10.16%. The 1 kHz impedance shows a smaller increased rate, and it indicates less susceptibility to changes in cell DoD. Therefore, the 1 kHz impedance is more suitable to be used for SoH estimation of cells for which the cell SoC is unknown. Contrastively, 1 Hz impedance shows larger increase rate, which means that it is more suitable for cell SoC estimation.

As shown in
For the calculation of the increase rate, 1 Hz impedance between DoD 5% and 95% in Figure 3 are used (if the full range of DoD 0% to 100% is considered, the difference in impedance deviation at 1 Hz and 1 kHz will be more prominent.). The reason is that at the beginning and end of cell discharge, the 1 Hz impedance is not measured correctly and cannot be used for cell SoC estimation. The voltage response of an electrochemical cell is not only affected by the cell impedance. The cell discharge process causes an overpotential, which is the voltage loss described by polarization. Polarization refers to a phenomenon in which the electrode potential becomes excessive or insufficient in an equilibrium state. During the reaction process in a battery cell, the rate of charge transfer process occurring in each cell component is not the same. If this rate is relatively slow for a particular process, it becomes the rate-limiting process for the entire reaction of the cell. As the cell discharges, current flows between both terminals, causing the voltage to be measured below its equilibrium potential. In this case, the difference between the voltage at both terminals and the equilibrium voltage is called overpotential, indicating the degree of polarization. Figure 4 shows the three polarization regions appearing in the typical discharge curve of a Li-ion battery cell. In the cell discharge process, polarization can be classified into three categories: activation polarization, concentration polarization, and ohmic polarization [58]. Ohmic polarization, also referred to as ohmic loss, is caused by the current flowing through the internal resistance of the battery cell. The greater the cell internal impedance, the greater the operating voltage drop as the ohmic losses increase proportional to the current density. Therefore, in this region, it has a slope related to the cell impedance and shows the most linearity when charging and discharging the cell. Due to this linearity, impedance can be measured accurately; hence, it is the most suitable region for cell state estimation through impedance measurement. Meanwhile, the cell output voltage in the Ohmic polarization, also referred to as ohmic loss, is caused by the current flowing through the internal resistance of the battery cell. The greater the cell internal impedance, the greater the operating voltage drop as the ohmic losses increase proportional to the current density. Therefore, in this region, it has a slope related to the cell impedance and shows the most linearity when charging and discharging the cell. Due to this linearity, impedance can be measured accurately; hence, it is the most suitable region for cell state estimation through impedance measurement. Meanwhile, the cell output voltage in the activation and concentration polarization regions has nonlinearity. Activation polarization is due to various delay factors inherent in the dynamics of electrochemical reactions, such as the work function that ions must overcome at the junction between the electrode and the electrolyte. It has a dominant effect at the beginning of cell discharge. Concentration polarization takes into account the resistance that ions face in the process of mass transfer (e.g., diffusion) as they move across the electrolyte and from one electrode to another. This polarization has a dominant effect at high cell DoD. In conclusion, the cell voltage drops significantly non-linearly as the cell discharges in the region of activation and concentration polarization. Because this nonlinearity is based on slow response, especially in battery cells, the lower the frequency, the more difficult it is to obtain an accurate cell impedance. As shown in Figure 3, the effect of activation and concentration polarization at a frequency of 1 Hz is more pronounced than at 1 kHz, because slow transport processes are well represented at lower frequencies. On the other hand, slow transport reactions cannot be detected at higher frequencies, but the effects of fast kinetics can be revealed.

Cell SoH Estimation
The measured 1 kHz impedance in Figure 3 is used to estimate the cell SoH. Figure 5 shows the average of the 1 kHz impedance at each SoH and the fitted line. and concentration polarization. Because this nonlinearity is based on slow response, especially in battery cells, the lower the frequency, the more difficult it is to obtain an accurate cell impedance. As shown in Figure 3, the effect of activation and concentration polarization at a frequency of 1 Hz is more pronounced than at 1 kHz, because slow transport processes are well represented at lower frequencies. On the other hand, slow transport reactions cannot be detected at higher frequencies, but the effects of fast kinetics can be revealed.

Cell SoH Estimation
The measured 1 kHz impedance in Figure 3 is used to estimate the cell SoH. Figure 5 shows the average of the 1 kHz impedance at each SoH and the fitted line. The fitted line in Figure 5 is obtained through linear regression of these 1 kHz impedances measured at each SoH, and it is expressed as Equation (7).
The measured 1 kHz impedance can be used to estimate the cell SoH using Equation (7). Here, root mean squared error (RMSE) is 0.21 mΩ and R is 0.99. R represents the goodness of fit of the estimated model. The closer this value is to 1, the better the estimate matches the target value.

Cell SoC Estimation
The cell 1 Hz impedance is greatly influenced by not only the cell SoC but also the cell SoH and cell temperature. Therefore, to be used for cell SoC estimation, the 1 Hz impedance is adjusted as cell temperature in Section 3.2.1 and normalized as cell SoH in Section 3.2.2.

Consideration of Cell Temperature
The change in temperature of a battery cell mainly depends on the C-rate, discharge time, relaxation time, and ambient temperature [59]. Figure 6 shows the cell temperature that rises when the cell is discharged at 1C current from DoD 0% to 100% in a temperature chamber set to 25 °C (same experimental conditions as shown in Table 2). The fitted line in Figure 5 is obtained through linear regression of these 1 kHz impedances measured at each SoH, and it is expressed as Equation (7).
The measured 1 kHz impedance can be used to estimate the cell SoH using Equation (7). Here, root mean squared error (RMSE) is 0.21 mΩ and R 2 is 0.99. R 2 represents the goodness of fit of the estimated model. The closer this value is to 1, the better the estimate matches the target value.

Cell SoC Estimation
The cell 1 Hz impedance is greatly influenced by not only the cell SoC but also the cell SoH and cell temperature. Therefore, to be used for cell SoC estimation, the 1 Hz impedance is adjusted as cell temperature in Section 3.2.1 and normalized as cell SoH in Section 3.2.2.

Consideration of Cell Temperature
The change in temperature of a battery cell mainly depends on the C-rate, discharge time, relaxation time, and ambient temperature [59]. Figure 6 shows the cell temperature that rises when the cell is discharged at 1C current from DoD 0% to 100% in a temperature chamber set to 25 • C (same experimental conditions as shown in Table 2). To measure the cell impedance at each constant temperature, cell impedances are measured through an EIS instrument (IM6ex from Zahner-Elektrik GmbH & CoKG) in a temperature chamber set to 20, 25, 30, and 35 °C. The measurement conditions are shown in Table 3.  To measure the cell impedance at each constant temperature, cell impedances are measured through an EIS instrument (IM6ex from Zahner-Elektrik GmbH & CoKG) in a temperature chamber set to 20, 25, 30, and 35 • C. The measurement conditions are shown in Table 3. The measurement results are shown in Figure 7 as a Nyquist plot and a Bode plot. To measure the cell impedance at each constant temperature, cell impedances are measured through an EIS instrument (IM6ex from Zahner-Elektrik GmbH & CoKG) in a temperature chamber set to 20, 25, 30, and 35 °C. The measurement conditions are shown in Table 3. The measurement results are shown in Figure 7 as a Nyquist plot and a Bode plot. As shown in Figure 7, the higher the cell temperature, the lower the cell impedance, and the 1 Hz impedance is more affected by temperature than the 1 kHz impedance.  As shown in Figure 7, the higher the cell temperature, the lower the cell impedance, and the 1 Hz impedance is more affected by temperature than the 1 kHz impedance. Figure 8 shows the 1 Hz impedance measured at each temperature and a curve fitted with a quadratic equation.
where T is the cell temperature, and Z(T) is the impedance obtained through the equation at the cell temperature T.
The measured 1 Hz impedance at various cell temperatures is uniformly adjusted to the cell impedance at 25 °C through Equation (9).
where Z(T) is the adjusted impedance to 25 °C and Z(T) is the measured impedance at temperature T.
The result of Equation (9) applied to the 1 Hz impedance in Figure 3 is shown in Figure 9.  Figure 8, where R 2 is 1.00 and RMSE is 0.19 mΩ.

Equation (8) is the quadratic equation fitted in
where T is the cell temperature, and Z(T) fitted is the impedance obtained through the equation at the cell temperature T. The measured 1 Hz impedance at various cell temperatures is uniformly adjusted to the cell impedance at 25 • C through Equation (9).
where Z(T) adjusted is the adjusted impedance to 25 • C and Z(T) measured is the measured impedance at temperature T.
The result of Equation (9) applied to the 1 Hz impedance in Figure 3 is shown in Figure 9.
where T is the cell temperature, and Z(T) is the impedance obtained through the equation at the cell temperature T.
The measured 1 Hz impedance at various cell temperatures is uniformly adjusted to the cell impedance at 25 °C through Equation (9).
where Z(T) is the adjusted impedance to 25 °C and Z(T) is the measured impedance at temperature T.
The result of Equation (9) applied to the 1 Hz impedance in Figure 3 is shown in Figure 9. In Figure 9, the 1 Hz impedance increases with increasing cell DoD. The higher the cell DoD, the higher the cell temperature ( Figure 6), and as the cell temperature increases, the cell impedance decreases (Figure 8). Figure 9 shows the impedance when the impedance of Figure 3 is adjusted to a lower temperature of 25 °C. Therefore, the temperatureadjusted 1 Hz impedance in Figure 9 becomes higher than the 1 Hz impedance in Figure  3.

Consideration of Cell SoH
As mentioned above, cell 1 Hz impedance is affected by cell SoH as well as cell temperature. The cell SoH can be estimated via the 1 kHz impedance even if the cell SoC is unknown (Section 3.1). The cell SoH obtained in Equation (7)  In Figure 9, the 1 Hz impedance increases with increasing cell DoD. The higher the cell DoD, the higher the cell temperature ( Figure 6), and as the cell temperature increases, the cell impedance decreases (Figure 8). Figure 9 shows the impedance when the impedance of Figure 3 is adjusted to a lower temperature of 25 • C. Therefore, the temperature-adjusted 1 Hz impedance in Figure 9 becomes higher than the 1 Hz impedance in Figure 3.

Consideration of Cell SoH
As mentioned above, cell 1 Hz impedance is affected by cell SoH as well as cell temperature. The cell SoH can be estimated via the 1 kHz impedance even if the cell SoC is unknown (Section 3.1). The cell SoH obtained in Equation (7) is used for normalizing the 1 Hz impedance. The measured cell 1 Hz impedance (Z measured ) can be normalized to Z norm by Equation (10).
In Equation (10), Z max and Z min represent the maximum and minimum values of the impedance between DoD 10% and 90% and can be obtained by Equations (11) and (12), respectively. In Equation (11), R 2 is 1.00 and RMSE is 0.12 mΩ, and in Equation (12), R 2 is 1.00 and RMSE is 0.26 mΩ. Z max = −0.00593·SoH 2 + 0.79922 SoH + 63.09 (11) Z min = −0.2969·SoH + 103.1 (12) The result of applying Equation (10) to the 1 Hz impedance in Figure 9 is shown in Figure 10. In Equation (10), Z and Z represent the maximum and minimum values of the impedance between DoD 10% and 90% and can be obtained by Equations (11) and (12), respectively. In Equation (11), R is 1.00 and RMSE is 0.12 mΩ, and in Equation (12) The result of applying Equation (10) to the 1 Hz impedance in Figure 9 is shown in Figure 10. The average of the normalized impedance at each SoH is fitted with a cubic equation and plotted as a curve in Figure 10. This fitted curve is expressed by Equation (13) and has an R of 1.00 and an RMSE of 1.51. Table 4 shows each R and RMSE when Equation (13) is applied to cells with differ- The average of the normalized impedance at each SoH is fitted with a cubic equation and plotted as a curve in Figure 10. This fitted curve is expressed by Equation (13) and has an R 2 of 1.00 and an RMSE of 1.51. Table 4 shows each R 2 and RMSE when Equation (13) is applied to cells with different SoHs. Because of the cell nonlinear response mentioned above, SoC estimation using 1 Hz impedance is made here between DoD 10% and 90%.

Battery Cell SoC Estimation at Different Initial SoCs
Section 3.2.2 shows the results of SoC estimation when a fully charged cell is fully discharged. However, battery cells are not always operated in a fully charged state. Here, battery cell SoC estimation at different cell initial SoCs is emulated. A total of 10 min of cell discharge and 60 min of cell relaxation are repeated until the cell is completely discharged. The cell impedance is measured every second while the cell is discharging, and each cell state after the relaxations represents the cell state at different initial SoCs. The experimental conditions are shown in Table 5.  Figure 11a shows the cell voltage and cell temperature for the cell DoD, and Figure 11b shows the 1 Hz and 1 kHz impedances.
As shown in Figure 11a, the cell voltage drops while the cell is discharging and increases during each relaxation. At the same time, the cell temperature increases while the cell is discharged and decreases during the relaxation time. Figure 11b shows that cell 1 Hz impedance is more affected by DoD than 1 kHz impedance. The cell 1 kHz impedance is relatively constant for DoD changes and can be used for cell SoH estimation as shown in Section 3.1, while the 1 Hz impedance is highly influenced by DoD, and it is used for SoC estimation. As shown in Figure 11a, the cell voltage drops while the cell is discharging and increases during each relaxation. At the same time, the cell temperature increases while the cell is discharged and decreases during the relaxation time. Figure 11b shows that cell 1 Hz impedance is more affected by DoD than 1 kHz impedance. The cell 1 kHz impedance is relatively constant for DoD changes and can be used for cell SoH estimation as shown in Section 3.1, while the 1 Hz impedance is highly influenced by DoD, and it is used for SoC estimation.
In Figure 11b, the DoD ranges where the 1 Hz impedance cannot be measured correctly are shown in gray. As explained earlier in Chapter 3, the cell 1 Hz impedance is not measured correctly at the beginning of each discharge, especially due to the nonlinearity of the cell discharge curve due to activation polarization shown in Figure 4. Each of these gray ranges in Figure 11b corresponds to 265 s. In addition, the range from DoD 95% in Figure 11b is also included in these gray ranges, as this range corresponds to the concentration polarization region in Figure 4. However, it should be noted that in most cell operations where cells are rarely fully discharged, a range that cell impedance cannot be In Figure 11b, the DoD ranges where the 1 Hz impedance cannot be measured correctly are shown in gray. As explained earlier in Chapter 3, the cell 1 Hz impedance is not measured correctly at the beginning of each discharge, especially due to the nonlinearity of the cell discharge curve due to activation polarization shown in Figure 4. Each of these gray ranges in Figure 11b corresponds to 265 s. In addition, the range from DoD 95% in Figure 11b is also included in these gray ranges, as this range corresponds to the concentration polarization region in Figure 4. However, it should be noted that in most cell operations where cells are rarely fully discharged, a range that cell impedance cannot be measured correctly only appear once in the beginning of cell discharge. In Figure 11b, multiple gray ranges are shown because the impedance at the beginning of the discharge in different initial SoCs is shown in a single figure.
The 1 Hz impedance of Figure 11b adjusted to the cell temperature is shown in Figure 12a. Equation (9) is used to take into account the effects of cell temperature. The temperature-adjusted 1 Hz impedance is then normalized to the cell SoH, and the result is shown as a dotted line in Figure 12b. Equation (10) is used to normalize the temperatureadjusted impedance to the cell SoH.
Batteries 2021, 7, x FOR PEER REVIEW 15 of 18 measured correctly only appear once in the beginning of cell discharge. In Figure 11b, multiple gray ranges are shown because the impedance at the beginning of the discharge in different initial SoCs is shown in a single figure.
The 1 Hz impedance of Figure 11b adjusted to the cell temperature is shown in Figure  12a. Equation (9) is used to take into account the effects of cell temperature. The temperature-adjusted 1 Hz impedance is then normalized to the cell SoH, and the result is shown as a dotted line in Figure 12b. Equation (10) is used to normalize the temperature-adjusted impedance to the cell SoH. The cell SoC can be estimated by applying Equation (13) to Figure 12b, where R is 0.99 and RMSE is 3.07 mΩ. Cell impedance in grayed-out ranges is excluded from SoC estimation. These ranges include a range of 265 s after the start of the discharge and a range of cell DoD exceeding 95%. The purpose of this experiment is to show that SoC estimation is possible even when discharges are initiated from different cell SoCs, so it can lead to misunderstanding that there are too many grayed-out ranges where SoC estimation is impossible. Again, in most cases when battery cells are used, this grayed-out range The cell SoC can be estimated by applying Equation (13) to Figure 12b, where R 2 is 0.99 and RMSE is 3.07 mΩ. Cell impedance in grayed-out ranges is excluded from SoC estimation. These ranges include a range of 265 s after the start of the discharge and a range of cell DoD exceeding 95%. The purpose of this experiment is to show that SoC estimation is possible even when discharges are initiated from different cell SoCs, so it can lead to misunderstanding that there are too many grayed-out ranges where SoC estimation is impossible. Again, in most cases when battery cells are used, this grayed-out range appears only once a battery cell begins to discharge and rarely twice when the cell is completely discharged.

Conclusions and Discussion
This paper introduces a SoC monitoring method of Li-ion battery cells using impedance measurement. For accurate SoC estimation, the estimated cell SoH is considered with the measured temperature. Unlike traditional EIS measurement methods, the proposed method does not require impedance at wide frequencies, saving measurement time and simplifying measurement saves hardware costs. Especially, a multi-sine signal is applied to measure the cell impedance at two frequencies simultaneously. The cell impedance at 1 kHz for SoH estimation and the cell impedance at 1 Hz for SoC estimation are used. As a result, the proposed cell SoC monitoring method enables simultaneous estimation of unknown cell SoH and SoC. This is verified through an experiment in several different initial SoCs in a cell as well.
One problem to be pointed out is that at the start of battery cell discharge (ca. 265 s at 1C) and in the high DoD range (over ca. 95%), accurate cell impedance measurement is not possible due to the large nonlinearity in the cell voltage response; hence, the SoC cannot be correctly estimated. Nonetheless, the proposed SoC estimation method using impedance can be used together with the SoC estimation method using cell voltage without any additional hardware or measurement. This is because the proposed method already measures the cell voltage response to obtain cell impedance. At the same time, this SoC estimation method using cell impedance can compensate for the weaknesses of the existing SoC estimation methods using cell voltage. The voltage of the Li-ion battery cell drops significantly at the beginning and end of the discharge. This enables SoC estimation as a simple method through cell voltage measurement. However, the decrease in voltage in the middle region of the cell discharge curve, the so-called "flat plateau", is not significantly noticeable, and this is one of the factors that makes it difficult to estimate SoC simply with the measured cell voltage. In particular, this flat plateau is notorious for estimating the SoC of lithium iron phosphate (LFP) battery cells. Contrariwise, this flat plateau makes impedance measurements more accurate. This is because accurate impedance measurement is possible when the target system is linear. That is, when the proposed method is used, SoC and SoH estimation is possible on this flat plateau and even more accurate.
In this paper, only the experimental results of lithium nickel manganese cobalt oxide (NMC) type cells are shown. The proposed method is based on EIS, which has already been validated in other types of battery cells through numerous literatures. Thus, it can be expected to be applied to other types of Li-ion cells. Nevertheless, it is worth comparing the results of applying this method to other types of cells in future studies. This paper only deals with SoC estimation during battery cell discharge. Accurate SoC estimation during cell charging is also worth studying. However, there is the problem that this method is only possible during CC charging. This is because the offset current changes frequently during CV charging, making it impossible to measure the impedance correctly. If SoC estimation is possible only during CC charging, the range of SoC that can be estimated is limited. The authors evaluate that monitoring SoC under this limited condition does not have a definite advantage; hence, this is not covered in this paper.
Lastly, the experiments in this paper consider the estimation of the state of a cell discharged at room temperature. Therefore, the experiment is conducted between 20 and 35 • C. However, since the cell state estimation method in this paper has the advantage of considering the effect of cell temperature on impedance, experiments will be conducted in a wider temperature range for wider application of this method in the future.