Oxygen-Terminated (1 × 1) Reconstruction of Reduced Magnetite Fe3O4(111).

The (111) facet of magnetite (Fe3O4) 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.

Analysis of Temperature-Programmed Desorption via Equilibrium Thermodynamics.

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. (Surf. Sci. 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.

Resolving the intrinsic short-range ordering of K+ ions on cleaved muscovite mica.

Muscovite mica, KAl2(Si3Al)O10(OH)2, is a common layered phyllosilicate with perfect cleavage planes. The atomically flat surfaces obtained through cleaving lend themselves to scanning probe techniques with atomic resolution and are ideal to model minerals and clays. Despite the importance of the cleaved mica surfaces, several questions remain unresolved. It is established that K+ ions decorate the cleaved surface, but their intrinsic ordering - unaffected by the interaction with the environment - is not known. This work presents clear images of the K+ distribution of cleaved mica obtained with low-temperature non-contact atomic force microscopy (AFM) under ultra-high vacuum (UHV) conditions. The data unveil the presence of short-range ordering, contrasting previous assumptions of random or fully ordered distributions. Density functional theory (DFT) calculations and Monte Carlo simulations show that the substitutional subsurface Al3+ ions have an important role for the surface K+ ion arrangement.

Water Structures Reveal Local Hydrophobicity on the In2O3(111) Surface.

Clean oxide surfaces are generally hydrophilic. Water molecules anchor at undercoordinated surface metal atoms that act as Lewis acid sites, and they are stabilized by H bonds to undercoordinated surface oxygens. The large unit cell of In2O3(111) provides surface atoms in various configurations, which leads to chemical heterogeneity and a local deviation from this general rule. Experiments (TPD, XPS, nc-AFM) agree quantitatively with DFT calculations and show a series of distinct phases. The first three water molecules dissociate at one specific area of the unit cell and desorb above room temperature. The next three adsorb as molecules in the adjacent region. Three more water molecules rearrange this structure and an additional nine pile up above the OH groups. Despite offering undercoordinated In and O sites, the rest of the unit cell is unfavorable for adsorption and remains water-free. The first water layer thus shows ordering into nanoscopic 3D water clusters separated by hydrophobic pockets.

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.

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.

Why and How Savitzky-Golay Filters Should Be Replaced.

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.

CO oxidation by Pt2/Fe3O4: Metastable dimer and support configurations facilitate lattice oxygen extraction.

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/Fe3O4(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 CO2. Pt2 dimers produce one CO2 molecule before falling apart into two adatoms, releasing the second CO. Olattice extraction only becomes facile when both the Pt-dimer and the Fe3O4 support can access metastable configurations, suggesting that substantial, concerted rearrangements of both cluster and support must be considered for reactions occurring at elevated temperature.

Structure of an Ultrathin Oxide on Pt3 Sn(111) Solved by Machine Learning Enhanced Global Optimization.

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 Pt3 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.

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 parameterUused to correct the electronic correlation in the DFT +Uformalism 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.

Single Rh Adatoms Stabilized on α-Fe2O3(11Ì 02) by Coadsorbed Water.

Oxide-supported single-atom catalysts are commonly modeled as a metal atom substituting surface cation sites in a low-index surface. Adatoms with dangling bonds will inevitably coordinate molecules from the gas phase, and adsorbates such as water can affect both stability and catalytic activity. Herein, we use scanning tunneling microscopy (STM), noncontact atomic force microscopy (ncAFM), and X-ray photoelectron spectroscopy (XPS) to show that high densities of single Rh adatoms are stabilized on α-Fe2O3(11̅02) in the presence of 2 × 10-8 mbar of water at room temperature, in marked contrast to the rapid sintering observed under UHV conditions. Annealing to 50 °C in UHV desorbs all water from the substrate leaving only the OH groups coordinated to Rh, and high-resolution ncAFM images provide a direct view into the internal structure. We provide direct evidence of the importance of OH ligands in the stability of single atoms and argue that their presence should be assumed when modeling single-atom catalysis systems.

Rapid oxygen exchange between hematite and water vapor.

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" ([Formula: see text]) surface of hematite (α-Fe2O3) 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.

Direct assessment of the acidity of individual surface hydroxyls.

The state of deprotonation/protonation of surfaces has far-ranging implications in chemistry, from acid-base catalysis1 and the electrocatalytic and photocatalytic splitting of water2, to the behaviour of minerals3 and biochemistry4. An entity's acidity is described by its proton affinity and its acid dissociation constant pKa (the negative logarithm of the equilibrium constant of the proton transfer reaction in solution). The acidity of individual sites is difficult to assess for solids, compared with molecules. For mineral surfaces, the acidity is estimated by semi-empirical concepts, such as bond-order valence sums5, and increasingly modelled with first-principles molecular dynamics simulations6,7. At present, such predictions cannot be tested-experimental measures, such as the point of zero charge8, integrate over the whole surface or, in some cases, individual crystal facets9. Here we assess the acidity of individual hydroxyl groups on In2O3(111)-a model oxide with four different types of surface oxygen atom. We probe the strength of their hydrogen bonds with the tip of a non-contact atomic force microscope and find quantitative agreement with density functional theory calculations. By relating the results to known proton affinities of gas-phase molecules, we determine the proton affinity of the different surface sites of In2O3 with atomic precision. Measurements on hydroxylated titanium dioxide and zirconium oxide extend our method to other oxides.

Unraveling CO adsorption on model single-atom catalysts.

Understanding how the local environment of a "single-atom" catalyst affects stability and reactivity remains a challenge. We present an in-depth study of copper1, silver1, gold1, nickel1, palladium1, platinum1, rhodium1, and iridium1 species on Fe3O4(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 Ag1-CO) or weaken it (as for Ni1-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.

IrO_{2} Surface Complexions Identified through Machine Learning and Surface Investigations.

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.

Fast low-noise transimpedance amplifier for scanning tunneling microscopy and beyond.

A transimpedance amplifier has been designed for scanning tunneling microscopy (STM). The amplifier features low noise (limited by the Johnson noise of the 1 GΩ feedback resistor at low input current and low frequencies), sufficient bandwidth for most STM applications (50 kHz at 35 pF input capacitance), a large dynamic range (0.1 pA-50 nA without range switching), and a low input voltage offset. The amplifier is also suited for placing its first stage into the cryostat of a low-temperature STM, minimizing the input capacitance and reducing the Johnson noise of the feedback resistor. The amplifier may also find applications for specimen current imaging and electron-beam-induced current measurements in scanning electron microscopy and as a photodiode amplifier with a large dynamic range. This paper also discusses the sources of noise including the often neglected effect of non-balanced input impedance of operational amplifiers and describes how to accurately measure and adjust the frequency response of low-current transimpedance amplifiers.

Movable holder for a quartz crystal microbalance for exact growth rates in pulsed laser deposition.

Controlling the amount of material deposited by pulsed laser deposition (PLD) down to fractions of one atomic layer is crucial for nanoscale technologies based on thin-film heterostructures. Albeit unsurpassed for measuring growth rates with high accuracy, the quartz crystal microbalance (QCM) suffers from some limitations when applied to PLD. The strong directionality of the PLD plasma plume and its pronounced dependence on deposition parameters (e.g., background pressure and fluence) require that the QCM is placed at the same position as the substrate during growth. However, QCM sensors are commonly fixed off to one side of the substrate. This also entails fast degradation of the crystal, as it is constantly exposed to the ablated material. The design for a movable QCM holder discussed in this work overcomes these issues. The holder is compatible with standard transfer arms, enabling easy insertion and transfer between a PLD chamber and other adjoining vacuum chambers. The QCM can be placed at the same position as the substrate during PLD growth. Its resonance frequency is measured in vacuum at any location where it can be in contact with an electrical feedthrough, before and after deposition. We tested the design for the deposition of hematite (Fe2O3), comparing the rates derived from the QCM and from reflection high-energy electron diffraction oscillations during homoepitaxial growth.

Atomic-Scale Studies of Fe3 O4 (001) and TiO2 (110) Surfaces Following Immersion in CO2 -Acidified Water.

Difficulties associated with the integration of liquids into a UHV environment make surface-science style studies of mineral dissolution particularly challenging. Recently, we developed a novel experimental setup for the UHV-compatible dosing of ultrapure liquid water and studied its interaction with TiO2 and Fe3 O4 surfaces. Herein, we describe a simple approach to vary the pH through the partial pressure of CO2 ( p C O 2 ) in the surrounding vacuum chamber and use this to study how these surfaces react to an acidic solution. The TiO2 (110) surface is unaffected by the acidic solution, except for a small amount of carbonaceous contamination. The Fe3 O4 (001)-( 2 × 2 )R45° surface begins to dissolve at a pH 4.0-3.9 ( p C O 2 =0.8-1 bar) and, although it is significantly roughened, the atomic-scale structure of the Fe3 O4 (001) surface layer remains visible in scanning tunneling microscopy (STM) images. X-ray photoelectron spectroscopy (XPS) reveals that the surface is chemically reduced and contains a significant accumulation of bicarbonate (HCO3 - ) species. These observations are consistent with Fe(II) being extracted by bicarbonate ions, leading to dissolved iron bicarbonate complexes (Fe(HCO3 )2 ), which precipitate onto the surface when the water evaporates.

Electrochemical Stability of the Reconstructed Fe3 O4 (001) Surface.

Establishing the atomic-scale structure of metal-oxide surfaces during electrochemical reactions is a key step to modeling this important class of electrocatalysts. Here, we demonstrate that the characteristic (√2×√2)R45° surface reconstruction formed on (001)-oriented magnetite single crystals is maintained after immersion in 0.1 M NaOH at 0.20 V vs. Ag/AgCl and we investigate its dependence on the electrode potential. We follow the evolution of the surface using in situ and operando surface X-ray diffraction from the onset of hydrogen evolution, to potentials deep in the oxygen evolution reaction (OER) regime. The reconstruction remains stable for hours between -0.20 and 0.60 V and, surprisingly, is still present at anodic current densities of up to 10 mA cm-2 and strongly affects the OER kinetics. We attribute this to a stabilization of the Fe3 O4 bulk by the reconstructed surface. At more negative potentials, a gradual and largely irreversible lifting of the reconstruction is observed due to the onset of oxide reduction.