Insights in the Oxygen Reduction Reaction: From Metallic Electrocatalysts to Diporphyrins

Tuning the product selectivity of single-atom catalysts for CO2 reduction beyond CO formation by orbital engineering

Electrochemical CO2 reduction (CO2R) is one of the promising strategies for developing sustainable energy resources. Single-atom catalysts (SACs) have emerged as efficient catalysts for CO2R. However, the efficiency of SACs for the formation of reduction products beyond two-step CO formation is low due to the lower binding strength of the CO intermediate. In this study, we present an orbital engineering strategy based on density functional theory calculations and the fragment molecular orbital approach to tune product selectivity for the CO2R reaction on macrocycle based molecular catalysts (porphyrin and phthalocyanine) and extended SACs (graphene and covalent organic frameworks) with Fe, Co, and Ni dopants. The introduction of neutral axial ligands such as imidazole, pyridine, and trimethyl phosphine to the metal dopants enhances the binding affinity of the CO intermediate. The stability of the catalysts is investigated through the thermodynamic binding energy of the axial ligands and ab initio molecular dynamics simulations (AIMD). The grand canonical potential method is used to determine the reaction free energy values. Using a unified activity volcano plot based on the reaction free energy values, we investigated the catalytic activity and product selectivity at an applied potential of -0.8 V vs. SHE and a pH of 6.8. We found that with the imidazole and pyridine axial ligands, the selectivity of Fe-doped SACs towards the formation of the methanol product is improved. The activity volcano plot for these SACs shows a similar activity to that of the Cu (211) surface. The catalytic activity is found to be directly proportional to the sigma-donating ability of the axial ligands.

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

What Is the Rate-Limiting Step of Oxygen Reduction Reaction on Fe-N-C Catalysts?

Oxygen reduction reaction (ORR) is essential to various renewable energy technologies. An important catalyst for ORR is single iron atoms embedded in nitrogen-doped graphene (Fe-N-C). However, the rate-limiting step of the ORR on Fe-N-C is unknown, significantly impeding understanding and improvement. Here, we report the activation energies of all of the steps, calculated by ab initio molecular dynamics simulations under constant electrode potential. In contrast to the common belief that a hydrogenation step limits the reaction rate, we find that the rate-limiting step is oxygen molecule replacing adsorbed water on Fe. This occurs through concerted motion of H2O desorption and O2 adsorption, without leaving the site bare. Interestingly, despite being an apparent "thermal" process that is often considered to be potential-independent, the barrier reduces with the electrode potential. This can be explained by stronger Fe-O2 binding and weaker Fe-H2O binding at a lower potential, due to O2 gaining electrons and H2O donating electrons to the catalyst. Our study offers new insights into the ORR on Fe-N-C and highlights the importance of kinetic studies in heterogeneous electrochemistry.

Activity And Stability of Single- And Di-Atom Catalysts for the O2 Reduction Reaction

Efficient and inexpensive catalysts for the O2 reduction reaction (ORR) are needed for the advancement of renewable energy technologies. In this study, we designed a computational catalyst-screening method to identify single and di-atom metal dopants from first-row transition elements supported on defect-containing nitrogenated graphene surfaces for the ORR. Based on formation-energy calculations and micro-kinetic modelling of reaction pathways using intermediate binding free energies, we have identified four potentially interesting single-atom catalysts (SACs) and fifteen di-atom catalysts (DACs) with relatively high estimated catalytic activity at 0.8 V vs RHE. Among the best SACs, MnNC shows high stability in both acidic and alkaline media according to our model. For the DACs, we found four possible candidates, MnMn, FeFe, CoCo, and MnNi doped on quad-atom vacancy sites having considerable stability over a wide pH range. The remaining SACs and DACs with high activity are either less stable or show a stability region at an alkaline pH.

Atomic metal coordinated to nitrogen-doped carbon electrocatalysts for proton exchange membrane fuel cells: a perspective on progress, pitfalls and prospectives

Proton exchange membrane fuel cells require reduced construction costs to improve commercial viability, which can be fueled by elimination of platinum as the O2 reduction electrocatalyst. The past 10 years has seen significant developments in synthesis, characterisation, and electrocatalytic performance of the most promising alternative electrocatalyst; single metal atoms coordinated to nitrogen-doped carbon (M-N-C). In this Perspective we recap some of the important achievements of M-N-Cs in the last decade, as well as discussing current knowledge gaps and future research directions for the community. We provide a new outlook on M-N-C stability and atomistic understanding with a set of original density functional theory simulations.

Origin of the superior oxygen reduction activity of zirconium nitride in alkaline media

The anion exchange membrane fuel cell (AEMFC), which can operate in alkaline media, paves a promising avenue for the broad application of earth-abundant element based catalysts. Recent pioneering studies found that zirconium nitride (ZrN) with low upfront capital cost can exhibit high activity, even surpassing that of Pt in alkaline oxygen reduction reaction (ORR). However, the origin of its superior ORR activity was not well understood. Herein, we propose a new theoretical framework to uncover the ORR mechanism of ZrN by integrating surface state analysis, electric field effect simulations, and pH-dependent microkinetic modelling. The ZrN surface was found to be covered by ∼1 monolayer (ML) HO* under ORR operating conditions, which can accommodate the adsorbates in a bridge-site configuration for the ORR. Electric field effect simulations demonstrate that O* adsorption on a 1 ML HO* covered surface only induces a consistently small dipole moment change, resulting in a moderate bonding strength that can account for the superior activity. Based on the identified surface state of ZrN and electric field simulations, pH-dependent microkinetic modelling found that ZrN reaches the Sabatier optimum of the kinetic ORR volcano model in alkaline media, with the simulated polarization curves being in excellent agreement with the experimental data of ZrN and Pt/C. Finally, we show that this theoretical framework can lead to a good explanation for the alkaline oxygen electrocatalysis of other transition metal nitrites such as Fe3N, TiN, and HfN. In summary, this study proposes a new framework to rationalize and design transition metal nitrides for alkaline ORR.

Theoretical Study on B-doped FeN4 Catalyst for Potential-Dependent Oxygen Reduction Reaction

Electrochemical reactions mostly take place at a constant potential, but traditional DFT calculations operate at a neutral charge state. In order to really model experimental conditions, we developed a fixed-potential simulation framework via the iterated optimization and self-consistence of the required Fermi level. The B-doped graphene-based FeN4 sites for oxygen reduction reaction were chosen as the model to evaluate the accuracy of the fixed-potential simulation. The results demonstrate that *OH hydrogenation gets facile while O2 adsorption or hydrogenation becomes thermodynamically unfavorable due to the lower d-band center of Fe atoms in the constant potential state than the neutral charge state. The onset potential of ORR over B-doped FeN4 by performing potential-dependent simulations agree well with experimental findings. This work indicates that the fixed-potential simulation can provide a reasonable and accurate description on electrochemical reactions.

Surface Curvature Effect on Dual-Atom Site Oxygen Electrocatalysis

Improved oxygen electrocatalysis is crucial for the ever-growing energy demand. Metal-nitrogen-carbon (M-N-C) materials are promising candidates for catalysts. Their activity is tunable via varying electronic and geometric properties, such as porosity. Because of the difficulty in modeling porosity, M-N-Cs with variable surface curvature remained largely unexplored. In this work, we developed a realistic in-pore dual-atom site M-N-C model and applied density functional theory to investigate the surface curvature effect on oxygen reduction and evolution reactions. We show that surface curving tailors both scaling relations and energy barriers. Thus, we predict that adjusting the surface curvature can improve the catalytic activity toward mono- and bifunctional oxygen electrocatalysis.

Surface states of dual-atom catalysts should be considered for analysis of electrocatalytic activity

Experimentally well-characterized dual-atom catalysts (DACs), where two adjacent metal atoms are stably anchored on carbon defects, have shown some clear advantages in electrocatalysis compared to conventional catalysts and emerging single-atom catalysts. However, most previous theoretical studies directly used a pristine dual-atom site to analyze the electrocatalytic activity of a DAC. Herein, by analyzing 8 homonuclear and 64 heteronuclear DACs structures with ab initio calculations, our derived surface Pourbaix diagrams show that the surface states of DACs generally differ from a pristine surface at electrocatalytic operating conditions. This phenomenon suggests that the surface state of a DAC should be considered before analyzing the catalytic activity in electrocatalysis, while the electrochemistry-driven pre-adsorbed molecules generated from the liquid phase may either change the electronic properties or even block the active site of DACs. Based on these results, we provide a critical comment to the catalyst community: before analyzing the electrocatalytic activity of a DAC, its surface state should be analyzed beforehand.

Limitations of Electrochemical Nitrogen Oxidation toward Nitrate

The electrocatalytic N₂ oxidation reaction (NOR) using renewable electricity is a promising alternative to the industrial synthesis of nitrate from NH₃ oxidation. However, breaking the triple bond in the nitrogen molecule is one of the most essential challenges in chemistry. In this work, we use density functional theory simulations to investigate the plausible reaction mechanisms of electrocatalytic NOR and its competition with oxygen evolution reaction (OER) at the atomic scale. We focus on the electrochemical conversion of inert N₂ to active *NO during NOR. We propose formation of *N₂O from *N₂ and *O as the rate-determining step (RDS). Following the RDS, a microkinetic model is utilized to study the rate of NOR on metal oxides. Our results demonstrate that a lower activation energy is obtained when a catalyst binds *O weakly. We show that the reaction is extremely challenging but also that design strategies have been suggested to promote electrochemical NOR.

Rare-Earth Single-Atom Catalysts: A New Frontier in Photo/Electrocatalysis

Single-atom catalysts (SACs) provide well-defined active sites with 100% atom utilization, and can be prepared using a wide range of support materials. Therefore, they are attracting global attention, especially in the fields of energy conversion and storage. To date, research has focused on transition-metal and precious-metal-based SACs. More recently, rare-earth (RE)-based SACs have emerged as a new frontier in photo/electrocatalysis owing to their unique electronic structure arising from the spin-orbit coupling of the 4f and valence orbitals, unsaturated coordination environment, and unique behavior as charge-transport bridges. However, a systematic review on the role of the RE active sites, catalytic mechanisms, and synthetic methods for RE SACs is lacking. Therefore, in this review, the latest developments in RE SACs having applications in photo/electrocatalysis are summarized and discussed. First, the theoretical advantages of RE SACs for photo/electrocatalysis are briefly introduced, focusing on the roles of the 4f orbitals and coupled energy levels. In addition, the most recent research progress on RE SACs is summarized for several important photo/electrocatalytic reactions and the corresponding catalytic mechanisms are discussed. Further, the synthetic strategies for the production of RE SACs are reported. Finally, challenges for the development of RE SACs are highlighted, along with future research directions and perspectives.

Computational Design of Gas Sensors Based on V3S4 Monolayer

Novel magnetic gas sensors are characterized by extremely high efficiency and low energy consumption, therefore, a search for a two-dimensional material suitable for room temperature magnetic gas sensors is a critical task for modern materials scientists. Here, we computationally discovered a novel ultrathin two-dimensional antiferromagnet V3S4, which, in addition to stability and remarkable electronic properties, demonstrates a great potential to be applied in magnetic gas sensing devices. Quantum-mechanical calculations within the DFT + U approach show the antiferromagnetic ground state of V3S4, which exhibits semiconducting electronic properties with a band gap of 0.36 eV. A study of electronic and magnetic response to the adsorption of various gas agents showed pronounced changes in properties with respect to the adsorption of NH3, NO2, O2, and NO molecules on the surface. The calculated energies of adsorption of these molecules were -1.25, -0.91, -0.59, and -0.93 eV, respectively. Obtained results showed the prospective for V3S4 to be used as effective sensing materials to detect NO2 and NO, for their capture, and for catalytic applications in which it is required to lower the dissociation energy of O2, for example, in oxygen reduction reactions. The sensing and reducing of NO2 and NO have great importance for improving environmental protection and sustainable development.

Study on the Structure-Activity Relationship Between Single-Atom, Cluster and Nanoparticle Catalysts in a Hierarchical Structure for the Oxygen Reduction Reaction

Literature reports have shown that in primary structures, single-atom catalysts exhibit better performance than cluster and nanoparticles due to their maximum atom utilization and the fine-tuning of the electronic structure of the active sites. Hierarchical structures have recently been extensively studied because of increased active sites and orderliness of channels significantly improves the catalytic performance compare to primary structures especially in nanoparticles, however, the different sized effect of catalysts research has not been reported. Herein, a unique hollow double-shell structure (a distinct cavity-containing) is used as a hierarchical model to study the possible difference between single atom, cluster, and nanoparticle and to establish the corresponding structure-activity relationship. Three Co catalysts are prepared: single atoms (Co-Catalyst-1), clusters (Co-Catalyst-2, 0.5-1 nm), and nanoparticles (Co-Catalyst-3, ≈5 nm) and their oxygen-reduction capacity is evaluated. The unique electronic interactions, the strong electron-withdrawing ability of N in Co-N₄ (Co-Catalyst-1), attract electrons from the electrode to Co, specifically by expediting the generation and transformation of the rate-determining step intermediates *OOH. The variant spatial structure which is caused by Co atom aggregation, and led to surface area, pore size, and carbon disorder, is a distinct, therefore significant variation in mass and charge transport efficiency, and activities.

Identifying the impact of the covalent-bonded carbon matrix to FeN4 sites for acidic oxygen reduction

The atomic configurations of FeN_x moieties are the key to affect the activity of oxygen rection reaction (ORR). However, the traditional synthesis relying on high-temperature pyrolysis towards combining sources of Fe, N, and C often results in the plurality of local environments for the FeN_x sites. Unveiling the effect of carbon matrix adjacent to FeN_x sites towards ORR activity is important but still is a great challenge due to inevitable connection of diverse N as well as random defects. Here, we report a proof-of-concept study on the evaluation of covalent-bonded carbon environment connected to FeN₄ sites on their catalytic activity via pyrolysis-free approach. Basing on the closed π conjugated phthalocyanine-based intrinsic covalent organic polymers (COPs) with well-designed structures, we directly synthesized a series of atomically dispersed Fe-N-C catalysts with various pure carbon environments connected to the same FeN₄ sites. Experiments combined with density functional theory demonstrates that the catalytic activities of these COPs materials appear a volcano plot with the increasement of delocalized π electrons in their carbon matrix. The delocalized π electrons changed anti-bonding d-state energy level of the single FeN₄ moieties, hence tailored the adsorption between active centers and oxygen intermediates and altered the rate-determining step.

Unveiling the effect of carbon matrix adjacent to Fe-N towards oxygen reduction reaction is important yet challenging. Here the authors investigate the carbon environment covalent-connected to FeN₄ sites on their catalytic activity using models prepared by pyrolysis-free approach.

Multifunctional Electrocatalysis on Single-Site Metal Catalysts: A Computational Perspective

Multifunctional electrocatalysts are vastly sought for their applications in water splitting electrolyzers, metal-air batteries, and regenerative fuel cells because of their ability to catalyze multiple reactions such as hydrogen evolution, oxygen evolution, and oxygen reduction reactions. More specifically, the application of single-atom electrocatalyst in multifunctional catalysis is a promising approach to ensure good atomic efficiency, tunability and additionally benefits simple theoretical treatment. In this review, we provide insights into the variety of single-site metal catalysts and their identification. We also summarize the recent advancements in computational modeling of multifunctional electrocatalysis on single-site catalysts. Furthermore, we explain each modeling step with open-source-based working examples of a standard computational approach.

Spectroscopic discernibility of dopants and axial ligands in pyridinic FeN4 environments relevant to single-atom catalysts

Single-atom catalysts (SACs) activate small molecules, e.g. the oxygen reduction reaction is catalysed by FeNC materials. Because the nature of active site(s) in this type of SAC is unclear, spectroscopic and computational insights are needed to clarify the atomistic composition and electronic structure. Using quantum chemistry, we show that key features of [Fe{phen2A2}L]n+ complexes (A = CH, N with n = 0, A = O with n = 0, 2; L = OH-, Cl-) can be differentiated spectroscopically.

Porphyrin-based frameworks for oxygen electrocatalysis and catalytic reduction of carbon dioxide

The recent progress made on porphyrin-based frameworks and their applications in energy-related conversion technologies (e.g., ORR, OER and CO₂RR) and storage technologies (e.g., Zn–air batteries).

Circumventing Scaling Relations in Oxygen Electrochemistry Using Metal-Organic Frameworks

It has been well-established that unfavorable scaling relationships between *OOH, *OH, and *O are responsible for the high overpotentials associated with oxygen electrochemistry. A number of strategies have been proposed for breaking these linear constraints for traditional electrocatalysts (e.g., metals, alloys, metal-doped carbons); such approaches have not yet been validated experimentally for heterogeneous catalysts. Development of a new class of catalysts capable of circumventing such scaling relations remains an ongoing challenge in the field. In this work, we use density functional theory (DFT) calculations to demonstrate that bimetallic porphyrin-based MOFs (PMOFs) are an ideal materials platform for rationally designing the 3-D active site environments for oxygen reduction reaction (ORR). Specifically, we show that the *OOH binding energy and the theoretical limiting potential can be optimized by appropriately tuning the transition metal active site, the oxophilic spectator, and the MOF topology. Our calculations predict theoretical limiting potentials as high as 1.07 V for Fe/Cr-PMOF-Al, which exceeds the Pt/C benchmark for 4e ORR. More broadly, by highlighting their unique characteristics, this work aims to establish bimetallic porphyrin-based MOFs as a viable materials platform for future experimental and theoretical ORR studies.

X-ray Absorption Spectroscopy Investigation of Platinum–Gadolinium Thin Films with Different Stoichiometry for the Oxygen Reduction Reaction

Alloys of platinum and lanthanides present a remarkable activity for the oxygen reduction reaction—both in the form of extended surfaces and nanoparticulate catalysts. Co-sputter-deposited thin film catalysts based on platinum and gadolinium show great oxygen reduction activity improvement over pure Pt. The sputter-deposition technique represents a viable and versatile approach for investigating model catalyst systems with different compositions. In this work, co-sputtered Pt5Gd and Pt7.5Gd thin films were investigated using X-ray absorption spectroscopy as well as standardized electrochemical techniques. These investigations revealed the importance of forming alloys with specific stoichiometry, supporting the need of forming compressively strained Pt overlayers in order to achieve optimum catalytic performances.