Structural and Electronic Descriptors of Catalytic Activity of Graphene-Based Materials: First-Principles Theoretical Analysis

Ultrathin, large area β-Ni(OH)2 crystalline nanosheet as bifunctional electrode material for charge storage and oxygen evolution reaction

Bifunctional electrode materials are highly desirable for meeting increasing global energy demands and mitigating environmental impact. However, improving the atom-efficiency, scalability, and cost-effectiveness of storage systems, as well as optimizing conversion processes to enhance overall energy utilization and sustainability, remains a significant challenge for their application. Herein, we devised an optimized, facile, economic, and scalable synthesis of large area (cm2), ultrathin (∼2.9 ± 0.3 nm) electroactive nanosheet of β-Ni(OH)2, which acted as bifunctional electrode material for charge storage and oxygen evolution reaction (OER). The β-Ni(OH)2 nanosheet electrode shows the volumetric capacity of 2.82 Ah.cm-3(0.82 µAh.cm-2) at the current density of 0.2 mA.cm-2. The device shows a high capacity of 820 mAh.cm-3 with an ultrahigh volumetric energy density of 0.33 Wh.cm-3 at 275.86 W.cm-3 along with promising stability (30,000 cycles). Furthermore, the OER activity of ultrathin β-Ni(OH)2 exhibits an overpotential (η10) of 308 mV and a Tafel value of 42 mV dec-1 suggesting fast reaction kinetics. The mechanistic studies are enlightened through density functional theory (DFT), which reveals that additional electronic states near the Fermi level enhance activity for both capacitance and OER.

Computational chemistry for water-splitting electrocatalysis

Electrocatalytic water splitting driven by renewable electricity has attracted great interest in recent years for producing hydrogen with high-purity. However, the practical applications of this technology are limited by the development of electrocatalysts with high activity, low cost, and long durability. In the search for new electrocatalysts, computational chemistry has made outstanding contributions by providing fundamental laws that govern the electron behavior and enabling predictions of electrocatalyst performance. This review delves into theoretical studies on electrochemical water-splitting processes. Firstly, we introduce the fundamentals of electrochemical water electrolysis and subsequently discuss the current advancements in computational methods and models for electrocatalytic water splitting. Additionally, a comprehensive overview of benchmark descriptors is provided to aid in understanding intrinsic catalytic performance for water-splitting electrocatalysts. Finally, we critically evaluate the remaining challenges within this field.

Modulated Ultrathin NiCo-LDH Nanosheet-Decorated Zr3+-Rich Defective NH2-UiO-66 Nanostructure for Efficient Photocatalytic Hydrogen Evolution

Defect engineering through modification of their surface linkage is found to be an effective pathway to escalate the solar energy conversion efficiency of metal-organic frameworks (MOFs). Herein, defect engineering using controlled decarboxylation on the NH2-UiO-66 surface and integration of ultrathin NiCo-LDH nanosheets synergizes the hydrogen evolution reaction (HER) under a broad visible light regime. Diversified analytical methods including positron annihilation lifetime spectroscopy were employed to investigate the role of Zr3+-rich defects by analyzing the annihilation characteristics of positrons in NH2-UiO-66, which provides a deep insight into the effects of structural defects on the electronic properties. The progressively tuned photophysical properties of the NiCo-LDH@NH2-UiO-66-D-heterostructured nanocatalyst led to an impressive rate of HER (∼2458 μmol h-1 g-1), with an apparent quantum yield of ∼6.02%. The ultrathin NiCo-LDH nanosheet structure was found to be highly favored toward electrostatic self-assembly in the heterostructure for efficient charge separation. Coordination of Zr3+ on the surface of the NiCo-LDH nanosheet support through NH2-UiO-66 was confirmed by X-ray absorption spectroscopy and electron paramagnetic resonance spectroscopy techniques. Femtosecond transient absorption spectroscopy studies unveiled a photoexcited charge migration process from MOF to NiCo-LDH which favorably occurred on a picosecond time scale to boost the catalytic activity of the composite system. Furthermore, the experimental finding and HER activity are validated by density functional theory studies and evaluation of the free energy pathway which reveals the strong hydrogen binding over the surface and infers the anchoring effect of the ultrathin layered double hydroxide (LDH) in the vicinity of the Zr cluster with a strong host-guest interaction. This work provided a novel insight into efficient photocatalysis via defect engineering at the linker modulation in MOFs.

Theoretical verification on adsorptive removal of caffeine by carbon and nitrogen-based surfaces: Role of charge transfer, π electron occupancy, and temperature

Eliminating an emerging water pollutant, caffeine molecules, from an aqueous solution using carbon and nitrogen-based adsorbents is of significant interest to public health. These adsorbents have been shown to have decent adsorption capacity toward caffeine due to their surface functionality. Therefore, screening various carbon and nitrogen-based surfaces can be a better option for high-performance adsorbents to remove caffeine efficiently from wastewater. Herein, we present combined first principles and molecular dynamics quantification of the adsorption enthalpies of caffeine molecules on the possible active sites of carbon and nitrogen-based adsorbents (graphene, phagraphene, graphdiyne, single-wall carbon nanotube, fullerene, and graphitic carbon nitride) with the incorporation of Van der Waals interactions. From the DFT calculations, N-doped carbon surfaces show the highest adsorption energies of single and dimer CAF compared to pristine carbon-based adsorbents. A charge density difference and Bader charge analysis display that high charge transfer occurs between the caffeine's oxygen and the surface's nitrogen atoms. An abundance of π-electrons from the nitrogen atoms, composed of large electron clouds of aromatic rings on the graphitic carbon surface, tends to favor extensive π-π interactions with the caffeine molecule. The high value of pz electron occupancy (1.445) of N in the hexagonal ring of the graphitic surface transfers additional charge transfer, which leads to strong adsorption energy of CAF than pristine surfaces. Also, the g-C3N4 surface adsorbs the CAF molecule with higher adsorption than other N-doped carbon surfaces due to the high pz_eo (1.5448) of N atoms on the surface. At 310 K, the water molecules' kinetics aids the single and dimer caffeine molecules to adsorb with the highest adsorption energies on the active sites of g-C3N4 surfaces than graphene adsorbent.

Characteristics of the Frustrated Lewis Pairs (FLPs) on the Surface of Albite and the Corresponding Mechanism of H2 Activation

The characteristics of frustrated Lewis pairs (FLPs) on albite surfaces were analyzed with density functional theory, and the reaction mechanism for H2 activation by the FLPs was studied. The results show that albite is an ideal substrate material with FLPs, and its (001) and (010) surfaces have the typical characteristics of FLPs. In the case of H2 activation, the interaction between the HOMO of H2 and the SOMO of the Lewis base and the electron acceptance characteristics of the Lewis acid are the key factors. In fact, the activation energy of H2 is the required activation energy from the ground state to the excited state, and once the excited state is produced, the dissociative adsorption of H2 will occur directly. This study provides a new ideas and a reference for research on the construction of novel solid FLPs catalysts using ultramicro channel materials.

Computational Insight into TM-N x Embedded Graphene Bifunctional Electrocatalysts for Oxygen Evolution and Reduction Reactions

Due to the energy crisis, development of bifunctional electrocatalysts for both oxygen evolution and reduction reactions is highly demanding. In this study, we have systematically investigated the bifunctional activity of metal (Co/Rh/Ir) and N co-doped graphene systems with varying N-dopant concentrations (TM-N _x @G, x = 0, 2, 4) using first-principles calculations. Charge transfer from the metal sites to the adsorbed intermediates and the adsorption free energy of the intermediates play important roles to help understand the potential-determining step and overpotential values for oxygen evolution reaction (OER)/oxygen reduction reaction (ORR). A dual volcano plot for all the systems using a common descriptor ΔG _OH* has been constructed. We find that the systems having ΔG _OH* values in the range of 0.40-0.70 eV can act as bifunctional electrocatalysts. Our study not only highlights the importance of metal and non-metal co-doped graphene as bifunctional catalysts but also can serve as a promising strategy for the design of efficient OER/ORR electrocatalysts.

Designing 3d metal oxides: selecting optimal density functionals for strongly correlated materials

Transition metal oxides (TMOs) have remarkable physicochemical properties, are non-toxic, and have low cost and high annual production, thus they are commonly studied for various technological applications. Density functional theory (DFT) can help to optimize TMO materials by providing insights into their electronic, optical and thermodynamic properties, and hence into their structure-performance relationships, over a wide range of solid-state structures and compositions. However, this is underpinned by the choice of the exchange-correlation (XC) functional, which is critical to accurately describe the highly localized and correlated 3d-electrons of the transition metals in TMOs. This tutorial review presents a benchmark study of density functionals (DFs), ranging from generalized gradient approximation (GGA) to range-separated hybrids (RSH), with the all-electron def2-TZVP basis set, comparing magneto-electro-optical properties of 3d TMOs against experimental observations. The performance of the DFs is assessed by analyzing the band structure, density of states, magnetic moment, structural static and dynamic parameters, optical properties, spin contamination and computational cost. The results disclose the strengths and weaknesses of the XC functionals, in terms of accuracy, and computational efficiency, suggesting the unprecedented PBE0-1/5 as the best candidate. The findings of this work contribute to necessary developments of XC functionals for periodic systems, and materials science modelling studies, particularly informing how to select the optimal XC functional to obtain the most trustworthy description of the ground-state electron structure of 3d TMOs.

Strategic Modulation of Target-Specific Isolated Fe,Co Single-Atom Active Sites for Oxygen Electrocatalysis Impacting High Power Zn-Air Battery

An effective modulation of the active sites in a bifunctional electrocatalyst is essentially desired, and it is a challenge to outperform the state-of-the-art catalysts toward oxygen electrocatalysis. Herein, we report the development of a bifunctional electrocatalyst having target-specific Fe-N₄/C and Co-N₄/C isolated active sites, exhibiting a symbiotic effect on overall oxygen electrocatalysis performances. The dualism of N-dopants and binary metals lower the d-band centers of both Fe and Co in the Fe,Co,N-C catalyst, improving the overpotential of the overall electrocatalytic processes (ΔE_ORR-OER = 0.74 ± 0.02 V vs RHE). Finally, the Fe,Co,N-C showed a high areal power density of 198.4 mW cm^-2 and 158 mW cm^-2 in the respective liquid and solid-state Zn-air batteries (ZABs), demonstrating suitable candidature of the active material as air cathode material in ZABs.

N-Doped Carbon Electrocatalyst: Marked ORR Activity in Acidic Media without the Contribution from Metal Sites?

Fe-N-C electrocatalysts have been demonstrated to be the most promising substitutes for benchmark Pt/C catalysts for the oxygen reduction reaction (ORR). Herein, we report that N-doped carbon materials with trace amounts of iron (0-0.08 wt. %) show excellent ORR activity and durability comparable and even superior to those of Pt/C in both alkaline and acidic media without significant contribution by the metal sites. Such an N-doped carbon (denoted as N-HPCs) features a hollow and hierarchically porous architecture, and more importantly, a noncovalently bonded N-deficient/N-rich heterostructure providing the active sites for oxygen adsorption and activation owing to the efficient electron transfer between the layers. The primary Zn-air battery using N-HPCs as the cathode delivers a much higher power density of 158 mW cm^-2 , and the maximum power density in the H₂ -O₂ fuel cell reaches 486 mW cm^-2 , which is comparable to and even better than those using conventional Fe-N-C catalysts at cathodes.

Descriptors and graphical construction for in silico design of efficient and selective single atom catalysts for the eNRR

Outline a screening protocol that uses density functional theory calculations to simultaneously optimize with respect to multiple criteria, thereby successfully identifying catalysts that are highly selective and also result in low overpotentials for ammonia production through eNRR.

Highly efficient catalysts of ruthenium clusters on Fe3O4/MWCNTs for the hydrogen evolution reaction

Anchoring of Ru onto the Fe3O4/MWCNTs composite, its HER mechanism and DFT calculated change in Gibbs free energy of adsorbed hydrogen are described.

Water dissociation and association on mirror twin boundaries in two-dimensional MoSe2: insights from density functional theory calculations

The adsorption and dissociation of water molecules on two-dimensional transition metal dichalcogenides (TMDs) is expected to be dominated by point defects, such as vacancies, and edges. At the same time, the role of grain boundaries, and particularly, mirror twinboundaries (MTBs), whose concentration in TMDs can be quite high, is not fully understood. Using density functional theory calculations, we investigate the interaction of water, hydroxyl groups, as well as oxygen and hydrogen molecules with MoSe₂ monolayers when MTBs of various types are present. We show that the adsorption of all species on MTBs is energetically favorable as compared to that on the basal plane of pristine MoSe₂, but the interaction with Se vacancies is stronger. We further assess the energetics of various surface chemical reactions involving oxygen and hydrogen atoms. Our results indicate that water dissociation on the basal plane should be dominated by vacancies even when MTBs are present, but they facilitate water clustering through hydroxyl groups at MTBs, which can anchor water molecules and give rise to the decoration of MTBs with water clusters. Also, the presence of MTBs affects oxygen reduction reaction (ORR) on the MoSe₂ monolayer. Unlike Se vacancies which inhibit ORR due to a high overpotential, it is found that the ORR process on MTBs is more efficient, indicating their important role in the catalytic activity of MoSe₂ monolayer and likely other TMDs.

Metal-Free Triazine-Based 2D Covalent Organic Framework for Efficient H2 Evolution by Electrochemical Water Splitting

Hydrogen evolution reaction (HER) by electrochemical water splitting is one of the most active areas of energy research, yet the benchmark electrocatalysts used for this reaction are based on expensive noble metals. This is a major bottleneck for their large-scale operation. Thus, development of efficient metal-free electrocatalysts is of paramount importance for sustainable and economical production of the renewable fuel hydrogen by water splitting. Covalent organic frameworks (COFs) show much promise for this application by virtue of their architectural stability, nanoporosity, abundant active sites located periodically throughout the framework, and high electronic conductivity due to extended π-delocalization. This study concerns a new COF material, C₆ -TRZ-TFP, which is synthesized by solvothermal polycondensation of 2-hydroxybenzene-1,3,5-tricarbaldehyde (TFP) and 4,4',4''-(1,3,5-triazine-2,4,6-triyl)tris[(1,1'-biphenyl)-4-amine]. C₆ -TRZ-TFP displayed excellent HER activity in electrochemical water splitting, with a very low overpotential of 200 mV and specific activity of 0.2831 mA cm^-2 together with high retention of catalytic activity after a long duration of electrocatalysis in 0.5 m aqueous H₂ SO₄ . Density functional theory calculations suggest that the electron-deficient carbon sites near the π electron-donating nitrogen atoms are more active towards HER than those near the electron-withdrawing nitrogen and oxygen atoms.

Assessment of the Accuracy of Density Functionals for Calculating Oxygen Reduction Reaction on Nitrogen-Doped Graphene

Experimental studies of the oxygen reduction reaction (ORR) at nitrogen-doped graphene electrodes have reported a remarkably low overpotential, on the order of 0.5 V, similar to Pt-based electrodes. Theoretical calculations using density functional theory have lent support to this claim. However, other measurements have indicated that transition metal impurities are actually responsible for the ORR activity, thereby raising questions about the reliability of both the experiments and the calculations. To assess the accuracy of the theoretical calculations, various generalized gradient approximation (GGA), meta-GGA, and hybrid functionals are employed here and calibrated against high-level wave-function-based coupled-cluster calculations (CCSD(T)) of the overpotential as well as self-interaction corrected density functional calculations and published quantum Monte Carlo calculations of O adatom binding to graphene. The PBE0 and HSE06 hybrid functionals are found to give more accurate results than the GGA and meta-GGA functionals, as would be expected, and for a low dopant concentration, 3.1%, the overpotential is calculated to be 1.0 V. The GGA and meta-GGA functionals give a lower estimate by as much as 0.4 V. When the dopant concentration is doubled, the overpotential calculated with hybrid functionals decreases, while it increases in GGA functional calculations. The opposite trends result from different potential-determining steps, the *OOH species being of central importance in the hybrid functional calculations, while the reduction of *O determines the overpotential obtained in GGA and meta-GGA calculations. The results presented here are mainly based on calculations of periodic representations of the system, but a comparison is also made with molecular flake models that are found to give erratic results due to finite size effects and geometric distortions during energy minimization. The presence of the electrolyte has not been taken into account explicitly in the calculations presented here but is estimated to be important for definitive calculations of the overpotential.

Main Descriptors To Correlate Structures with the Performances of Electrocatalysts

Traditional trial and error approaches to search for hydrogen/oxygen redox catalysts with high activity and stability are typically tedious and inefficient. There is an urgent need to identify the most important parameters that determine the catalytic performance and so enable the development of design strategies for catalysts. In the past decades, several descriptors have been developed to unravel structure-performance relationships. This Minireview summarizes reactivity descriptors in electrocatalysis including adsorption energy descriptors involving reaction intermediates, electronic descriptors represented by a d-band center, structural descriptors, and universal descriptors, and discusses their merits/limitations. Understanding the trends in electrocatalytic performance and predicting promising catalytic materials using reactivity descriptors should enable the rational construction of catalysts. Artificial intelligence and machine learning have also been adopted to discover new and advanced descriptors. Finally, linear scaling relationships are analyzed and several strategies proposed to circumvent the established scaling relationships and overcome the constraints imposed on the catalytic performance.

Molecular Dynamics and Machine Learning in Catalysts

Given the importance of catalysts in the chemical industry, they have been extensively investigated by experimental and numerical methods. With the development of computational algorithms and computer hardware, large-scale simulations have enabled influential studies with more atomic details reflecting microscopic mechanisms. This review provides a comprehensive summary of recent developments in molecular dynamics, including ab initio molecular dynamics and reaction force-field molecular dynamics. Recent research on both approaches to catalyst calculations is reviewed, including growth, dehydrogenation, hydrogenation, oxidation reactions, bias, and recombination of carbon materials that can guide catalyst calculations. Machine learning has attracted increasing interest in recent years, and its combination with the field of catalysts has inspired promising development approaches. Its applications in machine learning potential, catalyst design, performance prediction, structure optimization, and classification have been summarized in detail. This review hopes to shed light and perspective on ML approaches in catalysts.

Fast prediction of oxygen reduction reaction activity on carbon nanotubes with a localized geometric descriptor

By using the pyramidalization angle as a localized geometric descriptor for oxygen reduction reaction (ORR) activity of carbon nanotubes (CNTs), we show the ORR activity of these systems can be readily predicted with mere structural optimization of CNTs.

Contrasting Oxygen Reduction Reactions on Zero- and One-Dimensional Defects of MoS2 for Versatile Applications

Oxygen reduction reaction (ORR) is a key microscopic process in many electrochemical applications of materials, where the requirements of their ORR performances may vary strikingly, for example, during the uses of MoS₂ as an electrocatalyst and anticorrosion/lubricating coating in aqueous/humid environments, ORR should be activated and inhibited, respectively. To reveal a complete ORR profile of MoS₂, using first-principles calculations, we examine the stabilities of various possible zero-dimensional point defects on the surface and one-dimensional edge defects and comprehensively explore the ORR activities on pristine MoS₂ surface and those defects in acid/alkaline solutions. It is found that the ORRs on the pristine surface and surfaces with point defects always require large overpotentials (>1.0 V), indicating a defect-immune resistance of the planar MoS₂ surface against the ORR. However, the ORR overpotentials on edge defects can reach as low as 0.66 V, corresponding to a relatively high activity close to that of the prototypical catalyst Pt (overpotential ∼0.45 V). Such contrasting ORR behaviors of point and edge defects are also understood in depth by analyzing the underlying thermodynamic and electronic-structure mechanisms. This work not only quantitatively explains the performances of MoS₂ in both galvanic corrosion and electrochemical catalysis but also provides a useful structure-ORR map that can facilitate adapting the realistic MoS₂ to versatile electrochemical applications.

Catalytic trends of nitrogen doped carbon nanotubes for oxygen reduction reaction

Replacing the state-of-the-art fuel cell catalyst platinum for a cheaper and abundant alternative would make the hydrogen economy viable. Both nitrogen-doped graphene and nitrogen-doped carbon nanotubes (N-CNT) have been shown to be capable of acting as a metal-free catalyst for the oxygen reduction reaction (ORR). Until now, most of the research has been focused on the nitrogen doping and less on the structure of the nanotubes. Here, density functional theory calculations are used to calculate trends in ORR catalytic activity of graphitic-N-doped CNTs of different sizes and chirality of selected tubes between (4,0) and (20,10). This includes 13 armchair tubes, 17 zig-zag tubes and 42 chiral tubes, or 72 N-CNTs in total. 22 tubes are predicted to have a lower overpotential than the platinum catalyst and 46 tubes have lower overpotential than nitrogen doped graphene. The most active tubes are (14,7), (12,6), and (8,8), and display an overpotential of around 0.35 V, or 0.1 V lower overpotential than predicted on Pt(111) with the same level of theory.

Guiding Principles for Designing Highly Efficient Metal-Free Carbon Catalysts

Carbon nanomaterials are promising metal-free catalysts for energy conversion and storage, but the catalysts are usually developed via traditional trial-and-error methods. To rationally design and accelerate the search for the highly efficient catalysts, it is necessary to establish design principles for the carbon-based catalysts. Here, theoretical analysis and material design of metal-free carbon nanomaterials as efficient photo-/electrocatalysts to facilitate the critical chemical reactions in clean and sustainable energy technologies are reviewed. These reactions include the oxygen reduction reaction in fuel cells, the oxygen evolution reaction in metal-air batteries, the iodine reduction reaction in dye-sensitized solar cells, the hydrogen evolution reaction in water splitting, and the carbon dioxide reduction in artificial photosynthesis. Basic catalytic principles, computationally guided design approaches and intrinsic descriptors, catalytic material design strategies, and future directions are discussed for the rational design and synthesis of highly efficient carbon-based catalysts for clean energy technologies.

Zigzag carbon as efficient and stable oxygen reduction electrocatalyst for proton exchange membrane fuel cells

Non-precious-metal or metal-free catalysts with stability are desirable but challenging for proton exchange membrane fuel cells. Here we partially unzip a multiwall carbon nanotube to synthesize zigzag-edged graphene nanoribbons with a carbon nanotube backbone for electrocatalysis of oxygen reduction in proton exchange membrane fuel cells. Zigzag carbon exhibits a peak areal power density of 0.161 W cm⁻² and a peak mass power density of 520 W g⁻¹, superior to most non-precious-metal electrocatalysts. Notably, the stability of zigzag carbon is improved in comparison with a representative iron-nitrogen-carbon catalyst in a fuel cell with hydrogen/oxygen gases at 0.5 V. Density functional theory calculation coupled with experimentation reveal that a zigzag carbon atom is the most active site for oxygen reduction among several types of carbon defects on graphene nanoribbons in acid electrolyte. This work demonstrates that zigzag carbon is a promising electrocatalyst for low-cost and durable proton exchange membrane fuel cells.

Cost and stability of catalysts hinder widespread use of proton exchange membrane fuel cells. Here the authors synthesize zigzag-edged graphene nanoribbons for electrocatalysis of oxygen reduction. Employment of such a metal-free catalyst in a fuel cell yields remarkable power density and durability.