Stability of Excess Electrons Introduced by Ti Interstitial in Rutile TiO2

Dec 27, 2016 - The surface chemistry of rutile TiO2(110) has long been ascribed to surface bridging oxygen vacancies (VO) and accompanying excess elec...
2 downloads 12 Views 4MB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

Stability of Excess Electrons Introduced by Ti Interstitial in Rutile TiO(110) Surface 2

Kazuki Morita, Taizo Shibuya, and Kenji Yasuoka J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b09669 • Publication Date (Web): 27 Dec 2016 Downloaded from http://pubs.acs.org on January 1, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Stability of Excess Electrons Introduced by Ti Interstitial in Rutile TiO2(110) Surface Kazuki Morita,† Taizo Shibuya,‡ and Kenji Yasuoka∗,† †Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan ‡IoT Devices Research Laboratories, NEC Corporation, Tsukuba, Ibaraki 305-8501, Japan E-mail: [email protected]

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract The surface chemistry of rutile TiO2 (110) has long been ascribed to surface bridging oxygen vacancies (VO ) and accompanying excess electrons. However, recently there has been debate whether titanium interstitials (Tiint ), a subsurface defect, participates in the surface reactions of TiO2 (110). We used a combination of ab initio molecular dynamics and static density functional theory calculations to systematically investigate the spatial distribution of excess electrons introduced by Tiint . We found that these excess electrons form polarons and that the most stable structure had one polaron located at the Tiint site, and the other three at second-layer Ti sites below five-coordinate Ti sites. This behavior of Tiint providing excess electrons to the near-surface region is similar to that of VO , which suggests that Tiint may contribute to the surface chemical reactions of TiO2 (110).

INTRODUCTION The surfaces of transition metal oxides have attracted much attention because of their suitability for use in a wide variety of applications. 1 In particular, rutile TiO2 (110) has been studied intensively as a model surface. 2 Rutile TiO2 , which is used in photocatalysis, dyesensitized solar cells, and biocompatible materials, 3–5 is an insulator with a band gap of 3 eV. The (110) surface of rutile TiO2 consists of alternating rows of bridging oxygen atoms and five-coordinate titanium atoms. 2 The Ti sites are not equivalent because of the broken symmetry on the rutile TiO2 (110) surface, so we have labeled each site, as shown in Figure 1. A successful model used to describe the chemistry of surfaces is the bridging oxygen vacancy (VO ) model established by Henrich and co-workers. 6 In the VO model, the chemistry of a reduced surface is explained by VO and accompanying excess electrons attributed to the band gap states observed 0.7 eV below the conduction band edge in photoelectron spectroscopy (PES). Despite the simplicity of this model, it has been successfully used to explain many surface chemical reactions. 7–10 2 ACS Paragon Plus Environment

Page 2 of 18

Page 3 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Recent advances in experiment and theory have provided new evidence that the excess electrons introduced by VO form a polaron in the second layer below the Ti5c row (Ti5c-2nd ), 10–15 and that other defects such as the Ti interstitial (Tiint ) and the surface OH can also introduce excess electrons. 11,16–27 Wendt et al. 16 used PES with scanning tunneling microscopy to observe a band gap state on a rutile TiO2 (110) surface without VO . However, if all excess electrons are introduced by the VO , this result cannot be explained. They suggested that excess electrons are mainly withdrawn from Tiint in the near-surface region. More recently, Lira and co-workers used temperature-programmed desorption to measure the change in the amount of adsorbed O2 as they reduced the surface to increase the amount of VO . 22,23 The surface properties were explained by excess electrons introduced by these donors. However, the contribution of Tiint is under debate. 28 The existence of Tiint that diffuse on the rutile TiO2 (110) surface was confirmed more than a decade ago. 29 Recently, it was reported that one Tiint introduces four excess electrons. 20,21,30–34 Thus, both VO and Tiint can provide excess electrons. It is known that the former provides excess electrons to Ti5c-2nd and those excess electrons contribute to the surface chemistry of rutile TiO2 (110). In contrast, the behavior of excess electrons on the rutile TiO2 (110) surface introduced by the latter is yet unclear. Because TiO2 is largely ionic, a growing consensus is that its characteristics can often be explained using electrostatics. 13,35 The question here is: will excess electrons introduced by Tiint overcome the against attraction force and be provided to the surface? In this paper, we use density functional theory (DFT) to calculate the stability of excess electrons at different Ti sites and compare our results with the known behavior of excess electrons introduced by VO . First, we use ab initio molecular dynamics (AIMD) to generate possible polaron configurations and their probability at a finite temperature, and identify Ti sites that excess electrons are more likely to localize on. Next, we focus on these Ti sites and perform structure relaxation using static DFT to closely compare stability at low temperature by removing the effect of thermal fluctuations. We show that, with the exception

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of some excess electrons attracted by Tiint , most excess electrons introduced by Tiint behave in a similar manner to those introduced by VO , which suggests that Tiint contributes to the surface chemical reactions of TiO2 (110).

METHODOLOGY Calculations were performed with a plane-wave basis set using the Vienna Ab initio Simulation Package (VASP). 36 The Perdew–Burke–Ernzerhof functional was employed as an exchange and correlation functional. 37 An on-site coulomb correction, Hubbard U , was applied with the widely used U value of 4.2 eV. 38 Core-valence interaction was described using projector-augmented waves. The Brillouin zone was sampled at the gamma point with a plane-wave expansion cutoff of 450 eV. We adopted double-sided surface slab model to eliminate the effect from broken symmetry on the bulk terminated side. The slab thickness chosen was nine layers to account for the convergence of surface energy and bond length in the deepest layer. 39,40 A (6 × 2) slab model with a 15 Å vacuum between the surfaces was prepared with a Tiint at an octahedral site between the fourth and fifth layers (217 Ti and 432 O atoms in total). All atoms were allowed to move throughout all of the calculations. The AIMD simulation was performed using Born–Oppenheimer molecular dynamics, which strictly followed the Born–Oppenheimer approximation by minimizing the electronic energy for every ionic step. We raised the simulation temperature to 1000 K and enhanced thermal fluctuation to sample enough electron configurations, as conducted by Kowalski et al. 12 An AIMD simulation with a 3 fs time step was performed with a canonical ensemble for 6 ps with temperature control using a Nos´ e–Hoover thermostat. When an electron occupies an empty Ti 3d orbital, it induces a magnetic moment of 1.0 µB . Using this property, we determined the location of excess electrons from the difference between the spin up and spin down components of the projected charge density associated with each Ti site. When this difference exceeds 0.5 µB , at a particular Ti site, an excess electron is considered to be there. The

4 ACS Paragon Plus Environment

Page 4 of 18

Page 5 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

lifetime of excess electrons at a certain Ti site can be associated with the qualitative stability of excess electrons at that site. With the grouping shown in Figure 1, the localization times of excess electrons were summed over time and normalized to give a distribution function. Because our interest was the effect of the surface on the excess electrons, we grouped Ti sites according to their depth from the surface and coordinate of Ti the site was under (Ti5c or Ti6c ). The four rows in the [001] direction closest to Tiint were grouped as “Near Tiint .” Using the calculated distribution function, configurations were chosen and calculated with the static DFT. We prepared different initial structures and compared the final structures obtained following the structure relaxation. To create different terminal structures, we lengthened the bonds around the target Ti sites in the initial structures. It should be noted that changing the initial structures was only to trigger relaxed structures into certain electron structures, and ions were not constrained during the relaxation procedure.

RESULTS AND DISCUSSION We first describe the results of the AIMD simulation, emphasizing similarities and differences between VO studies and our results. Figure 3 depicts the change in location of excess electrons over time. The excess electrons stayed at the initial position near Tiint for first ps (Figure 2(a)), and after reaching equilibrium, they started to localize on various Ti sites. Each excess electron occupied a Ti3d orbital, formally creating Ti3+ . This behavior of excess electrons is in agreement with a previous Car–Parrinello molecular dynamics study performed for the VO . 12 The average lifetime of an excess electron at a Ti site was around 0.1 ps. Together with the results of an additional calculation at 700 K that gave an average lifetime of 0.2 ps, we used the Arrhenius equation to estimate that the activation energy for excess electron hopping was 0.18 eV. This is consistent with the previously calculated activation energy in rutile TiO2 of ∼ 0.3 eV. 41,42 A difference between our results and those of VO studies is the distribution function

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of excess electrons, which is shown on the right-hand side of Figure 2(b). The highest intensity in the distribution function is found for the “Near Tiint ” group followed by Ti5c-2nd and then Ti5c-1st . Although this distribution function mostly agrees with the case for VO , differences arise in the near-field of the defects. The distribution function in the case of VO shows higher intensity near VO ; this is because VO attracts excess electrons. 15,43 In contrast, our distribution function showed substantial population of the deeper “Near Tiint ” region, suggesting that Tiint can also attract some excess electrons. To further examine the stability of excess electrons in the “Near Tiint ” region, we compared the possible configurations. The removal of thermal fluctuation allowed us to simulate the behavior of excess electrons at low temperature, where many experimental results are available. We created different structures in static DFT with excess electrons localized at “Near Tiint ”, Ti5c-2nd , and Ti5c-1st , and then performed structure relaxation. Figure 4 shows four representative relaxed structures A to D, plotted together with the charge density of the band gap states. Figure 4 reveals that all the excess electrons localized in a Ti site are accompanied with lattice distortion, forming polarons. The stability of these structures was calculated and is presented in Table 1. The same configuration of excess electrons was calculated with different values of Hubbard U that were used in previous studies. 44,45 Different U values only resulted a slight change in the relative stability of structure C; this is discussed in detail in the Supporting Information. The difference in stability between “Near Tiint ”, Ti5c-2nd , and Ti5c-1st can be determined by comparing structure A, C, and D, because in all of these structures, one polaron is localized at Tiint and the remaining polarons are at Ti5c-2nd , “Near Tiint ”, and Ti5c-1st , respectively. The energy differences indicate that Ti5c-2nd is more stable than “Near Tiint ” and Ti5c-1st . An electron having more stability at Ti5c-2nd than at “Near Tiint ” seems to contradict the distribution function calculated from the AIMD simulation; however, this can be explained by comparing structure A and B. The difference between structure A and B, which both have excess electrons at Ti5c-2nd , is that A has a polaron at Tiint and B does not. Structure A is more stable than structure B (Table 1), indicating that it is energetically

6 ACS Paragon Plus Environment

Page 6 of 18

Page 7 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

more favorable for Tiint to trap a polaron than to release all four excess electrons. This result implies that the high population of excess electrons in the “Near Tiint ” region observed in the distribution function can be ascribed to trapping of excess electrons by Tiint . It is suggested that the other three excess electrons introduced to the surface behave similarly to the ones introduced by VO and are attracted to VO if it is present. This supports the Tiint -assisted surface chemical reaction model introduced by Lira et al. 22,23 In their model, Tiint provides excess electrons to the surface VO , because their results have shown that the excess electrons solely introduced by VO are insufficient to promote O2 adsorption. Table 1: Total Energies of Different Electronic Structures

Structure A B C D

Polaron location and number near Tiint Ti5c-1st Ti5c-2nd 1 0 3 0 0 4 4 0 0 1 3 0

Relative energy (eV) 0 0.148 0.179 0.291

Our results support those of Wendt et al. 16 They explained the dissociation rate of an O adatom by considering the diffusion of Tiint to the surface. We propose that Tiint may contribute more to surface chemistry if it is closer to the surface. Comparison of different electronic structures in our study revealed that Tiint can trap one excess electron, and provides only three excess electrons to the surface. This implies that with same reduction rate, two VO , which provide a total of four excess electrons to the surface, have more potential to contribute to surface chemical reactions than one Tiint . If the surface preparation and temperature conditions are met, Tiint will diffuse to the surface, allowing trapped excess electrons to participate in surface chemical reactions and increasing the role of Tiint . Yoon et al. 35 reported that Tiint may not exist until the third layer. Thus, it would be illuminating to study the stability of excess electrons introduced by a Tiint located deeper in the slab. Our results are also relevant to H diffusion on the surface. Li et al. proposed a model of H diffusion involving complex formation with an excess electron. 46 In their model, H diffusion

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

is explained by excess electrons moving in the [110] direction, and H on bridging oxygen hopping to the neighboring site following the excess electrons. The H hopping rate was underestimated in their DFT calculation. This is probably because in real systems, more excess electrons are provided from VO or Tiint than were considered in the calculation, which enhances the hopping rate. Although the distribution function indicates that it is rare for excess electrons to visit Ti6c-1st (Ti site under bridging oxygen), the presence of positively charged H+ may attract excess electrons to its nearby sites. Meanwhile, we still consider that VO strongly influence surface chemistry. If a surface is prepared to suppress the formation of Tiint , as done by Yim et al., 9 the chemistry of the resulting surface will be explained by VO . Results for the varying dominance of different defects suggest that the ratio of Tiint to VO in a sample depends on its preparation conditions. Modification to the VO model to consider excess electron donors and excess electrons will expand our knowledge of the surface chemistry of rutile TiO2 (110) surfaces in a more unified matter. Accounting for Tiint as an electron donor will allow us to understand surface chemical reactions on TiO2 more accurately, representing a step towards the development of TiO2 surfaces for specific applications such as photocatalysis.

CONCLUSIONS Explicit calculations with AIMD and static DFT showed that excess electrons introduced by Tiint localize on Ti sites and form polarons on rutile TiO2 (110) surfaces. Out of the four polarons introduced by one Tiint , the most stable structure was when one polaron was at Tiint and other three were at Ti5c-2nd . This result indicates that three quarters of excess electrons are trapped by Ti5c-2nd rather than stabilized near the Tiint region because of the electrostatic attraction of Tiint even without VO . It suggests that the behavior of excess electrons is influenced more strongly by the location, rather than the type of defect it originates from. Excess electrons associated with Tiint could contribute to surface reactions

8 ACS Paragon Plus Environment

Page 8 of 18

Page 9 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

such as O2 adsorption or O adatom diffusion. 16,22 Our results support the modification of the classical VO model to explain surface chemistry to also consider the excess electrons introduced by various electron donors, rather than directly relating band gap states and VO . By generalizing the VO model, a greater variety of surface chemical reactions may be explained for not only TiO2 but also other metal oxides such as ZnO, CeO2 , and SnO2 .

Supporting Information Available A comparison of calculation results obtained using different values of Hubbard U is presented in the Supporting Information. This material is available free of charge via the Internet at http://pubs.acs.org/.

Acknowledgement The authors thank the Supercomputer Center, Institute for Solid State Physics, The University of Tokyo for use of their facilities.

References (1) Diebold, U.; Li, S.-C.; Schmid, M. Oxide Surface Science. Annu. Rev. Phys. Chem. 2010, 61, 129–148. (2) Diebold, U. The Surface Science of Titanium Dioxide. Surf. Sci. Rep. 2003, 48, 53–229. (3) Fujishima, A.; Honda, K. Electrochemical Photolysis of Water at a Semiconductor Electrode. Nature 1972, 238, 37–38. (4) O’Regan, B.; Grätzel, M. A Low-Cost, High-Efficiency Solar Cell Based on DyeSensitized Colloidal TiO2 Films. Nature 1991, 353, 737–740.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(5) Albrektsson, T.; Hansson, H.; Ivarsson, B. Interface Analysis of Titanium and Zirconium Bone Implants. Biomaterials 1985, 6, 97–101. (6) Henrich, V. E.; Dresselhaus, G.; Zeiger, H. J. Observation of Two-Dimensional Phases Associated with Defect States on the Surface of TiO2 . Phys. Rev. Lett. 1976, 36, 1335– 1339. (7) Kimmel, G. A.; Petrik, N. G. Tetraoxygen on Reduced TiO2 (110): Oxygen Adsorption and Reactions with Bridging Oxygen Vacancies. Phys. Rev. Lett. 2008, 100, 196102. (8) Petrik, N. G.; Zhang, Z.; Du, Y.; Dohnálek, Z.; Lyubinetsky, I.; Kimmel, G. A. Chemical Reactivity of Reduced TiO2 (110): The Dominant Role of Surface Defects in Oxygen Chemisorption. J. Phys. Chem. C 2009, 113, 12407–12411. (9) Yim, C.; Pang, C.; Thornton, G. Oxygen Vacancy Origin of the Surface Band-Gap State of TiO2 (110). Phys. Rev. Lett. 2010, 104, 036806. (10) Krüger, P.; Jupille, J.; Bourgeois, S.; Domenichini, B.; Verdini, A.; Floreano, L.; Morgante, A. Intrinsic Nature of the Excess Electron Distribution at the TiO2 (110) Surface. Phys. Rev. Lett. 2012, 108, 126803. (11) Di Valentin, C.; Pacchioni, G.; Selloni, A. Electronic Structure of Defect States in Hydroxylated and Reduced Rutile TiO2 (110) Surfaces. Phys. Rev. Lett. 2006, 97, 166803. (12) Kowalski, P. M.; Camellone, M. F.; Nair, N. N.; Meyer, B.; Marx, D. Charge Localization Dynamics Induced by Oxygen Vacancies on the TiO2 (110) Surface. Phys. Rev. Lett. 2010, 105, 146405. (13) Deskins, N. A.; Rousseau, R.; Dupuis, M. Distribution of Ti3+ Surface Sites in Reduced TiO2 . J. Phys. Chem. C 2011, 115, 7562–7572. (14) Shibuya, T.; Yasuoka, K.; Mirbt, S.; Sanyal, B. A Systematic Study of Polarons Due to

10 ACS Paragon Plus Environment

Page 10 of 18

Page 11 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Oxygen Vacancy Formation at the Rutile TiO2 (110) Surface by GGA+U and HSE06 Methods. J. Phys. : Condens. Matter 2012, 24, 435504. (15) Shibuya, T.; Yasuoka, K.; Mirbt, S.; Sanyal, B. Bipolaron Formation Induced by Oxygen Vacancy at Rutile TiO2 (110) Surfaces. J. Phys. Chem. C 2014, 118, 9429–9435. (16) Wendt, S.; Sprunger, P. T.; Lira, E.; Madsen, G. K. H.; Li, Z.; Hansen, J. O.; Matthiesen, J.; Blekinge-Rasmussen, A.; Laegsgaard, E.; Hammer, B. et al. The Role of Interstitial Sites in the Ti 3d Defect State in the Band Gap of Titania. Science 2008, 320, 1755–1759. (17) Bowker, M.; Bennett, R. A. The Role of Ti3+ Interstitials in TiO2 (110) Reduction and Oxidation. J. Phys. : Condens. Matter 2009, 21, 474224. (18) Deskins, N. A.; Rousseau, R.; Dupuis, M. Localized Electronic States from Surface Hydroxyls and Polarons in TiO2 (110). J. Phys. Chem. C 2009, 113, 14583–14586. (19) Deskins, N. A.; Rousseau, R.; Dupuis, M. Defining the Role of Excess Electrons in the Surface Chemistry of TiO2 . J. Phys. Chem. C 2010, 114, 5891–5897. (20) Mulheran, P. A.; Nolan, M.; Browne, C. S.; Basham, M.; Sanville, E.; Bennett, R. A. Surface and Interstitial Ti Diffusion at the Rutile TiO2 (110) Surface. Phys. Chem. Chem. Phys. 2010, 12, 9763–9771. (21) Papageorgiou, A. C.; Beglitis, N. S.; Pang, C. L.; Teobaldi, G.; Cabailh, G.; Chen, Q.; Fisher, A. J.; Hofer, W. A.; Thornton, G. Electron Traps and Their Effect on the Surface Chemistry of TiO2 (110). Proc. Natl. Acad. Sci. 2010, 107, 2391–2396. (22) Lira, E.; Wendt, S.; Huo, P.; Hansen, J.; Streber, R.; Porsgaard, S.; Wei, Y.; Bechstein, R.; Lægsgaard, E.; Besenbacher, F. The Importance of Bulk Ti3+ Defects in the Oxygen Chemistry on Titania Surfaces. J. Am. Chem. Soc. 2011, 133, 6529–6532.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(23) Lira, E.; Huo, P.; Hansen, J. Ø.; Rieboldt, F.; Bechstein, R.; Wei, Y.; Streber, R.; Porsgaard, S.; Li, Z.; Lægsgaard, E. et al. Effects of the Crystal Reduction State on the Interaction of Oxygen with Rutile TiO2 (110). Catal. Today 2012, 182, 25–38. (24) Mitsuhara, K.; Okumura, H.; Visikovskiy, A.; Takizawa, M.; Kido, Y. The Source of the Ti 3d Defect State in the Band Gap of Rutile Titania (110) Surfaces. J. Chem. Phys. 2012, 136, 124707. (25) Walle, L.; Borg, A.; Uvdal, P.; Sandell, A. Probing the Influence from Residual Ti Interstitials on Water Adsorption on TiO2 (110). Phys. Rev. B 2012, 86, 205415. (26) Deák, P.; Aradi, B.; Frauenheim, T. Oxygen Deficiency in TiO2 : Similarities and Differences Between the Ti Self-Interstitial and the O Vacancy in Bulk Rutile and Anatase. Phys. Rev. B 2015, 92, 045204. (27) Moses, P.; Janotti, A.; Franchini, C.; Kresse, G.; Van de Walle, C. Donor Defects and Small Polarons on the TiO2 (110) Surface. J. Appl. Phys. 2016, 119, 181503. (28) Pang, C. L.; Lindsay, R.; Thornton, G. Structure of Clean and Adsorbate-Covered Single-Crystal Rutile TiO2 Surfaces. Chem. Rev. 2013, 113, 3887–3948. (29) Henderson, M. A. A Surface Perspective on Self-Diffusion in Rutile TiO2 . Surf. Sci. 1999, 419, 174–187. (30) Cho, E.; Han, S.; Ahn, H.-S.; Lee, K.-R.; Kim, S. K.; Hwang, C. S. First-Principles Study of Point Defects in Rutile TiO2-x . Phys. Rev. B 2006, 73, 193202. (31) Finazzi, E.; Di Valentin, C.; Pacchioni, G. Nature of Ti Interstitials in Reduced Bulk Anatase and Rutile TiO2 . J. Phys. Chem. C 2009, 113, 3382–3385. (32) Morgan, B. J.; Watson, G. W. Intrinsic n-Type Defect Formation in TiO2 : A Comparison of Rutile and Anatase from GGA+U Calculations. J. Phys. Chem. C 2010, 114, 2321–2328. 12 ACS Paragon Plus Environment

Page 12 of 18

Page 13 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(33) Lee, H.-Y.; Clark, S. J.; Robertson, J. Calculation of Point Defects in Rutile TiO2 by the Screened-Exchange Hybrid Functional. Phys. Rev. B 2012, 86, 075209. (34) Bjørheim, T. S.; Kuwabara, A.; Norby, T. Defect Chemistry of Rutile TiO2 from First Principles Calculations. J. Phys. Chem. C 2013, 117, 5919–5930. (35) Yoon, Y.; Du, Y.; Garcia, J. C.; Zhu, Z.; Wang, Z.-T.; Petrik, N. G.; Kimmel, G. A.; Dohnalek, Z.; Henderson, M. A.; Rousseau, R. et al. Anticorrelation Between Surface and Subsurface Point Defects and the Impact on the Redox Chemistry of TiO2 (110). ChemPhysChem 2015, 16, 313–321. (36) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Open-Shell Transition Metals. Phys. Rev. B 1993, 48, 13115. (37) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (38) Dohnálek, Z.; Lyubinetsky, I.; Rousseau, R. Thermally-Driven Processes on Rutile TiO2 (110)-(1× 1): A Direct View at the Atomic Scale. Prog. Surf. Sci. 2010, 85, 161–205. (39) Kiejna, A.; Pabisiak, T.; Gao, S. W. The Energetics and Structure of Rutile TiO2 (110). J. Phys. : Condens. Matter 2006, 18, 4207–4217. (40) Bredow, T.; Giordano, L.; Cinquini, F.; Pacchioni, G. Electronic Properties of Rutile TiO2 Ultrathin Films: Odd-Even Oscillations with the Number of Layers. Phys. Rev. B: Condens. Matter Mater. Phys. 2004, 70, 1–6. (41) Deskins, N. A.; Dupuis, M. Electron Transport via Polaron Hopping in Bulk TiO2 : A Density Functional Theory Characterization. Phys. Rev. B 2007, 75, 195212. (42) Janotti, A.; Franchini, C.; Varley, J.; Kresse, G.; Van de Walle, C. Dual Behavior of Excess Electrons in Rutile TiO2 . Phys. Status Solidi RRL 2013, 7, 199–203. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(43) Yim, C.; Watkins, M.; Wolf, M.; Pang, C.; Hermansson, K.; Thornton, G. Engineering Polarons at a Metal Oxide Surface. Phys. Rev. Lett. 2016, 117, 116402. (44) Mattioli, G.; Filippone, F.; Alippi, P.; Bonapasta, A. A. Ab Initio Study of the Electronic States Induced by Oxygen Vacancies in Rutile and Anatase TiO2 . Phys. Rev. B 2008, 78, 241201. (45) Setvin, M.; Franchini, C.; Hao, X.; Schmid, M.; Janotti, A.; Kaltak, M.; Van de Walle, C. G.; Kresse, G.; Diebold, U. Direct View at Excess Electrons in TiO2 Rutile and Anatase. Phys. Rev. Lett. 2014, 113, 086402. (46) Li, S.-C.; Zhang, Z.; Sheppard, D.; Kay, B. D.; White, J.; Du, Y.; Lyubinetsky, I.; Henkelman, G.; Dohnálek, Z. Intrinsic Diffusion of Hydrogen on Rutile TiO2 (110). J. Am. Chem. Soc. 2008, 130, 9080–9088.

14 ACS Paragon Plus Environment

Page 14 of 18

Page 15 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1: Schematic image of the surface slab model used in simulations. Small blue spheres represent titanium and large red spheres represent oxygen. Grouping of Ti sites was performed according to the row and depth of each Ti site. The “Near Tiint ” group was defined as the four rows near Tiint . Ti sites that did not belong to any group were denoted as the “other” group.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

(b)

(a) Ti5c-1st Ti6c-1st Ti5c-2nd Ti6c-2nd Ti5c-3rd Ti6c-3rd Near Tiint Other 0

1000

2000

3000

4000

5000

6000 0.0

0.5

1.0

1.5

Time (fs)

Distribution function of population

Figure 2: (a): Location of excess electrons during the simulation. The size of the circle indicates the population of excess electrons at specific Ti sites. Red, blue, dark green, and light green circles represent four, three, two, and one electron, respectively. Population of excess electrons at specific Ti sites. Grouping was performed as described in Figure 1. (b): Distribution function of the population calculated for each group.

1.2 1 0.8 Population

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 18

0.6 0.4 0.2 0 3000

3200

3400

3600

3800

4000

Time (fs) NearTiint Ti6c−1st Ti5c−1st

Ti5c−2nd Ti6c−2nd Ti5c−3rd

Ti6c−3rd Other

Figure 3: The change in population of excess electrons between 3000 and 4000 fs. The duration of population being around 1 represents the lifetime of excess electrons at that site. Grouping of Ti sites was conducted as defined in Figure 1.

16 ACS Paragon Plus Environment

2.0

Page 17 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4: Representative electronic structures calculated from different initial structures. Electron density corresponding to the band gap state is shown in yellow. Total energies of these structures are given in Table 1. The figures are ordered from the most stable structure A to the least stable structure D.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Graphical TOC Entry

18 ACS Paragon Plus Environment

Page 18 of 18