Asato
Inoue
a,
Takashi
Harada
ab,
Shuji
Nakanishi
*ab and
Kazuhide
Kamiya
*ab
aResearch Center for Solar Energy Chemistry, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan. E-mail: nakanishi.shuji.es@osaka-u.ac.jp; kamiya.kazuhide.es@osaka-u.ac.jp
bInnovative Catalysis Science Division, Institute for Open and Transdisciplinary Research Initiatives (ICS-OTRI), Osaka University, Suita, Osaka 565-0871, Japan
First published on 3rd November 2022
CO2 electrolysis to value-added products is a promising technology to close the carbon cycle and sequester anthropogenic CO2 into chemical feedstocks; an increase of the current density for multicarbon products is one of the requirements for practical implementation. We have successfully increased the partial current density for gaseous CO2 reduction reactions to multicarbon products (C2+) over Cu nanoparticles on gas diffusion electrodes in neutral electrolytes to a record value of 1.7 A cm−2. The faradaic efficiency for multicarbon products increased with the current density below total current density of 2000 mA cm−2 and reached 76% at a total current density of 1600 mA cm−2. The turnover frequency for the production of C2+ per Cu atoms exceeded 1.1 s−1. Optimizing the standard components and their assembly as the cathode elicits the high-turnover frequency of oxide-derived Cu catalysts, resulting in the record partial current density for C2+. Especially, we demonstrated that the thickness of catalyst layers was one highly sensitive factor in determining the maximum current density for C2+.
Broader contextThe electrochemical reduction of carbon dioxide (CO2) has attracted considerable attention to maintain the carbon cycle. An increase in the current density for value-added multicarbon products (C2+), such as ethylene and ethanol, is required for practical implementation of this technology. Herein, we successfully obtained the partial current density for C2+ over cupric oxide (CuO) nanoparticles on gas diffusion electrodes in neutral electrolytes to a value of 1.7 A cm−2. Additionally, the faradaic efficiency for multicarbon products increased with the current density below Jtotal = 2000 mA cm−2 and reached up to 76%. Our ultra-high-rate CO2 electrolysis was achieved by optimizing the ordinary components, including oxide-derived Cu nanoparticles and carbon-based gas-diffusion electrodes, and their assembly as the cathode. In the other word, the present study succeeded in eliciting the potential performance of an ordinary combination. Although more detailed analyses are needed, our electrodes will teach us essential requirements for ultra-high-rate electrolysis. Additionally, the obtained knowledge through the detailed analyses of our electrodes would be widely utilized to develop novel materials, such as catalysts and electrodes. |
The use of gas diffusion electrodes (GDEs) that allow the CO2RR to occur at the solid-catalyst/liquid-electrolyte/gaseous-CO2 interface, effectively accelerate the CO2RR by solving the problem of the mass transport limitation due to the inherently low diffusion and solubility of CO2 in water.15,16 A catalyst layer (CL) composed of metal Cu particles and ionomers on carbon-based microporous layers has been widely studied as the standard cathode for gaseous CO2RR to produce C2+ organics, such as C2H4, C2H5OH, C3H7OH and CH3COOH.8,17–19 Although the partial current density for C2+ (jC2+) has still remained below 500 mA cm−2 in many studies,18,20–27 a few studies have developed novel cathode assemblies or components toward an ultra-high-rate (UHR-) CO2RR with over jC2+ = 1 A cm−2. Sargent et al. successfully achieved a jC2+ over 1.2 A cm−2 in 7 M KOH using an electrode fabricated by dropping Cu nanoparticles on a sputtered Cu film/porous polytetrafluoroethylene (PTFE).28 Additionally, in this work, the bulk heterojunction of catalysts and ionomers decoupled the pathways for gases, electrons and ions, which resulted in a jC2+ over 1 A cm−2. Wang and co-authors reported that a fluorine-modified copper catalyst deposited on GDEs exhibited a total current density (Jtotal) of 1600 mA cm−2 with a faradaic efficiency for C2+ (FEC2+) of over 80% (i.e., jC2+ = over 1.2 A cm−2) in 0.75 M KOH.29 Very recently, Zheng et al. obtained jC2+ = 0.91 A cm−2 using nitrogen-engineered Cu nanoparticles in 1 M KOH.10
In these reports, UHR-CO2RR has been achieved using high alkaline solutions to suppress H2 evolution, which competes with the CO2RR. However, in addition to the toxicity and corrosivity, alkaline solutions quickly absorb CO2 and convert it to electrochemically inert (bi)carbonate.30–33 Based on these drawbacks, UHR electrolysis in neutral electrolytes is thus required from a practical viewpoint. Furthermore, these previous studies succeeded in achieving UHR-CO2RR by developing novel components in the three-phase interfaces (i.e., catalysts or electrodes). However, the obtainable maximum jC2+ value is still unknown when optimizing the standard components; CuO nanoparticle catalysts (CuONPs), which are known to serve as oxide-derived Cu catalysts for effective C2+ production,34–36 and carbon-based gas diffusion layers. An electrochemical system that can elicit the maximum jC2+ using only standard components will become the primary platform for future studies to design materials and interfaces for UHR-CO2RR. In other words, the design guidelines obtained from the UHR-CO2RR with only ordinary materials would be quite general, and thus, they are expected to be widely applicable to novel materials and systems.
In the present work, we attempt to pursue a high partial current density for C2+ (jC2+) using neutral electrolytes by optimizing the standard components and properly assembling them as the cathode. In particular, we achieved a record current density of 1.7 A cm−2 for multicarbon products in neutral electrolytes using synthetic or commercial CuONPs and carbon-based gas diffusion layers. Firstly, we will explain the details in characterizations and CO2RR activity of our optimized electrode, and let us discuss the mechanism of UHR-CO2RR. Then, the requirements for our UHR electrolysis will be discussed with varying components and conditions.
Fig. 1 (a) Representative TEM image of the synthesized CuONPs. (b) Top-view and (c and d) cross-sectional SEM micrographs of CuONPs-1.7/GDE. |
We measured the CO2RR activity of the Cu-based GDE catalysts. An in-house-built three-compartment electrochemical cell was used for CO2RR evaluation (Fig. S12, ESI†).37,38Jtotalvs. potential (U) curves were obtained in 1 M KCl under CO2 and argon conditions, as shown in Fig. S13 (ESI†). The onset potential of the cathodic current for bare and Cu-modified GDEs were −1.6 and −1.2 V vs. Ag/AgCl under CO2, respectively. When the Cu amount was increased from CuONPs-0.34/GDE to CuONPs-3.1/GDE, the cathodic current under CO2 gradually increased.
The CO2RR products in outlet gas streams from the gas chamber and electrolytes from the catholyte room were analysed using gas chromatography (Fig. S14, ESI†) and 1H nuclear magnetic resonance (NMR) spectroscopy (Fig. S15, ESI†) after constant-current electrolysis, respectively. Fig. 2a, b and Fig. S16 (ESI†) show the FE and partial current density of CO2RR products against the applied potential for CuONPs-1.7/GDE in 1 M KCl, respectively. Under low current conditions (Jtotal = 200 mA cm−2) at −1.36 V vs. Ag/AgCl electrode, CO was the main product (FECO > 50%), and the FEC2H4 was below 20%. The FE for C2+ products increased and reached 76% at −1.69 V vs. Ag/AgCl electrode (Jtotal = 1600 mA cm−2) when the applied potential was shifted to the negative side. In particular, the FEs for C2H4 and C2H5OH were 46% and 22% at −1.76 V vs. Ag/AgCl electrode, respectively. The partial current density for C2+ products (C2H4, C2H5OH, C3H7OH, and CH3COOH) was over jC2+ = 1.7 A cm−2 at −1.76 V vs. Ag/AgCl electrode. On the contrary, the FEs for C1 products (CO, CH4 and HCOOH) and H2 were 8.8 and 14% at −1.76 V vs. Ag/AgCl electrode, respectively. Fig. 2c and Table S2 (ESI†) show a comparison of jC2+ and FEC2+ in neutral electrolytes among the reported values. The partial current density for C2+ production (jC2+ = over 1.7 A cm−2) is a record value, and the FEC2+ is also one of the highest. Therefore, the present study succeeded in eliciting the potential performance of an ordinary combination for UHR-CO2RR (i.e., CuO-derived Cu nanoparticles and carbon-based GDEs). Significantly, such a high jC2+ value for our assembly was remarkably reproducible even using different batches of electrodes, as shown in Fig. S17 and Table S3 (ESI†). When the applied potential was changed to −1.79 V vs. Ag/AgCl electrode (Jtotal = 2800 mA cm−2), the FE for C2+ products suddenly decreased, and the FE for H2 increased. The total FEs at −1.79 V vs. Ag/AgCl electrode did not reach 100% (ca. 80%) because the bubbles contained CO2 substrate, CO2RR products, and H2 formed at the cathode surface and diffused into the catholyte chamber. The applied potential was compensated using the current interrupt method as written in the experimental section39–41 (see Fig. S18 for the values versus reversible hydrogen electrode, ESI†). However, as accurate potential corrections are difficult under ultra-high-current conditions, we must note that the applied potentials are still referenced values.
Fig. 2 (a) FEs of the CO2RR products under various applied potentials on CuONPs-1.7/GDE in 1 M KCl. Reaction time: 30 min. (b) Partial current density for C2+, C1 and H2 on CuONPs-1.7/GDE in 1 M KCl. (c) Comparison of FEC2+ and jC2+ for CuONPs-1.7/GDE with reported Cu-based catalysts in neutral solution (for the detailed references, see Table S2, ESI†). (d) FEs for the CO2RR products under various applied potentials on CuONPs-1.7/GDE in 1 M KOH. Reaction time: 30 min. (e) Comparison of FEC2+ in 1 M KOH and in 1 M KCl under various Jtotal over CuONPs-1.7/GDE. (f) Comparison of FEC2+ and jC2+ for CuONPs-1.7/GDE with reported Cu-based catalysts in alkaline solution (for the detailed references, see Table S5, ESI†). |
Even when using neutral solutions, the local pH at the surface of the catalyst increases significantly during UHR-CO2RR due to the consumption of protons with an increase in the CO2RR current density. Our rough calculation and previous reports showed that the local pH becomes almost the same in neutral and alkaline electrolytes during high-current electrolysis (Fig. S19 and Table S4, ESI†). Therefore, the cathode reaction proceeded under alkaline environments for UHR-CO2RR even with neutral solutions; thus, the fabricated cathode was expected to be suitable for UHR-CO2RR even in bulk alkaline solutions. Fig. 2d shows the FEs of CO2RR products against the applied potential for CuONPs-1.7/GDE in 1 M KOH alkaline electrolyte, and Fig. 2e shows a comparison of the FEs of C2+ products for CuONPs-1.7/GDE in 1 M KCl and 1 M KOH solutions. The FEs for CO2RR activity in alkaline and neutral solutions were almost the same in the range of Jtotal = 200–2400 mA cm−2, which indicates that the local pH around the CL in 1 M KOH and 1 M KCl solutions was not dependent on the bulk pH during electrolysis at Jtotal = 200 mA cm−2. The partial current density for C2+ products in alkaline solutions reached jC2+ = 1.8 A cm−2, which is a record value in alkaline solution, as compared and listed in Fig. 2f and Table S5 (ESI†) respectively, i.e., this electrode can exhibit a record jC2+ in both neutral and alkaline solutions. The mass activity, turnover frequency (TOF) and energy efficiency values are also summarized in Fig. S20 and Table S6 (ESI†) as reference data.
The volcano-type relation between FEC2+and current density was observed in the present UHR-CO2RR system. Especially, the FEC2+ increased with increasing current density below Jtotal = 2000 mA cm−2, and thus, a trade-on relationship between rate and selectivity was exhibited in this current region, as shown in Fig. 2a. Let us discuss the mechanism for the increase of FEC2+ with the current density. It should be noted here that it is difficult to strictly distinguish whether the selectivity toward the product depends on the current density (i.e., reaction rate) or on the applied potential because we cannot control each parameter independently. We here explain the product selectivity from the viewpoint of current density: the effect of (1) CO partial pressure, (2) local pH,42 and (3) in situ formation of the three-phase interface. The CO partial pressure in the CL became high with an increase in the current density, which facilitated CO dimerization. Jun Li et al. controlled local CO availability at the catalyst–electrolyte interface and demonstrated that a high concentration of CO accelerated CO dimerization, and caused a higher current density with C2+ production.43 Jing Li and co-workers reported that a high CO partial pressure suppresses H* adsorption, which leads to lowering of the hydrogen evolution reaction (HER).44 It should also be noted here that the FE for C2H5OH against C2H4 increased with increasing Jtotal (Fig. S21, ESI†). It has been reported that high CO pressure on a Cu surface favours C2H5OH production over C2H4 because the *CHCHOH intermediate for C2H5OH becomes more stable than *CCH for C2H4 under high CO coverage.43,45 This is because CCH* and *CHCHOH are adsorbed on Cu sites with side-on and end-on configurations, respectively, so that CCH* requires a larger occupied molecular area than *CHCHOH. This change in selectivity between C2H4 and C2H5OH supports that the CO local pressure depended on the current density. The increase in the local pH at the surface of the CL due to proton consumption by the CO2RR suppressed HER (vide supra). The pH at the electrode surfaces is known to reach up to 11.9 under electrolysis at just Jtotal = 50 mA cm−2 in neutral electrolytes.15,46 Considering that the almost same product selectivity was obtained in neutral (1 M KCl) and alkaline (1 M KOH) solutions between Jtotal = 200–2400 mA cm−2, the local pH in KCl solution became similar to that in KOH solution under electrolysis with Jtotal of just Jtotal = 200 mA cm−2. Further increase in the current density (Jtotal > 1000 mA cm−2) would lead to a pH of over 15 at the surface of the catalyst.7 Moreover, we propose that the in situ formation of the gas-transport pathway depending on the current density occurs during electrolysis. Evolved gases, such as H2, would migrate through gas-transport channels in the CL; therefore, the area of the CO2RR active interface (three-phase interfaces) could be expanded with an increase in the current density (see Fig. S22, ESI†).
It is important to identify the origin of the decrease in the CO2RR activity by application of Jtotal = 2800 mA cm−2 for CuONPs-1.7/GDE (Fig. 2a and d) toward a further increase in the production rate. It was assumed that the three-phase interface was broken over Jtotal = 2800 mA cm−2, which resulted in a severe decrease in the FE for C2+ products. Fig. 3 shows cross-sectional SEM images after application of Jtotal = 1000, 2000, and 2800 mA cm−2 (30 min) in 1 M KCl, respectively. The particle size of the CuONPs became about 10% smaller after electrolysis due to the reduction of CuONPs to Cu(I)/Cu(0)NPs (Fig. S23, ESI†). Cross-sectional elemental energy dispersive X-ray spectroscopy (EDX) mapping show K+ derived from electrolytes entered the MPL layer. One-tenth and one-fourth of the MPL layers were flooded after application of Jtotal = 1000 mA cm−2 and 2000 mA cm−2, respectively. In contrast, after the degradation with application of Jtotal = 2800 mA cm−2 (30 min), the K+ distribution was shown to reach half of the MPL layer (Fig. 3e). In addition, a residue of K+ crystals (KHCO3 and K2CO3, KCl and KOH) was observed throughout the MPL after application of Jtotal = 2800 mA cm−2 (Fig. S24, ESI†). Electrowetting effects, which decrease the solid–liquid interfacial tension and facilitate electrolyte penetration into the CL and MPL, are known to be emphasized under application of highly negative potentials.47–49 Although the evolved gases may push out and expand the gas-transport pathways (as the in situ formation of the three-phase interface), flooding of the MPL due to electrowetting may exceed the rate of pushing out the electrolyte. The flooded CL (and MPL) produced only H2 due to a lack of CO2 transportation, and thus, H2 production should become dominant at Jtotal = 2800 mA cm−2. This assumption is also supported by the diffusion of evolved gases into the catholyte chamber at Jtotal = 2800 mA cm−2 (vide supra).
Fig. S25 and Table S7 (ESI†) show the FE of CO2RR products under Jtotal = 2000 mA cm−2 when changing catalysts, the thickness of CL, electrolytes, and GDEs. One can notice that several kinds of commercial GDEs are available for UHR-CO2RR. Additionally, the high jC2+ of over 1.3 A cm−2 was obtained using 1 M KHCO3, which has been widely used as a neutral electrolyte for CO2RR. The slight decrease of jC2+ in 1 M KHCO3 is likely due to an acceleration of salt precipitation in MPL in the bicarbonate-rich electrolytes.
Surprisingly, UHR-CO2RR with jC2+ of over 1.4 A cm−2 was realized even when using commercial CuONPs with the size of 50 nm. These results indicated that the GDE and catalyst used in this study are not unique, and the optimized assembly elicits their undermined potential. Further studies are currently in progress to determine the detailed requirements for catalysts in our laboratory.
Among these parameters, Fig. S25 (ESI†) reveals that the thickness of CL is one high-sensitivity factor in determining maximum jC2+. The optimal thickness of CL was required for UHR-CO2RR, not too thin and not too thick. Thus, we investigated the detail of the dependence of CL thickness. The CO2RR was investigated for CuONPs-3.1/GDE and CuONPs-0.34/GDE (detailed product selectivities are shown in Fig. S26 and S27, ESI†) and the standard CuONPs-1.7/GDE (see the FEs in Fig. 2a) in 1 M KCl. The maximum jC2+ for CuONPs-3.1/GDE was achieved at Jtotal = 2400 mA cm−2, which was the same as that for CuONPs-1.7/GDE. The FEs for C2+ over CuONPs-3.1/GDE were 7% lower than those over CuONPs-1.7/GDE, although there was a similar tendency of the potential dependence (Fig. 4a). The FEH2 over CuONPs-3.1/GDE was higher than that over CuONPs-1.7/GDE for all current density regions (Fig. S28, ESI†). These results indicated that a lack of Cu active sites did not cause a decrease in FEC2+ over Jtotal = 2400 mA cm−2 for CuONPs-1.7/GDE. It was presumed that the Cu sites immersed in the electrolyte cannot serve as CO2RR sites but only as HER sites, and thus, selectivity toward H2 production increased with the amount of Cu (Fig. 4d). Next, the Cu amount on the GDE was decreased to 0.34 mg cm−2. Although the FEC2+ over CuONPs-0.34/GDE was greater than 80% at Jtotal = 1000 mA cm−2, it decreased to below 60% at Jtotal = 1600 mA cm−2. In contrast, when Jtotal was varied from 1000 mA cm−2 to 1600 mA cm−2, the FEs for CH4 and H2 increased from 1.6% to 14%, and from 5.4% to 17%, respectively. This is because the lack of catalytic Cu sites would result in relative enhancement of the HER on carbon particles that compose the MPL. In other words, the three-phase interface formed in CuONPs-0.34/GDE was thinner than that in CuONPs-1.7/GDE, resulting that the current density with the highest FEC2+ was smaller than the optimal for CuONPs-1.7/GDE (Fig. 4b and c). A significant increase in FE for CH4 with an increase in current density was observed on CuONPS-0.34/GDE. This tendency for CH4 was basically consistent with the reported potential dependence of the selectivity of the CO2RR on Cu-based catalysts.50–52 The same trends in the catalytic activity under various amounts of catalyst loadings were observed in 1 M KOH electrolyte (Fig. S29–S31, ESI†).
Fig. S32 (ESI†) shows the FEs depending on CL thickness at 500, 1000, 1600 and 2000 mA cm−2, to further demonstrate the importance of CL thickness for maximum jC2+. The thin CL (CuONPs-0.34) showed the best FE for C2+, and FEC2+ monotonically decreased with increasing the thickness of CL at Jtotal = 500 and 1000 mA cm−2. In contrast to FEC2+, the total FE for CO2RR was almost constant (85–90%) at these low current density regions, even when varying CL thickness. This is likely because CO partial pressure was not enough to form C–C bonds for thicker CL at the low current density regions. When increasing current density to Jtotal = 1600 mA cm−2, the FEC2+ for CuONPs-1.7 became the best in these three electrodes. Thus, our results clearly indicated that an optimal CL thickness for C2+ evolution depends on Jtotal.
Although our study is not the first paper to optimize and discuss each parameter, including the size of catalysts, porosity, and thickness of CL, we think that these interdependent parameters were optimized simultaneously, likely resulting in the highest jC2+ in this work. However, we must note that some requirements, especially for catalysts, toward UHR-CO2RR are still veiled. This work is focused on demonstrating the potential of ordinary catalysts; providing more detailed requirements will be left for a future paper.
Fig. 5b and Fig. S33 (ESI†) show the FEs for CO2RR products at Jtotal = 1600 mA cm−2 in 1 M KCl as a function of the electrolysis time. The change in the FEs for CO2RR products under electrolysis at Jtotal = 1600 mA cm−2 was almost negligible for over 2 h. Although the FE for C2+ and H2 then began to decrease and increase from 2 h, respectively, the FEC2+ production remained over 75% for 2.5 h and 60% for over 6 h. The electrode was analysed after the decrease of FEs for C2+ to clarify the origin of activity decrease. SEM observation (Fig. S34, ESI†) after electrolysis for 6 h at 1600 mA cm−2 showed that the morphology of CLs did not significantly vary, and most of the MPL was flooded in a similar manner to that described in previous reports.53,54 The recycling test was also conducted. Fig. S35 (ESI†) shows that the FEC2+ value was recovered by washing the electrode with pure water after the decrease. This result also indicates that the decrease in FEC2+ was due not to the degradation of catalysts but to the salt precipitation induced by flooded MPL. The stability dependence on the applied current density was investigated next (Fig. 5a–c and Fig. S33, ESI†). The FEC2+ of 70–80% almost unvaried for over 3 h and 1.5 h at Jtotal = 1000 mA cm−2 and Jtotal = 2000 mA cm−2, respectively. These results indicated a trade-off relationship between stability and current density; however, the decrease in FEC2+ within a short timescale was due to flooding of the CL and MPLs. Therefore, we can expect that the stability could be improved by the addition of hydrophobic materials and modulation of the porosity in the CL and MPL. The present electrochemical system and conditions will become one of the standard platforms toward these future studies.
Fig. 5 FEC2+ stability tests with respect to the current density over CuONPs-1.0/GDE in 1 M KCl under (a) Jtotal = 1000 mA cm−2, (b) 1600 mA cm−2 and (c) 2000 mA cm−2. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2ey00035k |
This journal is © The Royal Society of Chemistry 2023 |