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: In this article, a PGM-free nickel and nitrogen-doped porous carbon catalyst (Ni-N-C) was proposed to replace PGM catalysts in CO2-to-CO electrolyzers.
Abstract: The electrochemical CO2 reduction reaction (CO2RR) to pure CO streams in electrolyzer devices is poised to be the most likely process for near-term commercialization and deployment in the polymer industry. The reduction of CO2 to CO is electrocatalyzed under alkaline conditions on precious group metal (PGM) catalysts, such as silver and gold, limiting widespread application due to high cost. Here, we report on an interesting alternative, a PGM-free nickel and nitrogen-doped porous carbon catalyst (Ni–N–C), the catalytic performance of which rivals or exceeds those of the state-of-the-art electrocatalysts under industrial electrolysis conditions. We started from small scale CO2-saturated liquid electrolyte H-cell screening tests and moved to larger-scale CO2 electrolyzer cells, where the catalysts were deployed as Gas Diffusion Electrodes (GDEs) to create a reactive three-phase interface. We compared the faradaic CO yields and CO partial current densities of Ni–N–C catalysts to those of a Ag-based benchmark, and its Fe-functionalized Fe–N–C analogue under ambient pressures, temperatures and neutral pH bicarbonate flows. Prolonged electrolyzer tests were conducted at industrial current densities of up to 700 mA cm−2. Ni–N–C electrodes are demonstrated to provide CO partial current densities above 200 mA cm−2 and stable faradaic CO efficiencies around 85% for up to 20 hours (at 200 mA cm−2), unlike their Ag benchmarks. Density functional theory-based calculations of catalytic reaction pathways help offer a molecular mechanistic basis of the observed selectivity trends on Ag and M–N–C catalysts. Computations lend much support to our experimental hypothesis as to the critical role of N-coordinated metal ion, Ni–Nx, motifs as the catalytic active sites for CO formation. Apart from being cost effective, the Ni–N–C powder catalysts allow flexible operation under acidic, neutral, and alkaline conditions. This study demonstrates the potential of Ni–N–C and possibly other members of the M–N–C materials family to replace PGM catalysts in CO2-to-CO electrolyzers.
300 citations
TL;DR: A high-yield, flexible and self-supported single-atom nickel-decorated porous carbon membrane catalyst that possesses interconnected nanofibers and hierarchical pores, affording abundant effective nickel single atoms that participate in carbon dioxide reduction.
Abstract: Carbon dioxide electroreduction provides a useful source of carbon monoxide, but comparatively few catalysts could be sustained at current densities of industry level. Herein, we construct a high-yield, flexible and self-supported single-atom nickel-decorated porous carbon membrane catalyst. This membrane possesses interconnected nanofibers and hierarchical pores, affording abundant effective nickel single atoms that participate in carbon dioxide reduction. Moreover, the excellent mechanical strength and well-distributed nickel atoms of this membrane combines gas-diffusion and catalyst layers into one architecture. This integrated membrane could be directly used as a gas diffusion electrode to establish an extremely stable three-phase interface for high-performance carbon dioxide electroreduction, producing carbon monoxide with a 308.4 mA cm−2 partial current density and 88% Faradaic efficiency for up to 120 h. We hope this work will provide guidance for the design and application of carbon dioxide electro-catalysts at the potential industrial scale. Here the authors deploy Ni single atom-decorated carbon membranes as integrated gas diffusion electrodes to construct an extremely stable three-phase interface for CO2 electroreduction, producing CO with a partial current density of 308.4 mA cm–2 and a Faradaic efficiency of 88% for up to 120 h.
294 citations
TL;DR: It is reported that catalysts with nanocavities can confine carbon intermediates formed in situ, which in turn covers the local catalyst surface and thereby stabilizes Cu+ species, which leads to the marked C2+ selectivity at large conversion rate.
Abstract: Selective and efficient catalytic conversion of carbon dioxide (CO2) into value-added fuels and feedstocks provides an ideal avenue to high-density renewable energy storage. An impediment to enabling deep CO2 reduction to oxygenates and hydrocarbons (e.g., C2+ compounds) is the difficulty of coupling carbon-carbon bonds efficiently. Copper in the +1 oxidation state has been thought to be active for catalyzing C2+ formation, whereas it is prone to being reduced to Cu0 at cathodic potentials. Here we report that catalysts with nanocavities can confine carbon intermediates formed in situ, which in turn covers the local catalyst surface and thereby stabilizes Cu+ species. Experimental measurements on multihollow cuprous oxide catalyst exhibit a C2+ Faradaic efficiency of 75.2 ± 2.7% at a C2+ partial current density of 267 ± 13 mA cm-2 and a large C2+-to-C1 ratio of ∼7.2. Operando Raman spectra, in conjunction with X-ray absorption studies, confirm that Cu+ species in the as-designed catalyst are well retained during CO2 reduction, which leads to the marked C2+ selectivity at a large conversion rate.
289 citations
TL;DR: Zhang et al. as discussed by the authors designed a series of nickel phthalocyanine molecules supported on carbon nanotubes as molecularly dispersed electrocatalysts (MDEs), achieving CO2 reduction performances that are superior to aggregated molecular catalysts in terms of stability, activity and selectivity.
Abstract: Electrochemical reduction of CO2 is a promising route for sustainable production of fuels. A grand challenge is developing low-cost and efficient electrocatalysts that can enable rapid conversion with high product selectivity. Here we design a series of nickel phthalocyanine molecules supported on carbon nanotubes as molecularly dispersed electrocatalysts (MDEs), achieving CO2 reduction performances that are superior to aggregated molecular catalysts in terms of stability, activity and selectivity. The optimized MDE with methoxy group functionalization solves the stability issue of the original nickel phthalocyanine catalyst and catalyses the conversion of CO2 to CO with >99.5% selectivity at high current densities of up to −300 mA cm−2 in a gas diffusion electrode device with stable operation at −150 mA cm−2 for 40 h. The well-defined active sites of MDEs also facilitate the in-depth mechanistic understandings from in situ/operando X-ray absorption spectroscopy and theoretical calculations on structural factors that affect electrocatalytic performance. Widespread deployment of electrochemical CO2 reduction requires low-cost catalysts that perform well at high current densities. Zhang et al. show that methoxy-functionalized nickel phthalocyanine molecules on carbon nanotubes can operate as high-performing molecularly dispersed electrocatalysts at current densities of up to −300 mA cm–2.
287 citations
TL;DR: In this article, an electrocatalyst with confined reaction volume by coating Cu catalysts with nitrogen-doped carbon layers was developed, which achieved an ethanol Faradaic efficiency of (52 ± 1)% and an ethanol cathodic energy efficiency of 31%.
Abstract: The carbon dioxide electroreduction reaction (CO2RR) provides ways to produce ethanol but its Faradaic efficiency could be further improved, especially in CO2RR studies reported at a total current density exceeding 10 mA cm−2. Here we report a class of catalysts that achieve an ethanol Faradaic efficiency of (52 ± 1)% and an ethanol cathodic energy efficiency of 31%. We exploit the fact that suppression of the deoxygenation of the intermediate HOCCH* to ethylene promotes ethanol production, and hence that confinement using capping layers having strong electron-donating ability on active catalysts promotes C–C coupling and increases the reaction energy of HOCCH* deoxygenation. Thus, we have developed an electrocatalyst with confined reaction volume by coating Cu catalysts with nitrogen-doped carbon. Spectroscopy suggests that the strong electron-donating ability and confinement of the nitrogen-doped carbon layers leads to the observed pronounced selectivity towards ethanol. The electroreduction of CO2 to ethanol could enable the clean production of fuels using renewable power. This study shows how confinement effects from nitrogen-doped carbon layers on copper catalysts enable selective ethanol production from CO2 with a Faradaic efficiency of up to 52%.
286 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