
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39227576
51928
10.1038/s41467-024-51928-z
Article
Multivalent Cu sites synergistically adjust carbonaceous intermediates adsorption for electrocatalytic ethanol production
Chen Xiao 12
Jia Shuaiqiang sqjia@chem.ecnu.edu.cn

12
Zhai Jianxin 12
Jiao Jiapeng 12
Dong Mengke 12
Xue Cheng 12
Deng Ting 12
Cheng Hailian 12
Xia Zhanghui 12
Chen Chunjun 12
Xing Xueqing 3
Zeng Jianrong 4
http://orcid.org/0000-0001-6266-8290
Wu Haihong hhwu@chem.ecnu.edu.cn

12
http://orcid.org/0009-0006-9078-8322
He Mingyuan 12
http://orcid.org/0000-0003-0440-809X
Han Buxing hanbx@iccas.ac.cn

125
1 https://ror.org/02n96ep67 grid.22069.3f 0000 0004 0369 6365 Shanghai Key Laboratory of Green Chemistry and Chemical Processes, State Key Laboratory of Petroleum Molecular & Process Engineering, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
2 grid.22069.3f 0000 0004 0369 6365 Institute of Eco-Chongming, 20 Cuiniao Road, Chenjia Town, Chongming District, Shanghai, China
3 grid.9227.e 0000000119573309 Beijing Synchrotron Radiation Facility, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, China
4 grid.9227.e 0000000119573309 Shanghai Synchrotron Radiation Facility, Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai, China
5 grid.9227.e 0000000119573309 Beijing National Laboratory for Molecular Sciences, CAS Key Laboratory of Colloid and Interface and Thermodynamics, CAS Research/Education Center for Excellence in Molecular Sciences, Center for Carbon Neutral Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
3 9 2024
3 9 2024
2024
15 769124 4 2024
20 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Copper (Cu)-based catalysts show promise for electrocatalytic CO2 reduction (CO2RR) to multi-carbon alcohols, but thermodynamic constraints lead to competitive hydrocarbon (e.g., ethylene) production. Achieving selective ethanol production with high Faradaic efficiency (FE) and current density is still challenging. Here we show a multivalent Cu-based catalyst, Cu-2,3,7,8-tetraaminophenazine-1,4,6,9-tetraone (Cu-TAPT) with Cu2+ and Cu+ atomic ratio of about 1:2 for CO2RR. Cu-TAPT exhibits an ethanol FE of 54.3 ± 3% at an industrial-scale current density of 429 mA cm−2, with the ethanol-to-ethylene ratio reaching 3.14:1. Experimental and theoretical calculations collectively unveil that the catalyst is stable during CO2RR, resulting from suitable coordination of the Cu2+ and Cu+ with the functional groups in TAPT. Additionally, mechanism studies show that the increased ethanol selectivity originates from synergy of multivalent Cu sites, which can promote asymmetric C–C coupling and adjust the adsorption strength of different carbonaceous intermediates, favoring hydroxy-containing C2 intermediate (*HCCHOH) formation and formation of ethanol.

This study proposes a multivalent Cu-based catalyst (Cu-TAPT) for electrocatalytic CO2RR, which promotes specific intermediate formation by the synergy of Cu2+ and Cu+ sites, thereby achieving a 54.3% Faradaic efficiency of ethanol at 429 mA cm−2.

Subject terms

Electrocatalysis
Sustainability
the National Key R&D Program of China (2023YFA1507901, 2020YFA0710201), the National Natural Science Foundation of China (22293015, 22121002), the China Postdoctoral Science Foundation (2023M731096), and the Fundamental Research Funds for the Central Universities.issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The electrocatalytic carbon dioxide reduction (CO2RR) to produce high-value chemicals and fuels using renewable electricity provides a promising carbon-neutral route1–3. Multi-carbon (C2) products, such as ethylene and ethanol, have the prospect of higher energy density and market value than C1 products4. Notably, significant progress has been made in producing gaseous ethylene via the CO2RR5,6. In contrast, there has been limited progress to date in the electrochemical conversion of CO2 into liquid ethanol with high volumetric energy density7. Due to the dual constraints of thermodynamics and kinetics, CO2 electrocatalysts do not provide the majority of ethanol production, preferring ethylene instead8,9. Therefore, increasing ethanol selectivity relative to hydrocarbon is highly desired, but remains a challenge.

Previous studies have suggested that ethanol production can be optimized by modulating and stabilizing the key intermediate (*HCCOH) to undergo the hydrogenation pathway of *HCCOH to *HCCHOH (ethanol pathway) instead of dehydroxylation way to *HCC (ethylene pathway) over the copper (Cu)-based catalysts10–13. This optimization is achieved mainly through strategies such as tandem catalysis14, metal-support interactions15, multi-element doping12,13,16–18, interfacial interaction11 and various surface binding sites19,20. However, such Cu-based catalysts usually undergo electrochemical reduction under negative potential, which hinders their application in ethanol production. We hypothesize that ethanol production can be enhanced by synergistic effect of multivalent Cu sites with different adsorption capacities for carbonaceous intermediates, and improving their electrochemical stability under CO2RR conditions via suitable coordination interaction.

In this study, based on the strategy proposed above, we constructed a multivalent Cu-based coordination polymer catalyst Cu-2,3,7,8-tetraaminophenazine-1,4,6,9-tetraone (Cu-TAPT) with stable Cu2+ and Cu+ sites, and applied it as the catalyst in the ethanol synthesis via CO2RR. With the cooperative promotion of Cu2+ and Cu+ sites, the Faradaic efficiency (FE) ratio of ethanol to ethylene over the Cu-TAPT catalyst was 3.14:1, significantly outperforming the reference catalyst Cu-2,3,5,6-tetraaminobenzoquinone (Cu-TABQ). In situ X-ray absorption and Raman spectroscopic analysis suggested that the multivalent Cu sites are stabilized by the double confinement of N and O. Density functional theory (DFT) calculations further revealed that the synergy of multivalent Cu sites played a crucial role in hindering the breakage of the C–O bond in *HCCOH and facilitating the hydrogenation of C atoms. This process effectively retains the OH group, thereby promoting the selectivity of ethanol in CO2RR.

Results

Material synthesis and characterization

Cu-TAPT was obtained through the coordination reaction involving 2,3,5,6-Tetraaminobenzoquinone (TABQ) ligands and equimolar Cu2+ ions in deionized water at 120 °C. Cu-TABQ was prepared by substituting deionized water with dimethyl sulfoxide, following the similar procedure (Fig. 1a and Supplementary Fig. 1). Due to the solvent effect, the dimerization of the ligands in deionized water caused the partial reduction of Cu2+ into Cu+, resulting in structural differences. Powder X-ray diffraction (XRD) showed a broad diffraction peak of Cu-TAPT, suggesting an amorphous nature (Supplementary Fig. 2), which facilitates exposure of more active sites, as confirmed by the Brunauer–Emmett–Teller (BET) results (Supplementary Fig. 3)21. Scanning electron microscopy (SEM) images revealed distinct granule morphology of Cu-TAPT, showing the difference of Cu-TABQ and TABQ (Supplementary Fig. 4). Energy-dispersive X-ray (EDX) element mappings (Supplementary Fig. 5 and Supplementary Table 1) confirmed that Cu, N and O atoms were homogeneous distribution throughout Cu-TAPT. The contents of Cu in Cu-TAPT and Cu-TABQ obtained from inductively coupled plasma optical emission spectrometer (ICP-OES) were 27.6 and 28.2 wt%, respectively (Supplementary Table 2). Transmission electron microscopy (TEM) image showed the amorphous morphology and the absence of nanoparticles or clusters in Cu-TAPT (Supplementary Fig. 6), which is consistent with the XRD results. According to the high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images of Cu-TAPT and Cu-TABQ, ultrahigh-density bright dots were observed, suggesting atomically dispersion of Cu sites due to coordination (Fig. 1b and Supplementary Fig. 7).Fig. 1 Material synthesis and fine structural characterization.

a Schematic illustration of the Cu-TAPT and Cu-TABQ synthetic process. b HAADF-STEM image of Cu-TAPT catalyst. c, d N 1 s spectra c and O 1 s spectra d of the TABQ, Cu-TABQ and Cu-TAPT. e Cu 2p spectra of Cu-TAPT catalyst. f The first derivatives of the XANES region of Cu foil, Cu-TAPT, Cu-TABQ, Cu2O, CuO, CuPc and the inset is the content of Cu+ and Cu2+ in Cu-TAPT from first-order derivative of XANES. g Wavelet transforms (WT) for the k3-weighted Cu K-edge EXAFS signals of Cu-TAPT.

According to high-resolution X-ray photoelectron spectroscopy (XPS), the shifted binding energy of N and O in Cu-TAPT suggested that the electrons of N and O atoms were delocalized and were coordinated with Cu (Fig. 1c, d)22, which was also confirmed by the Fourier transform infrared spectroscopy (FT-IR) results (Supplementary Fig. 8). Different from the typical satellite peaks of Cu2+ in Cu-TABQ (Supplementary Fig. 9)23, the ratio of Cu2+ and Cu+ in Cu-TAPT was 1:1.84 (Fig. 1e, Supplementary Fig. 10 and Supplementary Table 3), agreed well with the theoretically-predicted value 1:2. The absorption threshold position in X-ray absorption near edge structure (XANES) of Cu-TABQ was basically the same as that of CuO and copper phthalocyanine (CuPc), indicating the similar oxidation state of +224. Interestingly, Cu-TAPT exhibited an absorption edge position between CuO and Cu2O, suggesting the coexistence of bivalent and monovalent of Cu and the ratio was 1:1.85 (Fig. 1f and Supplementary Fig. 11a), in accordance with the XPS results. The Fourier transforms extended X-ray absorption fine structure (EXAFS) spectra of Cu-TAPT displayed a prominent peak that is different from Cu–O and Cu–N and has a higher coordination number. Moreover, no peaks attributed to metal Cu–Cu bonds were observed around 2.2 Å, indicating that Cu sites are atomically distributed in Cu-TAPT (Supplementary Fig. 11b)25. The fitting results of EXAFS spectra (Supplementary Fig. 12 and Supplementary Table 4) showed that the average coordination number (CN) of Cu–O and Cu–N in Cu-TAPT was 4.10 and 4.22, respectively, which was different from that of Cu-TABQ (CNCu–N = 4.04) (Supplementary Fig. 13 and Supplementary Table 4). The differences in these structures can also be clearly observed from wavelet transform (WT) spectra (Fig. 1g and Supplementary Fig. 14). These characterizations agree with the building unit of the Cu-TAPT catalyst, which comprises Cu in different valence states and TAPT monomers extending to both ends through Cu–N and Cu–O coordination bonds. Such adjacent Cu sites with different valence states are expected to have synergistic effect for CO2RR, as substantiated in the subsequent section.

Electrocatalytic CO2RR performance

The electrocatalytic CO2RR activity and selectivity of the as-synthesized catalysts were tested in a flow cell reactor with 1 M KOH electrolyte, and the product was detected by gas chromatography (GC) and 1H nuclear magnetic resonance (NMR) spectroscopy. Cu-TAPT exhibited a lower onset potential than that of Cu-TABQ in linear sweep voltammetry (LSV) curves, suggesting that the multivalent catalyst could effectively reduce the overpotential of CO2RR (Supplementary Fig. 15). FEethanol increased and then decreased with applied potentials over Cu-TAPT in the range from −0.70 V to −0.85 V vs. RHE (Fig. 2a and Supplementary Figs. 16,17). At −0.80 V vs. RHE, FEethanol and FEC2+ reached its maxima of 54.3 ± 3% and 84.9 ± 2%, respectively (Supplementary Figs. 17). Cu-TAPT showed 2.31 times increase in FEethanol and 1.43 times increase in FEC2+ compared to the Cu-TABQ, which only reached the highest FEethanol and FEC2+ of 23.5% and 60%, respectively, at a more negative potential of −0.98 V vs. RHE (Supplementary Figs. 17-19). In particular, the FE ratio of ethanol to ethylene was almost more than 2 under the applied potential range, and even exceeds 3.1 times at −0.8 V vs. RHE (Fig. 2b), suggesting a preference for ethanol on Cu-TAPT. Compared with Cu-TABQ, Cu-TAPT exhibited greater current density, and the total current density and ethanol partial current density could reach 429 mA cm−2 and 233 mA cm−2 at −0.8 V vs. RHE, respectively (Supplementary Figs. 20,21). The C2+ and ethanol formation rates over Cu-TAPT increased with applied potentials and reached the maximum of 1,205 and 724 μmol h−1 cm−2 at −0.8 V vs. RHE, which were about 2.53 times and 4.16 times higher than that on Cu-TABQ at −0.98 V vs. RHE, respectively (Fig. 2c). The Cu-TAPT catalyst outperformed most of start-of-the-art Cu-based catalysts in FEethanol or current density or formation rate (Supplementary Table 5). Moreover, the mass activity, specific activity and turnover frequency (TOF) of Cu-TAPT were significantly better than those of Cu-TABQ (Supplementary Figs. 22,23). Based on the results above, we propose a mechanism involving concerted catalysis by different valence Cu sites to elucidate the variations in catalytic performance between Cu-TAPT and Cu-TABQ (Fig. 2d). Cu sites in different valence states exhibit distinct adsorption capacities for intermediates, thereby facilitating coupling process, promoting ethanol production, and inhibiting the formation of hydrocarbon products. Further specific details are discussed in the following sections.Fig. 2 Electrocatalytic CO2RR performance.

a FEs of all products under different applied potentials over Cu-TAPT in 1 M KOH electrolyte. b FE ratio of ethanol/ethylene at different potentials over Cu-TAPT. c, The formation rates of C2+ products and ethanol under different applied potentials over Cu-TAPT and Cu-TABQ in 1 M KOH electrolyte. d, Schematic illustration of C − C coupling formation process on Cu-TAPT based on the synergistic effect of Cu2+ and Cu+. The ball of dark blue, red, gray, azure, pink and white denote Cu2+, Cu+, C, N, O and H, respectively. e Time-dependent total current density, FEethanol and FEC2+ on Cu-TAPT were recorded at −0.8 V vs. RHE for at least 30 h. Error bars in accordance with the standard deviation of at least three independent measurements. All the potentials involved in this figure have been iR corrected (R = 1.89 Ω). The potential without iR is shown in Supplementary Fig. 17.

The properties of the electrode and electrolyte interfaces were further investigated via electrochemical double-layer capacitance and electrochemical impedance spectroscopy (EIS). The electrochemical active surface area (ECSA) of Cu-TAPT was 2.52 times of that of Cu-TABQ (Supplementary Fig. 24 and Table 6)18. The Nyquist plots showed that Cu-TAPT had a smaller charge transfer resistance (Rct) than Cu-TABQ (Supplementary Fig. 25 and Supplementary Table 7), which further reflected the more favorable charge transfer kinetics over Cu-TAPT for CO2RR26. Similarly, Cu-TAPT exhibited greater ionic conductivity, which is favorable to electrocatalytic reaction (Supplementary Table 7)27. In addition, the long-term durability of Cu-TAPT was also evaluated by a continuous CO2RR at −0.8 V vs. RHE for 30 h, and the FEethanol and FEC2+ products did not decrease obviously (Fig. 2e). As a comparison, the performance results for Cu-TABQ revealed a rapid decline in the FE for C2+ products within 7 h (Supplementary Fig. 26).

Various characterizations were performed to ulteriorly unveil the structural stability after CO2RR. The XRD pattern of Cu-TAPT indicated no metallic Cu appeared after CO2RR (Supplementary Fig. 27). Meanwhile, the XPS results suggested that the ratio of Cu2+ and Cu+ and N 1 s and O 1 s spectra in Cu-TAPT remains basically unchanged before and after the CO2RR process, whereas the molecular structure of Cu-TABQ could not be stabilized, and almost decomposes (Supplementary Fig. 28). Cu-TAPT could maintain the coexistence of Cu2+ and Cu+ with a ratio of 1:1.78 even after a long electrolysis process for 30 h (Supplementary Fig. 29). Additionally, the morphology of Cu-TAPT generally remained unchanged after CO2RR (Supplementary Fig. 30), which further indicated the excellent stability. Taken together, the Cu-TAPT catalyst has promising potential of application for the conversion of CO2 to ethanol products.

In situ characterization

To investigate the oxidation state and structural evolutions of catalysts during the CO2RR process, in situ XAS studies were carried out via a customized flow cell (Supplementary Fig. 31). The in situ XANES spectra of Cu-TAPT revealed an unchanged Cu oxidation state (Cu2+ and Cu+) for 2 h under −0.98 V vs. RHE (Fig. 3a, b), indicating that the active Cu sites in Cu-TAPT were stable during CO2RR process. Consistent with XANES results, no Cu–Cu bond at around 2.24 Å was discerned in the EXAFS (Supplementary Fig. 32), demonstrating that Cu remained as coordinated and atomically dispersed state of Cu-TAPT under CO2RR condition28. These results were also confirmed by the potential-dependence results (Fig. 3c and Supplementary Fig. 33). On the contrary, the oxidation state of Cu decreases with decreasing potential from OCP to −0.98 V vs. RHE of Cu-TABQ (Supplementary Fig. 34). The major peak ascribed to Cu–Cu bonds at approximately 2.2 Å emerged, implying that the catalysts underwent severe electrochemical reconstruction and transformed into metallic Cu particles. Linear combination fitting (LCF) of the XANES spectra was further conducted to calculate the chemical valence of Cu-TAPT, finding no metallic Cu (Fig. 3d and Supplementary Fig. 35). Likewise, in situ Raman spectroscopy of Cu-TAPT was conducted at the same constant potential of −0.98 V vs. RHE under CO2 atmosphere for 2 h (Fig. 3e). No considerable signal change was discovered for Cu-TPAT throughout extended operation, implying that the coordination structure of Cu-TAPT had high stability. We attribute these results to mutual constraints between Cu, N, and O atoms arising from asymmetric coordination29–32.Fig. 3 In situ characterization of Cu-TAPT during CO2RR.

a, b Time-dependent in situ Cu K-edge XANES a and first-order derivatives of the XANES region b of Cu-TAPT catalyst recorded at −0.98 V vs. RHE under CO2 atmosphere. c Potential-dependence of in situ first-order derivatives of the XANES spectra of the Cu-TAPT under CO2RR test. d Variation in the valence state of Cu species in Cu-TAPT catalyst with time. e In situ Raman spectra of Cu-TAPT catalyst operating at −0.98 V vs. RHE under CO2 atmosphere with a 532 nm laser. f Total energy as a function of simulation time during the AIMD simulations under 300 K and the electric field of −2.3 V/Angstorm. All the potentials involved in this figure are not iR corrected.

In addition, DFT and ab initio molecular dynamics (AIMD) calculations were used to further interpret the enhanced stability of Cu-TAPT (Model information in Supplementary Data 1). The crystal orbital Hamilton populations (COHP) were investigated for Cu and N or O atom bond strength of Cu-TAPT, Cu-TABQ and CuO (Supplementary Figs. 36,37 and Supplementary Table 8). According to the COHP theory, the negative and positive values of integral crystal orbital Hamilton populations (ICOHP) represent bonding and anti-bonding interactions respectively, and the zero value represents non-bonding interactions33. The more negative the ICOHP value, the stronger the bond strength. The energy-weighted ICOHP sums up to the fermi level of Cu–N and Cu–O bonds of Cu-TAPT were −1.18 and −1.09 eV, respectively, better than Cu-TABQ (Cu–N) and CuO (Cu–O). Therefore, this asymmetric coordination can effectively improve the stability of Cu–N or Cu–O, thereby improving the electrochemical stability during the CO2RR process, which is consistent with the results of in-situ XAS experiments. Subsequently, AIMD simulations were performed in the canonical ensemble (NVT) (Fig. 3f and Supplementary Fig. 38). The system energy remains basically unchanged under 300 K and electric field, indicating that the structure has excellent stability. Combined with the aforementioned in situ XAS and Raman results, as well as DFT calculations, we can conclude that the strong asymmetric coordination has ultra-high electrochemical stability and resistance to reconstruction over a broad potential window and long-term electrolysis during the CO2RR process, allowing Cu-TAPT to maintain excellent catalytic performance during the electrolysis.

Furthermore, the difference in ion transport capabilities based on Cu-TAPT and Cu-TABQ was echoed in Bode plots, obtained by in situ EIS measurements during electrically driven transport of single-component solutions in the flow cell (Supplementary Fig. 39). The characteristic frequency (f) which corresponds to the local maxima in the Bode phase plot has an inverse correlation to the overall impedance of the electrode34,35. Under the applied potential of OCP, the fCu-TAPT was higher than fCu-TABQ, indicating that lower impedance and better ion transport capability of Cu-TAPT (Supplementary Fig. 40). Likewise, as the potential increased, the fCu-TAPT gradually increased and a pronounced decrease in the low-frequency peak in the corresponding phase angle, indicating that the ion transport impedance of Cu-TAPT further decreased under the action of electric field and was better than Cu-TABQ (Supplementary Fig. 39)36. As expected, similar behaviors were also reflected in the corresponding Nyquist plots (Supplementary Fig. 41).

Reaction mechanism investigation

As is well-known that the formation of multi-carbon products undergoes the C–C coupling process and then further reduction in CO2RR37. To gain insight into the C–C coupling of multivalent Cu sites on the Cu-TAPT electrode, the Raman measurements under the potential of OCP and −0.78 V vs. RHE (optimum potential) were conducted to explore the adsorption behavior of the chemical intermediates (Supplementary Fig. 42). The observed peaks in the range 1800–2200 cm−1 can be ascribed to the *CO, including *COatop (2094 cm−1) and *CObridge (1978 cm−1) (Supplementary Fig. 43)14. Furthermore, the hydrocarbonate (HCO3−) peak at 1345 cm−1 on Cu-TAPT was almost absent in Raman spectra under −0.78 V vs. RHE (Supplementary Fig. 44)38, which indicates that the local pH value on the electrode during CO2RR is high. The high pH value of the electrode surface is believed to be beneficial to C–C coupling, lower the energy barrier for CO2 activation and inhibit the evolution of H2, and help improve the activity of CO2RR39.

To complement the above findings, we further performed in situ attenuated total reflection Fourier-transform infrared spectroscopy (ATR − FTIR) to monitor the co-adsorption behavior on the Cu-TAPT and Cu-TABQ catalyst (Supplementary Figs. 45,46). Previous studies have indicated that *CObridge protonation is energetically more favorable than *COatop40. Therefore, C–C coupling can be promoted under the asymmetry between *CO and *CHO (or *COH) after the *CObridge protonation step, and the result is also confirmed by Gibbs free energy calculations (Supplementary Fig. 47, Supplementary Table 9 and Supplementary Data 1). Notably, the adsorbed *CHO (or *COH) intermediate on Cu-TAPT was also observed from the ATR − FTIR spectrum and increased with the applied potential, strongly supporting this process (Supplementary Fig. 46). Therefore, due to the adsorption capacity difference of different valence Cu sites for CO, the energy barrier for the formation of asymmetric C–C coupling is low based on the synergistic effect of Cu2+ and Cu+, which helps to the generation of C2 products41. The peaks at 1339 and 1542 cm−1 are indexed to the absorbed *OC2H5 and *OCCOH on the catalysts, respectively (Fig. 4a, b). A ratio value of *OC2H5/*OCCOH on Cu-TAPT was larger than that on Cu-TABQ (Fig. 4c)12. This is ascribed to that the multivalent Cu sites induced unbalanced coupling environment, which accelerates the breakage of the Cu–COH bond, thus enhancing the production of ethanol intermediates, which agrees with the report that diversity in *CO binding sites enhances the formation of ethanol product19.Fig. 4 Reaction mechanism investigation.

a, b In situ ATR–FTIR spectra obtained under different potentials over Cu-TAPT a and Cu-TABQ b catalyst under CO2RR. c Potential dependence of *OC2H5/*OCCOH ratio obtained for Cu-TAPT and Cu-TABQ. d Gibbs free energy diagram for the hydrogenation (*HCCOH  to  *HCCHOH) and dehydroxylation (*HCCOH to *HCC) steps on Cu-TAPT and Cu-TABQ. e The COHP of C–OH and Cu–COH in *HCCOH over Cu-TAPT. f Charge transfer in C–OH and Cu–COH in *HCCOH adsorbed on Cu-TAPT and Cu-TABQ. All the potentials involved in this figure are not iR corrected.

Subsequently, DFT was carried out to study the intrinsic role of synergistic catalysis in promoting the selectivity of ethanol (Model information in Supplementary Data 1). Most previous mechanistic studies have suggested that formations of both ethylene and ethanol proceeds via further transformation of the *HCCOH intermediate. The *HCCOH evolves into a hydrocarbon intermediate (*HCC) or a hydroxyl-containing intermediate (*HCCHOH), serving as the formation of ethylene and ethanol. When the *HCCOH intermediate is hydrogenated and dehydrated to generate *HCC, hydrocarbons (ethylene) are generated, while when the hydroxyl group is retained to generate the *HCCHOH intermediate, ethanol is finally generated. Therefore, the distribution of ethylene and ethanol products is mainly regulated by regulating the strength of the Cu–COH and C–OH bonds on the intermediate11,15,42. We further calculated the Gibbs free energy of *HCCOH converted into *HCC and *HCCHOH to examine the impact of cooperative catalysis on C2 product distribution (Fig. 4d and Supplementary Table 9). Compared with Cu-TABQ, the synergistic effect of Cu2+ and Cu+ over Cu-TAPT can promote the hydrogenation of *HCCOH into hydroxyl-containing intermediates, reducing the reaction energy by 1.15 eV. The proton-electron coupling step of *HCCHOH is far superior to *HCC, providing an explanation for enhanced selectivity to ethanol on Cu-TAPT (Supplementary Fig. 48).

It is believed that the hydrogenation step involves the cleavage of the Cu–COH bond to form a C–H bond, while the dehydroxylation step requires the cleavage of the C–OH bond. COHP and charge density analysis were therefore used to quantify the relative strength of the Cu–COH and C–OH bonds in *HCCOH (Fig. 4e,f and Supplementary Table 10). Compared with Cu–COH, the energy-weighted ICOHP summed up to the fermi level for C–OH is more negative, and more charge transfer from C to OH group in Cu-TAPT on account of the synergistic effect of Cu2+ and Cu+, indicating that the C–OH bond is stronger than Cu–COH bond, thus promoting ethanol production.

Discussion

In summary, we have achieved highly efficient synthesis of ethanol in CO2RR by the synergy of Cu2+ and Cu+ sites on polymer catalyst Cu-TAPT, which suppresses the production of hydrocarbons and promotes the generation of ethanol. Cu-TAPT achieved a FE of 54.3 ± 3% towards ethanol at the current density of 429 mA cm−2, and the ratio of ethanol to ethylene reached 3.14:1. The coordination between Cu2+ and Cu+ with the functional groups in Cu-TAPT makes the catalyst very stable in CO2RR. Additionally, cooperative catalysis of the multivalent Cu sites can stabilize the hydroxyl-containing C2 intermediate (*HCCOH) and promotes the hydrogenation of C atoms (*HCCHOH), enhancing the selectivity to ethanol. We believe that of design of catalysts with multivalent active sites with synergistic effect can be used to prepared other efficient and selective catalysts for CO2RR.

Methods

Chemicals

Copper nitrate trihydrate (Cu(NO3)2·3H2O, 99%) was obtained from Beijing Innochem Science & Technology Co., Ltd. 2,3,5,6-Tetraaminobenzoquinone (TABQ, 97%) was purchased from Energy Chemical. Deuterated water (D2O, AR) was provided by Shanghai Meryer Biochemical Technology Co., LTD. Potassium hydroxide (KOH, AR), isopropanol (i-PrOH, AR), dimethylsulfoxide (DMSO, AR) and phenol (C6H5OH, AR) were all purchased from Sinopharm Chemical Reagent Co. Ltd. (Shanghai, China). CO2 and Ar (Shanghai Chemistry Industrial Zone Pujiang Special Type Gas Co., Ltd.) had a purity of 99.999% and were used as received. Anion exchange membrane FUMA–FAA–3–PK-130 (0.11–0.13 mm thick, 1.1–1.4 meg g−1 exchange capacity) was obtained from FuMA-Tech. Perfluorinated resin solution containing NafionTM 1100 W ( ~ 5 wt.% in lower aliphatic alcohols and water) was purchased from Sigma-Aldrich. Carbon paper-based gas diffusion layers (GDLs, Sigracet 29 BC) and nickel foam were provided by the Fuel Cell Store. Deionized water (resistivity 18.2 MΩ cm) was used for all experiments.

Synthesis of Cu-TAPT and Cu-TABQ

0.15 mmol of Cu(NO3)2·3H2O and 0.05 mmol TABQ monomer were added into 20.0 mL H2O and heated for 72 h at 120 °C. After cooling to room temperature, the precipitate was filtered and alternately washed with deionized water and ethanol. Subsequently, the Cu-TAPT was obtained by drying at 60 °C. Cu-TABQ was prepared with a similar procedure, and the only difference was that 20 mL DMSO was used as the solvent instead of water. The contents of Cu in Cu-TAPT and Cu-TABQ obtained from ICP-OES were 27.6 and 28.2 wt%, respectively.

Characterizations

XRD patterns were recorded with an X-ray diffractometer (Rigaku Ultima VI X-ray) with Cu-Kα radiation, and the scan speed was 5°/min. The morphology of catalysts was characterized by FESEM (JSM IT-800) equipped with EDS operated at 5 kV. High resolution transmission electron microscopy (HRTEM) images were obtained on Tecnai G2 F30 operated at 300 kV. HAADF-STEM images were achieved on an aberration-corrected electron projection microscope (FEI Titan G2 60-300). XPS spectra were obtained on an AXIS Supra surface analysis instrument using a monochromatic Al Kα X-ray beam (1,486.6 eV). FTIR spectra was determined utilizing an IRTracer-100 infrared spectrometer with a broad range of 400–4,000 cm−1. The Cu content was quantitatively determined by the ICP-OES (Agilent-5110). Nitrogen adsorption-desorption isotherms were measured using an automated gas sorption analyzer (AutoChem Discovery). X-ray absorption spectroscopy (XAS) measurements were carried out at the 4B9A beamline at Beijing Synchrotron Radiation Facility (BSRF) and the BL13SSW beamline at Shanghai Synchrotron Radiation Facility (SSRF), China. The Cu K-edge X-ray absorption spectra of Cu-TAPT and Cu-TABQ electrodes were obtained, and Cu foil, CuO and Cu2O and CuPc were used as references. The EXAFS oscillations χ(k) as a function of photoelectron wave number k was extracted by following standard procedures. The theoretical paths were generated using FEFF and the models were completed in a conventional way using the fitting program ARTEMIS. Part of the computational work was provided by “Ceshigo Research Service (www.ceshigo.com)” for technical support.

In situ Raman measurements

Stability measurements: In-situ Raman measurements were performed on a Renishaw inVia Reflex Raman microscope in a modified flow cell under CO2 atmosphere by using a water immersion objective lens (63x) with a 532 nm laser, with 5 s integration and averaging 20 scans per region. The modified flow-cell setup consisted of a three-electrode setup, in which platinum wire and Ag/AgCl electrodes were used as the counter and reference electrodes, respectively.

Reaction mechanism study: In situ Raman study was conducted on a Renishaw inVia Reflex Raman spectrometer with a 785 nm solid laser as an excitation source. The measurements were performed in a custom-built flow cell (Gaossunion) using graphite rods and Ag/AgCl electrodes as the anode electrode and reference electrode, respectively. A working electrode was prepared by spray brushing the catalyst onto the GDL. The electrolyte (0.5 M KOH was used as electrolytes to avoid scattering of the Raman signal by bubbles at high currents) was pumped into the GDE using a syringe pump. CO2 was passed into the back of the GDE.

In situ Attenuated Total Reflection Fourier-transform Infrared Spectroscopy (ATR–FTIR) measurements

In situ electrochemical ATR–FTIR measurements were carried out by utilizing an IRTracer-100 infrared spectrometer equipped with a liquid nitrogen-cooled MCT detector with a broad range of 1000–4000 cm−1. A hemisphere silicon prism was polished with diamond abrasion paste and sonicated with deionized water. The prim then was immersed in a piranha solution (98% H2SO4: 30% H2O2, v/v = 3/1) for 2 h. Electrocatalyst ink was consistent with electrochemical tests. 50 μL of the catalyst ink was dropped onto the prism and left to dry slowly. The above prism was mounted in a homemade electrochemical three electrode cell with an Ag/AgCl reference electrode and a platinum-wire counter electrode. The catalyst electrode was tightly pressed against the hemisphere silicon prism crystal window with a micrometre-scale gap to reduce the loss of infrared light. The electrolyte was 0.5 M KHCO3, which was constantly purged with CO2 during the experiment. The spectra under OCP were recorded for comparison.

In situ X-ray absorption spectroscopy (XAS) measurements

In situ XAS spectra were measured at the 4B9A station of the Beijing Synchrotron Radiation Facility based on CO2RR measurement conditions. A custom-built in situ electrochemical cell was used with Pt mesh counter electrode and Ag/AgCl reference electrode. The acquisition of each spectrum took around 18 minutes. Data analysis was performed using Athena and Artemis software.

XAFS analysis

XAS data were analyzed using ATHENA and ARTEMIS from the Demeter package software. The normalized EXAFS function was transformed, χ(E), in energy space to χ(k), where k is the photoelectron wave vector. χ(k) was multiplied by k3 to amplify the EXAFS oscillations in the mid-k region for assessing the interatomic interaction of Cu atoms. To differentiate the Cu K-edge EXAFS oscillation from different coordination shells, we applied Fourier transformation (from k space to R space) of the k3-weighted χ(k) with a k range of 2.5–12.4 Å−1. The Cu K edge EXAFS data were extracted from 2.5–12.4 Å−1 in k-space, and the Fourier transformed data were analyzed in the range of 1.0–2.8 Å for the prepared materials and with a k-space of 2.5–13.5 Å−1 and R-range of 1.0–3.0 Å for the reference materials. The experimental data were fitted with theoretical backscattering paths calculated using FEFF6, based upon CIF files from the Materials Project. The fits were performed based on k3 weighting. The amplitude reduction factor, S02, was determined by fitting a reference copper foil, using the first single-scattering paths with tabulated coordination numbers of bulk copper (fcc–Cu metal: mp–30, symmetrized). The value of S02 was determined to be 0.82. For Wavelet Transform analysis, the χ(k) exported from Athena was performed to provide the overall distribution.

Electrocatalytic CO2 reduction reaction (CO2RR) measurements

All the electrocatalytic experiments were performed in a flow cell (101017, Gaossunion) on the electrochemical workstation (CHI 660E, Shanghai CH Instruments Co., China) at ambient pressure and temperature. Typically, 6 mg of catalyst and 60 μL of 5 wt% Nafion 1100 W were dispersed in 1 mL i-PrOH and sonicated for 20 min. Then, the resulting ink was dropped onto a piece of 1.5 cm × 2.5 cm GDL to prepare gas diffusion electrodes (GDEs). The prepared electrodes were dried under an infrared lamp around ambient environment for 10 min. In the flow cell, the prepared GDEs, Ag/AgCl electrode (saturated KCl) and nickel foam electrode were employed as the working, reference and counter electrodes, respectively. The electrolyte was 1 M KOH aqueous solution. Anode chamber and cathode chamber were separated by an anionic membrane. 40 mL 1 M KOH was circulated in cathode chamber and anode chamber at a constant rate of 20 mL min−1 by peristaltic pump. The gaseous CO2 was flowed by the backside of the GDE at a rate of 20 SCCM, controlled by a mass flow meter. The geometric surface area of the catalyst was 1 cm2. All the potentials were reported with respect to the RHE using the formula14,25:1 E(vs.RHE)=E(vs.Ag/AgCl)+0.197V+0.059×pH+IR×70%

The ohmic losses was used to measure between the working and reference electrodes via EIS at the beginning of the electrolysis, with 70% IR compensation (i, current; R, uncompensated resistance) for manual correction of the potential.

During the stability test, the violent fluctuation of current was ascribed to the device itself, which would cause electrolyte flooding and carbonate precipitation in strong alkali environment to block the gas channel. Therefore, the flow cell needed to be disassembled and reassembled. Also, the electrolyte was also replaced in total circulation volume at each orange arrow, leading to the recovery of the catalytic performance.

Electrochemical measurements

All the electrocatalytic tests were conducted at ambient pressure and temperature in a flow cell. LSV curves were conducted with scan rate of 100 mV s−1. The Cdl value was obtained by recording cyclic voltammetry (CV) at various scan rates in the non-Faradaic region by flow cell.2 ECSA=Cdl/Cs

where Cdl and Cs are the double layer capacitance and specific capacitance, respectively.

EIS measurement was carried out on the electrochemical analyzer in a frequency range from 105 Hz to 10−2 Hz with an amplitude of 5 mV at an open circuit potential and equilibrium potential. The frequency and overall impedance as related by the following Equation35:3 Z=1/(2πf)n

The data obtained from the EIS measurements were fitted by the Zview software (Version 3.1, Scribner Associates, USA). From EIS data, it is straightforward to evaluate the electrical parameters. The intercept corresponding to the real impedance axis gives the Rs (also electrochemical cell resistance) and then the ionic conductivity can be obtained using the following Equation36:4 σion=L/(Rs*S)

where L and S are the sample thickness and contact area, respectively.

CO2RR product analysis

The gas products were analyzed by the GC (Agilent-8890), which was equipped with the flame ionization detectors (FID) and thermal conductivity detectors (TCD). The liquid products were quantified by NMR spectrometer (Bruker; Ascend 500 MHz) in D2O with C6H5OH (200 mM) and DMSO (5 mM) as internal standards (Supplementary Figs. 16,19). Solvent pre-saturation technique was implemented to suppress the water peak. The moles of a liquid product were calculated from integral areas and calibration curves. After the quantification, the FE of each product was calculated as follows:5 FE=(n×F×C×V)/Q×100%

Where Q is the amount of charge passed through the working electrode; F is the Faraday constant (96,485 C mol−1); C is the concentration of product; V is the volume of the electrolyte; n is the number of the electrons transfer for 1 mole product in electrocatalysis.

The calculation of TOF:6 TOF=(Iproduct×MCu)/(n×F×mcat×ω)×3600

Where Iproduct is the partial current for certain product (A); n is the number of electrons transferred for product formation (which is 12 for ethanol production.); mcat is the catalyst mass in the electrode (g); ω is in loading in the catalyst. MCu is the atomic mass of Cu, 63.55 g mol−1.

Computational details

DFT calculations were performed on the Dmol3 module by the Materials Studio 202043,44. The structures of intermediates in electrocatalysis were optimized, followed by the calculation of the frequency. The generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) function and TS for DFT-D correction were employed in the calculations45. The convergence tolerance of energy, force and displacement convergence were set as 1 × 10−5 Ha, 2 × 10−3 Ha/Å and 5 × 10−3 Å, respectively. The bond configuration analysis is based on the COHP theory46. The energy contribution of electrons filling the electron orbitals of a crystal is divided into the energy contribution of the overlapping atomic/fragment (molecule) orbital pairs. The time step for the AIMD simulation is 2 fs, and the canonical ensemble condition (NVT) and the electric field (Z-axis) was imposed with a target temperature of 300 K47. The smearing function was employed to improve the convergence of states near the Fermi level with a smearing width of 0.005 eV.

The free energy of a species was calculated according to the standard formula:7 ΔG=ΔE+ΔZPE−TΔS

where ΔE is the DFT total energy, ΔZPE stands for the correction of zero-point energy, T denotes the thermodynamic temperature (298.15 K), and ΔS is the variation of entropy. ΔZPE and ΔS are calculated with the vibrational frequencies of adsorbates and molecules as calculated within DFT.

Supplementary information

Supplementary Information

Peer Review File

Supplementary Dataset 1

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51928-z.

Acknowledgements

Han acknowledges the National Key R&D Program of China (2023YFA1507901, 2020YFA0710201), the National Natural Science Foundation of China (22293015, 22121002), and the Fundamental Research Funds for the Central Universities. Jia acknowledges the China Postdoctoral Science Foundation (2023M731096). We thank the Beijing Synchrotron Radiation Facility (BSRF, 4B9A beamline) and the Shanghai Synchrotron Radiation Facility (SSRF, BL13SSW beamline) for providing the XAS measurements.

Author contributions

X.C., S.Q.J., H.H.W. and B.X.H. proposed the project, designed the experiments, and wrote the manuscript. X.C. performed most of the experiments and analyzed the experimental data. X.C., J.X.Z. and J.P.J. carried out DFT calculations and analyzed the computational data. M.K.D., C.X., T.D. and H.L.C. provided help in materials synthesis and characterizations. X.C., Z.H.X., C.J.C., S.Q.J., X.Q.X and J.R.Z conducted XANES measurements and analyzed the results. H.H.W., M.Y.H. and B.X.H. co-supervised the whole project. All authors discussed the results and commented on the manuscript.

Peer review

Peer review information

Nature Communications thanks Jin Hyuk Cho and the other, anonymous, reviewer for their contribution to the peer review of this work. A peer review file is available.

Data availability

Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Birdja YY Advances and challenges in understanding the electrocatalytic conversion of carbon dioxide to fuels Nat. Energy 2019 4 732 745 10.1038/s41560-019-0450-y
Birdja, Y. Y. et al. Advances and challenges in understanding the electrocatalytic conversion of carbon dioxide to fuels. Nat. Energy 4, 732–745 (2019).10.1038/s41560-019-0450-y
2. Lees EW Mowbray BAW Parlane FGL Berlinguette CP Gas diffusion electrodes and membranes for CO2 reduction electrolysers Nat. Rev. Mater. 2021 7 55 64 10.1038/s41578-021-00356-2
Lees, E. W., Mowbray, B. A. W., Parlane, F. G. L. & Berlinguette, C. P. Gas diffusion electrodes and membranes for CO2 reduction electrolysers. Nat. Rev. Mater. 7, 55–64 (2021).10.1038/s41578-021-00356-2
3. Ager JW Lapkin AA Chemical storage of renewable energy Science 2018 360 707 708 10.1126/science.aat7918 29773731
Ager, J. W. & Lapkin, A. A. Chemical storage of renewable energy. Science 360, 707–708 (2018).29773731 10.1126/science.aat7918
4. Jouny M Luc W Jiao F General techno-economic analysis of CO2 electrolysis systems Ind. Eng. Chem. Res. 2018 57 2165 2177 10.1021/acs.iecr.7b03514
Jouny, M., Luc, W. & Jiao, F. General techno-economic analysis of CO2 electrolysis systems. Ind. Eng. Chem. Res. 57, 2165–2177 (2018).10.1021/acs.iecr.7b03514
5. Zhang B Highly electrocatalytic ethylene production from CO2 on nanodefective Cu nanosheets J. Am. Chem. Soc. 2020 142 13606 13613 10.1021/jacs.0c06420 32658474
Zhang, B. et al. Highly electrocatalytic ethylene production from CO2 on nanodefective Cu nanosheets. J. Am. Chem. Soc. 142, 13606–13613 (2020).32658474 10.1021/jacs.0c06420
6. Chen Y Efficient multicarbon formation in acidic CO2 reduction via tandem electrocatalysis Nat. Nanotechnol. 2023 19 311 318 10.1038/s41565-023-01543-8 37996517
Chen, Y. et al. Efficient multicarbon formation in acidic CO2 reduction via tandem electrocatalysis. Nat. Nanotechnol. 19, 311–318 (2023).37996517 10.1038/s41565-023-01543-8
7. Spurgeon JM Kumar B A comparative technoeconomic analysis of pathways for commercial electrochemical CO2 reduction to liquid products Energy Environ. Sci. 2018 11 1536 1551 10.1039/C8EE00097B
Spurgeon, J. M. & Kumar, B. A comparative technoeconomic analysis of pathways for commercial electrochemical CO2 reduction to liquid products. Energy Environ. Sci. 11, 1536–1551 (2018).10.1039/C8EE00097B
8. Calle‐Vallejo F Koper MTM Theoretical considerations on the electroreduction of CO to C2 species on Cu(100) electrodes Angew. Chem. Int. Ed. 2013 52 7282 7285 10.1002/anie.201301470
Calle‐Vallejo, F. & Koper, M. T. M. Theoretical considerations on the electroreduction of CO to C2 species on Cu(100) electrodes. Angew. Chem. Int. Ed. 52, 7282–7285 (2013).10.1002/anie.201301470
9. Nguyen TN Electrochemical CO2 reduction to ethanol: from mechanistic understanding to catalyst design J. Mater. Chem. A 2021 9 12474 12494 10.1039/D1TA01115D
Nguyen, T. N. et al. Electrochemical CO2 reduction to ethanol: from mechanistic understanding to catalyst design. J. Mater. Chem. A 9, 12474–12494 (2021).10.1039/D1TA01115D
10. Zang Y Selective CO2 electroreduction to ethanol over a carbon‐coated CuOx catalyst Angew. Chem. Int. Ed. 2022 61 e202209629 10.1002/anie.202209629
Zang, Y. et al. Selective CO2 electroreduction to ethanol over a carbon‐coated CuOx catalyst. Angew. Chem. Int. Ed. 61, e202209629 (2022).10.1002/anie.202209629
11. Xu A Copper/alkaline earth metal oxide interfaces for electrochemical CO2-to-alcohol conversion by selective hydrogenation Nat. Catal. 2022 5 1081 1088 10.1038/s41929-022-00880-6
Xu, A. et al. Copper/alkaline earth metal oxide interfaces for electrochemical CO2-to-alcohol conversion by selective hydrogenation. Nat. Catal. 5, 1081–1088 (2022).10.1038/s41929-022-00880-6
12. Zhang L Oxophilicity-controlled CO2 electroreduction to C2+ alcohols over lewis acid metal-doped Cuδ+ catalysts J. Am. Chem. Soc. 2023 145 21945 21954 10.1021/jacs.3c06697 37751566
Zhang, L. et al. Oxophilicity-controlled CO2 electroreduction to C2+ alcohols over lewis acid metal-doped Cuδ+ catalysts. J. Am. Chem. Soc. 145, 21945–21954 (2023).37751566 10.1021/jacs.3c06697
13. Peng C Surface Co‐modification of halide anions and potassium cations promotes high‐rate CO2‐to‐ethanol electrosynthesis Adv. Mater. 2022 34 2204476 10.1002/adma.202204476
Peng, C. et al. Surface Co‐modification of halide anions and potassium cations promotes high‐rate CO2‐to‐ethanol electrosynthesis. Adv. Mater. 34, 2204476 (2022).10.1002/adma.202204476
14. Wang P Boosting electrocatalytic CO2-to-ethanol production via asymmetric C–C coupling Nat. Commun. 2022 13 3754 10.1038/s41467-022-31427-9 35768462
Wang, P. et al. Boosting electrocatalytic CO2-to-ethanol production via asymmetric C–C coupling. Nat. Commun. 13, 3754 (2022).35768462 10.1038/s41467-022-31427-9
15. Wang X Efficient electrically powered CO2-to-ethanol via suppression of deoxygenation Nat. Energy 2020 5 478 486 10.1038/s41560-020-0607-8
Wang, X. et al. Efficient electrically powered CO2-to-ethanol via suppression of deoxygenation. Nat. Energy 5, 478–486 (2020).10.1038/s41560-020-0607-8
16. Feng J Improving CO2-to-C2+ product electroreduction efficiency via atomic lanthanide dopant-induced tensile-strained CuOx catalysts J. Am. Chem. Soc. 2023 145 9857 9866 10.1021/jacs.3c02428 37092347
Feng, J. et al. Improving CO2-to-C2+ product electroreduction efficiency via atomic lanthanide dopant-induced tensile-strained CuOx catalysts. J. Am. Chem. Soc. 145, 9857–9866 (2023).37092347 10.1021/jacs.3c02428
17. Zhang T Yuan B Wang W He J Xiang X Tailoring *H intermediate coverage on the CuAl2O4/CuO catalyst for enhanced electrocatalytic CO2 reduction to ethanol Angew. Chem. Int. Ed. 2023 62 e202302096 10.1002/anie.202302096
Zhang, T., Yuan, B., Wang, W., He, J. & Xiang, X. Tailoring *H intermediate coverage on the CuAl2O4/CuO catalyst for enhanced electrocatalytic CO2 reduction to ethanol. Angew. Chem. Int. Ed. 62, e202302096 (2023).10.1002/anie.202302096
18. Sun W V‐doped Cu2Se hierarchical nanotubes enabling flow‐cell CO2 electroreduction to ethanol with high efficiency and selectivity Adv. Mater. 2022 34 2207691 10.1002/adma.202207691
Sun, W. et al. V‐doped Cu2Se hierarchical nanotubes enabling flow‐cell CO2 electroreduction to ethanol with high efficiency and selectivity. Adv. Mater. 34, 2207691 (2022).10.1002/adma.202207691
19. Li YC Binding site diversity promotes CO2 electroreduction to ethanol J. Am. Chem. Soc. 2019 141 8584 8591 10.1021/jacs.9b02945 31067857
Li, Y. C. et al. Binding site diversity promotes CO2 electroreduction to ethanol. J. Am. Chem. Soc. 141, 8584–8591 (2019).31067857 10.1021/jacs.9b02945
20. Luo M Hydroxide promotes carbon dioxide electroreduction to ethanol on copper via tuning of adsorbed hydrogen Nat. Commun. 2019 10 5814 10.1038/s41467-019-13833-8 31862886
Luo, M. et al. Hydroxide promotes carbon dioxide electroreduction to ethanol on copper via tuning of adsorbed hydrogen. Nat. Commun. 10, 5814 (2019).31862886 10.1038/s41467-019-13833-8
21. Zhang R Molecular engineering of a metal-organic polymer for enhanced electrochemical nitrate-to-ammonia conversion and zinc nitrate batteries Angew. Chem. Int. Ed. 2023 62 e202309930 10.1002/anie.202309930
Zhang, R. et al. Molecular engineering of a metal-organic polymer for enhanced electrochemical nitrate-to-ammonia conversion and zinc nitrate batteries. Angew. Chem. Int. Ed. 62, e202309930 (2023).10.1002/anie.202309930
22. Luo M Coordination polymer electrocatalysts enable efficient CO-to-acetate conversion Adv. Mater. 2023 35 e2209567 10.1002/adma.202209567 36584285
Luo, M. et al. Coordination polymer electrocatalysts enable efficient CO-to-acetate conversion. Adv. Mater. 35, e2209567 (2023).36584285 10.1002/adma.202209567
23. Fan K Single crystals of a highly conductive three-dimensional conjugated coordination polymer J. Am. Chem. Soc. 2023 145 12682 12690 10.1021/jacs.3c02378 37204114
Fan, K. et al. Single crystals of a highly conductive three-dimensional conjugated coordination polymer. J. Am. Chem. Soc. 145, 12682–12690 (2023).37204114 10.1021/jacs.3c02378
24. Chen X Highly stable layered coordination polymer electrocatalyst towards efficient CO2‐to‐CH4 conversion Adv. Mater. 2024 36 2310273 10.1002/adma.202310273
Chen, X. et al. Highly stable layered coordination polymer electrocatalyst towards efficient CO2‐to‐CH4 conversion. Adv. Mater. 36, 2310273 (2024).10.1002/adma.202310273
25. Liang Y Stabilizing copper sites in coordination polymers toward efficient electrochemical C–C coupling Nat. Commun. 2023 14 474 10.1038/s41467-023-35993-4 36710270
Liang, Y. et al. Stabilizing copper sites in coordination polymers toward efficient electrochemical C–C coupling. Nat. Commun. 14, 474 (2023).36710270 10.1038/s41467-023-35993-4
26. Liao L MoS2 formed on mesoporous graphene as a highly active catalyst for hydrogen evolution Adv. Funct. Mater. 2013 23 5326 5333 10.1002/adfm.201300318
Liao, L. et al. MoS2 formed on mesoporous graphene as a highly active catalyst for hydrogen evolution. Adv. Funct. Mater. 23, 5326–5333 (2013).10.1002/adfm.201300318
27. Watanabe M Application of ionic liquids to energy storage and conversion materials and devices Chem. Rev. 2017 117 7190 7239 10.1021/acs.chemrev.6b00504 28084733
Watanabe, M. et al. Application of ionic liquids to energy storage and conversion materials and devices. Chem. Rev. 117, 7190–7239 (2017).28084733 10.1021/acs.chemrev.6b00504
28. Tan X Stabilizing copper by a reconstruction-resistant atomic Cu–O–Si interface for electrochemical CO2 reduction J. Am. Chem. Soc. 2023 145 8656 8664 10.1021/jacs.3c01638
Tan, X. et al. Stabilizing copper by a reconstruction-resistant atomic Cu–O–Si interface for electrochemical CO2 reduction. J. Am. Chem. Soc. 145, 8656–8664 (2023).10.1021/jacs.3c01638
29. Liu H Asymmetric coordination of single‐atom Co sites achieves efficient dehydrogenation catalysis Adv. Funct. Mater. 2022 32 2207408 10.1002/adfm.202207408
Liu, H. et al. Asymmetric coordination of single‐atom Co sites achieves efficient dehydrogenation catalysis. Adv. Funct. Mater. 32, 2207408 (2022).10.1002/adfm.202207408
30. Sun K Co(CN)3 catalysts with well-defined coordination structure for the oxygen reduction reaction Nat. Catal. 2023 6 1164 1173 10.1038/s41929-023-01047-7
Sun, K. et al. Co(CN)3 catalysts with well-defined coordination structure for the oxygen reduction reaction. Nat. Catal. 6, 1164–1173 (2023).10.1038/s41929-023-01047-7
31. Huang S Well-defined N3C1-anchored single-metal-sites for oxygen reduction reaction Angew. Chem. Int. Ed. 2024 63 e202314833 10.1002/anie.202314833
Huang, S. et al. Well-defined N3C1-anchored single-metal-sites for oxygen reduction reaction. Angew. Chem. Int. Ed. 63, e202314833 (2024).10.1002/anie.202314833
32. Song P Modulating the asymmetric atomic interface of copper single atoms for efficient CO2 electroreduction ACS Nano 2023 17 4619 4628 10.1021/acsnano.2c10701 36815694
Song, P. et al. Modulating the asymmetric atomic interface of copper single atoms for efficient CO2 electroreduction. ACS Nano 17, 4619–4628 (2023).36815694 10.1021/acsnano.2c10701
33. Duan X Stabilizing single-atomic ruthenium by ferrous ion doped NiFe-LDH towards highly efficient and sustained water oxidation Chem. Eng. J. 2022 446 136962 10.1016/j.cej.2022.136962
Duan, X. et al. Stabilizing single-atomic ruthenium by ferrous ion doped NiFe-LDH towards highly efficient and sustained water oxidation. Chem. Eng. J. 446, 136962 (2022).10.1016/j.cej.2022.136962
34. Mani AM Kumar A Chaudhury S Decoding transport selectivity of ions in polymer membranes by in-situ impedance spectroscopy Sep. Sci. Technol. 2023 58 2411 2421 10.1080/01496395.2023.2219377
Mani, A. M., Kumar, A. & Chaudhury, S. Decoding transport selectivity of ions in polymer membranes by in-situ impedance spectroscopy. Sep. Sci. Technol. 58, 2411–2421 (2023).10.1080/01496395.2023.2219377
35. Abdalrahman AA Aziz SB Karim WO EIS and FTIR approaches to study the ion transport parameters and relaxation dynamics of Na+1 ion in SPE based on MC polymer inserted with sodium salt Results Phys. 2022 36 105439 10.1016/j.rinp.2022.105439
Abdalrahman, A. A., Aziz, S. B. & Karim, W. O. EIS and FTIR approaches to study the ion transport parameters and relaxation dynamics of Na+1 ion in SPE based on MC polymer inserted with sodium salt. Results Phys. 36, 105439 (2022).10.1016/j.rinp.2022.105439
36. Kim HJ Chen B Suo Z Hayward RC Ionoelastomer junctions between polymer networks of fixed anions and cations Science 2020 367 773 776 10.1126/science.aay8467 32054759
Kim, H. J., Chen, B., Suo, Z. & Hayward, R. C. Ionoelastomer junctions between polymer networks of fixed anions and cations. Science 367, 773–776 (2020).32054759 10.1126/science.aay8467
37. Dinh C-T CO2 electroreduction to ethylene via hydroxide-mediated copper catalysis at an abrupt interface Science 2018 360 783 787 10.1126/science.aas9100 29773749
Dinh, C.-T. et al. CO2 electroreduction to ethylene via hydroxide-mediated copper catalysis at an abrupt interface. Science 360, 783–787 (2018).29773749 10.1126/science.aas9100
38. Li H Wei P Gao D Wang G In situ Raman spectroscopy studies for electrochemical CO2 reduction over Cu catalysts Curr. Opin. Green Sustain. Chem. 2022 34 100589 10.1016/j.cogsc.2022.100589
Li, H., Wei, P., Gao, D. & Wang, G. In situ Raman spectroscopy studies for electrochemical CO2 reduction over Cu catalysts. Curr. Opin. Green Sustain. Chem. 34, 100589 (2022).10.1016/j.cogsc.2022.100589
39. Chen X Electrochemical CO2-to-ethylene conversion on polyamine-incorporated Cu electrodes Nat. Catal. 2020 4 20 27 10.1038/s41929-020-00547-0
Chen, X. et al. Electrochemical CO2-to-ethylene conversion on polyamine-incorporated Cu electrodes. Nat. Catal. 4, 20–27 (2020).10.1038/s41929-020-00547-0
40. Shen H Asymmetrical C–C Coupling for electroreduction of CO on bimetallic Cu–Pd catalysts ACS Catal. 2022 12 5275 5283 10.1021/acscatal.2c00646
Shen, H. et al. Asymmetrical C–C Coupling for electroreduction of CO on bimetallic Cu–Pd catalysts. ACS Catal. 12, 5275–5283 (2022).10.1021/acscatal.2c00646
41. Ma W Electrocatalytic reduction of CO2 to ethylene and ethanol through hydrogen-assisted C–C coupling over fluorine-modified copper Nat. Catal. 2020 3 478 487 10.1038/s41929-020-0450-0
Ma, W. et al. Electrocatalytic reduction of CO2 to ethylene and ethanol through hydrogen-assisted C–C coupling over fluorine-modified copper. Nat. Catal. 3, 478–487 (2020).10.1038/s41929-020-0450-0
42. Zhang J Steering CO2 electroreduction pathway toward ethanol via surface-bounded hydroxyl species-induced noncovalent interaction Proc. Natl. Acad. Sci. USA 2023 120 e2218987120 10.1073/pnas.2218987120 36877842
Zhang, J. et al. Steering CO2 electroreduction pathway toward ethanol via surface-bounded hydroxyl species-induced noncovalent interaction. Proc. Natl. Acad. Sci. USA 120, e2218987120 (2023).36877842 10.1073/pnas.2218987120
43. Delley B An all-electron numerical method for solving the local density functional for polyatomic molecules J. Chem. Phys. 1990 92 508 517 10.1063/1.458452
Delley, B. An all-electron numerical method for solving the local density functional for polyatomic molecules. J. Chem. Phys. 92, 508–517 (1990).10.1063/1.458452
44. Delley B From molecules to solids with the DMol3 approach J. Chem. Phys. 2000 113 7756 7764 10.1063/1.1316015
Delley, B. From molecules to solids with the DMol3 approach. J. Chem. Phys. 113, 7756–7764 (2000).10.1063/1.1316015
45. Perdew JP Burke K Ernzerhof M Generalized gradient approximation made simple Phys. Rev. Lett. 1996 77 3865 3868 10.1103/PhysRevLett.77.3865 10062328
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).10062328 10.1103/PhysRevLett.77.3865
46. Dronskowski R Bloechl PE Crystal orbital Hamilton populations (COHP): energy-resolved visualization of chemical bonding in solids based on density-functional calculations J. Phys. Chem. 1993 97 8617 8624 10.1021/j100135a014
Dronskowski, R. & Bloechl, P. E. Crystal orbital Hamilton populations (COHP): energy-resolved visualization of chemical bonding in solids based on density-functional calculations. J. Phys. Chem. 97, 8617–8624 (1993).10.1021/j100135a014
47. Fischer SA Duncan WR Prezhdo OV Ab initio nonadiabatic molecular dynamics of wet-electrons on the TiO2 surface J. Am. Chem. Soc. 2009 131 15483 15491 10.1021/ja906599b 19780540
Fischer, S. A., Duncan, W. R. & Prezhdo, O. V. Ab initio nonadiabatic molecular dynamics of wet-electrons on the TiO2 surface. J. Am. Chem. Soc. 131, 15483–15491 (2009).19780540 10.1021/ja906599b
