EIS Tutorial Ver2b
EIS Tutorial Ver2b
Keywords Transfer Functions · Corrosion · Fuel Fig. 1 Representation of a black box. Taken from Orazem
Cells · Oxide Films and Tribollet.[1]
The transfer function provides a compact description Thus, for a time constant τC = 0.04 ms, the character-
of the input-output relation for a linear time-invariant istic frequency is 4.1 kHz.
(LTI) system. Due to the fact that most signals can The time constant for a faradaic reaction is given
be decomposed into summation of sinusoids via Fourier by
series, the response of a system is characterized by the
τ t = C 0 Rt (6)
frequency response of the system. A generalized system
is illustrated in Figure 2(a). The response to a step where Rt is the charge-transfer resistance, given, for
input signal X(t) shows different long and short time linear kinetics on a disk electrode, by
behaviors that can be represented as the dependence of
the transfer function on frequency. RT
Rt = (7)
The short-time behavior corresponds to high fre- nFi0
quencies, and the long-time behavior corresponds to Equations (6) and (7) yield a time constant of 0.51 ms
low frequencies. For an electrochemical system, charg- and a characteristic frequency of 310 Hz.
ing of the electrode-electrolyte interface occurs rapidly The time constant for diffusion to a rotating disk
and is associated with the high-frequency or short-time electrode is given by
response. Diffusion is a slower process with a large time
2
constant and correspondingly a smaller characteristic δN
τD = (8)
frequency. Di
The appropriate frequency range for the measure- where Di is the diffusion coefficient for the reacting
ment is determined in accordance to the system un- species and δN is the diffusion layer thickness given as
der study. For example, the phenomena associated with a function of rotation speed by
a typical electrochemical system include double-layer
charging, faradaic reactions, and diffusion. 1/3
3 1
r
ν 1
r
ν
Consider, for example, a disk electrode with radius δN,i = Γ (4/3) = 1.61
a Sci 1/3 Ω Sci 1/3 Ω
r0 = 0.25 cm and a capacity C0 = 20 µF/cm2 . The disk
(9)
is immersed in an electrolyte of resistivity ρ = 10 Ωcm.
Assume that the faradaic reaction has an exchange cur- where Sci = ν/Di is the Schmidt number which has
rent density i0 = 1 mA/cm2 and that the disk is rotat- a value Sci = 1000 for ν = 10−2 cm2 /s and a typical
ing at Ω = 400 rpm. The kinematic viscosity of the diffusivity Di = 10−5 cm2 /s. The time constant for dif-
electrolyte is ν = 10−2 cm2 /s. fusion of a species with a diffusivity of 10−5 cm2 /s is
The time constant for charging the electrode surface equal to 0.41 s. The corresponding characteristic fre-
is given by quency is 0.4 Hz.
Thus, the frequency range for an impedance mea-
τC = C0 Re (1) surement for disk electrode can be from 10 kHz to
A Tutorial on Electrochemical Impedance Spectroscopy 3
Y (t) = Y + |∆Y | cos(ωt + ϕ) (11) For an electrical or an electrochemical system, the input
is usually a potential, the output is a current, and the
where ϕ is the phase lag between the input and out- transfer function is called admittance. In the particular
put signals. An alternative representation of the time- case where the input is a current and the output is a
domain expressions is potential, the transfer function is an impedance. The
transfer function is, however, a property of the system
X(t) = X + Re{X
e exp(jωt)} (12) that is independent of the input signal. As the admit-
tance is the inverse of the impedance,
and
Ve (ω)
Y (t) = Y + Re{Ye exp(jωt)} (13) Z(ω) = (17)
ei(ω)
respectively, where X e and Ye are complex quantities Generally only the impedance is considered even if the
called phasors that are functions of frequency but are measurement corresponds to an admittance. The mea-
independent of time. The transfer function is a func- sured impedance can have a strong dependence on the
tion of frequency and is independent of both time and applied frequency. By analyzing the impedance as a
the magnitude of the input signal. While the measure- function of frequency, a transfer-function model could
ments are made in time domain, the determination of be defined which takes into account all time constants
the transfer function is obtained from subsequent anal- of the corresponding system.
ysis. Impedance spectroscopy is much more than a simple
The calculation of the transfer function at a given application of Ohm’s law, i.e., V = IR. The expres-
frequency ω is presented schematically in Figure 3. The sion V = IR represents a steady-state measurement.
ratio of the amplitudes of the output and input signals For a system consisting of a resistor, shown in Figure
4 Mark E. Orazem*, Bernard Tribollet
Ve Rt
Z= = Re + (24)
I
e 1 + jωR t Cdl
(a) from which the parameters Re , Rt , and Cdl may be
obtained easily by using graphical methods.[1, 3]
(a)
on steady-state potential as
∂iF ∂ γ̇
A= (28)
∂γ ∂V
(b)
where
Fig. 5 The polarization curve for a stationary iron electrode
in 0.5 M H2 SO4 : a) a scan including both the active and ∂γ i(25) − i(26)
passivated regions; b) zoomed portion of part (a) showing γ̇ = = (29)
∂t ΓF
the potential at which the impedance measurement was per-
formed. Data taken from Wu et al.[4].
When A = 0, the time-rate-of-change of fractional cov-
erage is independent of steady-state potential, and only
a single capacitive loop is seen in the impedance re-
The iron first oxidizes and forms a monovalent inter-
sponse. When A > 0, a low-frequency inductive loop is
mediate adsorbed on the electrode surface. This re-
visible, as is shown in Figure 6. The global impedance
action is followed by the oxidation of the intermedi-
response given in Figure 6 shows a low-frequency induc-
ate. The ferrous ion is soluble and diffuses away from
tive loop which may be associated with the formation
the electrode. Epelboin and Keddam[6], Armstrong et
of Fe(I) species adsorbed on the electrode surface.
al.[7], and Epelboin et al.[8] showed that the induc-
tive loops observed in the low-frequency impedance re- While the mathematical development is found else-
sponse could be attributed to partial coverage of the where, the results presented here show, from a qualita-
iron surface by an adsorbed intermediate. Such a film tive perspective, that the low-frequency inductive loop
is too thin to be observed directly. may be attributed to the coupling of electrochemical
reactions through an intermediate. In the case of the
The faradaic impedance for this system can be ex-
iron system, the reactions (25) and (26) are coupled
pressed as
through an adsorbed intermediate. A similar coupling
of elementary reactions has been proposed for corrosion
1 of magnesium.[9]
ZF = (27)
1 A Other forms of coupling may also be envisioned. Roy
+
Rt,eff B + jω et al.[10] suggested that low-frequency inductive loops,
often observed in proton-exchange-membrane (PEM)
where Rt,eff represents the effective charge-transfer re- fuel cells, may be attributed to oxidation of the plat-
sistance for reactions (25) and (26). The inductive term inum catalyst. The platinum dissolution was assumed
is scaled by a constant A that can be expressed in terms to occur by a reaction scheme similar to that reported
of the dependence of the steady-state faradaic current by Darling et al.,[11] by an electrochemical reaction
on the fractional coverage of the surface and the depen-
dence of the time-rate-of-change of fractional coverage Pt + H2 O PtO + 2H+ + 2e− (30)
6 Mark E. Orazem*, Bernard Tribollet
where KPt is the rate constant on a platinum site and Fig. 7 Electrochemical impedance response for a 5 cm2 nom-
KPtO is the rate constant on a platinum oxide site. The inal surface area PEM fuel cell at a current of 0.2 A/cm2 .
The line represents a model calculation based on the cou-
oxygen-reduction reaction was assumed to take place
pling of electrochemical reactions through the oxidation of a
according to platinum catalyst. The diamond is the polarization resis-
tance obtained from the slope of the polarization curve. Data
O2 + 4H+ + 4e− → 2H2 O (33) taken from Roy et al.[10].
In the special case that both the anodic and the ca-
thodic reactions are controlled by kinetics, the rate of
corrosion may be determined from impedance measure-
ments made under the open-circuit condition, or the
corrosion potential. The corrosion current may be ex-
tracted from the polarization or effective charge-transfer
resistance using the Stern–Geary relation[16]
Fig. 9 Pt 4f high-resolution spectra obtained from un-used
and used (after 600 hours of operation) MEA cathode catalyst −βc βa
layers. Taken from Roy et al.[12] icorr = (35)
2.303 × 103 Rt,eff (βa − βc )
where βa and βc are the Tafel slopes for anodic and
cathodic reactions, respectively, and
Nafionr membrane due to membrane degradation or to
carbon support corrosion. Rt,a Rt,c
Rt,eff = (36)
The high-resolution spectra of the Pt 4f region are Rt,a + Rt,c
presented in Figure 9. The binding energies of the Pt is the effective resistance associated with the parallel
4f peaks (∼ 71 eV) were consistent with those reported contributions of the charge-transfer resistances for the
for Pt metal [13]. Since the spectra have been normal- anodic and cathodic reactions. Equation (35) applies
ized, the lower Pt content on the un-used catalyst layer for systems such as the corrosion of iron in anaerobic
is evident in the lower signal-to-noise ratio of the spec- electrolytes, but does not apply for more complicated
trum obtained from this sample. In addition to a higher systems, for example, in aerated electrolytes, where the
Pt content on the used cathode catalyst layer, the Pt reactions are not strictly kinetically controlled. In more
4f peaks are broader compared to the un-used cathode complex systems, the steady-state corrosion current may
catalyst layer. There is a significant shoulder at higher be obtained from the charge-transfer resistance of the
binding energies compared to Pt metal. This is indica- corrosion reaction, i.e.,
tive of Pt oxidation. Roy et al.[12] observed that the
Pt 4f spectrum obtained from 3 monolayers of PtOx on βa
ia = (37)
Pt(100), prepared using atomic oxygen, was very sim- 2.303 × 103 Rt,a
ilar to the spectrum obtained from the used catalyst
A sophisticated model may be needed to differentiate
layer.[14] They concluded that, after extended opera-
the corrosion process from the other phenomena and
tion in the fuel cell, some of the platinum in the catalyst
reactions in the system.
layer was oxidized.
An example is provided by Frateur et al.[17], who
The attribution of the low-frequency inductive loop explored the corrosion of cast iron in EvianTM drinking
to PtOx formation is sensitive to the reversibility of water. A film is formed on the cast iron electrodes that
PtOx formation as a function of potential. While the includes a porous electrode and porous ferrous films.
reversibility of PtOx formation was not addressed by As the electrolyte is aerated, the cathodic reaction is
Roy et al.[12], Mathias et al.[15] reported that the re- reduction of oxygen, which is influenced by mass trans-
action is reversible. fer to the pore walls. The anodic reaction is dissolution
8 Mark E. Orazem*, Bernard Tribollet