Cation-induced changes in the inner- and outer-sphere mechanisms of electrocatalytic CO2 reduction.

The underlying mechanism of cation effects on CO2RR remains debated. Herein, we study cation effects by simulating both outer-sphere electron transfer (OS-ET) and inner-sphere electron transfer (IS-ET) pathways during CO2RR via constrained density functional theory molecular dynamics (cDFT-MD) and slow-growth DFT-MD (SG-DFT-MD), respectively. Our results show without any cations, only OS-ET is feasible with a barrier of 1.21 eV. In the presence of K+ (Li+), OS-ET shows a very high barrier of 2.93 eV (4.15 eV) thus being prohibited. However, cations promote CO2 activation through IS-ET with the barrier of only 0.61 eV (K+) and 0.91 eV (Li+), generating the key intermediate (adsorbed CO[Formula: see text]). Without cations, CO2-to-CO[Formula: see text](ads) conversion cannot proceed. Our findings reveal cation effects arise from short-range Coulomb interactions with reaction intermediates. These results disclose that cations modulate the inner- and outer-sphere pathways of CO2RR, offering substantial insights on the cation specificity in the initial CO2RR steps.

Structure and energetics of liquid water-hydroxyl layers on Pt(111).

The interactions between liquid water and hydroxyl species on Pt(111) surfaces have been intensely investigated due to their importance to fuel cell electrocatalysis. Here we present a molecular dynamics study of their structure and energetics using an ensemble of neural network potentials, which allow us to obtain unprecedented statistical sampling. We first study the energetics of hydroxyl formation, where we find a near-linear adsorption energy profile, which exhibits a soft and gradual increase in the differential adsorption energy at high hydroxyl coverages. This is strikingly different from the predictions of the conventional bilayer model, which displays a kink at 1/3ML OH coverage indicating a sizeable jump in differential adsorption energy, but within the statistical uncertainty of previously reported ab initio molecular dynamics studies. We then analyze the structure of the interface, where we provide evidence for the water-OH/Pt(111) interface being hydrophobic at high hydroxyl coverages. We furthermore explain the observed adsorption energetics by analyzing the hydrogen bonding in the water-hydroxyl adlayers, where we argue that the increase in differential adsorption energy at high OH coverage can be explained by a reduction in the number of hydrogen bonds from the adsorbed water molecules to the hydroxyls.

Improving the Activity of M-N4 Catalysts for the Oxygen Reduction Reaction by Electrolyte Adsorption.

Metal and nitrogen codoped carbons (M-N/Cs) have emerged as promising alternatives to platinum-based catalysts for the oxygen reduction reaction (ORR). DFT calculations are used to investigate the adsorption of anions and impurities from the electrolyte on the active site, modeled as an M-N4 motif embedded in a planar carbon sheet (M=Cr, Mn, Fe, Co). The two-dimensional catalyst structure implies that each metal atom has two potential active sites, one on each side of the sheet. Adsorption of anions or impurities on both sites results in poisoning, but adsorption on one of the sites leads to a modified ORR activity on the remaining site. The calculated adsorption energies show that a number of species adsorb only on one of the two sites under realistic experimental conditions. Furthermore, a few of these adsorbates modify the adsorption energies of the ORR intermediates on the remaining site, in such a way that the limiting potential is improved.

Combined DFT and Differential Electrochemical Mass Spectrometry Investigation of the Effect of Dopants in Secondary Zinc-Air Batteries.

Zinc-air batteries offer the potential of low-cost energy storage with high specific energy, but at present secondary Zn-air batteries suffer from poor cyclability. To develop economically viable secondary Zn-air batteries, several properties need to be improved: choking of the cathode, catalyzing the oxygen evolution and reduction reactions, limiting dendrite formation and suppressing the hydrogen evolution reaction (HER). Understanding and alleviating HER at the negative electrode in a secondary Zn-air battery is a substantial challenge, for which it is necessary to combine computational and experimental research. Here, we combine differential electrochemical mass spectrometry (DEMS) and density functional theory (DFT) calculations to investigate the fundamental role and stability when cycling in the presence of selected beneficial additives, that is, In and Bi, and Ag as a potentially unfavorable additive. We show that both In and Bi have the desired property for a secondary battery, that is, upon recharging they will remain on the surface, thereby retaining the beneficial effects on Zn dissolution and suppression of HER. This is confirmed by DEMS, where it is observed that In reduces HER and Bi affects the discharge potential beneficially compared to a battery without additives. Using a simple procedure based on adsorption energies calculated with DFT, it is found that Ag suppresses OH adsorption, but, unlike In and Bi, it does not hinder HER. Finally, it is shown that mixing In and Bi is beneficial compared to the additives by themselves as it improves the electrochemical performance and cyclic stability of the secondary Zn-air battery.

Giant onsite electronic entropy enhances the performance of ceria for water splitting.

Previous studies have shown that a large solid-state entropy of reduction increases the thermodynamic efficiency of metal oxides, such as ceria, for two-step thermochemical water splitting cycles. In this context, the configurational entropy arising from oxygen off-stoichiometry in the oxide, has been the focus of most previous work. Here we report a different source of entropy, the onsite electronic configurational entropy, arising from coupling between orbital and spin angular momenta in lanthanide f orbitals. We find that onsite electronic configurational entropy is sizable in all lanthanides, and reaches a maximum value of ≈4.7 k B per oxygen vacancy for Ce4+/Ce3+ reduction. This unique and large positive entropy source in ceria explains its excellent performance for high-temperature catalytic redox reactions such as water splitting. Our calculations also show that terbium dioxide has a high electronic entropy and thus could also be a potential candidate for solar thermochemical reactions.Solid-state entropy of reduction increases the thermodynamic efficiency of ceria for two-step thermochemical water splitting. Here, the authors report a large and different source of entropy, the onsite electronic configurational entropy arising from coupling between orbital and spin angular momenta in f orbitals.

Descriptors and Thermodynamic Limitations of Electrocatalytic Carbon Dioxide Reduction on Rutile Oxide Surfaces.

A detailed understanding of the electrochemical reduction of CO2 into liquid fuels on rutile metal oxide surfaces is developed by using DFT calculations. We consider oxide overlayer structures on RuO2 (1 1 0) surfaces as model catalysts to elucidate the trends and limitations in the CO2 reduction reaction (CO2RR) based on thermodynamic analysis. We aim to specify the requirements for CO2RR catalysts to establish adsorbate scaling relations and use these to derive activity volcanoes. Computational results show that the OH* binding free energy is a good descriptor of the thermodynamic limitations and it defines the left leg of the activity volcano for CO2RR. HCOOH* is a key intermediate for products formed through further reduction, for example, methanediol, methanol, and methane. The surfaces that do not bind HCOOH* are selective towards formic acid (HCOOH) production, but hydrogen evolution limits their suitability. We determine the ideal binding free energy for H* and OH* to facilitate selective CO2RR over H2 /CO evolution to be ΔGB [H]>0.5 eV and -0.5 eV<ΔGB [OH]<0.1 eV. The Re-containing overlayers considered in this work display excellent promise for selectivity, although they are active at a highly reducing potential.

Decoupling strain and ligand effects in ternary nanoparticles for improved ORR electrocatalysis.

Ternary Pt-Au-M (M = 3d transition metal) nanoparticles show reduced OH adsorption energies and improved activity for the oxygen reduction reaction (ORR) compared to pure Pt nanoparticles, as obtained by density functional theory. The strain and ligand effects in nanoparticles are decoupled and correlated with the extended Pt(111) surface for benchmarking. The ternary metal in the core allows for tuning the catalytic activity through strain effects. Pt-Au-M for M = Cr, Mn, Co, Cu, Zn nanoparticles are of particular interest as they exhibit an optimal contribution of strain, ligand effects and stability. Good agreement is found with experimental studies showing increased activity of Pt-Au-Fe/Ni nanoparticles, and mid to late 3d transition metals are predicted to exhibit enhance activity and stability with respect to pure Pt nanoparticles.

Selective Electrochemical Generation of Hydrogen Peroxide from Water Oxidation.

Water is a life-giving source, fundamental to human existence, yet over a billion people lack access to clean drinking water. The present techniques for water treatment such as piped, treated water rely on time and resource intensive centralized solutions. In this work, we propose a decentralized device concept that can utilize sunlight to split water into hydrogen and hydrogen peroxide. The hydrogen peroxide can oxidize organics while the hydrogen bubbles out. In enabling this device, we require an electrocatalyst that can oxidize water while suppressing the thermodynamically favored oxygen evolution and form hydrogen peroxide. Using density functional theory calculations, we show that the free energy of adsorbed OH* can be used to determine selectivity trends between the 2e(-) water oxidation to H2O2 and the 4e(-) oxidation to O2. We show that materials which bind oxygen intermediates sufficiently weakly, such as SnO2, can activate hydrogen peroxide evolution. We present a rational design principle for the selectivity in electrochemical water oxidation and identify new material candidates that could perform H2O2 evolution selectively.

A DFT-based genetic algorithm search for AuCu nanoalloy electrocatalysts for CO2 reduction.

Using a DFT-based genetic algorithm (GA) approach, we have determined the most stable structure and stoichiometry of a 309-atom icosahedral AuCu nanoalloy, for potential use as an electrocatalyst for CO2 reduction. The identified core-shell nano-particle consists of a copper core interspersed with gold atoms having only copper neighbors and a gold surface with a few copper atoms in the terraces. We also present an adsorbate-dependent correction scheme, which enables an accurate determination of adsorption energies using a computationally fast, localized LCAO-basis set. These show that it is possible to use the LCAO mode to obtain a realistic estimate of the molecular chemisorption energy for systems where the computation in normal grid mode is not computationally feasible. These corrections are employed when calculating adsorption energies on the Cu, Au and most stable mixed particles. This shows that the mixed Cu135@Au174 core-shell nanoalloy has a similar adsorption energy, for the most favorable site, as a pure gold nano-particle. Cu, however, has the effect of stabilizing the icosahedral structure because Au particles are easily distorted when adding adsorbates.

Elucidating the activity of stepped Pt single crystals for oxygen reduction.

The unexpectedly high measured activity of Pt[n(111) × (111)] and Pt[n(111) × (100)] stepped single crystal surfaces towards the oxygen reduction reaction (ORR) is explained utilizing the hydroxyl binding energy as the activity descriptor. Using this descriptor (estimated using experimental data obtained by different groups), a well-defined Sabatier-type volcano is observed for the activities measured for the Pt[n(111) × (111)] and Pt[n(111) × (100)] stepped single crystals, in remarkable agreement with earlier theoretical studies. We propose that the observed destabilisation of *OH species at these surfaces is due to the decreased solvation of the adsorbed hydroxyl intermediates on adjacent terrace sites.

Trends in electrochemical CO2 reduction activity for open and close-packed metal surfaces.

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.

Modeling CO2 reduction on Pt(111).

Density functional theory was used to model the electrochemical reduction of CO2 on Pt(111) with an explicit solvation layer and the presence of extra hydrogen atoms to represent a negatively charged electrode. We focused on the electronic energy barriers for the first four lowest energy proton-electron transfer steps for reducing CO2 on Pt(111) beginning with adsorbed *CO2 and continuing with *COOH, *CO + H2O, *COH, and ending with *C + H2O. We find that simple elementary steps in which a proton is transferred to an adsorbate (such as the protonation of *CO to *COH) have small barriers on the order of 0.1 eV. Elementary steps in which a proton is transferred and a C-O bond is simultaneously cleaved show barriers on the order of 0.5 eV. All barriers calculated for these steps show no sign of being insurmountable at room temperature. To explain why these barriers are so small, we analyze the charge density and the density of states plots to see that first, the electron transfer is decoupled from the proton transfer so that in the initial state, the surface and adsorbate are already charged up and can easily accept the proton from solution. Also, we see that in the cases where barriers are on the order of 0.1 eV, electron density in the initial state localizes on the oxygen end of the adsorbate, while electron density is more spread out on the surface for initial states of the C-O bond cleaving elementary steps.

First principles investigation of zinc-anode dissolution in zinc-air batteries.

With surging interest in high energy density batteries, much attention has recently been devoted to metal-air batteries. The zinc-air battery has been known for more than a hundred years and is commercially available as a primary battery, but recharging has remained elusive, in part because the fundamental mechanisms still remain to be fully understood. Here, we present a density functional theory investigation of the zinc dissolution (oxidation) on the anode side in the zinc-air battery. Two models are envisaged, the most stable (0001) surface and a kink surface. The kink model proves to be more accurate as it brings about some important features of bulk dissolution and yields results in good agreement with experiments. From the adsorption energies of hydroxyl species and experimental values, we construct a free energy diagram and confirm that there is a small overpotential associated with the reaction. The applied methodology provides new insight into computational modelling and design of secondary metal-air batteries.

Understanding Trends in the Electrocatalytic Activity of Metals and Enzymes for CO2 Reduction to CO.

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.

Solar hydrogen production with semiconductor metal oxides: new directions in experiment and theory.

An overview of a collaborative experimental and theoretical effort toward efficient hydrogen production via photoelectrochemical splitting of water into di-hydrogen and di-oxygen is presented here. We present state-of-the-art experimental studies using hematite and TiO(2) functionalized with gold nanoparticles as photoanode materials, and theoretical studies on electro and photo-catalysis of water on a range of metal oxide semiconductor materials, including recently developed implementation of self-interaction corrected energy functionals.

The Pt(111)/electrolyte interface under oxygen reduction reaction conditions: an electrochemical impedance spectroscopy study.

The Pt(111)/electrolyte interface has been characterized during the oxygen reduction reaction (ORR) in 0.1 M HClO(4) using electrochemical impedance spectroscopy. The surface was studied within the potential region where adsorption of OH* and O* species occur without significant place exchange between the adsorbate and Pt surface atoms (0.45-1.15 V vs RHE). An equivalent electric circuit is proposed to model the Pt(111)/electrolyte interface under ORR conditions within the selected potential window. This equivalent circuit reflects three processes with different time constants, which occur simultaneously during the ORR at Pt(111). Density functional theory (DFT) calculations were used to correlate and interpret the results of the measurements. The calculations indicate that the coadsorption of ClO(4)* and Cl* with OH* is unlikely. Our analysis suggests that the two-dimensional (2D) structures formed in O(2)-free solution are also formed under ORR conditions.

Electrochemical chlorine evolution at rutile oxide (110) surfaces.

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.

Surface Pourbaix diagrams and oxygen reduction activity of Pt, Ag and Ni(111) surfaces studied by DFT.

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.