Electrochemical conversion of CO 2 to carbon-based products, which can be used directly as fuels or indirectly as fuel precursors, is suggested as one of the promising solutions for sustainability. Not only does this process allow using renewables such as solar electricity as energy input, but CO 2 emitted from the consumption process can be recycled back into fuels. The success of this technology depends on the value added to the product that forms from CO 2 , and therefore it is important to facilitate multicarbon product generation. This work presents a copper-based catalyst, formed in situ from an ensemble of nanoparticles, that is able to selectively generate C 2 –C 3 products at low overpotentials with good stability, where their efficient formation has been difficult to achieve.

Direct conversion of carbon dioxide to multicarbon products remains as a grand challenge in electrochemical CO 2 reduction. Various forms of oxidized copper have been demonstrated as electrocatalysts that still require large overpotentials. Here, we show that an ensemble of Cu nanoparticles (NPs) enables selective formation of C 2 –C 3 products at low overpotentials. Densely packed Cu NP ensembles underwent structural transformation during electrolysis into electrocatalytically active cube-like particles intermixed with smaller nanoparticles. Ethylene, ethanol, and n-propanol are the major C 2 –C 3 products with onset potential at −0.53 V (vs. reversible hydrogen electrode, RHE) and C 2 –C 3 faradaic efficiency (FE) reaching 50% at only −0.75 V. Thus, the catalyst exhibits selective generation of C 2 –C 3 hydrocarbons and oxygenates at considerably lowered overpotentials in neutral pH aqueous media. In addition, this approach suggests new opportunities in realizing multicarbon product formation from CO 2 , where the majority of efforts has been to use oxidized copper-based materials. Robust catalytic performance is demonstrated by 10 h of stable operation with C 2 –C 3 current density 10 mA/cm 2 (at −0.75 V), rendering it attractive for solar-to-fuel applications. Tafel analysis suggests reductive CO coupling as a rate determining step for C 2 products, while n-propanol (C 3 ) production seems to have a discrete pathway.

With rising concerns about the anthropogenic impacts of current trends in energy use, as well as the prospect of continuing these trends to meet future needs (1), we are at a stage where revolutionary change to our energy paradigm is a must. Various methods for effectively using solar energy are being developed to power and support the global population (2⇓–4). Among them, artificial photosynthesis is considered vital to meeting our goal toward long-term global sustainability (5). The successful introduction of artificial photosynthesis technology will highly depend on the development of every functional component essential to the efficient operation of the overall system.

As energy from sunlight eventually ends up in chemical bonds by the photocatalytic or electrocatalytic component, development of an effective catalytic material to facilitate the conversion process becomes important. Over the past several decades, the focus has been on using water as the starting substrate and converting it to hydrogen gas (6). More recently, carbon dioxide has been considered as a promising substrate, and many efforts have been underway to find efficient electrocatalysts that can selectively operate for reducing CO 2 in aqueous solutions against the competing hydrogen evolution (7⇓⇓⇓⇓⇓⇓⇓⇓–16). However, major progress has been limited to two-electron reduced products of CO and formate. Still, the formation of multicarbon products involving multiple proton and electron transfers remains as one of the biggest scientific challenges to be addressed.

Starting from the idea that element copper is a key component to forming multicarbon products (17, 18), there have been various studies so far where formation of products such as C 2 H 4 , C 2 H 6 , and C 2 H 5 OH has been observed often with the requirement of large overpotentials (potential applied ≤−1 V vs. RHE) (19⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–35). These methods mostly rely on reducing certain forms of oxidized copper (19⇓⇓⇓⇓⇓–25, 27⇓⇓–30, 32⇓–34) (either oxides or halides) and even this approach has been extended to reduce carbon monoxide instead (36), a common intermediate for CO 2 reduction, to circumvent difficulties associated with C–C coupling starting from CO 2 . Furthermore, to instead create a favorable environment for multicarbon product formation, there have been attempts to use gas-diffusion electrodes with alkaline electrolytes (37). It would certainly be desirable to discover an electrocatalyst that can directly reduce CO 2 to multicarbon products with high selectivity and energy efficiency (i.e., minimal energy loss from low overpotentials).

Here, we show that an ensemble of densely packed copper nanoparticles (NPs) could enable selective conversion of CO 2 to multicarbon products, while significantly suppressing C 1 formation. Catalytically active cube-like structures, capable of forming ethylene, ethanol, and n-propanol, are formed during electrolysis by the structural transformation of the Cu NP ensemble. These structures can selectively generate C 2 and C 3 products together at low overpotentials in neutral pH aqueous media, illustrating the importance of in situ structural evolution in CO 2 electrocatalysis. We also find that the catalyst support plays an important role for high multicarbon selectivity. This work suggests an alternative route to development of catalysts for multicarbon products and understanding of their formation, where the field has been heavily reliant on using oxidized copper as starting materials.

Results and Discussion

Monodisperse Cu NPs (size 6.7 nm) prepared (SI Appendix, Fig. S1) were directly deposited onto carbon paper support (1 cm2 geo ) to form densely packed NP ensembles. Cu NP loading was systematically increased (SI Appendix, Table S1) starting from the lowest loading of ∼2 μg of Cu (×1). Number density of NPs was determined based on the estimated surface area of the carbon paper support (24, 26) (SI Appendix, Fig. S2), which was at ∼5.9 cm2 real /cm2 geo (roughness factor ∼5.9). Most of the NPs are isolated at the lowest loading condition, and increased loadings resulted in densely packed arrangements of Cu NPs (Fig. 1A). In the case of ×22.5 loading, the surface was mostly covered with closely packed Cu NPs (SI Appendix, Fig. S3).

Fig. 1. Varied density of Cu NP ensembles and their electrocatalytic activity. (A) SEM images of Cu NPs loaded on carbon-paper support at ×1 loading, ×5 loading, and ×22.5 loading. (B) FEs (%) for C 1 , C 2 , and C 3 products. (C) Relative ratio of the FEs. (D) Ethylene, ethanol, and n-propanol FE with the dotted line showing the overall C 2 –C 3 FE. Activity measured at −0.81 V vs. RHE, using 0.1 M KHCO 3 saturated under 1 atm CO 2 . Error bars shown in B–D are 1 SD from three independent measurements.

Cu NP ensembles with varied loading densities were tested for their electrocatalytic CO 2 reduction activity, under identical conditions of 0.1 M KHCO 3 at 1 atm CO 2 . From product analysis (SI Appendix, Fig. S4), we found that increased loadings resulted in a drastic rise of the C 2 –C 3 faradaic efficiency (FE) (Fig. 1B and SI Appendix, Fig. S5 and Table S2). This trend is consistent with the observed loss of C 1 products, indicating that carbon-based intermediates could be effectively coupled to yield multicarbon products. When plotting the relative ratio of the C 2 –C 3 FE to C 1 FE (Fig. 1C), charge consumed to reduce CO 2 was mainly from the reaction pathways to C 2 –C 3 products at increased loading conditions, reaching up to 76% out of the total CO 2 reduction products. Similar trends can be seen with CO 2 reduction dominating over the competing H 2 evolution (Fig. 1C). Among various C 2 –C 3 products, ethylene (C 2 H 4 ), ethanol (EtOH), and n-propanol (n-PrOH) were the majority, constituting 94% out of the total C 2 –C 3 products generated (Fig. 1D).

When probing the product distribution over time for the ×22.5 loading condition, an abrupt change occurred during the initial period (Fig. 2B). Hydrogen was the dominant product when gas products were measured 3 min after the start of electrolysis. Selectivity for C 2 H 4 increased afterward. A similar trend was found for the liquid products as well (SI Appendix, Fig. S6), where liquids analyzed for the first 7 min had less multicarbon products relative to formate. Visual inspection of the electrode (SI Appendix, Fig. S7) also supported the fact that product distribution was shifting during its initial electrochemical testing, as more gas bubbles were observed at the beginning of electrolysis, probably due to the majority of charge being consumed for two-electron transfer products, such as hydrogen.

Fig. 2. Structural transformation of Cu NP ensembles. (A) Schematic illustrating the transformation process of Cu NP ensembles to an active catalyst for C 2 –C 3 product formation. (B) Total current density (based on geometric area) versus time plot for ×22.5 loading condition at −0.81 V vs. RHE. FE of gas products are shown at the time point of measurement. FE of CH 4 and C 2 H 6 are omitted because of their low values. (C) SEM image of ×22.5 loaded carbon-paper support electrode after 7 min of electrolysis at −0.81 V vs. RHE. (D) Investigation of the parameters affecting structural transformation of Cu NP ensembles and their catalytic activity. Three different conditions have been tested: (i) separation of the NPs from their initial densely packed assembly, (ii) use of Cu nanocubes as starting materials, and (iii) change of support to a low surface area carbon plate. FE of C 1 , C 2 , and C 3 products obtained from trans-CuEn (left column, shaded, at −0.81 V vs. RHE) are compared with the activity measured for three different conditions [at (i) −0.84 V, (ii) −0.86 V, and (iii) −0.81 V vs. RHE, respectively]. Electrochemical tests were conducted using 0.1 M KHCO 3 solution at 1 atm CO 2 .

This observation indicated that the NP ensemble may go through a structural transformation process during initial electrolysis. Instead of the starting densely packed Cu NPs (×22.5 loading, Fig. 1A), cube-like particles (10 ∼ 40 nm) mixed together with smaller NPs were observed on the carbon support after electrolysis (SI Appendix, Fig. S8). Carbon-paper supports with lower loading were also checked after electrolysis and a trend could be observed: the more densely packed the Cu NPs, the more likely the formation of cube-like particles (SI Appendix, Fig. S9). When Cu NPs were sparsely covering the support, random aggregates of NPs together with the pristine NPs could be found after electrolysis. Surface-area analysis of Cu NP ensembles after electrolysis also showed that the densely packed NPs transformed to larger particles (SI Appendix, Fig. S10). As expected, loss of surface ligands during electrolysis and structural transformation was confirmed by elemental analysis (SI Appendix, Figs. S11 and S12). Structural transformation of the NP ensemble (Fig. 2A) occurred during the initial stage of electrolysis. This was confirmed from observation of the electrode 7 min after the start of electrolysis (Fig. 2C), which coincided with the shift in catalytic activity (Fig. 2B). Negligible catalytic activity of the underlying carbon paper (SI Appendix, Fig. S13) further supports that the structure derived from Cu NP ensembles is responsible for enhanced C 2 –C 3 formation. This catalytically active structure formed starting from densely packed Cu NP ensembles (×22.5 loading on carbon-paper support in 0.1 M KHCO 3 ), hereafter referred to as trans-CuEn, was further investigated.

As the initial loading density of Cu NP ensembles (and their densely “packed-ness”) tends to govern their structural transformation during electrolysis and resulting electrochemical activity, we tried to intentionally separate the Cu NPs in the precursor state to trans-CuEn. We expected the transformation process to cube-like structures to be disrupted, leading to diminished C 2 –C 3 selectivity. Cu NPs (×22.5 loading) were mixed with carbon black before depositing on carbon-paper support (SI Appendix, Fig. S14), which led to NPs being spatially separated from each other. Under this condition, substantial loss of C 2 –C 3 product selectivity (FE from 49 to 17%) was observed (Fig. 2D), while CO and HCOO− became major products. When particles were examined after electrolysis, the structure more resembled what would be observed for low-density conditions (SI Appendix, Figs. S9 and S14). Cu NPs have been observed to electrically fuse into irregularly shaped large crystals under strong bias conditions (<−1.25 V vs. RHE), reaching a similar state irrespective of the initial conditions (38). Here, we find that structural transformation can be caused not only under low bias conditions, but controlled by the initial arrangement of NPs, and consequently catalytic behavior for multicarbon products can be significantly improved.

As trans-CuEn displays cubic-shaped particles, copper nanocubes loaded onto carbon-paper support were tested under identical conditions for comparison. We used Cu nanocubes that have been previously studied for CO 2 reduction (26) (SI Appendix, Fig. S15). Specifically, cubes with edge length 25 nm were used (with copper loading mass identical to trans-CuEn) to approximately match the cubic-shaped particles that vary in size (10 ∼ 40 nm) for trans-CuEn. In contrast to trans-CuEn, observed structural changes were minimal where the cubes seem to have sintered or roughened (SI Appendix, Fig. S16). Furthermore, only small amounts of multicarbon products were detected (Fig. 2D). The result is consistent with the earlier report of Cu nanocubes, which claims multicarbon product formation at high overpotentials (<−1 V vs. RHE) (26). Therefore, we find that simple reproduction of the key morphological feature present in trans-CuEn is insufficient to reach high multicarbon selectivity.

This leads to the possibility of cube-like particles derived in situ during electrolysis featuring unique active sites for C 2 –C 3 formation. Recently, scanning tunneling microscopy investigation of copper for carbon monoxide reduction has shown not only the reconstruction of a polycrystalline surface to a (100) surface, but also the additional structural transformation unique to the (100) reconstructed copper, leading to stepped surfaces which selectively generate ethanol (39, 40). While this observation may have been for reducing CO, together with the findings here, it brings into attention the importance of in situ structural transformation for multicarbon product formation in copper-based catalysts. In addition, we would like to point out that while the vast majority of research has been to use oxide-derived structures, with even some reports claiming the importance of remaining oxidized copper (29, 41), the catalytically active structure derived here is from pristine Cu NPs (with a thin layer of surface oxide naturally present). Furthermore, we find that the structural transformation observed is unique to the original Cu NPs (SI Appendix, Fig. S17). Therefore, it would be important to understand how this structural transformation proceeds and what type of active site motifs are present under working conditions. This is especially the case for copper, which may oxidize after electrolysis, possibly leading to the loss of surface atomic information (SI Appendix, Fig. S12). However, we also cannot rule out the possibility that high multicarbon selectivity stems from having a mixture of particles (42), which are the cube-like ones together with smaller particles. With all of these taken into consideration, further investigation into the structural origin of high multicarbon selectivity from Cu NP ensembles is needed.

Furthermore, we investigated the role of the catalyst support by depositing Cu NPs onto a highly polished graphite plate (1 cm2 real , roughness factor ∼1), while keeping the NP density (/cm2 real ) identical to that of trans-CuEn. Structural transformation occurred in a similar way resulting in cubic-shaped particles (SI Appendix, Fig. S18). However, H 2 and C 1 products were the major products (Fig. 2D and SI Appendix, Fig. S18). We speculate that this difference is due to local pH effects discussed in earlier reports (32, 43⇓⇓–46), as the loss in the real surface area of the underlying support led to a sharp decrease of the geometric current density (lowered to ≈1/5 of the original). The increased local pH by large current density of trans-CuEn (on carbon-paper support) could play a role in determining its catalytic behavior. Therefore, it seems that it is important to not only start from a high density of closely packed Cu NPs to facilitate the structural transformation, but also have the underlying support provide sufficient surface area. This shows why high C 2 –C 3 selectivity was not observed from the previous report of Cu NP monolayers (16).

Catalytic activity of trans-CuEn was further probed at various potentials (Fig. 3 and SI Appendix, Fig. S19) in 0.1 M KHCO 3 . The onset of C 2 –C 3 formation was observed at only −0.57 V vs. RHE, with products mainly comprising C 2 H 4 , EtOH, and n-PrOH. Compared with that of the pristine copper foil (18), overpotentials were lowered by 180 mV for C 2 H 4 and 390 mV for EtOH and n-PrOH, respectively (SI Appendix, Table S3). Beyond this potential, a substantial rise in C 2 –C 3 FE was observed (SI Appendix, Table S4), with the highest selectivity toward C 2 –C 3 products (55%) achieved at −0.86 V vs. RHE. The high selectivity for C 2 –C 3 products, including oxygenates, is quite significant, compared with previously reported catalysts for C 2 –C 3 product formation around similar overpotentials applied in neutral pH aqueous media (SI Appendix, Table S5). So far, catalysts for multicarbon products have been Cu-based (mostly derived from oxides or halides) and require bias applied close to and beyond −1 V vs. RHE (SI Appendix, Table S6), where even only some of them reach product distributions dominated by C 2 –C 3 products (C 2 –C 3 > C 1 + H 2 ). Furthermore, with major efforts in the field toward using oxidized Cu as a starting template, the discovery of this catalyst presents an approach to achieving high C 2 –C 3 selectivity for electrochemical CO 2 reduction. In contrast, FE for two-electron reduced products (CO and formate) could be lowered to ∼1%, implying that almost all of the CO 2 interacting with the catalyst could undergo C–C coupling to yield more complex products (Fig. 3B). In assessing catalytic performance for multicarbon product formation, earlier reports have been using C 2 H 4 /CH 4 FE ratio as a figure of merit and trans-CuEn exhibits significantly high values at low overpotentials (∼252 at −0.78 V vs. RHE) that are comparable or better than previous catalysts reported for selective formation of C 2 H 4 (SI Appendix, Fig. S20 and Table S7). More negative bias applied leads to increase in CH 4 formation and C 1 FE.

Fig. 3. Electrochemical CO 2 reduction activity of trans-CuEn. (A) FE of C 1 , C 2 , and C 3 products at various potentials for trans-CuEn. (B) FE of major products at various potentials for trans-CuEn. Electrochemical tests were conducted using 0.1 M KHCO 3 solution at 1 atm CO 2 . Error bars shown are 1 SD from three independent measurements.

It has been suggested that larger cations promote higher concentrations of CO 2 near the catalyst surface, leading to increased activity (17, 47). For further optimization, Cu NP ensembles were tested in 0.1 M CsHCO 3 aqueous electrolyte saturated with 1 atm CO 2 and a similar trend was observed where increased loading densities resulted in higher C 2 –C 3 selectivity (SI Appendix, Fig. S21). Transformation of Cu NP ensembles (at optimized condition of ×32.5 loading in 0.1 M CsHCO 3 ) consistently resulted in cube-shaped particles mixed together with smaller NPs (SI Appendix, Fig. S22), hereafter referred to as trans-CuEn 2. Activity of trans-CuEn 2 was measured at various potentials (Fig. 4 and SI Appendix, Fig. S23) and high C 2 –C 3 selectivity was observed at more positive potentials with the onset of C 2 –C 3 formation at only −0.53 V vs. RHE (SI Appendix, Table S3), which is 40 mV less of applied overpotential compared with that observed in 0.1 M KHCO 3 . Highest C 2 –C 3 selectivity (∼50%) was observed at −0.75 V vs. RHE, shifting the potential 110 mV more positive relative to the point of maximum C 2 –C 3 FE in 0.1 M KHCO 3 . Therefore, with this catalytic structure, selective electrocatalytic conversion CO 2 to C 2 –C 3 hydrocarbons and oxygenates could be achieved at significantly reduced overpotentials, compared with what have been demonstrated up to now (SI Appendix, Tables S5 and S6). Similarly, the main products were C 2 H 4 , EtOH, and n-PrOH (Fig. 4B and SI Appendix, Table S8) constituting up to 95% of total C 2 –C 3 products (SI Appendix, Fig. S24). In addition, not only were FEs of CO and formate decreased to very low levels (1∼2%), but CH 4 formation could also be suppressed (<1%) across the entire potential region, resulting in a significantly high C 2 H 4 /CH 4 ratio (∼2,133 at −0.73 V vs. RHE) at low overpotentials (SI Appendix, Fig. S24 and Table S7). Owing to its high C 2 –C 3 selectivity in 0.1 M CsHCO 3 , the proportion of C 2 –C 3 products among the total CO 2 reduced products reached up to 90% (Fig. 4C).

Fig. 4. Electrocatalytic behavior of trans-CuEn 2 (×32.5 loading in 0.1 M CsHCO 3 ). (A) FE of C 1 , C 2 , and C 3 products at various potentials. (B) FE of major products at various potentials. (C) Relative ratio of the FE. (D) Logarithmic specific current density (corrected by the real surface area of the catalyst) plots for C 2 H 4 , EtOH, and n-PrOH. Electrochemical tests were conducted in 0.1 M CsHCO 3 solution at 1 atm CO 2 . Error bars shown are 1 SD from three independent measurements. (E) Long-term electrolysis at −0.75 V vs. RHE with gas products measured every hour. Column graph on the right shows FE of EtOH and n-PrOH measured after electrolysis and C 1 , C 2 , and C 3 product FEs for the overall run.

With the real surface area of trans-CuEn 2 measured (SI Appendix, Fig. S25), specific current density plots and Tafel slopes of the three major products could be obtained (Fig. 4D). Both C 2 H 4 and EtOH exhibit similar slopes (∼120 mV/dec), indicative of a rate-determining step with a common intermediate. Furthermore, C 2 H 4 and EtOH start forming in the potential region where CO evolution is dominant and increase while CO diminishes (Figs. 3B and 4B), suggesting that formation of these C 2 products is essentially limited by the coupling of major C 1 intermediates. It has been also shown that higher coverages of *CO can be expected in the region where CO formation is majorly observed (48). Therefore, with a slope close to 120 mV/dec suggesting a single electron transfer step, we expect the rate-determining step to be a reductive coupling (i.e., dimerization) step of adsorbed CO intermediates, predicted from theory and carbon monoxide reduction experiments on copper (43, 49⇓–51): * CO + * CO + e − → * C 2 O 2 − . On the other hand, n-PrOH exhibits a different slope, suggesting a distinct rate-determining step from that of C 2 products. Estimated value is rather close to that observed for CH 4 on copper foil (86 mV/dec) (38). In addition, it has been reported that n-PrOH formation only occurs when reactants include both CO (carbon monoxide) and C 2 H 4 , while CO reduction solely leads to EtOH (28). If C 3 products followed the same pathway as C 2 products, n-PrOH should have been observed upon CO reduction. Instead, it may be that n-PrOH formation requires coupling between CO and hydrogenated carbon [e.g., carbene (*CH 2 )], which is a major intermediate in the pathway to CH 4 (50). CH 4 formation activity of trans-CuEn and -CuEn 2 supports this idea as well (Figs. 3B and 4B and SI Appendix, Tables S4 and S8). In contrast to C 2 H 4 and EtOH, n-PrOH reaches peak selectivity at a more positive potential and the potential in which n-PrOH FE drops coincides well with the point where CH 4 FE starts to rise. However, it is still unclear how formation of C 3 products occur and an in-depth study of the mechanistic pathways to these products is needed.

Long-term stability was demonstrated by 10 h electrolysis of trans-CuEn 2 at −0.75 V vs. RHE (Fig. 4E). Average C 2 –C 3 FE reached ∼50% for the overall run and structural features of trans-CuEn 2 were maintained as well (SI Appendix, Fig. S26). Furthermore, stable C 2 –C 3 product current density of 10 mA/cm2 geo was achieved, which is potentially attractive for solar-to-fuel applications. As long-term electrolysis accumulates significant amounts of liquid products, propionaldehyde, likely to be the precursor to n-PrOH, was detected (SI Appendix, Fig. S27).