Bifunctional alloys for the electroreduction of CO2and CO

We use density functional theory to study the reduction of CO₂and CO to hydrocarbons through a formyl pathway on (111) and (211) facets of L1₂alloys with an A₃B composition.

Mechanistic insights into nitrogen fixation by nitrogenase enzymes

Biological nitrogen fixation by nitrogenase enzymes is a process that activates dinitrogen (N2) one of the most inert molecules in nature, within the confines of a living organism and at ambient conditions. Despite decades of study, there are still no complete explanations as to how this is possible. Here we describe a model of N2 reduction using the Mo-containing nitrogenase (FeMoco) that can explain the reactivity of the active site via a series of electrochemical steps that reversibly unseal a highly reactive Fe edge site. Our model can explain the 8 proton-electron transfers involved in biological ammonia synthesis within the kinetic scheme of Lowe and Thorneley, the obligatory formation of one H2 per N2 reduced, and the behavior of known inhibitors.

Strong Influence of Coadsorbate Interaction on CO Desorption Dynamics on Ru(0001) Probed by Ultrafast X-Ray Spectroscopy and Ab Initio Simulations

We show that coadsorbed oxygen atoms have a dramatic influence on the CO desorption dynamics from Ru(0001). In contrast to the precursor-mediated desorption mechanism on Ru(0001), the presence of surface oxygen modifies the electronic structure of Ru atoms such that CO desorption occurs predominantly via the direct pathway. This phenomenon is directly observed in an ultrafast pump-probe experiment using a soft x-ray free-electron laser to monitor the dynamic evolution of the valence electronic structure of the surface species. This is supported with the potential of mean force along the CO desorption path obtained from density-functional theory calculations. Charge density distribution and frozen-orbital analysis suggest that the oxygen-induced reduction of the Pauli repulsion, and consequent increase of the dative interaction between the CO 5σ and the charged Ru atom, is the electronic origin of the distinct desorption dynamics. Ab initio molecular dynamics simulations of CO desorption from Ru(0001) and oxygen-coadsorbed Ru(0001) provide further insights into the surface bond-breaking process.

Selective ultrafast probing of transient hot chemisorbed and precursor states of CO on Ru(0001)

We have studied the femtosecond dynamics following optical laser excitation of CO adsorbed on a Ru surface by monitoring changes in the occupied and unoccupied electronic structure using ultrafast soft x-ray absorption and emission. We recently reported [M. Dell'Angela et al. Science 339, 1302 (2013)] a phonon-mediated transition into a weakly adsorbed precursor state occurring on a time scale of >2 ps prior to desorption. Here we focus on processes within the first picosecond after laser excitation and show that the metal-adsorbate coordination is initially increased due to hot-electron-driven vibrational excitations. This process is faster than, but occurs in parallel with, the transition into the precursor state. With resonant x-ray emission spectroscopy, we probe each of these states selectively and determine the respective transient populations depending on optical laser fluence. Ab initio molecular dynamics simulations of CO adsorbed on Ru(0001) were performed at 1500 and 3000 K providing insight into the desorption process.

Li-O2 Kinetic Overpotentials: Tafel Plots from Experiment and First-Principles Theory

We report the current dependence of the fundamental kinetic overpotentials for Li-O2 discharge and charge (Tafel plots) that define the optimal cycle efficiency in a Li-air battery. Comparison of the unusual experimental Tafel plots obtained in a bulk electrolysis cell with those obtained by first-principles theory is semiquantitative. The kinetic overpotentials for any practical current density are very small, considerably less than polarization losses due to iR drops from the cell impedance in Li-O2 batteries. If only the kinetic overpotentials were present, then a discharge-charge voltaic cycle efficiency of ∼85% should be possible at ∼10 mA/cm(2) superficial current density in a battery of ∼0.1 m(2) total cathode area. We therefore suggest that minimizing the cell impedance is a more important problem than minimizing the kinetic overpotentials to develop higher current Li-air batteries.

Twin Problems of Interfacial Carbonate Formation in Nonaqueous Li-O2 Batteries

We use XPS and isotope labeling coupled with differential electrochemical mass spectrometry (DEMS) to show that small amounts of carbonates formed during discharge and charge of Li-O2 cells in ether electrolytes originate from reaction of Li2O2 (or LiO2) both with the electrolyte and with the C cathode. Reaction with the cathode forms approximately a monolayer of Li2CO3 at the C-Li2O2 interface, while reaction with the electrolyte forms approximately a monolayer of carbonate at the Li2O2-electrolyte interface during charge. A simple electrochemical model suggests that the carbonate at the electrolyte-Li2O2 interface is responsible for the large potential increase during charging (and hence indirectly for the poor rechargeability). A theoretical charge-transport model suggests that the carbonate layer at the C-Li2O2 interface causes a 10-100 fold decrease in the exchange current density. These twin "interfacial carbonate problems" are likely general and will ultimately have to be overcome to produce a highly rechargeable Li-air battery.

SURFACE SCIENCE: How to Power a Nanomotor

A successful nanotechnology will require nanomotors that can perform functions from switches to pumps and actuators. In their Perspective, Besenbacher and Nørskov discuss the study by Schmid et al., who show that tin islands on a copper surfaces propel themselves forward on the surface, drawing energy from alloy formation with the surface below. The system has about the same power-to-weight ratio as a car and may provide a paradigm for future nanomotors.

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.

Using scaling relations to understand trends in the catalytic activity of transition metals

A method is developed to estimate the potential energy diagram for a full catalytic reaction for a range of late transition metals on the basis of a calculation (or an experimental determination) for a single metal. The method, which employs scaling relations between adsorption energies, is illustrated by calculating the potential energy diagram for the methanation reaction and ammonia synthesis for 11 different metals on the basis of results calculated for Ru. It is also shown that considering the free energy diagram for the reactions, under typical industrial conditions, provides additional insight into reactivity trends.

Cyclic voltammograms for H on Pt(111) and Pt(100) from first principles

Cyclic voltammetry is a fundamental experimental method for characterizing electrochemical surfaces. Despite its wide use, a way to quantitatively and directly relate cyclic voltammetry to ab initio calculations has been lacking. We derive the cyclic voltammogram for H on Pt(111) and Pt(100), based solely on density functional theory calculations and standard molecular tables. By relating the gas phase adsorption energy to the electrochemical electrode potential, we provide a direct link between surface science and electrochemistry.

Scaling properties of adsorption energies for hydrogen-containing molecules on transition-metal surfaces

Density functional theory calculations are presented for CHx, x=0,1,2,3, NHx, x=0,1,2, OHx, x=0,1, and SHx, x=0,1 adsorption on a range of close-packed and stepped transition-metal surfaces. We find that the adsorption energy of any of the molecules considered scales approximately with the adsorption energy of the central, C, N, O, or S atom, the scaling constant depending only on x. A model is proposed to understand this behavior. The scaling model is developed into a general framework for estimating the reaction energies for hydrogenation and dehydrogenation reactions.

Estimations of electric field effects on the oxygen reduction reaction based on the density functional theory

By varying the external electric field in density functional theory (DFT) calculations we have estimated the impact of the local electric field in the electric double layer on the oxygen reduction reaction (ORR). Potentially, including the local electric field could change adsorption energies and barriers substantially, thereby affecting the reaction mechanism predicted for ORR on different metals. To estimate the effect of local electric fields on ORR we combine the DFT results at various external electric field strengths with a previously developed model of electrochemical reactions which fully accounts for the effect of the electrode potential. We find that the local electric field only slightly affects the output of the model. Hence, the general picture obtained without inclusion of the electric field still persists. However, for accurate predictions at oxygen reduction potentials close to the volcano top local electric field effects may be of importance.

Predicting Catalysis: Understanding Ammonia Synthesis from First-Principles Calculations

Here, we give a full account of a large collaborative effort toward an atomic-scale understanding of modern industrial ammonia production over ruthenium catalysts. We show that overall rates of ammonia production can be determined by applying various levels of theory (including transition state theory with or without tunneling corrections, and quantum dynamics) to a range of relevant elementary reaction steps, such as N(2) dissociation, H(2) dissociation, and hydrogenation of the intermediate reactants. A complete kinetic model based on the most relevant elementary steps can be established for any given point along an industrial reactor, and the kinetic results can be integrated over the catalyst bed to determine the industrial reactor yield. We find that, given the present uncertainties, the rate of ammonia production is well-determined directly from our atomic-scale calculations. Furthermore, our studies provide new insight into several related fields, for instance, gas-phase and electrochemical ammonia synthesis. The success of predicting the outcome of a catalytic reaction from first-principles calculations supports our point of view that, in the future, theory will be a fully integrated tool in the search for the next generation of catalysts.

Multidimensional effects on dissociation of N2 on Ru(0001)

The applicability of the Born-Oppenheimer approximation to molecule-metal surface reactions is presently a topic of intense debate. We have performed classical trajectory calculations on a prototype activated dissociation reaction, of N2 on Ru(0001), using a potential energy surface based on density functional theory. The computed reaction probabilities are in good agreement with molecular beam experiments. Comparison to previous calculations shows that the rotation of N2 and its motion along the surface affect the reactivity of N2 much more than nonadiabatic effects.

Bayesian error estimation in density-functional theory

We present a practical scheme for performing error estimates for density-functional theory calculations. The approach, which is based on ideas from Bayesian statistics, involves creating an ensemble of exchange-correlation functionals by comparing with an experimental database of binding energies for molecules and solids. Fluctuations within the ensemble can then be used to estimate errors relative to experiment on calculated quantities such as binding energies, bond lengths, and vibrational frequencies. It is demonstrated that the error bars on energy differences may vary by orders of magnitude for different systems in good agreement with existing experience.

Ammonia Synthesis from First-Principles Calculations

The rate of ammonia synthesis over a nanoparticle ruthenium catalyst can be calculated directly on the basis of a quantum chemical treatment of the problem using density functional theory. We compared the results to measured rates over a ruthenium catalyst supported on magnesium aluminum spinel. When the size distribution of ruthenium particles measured by transmission electron microscopy was used as the link between the catalyst material and the theoretical treatment, the calculated rate was within a factor of 3 to 20 of the experimental rate. This offers hope for computer-based methods in the search for catalysts.

Role of strain and ligand effects in the modification of the electronic and chemical properties of bimetallic surfaces

Periodic density functional calculations are used to illustrate how the combination of strain and ligand effects modify the electronic and surface chemical properties of Ni, Pd, and Pt monolayers supported on other transition metals. Strain and the ligand effects are shown to change the width of the surface d band, which subsequently moves up or down in energy to maintain a constant band filling. Chemical properties such as the dissociative adsorption energy of hydrogen are controlled by changes induced in the average energy of the d band by modification of the d-band width.

Spectroscopic link between adsorption site occupation and local surface chemical reactivity

In this Letter we show that sequences of adsorbate-induced shifts of surface core level (SCL) x-ray photoelectron spectra contain profound information on surface changes of electronic structure and reactivity. Energy shifts and intensity changes of time-lapsed spectral components follow simple rules, from which adsorption sites are directly determined. Theoretical calculations rationalize the results for transition metal surfaces in terms of the energy shift of the d-band center of mass and this proves that adsorbate-induced SCL shifts provide a spectroscopic measure of local surface reactivity.

Modification of the surface electronic and chemical properties of Pt(111) by subsurface 3d transition metals

The modification of the electronic and chemical properties of Pt(111) surfaces by subsurface 3d transition metals was studied using density-functional theory. In each case investigated, the Pt surface d-band was broadened and lowered in energy by interactions with the subsurface 3d metals, resulting in weaker dissociative adsorption energies of hydrogen and oxygen on these surfaces. The magnitude of the decrease in adsorption energy was largest for the early 3d transition metals and smallest for the late 3d transition metals. In some cases, dissociative adsorption was calculated to be endothermic. The surfaces investigated in this study had no lateral strain in them, demonstrating that strain is not a necessary factor in the modification of bimetallic surface properties. The implications of these findings are discussed in the context of catalyst design, particularly for fuel cell electrocatalysts.

Bonding of gold nanoclusters to oxygen vacancies on rutile TiO2(110)

Through an interplay between scanning tunneling microscopy (STM) and density functional theory (DFT) calculations, we show that bridging oxygen vacancies are the active nucleation sites for Au clusters on the rutile TiO2(110) surface. We find that a direct correlation exists between a decrease in density of vacancies and the amount of Au deposited. From the DFT calculations we find that the oxygen vacancy is indeed the strongest Au binding site. We show both experimentally and theoretically that a single oxygen vacancy can bind 3 Au atoms on average. In view of the presented results, a new growth model for the TiO2(110) system involving vacancy-cluster complex diffusion is presented.

Combined electronic structure and evolutionary search approach to materials design

We show that density functional theory calculations have reached an accuracy and speed making it possible to use them in conjunction with an evolutionary algorithm to search for materials with specific properties. The approach is illustrated by finding the most stable four component alloys out of the 192 016 possible fcc and bcc alloys that can be constructed out of 32 different metals. A number of well known and new "super alloys" are identified in this way.

Oxygen vacancies as active sites for water dissociation on rutile TiO(2)(110)

Through an interplay between scanning tunneling microscopy experiments and density functional theory calculations, we determine unambiguously the active surface site responsible for the dissociation of water molecules adsorbed on rutile TiO(2)(110). Oxygen vacancies in the surface layer are shown to dissociate H(2)O through the transfer of one proton to a nearby oxygen atom, forming two hydroxyl groups for every vacancy. The amount of water dissociation is limited by the density of oxygen vacancies present on the clean surface exclusively. The dissociation process sets in as soon as molecular water is able to diffuse to the active site.

One-dimensional metallic edge states in MoS2

By the use of density functional calculations it is shown that the edges of a two-dimensional slab of insulating MoS2 exhibit several metallic states. These edge states can be viewed as one-dimensional conducting wires, and we show that they can be observed directly using scanning tunneling microscopy for single-layer MoS2 nanoparticles grown on a support.

Adsorption-induced step formation

Through an interplay between density functional calculations, Monte Carlo simulations and scanning tunneling microscopy experiments, we show that an intermediate coverage of CO on the Pt(110) surface gives rise to a new rough equilibrium structure with more than 50% step atoms. CO is shown to bind so strongly to low-coordinated Pt atoms that it can break Pt-Pt bonds and spontaneously form steps on the surface. It is argued that adsorption-induced step formation may be a general effect, in particular at high gas pressures and temperatures.