Steady-State and Transient Response Analysis of DSSC Based on Electron Diffusion Coefficient and Chemical Capacitance

Article information

J. Electrochem. Sci. Technol. 2024;15(2):276-290
Publication date (electronic) : 2024 February 6
doi : https://doi.org/10.33961/jecst.2023.00808
1Centro de Física Aplicada y Tecnología Avanzada, Universidad Nacional Autónoma de México, Boulevard Juriquilla 3001, Querétaro 76230, México
2Centro de Investigación y Desarrollo Tecnológico en Electroquímica S.C., Parque Tecnológico Querétaro, Sanfandila 76703, Pedro Escobedo, Querétaro, México
*E-mail address: ranged@fata.unam.mx
Received 2023 September 1; Accepted 2023 December 12.

Abstract

A study of the transition from transitory state to steady state in DSSCs based on natural dyes is presented; cochineal was used as dye and Li+, Na+, and K+ were the ions added to the electrolyte. The photocurrent profiles were obtained as a function of time. Several DSSCs were prepared with different cations and their role and the transitory-to-steady transition was determined. A novel hybrid charge carrier source model based on the Heaviside function H(t) and the Lambert-Beer law, was developed and applied to analysis of the transient response of the output photocurrent. Additionally, the maximum effective light absorption coefficient α and the electronic extraction rate κ for each ion were determined: αLi+, Na+, K+ = (0.486, 0.00085, 0.1126) cm−1, and also the electronic extraction rate κext.Li+, Na+, K+ = (1410, 19.07, 19.69) cm s−1. The impedance model using Ficks second law was developed for carrier recombination to characterize the photocurrent.

1. Introduction

Different electrochemical devices have been developed to benefit from the energy and responsiveness of the redox reactions of different compounds. Devices such as organic LEDs have been developed based on the ability of certain polymers to emit light after certain redox reactions at the metal/material-organic interface [1]. In recent works, thin films composed of mixed metal oxides such as vanadium V and tungsten W have been deposited and electrochemically characterized to study electrochemical performances for applications in electrochromic and electrochemical energy storage [2]. Due to the nanostructured shape of such oxides, ion diffusion into and out of the films is promoted, which can facilitate charge storage during polarization. In addition to the fact that nanostructured oxide materials have turned out to be a new energy storage medium due to their electrochemical properties, they have also proved to be useful for generating electrical energy due to the high amount of surface area due to the inherent porosity. For this purpose, solar cells sensitized by pigment (DSSC) or quantum dots (QDSSC) have been developed. It should be noted that in other recent works, the performance of a QDSSC and a supercapacitor (built using electrochromic oxides) have been optimally coupled, resulting in a composite device capable of transforming solar energy, storing charge and generating charge indicators due to the color of the electrochromic material [3]. DSSCs have a lower efficiency than silicon cells, but the materials for their development are inexpensive and environmentally friendly [4,5]. Ruthenium complexes were commonly used as sensitizers in DSSCs leading to highest reported efficiencies. However, due to concerns about their environmental impact and potential toxicity, alternative approaches had been developed using natural organic pigments (NOPs) from leaves, seeds, flowers, insects, etc. [6,7]. Although NOPs-based DSSCs are less efficient, they are environmentally friendly. Platinum is the most efficient catalyst in the counter electrode due to its inherent catalytic properties; however, the cathode is covered with carbon since this also has catalytic properties but at lower cost.

For the reasons mentioned above, a device such as the NOPs-based DSSC and carbon-based electrode offers a promising eco-friendly inexpensive solution, aiming to enhance the sustainability and safety of DSSCs as an alternative to inorganic and fossil fuel solar cells, striving for comparable performance. Studies carried out to determine the performance and efficiency of synthetic pigment-based DSSCs reveal that they depend on several factors such as the molecular structure of pigment [811], the type of ions in the redox couple [12], the grain boundary conditions and the TiO2 film [13,14], etc.

The different techniques used for characterization enhance the different electrochemical parameters involved in the operation of a DSSC. The photoelectrochemical impedance technique (PIT), is one of the most used techniques in DSSCs characterization, as it allows to obtain information on the electrochemical parameters of each electrode and electrolyte of the DSSCs in the illumination state, this is done by imposing a sinusoidal signal in the steady state of illumination on the carriers:: n(x,t) = n + ñ, where n is the steady state electronic concentration and ñ is the sinusoidal excitation of amplitude 10 mV approx. [15]. The technique is mainly τn = 1/krec, based on Ficks laws with a certain recombination term related to the carrier lifetime. The PIT technique also allows determining the electronic concentration in the film by the chemical capacitance parameter, typical of a porous electrode system [16,17]: Cμ = Cμ, cb + Cμ, ss. This is a consequence of the distribution of the trap states and conduction band states from the oxygen vacancies at the grain boundaries (Ti3+).

Fig. 1 shows the electronic distribution in the conduction band in dark stage A and in light stage B. Two different trap state distributions are also shown: The first one with a lower slope, and the second one with a higher slope in the exponential distribution. When a DSSC is operating in dark (distribution A Fig. 1) and light (distribution B Fig. 1), electrons are trapped and released with certain velocities (average trapping velocity: vtrapp and average de-trapping velocity: vdetrapp) until they reach the collector.

Fig. 1

Schematic representation of the trapping-detrapping process for photogenerated electrons and insertion-releasing process of cations occurring at the dyed TiO2 film-electrolyte junction where x = L stands for the thickness film.

Such rate constants will depend on the ability to neutralize the charge in the trap states by the cations (Na+, K+, Li+) in the electrolyte. Therefore, a factor that will determine the average detrapping rate is the ease of the counterion to penetrate the TiO2 film and neutralize the trap states by forming an ionic pair in the pores and immobilizing the charge, which is why the size of the cation has an influence. As can be seen in Fig. 1, another factor that also influences the average rate of detrapping is the distribution of trap states, the average energies or depth of the trap states mc are higher in distribution 1 than in distribution 2, thus it will be easier to release an electron in a semiconductor oxide film that presents a distribution of trap states of type 1. When the trapped electrons are affected by either of the two factors mentioned above, vdetrapp tend to decrease, which causes the carriers to increase the time to reach the collector (TCO) τc.

In order to study separately the factor that represents the neutralization of charge in the trap states through the insertion of cations (Na+, K+, Li+) in the film and the consequent decrease in the rate of detrapping, electrodes with the same particle size and heat treatment at the same temperature have been prepared, therefore it can be considered that in the ratio of electron collection time τc and short circuit current density Jsc, mc remains constant in each electrode [18]:

(1) τcJscα-1

Where a=kBTmc is the slope of the exponential distribution of the trap states given in the interval 0 ≤ α ≤ 1; kBT is the electron energy, which, as can be seen from the distribution of the trap states, conditions the value of the short-circuit current density. On the other hand, in Bisquerts model an approximation is found in the dependence of τc and the electron diffusion coefficient Dn given by [19]:

(2) τc(1Dn)1mc-1

It should be noted that from Eq. (2), it is understood that when Dn is maximum, τc is minimum and Jsc is maximum. Moreover, if mc is constant, the short circuit current will depend on the diffusion coefficient Dn or the adsorption and desorption velocities (va and vd) of the counterions in the film.

The steady-state chemical capacitance is crucial since it is related to the effective electron diffusion coefficient Dn by the generalized Einstein relation. There are reports that the concentration of trap states in a porous TiO2 anatase electrode is of the order of 1015 cm−3 [20]; and there is large distribution throughout the film, both indicating that the effective electron diffusion coefficient Dn(n), is dependent on the electronic concentration in the film, and therefore on its illumination or dark state. As mentioned before, the PIT technique also provides information about the impedance of the electrolyte and the counter electrode, which are very often modeled using Warburg impedance ZW and a circuit CCERCE, respectively [21].

The photoelectrochemical impedance technique, by providing a frequency sweep, imposes a steady state for each frequency that allows obtaining the electrochemical parameters in steady-state illumination. The electrochemical parameters during the transition from the dark to the illuminated state are crucial in the full operation of DSSCs. In that transient, the current or voltage delivered by the cell can be affected by the connected load, due to the impedance coupling related to the maximum power transfer to the load [22]. The impedance coupling between the DSSC and the load can, at times, lead to inefficiencies or even prevent the cell from functioning optimally due to poor connections. To address this issue, our study focuses on understanding and analyzing the transition from the dark state to the illumination state of a NOPs-based DSSC with a carbon-based counter electrode. The transient phenomenon is studied by measuring the photocurrent density during the incidence of square pulses of light applied to the DSSC.

In order to know the electron diffusion coefficients in the light-dark transition and the behavior of the chemical capacitance with ions added to the electrolyte, Ficks second law was applied and adjusted to a system of this nature, obtaining a model for the determination of the diffusion coefficient given by a nonlinear partial equation of second order, which tends to be complex to solve in comparison with the model of Dn = const [23]. To avoid the nonlinearity of the equation, the current dark stage and the illumination stage have been treated separately, assuming different electron diffusion coefficients in both stages.

Likewise, a hybrid carrier source model was considered for t ≥ 0 and t ≤ 0, where t = 0 s the onset of the luminous part of the pulse.

2. Methodology

When a DSSC undergoes an abrupt dark-to-light transition, the performance of the cell is influenced by the electron diffusion coefficient, Dn(n). It depends on the electron density and the characteristics and concentration of the added ions; these produce important variations on the magnitude and shape of the photocurrent profiles. The diffusion coefficients of both, the transient and steady states in dark and illumination conditions, were obtained by Cyclic Voltammetry (CV). The chemical capacitance for each ion in illumination condition, Cμion(CμLi+,CμNa+,CμK+), was first determined; these can be associated with the trap states of the TiO2 film; this information allows to determine the electron density in that state; it is possible to mention that the characteristics and concentration of added ions control in some way affect the electron distribution on the mesoporous surface.

The photoelectrochemical impedance technique was used to determine the diffusion coefficients for each added ion: Dnss,.lightLi+=5.59×10-4cm2s-1,Dnss,.lightNa+=3.46×10-4cm2s-1 and Dnss,.lightK+=1.33×10-4cm2s-1. To determine the dark and transient state diffusion coefficients, an dark-light (8 seconds, total time) profile with an intensity 100mW cm−2 (7 seconds-light) was used for each DSSC allowing to obtain the dark diffusion coefficients for each ion: Dnss,.darkLi+=38.1cm2s-1×Dnss,.darkNa+=9.86cm2s-1 and Dnss,.darkK+=12.41cm2s-1. For the transient state: Dntrans.Li+=1.86×10-5cm2s-1,Dntrans.Na+=2.494×10-5cm2s-1 and Dntrans.K+=6.598×10-3cm2s-1.

2.1 Theoretical considerations

Due to the morphology of the photoelectrode, it is possible to neglect drift currents because the particle size (assuming spherical) of the anatase used is ɛr ≈ 102. This implies that the potential difference is less than the energy of the thermal potential kBT=45oq, where q is the charge. Therefore, the diffusive currents are those that control the system and can be studied under Ficks theory.

A photoelectrode in a DSSC can be considered as a mixed-phase electrode, where a pore flooded with electrolyte represents phase (1) liquid and phase (2) solid interconnected TiO2 anatase molecules, Fig. 2.

Fig. 2

Two-phase mixed electrode model based on the cylindrical pore model for dyed TiO2 film-electrolyte junctions where the main electrical impedances of a DSSC are shown.

In our study, we consider the model of Fig. 2, which shows De Levines cylindrical pore impedance model, flooded from x = L to the bottom of the phase-mixed electrode x = 0 [24]. Additionally, the impedance of both phases in the electrode, interface and film boundaries are observed.

2.2 Modeling the electron diffusion coefficient

As it is known, in a DSSC, electrons are generated through the dye molecules and diffuse throughout the film via the surface trap states and conduction band, moreover, the generated electrons can recombine with the oxidized species at the TiO2/E interface [16,25]. Such processes can be modeled by the electrochemical impedance proposed by Bisquert et al. [15]. For the study of the transition from the dark state to the light state, a hybrid model of the carrier source is given by:

(3) G(x,t)={0t0αI0e-αxt0

Such a model consists of a Heaviside step function H(t) in time for which a constant function in time has been taken for t ≥ 0 in the illumination state, but dependent on the position given by the Lambert-Beer absorption law. Here α represents the effective adsorption coefficient of the pigment molecules distributed on the film of thickness L, and I0 is the irradiance or the coefficient of the effective flux of incident photons. Thus, the model describing the system is given by:

(4) C2t=Dn2C2x2+k(C2-C20)+G(x,t)t

Where C2 represents the concentration of carriers in the solid phase (2), (electrons), C20 is the initial concentration of the same species in the same phase, k is the recombination frequency, and Dn the effective electron diffusion coefficient in the film.

However, during the dark-to-light transition, which is when the electronic charge is injected in an order of femtosecond [26], the recombination can be neglected since it is of the order of 101–102 Hz in this case. Therefore, Eq. (4) can be rewritten in the transient state as:

(5) C2t=Dn2C2x2+G(x,t)t

In reported studies [27], it has been shown that the electron diffusion coefficient is highly dependent on the same electronic concentration C2:

(6) Dn~C2l-1

Where l ≈2–3, is a parameter related to the distribution of the trap states [28]. Then, the model is a nonlinear partial differential equation with no analytical solutions. At Ref. [23], solutions were determined by numerical methods for illum,ination pulses of Eq. (5), but those results do not fit the experimental results. Therefore, in the study of dark-illumination transients proposed in this work, an approximation is used in which Dn is time-dependent rather than electron concentration-dependent:

(7) D(t)={Dnt0Dnt0

Where Dn is the effective dark steady-state diffusion coefficient (DSDC) in t ≤ 0, or when C2 = C20 and Dn is the effective diffusion coefficient in the illuminated state (EDCI) in t ≥ 0 of the transient part. This allows for avoiding the nonlinearity of the model as it is shown in the following equation for the approximate model in the transient:

(8a) C2t=Dn2C2x2t0
(8b) C2t=Dn
(8c) 2C2x2+αI0e-αxt0

Where the carrier source given by Eq. (3) has also been used. On the other hand, the boundary conditions when measuring the photocurrent generated by the DSSC in short circuit will be given by:

(9a) kext.C2(0,t)=DnC2x|x=0in the TCO
(9b) Jn=e-DnC2x|x=L=0in the MPEB/E interface

The first condition implies that in TCO the concentration is only time-dependent and is proportional to the concentration gradient evaluated in x = 0, where the proportionality constant is given by the ratio of the electron extraction rate kext. at the TCO contact.

At x = 0, due to the low charge generated by the NOP used [29], it is considered that there is no electrical or transport impediment for the electrons to reach the collector, this is reflected due to the low screening of the electrons arriving at x = 0 and that there is a maximum in the concentration gradient in that region. Therefore, at this boundary, it can be considered that the charge control is first-order kinetic (Arrhenius type) and consequently the diffusion coefficient in that region will be given by

Dn=kext.C2(0,t)C2x|x=0.

At the other boundary, in x = L, at the mixed-phase electrode boundary (MPEB) and the bulk of the electrolyte was determined with the electron flow completely blocked, thus the concentration gradient is minimal and equal to zero and hence there is no electron diffusion in that zone.

The solution for the model C2 (x, t) is given by:

(10a) C2(x,t)=c3exp(-kext.Dnx)exp(kext.2Dnt)+c4t+c5         for t0
(10b) =n=1Bnsin(nπxL)exp[-{n2π2DnL2t}]-I0αDne-αx+c1x+c2t0

That through the continuity conditions of the carrier density: nt0(0,0) = nt0(0,0) and the current density in the whole film: Jt0(x,0) = Jt0(x,0), it is obtained the parameter relationship with the coefficients: c1, c2, c3 and Bn

c1=-I0e-αLDnc2=I0(1-e-αL)kextc3=c2-I0αDnBn=2πn0L[c3kext.Dnexp(kext.Dnx)-c1-I0Dnexp(-αL)]cos(nπcL)dx

Parameters of model the DSSC for Li+; c4 and c5 are parameters determined by fitting experimental results.

As it can be observed Eq. (10b) for t ≥ 0 is a summation of different exponentials, which could represent the contribution of each of the charge species contributing to the current in phase (2). Each species is identified by its diffusion frequency given by:

(11) n2π2DnL2=ωn

Where ωn is the diffusion frequency of the umpteenth component of the current.

Up to this point, it has been considered in the model that the charge transfer at the bottom, i.e., at the phase (1) or TCO/electrolyte (TCO/E) interface of the pore is zero, however, sometimes the photocurrent measured at the TCO (x = 0) is affected by the charges accumulated at the bottom of the pores of the mixed electrode and in the steady state can be taken into account by means of the impedance ZB of Fig. 2. When this area is exposed directly with the electrolyte both recombination and charge accumulation can arise [13]. Since both processes demand charge in that region, they can modify the photocurrent profile measured in the transient.

In these cases, the model is modified if it is initially assumed that a quantity of carriers: δ(x)[C2|t 0(x = 0, t = 0)], is injected into the TCO from the sensitized TiO2 layers at the TCO/TiO2 + dye interface; where C2|t 0 indicates that the evaluation is done in part of the solution when t ≤ 0. Once the charge is injected into the TCO, some carriers tend to return to the mixed-phase electrode by a new reinjection at the TCO/E interface but in a different direction (from the TCO to the electrolyte). As discussed above, the movement of carriers in the opposite direction is due to the accumulation of ions of different signs at the TCO/E interface and may cease until their steady state impedance ZB in direct current is obtained. Since charge neutrality tends to be preserved in the film, it is equivalent to thinking that in x = L, a number of holes has been injected equal to the quantity of electrons injected in x = 0, and that the concentration in the same region is zero. One way to guarantee this boundary condition is to establish a charge of equal magnitude, but different sign at x = −L, the other condition will be given by the blocking of electrons at x = L. Considering that the same model of Eq. (8) for t ≥ 0 and the initial condition mentioned above we have that the electron concentration will be given by:

(12) C2(x,t)=N4πDhte-(x-L)24Dht-N4πDhte-(x+L)24Dht+c1x+c2-I0αDne-αx

Since the measured current density is given by: J=qDnC2x|x=0, where q is the electron charge. In addition, it is noted that the current density profile will depend on the hole diffusion coefficient Dh and the electron diffusion coefficient Dn of the mixed phase electrode:

(13) J=-DnqNL2DhπDht3exp(-L24Dht)+c1+qI0for t0

For the part of t ≤ 0 the solution will be given in the same way by Eq. (10) respectively.

3. Experimental

3.1 DSSCs construction

The semiconductor oxide film was prepared based on a clay of TiO2 anatase at 99.99% particle size < 10 nm (Sigma-Aldrich) and 99.99% glacial acetic acid (Sigma-Aldrich), the film was deposited by the doctor-Blade technique on a conductive FTO glass from 8-10Ωsquare (Sigma-Aldrich). The film received a heat treatment of 35 min at 450°C under normal atmospheric conditions. The thickness achieved in the TiO2 films is approximately 34 μm, measured with a profilometer 3D Bruker, Contour GT In Motion. The dye used was “Dactylopius coccus” known as “cochineal”, due to its low toxicity, its large amount of OH groups, and its abundance in Mexico. Its preparation was based on a solution of 0.5 g of cochineal powder from the Sonora-Mexico region, adding 20 mL of distilled water, 10 mL of 99.99% acetyl acetone (Sigma-Aldrich), and 0.25 g of alum (alumina sulfate and potash). It should be noted that a gap of 53 μm was established between the electrode and the counter electrode. The latter was prepared by a clay of 0.15 g of activated carbon ground and sieved to 74 μm, 0.13 g of polyurethane resin and 0.15 g of acetyl-acetone. The film was deposited by the same technique of roto-etching and a thickness of 23 μm.

Complementarily, the electrolyte employed for filled the three studied DSSCs was a 5:95% v/v solution of 3-methoxypropionitrile (Sigma Aldrich, 99%) and acetonitrile (Sigma Aldrich, 99%) containing 1 mM of LiI, KI, and NaI, in the presence of I2 to promote the formation of triiodide ions I3-.

3.2 Experimental model

As it has been determined, the behavior of the issued photocurrent depends on the diffusion coefficient Dn, which, as stated above, depends on the electronic concentration in the mesoporous given by: n = ncb + nt, the electronic concentration in the conduction band states and the electronic concentration in the trap states. Therefore, it is crucial to determine both the electronic structure and the presence and distribution of trap states in the TiO2 film used as electrodes.

Micro-Raman experiments were performed using a DXR Thermo Scientific dispersive Raman microscope to confirm the anatase structure of the TiO2 films which were covering the surface of the FTO photoanodes, as well as the graphitic (G band) and the amorphous (D band) composition of the carbon-based counter electrodes (see Fig. 3).

Fig. 3

Micro-Raman spectra obtained for (a) Photoelectrode spectra of TiO2 anatase film 34 μm thick on FTO. (b) Spectra of activated carbon film counter electrode and polyurethane resin as binder, thickness 24 μm.

In order to gain an understanding on the presence and distribution of trap sites in the anatase TiO2 films, cyclic voltammetry (CV) experiments were performed using an Epsilon BAS potentiostat-galvanostat which was connected with a three-electrode cell (see Fig. 4) where a silver (Ag) and a platinum (Pt) wire were used as pseudo-reference and counter-electrodes, respectively. Complementarily, FTO plates coated by anatase TiO2 films were employed as working electrodes. The electrochemical cell was filled by a 5:95% v/v solution of 3-methoxypropionitrile (Sigma Aldrich, 99%) and acetonitrile (Sigma Aldrich, 99%) containing 1 mM of LiI, KI, and NaI, in the absence of I2 to inhibit the formation of triiodide ions I3.

Fig. 4

Schematic representation of the three-electrode electrochemical cell employed for carrying out the cyclic voltammetry experiments.

The photoelectrochemical response of the DSSC was studied under visible light illumination (10 mW cm−2 - light power provided by an MR16 GE 12V–50W halogen lamp equipped with a dichroic reflector) by performing photo-electrochemical impedance spectroscopy (PEIS) experiments using an IM6 BAS-Zahner potentiostat-galvanostat where the frequency was scanned from 1 MHz to 50 mHz, whereas D.C. (equals to the open-circuit potential of the photocells) and A.C. signals (amplitude of ±10 mV) were applied under the steady-state of the DSSC.

For the study of the dark-light transient state, a dark-light time-varying amplitude step with an intensity of 100 mW cm−2 and a rising slope of 250 mV s−1 was applied on the electrode side of each DSSC using a 50 W halogen lamp, overlying a hot mirror lens (Edmund Optics model 47303) between the DSSC and the lamp to reduce heat transfer. The photocurrent obtained from each DSSC as a function of time was measured by the Keithley 2400 Source-Meter, using a KUSB-488B data acquisition interface and LabView v. 2014 based Software.

4. Results and Discussion

4.1 Micro-Raman spectra

The composition of both films was determined by the micro-Raman spectra of both electrodes. Fig. 3a shows the spectrum of the photoelectrode, which presents a greater orientation in the crystallographic plane (101) anatase. On the other hand, the counter electrode shows a concentration of graphite, identified by means of the G band, Fig. 3b.

In the mesoporous film, a higher concentration of crystallographic plates (101) has been determined, which is a crucial factor in the injection of electrons and the presence of trap states in the gap. In the activated carbon film, a concentration of graphite is found which allows catalysis at the counter electrode.

4.2 Trap states for each ion

Trap states are a determining factor in the variations of the diffusion coefficient Dn during a transient. Another important factor is the ions present in the electrolyte that can recombine and form ionic pairs with the electrons of the occupied trap states. Therefore, the presence and distribution of these ions were determined. In Fig. 512, the trap states that can be neutralized by different ions Li+, Na+ and K+ using the cyclic voltammetry (CV) technique.

Fig. 5

CV response for a TiO2-based working electrode, Pt auxiliary electrode, Ag reference electrode, and 1 mM LiI solution.

Fig. 6

Distribution of surface trap states neutralized by the Li+ ion on the photoelectrode.

Fig. 7

Distribution of deep trap states neutralized by Li+ ion on the photoelectrode.

Fig. 8

CV response for a TiO2-based working electrode TiO2, Pt auxiliary electrode, Ag reference electrode, and 1 mM NaI solution.

Fig. 9

Distribution of shallow trap states neutralized by Na+ ion on the photoelectrode.

Fig. 10

Distribution of deep trap states neutralized by Na+ ion on the photoelectrode.

Fig. 11

CV response for a TiO2-based working electrode TiO2, Pt auxiliary electrode, Ag reference electrode, and 1 mM KI solution.

Fig. 12

Distribution of deep trap states neutralized by the K+ ion on the photoelectrode.

Fig. 5 shows the voltammogram that allows determining the presence and location of the trap states that can neutralize the Li+ ion; the peaks in the voltagram in the interval [−0.579, −0.335] Potential vs. Fc+|Fc are evidence of the presence of surface trap states on the photoelectrode [30]. The presence of trap states at deeper levels neutralized by the same ion is also determined in the interval [−1.469, −0.986] Potential vs. Fc+|Fc.

Fig. 6 shows the distribution of the surface trap states neutralized by the ion Li+. These states are due to the crystallographic concentration (101) in the film which gives rise to states at the grain boundaries, then, it is also an important factor in the chemical capacitance. Ti3+ at the grain boundaries, and it is also an important factor in the chemical capacitance Cμ, cb at the TiO2/E interface [31]. Since the pores in the photoelectrode allow the passage of a certain number of cations.

The distribution of the deepest trap states is shown in Fig. 7, it is also an important factor in the variations of the electron diffusion coefficient and chemical capacitance Cμ, ss of the surface states in the solid phase of the mixed-phase electrode.

In Fig. 8, the voltammogram that allows determining the presence and location of the trap states that can neutralize the Na+ ion; in the interval [−0.672, −0.539] Potential vs. Fc+|Fc the presence of shallow trap states has been determined and in the interval [−1.548, −1.118] Potential vs. Fc+|Fc the deep trap states have been determined.

In this case, as can be seen in Fig. 9, the surface trap states neutralized by the cation Na+ are less than those neutralized by the cation Li+ shown in Fig. 6.

For the deep trap states neutralized by the same cation, a lower population is determined than those neutralized by the ion Li+, which reflects a lower diffusion of ions into Na+ in the pores of the photoelectrode, this could be attributed to its size, Fig. 10.

As can be seen in Fig. 11, the ions of K+ fail to diffuse in the film in the same way as ions of Li + or Na+, this is reflected in the null neutralization of shallow trap states. The states that do manage to neutralize are the deep states in the interval [−1.829, −1.304] Potential vs. Fc+|Fc The Fig. 12 shows the distribution of the states neutralized by the ion K+, which, as can be seen, neutralizes a larger number of deep trap states. It should be noted that the cyclic voltammetry (CV) technique was performed with a scanning speed of 25 mV s1, a speed with which the system is polarized, which allowed the distinction and distribution of trap states that manage to neutralize each ion.

4.3 Photoelectrochemical impedance spectroscopy (PEIS)

The trap states occupied during DSSC operation, in the illuminated state, as well as other parameters were estimated by PEIS, which also allows for estimating the steady state illumination diffusion coefficient. Fig. 13 shows the impedance behavior Z(ω) of the different DSSCs with different cations. Here, are determined: ωK+ = 89.047 s−1, ωNa+ = 186.74 s−1 and ωLi+ =353.3568 s−1 additionally, diffusion frequencies were estimated in the TiO2 film anatase of ωK+ = 689.6551 s−1, ωNa+ = 988.0495 s−1 and yωLi+ =2,500 s−1.

Fig. 13

Nyquist plots, for different DSSCs in which the different recombination frequencies were determined.

Using the Software Zview v.3.2b, an adjustment of Bisquert model DX1-Bisquert # 2 was performed [13,25], electrochemical parameters in the illuminated state of the different DSSC were estimated and are listed in Table 1.

Electrochemical parameters on light steady state

Where Rpore = rpore × L is the resistance to the ionic transport in the phase (1), Ret = ret × L−1 is the resistance to the electronic recombination at the TiO2/E interface, Cμ, cb = cμ, cb × L is the chemical capacitance due to the occupied states of the conduction band, Cμ, ss = cμ, ss × L−1 is the chemical capacitance due to the occupied trap states, Rtr = rtr × L is the electron transmission resistance in the phase (2). In addition, it was determined Rt = Rsubs + Rce that is the total resistance of the counter electrode and Rsubs = 196.1 Ω, 146 Ω and 207.6 Ω for the DSSC with Li+, Na+ and K+ respectively, the substrate resistance or the measured FTO.

On the other hand, it was also determined the electron diffusion length LD given by the relation [32]:

(14) LDL=(RetRtr)12

Table 2 shows the values of steady state diffusion coefficients in illumination for the different ions added to the DSSCs:

Length diffusion and steady state coefficients

4.4 Light-dark transients

To determine the diffusion coefficient in the transition from the dark steady state to the illumination steady state, the photocurrent signal emitted by each of the DSSCs was fitted to the model expressed by Eq. (10a,b) for the DSSC with the Li+, Na+ and K+ ions. By taking in each part of the model eC2(0,0) and the boundary conditions, the DSSC with Li+ ions yielded the following values of the fit parameters, shown in Tables 3 and 4:

Values of the fit parameters t ≤ 0

Values of the fit parameters t ≥ 0

For the case of DSSCs with Na+ and K+ ions, she density of injected electrons was considered to be given by C2t0(0,0), and the following parameters were obtained in Tables 5 and 6:

Model parameters for K+ ion for t ≤ 0

Model parameters for K+ ion for t ≥ 0

In Fig. 14, signal of the transients in each DSSC is observed, as can be seen, the profile depends both on the ions used that can modify the model and the different electrochemical parameters.

Fig. 14

Obtained photocurrents profiles of Li+, Na+ and K+ cations during a light step.

Tables 7 and 8 were obtained the model parameters for Na+ ion at the same dark-light, transitory and steady state conditions.

Model parameters for Na+ ion for t ≤ 0

Model parameters for Na+ ion for t ≥ 0

The revealed trap states can be observed in the voltammograms, since Li+ are the smallest ions with access to shallow and deep trap states, therefore, they tend to neutralize those states beneficially promoting a decrease in the available charge at the TiO2/E interface and of the shallow trap states as can be seen in the Cμ,cb and Cμ,ss values for each DSSC. In the case of the DSSC with K+ ion due to null neutralization of the shallow states, presents a higher charge concentration at the TiO2/E interface than the DSSC with Li+ but not higher than the DSSC with Na+, because Na+ ions are smaller and diffuse in greater quantity in the pores of the film.

From the Table 2 to Table 4 for DSSC with Li+ ions, the change in the electron diffusion coefficient during a light step can be determined. It is observed that in the dark state, the electron diffusion coefficient is higher, this is due to the majority of unoccupied states in the gap and conduction band, which facilitates the movement of carriers trapping and releasing from one state to another, which as can be seen from the same tables in the dark stage the concentration is 1.36699 × 107 C cm3. This relatively large value of the diffusion coefficient also explains why a large amount of charge is incorporated during injection. During the dark phase due to a relatively large electron diffusion coefficient, the dark current density takes a value of: 1.9852 × 107 C s1cm2, this is due to the low population of carriers in such a state. On the other hand, it is observed that during the transient, the electron diffusion coefficient is lower, this is due to the charge that has been injected in the trap and conduction states. In addition, the electron diffusion coefficient has decreased since the ions of phase (1) of the mixed electrode perceive the electric charge injected in x = 0.

For DSSCs with Na+ and K+ ions that partially and completely neutralize the shallow states, a differently modeled photocurrent profile is presented, due to a part of the ions are not in the TiO2/E interface, and remain active in phase (1) and being in contact at the bottom of the pore with the FTO promotes the diffusion of holes in the photoelectrode with a diffusion coefficient Dh respectively for each ion. This causes that when electrons are injected in x = 0, an instantaneous displacement of the holes occurs by diffusion in the photoelectrode to the same position causing the photocurrent to drop by recombination, which gives rise to the profile of the photocurrent signal in both cells. Moreover, due to the ion concentration at the TCO/E interface, the overvoltage in that region can be modified, so the extraction rate kext. can be affected since it depends on the voltage [33].

In addition, an effective absorption coefficient was estimated for each DSSC, which, as can be seen in the tables above, has different values for each DSSC. In the tables, the lowest value is for the DSSC with Na+ ion, which could explain why it emits lower photocurrent density as can be seen in Fig 14. This could be due to the interaction of the Na+ ion with the carminic acid molecules of the pigment.

Fig. 15 shows that the size of the ion negatively influences the diffusion coefficient Dn, which could be due to the variation of the velocities va and vd in each ion according to the following situations:

Fig. 15

Relationship between the diffusion coefficient on the photogenerated electrons (Dn), and the ion radii for K+, Na+, and Li+ cations contained in the studied DSSC.

  1. Jsc could decrease if va is large, which could be happening with larger cations adsorbing faster (easier) than smaller ones, leading to a decrease in the rate of electron detrapping from the trap states.

  2. Jsc could decrease if vd is small, which could be due to the fact that once the cations are adsorbed it is very difficult to desorb them, which causes the electrons to stay longer in the trap states and requires more energy to be released, also causing a decrease in the detrapping rate of the electrons from the trap states.

Therefore, if the open circuit voltage Voc is considered to remain constant, a variation in the short circuit current density Jsc or the diffusion coefficient Dn will be reflected in a variation of the DSSC efficiency: η=JscAVocPlight, where A is the cell area. It should be noted that it is necessary to consider the solvation energy of each ion since solvents such as acetonitrile are present and could influence their mobility [34].

5. Conclusions

In this study, it has been possible to model and experimentally determine the evolution of the electron diffusion coefficient from the dark steady state to the light steady state in DSSCs sensitized with carminic acid and a counter electrode based on activated carbon. The different DSSCs used in this work had different cations in the electrolyte solution: Li+, Na+, and K+ which allow revealing the trap states used to quantify the chemical capacitance in the pores of mesoporous TiO2. A strong dependence of the magnitude and shape of the profiles of the photocurrent sent by the DSSCs during the same transient on the cations used has been established since the trap states neutralized by the cations affect the electrochemical parameters in the DSSCs such as the diffusion coefficient itself and the rate of extraction kext. at the contact x = 0.

In addition, it has been possible to estimate an effective value of the absorption coefficient α for each incorporated cation, using a hybrid model for the source of charge carriers in DSSCs excited with a light step, thus having a better understanding of the electrochemical parameters involved in the dark-to-light transition in DSSCs.

Acknowledgments

To the personnel of the Light Scattering Laboratory and the Instrumentation Laboratory of the Center for Applied Physics and Advanced Technology of UNAM campus Juriquilla, Queretaro for facilitating the extraction of natural organic pigments and the performance of electrical measurements on the DSSCs, respectively. Also, to the staff of the Photoelectrochemistry Laboratory at the Center for Research and Technological Development in Electrochemistry in Queretaro, for facilitating the electrochemical tests. Juan Carlos Gallegos Mendoza would like to thank the “Consejo Nacional de Ciencia y Tecnología” for the grant awarded for the development of this project.

References

1. Bisquert J, Garcia-Belmonte G, Pitarch Á, Bolink H. J. Chem. Phys. Lett 2006;422(1–3):184–191.
2. Prasad A. K, Park J. Y, Kang S. H, Ahn K. S. Electrochim. Acta 2022;422:140340.
3. Prasad A. K, Park J.-Y, Jung H. Y, Kang J. W, Kang S.-H, Ahn K.-S. J. Ind. Eng. Chem 2023;117:500–509.
4. Kay A, Grätzel M. Sol. Energ. Mat. Sol. C 1996;44(1):99–117.
5. Ren Y, Zhang D, Suo J, Cao Y, Eickemeyer F. T, Vlachopoulos N, Zakeeruddin S. M, Hagfeldt A, Grätzel M. Nature 2023;613:60–65.
6. Park K.-H, Kim T.-Y, Han S, Ko H.-S, Lee S.-H, Song Y.-M, Kim J.-H, Lee J.-W. Spectrochim. Acta - A: Mol. Biomol. Spectrosc 2014;128:868–873.
7. Rangel D, Gallegos J. C, Vargas S, García F, Rodríguez R. Results Phys 2019;12:2026–2037.
8. Yao Z, Zhang M, Wu H, Yang L, Li R, Wang P. J. Am. Chem. Soc 2015;137(11):3799–3802.
9. Kim B.-G, Chung K, Kim J. Chem. - Eur. J 2013;19(17):5220–5230.
10. He J, Liu Y, Gao J, Han L. Photochem. Photobiol. Sci 2017;16(7):1049–1056.
11. Qian X, Wang X, Shao L, Li H, Yan R, Hou L. J. Power Sources 2016;326:129–136.
12. Listorti A, O’Regan B, Durrant J. R. Chem. Mater 2011;23(15):3381–3399.
13. Bisquert J. Phys. Chem. Chem. Phys 2000;2(18):4185–4192.
14. Fabregat-santiago F, Mora-Sero I, Garcia-Belmonte G, Bisquert J. J. Phys. Chem. B 2003;107(3):758–768.
15. Bisquert J, Garcia-Belmonte G, Fabregat-Santiago F, Bueno P. R. J. Electroanal. Chem 1999;475(2):152–163.
16. Bisquert J. Phys. Chem. Chem. Phys 2003;5(24):5360–5364.
17. Fleig J, Schmid A, Rupp G. M, Slouka C, Navickas E, Andrejs L, Hutter H, Volgger L, Nenning A. Acta Chim. Slov 2016;63(3):509–518.
18. Kopidakis N, Benkstein K. D, van de Lagemaat J, Frank A. J. J. Phys. Chem. B 2003;107(41):11307–11315.
19. Docampo P, Guldin S, Steiner U, Snaith H. J. J. Phys. Chem. Lett 2013;4(5):698–703.
20. Wang Y, Wu D, Fu L.-M, Ai X.-C, Xu D, Zhang J.-P. Phys. Chem. Chem. Phys 2014;16(23):11626–11632.
21. Fabregat-Santiago F, Bisquert J, Garcia-Belmonte G, Boschloo G, Hagfeldt A. Sol. Energy Mater. Sol. Cells 2005;87(1–4):117–131.
22. Zaky A. A, Sergeant P, Stathatos E, Falaras P, Ibrahim M. N. Machines 2022;10(10):882.
23. Cao F, Oskam G, Meyer G. J, Searson P. C. J. Phys. Chem 1996;100(42):17021–17027.
24. de Levie R. Electrochim. Acta 1964;9(9):1231–1245.
25. Bisquert J, Garcia-Belmonte G, Fabregat-Santiago F, Compte A. Electrochem. Commun 1999;1(9):429–435.
26. Sulzer D, Yasuda K. J. Chem. Theory Comput 2018;14(10):5090–5104.
27. Maldon B, Thamwattana N. J. Phys. Commun 2020;4(8):082001.
28. Barnes P. R. F, Miettunen K, Li X, Anderson A. Y, Bessho T, Gratzel M, O’Regan B. C. Adv. Mater 2013;25(13):1881–1922.
29. Bartolotta A, Calogero G. Dye-sensitized solar cells: From synthetic dyes to natural pigments. In : Enrichi F, Righini G. C, eds. Solar Cells and Light Management Elsevier; 2020. p. 107–161.
30. Peiris T. A. N, Sagu J. S, Wijayantha K. G. U, García-Cañadas J. ACS Appl. Mater. Interfaces 2014;6(17):14988–14993.
31. Hengerer R, Kavan L, Krtil P, Grätzel M. J. Electrochem. Soc 2000;147(4):1467.
32. Wang Q, Ito S, Grätzel M, Fabregat-Santiago F, Mora-Seró I, Bisquert J, Bessho T, Imai H. J. Phys. Chem. B 2006;110(50):25210–25221.
33. Dloczik L, Ileperuma O, Lauermann I, Peter L. M, Ponomarev E. A, Redmond G, Shaw N. J, Uhlendorf I. J. Phys. Chem. B 1997;101(49):10281–10289.
34. Duignan T. T, Zhao X. S. Phys. Chem. Chem. Phys 2020;22(43):25126–25135.

Article information Continued

Fig. 1

Schematic representation of the trapping-detrapping process for photogenerated electrons and insertion-releasing process of cations occurring at the dyed TiO2 film-electrolyte junction where x = L stands for the thickness film.

Fig. 2

Two-phase mixed electrode model based on the cylindrical pore model for dyed TiO2 film-electrolyte junctions where the main electrical impedances of a DSSC are shown.

Fig. 3

Micro-Raman spectra obtained for (a) Photoelectrode spectra of TiO2 anatase film 34 μm thick on FTO. (b) Spectra of activated carbon film counter electrode and polyurethane resin as binder, thickness 24 μm.

Fig. 4

Schematic representation of the three-electrode electrochemical cell employed for carrying out the cyclic voltammetry experiments.

Fig. 5

CV response for a TiO2-based working electrode, Pt auxiliary electrode, Ag reference electrode, and 1 mM LiI solution.

Fig. 6

Distribution of surface trap states neutralized by the Li+ ion on the photoelectrode.

Fig. 7

Distribution of deep trap states neutralized by Li+ ion on the photoelectrode.

Fig. 8

CV response for a TiO2-based working electrode TiO2, Pt auxiliary electrode, Ag reference electrode, and 1 mM NaI solution.

Fig. 9

Distribution of shallow trap states neutralized by Na+ ion on the photoelectrode.

Fig. 10

Distribution of deep trap states neutralized by Na+ ion on the photoelectrode.

Fig. 11

CV response for a TiO2-based working electrode TiO2, Pt auxiliary electrode, Ag reference electrode, and 1 mM KI solution.

Fig. 12

Distribution of deep trap states neutralized by the K+ ion on the photoelectrode.

Fig. 13

Nyquist plots, for different DSSCs in which the different recombination frequencies were determined.

Fig. 14

Obtained photocurrents profiles of Li+, Na+ and K+ cations during a light step.

Fig. 15

Relationship between the diffusion coefficient on the photogenerated electrons (Dn), and the ion radii for K+, Na+, and Li+ cations contained in the studied DSSC.

Table 1

Electrochemical parameters on light steady state

Parameter/Ion Li+ Na+ K+
Rpore (Ω) 430 453.02 1601
Ret (Ω) 179.9 357.23 337.9
Cμ,cb (μF) 12.24 507 50.63
Rtr (Ω) 1248 1888.03 2490
Cμ,ss (μF) 15.78 15.5 33.25
Cce (μF) 59.7 5.02 59.7
Rce (Ω) 3130 10352 14408

Table 2

Length diffusion and steady state coefficients

Ion Li+ Na+ K+

Parameter
LD (cm) 1.258 × 10−3 1.363 × 10−3 1.224 × 10−3
Dnlight, ss (cm2 s−1) 5.59 × 10−4 3.46 × 10−4 1.33 × 10−4

Table 3

Values of the fit parameters t ≤ 0

Model for Li+
e-C2(x,t)=e-[c3exp(-kext.Dnx)exp(-kext.2Dnt)+c4t+c5]for t0
Parameter Setting value
c3 × e 1.45576 × 10−7 C cm−3
kext. 1.41 × 103 cm s−1
Dn=Dnst.darkLi+ 38.1 cm2 s−1
c4 × e 1.36996 × 10−7 C s−1cm−3
c5 × e 1.36699 × 10−7 C cm−3

Table 4

Values of the fit parameters t ≥ 0

Model for Li+
Irradiance
I0 = 3.12 × 1017 cm−2 s−1
C2(x,t)=n=1Bnsin(nπxL)exp [-{n2π2DnL2t}]-I0αDne-αx+c1x+c2for t0
Parameter Setting value
ec2 2.7191 × 10−7 C cm−3
ec1 −2.6756 × 103 C cm−4
Dn=Dst.transLi+ 1.865 × 10−5 cm2 s−1
e-C2t0(0,0)=e-(c3+c5)=e-N 2.645 × 10−7 C cm−3
Dh 4.8229 × 10−6 cm2 s−1
B1πL 1.03205 × 10−4
α 0.486 cm−1

Table 5

Model parameters for K+ ion for t ≤ 0

Model for K+
e-C2(x,t)=e-[c3exp(-kext.Dnx)exp(-kext.2Dnt)+c4t+c5]for t0
Parameter Setting value
c3 × e 4.7052 × 10−8 C cm−3
kext. 19.6966 cm s−1
Dn=Dnst.darkKi+ 12.4133 cm2 s−1
c4 × e 2.2393 × 10−7 C s−1cm−3
c5 × e 2.24085 × 10−7 C cm−3

Table 6

Model parameters for K+ ion for t ≥ 0

Model for K+
Irradiance
I0 = 3.12 × 1017 cm−2 s−1
e-C2(x,t)=e-[N4πDhte-(x-L)24Dht-N4πDhte-(x+L)24Dht+c1x+c2-I0αDne-αx]for t0
Parameter Setting value
ec2 4.6707 × 10−6 C cm−3
ec1 −7.465 C cm−4
Dn=Dst.transK+ 6.5984 × 10−3 cm2 s−1
e-C2t0(0,0)=e-(c3+c5)=e-N 2.645 × 10−7 C cm−3
Dh 4.8229 × 10−6 cm2 s−1
α 0.1126 cm−1

Table 7

Model parameters for Na+ ion for t ≤ 0

Model for Na+
Irradiance
I0 = 3.12 × 1017 cm−2 s−1
e-C2(x,t)=e-[c3exp(-kext.Dnx)exp(-kext.2Dnt)+c4t+c5]for t0
Parameter Setting value
c3 × e 1.7216 × 10−8 C cm−3
kext. 19.07 cm s−1
Dn=Dnst.darkNa+ 9.8619 cm2 s−1
c4 × e 5.0514 × 10−8 C s−1cm−3
c5 × e 5.0415 × 10−8 C cm−3

Table 8

Model parameters for Na+ ion for t ≥ 0

Model for Na+
Irradiance I0 = 3.12 × 1017 cm−2 s−1 e-C2(x,t)=e-[N4πDhte-(x-L)24Dht-N4πDhte-(x+L)24Dht+c1x+c2-I0αDne-αx]for t0
Parameter Setting value
ec2 6.3105 × 10−7 C cm−3
ec1 −1.9946 C cm−4
Dn=Dst.transNa+ 2.4942 × 10−5 cm2 s−1
e-C2t0(0,0)=e-(c3+c5)=e-N 6.7631 × 10−8 C cm−3
Dh 2.3088 × 10−6 cm2 s−1
α 8.5 × 10−3 cm−1