Search

Oxygen reduction reaction (ORR) corresponds to the limiting process in fuel cells to convert gaseous hydrogen and oxygen from the air into electricity. Most commonly, electronic structure calculations in the density functional theory (DFT) approximation are applied to comprehend the elementary reaction steps on the atomic scale, including the identification of limiting reaction steps and approximation of electrocatalytic activity. A major challenge in the modeling of the ORR refers to the complexity of the mechanistic pathways given that a plethora of different descriptions, ranging from the mononuclear to OOH dissociation mechanisms, dissociative and oxide pathways, have been reported in the literature. In the present work, the adsorption free energy of the *OH intermediate, ΔG₁, is introduced as an effective descriptor for rapid screening of ORR mechanisms. By applying a rigorous thermodynamic analysis in conjunction with a descriptor-based approach, it is demonstrated how the calculation of a single binding energy in terms of ΔG₁ can be used to identify the energetically favored mechanistic description for any ORR catalyst by analyzing the location of the material in the volcano plot. The introduced procedure may support future DFT-based studies in that the preferred pathway is not missed when modeling the four-electron ORR by a binding-energy approach.

Electrochemical nitrogen reduction reaction (NRR) is mainly hampered by two facts. On the one hand, the hydrogen evolution reaction competes with the NRR under cathodic reaction conditions. On the other hand, the sluggish reaction kinetics of the NRR, due to the transfer of six proton-electron pairs in the reaction mechanism, causes large overpotentials, and thus results in humble intrinsic activity for ammonia formation. Optimization strategies for NRR electrocatalysts are ultimately called for to overcome these issues. While breaking scaling relation is considered as a universal remedy to enhance turnover for any electrocatalytic process, recently, the concept of catalytic resonance theory has been put forth as a second option to obtain reaction rates beyond the limiting Sabatier volcano. Yet, a comparison of these two concepts is missing but urgently needed to comprehend design principles for electrocatalyst optimization on the atomic scale. In the present work, breaking scaling relation and catalytic resonance by programmable catalysis are compared on the example of the NRR over transition-metal oxides. It is demonstrated that a fine tuning rather than a full breaking of nitrogen-containing scaling relations is required to enhance electrocatalytic activity of electrocatalysts at the volcano apex. The success of catalytic resonance theory is based on the occurrence of negative scaling correlations whereas in the case of positive scaling correlations, the impact of this approach on the optimization of electrocatalysts strongly depends on the location of the material in the volcano plot. The obtained insight may spur the development of NRR electrocatalysts with enhanced intrinsic activity by following the outlined design principles. in press

Metal-air batteries are encountered as a promising solution for energy storage due to their high energy density, cost effectiveness, and environmental benefits. Yet, the application of metal-air batteries in practice is still not mature, which is also related to the bifunctional oxygen electrocatalysis at the cathode, comprising the oxygen reduction (ORR) and oxygen evolution (OER) reactions during discharge and charge of the battery, respectively. Experimentally, the performance of electrocatalysts in the OER and ORR is described by bifunctional index (BI), but, so far, there is no direct approach to capture the BI on the atomic scale. Herein, we present a method to ascertain the BI from ab initio theory, thereby combining a data-driven methodology with thermodynamic considerations and microkinetic modeling as a function of the applied overpotential. Our approach allows deriving the BI from simple adsorption free energies, which are easily accessible to electronic structure theory in the density functional theory (DFT) approximation. We outline how our methodology may steer the design of efficient bifunctional catalysts on the atomic scale.

cited By 13

Electrochemical water splitting is a key technology for moving toward a promising energy scenario based on renewable (regenerative) energy resources in that wind and solar energy can be stored and ...

The computational hydrogen electrode (CHE) approach has spurred ab initio investigations in the field of electrocatalysis, since the underlying concept enables to quantify free energy changes, ΔG (...

Understanding the complex interactions of different building blocks within a sophisticated drug-delivery system (DDS), aimed at targeted transport of the drug to malignant cells, requires modeling techniques on different time and length scales. On the example of the anthracycline antibiotic doxorubicin (DOX), we investigate a potential DDS component, consisting of a gold nanoparticle and a short peptide sequence as carriers of DOX. The combination of atomistic molecular dynamics simulations and density functional theory calculations facilitates compiling a volcano plot, which allows deriving general conclusions on DDS constituents for chemotherapeutic agents within the class of anthracycline antibiotics: the nanoparticle and peptide carrier moieties need to be chosen in such a way that the anthracycline body of the drug is able to intercalate between both entities or between two (π-stacking) residues of the peptide. Using the popular volcano framework as a guideline, the present article connects the catalysis and biosimulation communities, thereby identifying a strategy to overcome the limiting volcano relation by tuning the coordination number of the drug in the DDS component.

Hydrogen electrocatalysis has been spurred by theoretical predictions, using simple ab initio thermodynamic considerations, in that the free-binding energy of adsorbed hydrogen has been applied in a heuristic fashion to search for sustainable electrocatalysts as a replacement for scarce platinum in electrolyzers and fuel cells. The original volcano model of Norskov et al. is given in [14] purports that the optimum hydrogen-evolution catalyst binds adsorbed hydrogen thermoneutrally at zero overpotential, a paradigm based on pure thermodynamic considerations. Recently, the Sabatier principle was revisited by factoring the applied overpotential and kinetics into the analysis. The extended Sabatier principle suggests that the optimum hydrogen-evolution catalyst binds adsorbed hydrogen weakly rather than thermoneutrally. This notion is corroborated by the fact that the most active hydrogen-evolution catalysts, Pt, MoS2, or Mo2C, indeed bind hydrogen weakly by about (100–200) meV rather than thermoneutrally at zero overpotential.

The oxygen evolution reaction (OER) is often designated as the enigma in water electrolysis because the development of active and stable OER catalysts is a challenging and formidable task. While ab initio theory in the density functional theory approximation initially focused on the mechanistic description via the OH, O, and OOH adsorbates, in recent years the lattice oxygen evolution reaction (LOER) mechanism attracted increasing attention, given that the LOER is seen as the main reason for catalyst instability under anodic potential conditions. The present concept article critically analyzes the LOER and indicates pitfalls in the interpretation of this mechanistic pathway. A method to assess the energetics of the LOER in relation to conventional OER mechanisms by the compilation of free-energy diagrams is introduced, which may contribute to enhance our understanding of the competing LOER and OER on the atomic scale. Further works are urgently needed to comprehend the interrelationship for the evolution of gaseous oxygen from the electrolyte or the crystal lattice.

Extended Tafel plots at various temperatures for an electrocatalyzed reaction and (possibly) its reversed reaction on single-crystalline model electrodes allow for constructing the (essential part of the) free energy surface, in particular the free energies of the transition states (TS). Free energies of the reaction intermediates (RIs) including the chemical nature of active surface sites (S) are hardly accessible to experiment and need therefore to be taken from constrained ab initio thermodynamics calculations. The compact compilation of experimental kinetic data in the form of a free energy diagram enables a critical assessment and validation of theoretical free energy landscapes based on first-principles kinetics. For three prototypical electrocatalyzed reactions, namely the chlorine evolution reaction (CER) and oxygen evolution reaction (OER) over RuO2(110) as well as hydrogen evolution reaction (HER) on Pt(111), we exemplify this universal approach and discuss potential benefits for theoretical mod...

Ultrathin single-crystalline RuO 2(110) films supported on Ru(0001) are employed as model electrodes to extract kinetic information about the industrially important chlorine evolution reaction (CER) in a 5M concentrated NaCl solution under well-defined electrochemical conditions and variable temperatures. A combination of chronoamperometry (CA) and online electrochemical mass spectrometry (OLEMS) experiments provides insight into the selectivity issue: At pH = 0.9, the CER dominates over oxygen evolution, whereas at pH = 3.5, oxygen evolution and other parasitic side reactions contribute mostly to the total current density. From temperature-dependent CA data for pH = 0.9, we determine the apparent free activation energy of the CER over RuO 2(110) to be 0.91 eV, which compares reasonably well with the theoretical value of 0.79 eV derived from first-principles microkinetics. The experimentally determined apparent free activation energy of 0.91 eV is considered as a benchmark for assessing future improved theoretical modeling from first principles.

Volcano plots are a powerful tool to screen electrode materials in the catalysis and battery science communities. Commonly, simple binding energies are analyzed by the concept of linear scaling relationships to describe activity trends in a homologous series of materials, putting forward the picture that an optimum electrode material in the hydrogen evolution reaction (HER) binds the reaction intermediate (RI) thermoneutrally at zero overpotential. This approach, however, consists of various oversimplifications since the applied overpotential and kinetics are not accounted for in the evaluation. In the present article, the apex of the HER volcano is modeled by microkinetics. It is demonstrated that the volcano's top shifts to weak bonding of the RI with increasing driving force as soon as kinetic effects are factored in the analysis. This paradigm change is corroborated by the fact that the constructed volcano plots, using microkinetics and scaling relations for the apex and legs of the volcano respectively, reproduce the high activities of Pt in the HER and RuO2 in the chlorine evolution reaction.

Volcano analyses have been established as a standard tool in the field of electrocatalysis for assessing the performance of electrodes in a class of materials. The apex of the volcano curve, where the most active electrocatalysts are situated, is commonly defined by a hypothetical ideal material that binds its reaction intermediates thermoneutrally at zero overpotential, in accordance with Sabatier's principle. However, recent studies report a right shift of the apex in a volcano curve, in which the most active electrocatalysts bind their reaction intermediates endergonically rather than thermoneutrally at zero overpotential. Focusing on two-electron process, this Viewpoint addresses the question of how the definition of an optimum catalyst needs to be modified with respect to the requirements of Sabatier's principle when kinetic effects and the applied overpotential are included in the analysis.

In this perspective we discuss how an intimate interaction of experiments with theory is able to deepen our insight into the catalytic reaction system on the molecular level. This strategy is illustrated by discussing various examples from our own research of surface chemistry and model catalysis. The particular examples were carefully chosen to balance the specific strength of both approaches – theory and experiment – and emphasize the benefit of this combined approach. We start with the determination of complex surface structures, where diffraction techniques in combination with theory are clear-cut. The promoter action of alkali metals in heterogeneous catalysis is rationalized with theory and experiment for the case of CO coadsorption. Predictive power of theory is limited as demonstrated with the apparent activity of chlorinated TiO2(110) in the oxidation of HCl: Even if we know all elementary reaction steps of a catalytic reaction mechanism, the overall kinetics may remain elusive and require the application kinetic Monte Carlo simulations. Catalysts are not always stable under reaction conditions and may chemically transform as discussed for the CO oxidation reaction over ruthenium. Under oxidizing reaction conditions ruthenium transforms into RuO2, a process which is hardly understood on the molecular level. Lastly we focus on electrochemical reactions. Here theory is clearly ahead since spectroscopic methods are not available to resolve the processes at the electrode surface.

Multielectron processes in electrochemistry require the stabilization of reaction intermediates (RI) at the electrode surface after every elementary reaction step. Accordingly, the bond strengths of these intermediates are important for assessing the catalytic performance of an electrode material. Current understanding of microscopic processes in modern electrocatalysis research is largely driven by theory, mostly based on ab initio thermodynamics considerations, where stable reaction intermediates at the electrode surface are identified, while the actual free energy barriers (or activation barriers) are ignored. This simple approach is popular in electrochemistry in that the researcher has a simple tool at hand in successfully searching for promising electrode materials. The ab initio TD approach allows for a rough but fast screening of the parameter space with low computational cost. However, ab initio thermodynamics is also frequently employed (often, even based on a single binding energy only) to comprehend on the activity and on the mechanism of an electrochemical reaction. The basic idea is that the activation barrier of an endergonic reaction step consists of a thermodynamic part and an additional kinetically determined barrier. Assuming that the activation barrier scales with thermodynamics (so-called Brønsted-Polanyi-Evans (BEP) relation) and the kinetic part of the barrier is small, ab initio thermodynamics may provide molecular insights into the electrochemical reaction kinetics. However, for many electrocatalytic reactions, these tacit assumptions are violated so that ab initio thermodynamics will lead to contradictions with both experimental data and ab initio kinetics. In this Account, we will discuss several electrochemical key reactions, including chlorine evolution (CER), oxygen evolution reaction (OER), and oxygen reduction (ORR), where ab initio kinetics data are available in order to critically compare the results with those derived from a simple ab initio thermodynamics treatment. We show that ab initio thermodynamics leads to erroneous conclusions about kinetic and mechanistic aspects for the CER over RuO

State-of-the-art material screening in the field of electrocatalysis mainly uses the concept of linear scaling relationships in order to express the (free) adsorption energies of different reaction intermediates, adsorbed on the surface of a solid-state electrocatalyst, as function of a descriptor. This thermodynamic analysis, based on the application of the computational hydrogen electrode approach (CHE), ultimately results in the construction of a Volcano plot, which facilitates identifying promising catalysts within a class of materials. The conventional ab initio Volcano concept, however, lacks of two critical aspects: on the one hand the applied overpotential, which constitutes the driving force of an electrocatalytic reaction, is not included in the underlying approach, since the thermodynamic analysis refers to the standard equilibrium potential of the electrocatalytic process; on the other hand, the kinetics is not accounted for. Herein, an alternate material-screening concept is presented, which promotes a discussion of the catalytic performance within a class of materials by explicitly including both the kinetic description and applied overpotential: kinetic scaling relations enable resolving the rate-determining reaction step in a homologous series of single-crystalline electrocatalysts in the overpotential regime of interest for practical applications. The proposed methodology is exemplified by the chlorine evolution reaction over transition-metal oxides, which corresponds to the anode reaction in the industrially relevant chlor-alkali process for the production of gaseous chlorine as basic chemical.

In the industrially important Chlor-Alkali process, the chlorine evolution reaction (CER) over a ruthenium dioxide (RuO2) catalyst competes with the oxygen evolution reaction (OER). This selectivity issue is elucidated on the microscopic level with the single-crystalline model electrode RuO2(110) by employing density functional theory (DFT) calculations in combination with the concept of volcano plots. We demonstrate that one monolayer of TiO2(110) supported on RuO2(110) enhances the selectivity towards the CER by several orders of magnitudes, while preserving the high activity for the CER. This win-win situation is attributed to the different slopes of the volcano curves for the CER and OER.

Constrained ab initio thermodynamics in the form of a Pourbaix diagram can greatly assist kinetic modeling of a particular electrochemical reaction such as the chlorine evolution reaction (CER) over RuO 2 (110). Pourbaix diagrams reveal stable surface structures, as a function of pH and the potential. The present DFT study indicates that the Pourbaix diagram in the CER potential region above 1.36 V and pH values around zero is dominated by a stable surface structure in which all coordinatively undercoordinated Ru sites (Ru cus ) are capped by on-top oxygen (O ot ). This oxygen saturated RuO 2 (110) surface is considered to serve as the catalytically active phase in the CER, quite in contrast to the heterogeneously catalyzed HCl oxidation (Deacon process), for which the active RuO 2 (110) surface is mainly covered by on-top chlorine. The active sites in the CER are suggested to be Ru cus O ot surface complexes, while in the Deacon process both undercoordinated surface Ru and oxygen sites must be available for the activation of HCl molecules.

Current progress in modern electrocatalysis research is spurred by theory, frequently based on ab initio thermodynamics, where the stable reaction intermediates at the electrode surface are identified, while the actual energy barriers are ignored. This approach is popular in that a simple tool is available for searching for promising electrode materials. However, thermodynamics alone may be misleading to assess the catalytic activity of an electrochemical reaction as we exemplify with the chlorine evolution reaction (CER) over a RuO2 (110) model electrode. The full procedure is introduced, starting from the stable reaction intermediates, computing the energy barriers, and finally performing microkinetic simulations, all performed under the influence of the solvent and the electrode potential. Full kinetics from first-principles allows the rate-determining step in the CER to be identified and the experimentally observed change in the Tafel slope to be explained.

Five and six electrons delocalized in the N,N,N,N plane are characteristic of the highly persistent 4N/5e radical anions (deep green) and σ-bis-homoaromatic 4N/6e dianions (red), respectively, which were generated by one- and two-electron reduction of rigid, proximate bis-diazenes (see reaction below). The longest wavelength UV/Vis absorptions of the radical anion and dianion are strongly dependent on the counterion.

In a specifically designed proximate and almost perfectly synperiplanar diazene–diazene oxide substrate, metathesis is the exclusive photoreaction and occurs with retention of optical purity, providing support for the [π2 + π2]photocycloaddition pathway (tetrazetidine oxide intermediate).