Subscriber access provided by UNIV OF ARIZONA
Article
Facile CO2 electro-reduction to formate via oxygen bidentate intermediate stabilized by high-index planes of Bi dendrite catalyst Jai Hyun Koh, Da Hye Won, Taedaehyeong Eom, Nak-Kyoon Kim, Kwang Deog Jung, Hyungjun Kim, Yun Jeong Hwang, and Byoung Koun Min ACS Catal., Just Accepted Manuscript • Publication Date (Web): 20 Jun 2017 Downloaded from http://pubs.acs.org on June 20, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
ACS Catalysis is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
Facile CO2 electro-reduction to formate via oxygen bidentate intermediate stabilized by high-index planes of Bi dendrite catalyst Jai Hyun Koh,a‡ Da Hye Won, a‡ Taedaehyeong Eom,b‡ Nak-Kyoon Kim,c Kwang Deog Jung, a Hyungjun Kim,b* Yun Jeong Hwang, a,d,* and Byoung Koun Mina,e,* a
Clean Energy Research Center, Korea Institute of Science and Technology, 5 Hwarang-ro 14-gil,
Seongbuk-gu, Seoul 02792, Republic of Korea b
Graduate School of Energy, Environment, Water, and Sustainability (EEWS), Korea Advanced
Institute of Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea c
Advanced Analysis Center, Korea Institute of Science and Technology, 5 Hwarang-ro 14-gil,
Seongbuk-gu, Seoul 02792, Republic of Korea d
Korea University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113,
Republic of Korea e
Green School, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
‡These authors contributed equally. *Corresponding Authors
1 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 24
ABSTRACT. Electrochemical CO2 conversion to chemical-products is a promising strategy for sustainable industrial development. However, the success of this approach requires an in-depth understanding of catalysis because it involves highly complex multi-step reactions. Herein, we suggest a rational design of a hierarchical Bi dendrite catalyst for an efficient conversion of CO2-to-formate. A high selectivity (~89% at –0.74 VRHE) and, more importantly, a stable performance during long-term operation (~12 h) were achieved with the Bi dendrite. Density functional theory (DFT) is used to investigate three possible reaction pathways in terms of surface intermediate, and the one via *OCOH surface intermediate is calculated to be the most energetically feasible. DFT calculations further elucidate the plane-dependent catalytic activity, and conclude that the high-index planes developed on the Bi dendrite are responsible for the sustainable performance of Bi dendrite. We expect that our experimental and theoretical study will provide a fundamental guideline to the CO2-to-formate conversion pathway as well as design principles for enhancing the catalytic performance.
KEYWORDS Electrocatalyst · CO2 Reduction · Bismuth · Dendrite · Formate
2 ACS Paragon Plus Environment
Page 3 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
INTRODUCTION Mankind has been heavily dependent on fossil fuels to meet his energy needs, and this reliance has led to various environmental problems such as global warming and the energy crisis.1 Therefore, the possibility of an artificial photosynthesis system that uses solar energy to synthesize chemical products from CO2 and water has attracted much attention as a potential breakthrough for sustainable development.2 One of the biggest challenges facing such a system is CO2 conversion because it has sluggish kinetics and is composed of highly complex multi-step reactions. To overcome this hurdle, various metal candidates (e.g., Ag3-11, Au11-20, Cu21-25, Sn26-29) have been experimentally screened for their catalytic performance on CO2 reduction and the mechanisms of the more promising ones have been extensively investigated. However, a catalyst which satisfies the requirements for practical applications, including low material cost, earthabundance, environmental friendliness and safety, is still being sorted. Bi has received attention as a catalyst of CO2 reduction due to advantages such as its low and stable cost and low-toxicity. In previous studies, almost all Bi catalysts produced CO from CO2 with high Faradaic efficiencies of 74–96% in an aprotic electrolyte (e.g., acetonitrile) containing ionic liquids.30-33 Interestingly, a few studies have recently reported that Bi-based electrodes could selectively reduce CO2 to formate in an aqueous electrolyte, which is more practical catalytic environment in terms of cost-effectiveness than an aprotic electrolyte with ionic liquid.34-36 Also, formic acid (HCOOH) and formate salts (HCOO–) are relatively highvalue commodity chemicals that are widely used in cleaning, feed preservation, and leather processing.37 A similar situation was also observed with Sn electrodes which are well-known for producing formate from CO2 in aqueous systems but in aprotic electrolytes catalyze the selective production of CO.32 These phenomena imply that the catalytic environment can tune the product
3 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 24
selectivity of the reaction by influencing the mechanism. Therefore, to design a catalyst rationally, it is important not only to modify the intrinsic properties of a catalyst experimentally, but also to understand and consider the specific mechanisms of catalysis under different operating conditions. Empirical and theoretical studies can be usefully integrated with the assistance of the density functional theory (DFT) method which can either suggest the origin of a catalytic activity based on experimental results or predict an appropriate structure for a catalyst. For instance, Zhu et al. experimentally found the optimum size of monodispersed Au nanoparticles, and their DFT study suggested that the origin of the superior CO2 electro-reduction activity was the numerous edge sites, which stabilize the reaction intermediates.13 On the contrary to this, an Au nanowire structure with a high edge-to-corner ratio was suggested based on the DFT calculations and the following experiments verified that those predicted structure did indeed exhibit the outstanding performance for electro-reduction of CO2.14 Collaborative studies integrating experimental and theoretical approaches have been conducted for the CO2 reduction reaction, and consequently, a rational design of a CO2-reducing catalyst has been proposed: the coordinately unsaturated site can effectively stabilize the reaction intermediate by lowering the energy barrier for its binding to the site.8,18,22,38-40 Based on the proposed strategy, we developed a nanostructured Bi electrode with a high density of low-coordinated sites for CO2 reduction in an aqueous solution. To enhance the product selectivity in a controlled manner, dendritic nanostructures were carefully incorporated to the surface of the Bi electrode. A dendritic morphology is promising for electrocatalysts due to their large surface area, and high degree of connectivity accompanied by a large number of active local sites.8,29 For example, these structural advantages of dendrites have been
4 ACS Paragon Plus Environment
Page 5 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
demonstrated by Sn29 and Zn8 with remarkable catalytic activities. The hierarchical Bi dendrite demonstrated a promising potential as a catalyst for CO2 reduction, showing the highest formate production rate among the Bi catalysts reported in an aqueous solution. The accompanying DFT study explored three possible reaction pathways and calculated that the most energetically feasible path was via *OCOH surface intermediate formation on the prepared Bi electrode. Furthermore, the plane dependence of the Bi catalyst proposed by the DFT calculation suggests fundamental guidelines and design principles for the development of a high performing CO2 electro-reduction catalyst.
EXPERIMENTAL SECTION Preparation of Bi dendrite electrode. Bi dendrite electrodes were deposited by a modified literature procedure.41 An electrodeposition of Bi3+ precursor was conducted in one-compartment cell using a potentiostat (PGSTAT128N, Metrohm Autolab). A three electrodes setup was employed with an Ag/AgCl (3 M NaCl, CH Instruments) and a platinum gauze (100 mesh, 99.9% trace metals basis, Sigma-Aldrich) as a reference electrode and a counter electrode, respectively. Bi thin film of 100 nm prepared by e-beam evaporation on a Cu substrate (99.98%, trace metals basis, Sigma-Aldrich) was used as a working electrode in order to minimize an amount of Bi used for the electrode. An electrolyte was 20 mM of Bi(NO3)3·5H2O (98%, Sigma-Aldrich) in ethylene glycol (anhydrous, 99.8%, Sigma-Aldrich) solution. Bi dendrite was deposited on the working electrode by repeated cycles of an electrodeposition and a resting step. A single cycle of deposition consists of a constant-potential electrodeposition at ‒1.8 V vs. Ag/AgCl for 60 s and a subsequent resting step for 2 s (Figure S1). This cycle was repeated for 10, 20, and 30 times to prepare three different Bi dendrite electrodes for electrochemical CO2 reduction. After the
5 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 24
deposition, the Bi dendrite electrodes were gently dipped in an ethanol to wash out the residual electrolytes on the surface and naturally dried at room temperature. Physical Characterization. The physical characteristics of the prepared Bi electrodes (i.e. Bi foil, e-beam Bi and Bi dendrite) were analyzed with a scanning electron microscopy (SEM), an X-ray diffraction (XRD), a high-resolution transmission electron microscopy (HR-TEM). The surface morphology of prepared Bi electrodes were confirmed by a SEM (Inspect F, FEI) at an acceleration voltage of 15.0 kV and an average working distance of 10.0 mm. The crystal structure of the prepared Bi electrodes was investigated by XRD (XRD-6000, Shimadzu) using a radiation from a Cu Kα source at 40 kV and 30 mA. Figure S2 shows the XRD patterns of Bi dendrite, e-beam Bi and polycrystalline Cu, which was used as the substrate. The morphology and crystal structure were also observed by HR-TEM (Titan 80-300TM, FEI) images and selected area diffraction (SAED) patterns. Electrochemical CO2 reduction. All electrochemical measurements were conducted in a gastight two-compartment cell of poly(ether ether ketone) (PEEK), and a proton exchange membrane (Nafion® 117) was placed in the middle of the electrochemical cell to separate catholyte and anolyte. Electrochemical CO2 reduction experiments were carried out in a CO2saturated 0.5 M KHCO3 (≥99.95%, Sigma-Aldrich) with pH of 7.3 by using a potentiostat (PGSTAT128N, Metrohm Autolab). A platinum gauze (100 mesh, 99.9% trace metals basis, Sigma-Aldrich) and the prepared Bi dendrite electrodes were used as a counter and working electrodes, respectively. A typical constant-potential electrolysis (chronoamperometry) of CO2 was conducted for an hour and steady-state current density was plotted against applied potential. All potentials were measured with a Ag/AgCl reference electrode (3 M NaCl, CH Instruments) and converted to reversible hydrogen electrode (RHE) scale using E ( . RHE) = E ( . Ag/
6 ACS Paragon Plus Environment
Page 7 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
AgCl) + 0.21 V + 0.059 V × pH . Prior to each electrochemical CO2 reduction experiment, solution resistance (Rs) was measured by electrochemical impedance spectroscopy taken at various potentials (Figure S3). The potential converted to RHE scale was further corrected for ohmic losses (iRs). Quantification of CO2 reduction products. The gas and liquid products from CO2 reduction were analyzed by using a gas chromatography (GC) and a nuclear magnetic resonance (NMR) spectrometer, respectively. Gas products as H2 and CO in the headspace of the electrochemical cell were quantitatively analyzed by GC (Younglin 6500 GC system) equipped with a 80/100 Carbosieve SII column (Supleco Analytical), a thermal conductivity detector (TCD), and a flame ionization detector (FID). The quantification of gas products were performed at an interval of ten minutes during the electrolysis for an hour. Ultra-high purity Ar (99.9999%) was used as a carrier gas. GC was directly connected to the cathodic side of the electrochemical cell and the gas products were injected through a six-port valve. Each compartment of the electrochemical cell contained 38 mL of electrolyte and 34 mL of headspace, and CO2 gas was purged through the electrolyte at 20 sccm. The Faradaic efficiency of gas products were evaluated from the area of each peak in the gas chromatogram by using the following equation: Faradaic efficiency = × ×
2
where Vi is the volume concentration of H2 or CO based on a calibration of the GC, F = 96,485 C mol , = 1.013 bar , R = 8.314 J mol K , T = 298.15 K , and itotal is total steady-state current measured during a constant-potential electrolysis. The quantification of liquid products were performed after the electrolysis for an hour on an Agilent 600 MHz NMR spectroscopy. 450 μL of electrolyte after CO2 reduction electrolysis was mixed with an internal standard, 50 μL of 1% 3-(trimethylsilyl)-1-propanesulfonic acid
7 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 24
(DSS, Sigma-Aldrich) in D2O, for NMR analysis. The H2O peak was suppressed by presaturation method. The ratio of the area of each liquid product peak to that of DSS was acquired to quantify the amount of the liquid product of CO2 reduction. The amount of charge passed to produce each liquid product was divided by the amount of total charge passed during the electrolysis to evaluate each Faradaic efficiency. The sum of Faradaic efficiencies of all gaseous and liquid products was normalized to 100%. Computational details & reduction potential calculation for CO2 reduction. We performed spin-polarized DFT calculations using Vienna Ab initio simulation package (VASP)42 with Perdew-Burke-Ernzerhof (PBE) exchange correlational functional43 and plane-wave basis sets. We first optimized the bulk rhombohedral Bi structure, yielding the lattice parameters of a = 4.809 Å and a = 57.03°, which are comparable to the experimental lattice parameters of a = 4.746 Å and a = 57.24° (Table S1). We then built slap models by cleaving the optimized bulk structure of Bi with the surface orientation of (003), (012), (110), and (104) (Figure S4), where the bottom layer is fixed and ~15 Å of vacuum layer is additionally included to avoid the artificial inter-slap interaction. For the reciprocal space sampling, we used 9×9×9 mesh for the optimization of primitive Bi bulk structure, 2×2×1 mesh for Bi (003) and (012) slab models, 2×3×1 mesh for (110) slab model, and 3×2×1 mesh for (104) slab model. We set the kinetic energy cutoff as 600 eV, and used the projector augmented wave (PAW) method44 to describe the core electrons. We computed the binding free energies (ΔGB) of possible intermediate states of H, COOH, OCOH, and HCOOH by diagonalizing the partial Hessian matrix to include the vibrational entropy and enthalpy to the DFT self-consistent-field (SCF) energy. We note that our calculation yields the thermodynamic potential of CO2 reduction to HCOOH as –0.22 V vs.
8 ACS Paragon Plus Environment
Page 9 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
standard hydrogen electrode (SHE), which is comparable to the experimental value of –0.20 V vs. SHE.45
9 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 24
RESULTS AND DISCUSSION
Scheme 1. Schematic illustration for the preparation of Bi dendrite electrode. Bi dendrite was prepared by e-beam evaporation of Bi source on a Cu foil, followed by electrodeposition of Bi3+ precursor in ethylene glycol solution. The Bi dendrite electrode was prepared by multi-step electrodeposition of Bi3+ precursor in ethylene glycol solution onto a 100 nm-thick Bi thin film substrate using different step cycles (Scheme 1). A Bi thin film substrate was prepared by e-beam deposition on a Cu substrate and had a rough surface comprised of hundreds-of-nanometer-sized particles (Figure 1b), which could facilitate deposition by improving the adhesion of deposits. According to the SEM images, the Bi dendrites deposited over 30 cycles showed a well-defined hierarchical structure in which smaller branches grew on bigger branches (Figure 1c and S5). The CO2 reduction activity of Bi dendrite electrodes increased with the number of deposition cycles used to create them and leveled at 30 cycles (Figure S6). Therefore, 30 cycles was determined to be the optimum number of cycles for Bi dendrite deposition and the subsequent CO2 reduction experiments were conducted on the corresponding electrode. The crystalline structure of the prepared Bi dendrite was examined by XRD. As shown in Figure 1d, the XRD patterns of the Bi dendrite wellmatched rhombohedral polycrystalline Bi (JCPDS No. 85-1330). Compared to a pristine Bi foil, relatively higher index planes such as (012), (107), (116), (214) and (009) were developed on the Bi dendrite surface, suggesting a well-formed dendrite structure (Figure S7).
10 ACS Paragon Plus Environment
Page 11 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
Figure 1. Characterization of Bi dendrite electrode. Top-down SEM images of (a) pristine Bi foil, (b) e-beam deposited Bi, and (c) Bi dendrite. (d) XRD patterns of as-prepared Bi dendrite and the pristine Bi foil. To investigate the catalytic activity and the selectivity on the prepared Bi dendrite in the aqueous solution, the electrochemical CO2 reduction was conducted at various constant potentials from ‒0.56 to ‒1.19 V [vs. RHE; all potentials are given with respect to RHE] for an hour in a CO2-saturated 0.5 M KHCO3 aqueous electrolyte (Figure S8). A commercial pristine Bi foil was employed as a control group in terms of CO2 reduction activity. The liquid and gas products of the CO2 reduction were analyzed by a NMR spectrometer and a GC, respectively, and the resulting data were used to calculate Faradaic efficiencies and partial current densities (Figure S9). The primary product of CO2 reduction on the prepared Bi electrodes was formate (HCOO–), and a trace of CO was also detected. H2 was produced as a byproduct by the hydrogen evolution reaction (HER), which is a competitive reaction with CO2 reduction. During CO2 electrolysis, the Bi dendrite showed a much higher current density with an improved Faradaic efficiency for formate compared to the pristine Bi foil (Figure 2a‒b). As shown in Figure 2b, the
11 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 24
Bi dendrite exhibited a stable current density of ‒2.7 mA cm‒2 with a Faradaic efficiency of ~89% for formate at ‒0.74 V whereas the pristine Bi foil exhibited a current density of only ‒0.4 mA cm‒2 with a Faradaic efficiency of ~56% at the same potential. Compared to the pristine Bi foil, the Bi dendrite electrode had a positively shifted onset potential, i.e., the potential at the start of the CO2 reduction reaction (Figure 2c‒d). Similarly, the overpotential of the dendritic electrode at the maximum formate production Faradaic efficiency was reduced by 180 mV. While the pristine Bi foil required an overpotential of 890 mV to reach its maximum formate Faradaic efficiency of ~79%, the Bi dendrite showed an anodically shifted overpotential of 710 mV to achieve its maximum Faradaic efficiency of ~89%. Furthermore, due to outstanding current densities and Faradaic efficiencies, the production rates of formate on the Bi dendrite surpassed not only those on the Bi foil, but also those on state-of-the-art Bi-based catalysts and on other formate-producing catalysts, such as Sn and Pb (Figure 2e and Table S3). Although the HER rate on the Bi dendrite exceeded the CO2 reduction rate at a high overpotential region over ‒1.04 V (Figure S10), this could be attributed to insufficient mass transport of CO2 into the layers of hierarchical dendritic structures.8,25
12 ACS Paragon Plus Environment
Page 13 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
Figure 2. Catalytic activities of pristine Bi foil and Bi dendrite. (a) Current densities during electrochemical CO2 reduction at various constant potentials in a CO2-saturated 0.5 M KHCO3. (b) Comparison of performance in terms of total current density (line; left axis) and Faradaic efficiency of formate (dot; right axis) at ‒0.74 V during 1 h. Faradaic efficiencies of (c) Bi foil and (d) Bi dendrite electrodes calculated after operation for 1 h of CO2 reduction at various applied potentials. (e) Formate production rate at various applied potentials on Bi foil and Bi dendrite electrodes. Since a durability of a catalyst is crucial for practical applications, the performance profile of the Bi dendrite during a long-term operation was examined by collecting the current density and Faradaic efficiency at ‒0.73 V (Figure 3). During constant-potential electrolysis, the Bi dendrite showed an outstanding stability over 12 h without any significant performance degradation. Faradaic efficiency, as well as current density, of formate production were also
13 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 24
maintained above 90% of their initial performances. In addition, negligible change in morphology (Figure 3b–c) and XRD patterns (Figure S11) implies that the Bi dendrite morphology is sufficiently stable during the long-term electrolysis. This result is particularly promising because Bi-based electrodes including their durability issues have not yet been extensively investigated in previous studies. Thus, our Bi dendrite electrode demonstrated significant potential as an electrocatalyst for CO2 reduction with high catalytic activity, selectivity, and stability.
Figure 3. (a) Current density (line; left axis) and Faradaic efficiency of formate (dot; right axis) on the Bi dendrite electrode during long-term operation of 12 h at ‒0.73 V. Morphology of the Bi dendrite electrode (b) before and (c) after CO2 reduction reaction for 12 h. It should be noted that the Bi dendrite exhibited 11 times higher CO2 reduction activity than that of the pristine Bi foil in terms of formate production rate at ‒0.74 V (Figure 2e and S10e). This 11 times higher CO2 reduction activity, despite the fact that Bi dendrite having only
14 ACS Paragon Plus Environment
Page 15 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
a 3.1-fold higher surface area than Bi foil according to the electrochemical surface area measurements (Figure S12), implies that the intrinsic catalytic activity of Bi dendrite was remarkably enhanced, perhaps mechanistically. To understand the mechanism of CO2 reduction on the electrode surface, Tafel plots were acquired at regions of low current density, where the reaction is mainly limited by electrokinetics (Figure S13). The dendritic Bi electrode had a Tafel slope of 92 mV dec‒1, smaller than that of pristine Bi foil of 150 mV dec‒1, and both values differed from 118 mV dec‒1. It is generally accepted that a reaction having a Tafel slope near 118 mV dec‒1 involves an initial one electron transfer step from catalyst to CO2 as the ratedetermining step (RDS), but neither electrode showed this value.12 Therefore, we can conclude that the RDSs on the Bi surfaces differ from the typical one electron transfer step, and that dendritic surfaces can lead to different mechanistic pathways from the pristine Bi foil as high index surfaces are exposed. In general, it is a well-known strategy to develop nanostructures on a catalyst surface to achieve enhanced catalytic performance.8,11,17,25,46 To gain a deeper understanding on the origin of the improved catalytic performance other than the effect of nanostructuring, DFT calculations were conducted considering possible reaction paths under the operating conditions (Figure 4). Three different possible reaction paths were suggested, as follows (an asterisk * indicates the catalytic surface): Path 1 via *COOH * + CO2 (g) + [H+ + e‒] ® *COOH (aq)
(R1-1; slow)
*COOH (aq) + [H+ + e‒] ® *HCOOH (aq)
(R1-2; fast)
*HCOOH (aq) ® * + HCOOH (aq)
(R1-3)
15 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 24
Path 2 via *OCOH * + CO2 (g) + [H+ + e‒] ® *OCOH (aq)
(R2-1; slow)
*OCOH (aq) + [H+ + e‒] ® *HCOOH (aq)
(R2-2; fast)
*HCOOH (aq) ® * + HCOOH (aq)
(R2-3)
Path 3 via *H * + [H+ + e‒] ® *H (aq)
(R3-1; slow)
*H (aq) + CO2 (g) + [H+ + e‒] ® *HCOOH (aq)
(R3-2; fast)
*HCOOH (aq) ® * + HCOOH (aq)
(R3-3)
They are schematically shown in Figure 4a. The difference between path 1 and path 2 is merely whether the COOH intermediate is adsorbed via a Bi-C interaction (i.e. *COOH) or a Bi-O interaction (i.e. *OCOH), respectively. In path 3 dissolved CO2 is added to surface hydrogen after a reversible hydrogen adsorption/desorption on the electrode surface, a mechanistic path which, as reported in previous research, could be possible in the aqueous electrolyte.47
16 ACS Paragon Plus Environment
Page 17 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
Figure 4. DFT calculation results on the CO2-to-formate reduction reaction pathways for various possible Bi planes of (003), (012), (110), and (104). (a) Schematic diagram of three possible reaction pathways; path 1 via the formation of *COOH intermediate, path 2 via the formation of *OCOH intermediate, and path 3 via the formation of *H. Reaction free energy diagrams for the reaction (b) path 1, (c) path 2, and (d) path 3, when zero overpotential is applied (bias potential U = ‒0.21 VSHE). Noting that (012), (110), and (104) planes are preferentially developed for Bi dendrite compared with the close-packed surface of (003) that is dominant in Bi foil. Using DFT, we calculated the reaction free energies of the above three reaction paths by diagonalizing the partial Hessian to evaluate zero point energies and vibrational free energies.
17 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 24
The computational hydrogen electrode method was used to determine the chemical potential of the proton and electron pair (H+ + e‒)48. We considered different planes of Bi surfaces, namely the close-packed surface of (003) and high-index surfaces of (012), (110), and (104) planes. For all three possible reaction paths on the close-packed surface of (003) plane, DFT results show that the first electron transfer step (R1-1, R2-1, or R3-1) is the most energetically difficult step. We also find that path 2, which is via the formation of *OCOH surface intermediate, provides the energetically most favorable path regardless of the choice of surface plane. It is of noted that the water participation in the reaction (which is not considered in our calculations) might kinetically prefer the O-H bond formation (i.e., *COOH path) over the C-H bond formation (i.e., *OCOH path)49-52. However, due to the strong thermodynamic preference toward *OCOH path, we still expect the dominant contribution of *OCOH path. Given that the dominant reaction pathway is path 2, we find that the *OCOH surface intermediate can be substantially more stabilized on all the explored high-index planes, (012), (110), and (104) than on the (003) plane. For the other possible pathways, namely path 1 and path 3, which could have a non-negligible contribution particularly under a high overpotential regime, we find that the (110) plane can serve as a better catalytic active site than the (003) plane. As our discussion being solely based on thermodynamic barriers, we estimated the potential dependent activation barriers (Eact)52-53 for the 1st and 2nd elementary steps of path 2 on the (012) surface as a representative example (Figures S14 and S15, and Table S2), leading to no barrier for the 1st step and only Eact ~ 0.2 eV for the 2nd step under the bias potential of –0.75 V (experimental condition). This shows that the barriers become surmountable at the applied potential, and thus confirms the validity of our discussion based on thermodynamic barriers.
18 ACS Paragon Plus Environment
Page 19 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
Indeed, as shown in the XRD pattern of Figure 1d, the (012), (110) and (104) planes were preferentially developed on the Bi dendrite, while the (003) plane was preferentially developed on the pristine Bi foil. In conjunction with DFT results, this suggests that the morphological deviation of the Bi dendrite (exposing high-index planes) contributes to the superior CO2 reduction to formate. Although we were not able to directly perform further DFT calculations due to limited computational resources, we further conceive that additional high index surfaces only present on the Bi dendrite electrode, such as (107), (116), (214) and (009) amongst others, could also positively affect CO2 reduction by providing low-coordinated sites with appropriate binding energy.
CONCLUSIONS In summary, a hierarchical Bi dendrite electrode was developed for the selective conversion of CO2 to formate in an aqueous electrolyte. The Bi dendrite exhibits high catalytic activity and selectivity with a Faradaic efficiency of ~89% at an overpotential of 710 mV as well as stable performance during 12 h of operation. Using DFT calculation, we investigated three possible CO2 reduction pathways, among which the path via the formation of *OCOH surface intermediate is the most favorable in terms of reaction free energetics. DFT results further elucidated that the high-index planes can efficiently stabilize the *OCOH intermediate to enhance CO2 reduction activity, which successfully explains the origin of the superior performance of the Bi dendrite compared to that of the Bi foil. Based on our structure-DFT corroborated study, we provide a greater understanding of the governing factors of CO2-toformate reduction catalysis, which can be utilized in the further development of various advanced metal electrodes for CO2 reduction with high activity, selectivity, and stability.
19 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 24
AUTHOR INFORMATION Corresponding Author *Email:
[email protected], Phone: +82-42-350-1725 *Email:
[email protected], Phone: +82-2-958-5227, Fax: +82-2-958-5809 *Email:
[email protected], Phone: +82-2-958-5853, Fax: +82-2-958-5809 Author Contributions ‡These authors contributed equally. Notes The authors declare no competing financial interests. ASSOCIATED CONTENT Supporting Information. Additional information about electrode preparation and DFT calculations, physical characterizations of the electrode, detail electrochemical results of CO2 reduction, and a supplementary table for performance comparison. ACKNOWLEDGMENT This work was supported by the program of the Korea Institute of Science and Technology (KIST) and partly by the KU-KIST program by the Ministry of Science, ICT and Future Planning.
20 ACS Paragon Plus Environment
Page 21 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
REFERENCES (1)
Olah, G. A.; Prakash, G. K. S.; Goeppert, A. J. Am. Chem. Soc. 2011, 133, 12881-12898.
(2)
Centi, G.; Quadrelli, E. A.; Perathoner, S. Energy Environ. Sci. 2013, 6, 1711-1731.
(3)
Rosen, B. A.; Salehi-Khojin, A.; Thorson, M. R.; Zhu, W.; Whipple, D. T.; Kenis, P. J.
A.; Masel, R. I. Science 2011, 334, 643-644. (4)
Salehi-Khojin, A.; Jhong, H.-R. M.; Rosen, B. A.; Zhu, W.; Ma, S.; Kenis, P. J. A.;
Masel, R. I. J. Phys. Chem. C 2012, 117, 1627-1632. (5)
Lim, H.-K.; Shin, H.; Goddard, W. A.; Hwang, Y. J.; Min, B. K.; Kim, H. J. Am. Chem.
Soc. 2014, 136, 11355-11361. (6)
Hsieh, Y.-C.; Senanayake, S. D.; Zhang, Y.; Xu, W.; Polyansky, D. E. ACS Catal. 2015,
5, 5349-5356. (7)
Kim, C.; Jeon, H. S.; Eom, T.; Jee, M. S.; Kim, H.; Friend, C. M.; Min, B. K.; Hwang, Y.
J. J. Am. Chem. Soc. 2015, 137, 13844-13850. (8)
Rosen, J.; Hutchings, G. S.; Lu, Q.; Rivera, S.; Zhou, Y.; Vlachos, D. G.; Jiao, F. ACS
Catal. 2015, 5, 4293-4299. (9)
Jee, M. S.; Jeon, H. S.; Kim, C.; Lee, H.; Koh, J. H.; Cho, J.; Min, B. K.; Hwang, Y. J.
Appl. Catal., B 2016, 180, 372-378. (10) Kim, C.; Eom, T.; Jee, M. S.; Jung, H.; Kim, H.; Min, B. K.; Hwang, Y. J. ACS Catal. 2017, 7, 779-785. (11) Kuhl, K. P.; Hatsukade, T.; Cave, E. R.; Abram, D. N.; Kibsgaard, J.; Jaramillo, T. F. J. Am. Chem. Soc. 2014, 136, 14107-14113. (12) Chen, Y.; Li, C. W.; Kanan, M. W. J. Am. Chem. Soc. 2012, 134, 19969-19972. (13) Zhu, W.; Michalsky, R.; Metin, Ö.; Lv, H.; Guo, S.; Wright, C. J.; Sun, X.; Peterson, A. A.; Sun, S. J. Am. Chem. Soc. 2013, 135, 16833-16836. (14) Zhu, W.; Zhang, Y.-J.; Zhang, H.; Lv, H.; Li, Q.; Michalsky, R.; Peterson, A. A.; Sun, S. J. Am. Chem. Soc. 2014, 136, 16132-16135. (15) Mistry, H.; Reske, R.; Zeng, Z.; Zhao, Z.-J.; Greeley, J.; Strasser, P.; Cuenya, B. R. J. Am. Chem. Soc. 2014, 136, 16473-16476. (16) Koh, J. H.; Jeon, H. S.; Jee, M. S.; Nursanto, E. B.; Lee, H.; Hwang, Y. J.; Min, B. K. J. Phys. Chem. C 2014, 119, 883-889.
21 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 24
(17) Feng, X.; Jiang, K.; Fan, S.; Kanan, M. W. J. Am. Chem. Soc. 2015, 137, 4606-4609. (18) Liu, M.; Pang, Y.; Zhang, B.; De Luna, P.; Voznyy, O.; Xu, J.; Zheng, X.; Dinh, C. T.; Fan, F.; Cao, C.; de Arquer, F. P. G.; Safaei, T. S.; Mepham, A.; Klinkova, A.; Kumacheva, E.; Filleter, T.; Sinton, D.; Kelley, S. O.; Sargent, E. H. Nature 2016, 537, 382-386. (19) Kim, K.-S.; Kim, W. J.; Lim, H.-K.; Lee, E. K.; Kim, H. ACS Catal. 2016, 6, 4443-4448. (20) Kim, H.; Jeon, H. S.; Jee, M. S.; Nursanto, E. B.; Singh, J. P.; Chae, K.; Hwang, Y. J.; Min, B. K. ChemSusChem 2016, 9, 2097-2102. (21) Tang, W.; Peterson, A. A.; Varela, A. S.; Jovanov, Z. P.; Bech, L.; Durand, W. J.; Dahl, S.; Norskov, J. K.; Chorkendorff, I. Phys. Chem. Chem. Phys. 2012, 14, 76-81. (22) Reske, R.; Mistry, H.; Behafarid, F.; Roldan Cuenya, B.; Strasser, P. J. Am. Chem. Soc. 2014, 136, 6978-6986. (23) Roberts, F. S.; Kuhl, K. P.; Nilsson, A. Angew. Chem., Int. Ed. 2015, 54, 5179-5182. (24) Varela, A. S.; Ju, W.; Reier, T.; Strasser, P. ACS Catal. 2016, 2136-2144. (25) Sen, S.; Liu, D.; Palmore, G. T. R. ACS Catal. 2014, 4, 3091-3095. (26) Chen, Y.; Kanan, M. W. J. Am. Chem. Soc. 2012, 134, 1986-1989. (27) Zhang, S.; Kang, P.; Meyer, T. J. J. Am. Chem. Soc. 2014, 136, 1734-1737. (28) Baruch, M. F.; Pander, J. E.; White, J. L.; Bocarsly, A. B. ACS Catal. 2015, 5, 3148-3156. (29) Won, D. H.; Choi, C. H.; Chung, J.; Chung, M. W.; Kim, E.-H.; Woo, S. I. ChemSusChem 2015, 8, 3092-3098. (30) DiMeglio, J. L.; Rosenthal, J. J. Am. Chem. Soc. 2013, 135, 8798-8801. (31) Medina-Ramos, J.; DiMeglio, J. L.; Rosenthal, J. J. Am. Chem. Soc. 2014, 136, 83618367. (32) Medina-Ramos, J.; Pupillo, R. C.; Keane, T. P.; DiMeglio, J. L.; Rosenthal, J. J. Am. Chem. Soc. 2015, 137, 5021-5027. (33) Zhang, Z.; Chi, M.; Veith, G. M.; Zhang, P.; Lutterman, D. A.; Rosenthal, J.; Overbury, S. H.; Dai, S.; Zhu, H. ACS Catal. 2016, 6, 6255-6264. (34) Zhang, H.; Ma, Y.; Quan, F.; Huang, J.; Jia, F.; Zhang, L. Electrochem. Commun. 2014, 46, 63-66. (35) Zhong, H.; Qiu, Y.; Zhang, T.; Li, X.; Zhang, H.; Chen, X. J. Mater. Chem. A 2016, 4, 13746-13753.
22 ACS Paragon Plus Environment
Page 23 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Catalysis
(36) Lv, W.; Zhou, J.; Bei, J.; Zhang, R.; Wang, L.; Xu, Q.; Wang, W. Appl. Surf. Sci. 2017, 393, 191-196. (37) Aresta, M.; Dibenedetto, A.; Angelini, A. Chem. Rev. 2014, 114, 1709-1742. (38) Gao, D.; Zhou, H.; Wang, J.; Miao, S.; Yang, F.; Wang, G.; Wang, J.; Bao, X. J. Am. Chem. Soc. 2015, 137, 4288-4291. (39) Rasul, S.; Anjum, D. H.; Jedidi, A.; Minenkov, Y.; Cavallo, L.; Takanabe, K. Angew. Chem., Int. Ed. 2015, 54, 2146-2150. (40) Won, D. H.; Shin, H.; Koh, J.; Chung, J.; Lee, H. S.; Kim, H.; Woo, S. I. Angew. Chem., Int. Ed. 2016, 55, 9297-9300. (41) Kang, D.; Park, Y.; Hill, J. C.; Choi, K.-S. J. Phys. Chem. Lett. 2014, 5, 2994-2999. (42) Kresse, G.; Furthmüller, J. Phys. Rev. B 1996, 54, 11169-11186. (43) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1997, 78, 1396-1396. (44) Blöchl, P. E. Phys. Rev. B 1994, 50, 17953-17979. (45) Allen J. Bard, R. P., Joseph Jordan, Standard Potentials in Aqueous Solution. CRC Press: 1985; p 195. (46) Rosen, J.; Hutchings, G. S.; Lu, Q.; Forest, R. V.; Moore, A.; Jiao, F. ACS Catal. 2015, 5, 4586-4591. (47) Min, X.; Kanan, M. W. J. Am. Chem. Soc. 2015, 137, 4701-4708. (48) Norskov, J. K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J. R.; Bligaard, T.; Jonsson, H. J. Phys. Chem. B 2004, 108, 17886-17892. (49) Michel, C.; Auneau, F.; Delbecq, F.; Sautet, P. ACS Catal. 2011, 1, 1430-1440. (50) Zhao, Y. F.; Yang, Y.; Mims, C.; Peden, C. H. F.; Li, J.; Mei, D. H. J. Catal. 2011, 281, 199-211. (51) Yang, Y.; Mims, C. A.; Mei, D. H.; Peden, C. H. F.; Campbell, C. T. J. Catal. 2013, 298, 10-17. (52) Nie, X. W.; Luo, W. J.; Janik, M. J.; Asthagiri, A. J. Catal. 2014, 312, 108-122. (53) Nie, X.; Esopi, M. R.; Janik, M. J.; Asthagiri, A. Angew. Chem., Int. Ed. 2013, 52, 24592462.
23 ACS Paragon Plus Environment
ACS Catalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 24 of 24
TOC GRAPHICS
24 ACS Paragon Plus Environment