# Density Functional Theory Study on the Electronic Structures of Oxadiazole Based Dyes as Photosensitizer for Dye Sensitized Solar Cells

Authors: Umer Mehmood; Ibnelwaleed A. Hussein; Khalil Harrabi; Shakeel Ahmed
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

The molecular structures and UV-visible absorption spectra of complex photosensitizers comprising oxadiazole isomers as theπ-bridges were analyzed by density functional theory (DFT) and time-dependent DFT. The ground state and excited state oxidation potentials, HOMOs and LUMOs energy levels, and electron injection from the dyes to semiconductor TiO2 have been computed in vacuum here. The results show that  all of the dyes may  potentially be good photosensitizers in DSSC. To justify the simulation basis, N3 dye was also simulated under the similar conditions. Simulated absorption spectrum, HOMO, LUMO, and band gap values of N3 were compared with the experimental values. We also computed the electronic structure properties and  absorption spectra  of dye/(TiO2)8 systems  to elucidate the electron injection efficiency at the interface. This work is expected to give proper orientation for experimental synthesis.


---

## Body

## 1. Introduction

The conversion of photoenergy into electrical energy is generally considered as the most potential way to resolve the world energy crisis, owing to its huge reserves. The inorganic silicon based solar cells are being currently used on a commercial scale because of their high efficiency. However, the need of highly purified silicon, use of toxic chemicals in their manufacture, and the high cost have restricted their worldwide use. These constraints encouraged the search for low cost and environmentally friendly solar cells. In this context, dye sensitized solar cells (DSSCs) have received widespread attention in recent years because of their easy processing and low cost [1–6].The working principle of the DSSCs is based on light absorption in a photosensitizer anchored on TiO2 semiconductor. Upon absorption of photoenergy, the electrons in the sensitizer become exited from ground state to the excited state. Due to the difference in energy levels of the electronic states, electrons from the excited state are injected to the conduction band of the semiconductor. As a result the dye becomes oxidized. The electrolyte, which is in contact with the dye, then donates electrons to the dye for the regeneration of sensitizer. Electrolyte then diffuses towards the counter electrode where the reduction reaction takes place and electrolyte restores its initial state by accepting electrons at the counter electrode.The overall efficiency of the DSSCs depends on the photosensitizer. Photosensitizers are categorized into metal complex and metal-free organic sensitizers. But metal-free organic photosensitizers are preferred over ruthenium based sensitizers because of their low cost and good transport properties. The basic structural unit of metal-free dyes is donor-pi-spacer-acceptor (D-π-A). The photovoltaic properties of such dyes can be finely tuned by selecting suitable groups within the D-π-A structure [7]. Many organic groups such as porphyrins [8], styrylarylenes [9], perylenes [10], benzofurans [11], indoles [12], and oxadiazole [13] have been used to tune optical properties of photosensitizers. But the oxadiazole derivatives are considered as the most efficient electron transport materials owing to their good thermal and chemical stabilities and high quantum yield [14]. Many studies have been presented in which 2-(4-biphenylyl)-5-(4-tertbutylphenyl)-(1,3,4) oxadiazole exhibits excellent charge transport properties [13, 15–18]. Zhu et al. successfully synthesized the naphthalimide derivatives containing oxadiazole moiety. They found that oxadiazole moiety increases the electron injection properties and minimizes the carrier recombination [15].DFT/TD-DFT is an effective tool in investigating the ground and excited state properties of photosensitizer complexes as compared to other high level quantum approaches because the computed orbitals are suitable for the typical MO-theoretical analyses and interpretations [19]. Many theoreticians have successfully applied this approach in the designing of photosensitizers [20–25].


## 2. Theoretical Background

The performance of DSSCs is evaluated by incident photon to conversion efficiency (IPCE). It is associated with charge collection efficiency (η
c), electron injection efficiency (ϕ
injc), and light harvesting efficiency (LHE), as (1)
IPCE
=
LHE
×
ϕ
injc
×
η
c
.LHE can be calculated in the following way [26]:(2)
LHE
=
1
-
10
-
f
,where “f” is the absorption of dye associated with maximum absorption also called oscillator strength and ϕ
injc is directly proportional to the free energy of electron injection as [20](3)
ϕ
injc
α
f
-
Δ
G
inject
.Equation (3) shows that the more positive is Δ
G
inject the greater will be the electron injection efficiency, while Δ
G
inject is the difference between oxidation potential energy of the excited state (Eox
dye
∗) and the reduction potential energy of TiO2 conduction band (E
C
B), which can be described as [24](4)
Δ
G
inject
=
Eox
dye
∗
-
E
C
B
.TiO2 cluster was also simulated and its E
C
B was found to be −4.01 eV (in vacuum). Similarly, Eox
dye
∗ can be calculated by using the following equation [21]:(5)
Eox
dye
∗
=
Eox
dye
-
Δ
E
,where Eoxdye (-HOMO) [22] is the dye’s ground state oxidation potential and Δ
E is the lowest absorption energy associated with λ
max
⁡.


## 3. Computational Details

All the DFT/TD-DFT calculations were executed using Amsterdam Density Functional (ADF) program (2013.01) [19]. The ground state geometries of oxadiazole dyes, dye/(TiO2)8, and standard dyes were optimized by applying hybrid B3LYP level together with triple-ζ polarization basis function. (TiO2)8 nanoparticles cluster was also simulated by considering generalized gradient approximation (GAD) at BYLP level and triple-ζ polarization basis function. UV-Vis spectra of oxadiazole dyes were simulated in ethanol solvent. Here the conductor-like screening model (COSMO) was selected to take the solvent effects into account [20], while the excitation energies were examined using TD-DFT and statistical average of orbital potentials (SAOP) model including the salvation effects [27]. Seventy singlet-singlet transitions were calculated, which are sufficient to fully define the whole absorption spectrum. In all the calculations, the relativistic effects were taken into account by the zero-order regular approximation (ZORA) Hamiltonian in its scalar approximation [28].


## 4. Results and Discussion

### 4.1. Designed Systems

The structures and names of new class of dyes are shown in Figure1. In these structures, biphenyl unit was used as electron-donating moiety and carboxyl and cyano groups (-COOH and -CN) were introduced as the electron acceptor and the anchor groups because of their high ability of electron-withdrawing and bonding to semiconductor. While oxadiazole isomers were introduced as π-conjugation to bridge the donor-acceptor systems, a double bond and a thiophene unit were also introduced to the pi-conjugation system for the fine tuning of molecular planar configurations and to broaden the absorption spectra. Here the performance of dyes has been tested by introducing different isomers of oxadiazoles as π-conjugated bridges.Figure 1
Chemical structures of organic dyes.


### 4.2. Energy Levels

The HOMOs, LUMOs, and band gap energies of photosensitizers play an important role in providing the thermodynamic driving force for the electrons injection. For efficient charge transfer, the LUMOs of dyes must be more negative (-LUMOs of dyes must be less negative) than the conduction band of the semiconductor while HOMO levels must be more positive (-HOMO levels must be more negative) than the redox potential of electrolyte. The HOMOs, LUMOs, and LUMO-HOMO energy gap of N3 computed at hybrid B3LYP level are in good agreement with the experimental values as shown in Table1(a). Thus, hybrid B3LYP function with triple-ζ polarization basis function can be a valid function for the geometry optimization of oxadiazole based photosensitizers.Table 1
(a) The FMOs, H-Lgap energies, and absorption peaks of N3. (b) The FMO (eV) and H-Lgap (eV) energies of systems 1–4.
(a)
Dye
LUMO (eV)
HOMO (eV)
H-Lgap (eV)
Wavelength (nm) of absorption peaks
N
3
+ (experimental)
−3.47
−5.36
1.89
311, 385, 530
N3 (simulated)
−3.67
−5.63
1.96
322, 405, 542
+Experimental [29].
(b)
System
LUMOs (eV)
HOMOs (eV)
H-Lgap (eV)
1
−3.867
−5.682
1.815
2
−3.630
−5.508
1.878
3
−3.762
−5.865
2.103
4
−3.882
−5.919
2.037Table1(b) shows the HOMOs, LUMOs, and HOMO-LUMO energy gaps of oxadiazole photosensitizers. The HOMO levels of the systems are in the order of system 4 (−5.919) < system 3 (−5.865) system 1 (−5.682) < system 2 (−5.508). The LUMO energy levels are in the order of system 4 (−3.882) < system 1 (−3.867) < system 3 (−3.762) < system 2 (−3.630). It can be seen that the insertion of oxadiazole isomers as π-conjugation significantly influences the HOMO and LUMO energy levels. Similarly, the H-L
g
a
p of systems are in the order of system 1 (1.815) < system 2 (1.878) < system 4 (2.037) < system 3 (2.103). These results suggest that systems 1–4 should be capable of injecting electrons into the conduction bands of TiO2.


### 4.3. Absorption Spectra of Dyes

An efficient photosensitizer should show intense absorption in the visible region (400 nm to 700 nm). The absorption peaks of simulated N3 sensitizer as shown in Figure2(a) are closer to the experimental values shown in Table 1(a). Thus COSMO solvation model along with SAOP function can be a valid choice for computing the excited state properties of oxadiazole based photosensitizers.Figure 2
Simulated absorption spectra of (a) N3 and (b) systems 1–4.
(a)
(b)Similarly, the simulated absorption spectra of systems 1–4 are shown in Figure2(b). All the absorption spectra of dyes can be clearly divided into two regions with first intense peak in the region of 326–500 nm and second peak in the region of 472–700 nm. The band in the UV-region is probably emerging from the transitions of electrons localized within the biphenyl unit. The absorption occurring in the visible region is due to transfer of electrons from donor to acceptor, which is sensitive to the nature of the conjugation pathway and red shifts on progressive addition of thiophene units. System 2 with 1,2,3-oxadiazole as pi-bridge exhibits a significant red shift as compared to other systems. It can also be observed that system 2 with 1,2,3-oxadiazole bridge shows a broader peak.Ground and excited state potential (E
o
x
d
y
e), maximum absorption (λ
max
⁡), oscillation strength (f), and main transitions are presented in Table 2. The main molecular orbitals (MOs) involved in the dominant electron transitions of systems 1–4 are shown in Figure 3. The major transitions may be attributed to the transfer of charge from phenyl group to -COOH and -CN through the thiophene unit and π-conjugation.Table 2
Optical, redox potentials, oscillator strengths, and major transitions of oxadiazole based photosensitizers.
System
E
ox
dye  (eV)
E
ox
dye
∗   (nm)
Δ
E (nm)
λ
max
⁡  (nm)
f
Main transition
1
5.682
3.482
2.3
540
0.59
H → L (96%)
2
5.508
3.308
2.2
565
0.58
H → L (91%)
3
5.865
3.765
2.2
561
0.72
H-1 → L (86%)
4
5.919
3.719
2.2
563
0.71
H → L (95%)Figure 3
Major transition of systems 1–4.


### 4.4. UV-Vis Absorption Spectra of Dyes on TiO2 Surface

The interface between sensitizers and TiO2plays an important role in the electron injection efficiency. But (TiO2)8 model was selected because of its balanced electronic properties (conduction band of ~−4.0 eV and band gap of 3.18 eV) [27, 28]. BAND mode was used to simulate the anatase TiO2 cluster. Here, we selected tetragonal anatase crystal structure with single-layer (001) surface slab. Then, a 4 × 1 supercell was created from this slab. All atoms were mapped within the unit cell [5, 27, 28]. The models of dyes linked on (TiO2)8 are shown in Figure 4, while the simulated UV-Vis spectra of dyes/(TiO2)8 are shown in Figure 5. The spectra show a dramatic red shift as compared to those in solution. The red shifts values of system 1, system 2, system 3, and system 4 are 20, 25, 15, and 19, respectively. This may be due to the increased delocalization of the π
∗ orbital of the conjugated framework caused by the interaction between the carboxylate group and the Ti+4 ions that directly decreases the energy of the π
∗ level. It can also be observed that dye (3)/TiO2 show higher red shift value as compared to other systems. The calculated excitation energies, oscillation strength, λ
max
⁡ (nm), and dominant transitions are shown in Table 3.Table 3
Lowest absorption energyΔ
E, maximum absorption λ
max
⁡, and oscillating strength (f).
System
λ
max
⁡ (nm)
Δ
E (eV)
f
1
560
2.217
0.31
2
690
2.104
0.34
3
576
2.140
0.36
4
582
2.159
0.30Figure 4
The model of systems 1–4 with (TiO2)8.Figure 5
Simulated absorption spectra of (TiO2)8/systems 1–4.


### 4.5. Free Energy of Electron Injection and Light Harvesting Efficiency

LHE andΔ
G
inject can be calculated by using (2) and (4), respectively. The results are given in Table 4. It can be found from Table 4 that all the calculated Δ
G
inject values are positive, which indicates that the conduction band edge of TiO2 lies blow the excited state of the dyes, thus favoring electron injection. The Δ
G
inject of systems 1–4 increases in the order of system 2 (0.692) > system 1 (0.618) > system 3 (0.335) > system 4 (0.281). Therefore, the order of driving force of the systems will be system 2 > system 1 > system 3 > system 4. There is only slight difference of LHE for systems 1–4. It demonstrates that all the photosensitizers will give comparable photocurrents. These results indicate that system 2 can be considered as plausible due to more positive Δ
G
inject (0.632 eV) and larger LHE value (0.74), which results in higher IPCE.Table 4
Free energy of electron injectionΔ
G
inject (eV) and light harvesting efficiency (LHE).
System
Δ
G
inject
LHE
1
0.613
0.738
2
0.692
0.740
3
0.335
0.801
4
0.281
0.796


## 4.1. Designed Systems

The structures and names of new class of dyes are shown in Figure1. In these structures, biphenyl unit was used as electron-donating moiety and carboxyl and cyano groups (-COOH and -CN) were introduced as the electron acceptor and the anchor groups because of their high ability of electron-withdrawing and bonding to semiconductor. While oxadiazole isomers were introduced as π-conjugation to bridge the donor-acceptor systems, a double bond and a thiophene unit were also introduced to the pi-conjugation system for the fine tuning of molecular planar configurations and to broaden the absorption spectra. Here the performance of dyes has been tested by introducing different isomers of oxadiazoles as π-conjugated bridges.Figure 1
Chemical structures of organic dyes.


## 4.2. Energy Levels

The HOMOs, LUMOs, and band gap energies of photosensitizers play an important role in providing the thermodynamic driving force for the electrons injection. For efficient charge transfer, the LUMOs of dyes must be more negative (-LUMOs of dyes must be less negative) than the conduction band of the semiconductor while HOMO levels must be more positive (-HOMO levels must be more negative) than the redox potential of electrolyte. The HOMOs, LUMOs, and LUMO-HOMO energy gap of N3 computed at hybrid B3LYP level are in good agreement with the experimental values as shown in Table1(a). Thus, hybrid B3LYP function with triple-ζ polarization basis function can be a valid function for the geometry optimization of oxadiazole based photosensitizers.Table 1
(a) The FMOs, H-Lgap energies, and absorption peaks of N3. (b) The FMO (eV) and H-Lgap (eV) energies of systems 1–4.
(a)
Dye
LUMO (eV)
HOMO (eV)
H-Lgap (eV)
Wavelength (nm) of absorption peaks
N
3
+ (experimental)
−3.47
−5.36
1.89
311, 385, 530
N3 (simulated)
−3.67
−5.63
1.96
322, 405, 542
+Experimental [29].
(b)
System
LUMOs (eV)
HOMOs (eV)
H-Lgap (eV)
1
−3.867
−5.682
1.815
2
−3.630
−5.508
1.878
3
−3.762
−5.865
2.103
4
−3.882
−5.919
2.037Table1(b) shows the HOMOs, LUMOs, and HOMO-LUMO energy gaps of oxadiazole photosensitizers. The HOMO levels of the systems are in the order of system 4 (−5.919) < system 3 (−5.865) system 1 (−5.682) < system 2 (−5.508). The LUMO energy levels are in the order of system 4 (−3.882) < system 1 (−3.867) < system 3 (−3.762) < system 2 (−3.630). It can be seen that the insertion of oxadiazole isomers as π-conjugation significantly influences the HOMO and LUMO energy levels. Similarly, the H-L
g
a
p of systems are in the order of system 1 (1.815) < system 2 (1.878) < system 4 (2.037) < system 3 (2.103). These results suggest that systems 1–4 should be capable of injecting electrons into the conduction bands of TiO2.


## 4.3. Absorption Spectra of Dyes

An efficient photosensitizer should show intense absorption in the visible region (400 nm to 700 nm). The absorption peaks of simulated N3 sensitizer as shown in Figure2(a) are closer to the experimental values shown in Table 1(a). Thus COSMO solvation model along with SAOP function can be a valid choice for computing the excited state properties of oxadiazole based photosensitizers.Figure 2
Simulated absorption spectra of (a) N3 and (b) systems 1–4.
(a)
(b)Similarly, the simulated absorption spectra of systems 1–4 are shown in Figure2(b). All the absorption spectra of dyes can be clearly divided into two regions with first intense peak in the region of 326–500 nm and second peak in the region of 472–700 nm. The band in the UV-region is probably emerging from the transitions of electrons localized within the biphenyl unit. The absorption occurring in the visible region is due to transfer of electrons from donor to acceptor, which is sensitive to the nature of the conjugation pathway and red shifts on progressive addition of thiophene units. System 2 with 1,2,3-oxadiazole as pi-bridge exhibits a significant red shift as compared to other systems. It can also be observed that system 2 with 1,2,3-oxadiazole bridge shows a broader peak.Ground and excited state potential (E
o
x
d
y
e), maximum absorption (λ
max
⁡), oscillation strength (f), and main transitions are presented in Table 2. The main molecular orbitals (MOs) involved in the dominant electron transitions of systems 1–4 are shown in Figure 3. The major transitions may be attributed to the transfer of charge from phenyl group to -COOH and -CN through the thiophene unit and π-conjugation.Table 2
Optical, redox potentials, oscillator strengths, and major transitions of oxadiazole based photosensitizers.
System
E
ox
dye  (eV)
E
ox
dye
∗   (nm)
Δ
E (nm)
λ
max
⁡  (nm)
f
Main transition
1
5.682
3.482
2.3
540
0.59
H → L (96%)
2
5.508
3.308
2.2
565
0.58
H → L (91%)
3
5.865
3.765
2.2
561
0.72
H-1 → L (86%)
4
5.919
3.719
2.2
563
0.71
H → L (95%)Figure 3
Major transition of systems 1–4.


## 4.4. UV-Vis Absorption Spectra of Dyes on TiO2 Surface

The interface between sensitizers and TiO2plays an important role in the electron injection efficiency. But (TiO2)8 model was selected because of its balanced electronic properties (conduction band of ~−4.0 eV and band gap of 3.18 eV) [27, 28]. BAND mode was used to simulate the anatase TiO2 cluster. Here, we selected tetragonal anatase crystal structure with single-layer (001) surface slab. Then, a 4 × 1 supercell was created from this slab. All atoms were mapped within the unit cell [5, 27, 28]. The models of dyes linked on (TiO2)8 are shown in Figure 4, while the simulated UV-Vis spectra of dyes/(TiO2)8 are shown in Figure 5. The spectra show a dramatic red shift as compared to those in solution. The red shifts values of system 1, system 2, system 3, and system 4 are 20, 25, 15, and 19, respectively. This may be due to the increased delocalization of the π
∗ orbital of the conjugated framework caused by the interaction between the carboxylate group and the Ti+4 ions that directly decreases the energy of the π
∗ level. It can also be observed that dye (3)/TiO2 show higher red shift value as compared to other systems. The calculated excitation energies, oscillation strength, λ
max
⁡ (nm), and dominant transitions are shown in Table 3.Table 3
Lowest absorption energyΔ
E, maximum absorption λ
max
⁡, and oscillating strength (f).
System
λ
max
⁡ (nm)
Δ
E (eV)
f
1
560
2.217
0.31
2
690
2.104
0.34
3
576
2.140
0.36
4
582
2.159
0.30Figure 4
The model of systems 1–4 with (TiO2)8.Figure 5
Simulated absorption spectra of (TiO2)8/systems 1–4.


## 4.5. Free Energy of Electron Injection and Light Harvesting Efficiency

LHE andΔ
G
inject can be calculated by using (2) and (4), respectively. The results are given in Table 4. It can be found from Table 4 that all the calculated Δ
G
inject values are positive, which indicates that the conduction band edge of TiO2 lies blow the excited state of the dyes, thus favoring electron injection. The Δ
G
inject of systems 1–4 increases in the order of system 2 (0.692) > system 1 (0.618) > system 3 (0.335) > system 4 (0.281). Therefore, the order of driving force of the systems will be system 2 > system 1 > system 3 > system 4. There is only slight difference of LHE for systems 1–4. It demonstrates that all the photosensitizers will give comparable photocurrents. These results indicate that system 2 can be considered as plausible due to more positive Δ
G
inject (0.632 eV) and larger LHE value (0.74), which results in higher IPCE.Table 4
Free energy of electron injectionΔ
G
inject (eV) and light harvesting efficiency (LHE).
System
Δ
G
inject
LHE
1
0.613
0.738
2
0.692
0.740
3
0.335
0.801
4
0.281
0.796


## 5. Conclusion

Four photosensitizers based on oxadiazole isomers as pi-conjugated spacer were designed and simulated using DFT/TD-DFT. The UV-Vis spectroscopic properties as well as the driving force (Δ
G
inject) of systems 1–4 show that these dyes are potential to be a good photosensitizer in DSSC. The results show that LUMO of the dyes is greater than the conduction band of TiO2 indicating that these dyes are thermodynamically favorable in charge transfer into the conduction band of TiO2. Similarly, from system 1 to system 4, the absorption spectra are red-shifted due to the increased delocalization of the π
∗ orbital of the conjugated framework caused by the interaction between the carboxylate group and the Ti+4 ions that directly decreases the energy of the π
∗ level. Moreover, there is no distinct difference in the LHE of systems 1–4. Therefore, it can be concluded that system 2 with 1,2,3-oxadiazole shows a balance among the different crucial parameters and is expected to be a promising sensitizer in the DSSC field.



---

*Source: 286730-2015-06-28.xml*