Minimum conditions for accurate modeling of urea production via co-electrolysis.

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

Evaluating Adsorbate-Solvent Interactions: Are Dispersion Corrections Necessary?

Incorporating solvent-adsorbate interactions is paramount in models of aqueous (electro)catalytic reactions. Although a number of techniques exist, they are either highly demanding in computational terms or inaccurate. Microsolvation offers a trade-off between accuracy and computational expenses. Here, we dissect a method to swiftly outline the first solvation shell of species adsorbed on transition-metal surfaces and assess their corresponding solvation energy. Interestingly, dispersion corrections are generally not needed in the model, but caution is to be exercised when water-water and water-adsorbate interactions are of similar magnitude.

Anodic and Cathodic Platinum Dissolution Processes Involve Different Oxide Species.

The degradation of Pt-containing oxygen reduction catalysts for fuel cell applications is strongly linked to the electrochemical surface oxidation and reduction of Pt. Here, we study the surface restructuring and Pt dissolution mechanisms during oxidation/reduction for the case of Pt(100) in 0.1 M HClO4 by combining operando high-energy surface X-ray diffraction, online mass spectrometry, and density functional theory. Our atomic-scale structural studies reveal that anodic dissolution, detected during oxidation, and cathodic dissolution, observed during the subsequent reduction, are linked to two different oxide phases. Anodic dissolution occurs predominantly during nucleation and growth of the first, stripe-like oxide. Cathodic dissolution is linked to a second, amorphous Pt oxide phase that resembles bulk PtO2 and starts to grow when the coverage of the stripe-like oxide saturates. In addition, we find the amount of surface restructuring after an oxidation/reduction cycle to be potential-independent after the stripe-like oxide has reached its saturation coverage.

A general but still unknown characteristic of active oxygen evolution electrocatalysts.

The unsatisfactory electrocatalysis of the oxygen evolution reaction (OER) is a major hurdle for the sustainable production of hydrogen using water electrolyzers. Besides, most state-of-the-art catalysts are based on expensive and scant elements such as Ru and Ir. Hence, it is paramount to establish the features of active OER catalysts to make well-informed searches. Here, an affordable statistical analysis exposes a general yet unnoticed characteristic of active materials for the OER: they frequently have three out of four electrochemical steps with free energies above 1.23 eV. For such catalysts, the first three steps (abbreviated as: H2O → *OH, *OH → *O, *O → *OOH) are statistically prone to be over 1.23 eV, and the second step is often potential limiting. Finally, "electrochemical symmetry", a recently introduced concept, is shown to be a simple and convenient criterion for the in silico design of enhanced OER catalysts, as materials with three steps over 1.23 eV tend to be highly symmetric.

The ABC of Generalized Coordination Numbers and Their Use as a Descriptor in Electrocatalysis.

The quest for enhanced electrocatalysts can be boosted by descriptor-based analyses. Because adsorption energies are the most common descriptors, electrocatalyst design is largely based on brute-force routines that comb materials databases until an energetic criterion is verified. In this review, it is shown that an alternative is provided by generalized coordination numbers (denoted by CN ¯ $\overline {{\rm{CN}}} $ or GCN), an inexpensive geometric descriptor for strained and unstrained transition metals and some alloys. CN ¯ $\overline {{\rm{CN}}} $ captures trends in adsorption energies on both extended surfaces and nanoparticles and is used to elaborate structure-sensitive electrocatalytic activity plots and selectivity maps. Importantly, CN ¯ $\overline {{\rm{CN}}} $ outlines the geometric configuration of the active sites, thereby enabling an atom-by-atom design, which is not possible using energetic descriptors. Specific examples for various adsorbates (e.g., *OH, *OOH, *CO, and *H), metals (e.g., Pt and Cu), and electrocatalytic reactions (e.g., O2  reduction, H2  evolution, CO oxidation, and reduction) are presented, and comparisons are made against other descriptors.

Extracting Features of Active Transition Metal Electrodes for NO Electroreduction with Catalytic Matrices.

Electrocatalytic reduction of oxidized nitrogen compounds (NOx) promises to help rebalance the nitrogen cycle. It is widely accepted that nitrate reduction to NH4+/NH3 involves NO as an intermediate, and NO hydrogenation is the potential-limiting step of NO reduction. Whether *NO hydrogenates to *NHO or *NOH is still a matter of debate, which makes it difficult to optimize catalysts for NOx electroreduction. Here, "catalytic matrices" are used to swiftly extract features of active transition metal catalysts for NO electroreduction. The matrices show that active catalysts statistically stabilize *NHO over *NOH and have undercoordinated sites. Besides, square-symmetry active sites with Cu and other elements may prove active for NO electroreduction. Finally, multivariate regressions are able to reproduce the main features found by the matrices, which opens the door for more sophisticated machine-learning studies. In sum, catalytic matrices may ease the analysis of complex electrocatalytic reactions on multifaceted materials.

Active and durable R2MnRuO7 pyrochlores with low Ru content for acidic oxygen evolution.

The production of green hydrogen in water electrolyzers is limited by the oxygen evolution reaction (OER). State-of-the-art electrocatalysts are based on Ir. Ru electrocatalysts are a suitable alternative provided their performance is improved. Here we show that low-Ru-content pyrochlores (R2MnRuO7, R = Y, Tb and Dy) display high activity and durability for the OER in acidic media. Y2MnRuO7 is the most stable catalyst, displaying 1.5 V at 10 mA cm-2 for 40 h, or 5000 cycles up to 1.7 V. Computational and experimental results show that the high performance is owed to Ru sites embedded in RuMnOx surface layers. A water electrolyser with Y2MnRuO7 (with only 0.2 mgRu cm-2) reaches 1 A cm-2 at 1.75 V, remaining stable at 200 mA cm-2 for more than 24 h. These results encourage further investigation on Ru catalysts in which a partial replacement of Ru by inexpensive cations can enhance the OER performance.

Energetics and Kinetics of Hydrogen Electrosorption on a Graphene-Covered Pt(111) Electrode.

The Angstrom-scale space between graphene and its substrate provides an attractive playground for scientific exploration and can lead to breakthrough applications. Here, we report the energetics and kinetics of hydrogen electrosorption on a graphene-covered Pt(111) electrode using electrochemical experiments, in situ spectroscopy, and density functional theory calculations. The graphene overlayer influences the hydrogen adsorption on Pt(111) by shielding the ions from the interface and weakening the Pt-H bond energy. Analysis of the proton permeation resistance with controlled graphene defect density proves that the domain boundary defects and point defects are the pathways for proton permeation in the graphene layer, in agreement with density functional theory (DFT) calculations of the lowest energy proton permeation pathways. Although graphene blocks the interaction of anions with the Pt(111) surfaces, anions do adsorb near the defects: the rate constant for hydrogen permeation is sensitively dependent on anion identity and concentration.

Computational description of surface hydride phases on Pt(111) electrodes.

Surface platinum hydride structures may exist and play a potentially important role during electrocatalysis and cathodic corrosion of Pt(111). Earlier work on platinum hydrides suggests that Pt may form clusters with multiple equivalents of hydrogen. Here, using thermodynamic methods and density functional theory, we compared several surface hydride structures on Pt(111). The structures contain multiple monolayers of hydrogen in or near the surface Pt layer. The hydrogen in these structures may bind the subsurface or reconstruct the surface both in the set of initial configurations and in the resulting (meta)stable structures. Multilayer stable configurations share one monolayer of subsurface H stacking between the top two Pt layers. The structure containing two monolayers (MLs) of H is formed at -0.29 V vs normal hydrogen electrode, is locally stable with respect to configurations with similar H densities, and binds H neutrally. Structures with 3 and 4 ML H form at -0.36 and -0.44 V, respectively, which correspond reasonably well to the experimental onset potential of cathodic corrosion on Pt(111). For the 3 ML configuration, the top Pt layer is reconstructed by interstitial H atoms to form a well-ordered structure with Pt atoms surrounded by four, five, or six H atoms in roughly square-planar and octahedral coordination patterns. Our work provides insight into the operando surface state during low-potential reduction reactions on Pt(111) and shows a plausible precursor for cathodic corrosion.

Using micro-solvation and generalized coordination numbers to estimate the solvation energies of adsorbed hydroxyl on metal nanoparticles.

Solvent-adsorbate interactions have a great impact on catalytic processes in aqueous systems. Implicit solvent calculations are inexpensive but inaccurate toward hydrogen bonds, while a full incorporation of explicit solvation is computationally demanding. Micro-solvation attempts to break this dilemma by including only those solvent molecules directly interacting with the solute and any nearby interfaces, thereby providing a compromise between accuracy and computational expenses. Here, we show that micro-solvation of *OH and its relation to adsorption sites is largely transferable across late transition metal nanoparticles. Solvation energies for *OH on nanoparticles of Ir, Pd, and Pt range from -0.63 ± 0.04 eV to -0.67 ± 0.12 eV, while those on Au and Ag are -0.75 ± 0.07 eV and -1.01 ± 0.05 eV, respectively. These results enable the use of average solvation corrections for *OH on late transition metal nanostructures.

Mechanistic insight into electrocatalytic glyoxal reduction on copper and its relation to CO2 reduction.

Copper electrodes produce several industrially relevant chemicals and fuels during the electrochemical CO2 reduction reaction (CO2RR). Knowledge about the reaction pathways can help tune the reaction selectivity toward higher-value products. To probe the uncertain role of the C2 molecule glyoxal, we electrochemically reduced it on polycrystalline Cu and quantified its liquid-phase products, namely, ethanol, ethylene glycol, and acetaldehyde. The gas phase contained hydrogen and traces of ethylene. In contrast with previous hypothesis, a one-to-one comparison with CO2RR on Cu indicates that glyoxal is neither a major intermediate in the pathway toward ethylene nor in the pathway toward ethanol. In addition, great possibilities for the selective, low-temperature production of ethylene glycol are open, as computational modelling shows that ethylene glycol and ethanol are produced on different active sites. Thus, apart from the mechanistic insight into CO2RR, this study gives new directions to facilitate the electrification of chemical processes at refineries.

Automated versus Chemically Intuitive Deconvolution of Density Functional Theory (DFT)-Based Gas-Phase Errors in Nitrogen Compounds.

Catalysis models involving metal surfaces and gases are regularly based on density functional theory (DFT) calculations at the generalized gradient approximation (GGA). Such models may have large errors in view of the poor DFT-GGA description of gas-phase molecules with multiple bonds. Here, we analyze three correction schemes for the PBE-calculated Gibbs energies of formation of 13 nitrogen compounds. The first scheme is sequential and based on chemical intuition, the second one is an automated optimization based on chemical bonds, and the third one is an automated optimization that capitalizes on the errors found by the first scheme. The mean and maximum absolute errors are brought down close to chemical accuracy by the third approach by correcting the inaccuracies in the NNO and ONO backbones and those in N-O and N-N bonds. This work shows that chemical intuition and automated optimization can be combined to swiftly enhance the predictiveness of DFT-GGA calculations of gases.

Revealing the Nature of Active Sites on Pt-Gd and Pt-Pr Alloys during the Oxygen Reduction Reaction.

For large-scale applications of hydrogen fuel cells, the sluggish kinetics of the oxygen reduction reaction (ORR) have to be overcome. So far, only platinum (Pt)-group catalysts have shown adequate performance and stability. A well-known approach to increase the efficiency and decrease the Pt loading is to alloy Pt with other metals. Still, for catalyst optimization, the nature of the active sites is crucial. In this work, electrochemical scanning tunneling microscopy (EC-STM) is used to probe the ORR active areas on Pt5Gd and Pt5Pr in acidic media under reaction conditions. The technique detects localized fluctuations in the EC-STM signal, which indicates differences in the local activity. The in situ experiments, supported by coordination-activity plots based on density functional theory calculations, show that the compressed Pt-lanthanide (111) terraces contribute the most to the overall activity. Sites with higher coordination, as found at the bottom of step edges or concavities, remain relatively inactive. Sites of lower coordination, as found near the top of step edges, show higher activity, presumably due to an interplay of strain and steric hindrance effects. These findings should be vital in designing nanostructured Pt-lanthanide electrocatalysts.

Finding Key Factors for Efficient Water and Methanol Activation at Metals, Oxides, MXenes, and Metal/Oxide Interfaces.

Activating water and methanol is crucial in numerous catalytic, electrocatalytic, and photocatalytic reactions. Despite extensive research, the optimal active sites for water/methanol activation are yet to be unequivocally elucidated. Here, we combine transition-state searches and electronic charge analyses on various structurally different materials to identify two features of favorable O-H bond cleavage in H2O, CH3OH, and hydroxyl: (1) low barriers appear when the charge of H moieties remains approximately constant during the dissociation process, as observed on metal oxides, MXenes, and metal/oxide interfaces. Such favorable kinetics is closely related to adsorbate/substrate hydrogen bonding and is enhanced by nearly linear O-H-O angles and short O-H distances. (2) Fast dissociation is observed when the rotation of O-H bonds is facile, which is favored by weak adsorbate binding and effective orbital overlap. Interestingly, we find that the two features are energetically proportional. Finally, we find conspicuous differences between H2O/CH3OH and OH activation, which hints toward the use of carefully engineered interfaces.

Interplaying coordination and ligand effects to break or make adsorption-energy scaling relations.

The linear relations between adsorption energies are one of the cornerstones of contemporary catalysis in view of the simplicity and predictive power of the computational models built upon them. Despite their extensive use, the exact nature of scaling relations is not yet fully understood, and a comprehensive theory of scaling relations is yet to be elaborated. So far, it is known that scalability is dictated by the degree of resemblance of the adsorbed species. In this work, density functional theory calculations show that CO and OH, two dissimilar species, scale or not depending on the surface facet where they adsorb at Pt alloys. This peculiar behavior arises from an interplay of ligand and geometric effects that can be used to modulate adsorption-energy scalability. This study opens new possibilities in catalysis, as it shows that surface coordination is a versatile tool to either balance or unbalance the similarities among adsorbates at alloy surfaces.

Monitoring the active sites for the hydrogen evolution reaction at model carbon surfaces.

Carbon is ubiquitous as an electrode material in electrochemical energy conversion devices. If used as a support material, the evolution of H2 is undesired on carbon. However, recently, carbon-based materials have aroused significant interest as economic and eco-conscious alternatives to noble metal catalysts. The targeted design of improved carbon electrode materials requires atomic scale insight into the structure of the sites that catalyse H2 evolution. This work shows that electrochemical scanning tunnelling microscopy under reaction conditions (n-EC-STM) can be used to monitor the active sites of highly oriented pyrolytic graphite for the hydrogen evolution reaction. With down to atomic resolution, the most active sites in acidic medium are pinpointed near edge sites and defects, whereas the basal planes remain inactive. Density functional theory calculations support these findings and reveal that only specific defects on graphite are active. Motivated by these results, the extensive usage of n-EC-STM on doped carbon-based materials is encouraged to locate their active sites and guide the synthesis of enhanced electrocatalysts.

Selectivity Map for the Late Stages of CO and CO2 Reduction to C2 Species on Copper Electrodes.

The electrochemical CO and CO2 reduction reactions (CORR and CO2 RR) using copper catalysts and renewable electricity hold promise as a carbon-neutral route to produce commodity chemicals and fuels. However, the exact mechanisms and structure sensitivity of Cu electrodes toward C2 products are still under debate. Herein, we investigate ethylene oxide reduction (EOR) as a proxy to the late stages of CORR to ethylene, and the results are compared to those of acetaldehyde reduction to ethanol. Density functional theory (DFT) calculations show that ethylene oxide undergoes ring opening before exclusively reducing to ethylene via *OH formation. Based on generalized coordination numbers (CN), a selectivity map for the late stages of CORR and CO2 RR shows that sites with moderate coordination (5.9 < CN < 7.5) are efficient for ethylene production, with pristine Cu(100) being more active than defective surfaces such as Cu(311). In contrast, kinks and edges are more active for ethanol production, while (111) terraces are relatively inert.

Designing water splitting catalysts using rules of thumb: advantages, dangers and alternatives.

Thermodynamic analysis of the oxygen evolution reaction (OER) hints toward an intrinsic overpotential caused by the nonoptimal adsorption-energy scaling relation between OH and OOH. Consequently, nowadays it is a widely accepted yet unverified rule of thumb that breaking such a scaling relation results in enhanced catalytic activity. In this perspective, we show that breaking the OH-OOH scaling relation does not per se lower the OER overpotential. Instead, electrocatalytic symmetry and ease of optimization are shown to be key factors when screening for enhanced OER catalysts. The essence of electrocatalytic symmetry is captured by a descriptor called the electrochemical-step symmetry index (ESSI). In turn, the ease of optimization and whether it should be scaling-based or scaling-free is provided by a procedure called δ-ε optimization. Finally, taking the search for bifunctional catalysts for oxygen electrocatalysis as an example, we show that the alternative analysis can be straightforwardly extended to other electrocatalytic reactions.

Trends in C-O and N-O bond scission on rutile oxides described using oxygen vacancy formation energies.

Reactivity trends on transition metals can generally be understood through the d-band model, but no analogous theory exists for transition metal oxides. This limits the generality of analyses in oxide-based catalysis and surface chemistry and has motivated the appearance of numerous descriptors. Here we show that oxygen vacancy formation energy (ΔE Vac) is an inexpensive yet accurate and general descriptor for trends in transition-state energies, which are usually difficult to assess. For rutile-type oxides (MO2 with M = 3d metals from Ti to Ni), we show that ΔE Vac captures the trends in C-O and N-O bond scission of CO2, CH3OH, N2O, and NH2OH at oxygen vacancies. The proportionality between ΔE Vac and transition-state energies is rationalized by analyzing the oxygen-metal bonds, which change from ionic to covalent from TiO2 to NiO2. ΔE Vac may be used to design oxide catalysts, in particular those where lattice oxygen and/or oxygen vacancies participate in the catalytic cycles.