With atmospheric CO
2 concentrations in the atmosphere reaching unacceptably high levels and continuing to rise, capturing atmospheric CO
2 and converting it into value-added fuels has become both essential and sustainable [
1–
3]. The electrochemical CO
2 reduction reaction (CO
2RR) stands out as an effective and promising strategy for transforming CO
2 into high-value chemicals and liquid fuels [
4]. Among the various products of electrochemical CO
2RR, CO serves as a fundamental feedstock for methanol synthesis, Fischer-Tropsch synthetic oils, and diverse carbonylation reactions [
5]. However, this conversion process involves two proton-coupled electron transfer steps, posing challenges such as identifying active sites and competing with the hydrogen evolution reaction (HER) [
6–
8]. In recent years, single-atom catalysts (SACs) have attracted significant attention due to their well-defined active sites and excellent selectivity for C
1 products [
6,
9–
11].
Carbon-based SACs play a pivotal role in electrocatalysis due to their versatile structures, exceptional stability, and high conductivity [
12–
14]. Their high activity in CO
2RR systems is typically associated with symmetric M-N
4 sites, such as Cu-N
4 [
15], Ni-N
4 [
16], or Co-N
4 [
17]; yet, challenges regarding activity, stability, and product selectivity persist. For instance, excessively strong
*CO adsorption energy or a high
*COOH energy barrier can result in unsatisfactory selectivity and stability [
18,
19]. To address these issues, dual-atom catalysts (DACs) with diverse metal combinations have been developed, leveraging interactions between the two metals to effectively lower reaction barriers, accelerate reaction rates, and thus achieve superior catalytic performance [
20–
22]. However, to the best of our knowledge, recent research has predominantly focused on comparing DACs with SACs, with limited studies investigating how intermetallic interactions influence CO
2RR performance in DACs systems [
23,
24].
This study explores the impact of intermetallic interactions on the CO2RR performance of Cu/Ni-NC DACs. Theoretical simulations were conducted on Cu/Ni-NC DACs models to assess how intermetallic interactions affect electronic structure, CO2RR activity, and selectivity, ultimately identifying the reason for the selectivity shift in Cu-N-C catalysts after Ni doping. These findings highlight that intermetallic interactions are a critical and reliable parameter influencing catalytic activity in DACs systems. Aberration-corrected high-angle annular dark-field scanning transmission electron microscopy (AC—HAADF-STEM) and X-ray absorption fine structure (XAFS) analyses further confirmed that the synthesized Cu/Ni-N-C catalyst features an atomic configuration with neighboring Cu/Ni sites coordinated with four nitrogen atoms, respectively. The resulting catalyst efficiently facilitates the selective conversion of CO2 to CO, achieving nearly 100% CO Faradaic efficiency (FECO) at −1.3 V vs. RHE in an H-type cell. Density functional theory (DFT) calculations reveal that the Cu/Ni dual-atom catalyst synergistically electroreduces CO2 to CO: The Cu atom, with low charge density, adsorbs the oxygen atom of the CO2 molecule; while the Ni atom, with high charge density, adsorbs the carbon atom of the CO2 molecule, collectively enhancing overall CO2RR activity. This study underscores the importance of tuning intermetallic interactions in DACs to maximize synergistic effects and optimize CO2RR performance.
The fabrication of Cu-Ni dual-atom catalysts is schematically illustrated in
Fig. 1. 1,3,5-Benzenetricarboxylic acid, copper(Ⅱ) nitrate trihydrate, and PVP were used as starting materials to synthesize HKUST-1 metal-organic framework (MOF)
via a hydrothermal reaction (see Experimental Section in Supporting information). When copper(Ⅱ) nitrate trihydrate was partially substituted with nickel nitrate hexahydrate, with all other experimental procedures kept unchanged, (CuNi)HKUST-1 MOFs with different Cu/Ni ratios were successfully obtained. The crystal phase and morphology of the as-synthesized samples were characterized by X-ray diffraction (XRD, Fig. S1 in Supporting information) and scanning electron microscopy (SEM,
Figs. 2a-c and Fig. S2 in Supporting information). Results confirmed the successful synthesis of HKUST-1 MOF and showed no changes in the structure or morphology of HKUST-1 after partial substitution of Cu atoms with Ni atoms. Finally, HKUST-1 and (Cu
xNi
y)HKUST-1 MOFs were mixed with dicyandiamide and then calcined; the resulting products were designated as Cu
xNi
y-N-C-BL (where x:y = 1:0, 0.8:0.2, 0.5:0.5 and 0:1; BL means before acid leaching). After acid leaching under an oxygen atmosphere, metal particles in the Cu
xNi
y-N-C-BL samples were removed, yielding target products of Cu
xNi
y-N-C. For comparison, nitrogen-doped carbon (N—C) was also synthesized without adding any metal components.
Scanning electron microscopy (SEM) images show that all Cu
xNi
y-N-C-BL samples consist of interconnected foam structures with wrinkled, entangled flakes (Fig. S7 in Supporting information), which is similar to the Cu-N-C sample (
Fig. 2e). The low-resolution transmission electron microscopy (TEM) image of the Cu-N-C catalyst reveals cross-linked, ultrathin folded-sheet features (
Fig. 2f) and the corresponding high-resolution TEM (HRTEM) image confirms no observable copper particles in the Cu-N-C catalyst (
Fig. 2g). The selected-area electron diffraction (SAED) pattern (inset of
Fig. 2g) further confirms the sample’s non-crystalline nature. Energy-dispersive spectroscopy (EDS) elemental mapping shows the existence of Cu, N, and C elements and the uniform distribution of N and Cu across the entire carbon matrix (
Fig. 2h). As presented in
Figs. 2i-k and Figs. S8-S11 (Supporting information), the Cu
0.5Ni
0.5—N-C, Cu
0.8Ni
0.2—N-C and Ni-N-C samples display a morphology analogous to that of the Cu-N-C sample, with no distinct metal particles observed on their surfaces. EDS elemental mapping further verifies the homogeneous distribution of N, Cu, and Ni across the carbon support (
Figs. 2l-m and Figs. S9-S11). As illustrated in the AC—HAADF-STEM image (
Fig. 2n), numerous bright spot pairs are distinctly observed, which are indicative of potential atomic-scale active sites. Furthermore, the 3D surface plot derived from region 1 of
Fig. 2n offers compelling evidence to corroborate the existence and structural features of these atomic pairs (
Fig. 2o). In addition, statistical analysis presented in
Fig. 2p clearly demonstrates that over 88.73% of the detected sites correspond to isolated dual-atom configurations (circled with red dashed lines) with a distance around 2.81–2.92 Å (
Fig. 2q). The electron energy-loss spectroscopy (EELS) diagram in Fig. S12 (Supporting information) further confirmed the co-existence of Ni and Cu atomic pairs. The small interatomic distance between Cu and Ni induces intense through-space and through-bond electronic coupling interactions, effectively expediting the migration of electrons between Cu and Ni. Collectively, these characterizations confirm that Cu/Ni dual-atom sites are not only the dominant active species but also uniformly distributed throughout the Cu
0.5Ni
0.5—N-C samples, laying a solid foundation for their enhanced catalytic performance. The Raman spectra of Cu-N-C, Cu
0.8Ni
0.2—N-C, and Cu
0.5Ni
0.5—N-C samples exhibit two prominent characteristic peaks at approximately 1355 and 1570 cm
-1 (Fig. S13 in Supporting information), which are ascribed to the d-band and G-band, respectively [
25,
26]. Defect concentrations were quantified by normalizing the intensity of the d-band relative to that of the G-band (
i.e., the
ID/
IG ratio, Table S1 in Supporting information) [
25]. Notably, the
ID/
IG values were calculated as 1.09, 1.08 and 1.07 for Cu-N-C, Cu
0.5Ni
0.5—N-C and Cu
0.8Ni
0.2—N-C samples, respectively, which are comparable and indicate that the defect content in these Cu
xNi
y-N-C samples is not the critical factor governing their catalytic performance. N
2 adsorption-desorption isotherms and corresponding pore size distribution analyses (Fig. S14 in Supporting information) confirm that all samples share comparable Brunauer-Emmett-Teller (BET) specific surface areas and total pore volumes. For example, the Cu
0.5Ni
0.5—N-C sample achieves a BET specific surface area of 137.71 m
2/g and a total pore volume of 0.62 cm
3/g, which is similar to that of the Cu-N-C and Ni-N-C sample. Such favorable textural characteristics are capable of optimizing the mass transfer behavior of CO
2 and electrolyte at the triple-phase interface, thereby creating favorable conditions for the smooth progression of the CO
2 electroreduction reaction.
X-ray photoelectron spectroscopy (XPS) measurements were conducted to systematically investigate the chemical composition and elemental valence states of the Cu
xNi
y-N-C and Cu-N-C samples. As revealed by the XPS survey spectra (Fig. S15 in Supporting information), Cu
xNi
y-N-C contained four main elements,
i.e., C, N, Cu and Ni, confirming the successful incorporation of N, Cu and Ni into the carbon sheet matrix. Notably, the N dopant is inherently derived from the dicyandiamide precursor. As shown in Fig. S16 (Supporting information), the high-resolution C 1s spectrum of the Cu
0.5Ni
0.5—N-C sample deconvolutes into four major peaks centered at 284.6, 286, 287.3 and 289.1 eV, which correspond to sp
2-bonded C=C, sp
3-bonded C—N, C—O, and C=O functional groups, respectively [
27].
Fig. 3a presents the high-resolution N 1s XPS spectrum of Cu
0.5Ni
0.5—N-C, revealing five distinct sub-peaks at 398.7, 399.2, 400.1, 401.1 and 403.1 eV. These peaks are assigned to pyridinic N, Cu/Ni-N
x, pyrrolic N, graphitic N, and oxidized N, respectively. The presence of the Cu/Ni-N
x peak in the deconvoluted N 1s spectrum provides direct evidence for the formation of Cu/Ni-N
x coordination groups in the sample [
27–
29]. For the Cu 2p region (
Fig. 3b), the peaks at 932.50 and 934.8 eV are attributed to the 2p
3/2 orbitals of Cu(Ⅰ) and Cu(Ⅱ) species, while those at 952.2 and 954.7 eV correspond to the 2p
1/2 orbitals of Cu(Ⅰ) and Cu(Ⅱ), respectively. Additional satellite peaks at 942.8 and 944.9 eV further confirm the presence of Cu oxide phases in the carbon matrix, which is further confirmed by the Cu Auger spectra in Fig. S17 (Supporting information) [
27,
30]. In the high-resolution Ni 2p spectrum of Cu
0.5Ni
0.5—N-C (
Fig. 3c), the peaks at 855.29 and 872.78 eV are assigned to Ni 2p
3/2 and Ni 2p
1/2, respectively. The absence of the metallic Ni 2p
3/2 peak (at 852.6 eV) indicates that the dispersed Ni species observed in HAADF-STEM images exist as Ni
+ or Ni
2+ cations [
31]. For comparison, XPS measurements were conducted on Cu-N-C, Ni-N-C and Cu
0.8Ni
0.2—N-C under the same conditions (Figs. S18-S20 in Supporting information). The results indicated that the chemical states of the major elements are analogous to those of Cu
0.5Ni
0.5—N-C. Notably, the divergent orbital hybridization behaviors and electron-donating properties of Cu and Ni in their N-coordination interactions trigger an asymmetric charge redistribution in the Cu-Ni diatomic pair. As observed, the binding energy of Cu 2p
3/2 in Cu-N-C is 932.35 eV, and that of Ni 2p
3/2 in Ni-N-C is 855.47 eV. While the Cu 2p
3/2 binding energy in Cu
0.5Ni
0.5—N-C shifts to 932.50 eV (an increase of 0.15 eV), and the Ni 2p
3/2 binding energy shifts to 855.29 eV (a decrease of 0.18 eV) (Fig. S21 in Supporting information). These shifts directly confirm the proposed electron transfer from Cu to Ni in the dual-atom catalyst.
X-ray absorption near-edge spectroscopy (XANES) and Fourier-transformed extended X-ray absorption fine structure (FT-EXAFS) were employed to characterize the structural and electronic properties of Cu and Ni sites in the Cu
0.5Ni
0.5—N-C sample. As shown in
Fig. 3d, the Cu K-edge XANES spectrum of the Cu
0.5Ni
0.5—N-C catalyst places the Cu absorption peak between those of CuPc and Cu foil, confirming the oxidized state of Cu. Similarly,
Fig. 3e displays the Ni K-edge XANES spectrum of the Cu
0.5Ni
0.5—N-C sample, where the Ni absorption peak is positioned between those of NiO/NiPc and a Ni foil, indicating that Ni is present in an oxidized state. For the FT-EXAFS analysis,
Fig. 3f displays the Cu K-edge FT-EXAFS spectrum of the Cu
0.5Ni
0.5—N-C sample, which is characterized by a dominant peak at 1.54 Å and can be assigned to Cu-N coordination environments. Notably, a small peak positioned at 2.12 Å was observed, which is due to the weak interaction between Cu and adjacent Ni atoms. Correspondingly, the Ni K-edge FT-EXAFS spectrum (
Fig. 3g) shows a major peak at 1.56 Å for the sample, associated with Ni-N bonding. To determine the detailed structural configuration of Cu and Ni in Cu
0.5Ni
0.5—N-C sample, the FT-EXAFS spectra were further fitted (
Figs. 3h and i, Figs. S22 and S23 in Supporting information). The fitting results (Tables S2 and S3 in Supporting information) reveal that the coordination numbers (CN) of Cu-N and Ni-N in the catalyst are 4.03 and 4, respectively. Wavelet transforms (WT) were used to validate the above conclusions, as they provide resolution in both R and k space. The WT contour plots of the Cu
0.5Ni
0.5—N-C catalyst exhibit a strong signal at 2.68 Å
-1, which corresponds to Cu-N contributions and is consistent with the signal of CuPc (
Fig. 3j and Fig. S24 in Supporting information). A comparable finding was obtained from the WT-EXAFS analysis of the Ni K-edge (
Fig. 3k and Fig. S25 in Supporting information).
To evaluate the performance of Cu
xNi
y-N-C catalysts, each catalyst powder was drop-coated onto a carbon paper with a Nafion binder to fabricate the working electrode. CO
2 reduction reaction tests were conducted in a three-electrode H-cell configuration using 0.1 mol/L KHCO
3 aqueous electrolyte saturated with CO
2. Linear sweep voltammetry (LSV) curves were recorded for all samples in both Ar- and CO
2-saturated 0.1 mol/L KHCO
3 solutions. As shown in
Fig. 4a, the current densities were notably higher under CO
2 atmosphere compared to Ar, indicating their strong catalytic activities toward CO
2RR. The potential-dependent gas product distribution of all samples was measured between −1 V and −1.5 V
vs. RHE, with in-line gas chromatography (GC) used for gas product quantification. For the Cu-N-C sample, CO
2RR yielded main products including CO and CH
4 with maximum Faradaic efficiency (FE) of CO 38.35% at −1 V
vs. RHE and maximum FE
CH4 21.89% at −1.4 V
vs. RHE (
Fig. 4b). In contrast, Cu
0.8Ni
0.2—N-C and Cu
0.5Ni
0.5—N-C exhibited distinct catalytic selectivity. The Cu
0.8Ni
0.2—N-C catalyst achieved a maximum FE
CO of 66.69% at −1.3 V
vs. RHE, corresponding to a partial current density of 10.06 mA/cm
2 (Fig. S26 in Supporting information).
Interestingly, increasing the Ni content enhanced CO generation, making it the dominant product. For Cu
0.5Ni
0.5—N-C, the maximum FE
CO reached 96.78% with a partial current density of 24.80 mA/cm
2 at −1.3 V
vs. RHE and negligible CH
4 and liquid products were detected (
Fig. 4c, Figs. S27 and S28 in Supporting information), exhibiting a different product selectivity from that of Cu-N-C catalyst. Notably, the optimal Cu
0.5Ni
0.5—N-C dual-atom catalyst displays a superior ability for CO formation relative to the Ni-N-C single-atom catalyst, where the latter yields a maximum FE
CO of 86% at −1.3 V
vs. RHE with a partial current density of 20.59 mA/cm
2 (Fig. S29 in Supporting information). The distinct shift in selectivity arises from the formation of Cu-Ni atomic pairs after Ni doping into the Cu-N-C catalyst. As Ni has a higher electronegativity than copper, the electron cloud density of Cu atoms is reduced, which weakens the adsorption of intermediate species during the CO
2 reduction reaction. This effect hinders the subsequent protonation process of the CO
2RR intermediates (
e.g.,
*CO), thereby suppressing CH
4 formation and elevating CO selectivity. The reduced electron cloud density around Cu atoms in the Cu/Ni-N-C catalyst, in comparison with the Cu-N-C catalyst, is further verified by the calculations of charge density difference and electron localization function (
Figs. 5a-c, Figs. S36 and S37 in Supporting information). For comparison, the N—C catalyst was tested and showed negligible CO
2RR activity. It exhibited nearly full Faradaic efficiency toward HER across various applied potentials (
Fig. 4d). These results confirm that Ni and Cu species serve as the main active sites for CO
2RR in Cu
xNi
y-N-C catalysts. To gain deeper insight into reaction kinetics, Tafel analysis was performed as presented in
Fig. 4e and the potential window adopted for Tafel fitting ranges from −0.75 V to −1.3 V
vs. RHE. This range was chosen because FE
CO stays above 90% to guarantee CO
2RR dominance, current density increases linearly with overpotential (a signature of kinetic control free of mass transfer constraints), and negligible side reactions like hydrogen evolution take place. The Cu
0.5Ni
0.5—N-C catalyst had a Tafel slope of 67.43 mV/dec, lower than those of Ni-N-C (77.57 mV/dec), Cu-N-C (122.41 mV/dec), and Cu
0.8Ni
0.2—N-C (97.60 mV/dec). This indicates that Cu
0.5Ni
0.5—N-C effectively promotes the first electron transfer to form the adsorbed
*COOH intermediate, accelerating overall CO
2RR dynamics. Electrochemical impedance spectroscopy (EIS) measurements further validated the reaction kinetics. As shown in
Fig. 4f, the
Rct values of N—C, Cu-N-C, Ni-N-C, Cu
0.8Ni
0.2—N-C, and Cu
0.5Ni
0.5—N-C are approximately 3055, 368.5, 265.7, 279.5 and 167.7 Ω, respectively, and the Cu
0.5Ni
0.5—N-C displays the smallest semicircle, which means faster electron transfer kinetics compared to the other catalysts, demonstrating efficient interfacial charge transfer.
The electrochemically active surface area (ECSA), a key factor influencing intrinsic catalytic activity, is proportional to the double-layer capacitance (
Cdl). Cyclic voltammetry (CV) measurements were used to determine
Cdl values for N—C, Cu-N-C, Ni-N-C, Cu
0.8Ni
0.2—N-C, and Cu
0.5Ni
0.5—N-C (Fig. S30). The calculated
Cdl values were 2.11 mF/cm
2 (N—C), 6.57 mF/cm
2 (Cu-N-C), 12.07 mF/cm
2 (Ni-N-C), 9.94 mF/cm
2 (Cu
0.8Ni
0.2—N-C), and 13.63 mF/cm
2 (Cu
0.5Ni
0.5—N-C) (
Fig. 4g). These results further support the significantly enhanced CO
2RR activity of Cu
0.5Ni
0.5—N-C. Furthermore, the superior intrinsic activity of Cu
0.5Ni
0.5—N-C is further highlighted by calculating the turnover frequency (TOF). The developed Cu
0.5Ni
0.5—N-C exhibits an exceptionally high TOF of 15,372 h
-1 at −1.3 V
vs. RHE, outperforming most of the reported state-of-the-art single atom catalysts evaluated under similar conditions (Table S4 in Supporting information). Cu
0.5Ni
0.5—N-C catalyst shows excellent durability with negligible decrease in current density, and FE
CO approximately remains about 90% after continuous electrocatalysis for 10 h (
Fig. 4h). Post-stability characterizations
via XRD, SEM, AC—HAADF-STEM image and XPS (
Figs. 4i-k and Figs. S31-S35 in Supporting information) revealed no significant changes in the catalyst’s phase structure, morphology, or elemental composition, confirming its good stability. This outstanding durability is due to the strong coupling interaction of N-coordinated dual atomic sites in the carbon support. This interaction reduces
*CO poisoning and accelerates CO desorption from active sites, thereby enhancing stability during CO
2RR.
To understand the origin of the enhancement in the catalytic activity of Cu/Ni-N-C catalyst, DFT calculations were conducted to analyze its electronic properties. Charge density difference (CDD) analysis reveals pronounced electron accumulation (yellow region) on Ni atoms and electron depletion (cyan region) on Cu atoms. This directly visualizes the electron density bias toward Ni and thus indicates a stronger interaction between Ni atoms and the intermediates of the electrochemical CO
2 reduction reaction (
Fig. 5a). More electron accumulation around Ni active centers than around Cu atoms is also observed in both Cu-N-C and Ni-N-C configurations (
Fig. 5b and Fig. S36 in Supporting information), further verifying that the catalytic activity of Ni-N sites is higher than that of Cu-N sites. This observation is further corroborated by electron localization function (ELF) analysis, which reveals localized electron density between the Cu-N and Ni-N moieties and demonstrates that Ni-N sites exhibit higher electron density than Cu-N sites (
Fig. 5c and Fig. S37 in Supporting information). This electron distribution is likely to have a substantial influence on the adsorption behavior of the intermediates and the reaction selectivity. When CO
2 is polarized, its lowest unoccupied molecular orbital is concentrated mainly on the carbon atom, making it vulnerable to nucleophilic attacks. Nevertheless, the presence of regions with localized positive charge on the catalyst surface can impede nucleophilic attacks on the carbon atom. This, in turn, raises the thermodynamic barrier for
*COOH adsorption. As a result, for the Cu/Ni-N-C catalyst, the carbon atom of polarized CO
2 molecules preferentially adsorbs onto the electron-enriched Ni active sites, while the oxygen atoms of CO
2 tend to bind more strongly to the electron-deficient Cu active sites (
Fig. 5d). This unique adsorption mode is a direct consequence of the electron density bias toward Ni and serves as theoretical evidence for the charge transfer from Cu to Ni, which is consistent with the XPS results.
The shift from monometallic Ni-N-C and Cu-N-C to bimetallic Cu/Ni-N-C brings about remarkable alterations to the electronic structure of the catalysts. To establish how such electronic structural changes are associated with CO
2RR catalytic activity, orbital analysis was further performed in this work. As Cu/Ni-N-C delivered boosted catalytic performance with metal sites acting as the active sites, this specific configuration was prioritized for systematic investigation and compared against Cu-N-C, aiming to unravel the underlying reasons for the disparities in their catalytic behaviors. Nickel possesses a valence electron configuration of 3d
84s
2, whereas copper has that of 3d
104s
1. Within the Cu-Ni dual-atom architecture, the unoccupied 3d orbitals of Ni (3d
8) are able to accept electron density from the fully filled 3d orbitals of Cu (3d
10), thereby enabling orbital hybridization-driven interatomic electron transfer. PDOS calculations demonstrate that the HOMO-LUMO band gap is narrowed in Cu/Ni-N-C, a feature that points to improved electronic conductivity and consequently higher catalytic activity potential, which stems from more efficient electron transfer throughout the reaction process (
Fig. 5e and Fig. S38 in Supporting information). In addition, the introduction of Ni into the Cu-N-C matrix causes the d-band center of Cu 3d to positively shift from −3.396 eV to −2.844 eV (
Fig. 5e), which infers stronger binding interactions with reaction intermediates. An analogous positive shift is also detected in the PDOS spectrum for the Ni 3d d-band center (Fig. S38).
The surface coverage of
*H species over catalysts with distinct hybrid heterostructures was quantified by evaluating the free energy variations related to water dissociation (Δ
GH2O) and hydrogen adsorption-desorption (Δ
GH) on Cu-N-C, Ni-N-C and Cu/Ni-N-C catalysts (
Figs. 5f and g). Specifically, Cu/Ni-N-C exhibits a Δ
GH2O value of 0.61 eV, which is considerably lower than those of Cu-N-C (1.65 eV) and Ni-N-C (1.20 eV) (
Fig. 5f). This result manifests that the incorporation of Ni endows Cu/Ni-N-C with more thermodynamically favorable water dissociation kinetics, thus yielding an elevated
*H concentration that is available for the subsequent electrocatalytic reactions. On the contrary, the absolute Δ
GH value of Cu/Ni-N-C is 1.08 eV, which is markedly lower than the corresponding values of Cu-N-C (1.96 eV) and Ni-N-C (1.86 eV) (
Fig. 5g). To further clarify the CO formation mechanism over Cu/Ni-N-C, DFT calculations were utilized to simulate the Gibbs free energy profiles of the complete reaction sequence. As illustrated in
Fig. 5h, the
*COOH formation on Cu/Ni-N-C is energetically more favorable (Δ
G = 1.78 eV) than on Cu-N-C (Δ
G = 2.52 eV) and Ni-N-C (Δ
G = 2.24 eV). This suggests that Cu/Ni-N-C was more conducive to CO formation than the Cu-N-C and Ni-N-C catalysts. Furthermore, the thermodynamic barrier for
*COOH formation is found to be higher on Cu-N-C (Δ
G = 2.52 eV) than on Ni-N-C (Δ
G = 2.24 eV), a result that verifies Ni (rather than Cu) as the preferred metal site for anchoring the C atom of
*COOH moieties in Cu/Ni-N-C, which aligns with the findings in
Fig. 5d.
To clarify the fundamental reason for this catalytic selectivity, we analyzed the structural configurations of the key
*COOH intermediate on the surfaces of Cu-N-C, Ni-N-C and Cu/Ni-N-C catalysts (
Fig. 5i). On Cu-N-C, the Cu-C bond length in the
*COOH intermediate is 2.027 Å, and its C—OH bond length is 1.376 Å. In the same manner, the Ni-N-C surface exhibits a Ni-C bond length of 2.020 Å and a C—OH bond length of 1.360 Å for the adsorbed
*COOH. In sharp contrast, the Cu/Ni-N-C surface shows a shortened Ni-C bond length of 1.921 Å for
*COOH, as well as a corresponding C—OH bond length of 1.388 Å. Such structural changes are of great significance to the underlying reaction mechanism of CO
2RR. A longer C—OH bond on Cu/Ni-N-C points to a weaker interaction between C and O atoms, thus making the C—OH bond easier to break. On the other hand, the shorter Ni-C bond on Cu/Ni-N-C means a stronger binding between the metal site and carbon atom, which is more favorable for the formation of the
*CO intermediate. This analysis of intermediate configurations offers a molecular-level interpretation for the improved
*CO formation efficiency observed on the Cu/Ni-N-C catalyst, and these structural understandings also confirm the validity of our earlier energy calculations.
This study systematically investigates the structure and CO2 electroreduction performance of Cu/Ni dual-atom catalysts (Cu/Ni-N-C) through a combination of experiments and theoretical calculations, clarifying the regulatory mechanism of intermetallic interactions on catalytic selectivity. Experimental results confirm that Cu and Ni atoms in Cu/Ni-N-C are each coordinated with four nitrogen atoms, forming uniformly distributed dual-atom active sites. The obtained Cu-N-C catalyst exhibits maximum FECO of 38.35% at −1 V vs. RHE and maximum FECH4 of 21.89% at −1.4 V vs. RHE, while the Cu0.5Ni0.5—N-C catalyst exhibits exceptional CO2RR performance: At −1.3 V vs. RHE, the CO Faradaic efficiency (FECO) reaches 96.78% with a partial current density of 24.80 mA/cm2. The observed selectivity variation is caused by the electron-withdrawing effect of Ni atoms on Cu atoms in the Cu/Ni-N-C catalyst, which lowers the electron cloud density around Cu atoms and thus weakens the adsorption of CO2RR intermediates on Cu sites. DFT calculations reveal the core mechanism underlying the enhanced catalytic performance for CO formation: Charge regulation between Cu and Ni forms an asymmetric adsorption mode. Electron-deficient Cu atoms act as O-terminal adsorption sites to bind the oxygen atoms of CO2; while electron-enriched Ni atoms serve as C-terminal adsorption sites to anchor the carbon atoms. This synergistic adsorption promotes CO2 polarization and activation, narrows the HOMO-LUMO band gap for efficient electron transfer, and reduces the thermodynamic barrier for *COOH intermediate formation. The Cu/Ni-N-C not only strengthens the binding between the metal site and carbon atom of CO2 but also weakens the C—OH interaction, facilitating CO generation. This work provides new insights into the rational design of dual-atom catalysts and highlights the potential of regulating intermetallic interactions to optimize CO2RR performance, advancing their industrial applications.