23. 실리콘·복합 음극

학습 목표

실리콘이 흑연과 달리 역학을 필요로 하는 이유를 설명하고, 응력 받는 고체의 화학 퍼텐셜에서 응력-OCP 이동 $\Delta U_\mathrm{mech} = \sigma_h \Omega / F$ 를 쓴다.
흑연과 실리콘을 한 전극 안의 두 병렬 입자로 결합하고, 복합 셀의 이동을 두 입자별 이동의 전류 가중 혼합으로 유도한다.
실리콘 고체 확산계수 $D_s$ 를 민감도 파라미터로 읽고, 그 값이 1C에서 입자 내부가 얼어붙는지를 정하는 이유를 말한다.

요약

실리콘은 흑연보다 훨씬 많은 리튬을 저장하지만 리튬화 시 부피가 약 3.7배 부푼다. 이 팽창이 개방전압을 이동시키며, 흑연의 ~10 % 변화는 무시할 수 있게 해 준다.
기계적 이동 $\Delta U_\mathrm{mech} = \sigma_h \Omega / F$ 는 경로에 의존하는(소성) 실제 기여이므로, $U_\mathrm{eq}$ 에 접지 않고 그 위에 얹어 모델링한다.
복합 Gr+Si 음극은 하나의 고체 전위를 공유하는 두 입자이며, 관측되는 셀 이동은 둘의 전류 가중 혼합이고, 실리콘의 확산계수가 실리콘의 참여 여부를 정한다.

실리콘은 에너지 밀도를 위해 이 분야가 당기는 지렛대다. 순수 실리콘이나 실리콘-흑연 복합 음극은 흑연 전극 용량을 대략 두 배로 올리지만, 흑연은 결코 풀라고 하지 않던 역학 문제를 가져온다. 이 장은 이 플랫폼의 실리콘 모델을 처음부터 끝까지 다룬다.

그림 · 실리콘의 팽창. 완전 리튬화된 Li₃.₇₅Si는 부피가 원래의 약 3.7배로 부푼다(반경 ~1.5배, 점선이 원래 크기). 이 큰 변형이 표면 응력을 만들고, 그 응력이 개방전압을 이동시켜 11장의 응력-OCP 히스테리시스를 낳는다. 흑연의 ~10% 팽창과 대조된다.Gemini로 생성한 그림

실리콘이 다른 이유

흑연은 리튬을 삽입해도 부피가 약 10 % 변하지만, 실리콘은 리튬과 합금화되며 완전 리튬화 시 부피가 약 3.7배로 부푼다(일반 응력-OCP 기구는 11. 역학 (응력-전위) 모델에 있다). 이 큰 변형이 흑연 전용 모델이 무시할 수 있는 두 가지를 만든다: 리튬화된 껍질이 팽창에 저항하며 내부 응력이 쌓이고, 그 응력이 전극의 평형 전위를 이동시킨다 .

응력-OCP 이동 ΔU_mech = σ_h Ω / F유도 보기

이 유도가 딛는 것

역학은 화학 퍼텐셜을 통해 들어온다. 정수압 응력 $\sigma_h$ 아래의 고체에 대해 Larché-Cahn 결과는 리튬 화학 퍼텐셜에 기계적 항을 더한다:

$\mu_\mathrm{Li} = \mu_\mathrm{Li}^0 - \sigma_h \Omega,$

여기서 $\Omega$ 는 host 속 리튬의 부분 몰 부피(리튬 1몰 삽입당 격자가 얼마나 팽창하는지)다. 전극 전위는 리튬 화학 퍼텐셜을 전압으로 잰 값이므로 $U = -\mu_\mathrm{Li}/F$ 이다 (와 2장의 Nernst 사슬이 이 부호를 준다). $\mu_\mathrm{Li}$ 을 응력 없는 부분과 기계적 부분으로 나누면, 일반 OCV에서 벗어난 이동은

$\Delta U_\mathrm{mech} = \frac{\sigma_h \Omega}{F}.$

소성 항복 아래의 응력은 변형 이력에 의존하므로, 이 이동은 경로에 의존한다: 충전과 방전이 서로 다른 곡선을 그려 히스테리시스 루프를 연다. 그래서 이것을 평형 $U_\mathrm{eq}$ 에 접지 않고 그 위에 더한다(상태함수는 이력을 지닐 수 없다). ∎

뿌리 사슬Larché-Cahn:

\mu_\mathrm{Li} = \mu_\mathrm{Li}^0 - \sigma_h\Omega

U = -\mu_\mathrm{Li}/F

\Delta U_\mathrm{mech} = \sigma_h\Omega/F

그림 · 실리콘 응력-OCP 히스테리시스 루프. 리튬화 시 팽창이 만드는 응력이 개방전압을 올리고, 탈리튬화 시 내려 충전과 방전이 서로 다른 곡선을 그린다. 변형이 소성이라 경로에 의존하고, 두 곡선 사이 넓이가 실제로 소산되는 열이다. 이력에 의존하므로 평형 OCV에 접지 않고 따로 모델링한다.

실전 예제· 실리콘 이동과 확산 시계의 크기 어림

중간 범위의 부분 몰 부피와 완만한 탄성 응력을 잡고, 그 스윙을 견디는 나노 실리콘 입자를 크기 매긴다.

양	값
$\Omega$ (부분 몰 부피)	$1\times10^{-5}~\mathrm{m^3/mol}$
$\sigma_h$ (정수압 응력, 탄성)	$50~\mathrm{MPa}$
$\Delta U_\mathrm{mech} = \sigma_h\Omega/F$	$(5\times10^{7}\times10^{-5})/96485 \approx 5~\mathrm{mV}$
$r_s$ (나노 Si 반경)	$1\times10^{-7}~\mathrm{m}$
$D_s$ (Si 고체 확산계수)	$1\times10^{-14}~\mathrm{m^2/s}$
$\tau_\mathrm{Si}=r_s^2/D_s$	$(10^{-7})^2/10^{-14} = 1~\mathrm{s}$

50 MPa에서 입자별 이동은 수 mV이고, 실리콘이 완전 리튬화에 이르는 수백 MPa에서는 같은 식이 수십 mV를 준다. 확산 시계는 $D_s = 10^{-14}$ 에서 100 nm 입자에 대해 약 1 s뿐이라 실리콘이 방전을 따라간다. $D_s$ 를 측정된 stiff end( $10^{-19}$ )로 내리면 $\tau$ 가 $\sim 10^5$ s로 불어나 내부가 얼어붙는다.

복합: 두 입자, 하나의 전위

실리콘-흑연 복합은 두 입자 집단을 병렬로 담은 하나의 전극이다: 같은 고체 전위 $\phi_s$ 를 공유하되 각자 고유의 반응속도·확산·OCV를 가진다. 전류가 둘 사이에서 갈리고, 관측되는 기계적 이동은 실리콘만의 것이 아니다.

복합 이동은 전류 가중 혼합유도 보기

이 유도가 딛는 것

두 입자가 하나의 고체 전위에 있으므로, 전극의 총 반응 전류는 두 가지 전류의 합 $j_\mathrm{tot} = j_\mathrm{Gr} + j_\mathrm{Si}$ 이다(공유 노드에서 전하 보존). 각 입자는 자기 기계적 이동을 지닌다: 흑연은 $\Delta U_1$ (작음), 실리콘은 $\Delta U_2$ (큼). 셀 단자가 실제로 보는 이동은 한 쿨롱의 전하가 평균적으로 겪는 것, 즉 각 가지가 나르는 전류로 가중한 값이다:

$\Delta U_\mathrm{cell} = \frac{j_\mathrm{Gr}\,\Delta U_1 + j_\mathrm{Si}\,\Delta U_2}{j_\mathrm{tot}}.$

실리콘 분율이 낮으면 방전 초기에 흑연이 대부분의 전류를 나르므로, 실리콘의 큰 이동을 완만한 셀 수준 값으로 완충한다(이 플랫폼 복합 프리셋에서 C/10 중앙값 약 80 mV). N/P 균형이 이미 두 입자의 용량을 합하므로, 실리콘 슬롯은 음극 예산의 일부이지 덧붙임이 아니다. ∎

뿌리 사슬

j_\mathrm{tot} = j_\mathrm{Gr} + j_\mathrm{Si}

입자별 이동

\Delta U_1, \Delta U_2

\Delta U_\mathrm{cell} = (j_\mathrm{Gr}\Delta U_1 + j_\mathrm{Si}\Delta U_2)/j_\mathrm{tot}

입자 확산 탐색기 (실시간). D_s와 C-rate를 끌면 확산 시계 τ=r²/D와 반경 농도 프로파일이 갱신됩니다 (r_s = 5 µm 고정).

τ ≈ 42 min · Fo ≈ 1.44

D_s (m²/s)1e-14.0C-rate1.0C

확산 영역기울기 형성

확산 시계 τ = r_s²/D_s를 방전 시간과 비교: Fo = t_방전/τ 가 크면 표면과 중심이 같이 차고(잘 섞임), 작으면 표면만 차고 중심은 얼어붙는다. c(r) 형상은 정성적 스케일링이며, 실제 프로파일은 스튜디오 P2D 솔버가 적분한다.

실리콘 확산계수는 민감도 파라미터

비정질 실리콘의 $D_s$ 문헌값은 약 다섯 자릿수( $10^{-20}$ ~ $10^{-14}~\mathrm{m^2/s}$ )에 걸쳐 있고, 어느 끝도 "틀린" 것이 아니다: 하나의 상수가 아니라 민감도 파라미터다 . 이 플랫폼의 material_db는 그 범위를 싣되 loose end( $10^{-14}$ )를 기본값으로 둔다. 그 값이 실리콘을 방전 안에서 실제로 활용되게 하고, full Newman P2D가 축소 차수 SPMe와 곡선 본체에서 수 mV 이내로 일치하게 하기 때문이다. stiff 측정 끝( $10^{-19}$ )은 1C에서 내부를 얼린다( $\tau \sim 10^5$ s). 연구용으로는 유효한 경계지만 전형적 셀 거동은 아니다. 끝을 고르는 것은 수치 tractability가 아니라 셀별 물리 충실도이며, 복합 기계 솔버는 전 범위에서 수렴한다.

이 장이 다루지 않는 것

수명. 히스테리시스 루프를 여는 그 부피 스윙이 입자를 부수고 매 사이클 SEI를 다시 만들며, 이것이 실리콘의 지배적 열화 방식이자 열화 장의 주제다. 이 장은 한 번의 충전 또는 방전을 모델링할 뿐, 전극을 닳게 하는 천 번을 다루지 않는다.

참고문헌

[1] G. Bucci, S. P. V. Nadimpalli, V. A. Sethuraman, A. F. Bower, P. R. Guduru (2014). Measurement and Modeling of the Mechanical and Electrochemical Response of Amorphous Si Thin Film Electrodes J. Mech. Phys. Solids 62, 276. doi:10.1016/j.jmps.2013.10.005 ↗

[2] C.-H. Chen, F. Brosa Planella, K. O'Regan, D. Gastol, W. D. Widanage, E. Kendrick (2020). Development of Experimental Techniques for Parameterization of Multi-scale Lithium-ion Battery Models J. Electrochem. Soc. 167, 080534. doi:10.1149/1945-7111/ab9050 ↗

REAL 에너지 디바이스 가이드북