CO2 electroreduction to ethylene via hydroxide-mediated copper catalysis at an abrupt interface
18 May 2018-Science (American Association for the Advancement of Science)-Vol. 360, Iss: 6390, pp 783-787
TL;DR: A copper electrocatalyst at an abrupt reaction interface in an alkaline electrolyte reduces CO2 to ethylene with 70% faradaic efficiency at a potential of −0.55 volts versus a reversible hydrogen electrode (RHE).
Abstract: Carbon dioxide (CO 2 ) electroreduction could provide a useful source of ethylene, but low conversion efficiency, low production rates, and low catalyst stability limit current systems. Here we report that a copper electrocatalyst at an abrupt reaction interface in an alkaline electrolyte reduces CO 2 to ethylene with 70% faradaic efficiency at a potential of −0.55 volts versus a reversible hydrogen electrode (RHE). Hydroxide ions on or near the copper surface lower the CO 2 reduction and carbon monoxide (CO)–CO coupling activation energy barriers; as a result, onset of ethylene evolution at −0.165 volts versus an RHE in 10 molar potassium hydroxide occurs almost simultaneously with CO production. Operational stability was enhanced via the introduction of a polymer-based gas diffusion layer that sandwiches the reaction interface between separate hydrophobic and conductive supports, providing constant ethylene selectivity for an initial 150 operating hours.
Citations
More filters
TL;DR: The results of this work highlight the importance of carrying out a carbon balance, in addition to traditional measurements of activity and selectivity, to adequately assess the performance of CO2 reduction devices at high current densities, and inform future efforts aimed at mitigating membrane crossover in MEA-type electrolyzers forCO2 reduction.
Abstract: Cell designs that integrate membrane-electrode assemblies (MEAs) with highly selective catalysts are a promising route to reduce ohmic losses and achieve high energy efficiency in CO2 reduction at industrially relevant current densities. In this work, porous silver filtration membranes are demonstrated as simple and efficient gas-diffusion electrodes for CO2 reduction to CO at high current densities in an MEA-type device. A partial current density for CO of up to ca. 200 mA cm-2 was achieved at a cell voltage of ca. 3.3 V, in tandem with minimal H2 production. However, the analysis of cathodic and anodic outlet streams revealed that CO2 cross-over across the anion-exchange membranes, mostly in the form of CO32- but partially as HCOO- generated over the cathode, actually exceeds the amount of CO2 converted to the target product, resulting in a poor utilization of the reactant and in the early onset of mass transfer limitations. In addition, CO2 cross-over leads to a nonstoichiometric decrease of the outlet flow rate from the cathodic compartment. This effect can lead to a substantial overestimation of catalytic performance if the inlet flow rate of CO2 is used as reference for calculating partial current densities and Faradaic efficiencies. The results of this work highlight the importance of carrying out a carbon balance, in addition to traditional measurements of activity and selectivity, to adequately assess the performance of CO2 reduction devices at high current densities, and inform future efforts aimed at mitigating membrane cross-over in MEA-type electrolyzers for CO2 reduction.
140 citations
TL;DR: This study presents a new perovskite-type copper(I) nitride (Cu3N) nanocube (NC) catalyst for selective CO2RR, and presents a good example of utilizing metal nitrides as highly efficient nanocatalysts for selectiveCO2RR to hydrocarbons that will be important for sustainable chemistry/energy applications.
Abstract: Understanding the Cu-catalyzed electrochemical CO2 reduction reaction (CO2RR) under ambient conditions is both fundamentally interesting and technologically important for selective CO2RR to hydroca...
139 citations
TL;DR: In this paper, two stable copper(I)-based coordination polymer (NNU-32 and NNU-33(S)) catalysts are synthesized and integrated into a CO2 flow cell electrolyzer, which exhibited very high selectivity for electrocatalytic CO2-to-CH4 conversion due to clearly inherent intramolecular cuprophilic interactions.
Abstract: Cu(I)-based catalysts have proven to play an important role in the formation of specific hydrocarbon products from electrochemical carbon dioxide reduction reaction (CO2RR). However, it is difficult to understand the effect of intrinsic cuprophilic interactions inside the Cu(I) catalysts on the electrocatalytic mechanism and performance. Herein, two stable copper(I)-based coordination polymer (NNU-32 and NNU-33(S)) catalysts are synthesized and integrated into a CO2 flow cell electrolyzer, which exhibited very high selectivity for electrocatalytic CO2-to-CH4 conversion due to clearly inherent intramolecular cuprophilic interactions. Substitution of hydroxyl radicals for sulfate radicals during the electrocatalytic process results in an in situ dynamic crystal structure transition from NNU-33(S) to NNU-33(H), which further strengthens the cuprophilic interactions inside the catalyst structure. Consequently, NNU-33(H) with enhanced cuprophilic interactions shows an outstanding product (CH4) selectivity of 82% at -0.9 V (vs reversible hydrogen electrode, j = 391 mA cm-2), which represents the best crystalline catalyst for electrocatalytic CO2-to-CH4 conversion to date. Moreover, the detailed DFT calculations also prove that the cuprophilic interactions can effectively facilitate the electroreduction of CO2 to CH4 by decreasing the Gibbs free energy change of potential determining step (*H2COOH → *OCH2). Significantly, this work first explored the effect of intrinsic cuprophilic interactions of Cu(I)-based catalysts on the electrocatalytic performance of CO2RR and provides an important case study for designing more stable and efficient crystalline catalysts to reduce CO2 to high-value carbon products.
139 citations
TL;DR: In this article, the authors investigated the reactivities and coverages of atop-and bridge-bound CO on a polycrystalline Cu electrode in contact with alkaline electrolytes and showed that a fraction of atopbound CO converts to bridge-bonded CO when the total CO coverage drops below the saturation coverage.
Abstract: Surface-adsorbed CO is generally considered a reactive on-pathway intermediate in the aqueous electrochemical reduction of CO2 on Cu electrodes. Though CO can bind to a variety of adsorption sites (e.g., atop or bridge), spectroscopic studies of the Cu/electrolyte contact have mostly been concerned with atop-bound CO. Using surface-selective infrared (IR) spectroscopy, we have investigated the reactivities and coverages of atop- and bridge-bound CO on a polycrystalline Cu electrode in contact with alkaline electrolytes. We show here that (1) a fraction of atop-bound CO converts to bridge-bonded CO when the total CO coverage drops below the saturation coverage and (2) unlike atop-bound CO, bridge-bonded CO is an unreactive species that is not reduced at a potential of −1.75 V vs SHE. Our results suggest that bridge-bonded CO is not an on-pathway intermediate in CO reduction. Using density functional theory (DFT) calculations, we further reveal that the activation barrier for the hydrogenation of bridge-bon...
136 citations
01 Aug 2019
TL;DR: In this article, the roughness factor of an electrode has been used to increase total rates of production, although rarely as a means to improve selectivity, which is an effective design principle to direct the selectivity of CO reduction towards multicarbon oxygenates at low overpotentials and concurrently suppressing hydrocarbon and hydrogen production.
Abstract: Using renewable electricity to convert CO/CO2 into liquid products is touted as a sustainable process to produce fuels and chemicals, yet requires further advances in electrocatalyst understanding, development and device integration. The roughness factor of an electrode has generally been used to increase total rates of production, although rarely as a means to improve selectivity. Here we demonstrate that increasing the roughness factor of Cu electrodes is an effective design principle to direct the selectivity of CO reduction towards multicarbon oxygenates at low overpotentials and concurrently suppressing hydrocarbon and hydrogen production. The nanostructured Cu electrodes are capable of achieving almost full selectivity towards multicarbon oxygenates at an electrode potential of only –0.23 V versus the reversible hydrogen electrode. The successful implementation of this catalytic system has enabled an excellent CO reduction performance and elucidated viable pathways to improve the energy efficiency towards liquid fuels in high-power conversion electrolysers. The roughness factor of an electrode has been generally used to increase total rates of production, though rarely as a means to improve selectivity. Now, Jaramillo, Hahn and co-workers direct the selectivity of CO reduction to multicarbon oxygenates at low overpotentials by increasing the roughness factor of nanostructured Cu electrodes.
136 citations
References
More filters
TL;DR: A simple derivation of a simple GGA is presented, in which all parameters (other than those in LSD) are fundamental constants, and only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked.
Abstract: Generalized gradient approximations (GGA’s) for the exchange-correlation energy improve upon the local spin density (LSD) description of atoms, molecules, and solids. We present a simple derivation of a simple GGA, in which all parameters (other than those in LSD) are fundamental constants. Only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked. Improvements over PW91 include an accurate description of the linear response of the uniform electron gas, correct behavior under uniform scaling, and a smoother potential. [S0031-9007(96)01479-2] PACS numbers: 71.15.Mb, 71.45.Gm Kohn-Sham density functional theory [1,2] is widely used for self-consistent-field electronic structure calculations of the ground-state properties of atoms, molecules, and solids. In this theory, only the exchange-correlation energy EXC › EX 1 EC as a functional of the electron spin densities n"srd and n#srd must be approximated. The most popular functionals have a form appropriate for slowly varying densities: the local spin density (LSD) approximation Z d 3 rn e unif
146,533 citations
TL;DR: An efficient scheme for calculating the Kohn-Sham ground state of metallic systems using pseudopotentials and a plane-wave basis set is presented and the application of Pulay's DIIS method to the iterative diagonalization of large matrices will be discussed.
Abstract: We present an efficient scheme for calculating the Kohn-Sham ground state of metallic systems using pseudopotentials and a plane-wave basis set. In the first part the application of Pulay's DIIS method (direct inversion in the iterative subspace) to the iterative diagonalization of large matrices will be discussed. Our approach is stable, reliable, and minimizes the number of order ${\mathit{N}}_{\mathrm{atoms}}^{3}$ operations. In the second part, we will discuss an efficient mixing scheme also based on Pulay's scheme. A special ``metric'' and a special ``preconditioning'' optimized for a plane-wave basis set will be introduced. Scaling of the method will be discussed in detail for non-self-consistent and self-consistent calculations. It will be shown that the number of iterations required to obtain a specific precision is almost independent of the system size. Altogether an order ${\mathit{N}}_{\mathrm{atoms}}^{2}$ scaling is found for systems containing up to 1000 electrons. If we take into account that the number of k points can be decreased linearly with the system size, the overall scaling can approach ${\mathit{N}}_{\mathrm{atoms}}$. We have implemented these algorithms within a powerful package called VASP (Vienna ab initio simulation package). The program and the techniques have been used successfully for a large number of different systems (liquid and amorphous semiconductors, liquid simple and transition metals, metallic and semiconducting surfaces, phonons in simple metals, transition metals, and semiconductors) and turned out to be very reliable. \textcopyright{} 1996 The American Physical Society.
81,985 citations
TL;DR: In this paper, the formal relationship between US Vanderbilt-type pseudopotentials and Blochl's projector augmented wave (PAW) method is derived and the Hamilton operator, the forces, and the stress tensor are derived for this modified PAW functional.
Abstract: The formal relationship between ultrasoft (US) Vanderbilt-type pseudopotentials and Bl\"ochl's projector augmented wave (PAW) method is derived. It is shown that the total energy functional for US pseudopotentials can be obtained by linearization of two terms in a slightly modified PAW total energy functional. The Hamilton operator, the forces, and the stress tensor are derived for this modified PAW functional. A simple way to implement the PAW method in existing plane-wave codes supporting US pseudopotentials is pointed out. In addition, critical tests are presented to compare the accuracy and efficiency of the PAW and the US pseudopotential method with relaxed core all electron methods. These tests include small molecules $({\mathrm{H}}_{2}{,\mathrm{}\mathrm{H}}_{2}{\mathrm{O},\mathrm{}\mathrm{Li}}_{2}{,\mathrm{}\mathrm{N}}_{2}{,\mathrm{}\mathrm{F}}_{2}{,\mathrm{}\mathrm{BF}}_{3}{,\mathrm{}\mathrm{SiF}}_{4})$ and several bulk systems (diamond, Si, V, Li, Ca, ${\mathrm{CaF}}_{2},$ Fe, Co, Ni). Particular attention is paid to the bulk properties and magnetic energies of Fe, Co, and Ni.
57,691 citations
TL;DR: An improved way of estimating the local tangent in the nudged elastic band method for finding minimum energy paths is presented, and examples given where a complementary method, the dimer method, is used to efficiently converge to the saddle point.
Abstract: An improved way of estimating the local tangent in the nudged elastic band method for finding minimum energy paths is presented. In systems where the force along the minimum energy path is large compared to the restoring force perpendicular to the path and when many images of the system are included in the elastic band, kinks can develop and prevent the band from converging to the minimum energy path. We show how the kinks arise and present an improved way of estimating the local tangent which solves the problem. The task of finding an accurate energy and configuration for the saddle point is also discussed and examples given where a complementary method, the dimer method, is used to efficiently converge to the saddle point. Both methods only require the first derivative of the energy and can, therefore, easily be applied in plane wave based density-functional theory calculations. Examples are given from studies of the exchange diffusion mechanism in a Si crystal, Al addimer formation on the Al(100) surfa...
6,825 citations
TL;DR: This paper describes how accurate off-lattice ascent paths can be represented with respect to the grid points, and maintains the efficient linear scaling of an earlier version of the algorithm, and eliminates a tendency for the Bader surfaces to be aligned along the grid directions.
Abstract: A computational method for partitioning a charge density grid into Bader volumes is presented which is efficient, robust, and scales linearly with the number of grid points. The partitioning algorithm follows the steepest ascent paths along the charge density gradient from grid point to grid point until a charge density maximum is reached. In this paper, we describe how accurate off-lattice ascent paths can be represented with respect to the grid points. This improvement maintains the efficient linear scaling of an earlier version of the algorithm, and eliminates a tendency for the Bader surfaces to be aligned along the grid directions. As the algorithm assigns grid points to charge density maxima, subsequent paths are terminated when they reach previously assigned grid points. It is this grid-based approach which gives the algorithm its efficiency, and allows for the analysis of the large grids generated from plane-wave-based density functional theory calculations.
5,417 citations