Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Coordination tailoring of Cu single sites on C3N4 realizes selective CO2 hydrogenation at low temperature

## Abstract

CO2 hydrogenation has attracted great attention, yet the quest for highly-efficient catalysts is driven by the current disadvantages of poor activity, low selectivity, and ambiguous structure-performance relationship. We demonstrate here that C3N4-supported Cu single atom catalysts with tailored coordination structures, namely, Cu–N4 and Cu–N3, can serve as highly selective and active catalysts for CO2 hydrogenation at low temperature. The modulation of the coordination structure of Cu single atom is readily realized by simply altering the treatment parameters. Further investigations reveal that Cu–N4 favors CO2 hydrogenation to form CH3OH via the formate pathway, while Cu–N3 tends to catalyze CO2 hydrogenation to produce CO via the reverse water-gas-shift (RWGS) pathway. Significantly, the CH3OH productivity and selectivity reach 4.2 mmol g–1 h–1 and 95.5%, respectively, for Cu–N4 single atom catalyst. We anticipate this work will promote the fundamental researches on the structure-performance relationship of catalysts.

## Introduction

CO2 hydrogenation has attracted great attention because it can simultaneously reduce the CO2 emission1,2,3,4 and produce high value-added chemicals5,6,7,8, yet it suffers from low activity and selectivity due to the chemical inertness of CO2 and the complicated process of CO2 hydrogenation9,10,11. To promote CO2 hydrogenation, one needs to raise the temperature or pressure, which however results in increased selectivity of CO and decreased selectivities of value-added chemicals. For instance, Cu-based catalysts have been widely used for CO2 hydrogenation, which are generally operated at temperature and pressure ranges of 220–300 °C and 50–100 atm, respectively10,11. Over the past decades, great efforts have been devoted to developing efficient catalysts for CO2 hydrogenation, such as compositing Cu sites with metal oxides (e.g., ZnO, ZrO2, and TiO2)12,13,14,15,16,17,18,19,20. However, the request for efficient CO2 hydrogenation catalysts is driven by the current drawbacks including poor activity and/or CH3OH selectivity, as well as the ambiguous mechanism.

Recently, single atom catalysts (SACs) have attracted great attention in catalysis due to their properties including 100% atom utilization efficiency, low coordination environment, confinement effect, etc21,22,23,24,25,26,27,28,29,30. Typical for CO2 hydrogenation, SACs are regarded as promising catalysts for producing high value-added chemicals. For instance, Ye et al. synthesized Ir1–In2O3 SAC for CO2 hydrogenation with a high C2H5OH selectivity of 99% and a turnover frequency (TOF) value of 481 h–1 31. It is found that CO2 can be activated into CO* on the isolated Ir atom, and then couples with the methoxide adsorbed on the In2O3 to form C–C bond. Zhang and co-workers fabricated a stable Ir-based SAC for CO2 hydrogenation to formate. Detailed characterizations and theoretical calculations indicate that such Ir-based SAC has a similar structure with a homogeneous mononuclear Ir pincer complex catalyst, which can serve as an efficient catalyst for the conversion of CO2 to formate32. Despite this progress, CO2 hydrogenation to value-added liquid products over SACs is still far away from the satisfaction of chemical industry because of the following issues: (1) it is challenging to precisely regulate the selectivities of products during CO2 hydrogenation over SACs, (2) the poor stability of SACs at high temperature and pressure needs to be overcome, and (3) the intrinsic structure-performance relationship of SACs has not been clearly understood yet. It is thus highly desired to develop efficient SACs with enhanced activity, selectivity and stability for CO2 hydrogenation under mild conditions.

In this work, we used C3N4, a functional support with abundant N sites, to anchor Cu single sites for CO2 hydrogenation at a low temperature range of 70–150 °C. By simply altering the treatment conditions, the coordination structures of Cu single atoms can be regulated, for instance, Cu–N4 and Cu–N3 structures. Experimental results imply that Cu–N3 SAC displays a CO selectivity of 94.3%, while Cu–N4 SAC exhibit a CH3OH selectivity of 95.5% at a CH3OH productivity of 4.2 mmol g–1 h–1 at 150 °C. Impressively, the CH3OH productivity for Cu–N4 SAC has surpassed the state-of-the-art Cu-ZnO/Al2O3 by 3.2 times (1 mmol g–1 h–1). Moreover, Cu–N4 SAC shows promising stability without decay of CH3OH productivity and changes in the structure after 5 consecutive cycles. Further studies imply that CO2 hydrogenation follows the formate pathway to form CH3OH on Cu–N4 SAC and the RWGS pathway to form CO on Cu–N3 SAC, respectively. This work not only provides an efficient SAC for CO2 hydrogenation to form CH3OH under mild conditions, but also sheds interesting light on studying the structure-performance relationship of catalyst.

## Results and discussion

### Structure analysis of catalysts

C3N4 supported Cu SACs were synthesized via the thermal pyrolysis of melamine (see experimental section for details). The physicochemical properties of catalysts including the compositions (determined by inductively coupled plasma-atomic emission spectroscopy (ICP-AES) measurement) are listed in Supplementary Table 1. Results from the N2 adsorption–desorption measurement suggest the generation of mesopores in Cu–N4 and Cu–N3 SACs. Consequently, the surface areas of pristine C3N4, Cu–N4 SAC and Cu–N3 SAC are 0.70, 22.5, and 20.8 m2/g, respectively (Supplementary Fig. 1). Transmission electron microscopy (TEM) image shows that the as-prepared C3N4 displays a morphology of nanosheet (Supplementary Fig. 2). C3N4 supported Cu SACs with tailored coordination environment, namely, Cu–N4 and Cu–N3 SACs, were generated by depositing CuCl2 onto C3N4 under specific conditions (Fig. 1a, b). After loading Cu on C3N4, no obvious changes in the morphologies of Cu–N4 and Cu–N3 SACs are observed (Fig. 1c, f). Scanning transmission electron microscopy energy-dispersive X-ray spectroscopy (STEM-EDS) images indicate that Cu, C, and N are well dispersed in the Cu–N4 and Cu–N3 SACs (Fig. 1d, g). Moreover, the atomic ratios of Cu in Cu–N4 and Cu–N3 SACs are 12.1% and 13.1%, respectively (Supplementary Fig. 3), which are close to those values from ICP-AES measurement (Supplementary Table 1). Atomic force microscopy (AFM) images demonstrate that Cu–N4 and Cu–N3 SACs have similar thicknesses with C3N4 (~1.6 nm, Supplementary Fig. 4). Note that the decrease of thickness of Cu–N4 and Cu–N3 SACs compared with the pristine C3N4 might be attributed to the delamination of C3N4 during the synthesis of SACs. Moreover, spherical aberration TEM (SA-TEM) measurement was employed to study the states of Cu atoms in different samples. It is found that Cu atoms, which are presented as dense white dots, display as single state in the Cu–N4 and Cu–N3 SACs (Fig. 1e, h). On the other hand, Cu NPs with a mean size of 3 nm are observed in the C3N4 supported Cu nanoparticles (Cu NPs/C3N4) (Supplementary Fig. 2b). The structures of Cu–N4 and Cu–N3 SACs were studied by X-ray diffraction (XRD) measurement. In the XRD pattern of C3N4, typical peaks appear at 2θ of 13.1° and 27.3°, which are ascribed to (100) and (002) facet, respectively (black curve in Supplementary Fig. 5). No obvious differences in the XRD patterns of Cu–N4 and Cu–N3 SACs are attributed to the isolated state of Cu atom (red and blue curves in Supplementary Fig. 5). By contrast, the typical peaks of Cu NPs appear at 43.3°, 50.4°, and 74.1° in the XRD pattern of Cu NPs/C3N4, which correspond to the characteristic peaks of Cu (111), (200), and (220) facets, respectively (Supplementary Fig. 6). Moreover, the peaks at 3145, 1639, 1145, and 798 cm–1 in the Fourier transform infrared spectroscopy (FTIR) spectrum of C3N4 can be indexed as the stretching vibrations of triazine ring, –NH2, –OH, and the bending vibration of triazine ring, respectively, and the presence of characteristic peaks of C3N4 in the FTIR spectra of Cu–N4 and Cu–N3 SACs suggested that the structure of C3N4 is largely maintained after Cu loading (Supplementary Fig. 7).

X-ray photoelectron spectroscopy (XPS) measurement was further used to reveal the surface properties of SACs. It is found that the surface of Cu NPs is partially oxidized, as presented by the shoulders at 934.7 and 955.0 eV in the XPS spectrum of Cu NPs/C3N4 (Fig. 2a). Considering the appearance of Cu peaks in the XRD pattern of Cu NPs/C3N4 (Supplementary Fig. 6), the peaks at 932.5 and 952.2 eV are ascribed to 2p1/2 and 2p3/2 peaks of Cu0. For Cu–N4 SACs, the characteristic peaks of Cu+ 2p2/3 and 2p1/2 appear at 932.5 and 952.2 eV, respectively (Note that the peaks of Cu0 and Cu+ locates at similar positions in XPS spectrum). For Cu–N3 SACs, in addition to the features of Cu+, two peaks corresponding to Cu2+ appear at 935.0 and 955.0 eV in the Cu 2p XPS spectrum with a Cu2+/Cu+ ratio of ~2 (Fig. 2a). For N 1s XPS spectra of Cu–N4 and Cu–N3 SACs, the features of pyridinic N, pyrrolic-N and graphic-N are observed at 398.6, 400.2 and 401.3 eV, respectively, which are similar to those of C3N4 (Supplementary Fig. 8). Moreover, the oxidation states of Cu in Cu–N4 and Cu–N3 SACs have been further confirmed by Auger Electron Spectroscopy (AES) spectra, in which the peaks negatively shift to lower binding energies compared with Cu NPs/C3N4 (Supplementary Fig. 9). To further reveal the chemical states of Cu atoms in Cu–N4 and Cu–N3 SACs, diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) measurement using CO as the probe molecule was performed. Two weak bands appear at 2173 and 2119 cm–1 in the CO-DRIFTS spectra of Cu–N4 and Cu–N3 SACs, which are ascribed to CO linear adsorption on Cu atoms (Fig. 2b). Additionally, X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) measurement were performed at Cu K-edge to investigate the local structures of Cu atoms in catalysts. Compared to Cu foil, the absorption edges for Cu–N4 and Cu–N3 SACs shift to higher energy, suggesting the oxidation states of Cu in Cu–N4 and Cu–N3 SACs (Fig. 2c). Moreover, we estimated the mean chemical valences of Cu in Cu–N4 and Cu–N3 SACs based on the XANES spectra of Cu–N4 SAC, Cu–N3 SAC, Cu foil, Cu2O, and CuO. As shown in Supplementary Fig. 10, the mean chemical valences of Cu in Cu–N4 and Cu–N3 SACs are +1.05 and +1.64, respectively. In the EXAFS spectra of Cu2O and CuO, the peaks at 1.55 and 2.23 Å are ascribed to the coordination of Cu–O and Cu–Cu, respectively (Fig. 2d). The absences of Cu–Cu coordination in the EXAFS spectra of Cu–N4 and Cu–N3 SACs suggest that Cu atoms are presented as isolated state (Fig. 2d and Supplementary Table 2). Additionally, wavelet transform (WT) at the Cu K-edge EXAFS analysis further confirms the isolated state of Cu atoms in Cu–N4 and Cu–N3 SACs. As shown in Fig. 2e, the intense Cu–N coordination, which is illustrated by red color, is observed in the WT contour plots of Cu–N4 SAC and Cu–N3 SAC at the radial distances of ~1.45 Å (Y axis). Also, the strong Cu–Cu (~2.2 Å) and Cu–O (1.5 Å) coordination can be observed in the WT contour plots of Cu foil and CuO, which agrees with the results from EXAFS spectra (Fig. 2e).

### CO2 hydrogenation

The performance of CO2 hydrogenation was evaluated in a batch stainless-steel reactor with 10 mg catalyst and 10 mL water. Commercial Cu-ZnO/Al2O3, C3N4, and Cu NPs/C3N4 were used as references. Prior to catalytic tests, the reactor was pre-flushed with N2 for five times. Afterwards, a mixture of H2, CO2, and N2 (72:24:4 vol.%) was pressurized into the reactor with a total pressure of 3.2 MPa. As shown in Fig. 3a, no gaseous (e.g., CO and CH4) and liquid products (e.g., CH3OH) were detected at 150 °C when C3N4 was as catalyst, indicating that C3N4 is inactive under the indicated conditions. When Cu NPs/C3N4 was used as catalyst, CH3OH productivity is 1.8 mmol g–1 h–1 at a CH3OH selectivity of 78.6%. CH3OH productivity significantly increases to 4.2 mmol g–1 h–1 at a CH3OH selectivity of 95.5% for Cu–N4 SAC. In contrast, CO selectivity reaches 94.3% with a CO productivity of 2.5 mmol g–1 h–1 when Cu–N3 SAC was used as catalyst. It is noted that the productivity and selectivity of CH3OH for Cu–N4 SAC are much higher than those of Cu NPs/C3N4 (Supplementary Fig. 11) and commercial Cu-ZnO/Al2O3 (Supplementary Fig. 12). Moreover, CH3OH and CO productivities reach 4.2 and 2.5 mmol g–1 h–1, respectively, when the reaction time was extended to 3 h. As shown in Fig. 3c, when Cu–N4 SAC was used as catalyst, the CH3OH productivity increases from 0.6 to 4.2 mmol g–1 h–1 with the increased temperature from 70 to 150 °C, whereas the selectivity of CH3OH is kept at ~95%. On the other hand, the CO productivity increases from 0.2 to 2.6 mmol g–1 h–1 for Cu–N3 SAC when the reaction temperature was elevated from 70 to 150 °C (Fig. 3c). We also evaluated the catalytic performance at the temperature range from 70 to 150 °C (Supplementary Figs. 13, 14). It is found high temperature favors CO2 hydrogenation over Cu–N4 (to produce CH3OH) and Cu–N3 (to produce CO) SACs. The above results indicate that Cu–N4 SAC can serve as highly active and selective catalyst for CO2 hydrogenation to CH3OH, while Cu–N4 SAC favors for CO formation (Fig. 3b). Moreover, the kinetics for CO2 hydrogenation were studied for further evaluating the catalytic performance. Arrhenius plots indicate that the apparent activation energies (Ea) for the commercial Cu-ZnO/Al2O3, Cu NPs/C3N4, Cu–N4 SAC, and Cu–N3 SAC are 130.6, 130.3, 50.2, and 86.6 kJ mol–1, respectively (Fig. 3d). The smallest value of Ea for Cu–N4 SAC indicates that Cu–N4 SAC can serve as the most active catalyst for CO2 hydrogenation, which is consistent with our experimental observations. We further calculated the reaction orders with respect to CO2 and H2 for Cu–N4 SAC, and the reaction orders of CO2 and H2 are 1.18 and 1.19, respectively (Supplementary Fig. 15). Note that the catalytic performance of Cu–N4 SAC has surpassed most of the reported Cu-based catalysts for CO2 hydrogenation (Supplementary Table 3).

Additionally, the stability of Cu–N4 SAC was evaluated in five consecutive cycles at 150 °C. No obvious decay of the CH3OH productivity indicates the promising stability of Cu–N4 SAC for CO2 hydrogenation (Fig. 3e). Besides, the local structures of Cu atoms are largely maintained for the spent Cu–N4 SAC (Fig. 3f, Supplementary Figs. 16, 17). ICP-AES results indicate that the content of Cu in spent Cu–N4 SAC is close to that of fresh Cu–N4 SAC (10.2%), confirming the enhanced stability of Cu–N4 SACs for CO2 hydrogenation. Additionally, in situ IR (Supplementary Fig. 18) and XRD measurements were performed for Cu–N4 and Cu–N3 SACs for further studying the stability of Cu–N4 SAC for CO2 hydrogenation. No obvious changes are observed in the IR spectra collected at 3 MPa and 150 °C with different reaction times, suggesting the promising stabilities of catalysts for CO2 hydrogenation (Supplementary Fig. 19). In particular, the peak at 3170, 1650, 1175 and 808 cm-1 are ascribed to the vibrations of triazine ring stretching, –NH2, –OH stretching, and triazine ring bending, respectively. Besides, in situ XRD patterns were collected at different temperature and in the presence of H2/CO2/N2 (72:24:4 vol.%) with a total pressure of 0.03 MPa. As shown in Supplementary Fig. 20, the peak at 27.3° can be indexed as the (002) facet of C3N4, and no obvious changes were observed when the temperature was increased from 40 to 150 °C, further confirming the enhanced stability of Cu–N4 SAC. Note that the slight decrease of productivity in 5 cycles might be attributed to the catalyst loss during centrifugation.

### Catalysis mechanism studies

To understand the distinct performance and catalysis mechanism of Cu–N4 and Cu–N3 SACs towards CO2 hydrogenation, in situ DRIFTS measurement was performed by exposing different catalysts to a mixture of CO2 and H2 (1:3 vol.%) at 0.1 MPa and 150 °C. No peaks were observed in the DRIFTS spectra before introducing CO2 and H2 onto Cu–N4 and Cu–N3 SACs (red curves in Supplementary Fig. 21a, b). After the introduction of CO2 and H2, several peaks appear at 3265, 2964, 2848, 2514, 2014, and 1917 cm–1 in the DRIFTS spectrum of Cu–N4 SAC, which can be indexed as the characteristic peaks of CH3O* (Supplementary Fig. 21a). In contrast, two weak peaks are observed at 2114 and 2176 cm–1 in the DRIFTS spectrum of Cu–N3 SAC, which are ascribed to CO* adsorption (Supplementary Fig. 21b). Additionally, given the significant differences in the coordination environment of Cu–N4 SAC (consists of Cu+) and Cu–N3 SAC (consists of Cu+ and Cu2+), we further studied the effects of Cu2+/Cu+ ratio on the performance of CO2 hydrogenation. Considering H2 can reduce Cu2+ into Cu+, we thus treated Cu–N3 SAC at 500 °C for different times. In particular, the Cu2+/Cu+ ratios are 1 and 0.5 when Cu–N3 SAC was treated at 500 °C for 15 and 30 min, respectively. When the Cu2+/Cu+ ratio is increased from 0.5 to 2 (Supplementary Fig. 22), CH3OH productivity significantly decreases from 2.8 to 0.25 mmol g–1 h–1, while CO productivity significantly increases from 0.5 to 2.6 mmol g–1 h–1 (Supplementary Fig. 23). The aforementioned results imply that the catalytic performance of CO2 hydrogenation was strongly depended on the coordination environment of Cu atoms.

Density functional theory (DFT) calculations were further performed to confirm the stable configurations of Cu single atoms on C3N4 and the corresponding mechanism for catalysis. We choose the buckled C3N4 structure as the substrate, because it is more stable than the planar one, which has been fully discussed in the previous studies33,34,35. EXAFS analysis indicates that the coordination numbers of Cu–N for Cu–N3 and Cu–N4 SACs are ~3.1 and 3.7, respectively (Supplementary Table S2), we thus selected two models to mimic the structures of Cu single atoms on C3N4: 1) Cu atom locates at the hole site of C3N4 (Supplementary Fig. 24). After structural relaxation, Cu atom will not locate at the center of the hole but slightly moves to one side. Cu atom strongly coordinates with two N atoms with the bond length of 2.18 Å and 2.16 Å, and two longer Cu–N bonds of 2.28 Å and 2.38 Å are also formed. In this model, Cu atom is coordinated with 4 N atoms to form structure of Cu–N4, which is consistent with the results from EXAFS measurement. The adsorption energy of Cu is 2.98 eV, which is higher than the adsorption energy of Pd single atoms on C3N4 (2.17 eV) in previous work36, indicating that Cu atoms can be firmly anchored on C3N4 in single state; 2) Cu atom replaces one of C atoms in C3N4 to form three Cu–N bonds (Supplementary Fig. 25). In this model, Cu is coordinated with 3 N atoms, and the length of Cu–N bond is as short as ~1.80 Å, indicating a strong coordination between Cu and N atoms. To confirm the thermodynamic stability of our two models, we have performed the molecular dynamics (MD) calculations at the temperature of 450 K for totally 10 ps with a time step of 1 fs for both Cu–N3 and Cu–N4, respectively. From our results, we can observe the fluctuation of the total energy with the time evolution, as shown in Supplementary Fig. 26. The energies are oscillating near the equilibrium state and there is no structural disruption throughout simulations, indicating the two models are thermodynamically stable. Furthermore, we estimated the Bader charges of Cu in the above two models. It is found that the Bader charge of Cu in Cu–N3 model is +0.95e, which is larger than that of Cu–N4 model (+0.73e), suggesting that more electrons transfer from Cu to neighboring N atoms in Cu–N3 model compared to Cu–N4 model (Supplementary Fig. 27). Consequently, the valence of Cu in Cu–N3 is higher than that in Cu–N4, which is consistent with the results from XANES measurement. Additionally, we analysed the projected density of states (PDOS) of Cu in Cu–N3 and Cu–N4 models to study the electronic structures. It is found that the 3d orbital of Cu atom strongly couples with the 2p orbital of N atom, which confirms the strong coordination between Cu single atoms and C3N4 (Fig. 4a, b). Compared to Cu–N4 ( + 0.73e), the orbital coupling between Cu and N atoms for Cu–N3 (+0.95e) shifts to lower energy levels (far away from the Fermi level), indicating a relatively higher valence state of Cu, which is in agreement with characterizations and Bader charge analysis. We finally simulated the K-edge XANES spectra of Cu–N4 and Cu–N3 SACs based on the selected models. As shown in Supplementary Fig. 28, it is found that both the simulated XANES spectra of Cu–N4 and Cu–N3 models display similar features with those of the original ones, further indicating that the selected Cu–N4 and Cu–N3 models are representative for investigating the mechanism for catalysis.

Moreover, DFT calculations were conducted to investigate the reaction pathways of CO2 hydrogenation on Cu–N4 and Cu–N3 SACs. Based on the previous reports19,20, the major reaction pathway for CO formation is the reverse-water-gas-shift (RWGS) reaction via carboxyl (*HOCO) intermediate, while the formation of CH3OH involves a formate pathway. We then systematically investigated all these possible pathways on Cu–N4 and Cu–N3 SACs. On Cu–N4, DFT calculations imply that the first H prefers to bond with C atom of CO2 to form *HCOO, instead of forming the *HOCO intermediate (the free energy for *HOCO formation is higher than that for *HCOO formation by 0.30 eV at the experimental temperature). *HCOO intermediate then undergoes hydrogenation to form *HCOOH and *H2COOH. The formation of *H2CO is energetically allowed through the cleavage of C − O bond in *H2COOH with the free energy decrease by 0.36 eV (Fig. 4c). In contrast, DFT calculations indicate that the formation of *HCOO and *HOCO completely occurs on Cu–N3 SAC, and the free energies for *HCOO and *HOCO are 0.59 and 0.63 eV, respectively (Fig. 4d). However, the formation of *H2COOH via *HCOOH hydrogenation is not allowed by giving a high free energy barrier of 0.75 eV. Consequently, *HOCO is converted into CO via the RWGS pathway. On the other hand, the adsorption energy of CO on Cu–N3 SAC is only –0.36 eV, indicating that CO can desorb from the surface of catalyst with a free energy decrease of 0.66 eV (Fig. 4d). Based on the aforementioned results, we conclude that CO2 hydrogenation follows the formate pathway to form CH3OH on Cu–N4 SAC (Fig. 4e) and the RWGS pathway to form CO in Cu–N3 SAC (Fig. 4f).

In summary, we have synthesized C3N4 supported Cu SAC with a high loading amount of 12 wt.% for low temperature CO2 hydrogenation. By simply altering the parameters for treatment, the coordination structure of Cu atoms in SAC can be systematically tuned by altering the conditions for treatment (e.g., Cu–N4 and Cu–N3), which leads to the significant differences in the CO2 hydrogenation. Mechanism studies suggest that the reaction pathways are significantly related to the coordination structure of Cu atoms. In particular, Cu–N4 SAC exhibits a CH3OH productivity of 4.2 mmol g–1 h–1 at a CH3OH selectivity of 95.5% at 150 °C, which has surpassed that of commercial Cu/ZnO/Al2O3 by 3.2 times (1 mmol g–1 h–1). Detailed investigations show that CO2 hydrogenation follows the formate pathway to form CH3OH on Cu–N4 SAC and the RWGS pathway to form CO on Cu–N3 SAC, respectively. This work not only provides an active and selective SAC for low temperature CO2 hydrogenation, but also sheds lights on studying the structure-performance relationship of catalyst.

## Methods

### Chemicals

Copper (II) chloride dihydrate (CuCl2.2H2O, > 99.9%) and melamine were purchased from Shanghai Chemical Reagent company. All chemical reagents were used without further purification. All aqueous solutions were prepared using deionized water with a resistivity of 18.2 MΩ cm–1.

### Synthesis of C3N4

In a typical preparation of C3N4, 500 mg melamine was heated to 550 °C with a rate of 8 °C min–1 in a tube furnace under 5% H2/Ar atmosphere for 4 h. After cooling to room temperature, the products were grinded and collected.

### Synthesis of Cu–N4 SAC

In a typical synthesis of Cu–N4 SAC, 500 mg melamine, 200 mg CuCl2, 5 mL ethanol and 5 mL water were added into a 30 mL vial. The mixture was stirred for 24 h and then centrifuged. The obtained sample was dried in a vacuum oven at 60 °C for 24 h and was thoroughly grinded. Finally, the sample was heated to 550 °C at 8 °C min–1 and maintained at the temperature for 4 h in H2/Ar flow (5 vol.% H2 in Ar). Afterwards, the sample was collected after cooling to room temperature.

### Synthesis of Cu–N3 SAC

For the synthesis of Cu–N3 SAC, 100 mg C3N4 and 23 mg CuCl2 were added into 20 mL ethanol under moderate stirring for 12 h at 80 °C. After that, the product was centrifuged and washed with distilled water and ethanol for three times, followed by drying in vacuum at 70 °C for 24 h. Finally, the sample was heated at 500 °C for 2 h in Ar flow to obtain Cu–N3 SAC.

### Synthesis of Cu NPs/C3N4

For the synthesis of Cu NPs/C3N4, the parameters were same with those of C3N4 supported Cu SACs, except for the final treatment in 0.5 M NaBH4 for 30 min at room temperature.

### Characterizations

X-ray diffraction (XRD) measurement was conducted on an X’Pert-Pro X-ray powder diffractometer equipped with a Cu radiation source (λ = 0.15406 nm). High resolution transmission electron microscopy and elemental line scans were performed on a JEOL-2100F transmission electron microscope at an acceleration voltage of 200 kV. High-angle annular dark-field STEM energy-dispersive X-ray spectroscopy (HAADF-STEM-EDS) were conducted on a JEOL GrandARM300F scanning transmission electron microscope with double Cs correctors at an acceleration voltage of 300 kV. Raman spectra were recorded on a HORIBA HR800 instrument. AFM was collected on a Bruker Dimension. Attenuated total reflectance infrared (ATR-IR) spectra were recorded on a Bruker Vertex 70 spectrometer in the spectral range of 4000–400 cm–1. All the XPS spectra of the Cu NPs were collected by XPS (Thermo Scientific, ESCALAB 250 XI). The concentrations of all the catalysts were determined by the inductively coupled plasma-atomic emission spectroscopy (ICP-AES) (710-ES, Varian). The surface properties like Brunauer–Emmett–Teller surface area, Barrett-Joyner-Halenda (BJH) porosity, pore volume of materials was examined by N2 adsorption-desorption isotherms at 77 K by using Autosorb Quantachrome 1MP instrument. The in situ Fourier transform infrared (FTIR) spectra of CO adsorption were measured by a Nicolet 6700 equipped with the liquid nitrogen cooled mercury-cadmium-telluride detector (MCT). Quick-scanning X-ray absorption spectroscopy beamline using bending magnet at Taiwan Photon Source (TPS).

### Catalytic tests

The CO2 hydrogenation was performed in a 60 mL stainless-steel autoclave. The commercial Cu-ZnO/Al2O3 was used as a reference, which was pre-treated in 10%H2/Ar at 250 °C for 4 h. After the additions of 10 mL H2O and 10 mg catalysts into a Teflon inlet, the autoclave was pressurized with CO2 (0.8 MPa) and H2 (2.4 MPa). The reaction was performed at 150 °C with moderate magnetic stirring for 3 h. After completion of the reaction, the gaseous mixture was analyzed using a gas chromatograph (Shiweipx GC-7806) equipped with a GDX-502 column connected to a thermal conductivity detector. The liquid mixture was collected by centrifugation at 14000 × g for 3 min. 10 μL isopropanol was introduced into 1 mL reaction mixture as an internal standard. The liquid mixture was analysed using a gas chromatograph (Persee G5) equipped with a KB-5 column connected to a flame ionization detector. The tests were repeated three times for each catalyst.

### Theoretical methods

DFT calculations were performed by using the Perdew–Burke–Ernzerhof (PBE) functional37,38 within the generalized gradient approximation implemented in Vienna Ab Initio Simulation Package39,40. The projector-augmented wave method41,42 was applied to describe the electron–ion interactions, and the D3 Grimme’s method was employed to correct van der Waals interaction43. We used a plane-wave cutoff energy of 520 eV and the Gaussian smearing with a width of 0.05 eV. Periodic boundary conditions were applied, and more than 15 Å of vacuum space was used to avoid the interaction of the adjacent images. A Γ-centered (3, 3, 1) k-point grid was adopted to sample the Brillouin zone of a 2 × 2 buckled C3N4 supercell with the lattice constant of 6.94 Å. All the structures were fully relaxed until the force components were < 0.02 eV Å–1.

The free energy (G) for each species is expressed as Eq. (1):

$${{{{{\rm{G}}}}}}\,=\,{{{{{{\rm{E}}}}}}}_{{{{{{\rm{DFT}}}}}}}\,+\,{{{{{{\rm{E}}}}}}}_{{{{{{\rm{ZPE}}}}}}}\,-\,{{{{{\rm{TS}}}}}}$$
(1)

where EDFT, EZPE, S, and T represents electronic energy, zero-point energy, entropy and temperature (423 K), respectively. The ideal gas approximation and the harmonic approximation were used for H2O, H2, CO2, CO and CH3OH molecules, and all atomic nuclear motions for adsorbates were considered as harmonic oscillators. The relevant contributions to G were listed in Supplementary Table 4. A series of gas-phase thermochemical reaction enthalpies (Supplementary Table 5) were tested to correct the DFT energies of CO2 and CO due to the inaccuracy of the PBE functional to describe those molecules44. The correction energies of +0.11 and −0.34 eV were applied for CO2 and CO, respectively.

XANES simulations. Cu K-edge XANES spectra were simulated on FDMNES45 based on the cell structure of Cu–N4 (Supplementary Fig. 24) and Cu–N3 models (Supplementary Fig. 25) used for DFT calculations. Cu reference (ICSD no. 64699) was used for calibrating energy. Hubbard energy of 5 eV was used for all the Cu K-edge simulations46, and the Green function employing Muffin-tin potential was applied for the XANES calculations with Athena 0.9.2647.

### Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

All relevant data are available from the authors on request. The theoretical calculation data generated in this study are provided in the Supplementary InformationSource data are provided with this paper.

## References

1. 1.

Schoedel, A., Ji, Z. & Yaghi, O. M. The role of metal-organic frameworks in a carbon-neutral energy cycle. Nat. Energy 1, 16034 (2016).

2. 2.

Jiao, Y., Zheng, Y., Jaroniec, M. & Qiao, S. Z. Design of electrocatalysts for oxygen- and hydrogen-involving energy conversion reactions. Chem. Soc. Rev. 44, 2060–2086 (2015).

3. 3.

Arakawa, H. et al. Catalysis research of relevance to carbon management? Progress, challenges, and opportunities. Chem. Rev. 101, 953–996 (2001).

4. 4.

Nitopi, S. et al. Progress and perspectives of electrochemical CO2 reduction on copper in aqueous electrolyte. Chem. Rev. 119, 7610–7672 (2019).

5. 5.

Mei, Q. et al. Selective utilization of methoxyl group in lignin to produce acetic acid. Angew. Chem. Int. Ed. 56, 14868–14872 (2017).

6. 6.

Haisch, T. et al. Quantification of formaldehyde production during alkaline methanol electrooxidation. Electrochem. Commun. 102, 57–62 (2019).

7. 7.

Raun, K. V. et al. Modeling of the molybdenum loss in iron molybdate catalyst pellets for selective oxidation of methanol to formaldehyde. Chem. Eng. J. 361, 1285–1295 (2019).

8. 8.

Trimpalis, A. et al. NiAu single atom alloys for the selective oxidation of methacrolein with methanol to methyl methacrylate. Catal. Today 355, 804–814 (2019).

9. 9.

Peppley, B. A. et al. Methanol-steam reforming on Cu/ZnO/Al2O3. Part 1: the reaction network. Appl. Catal. A Gen. 179, 21–29 (1999).

10. 10.

Behrens, M. et al. The active site of methanol synthesis over Cu/ZnO/Al2O3 industrial catalysts. Science 336, 893–897 (2012).

11. 11.

Behrens, M. Methanol synthesis over Cu/ZnO/Al2O3: the active site in industrial catalysis. Science 336, 893–897 (2012).

12. 12.

Wu, P. & Yang, B. Significance of surface formate coverage on the reaction kinetics of methanol synthesis from CO2 hydrogenation over Cu. ACS Catal. 7, 7187–7195 (2017).

13. 13.

Kattel, S. et al. Active sites for CO2 hydrogenation to methanol on Cu/ZnO catalysts. Science 355, 1296–1299 (2017).

14. 14.

An, B. et al. Confinement of ultrasmall Cu/ZnOx nanoparticles in metal-organic frameworks for selective methanol synthesis from catalytic hydrogenation of CO2. J. Am. Chem. Soc. 139, 3834–3840 (2017).

15. 15.

Rungtaweevoranit, B. et al. Copper nanocrystals encapsulated in Zr-based metal-organic frameworks for highly selective CO2 hydrogenation to methanol. Nano Lett. 16, 7645–7649 (2016).

16. 16.

Larmier, K. et al. CO2 to methanol hydrogenation on zirconia-supported copper nanoparticles: Reaction intermediates and the role of the metal-support interface. Angew. Chem. Int. Ed. 56, 2318–2323 (2017).

17. 17.

Bando, K. K. et al. In-situ FT-IR study on CO2 hydrogenation over Cu catalysts supported on SiO2, Al2O3, and TiO2. Appl. Catal. A Gen. 165, 391–409 (1997).

18. 18.

Zall, C. M., Linehan, J. C. & Appel, A. M. Triphosphine-ligated copper hydrides for CO2 hydrogenation: structure, reactivity, and thermodynamic studies. J. Am. Chem. Soc. 138, 9968–9977 (2016).

19. 19.

Kattel, S., Yan, B., Yang, Y., Chen, J. G. & Liu, P. Optimizing binding energies of key intermediates for CO2 hydrogenation to methanol over oxide-supported copper. J. Am. Chem. Soc. 138, 12440–12450 (2016).

20. 20.

Kattel, S., Liu, P. & Chen, J. G. Tuning selectivity of CO2 hydrogenation reactions at the metal/oxide interface. J. Am. Chem. Soc. 139, 9739–9754 (2017).

21. 21.

Peng, Y. et al. Molecular-level insight into how hydroxyl groups boost catalytic activity in CO2 hydrogenation into methanol. Chem 4, 613–625 (2018).

22. 22.

Chen, Y. et al. Optimizing reaction paths for methanol synthesis from CO2 hydrogenation via metal-ligand cooperativity. Nat. Commun. 10, 1885 (2019).

23. 23.

Wang, A., Li, J. & Zhang, T. Heterogeneous single-atom catalysis. Nat. Rev. Chem. 2, 65–81 (2018).

24. 24.

Qiao, B. et al. Single-atom catalysis of CO oxidation using Pt1/FeOx. Nat. Chem. 3, 634–641 (2011).

25. 25.

Liu, P. et al. Photochemical route for synthesizing atomically dispersed palladium catalysts. Science 352, 797–800 (2016).

26. 26.

Zhang, G. et al. A general route via formamide condensation to prepare atomically dispersed metal-nitrogen-carbon electrocatalysts for energy technologies. Energy Environ. Sci. 12, 1317–1325 (2019).

27. 27.

Li, H. et al. Synergetic interaction between neighbouring platinum monomers in CO2 hydrogenation. Nat. Nanotechnol. 13, 411–417 (2018).

28. 28.

Yong, X. et al. Revealing the correlation between catalytic selectivity and the local coordination environment of Pt single atom. Nano Lett. 33, 6865–6872 (2020).

29. 29.

Zhang, X. et al. Product selectivity in plasmonic photocatalysis for carbon dioxide hydrogenation. Nat. Commun. 8, 14542 (2017).

30. 30.

An, B. et al. Cooperative copper centres in a metal-organic framework for selective conversion of CO2 to ethanol. Nat. Catal. 2, 709–717 (2019).

31. 31.

Ye, X. et al. Highly selective hydrogenation of CO2 to ethanol via designed bifunctional Ir1–In2O3 single-atom catalyst. J. Am. Chem. Soc. 142, 19001–19005 (2020).

32. 32.

Shao, X. Z. et al. Iridium single-atom catalyst performing a quasi-homogeneous hydrogenation transformation of CO2 to formate. Chem 5, 693–705 (2019).

33. 33.

Kong, N. et al. Single vanadium atoms anchored on graphitic carbon nitride as a high-performance catalyst for non-oxidative propane dehydrogenation. ACS Nano 14, 5772–5779 (2020).

34. 34.

Ji, Y. et al. Heptazine-based graphitic carbon nitride as an effective hydrogen purification membrane. RSC Adv. 6, 52377–52383 (2016).

35. 35.

Wu, H. Z., Liu, L. M. & Zhao, S. J. The effect of water on the structural, electronic and photocatalytic properties of graphitic carbon nitride. Phys. Chem. Chem. Phys. 16, 3299 (2014).

36. 36.

Gao, G., Jiao, Y., Waclawik, E. R. & Du, A. Single atom (Pd/Pt) supported on graphitic carbon nitride as an efficient photocatalyst for visible-light reduction of carbon dioxide. J. Am. Chem. Soc. 138, 6292–6297 (2016).

37. 37.

Perdew, J. P. & Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 45, 13244–13249 (1992).

38. 38.

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).

39. 39.

Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).

40. 40.

Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).

41. 41.

Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558–561 (1993).

42. 42.

Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).

43. 43.

Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010).

44. 44.

Peterson, A. A., Abild-Pedersen, F., Studt, F., Rossmeisl, J. & Nørskov, J. K. How copper catalyzes the electroreduction of carbon dioxide into hydrocarbon fuels. Energy Environ. Sci. 3, 1311 (2010).

45. 45.

Bunau, O. & Joly, Y. Self-consistent aspects of x-ray absorption calculations. J. Phys. Condens. Matter 21, 345501 (2009).

46. 46.

Schnell, I., Czycholl, G. & Albers, R. C. Hubbard-Ucalculations for Cu from first-principle Wannier functions. Phys. Rev. B 65, 075103 (2002).

47. 47.

Ravel, B. & Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat. 12, 537–541 (2005).

## Acknowledgements

This work was financially supported by Ministry of Science and Technology (2017YFA0208200, 2016YFA0204100), the National Natural Science Foundation of China (22025108, 51802206), Natural Science Foundation of Jiangsu Province (BK20180846), Guangdong Provincial Natural Science Fund for Distinguished Young Scholars (2021B1515020081), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and the start-up supports from Xiamen University and Guangdong University of Technology.

## Author information

Authors

### Contributions

X.H. and Y.X. conceived and supervised the research. T.Y., X.W., M.C., and Y.Z. performed the synthesis, most of the structural characterizations, and CO2 hydrogenation tests. X.M. and L.W. performed the DFT simulations. S.Y. helped test spherical aberration electron microscope. C.W.P. helped test EXAFS and XANES data. X.H., T.Y., and Y.X. wrote the paper. T.Y. and X.M. contributed equally to this work. All authors discussed the results and commented on the paper.

### Corresponding authors

Correspondence to Lu Wang, Yong Xu or Xiaoqing Huang.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Peer review information Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Yang, T., Mao, X., Zhang, Y. et al. Coordination tailoring of Cu single sites on C3N4 realizes selective CO2 hydrogenation at low temperature. Nat Commun 12, 6022 (2021). https://doi.org/10.1038/s41467-021-26316-6

• Accepted:

• Published: