## Slow Voltage Relaxation of Silicon Nanoparticles with a Chemo-Mechanical Core-Shell Model

∗

Lukas K¨ obbing, Yannick Kuhn, and Birger Horstmann Institute of Engineering Thermodynamics, German Aerospace Center (DLR), Wilhelm-Runge-Straße 10, 89081 Ulm, Germany Helmholtz Institute Ulm (HIU), Helmholtzstraße 11, 89081 Ulm, Germany and Faculty of Natural Sciences, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany

Silicon presents itself as a high-capacity anode material for lithium-ion batteries with a promising future. The high ability for lithiation comes along with massive volume changes and a problematic voltage hysteresis, causing reduced efficiency, detrimental heat generation, and a complicated stateof-charge estimation. During slow cycling, amorphous silicon nanoparticles show a larger voltage hysteresis than after relaxation periods. Interestingly, the voltage relaxes for at least several days, which has not been physically explained so far. We apply a chemo-mechanical continuum model in a core-shell geometry interpreted as a silicon particle covered by the solid-electrolyte interphase to account for the hysteresis phenomena. The silicon core (de)lithiates during every cycle while the covering shell is chemically inactive. The visco-elastoplastic behavior of the shell explains the voltage hysteresis during cycling and after relaxation. We identify a logarithmic voltage relaxation, which fits with the established Garofalo law for viscosity. Our chemo-mechanical model describes the observed voltage hysteresis phenomena and outperforms the empirical Plett model. In addition to our full model, we present a reduced model to allow for easy voltage profile estimations. The presented results support the mechanical explanation of the silicon voltage hysteresis with a coreshell model and encourage further efforts into the investigation of the silicon anode mechanics.

## I. INTRODUCTION

For the enhancement of next-generation lithium-ion batteries, research and industry consider the application of pure silicon anodes [1-3]. Silicon is a popular choice as it is an abundant and cheap material. Anodes made from silicon possess a high theoretical capacity, leading to a massive volume expansion of up to 300 % during lithiation and respective shrinkage during delithiation [4]. The massive deformations induce significant stresses inside the anode material, causing fracture of large silicon particles above a critical diameter of 150 nm [5]. Larger silicon particles suffer from cracks, particle pulverization, and are prone to losing contact with the current collector [6]. Anodes made from silicon nanoparticles promise a higher stability and cycle life compared to anodes with larger silicon particles. Thus, research and industry focus on the behavior of nano-structured silicon anodes [7].

tion gradients due to slow diffusion in micrometer-sized particles [16-18], phase transformation in the very first cycle [19], slow reaction kinetics [20], and mechanical impact of the solid-electrolyte interphase on the lithiation behavior of silicon nanoparticles [21].

A severe challenge for the commercialization of silicon anodes is the handling and possible reduction of the voltage hysteresis observed in various experiments [8-11]. Silicon anodes reveal a different voltage during slow lithiation compared to delithiation reducing efficiency and causing detrimental heat generation [12, 13]. Experiments observe this hysteresis phenomenon of amorphous silicon anodes in thin-film geometries, micron-sized particles, and nanoparticles. Literature discusses different reasons for the voltage hysteresis: mechanics and plastic flow of silicon in thin-film geometries [14, 15], concentra-

∗ birger.horstmann@dlr.de

The solid-electrolyte interphase (SEI) forms on the surface of the silicon anode particles resulting from the electrochemical instability of the electrolyte in contact with the anode [22-24]. This process leads to a significant capacity loss, even more pronounced for silicon nanoparticle anodes due to the high surface-to-volume ratio. The SEI film is only several nanometers thick but can effectively limit further electrolyte decomposition [25, 26]. Nevertheless, a leakage of electrons through the SEI causes a continued, long-term growth of the SEI during storage and cycling [27, 28]. Concerning the massive volume changes of silicon particles during cycling, it is important to consider the mechanics of the SEI, namely elastic and plastic deformations as well as fracture [29]. However, experimental and theoretical results indicate that the innermost layer of the SEI withstands the large stresses due to plastic flow and fast self-healing [29-31]. In addition to the SEI, a natural silicon oxide layer usually covers the surface of silicon particles and contributes to the innermost SEI layer [32].

For silicon nanoparticles, the SEI has a beneficial effect on their mechanical integrity [1, 33, 34]. Furthermore, carbon-coated silicon particles show a reduced voltage hysteresis [10]. Both observations support our recent explanation of the voltage hysteresis for silicon nanoparticles based on SEI mechanics [21]. Due to the reported pulverization of large particles [6] and growth of SEI into the interior of silicon anodes [35], the mechanical expla-

nation of the hysteresis for silicon nanoparticles covered by a surface layer also has the potential to reason the hysteresis observed for larger geometries [21]. Nonetheless, recent experimental results reveal a slow voltage relaxation behavior for at least 300 h and a significantly enlarged cycling voltage hysteresis even at C / 100 not explained so far [11]. The slow voltage relaxation is in line with the experimental findings of Sethuraman et al. [20]. However, their theoretical explanation with reaction kinetics in the Tafel regime requires extreme values for the exchange current density and the transfer coefficients. Here, we present a consistent picture of the voltage relaxation and other voltage hysteresis phenomena based on mechanics.

This manuscript builds on our previous explanation of the voltage hysteresis of silicon nanoparticles by the mechanical interaction of silicon and SEI [21]. We explain the basic principles of our chemo-mechanical coreshell model in Section II. Furthermore, we introduce the Garofalo viscosity model necessary because of the large stresses arising inside the SEI shell and discuss its behavior in the core-shell system with an analytical approximation and a reduced model. In Section III, we describe the recent experiments performed by Wycisk et al. [11], which we analyze in detail in Section IV. In conclusion, we present a consistent description of the observed slow voltage relaxation, hysteresis shape, C-rate dependence, and voltage transition profiles.

## II. THEORY

Our theoretical framework describes the behavior of a core-shell system, where the silicon particle as core can lithiate and delithiate while the shell is chemically inactive and deforms only mechanically as illustrated in Fig. 1. The core-shell model can in principle relate to different scenarios: (i) SEI on silicon nanoparticle or pulverized silicon micro-particle. Literature reports that the inner SEI on silicon is robust against the severe volume changes [29-31]. (ii) Silicon oxide layer on silicon due to exposure to air or electrolyte, which is reported to affect the silicon lithiation behavior [36, 37]. This effect potentially combines with the influence of the SEI [32]. (iii) Silicon nanodomains in silicon particles surrounded by chemically inactive regions. Literature reports the existence of silicon nanodomains for amorphous silicon under high pressure [38], for crystalline silicon [39], and generically for silicon oxide particles [40, 41]. We interpret the core-shell behavior as a particle-SEI system in the following.

We have presented the foundations of the chemomechanical particle-SEI model used in this study in our previous publications [21, 29]. In the following, we summarize the most important assumptions and equations. Further, we highlight advancements compared to our previous works.

## A. Chemo-Mechanical Core-Shell Model

The silicon particle core deforms due to the chemical lithiation and delithiation F core ch , , elastic deformation F core el , , and plastic deformation F core pl , when reaching the yield criterion. The large deformation approach determines the total deformation F core as

$$F _ { \text{core} } = F _ { \text{core,pl} } F _ { \text{core,el} } F _ { \text{core,ch} }. \quad \quad ( 1 )$$

The concentration of lithium atoms c Li 0 , inside the silicon particle expressed in the undeformed Lagrangian frame determines the chemical deformation

$$\mathbf F _ { \text{core,ch} } = \lambda _ { \text{ch} } \mathbf l d = ( 1 + v _ { \text{Li} } c _ { \text{Li}, 0 } ) ^ { 1 / 3 } \mathbf l d \quad ( 2 )$$

with v Li the molar volume of lithium inside silicon.

The strain tensors E core k , read

$$E _ { \text{core,k} } = \frac { 1 } { 2 } \left ( \mathbf F _ { \text{core,k} } ^ { \text{T} } \mathbf F _ { \text{core,k} } - \text{Id} \right ), \quad \ \ ( 3 )$$

where the subscript k indicates the kind of deformation, which is either the total deformation or one of the mentioned deformation contributions from Eq. (1).

The Cauchy stress σ describes the stress in the deformed Euler frame and the Piola-Kirchhoff stress P = Jσ F -T describes the stress in the undeformed Lagrangian frame with J = det F . The Piola-Kirchhoff stress due to elastoplastic deformation reads

$$\mathbf P _ { \text{core} } = \lambda _ { \text{ch} } ^ { - 2 } \mathbf F _ { \text{core} } \mathbf F _ { \text{core,pl} } ^ { - T } \mathbf F _ { \text{core,pl} } ^ { - 1 } \\ ( \lambda _ { \text{core} } \, \text{tr} ( \mathbf E _ { \text{core,el} } ) \text{Id} + 2 G _ { \text{core} } \mathbf E _ { \text{core,el} } )$$

with the first Lam´ constant e λ core and the second Lam´ e constant G core .

Due to the chemo-mechanical coupling, the stress inside the particle affects the voltage as

$$U = U _ { 0 } + \frac { v _ { \text{Li} } } { 3 F \lambda _ { \text{ch} } ^ { 3 } } \mathbf P _ { \text{core} } \colon \mathbf F _ { \text{core} }$$

with the true open-circuit voltage (OCV) of silicon U 0 and the Faraday constant F .

The differential equations of interest inside the particle are the continuity equation for the time derivative of the lithium concentration ˙ c Li 0 , , the momentum balance, and the equation for the plastic flow rate ˙ F core pl , ,

$$\dot { c } _ { \text{Li}, 0 } = - \nabla _ { 0 } \cdot \vec { N } _ { \text{Li}, 0 } \text{ \quad \ \ } ( 6 )$$

$$0 = \nabla _ { 0 } \cdot { \mathbf P } _ { \text{core} } \quad \quad ( 7 )$$

$$\dot { \mathbf F } _ { \text{core,pl} } \mathbf F _ { \text{core,pl} } ^ { - 1 } = \phi _ { \text{core} } \frac { \partial f _ { \text{core} } } { \partial \mathbf M _ { \text{core} } }. \quad \quad ( 8 )$$

For the lithiation equation, we define the lithium flux ⃗ N Li 0 , = - ∇ L 0 µ Li , the electro-chemo-mechanical potential µ Li = -FU , the mobility L = D Li ( ∂µ Li /∂c Li 0 , ) -1 , and the diffusion coefficient D Li . At the particle boundary, the (de)lithiation rate determines the lithium flux

FIG. 1. Scheme of volume change and SEI stress during lithiation, delithiation, and relaxation periods.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-091d7c29-7169-44d2-b952-e8dea98f588f.jpg)

⃗ N Li 0 , ( R ). For the plastic flow, the von Mises yield criterion f core = 3 2 | M dev core | 2 /σ 2 Y core , -1 ≤ 0 determines plasticity with M dev core = M core -1 / 3 tr M core the deviatoric part of the adapted Mandel stress M core = F T core rev , σ core F -T core rev , , the reversible deformation F core rev , = F core el , F core ch , , and the yield stress σ Y core , . The consistency condition ˙ f core = 0 determines the plastic multiplier ϕ core .

dependence of the Garofalo law is reasoned in Ref. [43] by a change in the energy landscape due to mechanical deformations and lattice distortions. Furthermore, positive entropy production is guaranteed analogously to the derivation in Ref. [21] as the inverse hyperbolic sine is positive for positive arguments and negative for negative ones.

For the shell behavior, we assume that the SEI shell deforms only mechanically, namely elastically and plastically,

$$F _ { \text{shell} } = F _ { \text{shell}, \text{pl} } F _ { \text{shell}, \text{el} }$$

leading to massive mechanical strains and stresses when experiencing the significant volume change of the silicon particle during cycling. Analogous to the particle, Eq. (3) determines the strain tensors E shell k , inside the SEI.

In addition to the elastoplastic stress P shell el , determined analogously to Eq. (4), we consider the viscous behavior of the SEI. To describe large viscous stresses during cycling on the one hand and small viscous stresses during relaxation on the other hand, we use the Garofalo law or inverse hyperbolic sine law

$$\mathbf P _ { \text{shell,visc} } = J _ { \text{shell} } \sigma _ { \text{ref} } \cdot \text{asinh} \left ( \tau \dot { \mathbf E } _ { \text{shell} } \right ) \mathbf F _ { \text{shell} } ^ { - T } \quad ( 1 0 ) \quad \text{ multi} \quad$$

calculated component-wise and presented initially in Ref. [42]. The parameter σ ref describes as a reference stress the magnitude of viscous stress at a certain strain rate. The parameter τ describes the time constant of the system and the dependence on the strain rate. In this study, we use the Garofalo viscosity model stated in Eq. (10) instead of a standard Newtonian model or a shear-thinning model [21] to account more adequately for the complexity of the mechanical behavior. The particular functional

The differential equations of interest inside the SEI shell are the momentum balance and the equation for plastic flow,

$$0 = \nabla _ { 0 } \cdot ( \mathbf P _ { \text{shell}, \text{el} } + \mathbf P _ { \text{shell}, \text{visc} } ) \quad ( 1 1 )$$

$$\dot { \mathbf F } _ { \text{shell,pl} } \mathbf F _ { \text{shell,pl} } ^ { - 1 } = \phi _ { \text{shell} } \frac { \partial f _ { \text{shell} } } { \partial \mathbf M _ { \text{shell,el} } }.$$

The yield criterion f shell = 3 2 | M dev shell el , | 2 /σ 2 Y shell , -1 ≤ 0 is determined by the deviatoric part M dev shell el , of the adapted elastic Mandel stress M shell el , = F T shell el , σ shell el , F -T shell el , and the plastic multiplier ϕ shell results from the consistency condition ˙ f shell = 0.

Note that we model the mechanical deformations on a continuum scale. Thus, the visco-elastoplastic behavior is not necessarily an intrinsic property of a single material domain. Instead, interfaces and grain boundaries of multiple crystal domains can determine the continuum mechanics. Hence, the described visco-elastoplasticity can be a consequence of repeated partial cracking and healing, as discussed for the SEI in Ref. [29]. This description is reasonable because the literature does not observe significant fracture of the inner SEI layer on silicon [29-31].

We assume that the surfaces of the silicon particle and SEI stick tightly together, meaning that the radial part of the stress coincides

$$\mathbf P _ { \text{core,rr} } \Big | _ { r = R } = \mathbf P _ { \text{shell,rr} } \Big | _ { r = R } \quad \quad ( 1 3 )$$

when evaluated at the particle-SEI interface r = R . Due to the merely mechanical deformation of the SEI, significant stresses arise inside the SEI impacting the silicon particle stress and voltage.

As presented in Ref. [21], the expansion of the silicon particle during lithiation leads to a mechanical reaction of the SEI, namely, first elastic and then plastic deformation. The strains inside the SEI generate significant compressive stress acting on the silicon particle as visualized in Fig. 1. Additionally, viscous behavior increases the total compressive stress during lithiation depending on the strain rate. During the subsequent delithiation, tensile stress originates from elastic and plastic deformations as well as viscosity, which causes a stress hysteresis inside the SEI, impacting the voltage of the silicon particle according to Eq. (5). Hence, the visco-elastoplastic behavior of the SEI describes the voltage hysteresis observed for silicon nanoparticles.

## B. Analytical Approximation for the Voltage Relaxation

To gain an analytical approximation for the voltage relaxation, we investigate the behavior of the presented chemo-mechanical core-shell model in a simplified setup. Thus, we analyze all local variables at the interface accounting for the central role of the interface coupling. In the following, we discuss several assumptions paving the way to a simplified analytical expression.

First, we choose the simplified description that during relaxation the silicon particle behaves purely elastically according to Hooke's law

$$\sigma _ { \text{ev} } = E _ { \text{core} } \cdot \mathbf E _ { \text{core,ev}, \text{rr} } \quad \quad ( 1 4 ) \quad \Delta$$

with Young's modulus E core and the elastic strain of the core E core ev rr , , due to viscous stress of the shell.

Furthermore, we consider only the viscous stress contribution inside the shell as the elastic stress of the shell stays constant, i.e.,

$$\sigma _ { \text{shell}, \text{visc} } = \sigma _ { \text{ref} } \cdot \sinh \left ( \tau \dot { \mathbf E } _ { \text{shell} } \right ). \quad \ \ ( 1 5 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The time evolution of the radial stress component in the silicon particle resulting from the time derivative of Eq. (14) states

$$\frac { d \sigma _ { \text{ev} } } { d t } = E _ { \text{core} } \cdot \dot { \mathbf E } _ { \text{core,ev}, \text{rr} }. \quad \quad ( 1 6 )$$

The silicon core deforms only elastically during relaxation of viscous shell stress and isotropically, thus

$$\dot { \mathbf E } _ { \text{core,ev} } \approx \dot { \mathbf F } _ { \text{core,ev} } \approx \frac { \dot { \mathbf F } _ { \text{core} } } { \lambda _ { \text{ch} } } \approx \frac { \dot { \mathbf E } _ { \text{core} } } { \lambda _ { \text{ch} } ^ { 2 } } \quad ( 1 7 ) \quad T _ { \text{visc} }$$

and

$$\dot { E } _ { \text{core,ev,rr} } = \dot { E } _ { \text{core,ev,\varphi\varphi} } = \frac { \dot { E } _ { \text{core,\varphi\varphi} } } { \lambda _ { \text{ch} } ^ { 2 } }. \quad ( 1 8 ) \quad \text{reve}$$

The radial and tangential stresses are related by the momentum balance as

$$\sigma _ { \text{shell}, \varphi \varphi } = - \alpha \lambda _ { \text{ch} } ^ { 3 } \sigma _ { \text{shell}, \text{rr} } \quad \quad ( 1 9 )$$

with the parameter α = 1 2 ( R core L shell -1 ) defined by the core radius R core and the shell thickness L shell .

Using equations (13), (15), (16), and (19), we find the differential equation for the radial stress component

$$\frac { d \sigma _ { e v } } { d t } = \frac { E _ { \text{core} } \cdot \frac { \dot { \mathbf E } _ { \text{core}, \varphi \varphi } } { \lambda _ { \text{ch} } ^ { 2 } } } { \lambda _ { \text{ch} } ^ { 2 } } \quad \quad ( 2 0 )$$

$$= \frac { E _ { \text{core} } } { \tau \lambda _ { \text{ch} } ^ { 2 } } \sinh \overset { \dots } { \longrightarrow } \overset { \sigma _ { \text{shell}, \text{visc}, \varphi \varphi } } { \sigma _ { \text{ref} } } \right ) \quad \ \ ( 2 1 )$$

$$= - \frac { \stackrel { \dots } { E _ { \text{core} } } } { \tau \lambda _ { \text{ch} } ^ { 2 } } \sinh \left ( \frac { \alpha \lambda _ { \text{ch} } ^ { 3 } \sigma _ { \text{ev} } } { \sigma _ { \text{ref} } } \right ) ^ { \cdot }. \quad \ \ ( 2 2 )$$

We solve the simplified differential equation in Eq. (22) analytically to describe the whole time dependence with a single analytical solution

$$\sigma _ { \text{ev} } = \frac { 2 \sigma _ { \text{ref} } } { \alpha \lambda _ { \text{ch} } ^ { 3 } } \cdot \text{atanh} \left ( C \cdot \exp \left ( - \frac { E _ { \text{core} } \alpha \lambda _ { \text{ch} } } { \tau \sigma _ { \text{ref} } } t \right ) \right ), \ \ ( 2 3 )$$

where the constant C can be determined from the boundary condition at time t = 0 with σ ev ( t = 0) = σ 0 .

For the calculation of the stress effect on the silicon voltage according to Eq. (5), we approximate the deformation of the silicon particle core as purely chemical, F core = F core ch , = λ ch Id , and we assume isotropic stress distribution inside the particle P core ev , = λ 2 ch σ ev Id . Therefore, the impact of the stress on the voltage according to Eq. (5) simplifies to ∆ U ev = v Li σ ev /F in the reduced model and the voltage relaxation reads

$$( 1 4 ) \quad \Delta U _ { \text{ev} } = \frac { 2 v _ { \text{Li} } \sigma _ { \text{ref} } } { \alpha F \lambda _ { \text{ch} } ^ { 3 } } \, \text{atanh} \left ( C \exp \left ( - \frac { E _ { \text{core} } \alpha \lambda _ { \text{ch} } } { \tau \sigma _ { \text{ref} } } t \right ) \right ). \ ( 2 4 )$$

To understand the origin and the regimes of the convoluted functional behavior in Eq. (23), we analyze the relaxation behavior in the limits of large and low stress magnitudes in section SI in the Supporting Information. Due to the importance of the long-term relaxation, here we present only the large stress limit. In the limit of large compressive stress, the differential equation (22) simplifies to

$$\cdot \quad \frac { \text{d} \sigma _ { \text{ev} } } { \text{d} t } = - \frac { E _ { \text{core} } } { \tau \lambda _ { \text{ch} } ^ { 2 } } \cdot \left ( - \frac { 1 } { 2 } \right ) \exp \left ( - \frac { \alpha \lambda _ { \text{ch} } ^ { 3 } \sigma _ { \text{ev} } } { \sigma _ { \text{ref} } } \right ). \quad ( 2 5 )$$

The analytical solution for this differential equation is

$$\sigma _ { \text{ev} } = \frac { \sigma _ { \text{ref} } } { \alpha \lambda _ { \text{ch} } ^ { 3 } } \cdot \ln \left ( \frac { E _ { \text{core} } \alpha \lambda _ { \text{ch} } } { 2 \tau \sigma _ { \text{ref} } } t + C _ { \text{exp} } \right ) \quad ( 2 6 )$$

with the integration constant C exp determined from the boundary condition.

Thus, the voltage relaxation according to the Garofalo viscosity

$$\Delta U _ { e v } = \frac { v _ { L i } \sigma _ { r e f } } { \alpha F \lambda _ { \text{ch} } ^ { 3 } } \cdot \ln \left ( \frac { E _ { \text{core} } \alpha \lambda _ { \text{ch} } } { 2 \tau \sigma _ { r e f } } t + C _ { \text{exp} } \right ) \quad ( 2 7 )$$

reveals logarithmic behavior in the large stress limit.

## C. Reduced Model Equations

Complementary to our full model presented in Section II A, we derive a reduced model with the key features in section SII in the Supporting Information. The reduced model describes the elastic stress contribution of the core at the interface between core and shell due to elastoplastic behavior of the shell σ ee and due to viscous behavior of the shell σ ev .

The system of equations defining the reduced chemomechanical hysteresis model reads

$$\frac { d \, S O C } { d t } = \frac { \dot { c } _ { \text{Li}, 0 } } { c _ { \text{Li}, \max } } = \pm \frac { C _ { \text{rate} } } { 3 6 0 0 } \, \frac { 1 } { s } \\ \hat { \ }.$$

$$\frac { \partial u } { d t } & = \begin{cases} - E _ { \text{shell} } \frac { 2 v _ { L i } ^ { 2 } } { 3 F \lambda _ { \text{ch} } ^ { 7 } } \dot { c } _ { L i, 0 }, & \text{if $f_{red}<0$} \\ \frac { \alpha \sigma _ { \text{Y,shell} } v _ { L i } ^ { 2 } } { F \left ( 1 + \alpha \lambda _ { \text{ch} } ^ { 3 } \right ) ^ { 2 } } \left | \dot { c } _ { L i, 0 } \right |, & \text{otherwise} \\ d \Delta U _ { \text{...} } & F _ { \text{...} } \dots v _ { r } \colon & / \alpha \lambda ^ { 3 }, \, F \Delta U _ { \text{...} } \dots \, & F _ { \text{...} } \dots v _ { r } ^ { 2 }. & \text{The d} \end{cases}$$

$$\frac { \text{d} \, \Delta U _ { e v } } { \text{d} t } = - \frac { \text{e} _ { \text{core} } v _ { \text{Li} } } { \tau F \lambda _ { \text{ch} } ^ { 2 } } \sinh \left ( \frac { \alpha \lambda _ { \text{ch} } ^ { 3 } F \Delta U _ { e v } } { \sigma _ { \text{ref} } v _ { \text{Li} } } \right ) - \frac { E _ { \text{core} } v _ { \text{Li} } ^ { 2 } } { 3 F \lambda _ { \text{ch} } ^ { 3 } } \dot { c } _ { \text{Li}, 0 } } \quad \text{uring $\mathfrak{s.d}} } { \text{The dash} }$$

with the parameter α = 1 2 ( R core L shell -1 ) and the yield condition for plastic flow for the reduced model

$$f _ { \text{red} } = - \text{sgn} \left ( \dot { c } _ { L i, 0 } \right ) \left ( 1 + \alpha \lambda _ { \text{ch} } ^ { 3 } \right ) \frac { F \Delta U _ { \text{ee} } } { v _ { L i } \sigma \text{Y,shell} } - 1 < 0. \ ( 3 1 ) \quad \text{silicor} \\ \quad \text{defor}.$$

The equations defining the reduced model describe the silicon anode voltage as U = U mean +∆ U ee +∆ U ev . Equation (28) states the change of SOC for lithiation (+) and delithiation ( -). The upper case in Eq. (29) describes the voltage evolution caused by elastic behavior of the silicon core due to elastic behavior of the shell. The lower case describes elastic core stress due to plastic behavior of the shell. The first term in Eq. (30) considers the viscous shell stress relaxation. The second term considers viscous shell stress increase because of silicon volume changes.

In Fig. 2, we depict the voltage profile predicted by the reduced model for a GITT procedure with (de)lithiation steps of ∆SOC = 0 02% with C 20 and relaxation pe-. / riods of 3 h. Furthermore, the figure shows the voltage during C 20 cycling and after / 12 h relaxation periods. The dashed black line depicts the fitted mean OCV curve U mean between the measured lithiation and delithiation voltage after 3 h rest period for a pure silicon anode from Ref. [8] used as true OCV curve for the simulations.

## III. COMPUTATIONAL &amp; EXPERIMENTAL DETAILS

## A. Simulation Setup

Our simulations describe the behavior of a silicon nanoparticle anode with a single-particle model. We im-

FIG. 2. Voltages according to the presented reduced model during GITT, C / 20 cycling, and after 12 h relaxation periods. The dashed black line depicts the mean OCV measured for a silicon anode in Ref. [8].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-18923585-f967-4676-b2e8-6da770167fb4.jpg)

plement our model in MATLAB using a finite-difference approach by discretizing the radial dimension. To solve the set of differential equations (6), (7), (8), (11), and (12), we use the solver ode15i. The variables inside the silicon core are the concentration of lithium c Li 0 , , the deformed radius r core , and the radial component of the plastic deformation F core pl rr , , of each silicon core element. The variables inside the SEI shell are the deformed radius r shell and the radial component of the plastic deformation F shell pl rr , , of each SEI shell element.

## B. Material Parameters

We adopt the parameters from our previous publication [21] and adapt where necessary. Particularly, we consider the stiff, inorganic SEI shell with Young's modulus of E shell = 100GPa compatible with experiments [44, 45]. The viscosity of the inner SEI shell is considered as a fit value and may range from η = 10 7 Pas for a highly viscous polymer [46] to η = 10 15 Pas for silicon oxide [47-49].

## C. Experimental Setup

The experiments analyzed in this study have been performed and published by Wycisk et al. [11] at Mercedes following discussions with the authors of this manuscript. The publication discusses full-cell voltage measurements with an NMC811 cathode and anodes with varying contents of silicon active material. Here, we constrain solely to the experimental results discussing anodes with pure silicon active material. The silicon anode consists of silicon nanoparticles attached to a conductive carbon net-

FIG. 3. (a) Experimental voltage relaxation of silicon at SOC = 0 3 over 300 hours after a charge and discharge period (protocol . SIII A) [11]. (b) The semi-logarithmic plot unveils the logarithmic voltage relaxation behavior.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-ad1db87c-5ad9-4c32-b717-726baef6e4bc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-071de362-b095-4d4a-a329-49e3a4d5572d.jpg)

work discussed as 'silicon-carbon composite granules' in Ref. [50]. We summarize the experimental and our simulation protocols in section SIII in the Supporting Information but refer to the experimental publication for the experimental details [11].

Throughout this manuscript, we consider voltages from the anode perspective and calculate voltage differences to the mean OCV, U -U mean . For comparison, the voltage difference for the full-cell measurements is calculated as U -U mean = -( U full -U full mean , ).

## IV. RESULTS AND DISCUSSION

## A. Experimental Results: Logarithmic Voltage Relaxation

over time as a semi-logarithmic plot. Apparently, the voltage relaxation profile does not follow a typical exponential relaxation behavior as illustrated in yellow. We identify a linear regime in the semi-logarithmic plot and fit a logarithmic function to the experimental data. The logarithmic fit agrees with the experimental data in a wide range of times t &lt; 20 h. Only for times larger than 20 h, the voltage relaxation slightly diminishes leaving the logarithmic regime. This is expected as logarithmic behavior would diverge for large times. The logarithmic voltage relaxation found in the experiment agrees with the experimentally observed voltage relaxation of silicon thin-film electrodes in Ref. [20]. It is in stark contrast to diffusion limitation with exponential long-term behavior and supports the theory of a mechanical origin.

First, we analyze the long-time relaxation experiment performed by Wycisk et al. [11] following the protocol described in section SIII A in the Supporting Information. In Fig. 3, we depict the voltage relaxation at the same SOC measured once in charge and once in discharge direction.

Interestingly, the authors of Ref. [11] find that even after 300 h of rest, the voltage depicted in Fig. 3(a) is not completely relaxed. Therefore, the true OCV value deviates from the relaxed voltage after 300 h and strongly deviates from standard GITT measurements with only a few hours of voltage relaxation. The authors of Ref. [11] exclude degradation or self-discharge due to the similar voltage relaxation profiles after the charge and discharge period. However, the mean value of the relaxed voltage after 300 h varies from the mean OCV measured with GITT for C 20 and 12 h rest period due to deviations of / the experimental cells.

Here, we investigate the voltage relaxation profile in detail again. In Fig. 3(b), we show the voltage relaxation

## B. Simulation Results: Slow Voltage Relaxation

As discussed in Section II A, the silicon OCV hysteresis results from elastoplastic stress generated by the shell, and the enlarged voltage hysteresis during cycling results from viscous shell stress acting on the particle core [21]. A simple Newtonian viscosity model, σ shell = η shell ˙ E shell , with constant viscosity η shell would imply exponential voltage relaxation behavior during rest contrasting the experimental observations. Due to the large stresses inside the SEI shell, the Newtonian model is not suitable for describing the viscous behavior. Instead, for large stresses, the strain rate is known to depend exponentially on the stress, leading to a logarithmic stress relaxation behavior. Therefore, we use the established Garofalo law given in Eq. (10) to describe both regimes.

Using the Garofalo model, Fig. 4 depicts our simulation results in comparison to the experimental data. The parameters are given in the Supporting Information in Table S1. We shift our simulations to match the observed

FIG. 4. (a) Voltage relaxation of silicon at SOC = 0 3 over 300 hours (protocol SIII A). Comparison of simulation, experiment . [11], and the analytical Garofalo approximation. (b) The semi-logarithmic plot shows agreement of the various curves.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-cfcd8581-b599-4a81-8916-7ed26483aa49.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-05597d5d-58c9-4d0c-bc5c-92bbccb1b3e0.jpg)

voltage after relaxation. The simulations reproduce the voltage relaxation profiles after the charge and discharge period. In particular, the simulation using the Garofalo law describes both the logarithmic relaxation regime as well as the decreasing relaxation after 20 h. The agreement confirms the explanation of the silicon voltage hysteresis by a visco-elastoplastic SEI behavior.

sulting from the elastoplastic contribution. Furthermore, the simulation shows an enlarged hysteresis during cycling caused by viscous stress.

To validate our simulation results, Fig. 4 compares our simulation and the experiment to the analytical approximation presented in Section II B. The analytical approximation for the voltage relaxation with Garofalo law viscosity reveals a similar logarithmic relaxation regime followed by a slowed relaxation. Thus, the specific trends observed for our simulation and the analytical approximation agree while the actual values deviate slightly. Nevertheless, as the analytical approach relies on several assumptions and approximations, the similarity of the voltage profile supports our simulation results.

## C. OCV and Cycling Voltage Hysteresis

Silicon anodes are generally known to show a significant voltage hysteresis. In Fig. 5, we depict the experimental OCV hysteresis after long relaxation and the enlarged voltage hysteresis during slow cycling [11]. We describe the protocol in SIII B in the Supporting Information. To check the consistency of our model with the experimental voltage hysteresis, Fig. 5 shows the simulation of the anode voltage during slow cycling and the OCV after long relaxation depending on the SOC for the parameters obtained from the voltage relaxation behavior. The illustrated voltages describe the influence on the silicon anode voltage. Hence, the voltage decreases during lithiation due to compressive stress and increases during delithiation due to tensile stress. The simulation results in Fig. 5 reveal a significant OCV hysteresis re-

The comparison of the cycling and relaxed voltages reveals a good agreement between simulation and experiment in a wide SOC regime. However, our simulation and the experiment deviate slightly at both extremes, SOC &lt; 0 2 and SOC . &gt; 0 8. . This disagreement results at least partially from the determination of the true, stress-free OCV curve as the mean between lithiation and delithiation OCV. At very high SOC, the compressive stress during lithiation is fully developed, while the tensile stress during the following delithiation has to build up gradually after the change of direction. Analogously, the tensile stress during delithiation is fully developed, while the compressive stress during the following lithiation has to build up gradually after the change of direction at low SOC. Therefore, the mean value between the lithiation and delithiation OCV at both extremes is not stress-free. Its consideration as true, stress-free OCV in the simulation leads to an apparent deviation. In the Supporting Information, we discuss a corrected OCV curve assuming a constant hysteresis size in the extreme SOC regimes. Fig. S3 reveals a better agreement between simulation and experiment in the extreme SOC regimes.

In our previous publication [21], we compared our simulation to the GITT measurement performed for a silicon half cell by Pan et al. [8, 9]. The cells differ significantly from the cells investigated by Wycisk et al. [11] due to a presumably different silicon raw material and electrolyte composition. Nevertheless, we compare our new model and the parameters obtained from the voltage relaxation [11] to the GITT measurement [8, 9] in section SV in the Supporting Information. Fig. S4 shows a reasonable match of simulation and experiment considering the full GITT procedure as well as a single GITT pulse. The agreement confirms the applicability of our chemo-

FIG. 5. C 20 and open-circuit voltage hysteresis after 12 h / relaxation in simulation and experiment (protocol SIII B) [11].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-e5459cef-2655-47f1-818e-02049b6d8351.jpg)

mechanical model to GITT measurements with different cells.

## D. C-Rate Dependence of Voltage Hysteresis

The experimental data obtained by Wycisk et al. [11] also cover the C-rate dependence of the voltage difference between the cycling voltage and the relaxed voltage after 12 h at SOC = 0 5 following the protocol given in . SIII C in the Supporting Information. As displayed in Fig. 6, the data reveal a linear dependence of the voltage on the C-rate. However, extrapolating this linear dependence to zero current results in a significant voltage offset compared to the OCV after infinite relaxation time. This offset would imply an enlarged hysteresis even for infinitely slow cycling, which is unexpected. Therefore, the authors conclude that the voltage will depart from the linear trend at particularly low C-rates.

The Newtonian viscosity model has a linear relation between the strain rate and the viscous stress. Hence, the size of the additional voltage hysteresis is linearly dependent on the C-rate as illustrated in yellow in Fig. 6. However, the Newtonian model explains no voltage offset, and the slope disagrees with the experiment when matching the hysteresis size at C / 10.

In comparison to the experimental and the Newtonian C-rate dependence, Fig. 6 also depicts the simulated Crate dependence. The inverse hyperbolic sine in Eq. (10) determines the C-rate dependence of the viscous stress and, consequently, the C-rate dependence of the additional voltage hysteresis during cycling. Thus, the simulation reveals a non-linear dependence of the voltage on the current. Nonetheless, after a swift increase of the voltage at current rates smaller C / 100, the increase slows down, approaching an almost linear trend with small curvature. Although the three experimental data points follow the linear trend exactly, we assume that our simula-

FIG. 6. C-rate dependence of voltage hysteresis at SOC = 0 5 . in simulation and experiment (protocol SIII C) [11].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-a5dcdb24-0d8a-4c76-a250-9a6eff4f7189.jpg)

tion is in reasonable agreement with the experiment and additionally describes the transition to vanishing voltage at zero current. We expect that more experimental data points particularly at low C-rates might indicate a curvature and deviation from the linear trend.

## E. Voltage Transition Profiles

Another interesting behavior is the silicon anode voltage profile of transitions between cycling and rest periods. In the following, we discuss the features of different transitions and compare our simulation to the experimental data from Ref. [11] wherever possible.

First, we investigate the transition profile between lithiation and delithiation according to protocol SIII D in the Supporting Information. In Fig. 7, we show the delithiation with either C / 10, C / 20, or GITT procedure after a continuous lithiation and rest period. For reference, the figure also includes the simulated and measured lithiation and delithiation OCV curves from Fig. 5, which almost coincide in the depicted regime 0 3 . &lt; SOC &lt; . 0 5. All experimental data [11] reveal a smooth transition between the lithiation and delithiation voltage. The slope of the voltage profiles is large directly after the change of direction and slows down gradually when approaching the delithiation voltage.

The numerical results are depicted in Fig. 7 compared to the experiment. When switching the current direction from lithiation to delithiation, the simulated voltage profiles for C / 10 (yellow) and C / 20 (purple) currents reveal three regimes. Immediately after the change of direction, the voltage shows a steep increase for a small span of ∆SOC ≈ 0 01 attributed to the rapid build-up of vis-. cous stress. Afterward, for a range of ∆SOC ≈ 0 1, . a constant, moderate voltage slope demonstrates the decrease of compressive elastic stress and the subsequent

FIG. 7. Voltage transition from lithiation to delithiation in simulation and experiment (protocol SIII D) [11].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-3a6f5a1c-321a-4347-96ca-c59faee3b45f.jpg)

increase in tensile elastic stress. In the third regime, the slope slows down, and the voltage approaches a maximum value when reaching the yield criterion for plasticity. The higher current C / 10 shows a slightly faster voltage transition compared to the lower current C / 20. For the GITT transition curve (green), the relaxation of viscous stress during the rest periods suppresses the viscous regime after the change of direction. Contrary to the simulation, the experimental curves do not reveal clearly defined regimes but are in line with the general trend of a rapid voltage increase after the change of direction followed by an attenuated transition to the delithiation voltage curve. The much smoother experimental results compared to our simulation are expected as we consider only a single-particle model but the detailed features average out in the multi-particle experiment. Thus, we conclude that our simulation result agrees reasonably with the experimental measurement.

In the Supporting Information in section SVI, we evaluate the behavior of an interrupted lithiation pulse for different C-rates and at different SOC values. All voltage profiles in Figures S5 and S6 show a steep slope at the beginning of the pulses, revealing the increase in viscous stress followed by a slower convergence to the lithiation voltage, indicating elastoplastic behavior. The similar voltage profiles for different C-rates indicate that the voltage transition needs a certain amount of charge throughput or SOC change ∆SOC in accordance with the experimental results from Ref. [11] for a blended graphite-silicon anode. Additionally, the voltage profiles at different SOC values in Fig. S6 show that the general trends of the chemo-mechanical simulation agree with the ones of the experiment. However, all experimental curves show an overshoot instead of a smooth convergence to the lithiation voltage, which is not visible in our simulations. In terms of mechanics, this overshoot might result from a thixotropic behavior of the SEI shell as discussed in the Supporting Information.

Another voltage hysteresis effect measured for silicon anodes is a pronounced relaxation during rest observed for higher applied currents [11, 51]. Higher C-rates show an increased voltage hysteresis during cycling in agreement with our viscosity model. However, this dependence surprisingly inverts after relaxation. This phenomenon is not captured in our chemo-mechanical single-particle model. Therefore, we support the interpretation as a multi-particle effect [11] and add a mechanical explanation. For fast charging, the silicon particles inside the anode will lithiate more inhomogeneously, causing enhanced plastic flow of the shell around particles with a higher lithiation level. During the subsequent rest period, the silicon particles with initially higher lithiation degrees delithiate slightly. The shrinkage of those particles reduces the remaining compressive stress, while the stress in the particles with initially lower lithiation levels can not exceed the yield stress for plastic flow. Hence, this multi-particle effect can reduce the mean stress hysteresis inside the silicon anode and, consequently, the voltage hysteresis after relaxation.

Finally, we estimate the voltage transition behavior for alternating short lithiation and delithiation pulses following protocol SIII F in the Supporting Information. The silicon voltage hysteresis is often described empirically with the Plett model presented in section SVII in the Supporting Information [52-54]. In Fig. 8(a), we depict the behavior for alternating pulses with ∆SOC = 0 01 . predicted by the empirical Plett model with the parameters adjusted to fit the experimental voltage hysteresis. The Plett model does not reveal a constant hysteresis behavior during 10 subsequent cycles but rather approaches the mean OCV within the first cycles and then describes a hysteresis around it. Additionally, the Plett model is not able to account for a relaxation phase without a change in SOC. In contrast, Fig. 8(b) shows the simulation of alternating pulses, which reveal a permanent hysteresis during 10 subsequent cycles. Only the very first pulse initially shows a slightly different behavior with an enlarged hysteresis size because of a different stress state in the initial situation after the 12 h relaxation period. We know that experiments show a permanent hysteresis behavior upon alternate lithiation and delithiation pulses in line with our simulation result. Thus, we conclude that our chemo-mechanical core-shell model outperforms the empirical Plett model in the description of voltage hysteresis phenomena.

## V. CONCLUSIONS

Detailed analysis of the silicon voltage hysteresis experiments performed by Wycisk et al. [11] reveals a slow, non-exponential voltage relaxation. We identify a logarithmic voltage relaxation for a wide range of times and a transition to exponential relaxation for larger times due to the divergence of the logarithmic behavior. With a chemo-mechanical core-shell model, we have illustrated

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-348aed68-2be7-411e-a37b-340d024c36da.jpg)

FIG. 8. Voltage for alternating lithiation and delithiation pulses with ∆SOC = 1% (protocol SIII F) for (a) the phenomenological Plett model and (b) our chemo-mechanical simulation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01106v1-figure-6f9c827a-a0a6-43cc-9318-985a423c1801.jpg)

that the visco-elastoplastic SEI shell behavior following the Garofalo law or inverse hyperbolic sine law for viscosity can accurately describe the voltage relaxation of a silicon anode over the whole time span. Our simulations also reproduce the observed voltage hysteresis and GITT measurement with the parameters obtained from the relaxation experiment.

Plett model regarding physical understanding as well as the description of the various features of the hysteresis phenomenon.

Additionally, the Garofalo viscosity model can approach the experimentally observed C-rate dependence of the cycling voltage hysteresis. The inverse hyperbolic sine behaves approximately linear in a wide span of Crates but shows a kink and reveals vanishing additional voltage hysteresis at zero current. Therefore, the Garofalo law viscosity model fits much better to the C-rate dependence than Newtonian viscosity, which reveals a proportional relation between the voltage and the applied C-rate.

With a focus on the voltage transition behavior between lithiation and delithiation, the presented chemomechanical model can adequately describe the general trends of an initially fast voltage transition followed by an attenuated convergence to the delithiation voltage curve. The interplay of viscous, elastic, and plastic contributions to the simulated voltage explains this voltage profile. Furthermore, our model reasonably describes the lithiation behavior after a rest period. Thus, our chemomechanical core-shell model outperforms the empirical

[1] L. Sun, Y. Liu, R. Shao, J. Wu, R. Jiang, and Z. Jin, Recent progress and future perspective on practical silicon anode-based lithium ion batteries, Energy Storage Materials 46 , 482 (2022).

[2] X. Zuo, J. Zhu, P. M¨ller-Buschbaum, and Y.-J. Cheng, u Silicon based lithium-ion battery anodes: A chronicle perspective review, Nano Energy 31 , 113 (2017).

[3] K. Feng, M. Li, W. Liu, A. G. Kashkooli, X. Xiao, M. Cai,

The overall accordance of our simulations to experimental results supports our chemo-mechanical explanation of the voltage hysteresis presented initially in Ref. [21]. The description of the viscous behavior using the Garofalo law is more suitable than linear Newtonian viscosity because of the large stresses reached inside the SEI shell. In conclusion, we have demonstrated that our physical model presents a consistent picture of the various features of the silicon voltage hysteresis phenomenon.

## ACKNOWLEDGMENTS

Lukas K¨ obbing gratefully acknowledges funding and support by the European Union's Horizon Europe within the research initiative Battery 2030+ via the OPINCHARGE project under the grant agreement number 101104032 and by the German Research Foundation (DFG) within the research training group SiMET under project number 281041241/GRK2218. The authors appreciate fruitful discussions with Dominik Wycisk, Felix Schwab, Martin Werres, Raphael Schoof, and Arnulf Latz.

and Z. Chen, Silicon-Based Anodes for Lithium-Ion Batteries: From Fundamentals to Practical Applications, Small 14 , 10.1002/smll.201702737 (2018).

[4] L. Y. Beaulieu, K. W. Eberman, R. L. Turner, L. J. Krause, and J. R. Dahn, Colossal Reversible Volume Changes in Lithium Alloys, Electrochemical and SolidState Letters 4 , A137 (2001).

[5] X. H. Liu, L. Zhong, S. Huang, S. X. Mao, T. Zhu,

- and J. Y. Huang, Size-Dependent Fracture of Silicon Nanoparticles During Lithiation, ACS Nano 6 , 1522 (2012).
- [6] M. Wetjen, S. Solchenbach, D. Pritzl, J. Hou, V. Tileli, and H. A. Gasteiger, Morphological Changes of Silicon Nanoparticles and the Influence of Cutoff Potentials in Silicon-Graphite Electrodes, Journal of The Electrochemical Society 165 , A1503 (2018).
- [7] F. Kilchert, M. Schammer, A. Latz, and B. Horstmann, Silicon Nanowires as Anodes for Lithium-Ion Batteries: Full Cell Modeling, Energy Technology 2400206 , 10.1002/ente.202400206 (2024), 2401.16125.
- [8] K. Pan, F. Zou, M. Canova, Y. Zhu, and J.-H. Kim, Systematic electrochemical characterizations of Si and SiO anodes for high-capacity Li-Ion batteries, Journal of Power Sources 413 , 20 (2019).
- [9] K. Pan, Dissertation , Ph.D. thesis, Ohio State University (2020).
- [10] P. Bernard, J. P. Alper, C. Haon, N. Herlin-Boime, and M. Chandesris, Electrochemical analysis of silicon nanoparticle lithiation - Effect of crystallinity and carbon coating quantity, Journal of Power Sources 435 , 226769 (2019).
- [11] D. Wycisk, G. K. Mertin, M. Oldenburger, O. von Kessel, and A. Latz, Challenges of open-circuit voltage measurements for silicon-containing Li-Ion cells, Journal of Energy Storage 89 , 111617 (2024).
- [12] M. T. McDowell, S. W. Lee, W. D. Nix, and Y. Cui, 25th Anniversary Article: Understanding the Lithiation of Silicon and Other Alloying Anodes for Lithium-Ion Batteries, Advanced Materials 25 , 4966 (2013).
- [13] D. Wycisk, G. K. Mertin, M. Oldenburger, and A. Latz, Analysis of heat generation due to open-circuit voltage hysteresis in lithium-ion cells, Journal of Energy Storage 61 , 106817 (2023).
- [14] V. A. Sethuraman, M. J. Chon, M. Shimshak, V. Srinivasan, and P. R. Guduru, In situ measurements of stress evolution in silicon thin films during electrochemical lithiation and delithiation, Journal of Power Sources 195 , 5062 (2010).
- [15] B. Lu, Y. Song, Q. Zhang, J. Pan, Y.-T. Cheng, and J. Zhang, Voltage hysteresis of lithium ion batteries caused by mechanical stress, Physical Chemistry Chemical Physics 18 , 4721 (2016).
- [16] R. Chandrasekaran, A. Magasinski, G. Yushin, and T. F. Fuller, Analysis of Lithium Insertion/Deinsertion in a Silicon Electrode Particle at Room Temperature, Journal of The Electrochemical Society 157 , A1139 (2010).
- [17] K. Zhao, M. Pharr, S. Cai, J. J. Vlassak, and Z. Suo, Large Plastic Deformation in High-Capacity Lithium-Ion Batteries Caused by Charge and Discharge, Journal of the American Ceramic Society 94 , s226 (2011).
- [18] Z. Cui, F. Gao, and J. Qu, A finite deformation stressdependent chemical potential and its applications to lithium ion batteries, Journal of the Mechanics and Physics of Solids 60 , 1280 (2012).
- [19] M. T. McDowell, S. W. Lee, J. T. Harris, B. A. Korgel, C. Wang, W. D. Nix, and Y. Cui, In Situ TEM of TwoPhase Lithiation of Amorphous Silicon Nanospheres, Nano Letters 13 , 758 (2013).
- [20] V. A. Sethuraman, V. Srinivasan, and J. Newman, Analysis of Electrochemical Lithiation and Delithiation Kinetics in Silicon, Journal of The Electrochemical Society 160 , A394 (2013).
- [21] L. K¨ obbing, A. Latz, and B. Horstmann, Voltage Hysteresis of Silicon Nanoparticles: Chemo-Mechanical Particle-SEI Model, Advanced Functional Materials 34 , 2308818 (2024).
- [22] B. Horstmann, F. Single, and A. Latz, Review on multiscale models of solid-electrolyte interphase formation, Current Opinion in Electrochemistry 13 , 61 (2019).
- [23] X. Zhang, S. Weng, G. Yang, Y. Li, H. Li, D. Su, L. Gu, Z. Wang, X. Wang, and L. Chen, Interplay between solidelectrolyte interphase and (in)active LixSi in silicon anode, Cell Reports Physical Science 2 , 100668 (2021).
- [24] M. Nie, D. P. Abraham, Y. Chen, A. Bose, and B. L. Lucht, Silicon Solid Electrolyte Interphase (SEI) of Lithium Ion Battery Characterized by Microscopy and Spectroscopy, The Journal of Physical Chemistry C 117 , 13403 (2013).
- [25] P. Verma, P. Maire, and P. Nov´k, A review of the feaa tures and analyses of the solid electrolyte interphase in Li-ion batteries, Electrochimica Acta 55 , 6332 (2010).
- [26] E. Peled and S. Menkin, Review-SEI: Past, Present and Future, Journal of The Electrochemical Society 164 , A1703 (2017).
- [27] L. K¨ obbing, A. Latz, and B. Horstmann, Growth of the solid-electrolyte interphase: Electron diffusion versus solvent diffusion, Journal of Power Sources 561 , 232651 (2023).
- [28] L. von Kolzenberg, A. Latz, and B. Horstmann, Solid-Electrolyte Interphase During Battery Cycling: Theory of Growth Regimes, ChemSusChem 13 , 3901 (2020).
- [29] L. von Kolzenberg, A. Latz, and B. Horstmann, ChemoMechanical Model of SEI Growth on Silicon Electrode Particles, Batteries &amp; Supercaps 5 , 1 (2022).
- [30] K. Guo, R. Kumar, X. Xiao, B. W. Sheldon, and H. Gao, Failure progression in the solid electrolyte interphase (SEI) on silicon electrodes, Nano Energy 68 , 104257 (2020).
- [31] J. Lee, J.-y. Jeong, J. Ha, Y.-t. Kim, and J. Choi, Understanding solid electrolyte interface formation on graphite and silicon anodes in lithium-ion batteries: Exploring the role of fluoroethylene carbonate, Electrochemistry Communications 163 , 107708 (2024).
- [32] C. Cao, I. I. Abate, E. Sivonxay, B. Shyam, C. Jia, B. Moritz, T. P. Devereaux, K. A. Persson, H.-G. Steinr¨ck, and M. F. Toney, Solid Electrolyte Interphase u on Native Oxide-Terminated Silicon Anodes for Li-Ion Batteries, Joule 3 , 762 (2019).
- [33] Q. Li, X. Liu, X. Han, Y. Xiang, G. Zhong, J. Wang, B. Zheng, J. Zhou, and Y. Yang, Identification of the Solid Electrolyte Interface on the Si/C Composite Anode with FEC as the Additive, ACS Applied Materials &amp; Interfaces 11 , 14066 (2019).
- [34] J. Chen, X. Fan, Q. Li, H. Yang, M. R. Khoshi, Y. Xu, S. Hwang, L. Chen, X. Ji, C. Yang, H. He, C. Wang, E. Garfunkel, D. Su, O. Borodin, and C. Wang, Electrolyte design for LiF-rich solid-electrolyte interfaces to enable high-performance microsized alloy anodes for batteries, Nature Energy 5 , 386 (2020).
- [35] Y. He, L. Jiang, T. Chen, Y. Xu, H. Jia, R. Yi, D. Xue, M. Song, A. Genc, C. Bouchet-Marquis, L. Pullan, T. Tessner, J. Yoo, X. Li, J.-G. Zhang, S. Zhang, and C. Wang, Progressive growth of the solid-electrolyte interphase towards the Si anode interior causes capacity fading, Nature Nanotechnology 16 , 1113 (2021).

- [36] M. Schnabel, S. P. Harvey, E. Arca, C. Stetson, G. Teeter, C. Ban, and P. Stradins, Surface SiO 2 Thickness Controls Uniform-to-Localized Transition in Lithiation of Silicon Anodes for Lithium-Ion Batteries, ACS Applied Materials &amp; Interfaces 12 , 27017 (2020).
- [37] K. W. Schroder, A. G. Dylla, S. J. Harris, L. J. Webb, and K. J. Stevenson, Role of Surface Oxides in the Formation of Solid-Electrolyte Interphases at Silicon Electrodes for Lithium-Ion Batteries, ACS Applied Materials &amp; Interfaces 6 , 21510 (2014).
- [38] V. L. Deringer, N. Bernstein, G. Cs´nyi, C. Ben Maha moud, M. Ceriotti, M. Wilson, D. A. Drabold, and S. R. Elliott, Origins of structural and electronic transitions in disordered silicon, Nature 589 , 59 (2021).
- [39] V. L. Chevrier, L. Liu, D. B. Le, J. Lund, B. Molla, K. Reimer, L. J. Krause, L. D. Jensen, E. Figgemeier, and K. W. Eberman, Evaluating Si-Based Materials for Li-Ion Batteries in Commercially Relevant Negative Electrodes, Journal of The Electrochemical Society 161 , A783 (2014).
- [40] K. Kitada, O. Pecher, P. C. M. M. Magusin, M. F. Groh, R. S. Weatherup, and C. P. Grey, Unraveling the Reaction Mechanisms of SiO Anodes for Li-Ion Batteries by Combining in Situ 7 Li and ex Situ 7 Li/ 29 Si Solid-State NMR Spectroscopy, Journal of the American Chemical Society 141 , 7014 (2019).
- [41] J. Wang, X. Wang, B. Liu, H. Lu, G. Chu, J. Liu, Y.G. Guo, X. Yu, F. Luo, Y. Ren, L. Chen, and H. Li, Size effect on the growth and pulverization behavior of Si nanodomains in SiO anode, Nano Energy 78 , 105101 (2020).
- [42] F. Garofalo, O. Richmond, W. F. Domis, and F. von Gemmingen, Strain-Time, Rate-Stress and RateTemperature Relations during Large Deformations in Creep, Proceedings of the Institution of Mechanical Engineers, Conference Proceedings 178 , 31 (1963).
- [43] E. T. Stang, Constitutive Modeling of Creep in Leaded and Lead-Free Solder Alloys Using Constant Strain-Rate Tensile Testing , Master of science in mechanical engineering, Wright State University (2018).
- [44] H. Shin, J. Park, S. Han, A. M. Sastry, and W. Lu,
- Component-/structure-dependent elasticity of solid electrolyte interphase layer in Li-ion batteries: Experimental and computational studies, Journal of Power Sources 277 , 169 (2015).
- [45] Y. Chai, W. Jia, Z. Hu, S. Jin, H. Jin, H. Ju, X. Yan, H. Ji, and L.-J. Wan, Monitoring the mechanical properties of the solid electrolyte interphase (SEI) using electrochemical quartz crystal microbalance with dissipation, Chinese Chemical Letters 32 , 1139 (2021).
- [46] R. Edgeworth, B. J. Dalton, and T. Parnell, The pitch drop experiment, European Journal of Physics 5 , 198 (1984).
- [47] P. Sutardja and W. Oldham, Modeling of stress effects in silicon oxidation, IEEE Transactions on Electron Devices 36 , 2415 (1989).
- [48] V. Senez, D. Collard, B. Baccus, M. Brault, and J. Lebailly, Analysis and application of a viscoelastic model for silicon oxidation, Journal of Applied Physics 76 , 3285 (1994).
- [49] M. I. Ojovan, Viscosity and Glass Transition in Amorphous Oxides, Advances in Condensed Matter Physics 2008 , 1 (2008).
- [50] J. Schwan, G. Nava, and L. Mangolini, Critical barriers to the large scale commercialization of silicon-containing batteries, Nanoscale Advances 2 , 4368 (2020).
- [51] A. Durdel, S. Friedrich, L. H¨sken, and A. Jossen, Modelu ing Silicon-Dominant Anodes: Parametrization, Discussion, and Validation of a Newman-Type Model, Batteries 9 , 558 (2023).
- [52] G. L. Plett, Extended Kalman filtering for battery management systems of LiPB-based HEV battery packs, Journal of Power Sources 134 , 262 (2004).
- [53] C. P. Graells, M. S. Trimboli, and G. L. Plett, Differential hysteresis models for a silicon-anode Li-ion battery cell, in 2020 IEEE Transportation Electrification Conference &amp; Expo (ITEC) , Vol. 1 (IEEE, 2020) pp. 175-180.
- [54] D. Wycisk, M. Oldenburger, M. G. Stoye, T. Mrkonjic, and A. Latz, Modified Plett-model for modeling voltage hysteresis in lithium-ion cells, Journal of Energy Storage 52 , 105016 (2022).
Lukas Köbbing Yannick Kuhn Birger Horstmann
Slow Voltage Relaxation of Silicon Nanoparticles with a Chemo-Mechanical Core-Shell Model
Silicon presents itself as a high-capacity anode material for lithium-ion
batteries with a promising future. The high ability for lithiation comes along
with massive volume changes and a problematic voltage hysteresis, causing
reduced efficiency, detrimental heat generation, and a complicated
state-of-charge estimation. During slow cycling, amorphous silicon
nanoparticles show a larger voltage hysteresis than after relaxation periods.
Interestingly, the voltage relaxes for at least several days, which has not
been physically explained so far. We apply a chemo-mechanical continuum model
in a core-shell geometry interpreted as a silicon particle covered by the
solid-electrolyte interphase to account for the hysteresis phenomena. The
silicon core (de)lithiates during every cycle while the covering shell is
chemically inactive. The visco-elastoplastic behavior of the shell explains the
voltage hysteresis during cycling and after relaxation. We identify a
logarithmic voltage relaxation, which fits with the established Garofalo law
for viscosity. Our chemo-mechanical model describes the observed voltage
hysteresis phenomena and outperforms the empirical Plett model. In addition to
our full model, we present a reduced model to allow for easy voltage profile
estimations. The presented results support the mechanical explanation of the
silicon voltage hysteresis with a core-shell model and encourage further
efforts into the investigation of the silicon anode mechanics.