Expanding Multinary Selenide Based High ... - ACS Publications

Apr 19, 2019 - Emily Johnson,. ‡ and Manashi Nath*,†. †. Department of Chemistry, Missouri University of Science & Technology, 400 W. 11th Stree...
0 downloads 0 Views 2MB Size
Subscriber access provided by OCCIDENTAL COLL

Article

Expanding Multinary Selenide based High-Efficiency Oxygen Evolution Electrocatalysts through Combinatorial Electrodeposition: Case Study with Fe-Cu-Co Selenides Xi Cao, Emily Johnson, and Manashi Nath ACS Sustainable Chem. Eng., Just Accepted Manuscript • DOI: 10.1021/ acssuschemeng.9b01095 • Publication Date (Web): 23 Apr 2019 Downloaded from http://pubs.acs.org on April 24, 2019

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

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 43 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

ACS Sustainable Chemistry & Engineering

Expanding Multinary Selenide based HighEfficiency Oxygen Evolution Electrocatalysts through Combinatorial Electrodeposition: Case Study with Fe-Cu-Co Selenides Xi Caoa, Emily Johnsonb and Manashi Nath*a a

Department of Chemistry, Missouri University of Science & Technology, 400 W 11th Street,

Rolla, Missouri 65409, United States. b

Department of Chemical and Biochemical Engineering, Missouri University of Science &

Technology, 1101 N. State Street, Rolla, Missouri 65409, United States. *

Email: [email protected]

KEYWORDS (Mixed metal selenides, Oxygen evolution electrocatalyst, Water splitting, Electrocatalyst)

ACS Paragon Plus Environment

1

ACS Sustainable Chemistry & Engineering 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 43

ABSTRACT

Developing low-cost electrocatalysts with high-efficiency for water splitting is a critical task to make this technology viable for large-scale clean energy generation. Transition metal selenides, comprising of earth abundant elements, such as Fe, Co and Cu, have gained attention as superior electrocatalysts for oxygen evolution reaction (OER) in the alkaline medium. In this article, we have systematically investigated the evolution of OER catalytic activity as a function of composition for a series of Fe-Co-Cu quaternary selenides by exploring a trigonal phase diagram. The OER activity was dependent on the quantity of Cu and Fe in the Co-Fe-Cu-Se quaternary selenide electrocatalysts, while surprisingly, Fe-Cu ternary selenides exhibit reduced OER activity in comparison to their pure parent compounds FeSe and Cu3Se2. Quaternary selenides exhibited more efficient catalytic activity with increasing amount of Fe or Cu in the catalysts, and the quaternary mixed metal selenide thin film of composition (Fe0.48Co0.38Cu0.14)Se showed the best catalytic performance with a small overpotential of 256 mV at 10 mA cm-2 and a low Tafel slope of 40.8 mV dec-1 in N2-saturated 1.0 M KOH solution. The outstanding catalytic performance of quaternary selenides may be explained by the possible electron cloud delocalization among the transition metal sites in the catalytic system through d-bands, leading to lower charge transport resistance at the catalyst-electrolyte interface as well better film conductivity, as has also been observed through electrochemical impedance spectroscopy (EIS). Such enhanced charge transfers eventually facilitate the rate of O2 release from the catalyst surface leading to enhanced activity.

ACS Paragon Plus Environment

2

Page 3 of 43 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

ACS Sustainable Chemistry & Engineering

Introduction Hydrogen has been widely considered as an ideal alternative for sustainable energy resource to diminish the utilization of fossil fuels, especially in recent times when fossil fuel depletion has reached an all-time high along with the environmental pollution caused by carbon dioxide, a byproduct of fossil fuel combustion.1,2,3 While hydrogen fuel holds tremendous promise as clean source of energy, the production of pure hydrogen is a challenge that must be met before hydrogen fuel economy can replace the mature conventional technologies. Alkaline water electrolysis is a well-established technology to produce high-purity hydrogen in copious quantities for use in fuel cells and other energy storage devices.4,5 Water splitting includes two critical half-cell reactions: hydrogen evolution reaction (HER) at the cathode (2H2O + 2e- → H2 + 2OH-) and oxygen evolution reaction (OER) at the anode (4OH- → O2 + 4e- + 2H2O).6,7,8,9 Among these, OER is the more challenging kinetically sluggish process with multistep four-electron transfers involved which has greatly reduced the efficiency of

water splitting reaction thereby, limiting the

applicability of this technology to produce hydrogen on a mass scale. An efficient electrocatalyst is required to speed up OER by lowering the activation barrier and reducing the overpotential required for water oxidation. In this context, large number of research endeavors have been devoted to designing high-efficiency and low-cost OER electrocatalysts.10,11,12,13,14 Although research over the last several decades have identified some state-of-the-art OER electrocatalysts based on noble-metal-based oxides, such as IrO2 and RuO2 exhibiting low overpotentials and excellent stability, the scarcity and high cost of these precious metals has hindered their large scale industrial applications.15,16 Consequently, over the past several years, researchers have focused on exploring earth abundant transition metal based chemistry to identify and develop low-cost electrocatalysts comprising first row transition metals (primarily Fe, Ni, and Co). The quest for

ACS Paragon Plus Environment

3

ACS Sustainable Chemistry & Engineering 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 4 of 43

these transition metal electrocatalysts has led to oxides and (oxy)hydroxides, that could potentially replace precious metals due to their comparable activity and durability.17,18,19,20,21,22,23,24 Continued research conducted by our group as well as by other groups over the last couple of years has made it quite apparent that transition-metal based selenides are more likely to exhibit an improved OER catalytic activity and a reduced overpotential compared to oxides. Such enhancement of OER catalytic activity is attributed to the increased covalency of the metal-ligand bond as well as lowering the electrochemical redox potential of the catalytically active transition metal site.25,26,27,28,29,30 It was also observed that the smaller electronegativity of selenium compared to oxygen led to enhancement of preferential attachment of hydroxyl ions on the catalytically active metal sites and facilitated further steps of the OER process.31,32,33 Accordingly, binary metal selenides (Ni3Se2, NiSe2, CoSe, and Co7Se8) have been widely investigated as excellent OER electrocatalysts with impressively low overpotentials.26,34,35,36 Ternary mixed metal selenides such as NixFe1-xSe2, Ni0.89Co0.11Se2, FeNi2Se4 and NiCoSe2 exhibits more efficient catalytic activity suggesting that introduction of transition metal dopants improves catalytic activity.37,38,39,40 Further studies revealed that transition metal doping increased the number of actual catalytically active sites and accelerated the rate determining steps by modulating the OH- adsorption kinetics on the catalyst surface through tuning the local electron density around the catalytic site.41,42 Similar effects can also be observed in the transition metal oxides/(oxy)hydroxides, such as Fe or Co doped NiO/NiOOH exhibiting a much higher catalytic activity than NiO/NiOOH itself.43,44,45,46 More recently, we have systematically studied how transition metal doping in catalytically active selenide compositions effect OER activity through the exploration of a trigonal phase diagram of Ni-Fe-Co selenides and demonstrated that quaternary mixed metal selenides exhibited better OER activity than binary or even ternary selenides.33 In the Ni-Fe-Co selenide system, quaternary

ACS Paragon Plus Environment

4

Page 5 of 43 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

ACS Sustainable Chemistry & Engineering

selenides with progressively increasing amount of Fe showed higher catalytic activity and the optimal composition was identified as (Ni0.25Fe0.68Co0.07)3Se4 which could achieve 10 mA cm-2 at a small overpotential of 230 mV and showed remarkable durability in alkaline medium.33 While Ni and Co have traditionally shown excellent OER activity and have been given primary attention due to their relative abundance on the earth’s crust compared to precious metals, another abundant element with high d-electron occupancy, namely copper, has been explored far less for electrocatalytic water splitting. Cu-based OER electrocatalysts will be even more impactful since copper is one of the most abundant elements on the earth’s crust found in lots of minerals.47 Very recently there has been a couple of reports where Cu2O has been reported as an OER active catalyst.48,49 We have synthesized Cu2Se via electrodeposition as well as hydrothermal and high temperature CVD methods and have shown that Cu2Se exhibited better electrocatalytic activity compared to the oxide analogue.50 Additionally, introduction of Cu with highly occupied d-levels impart enhanced conductivity to the matrix which will improve charge transport on the catalyst surface. Hence, in order to understand the effect of transition metal doping in copper selenide on the OER catalytic activity, we have performed a systemic study of the phase diagram of Co-Ni-Cu selenides.51 These studies reveal that OER catalytic activity is influenced by the amount of Cu introduced into the catalytic system, specifically for quaternary selenides, and the best-performing catalyst composition in this case was identified as (Co0.21Ni0.25Cu0.54)3Se2 thin film which required an overpotential of 272 mV to reach 10 mA cm-2.51 Fe has been shown to act as an enhancer for OER catalytic activity in oxides and selenides,33,52,53 and quaternary selenides incorporating both Cu and Fe are expected to have better OER catalytic activity. Inspired by this idea, we have systematically investigated a series of mixed metal selenides comprising various amounts of Fe, Co and Cu by exploring a trigonal phase diagram through combinatorial electrodeposition as has

ACS Paragon Plus Environment

5

ACS Sustainable Chemistry & Engineering 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 6 of 43

been reported in this article. Surprisingly, Fe-Cu ternary selenides showed reduced OER activity in comparison to their pure parent compounds FeSe and Cu3Se2, indicating that neither Fe nor Cu doping can improve the activity of Cu3Se2 or FeSe, respectively. On the other hand, with introduction of Co dopant into the composition, all Fe-Co-Cu quaternary selenides showed enhanced catalytic activity with low overpotential as well as lower Tafel slopes. The catalytic activity of the quaternary selenides improved as the relative amounts of either Fe or Cu increased in the catalysts. The optimal catalyst composition in this phase space has been identified as (Fe0.48Co0.38Cu0.14)Se which required an overpotential of 256 mV to achieve 10 mA cm-2 and a showed Tafel slope of 40.8 mV dec-1. The low value of Tafel slope indicates that the catalyst has favorable kinetics for OER in alkaline medium. It is interesting to note that even in the absence of Ni, which has been recognized as one of the most active catalytic sites, it is possible to have highefficiency catalytic activity for OER through transition metal doping. The enhanced catalytic activity can be attributed to the possible electron cloud delocalization through formation of dbands between the different transition metal centers within the catalytic system.7,54 This systematic exploration of mixed-metal selenides comprised of iron, cobalt and copper offers more opportunities to understand the evolution of catalytic property as a function of covalency and delectron occupancy of the transition metals and will have far-reaching implications in practical application of electrocatalytic water splitting. Experimental Section Materials: The metal precursors comprising iron (II) sulfate [FeSO4·7H2O, 99.0 %], cobalt (II) sulfate [CoSO4·7H2O, 99.0 %], copper (II) sulfate [CuSO4·5H2O, 99.0 %], selenium dioxide [SeO2], ammonium sulfate [(NH4)2SO4, 99 %], potassium hydroxide [KOH, 85.0 %], were obtained from Fisher Scientific and used as is. All solutions were prepared in deionized (DI) water

ACS Paragon Plus Environment

6

Page 7 of 43 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

ACS Sustainable Chemistry & Engineering

(18 MΩ / cm). The Au-coated glass substrates used in electrodeposition were purchased from Deposition Research Lab Incorporated (DRLI), Lebanon, Missouri. Combinatorial electrodeposition of thin films. Systematic exploration of a trigonal phase diagram was carried out to synthesize a series of transition metal doped ternary and quaternary selenide films through combinatorial electrodeposition as shown in Figure S1. This method is very similar to our previously reported work featuring Co-Fe-Ni selenides, and the process has been described in details in the previous publication.33 Specifically, the electrolyte contained a combination of freshly prepared 10 mM of FeSO4·7H2O, CoSO4·7H2O, and CuSO4·5H2O solutions along with 10 mM SeO2 and 25 mM (NH4)2SO4. Before electrodeposition, all electrolytes were purged in N2 for 30 mins to remove dissolved oxygen. Prior to electrodeposition, Au-coated glass substrates were sonicated in micro-90, isopropanol followed by rinsing with deionized water to obtain clean surface.33,51 To ensure same geometric area for all electrodeposited thin films studied in this article, the cleaned substrates were confined with Teflon tape punched with a hole of known geometric area. An IviumStat potentiostat was used for all the electrodeposition experiments. The films were deposited on Au-coated glass working electrode using a threeelectrode electrochemical cell consisting of Ag|AgCl (KCl saturated) reference electrode and glassy carbon (GC) counter electrode. The electrodeposition potential was maintained at –0.8 V 𝑣𝑠 Ag|AgCl (KCl saturated) for 300 s at room temperature. Characterization. KRATOS AXIS 165 X-ray Photoelectron Spectrometer along with the monochromatic Al X-ray source were used to conduct X-ray photoelectron spectroscopy (XPS) of the catalyst. All XPS analysis was collected on the as-prepared sample surface without any sputtering. The morphologies of the thin film before and after electrochemical activity were observed with scanning electron microscopy (SEM) using the FEI Helios NanoLab 600

ACS Paragon Plus Environment

7

ACS Sustainable Chemistry & Engineering 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 8 of 43

FIB/FESEM at an accelerating voltage of 15 kV and a working distance of around 5.0 mm. Energy dispersive spectroscopy (EDS) from the same SEM microscope was used to evaluate the elemental compositions of all electrodeposited films. In this work, the elemental analysis was obtained at three to five different spots on the surface of each compound and the average EDS data was considered for obtaining the relative atomic ratio of the constituent elements. For all compounds, the molecular formula was written as (FexCoyCuz)Sen, where x, y, z, and n were the EDS atomic ratios of the transition elements to total amount of selenium, respectively. Electrochemical Measurements A typical three-electrode electrochemical cell connected to an IviumStat potentiostat was used for measuring the electrocatalytic performances. Linear sweep voltammetry (LSV) was carried out to study the OER catalytic activity, while the catalyst stability was investigated by carrying out constant potential chronoamperometry for an extended period of time in N2-saturated 1.0 M KOH at room temperature. In all the electrochemical experiments, the deposited thin films served as the working electrode, while GC was used as the counter electrode and Ag|AgCl (KCl saturated) as the reference electrode. Nernst equation as shown in eq. 1 was used to convert the experimentally measured potentials vs Ag|AgCl (KCl saturated) to the reversible hydrogen electrode (RHE) . 𝑜 𝐸𝑅𝐻𝐸 = 𝐸𝐴𝑔|𝐴𝑔𝐶𝑙 + 0.059𝑝𝐻 + 𝐸𝐴𝑔|𝐴𝑔𝐶𝑙

(1)

The standard potential of Ag|AgCl (KCl saturated) at 25.1 ºC was 0.197 V. All the potentials reported in this article have been iR corrected. Calculation of Tafel plot Tafel equation relating the overpotential 𝜂 with the current density 𝑗 as shown in eq. 2 was used to study the electrochemical kinetics of the (FexCoyCuz)Sen thin films for OER . 𝜂=𝑎+

2.3𝑅𝑇 𝛼𝑛𝐹

log(𝑗)

(2)

ACS Paragon Plus Environment

8

Page 9 of 43 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

ACS Sustainable Chemistry & Engineering

Here 𝜂 refers to the overpotential, 𝑗 represents the current density. The Tafel slope is given by2.3𝑅𝑇/𝛼𝑛𝐹. The Tafel plots in this work were calculated from the LSV collected at a scan rate of 2 𝑚𝑉 𝑠 −1 in a non-stirred N2-saturated 1.0 M KOH solution. Electrochemically active surface area (ECSA) Electrochemical double layer capacitance was used to measure the ECSA of these catalysts following standard methods described in detail in previous reports.17,33 In the measurement of ECSA through double layer capacitance it is assumed that the current obtained in the non-Faradaic region is caused by double-layer charging instead of electrochemical reactions or charge transfer. The double layer current (𝑖𝐷𝐿 ) was obtained by performing cyclic voltammograms (CVs) with various scan rates in a non-stirred N2-saturated 1.0 M KOH solution. The ratio of double-layer current 𝑖𝐷𝐿 and the scan rate (ν) of CV yielded specific electrochemical double-layer capacitance (𝐶𝐷𝐿 ) as shown in eq. 3. The cathodic and anodic charging currents were a function of scan rates and the 𝐶𝐷𝐿 can be calculated by averaging the absolute values of cathodic and anodic slopes. Eq. 4 can be used to calculate ECSA. 𝑖𝐷𝐿 = 𝐶𝐷𝐿 × 𝑣

(3)

𝐸𝐶𝑆𝐴 = 𝐶𝐷𝐿 ⁄𝐶𝑠

(4)

Here, double-layer capacitance (𝐶𝐷𝐿 ) is in units of mF. Based on other reported values, the specific capacitance (Cs) was assumed to be 0.040 mF cm-2 in the alkaline solution.17,26 Roughness factor (RF) was estimated from the ratio of ECSA and the geometric area (0.07 cm2). Results and Discussion Structural and Elemental Compositions of Electrodeposited Thin Films. The trigonal phase diagram as shown in Figure S1 was systematically explored to electrodeposit series of mixed metal selenides. This method referred to as combinatorial electrodeposition is similar to our previous

ACS Paragon Plus Environment

9

ACS Sustainable Chemistry & Engineering 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 43

published reports where the process has been described in detail with different transition metals.33,51 In the trigonal phase diagram, the vertices represented three binary selenides (CoSe, FeSe and Cu3Se2), while the sides of the triangle referred to the ternary selenides FexCoySen, FexCuySen and CoxCuySen, respectively. The interior of this trigonal phase diagram represented the quaternary mixed-metal selenides (FexCoyCuzSen). It should be noted here that all the electrodeposited thin films of ternary and quaternary compositions were amorphous in nature and did not yield clear pxrd patterns making proper phase determination from pxrd alone less probable. The elemental composition of the electrodeposited films on the other hand, were determined from the average relative atomic ratio obtained from EDS analysis as listed in the Table S1, juxtaposed with the precursor ratio in the electrolyte that produced the respective film. Typically, EDS was collected at several regions on the sample surface and interestingly, it showed similar quantification results for the elemental ratio, indicating that while the films were not crystalline, they were indeed of uniform composition throughout each catalytic surface. As can be seen in Table S1, three as-prepared binary thin films were composed of Fe, Co or Cu and Se with approximate atomic ratio of 1 : 1 (Fe : Se), 1 : 1 (Co : Se) and 3 : 2 (Cu : Se) respectively, hence the molecular formulae of three binaries were written as FeSe, CoSe and Cu3Se2 as shown in Table S2, respectively. The molecular formulae identified in the various current voltage plots of thr ternary and quaternary selenides can be found in Table S2. Careful observation of Table S1 revealed that increasing (or decreasing) amounts of the metal precursor in the electrolytic bath led to an increase (or decrease) of the relative atomic percentage of the corresponding element in the electrodeposited film. This demonstrates the ability to control thin film composition through combinatorial electrodeposition.

ACS Paragon Plus Environment

10

Page 11 of 43 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

ACS Sustainable Chemistry & Engineering

Electrochemical Characterization. Linear sweep voltammetry (LSV) was carried out in N2saturated 1.0 M KOH at a scan rate of 10 mV s-1 to measure OER catalytic activity of the deposited films . The electrocatalytic performance was benchmarked in comparison with other highly active catalysts by comparing the onset overpotential (onset η) and the overpotential required to reach a current density of 10 mA cm-2 (η at 10 mA cm-2) obtained from the LSV data.17 Table S1 lists the onset overpotential, overpotential at 10 mA cm-2 and the corresponding EDS atomic ratios of all the electrodeposited selenides reported in this work. The OER polarization curves of the three binary selenide films can be found in Figure S2, demonstrating the OER catalytic activity trend of binary films was CoSe > Cu3Se2 > FeSe. The LSV plots for all the ternary selenide films as shown in Figure 1. In Figure 1a, Co-Cu selenides were obviously the best group in comparison to the other two, and the optimal composition in this group was identified as (Co0.93Cu0.07)3Se2 which onset overpotential was 235 mV and overpotential was 272 mV at 10 mA cm-2. Interestingly, the LSV plots of Fe-Cu selenide group in Figure 1b displayed the worst catalytic activity among the three ternary groups. The ternary composition in this group exhibiting the best catalytic activity with an onset overpotential of 307 mV and an overpotential of 347 mV to reach a current density of 10 mA cm-2 was identified as (Fe0.78Cu0.22)3Se4. Some of the LSV plots from the Co-Fe selenide group shown in Figure 1c was also included in our previous work on combinatorial electrodeposition of Ni-Co-Fe multinary selenides.33 In this group, (Co0.90Fe0.10)3Se4 exhibited the lowest onset overpotential of 257 mV and overpotential of 295 mV at 10 mA cm-2.33 It is worth mentioning that several pre-oxidation peaks were observed in these LSVs prior to OER process corresponding to oxidation of transition metal ions within the catalyst composite as shown in Figures 1a and 1d. Pre-oxidation peaks are commonly noted in the Ni- and Co- based electrocatalysts, which result from the change in oxidation state of transition metal during the

ACS Paragon Plus Environment

11

ACS Sustainable Chemistry & Engineering 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 43

anodic scan.55,56,57,58 As can be seen in Figures 1a and 1d, the peaks observed around 1.10 V vs RHE were attributed to the conversion of Co2+ to Co3+.33,57,59 Interestingly, the intensity of Co preoxidation peaks increased with increasing amount of Co in the electrocatalysts, which further confirmed combinatorial electrodeposition as a feasible method for compositional control. The polarization curves of the most efficient ternary selenide composition from each group has been compared with those of the binaries in Figure 1d to illustrate the effect of transition metal doping on OER catalytic activity. It was observed that doping CoSe with 7 % Cu improved the OER catalytic activity by decreasing η at 10 mA cm-2 by about 26 mV. Similarly, 22 % Cu dopant in FeSe lowered the overpotential by 46 mV to reach a current density of 10 mA cm-2. However, it should be noted that the catalytic activity of (Fe0.78Cu0.22)3Se4 was still observed to be worse than the binary copper selenide, Cu3Se2. Doping of of 10 % Fe into CoSe did not lead to any observable change in the OER catalytic activity. It has been widely reported that incorporation of transitionmetal dopants improved the OER catalytic activity, for example, small amount of Fe doped into Ni- or Co- based electrocatalysts can radically ameliorate the OER activity.43,60,61,62,63,64,65,66 It is interesting to note that in these ternary metal selenides, transition metal doping forces the formation of a different stoichiometry and structure type compared to the binary selenide. The ternary films of Co-Cu selenides formed the M3Se2 composition, and accordingly they could be described as Co-doped Cu3Se2. On the other hand, the Co-Fe selenide group formed the structure type M3Se4. Interestingly, the Fe-Cu group exhibits different stoichiometries at different relative ratios of Fe and Cu. At high Cu concentration with minimal Fe doping, the deposition adopts M3Se2 composition and can be described as doped Fe-doped Cu3Se2. At high relative ratio of Fe with minimal Cu doping, the structure can be described as M3Se4 (Cu-doped Fe3Se4). At comparable relative ratios of Fe:Cu, the stoichiometries obtained were closer to monoselenides, and can be

ACS Paragon Plus Environment

12

Page 13 of 43 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

ACS Sustainable Chemistry & Engineering

described as MSe. It should also be noted that each of these stoichiometries are distinctly different from each other with respect to crystal structure and may also contain different coordination geometry around the transition metal ions, as discussed later. Different stoichiometries of the metal selenides also give rise to different catalytic properties as has been observed with the Ni-selenide,

(a)

(c)

(b)

(d)

Figure 1. OER electrocatalytic performances of the various ternary metal selenides obtained along the edges of the trigonal phase diagram measured through linear sweep voltammetries (LSVs). (a) Co-Cu group (b) Fe-Cu group, and (c) Co-Fe group. (d) OER activity of the most active ternary composition (solid line) from each group was compared with the three binary compositions (dotted line). The dashed black line in each panel indicates the current density at 10 mA cm-2.

ACS Paragon Plus Environment

13

ACS Sustainable Chemistry & Engineering 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 14 of 43

Ni3Se2, NiSe, NiSe2 and Ni3Se4.26,27,67,68 Additionally, for mixed metal catalysts, it has been shown that the OER activity depends not only on the quantity of transition metal dopants but it is also strongly related to the location of metal dopants in the catalytic system. D. Friebel et al. have reported that Fe in Ni1-xFexOOH showed 500-fold increased OER activity when Fe was doped into the octahedral sites where the edge-sharing of [NiO6] and [FeO6] occurred affecting the local electronic structure of the catalyst.69 We have also previously reported through DFT calculations that introducing dopants into the surface layers altered the adsorption energy of hydroxyl species on the active metal sites thereby lowering the potential to initiate OER and leading to enhancement of the catalytic activity.33 Hence, it is not surprising to see Fe-Cu selenides showing less efficient OER catalytic activity, which can be a reflection of the inherent property of the structure type and composition. The trend of OER catalytic activity in the family of Co-Fe-Cu-Se quaternary selenides comprising of 66 individual compositions, has been shown in the contour plots of Figure 2, created from the two overpotential values (onset and η at 10 mA.cm-2) in Table S1 and the EDS atomic ratios. Contour plots were constructed to investigate the trend of OER onset potential and overpotential at 10 mA cm-2 as a function of composition as shown in Figures 2 a and 2b, respectively. The range of overpotentials obtained has been represented with a color gradient with the catalytic activity becoming better from bottom to top of the gradient scale corresponding to lower value of the overpotential. Typically, regions with red color showed better catalytic activity compared to the blue colored regions which showed the highest overpotentials. The very interior region in the trigonal phase diagram showed the lowest overpotential as shown in Figure 2a, demonstrating the best catalytic activity, while the areas adjacent to three axes exhibited less efficient catalytic activity. This demonstrated that heavily doped quaternary selenides were most

ACS Paragon Plus Environment

14

Page 15 of 43 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

ACS Sustainable Chemistry & Engineering

likely to exhibit higher OER catalytic activity in comparison to binary or ternary selenides in this Cu-Co-Fe phase space. As can be seen, the best quaternary selenides were confined to the region where Cu was around 10 – 50 %, Co was about 15 – 50 %, and Fe was roughly 30 – 70 %. The trends indicated that the quaternary selenides showed the best catalytic performance when the level of Fe was higher, suggesting that while Fe played an important role in improving the catalytic activity, it probably was not the main catalytic site, since Fe-rich ternary phases showed poor catalytic activity as can be observed in the Co-Fe group (Figure 1c). Meanwhile, amount of Cu effected the catalytic activity as well with the OER activity being enhanced by modest amount of Cu. It should be noted that a relatively higher amount of Fe or Cu compared to that of Co in the catalyst is more likely to produce better activity with lower overpotentials. However, as explained above, Fe-Cu-Se ternary selenide still showed poor performance indicating that Fe and Cu sites by (a)

(b)

Figure 2. The trend of OER catalytic activity within the entire Fe-Co-Cu trigonal phase space. Mapping of overpotential η (in units of V) (a) at the onset of OER activity and (b) at a current density of 10 mA cm-2. These contour plots were created by using the actual relative atomic ratio of the elements as obtained from EDS. The color gradient represents the overpotential measured in volts.

ACS Paragon Plus Environment

15

ACS Sustainable Chemistry & Engineering 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 16 of 43

themselves do not have higher catalytic activity, but the inductive effect of these ions on the neighboring Co sites increased the catalytic activity of Co thereby lowering the overpotentials. These observations further confirmed the benefits of such systematic studies in discovering untargeted mixed-metal based electrocatalysts. Figure 3 indicated the LSV polarization curves of the six best performing quaternary mixed metal selenides measured in N2-saturated 1.0 M KOH at scan rate of 10 mV s-1. These six compounds showed excellent catalytic activity for the OER process. The onset η of (Fe0.09Co0.29Cu0.62)3Se2 was 248 mV and it needed a η of 272 mV to reach a current density of 10 mA cm-2. The onset η and the η at 10 mA cm-2 of (Fe0.25Co0.61Cu0.14)Se were 258 mV and 278 mV, respectively. The onset η of (Fe0.29Co0.51Cu0.2)Se was 256 mV and it generated a current density of 10 mA cm-2 and a η of 278 mV. The onset η of (Fe0.44Co0.41Cu0.15)Se was 250 mV and required 271 mV to achieve a current density of 10 mA cm-2. The onset η and the η at 10 mA cm-2 of (Fe0.55Co0.26Cu0.19)Se were 255 mV and 276 mV respectively. The thin film (Fe0.48Co0.38Cu0.14)Se displayed the best catalytic activity in this work, with an onset η of 235 mV and η of 256 mV at 10 mA cm-2. It is important to note that most of the quaternary compounds adopted the monoselenide MSe structure type except for the ones with high Cu concentration, which adopted the M3Se2 structure type. The Tafel slopes of all electrocatalysts were obtained and has been listed in Table S1. The Tafel plots (η vs log j) and slopes of the six best performing quaternary selenides has been graphed in Figure

3b.

The

Tafel

slopes

of

(Fe0.55Co0.26Cu0.19)Se,

(Fe0.09Co0.29Cu0.62)3Se2,

(Fe0.44Co0.41Cu0.15)Se, (Fe0.48Co0.38Cu0.14)Se, (Fe0.29Co0.51Cu0.2)Se and (Fe0.25Co0.61Cu0.14)Se were 45.2, 59.5, 41.1, 40.8, 41.8, and 41.3 mV dec-1, respectively. Tafel slopes are closely related to the number of electrons involved in the rate determining step of the OER process. For the six

ACS Paragon Plus Environment

16

Page 17 of 43 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

ACS Sustainable Chemistry & Engineering

quaternaries with small Tafel slopes, it is likely that the second electron transfer step of OER is the rate-determining step suggesting that they have faster kinetics for oxygen evolution.7,70 (Fe0.48Co0.38Cu0.14)Se showed the lowest overpotential as well as lowest Tafel slope which further confirmed that it is indeed the best catalyst composition in this quaternary phase space. Electrochemical double-layer capacitance of (Fe0.48Co0.38Cu0.14)Se thin film was measured to estimate its electrochemically active surface area (ECSA). Double-layer charging current was measured in the non-Faradaic region at various scan rates using cyclic voltammograms (CVs). In this study, (Fe0.48Co0.38Cu0.14)Se thin film as the working electrode was suspended in N2-saturated 1.0 M KOH and the CVs, found in Figure 4, were measured from -0.16 to -0.05 V vs Ag|AgCl (KCl saturated) at scan rates of 2.5 to 20 mV s-1. The ECSA was calculated using eq. 4, and the value was estimated to be 1.38 cm2 with a roughness factor (RF) of 19.7. The high value of RF confirmed the highly granular surface of the film as seen in the SEM image of Figure 5b. Since rough surfaces allow for amplified contact of the electrolyte with the catalytically active sites, catalytic activity tends to increase with increasing surface roughness. (a)

(b)

Figure 3. (a) OER activity of the six best quaternary compositions measured by conducting LSVs in N2-saturated 1.0 M KOH solution at a scan rate of 10 mV s-1. (b) Tafel plots of these catalysts.

ACS Paragon Plus Environment

17

ACS Sustainable Chemistry & Engineering 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 18 of 43

The catalytic activity of the best performing quaternary (Fe0.48Co0.38Cu0.14)Se was compared with those of CoSe, FeSe, and Cu3Se2 in N2-saturated 1.0 M KOH at a scan rate of 10 mV s-1 to understand the effect of transition metal doping on catalytic activity as shown in Figure 5a. The onset overpotentials of the binary selenides were 263, 346 and 270 mV for CoSe, FeSe, and Cu3Se2 respectively. To achieve an OER current density of 10 mA cm-2, the overpotentials were 298 mV, 393 mV and 330 mV for CoSe, FeSe and Cu3Se2, respectively. As a result of transition-metal doping, the onset η and η at 10 mA cm-2 of (Fe0.48Co0.38Cu0.14)Se were decreased to 235 mV and 256 mV, respectively which were lower than any of the binary selenides. The Tafel slope of (Fe0.48Co0.38Cu0.14)Se (40.8 mV dec-1) was also lower than that of CoSe (62.9 mV dec-1), FeSe (103.5 mV dec-1) and Cu3Se2 (93.4 mv dec-1) as shown in Table S1 suggesting the OER activity

Figure 4. ECSA for the (Fe0.48Co0.38Cu0.14)Se thin film was measured by performing cyclic voltammograms in N2-saturated 1.0 M KOH solution at different scan rates from 2.5 to 20 mV s-1. The inset shows the plot of both anodic and cathodic currents measured at -0.10 V vs Ag|AgCl (KCl saturated) as a function of scan rate.

ACS Paragon Plus Environment

18

Page 19 of 43 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

ACS Sustainable Chemistry & Engineering

on the quaternary selenide surface is increased kinetic favorability compared to that on the binaries. This can be attributed to transition-metal doping, which increases the amount of actual catalytically active sites as well as influence the energetics of hydroxyl attachment, thereby facilitating the OER reaction rate.7,14,71,72,73 The LSV curve of (Fe0.48Co0.38Cu0.14)Se was also compared with the catalytic performance of the best ternary phases from each group, viz. (Fe0.78Cu0.22)3Se4, (Co0.93Cu0.07)3Se2, (Co0.9Fe0.1)3Se4 shown in Figure S3. As can be seen, (Fe0.48Co0.38Cu0.14)Se exhibited ameliorated OER activity compared to the best ternary composition (Co0.93Cu0.07)3Se2 with a decrease in η at 10 mA cm-2 by 16 mV. Overall, the OER catalytic activity was decidedly enhanced in mixed-metal selenides in terms of both overpotential and Tafel slope. The morphology and composition of the best performing quaternary selenide was further studied through surface analytical techniques. Detailed SEM imaging of the (Fe0.48Co0.38Cu0.14)Se film as shown in Figure 5b suggested that the composition of the thin film surface morphology was composed primarily of cauliflower-like microstructures with a wide size distribution of 0.5 to 2.5 µm. X-ray photoelectron spectroscopy (XPS) was carried out to further study the elemental composition and identify oxidation states of the constituent elements in (Fe0.48Co0.38Cu0.14)Se thin film. The XPS spectra of Fe 2p, Cu 2p, Co 2p and Se 3d have been shown in the Figure 5c-f. It should be noted that the binding energy in all XPS spectra has been calibrated with respect to the standard C 1s (284.5 eV) reference spectra. The spectrum of Fe 2p in Figure 5c showed peaks at 708.45 and 722.27 eV corresponding to Fe2+ while peaks at 710.26 and 724.75 eV were assigned to Fe3+ with satellite peaks at 716.08 and 730.57 eV.74,75,76 In the Cu 2p spectrum as shown in Figure 5d, peaks at 930.32 and 950.12 eV were attributed to Cu+ and those at 932.16 and 952.23 eV corresponded to Cu2+, respectively, while peaks at 939.75 and 960.15 eV were assigned to the

ACS Paragon Plus Environment

19

ACS Sustainable Chemistry & Engineering 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 20 of 43

satellite peaks.77,78 In the Co 2p spectrum shown in Figure 5e, peaks at 779.14 and 794.48 eV (a)

(b)

500 nm

3 µm

(c)

(d)

(e)

(f)

Figure 5. (a) Polarization curves of (Fe0.48Co0.38Cu0.14)Se compared with the three binary catalyst films. (b) SEM image of (Fe0.48Co0.38Cu0.14)Se film deposited for 300 s. (c) – (f) XPS spectra of the as-deposited (Fe0.48Co0.38Cu0.14)Se showing, (c) Fe 2p (d) Cu 2p (e) Co 2p and (f) Se 3d.

ACS Paragon Plus Environment

20

Page 21 of 43 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

ACS Sustainable Chemistry & Engineering

referred to Co2+, and those at 783.80 and 800.61 eV corresponded to its satellite peaks, respectively.79,80 In the Se 3d spectrum shown in Figure 5f, peaks at 55.24 and 56.22 eV were assigned to Se 3d5/2 and Se 3d3/2 respectively, while the big sharp peak around 59.0 eV was attributed to SeOx which might be a result of the surface oxidation of selenide.81,82 The deconvoluted spectra revealed that both Fe 2p and Cu 2p were mixed-valence metal cations while Co 2p was present as only divalent ions. The XPS spectra also confirmed that the film composition was exclusively selenide with no indication of presence of other metal oxide/oxyhydroxide impurity. The Faradaic efficiency of the electrodeposited (Fe0.48Co0.38Cu0.14)Se thin film for OER was assessed with the water displacement method,83 and the measured quantity of O2 evolved was compared to the theoretical yield as seen in Figure 6. A Faradaic efficiency of nearly 100% was measured for this electrocatalyst.

Figure 6. The amount of theoretically calculated oxygen (black dots) and experimentally measured oxygen (red dots) obtained for (Fe0.48Co0.38Cu0.14)Se versus time at 0.7 V vs Ag|AgCl (KCl saturated).

ACS Paragon Plus Environment

21

ACS Sustainable Chemistry & Engineering 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 22 of 43

As mentioned above, transition metal doping is expected to influence electron density through inductive effect around the catalytic centers, and hence will influence charge transport at the catalyst-electrolyte interface as well as through the catalyst composite. Electrochemical impedance spectroscopy (EIS) was executed on the best performing quaternary composition along with the three binary selenide films to examine the result of transition metal doping on the various charge transfer pathways. Figure 7 shows the Nyquist plots of the best-performance quaternary selenide, (Fe0.48Co0.38Cu0.14)Se, compared with those of three binary selenides (FeSe, CoSe and Cu3Se2), which were collected at an applied potential of 0.5 V vs Ag|AgCl (KCl saturated) in N2-saturated 1.0 M KOH solution. The Nyquist plots were fitted to an equivalent circuit as shown in insert of Figure 7 which included the electrolyte resistance (Rs), interfacial electron charge transfer resistance (Rct), and the constant phase element of double layer (CPEdl). The values of these parameters have been listed in Table 1. In the higher frequency region, the low value of Rs of ca. 5.0 Ω was attributed to the uncompensated solution resistance of 1.0 M KOH (Figure S4). Derived from the equivalent circuit, the charge transfer resistance (Rct) and constant phase element of double layer (CPEdl) were obtained in the lower-frequency region. The charge transfer resistance Rct is related to the electron transfer across the electrode-electrolyte interface and hence represents the kinetics of the electrocatalysis and its value is inversely proportional to the reaction rate.84,85,86 From Table 1, it can be seen that (Fe0.48Co0.38Cu0.14)Se showed the smallest Rct value, indicating that, at the film-electrolyte interface, the charge transfer rate was kinetically faster for (Fe0.48Co0.38Cu0.14)Se than that of binaries (CoSe, Cu3Se2 and FeSe). Fast interfacial charge transfer leads to high catalytic efficiency of (Fe0.48Co0.38Cu0.14)Se by facilitating the rate of intermediate formation on the catalyst surface and decreasing the OER onset overpotential as a

ACS Paragon Plus Environment

22

Page 23 of 43 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

ACS Sustainable Chemistry & Engineering

result.87,88 Likewise, the catalytic activity can be influenced by the thin film resistance (Rf) by enhancing electron transfer within the bulk of the catalyst composite leading to increased current density. It can be observed that (Fe0.48Co0.38Cu0.14)Se showed the smallest Rf value in comparison to that of three binaries in Table 1, indicating the quaternary selenide had an enhanced conductivity and helped the local electron transfer in the catalytic system. From these studies it is very apparent

Figure 7. EIS measurements and the corresponding Nyquist plots for (Fe0.48Co0.38Cu0.14)Se, CoSe, Cu3Se2 and FeSe catalysts obtained at an applied potential of 0.5 V vs Ag|AgCl (KCl saturated) in N2-saturated 1.0 M KOH solution. Symbols indicate the raw data, while solid lines represent the corresponding fit to equivalent circuit model (inset).

ACS Paragon Plus Environment

23

ACS Sustainable Chemistry & Engineering 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 24 of 43

that transition metal doping improves surface electronic structure and thereby enhancing various charge transport pathways for OER catalytic activity. Table 1. Equivalent circuit parameters obtained from fitting of EIS experimental data. Catalysts

Rs/Ω

Rct/Ω

CPEdl/F

Rf /Ω

CPEf /F

(Fe0.48Co0.38Cu0.14)Se

5.0

10.9

0.24

113.1

0.80

CoSe

5.0

13.2

0.30

248.9

0.83

Cu3Se2

5.0

18.4

0.48

1022.3

0.87

FeSe

5.0

291.9

0.75

2040.0

0.93

Rs is the resistance of the electrolyte; Rct is the electron transfer resistance; CPEdl is the constant phase element of double-layer nonideal capacitance; Rf is resistance of the catalyst layer; and CPEf is the capacitance.

The stability of any OER electrocatalyst is a critical factor to evaluate its potential in large-scale industrial application. Hence, the stability of (Fe0.48Co0.38Cu0.14)Se thin film under continuous oxygen evolution for extended period of time (18 h) was estimated through chronoamperometric study as shown in the inset of Figure 8a. The chronoamperometric study was performed at a constant potential of 1.50 V vs RHE, where the thin film reached a current density of 10 mA cm-2 in 1.0 M KOH solution. It can be observed that (Fe0.48Co0.38Cu0.14)Se catalyst showed outstanding stability without any current density loss even after 18 h in the harsh alkaline medium. To further confirm

the

stability,

LSV

of

(Fe0.48Co0.38Cu0.14)Se

was

re-measured

after

18

h

chronoamperometry and was compared with that of the as-deposited catalyst shown in Figure 8a, which revealed the quaternary selenide showed almost the same catalytic activity before and after stability study with similar onset potential and overpotential at 10 mA cm-2. The SEM shown in Figure 8b exhibited that the surface morphology of (Fe0.48Co0.38Cu0.14)Se had been slightly modified by making the surface rougher possibly due to the evolution of large amounts of O2 on

ACS Paragon Plus Environment

24

Page 25 of 43

catalyst surface, resulting in the production of pores. The compositional stability of the catalyst 30

(a)

(b)

j / mA cm-2

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

ACS Sustainable Chemistry & Engineering

20 10 0 0

2

4

6

8 10 12 14 16 18

Time /h

2 µm

(c)

(d)

(e)

(f)

Figure 8. (a) Comparison of LSVs before and after 18 h in N2-saturated 1.0 M KOH solution. The inset is the chronoamperometry plot of (Fe0.48Co0.38Cu0.14)Se for 18 h at an applied potential of 1.50 V vs RHE. (b) SEM image of (Fe0.48Co0.38Cu0.14)Se after electrochemical measurements. (c−f) XPS spectra of (Fe0.48Co0.38Cu0.14)Se after chronoamperometry.

ACS Paragon Plus Environment

25

ACS Sustainable Chemistry & Engineering 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 26 of 43

after electrochemical measurement was corroborated using XPS spectra, seen in Figure 8c-f, which revealed that the catalyst was predominantly comprised of Fe, Co, Cu and Se and they were still in the mixed oxidation states. The XPS peak positions for Fe, Co, and Cu did not show any changes indicating that the composition of the film was still selenide-based. The comparison of EDS before and after chronoamperometry in Table S3 further demonstrated that the relative elemental ratio was preserved even after prolonged catalytic activity signifying compositional stability of the catalyst. Interestingly, while the ternary selenides adopt various non-stoichiometric compositions including M3Se4, MSe, and M3Se2 depending on the transition metal doping, the quaternary selenides are mostly forced to adopt the monoselenide, MSe type phase. The primary difference between the M3Se2 and MSe structure type can be illustrated with Cu3Se2 and CoSe as shown in Figure S5 in supporting information. In Cu3Se2, the transition metal is tetrahedrally coordinated to 4 Se atoms, and the CuSe4 tetrahedra are linked further through edge sharing of the polyhedra. CoSe on the other hand represents the fully filled NiAs structure type containing octahedrally coordinated Co atoms, while the Se atoms show tetrahedral coordination. On doping in CoSe lattice, it can be expected that the transition metal dopant atoms will occupy the octahedral sites, while similar dopants in Cu3Se2 will occupy the tetrahedral sites. From the trend in catalytic activity observed here, it is apparent that for Co and Fe, the transition metal in octahedral coordination shows better catalytic activity than that in tetrahedral coordination. It has been reported previously also that Fe in the octahedral site shows better enhancement of catalytic property.60,64,89 Cu on the other hand shows better activity in tetrahedral coordination. Apart from influence of structure type and coordination geometry, other factors can also play a crucial role in determining the OER catalytic property. For example, it was reported by previous researchers that

ACS Paragon Plus Environment

26

Page 27 of 43 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

ACS Sustainable Chemistry & Engineering

a catalyst surface containing anion vacancies and defects was favorable for the adsorption and desorption of hydroxyl species during the OER process, leading to small charge transfer resistance at the electrode-electrolyte interface and which further enhances the OER activity.90,91,92 Hence, the structural defects might also play a crucial role in affecting the catalytic activity of these amorphous electrocatalysts reported in this article. However, to comprehensively analyze the role of structural defects on the OER catalytic activity, one needs to perform detailed DFT studies for each stoichiometry and related structure type. We have initiated such DFT studies and have recently observed that for one particular stoichiometry based on Cu-Co-Se, identified from this phase diagram, the structural defects on different lattice planes indeed has a significant influence on the OER catalytic activity. The results of these DFT studies will be published separately. Conclusions The ternary phase space of Co-Fe-Cu has been successfully explored through combinatorial electrodeposition and have identified the optimal doping compositions that yield maximum catalytic efficiency in this family of compounds. EDS characterization of all selenides films revealed the ability of compositional control through combinatorial electrodeposition as the relative elemental composition of the thin films were proportional to the corresponding precursor ratios in the electrolytes. A full investigation of OER kinetic parameters from sixty-six compounds revealed the catalytic trends as follows: (i) the catalytic activity is sensitive to the amount of transition metal dopants, and the presence of relatively higher amount of Fe or Cu and lower amount of Co leads to small overpotentials, while overpotentials tend to decrease at Fe concentration 30 – 70 %; (ii) quaternary selenides are more likely to show lower overpotentials as well as Tafel slopes than ternary and binary selenides indicating the inductive effect of heteroatoms play a crucial role in influencing the catalytic activity; (iii) Co might still be better catalytic site

ACS Paragon Plus Environment

27

ACS Sustainable Chemistry & Engineering 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 28 of 43

although Fe and Cu have a more positive influence on increasing its activity; (iv) the enhancement in catalytic activity is supported by increased charge transfer at the catalyst-electrolyte interface and within the catalyst composite. Specifically, the best composition in this work has been identified as (Fe0.48Co0.38Cu0.14)Se thin film with a low onset overpotential and overpotential to attain a current density of 10 mA cm-2 of 235 mV and 256 mV, respectively in N2-saturated 1.0 M KOH demonstrating promising catalytic activity. This study brings out a new avenue for systematic study of mixed-transition-metal selenides as OER electrocatalysts. ASSOCIATED CONTENT AUTHOR INFORMATION Corresponding Author *

Email: [email protected]

Author Contributions All authors have given approval to the final version of the manuscript. Notes The authors declare no competing financial interest. Supporting Information. The Supporting Information includes trigonal phase diagram, LSV plots, summary table, EDS table before and after chronoamperometry, atomic percentage table and crystal structures. The Supporting Information is available free of charge on the ACS Publications website. ACKNOWLEDGMENT

ACS Paragon Plus Environment

28

Page 29 of 43 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

ACS Sustainable Chemistry & Engineering

The authors would like to acknowledge financial support from National Science Foundation (DMR 1710313), American Chemical Society Petroleum Research Fund (54793-ND10), and Energy Research and Development Center (ERDC) Missouri S&T. The authors would also like to acknowledge Dr. Wenyan Liu for help with the ICP-AES measurements and Materials Research Center (MRC) for equipment usage.

REFERENCES (1)

Chow, J.; Kopp, R. J.; Portney, P. R. Energy Resources and Global Development. Science 2003, 302, 1528–1531.

(2)

Turner, J. A. Sustainable Hydrogen Production. Science 2004, 305, 972–974.

(3)

Dresselhaus, M. S.; Thomas, I. L. Alternative Energy Technologies. Nature 2001, 414, 332– 337.

(4)

Lewis, N. S. Toward Cost-Effective Solar Energy Use. Science 2007, 315, 798–801.

(5)

Walter, M. G.; Warren, E. L.; Mckone, J. R.; Boettcher, S. W.; Mi, Q.; Santori, E. A.; Lewis, N. S. Solar Water Splitting Cells. Chem. Rev. 2010, 110, 6446–6473.

(6)

Rossmeisl, J.; Logadottir, A.; Nørskov, J. K. Electrolysis of Water on (Oxidized) Metal Surfaces. Chem. Phys. 2005, 319, 178–184.

(7)

Suen, N.-T.; Hung, S.-F.; Quan, Q.; Zhang, N.; Xu, Y.-J.; Chen, H. M. Electrocatalysis for the Oxygen Evolution Reaction: Recent Development and Future Perspectives. Chem. Soc. Rev. 2017, 46, 337–365

ACS Paragon Plus Environment

29

ACS Sustainable Chemistry & Engineering 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

(8)

Page 30 of 43

Dau, H.; Limberg, C.; Reier, T.; Risch, M.; Roggan, S.; Strasser, P. The Mechanism of Water Oxidation: From Electrolysis via Homogeneous to Biological Catalysis. ChemCatChem 2010, 2, 724–761.

(9)

Yan, Y.; Xia, B. Y.; Zhao B.; Wang, X. A Review on Noble-Metal-Free Bifunctional Heterogeneous Catalysts for Overall Electrochemical Water Splitting. J. Mater. Chem. A 2016, 4, 17587–17603.

(10) Görlin, Y.; Jaramillo, T. F. A Bifunctional Nonprecious Metal Catalyst for Oxygen Reduction and Water Oxidation. J. Am. Chem. Soc. 2010, 132, 13612–13614. (11) Chen, Q.; Switzer, J. A. Photoelectrochemistry of Ultrathin, Semitransparent, and Catalytic Gold Films Electrodeposited Epitaxially onto n-Silicon (111). ACS Appl. Mater. Interfaces 2018, 10, 21365–21371. (12)

Mahenderkar, N. K.; Chen, Q.; Liu, Y. C.; Duchild, A. R.; Hofheins, S.; Chason, E.; Switzer, J. A. Epitaxial Lift-off of Electrodeposited Single-Crystal Gold Foils for Flexible Electronics. Science 2017, 355, 1203–1206.

(13)

Wang, H.; Lee, H.; Deng, Y.; Lu, Z.; Hsu, P.; Liu, Y.; Lin, D.; Cui, Y. Bifunctional NonNoble Metal Oxide Nanoparticle Electrocatalysts through Lithium-Induced Conversion for Overall Water Splitting. Nat. Commun. 2015, 6, 7261.

(14)

Suntivich, J.; May, K. J.; Gasteiger, H. A.; Goodenough, J. B.; Shao-Horn, Y. A Perovskite Oxide Optimized for Oxygen Evolution Catalysis from Molecular Orbital Principles. Science 2011, 334, 1383–1385.

(15)

Over, H. Surface Chemistry of Ruthenium Dioxide in Heterogeneous Catalysis and

ACS Paragon Plus Environment

30

Page 31 of 43 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

ACS Sustainable Chemistry & Engineering

Electrocatalysis: From Fundamental to Applied Research. Chem. Rev. 2012, 112, 3356– 3426. (16)

Reier, T.; Teschner, D.; Lunkenbein, T.; Bergmann, A.; Selve, S.; Kraehnert, R.; Schlogl, R.; Strasser, P. Electrocatalytic Oxygen Evolution on Iridium Oxide: Uncovering CatalystSubstrate Interactions and Active Iridium Oxide Species. J. Electrochem. Soc. 2014, 161, 876–882.

(17)

McCrory, C. C. L.; Jung, S.; Peters, J. C.; Jaramillo, T. F. Benchmarking Heterogeneous Electrocatalysts for the Oxygen Evolution Reaction. J. Am. Chem. Soc. 2013, 135, 16977−16987.

(18)

Gao, X.; Zhang, H.; Li, Q.; Yu, X.; Hong, Z.; Zhang, X.; Liang, C.; Lin, Z. Hierarchical NiCo2O4 Hollow Microcuboids as Bifunctional Electrocatalysts for Overall WaterSplitting. Angew. Chem. Int. Ed. 2016, 55, 6290–6294.

(19)

Görlin, M.; Chernev, P.; De Araújo, J. F.; Reier, T.; Dresp, S.; Paul, B.; Krähnert, R.; Dau, H.; Strasser, P. Oxygen Evolution Reaction Dynamics, Faradaic Charge Efficiency, and the Active Metal Redox States of Ni-Fe Oxide Water Splitting Electrocatalysts. J. Am. Chem. Soc. 2016, 138, 5603−5614.

(20)

Chen, Q.; Switzer, J. A. Electrodeposition of Nanometer-Thick Epitaxial Films of Silver onto Single-Crystal Silicon Wafers. J. Mater. Chem. C 2019, 7, 1720-1725

(21)

Gong, M.; Li, Y.; Wang, H.; Liang, Y.; Wu, J. Z.; Zhou, J.; Wang, J.; Regier, T.; Wei, F.; Dai, H. An Advanced Ni-Fe Layered Double Hydroxide Electrocatalyst for Water Oxidation. J. Am. Chem. Soc. 2013, 135, 8452–8455.

ACS Paragon Plus Environment

31

ACS Sustainable Chemistry & Engineering 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

(22)

Page 32 of 43

Kuang, M.; Han, P.; Wang, Q.; Li, J.; Zheng, G. CuCo Hybrid Oxides as Bifunctional Electrocatalyst for Efficient Water Splitting. Adv. Funct. Mater. 2016, 26, 8555–8561.

(23)

Chemelewski, W. D.; Lee, H. C.; Lin, J. F.; Bard, A. J.; Mullins, C. B. Amorphous FeOOH Oxygen Evolution Reaction Catalyst for Photoelectrochemical Water Splitting. J. Am. Chem. Soc., 2014, 136, 2843–2850.

(24)

Deng, X.; Tüysüz, H. Cobalt-Oxide-Based Materials as Water Oxidation Catalyst: Recent Progress and Challenges. ACS Catal. 2014, 4, 3701–3714.

(25)

Anantharaj, S.; Ede, S. R.; Sakthikumar, K.; Karthick, K.; Mishra, S.; Kundu, S. Recent Trends and Perspectives in Electrochemical Water Splitting with an Emphasis on Sulfide, Selenide, and Phosphide Catalysts of Fe, Co, and Ni: A Review. ACS Catal. 2016, 6, 8069– 8097.

(26)

Swesi, A. T.; Masud, J.; Nath, M. Nickel Selenide as a High-Efficiency Catalyst for Oxygen Evolution Reaction. Energy Environ. Sci. 2016, 9, 1771–1782.

(27)

Tang, C.; Cheng, N.; Pu, Z.; Xing, W.; Sun, X. NiSe Nanowire Film Supported on Nickel Foam: An Efficient and Stable 3D Bifunctional Electrode for Full Water Splitting. Angew. Chem. Int. Ed. 2015, 54, 9351–9355.

(28)

Kwak, I. H.; Im, H. S.; Jang, D. M.; Kim, Y. W.; Park, K.; Lim, Y. R.; Cha, E. H.; Park, J. CoSe2 and NiSe2 Nanocrystals as Superior Bifunctional Catalysts for Electrochemical and Photoelectrochemical Water Splitting. ACS Appl. Mater. Interfaces 2016, 8, 5327–5334.

(29)

Li, W.; Gao, X.; Xiong, D.; Wei, F.; Song, W. G.; Xu, J.; Liu, L. Hydrothermal Synthesis of Monolithic Co3Se4 Nanowire Electrodes for Oxygen Evolution and Overall Water

ACS Paragon Plus Environment

32

Page 33 of 43 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

ACS Sustainable Chemistry & Engineering

Splitting with High Efficiency and Extraordinary Catalytic Stability. Adv. Energy Mater. 2017, 7, 1602579. (30)

De Silva, U.; Masud, J.; Zhang, N.; Hong, Y.; Liyanage, W. P. R.; Asle Zaeem, M.; Nath, M. Nickel Telluride as a Bifunctional Electrocatalyst for Efficient Water Splitting in Alkaline Medium. J. Mater. Chem. A 2018, 6, 7608–7622.

(31) Xu, K.; Ding, H.; Lv, H.; Tao, S.; Chen, P.; Wu, X.; Chu, W.; Wu, C.; Xie, Y. Understanding Structure-Dependent Catalytic Performance of Nickel Selenides for Electrochemical Water Oxidation. ACS Catal. 2017, 7, 310–315. (32)

Carlson, R. R.; Meek, D. W. Dipole Moments of Several Tertiary Phosphine Oxides, Sulfides, and Selenides and of Some Tertiary Arsine Oxides and Sulfides. Inorg. Chem. 1974, 13, 1741–1747.

(33)

Cao, X.; Hong, Y.; Zhang, N.; Chen, Q.; Masud, J.; Asle Zaeem, M.; Nath, M. Phase Exploration and Identification of Multinary Transition-Metal Selenides as High-Efficiency Oxygen Evolution Electrocatalysts through Combinatorial Electrodeposition. ACS Catal. 2018, 8, 8273–8289.

(34)

Swesi, A. T.; Masud, J.; Liyanage, W. P. R.; Umapathi, S.; Bohannan, E.; Medvedeva, J.; Nath, M. Textured NiSe2 Film: Bifunctional Electrocatalyst for Full Water Splitting at Remarkably Low Overpotential with High Energy Efficiency. Sci. Rep. 2017, 7, 2401.

(35)

Gao, M. R.; Cao, X.; Gao, Q.; Xu, Y. F.; Zheng, Y. R.; Jiang, J.; Yu, S. H. Nitrogen-Doped Graphene Supported CoSe2 Nanobelt Composite Catalyst for Efficient Water Oxidation. ACS Nano 2014, 8, 3970–3978.

ACS Paragon Plus Environment

33

ACS Sustainable Chemistry & Engineering 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

(36)

Page 34 of 43

Masud, J.; Swesi, A. T.; Liyanage, W. P. R.; Nath, M. Cobalt Selenide Nanostructures: An Efficient Bifunctional Catalyst with High Current Density at Low Coverage. ACS Appl. Mater. Interfaces 2016, 8, 17292–17302.

(37)

Liu, B.; Zhao, Y. F.; Peng, H. Q.; Zhang, Z. Y.; Sit, C. K.; Yuen, M. F.; Zhang, T. R.; Lee, C. S.; Zhang, W. J. Nickel–Cobalt Diselenide 3D Mesoporous Nanosheet Networks Supported on Ni Foam: An All-pH Highly Efficient Integrated Electrocatalyst for Hydrogen Evolution. Adv. Mater. 2017, 29, 1606521.

(38)

Xu, X.; Song, F.; Hu, X. A Nickel Iron Diselenide-Derived Efficient Oxygen-Evolution Catalyst. Nat. Commun. 2016, 7, 12324.

(39)

Umapathi, S.; Masud, J.; Swesi, A. T.; Nath, M. FeNi2Se4–Reduced Graphene Oxide Nanocomposite: Enhancing Bifunctional Electrocatalytic Activity for Oxygen Evolution and Reduction through Synergistic Effects. Adv. Sustainable Syst. 2017, 1, 1700086.

(40)

Akbar, K.; Jeon, J. H.; Kim, M.; Jeong, J.; Yi, Y.; Chun, S.-H. Bifunctional Electrodeposited 3D NiCoSe2/Nickle Foam Electrocatalysts for Its Applications in Enhanced Oxygen Evolution Reaction and for Hydrazine Oxidation. ACS Sustainable Chem. Eng. 2018, 6, 7735-7742.

(41)

Suen, N. T.; Hung, S. F.; Quan, Q.; Zhang, N.; Xu, Y. J.; Chen, H. M. Electrocatalysis for the Oxygen Evolution Reaction: Recent Development and Future Perspectives. Chem. Soc. Rev. 2017, 46, 337–365.

(42)

Maiyalagan, T.; Jarvis, K. A.; Therese, S.; Ferreira, P. J.; Manthiram, A. Spinel-Type Lithium Cobalt Oxide as a Bifunctional Electrocatalyst for the Oxygen Evolution and

ACS Paragon Plus Environment

34

Page 35 of 43 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

ACS Sustainable Chemistry & Engineering

Oxygen Reduction Reactions. Nat. Commun. 2014, 5, 3949. (43) Gong, M.; Dai, H. A Mini Review on NiFe-Based Materials as Highly Active Oxygen Evolution Reaction Electrocatalysts. Nano Res. 2015, 8, 23–39. (44)

Enman, L. J.; Burke, M. S.; Batchellor, A. S.; Boettcher, S. W. Effects of Intentionally Incorporated Metal Cations on the Oxygen Evolution Electrocatalytic Activity of Nickel (Oxy) Hydroxide in Alkaline Media. ACS Catal. 2016, 6, 2416–2423.

(45)

Görlin, M.; Ferreira de Araujo, J..; Schmies, H.; Bernsmeier, D.; Dresp, S.; Gliech, M.; Jusys, Z.; Chernev, P.; Kraehnert, R.; Dau, H.; Strasser, P. Tracking Catalyst Redox States and Reaction Dynamics in Ni-Fe Oxyhydroxide Oxygen Evolution Reaction Electrocatalysts: The Role of Catalyst Support and Electrolyte PH. J. Am. Chem. Soc. 2017, 139, 2070–2082.

(46)

Yang, Y.; Fei, H.; Ruan, G.; Xiang, C.; Tour, J. M. Efficient Electrocatalytic Oxygen Evolution on Amorphous Nickel-Cobalt Binary Oxide Nanoporous Layers. ACS Nano 2014, 8, 9518–9523.

(47)

Simon, G.; Kresler, S. E.; Essene, E. J. Phase Relations among Selenides, Sulfides, Tellurides, and Oxides: II. Applications to Selenide-Bearing Ore Deposits. Econ. Geol. 1997, 92, 468–484.

(48)

Xu, H.; Feng, J. X.; Tong, Y. X.; Li, G. R. Cu2O-Cu Hybrid Foams as High-Performance Electrocatalysts for Oxygen Evolution Reaction in Alkaline Media. ACS Catal. 2017, 7, 986–991.

(49)

Zhang, H.; Zhang, Z.; Li, N.; Yan, W.; Zhu, Z. Cu2O@C Core/Shell Nanoparticle as an

ACS Paragon Plus Environment

35

ACS Sustainable Chemistry & Engineering 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 36 of 43

Electrocatalyst for Oxygen Evolution Reaction. J. Catal. 2017, 352, 239–245. (50)

Masud, J.; Liyanage, W. P. R.; Cao, X.; Saxena, A.; Nath, M. Copper Selenides as HighEfficiency Electrocatalysts for Oxygen Evolution Reaction. ACS Appl. Energy Mater. 2018, 1, 4075–4083.

(51)

Cao, X.; Johnson, E.; Nath, M. Identifying High-efficiency Oxygen Evolution Electrocatalysts from Co-Ni-Cu Based Selenides through Combinatorial Electrodeposition. Submitted to J. Mater. Chem. A. 2019, doi: 10.1039/c9ta00863b.

(52)

Batchellor, A. S.; Kwon, G.; Laskowski, F. A. L.; Tiede, D. M.; Boettcher, S. W. Domain Structures of Ni and NiFe (Oxy)Hydroxide Oxygen-Evolution Catalysts from X-Ray Pair Distribution Function Analysis. J. Phys. Chem. C 2017, 121, 25421–25429.

(53)

Zhuang, L.; Ge, L; Yang, Y.; Li, M.; Jia, Y.; Yao, X.; Zhu, Z. Ultrathin Iron‐Cobalt Oxide Nanosheets with Abundant Oxygen Vacancies for the Oxygen Evolution Reaction. Adv. Mater. 2017, 29, 1606793.

(54)

Santos, E.; Schmickler, W. d-Band Catalysis in Electrochemistry. ChemPhysChem 2006, 7, 2282–2285.

(55)

Stevens, M. B.; Trang, C. D. M.; Enman, L. J.; Deng, J.; Boettcher, S. W. Reactive Fe-Sites in Ni/Fe (Oxy)Hydroxide Are Responsible for Exceptional Oxygen Electrocatalysis Activity. J. Am. Chem. Soc. 2017, 139, 11361−11364.

(56)

Xu, D.; Stevens, M. B.; Rui, Y.; DeLuca, G.; Boettcher, S. W.; Reichmanis, E.; Li, Y.; Zhang, Q.; Wang, H. The Role of Cr Doping in Ni–Fe Oxide/(Oxy)Hydroxide Electrocatalysts for Oxygen Evolution. Electrochim. Acta 2018, 265, 10–18.

ACS Paragon Plus Environment

36

Page 37 of 43 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

ACS Sustainable Chemistry & Engineering

(57)

Burke, M. S.; Kast, M. G.; Trotochaud, L.; Smith, A. M.; Boettcher, S. W. Cobalt-Iron (Oxy)Hydroxide Oxygen Evolution Electrocatalysts: The Role of Structure and Composition on Activity, Stability, and Mechanism. J. Am. Chem. Soc. 2015, 137, 3638−3648.

(58)

Xu, X.; Liang, H.; Ming, F.; Qi, Z.; Xie, Y.; Wang, Z. Prussian Blue Analogues Derived Penroseite (Ni,Co)Se2 Nanocages Anchored on 3D Graphene Aerogel for Efficient Water Splitting. ACS Catal. 2017, 7, 6394–6399.

(59)

Zhao, X.; Zhang, H.; Yan, Y.; Cao, J.; Li, X.; Zhou, S.; Peng, Z.; Zeng, J. Engineering the Electrical Conductivity of Lamellar Silver-Doped Cobalt (II) Selenide Nanobelts for Enhanced Oxygen Evolution. Angew. Chem. Int. Ed. 2017, 56, 328–332.

(60)

Trotochaud, L.; Young, S. L.; Ranney, J. K.; Boettcher, S. W. Nickel-Iron Oxyhydroxide Oxygen-Evolution Electrocatalysts: The Role of Intentional and Incidental Iron Incorporation. J. Am. Chem. Soc. 2014, 136, 6744–6753.

(61)

Qi, J.; Zhang, W.; Xiang, R.; Liu, K.; Wang, H. Y.; Chen, M.; Han, Y.; Cao, R. Porous Nickel–Iron Oxide as a Highly Efficient Electrocatalyst for Oxygen Evolution Reaction. Adv. Sci. 2015, 2, 1500199.

(62)

Wang, Z.; Li, J.; Tian, X.; Wang, X.; Yu, Y.; Owusu, K. A.; He, L.; Mai, L. Porous NickelIron Selenide Nanosheets as Highly Efficient Electrocatalysts for Oxygen Evolution Reaction. ACS Appl. Mater. Interfaces 2016, 8, 19386–19392.

(63)

Li, Y. F.; Selloni, A. Mechanism and Activity of Water Oxidation on Selected Surfaces of Pure and Fe-Doped NiOx. ACS Catal. 2014, 4, 1148–1153.

ACS Paragon Plus Environment

37

ACS Sustainable Chemistry & Engineering 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

(64)

Page 38 of 43

Landon, J.; Demeter, E.; Inoǧlu, N.; Keturakis, C.; Wachs, I. E.; Vasić, R.; Frenkel, A. I.; Kitchin, J. R. Spectroscopic Characterization of Mixed Fe-Ni Oxide Electrocatalysts for the Oxygen Evolution Reaction in Alkaline Electrolytes. ACS Catal. 2012, 2, 1793–1801.

(65)

Enman, L. J.; Stevens, M. B.; Dahan, M. H.; Nellist, M. R.; Toroker, M. C.; Boettcher, S. W. Operando X-Ray Absorption Spectroscopy Shows Iron Oxidation Is Concurrent with Oxygen Evolution in Cobalt–Iron (Oxy)Hydroxide Electrocatalysts. Angew. Chem. Int. Ed. 2018, 57, 12840–12844.

(66)

Zhang, X.; Zhang, X.; Xu, H.; Wu, Z.; Wang, H.; Liang, Y. Iron-Doped Cobalt Monophosphide Nanosheet/Carbon Nanotube Hybrids as Active and Stable Electrocatalysts for Water Splitting. Adv. Funct. Mater. 2017, 27, 1606635.

(67)

Zhang, J.; Wang, Y.; Zhang, C.; Gao, H.; Lv, L.; Han, L.; Zhang, Z. Self-Supported Porous NiSe2 Nanowrinkles as Efficient Bifunctional Electrocatalysts for Overall Water Splitting. ACS Sustain. Chem. Eng. 2018, 6, 2231–2239.

(68)

Anantharaj, S.; Kennedy, J.; Kundu, S. Microwave-Initiated Facile Formation of Ni3Se4 Nanoassemblies for Enhanced and Stable Water Splitting in Neutral and Alkaline Media. ACS Appl. Mater. Interfaces 2017, 9, 8714–8728.

(69)

Friebel, D.; Louie, M. W.; Bajdich, M.; Sanwald, K. E.; Cai, Y.; Wise, A. M.; Cheng, M. J.; Sokaras, D.; Weng, T. C.; Alonso-Mori, R.; Davis, R. C.; Bargar, J. R.; Nørskov, J.K.; Nilsson, A.; Bell, A. T. Identification of Highly Active Fe Sites in (Ni,Fe)OOH for Electrocatalytic Water Splitting. J. Am. Chem. Soc. 2015, 137, 1305–1313.

(70)

Shinagawa, T.; Garcia-Esparza, A. T.; Takanabe, K. Insight on Tafel Slopes from a

ACS Paragon Plus Environment

38

Page 39 of 43 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

ACS Sustainable Chemistry & Engineering

Microkinetic Analysis of Aqueous Electrocatalysis for Energy Conversion. Sci. Rep. 2015, 5, 13801. (71)

Chi, J. Q.; Yan, K. L.; Xiao, Z.; Dong, B.; Shang, X.; Gao, W. K.; Li, X.; Chai, Y. M.; Liu, C. G. Trimetallic Ni-Fe-Co Selenides Nanoparticles Supported on Carbon Fiber Cloth as Efficient Electrocatalyst for Oxygen Evolution Reaction. Int. J. Hydrogen Energy 2017, 42, 20599–20607.

(72)

Smith, R. D. L.; Prévot, M. S.; Fagan, R. D.; Trudel, S.; Berlinguette, C. P. Water Oxidation Catalysis: Electrocatalytic Response to Metal Stoichiometry in Amorphous Metal Oxide Films Containing Iron, Cobalt, and Nickel. J. Am. Chem. Soc. 2013, 135, 11580–11586.

(73)

Grimaud, A.; May, K. J.; Carlton, C. E.; Lee, Y. L.; Risch, M.; Hong, W. T.; Zhou, J.; ShaoHorn, Y. Double Perovskites as a Family of Highly Active Catalysts for Oxygen Evolution in Alkaline Solution. Nat. Commun. 2013, 4, 2439.

(74)

Shen, M.; Ruan, C.; Chen, Y.; Jiang, C.; Ai, K.; Lu, L. Covalent Entrapment of Cobalt-Iron Sulfides in N-Doped Mesoporous Carbon: Extraordinary Bifunctional Electrocatalysts for Oxygen Reduction and Evolution Reactions. ACS Appl. Mater. Interfaces 2015, 7, 1207– 1218.

(75)

Zhang, P.; Li, L.; Nordlund, D.; Chen, H.; Fan, L.; Zhang, B.; Sheng, X.; Daniel, Q.; Sun, L. Dendritic Core-Shell Nickel-Iron-Copper Metal/Metal Oxide Electrode for Efficient Electrocatalytic Water Oxidation. Nat. Commun. 2018, 9, 381.

(76)

Liu, Y.; Niu, Z.; Lu, Y.; Zhang, L.; Yan, K. Facile Synthesis of CuFe2O4 crystals Efficient for Water Oxidation and H2O2 reduction. J. Alloys Compd. 2018, 735, 654–659.

ACS Paragon Plus Environment

39

ACS Sustainable Chemistry & Engineering 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

(77)

Page 40 of 43

Chen, H.; Gao, Y.; Ye, L.; Yao, Y.; Chen, X.; Wei, Y.; Sun, L. A Cu2Se-Cu2O Film Electrodeposited on Titanium Foil as a Highly Active and Stable Electrocatalyst for the Oxygen Evolution Reaction. Chem. Commun. 2018, 54, 4979–4982.

(78)

Niu, Y.; Qian, X.; Xu, C.; Liu, H.; Wu, W.; Hou, L. Cu-Ni-CoSex Quaternary Porous Nanocubes as Enhanced Pt-Free Electrocatalysts for Highly Efficient Dye-Sensitized Solar Cells and Hydrogen Evolution in Alkaline Medium. Chem. Eng. J. 2019, 357, 11–20.

(79)

Ali, Z.; Asif, M.; Huang, X.; Tang, T.; Hou, Y. Hierarchically Porous Fe 2CoSe4 BinaryMetal Selenide for Extraordinary Rate Performance and Durable Anode of Sodium-Ion Batteries. Adv. Mater. 2018, 30, 1802745.

(80)

Zhao, Y.; Jin, B.; Zheng, Y.; Jin, H.; Jiao, Y.; Qiao, S. Z. Charge State Manipulation of Cobalt Selenide Catalyst for Overall Seawater Electrolysis. Adv. Energy Mater. 2018, 8, 1801926.

(81)

Yu, J.; Cheng, G.; Luo, W. 3D Mesoporous Rose-like Nickel-Iron Selenide Microspheres as Advanced Electrocatalysts for the Oxygen Evolution Reaction. Nano Res. 2018, 11, 2149–2158.

(82)

Quan, L.; Liu, T.; Yi, M.; Chen, Q.; Cai, D.; Zhan, H. Construction of Hierarchical Nickel Cobalt Selenide Complex Hollow Spheres for Pseudocapacitors with Enhanced Performance. Electrochim. Acta 2018, 281, 109–116.

(83)

Jia, J.; Seitz, L. C.; Benck, J. D.; Huo, Y.; Chen, Y.; Ng, J. W. D.; Bilir, T.; Harris, J. S.; Jaramillo, T. F. Solar Water Splitting by Photovoltaic-Electrolysis with a Solar-toHydrogen Efficiency over 30 %. Nat. Commun. 2016, 7, 13237.

ACS Paragon Plus Environment

40

Page 41 of 43 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

ACS Sustainable Chemistry & Engineering

(84)

Xiao, P.; Sk, M. A.; Thia, L.; Ge, X.; Lim, R. J.; Wang, J. Y.; Lim, K. H.; Wang, X. Molybdenum Phosphide as an Efficient Electrocatalyst for the Hydrogen Evolution Reaction. Energy Environ. Sci. 2014, 7, 2624–2629.

(85)

Merki, D.; Vrubel, H.; Rovelli, L.; Fierro, S.; Hu, X. Fe, Co, and Ni Ions Promote the Catalytic Activity of Amorphous Molybdenum Sulfide Films for Hydrogen Evolution. Chem. Sci. 2012, 3, 2515–2525.

(86)

Sivanantham, A.; Ganesan, P.; Shanmugam, S. Hierarchical NiCo2S4 Nanowire Arrays Supported on Ni Foam: An Efficient and Durable Bifunctional Electrocatalyst for Oxygen and Hydrogen Evolution Reactions. Adv. Funct. Mater. 2016, 26, 4661–4672.

(87)

Chen, P.; Xu, K.; Zhou, T.; Tong, Y.; Wu, J.; Cheng, H.; Lu, X.; Ding, H.; Wu, C.; Xie, Y. Strong-Coupled Cobalt Borate Nanosheets/Graphene Hybrid as Electrocatalyst for Water Oxidation under Both Alkaline and Neutral Conditions. Angew. Chem. Int. Ed. 2016, 55, 2488–2492.

(88)

Jia, X.; Zhao, Y.; Chen, G.; Shang, L.; Shi, R.; Kang, X.; Waterhouse, G. I. N.; Wu, L. Z.; Tung, C. H.; Zhang, T. Ni3FeN Nanoparticles Derived from Ultrathin NiFe-Layered Double Hydroxide Nanosheets: An Efficient Overall Water Splitting Electrocatalyst. Adv. Energy Mater. 2016, 6, 1502585.

(89)

Schwanke, C.; Stein, H. S.; Xi, L.; Sliozberg, K.; Schuhmann, W.; Ludwig, A.; Lange, K. M. Correlating Oxygen Evolution Catalysts Activity and Electronic Structure by a HighThroughput Investigation of Ni1-y-zFeyCrzOx. Sci. Rep. 2017, 7, 44192.

(90)

Yan, X.; Tian, L.; Li, K.; Atkins, S.; Zhao, H.; Murowchick, J.; Liu, L.; Chen, X.

ACS Paragon Plus Environment

41

ACS Sustainable Chemistry & Engineering 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 42 of 43

FeNi3/NiFeOx Nanohybrids as Highly Efficient Bifunctional Electrocatalysts for Overall Water Splitting. Adv. Mater. Interfaces 2016, 3, 1600368. (91)

Yan, X.; Tian, L.; Atkins, S.; Liu, Y.; Murowchick, J.; Chen, X. Converting CoMoO4 into CoO/MoOx for Overall Water Splitting by Hydrogenation. ACS Sustain. Chem. Eng. 2016, 4, 3743–3749.

(92)

Liu, L.; Niu, D.; Song, F.; Lyu, L.; He, J.; Yan, X.; Li, K.; Hu, X.; Chen, X. From Water Oxidation to Reduction: Transformation from NixCo3– xO4 Nanowires to NiCo/NiCoOx Heterostructures. ACS Appl. Mater. Interfaces 2016, 8, 3208–3214.

ACS Paragon Plus Environment

42

Page 43 of 43 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

ACS Sustainable Chemistry & Engineering

TOC Graphic Combinatorial electrodeposition in Fe-Co-Cu phase-space has led to discovery of high-efficiency catalyst, (Fe0.48Co0.38Cu0.14)Se showing overpotential of 256 mV at 10mA.cm-2.

η (V / at 10 mA cm-2)

ACS Paragon Plus Environment

43