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

3. 반응 동역학: Butler-Volmer

학습 목표
  • Butler-Volmer 식을 전이상태 이론에서 정직하게 유도하고, 전위 이동이 두 활성화 장벽에 나뉘며 전달계수의 합이 1이 됨을 (선언이 아니라) 증명한다.
  • 식에서 두 극한을 읽어낸다: 작은 과전압의 선형(저항 같은) 영역과 전하전달저항 RctR_\mathrm{ct}, 그리고 큰 과전압의 Tafel(로그 전류) 영역과 고정된 기울기.
  • 하나의 구체적 계면에 대해 과전압에서 전류밀도를, 같은 파라미터에서 RctR_\mathrm{ct}나 Tafel 기울기를 계산한다.
요약
  • 동역학은 속도, 곧 주어진 전류로 반응을 밀 때 드는 과전압 η\eta의 가격을 정합니다.
  • Butler-Volmer 식 j=j0[exp(αaFη/RT)exp(αcFη/RT)]j=j_0[\exp(\alpha_a F\eta/RT)-\exp(-\alpha_c F\eta/RT)]η\eta에 따라 반대로 움직이는 두 활성화 장벽에서 유도됩니다.
  • 두 극한이 전부입니다: 작은 η\eta선형 영역(저항처럼 행동), 큰 η\etaTafel 영역(전류 10배마다 고정 전압).

열역학은 반응이 일어날 수 있는지를 말하고, 동역학은 주어진 속도로 일어나게 하려면 얼마나 세게 밀어야 하는지를 말합니다. 그 밀기가 η\eta이고, 가격표가 Butler-Volmer 방정식입니다.

양방향 교통 그림

평형 상태의 전하 전달 계면은 놀고 있지 않습니다: 정방향·역방향 반응이 같은 속도로, 각각 교환 전류밀도 j0j_0만큼 돌고 있습니다. 과전압을 걸면 에너지 지형이 기울어집니다: 한 방향은 쉬워지고 반대는 어려워지는데, 지수적으로 그렇습니다:

j=j0[exp ⁣(αaFηRT)exp ⁣(αcFηRT)]j = j_0\left[\exp\!\left(\frac{\alpha_a F \eta}{RT}\right) - \exp\!\left(-\frac{\alpha_c F \eta}{RT}\right)\right]

αa,αc\alpha_a, \alpha_c(전달계수, 보통 ≈ 0.5)는 기울임이 각 방향을 얼마나 돕는지를 정합니다 .

두 활성화 장벽에서 유도하는 Butler-Volmer유도 보기
이 유도가 딛는 것

반응이 자유에너지 장벽을 넘는 그림을 떠올려 보세요: 정방향(산화)으로 가려면 활성화 자유에너지 ΔGa\Delta G^\ddagger_a를, 역방향(환원)으로 가려면 장벽 ΔGc\Delta G^\ddagger_c를 넘어야 합니다. 전이상태 이론은 각 속도를 Arrhenius 인자로 만들므로, 두 부분 전류밀도는

ja=Fka0cRexp ⁣(ΔGaRT),jc=Fkc0cOexp ⁣(ΔGcRT)j_a = F k_a^0 c_R\,\exp\!\left(-\frac{\Delta G^\ddagger_a}{RT}\right), \qquad j_c = F k_c^0 c_O\,\exp\!\left(-\frac{\Delta G^\ddagger_c}{RT}\right)

평형(η=0\eta = 0)에서 둘은 크기가 같고 방향이 반대이며, 그 공통 크기가 교환 전류밀도 j0j_0입니다.

이제 전극 전위를 η\eta만큼 올리면 전달되는 전하의 전기 에너지가 이동하고, 그 이동분이 두 장벽에 나뉩니다: 분율 αa\alpha_a는 정방향 장벽을 낮추고 αc\alpha_c는 역방향 장벽을 높입니다,

ΔGaΔGa,0αaFη,ΔGcΔGc,0+αcFη\Delta G^\ddagger_a \to \Delta G^\ddagger_{a,0} - \alpha_a F\eta, \qquad \Delta G^\ddagger_c \to \Delta G^\ddagger_{c,0} + \alpha_c F\eta

두 분율의 합이 1이라는 것은 편의가 아니라 강제된 사실이며, 꼼꼼한 서술이라면 선언하지 않고 증명합니다.

보조정리왜 분배가 α_a + α_c = 1 을 따르는가보조정리 보기
이 유도가 딛는 것

전극 전위를 η\eta만큼 올립니다. 패러데이 법칙에 의해 전달되는 전하의 정전기 에너지는 전자 1몰당 정확히 FηF\eta만큼 바뀝니다: 이것이 기울임이 쓸 수 있는 예산 전부이며, 그 이상도 이하도 아닙니다. 전이상태는 반응물 우물과 생성물 우물 사이 반응 좌표의 어떤 분율 위치에 놓입니다. 정방향 쪽에서 거기로 내려가며 그 예산의 분율 αa\alpha_a를 차지하므로 정방향 장벽이 αaFη\alpha_a F\eta만큼 내려가고, 반대편에서 올라가는 역방향 장벽은 나머지 몫 αcFη\alpha_c F\eta만큼 올라갑니다. 반응물에서 생성물까지 자유에너지 차의 변화는 이 둘의 합이므로,

αaFη+αcFη=Fηαa+αc=1\alpha_a F\eta + \alpha_c F\eta = F\eta \quad\Rightarrow\quad \alpha_a + \alpha_c = 1

단일 단계·단일 전자 전달에서 성립합니다. 구동력 FηF\eta 전체는 나뉠 뿐 새로 생기지 않습니다: 그 보존이 의 내용입니다. ∎

바뀐 장벽을 대입하고 j0j_0로 묶어내면:

j=jajc=j0[exp ⁣(αaFηRT)exp ⁣(αcFηRT)]j = j_a - j_c = j_0\left[\exp\!\left(\frac{\alpha_a F\eta}{RT}\right) - \exp\!\left(-\frac{\alpha_c F\eta}{RT}\right)\right]

지수 두 개는 계면을 기울일 때 손 닿는 거리로 들고 나는 바로 그 장벽들입니다. 동역학 법칙 전체가 활성화 장벽 하나에 대한 장부 정리인 셈입니다. ∎

뿌리 사슬TST 부분 전류평형이 j0j_0를 정함장벽을 αFη\mp\alpha F\eta 이동j=j0[eαaFη/RTeαcFη/RT]j = j_0[e^{\alpha_a F\eta/RT}-e^{-\alpha_c F\eta/RT}]
Butler–Volmer 과전압 특성 (실시간). 슬라이더로 j₀·α·T를 바꾸면 비선형 곡선(j)과 선형화 근사(jLin = η/R_ct)가 즉시 갱신됩니다.
R_ct ≈ 0.0257 Ω·m² · j(50 mV) ≈ 2.268 A/m²
-0.12-0.0600.060.12η (V)-12.0-6.00.06.012.0j (A/m²)

비선형 BV (실선)와 저과전압 선형 근사 (점선). α_a ≈ α_c = 0.5, T 낮을수록 R_ct ↑, j₀ ↑일수록 곡선이 가팔라집니다. 스튜디오 anode/cathode kinetics 패널의 k_0 값과 비교해 보세요.

위 탐색기를 움직여 보세요: 작은 η\eta에서는 두 지수가 거의 상쇄되어 곡선이 선형입니다. 계면이 전하전달저항 Rct=RT/(j0(αa+αc)F)R_\mathrm{ct} = RT / (j_0(\alpha_a+\alpha_c)F)짜리 저항처럼 행동합니다. 큰 η\eta에서는 한 지수가 지배하는 Tafel 영역: 전류 10배마다 고정된 밀리볼트를 냅니다.

Tafel 플롯과전압에 대한 log 전류밀도 플롯. 원점에서 교환 전류밀도로 내려가고, 높은 과전압에서 기울기 b의 직선이 된다.과전압 ηlog |j|0log j₀기울기 b양극 가지음극 가지선형 근처(양쪽 지수 균형)
그림 · Tafel 플롯. 원점(η=0) 근처에서는 정·역 지수가 균형을 이뤄 log|j|가 교환 전류밀도 j₀로 내려간다. 과전압이 커지면 한 지수가 지배해 직선이 되고, 그 기울기가 Tafel 기울기 b = 2.303RT/(αF)이다.
Tafel 극한: 큰 η, 한 지수가 이긴다유도 보기
이 유도가 딛는 것

계면을 한 방향으로 세게 밀어 ηRT/F\eta \gg RT/F(상온에서 약 26 mV)로 가면, 역방향 지수 exp(αcFη/RT)\exp(-\alpha_c F\eta/RT)가 0으로 붕괴하고 Butler-Volmer는 정방향 가지만 남깁니다:

jj0exp ⁣(αaFηRT)j \approx j_0\exp\!\left(\frac{\alpha_a F\eta}{RT}\right)

자연로그를 취해 과전압에 대해 풀면:

η=RTαaFlnjj0\eta = \frac{RT}{\alpha_a F}\ln\frac{j}{j_0}

상용로그로 바꾸면 모든 전기화학자가 그리는 직선 η=a+blog10j\eta = a + b\log_{10} j가 되고, 기울기는

b=2.303RTαaF118 mV/decade(αa=0.5, 25C)b = \frac{2.303\,RT}{\alpha_a F} \approx 118\ \mathrm{mV/decade} \quad (\alpha_a = 0.5,\ 25\,^\circ\mathrm{C})

이제 전류 10배마다 고정 전압 bb를 냅니다. 굼뜬 반응(작은 αa\alpha_a)은 더 가파른 기울기를 걸치고, 이것이 바로 분극 곡선에서 동역학을 읽어내는 방법입니다 . ∎

뿌리 사슬ηRT/F\eta \gg RT/F역방향 가지를 버림η=RTαaFln(j/j0)\eta = \tfrac{RT}{\alpha_a F}\ln(j/j_0)b=2.303RT/αaFb = 2.303\,RT/\alpha_a F
선형 극한: 작은 η, Ohm의 법칙이 돌아온다유도 보기
이 유도가 딛는 것

이번엔 부드러운 경우, ηRT/F|\eta| \ll RT/F. 두 지수를 1차까지 전개하면 (exp(x)1+x\exp(x)\approx 1+x):

jj0[(1+αaFηRT)(1αcFηRT)]j \approx j_0\left[\left(1 + \frac{\alpha_a F\eta}{RT}\right) - \left(1 - \frac{\alpha_c F\eta}{RT}\right)\right]

상수 11은 상쇄되고 두 선형 항은 더해집니다:

jj0(αa+αc)FRTηj \approx j_0\,\frac{(\alpha_a + \alpha_c)F}{RT}\,\eta

이것이 Ohm의 법칙 η=jRct\eta = j\,R_\mathrm{ct}이고, 여기서

Rct=RTj0(αa+αc)FR_\mathrm{ct} = \frac{RT}{j_0(\alpha_a + \alpha_c)F}

1/j01/j_0에 비례하는 면적당 저항입니다: 5장의 EIS 반원이 직접 재는 바로 그 숫자이자, 어떤 분극 곡선에서든 첫 전압 강하의 미시적 뿌리입니다. ∎

뿌리 사슬\eta\ll RT/Fex1+xe^x \approx 1+xjj0(αa+αc)FRTηj \approx j_0\tfrac{(\alpha_a+\alpha_c)F}{RT}\etaRct=RT/j0FR_\mathrm{ct} = RT/j_0 F

j0j_0가 기기의 성격인 이유

j0j_0는 반응이 애초에 얼마나 기꺼이 도는지의 척도입니다. 비교:

  • 리튬이온 삽입: j0j_0는 전해질과 고체 표면 농도 모두에 의존합니다. j0kce0.5cs0.5(cs,maxcs)0.5j_0 \propto k\, c_e^{0.5} c_s^{0.5} (c_{s,\max}-c_s)^{0.5}가 P2D 솔버가 모든 입자 표면에서 평가하는 바로 그 형태입니다(8. P2D 지배방정식) . 표면이 포화되거나 비면 0으로 사라집니다. 이것이 스튜디오에서 보는 확산 한계 컷오프의 한 뿌리입니다.
  • Pt 위 수소 산화: j0j_0가 거대합니다. PEMFC의 음극은 거의 공짜입니다.
  • Pt 위 산소 환원: j0j_0가 5에서 6자릿수 작음. 이 숫자 하나가 연료전지 분극 곡선(14. 분극 곡선 읽기)이 시작하자마자 ~0.3 V를 잃는 이유이고, 백금 로딩이 원가 이야기가 되는 이유입니다.

온도

j0j_0에는 Arrhenius 인자 j0exp(Ea/RT)j_0 \propto \exp(-E_a/RT)가 붙습니다. 이는 스튜디오 스키마의 k_Ea 파라미터입니다. 차가운 동역학은 저온 용량 붕괴의 큰 조각입니다 (나머지는 수송, 4장).

실전 예제· Butler-Volmer로 과전압에서 전류밀도 구하기

대칭 계면 αa=αc=0.5\alpha_a = \alpha_c = 0.5, 교환 전류밀도 j0=5 A/m2j_0 = 5~\mathrm{A/m^2}, T=298 KT = 298~\mathrm{K}을 잡습니다. η=0.10 V\eta = 0.10~\mathrm{V}을 걸었을 때 순 전류밀도를 구합니다.

F/RTF/RT (298 K)96485/(8.314×298)=38.9 V196485 / (8.314 \times 298) = 38.9~\mathrm{V^{-1}}
αaFη/RT\alpha_a F\eta/RT0.5×38.9×0.10=1.950.5 \times 38.9 \times 0.10 = 1.95
정방향 항 exp(1.95)\exp(1.95)7.07.0
역방향 항 exp(1.95)\exp(-1.95)0.140.14
j=j0(7.00.14)j = j_0(7.0 - 0.14)5×6.87=34 A/m25 \times 6.87 = 34~\mathrm{A/m^2}

역방향 항은 이미 정방향의 약 2 %에 불과하므로, 100 mV에서 이 계면은 사실상 Tafel 영역에 있습니다: 환원 가지는 거의 꺼졌고 jj는 하나의 상승 지수를 따라갑니다.

전류에 따라 커지는 동역학 과전압 보기

NMC811/Gr 프리셋을 C/10으로 방전한 뒤 2C로 다시 돌려, 같은 중간 SOC 지점에서 단자전압 곡선을 비교해 보세요.

  • C/10에서는 전류밀도가 작아 η\eta가 선형 영역에 있고 동역학 강하는 수 mV 수준입니다(Ohm 같은 ηjRct\eta \approx j R_\mathrm{ct}).
  • 2C에서는 벌어지는 간격이 일부는 동역학(Butler-Volmer), 일부는 수송입니다. Tafel 영역에 들어서면 동역학 몫은 전류에 대해 선형이 아니라 대략 로그로 자랍니다.
스튜디오 열기 · 설정 자동 입력

연습문제

Q1

워크드 예제의 계면(j0=5 A/m2j_0 = 5~\mathrm{A/m^2}, αa+αc=1\alpha_a + \alpha_c = 1, T=298 KT = 298~\mathrm{K})에 대해 소신호 전하전달저항 RctR_\mathrm{ct}를 계산하고, 어떤 측정이 그것을 직접 읽는지 말하라.

풀이 보기

선형 극한에서 Rct=RT/(j0(αa+αc)F)=RT/(j0F)R_\mathrm{ct} = RT / (j_0(\alpha_a+\alpha_c)F) = RT/(j_0 F). 수치로,

Rct=8.314×2985×96485=24784824255.1×103 Ωm2.R_\mathrm{ct} = \frac{8.314 \times 298}{5 \times 96485} = \frac{2478}{482425} \approx 5.1 \times 10^{-3}~\Omega\,\mathrm{m^2}.

1/j01/j_0에 비례하므로 교환 전류가 클수록 저항이 작습니다. 이 면적당 RctR_\mathrm{ct}는 정확히 5장 EIS 반원의 지름이며, 선형 영역을 벗어나지 않고 j0j_0를 재는 방법입니다.

Q2

차가운 동역학. 냉각으로 j0j_0가 절반이 되는 동안 계면이 고정된 전류밀도에서 Tafel 영역에 머문다고 하자. j0j_0가 만드는 부분에 대해 b=118 mV/decadeb = 118~\mathrm{mV/decade}를 써서 이때 늘어나는 동역학 과전압을 추정하라.

풀이 보기

Tafel 영역에서 η=(RT/αaF)ln(j/j0)=(b/2.303)ln(j/j0)\eta = (RT/\alpha_a F)\ln(j/j_0) = (b/2.303)\ln(j/j_0)입니다. 고정된 jj에서 j0j_0를 절반으로 줄이면 η\eta

Δη=b2.303lnj0,warmj0,cold=0.1182.303ln2=0.0512×0.69336 mV\Delta\eta = \frac{b}{2.303}\ln\frac{j_{0,\text{warm}}}{j_{0,\text{cold}}} = \frac{0.118}{2.303}\ln 2 = 0.0512 \times 0.693 \approx 36~\mathrm{mV}

만큼 오릅니다. 즉 j0j_0가 2배 줄면 RT/αaFRT/\alpha_a F 앞인자와 수송 항이 더하는 몫 위에 약 36 mV의 동역학 페널티가 더해집니다. 이것이 저온 전압 붕괴의 한 가닥이고, 나머지는 4장의 수송 이야기입니다.

이 모델이 다루지 못하는 것

Butler-Volmer는 계면을 단일 이상 단계로 취급합니다. 실제 계면에는 막 저항(SEI, film_r 파라미터), 흡착 단계, 이중층 충전, 촉매 열화가 얹힙니다. 그것들은 핵심 방정식이 아니라 보정들의 몫입니다.

참고문헌

[1] A. J. Bard, L. R. Faulkner (2001). Electrochemical Methods: Fundamentals and Applications, 2nd ed. Wiley. 원문 링크 ↗

[2] J. Newman, K. E. Thomas-Alyea (2004). Electrochemical Systems, 3rd ed. Wiley. 원문 링크 ↗

함께 보기

2. 열역학: Nernst와 OCV · 4. 수송: 확산과 이동 · 8. P2D 지배방정식 · Bard & Faulkner 3장; Newman 8장.

REAL

Renewable Engine Analysis Laboratory

Hanyang University ERICA

경기도 안산시 상록구 한양대학로 55 공학관V 502호

Rm #502, Engineering Bldg. V, 55 Hanyangdaehak-ro, Sangnok-gu, Ansan, Gyeonggi-do, Korea

© 2026 REAL — Department of Mechanical Engineering, Hanyang University ERICA. All rights reserved.