Numerical Modeling of a Low-Cobalt All-Solid-State Cell with Ceramic Electrolyte Using a Deformable Geometry

: All-solid-state batteries with a lithium negative electrode and a ceramic electrolyte are key toward high energy density. To ensure a safe, fast, accurate, and cost-effective development of this technology, the experimental methodology must be supported by the numerical modeling approach. This work proposes and describes an electrochemical model of a Li 7 La 3 Zr 2 O 12 (LLZO) and Ni-rich NMC-based lithium cell with a deformable lithium negative electrode. Simulations were computed using the finite element method at different operating conditions to demonstrate the scope of the modeling work. Discharge rate tests, deformation tracking, geometric defect investigation, and polarization decomposition are described. Theoretical validation of the mass balance, the stripping rate, the ohmic polarization, and the mesh deformation demonstrated the consistency of the volumetric deformation strategy. We demonstrated in this study a deformable modeling strategy, which was found to be useful for the electrostripping analysis of anodic geometry defects during discharge. Non-uniformity in the lithium stripping rate was found along the anodic interface with defects, and this non-uniformity was accentuated with a higher discharge rate. The cell’s discharge potential was decomposed by considering the equilibrium potential and the polarizations of the main components of the cell. This post-processing was found to be useful for the understanding of the cell’s behavior.


Introduction
One way of improving lithium secondary batteries is to use lithium metal instead of conventional graphite as the negative electrode material, which increases the theoretical specific capacity of the negative electrode by a factor of ten (3861 mAh g −1 versus 372 mAh g −1 for a lithiated graphite electrode).A further way of improving lithium metal batteries is to use a solid electrolyte for enhanced safety and capacity retention.On one hand, this eliminates the need for flammable liquid electrolytes, and on the other, it mechanically impedes lithium dendrite growth [1][2][3][4][5][6][7][8].
The technological development approach to new batteries is generally experimental, which comes with a few shortcomings.First, numerous preparation steps and multiple tests are required to research, develop, validate, or optimize a new generation of materials, processes, or cells.Second, technology scale-up, for example, from button cell to pouch cell, can be a limitation when performing only experiments.Third, common characterization techniques are not very informative about the internal behavior of cells.Ultimately, a deeper understanding of electrochemistry through numerical modeling would ensure safer, faster, more accurate, and more cost-effective technological development [9][10][11].
In the 1990s, Doyle, Fuller, and Newman simulated the charging and discharging of Li-ion cells [12,13].Lithium insertion into the active material was described following Newman's work on porous electrode theory [14].In 2010, Danilov et al. presented the modeling of an all-solid-state Li/Li 3 PO 4 /LiCoO 2 cell and the decomposition of the total polarization.Three polarizations were evaluated: mass transfer polarization for the electrolyte and charge transfer and diffusion polarizations for the positive electrode [15].In 2021, Chen et al. presented a detailed breakdown of the total polarization of a graphite electrode in a half-cell with a mobile anions electrolyte and a conductive additive in the negative electrode.Four polarizations were evaluated: mass transfer polarization for the electrolyte and also charge transfer, diffusion, and ohmic resistance polarizations for the electrode [16].In 2023, Wang et al. presented a detailed analysis of resistances in current collectors [17].In 2023, Xiong et al. presented a detailed polarization decomposition method for lithium-ion battery porous electrodes [18].Several other authors have published electrochemical modeling studies of lithium batteries with liquid or solid electrolytes [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38].Among others, they have studied the thermo-mechanical behavior, predicted the performance, simplified the multiphysics models, and computed the Li transports for all-solid-state and Li-ion battery technologies.However, few authors have focused on modeling all-solid-state batteries with ceramic electrolytes, and only a little work was done on the decomposition of the total polarization of these batteries.Also, these studies do not account for the common geometric defects on the surface of the anode or solid electrolyte.
Currently, the transport electrification industry is considering using all-solid-state batteries and is looking for low-cost technologies that minimize the content of Co in the positive electrode.In this context, the present modeling work was focused on a Li/LLZO/(LLZO-NMC811) cell, where NMC811 is a Ni-rich layered lithiated nickelmanganese-cobalt oxide (LiNi 0.8 Mn 0.1 Co 0.1 O 2 ).The layered lithiated cobalt oxide or LCO (LiCoO 2 ) was also considered a comparative material for the cell's voltage curves.The deformation of the negative electrode was tracked through the discharge in a ceramic electrolyte all-solid-state cell.The deformable geometry strategy was applied to investigate the electrostripping of anodic interface with common geometric defects.A decomposition of polarizations is also presented.This study describes in detail an innovative approach to predict the behavior of all-solid-state batteries (ASSBs) and a new approach that can support the development of technologies incorporating metallic lithium and solid electrolytes.

Methods
The all-solid-state cell, schemed in Figure 1, was divided into three main components: the lithium metal (negative electrode), the solid electrolyte, and the positive electrode consisting of a mixture of active material and solid electrolyte.During discharge, lithium was stripped from the negative electrode and transported to the positive electrode, where it became inserted in particles of active material.

Governing Equations
A numerical model was developed based on the following assumptions: (1) a fully reversible electrochemical cell, (2) isothermal, (3) a variable-volume metal negative electrode, (4) a solid electrolyte with immobile anions, (5) a positive electrode of invariable volume in which spherical active particles are homogeneously mixed with solid electro-

Governing Equations
A numerical model was developed based on the following assumptions: (1) a fully reversible electrochemical cell, (2) isothermal, (3) a variable-volume metal negative electrode, (4) a solid electrolyte with immobile anions, (5) a positive electrode of invariable volume in which spherical active particles are homogeneously mixed with solid electrolyte, (6) with perfectly conducting current collectors, and (7) we considered a perfect interface with no interfacial resistance at both LLZO/Li and LLZO/NMC boundaries, as other studies did [36,49].This assumption was valid as we simulated a cell with a non-degraded interface (without SEI, without additives, and so on).Additionally, LLZO-based solid electrolytes are often prepared using methods similar to those of the semiconductor industry [36], where there is excellent contact between the electrode and electrolyte surfaces.The symbols used in this model are shown in Table 1.The governing equations of the multiphysics numerical model are presented in Table 2.

Main Component Governing Equations Comments
Negative electrode Fick's first law of diffusion ∂c Li + ∂r r=0 = 0 Centrale symmetry condition Lithium insertion in oxide (II) Butler-Volmer equation

Adjustment of exchangecurrent density
For the negative electrode, since the material was solid, electronically conductive, and ionically insulating, the only charge transport mechanism was the electronic conduction described by Ohm's law.The conservation of electric charge was considered stationary, given the speed order of electronic conduction relative to the cell's cycling timescales.The redox half-reaction was that of lithium, and its kinetic was described by the Butler-Volmer relation.
For the solid electrolyte, the lithium cation was the only mobile ion.Oxygen anions, on the other side, were immobile due to their high activation barrier.The dominant mass transport mechanism of lithium cations was migration [50].There was neither polarization nor limiting current, as is the case for electrolytes with diffusive mass transport.The only polarization was ohmic.Since ionic conduction was much higher than electronic conduction, electronic conduction was negligible.The governing equations were the same as for electronic conduction in the negative electrode but this time considering the ionic conductivity of the solid electrolyte.
For the positive electrode, solid electrolyte was uniformly mixed with spherical particles of active material.Three sets of governing competing equations must be considered and applied over the entire electrode domain: one set of equations associated with electronic conduction in the active material, one with lithium migration into the solid electrolyte, and one with lithium insertion into the active material.Bruggeman's relationship applied to the efficient electronic conductivity of the active material of the positive electrode and the efficient ionic conductivity of the electrolyte.Within the active material particles, the lithium-ion transport mechanism was diffusion.The quantity conserved was the mass of lithium cations.Fick's first law describes the internal molar flux.Two boundary conditions were required.At the center of the particle, i.e., at r = 0 in the spherical coordinate system, the condition of central symmetry applied.At the surface of the particle, i.e., at r = r p , i.e., at the interface between the active material and the solid electrolyte, the diffusion-insertion flux equivalence condition applied.The Butler-Volmer relationship described the chemical kinetics of lithium insertion into the active material.The exchange current density was adjusted according to the local lithiation state.

Deformation of the Negative Electrode
During discharge, the lithium metal negative electrode thinned; during charging, it thickened.This volume change was considered in the numerical model by a deformable negative electrode domain and a deformable negative electrode surface.The strategy for deforming the geometry is illustrated in Figure 2. The absolute plating or stripping rate was proportional to the anodic interface current density (Equation ( 1)), where the proportionality constant was κ = M Li /(nFρ Li ).There were a few other mechanical considerations within the adopted model that should be noted.Mechanical stresses were not considered while the cell was subjected to volumetric and geometric change upon charge and discharge.An external stack pressure is generally required to avoid void formation between the lithium and the electrolyte films.This void formation was not considered in the present study.The volume change in the positive electrode was not considered either.
Batteries 2024, 10, x FOR PEER REVIEW 6 of 16 Since the mass transport mechanism in the solid electrolyte was migration, when the cell was discharged, the amount of lithium stripped from the negative electrode must equal the amount of lithium inserted into the positive electrode.To validate the movingboundary strategy, these two quantities were used to post-process the areal discharge capacity ([=] mAh cm −2 ) in two ways.The mass balance was validated by the equality of the two calculated capacity values.The surface discharge capacity according to the lithium stripped from the negative electrode and the surface discharge capacity according to the Since the mass transport mechanism in the solid electrolyte was migration, when the cell was discharged, the amount of lithium stripped from the negative electrode must equal the amount of lithium inserted into the positive electrode.To validate the movingboundary strategy, these two quantities were used to post-process the areal discharge capacity ([=] mAh cm −2 ) in two ways.The mass balance was validated by the equality of the two calculated capacity values.The surface discharge capacity according to the lithium stripped from the negative electrode and the surface discharge capacity according to the lithium inserted in the positive electrode were, respectively,

Polarization Decomposition and Evaluation
The polarizations observed in an electrochemical cell were linked either to reaction activation, electrical or ionic resistance, or charge carrier depletion.The total polarization represents the potential difference between the electrochemical reaction's ideal state and the current state.Furthermore, the equilibrium potential of the positive electrode varies according to its state of charge and temperature.The more lithiated the positive electrode, or the higher the temperature, the lower its oxidizing potential.
Two activation polarizations were observed in the cell under study.According to the Butler-Volmer equation, a kinetic polarization at the negative electrode, η ct − , varied according to the imposed current.It was evaluated at the interface of the negative electrode and the electrolyte.At the positive electrode, also following the Butler-Volmer equation, a kinetic polarization, η ct + , varied according to the imposed current, but this time the exchange current varied according to the local lithiation state of the electrode.It was evaluated at the interface of the electrolyte and the positive electrode.
Four ohmic polarizations were observed at (i) the negative electrode, (ii) the electrolyte, (iii) the active section of the positive electrode, and (iv) the electrolytic section of the positive electrode.The first two were evaluated by the potential difference in their domain.The third could be evaluated by ∆ϕ s + .However, the fourth could not be evaluated by ∆ϕ l + .The effective oxidation potential of the active material was that at the particle surface.Here, the concentration of lithium inserted differed from the concentration averaged across the whole electrode.Thereby, the ohmic polarization of the electrolytic section of the positive electrode could be evaluated as follows.
A charge carrier depletion polarization was observed.At the positive electrode, inside the active material particles, lithium diffused according to Fick's second law.If the flux and particle radius were too high compared with the diffusivity of lithium in the active material, a concentration gradient set in.Here, the concentration of inserted lithium was different from the concentration averaged inside the particle.Thereby, the concentration polarization associated with this phenomenon was as follows.
The total polarization, η tot , was equal to the sum of the absolute polarizations.
The cell potential was the difference between equilibrium potential and total polarization.
The total difference potential was the difference between the actual potential and the fully charged cell potential.

Simulation
The numerical model was solved following the finite element approach available in Comsol Multiphysics 6.1 [51].The following considerations were taken.For the cell boundary conditions, a current density was imposed at one electrode, and an electrical ground was imposed at the other electrode.In the negative electrode, the geometry and the mesh were deforming through the discharge; in the electrolyte film and in the positive electrode, the geometry and the mesh were kept constant.The current distribution was initialized with a primary-type current distribution strategy.The time-dependent study was solved with the fully coupled parallel direct solver PARDISO.The LiionDB database was consulted to assist in selecting the appropriate parameter values [52].Some parameters could have a range of values, e.g., ionic conductivity of LLZO varying from 0.1 and 0.01 S m −1 .To illustrate the effect of the ionic conductivity of LLZO, we simulated the discharge potential curve at 1 C in Figure S1 of the Supplementary Materials.For this study, we chose a typical value of 0.1 S m −1 , ensuring this aligned with practical applications.The parameters used in simulations are presented in Table 3.Since the cell's capacity was limited by polarizations, a stop condition ensured that the discharge ended when the potential fell below 2.5 V.

Discharge Rate
We first focused on the discharge curves of the modeled all-solid-state cell with NMC811 as the positive electrode active material at different discharge rates, and we compared this result with a modeled cell with LCO as the positive electrode active material.The governing equations (Table 2) were solved with nominal parameters (Table 3) at discharge rates of 0.2 C, 0.5 C, 1 C, 2 C, and 5 C. The simulated discharge voltage profiles are presented in Figure 3.

𝜀
Volume fraction of solid electrolyte in the positive electrode 0.

Discharge Rate
We first focused on the discharge curves of the modeled all-solid-state cell with NMC811 as the positive electrode active material at different discharge rates, and we compared this result with a modeled cell with LCO as the positive electrode active material.The governing equations (Table 2) were solved with nominal parameters (Table 3) at discharge rates of 0.2 C, 0.5 C, 1 C, 2 C, and 5 C. The simulated discharge voltage profiles are presented in Figure 3.The equilibrium potential, which indicated the potential of an ideal cell, thus without polarization or at a very low discharge rate, was also calculated and presented in the same figure [56,57].Here, we see that even polarization-free or at very low discharge rates, the capacity and the cutoff voltage of a cell with an NMC-based positive electrode were lower than those with an LCO-based positive electrode.The theoretical areal capacity was calculated based on the lithium insertion in the positive electrode (NMC811) and was limited by polarizations, which will be discussed further.

Deformation of the Negative Electrode
During a discharge of the modeled all-solid-state cell, we monitored the deformation of the negative electrode with a plane interface while using NMC811 as the active material The equilibrium potential, which indicated the potential of an ideal cell, thus without polarization or at a very low discharge rate, was also calculated and presented in the same figure [56,57].Here, we see that even polarization-free or at very low discharge rates, the capacity and the cutoff voltage of a cell with an NMC-based positive electrode were lower than those with an LCO-based positive electrode.The theoretical areal capacity was calculated based on the lithium insertion in the positive electrode (NMC811) and was limited by polarizations, which will be discussed further.

Deformation of the Negative Electrode
During a discharge of the modeled all-solid-state cell, we monitored the deformation of the negative electrode with a plane interface while using NMC811 as the active material for the positive electrode.The governing equations (Table 2) were solved with nominal parameters (Table 3) at a discharge rate of 1 C.The simulation results are presented in Figure 4.A video of the deforming mesh coupled with the voltage profile during the discharge is available in Supplementary Materials.
To validate the deforming geometry strategy, four elements are observed: (i) the mass balance, (ii) the stripping rate, (iii) the ohmic polarization, and (iv) the mesh deformation.Due to the immediacy of the migration mass transfer mechanism, there was no accumulation of lithium in the electrolyte.Hence, the quantity of stripped metallic lithium from the anode surface equaled the quantity of inserted ionic lithium in the cathode active material.The two methods proposed for calculating the areal capacity during the discharge (Equations ( 2) and (3)) were followed (Figure 4a).They were in line with each other and validated the mass balance.The stripping rate (Figure 4b) and the anodic ohmic polarization (Figure 4c) were probed during the discharge.The stripping rate was positive and constant.Also, the anodic ohmic polarization was negligible, positive, and decreasing.These observations aligned with theoretical considerations of Equation ( 1) and Ohm's law.The mesh was actualized at each time step from the beginning of the discharge (Figure 4d) until its end (Figure 4e).The negative electrode was contracting while the rest of the cell remained unchanged.At the end of the discharge, the cell was almost anode-free with only a 10 µm film of lithium negative electrode left as a cathodic reactive surface for the charging step.The mesh contraction at the end of the discharge, when there was a sudden drop of voltage, was consistent with the mass of lithium migrated from the anode to the cathode.
Batteries 2024, 10, x FOR PEER REVIEW 9 of 16 for the positive electrode.The governing equations (Table 2) were solved with nominal parameters (Table 3) at a discharge rate of 1 C.The simulation results are presented in Figure 4.A video of the deforming mesh coupled with the voltage profile during the discharge is available in Supplementary Materials.To validate the deforming geometry strategy, four elements are observed: (i) the mass balance, (ii) the stripping rate, (iii) the ohmic polarization, and (iv) the mesh deformation.Due to the immediacy of the migration mass transfer mechanism, there was no accumulation of lithium in the electrolyte.Hence, the quantity of stripped metallic lithium from the anode surface equaled the quantity of inserted ionic lithium in the cathode active material.The two methods proposed for calculating the areal capacity during the discharge (Equations ( 2) and (3)) were followed (Figure 4a).They were in line with each other and validated the mass balance.The stripping rate (Figure 4b) and the anodic ohmic polarization (Figure 4c) were probed during the discharge.The stripping rate was positive and constant.Also, the anodic ohmic polarization was negligible, positive, and decreasing.These observations aligned with theoretical considerations of Equation ( 1) and Ohm's law.The mesh was actualized at each time step from the beginning of the discharge (Figure 4d) until its end (Figure 4e).The negative electrode was contracting while the rest of the cell remained unchanged.At the end of the discharge, the cell was almost anode-free with only a 10 µm film of lithium negative electrode left as a cathodic reactive surface for the charging step.The mesh contraction at the end of the discharge, when there was a sudden

Electrostripping of Geometric Defects in the Anodic Interface
There would be no significant difference between the results obtained from a 1D or 2D model regarding the ohmic potential drop when using a planar anodic interface.However, one insight of the 2D deformable geometry strategy was the capability to study the transient electrostripping of defects in the anodic interface.
The deformation strategy presented in Figure 2 was extended with the use of a volumeconstant deformable single-ion electrolyte.The electrostripping of a lithium anode with a peak and a valley, two common anodic interface defects, was investigated.It was assumed that the submicron defects were thin enough to allow perfect interfacial contact during the simulated discharges.The governing equations (Table 2) with NMC811 were solved with nominal parameters (Table 3) at a discharge rate of 1 C. Figure 5 shows the interface deformation progression.The anodic interface x-position is shown in time with three proportional variables, which we identify.These are the electrolyte film potential, the normal current density at the interface, and normal stripping rate.
during the simulated discharges.The governing equations (Table 2) with NMC811 were solved with nominal parameters (Table 3) at a discharge rate of 1 C. Figure 5 shows the interface deformation progression.The anodic interface x-position is shown in time with three proportional variables, which we identify.These are the electrolyte film potential the normal current density at the interface, and normal stripping rate.It was demonstrated that on a lithium peak, the electrolyte film potential, the norma interface current density, and the normal stripping rate were maximized on the tip and minimized on the edges.It was also demonstrated that the opposite phenomenon occurred in a lithium valley: each of these proportional variables was minimized in the base and maximized on the edges.The geometric defects on the lithium made the stripping rate non-uniform, and this non-uniformity along the interface, while reduced during the discharge, was maximal at the initial time step.
Figure 6 shows the initial value of the stripping rate at different discharge rates.It is shown that the non-uniformity of the initial stripping rate along the interface was accentuated at higher discharge rates.The locations around a geometric defect with the highes stripping rate were on the edges for a valley and on the tip for a peak.If the interfacia contact was not maintained, it was at these locations that voids would be inclined to form This could occur if the lithium flow from inside the electrode was locally insufficient or i It was demonstrated that on a lithium peak, the electrolyte film potential, the normal interface current density, and the normal stripping rate were maximized on the tip and minimized on the edges.It was also demonstrated that the opposite phenomenon occurred in a lithium valley: each of these proportional variables was minimized in the base and maximized on the edges.The geometric defects on the lithium made the stripping rate nonuniform, and this non-uniformity along the interface, while reduced during the discharge, was maximal at the initial time step.
Figure 6 shows the initial value of the stripping rate at different discharge rates.It is shown that the non-uniformity of the initial stripping rate along the interface was accentuated at higher discharge rates.The locations around a geometric defect with the highest stripping rate were on the edges for a valley and on the tip for a peak.If the interfacial contact was not maintained, it was at these locations that voids would be inclined to form.This could occur if the lithium flow from inside the electrode was locally insufficient or if the surface of the electrolyte was strictly non-deformable.Recently, Vishnugopi et al. have also reported that a contact loss might appear during stripping [58].This is more inclined to occur at higher discharge rates, i.e., when the stripping rate variability along the interface is accentuated.These results on the effect of different discharge rates are in accordance with the predictions recently presented by Carmona and Albertus in their work on stationary current distribution at the Li/LLZO interface [59].
the surface of the electrolyte was strictly non-deformable.Recently, Vishnugopi et al. have also reported that a contact loss might appear during stripping [58].This is more inclined to occur at higher discharge rates, i.e., when the stripping rate variability along the interface is accentuated.These results on the effect of different discharge rates are in accordance with the predictions recently presented by Carmona and Albertus in their work on stationary current distribution at the Li/LLZO interface [59].

Decomposition of Polarizations
We then decomposed the polarizations of a general case Li/LLZO/(LLZO-NMC811) cell during a discharge.The polarizations were extracted from the preceding section's simulation where the governing equations (Table 2) were solved with nominal parameters (Table 3) at a discharge rate of 1 C.The polarizations of main components are presented in Figure 7, and the order of magnitudes of all polarizations are shown in Table 4.

Decomposition of Polarizations
We then decomposed the polarizations of a general case Li/LLZO/(LLZO-NMC811) cell during a discharge.The polarizations were extracted from the preceding section's simulation where the governing equations (Table 2) were solved with nominal parameters (Table 3) at a discharge rate of 1 C.The polarizations of main components are presented in Figure 7, and the order of magnitudes of all polarizations are shown in Table 4. also reported that a contact loss might appear during stripping [58].This is more inclined to occur at higher discharge rates, i.e., when the stripping rate variability along the interface is accentuated.These results on the effect of different discharge rates are in accordance with the predictions recently presented by Carmona and Albertus in their work on stationary current distribution at the Li/LLZO interface [59].

Decomposition of Polarizations
We then decomposed the polarizations of a general case Li/LLZO/(LLZO-NMC811) cell during a discharge.The polarizations were extracted from the preceding section's simulation where the governing equations (Table 2) were solved with nominal parameters (Table 3) at a discharge rate of 1 C.The polarizations of main components are presented in Figure 7, and the order of magnitudes of all polarizations are shown in Table 4.The sum of polarizations with the equilibrium potential equaled the potential difference.The polarization of the negative electrode, especially the ohmic contribution, was several orders of magnitude lower than the total polarization.The ohmic polarization of the electrolyte between the electrodes was also orders of magnitude lower that the total polarization.The cell's discharge concluded not at 3600 s but at 3595 s with a sudden and significant rise of the polarization in the positive electrode.The capacity was limited by the charge transfer polarization that increased when the surface insertion concentration approached the maximal concentration.
The charge transfer contribution was dominant, showing that the cell was reaction limited.This limitation was linked to the 0.942 maximal state of the charge coefficient of the NMC811, which highlighted the charge transfer polarization.This effect was unique to active materials with a maximal state of charge coefficient lower than one.On an equivalent cell with LCO as a positive active material, with a maximal state of the charge coefficient of 1, the final voltage drop would have a significant diffusion polarization contribution together with an equilibrium potential drop.
One insight of this decomposition of polarization was to examine the impact and the tuning of key design parameters on the cell's voltage profile.Here, we repeated the preceding section's simulation but here with (i) a thicker film of electrolyte between the electrodes and (ii) a larger radius of active particles in the positive electrode.Figure 8 shows that a thicker electrolyte film would lead to an undesired decrease in voltage by an ohmic drop.However, a thinner film of electrolyte, under 50 µm for this simulation case, would not be a key to enhancing the cell's capacity or voltage.The aim of achieving ultrathin ceramic electrolyte film was to increase energy density and specific capacity of the cell due to a weight and volume reduction.The sum of polarizations with the equilibrium potential equaled the potential difference.The polarization of the negative electrode, especially the ohmic contribution, was several orders of magnitude lower than the total polarization.The ohmic polarization of the electrolyte between the electrodes was also orders of magnitude lower that the total polarization.The cell's discharge concluded not at 3600 s but at 3595 s with a sudden and significant rise of the polarization in the positive electrode.The capacity was limited by the charge transfer polarization that increased when the surface insertion concentration approached the maximal concentration.
The charge transfer contribution was dominant, showing that the cell was reaction limited.This limitation was linked to the 0.942 maximal state of the charge coefficient of the NMC811, which highlighted the charge transfer polarization.This effect was unique to active materials with a maximal state of charge coefficient lower than one.On an equivalent cell with LCO as a positive active material, with a maximal state of the charge coefficient of 1, the final voltage drop would have a significant diffusion polarization contribution together with an equilibrium potential drop.
One insight of this decomposition of polarization was to examine the impact and the tuning of key design parameters on the cell's voltage profile.Here, we repeated the preceding section's simulation but here with (i) a thicker film of electrolyte between the electrodes and (ii) a larger radius of active particles in the positive electrode.Figure 8 shows that a thicker electrolyte film would lead to an undesired decrease in voltage by an ohmic drop.However, a thinner film of electrolyte, under 50 µm for this simulation case, would not be a key to enhancing the cell's capacity or voltage.The aim of achieving ultrathin ceramic electrolyte film was to increase energy density and specific capacity of the cell due to a weight and volume reduction.The upper part of Figure 9 shows that a larger radius of the active particles in the positive electrode would lead to a deteriorated voltage profile.The lower part of Figure 9 shows that this was more precisely due to charge transfer activation limitation and also shows that the other polarization components in the positive electrode were not significantly affected.A larger radius of active particles led to undesired decreased capacity and voltage of the cell.A radius smaller than 10 µm for this simulation case did not lead to improvement in capacity but in voltage only.The limitation in capacity occurred due to the rise in charge transfer polarization as the surface insertion concentration drew closer to its maximum level.The lower voltage from the beginning of the discharge was caused by the increased surface area of the active materials.This indicated that in a positive electrode with a particle-sized distributed active material, larger particles could contribute more to the voltage drop than smaller ones.
to improvement in capacity but in voltage only.The limitation in capacity occurred due to the rise in charge transfer polarization as the surface insertion concentration drew closer to its maximum level.The lower voltage from the beginning of the discharge was caused by the increased surface area of the active materials.This indicated that in a positive electrode with a particle-sized distributed active material, larger particles could contribute more to the voltage drop than smaller ones.

Conclusions
This research paper presents a numerical electrochemical model of an all-solid-state cell with lithium metal electrode, a ceramic solid electrolyte, and a low-cobalt positive electrode.The studied ceramic solid electrolyte was LLZO, and the studied positive electrode active material was NMC811.The mobility of the lithium ions in the electrolyte was modeled as migrating in an immobile anions medium.
We proposed a way to consider the negative electrode as deformable and to track its deformation in time.We calculated the areal capacity through discharge from the deformation of the negative electrode and from insertion in the positive electrode.We showed

Conclusions
This research paper presents a numerical electrochemical model of an all-solid-state cell with lithium metal electrode, a ceramic solid electrolyte, and a low-cobalt positive electrode.The studied ceramic solid electrolyte was LLZO, and the studied positive electrode active material was NMC811.The mobility of the lithium ions in the electrolyte was modeled as migrating in an immobile anions medium.
We proposed a way to consider the negative electrode as deformable and to track its deformation in time.We calculated the areal capacity through discharge from the deformation of the negative electrode and from insertion in the positive electrode.We showed that the deformable negative electrode strategy was consistent from viewpoints of (i) the mass balance, (ii) the stripping rate, (iii) the ohmic polarization, and (iv) the mesh deformation.
We showed that this deformable modeling strategy did not affect the voltage profile predictions since the lithium electrode's ohmic polarization was negligible.However, it could be used to illuminate the behavior of geometric defects in the anodic interface.We investigated the cases of a valley and a peak and illustrated the electrostripping of geometric defects in the anodic interface.This simulation case demonstrated that the stripping rate was non-uniform along the interface, and this phenomenon was accentuated at high discharge rates.This work in finite element modeling of the lithium-electrolyte interface, which moves over time, is a step toward a powerful tool for the visualization and optimization of Li/LLZO interfaces.This would be useful to developers of ASSBs looking for robust tools for predicting the performance of this promising technology.

Figure 1 .
Figure 1.The 2D scheme of the all-solid-state cell.

Figure 2 .
Figure 2. The 2D scheme of the geometry deformation strategy.The negative electrode domain and surface were deformable.The domain of the solid electrolyte and that of the positive electrode were constant.

Figure 2 .
Figure 2. The 2D scheme of the geometry deformation strategy.The negative electrode domain and surface were deformable.The domain of the solid electrolyte and that of the positive electrode were constant.

Figure 3 .
Figure 3. Discharge potential curves at different discharge rates with (a) NMC811 and (b) LCO positive electrode active material.

Figure 3 .
Figure 3. Discharge potential curves at different discharge rates with (a) NMC811 and (b) LCO positive electrode active material.

Figure 4 .
Figure 4. Deformation of the negative electrode with a plane interface.(a) Voltage profile calculated from the deformation of the negative electrode (dashed green curve) and from insertion in the positive electrode (dashed blue curve) compared with the profile of an ideal cell.(b) Stripping rate and (c) ohmic negative polarization.(d) Initial state and (e) final state of the deformable mesh.

Figure 4 .
Figure 4. Deformation of the negative electrode with a plane interface.(a) Voltage profile calculated from the deformation of the negative electrode (dashed green curve) and from insertion in the positive electrode (dashed blue curve) compared with the profile of an ideal cell.(b) Stripping rate and (c) ohmic negative polarization.(d) Initial state and (e) final state of the deformable mesh.

Figure 5 .
Figure 5. Progression of the electrostripping of a lithium anode with two geometric defects at the interface: a lithium peak and a lithium valley.Contour lines represent time from 0 to 1 h (1 C) with 30 s time steps.

Figure 5 .
Figure 5. Progression of the electrostripping of a lithium anode with two geometric defects at the interface: a lithium peak and a lithium valley.Contour lines represent time from 0 to 1 h (1 C) with 30 s time steps.

Figure 6 .
Figure 6.Initial stripping rate of anodic interface with defects at different discharge rates.

Figure 6 .
Figure 6.Initial stripping rate of anodic interface with defects at different discharge rates.

Figure 6 .
Figure 6.Initial stripping rate of anodic interface with defects at different discharge rates.

Table 4 .
Order of magnitude of all polarizations.

Table 1 .
Symbols used in the model.Li + ,min Minimal lithium concentration in the active material of the positive electrode mol m −3 c Li + ,max Maximal lithium concentration in the active material of the positive electrode mol m −3 c Li + ,re f Reference lithium concentration in the active material of the positive electrode mol m −3 D Li + Diffusion coefficient of lithium in the active material of the positive electrode m 2 s −1

Table 2 .
Governing equations of the numerical model.

Table 3 .
List of parameters used for simulations.Li + ,max − c Li + ,min nFε s + A / c