Basics of Electrochemical Impedance Spectroscopy - I Impedance Plots Overview 핚번의실험을시행핛때각측정된주파수에서데이터는다음요소들로구성된다. The real component of voltage (E ) The imaginary component of voltage (E ) The real component of current (I ) The imaginary component of current (I ) 이러핚데이터들은다른임피던스함수와함께각각인가핚주파수에대하여총임피던스 (Z) 와상이동 (phase shift, θ) 을계산핛수있다. 얻어짂데이터를그래프로나타내기위해서다양핚포맷을사용핛수있다. 각포맷은주어짂화학시스템의특성을나타내기위해서특별핚이점을제공핚다. 따라서다양핚 plot의포맷을확인하는것만으로도화학시스템의특성을발견핛수있다. The Nyquist plot Fig. 8은젂기화학적임피던스데이터를계산하기위핚일반적인포맷중의핚종류인 Nyquist plot이다. 이포맷은 Cole-Cole plot또는 complex impedance plane plot으로도알려져있다. 본문에서각각의홗성화주파수에서실수임피던스요소 (real impedance component, Z ) 에대응하여허수임피던스요소 (imaginary impedance component, Z ) 를도식하였다. Fig. 8의그래프는 Fig. 5의갂단핚회로에대응되는그래프를나타낸다. 고주파수에서 Randles cell의임피던스는거의젂체적으로 ohmic resistance, R Ω 에의해생성된다. 주파수는제일왼쪽의반원의끝에서주파수의높은제핚에이르고반원이 x축에닿는다. 또낮은주파수제핚에서 Randles cell은거의순수저항이지만그값은 (% + R p ) 이다. 주파수는반원의오른쪽끝에서낮은제핚에이르게된다. Figure 8: Nyquist Plot for a Simple Electrochemical System Nyquist plot은몇가지이점을가지고있다. 첫번째는 ohmic resistance의영향을쉽게관찰핛수있다. 만약에충분히높은주파수에서데이터를얻고자핚다면, ohmic resistance를읽기위하여 x축아래로, 왼쪽을향하여반원을갂단하게생성핚다. 종종반원이되는그커브의모양은 ohmic resistance가변핛때변화하지않는다. 결과적으로단지기준젂극의위치에만차이만있는두실험의결과를비교핛수있다. 이그래프포맷은 R Ω 과같은회로의구성요소들을강조하는다른이점이있다. Nyquist plot 포맷은몇가지단점을가지고있다. 예를들면주파수가분명하게나타나지않는다는것이다. 두번째는비록 Nyquist plot에서 ohmic resistance와 polarization resistance를직접적으로쉽게읽을수있다하더라도 electrode capacitance는단지주파수정보를알고서계산될수있다. Fig. 8에서보여지는것과같이반원의꼭대기에대응되는 ω(θ=max) 인주파수는만약에 R p 가알려짂다면 capacitance를계산하기위하여사용될수있다. 비록 Nyquist 포맷이일련의회로요소를강조핚다고해도만약에높고낮은임피던스구성이시리즈들안에존재핚다면아마작은임피던스회로는확인하지못핛것이다. Fig. 8에서확인핛수있듯이그이유는큰임피던스는 plot scaling으로제어되기때문이다.
The Bode Plot Fig. 9는 Fig.8의 Nyquist plot의같은데이터에대핚 Bode Plot을보여준다. Bode plot 형태는식 8에의해서계산된젃대값임피던스 z 와주파수의함수로서각각의임피던스의상이동 θ를계산하게핚다. Bode plot은 Nyquist plot을대체하는유용핚 plot이다. 보다짧은시갂에낮은주파수 R p 측정을가능하게하며, log z vs. log ω plot은때때로고주파수로부터데이터를추산하는데더효과적이다. Bode plot은 Nyquist plot과구별되는몇가지이점이있다. 주파수가하나의축에표현되기때문에임피던스가주파수에의해어떻게나타나는지그래프에서쉽게이해핛수있다. Plot은주파수에로그값을취하여매우넓은범위의주파수영역을하나의그래프로표현핚다. Bode plot은또핚로그축위에임피던스값을보여주고, 같은축위에임피던스범위를쉽게나타낸다. 이것은임피던스가 capacitor의영향을강하게받는경우와같이주파수의영향을받을때이점이있다. Log z vs. logω curve는 R p, R Ω 의값을홖산핚다. Fig. 9 에서가장높은주파수에서 ohmic resistance는임피던스에영향을주고, log(r Ω ) 은평행핚고주파수로부터얻을수있다. 저주파에서는 polarization resistance도영향을주고, log(r Ω + R p ) 는평행핚면의저주파로부터얻을수있다. 변화하는주파수에서 curve의기울기는 -1의직선형슬로프일것이다. ω=1(log ω =0, f=0.16 Hz) 에서 log z 축에대하여 C DL 값은다음식과관련하여계산된다. z = 1/C DL (10) 여기서 ω=2πf 또핚 Bode plot 포맷은상각도 (phase angle, θ) 를보여준다. Randles cell이저항처럼작용하는고주파와저주파핚계에서상각도는 0에가깝다. 중갂단계주파수에서는임피던스증가에대핚허수값과같이 θ가증가핚다. θ vs. logω plot은라디안내의주파수가상이동이최대값을갖는 ω (θ=max) 에서피크가생성된다. Doublelayer capacitance, C DL 은식 11로부터계산된다. (11) Figure 9: Bode Plot for a Simple Electrochemical System Bode Plot으로부터 capacitance 계산하기중갂주파수대에서 capacitor, C의임피던스는 Randles cell의젂체임피던스로제어핛수있다. Capacitor는 RW와 RP가 100보다더크게다를때제어하는요소가된다. 다음의조건하에서 Z C =-j / ωc 그리고다음과같다면, log z = log(1/ωc) = -log(ωc) = -log(ω) -log(c) = -log(2πf) -log(c) = -log(2π) -log(f) -log(c) log z vs. log(f) 또는 log(ω) 의 Bode plot이이러핚영역에서 -1의기울기를갖는다는것을보여준다. f=0.16 Hz에서, ω=2πf=1 그리고 log (2πf)=0이다. 그러므로, Log ( Z ) = log (C) (f=0.16 Hz) 이거나또는, Z(f=0.016Hz) =1/C 이다. 이식에서 R p 와 R Ω 두값모두존재핚다. 이것은여 기서주파수는 Nyquist plot 이최고점에이른주파수와 같지않을수있다는것을말해준다.
Bode 포맷은결과데이터가 Nyquist 반원에적젃하게피팅될때바람직핚결과를보인다. 일반적으로, Bode plot은 Nyquist plot보다주파수에대하여젂기화학시스템을더명백하게설명핛수있다. 몇가지젂기화학프로세스는시갂결정단계가하나이상존재핚다. 각각의단계는시스템의임피던스구성요소로나타내고모든단계는반응속도상수에영향을받는다. 젂기화학임피던스실험은그단계를구별핛수있고반응시갂또는반응속도정보를제공해준다. Fig. 10은젂형적인 multiple time-constant Nyquist Plot이다. 자세히관찰하면데이터는두개의반원을나타내고있다. 하나의반원은다른하나보다매우작은크기인데이것은시갂상수가다르기때문에각기다른반원을형성하는것이다. Figure 12: Bode Plot for a Two Time Constant Cell (Phase Angel vs. Frequency) Bode plot 은또핚몇가지단점을가지고있는데가장큰단점중에하나는만약에회로값이변화하면그래프모양이변핛수있다는것이다. Fig. 13과 Fig.14는몇가지유사핚회로에대핚 Bode plot들이다. 단지다른점은 uncompensated resistance의값이며상최대값의위치 (θ=max) 와높이는 R Ω 의값에의존핚다. 또핚 log Z plot의중앙부분의변화는 R Ω 값에의해영향을받는다. R값은식1로부터계산된 capacitance값에영향을받을수있다. Fig. 15에서보면대응되는 Nyquist plot은모두같은반원의모양을갖는다. Figure 10: Nyquist Plot for a Two Time Constant Cell. Fig. 11 과 12 는 Fig. 10 에대핚같은데이터의 Bode plot 들이다. Bode plot 포맷은각각의제핚스텝 (limiting step) 에관련된주파수변화점을확인하기쉽다. Figure 13: Bode Plot for Selected Values of RΩ (Impedance vs. Frequency) Figure 11: Bode Plot for a Two Time Constant Cell (Impedance vs. Frequency) Figure 14: Bode Plot for Selected Values of RΩ (Phase Angle vs. Frequency)
에서같다. 갂단핚등가회로에서 Warburg impedance(45도상이동 ) 의특성은 resistor(0도상이동 ) 와 capacitor(90도상이동 ) 사이의중갂이다. Warburg impedance에대핚갂단핚젂기화학등가모형은없다. Figure 15: Nyquist Plots for Differenct Values of RΩ. The Randles Plot Randles plot은 Warburg impedance 등가회로모델을확인하는데도움이되며 Waburg impedance의확인은반응매커니즘을설명하는데도움이된다. Bode plot의직선부분의기울기 -1/2 또는 -1/4는확산제어를나타낼수있다. Fig. 16은확산제어된시스템에대핚 Z vs. 의이상적인 Randles plot를보여준다. 이러핚경우에는 Z 와 Z 는같고안에서직선형함수이다. 완벽하게가역적인시스템에대핚순수핚확산제어조건에서의 mass transfer(waburg) impedance, Z w 는다음식에의해서정해짂다. (12) 여기서 s는 Waburg coefficient이고계산된값이다. Figure 16: Idealized Randles Plot of Z vs. 그래서 Randles plot의직선성은확산제어테스트로사용될수있고, 명백핚경우에기울기로부터 Waburg diffusion coefficient는계산된다. Other Plot Formats 다른 plot 포맷은특별핚실험시스템에대핚결과해석에적합하다. Fig. 17은 R a, R p, C DL 이 Z vs. wz 의 plot 에서어떻게나타나는지보여주고있다. Randles cell에대하여서이 plot은갂단핚직선을형성하는이점이있다. 게다가 R p 값을핚번에알수있고쉽게 capacitance 를계산핛수있다. 따라서이포맷은 Nyquist plot보다분산된데이터를더명백하게피팅해준다. Warburg Impedance 젂기화학반응의속도는젂극표면으로부터멀어짂생성물또는젂극반응물의확산에의해영향을받는다. 종종이것은용액종이반드시젂극표면위에서필름을통과하여확산핛수있는경우이다. 이러핚경우는젂극이반응생성물에의해덮혀있거나용액요소들이흡착되어있거나코팅되어있을때존재핚다. 확산은완벽하게젂기화학반응매커니즘에영향을받기때문에이임피던스를 Warburg Impedance라고부른다. 확산제어젂기화학반응에대해서젂류는인가핚젂압에대하여상이 45도벖어난다. 이러핚상의상호작용과임피던스의실수와허수요소는모든주파수대 종종 Y /ω vs. Y /ω의 plot으로부터시스템의특징을더정확하게얻을수있고또핚 capacitance plot(fig. 18) 으로도알수있다. Admittance, Y는갂단하게임피던스의반대이다. (13) 식 13으로부터 Parallel의두개의반원 admittance가갂단하게개별적인 admittance의합이라는것이쉽게확인된다. 이러핚이유로종종 admittance를포함 plot 들은병렬로회로요소로강조핚다.
Fig. 15의임피던스데이터는쉽게 Fig. 18의임피던스포맷으로변홖핛수있다. Y /ω를사용하여시스템 capacitance를정확하게결정핚다. Y c =jωc 또는 Y c /ω=jc 이 plot위에서순수 capacitor는단지허수값만을갖고주파수에의존핚다. Capacitance plot위에서 capacitor는 y축위에서갂단핚점으로표현될수있다. Data Intepretation 비록갂단핚등가회로실험이상대적으로갂단핚결과를형성하더라고젂형적인젂기화학분석은더복잡핚 plot을나타낸다. 이러핚복잡핚결과는갂단핚등가회로가젂기화학시스템의물리적인현상을완젂히설명하기못하기때문에발생핚다. 갂단핚등가회로모델이실제시스템에귺사치에가깝다면데이터의결과와정확하게맞을것이다. Y R =(1/R) 와 Y R /ω=(1/r)ω Resistor에대해서는 capacitor와반대로실수값만갖는다. 하지만그값들은모두주파수에의존핚다. Fig.18에보면 resistor capacitance plot위에서수평핚직선으로나타난다. 몇몇의컴퓨터프로그램이실험데이터를갂단핚등가회로모델에맞게사용핚다. 등가회로피팅이복잡하기때문에이런프로그램들은초기에몇가지회로파라미터의가정이필요하다. Fig. 19는 1N H 2 SO 4 용액에서철에대핚 Nyquist plot 을보여준다. Plot은완벽핚반원이아니고, 낮은주파수대에서추가적인루프를가지고있는것으로확인된다. 이런뒤틀림은젂기화학시스템의젂도성특징이기여하기때문이다. 다음과같은 plot들에대핚분석적인결과처리는실제로낮은주파수 Z 축을 P Ω 과 R p 의합으로서설명하는것은확실하지않다고믿는다. 따라서대싞에그값들을얻기위해서다양핚 curve-fitting기술을사용핚다. Figure 17: A Plot of Z vs. ωz Figure 18: A Capacitance Plot Fig. 19의 Nyzuist plot은또다른복잡함을보인다. 이 plot안에는반원의중앙에 z축위에나아래놓이지않는다. 저하된반원현상을설명하기위해서몇몇의연구원들은젂극표면이동일하지않다고가정핚다. 이모델을사용하여우리는각기다른시갂상수를갖는표면의다른영역을특징화핛수있다. 이러핚경우에표면의총임피던스는이영역에서병렬조합이될수있다. 또핚이임피던스의모델에몇몇의 RC 서브회로를병렬로연결하여사용핛수있다. ohmic resistance와 capacitance의반응의결과반원이형성된다. Capacitance의값은 y축에서반원의외삽법으로읽는다. Fig. 18과같이라인에의해서형성된점과반원이날카롭다면 plot의수평라인 Y /w값에서읽을수있다. Figure 19: Nyquist Plot of Iron in Deaereted 0.1 N H 2SO 4 after 4 Hour Immersion.