Bayesian regression for
multi-level machine-learned potentials
Subproject P03
The first-principles description of the properties of multi-component metal oxides is an exceedingly challenging problem. The reasons are that the configurational space grows exponentially with the number of species and standard Density Functional Theory (DFT) is often not accurate enough. The long-term objective of P03 is to accelerate first-principles calculations by developing machine-learning approaches for the description of the interatomic forces, Born effective charges, and other tensorial properties of multivalent oxides. The project will rely on kernel-based methods and Bayesian inference to implement fully automatic “on-the-fly” learning.
In the first project period, we will develop machine-learned force fields (MLFF) for DFT and DFT+U, whereby the number of components in the FF will be gradually increased. A concise framework for learning tensorial properties will be implemented. We will use this to simulate infrared spectra of oxide materials, which can be readily compared to the finite-temperature spectra measured by the experimental groups.
The difference between DFT and hybrid functionals will be machine-learned to go beyond semi-local functionals (Delta-learning). The long-term perspective is to extend this approach to highly accurate beyond-DFT methods, such as the random phase approximation and quantum chemistry (coupled cluster) methods. Although kernel-based methods are exceedingly accurate, they are often less efficient than NN. We will collaborate with other projects to recast the on-the-fly trained FF into NN potentials to address this issue.
Expertise
The main research efforts of the group are directed towards the development of quantum-mechanical tools for atomic-scale simulations of properties and processes in materials and the application of these methodologies to key areas of condensed matter physics and materials research. An important pillar of the research is the Vienna Ab initio Simulation Package (VASP), a general-purpose ab initio code for solving the many-electron Schrödinger equation. The code is among the world leaders in its field, with more than 3500 licensees worldwide. We have expertise with simulations for a vast number of properties using many different techniques:
- Density functional theory (DFT), including spin and non-collinear DFT
- Linear response theory to calculate phonons and dielectric properties
- Hartree-Fock techniques and many flavors of hybrid functionals
- Many-body perturbation theory, including GW and Bethe-Salpeter
- Wavefunction-based correlated methods (Møller-Plesset perturbation theory)
- Surface science, including growth and oxide formation
- Simulation of nanostructures
- Semiconductor physics: charge trapping, polarons
- Electronic excitations
- Defect energies in extended systems
For TACO, we will adapt our machine-learning techniques to tensorial properties and correlated wavefunction techniques. These techniques are directly integrated into VASP and allow to accelerate finite-temperature simulations by many orders of magnitudes.
Team
Former Members
Publications
2024
Rath, David; Mikerásek, Vojtěch; Wang, Chunlei; Eder, Moritz; Schmid, Michael; Diebold, Ulrike; Parkinson, Gareth S.; Pavelec, Jiří
Journal ArticleOpen AccessSubmittedarXivIn: arXiv, 2024.
Abstract | Links | BibTeX | Tags: P02, P04
@article{Rath2024,
title = {Infrared Reflection Absorption Spectroscopy Setup with Incidence Angle Selection for Surfaces of Non-Metals},
author = {David Rath and Vojtěch Mikerásek and Chunlei Wang and Moritz Eder and Michael Schmid and Ulrike Diebold and Gareth S. Parkinson and Jiří Pavelec},
url = {https://arxiv.org/abs/2403.19263},
doi = {10.48550/ARXIV.2403.19263},
year = {2024},
date = {2024-03-28},
journal = {arXiv},
publisher = {arXiv},
abstract = {Infrared Reflection Absorption Spectroscopy (IRAS) on dielectric single crystals is challenging because the optimal incidence angles for light-adsorbate interaction coincide with regions of low IR reflectivity. Here, we introduce an optimized IRAS setup that maximizes the signal-to-noise ratio for non-metals. This is achieved by maximizing light throughput, and by selecting optimal incidence angles that directly impact the peak heights in the spectra. The setup uses a commercial FTIR spectrometer and is usable in ultra-high vacuum (UHV). Specifically, the design features sample illumination and collection mirrors with a high numerical aperture inside the UHV system, and an adjustable aperture to select the incidence angle range on the sample. This is important for p-polarized measurements on dielectrics, because the peaks in the spectra reverse direction at the Brewster angle (band inversion). The system components are connected precisely via a single flange, ensuring long-term stability. We studied the signal-to-noise (SNR) variation in p-polarized IRAS spectra for one monolayer of CO on TiO_{2}(110) as a function of incidence angle range, where a maximum signal-to-noise ratio of 70 was achieved at 4 cm^{-1} resolution in five minutes measurement time. The capabilities for s-polarization are demonstrated by measuring one monolayer D_{2}O adsorbed on a TiO_{2}(110) surface, where a SNR of 65 was achieved at a delta_R/R0 peak height of 1.4x10-4 in twenty minutes.},
keywords = {P02, P04},
pubstate = {published},
tppubtype = {article}
}
Wang, Chunlei; Sombut, Panukorn; Puntscher, Lena; Jakub, Zdenek; Meier, Matthias; Pavelec, Jiri; Bliem, Roland; Schmid, Michael; Diebold, Ulrike; Franchini, Cesare; Parkinson, Gareth S.
CO‐Induced Dimer Decay Responsible for Gem‐Dicarbonyl Formation on a Model Single‐Atom Catalyst
Journal ArticleOpen AccessIn PressIn: Angewandte Chemie - International Edition, no. e202317347, 2024, ISSN: 1521-3773.
Abstract | Links | BibTeX | Tags: P02, P04, P07
@article{Wang2024,
title = {CO‐Induced Dimer Decay Responsible for Gem‐Dicarbonyl Formation on a Model Single‐Atom Catalyst},
author = {Chunlei Wang and Panukorn Sombut and Lena Puntscher and Zdenek Jakub and Matthias Meier and Jiri Pavelec and Roland Bliem and Michael Schmid and Ulrike Diebold and Cesare Franchini and Gareth S. Parkinson},
doi = {10.1002/anie.202317347},
issn = {1521-3773},
year = {2024},
date = {2024-01-31},
journal = {Angewandte Chemie - International Edition},
number = {e202317347},
publisher = {Wiley},
abstract = {The ability to coordinate multiple reactants at the same active site is important for the wide-spread applicability of single-atom catalysis. Model catalysts are ideal to investigate the link between active site geometry and reactant binding, because the structure of single-crystal surfaces can be precisely determined, the adsorbates imaged by scanning tunneling microscopy (STM), and direct comparisons made to density functional theory. In this study, we follow the evolution of Rh_{1} adatoms and minority Rh_{2} dimers on Fe_{3}O_{4}(001) during exposure to CO using time-lapse STM at room temperature. CO adsorption at Rh_{1} sites results exclusively in stable Rh_{1}CO monocarbonyls, because the Rh atom adapts its coordination to create a stable pseudo-square planar environment. Rh_{1}(CO)_{2} gem-dicarbonyl species are also observed, but these form exclusively through the breakup of Rh_{2} dimers via an unstable Rh_{2}(CO)_{3} intermediate. Overall, our results illustrate how minority species invisible to area-averaging spectra can play an important role in catalytic systems, and show that the decomposition of dimers or small clusters can be an avenue to produce reactive, metastable configurations in single-atom catalysis.},
keywords = {P02, P04, P07},
pubstate = {published},
tppubtype = {article}
}
Birschitzky, Viktor; Sokolovic, Igor; Prezzi, Michael; Palotas, Krisztian; Setvin, Martin; Diebold, Ulrike; Reticcioli, Michele; Franchini, Cesare
Machine Learning Based Prediction of Polaron-Vacancy Patterns on the TiO2(110) Surface
Journal ArticleOpen AccessAccepted ArticlearXivIn: npj Computational Materials, 2024.
Abstract | Links | BibTeX | Tags: P02, P07
@article{Birschitzky_2024a,
title = {Machine Learning Based Prediction of Polaron-Vacancy Patterns on the TiO2(110) Surface},
author = {Viktor Birschitzky and Igor Sokolovic and Michael Prezzi and Krisztian Palotas and Martin Setvin and Ulrike Diebold and Michele Reticcioli and Cesare Franchini},
doi = { https://doi.org/10.48550/arXiv.2401.12042},
year = {2024},
date = {2024-01-22},
journal = {npj Computational Materials},
abstract = {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 (VO) and induced small polarons on rutile TiO2(110), effectively disentangling the interactions between polarons and defects. By combining neural-network supervised learning and simulated annealing, we elucidate the inhomogeneous VO 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 VO-configurations are identified, which could have consequences for surface reactivity.},
keywords = {P02, P07},
pubstate = {published},
tppubtype = {article}
}
2023
Gamba, Oscar; Eder, Moritz; Poglitsch, Matthias; Pavelec, Jiri; Sombut, Panukorn; Meier, Matthias; Diebold, Ulrike; Schmid, Michael; Parkinson, Gareth S.
Formation and stability of Fe-rich terminations of the Fe3O4(001) surface
Journal ArticleOpen AccessIn: Materials Research Express, vol. 10, iss. 44, no. 116517, 2023.
Abstract | Links | BibTeX | Tags: P02, P04
@article{Gamba2023,
title = {Formation and stability of Fe-rich terminations of the Fe_{3}O_{4}(001) surface},
author = {Oscar Gamba and Moritz Eder and Matthias Poglitsch and Jiri Pavelec and Panukorn Sombut and Matthias Meier and Ulrike Diebold and Michael Schmid and Gareth S. Parkinson},
doi = {10.1088/2053-1591/ad0ac5},
year = {2023},
date = {2023-11-22},
urldate = {2023-11-22},
journal = {Materials Research Express},
volume = {10},
number = {116517},
issue = {44},
publisher = {IOP Publishing},
abstract = {Understanding how the structure of iron oxide surfaces varies with their environment is essential for rationalizing their role in (geo-)chemistry and optimizing their application in modern technologies. In this paper, we create Fe-rich terminations of Fe_{3}O_{4}(001) by depositing iron directly onto the 'subsurface cation vacancy'-reconstructed surface, which is the most stable surface under ultrahigh vacuum conditions. Scanning tunneling microscopy and x-ray photoelectron spectroscopy data reveal that the excess iron is initially accommodated as two-fold coordinated adatoms and later incorporates into the subsurface cation vacancies. As the coverage increases, small patches of the octahedral pair termination (also known as the 'Fe dimer' termination) nucleate, eventually covering the entire surface after the deposition of 2 iron atoms per (√2×√2)R45° unit cell. This conclusion effectively rules out some existing models for the termination and provides support for the model proposed by Rustad \textit{et al.} (Surface Science 432, L583-L588, 1999), highlighting the need for further theoretical work to complete the Fe_{3}O_{4}(001) surface phase diagram. The octahedral pair termination is found to be unstable above 523 K and upon exposure to molecular O2 because the excess iron atoms agglomerate to form small FeO_{x} clusters.},
keywords = {P02, P04},
pubstate = {published},
tppubtype = {article}
}
Puntscher, Lena; Sombut, Panukorn; Wang, Chunlei; Ulreich, Manuel; Pavelec, Jiri; Rafsanjani-Abbasi, Ali; Meier, Matthias; Lagin, Adam; Setvin, Martin; Diebold, Ulrike; Franchini, Cesare; Schmid, Michael; Parkinson, Gareth S.
A Multitechnique Study of C2H4 Adsorption on Fe3O4(001)
Journal ArticleOpen AccessIn: Journal of Physical Chemistry C, vol. 127, iss. 37, pp. 18378–18388, 2023.
Abstract | Links | BibTeX | Tags: P02, P04, P07
@article{Puntscher2023,
title = {A Multitechnique Study of C_{2}H_{4} Adsorption on Fe_{3}O_{4}(001)},
author = {Lena Puntscher and Panukorn Sombut and Chunlei Wang and Manuel Ulreich and Jiri Pavelec and Ali Rafsanjani-Abbasi and Matthias Meier and Adam Lagin and Martin Setvin and Ulrike Diebold and Cesare Franchini and Michael Schmid and Gareth S. Parkinson},
doi = {10.1021/acs.jpcc.3c03684},
year = {2023},
date = {2023-09-11},
urldate = {2023-09-11},
journal = {Journal of Physical Chemistry C},
volume = {127},
issue = {37},
pages = {18378--18388},
publisher = {American Chemical Society (ACS)},
abstract = {The adsorption/desorption of ethene (C_{2}H_{4}), also commonly known as ethylene, on Fe_{3}O_{4}(001) was studied under ultrahigh vacuum conditions using temperature-programmed desorption (TPD), scanning tunneling microscopy, X-ray photoelectron spectroscopy, and density functional theory (DFT)-based computations. To interpret the TPD data, we have employed a new analysis method based on equilibrium thermodynamics. C_{2}H_{4} adsorbs intact at all coverages and interacts most strongly with surface defects such as antiphase domain boundaries and Fe adatoms. On the regular surface, C_{2}H_{4} binds atop surface Fe sites up to a coverage of 2 molecules per (√2 × √2)R45° unit cell, with every second Fe occupied. A desorption energy of 0.36 eV is determined by analysis of the TPD spectra at this coverage, which is approximately 0.1–0.2 eV lower than the value calculated by DFT + U with van der Waals corrections. Additional molecules are accommodated in between the Fe rows. These are stabilized by attractive interactions with the molecules adsorbed at Fe sites. The total capacity of the surface for C_{2}H_{4} adsorption is found to be close to 4 molecules per (√2 × √2)R45° unit cell.},
keywords = {P02, P04, P07},
pubstate = {published},
tppubtype = {article}
}
Redondo, Jesús; Michalička, Jan; Kraushofer, Florian; Franceschi, Giada; Šmid, Břetislav; Kumar, Nishant; Man, Ondřej; Blatnik, Matthias; Wrana, Dominik; Mallada, Benjamin; Švec, Martin; Parkinson, Gareth S.; Setvin, Martin; Riva, Michele; Diebold, Ulrike; Čechal, Jan
Journal ArticleOpen AccessIn: Advanced Materials Interfaces, no. 2300602, 2023.
Abstract | Links | BibTeX | Tags: P02, P04
@article{Redondo2023,
title = {Hematite α-Fe_{2}O_{3}(0001) in Top and Side View: Resolving Long-Standing Controversies about Its Surface Structure},
author = {Jesús Redondo and Jan Michalička and Florian Kraushofer and Giada Franceschi and Břetislav Šmid and Nishant Kumar and Ondřej Man and Matthias Blatnik and Dominik Wrana and Benjamin Mallada and Martin Švec and Gareth S. Parkinson and Martin Setvin and Michele Riva and Ulrike Diebold and Jan Čechal},
doi = {10.1002/admi.202300602},
year = {2023},
date = {2023-08-18},
urldate = {2023-08-18},
journal = {Advanced Materials Interfaces},
number = {2300602},
publisher = {Wiley},
abstract = {Hematite is a common iron oxide found in nature, and the α-Fe_{2}O_{3}(0001) plane is prevalent on the nanomaterial utilized in photo- and electrocatalytic applications. The atomic-scale structure of the surface remains controversial despite decades of study, partly because it depends on sample history as well as the preparation conditions. Here, a comprehensive study is performed using an arsenal of surface techniques (non-contact atomic force microscopy, scanning tunneling microscopy, low-energy electron diffraction, and X-ray photoemission spectroscopy) complemented by analyses of the near surface region by high-resolution transmission electron microscopy and electron energy loss spectroscopy. The results show that the so-called “bi-phase” termination forms even under highly oxidizing conditions; a (1 × 1) surface is only observed in the presence of impurities. Furthermore, it is shown that the biphase is actually a continuous layer distorted due to a mismatch with the subsurface layers, and thus not the proposed mixture of FeO(111) and α-Fe_{2}O_{3}(0001) phases. Overall, the results show how combining surface and cross-sectional imaging provides a full view that can be essential for understanding the role of the near-surface region on oxide surface properties.},
keywords = {P02, P04},
pubstate = {published},
tppubtype = {article}
}
Kraushofer, Florian; Meier, Matthias; Jakub, Zdeněk; Hütner, Johanna; Balajka, Jan; Hulva, Jan; Schmid, Michael; Franchini, Cesare; Diebold, Ulrike; Parkinson, Gareth S.
Oxygen-Terminated (1 × 1) Reconstruction of Reduced Magnetite Fe3O4(111)
Journal ArticleOpen AccessIn: vol. 14, no. 13, pp. 3258–3265, 2023.
Abstract | Links | BibTeX | Tags: P02, P04, P07
@article{Kraushofer2023,
title = {Oxygen-Terminated (1 × 1) Reconstruction of Reduced Magnetite Fe_{3}O_{4}(111)},
author = {Florian Kraushofer and Matthias Meier and Zdeněk Jakub and Johanna Hütner and Jan Balajka and Jan Hulva and Michael Schmid and Cesare Franchini and Ulrike Diebold and Gareth S. Parkinson},
doi = {10.1021/acs.jpclett.3c00281},
year = {2023},
date = {2023-03-28},
urldate = {2023-03-28},
volume = {14},
number = {13},
pages = {3258--3265},
publisher = {American Chemical Society (ACS)},
abstract = {The (111) facet of magnetite (Fe_{3}O_{4}) has been studied extensively by experimental and theoretical methods, but controversy remains regarding the structure of its low-energy surface terminations. Using density functional theory (DFT) computations, we demonstrate three reconstructions that are more favorable than the accepted Feoct2 termination under reducing conditions. All three structures change the coordination of iron in the kagome Feoct1 layer to be tetrahedral. With atomically resolved microscopy techniques, we show that the termination that coexists with the Fetet1 termination consists of tetrahedral iron capped by 3-fold coordinated oxygen atoms. This structure explains the inert nature of the reduced patches.},
keywords = {P02, P04, P07},
pubstate = {published},
tppubtype = {article}
}
Corrias, Marco; Papa, Lorenzo; Sokolovíc, Igor; Birschitzky, Viktor; Gorfer, Alexander; Setvin, Martin; Schmid, Michael; Diebold, Ulrike; Reticcioli, Michele; Franchini, Cesare
Automated Real-Space Lattice Extraction for Atomic Force Microscopy Images
Journal ArticleOpen AccessIn: Machine Learning: Science and Technology, vol. 4, pp. 015015, 2023.
Abstract | Links | BibTeX | Tags: P02, P07
@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}
}
2022
Schmid, Michael; Parkinson, Gareth S.; Diebold, Ulrike
Analysis of Temperature-Programmed Desorption via Equilibrium Thermodynamics
Journal ArticleOpen AccessIn: ACS Physical Chemistry Au, vol. 3, iss. 1, pp. 44–62, 2022.
Abstract | Links | BibTeX | Tags: P02, P04
@article{Schmid2022,
title = {Analysis of Temperature-Programmed Desorption via Equilibrium Thermodynamics},
author = {Michael Schmid and Gareth S. Parkinson and Ulrike Diebold},
doi = {10.1021/acsphyschemau.2c00031},
year = {2022},
date = {2022-11-15},
journal = {ACS Physical Chemistry Au},
volume = {3},
issue = {1},
pages = {44--62},
publisher = {American Chemical Society (ACS)},
abstract = {Temperature-programmed desorption (TPD) experiments in surface science are usually analyzed using the Polanyi–Wigner equation and/or transition-state theory. These methods are far from straightforward, and the determination of the pre-exponential factor is often problematic. We present a different method based on equilibrium thermodynamics, which builds on an approach previously used for TPD by Kreuzer et al. (\textit{Surf. Sci.}\textbf{1988}). Equations for the desorption rate are presented for three different types of surface–adsorbate interactions: (i) a 2D ideal hard-sphere gas with a negligible diffusion barrier, (ii) an ideal lattice gas, that is, fixed adsorption sites without interaction between the adsorbates, and (iii) a lattice gas with a distribution of (site-dependent) adsorption energies. We show that the coverage dependence of the sticking coefficient for adsorption at the desorption temperature determines whether the desorption process can be described by first- or second-order kinetics. The sticking coefficient at the desorption temperature must also be known for a quantitative determination of the adsorption energy, but it has a rather weak influence (like the pre-exponential factor in a traditional TPD analysis). Quantitative analysis is also influenced by the vibrational contributions to the energy and entropy. For the case of a single adsorption energy, we provide equations to directly convert peak temperatures into adsorption energies. These equations also provide an approximation of the desorption energy in cases that cannot be described by a fixed pre-exponential factor. For the case of a distribution of adsorption energies, the desorption spectra cannot be considered a superposition of desorption spectra corresponding to the different energies. Nevertheless, we present a method to extract the distribution of adsorption energies from TPD spectra, and we rationalize the energy resolution of TPD experiments. The analytical results are complemented by a program for simulation and analysis of TPD data.},
keywords = {P02, P04},
pubstate = {published},
tppubtype = {article}
}
Wang, Zhichang; Reticcioli, Michele; Jakub, Zdenek; Sokolović, Igor; Meier, Matthias; Boatner, Lynn A; Schmid, Michael; Parkinson, Gareth S.; Diebold, Ulrike; Franchini, Cesare; Setvin, Martin
Surface chemistry on a polarizable surface: Coupling of CO with KTaO 3(001)
Journal ArticleOpen AccessIn: Science Advances, vol. 8, iss. 33, 2022.
Abstract | Links | BibTeX | Tags: P02, P04, P07
@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}
}