Scholarly article on topic 'Enhanced hydrogen evolution reaction on hybrids of cobalt phosphide and molybdenum phosphide'

Enhanced hydrogen evolution reaction on hybrids of cobalt phosphide and molybdenum phosphide Academic research paper on "Nano-technology"

Share paper
Academic journal
R. Soc. open sci.
OECD Field of science

Academic research paper on topic "Enhanced hydrogen evolution reaction on hybrids of cobalt phosphide and molybdenum phosphide"



VJ Check for ^^ updates

Cite this article: Fang S-L, Chou T-C, Samireddi S, Chen K-H, Chen L-C, Chen W-F. 2017 Enhanced hydrogen evolution reaction on hybrids of cobalt phosphide and molybdenum phosphide. R. Soc. openscl. 4:161016.

Received: 8 December 2016 Accepted: 30 January 2017

Subject Category:


Subject Areas:

materials science/inorganic chemistry Keywords:

hydrogen evolution reaction, electrocatalyst, ternary, cobalt phosphide, molybdenum phosphide

Author for correspondence:

Wei-Fu Chen


This article has been edited by the Royal Society of Chemistry, including the commissioning, peer review process and editorial aspects up to the point of acceptance.

Electronic supplementary material is available online at figshare.c.3691936.

ROY %t




Enhanced hydrogen evolution reaction on hybrids of cobalt phosphide and molybdenum phosphide

Si-Ling Fang1, Tsu-Chin Chou1, Satyanarayana Samireddi2,3, Kuei-Hsien Chen1,2, Li-Chyong Chen1 and Wei-Fu Chen1

1Center for Condensed Matter Sciences, National Taiwan University, Taipei 10617, Taiwan, Republic of China

institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan, Republic of China

department of Chemistry, National Tsing Hua University, Hsinchu 30012, Taiwan, Republic of China

<2> W-FC, 0000-0002-3009-0196

Production of hydrogen from water electrolysis has stimulated the search of sustainable electrocatalysts as possible alternatives. Recently, cobalt phosphide (CoP) and molybdenum phosphide (MoP) received great attention owing to their superior catalytic activity and stability towards the hydrogen evolution reaction (HER) which rivals platinum catalysts. In this study, we synthesize and study a series of catalysts based on hybrids of CoP and MoP with different Co/Mo ratio. The HER activity shows a volcano shape and reaches a maximum for Co/Mo = 1. Tafel analysis indicates a change in the dominating step of Volmer-Hyrovsky mechanism. Interestingly, X-ray diffraction patterns confirmed a major ternary interstitial hexagonal CoMoP2 crystal phase is formed which enhances the electrochemical activity.

1. Introduction

Generation of hydrogen fuel from water as alternative to fossil oils without releasing carbonaceous gases, such as carbon monoxide and carbon dioxide [1], has been considered as a promising green technology. Electrolytic hydrogen evolution reaction (HER) has received great attention, because it can be activated from renewable sources of energy like wind and solar

2017 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution License, which permits unrestricted use, provided the original author and source are credited.

[2,3]. However, water electrolysers are hampered by high costs and limited abundance of electrode materials owing to the use of noble metals like platinum. This makes water electrolysers unfavourable to compete with natural gas reforming. It is of high priority to research earth-abundant catalysts as possible alternatives to noble metals [4-6]. Unfortunately, transition metal-based electrocatalysts have suffered from high overpotentials and corrosion problems in acidic media [7].

Transition metal-based pnictides [8,9] (nitrogen-group elements) are potential electrocatalysts for the HER, as they possess excellent corrosion resistance in the HER condition and good electrical conductivity as an electrode. Very recently, metal phosphides have received great attention owing to their superior catalytic activity towards the HER which rivals platinum catalysts [10-12]. It has previously been reported by Schaak and co-workers [13-15] that orthorhombic CoP and Co2P nanocrystallines produced cathodic current density of 20mAcm-2 at overpotentials ranging from 85 to 117mV and stable over 24 h of operation. Bulk hexagonal molybdenum phosphide (MoP) has shown stable HER activity in acidic solution with overpotentials of 140-246mV for driving a current density of 10mAcm-2 [16,17]. Even lower overpotentials were reported on amorphous or nanosized MoP HER when compared with bulk MoP [18,19]. It is established that doping Co to MoP increases the intrinsic activity of MoP (so-called Co-promoted MoP); however, the positioning of Co and the exact structural polytype are still not resolved [20]. A decent comprehension of the correspondence between the HER activity and the crystal structure of these materials is not yet thoroughly developed. For bimetal alloys, material design methodologies for optimizing the electronic structure and electrochemical activity are well investigated [21-25]; however, only a few studies were focused on the modification of structure of mixed metal pnictides to accelerate the electrolytic reaction. Cao et a/. [26] reported a layered Co0.6Mo1.4N2 allows the Co to tune the electronic states of molybdenum at the catalyst surface without disrupting the catalytic activity. Staszak-Jirkovsky et a/. [27] demonstrated Con+ cations in a CoMoS^ chalcogel structure helped accelerate the rate-determining Volmer step. Ternary Chevrel-phase NiMo3S4 was reported by Jiang et a/. [28] that the interconnected [Mo6S8]2 cluster units allow faster charge transfer.

Herein, we show that this research improves upon previous Co-promoted MoP catalysts towards the HER by optimizing the CoP and MoP compositions. We have used a solid-state synthesis route to prepare hybrids of CoP and MoP. Structural studies indicate that at the ratio of Co/Mo = 1, a new crystalline hexagonal CoMoP2 phase is formed in addition to CoP and MoP phases. The hybrid with Co/Mo = 1 exhibits good HER performance and shows promise as an efficient cost-effective cathode material for water splitting.

2. Material and methods

2.1. Synthesis of hybrids of CoP and MoP

Hybrids of cobalt phosphide and molybdenum phosphide were prepared by a two-step thermal treatment of mixtures of cobalt nitrate (Co(NO3)2.6H2O, Acros), ammonium molybdate ((NH4)6Mo7 O24.4H2O, Aldrich) and ammonium dihydrogen phosphate (NH4H2PO4, Aldrich). A typical procedure is shown as follows. For preparing Co0.5Mo0.5P catalyst, 10pmol ammonium molybdate, 70pmol cobalt nitrate and 140pmol ammonium dihydrogen phosphate were mixed in 50 ml water and were ultrasonicated until all salts were dissolved. The ratio of cobalt nitrate to ammonium molybdate was varied to create molar ratios [Co]/[Mo] of 0.1/0.9,0.3/0/7,0.5/0.5,0.7/0.3 and 0.9/0.1. After mixing, the transparent solution was dried at 80°C in an oven. The solid mixture was then annealed in a quartz tube furnace with a 100 sccm Ar flow from ambient to 350°C at a rate of 10°C min-1 and then held at 300°C for 1 h. Then the gas flow was switched to H2/ Ar mix flow (H2 50 sccm; Ar 50 sccm). The temperature was increased to 800°C at a ramping rate of 15°Cmin-1 and held at 800°C for 2h. Bulk CoP and MoP were prepared by the same procedure using corresponding precursors with [Co or Mo]/[P] = 1. To study the effect of annealing temperature, the precursors of Co°.5Mo0 5P were treated at 650°C and 1000°C under H2/Ar mix flow for 2 h after the annealing at 350°C under Ar.

2.2. Structural characterization

The micromorphology of hybrid catalysts were observed on a field emission scanning electron microscope (JOEL JSM-6700F) and on a JEOL JEM-2100 transmission electron microscope. The crystalline compositions of the Co^ MoTP hybrids were verified by X-ray diffraction (XRD) patterns with a Bruker D2 PHASER diffractometer. X-ray photoelectron spectra (XPS) were collected on a VG Scientific ESCALAB 250 and the binding energy was calibrated using the C 1s peak at 284.6 eV.

2.3. Electrochemical measurements

Electrochemical measurements were performed in a three-electrode electrochemical cell using a Zahner ZENNIUM E electrochemical workstation. The electrodes for the electrochemical measurements were fabricated as follows. Catalyst ink was prepared by mixing catalyst powder with milli-Q water solution (1ml for 10 mg of electrocatalyst). In total, 5% Nafion dispersion (DuPont) was added (50 mg solid Nafion for 100 mg of catalyst) to the catalyst slurry. Catalyst coating on glassy carbon electrode with 0.196 cm2 active area was fabricated by drop-casting the catalyst ink on with a micropipette. The catalyst loading was 0.4 mg catalyst cm-2. The electrode was then dried under vacuum at room temperature. The electrolytes for electrochemical measurements were prepared with perchloric acid (Aldrich) and Milli-Q water (Millipore). Pt foil (purity 99.999%) was purchased from Aldrich. All the electrochemical measurements were performed in 0.1 M HClO4 (aq) electrolyte, which was deaerated with hydrogen gas before use. A platinum wire, and an Ag/AgCl reference electrode (3M KCl) were used as the counter electrode and reference electrode, respectively. All potentials, E, are quoted with respect to reversible hydrogen electrode (RHE). The electrocatalysis was studied using linear sweeping voltage in the range of +0.2 V to -0.5 V (versus RHE). Electrochemical impedance spectroscopy (EIS) was performed in potentiostatic mode at an applied overpotential of 100 mV with frequency from 10 mHz to 0.1 MHz and amplitude of 5 mV. The accelerated deterioration test of electrocatalysts was examined by applying cyclic voltammetry in the potential range from +0.1 to -0.5 V versus RHE for 1000 cycles.

3. Results and discussion

The goal of this investigation is combining Mo with Co in order to tune the properties of their phosphides. Early attempts found enhancement of the HER activity on MoP by the doping of small amounts of Co onto its hexagonal structure [20]. The synthesis was reported to produce hexagonal MoP, while the state of cobalt was not known. In this study, a set of catalysts was prepared with different cobalt-to-all-metal atomic ratios, y (y = [Co]/[Co + Mo], t = [Mo]/[Co + Mo]). Co^MoTP (e.g. Co01Mo0 9P) is defined as the sample code of the obtained catalysts; it does not represent any crystalline phase. The Co^MoTP catalysts investigated were Co01Mo0 9P, Co0 3Mo0 7P, Co05Mo05P, Co0 7Mo0 3P and Co0 9Mo01P as shown in scheme 1.

We investigated the HER activities of the Co9MoTP catalysts, MoP (the same process for making Co^MoTP) and CoP in 0.1 M HQO4 solution using a typical three-electrode electrochemical cell. The HER activity is compared using the overpotential (n) at 10 mA cm-2 of cathodic current density (^10) and the current density at n = 200 mV. For assessing the electrochemical analysis, in figure 1a, a Pt foil exhibits HER activity with n10 of 50 mV which is comparable to other studies. The polarization curve recorded with Co9MoTP showed good activity for the HER. It is observed clearly in the electrochemical cell that the hydrogen bubbled more vigorously with increasing overpotentials. The n10 of the MoP and the CoP catalysts were read at 250 mV and 283 mV, respectively, which are similar to the overpotentials reported previously [29]. For Co^MoTP catalysts, the polarization curves for the HER in figure 1a changed with the composition. In figure 1b, the overpotential, n10, is plotted versus the cobalt content, y. The overpotential data of bulk MoS2 are collected from the literature [30]. It is observed that both bulk CoP and MoP show better HER activity than bulk MoS2. MoS2 is known as an edge-active material owing to its unique two-dimensional structure [31-34], while CoP and MoP are active in most facets. In figure 1b, it is shown that on adding Co to the MoP crystal (y = 0 — 0.5), the n10 decreases (better activity) from 250 to 165 mV, and then it forms a volcano shape with maxima at y = 0.5 as y increases to 1.0. This result indicates that when mixing CoP and MoP, the HER activity is enhanced. The Co0 5Mo0 5P catalyst showed an n10 at 165 mV which rivals other bulk catalysts reported. Semimetallic MoP2 nanoparticles showed an n10 of 143mV [35]. On CoS2 thin film, an n10 of 192mV was observed by Faber et al. [36]. Very recently, Jiang et al. [28] demonstrated that Chevrel-phase bimetal sulfide NiMo3S4 possesses good activity with an n10 of 257mV. Thus, the Co05Mo05P catalyst presented here is one of the top non-Pt HER catalysts in the bulk form. In figure 1c, the current density at n = 200 mV of the Co9MoTP catalysts showed a volcano shape as a function of the Co content. The Co0.5Mo0.5P catalyst showed an activity of 34.3 mA cm-2 which is 5.9 and 16.3 times the activity of MoP and CoP catalysts, respectively.

To understand how the composition affects the HER activity, XRD patterns were collected from these catalysts. The powder XRD studies in figure 2 evidence the presence of cobalt phosphide and molybdenum phosphide, and their relative peak intensities differ with the various ratios of Co/Mo. For y from 0 to 0.3, the patterns suggest that the major crystal phase is WC-type MoP (ICSD No: 186 874), with a = 3.25 Á, b = 3.25 Á, c = 3.24 Á, and space group symmetry P6m2 (#187). For y = 0.1,

Co(NO3)2 •«^fc*»®


mix in wat and dr'


Cno.')Mflo.ip Co"7MoIL'P

( <i" fMit"^l> 350°C, air



Mo:P= 1:1

800°C, H2/Ar

«x » •x •

Scheme 1. Schematic depictionsof the synthesis of Co^ MoTP through treating a mixtureof metal precursorsand ammonium hydrogen phosphate with a two-steps solid-state process. The difference in colour of these precursor mixtures with different Co: Mo: P ratios is shown as the inset.

new diffraction peaks appeared at 20 = 31.3, 39.7, 42.0, and 45.2° (marked as asterisk) with a very low intensity. The intensity of these peaks increases as p increases and reaches maximum at p = 0.5. These peaks were confirmed as hexagonal CoMoP2 structure (hereinafter h-CoMoP2 is used to represent this crystal phase in order to distinguish the difference from the hybrid composition Co0.5Mo0.5P) with a = 3.33 A, c = 11.22 A, and a space group P63/mmc (ICSD No: 624219, see the electronic supplementary material, figure S1). In addition to MoP and h-CoMoP2 phases at p = 0.5, minor signals related to MnP-type orthorhombic CoP (ICSD No. 43249) also exist. For spectra with p = 0.5-1.0, the intensity of CoP peaks increase, but h-CoMoP2 peaks decrease and MoP peaks disappear completely. The XRD pattern of the sample with p = 0.9 is not shown here because it showed a similar pattern to CoP. The formation of ternary h-CoMoP2 phase was previously explained by the composition of trigonal-prismatic MoP6 and octahedral CoP6 prism to form intermetallic linear -Mo-Co-Mo- chains along the c-axis by Guerin & Sergent [37]. Interestingly, the h-CoMoP2 crystal possesses a similar trigonal-prismatic MoP6 prism to MoP crystal. X-ray photoelectron spectroscopy (XPS) was used to probe the surface electronic properties in the Co0.5Mo0.5P hybrid. The Co 2p spectrum (see the electronic supplementary material, figure S2a) showed two dominating peaks appeared at 777.0 and 791.8 eV corresponding to the Co 2p3/2 and Co 2p1/2 signals, respectively. The minor broad peak at 780 eV indicates the presence of surface oxidized Co species resulting from the contact with air. The Mo 3d spectrum, as shown in electronic supplementary material, figure S2b, showed a strong 3d5/2 peak at 227.8 eV and 3d3/2 peak at 231.1 eV. These peaks are assigned to zero valence state metallic Mo which indicates no or very little oxidation occurs on the Mo site. In the P 2p XPS spectrum shown in the electronic supplementary material, figure S2c, the strong peak at 129.4eV is assigned to negative charged (metal-P6-) phosphide. The minor peak at 133.4eV is referred to the P-O species. The crystal field model describes that the strong P ligands split the Co 3d states in the octahedral CoP6 prism into high-spin t2g5eg2 ground state. Kibsgaard et a/. [11] have shown that smaller differential hydrogen adsorption-free energies AGh of the Co-bridge phosphorus site on CoP surface at a low coverage than that of the P site on MoP surface at a high coverage makes the CoP6 prism a high turnover-frequency towards the HER.

The reaction between the cobalt and/or molybdenum precursors with ammonium hydrogen phosphate under hydrogen environment at 800°C generated uniform phosphide microparticles. The SEM image in figure 3 shows that the CoP and the Co0.5Mo0.5P samples are composed of primary particles, but the MoP contains mostly secondary aggregated particles. These results indicate that the ratio of Co and Mo precursors affects the basic morphology of the catalysts. This can be ascribed to the complicated decomposition process regarding the reduction of MoOx in the high-temperature reduction reaction. The Co0.5Mo0.5P hybrid was characterized by transmission electron microscopy (TEM).

¿T I -5

£ -15

si -20

n e -a -25

nt e r -30

r u u -35

-40 -0

(mV) 50

CM 'S 100

ial 250

nt te ot 300

rp er 350

..jT"/ • .üf

vv ■ ▼*

■ T*

: t* ■

Co0'9Mo0'1P Co0'7Mo0'3P Co05Mo05P Co0'3Mo0'7P Co0'1Mo0'9P Pt

-0.3 -0.2 -0.1 0 0.1 potential (V versus RHE)

0.4 0.6 Co/(Co + Mo), j

< 25 S

^ 20 si

10 5 0

0 0.1 0.3 0.5 0.7 Co/(Co + Mo), j

0.9 1.0

Figure 1. (a) Polarization curves of the Co ^ MoTP catalysts annealed at 800°Cand aPtfoil inhydrogen-purged 0.1 M HClO4.(6) Activity volcano for the HER showing the overpotentials from (a)at current density of10 mA cm-2 asafunctionofcobaltcontent ^.The data of Pt and MoS2 from the literature are also included. (c) The dependence of current density at an overpotential of 200 mV from (a) on the cobalt content.

The TEM image (see the electronic supplementary material, figure S3) showed interconnected structure composed of large particles (0.3-0.6 pm) which is in good agreement with the SEM result. The energy dispersive X-ray (EDX) spectra (as shown in electronic supplementary material, figure S4) collected on the TEM showed signals from Co, Mo and P elements. The Cu signals resulted from the Cu grid. The atomic percentages of Co, Mo and P elements obtained are 26.9, 24.2 and 48.9%, respectively, which answers to the composition of Co0.5Mo0.5P well.

The effect of annealing temperature on the crystal structure and the HER activity of the Co0.5Mo0.5P catalyst are studied as shown in figure 4. The Co0 5Mo0 5P catalyst annealed at 650oC (Co0 5Mo05P-650, pink line in figure 4a) presents diffraction signals composed of h-CoMoP2 (black asterisk) and CoP (red hollow square) crystals. After annealed at 1000°C (Co05Mo05P-1000, black line in figure 4a),

Figure 2. X-ray diffraction patterns corresponding to the Co^MoTP catalysts annealed at 800°Cas well as the single-phase MoP and CoP powder. The diffraction peaks are assigned to specific crystalline phases as follows. Black hollow circle, hexagonal MoP; red hollow square, orthorhombic CoP; black asterisk, hexagonal CoMoP2.

Figure 3. SEM images of (a) CoP, (b) Co05 Mo05 Pa nd (c) MoP crystals taken at a magnification of 100 000x


H,AAA a a.


* Co0-5Mo0-5P-650

-0.4 -0.3 -0.2 potential (V versus

-0.1 RHE)

Figure 4. (a) X-ray diffraction patterns and (b) HER polarization curves of the Co05 Mo05 P samples annealed at 650°C, 800°C and 1000°C. The assignments of the diffraction peaks are included: black hollow circle, hexagonal MoP; red hollow square, orthorhombic CoP; black asterisk, hexagonal CoMoP2;green hollowtriangle, orthorhombic CoMoP. The HER polarization curves were collected in hydrogen-purged 0.1 M HClO4 solution.

a new orthorhombic CoMoP with a Pnma space group (ICSD No. 2421) is formed (marked as green hollow triangle). The result indicates that increasing annealing temperature to 1000°C induced a phase transformation from h-CoMoP2 to CoMoP, and the excess phosphorus was evaporated. The polarization curves for the HER of Coa5Moa5P-650, Coa5Moa5P-800 and Co°.5Mo°.5P-1000 are compared in figure 4b.

§ -0.15

v -0.20

| -0.25

■ Co01Mo09P

• Coa3Moa7P

★ Co05Mo05P ▲ Coa7Moa3P □ Co09Mo01P

-0.2 0 0.2 0.4 0.6 0.8 1.0 log (current density (mA cm-2))

120 100 80

a 40 H

Volmer step dominated

□ B*

Heyrovsky step dominated Pt Tafel step determined

0 0.2 0.4 0.6 0.8 1.0 Co/(Co + Mo), j

Figure 5. (a) Tafel plots ofthe CotpMoTP catalysts derived from figure 1a. (b) The dependence of average Tafel slopes on the cobalt content, y. The Tafel slope of a Pt electrode is also included.

The HER activities of these catalysts are in series of Co0 5Mo0 5P-1000 < Coa5Moa5P-650 < Coa5Moa5P-800. The Co°.5Mo°.5P-1000 catalyst showed a poor HER activity with nW of 338 mV, which indicates that CoMoP structure is less active than the P-rich h-CoMoP2. This can be explained by the higher phosphorus content in the h-CoMoP2 (28.6%) than that in CoMoP (16.7%). Callejas et al. [ 15] reported that CoP showed significantly lower overpotential than Co2P to produce the same current density. Xiao et al. [ L7] compared the HER activity between bulk MoP and Mo3P and concluded that MoP is more active than Mo3 P. Thus, in this study, the h-CoMoP2 with higher P content showed better HER activity.

The dependence of the gradient of logarithm of current density versus potential (so-called Tafel plot) corresponds to the reaction mechanism of the HER. The Tafel curves recorded on the Co9MoTP catalysts (annealed at 800°C) as shown in figure 5a proceed in accordance with the classical two-electron-reaction model. The cathodic hydrogen evolution in acidic conditions involves either Volmer-Tafel (proton discharge followed by the recombination of two bound hydrogen) or Volmer-Heyrovsky mechanism (proton discharge followed by electrochemical hydrogen desorption). Tafel slopes were obtained by linear fitting of the low current density region (-0.1 to +0.5log[mAcm-2]). Figure 5b shows the dependence of the obtained Tafel slopes on the cobalt content, y. For Pt foil, the Tafel slope of 30.1 mVdec-1 suggests a Volmer-Tafel mechanism that the recombination of two Hads is rate-determining. The data in figure 5b showed decreasing Tafel slopes from 84.0 to 60.5 mV dec-1 for samples from y = 0 to 0.5. These results suggest that at a small y (y = 0, high Tafel slope) hydrogen evolution occurs via a Volmer-Heyrovsky mechanism in which slow adsorption of proton dominates the kinetics. For y = 0.5, the small Tafel slope of 60.5mVdec-1 indicates that the electrochemical desorption of hydrogen influences the reaction kinetics through the Volmer-Heyrovsky mechanism much more than the case of MoP catalyst. Comparing the electrochemical properties of Co0.5Mo0.5P catalyst to bulk CoP and MoP catalysts, the small overpotential and the low Tafel slope of Co^Mo^P demonstrate that the presence of h-CoMoP2 crystal phase in the catalyst favours proton adsorption kinetics. To further ensure that h-CoMoP2 crystal phase enhances HER activity, we collected polarization curves from three 50/50 mixtures of CoP/MoP powder. The obtained polarization curves of the physical mixtures (see the electronic supplementary material, figure S5) showed HER activity ranging between pure CoP and MoP catalysts. This result indicates clearly that the physical mixture of CoP and MoP powder does not enhance the HER activity, while the h-CoMoP2 crystal phase does.

We carried out EIS analyses on these catalysts to study the charge transfer in the catalytic turnover. The measurements were conducted at overvoltage of 100 mV. The Nyquist plots of the Co^Mo^P, CoP and MoP samples as shown in figure 6 revealed classical two time-constants characteristics that correspond to a combination of kinetic and diffusion response of the HER on rough electrode surfaces. This two time-constants equivalent circuit is included in the electronic supplementary material, figure S6. The series resistance, Rs, is a sum of all ohmic resistances that the electron passed through from wires to the catalyst surface. The double-layer capacitance, C^i, in parallel with the charge transfer resistance, Rct, contributes to the charge-transfer process; the capacitance of the catalyst coating, C^2, in parallel with the resistance of ion conducting paths that develops in the porous surface, Rp. The curves in figure 6 were fitted with the above two time-constant equivalent circuit. The fitting agrees well with the experimental results (see the electronic supplementary material, figure S7). These electrodes showed small series resistances between 2.5 and 4 O which indicates good adhesion of the powder on glassy

r 30 ^

20 10 0

★ ★ ★

Coa5Moa5P £

0 20 40 60 80 100 120 140 160 Z (Ohm)

Figure 6. Comparison of the Nyquist plotsof the Co05 Mo05 P, CoP and MoP catalysts at an overpotential of100 mV in 0.1 M HClO4 solution.

-30 1 -0.3

gdCI 0*0«D CD

o" a" 0»

o* d* 0» Q»

• initial

o after 1000 cycles

-0.2 -0.1 0 potential (V versus RHE)

time (min)

100 90 £ 80 * 70 g 60 -g 50 £

40 ic e 30 | 20 | 10 p-0

Figure7. (a) HER polarization curves of Co05Mo0 5 P electrode before and after potential sweeps (-0.5 ~ +0.1 V versus RHE) for 1000 cycles in 0.1 M HClO4 solution (scan rate 2 mV s-1) (b). Time dependence of the quantity of hydrogen produced experimentally (dashed line) and the Faradaic efficiency of the Co05Mo05P catalyst at a fixed cathodic current density of 10 mA cm-2 in 0.1 M HClO4 solution 60 min.

carbon electrode. In figure 6, the pronounced semicircle at low frequencies (high Z') returns estimates of the charge-transfer resistance, Rct. The Rct of the Co^Mo^P catalyst (92.0 O) is found much lower than CoP (154.1 O) and MoP catalysts (107.5 O). As reported recently, the Rct of bulk MoS2 was found to be 150 O at n = 150 mV [30]. MoS2/RGO has an Rct of 250 O at n = 120 mV [38]; Mo2C nanowires showed a low Rct of 90 O at n = 150 mV [39]. Thus, the present Co^Mo^P catalyst is regarded as a catalyst with efficient charge transfer property. Efficient charge transfer kinetics reflects to the accelerating of the proton discharge step in the Volmer-Heyrovsky mechanism and thus lowers its Tafel slope as aforementioned. The enhancement in Rct can be ascribed to the modification of d-band structure owing to the formation of Co-P-Mo linkage.

To assess the stability of the Co0.5Mo0.5P electrodes, we performed the accelerated deterioration experiment by sweeping the applied potential from -0.5 to + 0.2 V versus RHE in 0.1 M HQO4 solution. After 1000 cycles on the Co^Mo^P electrode, only a slight shift (11 mV increase in n10) in the polarization curve is observed as shown in figure 7a, indicating a good stability of the Co0.5Mo0.5P catalyst in the operating conditions. Faradaic efficiency was determined by performing a controlled potential electrolysis experiment using an H-type cell. The amount of hydrogen gas produced was calculated from the volume of gas at a fixed cathodic current density of 10 mA cm-2 in 0.1 M HClO4 solution for 60 min. The experimentally determined hydrogen quantity was compared to the calculated amount based on the charge consumed as shown in figure 7b. The Faradaic efficiency reached 98% after 60 min of operation. The cathodic current density recorded for driving the HER at overpotential of 165 mV as a function of time is plotted in the electronic supplementary material, figure S8. The current density slightly decreased about 5% after 6 h of electrolysis. These results imply that the Co0.5Mo0.5P catalyst is a highly efficient and cathodically stable HER catalyst in the acidic environment.

4. Conclusion

We have studied the relationship between the composition of Co^MoTP hybrids and their activity in the HER. The HER activity of the Co9MoTP hybrids showed a volcano shape with a maxima at Co content of 0.5. The Co0 5Mo0 5P catalyst possesses a smaller overpotential (165 mV for driving 10 mA cm-2 of current density) when compared with CoP and MoP and is comparable to bulk non-precious catalysts reported. The hexagonal CoMoP2 phase was found to be responsible for the HER activity in the Co05Mo05P catalyst. The Tafel analysis showed a change in the dominating step in Volmer-Heyrovsky mechanism. The low Tafel slope of the Co0 5Mo0 5P catalyst suggests that the electrochemical desorption of adsorbed hydrogen is the rate-determining step in the HER. The Co0 5Mo0 5P catalyst showed low charge transfer resistance, high stability and high Faradaic efficiency; all demonstrate that the Co0 5Mo0 5P catalyst is a potential, high-performance electrocatalyst for water electrolysis in acidic environments.

Data accessibility. Electronic supplementary material including models and fitting of the EIS, XRD, XPS, EDX spectra and polarization curves.

Authors' contributions. W.-F.C. conceived and designed the experiments. S.-L.F. carried out the synthesis and characterization. S.S. contributed to the SEM analysis and the discussion. L.-C.C. and K.-H.C. contributed to the discussion of the results. W.-F.C. and T.-C.C. contributed to the analysis of the data and the writing of the manuscript. Competing interests. The authors declare that they have no competing interests.

Funding. This work was supported by Ministry of Science and Technology (104-2113-M-002-013-MY2). S.-L.F., T.-C.C. and L.-C.C. are grateful for the support by Academia Summit Program from Ministry of Science and Technology (104-2745-M-002-004-ASP; 105-2745-M-002-003-ASP).


1. Vesborg PCK, Seger B, Chorkendorff 1.2015 Recent development in hydrogen evolution reaction catalysts and their practical implementation.

J. Phys. Chem. Lett. 6,951-957. (doi:10.1021/ acs.jpclett.5b00306)

2. May MM, Lewerenz H-J, Lackner D, Dimroth F, Hannappel T. 2015 Efficient direct solar-to-hydrogen conversion by in situ interface transformation of a tandem structure. Nat. Commun. 6,8286. (doi:10.1038/ncomms

3. Burkhardt J, Patyk A, Tanguy P, Retzke C. 2016 Hydrogen mobility from wind energy—a life cycle assessment focusing on the fuel supply. Appl. Energ. 181,54-64. doi:10.1016/j.apenergy.2016.07.104)

4. Zeng M, Li Y. 2015 Recent advances in heterogeneous electrocatalysts for the hydrogen evolution reaction. J. Mater. Chem. A. 3,

14 942-14 962. (doi:10.1039/c5ta02974k)

5. Morales-Guio CG, Stern L-A, Hu X. 2014 Nanostructured hydrotreating catalysts for electrochemical hydrogen evolution. Chem. Soc. Rev. 43,6555-6569. (doi:10.1039/c3cs60468c)

6. LeeC-P etal. 2016 Beaded stream-like CoSe2 nanoneedle array for efficient hydrogen evolution electrocatalysis. J. Mater. Chem. A 4,4553-4561. (doi:10.1039/c6ta00464d)

7. Zou X, Zhang Y.2015 Noble metal-free hydrogen evolution catalysts for water splitting. Chem. Soc. Rev. 44,5148-5180. (doi:10.1039/c4cs00448e)

8. Chen WF, Sasaki K, MaC, Frenkel AI, MarinkovicN, Muckerman JT, Zhu YM, Adzic RR. 2012 Hydrogen-evolution catalysts based on non-noble metal nickel-molybdenum nitride nanosheets. Angew. Chem. Int. Ed. SI, 6131-6135. (doi:10.1002/ anie.201200699)

9. Chen W-F, Muckerman JT, Fujita E. 2013 Recent developments in transition metal carbides and nitrides as hydrogen evolution electrocatalysts.

Chem. Commun. 49,8896-8909. (doi:10.1039/ C3CC44076A)

10. Shi Y, Zhang B.2016 Recent advances in transition metal phosphide nanomaterials: synthesis and applications in hydrogen evolution reaction. Chem. Soc. Rev. 45,1529-1541. (doi:10.1039/c5cs00434a)

11. KibsgaardJ,Tsai C,Chan K, BenckJD, NorskovJK, Abild-Pedersen F, Jaramillo TF. 2015 Designing an improved transition metal phosphide catalyst for hydrogen evolution using experimental and theoretical trends. Energy Environ. Sci. 8, 3022-3029. (doi:10.1039/c5ee02179k)

12. ZhuoJ, Caban-AcevedoM, Liang H,Samad L, Ding Q, FuY, Li M, Jin S.2015 High-performance electrocatalysis for hydrogen evolution reaction using Se-doped pyrite-phase nickel diphosphide nanostructures. ACSCatal. 5,6355-6361. (doi:10.1021/acscatal.5b01657)

13. Popczun EJ, Read CG, Roske CW, Lewis NS, Schaak RE. 2014 Highly active electrocatalysis of the hydrogen evolution reaction by cobalt phosphide nanoparticles. Angew. Chem. Int. Ed. 53,5427-5430. (doi:10.1002/anie.201402646)

14. Popczun EJ, Roske CW, Read CG, Crompton JC, McEnaney JM, Callejas JF, Lewis NS, Schaak RE. 2015 Highly branched cobalt phosphide nanostructures for hydrogen-evolution electrocatalysis. J. Mater. Chem. A 3,5420-5425. (doi:10.1039/c4ta06642a)

15. Callejas JF, Read CG, Popczun EJ, McEnaney JM, Schaak RE. 2015 Nanostructured Co2P electrocatalyst for the hydrogen evolution reaction and direct comparison with morphologically equivalent CoP. Chem. Mater. 27,3769-3774. (doi:10.1021/acs.chemmater.5b01284)

16. ChenX, WangD, WangZ, ZhouP, WuZ, Jiang F. 2014 Molybdenum phosphide: a new highly efficient catalyst for the electrochemical hydrogen evolution reaction. Chem. Commun. 50,

11683-11 685. (doi:10.1039/c4cc05936k)

17. Xiao P, Sk MA, Thia L, Ge X, Lim RJ, Wang J-Y, Lim KH, Wang X. 2014 Molybdenum phosphide as an efficient electrocatalyst for the hydrogen evolution reaction. Energy Environ. Sci. 7,2624-2629. (doi:10.1039/c4ee00957f)

18. McEnaney JM, Crompton JC, Callejas JF, Popczun EJ, Biacchi AJ, Lewis NS, Schaak RE. 2014 Amorphous molybdenum phosphide nanoparticles for electrocatalytic hydrogen evolution. Chem. Mater. 26,4826-4831. (doi:10.1021/cm502035s)

19. Xing Z, Liu Q, Asiri AM, Sun X. 2014 Closely interconnected network of molybdenum phosphide nanoparticles: a highly efficient electrocatalyst for generating hydrogen from water. Adv. Mater. 26, 5702 -5707. (doi:10.1002/adma.201401692)

20. Wang D, Zhang D, Tang C, Zhou P, Wu Z, Fang B. 2016 Hydrogen evolution catalyzed by cobalt-promoted molybdenum phosphide nanoparticles. Catal.Sci. Technol. 6,1952-1956. (doi:10.1039/c5cy01457c)

21. Tymoczko J, Calle-Vallejo F, Schuhmann W, Bandarenka AS. 2016 Making the hydrogen evolution reaction in polymer electrolyte membrane electrolysers even faster. Nat. Commun. 7,10990. (doi:10.1038/ncomms10990)

22. Jeremiasse AW, Bergsma J, Kleijn JM, Saakes M, BuismanCJN, CohenStuart M, Hamelers HVM. 2011 Performance of metal alloys as hydrogen evolution reaction catalysts in a microbial electrolysis cell. Int. J. Hydrogen Energy 36,10 482-10 489. (doi:10.1016/ j.ijhydene.2011.06.013)

23. Lu Q etal. 2015 Highly porous non-precious bimetallic electrocatalysts for efficient hydrogen evolution. Nat. Commun. 6,6567. (doi:10.1038/ ncomms7567]

24. Koboski KR, Nelsen EF, HamptonJR.2013 Hydrogen evolution reaction measurements of dealloyed porous NiCu. Nanoscale Res. Lett. 8,528. (doi:10.1186/1556-276x-8-528)

25. Tang MH,HahnC, KlobucharAJ, NgJWD, Wellendorff J, Bligaard T,Jaramillo TF.2014 Nickel-silver alloy electrocatalysts for hydrogen evolutionandoxidationinanalkalineelectrolyte. Phys. Chem. Chem. Phys. 16,19 250-19 257. (doi:10.1039/c4cp01385a)

26. Cao B, Veith GM, Neuefeind JC, Adzic RR, Khalifah PG.2013 Mixed close-packed cobalt molybdenum nitrides as non-noble metal electrocatalysts for the hydrogen evolution reaction. J. Am. Chem. Soc. 135, 19 1 86-19192. (doi:10.1021/ja4081056)

27. Staszak-Jirkovsky J etal. 2016 Design of active and stable Co-Mo-Sxchalcogels as pH-universal catalysts for the hydrogen evolution reaction. Nat. Mater. 15,197-203. (doi:10.1038/nmat

28. Jiang J, Gao M, Sheng W, Yan Y. 2016 Hollow Chevrel-phase NiMo3S4 for hydrogen evolution in alkaline electrolytes. Angew. Chem. Int. Ed. 55, 15 240-15 245. (doi:10.1002/anie .201607651)

29. PanY, LinY, ChenY, LiuY, LiuC. 2016Cobalt phosphide-based electrocatalysts: synthesis and phase catalytic activity comparison for hydrogen evolution. J. Mater. Chem. A 4,4745-4754. (doi:10.1039/C6TA00575F)

30. Behranginia A et al. 2016 Highly efficient hydrogen evolution reaction using crystalline layered three-dimensional molybdenum disulfides grown on graphene film. Chem. Mater. 28,549-555. (doi:10.1021/acs.chemmater.5b03997)

31. JaramilloTF,j0rgensen KP, BondeJ, Nielsen JH, HorchS, Chorkendorff 1.2007 Identificationofactive edge sites for electrochemical H2 evolution from MoS2 nanocatalysts. Science 317,100-102. (doi:10.1126/science.1141483)

32. Kibsgaard J, Chen ZB, Reinecke BN, Jaramillo TF. 2012 Engineering the surface structure of MoS2 to preferentially expose active edge sites for electrocatalysis. Nat. Mater. 11, 963-969. (doi:10.1038/nmat3439)

33. Xie J, Zhang H, Li S, Wang R, Sun X, Zhou M, Zhou J, Lou XW, Xie Y. 2013 Defect-rich MoS2 ultrathin nanosheets with additional active edge sites for enhanced electrocatalytic hydrogen evolution. Adv. Mater. 25,5807-5813. (doi:10.1002/adma. 201302685)

34. Karunadasa HI,Montalvo E,SunY,Majda M,Long JR, Chang CJ. 2012 A molecular MoS2 edge site mimicfor catalytic hydrogen generation. Science 335,698. (doi:10.1126/science.1215868)

35. PuZ,SaanaAmiinuI,WangM,YangY,MuS.2016 Semimetallic MoP2: an active and stable hydrogen evolution electrocatalyst over the whole pH range. Nanoscale 8,8500-8504. (doi:10.1039/c6nr00 820h)

36. FaberMS, Lukowski MA, Ding Q, Kaiser NS,JinS. 2014 Earth-abundant metal pyrites (FeS2,CoS2, NiS2, and their alloys) for highly efficient hydrogen evolution and polysulfide reduction electrocatalysis. J. Phys. Chem. C. 118,21 347-21356. (doi:10.1021/jp506288w)

37. Guérin R, Sergent M. 1976 Nouveau phosphure ternaire: NiMoP2, a chaine linéaire infinie -Mo-Ni-Mo- et composés isotypes: NiWP2, CoMoP2 et CoWP2. J. Solid State Chem. 18,317-323. (doi:10.1016/0022-4596(76)90113-4)

38. Li Y, WangH, Xie L, LiangY, HongG, Dai H. 2011 MoS2 nanoparticles grown on graphene: an advanced catalyst for the hydrogen evolution reaction. J. Am. Chem. Soc. 133,7296-7299. (doi:10.1021/ja201269b)

39. LiaoL etal.2014Ananoporousmolybdenum carbide nanowire as an electrocatalyst for hydrogen evolution reaction. Energy Environ. Sci. 7,387-392. (doi:10.1039/c3ee42441c)