Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe-N-C Materials

Unveiling the Composition-Dependent Catalytic Mechanism of Pt-Ni Alloys for Oxygen Reduction: A First-Principles Study

Unveiling the composition-dependent catalytic mechanism of Pt-based alloy cathodes for the oxygen reduction reaction (ORR) helps improve the proton exchange membrane fuel cells. Using density functional theory calculations, this study investigates the ORR catalytic performance of the Pt-Ni system with various compositions (1.00, ∼0.99, 0.75, 0.50, 0.25, ∼0.01, and 0.00). The ordered solid solution PtNi3(111) system shows activity comparable to Pt(111) and is cost-effective. The Ni1/Pt(111) system, featuring a single Ni atom on the Pt(111) surface as a surface single-atom alloy (SSAA), demonstrates the highest activity with an overpotential of only 0.28, which could be further reduced to 0.21 V by decreasing the surface Ni concentration to 1/16 monolayer coverage. The predicted high activity of Ni1/Pt(111) is confirmed when considering factors such as the implicit solution environment, constant potential conditions, and protonation capability. Moreover, surface-adsorbed oxygen species driven by reaction conditions stabilize these single Ni atoms of Ni1/Pt(111) by preventing segregation and dissolution processes, thereby exhibiting a dual functionality. This study reveals the composition dependence of Pt-based alloys and highlights the stability mechanisms of SSAA catalysts during the ORR.

Leveraging Curvature on N-Doped Carbon Materials for Hydrogen Storage

Carbon sorbent materials have shown great promise for solid-state hydrogen (H2) storage. Modification of these materials with nitrogen (N) dopants has been undertaken to develop materials that can store H2 at ambient temperatures. In this work density functional theory (DFT) calculations are used to systematically probe the influence of curvature on the stability and activity of undoped and N-doped carbon materials toward H binding. Specifically, four models of carbon materials are used: graphene, [5,5] carbon nanotube, [5,5] D5d-C120, and C60, to extract and correlate the thermodynamic properties of active sites with varying degrees of sp2 hybridization (curvature). From the calculations and analysis, it is found that graphitic N-doping is thermodynamically favored on more pyramidal sites with increased curvature. In contrast, it is found that the hydrogen binding energy is weakly affected by curvature and is dominated by electronic effects induced by N-doping. These findings highlight the importance of modulating the heteroatom doping configuration and the lattice topology when developing materials for H2 storage.

A General Descriptor for Single-Atom Catalysts with Axial Ligands

Decoration of an axial coordination ligand (ACL) on the active metal site is a highly effective and versatile strategy to tune activity of single-atom catalysts (SACs). However, the regulation mechanism of ACLs on SACs is still incompletely known. Herein, we investigate diversified combinations of ACL-SACs, including all 3d-5d transition metals and ten prototype ACLs. We identify that ACLs can weaken the adsorption capability of the metal atom (M) by raising the bonding energy levels of the M-O bond while enhancing dispersity of the d orbital of M. Through examination of various local configurations and intrinsic parameters of ACL-SACs, a general structure descriptor σ is constructed to quantify the structure-activity relationship of ACL-SACs which solely based on a few key intrinsic features. Importantly, we also identified the axial ligand descriptor σACL, as a part of σ, which can serve as a potential descriptor to determine the rate-limiting steps (RLS) of ACL-SACs in experiment. And we predicted several ACL-SACs, namely, CrN4-, FeN4-, CoN4-, RuN4-, RhN4-, OsN4-, IrN4- and PtN4-ACLs, that entail markedly higher activities than the benchmark catalysts of Pt and IrO2 for oxygen reduction reaction and oxygen evolution reaction, respectively, thereby supporting that the general descriptor σ can provide a simple and cost-effective method to assess efficient electrocatalysts.

Metal-Free Activation of Molecular Oxygen by Quaternary Ammonium-Based Ionic Liquid: A Detail Mechanistic Study

Most oxidation processes in common organic synthesis and chemical biology require transition metal catalysts or metalloenzymes. Herein, we report a detailed mechanistic study of a metal-free oxygen (O2) activation protocol on benzylamine/alcohols using simple quaternary alkylammonium-based ionic liquids to produce products such as amide, aldehyde, imine, and in some cases, even aromatized products. NMR and various control experiments established the product formation and reaction mechanism, which involved the conversion of molecular oxygen into a hydroperoxyl radical via a proton-coupled electron transfer process. Detection of hydrogen peroxide in the reaction medium using colorimetric analysis supported the proposed mechanism of oxygen activation. Furthermore, first-principles calculations using density functional theory (DFT) revealed that reaction coordinates and transition state spin densities have a unique spin conversion of triplet oxygen leading to formation of singlet products via a minimum energy crossing point. In addition to DFT, domain-based local pair natural orbital coupled cluster, (DLPNO-CCSD(T)), and complete active space self-consistent field, CASSCF(20,14) methods complemented the above findings. Partial density of states analysis showed stabilization of π* orbital of oxygen in the presence of ionic liquid, making it susceptible to hydrogen abstraction in a mild, metal-free condition. Inductively coupled plasma atomic emission spectroscopic (ICP-AES) analysis of reactant and ionic liquids clearly showed the absence of any significant transition metal contamination. The current results described the origin of O2 activation within the context of molecular orbital (MO) theory and opened up a new avenue for the use of ionic liquids as inexpensive, multifunctional and high-performance alternative to metal-based catalysts for O2 activation.

Inter-site structural heterogeneity induction of single atom Fe catalysts for robust oxygen reduction

Metal-nitrogen-carbon catalysts with hierarchically dispersed porosity are deemed as efficient geometry for oxygen reduction reaction (ORR). However, catalytic performance determined by individual and interacting sites originating from structural heterogeneity is particularly elusive and yet remains to be understood. Here, an efficient hierarchically porous Fe single atom catalyst (Fe SAs-HP) is prepared with Fe atoms densely resided at micropores and mesopores. Fe SAs-HP exhibits robust ORR performance with half-wave potential of 0.94 V and turnover frequency of 5.99 e-1s-1site-1 at 0.80 V. Theoretical simulations unravel a structural heterogeneity induced optimization, where mesoporous Fe-N4 acts as real active centers as a result of long-range electron regulation by adjacent microporous sites, facilitating O2 activation and desorption of key intermediate *OH. Multilevel operando characterization results identify active Fe sites undergo a dynamic evolution from basic Fe-N4 to active Fe-N3 under working conditions. Our findings reveal the structural origin of enhanced intrinsic activity for hierarchically porous Fe-N4 sites.

Modeling Interfacial Dynamics on Single Atom Electrocatalysts: Explicit Solvation and Potential Dependence

ConspectusSingle atom electrocatalysts, with noble metal-free composition, maximal atom efficiency, and exceptional reactivity toward various energy and environmental applications, have become a research hot spot in the recent decade. Their simplicity and the isolated nature of the atomic structure of their active site have also made them an ideal model catalyst system for studying reaction mechanisms and activity trends. However, the state of the single atom active sites during electrochemical reactions may not be as simple as is usually assumed. To the contrary, the single atom electrocatalysts have been reported to be under greater influence from interfacial dynamics, with solvent and electrolyte ions perpetually interacting with the electrified active center under an applied electrode potential. These complexities render the activity trends and reaction mechanisms derived from simplistic models dubious.In this Account, with a few popular single atom electrocatalysis systems, we show how the change in electrochemical potential induces nontrivial variation in the free energy profile of elemental electrochemical reaction steps, demonstrate how the active centers with different electronic structure features can induce different solvation structures at the interface even for the same reaction intermediate of the simplest electrochemical reaction, and discuss the implication of the complexities on the kinetics and thermodynamics of the reaction system to better address the activity and selectivity trends. We also venture into more intriguing interfacial phenomena, such as alternative reaction pathways and intermediates that are favored and stabilized by solvation and polarization effects, long-range interfacial dynamics across the region far beyond the contact layer, and the dynamic activation or deactivation of single atom sites under operation conditions. We show the necessity of including realistic aspects (explicit solvent, electrolyte, and electrode potential) into the model to correctly capture the physics and chemistry at the electrochemical interface and to understand the reaction mechanisms and reactivity trends. We also demonstrate how the popular simplistic design principles fail and how they can be revised by including the kinetics and interfacial factors in the model. All of these rich dynamics and chemistry would remain hidden or overlooked otherwise. We believe that the complexity at an electrochemical interface is not a curse but a blessing in that it enables deeper understanding and finer control of the potential-dependent free energy landscape of electrochemical reactions, which opens up new dimensions for further design and optimization of single atom electrocatalysts and beyond. Limitations of current methods and challenges faced by the theoretical and experimental communities are discussed, along with the possible solutions awaiting development in the future.

Gaining insight into the impact of electronic property and interface electrostatic field on ORR kinetics in alloy engineering via theoretical prognostication and experimental validation

Alloy engineering has been utilized as a potent strategy to modulate the oxygen reduction reaction （ORR） activity. However, the regulatory mechanism underpinning the ORR kinetics by means of alloy engineering is still shrouded in ambiguity. This work places emphasis on the kinetics of the ORR concerning Pt3M (M = Cr, Co, Cu, Pd, Sn, and Ir) catalysts, and integrates theoretical prognostication and experimental validation to illuminate the fundamental principles of alloy engineering. The ORR kinetic activity, as prognosticated by theory, shows significant agreement with experimental results, provided that the rate-determining step (RDS) accounts for a dominant role in the potential-independent kinetic mechanism. In essence, alloy engineering manipulates electronic properties through electron transfer to modulate intermediate adsorption and adjusts the interface electric field (Efield) to regulate hydrogen atom transport, ultimately influencing kinetics. The Efield holds greater significance in ORR kinetics compared to the intermediate adsorption (EadsO), the corresponding degrees of correlation with free energy barriers (Ea) of RDS are -0.89, and 0.75, respectively. This work highlights the nature of alloy engineering for ORR kinetics modulation and assists in the design of efficient catalysts.

Molecular Recognition Regulates Coordination Structure of Single-Atom Sites

Coordination engineering for single-atom sites has drawn increasing attention, yet its chemical synthesis remains a tough issue, especially for tailorable coordination structures. Herein, a molecular recognition strategy is proposed to fabricate single-atom sites with regulable local coordination structures. Specifically, a heteroatom-containing ligand serves as the guest molecule to induce coordination interaction with the metal-containing host, precisely settling the heteroatoms into the local structure of single-atom sites. As a proof of concept, thiophene is selected as the guest molecule, and sulfur atoms are successfully introduced into the local coordination structure of iron single-atom sites. Ultrahigh oxygen reduction electrocatalytic activity is achieved with a half-wave potential of 0.93 V versus reversible hydrogen electrode. Furthermore, the strategy possesses excellent universality towards diversified types of single-atom sites. This work makes breakthroughs in the fabrication of single-atom sites and affords new opportunities in structural regulation at the atomic level.

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.

Electrochemical Potential-Driven Shift of Frontier Orbitals in M-N-C Single-Atom Catalysts Leading to Inverted Adsorption Energies

Electronic structure is essential to understanding the catalytic mechanism of metal single-atom catalysts (SACs), especially under electrochemical conditions. This study delves into the nuanced modulation of "frontier orbitals" in SACs on nitrogen-doped graphene (N-C) substrates by electrochemical potentials. We observe shifts in Fermi level and changes of d-orbital occupation with alterations in electrochemical potentials, emphasizing a synergy between the discretized atomic orbitals of metals and the continuous bands of the N-C based environment. Using O2 and CO2 as model adsorbates, we highlight the direct consequences of these shifts on adsorption energies, unveiling an intriguing inversion of adsorption energies on Co/N-C SAC under negative electrochemical potentials. Such insights are attributed to the role of the dxz and dz2 orbitals, pivotal for stabilizing the π* orbitals of O2. Through this exploration, our work offers insights on the interplay between electronic structures and adsorption behaviors in SACs, paving the way for enhanced catalyst design strategies in electrochemical processes.

Leveraging Interlayer Interaction in M-N-C Catalysts for Enhanced Activity in Oxygen Reduction Reactions

Atomically dispersed metal-nitrogen-carbon (M-N-C) materials are deemed promising catalysts for the oxygen reduction reaction (ORR) in fuel cells. Yet the multilayer nature of M-N-C has been largely neglected in computational analysis. To bridge the gap, we conducted a first-principles investigation using bilayer M-N-C models (TMNx/G-TMNy/G, TM = Mn, Fe, Co, Ni, Cu, G = graphene, x, y = 3 or 4), where the TMs on the top serves as the active center. While in-plane TMN4 at the bottom has a minimal impact on the ORR, out-of-plane TMN3 substantially influences the adsorption free energy of OH through a strong interlayer bonding interaction. By leveraging interlayer interactions, we appreciably lowered the overpotential of selected TMN4 (TM = Co, Ni, Cu) and achieved a minimum of 0.40 V on CoN4/G-CuN3/G. Constant potential calculations revealed weak dependence of OH binding energy on external voltage and obtained results comparable to constant charge calculation. This study provided new physical insight into modulating naturally occurring multilayer M-N-C catalysts.

Targeted Spin-State Regulation to Boost Oxygen Reduction Reaction

Catalytic reactions are known to be significantly affected by spin states and their variations during reaction processes, yet the mechanisms behind them remain not fully understood, thus preventing the rational optimization of catalysis. Here, we explore the relationship between the spin states of active sites and their catalytic performance, taking the oxygen reduction reaction as an example. We demonstrate that the catalytic performance is spin-state-dependent and can be improved by adjusting spin states during the catalytic process. To this end, we further investigate the possibility of altering the spin states of transition metals through the application of external fields, such as adsorbed species. By studying the influence of the strength of adsorbed ligands on spin states and its impact on catalytic performance, our results show that optimal catalytic performance is achieved when the strength of the external field is neither too strong nor too weak, forming a volcano-like relationship between the catalytic performance and the external field strength. Our findings can have far-reaching implications for the rational design of high-performance catalysis.

Collective Behavior of Single-Atom Catalysts: A Synergistic Effect between Strain and Site Configuration

Fe single-atom catalysts on N-doped graphene (Fe-NC) exhibit good and variable catalytic activity linked to the active site density and configuration. Here, we comprehensively investigate the Fe-NC catalysts under various strained states and site densities to address the interplay between the active site density, local strain, site geometry, and oxygen evolution reaction (OER) activity. It is found that the active site density is closely associated with in-plane strain, which can be tuned by popping up Fe single atoms from the graphene film and, thereby, modulating the OER catalytic activity. Further analysis indicates that there exist three orientations of the FeN4 active site, each introducing specific anisotropic strain. As a result, the in-plane strain correlates with both the orientation and density of the active site, ultimately influencing catalytic activity. Our findings demonstrate the synergistic effects of multiple factors in single-atom catalysts, providing new insights into the rational design and fine tuning single-atom catalysts via collective interactions.

Bifunctional Single Atom Catalysts for Rechargeable Zinc-Air Batteries: From Dynamic Mechanism to Rational Design

Ever-growing demands for rechargeable zinc-air batteries (ZABs) call for efficient bifunctional electrocatalysts. Among various electrocatalysts, single atom catalysts (SACs) have received increasing attention due to the merits of high atom utilization, structural tunability, and remarkable activity. Rational design of bifunctional SACs relies heavily on an in-depth understanding of reaction mechanisms, especially dynamic evolution under electrochemical conditions. This requires a systematic study in dynamic mechanisms to replace current trial and error modes. Herein, fundamental understanding of dynamic oxygen reduction reaction and oxygen evolution reaction mechanisms for SACs is first presented combining in situ and/or operando characterizations and theoretical calculations. By highlighting structure-performance relationships, rational regulation strategies are particularly proposed to facilitate the design of efficient bifunctional SACs. Furthermore, future perspectives and challenges are discussed. This review provides a thorough understanding of dynamic mechanisms and regulation strategies for bifunctional SACs, which are expected to pave the avenue for exploring optimum single atom bifunctional oxygen catalysts and effective ZABs.

Electronic Structure Modulation of Fe-N4-C for Oxygen Evolution Reaction via Transition Metal Dopants and Axial Ligands

The popular single-atom catalyst (SAC) Fe-N4 is generally believed to be an excellent oxygen reduction reaction (ORR) electrocatalyst, which is less active in the oxygen evolution reaction (OER). Herein, FeM-N6 configuration catalysts (M = Fe, Co, Ni, Cu, Ag, and Au) were constructed for the oxygen evolution reaction by embedding M dopants on Fe-N4 systems based on the density functional theory. The electronic structure analysis reveals that the Fe-M metal interactions play dominant roles in regulating the d orbital distributions of Fe sites, which in turn alter the catalytic OER performance. Subsequent thermodynamic results indicate that the potential-determining step (PDS) for all catalysts is the formation of OOH*, which exhibits a tendency of decreased overpotentials with enhanced metal interactions. Apart from these, the effects of axial ligands on the OER activity of the catalysts in practical conditions were considered. Generally, most of the axial ligands are found to be thermodynamically favorable for the OER process. Interestingly, a competitive relationship of the electrons from the d orbital of Fe sites was found between the axial ligand and the adsorbed intermediate species during the reaction, which raises the energy barrier for OH* to O* conversion and can even alter the PDS in certain cases. The present work sheds new light on the design of future high-performance OER catalysts.

Iron clusters regulate local charge distribution in Fe-N4 sites to boost oxygen electroreduction

The atomically-dispersed and nitrogen-coordinated iron (FeNC) on a carbon catalyst is a potential non-noble metal catalyst that can replace precious metal electrocatalysts. However, its activity is often unsatisfactory owing to the symmetric charge distribution around the iron matrix. In this study, atomically- dispersed Fe-N₄ and Fe nanoclusters loaded with N-doped porous carbon (Fe_NCs/Fe_SAs-NC-Z8@34) were rationally fabricated by introducing homologous metal clusters and increasing the N content of the support. Fe_NCs/Fe_SAs-NC-Z8@34 exhibited a half-wave potential of 0.918 V, which exceeded that of the commercial benchmark Pt/C catalyst. Theoretical calculations verified that introducing Fe nanoclusters can break the symmetric electronic structure of Fe-N₄, thus inducing charge redistribution. Furthermore, it can optimize a part of Fe 3d occupancy orbitals and accelerate OO fracture in OOH* (rate-determining step), thus significantly improving oxygen reduction reaction activity. This work provides a reasonably advanced pathway to modulate the electronic structure of the single-atom center and optimize the catalytic activity of single-atom catalysts.

Exploring the catalytic activity of graphene-based TM-NxC4-x single atom catalysts for the oxygen reduction reaction via density functional theory calculation

Electrocatalysts for the oxygen reduction reaction (ORR) are extremely crucial for advanced energy conversion technologies, such as fuel cell batteries. A promising ORR catalyst usually should have low overpotentials, rich catalytic sites and low cost. In the past decade, single-atom catalyst (SAC) TM-N₄ (TM = Fe, Co, etc.) embedded graphene matrixes have been widely studied for their promising performance and low cost for ORR catalysis, but the effect of coordination on the ORR activity is not fully understood. In this work, we will employ density functional theory (DFT) calculations to systematically investigate the ORR activity of 40 different 3d transition metal single-atom catalysts (SACs) supported on nitrogen-doped graphene supports, ranging from vanadium to zinc. Five different nitrogen coordination configurations (TM-N_xC_4-x with x = 0, 1, 2, 3, and 4) were studied to reveal how C/N substitution affects the ORR activity. By looking at the stability, free energy diagram, overpotential, and scaling relationship, our calculation showed that partial C substitution can effectively improve the ORR performance of Mn, Co, Ni, and Zn-based SACs. The volcano plot obtained from the scaling relationship indicated that the substitution of N by C could distinctively affect the potential-limiting step in the ORR, which leads to the enhanced or weakened ORR performance. Density of states and d-band center analysis suggested that this coordination-tuned ORR activity can be explained by the shift of the d-band center due to the coordination effect. Finally, four candidates with optimal ORR activity and dynamic stability were proposed from the pool: NiC₄, CoNC₃, CrN₄, and ZnN₃C. Our work provides a feasible designing strategy to improve the ORR activity of graphene-based TM-N₄ SACs by tuning the coordination environment, which may have potential implication in the high-performance fuel cell development.

Strongly bound Wannier-Mott exciton in pristine (LaO)MnAs and origin of ferrimagnetism in F-doped (LaO)MnAs

We study the electronic, magnetic, and optical properties of (LaO1-xFx)MnAs (x = 0, 0.0625, 0.125, 0.25) systems, calculated using the generalized gradient approximation (GGA) corrected by Hubbard energy (U) = 1 eV. For x = 0, this system shows equal bandgap (Eg) values for spin-up and spin-down of 0.826 eV, with antiferromagnetic (AFM) properties and local magnetic moment in the Mn site of 3.86 μB per Mn. By doping F with x = 0.0625, the spin-up and spin-down Eg values decrease to 0.778 and 0.798 eV, respectively. This system, along with antiferromagnetic properties, also has a local magnetic moment in the Mn site of 3.83 μB per Mn. Increasing doping F to x = 0.125 induces increases of Eg to 0.827 and 0.839 eV for spin-up and spin-down. However, the AFM remains, where μMn slightly decreases to 3.81 μB per Mn. Furthermore, the excess electron from the F ion induces the Fermi level to move toward the conduction band and changes the bandgap type from indirect bandgap (Γ → M) to direct bandgap (Γ → Γ). Increasing x to 25% induces the decrease of spin-up and spin-down Eg to 0.488 and 0.465 eV, respectively. This system shows that the AFM changes to ferrimagnetism (FIM) for x = 25%, with a total magnetic moment of 0.78 μB per cell, which is mostly contributed by Mn 3d and As 4p local magnetic moments. The change from AFM to FIM behavior results from competition between superexchange AFM ordering and Stoner's exchange ferromagnetic ordering. Pristine (LaO)MnAs exhibits high excitonic binding energy (∼146.5 meV) due to a flat band structure. Our study shows that doping F in the (LaO)MnAs system significantly modifies the electronic, magnetic, and optical properties for novel advanced device applications.

Potential-Dependent Active Moiety of Fe-N-C Catalysts for the Oxygen Reduction Reaction

The real active moiety of Fe-N-C single-atom catalysts (SACs) during the oxygen reduction reaction (ORR) depends on the applied potential. Here, we examine the ORR activity of various SAC active moieties (Fe-N₄, Fe-(OH)N₄, Fe-(O₂)N₄, and Fe-(OH₂)N₄) over a wide potential window ranging from -0.8 to 1.0 V (vs. SHE) using constant potential density functional theory calculations. We show that the ORR activity of the Fe-N₄ moiety is hindered by the slow *OH protonation, while the Fe-(OH₂)N₄ (0.4 V ≤ U ≤ 1.0 V), *O₂-assisted Fe-N₄ (-0.6 V ≤ U ≤ 0.2 V), and Fe-(OH)N₄ (U = -0.8 V) moieties dominate the ORR activity of the Fe-N-C catalysts at different potential windows. These oxygenated species modified the single-atom Fe sites and can promote *OH protonation by regulating the electron occupancy of the Fe 3d_z² (spin-up) and Fe 3d_xz (spin-down) orbitals. Overall, our findings provide guidance for understanding the active moieties of SACs.

FeNC Oxygen Reduction Electrocatalyst with High Utilization Penta-Coordinated Sites

Atomic Fe in N-doped carbon (FeNC) electrocatalysts for oxygen (O₂ ) reduction at the cathode of proton exchange membrane fuel cells are the most promising alternative to platinum-group-metal catalysts. Despite recent progress on atomic FeNC O₂  reduction, their controlled synthesis and stability for practical applications remain challenging. A two-step synthesis approach has recently led to significant advances in terms of Fe-loading and mass activity; however, the Fe utilization remains low owing to the difficulty of building scaffolds with sufficient porosity that electrochemically exposes the active sites. Herein, this issue is addressed by coordinating Fe in a highly porous nitrogen-doped carbon support (≈3295 m²  g^-1 ), prepared by pyrolysis of inexpensive 2,4,6-triaminopyrimidine and a Mg²⁺ salt active site template and porogen. Upon Fe coordination, a high electrochemical active site density of 2.54 × 10¹⁹  sites g_FeNC ^-1  and a record 52% FeN_x electrochemical utilization based on in situ nitrite stripping are achieved. The Fe single atoms are characterized pre- and post-electrochemical accelerated stress testing by aberration-corrected high-angle annular dark field scanning transmission electron microscopy, showing no Fe clustering. Moreover, ex situ X-ray absorption spectroscopy and low-temperature Mössbauer spectroscopy suggest the presence of penta-coordinated Fe sites, which are further studied by density functional theory calculations.

Inductive Effect Alone Cannot Explain Lewis Adduct Formation and Dissociation at Electrode Interfaces

Understanding breaking and formation of Lewis bonds at an electrified interface is relevant to a large range of phenomena, including electrocatalysis and electroadsorption. The complexities of interfacial environments and associated reactions often impede a systematic understanding of this type of bond at interfaces. To address this challenge, we report the creation of a main group classic Lewis acid-base adduct on an electrode surface and its behavior under varying electrode potentials. The Lewis base is a self-assembled monolayer of mercaptopyridine and the Lewis acid is BF₃, forming a Lewis bond between nitrogen and boron. The bond is stable at positive potentials but cleaves at potentials more negative of approximately -0.3 V vs Ag/AgCl without an associated current. We also show that if the Lewis acid BF₃ is supplied from a reservoir of Li⁺BF₄^- electrolyte, the cleavage is completely reversible. We propose that the N-B Lewis bond is affected both by the field-induced intramolecular polarization (electroinduction) and by the ionic structures and ionic equilibria near the electrode. Our results indicate that the second effect is responsible for the Lewis bond cleavage at negative potentials. This work is relevant to understanding the fundamentals of electrocatalytic and electroadsorption processes.

Candied Haws-Like Fe-N-C Catalysts with Broadened Carbon Interlayer Spacing for Efficient Zinc-Air Battery

As efficient nonprecious metal catalysts for oxygen reduction reaction (ORR), Fe-N-C materials are one of the most promising alternatives to Pt-based catalysts for fuel cells and metal-air batteries. However, the intrinsically low density of key active sites like FeN₄ moieties hampers their commercial applications. Herein, we provide a smart strategy to construct a candied haws-like Fe-N-C catalyst (CH-FeNC) with broadened carbon interplanar spacing (>4 Å), starting with trehalose as a structure-built brick coupled with a zinc-zeolite imidazole framework (ZIF-8) and polyaniline (PANI) and then followed by copyrolysis carbonization of them. The obtained CH-FeNC exhibits half-wave potentials of 0.92 and 0.90 V (vs RHE) before and after 10,000 cycles in 0.1 M KOH, which are superior to the 0.90 and 0.85 V obtained by commercial Pt/C for ORR. The power density of a homemade zinc-air battery equipped with the catalyst is up to 131 mW cm^-2, greater than that of Pt/C (124 mW cm^-2). The extended X-ray absorption fine structure (EXAFS) results and density functional theory (DFT) theoretical calculations reveal that there exists enriched zigzag or armchair edge-hosted FeN₄ active sites, located at the abundant interface between carbon components in this composite. Furthermore, the unique broadened carbon interlayer spacing plays a key role in deciding the ORR rate in alkaline but not in acidic environments because there exists a fifth ligand of active Fe in the form of FeN₄ centers coupled with SO₄^2- and ClO₄^- from acids.

Exploring the Potential Energy Surface of Pt6 Sub-Nano Clusters Deposited over Graphene

Catalytic systems based on sub-nanoclusters deposited over different supports are promising for very relevant chemical transformations such as many electrocatalytic processes as the ORR. These systems have been demonstrated to be very fluxional, as they are able to change shape and interconvert between each other either alone or in the presence of adsorbates. In addition, an accurate representation of their catalytic activity requires the consideration of ensemble effects and not a single structure alone. In this sense, a reliable theoretical methodology should assure an accurate and extensive exploration of the potential energy surface to include all the relevant structures and with correct relative energies. In this context, we applied DFT in conjunction with global optimization techniques to obtain and analyze the characteristics of the many local minima of Pt₆ sub-nanoclusters over a carbon-based support (graphene)-a system with electrocatalytic relevance. We also analyzed the magnetism and the charge transfer between the clusters and the support and paid special attention to the dependence of dispersion effects on the ensemble characteristics. We found that the ensembles computed with and without dispersion corrections are qualitatively similar, especially for the lowest-in-energy clusters, which we attribute to a (mainly) covalent binding to the surface. However, there are some significant variations in the relative stability of some clusters, which would significantly affect their population in the ensemble composition.

Correlation between Electronic Descriptor and Proton-Coupled Electron Transfer Thermodynamics in Doped Graphite-Conjugated Catalysts

Graphite-conjugated catalysts (GCCs) provide a powerful framework for investigating correlations between electronic structure features and chemical reactivity of single-site heterogeneous catalysts. GCC-phenazine undergoes proton-coupled electron transfer (PCET) involving protonation of phenazine at its two nitrogen atoms with the addition of two electrons. Herein, this PCET reaction is investigated in the presence of defects, such as heteroatom dopants, in the graphitic surface. The proton-coupled redox potentials, E_PCET, are computed using a constant potential periodic density functional theory (DFT) strategy. The electronic states directly involved in PCET for GCC-phenazine exhibit the same nitrogen orbital character as those for molecular phenazine. The energy ε_LUS of this phenazine-related lowest unoccupied electronic state in GCC-phenazine is identified as a descriptor for changes in PCET thermodynamics. Importantly, ε_LUS is obtained from only a single DFT calculation but can predict E_PCET, which requires many such calculations. Similar electronic features may be useful descriptors for thermodynamic properties of other single-site catalysts.

Metal-Organic Framework-Derived Atomically Dispersed Co-N-C Electrocatalyst for Efficient Oxygen Reduction Reaction

In this work, an atomically dispersed cobalt-nitrogen-carbon (Co-N-C) catalyst is prepared for the oxygen reduction reaction (ORR) by using a metal-organic framework (MOF) as a self-sacrifice template under high-temperature pyrolysis. Spherical aberration-corrected electron microscopy is employed to confirm the atomic dispersion of high-density Co atoms on the nitrogen-doped carbon scaffold. The X-ray photoelectron spectroscopy results verify the existence of Co-N-C active sites and their content changes with the Co content. The electrochemical results show that the electrocatalytic activity shows a volcano-shaped relationship, which increases with the Co content from 0 to 0.99 wt.% and then decreases when the presence of Co nanoparticles at 1.61 wt.%. The atomically dispersed Co-N-C catalyst with Co content of 0.99 wt.% shows an onset potential of 0.96 V vs. reversible hydrogen electrode (RHE) and a half-wave potential of 0.89 V vs. RHE toward ORR. The excellent ORR activity is attributed to the high density of the Co-N-C sites with high intrinsic activity and high specific surface area to expose more active sites.

What is the Real Origin of the Activity of Fe-N-C Electrocatalysts in the O2 Reduction Reaction? Critical Roles of Coordinating Pyrrolic N and Axially Adsorbing Species

Fe-N-C electrocatalysts have emerged as promising substitutes for Pt-based catalysts for the oxygen reduction reaction (ORR). However, their real catalytic active site is still under debate. The underlying roles of different types of coordinating N including pyridinic and pyrrolic N in catalytic performance require thorough clarification. In addition, how to understand the pH-dependent activity of Fe-N-C catalysts is another urgent issue. Herein, we comprehensively studied 13 different N-coordinated FeN_xC configurations and their corresponding ORR activity through simulations which mimic the realistic electrocatalytic environment on the basis of constant-potential implicit solvent models. We demonstrate that coordinating pyrrolic N contributes to a higher activity than pyridinic N, and pyrrolic FeN₄C exhibits the highest activity in acidic media. Meanwhile, the in situ active site transformation to *O-FeN₄C and *OH-FeN₄C clarifies the origin of the higher activity of Fe-N-C in alkaline media. These findings can provide indispensable guidelines for rational design of better durable Fe-N-C catalysts.