Abstract
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.
1. Introduction
Impedance spectroscopy is a powerful method for characterizing the electrical properties of materials and their interfaces [1–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–17]. Typically, a dye-sensitized solar cell (DSSC) is composed of a ruthenium dye loaded mesoporous film of nanocrystalline TiO2 on fluorine-doped tin oxide (FTO) glass substrate as photoelectrode (PE), an iodide/triiodide based redox electrolyte solution, and a Pt coated FTO glass substrate as counter electrode (CE) [14, 15, 18–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 TiO2 and ionic diffusion in the electrolyte solution; (ii) charge transfer due to electron back reaction at the FTO/electrolyte interface and recombination at the TiO2/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 TiO2 [2, 7, 15, 21–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. 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 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 , 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 complex number, a small-amplitude AC voltage can be described as and response to this potential is the AC current , where is the phase difference between and and . Therefore, (2) can be written as [1, 4] Again, (4) can be rewritten in terms of magnitude as [1, 4] Applying Euler’s relationship and replacing with , (5) can be expressed as [1, 4] In general, impedance is expressed as [1, 4] or more simply as [1, 4] where and 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 and are plotted against . Sometimes, it is helpful to plot against [1, 2].
In the frequency domain, current-voltage relations can be rearranged as (2). If a purely sinusoidal voltage is applied across a resistor with resistance then the current that flows through the resistor will be , which can be written as . 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 , where and . The above expression for the current passing through the capacitor can be written as or , where . The impedance of the capacitor, , is thus [1] where or in complex notation is the reactance of a capacitor and 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 . 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.
2.2. Equivalent Circuit of Some Electrochemical Systems and Their Impedance
2.2.1. 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 is the double layer capacitance then the total impedance of the system becomes where and are the impedance for and , 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 , 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 Ω, while the capacitive impedance contributes to the negative imaginary part of the impedance. At the low frequency limit the capacitive impedance is so large that the total impedance is infinity. Therefore, the dc resistance, , 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)), and are plotted against . The plot of impedance (red circle) versus frequency has a breakpoint, which corresponds to the characteristic frequency or characteristic time constant 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.
(a)
(b)
(c)
(d)
(e)
2.2.2. Nonpolarizable Electrode in Contact with Electrolyte
If the electrode is nonpolarizable, then the system can be modeled by introducing a resistance parallel to the capacitance 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 and a double layer capacitance . The impedance of the system can be written as where and is the impedance of the parallel combination of the and the .
Thus, (16) can be written in terms of reactance as [1] Here, and 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 . 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 . So, the dc resistance of this system is . It can be noticed from (17) that the maximum of the occurs at , which corresponds to the characteristic frequency of the charge transfer process .
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 : and the second breakpoint corresponds to the time constant : Here, the frequency in the Bode magnitude plot (red circle, Figure 2(d)) can be calculated from (18) as Hz. On the other hand, is calculated to be 15.92 Hz for and μF. The Bode phase plot (blue square, Figure 2(d)) has a maximum at around the frequency 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 , 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 . The semicircle progressively increased as the value of increased from 50 to 100 while 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 with the increase of . On the other hand, the complex plane plot (Figure 4(a)) remained unchanged for a fixed value of and different values of . Thus, changes as the changes while it remains fixed for all values of . For both cases, the Bode magnitude and phase plots depicted in Figures 4(b) and 4(c) clearly show that shifts towards the low frequency edge for increasing either or . It is to be noted here that phase angle at the maxima decreases with the decrease of while it is the same for different values of .
(a)
(b)
(c)
(a)
(b)
(c)
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].
(a)
(b)
2.2.4. Constant Phase Element
In equivalent circuit model of an electrochemical system, the capacitance is often replaced by a constant phase element (CPE) to account for the deviation of the 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 is expressed as [2] 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 of the electrochemical interface corresponding to the parallel combination of and of Figure 6(a) can be calculated by comparing (22) with (19) as
(a)
(b)
(c)
(d)
2.2.5. 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 , (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.
(a)
(b)
(c)
(d)
(e)
2.2.6. 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 of the faradic or charge transfer kinetics is too fast compared to the time constant 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].
2.2.7. 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 .
(a)
(b)
(c)
(d)
(e)
(f)
The diffusion impedance for a reflecting boundary condition is expressed as [26] where and 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 is expressed as [26] and the impedance for the absorbing boundary condition is expressed as [26] 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] 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 diffusion-reaction 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 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, 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).
3. EIS Spectra of DSSCs
The charge transfer kinetics, involved in working DSSCs based on liquid electrolyte containing 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 steady-state, at around open-circuit voltage 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 , the second or middle semicircle corresponds to the electron diffusion in the TiO2 film and electron back reaction with oxidized redox species at the TiO2/electrolyte interface, and the third semicircle at the low frequency region corresponds to the diffusion of 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).
(a)
(b)
(c)
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 is the charge transfer resistance of the charge recombination process at the in electrolyte; is the chemical capacitance of the TiO2 film; is the transport resistance of electrons in TiO2 film; is the Warburg element showing the Nernst diffusion of in electrolyte; and are the charge transfer resistance and double-layer capacitance at the Pt CE; and are the charge transfer resistance and the corresponding double-layer capacitance at exposed transparent conducting oxide (TCO)/electrolyte interface; and are the resistance and the capacitance at TCO/TiO2 contact; is the series resistance; and is the thickness of the mesoscopic TiO2 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 ) 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.
(a)
(b)
3.1. Ohmic Series Resistance
The sheet resistance of electrode substrate and the resistance of electrolyte solution are the main contributor to the Ohmic series resistance in DSSCs. The impedance for the is
3.2. Charge Transfer at the CE
The charge transfer resistance at the Pt CE is associated with the redox reaction involving and . The exchange current density of the reaction is related to by Buttler-Volmer equation as [13] 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 can be expressed in terms of CPE as So, the characteristic frequency of the charge transfer process can be calculated as and the equivalent capacitance of can be calculated as
3.3. Electron Diffusion and Recombination at the PE
In DSSCs, electron transport through diffusion in the TiO2 is coupled with electron back reaction, generally termed as recombination, at the TiO2/electrolyte interface. The impedance of diffusion and recombination of electrons at the PE of DSSCs has been extensively studied by several research groups [5, 12, 26, 31]. Impedance of this diffusion-recombination process appears in the middle semicircle of EIS spectra of DSSCs (Figure 9(b)) with characteristic frequencies and . The PE permeated with liquid electrolyte clearly resembles the electrochemical system with reflecting boundary as shown in Figure 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] Similarly, the characteristic frequency can be written as [32, 33] 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 TiO2 film permeated with electrolyte can be calculated from CPE as According to (35) and (46), (48) can be rearranged as
3.4. Diffusion of in the Electrolyte Solution
In the electrolyte solution, concentration of is much higher than that of . As a result, only contributes to the diffusion impedance that appears at the low frequency region [6, 12, 13]. The impedance of diffusion can be modelled as FTL with short-circuit terminus and without (Figure 8(b)) and the corresponding finite-length Warburg impedance can be expressed as where is the diffusion resistance, is the diffusion coefficient of , 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 where is the characteristic frequency of the diffusion. The frequency maxima are related to as .
3.5. Constructing EIS Spectra of Complete DSSCs
According to Figures 9(a) and 10(a), a simple electrical equivalent circuit of DSSCs can be constructed by combining the elements that are involved in the impedances , , , and [7, 13]. Thus, the impedance of complete DSSCs can be calculated by summing up (41), (43), (50), and (51) as 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 TiO2 as a PE, 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.
4. 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 as well as to develop a basic understanding of EIS in general from scratch.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by National Research Foundation of Korea (NRF) Grants (NRF-2009-C1AAA001-2009-0093168 and 2012-014844) funded by the Ministry of Education, Science and Technology (MEST). Also, this work was partially supported by the NRF Grant 2011-0024237 funded by MEST through the Basic Science Research Program.
Supplementary Materials
The supplementary material is an MS Excel 2010 file that demonstrates spreadsheet calculation of EIS spectra. In spreadsheet, a complex number can be constructed using the built-in function and the number can be operated with all the basic mathematical operators available in the spreadsheet as functions for complex numbers. Based on the above concept, all the EIS plots discussed in the present paper are calculated on spreadsheets of this workbook.