Surface structure and reactivity of
multi-component oxides at the atomic scale

@article{ANGEWANDTE2022,

title = {Structure of an Ultrathin Oxide on Pt_{3}Sn(111) Solved by Machine Learning Enhanced Global Optimization},

author = {Lindsay R. Merte and Malthe Kjær Bisbo and Igor Sokolović and Martin Setvín and Benjamin Hagman and Mikhail Shipilin and Michael Schmid and Ulrike Diebold and Edvin Lundgren and Bjørk Hammer},

url = {https://onlinelibrary.wiley.com/doi/10.1002/anie.202204244},

doi = {10.1002/anie.202204244},

year  = {2022},

date = {2022-04-05},

urldate = {2022-04-05},

journal = {Angewandte Chemie - International Edition},

abstract = {Determination of the atomic structure of solid surfaces typically depends on comparison of measured properties with simulations based on hypothesized structural models. For simple structures, the models may be guessed, but for more complex structures there is a need for reliable theory-based search algorithms. So far, such methods have been limited by the combinatorial complexity and computational expense of sufficiently accurate energy estimation for surfaces. However, the introduction of machine learning methods has the potential to change this radically. Here, we demonstrate how an evolutionary algorithm, utilizing machine learning for accelerated energy estimation and diverse population generation, can be used to solve an unknown surface structure—the (4×4) surface oxide on Pt_{3}Sn(111)–based on limited experimental input. The algorithm is efficient and robust, and should be broadly applicable in surface studies, where it can replace manual, intuition based model generation.},

keywords = {P02},

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

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{Franceschi2022,

title = {Reconstruction changes drive surface diffusion and determine the flatness of oxide surfaces},

author = {Giada Franceschi and Michael Schmid and Ulrike Diebold and Michele Riva},

doi = {10.1116/6.0001704},

year  = {2022},

date = {2022-02-22},

urldate = {2022-02-22},

journal = {Journal of Vacuum Science & Technology A},

volume = {40},

number = {2},

pages = {023206},

publisher = {American Vacuum Society},

abstract = {Surface diffusion on metal oxides is key in many areas of materials technology, yet it has been scarcely explored at the atomic scale. This work provides phenomenological insights from scanning tunneling microscopy on the link between surface diffusion, surface atomic structure, and oxygen chemical potential based on three model oxide surfaces: Fe}2}O_{3}(1-102), La_{1−x}Sr_{x}MnO_{3}(110), and In_{2}O_{3}(111). In all instances, changing the oxygen chemical potential used for annealing stabilizes reconstructions of different compositions while promoting the flattening of the surface morphology—a sign of enhanced surface diffusion. It is argued that thermodynamics, rather than kinetics, rules surface diffusion under these conditions: the composition change of the surface reconstructions formed at differently oxidizing conditions drives mass transport across the surface.},

keywords = {P02},

pubstate = {published},

tppubtype = {article}

}

@article{ACSMEASURE2022,

title = {Why and How Savitzky–Golay Filters Should Be Replaced},

author = {Michael Schmid and David Rath and Ulrike Diebold},

url = {https://pubs.acs.org/doi/10.1021/acsmeasuresciau.1c00054},

doi = {10.1021/acsmeasuresciau.1c00054},

year  = {2022},

date = {2022-02-17},

urldate = {2022-02-17},

journal = {ACS Measurement Science Au},

volume = {2},

number = {2},

pages = {185--196},

abstract = {Savitzky–Golay (SG) filtering, based on local least-squares fitting of the data by polynomials, is a popular method for smoothing data and calculations of derivatives of noisy data. At frequencies above the cutoff, SG filters have poor noise suppression; this unnecessarily reduces the signal-to-noise ratio, especially when calculating derivatives of the data. In addition, SG filtering near the boundaries of the data range is prone to artifacts, which are especially strong when using SG filters for calculating derivatives of the data. We show how these disadvantages can be avoided while keeping the advantageous properties of SG filters. We present two classes of finite impulse response (FIR) filters with substantially improved frequency response: (i) SG filters with fitting weights in the shape of a window function and (ii) convolution kernels based on the sinc function with a Gaussian-like window function and additional corrections for improving the frequency response in the passband (modified sinc kernel). Compared with standard SG filters, the only price to pay for the improvement is a moderate increase in the kernel size. Smoothing at the boundaries of the data can be improved with a non-FIR method, the Whittaker–Henderson smoother, or by linear extrapolation of the data, followed by convolution with a modified sinc kernel, and we show that the latter is preferable in most cases. We provide computer programs and equations for the smoothing parameters of these smoothers when used as plug-in replacements for SG filters and describe how to choose smoothing parameters to preserve peak heights in spectra.},

keywords = {P02},

pubstate = {published},

tppubtype = {article}

}

@article{Franceschi2021,

title = {Two-dimensional surface phase diagram of a multicomponent perovskite oxide: La_{0.8}Sr_{0.2}MnO_{3} (110)},

author = {Giada Franceschi and Michael Schmid and Ulrike Diebold and Michele Riva},

doi = {10.1103/physrevmaterials.5.l092401},

year  = {2021},

date = {2021-09-24},

urldate = {2021-09-24},

journal = {Physical Review Materials},

volume = {5},

number = {9},

pages = {L092401},

publisher = {American Physical Society (APS)},

abstract = {The many surface reconstructions of (110)-oriented lanthanum strontium manganite (\textbf{La_{0.8}Sr_{0.2}MnO}3}}, LSMO) were followed as a function of the oxygen chemical potential (\textit{\textbf{μ}_{O}}) and the surface cation composition. Decreasing \textit{\textbf{μ}_{O}} causes Mn to migrate across the surface, enforcing phase separation into \textit{\textbf{A}}-site-rich areas and a variety of composition-related, structurally diverse \textit{\textbf{B}}-site-rich reconstructions. The composition of these phase-separated structures was quantified with scanning tunneling microscopy, and these results were used to build a two-dimensional phase diagram of the LSMO(110) equilibrium surface structures.},

keywords = {P02},

pubstate = {published},

tppubtype = {article}

}

@article{Mirabella2021,

title = {Ni-modified Fe_{3}O_{4}(001) surface as a simple model system for understanding the oxygen evolution reaction},

author = {Francesca Mirabella and Matthias Müllner and Thomas Touzalin and Michele Riva and Zdenek Jakub and Florian Kraushofer and Michael Schmid and Marc T M Koper and Gareth S Parkinson and Ulrike Diebold},

doi = {10.1016/j.electacta.2021.138638},

year  = {2021},

date = {2021-09-01},

urldate = {2021-09-01},

journal = {Electrochimica Acta},

volume = {389},

pages = {138638},

publisher = {Elsevier BV},

abstract = {Electrochemical water splitting is an environmentally friendly technology to store renewable energy in the form of chemical fuels. Among the earth-abundant first-row transition metal-based catalysts, mixed Ni-Fe oxides have shown promising performance for effective and low-cost catalysis of the oxygen evolution reaction (OER) in alkaline media, but the synergistic roles of Fe and Ni cations in the OER mechanism remain unclear. In this work, we report how addition of Ni changes the reactivity of a model iron oxide catalyst, based on Ni deposited on and incorporated in a magnetite Fe_{3}O_{4}(001) single crystal, using a combination of surface science techniques in ultra-high vacuum such as low energy electron diffraction (LEED), x-ray photoelectron spectroscopy (XPS), low-energy ion scattering (LEIS), and scanning tunneling microscopy (STM), as well as atomic force microscopy (AFM) in air, and electrochemical methods such as cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) in alkaline media. A significant improvement in the OER activity is observed when the top surface presents an iron fraction among the cations in the range of 20-40%, which is in good agreement with what has been observed for powder catalysts. Furthermore, a decrease in the OER overpotential is observed following surface aging in electrolyte for three days. At higher Ni load, AFM shows the growth of a new phase attributed to an (oxy)-hydroxide phase which, according to CV measurements, does not seem to correlate with the surface activity towards OER. EIS suggests that the OER precursor species observed on the clean and Ni-modified surfaces are similar and Fe-centered, but form at lower overpotentials when the surface Fe:Ni ratio is optimized. We propose that the well-defined Fe_{3}O_{4}(001) surface can serve as a model system for understanding the OER mechanism and establishing the structure-reactivity relation on mixed Fe-Ni oxides.},

keywords = {P02, P04, pre-TACO},

pubstate = {published},

tppubtype = {article}

}

@article{Franchini2021,

title = {Polarons in materials},

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

doi = {10.1038/s41578-021-00289-w},

year  = {2021},

date = {2021-03-19},

journal = {Nature Reviews Materials},

publisher = {Springer Science and Business Media LLC},

abstract = {Polarons are quasiparticles that easily form in polarizable materials due to the coupling of excess electrons or holes with ionic vibrations. These quasiparticles manifest themselves in many different ways and have a profound impact on materials properties and functionalities. Polarons have been the testing ground for the development of numerous theories, and their manifestations have been studied by many different experimental probes. This Review provides a map of the enormous amount of data and knowledge accumulated on polaron effects in materials, ranging from early studies and standard treatments to emerging experimental techniques and novel theoretical and computational approaches.},

keywords = {P02, P07, pre-TACO},

pubstate = {published},

tppubtype = {article}

}

@article{Hulva2021,

title = {Unraveling CO adsorption on model single-atom catalysts},

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

doi = {10.1126/science.abe5757},

year  = {2021},

date = {2021-01-22},

journal = {Science},

volume = {371},

number = {6527},

pages = {375--379},

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

abstract = {Understanding how the local environment of a “single-atom” catalyst affects stability and reactivity remains a challenge. We present an in-depth study of copper_{1}, silver_{1}, gold_{1}, nickel_{1}, palladium_{1}, platinum_{1}, rhodium_{1}, and iridium_{1} species on Fe_{3}O_{4}(001), a model support in which all metals occupy the same twofold-coordinated adsorption site upon deposition at room temperature. Surface science techniques revealed that CO adsorption strength at single metal sites differs from the respective metal surfaces and supported clusters. Charge transfer into the support modifies the d-states of the metal atom and the strength of the metal–CO bond. These effects could strengthen the bond (as for Ag_{1}–CO) or weaken it (as for Ni_{1}–CO), but CO-induced structural distortions reduce adsorption energies from those expected on the basis of electronic structure alone. The extent of the relaxations depends on the local geometry and could be predicted by analogy to coordination chemistry.},

keywords = {P02, P04, P07, pre-TACO},

pubstate = {published},

tppubtype = {article}

}

@article{Timmermann2020,

title = {IrO_{2} Surface Complexions Identified through Machine Learning and Surface Investigations},

author = {Jakob Timmermann and Florian Kraushofer and Nikolaus Resch and Peigang Li and Yu Wang and Zhiqiang Mao and Michele Riva and Yonghyuk Lee and Carsten Staacke and Michael Schmid and Christoph Scheurer and Gareth S Parkinson and Ulrike Diebold and Karsten Reuter},

doi = {10.1103/physrevlett.125.206101},

year  = {2020},

date = {2020-11-10},

journal = {Physical Review Letters},

volume = {125},

number = {20},

pages = {206101},

publisher = {American Physical Society (APS)},

abstract = {A Gaussian approximation potential was trained using density-functional theory data to enable a global geometry optimization of low-index rutile IrO_{2} facets through simulated annealing. Ab initio thermodynamics identifies (101) and (111) (1×1) terminations competitive with (110) in reducing environments. Experiments on single crystals find that (101) facets dominate and exhibit the theoretically predicted (1×1) periodicity and x-ray photoelectron spectroscopy core-level shifts. The obtained structures are analogous to the complexions discussed in the context of ceramic battery materials.},

keywords = {P02, P04, pre-TACO},

pubstate = {published},

tppubtype = {article}

}