Polaron pattern recognition
in correlated oxide surfaces

@article{Ranalli2023,

title = {Temperature-dependent anharmonic phonons in quantum paraelectric KTaO_{3} by first principles and machine-learned force fields},

author = {Luigi Ranalli and Carla Verdi and Lorenzo Monacelli and Georg Kresse and Matteo Calandra and Cesare Franchini},

doi = {10.1002/qute.202200131},

year  = {2023},

date = {2023-02-22},

urldate = {2023-02-22},

journal = {Advanced Quantum Technology},

volume = {6},

issue = {4},

abstract = {Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anharmonic phonon effects. The temperature-dependent evolution of the soft ferroelectric phonon mode in the quantum paraelectric KTaO_{3} in the range 0–300 K is modeled by combining density functional theory (DFT) calculations with the stochastic self-consistent harmonic approximation assisted by an on-the-fly machine-learned force field. The calculated data show that including anharmonic terms is essential to stabilize the spurious imaginary ferroelectric phonon predicted by DFT in the harmonic approximation, in agreement with experiments. Augmenting the DFT workflow with machine-learned force fields allows for efficient stochastic sampling of the configuration space using large supercells in a wide temperature range, inaccessible to conventional ab initio protocols. This work proposes a robust computational workflow capable of accounting for collective behaviors involving different degrees of freedom and occurring at large time/length scales, paving the way for precise modeling and control of quantum effects in materials.},

keywords = {P03, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Corrias2023,

title = {Automated Real-Space Lattice Extraction for Atomic Force Microscopy Images},

author = {Marco Corrias and Lorenzo Papa and Igor Sokolovíc and Viktor Birschitzky and Alexander Gorfer and Martin Setvin and Michael Schmid and Ulrike Diebold and Michele Reticcioli and Cesare Franchini},

doi = {10.1088/2632-2153/acb5e0},

year  = {2023},

date = {2023-01-24},

urldate = {2023-01-24},

journal = {Machine Learning: Science and Technology},

volume = {4},

pages = {015015},

abstract = {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 TiO_{2}(101), oxygen deficient rutile TiO_{2}(110) with and without CO adsorbates, SrTiO_{3}(001) with Sr vacancies and graphene with C vacancies. The code delivers excellent results and is tested against atom misclassification and artifacts, thereby facilitating the interpretation of scanning probe microscopy images.},

keywords = {P02, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Wang2022,

title = {Surface chemistry on a polarizable surface: Coupling of CO with KTaO _{3}(001)},

author = {Zhichang Wang and Michele Reticcioli and Zdenek Jakub and Igor Sokolović and Matthias Meier and Lynn A Boatner and Michael Schmid and Gareth S. Parkinson and Ulrike Diebold and Cesare Franchini and Martin Setvin},

url = {https://www.science.org/doi/10.1126/sciadv.abq1433},

doi = {10.1126/sciadv.abq1433},

year  = {2022},

date = {2022-08-19},

urldate = {2022-08-19},

journal = {Science Advances},

volume = {8},

issue = {33},

publisher = {American Association for the Advancement of Science (AAAS)},

abstract = {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 KTaO_{3}. This intrinsically polar material cleaves along (001) into KO- and TaO_{2}-terminated surface domains. At TaO_{2} terraces, the polarity-compensating excess electrons form a two-dimensional electron gas and can also localize by coupling to ferroelectric distortions. TaO_{2} 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.},

keywords = {P02, P04, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Sombut2022,

title = {Role of Polarons in Single-Atom Catalysts: Case Study of Me_{1}[Au_{1},Pt_{1} and Rh_{1}] on TiO_{2}(110)},

author = {Panukorn Sombut and Lena Puntscher and Marlene Atzmüller and Zdenek Jakub and Michele Reticcioli and Matthias Meier and Gareth S. Parkinson and Cesare Franchini},

doi = {10.1007/s11244-022-01651-0},

year  = {2022},

date = {2022-07-25},

journal = {Topics in Catalysis},

volume = {65},

pages = {1620--1630},

abstract = {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 TiO_{2}(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.},

keywords = {P04, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Reticcioli2022,

title = {Competing electronic states emerging on polar surfaces},

author = {Michele Reticcioli and Zhichang Wang and Michael Schmid and Dominik Wrana and Lynn A. Boatner and Ulrike Diebold and Martin Setvin and Cesare Franchini},

url = {https://www.nature.com/articles/s41467-022-31953-6},

doi = {10.1038/s41467-022-31953-6},

year  = {2022},

date = {2022-07-25},

urldate = {2022-07-25},

journal = {Nature Communications},

volume = {13},

number = {4311},

publisher = {Springer Science and Business Media LLC},

abstract = {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.},

keywords = {P02, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Birschitzky2022,

title = {Machine learning for exploring small polaron configurational space},

author = {Viktor C Birschitzky and Florian Ellinger and Ulrike Diebold and Michele Reticcioli and Cesare Franchini},

url = {https://www.nature.com/articles/s41524-022-00805-8},

doi = {10.1038/s41524-022-00805-8},

year  = {2022},

date = {2022-06-06},

urldate = {2022-06-06},

journal = {npj Computational Materials},

volume = {8},

number = {125},

publisher = {Springer Science and Business Media LLC},

abstract = {Polaron defects are ubiquitous in materials and play an important role in many processes involving carrier mobility, charge transfer and surface reactivity. Determining small polarons’ spatial distributions is essential to understand materials properties and functionalities. However, the required exploration of the configurational space is computationally demanding when using first principles methods. Here, we propose a machine-learning (ML) accelerated search that determines the ground state polaronic configuration. The ML model is trained on databases of polaron configurations generated by density functional theory (DFT) via molecular dynamics or random sampling. To establish a mapping between configurations and their stability, we designed descriptors modelling the interactions among polarons and charged point defects. We used the DFT+ML protocol to explore the polaron configurational space for two surface-systems, reduced rutile TiO_{2}(110) and Nb-doped SrTiO_{3}(001). The ML-aided search proposes additional polaronic configurations and can be utilized to determine optimal polaron distributions at any charge concentration.},

keywords = {P02, P07},

pubstate = {published},

tppubtype = {article}

}

@article{SCIADV2022,

title = {CO oxidation by Pt_{2}/Fe_{3}O_{4}: Metastable dimer and support configurations facilitate lattice oxygen extraction},

author = {Matthias Meier and Jan Hulva and Zdenek Jakub and Florian Kraushofer and Mislav Bobić and Roland Bliem and Martin Setvin and Michael Schmid and Ulrike Diebold and Cesare Franchini and Gareth S. Parkinson},

url = {https://www.science.org/doi/10.1126/sciadv.abn4580},

doi = {10.1126/sciadv.abn4580},

year  = {2022},

date = {2022-04-01},

urldate = {2022-04-01},

journal = {Science Advances},

volume = {8},

issue = {13},

pages = {eabn4580},

abstract = {Heterogeneous catalysts based on subnanometer metal clusters often exhibit strongly size-dependent properties, and the addition or removal of a single atom can make all the difference. Identifying the most active species and deciphering the reaction mechanism is extremely difficult, however, because it is often not clear how the catalyst evolves in operando. Here, we use a combination of atomically resolved scanning probe microscopies, spectroscopic techniques, and density functional theory (DFT)–based calculations to study CO oxidation by a model Pt/Fe_{3}O_{4}(001) “single-atom” catalyst. We demonstrate that (PtCO)_{2} dimers, formed dynamically through the agglomeration of mobile Pt-carbonyl species, catalyze a reaction involving the oxide support to form CO_{2}. Pt_{2} dimers produce one CO_{2} molecule before falling apart into two adatoms, releasing the second CO. Olattice extraction only becomes facile when both the Pt-dimer and the Fe_{3}O_{4} support can access metastable configurations, suggesting that substantial, concerted rearrangements of both cluster and support must be considered for reactions occurring at elevated temperature.},

keywords = {P02, P04, P07},

pubstate = {published},

tppubtype = {article}

}

@article{JPCM2022,

title = {Modeling polarons in density functional theory: lessons learned from TiO_{2}},

author = {Michele Reticcioli and Ulrike Diebold and Cesare Franchini},

url = {https://iopscience.iop.org/article/10.1088/1361-648X/ac58d7},

doi = {10.1088/1361-648X/ac58d7},

year  = {2022},

date = {2022-03-14},

urldate = {2022-03-14},

journal = {Journal of Physics: Condensed Matter},

volume = {34},

number = {20},

pages = {204006},

abstract = {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 TiO_{2}, 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.},

keywords = {P02, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Jakub2021,

title = {Rapid oxygen exchange between hematite and water vapor},

author = {Zdenek Jakub and Matthias Meier and Florian Kraushofer and Jan Balajka and Jiri Pavelec and Michael Schmid and Cesare Franchini and Ulrike Diebold and Gareth S. Parkinson},

doi = {10.1038/s41467-021-26601-4},

year  = {2021},

date = {2021-11-10},

journal = {Nature Communications},

volume = {12},

number = {6488},

issue = {1},

publisher = {Springer Science and Business Media LLC},

abstract = {Oxygen exchange at oxide/liquid and oxide/gas interfaces is important in technology and environmental studies, as it is closely linked to both catalytic activity and material degradation. The atomic-scale details are mostly unknown, however, and are often ascribed to poorly defined defects in the crystal lattice. Here we show that even thermodynamically stable, well-ordered surfaces can be surprisingly reactive. Specifically, we show that all the 3-fold coordinated lattice oxygen atoms on a defect-free single-crystalline “r-cut” (1-102) surface of hematite (α-Fe_{2}O_{3}) are exchanged with oxygen from surrounding water vapor within minutes at temperatures below 70 °C, while the atomic-scale surface structure is unperturbed by the process. A similar behavior is observed after liquid-water exposure, but the experimental data clearly show most of the exchange happens during desorption of the final monolayer, not during immersion. Density functional theory computations show that the exchange can happen during on-surface diffusion, where the cost of the lattice oxygen extraction is compensated by the stability of an HO-HOH-OH complex. Such insights into lattice oxygen stability are highly relevant for many research fields ranging from catalysis and hydrogen production to geochemistry and paleoclimatology.},

keywords = {P02, P04, P07},

pubstate = {published},

tppubtype = {article}

}