Bjørk Hammer

Electronic factors determining the reactivity of metal surfaces

Abstract Based on density functional theory calculations of H 2 dissociation on Al(111), Cu(111), Pt(111) and Cu 3 Pt(111) we present a consistent picture of some key physical properties determining the reactivity of metal and alloy surfaces. The four metal surfaces are chosen to represent metals with no t -bands, with filled d -bands and with d -states at the Fermi level. We show that electronic states in the entire valence band of the metal surface are responsible for the reactivity, which consequently cannot be understood solely in terms of the density of states at the Fermi nor in terms d -states above it. Rather we suggest that trends in reactivities can be understood in terms of the hybridization energy between the bonding and anti-bonding adsorbate states and the metal d -bands (when present), and we demonstrate that a simple frozen potential based estimate of the hybridization energy correlates well with the calculated variation of the barrier height for the different metal surfaces.

Theoretical Surface Science and Catalysis: Calculations and Concepts

Abstract The application of density functional theory to calculate adsorption properties, reaction pathways, and activation energies for surface chemical reactions is reviewed. Particular emphasis is placed on developing concepts that can be used to understand and predict variations in reactivity from one transition metal to the next or the effects of alloying, surface structure, and adsorbate-adsorbate interactions on the reactivity. Most examples discussed are concerned with the catalytic properties of transition metal surfaces, but it is shown that the calculational approach and the concepts developed to understand trends in reactivity for metals can also be used for sulfide and oxide catalysts.

Bandgap opening in graphene induced by patterned hydrogen adsorption

Graphene, a single layer of graphite, has recently attracted considerable attention owing to its remarkable electronic and structural properties and its possible applications in many emerging areas such as graphene-based electronic devices. The charge carriers in graphene behave like massless Dirac fermions, and graphene shows ballistic charge transport, turning it into an ideal material for circuit fabrication. However, graphene lacks a bandgap around the Fermi level, which is the defining concept for semiconductor materials and essential for controlling the conductivity by electronic means. Theory predicts that a tunable bandgap may be engineered by periodic modulations of the graphene lattice, but experimental evidence for this is so far lacking. Here, we demonstrate the existence of a bandgap opening in graphene, induced by the patterned adsorption of atomic hydrogen onto the Moiré superlattice positions of graphene grown on an Ir(111) substrate.

Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method

Electronic structure calculations have become an indispensable tool in many areas of materials science and quantum chemistry. Even though the Kohn-Sham formulation of the density-functional theory (DFT) simplifies the many-body problem significantly, one is still confronted with several numerical challenges. In this article we present the projector augmented-wave (PAW) method as implemented in the GPAW program package (https://wiki.fysik.dtu.dk/gpaw) using a uniform real-space grid representation of the electronic wavefunctions. Compared to more traditional plane wave or localized basis set approaches, real-space grids offer several advantages, most notably good computational scalability and systematic convergence properties. However, as a unique feature GPAW also facilitates a localized atomic-orbital basis set in addition to the grid. The efficient atomic basis set is complementary to the more accurate grid, and the possibility to seamlessly switch between the two representations provides great flexibility. While DFT allows one to study ground state properties, time-dependent density-functional theory (TDDFT) provides access to the excited states. We have implemented the two common formulations of TDDFT, namely the linear-response and the time propagation schemes. Electron transport calculations under finite-bias conditions can be performed with GPAW using non-equilibrium Green functions and the localized basis set. In addition to the basic features of the real-space PAW method, we also describe the implementation of selected exchange-correlation functionals, parallelization schemes, ΔSCF-method, x-ray absorption spectra, and maximally localized Wannier orbitals.

Surface electronic structure and reactivity of transition and noble metals

We present self-consistent density functional calculations using the LMTO-ASA method of the variations in the surface electronic structure for pseudomorfic overlayers and impurities of Fe, Co, Ni, Cu, Ru, Rh, Pd, Ag, Ir, Pt, and Au on the other metals. Knowledge of these variations is of importance in understanding trends in the reactivity of metal surfaces. A simple model is presented which gives a description of the overall trends in the self-consistently calculated results.

The Role of Interstitial Sites in the Ti3d Defect State in the Band Gap of Titania

Titanium dioxide (TiO2) has a number of uses in catalysis, photochemistry, and sensing that are linked to the reducibility of the oxide. Usually, bridging oxygen (Obr) vacancies are assumed to cause the Ti3d defect state in the band gap of rutile TiO2(110). From high-resolution scanning tunneling microscopy and photoelectron spectroscopy measurements, we propose that Ti interstitials in the near-surface region may be largely responsible for the defect state in the band gap. We argue that these donor-specific sites play a key role in and may dictate the ensuing surface chemistry, such as providing the electronic charge required for O2 adsorption and dissociation. Specifically, we identified a second O2 dissociation channel that occurs within the Ti troughs in addition to the O2 dissociation channel in Obr vacancies. Comprehensive density functional theory calculations support these experimental observations.

Enhanced Bonding of Gold Nanoparticles on Oxidized TiO2(110)

We studied the nucleation of gold clusters on TiO2(110) surfaces in three different oxidation states by high-resolution scanning tunneling microscopy. The three TiO2(110) supports chosen were (i) reduced (having bridging oxygen vacancies), (ii) hydrated (having bridging hydroxyl groups), and (iii) oxidized (having oxygen adatoms). At room temperature, gold nanoclusters nucleate homogeneously on the terraces of the reduced and oxidized supports, whereas on the hydrated TiO2(110) surface, clusters form preferentially at the step edges. From interplay with density functional theory calculations, we identified two different gold-TiO2(110) adhesion mechanisms for the reduced and oxidized supports. The adhesion of gold clusters is strongest on the oxidized support, and the implications of this finding for catalytic applications are discussed.

Structure sensitivity in adsorption: CO interaction with stepped and reconstructed Pt surfaces

The structure sensitivity of CO adsorption on different flat, stepped, kinked and reconstructed Pt surfaces is studied using large-scale density-functional calculations. We find an extremely strong structure sensitivity in the adsorption energy with variations up to 1 eV (or 100%) from one structure to the next. We propose a model to explain this behavior, and use it to discuss more generally the origin of structure sensitivity in heterogeneous catalysis.

A theoretical study of CH4 dissociation on pure and gold‐alloyed Ni(111) surfaces

We present a density functional theory study of the first step of CH4 adsorption on the Ni(111) surface, dissociation into adsorbed CH3 and H. The rupture of the C–H bond occurs preferentially on top of a Ni atom, with a dissociation barrier of about 100 kJ/mol (including zero point corrections). The transition state involves considerable internal excitation of the molecule. The active C–H bond is both stretched to 1.6 A and tilted relative to the methyl group. A normal mode analysis shows that the reaction coordinate is mainly a C–H stretch, while the orientation of the C–H bond relative to the surface is responsible for the highest real mode. Alloying the surface with gold also affects the reactivity of the Ni atoms on adjacent surface sites. The dissociation barrier is increased by 16 and 38 kJ/mol for a Ni atom with one or two gold neighbors, respectively. We attribute these changes to a shift in the local density of d states at the nickel atoms in the neighborhood of gold.

Special Sites at Noble and Late Transition Metal Catalysts

An overview of recent advancements in density functional theory modeling of particularly reactive sites at noble and late transition metal surfaces is given. Such special sites include sites at the flat surfaces of thin metal films, sites at stepped surfaces, sites at the metal/oxide interface boundary for oxide-supported metal clusters, and sites at the perimeter of oxide islands grown on metal surfaces. The Newns–Anderson model of the electronic interaction underlying chemisorption is described. This provides the grounds for introducing the Hammer–Nørskov d-band model that correlates changes in the energy center of the valence d-band density of states at the surface sites with their ability to form chemisorption bonds. A reactivity change described by this model is characterized as an electronic structure effect. Brønsted plots of energy barriers versus reaction energies are discussed from the surface reaction perspective and are used to analyze the trends in the calculated changes. Deviations in the relation between energy barriers and reaction energies in Brønsted plots are identified as due to atomic structure effects. The reactivity change from pure Pd surfaces to Pd thin films supported on MgO can be assigned to an electronic effect. Likewise for the reactivity change from flat Au surfaces, over Au thin films to Au edges and the Au/MgO interface boundary. The reactivity enhancement at atomic step sites is of both electronic and atomic structure nature for NO dissociation at Ru, Rh and Pd surfaces. The enhancement of the CO oxidation reactivity when moving from a CO+O coadsorption structure on Pt(111) to the PtO2 oxide island edges supported by Pt(111) is, however, identified as mainly an atomic structure effect. As such, it is linked to the occurrence of favorable pathways at the oxide island edges and is occurring despite of stronger adsorbate binding of the oxygen within the oxide edge, i.e. despite of an opposing electronic effect. As a final topic, a discussion is given of the accuracy of density functional theory in conjunction with surface reactions; adsorption, desorption, diffusion, and dissociation. Energy barriers are concluded to be more robust with respect to changes in the exchange-correlation functional than are molecular bond and adsorption energies.