Interplay between adsorbates and polarons: CO on rutile TiO2(110)

Michele Reticcioli,1 Igor Sokolovi´c,2 Michael Schmid,2 Ulrike Diebold,2 Martin Setvin,2, ∗ and Cesare Franchini1, † 1University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria 2Institute of Applied Physics, Technische Universit¨atWien, Vienna, Austria (Dated: July 17, 2018) Polaron formation plays a major role in determining the structural, electrical and chemical prop- erties of ionic crystals. Using a combination of first principles calculations and scanning tunneling microscpoy/atomic force microscopy (STM/AFM), we analyze the interaction of polarons with CO molecules adsorbed on the rutile TiO2(110) surface. Adsorbed CO shows attractive coupling with polarons in the surface layer, and repulsive interaction with polarons in the subsurface layer. As a result, CO adsorption depends on the reduction state of the sample. For slightly reduced surfaces, many adsorption configurations with comparable adsorption energies exist and polarons reside in the subsurface layer. At strongly reduced surfaces, two adsorption configurations dominante: either inside an oxygen vacancy, or at surface Ti5c sites, coupled with a surface polaron.

A wide range of materials form polaronic in-gap states unfavorable but might eventually occur at elevated tem- upon injection of extra charge, as the excess electrons peratures [4,8, 24, 26]. or holes couple to the lattice phonon field. The charge The effect of polarons is usually not considered in ad- carrier generated by defects [1–6], doping [7–9], ad- sorption studies. CO adsorption on the rutile (110) sur- sorbates [10–12], or irradiation [13–16], interact with face is a well-studied phenomenon from both theoretical the lattice field to different extents depending on the and experimental points of view [37–43], yet controversies electron-phonon coupling, which is strongly material- appear even in elementary issues. Beyond a general con- dependent [4, 17, 18]. The formation of polarons prevents sensus on the local geometric properties (CO molecules doping-driven insulator-to-metal transition and dramati- adsorb vertically at Ti5c sites at low coverage) [37, 38], cally alters the properties of the system and its function- conflicting outcomes have been reported, which either alities [17, 19, 20]. In the strong short-range coupling suggest [39] or exclude [40–42] the possibility of CO ad- limit, localized (so called small) polarons form; they lo- sorption at VO sites, in the latter case endorsing a major cally distort the lattice and lead to the formation of sharp role played by the fivefold-coordinated Ti5c sites at S0. in-gap states [21, 22]. At low temperature, the ground We show that the apparent disagreements in the litera- state is determined by the polaronic configuration that ture can be resolved by a proper treatment of polarons minimizes the energy of the system [23, 24], but even and their interaction with adsorbates, by combining den- small thermal energies can activate polaron hopping to sity functional theory (DFT) simulations [44, 45] with different hosting sites, thereby changing the nature and scanning tunneling microscopy (STM) and atomic force properties of the polaronic state [9, 25–27]. microscopy (AFM). Formation of polarons is particularly favorable in We start by showing how the presence of CO alters the transition-metal oxides and is further promoted in the stability and orbital topology of polarons. We consider vicinity of the surface, where the crystal lattice is more one CO molecule adsorbed at the Ti5c site next nearest flexible [27–30]. Here, we address the interactions be- neighbor to the VO (NNN-Ti5c) for a VO concentration tween electron polarons and adsorbates. We show of 5.6% (see Methods [46] for details). The spatial ex- that adsorbates are able to alter the stability of po- tension of the S1 and the S0 polaron electronic charge is larons and, in turn, the polarons affect the energet- shown in Fig.1(a) and1(b), respectively. The S1 polaron ics and configuration of the adsorbates. We considered retains the same characteristic spatial distribution as in CO molecules adsorbed on rutile titanium dioxide, an the absence of adsorbates [24, 47], with a dz2 -dx2−y2 -like archetypal polaronic material [31]. In clean (i.e. with- orbital character at the hosting Ti6c site and about 1/3 out CO molecules) TiO2(110) the creation of oxygen va- of the charge distributed on the surrounding atoms. The cancies (VO) promotes the formation of polarons at the CO does not alter this situation, and very little (0.1%) arXiv:1807.05859v1 [cond-mat.mtrl-sci] 16 Jul 2018 sixfold-coordinated Ti6c atoms lying in the subsurface polaronic charge is transfered to the CO molecule. This layer (S1) [23, 24, 32–34]. Due to the attractive interac- can be observed in filled-state STM images [see the weak tion, polarons tend to populate sites near VO [34–36]. Po- circular spots in the inset of Fig.1(a)]. laron hopping from S1 to the surface layer (S0) is highly Conversely, a polaron in the S0 layer is strongly af- fected by the CO molecule [Fig.1(b)]. The S0 polaron acquires a stronger dyz (54%) character as compared to the case with no adsorbates (43%) [24], and a non negli- ∗ [email protected] gible portion of the polaronic charge (1%) is transfered to † [email protected] the 2π* antibonding orbital of the adsorbed CO [48, 49]. 2

FIG. 1. Effects of a CO molecule adsorbed at NNN-Ti5c site on polaronic states at low reduction level (5.6%, i.e., one VO in a 9×2 two-dimensional unit cell). (a,b): electronic charge density of the S1 (a) and S0 (b) polarons in presence of CO. Atoms at the back are depicted by blurred spheres. A top view of the considered configuration is also sketched in each panel. The insets represent the experimental and simulated STM images. (c): Polaron formation energy of S0 and S1 polarons, in case of a TiO2(110) surface with and without adsorbed CO.

This causes the formation of a double-lobed polaronic age), and less frequently at VO sites (marked as CO+VO). cloud above the CO molecule [see the STM images in the The filled-state STM image of the same region [Fig.2(b)] inset of Fig.1(b)], which has a distinctly different shape shows mostly weak circular spots at the CO molecules on from the one shown in Fig.1(a). Remarkably, the CO- Ti5c sites, which we attribute to the electronic cloud of polaron interaction affects strongly the polaron forma- S1 polarons in the vicinity of the CO, similar to those tion energies (EPOL), as shown in Fig.1(c). While in the predicted in the inset of Fig.1(a). There is only one S1 case the polaron formation energy is only marginally intense double-lobed feature in the STM image, sand- S1 destabilized by the CO (∆EPOL = +23 meV), the for- wiched by two oxygen vacancies. We attribute this to a mation of an S0 polaron becomes much more favorable CO at the NNN-Ti5c site, coupled to an S0 polaron (CO- S0 (∆EPOL = −161 meV): The adsorbate changes the po- S0-polaron complex), as predicted in Fig.1(b). We note laronic ground state of the system, an effect that was that the CO adsorbed at VO sites observed by AFM car- overlooked in previous studies [50, 51]. By inspecting ries no in-gap state, and is transparent to the filled-state the contributions to the polaron formation energy (deter- STM imaging. mined by the balance between the electronic energy gain At a more reduced surface (cVO = 14.5% and θCO= due to the electron-phonon coupling and the energy cost 0.15 ML), the AFM image shows that all VO sites are to locally distort the lattice [8]) we find that the presence occupied by CO molecules [Fig.2(c)]. The filled-states of the CO reduces significantly the structural energy cost STM images of the same region show many double-lobed in the S0 polaron case. Moreover, we verified that the CO molecules at Ti5c sites, indicating the formation of formation of an S0 polaron at the CO-adsorption NNN- CO+S0-polaron complexes [Figs.2(d,e)]. Figs.2(d,e) Ti5c site is more favorable than any other configuration, show two consecutive STM images of the same area, by considering different adsorption and polaronic sites. showing diffusion of several CO molecules, associated However, although the CO adsorption reduces the energy with the change of their coupling to S0 and S1 polarons barrier for the S1-to-S0 polaron hopping by 70 meV, the (indicated by arrows in the figure). The increased num- S1-to-S1 hopping remains predominant (see Fig. S1 [46]). ber of the CO+S0-complexes on the strongly reduced Next we focus on the evolution of the CO adsorption surface is in line with the previously reported polaron process as a function of the VO concentration (cVO ) and dynamics [27]: The S1-to-S0 polaron hopping is more CO coverage (θCO), see Figure2. We use AFM imag- favorable due to the repulsive polaron-polaron interac- ing with a CO-terminated tip to locate the adsorbed CO tions in the S1 layer and attraction of polarons to the molecules. Compared to more commonly used empty- VO centers [24]. state STM imaging, AFM disturbs the system less (tip- Finally, by further increasing the CO concentration induced CO hopping can be avoided) and provides better up to θCO = 0.7 ML [52] at highly reduced condition resolution of the adsorption geometry. Complementary (cVO = 14.5%), the AFM image clearly shows that CO to the AFM, filled-state STM imaging provides informa- tends to avoid adsorption sites nearest neighbor to VO tion about the polarons. [Fig.2(f)], and preferentially populates NNN-Ti 5c sites in Figure2(a,b) shows AFM and STM images of a slightly combination with the S0 polaron [double-lobed features reduced (cVO = 5.8%) surface with a low CO coverage in Fig.2(f), particularly intense when the CO is adsorbed (θCO=0.09 ML). The CO molecules are predominantly at a NNN-Ti5c site sandwiched between two VO’s]. CO adsorbed on Ti5c sites (brighter spots in the AFM im- molecules adsorb also at other Ti5c sites, characterized 3

-27 Hz -24 Hz -12 pA 0 pA a Vs = -1 V 5x5 nm2 b

O br Ti 5c

]

0 - 1 1 [ [ 001]

Vs = -1 V 2870 Vs = -1 V 2870 c Vs = 0 V 4x4 nm2 d Vs = -0.8 V e Vs = -0.9 V O Ti br 5c

FIG. 3. Site dependent adsorption energy. A CO molecule explores Ti5c sites at various distances from the oxygen va-

2048 2047 2046 cancy (down pointing triangles), and the VO site (circle), in f 8x8 nm2 g presence of polarons localized at the S1 layer in a reduced slab with cVO = 5.6% (a) and cVO = 17.6% (b). The inset sketches the considered configurations. In addition, we report the case of adsorption at the NNN-Ti5c hosting an S0 polaron O br (up pointing triangle). Ti 5c

from the VO, in accordance with the experiment [see CO+S1 circular spots in Fig.2(b)]. CO adsorptions at VO (CO+VO) or at polaronic NNN-Ti5c (CO+S0- Vs = -1 V Vs = -1 V 1943 1943 polaron complex) sites are essentially degenerate in en- ergy, comparable to the most favorable CO+S1-polaron configurations. In the experiment, the rare occurence FIG. 2. Experimental AFM (blue-white color gradient) and of the double-lobed CO+S0 spots and CO+VO features filled state STM images (yellow-black color) of CO adsorbed [Figs.2(a,b)] originate from a smaller number of avail- on the rutile (110) surfaces at different surface reduction able adsorption sites, as compared to the non-polaronic states. Dashed circles show the positions of VOs, solid cir- cles show VOs with an adsorbed CO molecule, triangles show Ti5c sites [39].

CO molecules adsorbed on Ti5c atoms. Down-pointing tri- For a strongly reduced surface [cVO = 16.7%, angles indicate coupling with S1 polarons, up-pointing tri- Fig.3(b)], adsorption at Ti 5c sites in combination with angles indicate coupling with S0 polarons. (a,b): Low CO S1 polarons becomes significantly less favorable, mainly coverage and low reduction level (θ = 0.09 ML and c = CO VO due to the absence of Ti sites at large distances from the 5.8%). (c,d,e): moderate CO coverage and high reduction 5c VO’s, whereas the CO+S0 and CO+VO configurations level (θCO = 0.15 ML and cVO = 14.5%). The STM im- ages (d,e) were measured sequentially, and the arrows indicate retain their high stability and represent the most stable the diffusion of CO along Ti5c sites, accompanied by polaron solutions. This is in excellent agreement with the STM hopping. (f,g): high CO coverage and high reduction level measurements that show a progressive increase of double-

(θCO = 0.7 ML and cVO = 14.5%). lobed spots arising from CO+S0-polaron complexes with

increasing cVO [see Figs.2(c-g)] combined with a large density of CO+VO features [see Figs.2(c,f)]. by weaker STM signals, due to the long-range tail of the The interplay between polarons and adsorbed CO electronic clouds of CO-S0-polaron complexes at neigh- molecules significantly affects the adsorption scheme; the bor NNN-Ti5c sites. All oxygen vacancies are occupied different CO-polaron couplings affect the adsorption en- by CO molecules, again with no in-gap state found in the ergies, the bonding distances from the surface as well as STM images. the C-O bond length (Fig. S2 [46]). The various polaron- The interpretation of the experimental data is sup- CO coupling schemes reported here are consistent with ported by calculated site-dependent CO adsorption en- reported experimental data on the CO adsorption: In ergies (Eads), see Fig.3. At low reduction [ cVO = 5.6%, temperature programmed desorption (TPD), the TPD Fig.3(a)], the stability of the CO adsorption at non- spectra measured on the rutile (110) surface [43] are com- polaronic Ti5c sites (i.e. CO+S1 configurations, down plex, with multiple desorption peaks. In contrast, the pointing triangles) increases with an increasing distance same experiment performed on the anatase TiO2 (101) 4 surface [53] shows a single desorption peak only. This [4] M. Setvin, C. Franchini, X. Hao, M. Schmid, A. Jan- can be associated with the absence of small polarons at otti, M. Kaltak, C. G. Van De Walle, G. Kresse, and the anatase (101) surface [4], which simplifies the adsorp- U. Diebold, Physical Review Letters 113, 1 (2014). tion scheme in comparison to the polaronic rutile. Sim- [5] P. De´ak,J. Kullgren, and T. Frauenheim, Physica Status Solidi - Rapid Research Letters 8, 583 (2014). ilarly, infrared absorption spectra of CO on the anatase [6] C. Verdi, F. Caruso, and F. Giustino, Nature Commu- (101) surface always exhibit a single absorption peak [53], nications 8, 1 (2017). while the rutile (110) surface shows either one or two vi- [7] C. Crevecoeur and H. De Wit, Journal of Physics and brational states, depending on the reduction level of the Chemistry of Solids 31, 783 (1970). crystal [41]. [8] A. Janotti, C. Franchini, J. B. Varley, G. Kresse, and In summary, by combining first principles calculations C. G. Van de Walle, physica status solidi (RRL) - Rapid Research Letters 7, 199 (2013). and surface sensitive techniques we have elucidated the [9] X. Hao, Z. Wang, M. Schmid, U. Diebold, and C. Fran- key role of the interaction between polarons and CO ad- chini, Physical Review B 91, 085204 (2015). sorbates on rutile TiO2(110). We have shown that CO [10] C. Di Valentin, G. Pacchioni, and A. Selloni, adsorption promotes polaron transfer from subsurface Physical Review Letters 97 (2006), 10.1103/Phys- to surface sites, in particular at highly reduced TiO2 RevLett.97.166803. [11] N. A. Deskins, R. Rousseau, and M. Dupuis, Journal of samples, thus reinforcing the activity of surface Ti5c sites. We have identified three distinct adsorption con- Physical Chemistry C 113, 14583 (2009). [12] H. Sezen, M. Buchholz, A. Nefedov, C. Natzeck, figurations with well-defined characteristics: CO at V O S. Heissler, C. Di Valentin, and C. W¨oll, Scientific Re- sites, CO at Ti5c sites weakly coupled with polarons in ports 4, 3808 (2015). the subsurface (manifested by weak circular features in [13] H. Sezen, H. Shang, F. Bebensee, C. Yang, M. Buch- STM), and strongly coupled CO+S0-polaron complexes holz, A. Nefedov, S. Heissler, C. Carbogno, M. Scheffler, at NNN-Ti5c sites (appearing as double-lobes in STM P. Rinke, and C. W¨oll, Nature Communications 6, 6901 images). The coupling between CO and polarons and (2015). its interaction with the defective lattice effectively in- [14] F. Freytag, G. Corradi, and M. Imlau, Scientific Reports 6, 36929 (2016). fluences the dynamics of the CO+polarons complexes, [15] L. M. Carneiro, S. K. Cushing, C. Liu, Y. Su, P. Yang, enabling breaking and recombination of the CO-S0/S1 A. P. Alivisatos, and S. R. Leone, Nature Materials 16, complexes. Our study delivers a consistent and com- 819 (2017). prehensive picture of the CO adsorption process in an [16] K. Miyata, D. Meggiolaro, M. T. Trinh, P. P. Joshi, archetypal polaronic material, solves long-standing am- E. Mosconi, S. C. Jones, F. De Angelis, and X.-Y. Zhu, biguities and conflicting interpretations of experimental Science Advances 3, e1701217 (2017). results, and sets the path for revisiting the interpreta- [17] A. S. Alexandrov and J. T. Devreese, Advances in Po- laron Physics, Springer Series in Solid-State Sciences tion of adsorption processes in polar semiconductors and (Springer-Verlag, Berlin, 2010). transition metal oxides. [18] F. Giustino, Rev. Mod. Phys. 89, 015003 (2017). [19] I. G. Austin and N. F. Mott, Advances in Physics 50, 757 (2001). [20] A. M. Stoneham, J. Gavartin, A. L. Shluger, A. V. Kim- ACKNOWLEDGEMENTS mel, D. M. Ramo, H. M. Rønnow, G. Aeppli, and C. Ren- ner, Journal of Physics: Condensed Matter 19, 255208 This work was supported by the Austrian Science (2007). Fund (FWF) projects ViCom (F4109-N28) and POLOX [21] A. M. Stoneham, Journal of the Chemical Society, Fara- (Grant No. I 2460-N36), by the ERC Advanced Research day Transactions 2 85, 505 (1989). Grant ‘OxideSurfaces’ and by the FWF Wittgenstein- [22] A. L. Shluger and A. M. Stoneham, Journal of Physics: Condensed Matter 5, 3049 (1993). prize (Z250-N16). The computational results were [23] N. A. Deskins, R. Rousseau, and M. Dupuis, Journal of achieved by using the Vienna Scientific Cluster (VSC). Physical Chemistry C 115, 7562 (2011). [24] M. Reticcioli, M. Setvin, M. Schmid, U. Diebold, and C. Franchini, arXiv 2 (2018). [25] Z. Dohn´alek,I. Lyubinetsky, and R. Rousseau, Progress in Surface Science 85, 161 (2010). [1] E. Finazzi, C. D. Valentin, G. Pacchioni, C. Di Valentin, [26] P. M. Kowalski, M. F. Camellone, N. N. Nair, B. Meyer, and G. Pacchioni, The Journal of Physical Chemistry C and D. Marx, Physical Review Letters 105, 5 (2010). 113, 3382 (2009). [27] M. Reticcioli, M. Setvin, X. Hao, P. Flauger, G. Kresse, [2] X. Mao, X. Lang, Z. Wang, Q. Hao, B. Wen, Z. Ren, M. Schmid, U. Diebold, and C. Franchini, Physical Re- D. Dai, C. Zhou, L.-M. Liu, and X. Yang, The Journal view X 7, 031053 (2017). of Physical Chemistry Letters 4, 3839 (2013). [28] P. Kr¨uger,S. Bourgeois, B. Domenichini, H. Magnan, [3] S. Moser, L. Moreschini, J. Ja´cimovi´c, O. S. Bariˇsi´c, D. Chandesris, P. Le F`evre, A. M. Flank, J. Jupille, H. Berger, A. Magrez, Y. J. Chang, K. S. Kim, A. Bost- L. Floreano, A. Cossaro, A. Verdini, and A. Morgante, wick, E. Rotenberg, L. Forr´o, and M. Grioni, Physical Physical Review Letters 100, 2 (2008). Review Letters 110, 1 (2013). 5

[29] Y. Yoon, Y. Du, J. C. Garcia, Z. Zhu, Z.-t. T. Wang, 4731 (2012). N. G. Petrik, G. A. Kimmel, Z. Dohnalek, M. A. Hen- [42] Y. Cao, M. Yu, S. Qi, S. Huang, T. Wang, M. Xu, S. Hu, derson, R. Rousseau, N. A. Deskins, and I. Lyubinetsky, and S. Yan, Scientific Reports 7, 6148 (2017). ChemPhysChem 16, 313 (2015). [43] N. G. Petrik and G. A. Kimmel, The Journal of Physical [30] T. Shibuya, K. Yasuoka, S. Mirbt, and B. Sanyal, The Chemistry Letters 3, 3425 (2012). Journal of Physical Chemistry C 121, 11325 (2017). [44] G. Kresse and J. Furthm¨uller, Physical Review B 54, [31] B. V. N. and M. D. N., physica status solidi (b) 27, 443 11169 (1996). (1968). [45] G. Kresse and J. Furthm¨uller, Computational Materials [32] U. Diebold, S.-c. Li, and M. Schmid, Annual Review of Science 6, 15 (1996). Physical Chemistry 61, 129 (2010). [46] “Supplementary material,”. [33] M. A. Henderson, Surface Science Reports 66, 185 [47] T. Shibuya, K. Yasuoka, S. Mirbt, and B. Sanyal, Jour- (2011). nal of Physics: Condensed Matter 24, 435504 (2012). [34] T. Shibuya, K. Yasuoka, S. Mirbt, and B. Sanyal, The [48] P. S. Bagus, C. J. Nelin, and C. W. Bauschlicher, Phys- Journal of Physical Chemistry C 118, 9429 (2014). ical Review B 28, 5423 (1983). [35] C. M. Yim, M. B. Watkins, M. J. Wolf, C. L. Pang, [49] K. I. Hadjiivanov and G. N. Vayssilov, in Advances in K. Hermansson, and G. Thornton, Physical Review Let- Catalysis, Vol. 47 (Academic Press, 2002) pp. 307–511. ters 117, 1 (2016). [50] M. Farnesi Camellone, P. M. Kowalski, and D. Marx, [36] P. G. Moses, A. Janotti, C. Franchini, G. Kresse, and Physical Review B 84, 035413 (2011). C. G. Van De Walle, Journal of Applied Physics 119 [51] Y.-Y. Yu and X.-Q. Gong, ACS Catalysis 5, 2042 (2015). (2016), 10.1063/1.4948239. [52] “The high coverage was determined by counting the total [37] U. Diebold, Surface Science Reports 48, 53 (2003). number of Ti5c sites (415) and vacancies (14.5%) occu- [38] J. P. Prates Ramalho, F. Illas, and J. R. B. Gomes, pied by co.”. Physical Chemistry Chemical Physics 19, 2487 (2017). [53] M. Setvin, M. Buchholz, W. Hou, C. Zhang, B. St¨oger, [39] R. Mu, A. Dahal, Z.-T. Wang, Z. Dohn´alek,G. A. Kim- J. Hulva, T. Simschitz, X. Shi, J. Pavelec, G. S. Parkin- mel, N. G. Petrik, and I. Lyubinetsky, The Journal of son, M. Xu, Y. Wang, M. Schmid, C. W¨oll,A. Selloni, Physical Chemistry Letters 8, 4565 (2017). and U. Diebold, The Journal of Physical Chemistry C [40] Y. Zhao, Z. Wang, X. Cui, T. Huang, B. Wang, Y. Luo, 119, 21044 (2015). J. Yang, and J. Hou, Journal of the American Chemical Society 131, 7958 (2009). [41] M. Xu, H. Noei, K. Fink, M. Muhler, Y. Wang, and C. W¨oll, Angewandte Chemie International Edition 51,