An important feature of the curve shown in Fig. 4a is that,
along with the main peak falling on the short-time domain,
one can clearly see a portion of the curve of a relatively small
slope in the time domain from 50 to 125 s. An ideal uniform
film should have only one Cottrell domain (i.e., a single value
of td), and thus the above inflection on the It1/2 vs. log t curve
implies a secondary, less expressed Cottrell domain. It can be
formally related to a higher value of td due to the film’s
non-uniformity.21 This non-uniformity may originate from
either a distribution of CP fibrils thickness (although Di
remains the same), or, alternatively, from a distribution of Di
due to the different extent of swelling of different parts of the
PTFPT film. The latter factor of the film’s non-uniformity
seems to be in line with the relatively slow kinetics of its swel-
ling in the SF solution, as verified by CV. Note that the shape
of the It1/2 vs. log t curve (i.e., the presence of a major and a
secondary Cottrell domain) was qualitatively independent of
the equilibrium potential of the film (or the related doping
level). This latter observation substantiates, as a first approxi-
mation, the use of finite-space diffusion models (elaborated
for uniform intercalation electrodes)17,18 for the determination
of td related to a specific, particular domain of the non-
uniform intercalation electrode.
on the determination of td .23 As seen from the solid line in
Fig. 4d, the transition from the Cottrell domain to the finite-
space diffusion domain occurs within a narrow time range near
t/td ¼ 0.2.23 In the short-time domain, (dlog|I(t)|/dlog t) is
equal to ꢀ0.5 (Cottrell domain), whereas in the long-time
domain, (dlog|I(t)|/dlogt) ¼ ꢀ(p2t/4td).23 Aninterestingproperty
of the plot dlog|I(t)|/dlog t vs. t/td for the case of a uniform
electrode is that it does not depend on td . In view of the con-
siderable Ohmic potential drops and the kinetic contributions
to the current transient of the PTFPT film under considera-
tion, the horizontal plateau corresponding to (dlog|I(t)|/dlog
t) ¼ ꢀ0.5 at (t/td) < 0.2, is absent in the experimental curve
in Fig. 4d. In addition, in the range of (t/td) > 0.2, the experi-
mental curve does not approach to the limiting straight line,
and clearly demonstrates a concave shape. Comparison
between Figs. 4d and c shows that the concave curve in the for-
mer figure corresponds to a gradual increase in the local value
of td . Thus, such an analysis of the shape of the dlog|I(t)|/dlog
t vs. t/td plot appears to be very useful in identifying the
non-uniform character of doping of the PTFPT film under
consideration.
It was significant to study how the non-uniform character of
doping of the PTFPT film is reflected by parallel EIS measure-
ments. Fig. 5a shows typical Nyquist plots for the PTFPT film
during its n-doping at five selected potentials in the SF solu-
tion. In the high-frequency domain, a non-closed semicircle
(HFS) appears as a result of the relatively high bulk polymer
film resistance in the SF-based solution (see insert in Fig.
5a). This semicircle was absent in the impedance spectra of
the same electrodes in the AN- and the PC-based solutions.
Thus, the above high bulk film resistance deduced from the
Nyquist plots for the PTFPT film in the SF solution could
be the major reason for deviation of the current transients
described above, from the Cottrell behavior in the short-time
domain (see Fig. 4). This high, potential-dependent bulk film
resistance probably arises because of insufficient swelling of
the film in the SF solution. Thus, a correlation between the
related PITT, CV, and EIS characteristics of the PTFPT film
in the SF solution has been obtained.
Our major focus is on the portions of the Nyquist plots
related to the medium-to-low frequency domain, because they
reflect the impact of diffusion. The insert in Fig. 5a shows that
a potential-dependent, medium-frequency semicircle (MFS)
appears on the right-hand side of the HFS. Fitting the MFS
using a simple equivalent circuit analog shows that the capaci-
tance in its maximum is of the order of 6 mF cmꢀ2, i.e., close to
a typical double-layer capacitance. Taking into account the
fact that the coupled resistance decreases as the doping level
of the film increases, the MFS was ascribed to slow, interfacial
ion-transfer kinetics. At the beginning of the n-doping (E ¼
1.6 V vs. Li/Li+), as is seen from Fig. 5a, an extensive,
depressed, HFS overlaps with the MFS. This latter semicircle,
in turn, overlaps with the very narrow (in terms of the fre-
quency range) Warburg-type response, which approaches a
limiting, distributed capacitance behavior (constant-phase
element, CPE) in the limit of the very low frequencies.
Fig. 5b presents a detailed Nyquist plot measured with the
PTFPT film during its n-doping at 1.6 V (vs. a Li ref. elec-
trode). This plot is used as an example for the analysis of td
obtained by EIS in the high- and the low-frequency domains,
and for a comparison of the td values obtained by EIS and
PITT.
The potential dependence of td related to the short time
Cottrell domains of the PTFPT film can be conveniently calcu-
lated by a simple expression:20
2
tdðEÞ ¼ ½Cdif =ðp1=2It1=2=DEÞꢂ ;
ð2Þ
where td(E) is the characteristic diffusion time constant of
a specific part of the film from which the chemical diffusion
coefficient Di (related to this specific part) can be calculated:
Di(E) ¼ l2/td(E). It is assumed that the characteristic diffusion
length l is the film’s thickness. Note that this assumption affects
the absolute value of Di rather than the plot of log Di vs. poten-
tial, since the characteristic diffusion length is most probably
invariant with the potential. The calculated value, td ¼ 12.3 s,
is marked in Fig. 4a.
Instead of determining the Cottrell parameter as minimum
in Fig. 4a, one can use the conventional I vs. t1/2 plot (see Fig.
4b). A straight line tangential to the short-time domain of the
experimental curve, crossing the origin of the coordinates,
results in exactly the same value of td as that obtained from
Fig. 4a. Involvement of slow kinetics and the secondary Cot-
trell domain are indicated in Fig. 4b, and can be compared
with those marked in Fig. 4a.
It was interesting to examine the long-time domain of the
response. This can be done using different current–time plots.
In Fig. 4c we present the same experimental data in the form
of a log|I| vs. t curve, appropriate for the identification of an
exponential decay of the diffusion current with time, in the
finite-space (long-time) domain.17,18,20,23 The diffusion time
constants td can be found from the linear portions of this plot,
whose slope is ꢀp2/4td . From Fig. 4c it is seen that in the
range of 50 < t/s < 125 the plot is linear, and the related td is
120 s. However, in the range of 10 < t/s < 13, the effective
value of td is considerably smaller, about 53 s. Fig. 4c actually
provides evidence that in the range of 10 < t/s < 50, the slope
of the log I vs. t curve gradually decreases, and thus the effec-
tive values of td increase accordingly.
In order to complete the comparison of the shape of the cur-
rent transients for uniform and non-uniform films, we used the
plot of the differential quantity, dlog|I(t)|/dlog t vs. dimension-
less time, t/td proposed by Montella.23 The solid line in Fig. 4d
relates to the case of semi-infinite (Cottrell) diffusion in the
short-time domain, and to finite-space diffusion in the long-
time domain (uniform electrode). The current was calculated
with the use of the exponent series sum,21 which was previously
used by us for modeling Li-ion diffusion in composite graphite
electrodes,24 and by Montella in his detailed analysis of the
effect of charge-transfer kinetics and Ohmic potential drops
Similar to the short-time domain of the PITT response eqn.
(2), the high-frequency Warburg-type EIS response can be
conveniently treated in terms of the differential capacitance,
Cdif ¼ Qmdy/dE (equilibrium characteristic), and the Warburg
slope, Aw (a kinetic characteristic), where Aw ¼ DRe/Doꢀ1/2
¼
DIm/Doꢀ1/2 (DRe and DIm are the differences in the real
and the imaginary components of the impedance, respec-
tively, corresponding to a finite variation in the angular
2890
Phys. Chem. Chem. Phys., 2003, 5, 2886–2893