**Theoretical Studies of Titanium Dioxide for Dye-Sensitized Solar Cell and Photocatalytic Reaction**

Fu-Quan Bai, Wei Li and Hong-Xing Zhang

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/intechopen.68745

#### Abstract

This chapter aims to provide researchers in the field of photovoltaics with the valuable information and knowledge needed to understand the physics and modeling of titanium dioxide for dye-sensitized solar cell and photocatalytic reaction. The electronic band structure of titanium dioxide, the treatment of the excited state of titanium dioxide, the molecular dynamics and ultrafast quantum dynamics simulations, and several promising photocatalytic schemes and important considerations for theoretical study are addressed and reviewed. The advanced computational strategies and methods and optimized models to achieve exact simulation are described and discussed, including first principle calculations, nonadiabatic molecular and quantum dynamics, wave function propagation methods, and surface construction of titanium dioxide. These advanced theoretical investigations have become highly active areas of photovoltaics research and powerful tools for the supplement and prediction of related experimental efforts.

Keywords: electronic band structure, first principle calculation, excited state, molecular modeling, photophysical process, photochemistry and photocatalysis

### 1. Introduction

With the increasing concerns about energy issues and the associated environmental pollution, renewable energy technologies are devoted to being the most promising strategy for sustainable energy supply [1, 2]. Among all of these technologies, including hydro, solar, wind, geothermal heat, and biomass, photovoltaic (PV) technology that converts solar energy into electricity have attracted considerable attention. At present, the solid-state junction devices based on silicon material have dominated PV solar energy converters. However, the light-induced degradation of silicon materials limits the device stability and hinders its effective. Third-generation photovoltaics are

© 2017 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

able to produce high-efficiency photon to electricity conversion devices at a cheaper production cost. This is a consensus that the price is decreased, and the cell efficiency is maintained for the third-generation solar cell. There are unique features for this kind of cells; the cell should be sensitized by a photoactive material among dye, quantum dots (QDs), and perovskite system [3]. Work on sensitized photovoltaics started during the 1970s with the use of organic dyes as the sensitizer [4]. The similar structure of these solar cells is that in the case of the original design: The cell has three primary parts, glass sheet with transparent conducting oxide coating (ITO or FTO) as anode and the counter electrode (CE) on top and bottom, respectively; and TiO2 film deposits on the conductive side of the glass sheet which is then immersed in a mixture of a photosensitizer and a solvent for charge transport. The working principle looks also dynamically analogical. But in quantum dot-sensitized solar cells (QDSSCs), that replaced organic dyes in dye-sensitized solar cells (DSSCs) with inorganic sensitizers, the total performance is promoted by the utilization of nano-sized crystals with a short band gap and a high extinction coefficient. Along with this line of thought, since 2009, the perovskite materials as sensitizers have been used initially and works very well with the solid-state hole transfer material. Until now, its efficiency has touched over 21% [5, 6]. For DSSCs, the photoanode components are the dye sensitizer, a mesoporous semiconducting oxide layer and a transparent conducting oxide (TCO). In DSSCs, the transport of charges (electrons) to the external circuit begins when electrons exit the semiconducting network layer and ends when the redox mediator in the charge transport medium returns them to the sensitizers [7]. The photoelectric anode consists of a TCO substrate with a semiconductor oxide layer (usually TiO2) and a dye sensitizer. In fact, there are two TiO2 layers. The first layer is a slender and blocking layer to inhibit electron recombination with ionized dyes and/or mediators. The second layer is more than 20–30 nm thickness of mesoporous TiO2. The larger surface area of the mesoporous TiO2 area allows a greater amount of dye to be adsorbed on its surface. An electrolyte usually with an iodide/triiodide couple is needed for DSSC. The electrolyte can be in liquid/gel or solid form.

### 2. The electronic band structure of titanium dioxide

A typical DSSC consists of the wide band gap TiO2 semiconductor, the dye is employed to absorb light radiation, and the photo-excited electrons are subsequently transferred to the TiO2 [8]. The fast interfacial electron transfer between the dye and TiO2 could, in principle, minimize the loss of the utilized photon energy as heat through the electron-phonon relaxation pathway. Through this process, the absorption range can be efficiently extended to the visible region, which dominates the solar spectrum. In the many semiconductor metal oxides, titanium dioxide (TiO2) is expected to play an import role in PV devices due to its high chemical and optical stability, nontoxicity, and corrosion resistance, low-cost [9]. TiO2 is a semiconductor material with a band gap of 3.2 eV, corresponding to a wavelength of 390 nm. The anatase structure is preferred over other polymorphs for solar cell applications because of its potentially higher conduction band edge energy and lower recombination rate of electron–hole pairs [10]. This is not the case for excellent UV photocatalysts such as rutile or anatase polymorphs of TiO2 due to the rather high band gap that these materials exhibit (Figure 1). This has prompted chemists and material scientists to search for modifications of

Theoretical Studies of Titanium Dioxide for Dye-Sensitized Solar Cell and Photocatalytic Reaction http://dx.doi.org/10.5772/intechopen.68745 231

able to produce high-efficiency photon to electricity conversion devices at a cheaper production cost. This is a consensus that the price is decreased, and the cell efficiency is maintained for the third-generation solar cell. There are unique features for this kind of cells; the cell should be sensitized by a photoactive material among dye, quantum dots (QDs), and perovskite system [3]. Work on sensitized photovoltaics started during the 1970s with the use of organic dyes as the sensitizer [4]. The similar structure of these solar cells is that in the case of the original design: The cell has three primary parts, glass sheet with transparent conducting oxide coating (ITO or FTO) as anode and the counter electrode (CE) on top and bottom, respectively; and TiO2 film deposits on the conductive side of the glass sheet which is then immersed in a mixture of a photosensitizer and a solvent for charge transport. The working principle looks also dynamically analogical. But in quantum dot-sensitized solar cells (QDSSCs), that replaced organic dyes in dye-sensitized solar cells (DSSCs) with inorganic sensitizers, the total performance is promoted by the utilization of nano-sized crystals with a short band gap and a high extinction coefficient. Along with this line of thought, since 2009, the perovskite materials as sensitizers have been used initially and works very well with the solid-state hole transfer material. Until now, its efficiency has touched over 21% [5, 6]. For DSSCs, the photoanode components are the dye sensitizer, a mesoporous semiconducting oxide layer and a transparent conducting oxide (TCO). In DSSCs, the transport of charges (electrons) to the external circuit begins when electrons exit the semiconducting network layer and ends when the redox mediator in the charge transport medium returns them to the sensitizers [7]. The photoelectric anode consists of a TCO substrate with a semiconductor oxide layer (usually TiO2) and a dye sensitizer. In fact, there are two TiO2 layers. The first layer is a slender and blocking layer to inhibit electron recombination with ionized dyes and/or mediators. The second layer is more than 20–30 nm thickness of mesoporous TiO2. The larger surface area of the mesoporous TiO2 area allows a greater amount of dye to be adsorbed on its surface. An electrolyte usually with an iodide/triiodide couple is needed for DSSC. The electrolyte can be

in liquid/gel or solid form.

230 Titanium Dioxide

2. The electronic band structure of titanium dioxide

A typical DSSC consists of the wide band gap TiO2 semiconductor, the dye is employed to absorb light radiation, and the photo-excited electrons are subsequently transferred to the TiO2 [8]. The fast interfacial electron transfer between the dye and TiO2 could, in principle, minimize the loss of the utilized photon energy as heat through the electron-phonon relaxation pathway. Through this process, the absorption range can be efficiently extended to the visible region, which dominates the solar spectrum. In the many semiconductor metal oxides, titanium dioxide (TiO2) is expected to play an import role in PV devices due to its high chemical and optical stability, nontoxicity, and corrosion resistance, low-cost [9]. TiO2 is a semiconductor material with a band gap of 3.2 eV, corresponding to a wavelength of 390 nm. The anatase structure is preferred over other polymorphs for solar cell applications because of its potentially higher conduction band edge energy and lower recombination rate of electron–hole pairs [10]. This is not the case for excellent UV photocatalysts such as rutile or anatase polymorphs of TiO2 due to the rather high band gap that these materials exhibit (Figure 1). This has prompted chemists and material scientists to search for modifications of

Figure 1. Band structure of rutile and anatase calculated by PBEsol with the respectively relaxed crystal structures. (Copyright 2016 Royal Society of Chemistry).

TiO2 which could exhibit photocatalytic activity under visible light. Doping seemed to be one of the obvious options.

### 3. First principle calculations and modeling simulation

In this section, we will provide a comprehensive insight on the different levels of theoretical methods used in the description of the electronic structures, and particularly attention was paid to methods involved in the excited state in the metal oxide. At present, the wave function-based ab initio method and density functional theory (DFT) are the two main types of theoretical method. However, both methods have drawbacks, which are Born-Oppenheimer (BO) approximation break down when treating the excited state potential energy surface (PES). Of course, the BO approximation is good for description of the ground state properties, such as the geometry optimization, electronic structure calculation, since the nuclei are fixed for the ground state. Stochastic methods need to be developed to properly account for the excited state, but this kind of method is usually time-consuming since it requires the evaluation of the second-derivative of Hessian matrix. In this respect, choosing a particular computational method must take a compromise between accuracy and its feasibility.

#### 3.1. The DFT theory

Quantum mechanics calculations have widely used to address the physical and chemical properties of metal oxide. The level of the theoretical method includes Hartree-Fock (HF) selfconsistent field (SCF), second-order Møller-Plesset (MP2) perturbation theory, coupled cluster (CC) theory, as well as the most popular DFT method. DFT has been the standard "answer" for the electronic structure calculations of a series of material systems, from molecules to solids, and from clusters to periodic systems. The one-electron density matrices can be expressed as follows:

$$\rho\_0(\vec{\mathbf{X}}) = \mathcal{N}\left[\Psi\_0^\*(\vec{\mathbf{X}}\_1, \vec{\mathbf{X}}\_2, \dots, \vec{\mathbf{X}}\_N)\Psi(\vec{\mathbf{X}}\_1, \vec{\mathbf{X}}\_2, \dots, \vec{\mathbf{X}}\_N)d\_{\underset{\mathbf{X}\_2}{\rightarrow}}\dots d\_{\underset{\mathbf{X}\_N}{\rightarrow}}\right] \tag{1}$$

where the subindex zero indicates density originates from the ground state wave function. For one N-electron system, DFT approximates the ground state energy of the system as a function of the electron density, that is to say, the ground state energy can be determined as long as the density function is given. However, the exact expression of density functional is still unknown. There are substantial works devoted to exploring the accurate density functionals. Generally, the density functional is given by:

$$E[\rho] = T\_s[\rho] + V\_{\text{erf}}[r] + V\_{\text{Coulomb}}[\rho] + V\_{\text{XC}}[\rho] \tag{2}$$

Ts½ρ� is the kinetic energy term of the noninteracting electrons, the Vert½r� is the external potential, and VCoulomb represents the Coulomb repulsion. The last term VXC[ρ] accounts for the exchange-correlation effect, which is difficult to get the exact expression. The success of DFT depends on whether we can find out the appropriate approximations to the exchangecorrelation term. Generalized gradient approximation (GGA) and local density approximation (LDA) are mostly used for exchange-correlation functionals in condensed matter systems, and meta-GGA and hybrid functionals are often used in isolated molecule. DFT in the local density approximation (LDA) and various GGA consistently predicts that anatase is more stable than rutile [11], and no significant improvement can be achieved when using more sophisticated DFT approaches, including hybrid functionals (PBE0 and HSE06) [12] and DFT plus the Hubbard U correction (DFT + U) [13, 14], and various dispersion-corrected DFT methods. A similar situation also occurs for the DFT + U method, which predicts rutile to be more stable only if an unphysically large U value is used.

#### 3.2. The treatment of excited state of titanium dioxide

The relative stability of different polymorphic phases of a material is often addressed in terms of the ground state total energy [15, 16], but generally speaking, other degrees of freedom, including, in particular, those of nuclei vibration, may also play a significant role. The contributions of vibration, including the zero-point energy (ZPE) and finite-temperature vibration enthalpy and entropy, can be obtained from the phonon spectrum calculated in the harmonic approximation. We have found the resort that by considering high-order correlations in the adiabatic connection fluctuation–dissipation theory with the random phase approximation (ACFDT-RPA), rutile is correctly predicted to be more stable than anatase, which can be physically attributed to different characters in the electronic band structure of rutile and anatase, including, in particular, that rutile has a smaller band gap than anatase [17]. In the framework of DFT, the most efficient way to approach the excitation problem is its time-dependent formalism (TD-DFT) based on the Runge-Gross theorems, which may be seen as a consequence of the time-dependent Schrödinger equation (TD-SE) to describe the time-dependent phenomena [18]. TD-DFT can be viewed an alternative formulation of time-dependent quantum mechanics but, in contrast to the normal approach that relies on wave functions and on the many-body Schrödinger equation, its basic variable is the one-body electron density. The key exceptions so far that should be mentioned here are the development of extending CIS, TDHF, and TD-DFT to one-dimensional system like a one-dimensional chain of hydrogen molecules by Hirata et al. [19, 20], and the recent implementation of TD-DFT to periodic systems in the CRYSTAL code by Bernasconi et al. [21]. Walker et al. also developed the Turbo TDDFT code to make it is possible to carry out TD-DFT calculations using plane wave basis sets for the large condensed system [22]. However, in spite of using periodic boundary conditions, which is a necessary requirement to date, this formalism is only implemented for finite systems. The introduction of a portion of the Fock exchange is not the only strategy to improve the description of the band gap in the oxide and related materials, where LDA and GGA are even qualitatively ineffective. One common method is to add two parameters to the LDA or GGA exchange-correlation function, as in the Hubbard model, an exchange (J) term between electrons with the same orbital angular momentum, to mimic the effective field Coulomb (U) term [23].

Hessian matrix. In this respect, choosing a particular computational method must take a compro-

Quantum mechanics calculations have widely used to address the physical and chemical properties of metal oxide. The level of the theoretical method includes Hartree-Fock (HF) selfconsistent field (SCF), second-order Møller-Plesset (MP2) perturbation theory, coupled cluster (CC) theory, as well as the most popular DFT method. DFT has been the standard "answer" for the electronic structure calculations of a series of material systems, from molecules to solids, and from clusters to periodic systems. The one-electron density matrices can be expressed as

> , …, XN !

ÞΨðX<sup>1</sup> ! , X<sup>2</sup> !

where the subindex zero indicates density originates from the ground state wave function. For one N-electron system, DFT approximates the ground state energy of the system as a function of the electron density, that is to say, the ground state energy can be determined as long as the density function is given. However, the exact expression of density functional is still unknown. There are substantial works devoted to exploring the accurate density functionals. Generally,

Ts½ρ� is the kinetic energy term of the noninteracting electrons, the Vert½r� is the external potential, and VCoulomb represents the Coulomb repulsion. The last term VXC[ρ] accounts for the exchange-correlation effect, which is difficult to get the exact expression. The success of DFT depends on whether we can find out the appropriate approximations to the exchangecorrelation term. Generalized gradient approximation (GGA) and local density approximation (LDA) are mostly used for exchange-correlation functionals in condensed matter systems, and meta-GGA and hybrid functionals are often used in isolated molecule. DFT in the local density approximation (LDA) and various GGA consistently predicts that anatase is more stable than rutile [11], and no significant improvement can be achieved when using more sophisticated DFT approaches, including hybrid functionals (PBE0 and HSE06) [12] and DFT plus the Hubbard U correction (DFT + U) [13, 14], and various dispersion-corrected DFT methods. A similar situation also occurs for the DFT + U method, which predicts rutile to be more stable

The relative stability of different polymorphic phases of a material is often addressed in terms of the ground state total energy [15, 16], but generally speaking, other degrees of freedom, including, in particular, those of nuclei vibration, may also play a significant role. The contributions of vibration, including the zero-point energy (ZPE) and finite-temperature

, …, XN ! Þd X2 ! …d XN

E½ρ� ¼ Ts½ρ� þ Vert½r� þ VCoulomb½ρ� þ VXC½ρ� ð2Þ

! ð1Þ

mise between accuracy and its feasibility.

ρ0ðX ! Þ ¼ N ð Ψ� <sup>0</sup>ðX<sup>1</sup> ! , X<sup>2</sup> !

the density functional is given by:

only if an unphysically large U value is used.

3.2. The treatment of excited state of titanium dioxide

3.1. The DFT theory

follows:

232 Titanium Dioxide

Among the various possible approaches going beyond the DFT band structure calculations to predict the band gap of oxides, a growing interest is nowadays dedicated to the so-called GW quasiparticle approach (Green's function G and the screened Coulomb interaction W). The quasiparticle approach permits simulation of the ejection or the absorption of an electron from an N-electron state to an N 1 electron state. And then, the direct comparison of the calculated results with photoemission or inverse-photoemission experiments can be performed based on this in principle. The first formulation of the GW method dates back to the work from 1965. However, the approach was not applied to large-scale, numerical electronic structure calculations before the mid-eighties. The resulting band gaps compare much better with experiment. The GW approximation is an approximation made in order to calculate the self-energy of a many-body system of electrons. The approximation is the expansion of the self-energy in terms of the single particle Green's function G and the screened Coulomb interaction W [24]. The whole formalism is rather complex, and a more detailed description can be found in the review paper of Huang and Carter [25]. Nevertheless, the ever-growing computational facilities have made GW calculations indeed attractive in solid-state computational simulations. The GW approach consists of updating both Green's function and the dielectric potential at every iteration. The amplitudes and the excitation energies, E, can be obtained from the Bethe-Salpeter equation (BSE), which can be written as an eigenvalue equation according to

$$E\left[\left(\varepsilon\_p^{QP} - \varepsilon\_h^{QP}\right) - E\right]\mathbf{c}\_{ph} + \sum\_{p'h'} \mathbf{K}\_{ph, p'h'} \mathbf{c}\_{p'h'} = \mathbf{0} \tag{3}$$

where εQP <sup>p</sup> and <sup>ε</sup>QP <sup>h</sup> are the quasiparticle eigenvalues of the particle and the hole, respectively. The sum in Eq. (3) runs over the kernel Kph, <sup>p</sup><sup>0</sup> <sup>h</sup><sup>0</sup> , which can be split into two terms:

$$\mathbf{K}\_{ph,p'h'} = \mathbf{K}\_{ph,p'h'}^{\mathrm{exc}} + \mathbf{K}\_{ph,p'h'}^{\mathrm{dir}} \tag{4}$$

Kexc ph,p<sup>0</sup> <sup>h</sup><sup>0</sup> is the particle�hole exchange operator, which has the usual form:

$$\mathbf{K}\_{ph, p'h'}^{\mathrm{exc}} = \int d\mathbf{x} d\mathbf{x}' \Psi\_p^\*(\mathbf{x}) \Psi\_h(\mathbf{x}) \frac{1}{|\mathbf{x} - \mathbf{x}'|} \Psi\_{p'}(\mathbf{x}) \Psi\_{h'}^\*(\mathbf{x}') \tag{5}$$

Kdir ph,p<sup>0</sup> <sup>h</sup><sup>0</sup> is the screened Coulomb operator, where a term depending on the inverse screened dielectric tensor is applied, similarly to the case for GW:

$$\mathbf{K}\_{ph, p'h'}^{\mathrm{dir}} = \int d\mathbf{x} d\mathbf{x'} \Psi\_p^\*(\mathbf{x}) \Psi\_h(\mathbf{x'}) \eta^{-1}(\mathbf{x}, \mathbf{x'}) \frac{1}{|\mathbf{x} - \mathbf{x'}|} \Psi\_{p'}(\mathbf{x}) \Psi\_{h'}^\*(\mathbf{x'}) \tag{6}$$

In the same framework, neutral excitations (e.g., optical and energy-loss spectra) are also well described today through the solution of the Bethe-Salpeter equation. A self-contained method is able to perform both the starting self-consistent field calculation in the generalized Kohn-Sham (gKS) framework and the many-body perturbation theory (MBPT) post-treatment, and then that can produce GW quasiparticle energies, optical spectra and the dielectric screening, prior to Eqs. (3)–(6) through a solution of the BSE. Entertainingly, Onida et al. have shown that the BSE is intimately related to the linear response formulation of TDDFT [26].

#### 3.3. The molecular dynamics simulation

The extremely small size of these particles makes them ideal candidates for investigation with molecular dynamics (MD) simulation, which has been a long-established simulation technique to extensively investigate structural and dynamic properties of solids and fluids at the atomic scale. A large number of force fields have been reported in the literature for modeling TiO2 polymorphs [27–29]. For example, based on the Hubbard-corrected DFT + U and first-principles molecular dynamics (FPMD) simulations, facet-dependent trapping and dynamics of excess electrons at anatase TiO2 surfaces and aqueous interfaces could be achieved. Whereas no electron trapping is observed on the (101) surface in vacuo, an excess electron at the aqueous (101) interface can trigger water dissociation and become trapped into a stable surface Ti3+-bridging OH complex [30].

#### 4. The titanium dioxide used in DSSC

TiO2 semiconductor is the widely used electron conductor in Gråtzel-type photovoltaic DSSC, because of its low-cost and ease of synthesis. A typical DSSC consists of the wide band gap TiO2 semiconductor, the dye is employed to absorb light radiation, and the photo-excited electrons are subsequently transferred to the TiO2. Two major crystal phases of TiO2: rutile (a) and anatase (b), are shown in Figure 2. The rutile phase is energetically stable since the crystal structure is symmetric. Each Ti atom is sixfold, and coordinated with threefold O atom. The anatase phase has the similar coordination number as rutile but with asymmetrical structure. DFT calculation suggests that the anatase (101) surface is the most active surface since it has unpassivated Ti and O atoms, which has been widely used to simulate the interface electron transfer process [31, 32].

#### 4.1. The rule of titanium dioxide

εQP <sup>p</sup> � <sup>ε</sup>QP h � �

where εQP

234 Titanium Dioxide

Kexc ph,p<sup>0</sup>

Kdir ph,p<sup>0</sup> <sup>p</sup> and <sup>ε</sup>QP

The sum in Eq. (3) runs over the kernel Kph, <sup>p</sup><sup>0</sup>

Kexc ph, p<sup>0</sup> <sup>h</sup><sup>0</sup> ¼ ð dxdx<sup>0</sup> Ψ�

Kdir ph,p<sup>0</sup> <sup>h</sup><sup>0</sup> ¼ ð dxdx<sup>0</sup> Ψ�

3.3. The molecular dynamics simulation

into a stable surface Ti3+-bridging OH complex [30].

4. The titanium dioxide used in DSSC

dielectric tensor is applied, similarly to the case for GW:

h i

Kph,p<sup>0</sup>

<sup>h</sup><sup>0</sup> is the particle�hole exchange operator, which has the usual form:

pðxÞΨhðx<sup>0</sup>

the BSE is intimately related to the linear response formulation of TDDFT [26].

� E

<sup>h</sup><sup>0</sup> <sup>¼</sup> <sup>K</sup>exc ph, p<sup>0</sup>

<sup>p</sup>ðxÞΨhðx<sup>Þ</sup> <sup>1</sup>

<sup>h</sup><sup>0</sup> is the screened Coulomb operator, where a term depending on the inverse screened

<sup>Þ</sup>η�<sup>1</sup> ðx, x<sup>0</sup>

In the same framework, neutral excitations (e.g., optical and energy-loss spectra) are also well described today through the solution of the Bethe-Salpeter equation. A self-contained method is able to perform both the starting self-consistent field calculation in the generalized Kohn-Sham (gKS) framework and the many-body perturbation theory (MBPT) post-treatment, and then that can produce GW quasiparticle energies, optical spectra and the dielectric screening, prior to Eqs. (3)–(6) through a solution of the BSE. Entertainingly, Onida et al. have shown that

The extremely small size of these particles makes them ideal candidates for investigation with molecular dynamics (MD) simulation, which has been a long-established simulation technique to extensively investigate structural and dynamic properties of solids and fluids at the atomic scale. A large number of force fields have been reported in the literature for modeling TiO2 polymorphs [27–29]. For example, based on the Hubbard-corrected DFT + U and first-principles molecular dynamics (FPMD) simulations, facet-dependent trapping and dynamics of excess electrons at anatase TiO2 surfaces and aqueous interfaces could be achieved. Whereas no electron trapping is observed on the (101) surface in vacuo, an excess electron at the aqueous (101) interface can trigger water dissociation and become trapped

TiO2 semiconductor is the widely used electron conductor in Gråtzel-type photovoltaic DSSC, because of its low-cost and ease of synthesis. A typical DSSC consists of the wide band gap

cph <sup>þ</sup><sup>X</sup> p0 h0 Kph,p<sup>0</sup> <sup>h</sup><sup>0</sup> cp<sup>0</sup>

<sup>h</sup> are the quasiparticle eigenvalues of the particle and the hole, respectively.

<sup>h</sup><sup>0</sup> <sup>þ</sup> <sup>K</sup>dir ph,p<sup>0</sup>

> jx � x<sup>0</sup> j

> > <sup>Þ</sup> <sup>1</sup> jx � x<sup>0</sup> j

<sup>h</sup><sup>0</sup> , which can be split into two terms:

Ψp0ðxÞΨ�

h0ðx<sup>0</sup>

Ψp0ðxÞΨ�

h0ðx<sup>0</sup>

<sup>h</sup><sup>0</sup> ¼ 0 ð3Þ

<sup>h</sup><sup>0</sup> ð4Þ

Þ ð5Þ

Þ ð6Þ

So far, solid-state devices based on silicon have dominated photovoltaic solar converters. However, light-induced degradation of silicon materials limits device stability, while the relatively high cost (high purity) of photovoltaic power generation hinders its effective competition with fossil energy. Pure silicon has a little use, but doped silicon is the basis for most semiconductors. An attractive alternative to crystalline silicon photovoltaic devices is a unit made of mesoscopic inorganic semiconductors that can be easily prepared and provide a very low-cost manufacturing prospect. In these devices, inorganic semiconductors are primarily used as electron acceptors and provide direct or tortuous paths for electron transport, and in some cases, they can act as scaffolds for absorbing light collectors. Holes reside in the valence band, a level below the conduction band. Doping with an electron acceptor, an atom which may accept an electron, creates a deficiency of electrons, the same as an excess of holes. Among many semiconductor metal oxides, TiO2 nanomaterials appear to be differentiated candidates because of their high chemical and optical stability, non-toxic, low-cost and corrosion resistance. Anatase TiO2 has a crystal structure corresponding to the tetragonal system, but the distortion of the TiO6 octahedron is slightly larger for the anatase phase. The anatase structure with improved electrical properties is superior to other polymorphs used in solar cells because of its potentially higher conduction band edge energy, and it enables

Figure 2. The crystal structure of rutile (a) and anatase (b) TiO2. (Copyright 2013 American Chemistry Society).

remarkably efficient spatial separation of electron–hole pairs on the submicrosecond time scale. The physical and chemical properties of TiO2 nanocrystals are affected not only by the intrinsic electronic structure but also by their size, shape, texture, doping and surface properties. The smaller effective mass in anatase is consistent with the high mobility, band like conduction observed in anatase crystals. It is also responsible for the very shallow donor energies in anatase. In addition to the large surface area, TiO2 nanomaterials should also have high electron mobility, making it possible to efficiently collect electrons injected into TiO2. Since the defects in the TiO2 nanomaterials can be used as electron captures and present in the grain boundaries at the junctions between the nano-sized particles, it is expected to use a network structure consisting of TiO2 monocrystalline nanowires rather than TiO2 nanoparticles rapid electronic transmission. Besides, the properties of TiO2 also strongly relay on the modifications of the TiO2 material host and on the interactions of TiO2 materials with the environment. The surface modification not only affects the interfacial energy offset but also has significant impact on the charge separation, transport, and recombination processes.

#### 4.2. The charge transfer process and Marcus theory

The Marcus theory is the most fundamental theory for a description of the electron transfer (ET) process in condensed material. Its formula can be given by [33–35]:

$$k(T) = \kappa A \sigma^2 \exp\left\{-\frac{w^r + \frac{\Lambda}{4}\left(1 + \frac{\Delta \tilde{G}^{\circ}}{\Lambda}\right)^2}{k\_B T}\right\} \tag{7}$$

$$
\Delta \tilde{\mathbf{G}}^o = \Delta \mathbf{G}^o + \mathbf{u}^p - \mathbf{u}^r \tag{8}
$$

The above equation explicitly takes various contributions to the free energy into account, where Aσ2represents the collision frequency, and σ is the distance between the reacting species. <sup>Δ</sup>G<sup>ο</sup> is the reaction free energy, <sup>Δ</sup>G<sup>~</sup> <sup>ο</sup> is the free energy after work correction. wr and wp are the work corresponds to the reactant and product state. λ is the reorganization energy, contains the contributions coming from the reactant and solvent.

$$
\lambda = \lambda\_r + \lambda\_{sd} \tag{9}
$$

The reactant and solvent contributions of the reorganization energy are given by:

$$
\lambda\_r = \sum\_i \frac{f\_i^{(r)} f\_i^{(p)}}{f\_i^{(r)} + f\_i^{(p)}} (\Delta q\_i)^2 \tag{10}
$$

$$
\lambda\_{sol} = (\Delta e)^2 \left( \frac{1}{2R\_1} + \frac{1}{2R\_2} - \frac{1}{\sigma} \right) \left( \frac{1}{\varepsilon\_s} - \frac{1}{\varepsilon\_{op}} \right) \tag{11}
$$

where f ðrÞ <sup>i</sup> and f ðpÞ <sup>i</sup> are the force constants of the given normal mode for reactant and product states, respectively; Δq is the changes of the ith nuclear coordinate; and Δe is the amount of transferred charge from reactant to product. ε<sup>s</sup> and εop are dielectric constants of the solvent for static and optical. The original Marcus theory assumes that diabetic PES can correspond to the reactant and product states, and this simplifies the description of the ET.

#### 4.3. The ultrafast quantum dynamics

remarkably efficient spatial separation of electron–hole pairs on the submicrosecond time scale. The physical and chemical properties of TiO2 nanocrystals are affected not only by the intrinsic electronic structure but also by their size, shape, texture, doping and surface properties. The smaller effective mass in anatase is consistent with the high mobility, band like conduction observed in anatase crystals. It is also responsible for the very shallow donor energies in anatase. In addition to the large surface area, TiO2 nanomaterials should also have high electron mobility, making it possible to efficiently collect electrons injected into TiO2. Since the defects in the TiO2 nanomaterials can be used as electron captures and present in the grain boundaries at the junctions between the nano-sized particles, it is expected to use a network structure consisting of TiO2 monocrystalline nanowires rather than TiO2 nanoparticles rapid electronic transmission. Besides, the properties of TiO2 also strongly relay on the modifications of the TiO2 material host and on the interactions of TiO2 materials with the environment. The surface modification not only affects the interfacial energy offset but also has significant impact on the charge separation, transport, and recom-

The Marcus theory is the most fundamental theory for a description of the electron transfer

8 >>><

>>>:

The above equation explicitly takes various contributions to the free energy into account, where Aσ2represents the collision frequency, and σ is the distance between the reacting species.

work corresponds to the reactant and product state. λ is the reorganization energy, contains

f ðrÞ i f ðpÞ i

� � 1

εs � 1 εop � �

f ðrÞ <sup>i</sup> þ f ðpÞ i ðΔqi Þ

The reactant and solvent contributions of the reorganization energy are given by:

<sup>λ</sup><sup>r</sup> <sup>¼</sup> <sup>X</sup> i

> <sup>2</sup> 1 2R<sup>1</sup> þ 1 2R<sup>2</sup> � 1 σ

λsol ¼ ðΔeÞ

wr <sup>þ</sup> <sup>λ</sup>

<sup>4</sup> <sup>1</sup> <sup>þ</sup> <sup>Δ</sup>G<sup>~</sup> <sup>ο</sup> λ � �<sup>2</sup>

9 >>>=

>>>;

<sup>Δ</sup>G<sup>~</sup> <sup>ο</sup> <sup>¼</sup> <sup>Δ</sup>G<sup>ο</sup> <sup>þ</sup> <sup>w</sup><sup>p</sup> � <sup>w</sup><sup>r</sup> <sup>ð</sup>8<sup>Þ</sup>

is the free energy after work correction. wr and wp are the

λ ¼ λ<sup>r</sup> þ λsol ð9Þ

<sup>2</sup> <sup>ð</sup>10<sup>Þ</sup>

ð7Þ

ð11Þ

kBT

bination processes.

236 Titanium Dioxide

4.2. The charge transfer process and Marcus theory

<sup>Δ</sup>G<sup>ο</sup> is the reaction free energy, <sup>Δ</sup>G<sup>~</sup> <sup>ο</sup>

the contributions coming from the reactant and solvent.

(ET) process in condensed material. Its formula can be given by [33–35]:

<sup>k</sup>ðTÞ ¼ <sup>κ</sup>Aσ<sup>2</sup> exp �

Nonadiabatic ET involves the participation of the nuclear, the Fermi's Golden rule is usually employed to evaluate the electron transfer rate KI->F of nonadiabatic ET.

$$k\_{I \to F} = \frac{2\pi}{\hbar} |\langle F|H'|I\rangle|^2 \delta(E\_F - E\_I) \tag{12}$$

This formula originates from the first-order time-dependent perturbation theory (TD-PT). Where I〉 is the initial state, F〉 is the final state, and H<sup>0</sup> is the time-dependent perturbation Hamiltonian due to the nuclei vibration. The advantage of FGR is that ET rate can be determined from electronic structure calculation, that is, the coupling between initial and final state depends on the atomic orbital coefficient. This greatly simplifies the calculation of electron transfer for the chromophore-semiconductor system.

A conceptually similar but rougher approach is proposed by May, in which the rate constant derivates to be:

$$k\_{i \to \{f\}} = 2\overline{I} = \frac{2\pi}{\hbar^2} \overline{N} |\overline{V\_e}|^2 \tag{13}$$

This equation is based on linear response approaches, in which Γ is the absorption line broadening. Ve is average density of states. Labat and coworkers derivate an even more simple expression for the electron transfer time based on the Newns�Anderson approach [36]:

$$
\pi = \frac{658}{\Delta} \tag{14}
$$

Δ is the broadening of the dye donor orbital. Generally, this method gives the electron transfer timescale half of the experimental value [37].

The electron transfer rate constant estimated using Marcus and Fermi's golden rule can be calculated at different levels of theory method. But these methods just focus on the description of the single process. However, the electron transfer process might be subject to the influence of environment. Accurate estimation of the electron transfer rate depends on the exact solution of TD-SE, especially the nuclei part at the quantum mechanics level.

$$i\hbar \frac{\partial \check{\chi}\_i(t, R(t))}{\partial t} = \sum\_j \left\{ \left( T\_{mcl} + E\_j(\mathbb{R}) \right) \delta\_{\vec{\eta}} - \hbar^2 \frac{d\_{\vec{\eta}}^{(2)}}{2M} \right\} \check{\chi}\_j(t, R(t)) \tag{15}$$

It is challenging to treat nuclei quantum mechanics. This is because a given system can contain a nucleus with different mass and speed. Therefore, it is difficult to developing a universe density functional to treat nuclei vibration. The accurate solution of Eq. (15) only fed with the small system based on the grid method. But as the size of the system increases, the above method becomes impractical. The alternative is to treat the nuclei classically, and then a semiclassical expression can be obtained:

$$i\hbar \frac{\partial \check{\chi}\_i\left(t, R(t)\right)}{\partial t} = \sum\_j \left\{ E\_j(\mathbf{R}) \delta\_{ij} - i\hbar^2 \frac{p}{M} d\_{ij}^{(1)} \right\} \check{\chi}\_j\left(t, R(t)\right) \tag{16}$$

In the above equation, the classical momenta p is introduced. Surface hopping method is suitable for mimic the wave-packets branching using a classical trajectory. Fewest switches surface hopping (FSSH) is the most popular semi-classical method for nonadiabatic molecular dynamics. The electron transition rate from state i to state j at the small time interval g(t) can be expressed as follows:

$$\mathbf{g}\_{i \to j}(t) = \max(\mathbf{0}, P\_{i \to j}(t))\tag{17}$$

where

$$P\_{i \to j}(t) = 2 \frac{\int\_{t}^{t + \Delta t} \text{Re}\left(c\_i^\*(t')c\_j(t')\frac{\overrightarrow{P}}{M}\overrightarrow{d}\_{\vec{\eta}}(t')\right)dt'}{c\_i^\*(t)c\_i(t)} \approx 2 \frac{\text{Re}\left(c\_i^\*(t)c\_j(t)\frac{\overrightarrow{P}}{M}\overrightarrow{d}\_{\vec{\eta}}(t')\right)\Delta t}{c\_i^\*(t)c\_j(t)}\tag{18}$$

Here, ci(t) is the time-dependent expansion coefficient. The hopping probability is compared with a random number to determine whether it is hopped to a new state or remain in the current state. Nuclear velocity is rescaled to conserve the total energy. More details of the implementation of FSSH can be found in the literature [38, 39]. It is worthy of noting that classical path approximation (CPA) is employed for treating the large system. CPA assumes that the electron depends on nuclei motion, but a nucleus is independent of the electron, and this is the so-called quantum backreaction problem. CPA is valid in case the heavy nuclei has the large energy than a lighter electron, and this could enable the nuclear evolution is not affected by the electronnuclear interaction. The use of CPA could lead to the large computational saving since the many stochastic surface hops can be achieved along the precomputed MD trajectory.

The original FSSH formulas take the hop-rejection into consideration, and this is given by the following:

$$\mathbf{g}\_{i \to j}(\mathbf{t}) \to \mathbf{g}\_{i \to j}(\mathbf{t}) b\_{i \to j}(\mathbf{t}) \tag{19}$$

$$b\_{i \to j}(t) = \begin{cases} \exp\left(-\frac{E\_j - E\_i}{k\_B T}\right), E\_j > E\_i \\ 1, E\_j \le E\_i \end{cases} \tag{20}$$

The transition probabilities gi!<sup>j</sup> ðtÞ in the above equation are rescaled by the Boltzmann factor, taking the detailed balance between transition up and transition down in energy.

#### 4.4. Wavefunction propagation methods

It is challenging to treat nuclei quantum mechanics. This is because a given system can contain a nucleus with different mass and speed. Therefore, it is difficult to developing a universe density functional to treat nuclei vibration. The accurate solution of Eq. (15) only fed with the small system based on the grid method. But as the size of the system increases, the above method becomes impractical. The alternative is to treat the nuclei classically, and then a semi-

EjðRÞδij � <sup>i</sup>ħ<sup>2</sup> <sup>p</sup>

In the above equation, the classical momenta p is introduced. Surface hopping method is suitable for mimic the wave-packets branching using a classical trajectory. Fewest switches surface hopping (FSSH) is the most popular semi-classical method for nonadiabatic molecular dynamics. The electron transition rate from state i to state j at the small time interval g(t) can be

n o

dt<sup>0</sup>

Re c�

M d ð1Þ ij

χ~j � t, RðtÞ �

ðtÞ ¼ maxð0, Pi!<sup>j</sup>ðtÞÞ ð17Þ

<sup>i</sup> <sup>ð</sup>tÞcjðt<sup>Þ</sup> <sup>P</sup> ! <sup>M</sup> dij ! ðt 0 Þ

c�

� �

ðtÞbi!<sup>j</sup>ðtÞ ð19Þ

Δt

<sup>i</sup> <sup>ð</sup>tÞcjðt<sup>Þ</sup> <sup>ð</sup>18<sup>Þ</sup>

ð16Þ

ð20Þ

j

gi!<sup>j</sup>

� �

<sup>i</sup> <sup>ð</sup>tÞciðt<sup>Þ</sup> <sup>≈</sup> <sup>2</sup>

Here, ci(t) is the time-dependent expansion coefficient. The hopping probability is compared with a random number to determine whether it is hopped to a new state or remain in the current state. Nuclear velocity is rescaled to conserve the total energy. More details of the implementation of FSSH can be found in the literature [38, 39]. It is worthy of noting that classical path approximation (CPA) is employed for treating the large system. CPA assumes that the electron depends on nuclei motion, but a nucleus is independent of the electron, and this is the so-called quantum backreaction problem. CPA is valid in case the heavy nuclei has the large energy than a lighter electron, and this could enable the nuclear evolution is not affected by the electronnuclear interaction. The use of CPA could lead to the large computational saving since the many

The original FSSH formulas take the hop-rejection into consideration, and this is given by the

kBT � �

, Ej > Ei

ðtÞ in the above equation are rescaled by the Boltzmann factor,

<sup>ð</sup>tÞ ! gi!<sup>j</sup>

1, Ej ≤ Ei

bi!<sup>j</sup>ðtÞ ¼ exp � Ej � Ei

classical expression can be obtained:

Pi!<sup>j</sup>ðtÞ ¼ 2

The transition probabilities gi!<sup>j</sup>

ð<sup>t</sup>þΔ<sup>t</sup> t

Re c� i ðt 0 Þcjðt 0 Þ P ! <sup>M</sup> dij ! ðt 0 Þ

c�

stochastic surface hops can be achieved along the precomputed MD trajectory.

gi!<sup>j</sup>

8 ><

>:

taking the detailed balance between transition up and transition down in energy.

expressed as follows:

where

238 Titanium Dioxide

following:

iℏ ∂χ~<sup>i</sup> � t, RðtÞ � <sup>∂</sup><sup>t</sup> <sup>¼</sup> <sup>X</sup> The previous section discussed the semi-classical method for nonadiabatic molecular dynamics. That semi-classical method does treat the nuclear classically. However, in some cases, the nuclear quantum effect cannot be neglected, for example, the evaluation of Franck-Condon factor involved in the electron excitation. Explicitly propagation of the nuclear wave function can achieve this goal. Methods like multiconfiguration time-dependent Hartree (MCTDH) method and quantized Hamiltonian dynamics (QHD) are developed for a different level of accuracy. These methods are based on the propagation of Gaussian wavepacket, which has been proposed to study the quantum effects of nuclear dynamics. The MCTDH method is one of the most accurate methods for wave-packet propagation. The nuclei wave function is expressed as follows:

$$\begin{aligned} \left| \ddot{\chi}\_i(t, \mathsf{R}\_1, \mathsf{R}\_2, \dots \mathsf{R}\_f) \right\rangle \\ = \sum\_{j\_1}^{n\_1} \dots \sum\_{j'}^{n\_{\tilde{f}}} A^1\_{j\_1 \dots j'}(t) |\chi^1\_{j\_1}(t, \mathsf{R}\_1)\rangle \dots |\chi^1\_{j\_{\tilde{f}}}(t, \mathsf{R}\_{\tilde{f}}) \rangle \end{aligned} \tag{21}$$

MCDTH provides the exact and accurate solution to the quantum description of nuclear behaviors. Therefore, it is extremely time-demanding in computation. Its complexity scales as Nf , where fis the number of freedom of the total system, while Nf is the number of the atomic orbital basis function. However, this method is only suitable for the system with limited atoms and has been improved using multilayer MCTDH (ML-MCTDH) technique:

$$|\chi\_i(t, R\_1, R\_2 \dots R\_f)\rangle = \sum\_{j\_1}^{n\_1} \dots \sum\_{j\_p}^{n\_p} A^1\_{j\_1 \dots j\_p}(t) |\chi^1\_{j\_1}(t, R\_1^{col, 1})\rangle \dots |\chi^1\_{j\_p}(t, R\_p^{col, 1})\rangle \tag{22}$$

The concept of collective coordinates (Rc�ll) is introduced. In principle, ML-MCTDH theory can be fed with the system contains several thousands of degrees of freedom [33].

In addition to the MCTDH method, quantum dynamics method based on the propagation of electron wave function is also of interest. Time-dependent electron wave function can be expressed as linear combination of atomic orbitals

$$|\mathfrak{D}(t)\rangle = \sum\_{i,a} B\_{i,a}(t)|i,a\rangle \tag{23}$$

ji, α〉 is the corresponding atomic orbital α of atom i. The expansion coefficient can be expressed as follows:

$$B\_{i,a}(t) = \sum\_{q} Q\_{i,a}^{q} \mathbb{C}\_{q} \exp\left(-\frac{i}{\hbar} E\_{q}t\right) \tag{24}$$

Solving the generalized eigenvalue equation:

$$
\delta H \mathbb{Q}^q = E\_q S \mathbb{Q}^q \tag{25}
$$

S is the overlap matrix in the atomic orbital basis, H is the Hamiltonian matrix which is determined from the extended Hückel theory (EHT) calculations. The extended Hückel calculation is a semi-classical quantum chemistry method. Within the scheme of EHT, the off-diagonal Hamiltonian matrix can be given as follows:

$$H\_{\vec{\eta}} = \mathcal{K} \mathcal{S}\_{\vec{\eta}} \frac{H\_{\vec{\imath}\vec{\imath}} + H\_{\vec{\jmath}\vec{\jmath}}}{2} \tag{26}$$

where K is the Wolfsberg-Helmholtz constant, ~1.75. In EHT, only valence electrons are included; the core electron energies are assumed to be constant.

The population of the electron wave function can be determined from the overlap matrix and atomic orbital basis function,

$$P(t) = \left| \sum\_{i,a}^{MOL} \sum\_{i,\beta} B\_{i,\alpha}^\*(t) B\_{\dot{\beta},\beta}(t) S\_{\alpha,\beta}^{i,j} \right| \tag{27}$$

The combination of theoretical methodology enables the realistic description of the ET process in the large system, and we have this method to investigate the ET in dye/TiO2 interface. We have established the detailed mechanism of electron transfer in DSSC [40].

#### 5. The titanium dioxide used as photocatalyst

Another important application of TiO2 is a photocatalytic reaction. The titanium dioxide (TiO2) has been also widely used as a photocatalyst. Numerous studies have been carried out on the popular photocatalysts of TiO2. The reaction turns out to be very complicated, and its mechanism is still not well understood. It is time to use the theoretical tools to find some important insights out. Photoexcitation of TiO2 leads to the formation of free electrons and holes in CB and VB, respectively. There are several carrier dynamic pathways that might happen. First, an electron in CB recombines non-radiatively with a hole in VB in the form of heat. Second, electron acceptors are reduced by photogenerated electrons; and third, electron donors are oxidized by photogenerated holes. The photochemical reaction may happen once the electron and hole diffuse to the TiO2 crystal surface, and then the photogenerated electron and hole can interact with the O2 or H2O molecules to achieve the oxidization or reduction process. Some typical oxidization and reduction process are shown in Figure 3.

The elementary reaction steps can be expressed as follows:

$$\begin{aligned} \mathrm{TiO}\_{2} &+ hv(\lambda < 390nm) \rightarrow e^{-} + h^{+} \\ h^{+} + H\_{2}O &\rightarrow OH^{\bullet} + H^{+} \\ e^{-} + O\_{2} &\rightarrow O\_{2}^{-} \end{aligned} \tag{28}$$

#### 5.1. The surface of titanium dioxide

The performance of the TiO2-based photocatalyst device can be regulated by the surface molecular adsorption, and particularly, the adsorption of O2 and H2O on TiO2 has been widely discussed. Investigation of adsorption of H2O onto TiO2 mainly focuses on which sites water molecular adsorption is favored, and also the form of water adsorption, molecular or dissociative. Rutile

Theoretical Studies of Titanium Dioxide for Dye-Sensitized Solar Cell and Photocatalytic Reaction http://dx.doi.org/10.5772/intechopen.68745 241

Figure 3. Schematic representation of typical photocatalyst process. (Copyright 2013 American Chemistry Society).

TiO2 (110) surface, the most stable TiO2 polymorph, is usually employed to model the surface adsorption. Experimentally, a water molecule was predicted on being molecularly adsorption onto the stoichiometric rutile TiO2 (110) surface under the low coverage [41–44]. Theoretically, both molecularly and dissociative absorptions have been suggested by the DFT calculation [45–47]. The difference between theory and experiment arising from the exact parameters used in DFT calculation, such as the exchange-correlation function, and from the limitation in the simulation scales, for example, experimental system usually has the diameter up to tens of hundreds of angstrom, which is far beyond the capability DFT can handle at the current stage.

Band-gap engineering is an efficient way to tune the photoactive of the semiconductor. Different heteroatoms like anion nonmetals and cation-metals as well as cation-nonmetals have been used to doping in the semiconductor. Generally, semiconductor doping can introduce gap-state in the band gap, and this leads to the reduction of the optical band gap, improving the photoactive. Electronic structure calculation based on DFT has been used to elucidate the possible effect caused by the dopant. For example, substitute the Ti atom of bulk TiO2 with metal cations, such as Cr and V, is an n-type doping. This will introduce the electron donor states within the band gap, which allows for the reduction of the optical band gap, thus the better performance for utilizing the photon energy, extending the photocatalytic activity into visible part of the solar spectrum [48– 51]. Replace the O with N or C atoms belongs to the p-doping, which introduce the electron acceptor states. Therefore, dopant states have the large extent of overlap between the electronic states of TiO2, this minimizing the problems arising with the localized states in n-type doping.

#### 5.2. The CO2 photocatalytic reaction

calculation is a semi-classical quantum chemistry method. Within the scheme of EHT, the

where K is the Wolfsberg-Helmholtz constant, ~1.75. In EHT, only valence electrons are included;

The population of the electron wave function can be determined from the overlap matrix and

The combination of theoretical methodology enables the realistic description of the ET process in the large system, and we have this method to investigate the ET in dye/TiO2 interface. We

Another important application of TiO2 is a photocatalytic reaction. The titanium dioxide (TiO2) has been also widely used as a photocatalyst. Numerous studies have been carried out on the popular photocatalysts of TiO2. The reaction turns out to be very complicated, and its mechanism is still not well understood. It is time to use the theoretical tools to find some important insights out. Photoexcitation of TiO2 leads to the formation of free electrons and holes in CB and VB, respectively. There are several carrier dynamic pathways that might happen. First, an electron in CB recombines non-radiatively with a hole in VB in the form of heat. Second, electron acceptors are reduced by photogenerated electrons; and third, electron donors are oxidized by photogenerated holes. The photochemical reaction may happen once the electron and hole diffuse to the TiO2 crystal surface, and then the photogenerated electron and hole can interact with the O2 or H2O molecules to achieve the oxidization or reduction process. Some typical oxidization and

TiO<sup>2</sup> þ hvðλ < 390nmÞ ! e� þ h<sup>þ</sup>

The performance of the TiO2-based photocatalyst device can be regulated by the surface molecular adsorption, and particularly, the adsorption of O2 and H2O on TiO2 has been widely discussed. Investigation of adsorption of H2O onto TiO2 mainly focuses on which sites water molecular adsorption is favored, and also the form of water adsorption, molecular or dissociative. Rutile

<sup>h</sup><sup>þ</sup> <sup>þ</sup> <sup>H</sup>2<sup>O</sup> ! OH• <sup>þ</sup> <sup>H</sup><sup>þ</sup>

2

e� þ O<sup>2</sup> ! O�

X i,β B�

Hii þ Hjj

i,αðtÞBj,βðtÞSi,j

α,β � �

<sup>2</sup> <sup>ð</sup>26<sup>Þ</sup>

� <sup>ð</sup>27<sup>Þ</sup>

ð28Þ

Hij ¼ KSij

off-diagonal Hamiltonian matrix can be given as follows:

the core electron energies are assumed to be constant.

PðtÞ ¼ � � � XMOL i,α

5. The titanium dioxide used as photocatalyst

The elementary reaction steps can be expressed as follows:

reduction process are shown in Figure 3.

5.1. The surface of titanium dioxide

have established the detailed mechanism of electron transfer in DSSC [40].

atomic orbital basis function,

240 Titanium Dioxide

From the viewpoints of a theoretical worker, one can determine the every property of a given materials using just its coordinates, and this is exactly the idea of first-principles calculation. The properties include the geometrical structure, the ground state energy, the work function, and the band structure. In order to design the more efficient photocatalyst, it is necessary for theoretical worker knows the basic physics prior to performing first principle calculations.

First, the photocatalyst performance depends on the chosen of the TiO2 surface. Asymmetric anatase (101) surface has the most uncoordinated Ti and O atoms, and thus, it is suitable for the molecular adsorption of CO2, thus the better photocatalyst performance (Figure 4). Second, proper energy alignment between TiO2 and the adsorbed molecules is required. This is due to the built-in driving force largely affect the direction of the ET, and a negative driving force should lead to the inefficient ET process, in order to avoid this issue. The TiO2 CB should be higher than the molecular LUMO orbital, whereas the TiO2 VB should be lower than the molecular HOMO orbital. Third, the interface details between TiO2 and adsorbed molecule affect the time scale of interface charge separation. This is because interface charge transfer can be described using the Fermi's Golden rule if the coupling between the electron/hole donor orbital in TiO2 and acceptor orbital in adsorbed molecule is large, then the ET will be fast. The interface interaction depends on the adsorption details of molecular onto TiO2 surfaces.

The last, the work function is defined as the minimum energy required transferring an electron from the highest filled level of a solid to a point in the vacuum outside the solid surface. The work function is a property of the surface, which has been widely used to measure how difficult the charge transfer is. Screen the materials with the proper work function is required.

#### 5.3. Other main photocatalytic reactions

The (101)-(001) surface heterojunction constructed on polyhedral TiO2 nanocrystals has recently been proposed to be favorable for the efficient electron-hole spatial separation due

Figure 4. Photocatalytic reduction of CO2 on the anatase TiO2 (101) surface (Copyright 2016 American Chemical Society).

to the preferential accumulation of electron and hole on (101) and (001) facets, respectively. The formed free electron and hole can promote reactive oxygen species (ROS) production, which potentially can be used for inactivation of bacteria. For the biological system, the theoretical investigation with electronic level is still a challenging for researchers and here only an example for readers a general understanding. A series of truncated octahedral double cone TiO2 nanocrystals co-planar with (101) and (001) facets are prepared to form (101) to (100) of various ratios for optimizing the electron-(001) surface. All of these polyhedral TiO2 nanocrystals can produce ROS more prominently than spherical TiO2 nanocrystals, exhibiting more bacteria against gram-negative (E. coli) and Gram-positive (Staphylococcus aureus) bacteria under simulated sunlight high antibacterial activity. Surface heterojunction as a new strategy is expected to develop effective antimicrobial nanomaterials [52].
