Zhengyuan
Li
a,
Yanbo
Fang
b,
Jianfang
Zhang
a,
Tianyu
Zhang
a,
Juan D.
Jimenez
c,
Sanjaya D.
Senanayake
c,
Vesselin
Shanov
ab,
Shize
Yang
*d and
Jingjie
Wu
*a
aDepartment of Chemical and Environmental Engineering, University of Cincinnati, Cincinnati, OH 45221, USA. E-mail: jingjie.wu@uc.edu
bDepartment of Mechanical and Materials Engineering, University of Cincinnati, Cincinnati, OH 45221, USA
cChemistry Division, Brookhaven National Laboratory, Upton, NY 11973, USA
dEyring Materials Center, Arizona State University, Tempe, AZ 85287, USA. E-mail: shize.yang@asu.edu
First published on 15th June 2021
The selectivity towards a specific C2+ product, such as ethylene (C2H4), is sensitive to the surface structure of copper (Cu) catalysts in carbon dioxide (CO2) electro-reduction. The fundamental understanding of such sensitivity can guide the development of advanced electrocatalysts, although it remains challenging at the atomic level. Here we demonstrated that planar defects, such as stacking faults, could drive the electrocatalysis of CO2-to-C2H4 conversion with higher selectivity and productivity than Cu(100) facets in the intermediate potential region (−0.50 ∼ −0.65 V vs. RHE). The unique right bipyramidal Cu nanocrystals containing a combination of (100) facets and a set of parallel planar defects delivered 67% faradaic efficiency (FE) for C2H4 and a partial current density of 217 mA cm−2 at −0.63 V vs. RHE. In contrast, Cu nanocubes with exclusive (100) facets exhibited only 46% FE for C2H4 and a partial current density of 87 mA cm−2 at an identical potential. Both ex situ CO temperature-programmed desorption and in situ Raman spectroscopy analysis implied that the stronger *CO adsorption on planar defect sites facilitates CO generation kinetics, which contributes to a higher surface coverage of *CO and in turn an enhanced reaction rate of C–C coupling towards C2+ products, especially C2H4.
It has been demonstrated that the coupling of the key intermediate (e.g., surface adsorbed *CO) is a critical step in CO2-to-C2+ product conversion.5–8 Therefore, numerous strategies have been developed to maximize the utilization of *CO. Cascade catalysis involving tandem catalysts,9,10 tandem electrodes,11,12 and tandem reactors13 has been proposed to supply extra CO and thus to raise local CO concentration near the Cu surface for facilitating C–C coupling. Another widely investigated approach is to properly increase the *CO binding energy via tuning the surface structures of Cu, such as alloying,14 surface doping,15,16 crystal faceting,17,18 subsurface oxygen engineering,19 and surface reconstruction.20 It is worth noting that the CO generation rate can also be enhanced with the *CO binding energy since the activation energy for CO formation decreases according to the scaling relations and Brønsted–Evans–Polanyi (BEP) relations.21 Subsequently, the surface coverage of *CO on Cu can be simultaneously enhanced as *CO binding energy increases.5,22,23
It is well known that the Cu(100) facet with square atomic configurations favors the formation of C2H4, which benefits from the stronger *CO binding energy and subsequently lower *CO dimerization barrier compared to that on the Cu(111) facet.8 Meanwhile, low-coordinated defect sites, such as vacancies, steps, twin boundaries, and grain boundaries, arisen from in situ or ex situ growth on the metallic Cu surface were proposed to be responsible for the enhanced catalytic reactivity.23–26 Indeed, research on twin boundaries and grain boundaries has pointed out that the high activity and selectivity are correlated with the sites binding CO more strongly than low-index Cu facets, and the reactivity is linearly proportional to the density of such active sites.26 But most of these two-dimensional defects are associated with the Cu(111) facets or interfaces between (111) and (100) facets.23,27 A comprehensive comparison between (100) facets and planar defects, which is missing, could assist in designing electrocatalysts for further improvement of CO2-to-C2+ conversion.
Herein, we demonstrated that planar defects could drive the electrocatalysis of CO2-to-C2H4 conversion with higher selectivity and productivity than Cu low-index facets. To this end, we prepared a right bipyramidal Cu (Rbp-Cu) electrocatalyst with both exposed Cu(100) facets and a set of parallel planar defects (e.g., stacking faults). Rbp-Cu achieved a Faradaic efficiency (FE) of 67% and 89% for C2H4 and C2+ products, respectively, at a current density of over 325 mA cm−2 and a corresponding potential of −0.63 V vs. RHE. In contrast, Cu nanocubes (Cube-Cu) with exclusive Cu(100) facets exhibited only 46% FE for C2H4 and a partial current density of 87 mA cm−2 at an identical potential. The planar defect was discovered to show superior activity and selectivity towards C2H4 formation to Cu(100) by comparing the performance between Rbp-Cu and Cube-Cu. The reaction mechanism study using CO temperature-programmed desorption and in situ Raman spectroscopy implied that the enhancement of reactivity for CO2-to-C2H4 conversion originates from the stronger *CO adsorption energy on planar defects than Cu(100).
CO temperature-programmed desorption (CO-TPD) was performed under a 10% CO/He atmosphere. The α-alumina support was applied to disperse 10 wt% Cu NCs by a simple impregnation method. 50 mg of Cu samples were first purged with 10% H2/He at 200 °C for 30 min to remove any impurity on the surface, and then adsorbed CO gas at 50 °C for 60 min. After purging again in pure He at 50 °C to remove physisorbed CO on the sample surface for 60 min, the temperature was raised to 500 °C at a ramp rate of 10 °C min−1. The gas flow rate was fixed at 50 sccm for all steps. The resulting species were traced by on-line mass spectrometry (SRS, RGA 100 Headgas Analyzer). The mass-to-charge (m/z) signal of 28 (CO) was monitored with a time interval of 2 s.
For each applied potential, an on-line gas chromatograph (GC, SRI Instruments MultipleGas#5) equipped with both a thermal conductivity detector and a flame ionization detector was used to monitor the gas products. To calibrate the outlet gas flow rate of CO2, Ar, as an internal standard, was fed at 10 sccm and mixed with the outlet gas stream from the flow cell before looping to the GC.11 The FEs for gas products were calculated using the following equation (eqn (1)):
(1) |
Meanwhile, the catholyte was collected after electrolysis followed by quantifying the liquid products via1H NMR (Bruker AV 400 MHz spectrometer). 500 μl of the catholyte was mixed with a 100 μl internal standard of 5 mM 3-(trimethylsilyl)propionic-2,2,3,3-d4 acid sodium salt in D2O. The standard deviations were calculated based on the measurements of three independent electrodes.
The total CO generation and dimerization rates were calculated according to the following two equations (eqn (2) and (3)):
(2) |
(3) |
We used Pb underpotential deposition (UPD), a surface-sensitive electrochemical technique, to further probe the surface structures of Cu NCs. Different Cu facets exhibit distinguished Pb-UPD peaks in the CV curves.17,33 As displayed in Fig. 1d, the peaks located at ca. −0.10 V and −0.08 V can be assigned to Cu(100) and Cu(111), respectively.17 The peaks of Pb electrodeposition on Cu(100) are more prominent for both Rbp-Cu and Cube-Cu samples than those on Cu(111). An additional peak centered at a less cathodic potential (ca. −0.06 V) is only observed for Rbp-Cu, which is ascribed to the Pb-UPD on the low-coordinated sites from surface defects.33,34 Such a kind of phenomenon that the Pb over-layer on the face-centered cubic metal preferentially nucleates at the edges or steps and then expands to the terraces has been demonstrated by in situ scanning probe microscopy studies.35 OH− electro-adsorption was also performed to verify the predominant (100) facets for both Cu NCs (Fig. S4†). An additional shoulder peak at ca. 0.34 V was observed for Rbp-Cu, which is associated with surface defects.36 The percent of high-index surface defects on Rbp-Cu is estimated to be 17.3% according to the integral peak areas,37 which is close to the value obtained from TEM aforementioned. In agreement with TEM, Pb-UPD and OH− electro-adsorption results confirm the presence of defect sites on Rbp-Cu, although the exact facet index is hard to identify. XRD patterns show that both types of Cu NCs exhibit a much more pronounced (100) peak and a weaker (111) peak (Fig. 1e), which is consistent with electrochemical characterization results and confirms the uniform shape of the NCs. Moreover, other than Cu peaks, no other peaks associated with copper oxide or other impurities were observed. XPS analysis reveals that Rbp-Cu and Cube-Cu both present two main peaks at 932.8 eV and 952.6 eV in the Cu 2p spectrum (Fig. 1f), corresponding to Cu0 or Cu+ species.38 The Auger Cu LMM spectrum (Fig. S5†) shows that both Rbp-Cu and Cube-Cu mainly consist of Cu0 (918.6 eV) at the surface, with a minimal amount of Cu+ (916.6 eV). The formation of Cu+ is caused by the rapid oxidation of Cu0 in air and would be quickly in situ reduced to metallic Cu under CO2 electrolysis.22,39 Due to the introduction of a trace amount of Pd (<0.55 at%) as the seed during the synthesis of Rbp-Cu, Pd 3d XPS was carried out. There is no Pd element on the Rbp-Cu surface (Fig. S6†), suggesting the absence of Pd dopants or Pd–Cu bimetallic structures.
The CO2RR performance on Rbp-Cu and Cube-Cu based GDEs was measured in a customized flow cell using 1 M KOH as the catholyte (Fig. S7 and S8†). Rbp-Cu exhibited a higher selectivity towards C2H4 than Cube-Cu. The maximum FE of C2H4 for Rbp-Cu reached 67% at −0.63 V (vs. RHE, thereafter), at which the partial current density (jC2H4) reached 217 mA cm−2 (Fig. 2a and b). The overall performance of Rbp-Cu for CO2-to-C2H4 conversion is superior to that of most copper-based catalysts tested under identical conditions (Table S1†). Combined with C2+ liquid products analyzed by 1H NMR, including ethanol (C2H5OH), acetate (CH3COO−), and n-propone (C3H7OH), the FE and partial current density of C2+ products (jC2+) on Rbp-Cu reached 89% and 289 mA cm−2, respectively, at −0.63 V (Fig. 2c and d). In contrast, Cube-Cu showed only 46% FE for C2H4 and 68% FE for C2+ products at the same overpotential. Moreover, the FE ratio of C2H4 over C2H5OH was enhanced from 3.10 for Cube-Cu to 3.77 for Rbp-Cu at −0.63 V (Fig. S9†). Compared to Rbp-Cu, Cube-Cu manifested lower jC2H4 and jC2+ over the applied potential range (Fig. 2b and d). Interestingly, we found that Rbp-Cu showed a relatively much lower FE for CO than Cube-Cu (Fig. S10†), indicating the more efficient conversion of the *CO intermediate to hydrocarbons and oxygenates on Rbp-Cu.10 The improvement of activity and selectivity towards C2H4 and C2+ products suggests that the planar defect sites on Rbp-Cu may be responsible for the enhanced C–C coupling rate.
To exclude the possible effect of surface roughness, the electrochemical surface area (ECSA) was determined by using the double-layer capacitance (Cdl) method (Fig. S11†). The ECSA normalized current density of each product followed the same trend as the geometric current density for both Rbp-Cu and Cube-Cu (Fig. S12†). At an optimal potential of −0.63 V for C2H4 and C2+ product formation, the normalized jC2H4 and jC2+ for Rbp-Cu were 1.9-fold and 1.7-fold higher than those for Cube-Cu, respectively, indicating the excellent intrinsic activity of Rbp-Cu towards C2H4 and C2+ products. Considering a lower surface-to-volume ratio of Cube-Cu related to Rbp-Cu,28 the loading of Cube-Cu was increased from 0.30 to 0.42 mg cm−2, which increased the Cdl (proportional to the ECSA) to 0.57 mF cm−2 accordingly, similar to that of Rbp-Cu (0.56 mF cm−2) with a loading of 0.30 mg cm−2 (Fig. S13†). Although the current density increased with the catalyst loading, a similar product distribution was obtained for Cube-Cu with two different loadings. Especially, the FEs for C2H4 were almost identical in the investigated overpotential range (Fig. S13e†). Therefore, we infer that a higher ECSA is not the main factor for promoting selectivity towards C2H4.
Notably, Rbp-Cu exhibited activity towards C2H4 and C2+ products similar to Cube-Cu at potentials between −0.35 and −0.50 V (Fig. 2 and S12†), likely owing to the Cu(100) facet on both Cu NCs which facilitates the *COCO dimerization in such a low overpotential region.7,39 In the potential region of −0.50 ∼ −0.65 V, the planar defect sites on the Rbp-Cu surface play a more critical role in CO2-to-C2+ product conversion. The possible dimerization pathways in this overpotential range may involve coupling between *CO and reduced *CO species (i.e., *CHO and/or *COH),8,40 the reaction barrier of which could be decreased on low-coordinated defects.37,41,42 However, the selectivity towards C2H4 for Rbp-Cu decreased dramatically at higher cathodic potentials (−0.65 ∼ −0.70 V). Simultaneously, the activity and selectivity towards hydrogen (H2) and methane (CH4) rose accordingly (Fig. S10 and S12†). The upsurge of reactivity towards H2 arises from the stronger competitive hydrogen evolution reaction (HER) occurring on low-coordinated sites in the high cathodic potential region.30,43 To verify this, the HER was performed under an Ar atmosphere in a flow cell. Linear sweep voltammetry (LSV) showed higher electrocatalytic activity of Rbp-Cu for the HER than Cube-Cu as the potential was swept beyond −0.50 V (Fig. S14†), demonstrating that planar defects are more active for the HER as well as the CO2RR in the high cathodic potential region. At higher overpotentials, increased adsorption of hydrogen would suppress CO2-to-C2+ product conversion.44
To demonstrate a stable CO2-to-C2+ product conversion, a durability test for Rbp-Cu was conducted at −0.63 V, where optimal selectivity towards C2H4 and C2+ products was achieved. By using polytetrafluoroethylene-treated GDL to improve the hydrophobicity, a relatively stable FE of (61 ± 5)% for C2H4 was recorded while FE for C2+ products was maintained at (82 ± 4)% with a total current density of (330 ± 30) mA cm−2 during a 50 hour continuous test (Fig. 2e). As examined by TEM, the post-electrolysis Rbp-Cu retained the planar defects (Fig. S15†), although the corners or edges of the well-defined morphology were degraded to some extent.45,46 Surface adsorbed *CO, the key intermediate for CO2 reduction to hydrocarbons and oxygenates, may stabilize Cu defects in line with prior density functional theory (DFT) calculations.22
Given that the same Cu(100) facets were exposed for both Rbp-Cu and Cube-Cu samples, we envisioned that the superior CO2-to-C2H4 conversion is attributed to the introduction of the planar defects on the Rbp-Cu surface. These planar defects with a low coordination number enable stronger *CO adsorption that assists in building a sufficiently high surface coverage of *CO to direct the C–C coupling.6,23 To prove this hypothesis, an ex situ CO-TPD experiment and in situ electrochemical Raman spectroscopy measurements were conducted.
The monitored CO-TPD traces, which can give the amount and strength of *CO adsorption, are shown in Fig. 3a and S16.† For both types of Cu NCs, the onset CO desorption began at approximately 60 °C. The observation of a later CO desorption peak at 90 °C for Rbp-Cu indicates a stronger adsorption affinity for CO and CO-related intermediates, resulting from the presence of planar defects.36,47 Additionally, the larger CO desorption peak area for Rbp-Cu implies a larger surface coverage of *CO on Rbp-Cu than Cube-Cu.48 These CO-TPD results preliminarily support that both the adsorption energy and surface coverage of *CO are boosted on Rbp-Cu, benefiting from planar defect sites.
To evaluate the *CO adsorption behaviors under CO2 electrocatalysis conditions, in situ Raman spectroscopy measurements were further carried out in a modified flow cell with a quartz window in front of the cathode GDE (Fig. 3b). The Raman spectra were recorded in a range of applied potentials from the open circuit potential (OCP) to −0.70 V for both types of Cu NCs (Fig. 3c and d). For Rbp-Cu, the low-frequency bands of *CO at ca. 282 and 371 cm−1 are characteristic of Cu–CO frustrated rotation and Cu–CO stretch, respectively (Fig. 3c).49 The high-frequency bands are associated with CO stretching vibrations with different adsorption configurations, including bridge-bound CO (CObridge, 1850–1880 cm−1) and atop-bound CO (COatop, 2000–2090 cm−1).49,50 Notably, the asymmetric broad band for COatop is likely due to the *CO adsorption on various Cu surface sites, which, in turn, may affect the activity and selectivity of the CO2RR.49,51 The difference in *CO adsorption behaviors between Rbp-Cu and Cube-Cu was assessed by comparing the in situ Raman spectra collected at the same applied potential of −0.63 V. First, a blue shift of the COatop sharp band suggests the stronger binding of *CO at low-coordinated defect sites for Rbp-Cu (2084 cm−1) relative to Cube-Cu (2067 cm−1) (Fig. 3f).49 Meanwhile, the band for the Cu–CO stretch also exhibited a blue shift on Rbp-Cu (371 cm−1) compared to that on Cube-Cu (365 cm−1) (Fig. 3e), confirming a stronger binding of *CO to the Rbp-Cu surface.52 Second, the integrated areas of bands (e.g., Cu–CO and COatop bands), which are proportional to the *CO surface coverage,49 increase for Rbp-Cu compared with Cube-Cu (Fig. 3e and f). Taken together, ex situ CO-TPD and the in situ Raman spectroscopy results unravel that Rbp-Cu promotes the adsorption and surface coverage of *CO, which could enhance the subsequent C–C coupling kinetics and thus the production of C2+ products.
To further explore the impact of the enhancement of *CO adsorption energy on the C–C coupling kinetics, we analyzed and compared the CO generation and dimerization rates (Fig. 4 and S17†).12 The CO generation rate was referred to as the summary of the normalized production rates of CO, CH4, and C2+ products (eqn (2) in the Experimental section). According to the scaling relations of adsorption energy between intermediates and BEP relations between the intermediate adsorption energy and activation barrier,21 the enhancement of *CO binding indicates stronger *COOH binding and a lower activation barrier in the elementary step of CO2-to-*COOH conversion, which leads to an accelerated CO generation rate. Fig. 4a shows that Rbp-Cu indeed provided more CO than Cube-Cu, manifesting that planar defects facilitate the CO generation rate compared to terraces such as (100) facets in our case. Meanwhile, the CO dimerization rate was derived based on the normalized production rates of C2+ products (eqn (3) in the Experimental section). Again, Rbp-Cu exhibited a much faster CO dimerization rate than Cube-Cu due partly to the increased surface coverage of *CO (Fig. 4b). For example, the CO dimerization rate of Rbp-Cu reached >2.2-times as high as that of Cube-Cu at around −0.65 V. It is also reasonable to propose that C–C coupling intermediates, such as *COCO, *COCHO, or *COCOH, can be efficiently stabilized on defect sites, which further facilitates the dimerization reaction.40,41 In short, compared to (100) facets, the stronger *CO adsorption on planar defect sites facilitates CO generation kinetics in the intermediate overpotential region, which contributes to a higher surface coverage of *CO and in turn an enhanced reaction rate of C–C coupling towards C2+ products, especially C2H4.
Fig. 4 Comparison of the (a) CO generation rate and (b) CO dimerization rate between Rbp-Cu and Cube-Cu at different applied potentials. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ta02565a |
This journal is © The Royal Society of Chemistry 2021 |