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

5. 임피던스(EIS)의 언어

학습 목표
  • Nyquist 플롯을 시간대 이야기로 읽는다: 고주파 절편에서 옴 저항을, 반원의 지름과 꼭대기에서 전하전달저항 RctR_\mathrm{ct}과 이중층 커패시턴스 CdlC_\mathrm{dl}을 뽑아낸다.
  • RctR_\mathrm{ct}을 Butler-Volmer 법칙의 소신호 기울기로 유도하고, 병렬 RctCdlR_\mathrm{ct}\parallel C_\mathrm{dl} 회로가 지름 RctR_\mathrm{ct}인 반원을 그림을 보인다.
  • 반무한 확산이 45° Warburg 직선을 만드는 이유를 설명하고, 측정 스펙트럼을 믿을 수 있게 하는 세 조건(선형성, 정상성, 인과성)을 말한다.
요약
  • EIS는 셀을 한 번에 한 주파수씩 흔들어, 모든 공정이 한꺼번에 답하는 대신 각자의 시간대에서 답하게 합니다.
  • Nyquist 플롯은 왼쪽에서 오른쪽으로 읽는 이야기입니다: 옴 절편, 전하전달 반원(RctCdlR_\mathrm{ct}\parallel C_\mathrm{dl}), 이어 45° Warburg 확산 꼬리.
  • 둘 다 유도됩니다: 반원은 병렬 RC에서(그 저항은 선형화한 Butler-Volmer 기울기), 45° 직선은 반무한 확산의 1/jω1/\sqrt{j\omega} 스케일링에서.

전기화학 임피던스 분광법은 가장 정중한 질문을 던집니다: 한 주파수에서 전류를 아주 조금 흔들고, 전압이 어떻게 흔들려 돌아오는지 본 뒤, 주파수를 바꿔가며 반복합니다. 그 정중함의 보상으로 셀은 한 번에 한 공정씩 대답합니다 .

왜 주파수가 공정을 분리하는가

물리 공정마다 고유 시간이 있습니다: 이중층 충전(µsms), 전하 전달(ms), 고체 확산(smin). 주어진 가진 주파수에서, 흔들림보다 훨씬 빠른 공정은 즉시 따라오고(저항처럼 보임), 훨씬 느린 공정은 아예 못 따라오며(열린 회로처럼 보임), 그 시간대 근처의 공정이 EIS가 측정하는 위상 지연을 만듭니다. 주파수를 쓸면 시간 사다리를 한 칸씩 내려가는 셈입니다.

Nyquist 플롯 읽기

Nyquist 선도 해부실수 임피던스 축 위에 옴 절편, 전하전달 반원, 45도 Warburg 확산 꼬리가 표시된 Nyquist 선도.Re(Z)−Im(Z)R_sR_ct정점 f = 1/2πR_ctCWarburg 45° (확산)고주파저주파
그림 · Nyquist 선도 해부. 고주파 절편은 옴 저항 R_s, 반원의 지름은 전하전달 저항 R_ct, 반원 정점 주파수는 1/(2πR_ctC_dl), 저주파 45° 꼬리는 확산(Warburg). 각 특징이 물리적 주소 하나에 대응.

Im(Z)-\mathrm{Im}(Z) vs Re(Z)\mathrm{Re}(Z)를 주파수마다 한 점씩 그리면:

  • 실수축과 만나는 고주파 절편은 순수 저항(전해질, 멤브레인, 접촉)입니다. PEMFC에서는 멤브레인 가습 상태를 재는 표준 게이지입니다.
  • 반원(들)은 이중층 커패시턴스와 병렬인 전하전달저항 RctR_\mathrm{ct}입니다. 3장의 BV 탐색기가 보여줬듯 Rct1/j0R_\mathrm{ct} \propto 1/j_0: 큰 반원은 굼뜬 동역학이며, 연료전지 스펙트럼을 산소 전극이 지배하는 이유입니다.
  • 저주파의 45° 꼬리(Warburg)가 용량성 직선으로 꺾이는 것은 고체 확산입니다. 시작 주파수에 rs2/Dsr_s^2/D_s가 담겨 있습니다. 4장의 고체 시계를 전기적으로 읽는 셈입니다.
전하전달 반원 (병렬 RC)유도 보기
이 유도가 딛는 것

계면을 미소 교류 전류가 흐르는 두 병렬 경로로 봅니다: RctR_\mathrm{ct}짜리 패러데이 가지와 이중층을 충전하는 커패시턴스 가지 CdlC_\mathrm{dl}. 이 패러데이 가지의 저항은 자유 매개변수가 아니라, 전류-과전압 법칙의 영점 바이어스 기울기이며, 이를 먼저 증명합니다.

보조정리R_ct, 선형화한 Butler-Volmer 기울기로서보조정리 보기
이 유도가 딛는 것

Butler-Volmer 법칙은 패러데이 전류밀도를 η\eta에 대한 두 지수함수의 차로 씁니다:

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]

EIS는 η\eta를 아주 작게 유지하므로 η=0\eta=0 근처에서 선형화합니다. 미분하여 η=0\eta=0에서 평가하면(두 지수가 모두 1),

djdηη=0=j0FRT(αa+αc)\left.\frac{dj}{d\eta}\right|_{\eta=0} = j_0\frac{F}{RT}\left(\alpha_a + \alpha_c\right)

전하전달저항은 그 역기울기 Rct(dη/dj)η=0R_\mathrm{ct} \equiv (d\eta/dj)|_{\eta=0}입니다:

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

따라서 교환전류밀도 j0j_0가 작으면(굼뜬 동역학) 저항이 커집니다 . ∎

RctR_\mathrm{ct}이 정해지면 두 가지가 교류 전류를 병렬로 나르므로 어드미턴스가 더해지고(j=1j=\sqrt{-1}),

Z(ω)=(1Rct+jωCdl)1=Rct1+jωRctCdlZ(\omega) = \left(\frac{1}{R_\mathrm{ct}} + j\omega C_\mathrm{dl}\right)^{-1} = \frac{R_\mathrm{ct}}{1 + j\omega R_\mathrm{ct}C_\mathrm{dl}}

τ=RctCdl\tau = R_\mathrm{ct}C_\mathrm{dl}로 두고 켤레복소수를 곱해 실부와 허부로 나누면:

Z=Rct1+(ωτ)2,Z=Rctωτ1+(ωτ)2Z' = \frac{R_\mathrm{ct}}{1+(\omega\tau)^2}, \qquad -Z'' = \frac{R_\mathrm{ct}\,\omega\tau}{1+(\omega\tau)^2}

둘에서 ω\omega를 소거하면 짧은 대수로

(ZRct2)2+(Z)2=(Rct2)2\left(Z' - \tfrac{R_\mathrm{ct}}{2}\right)^2 + (Z'')^2 = \left(\tfrac{R_\mathrm{ct}}{2}\right)^2

즉 중심 (Rct/2,0)(R_\mathrm{ct}/2,\,0), 반지름 Rct/2R_\mathrm{ct}/2인 원, 곧 Nyquist 반원입니다. 지름은 정확히 RctR_\mathrm{ct}이고 꼭대기는 ω=1/(RctCdl)\omega = 1/(R_\mathrm{ct}C_\mathrm{dl})에 있어서, 호의 폭은 동역학을, 꼭대기 주파수는 이중층을 읽어냅니다. ∎

뿌리 사슬BV 선형화 Rct\to R_\mathrm{ct}병렬 RC 어드미턴스Z=Rct/(1+jωτ)Z=R_\mathrm{ct}/(1+j\omega\tau)지름 RctR_\mathrm{ct}인 원
Warburg 45° 직선 (반무한 확산)유도 보기
이 유도가 딛는 것

저주파에서는 전류가 표면으로 오는 반응물의 확산에 막힙니다. 농도를 c(x,t)=c0+Δc(x)ejωtc(x,t) = c_0 + \Delta c(x)\,e^{j\omega t}로 섭동해 Fick 제2법칙에 넣으면, 시간미분이 jωj\omega를 끌어내리고 xx에 대한 상미분방정식만 남습니다:

jωΔc=Dd2Δcdx2d2Δcdx2=jωDΔcj\omega\,\Delta c = D\frac{d^2 \Delta c}{dx^2} \quad\Rightarrow\quad \frac{d^2\Delta c}{dx^2} = \frac{j\omega}{D}\,\Delta c

반무한 전극 안으로 감쇠하는 해는

Δc(x)=Δc(0)exp ⁣(jωD  x)\Delta c(x) = \Delta c(0)\,\exp\!\left(-\sqrt{\tfrac{j\omega}{D}}\;x\right)

표면 교류 플럭스는 Fick 제1법칙에서 N~=DΔc(0)=jωD  Δc(0)\tilde{N} = -D\,\Delta c'(0) = \sqrt{j\omega D}\;\Delta c(0)이고, 전류는 I~N~\tilde{I}\propto\tilde{N}입니다. 표면 전위가 표면 농도에 실수 계수 E/c\partial E/\partial c로 반응하므로, 임피던스 ZW=E~/I~Z_W = \tilde{E}/\tilde{I}는 그 주파수 의존성을 물려받습니다:

ZW1jω=1ωejπ/4=12ω(1j)Z_W \propto \frac{1}{\sqrt{j\omega}} = \frac{1}{\sqrt{\omega}}\,e^{-j\pi/4} = \frac{1}{\sqrt{2\omega}}\,(1 - j)

실부와 허부가 모든 주파수에서 같으므로, Nyquist 평면에서 임피던스는 4545^\circ 직선을 그립니다. 저항성과 용량성의 정확히 중간에 기운 확산의 지문입니다. ∎

뿌리 사슬Fick 제2법칙 섭동Δcejω/Dx\Delta c \propto e^{-\sqrt{j\omega/D}\,x}ZW1/jωZ_W\propto 1/\sqrt{j\omega}45° 직선
실전 예제· 반원 하나에서 R_ct, C_dl, 꼭대기 주파수 읽기

어떤 산소 전극이 25 °C에서 교환전류밀도 j0=5 A/m2j_0 = 5~\mathrm{A/m^2}, 대칭 전달 (αa=αc=0.5\alpha_a=\alpha_c=0.5이라 αa+αc=1\alpha_a+\alpha_c=1), 이중층 커패시턴스 Cdl=0.20 F/m2C_\mathrm{dl}=0.20~\mathrm{F/m^2}을 가진다고 하자. 반원의 지름과 꼭대기 주파수를 구하라.

RT/FRT/F (25 °C)0.0257 V0.0257~\mathrm{V}
Rct=RT/(j0(αa+αc)F)R_\mathrm{ct} = RT/(j_0(\alpha_a+\alpha_c)F)0.0257/5=5.1 mΩm20.0257/5 = 5.1~\mathrm{m\Omega\,m^2}
반원 지름Rct=5.1 mΩm2R_\mathrm{ct} = 5.1~\mathrm{m\Omega\,m^2}
τ=RctCdl\tau = R_\mathrm{ct}C_\mathrm{dl}5.1×103×0.20=1.0 ms5.1\times10^{-3}\times0.20 = 1.0~\mathrm{ms}
꼭대기 주파수 f=1/(2πτ)f = 1/(2\pi\tau)155 Hz\approx 155~\mathrm{Hz}

j0j_0을 절반으로 줄이면(예: 셀을 식히면) RctR_\mathrm{ct}이 두 배가 되어 반원 폭도 두 배가 되지만, 옴 절편은 건드리지 않습니다: 손실을 그 물리적 주소에서 곧바로 읽어낸 것입니다.

이 플랫폼의 EIS

Randles 등가회로옴 저항 R_s가 직렬로, 전하전달 저항 R_ct와 Warburg가 이중층 용량 C_dl과 병렬로 연결된 Randles 회로.R_sC_dlR_ctWZ(ω)
그림 · Randles 등가회로. 전해질·접촉의 옴 저항 R_s가 직렬로 오고, 반응 계면에서는 이중층 용량 C_dl이 전하전달 저항 R_ct(+저주파 Warburg 확산 W)와 병렬로 놓인다. 이 회로가 앞의 Nyquist 반원 하나와 45° 꼬리를 정확히 만든다.

Battery 1D에는 가상 EIS 벤치가 있습니다: 솔버가 P2D 전체 모델에 미소 정현 섭동을 가해 스펙트럼을 내놓으므로, 셀을 식히거나 j0j_0를 줄이면 반원이 자라는 걸 볼 수 있습니다 . 조작법은 EIS 분석 가이드, 실행은 Battery 1D에서 (Analysis mode → EIS).

Nyquist 탐색기 (실시간). R_s·R_ct·C_dl을 끌면 반원이 즉시 갱신됩니다: R_s는 절편, R_ct는 지름, C_dl은 정점 주파수.
정점 f ≈ 398 Hz
510152025Re(Z) (mΩ·cm²)036912−Im(Z)

단일 RC 반원(Warburg 꼬리 없음). 고주파 절편 = R_s, 지름 = R_ct, 정점 주파수 = 1/(2πR_ctC_dl). 실제 스펙트럼은 저주파에 확산(Warburg) 가지가 더 붙습니다(위 그림).

작은 글씨

EIS는 세 조건에서만 정직합니다: 선형성(작은 진폭), 정상성(스윕 중 셀이 표류하면 안 됨; 천천히 방전 중인 배터리는 저주파에서 이를 위반), 인과성. Kramers-Kronig 관계식은 측정 스펙트럼을 정확히 이 가정들에 대해 감사하기 위해 존재합니다 . 그리고 Nyquist 플롯은 증명이 아닙니다: 서로 다른 회로가 같은 호를 만들 수 있습니다. 등가회로 피팅은 식별이 아니라 가설 검정입니다.

연습문제

Q1

측정된 Nyquist 호가 실수축을 8 mΩm28~\mathrm{m\Omega\,m^2}28 mΩm228~\mathrm{m\Omega\,m^2}에서 지나고, 꼭대기(Z-Z'' 최대)가 f=80 Hzf=80~\mathrm{Hz}에 있다. 옴 저항, 전하전달저항, 이중층 커패시턴스를 읽어내라.

풀이 보기

고주파 절편이 옴 저항으로 RΩ=8 mΩm2R_\Omega = 8~\mathrm{m\Omega\,m^2}입니다. 반원의 지름이 전하전달저항으로 Rct=288=20 mΩm2R_\mathrm{ct} = 28 - 8 = 20~\mathrm{m\Omega\,m^2}입니다. 꼭대기는 ωτ=1\omega\tau = 1이므로 τ=1/(2πf)=1/(2π80)1.99 ms\tau = 1/(2\pi f) = 1/(2\pi\cdot80) \approx 1.99~\mathrm{ms}이고, Cdl=τ/Rct=1.99×103/0.0200.10 F/m2C_\mathrm{dl} = \tau/R_\mathrm{ct} = 1.99\times10^{-3}/0.020 \approx 0.10~\mathrm{F/m^2}입니다.

Q2

셀을 25 °C에서 5 °C로 식히면 교환전류밀도 j0j_0이 3배 작아진다. Rct1/j0R_\mathrm{ct}\propto 1/j_0을 써서, 반원 지름은 몇 배로 변하고, Warburg 꼬리의 45° 기울기는 바뀌는가? 위 두 유도에서 근거를 대라.

풀이 보기

Rct=RT/(j0(αa+αc)F)R_\mathrm{ct}=RT/(j_0(\alpha_a+\alpha_c)F)1/j01/j_0으로 스케일하므로(작은 RT/FRT/F 온도 앞계수는 20 K 동안 7 %만 변해 j0j_0 변화에 압도됨) 지름은 대략 3배로 커집니다: 동역학이 느려지며 반원이 부풀어 오릅니다. Warburg 기울기는 정확히 4545^\circ로 유지되는데, 그 각도는 1/jω1/\sqrt{j\omega} 스케일링만으로 나오고 DD와 무관하기 때문입니다. 식히면 꼬리의 시작 주파수(rs2/Dsr_s^2/D_s를 담은)만 이동할 뿐 기울기는 그대로입니다.

이 장이 다루지 못하는 것

비선형(대진폭) EIS, 부반응·흡착 중간체가 만드는 저주파 유도성 루프(실제 PEMFC 스펙트럼에 흔함), 이완시간분포(DRT) 분석.

참고문헌

[1] M. E. Orazem, B. Tribollet (2017). Electrochemical Impedance Spectroscopy, 2nd ed. Wiley. doi:10.1002/9781119363682 ↗

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

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

함께 보기

3. 반응 동역학: Butler-Volmer · 4. 수송: 확산과 이동 · 14. 분극 곡선 읽기 · Bard & Faulkner 10장; Orazem & Tribollet, Electrochemical Impedance Spectroscopy.

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.