リチウムイオン電池Siナノ粒子電極高速変形向け改良化学ポテンシャルモデル

概要

低変形速度仮定排除した改良化学ポテンシャル式提示、Siナノ粒子電極拡散誘起応力とLi分布精密モデル化可能。定電位・定電流サイクル下で新旧式比較、定電位操作で顕著乖離示す一方定電流結果ほぼ不変。Liフラックスと静水圧応力が両モデル本質的に同一となる臨界半径初同定。

はじめに

リチエーション時体積変化最大400%示すSi等高容量電極は携帯電子機器、EV、グリッド貯蔵に不可欠。しかし均質体積膨張生成拡散誘起応力が脆性破壊誘発、サイクル寿命短縮と容量低下招く。複合電極複雑形態が純粋理論記述困難、多くモデル材料特性空間変化扱うも界面効果無視。理論考慮典型電極幾何は球形、円柱、平板、1-2次元モデル対応。

Siナノ粒子電極最近研究でリチエーション誘起相変態、形態進化、応力生成相互作用浮上。これら研究固体熱力学平衡依存、原子拡散応力影響受け逆も真。LarcheとCahn枠組みで小等方変形仮定応力依存化学ポテンシャル導入、Wu式でEshelby運動量テンソル組込み。Cuiらは有限変形拡張も導出は低変形速度（拡散相対）でのみ精度保持。Si急速充電時大体積膨張でこの仮定大き誤差招く可能性。

本研究では低速度仮定なし高速変形妥当化学ポテンシャル式導出、電極形態非依存。新モデルSiナノ粒子電極定電位・定電流制御下適用、応力とLi分布影響定量化。また両モデル収束臨界半径同定、応力影響空間限界洞察提供。

実験方法

力学方程式

Li挿入体積変化誘起。ラグランジュ座標とオイラー座標両用い変形記述：$$ \mathbf{U} = \mathbf{x} - \mathbf{X} \quad (1) $$ ここでU変位場。変形勾配テンソル：$$ \mathbf{F} = \frac{\partial \mathbf{x}}{\partial \mathbf{X}} = \mathbf{I} + \mathrm{Grad}\,\mathbf{U} \quad (2) $$

球形粒子座標(R,Θ,Φ)（ラグランジュ）、(r,θ,φ)（オイラー）で勾配テンソル簡略化：$$ \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) $$

変形をリチエーション誘起形状変化と可逆弾性変形分解：$$ \mathbf{F} = \mathbf{F}^e \mathbf{F}^c \quad (4) $$ Fᶜ等方的：$$ \mathbf{F}^c = (1+\Omega C)^{1/3}\mathbf{I} \quad (5) $$ ここでΩ偏モル体積。従って $$ \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) $$

Green-Lagrange歪みテンソルE：$$ \mathbf{E} = \frac{1}{2}(\mathbf{F}^T\mathbf{F} - \mathbf{I}) \quad (7) $$ 径方向・接線方向成分：$$ E_R^e = \tfrac{1}{2}\left[\tfrac{(1+\partial u/\partial R)^2}{(1+\Omega(R)C)^{2/3}} - 1\right] \quad (9) $$ $$ E_\Theta^e = E_\Phi^e = \tfrac{1}{2}\left[\tfrac{(1+u/R)^2}{(1+\Omega(R)C)^{2/3}} - 1\right] \quad (10) $$

線形弾性仮定で歪みエネルギー密度：$$ W = \tfrac{1}{2}\mathbf{E}^e:\mathbf{C}:\mathbf{E}^e = \det(\mathbf{F}^c)\frac{E_h}{2(1+\nu)}\left[\tfrac{\nu}{1-2\nu}(tr\,\mathbf{E}^e)^2 + tr\,\mathbf{E}^e\mathbf{E}^e\right] \quad (12) $$ ここでEₕ、νヤング率、ポアソン比。第1ピオラ-キルヒホフ応力：$$ \mathbf{P} = \det(\mathbf{F}^c)\frac{E_h}{1+\nu}\left[\tfrac{\nu}{1-2\nu}tr\,\mathbf{E}^e + \mathbf{E}^e\right]\mathbf{F}^e(\mathbf{F}^c)^{-1} \quad (13) $$

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

体力なし平衡条件：$$ \frac{\partial P_R}{\partial R} + \frac{2(P_R - P_\Phi)}{R} = 0 \quad (16) $$ 境界条件：$$ u(0,t) = 0,\; P_R(R_0,t) = 0 \quad (17) $$

物質輸送方程式

Li濃度C(X,t)とフラックスJ(X,t)満たす：$$ \frac{\partial C}{\partial t} + \mathrm{Div}\,\mathbf{J} = 0 \quad (18) $$ 球対称性で：$$ \frac{\partial C(R,t)}{\partial t} + \frac{1}{R^2}\frac{\partial(R^2 J(R,t))}{\partial R} = 0 \quad (19) $$ フラックスは化学ポテンシャル駆動力フィック則従う：$$ J = -\frac{CD}{R_g T F_{11}F_{11}}\frac{\partial \mu}{\partial R} \quad (20) $$

化学ポテンシャル全内部エネルギー密度由来：$$ \Pi = \varphi(C) + W(C,\mathbf{E}^e) \quad (22) $$ より：$$ \mu = \frac{\partial \Pi}{\partial C} = \mu_0(C) + \tau(\mathbf{E}^e,C) \quad (23) $$ ここでτ応力依存部。希薄溶液で：$$ \mu_0(C) = \mu_0 + R_g T \ln(\gamma C) \quad (24) $$ γ=1。

従来ヘルムホルツ自由エネルギーアプローチで：$$ \tau_H = \frac{\partial W}{\partial C} = -\det(\mathbf{F}^e)\sigma_m\Omega \quad (25) $$ より：$$ \mu = \mu_0 + R_g T \ln C - \det(\mathbf{F}^e)\Omega\sigma_m \quad (26) $$ ここでσₘコーシー静水圧応力：$$ \sigma_m = \tfrac{1}{3}tr(\det^{-1}\mathbf{F}\mathbf{P}\mathbf{F}^T) \quad (27) $$ 式(26)が従来化学ポテンシャルモデル。

一方ギブス自由エネルギー扱いで：$$ \tau_G = \frac{\partial W}{\partial C} \quad (28) $$ より：$$ \mu = \mu_0 + R_g T \ln C - \frac{\partial W}{\partial C} \quad (29) $$ 開発モデル。質量輸送系は(19)、(20)、(26)、(29)連結。定電位（C(R₀,t)=C_max）と定電流（J(R₀,t)=j₀(1+u/R)²）境界条件：$$ C(R_0,t) = C_{\max},\; J(R_0,t) = j_0(1+u/R)^2 \quad (30-31) $$ 初期条件：$$ C(R,0) = 0 \quad (32) $$

数値実装

解析解困難、COMSOL Multiphysicsで連結PDEs解く、Table 1材料パラメータ採用。無次元変数使用で明確化。充電状態（SOC）計算：$$ SOC = \frac{\int_0^{R_0}C(R,t)R^2 dR}{\int_0^{R_0}C_{\max}R^2 dR} \quad (33) $$ 2化学ポテンシャルモデル由来フラックス：$$ J_H = \frac{\partial \tau_H}{\partial R},\; J_G = \frac{\partial \tau_G}{\partial R} \quad (34) $$

結果と考察

Figure 1で球形Siナノ粒子定電流サイクル下Li濃度、径方向応力、箍応力プロファイル示す。開発モデル（実線）と従来モデル（三角）ほぼ同一濃度分布全SOCs。応力プロファイル微差のみ、46.7%と65.5% SOCで中心部やや上昇。

Figure 2で定電位サイクル類似データ。ここ2モデル乖離顕著、特に中間SOCs中心応力。定電位操作は表面濃度一定課すため充電速度、結果変形速度が定電流より高い。与えられたSOC全体変形類似、鍵因子は絶対変形でなく変形速度示唆。

Figure 3（Ref.[17]）で620 nm Siナノ粒子が2 V vs Liで1分以内完全膨張可能示す、高速変形モデリング適合性強調。かかる領域で2化学ポテンシャルモデル著しく異なる応力状態予測、実験応力データ未入手で定量的検証不能。

Figure 4で定電流サイクル様々充電電流下フラックス比Jₕ/J₉検証。比は粒子全域で1近傍、中心付近わずか増大、表面付近減少。電流増大で広がり拡大、変形速度感度強化示唆。全曲線臨界半径付近交差、化学ポテンシャル非依存領域（CIR）命名、モデル選択拘わらずフラックス収束。CIRは電流増大で表面側シフト。

Figures 5、6で2モデルコーシー静水圧応力分布比較。CIR内σₘリチエーション大半期間ほぼゼロ、モデル一致理由説明。CIR外σₘ発散、化学ポテンシャル差拡大、フラックス差増幅。

結論

高速変形妥当化学ポテンシャル式導入、従来モデル固有低速度仮定排除。従来式比較で定電流サイクルでは差微小、定電位制御で顕著乖離、変形速度駆動要因。臨界半径（CIR）出現、フラックスと応力実質同一、高電流で表面側移動。Siナノ粒子電極急速充電時応力進化精密予測に開発モデル必須性強調。