Abstract
Designing highly durable and active electrocatalysts applied in polymer electrolyte membrane (PEM) electrolyzer for the oxygen evolution reaction remains a grand challenge due to the high dissolution of catalysts in acidic electrolyte. Hindering formation of oxygen vacancies by tuning the electronic structure of catalysts to improve the durability and activity in acidic electrolyte was theoretically effective but rarely reported. Herein we demonstrated rationally tuning electronic structure of RuO2 with introducing W and Er, which significantly increased oxygen vacancy formation energy. The representative W0.2Er0.1Ru0.7O2-δ required a super-low overpotential of 168âmV (10âmA cmâ2) accompanied with a record stability of 500âh in acidic electrolyte. More remarkably, it could operate steadily for 120âh (100âmA cmâ2) in PEM device. Density functional theory calculations revealed co-doping of W and Er tuned electronic structure of RuO2 by charge redistribution, which significantly prohibited formation of soluble Rux>4 and lowered adsorption energies for oxygen intermediates.
Similar content being viewed by others
Introduction
Electrochemical splitting water into H2 is considered as a promising technology for storing renewable and sustainable energy1,2,3,4,5,6,7. Acidic proton exchange membrane (PEM) water electrolyzers is highly desirable because the reaction rate in acidic electrolyte is three or more orders of magnitudes speedier than pH-neutral and alkaline electrolyzers owing to its higher voltage efficiency, more compact designing system, lower ohmic loss, and wider partial load range8,9,10,11. For decades, anode catalyst for PEM is mainly concentrated on IrO2 with the compromise between stability and activity in acidic electrolyte12,13,14,15,16,17. However, the low mass activity and the high cost extremely hindered its practical applications at large scale. Compared with IrO2, Ru-based materials usually exhibit the relatively higher activity owing to the suitable binding ability of oxygenated intermediate species (*OH, *O, and *OOH) with the surface active sites18,19,20,21,22, but with a poor stability. Recently, a series of Ru-based oxides, including oxide perovskites23,24, 3d metals doped RuO225, heterostructured Ru@IrOx3 and Ru1âPt3Cu18 with compressive strain, were explored to simultaneously enhance their activities and stabilities in acidic electrolyte by regulating the compositions and surface strain. Although the activities of Ru-based catalysts are attractive, the stability is still only dozens of hours in acidic electrolyte. It has been proven that the over-oxidation of Ru-based materials results in soluble Rux>4 (e.g. RuO4) derivatives14,26,27, which seriously hinders the application of RuO2 in the commercial PEM devices.
The over-oxidation of Ru-based catalysts in acidic electrolyte is mainly caused by the oxidation of lattice oxygen, resulting in oxygen vacancies (VO) during OER9,18,28. The generation of VO would expose Ru atoms on the surface3,29,30,31,32, which will be over-oxidized to soluble high valence Rux>4 (e.g. RuO4) materials. The over-oxidation of Ru unavoidably leads to the collapse of the crystal structure, damaging the stability28,29. Therefore, it is reasonable to speculate that the over-oxidation of Ru-based catalysts could be thermodynamically hindered during OER if we can tune the electronic structure of RuO2 to make the formation energy of VO much higher than the redox H2O/O2 energy. Foreign elements doping is the classical and effective strategy to enlarge the localized gap between O 2p band centers and Fermi level28,33,34. This would enhance the energy barrier for oxidation of lattice oxygen, prohibiting the formation of VO.
Herein, W and Er are co-doped into the lattice of RuO2 to modify its electronic structure, avoiding over-oxidation of Ru. By down-shifting O 2p-band centers, the energy (ÎG) for the VO formation in WmErnRu1âmânO2âδ is significantly increased, preventing the generation of soluble high valence Rux>4. As a result, the representative W0.2Er0.1Ru0.7O2âδ offers a record low overpotential of 168âmV to achieve 10âmAâcmâ2 accompanied with a stability at least 500âh in 0.5âM H2SO4 electrolyte. It can also be applied as anode catalyst in acidic PEM with a high current of 100âmAâcmâ2 for over 120âh. Thus, tuning the electronic structure of RuO2 to prevent the over-oxidation and dissolution of Ru opens a feasible route to keep Ru4+ active sites for acidic OER. This work provides the insight to understand the stability of metal oxide catalysts in acidic electrolyte and a strategy to promote their activity and stability by fixing the catalysts.
Results and discussion
Mechanisms of W0.2Er0.1Ru0.7O2âδ toward acidic OER
As for the metal dissolution in acidic electrolyte, it is mainly correlated with the adsorbate evolution and lattice oxygen oxidation mechanisms9,18,35. For lattice oxygen oxidation, there are four electrochemical steps as well as a non-electrochemical step for desorption of O2. In lattice oxygen oxidation route, the origin for releasing oxygen was occurred in step III (O*â+âH2Oâ+âOLâ+â2(H+â+âeâ)âOVâ+âO2â+âH2Oâ+â2(H+â+âeâ), L represents lattice). With releasing oxygen from the lattice of electrocatalysts, the oxygen vacancies were created, in turn beneficially enhancing the OER activity (Fig. 1a, c)28. However, in acidic electrolyte, the existing oxygen vacancies would accelerate the over-oxidation of the exposed Ru, leading to the highly soluble Rux>4 (e.g. RuO4) derivatives. Therefore, the structure of Ru-based materials on anodic side would be destroyed with releasing lattice oxygen from catalysts. As compared with lattice oxygen oxidation, there are only four electrochemical steps in the adsorbate evolution way toward OER (Fig. 1b)9,36. Simultaneously, in the adsorbate evolution route, OâO coupling from H2O in steps I and III was the origin for releasing O2 (step I: W0.2Er0.1Ru0.7O2âδâ+âH2OâââW0.2Er0.1Ru0.7O2âδâOH*â+âH+â+âeâ, step III: W0.2Er0.1Ru0.7O2âδâO*â+âH2OâââW0.2Er0.1Ru0.7O2âδâOOH*â+âH+â+âeâ) (Fig. 1d). The structures of catalysts would not be destroyed without releasing lattice oxygen. Therefore, effectively enhancing the energy for formation of oxygen vacancies and inhibiting the direct OâO coupling together from the lattice of W0.2Er0.1Ru0.7O2âδ would be beneficial for maintaining the stability of Ru-based materials in acidic electrolyte (Fig. 1b, d).
Density functional theory (DFT) calculations
In light of the reacted ways toward acidic OER, we carried out DFT studies aiming at rationally tuning the electronic of RuO2. To enhance the acidic OER stability, the energy for VO formation should be thermodynamically suppressed. Density of states (DOS) from DFT was established to investigate how W and Er synergistically rationalize the Ru 4d and O 2p orbitals in RuO2 (Fig. 2a and Supplementary Fig. 1). The DOS spectra shows that O 2p band center (εp) moves from â3.31âeV (RuO2) to â4.12âeV with W and Er introduction. And the gap between εp and Fermi level is obviously enlarged28. This result indicates the covalency of RuâO bond is decreased by introduction of W and Er37. Therefore, in W0.2Er0.1Ru0.7O2âδ, the direct OâO coupling on O 2p states above the Fermi level will be not thermodynamically favored, which could be seen from ÎG2 for formation of OH* (Fig. 2b, right). This result will make ÎG for formation of VO higher than that of RuO2 (Fig. 2b, left). Further, the DFT results also reveal ÎG for formation of VO at the pyramid vertex of W0.2Er0.1Ru0.7O2âδ shifts from lower energy band of 0.67âeV (RuO2) to higher energy band of 2.29âeV (Fig. 2c). These results indicate W and Er could effectively enhance the energy for formation of VO in W0.2Er0.1Ru0.7O2âδ toward acidic OER. Therefore, the dissolution of Ru in acidic electrolyte will be suppressed. Moreover, Bader charge analyses also prove the downshift of p band centers, resulting in the negative charge of the lattice O (Supplementary Fig. 2). Consequently, this electronic state could effectively enhance ÎG for the formation of VO and suppress the lattice O binding with oxygen intermediates to release O2. Besides that, the corresponding model structures were established for oxidation of lattice oxygen on surface of (110) in W0.2Er0.1Ru0.7O2âδ (Supplementary Fig. 3). However, the OH* adsorbed on W0.2Er0.1Ru0.7O2âδ was not stable and transferred to Ru, W, or Er sites after optimization of the models (step II in Fig. 1c and Supplementary Fig. 3). Therefore, the lattice oxygen in W0.2Er0.1Ru0.7O2âδ will not participate in acidic OER process. This result also confirmed that introducing W and Er into RuO2 could suppress the dissolution rate of Ru.
DFT calculations were then performed to rationalize the acidic OER performance (Fig. 2d and Supplementary Figs. 4â12). As also motived by the experimental results of synthesized samples, all the models were established based on (110) surface. Besides that, the solvation was not considered in our calculation because the DFT was intended to understand trends when doping the RuO2 surface with foreign metal atoms. All potential determining steps (PDS) for these designed samples were between OOH* and O*. The calculated ÎG for PDS followed the order of W0.2Er0.1Ru0.7O2âδâ1 (0.53âeV)â<âW0.2Ru0.8O2âδâ1 (0.6âeV)â<âEr0.1Ru0.9O2âδâ1 (0.72âeV)â<âRuO2 (0.79âeV) (Fig. 2d and Supplementary Figs. 4â6), revealing that simultaneously doping W and Er could reduce energy barriers for boosting activities. This phenomenon also agreed well with the upshift of Ru 4d band centers, which tuned the adsorption energy of oxygen intermediates on active sites. Moreover, the surface of the prepared and established models for these RuO2-based materials was not bare anymore at Uâ>â1.23âV vs. RHE. And there should be neighboring intermediates around active sites, which influenced on the energetics of the elementary processes in the OER (Supplementary Figs. 9â11). According to the calculations, the energy barrier for PDS significantly decreased (Supplementary Figs. 10 and 11), compared with the traditional models without considering the neighboring intermediates (Supplementary Figs. 4â6). Therefore, these results also indicated that the neighboring intermediates around active sites could also contribute to enhancing the activity of W0.2Er0.1Ru0.7O2âδâ1 toward OER. Besides that, ÎG2 for formation of OH* has always acted as a descriptor for indicating the activity trends of the modeled materials38. The scaling relation of ÎG2 and ÎG3 for formation of OH* and OOH* should be also coupled: 2.6âeVââ¤âÎG2â+âÎG3ââ¤â3.6âeV39. As seen from the relationship of overpotential and ÎG2, the calculated ÎG2 for W0.2Er0.1Ru0.7O2âδâ1 (1.01âeV) was the apex of the volcano (Supplementary Fig. 12). Additionally, the calculated ÎG2 for W0.2Er0.1Ru0.7O2âδâ1 (1.01âeV) was also close to the theoretical ÎG2 (1.23âeV)22. These results revealed that W0.2Er0.1Ru0.7O2âδâ1 possessed the moderate ÎG2 among these structures toward acidic OER. According to above results, the co-doping of W and Er into RuO2 could effectively reduce generation of VO and soluble Rux>4, which sufficiently enhance the acidic OER stabilities and activities.
Synthesis and characterization
To experimentally probe the predictions, WmErnRu1âmânO2âδ nanosheets with various ratios were synthesized by a hydrothermal method, in which ethylenediaminetetraacetic acid (EDTA), ruthenium(III) acetate, WCl6, ErCl3, and citric acid were used as precursors (Fig. 3a). Supplementary Fig. 13 shows that these as-synthesized materials are of the rutile ruthenium oxides (RuO2, PDF#:40-1290) along with the dominant (110) peak. Transmission electron microscopy (TEM) exhibited that the prepared WmErnRu1âmânO2âδ samples possessed the ultrathin nanosheet morphology (Fig. 3b, c, and Supplementary Figs. 14â19). The corresponding elemental mapping confirmed that Er, W, and Ru were uniformly distributed in these nanosheets (Fig. 3d). High-resolution TEM (HR-TEM), fast Fourier transform (FFT), and selected area electron diffraction (SAED) revealed that W0.2Er0.1Ru0.7O2âδ nanosheets possessed a face centered cubic (FCC) crystal face along (110) direction (Fig. 3e, f). The atomic ratios of Ru, W, and Er in these as-prepared samples were analyzed by inductively coupled plasma mass spectrometry (ICP-MS). The results revealed the atomic ratio of these elements was close to the expected values (Supplementary Tables 1â6).
Oxidation state analysis
X-ray photoelectron spectra (XPS) was carried out on these prepared Ru-based nanosheet catalysts to analyze their compositions and chemical states (Supplementary Fig. 20). From Fig. 4a, the Ru 3d spectra for RuO2âδ, W0.2Ru0.8O2âδ, Er0.1Ru0.9O2âδ, and W0.2Er0.1Ru0.7O2âδ could be deconvoluted into two kinds of doublets. The binding energies of lower energy couples situated at 280.6 (Ru 3d5/2) and 284.8âeV (Ru 3d3/2) were, respectively, attributed to Rux<4 and Ru4+ 3,40. Simultaneously, the higher energy couples situated at 282.2 and 286.7âeV were, respectively, ascribed to the satellite peaks. The average area and intensity ratios for Ru 3d3/2 and Ru 3d5/2 were, respectively, higher (2:1.2) and lower (2:3.3) than the expected value of 2:3. This should be ascribed to the fact that the C 1s peak was coincided with Ru 3d41. The Ru 3d3/2 for W0.2Er0.1Ru0.7O2âδ was negatively shifted (0.15âeV) compared with those of RuO2âδ, W0.2Ru0.8O2âδ, and Er0.1Ru0.9O2âδ, demonstrating that Ru in W0.2Er0.1Ru0.7O2âδ was in a little lower valence state than Ru4+. For W 4f, the binding energies were deconvoluted into W 4f5/2 (37.4âeV) and W 4f7/2 (35.5âeV), indicating the element of W was in an valence state of W6+ (Fig. 4b)42. Similarly, the average area and intensity ratios for W 4f5/2 and W 4f7/2 doublet were, respectively, 15:17 and 9:10, which were close to the expected value of 3:4. Moreover, compared with W0.2Ru0.8O2âδ, the binding energies of W 4f5/2 and W 4f7/2 for W0.2Er0.1Ru0.7O2âδ were, respectively, negatively shifted with 0.74 and 0.60âeV (Fig. 4b), confirming that W in W0.2Er0.1Ru0.7O2âδ exhibited a lower valence state than W6+. From Fig. 3c, the binding energies for Er 4d could be fitted with two peaks, which should be, respectively, assigned to Er 4d (170.0âeV) and Er 4d5/2 (168.5âeV), revealing the Er element in these samples were in Er3+ valence state43. Interestingly, as compared with Er0.1Ru0.9O2âδ, the Er 4d5/2 and Er 4d for W0.2Er0.1Ru0.7O2âδ were, respectively, positively shifted with 0.20 and 0.47âeV (Fig. 4c). This indicates that Er3+ in W0.2Er0.1Ru0.7O2âδ was a little more positively charged. Besides that, O 1s peaks around at 532, 530.5, and 529.3âeV are, respectively, confirmed as hydroxyl groups, lattice oxygen, and MâO bonds. The lattice oxygen and MâO bonds in W0.2Er0.1Ru0.7O2âδ, W0.2Ru0.8O2âδ, and Er0.1Ru0.9O2âδ nanosheets were slightly positive shifted compared with RuO2âδ, respectively (Fig. 4d), indicating O 1s was positively charged44. Through above analyses, it could be concluded that some electrons were transferred from Er to Ru and W through O in W0.2Er0.1Ru0.7O2âδ, confirming the electronic interactions among Ru, W, Er, and O45.
Electrochemical OER in acid electrolyte
To establish correlation between the electrocatalytic performances and the specific electronic structure, the OER activities of the prepared WmErnRu1âmânO2âδ nanosheets and C-RuO2 were investigated in 0.5âM H2SO4 (Supplementary Figs. 21 and 22). The linear sweep voltammetry (LSV, scan rate: 5âmVâsâ1) was corrected with iR and normalized with geometrical area, respectively. It was seen that the performances of prepared electrocatalysts followed an order of W0.2Er0.1Ru0.7O2âδâ>âEr0.1Ru0.9O2âδâ>âW0.2Ru0.8O2âδ nanosheetsâ>âC-RuO2 NPs as shown in Fig. 6a. Simultaneously, the representative W0.2Er0.1Ru0.7O2âδ nanosheet exhibited a record performance toward OER with an overpotential (ɳ) as low as 168âmV at 10âmAâcmâ2 (Fig. 6a), remarkably preceding other prepared Ru-based nanosheets. Simultaneously, W0.2Er0.1Ru0.7O2âδ also exhibited the fastest reaction rate among these prepared samples with a small Tafel slope of 66.8âmVâdecâ1 (Fig. 6b). The current density of W0.2Er0.1Ru0.7O2âδ nanosheets could achieve 500âmAâcmâ2 at ɳâ=â275âmV, which was 28.6-folds higher than that of C-RuO2. It should be noted that the performance for the representative W0.2Er0.1Ru0.7O2âδ was also superior to those of the catalysts reported recently (Supplementary Table 9).
To reveal the origin of the performances for prepared samples, electrochemical double-layer capacitance (Cdl) was performed to measure the surface area and roughness factor (Rf) (Fig. 6c and Supplementary Figs. 23â27). Moreover, CV curves for these materials supported on glass carbon electrode were also measured to confirm Cdl (Supplementary Figs. 26 and 27). As listed in Supplementary Table 8, the calculated Cdl, Rf, and surface area exhibited a similar trend as activity. W0.2Ru0.8O2âδ, Er0.1Ru0.9O2âδ, and W0.2Er0.1Ru0.7O2âδ nanosheet structures revealed successively increased roughness compared with RuO2âδ nanosheets. The representative W0.2Er0.1Ru0.7O2âδ (1231.8âm2âgâ1) exhibited much larger surface area than Er0.1Ru0.9O2âδ (494.9âm2âgâ1), W0.2Ru0.8O2âδ (326.3âm2âgâ1), and RuO2âδ (166.7âm2âgâ1), probably originating from more active sites exposure with co-doping of W and Er into RuO2 nanosheets (Supplementary Fig. 21d). The calculated specific activity followed an order of W0.2Er0.1Ru0.7O2âδ (1.23âmAâcmâ2ox)â>âW0.2Ru0.8O2âδ (1.16âmAâcmâ2ox)â>âEr0.1Ru0.9O2âδ (1.15âmAâcmâ2ox)â>âRuO2âδ (0.98âmAâcmâ2ox)â>âC-RuO2 (0.50âmAâcmâ2ox). Simultaneously, W0.2Er0.1Ru0.7O2âδ also possessed the highest mass activity of 1518.6âAâgâ1ox, which was 2.7, 5.8, and 28.5-folds higher than Er0.1Ru0.9O2âδ, W0.2Ru0.8O2âδ, and C-RuO2, respectively (Supplementary Fig. 21d). The highest specific and mass activities for W0.2Er0.1Ru0.7O2âδ among these samples indicated the highest intrinsic activity of W0.2Er0.1Ru0.7O2âδ.
Nyquist plots in Fig. 6d were applied to fit the equivalent circuit diagram for the electrochemical impedance spectroscopy, revealing the reaction impendence between the electrode and solution. Figure 6d shows that the charge-transfer resistance (Rct) for W0.2Er0.1Ru0.7O2-δ was 3.3âΩ, which was much smaller than those of Er0.1Ru0.9O2âδ (13.8âΩ), W0.2Ru0.8O2âδ (20.2âΩ), and RuO2âδ (50.8âΩ). These results further confirmed the co-doping of W and Er into RuO2âδ nanosheets enhanced the intrinsic OER property of the samples. Additionally, Faraday efficiency for OER on W0.2Er0.1Ru0.7O2âδ could reach 99.5%, indicating the current mainly came from OER as shown in Supplementary Fig. 28.
Stability for acidic OER and PEM
The stability in acidic electrolyte is a critical factor to evaluate the performances of OER electrocatalysts due to highly oxidative operating conditions and corrosive electrolytes. The stability of W0.2Er0.1Ru0.7O2âδ could be assessed by âɳâ=âɳfinalâɳinitial, the gap between the final and initial overpotential (ɳ) can be calculated from a chronopotentiometric test of 10âmAâcmâ2. From Fig. 7a, the chronopotentometric line exhibited a highly durable stability for 250âh with a small increase of âɳâ=â83âmV, indicating the excellent stability of W0.2Er0.1Ru0.7O2âδ (Fig. 7b and Supplementary Table 9). Remarkably, from 250 to 500âh during the stability test, the operation ran stably with a tiny increase of âɳâ=â5âmV, confirming the durable properties of W0.2Er0.1Ru0.7O2âδ in acidic electrolyte. Besides that, the stability of W0.2Er0.1Ru0.7O2âδ loaded on glass carbon electrode also revealed the durable stability in acidic electrolyte (Supplementary Fig. 29). As a contrast, the C-RuO2 could only run 24âh with the large loss of âɳâ=â615âmV (Supplementary Fig. 30). To confirm the fixation of Ru in W0.2Er0.1Ru0.7O2âδ, ICP-MS was performed to analyze the concentration of the dissolved W, Er, and Ru in the solution during the 500âh stability test at an interval of 100âh (Fig. 7c). The results of ICP-MS revealed that the Ru cations in the solution were <12âppb after 500âh (Supplementary Tables 10â18). More significantly, W0.2Er0.1Ru0.7O2âδ as anodic electrocatalyst could also operate steadily for 120âh at high current of 100âmAâcmâ2 in acidic PEM (Fig. 7d, e).
Charge redistribution of W0.2Er0.1Ru0.7O2âδ for the long-durable acidic OER stability
To prove the highly improved dissolution and oxidation resistance of W0.2Er0.1Ru0.7O2âδ toward OER in acidic electrolyte, the chemical states for Ru, W, Er, and O in W0.2Er0.1Ru0.7O2âδ after the acidic OER measurement were further investigated and compared with these elements before OER, respectively (Fig. 8). The states of W0.2Er0.1Ru0.7O2âδ after the 500âh stability test were analyzed by XPS. For Ru 3d, the peaks at 284.8âeV for Ru 3d3/2 was not obviously positively or negatively shifted as compared with that of the sample before OER test, suggesting that the valence of active sites Ru remained Ru4+ during OER test (Fig. 8a). Simultaneously, the peak at 280.6âeV (Ru 3d5/2) disappeared, indicating that the Rux<4 were all oxidized to Ru4+.
Besides that, Ru valence state in W0.2Er0.1Ru0.7O2âδ after OER was also characterized with both XANES and EXAFS (Supplementary Fig. 31). W0.2Er0.1Ru0.7O2âδ showed similar adsorption energy with that for C-RuO2, indicating Ru stayed Ru4+. Simultaneously, compared with adsorption energy of W0.2Er0.1Ru0.7O2âδ before OER, adsorption energy of W0.2Er0.1Ru0.7O2âδ after OER shifted to higher energy, implying Rux<4 were all oxidized to Ru4+. This result also implied that there was no RuO4 generation during acidic OER. Moreover, EXAFS spectra exhibited that the RuâRu bond (2.71âà ) in W0.2Er0.1Ru0.7O2âδ disappeared after OER, indicating that the Rux<4 were oxidized to Ru4+ (Supplementary Table 19). Simultaneously, the RuâRu, RuâO, RuâW/Er bonds in W0.2Er0.1Ru0.7O2âδ after OER was very close to the W0.2Er0.1Ru0.7O2âδ, indicating the crystal structure was kept almost unchangeable, except for the lower Rux<4.
The binding energies for W and Er in W0.2Er0.1Ru0.7O2âδ after OER were positively shifted compared with those before OER, respectively. This suggested that W and Er in W0.2Er0.1Ru0.7O2âδ after OER were positively charged (Fig. 8b, c). Compared with O 1s in W0.2Er0.1Ru0.7O2âδ before OER, the lattice O and MâO groups in W0.2Er0.1Ru0.7O2âδ after OER were slightly positively shifted, indicating O 1s were slightly positively charged after OER (Fig. 8d). Simultaneously, the ratio between the area of the lattice oxygen and MâO groups of W0.2Er0.1Ru0.7O2âδ after OER stability became smaller, compared with that for W0.2Er0.1Ru0.7O2âδ before OER. Additionally, the hydroxyl groups disappeared in W0.2Er0.1Ru0.7O2âδ after OER, suggesting the hydroxyl groups were oxidized during OER. The surface rearrangement in W0.2Er0.1Ru0.7O2âδ during acidic OER should be responsible for the shift in these electron states (Fig. 8aâd). Simultaneously, after the surface rearrangement, partial electrons from W and Er were transferred to Ru atoms, achieving a charge redistribution among Ru, W, Er, and O. The similar redistribution was also observed in IrOx/SrIrO3 for acidic OER46. Remarkably, after electron charge redistribution, Ru 4d still keeps Ru4+ without the generation of soluble Rux+ (xâ>â4) derivatives. Therefore, charge redistribution among W, Er, Ru, and O is one reason for achieving the long-durable stability of W0.2Er0.1Ru0.7O2âδ in acidic OER.
Characterization of W0.2Er0.1Ru0.7O2âδ after acidic OER
To prove there was no generation of oxygen vacancies on W0.2Er0.1Ru0.7O2âδ during acidic OER, electron spin resonance (ESR) was further applied to characterize W0.2Er0.1Ru0.7O2âδ after the stability test (Fig. 8e). ESR spectra revealed that there were no obvious peaks around 3275âG (oxygen vacancies) in W0.2Er0.1Ru0.7O2âδ, confirming that there were no generation of oxygen vacancies during acidic OER, directly proving co-doping of W and Er could enhance the energy formation of oxygen vacancies (Supplementary Fig. 2). Simultaneously, it also proved that the OâO coupling was from H2O, rather than from lattice oxygen (Fig. 1). Consequently, there was no generation of high-valance RuO4, avoiding collapsing the crystal structure of W0.2Er0.1Ru0.7O2âδ. Additionally, the morphology of the W0.2Er0.1Ru0.7O2âδ remained intact after 500âh operation under acidic electrolyte (Fig. 8f, g). The fractional lattice spacing of W0.2Er0.1Ru0.7O2âδ also kept at ~0.356ânm, which may be due to the surface rearrangement during acidic OER (Fig. 8f and Supplementary Fig. 32). The elemental mappings of W0.2Er0.1Ru0.7O2âδ revealed that W, Er, Ru, and O were still distributed uniformly in the sample (Fig. 8g), suggesting the unchanged composition. Moreover, XRD pattern of W0.2Er0.1Ru0.7O2âδ after OER did not change obviously, implying that the structure kept unchangeable (Supplementary Fig. 33).
In summary, we applied doping engineering with W and Er for simultaneously boosting the OER activities and stabilities of RuO2 in acidic electrolyte by tuning electronic structure of RuO2. The generation of oxygen vacancies during OER was thermodynamically difficult due to the obviously downshift O 2p band centers, which improved the dissolution and oxidation resistance of Ru. Remarkably, the representative W0.2Er0.1Ru0.7O2âδ OER catalyst only required a 168âmV overpotential to reach 10âmAâcmâ2 accompanied with long-durable stability of 500âh in acidic electrolyte, overwhelmingly outperforming other Ru-based electrocatalysts reported so far. Particularly, the obtained sample could also operate stably for 120âh at a high current of 100âmAâcmâ2 in acidic PEM, further pushing its applications in industrial field. This study not only offers deep insights into tuning electronic properties of catalysts to enhance stabilities, but also opens a new-record horizon for designing electrocatalysts with super activity applying in acidic PEM.
Methods
Materials: Ruthenium(III) acetate (C14H21O14Ru3), Tungsten(VI) chloride (WCl6), and Erbium(III) chloride (ErCl3) were all directly obtained from Aladdin Reagent (Shanghai) Co., Ltd. HCl and H2SO4 was purchased from Sinopharm Co., Ltd. EDTA and citric acid were provided with Shanghai Macklin Biochemical Co., Ltd. Carbon paper (CP) was supplied with Shanghai Hesen Electric Co., Ltd. Ultra-pure (18.2âMΩâcmâ1) water was applied to deal with the prepared electrodes and electrolytes.
Preparation of RuO2âδ, W0.1Ru0.9O2âδ, W0.2Ru0.8O2âδ, W0.3Ru0.7O2âδ, Er0.1Ru0.9O2âδ, W0.2Er0.1Ru0.7O2âδ, and W0.2Er0.2Ru0.6O2âδ nanosheets
All samples were prepared by the hydrothermal method. Firstly, C14H21O14Ru3 (0.07âmol, 5.02âg), WCl6 (0.002âmol, 0.79âg), and ErCl3 (0.001âmol, 0.27âg) were added into ultra-pure water with vigorous stirring. Then, EDTA (0.01âmol, 2.92âg) and citric acid (0.01âmol, 1.92âg) was added into ultra-pure water with regulating pH to 9 by NH3·H2O. Finally, the solution containing EDTA and critic acid was added drop-wise into the metal salt solution for reaction for 12âh at 353âK. The obtained precursors were tiled in a treated porcelain boat and heated to 573âK for 4âh and 673âK for 2âh with a heating rate of 5âK/min in muffle furnace, respectively. The porous W0.2Er0.1Ru0.7O2âδ nanosheets were obtained after the muffle furnace cooling to 303âK. The methods for preparation of the RuO2âδ, W0.1Ru0.9O2âδ, W0.2Ru0.8O2âδ, W0.3Ru0.7O2âδ, Er0.1Ru0.9O2âδ, and W0.2Er0.2Ru0.6O2âδ nanosheets were the same as that of W0.2Er0.1Ru0.7O2âδ nanosheets.
Characterization
The XRD patterns for prepared electrocatalysts were all obtained on a X-pert Powder (PANalytical B.V., Netherlands) with PIXcel 1d detector and Cu-Kα (λâ=â1.54178âà ) radiation. Additionally, Tecnai G2 F20 S-TWIN (FEI, America) was used to obtain the TEM, HR-STEM, and elemental mapping at a voltage of 200âkV. For XPS analysis, X-ray photoelectron spectrometer (Thermo Scientific Escalab Xi+, England) was used to characterize these representative electrocatalysts. For ESR, JEOL FA200 was applied to characterize these samples. ICP-MS (Agilent ICP-OES730) was employed to show the ratios of these elements in W0.2Er0.1Ru0.7O2âδ. BL01C1 Beamline in NSRRC was employed to collect X-ray absorption spectra. Simultaneously, a double-crystal monochromator (Si, 111) was monochromatized for the radiation. Athena software was employed to analyze the XANES, EXAFS data.
Electrochemical measurements
These prepared electrodes were characterized by Bio-Logic VSP potentiostat. The mass-loading for prepared electrocatalysts on CP was 0.33âmgâcmâ2. The prepared C-RuO2/CP, RuO2âδ/CP, W0.1Ru0.9O2âδ/CP, W0.2Ru0.8O2âδ/CP, W0.3Ru0.7O2âδ/CP, Er0.1Ru0.9O2âδ/CP, and W0.2Er0.1Ru0.7O2âδ/CP, W0.2Er0.2Ru0.6O2âδ/CP (1âÃâ0.5âcm2) were, respectively, applied as working electrodes. Meanwhile, the platinum column electrode and Hg/Hg2SO4 (0.645âV) were, respectively, employed as counter and reference electrodes. Besides that, LSV with iR correction was applied to characterize activities of these samples in electrolytic cell (40âmL, 0.5âM H2SO4) with purging O2 about 30âmin. The EIS was tested in a range of 0.01â100âkHz at onset potential and open circuit voltage. The stabilities of C-RuO2 and W0.2Er0.1Ru0.7O2âδ were carried out at 10âmAâcmâ2 toward OER. Additionally, the catalysts loaded on new glass carbon electrodes were employed for the stability test at 10âmAâcmâ2 or 1.4âV vs. RHE.
The measurement of PEM was performed on the self-made cell, which was composed of two SS316 steel plates. In PEM, the anodic W0.2Er0.1Ru0.7O2âδ catalyst and cathodic Pt/C (20%, commercial) were, respectively, sprayed on Ti foam and CP (2âÃâ2âcm2), respectively. Additionally, the Ti foam and CP-loaded electrocatalysts were fixed with N117 membrane together using hot pressing at 363âK for 30âmin. 0.5âM H2SO4 (3âmLâminâ1) was applied as electrolyte during the test.
Double layered capacitances (Cdl) for these samples were assessed by CV with a san rate from 10 to 30âmVâsâ1 at a window of 1.01â1.11âV vs. RHE. The Cdl can be calculated from the scan rates with current densities (j) obtaining from CV curves. Simultaneously, specific capacitance of 40âμFâcmâ2 was applied to assess the ECSA. Additionally, the catalysts loaded on new glass carbon electrodes were also employed for the CV test.
Faradaic efficiency (FE)
FE for W0.2Er0.1Ru0.7O2âδ toward acidic OER (0.5âM H2SO4) was measured with the three-electrode configuration as well as detected by a GC-9790II gas chromatography. The column is 5âA molecular sieve (length: 2âm, diameter: 3âmm). The detector is thermal conductivity detector (TCD). The electrolyte under stirring was firstly degassed by Ar gas for half an hour. Then, the gaseous sample was taken out with a gas tight syringe every 20âmin at 10âmA. Additionally, the samples were analyzed via a GC calibrated O2. Each injection into GC was repeated to eliminate error.
FE could be calculated through Eq.(1):
where VO2 represents the volume of O2, F represents Faraday constant (96,485.33289âC/mol), i stands for the applying current, t is the total time for OER, and Vm represents the molar volume of the generated O2.
Theoretical calculations
Spin-polarized DFT calculations were conducted on projector-augmented wave (PAW)47 in the Vienna Ab initio Simulation Package (VASP)48. The generalized gradient approximation (GGA) of PerdewâBurkeâErnzerhof (PBE)49 exchange-functional functional was applied. The cut-off energy for plane-wave basis was set as 450âeV. A 20âà vacuum slab in a direction perpendicular to the surface of catalyst was adopted to avoid periodic interactions. The Brillouin zone integration was performed with 3âÃâ3âÃâ1 MonkhorstâPack k-point sampling for geometry relaxation50. For the calculations of DOS, the k-point mesh was increased to 6âÃâ6âÃâ1. The convergence threshold for force and energy during optimization were set as 0.03âeV/à , 10â4âeV, respectively. A p(4âÃâ4) unit cell of RuO2 (110) surface with three layers in their slab was comprised. In the maintext and supplementary information, the influence of neighboring intermediates on the energetics was not considered when free-binding energies were calculated without special instructions. Additonally, the influence of aqueous solvent in calculation were not considered. The first two layers as well as the adsorbate intermediates were relaxed with the fixation of other atoms for geometry optimization. For Er and W doping, three different initial structures are considered, +3 and +6 valences of Er and W are taken into account, respectively. It is worth mentioning that the model of Er-doped exhibited a relatively lower energy as compared with Er and W co-doped, because of the larger distortion of the adjacent O at the Er defect position, another more reasonable structural model is selected. On this basis, we simulated the adsorption behavior of *O, *OH, and *OOH intermediates for each catalyst, and each model was optimized to convergence. ÎG for each OER step was calculated through the model of computational hydrogen electrode51 along with the equation as following:
The formula for calculation of VO:
where ÎE refers to DFT energy difference; ÎS and ÎEZPE refer to corrections with entropy through vibrational analysis and zero point energy at 300âK, respectively; G(VO) refers to the energy of the structure after leaving a vacancy; G(*O) refers to the energy that adsorbs the structure of *O intermediate; GH2O and GH2 are the energy for water molecules and hydrogen molecules, respectively.
Data availability
All of the authors claim that the presented data in this work will be available for the contacts from corresponding author with reasonable request.
References
Zhang, B. et al. Homogeneously dispersed multimetal oxygen-evolving catalysts. Science 352, 333 (2016).
Sun, H. et al. Self-supported transition-metal-based electrocatalysts for hydrogen and oxygen evolution. Adv. Mater. e1806326 (2019).Â
Shan, J. et al. Charge-redistribution-enhanced nanocrystalline Ru@IrOx electrocatalysts for oxygen evolution in acidic media. Chem 5, 445â459 (2019).
Blasco-Ahicart, M., Soriano-López, J., Carbó, J. J., Poblet, J. M. & Galan-Mascaros, J. R. Polyoxometalate electrocatalysts based on earth-abundant metals for efficient water oxidation in acidic media. Nat. Chem. 10, 24 (2017).
Kasian, O., Grote, J. P., Geiger, S., Cherevko, S. & Mayrhofer, K. J. J. The common intermediates of oxygen evolution and dissolution reactions during water electrolysis on iridium. Angew. Chem. Int. Ed. 57, 2488â2491 (2018).
Geiger, S. et al. The stability number as a metric for electrocatalyst stability benchmarking. Nat. Catal. 1, 508â515 (2018).
Nong, H. N. et al. A unique oxygen ligand environment facilitates water oxidation in hole-doped IrNiOx coreâshell electrocatalysts. Nat. Catal. 1, 841â851 (2018).
King, L. A. et al. A non-precious metal hydrogen catalyst in a commercial polymer electrolyte membrane electrolyser. Nat. Nanotechnol. 14, 1071â1074 (2019).
Tobias Reier, H. N. N., Teschner, D., Schlögl, R. & Strasser, P. Electrocatalytic oxygen evolution reaction in acidic environments-reaction mechanisms and catalysts. Adv. Energy Mater. 7, 1601275 (2017).
Kibsgaard, J. & Chorkendorff, I. Considerations for the scaling-up of water splitting catalysts. Nat. Energy 4, 430â433 (2019).
Chatti, M. et al. Intrinsically stable in situ generated electrocatalyst for long-term oxidation of acidic water at up to 80â°C. Nat. Catal. 2, 457â465 (2019).
Lei Wang, Z. Z. et al. Tunable intrinsic strain in two-dimensional transition metal electrocatalysts. Science 363, 870â874 (2019).
Cao, L. et al. Dynamic oxygen adsorption on single-atomic ruthenium catalyst with high performance for acidic oxygen evolution reaction. Nat. Commun. 10, 4849 (2019).
Kim, J. et al. High-performance pyrochlore-type yttrium ruthenate electrocatalyst for oxygen evolution reaction in acidic media. J. Am. Chem. Soc. 139, 12076â12083 (2017).
Gao, J. et al. Breaking long-range order in iridium oxide by alkali ion for efficient water oxidation. J. Am. Chem. Soc. 141, 3014â3023 (2019).
Shang, C. et al. Electron correlations engineer catalytic activity of pyrochlore iridates for acidic water oxidation. Adv. Mater. 31, e1805104 (2019).
Zhao, M. et al. Ru octahedral nanocrystals with a face-centered cubic structure, {111} facets, thermal stability up to 400 degrees C, and enhanced catalytic activity. J. Am. Chem. Soc. 141, 7028â7036 (2019).
Yao, Y. et al. Engineering the electronic structure of single atom Ru sites via compressive strain boosts acidic water oxidation electrocatalysis. Nat. Catal. 2, 304â313 (2019).
Shan, J., Zheng, Y., Shi, B., Davey, K. & Qiao, S. Z. Regulating electrocatalysts via surface and interface engineering for acidic water electrooxidation. ACS Energy Lett. 4, 2719â2730 (2019).
Stoerzinger, K. A. et al. Orientation-dependent oxygen evolution on RuO2 without lattice exchange. ACS Energy Lett. 2, 876â881 (2017).
Rossmeisl, J., Qu, Z. W., Zhu, H., Kroes, G. J. & Nørskov, J. K. Electrolysis of water on oxide surfaces. J. Electroanal. Chem. 607, 83â89 (2007).
Exner, K. S. Design criteria for oxygen evolution electrocatalysts from First Principles: introduction of a unifying material-screening approach. ACS Appl. Energy Mater. 2, 7991â8001 (2019).
Retuerto, M. et al. Na-doped ruthenium perovskite electrocatalysts with improved oxygen evolution activity and durability in acidic media. Nat. Commun. 10, 2041 (2019).
Miao, X. et al. Quadruple perovskite ruthenate as a highly efficient catalyst for acidic water oxidation. Nat. Commun. 10, 3809 (2019).
Su, J. et al. Assembling ultrasmall copper-doped ruthenium oxide nanocrystals into hollow porous polyhedra: highly robust electrocatalysts for oxygen evolution in acidic media. Adv. Mater. e1801351 (2018).
Ge, R. et al. Ultrafine defective RuO2 electrocatayst integrated on carbon cloth for robust water oxidation in acidic media. Adv. Energy Mater. 9, 1901313 (2019).
Lin, Y. et al. Chromiumâruthenium oxide solid solution electrocatalyst for highly efficient oxygen evolution reaction in acidic media. Nat. Commun. 10, 162 (2019).
Grimaud, A. et al. Activating lattice oxygen redox reactions in metal oxides to catalyse oxygen evolution. Nat. Chem. 9, 457â465 (2017).
Teliska, M., OâGrady, W. E. & Ramaker, D. E. Determination of O and OH adsorption sites and coverage in situ on Pt electrodes from Pt L23 X-ray absorption spectroscopy. J. Phys. Chem. B 109, 8076â8084 (2005).
Fang, Y. H. & Liu, Z. P. Mechanism and Tafel lines of electro-oxidation of water to oxygen on RuO2 (110). J. Am. Chem. Soc. 132, 18214â18222 (2010).
Wohlfahrt-Mehrens, M. & Heitbaum, J. Oxygen evolution on Ru and RuO2 electrodes studied using isotope labelling and on-line mass spectrometry. J. Electroanal. Chem. Interfacial Electrochem. 237, 251â260 (1987).
Macounova, K., Makarova, M. & Krtil, P. Oxygen evolution on nanocrystalline RuO2 and Ru0.9Ni0.1O2âδ electrodesâDEMS approach to reaction mechanism determination. Electrochem. Commun. 11, 1865â1868 (2009).
Mefford, J. T. et al. Water electrolysis on La(1-x)Sr(x)CoO(3-delta) perovskite electrocatalysts. Nat. Commun. 7, 11053 (2016).
Grimaud, A., Hong, W. T., Shao-Horn, Y. & Tarascon, J. M. Anionic redox processes for electrochemical devices. Nat. Mater. 15, 121â126 (2016).
Zagalskaya, A. & Alexandrov, V. Role of defects in the interplay between adsorbate evolving and lattice oxygen mechanisms of the oxygen evolution reaction in RuO2 and IrO2. ACS Catal. 10, 3650â3657 (2020).
Halck, N. B., Petrykin, V., Krtil, P. & Rossmeisl, J. Beyond the volcano limitations in electrocatalysisâoxygen evolution reaction. Phys. Chem. Chem. Phys. 16, 13682â13688 (2014).
Huang, Z. F. et al. Chemical and structural origin of lattice oxygen oxidation in CoâZn oxyhydroxide oxygen evolution electrocatalysts. Nat. Energy 4, 329â338 (2019).
Man, I. C. et al. Universality in oxygen evolution electrocatalysis on oxide surfaces. ChemCatChem 3, 1159â1165 (2011).
Exner, K. S. Universality in oxygen evolution electrocatalysis: high-throughput screening and a priori determination of the rate-determining reaction step. ChemCatChem 12, 2000â2003 (2020).
Shen, J. Y., Adnot, A. & Kaliaguine, S. An ESCA study of the interaction of oxygen with the surface of ruthenium. Appl. Surf. Sci. 51, 47â60 (1991).
Petrykin, V. et al. Local structure of nanocrystalline Ru1âxNixO2âδ dioxide and its implications for electrocatalytic behaviorâan XPS and XAS study. J. Phys. Chem. C 113, 21657â21666 (2009).
Kumari, S. et al. A low-noble-metal W1âxIrxO3âδ water oxidation electrocatalyst for acidic media via rapid plasma synthesis. Energy Environ. Sci. 10, 2432â2440 (2017).
ReszczyÅska, J. et al. Visible light activity of rare earth metal doped (Er3+, Yb3+ or Er3+/Yb3+) titania photocatalysts. Appl. Catal. BâEnviron. 163, 40â49 (2015).
Shan, J., Ling, T., Davey, K., Zheng, Y. & Qiao, S. Z. Transition-metal-doped RuIr bifunctional nanocrystals for overall water splitting in acidic environments. Adv. Mater. 31, e1900510 (2019).
Jiang, J. et al. Atomic-level insight into super-efficient electrocatalytic oxygen evolution on iron and vanadium co-doped nickel (oxy)hydroxide. Nat. Commun. 9, 2885 (2018).
Seitz, L. C. et al. A highly active and stable IrOx/SrIrO3 catalyst for the oxygen evolution reaction. Science 353, 1011â1014 (2016).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953â17979 (1994).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169â11186 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865â3868 (1996).
Weiler, M., Aggarwal, R. & Lax, B. Interband magnetoreflectance in semiconducting Hg1-xCdxTe alloys. Phys. Rev. B 16, 3603â3607 (1977).
Nørskov, J. et al. Origin of the overpotential for oxygen reduction at a fuel-cell cathode. J. Phys. Chem. B 108, 17886â17892 (2004).
Acknowledgements
Zhejiang Provincial Natural Science Foundation of China (Grant No. LR17B060003), Natural Science Foundation of China (Project Nos. 21872174, U1932148, 21776248, and 21676246), and International Science and Technology Cooperation Program (Grant Nos. 2017YFE0127800 and 2018YFE0203400) supported this work. This study was also supported by the Fundamental Research Funds for the Central Universities (2020XZZX002-07).
Author information
Authors and Affiliations
Contributions
S.Y.H. synthesized the catalysts. S.Y.H. analyzed the data and wrote the manuscript. M.L. and X.W.Z. revised the manuscript. J.J.P., X.N.L., X.L.T., Y.H. and N.X. conducted the theoretical calculations. L.C.L. helped revise the manuscript. All the authors commented the results and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Qianwang Chen, Kai Exner, and Maria Retuerto for their contributions to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hao, S., Liu, M., Pan, J. et al. Dopants fixation of Ruthenium for boosting acidic oxygen evolution stability and activity. Nat Commun 11, 5368 (2020). https://doi.org/10.1038/s41467-020-19212-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-19212-y