Catalyst design and DFT prediction. Most of Ni-N-C catalysts have been revealed to show very weak binding to COOH* intermediate upon CO2-to-CO coversion.13,16 Therefore, injecting foreign electrons into Ni single atom to tune its d band structure is desirable for enhancing the adsorption of COOH*. Based on this principle, previous studies have found that Ni(i)-N-C shows significantly enhanced CO2RR activity compared to Ni(ii)-N-C.17 Note that metal nanoparticles can be easily formed with single atoms because of the aggregation effect, and the formation of nanoparticles sometimes is inevitable if a high loading of single atom is pursued. In this case, the electron-rich structure of metal nanoparticles can potentially tune the electronic configuration of nearby metal-N-C,18 and enhance its electrocatalytic performance. Given that Ni is a strongly interacting metal with carbon,19 we propose a Ni-N-C on Ni nanoparticles design, and firstly use DFT calculations to predict the electron distribution and charge transfer between Ni nanoparticles and a layer of carbon via the deformation electronic density calculation and Bader analysis. Here face-centered cubic Ni(111) is used as the targeted model because of its lowest formation energies among the different phases.20 As shown in Fig. 1a,b, the first layer of Ni(111) can easily bind with the adjacent carbon, which leads to a total transfer of 2.45 e− to the carbon surface.
We next shift our attention to study the effects of extra electrons on the Ni-N-C structures for CO2RR. Here, two representative structures including coordinatively saturated NiN4 (Fig. S1) and unsaturated NiN3V are considered (Fig. S2),13,21,22 where V stands for vacancy and Ni single atom is the active site.23 COOH* and CO* are intermediates for electrochemical CO2-to-CO conversion.24 On neutral Ni-N-C catalysts, as shown in Fig. 1c,d, the rate-limiting steps are CO2(g)→COOH* on both NiN4 and NiN3V with the overpotential of 1.50 V and 0.84 V, respectively, attributed to the weak binding of COOH* on NiN4 and NiN3V. After injecting extra electrons into the Ni-N-C structure, the binding strength of COOH* and CO* on NiN4 and NiN3V are greatly enhanced with the increase of extra electrons. In comparison with the neutral cases, the stronger binding of COOH* on negatively charged NiN4 and NiN3V reduces the overpotential of CO2RR and enhances their activities. To understand the effects of extra electrons on the enhanced binding strength of COOH* and CO* on Ni-N-C, we then calculated the number of electrons localized on Ni single atom (Fig. 1e). The results show that the number of electrons on Ni atom increases with the injected extra electrons on Ni-N-C. More electrons localized on Ni single atom results in stronger binding of intermediates compared to those on neutral cases.
Given that hydrogen evolution reaction (HER) is a dominant side reaction that is competitive with the CO2RR, we calculate the difference between thermodynamic limiting potentials for CO2RR and HER (denoted as UL(CO2)-UL(H2)), which can reflect the selectivity (Fig. S4).25,26 Figure 1f shows the UL(CO2)-UL(H2) on negatively charged NiN4 and NiN3V, in comparison with the neutral cases. Clearly, negatively charged NiN4 and NiN3V show more positive values of UL(CO2)-UL(H2), corresponding to their higher selectivity for CO production.
Catalysts synthesis and characterizations. We pursued the nanoparticle co-catalysis effect experimentally by one-step solid-state pyrolysis with excessive Ni sources. The schematic diagram of the catalyst design and formation mechanism is shown in Fig. S5a. First, the melamine was decomposed into carbon nitride at ~400°C.27 With the increase of annealing temperature, the Ni nanoparticles catalyzed the carbon nitrite to generate the carbon nanotubes (CNTs),28 meanwhile Ni single atoms were trapped by the nitrogen-doped CNTs via thermal emitting.29,30 The NiSA/NP was obtained after the heat treatment followed by the acid washing to remove the Ni nanoparticles outside of CNTs. Figure 2a shows the morphology of the NiSA/NP catalyst, where large amounts of Ni nanoparticles are confined inside of the nitrogen-doped CNTs. The EDS mapping confirms the presence of Ni and N in the carbon matrix, and evenly distributed Ni nanoparticles at micron-scale (Fig. S5d). The bright dots displayed in Fig. 2b are assigned to the Ni single atoms anchored on CNTs. For comparison, pure NiSA was prepared by the treatment of NiSA/NP with NH4Cl to remove Ni nanoparticles inside of the CNTs.31 As confirmed by TEM (Fig. 2c and Fig. S6) and X-ray diffraction (XRD, Fig. 2d), most of the Ni nanoparticles inside of CNTs were removed after the NH4Cl treatment while Ni single atoms remained. XRD also shows a dominant Ni(111) peak of the catalysts which is consistent with the DFT modelling. The XPS N 1s spectra of NiSA and NiSA/NP exhibit similar pyridinic (398.5 eV), pyrrolic (401.2 eV), graphitic (402.7 eV), and Ni-N (399.4 eV) peaks (Fig. S7b,c), indicative of the identical Ni-N-C structure after the removal of Ni nanoparticles. According to the ICP-OES, the Ni contents in NiSA/NP and NiSA are 8.8 wt% and 4.6 wt%, respectively. The high content of Ni single atom can be attributed to (i) the in-situ formation of Ni-N-C during the growth of CNTs and (ii) the excessive Ni sources provided by the Ni nanoparticles.
The unpaired electron in catalysts was directly observed by EPR spectroscopy (Fig. 2e), where a g-value of 2.17 is obtained for both samples and could be assigned to the unpaired electron in the \(3{d}_{{x}^{2}-{y}^{2}}\) orbital of Ni(i).32,33 The EPR signal intensity of NiSA/NP is above 2 times higher than NiSA, corresponding to the electron-rich structure due to the presence of Ni nanoparticles, consistent with the Bader analysis. We further examined the d band structure of Ni using ultraviolet photoelectron spectroscopy (UPS, Fig. 2f and Fig. S8). No obvious metallic Ni peak can be observed across the Fermi level, which can be attributed to the fully confined Ni nanoparticles inside of CNTs that exceeded the detection limit of UPS (2-3 nm), in line with XPS results (Fig. S7a). Therefore, the shift of near Fermi level signals (0-1 eV) can be attributed to the change of electronic structure of Ni single atoms on the surface of CNTs. The position of NiSA increases in binding energy compared to NiSA/NP (away from the Fermi level). According to d band theory, the downshifted d band of NiSA (grey curve) corresponds to a weakened adsorption strength,34 which is in good accordance with DFT calculation.
The Ni K-edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) measurements were carried out to identify the electronic structure and coordination environment of catalysts. Ni foil and NiPc were used as references for metallic Ni and Ni2+, respectively. The pre-edge curve of NiSA is in the middle of Ni foil and NiPc (Fig. 2g), suggesting a valence state of Ni close to +1.35 Ni(i)-N-C has been revealed as a highly active site for CO2-to-CO conversion.32 More importantly, the pre-edge position of NiSA/NP is between the Ni foil ad NiSA, corresponding to an even lower average valence state than Ni+. The Fourier-transformed k3-weighted spectra of the samples exhibit two major peaks at 1.45 Å and 2.1 Å, which are assigned to Ni-N and Ni-Ni coordination shells, respectively (Fig. 2h). The significantly decreased Ni-Ni peak intensity in NiSA compared with NiSA/NP corresponds to the removal of Ni nanoparticles. Note that the Ni-Ni peak can be easily observed even there are only residual Ni NPs in the catalysts.31 Based on the EXAFS fitting of Ni-N shell on NiSA, the coordination number is determined to be 2.7 which is close to the coordinative unsaturated NiN3V structure (Fig. 2i and Table S2).
Electrochemical CO 2 conversion. The electrochemical testing was firstly conducted in an H-cell using 0.5 M KHCO3 electrolyte and the corresponding calibration curve is shown in Fig. S9. Figure 3a displays the partial CO currents at different potentials, where the NiSA/NP outperforms the NiSA in a wide potential window, achieving a high current of 131 mA cm−2 at -1.0 V vs RHE which is almost doubled compared to NiSA. To the best of our knowledge, this result outperforms most, if not all, CO2-to-CO catalysts in H-cells reported to date (Table S3). The Faradic efficiency of CO (FECO) is shown in Fig. S10, where NiSA/NP can achieve an ultra-high FECO above 99% at -0.8 V. Besides, the onsite potential of NiSA/NP is at around -0.5 V vs RHE, which is approaching to the values predicted from DFT calculation on 2 e− negatively charged NiN3V model. The electrochemically active surface area (ECSA) is evaluated from double-layer capacitance (Cdl) to elucidate the origin of the enhanced catalytic activity on NiSA/NP (Fig. S11).36 The NiSA/NP shows slightly higher Cdl compared with NiSA (1.4 mF vs 1.2 mF), suggesting that the improved catalytic activity is primarily attributed to the intrinsic reactivity of each site, rather than a surface area effect (Fig. 3b). Stability testing is then performed at -0.75 V vs RHE, NiSA/NP shows robust durability for CO2RR during 30 hours of continuous electrolysis, maintaining constant current densities and 99% of the initial FE for CO production (Fig. 3c). This can be attributed to the strong immobilization and close proximity between catalysts and CNTs, which reduce aggregation and bypass undesired side reactions to suppress catalyst deactivation.37
To assess the performance of catalysts at industrial-relevant devices and current densities, we further evaluate the CO2RR in flow cells using gas diffusion electrodes (Fig. S12). Figure 4a displays the CO FEs and partial current densities plotted against the iR-corrected potentials in 1 M KOH flow cells. The NiSA/NP shows a low onsite potential at around -0.25 V vs RHE. Then, the CO partial currents quickly increase to above 200 mA cm−2 at -0.4 V with FECO of 96%, and further reach 346 mA cm−2 at -0.5 V with FECO of 98%. As a comparison, the highest current obtained with NiSA is 263 mA cm−2 at -0.65 V (Fig. S13). The superior CO2RR performance of NiSA/NP in the flow cell can be attributed to (i) the significantly enhanced CO2 mass transport based on gas diffusion electrode and (ii) the catalytic promotion effect in alkaline electrolytes.43
The full cell performance was then tested based on the zero-gap MEA device, which is regarded as a promising prototype for practical application because of the very low cell resistance and robust structure. Fig. 4b exhibits the LSV curves of two different full cells by using IrO2 and NiFe anodes, respectively. With 1 M KOH as the anolyte, the NiSA/NPǀǀNiFe cell shows even better performance than the widely used IrO2 anode. It is worth noting that the Ir-based anode can cost more than 60% of the total price of MEA device,7 so that the development of the low-cost substitution is highly desirable. As shown in Fig. 4c, the NiSA/NPǀǀNiFe cell delivers a jCO of 163 mA cm−2 with FECO of 92% at -2.0 V without iR compensation. With further increase of cell voltage to -2.3 V, the jCO can reach up to 310 mA cm−2 with FECO of ~99%, corresponding to an outstanding full cell energy efficiency of 57% based on the calculation: EE = ΔE0/ΔEApplied × FECO, where EE is the energy efficiency, ΔE0 is the equilibrium full cell voltage (E0CO- E0OER = -0.11 V - 1.23 V = -1.34 V), and ΔEApplied is the applied full cell voltage without iR-correction. Note that certain Ag-based MEA devices can achieve current densities above 600 mA cm−2 but at a much higher cell voltage of -3.5 V with a low energy efficiency of ~35%.38 In comparison, our NiSA/NPǀǀNiFe cell shows the highest energy efficiency and the lowest cell voltage ever reported for CO production at an industrially relevant current of ~300 mA cm−2 (Fig. 4d and Table S4).