Yu-Xin
Luan‡
,
De-Kun
Liu‡
,
Yu-Xiang
Chen
,
Xue-Zhi
Song
*,
Jing
Liang
*,
Xiao-Bing
Wang
,
Wen-Qiang
Sun
,
Yu-Lan
Meng
,
Chen-Yu
Shen
,
Ruichen
Zhou
,
Jinxuan
Liu
,
Xiao-Feng
Wang
* and
Zhenquan
Tan
*
Leicester International Institute, School of Chemical Engineering, Ocean and Life Sciences, School of General Education, State Key Laboratory of Fine Chemicals, Dalian University of Technology, 2 Dagong Road, Liaodongwan New District, Panjin, Liaoning 124221, China. E-mail: songxz@dlut.edu.cn; jingliang@mail.dlut.edu.cn; wangxf@dlut.edu.cn; tanzq@dlut.edu.cn
First published on 4th July 2025
Developing highly efficient and durable electrocatalysts for the oxygen evolution reaction (OER) holds great promise in revolutionizing the sustainable energy-conversion technologies, which greatly rely on the surface electronic structures. However, constructing catalytically active medium-entropy materials encounters great challenges due to the complex composition and the unsolved electronic structure–performance relationship. In this study, medium-entropy metal selenide (CoFeNiSe) with a cumulus-like architecture is fabricated through the selenylation of trimetallic hydroxide precursors. The medium-entropy CoFeNiSe manifests outstanding electrocatalytic OER activity and kinetics, that is, a low overpotential of 268 mV at 10 mA cm−2 and a Tafel slope of 53.33 mV dec−1 in 1 M KOH, outperforming the entropy-poor binary CoFeSe and CoNiSe, as well as its unary CoSe counterparts. It has been found experimentally and theoretically that the effective modulations of the valency and the d-band center at an optimal level regulate the adsorption/desorption ability of reaction intermediates via multimetallic electron interactions. This work provides valuable insight into the electronic structure regulation in entropy-rich materials, shedding light on the development of current green energy conversion technologies.
Numerous studies have witnessed that the design and preparation of Earth-abundant 3d-metal-based materials via controlling morphology and varying compositions, including metal hydroxides, metal phosphides, metal–organic frameworks, nitrides, etc., can regulate their activity towards the OER.7–15 Among these candidates, metal selenides show great superiority to other metal chalcogenides (e.g., sulfides and oxides) in improving the catalytic activity for the OER due to the higher intrinsic electronic conductivity.16–18 Researchers have developed various strategies to modulate the structures for further enhancing electrocatalytic performance. For example, seminal investigations demonstrated that the modulation of the metal to Se ratio in NiSex and CoSex would affect catalytic activity.19,20 However, mono-metallic selenides face the dilemma of simple structures, untunable morphologies and limited ability to adjust the electronic structure. These issues can be overcome to some degree by constructing bimetallic selenides or heterostructures, such as Fe doping in NiSe or MoSe2@NiCo2Se4 because of the synergistic effect between different metallic constituents.21,22 Notwithstanding these achievements, it is still challenging to explore novel metal selenide-based electrocatalysts by exquisitely controlling the composition to engineer binding energy with intermediates accurately for efficient and durable electrocatalysis.
Since their discovery in 2004, high-entropy materials have emerged into a hotspot of electrocatalysis and shown great potential propelled by their core effects, that is, the high-entropy effect, the lattice distortion effect, the slow diffusion effect, and the cocktail effect.23,24 These features and broad metal elements endow high-entropy materials with distorted lattice structures and excellent stability. Specifically, the random arrangement of multi-element atoms and the enhanced d–d interaction allow high-entropy materials for fine-tuning the electronic structure of active sites and desirable binding energy, providing a new solution to active electrocatalysis.25–28 Under such circumstances, medium-entropy materials (MEMs) constituting four or fewer metallic components also have the above-mentioned merits, and MEMs will deliver enhanced catalytic activity when the size comes to the nanoscale.29 Moreover, MEMs have great ease in optimizing the metal component and recovering metal species in comparison with their high-entropy counterparts, which is beneficial for industrial applications. Previous studies have introduced high entropy or medium entropy into metal selenides, including an FeCoMnNiSe2 electrocatalyst on Ni foam, (Ni,Fe,Co)Se2 nanosheets on a stainless steel mesh substrate and Ni–Co–Fe–Se@NiCo-LDH nanoarrays on Ni foam.30–33 However, the exploration of advanced unsupported medium-entropy metal selenides comes across two obstacles: the miscibility of multi-elements partially inhibiting multiplicity synergies and the difficulty in deciphering the relationship between the electronic structure and activity. As a consequence, it is highly desirable to develop powerful methods to uniformly distribute metal components and investigate their electronic structure in depth.
In this work, one easily accessible approach was developed to fabricate homogeneous tri-metallic precursors. Then, the selenylation step was conducted under hydrothermal conditions to synthesize the medium-entropy CoFeNi selenide (CoNiFeSe). Significantly, the incorporation of Co, Ni and Fe elements adjusted the d-band center's position, thereby enhancing the catalytic OER activity of the newly designed medium-entropy CoNiFeSe material with an impressively low overpotential of 268 mV at 10 mA cm−2. This work provides some fundamental insight into the electronic structure regulation for the electrocatalysts applied in water electrolysis.
In the second step, for solution C, 100 mg of NaBH4 was dissolved in 10 mL of deionized water. Solution C was slowly added dropwise into the mixed solution and continuously stirred for 30 min. Then the final suspension was centrifuged and the products were washed with deionized water twice and anhydrous ethanol once and then dried at 60 °C.
For the preparation of comparable CoNi-LDH, CoFe-LDH, and Co-LDH materials, the initial solutes in solution A were different from those in CoNiFe-LDH, which were Co(NO3)2·6H2O (1.00 mmol) and Ni(NO3)2·6H2O (1.00 mmol) for CoNi-LDH, Co(NO3)2·6H2O (1.00 mmol) and Fe(NO3)3·9H2O (1.00 mmol) for CoFe-LDH, and Co(NO3)2·6H2O (2.00 mmol) for Co-LDH. The other procedures are identical to the preparation of CoNiFe-LDH.
The morphologies of the as-synthesized LDH-based precursors and metal selenide materials were characterized using SEM. The as-synthesized CoNiFe-LDH displayed uniform sakura-like morphology with abundant nanosheets as petals assembling together with a diameter of about 100–200 nm (Fig. 1a and b). Moreover, the metallic compositions in different LDH-based precursors also exert a profound effect on their morphologies, in which the morphology of CoNi-LDH (Fig. S1a and b†) is similar to that of CoNiFe-LDH. For Co-LDH, the lateral size becomes much larger, which is similar to the surface of a cocoon (Fig. S2a and b†). But the CoFe-LDH shows coralline-like morphology (Fig. S3a and b†). After the selenization treatment of CoNiFe-LDH, CoNiFeSe basically retains the initial state of CoNiFe-LDH and further reveals a cumulus-like architecture (Fig. 1c and d). The cumulus-like architecture features uniform particle morphology with visible interconnections originating from the flowerlike CoNiFe-LDH. The TEM image in Fig. 1e reveals that CoNiFeSe exhibits a clustered structure formed by ultrathin nanosheets as building blocks, leaving some pores. These special morphological features, such as nanosheets and pores, would expose more active catalytic sites on the material's surface and enable better electrolyte diffusion, faster charge transfer, and efficient removal of gases generated during the catalytic reaction, ultimately contributing to superior OER catalytic performance. The high-resolution TEM (HRTEM) image in Fig. 1f displays lattice spacings of 0.203 nm, corresponding to the (102) plane of FeCoNiSe (PDF#52-1008). The SAED image showing a set of diffraction rings suggests a polycrystalline diffraction pattern (Fig. S4†).36 The prominent diffraction rings observed in the SAED pattern can be assigned to the (101), (103), and (202) planes of Co0.85Se (JCPDS: 52-1008), supporting the successful formation of CoNiFeSe with homogeneous mixed metallic elements. Furthermore, high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images and EDS elemental mapping images illustrate a uniform distribution of Fe, Co, Ni and Se throughout the CoNiFeSe structure (Fig. 1g). The uniform distribution with an almost equal atomic content of Co, Ni and Fe (Fig. S5†) fosters a synergistic interaction among the multiple elements to regulate the electronic structure, further boosting conductivity and improving oxygen adsorption/desorption capabilities, both of which positively contribute to the electrochemical catalytic performance.37 According to the EDS results, the calculated ΔSmix of the CoNiFeSe material is 1.08R, which confirms its medium-entropy characteristic.
To further confirm the crystallographic feature of the as-synthesized metal selenides, powder XRD measurement was performed and the corresponding patterns are compiled in Fig. 2a. The diffractions of CoNiFeSe, CoNiSe, CoFeSe and CoSe basically match the simulated ones in Co0.85Se (JCPDS 52-1008). In the case of CoNiFeSe, it is clear that the reflections at 2θ being 33.26°, 44.74°, and 50.56° were attributed to the (101), (102), and (110) crystal planes of Co0.85Se (JCPDS 52-1008), which agrees with the SAED result above. In comparison with unary CoSe, a slight negative shift to the low diffraction angle direction is observed in binary CoFeSe. This shift is likely caused by the replacement of Co ions with larger Fe ions, increasing the lattice constant and decreasing the 2θ value according to the Bragg equation (2dsin
θ = nλ). The Ni-containing materials, such as CoNiSe and CoNiFeSe, show enhanced crystallinity, which may illustrate Ni doping significantly endowing the strong synergistic effect among different centers in ternary metallic selenides.38
![]() | ||
Fig. 2 (a) XRD of CoNiFeSe, CoFeSe, CoNiSe and CoSe, (b) full XPS spectrum, (c) Co 2p spectrum, (d) Ni 2p spectrum, (e) Fe 2p spectrum, and (f) Se 3d spectrum of CoNiFeSe. |
The full XPS spectrum of CoNiFeSe shown in Fig. 2b shows the presence of Ni 2p, Co 2p, Fe 2p, O 1s and Se 3d. The Co 2p XPS spectrum reveals four distinct 2p3/2 peaks and four 2p1/2 peaks (Fig. 2c). These peaks in the 2p3/2 region include the peak at 778.75 eV which can be attributed to the Co–Se bond, the peak at 781.45 eV which can be assigned to high-valence Co3+, the intermediate state Co2+and the satellite peaks located at 785.45 eV and 788.52 eV. There are also four corresponding peaks in the 2p1/2 region.39,40 The Ni 2p3/2 spectrum of Fig. 2d can be deconvoluted into two components: one at 853.5 eV, accompanied by a satellite peak, corresponding to the Ni–Se bond in CoNiFeSe, and the other at 856.5 eV, with a satellite peak being attributed to Ni2+.41 In the Fe 2p region of CoNiFeSe, as depicted in Fig. 2e, the binding energy at 707.20 eV is assigned to the Fe–Se bond, aligning with previous studies.38 Additionally, two peaks at 711.52 eV and 724.80 eV, along with their satellite peaks, correspond to Fe 2p3/2 and Fe 2p1/2 of the Fe–O bond, respectively.38,42 The high-resolution Se 3d XPS spectrum can be deconvoluted into three distinct peaks (Fig. 2f). The two primary peaks at binding energies of 53.85 eV and 55.35 eV correspond to the Se 3d5/2 and Se 3d3/2 spin–orbit components, respectively. These binding energies are characteristic of selenium in the −2 oxidation state, indicating that the selenium in the CoNiFeSe compound exists predominantly as Se2−. Additionally, the broad peak observed at 59.40 eV can be attributed to the presence of the oxidized Se species, indicating that the surface of selenides may be oxidized in the air.43,44
The electrocatalytic performance of the OER of the selenides was studied by linear sweep voltammetry (LSV) under the conditions of a typical three-electrode system using 1.0 M KOH as the electrolyte. Fig. 3a shows the LSV curves of CoNiFeSe, CoFeSe, CoNiSe and CoSe, respectively. The CoNiFeSe material can drive an earliest current response at 1.33 V (vs. RHE), and the current will dramatically increase when the anodic applied potential is increased. The overpotentials to attain a density of 10 mA cm−2 were 268 mV for CoNiFeSe, 307 mV for CoFeSe, 356 mV for CoNiSe and 338 mV for CoSe, respectively (Fig. 3b). Clearly, the overpotential of CoNiFeSe is much smaller than others, which proves that this medium-entropy catalyst has excellent electrocatalytic activity. OER kinetics was studied by extracting the Tafel slope from the LSV curve, since the value of the Tafel slope is a key parameter to reveal the OER kinetics. As illustrated in Fig. 3c, the Tafel slope value of CoNiFeSe is 53.33 mV dec−1, which is smaller than CoFeSe (57.04 mV dec−1), CoNiSe (89.26 mV dec−1) and CoSe (88.55 mV dec−1), revealing the favorable kinetics accelerated by trimetallic components. We also evaluated the impressive CoNiFeSe material against other reported materials, including but not limited to high/medium-entropy materials and metal selenides, based on overpotential at a current of 10 mA cm−2 and the Tafel slope (Table S1†). Notably, its performance either matches or exceeds that of many reported materials, such as Ni(CN)2/NiSe2, (NiFeCoMn)3S4 and Ni@Co/Se2 or NiFe-Se/CFP.45–48
The electrochemically active surface area (ECSA) indicates the density of electrochemically active sites present per unit area of the catalyst material. The double-layer capacitance (Cdl), derived from cyclic voltammetry (CV) curves at varying scan rates, is directly proportional to the ECSA.49 The magnitude of Cdl serves as an indicator of the quantity of electrochemically active sites available. CVs were collected at sweep rates of 20 mV s−1, 40 mV s−1, 60 mV s−1, 80 mV s−1, 100 mV s−1 and 120 mV s−1 for 40 consecutive cycles and are shown in Fig. S6.† The Cdl is calculated by the capacitive current difference (Δj) at different scan rates in the non-faradaic region. The Cdl of CoNiFeSe is 1.70 mF cm−2, CoFeSe is 1.34 mF cm−2, CoNiSe is 4.32 mF cm−2 and CoSe is 12.93 mF cm−2, revealing that uni-metallic CoSe possesses a large number of active sites (Fig. 3d). And the decreased active site number may be related to the micromorphology change caused by the introduction of one or two extra metallic components. The distinguished electrocatalytic activity and the limited active site of CoNiFeSe provide a hint that other factors may determine the final activity. To further understand the OER performance of the polymetallic selenides, electrochemical impedance spectroscopy (EIS) was conducted to analyze the charge transfer behavior at the interface between the electrolyte and the catalyst surface, covering a wide frequency range from 100000 Hz to 0.01 Hz. The radius of curvature in this context is directly related to the charge transfer resistance and inversely related to the charge transfer capacity. The Nyquist plots recorded for the CoNiFeSe, CoFeSe, CoNiSe and CoSe catalysts are shown in Fig. 3e with an equivalent circuit. The CoNiFeSe electrode has an Rct value of 30.52 Ω, smaller than CoSe (110.70 Ω), CoFeSe (92.23 Ω), and CoNiSe (202.50 Ω). It proves that CoNiFeSe has the lowest charge transfer resistance, which is in good agreement with the LSV trend above. The exceptional durability of the polymetallic catalyst is essential for sustainable energy-related electrocatalytic applications. The CoNiFeSe catalyst remained stable for 96 hours (Fig. 3f)with a good retention of the initial current density of 10 mA cm−2, which underscores its excellent stability. After the electrocatalytic stability test, the recorded XPS spectra of Co, Ni, and Fe changed significantly in comparison with their initial states (Fig. S9–11†). After the stability test, the Co 2p XPS spectrum showed the disappearance of Co–Se bonds and the formation of CoOOH, which may serve as the primary active phase.44,50 The HRTEM image of CoNiFeSe after the stability test reveals the formation of amorphous CoOOH on the outer layer with the retention of metal selenides inside (Fig. S12†). The Raman spectrum of CoNiFeSe after the electrochemical stability test also indicates the formation of an active MOOH phase (Fig. S13†).
To further elucidate the mechanism of the activity enhancement in medium-entropy CoNiFeSe, a set of analyses and investigations were conducted. The ECSA-normalized current curves presented in Fig. S14† clearly demonstrate the much larger current density of CoNiFeSe than other counterparts, which reveals that the medium-entropy CoNiFeSe has the largest intrinsic activity at each active site. To the best of our knowledge, the intrinsic per-site activity is relevant to the electronic structure of the active site. Then, the XPS spectra of different materials have been compiled and compared in Fig. 4a to clarify the electronic structure modulated effect. The four materials show similar deconvoluted peaks, with Co–Se peaks, Co3+/Co2+ peaks and satellite peaks. The area ratio of the Co–Se bond basically increases from unary CoSe to the target ternary CoNiFeSe material, indicating the stabilization effect on Co–Se species induced by the entropy increase (Fig. 4b). And the content of the easily oxidized species including Co–Se and Co2+ may endow the material with more active sites, possibly improving the activity for the OER.51,52
Based on the well-established d-band theory, the d-band center (Ed) serves as a valuable descriptor for examining the interactions between metal sites on the catalyst surface and oxygen-related intermediates.53 On the basis of valence band spectrum analysis among the materials studied (Fig. 4c), the Ed energy levels were calculated. The introduction of individual Ni species would downshift the d-band from CoSe (4.95 eV) to CoNiSe (5.05 eV). This d-band shift would lead to weak bonding to the reaction intermediates, hindering their adsorption on the CoNiSe and CoSe. Conversely, the d-band manifests an upward shift when the single Fe element is introduced to CoSe, with a d-band center of 4.46 eV in CoFeSe. The elevated Ed energy level relative to the Fermi level raises the antibonding states and reduces the occupancy of the antibonding state, resulting in overly tight bonding between the metal sites and intermediates.54 The medium-entropy CoNiFeSe material possesses an optimized Ed energy level (4.62 eV), which allows it to effectively balance the adsorption and desorption processes of OH*, O*, and OOH* intermediates on its metal sites, thereby improving its electrocatalytic OER performance.55,56 The electronic structure modulation is directly related to the charge transfer and electronic interactions, which can be elucidated in terms of the valence electron structure on the basis of XPS results. In this context, the M–Se–M bonding (where M = Co or Fe) is used to describe the electron interactions between transition metal elements. As illustrated in Fig. 4d, the unpaired electron in the π-symmetry (t2g) d-orbitals of Co and Fe cations can only achieve a weak electron interplay via an Se bridge with π donation features. The two Ni2+ (t2g6eg2) ions have fully occupied π-symmetry (t2g) d-orbitals, leading to strong electron–electron repulsion between Se2− and Ni2+ ions.57,58 After coupling Co, Ni and Fe cations, the simultaneous e−–e− repulsion and π-donation would trigger the modulation of the electronic structure of the metal center, being conducive to efficient OER elelctrocatalysis.59
To further explore the inherent electronic effect in medium-entropy CoNiFeSe, density functional theory (DFT) calculations were performed in this work. The structure of medium-entropy CoNiFeSe (102) was constructed for calculation in Fig. 4e, while CoNiSe (102) and CoFeSe (102) structures are shown in Fig. S16 and S17.† The projected density of states (PDOS) for different materials was calculated and is compiled in Fig. 4f. Then, the d-band centers of CoNiFeSe, CoNiSe and CoFeSe are determined to be −1.271 eV, −1.344 eV and −1.163 eV, respectively. The theoretical d-band center order is in good agreement with the experimental XPS valence band spectrum analysis. These results further demonstrate that the medium-entropy CoNiFeSe material exhibits moderate intermediate adsorption strength by the synergistic Ni and Fe electronic modulated effect, with the adsorption strength order of CoFeSe > CoNiFeSe > CoNiSe.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5dt00739a |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2025 |