Ruthenium-cobalt nanoalloys encapsulated in nitrogen-doped graphene as active electrocatalysts for producing hydrogen in alkaline media.

The scalable production of hydrogen could conveniently be realized by alkaline water electrolysis. Currently, the major challenge confronting hydrogen evolution reaction (HER) is lacking inexpensive alternatives to platinum-based electrocatalysts. Here we report a high-efficient and stable electrocatalyst composed of ruthenium and cobalt bimetallic nanoalloy encapsulated in nitrogen-doped graphene layers. The catalysts display remarkable performance with low overpotentials of only 28 and 218 mV at 10 and 100 mA cm-2, respectively, and excellent stability of 10,000 cycles. Ruthenium is the cheapest platinum-group metal and its amount in the catalyst is only 3.58 wt.%, showing the catalyst high activity at a very competitive price. Density functional theory calculations reveal that the introduction of ruthenium atoms into cobalt core can improve the efficiency of electron transfer from alloy core to graphene shell, beneficial for enhancing carbon-hydrogen bond, thereby lowing ΔGH* of HER.

Hydrogen energy, a renewable energy resource, has been intensely investigated as an ideal alternative to the n class="Chemical">conservative fossil fuels for its high gravimetric energy density, zero-emission and earth-abundance123. Despite its promising application prospect, the current industrial route of hydrogen production not only consumes the fossil fuels but also brings on the greenhouse gas CO2 emission45. Thus, exploring a cleaner, more sustainable and efficient strategy is urgently needed. Currently, electrochemical water splitting, a non-fossil fuel-based technology, is evoking increasing interests and stimulating intense investigations to produce low-costing and high-pure hydrogen67. As the heart of the scalable hydrogen production, it is of great importance to develop highly efficient electrocatalysts to drive the hydrogen evolution reaction (HER)8.

It is well known, electrochemical water splitting can be carried out in acidic or alkaline media24. Unfortunately, the acid electrolyzers are technologically and commercially hindered by the lacking of efficient and low-cost counter electrode catalysts in acidic electrolytes2910. Consequently, tremendous efforts have been devoted to developing HER catalysts with high activity and stability in basic media on the basis of the available alkaline oxygen evolution electrocatalysts, aiming at accelerating commercialization of the basic electrolyzers for H2-production.

Platinum, as ‘the Holy Grail’ of HER electrocatalysts, remains the best HER catalysts with nearly zero overpotential and excellent long-term durability3111213. Unfortunately, the widespread commercialization of Pt-based electrocatalysts are hindered by their scarcity and expensive price1415. Therefore, it is of great significance to explore inexpensive alternatives for Pt electrocatalysts. To date, plenty of robust and efficient alternative catalysts have been reported aiming at replacing Pt-based electrocatalysts in alkaline media. Among these materials, transition metal (TM) -based catalysts, including Mo-based catalysts2491617, Ni-based catalysts6101819 and Co-based catalysts14202122, have been proven to be competitive eletrocatalysts as a result of their high efficiency and low cost. In particular, the biphasic nanocrystalline Ni–Mo–N catalyst recently reported by Li’s group exhibited amazing HER activity with low overpotentials of 43 mV in 1 M KOH and 53 mV in 0.5 M H2SO4, which is close to commercial Pt/C catalyst at the same mass loading of 1 mg cm−2 (ref. 2). Although, some of the TM-based catalysts show high HER performance, but still inferior to the Pt-based catalysts in overpotential and durability, which would increase energy consumption and hence decreasing economic competitiveness2324. Currently, the pure TM-based catalysts are unable to meet the requirement of replacing Pt-based electrocatalysts.

For decades, other cheaper platinum-group n class="Chemical">metals have been investigated for HER in view of their high similarity to Pt in chemical inertness13252627. In particular, ruthenium (42 $ per oz) is more economically advantageous in price than the rest of Pt-group metals, such as Pt (992 $ per oz), Pd (551 $ per oz), Ir (500 $ per oz) and so on28. As a matter of fact, Ru has evoked special attention as a top oxygen evolution electrocatalytic material29. However, it is not very active for HER in basic media in previous works30313233. It is shown that alloying noble metals with other TMs, with the amount of noble metal even decreased by up to an order of magnitude, is a major route to prepare highly efficient catalysts with balance of good cost-competitiveness. Moreover, previous studies have proven that the chemical properties of bimetallic surfaces could be modified by the combined changes in the average energy of the surface d-band and in the width of the d-band due to the cumulative strain and ligand effects originated from the formation of heteroatom bonds and the alteration of the bond length834. Therefore, in addition to lowering the material cost, the noble-transition bimetallic alloys could remarkably boost the electrocatalytic activities contributed by the shift of charge distributions and the resulting modification of surface properties during the formation of alloys263536.

In this work, a Ru-based electrochemical catalyst for HER was developed by alloy n class="Chemical">Ru with TM Co. Pervious work of our group and Bao’s group have shown that a metal core coated with a carbon shell, especially N-doped graphene, can simultaneously promote the HER activity through the synergism and enhance the stability due to the protection from carbon cage837. Inspired by this, a metal-organic frameworks (MOFs) -assisted strategy was adopted for the preparation of RuCo nanoalloys encapsulated in nitrogen-doped graphene layers (RuCo@NC). We adopted an in situ method that consists of the one-step annealing of Ru-doped Prussian blue analogues, which are ideal nitrogen-rich precursors composed of metals as nodes and CN- groups as linkers, for the fabrication of alloy materials wrapped in N-doped carbon82638394041. The RuCo@NC hybrid material shows an unprecedented high electrocatalytic performance towards HER in basic conditions, even superior to the commercial Pt/C catalysts.

Results

Synthesis and characterization of RuCo@NC catalyst

Obviously, Fig. 1 illustrates the synthetic route and model of the Run class="Chemical">Co@NC hybrids. The Co3[Co(CN)6]2 MOF precursor particles, designated as S-0-MOF, were synthesized according to our previous studies4142. As revealed by the field-emission scanning electron microscopy (FESEM) and transmission electron microscopy (TEM) images (Supplementary Fig. 1a,b), the morphologies of the S-0-MOF particles were truncated nanocubes which have very narrow diameter distributions with a mean diameter of ∼100 nm. In Fig. 1b, the Ru-doped Co3[Co(CN)6]2 precursor were obtained via an ion-exchange reaction in the liquid phase. In brief, various RuCl3 solution was added into the Co3[Co(CN)6]2 solution, and the corresponding products were hereinafter designated as S-1-MOF, S-2-MOF, S-3-MOF, S-4-MOF, S-5-MOF and S-6-MOF, respectively. After adding Ru source, Ru3+ ion diffused into the open framework of Co3[Co(CN)6]2, which triggered an ion-exchange reaction between Ru3+ and Co3+ while maintaining intrinsic framework structure4344. The FESEM and TEM images of the obtained Ru-doped precursor particles were illustrated respectively in Fig. 2a,b and Supplementary Fig. 2–6a,b, exhibiting that the as-prepared noble metal-doped Co3[Co(CN)6]2 inherited the nanocubic morphology with an average diameter of ∼100 nm. The corresponding X-ray diffraction patterns are shown in Fig. 2c. The patterns of the samples all showed no additional reflections except a series of Bragg reflections corresponding to the diffractions from the Co3[Co(CN)6]2 (JCPDS No. 77−1161). The unchanged X-ray diffraction patterns of Ru-doped MOFs suggest that Co3[Co(CN)6]2 could maintain the intrinsic framework structure during the Ru doping process, which is in good agreement with the unchanged morphologies of the precursors discussed above. The X-ray photoelectron spectroscopy (XPS) spectra of S-4-MOF were revealed in Fig. 2d and Supplementary Fig. 9. Obviously, the Ru 3p peaks at 463.3 and 485.8 eV were, respectively, assigned to Ru 3p 3/2 and Ru 3p 1/2 for Ru(III) species, conforming that Ru3+ have taken the place of some Co3+ ions. Essentially, the Ru doping process was carried out in the MOF framework structure via an ion-exchange reaction43.

Figure 1

Schematic illustration of the synthetic route and model of the RuCo nanoalloys encapsulated in nitrogen-doped graphene layers.

(a) a nanocube of Co3[Co(CN)6]2 precursor, (b) a nanocube of Ru-doped Co3[Co(CN)6]2, (c) an aggregate of RuCo alloys encapsulated in graphene layers and (d) enlarged model of RuCo nanoalloy encapsulated in nitrogen-doped graphene layers as an electrocatalyst towards hydrogen evolution reaction in alkaline media.

Figure 2

Ru-doped MOF Co3[Co(CN)6]2 MOF precursor particles.

(a,b) FESEM and TEM images of the as-prepared S-4-MOF. Scale bars, 100 nm. (c) The X-ray diffraction patterns of S-0-MOF, S-1-MOF, S-2-MOF, S-3-MOF, S-4-MOF, S-5-MOF and S-6-MOF. (d) The XPS result of the Ru 3p spectrum enlarged from Supplementary Fig. 9a.

In this work, the electrocatalysts were synthesized via an ‘in n class="Chemical">situ’ formation process without adding any other carbon sources. Co3[Co(CN)6]2 is composed of CN− groups as linkers and TM Co as metallic nodes, which is an ideal precursor for the preparation of composite materials of TM-based material with highly N-doped carbon39. Specifically, the corresponding RuCo@NC hybrids were synthesized via one-step annealing of the obtained Ru-doped MOFs, hereinafter marked as S-0, S-1, S-2, S-3, S-4, S-5 and S-6, respectively. As illustrated in Fig. 1c,d, the nanocubic MOF precursors were directly carbonized at 600 °C under a nitrogen flow without adding any other carbon sources. The resulting product is composed of the bimetallic RuCo nanoalloys encapsulated in nitrogen-doped graphene layers. During the annealing process, Ru and Co atoms from the precursor would form the bimetallic RuCo nanocrystals, meanwhile some remaining CN− group linkers would transfer into nitrogen-doped graphene layers. As can be seen in Fig. 1d, the formation of RuCo nanoalloys will be coated speedily by the in situ formed N-doped graphene layers, which is able to efficiently avoid the agglomeration of the inside alloy particles to provide more eletrocatalytic active sites which are beneficial to HER activity and long-term corrosion protection to enhance the catalytic durability.

All the obtained samples embraced an irregular particle-like morphology according to the FESEM results (Fig. 3a and Supplementary Fig. 2–6d). As revealed in the TEM images (Fig. 3b, Supplementary Figs 2–6c,e,f and 7 and Supplementary Table 1), the irregular products were composed of small alloy particles with a mean diameter of ∼30 nm encapsulated carbon layers. The high-resolution TEM (HRTEM) images (Fig. 3c,d) showed that the small alloy particles were coated with N-doped graphene layers. Most of the graphene layers were ∼6–15 layers thick. However, as shown in Supplementary Fig. 8, a large proportion of graphene shells consisting of 1–5 layers still existed, which are good active sites for HER. The thin graphene shells were beneficial in electron transfer from the alloy core to the shell and hence improved the catalytic activity on the graphene surface837. Besides, Fig. 3e clearly revealed that lattice fringe spaces of 0.251 and 0.204 nm are, respectively, consistent with the (111) and (220) planes of the cubic Co phase or the (002) and (220) planes of the hexagonal Co phase. While, no lattice fringe spaces of metallic Ru were found, implying the formation of RuCo alloy which kept the crystal structure of metallic Co. As shown in Fig. 3f–j, the images of elemental mapping from energy filtered TEM showed that the Co and Ru elements (5.69 wt.% Ru and 94.31 wt.% Co) were uniformly distributed in the inner particles and surrounded by the C and N elements, further confirming the formation of RuCo nanoalloys encapsulated in nitrogen-doped graphene layers.

Figure 3

SEM and TEM characterization of RuCo@NC hybrid.

(a,b) The FESEM and TEM images of S-4. Scale bars, 100 nm. (c–e) HRTEM images of S-4. Scale bars, 10 nm. (f–j) HAADF-STEM image and corresponding EDX maps of S-4 for Co (g), Ru (h), C (i) and N (j), respectively. Scale bars, 20 nm.

The corresponding X-ray diffraction patterns of n class="Chemical">RuCo@NC were shown in Fig. 4a. Similar diffraction features of a broad and weak peak of C (002) were detected, confirming the existence of carbon layers. Besides, the patterns of all the samples showed no additional reflections except a series of Bragg reflections corresponding to the diffractions from the hexagonal structure Co (JCPDS card no. 05−0727) phase and face centred cubic (FCC) Co (JCPDS card no. 15-0806) (refs 7, 26), suggesting that Ru atoms were in situ dissolved in metallic Co to form RuCo alloy, which is in good agreement with the value determined by the HRTEM observations discussed above.

Figure 4

Structural analysis of RuCo@NC hybrid.

(a) The X-ray diffraction patterns of S-0, S-1, S-2, S-3, S-4, S-5 and S-6, respectively. (b) Raman spectrum of S-4, (c) The XPS result of the N1s spectrum enlarged from Supplementary Fig. 10a, (d) N2 adsorption–desorption isotherm and pore size distribution plot (inset) of S-4.

The Raman spectra of S-4 were shown in Fig. 4b. The product displayed three Raman peaks locating at ∼1,349, 1,583 and 2,703 cm−1, which corresponded to the D, G and 2D bands, respectively. The high ID/IG band intenn class="Chemical">sity ratio of S-4 indicated the generation of large amounts of defects, suggesting that a large amount of N atoms were doped in the graphitic carbon layers. Moreover, the second-order band is broad and weak, implying that the coated carbon is thin graphene with several layers45.

The XPS spectra of S-4 were illustrated in Supplementary Fig. 10a. The nitrogen n class="Chemical">content of S-4 was ∼3.51 atom%. The high-resolution N1s spectrum of S-650 (as shown in Fig. 4c) can be deconvoluted into four individual peaks assigned to pyridinic-N (398.5 eV), pyrrolic-N (399.6, 400.9 eV) and quaternary-N (401.6 eV), respectively846. Importantly, two pyrrolic-N-binding energies are observed here. This phenomenon might result from the energy shift which was induced by the interaction between some pyrrolic-N and metal atoms47. Therefore, the percentages of pyridinic-N, metal-N, pyrrolic-N and quaternary-N were 40%, 17%, 27% and 16%, respectively. Besides, XPS investigation of the Co 2p spectrum revealed the presence of two distinct chemical species: Co0 and CoII species peaks, indicating the existence of CoII originated from surface oxidation of metallic Co. The existence of metallic Ru was also confirmed by the weak peak at 462.2 eV (Supplementary Fig. 10c), suggesting that the Ru content was very tiny in the RuCo alloys with good balance of its cost.

The specific surface area and pore size distribution were obtained by n class="Chemical">N2 adsorption/desorption isotherms. As can be seen in Fig. 4d, a type-II isotherm with a H3-type hysteresis loop is obtained, which is characteristic of mesoporous non-rigid aggregates48. Obviously, S-4 have a specific surface area of 41.365 m2 g−1. The pore size distribution (insert image) revealed a strong and narrow distribution centred at ∼3 nm. The other samples showed the similar adsorption/desorption isotherm curves and pore size distribution in Supplementary Fig. 11. The BET surface area were listed in Supplementary Table 1. These results were in good agreement with the value determined by the morphology observations discussed above.

Electrochemical characterization for HER catalysis

The electrochemical catalytic activities of the annealed samples for HER were recorded by a typical three-electrode electrochemical cell in N2 saturated 1 M KOH electrolyte. The HER performance was evaluated by observing the overpotential versus the reversible hydrogen electrode (RHE) at 10 mA cm−2, which is the current density expected for a 12.3% efficient solar water-splitting device15. As revealed in Fig. 5a, the undoped S-0 catalyst exhibited an excellent HER activity but far inferior to the doped counterparts of RuCo alloy encapsulated in N-doped graphene layers, proving that alloying TMs with less noble metals played a key role in boosting the electrocatalysis activity with good balance of the cost. The polarization curves also showed that the S-4 catalyst had the highest HER activity among seven catalysts with an overpotential of only 28 mV without IR-correction. As illustrated in Table 1, the trend in the overpotentials was found to be S-4 (28 mV)sing amount of Ru. Interestingly, the activity would decrease with the further increase of Ru amount. Besides, the similar trend of overpotentials could also be observed reaching a higher current density of 100 mA cm−2 (Table 1 and Supplementary Fig. 12b). To gain further insight into the activity of S-4, the inner metal of S-4 was etched via 1 M HCl solution. Obviously, S-4 retains the same morphology, but a part of its metallic cores are removed after the acid leaching (Supplementary Fig. 13a). Besides, the HCl solution became pink in Supplementary Fig. 13c, suggesting that just a part of metallic Co was corroded and dissolved into the solution. It was also testified by the neglectable Ru concentration (0.46 wt.% Ru in Co and Ru) of ICP. It was shown in Fig. 5d and Supplementary Fig. 13d that both of the overpotentials reaching current densities of 10 and 100 mA cm−2 increased after etching, revealing that the metallic Co in S-4 was indispensable and had a key role in HER activity. It was also found that S-4 with lower Ru content (3.58 wt.% Ru in RuCo alloy, obtained from ICP) showed better activity than the etched counterpart with higher Ru content (17.7 wt.% Ru in RuCo alloy, obtained from ICP), which was in good agreement with the above- mentioned trend of overpotentials. All these findings revealed that S-4 (Ru, 3.58 wt. %, obtained from ICP) hold the best alloy structure for HER among RuCo alloys.

Figure 5

Electrocatalytic HER performance test of catalysts in N2 saturated 1 M KOH solution.

(a) HER polarization curves of RuCo@NC samples with the same mass loading, (b) HER polarization curves of S-4, Ru and Pt/C with the same mass loading and durability test after 10,000th cycles, (c) the Tafel plots of S-4, Ru and Pt/C, (d) HER polarization curve of etched S-4 by 1 M HCl.

Table 1

Comparison of HER activity of different catalysts.

Sample	S-0	S-1	S-2	S-3	S-4	S-5	S-6	Pt/C	Ru
η@10 (mV)	300	91	83	67	28	91	85	40	58
η@100 (mV)		338	322	291	218	281	284	287	462
Ru (wt.%)	0	1.858	2.316	3.074	3.58	4.00	4.234	20 (Pt )	99.9

HER, hydrogen evolution reaction.

The electrochemically active surface area of samples was estimated un class="Chemical">sing a simple cyclic voltammetry (CV) method949. Due to the unknown capacitive behaviour (Cs) of the RuCo alloy electrode especially with N-doped graphene shell, we can safely estimate relative surface areas of seven samples, since the double-layer capacitance (Cdl) is expected to be linearly proportional to effective active surface area for samples with similar composition and this method was also employed in previous study8. The results in Supplementary Fig. 14 suggested that with the increasing density of catalytically active sites of samples also follow the similar trend of the catalytic performance. A series of activity normalization of special activity (SA) and active site activity (ASA), taking into account of the surface area and active site concentration respectively at overpotential of 100 mV, were shown in Supplementary Fig. 16 (ref. 50). Especially, the S-4 catalyst exhibited a SA of 0.707 mA cm−2, which are much higher than other catalysts. Besides, the result also suggests that the SA of samples increased with the increasing catalytic active in general, suggesting the effect of surface area is a negligible factor on activity. Due to the unknown capacitive behaviour (Cs) of the RuCo alloy electrode especially with N-doped graphene shell, we carried out activity normalization of ASA*CS. Especially, the S-4 catalyst exhibited an ASA*CS of 1.242 mA cm−2, which are much higher than other catalysts with the similar value of ASA*CS. Besides, the electrochemical impedance spectroscopy images (Supplementary Fig. 15) showed that the Co@NC has a biggest semicircle radius, indicative of a slightly higher charge-transfer impedance than of S-1 (8.53 Ω), of S-2 (8.55 Ω), of S-3 (8.58 Ω), of S-4, of S-5 (8.53 Ω) and of S-6 (8.47 Ω) (refs 51, 52). Therefore, alloying Co with small amount of Ru could lead to a higher charge-transfer rate and more facile catalytic kinetics toward HER.

The commercial n class="Chemical">Ru powders catalyst (Ru: 99.9 wt.%) and the commercial Pt/C catalyst (Pt: 20 wt.%) were also measured as references. Surprisingly, the obtained S-4 even exhibited better activity than Ru catalyst (58 mV) and Pt/C catalyst (40 mV) (Fig. 5b). Moreover, the fast HER kinetics also enabled the S-4 catalyst to reach a high current density of 100 mA cm−2 at η as low as 0.218 V, outperforming the other HER catalysts including Pt/C (Supplementary Fig. 12a and Table 1). It is also the best one among some catalysts in the recent reports towards HER in basic media (Table 2)1246917181921222553545556. The Tafel slope (Fig. 5c) of the S-4 sample was 31 mV dec−1, which was as same as the value of Pt/C, implying a rapid HER rate and a Tafel–Volmer mechanism with electrochemical desorption of H2 as the rate-determining step in the HER process57. While, the Tafel slope was 55 mV for the Ru powder catalyst, which was much larger and indicative of a typical Volmer−Heyrovsky mechanism with the Volmer step as rate-limiting step for HER653. The exchange current density (j0) values were obtained from Tafel plots29. As read from Supplementary Fig. 17, the catalysts showed an increase in jo in the following order: S-4 (10−2.48)> Ru powder (10−2.74)> Pt/C (10−2.94). Durability was another critical parameter to assess the electrocatalytic performance. The durability of S-4 was also evaluated by measuring polarization curves after 10,000 CV sweeps between −1.1 and −0.9 V (versus Ag/AgCl) at 100 mV s−1. As illustrated in Fig. 5b, the polarization curve of S-4 after 10,000 cycles retained an almost similar performance to the initial test, just with the overpotential increased by 4 mV, which is also slightly superior to the durability of Pt/C (8 mV) as well as much better than Ru powder. The graphene shell and the alloying of Ru and Co made great contributions to protect the catalysts from corrosion during cycling837. Besides, Ru (42 $ per oz) is further more economically advantageous in price than Pt (992 $ per oz)28. The Ru content of S-4 was 3.58 wt. % obtained from ICP in Table 1, which was much lower than the Pt content of 20 wt. % in Pt/C catalysts. Specific to the field of noble metal’s cost, the S-4 was only 0.76% of commercial Pt/C in price. Therefore, the S-4 sample performed higher activity, durability and economic competitiveness than the commercial Pt/C catalysts, demonstrating that it had a great potential to be a substitution of Pt/C catalyst for HER in alkaline media.

Table 2

HER electrocatalysts in alkalic media reported recently.

Catalyst	electrode	Loading amount (mg cm−2)	Electrolyte	Overpotential at 10 mA cm−2 (mV)	Tafel plots (mV dec−1)	Reference
S-4	GHE	0.275	1 M KOH	28	31	This work
NiO/Ni-CNT	GHE	0.28	1 M KOH	80	82	Nat. Commun., 20146
Ni/NixPy	Ni foam	N/A	1 M KOH	130 (with iR-correction)	58.5	Adv. Funct. Mater., 201618
NF-Ni3Se2/Ni	Ni foam	8.87	1 M KOH	203 (with iR-correction)	79	Nano Energy, 201619
Pt13Cu73Ni14/CNF@CF	CNF@CF	N/A	1 M KOH	148 (η5)	54	ACS Appl Mater Interfaces, 201621
Pd-CNx	GCE	0.28	0.5 M KOH	180 (η5)	150	ACS Catalysis, 201625
Mo2C@N-C	GCE	0.28	1 M KOH	60	N/A	Angew Chem Int Ed Engl, 20154
Ni–Mo–N	GCE	1	1 M KOH	43 (η20)	40	Nano energy, 20162
MoO2/CC	Carbon paper	2.9	1 M KOH	100 (with iR-correction)	41	Adv. Mater., 20169
MoOx/Ni3S2 /NF	Ni foam	12	1 M KOH	110 (η15)	90	Adv. Funct. Mater., 201617
CoOx@CN	Ni foam	0.42	1 M KOH	232	N/A	J. Am. Chem. Soc., 201556
CoP2/RGO	GCE	0.285	1 M KOH	88(with iR-correction)	50	J. Mater. Chem. A, 201653
CoP/rGO-400	GCE	0.28	1 M KOH	150	38	Chem. Sci., 201654
c-CoSe2/CC	Carbon cloth	N/A	1 M KOH	190(with iR-correction)	85	Adv. Mater., 201622
N-Co@G	GCE	0.285	0.1 M NaOH	337	N/A	ACS Appl Mater Interfaces, 201555
WC-CNTs	Si wafer	N/A	0.1 M KOH	137	106	ACS Nano, 20151

HER, hydrogen evolution reaction.

Some of the information was not specified in the literature and was estimated according to the data graphs.

HER enhancement mechanism

To investigate the origin of excellent activity of RuCo alloy encapsulated in nitrogen class="Chemical">n-doped graphene layers for HER, density functional theory calculations were carried out using software of Vienna Ab Initio Simulation Package. A graphitic carbon cage C240 encapsulated 55 metal atoms was used as the basic model of graphene-encapsulated alloys, which worked well in previous studies3758. The detailed calculation information could be read in the Supplementary Note 1. In general, for HER performed in both acid and alkaline electrolyte, ΔGH* is one of the key descriptors in theoretical prediction of the activity for HER. Previous studies have used |ΔGH*| as a catalytic descriptor for HER and proposed the optimal value should be close to 0 (refs 22, 28). Therefore, the ΔGH* of different models were calculated, including pure graphene (C240), graphene doped with nitrogen atoms (C239N1), graphene-encapsulated Co (C239N1Co55) as well as graphene-encapsulated RuCo alloys. In order to find the influence of Ru content in RuCo alloys, three established different models were abbreviated to Ru1Co (C239N1Ru1Co54), Ru2Co (C239N1Ru2Co53) and Ru3Co (C239N1Ru3Co52), in which the metal atom ratio is comparable to experimental materials measured by ICP. The optimized adsorption structures are shown in Fig. 6, which exhibited the C atom next to the doping N atom in graphene shell is the active adsorption site for H*, which is also consistent with previous work83753. As illustrated in Fig. 7, the nitrogen doping and the combination of metal Co with graphene can significantly reduce ΔGH*, this result was consistent with our previous work8. Moreover, the introduction of Ru atoms into Co metal can further decrease ΔGH* and properly increasing Ru content will further reduce the value of ΔGH*. As a result, Ru3Co model exhibit the lowest ΔGH* (0.31 eV) among graphene-encapsulated models. We also calculated H* adsorbed on Pt55 and Ru55 metal cluster models. Besides, some possible H* adsorption sites such as top, bridge, hcp and hollow sites are taken into considerations in our cluster models (Supplementary Figs 19 and 20). It should be noted that, during optimization process the hollow site of Pt55 and top site of Ru55 are tended to be bridge and hollow sites, respectively, indicating the latter are more moderate adsorption sites. The calculated ΔGH* of different models and sites are illustrated in Supplementary Tables 2 and 4.

Figure 6

Optimized structures of H* adsorbed on nitrogen-doped graphene-encapsulated Co and RuCo alloys.

(a) Pure Co model (C239N1Co55), (b) Ru1Co alloy model (C239N1Ru1Co54) and (c) Ru2Co alloy model (C239N1Ru2Co53) (d) Ru3Co alloy model (C239N1Ru3Co52). The black, blue, pink, brown and white balls refer to C, N, Co, Ru and H atoms, respectively.

Figure 7

HER-free energy diagram.

ΔGH* calculated at the equilibrium potential of different models.

The calculated results show n class="Chemical">Ru55 cluster exhibited even better ΔGH* than RuCo and Pt55, which is not in good agreement with our previous experimental results based on purchased Ru powder. Unlike bulk metal, metal cluster and nano-sized metal particles have more high active sites such as facets, corners, edges. However, the Ru powder measured before was purchased from manufacturer and the size of them was in micron order as illustrated in Supplementary Fig. 18a. Besides, recent researches report nano metallic Ru-based electrocatalysts with very active performance in basic media5960. Therefore, we synthesized nano-sized Ru nanoparticles (Supplementary Fig. 18b) by ourselves through RuCl3 and NaBH4 as Precursors. The performance of nano-sized Ru catalyst was even better than our S-4 at high overpotential (Supplementary Fig. 18d–f), which was in good agreement with our calculated cluster result. For Pt cluster, even its calculated result was as good as our CoRu models; however, the measured activity was inferior to our S-4 sample. Similar results also existed in a recently published research59, Qiao and his co-workers have proved that, apart from the ΔGH*, the water dissociation kinetics would also affect the overall reaction rate, especially under basic electrolyte. According to their calculated results, when the kinetics of water dissociation from the Volmer step is considered, Pt exhibits a significant higher energy barrier than Ru and others, indicating sluggish water dissociation during catalytic process. Therefore, from the kinetic viewpoint, conversely, Pt did not demonstrate as good catalytic performance as depicted by calculated ΔGH*, which might explain the inferior measured activity than our CoRu alloys. We also tried to calculate energy barrier of water dissociation process for our CoRu alloys; however, due to the very large models (C239N1Co52Ru3 containing 295 atoms), it’s very difficult and time-consuming to search for the exact transition state and energy barrier of ours. However, combing our calculated ΔGH* results with Qiao’s result, it’s sufficient to give a reasonable explanation. Therefore, both the experimental and the calculated results indicated that alloying Co with small amount of Ru could obtain highly active eletrocatalyst comparable to pure Ru catalyst.

The charge-density difn class="Chemical">ference of model Co and Ru3Co was also calculated and illustrated in Fig. 8. It is shown that electrons transferred from alloy core to C atoms of graphene shell, which was beneficial to enhance C–H bond thereby lowing ΔGH* of the model as a whole37. As shown in Supplementary Table 3, the exact number of transferred electron was calculated using bader charge analysis. Ru3Co model showed the largest number of transferred electrons (5.91e−) among graphene-encapsulated models, therefore, ΔGH* of it was also the lowest among them. The calculated result is consistent with the experimental results, proving the excellent HER activity originates from Ru alloying with Co atoms as well as the unique structure derived from MOFs.

Figure 8

Calculated charge-density differences of different models.

(a) Co and (b) Co3Ru models. The isosurface value of the colour region is 0.01 e Å−3. The yellow and cyan regions refer to increased and decreased charge distributions, respectively.

Discussion

In summary, a novel RuCo alloy catalysts with high catalytic activity and stability for HER in ban class="Chemical">sic solutions was prepared, which exceeds almost all the documented electrocatalysts including 20 wt.% Pt/C catalysts, hence enable it a cheaper alternative to Pt-based electrocatalysts for HER in basic media. The RuCo alloy structures with a low content of Ru provide it a huge economical advantage in price over commercial Pt/C catalysts. The density functional theory calculations indicate that RuCo alloy core can transfer more number of electrons to the graphene shell than pure Co metal core so as to enhance C–H bond, which will significantly decrease ΔGH*, and thereby improve electrocatalysis activity. The results here show a new way for the development of high-performance HER electrocatalysts in alkaline media while reducing the cost of noble metal electrocatalysts.

Methods

RuCo@NC hybrids synthesis

The Co3[Co(CN)6]2 MOF was prepared acn class="Chemical">cording to our previous researchs4142. The as-prepared Co3[Co(CN)6]2 nanoparticles (25 mg) were dispersed in a 20 ml distilled water system under agitated stirring to get an absolutely homogeneous mixed solution, followed by the addition of 0.55, 1.1, 1.65, 2.2, 2.75 or 3.3 ml RuCl3 solution (0.01 g ml−1). After agitated stirring for 10 h in dark, the brown products were collected and rinsed several times by distilled water, and finally dried under oven at 60 °C. The thermal decomposition of the MOF precursor was performed at 600 °C for 4 h under nitrogen atmosphere in the oven with a heating rate of 10 °C min−1. The obtained sample should be kept in the vacuum drying oven.

Nano Ru synthesis

Briefly, 100 mg sodium borohydride (n class="Chemical">NaBH4) was dissolved in a 20 ml distilled water system under agitated stirring to get an absolutely transparent solution. And then the above solution was added into 20 ml RuCl3 solution (5 mg ml−1) slowly. After agitated stirring for 20 h, the black products were collected and rinsed several times by distilled water and ethanol, and finally dried under oven at 60 °C in the vacuum drying oven.

Characterization

The powder X-ray diffraction patterns of the samples were analysed with an X-ray diffractometer (Japan Rigaku D/MAX-γA) un class="Chemical">sing Cu-Kα radiation (λ=1.54178 Å) with 2θ range of 20–80°. The morphology and size of all as-synthesized samples were characterized by a JEOL JSM-6700 M field-emission scanning electron microscopy and a scanning TEM (STEM; Talos F200X) with energy dispersive X-ray (EDX) spectroscopy. Raman spectrum was carried out using a LabRAM HR Evolution ranging from 1,000 to 3,000 cm−1. XPS was measured on an ESCALAB 250 X-ray photoelectron spectrometer using Al Ka radiation. The specific surface and pore diameters were obtained from the results of N2 physisorption at 77 K (Micromeritics ASAP 2020) by using the BET (Brunauer-Emmet-Teller) and BJH (Barrett-Joyner-Halenda) methods, respectively. The content of Co and Ru were obtained via the inductively coupled plasma-atomic emission spectrometer (ICP-AES) (Optima 7300 DV).

HER electrochemical measurements

The HER electrochemical perform measurements were performed in a three-electrode system on an electrochemical workstation (CHI 660D) in 1 M KOH electrolyte. Typically, 4 mg of catalyst and 30 μl n class="Chemical">Nafion solution (Sigma Aldrich, 5 wt.%) were dispersed in 1 ml of 3:1 v/v water/isopropanol mixed solvent by at least 30 min sonication to form a homogeneous ink. Then 5 μl of the dispersion (containing ∼19.4 μg of catalyst) was loaded onto a glassy carbon electrode with 3 mm diameter (loading ∼0.275 mg cm−2). An Ag/AgCl (filled with 3 M KCl solution) electrode and a platinum wire were served as the reference electrode and counter electrode, respectively. All of the potentials were calibrated to a RHE. The working electrode was polished with Al2O3 powders with size ranging from 1 to 0.05 μm. Commercial Ru powders catalysts (99.9%, Aladdin) and commercial Pt/C catalysts (20%, Alfa Aesar) were used as a reference to evaluate the electrocatalytic performance of various samples. Linear sweep voltammetry with a scan rate of 2 mV s−1 was conducted between −1.1 V and −0.9 V versus AgCl/Ag electrode into N2 saturated 1 M KOH electrolyte. CV was conducted in 1 M KOH solution in the potential region from −1.1 V to −0.9 V versus AgCl/Ag electrode at a sweep rate of 100 mV s−1 for 10,000 times to investigate the cycling stability.

Data availability

The authors declare that the data supporting data supporting the findings of this study are available within the article and its Supplementary Information files.

Additional information

How to cite this article: Su, J. et al. Ruthenium-Cobalt nanon class="Chemical">alloys encapsulated in nitrogen-doped graphene as active electrocatalysts for producing hydrogen in alkaline media. Nat. Commun. 8, 14969 doi: 10.1038/ncomms14969 (2017).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.