Computational Discovery of Nickel-Based Catalysts for CO2

Sep 7, 2017 - Electrochemical reduction of CO2 into chemical fuels is crucial to clean energy production and environment remediation. First-principles...
0 downloads 11 Views 2MB Size
Subscriber access provided by GRIFFITH UNIVERSITY

Article

Computational Discovery of Ni-Based Catalysts for CO Reduction to Formic Acid 2

Zhonglong Zhao, Zhengzheng Chen, and Gang Lu J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b06895 • Publication Date (Web): 07 Sep 2017 Downloaded from http://pubs.acs.org on September 9, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Computational Discovery of Ni-based Catalysts for CO2 Reduction to Formic Acid Zhonglong Zhao, Zhengzheng Chen, and Gang Lu* Department of Physics and Astronomy, California State University Northridge, California, 91330, United States ABSTRACT: Electrochemical reduction of CO2 into chemical fuels is crucial to clean energy production and environment remediation. First-principles calculations are performed to elucidate reaction mechanism of CO2 reduction to formic acid on Ni-based catalysts. The origin of CO poisoning is examined and a novel design strategy is proposed to eliminate CO poisoning. Three design criteria are derived based on which computational screening is performed to identify several Ni-based near-surface-alloys (NSAs) with both high selectivity and reactivity. The effect of elastic strain on CO2 reduction is studied on these NSAs. We predict that Ni/Ti, Cu/Ni, and strained Cu/Ni NSAs could lead to highly selective and efficient production of formic acid.

1 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 25

INTRODUCTION Electrochemical reduction of CO2 into chemical fuels, such as hydrocarbons, using renewable energy sources (solar, hydro, and wind) represents a promising route to mitigate our dependence on fossil fuels and to ease the threat of global warming.1-4 However, despite intense research in recent years, the discovery of efficient, selective and stable catalysts for CO2 reduction (CO2RR) still eludes us. For example, although Cu is regarded as the best monometallic catalyst, capable of producing significant quantities of hydrocarbons, it does so inefficiently with a high overpotential (~1 V) and a poor selectivity (16 major and minor products).5-7 Formic acid (HCOOH) is one of the main products of CO2RR, and with only two proton-electron pairs in the reaction, it is also one of the simplest products. Despite its simplicity, formic acid has been found as an important fuel in fuel cells,8 a promising feedstock for fine chemicals,9 a preservative and antibacterial agent,10 and a suitable material for hydrogen storage.11 Hence it is of broad scientific and technological importance to develop catalysts that can efficiently and selectively reduce CO2 to useful products such as formic acid at low overpotentials.12-22 There are two competing reaction pathways for CO2RR to formic acid (HCOOH):13, 18 +



+



H +e H +e CO 2  → COOH*  → HCOOH

+



+

(1)



H +e H +e CO 2  → HCOO*  → HCOOH

(2)

In the first pathway, carboxyl (COOH*, * indicates adsorbed species) is the intermediate and in the second pathway, formate (HCOO*) is the intermediate. Among the factors that adversely affect the activity and the selectivity of the catalysts, CO poisoning stands out, particularly on transition metals, such as Fe, Co, Ni, Pd, and Pt.6,

14, 23-24

The poisoning results from the

following reaction (pathway 3) that shares the same intermediate of carboxyl as in the pathway 1:13 2 / 21 ACS Paragon Plus Environment

Page 3 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

+



+



H +e H +e CO2  → COOH*  → CO* + H 2 O  → CO

(3)

Due to its strong binding to the metal surface, CO*, once formed, can block active sites and hinder the reaction to HCOOH. To remove CO* from the surface, a high voltage would have to be applied, which may but not always produce hydrogen as the main product. The CO poisoning has been observed by Hori et al.23 and Kuhl et al.6 where the Faradaic yields of H2 on Ni, Pt, and Fe electrodes can reach 90% and only 2% of hydrocarbons is produced at an applied potential of ~1 V vs. RHE. To alleviate CO poisoning on the transition metals, research effort has been devoted to weakening CO affinity to the surface.14, 17 For example, Min et al.14 demonstrated that with a short oxidative treatment to the poisoned Pd surface, formate production can be reactivated. However, such treatment has to be performed repeatedly or the poisoning will return. Guided by theoretical predictions, Klinkova et al.17 synthesized Pd nanoparticles with high-index facets which could weaken the CO binding affinity and in turn enhance the catalytic activity. However, CO poisoning persisted and periodic surface treatments were required for a long-term performance of the catalyst. In this work, we propose an alternative strategy which deactivates the first pathway involving carboxyl. If carboxyl is no longer produced, the competing CO production (pathway 3) is no longer possible, thus CO poisoning is completely eliminated. By means of first-principles calculations, we show that this strategy can be realized in practice by modulating the free energies of the key intermediates (COOH* and HCOO*). Following a computational screening, we predict several inexpensive Ni-based near-surface alloys (NSAs) to be more active and selective than the best performing metal catalysts for the production of formic acid. We should point out that several Ni-based catalysts have been proposed as catalysts for CO2 reduction to other products than formic acid, for example, Ni-Ga alloys for the

3 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

production of methane, ethylene, and ethane,25 Ni-In alloys for CO production,26 atomic Ni on defective graphene support for methanol production,27 and atomic Ni on metal-organic frameworks support for CO production.28 COMPUTATIONAL DETAILS The Ni (111) surface and related NSAs were simulated with four atomic layers and a 3×3 inplane supercell. A thicker six-layer slab was also used and found to yield converged results as the four-layer model (Figure S1). The adjacent slabs were separated by a 17 Å vacuum in the normal direction of the surface. The top three atomic layers were allowed to relax while the bottom layer was fixed. The equilibrium lattice parameter of Ni was calculated as 3.526 Å. The spin-polarized DFT calculations were performed using Vienna ab initio simulation package (VASP).29 The ion-electron interaction was described by the projector-augmented wave method.30 The revised Perdew-Burke-Ernzerhof (RPBE) functional was used for the exchangecorrelation potential.31-32 RPBE+U (U = 1 eV) method33 was employed to minimize the selfinteraction error and yielded the magnetic moment of Ni as 0.66 µB, same as the experimental value.34-35 A plane-wave energy cutoff of 400 eV was used for all calculations and the Brillouinzone was sampled with a 3 × 3 × 1 k-mesh according to the Monkhorst-Pack scheme.36 Ab initio MD simulations were performed using NVT (Number of particles, Volume, Temperature) ensemble37-39 and the Brillouin zone integration was restricted to the Γ point. The computational hydrogen electrode (CHE) model40 was used to calculate the free energy diagrams. In the CHE model, the chemical potential of a proton-electron pair is defined in equilibrium with half of that of gaseous H2 at U = 0 V, 101325 Pa, any pH values and temperatures. The chemical potential is shifted by –eU (e is the elementary positive charge)

4 / 21 ACS Paragon Plus Environment

Page 4 of 25

Page 5 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

when an external potential U is applied to the system. Details about the free energy calculations can be found in the Supporting Information. RESULTS AND DISCUSSION In this work, we choose Ni as a representative of transition metal catalysts to illustrate our design strategy for CO2RR to formic acid. Figure 1a depicts the reaction pathways for HCOOH production on Ni (111) surface. The reaction pathway consisted of COOH* (HCOO*) as the intermediate is shown in black (blue) color, and labeled as route 1 (route 2) in the figure. In addition, the undesirable pathway to CO poisoning (route 3) and two branched routes to COHOH* and H2COO* are also included in Figure 1a. Although the intermediates COHOH* and H2COO* can be further reduced to hydrocarbons,13, 18, 41 the corresponding overpotentials are much higher, thus not examined here. The further reduction of CO* is not considered here because the Faradaic yields of its reduced products on Ni are less than 2% as reported in experiment.6

Figure 1. (a) Reaction pathways for CO2 reduction to formic acid (HCOOH). (b) The free energy diagram for CO2 reduction to HCOOH on Ni (111) surface at 0 V vs. RHE. Relevant 5 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

free energy differences (∆G1, ∆G2, and ∆G3) are labeled. The same color scheme is used for (a) and (b). Figure 1b displays the free energy diagram for HCOOH production on Ni (111) surface at 0 V vs. RHE from first-principles calculations. The limiting potential required to activate each reaction step can be derived as UL = –∆G/e according to the Computational Hydrogen Electrode (CHE) model40 where ∆G represents the free energy difference between the two adjacent reaction steps. The calculation of the free energy can be found in Supporting Information. The most negative UL along each pathway is defined as its overpotential. Clearly, on Ni (111)surface, the route 1 (black) is activated while the route 2 (blue) is deactivated because the overpotential of the route 1 (–∆G1/e = –0.16 V) is less negative than that (–∆G3/e = –0.33 V) of the route 2. However, for the route 1, the intermediate COOH* is reduced to CO* and H2O instead of HCOOH, owing to a large free energy drop of the former. Since CO* is strongly bound to the Ni (111) surface, its desorption would require a very high energy input (Figure 1b). Hence, the Ni surface is poisoned by CO* and no formic acid can be produced. These results are consistent with the experimental finding that no significant amount of hydrocarbons can be detected on pure Ni electrode even at an elevated potential above 1 V.6, 23 Since the route 1 leads to CO poisoning, one turns to the route 2 for a possible solution. Evidently, the route 2 has one crucial advantage: there are no viable byproducts in this route because the only competing reaction HCOO* → H2COO* would require a very high limiting potential of 1.5 V. Thus, if the route 1 is deactivated and the route 2 is activated, CO poisoning would be eliminated and formic acid would be the only product. More specifically, we require that the first step in the route 1, i.e., the formation of COOH* be turned off, while both steps in the route 2 be turned on. To satisfy both requirements, we must have 6 / 21 ACS Paragon Plus Environment

Page 6 of 25

Page 7 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

∆G1 − ∆G 2 > 0

(4)

∆G1 − ∆G 3 > 0

(5)

Catalysts satisfying both conditions (selectivity criteria) are predicted to produce formic acid without CO poisoning. Now that the route 2 is activated, the next step is to minimize its overpotential. This can be achieved by meeting the following condition (activity criterion): ∆G 2 = ∆G 3

(6)

This equation is derived by noting that the overpotential is the larger of ∆G2 and ∆G3. The three criteria provide a guidance to computationally screen Ni-based catalysts for HCOOH production with a low overpotential and a high selectivity.

Figure 2. Free energy change between the reaction steps (∆G1: CO2 → COOH*, ∆G2: CO2 → HCOO*, ∆G3: HCOO* → HCOOH) on (a) Ni-based NSAs, Cu (211) and Ni (111)

7 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

surfaces; (c) strained Cu*/Ni as a function of the free energy of HCOO*. The NSAs above zero meet the two selectivity criteria as indicated by solid symbols. The corresponding volcano plots via the route 2 are shown in (b) and (d). Schematics for the overlayer (M*/Ni) and the subsurface (Ni/M) NSAs are shown in the insets of (b). Blue and gray spheres denote the M atoms and Ni, respectively. In the following, we focus on Ni-based NSAs that meet the selectivity and the activity criteria. Twenty-seven alloy elements M (M = Sc, Ti, V, Cr, Mn, Fe, Co, Cu, Zn, Y, Zr, Nb, Mo, Tc, Ru, Rh, Pd, Ag, Cd, Hf, Ta, W, Re, Os, Ir, Pt, Au) have been considered in this work. A monolayer of M is introduced either to the surface layer (M*/Ni) or the subsurface layer (Ni/M) on Ni (111) surface. We note that bimetallic NSAs composed of diverse metal elements have been previously synthesized. For example, electrodeposition of Cu on Pt (111) could yield both Cu*/Pt and Pt/Cu NSAs depending on annealing temperature.42 Similarly, electrodeposition of Au on polycrystalline W produced Au*/W NSAs with controllable Au layers.43 Hence our design strategy can be realized in practice. The relative stability of M*/Ni and Ni/M NSAs can be evaluated based on segregation energy (Eseg), defined as the energy difference between NSAs with an M atom substituting for a Ni atom in the subsurface and the surface.44 If Eseg is negative, the M atom would prefer to segregate to the surface (M*/Ni); conversely, it would prefer to segregate to the subsurface (Ni/M). Based on Eseg, we can narrow down Ni-based NSAs which are energetically stable (Table S1). Furthermore, it has been suggested that the presence of strongly adsorbed species could destabilize the NSA surface.42 Hence, we perform additional screening by fully relaxing the atomic geometries of the NSAs with HCOO* and COOH* on the surfaces. We have thus identified six stable Ni-based NSAs (Ni/Ti, Ni/V, Ni/Cr, Ni/Co, Cu*/Ni, and Ni/Mo) which show 8 / 21 ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

no substantial surface reconstruction. Lastly, we carry out ab initio molecular dynamics (MD) simulations at 300 K to examine dynamical stability of these six NSAs with HCOO* on the surfaces. We find that the NSAs remain stable during the MD simulations, as evidenced by negligible changes in the radial distribution functions (Figure S2). After establishing the thermodynamic stability of the NSAs, we proceed to examine their selectivity and activity. As shown in Figure 2a, we find that Ni/Cr, Ni/Mo, Ni/Ti, and Cu*/Ni meet the selectivity criteria as their corresponding data points are above zero. The remaining alloys satisfy only one of the two criteria. The overpotentials of the NSAs are displayed in Figure 2b. Among the four NSAs meeting the selectivity criteria, Cu*/Ni and Ni/Ti have the lowest overpotentials, thus are predicted as the most promising catalysts for the production of formic acid. In particular, Cu*/Ni is the most selective catalyst with the maximum free energy differences, while Ni/Ti is the most active catalyst with the lowest overpotential. It is well known that elastic strains can influence adsorbate binding energy on a catalyst surface, thus change its catalytic activity.45-47 Recently, there is a surge of interest to optimize the performance of metal catalysts by tuning the strain on the catalyst surface.48-57 The elastic strain can be applied to the catalysts directly via mechanical loading48-49 or indirectly via lattice mismatch in core/shell nanostructures (e.g., nanoparticles and nanoplates).51-57 Although CO2RR on the core/shell nanostructures is beyond the scope of this work, the analogy between NSAs and multimetallic core/shell nanostructures prompts us to examine how strain may modulate the selectivity and the reactivity of the NSAs. To this end, biaxial tensile and compressive strains, εXY, are applied to Cu*/Ni surface in and directions. As shown in Figure 2d, we find that +3% biaxial tensile strain leads to the lowest overpotential, while the selectivity drops slightly (Figure 2c). Hence, the biaxial tensile strain is beneficial to Cu*/Ni catalyst. Since it is

9 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 25

generally believed that the catalysts on the opposite legs of the volcano plot tend to have opposite responses to strain, we speculate that Ni/Mo may become more active under appropriate compressive strains.

Figure 3. Free energy diagrams for CO2 electroreduction to HCOOH on (a) Ni/Ti, (b) Cu*/Ni, and (c) strained (εXY = +3%) Cu*/Ni NSAs at 0 V vs. RHE. In Figure 3, we summarize the free energy diagrams for the production of formic acid on strain-free Ni/Ti and Cu*/Ni surfaces, as well as strained (εXY = +3%) Cu*/Ni surface. Clearly, the route 2 is now strongly preferred and the route 1 is effectively deactivated. As a result, CO poisoning is eliminated on the three surfaces. It is evident from Figure 3 that the deactivation of the route 1 on the proposed NSAs is due to an increase of COOH* free energy, which corresponds to a weakened adsorption on the surfaces. The weaker adsorption can be understood from a downshift of the d-band center in the local density of states (dxz, dyz, and dz2 components) 10 / 21 ACS Paragon Plus Environment

Page 11 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of the surface atoms in the NSAs as compared to that in pure Ni (Figure S3). Note that this correlation between the adsorption strength and the d-band center has also been observed in other metal catalysts.50,

58

The overpotentials for the route 2 on Ni/Ti and strained Cu*/Ni are

calculated as –0.15 and –0.14 V, respectively. These NSAs are predicted to be more active than some of the most efficient catalysts for HCOOH production, including partially oxidized and atomically thin Co layers (–0.24 V),16 Fe carbonyl cluster (–0.23 ~ –0.44 V),15 Pd nanoparticles (NPs) with high-index facets (–0.20V),17 and Pd NPs dispersed on a carbon support (–0.20 V).14

Figure 4. Free energy diagram for HER on various NSA surfaces at 0 V vs. RHE. The adsorption geometries are shown in the insets. The pink and grey spheres represent H and metals, respectively. Finally, we examine hydrogen evolution reaction (HER) on these NSAs. A superb catalyst for CO2RR should possess a high overpotential for HER to minimize the consumption of protons and electrons. Hence, the catalyst should have either a very strong or very weak binding to H*. Between the two, the weak binding of H* is more desirable so that H* will not poison the surface.59 As shown in Figure 4, we find that the predicted NSAs exhibit weak binding to H* (or 11 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 25

large positive free energy), leading to suppressed HER. More specifically, Ni/Ti and Cu*/Ni are predicted to have the most negative overpotentials for HER (–0.08 V and –0.12 V). In comparison, the HER overpotential on Cu (211) surfaces is only –0.03 V. However, we note that the overpotentials for HER on the proposed NSAs are slightly (~0.08 V) less negative than those for HCOOH production, which means that HER is a still major competing reaction. The weak H* binding on the proposed NSA surfaces can also be understood from the downshift of the dband center as mentioned above. A previous work by Liu et al.60 showed that the overpotential for HER on small metal clusters could be increased by employing a defective graphene support which modified the electronic structure of the active sites. CONCLUSIONS In summary, we have performed first-principles calculations to elucidate the reaction mechanism of CO2RR to formic acid on Ni-based catalysts. We revealed the origin of CO poisoning on the metal surface and proposed a strategy to eliminate CO poisoning by activating the HCOO* pathway and deactivating the competing COOH* pathway. Three design criteria were derived based on which computational screening was performed to identify the Ni-based NSAs with high selectivity and reactivity. The effect of elastic strain on CO2RR was examined on these NSAs. We predicted that Ni/Ti, Cu*/Ni, and strained Cu*/Ni surfaces could lead to highly selective and efficient production of formic acid. It is hoped that this work could inspire future discovery of superior catalysts for CO2RR. ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website.

12 / 21 ACS Paragon Plus Environment

Page 13 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Detailed information regarding the free energy calculations, NSA segregation energies, and supportive figures. AUTHORINFORMATION Corresponding Author *E-mail: [email protected] Notes The authors declare no competing financial interest. ACKNOWLEDGMENTS This work was supported by the Office of Naval Research (N00014-15-1-2092) and National Science Foundation (DMR-1205734). REFERENCES (1)

Lewis, N. S.; Nocera, D. G. Powering the planet: chemical challenges in solar energy

utilization. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 15729-15735. (2)

Olah, G. A.; Prakash, G. K. S.; Goeppert, A. Anthropogenic chemical carbon cycle for a

sustainable future. J. Am. Chem. Soc. 2011, 133, 12881-12898. (3)

Jhong, H.-R. M.; Ma, S.; Kenis, P. J. A. Electrochemical conversion of CO2 to useful

chemicals: current status, remaining challenges, and future opportunities. Curr. Opin. Chem. Eng. 2013, 2, 191-199. (4)

Wilcox, J. Carbon capture; Springer: New York, 2012.

(5)

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.

13 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(6)

Page 14 of 25

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. (7)

Peterson, A. A.; Nørskov, J. K. Activity descriptors for CO2 electroreduction to methane

on transition-metal catalysts. J. Phys. Chem. Lett. 2012, 3, 251-258. (8)

Yu, X.; Pickup, P. G. Recent advances in direct formic acid fuel cells (DFAFC). J. Power

Sources 2008, 182, 124-132. (9)

Mura, M. G.; Luca, L. D.; Giacomelli, G.; Porcheddu, A. Formic acid: a promising bio-

renewable feedstock for fine chemicals. Adv. Synth. Catal. 2012, 354, 3180-3186. (10)

Hietala, J.; Vuori, A.; Johnsson, P.; Pollari, I.; Reutemann, W.; Kieczka, H. Formic acid.

Ullmann's encyclopedia of industrial chemistry; Wiley: New York, 2016; pp 1-22. (11)

Grasemann, M.; Laurenczy, G. Formic acid as a hydrogen source - recent developments

and future trends. Energy Environ. Sci. 2012, 5, 8171-8181. (12)

Sen, S.; Liu, D.; Palmore, G. T. R. Electrochemical reduction of CO2 at copper

nanofoams. ACS Catal. 2014, 4, 3091-3095. (13)

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. (14)

Min, X.; Kanan, M. W. Pd-catalyzed electrohydrogenation of carbon dioxide to formate:

high mass activity at low overpotential and identification of the deactivation pathway. J. Am. Chem. Soc. 2015, 137, 4701-4708.

14 / 21 ACS Paragon Plus Environment

Page 15 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(15)

Taheri, A.; Thompson, E. J.; Fettinger, J. C.; Berben, L. A. An iron electrocatalyst for

selective reduction of CO2 to formate in water: including thermochemical insights. ACS Catal. 2015, 5, 7140-7151. (16)

Gao, S.; Lin, Y.; Jiao, X.; Sun, Y.; Luo, Q.; Zhang, W.; Li, D.; Yang, J.; Xie, Y. Partially

oxidized atomic cobalt layers for carbon dioxide electroreduction to liquid fuel. Nature 2016, 529, 68-71. (17)

Klinkova, A.; Luna, P. D.; Dinh, C.-T.; Voznyy, O.; Larin, E. M.; Kumacheva, E.;

Sargent, E. H. Rational design of efficient palladium catalysts for electroreduction of carbon dioxide to formate. ACS Catal. 2016, 6, 8115-8120. (18)

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. (19)

Zhang, S.; Kang, P.; Meyer, T. J. Nanostructured tin catalysts for selective

electrochemical reduction of carbon dioxide to formate. J. Am. Chem. Soc. 2014, 136, 1734-1737. (20)

Studt, F.; Sharafutdinov, I.; Abild-Pedersen, F.; Elkjær, C. F.; Hummelshøj, J. S.; Dahl,

S.; Chorkendorff, I.; Nørskov, J. K. Discovery of a Ni-Ga catalyst for carbon dioxide reduction to methanol. Nat. Chem. 2014, 6, 320-324. (21)

Li, C. W.; Ciston, J.; Kanan, M. W. Electroreduction of carbon monoxide to liquid fuel

on oxide-derived nanocrystalline copper. Nature 2014, 508, 504-507. (22)

Chen, Z.; Zhang, X.; Lu, G. Overpotential for CO2 electroreduction lowered on strained

penta-twinned Cu nanowires. Chem. Sci. 2015, 6, 6829-6835.

15 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(23)

Page 16 of 25

Hori, Y.; Wakebe, H.; Tsukamoto, T.; Koga, O. Electrocatalytic process of CO selectivity

in electrochemical reduction of CO2 at metal electrodes in aqueous media. Electrochim. Acta 1994, 39, 1833-1839. (24)

Akhade, S. A.; Luo, W.; Nie, X.; Bernstein, N. J.; Asthagiri, A.; Janik, M. J. Poisoning

effect of adsorbed CO during CO2 electroreduction on late transition metals. Phys. Chem. Chem. Phys. 2014, 16, 20429-20435. (25)

Torelli, D. A.; Francis, S. A.; Crompton, J. C.; Javier, A.; Thompson, J. R.; Brunschwig,

B. S.; Soriaga, M. P.; Lewis, N. S. Nickel–gallium-catalyzed electrochemical reduction of CO2 to highly reduced products at low overpotentials. ACS Catal. 2016, 6, 2100-2104. (26)

He, J.; Dettelbach, K. E.; Salvatore, D. A.; Li, T.; Berlinguette, C. P. High-throughput

synthesis of mixed-metal electrocatalysts for CO2 reduction. Angew. Chem., Int. Ed. 2017, 56, 6068-6072. (27)

Back, S.; Lim, J.; Kim, N.-Y.; Kim, Y.-H.; Jung, Y. Single-atom catalysts for CO2

electroreduction with significant activity and selectivity improvements. Chem. Sci. 2017, 8, 1090-1096. (28)

Zhao, C.; Dai, X.; Yao, T.; Chen, W.; Wang, X.; Wang, J.; Yang, J.; Wei, S.; Wu, Y.; Li,

Y. Ionic exchange of metal–organic frameworks to access single nickel sites for efficient electroreduction of CO2. J. Am. Chem. Soc. 2017, 139, 8078-8081. (29)

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. (30)

Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953-17979.

16 / 21 ACS Paragon Plus Environment

Page 17 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(31)

Hammer, B.; Hansen, L. B.; Nørskov, J. K. Improved adsorption energetics within

density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 1999, 59, 7413-7421. (32)

Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple.

Phys. Rev. Lett. 1996, 77, 3865-3868. (33)

Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P. Electron-

energy-loss spectra and the structural stability of nickel oxide: an LSDA+U study. Phys. Rev. B 1998, 57, 1505-1509. (34)

Chen, C. T.; Smith, N. V.; Sette, F. Exchange, spin-orbit, and correlation effects in the

soft-x-ray magnetic-circular-dichroism spectrum of nickel. Phys. Rev. B 1991, 43, 6785-6787. (35)

Vogel, J.; Sacchi, M. Polarization and angular dependence of the L2,3 absorption edges in

Ni(110). Phys. Rev. B 1994, 49, 3230-3234. (36)

Monkhorst, H. J.; Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B

1976, 13, 5188-5192. (37)

Nosé, S. A unified formulation of the constant temperature molecular dynamics methods.

J. Chem. Phys. 1984, 81, 511-519. (38)

Shuichi, N. Constant temperature molecular dynamics methods. Prog. Theor. Phys. 1991,

103, 1-46. (39)

Bylander, D. M.; Kleinman, L. Energy fluctuations induced by the Nosé thermostat. Phys.

Rev. B 1992, 46, 13756-13761. (40)

Nørskov, J. K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J. R.; Bligaard, T.;

Jónsson, H. Origin of the overpotential for oxygen reduction at a fuel-cell cathode. J. Phys. Chem. B 2004, 108, 17886-17892.

17 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(41)

Page 18 of 25

Luo, W.; Nie, X.; Janik, M. J.; Asthagiri, A. Facet dependence of CO2 reduction paths on

Cu electrodes. ACS Catal. 2016, 6, 219-229. (42)

Bandarenka, A. S.; Varela, A. S.; Karamad, M.; Calle-Vallejo, F.; Bech, L.; Perez-Alonso,

F. J.; Rossmeisl, J.; Stephens, I. E. L.; Chorkendorff, I. Design of an active site towards optimal electrocatalysis: overlayers, surface alloys and near-surface alloys of Cu/Pt(111). Angew. Chem., Int. Ed. 2012, 51, 11845-11848. (43)

Javier, A.; Baricuatro, J. H.; Kim, Y. G.; Soriaga, M. P. Overlayer Au-on-W near-surface

alloy for the selective electrochemical reduction of CO2 to methanol: empirical (DEMS) corroboration of a computational (DFT) prediction. Electrocatalysis 2015, 6, 493-497. (44)

Greeley, J.; Mavrikakis, M. Alloy catalysts designed from first principles. Nat. Mater.

2004, 3, 810-815. (45)

Mavrikakis, M.; Hammer, B.; Nørskov, J. K. Effect of strain on the reactivity of metal

surfaces. Phys. Rev. Lett. 1998, 81, 2819-2822. (46)

Gsell, M.; Jakob, P.; Menzel, D. Effect of substrate strain on adsorption. Science 1998,

280, 717-720. (47)

Kitchin, J. R.; Nørskov, J. K.; Barteau, M. A.; Chen, J. G. Role of strain and ligand

effects in the modification of the electronic and chemical properties of bimetallic surfaces. Phys. Rev. Lett. 2004, 93, 156801. (48)

Yan, K.; Kim, S. K.; Khorshidi, A.; Guduru, P. R.; Peterson, A. A. High elastic strain

directly tunes the hydrogen evolution reaction on tungsten carbide. J. Phys. Chem. C 2017, 121, 6177-6183.

18 / 21 ACS Paragon Plus Environment

Page 19 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(49)

Yan, K.; Maark, T. A.; Khorshidi, A.; Sethuraman, V. A.; Peterson, A. A.; Guduru, P. R.

The influence of elastic strain on catalytic activity in the hydrogen evolution reaction. Angew. Chem., Int. Ed. 2016, 55, 6175-6181. (50)

Kibler, L. A.; El-Aziz, A. M.; Hoyer, R.; Kolb, D. M. Tuning reaction rates by lateral

strain in a palladium monolayer. Angew. Chem., Int. Ed. 2005, 44, 2080-2084. (51)

Strasser, P.; Koh, S.; Anniyev, T.; Greeley, J.; More, K.; Yu, C.; Liu, Z.; Kaya, S.;

Nordlund, D.; Ogasawara, H.; et al. Lattice-strain control of the activity in dealloyed core–shell fuel cell catalysts. Nat. Chem. 2010, 2, 454-460. (52)

Greeley, J.; Nørskov, J. K.; Kibler, L. A.; El-Aziz, A. M.; Kolb, D. M. Hydrogen

evolution over bimetallic systems: understanding the trends. ChemPhysChem 2006, 7, 10321035. (53)

Tang, W.; Henkelman, G. Charge redistribution in core-shell nanoparticles to promote

oxygen reduction. J. Chem. Phys. 2009, 130, 194504. (54)

Bu, L.; Zhang, N.; Guo, S.; Zhang, X.; Li, J.; Yao, J.; Wu, T.; Lu, G.; Ma, J.-Y.; Su, D.;

et al. Biaxially strained PtPb/Pt core/shell nanoplate boosts oxygen reduction catalysis. Science 2016, 354, 1410-1414. (55)

Zhang, X.; Lu, G. Computational design of core/shell nanoparticles for oxygen reduction

reactions. J. Phys. Chem. Lett. 2014, 5, 292-297. (56)

Guo, S.; Zhang, X.; Zhu, W.; He, K.; Su, D.; Mendoza-Garcia, A.; Ho, S. F.; Lu, G.; Sun,

S. Nanocatalyst superior to Pt for oxygen reduction reactions: the case of core/shell Ag(Au)/CuPd nanoparticles. J. Am. Chem. Soc. 2014, 136, 15026-15033.

19 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(57)

Page 20 of 25

Zhang, S.; Zhang, X.; Jiang, G.; Zhu, H.; Guo, S.; Su, D.; Lu, G.; Sun, S. Tuning

nanoparticle structure and surface strain for catalysis optimization. J. Am. Chem. Soc. 2014, 136, 7734-7739. (58)

Rochana, P.; Wilcox, J. A theoretical study of CO adsorption on FeCo(100) and the effect

of alloying. Surf. Sci. 2011, 605, 681-688. (59)

Mistry, H.; Reske, R.; Zeng, Z.; Zhao, Z.-J.; Greeley, J.; Strasser, P.; Cuenya, B. R.

Exceptional size-dependent activity enhancement in the electroreduction of CO2 over Au nanoparticles. J. Am. Chem. Soc. 2014, 136, 16473-16476. (60)

Liu, C.; He, H.; Zapol, P.; Curtiss, L. A. Computational studies of electrochemical CO2

reduction on subnanometer transition metal clusters. Phys. Chem. Chem. Phys. 2014, 16, 2658426599.

20 / 21 ACS Paragon Plus Environment

Page 21 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC Graphic

21 / 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. (a) Reaction pathways for CO2 reduction to formic acid (HCOOH). (b) The free energy diagram for CO2 reduction to HCOOH on Ni (111) surface at 0 V vs. RHE. Relevant free energy differences (∆G1, ∆G2, and ∆G3) are labeled. The same color scheme is used for (a) and (b). 73x36mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 25

Page 23 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Free energy change between the reaction steps (∆G1: CO2 → COOH*, ∆G2: CO2 → HCOO*, ∆G3: HCOO* → HCOOH) on (a) Ni-based NSAs, Cu (211) and Ni (111) surfaces; (c) strained Cu*/Ni as a function of the free energy of HCOO*. The NSAs above zero meet the two selectivity criteria as indicated by solid symbols. The corresponding volcano plots via the route 2 are shown in (b) and (d). Schematics for the overlayer (M*/Ni) and the subsurface (Ni/M) NSAs are shown in the insets of (b). Blue and gray spheres denote the M atoms and Ni, respectively. 109x85mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Free energy diagrams for CO2 electroreduction to HCOOH on (a) Ni/Ti, (b) Cu*/Ni, and (c) strained (εXY = +3%) Cu*/Ni NSAs at 0 V vs. RHE. 105x74mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 24 of 25

Page 25 of 25

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4. Free energy diagram for HER on various NSA surfaces at 0 V vs. RHE. The adsorption geometries are shown in the insets. The pink and grey spheres represent H and metals, respectively. 76x58mm (300 x 300 DPI)

ACS Paragon Plus Environment