Abstract

We propose a numerical procedure consisting of a simplified physical model and a numerical method with the aim of optimizing the performance parameters of dye-sensitized solar cells (DSSCs). We calculate the real rate of absorbed photons (in the dye spectral range) by introducing a factor in order to simplify the light absorption and reflection on TCO electrode. We consider the electrical transport to be purely diffusive and the recombination process only to occur between electrons from the TiO2 conduction band and anions from the electrolyte. The used numerical method permits solving the system of differential equations resulting from the physical model. We apply the proposed numerical procedure on a classical DSSC based on Ruthenium dye in order to validate it. For this, we simulate the J-V characteristics and calculate the main parameters: short-circuit current density , open circuit voltage , fill factor FF, and power conversion efficiency . We analyze the influence of the nature of semiconductor (TiO2) and dye and also the influence of different technological parameters on the performance parameters of DSSCs. The obtained results show that the proposed numerical procedure is suitable for developing a numerical simulation platform for improving the DSSCs performance by choosing the optimal parameters.

1. Introduction

The technology and materials used for the third generation solar cells give the opportunity to obtain cells with high efficiency [15]. The solar cells based on dye-sensitized nanostructure with mesoporous metal oxides (DSSCs) have attracted considerable attention since the work of O’Regan and Grätzel [6], their manufacturing being environment-friendly and energy-efficient [1, 7]. Up to now, certified efficiencies over 10% under standard conditions or even higher (12,4%) at the laboratory scale were reported [710]. Based on the low cost of materials and the simplicity of fabrication process, DSSC can have lower fabrication costs than conventional silicon-based solar cells. Taking into account this advantage, the improvement of DSSC parameters for making them widely used appears as a strong necessity.

Further optimization of the DSSC parameters requires a better correlation between interrelated processes of transport and accumulation of electrons in the mesoporous oxide phase and recombination of electrons with electron acceptors [11]. In order to understand the different processes governing the DSSC’s mode of operation and to enhance the DSSCs performance, modeling of processes and numerical simulation of the cells were carried out [1223]. The thickness, morphology (particle diameter and porosity) for the TiO2 layer, electron mobility and recombination rate (electron lifetime), absorption spectrum of the dye, thickness and material quality of transparent conductive oxide (TCO) layers and the used electrolyte determine the characteristics [17, 18, 2022, 24]. Consequently, the extracted parameters from simulation, the short-circuit current density , the open circuit voltage , the fill factor FF, and the power conversion efficiency were analyzed.

In the present paper, we propose a numerical procedure for optimizing the DSSCs consisting of a simplified physical model and a numerical method capable of solving the system of differential equations resulted from the physical model. We applied this proposed procedure on a DSSC with Ru dye in order to show its validity. For this, the current density-voltage characteristics were simulated, and the main parameters  , , FF, and were obtained and discussed. The influence of the thickness of TiO2 film and electron lifetime on the cell performance parameters was analyzed aiming to be used in the cell manufacturing.

2. Physical Model and Numerical Method

We use a classical DSSC which is presented in Figure 1. This device contains two electrodes made of TCO glass, commonly used being fluorine-doped tin oxide (FTO). The illuminated electrode is coated with a nanoporous TiO2 layer which is also coated with a monolayer of the Ru-complex dye. On the counter electrode of TCO glass, a thin layer of Pt (~5 nm) which acts as a catalyst for the redox reaction is deposited [25]. The space between the two electrodes is filled with an electrolyte (propyl methyl imidazolium iodide) containing an iodide/triiodide () redox couple.

In a DSSC under visible light, an electron from the dye molecule is excited, then injected into the conduction band of TiO2 nanostructured semiconductor, and finally collected by the TCO glass electrode and transported in the external load. The positive charge (dye cations) is reduced by receiving an electron from the iodide ion which will be subsequently regenerated by reducing triiodide ions () at the platinized counter electrode.

By using the nanostructured TiO2, the depletion process is less pronounced, so that under reverse bias the electric field is negligible [14, 26]. The origin of the photovoltage is explained by a built-in potential barrier between TiO2 semiconductor and the TCO electrode and/or by changing the position of Fermi level due to the electron injection [26]. The most probable recombination process takes place between electrons from TiO2 and ions [18, 19].

In the ideal case, the rate of the absorbed photons in the volume unit can be written as where describes the position inside the TiO2 film (including the dye) with a thickness , is the wavelength of the radiation, limits and are imposed by the absorption characteristics of the dye through the absorption coefficient , and is the spectral incident photon flux density.

In the real case, the rate of the absorbed photons in the volume unit is reduced to due to the light absorption in TCO glass electrodes and to the reflection onto [16]. We describe the contributions of these processes by a factor aiming to propose a simplified model to be used for designing and optimizing any type of DSSCs. Therefore, is given by

The rate of electron injection should take into account the injection efficiency [11], so that

In our simplified model we consider the electrical transport to be purely diffusive, and therefore it is described by a diffusive transport equation for the electrical current density [17, 19]. We neglect the internal electric field generated by unbalanced local charge because of its very low value [12, 17]. The electron density in the conduction band of the TiO2 layer is described by the continuity equation in quasi-stationary regime. In this case, the system of equations has the following form: where is the elementary charge, describes the position inside the TiO2 film , is the recombination rate, and is the electron diffusion coefficient.

We neglect the trapping/detrapping processes, so that in the continuity equation (see (4)) the corresponding terms do not appear [15]. Also we assume that the recombination process takes place only between electrons from conduction band of TiO2 layer and anions present in the electrolyte [18]. Consequently, the recombination rate is proportional to , so that where is the electron density at equilibrium (in the dark) and is the electron lifetime.

The boundary conditions are the following.(1)Assuming that electrons are extracted by application of an external voltage , the boundary condition at is where is the density of states in the conduction band of TiO2 semiconductor, given by the formula where is the effective mass of the electron, is Boltzmann constant, is Planck constant, is the cell temperature, and is the energy between the conduction band edge and the Fermi quasi-level .(2)Assuming that all the electrons are collected at (there is no recombination process inside TiO2 film meaning that TiO2 film has a good quality from the technological point of view), the boundary condition at is

The open circuit voltage is derived from the condition and is given by the expression where (0) is the rate of electron injection at , while the short-circuit current density is obtained for the condition .

The system of differential equations (4) and (5) together with boundary conditions (see (7) and (9)) has no enough conditions at but has the boundary condition at (see (9)). This problem is known as “the boundary value problem,” being defined by differential equations in which some conditions are specified at the initial point while the others are specified at the end point. In order to solve this system, we used a numerical method based on the “shooting method” [27] which generates different values for the initial condition at until the final condition at is satisfied. We used this method in a program developed in Mathcad.

Consequently, the numerical procedure we propose for designing and optimizing any type of DSSCs considers the essential processes which occur in a real DSSC and its technological parameters and simplifies the others, for example, the contribution of the light absorption and reflection on the TCO electrodes by introducing for obtaining .

3. Results and Discussion

In order to validate our numerical procedure (simplified model and numerical method) we apply it on a DSSC based on Ru dye and simulate the characteristics. For this, we took into account the absorption coefficient of the Ru dye between the spectral limits and where its absorption is very good (~103 M−1 cm−1) [28]. The numerical procedure can be applied on a DSSC based on a different dye by changing the spectral limits and the absorption coefficient. For the incident flux density we considered the AM1.5 G spectrum [29]. The other constants and parameters used in the simulation were taken from the literature [11, 12, 21] and are presented in Table 1.

Figure 2 shows the simulated characteristics obtained for and and different values for the electron lifetime in the conduction band.

The values obtained for , , FF, and using the parameters and and different values for are listed in Table 2. One can see that an increase of with 20% (by using suitable TCO electrodes) produces an increase of both and with about 24%, while and FF remain practically constant.

The simulated values for , , FF, and as a function of the TiO2 film thickness, using , are presented in Table 3. According to this, Figures 3(a) and 3(b) present the thickness dependence of the short-circuit current density and the power conversion efficiency , the TiO2 film thickness being an important technological parameter for designing and optimizing the cell. One can see that both and increase with the film thickness up to maximum values of and , respectively, reached for . By further increasing the film thickness up to 50 μm, they insignificantly decrease (-). Consequently, a thicker layer will not bring any improvement to or (helpful information for design and technological processes). Also, one can observe that the increase of film thickness from 2 to 10 μm does not influence the open circuit voltage (remains constant) and the fill factor FF (decreases with about ).

In Table 4 the values for , , FF, and calculated for different values of the lifetime using are given. Figures 4(a) and 4(b) present the corresponding curves, that is, the lifetime dependence of the short-circuit current and the power conversion efficiency . They show that both and significantly increase when increases up to 50 ms. This means that if other materials instead of TiO2 ( takes values up to 50 ms) are used, the performance parameters of DSSC can be controlled and consequently optimized.

The values for all parameters we obtained from simulation are in good agreement with those reported in the literature for this type of DSSC [14, 19, 21, 25, 3032] which demonstrate that our proposed numerical procedure is valid and versatile. Therefore, it can be used for designing and optimizing DSSCs.

4. Conclusions

In this paper we propose a numerical procedure for optimizing the parameters of any type of DSSC, consisting of a simplified physical model and a numerical method capable of solving the system of differential equations resulted from the model.

We applied this procedure on a classical DSSC based on Ru dye by simulating the characteristic and calculating all main parameters, , , FF, and , and the obtained results are in good agreement with those reported in the literature, which validate our proposed procedure.

For this, we analyzed the influence of the nature of semiconductor and dye and also the influence of different technological parameters on the performance parameters of DSSCs. Thus, for obtaining the real rate of the absorbed photons we simplified the light absorption and reflection of TCO by introducing the factor and we considered the nature of dye by using the specific absorption coefficient . We analyzed the influence of TiO2/semiconductor thickness and electron lifetime in the conduction band on the main cell parameters. We showed that the increase of TiO2 semiconductor thickness over the optimal thickness of 15 μm does not influence the and values ( = 12.601 mA/cm2 and ). We obtained a similar result for the electron lifetime higher than 40–50 ms. The obtained results in good agreement with the literature validate our numerical procedure. The analysis made on this classical DSSC based on TiO2 with Ru dye can be extended to any other DSSCs.

The proposed numerical procedure allows a comprehensive analysis of the performance parameters of DSSCs and permits the development of a numerical simulation platform for designing and optimizing any type of DSSCs.

Acknowledgment

This work was partially supported by the Romanian National Authority for Scientific Research through the CNCS-UEFISCDI under Contract no. PN II-PT-PCCA-9/2012.