Yimeng Sun
a,
Lin Tao
*a,
Yaqiong Su
*b,
Davoud Dastanc,
Han Zhang
a,
Hongwei Zhaoa,
Lixiang Li
a and
Baigang An
*a
aSchool of Chemical Engineering, University of Science and Technology Liaoning, Anshan 114051, China. E-mail: taolin@ustl.edu.cn; bgan@ustl.edu.cn
bSchool of Chemistry, Engineering Research Center of Energy Storage Materials and Devices of Ministry of Education, National Innovation Platform (Center) for Industry-Education Integration of Energy Storage Technology, Xi'an Jiaotong University, Xi'an, 710049, China. E-mail: yqsu1989@xjtu.edu.cn
cDepartment of Materials Science and Engineering, Cornell University, Ithaca, NY 14850, USA
First published on 27th January 2026
In recent years, the electrocatalytic carbon dioxide reduction reaction (CO2RR), driven by renewable energy and operated under mild conditions with controllable reaction pathways, has emerged as a promising route for carbon-neutral energy conversion. Two-dimensional supported catalysts have attracted particular interest owing to their tunable electronic structures and well-defined active sites. However, how the number and spatial configuration of active centers govern CO2RR activity and selectivity remains insufficiently understood, limiting the rational design of efficient catalysts. This review provides a comprehensive overview of recent experimental and theoretical advances in two-dimensional supported catalysts for the CO2RR, including single-atom (SACs), double-atom (DACs), and three-atom (TACs) catalysts, and metal clusters, with an emphasis on insights obtained from density functional theory (DFT). The fundamental reaction pathways of the CO2RR are first summarized, highlighting structure–activity relationships between active-site characteristics and catalytic performance. Subsequently, the advantages and limitations of different catalyst architectures are critically compared, and the mechanisms of CO2 reduction to C1 products such as CO, HCOOH, and CH4 are systematically analyzed. Particular attention is given to the role of DFT in elucidating reaction pathways, charge transfer, and adsorption energetics, thereby revealing key descriptors governing activity and selectivity. By integrating experimental observations with theoretical insights, this review aims to provide a mechanistic framework and design principles for the development of advanced two-dimensional supported catalysts for efficient CO2RRs.
![]() | ||
| Fig. 1 The spatial distribution of (a) the posteriori emissions, (b) the absolute discrepancies between the a posteriori and MEIC inventories (Mg km−2 year−1) for total annual CO2 emissions across China, and (c) a provincial-level comparison (TgCO2 per year) of the two emission datasets. Reproduced with permission from ref. 3, Copyright 2024 by American Chemical Society. | ||
The utilization of CO2 serves as a fundamental prerequisite for establishing a sustainable artificial carbon cycle. As a key feedstock in C1 chemistry, CO2 can currently be converted through three primary pathways. The first involves reforming CO2 with CH4 to produce syngas.6 The second entails transforming CO2 into fuels via electrochemical, photochemical, or thermochemical reactions.7 The third pathway focuses on the production of high-value chemicals through catalytic hydrogenation.8 CO2 transformation not only reduces fossil fuel consumption but also effectively alleviates the environmental impact of CO2 generated from human production and daily life,9 making it one of the important approaches to tackle excessive CO2 emissions. However, the CO2 molecule is linearly symmetrical in shape with the carbon atom in the highest oxidation state. The C–O bond length is only 0.116 nm and exhibits certain triple bond characteristics, and has a bond energy as high as 750 kJ mol−1, which indicates that breaking the C
O double bond requires great energy.10 In general, CO2 molecules exhibit chemical stability, thermodynamic stability, and kinetic inertness, making their activation and transformation highly challenging.11,12 Additionally, the CO2RR involves complex processes entailing multiple electron transfer steps, hydrogenation, and C–C bond coupling, which result in CO2 conversion with high energy consumption and low efficiency.
Light, electric, and thermal energy act as the requisite energy sources in the photocatalytic, electrocatalytic, and thermal catalytic reductions of CO2, respectively.13 Among these approaches, the thermal catalytic CO2RR is capable of generating a wide range of products, including various liquid fuels. However, it typically requires elevated temperatures and pressures and often depends on fossil fuel-derived energy sources.14 In contrast, photocatalysis and electrocatalysis can harness solar energy and electrical energy to drive the transformation of CO2 into C1 or C2+ products.15,16 Consequently, photocatalytic and electrocatalytic CO2RRs represent key pathways for achieving the sustainable utilization of new energy and are important strategies for attaining carbon neutrality. Notably, the CO2RR stands out as a pivotal technology for transforming CO2 into valuable energy products such as CO, HCOOH, and other C1 compounds, owing to its efficiency and controllability.17–22
The core function of an electrocatalyst lies in reducing the activation overpotential of the reaction. In addition, the active sites of catalysts can optimize the reaction process and decrease the activation energy by opening up an alternative reaction pathway. That is, catalysts can adsorb CO2 molecules and generate new intermediates. By controlling the adsorption behavior of the above new intermediates and weakening the interaction force between them, a more favorable condition for the reaction is obtained. Nevertheless, the diverse range of CO2RR products (including CO, CH4, CH3OH, C2H4, C2H5OH, etc.) and their similar reaction potentials23 tend to result in the formation of mixed products. Furthermore, the hydrogen evolution reaction (HER) is a major competitive side reaction in the CO2RR, leading to extra energy consumption and reduction of reaction efficiency. Therefore, exploring novel, high-performance CO2RR electrocatalysts with low overpotential and the ability to inhibit the HER is of great significance.
The unique quantum confinement effects of two-dimensional (2D) materials have positioned them as a rapidly emerging and highly regarded platform in catalytic research, adjustable electronic structures, abundant active sites, large specific surface areas, and cost-effectiveness.24–31 In the selection of carriers for atomic dispersion catalysts, apart from two-dimensional materials, porous crystal materials such as metal–organic frameworks (MOFs) and zeolites have also attracted extensive attention due to their regular and ordered pore structures, extremely high specific surface areas, and precisely modifiable coordination sites. For instance, MOFs provide an ideal space for anchoring and confining metal precursors, which is conducive to the preparation of SACs with high loading capacity.32–34 The rigid framework and ion-exchange sites of zeolites can effectively stabilize metal cations and prevent their migration and agglomeration.35,36 These carrier systems have demonstrated unique potential in the CO2RR. Graphene, with its two-dimensional layered architecture37–39 and remarkable properties, has sustained global focus. This material exhibits a high visible-light transmittance of 97.7%, a heat conduction capability up to 3 × 103 W m−1 K−1 at 298 K, an electrical conduction of approximately 104 Ω−1 cm−1, along with a specific surface area of 2630 m2 g−1 as well as a Young's modulus of 1.1 TPa.40,41 Owing to these properties, it demonstrates significant potential for applications spanning energy (e.g., storage, solar cells), electronics (e.g., transparent electrodes), aerospace, and advanced materials (e.g., composites, biomedicine, catalysis, HER),42–45 holding the potential to replace numerous traditional materials. However, its zero-bandgap nature restricts its application in the field of logic electronic devices.46,47 Beyond graphene, a multitude of two-dimensional materials have achieved successful prediction and laboratory-based synthesis in recent years. Of the hundreds of existing 2D materials, there are single-atom48–52 and double-atom53–57 catalysts where metal atoms serve as active centers embedded in a 2D substrate, as well as three-atom58–62 and metal cluster catalysts,63–67 transition metal-supported two-dimensional carbon nitrides,68–72 and so on. The combination of high charge carrier mobility, efficient heat dissipation, and mechanical stability in these standalone materials makes them promising for high-performance CO2RR electrocatalysis.
Since the 20th century, as research into material structures has deepened, traditional research methods have struggled to meet specific investigative demands. However, the rapid advancements in computing technologies and quantum physical theories have enabled the computer-based simulation of material characteristics. Researchers have successively created a variety of computational methods and theoretical frameworks, including quantum mechanical first-principles, molecular dynamics, and machine learning.73,74 Among these, quantum mechanical first-principles simulate material properties based on the atomic spatial arrangement and electronic structures.75 The scale of the systems it can handle typically ranges from several to hundreds of atoms, and it exhibits remarkable strengths in calculating electronic structures and optomagnetic properties. Quantum chemistry addresses chemical problems by applying the classical theories of quantum mechanics, allowing the acquisition of data including molecular orbital energies and the total energy of chemical systems.76 Density functional theory (DFT), a major breakthrough during the evolution of electronic structure theory, establishes a correspondence between the ground-state wave function and the ground-state electron density.77,78 This approach significantly reduces computational complexity while maintaining accuracy, proving capable of conducting in-depth studies on material reaction mechanisms thus enabling the rational design of materials and screening.
However, the reduction reaction of CO2 still faces several fundamental challenges at the mechanism level in theoretical calculation research. Firstly, the initial activation path of the reaction remains controversial, and the adsorption strength of the intermediate has a significant impact on the selectivity of the main product. The identification of key intermediates (such as COOH and OCHO) and their formation sequence are difficult to define clearly. The reason why different metal and non-metal catalysts can generate a variety of CO2RR products is mainly because they have unique adsorption capacity and reactivity for intermediates. Secondly, the carbon–carbon coupling mechanism for generating C2+ products has not been fully elucidated. Multiple pathways such as CO dimerization and Co–*CHO coupling are relatively close on the reaction energy barrier, leading to disputes over the rate-determining steps of the reaction and the sources of product selectivity.79,80 The essence of these problems lies in the fact that the current theoretical models have limitations in accurately describing the real electrochemical interface environment, including the synergistic effects of multiple factors such as solvation and interface electric fields, which have not been fully quantified. Furthermore, the dynamic structural evolution of the catalyst surface at the working potential leads to significant differences between the calculation results based on the ideal static model and the structure and properties of the actual active sites, thereby severely restricting the ability to precisely predict reaction pathways and product selectivity at the atomic scale. To effectively overcome the energy barrier in the electron-transfer-proton coupling process, accelerate the catalytic reaction rate and suppress the occurrence of unwanted side reaction pathways, it is crucial to develop and design ideal catalysts with high selectivity and low energy consumption.
This review summarizes the recent research progress in CO2RR catalysts for converting CO2 into one-carbon products. The catalysts include two-dimensional materials such as single-atom, diatomic, and triatomic catalysts, and metal cluster catalysts embedded in two-dimensional substrates with metal active sites, and two-dimensional carbon nitrides supported by transition metals. It reveals the structure–activity relationship between the material properties and performance. In addition, it not only elaborates on the advantages and disadvantages of different catalysts and details the mechanism of CO2 reduction to one-carbon products but also reviews the relevant content of density functional theory (DFT) calculations and electrocatalysis. This review aims to provide a basis for more suitable catalysts and accurate mechanisms for the CO2RR.
Compared with traditional three-dimensional porous carbon or oxide carriers, two-dimensional carriers (such as graphene, carbon nitride, MXene, etc.) have clear surface geometric structures and tunable electronic properties, providing an ideal platform for constructing and stabilizing specific coordination structures of SACs.94–97 A large number of studies have shown that by precisely regulating the coordination number, types of coordination atoms (such as N, O, S, and Cl), and coordination configurations of the metal center on the two-dimensional carrier, the physical and chemical state can be effectively optimized, thereby directing the improvements in CO2RR performance.
Specifically, coordination number regulation is a key means to adjust the d-band center of the catalyst and change the adsorption strength of key intermediates (such as COOH). For example, compared with the common M–N4 configuration, low coordination number Co–N2 sites have been proven to effectively reduce the formation energy barrier of COOH, achieving up to 94% CO selectivity under appropriate overpotential.98 Coordination atom engineering involves introducing different electronegative heteroatoms to adjust the charge distribution of the metal center. For example, pyrrole N species can stabilize the Fe3+ valence state and optimize the reaction pathway;99 while axial coordination of Cl (such as Fe–N4Cl) can significantly change the adsorption behavior of intermediates.100 Moreover, the defects and pore engineering of the two-dimensional carrier (such as constructing nitrogen vacancy Ni–N3–V structure) not only provide abundant anchoring sites for metal atoms but also enhance the mass transfer of CO2 and electrolyte, raising the current density to industrial relevant levels (such as >100 mA cm−2).101
Further strategies include constructing double-atom sites (such as Ni–Fe and Co–Cu) to achieve dual-center cooperative catalysis, and using tandem catalytic design (coupling SACs with Cu-based catalysts) to expand the product range from C1 (CO, HCOOH) to high-value C2+ products (ethylene, ethanol).102–104 It is worth noting that advanced in situ characterization techniques have revealed that SACs may undergo dynamic structural evolution under reaction conditions (such as Fe3+/Fe2+ valence state changes), which is closely related to their final activity and selectivity.105 More cutting-edge research even utilizes atomic-level-resolved in situ electrochemical STM to directly observe the dynamic migration and reorganization process of Cu single-atom sites at the reaction interface and correlate it with the selectivity switch of products.106 In summary, SAC systems based on two-dimensional carriers achieve precise control of the CO2RR pathway through multi-dimensional and cross-scale structural design, demonstrating great potential from fundamental understanding to practical applications. Since Zhang's research group published their pioneering work in 2011,107 SACs have developed in a crucial direction for overcoming bottlenecks within the realm of heterogeneous catalysis. Their research demonstrated that the activity of the Pt1/FeOx SACs in the CO oxidation reaction exceeds that of nano-platinum catalysts by three times.108,109 With an atomic architecture where metal centers are coordinated by non-metal atoms, this type of catalyst offers unique advantages by maximizing the exposure of active sites. Its distinctive coordination structure and electronic properties can significantly accelerate the kinetic process of the CO2RR, thus exhibiting great potential in this field.110–113 Notably, the agglomeration of single atoms on the substrate continues to be a challenging issue in the preparation of SACs, as it can diminish the catalyst's performance. Therefore, a suitable support is required to form strong bonds with single atoms, preventing metal aggregation. Typical supports include chalcogenides, oxides, and carbon materials.114 While metal oxide supports can enhance CO2RR activity through strong interactions with single atoms,115–117 the oxygen species on their surfaces may cause the oxidation of metal single atoms, thereby reducing catalytic efficiency. In contrast, two-dimensional carbon-centric materials possess considerable application prospects, owing to their outstanding stability in reaction environments.118–122
With its unique structure, two-dimensional graphitic carbon nitride (g-C3N4) has proven effective for anchoring Cu, Pd, and Pt single atoms, serving as an excellent SAC support. This prominence stems from its cost-effectiveness, robust physicochemical stability, and facile synthesis protocols.123–126 Sun and colleagues69 elucidated a general mechanistic principle governing g-C3N4-catalyzed CO2 conversion (Fig. 2a), non-spontaneous processes are characterized by positive reaction energies, and such processes arise when CO2 is adsorbed on the surface or undergoes transformation via an odd count of H+/e− pairs to form free radicals. Conversely, spontaneous processes with negative reaction energies are associated with the release of products involving even H+/e− pairs, including CO, H2CO, and CH3OH. Transition metal dichalcogenides (TMDs) have also garnered growing interest, with their asymmetric layered structures127–129 conferring unique physical properties that theoretically and experimentally enhance CO2RR efficiency.130–134
![]() | ||
| Fig. 2 (a) Optimized configurations of planar and corrugated 2 × 2 g-C3N4 sheets, showing their resulting electronic DOS. Reproduced with permission from ref. 69, Copyright 2016 by the Royal Society of Chemistry. (b) Schematic of the TM@J-WSSe structure with 28 transition metal (TM) atoms anchored onto the substrate. Reproduced with permission from ref. 136, Copyright 2025 by Elsevier. (c) Stable structures of Ge-N4G and Tl-N4G. Reproduced with permission from ref. 137, Copyright 2025 by American Chemical Society. (d) Schematic of first-row transition metal single atoms anchored on a g-C3N4 monolayer for the CO2RR. Reproduced with permission from ref. 68, Copyright 2025 by American Chemical Society. (e) Optimized atomic structures of the investigated catalyst series. Reproduced with permission from ref. 70, Copyright 2024 by American Chemical Society. (f) Schematic diagram of the structure of M–N4@Gr (M = Fe/Co/Ni). Reproduced with permission from ref. 110, Copyright 2022 by American Chemical Society. (g) Top and side views of the optimized geometric structures for the Fe–Nx–C catalyst series (models I–IV). Reproduced with permission from ref. 140, Copyright 2023 by American Chemical Society. | ||
The Janus WSSe (J-WSSe) monolayer, a novel two-dimensional TMD, is synthesized via selective epitaxial atomic replacement (SEAR), where one Se layer in WSe2 is substituted with S.135 This architecture—with W bonded to S on one side and Se on the other—breaks inversion symmetry. Transition metal J-WSSe SACs with anchored (TM) show outstanding electrocatalytic activity, yet the potential of J-WSSe with TM anchoring and Se incorporation vacancies as high-performance CO2RR catalysts warrants systematic investigation. Chen et al.136 explored 2d 3d, 4d, and 5d transition metals anchored on J-WSSe via anchoring, designing and screening a library of SACs (Fig. 2b). Their findings identified Os@J-WSSe, Cu@J-WSSe, Mn@J-WSSe, and Cr@J-WSSe as optimal catalysts for CO, HCOOH, CH3OH, and CH4 production, respectively. Notably, Cr@J-WSSe displayed unique dynamic stability, maintaining consistent active sites pre- and post-reaction. Tan and coworkers137 reported CO2RR performance of p-block metal SACs supported on N-doped graphene (M–N4G, where M = In, Tl, Ge, Sn, Pb, Sb, Bi. Fig. 2c). Experimental results showed most M–N4G catalysts exhibit superior CO2RR selectivity over the competing HER. As a key intermediate species, *HCOO's adsorption energy has stood out as a trustworthy indicator to evaluate activity. Among these, Ge/Pb–N4G demonstrated exceptional electrocatalytic efficacy for HCOOH synthesis, offering insights for rational SAC design.
Anna and colleagues68 studied the stability of first-row transition metal atoms anchored on g-C3N4 monolayers in the CO2RR. As displayed in Fig. 2(d), Sc, Ti, V, and Zn have large proton affinity (ΔG(H*) < −0.3 eV), which will result in active site blockage and possibly inhibit CO2 conversion. In contrast, Cr, Mn, Fe, Co, and Ni have moderately positive energies of proton binding (ΔG(H*) = 0.1–0.4 eV), enabling suppression of the competing HER. Among systems with relatively strong binding energies showing repulsion (Ni, Mn, Co) Ni1/C3N4 stands out, as its moderate *CO adsorption energy (−0.52 eV), optimized d-band center position (−2.3 eV vs. Fermi level), and charge modulation of the metal by coordinated nitrogen sites (Δq = 0.53 e−) achieve an optimal balance between overpotential (ηCO2RR = 0.38 V) and selectivity (78% CO Faraday efficiency). Liu et al.70 systematically investigated CO2RR over metal-modified C3N4 nanotubes (Mn/CNNTs) in Fig. 2e for CO and formic acid production. Detailed analysis of geometric and electronic structures revealed that CNNTs’ unique adsorption sites efficiently activate CO2 and preferentially direct the CO2RR toward the CO pathway, enhancing selectivity. Single-atom-modified armchair CNNTs (M/aN-6) exhibit the highest CO selectivity in the CO2RR.
Furthermore, extensive research confirms that the local coordination environment dictates the electrocatalytic activity and selectivity of SACs, serving as a cornerstone for efficient electrocatalysis. Typically, SACs’ active centers adopt coordination frameworks, with catalytic efficiency hinging on the central metal, coordination atom types/numbers, support properties, and functional group modifications.71,72 Even subtle structural perturbations can markedly alter SACs’ performance.138 In M–N–C systems, for instance, modifying the carbon support morphology and generating nearby defects serves to both expose active sites and regulate the metal's electronic states, leading to enhanced catalytic efficacy.139 Yang et al.110 confirmed that the coordination environment modulates d-orbital electron occupancy in M–N4@Gr (M = Fe/Co/Ni) SACs (Fig. 2f), inducing distinct electron spin polarization that influences catalytic selectivity and performance. Fe–N4@Gr exhibits significant spin polarization, with decreasing polarization following the sequence of Fe > Co > Ni.
In the deformable coordination field, metal–ligand bond distortion alters d-orbital splitting, regulating intermediate adsorption strength. Ni's dx2−y2 and dxy orbitals bond with N 2p orbitals, leaving few empty orbitals to accept electrons and resulting in weak CO2 adsorption. In contrast, Fe's dz2/dzx orbitals and Co's d orbitals remain partially occupied, endowing Fe/Co single atoms with both electron-donor and acceptor capabilities. This work demonstrates optimized catalytic performance via precise regulation of the coordination environment around non-precious metals, underscoring its significance for SACs. Wang et al.140 elucidated the structure–activity relationship of Fe–NxCy–C (x = 2–4, y = 0–2) catalysts in the CO2RR in Fig. 2g, that spin characteristics of Fe sites can induce modulation of the nearby coordination environment and density of active sites. The Fe–N3P–C structure can sustain high CO yields over a wide potential range. This study demonstrates the intrinsic structure–performance relationship revealed by the reaction itself, and suggests that the synergism between coordination environment and the density of active sites is the key to understanding and enhancing the performance of SACs.
Although SACs have great potential, the homogeneity of their active sites still requires in-depth research, and achieving their controllable synthesis under different coordination environments remains challenging. Importantly, current SAC-derived CO2RR products are predominantly CO, with limited reports of CH4, HCOOH, etc. There is a paucity of research on generating high-value C2+ hydrocarbon fuels by tuning single-atom-support interactions, emphasizing the need to develop new SAC materials for highly selective hydrocarbon synthesis.
As illustrated in Fig. 3(a), Liu et al.152 investigated CeO2(110)-supported diatomic catalysts (Mα–Mβ/CeO2, where Mα–Mβ are 3d/4d/5d transition metals). Analysis of binding and aggregation energies revealed that metals with lower electronegativity and higher d-orbital occupancy stabilize supported DACs. High-throughput descriptor calculations identified a series of potential DACs with both high activity (turnover frequency, TOF > 0.1 s−1) and selectivity (>50%). Among these, heterogeneous Au-based DACs (e.g., Au–Fe) exhibited good performance. The introduction of Au further activated adsorbed oxygen and weakened CO binding. This theoretical study establishes a foundational framework for the rational design of dual-site catalysts. Current DAC research remains in its infancy. DACs are classified as homonuclear or heteronuclear based on the identity of the two atoms in the active center.
![]() | ||
| Fig. 3 (a) Atomic structure of CeO2 (1 × 1 × 1) with Ce (green) and O (pink) atoms, and schematic search diagram of DACs/CeO2(110) for CO PROX. Reproduced with permission from ref. 152, Copyright 2025 by American Chemical Society. (b) DFT study of bimetallic Ni–Fe–N–C catalysts, structural configurations (Ni–N4, Fe–N4, FeNi–NC, FeNi–NSC with C-gray, N-blue, S-yellow, Ni-green, Fe-cyan) and their catalytic reaction energetics, including free energy profiles and stepwise energy comparisons. Reproduced with permission from ref. 153, Copyright 2024 by American Chemical Society. (c) Computational characterization of bimetallic TM1/TM2-N@Gra systems showing geometric configurations, thermodynamic stability, and CO2 adsorption geometry. Reproduced with permission from ref. 57, Copyright 2022 by Elsevier. (d) Structural and electronic properties of pristine g-CN and Fe2-doped g-CN systems. Illustrated are the optimized geometries (top and side views), TDOS, and comparative binding energy analysis. Reproduced with permission from ref. 154, Copyright 2021 by the Royal Society of Chemistry. (e) Proposed reaction pathway for CO2 to C2H5OH, C2H4 on N-graphene based dual-atom catalysts, geometric configurations and overall reaction network encompassing both C1 (CH3OH, CH4) and C2 products, and two distinct reaction pathways from the *CO intermediate to C2-coupled species. Reproduced with permission from ref. 155, Copyright 2020 by the Royal Society of Chemistry. (f) DFT study of Ni-based SACs for the CO2RR, showing COOH* adsorption configurations on NiMN6 and NiN4 sites. Reproduced with permission from ref. 56, Copyright 2022 by American Chemical Society. | ||
As shown in Fig. 3(b), Huang's group153 demonstrated that modulation via diatomic catalysts or heteroatomic ligands is a promising strategy to optimize SACs for more efficient CO2 conversion to high-value chemicals. Their study focused on electrocatalytic CO2 reduction to CO using an N/S-coordinated Fe–Ni diatomic pair embedded in a carbon matrix (FeNi–NSC). Gibbs free energy analysis of CO2RR steps revealed that S coordination weakens Fe site activity and collaborates with adjacent Ni sites to regulate bidentate adsorption of the COO* intermediate, significantly reducing the CO desorption barrier and altering the rate-determining step. Wei et al.57 proposed that the catalytic property of bimetallic atomic sites may be enhanced by regulating the binding mode of crucial intermediates supported by bonding interaction and electronic structure analyses. As shown in Fig. 3(c), their study screened four typical transition metals (TM = Fe, Co, Ni, Cu) and then constructed ten TM1/TM2-N@Gra to systematically assess homonuclear and heteronuclear DACs for CO2 to C1 products. Among them, Fe/Co and CO2 displayed optimal limiting potentials (UL) for HCOOH formation (−0.37 V and −0.34 V, correspondingly), indicating the best catalytic performance. The enhanced activity for CO2 reduction was attributed to the optimal adsorption energy of the HCOOH intermediate on Fe/Co sites, as revealed by electronic structure analysis.
Bui's team154 proposed a new strategy to boost CO2RR performance by integrating well-defined structural clusters. As shown in Fig. 3(d), they designed g-CN molecular skeletons hosting various transition metal (TM) dimers for the CO2RR, ultimately identifying five atomic pairs (Cr2, CrFe, Mn2, MnFe, Fe2) with excellent active performance. The performance of Fe2@g-CN in the CO2RR was further optimized by regulating the integration of iron atomic clusters. It is worth noting that the introduction of Fe13 atomic clusters enhanced the CO2 adsorption capacity while effectively inhibiting hydrogen activation. In addition, it also breaks the scale relationship between the intermediates (CO and CHO), thereby significantly enhancing the performance of CO2RR and reducing the limiting potential of the C1 path to −0.45V. This work broadens understanding of DAC mechanisms in the CO2RR and offers insights for rational design and integration of heterogeneous DACs with specific atomic clusters for other applications. Chen et al.155 screened bimetallic dimers supported on nitrogen-doped graphene (including homonuclear and heteronuclear transition metals). As shown in Fig. 3(e), C–C coupling analysis revealed that six catalysts (Cr–Cu, Mn–Cu, Co–Co, Co–Ni, Co–Cu, and Ni–Cu) showed coupling energy barriers under the threshold (0.75 eV), with the top five showing notably lower barriers. All six DACs demonstrated superior C2 pathway selectivity in competitive HER and CO2RR.
Despite atomic vibrations, the bimetallic atoms and their six coordinated nitrogen atoms maintained strong bonding without obvious structural distortion, indicating excellent heat stability at 500 K. This suggests these DACs are able to stably retain their configuration under both mild and high-temperature conditions. As displayed in Fig. 3(f), Zhu et al.56 studied CO2 activation based on the Ni–Cu catalyst. As shown, the secondary Cu could tune the energy level of the Ni 3d orbital into the Fermi level and then enhance the emergence of the rate-determining *COOH intermediate. Remarkably, a CO TOF of 20
695 h−1 was obtained at −0.6 V (vs RHE) with the maximum CO Faraday efficiency of 97.7%. In situ XAFS monitoring of dynamic structural evolution further clarified interactions between Ni centers and CO2 molecules, along with the Ni–Cu synergistic mechanism enhancing CO2RR activity.
Constructing structures where one metal atom binds exclusively to another (especially in dimer form) remains a major challenge. Inhomogeneity in support pores and defects hinders precise fabrication of diatomic sites. In contrast, nano synthetic catalysts or SACs outperform DACs in terms of controllability and performance.156 Excessive metal precursor usage can form metal nanoparticles, while reducing dosage may promote single-atom dispersion and isolation. However, DACs exist in an intermediate state, making precise control of precursor dosage difficult and often leading to uneven metal distribution in catalytic centers. Due to imprecise metal loading control, coexistence of single atoms, atomic clusters, and even larger aggregates is inevitable, triggering active site aggregation. When catalysts contain ill-defined clusters or particles, clarifying the structure–performance relationship of dual-site metal catalysts becomes challenging. Most current DACs use carbon-based supports, with relatively few studies concerning other substrates. Thus, the feasibility and performance of DACs on non-carbon supports require further exploration.
Jia et al.163 reported that TACs behave as high current electrocatalysts for the CO2RR because the large number of active sites in TACs benefit from both the high electrochemical surface area offered by the nanocubes and fast charge migration. The electronic properties and stability of TACs are also related to their carrier materials. By anchoring the three atoms on appropriate carriers and taking advantage of their interfacial interactions, the agglomeration of clusters can be effectively avoided,164 and an optimized triatomic coordination environment can effectively enhance catalytic performance.165 For example, when three metal atoms are atomically dispersed on supports such as nitrogen-doped graphene or g-C3N4, TACs with either random166–168 or triangular169,170 structures can form. Such catalysts actively participate in electrochemical reactions—including N–N bond/O–O bond cleavage or C–C coupling—via orbital coupling or non-bonding interactions between metal atoms.171
Zhou's group60 explored the utility of graphdiyne-supported 3TM-GY (TM = Cu, Fe, Co) applied in the CO2RR to generate C1 products, as shown in Fig. 4(a). The singular triangular TM trimer anchored on graphdiyne, characterized by strong adsorption energy and a distinct curved configuration, creates a beneficial setting for CO2 adsorption and initial activation. 3Cu-GY preferentially reduces CO2 to CH4, exhibiting higher selectivity than other C1 products and surpassing most atomically dispersed electrocatalysts. At increased applied potentials, the CO2RR on 3Cu-GY favors high-throughput production of HCOOH and CH4, while 3Fe-GY/3Co-GY favor deeply reduced products. Hydrogen evolution is strongly inhibited, rendering CO2RR dominant. This performance enhancement is attributed to synergies among multiple metal centers, which promote adsorption of critical species, notably HCOO* and CH* and facilitate pathways conducive to CO2 reduction.172 Sun et al.173 compared SACs and TACs (Fig. 4(b)), and showed that SACs maximize atomic utilization and provide uniform active sites, while TACs precisely regulate intermediate adsorption via interatomic synergies, overcoming traditional catalytic limitations. Among their findings, 3Mo-C2N1 and 3Ti-C2N1 exhibited exceptional catalytic activity for CO and HCOOH production, respectively. Swapnil et al.174 investigated CO2 adsorption via TMn-doped C2N monolayers, revealing that TM3-C2N achieves a maximum adsorption capacity of 6 CO2 molecules per trimer—significantly outperforming TM1-C2N and TM2-C2N systems.
![]() | ||
| Fig. 4 (a) Optimized atomic structures and charge density difference distributions of the 3TM-GY catalysts. Reproduced with permission from ref. 60, Copyright 2021 by Elsevier. (b) A schematic diagram showing the stability and conductivity of TM-C2N1. Reproduced with permission from ref. 173, Copyright 2024 by the Royal Society of Chemistry. (c) Design strategy and synthesis implementation. Reproduced with permission from ref. 176, Copyright 2024 by American Chemical Society. (d) Top-view structure of the trimetallic TM1TM2TM3@GY catalyst. Reproduced with permission from ref. 58, Copyright 2024 by Elsevier. (e) Optimized CO* adsorption configurations on TM-NiFe catalysts, with corresponding Mulliken charges and projected density of states (PDOS) of the transition metal (TM) atoms. Reproduced with permission from ref. 175, Copyright 2023 by Elsevier. (f) Structural models of SACs, DACs, and TACs on Mo2CO2 MXene, with their corresponding transition metal constituents, formation energies, and dissolution potentials. Reproduced with permission from ref. 172, Copyright 2023 by the Royal Society of Chemistry. | ||
Copper demonstrates unique efficacy in catalyzing the CO2 reduction reaction toward hydrocarbon products, owing to its unique ability to adsorb *CO intermediates. However, traditional Cu catalysts suffer from deficiencies in CH4 selectivity and activity. As shown in Fig. 4(c), Pan's team59 addressed this atomic-level engineering of Cu sites and mesoscale control of support porosity, enabling dynamic transformation of single-atom Cu–N3 into N- and OH−-anchored Cu3 clusters under reaction conditions. The N, OH-Cu3 site stabilizes *CO at −0.76 eV. The kinetic barrier for its hydrogenation to *CHO is 0.85 eV, markedly lower than the C–C coupling barrier (2.03 eV), suggesting that formate formation is kinetically preferential, selectively promoting CH4 formation. Its moderate adsorption of CO capacity and low barriers for CO hydrogenation are key to high selectivity. Additionally, dual continuous mesoporous channels create a water-deficient local microenvironment, inhibiting hydrogen evolution and prolonging *CO intermediate retention, thereby promoting deep CO2 reduction.
Similarly, as presented in Fig. 4(d), Li et al.58 developed a series of high-performance non-noble triatomic catalysts (TM1TM2TM3@GY, where TM = Mn, Fe, Co, Ni, Cu, Mo) supported on the 2D material graphdiyne, investigating their CO2RR selectivity. MnMoFe@GY and MnMoCu@GY displayed optimal CO2RR performance, achieving the highest activity and selectivity among all tested catalysts, with effective hydrogen evolution inhibition. They demonstrated good experimental feasibility, stable catalytic performance, and excellent CO2RR selectivity—confirming TACs’ potential in CO2RR catalysis. Liu et al.175 designed a series of trimetallic TM-NiFe catalysts as presented in Fig. 4(e), finding that Zn–NiFe performs exceptionally in the CO2RR-to-CO pathway. Cr–NiFe prefers the COOH* pathway over HCOO* for HCOOH production, enabling efficient HCOOH synthesis while significantly inhibiting competing reactions (hydrogen evolution and CO2RR-to-CO). These results highlight the abundance of CO2RR active sites in TACs. Xiao et al.172 observed that homonuclear diatomic and triatomic catalysts supported on 2D Mo2CO2 materials exhibit superior CO2RR performance compared to single-atom systems. As shown in Fig. 4(f), multinuclear reaction sites in multi-atom catalysts significantly enhance adsorption of crucial intermediates (e.g., HCOO* and CH*), which facilitates the selective reduction to CH4 at ultra-low overpotentials. Additionally, multinuclear sites effectively promote C–C coupling, creating conditions for efficient C2H5OH production.
Benefiting from multi-metal synergy, high selectivity, and cost-effectiveness, TACs have attracted significant research interest in CO2 electroreduction. However, rational TAC design faces fundamental challenges. Future efforts must address precise synthesis, stability enhancement, and industrial adaptation. Integrating advanced characterization and computational simulation will be critical to advancing TACs from laboratory research to large-scale applications, contributing to global carbon neutrality goals.
![]() | ||
| Fig. 5 (a) Material classification by size scale. Reproduced with permission from ref. 176, Copyright 2024 by the Royal Society of Chemistry. (b) Schematic of the synthesis of Fe–NC–SBA-15. Reproduced with permission from ref. 185, Copyright 2020 by the Royal Society of Chemistry. (c) Structural model of the CoNiCuRh@FeN4 high-entropy catalyst (ΔSmix = 0.14 meV K−1), and the correlation between its dissolution potential (Udiss) and formation energy (Eformation). Reproduced with permission from ref. 187, Copyright 2025 by Elsevier. (d) Molecular structures of three copper-bromide-phosphine complexes: [Cu3(μ-Br)2(bis-methyl-bisphenylphosphine)3] Br-(CuBr–BisM), Cu2(μ-Br)2(1,2-phenyl-bisphenylphosphine)2 (CuBr–12B) and Cu2(μ-Br)2(triphenylphosphine)2(4-phenylpyridine)2 (CuBr–4PP). Reproduced with permission from ref. 188, Copyright 2024 by Nature. (e) Most stable CO2 adsorption configurations, charge density difference isosurfaces, and projected density of states (PDOS) for CO2 on the FT Cu3(Se)@MoSSe, FT Cu3(S)@MoSSe, FT Cu3(Te)@MoSTe, and FT Cu3(S)@MoSTe catalysts. Reproduced with permission from ref. 196, Copyright 2024 by American Chemical Society. (f) Theoretical simulations of the CO/HCOOH formation pathways, electronic structure analysis (charge density and PDOS of d-/p-bands) of the ZnO, Ag, and ZnO–Ag systems, and the corresponding HER processes. Reproduced with permission from ref. 200, Copyright 2021 by American Chemical Society. | ||
As illustrated in Fig. 5(b), Wu et al.186 demonstrated that iron cluster sites anchored on mesoporous carbon materials exhibit enhanced electrocatalytic performance in the CO2RR. Similarly, a study on CuAg bimetallic clusters187 utilized DFT to analyze relationships between catalyst structure, composition, and activity, revealing that specific cluster configurations favor the formation of target products such as CH3OH. These findings underscore the critical role of atomic arrangement within clusters in optimizing CO2RR efficiency. Tamtaji's group188 reported five configurations comprising five metals Fe, Co, Ni, Cu, and Rh. In Fig. 5(c), CoNiCuRh@FeN4, FeN4 serves as the active center, while the CoNiCuRh high-entropy cluster serves as a regulatory unit, achieving a low overpotential for CO production in the CO2RR. In FeNiCuRh@CoN4, CoN4 synergizes with the FeNiCuRh high-entropy cluster as the active site, enabling selective reduction of CO2 to methanol at an overpotential—outperforming traditional SACs. Fig. 5(d) highlights that binuclear Cu(I) complexes serve as efficient catalysts for CO2-to-C3 product conversion via C–C coupling, indicating that multinuclear arrangements can promote C2+ formation by offering a favorable microenvironment for CO2 activation and subsequent coupling steps.189
DFT-based explorations of alloy and bimetallic clusters are also prominent. Mechanistic analyses indicate that C–C coupling is influenced by the catalyst's ability to manage competing reactions (e.g., hydrogen evolution). Multinuclear catalysts can regulate these kinetics to favor CO2 reduction pathways. For example, a study on AuCu clusters190 employed first-principles calculations to simulate catalytic behavior, emphasizing how alloying modulates adsorption energies and reaction pathways. Similarly, investigations of Cu@CNT catalysts191 established that certain Ni-based clusters preferentially convert CO2 to CO, demonstrating the utility of computational methods in predicting catalytic performance from cluster composition. Studies on defects and near-surface atomic structures have further elucidated their impact on catalytic activity. Zhang et al.192 used DFT to analyze defective Cu4Xn clusters, finding that defects alter adsorption sites and lower activation barriers, thereby enhancing the CO2RR. Jo et al.193 explored surface FCC dependence, using DFT to clarify how different Cu surface facets influence HCOOH formation pathways. They found that electronic properties related to surface coordination and roughness significantly affect catalytic efficiency, emphasizing the critical role of surface engineering for nanoparticle catalysts. Xu et al.194 synthesized CuO nanosheets with controllable nanostructures as effective CO2RR catalysts, noting that surface oxygen and Cu+ species are not essential for C2+ product formation, suggesting catalytic activity is primarily determined by metal surface structure and morphology. Deng et al.195 further supported size-dependent activity, combining theoretical and experimental methods to identify an optimal silver nanoparticle size (8–10 nm) balancing activity and selectivity—emphasizing the need for precise size control in CO2RR nanoparticle catalysts.
Fig. 5e demonstrates the outstanding CO2RR performance of the FT-Cu3@MoSX (X = Se, Te) catalyst. CO2 is stably adsorbed on the catalyst surface, and the charge transfer (0.59–0.97e) and orbital hybridization are significant, converting linear CO2 into V-shaped CO2δ− active species (OCO bond angle 121.7°–136.8°, C
O bond length increasing to 1.30 Å). Moreover, the electron transfer efficiency at the S-terminated surface is better, with an adsorption energy (−0.51 to −0.55 eV) superior to that of the non-polar MoS2 system. This effectively activates CO2 and lowers the reaction energy barrier, demonstrating its high catalytic potential.196 Structural engineering of atomic catalysts—including multinuclear arrangements—has been shown to provide multiple reaction pathways and synergistic effects, enhancing overall electrocatalytic activity.188 Innovative catalyst design strategies include Zhang et al.'s197 Zn–Ag–O catalytic materials in Fig. 5f, featuring nanoscale nanoparticles within nanopores to achieve high performance and durability, illustrating how nanoscale structural design boosts CO2RR efficiency. Gang et al.198 explored the scalability and practicality of nanoparticle catalysts, developing a one-step pyrolysis method to prepare M–N–C catalysts synthesized from carbon nanotubes and nitrogen precursors—an economical, effective approach enabling scalable CO2RR catalysts and expanding practical application potential.
It is worth noting that, based on multi-metal cluster catalysts, high-entropy alloy clusters (HEA Clusters) are emerging as an emerging frontier. HEA clusters are typically composed of five or more principal metal elements in a nearly equimolar ratio, and their extremely high configurational entropy is expected to confer them unique electronic structures and outstanding structural stability. For instance, as mentioned earlier in this review, the CoNiCuRh@FeN4 system188,199 already possesses characteristics of being multi-component and having high entropy, demonstrating its unique control ability over different product pathways in the CO2RR. Theoretically, the vast component space of HEA clusters contains infinite possibilities for electronic structures, providing an unprecedented platform for 'tailoring' CO2RR catalysts with high activity, high selectivity, and ultra-high stability. However, the controllable synthesis, precise characterization, and elucidation of the relationship between structure and function remain significant challenges. In the future, by leveraging high-throughput computing and machine learning, the optimal combinations from the vast ‘high-entropy chemical space’ will be selected, and in situ characterization techniques for sub-nanometer-scale high-entropy systems will be developed, which will be the key to achieving breakthroughs in this field.
Theoretical and spectroscopic techniques, such as multinuclear magnetic resonance, have been used to study hyperpolarized nuclear spin behavior in multinuclear systems, informing the design of more effective catalysts. Additionally, integrating multinuclear complexes into molecular junctions to enhance thermoelectric performance highlights the broader applicability of multinuclear structures in energy-related catalysis. In summary, these studies collectively emphasize the significance of multinuclear metal complexes in advancing CO2RR technology. Incorporating multinuclear structures into catalyst design provides a universal strategy to overcome current CO2RR limitations, enabling more efficient, selective, and stable catalytic systems. Sustained mechanistic research and structural innovation are critical to fully harnessing the viability of multinuclear catalysts for sustainable energy CO2 conversion technologies.
The mechanism for CO2-to-CO conversion involves a process of transferring 2 electrons (2e−) and 2 protons (2H+). The CO2 reduction process commences with the efficient adsorption of CO2 molecules on the surface (Fig. 6a). The adsorbed CO2 is reduced to the *COOH intermediate. Subsequent attack by an additional H+ and e− converts *COOH to H2O and *CO, with the final step involving *CO desorption from the electrode surface to release CO.202 Current research has confirmed that the two key bottlenecks restricting the production efficiency of carbon monoxide lie in the fact that too weak COO− adsorption energy hinders the formation of COOH intermediates, while too strong *CO adsorption inhibits the desorption of the products.204,205 Multiple studies corroborate that the *CO binding energy on Cu surfaces acts as a key activity descriptor to characterize CO2RR products beyond 2e− pathways, aligning with the Sabatier principle. Overly strong *CO binding results in catalyst poisoning, whereas weak binding causes premature CO desorption, preventing further reduction. Optimal *CO adsorption balances high CO uptake with HER suppression, thereby maintaining high CO2RR Faraday efficiency.
![]() | ||
| Fig. 6 (a) Proposed reaction pathways and key intermediates for the CO2RR to CO products. Reproduced with permission from ref. 201, Copyright 2025 by Elsevier. (b) Depicted are multiple aspects of the catalytic system, including the model structures, computed Gibbs free energy profiles for CO2-to-CO conversion, the disparity between CO2RR and HER limiting potentials, potential-dependent in situ ATR-IR spectra in CO2-saturated 0.5 M KHCO3, and the proposed CO2RR mechanism on Ni-SAs@FNC. Reproduced with permission from ref. 99, Copyright 2021 by Elsevier. (c) Gibbs free energy diagram for the electrochemical CO2RR process. Reproduced with permission from ref. 207, Copyright 2022 by American Chemical Society. (d) XAS characterization and theoretical computation of the single-atom Fe catalyst. Reproduced with permission from ref. 208, Copyright 2022 by American Chemical Society. | ||
The coordinative environment and electronic structure are pivotal in determining the catalytic efficacy. Wang's group212 developed a catalyst featuring bismuth atomic sites co-coordinated by nitrogen and sulfur. Synthesized via a cation–anion diffusion strategy, this catalyst exhibits enhanced electrocatalytic activity.213 Li et al.214 explored low-valent Znδ+ single-atom sites, which demonstrated high current density and promising industrial application potential. Moreover, Song et al.215 introduced a boron-doped Ni–N4 SAC, where boron atoms collaborate with nickel to lower the energy barrier of *COOH generation, thereby realizing excellent performance in CO2 reduction. Guo et al.216 revealed a correlation between the tailored atomic environment of SnN3O1 sites and enhanced transformation of CO2 into CO via regulating the intermediates’ binding energy, while simultaneously inhibiting competitive pathways such as HCOOH generation. In a related vein, Chaipraditgul et al.217 investigated transition metal modifications on alumina supports to modulate surface interactions, though their research extended to CO2 hydrogenation for alkene production.
Sun et al.221 showed that regulating the microenvironment of Ni–Fe binary SACs using ionic liquids can further enhance CO2 electroreduction, indicating that microenvironment engineering is pivotal for the performance of diatomic catalysts. The synergistic effects observed in these diatomic systems are supported by in situ spectroscopic and theoretical studies, which confirm that such configurations can reduce activation energy and promote the selective generation of CO. Additionally, Sun's team173 conducted a comparative study between single-atom and multi-atom catalysts. By designing transition metal catalysts (TM-C2N1) supported on single-layer nitrogen-doped graphene, they achieved CO2 reduction to CO using both single-atom and multi-atom catalyst systems. Gibbs free energy and electronic structure studies on 3TM-C2N1 indicated their superior catalytic performance for CO during the CO2RR, with relatively low limiting potentials, which offer a novel approach to promoting sustainable CO2 conversion.
The Huang team222 reported the Fe–Se bimetallic single-atom catalyst FeSe–NC, prepared by a selenic acid etching-assisted method, achieves efficient CO2-to-CO conversion through the synergistic effect of Fe–N4 and Se–C2 sites. The core mechanism was characterized and verified by Fig. 7a and theoretical calculations. In situ ATR-FTIR showed that the introduction of Se atoms significantly weakened the adsorption strength of *CO at the Fe sites and accelerated the desorption of *CO. KIE analysis indicated an improvement in proton transfer kinetics. DFT calculations confirmed that the Fe-Se synergy not only lowered the formation energy barrier of *COOH but also optimized the thermodynamics of *CO desorption. Moreover, the negative shift of the d-band center of Fe (−0.61 to −0.88 eV) and the increase in the difference in the limit potential between the CO2RR/HER enhanced the inhibition of the HER side reaction. Combined with the mass transfer advantage of the hierarchical porous structure, the catalyst was stably generated as CO at high FECO (up to 97.7%) and industrial-grade jCO (228 mA cm−2).
![]() | ||
| Fig. 7 (a) In situ ATR-FTIR spectra, CO partial current densities in H2O/D2O, kinetic isotope effects, and theoretical calculations comparing the CO2 reduction reaction performance and electronic structure of Fe–NC and FeSe–NC catalysts. Reproduced with permission from ref. 222, Copyright 2024 by American Chemical Society. (b) Performance evaluation of the three ML algorithms. Reproduced with permission from ref. 224, Copyright 2023 by American Chemical Society. (c) Gibbs free energy profiles for CO2-to-CO and the HER on various structural models, a schematic of the in situ ATR-SEIRAS setup, and the corresponding spectra for the ZnNi-1000, Zn-SAC, and Ni-SAC catalysts. Reproduced with permission from ref. 225, Copyright 2024 by Elsevier. (d) Electrochemical CO2 reduction performance in a GDE flow cell and a zero-gap MEA electrolyzer using a NiFe-DAC catalyst. Metrics include CO faradaic efficiency (FECO), partial current density, and comparative performance analysis against state-of-the-art catalysts across multiple scales. Reproduced with permission from ref. 223. Copyright 2024 by Elsevier. | ||
Similarly, Han et al.223 investigated that the high CO selectivity of NiFe-DAC is attributed to the cooperative electronic regulation of Fe–N5and Ni–N4Ni–N4 sites. As shown in Fig. 7d, the Fe–N5 site preferentially adsorbs activated CO2 to form the *COOH intermediate (with a binding energy similar to that of Fe-SAC), and the adjacent Ni–N4 site weakens the Fe's adsorption of *CO through electron coupling (the desorption free energy decreases from 1.1 eV to 0.1 eV), avoiding site poisoning and accelerating the desorption of *CO, moreover, the Ni site has a stronger adsorption of H2O, inhibiting its own adsorption of CO2, thus enabling the CO2 to selectively accumulate at the Fe site.
Xiong et al.224 utilized machine learning (support vector regression algorithm) to predict that the overpotential of this catalyst is 0.11 V. As shown in Fig. 7b, the decisive step for the reduction of CO2 to CO is CO2 → COOH*, and its overpotential is only 0.09 V. The results of the two are highly consistent, confirming the reliability of the catalytic activity. Electrochemical tests showed that at an electrode potential of −0.8 V vs RHE, the CO faradaic efficiency of this catalyst is as high as 82.12%, significantly higher than that of pure Cu (57.04%), Pd@Cu (79.25%), and Pt@Cu (71.24%). Meanwhile, the faradaic efficiency of H2 is lower than those of the above control catalysts. This effectively inhibits the HER. In the stability experiment, after continuous electrolysis in the H-cell for 12 hours, the current density and CO faradaic efficiency did not show significant decline, and XRD characterization confirmed that its crystal structure remained stable with no significant changes. This verified the core value of PdPt@Cu as an efficient catalyst for CO2 reduction to CO.
Furthermore, studies have shown that the performance of the ZnNi-TACs catalyst in converting CO2 to CO is superior to that of SACs.225 As confirmed by Fig. 7c, the free energy curve of the CO2RR for its core triatomic structure (Zn–Ni1–Ni2, Ni1–Zn–Ni2) is significantly lower than that of SACs, significantly reducing the activation energy barrier for the formation of COOH, which is more favorable for subsequent CO hydrogenation and desorption thermodynamics, and effectively inhibiting the HER. The in situ ATR-SEIRAS spectroscopy further captured the evolution of key intermediates. The peaks at 3365 cm−1 and 1623 cm−1 confirmed the adsorption and dissociation of H2O on the Zn atoms, the peak at 3706 cm−1 and the peak at 1482 cm−1 indicated that *H was stabilized by the N atom; the intensity of the characteristic peaks at 1396 cm−1 (*COOH) and 1933 cm−1 (*CO) increased with the negative shift of the potential, and the peak positions and intensities were higher than those of Zn-SAC and Ni-SAC, confirming that the triatomic structure accelerates the formation and transformation of intermediates. This enables efficient and highly selective CO generation.
:
1) samples, which enable CO2-to-CO conversion via the CH4 + CO2/CO2 redox cycle. The incorporation of Ni enhances the material's activity in the temperature range of 873–973 K. During the reductive step CH4 reacts with CO2 on Ni sites to diminish the oxygen storage component, while in the oxidation step, the oxygen storage component is further oxidized by CO2 to form CO. Similarly, Zhanaidarova et al.231 investigated the anchoring of Re(tBu-bpy)(CO)3Cl on MWCNTs, which increases the current density, lowers the overpotential, and maintains CO selectivity under aqueous conditions. This highlights that polyatomic catalysts supported on conductive substrates can boost CO2 electroreduction efficiency.
![]() | ||
| Fig. 8 (a) Electrochemical performance of the CO2RR on the catalyst structure. Reproduced with permission from ref. 226, Copyright 2022 by Wiley. (b) Depicted are the predicted *CO adsorption energies and DOS at different sites, the optimized adsorption geometries of key intermediates on N, OH-Cu3, and the free energy profiles (at 0 V vs. RHE) for CO2 reduction to CH4 and for *COCO coupling on the same active site. Reproduced with permission from ref. 59, Copyright 2024 by American Chemical Society. (c) Energy profiles along the pathway for photochemical CO2-to-CO conversion. Reproduced with permission from ref. 227, Copyright 2021 by Elsevier. (d) The catalyst library of Pd–Fe/Co/Ni/Cu bimetallic systems with different doping ratios, applied in the ethane/CO2 reaction. Reproduced with permission from ref. 228, Copyright 2021 by Elsevier. | ||
The role of metal/oxide interactions in regulating catalytic activity was further exemplified by Huo et al.,232 who studied Cu/SnOx heterostructure CNT-supported nanoparticles. Their systematic analysis revealed that component-dependent metal-oxide interplay significantly affects CO2 electroreduction, suggesting that multi-atomic configurations can optimize the selectivity toward target products. Nanostructured Au catalysts are also employed in electrochemical CO2 reduction. Pan et al.59 obtained dynamic reconfiguration of N,OH–Cu3 clusters originated from atomically single-metal sites. Fig. 8b shows Cu–N3 sites by regulating the atomic-level structure of Cu active sites and designing the engineered macro-morphology of the carbon support. Benefiting from the moderate CO adsorption ability and low CO hydrogenation energy barrier, the N,OH–Cu3 site attains an unprecedented Faraday efficiency of 74.2% for the CO2-to-CH4 conversion at an industrial-grade current density of 300 mA cm−2. Nursanto et al.233 found that the CO selectivity exhibited a positive dependence on Au loading, plateauing at 78% as the morphology evolved from clustered to layered structures. The superior mass activity of 4 nm nanoparticles definitively confirmed that the nanoscale architecture dictates the CO production reactivity and selectivity.
Alfonso et al.234 further analyzed ligand-protected Au25 clusters, noting that fully ligand-protected clusters lack efficiency for CO2 reduction due to the high electrocatalytic potential required to form the key carboxyl intermediate, emphasizing the importance of surface and ligand engineering in cluster catalysts. Au25 nanoclusters reported by Wu's team235 are equipped with thiol ligands, one hand to maintain the Au25(SR)18 structure, but cover all Au sites to CO adsorption and make the intact catalyst dead, while the other to detach from the Au25(SR)18/CeO2 interface at ≥423 K to activate CO and induce low-temperature CO oxidation on the Au25(SR)18/CeO2 rod-shaped catalyst. Xiao et al.236 highlighted the significance of metal–ligand cooperation in dinickel complexes, showing that bimetallic synergy significantly enhances catalytic reactivity—their NiIINiII(bphpp)(AcO)2 complex exhibits 5-fold higher reactivity than mononuclear analogs.
As indicated in Fig. 8(c), Xu et al.227 found epitaxially grown Co9S8 embedded with CoO, which strengthens CO2 adsorption affinity and activates the molecule via bond length extension, facilitating reduction to CO. This underscores the value of engineering multi-atomic structures to promote CO2 activation. Additionally, doping the catalyst matrix can regulate activity and selectivity. Han's team237 reported a high-performance atomic Fe–In–NC catalyst for CO2 reduction, which achieved a 95% CO faradaic efficiency and surpassed Fe–NC across a broad potential range. This performance was attributed to the Fe–In d–p orbital hybridization, which concurrently promotes *CO desorption and reduces the *COOH formation barrier. Li et al.228 explored Pd doping in Fe/Ni/CeO2 catalysts (Fig. 8(d)), where Ni6/CeO2 achieved 96.5% CO selectivity in ethane-CO2 dry reforming, demonstrating excellent performance. However, Pd doping in Ni-based catalysts reduces syngas selectivity due to enhanced cleavage of C–C bonds and increased activation of CO2 energy barriers, which is unfavorable for CO generation. Fe-based and Pd-doped Fe-based catalysts preferentially produce ethylene, with minimal impact on CO.
In summary, the production of *COOH is considered to be the rate-determining step in the conversion of CO2-to-CO.238 Therefore, the final selectivity and energy yield of this reaction are greatly limited by the *COOH intermediate adsorption on the effective sites of the surface. Furthermore, it is reasonable to construct a strong atomic local electric field by heteroatom doping to improve the *COOH adsorption.
To date, extensive efforts have focused on developing catalysts for HCOOH/HCOO− production in the CO2RR, primarily involving p-block metals.241 This 2-electron transfer process has a complex reaction mechanism. After decades of research, the potential pathways have been identified. CO2 is adsorbed at the active sites of the catalyst, undergoes activation and reduction processes to form intermediates such as *HCOO, *COOH or *OCHO, and eventually generates HCOOH/HCOO−.
C double bond vacancies of C3N monolayers—alongside TM-N-C (TM-N co-doped carbon) materials. M-CC catalysts facilitate electron transport in electrocatalysis, with all M-CCs showing superior selectivity for the CO2RR over the HER. Notably, Cu-CC exhibits a low UL of 0.68 V toward HCOOH production. Song et al.244 explored SACs embedded in antimonene monolayers for the CO2RR, finding that 3d non-precious TMs anchored on antimonene exhibit superior CO2RR selectivity over the HER. Among these, Zn-based SACs uniquely produce HCOOH, while others primarily generate CH4. Co@antimonene exhibits a low overpotential (0.50 V), which rivals that of state-of-the-art electrocatalysts. The interaction between TMs and antimonene modulates SACs intrinsic activity, while TM binding to potential-determining step (PDS) intermediates dictates CO2RR overpotential and product selectivity.The He team245 investigated the single-atom Cu adsorbed Janus MoSSe monolayer catalyst (MoSSe-Cu). Fig. 9a clearly presents the complete reaction pathway of CO2 → HCOOH. The rate-determining step is *OCHO → *HCOOH, and the corresponding free energy is only 0.49 eV, which is significantly lower than that of g-C3N4 and C2N monolayer-loaded transition metal trimer catalysts (0.86 eV, 0.57 eV). The thermodynamic feasibility is better. From the free energy change in Fig. 9a, it can be seen that the step of *CO2 combining with the proton-electron pair to form *OCHO has a free energy reduction of 0.78 eV, which is a spontaneous reaction, while the free energy for generating the competing intermediate *COOH is only reduced by 0.10 eV. Thermodynamically, the *OCHO path is dominant, ensuring the high selectivity of HCOOH. Chen et al.246 found that introducing sulfur atoms into the NiN4 coordination structure specifically formed NiN2S2 (Fig. 9(b)). Their work highlights the critical role of heteroatom coordination in modulating the performance of Ni SACs, as the sulfur dopants regulate the electronic properties environment of Ni active centers to favor *OCHO adsorption and subsequent formic acid formation.
![]() | ||
| Fig. 9 (a) The reaction pathways and Gibbs free energy of CO2 reduction to HCOOH on MoSSe-Cu. Reproduced with permission from ref. 245, Copyright 2023 by Elsevier. (b) Gibbs free energy profiles for the CO2RR pathways to CO and HCOOH, along with the adsorption configurations of key intermediates on NiN4, NiN3S1, and NiN2S2. Reproduced with permission from ref. 246, Copyright 2023 by the Royal Society of Chemistry. (c) The full catalytic cycle depicting the reduction of two CO2 molecules by one equivalent of Fe@B10H14, computed at the B3LYP-D3 level. Reproduced with permission from ref. 250, Copyright 2017 by the Royal Society of Chemistry. (d) Potential reaction pathways from *HCOOH and comparative limiting potentials for CO2-to-HCOOH conversion on single-atom transition metal catalysts supported on TMDs. Reproduced with permission from ref. 249, Copyright 2025 by American Chemical Society. | ||
Similarly, Zeng's team247 explored MOF-derived 2D SACs based on transition metal-tetra hydroxybenzoquinone frameworks. These materials offer abundant accessible active sites and versatile functionalization capabilities, enabling precise tuning for specific reaction pathways—including selective reduction of CO2. The MOF-based platform offers a flexible avenue for optimizing catalytic selectivity through structural design. Qian et al.248 demonstrated that doping iron into a boron cage (B10H14) yields an innovative catalyst (Fe@B10H14) capable of catalyzing CO2 hydrogenation by means of a two-step reduction pathway. As presented in Fig. 9(c), their quantum mechanical studies revealed that this catalyst can sequentially reduce two CO2 molecules to produce HCOOH, highlighting the potential of metal-doped boron-based frameworks in CO2 reduction. Further research into nickel-based SACs has also been conducted to explore their catalytic performance in formic acid generation, expanding the portfolio of promising catalyst systems for this reaction. As shown in Fig. 9d, Li et al.249 reported the further protonation behavior of the HCOOH intermediate on the surface of Mn@NbS2 indicates that the energy barriers required for its conversion to HOCHOH and *OCH intermediates are as high as 2.05 eV and 1.91 eV respectively, significantly higher than the energy barrier of 0.53 eV required for direct desorption to form HCOOH. Thermodynamically, this effectively inhibits the formation of by-products such as CH4 and CH3OH, ensuring the high selectivity of the HCOOH product. In terms of catalytic activity, Fig. 9d compares the limiting potential of Mn@NbS2 with those of previously reported single-atom transition metal anchored two-dimensional transition metal chalcogenide (TMD) catalysts such as Ti@WTe2, Fe@MoTe2, and Ni@MoTe2. This shows that the limiting potential of Mn@NbS2 is as low as −0.53 V, which is within the optimal potential range for efficient catalytic reactions (absolute value <0.6 V), and only requires a lower applied driving voltage to achieve efficient reaction.
(from chemisorbed to physiosorbed states) drives the reaction to *OCHO
(rather than *COOH) intermediates as a result of asymmetric charge buildup at bimetallic centers, and correlates this charge buildup with a proposed charge aggregation intensity (CAI) reactivity descriptor. Screening catalysts according to the CAI descriptor identifies the NiSb diatomic catalyst as very promising for HCOOH conversion. During formic acid synthesis, as shown in Fig. 10(b), the Cr-NiFe catalyst exhibits a preference for the COOH* pathway rather than the HCOO* pathway, with surface Cr atoms' negative charge and low d-band center enabling efficient formic acid production at an exceptional overpotential as low as 0.080 V while markedly inhibiting competing reactions.
![]() | ||
| Fig. 10 (a) Calculated binding energies of NiSn DACs in various configurations, pCOHP between Ni and Sn atoms, reaction pathways for *CO2 protonation to COOH/OCHO, potential-dependent total energies of adsorption configurations, and free energy profiles for the CO2-to-CO reduction process. Reproduced with permission from ref. 226, Copyright 2025 by American Chemical Society. (b) Gibbs free energy comparisons for the CO2RR and HER, the correlation between ΔGh* and εd, the HER overpotential on TM-NiFe catalysts, and mechanistic analysis of CO/HCOOH selectivity via the *COOH intermediate. Reproduced with permission from ref. 175, Copyright 2023 by Elsevier. (c) Reaction pathways and corresponding Gibbs free energy diagrams for the HER and the conversion of CO2 to COOH/HCOO intermediates. Reproduced with permission from ref. 252, Copyright 2021 by Elsevier. (d) Correlation between the adsorption energy of products and their corresponding overpotential. Reproduced with permission from ref. 173, Copyright 2024 by the Royal Society of Chemistry. | ||
Similarly, Fig. 10(c) shows Yang et al.252 used DFT to screen M@2D-FeS2 catalysts, revealing Co@2D-FeS2 exhibits excellent activity and selectivity for HCOOH while strongly inhibiting the HER. Sun's team173 compared SACs and multi-atom catalysts (MACs), developing TM-C2N1 for the CO2RR. Gibbs free energy and electronic structure analyses on 3TM-C2N1 show favorable CO2RR catalytic performance with low limiting potentials. Notably, 3Ti-C2N1 achieves optimal HCOOH production corresponding to a UL of −0.42 V, positioned at the volcano plot vertex (Fig. 10(d)), offering a new route for sustainable CO2 conversion. Liu et al.175 designed trimetallic TM-NiFe catalysts, constructing 2D activity volcano curves to evaluate their performance in CO2-to-CO and HCOOH conversion. Additionally, Sun et al.253 engineered PtS2-based SACs with 3d TMs as candidate CO2RR electrocatalysts, systematically investigating pathways to C1 products. All TM-PtS2 SACs show higher selectivity for C1 products over the HER. For Sc, Ti, V, Cr, Mn, Fe, and Cu supported on PtS2-SV (sulfur-vacancy PtS2), HCOOH is the dominant product, while Co−, Ni−, and Zn–PtS2 generate diverse C1 products.
![]() | ||
| Fig. 11 (a) Comparative free energy diagrams for CO2 reduction to HCOOH and CO on the In (101) (red) versus CuO (111) (black) facets. Reproduced with permission from ref. 255, Copyright 2021 by Elsevier. (b) Displayed are the DFT-calculated ΔG for the CO2-to-formate pathways on both facet and edge sites of Bi (003) and (012) planes. A schematic diagram illustrates the mechanism for selective formate formation on Bi nanosheets. Reproduced with permission from ref. 256, Copyright 2018 by Elsevier. (c) The adsorption configurations of CO2, HCOO* and HCOOH in the Co–PbCO3 heterojunction and uniform metal doping, as well as the free energy diagram of the CO2RR on the surface. Reproduced with permission from ref. 238, Copyright 2023 by Elsevier. (d) The average current density obtained at different overpotentials, faradaic efficiencies, and the product selectivities of H2, HCOOH and CO generated at different overpotentials. Reproduced with permission from ref. 257, Copyright 2017 by American Chemical Society. (e) The faradaic efficiency, bias current density, specific current and conversion frequency of CO + HCOOH and H2 are presented. Reproduced with permission from ref. 258, Copyright 2021 by Elsevier. | ||
At the atomic scale, Jiang et al.256 observed that in CO2-saturated 0.1 M KHCO3 solution, Pd-B/C catalysts achieve a format Faraday efficiency (ηHCOO−) of 70% following 2 hours of electrolysis. Even at potentials more negative than the −0.5 V mark—where CO formation (ηCO) increases—the ηHCOO−/ηCO selectivity ratio remains significantly higher for Pd–B/C than Pd/C, indicating that subsurface B doping on Pd preferentially promotes adsorption and formation of the format-pathway intermediate *HCOO over the CO− pathway intermediate *COOH. Catalyst structural features, such as grain boundaries and nanostructures, also modulate CO2 reduction performance. Kumar et al.250 then demonstrated that SnO2 porous nanowires with diminished SnO2 grain boundaries, and thus increased active-site density, exhibit enhanced CO2-to-HCOOH conversion. This catalyst can initiate formic acid synthesis at a low overpotential (350 mV) and sustain a faradaic efficiency of 80% at −0.8 V vs. RHE, with higher energy conversion efficiency than similar systems.
Advances in catalyst design extend to ultrathin 2D nanostructures. Zhang et al.259 obtained ultrathin Bi nanosheets synthesized by liquid-phase exfoliation (Fig. 11(b)), where *OCHO intermediate formation preferentially occurs at edge sites (confirmed by lower Gibbs free energy) rather than basal planes. These nanosheets achieve an 86.0% formate Faraday efficiency at −1.1 V vs. RHE, where the current density reaches 16.5 mA cm−2, outperforming bulk Bi due to superior conductivity and plentiful edge sites. Atomically precise clusters offer tunable selectivity. Liu et al.260 found that double tetrahedral Cu8 clusters exhibit a formic acid Faraday efficiency of 92% at −1.0 V vs. RHE, surpassing cubic Cu8 clusters. Theoretical studies reveal weaker competition with hydrogen evolution and lower *HCOO adsorption free energy on the double tetrahedral structure, clarifying how intermediate interactions govern activity and selectivity. Zhang's group239 reported Co-modified PbCO3 electrocatalysts achieving a formate faradaic efficiency of 98.15% at −0.70 V vs. RHE (Fig. 11(c)), and this performance stems from reduced *HCOO formation barriers and enhanced adsorption at Co–Pb2+ sites, where Pb2+ charge enrichment boosts formate selectivity. Tang et al.257 highlighted the role of negatively charged hydrides in copper hydride nanoclusters (e.g., Cu32H20L12. Fig. 11(d)), which direct selectivity toward formic acid over CO at low overpotentials. DFT calculations further predict weaker hydrogen evolution competitiveness under these conditions. Wang et al.261 developed heterostructures of ultrasmall polymetallic sulfide clusters on MWCNTs (Fig. 11(e)), achieving >80% faradaic efficiency for CO2-to-C1 conversion (including formic acid) with a high turnover frequency (5974.62 h−1) and 12-hour stability.
Thus, designing efficient electrocatalysts remains central to advancing CO2RR for HCOOH production. Component engineering, defect engineering, and controlled morphology engineering are established strategies, but critical gaps persist including precise pore structure tuning for enhanced active surface area properties and mass transport capabilities. Future work still needs to address the issue of separating electrolytes from formic acid. The form of the product is closely related to the pH value of the system. Formic acid is formed under acidic conditions, while formate is formed under weakly alkaline conditions. Therefore, the separation process needs to be designed specifically based on this.258 Notably, suppressing the thermodynamically favored HER still stands as a key challenge, requiring synergistic optimization of catalyst and reactor design.
![]() | ||
| Fig. 12 (a) Structure and spin charge density analysis of Ru1/CeO2 and Ru2/CeO2 catalysts for selective CO2 hydrogenation, including adsorption configurations, reaction pathways, and calculated energy profiles. Reproduced with permission from ref. 262, Copyright 2025 by the Royal Society of Chemistry. (b) Free energy diagrams and corresponding stable intermediate structures for the electrocatalytic reduction of CO2 to CH4 on a Zn–N4-graphene catalyst. Reproduced with permission from ref. 263, Copyright 2020 by American Chemical Society. (c) Free energy for converting CO2 to CH4 via an 8-electron process on o-RuB2N2@TiN, p-RuB2N2@TiN and RuB3N1@TiN. Reproduced with permission from ref. 264, Copyright 2024 by Elsevier. | ||
Similarly, the Han team263 reported the SA-Zn/MNC single-atom catalyst, which achieves selective regulation of intermediates by the Zn-N4 active site (Zn is +2 valence), successfully completing the CO2 electroreduction to CH4. Fig. 12b clearly shows its 8-electron transfer free energy pathway. CO2 preferentially forms the OCHO intermediate at the Zn site (energy barrier 0.46 eV), rather than the key COOH (energy barrier 1.2 eV), blocking the competing pathways at the source. Fig. 12b also shows that OCHO forms bonds with the O atom through the O atom and carbon terminal protonation to form a stable configuration, with the rate-determining step being OCHOH → CHO (energy barrier 0.87 V). Each intermediate maintains stable coordination with Zn in Fig. 12b, and the OHCH3 protonation releases CH4, and the residual OH combines with H+ to form H2O, completing the cycle. In situ ATRSEIRAS further confirmed the existence of OCH2 and *O CH3 and no CO signal. Combined with the catalyst's microporous structure (525 m2 g−1, 1.04 nm) and the charge transport advantage of nitrogen-doped carbon support, SA-Zn/MNC achieved 85% CH4 faradaic efficiency at −1.8 V vs. SCE, −31.8 mA cm−2 partial current density, and stable operation for 35 h.
Pan et al.264 further optimized the active center through boron coordination, developing a TiN-supported Ru SAC (RuBxN4−x@TiN) for CO2 reduction to CH4 as another efficient solution. Among them, o-RuB2N2@TiN, p-Ru B2N2@TiN and Ru B3N1@TiN in Fig. 12c exhibited the best performance. Their formation energies were negative and the structure was stable at 500 K (the fluctuation of Ru–B bond length was only 0.02 Å). The B atom acted as a Lewis acid site, reducing the CO2 adsorption energy to −1.57 eV, extending the C–O bond to 1.3 Å, and the charge transfer amount was 0.52–0.64e−, achieving complete activation of CO2. Moreover, the rate-determining step energy barrier for CH4 generation was <0.7 eV, and the limiting potential was as low as 0.59–0.68 V, which was significantly superior to RuN4@TiN (1.28 V) and pure TiN (1.72 V).
The excellent activity and selectivity were demonstrated by the synergistic effect of the coordination environment and the carrier, which showed a key reference for the structural design and path optimization of SACs in multi-electron transfer CO2 reduction reactions.
![]() | ||
| Fig. 13 (a) Free energy profiles of the CO2RR to CH4 and corresponding adsorption configurations. Reproduced with permission from ref. 265, Copyright 2024 by Elsevier. (b) Possible C1 pathways of CO2 reduction on the Fe2Ir@NG, and the corresponding free energy diagrams. Data denote the ΔG of each elementary step. Reproduced with permission from ref. 266, Copyright 2022 by American Chemical Society. (c) Competition between the CO2RR and HER. Reproduced with permission from ref. 267, Copyright 2025 by Elsevier. | ||
The Fe2Ir@NG bimetallic catalyst developed by Han et al.266 exhibits distinct reaction pathway characteristics. Fig. 13b comprehensively depicts the key process of CO2 conversion to CH4 on its surface. As can be seen from the figure, after CO2 adsorption, a proton-electron pair transfer occurs to generate the CO + OH intermediate, which lowers the free energy by 0.74 eV. This step is a spontaneous reaction and does not require overcoming an energy barrier. Subsequently, CO gradually undergoes hydrogenation to form intermediates such as CHO (ΔG = 0.36 eV) and CHOH (ΔG = 0.49 eV), with the rate-determining step being CHO → CHOH. The corresponding limiting potential is only 0.49 V, significantly lower than that of the similar Fe3@NG catalyst (0.79 V). Comparing the free energy data of the CH3OH and CO generation pathways in the figure, the free energies of all intermediates in the CH4 generation pathway are at a lower level, and there are no high-energy barrier bottlenecks throughout the process. The thermodynamic advantage is significant, which is attributed to the synergistic effect of Fe and Ir atoms, optimizing the adsorption strength of key intermediates such as CO and *CHO, enabling the CH4 generation process to have a low limiting potential and high spontaneity characteristics.
The Ti3@NG catalyst studied by the Wang team267 exhibits extremely high structural stability, providing a stable active environment for catalytic reactions. The adsorption process of CO2 on its surface occurs spontaneously (with an adsorption free energy of negative value), and Ti3@NG transfers 1.69e of charge to the CO2 molecule, causing a significant reduction in the O–C–O bond angle and achieving efficient activation of the CO2 molecule. The reaction path calculation in Fig. 13c indicates that the Ti3@NG catalyst for CO2 reduction follows the optimal path of CO2* → COOH* → CO* → CHO* → HCHO* → H3CO* → CH4, with the rate-determining step being the process of combining CHO with proton-electron pairs to form HCHO. The corresponding limiting potential is only −0.53 V, which is superior to most reported catalysts for CH4 generation. The volcano curve analysis further reveals that the adsorption intensity of Ti3@NG for the key intermediate CO is moderate, precisely located at the summit active region of the CH4 generation reaction, avoiding the retention of intermediates due to excessive adsorption or the inhibition of protonation due to insufficient adsorption. At the same time, its adsorption free energy for CO2 is more negative than that of the H atom, and the limiting potential of the CO2RR is lower than that of the HER, effectively inhibiting the competition of side reactions.
The above study achieved efficient CO2 reduction to CH4 by designing dual-atom (Ru2-NC, Fe2Ir@NG) and metal cluster (Ti3@NG) catalysts, and optimizing the adsorption of intermediates and reaction pathways through the synergistic effect of metal sites. The core performance advantages were confirmed through free energy calculations, structural stability characterization, and selectivity analysis, providing important references for the design and application of multi-atom cooperative catalytic systems in multi-electron transfer CO2 reduction reactions.
![]() | ||
| Fig. 14 (a) Coordination environment modulation strategies for CO2 reduction and theoretical limiting potentials of CO2-to-CH4 conversion on Fe, Co, and Ni doped Cu55 clusters, along with the free energy diagram and optimized intermediate structures for CH4 formation on a Co666 cluster. Reproduced with permission from ref. 268, Copyright 2024 by Elsevier. (b) In situ ATR-SEIRAS spectra of Ce/CuOx-NPs and CuOx-NPs catalysts under varying potentials, proposed reaction pathways for electrochemical CO2 reduction to CH4, calculated Gibbs free energy diagrams, and projected density of states analyses for *CHO intermediate adsorption. Reproduced with permission from ref. 269, Copyright 2025 by Elsevier. (c) In situ ATR-SEIRAS spectra of Cu-CeO2 catalysts with varying compositions (Cu-CeO2-1.3, -3.0, -5.5), total peak area analysis of *CO and *CHO adsorbates, and the relative proportion of *CHO within their combined signal. Reproduced with permission from ref. 270, Copyright 2025 by American Chemical Society. | ||
Similarly, the atomic-level cerium-doped copper oxide catalyst (Ce/CuOx-NPs) reported by Chen's team269 achieved efficient conversion by regulating the electronic structure and the adsorption behavior of intermediates, as shown in Fig. 14b. It followed the reaction path of CO2 → *CO → *CHO → *CH2O → *CH3O → CH4. In situ ATR-SEIRAS characterization detected key intermediate signals of *CHO and *CH3O at 1490 cm−1 and 1397 cm−1, without the appearance of C2 product signals such as *OCCOH (1190 cm−1). This demonstrated high selectivity. DFT calculations indicated that cerium doping reduced the energy barrier of the rapid step *CO → *CHO to 0.33 eV, which was much lower than that of pure Cu (0.82 eV) and CuOx (1.10 eV), and simultaneously promoted the breakage of the C–O bond in *CH3O. Projection density of states (PDOS) analysis showed that the overlap degree of the d orbitals of Cu in Ce/CuO3 with the p orbitals of C in *CHO was significantly higher than that of pure Cu and CuO3, strengthening the electron interaction. In electrochemical tests, the catalyst achieved a CH4 faradaic efficiency (FECH4) of 67.4% at −1.6 V (vs. RHE), with a partial current density of 293 mA cm−2. The FE_CH4 remained above 62% within the voltage window of −1.4 to −1.75 V, and the 10-hour stability test showed no significant degradation in structure and performance, providing solid experimental and theoretical support.
In addition to atomic-level doping, the regulation of nanostructures has also achieved remarkable results. Xiong et al.270 developed mesoporous Cu-CeO2 nanospherical catalysts, which achieved nano-confined effect regulation through pore size engineering, providing a new path for reaction optimization. As shown in Fig. 14c, in situ ATR-SEIRAS characterization revealed that the intensity of the *CO and *CHO characteristic peaks at 2190–2110 cm−1 and 1480–1430 cm−1 increased with the reduction of pore size, confirming that small-sized pores can enhance the activation of CO2 and the accumulation of intermediates. Among them, Cu-mCeO2-3.0 with a pore size of 3.0 nm performed the best, and the peak area ratio of *CHO to *CO was the highest, indicating that the rate of *CO hydrogenation to *CHO (the rapid step of CH4 formation) was the fastest. This catalyst achieved a CH4 selectivity of 66.1 ± 2.9% at −1.4 V (vs. RHE), with a partial current density of 237.6 ± 14.5 mA cm−2, and maintained high selectivity within a wide potential window. Comparing Cu-mCeO2-1.3 with a pore size of 1.3 nm (where *CO and *CHO coverage was too high, resulting in C–C coupling) and Cu-mCeO2-5.5 with a pore size of 5.5 nm (with a weak confined effect), it was confirmed that CH4 selectivity has a volcano-like relationship with pore size.
The above study optimized the catalytic performance through two strategies, atomic-level doping (doping of transition metals in Cu clusters and Ce in CuOx) and nano-confined effect (control of mesopore diameters). It focused on three core dimensions: electronic structure, intermediate adsorption, and reaction pathway. All achieved efficient CO2 reduction to CH4. The common mechanism lies in precisely regulating the local environment and macroscopic structure of the catalyst, reducing the energy barriers of key steps, strengthening the conversion of target intermediates, and inhibiting side reactions. This provides a comprehensive reference for the multi-dimensional design of highly active and highly selective CO2-to-CH4 catalysts.
In summary, the core performance and characteristics of different types of 2D loadable catalysts in the CO2RR process for the production of C1 products can be summarized as shown in Table 1.
| Target product | Catalyst category | Representative catalyst | Key advantages | Critical performance metrics |
|---|---|---|---|---|
| CO | SACs | FeN5@graphene, F–Ni–N4, Fe–N3P–C | 100% atom utilization, high selectivity, tunable structure | FECO ≤ 98%, η ≤ 0.34V, h = 24 h |
| DACs | FeNi–NSC, PdPt@Cu, Ni–Cu@N–C | Dual-site synergy, optimized *CO desorption, suppressed HER | FECO ≤ 97.7%, TOF 20 695 h−1 |
|
| TACs | Zn–Ni1–Ni2@C, MnMoFe@graphyne | Multi-atom synergy, breaks scaling relationship | UL = 0.34 V, FECO >90% | |
| Clusters | Fe@meso-C, Fe–In–NC, AuCu alloycluster | Quantum-size effect, rich active sites, hybridized e-states | FECO ≤ 95%, enhanced *CO adsorption, durable | |
| HCOOH | SACs | MoSSe–Cu, Zr1@C2N, NiN2S2 | Pathway-specific, low η, robust structure | FEHCOOH ≤ 96%, E = 0.23 V |
| DACs | NiSb@C, Cr–NiFe@N–C | Charge accumulation, suppresses competing reactions | η ≤ 0.08 V, FEHCOOH >90% | |
| TACs | 3Ti-C2N1, MnMoCu@graphyne | Asymmetric active sites, favorable thermodynamics | UL = 0.42 V, strong HER suppression | |
| Clusters | Bi–Sn aerogel, In–Cu NP, Cu8 cluster | Prominent interface effect, fast mass transfer, dense sites | FEHCOOH ≤ 98.15%, j = 16.5 mA cm−2 | |
| CH4 | SACs | Ru1/CeO2, SA-Zn/MNC, RuB2N2@TiN | Single active center, precise intermediate tuning | FECH4 ≤ 85%, h = 35 h |
| DACs | Ru2-NC, Fe2Ir@NG | Bimetallic synergy, lowers multi-e transfer barrier | UL = 0.28 V, FECH4 >82% | |
| TACs | Ti3@NG, N,OH-Cu3 cluster | Stable framework, moderate *CO, suppressed C–C coupling | FECH4 ≤ 74.2%, j = 300 mA cm-2 | |
| Clusters | Co666@Cu55, Ce/CuOx NPs, Cu-mCeO2 | Tunable e-structure, nano-confinement, abundant sites | FECH4 ≤ 67.4%, j = 293 mA cm−2 |
Advanced associated probe characterization requires the development and application of associated probe spectroscopy and microscopy techniques, which can simultaneously provide chemical (e.g., XAS, Raman), structural (e.g., transmission electron microscopy at the same location, electrochemical scanning tunneling microscopy) and electronic information for the same catalyst under the same operating conditions. This will make it possible to directly correlate dynamic structural changes with catalytic activity and selectivity.
Despite progress, challenges remain in elucidating how electronic structure relates to electrochemical performance, critical for enhancing CO2RR efficiency and advancing practical applications. A primary hurdle is probing the dynamic surface states exhibited by catalysts and intermediates. For instance, while CO2RR studies often emphasize transitions between intermediates when analyzing energy barrier heights and rate-determining steps, conformational changes of intermediates—critical to overall reaction kinetics—are frequently overlooked. A single intermediate, such as *COOH, can adopt multi-site adsorption configurations involving different atoms, profoundly influencing subsequent pathways and final products.279–281
Moreover, intermediate transformations are not mere atomic additions, subtractions, or substitutions but complex dynamic processes involving bond formation/cleavage and interactions with both the catalyst and pre-adsorbed species. Concurrently, catalyst surfaces evolve significantly over reaction time—e.g., alloying in metal heterostructures—raising questions about the true active species and deactivation mechanisms under operational conditions. Understanding the evolution of actual active centers throughout long-term reaction processes is thus essential to ensuring catalyst stability and reliability. Current research largely confines catalyst-adsorbate interactions to single adsorbed molecules, neglecting cumulative effects of multiple intermediates or complex reaction patterns arising from interactions between numerous intermediates and the catalyst—phenomena common in practical scenarios. During prolonged CO2RR, such clustering becomes more pronounced, with diverse intermediates coexisting on the surface. This oversight also extends to the transformation mode of intermediates. Because when the surface is densely covered, its transformation behavior may be completely different from that of isolated molecules, which highlights the inherent limitations of guiding experiments by simulating single-molecule dynamics.282–284
Expanding research to encompass dynamic interactions between catalysts and multiple intermediates is therefore urgent to comprehensively understand the CO2RR under practical conditions. DFT calculations have successfully guided the screening of atomic catalysts,278,285–287 but computational cost and time restrict the exploration of coordination environments around central metal atoms, rendering traditional DFT unsuitable for large-scale systems. With advancing computing power, first-principles high-throughput screening is poised to guide the design of superior catalysts, when combined with in situ characterization, such methods can also screen reaction intermediate structures, further advancing mechanistic studies of dynamics. However, these calculations typically assume ideal environments, whereas electrocatalytic processes involve highly complex external conditions and internal environments, leading to discrepancies between experimental results and theoretical simulations. Larger models are needed to fully simulate electrocatalytic environments, demanding the timely development of new methodologies.
Artificial intelligence (AI) offers innovative approaches to catalyst design, integrating advanced computing technologies e.g., DFT and machine learning (ML) assisted models.285,288–292 Operating under a variety of conditions can make predictions more directly relevant to real-world use—traditional high-throughput screening typically operates only under “standard” conditions, but these do not account for temperature, pressure, or pH. Building and supporting data sets and open databases are key to improving the effectiveness of the training sets for AI. Catalyst activity originates from active atoms, support materials, coordination environments, and defect distributions, and establishing multi-factor structure–activity relationships models and multi-scale simulations enables effective integration of structural information to refine comprehension of structure–activity relationships. AI-enabled design is capable of enhancing model precision and the efficiency of experimental verification, which in turn speeds up the conversion of theoretical research into real-world applications.
The future of the CO2RR lies in the integration of precise experiments with complex theories. By resolving the intertwined experimental and computational challenges mentioned above, this field can progress from accidental discoveries to the rational design of next-generation catalysts. The ultimate goal is to develop reliable, highly selective, and energy-efficient catalytic systems that can be integrated into practical electrolyzers, converting carbon dioxide from waste into sustainable fuel and chemical sources, and making meaningful contributions to the circular carbon economy.
Among these catalyst systems, metal clusters stand out as a critical architecture for advancing the CO2RR, especially for complex products requiring multi-step electron transfer and C–C coupling. Composed of several to hundreds of atoms (1–3 nm in size), clusters integrate the advantages of atomic-level precision and nano-scale functionality, exhibiting unique quantum size effects and hybrid electronic structures that bridge molecular-like states and bulk material band structures. Their structural diversity including constituent atom count, elemental proportions, atomic arrangement, and geometric symmetry enables the creation of abundant active sites tailored to adsorb and convert specific reactants and intermediates in the CO2RR. Unlike isolated SACs/DACs/TACs, clusters provide a synergistic multi-nuclear environment that effectively modulates the adsorption energetics of key intermediates (e.g., *COOH, *OCHO, *CO) and reduces kinetic barriers for rate-determining steps, particularly C–C coupling for C2+ products and deep hydrogenation for CH4.
In the reduction of CO2 to C1 products, clusters demonstrate remarkable performance advantages. For CO production, Fe-based clusters anchored on mesoporous carbon, Fe–In–NC clusters, and AuCu alloy clusters leverage quantum-size effects and rich active sites to enhance *CO adsorption and suppress the HER, achieving CO faradaic efficiencies (FE) up to 95% and durable stability. In HCOOH synthesis, Bi–Sn bimetallic aerogels, In–Cu nanoparticles, and double tetrahedral Cu8 clusters exhibit prominent interface effects and fast mass transfer, with FEHCOOH reaching as high as 98.15% and current densities of 16.5 mA cm−2. For CH4 generation, Co666@Cu55 clusters, Ce/CuOx NPs, and mesoporous Cu–CeO2 nanospheres rely on tunable electronic structures and nano-confinement effects to optimize intermediate adsorption and hydrogenation pathways, achieving FECH4 up to 74.2% at industrial-grade current densities (300 mA cm−2). Notably, high-entropy alloy clusters (e.g., CoNiCuRh@FeN4) have emerged as a frontier, with their ultra-high configurational entropy endowing unique electronic structures and structural stability, enabling precise regulation of multiple product pathways.
This review highlights recent progress in classical two-dimensional material catalysts for the CO2RR and strategies to regulate the selectivity of C1+ products. By systematically examining the performance of diverse electrocatalysts, we elucidate reaction mechanisms and design principles at a fundamental level. The pursuit of efficient CO2RR catalysts reveals a fundamental design principle, in which the number and size of active sites determine the reaction pathway and product selectivity. From SACs to DACs and TACs, and finally to the development of atomic clusters, a clear scale-functional paradigm has been defined. SACs feature uniform and isolated metal centers, which exhibit excellent selectivity for simple two-electron transfer products such as CO while maximizing atomic utilization efficiency. However, their peculiar geometric structures inherently restrict the adsorption configuration of key intermediates, thereby hindering complex multi-step reactions, most crucially, C–C coupling.
To overcome this limitation, DACs and TACs introduced adjacent metal sites, creating a unique coordinated environment and thus achieving a synergy effect. This atomic-level partnership allows for the simultaneous stabilization of multiple intermediates or the cleavage of specific bonds, effectively circumventing the limitations imposed by linear scaling relationships in traditional heterogeneous catalysis. This development has opened up thermodynamic and kinetic pathways for more valuable multi-electron products such as HCOOH, C2H4 and CH4. When the catalytic target further shifts to high-value C2+ products, atomic clusters typically composed of several to hundreds of atoms become the preferred architecture. Their multi-nuclear structure provides a set of active sites with hybrid electronic states, which are particularly good at reducing the kinetic barrier of C-C coupling, a step that is often regarded as the bottleneck for the formation of C2+. Therefore, the complexity of the required CO2RR product combination directly determines the optimal scale of the active site, among which SACs exhibit excellent catalytic performance for the generation of CO. The synthesis of C2+ substance requires the synergy of multi-atom sets, although the latter faces greater challenges in terms of precise synthesis and structural stability.
2D Electron Gas Properties and a Route toward Graphene-based Nanoelectronics, J. Phys. Chem. B, 2004, 108, 19912–19916 Search PubMed.| This journal is © The Royal Society of Chemistry 2026 |