Identification of Fenton-like active Cu sites by heteroatom modulation of electronic density.

Developing heterogeneous catalysts with atomically dispersed active sites is vital to boost peroxymonosulfate (PMS) activation for Fenton-like activity, but how to controllably adjust the electronic configuration of metal centers to further improve the activation kinetics still remains a great challenge. Herein, we report a systematic investigation into heteroatom-doped engineering for tuning the electronic structure of Cu-N4 sites by integrating electron-deficient boron (B) or electron-rich phosphorus (P) heteroatoms into carbon substrate for PMS activation. The electron-depleted Cu-N4/C-B is found to exhibit the most active oxidation capacity among the prepared Cu-N4 single-atom catalysts, which is at the top rankings of the Cu-based catalysts and is superior to most of the state-of-the-art heterogeneous Fenton-like catalysts. Conversely, the electron-enriched Cu-N4/C-P induces a decrease in PMS activation. Both experimental results and theoretical simulations unravel that the long-range interaction with B atoms decreases the electronic density of Cu active sites and down-shifts the d-band center, and thereby optimizes the adsorption energy for PMS activation. This study provides an approach to finely control the electronic structure of Cu-N4 sites at the atomic level and is expected to guide the design of smart Fenton-like catalysts.

The Fenton-like process presents one of the most powerful water treatment technologies to tackle persistent organic pollutants resulting from rapid economic development and unsustainable industrial and agricultural expansion (1 –4). The peroxymonosulfate (PMS)-based advanced oxidation process has attracted extensive attention due to its high efficiency at a wide pH range and ease of transport and storage (5 –7). However, the sluggish kinetics of PMS activation during oxidation processes results in prohibitive costs and substantial chemical inputs (8, 9). Therefore, developing efficient catalysts to accelerate the reaction kinetics of PMS is crucial toward efficient catalytic oxidation of recalcitrant organics. Although homogeneous first-row transition metals (Co2+, Fe2+, Cu2+, and Mn2+) generally exhibit remarkable capabilities for PMS activation, they also suffer problems such as poor recyclability and accumulation of sludge (10 –12). Comparatively, heterogeneous catalysts [e.g., transition metal oxides (13, 14), supported nanoparticles (NPs) (15, 16), and carbon-based materials (17, 18)] can be readily recovered and regenerated and are recognized as promising candidates for PMS activation. Nevertheless, the heterogeneity of NPs results in lower utilization efficiency of surface atoms (with 81.6% atoms buried and unavailable for 6-nm nickel NPs) and generally slower reaction kinetics than their homogeneous counterparts (19).

Single-atom catalysts (SACs) featuring utmost atom-utilization efficiency and tunable electronic structure can break the limitations of heterogeneous catalysts in terms of the kinetics and catalytic activity (20, 21). Thus, SACs show a great potential to address the slow reaction kinetics of PMS for the Fenton-like process via maximizing the number of catalytic sites (22). For instance, a single-site Fe catalyst exhibited much faster reaction kinetics toward the degradation of phenol than the Fe NP catalyst, owing to the maximized atomic utilization (23). In addition, the synergetic effect between the atomic center and pyrrolic N site of supports endowed Co SACs with dual reaction sites and high activity for PMS-based oxidation (24). To further accelerate the reaction kinetics of PMS, various strategies have been developed to improve the intrinsic activity of single atomic sites. By controlling the configurations of single atomic sites, PMS was more favorable for adsorption and activation on the CoN2 + 2 site than the CoN4 site (25). Previous work shows that manipulating the electronic structure of single sites plays an essential role in mediating the intrinsic activity (26, 27). It is highly desirable to gain insights into tuning the electronic structure of single-atom sites to achieve superior PMS activation kinetics.

Recent studies demonstrate that the electronic structure of isolated metal sites can be directly modulated by altering the coordinated atom species of the metal centers, favorable for expediting catalytic activity (28, 29). Notably, controlling the long-range interactions with suitable functionalities on the substrate of SACs can be a promising approach for tuning the electronic structure of metal centers (30). Indeed, the kinetic activity of single atomic sites was successfully tuned by the introduction of electron-withdrawing oxidized S groups or electron-donating thiophene-like S species into carbon supports of SACs (31). To this end, nonmetallic heteroatoms offer a substantial potential to serve as electron-withdrawing/donating functionalities on the carbon plane by chemical substitution (32, 33). Specifically, boron (B) with a vacant 2p orbital conjugating with the carbon π system extracts the electrons, while phosphorus (P) with a readily available lone electron pair and low electronegativity is expected to donate electron in graphene (34, 35). With this strategy, incorporating particular heteroatoms (B/P) into the substrate is a possible route to deplete/enrich the electronic density of metal centers, tuning the electronic structure of single sites to promote PMS activation kinetics.

In this work, we designed a versatile strategy to systematically tune the electronic structure of Cu-N4 sites by integrating specific heteroatoms (B/P) into N-doped carbon substrates of Cu SACs. Subsequently, the effect of the controlled electronic features of Cu centers on facilitating PMS reaction kinetics was explored. Here, the heteroatom modified Cu-N4 catalysts were first prepared by using a hydrogen-bonding-assisted pyrolysis approach. Synchrotron X-ray adsorption spectroscopy and the projected density of states (PDOS) analysis verified the successful regulation of the electronic configuration of the Cu-N4 SACs by different heteroatom functionalities. Furthermore, electron paramagnetic resonance and Raman spectra were employed to elucidate the PMS activation mechanism in the Cu-N4/C-B/PMS system. This study opens an avenue to regulating the electronic structures of single active site of SACs to accelerate PMS activation kinetics for pollutant degradation.

Results and Discussion

Heteroatom-Doped Engineering of Cu SACs.

Single-atom Cu catalysts decorated with specific nonmetallic heteroatom (Cu-N4/C-B and Cu-N4/C-P) were prepared through an H-bonding-assisted pyrolysis strategy to study the electronic structure effect of single Cu sites for PMS activation, as illustrated in Fig. 1. Typically, boric/phosphoric acid and carbamide were adopted as building units for self-assembled architecture by H-bonding interaction. The Cu-contained precursor was homogeneously mixed with the supramolecular assembly. SACs consisting of Cu-N4 active sites anchored on heteroatom substituted carbon supports were obtained after pyrolyzing the mixed powder under an Ar atmosphere. The control sample (denoted as Cu-N4/C) without heteroatom was also prepared through a similar procedure but with some modifications (see ).

Fig. 1.

Synthetic illustration and morphology characterizations. (A) Schematic of the preparation strategy for Cu-N4/C-B and Cu-N4/C-P. The color bar indicates the electronic density of Cu-N4 site, electro-rich (blue) and electro-poor (red). (B) HAADF-STEM image and the corresponding EDS mapping images of Cu-N4/C-B. (C) AC HAADF-STEM image and (D) enlarged intensity image of Cu-N4/C-B. (E) Atom-overlapping Gaussian-function-fitting mapping of the square from D, intensity profile along X–Y in D.

The transmission electron microscopy (TEM) images illustrate that Cu-N4/C-B displays a thin graphene-like layered structure, without obvious agglomerates of Cu nanoparticles ( ). As revealed by , the ring-like selected-area electron diffraction (SAED) pattern of Cu-N4/C-B excludes the existence of crystalline Cu. The energy-dispersive spectroscopy (EDS) mapping images from Fig. 1 show that Cu, B, N, and C are distributed uniformly across the entire architecture. To investigate the dispersion properties of Cu species at the atomic level, the aberration-corrected high-angle dark-field scanning TEM (AC HAADF-STEM) measurement was carried out. As shown in Fig. 1 and , the single bright dots corresponding to isolated Cu atoms can be clearly distinguished from the heteroatom-doped carbon matrix. The atomic dispersion of Cu in Cu-N4/C-B is further confirmed by the atom-overlapping Gaussian-function fitting mapping and analysis of intensity profile (Fig. 1). The morphologies of Cu-N4/C and Cu-N4/C-P were also checked, revealing a similar thin graphene-like layered structure ( ). Additionally, the EDS mappings of Cu-N4/C-P reveal the presence of Cu, C, N, and P elements ( ).

The powder X-ray diffraction pattern of Cu-N4/C-B exhibits only one broad peak at approximately 25° indexed to the (002) plane of the graphitic carbon (36), which is similar to these of Cu-N4/C and Cu-N4/C-P ( ). No typical peaks ascribed to crystalline Cu species are observed, which is consistent with the observation from the SAED pattern. Moreover, the Raman spectra of the Cu-N4 catalysts reveal the characteristic D and G bands of conductive carbon materials with similar calculated I D/I G values ( ). Heteroatom modified SACs exhibit slightly larger values than the primary sample, indicating the formation of more defects and disordered structures in Cu-N4/C-B and Cu-N4/C-P. N2 physisorption analysis confirms that all the catalysts have comparable BET surface areas and pore sizes ( ). Furthermore, the Cu loadings, measured by inductively coupled plasma atomic emission spectroscopy, are close to each other in the Cu-N4 catalysts ( ). Considering all the above findings, there are no obvious differences in carbon crystallinity, surface area, or total Cu content for Cu-N4/C, Cu-N4/C-B, and Cu-N4/C-P, which ensures that the effect of heteroatom decoration on the electronic structure of the Cu-N4 catalysts could be investigated independently.

Atomic Structure and Chemical State Analysis of Cu SACs.

Synchrotron X-ray absorption spectra were used to clarify the fine structures of Cu-N4 catalysts. As shown in Fig. 2, the Fourier transform extended X-ray absorption fine structure (FT-EXAFS) spectra of the Cu-N4 catalysts all exhibit one major peak at about 1.5 Å, which can be attributed to Cu-N coordination (37). In contrast to Cu foil and CuO, the absence of a metallic Cu-Cu scattering path (2.2 Å) in the spectra of Cu-N4 catalysts further verifies the atomic dispersion of Cu. To better investigate the atomic configuration of Cu-N4 catalysts, wavelet transform (WT) analysis was performed, which can provide powerful resolution in both R and k spaces and discriminate the backscattering atoms. In line with the FT-EXAFS analysis, the WT contour plots of Cu-N4 catalysts exhibit only one intensity maximum at 4.8 Å−1, which is assigned to the Cu-N bond (Fig. 2 ). For comparison, both Cu foil and CuO counterparts show higher intensity maxima at 7.6 and 6.8 Å−1, corresponding to Cu-Cu and Cu-O-Cu coordination, respectively (Fig. 2 ). We performed EXAFS fitting to analyze the structural parameters of Cu atoms ( ), and the corresponding parameters are listed in . It is observed that the best-fit EXAFS results match quite well with the experiment spectra. Based on the results, Cu species in three single-atom samples exist as the form of isolated Cu-N4 sites, as illustrated by the local atomic structure model (insets in ). The above results strongly suggest that Cu atoms are atomically dispersed in Cu-N4/C, Cu-N4/C-B, and Cu-N4/C-P, without the existence of metal-derived crystalline structures.

Fig. 2.

Atomic local structure and chemical state of Cu-N4 catalysts. (A) Cu K-edge FT-EXAFS spectra of Cu-N4/C-B, Cu-N4/C-P, Cu-N4/C, and reference samples. (B) N K-edge XANES spectra of the Cu-N4 catalysts. (C) The normalized Cu K-edge XANES spectra of the Cu-N4 catalysts and the references (Cu foil and CuO). (D) Cu 2p XPS spectra of Cu-N4/C-B, Cu-N4/C-P, and Cu-N4/C. (E–I) WT-EXAFS plots of Cu foil, CuO, Cu-N4/C-B, Cu-N4/C, and Cu-N4/C-P, respectively.

The electronic structure of the prepared catalysts was first probed by soft X-ray absorption near-edge structure (XANES) analysis. As revealed in Fig. 2, the N K-edge spectra of Cu-N4 catalysts are dominated by three well-resolved resonance peaks attributed to the pyridinic π* (N1), graphitic π* (N2), and C–N σ* (N3) transitions, respectively (38). Notably, the typical peaks of Cu-N4 catalyst shift to lower energy with P incorporation. A similar shift, but in the opposite direction (i.e., to higher energy), was observed in the case of B-decorated Cu-N4 catalyst. The observed changes in the soft XANES spectra indicate a variation in the chemical environment around nitrogen atoms upon heteroatom doping (39, 40), strongly implying a possible change of electronic structure of Cu metal centers. Meanwhile, the B and P K-edge spectra confirm that heteroatoms (B/P) have been successfully doped into carbon substrate ( ).

The Cu K-edge XANES spectra of Cu-N4 catalysts were employed to survey the chemical state of Cu. The energy absorption thresholds of the three Cu-N4 catalysts are situated between that of Cu foil and CuO, indicating the atomically dispersed Cu atoms carry partially positive charges (Fig. 2). Of particular note is that the valence state of single Cu sites was found to rank in the order 2 > Cu-N4/C-B > Cu-N4/C > Cu-N4/C-P > 0, which means the electronic properties of the Cu-N4 sites can be tailored by controlling the heteroatom species in carbon substrate. This observation also shows a good agreement with the Cu 2p high-resolution X-ray photon spectroscopy (XPS) spectra. As shown in Fig. 2, the fitted ratios of Cu (0) to Cu (2+) valence state in three Cu-N4 catalysts show the same trend as the Cu K-edge XANES results, indicating that Cu-N4/C-B possesses a higher chemical valence compared with Cu-N4/C and Cu-N4/C-P. Moreover, N 1s XPS spectra of Cu-N4 catalysts exhibit three characteristic peaks derived from pyridinic N, pyrrolic N, and graphitic N, respectively ( ) (41, 42). These N species exhibit different composition distributions and consequently are associated with the heteroatom incorporation. The successful modification of Cu-N4 catalysts by heteroatoms was further confirmed by B 1s and P 2p XPS spectra of Cu-N4/C-B and Cu-N4/C-P, respectively.

Regulation of Cu Electronic Structure for PMS Activation.

To assess the effect of heteroatom-induced electronic structure regulation of Cu centers on tuning PMS activation activity, the Fenton-like catalytic performance was evaluated for bisphenol A (BPA) removal. As shown in Fig. 3, the pristine Cu-N4/C catalyst achieved less than 60% degradation of BPA in 5 min. Remarkably, the electron-depleted Cu-N4 sites drastically boosted the activation of PMS, with BPA almost completely decomposed using Cu-N4/C-B as the catalyst. However, only about 11% BPA was degraded by PMS when Cu-N4/C-P with electron-enriched Cu centers was employed. In order to provide a better comparison of the catalytic performance, the BPA removal kinetics were then fitted by the pseudo-first-order reactions. It is worth mentioning that the apparent rate constant (k) of Cu-N4/C-B was ∼70 and 5.5 times that of Cu-N4/C-P and Cu-N4/C, respectively (Fig. 3), being one of the most active Cu-based catalysts for PMS activation ( ). Homogeneous Cu catalyst [Cu(II)] showed very poor activity for the activation of PMS ( ). Furthermore, the Fenton-like performance of Cu-N4/C-B is superior to most of the state-of-the-art heterogeneous catalysts ( ). In addition, the comparison with the BPA degradation kinetic curve of BCN without Cu active centers reveals that the Fenton-like activity was mainly ascribed to Cu-N4 sites (Fig. 3).

Fig. 3.

Fenton-like performance for Cu-N4 catalysts. (A) Kinetics of BPA degradation by PMS catalyzed by Cu-N4/C-B, Cu-N4/C, and Cu-N4/C-P within 5 min. (B) Comparison of the rate constant of BPA removal by the three Cu-N4 catalysts. (C) Kinetics of BPA degradation by PMS catalyzed by Cu-N4/C-B and BCN. (D) The relationship between the rate constant and the Cu valence state in the prepared Cu-N4 catalysts. (E) Influence of pH on BPA degradation in the Cu-N4/C-B/PMS system. (F) EPR spectra in the activation of PMS in the presence of Cu-N4/C-B catalyst. (G) Comparison of degradation kinetics under different quenching conditions. (H) Current responses after the sequential injection of PMS and BPA at the Cu-N4/C-B, Cu-N4/C, and Cu-N4/C-P working electrodes. (I) Raman spectra of the Cu-N4/C-B/PMS and Cu-N4/C/PMS systems (ABS: pH buffer). Reaction condition: [BPA] = 20 mg⋅L−1, [PMS] = 0.2 g⋅L−1, catalyst = 0.1 g⋅L−1, T = 298 K, initial solution pH = 6.0.

The dynamic enhancement in the oxidation performance of Cu-N4/C-B indicates the critical role of electronic structure of Cu sites for catalytic PMS activation. Specifically, the reaction rate constant for the synthesized Cu-N4 catalysts was found to correlate well with the Cu valence state calculated from XPS results (Fig. 3). The B modified Cu SACs with a relatively high valence state led to a higher value of k in BPA removal. Cu-N4/C-B catalyst with optimal B content showed the best activity for pollutant oxidation. However, the incorporation of P heteroatom into the substrates of Cu SACs enriched the electronic density of Cu centers (relatively low valence state), resulting in loss of degradation activity. In this regard, electronic structure regulation of metal centers can be an approach to enhance PMS activation kinetics. Moreover, Cu-N4/C-B exhibited satisfying catalytic performance over a broad pH range (Fig. 3), which is beneficial for practical treatment of wastewater.

To identify the reactive oxygen species (ROS) generated during PMS activation, electron paramagnetic resonance (EPR) experiments were conducted. As shown in Fig. 3, the signals for the hydroxyl radical (HO•) and sulfate radical (SO4 •−) were barely detectable in the Cu-N4/C-B/PMS system. Specifically, the characteristic triplet signals of TMP-1O2 adducts for the three Cu-N4 catalysts exhibited enhanced intensity compared to PMS ( ), implying the existence of singlet oxygen (1O2) during PMS activation (43). The EPR analysis suggested that the nonradical ROS was the dominant oxidant in the Cu-N4/C-B/PMS system. Radical quenching experiments were carried out to further identify the active species during PMS activation. As shown in Fig. 3, excessive ethanol (EtOH) served as radical scavengers for SO4 •−/•OH (18, 44), with slight inhibition on BPA removal, indicating that both SO4 •− and •OH contributed little to BPA degradation. Both benzoquinone (BQ) and sodium azide (NaN3) demonstrated noticeable inhibition effect on BPA degradation, suggesting the potential role of 1O2 (45, 46). However, the solvent exchange (H2O to D2O) did not result in an acceleration of BPA decomposition ( ). This contradicts the previous findings that D2O is a promoter for singlet oxygenation because the lifetime of 1O2 in D2O is extended up to 10 times (47, 48). Thus, 1O2 was not primarily responsible for BPA degradation under such conditions, indicating the existence of a secondary nonradical PMS activation pathway (49, 50).

In fact, BQ and NaN3 can also be oxidized by the activated PMS complex via a nonradical pathway (17). Therefore, potassium dichromate (K2Cr2O7) was alternatively employed as an electron scavenger to investigate the electron transfer process. BPA decomposition was barely inhibited in the Cu-N4/C-B/PMS upon K2Cr2O7 addition ( ), suggesting that the system was likely based on direct electron transfer on the catalyst surface other than in solution phase (51). Moreover, chronoamperometry measurements were performed (Fig. 3) to verify PMS activation with Cu-N4 catalysts. The injection of PMS caused distinct current jumps for the three Cu SACs, which verifies electron transfer in the Cu SACs/PMS systems and most likely from Cu SACs to PMS. Impressively, Cu-N4/C-B exhibited the maximum intensity of current jump, demonstrating efficient electron transfer between Cu-N4/C-B and PMS. It is well documented that high-valent intermediates are formed through heterolytically cleaving the peroxide O–O bond for transition metal catalyzed systems (52, 53). In this case, after electron transfer from Cu-N4/C-B to PMS, high-valent copper-oxo species [Cu(III)-OH] were generated by heterolytic cleavage of the O–O bond and served as key intermediates in PMS activation. Furthermore, the subsequent injection of BPA led to the opposite direction of current change, implying electron transfer from BPA to Cu SACs, which resulted in BPA degradation.

To further verify the formation of high-valent copper during the PMS activation, Raman spectra were collected. As shown in Fig. 3, the new peak around 614 cm−1 in the Cu-N4/C-B/PMS system was likely due to high-valent copper-oxo species, because this is a characteristic Cu–O stretching vibration band of a metastable Cu(III) (7, 54). Moreover, this featured peak for Cu(III)-OH was also observed in the Cu-N4/C/PMS system, further confirming the presence of high-valent copper intermediate during PMS activation. The peak at 654 cm−1 was associated with the ABS buffer’s contribution. Moreover, the PMS decomposition rate was basically the same in the presence and absence of BPA, confirming that PMS was decomposed after being activated by Cu-N4/C-B. Therefore, the degradation of BPA in Cu-N4/C-B/PMS system displayed a nonradical pathway based on high-valent copper formation ( ). First, PMS was adsorbed and activated on Cu sites, forming Cu(III)-OH intermediate. Subsequently, the pollutant BPA was adsorbed and attacked by the reactive complex. Afterward, Cu-N4 sites were reactivated after the desorption of the oxidized BPA.

Theoretical Study of Electronic Structure of Cu SACs.

To gain a fundamental understanding of the heteroatom-doped engineering effect on Cu-N4 SACs, density functional theory (DFT) calculations were performed to shed light on regulating the electronic structure of Cu active sites for PMS activation. As shown in Fig. 4, the PDOS for the Cu centers, heteroatoms in the substrate and oxygen of PMS adsorbed on the Cu sites were calculated. In apparent contrast to primary Cu-N4/C, the PDOS of Cu sites in Cu-N4/C-B displays a negative shift, indicating a decrease in d-band center. Comparatively, the PDOS of Cu atoms in Cu-N4/C-P moves toward the opposite direction, with an increased d-band center ( ). Additionally, the PDOS for PMS adsorbed on Cu centers of the three SACs shows that Cu-N4/C-B has a relatively weak interaction with PMS compared to Cu-N4/C ( ). On the contrary, the interaction between Cu-N4/C-P and PMS is pretty strong. As evidenced from the charge density analysis (Fig. 4 and ), electron transfer between PMS and Cu-N4 is observed, suggesting the chemisorption of PMS on all the Cu SACs.

Fig. 4.

PDOS and charge density differences analyses. (A) PDOS of Cu atom, heteroatoms in the substrate, and oxygen of PMS adsorbed on the Cu center (EF is marked in each graph with the black dashed line). (B) The calculated electron density difference diagrams of Cu-N4/C-B with PMS adsorbed on the Cu-N4 site. (C) The adsorption energy values of PMS on Cu-N4/C-B, Cu-N4/C and Cu-N4/C-P catalysts. (D) Relationship between the d-band center and adsorption energy and reaction rate constant for three Cu-N4 samples. (E–G) The optimized structure and the corresponding electron density plots of Cu-N4/C-B, Cu-N4/C, and Cu-N4/C-P catalysts.

The corresponding adsorption energies (Eads) were further calculated, as shown in Fig. 4. Consistent with the PDOS analysis, Cu-N4/C-B possesses a small adsorption energy, while Cu-N4/C and Cu-N4/C-P display strong binding for PMS. It is inferred that Cu-N4/C-B catalyst corresponds to a moderate adsorption energy for PMS activation, which results in its superior Fenton-like activity (Fig. 4). On the other hand, both Cu-N4/C and Cu-N4/C-P might cause poison of the active sites. Moreover, Bader charge analysis (Fig. 4 ) shows that the incorporation of heteroatoms has appreciable influence on electron distribution, especially for the Cu centers. As revealed in , the calculated Cu valence state in Cu-N4/C-B increases markedly in contrast to that of Cu-N4/C, while Cu-N4/C-P displays the opposite behavior. The trend of Cu valence state is qualitatively consistent with the XANES and XPS results. Therefore, heteroatom functionalization can serve as an effective strategy to optimize the electronic structure of Cu-N4 sites, boosting the activation of PMS.

Conclusions

In summary, we systematically investigated the modulation of the electronic structure of Cu-N4 catalytic sites by heteroatom-doped engineering for Fenton-like oxidation at the atomic level. A series of atomically dispersed Cu-N4 catalysts were successfully prepared by a hydrogen-bonding-assisted pyrolysis strategy. Both DFT calculations and experimental investigations indicate that the electronic density of the active Cu centers is well controlled via the long-range interaction with heteroatoms. The modified catalyst can readily achieve high improvement of PMS activation kinetics. Importantly, the electron-depleted Cu-N4/C-B catalyst induces optimized adsorption energy for PMS with an increase of oxidation activity, whereas electron-rich Cu-N4/C and Cu-N4/C-P display strong binding of PMS, causing the poison of active sites. This work would provide a deep insight into the electronic structure regulation of metal centers and structure–activity relationship at the atomic level, which could be helpful to develop advanced Fenton-like catalysts.

Materials and Methods

Materials.

Copper (II) acetate [Cu(OAc)2], 2,2′:6′,2′′-terpyridine, boric acid, phosphoric acid, glucose, and urea were purchased from Sinopharm Chemical Reagent Beijing Co., Ltd. BPA and polyether F127 were obtained from Sigma-Aldrich Chemical Co., Ltd. KHSO5·0.5KHSO4·0.5K2SO4 (PMS), 5,5-dimethyl-1-pyrroline-N-oxide (DMPO), and 2,2,6,6-tetramethyl-4-piperidinol (TMP) were bought from Alfa Aesar Co. Inc. All the chemical reagents were used as received without any other purification.

Synthesis of Cu-N4/C-B Catalysts.

In a typical synthesis of Cu-N4/C-B catalyst, first the Cu precursor was prepared by mixing Cu(OAc)2 (1 mmol) and 2,2′:6′,2′′-terpyridine (1 mmol) in 10 mL of tetrahydrofuran under stirring for 24 h. The resultant precipitate was centrifuged, washed, and vacuum-dried. Subsequently, 50 mmol of carbamide and 1.6 mmol of boric acid were dispersed in 10 mL deionized (DI) water under ultrasonication for 15 min. Then, 0.034 mmol of Cu-contained precursor and 2 mmol of F127 in ethanol (10 mL) were injected into the above solution and the resulting mixture was stirred for 4 h. Then, the H-bonded assemblies were obtained by rotary evaporation of the solvent. Afterward, the sample were pyrolyzed in the Ar atmosphere, maintaining 800 °C for 2 h at a heating rate of 5 °C/min. After being naturally cooled to room temperature, a black powder sample can be obtained. The Cu SACs with B content gradient were prepared by a similar procedure except for changing the amount of boric acid to 0.8 and 2.4 mmol. They are denoted as Cu-N4/C-B-0.8 and Cu-N4/C-B-2.4, respectively.

Synthesis of BCN.

For the BCN catalyst preparation, 50 mmol of carbamide and 1.6 mmol of boric acid were dispersed in 10 mL DI water under ultrasonication for 15 min, then 2 mmol of F127 in ethanol (10 mL) were injected into the above solution and the resulting mixture was stirred for 4 h. Then, the H-bonded assemblies were obtained by rotary evaporation of the solvent. Afterward, the sample was pyrolyzed in the Ar atmosphere, maintaining 800 °C for 2 h at a heating rate of 5 °C/min. After being naturally cooled to room temperature, a black powder sample can be obtained.

Synthesis of Cu-N4/C Catalysts.

Typically, for the synthesis of Cu-N4/C catalyst, 50 mmol of carbamide and 0.56 mmol glucose were dispersed in 10 mL DI water under ultrasonication for 15 min. Then, 0.034 mmol of Cu-contained precursor and 2 mmol of F127 in ethanol (10 mL) were injected into the above solution and the resulting mixture was stirred for 4 h. Then, the H-bonded assemblies were obtained by rotary evaporation of the solvent. Afterward, the sample were pyrolyzed in the Ar atmosphere, maintaining 800 °C for 2 h at a heating rate of 5 °C/min. After being naturally cooled to room temperature, a black powder sample can be obtained.

Synthesis of Cu-N4/C-P Catalysts.

In brief, for the synthesis of Cu-N4/C-P catalyst, 50 mmol of carbamide and 100 μL (1.6 mmol) phosphoric acid were dispersed in 10 mL DI water under ultrasonication for 15 min. Then, 0.034 mmol of Cu-contained precursor and 2 mmol of F127 in ethanol (10 mL) were injected into the above solution and the resulting mixture was stirred for 4 h. Then, the H-bonded assemblies were obtained by rotary evaporation of the solvent. Afterward, the sample were pyrolyzed in the Ar atmosphere, maintaining 800 °C for 2 h at a heating rate of 5 °C/min. After being naturally cooled to room temperature, a black powder sample can be obtained. The Cu SACs with P content gradient were prepared by a similar procedure except for changing the amount of phosphoric acid to 0.8 and 2.4 mmol. They are denoted as Cu-N4/C-P-0.8 and Cu-N4/C-P-2.4, respectively.

Synthesis of PCN.

For the PCN catalyst preparation, 50 mmol of carbamide and 100 μL phosphoric acid were dispersed in 10 mL DI water under ultrasonication for 15 min, then 2 mmol of F127 in ethanol (10 mL) were injected into the above solution and the resulting mixture was stirred for 4 h. Then, the H-bonded assemblies were obtained by rotary evaporation of the solvent. Afterward, the sample were pyrolyzed in the Ar atmosphere, maintaining 800 °C for 2 h at a heating rate of 5 °C/min. After being naturally cooled to room temperature, a black powder sample can be obtained.