Theoretical insights into selective electrochemical conversion of carbon dioxide

Electrochemical conversion of CO2 and water to valuable chemicals and fuels is one of the promising alternatives to replace fossil fuel-based processes in realizing a carbon–neutral cycle. For practical application of such technologies, suppressing hydrogen evolution reaction and facilitating the activation of stable CO2 molecules still remain major challenges. Furthermore, high production selectivity toward high-value chemicals such as ethylene, ethanol, and even n-propanol is also not easy task to achieve. To settle these challenges, deeper understanding on underlying basis of reactions such as how intermediate binding affinities can be engineered at catalyst surfaces need to be discussed. In this review, we briefly outline recent strategies to modulate the binding energies of key intermediates for CO2 reduction reactions, based on theoretical insights from density functional theory calculation studies. In addition, important design principles of catalysts and electrolytes are also provided, which would contribute to the development of highly active catalysts for CO2 electroreduction.


Introduction
Climate change due to greenhouse gases has been a critical issue that must be solved for human being's sustainable life [1]. Indiscriminate exploitation of fossil fuel has accelerated the global warming issue associated with greenhouse gas emission [2], therefore, clean energy generation technologies should be developed to mostly or partly replace the fossil fuel-based processes. Electrochemically converting CO 2 and water to valuable chemicals and fuels, using renewable energy resources such as solar energy, has been one of the promising strategies to realize carbon-neutral energy cycle [3,4]. It can generate a variety of products including HCOOH, CO, CH 4 , C 2 H 4 , C 2 H 5 OH, and even C 3 -C 4 chemicals under ambient conditions in a sustainable fashion [5][6][7][8][9]. However, for practical implementation of the electrochemical CO 2 conversion technology, suppressing hydrogen evolution reactions (HER) and enhancing selectivity for specific products remain challenges. In general, the HER, a proton reduction reaction, easily occurs as a competing reaction of CO 2 reduction reactions (CO 2 RR) because it is kinetically more facile [10]. In addition, the CO 2 can be reduced to various intermediates and products via multiple reaction pathways, resulting in poor selectivity of the reaction [6].
In the electrochemical reaction system, product selectivity of CO 2 RR is deeply associated with the binding affinity of reaction intermediates on the catalyst surface. In specific, polycrystalline Au, Ag and Zn electrodes can evolve CO with high selectivity because these metals intrinsically possess low binding affinities for an adsorbed *CO intermediate [10][11][12][13][14][15]. In contrast, Pt, Ni, and Fe electrodes are poisoned by *CO owing to their strong binding affinities for *CO [16]. As a result, HER proceeds with high selectivity over CO 2 RR. Meanwhile, high selectivity for HCOO − production have been shown from the metals like Sn, Pb and Bi that relatively favor *OCHO over *COOH [17][18][19]. Therefore, engineering or controlling binding affinities of key intermediates is a critical issue in the product selectivity enhancement.
Recently, effective strategies to tune intermediate binding affinities have been suggested and demonstrated by several pioneering works. Representatively, it was reported that electrochemically reduced Cu oxides possess subsurface oxygens near the catalyst surface, which can increase the *CO binding affinity because the oxygens withdraw electrons from the Cu sp-band and reduce the σ-repulsion between a Cu atom and *CO molecule [20]. As a consequence, C-C coupling reactions by CO dimerization at the catalyst surface with high CO coverage can be kinetically promoted. Such a stronger CO binding through subsurface oxygens was also theoretically rationalized in Ag case [21]. In addition, it has been suggested that nanoneedle morphology can enhance electric field at the tip of nanostructure, thereby the alkali metal ions can be concentrated [22]. The solvated alkali metal cations selectively stabilize the polarizable intermediates like *CO 2 , *CO, *OCCO and *OCCHO due to the electrostatic interactions while *H is insensitive to the cations due to negligible polarizability [23]. As a result, higher CO 2 RR selectivity could be achieved.
Here, we first address intrinsic limitations of conventional pure metals in terms of binding affinity and product selectivity. Then, we review recent strategies to allow for the binding affinity modulation of CO 2 reduction intermediates. The strategies are classified into five groups: (1) oxygen incorporation, (2) electrolyte engineering, (3) nanostructuring, (4) introducing surface ligands, and (5) single-atom catalyst. In each subject, we cover theoretical insights into how the modifications can induce the tuning of binding affinity based on density functional theory (DFT) calculation studies. This review provides important design principles of electrocatalysts to promote CO 2 RR over HER, as well as C-C coupling via CO dimerization pathway.

Challenges of electrochemical CO reduction reactions
CO 2 is a stable molecule (ΔG f = − 394.4 kJ/mol) that has electrophilic carbon center due to neighboring oxygen atoms with a linear configuration. The one-electron transfer to CO 2 to form the CO 2 − anion, an initially created intermediate in most cases, requires a very negative potential of − 1.9 V vs. SHE in a neutral aqueous solution (Table 1) [24]. This thermodynamic limitation makes the overpotential for CO 2 RR much larger than HER, leading to low selectivity for CO 2 RR. Furthermore, although the CO 2 molecule can be activated at more favorable overpotentials by the concerted proton-electron transfer pathway or the promotional effect of alkali metal cations contained in electrolytes [25,26], a criticial issue is laid on controlling product selectivity. The multiple CO 2 reduction reactions that produce different chemicals including HCOO − , CO, CH 4 , C 2 H 4 and C 2 H 5 OH can occur concurrently as described in Table 1. In addition, HER acts as a major competing reaction. The HER is kinetically more facile than CO 2 RR because the number of electrons involved in the reaction is only two, and because the concentrations of proton sources such as water and bicarbonate ions are much higher than that of CO 2 in the neutral or basic aqueous solutions frequently used in CO 2 RR experiments [27]. Due to the thermodynamic as well as kinetic reasons, selectively converting CO 2 to desirable products at low overpotential is challenging. To overcome the issues, we need to deeply understand crucial factors that affect the CO 2 RR properties.
Previously, Hori et al. [16] have shown that the product selectivity changes with the type of transition metal  [28][29][30]. Norskov group explained that Au and Ag metals, representative catalysts for CO production, have weak *CO binding energy, which allows for *CO desorption rather than further hydrogenation toward *CHO and *COH as shown in Fig. 1 [30]. Additionally, the Au and Ag metal can show a low catalytic activity for HER because of the weak binding energy of *H. In contrast, the Ni, Pt, Pd and Rh metals which dominantly produce H 2 , have too strong *CO binding energy to evolve *CO, whereas HER is facilitated due to their ideal *H binding energy based on Sabatier principle. Exceptionally, the Cu metal has a medium *CO binding energy and relatively unfavorable limiting potential for HER. Thereby, further hydrogenation of *CO intermediates can proceed at the catalyst surface, generating various hydrocarbons and alcohols.
While the binding energy of the key reaction intermediates can give us important insights as a descriptor to understand intrinsic activity of various metals, enhancing their catalytic properties such as selectivity and overpotential is limited by scaling relations between the binding energies of various intermediates. For instance, the binding energies for *COOH and *CO intermediates are correlated on transition metals by the equation, ΔG CO = 1.31 × ∆G COOH − 0.47 [5]. Thus, the binding energies are hard to be controlled independently although stronger binding of *COOH and weaker binding of *CO is required to efficiently promote CO production. In addition, the scaling relation also exists between  [5]. which also makes it difficult to selectively suppress HER. In this regard, breaking these scaling relations is critical in enhancing the catalyst properties for CO 2 RR. In the following chapter, we outline recent strategies that allow for the independent modulation of the binding energy for specific intermediate, in terms of catalyst development and eletrolyte engineering.

Oxygen incorporation
Electrochemically reduced copper oxides, i.e. oxidederived Cu (OD-Cu), have been intensively studied as an electrocatalyst for the production of C 2+ chemicals from CO 2 RR [20,24,31], since Li et al. [32] first reported that the OD-Cu electrode can promote C-C coupling reactions with the suppression of HER in CO reduction reactions (CORR). It is known that the C-C coupling reactions in CO 2 RR actually proceed via similar pathways with those of CORR, such as CO dimerization [33,34], therefore the OD-Cu electrode can show the enhanced C-C coupling capability even in CO 2 RR. The promotional effect of the OD-Cu electrode has been explained based on several hypotheses. Feng et al. [35] demonstrated that surface density of grain boundary is correlated with the catalytic activity of CO reduction, suggesting that defective grain boundary formed by oxidation treatment is an active site for C-C coupling reactions. Meanwhile, Xiao et al. [36] showed from DFT calculations that the Cu 2 O-Cu interface can lower the energetic barrier of C-C bond formation by attractive electrostatics between two carbon atoms bound to Cu + and Cu 0 atoms, respectively. Furthermore, it has been proposed that subsurface oxygen (O sub ) plays an important role in modulating the binding strength of CO molecules which are known to be a key intermediate for C-C coupling reaction [20]. Here, we explain the promotional effect of OD-Cu electrodes focusing on the subsurface oxygen.
Eilert et al. [20] performed in situ ambient pressure X-ray photoelectron spectroscopy (APXPS) analysis to identify that oxygen species can be stably maintained at the surface even after applying cathodic potentials enough to reduce Cu oxides. Figure 2a shows the O 1 s APXPS spectra of the initial, oxidized, and reduced Cu foil. After oxidation, much thicker water overlayer was formed and oxidized compounds assignable to CuCO 3 , Cu(OH) 2 and Cu 2 O were created. After re-reduction, they found that an oxygen peak increased much more in comparison with the initial spectrum of Cu foil. They hypothesized that the enhanced oxygen peak can be attributed to the formation of O sub , based on the peak position and additional oxygen K-edge electron energyloss spectra (EELS) analysis.
In order to elucidate the influence of O sub in CO 2 RR properties, the CO binding energy was calculated on clean Cu (100) and on Cu (100) with subsurface  [20], copyright 2017 American Chemical Society. c Reproduced with permission from [37], copyright 2017 National Academy of Sciences) oxygens in different layers below surface. The CO binding became stronger in the presence of subsurface oxygens. Such enhancement of CO binding was explained by local charge interaction between O sub and a Cu atom as shown in Fig. 2b. The O sub could withdraw electrons from Cu sp-band, and lower the σ-repulsion between a Cu atom and CO molecule. Ultimately, it was suggested that stronger CO binding leads to a higher CO coverage at the catalyst surface, and kinetically promotes C-C coupling reactions over hydrogenation reactions such as CH 4

and H 2 evolutions.
It is also reported that the O sub can lower the activation barrier for the transition from a physisorbed linear CO 2 (l-CO 2 ) to a bent chemisorbed CO 2 (b-CO 2 ). Favaro et al. [37] predicted the binding structures of chemisorbed CO 2 and H 2 O on Cu (111) with different levels of O sub (Fig. 2c). Also, they identified that the free energies of the chemisorbed systems depend on O sub levels. When 0.25 mono-layer (ML) O sub was incorporated, the free energy was much lowered in comparison with unmodified Cu (111). The presence of O sub generated a heterovalent surface composed of Cu 0 and Cu + atoms, where the Cu 0 and Cu + centers were synergistically bound to the C and O atoms of the b-CO 2 , respectively. The electronic interactions between the binding centers and the adsorbed species stabilized b-CO 2 structure significantly.
However, there are still on-going discussions as to whether the lattice oxygen is stable under highly cathodic conditions. Mandal et al. synthesized Cu 2 O catalysts and investigated their reduction behavior using in situ Raman spectroscopy [38]. They found that Cu 2 O signals completely disappeared during initial 3 min at CO 2 RR conditions, and CO 2 reduction products were detected after complete reduction by selected-ion flow tube mass spectrometry (SIFT-MS). The experimental result was also supported by DFT calculation that Cu 2 O reduction is energetically preferred over CO 2 RR. Moreover, Lum et al. [39] also identified from 18 O labelling experiments where the oxygen species of OD-Cu catalysts were observed to survive less than 1% of the original oxygen content. In contrast, it has been proposed that the O sub can be stably sustained in (100) surfaces of Cu nanocubes and amorphous surface layers due to their disordered nature [40,41], although O sub is not stable in ordered Cu slab models. In fact, it has been reported that loss of Cu cube morphology by roughening during CO 2 RR caused the absence of Cu + signals and the increase in HER over CO 2 RR [42]. The aggregated amorphous nanoclusters of Cu 2 (OH) 3 Cl were also found to have the mixed metal valency of Cu + and Cu 0 at highly cathodic potentials [33].

Electrolyte engineering
The cations and anions contained in an aqueous electrolyte also play an important role in controlling the product selectivity. Resasco et al. [23] measured the partial current densities of CO 2 RR on Cu(111) using different alkali metal cations such as Li + , Na + , K + , Rb + and Cs + as shown in Fig. 3a. Interestingly, as the cation size became larger, the partial current densities for HCOO − , C 2 H 4 and C 2 H 5 OH production increased while the production rates of H 2 , CO and CH 4 were unchanged. Larger cation size facilitated the CO 2 activation over HER as well as the C 2 production over C 1 . The effect of cation on product distribution was explained by electrostatic interactions between the intermediates and the solvated cations at the outer Helmholtz plane. Based on DFT calculation, it was demonstrated that larger cations are more stabilized in the outer Helmholtz plane at the Cu (111) surface (Fig. 3b), and show higher coverage near the catalyst surface. The stabilized cations facilitated the adsorption of intermediates with large dipole moments (e.g., *CO 2 , *CO, *OCCO) over *H intermediate with a negligible dipole moment through an electric field (Fig. 3c).
It has been also suggested that the other role of hydrated cations is to control the local CO 2 concentration near the catalyst surface. Singh et al. [43] reported that cation hydrolysis occurs near the cathode surface during CO 2 reduction. As the cation migrates toward the negatively polarized electrode, an increasing electrostatic interaction can reduce the pK a for hydrolysis. If the local pH is higher than the pK a of the hydrated cations, the dissociation of a water in the hydration shell releases protons, and buffers the change of local pH. The buffering ability declines in the order of Cs + > Rb + > K + > Na + > Li + . As local CO 2 concentration competes with local pH due to the equilibrium reactions between CO 2 , HCO 3 − and CO 3 2− species, the pH decreases and the CO 2 concentration increases near the cathode with increasing cation size.
Furthermore, Dinh et al. [44] found that higher concentration of OH − ions can facilitate CO 2 activation on the Cu-coated gas diffusion layer (GDL) electrodes. When linear sweep voltammetry (LSV) curves under CO 2 and Ar flow were scanned in the electrolytes of various KOH concentrations, the onset of LSV curves recorded under CO 2 shifted markedly to more positive potentials with increasing KOH concentration, whereas the onset under Ar remained similar. Furthermore, at 10 M KOH concentration, the onset potentials for CO and C 2 H 4 production were similar, whereas the onset potential gap between CO and C 2 H 4 is observed at low KOH concentration because the build-up of CO coverage is required to proceed CO dimerization. Therefore, high KOH concentration made CO dimerization much easier. Based on DFT calculation, they proposed that the hydroxide adsorption increases the charge imbalance between carbon atoms in the *OCCO intermediate. As a result, the dipole interaction between the carbon atoms more stabilized the *OCCO, and reduced the activation barrier for CO dimerization step.

Nanostructuring
Sharpening metal can perturb the local electron density of metal surface and enhance the local electric field because free electrons migrate to sharp regions due to electrostatic repulsion. The field enhancement effect by tip electrode can be exploited for improving the CO 2 RR properties. Liu et al. [22] simulated the electron density distribution on various electrode surfaces as a function of tip radius (Fig. 4a). It was shown that the tip-concentrated electron density dramatically increases as the electrodes sharpen. The tip sharpening from a radius of 140 to 5 nm enhanced electrostatic field intensity from approximately 3 to 43 MV m −1 . The enhancement of electric field also affected the concentration of surrounding alkali metal cations. The adsorbed K + density at the tip elevated 20-fold due to the field enhancement (Fig. 4b). The highly concentrated cations can help to tune the free energy of adsorption of intermediates of CO 2 RR through the electrostatic interaction. According to DFT calculation study, the K + ions greatly lowered free energy of *COOH and *CO adsorption, in comparison with the case without K + ions (Fig. 4c). Based on these theoretical results, they attempted to measure the electrochemical CO 2 reduction properties of Au needles, rods and particles. Surprisingly, Au needle electrodes exhibited a dramatic enhancement in CO 2 RR selectivity, as well as the onset potential shift to lower overpotentials (Fig. 4d). From Kelvin probe atomic force microscopy and inductively coupled plasma optical emission spectrometry (ICP-OES), it was also confirmed that electric fields and adsorbed K + concentration are the highest for the needles.
Controlling the shape of catalysts has been reported to be effective even in C 2+ production. Jiang et al. [45] prepared a Cu nanocube catalyst with exposed (100) facets through Cu 2+ ion cycling method on Cu foil. When the cycle number increased from 0 to 100 cycles, the FE of C 2+ products was dramatically enhanced from 26.0 to 60.6%. Concurrently, the FE of C 1 products decreased from 59.1 to 22.7%. The enhanced C-C coupling capability was attributed to the exposure of (100) facets in addition to local pH effect resulting from high surface area because activity enhancement in the C 2+ production was observed in low overpotential regions where mass transport effects can be neglected. According to DFT calculation, the Cu (100) surface showed stronger CO binding as well as lower CO dimerization barrier compared to Cu (111) surface. As a result, higher selectivity toward C 2+ products was achieved. In addition to the facet-controlling strategy, constructing nanocavities at the catalyst surface can promote C 2 -C 1 coupling to form C 3 products. Zhuang et al. [34] simulated the spatial concentrations and flux distribution of CO, C 2 and C 3 species using a hollow spherical shell with a circular Fig. 4 Electric field enhancement and the impacts on CO 2 reduction properties on nanoneedle catalysts. a Free electron density distribution on the Au electrodes with different tip radius. The tip radius of the structure is 5, 60 and 140 nm for left, middle and right panels. b Adsorbed K + density and current density distributions on the surface of Au needles. c Free energy diagrams of the electrochemical CO production on Au (111) surface in the presence of adsorbed K + and in the absence of adsorbed K + . d CO FE on Au needles, rods and particles at different applied potentials (a-d Reproduced with permission from [22], copyright 2016 Nature publishing) opening. They found that the inner cavity confines the mass transport of produced C 2 species into exterior electrolyte, which results in higher local concentration of C 2 species inside the cavity. Ultimately, the cavity structure increased surface coverage of C 2 intermediate, and elevated the production rate of C 3 species. To experimentally demonstrate the cavity effect, Cu catalysts with the cavity structure were prepared through acidic etching of Cu 2 O nanoparticles and then electrochemical reduction. The nanocavity Cu catalyst exhibited the enhanced FE of 21.0% towards propanol production compared to closed Cu nanoparticles (NPs), which was in accordance with the simulation results.

Surface ligand effects
The ligand assembly can induce charge rearrangement, causing the difference in the intermediate stabilization and the adsorption configuration [46][47][48]. Kim et al. [46] synthesized Ag NP supported by carbon black in the presence of cysteamine as a ligand, and evaluated the electrocatalytic performance. The cysteamineanchored Ag NPs showed a higher maximum FE of 84.4% toward CO production in comparison with 70.5% FE of Ag foil. More importantly, the overpotential required to achieve the CO partial current density of 1 mA cm −2 was reduced as 300 mV by the ligand assembly (Fig. 5a). By DFT calculation study, it was found that the binding energy of *COOH intermediate monotonically increases as the cysteamine coverage increases, while the CO binding energy exhibited marginal changes from the ligand assembly (Fig. 5b). Specifically, the binding energy of *COOH was enhanced up to 0.26 eV as the cysteamine/Ag ratio increases from 0 to 2.7%. The preferential stabilization of *COOH over *CO led to the reduction of overpotential. Deeper mechanistic insight was taken from the spatial electron spin density analysis. When cysteamine molecules were attached, spatial localization of the unpaired electron was stabilized particularly at the Ag surface, whereas the unpaired electrons were delocalized over the entire Ag atoms without the cysteamine. The unpaired electrons made the Ag-COOH bond more covalent, and stabilized the *COOH as a radical form.
The ligand effect has been also demonstrated with the different type of ligand for Au NP. Cao et al. [48] successfully prepared the porphyrin ligand-capped and naked Au NPs (denoted as P1-AuNP and naked-AuNP, respectively). The P1-AuNP showed higher electrocatalytic activity toward CO production compared to the naked-AuNP. At − 0.45 V vs. RHE, the CO FE was measured to be 93% and 44% for P1-AuNP and naked-AuNP, respectively. Moreover, at − 0.45 V vs. RHE, the partial current density for CO production was approximately sixfold higher on P1-AuNP than on naked-AuNP. Regarding the effect of porphyrin ligand, DFT calculation study showed that the formation of *COOH and *CO intermediates is more favorable as 0.55 and 0.08 eV on P1-Au(111) surfaces than on Au(111) surfaces. The porphyrin ligand dramatically lowered the energetic barrier for *COOH formation step.

Single-atom catalyst
Recently, single-atom catalysts have attracted great attention in the research field of electrochemical CO 2 reduction due to their totally different activity compared to conventional bulk metals [49,50]. One of the representative single-atom catalysts is metal-incorporated N-doped Fig. 5 The effects of introducing ligand on the metal surfaces. a CO partial current density of cysteamine-anchored Ag/C, Ag foil, and cysteamine-C samples. b The binding energies of the COOH and CO intermediates examined from Ag (147-n) Cys n (n = 0, 1, 2, 4) models. Colored map represents the theoretical CO 2 reduction potential (E 0 ) as a function of COOH and CO binding free energies (ΔG B COOH and ΔG B CO ) (a-b Reproduced with permission from [46], copyright 2015 American Chemical Society) carbon (M-N-C) where a metal atom is coordinated by neighboring nitrogen atoms. Ju et al. [51] investigated the catalytic activity of M-N-C including various metals such as Mn, Co, Fe, Ni and Cu. They identified that the selectivity for CO production exhibit a strong dependence on the nature of metal center (Fig. 6a, b). The Co-N x sites selectively produced H 2 with the FE of 90% while the Ni-N x sites showed the highest selectivity of 85% for CO production. To understand the experimental trend, DFT calculation study was conducted (Fig. 6c). It was discovered that Co-N x sites have energetic downhill for H 2 production, but represent too strong *CO adsorption to desorb the *CO. In contrast, on the Ni-N x sites, the *H adsorption was energetically uphill reaction. Moreover, the formation and desorption process of the *CO intermediate had all-downhill energetics. In terms of the CO production rate, the Ni-NC catalyst outperformed the state-of-the-art Au catalysts. The CO 2 RR property of the Ni-NC catalyst is interesting, given that bulk Ni foil is specific for H 2 production.
The special catalytic performance of single-atom metal catalysts can be attributed to the unique electronic structure of the metal center affected by surrounding coordination atoms. According to a previous study [52], it was explained that the unstable Ni(I) atoms with an unpaired electron can be stabilized by pyridinic N with low electronegativity in an N-doped carbon matrix. The Ni(I) site can transfer charges spontaneously to the carbon 2p orbital in CO 2 , forming a CO 2 δ− species (Fig. 6d). Such an electronic interaction helps to reduce the energy barrier for electrochemical CO 2 reduction. However, it should be also considered that various N species present in the N-doped carbons can act as an active site, in addition to M-N x sites. Jung et al. investigated the catalytic performance of various M-N-C (M = Fe, Co, Cu) catalysts for CO 2 RR [53]. The maximum FE for CO production increased in the order of Cu-NC < Co-NC < Fe-NC. The maximum FE values were considered to be related with the atomic contents of pyridinic, pyrrolic, and graphitic N species based on XPS analysis. The catalysts containing higher pyridinic and graphitic N contents with lower pyrrolic N contents showed a higher selectivity toward CO production. Previous DFT calculation study also proposed that the N species in N-doped carbon nanotubes (NCNTs) can lower the energy barrier for *COOH formation because the N defects have lonepair electrons and can donate electrons to the intermediate with strong binding, although the promotional effect depends on the type of N species [54].

Conclusions
In summary, we have reviewed recent strategies for enhancing electrocatalytic properties of CO 2 reduction reactions, and discussed new insights in terms of intermediate binding energy based on the DFT calculation results. Theoretically, suppressing the H 2 evolution and improving the selectivity for a specific product on pure metal surfaces are limited by scaling relations between the binding energies of reaction intermediates, whereas it has been demonstrated that the binding energies can be tuned independently by implementing electronic perturbation and electrostatic interaction. Specifically, incoporating oxygens into the Cu surface stabilized the bent chemisorbed CO 2 and strenghthened the CO binding, which kinetically promoted C-C coupling reaction for the production of C 2 H 4 and C 2 H 5 OH. The utilization of larger alkali metal cations included in the electrolyte selectively enhanced the adsorption of intermediates with large dipole moments such as *CO 2 , *CO and *OCCO. In addition, it was revealed that the field enhancement by tip electrode increased the concentration of surrounding alkali metal cations, dramaticlaly lowering the adsorption free energies of *COOH and *CO. Furthermore, introducing surface ligand on Au and Ag NPs stabilized *COOH over *CO, which greatly elevated the selectivity and reduced the overpotential toward CO production. Also, the single-atom catalysts like Ni-N-C showed the high CO selectivity which cannot be achieved from the bulk Ni foil owing to strong *CO binidng. We believe that this review not only provide deep understandings of CO 2 reduction reactions, but also help to create new ideas for future electrocatalyst developments.
Authors' contributions CWL, CK and BKM wrote the manuscript. All authors read and approved the final manscript.