Bifunctional alloys for the electroreduction of CO2 and CO.

We use density functional theory to study the reduction of CO2 and CO to hydrocarbons through a formyl pathway on (111) and (211) facets of L12 alloys with an A3B composition. We find that several alloys may reduce the thermodynamic overpotential for CO reduction by more than 0.2 V compared to a copper step, however, these alloys are most often rather unstable in aqueous environment or have low alloy formation energies and may be susceptible to segregation destroying the active sites. Strategies to improve alloy stability against corrosion or segregation would likely be needed in order to realize the full potential of these alloys.

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.

Surface chemistry. Probing the transition state region in catalytic CO oxidation on Ru.

Femtosecond x-ray laser pulses are used to probe the carbon monoxide (CO) oxidation reaction on ruthenium (Ru) initiated by an optical laser pulse. On a time scale of a few hundred femtoseconds, the optical laser pulse excites motions of CO and oxygen (O) on the surface, allowing the reactants to collide, and, with a transient close to a picosecond (ps), new electronic states appear in the O K-edge x-ray absorption spectrum. Density functional theory calculations indicate that these result from changes in the adsorption site and bond formation between CO and O with a distribution of OC-O bond lengths close to the transition state (TS). After 1 ps, 10% of the CO populate the TS region, which is consistent with predictions based on a quantum oscillator model.

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.

Real-time observation of surface bond breaking with an x-ray laser.

We used the Linac Coherent Light Source free-electron x-ray laser to probe the electronic structure of CO molecules as their chemisorption state on Ru(0001) changes upon exciting the substrate by using a femtosecond optical laser pulse. We observed electronic structure changes that are consistent with a weakening of the CO interaction with the substrate but without notable desorption. A large fraction of the molecules (30%) was trapped in a transient precursor state that would precede desorption. We calculated the free energy of the molecule as a function of the desorption reaction coordinate using density functional theory, including van der Waals interactions. Two distinct adsorption wells-chemisorbed and precursor state separated by an entropy barrier-explain the anomalously high prefactors often observed in desorption of molecules from metals.

Theoretical evidence for low kinetic overpotentials in Li-O2 electrochemistry.

We develop a density functional theory model for the electrochemical growth and dissolution of Li(2)O(2) on various facets, terminations, and sites (terrace, steps, and kinks) of a Li(2)O(2) surface. We argue that this is a reasonable model to describe discharge and charge of Li-O(2) batteries over most of the discharge-charge cycle. Because non-stoichiometric surfaces are potential dependent and since the potential varies during discharge and charge, we study the thermodynamic stability of facets, terminations, and steps as a function of potential. This suggests that different facets, terminations, and sites may dominate in charge relative to those for discharge. We find very low thermodynamic overpotentials (<0.2 V) for both discharge and charge at many sites on the facets studied. These low thermodynamic overpotentials for both discharge and charge are in very good agreement with the low kinetic overpotentials observed in recent experiments. However, there are other predicted paths for discharge/charge that have higher overpotentials, so the phase space available for the electrochemistry opens up with overpotential.

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.

Electrical conductivity in Li2O2 and its role in determining capacity limitations in non-aqueous Li-O2 batteries.

Non-aqueous Li-air or Li-O(2) cells show considerable promise as a very high energy density battery couple. Such cells, however, show sudden death at capacities far below their theoretical capacity and this, among other problems, limits their practicality. In this paper, we show that this sudden death arises from limited charge transport through the growing Li(2)O(2) film to the Li(2)O(2)-electrolyte interface, and this limitation defines a critical film thickness, above which it is not possible to support electrochemistry at the Li(2)O(2)-electrolyte interface. We report both electrochemical experiments using a reversible internal redox couple and a first principles metal-insulator-metal charge transport model to probe the electrical conductivity through Li(2)O(2) films produced during Li-O(2) discharge. Both experiment and theory show a "sudden death" in charge transport when film thickness is ~5 to 10 nm. The theoretical model shows that this occurs when the tunneling current through the film can no longer support the electrochemical current. Thus, engineering charge transport through Li(2)O(2) is a serious challenge if Li-O(2) batteries are ever to reach their potential.

Universal transition state scaling relations for (de)hydrogenation over transition metals.

We analyse the transition state energies for 249 hydrogenation/dehydrogenation reactions of atoms and simple molecules over close-packed and stepped surfaces and nanoparticles of transition metals using Density Functional Theory. Linear energy scaling relations are observed for the transition state structures leading to transition state scaling relations for all the investigated reactions. With a suitable choice of reference systems the transition state scaling relations form a universality class that can be approximated with one single linear relation describing the entire range of reactions over all types of surfaces and nanoclusters.

Communications: Elementary oxygen electrode reactions in the aprotic Li-air battery.

We discuss the electrochemical reactions at the oxygen electrode of an aprotic Li-air battery. Using density functional theory to estimate the free energy of intermediates during the discharge and charge of the battery, we introduce a reaction free energy diagram and identify possible origins of the overpotential for both processes. We also address the question of electron conductivity through the Li(2)O(2) electrode and show that in the presence of Li vacancies Li(2)O(2) becomes a conductor.

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.

Modeling ethanol decomposition on transition metals: a combined application of scaling and Brønsted-Evans-Polanyi relations.

Applying density functional theory (DFT) calculations to the rational design of catalysts for complex reaction networks has been an ongoing challenge, primarily because of the high computational cost of these calculations. Certain correlations can be used to reduce the number and complexity of DFT calculations necessary to describe trends in activity and selectivity across metal and alloy surfaces, thus extending the reach of DFT to more complex systems. In this work, the well-known family of Brønsted-Evans-Polanyi (BEP) correlations, connecting minima with maxima in the potential energy surface of elementary steps, in tandem with a scaling relation, connecting binding energies of complex adsorbates with those of simpler ones (e.g., C, O), is used to develop a potential-energy surface for ethanol decomposition on 10 transition metal surfaces. Using a simple kinetic model, the selectivity and activity on a subset of these surfaces are calculated. Experiments on supported catalysts verify that this simple model is reasonably accurate in describing reactivity trends across metals, suggesting that the combination of BEP and scaling relations may substantially reduce the cost of DFT calculations required for identifying reactivity descriptors of more complex reactions.

Towards the computational design of solid catalysts.

Over the past decade the theoretical description of surface reactions has undergone a radical development. Advances in density functional theory mean it is now possible to describe catalytic reactions at surfaces with the detail and accuracy required for computational results to compare favourably with experiments. Theoretical methods can be used to describe surface chemical reactions in detail and to understand variations in catalytic activity from one catalyst to another. Here, we review the first steps towards using computational methods to design new catalysts. Examples include screening for catalysts with increased activity and catalysts with improved selectivity. We discuss how, in the future, such methods may be used to engineer the electronic structure of the active surface by changing its composition and structure.

Alloys of platinum and early transition metals as oxygen reduction electrocatalysts.

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.

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.

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.