Next Article in Journal
Comparative Studies of Fischer-Tropsch Synthesis on Iron Catalysts Supported on Al2O3-Cr2O3 (2:1), Multi-Walled Carbon Nanotubes or BEA Zeolite Systems
Previous Article in Journal
Catalytic Oxidation of Dimethyl Disulfide over Bimetallic Cu–Au and Pt–Au Catalysts Supported on γ-Al2O3, CeO2, and CeO2–Al2O3
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Graphene Supported Tungsten Carbide as Catalyst for Electrochemical Reduction of CO2

by
Sahithi Ananthaneni
,
Zachery Smith
and
Rees B. Rankin
*
Department of Chemical Engineering, Villanova University, 800 E Lancaster University, Villanova, PA 19085, USA
*
Author to whom correspondence should be addressed.
Catalysts 2019, 9(7), 604; https://doi.org/10.3390/catal9070604
Submission received: 14 June 2019 / Revised: 28 June 2019 / Accepted: 4 July 2019 / Published: 15 July 2019

Abstract

:
Electrochemical reduction of CO2 to useful chemical and fuels in an energy efficient way is currently an expensive and inefficient process. Recently, low-cost transition metal-carbides (TMCs) have been proven to exhibit similar electronic structure similarities to Platinum-Group-Metal (PGM) catalysts and hence, can be good substitutes for some important reduction reactions. In this work, we test graphene-supported WC (Tungsten Carbide) nanoclusters as an electrocatalyst for the CO2 reduction reaction. Specifically, we perform density functional theory (DFT) studies to understand various possible reaction mechanisms and determine the lowest thermodynamic energy landscape of CO2 reduction to various products, such as CO, HCOOH, CH3OH, and CH4. This in-depth study of reaction energetics could lead to improvements and development of more efficient electrocatalysts for CO2 reduction.

Graphical Abstract

1. Introduction

The carbon dioxide (CO2) is very stable under environmental conditions, and reduction to some hydrocarbon products is an endothermic (or endergonic) process. The process requires a moderate to highly negative potential combined with excellent catalyst kinetic barrier efficiency to minimize the combined energy barriers for this reduction reaction. Among various methods of CO2 reduction, photocatalytic and electrocatalytic processes are predominantly studied. In the case of electrocatalysis, the reaction rate and products can be controlled by tuning the external potential, whereas the driving force for photocatalysis is solar energy. Although this work does not focus on studying photocatalytic reduction of CO2, it is well worth understanding this process and reduction mechanism while discussing electrochemical CO2 reduction and its state-of-the-art catalysts.
Photocatalysis involves photon absorption by a semiconductor photocatalyst followed by electron-hole pair excitation, charge transfer, and surface (redox) chemistry. The surface reaction pathways in CO2 photocatalytic reduction are similar to electrocatalytic reduction. However, the oxidation and reduction reactions occur on different sides of the catalyst, and the driving force for this process is from solar energy. It still remains as a significant challenge to fabricate a semiconductor photocatalyst which can efficiently separate the electron-hole pairs generated and prevent recombination, transfer electrons to the surface, and provide a catalytic surface for the reaction to occur [1,2,3]. There have been continuous improvements in design and engineering towards the development of a novel photocatalyst in various ways hoping to overcome the above challenges. Some of them are by:
  • Adjusting bandgap and band position to capture solar energy effectively, thereby improving the efficiency of CO2 reduction. For example, the potential of the valence band of various metal oxide catalysts is around 3 eV vs. NHE (normal hydrogen electrode), utilizing only ultraviolet light. Replacing the O atom with N or C will narrow the bandgap thereby efficiently absorbing the solar radiations [4,5].
  • Modifying surface structure and electronic properties by addition of oxygen vacancies, which can lower the energy barrier due to the attraction between CO2 molecules and oxygen vacancy. DFT calculations on TiO2 surface with oxygen vacancies is more effective for CO2 activation than TiO2 [6,7].
  • The use of nanostructured semiconductors in the form of nanorods, nanowires, nanotubes, nanobelts, etc. decreases the electron-hole recombination rate as these are usually in a single crystalline phase that eliminates the possibility of grain boundaries and defects in the materials as they act as recombination sites for the electron-hole pair in polycrystalline materials. At the same time, one-dimensional nanostructures improve electron transport by improving the separation of electron-hole pairs. For instance, Zn2GeO4 nanoribbons are proven to show improved photocatalytic activity towards CH4 formation when compared to bulk Zn2GeO4 [8,9]
  • Use of co-catalysts in the process to promote the separation and movement of charge carriers. This will minimize the recombination of electron-hole pairs due to the barrier between semiconductor and co-catalyst. Metal nanoparticles, such as Pt, Rh, Pd, Cu, Ag, Au, supported on semiconductor are proven to work efficiently when compared to pure semiconducting photocatalysts [10,11].
Several electrocatalysts are being studied to convert CO2 to useful chemicals and fuels, but they are either limited by high overpotentials or poor product selectivity. In 1985, Hori et al. reported that Cu is a unique metal catalyst which can reduce CO2 to hydrocarbons efficiently, and further studies described that a Cu catalyst could electrochemically reduce CO2 to 16 different products among which CH4 and C2H4 showed higher current densities but at larger overpotentials of up to 1V [12,13,14]. To overcome this obstacle and make the CO2 reduction reaction more viable, we require more complex and tailored materials than simple, pure transition metals.
Recently, low-cost transition metal-carbides (TMCs) have received special attention as electrocatalysts as they have shown improved catalyst stability, activity, and selectivity when compared to their parent metal/elements [15,16]. Formation of carbides (with C in the crystal lattice or nanoparticle surface) modifies the metal–atom bonding, increasing the metal–metal distance, thereby causing a contraction in the metal atoms’d-band. These d-band contractions would give a better density of states near the Fermi level than their parent metal. In particular, non-noble metal-carbides, such as Mo and W, display a similar electronic structure to noble metals due to their metal-covalent binding and hence, can be significant substitutes for precious catalysts in various important catalytic reactions [16]. Additional to electronic and magnetic properties similar to transition metals, they exhibit high melting points as in ionic compounds and hardness similar to covalent solids. Therefore, the bonding in TMCs can be explained as a combination of metallic, covalent, and ionic components [17].
Theoretical and experimental studies proved that TMCs show better catalyst activity in“hydrogen-participating” reactions [18]. For example, metal terminated WC exhibited good activity for hydrogen evolution reaction (HER) and WC coated with Fe when tested for CO2 reduction are selective towards methane formation. This is because electronic properties of W atoms surrounding Fe are modified, thereby modifying the selectivity. Another interesting phenomenon of TMCs is they do not follow scaling relations that correlate binding energies of intermediates in a reaction network as with other transition metals mainly because of their oxygen affinity, i.e., their tendency to bind carbon-bound species weakly compared to oxygen-bound species. This will open up several possibilities to improve their catalyst activity, selectivity compared to metal catalysts [19,20]. In the same way, in the density functional theory (DFT) study of Fe and Co carbides for Fischer–Tropsch synthesis, it was found that FeC (iron carbide) is more active than pure Fe for CO hydrogenation and methane selectivity is higher on Co carbide compared to pure Co [21].
Among different metal-carbides, tungsten carbide (WC) based compounds are widely studied electrocatalysts [22]. They are investigated in various forms, such as:
  • Alloys to combine the electronic properties of WC with other metal(s), for example, Tantalum doped WC displayed better activity towards hydrogen evolution (HER) when compared to unmodified WC [23].
  • The specific shape and structural arrangement, such as core-–shell structures of WC with monolayer metal coatings which are stable against CO poisoning, thereby improving the activity of methanol electro-oxidation [24,25,26,27].
  • Catalyst support to increase the electrocatalytic activity leading to better performance of fuel cells, such as WC supported Pt which is found to be more thermally and electrochemically stable than Pt/C for oxygen reduction reaction [28,29].
  • Co-catalyst to the catalytic system where strong electronic interactions between them might modify (maximize) the electrocatalytic activity. As an example, Ni with WC nanocluster for urea electro-oxidation showed high tolerance towards CO poisoning, and high stability thereby enhancing catalyst activity [30].
All these studies explain that TMCs have the potential to work as better catalysts making them an attractive alternative for traditional metallic catalysts in some of the industrially relevant catalytic reactions. Nanocatalysts in the form of nanoclusters where the atoms are structurally quasi-defined to well-defined help in experimental and theoretical investigations of important electronic structure properties in CO2 electroreduction reactions (CO2RR). Additionally, graphene as a catalyst support further improves the active surface area for the catalyst systems by providing a minimal adsorption footprint for the TMC nanoparticles (NP). Other unique properties of this two-dimensional structure, such as high stability and electrical conductivity which can modify the TMC NP electronic structure, help in selectivity and cost reduction of catalysts playing a pivotal role in most of the heterogeneous catalyst systems [31].
Motivated by all the above aspects, in this study, we test graphene supported WC nanocluster as an electrocatalyst for CO2RR. This work determines the performance of WC/graphene as a catalyst system for CO2 reduction to various products, such as CO, HCOOH, CH3OH, and CH4, which possibly could be the platform for designing new and improved TMC electrocatalysts for this important reduction reaction. To achieve this, we have calculated binding free energies of all possible reaction intermediates and analyzed the reaction mechanisms in detail by focusing on understanding the effect solvation energies on product selectivity and catalyst activity and lastly determined the lowest energy pathways for all the products mentioned above. We inferred that CH4 is favored over CH3OH on WC/graphene at lower reducing potentials. From literature, CH3OH formation is less preferred on WC and metal-coated WC because of the stronger binding of O* and OH* bound species [20,27,32]. Results from the current work explain that the reduction to CH3OH on graphene supported WC is thermodynamically favorable at higher negative potentials compared to reduction to CH4, which is in qualitative agreement with the results from the literature.

2. Results and Discussion

We studied several possible reaction intermediates and plotted free energy diagrams for various possible reaction mechanisms to determine the lowest energy pathway for each product. The stability of these adsorbed species can be determined in terms of the difference in binding free energies between intermediates formed along the reaction pathway and reactants. That is, when the adsorption of an intermediate is unstable (ΔG is positive), the downhill reaction is possible by applying an external potential. Intermediate species with net positive ΔG formation in the current study are COH*, C*, CH*, CHOH*. At this point, it cannot be decoupled as to how much of this energy is purely attributed to the adsorbed intermediate and how much to the catalyst system. However, as we did not observe any significant changes to the catalyst morphology after adsorption of these intermediates, we assume that contribution to the net positive ΔG formation of adsorbed intermediates is mostly due to the net instability or metastability of the adsorbate/intermediate itself relative to the reactant. We expect to fully address this area in the follow-up work when we can locate experimental collaboration to further validate the optimization of this catalytic system.
As we apply higher (larger magnitude) negative voltages to different reaction pathways, the pathway with the smallest positive rate-limiting step will be the first pathway to become entirely exergonic across all steps in the pathway. This pathway will be the best lowest energy reaction pathway for CO2RR. These free energy diagrams (FEDs) provide an overall understanding of the reaction mechanism, electrocatalytic activity for CO2RR as well as selectivity of products. Here, we focus on the lowest ∆G pathway for CO2RR to CO, HCOOH, CH3OH, and CH4.

2.1. CO and HCOOH as Products

Figure 1 and Figure 2 show the lowest energy pathway for CO2 reduction to HCOOH and CO at 0 V vs. RHE (reversible hydrogen electrode) on graphene supported WC. These products are obtained by two proton–electron transfers along the reaction pathway. CO2 is first protonated to form either COOH* or OCHO*. The limiting potential for HCOOH and CO formation depends on how strongly or weakly COOH* and OCHO* bind to the surface. When we compare binding free energies, the formation of COOH* is less exergonic compared to OCHO*. A second proton–electron transfer results in the formation of HCOOH and CO. Although OCHO* is more stable compared to COOH*, the pathway via COOH* will minimize the rate-limiting step by around 1.9 V for HCOOH formation and 2.5 V for CO formation. This larger difference in binding free energies of COOH* and OCHO* is due to the fact that COOH tries to bind to the surface of the catalyst via C and O atom and OCHO binds to the surface of the catalyst via two oxygen atoms and WC has strong oxygen affinity, i.e., it binds O* and OH* bound intermediate species very strongly, creating significant energy barriers.
Figure 1 shows the lowest energy pathway for HCOOH formation. The calculated limiting potential for HCOOH formation is −2.45 V, and the potential determining step, i.e., the rate-limiting step is the formation of HCOOH from COOH*. Figure 2 shows the lowest energy pathway for CO formation. The calculated limiting potential for CO formation is −1.36 V, and the rate-limiting step is the formation of CO from COOH*. This is in agreement with previous experimental and theoretical studies showing the reduction potential for CO2 to CO is in the range of −0.72 V to −1.5 V on pure metal surfaces, such as Ag, Au, Zn. This comparison also confirms that WC/graphene can be used as an alternative for precious electrocatalysts for CO2 reduction to CO. Further reduction of CO* to CH3OH and CH4 as products is determined by the binding energy of CO*. Stronger binding of CO* results in hydrogen evolution due to CO poisoning and weaker binding of CO* results in CO desorption before further reduction to products. It is proven from previous work that metals, such as Ag, Au, and Zn, that bind CO weakly during CO2 reduction can further reduce to CH3OH and CH4 but at higher negative potentials. Therefore, in the next section, we will discuss CO2 reduction to CH3OH and CH4 on graphene supported WC.

2.2. CH3OH and CH4 as Products

This section discusses the lowest energy pathways for CO2 reduction to CH3OH and CH4 at 0 V on graphene supported WC. These products are obtained by six and eight proton–electron transfers, respectively, along the reaction pathway. The first two steps up to CO* formation are described in detail in the previous section. However, complexity arises when CO* is further protonated. As the number of protons–electrons transferred increase, the number of required intermediates in each reaction network increases. This is one of the reasons why the CO2 electrochemical reduction reaction mechanism and its thermodynamics is more complex to study when compared to oxygen reduction reaction (ORR), hydrogen evolution reaction (HER), and other reactions involving the transfer of fewer proton–electron pairs. Similar to CO and HCOOH formation, CH3OH and CH4 product formation follow a reaction pathway via COOH*. The protonation of CO* in electrochemical reduction to CH3OH and CH4 can follow either an oxophilic pathway (CO* is protonated at C atom to form HCO* and binds to the surface of the catalyst via O atom) or a carbophilic pathway (CO* is protonated at O atom to form COH* and binds to the surface of the catalyst via C atom).

2.3. CH3OH as Product

Figure 3 and Figure 4 display the lowest energy pathway for CH3OH formation with and without taking solvation energies into account. Apart from the reaction mechanism, these free energy diagrams also guide us in understanding the effect of solvation on the reaction thermodynamics in CH3OH formation. In Figure 3, we have shown the lowest energy pathway for CH3OH formation when solvation energies are excluded. Real aqueous phase electrochemical CO2RR has water molecules present around the TMC NP, and these water molecules can stabilize the reactants, products, and intermediates through H-bonding. The protonation of CO* in this product formation is via the oxophilic pathway, i.e., via HCO* species. If the solvation effect is not considered, the calculated limiting potential for the product formation is −0.56 V with the rate-limiting step of CH3O* to CH3OH formation. Figure 4 shows the lowest energy pathway with the addition of solvation energies. Inclusion of solvation energies in electrochemical reduction to CH3OH modified not only the reducing potential but also the lowest energy reaction pathway and the rate-limiting step. The limiting potential is increased from −0.56 V to −1.79 V, and the reaction pathway is shifted from an oxophilic to the carbophilic pathway, i.e., the pathway is shifted from HCO* to COH*. The new rate-limiting step is the formation of CHOH* from COH*. It is worth noting that these changes in the reaction pathway, rate-limiting step, and its corresponding potential is due to the following reason: WC has strong oxygen affinity; therefore, it strongly binds all O* and OH* bound intermediate species. Inclusion of the solvent effect has further stabilized these O*, and OH* bound intermediate species altering the electronic binding free energies, thereby creating huge energy barriers in the reaction network.
For additional comparison, in Figure 4, we have included the carbophilic pathway without solvation energies (red pathway). The calculated limiting potential for the product formation is −1.73 V with rate-limiting step COH* to CHOH*. Inclusion of solvation energies (blue pathway) has stabilized the adsorbate species and minimized the reducing potential by ~0.06 V without modifying the rate-limiting step of the carbophilic pathway.

2.4. CH4 as Product

Figure 5, red pathway, shows the lowest energy pathway without the solvation energy correction factor and the calculated limiting potential for the CH4 formation is −1.58 V with the same rate-limiting step, i.e., CO* to COH*. The blue pathway is the lowest energy pathway after adding the solvation correction factor and the calculated limiting potential for product formation is −0.84 V, and the rate-limiting step is the formation of CO* from COH*. The key point here is the inclusion of solvation energies minimized the reducing potential by around 0.75V but did not modify the lowest energy reaction pathway. This may be due to the fact that CH4 formation pathway is via the carbophilic pathway and carbon bound species. Therefore, the effect of solvation stabilized the COOH*, CO*, and COH* binding free energies modifying just the first four steps in the reaction pathway.
Figure 6 and Figure 7 show the calculated lowest energy profiles, including solvation energies for a complete series of elementary steps leading to CH3OH and CH4 formation, respectively. We have also included the free energy profile of product formation at the thermodynamic limiting potential/reducing potential, i.e., the potential at which all the elementary steps are downhill (exergonic) in free energy. This is used to determine the overpotentials of the reaction on a particular electrocatalyst.
Overpotentials can be directly related to the catalyst activity and energy efficiency and can be obtained by the difference between the equilibrium potential and limiting potential. On our catalyst system, CO2 can reduce to CH3OH with an overpotential of 1.81 V and CH4 with an overpotential of 0.67 V. For additional comparison we have shown free energy profile at the equilibrium potentials, and the maximum potential allowed by the thermodynamics.

2.5. Comparison of CO2 Reduction to CH4 and CH3OH on Graphene Supported WC Nanocluster and WC (0001)

In this section, we compare CO2 reduction on WC/ graphene to that on WC (0001). Our analysis and comparison are based on the reaction free energy profile. From the work of Wannakao S. et al., the limiting potential for CO2 reduction to CH4 on WC (0001) is −0.35 V (no solvation) which is 0.5 V (with solvation) lower than the potential achieved in the current work. Similarly, the limiting potential for CH3OH formation on WC (0001) is around −0.39 V (no solvation), which is around 1.5 V (with solvation) lower than the potential achieved in the current work. This difference in potential could be due to the following reasons:
  • Employing different functionals in the DFT calculations would result in differences in binding free energies of intermediates. In the work of Wannakao S. et al., it is proven that
  • Another reason is predicted to be due to the coverage of adsorbed intermediates on the surface of the catalyst. Our catalyst system is designed by placing only one adsorbate species on a single nanocluster which is approximately equal to 1/6 or 1/9 monolayer (ML) coverage of adsorbate species on the surface (assuming each side of the cluster mimics 3 × 2 or 3 × 3 slab surface). However, in the work of Wannakao S. et al., WC (0001) was modeled by 1/6 to 1/9 ML coverage of adsorbate species which means these are placed comparatively closer than our adsorbate species arrangement. We predict that the influence of lateral interactions between adsorbed intermediates could also lead to a difference in reaction free energies. To investigate this, we vary (increase) the surface coverage of intermediate species by placing two species instead of one in our catalyst system as neighboring atoms/moieties near the active site. In other words, this system is modified to try to approximate the effects of 1/6 to 1/9 ML coverage of adsorbate species in their work. Figure 8 shows how the reaction free energies vary when the proportion of adsorbate coverage on the catalyst system varies. Here, we have computed the free energies of initial steps in the reaction network (COOH*, CO*, COH*) as these are the pathway determining intermediates. Co-adsorption of these species generated an upward shift of binding free energies of all the initial three steps. Consequently, the rate-limiting step in the case of co-adsorbed species shifted to the COOH* protonation step from the CO* protonation step. This analysis also explains that the surface coverage of the adsorbed species plays an equally important role in determining the energetics of CO2 reduction reaction.

2.6. Role of Graphene

In the WC-graphene system, the energy barrier of the rate-limiting step is 0.85 V, whereas the energy barrier of the rate-limiting step for the WC nanocluster is 1.04 V. In other words, CO* and COH* are further stabilized on graphene supported WC. Therefore, graphene as a support for WC catalyst enhances the energy efficiency of CO2 reduction reaction by lowering the limiting potential by ~0.2 V. Figure 9 compares the free energies of initial steps in the reaction network (COOH*, CO*, COH*) on graphene supported WC and plain WC nanocluster. Therefore, instead of a numerical comparison of binding free energies and reduction potentials on graphene supported WC and WC (0001), we have compared the reaction pathways and potential determining steps. We have graphically represented the reaction pathways of CH4 formation on WC/graphene and WC (0001) in Figure 10 to show the similarities between them.
From literature, CH4 formation on WC (0001) follows a reaction pathway via COOH* and COH* with CO* protonation to COH* as the potential determining step. Similarly, current work shows that the lowest energy pathway for CO2 reduction follows the carbophilic pathway with an identical rate-limiting step and binding site (binds to W through C from COH*). In addition, while comparing COH* and HCO* binding energies, both on WC/graphene and WC (0001), HCO* (via C and O atoms) binds more strongly than COH* (via C atom) and in both the cases’ pathway via COH* minimizes the reducing potential when compared to HCO*.

2.7. HER

There are various studies explaining that the catalysts for CO2 electrochemical reduction are also good catalysts for HER as thermodynamically both HER and CO2RR require approximately 0 V (RHE) for the reaction to take place and this is one of the major challenges in CO2RR. As HER is an unwanted side reaction, we need to design a catalyst which can suppress this reaction and further reduce CO2 to various products. Hence, in this study, we have considered understanding HER in parallel. This is a straightforward reaction where 2 proton electrons are required to form H2. This is a two-step mechanism, and the lowest energy reaction pathway at U = 0 V is shown in Figure 11. The limiting potential required for HER to take place is −0.98 V, which is higher than the limiting potential required to form CH4. A recent theoretical study explained that CO* formed along the CO2 reduction can occupy the adsorption sites preventing H* contact to the surface and eventually H2 formation [33].
To further understand if WC/graphene system favors HER or CO2RR or both, we have compared binding free energies of H* and CO* at a range of applied potentials (from 0 V to 1 V vs. RHE) as shown in Figure 12. At the range of potentials studied, CO* formation is thermodynamically more preferred than H* formation at higher negative potentials. In other words, CO* becomes more stable (higher negative binding free energy) than H* at higher negative potentials preventing HER and at the same time opening up pathways for further reduction to products with C–H bonds.

3. (Computational) Materials and Methods

We used density functional theory (DFT) to understand catalyst surface reaction pathways in detail. We performed plane wave DFT calculations with VASP (Vienna Ab Initio Simulation Package) to find optimized surface structures and calculate electronic structure properties, such as ground state energies, binding energies, charge densities, and perform Bader charge analysis [34,35,36,37,38]. Throughout this work, all the electronic structure calculations were performed using Van Der Waals, opt-PBE functional as they are proven to show high accuracy to study adsorption properties [38,39,40]. A Fermi smearing of 0.2 eV is used, and calculations were performed with gamma centered k-points mesh of 2 × 2 × 1 with a convergence of ground state energies up to 10−5eV/mol-unit cell with respect to k-point sampling [41]. A vacuum space of 12 Å was defined to minimize the interactions between repeated structures in the direction parallel to the surface normal of the graphene plane. All the reaction energy calculations were completed using the lowest energy conformation of the intermediate species. Detailed images of the structure and some adsorbate snapshots are provided in the Supporting Information (Figure S1) submitted with this manuscript. We took advantage of the computational hydrogen electrode (CHE) approach for screening and designing electrocatalysts primarily to understand reaction mechanisms for CO2 electroreduction to CH4 and CH3OH [42,43,44]. Possible reaction pathways for electrochemical reduction of CO2 to CO, CH4, and CH3OH are shown in Figure 13. We used RHE (reversible hydrogen electrode) as a reference electrode throughout this work. Since RHE is the reference for all the reactions, it can be set to zero.
The overall reactions of CO2 reduction to products and their corresponding thermodynamic equilibrium potentials are shown in Table 1 below.
The Binding energy of each intermediate species in the reaction network can be calculated using Equation (1). This is the difference between the DFT energy of the adsorbed intermediate and the sum of the DFT energy of the bare surface and the formation energies of C, H, O w.r.t gaseous CO, H2, and H2O, where x,y, and z denote stoichiometries of each element.
Δ E binding = E DFT ( E surface + x E C + y E H + z E O )
Similarly, the free energy of each electrochemical step in the reaction pathway corrected by zero-point energies (ZPE) with enthalpy and entropy contributions at potentials other than 0 V can be estimated using Equation (2) given below.
Δ G ( U ) = Δ E rxn + Δ ZPE + C p dT T Δ S neU
U is the applied potential to make the elementary step exergonic (known as the limiting potential, UL) and n is the number of proton–electron pairs consumed in each step. The ZPE, enthalpy, and entropy of adsorbed species are obtained from previously determined values as they are assumed to be largely independent of catalyst surface and therefore can be approximated to be the same for all structures [42,46]. Since an explicit treatment of a number of water molecules is tedious to carry out using DFT methods, we have used the implicit method in VASP with the default dielectric constant of H2O to calculate solvation energies. These are generally much less computationally demanding than explicit methods but can reproduce significant results as with explicit methods for O*, OH* bound intermediate species [47,48].

4. Conclusions and Future Work

This work provides theoretical evidence that graphene supported WC nanoparticles could be a useful catalyst system for CO2 reduction to light hydrocarbons and fuels. We have presented an improved understanding of CO2 reduction reaction mechanisms and provided the lowest energy pathway for various products. The results imply that CH4 is favored over CH3OH on this catalyst system because of the strong oxygen affinity towards WC. It is also worth noting that the solvation effect plays an important role in determining the reaction pathway. We have shown how the addition of the solvation effect has shifted the rate-limiting step from CH3OH formation step to the COH* protonation step on the lowest energy pathway for CO2 reduction to CH3OH. At the same time, we have also shown how the binding free energies are overestimated in the CH4 reduction reaction when the effect of solvation is not considered, thus resulting in higher limiting potentials. The addition of the solvation effect to the CH4 reduction pathway has minimized the limiting potential by ~0.75 V. In our catalyst system, CO2 can reduce to CH4 with an overpotential of 0.67 V and to CH3OH at higher negative potentials, i.e., with an overpotential of 1.81 V. We have also qualitatively compared results from current work with the work of Wannakao S. et al. and explained the effect of adsorbate coverage in determining the energetics of CO2 reduction reaction. Furthermore, we have compared the binding free energies of H* and CO* to explain that HER is possible along the CO2RR, but at higher applied negative potentials, our catalyst system suppresses HER thus promoting CO2 reduction to CH4. Future work needs to focus on two main tasks: (i) further studies on CO2 reduction to C2 products and (ii) improving the efficiency of current catalyst system by further minimizing the reduction potential. This can be achieved by introducing catalytically active dopants, thereby tuning the electronic structure properties or by varying the WC nanoparticle size and interaction with the graphene support. Both these tasks are currently the subject of consideration for future work to build on the results from this manuscript.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4344/9/7/604/s1.

Author Contributions

S.A was responsible for running majority of calculations, majority of data analysis, and majority of writing. Z.S initiated the first DFT calculations that led to the total scope of the research in this manuscript. R.B.R was responsible for creating graphics in S.I, reviewing and editing the manuscript before submission, and conceiving the research project and mentoring the students involved in the research in this manuscript.

Funding

This research received no external funding.

Acknowledgments

R.B.R. gratefully acknowledges financial and research support from the Department of Chemical Engineering at Villanova University through his startup package fund. Z.S. acknowledges support from Villanova University through the VURF summer program. S.A. acknowledges financial support from the Department of Chemical Engineering and the College of Engineering at Villanova University as well.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Yang, Y.; Ajmal, S.; Zheng, X.; Zhang, L. Efficient nanomaterials for harvesting clean fuels from electrochemical and photoelectrochemical CO2 reduction. Sustain. Energy Fuels 2018, 2, 510–537. [Google Scholar] [CrossRef]
  2. Tu, W.; Zhou, Y.; Zou, Z. Photocatalytic Conversion of CO2 into Renewable Hydrocarbon Fuels: State-of-the-Art Accomplishment, Challenges, and Prospects. Adv. Mater 2014, 26, 4607–4626. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, J.; Huang, Y.; Ye, W.; Li, Y. CO2 reduction: From the electrochemical to photochemical approach. Adv. Sci. 2017, 4, 1700194. [Google Scholar] [CrossRef] [PubMed]
  4. Kudo, A.; Miseki, Y. Heterogeneous photocatalyst materials for water splitting. Chem. Soc. Rev. 2009, 38, 253–278. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, Q.; Zhou, Y.; Tian, Z.; Chen, X.; Gao, J.; Zou, Z. Zn2 GeO4 crystal splitting toward sheaf-like, hyperbranched nanostructures and photocatalytic reduction of CO2 into CH 4 under visible light after nitridation. J. Mater. Chem. 2012, 22, 2033–2038. [Google Scholar] [CrossRef]
  6. Indrakanti, V.P.; Kubicki, J.D.; Schobert, H.H. Photoinduced activation of CO2 on Ti-based heterogeneous catalysts: Current state, chemical physics-based insights and outlook. Energy Environ. Sci. 2009, 2, 745–758. [Google Scholar] [CrossRef]
  7. Indrakanti, V.P.; Schobert, H.H.; Kubicki, J.D. Quantum mechanical modeling of CO2 interactions with irradiated stoichiometric and oxygen-deficient anatase TiO2 surfaces: Implications for the photocatalytic reduction of CO2. Energy Fuels 2009, 23, 5247–5256. [Google Scholar] [CrossRef]
  8. Yan, S.; Wan, L.; Li, Z.; Zou, Z. Facile temperature-controlled synthesis of hexagonal Zn2GeO4 nanorods with different aspect ratios toward improved photocatalytic activity for overall water splitting and photoreduction of CO2. Chem. Commun. 2011, 47, 5632–5634. [Google Scholar] [CrossRef]
  9. Chen, X.; Li, C.; Grätzel, M.; Kostecki, R.; Mao, S.S. Nanomaterials for renewable energy production and storage. Chem. Soc. Rev. 2012, 41, 7909–7937. [Google Scholar] [CrossRef]
  10. Anpo, M.; Yamashita, H.; Ichihashi, Y.; Fujii, Y.; Honda, M. Photocatalytic reduction of CO2 with H2O on titanium oxides anchored within micropores of zeolites: Effects of the structure of the active sites and the addition of Pt. J. Phys. Chem. B 1997, 101, 2632–2636. [Google Scholar] [CrossRef]
  11. Tong, H.; Ouyang, S.; Bi, Y.; Umezawa, N.; Oshikiri, M.; Ye, J. Nano-photocatalytic materials: Possibilities and challenges. Adv. Mater. 2012, 24, 229–251. [Google Scholar] [CrossRef] [PubMed]
  12. Hori, Y.; Kikuchi, K.; Suzuki, S. Production of CO and CH4 in electrochemical reduction of CO2 at metal electrodes in aqueous hydrogencarbonate solution. Chem. Lett. 1985, 14, 1695–1698. [Google Scholar] [CrossRef]
  13. Kuhl, K.P.; Cave, E.R.; Abram, D.N.; Jaramillo, T.F. New insights into the electrochemical reduction of carbon dioxide on metallic copper surfaces. Energy Environ. Sci. 2012, 5, 7050–7059. [Google Scholar] [CrossRef]
  14. Gonçalves, M.R.; Gomes, A.; Condeço, J.; Fernandes, R.; Pardal, T.; Sequeira, C.A.C.; Branco, J.B. Selective electrochemical conversion of CO2 to C2 hydrocarbons. Energy Convers. Manag. 2010, 51, 30–32. [Google Scholar] [CrossRef]
  15. Tackett, B.M.; Sheng, W.; Chen, J.G. Opportunities and Challenges in Utilizing Metal-Modified Transition Metal Carbides as Low-Cost Electrocatalysts. Joule 2017, 1, 253–263. [Google Scholar] [CrossRef] [Green Version]
  16. Levy, R.B.; Boudart, M. Platinum-like behavior of tungsten carbide in surface catalysis. Science 1973, 181, 547–549. [Google Scholar] [CrossRef] [PubMed]
  17. Zhong, Y.; Xia, X.; Shi, F.; Zhan, J.; Tu, J.; Fan, H.J. Transition metal carbides and nitrides in energy storage and conversion. Adv. Sci. 2016, 3, 1500286. [Google Scholar] [CrossRef]
  18. Liang, C.; Ying, P.; Li, C. Nanostructured β-Mo2C Prepared by Carbothermal Hydrogen Reduction on Ultrahigh Surface Area Carbon Material. Chem. Mater. 2002, 14, 3148–3151. [Google Scholar] [CrossRef]
  19. Michalsky, R.; Zhang, Y.J.; Medford, A.J.; Peterson, A.A. Departures from the Adsorption Energy Scaling Relations for Metal Carbide Catalysts. J. Phys. Chem. C 2014, 118, 13026–13034. [Google Scholar] [CrossRef]
  20. Wannakao, S.; Artrith, N.; Limtrakul, J.; Kolpak, A.M. Engineering Transition-Metal-Coated Tungsten Carbides for Efficient and Selective Electrochemical Reduction of CO2 to Methane. ChemSusChem 2015, 8, 2745–2751. [Google Scholar] [CrossRef]
  21. Cheng, J.; Hu, P.; Ellis, P.; French, S.; Kelly, G.; Lok, C.M. Density functional theory study of iron and cobalt carbides for Fischer− Tropsch synthesis. J. Phys. Chem. C 2009, 114, 1085–1093. [Google Scholar] [CrossRef]
  22. Stottlemyer, A.L.; Weigert, E.C.; Chen, J.G. Tungsten Carbides as Alternative Electrocatalysts: From Surface Science Studies to Fuel Cell Evaluation. Ind. Eng. Chem. Res. 2011, 50, 16–22. [Google Scholar] [CrossRef]
  23. Hunt, S.T.; Kokumai, T.M.; Zanchet, D.; Román-Leshkov, Y. Alloying tungsten carbide nanoparticles with tantalum: Impact on electrochemical oxidation resistance and hydrogen evolution activity. J. Phys. Chem. C 2015, 24, 13691–13699. [Google Scholar] [CrossRef]
  24. Hunt, S.T.; Milina, M.; Alba-Rubio, A.C.; Hendon, C.H.; Dumesic, J.A.; Román-Leshkov, Y. Self-assembly of noble metal monolayers on transition metal carbide nanoparticle catalysts. Science 2016, 352, 974–978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Obradović, M.D.; Gojković, S.L.; Elezović, N.R.; Ercius, P.; Radmilović, V.R.; Vračar, L.D.; Krstajić, N.V. The kinetics of the hydrogen oxidation reaction on WC/Pt catalyst with low content of Pt nano-particles. J. Electroanal. Chem. 2012, 671, 24–32. [Google Scholar] [CrossRef]
  26. Esposito, D.V.; Hunt, S.T.; Kimmel, Y.C.; Chen, J.G. A New Class of Electrocatalysts for Hydrogen Production from Water Electrolysis: Metal Monolayers Supported on Low-Cost Transition Metal Carbides. J. Am. Chem. Soc. 2012, 134, 3025–3033. [Google Scholar] [CrossRef] [PubMed]
  27. Wannakao, S.; Artrith, N.; Limtrakul, J.; Kolpak, A.M. Catalytic Activity and Product Selectivity Trends for Carbon Dioxide Electroreduction on Transition Metal-Coated Tungsten Carbides. J. Phys. Chem. C 2017, 121, 20306–20314. [Google Scholar] [CrossRef]
  28. Chhina, H.; Campbell, S.; Kesler, O. Thermal and electrochemical stability of tungsten carbide catalyst supports. J. Power Sources 2007, 164, 431–440. [Google Scholar] [CrossRef]
  29. Liu, Y.; Kelly, T.G.; Chen, J.G.; Mustain, W.E. Metal carbides as alternative electrocatalyst supports. ACS Catal. 2013, 3, 1184–1194. [Google Scholar] [CrossRef] [PubMed]
  30. Wang, L.; Li, M.; Huang, Z.; Li, Y.; Qi, S.; Yi, C.; Yang, B. Ni–WC/C nanocluster catalysts for urea electrooxidation. J. Power Sources 2014, 264, 282–289. [Google Scholar] [CrossRef]
  31. Antolini, E. Graphene as a new carbon support for low-temperature fuel cell catalysts. Appl. Catal. B Environ. 2012, 123, 52–68. [Google Scholar] [CrossRef]
  32. Ross, P.N., Jr.; Stonehart, P. The relation of surface structure to the electrocatalytic activity of tungsten carbide. J. Catal. 1977, 48, 42–59. [Google Scholar] [CrossRef]
  33. Ooka, H.; Figueiredo, M.C.; Koper, M.T. Competition between hydrogen evolution and carbon dioxide reduction on copper electrodes in mildly acidic media. Langmuir 2017, 33, 9307–9313. [Google Scholar] [CrossRef] [PubMed]
  34. MedeA, M. Exploration and Design Analysis; Materials Design Inc.: New York, NY, USA, 2016. [Google Scholar]
  35. Kresse, G.; Furthmüller, J. Software VASP, vienna. Rev. B 1996, 54, 169. [Google Scholar]
  36. Hafner, J. Ab-initio simulations of materials using VASP: Density-functional theory and beyond. J. Comput. Chem. 2008, 29, 2044–2078. [Google Scholar] [CrossRef] [PubMed]
  37. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  38. Piotrowski, M.J.; Ungureanu, C.G.; Tereshchuk, P.; Batista, K.E.; Chaves, A.S.; Guedes-Sobrinho, D.; Da Silva, J.L. Theoretical Study of the Structural, Energetic, and Electronic Properties of 55-Atom Metal Nanoclusters: A DFT Investigation within van der Waals Corrections, Spin–Orbit Coupling, and PBE+U of 42 Metal Systems. J. Phys. Chem. C 2016, 120, 28844–28856. [Google Scholar] [CrossRef]
  39. Bucko, T.; Hafner, J.; Lebegue, S.; Angyán, J.G. Improved Description of the Structure of Molecular and Layered Crystals: Ab Initio DFT Calculations with van der Waals Corrections. J. Phys. Chem. A 2010, 114, 11814–11824. [Google Scholar] [CrossRef]
  40. Klimeš, J.; Michaelides, A. Perspective: Advances and challenges in treating van der Waals dispersion forces in density functional theory. J. Chem. Phys. 2012, 137, 120901. [Google Scholar] [CrossRef] [Green Version]
  41. Watson, G.W.; Wells, R.P.K.; Willock, D.J.; Hutchings, G.J. π adsorption of ethene on to the {111} surface of copper: A periodic ab initio study of the effect of k-point sampling on the energy, atomic and electronic structure. Surf. Sci. 2000, 459, 93–103. [Google Scholar] [CrossRef]
  42. Peterson, A.A.; Abild-Pedersen, F.; Studt, F.; Rossmeisl, J.; Nørskov, J.K. How copper catalyzes the electroreduction of carbon dioxide into hydrocarbon fuels. Energy Environ. Sci. 2010, 3, 1311–1315. [Google Scholar] [CrossRef]
  43. Nørskov, J.K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.R.K.J.; Kitchin, J.R.; Bligaard, T.; Jonsson, H. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108, 17886–17892. [Google Scholar] [CrossRef]
  44. Yoo, J.S.; Christensen, R.; Vegge, T.; Nørskov, J.K.; Studt, F. Theoretical Insight into the Trends that Guide the Electrochemical Reduction of Carbon Dioxide to Formic Acid. ChemSusChem 2016, 9, 358–363. [Google Scholar] [CrossRef] [PubMed]
  45. Kuhl, K.P.; Hatsukade, T.; Cave, E.R.; Abram, D.N.; Kibsgaard, J.; Jaramillo, T.F. Electrocatalytic Conversion of Carbon Dioxide to Methane and Methanol on Transition Metal Surfaces. J. Am. Chem. Soc. 2014, 136, 14107–14113. [Google Scholar] [CrossRef] [PubMed]
  46. Hirunsit, P. Electroreduction of Carbon Dioxide to Methane on Copper, Copper–Silver, and Copper–Gold Catalysts: A DFT Study. J. Phys. Chem. C 2013, 117, 8262–8268. [Google Scholar] [CrossRef]
  47. Roux, B.; Simonson, T. Implicit solvent models. Biophys. Chem. 1999, 78, 1–20. [Google Scholar] [CrossRef]
  48. Cramer, C.J.; Truhlar, D.G. Implicit solvation models: Equilibria, structure, spectra, and dynamics. Chem. Rev. 1999, 99, 2161–2600. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Lowest energy pathway of CO2 reduction to HCOOH at U = 0 V on WC/graphene. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 1. Lowest energy pathway of CO2 reduction to HCOOH at U = 0 V on WC/graphene. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g001
Figure 2. Lowest energy pathway of CO2 reduction to CO at U = 0 V on WC/graphene. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 2. Lowest energy pathway of CO2 reduction to CO at U = 0 V on WC/graphene. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g002
Figure 3. Lowest energy pathway of CO2 reduction to CH3OH on the WC/graphene excluding solvation effect. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 3. Lowest energy pathway of CO2 reduction to CH3OH on the WC/graphene excluding solvation effect. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g003
Figure 4. Lowest energy pathways of CO2 reduction to CH3OH on the WC/graphene. Blue pathway: including solvation energy correction factor. Red pathway: without solvation energy correction factor. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 4. Lowest energy pathways of CO2 reduction to CH3OH on the WC/graphene. Blue pathway: including solvation energy correction factor. Red pathway: without solvation energy correction factor. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g004
Figure 5. Lowest energy pathway of CO2 reduction to CH4 on the WC/graphene. Red: excluding the solvation effect. Blue: including solvation effect. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 5. Lowest energy pathway of CO2 reduction to CH4 on the WC/graphene. Red: excluding the solvation effect. Blue: including solvation effect. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g005
Figure 6. Reaction energy profiles for CO2 reduction to CH3OH at 0 V vs. RHE (blue), thermodynamic limiting potential of (UL) −1.79 V vs. RHE (black) and equilibrium potential (UE) of −0.03 V vs. RHE (red). Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 6. Reaction energy profiles for CO2 reduction to CH3OH at 0 V vs. RHE (blue), thermodynamic limiting potential of (UL) −1.79 V vs. RHE (black) and equilibrium potential (UE) of −0.03 V vs. RHE (red). Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g006
Figure 7. Reaction energy profiles for CO2 reduction to CH4 at 0V vs. RHE (blue), thermodynamic limiting potential (UL) of −0.84 V vs. RHE (black) and equilibrium potential (UE) of 0.17 V vs. RHE (red). Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 7. Reaction energy profiles for CO2 reduction to CH4 at 0V vs. RHE (blue), thermodynamic limiting potential (UL) of −0.84 V vs. RHE (black) and equilibrium potential (UE) of 0.17 V vs. RHE (red). Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g007
Figure 8. Role of surface coverage on reaction energetics. Blue: one adsorbed species per WC/graphene catalyst system. Red: Co-adsorbed species per WC/graphene catalyst system. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 8. Role of surface coverage on reaction energetics. Blue: one adsorbed species per WC/graphene catalyst system. Red: Co-adsorbed species per WC/graphene catalyst system. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g008
Figure 9. Comparison of reaction energetics of CO2 reduction on WC nanocluster (blue) and graphene supported WC nanocluster (red). Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 9. Comparison of reaction energetics of CO2 reduction on WC nanocluster (blue) and graphene supported WC nanocluster (red). Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g009
Figure 10. Lowest energy pathway for CO2 reduction to CH4. Red: on graphene supported WC. Blue: on WC (0001) from the work of Wannakao S. et al. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Figure 10. Lowest energy pathway for CO2 reduction to CH4. Red: on graphene supported WC. Blue: on WC (0001) from the work of Wannakao S. et al. Color code of atoms: red: oxygen, black: carbon, white: hydrogen.
Catalysts 09 00604 g010
Figure 11. Reaction energy profiles for hydrogen evolution reaction (HER) at 0 V vs. RHE (blue), thermodynamic limiting potential of (UL) −0.98 V vs. RHE (black).
Figure 11. Reaction energy profiles for hydrogen evolution reaction (HER) at 0 V vs. RHE (blue), thermodynamic limiting potential of (UL) −0.98 V vs. RHE (black).
Catalysts 09 00604 g011
Figure 12. Comparison of adsorption/binding free energies of formation of H* and CO* on WC/Graphene catalyst system at various applied potentials ranging from 0 to 1 V vs. RHE. Red: H* binding free energies Black: CO* binding free energies.
Figure 12. Comparison of adsorption/binding free energies of formation of H* and CO* on WC/Graphene catalyst system at various applied potentials ranging from 0 to 1 V vs. RHE. Red: H* binding free energies Black: CO* binding free energies.
Catalysts 09 00604 g012
Figure 13. Possible reaction energy mechanisms for electrochemical reduction of CO2 to CO, CH4, and CH3OH. OH, and H2O formed along the reaction pathways are omitted in the figure for clarity.
Figure 13. Possible reaction energy mechanisms for electrochemical reduction of CO2 to CO, CH4, and CH3OH. OH, and H2O formed along the reaction pathways are omitted in the figure for clarity.
Catalysts 09 00604 g013
Table 1. Overall reactions for CO2 reduction to different products and their equilibrium potentials (U, V vs. RHE) [45].
Table 1. Overall reactions for CO2 reduction to different products and their equilibrium potentials (U, V vs. RHE) [45].
ReactionU (V vs. RHE)
2(H+ + e) → H20
CO2 + 2(H+ + e−) → CO + H2O−0.10
CO2 + 2(H+ + e−) → HCOOH−0.20
CO2 + 6(H+ + e−) → CH3OH + H2O−0.03
CO2 + 8(H+ + e−) → CH4 + H2O0.17

Share and Cite

MDPI and ACS Style

Ananthaneni, S.; Smith, Z.; Rankin, R.B. Graphene Supported Tungsten Carbide as Catalyst for Electrochemical Reduction of CO2. Catalysts 2019, 9, 604. https://doi.org/10.3390/catal9070604

AMA Style

Ananthaneni S, Smith Z, Rankin RB. Graphene Supported Tungsten Carbide as Catalyst for Electrochemical Reduction of CO2. Catalysts. 2019; 9(7):604. https://doi.org/10.3390/catal9070604

Chicago/Turabian Style

Ananthaneni, Sahithi, Zachery Smith, and Rees B. Rankin. 2019. "Graphene Supported Tungsten Carbide as Catalyst for Electrochemical Reduction of CO2" Catalysts 9, no. 7: 604. https://doi.org/10.3390/catal9070604

APA Style

Ananthaneni, S., Smith, Z., & Rankin, R. B. (2019). Graphene Supported Tungsten Carbide as Catalyst for Electrochemical Reduction of CO2. Catalysts, 9(7), 604. https://doi.org/10.3390/catal9070604

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop