scispace - formally typeset
Search or ask a question
Journal ArticleDOI

CO2 electroreduction to ethylene via hydroxide-mediated copper catalysis at an abrupt interface

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
Journal ArticleDOI
01 Jul 2019
TL;DR: In this article, the authors describe progress and identify mechanistic questions and performance metrics for catalysts that can enable carbon-neutral renewable energy storage and utilization, and discuss design principles for improved activity and selectivity.
Abstract: Electrochemical carbon dioxide recycling provides an attractive approach to synthesizing fuels and chemical feedstocks using renewable energy. On the path to deploying this technology, basic and applied scientific hurdles remain. Integrating catalytic design with mechanistic understanding yields scientific insights and progresses the technology towards industrial relevance. Catalysts must be able to generate valuable carbon-based products with better selectivity, lower overpotentials and improved current densities with extended operation. Here, we describe progress and identify mechanistic questions and performance metrics for catalysts that can enable carbon-neutral renewable energy storage and utilization. Electrochemical carbon dioxide reduction is an attractive approach for obtaining fuels and chemical feedstocks using renewable energy. In this Review, the authors describe progress so far, identify mechanistic questions and performance metrics, and discuss design principles for improved activity and selectivity.

657 citations

Journal ArticleDOI
14 May 2020-Nature
TL;DR: Computational studies suggest that the Cu-Al alloys provide multiple sites and surface orientations with near-optimal CO binding for both efficient and selective CO2 reduction, and in situ X-ray absorption measurements reveal that Cu and Al enable a favourable Cu coordination environment that enhances C–C dimerization.
Abstract: The rapid increase in global energy demand and the need to replace carbon dioxide (CO2)-emitting fossil fuels with renewable sources have driven interest in chemical storage of intermittent solar and wind energy1,2. Particularly attractive is the electrochemical reduction of CO2 to chemical feedstocks, which uses both CO2 and renewable energy3–8. Copper has been the predominant electrocatalyst for this reaction when aiming for more valuable multi-carbon products9–16, and process improvements have been particularly notable when targeting ethylene. However, the energy efficiency and productivity (current density) achieved so far still fall below the values required to produce ethylene at cost-competitive prices. Here we describe Cu-Al electrocatalysts, identified using density functional theory calculations in combination with active machine learning, that efficiently reduce CO2 to ethylene with the highest Faradaic efficiency reported so far. This Faradaic efficiency of over 80 per cent (compared to about 66 per cent for pure Cu) is achieved at a current density of 400 milliamperes per square centimetre (at 1.5 volts versus a reversible hydrogen electrode) and a cathodic-side (half-cell) ethylene power conversion efficiency of 55 ± 2 per cent at 150 milliamperes per square centimetre. We perform computational studies that suggest that the Cu-Al alloys provide multiple sites and surface orientations with near-optimal CO binding for both efficient and selective CO2 reduction17. Furthermore, in situ X-ray absorption measurements reveal that Cu and Al enable a favourable Cu coordination environment that enhances C–C dimerization. These findings illustrate the value of computation and machine learning in guiding the experimental exploration of multi-metallic systems that go beyond the limitations of conventional single-metal electrocatalysts. Machine learning predicts Cu-Al electrocatalysts provide better efficiency and productivity than copper when using intermittent renewable electricity to convert carbon dioxide to useful chemicals and fuels.

656 citations

Journal ArticleDOI
TL;DR: A techno-economic analysis is presented with the goal of identifying maximally profitable products and the performance targets that must be met to ensure economic viability-metrics that include current density, Faradaic efficiency, energy efficiency, and stability.
Abstract: The electrochemical reduction of CO2 is a promising route to convert intermittent renewable energy to storable fuels and valuable chemical feedstocks. To scale this technology for industrial implementation, a deepened understanding of how the CO2 reduction reaction (CO2 RR) proceeds will help converge on optimal operating parameters. Here, a techno-economic analysis is presented with the goal of identifying maximally profitable products and the performance targets that must be met to ensure economic viability-metrics that include current density, Faradaic efficiency, energy efficiency, and stability. The latest computational understanding of the CO2 RR is discussed along with how this can contribute to the rational design of efficient, selective, and stable electrocatalysts. Catalyst materials are classified according to their selectivity for products of interest and their potential to achieve performance targets is assessed. The recent progress and opportunities in system design for CO2 electroreduction are described. To conclude, the remaining technological challenges are highlighted, suggesting full-cell energy efficiency as a guiding performance metric for industrial impact.

599 citations

Journal ArticleDOI
01 Jun 2020
TL;DR: In this article, a fluorine-modified copper catalyst was proposed for electrocatalytic CO2 electroreduction in a flow cell, achieving an ultrahigh current density of 1.6
Abstract: Electrocatalytic reduction of CO2 into multicarbon (C2+) products is a highly attractive route for CO2 utilization; however, the yield of C2+ products remains low because of the limited C2+ selectivity at high CO2 conversion rates. Here we report a fluorine-modified copper catalyst that exhibits an ultrahigh current density of 1.6 A cm−2 with a C2+ (mainly ethylene and ethanol) Faradaic efficiency of 80% for electrocatalytic CO2 reduction in a flow cell. The C2–4 selectivity reaches 85.8% at a single-pass yield of 16.5%. We show a hydrogen-assisted C–C coupling mechanism between adsorbed CHO intermediates for C2+ formation. Fluorine enhances water activation, CO adsorption and hydrogenation of adsorbed CO to CHO intermediate that can readily undergo coupling. Our findings offer an opportunity to design highly active and selective CO2 electroreduction catalysts with potential for practical application. Electrocatalytic reduction of CO2 into multicarbon (C2+) products is a highly attractive route for CO2 utilization. Now, a fluorine-modified copper catalyst is shown to achieve current densities of 1.6 A cm−2 with a C2+ Faradaic efficiency of 80% for electrocatalytic CO2 reduction in a flow cell.

591 citations

Journal ArticleDOI
TL;DR: In this article, the impact of reaction rate on catalytic behavior and the operation of gas-diffusion layers for electrocatalytic CO2 reduction is discussed. But, the authors focus on high current density (∼200 mA cm−2) and do not evaluate the performance of these catalysts.
Abstract: Electrocatalytic CO2 reduction has the dual-promise of neutralizing carbon emissions in the near future, while providing a long-term pathway to create energy-dense chemicals and fuels from atmospheric CO2. The field has advanced immensely in recent years, taking significant strides towards commercial realization. Catalyst innovations have played a pivotal role in these advances, with a steady stream of new catalysts providing gains in CO2 conversion efficiencies and selectivities of both C1 and C2 products. Comparatively few of these catalysts have been tested at commercially-relevant current densities (∼200 mA cm−2) due to transport limitations in traditional testing configurations and a research focus on fundamental catalyst kinetics, which are measured at substantially lower current densities. A catalyst's selectivity and activity, however, have been shown to be highly sensitive to the local reaction environment, which changes drastically as a function of reaction rate. As a consequence of this, the surface properties of many CO2 reduction catalysts risk being optimized for the wrong operating conditions. The goal of this perspective is to communicate the substantial impact of reaction rate on catalytic behaviour and the operation of gas-diffusion layers for the CO2 reduction reaction. In brief, this work motivates high current density catalyst testing as a necessary step to properly evaluate materials for electrochemical CO2 reduction, and to accelerate the technology toward its envisioned application of neutralizing CO2 emissions on a global scale.

581 citations

References
More filters
Journal ArticleDOI
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

Journal ArticleDOI
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

Journal ArticleDOI
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

Journal ArticleDOI
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

Journal ArticleDOI
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