Catalyst preparation and characterization
We synthesized CuBx samples and characterized their structures. To prepare the CuBx catalysts with varying B content, we employed a one-step wet chemical reduction method (Fig. 2a and S2-3).15 Varying the molar ratio of CuCl2 and NaBH4 solutions allowed us to control the B content. Inductively coupled plasma mass spectroscopy (ICP-MS) confirmed the presence of B and its tunability (Table S1). Interestingly, the B/Cu mass ratio (wt. %) increased in the bulk phase with increasing B content, while the B/Cu mass ratio in the X-ray photoelectron spectroscopy (XPS) fine spectra showed a different trend, showing a maximum of 5.4 wt. % at CuB2, and low to 0.6 wt. % at CuB5 (Fig. 2b and S4).12 In order to accurately measure the bulk-phase boron concentration of CuBx samples, we measured the dissolution time-dependent boron concentration, whose depth of detection was positively correlated with the dissolution time in dilute nitric acid solution indicating surface enrichment of B at lower B doping levels and subsidence into the Cu bulk phase at higher B levels (Fig. 2c).27
Similarly, XPS was used to study the oxidation state and chemical composition of the sample surface.28 In the corresponding B 1s spectra (Fig. 2e), an expected increase in B binding energy with increasing B content was observed. However, it was also accompanied by a relative decrease in the surface B peak with increasing B content. In the Cu 2p3/2 fine XPS, an attributed peak of CuO is observed at about 933.7 eV, indicating an enhancement of the electron-binding energy of Cu, i.e., an enhanced surface oxidation state. Similarly, the powder X-ray diffraction (XRD) pattern exhibits a transition from the primary Cu crystalline phase to the Cu2O crystalline phase and finally to the CuO crystalline phase with increasing B content (Fig. 2f), which is consistent with the trend of elevated Cu oxidation state. The magnified Cu(111) facet shows a shift Fig. 2. a) Schematic representation of the preparation and components of CuBx samples. b) B/Cu wt.% of bulk and surface of CuBx samples with different B contents measured by XPS and ICP-MS. c) Dissolution time-dependent boron concentration in CuBx samples measured by ICP-MS. d) k-edge XANES spectra of CuB2 and CuB5 samples. e) Boron XPS spectrum of CuBx samples. f) XRD spectrum of CuBx samples. (Inside the inset is a zoomed-in Cu(111) peak that is slightly cheaper as the B content changes) g-j) Wavelet transform (WT) –EXAFS Cu K-edge spectra for CuB5, CuB2, Cu2O and Cu foil as a reference. WT maximum positions are indicated as imaginary lines. ‘+α’ indicates that R is not phase-corrected.of the peak to a lower angle with increasing B content, i.e., accompanied by a stretching phenomenon of the Cu-Cu lattice, which can be inferred that B doping into the Cu in the lattice gap.29,30 We performed X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) analyses (Fig. 2d).31 The Cu K-edge white line peak positions confirm the overall Cu valence state between 0 and 1. However, the leading edge position indicates a higher oxidation state of CuB5 than Cu+, suggesting inhomogeneous oxidation in the structure and an enhanced local Cu oxidation state. This observation is consistent with the near-edge first-order derivative diagram of Cu K-edge (Figure S5), where the first-order derivative of CuB5 exhibits the highest peak at high energy levels. A high-resolution wavelet transform (WT)-EXAFS analysis was performed in both R-space and K-space (Fig. 2g-j). In the lower K-space, the 6.2 Å−1 K peak of the Cu-O coordination signal is not observed, so the O element is Fig. 3. SEM images of Ni-SAC (a), CuB2(b). c) EDS maps of the corresponding Cu, B. Ni elements of the catalysts after hybridization. d, e) TEM image of CuB2 and Hybridized catalyst, respectively. not present in the bulk phase of CuBx. Combined with the previous XRD patterns, we can infer that doping B into the crystalline phase of Cu occupies an O vacancy similar to that of Cu oxide crystals. Meanwhile, increasing the B content also leads to the appearance of new peaks of Cu bonded to lighter elements in the low K region. Moreover, as the peak height increased, it indicated an increase in Cu-B bonding. In the high K space, Cu-Cu coordination is observed, and the Cu-Cu bond in CuBx moves to higher energy space, indicating elevated Cu-Cu bond energy, electron transfer to the surrounding environment, and elevated oxidation state. The results of EXAFS curve fitting for CuB2 and CuB5 showed that the coordination numbers of metal Cu-Cu were 8.0 and 3.2, respectively, and the average bond lengths were 2.58 Å and 2.82 Å, respectively (Fig. 5b, S6-7 and Table s2-3).
We utilized a previously synthesized Ni-SAC catalyst with the "Ni-N2C2" configuration, known for its robust CO2 catalytic performance under neutral conditions.32 High-angle annular dark-field scanning transmission electron microscopy (HAADF-TEM) images provided confirmation of the monodispersion of Ni atoms in the Ni-SAC (see Figure S1). In contrast, CuB2 displayed a distinctive porous morphology, featuring a specific surface area (SBET) of 24.5 m3/g, as measured in Figure S8. This porous structure exhibited a prominent hysteresis loop curve, indicative of the presence of micropores and mesopores around 1.3 nm and 100 nm, respectively (refer to Fig. 3b). Notably, Ni-SAC displayed a thin sheet with an interlayer convolutional morphology (as observed in Fig. 3a). The hybrid catalysts were created by mixing CuBx and Ni-SAC, and energy-dispersive spectroscopy (EDS) elemental mapping confirmed the relative independence of Cu and B elements from C and Ni elements (Fig. 3c). The morphology of CuBx remained largely unaltered both before and after assembly (as depicted in Fig. 3d-e).
Performance for CO-to-ethanol and ethylene electroreduction
The CO2RR performance of the synthesized catalysts was tested in a flow cell in 0.5 M KHCO3 electrolyte (Figure S9). Gas chromatography (GC) analysis revealed the presence of H2, CO, methane (CH4), and ethylene (C2H4), while headspace injection GC analysis detected methanol (CH3OH), ethanol (C2H5OH), and n-propanol in the liquid phase products. Initially, we investigated the effects of different B doping levels on the CO2RR performance (Fig. 4a). Among the catalysts tested, CuB2 demonstrated the highest C2+ faradaic efficiency (FE) of 67% at a voltage of -1.2 V vs. RHE, with a current density exceeding 250 mA cm− 2. 33 CuB2 and CuB5 exhibited specific ethylene-ethanol product transition patterns, resulting in a reduction in FE of C2H4 from 46–34% and Figure 5. a-b) In situ Raman spectra of the CuB2 and Cu samples working at different potentials in 1 M KHCO3.c-d) In situ ATR-FTIR spectra of CuB2 and CuB5 from 2500–1000 cm− 1 working at different potentials in 1 M KOH. OCP means open circuit potential. an increase in ethanol selection by 10%. These observations echoed the differences in B-element content between the surface and bulk phases of CuB2 and CuB5, as characterized earlier. To enhance this effect, we prepared hybrid catalysts with different ratios of Ni-SAC to CuBx (Fig. 4b), achieving ethylene yields up to 72% at a Ni-SAC to CuB2 mass ratio of 0.05, at -1.2 V vs. RHE, with a maximum current density close to 450 mA cm− 2. Ni-SAC was used as a substrate as well as an indirect source of CO, and the separate reduction product monitoring of Ni-SAC showed only H2 and CO, the latter with a Faraday efficiency (FE) of 85 ± 1% at -1.2 V vs. RHE and a partial current density of 120 mA cm− 2 (Figure S10).
To study the reaction kinetics of the catalyst, Nyquist plots were measured at open-circuit potential. Electrochemical impedance spectra (EIS) analysis (Figure S16) revealed that the addition of Ni-SAC effectively reduced the charge transfer impedance (Rct) at the operational voltage. Considering the reactant mass transfer equation ZRe = RΩ + Rct − 2σ2Cd, where σ represents the mass transfer rate from the reactant to the reaction interface during the reaction process, the indirect CO generation by adding Ni-SAC improves the transfer efficiency of CO diffusion to the Cu reaction interface.34–36 Furthermore, we measured the electrochemically active surface area (ECSA) of the catalysts, which indicated the homogeneous dispersion of CuBx by the carbon carriers effectively enhanced the reactive area of the catalysts, leading to increased mass transfer rates and subsequent enlargement of the active surface area (Figure S17).
Despite a decrease in current density and FEC2+ yield with increased Ni-SAC content impregnated with CuB2, the biased generation trend for C2H5OH and CO products remained consistent. We conducted coordinated CO2RR testing of CuB1 − 5 and Ni-SAC ratios with different mass fractions of 0.05, 0.1, and 0.2. Surprisingly, the different Ni-SAC
ratios amplified the tendency for ethanol-ethanol product differentiation, with the highest Figure 6. a) Average bond energies of individual B atoms corresponding to the layers of the Cu lattice when different numbers of B atoms are inserted through them. b) Gibbs free energy step diagram for the transition of *co to ethylene ethanol differentiation intermediates (*HCCOH and *HCCHO). c) The Gibbs free energy of adsorption of *CO at the B site and the Cu site changes when a B atom is inserted into the interlayer of the crystalline phase of Cu as well as without the addition of B. d) Electron density difference plots for the Cu-B bonding configuration and when it adsorbs *CO. e) The Gibbs free energy step diagrams for final C2+ product differentiation at the B site of Cu-L0, and at the Cu site of Cu(111), and of Cu-L1. f) The optimal ethylenic ethanol product differentiation routes and atom migration steps. ethylene selectivity of 72% observed for CuB2-Ni0.05SAC and 46% for CuB5-Ni0.2SAC (Fig. 4d, 4f and S11-14). Further product distribution analysis indicated that in the hybridized Ni-SAC system, FEs of H2 were significantly lower compared to the unhybridized condition, while FEs of CO were elevated and the catalytic activity of CO2 was enhanced (Fig. 4c and S15). This suggests that the tandem system of indirectly generated CO inhibits the hydrogen evolution reaction (HER) and reactive energy barriers for CO generation. Moreover, the results of FEs ratios of C2+/CO indicate that a majority of the formed *CO is inclined to participate in the C-C coupling process, promoting sustained generation of ethylene products with a specific localized CO coverage. However, as the CO source hybridization ratio increased continuously, the CO conversion ability weakened, while CuB5 in a higher local CO environment favored the differentiation pathway towards ethanol products. The bias current densities of the best two samples at different bias voltages are displayed in Fig. 4f. To evaluate the CO2RR performance, we compared the performance of the CuBx-NiySAC catalyst with other top catalysts reported in the literature (Fig. 4e).11,14,18,25,37–45 The stabilized CO2RR properties of CuB2-Ni0.05SAC and CuB5-Ni0.2SAC were tested at a constant voltage of -1.2 V vs. RHE. As shown in Fig. 4g, the bias current densities of the products remained above 200 mA cm− 2 for 10 hours. The morphology of the catalysts remained largely unchanged after a long period of CO2 electrolysis (Figure. S18). In summary, the CuBx-NiySAC catalyst shows great promise for the conversion of CO2 to C2+ products.
In suit characterization of the CuB X
we utilized in situ Raman measurements to investigate the conformational and phase transitions of the CuB2 catalyst during the CO2 reduction reaction (CO2RR), as depicted in Fig. 5a.43,46–49 Initially, at the open-circuit potential (OCP) before CO2RR initiation, the Raman spectrum in Fig. 5b reveals two distinct peaks around 520 cm− 1, corresponding to the A1g mode of Cu2O. These peaks gradually decrease as a negative potential is applied, concomitant with the emergence of a rotational vibrational peak at approximately 280 cm− 1, indicating the involvement of CO in the reduction reaction. Further reduction of the applied potential to -1.0 V vs. RHE leads to nearly complete disappearance of the Cu2O peak on the copper catalyst. Combining this with our previous XPS and XRD characterizations, it is evident that B is present within the oxygen vacancies of CuBx. Additionally, under highly negative potentials, B plays a role in stabilizing the + 1 valence state of CuBx compared to the pure Cu catalyst. Solid-state Raman spectra for both Cu and CuB2, as shown in Figure S19, confirm distinctive Cu-B bond peaks at around 820 cm− 1. In contrast, in situ Raman spectra of Cu do not show any discernible presence of the Cu-B bond. During continuous CO2RR, two characteristic Raman peaks at approximately 1900 cm− 1 and 2130 cm− 1 become evident, corresponding to Cu-surface bridge adsorption and top-adsorption C = O stretching vibrations, respectively. Notably, on the CuB2 catalyst, a novel C = O adsorption site appears at around 2000 cm− 1, suggesting that B functions as an additional adsorption site to stabilize *CO. This intermediate, crucial for the coupling of C = C, is likely instrumental in the formation of C2+ products. Further in situ characterization is necessary to determine the influence of B content on the differentiation of C2+ products.
To comprehend the disparate CO2RR performance of CuB2 and CuB5 catalysts, we employed surface-enhanced Fourier Transform InfraRed-based spectroscopy to identify crucial adsorbed intermediates within the electrocatalytic process.16,50–52 In Fig. 5c and 5d, we observed the corresponding CO2 adsorption peaks at 2360 cm-1 for both samples, and the one for CuB5 is significantly stronger than that for CuB2, which indicates the rapid conversion of CO2 on CuB5. The absorption bands at 2050 cm− 1 correspond to surface-bound CO (*CO) and at 1723 cm− 1 to C = C-coupled *CHO material. Interestingly, the higher *CO peak mirrors that of CuB2 indicating additional *CO sites and that the conversion of *CO on CuB5 is more favorable to *CHO. This observation is consistent with the subsequent pathway of ethanol generation by the Cu-L1 model. The stretching vibrations (Vs) at 1610 cm− 1 *CO = CO, 1290 cm− 1 *CO = COH, and 1084 cm-1 *C-OH indicate that both samples have similar reaction intermediates in their respective reaction pathways. However, CuB5 exhibited more robust Vs for *C-OH, reflecting its enhanced ethanol selectivity. In addition, the bending vibration (Vb) of *C-O-H at 1390 cm− 1 in CuB5 showed significantly stronger absorption peaks compared to CuB2, indicating the importance of this intermediate for ethanol production. Notably, the trend of absorption peak intensities assigned to the *CHO and *COCOH species closely matches the trend of selectivity for the C2+ product, emphasizing the validity of our peak assignments.
Theoretical calculations
We employed first-principles Density Functional Theory (DFT) calculations to investigate the structural stability of B atoms within the bulk phase of Cu. These computations were based on experimentally determined variations in B content under strongly reducing conditions (Fig. 6a). Specifically, we introduced 1, 2, 8, and 16 B atoms into discrete layers of the Cu (4 × 4 × 6) crystalline phase and subsequently determined the optimized geometries, as depicted in Figure S20-21. Remarkably, the atomic radii of B are considerably smaller than those of Cu. Computational analyses reveal that B atoms are densely clustered and distributed within the layers of the Cu crystalline structure. Subsequent to these optimizations, we conducted bond energy calculations for the B atoms. The Fig. 6a showcases the average bond energies of individual B atoms in various layers within the crystalline phase and sheds light on the optimal structural configurations of B atoms at different B contents. Consistent with the experimental findings, at low B atom concentrations, the atoms predominantly adhere to the surface. Conversely, at higher B atom concentrations, clusters are distributed within the bulk phase of Cu.
Based on the experimental results of the changes in product selectivity associated with different levels of B distribution within the Cu crystalline phase, we designed three catalytic configurations; Cu(111), Cu-L0 (low B content) and Cu-L1 (high B content). The structural configurations of both Cu and Cu-L0-5 are listed in Figure S22. According to the free energies of individual B atoms calculated in the Fig. 6a, Cu-L0 is the introduction of individual B atoms on the surface of the crystalline phase of Cu(111), and Cu-L1 is the insertion of eight B atoms through the L1 layer of Cu(111). We then performed first-principles DFT calculations to explore the pathways of C2+ product formation. Previous DFT studies have emphasized the efficacy of the Ni-N2C2 SAC in converting CO2 to CO products.32 Therefore, we first discuss the unique adsorption properties of *CO on these three catalytic configurations. As shown in Fig. 6c, the Gibbs free energy of *CO adsorption on Cu(111) was − 1.605 eV. When single B atom is introduced on the Cu-based surface, the differential density of adsorption charge on the Cu site is concentrated on *CO, which is unfavorable for the retention of *CO, while the B site exhibits a strong adsorption of *CO, which is lower than the Gibbs free energy of *CO adsorption on Cu(111) of 0.65 eV (Fig. 6b). Therefore, the B site on Cu-L0 can be used as an additional adsorption site for *CO for subsequent C2+ product conversion.
To delve further into the impact of B on Cu sites within the Cu crystalline phase, we analyze the pathway from *CO to the intermediate *HCCHO (Fig. 6d). In previous calculations, we addressed the sequential problem of proton charge transfer, manifesting as the optimal reaction pathway involving H addition followed by transition state (TS) coupling.53 Notably, the C = C coupling process, being heat-absorbing and influenced by external factors, serves as the rate-determining step (RDS) in this system. The free energy of this coupling process on the Cu site changes from 0.43 eV to 0.42 eV with the introduction of higher B content, indicating that *CO is more likely to generate C2+ products on Cu sites with CuB5. Interestingly, when we employed the B site on CuB2 for additional *CO adsorption, it exhibited lower energy barriers (*COH,TS and coupling) due to stronger *COL adsorption on B. Experimental in-situ Raman results also indicated enhanced *CO adsorption. Furthermore, the stronger B-CO bonding raised the difficulty of hydrogenolysis of *COH to dehydrate. This resulted in an alternative pathway involving H addition, followed by hydrogenolysis to produce *HCCOH, which had lower energy compared to *HCCHO by 0.33 eV. This suggests the advantage of B as an adsorbate site for *CO in facilitating the differentiation of C2+ products, the corresponding divergent charge differential distributions of *CCO and *HOCCOH are also shown in Fig. 6b.
Next, we needed to continue exploring the intermediate pathways for C2+ product differentiation on three catalytic configurations. Specifically, we focused on the corresponding reaction energies for the generation of ethylene and ethanol products after initiation with *HCCOH and *HCCHO (as shown in Fig. 6e). In the Cu(111) and Cu-L1 configurations with *HCCHO as the promoter at the Cu site, *H tends to be added to C due to the relatively weak Cu-CO bond energy. This leads to an increase in the energy of *HCCHOH on O, and an increase in the free energy of the ethylene pathway by 0.15 eV compared to the ethanol pathway. This suggests that C2+ product partitioning at the Cu site facilitates the protonation of C and subsequent retention of -OH in Cu-C-O, which is essential for subsequent ethanol formation. In the Cu-L1 configuration, the presence of interlayer B elements leads to lattice expansion and the formation of electronically delocalized domains within the Cu crystalline phase, as observed by XRD and XPS experiments. Calculations show lower interstitial adsorption energies, indicating higher stability. In the presence of B, the free energy change of the C2+ product is 0.07 eV lower than that on Cu(111), indicating that B doping enhances the catalytic activity of CO2RR. Notably, the free energy change of the *HCCOH pathway is usually lower than that of *HCCHO, and the subsequent hydrogenation and dehydration processes become easier due to the addition of H on O, resulting in a continuous decrease in the overall stepwise trend. Since H-O- leaves easily with O, the -OH group cannot be favorably retained, ultimately leading to ethylene production.39 The optimal ethylenic ethanol product differentiation routes and atom migration steps are shown by us in Fig. 6f. Throughout the reaction pathway, when B is deposited on the C surface as an adsorption site for *CO, the total free energy change of *CO during ethylene production is minimized and the reaction step shows a decreasing trend. Our experiments confirmed that B was deposited on the surface only at low concentrations, indicating effective but limited reaction sites. Consequently, the C2+ selectivity efficiency increases significantly when an additional source of *CO generation is introduced.