Search

Polarons are crucial for charge transport in semiconductors, significantly impacting material properties and device performance. The dynamics of small polarons can be investigated using first-principles molecular dynamics. However, the limited timescale of these simulations presents a challenge for adequately sampling infrequent polaron hopping events. Here, we introduce a message-passing neural network that learns the polaronic potential energy surface by encoding the polaronic state, allowing for simulations of polaron hopping dynamics at the nanosecond scale. By leveraging the statistical significance of the long timescale, our framework can accurately estimate polaron (anisotropic) mobilities and activation barriers in prototypical polaronic oxides across different scenarios (hole polarons in rocksalt MgO and electron polarons in pristine and F-doped rutile TiO$_2$) in excellent agreement with experimental observations.

The multifaceted physics of oxides is shaped by their composition and the presence of defects, which are often accompanied by the formation of polarons. The simultaneous presence of polarons and defects, and their complex interactions, pose challenges for first-principles simulations and experimental techniques. In this study, we leverage machine learning and a first-principles database to analyze the distribution of surface oxygen vacancies (V$_{\rm O}$) and induced small polarons on rutile TiO$_2$(110), effectively disentangling the interactions between polarons and defects. By combining neural-network supervised learning and simulated annealing, we elucidate the inhomogeneous V$_{\rm O}$ distribution observed in scanning probe microscopy (SPM). Our innovative approach allows us to understand and predict defective surface patterns at previously inaccessible length scales, identifying the specific role of individual types of defects. Specifically, surface-polaron-stabilizing V$_{\rm O}$-configurations are identified, which could have consequences for surface reactivity.

The discovery of ambient superconductivity would mark an epochal breakthrough long-awaited for over a century, potentially ushering in unprecedented scientific and technological advancements. The recent findings on high-temperature superconducting phases in various hydrides under high pressure have ignited optimism, suggesting that the realization of near-ambient superconductivity might be on the horizon. However, the preparation of hydride samples tends to promote the emergence of various metastable phases, marked by a low level of experimental reproducibility. Identifying these phases through theoretical and computational methods entails formidable challenges, often resulting in controversial outcomes. In this paper, we consider N-doped LuH3 as a prototypical complex hydride: By means of machine-learning-accelerated force-field molecular dynamics, we have identified the formation of H2 molecules stabilized at ambient pressure by nitrogen impurities. Importantly, we demonstrate that this molecular phase plays a pivotal role in the emergence of a dynamically stable, low-temperature, experimental-ambient-pressure superconductivity. The potential to stabilize hydrogen in molecular form through chemical doping opens up a novel avenue for investigating disordered phases in hydrides and their transport properties under near-ambient conditions.

In polarizable materials, electronic charge carriers interact with the surrounding ions, leading to quasiparticle behaviour. The resulting polarons play a central role in many materials properties including electrical transport, optical properties, surface reactivity and magnetoresistance, and polaron properties are typically investigated indirectly through such macroscopic characteristics. Here, noncontact atomic force microscopy (nc-AFM) is used to directly image polarons in Fe2O3 at the single quasiparticle limit. A combination of Kelvin probe force microscopy (KPFM) and kinetic Monte Carlo (KMC) simulations shows that Ti doping dramatically enhances the mobility of electron polarons, and density functional theory (DFT) calculations indicate that a metallic transition state is responsible for the enhancement. In contrast, hole polarons are significantly less mobile and their hopping is hampered further by the introduction of trapping centres.

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.

The cubic perovsike strontium titanate SrTiO$_3$ (STO) is one of the most studied, polarizable transition metal oxides. When excess charge is introduced to this material, e.g., through doping or atomic defects, STO tends to host polarons: Quasi-particles formed by excess charge carriers coupling with the crystal phonon field. Their presence alters the materials properties, and is a key for many applications. Considering that polarons form preferentially on or near surfaces, we study small polaron formation at the TiO$_2$ termination of the STO(001) surface via density functional theory calculations. We model several supercell slabs of Nb-doped and undoped STO(001) surfaces with increasing size, also considering the recently observed as-cleaved TiO$_2$ terminated surface hosting Sr-adatoms. Our findings suggest that small polarons become less stable at low concentrations of Nb-doping, in analogy with polarons localized in the bulk. Further, we inspect the stability of different polaron configurations with respect to Nb- and Sr-impurities, and discuss their spectroscopic properties., Comment: 8 pages, 7 figures, 1 table

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 TaO$_2$ termination of KTaO$_3$(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. Controlling the degree of charge ordering and the transition from ferroelectric to paraelectric states could be of great benefit for the generation and transport of carriers in electronic applications.

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 Rh1, Pt1, and Au1 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 Pt1 and Au1 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 Rh1 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.

Polaron defects are ubiquitous in materials and play an important role in many processes involving carrier mobility, charge transfer and surface reactivity. Determining the spatial distribution of small polarons is essential to understand materials properties and functionalities. This requires an exploration of the configurational space, which is computationally demanding when using standard first principles methods, and technically prohibitive for many-polaron systems. Here, we propose a machine-learning (ML) accelerated search that compares the energy stability of different polaron patterns and determines the ground state configuration. The kernel-regression based ML model is trained on databases generated by density functional theory (DFT) calculations on a minimal set of initial polaron patterns, obtained by using either molecular dynamics simulations or a random sampling approach. To establish an efficient mapping between training data and configuration stability we designed simple descriptors that model the interactions among polarons and charged point defects. The proposed DFT+ML protocol is used here to explore millions of polaron configurations for two different systems, oxygen defective rutile TiO$_2$(110) and Nb-doped SrTiO$_3$(001). Our data shows that the ML-aided search correctly individuates the ground-state polaron patterns, proposes polaronic configurations not visited in the training and can be used to efficiently determine the optimal distribution of polarons at any charge concentration.

Modern computing facilities grant access to first-principles density-functional theory study of complex physical and chemical phenomena in materials, that require large supercell to properly model the system. However, supercells are associated to small Brillouin zones in the reciprocal space, leading to folded electronic eigenstates that make the analysis and interpretation extremely challenging. Various techniques have been proposed and developed in order to reconstruct the electronic band structures of super cells, unfolded into the reciprocal space of an ideal primitive cell. Here, we propose an efficient unfolding scheme embedded directly in the Vienna Ab-initio Simulation Package (VASP), that requires modest computational resources and allows for an automatized mapping from the reciprocal space of the supercell to primitive cell Brillouin zone. This algorithm can computes band structures, Fermi surfaces and spectral functions, by using an integrated post-processing tool (bands4vasp). The method is here applied to a selected variety of complex physical situations: the effect of doping on the band dispersion in the BaFe$_{\rm 2(1-x)}$Ru$_{\rm 2x}$As$_2$ superconductor, the interaction between adsorbates and polaronic states on the TiO$_2$(110) surface, and the band splitting induced by non-collinear spin fluctuations in EuCd$_2$As$_2$.

Growth of two-dimensional metals has eluded materials scientists since the discovery of the atomically thin graphene and other covalently bound 2D materials. Here, we report a two-atom-thick hexagonal copper-gold alloy, grown through thermal evaporation on freestanding graphene and hexagonal boron nitride. The structures are imaged at atomic resolution with scanning transmission electron microscopy and further characterized with spectroscopic techniques. Electron irradiation in the microscope provides sufficient energy for a phase transformation of the 2D structure--atoms are released from their lattice sites with the gold atoms eventually forming face-centered cubic nanoclusters on top of 2D regions during observation. The presence of copper in the alloy enhances sticking of gold to the substrate, which has clear implications for creating atomically thin electrodes for applications utilizing 2D materials. Its practically infinite surface-to-bulk ratio also makes the 2D CuAu particularly interesting for catalysis applications., Comment: 24 pages, including 4 figures in main text and 6 in the supplementary information

The formation of polarons is a pervasive phenomenon in transition metal oxide compounds, with a strong impact on the physical properties and functionalities of the hosting materials. In its original formulation the polaron problem considers a single charge carrier in a polar crystal interacting with its surrounding lattice. Depending on the spatial extension of the polaron quasiparticle, originating from the coupling between the excess charge and the phonon field, one speaks of small or large polarons. This chapter discusses the modeling of small polarons in real materials, with a particular focus on the archetypal polaron material TiO2. After an introductory part, surveying the fundamental theoretical and experimental aspects of the physics of polarons, the chapter examines how to model small polarons using first principles schemes in order to predict, understand and interpret a variety of polaron properties in bulk phases and surfaces. Following the spirit of this handbook, different types of computational procedures and prescriptions are presented with specific instructions on the setup required to model polaron effects., Comment: 36 pages, 12 figures

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 realistically within a supercell approach by replacing sodium with magnesium at different concentrations. Our data indicate that by increasing carrier concentration the system is subjected to two types of transition: (i) insulator to bad metal at low doping and low temperature and (ii) bad metal to metal at high doping and/or high-temperature. The predicted doping-induced insulator to metal transition (MIT) has similar traits with the temperature driven MIT reported in the undoped compound. Both develops in an itinerant background and exhibit a coupled electronic and magnetic behaviour characterized by the gradual quenching of the (pseudo)-gap associated with an reduction of the local spin moment. Unlike the temperature-driven MIT, chemical doping induces substantial modifications of the band structure and the MIT cannot be fully described as a Lifshitz process., Comment: 5 pages, 4 figures

Activator impurities and their distribution in the host lattice play a key role in scintillation phenomena. Here a combination of cross-sectional noncontact atomic force microscopy (nc-AFM), X-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) was used to study the distribution of Eu2+ dopants in a NaI scintillator activated by 3% of EuI2. Two types of precipitate structures were found. First, needle-shaped EuI2 precipitates with a layered structure are likely responsible for scattering the scintillation light. In transparent crystals with good scintillation properties, precipitates with a cubic crystal structure and a size below 4 nm were found. A surprisingly low concentration of point defects was detected in all of the investigated samples. Upon annealing, Eu segregates towards the surface, which results in the formation of an ordered hexagonal overlayer with the EuI2 composition and a pronounced, unidirectional moire pattern.

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.

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.

We study the effects of dilute La and Rh doping on the electronic structure of the relativistic Mott insulator Sr$_2$IrO$_4$ using fully relativistic and magnetically non-collinear density functional theory with the inclusion of an on-site Hubbard $U$. To model doping effects, we have adopted the supercell approach that allows for a realistic treatment of structural relaxations and electronic effects beyond a purely rigid band approach. By means of the band unfolding technique we have computed the spectral function and constructed the effective band structure and Fermi surface (FS), which are readily comparable with available experimental data. Our calculations clearly indicate that La and Rh doping can be interpreted as effective electron and (fractional) hole doping, respectively. In Sr$_{2-x}$La$_x$IrO$_4$ the IMT is accompanied by a moderate renormalization of the electronic correlation substantiated by a reduction of the effective on-site Coulomb repulsion $U-J$ from 1.6 eV ($x=0$) to 1.4 eV (metallic regime $x=12.5\%$). The progressive closing of the relativistic Mott gap leads to the emergence of connected elliptical electron pockets at ($\pi$/2,$\pi$/2) and less intense features at $X$ in the FS. The substitution of Ir with the nominally isovalent Rh is accompanied by a substantial hole transfer from the Rh site to the nearest neighbor Ir sites. This shifts down the chemical potential, creates almost circular disconnected hole pockets in the FS and establishes the emergence of a two-dimensional metallic state formed by conducting Rh-planes intercalated by insulating Ir-planes. Finally, our data indicate that hole doping causes a flipping of the in-plane net ferromagnetic moment in the Rh plane and induces a magnetic transition from the AF-I to the AF-II ordering.