Modelo Potencial Químico Mejorado para Deformación Rápida en Electrodos Nanopartículas Silicio de Baterías Li-Ion

Resumen

Presentamos una expresión potencial químico refinada que elimina la suposición baja tasa deformación, permitiendo modelado preciso estrés inducido difusión y distribución litio en electrodos nanopartículas silicio. Comparando nueva expresión con tradicional bajo ciclado potenciostático y galvanoestático, demostramos divergencia pronunciada durante operación potenciostática, mientras resultados galvanoestáticos permanecen largamente inalterados. Un radio crítico —donde flujo Li y tensión hidrostática son esencialmente idénticos para ambos modelos— se identifica por primera vez.

Introducción

Electrodos alta capacidad como silicio, que sufren cambios volumétricos hasta 400 % durante litio, son fundamentales para electrónica portátil, vehículos eléctricos y almacenamiento escala red [1–3]. Sin embargo, estrés inducido difusión generado por expansión volumétrica homogénea puede desencadenar fractura frágil, acortando vida ciclo y reduciendo capacidad [4]. Electrodos compuestos exhiben morfologías complejas que desafían descripciones puramente teóricas; la mayoría modelos tratan propiedades materiales como espacialmente variables mientras ignoran efectos interfaciales. Geometrías electrodo típicas consideradas en teoría son esféricas, cilíndricas y planas, con modelos uno o dos dimensiones accordingly.

Investigaciones recientes electrodos nanopartículas silicio han destacado interacción entre transformación fase inducida litio, evolución morfológica y generación estrés [5–7]. Estos estudios basan en equilibrio termodinámico sólidos, donde difusión atómica influida por estrés y viceversa. Marco Larche y Cahn [8] introdujo potencial químico dependiente estrés asumiendo deformaciones pequeñas isotrópicas, mientras formulación Wu [9] incorporó tensor momento Eshelby. Cui et al. [10] extendió esto a deformaciones finitas, yet su derivación permanece precisa solo para bajas tasas deformación relativas a difusión. Dada gran expansión volumétrica silicio durante carga rápida, esta suposición puede llevar errores significativos.

En este trabajo, derivamos expresión potencial químico válida deformación rápida sin suposición baja tasa, independiente morfología electrodo. Aplicando nuevo modelo electrodos nanopartículas silicio bajo control potenciostático y galvanoestático, cuantificamos impacto en estrés y distribución litio. También identificamos radio crítico donde dos modelos convergen, proporcionando insight límites espaciales influencia estrés.

Métodos

Ecuaciones Mecánicas

Inserción litio induce cambio volumétrico. Describimos deformación usando coordenadas Lagrangianas y Eulerianas:
$$\mathbf{U}=\mathbf{x}-\mathbf{X}\quad (1)$$ donde U campo desplazamiento. Tensor gradiente deformación:
$$\mathbf{F}=\frac{\partial \mathbf{x}}{\partial \mathbf{X}}=\mathbf{I}+\mathrm{Grad}\,\mathbf{U}\quad (2)$$

Para partícula esférica coordenadas (R,Θ,Φ) (Lagrangiana) y (r,θ,φ) (Euleriana), tensor gradiente simplifica a
$$\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)$$

Deformación descompone en cambio forma inducido litio y deformación elástica reversible:
$$\mathbf{F}=\mathbf{F}^e\mathbf{F}^c\quad (4)$$ con F^c isotrópico:
$$\mathbf{F}^c=(1+\Omega C)^{1/3}\mathbf{I}\quad (5)$$ donde Ω volumen molar parcial. Consecuentemente,
$$\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)$$

Tensor deformación Green-Lagrange E es
$$\mathbf{E}=\frac{1}{2}\left(\mathbf{F}^T\mathbf{F}-\mathbf{I}\right)\quad (7)$$ produciendo componentes radial y tangencial
$$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)$$

Asumiendo elasticidad lineal, densidad energía deformación es
$$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)$$ donde E_h y ν módulo Young y ratio Poisson. Primer esfuerzo Piola-Kirchhoff sigue
$$\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)$$

Componentwise,
$$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)$$

Equilibrio sin fuerzas cuerpo requiere
$$\frac{\partial P_R}{\partial R}+\frac{2(P_R-P_{\Phi})}{R}=0\quad (16)$$ con condiciones frontera
$$u(0,t)=0,\;P_R(R_0,t)=0\quad (17)$$

Ecuación Transporte Masa

Concentración litio C(X,t) y flujo J(X,t) satisfacen
$$\frac{\partial C}{\partial t}+\mathrm{Div}\,\mathbf{J}=0\quad (18)$$. Bajo simetría esférica, reduce a
$$\frac{\partial C(R,t)}{\partial t}+\frac{1}{R^2}\frac{\partial\big(R^2J(R,t)\big)}{\partial R}=0\quad (19)$$. Flujo sigue ley Fick con fuerza impulsora potencial químico:
$$J=-\frac{CD}{R_gTF_{11}F_{11}}\frac{\partial \mu}{\partial R}\quad (20)$$.

Potencial químico deriva densidad energía interna total
$$\Pi=\varphi(C)+W(C,\mathbf{E}^e)\quad (22)$$, dando
$$\mu=\frac{\partial\Pi}{\partial C}=\mu_0(C)+\tau(\mathbf{E}^e,C)\quad (23)$$ donde τ parte dependiente estrés. Para solución diluida,
$$\mu_0(C)=\mu_0+R_gT\ln(\gamma C)\quad (24)$$ con γ=1.

Tradicionalmente, enfoque energía libre Helmholtz produce
$$\tau_H=\frac{\partial W}{\partial C}= -\det(\mathbf{F}^e)\sigma_m\Omega\quad (25)$$ llevando a
$$\mu={\mu}_0+R_gT\ln C-\det(\mathbf{F}^e)\Omega\sigma_m\quad (26)$$. Aquí, σ_m esfuerzo hidrostático Cauchy
$$\sigma_m=\tfrac{1}{3}tr\,(\det^{-1}\mathbf{F}\mathbf{PF}^T)\quad (27)$$. Esta expresión (26) es modelo potencial químico tradicional.

Alternativamente, tratando energía libre Gibbs produce
$$\tau_G=\frac{\partial W}{\partial C}\quad (28)$$ y consecuentemente
$$\mu={\mu}_0+R_gT\ln C-\frac{\partial W}{\partial C}\quad (29)$$, nuestro modelo desarrollado. Sistema transporte masa acopla (19), (20), (26) y (29). Condiciones frontera operación potenciostática (C(R₀,t)=C_max) y galvanoestática (J(R₀,t)=j₀(1+u/R)²) son
$$C(R_0,t)=C_{\max},\;J(R_0,t)=j_0\big(1+u/R\big)^2\quad (30-31)$$ con inicial
$$C(R,0)=0\quad (32)$$.

Implementación Numérica

Soluciones analíticas intratables; resolvimos EDPs acopladas usando COMSOL Multiphysics, adoptando parámetros materiales listados Tabla 1. Variables adimensionales usadas para claridad. Estado carga (SOC) computado como
$$SOC=\frac{\int_0^{R_0}C(R,t)R^2dR}{\int_0^{R_0}C_{\max}R^2dR}\quad (33)$$. Flujos derivados dos modelos potencial químico son
$$J_H=\frac{\partial\tau_H}{\partial R},\;J_G=\frac{\partial\tau_G}{\partial R}\quad (34)$$.

Resultados y Discusión

Figura 1 ilustra concentración litio, esfuerzo radial y esfuerzo circunferencial perfiles nanopartícula Si esférica bajo ciclado galvanoestático. Modelos desarrollado (líneas sólidas) y tradicional (triángulos) producen distribuciones concentración casi idénticas across todos SOCs. Perfiles esfuerzo difieren solo ligeramente, con modesta elevación central 46,7 % y 65,5 % SOC.

Figura 2 presenta datos análogos ciclado potenciostático. Aquí, divergencia entre dos modelos más pronunciada, particularmente esfuerzos centrales SOCs intermedios. Porque operación potenciostática impone concentración superficial constante, tasa carga —y consecuentemente tasa deformación— es mayor que modo galvanoestático. Deformación general para dado SOC permanece similar, subrayando factor clave es tasa deformación no deformación absoluta.

Figura 3 (de Ref. [17]) demuestra nanopartícula Si 620-nm puede expandir completamente dentro un minuto a 2 V versus Li, destacando relevancia modelado deformación rápida. En tales regímenes, dos modelos potencial químico predicen estados estrés marcadamente diferentes, aunque datos experimentales estrés no disponibles para validación cuantitativa.

Figura 4 examina ratio flujo J_H/J_G bajo ciclado galvanoestático varias corrientes carga. Ratio permanece cercano unidad across partícula, con ligero aumento cerca centro y disminución cerca superficie. Dispersión se ensancha con mayor corriente, indicando mayor sensibilidad tasa deformación. Todas curvas intersectan cerca radio crítico —termed región independiente potencial químico (CIR)— donde flujos convergen independientemente modelo elegido. CIR se desplaza más cerca superficie al aumentar corriente.

Figuras 5 y 6 comparan distribuciones esfuerzo hidrostático Cauchy para dos modelos. En CIR, σ_m casi cero la mayor parte período litio, explicando por qué modelos acuerdan allí. Fuera CIR, σ_m diverge, amplificando diferencias potencial químico y así flujo.

Conclusiones

Hemos introducido expresión potencial químico válida deformación rápida, eliminando suposición baja tasa inherente modelos tradicionales. Comparación con formulación tradicional revela diferencias insignificantes bajo ciclado galvanoestático pero desviaciones significativas bajo control potenciostático, impulsadas por tasa deformación. Un radio crítico —el CIR— emerge donde flujos y esfuerzos virtualmente idénticos, moviéndose más cerca superficie a corrientes mayores. Estos hallazgos destacan necesidad usar modelo desarrollado para predicción precisa evolución estrés en electrodos nanopartículas silicio durante carga rápida.

Disponibilidad Datos y Materiales

Conjuntos datos analizados en este estudio disponibles del autor correspondiente bajo petición razonable.

Abreviaturas

CIR

Región donde flujos difusión predichos dos expresiones potencial químico son esencialmente idénticos.