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.

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.

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.

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.

Influence of Van der Waals Interactions on the Solvation Energies of Adsorbates at Pt-Based Electrocatalysts.

Solvation can significantly modify the adsorption energy of species at surfaces, thereby influencing the performance of electrocatalysts and liquid-phase catalysts. Thus, it is important to understand adsorbate solvation at the nanoscale. Here we evaluate the effect of van der Waals (vdW) interactions described by different approaches on the solvation energy of *OH adsorbed on near-surface alloys (NSAs) of Pt. Our results show that the studied functionals can be divided into two groups, each with rather similar average *OH solvation energies: (1) PBE and PW91; and (2) vdW functionals, RPBE, PBE-D3 and RPBE-D3. On average, *OH solvation energies are less negative by ∼0.14 eV in group (2) compared to (1), and the values for a given alloy can be extrapolated from one functional to another within the same group. Depending on the desired level of accuracy, these concrete observations and our tabulated values can be used to rapidly incorporate solvation into models for electrocatalysis and liquid-phase catalysis.

Enhanced Electroreduction of Carbon Dioxide to Methanol Using Zinc Dendrites Pulse-Deposited on Silver Foam.

The electrocatalytic CO2 reduction reaction (CO2 RR) can dynamise the carbon cycle by lowering anthropogenic CO2 emissions and sustainably producing valuable fuels and chemical feedstocks. Methanol is arguably the most desirable C1 product of CO2 RR, although it typically forms in negligible amounts. In our search for efficient methanol-producing CO2 RR catalysts, we have engineered Ag-Zn catalysts by pulse-depositing Zn dendrites onto Ag foams (PD-Zn/Ag foam). By themselves, Zn and Ag cannot effectively reduce CO2 to CH3 OH, while their alloys produce CH3 OH with Faradaic efficiencies of approximately 1 %. Interestingly, with nanostructuring PD-Zn/Ag foam reduces CO2 to CH3 OH with Faradaic efficiency and current density values reaching as high as 10.5 % and -2.7 mA cm-2 , respectively. Control experiments and DFT calculations pinpoint strained undercoordinated Zn atoms as the active sites for CO2 RR to CH3 OH in a reaction pathway mediated by adsorbed CO and formaldehyde. Surprisingly, the stability of the *CHO intermediate does not influence the activity.

Structure- and Potential-Dependent Cation Effects on CO Reduction at Copper Single-Crystal Electrodes.

The complexity of the electrocatalytic reduction of CO to CH4 and C2H4 on copper electrodes prevents a straightforward elucidation of the reaction mechanism and the design of new and better catalysts. Although structural and electrolyte effects have been separately studied, there are no reports on structure-sensitive cation effects on the catalyst's selectivity over a wide potential range. Therefore, we investigated CO reduction on Cu(100), Cu(111), and Cu(polycrystalline) electrodes in 0.1 M alkaline hydroxide electrolytes (LiOH, NaOH, KOH, RbOH, CsOH) between 0 and -1.5 V vs RHE. We used online electrochemical mass spectrometry and high-performance liquid chromatography to determine the product distribution as a function of electrode structure, cation size, and applied potential. First, cation effects are potential dependent, as larger cations increase the selectivity of all electrodes toward ethylene at E > -0.45 V vs RHE, but methane is favored at more negative potentials. Second, cation effects are structure-sensitive, as the onset potential for C2H4 formation depends on the electrode structure and cation size, whereas that for CH4 does not. Fourier Transform infrared spectroscopy (FTIR) and density functional theory help to understand how cations favor ethylene over methane at low overpotentials on Cu(100). The rate-determining step to methane and ethylene formation is CO hydrogenation, which is considerably easier in the presence of alkaline cations for a CO dimer compared to a CO monomer. For Li+ and Na+, the stabilization is such that hydrogenated dimers are observable with FTIR at low overpotentials. Thus, potential-dependent, structure-sensitive cation effects help steer the selectivity toward specific products.

Spectroscopic Observation of a Hydrogenated CO Dimer Intermediate During CO Reduction on Cu(100) Electrodes.

Carbon dioxide and carbon monoxide can be electrochemically reduced to useful products such as ethylene and ethanol on copper electrocatalysts. The process is yet to be optimized and the exact mechanism and the corresponding reaction intermediates are under debate or unknown. In particular, it has been hypothesized that the C-C bond formation proceeds via CO dimerization and further hydrogenation. Although computational support for this hypothesis exists, direct experimental evidence has been elusive. In this work, we detect a hydrogenated dimer intermediate (OCCOH) using Fourier transform infrared spectroscopy at low overpotentials in LiOH solutions. Density functional theory calculations support our assignment of the observed vibrational bands. The formation of this intermediate is structure sensitive, as it is observed only during CO reduction on Cu(100) and not on Cu(111), in agreement with previous experimental and computational observations.

Capturing Solvation Effects at a Liquid/Nanoparticle Interface by Ab Initio Molecular Dynamics: Pt201 Immersed in Water.

Solvation can substantially modify the adsorption properties of heterogeneous catalysts. Although essential for achieving realistic theoretical models, assessing such solvent effects over nanoparticles is challenging from a computational standpoint due to the complexity of those liquid/metal interfaces. This effect is investigated by ab initio molecular dynamics simulations at 350 K of a large platinum nanoparticle immersed in liquid water. The first solvation layer contains twice as much physisorbed water molecules above the terraces, than chemisorbed ones located only at edges and corners. The solvent stabilizes the binding energy of chemisorbates: 66% of the total gain comes from interactions with physisorbed molecules and 34% from the influence of bulk liquid.

Double-Stranded Water on Stepped Platinum Surfaces.

The interaction of platinum with water plays a key role in (electro)catalysis. Herein, we describe a combined theoretical and experimental study that resolves the preferred adsorption structure of water wetting the Pt(111)-step type with adjacent (111) terraces. Double stranded lines wet the step edge forming water tetragons with dissimilar hydrogen bonds within and between the lines. Our results qualitatively explain experimental observations of water desorption and impact our thinking of solvation at the Pt electrochemical interface.

Anisotropic etching of rhodium and gold as the onset of nanoparticle formation by cathodic corrosion.

Cathodic corrosion is a phenomenon in which negatively polarized metal electrodes are degraded by cathodic etching and nanoparticle formation. Though these changes are dramatic and sometimes even visible by eye, the exact mechanisms underlying cathodic corrosion are still unclear. This work aims to improve the understanding of cathodic corrosion by studying its onset on rhodium and gold electrodes, which are subjected to various constant cathodic potentials in 10 M NaOH. After this polarization, the electrodes are studied using cyclic voltammetry and scanning electron microscopy, allowing a corrosion onset potential of -1.3 V vs. NHE for rhodium and -1.6 V vs. NHE for gold to be defined. The mildness of the potentials on both metals suggests that cathodic corrosion is less extreme and more ubiquitous than expected. Furthermore, we are able to observe well-defined rectangular etch pits on rhodium. Combined with rhodium cyclic voltammetry, this indicates a strong preference for forming (100) sites during corrosion. In contrast, a (111) preference is indicated on gold by voltammetry and the presence of well-oriented quasi-octahedral nanoparticles. This different etching behavior is suggested to be caused by preferential adsorption of sodium ions to surface defects, as is confirmed by density functional theory calculations.

Initial stages of water solvation of stepped platinum surfaces.

Platinum is an active catalyst for a large number of (electro)chemical reactions in aqueous solution. The observed catalytic activities result from an interplay between the intrinsic adsorption properties of platinum surfaces and their interaction with the aqueous environment. Although water networks have been extensively studied on close-packed surfaces, little is known about high-coverage solvation environments around defects. Here, we report DFT calculations on medium- to high-coverage water adsorption structures near the (100) step edge on Pt(533). We find that isolated ring structures adjacent to step edges form hexagons or pentagons. For higher coverages, 6 possible adsorption structures with varying ring sizes along the step edge and almost identical adsorption energies are observed. From our results we conclude that the favorable interaction of the H-down oriented water molecules, adjacent to the step edge, with the step dipole plays an important role in the formation of these structures. Furthermore, our results explain why water networks on stepped surfaces originate at the step edges, and extend towards the adjacent terraces, in agreement with previous experiments. These results show how step edges act as anchoring points for water adsorption and suggest that solvation of defects might dominate water structures on real platinum surfaces.

Bond-making and breaking between carbon, nitrogen, and oxygen in electrocatalysis.

Many catalytic reactions involving small molecules, which are key transformations in sustainable energy and chemistry, involve the making or breaking of a bond between carbon, nitrogen and oxygen. It has been observed that such heterogeneously (electro)catalyzed reactions often exhibit remarkable and unusual structure sensitivity, in the sense that they take place preferentially on catalyst surfaces with a long-ranged two-dimensional (100) atomic structure. Steps and defects in this two-dimensional structure lower the catalytic activity. Such structure sensitivity must be due to the existence of a special active site on these two-dimensional (100) terraces. Employing detailed density functional theory calculations, we report here the identification of this special active site for a variety of catalytic reactions. The calculations also illustrate how this specific site breaks the well-known rule that under-coordinated surface atoms bind adsorbates stronger, thereby providing the atomic-level explanation for the lack of reactivity of steps and defects for the reactions under consideration. The breakdown of such rule results in significant deviations from commonly observed energetic scaling relations between chemisorbates. Thus, this work provides new design rules for the development of thermodynamically efficient catalysts for an important class of bond-making and bond-breaking reactions.

Density functional theory study of adsorption of H2O, H, O, and OH on stepped platinum surfaces.

We report on density functional theory (DFT)-GGA (generalized gradient approximation) computed adsorption energetics of water and the water-related fragments OH, O, and H on stepped Pt surfaces in the low coverage limit. The Pt(100) step edge as encountered on Pt(533) shows increased binding for all species studied, while the Pt(110) step edge, as found on Pt(553) shows only significantly enhanced binding for O and OH. Comparing these results to ultra high vacuum experiments reveals that DFT can explain the main experimental trends semiquantitatively.

Tailoring the catalytic activity of electrodes with monolayer amounts of foreign metals.

During the past decade, electrocatalysis has attracted significant attention primarily due to the increased interest in the development of new generations of devices for electrochemical energy conversion. This has resulted in a progress in both fundamental understanding of the complex electrocatalytic systems and in the development of efficient synthetic schemes to tailor the surface precisely at the atomic level. One of the viable concepts in electrocatalysis is to optimise the activity through the direct engineering of the properties of the topmost layers of the surface, where the reactions take place, with monolayer and sub-monolayer amounts of metals. This forms (bi)metallic systems where the electronic structure of the active sites is optimised using the interplay between the nature and position of the atoms of solute metals at the surface. In this review, we focus on recent theoretical and experimental achievements in designing efficient (bi)metallic electrocatalysts with selective positioning of foreign atoms to form a variety of active catalytic sites at the electrode surface. We summarize recent results published in the literature and outline challenges for computational and experimental electrocatalysis to engineer active and selective catalysts using atomic layers.

Fast prediction of adsorption properties for platinum nanocatalysts with generalized coordination numbers.

Platinum is a prominent catalyst for a multiplicity of reactions because of its high activity and stability. As Pt nanoparticles are normally used to maximize catalyst utilization and to minimize catalyst loading, it is important to rationalize and predict catalytic activity trends in nanoparticles in simple terms, while being able to compare these trends with those of extended surfaces. The trends in the adsorption energies of small oxygen- and hydrogen-containing adsorbates on Pt nanoparticles of various sizes and on extended surfaces were analyzed through DFT calculations by making use of the generalized coordination numbers of the surface sites. This simple and predictive descriptor links the geometric arrangement of a surface to its adsorption properties. It generates linear adsorption-energy trends, captures finite-size effects, and provides more accurate descriptions than d-band centers and usual coordination numbers. Unlike electronic-structure descriptors, which require knowledge of the densities of states, it is calculated manually. Finally, it was shown that an approximate equivalence exists between generalized coordination numbers and d-band centers.

Error awareness in the volcano plots of oxygen electroreduction to hydrogen peroxide.

Electrocatalysis holds the key to the decentralized production of hydrogen peroxide via the two-electron oxygen reduction reaction (ORR, ). However, cost-effective, active, and selective catalysts are still sought after. While density functional theory (DFT) has already led to the discovery of various enhanced catalysts, it has a severe yet often unnoticed drawback: the ill description of O2 and H2O2. Here, we analyze the impact of the errors in those two species on the most widespread activity plots in the literature, namely free-energy diagrams and Sabatier-type volcano plots. Uncorrected or partially corrected gas-phase energies lead to appreciably different activity plots that may provide inaccurate predictions. Indeed, we show for a variety of electrocatalysts that only when the errors in O2 and H2O2 are corrected can DFT mimic the experiments. In sum, this work provides concrete guidelines to avoid a common pitfall of computational models for electrocatalytic hydrogen peroxide production.

Cation Effects on the Adsorbed Intermediates of CO2 Electroreduction Are Systematic and Predictable.

The electrode-electrolyte interface, and in particular the nature of the cation, has considerable effects on the activity and product selectivity of the electrochemical reduction of CO2. Therefore, to improve the electrocatalysis of this challenging reaction, it is paramount to ascertain whether cation effects on adsorbed intermediates are systematic. Here, DFT calculations are used to show that the effects of K+, Na+, and Mg2+, on single carbon CO2 reduction intermediates can either be stabilizing or destabilizing depending on the metal and the adsorbate. Because systematic trends are observed, cation effects can be accurately predicted in simple terms for a wide variety of metals, cations and adsorbed species. These results are then applied to the reduction of CO2 to CO on four different catalytic surfaces (Au, Ag, Cu, Pd) and activation of weak-binding metals is consistently observed by virtue of the stabilization of the key intermediate *COOH.