Dynamics in the O(2 × 1) adlayer on Ru(0001): bridging timescales from milliseconds to minutes by scanning tunneling microscopy.

The dynamics within an O(2 × 1) adlayer on Ru(0001) is studied by density functional theory and high-speed scanning tunneling microscopy. Transition state theory proposes dynamic oxygen species in the reduced O(2 × 1) layer at room temperature. Collective diffusion processes can result in structural reorientations of characteristic stripe patterns. Spiral high-speed scanning tunneling microscopy measurements reveal this reorientation as a function of time in real space. Measurements, ranging over several minutes with constantly high frame rates of 20 Hz resolved the gradual reorientation. Moreover, reversible fast flipping events of stripe patterns are observed. These measurements relate the observations of long-term atomic rearrangements and their underlying fast processes captured within several tens of milliseconds.

Structure and registry of the silica bilayer film on Ru(0001) as viewed by LEED and DFT.

Silica bilayers are stable on various metal substrates, including Ru(0001) that is used for the present study. In a systematic attempt to elucidate the detailed structure of the silica bilayer film and its registry to the metal substrate, we performed a low energy electron diffraction (I/V-LEED) study. The experimental work is accompanied by detailed calculations on the stability, orientation and dynamic properties of the bilayer at room temperature. It was determined, that the film shows a certain structural diversity within the unit cell of the metal substrate, which depends on the oxygen content at the metal-bilayer interface. In connection with the experimental I/V-LEED study, it became apparent, that a high-quality structure determination is only possible if several structural motifs are taken into account by superimposing bilayer structures with varying registry to the oxygen covered substrate. This result is conceptually in line with the recently observed statistical registry in layered 2D-compound materials.

Adsorption of CH4 on the Pt(111) surface: Random phase approximation compared to density functional theory.

We investigate the adsorption of CH4 on the Pt(111) surface for two adsorption modes, hcp (hexagonal closed packed) hollow tripod and top monopod in a (√3 × âˆš3)R30° surface cell that corresponds to experimental surface coverage. Surface structures are optimized with density functional theory using the Perdew-Burke-Ernzerhof (PBE) functional augmented with the many-body dispersion scheme of Tkatchenko (PBE+MBD). Whereas the Random Phase Approximation (RPA) predicts a clear preference of about 5 kJ mol-1 for the hcp tripod compared to the top monopod structure, in agreement with vibrational spectra, PBE+MBD predicts about equal stability for the two adsorption structures. For the hcp tripod, RPA yields an adsorption energy of -14.5 kJ mol-1, which is converged to within 1.0 ± 0.5 kJ mol-1 with respect to the plane wave energy cutoff (500 eV), the k-point mesh (4 × 4 × 1), the vacuum layer (about 10.3 Å, with extrapolation to infinite distance), and the number of Pt layers (3). Increments for increasing the number of Pt layers to 4 (+1.6 kJ mol-1) and the k-point mesh to 6 × 6 × 1 (-0.6 kJ mol-1) yield a final estimate of -13.5 ± 2.1 kJ mol-1, which agrees to within 2.2 ± 2.1 kJ mol-1 with experiment (-15.7 ± 1.6), well within the chemical accuracy range.

Elucidating Surface Structure with Action Spectroscopy.

Surface Action Spectroscopy, a vibrational spectroscopy method developed in recent years at the Fritz Haber Institute is employed for structure determination of clean and H2O-dosed (111) magnetite surfaces. Surface structural information is revealed by using the microscopic surface vibrations as a fingerprint of the surface structure. Such vibrations involve just the topmost atomic layers, and therefore the structural information is truly surface related. Our results strongly support the view that regular Fe3O4(111)/Pt(111) is terminated by the so-called Fetet1 termination, that the biphase termination of Fe3O4(111)/Pt(111) consists of FeO and Fe3O4(111) terminated areas, and we show that the method can differentiate between different water structures in H2O-derived adsorbate layers on Fe3O4(111)/Pt(111). With this, we conclude that the method is a capable new member in the set of techniques providing crucial information to elucidate surface structures. The method does not rely on translational symmetry and can therefore also be applied to systems which are not well ordered. Even an application to rough surfaces is possible.

Vibrational properties of CO2 adsorbed on the Fe3O4 (111) surface: Insights gained from DFT.

By virtue of density functional theory calculations, this work discusses several carbonate, carboxylate, and bicarbonate species on two thermodynamically relevant metal terminations of the (111) surface of magnetite, Fe3O4. We present adsorption energies and vibrational wavenumbers and conclude in assigning the observed infrared reflection-absorption spectroscopy bands. CO2 prefers to adsorb molecularly on the Fetet1 terminated Fe3O4(111) surface, a finding consistent with observation. Calculations compared with the experiment lead to interpreting results in favor of the Fetet1 (single metal) terminated Fe3O4(111) surface as the regular surface termination. Formation of carbonate and bicarbonate requires metal impurities on that surface. Such impurities exist, for instance, on the Feoct2 (double metal) termination, which can thus be used as a model for "metal-rich terminations" of more complex surfaces.

Water adsorption on the Fe3O4(111) surface: dissociation and network formation.

We monitored adsorption of water on a well-defined Fe3O4(111) film surface at different temperatures as a function of coverage using infrared reflection-absorption spectroscopy, temperature programmed desorption, and single crystal adsorption calorimetry. Additionally, density functional theory was employed using a Fe3O4(111)-(2 × 2) slab model to generate 15 energy minimum structures for various coverages. Corresponding vibrational properties of the adsorbed water species were also computed. The results show that water molecules readily dissociate on regular surface Fetet1-O ion pairs to form "monomers", i.e., terminal Fe-OH and surface OH groups. Further water molecules adsorb on the hydroxyl covered surface non-dissociatively and form "dimers" and larger oligomers, which ultimately assemble into an ordered (2 × 2) hydrogen-bonded network structure with increasing coverage prior to the formation of a solid water film.

Cooperative Formation of Long-Range Ordering in Water Ad-layers on Fe3 O4 (111) Surfaces.

The initial stages of water adsorption on magnetite Fe3 O4 (111) surface and the atomic structure of the water/oxide interface remain controversial. Herein, we provide experimental results obtained by infrared reflection-absorption spectroscopy (IRAS) and temperature-programmed desorption (TPD), corroborated by density functional theory (DFT) calculations showing that water readily dissociates on Fetet sites to form two hydroxo species. These act as an anchor for water molecules to form a dimer complex which self-assembles into an ordered (2×2) structure. Water ad-layer ordering is rationalized in terms of a cooperative effect induced by a hydrogen-bonding network.

Toward an Understanding of Selective Alkyne Hydrogenation on Ceria: On the Impact of O Vacancies on H2 Interaction with CeO2(111).

Ceria (CeO2) has recently been found to be a promising catalyst in the selective hydrogenation of alkynes to alkenes. This reaction occurs primarily on highly dispersed metal catalysts, but rarely on oxide surfaces. The origin of the outstanding activity and selectivity observed on CeO2 remains unclear. In this work, we show that one key aspect of the hydrogenation reaction-the interaction of hydrogen with the oxide-depends strongly on the presence of O vacancies within CeO2. Through infrared reflection absorption spectroscopy on well-ordered CeO2(111) thin films and density functional theory (DFT) calculations, we show that the preferred heterolytic dissociation of molecular hydrogen on CeO2(111) requires H2 pressures in the mbar regime. Hydrogen depth profiling with nuclear reaction analysis indicates that H species stay on the surface of stoichiometric CeO2(111) films, whereas H incorporates as a volatile species into the volume of partially reduced CeO2-x(111) thin films (x ∼ 1.8-1.9). Complementary DFT calculations demonstrate that oxygen vacancies facilitate H incorporation below the surface and that they are the key to the stabilization of hydridic H species in the volume of reduced ceria.

Reduction and oxidation of Au adatoms on the CeO2(111) surface - DFT+U versus hybrid functionals.

Recently we showed that Au atoms may titrate Ce3+ ions in near-surface layers of reduced CeO2(111). This surface contained oxygen vacancies in subsurface position within the topmost O-Ce-O trilayer [Pan et al., Phys. Rev. Lett., 2013, 111, 206101.]. The present work builds upon these findings and discusses additional results obtained using PBE+U and hybrid functionals. These approaches do not predict the same relative stabilities for the various adsorption sites of a single Au adatom at an O-defect concentration of a » ML or 1.984 nm-2. We attribute this discrepancy to a different alignment within the O 2p-Ce 4f gap, i.e. a different order by energy of partially occupied Ce 4f and Au 6s orbitals. The energy offset of these orbitals matters, because the adsorption of Au0(6s1) atop Ce3+(4f1) or atop a subsurface oxygen atom in the first coordination shell of a Ce3+(4f1) involves creation of Au-(6s2) and Ce4+(4f0) ions. The electron transfer to Au is coupled to stabilizing ionic relaxation in the lattice, commonly known as polaronic distortion, reinforcing the Au-Ce bond. The order of 4f and 6s orbitals depends on the density functional approximation and is also strongly influenced by the oxygen defect concentration.

O2 Activation on Ceria Catalysts-The Importance of Substrate Crystallographic Orientation.

An atomic-level understanding of dioxygen activation on metal oxides remains one of the major challenges in heterogeneous catalysis. By performing a thorough surface-science study of all three low-index single-crystal surfaces of ceria, probably the most important redox catalysts, we provide a direct spectroscopic characterization of reactive dioxygen species at defect sites on the reduced ceria (110) and (100) surfaces. Surprisingly, neither of these superoxo and peroxo species was found on ceria (111), the thermodynamically most stable surface of this oxide. Applying density functional theory, we could relate these apparently inconsistent findings to a sub-surface diffusion of O vacancies on (111) substrates, but not on the less-closely packed surfaces. These observations resolve a long standing debate concerning the location of O vacancies on ceria surfaces and the activation of O2 on ceria powders.

Structural and Electronic Effects on the Properties of Fe2(dobdc) upon Oxidation with N2O.

We report electronic, vibrational, and magnetic properties, together with their structural dependences, for the metal-organic framework Fe2(dobdc) (dobdc(4-) = 2,5-dioxido-1,4-benzenedicarboxylate) and its derivatives, Fe2(O)2(dobdc) and Fe2(OH)2(dobdc)-species arising in the previously proposed mechanism for the oxidation of ethane to ethanol using N2O as an oxidant. Magnetic susceptibility measurements reported for Fe2(dobdc) in an earlier study and reported in the current study for Fe(II)0.26[Fe(III)(OH)]1.74(dobdc)(DMF)0.15(THF)0.22, which is more simply referred to as Fe2(OH)2(dobdc), were used to confirm the computational results. Theory was also compared to experiment for infrared spectra and powder X-ray diffraction structures. Structural and magnetic properties were computed by using Kohn-Sham density functional theory both with periodic boundary conditions and with cluster models. In addition, we studied the effects of different treatments of the exchange interactions on the magnetic coupling parameters by comparing several approaches to the exchange-correlation functional: generalized gradient approximation (GGA), GGA with empirical Coulomb and exchange integrals for 3d electrons (GGA+U), nonseparable gradient approximation (NGA) with empirical Coulomb and exchange integrals for 3d electrons (NGA+U), hybrid GGA, meta-GGA, and hybrid meta-GGA. We found the coupling between the metal centers along a chain to be ferromagnetic in the case of Fe2(dobdc) and antiferromagnetic in the cases of Fe2(O)2(dobdc) and Fe2(OH)2(dobdc). The shift in magnetic coupling behavior correlates with the changing electronic structure of the framework, which derives from both structural and electronic changes that occur upon metal oxidation and addition of the charge-balancing oxo and hydroxo ligands.

Erratum: Titration of Ce

This corrects the article DOI: 10.1103/PhysRevLett.111.206101.

Surface structure of V_{2}O_{3}(0001) revisited.

In a recent paper [A. J. Window et al., Phys. Rev. Lett. 107, 016105 (2011)], it was proposed that V_{2}O_{3}(0001) is terminated by the so-called O_{3} termination, a reconstruction with a terminating distorted hexagonal oxygen layer. We show that the surface is terminated by vanadyl (V═O) groups instead. This conclusion is based on quantitative low-energy electron diffraction combined with scanning tunneling microscopy, fast atom scattering, and density functional theory employing the Heyd-Scuseria-Ernzerhof functional. New insights into the subsurface sensitivity of ion beam triangulation show that results previously interpreted in favor of the O_{3} termination are reconcilable with vanadyl termination as well.

Water Interaction with Iron Oxides.

We present a mechanistic study on the interaction of water with a well-defined model Fe3O4(111) surface that was investigated by a combination of direct calorimetric measurements of adsorption energies, infrared vibrational spectroscopy, and calculations bases on density functional theory (DFT). We show that the adsorption energy of water (101 kJ mol(-1)) is considerably higher than all previously reported values obtained by indirect desorption-based methods. By employing (18)O-labeled water molecules and an Fe3 O4 substrate, we proved that the generally accepted simple model of water dissociation to form two individual OH groups per water molecule is not correct. DFT calculations suggest formation of a dimer, which consists of one water molecule dissociated into two OH groups and another non-dissociated water molecule creating a thermodynamically very stable dimer-like complex.

Support effect in oxide catalysis: methanol oxidation on vanadia/ceria.

Density functional theory is used for periodic models of monomeric vanadia species deposited on the CeO2(111) surface to study dissociative adsorption of methanol and its subsequent dehydrogenation to formaldehyde. Dispersion-corrected PBE+U calculations are performed and compared with HSE and B3LYP results. Dissociative adsorption of methanol at different sites on VO2·CeO2(111) is highly exothermic with adsorption energies of 1.8 to 1.9 eV (HSE+D). Two relevant pathways for desorption of formaldehyde are found with intrinsic barriers for the redox step of 1.0 and 1.4 eV (HSE+D). The calculated desorption temperatures (370 and 495 K) explain the peaks observed in temperature-programmed desorption experiments. Different sites of the supported catalyst system are involved in the two pathways: (i) methanol can chemisorb on the CeO2 surface filling a so-called pseudovacancy and the H atom is transferred to an V-O-Ce interphase bond or (ii) CH3OH may chemisorb at the V-O-Ce interphase bond and form a V-OCH3 species from which H is transferred to the ceria surface, providing evidence for true cooperativity. In both cases, ceria is directly involved in the redox process, as two electrons are accommodated in Ce f states forming two Ce(3+) ions whereas vanadium remains fully oxidized (V(5+)).

Sites for methane activation on lithium-doped magnesium oxide surfaces.

Density functional calculations yield energy barriers for H abstraction by oxygen radical sites in Li-doped MgO that are much smaller (12±6 kJ mol(-1)) than the barriers inferred from different experimental studies (80-160 kJ mol(-1)). This raises further doubts that the Li(+)O(˙-) site is the active site as postulated by Lunsford. From temperature-programmed oxidative coupling reactions of methane (OCM), we conclude that the same sites are responsible for the activation of CH4 on both Li-doped MgO and pure MgO catalysts. For a MgO catalyst prepared by sol-gel synthesis, the activity proved to be very different in the initial phase of the OCM reaction and in the steady state. This was accompanied by substantial morphological changes and restructuring of the terminations as transmission electron microscopy revealed. Further calculations on cluster models showed that CH4 binds heterolytically on Mg(2+)O(2-) sites at steps and corners, and that the homolytic release of methyl radicals into the gas phase will happen only in the presence of O2.

Hybrid RPA:DFT Approach for Adsorption on Transition Metal Surfaces: Methane and Ethane on Platinum (111).

The hybrid QM:QM approach is extended to adsorption on transition metal surfaces. The random phase approximation (RPA) as the high-level method is applied to cluster models and, using the subtractive scheme, embedded in periodic models which are treated with density functional theory (DFT) that is the low-level method. The PBE functional, both without dispersion and augmented with the many-body dispersion (MBD), is employed. Adsorption of methane and ethane on the Pt(111) surface is studied. For methane in a 2 × 2 surface cell, the hybrid RPA:PBE and RPA:PBE+MBD results, -14.3 and -16.0 kJ mol-1, respectively, are in close agreement with the periodic RPA value of -13.8 kJ mol-1 at significantly reduced computational cost (factor of ∼50). For methane and ethane, the RPA:PBE results (-14.3 and -17.8 kJ mol-1, respectively) indicate underbinding relative to energies derived from experimental desorption barriers for relevant loadings (-15.6 ± 1.6 and -27.2 ± 2.9 kJ mol-1, respectively), whereas the hybrid RPA:PBE+MBD results (-16.0 and -24.9 kJ mol-1, respectively) agree with the experiment well within experimental uncertainty limits (deviation of -0.4 ± 1.5 and +2.3 ± 2.9 kJ mol-1, respectively). Finding a cluster that adequately and robustly represents the adsorbate at the bulk surface is important for the success of the RPA-based QM:QM scheme for metals.

Titration of Ce3+ ions in the CeO2(111) surface by Au adatoms.

The role of surface and subsurface O vacancies for gold adsorption on crystalline CeO2(111) films has been investigated by scanning tunneling microscopy and density functional theory. Whereas surface vacancies serve as deep traps for the Au atoms, subsurface defects promote the formation of characteristic Au pairs with a mean atom distance of two ceria lattice constants (7.6 Å). Hybrid density functional theory calculations reveal that the pair formation arises from a titration of the two Ce3+ ions generated by a single O vacancy. The Au-Ce3+ bond forms due to a strain effect, as the associated charge transfer from the spacious Ce3+ into the adgold enables a substantial relaxation of the ceria lattice. Also the experimentally determined Au-pair length is reproduced in the calculations, as we find a Ce3+-Ce3+ spacing of two ceria lattice parameters to be energetically preferred. Single Au atoms can thus be taken as position markers for Ce3+ ion pairs in the surface, providing unique information on electron-localization phenomena in reduced ceria.

Van der Waals interactions in ionic and semiconductor solids.

van der Waals (vdW) energy corrected density-functional theory [Phys. Rev. Lett. 102, 073005 (2009)] is applied to study the cohesive properties of ionic and semiconductor solids (C, Si, Ge, GaAs, NaCl, and MgO). The required polarizability and dispersion coefficients are calculated using the dielectric function obtained from time-dependent density-functional theory. Coefficients for "atoms in the solid" are then calculated from the Hirshfeld partitioning of the electron density. It is shown that the Clausius-Mossotti equation that relates the polarizability and the dielectric function is accurate even for covalently-bonded semiconductors. We find an overall improvement in the cohesive properties of Si, Ge, GaAs, NaCl, and MgO, when vdW interactions are included on top of the Perdew-Burke-Ernzerhof or Heyd-Scuseria-Ernzerhof functionals. The relevance of our findings for other solids is discussed.