Practical identifiability of electrochemical P2D models for lithium-ion batteries

Electrochemical models play a significant role in today’s rapid development and enhancement of lithium-ion batteries. For instance, they are applied for design and process optimization. More recently, model and parameter identifiability are gaining interest as thorough model parameterization is key to reliable simulation results. Especially electrochemical models are often prone to unidentifiability and overfitting due to their high number of adjustable parameters. In this article, the most common electrochemical peudo-2D model of a lithium-ion battery is parameterized. A three-step procedure is applied which considers quasi-static 3-electrode measurements of the open-circuit potential, C-rate tests, and electrochemical impedance spectra. Identifiability of each step is discussed in-depth and a general guidance for future parameterizations is derived. The conducted study reveals the insufficiency of open-circuit potential and C-rate tests to fully parameterize the electrochemical model. Highly dynamic tests, e.g., impedance spectroscopy, are required to resolve the ambiguity of diffusive and electric processes under quasi-static conditions. Any parameterization of electrochemical models requires experimental data of electrode-resolved tests, as well as a combination of quasi-static and highly dynamic tests. The results of this study provide guidance for the use of electrochemical models in applied sciences and industry.


Introduction
Nowadays, the lithium-ion battery (LIB) is the dominant electrochemical storage technology for a large variety of applications. The further development of LIBs is based on enhancements of battery materials and the extended application of novel modeling approaches used for optimization.
Due to new applications for models and a significant increase of computational power, a wide variety of models has been developed and different applications have been addressed, e.g., models have been used to quantify solid-electrolyte interface (SEI) formation and aging [14,34,39]. The variety of models is summarized in reviews about modeling of LIBs with focus on systems engineering, multi-scale modeling, and state estimation in electric cars, respectively [12,19,28,33].
Besides the governing equations, model parameters are an essential part of every model. Their correct choice can be not only crucial but also challenging. While parameter estimation (PE) is commonly used to parameterize models, it is also a powerful tool to derive non-measurable parameters from cell performance. These yield insight, e.g., into the impact of production parameters [22]. Both applications require practical identifiability of the model parameters. In the following, parameter estimation approaches are briefly reviewed, and direct and sample-based approaches of proofing identifiability are summarized.
The most common parameters in electrochemical LIB models are reaction rate constants, diffusion coefficients of active materials and electrolyte, and conductivity of solid and electrolyte. Temperature dependencies are commonly modeled using Arrhenius' law. It has a pre-exponential coefficient and an activation energy which are estimated from temperature variations. Concentration dependencies can also be considered and may further increase the number of parameters. For instance, in Ref. [11], the state of charge (SOC) or concentration-dependent open-circuit potential (OCP) between both electrodes was estimated from electrochemical experiments alongside the other model parameters. Fitting empirical OCP equations for both electrodes together with kinetic constants. etc., led to a large number of parameters. To deconvolute the joint parameter estimation of OCP curves and other model parameters Refs. [20] and [22] first parameterized the thermodynamic, OCP-related parameters. Then, in the subsequent main parameter estimation kinetic constants, diffusion coefficients, etc., were estimated. For equivalent circuit models, resistances, capacitance, and OCP have to be estimated, and concentration dependencies are common [19]. The number of estimated parameters varies significantly between different publications. This is exemplary illustrated in Table 1. The difference in numbers is partially related to different model complexities and types, and partially to the application of different parameter estimation procedures. Further, the number of uncertain parameters is dependent on the model application, e.g., whether temperatures are varied in the experiments and its influence is modeled.
As Table 1 shows, the most common case is a P2D model without temperature dependency in combination with C-rate tests. This commonly leads to three to 14 parameters. The difference in total number of parameters is partially related to differences in parameterization approaches and research objectives. But this differences also raises the equation about the maximum number of independent uncertain parameters which can be estimated in one model. Different parameter procedures were introduced in the literature. For physics-based full-order models (P2D models and SPMs in Table 1), Ramadesigan et al. applied gradientbased, thus local, least-squares fitting to estimate parameter changes due to cell aging in combination with an uncertainty quantification to calculate the confidence intervals of the estimated parameters [32]. Vazquez-Arenas et al. coupled the parameter estimation with a sensitivity analysis to assess the estimation accuracy [43]. Lenze et al. estimated parameters of an electrochemical model manually adjusting parameters [22]. Chun and Han introduced a cascaded improved harmony search to estimate the parameters of the P2D model, wherein a partially random-based optimization approach was conducted repeatedly [4].
In the engineering and control community, there exists a variety of reduced-order models or equivalent circuit models (ECMs) [1,18,29,45]. For ECMs, parameter estimation is frequently conducted [1]. These models commonly focus on estimation of battery states like SOC, state of power, or state of health [16,25,29,42,44]. Dvorak et al. introduced a parameter estimation method for an ECM containing 2 RC elements. This method applies step-wise discharge at constant current with repetitions at different C-rates and temperatures [7]. A Kalman filter was applied to the ECM. It was shown that in this case, parameter estimation should facilitate electrochemical tests similar to the actual application for optimal accordance of the model [2].
Besides the choice of the experimental-electrical test procedure and optimization algorithms, in some publications, the objective function was adjusted to enhance identifiability. Different objective functions were applied to different parts of the experimental data. For instance, Li et al. applied three steps: high SOC and low current, low SOC and low current, and eventually considered a dynamic driving cycle [23]. Jin et al. considered thermodynamic parameters and kinetic parameters independently for a P2D model as well as linear or logarithmic scaling for different parameters [17].
All these publications have a strong theoretical focus. In contrast, Ecker et al. used a variety of different experimental analysis methods to measure model parameters like diffusion coefficients. They used electrochemical impedance spectroscopy in combination with equivalent circuit models to partially estimate the parameters of the P2D model [8,9]. Schmalstieg et al. presented a method to parameterize a model of bought-in commercial cells. Their method combined EIS, galvanostatic intermittent titration technique, porosity measurement by mercury intrusion, and ECM fitting [36].
Identifiability is basically a model property which describes whether its parameters can be identified from its output for a specified input [3]. Further, the terms structural identifiability [3], practical identifiability [27], and uniqueness are introduced. Structural identifiability focuses on the model function itself and assesses whether identifiability is given or not. In general, it refers to global identifiability. Practical identifiability also considers the sensitivity of the model to the parameters. Practical identifiability requires sensitivity and is related to local identifiability. Uniqueness of the estimated parameter set yields global identifiability.
Sharma and Fathi [38] and Lin and Stefanopoulou [24] used the Fisher information matrix to assess local identifiability. From the Fisher information matrix, the eigenvalues can be used to assess identifiability [10,37], and a singular value decomposition can be carried out. Bizeray et al. linearized a single-particle model to investigate its structural identifiability [3]. This approach has the advantage to be able to assess global identifiability and hence uniqueness but requires a closed formulation of the models output and linearization. Thus, it is not applicable to the commonly used P2D battery model which lacks the closed formulation of its output. Also, it is possible that linearization reduces the model's identifiability. In contrast, the Fisher information matrix could be applied to a P2D model, but computational cost would be enormous for a state-of-the-art desktop PC. For a SPM, the Fisher information matrix was successfully applied by Pozzi et al. to derive optimal experiments for parameterization [31].
Barcellona et al. published a review on parameter identification techniques for lithium-ion battery models [1]. Therein, they were classified in online methods, offline methods, or analytically/numerical calculation methods. Online methods were commonly applied to ECMs for state estimation [1]. Tian et al. estimated parameters of an ECM from constant-current discharge data and assessed local identifiability analyzing the rank of the sensitivity matrix [41]. Pozna et al. showed that even a first-order RC equivalent circuit battery model is not identifiable [30]. For fractional-order models (FOM) of lithium-ion batteries, parameter identifiability was assessed by Guo et al. and Li et al. [13,23]. FOM is based on Laplace transformation and provides lumped parameters which have to be estimated.
The review of relevant literature shows the variety of parameter estimation approaches in different scientific communities. The intention of this article is to interconnect those groups and to provide hands-on guidance for parameter estimation of electrochemical models, especially for applied sciences. The advantage of the introduced approach is that it enables estimation of a unique set of model parameters of a P2D model. All existing methods were either global or applicable to the P2D model. In contrast to the published literature of P2D models, it allows assessment of practical identifiability and uniqueness. Thus, it provides the most reliable base for any further physical or electrochemical model-based analysis.
The outline of this article is as follows: In Sect. 2, the model as well as the parameterization strategy is introduced. In Sect. 3, the model is parameterized in three steps. The second step is the common combination of P2D model with C-rate tests. Its practical identifiability is discussed in-depth and a mitigation strategy for its partial unidentifiability is introduced in Sect. 3.3. This article ends with conclusions which provide guidance for the use of electrochemical models in applied sciences and industry.

Electrochemical model
In this article, the common P2D model is applied to simulate C-rate tests. The model was first introduced by Doyle et al. [5,6]. In the following, the governing equations are briefly reviewed. For a detailed set of equations and boundary conditions, we refer to our previous publication [20]. The model is a set of partial differential and algebraic equations. All spatial derivatives are discretized applying finite volume method.
Intercalation reactions at both electrodes are described by Butler-Volmer kinetics: with a concentration-dependent exchange current density i 0 In Eqs. 1 and 2, j Li is the intercalation current density, a s the volume-specific active area, the symmetry coefficient of the reaction, is the electrochemical overpotential, k the reaction rate constant, c e Li + concentration in electrolyte, c s the Li concentration in active material, and c max is the maximal Li concentration in active material; thus, c max − c s are available lattice vacancies for Li in the active material. Diffusion of lithium in spherical particles is described by Fick's law in a radial coordinate r: wherein D s is the solid-phase diffusion coefficient. Further, initial values c a,0 and c c,0 are introduced for the solid-phase volume elements of anode and cathode at t = 0 at a cell voltage of 4.2 V , respectively. In the electrolyte, diffusion and migration are considered and the liquid-phase potential e is governed by In Eq. 4, e,eff is the effective electrolyte conductivity, e is potential, x is a linear coordinate orthogonal to the plain electrode area, and t p is the transference number. The effective diffusion coefficient is derived from porosity and tortuosity : Further, the model considers electrochemical double layers at both electrodes which was first introduced to this model type by Legrand et al. [21]: with double-layer capacitance C DL and At last, the Nernst-Einstein equation is applied to link electrolyte diffusivity D, a frequently adjusted model parameter, and electrolyte conductivity e to reduce the number of independent adjustable parameters: It assumes a linear concentration dependency of e and is, thus, a simplification of the model. However, its influence is assumed to be small at low and moderate concentrations. The open-circuit potential curves are estimated by empirical equations 10 and 11. The equations are motivated by the work of Smith and Wang [40] but are adjusted to reproduce the characteristics of the given materials more precisely.
The normalized concentration c is given by Parameters of Eqs. 10 and 11 are listed in Table 7 in the Online Appendix.
In the end, cell voltage can be derived from the solidphase potential different between the two active material-tocurrent collector interfaces at x = 0 and x = L: This potential difference includes all potential losses in solid and liquid phases, as well as the OCPs and electrochemical overpotentials of both electrodes.
In total, there are 14 uncertain parameters. For these, the parameter estimation procedure is introduced in the following section. Temperature effects are not considered.
Due to computational efficiency, a single-particle model is used to simulate EIS. Where applicable, the same parameters are used as in the P2D model. For the SEI, there are additional parameters which are taken from Ref. [14] as the SEI is beyond the scope of this study. The model simulated diffusion solely in a single particle per electrode.
The model considers an additional adsorption kinetic at the SEI surface: wherein a vacant SEI lattice space is v and a lithium-filled vacancy Li . Further diffusion and charge transport in the SEI are considered: The linear coordinate in the SEI is denoted . However, for details about the single-particle model, it is referred to Ref. [14], as it is only a tool to derive parameters from impedance spectra. It could be substituted by other tools such as equivalent circuit models.
The model was implemented in Matlab R2017b. Multistart fitting was conducted at a high-performance cluster due to its parallelization feasibility. To derive the impedance from simulation data, Matlab built-in fast Fourier transformation was used. For all least-square problems, a Matlab built-in Levenberg-Marquardt algorithm was used.

Parameterization procedure
Based on the literature reviewed above, a three-step parameter estimation procedure is applied in this article. For instance, see Ref. [23] for sequential parameter estimation procedures. The experiments are an OCP measurement for PE step 1, C-rate tests for kinetic parameters in PE step 2, and EIS data for improved identification of kinetic parameters in PE step 3. All experiments are conducted in a threeelectrode setup. The different steps of the applied parameter estimation procedure are summarized in Table 2. This approach uses the higher sensitivity of various experiments to certain parameters to gain more reliable results.
PE step 1 identifies static parameters, like the initial concentration of both electrodes, c a,0 and c c,0 , and the specific capacity of both active materials, c a,max and c c,max which are not affected by kinetics. PE step 2 is denoted as quasistatic. It is related to kinetic parameters, which are only sensitive at non-zero cell currents. This includes diffusion coefficients of both active materials, respectively, D s,a and D s,c , exchange current densities, i 0,a and i 0,c , electric conductivities, s,q and s,c , and electrode tortuosities, a and c , effecting, e.g., the ionic conductivity of the electrolyte. These parameters affect the cell performance at a time scale of seconds to hours.  The electrochemical cell response at small time scales of seconds and below is considered in PE step 3. It concerns especially reaction kinetics and electric transport. The dynamic cell behavior is analyzed, e.g., via impedance spectroscopy, and the double-layer capacitances, C DL,a and C DL,c , are estimated. Exchange current densities, i 0,a and i 0,c , and solid-phase conductivities, s,c and s,a , are recalculated, wherein the result of PE step 2 is used as starting value. Solid-phase diffusion coefficients are not recalculated in this step, as they can be estimated precisely in PE step 2. The choice of parameters is based on the literature review in Sect. 1 and especially on the sensitivity analysis conducted in our previous work [20].
In the following, the least-square formulations for all three steps are introduced. As the measured charges and voltages have different magnitudes, the difference between simulation and experiment in the least-square formulations F j ( j ), ∀j ∈ {1, 2, 3} in Eqs. 17a to 19b is normalized vs. the maximum of the respective experimental values. Further, the parameter vectors j , ∀j ∈ {1, 2, 3} are normalized vs. their starting values for numerical smoothness. For exchange current densities, a logarithmic scale is applied, as the sensitivity of the exchange current density is known to be smaller than the sensitivity of, e.g., diffusion coefficients. This approach is in accordance with the work of, e.g., Jin et al. [17].
As kinetic parameters do not affect the simulated open fullcell voltage curve E OCV vs. discharged charge Q, the static parameters * 1 are identified in PE step 1 using solely the OCP measurement for the two electrodes i ∈ {an, cath} and n sample points k: Both deviations in voltage and charge direction of the simulation from experimental points are considered to ensure accordance for all SOCs. E.g., at intermediate SOCs, the OCV curve is flat and a deviation in charge is large compared to a deviation in capacity. In contrast, at low SOC at a steep OCV curve, deviations in voltage are dominant. Summing up both deviations leads to a high accordance for the entire SOC range. Further, this allows to precisely identify the intercalation steps of the graphite anode. In PE step 2, the identified static parameters are used. Hence, a subset of kinetic parameters * 2 is estimated from the C-rate tests with m different C-rates j for the two electrodes, i ∈ {an, cath} , and n equidistant sample points k: In PE step 3, EIS is simulated at 50% SOC. This step is, e.g., in accordance with the approach in Refs. [9] and [8]. The least-square formulation of this step is as follows for the two electrodes i ∈ {an, cath} , and n frequencies in the experiment: EIS spans a broad frequency range and enables to consider the different time constants of the dynamic processes. This allows to distinguish between parameters related to fast electrochemical reactions and slow diffusion processes. Impedance data, thus, lead to further sensitive parameters like double-layer capacitances of both electrodes. These have a negligible impact in C-rate tests, but a significant impact on the impedance spectra.

Electrochemical reference system
The model is parameterized using experimental data of NMC111 vs. graphite cell. The electrodes are made at a pilot plant-scale production line in the Battery LabFactory Braunschweig. For a detailed description of electrode compositions and production processes, it is referred to our previous publication [20]. Geometric properties and transport process parameters of the electrochemical reference system are listed in Table 3.
All tests were conducted in a commercial three-electrode setup of EL-Cells GmbH, which was placed in a temperature chamber at 25 °C. Test protocols of C-rate tests and OCP measurements are given in Ref. [20]. EIS was conducted at 65 logarithmic-equidistant frequencies ranging from 20 mHz to 50 kHz.

Multi-start identifiability test
In the following, a sample-based identifiability test is introduced. In principle, this test can be applied to any of the three PE steps. Here however, it is solely applied to PE step 2 as this is the main step having the most adjustable parameters. This decision will be justified by the results of Sects.

and 3.3.
Direct approaches to access identifiability require a closed formulation of the model equation. The P2D model, however, does not satisfy this requirement. A solution to get a closed formulation could be a drastic simplification of the model and linearization. However, the nonlinearity bears a (19b) significant amount of information which enhances the identifiability of the model parameters significantly as shown by Lenze et al. [22]. Thus, an indirect, sample-based approach is chosen: multi-start parameter estimation. In this approach, different starting points m,0 , in the domain of model parameters m , are chosen for the parameter estimation algorithm. The domain m,j of parameter j is limited by physically plausible parameter values x lb,j at the lower boundary (lb) and x ub,j at the upper boundary (ub) provided by the literature. The choice of the various starting points can be guided by different approaches. The parameter space could be sampled in an equidistant mesh. This would lead to a large number of sample points and to many unlikely parameter combinations. Random sampling could reduce the number of sample points but would lead to non-deterministic results. Further, a small number of parameter vectors could be chosen by physical insight. This could reduce the number of chosen unlikely parameter combinations, such as that all parameters are at the lower boundary, which would lead to negligible dischargeable capacity. To minimize computational cost, this approach is applied as explained in the following.
Not all parameters show a similar sensitivity towards C-rate curves, so it is focused on the sensitive parameters, as identified in our previous study [20]. Namely this is solidphase diffusion coefficient of both electrodes, electric phase conductivity of both electrodes, and the exchange current densities of both electrodes. Starting values are varied in the following for those parameters.
For each dimension j of the parameter space m three values are defined: the boundaries x j,lb and x j,ub and the reference starting value x j,ref listed in Table 3. From those three values, three potential starting points x j,l,0 , ∀l ∈ {−, 0, +} are derived: The applied logarithm gives weight to the fact that literature values for some parameters range over several orders of magnitude, and thus, the starting values should do as well. The sample generator function, Eqs. 20a to 20c, should be adjusted to the individual parameter estimation problem. As three potential starting values for six parameters would lead to a total of 729 different starting points, promising starting points for PE step 2 are chosen by combination of following deterministic rules: -At least one parameter is set to x j,−,0 , -At least one parameter is set to x j,+,0 , -At least three parameters are set to x j,0,0 .
This reduces the number of starting points to 151. The reference starting point x j,0 = x j,0,0 , ∀j is added as starting point number 152. Further, starting points will be rejected if the discharge capacity at 0.5C is below 50% of the experimental value. If a starting point is not rejected, parameter estimation step 2 is conducted as introduced above.

Point estimate method
In Sect. 3.2, a sensitivity analysis is conducted. A nested point estimate method (PEM) is applied. This method is a tool for global sensitivity analysis which is based on deterministic sampling and Taylor series expansion. For details about the method, it is referred to Refs. [35] and [20]. Sensitivity analysis was shown to be feasible to assess accuracy of parameters, e.g., by Refs. [32] and [43].

Results and discussion
In this section, parameter estimation steps 1 to 3 are conducted. For steps 1 and 3, identifiabily is discussed briefly, while for PE, step 2 multi-start fitting is applied. The focus is on identifying a unique parameter set of the P2D model for simulation of quasi-static electrochemical tests, as this was a common combination as shown above.

Identifiability of full-cell open-circuit potentials
The essential part of PE step 1 is to quantify initial capacity and active material losses during the formation at the end of cell production. PE step 1 affects four parameters: c c,0 , c a,0 , c c,max , and c a,max . Only the measurement of halfcell open-circuit potentials is considered. As OCPs are formulated as functions vs. intercalation ratio, respectively, normalized concentration c s ∕ c max , c max can be used to stretch the half-cell voltage of one electrode compared to full-cell SOC. The initial concentration c s,0 of one electrode can be used to shift one electrode vs. full-cell SOC. However, PE step 1 is quite straightforward as the model interpolates the experimental OCP curve with solely an altered base. Basically, PE step 1 is a fine tuning of the altered base.
Results of PE step 1 are listed in Table 4. It is noteworthy that all four parameters describe the technical usable concentration range, analog to full-SOC, not the entire stoichiometric range of LiC 6 and LiNi 1∕3 Mn 1∕3 Co 1∕3 . For some active materials, the later can be significantly higher than the technical usable range.
To access uniqueness of the estimated parameter set, it is referred to the characteristic shape of both electrode OCP curves, respectively, their nonlinearity. Characteristics of any graphite-based anode's OCP are a wide flat range with discrete step and a steep section at high discharge capacities (low SOCs). In contrast, the NMC cathode's OPC is decreasing with increasing discharge capacity, while the slope is changing as well. For illustration, it is referred to the experimental OCP curves in Fig. 5 in the Online Appendix. The nonlinearity of the OCPs and the availability of half-cell potentials allow identification of c j,0 and Δc j,max for both electrodes. As this approach is based on the characteristics of the OCP curves, it is potentially only valid for the given well-balanced NMC vs. graphite cell. For cells with different OCP characteristics or a Fig. 1 Histogram of the normalized residuum F 2 of the entity of parameter estimations which fulfilled the initial requirement of at least 50% of the experimental capacity at 0.5C Table 5 Exemplary parameter sets from multi-start parameter estimation Green highlights a positive change of the parameter between set 1 and set 2, red highlights a negative change. Simulated C-rate tests for both parameter sets are shown in Fig. 6 in the Online Appendix different balancing, a different approach can be more suitable. Also, analytic approaches at distinct points may be possible.

Practical identifiability from C-rate tests
Multi-start parameter estimation is conducted for PE step 2 as described in Sect. 2.4.
Out of 152 initial points, 49 have been chosen by the algorithm due to the rejection criterion introduced in Sect. 2.4. Further, results with a normalized residuum below 6.0 × 10 −3 were chosen for the further analysis. Consequently, the number of parameter vectors is reduced to 35. Figure 1 shows the normalized residuum of the parameter estimation defined in Eqs. 18a and 18b. There is a narrow band of parameter sets with a residuum between 5.8 × 10 −3 and 6.0 × 10 −3 . Thus, they only vary about 3.3 % and all are equivalent to an excellent accordance between simulation and experimental C-rate test. Thus, the least-square solver is able to converge for a large number of initial points. In this context, solver convergence implies that any parameter set is found which leads to a minimal residuum of the objective function. In the following, it is assessed whether the solver converges to the same parameter set.
For further discussion of the P2D model's unidentifiability in PE step 2, two parameter sets are handpicked. Both sets lead to an excellent accordance between the experimental discharge curve and the simulation but show significant differences in their parameter values. The two parameter sets are listed in Table 5. The corresponding simulated discharge curves are shown in Fig. 6 in the Online Appendix. Differences between the curves are negligible and below the typical experimental accuracy.
Despite their equal discharge curves, the two different parameter sets vary significantly regarding the cathode properties. The cathode solid-phase diffusivity is decreased between the first and second parameter sets by a factor of two, solid-phase conductivity by a factor of 450, and the exchange current density is increased by a factor of 70. At the anode, only the solid-phase conductivity is changed significantly.
Further, a 2D parameter variation is conducted wherein cathode exchange current density is varied (first dimension) and cathode solid-phase conductivity and simultaneously cathode solid-phase diffusivity (second dimension) are varied simultaneously. It is ensured that both parameter sets in Table 5 are a sample of this parameter variation. In Fig. 2, they are highlighted by a circle and a square, respectively. Results of the parameter variation are shown in Fig. 2. It displays F 2 ( * 2 ) for all samples. The parameter variation reveals the unidentifiability of the P2D model with C-rate test data, as an infinite number of solutions is visible. Thus, the combined impact of diffusion and conduction in cathode particles is not distinguishable from the impact of the surface reaction. To overcome this limitation, EIS data are used for the model identification in the following section. The different time constants of reaction and transport processes are used for separation of   Table 5. Point estimate method was applied for six normally distributed uncertain parameters s,c in S m −1 0.18 0.11 [15] i 0,c in A m −2 7 1.20 [43] cathode exchange current density and cathode solid-phase conductivity.
To further illustrate the cause of unidentifiability, a sensitivity analysis is carried out for the two different parameter sets, both providing an optimal accordance of simulated and measured C-rate test. Sobol indices are determined for the sensitivity of the residuum applying PEM as introduced in Sect. 2.5. Results are displayed in Fig. 3.
The model is most sensitive to D s,a , followed by either s,c or i c,0 in dependence on the case. This reveals the transition between the two arcs of the L-shape of the minimal residuum in Fig. 2. One arc relates potential losses to the surface overpotential of the cathode, while the other arc relates them to ohmic losses in the cathode. Sobol indices of D s,c and i a,0 are marginal, and for s,a , they are zero for both cases. Due to that, the significant changes of anodic solid-phase conductivity between the two parameter sets do not affect the discharge curve of the cell, respectively, the residuum F 2 . Further, this leads to practical unidentifiability of the anode conductivity due to non-sensitivity. This further explains the large difference of anode conductivity in Table 5. To enable identifiability of the anode conductivity, a further experiment could be added or the conductivity could be measured ex situ.
Concluding PE step 2, out of the adjusted parameters, the anode solid-phase conductivity remains unidentified as it is insensitive for the investigated C-rates. Further, cathode exchange current density and cathode solid-phase conductivity remain unidentified due to their ambiguous impact on the discharge curve. To overcome this ambiguity, PE step 3 is conducted in the following. The reader is not commended to take PE step 2 as a step-by-step guide for parameter estimation as the optimal parameter estimation procedure may be dependent on applied active materials and electrode design. However, it is emphasized to the different tools, alongside others multi-start parameter estimation and sensitivity analysis, to assess reliability of the estimated parameter set.

Identifiability by addition of assessment of impedance spectra
Electrochemical impedance spectroscopy reveals real and imaginary parts of the battery impedance Z at different frequencies of the sinusoidal current or voltage input signal. At different frequencies, the impedance is governed by different processes and semicircles are related to electrochemical reactions at particle surfaces or in the SEI [15]. Figure 7 in the Online Appendix shows the simulated impedance spectrum of the lithium-ion battery at 50% SOC for illustration. Commonly, the cathodic impedance contains a single semicircle representing electrochemical reactions and charge of double layers, a low-frequency arc representing slow diffusion processes, and a positive non-zero real part at an imaginary part of zero at high frequencies. This point is related to the conductivity of electrolyte in cathode and separator. The full-cell impedance spectrum contains beside anode and cathode semicircle a further semicircle at high frequencies, which is related to the SEI [15]. Applying the SPM and PE step 3, solid-phase conductivity s,c and exchange current density i 0,c are estimated as listed in Table 6. Simulation and experimental data of the identified cathode are shown in Fig. 8 in the Online Appendix. The linear part of the spectra at high frequencies is an artifact from the lithium reference electrode. Also, a cathode double-layer capacity of 4.18 F m −2 was estimated.
Reference values from the literature are also listed in Table 6, which are qualitatively in good accordance to the value determined based on simulations. It should be pointed out that Heins and Schröder used exactly the same electrodes as in this work [15]. Vazquez-Arenas had different electrodes but at least the same active material [43]. The parameters derived from EIS are in good accordance with one of the infinite parameter sets derived from the C-rate test. Thus, this parameter set indicated by the black square is considered as the final parameter set of the P2D model after PE step 3. Cathode diffusivity is taken from PE step 2 from the parameter set matching s,c and i 0,c of PE step 3. This approach is based on the results of Beelen et al. [2]. They suggest that for parameterization, tests should be used which are similar to the later application. In this case, a strong focus lies on the P2D model and C-rate tests as the most common combination.
In this study, the objective of PE step 3 is to distinguish between cathode solid-phase conductivity s,c and cathode exchange current density j 0,c , as the subordinated objective is to parameterize the P2D model for C-rates. Further parameterization of the SPM model and its SEI parameters are beyond the scope of this work. However, for different application a wider parameter study on the impedance data may be beneficial.
To illustrate the feasibility of EIS to do this, the influence of parameter variations of s,c and j 0,c are shown in Fig. 4 for the cathode impedance spectrum.
With increasing solid-phase conductivity, the semicircle moves to lower real parts of impedance, while the size of semicircle stays the same. If the exchange current density is increased, the semicircle, related to cathode intercalation reaction, is decreasing in size, while the abscissa intercept at high frequencies remain the same. Thus, with position and size of the cathode semicircle, the two parameters s,c and j 0,c can be distinguished. It is noteworthy that the abscissa intercept also depends on the electrolyte conductivity. However, this is not an adjustable parameter of the applied model. The electrolyte diffusivity is identified in PE step 2 and the electrolyte conductivity is derived therefrom by the Nernst-Einstein equation.
Further, it is referred to studies on the structural identifiability of a linearized SPM from EIS data of Bizeray et al. They concluded that besides the two solid-phase diffusion coefficients, a lumped charge-transfer resistance could be identified [3]. However, it has to be noticed that their studies were based on full-cell data only. Thus, if half-cell potentials are available, as in our case, separate charge transfer resistances of both electrodes are identifiable.
In conclusion, a three-step parameter estimation procedure was conducted. Out of the 14 adjustable parameters of the P2D model only, the anode solid-phase conductivity remains unidentifiable as it is insensitive at the investigated C-rates. In this case, a reasonable value from literature can be chosen.

Conclusions
In this article, a P2D model was parameterized in a threestep procedure, and practical identifiability was assessed. Unique parameters of the P2D model were successfully identified through a combination of OCP, C-rate, and EIS data in a three-electrode setup. It was shown that the classical C-rate test in combination with OCP measurements does not provide sufficient information for P2D model identification. This result is remarkable as the combination of C-rate tests and P2D models is the most common case in literature. Specifically, reaction and conductivity coefficients cannot be distinguished. This may lead to wrong conclusions on limiting processes and losses in the cell. Further, it may mislead cell and electrode optimization to a potentially wrong direction of attention. Further in this study, the anodic solidphase conductivity is unidentifiable due to non-sensitivity in C-rate tests. This is related to the commonly high conductivity of graphite anodes. The conducted study leads to general guidelines for model parameterization concerning identifiability. First, parameterization of an electrochemical model requires a combination of dynamic and constantcurrent measurements. Second, half-cell and full-cell measurements are required. The parameterization procedure and insights in this work can provide guidance for researchers and developers to conduct a proper parameterization of their models. It gives helpful advice and sensitizes the issue of non-uniqueness.