Wenxiu
Ma‡
a,
Zhe
Chen‡
b,
Jun
Bu
a,
Zhenpeng
Liu
a,
Jinjin
Li
a,
Chen
Yan
a,
Lin
Cheng
d,
Lei
Zhang
a,
Hepeng
Zhang
a,
Jichao
Zhang
*c,
Tao
Wang
*b and
Jian
Zhang
*a
aKey Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, Xi'an, 710129, P. R. China. E-mail: zhangjian@nwpu.edu.cn
bCenter of Artificial Photosynthesis for Solar Fuels, School of Science, Westlake University, Hangzhou 310024, P. R. China. E-mail: twang@westlake.edu.cn
cShanghai Synchrotron Radiation Facility, Zhangjiang Laboratory, Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai, 201204, P. R. China. E-mail: zhangjichao@zjlab.org.cn
dSchool of Science, Xi'an Polytechnic University, Xi'an, 710048, P. R. China
First published on 15th November 2021
In pursuit of environmental friendliness and high economic efficiency, renewable energy-driven selective acetylene semihydrogenation under ambient conditions is urgently required but remains a grand challenge. Herein, we explore single-atom nickel (Ni) dispersed N-doped carbon (SA-Ni-NC) as a novel electrocatalyst for catalyzing selective acetylene semihydrogenation. In situ electrochemical Raman and theoretical investigations reveal that weak π-adsorption of ethylene on individual Ni sites in SA-Ni-NC facilitates its desorption and thus avoids its overhydrogenation. Eventually, under pure acetylene flow, SA-Ni-NC exhibits a high ethylene faradaic efficiency (FEethylene) of 91.3% and a large current density of −92.2 mA cm−2 at −0.6 V vs. the reversible hydrogen electrode (RHE). Even in crude ethylene stream containing 1% acetylene impurities (1 × 104 ppm), SA-Ni-NC still manifests a high acetylene conversion of 97.4% with a large space velocity (SV) of 2.4 × 104 mL gcat−1 h−1 and a high turnover frequency (TOF) of 22.9 h−1 for each Ni atom.
For enhancing ethylene selectivity, weakening ethylene adsorption to avoid its over-hydrogenation is crucial during acetylene semihydrogenation.17–19 Among different ethylene adsorption modes ethylidyne mode, di-σ-mode, and π-bonded mode, adsorption strength of the π-bonding configuration on one metal atom is the lowest.1,20,21 Meanwhile, isolated metal atoms can suppress the carbon–carbon coupling reaction between two neighboring acetylene molecules. Accordingly, electrocatalysts featuring single-atom active sites are promising for facilitating electrocatalytic acetylene semihydrogenation to ethylene but remain unexplored until now.
In this work, we synthesize single-atom Ni supported on N-doped carbon (SA-Ni-NC) by pyrolyzing nickel doped zeolitic imidazolate framework-8 (ZnNi-ZIF) precursors.22,23 When the SA-Ni-NC serves as an electrocatalyst for acetylene semihydrogenation, it shows an excellent electrocatalytic performance with an ethylene faradaic efficiency (FEethylene) of 91.3%, a partial ethylene current density of −92.2 mA cm−2, and a turnover frequency (TOF) of 3847.8 h−1 at −0.6 V vs. the reversible hydrogen electrode (RHE) under a pure acetylene flow. For crude ethylene flow containing 1 × 104 ppm acetylene, SA-Ni-NC presents a high conversion of 97.4% and a large space velocity of 2.4 × 104 mL gcat−1 h−1. Density functional theory calculations combined with in situ electrochemical Raman analyses demonstrate that the superior electrocatalytic performance of SA-Ni-NC is attributed to weak π-adsorption of ethylene on single Ni atoms, which favors the rapid ethylene desorption and eventually promotes selective semihydrogenation of acetylene to ethylene.
:
1, Ni nanoparticle immobilized N-doped carbon (Ni NP/N-C) was fabricated by increasing the Ni/Zn ratio of the ZnNi-ZIF to 9
:
1 (Fig. S6–S9, ESI†). The specific surface area of SA-Ni-NC was measured to be ∼727.9 m2 g−1 based on Brunauer–Emmett–Teller (BET) adsorption/desorption isotherms (Fig. S10, ESI†).
X-ray photoelectron spectroscopy (XPS) was employed to study the chemical structure of SA-Ni-NC (Fig. S11, ESI†). In the XPS spectra of Ni in SA-Ni-NC, two characteristic peaks at 855.3 and 872.5 eV correspond to Ni 2p1/2 and Ni 2p3/2, respectively.24 The high-resolution XPS spectrum of N 1s revealed the coexistence of pyridinic (398.6 eV), pyrrolic (400.8 eV) and graphitic (403.8 eV) nitrogen species.25 The Ni content in SA-Ni-NC was determined to be ∼2.39 wt% by utilizing inductively coupled plasma analysis (ICP, Table S1, ESI†). Transmission electron microscopy (TEM) was conducted to scrutinize the morphology of SA-Ni-NC and no Ni particles existed (Fig. 1b and S12, ESI†). Fig. 1c and S13d (ESI†) show the corresponding elemental mapping images of C, N, and Ni, evidencing their homogeneous distributions over SA-Ni-NC. The high-angle annular dark field-scanning transmission electron microscopy (HAADF-STEM) showed the atomic dispersion of Ni in the carbon matrix of SA-Ni-NC (Fig. 1d, e and S13a–c, ESI†). The electron energy loss spectroscopy (EELS) analysis was further conducted. As shown in Fig. S14,† the bright atoms are single-atom Ni in SA-Ni-NC.
To further unveil the bonding environment of Ni atoms in SA-Ni-NC, the X-ray absorption fine structure (XAFS) analysis was employed. With the Ni foil and Ni porphyrin complex (NiPc) as references, the Ni K-edge X-ray absorption near-edge structure (XANES) of SA-Ni-NC is depicted in Fig. 2a. Clearly, the XANES spectra of the Ni K-edge indicated that the energy absorption edge of SA-Ni-NC was located between those of Ni foil and NiPc, and close to NiPc, indicating that the valence state of Ni species in SA-Ni-NC was about +2. Furthermore, the Fourier transform (FT) k3-weighted extended X-ray absorption fine structure (EXAFS) curve of SA-Ni-NC is shown in Fig. 2b. The main peak of SA-Ni-NC at 1.41 Å corresponded to the first shell scattering of the Ni–N bonds rather than the Ni–Ni bonds at 2.18 Å. The EXAFS fitting analysis revealed that the isolated Ni atoms were in a Ni–N4 configuration with a Ni–N bond length of 1.83 Å (Fig. 2c and Table S2, ESI†).26–28 The wavelet transform (WT) contour plot of SA-Ni-NC had only one intensity maximum at 5.6 Å−1, which was assigned to the Ni–N bonds rather than Ni–Ni scattering (8.2 Å−1). These results clearly prove the dispersion of Ni–N4 species in SA-Ni-NC (Fig. 2d–f).29,30
For electrocatalytic acetylene semihydrogenation, SA-Ni-NC and Ni NP/N-C were tested using a three-electrode flow cell in 1 M KOH aqueous solution under pure acetylene flow with a flow rate of 20 mL min−1 (Fig. S15–S17, ESI†). As shown in the linear sweep voltammetry (LSV) curve in Fig. 3a, the cathodic current density of Ni NP/N-C at −0.5 V considerably decreases to −31.9 mA cm−2 from −47.8 mA cm−2 when the gas flow switched from Ar to pure acetylene. By contrast, SA-Ni-NC exhibited an increased current density of −55.3 mA cm−2 at −0.5 V in acetylene flow, which was 3.1 times higher than that under Ar. These results undoubtedly reflected the outstanding electrocatalytic activity of SA-Ni-NC toward acetylene. Based on gas chromatography (GC) analyses, the total faradaic efficiencies of gaseous products (H2 and ethylene) were close to 100% and no liquid products were detected (Fig. S18b, ESI†). Noticeably, as shown in Fig. 3b, SA-Ni-NC achieves a FEethylene of 91.3% at −0.6 V, which is dramatically higher than 71.7% for Ni NP/N-C. Meanwhile, SA-Ni-NC also exhibited an enhanced ethylene partial current density (jethylene) over all applied potentials (Fig. 3c). The jethylene of SA-Ni-NC reached −84.2 mA cm−2 at −0.6 V, which was approximately 1.8 times higher than that for Ni NP/N-C. Impressively, according to the inset in Fig. 3c, the calculated TOF of SA-Ni-NC was up to 3847.8 h−1 at a potential of −0.6 V. Then, long-term electrocatalytic acetylene semihydrogenation at −30 mA cm−2 was carried out (Fig. S18c, ESI†). Apparently, over 12 h, SA-Ni-NC presented a stable potential of about −0.5 V and an average FEethylene of ∼95%. For deeply probing the underlying mechanism of superior electrocatalytic acetylene semihydrogenation on SA-Ni-NC, electrochemical active surface area (ECSA) and electrochemical impedance spectroscopy (EIS) measurements were evaluated (Fig. S19, ESI†). The ECSA of SA-Ni-NC and Ni NP/N-C was determined by measuring the electrochemical double-layer capacitance (Cdl), which was derived from the cyclic voltammograms at different scan rates from 20 to 100 mV s−1. As indicated in Fig. S19c,† the Cdl of SA-Ni-NC was about 34.1 mF cm−2, which was much higher than 20.0 mF cm−2 for Ni NP/N-C, suggesting that SA-Ni-NC exposed more active sites. The EIS plot of SA-Ni-NC showed a smaller charge transfer resistance of ∼29 Ω than >120 Ω for Ni NP/N-C, indicating its fast electron transfer kinetics.
To assess the practical application feasibility of SA-Ni-NC, a two-electrode full cell was constructed by integrating NiFe-LDH on Ni foam as the anode for the oxygen evolution reaction. As shown in Fig. S20a,† in 1 M KOH aqueous solution, the current density of the full cell was as high as −100 mA cm−2 at a voltage of 3.25 V under pure acetylene flow. The FEethylene of SA-Ni-NC was ∼99% at a current density of −30 mA cm−2 (Fig. 3d). In particular, during a 45 h long-term durability at −30 mA cm−2, the FEethylene of SA-Ni-NC was steadily >90% (Fig. 3e).
Encouraged by such outstanding semihydrogenation performance of SA-Ni-NC under a pure acetylene atmosphere, its performance in the crude ethylene feedstock containing 1% acetylene impurities (1 × 104 ppm) was further assessed in a flow cell. Limited by the large input flow rate and short gas residence time, more than 70% acetylene impurities still existed in output ethylene flow. Therefore, we customized a flow cell with a 5 × 5 cm2 electrode area as a preliminary demonstration for further evaluating the application possibility of SA-Ni-NC (Fig. S21–S23, ESI†). Fig. 4a shows the contents of acetylene impurities at different current densities. The acetylene concentration in output ethylene flow was considerably decreased to 263.6 ppm at −40 mA (−1.6 mA cm−2), corresponding to an acetylene conversion of ∼97.4%. Accordingly, the space velocity (SV) of SA-Ni-NC was up to 2.4 × 104 mL gcat−1 h−1 with a high TOF of 22.9 h−1, which was much higher than those for the state-of-the-art thermocatalysts, as shown in Table S3,† for example, 3.26 × 102 mL gcat−1 h−1 for Pd/PPS,31 1.2 × 103 mL gcat−1 h−1 for Pd SAs,32 and 3 × 103 mL gcat−1 h−1 for Cu1/ND@G.5 Next, during a 500 min stability test in crude ethylene flow, acetylene conversion was stably above 97.0% (Fig. 4b). These results unambiguously demonstrate the promising applications of SA-Ni-NC for electrocatalytic acetylene semihydrogenation.
To reveal the nature of superior acetylene semihydrogenation performance of SA-Ni-NC, in situ electrochemical Raman (Fig. S24, ESI†) and density functional theory (DFT) investigations were employed. As illustrated in Fig. 4c, in comparison with that of SA-Ni-NC under open circuit potential (OCP), a new Raman peak representing the ν(C
C) stretching vibration appeared at 1961 cm−1 at 0 V under acetylene flow, verifying acetylene adsorption on SA-Ni-NC.33,34 When the potential increased from 0 to −0.6 V, the ν(C
C) stretching vibration of π-bound ethylene gradually emerged at 1546 cm−1.34–37 Weak π-adsorption of ethylene is beneficial for ethylene desorption, eventually facilitating acetylene semi-hydrogenation kinetics. As illustrated by the DFT simulations in Fig. 4d, the adsorbed acetylene molecules were first hydrogenated into *CHCH2 with a lower free energy change of 0.07 eV on SA-Ni-NC than that on Ni NP/N-C (0.37 eV) (the main crystal planes of Ni NP/N-C are (111)), implying a lower activity of Ni NP/N-C than SA-Ni-NC. More importantly, the positive free energy change (0.61 eV) for the *CH2CH3 further hydrogenation step over SA-Ni-NC indicates that deep hydrogenation can be effectively suppressed. Meanwhile, the weak desorption of the produced ethylene is exothermic by 0.23 eV on SA-Ni-NC, implying a rapid ethylene desorption. In addition, energy barriers for water dissociation on both catalysts were calculated to further evaluate the competitive hydrogen evolution reaction (HER) under alkaline conditions. As shown in Fig. S25a,† there is a kinetic energy barrier of 1.02 eV for water dissociation on the Ni (111) surface indicating a sluggish kinetics for water splitting. For Ni SAC shown in Fig. S25b,† the H2O molecule only shows very weak physical adsorption (−0.07 eV), which clearly indicates an even harsher kinetics for water splitting. Therefore, under alkaline conditions, the SA-Ni-NC catalyst exhibits a higher catalytic selectivity (i.e., faradaic efficiency) for acetylene semihydrogenation by suppressing the HER process. These DFT results clearly confirm the efficient occurrence of acetylene semihydrogenation to ethylene on SA-Ni-NC.
| E(vs. RHE) = E(vs. Hg/HgO) + 0.099 V + 0.059 V × pH |
Two-electrode flow cells were also constructed SA-Ni-NC loaded on a gas diffusion electrode as the cathode and NiFe-LDH loaded on nickel foam as the anode for electrocatalytic acetylene semihydrogenation in a 1 M KOH aqueous solution.
485 C mol−1. mcat: catalyst mass in the electrode, g. ω: Ni loading in the catalyst. MNi: atomic mass of Ni, 58.69 g mol−1.
The space velocity (SV) of the catalysts was calculated based on the amount of feed gas according to the following equation:
The SA-Ni-NC catalyst was constructed by embedding the NiN4 moiety in a (5 × 5) supercell of graphene. As for the Ni (111) surface, a (3 × 3) supercell composed of 4 layers of Ni atoms was used, in which the bottom two atomic layers were fixed, while atoms in the other layers were fully relaxed during geometry optimization. For both models, a vacuum layer of around 15 Å was added in the z direction to eliminate the interactions between periodic images.
The Gibbs free energy change (ΔG) of each elementary step during the acetylene hydrogenation process was calculated by using the computational hydrogen electrode (CHE) model proposed by Nørskov et al.43 The chemical potential of the proton–electron pair in aqueous solution is related to one-half of the chemical potential of an isolated hydrogen molecule. Based on this model, the ΔG value can be obtained using the formula ΔG = ΔE + ΔZPE − TΔS, where ΔE is the reaction energy of the reactant and product species adsorbed on the catalyst directly obtained from DFT calculations; ΔZPE and ΔS are the changes between the adsorbed species and the gas phase molecules in zero-point energies and entropy at 298.15 K, which can be calculated from the vibrational frequencies.
Footnotes |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ta08002d |
| ‡ These authors contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2022 |