ORIGINAL RESEARCH article
Electronic Structure and Band Alignments of Various Phases of Titania Using the Self-Consistent Hybrid Density Functional and DFT+U Methods
- 1CNRS, LPCT, UMR 7019, Université de Lorraine, Vandœuvre-lès-Nancy,France
- 2Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, South Korea
To understand, and thereby rationally optimize photoactive interfaces, it is of great importance to elucidate the electronic structures and band alignments of these interfaces. For the first-principles investigation of these properties, conventional density functional theory (DFT) requires a solution to mitigate its well-known bandgap underestimation problem. Hybrid functional and Hubbard U correction are computationally efficient methods to overcome this limitation, however, the results are largely dependent on the choice of parameters. In this study, we employed recently developed self-consistent approaches, which enable non-empirical determination of the parameters, to investigate TiO2 interfacial systems—the most prototypical photocatalytic systems. We investigated the structural, electronic, and optical properties of rutile and anatase phases of TiO2. We found that the self-consistent hybrid functional method predicts the most reliable structural and electronic properties that are comparable to the experimental and high-level GW results. Using the validated self-consistent hybrid functional method, we further investigated the band edge positions between rutile and anatase surfaces in a vacuum and electrolyte medium, by coupling it with the Poisson-Boltzmann theory. This suggests the possibility of a transition from the straddling-type to the staggered-type band alignment between rutile and anatase phases in the electrolyte medium, manifested by the formation of a Stern-like layer at the interfaces. Our study not only confirms the efficacy of the self-consistent hybrid functional method by reliably predicting the electronic structure of photoactive interfaces, but also elucidates a potentially dramatic change in the band edge positions of TiO2 in aqueous electrolyte medium which can extensively affect its photophysical properties.
Titania (TiO2) is one of the most prototypical materials, utilized in a wide range of photocatalytic and photovoltaic applications (Diebold, 2003; Thompson and Yates, 2006; Fujishima et al., 2008; Henderson, 2011; Schneider et al., 2014). Depending on octahedron connectivity, TiO2 has several polymorphs, among which the rutile phase is known to be the ground state and the anatase phase is usually regarded as a better photocatalyst (Luttrell et al., 2014). For mechanistic investigations of the photophysical processes, toward optimization of the photocatalytic or photovoltaic activity, it is important to understand the electronic structures, such as bandgap properties and band edge alignments of the various TiO2 phases. However, due to the sensitivity of electronic structures, to the surface state (e.g., defects or molecular adsorptions, particle size, dielectric environment, etc.) (Diebold, 2003; Yu et al., 2003; Stevanovic et al., 2012; Ronca et al., 2013; Amano et al., 2016; Macounová et al., 2017), a complete understanding of the electronic and optical properties of TiO2 remains elusive, despite a number of experimental and theoretical studies that have been conducted. For example, the relative band edge positions between the rutile and anatase phases are still under debate (Li and Gray, 2007; Scanlon et al., 2013; Zhang et al., 2014; Nolan et al., 2016).
First-principles based density functional theory (DFT) is a total energy theory, which provides a wealth of understanding on the electronic structures of materials, with a reasonable computational cost. However, the (semi-)local approximation of the Kohn-Sham (KS) DFT in describing the exchange-correlation (XC) energy invokes an inevitable problem, in which the bandgaps of semiconductors and insulators are significantly underestimated. The fundamental origin of this problem lies in the fact that the total energy vs. the number of electrons obtained with a (semi-)local XC functional, is not a series of linear segments between integer numbers (Sham and Schlüter, 1983; Anisimov et al., 1991).
To resolve the bandgap underestimation problem of KS-DFT, many theoretical advances have been achieved, including the weighted density approximation (Alonso and Girifalco, 1978), Hubbard U correction (DFT+U) method (Anisimov et al., 1991), self-interaction correction method (Perdew and Zunger, 1981), screened exchange approximation (Bylander and Kleinman, 1990), optimized effective potentials (Grüning et al., 2006), generalized Kohn-Sham (GKS) scheme (Seidl et al., 1996), meta-GGA potentials (Tao et al., 2003), and hybrid functionals (Becke, 1993). In addition to these methods, quasiparticle GW approximation (Hedin, 1965) is regarded as the most reliable and formally accurate method in terms of predicting exact band edge positions of the valence band maximum (VBM) and the conduction band minimum (CBM). However, it also negatively affected by expensive computational cost, which prohibits its routine application in large-scale systems such as surfaces or nanoparticles. On the other hand, the hybrid functional or DFT+U methods provide a more computationally efficient manner to accurately describe the electronic structure (Anisimov et al., 1991; Becke, 1993). By including a portion of the nonlocal Hartree-Fock (HF) exchange or an additional Hubbard-like term, the hybrid functional or DFT+U methods, respectively, reduce the self-interaction error of conventional DFT. However, these methods require an empirical parameter to determine either a HF mixing ratio or a Hubbard U parameter.
To avoid such empiricism in the hybrid functional or DFT+U methods, self-consistent approaches have recently been proposed for the non-empirical determination of the HF mixing ratio (Skone et al., 2014) or the Hubbard U parameter (Cococcioni and de Gironcoli, 2005). The self-consistent hybrid (sc-hybrid) functional method is based on the simplified form of the many-body self-energy under the static approximation, called the static Coulomb hole plus screened exchange (COHSEX), which enables a relation of the HF mixing ratio with the inverse macroscopic dielectric constant (Skone et al., 2014). Based on the linear-response property of the total energy with respect to the occupation number, the Hubbard U parameter can also be determined in a self-consistent manner (Cococcioni and de Gironcoli, 2005).
In this study, we investigated the electronic structures and band edge positions of the rutile and anatase phases of TiO2 by means of these recently developed self-consistent hybrid functional and DFT+U methods. In comparison with the conventional DFT results based on the generalized gradient approximation (GGA) of the XC energy, and with the GW results from previous studies, we assessed the reliability of the self-consistent methods by describing the structural, electronic, and optical properties of bulk rutile and anatase TiO2. We then investigated surface band alignment in rutile (110) and anatase (101) surfaces. Further, as most interesting photocatalytic reactions (e.g., water splitting) occurred in water, we examined the band edge positions in a vacuum and aqueous environment, where the solvation effect was implicitly modeled using the Poisson-Boltzmann theory (Mathew et al., 2014; Mathew and Hennig, 2016).
Self-Consistent Approaches in Hybrid Functional and GGA+U Methods
Hybrid functional (Becke, 1993) includes a portion of the HF exchange energy in addition to the (semi-)local DFT XC energy, where the mixing ratio is controlled by the parameter α. For the particular GGA functional chosen by Perdew, Burke, and Ernzerhof (PBE) (Perdew et al., 1996a), the hybrid PBE functional (PBEh) of the XC energy is written as (Perdew et al., 1996b):
The choice of α = 0.25 gives the standard PBE0 functional, which is best fitted to yield the atomization energies of typical molecules (Adamo and Barone, 1999). More recent studies have related the parameter α with the inverse static dielectric constant of the system (ε∞) by comparing the HF exchange with the screened Coulomb interaction (Marques et al., 2011). Then, since the ε∞ is dependent on the choice of α, a self-consistent scheme was suggested to determine α without empiricism (Skone et al., 2014). From the initial guess of αin = 0.25, self-consistency could be achieved by repeating the calculation of ε∞ using the PBEh functional with αin as an input parameter, until the difference between αin and the resultant is smaller than a specified threshold. We hereafter term the self-consistent PBE0 functional method scPBE0.
On the other hand, DFT+U (Anisimov et al., 1991) solves the problem of the unphysical curvature by adding the correction of a Hubbard-like interaction for the atom I with occupancy of its localized orbital nIσ:
where the last term, , is subtracted to avoid double counting since it is included in both the Hubbard term and the DFT energy. Usually, as a practical choice, a simplified rotationally invariant version introduced by Dudarev et al. is employed (Dudarev et al., 1998):
where the term Ueff = U − J accounts for on-site Coulomb repulsion and a mimicked effects of exchange interaction.
Cococcioni and de Gironcoli suggested a self-consistent approach to determine Ueff using linear response theory (Cococcioni and de Gironcoli, 2005). In this approach, the eigenvalue shift, αI, to localize (or delocalize) the occupation nI of the Hubbard site I is introduced to the total energy of a constrained system:
In the case of molecules or solid systems, the non-linear behavior of total energy as a function of the number of electrons is also induced by the rehybridization of localized orbitals. The calculation of the constrained total energy, by fixing the non-interacting KS potential, is therefore also needed, to subtract this rehybridization effect:
Then, Ueff can be determined as from the interacting and non-interacting density response functions of the system with respect to these constrained problems:
The response functions can be calculated as follows: we first obtained a well-converged self-consistent potential with zero perturbation (αI = 0), and then performed both non-self-consistent and self-consistent calculations under the KS potential by applying a small potential shift (non-zero αI). From the non-self-consistent and self-consistent calculations, we obtained the non-interacting and interacting response functions, respectively.
By iteratively calculating the Ueff and response functions until Ueff converges with a specified threshold, we obtained the self-consistent value of Ueff. In the latter part of this manuscript, we will drop the subscript from Ueff for brevity (hereafter referred to simply as U).
We performed DFT calculations of bulk rutile TiO2, bulk anatase TiO2, and their corresponding surface slab models using the projector augmented wave (PAW) method as implemented in the Vienna Ab-initio Simulation Package (VASP) (Kresse and Furthmüller, 1996; Kresse and Joubert, 1999). A plane-wave expansion of 500 eV was used. The electronic iterations were continued until the energy difference from the previous step became smaller than 10−7 eV, and the geometry optimization was iterated until all atomic forces became smaller than 0.02 eV/Å. We optimized the bulk structures using the PBE functional (Perdew et al., 1996a) with the Γ-centered k-mesh, with grid size of 0.4 Å−1 (6 × 6 × 4 and 5 × 5 × 2 for rutile and anatase phases, respectively). We then employed these PBE-relaxed geometries to obtain the α parameters of scPBE0 and the self-consistent U parameters.
To obtain the α of scPBE0 for each phase, we calculated the static dielectric constants using the perturbation expansion after discretization (PEAD) method (Nunes and Gonze, 2001; Souza et al., 2002). Since convergence of the dielectric constant is slower than that of total energy, we used the grid size of 0.3 Å−1 for the k-mesh in the PEAD calculations.
We employed Quantum Espresso (QE) code (Giannozzi et al., 2009, 2017) to obtain the self-consistent U parameters, while all other calculations were performed using VASP. This was due to the lack of implementation of the calculation of the self-consistent U parameters in the official VASP code. During QE calculations, we used PAW potentials and the PBE functional in a manner analogous to that used for the VASP calculations. We further confirmed the consistency of the calculation results from two different codes by comparing the equilibrium lattice parameters, where the QE calculation used a plane-wave kinetic energy cutoff of 100 Ry and the same grid spacing to the VASP calculations (see Table 1). To minimize the artifact due to the periodic images of the density perturbation, we used 2 × 2 × 2 supercells to obtain carefully converged response functions.
Table 1. Structural properties, bandgaps, and static dielectric constants of bulk rutile and anatase TiO.
Using the convergence criterion of 0.01, we determined both the self-consistent α of the scPBE0 and the self-consistent U parameter of the PBE+U. Using these self-consistent parameters, we then optimized the bulk structures and surface slab models of rutile and anatase TiO2 phases. Both rutile (110) and anatase (101) slabs were composed of 6 TiO2 layers, and the middle two layers were kept fixed during the geometry relaxation to mimic the bulk structure. The lattice parameters of the slab models were taken from those of bulk structures which were optimized using the corresponding level of theory. In addition to the scPBE0 and PBE+U, which uses the self-consistent U parameter (namely, scPBE+U), we further performed conventional PBE and PBE0 calculations for the sake of comparison.
Results and Discussion
Figure 1 shows the evolution of the self-consistent parameters α and U during iteration. For scPBE0, we obtained α values of 0.151 and 0.178 for the rutile and anatase TiO2, respectively. Our α value for the rutile phase is slightly different from the self-consistent value determined by He and Franchini, 0.142 (He and Franchini, 2017). The difference was ascribed to the use of slightly different geometries optimized using different XC functionals. Our value was further compared with the empirical value of 0.159, determined to reproduce the experimental bandgap.
Figure 1. Self-consistent convergence of (A) α of scPBE0 and (B) U of scPBE+U during iterations. Black circles and red triangles denote rutile and anatase phases, respectively. The converged values are indicated by the dashed lines.
We obtained the same self-consistent U values of 2.47 eV for both the rutile and anatase phases. Our value is smaller than the frequently used empirical value of 4.20 eV determined to reproduce experimental spectroscopic data (Morgan and Watson, 2007; Araujo-Lopez et al., 2016).
Structural and Electronic Properties of Rutile and Anatase TiO2
Using the self-consistently determined values of α and U, we investigated the equations of state (EoS) of the rutile and anatase phases of TiO2. Figure 2 shows the EoS of the rutile and anatase phases calculated using PBE, scPBE+U, PBE0, and scPBE0 methods. The lattice parameters, equilibrium volumes, and bulk moduli are shown in Table 1.
Figure 2. Equations of state (EoS) for (A) rutile and (B) anatase phases of TiO2 calculated using PBE (black circles), scPBE+U (red triangles), PBE0 (blue squares), and scPBE0 (green diamonds) methods. The vertical magenta lines represent the experimental volumes of the two phases.
We found that PBE and scPBE+U both overestimated the equilibrium volume and underestimated the bulk modulus, while the hybrid functionals showed better agreement with experimental values. Compared to the conventional PBE0, the self-consistent method of scPBE0 demonstrated a slight improvement.
Table 1 also shows the bandgaps and dielectric constants of rutile and anatase TiO2 calculated using various DFT methods. Unsurprisingly, PBE underestimated the bandgaps. However, even after the inclusion of the Hubbard U correction, the scPBE+U still significantly underestimated the bandgaps, implying that the self-consistently determined U value was not large enough to reproduce the experimental bandgaps. Indeed, a previous study has demonstrated that a fairly large U value of 4.2 eV (Morgan and Watson, 2007) was required to empirically fit the calculated TiO2 bandgap to the experimental spectroscopic data.
In contrast, PBE0 severely overestimated the bandgaps of both the rutile and anatase phases, leading to underestimation of the dielectric constants. The self-consistent approach, which gave smaller α values than the value of 0.25 given by PBE0, resulted in smaller bandgaps and higher dielectric constants. In the case of the rutile phase, the scPBE0-determined bandgap of 3.25 eV is close to the electronic bandgap of 3.30 eV measured experimentally using photoelectron and inverse-photoelectron spectroscopies (PES/IPES) (Tezuka et al., 1994). This value further agrees with the G0W0 calculation results, with values ranging from 3.30 to 3.59 eV (Chiodo et al., 2010; Kang and Hybertsen, 2010; Landmann et al., 2012; Zhu and Gao, 2014; Sun et al., 2015). For the anatase phase, the scPBE0-determined bandgap of 3.75 eV also agrees with the G0W0 calculation results ranging from 3.56 to 3.86 eV (Chiodo et al., 2010; Kang and Hybertsen, 2010; Landmann et al., 2012; Patrick and Giustino, 2012; Zhu and Gao, 2014; Sun et al., 2015). Unfortunately, to the best of our knowledge, no experimental value for the electronic bandgap of anatase TiO2 is available in the literature, although the optical bandgap (which is known to be smaller than the electronic bandgap) was experimentally determined to be 3.47–3.53 eV (Reyes-Coronado et al., 2008; Baldini et al., 2017). It was thus concluded that scPBE0 most reliably predicts the structural and electronic properties of both TiO2 phases.
Band Alignments of Rutile (110) and Anatase (101) Surfaces
For photocatalytic applications, an understanding of the surface band alignments in rutile and anatase TiO2 surfaces is important to elucidate the underlying physics of the charge carrier transfer between two phases. One open question in this regard relates to the relative band edge positions in rutile and anatase TiO2 surfaces (Scanlon et al., 2013). Two possibilities have been suggested, as illustrated in Figure 3: one is the straddling type (Kang et al., 2012), in which both the VBM and CBM of the rutile phase are located in the bandgap region of the anatase phase, and the second is the staggered type (Deák et al., 2011; Scanlon et al., 2013; Ju et al., 2014; Garcia et al., 2015), in which the VBM and CBM of the anatase phase are either up- or down-shifted from the VBM and CBM of the rutile phase.
Figure 3. The proposed types of band alignment between rutile and anatase phases of TiO2. In type-I, also called “straddling type,” the excited electrons and holes are accumulated in the rutile phase. On the other hand, in type-II (“staggered type”), the two charge carriers are preferentially accumulated in different phases according to the sign of the relative band edge shifts.
As such, we investigated the positions of the VBM and CBM of rutile and anatase TiO2 surfaces using various DFT methods. We investigated rutile (110) and anatase (101) surfaces, which are known to be the most stable surfaces for each phase. Before our main results, it is noteworthy to note about the effects of junction between the rutile and anatase TiO2. In its usual mixed-phase form of the TiO2 photocatalyst, the direct contact between two phases gives rise to the electron transfer from anatase to rutile phase, which affects to the band alignment and eventually to the photocatalytic activity (Kawahara et al., 2002). Not only the effects of the junction highly depend on the structure, especially at the interface, of the photocatalyst (Deák et al., 2011), but also, we need to model rutile-anatase interface to take account them. The latter, in terms of the computational cost, is a formidable task in our study with the hybrid functionals. As fundamental properties of their catalytic activities, we have focused on the band alignments between the rutile and anatase TiO2 without the junction effects by using the models of their surfaces separately. Figure 4A shows the surface band edge locations for each phase, aligned with respect to the vacuum level.
Figure 4. Band edge positions of rutile (110) and anatase (101) surfaces calculated using PBE, scPBE+U, PBE0, and scPBE0 methods. (A) Band edge positions in vacuum are aligned with respect to the vacuum potential, while (B) band edge positions in the aqueous electrolyte are aligned with respect to the normal hydrogen electrode (NHE) using the theoretical absolute electrode potential of 4.6 V. The reduction potential (H+/H2) and the oxidation potential (O2/H2O) for water splitting are presented as red and green dashed lines, respectively.
For the rutile (110) surface, PBE gave −7.21 and −5.46 eV for the VBM and CBM, respectively. When the U correction was included using scPBE+U, the VBM shifted downward slightly (−7.25 eV), whereas the CBM was up-shifted by a larger degree due to the increased band gap (−5.27 eV). The VBMs calculated using the hybrid functionals PBE0 (−8.69 eV) and scPBE0 (−8.10 eV) are much deeper than those obtained by PBE and scPBE+U, while the CBMs (−4.59 and −4.94 eV for PBE0 and scPBE0, respectively) are higher than those of the PBE and scPBE+U results.
Migani and co-workers calculated the band alignment of a rutile (110) surface using various approximated variants of quasiparticle GW methods (Migani et al., 2014). In their study, the positions of the VBMs were in the region between −7.3 and −8.8 eV, which is consistent with the experimentally derived VBM alignment. However, their band gaps were largely overestimated, resulting in a large up-shift of the CBMs in comparison to the experimentally derived energies. The position of the VBM from our scPBE0 method agrees reasonably well with the GW results of the same authors (Migani et al., 2014). Considering that the scPBE0-determined bandgap is also in good agreement with the experimentally determined electronic bandgap (as discussed above), it was concluded that the location of the CBM is also reliably predicted by the scPBE0 calculations.
Notably, the scPBE0 results support the straddling-type band edge alignment between the rutile (110) and anatase (101) surfaces. This is also consistent with the results of a previous DFT+U study by Zhang et al. (2015).
Motivated by the knowledge that most photocatalytic applications of TiO2 occur in aqueous medium (e.g., water splitting) (Ge et al., 2016), we further examined changes in the band edge positions in response to the effects of an aqueous electrolyte. We employed the Poisson-Boltzmann model to include the dielectric screening effect of the medium, as implemented in VASPsol code (Mathew et al., 2014; Mathew and Hennig, 2016). We considered an electrolyte consisting of an aqueous solution of monovalent anions and cations in 1M concentrations, by using a relative permittivity of 78.4 and a Debye length of 3 Å.
Figure 4B shows the band alignments of TiO2 aligned with respect to the theoretically determined absolute electrode potential of 4.6 V (Mathew and Hennig, 2016) instead of the vacuum level. We first compared the band edge positions with the reduction potential of H+/H2 (0.0 V vs. NHE) and the oxidation potential of O2/H2O (+1.23 V vs. NHE). From the PBE and scPBE+U results, the VBMs of the rutile (110) surface lie at a lower potential than the oxidation potential of O2/H2O, which is inconsistent with a previous experimental finding that TiO2 is one of the most-active water-splitting photocatalyst materials (Zhao and Liu, 2014; Ge et al., 2016; Gellé and Moores, 2017). In contrast, PBE0 and scPBE0 provided reasonable band edge positions that would energetically allow reduction and oxidation of the water molecule.
Interestingly, in comparison to the results in vacuum, we found that the CBM positions in the rutile and anatase surfaces became almost identical when the electrolyte solvation effect was taken into consideration. To elucidate the origin, we show the bound charge distributions at the water-rutile (110) and water-anatase (101) interfaces in Figure 5. We found that ions tended to accumulate at the water-TiO2 interface, forming a Stern-like layer at both interfaces. However, the water-rutile (110) interface showed a more profound tendency for Stern layer formation in comparison to the water-anatase (101) interface. Consequently, a stronger built-in potential imposed at the water-rutile (110) interface led to greater upward shifts of the band edge positions and CBM positions comparable between the rutile and anatase phases. This implies that the band alignments of TiO2 in an electrolyte can differ from those in vacuum, and that a transition from the straddling-type to the staggered-type band alignment may occur in an electrolyte medium, depending on the salt concentrations.
Figure 5. Bound charge distributions at (A) water-rutile (110) and (B) water-anatase (101) interfaces calculated using Poisson-Boltzmann theory coupled with scPBE0 calculations. The magenta lines depict the positions of the uppermost atoms in the surface slab of TiO2.
In this study, we investigated the bulk properties of rutile and anatase phases of TiO2 as well as their surface electronic properties, using recently developed self-consistent variations of the hybrid functional and DFT+U methods, and with the conventional GGA(PBE) and PBE0 methods. The self-consistent hybrid functional method (i.e., scPBE0) demonstrated the best performance in predicting the bulk structural, elastic, and electronic properties of both phases of TiO2. In particular, the bandgaps of both the rutile and anatase phases of TiO2 were accurately described using the scPBE0 method, and the results were in good agreement with the experimental and/or the highly accurate GW results.
Based on the success of the scPBE0 method, we applied this method to calculations of surface slab models. We investigated the surface band alignment of rutile (110) and anatase (101) surfaces in a vacuum and in aqueous electrolyte, described using the Poisson-Boltzmann theory. Our scPBE0 method not only demonstrated the most reliable calculation of band edge positions, but also elucidated that the two surfaces form a straddling-type alignment with respect to each other in a vacuum. However, in the aqueous electrolyte, the solid-liquid interfacial field differently shifts the band edge locations, resulting in the CBMs of the rutile and anatase surfaces being located at nearly the same level. The difference was ascribed to the different tendencies for Stern-like layer formation at the water-rutile and water-anatase interfaces.
Our study suggests that the scPBE0 method, as an approximation of the GW method, can be a practical replacement for the computationally demanding GW calculations. It also suggests that the surface band alignment in an electrolyte can be different from that in a vacuum, and even implies that a transition from the straddling-type to the staggered-type band alignment might be possible by changing the salt concentration of the electrolyte. Our current work will not only provide a methodological suggestion for theoretical investigations on the electronic structures and band edge positions of large-scale systems, but will also provide useful insight into band alignment of two different TiO2 phases which can be helpful in the design of photoactive interfaces.
WK and HK conceptualized and designed the research. WK and MH performed calculations. WK, MH, and HK analyzed the results. WK wrote the original draft. WK, SL, EL, and HK reviewed and edited the manuscript. SL, EL, and HK supervised the research. WK and MH contributed equally to this work.
This work was supported by the Creative Materials Discovery Program (Grant 2017M3D1A1039378) and the French-Korean STAR Program (No. NRF2017K1A3A1A21013734) through the National Research Foundation of Korea (NRF). We also acknowledge the grant funded by the Ministry of Science and ICT of Korea (MSIT) (No. NRF-2017R1A5A1015365).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are very grateful to the Guest Editors (Mariachiara Pastore, Stefano Caramori, and Marina Freitag) for their invitation to contribute to this article and to the reviewers for their insightful comments.
Alonso, J. A., and Girifalco, L. A. (1978). Nonlocal approximation to the exchange potential and kinetic energy of an inhomogeneous electron gas. Phys. Rev. B 17, 3735–3743. doi: 10.1103/PhysRevB.17.3735
Amano, F., Nakata, M., Yamamoto, A., and Tanaka, T. (2016). Effect of Ti3+ ions and conduction band electrons on photocatalytic and photoelectrochemical activity of rutile titania for water oxidation. J. Phys. Chem. C 120, 6467–6474. doi: 10.1021/acs.jpcc.6b01481
Araujo-Lopez, E., Varilla, L. A., Seriani, N., and Montoya, J. A. (2016). TiO2 anatase's bulk and (001) surface, structural and electronic properties: a DFT study on the importance of Hubbard and van der Waals contributions. Surface Sci. 653, 187–196. doi: 10.1016/j.susc.2016.07.003
Baldini, E., Chiodo, L., Dominguez, A., Palummo, M., Moser, S., Yazdi-Rizi, M., et al. (2017). Strongly bound excitons in anatase TiO2 single crystals and nanoparticles. Nat. Commun. 8:13. doi: 10.1038/s41467-017-00016-6
Burdett, J. K., Hughbanks, T., Miller, G. J., Richardson, J. W., and Smith, J. V. (1987). Structural-electronic relationships in inorganic solids: powder neutron diffraction studies of the rutile and anatase polymorphs of titanium dioxide at 15 and 295 K. J. Am. Chem. Soc. 109, 3639–3646. doi: 10.1021/ja00246a021
Chiodo, L., García-Lastra, J. M., Iacomino, A., Ossicini, S., Zhao, J., Petek, H., et al. (2010). Self-energy and excitonic effects in the electronic and optical properties of TiO2 crystalline phases. Phys. Rev. B 82:045207. doi: 10.1103/PhysRevB.82.045207
Cococcioni, M., and de Gironcoli, S. (2005). Linear response approach to the calculation of the effective interaction parameters in the LDA+U method. Phys. Rev. B 71:035105. doi: 10.1103/PhysRevB.71.035105
Dudarev, S. L., Botton, G. A., Savrasov, S. Y., Humphreys, C. J., and Sutton, A. P. (1998). Electron-energy-loss spectra and the structural stability of nickel oxide: an LSDA+U study. Phys. Rev. B 57, 1505–1509. doi: 10.1103/PhysRevB.57.1505
Garcia, J. C., Nolan, M., and Deskins, N. A. (2015). The nature of interfaces and charge trapping sites in photocatalytic mixed-phase TiO2 from first principles modeling. J. Chem. Phys. 142:024708. doi: 10.1063/1.4905122
Giannozzi, P., Andreussi, O., Brumme, T., Bunau, O., Buongiorno Nardelli, M., Calandra, M., et al. (2017). Advanced capabilities for materials modelling with Quantum ESPRESSO. J. Phys. Condens. Matter 29:465901. doi: 10.1088/1361-648X/aa8f79
Giannozzi, P., Baroni, S., Bonini, N., Calandra, M., Car, R., Cavazzoni, C., et al. (2009). QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phy. Condens. Matter 21:395502. doi: 10.1088/0953-8984/21/39/395502
Grüning, M., Marini, A., and Rubio, A. (2006). Density functionals from many-body perturbation theory: the band gap for semiconductors and insulators. J. Chem. Phys. 124:154108. doi: 10.1063/1.2189226
He, J., and Franchini, C. (2017). Assessing the performance of self-consistent hybrid functional for band gap calculation in oxide semiconductors. J. Phys. Condens. Matter 29:454004. doi: 10.1088/1361-648X/aa867e
Ju, M.-G., Sun, G., Wang, J., Meng, Q., and Liang, W. (2014). Origin of high photocatalytic properties in the mixed-phase TiO2: a first-principles theoretical study. ACS Appl. Mater. Interfaces 6, 12885–12892. doi: 10.1021/am502830m
Kang, J., Wu, F., Li, S.-S., Xia, J.-B., and Li, J. (2012). Calculating band alignment between materials with different structures: the case of anatase and rutile titanium dioxide. J. Phys. Chem. C 116, 20765–20768. doi: 10.1021/jp3067525
Kawahara, T., Konishi, Y., Tada, H., Tohge, N., Nishii, J., and Ito, S. (2002). A patterned TiO2(anatase)/TiO2(rutile) bilayer-type photocatalyst: effect of the anatase/rutile junction on the photocatalytic activity. Angew. Chem. 41, 2811–3. doi: 10.1002/1521-3773(20020802)41:15<2811::AID-ANIE2811>3.0.CO;2-%23
Landmann, M., Rauls, E., and Schmidt, W. G. (2012). The electronic structure and optical response of rutile, anatase and brookite TiO2. J. Phys. Condens. Matter 24:195503. doi: 10.1088/0953-8984/24/19/195503
Li, G., and Gray, K. A. (2007). The solid–solid interface: explaining the high and unique photocatalytic reactivity of TiO2-based nanocomposite materials. Chem. Phys. 339, 173–187. doi: 10.1016/j.chemphys.2007.05.023
Luttrell, T., Halpegamage, S., Tao, J., Kramer, A., Sutter, E., and Batzill, M. (2014). Why is anatase a better photocatalyst than rutile?—Model studies on epitaxial TiO2 films. Sci. Rep. 4:4043. doi: 10.1038/srep04043
Macounová, K. M., Klusáčková, M., Nebel, R., Zukalová, M., Klementová, M., Castelli, I. E., et al. (2017). Synergetic surface sensitivity of photoelectrochemical water oxidation on TiO2 (anatase) electrodes. J. Phys. Chem. C 121, 6024–6032. doi: 10.1021/acs.jpcc.6b09289
Mathew, K., and Hennig, R. G. (2016). Implicit Self-Consistent Description of Electrolyte in Plane-Wave Density-Functional Theory. ArXiv e-prints. Available online at: https://ui.adsabs.harvard.edu/#abs/2016arXiv160103346M (Accessed January 1, 2016).
Mathew, K., Sundararaman, R., Letchworth-Weaver, K., Arias, T. A., and Hennig, R. G. (2014). Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways. J. Chem. Phys. 140:084106. doi: 10.1063/1.4865107
Nolan, M., Deskins, N. A., Schwartzenberg, K. C., and Gray, K. A. (2016). Local interfacial structure influences charge localization in titania composites: beyond the band alignment paradigm. J. Phys. Chem. C 120, 1808–1815. doi: 10.1021/acs.jpcc.5b12326
Reyes-Coronado, D., Rodríguez-Gattorno, G., Espinosa-Pesqueira, M. E., Cab, C., de Coss, R., and Oskam, G. (2008). Phase-pure TiO2 nanoparticles: anatase, brookite and rutile. Nanotechnology 19:145605. doi: 10.1088/0957-4484/19/14/145605
Ronca, E., Pastore, M., Belpassi, L., Tarantelli, F., and De Angelis, F. (2013). Influence of the dye molecular structure on the TiO2 conduction band in dye-sensitized solar cells: disentangling charge transfer and electrostatic effects. Energy Environ. Sci. 6, 183–193. doi: 10.1039/C2EE23170K
Scanlon, D. O., Dunnill, C. W., Buckeridge, J., Shevlin, S. A., Logsdail, A. J., Woodley, S. M., et al. (2013). Band alignment of rutile and anatase TiO2. Nat. Mater. 12, 798–801. doi: 10.1038/nmat3697
Schneider, J., Matsuoka, M., Takeuchi, M., Zhang, J., Horiuchi, Y., Anpo, M., et al. (2014). Understanding TiO2 photocatalysis: mechanisms and materials. Chem. Rev. 114, 9919–9986. doi: 10.1021/cr5001892
Stevanovic, A., Büttner, M., Zhang, Z., and Yates, J. T. (2012). Photoluminescence of TiO2: effect of UV light and adsorbed molecules on surface band structure. J. Am. Chem. Soc. 134, 324–332. doi: 10.1021/ja2072737
Sun, H., Mowbray, D. J., Migani, A., Zhao, J., Petek, H., and Rubio, A. (2015). Comparing quasiparticle H2O level alignment on anatase and rutile TiO2. ACS Catal. 5, 4242–4254. doi: 10.1021/acscatal.5b00529
Tao, J., Perdew, J. P., Staroverov, V. N., and Scuseria, G. E. (2003). Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids. Phys. Rev. Lett. 91:146401. doi: 10.1103/PhysRevLett.91.146401
Tezuka, Y., Shin, S., Ishii, T., Ejima, T., Suzuki, S., and Sato, S. (1994). Photoemission and Bremsstrahlung Isochromat spectroscopy studies of TiO2 (rutile) and SrTiO3. J. Phys. Soc. Jpn. 63, 347–357. doi: 10.1143/jpsj.63.347
Yu, J., Yu, J. C., Leung, M. K. P., Ho, W., Cheng, B., Zhao, X., et al. (2003). Effects of acidic and basic hydrolysis catalysts on the photocatalytic activity and microstructures of bimodal mesoporous titania. J. Catal. 217, 69–78. doi: 10.1016/S0021-9517(03)00034-4
Zhang, J., Zhou, P., Liu, J., and Yu, J. (2014). New understanding of the difference of photocatalytic activity among anatase, rutile and brookite TiO2. Phys. Chem. Chem. Phys. 16, 20382–20386. doi: 10.1039/C4CP02201G
Keywords: photocatalyst, titania, band alignment, self-consistent GGA+U, self-consistent hybrid functional
Citation: Kim WJ, Han MH, Lebègue S, Lee EK and Kim H (2019) Electronic Structure and Band Alignments of Various Phases of Titania Using the Self-Consistent Hybrid Density Functional and DFT+U Methods. Front. Chem. 7:47. doi: 10.3389/fchem.2019.00047
Received: 29 November 2018; Accepted: 18 January 2019;
Published: 07 February 2019.
Edited by:Stefano Caramori, University of Ferrara, Italy
Reviewed by:Sergio Tosoni, Università degli Studi di Milano Bicocca, Italy
Zhimin Ao, Guangdong University of Technology, China
Copyright © 2019 Kim, Han, Lebègue, Lee and Kim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work