Verbeterd chemisch-potentieel model voor snelle deformatie in silicium nanodeeltje elektroden van Li-ionbatterijen

Samenvatting

We presenteren een verfijnde chemisch-potentieel expressie die de lage-deformatiesnelheid assumptie verwijdert, wat accurate modellering van diffusie-geïnduceerde spanning en lithiumdistributie in silicium nanodeeltje elektroden mogelijk maakt. Door het nieuwe model te vergelijken met het traditionele onder potentiostatische en galvanostatische cycli, tonen we een duidelijke divergentie tijdens potentiostatische operatie, terwijl galvanostatische resultaten grotendeels ongewijzigd blijven. Een kritieke straal—waar Li-flux en hydrostatische spanning essentieel identiek zijn voor beide modellen—wordt voor het eerst geïdentificeerd.

Inleiding

Hoog-capaciteit elektroden zoals silicium, die volumetrische veranderingen van tot 400% tijdens lithiatie ondergaan, zijn crucieel voor draagbare elektronica, elektrische voertuigen en net-schaal opslag. Echter, de diffusie-geïnduceerde spanning gegenereerd door homogene volumetrische expansie kan bros breuk activeren, wat cycluslevensduur verkort en capaciteit vermindert. Composiet elektroden vertonen complexe morfologieën die zuivere theoretische beschrijvingen uitdagen; de meeste modellen behandelen materiaaleigenschappen als ruimtelijk variërend terwijl interfaciale effecten genegeerd worden. Typische elektrodegeometrieën in theorie zijn sferisch, cylindrisch en planaar, met een- of tweedimensionale modellen dienovereenkomstig.

Recente onderzoeken op silicium nanodeeltje elektroden hebben de wisselwerking tussen lithiatie-geïnduceerde fasetransformatie, morfologische evolutie en spanningsgeneratie benadrukt. Deze studies vertrouwen op thermodynamisch evenwicht van vastestoffen, waar atomaire diffusie beïnvloed wordt door spanning, en vice versa. Larche en Cahn's framework introduceerde een spanning-afhankelijk chemisch potentieel aannemend kleine, isotrope deformaties, terwijl Wu's formulering de Eshelby-impulstensor incorporeerde. Cui et al. breidden dit uit tot eindige deformaties, doch hun herleiding blijft enkel accuraat voor lage deformatiesnelheden ten opzichte van diffusie. Gegeven silicium's grote volumetrische expansie tijdens snelle lading, kan deze assumptie tot significante fouten leiden.

In dit werk herleiden we een chemisch-potentieel expressie geldig voor snelle deformatie zonder de lage-snelheid assumptie, onafhankelijk van elektrodemorfologie. Door het nieuwe model toe te passen op silicium nanodeeltje elektroden onder zowel potentiostatische als galvanostatische controle, kwantificeren we de impact op spanning en lithiumdistributie. We identificeren ook een kritieke straal waar de twee modellen convergent, wat inzicht geeft in de ruimtelijke grenzen van spanningsinvloed.

Methoden

Mechanische vergelijkingen

Lithiuminvoeging induceert volumetrische verandering. We beschrijven deformatie met zowel Lagrangiaanse als Euleriaanse coördinaten: $$\mathbf{U}=\mathbf{x}-\mathbf{X}\quad (1)$$ waar U het verplaatsingsveld is. De deformatiegradiënttensor is $$\mathbf{F}=\frac{\partial \mathbf{x}}{\partial \mathbf{X}}=\mathbf{I}+\mathrm{Grad}\,\mathbf{U}\quad (2)$$

Voor een sferisch deeltje met coördinaten (R,Θ,Φ) (Lagrangiaans) en (r,θ,φ) (Euleriaans), vereenvoudigt de gradiënttensor tot $$\mathbf{F}=\begin{bmatrix}{F}_R&0&0\0&{F}_{\Theta}&0\0&0&{F}_{\Phi}\end{bmatrix}=\begin{bmatrix}1+\partial u/\partial R&0&0\0&1+u/R&0\0&0&1+u/R\end{bmatrix}\quad (3)$$

Deformatie wordt verdeeld in lithiatie-geïnduceerde vormverandering en reversibele elastische deformatie: $$\mathbf{F}=\mathbf{F}^e\mathbf{F}^c\quad (4)$$ met F^c isotrop: $$\mathbf{F}^c=(1+\Omega C)^{1/3}\mathbf{I}\quad (5)$$ waar Ω het partiële molaire volume is. Gevolglijk, $$\mathbf{F}^e=(1+\Omega(R)C)^{-1/3}\begin{bmatrix}1+\partial u/\partial R&0&0\0&1+u/R&0\0&0&1+u/R\end{bmatrix}\quad (6)$$

De Green-Lagrange spanningstensor E is $$\mathbf{E}=\frac{1}{2}\left(\mathbf{F}^T\mathbf{F}-\mathbf{I}\right)\quad (7)$$ wat radiale en tangentiële componenten oplevert $$E_R^e=\tfrac{1}{2}\Big[\tfrac{(1+\partial u/\partial R)^2}{(1+\Omega(R)C)^{2/3}}-1\Big]\quad (9)$$ $$E_{\Theta}^e=E_{\Phi}^e=\tfrac{1}{2}\Big[\tfrac{(1+u/R)^2}{(1+\Omega(R)C)^{2/3}}-1\Big]\quad (10)$$

Aannemende lineaire elasticiteit, is de spanningenergiedichtheid $$W=\tfrac{1}{2}\mathbf{E}^e:\mathbf{C}:\mathbf{E}^e=\det(\mathbf{F}^c)\frac{E_h}{2(1+\nu)}\Big[\tfrac{\nu}{1-2\nu}\big(tr\,\mathbf{E}^e\big)^2+tr\,\mathbf{E}^e\mathbf{E}^e\Big]\quad (12)$$ waar E_h en ν Young's modulus en Poisson's ratio zijn. De eerste Piola-Kirchhoff spanning volgt $$\mathbf{P}=\det(\mathbf{F}^c)\frac{E_h}{1+\nu}\Big[\tfrac{\nu}{1-2\nu}tr\,\mathbf{E}^e+\mathbf{E}^e\Big]\mathbf{F}^e(\mathbf{F}^c)^{-1}\quad (13)$$

Componentgewijs, $$P_R=\left(1+\Omega C\right)^{1/3}\frac{E_h}{2(1+\nu)(1-2\nu)}\left(1+\tfrac{\partial u}{\partial R}\right)\Big[(1-\nu)E_R^e+2\nu E_{\Theta}^e\Big]\quad (14)$$ $$P_{\Theta}=P_{\Phi}=\left(1+\Omega C\right)^{1/3}\frac{E_h}{2(1+\nu)(1-2\nu)}\left(1+\tfrac{u}{R}\right)\Big[\nu E_R^e+E_{\Theta}^e\Big]\quad (15)$$

Evenwicht zonder volumekrachten vereist $$\frac{\partial P_R}{\partial R}+\frac{2(P_R-P_{\Phi})}{R}=0\quad (16)$$ met randvoorwaarden $$u(0,t)=0,\;P_R(R_0,t)=0\quad (17)$$

Massa-transportvergelijking

Lithiumconcentratie C(X,t) en flux J(X,t) voldoen aan $$\frac{\partial C}{\partial t}+\mathrm{Div}\,\mathbf{J}=0\quad (18)$$. Onder sferische symmetrie reducert dit tot $$\frac{\partial C(R,t)}{\partial t}+\frac{1}{R^2}\frac{\partial\big(R^2J(R,t)\big)}{\partial R}=0\quad (19)$$. De flux volgt Fick's wet met chemisch-potentieel dryvende kracht: $$J=-\frac{CD}{R_gTF_{11}F_{11}}\frac{\partial \mu}{\partial R}\quad (20)$$.

Het chemisch potentieel afgeleid van totale interne energiedichtheid $$\Pi=\varphi(C)+W(C,\mathbf{E}^e)\quad (22)$$, geeft $$\mu=\frac{\partial\Pi}{\partial C}=\mu_0(C)+\tau(\mathbf{E}^e,C)\quad (23)$$ waar τ het spanning-afhankelijke deel is. Voor een verdunde oplossing, $$\mu_0(C)=\mu_0+R_gT\ln(\gamma C)\quad (24)$$ met γ=1.

Traditioneel levert de Helmholtz-vrije-energie benadering $$\tau_H=\frac{\partial W}{\partial C}= -\det(\mathbf{F}^e)\sigma_m\Omega\quad (25)$$ leidend tot $$\mu={\mu}_0+R_gT\ln C-\det(\mathbf{F}^e)\Omega\sigma_m\quad (26)$$. Hier is σ_m de Cauchy hydrostatische spanning $$\sigma_m=\tfrac{1}{3}tr\,(\det^{-1}\mathbf{F}\mathbf{PF}^T)\quad (27)$$. Deze expressie (26) is het traditionele chemisch-potentieel model.

Alternatief, het behandelen van de Gibbs-vrije-energie levert $$\tau_G=\frac{\partial W}{\partial C}\quad (28)$$ en gevolglik $$\mu={\mu}_0+R_gT\ln C-\frac{\partial W}{\partial C}\quad (29)$$, ons ontwikkelde model. Het massa-transportsysteem koppelt (19), (20), (26) en (29). Randvoorwaarden voor potentiostatische (C(R₀,t)=C_max) en galvanostatische (J(R₀,t)=j₀(1+u/R)^2) operaties zijn $$C(R_0,t)=C_{\max},\;J(R_0,t)=j_0\big(1+u/R\big)^2\quad (30-31)$$ met initiële $$C(R,0)=0\quad (32)$$.

Numerieke implementatie

Analytische oplossingen zijn intricater; we losten de gekoppelde PDE's op met COMSOL Multiphysics, materiaalparameters uit Tabel 1 aannemend. Dimensionloze variabelen werden gebruikt voor duidelijkheid. Laadtoestand (SOC) wordt berekend als $$SOC=\frac{\int_0^{R_0}C(R,t)R^2dR}{\int_0^{R_0}C_{\max}R^2dR}\quad (33)$$. Fluxen afgeleid van de twee chemisch-potentieel modellen zijn $$J_H=\frac{\partial\tau_H}{\partial R},\;J_G=\frac{\partial\tau_G}{\partial R}\quad (34)$$.

Resultaten en discussie

Figuur 1 illustreert lithiumconcentratie, radiale spanning en hoopspanning profielen voor een sferisch Si nanodeeltje onder galvanostatische cycli. Het ontwikkelde (volle lijnen) en traditionele (driehoeken) modellen leveren vrijwel identieke concentratiedistributies op over alle SOC's. Spanningprofielen verschillen slechts licht, met bescheiden centrale elevatie bij 46,7% en 65,5% SOC.

Figuur 2 presenteert analoge data voor potentiostatische cycli. Hier is de divergentie tussen de twee modellen duidelijker, met name in centrale spanningen bij intermediaire SOC's. Omdat potentiostatische operatie een constante oppervlakteconcentratie oplegt, is de laadsnelheid—en gevolglik de deformatiesnelheid—hoger dan in galvanostatische modus. De algehele deformatie voor een gegeven SOC blijft vergelijkbaar, wat benadrukt dat de sleutelfactor de deformatiesnelheid is en niet de absolute deformatie.

Figuur 3 (uit Ref. [17]) demonstreert dat een 620-nm Si nanodeeltje volledig kan expanderen binnen één minuut bij 2 V versus Li, wat de relevantie van snelle-deformatie modellering benadrukt. In dergelijke regimes voorspellen de twee chemisch-potentieel modellen merkbaar verschillende spanningsstromen, hoewel experimentele spanningdata onbeschikbaar blijft voor kwantitatieve validatie.

Figuur 4 onderzoekt de fluxratio J_H/J_G onder galvanostatische cycli bij diverse laadstromen. De ratio blijft dicht bij eenheid over het deeltje, met een lichte toename naar het centrum en afname naar het oppervlak. De spreiding verbreedt met hogere stroom, wat sterkere gevoeligheid voor deformatiesnelheid aangeeft. Alle krommen snijden zich nabij een kritieke straal—getermd de chemisch-potentieel-onafhankelijke regio (CPOR)—waar fluxes convergent zijn ongeacht het gekozen model. CPOR verschuift dichter bij het oppervlak naarmate de stroom toeneemt.

Figuren 5 en 6 vergelijken Cauchy hydrostatische spanningverdelingen voor de twee modellen. In de CPOR is σ_m bijna nul voor het grootste deel van de lithiatieperiode, wat verklaart waarom de modellen daar overeenkomen. Buiten CPOR divergeert σ_m, wat verschillen in chemisch potentieel en dus flux versterkt.

Conclusies

We hebben een chemisch-potentieel expressie geïntroduceerd die geldig is voor snelle deformatie, de lage-snelheid assumptie inherent in traditionele modellen eliminerend. Vergelijking met de traditionele formulering onthult verwaarloosbare verschillen onder galvanostatische cycli maar significante afwijkingen onder potentiostatische controle, gedreven door deformatiesnelheid. Een kritieke straal—de CPOR—ontstaat waar fluxes en spanningen vrijwel identiek zijn, bewegend dichter bij het oppervlak bij hogere stromen. Deze bevindingen benadrukken de noodzaak van het ontwikkelde model voor accurate voorspelling van spannings evolutie in silicium nanodeeltje elektroden tijdens snelle lading.