Titanium dioxide (TiO2) is one of the most used oxides in renewable energy applications, such as hydrogen production, photovoltaics, and light-emitting diodes. To further improve the efficiency of the devices, doping strategies are used to modify their fundamental properties. Here, we used density functional theory (DFT) simulations to explore the effect of all the halogen dopants on the structural, electronic, and optical properties of TiO2. We investigated both the interstitial and the oxygen substitutional positions, and for the optimized structures, we used hybrid DFT calculations to predict the electronic and optical properties. In all cases, we found that halogen dopants reduce the bandgap of the pristine TiO2 while gap states also arise. The halogen dopants constitute a single acceptor when they occupy interstitial sites, while when they are inserted in oxygen sites, they act as donors. This can be established by the states that form above the valence band. It is proposed that these states contribute to the significant changes in the optical and electronic properties of TiO2 and can be beneficial to the photovoltaic and photocatalytic applications of TiO2. Importantly, the iodine doping of TiO2 significantly reduces the bandgap of TiO2 while increasing its dielectric constant, making it suitable for light-harvesting applications.
INTRODUCTION
TiO2 is a semiconductor with n-type conductivity that crystallizes in anatase, rutile, and brookite structures.1–6 Doping is considered an appropriate and low-cost method to increase the electrical conductivity and the transparency of TiO2.7–9 It is mainly used in solar cell technology as an electron transport layer material,10 in photocatalysis for hydrogen production,11 and in sensing applications.12 Many experiments have investigated the halogen doping of TiO2,13–16 with fluorine being the most used dopant for optoelectrical applications.17,18 As it is stated in these studies, fluorine doping modifies the properties of the bulk and the surface.17,18 In particular, fluorine leads to highly crystalline and photoactive samples as it hampers the anatase to rutile phase transformation.13 Fluorine also enables the formation of defects on the surface of anatase, and this leads to an improved absorption in the near-ultra-violet (UV) region and favors the formation of long-living luminescent surface trapping sites.13 Wang et al.19 investigated the photocatalytic activity of Cl:TiO2 with various doping concentrations. Their results showed that the photocatalytic activity increases up to a level with respect to the Cl content. This is due to the additional energy levels formed in the bandgap when the surface Ti3+ exceeds a particular value. As a result, higher Cl dopant results in excessive surface Ti3+, which act as surface recombination centers for e−/h+ and accordingly suppress the photocatalytic activity.19,20
Regarding bromine doping, there are fewer reports that focus on the effect on TiO2, but it was predicted that it reduces the bandgap of TiO2 and the co-doping with chlorine enhances the photocatalytic activity.21 The electronic and optical properties of I:TiO2 were extensively investigated from an experimental point of view. In many works, it was determined that a profound gap reduction occurred due to iodine incorporation.16,22,23 For example, Hong et al.24 predicted that iodine doping of TiO2 that is prepared with a hydrolysis method exhibited an efficient decomposition of phenol under constant visible irradiation.
Apart from the experimental results, density functional theory (DFT) studies also examine the effect of halogen doping in TiO2. In their work, Yang et al.25 investigated the role of halogen dopants in substitutional positions in TiO2. They predicted that iodine reduces the bandgap from 3.16 to 2.50 eV. Notably, in their simulations, they used the “scissors operation” to tackle the well-known problem of gap underestimation.26 Apart from the scissors operation, a variety of methods have been developed to correct the bandgap of TiO2, such as the Hubbard + U method,27 and the more computationally expensive hybrid pseudopotentials that produce the most accurate bandgap values.
Here, we employed hybrid DFT calculations to extensively investigate the effect of halogens in the bulk properties of TiO2. Specifically, we examine the changes in the electronic and optical properties and the potential applications in photocatalysis and photovoltaics. In this perspective, we predicted a bandgap that reaches a value of 3.56 eV, agreeing with the experimental bandgap.28 Herein, we have examined interstitial and substitutional defect sites, and for the minimum energy systems, we predicted the changes in the electronic, electrical, and optical properties. Using the total and partial density of states (pDOS), we analyzed the minimum sites and their effect on the energy gap. Furthermore, we have predicted the absorption for all the doping cases and we discuss our results with the available literature data.
METHODOLOGY
For our study, we used the Cambridge Serial Total Energy Package (CASTEP).29,30 We employed the hybrid functional PBE0 to encounter the bandgap underestimation due to the effect of the localized electrons.31 We chose the cut-off energy at 800 eV for a supercell consisting of 48 atoms. In order to sample the first Brillouin zone during the geometry relaxation, we use a k-point mesh of 2 × 2 × 3, while for the DOS calculations, we choose a set of 4 × 4 × 4. Finally, the convergence criteria of our simulations were set at 2.0 · 10−5 eV/atom for the SCF tolerance, 0.05 eV/Å for the force tolerance, and 0.001 Å for the max displacement tolerance. The analysis of the electronic bands and the semi-classical Boltzmann theory, to calculate the electrical conductivity, were employed using the BoltzTrap software.32
RESULTS
Structural properties
Anatase TiO2 belongs to the I4/amd space group of the tetragonal system. The calculated lattice parameters of the unit cell are a = b = 3.74 Å and c = 9.45 Å, which are in good agreement with the experimental predicted values.33 In Figs. 1(a)–1(h), the optimized dopant position in TiO2 are presented. The introduction of dopants in the lattice creates local distortions, which translate to changes in the lattice parameters. The structural parameters are shown in Table I. In Fig. 2, we have plotted the change of the volume for each doping case, as a function of the radius of the dopant. The volume of the material increases with the increase of the dopant radius with I:TiO2 having the highest volume.
The optimized structure of (a) FO:TiO2, (b) Fi:TiO2, (c) ClO:TiO2, (d) Cli:TiO2, (e) BrO:TiO2, (f) Bri:TiO2, (g) IO:TiO2, and (h) Ii:TiO2.
The optimized structure of (a) FO:TiO2, (b) Fi:TiO2, (c) ClO:TiO2, (d) Cli:TiO2, (e) BrO:TiO2, (f) Bri:TiO2, (g) IO:TiO2, and (h) Ii:TiO2.
The lattice parameters after geometry relaxation for every examined dopant.
| . | a (Å) . | c (Å) . | Volume (Å3) . | |||
|---|---|---|---|---|---|---|
| Interstitial . | Substitutional . | Interstitial . | Substitutional . | Interstitial . | Substitutional . | |
| F:TiO2 | 3.72 | 3.76 | 9.56 | 9.40 | 132.29 | 132.54 |
| Cl:TiO2 | 3.75 | 3.84 | 9.56 | 9.38 | 134.07 | 138.31 |
| Br:TiO2 | 3.76 | 3.86 | 9.57 | 9.41 | 135.30 | 140.20 |
| I:TiO2 | 3.79 | 3.77 | 9.59 | 9.87 | 137.75 | 140.28 |
| Undoped | 3.74 | 9.45 | 132.18 | |||
| . | a (Å) . | c (Å) . | Volume (Å3) . | |||
|---|---|---|---|---|---|---|
| Interstitial . | Substitutional . | Interstitial . | Substitutional . | Interstitial . | Substitutional . | |
| F:TiO2 | 3.72 | 3.76 | 9.56 | 9.40 | 132.29 | 132.54 |
| Cl:TiO2 | 3.75 | 3.84 | 9.56 | 9.38 | 134.07 | 138.31 |
| Br:TiO2 | 3.76 | 3.86 | 9.57 | 9.41 | 135.30 | 140.20 |
| I:TiO2 | 3.79 | 3.77 | 9.59 | 9.87 | 137.75 | 140.28 |
| Undoped | 3.74 | 9.45 | 132.18 | |||
Electrical properties
In Fig. 3, the DOS results for the halogen doped and undoped TiO2 are presented. Generally, there are band structure changes in all the doping cases in interstitial and substitutional defects. In Fig. 3(i), we include as a reference the DOS of the undoped TiO2 with a predicted gap value of 3.56 eV. This value agrees with previous experiments28 and theoretical studies.34 In particular, we predict a reduction to the bandgap, and specifically in the interstitial defects, energy levels near the conduction and the valence band are created. In Fig. 3(a), the DOS of Fo:TiO2 is shown. It is seen that the F substitutional does not exhibit any energy states inside the gap, but it slightly decreases the gap to 3.43 eV. Conversely, interstitial doping [Fig. 3(b)] gives rise to gap states near the conduction and the valence band edges. In this case, the bandgap is significantly reduced to a value of 3.15 eV. For the Cl doping, we can see that in the substitution of oxygen the bandgap is again reduced to a value of 3.33 eV, while in the interstitial position, it reaches the value of 3.09 eV. Similar to the fluorine case, we can predict the formation of energy states at the edge of the conduction band. In previous work on F and Cl doping, we also observed a significant reduction in the case of chlorine doping even at lower dopant contents.27,35
The density of states graphs for (a) FO:TiO2, (b) Fi:TiO2, (c) ClO:TiO2, (d) Cli:TiO2, (e) BrO:TiO2, (f) Bri:TiO2, (g) IO:TiO2, (h) Ii:TiO2, and (i) the undoped.
The density of states graphs for (a) FO:TiO2, (b) Fi:TiO2, (c) ClO:TiO2, (d) Cli:TiO2, (e) BrO:TiO2, (f) Bri:TiO2, (g) IO:TiO2, (h) Ii:TiO2, and (i) the undoped.
Regarding bromine doping, we predict that the BrO:TiO2 has a bandgap value of 3.24 eV, while interstitial doping significantly reduces the gap at a value of 2.70 eV. Again, gap states are formed near the valence band, which are mainly attributed to bromine as well as near the conduction band with the contribution of oxygen and bromine. The case is quite different with iodine dopant. Interestingly, when inserted as a substitute for oxygen, it increases the bandgap to a value of 3.65 eV. This time, the Io:TiO2 exhibits additional energy levels near the valence band maximum (VBM). Such states near the oxide’s VBM and those near the conduction band can act as shallow and deep acceptors, respectively, hence contributing to the performance of solar-cell devices based on TiO2. Conversely, the interstitial doping reduces the bandgap to a value of 2.66 eV and mid-gap states are formed. These states can also act as traps for the photogenerated carriers leading to a decrease in the device photocurrent and efficiency. The present results considering the F, Cl, Br, and I substitutes are consistent with the work of Yang et al.25 regarding the decrease of the forbidden bandgap. Nevertheless, using the hybrid functionals, we obtain different results for I substitutional doping, where we see both the bandgap to increase and the formation of donor-type states.
Generally, we conclude that the substitutional doping of halogen atoms does not introduce any available energy levels, with iodine being an exception. This contrasts with the interstitial dopants that change the edges both in the valence and conduction band. This strongly suggests the presence of free electrons that enhance the n-type conductivity of TiO2 and has also been observed in other similar oxides, such as SnO2.36
In Table II, we have summarized the bandgap values for all the doping cases. The lowest bandgap is achieved with Ii:TiO2 and Bri:TiO2. In the case of Io:TiO2 at a doping percentage of 2.08%, we can see the highest gap value. Therefore, it is proposed that iodine doping at these doping percentages should be examined as an electron transport layer in photovoltaics, while bromine should also be examined as a dopant to develop better photocatalysts.
The bandgap values for every halogen doped TiO2.
| . | Bandgap (eV) . | |
|---|---|---|
| Interstitial . | Substitutional . | |
| F:TiO2 | 3.15 | 3.43 |
| Cl:TiO2 | 3.09 | 3.33 |
| Br:TiO2 | 2.70 | 3.24 |
| I:TiO2 | 2.66 | 3.65 |
| Undoped | 3.56 | |
| . | Bandgap (eV) . | |
|---|---|---|
| Interstitial . | Substitutional . | |
| F:TiO2 | 3.15 | 3.43 |
| Cl:TiO2 | 3.09 | 3.33 |
| Br:TiO2 | 2.70 | 3.24 |
| I:TiO2 | 2.66 | 3.65 |
| Undoped | 3.56 | |
In Fig. 4, the projected DOS for the interstitial doped cases are examined. Through the pDOS, we can decode the hybridization of the orbitals that create the gap states. As we show in Fig. 4(a), F doping has new additional energy levels near the conduction band edge. The states located at 2.5 eV arise due to the hybridization of O-2p with F-2p orbitals. In the case of chlorine, it is seen that the Cl-3p orbitals contribute to the gap states near the valence and the conduction band. We also observe the same in the bromine interstitial where the states occur because of the hybridization of Br-4p and O-2p orbitals. Finally, in the case of iodine interstitial, the inner gap states arise at 0.2 eV due to the hybridization of I-5p, Ti-3d, and O-2p. The present results agree with previous DFT work on halogen doped SnO2, which in turn indicates that the formation of gap states in halogen doped oxides may have a universal origin.37
The projected density of states graphs for (a) Fi:TiO2, (b) Cli:TiO2, (c) Bri:TiO2, and (d) Ii:TiO2.
The projected density of states graphs for (a) Fi:TiO2, (b) Cli:TiO2, (c) Bri:TiO2, and (d) Ii:TiO2.
To further investigate the electrical properties of halogen doped anatase TiO2, such as the electrical conductivity, the BoltzTrap code,32 which is based on the semiclassical Boltzmann transport equation, was used for each dopant in both interstitial and substitutional positions. To gain a better insight into the electrical conductivity changes, we predict the conductivity value at the Fermi level at 300 K. The electrical conductivity estimates the flow of free charge carriers (e−/h+ pairs) inside the bulk structure as a function of temperature and chemical potential. The calculated electrical conductivity to scattering time (σ/τ) of halogen doped TiO2 is presented in Fig. 5, while the results for the undoped agree with the available literature.38–40 We predict a reduction of the conductivity with the increase of the temperature for every case. We can see that both in the substitutional and interstitial doping, bromine and iodine have shown interesting results as they have improved the conductivity in the n-type area of the chemical potential. In Table III, we have collected all the conductivity values at the conduction band edge, where significant changes are predicted for each compound. This improvement of conductivity is beneficial for applying halogen doping on modern solar cells. The increase in the case of iodine interstitial is of great importance, and it opens the road to the application of I:TiO2 to photovoltaics as an electron transport layer or to sensors as a primary sensing material. Furthermore, improvements may be achieved by co-doping strategies and processing conditions as has been previously demonstrated in related materials.41,42
The conductivity graphs for (a) substitutional and (b) interstitial doping. The conductivity of the undoped is presented with purple for reference.
The conductivity graphs for (a) substitutional and (b) interstitial doping. The conductivity of the undoped is presented with purple for reference.
The conductivity to scattering time values for every halogen doped TiO2 in the conduction band edge.
| . | Conductivity to scattering time (S/ms) . | |
|---|---|---|
| Interstitial . | Substitutional . | |
| F:TiO2 | 0.20 × 105 | 0.74 × 105 |
| Cl:TiO2 | 0.66 × 105 | 0.78 × 105 |
| Br:TiO2 | 0.97 × 105 | 0.91 × 105 |
| I:TiO2 | 0.10 × 106 | 0.47 × 105 |
| Undoped | 0.79 × 105 | |
| . | Conductivity to scattering time (S/ms) . | |
|---|---|---|
| Interstitial . | Substitutional . | |
| F:TiO2 | 0.20 × 105 | 0.74 × 105 |
| Cl:TiO2 | 0.66 × 105 | 0.78 × 105 |
| Br:TiO2 | 0.97 × 105 | 0.91 × 105 |
| I:TiO2 | 0.10 × 106 | 0.47 × 105 |
| Undoped | 0.79 × 105 | |
Optical properties
Finally, in Fig. 6, the absorption coefficiency is presented for the undoped and doped samples. All the cases considered here have an absorption range in the ultraviolet light region. Compared to the undoped cases, we predicted that, in the interstitial doping [Fig. 6(a)], the absorbance is significantly enhanced with the Ii:TiO2 having the highest absorbance among the other dopants. In the case of substitutional doping [Fig. 6(b)], I and Br doping again show an improvement to the absorbance compared to the other halogens or the undoped TiO2. The present results for the absorbance of undoped TiO2 agree well with previous DFT43 and experimental studies.44 From the optical results, it can be concluded that when halogens incorporate oxygen sites, it is more beneficial for applying TiO2 to photovoltaic technologies. This can be achieved when the dopant is inserted in a small amount where it sitsite in an oxygen vacancy site. However, when the dopant is experimentally inserted in a large amount then due to the increase of the absorption it is more practical to be examined for alternative energy production, such as hydrogen production. Previous studies have shown that a high percentage of doping significantly improves the photocatalytic properties of TiO2 and thus the hydrogen production.28
The absorption coefficiency for (a) the undoped and (b) doped samples.
CONCLUSION
In this work, we used DFT to predict the changes of the structural, electrical, and optical properties of TiO2 due to halogen doping. In all the examined cases, we predicted that the halogen incorporation significantly reduces the bandgap, except for the case of iodine, where when inserted as an oxygen substitutional it increases the bandgap. We also predicted the electrical conductivity of interstitial and substitutional doping and we found that iodine and bromine significantly enhance it. Finally, we predicted that halogens improve the absorbance of TiO2. Interestingly, iodine doping has the highest absorbance among the other halogens, making it the most suitable dopant candidate for light-harvesting applications.
ACKNOWLEDGMENTS
P.-P.F., M.V., D.D., and A.C. are grateful for LRF ICON funding from the Lloyd’s Register Foundation, a charitable foundation helping to protect life and property by supporting engineering-related education, public engagement and the application of research. N.K. and A.C. acknowledges support from European Union’s H2020 Program under Grant Agreement No. 824072- HARVESTORE.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Petros-Panagis Filippatos: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal). Nikolaos Kelaidis: Investigation (equal); Methodology (equal); Visualization (equal). Maria Vasilopoulou: Methodology (equal); Resources (equal); Supervision (equal); Visualization (equal). Dimitris Davazoglou: Data curation (equal); Supervision (equal). Alexander Chroneos: Funding acquisition (equal); Methodology (equal); Software (equal); Supervision (equal).
DATA AVAILABILITY
The data that support the findings of this study are available within the article.





