역학 (응력-전위, Stress-OCP) 모델

학습 목표

응력을 받는 호스트 속 리튬 원자의 화학 퍼텐셜에서 응력-전위 이동 $\Delta U_\mathrm{mech} = \sigma_h\Omega/F$ 을 유도하고, 압축과 인장에서 그 부호를 읽어낸다.
부분 몰 부피 $\Omega$ 과 정수압 응력 $\sigma_h$ 에서 입자당 전압 이동량을 어림하고, 흑연이 이를 컴포지트 셀 단자에서 ~80 mV로 완충하는 이유를 설명한다.
가역 탄성 이동과 이력 의존 소성 히스테리시스 루프를 구별하고, 루프 폭이 $U_\mathrm{eq}$ 에 접어 넣으면 안 되는 실제 소산 발열인 이유를 설명한다.

요약

실리콘의 ~4배 리튬화 팽창은 자신의 평형 전위를 이동시킬 만큼 큰 내부 응력을 만든다. 흑연은 팽창이 작아 무시하는 결합이다.
응력은 정수압 응력을 받는 고체의 화학 퍼텐셜에서 유도되는 항 $\Delta U_\mathrm{mech} = \sigma_h\Omega/F$ 을 OCV에 더한다.
소성 항복은 응력을 이력 의존으로 만들어 충·방전 히스테리시스 루프를 열며, 그 폭은 수치 아티팩트가 아니라 실제 소산 발열이다.

왜 실리콘은 흑연과 다른 모델이 필요한가

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

흑연은 리튬을 삽입해도 부피가 ~10%만 변한다. 실리콘은 리튬과 합금화되며 훨씬 크게 팽창한다. 이 플랫폼의 Si 항목은 $V/V_0 = 1 + \text{expansion}\cdot x$ 에서 정확히 $3.11$ 의 expansion 계수를 쓴다(server/p2d-schema.js: "Silicon uses 3.11, giving the expected large alloying expansion"). 완전히 리튬화된 입자는 원래 부피의 약 4.1배가 된다. 이 정도의 부피 변형률은 흑연만 있는 모델이 무시할 수 있는 두 가지를 만든다: 리튬화된 영역이 주변의 팽창 저항과 맞서면서 내부 응력이 생기고, 그 응력이 입자의 전기화학적 평형 전위를 이동시킨다.

지배 관계식

$U_\mathrm{eff} = U_\mathrm{eq} + \Delta U_\mathrm{mech}$

$U_\mathrm{eq}$ 는 일반적인 무응력 (Electrode DB 페이지의 곡선)다. 이 플랫폼은 역학적 이동량을 다음과 같이 계산한다:

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

(스키마의 omega 파라미터 설명). 여기서 $\sigma_h$ 는 국소 정수압 응력, $\Omega$ (omega, 단위 m³/mol, 유효 범위 $[0, 10^{-4}]$ )는 호스트 내 리튬의 부분 몰 부피다. $\sigma_h$ 자체는 입자가 얼마나 팽창하려 하는지(위의 팽창식)와 주변이 그 팽창을 얼마나 강하게 저항하는지에 좌우된다. 강성 법칙은 young_type(범위 $[0, 3]$ )으로 선택한다. 스키마가 명시적으로 문서화한 것은 type 2: "실리콘 농도의존 Bucci 스타일 모델", 실리콘 재료 에세이가 확산계수 인용에 쓰는 바로 그 Bucci 2014 데이터셋과 같은 계열이다 . 나머지 선택값들은 솔버 내부에서 처리되는 더 단순한 상수/클램프 강성 법칙이다. 아래의 전극 전위 계산, 즉 전위를 리튬 화학 퍼텐셜을 $-F$ 로 나눈 값으로 두는 방식은 표준 서술을 따른다 .

응력을 받는 고체의 화학 퍼텐셜에서 응력-전위 이동 유도유도 보기

이 유도가 딛는 것

호스트에 리튬 1몰을 삽입한다. 무응력(unstressed) 고체에서 그 화학 퍼텐셜은 $\mu_\mathrm{Li}^0(\theta)$ , 즉 $U_\mathrm{eq}$ 뒤에 있는 평범한 조성 의존성이다. 응력을 받는(stressed) 고체에서는 추가 비용이 있다: 삽입된 리튬이 격자를 자신의 부분 몰 부피 $\Omega$ 만큼 팽창시키고, 그 팽창을 국소 응력장에 맞서 하면서 화학 퍼텐셜이 이동한다. 정수압 응력 $\sigma_h$ (인장이 양수)를 받는 고체에 대한 Larché-Cahn 결과는

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

이다. 압축( $\sigma_h < 0$ )에서는 $-\sigma_h\Omega > 0$ 이 화학 퍼텐셜을 올린다. 이미 짓눌린 격자에 리튬을 밀어 넣기가 더 어렵다는, 직관 그대로의 결과다. 전극 전위는 리튬 화학 퍼텐셜을 $-F$ 로 나눈 값(리튬당 전자 하나)이므로, $\mu$ 의 이동은 곧바로 전위의 이동으로 이어진다:

$\Delta U_\mathrm{mech} = -\frac{\Delta\mu_\mathrm{Li}}{F} = -\frac{-\sigma_h\Omega}{F} = \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

탄성 vs 소성 거동

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

작은 변형률에서는 입자가 탄성적으로 거동한다: 응력이 가역적으로 쌓이고 풀리며, $\Delta U_\mathrm{mech}$ 는 충전·방전에서 같은 곡선을 그릴 것이다. 실리콘의 팽창은 실용적인 Si 함량에서 순수 탄성으로 감당하기엔 너무 크다: 항복 응력(yield_stress, 단위 Pa, 범위 $[0, 2\times10^{10}]$ )을 넘으면 재료가 소성 변형한다(plastic 0/1: "OCV 히스테리시스 루프를 여는 소성 배후응력을 켠다"). 소성 변형은 비가역적이다: 역학 에너지를 회복 가능한 변형으로 저장하는 대신 열로 소산시킨다. 응력 상태가 현재 리튬 함량만이 아니라 어떻게 거기에 도달했는지에도 좌우되기 때문에, 충전과 방전은 서로 다른 단자 전압 곡선을 그린다, 곧 기계적 히스테리시스 루프다(C/10로, 충·방전 용량 축, mech ON/OFF). 흑연과 실리콘이 서로 다른 전위로 병렬 공존하므로 하나의 셀 OCV로 그릴 수 없어, 루프는 실제 단자 전압을 쓴다. 관측 가능한 기계적 효과는 mech ON 루프와 mech OFF 루프 사이의 간격이다(C/10에서 중앙값 ~80 mV, 흑연이 고갈되어 실리콘이 전류를 전담하는 저SOC에서 약 250에서 300 mV까지). 실리콘 입자 수준에서는 이동량이 몇 배 크지만 흑연이 셀 단자에서 완충한다.

크기 감(order-of-magnitude) 예시

실전 예제· 응력-전위 이동과 실리콘 확산 시계

슬라이더 중간값의 부분 몰 부피와, 아직 탄성 영역인 적당한 응력을 대입하고, 그런 입자가 팽창 스윙을 견디게 해 주는 나노 실리콘 확산 시간을 가늠한다.

양	값
$\Omega$ (부분 몰 부피)	$1\times10^{-5}~\mathrm{m^3/mol}$
$\sigma_h$ (정수압 응력, 탄성)	$50~\mathrm{MPa} = 5\times10^{7}~\mathrm{Pa}$
$F$ (패러데이 상수)	$96485~\mathrm{C/mol}$
$\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}$ (100 nm)
$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에서 낮은 GPa)에서는 같은 식이 수백 mV를 주고, 흑연이 이를 셀 단자에서 중앙값 ~80 mV 히스테리시스로 완충한다. 1초의 확산 시계는 실리콘을 ~100 nm로 분쇄하는 이유다. 마이크론 크기 입자라면 $\tau_\mathrm{Si}$ 가 $10^4$ 배 길어져, 고르게 리튬화되기 훨씬 전에 응력 구배로 파쇄될 것이다.

$\Omega \sim 10^{-5}\ \text{m}^3/\text{mol}$ ( $[0,10^{-4}]$ 슬라이더의 중간값)와 적당한 정수압 응력 $\sigma_h \sim 50\ \text{MPa}$ (항복 한계 $2\times10^{10}$ Pa보다 훨씬 낮음, 즉 아직 탄성 영역)를 대입하면: $\Delta U_\mathrm{mech} = \sigma_h\Omega/F \approx (5\times10^7 \times 10^{-5}) / 96485 \approx 5\ \text{mV}$ . 실리콘이 완전 리튬화에 가까워지며 실제로 도달하는 더 높은 응력(수백 MPa~낮은 GPa)에서는 같은 식이 입자 기준 수백 mV 이동량을 준다. 컴포지트 셀 단자에서는 이보다 훨씬 작게 나타나며(흑연이 완충), 관측 가능한 기계적 히스테리시스는 C/10에서 ~80 mV, 실리콘이 전류를 전담하는 저SOC 쪽에서 더 커진다.

토글 읽는 법

mech(0/1): 마스터 스위치, 기본값 꺼짐.
omega(Ω, m³/mol, $[0,10^{-4}]$ ): 단위 응력당 전위 이동량.
expansion: $V/V_0 = 1 + \text{expansion}\cdot x$ 의 기울기; Si = 3.11.
young_type(0 ~ 3): 영률 법칙 선택자; type 2 = 농도의존(Bucci 스타일, 실리콘).
plastic(0/1) + yield_stress(Pa, $[0, 2\times10^{10}]$ ): 입자가 언제 소성 유동하는지; 이것이 히스테리시스 루프를 실제로 여는 스위치다.

연습문제

실리콘 입자가 부분 몰 부피 $\Omega = 8\times10^{-6}~\mathrm{m^3/mol}$ 로 정수압 압축 $\sigma_h = -200~\mathrm{MPa}$ 을 받고 있다. $\Delta U_\mathrm{mech}$ 을 계산하고, 그 부호를 제시하며, 입자에 리튬을 밀어 넣는 관점에서 그 부호가 무엇을 뜻하는지 물리적으로 설명하라.

풀이 보기

$\Delta U_\mathrm{mech} = \sigma_h\Omega/F = (-2\times10^{8}\times8\times10^{-6})/96485 \approx -1600/96485 \approx -17~\mathrm{mV}$ . 이동량은 음수다. 압축은 리튬 화학 퍼텐셜을 올리고( $\mu_\mathrm{Li} = \mu_\mathrm{Li}^0 - \sigma_h\Omega$ 에서 $\sigma_h < 0$ 이면 $-\sigma_h\Omega > 0$ ), $U = -\mu_\mathrm{Li}/F$ 이므로 유효 전위는 내려간다. 물리적으로, 이미 짓눌린 격자는 리튬을 더 받기를 거부하므로, 리튬을 밀어 넣으려면 더 낮은 전위가 필요하다(충전은 더 어렵고 방전은 더 쉬움). 직관 그대로의 결과다.

스튜디오 실습. 컴포지트 Gr+Si 프리셋을 mech 켠 상태로 C/10에서 방전한 뒤 충전하라. plastic 토글을 이용해 가역 탄성 이동과 이력 의존 소성 히스테리시스 폭을 분리하는 절차를 설계하라.

풀이 보기

같은 omega에서 루프를 두 번 돌린다. 한 번은 plastic = 0, 한 번은 plastic = 1이다. plastic = 0이면 응력이 순수 탄성이라 $\Delta U_\mathrm{mech}$ 가 충전과 방전에서 같은 곡선을 되짚어 루프가 닫힌다. 남는 간격은 물리가 아니라 수치 오차다. plastic = 1이면 배후응력이 상태를 이력 의존으로 만들어 충전과 방전이 벌어지고, 둘러싸인 루프 면적이 실제 소산 발열이다. plastic을 켤 때만 나타나는 폭이 소성 기여분이고, 두 실행이 공유하는 것이 가역 탄성 이동이다. 이것이 상태함수인 $U_\mathrm{eq}$ 와 루프를 따로 모델링하는 이유이기도 하다.

이 모델이 다루지 못하는 것

이것은 집중형(lumped) 단일 입자 응력 모델이지, 파괴역학(fracture-mechanics) 모델이 아니다. 균열 개시, 입자 파쇄, 그로 인한 SEI 재성장(실제 실리콘의 수백 사이클 열화를 지배하는 현상)을 예측하지 않는다. $\Delta U_\mathrm{mech}$ 와 루프 폭을 1차(first-order) 전압·발열 보정치로 다루고, 수명 예측으로 오해하지 말 것. 이 빌드에서는 음극의 2번 슬롯(Si)에만 배선되어 있다. 양극 입자에는 응력 모델이 없다.

참고문헌

[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] J. Newman, K. E. Thomas-Alyea (2004). Electrochemical Systems, 3rd ed. Wiley. 원문 링크 ↗

함께 보기

SOC·stoichiometry · 결과 해석 · 재료 에세이: 실리콘 (Electrode DB 페이지). 이론: Larché & Cahn (1985), Acta Metall. 33, 331; Christensen & Newman (2006), J. Solid State Electrochem. 10, 293.