Modeling polarons in density functional theory: lessons learned from TiO2

Density functional theory (DFT) is nowadays one of the most broadly used and successful techniques to study the properties of polarons and their effects in materials. Here, we systematically analyze the aspects of the theoretical calculations that are crucial to obtain reliable predictions in agreement with the experimental observations. We focus on rutile TiO2, a prototypical polaronic compound, and compare the formation of polarons on the (110) surface and subsurface atomic layers. As expected, the parameter U used to correct the electronic correlation in the DFT + U formalism affects the resulting charge localization, local structural distortions and electronic properties of polarons. Moreover, the polaron localization can be driven to different sites by strain: due to different local environments, surface and subsurface polarons show different responses to the applied strain, with impact on the relative energy stability. An accurate description of the properties of polarons is key to understand their impact on complex phenomena and applications: as an example, we show the effects of lattice strain on the interaction between polarons and CO adsorbates

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

Polaron formation plays a major role in determining the structural, electrical and chemical properties 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 TiO$_2$(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 Ti$_{5c}$ sites, coupled with a surface polaron

Formation and dynamics of small polarons on the rutile TiO2_2(110) surface

Charge trapping and formation of polarons is a pervasive phenomenon in transition metal oxide compounds, in particular at the surface, affecting fundamental physical properties and functionalities of the hosting materials. Here we investigate via first-principle techniques the formation and dynamics of small polarons on the reduced surface of titanium dioxide, an archetypal system for polarons, for a wide range of oxygen vacancy concentrations. We report how the essential features of polarons can be adequately accounted in terms of few quantities: the local structural and chemical environment, the attractive interaction between negatively charged polarons and positively charged oxygen vacancies, and the spin-dependent polaron-polaron Coulomb repulsion. We combined molecular dynamics simulations on realistic samples derived from experimental observations with simplified static models, aiming to disentangle the various variables at play. We find that depending on the specific trapping site, different types of polarons can be formed, with distinct orbital symmetries and different degree of localization and structural distortion. The energetically most stable polaron site is the subsurface Ti atom below the undercoordinated surface Ti atom, owing to a small energy cost to distort the lattice and a suitable electrostatic potential. Polaron-polaron repulsion and polaron-vacancy attraction determine the spatial distribution of polarons as well as the energy of the polaronic in-gap state. In the range of experimentally reachable oxygen vacancy concentrations the calculated data are in excellent agreement with observations, thus validating the overall interpretation

Role of Polarons in Single-Atom Catalysts: Case Study of Me1 [Au1, Pt1, and Rh1] on TiO2(110)

The local environment of metal-oxide supported single-atom catalysts plays a decisive role in the surface reactivity and related catalytic properties. The study of such systems is complicated by the presence of point defects on the surface, which are often associated with the localization of excess charge in the form of polarons. This can affect the stability, the electronic configuration, and the local geometry of the adsorbed adatoms. In this work, through the use of density functional theory and surface-sensitive experiments, we study the adsorption of Rh-1, Pt-1, and Au-1 metals on the reduced TiO2(110) surface, a prototypical polaronic material. A systematic analysis of the adsorption configurations and oxidation states of the adsorbed metals reveals different types of couplings between adsorbates and polarons. As confirmed by scanning tunneling microscopy measurements, the favored Pt-1 and Au-1 adsorption at oxygen vacancy sites is associated with a strong electronic charge transfer from polaronic states to adatom orbitals, which results in a reduction of the adsorbed metal. In contrast, the Rh-1 adatoms interact weakly with the excess charge, which leaves the polarons largely unaffected. Our results show that an accurate understanding of the properties of single-atom catalysts on oxide surfaces requires a careful account of the interplay between adatoms, vacancy sites, and polarons

Competing electronic states emerging on polar surfaces

Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO2 termination of KTaO3(001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions

Doping-induced insulator-metal transition in the Lifshitz magnetic insulator NaOsO3.

By means of first principles schemes based on magnetically constrained density functional theory and on the band unfolding technique we study the effect of doping on the conducting behaviour of the Lifshitz magnetic insulator NaOsO3. Electron doping is treated within a supercell approach by replacing sodium with magnesium at different concentrations ([Formula: see text], [Formula: see text]). Undoped NaOsO3 is subjected to a temperature-driven Lifshitz transition involving a continuous closing of the gap due to longitudinal and rotational spin fluctuations (Kim et al 2016 Phys. Rev. B 94 241113). Here we find that Mg doping suppresses the insulating state, gradually drives the system to a metallic state (via an intermediate bad metal phase) and the transition is accompanied by a progressive lowering of the Os magnetic moment. We inspected the role of longitudinal spin fluctuations by constraining the amplitude of the local Os moments and found that a robust metal state can be achieved below a critical moment. In analogy with the undoped case we conjecture that the decrease of the local moment can be controlled by temperature effects, in accordance with the theory of itinerant electron magnetism

Automated Real-Space Lattice Extraction for Atomic Force Microscopy Images

Analyzing atomically resolved images is a time-consuming process requiring solid experience and substantial human intervention. In addition, the acquired images contain a large amount of information such as crystal structure, presence and distribution of defects, and formation of domains, which need to be resolved to understand a material's surface structure. Therefore, machine learning techniques have been applied in scanning probe and electron microscopies during the last years, aiming for automatized and efficient image analysis. This work introduces a free and open source tool (AiSurf: Automated Identification of Surface Images) developed to inspect atomically resolved images via Scale-Invariant Feature Transform (SIFT) and Clustering Algorithms (CA). AiSurf extracts primitive lattice vectors, unit cells, and structural distortions from the original image, with no pre-assumption on the lattice and minimal user intervention. The method is applied to various atomically resolved non-contact atomic force microscopy (AFM) images of selected surfaces with different levels of complexity: anatase TiO2(101), oxygen deficient rutile TiO2(110) with and without CO adsorbates, SrTiO3(001) with Sr vacancies and graphene with C vacancies. The code delivers excellent results and has proved to be robust against atom misclassification and noise, thereby facilitating the interpretation scanning probe microscopy images

Surface chemistry on a polarizable surface: Coupling of CO with KTaO3(001)

Polarizable materials attract attention in catalysis because they have a free parameter for tuning chemical reactivity. Their surfaces entangle the dielectric polarization with surface polarity, excess charge, and orbital hybridization. How this affects individual adsorbed molecules is shown for the incipient ferroelectric perovskite KTaO3. This intrinsically polar material cleaves along (001) into KO- and TaO2-terminated surface domains. At TaO2 terraces, the polarity-compensating excess electrons form a two-dimensional electron gas and can also localize by coupling to ferroelectric distortions. TaO2 terraces host two distinct types of CO molecules, adsorbed at equivalent lattice sites but charged differently as seen in atomic force microscopy/scanning tunneling microscopy. Temperature-programmed desorption shows substantially stronger binding of the charged CO; in density functional theory calculations, the excess charge favors a bipolaronic configuration coupled to the CO. These results pinpoint how adsorption states couple to ferroelectric polarization