Trends in electrocatalytic activity of the oxygen evolution reaction (OER) are investigated on the basis of a large database of HO* and HOO* adsorption energies on oxide surfaces. The theoretical overpotential was calculated by applying standard density functional theory in combination with the computational standard hydrogen electrode (SHE) model. We showed that by the discovery of a universal scaling relation between the adsorption energies of HOO* vs HO*, it is possible to analyze the reaction free energy diagrams of all the oxides in a general way. This gave rise to an activity volcano that was the same for a wide variety of oxide catalyst materials and a universal descriptor for the oxygen evolution activity, which suggests a fundamental limitation on the maximum oxygen evolution activity of planar oxide catalysts.

The widespread use of low-temperature polymer electrolyte membrane fuel cells for mobile applications will require significant reductions in the amount of expensive Pt contained within their cathodes, which drive the oxygen reduction reaction (ORR). Although progress has been made in this respect, further reductions through the development of more active and stable electrocatalysts are still necessary. Here we describe a new set of ORR electrocatalysts consisting of Pd or Pt alloyed with early transition metals such as Sc or Y. They were identified using density functional theory calculations as being the most stable Pt- and Pd-based binary alloys with ORR activity likely to be better than Pt. Electrochemical measurements show that the activity of polycrystalline Pt(3)Sc and Pt(3)Y electrodes is enhanced relative to pure Pt by a factor of 1.5-1.8 and 6-10, respectively, in the range 0.9-0.87 V.

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.

We develop a model based on density functional theory calculations to describe trends in catalytic activity for CO2 electroreduction to CO in terms of the adsorption energy of the reaction intermediates, CO and COOH. The model is applied to metal surfaces as well as the active site in the CODH enzymes and shows that the strong scaling between adsorbed CO and adsorbed COOH on metal surfaces is responsible for the persistent overpotential. The active site of the CODH enzyme is not subject to these scaling relations and optimizes the relative binding energies of these adsorbates, allowing for an essentially reversible process with a low overpotential.

Based on density functional theory calculations we investigate the electrochemically most stable surface structures as a function of pH and electrostatic potential for Pt(111), Ag(111) and Ni(111), and we construct surface Pourbaix diagrams. We study the oxygen reduction reaction (ORR) on the different surface structures and calculate the free energy of the intermediates. We estimate their catalytic activity for ORR by determining the highest potential at which all ORR reaction steps reduce the free energy. We obtain self-consistency in the sense that the surface is stable under the potential at which that particular surface can perform ORR. Using the self consistent surfaces, the activity of the very reactive Ni surface changes dramatically, whereas the activity of the more noble catalysts Pt and Ag remains unchanged. The reason for this difference is the oxidation of the reactive surface. Oxygen absorbed on the surface shifts the reactivity towards the weak binding region, which in turn increases the activity. The oxidation state of the surface and the ORR potential are constant versus the reversible hydrogen electrode (RHE). The dissolution potential in acidic solution, on the other hand, is constant vs. the standard hydrogen electrode (SHE). For Ag, this means that where the potential for dissolution and ORR are about the same at pH = 0, Ag becomes more stable relative to RHE as pH is increased. Hence the pH dependent stability offers an explanation for the possible use of Ag in alkaline fuel cell cathodes.

In this work, we extend the activity volcano for oxygen reduction from the face-centered cubic (fcc) metal (111) facet to the (100) facet. Using density functional theory calculations, we show that the recent findings of constant scaling between OOH* and OH* holds on the fcc metal (100) facet, as well. Using this fact, we show the existence of a universal activity volcano to describe oxygen reduction electrocatalysis with a minimum overpotential, η min = 0.37 ± 0.1 V. Specifically, we find that the (100) facet of Pt is found to bind oxygen intermediates too strongly and is not active for oxygen reduction reaction (ORR). In contrast, Au( 100) is predicted to be more active than Au(111) and comparable in activity to Pt alloys. Using this activity volcano, we further predict that Au alloys that bind OH more strongly could display improved ORR activity on the (100) facet. We carry out a computational search over candidate alloys and suggest that alloying Au with early transition metals could lead to materials that exhibit enhanced ORR activity. We also prescribe a simple principle in the choice of nanoparticle catalysts in which the chosen catalyst should have its dense facet on the weak binding leg of the volcano so that its more open facets and defects could also be active.

We present a theoretical analysis of trends in overpotentials for electrocatalytic CO2 reduction based on density functional theory calculations. The analysis is based on understanding variations in the free energy of intermediates and mapping out the potential at which different elementary steps are exergonic as a measure of the catalytic activity. We study different surface structures and introduce a simple model for including the effect of adsorbate-adsorbate interactions. We find that high coverages of CO under typical reaction conditions for the more reactive transition metals affect the catalytic activity towards the CO2 reduction reaction, but the ordering of metal activities is not changed. For the hydrogen evolution reaction, a high CO coverage shifts the maximum activity towards more reactive metals than Pt.

Based on density functional theory (DFT) calculations we study the electrochemical chlorine evolution reaction on rutile (110) oxide surfaces. First we construct the Pourbaix surface diagram for IrO(2) and RuO(2), and from this we find the chlorine evolution reaction intermediates and identify the lowest overpotential at which all elementary reaction steps in the chlorine evolution reaction are downhill in free energy. This condition is then used as a measure for catalytic activity. Linear scaling relations between the binding energies of the intermediates and the oxygen binding energies at cus-sites are established for MO(2) (M being Ir, Ru, Pt, Ti). The linear relations form the basis for constructing a generalized surface phase diagram where two parameters, the potential and the binding energy of oxygen, are needed to determine the surface composition. We calculate the catalytic activity as function of the oxygen binding energy, giving rise to a Sabatier volcano. By combining the surface phase diagram and the volcano describing the catalytic activity, we find that the reaction mechanism differs depending on catalyst material. The flexibility in reaction path means that the chlorine evolution activity is high for a wide range of oxygen binding energies. We find that the required overpotential for chlorine evolution is lower than the overpotential necessary for oxygen evolution.