Electrochemical Impedance Spectra of Dye-Sensitized Solar Cells : Fundamentals and Spreadsheet Calculation

Electrochemical impedance spectroscopy (EIS) is one of the most important tools to elucidate the charge transfer and transport processes in various electrochemical systems including dye-sensitized solar cells (DSSCs). Even though there are many books and reports on EIS, it is often very difficult to explain the EIS spectra of DSSCs. Understanding EIS through calculating EIS spectra on spreadsheet can be a powerful approach as the user, without having any programming knowledge, can go through each step of calculation on a spreadsheet and get instant feedback by visualizing the calculated results or plot on the same spreadsheet. Here, a brief account of the EIS of DSSCs is given with fundamental aspects and their spreadsheet calculation. The review should help one to develop a basic understanding about EIS of DSSCs through interacting with spreadsheet.


Introduction
Impedance spectroscopy is a powerful method for characterizing the electrical properties of materials and their interfaces [1][2][3][4].When applied to an electrochemical system, it is often termed as electrochemical impedance spectroscopy (EIS); examples of such systems are electrochemical cells such as fuel cells, rechargeable batteries, corrosion, and dye-sensitized solar cells (DSSCs) [2,3].Recently, EIS has become an essential tool for characterizing DSSCs [5][6][7][8][9][10][11][12][13][14][15][16][17].Typically, a dye-sensitized solar cell (DSSC) is composed of a ruthenium dye loaded mesoporous film of nanocrystalline TiO 2 on fluorine-doped tin oxide (FTO) glass substrate as photoelectrode (PE), an iodide/triiodide (I − /I 3 − ) based redox electrolyte solution, and a Pt coated FTO glass substrate as counter electrode (CE) [14,15,[18][19][20].Analysis of EIS spectrum of a DSSC provides information about several important charge transport, transfer, and accumulation processes in the cell.These are (i) charge transport due to electron diffusion through TiO 2 and ionic diffusion in the electrolyte solution; (ii) charge transfer due to electron back reaction at the FTO/electrolyte interface and recombination at the TiO 2 /electrolyte interface and the regeneration of the redox species at CE/electrolyte interfaces; and (iii) charging of the capacitive elements in the cells including the interfaces, the conduction band, and surface states of the porous network of TiO 2 [2,7,15,[21][22][23].Even though there are many books and reports on EIS, it is often very difficult to explain the EIS spectra of DSSCs.Moreover, the details of EIS calculation always remain under several layers of programming abstraction and thus cannot be accessed by the user.
Calculating EIS spectra on spreadsheet can be a powerful approach as the user, without having any programming knowledge, can go through each step of calculation on a spreadsheet and get instant feedback by visualizing the calculated results or plot on the same spreadsheet.From our experience of learning EIS of DSSCs from scratch, we found that it was far more easy and fun to learn EIS through spreadsheet calculation than trying to decipher the abstract ideas of EIS on books or papers.
Here, a brief account of the general aspects of EIS is given with mathematical expressions and their calculation on spreadsheet (see the interactive Microsoft Excel 2010 file in the Supplementary Material available online at http://dx.doi.org/10.1155/2014/851705).Most importantly, we summarize the fundamental charge transfer processes that take place in working DSSCs and how those processes give rise to EIS spectra.
2 International Journal of Photoenergy

Electrochemical Impedance Spectroscopy
2.1.Fundamentals.Let us begin with the notion of an ideal resistor having resistance .According to Ohm's law, current () flowing through the resistor and voltage () across the two terminals of the resistor is expressed by the following relation: On the other hand, impedance is a more general concept than resistance because it involves phase difference [4].During impedance measurement, a small-amplitude modulated voltage (, ) is applied over a wide range of frequency ( = /2) and the corresponding currents (, ) are recorded, or vice versa.The resultant impedance () of the system is calculated as [1,2,4]  () =  (, )  (, ) provided that (, ) is small enough to be linear with respect to (, ), or vice versa.At a certain frequency , (, ) may have different amplitude and phase than that of (, ) depending on the nature of the charge transfer processes in the system that results in impedance of the corresponding charge transfer process.When the frequency of the applied perturbation is very low, the system is said to be driven with dc current and the impedance of the system coincides with its dc resistance ( dc ), that is, impedance with zero phase difference [2,10]: It is to be noted here that there are other response quantities related to impedance such as admittance (), modulus function (), and complex dielectric constant or dielectric permittivity () [2,4].
In general, impedance is expressed as [1,4]  () =  Re +  Im (7) or more simply as [1,4] where  Re =   = || cos  and  Im =   = || sin  are the real and the imaginary parts of the impedance, respectively.The real and imaginary parts of the impedance are related to the phase angle  as and the magnitude || as EIS data can be displayed in different ways.In the complex plane,   is plotted against   .The complex plane plots are often termed as Nyquist plots [4].In Bode plot, both log || and  are plotted against log .Sometimes, it is helpful to plot log   against log  [1,2].
In the frequency domain, current-voltage relations can be rearranged as (2).If a purely sinusoidal voltage (, ) =  0 sin() is applied across a resistor with resistance  then the current that flows through the resistor will be (, ) = (, )/ =  0 sin()/, which can be written as (, ) =  0 sin().So, the impedance of the resistor,   (), is [1] In this case, the applied voltage and the resultant current are in phase.If the voltage is applied to a capacitor having capacitance  then the resultant current is (, ) = (, )/ =  0 cos(), where  = / and  = .The above expression for the current passing through the capacitor can be written as (, ) =  0 cos( − /2) or (, ) =  0 sin(), where  0 =  0 .The impedance of the capacitor, where 1/ or in complex notation 1/ is the reactance of a capacitor and −/2 is the phase difference.According to the above description, reactance for any electrical element can be deduced using fundamental relation between current and voltage for that element as summarized in Table 1 [2,4].Analysis of EIS data is central to the study of EIS of an electrochemical system.An overview of the system of interest facilitates the translation of the charge transfer, transport, and accumulation processes in the system to an electrical circuit composed of a lump of series and parallel combination of resistors, capacitors, inductors, and so forth.The equivalent model is used to deduce the physically meaningful properties of the system.Any equivalent circuit model can be constructed using Kirchoff 's rules [1,2].For example, if two elements are in series then the current passing through them are the same and if two elements are in parallel then the voltages across them are the same.In spreadsheet, a complex number can be constructed using built-in function and the number can be operated with all the basic mathematical operators available in the spreadsheet as functions for complex numbers.Figure 1 shows such calculation implemented for impedance of a capacitor ( dl ).Thus, spreadsheet enables one to calculate EIS in its user friendly interface.Based on the above concept, all the EIS plots discussed in the present paper are calculated on spreadsheet (see the Microsoft Excel 2010 file in the Supplementary Material) unless otherwise mentioned.

Equivalent Circuit of Some Electrochemical Systems and
Their Impedance

Ideally Polarizable Electrode in Contact with Electrolyte.
An ideally polarizable electrode behaves as an ideal capacitor because there is no charge transfer across the solution/electrode interface [1].Impedance of such system can be modeled as a series combination of a resistor and a capacitor as shown in the inset of Figure 2(a).If   is the solution resistance and  dl is the double layer capacitance then the total impedance of the system becomes where    () and   dl () are the impedance for   and  dl , respectively.Equation ( 13) can be written in terms of reactance as [1,24] Rearranging ( 14), one gets Here, the real (  ) and the imaginary (  ) parts of the impedance are   and −1/ dl , respectively.Figure 2(a) shows complex plane plot of the impedance as a straight line perpendicular to the real or -axis at   , in this case   = 50 Ω, while the capacitive impedance contributes to the negative imaginary part of the impedance.At the low frequency limit ( → 0) the capacitive impedance is so large that the total impedance is infinity.Therefore, the dc resistance, (0), of the system is infinity and there is no dc current to flow through the system.As the frequency increases the capacitive impedance decreases.At the limit of very high frequency ( → ∞), the capacitor becomes short-circuited and there remains the resistance   only.However, complex plane plot does not tell us about the corresponding frequency of the impedance explicitly.In the Bode plot (Figure 2(b)), log || and  are plotted against log .The plot of impedance (red circle) versus frequency has a breakpoint, which corresponds to the characteristic frequency  = 1/   dl or characteristic time constant  = 1/ =    dl = 0.005 s of the system.On the other hand, the Bode phase plot (blue square) shows that the phase angle changes from 0 ∘ at high frequency to −90 ∘ at low frequency.

Nonpolarizable Electrode in Contact with Electrolyte.
If the electrode is nonpolarizable, then the system can be modeled by introducing a resistance  ct parallel to the capacitance  dl as shown in the inset of Figure 2(c), which is known as simplified Randle's circuit [1,24].Eventually, the circuit consists of a series connection of a solution resistance   with a parallel combination of a charge transfer resistance  ct and a double layer capacitance  dl .The impedance of the system can be written as where   () =   and  pl () is the impedance of the parallel combination of the  ct and the  dl .
Thus, ( 16) can be written in terms of reactance as [1] Here, ) are the real and imaginary parts of the impedance, respectively.Figure 2(c) shows the impedance of the system in complex plane plot.The plot has a semicircle, which is typical for a kinetic control system.When  → ∞, the capacitive impedance is short-circuited, and this eventually shunts the  ct .Therefore, only the   remains at the high frequency intercept.As the frequency decreases the capacitive impedance increases.At the low frequency intercept the capacitive impedance is infinitely large but still there is the  ct .So, the dc resistance (0) of this system is (0) =   + ct .It can be noticed from ( 17) that the maximum of the   occurs at   =   +  ct /2, which corresponds to the characteristic frequency of the charge transfer process ( max ).
In Figure 2(d), the Bode magnitude plot (red circle) of the system has two breakpoints [1].From the high frequency edge, the first breakpoint corresponds to the time constant  1 : and the second breakpoint corresponds to the time constant  2 : Here, the frequency  1 in the Bode magnitude plot (red circle, Figure 2(d)) can be calculated from ( 18) as  1 =  1 /2 = 1/2 1 = 47.75 Hz.On the other hand,  max is calculated to be 15.92Hz for  ct = 100 Ω and  dl = 100 F.The Bode phase plot (blue square, Figure 2(d)) has a maximum at around the frequency  max and 0 ∘ phase shift at both the high and low frequency limit.However, the maximum of the phase angle appears at, somewhat, higher frequency than the actual  max , which appears at the maxima of the Bode imaginary plot (Figure 2(e)) [1]. Figure 3(a) shows EIS spectra in complex plane for different values of  ct .The semicircle progressively increased as the value of  ct increased from 50 to 100 Ω while  dl remained the same.The Bode magnitude and phase plots depicted in Figures 3(b) and 3(c) clearly show increase of magnitude and decrease of characteristic frequency ( max ) with the increase of  ct .On the other hand, the complex plane plot (Figure 4(a)) remained unchanged for a fixed value of  ct and different values of  dl .Thus, (0) changes as the  ct changes while it remains fixed for all values of  dl .For both cases, the Bode magnitude and phase plots depicted in Figures 4(b) and 4(c) clearly show that  max shifts towards the low frequency edge for increasing either  ct or  dl .It is to be noted here that phase angle at the maxima decreases with the decrease of  ct while it is the same for different values of  dl .

2.2.3.
Inductance.So far we have seen that the imaginary part of the impedances for different combination of resistances and capacitors showed negative values and the spectra appeared in the first quadrant of the complex plane.However, the imaginary parts sometimes take positive values and thus the spectra appear in both first and forth quadrants due to the inductance of the contact wire, which often produces a tail at high frequencies (Figure 5(a)) [2].On the other hand, impedances of several types of solar cells show similar phenomenon, however at low frequency region, as a loop that forms an arc in the fourth quadrant (Figure 5(b)), which is attributed to specific adsorption and electrocrystallization processes at the electrode [2,4].

Constant Phase Element.
In equivalent circuit model of an electrochemical system, the capacitance  dl is often replaced by a constant phase element (CPE) to account for the deviation of the  dl from an ideal capacitor.The impedance of the CPE is expressed as [2,4,25] where   and  are the CPE prefactor and index, respectively.If the index  is equal to 1.0 the CPE coincides with a pure capacitor.Generally,  varies from 1.0 to 0.5 to fit an experimental data.The impedance corresponding to the simplified Randle's circuit with CPE (Figure 6(a)) can be expressed as Figure 6(b) shows EIS spectra for the impedance corresponding to the equivalent circuit (Figure 6(a)) in complex plane for different values of CPE index .As the value of  decreases from 1.0 to 0.5 the semicircle deviates to a depressed semicircle.In this case, the characteristic frequency  max is expressed as [2]  max = 1 From ( 22), we can see that the CPE response decelerates with the decrease of , which is evident at the second breakpoint from high frequency end of Figure 6(c).Moreover, the phase angle at the maxima decreases as well (Figure 6(d)).The equivalent capacitance ( dl ) of the electrochemical interface corresponding to the parallel combination of  ct and   of Figure 6(a) can be calculated by comparing (22) with (19) as

Semi-Infinite Diffusion.
There is another important impedance element that accounts for the impedance of redox species diffuse to and from the electrode surface.The impedance is known as semi-infinite Warburg impedance and is expressed as [1] Since 1/√ = (1 − )/ √ 2, (24) can be written as The coefficient  is defined as [1,24] where  *  and  *  are the bulk concentration of oxidant and reductant, respectively;   and   are the diffusion coefficients of the oxidant and reductant, respectively;  is the surface area of the electrode; and  is the number of electrons involved.The semi-infinite diffusion impedance cannot be modeled by simply connecting resistor and capacitor because of square root of frequency (√) [1,24].A semi-infinite transmission line (TL) composed of resistors and capacitors (Figure 7(a)) describes the impedance as a distributed element.This impedance appears as a diagonal line with a slope of 45 ∘ in complex plane plot (Figure 7(b)).In the Bode plot (Figure 7(c)), the magnitude of the impedance (red circle) increases linearly from a very low value at high frequency limit to a high value at low frequency limit and the phase angle (blue square) always remains at 45 ∘ , which is the characteristic of a diffusion process.This kind of diffusion phenomenon is seen where diffusion layer has infinite thickness.

Randle's Circuit.
If the kinetic control process as discussed in Section 2.2.2 is coupled with mass transfer process then the simplified Randle's circuit can be modified by introducing Warburg impedance () as shown in the inset of Figure 7(d) to model the mixed control process [1,24].The model of this mixed control system is known as Randle's circuit.The complex plane plot (Figure 7(d)) shows that the impedance of faradic process appears as a semicircle at high frequency edge and the diffusion process appears as a diagonal line with a slope of 45 ∘ at the low frequency edge.The Bode magnitude plot (red circles, Figure 7(e)) of the same system has three breakpoints, in the order of decreasing frequency; the first two breakpoints are similar to that of the case for kinetic control process, which is modeled as simplified Randle's circuit, and the last one corresponds to the diffusion process.The Bode phase plot (blue squares, Figure 7(e)) is similar to the Bode phase plot for simplified Randle's circuit except at the low frequency region where phase angle gradually increases and at the limit of low frequency it reaches 45 ∘ due to diffusion process.If the time constant (  = 1/ max =  ct  dl ) of the faradic or charge transfer kinetics is too fast compared to the time constant (  =  2 ct /2 2 ) of diffusion process then the system is said to be under diffusion control.On the other hand, the system will be under kinetic control  if the time constant associated with the kinetics is relatively slower than that of diffusion [1,24].

Diffusion in a Thin
Film.Diffusion occurs in a thin film also, for example, triiodide diffusion in the electrolyte solution of DSSCs.Moreover, diffusion can be coupled with reaction such as the electron diffusion-recombination at the PE of DSSCs.Impedance of such diffusion is known as finite-length diffusion impedance.The impedance of the diffusion and recombination or diffusion and coupled reaction can be modeled as a finite-length transmission line (FTL) composed of distributed elements   ,   , and   as shown in Figures 8(a) and 8(b), where   is given by [26] In thin film diffusion, the diffusion layer is bounded and the impedance at lower frequencies no longer obeys the equation for semi-infinite Warburg diffusion [1,2,4].Professor Bisquert has modeled various aspects of diffusion of particles with diffusion coefficient  in a thin film of thickness , where the characteristic frequency   is [26] In a reflecting boundary condition, electrons, being injected at the interface between a conducting substrate and a porous semiconductor film, diffuse through the film to the outer edge of the film where electron transport is blocked.This diffusion phenomenon can be modeled as a FTL with short-circuit at the terminus similar to that in Figure 8(a), however, without   as the diffusion is not coupled with reaction.On the other hand, in an absorbing boundary condition, electrons are injected at - junction and are collected at the outer edge of the neutral  region of a semiconductor.The diffusion process can be modelled as a FTL with open-circuit at the terminus similar to that in Figure 8(b), of course, without   .The diffusion impedance ( , ) for a reflecting boundary condition is expressed as [26] where   (=   ) and   (= 1/    ) are the diffusion resistance and characteristic frequency of diffusion, respectively.Complex plane plot of this impedance shows a straight line with 45 ∘ at high frequency and then vertically goes up at the low frequency (Figure 8(c)).The high and the low frequency regions clearly show two distinct features separated by the characteristic frequency   .When  ≫   , the system behaves as a semi-infinite and (29) coincides with (24) as [26]  , () =   √    .
At the low frequency region, the impedance becomes [26]  , ( For absorbing boundary condition, the diffusion impedance ( , ) can be expressed as [26] The impedance in complex plane plot appears as an arc at the low frequency region and a straight line with 45 ∘ showing semi-infinite behavior at high frequency region that follows (30) as shown in Figure 8(d).
The impedance of the diffusion and recombination for the reflective boundary condition ( dr, ) is expressed as [26]  dr, ( and the impedance for the absorbing boundary condition ( dr, ) is expressed as [26]  dr, ( where   and   are the diffusion resistance and characteristic frequency for diffusion, respectively, as in ( 29) and (32).The additional terms   and   are the resistance corresponding to homogeneous reactions and the characteristic frequency of the reaction, respectively.Equations ( 33) and (34) have three independent parameters, for example,   ,   , and   .The relation among the physicochemical parameters is expressed as [26] International Journal of Photoenergy where  and   are the film thickness and the diffusion length, respectively.Comparing ( 28) and ( 35), one can write Figure 8(e) shows EIS spectra for impedance of diffusionreaction with reflective boundary condition in complex plane plot for different ratio of   /  .When   is very large (red circles, Figure 8(e)), (33) reduces to (30) of simple diffusion.
In this case, the reaction resistor   in the transmission line model (Figure 8(a)) is open circuit.For a finite   , the impedance takes two different shapes depending on the quotient of (35).If   >   (blue squares, Figure 8(e)), the impedance at high frequency region ( ≫   ) follows (30) and at the low frequency region ( ≪   ) the expression is Thus, the complex plane plot of the impedance has a small Warburg part at high frequency and a large arc at low frequency.In this case, the dc resistance is expressed as When   <   (green triangles, in the inset of Figure 8(e)), (33) gives the expression where the reaction time is shorter than the time for diffusion across the layer (  ≫   ).This is the case when diffusing species are lost before they reach the outer edge of the film.The model corresponding to (39) is called Gerischer's impedance and the dc resistance has the form If   >   (blue squares, Figure 8(f)), (34) approximates to (32); however, the dc resistance is slightly less than that of the case for very large value of   due to additional contribution of   's as in Figure 8(d).When   <   (green triangles, Figure 8(f)), (34) reduces to Gerischer's impedance of (39) and the dc resistance of the impedance is given by (40).

EIS Spectra of DSSCs
The charge transfer kinetics, involved in working DSSCs based on liquid electrolyte containing I 3 − /I − redox couple, are shown in Figure 9(a) with plausible time constants [19,27,28].Within the frequency range of EIS measurement, several time constants are well dispersed in the frequency domain and they give rise to three distinct semicircles in complex plane plot (Figure 9(b)) or three distinct peaks in Bode plot (Figure 9(c)) of EIS of a DSSC at a certain steadystate, at around open-circuit voltage ( oc ) under illumination or at high potential under dark, attained by applying a voltage and illumination.These semicircles in the EIS spectra have been assigned to corresponding charge transfer processes by means of theoretical and experimental approach [5,6,12,13,29].Among the three semicircles of the complex plane plot (Figure 9(b)), in the order of decreasing frequency, the first semicircle corresponds to the charge transfer processes at the Pt/electrolyte and uncovered FTO/electrolyte interfaces with a characteristic frequency  CE , the second or middle semicircle corresponds to the electron diffusion in the TiO 2 film and electron back reaction with oxidized redox species at the TiO 2 /electrolyte interface, and the third semicircle at the low frequency region corresponds to the diffusion of I 3 − in the electrolyte solution with a characteristic frequency   .The characteristic frequency for electron transport or diffusion (  ) appears at the high frequency region of the middle semicircle while the peak frequency (  ) of that semicircle corresponds to the electron back reaction.Similarly, the Bode plots (Figure 9(c)) show all characteristic frequencies except   , which may appear as a break point at the high frequency limit of second semicircle in complex plane plot at certain steady-states but not in Bode plot.The above description is consistent with the time constants shown in Figure 9(a).
Several research groups have already demonstrated systematic approach to characterize EIS of DSSCs [7,12,13,30].Determination of physical parameters from EIS spectra of DSSCs is often done by fitting the spectra to an equivalent circuit.The most widely used equivalent circuit of the complete DSSCs is a transmission line model as shown in Figure 10(a), where  ct is the charge transfer resistance of the charge recombination process at the TiO 2 /I 3 − in electrolyte;   is the chemical capacitance of the TiO 2 film;   is the transport resistance of electrons in TiO 2 film;   is the Warburg element showing the Nernst diffusion of I 3 − in electrolyte;  Pt and  Pt are the charge transfer resistance and double-layer capacitance at the Pt CE;  TCO and  TCO are the charge transfer resistance and the corresponding doublelayer capacitance at exposed transparent conducting oxide (TCO)/electrolyte interface;  CO and  CO are the resistance and the capacitance at TCO/TiO 2 contact;   is the series resistance; and  is the thickness of the mesoscopic TiO 2 film [7].At high illumination the equivalent circuit may be simplified to Figure 10(b).In addition to selecting an appropriate equivalent circuit, one must be able to estimate the parameters to a good approximation from the EIS spectra to initiate the fitting on a program that usually comes with every EIS workstation.Adachi et al. showed how to determine the parameters relating to charge (electrons and I 3 − ) transport in a DSSC from EIS spectra [6].The EIS spectra of DSSCs do not necessarily show three distinct arcs in the complex plane plot or three peaks in Bode plot; however, proper inspection of the experimental data may help to extract the important parameters efficiently.Even though the charge transfer processes in a working DSSCs are more complicated than the above description, we will mainly discuss most significant processes and how the impedance of those individual processes shapes the EIS spectra of complete DSSCs.

Ohmic Series Resistance.
The sheet resistance of electrode substrate and the resistance of electrolyte solution are the main contributor to the Ohmic series resistance ( OS ) in DSSCs.The impedance ( OS ) for the  OS is (3) (4) (5) (2)  Typical impedance spectra of a DSSC presented in complex plane and (c) Bode imaginary (blue solid line) and phase (green solid line) plot.In the order of decreasing frequency, the characteristic frequencies  CE ,   ,   , and   correspond to the charge transfer processes at the Pt/electrolyte interface, electron diffusion in the TiO 2 film, electron back reaction with oxidized redox species in the electrolyte, and diffusion of redox species in the electrolyte solution, respectively.
where  is the ideal gas constant,  is the Faraday constant,  is the temperature, and  is the number of electrons involved in the reaction.The charge transfer process at the CE can be modeled as a - parallel circuit and the corresponding impedance ( Pt ) can be expressed in terms of CPE as So, the characteristic frequency of the charge transfer process ( CE ) can be calculated as and the equivalent capacitance of  Pt ( Pt ) can be calculated as  8(a).Thus, the impedance of diffusion and recombination of electrons at the PE of DSSCs can be expressed by (33).In practice, the distributed capacitance   is replaced with distributed CPE to account for the nonideality in the diffusion-recombination processes.In this case, the characteristic frequency   can be expressed in terms of CPE as [32,33]   = 1 Similarly, the characteristic frequency   can be written as [32,33]   = 1 and the impedance can be expressed as [32,33] where   is the electron transport or diffusion resistance,   is the electron recombination resistance, and   and   are the CPE prefactor and index, respectively.The chemical capacitance (  ) of the TiO 2 film permeated with electrolyte can be calculated from CPE as  2. The blue circle shows the EIS spectra simulated on Zview software using the same parameters and according to the equivalent circuit as shown in the inset, where DX is the extended element 11: Bisquert number 2 that corresponds to the impedance of the diffusion-recombination process at the PE of DSSCs ( PE ).
According to (35) and ( 46), (48) can be rearranged as (50) 3.4.Diffusion of  3 − in the Electrolyte Solution.In the electrolyte solution, concentration of I − is much higher than that of I 3 − .As a result, only I 3 − contributes to the diffusion impedance that appears at the low frequency region [6,12,13].The impedance of I 3 − diffusion can be modelled as FTL with short-circuit terminus and without   (Figure 8(b)) and the corresponding finite-length Warburg impedance (  I ) can be expressed as where   I is the diffusion resistance,  I is the diffusion coefficient of I 3 − , and  is the diffusion length, which is considered as half of the thickness of the electrolyte film [12].Equation (51) is the same as (32), which expresses the impedance of finite-length diffusion with absorbing boundary condition provided that   =  I / 2 where   is the characteristic frequency of the diffusion.The frequency maxima ( max ) are related to   as  max = 2.5  .(52)

Constructing EIS Spectra of Complete
Figure 11 shows complex plane plot for the impedance of a DSSC showing individual components calculated through (52) using the parameters obtained from an EIS spectrum of a DSSC with N719 loaded TiO 2 as a PE, I − /I 3 − based liquid electrolyte, and a platinized CE measured at open-circuit voltage under 1 sun condition (Table 2).To compare the EIS spectrum calculated on spreadsheet (green solid line, Figure 11) with that obtained by commercially available software, EIS spectrum of DSSC (blue circle, Figure 11) was also calculated on Zview software (Zview version 3.1, Scribner Associates Inc., USA) according to the equivalent circuit shown in the inset.It is found that both spreadsheet calculation and Zview simulation generate exactly the same EIS spectrum of DSSC.

Conclusions
Spreadsheet calculation can successfully simulate EIS spectra of DSSCs.Calculation of EIS on spreadsheet allows one to get overall idea of how EIS spectra of DSSCs evolve from impedance response of individual components of DSSCs and how the properties of the EIS spectra are related to each other.Any kind of EIS spectra can be calculated on spreadsheet using the built-in function available in the spreadsheet provided that the corresponding impedance expression is known.This review should help one to learn EIS of DSSCs

Figure 1 :
Figure 1: Screenshots of the spreadsheet calculation of impedance of a capacitor ( dl ) with capacitance of 100 F at frequencies 10 mHz and 100 kHz showing formulas and corresponding results in MS Excel.

Figure 2 :
Figure 2: (a) Complex plane plot for the impedance corresponding to the equivalent circuit as shown in the inset with   = 50 Ω and  dl = 100 F and (b) Bode magnitude and phase plot of the impedance.(c) Complex plane plot for the impedance corresponding to the simplified Randle's circuit with   = 50 Ω,  ct = 100 Ω, and  dl = 100 F as shown in the inset; (d) Bode magnitude and phase; and (e) Bode imaginary and phase plot of the impedance.

Figure 3 :
Figure 3: (a) Complex plane plot for the impedance corresponding to the simplified Randle's circuit with   = 50 Ω,  dl = 100 F, and different values of  ct as mentioned.(b) Bode magnitude and (c) phase plot for the impedance.

Figure 4 :
Figure 4: (a) Complex plane plot for the impedance corresponding to the simplified Randle's circuit with   = 50 Ω,  ct = 100 Ω and different values of  dl as mentioned.(b) Bode magnitude and (c) phase plot for the impedance.

Figure 5 :
Figure 5: Complex plane plots for the impedances that show inductive effects at (a) the high frequency and (b) the low frequency regions.Inset shows the corresponding equivalent circuits with   = 10 Ω,  ct = 100 Ω,  1 = 160 Ω,  dl = 100 F, and  1 = 10 H for (a) and 100 H for (b).

Figure 6 :
Figure 6: (a) Equivalent circuit with   as CPE.(b) Complex plane, (c) Bode magnitude, and (d) phase plot for the impedance corresponding to the equivalent circuit with   = 50 Ω,  ct = 100 Ω, and   = 100 F⋅s −1 and different values of CPE index  as mentioned.

Figure 7 :
Figure 7: (a) Semi-infinite transmission line depicting diffusion process.(b) Complex plane and (c) Bode plot for the Warburg diffusion where the coefficient  = 150 Ω s −0.5 .(d) Complex plane and (e) Bode plot for the impedance corresponding to Randle's circuit with   = 50 Ω,  ct = 100 Ω, and  dl = 100 F and the Warburg coefficient  = 150 Ω s −0.5 .Inset (d) shows Randle's circuit and magnitude and phase plot for Randle's circuit.

Figure 8 :
Figure 8: Finite-length transmission line models of diffusion-reaction impedance with (a) reflective and (b) absorbing boundary condition.Complex plane plots of the impedance model for diffusion with the (c) reflective boundary condition and (d) the absorbing boundary condition.Complex plane plot of the impedance model for diffusion coupled with a homogeneous reaction with the (e) reflective boundary condition (inset shows magnified view of the high frequency region of the plot) and (f) the absorbing boundary condition.

Figure 8 (
Figure 8(f) shows the complex plane plot of the impedance for diffusion-reaction with the absorbing boundary condition for different cases of   /  .For a very large value of   (red circles, Figure 8(f)), (34) turns into (32) of simple diffusion as in Figure 8(d).The dc resistance of the impedance equals   .If   >   (blue squares, Figure8(f)), (34) approximates to(32); however, the dc resistance is slightly less than that of the case for very large value of   due to additional contribution of   's as in Figure8(d).When   <   (green triangles, Figure8(f)), (34) reduces to Gerischer's impedance of (39) and the dc resistance of the impedance is given by (40).

Figure 9 :
Figure 9: (a) Charge transfer kinetics involved in dye-sensitized solar cells where dark arrow shows loss mechanism: (1) injection of electrons, (2) diffusion of electrons in the TiO 2 , (3) regeneration of dye, (4) regeneration of redox mediator (I 3 − /I − ), (5) diffusion of I 3 − , (6) diffusion of I − , (7) back reaction of TiO 2 conduction band electrons with I 3 − , (8) recombination of electrons with oxidized dye, and (9) back reaction of electrons from FTO to I 3 − .(b)Typical impedance spectra of a DSSC presented in complex plane and (c) Bode imaginary (blue solid line) and phase (green solid line) plot.In the order of decreasing frequency, the characteristic frequencies  CE ,   ,   , and   correspond to the charge transfer processes at the Pt/electrolyte interface, electron diffusion in the TiO 2 film, electron back reaction with oxidized redox species in the electrolyte, and diffusion of redox species in the electrolyte solution, respectively.

3. 2 .
Charge Transfer at the CE.The charge transfer resistance ( Pt ) at the Pt CE is associated with the redox reaction involving I − and I 3 − .The exchange current density ( 0 ) of the reaction is related to  Pt by Buttler-Volmer equation as[13]

Figure 11 :
Figure 11: Complex plane plot for the impedance of a DSSC showing calculated impedance of individual components and complete DSSC using parameters as summarized in Table2.The blue circle shows the EIS spectra simulated on Zview software using the same parameters and according to the equivalent circuit as shown in the inset, where DX is the extended element 11: Bisquert number 2 that corresponds to the impedance of the diffusion-recombination process at the PE of DSSCs ( PE ).

Table 1 :
Basic electrical elements and their current-voltage relation.