Modeling Potential-Dependent Electrochemical Activation Barriers: Revisiting the Alkaline Hydrogen Evolution Reaction

Emerging Atomistic Modeling Methods for Heterogeneous Electrocatalysis

Heterogeneous electrocatalysis lies at the center of various technologies that could help enable a sustainable future. However, its complexity makes it challenging to accurately and efficiently model at an atomic level. Here, we review emerging atomistic methods to simulate the electrocatalytic interface with special attention devoted to the components/effects that have been challenging to model, such as solvation, electrolyte ions, electrode potential, reaction kinetics, and pH. Additionally, we review relevant computational spectroscopy methods. Then, we showcase several examples of applying these methods to understand and design catalysts relevant to green hydrogen. We also offer experimental views on how to bridge the gap between theory and experiments. Finally, we provide some perspectives on opportunities to advance the field.

Dual-strategy engineered nickel phosphide for achieving efficient hydrazine-assisted hydrogen production in seawater

Electrocatalytic hydrogen production in seawater to alleviate freshwater shortage pressures is promising, but is hindered by the sluggish oxygen evolution reaction and detrimental chloride electrochemistry. Herein, a dual strategy approach of Fe-doping and CeO2-decoration in nickel phosphide (Fe-Ni2P/CeO2) is rationally designed to achieve superior bifunctional catalytic performance for the hydrogen evolution reaction (HER) and hydrazine oxidation reaction (HzOR) in seawater. Notably, the two-electrode Fe-Ni2P/CeO2-based hybrid seawater electrolyzer realizes energy-efficient and chlorine-free hydrogen production with ultralow cell voltages of 0.051 and 0.597 V at 10 and 400 mA cm-2, which are significantly lower than those needed in the hydrazine-free seawater electrolyzer. Density functional theory calculations manifest that the combination of Fe doping and heterointerface construction between Fe-Ni2P and CeO2 can adjust the electronic structure of the Ni2P and optimize the water dissociation barrier and hydrogen adsorption free energy, leading to improvement of the intrinsic catalytic performance. This route affords a feasible solution for future large-scale hydrogen generation using abundant ocean water.

Bubble-water/catalyst triphase interface microenvironment accelerates photocatalytic OER via optimizing semi-hydrophobic OH radical

Photocatalytic water splitting (PWS) as the holy grail reaction for solar-to-chemical energy conversion is challenged by sluggish oxygen evolution reaction (OER) at water/catalyst interface. Experimental evidence interestingly shows that temperature can significantly accelerate OER, but the atomic-level mechanism remains elusive in both experiment and theory. In contrast to the traditional Arrhenius-type temperature dependence, we quantitatively prove for the first time that the temperature-induced interface microenvironment variation, particularly the formation of bubble-water/TiO2(110) triphase interface, has a drastic influence on optimizing the OER kinetics. We demonstrate that liquid-vapor coexistence state creates a disordered and loose hydrogen-bond network while preserving the proton transfer channel, which greatly facilitates the formation of semi-hydrophobic •OH radical and O-O coupling, thereby accelerating OER. Furthermore, we propose that adding a hydrophobic substance onto TiO2(110) can manipulate the local microenvironment to enhance OER without additional thermal energy input. This result could open new possibilities for PWS catalyst design.

Designing Efficient Nitrate Reduction Electrocatalysts by Identifying and Optimizing Active Sites of Co-Based Spinels

Cobalt-based spinel oxides (i.e., Co3O4) are emerging as low-cost and selective electrocatalysts for the electrochemical nitrate reduction reaction (NO3-RR) to ammonia (NH3), although their activity is still unsatisfactory and the genuine active site is unclear. Here, we discover that the NO3-RR activity of Co3O4 is highly dependent on the geometric location of the Co site, and the NO3-RR prefers to occur at octahedral Co (CoOh) rather than tetrahedral Co (CoTd) sites. Moreover, CoOhO6 is electrochemically transformed to CoOhO5 along with the formation of O vacancies (Ov) during the process of NO3-RR. Both experimental and theoretic results reveal that in situ generated CoOhO5-Ov configuration is the genuine active site for the NO3-RR. To further enhance the activity of CoOh sites, we replace inert CoTd with different contents of Cu2+ cations, and a volcano-shape correlation between NO3-RR activity and electronic structures of CoOh is observed. Impressively, in 1.0 M KOH, (Cu0.6Co0.4)Co2O4 with optimized CoOh sites achieves a maximum NH3 Faradaic efficiency of 96.5% with an ultrahigh NH3 rate of 1.09 mmol h-1 cm-2 at -0.45 V vs reversible hydrogen electrode, outperforming most of other reported nonprecious metal-based electrocatalysts. Clearly, this work paves new pathways for boosting the NO3-RR activity of Co-based spinels by tuning local electronic structures of CoOh sites.

The electrocatalytic N2 reduction activity of core-shell iron nanoalloy catalysts: a density functional theory (DFT) study

A molecular level understanding of the property evolution in binary nanoalloy catalysts is crucial for designing novel electrocatalysts for ammonia synthesis. In this regard, designing core-shell catalyst structures has been a versatile approach to achieve the product selectivity. Herein, we investigated the activity evolution of Fe-based core-shell (M15@Fe50) (M = Co, Ni, or Cu) clusters for the nitrogen reduction reaction (NRR). Nitrogen reduction following the associative mechanistic pathway is significantly activated over the Cu15@Fe50 cluster. The d-band center from the electronic structure analysis is found to be upshifted, justifying the activity towards the NRR. The reduction reaction occurs via the surface restructuring of the catalyst, in which the *NH2 formation is found to be the lowest endergonic potential determining step compared to pristine Fe(110). Based on this, the high NRR activity of the Cu15@Fe50 cluster has been proposed, which, we envision, will provide useful insights into the position and compositional effects of core-shell structures for the discovery of efficient NRR electrocatalysts.

Tailoring a local acid-like microenvironment for efficient neutral hydrogen evolution

Electrochemical hydrogen evolution reaction in neutral media is listed as the most difficult challenges of energy catalysis due to the sluggish kinetics. Herein, the Ir-H_xWO₃ catalyst is readily synthesized and exhibits enhanced performance for neutral hydrogen evolution reaction. H_xWO₃ support is functioned as proton sponge to create a local acid-like microenvironment around Ir metal sites by spontaneous injection of protons to WO₃, as evidenced by spectroscopy and electrochemical analysis. Rationalize revitalized lattice-hydrogen species located in the interface are coupled with H_ad atoms on metallic Ir surfaces via thermodynamically favorable Volmer-Tafel steps, and thereby a fast kinetics. Elaborated Ir-H_xWO₃ demonstrates acid-like activity with a low overpotential of 20 mV at 10 mA cm^-2 and low Tafel slope of 28 mV dec^-1, which are even comparable to those in acidic environment. The concept exemplified in this work offer the possibilities for tailoring local reaction microenvironment to regulate catalytic activity and pathway.

Interface Engineering on Amorphous/Crystalline Hydroxides/Sulfides Heterostructure Nanoarrays for Enhanced Solar Water Splitting

Developing highly efficient and stable noble-metal-free electrocatalysts for water splitting is critical for producing clean and sustainable energy. Here, we design a hierarchical transition metal hydroxide/sulfide (NiFe(OH)_x-Ni₃S₂/NF) electrode with dual heterointerface coexistence using a cation exchange-induced surface reconfiguration strategy. The electrode exhibits superior electrocatalytic activities, achieving low overpotentials of 55 mV for hydrogen evolution and 182 mV for oxygen evolution at 10 mA cm^-2. Furthermore, the assembled two-electrode system requires voltages as low as 1.55 and 1.62 V to deliver industrially relevant current densities of 500 and 1000 mA cm^-2, respectively, with excellent durability for over 200 h, which is comparable to commercial electrolysis. Theoretical calculations reveal that the hierarchical heterostructure increases the electronic delocalization of the Fe and Ni catalytic centers, lowering the energy barrier of the rate-limiting step and promoting O₂ desorption. Finally, by implementing the catalysts in a solar-driven water electrolysis system, we demonstrate a record and durable solar-to-hydrogen (STH) conversion efficiency of up to 20.05%. This work provides a promising strategy for developing low-cost and high-efficiency bifunctional catalysts for a large-scale solar-to-hydrogen generation.

Activation of H2O Tailored by Interfacial Electronic States at a Nanoscale Interface for Enhanced Electrocatalytic Hydrogen Evolution

Despite the fundamental and practical significance of the hydrogen evolution reaction (HER), the reaction kinetics at the molecular level are not well-understood, especially in basic media. Here, with ZIF-67-derived Co-based carbon frameworks (Co/NCs) as model catalysts, we systematically investigated the effects of different reaction parameters on the HER kinetics and discovered that the HER activity was directly dependent not on the type of nitrogen in the carbon framework but on the relative content of surface hydroxyl and water (OH-/H2O) adsorbed on Co active sites embedded in carbon frameworks. When the ratio of the OH-/H2O was close to 1:1, the Co/NC nanocatalyst showed the best reaction performance under the condition of high-pH electrolytes, e.g., an overpotential of only 232 mV at a current density of 10 mA cm-2 in the 1 M KOH electrolyte. We unambiguously identified that the structural water molecules (SWs) in the form of hydrous hydroxyl complexes absorbed on metal centers {OHad·H2O@M+} were catalytic active sites for the enhanced HER, where M+ could be transition or alkaline metal cations. Different from the traditional hydrogen bonding of water, the hydroxyl (hydroxide) groups and water molecules in the SWs were mainly bonded together via the spatial interaction between the p orbitals of O atoms, exhibiting features of a delocalized π-bond with a metastable state. These newly formed surface bonds or transitory states could be new weak interactions that synergistically promote both interfacial electron transfer and the activation of water (dissociation of O-H bonds) at the electrode surface, i.e., the formation of activated H adducts (H*). The capture of new surface states not only explains pH-, cation-, and transition-metal-dependent hydrogen evolution kinetics but also provides completely new insights into the understanding of other electrocatalytic reductions involving other small molecules, including CO2, CO, and N2.

Reconciling the Volcano Trend with the Butler-Volmer Model for the Hydrogen Evolution Reaction

The volcano trend has been widely utilized to forecast new optimum catalysts in computational chemistry while the Butler-Volmer relationship is the norm to explain current-potential characteristics from cyclic voltammetry in analytical chemistry. Herein, we develop an electrochemical model for hydrogen evolution reaction exchange currents that reconciles device-level chemistry, atomic-level volcano trend, and the Butler-Volmer relation. We show that the model is a function of the easy-to-compute hydrogen adsorption energy invariably obtained from first-principles atomic simulations. In addition, the model reproduces with high fidelity the experimental exchange currents for elemental metal catalysts over 15 orders of magnitude and is consistent with the recently proposed analytical model based on a data-driven approach. Our findings based on fundamental electrochemistry principles are general and can be applied to other reactions including CO2 reduction, metal oxidation, and lithium (de)intercalation reactions.

Modeling Operando Electrochemical CO2 Reduction

Since the seminal works on the application of density functional theory and the computational hydrogen electrode to electrochemical CO₂ reduction (eCO₂R) and hydrogen evolution (HER), the modeling of both reactions has quickly evolved for the last two decades. Formulation of thermodynamic and kinetic linear scaling relationships for key intermediates on crystalline materials have led to the definition of activity volcano plots, overpotential diagrams, and full exploitation of these theoretical outcomes at laboratory scale. However, recent studies hint at the role of morphological changes and short-lived intermediates in ruling the catalytic performance under operating conditions, further raising the bar for the modeling of electrocatalytic systems. Here, we highlight some novel methodological approaches employed to address eCO₂R and HER reactions. Moving from the atomic scale to the bulk electrolyte, we first show how ab initio and machine learning methodologies can partially reproduce surface reconstruction under operation, thus identifying active sites and reaction mechanisms if coupled with microkinetic modeling. Later, we introduce the potential of density functional theory and machine learning to interpret data from Operando spectroelectrochemical techniques, such as Raman spectroscopy and extended X-ray absorption fine structure characterization. Next, we review the role of electrolyte and mass transport effects. Finally, we suggest further challenges for computational modeling in the near future as well as our perspective on the directions to follow.

Fast-Decoding Algorithm for Electrode Processes at Electrified Interfaces by Mean-Field Kinetic Model and Bayesian Data Assimilation: An Active-Data-Mining Approach for the Efficient Search and Discovery of Electrocatalysts

The microscopic origins of the activity and selectivity of electrocatalysts has been a long-lasting enigma since the 19th century. By applying an active-data-mining approach, employing a mean-field kinetic model and a statistical approach of Bayesian data assimilation, we demonstrate here a fast decoding to extract key properties in the kinetics of complicated electrode processes from current-potential profiles in experimental and literary data. As the proof-of-concept, kinetic parameters on the four-electron oxygen reduction reaction in the 0.1 M HClO4 solution (ORR: O2 + 4e- + 4H+ → 2H2O) of various platinum-based single-crystal electrocatalysts are extracted from our own experiments and third-party literature to investigate the microscopic electrode processes. Furthermore, data assimilation of the mean-field ORR model and experimental data is performed based on Bayesian inference for the inductive estimation of kinetic parameters, which sheds light on the dynamic behavior of kinetic parameters with respect to overpotential. This work shows that a fast-decoding algorithm based on a mean-field kinetic model and Bayesian data assimilation is a promising data-driven approach to extract key microscopic features of complicated electrode processes and therefore will be an important method toward building up advanced human-machine collaborations for the efficient search and discovery of high-performance electrochemical materials.

High-Throughput Screening of Efficient Biatom Catalysts Based on Monolayer Carbon Nitride for the Nitric Oxide Reduction Reaction

Exploring efficient catalysts for the nitric oxide reduction reaction (NORR) toward NH₃ synthesis is becoming increasingly important for tackling both NH₃ synthesis and NO removal problems. Currently, only a few NORR catalysts have been proposed, which are exclusively concentrated on bulk metals or single-atom catalysts. Here, taking monolayer C₂N as an example, we explore the potential of biatom catalysts (BACs) for direct NO-to-NH₃ conversion by means of high-throughput first-principles calculations. According to a rational five-step screening strategy, a promising BAC of Cr₂-C₂N is successfully screened out, exhibiting high stability, activity, and selectivity and a low kinetic barrier for the NORR toward NH₃ synthesis. Importantly, the adsorption energy of N atoms (ΔE_*N) and the Gibbs free energy of NO adsorption (ΔG_*NO) are identified as effective descriptors for efficient NORR catalysts. In addition, through tuning the NO coverage, the NORR on Cr₂-C₂N could produce different products of NH₃ and N₂O, providing the possibility to realize controllable multiproduct BACs. These findings not only suggest the great potential of BACs for direct NO-to-NH₃ conversion but also help in rationally designing high-performance BACs.

Electrolyte pH-dependent hydrogen binding energies and coverages on platinum, iridium, rhodium, and ruthenium surfaces

The weakened Hatop binding strength and increased Hatop coverage are universal phenomena on Pt, Ir, Rh, and Ru surfaces from acidic to alkaline media, which are important factors in the pH-dependent hydrogen reaction kinetics.