Practical Applications of Grand-canonical Electronic Structure Calculations in Electrochemical Simulation

Modeling electrified interfaces has long been a great challenge in electrochemistry. In recent years, the grand-canonical treatment for electrons has gradually been developed, and its combination with density functional theory has been widely used to simulate electrochemical processes on an atomistic scale. In this Perspective, we aim to discuss several practical applications of this powerful technique after a short review of necessary fundamentals. We will begin with capacitor-based parametrization method of grand-canonical calculated results. If considering the electrodes under different applied potentials as different materials, the parametrization can be viewed as a kind of "quadratic scaling relation", which might reduce the overall computational costs by data postanalysis rather than algorithm development. Following an example of the abnormal potential-independent energetic curve within the bandgap area, we turn the topic to the semiconducting electrodes. Meanwhile, the specific behaviors of the bandgap also indicate that besides the reaction thermodynamics and kinetics, the detailed electronic structure of the system can also be well described by the grand-canonical treatment on electrons. Several possibilities for further applications are proposed correspondingly and summarized at the end of paper. We believe that the grand-canonical treatment for electronic structure calculations can greatly enrich our understanding of the fundamental mechanisms under electrochemical environments.

Mesoporous Single Atom-Cluster Fe-N/C Oxygen Evolution Electrocatalysts Synthesized with Bottlebrush Block Copolymer-Templated Rapid Thermal Annealing

Current electrocatalysts for oxygen evolution reaction (OER) are either expensive (such as IrO2, RuO2) or/and exhibit high overpotential as well as sluggish kinetics. This article reports mesoporous earth-abundant iron (Fe)-nitrogen (N) doped carbon electrocatalysts with iron clusters and closely surrounding Fe-N4 active sites. Unique to this work is that the mechanically stable mesoporous carbon-matrix structure (79 nm in pore size) with well-dispersed nitrogen-coordinated Fe single atom-cluster is synthesized via rapid thermal annealing (RTA) within only minutes using a self-assembled bottlebrush block copolymer (BBCP) melamine-formaldehyde resin composite template. The resulting porous structure and domain size can be tuned with the degree of polymerization of the BBCP backbone, which increases the electrochemically active surface area and improves electron transfer and mass transport for an effective OER process. The optimized electrocatalyst shows a required potential of 1.48 V (versus RHE) to obtain the current density of 10 mA/cm2 in 1 M KOH aqueous electrolyte and a small Tafel slope of 55 mV/decade at a given overpotential of 250 mV, which is significantly lower than recently reported earth-abundant electrocatalysts. Importantly, the Fe single-atom nitrogen coordination environment facilitates the surface reconstruction into a highly active oxyhydroxide under OER conditions, as revealed by X-ray photoelectron spectroscopy and in situ Raman spectroscopy, while the atomic clusters boost the single atoms reactive sites to prevent demetalation during the OER process. Density functional theory (DFT) calculations support that the iron nitrogen environment and reconstructed oxyhydroxides are electrocatalytically active sites as the kinetics barrier is largely reduced. This work has opened a new avenue for simple, rapid synthesis of inexpensive, earth-abundant, tailorable, mechanically stable, mesoporous carbon-coordinated single-atom electrocatalysts that can be used for renewable energy production.

Copper Single Atoms Chelated on Ligand-Modified Carbon for Ullmann-type C-O Coupling

Cross-coupling reactions are of great importance in chemistry due to their ability to facilitate the construction of complex organic molecules. Among these reactions, the Ullmann-type C-O coupling between phenols and aryl halides is particularly noteworthy and useful for preparing diarylethers. However, this reaction typically relies on homogeneous catalysts that rapidly deactivate under harsh reaction conditions. In this study, we introduce a novel heterogeneous catalyst for the Ullmann-type C-O coupling reaction, comprised of isolated Cu atoms chelated to a tetraethylenepentamine-pyrrole ligand that is immobilized on graphite nanoplatelets. The catalytic study reveals the recyclability of the material, and demonstrates the crucial role of the pyrrole linker in stabilizing the Cu sites. The work expands the potential of single-atom catalyst nanoarchitectures and underscores the significance of ligands in stabilizing metals in cationic forms, providing a novel, tailored catalyst for cross-coupling chemistries.

Modeling Single-Atom Catalysis

Electronic structure calculations represent an essential complement of experiments to characterize single-atom catalysts (SACs), consisting of isolated metal atoms stabilized on a support, but also to predict new catalysts. However, simulating SACs with quantum chemistry approaches is not as simple as often assumed. In this work, the essential factors that characterize a reliable simulation of SACs activity are examined. The Perspective focuses on the importance of precise atomistic characterization of the active site, since even small changes in the metal atom's surroundings can result in large changes in reactivity. The dynamical behavior and stability of SACs under working conditions, as well as the importance of adopting appropriate methods to solve the Schrödinger equation for a quantitative evaluation of reaction energies are addressed. The Perspective also focuses on the relevance of the model adopted. For electrocatalysis this must include the effects of the solvent, the presence of electrolytes, the pH, and the external potential. Finally, it is discussed how the similarities between SACs and coordination compounds may result in reaction intermediates that usually are not observed on metal electrodes. When these aspects are not adequately considered, the predictive power of electronic structure calculations is quite limited.

Elucidation of the Active Site for the Oxygen Evolution Reaction on a Single Pt Atom Supported on Indium Tin Oxide

Single-atom catalysts (SACs) have attracted attention for their high catalytic activity and selectivity, but the nature of their active sites under realistic reaction conditions, involving various ligands, is not well-understood. In this study, we use density functional theory calculations and grand canonical basin hopping to theoretically investigate the active site for the oxygen evolution reaction (OER) on a single Pt atom supported on indium tin oxide, including the influence of the electrochemical potential. We show that the ligands on the Pt atom change from Pt-OH in the absence of electrochemical potential to PtO(OH)₄ in electrochemical conditions. This change of the chemical state of Pt is associated with a decrease of 0.3 V for the OER overpotential. This highlights the importance of accurately identifying the nature of the active site under reaction conditions and the impact of adsorbates on the electrocatalytic activity. This theoretical investigation enhances our understanding of SACs for the OER.

Origin of the Rarely Reported High Performance of Mn-doped Carbon-based Oxygen Reduction Catalysts

Recent efforts to develop durable high-performance platinum-group metal (PGM)-free oxygen reduction reaction (ORR) electrocatalysts have focused on Fe- and Co-based molecular and pyrolyzed catalysts. While Mn-based catalysts have advantages of lower toxicity and higher durability, their activity has been generally poor. Nevertheless, several examples of high-performance Mn-based catalysts have been reported. Thus, it is necessary to understand why Mn-based materials much more rarely show high catalytic ORR performance and to determine the factors that can lead to the achievement of such high performance in these rare cases. We have studied the effects of the changes in the macrocycle structure, axial ligand, distance between the active sites, interactions with the dopant N atoms and the presence of an extended carbon network on the ORR catalysis of various Mn-, Fe-, and Co-based systems through the comparison of the adsorption energies of the ORR intermediates. We find that the sensitivity to the local environment changes is the largest for Mn and is the smallest for Co, with Fe between Mn and Co. Our results showed that the strong binding of OH by Mn and the strong sensitivity of the Mn to the modification of its environment necessitate a precise combination of local environment changes to achieve a high onset potential (V_onset ) in Mn-based catalysts. By contrast, the weaker binding of OH by Fe and Co and their weaker sensitivity to local environment changes lead to a wide variety of local environments with favorable catalytic activity (V_onset >0.7 V) for Co- and Fe-based systems. This explains the scarcity of reported Mn-based pyrolyzed catalysts and suggests that precise material synthesis and engineering of the active site can achieve high-performance Mn-based ORR electrocatalysts with high activity and durability.

Modeling Electrochemical Processes with Grand Canonical Treatment of Many-Body Perturbation Theory

Electrocatalysis plays a key role in sustainable energy conversion and storage. It is critical to model the grand canonical treatment of electrons, which accounts for the electrochemical potential explicitly, at the atomic scale and understand these reactions at electrified interfaces. However, such a grand canonical treatment for electrocatalytic modeling is in practice restricted to a treatment of electronic structure with density functional theory, and more accurate methods are in many cases desirable. Here, we develop an original workflow combining the grand canonical treatment of electrons with many-body perturbation theory in its random phase approximation (RPA). Using the potential dependent adsorption of carbon monoxide on the copper (100) facet, we show that the grand canonical RPA energetics provide the correct on-top Cu geometry for CO at reducing potential. The match with experimental results is significantly improved compared to the functionals at the generalized gradient approximation level, which is the most commonly used approximation for electrochemical applications. We expect this development to pave the way to further electrochemical applications using RPA.

Computational screening of transition-metal doped boron nanotubes as efficient electrocatalysts for water splitting

The search for efficient and low-cost electrocatalysts for the oxygen evolution reaction (OER) and the hydrogen evolution reaction (HER) is of utmost importance for the production of hydrogen and oxygen via water splitting. In this work, the catalytic performance of the OER and HER on transition metal doped boron nanotubes (BNTs) was investigated using density functional theory. It was found that the doped transition metal atoms determine the catalytic activity of the BNTs. Rhodium-doped BNTs exhibited excellent OER activity, while cobalt-doped BNTs displayed great catalytic activity toward the HER. Volcano relationships were found between the catalytic activity and the adsorption strength of reaction intermediates. Rhodium- and cobalt-doped BNTs exhibited great OER and HER catalytic activity due to the favorable adsorption strength of reaction intermediates. This work sheds light on the design of novel electrocatalysts for water splitting and provides helpful guidelines for the future development of advanced electrocatalysts.

Rhodium-doped BNTs demonstrated excellent OER activity, while cobalt-doped BNTs exhibited the best catalytic activity toward the HER among 12 different transition metal-doped BNTs.

Single-Atom (Iron-Based) Catalysts: Synthesis and Applications

Supported single-metal atom catalysts (SACs) are constituted of isolated active metal centers, which are heterogenized on inert supports such as graphene, porous carbon, and metal oxides. Their thermal stability, electronic properties, and catalytic activities can be controlled via interactions between the single-metal atom center and neighboring heteroatoms such as nitrogen, oxygen, and sulfur. Due to the atomic dispersion of the active catalytic centers, the amount of metal required for catalysis can be decreased, thus offering new possibilities to control the selectivity of a given transformation as well as to improve catalyst turnover frequencies and turnover numbers. This review aims to comprehensively summarize the synthesis of Fe-SACs with a focus on anchoring single atoms (SA) on carbon/graphene supports. The characterization of these advanced materials using various spectroscopic techniques and their applications in diverse research areas are described. When applicable, mechanistic investigations conducted to understand the specific behavior of Fe-SACs-based catalysts are highlighted, including the use of theoretical models.

Metal-Nitrogen-Carbon Catalysts of Specifically Coordinated Configurations toward Typical Electrochemical Redox Reactions

Metal-nitrogen-carbon (M-N-C) material with specifically coordinated configurations is a promising alternative to costly Pt-based catalysts. In the past few years, great progress is made in the studies of M-N-C materials, including the structure modulation and local coordination environment identification via advanced synthetic strategies and characterization techniques, which boost the electrocatalytic performances and deepen the understanding of the underlying fundamentals. In this review, the most recent advances of M-N-C catalysts with specifically coordinated configurations of M-N_x (x = 1-6) are summarized as comprehensively as possible, with an emphasis on the synthetic strategy, characterization techniques, and applications in typical electrocatalytic reactions of the oxygen reduction reaction, oxygen evolution reaction, hydrogen evolution reaction, CO₂ reduction reaction, etc., along with mechanistic exploration by experiments and theoretical calculations. Furthermore, the challenges and potential perspectives for the future development of M-N-C catalysts are discussed.

Carbon-Supported Single Metal Site Catalysts for Electrochemical CO2 Reduction to CO and Beyond

The electrochemical CO₂ reduction reaction (CO₂ RR) is a promising strategy to achieve electrical-to-chemical energy storage while closing the global carbon cycle. The carbon-supported single-atom catalysts (SACs) have great potential for electrochemical CO₂ RR due to their high efficiency and low cost. The metal centers' performance is related to the local coordination environment and the long-range electronic intercalation from the carbon substrates. This review summarizes the recent progress on the synthesis of carbon-supported SACs and their application toward electrocatalytic CO₂ reduction to CO and other C₁ and C₂ products. Several SACs are involved, including MN_x catalysts, heterogeneous molecular catalysts, and the covalent organic framework (COF) based SACs. The controllable synthesis methods for anchoring single-atom sites on different carbon supports are introduced, focusing on the influence that precursors and synthetic conditions have on the final structure of SACs. For the CO₂ RR performance, the intrinsic activity difference of various metal centers and the corresponding activity enhancement strategies via the modulation of the metal centers' electronic structure are systematically summarized, which may help promote the rational design of active and selective SACs for CO₂ reduction to CO and beyond.

Fe/N co-doped mesoporous carbon derived from cellulose-based ionic liquid as an efficient heterogeneous catalyst toward nitro aromatic compound reduction reaction

Iron and nitrogen-doped carbon substances with abundant active sites that related to dispersion of heteroatom species (Fe and N) on the surface of carbonous structure, are promising choice for eco-friendly catalytic reactions. Herein, cellulose-based ionic liquid (IL) derivative not only employed as the both nitrogen and iron heteroatom precursors, but also has been used as the green and biodegradable substrate. The non-noble Fe-NC@550, was successfully fabricated by convenient carbonization of cellulose-based IL. Further, the FeCl₄^- anion was used as the iron precursor and also it has been applied to elevate the SSA (specific surface area) of catalyst (from 40.96 to 160.42 m²/g) due to the presence of chlorine. On the basis of several pertinent physicochemical and experimental outcomes, the structure of the catalyst was successfully proved in different synthetic steps. As expected, the Fe-NC@550 exhibited the substantial efficiency toward hydrogenation of nitroarenes with high TOF value and also remarkable reusability.

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.

Transition-metal single atoms embedded into defective BC3 as efficient electrocatalysts for oxygen evolution and reduction reactions

Searching for high-activity, stable and low-cost catalysts toward oxygen evolution reaction (OER) and oxygen reduction reaction (ORR) are of significant importance to the development of renewable energy technologies. By using the computational screening method based on the density functional theory (DFT), we have systematically studied a wide range of transition metal (TM) atoms doped a defective BC3 monolayer (B atom vacancy VB and C atom vacancy VC), denoted as TM@VB and TM@VC (TM = Mn, Fe, Co, Ni, Cu, Ru, Rh, Pd, Ir and Pt), as efficient single atom catalysts for OER and ORR. The calculated results show that all the considered TM atoms can tightly bind with the defective BC3 monolayers to prevent the atomically dispersed atoms from clustering. The interaction strength between intermediates (HO*, O* and HOO*) and catalyst govern the catalytic activities of OER and ORR, which has a direct correlation with the d-band center (εd) of the TM active site that can be tuned by adjusting TM atoms with various d electron numbers. For TM@VB catalysts, it was found that the best catalyst for OER is Co@VB with an overpotential ηOER of 0.43 V, followed by Rh@VB (ηOER = 0.49 V), while for ORR, Rh@VB exhibits the lowest overpotential ηORR of 0.40 V, followed by Pd@VB (ηORR = 0.45 V). For TM@VC catalysts, the best catalyst for OER is Ni@VC (ηOER = 0.47 V), followed by Pt@VC (ηOER = 0.53 V), and for ORR, Pd@VC exhibits the highest activity with ηORR of 0.45 V. The results suggest that the high activity of the newly predicted well dispersed Rh@VB SAC is comparable to that of noble metal oxide benchmark catalysts for both OER and ORR. Importantly, Rh@VB may remain stable against dissolution at pH = 0 condition. The high energy barrier prevents the isolated Rh atom from clustering and ab initio molecule dynamic simulation (AIMD) result suggests that Rh@VB can remain stable under 300 K, indicating its kinetic stability. Our findings highlight a novel family of efficient and stable SAC based on carbon material, which offer a useful guideline to screen the metal active site for catalyst designation.

Demystifying the Atomistic Origin of the Electric Field Effect on Methane Oxidation

Understanding the role of an electric field on the surface of a catalyst is crucial in tuning and promoting the catalytic activity of metals. Herein, we evaluate the oxidation of methane over a Pt surface with varying oxygen coverage using density functional theory. The latter is controlled by the electrode polarization, giving rise to the non-Faradaic modification of catalytic activity phenomenon. At -1 V, the Pt(111) surface is present, while at 1 V, α-PtO₂ on Pt(111) takes over. Our results demonstrate that the alteration of the platinum oxide surface under the influence of an electrochemical potential promotes the catalytic activity of the metal oxides by lowering the activation energy barrier of the reaction.

Revisiting the Active Sites at the MoS2/H2O Interface via Grand-Canonical DFT: The Role of Water Dissociation

MoS₂ is a promising low-cost catalyst for the hydrogen evolution reaction (HER). However, the nature of the active sites remains a subject of debate. By taking the electrochemcal potential explicitly into account using grand-canonical density functional theory (DFT) in combination with the linearized Poisson-Boltzmann equation, we herein revisit the active sites of 2H-MoS₂. In addition to the well-known catalytically active edge sites, also specific point defects on the otherwise inert basal plane provide highly active sites for HER. Given that HER takes place in water, we also assess the reactivity of these active sites with respect to H₂O. The thermodynamics of proton reduction as a function of the electrochemical potential reveals that four edge sites and three basal plane defects feature thermodynamic overpotentials below 0.2 V. In contrast to current proposals, many of these active sites involve adsorbed OH. The results demonstrate that even though H₂O and OH block "active" sites, HER can also occur on these "blocked" sites, reducing protons on surface OH/H₂O entities. As a consequence, our results revise the active sites, highlighting the so far overlooked need to take the liquid component (H₂O) of the functional interface into account when considering the stability and activity of the various active sites.