
==== Front
RSC Adv
RSC Adv
RA
RSCACL
RSC Advances
2046-2069
The Royal Society of Chemistry

d4ra04011b
10.1039/d4ra04011b
Chemistry
Computational screening of transition metal atom doped ZnS and ZnSe nanostructures as promising bifunctional oxygen electrocatalysts†
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra04011b

https://orcid.org/0000-0001-9621-4502
Xia Feifei a
Shu Li a
Yang Fengli a
Wen Yingpin a
Zheng Chunzhi a
a School of Chemistry and Chemical Engineering, Jiangsu University of Technology Changzhou 213001 Jiangsu P. R. China ffxia@jsut.edu.cn

12 9 2024
12 9 2024
12 9 2024
14 40 2899829005
31 5 2024
30 8 2024
This journal is © The Royal Society of Chemistry
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by-nc/3.0/ This article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given and it is not used for commercial purposes.
The design of bifunctional oxygen electrocatalysts showing high catalytic performance for the oxygen evolution reaction (OER) and oxygen reduction reaction (ORR) is of great significance for developing new renewable energy storage and conversion technologies. Herein, based on the first principles calculations, we systematically explored the electrocatalytic activity of a series of transition metal atom (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe nanostructures for OER and ORR. The calculated results revealed that Ni- and Pt-doped ZnS and ZnSe nanostructures exhibit promising electrocatalytic performance for both OER and ORR in comparison to the pristine ZnS and ZnSe nanostructures. Especially, the OER/ORR overpotentials of Ni-doped ZnS and ZnSe nanostructures are estimated to be 0.28/0.30 and 0.31/0.31 V, respectively, disclosing their great potential as bifunctional oxygen electrocatalysts. Moreover, it is found that Ni-doped ZnS and ZnSe nanostructures for OER and ORR are on the top of the volcano plots, evincing promising catalytic performance. Our results provide theoretical insights into a feasible strategy to synthesize highly efficient ZnS- and ZnSe-based bifunctional oxygen electrocatalysts in the future.

The computed OER and ORR overpotentials of pristine and TM doped ZnX (X = S and Se) nanostructures.

National Natural Science Foundation of China 10.13039/501100001809 21703087 pubstatusPaginated Article
==== Body
pmc1. Introduction

Developing sustainable and green energy technologies is extremely urgent to solve the growing energy crisis and environmental pollution, which can ensure a safe and sustainable future.1–5 As two critical processes in the electrochemical energy conversion, oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) involve a four-electron process with complicated mechanisms and exhibit sluggish kinetics and high overpotentials,6–8 which requires highly active catalysts for the practical applications. Currently, RuO2 and IrO2 are popular as outstanding OER electrocatalysts,9–11 and noble metals such as Pt-based materials are the best ORR electrocatalysts.12,13 However, the high cost and low selectivity of these electrocatalysts considerably impede their large-scale commercial applications. Therefore, it is urgent to develop low-cost and high-selectivity OER/ORR electrocatalysts with promising catalytic performance.

It is reported that transition metal chalcogenides (TMCs) containing sulfide and selenide have considerable electrical conductivity and tunable electronic configurations, which is beneficial to their wide application in electrocatalysis.14–22 Moreover, these TMCs can draw considerable attention from researchers due to their natural abundance and intrinsic activity.17 For instance, Garaj et al.18 demonstrated that the Gibbs free energy of hydrogen adsorption on the edges of MoS2 is close to zero, leading to the high catalytic activity for hydrogen evolution reaction (HER). Zheng et al.19 explored that nickel selenide with different stoichiometric ratios (NiSe, Ni0.5Se, Ni0.75Se and Ni0.85Se) exhibit HER and OER catalytic activity in an alkaline medium. Moreover, the synthesized NiS and FeSe materials20 exhibited excellent catalytic activity with overpotentials of 83 and 171 mV for HER and OER, respectively. All these reported results show that TMCs exhibit the potential catalytic performance and can serve as electrocatalysts for HER and OER.

TMCs, ZnS and ZnSe have become the electrode materials for HER, OER, and ORR because of their natural abundance, higher electrical conductivity and tunable electronic configurations.23–28 For example, Chen et al.23 presented that Ni and Sn anchored on the surface of ZnS materials can obviously improve their HER catalytic activity, while Sn and Pt supported on ZnS surface display higher OER catalytic activity. Yang et al.25 demonstrated that the ORR catalytic performance of ZnSe was ameliorated by the introduction of C and N atoms, which tunes the electronic properties of ZnSe and leads to the lower reaction energy barrier for the *O to *OOH conversion. These works suggest that metal or non-metal atoms can modulate the electronic properties of ZnS and ZnSe, resulting in the improvement of catalytic performance for HER, OER and ORR. Moreover, elemental doping is a facile and effective method to improve the intrinsic activity of catalysts during the reaction process.29 The introduction of metal or non-metal atoms can not only agitate the crystal lattice to expose more catalytic centers but also regulate their intrinsic electronic structure, which can activate the inherent catalytic activity.29,30 Motivated by these advantages, it is anticipated that transition metal (TM) single atoms (Fe, Co, Ni, Cu, Pd and Pt) doping can tune the electronic properties of ZnS and ZnSe nanostructures, which is beneficial to their catalytic activity.

Herein, we systemically explored the OER and ORR catalytic activity of pristine and TM (Fe, Co, Ni, Cu, Pd and Pt) doped ZnS and ZnSe nanostructures using the first principles calculations based on density functional theory (DFT).31,32 Firstly, the stabilities and catalytic active sites of TM (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe nanostructures were investigated. Subsequently, the OER/ORR catalytic activity of TM (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe nanostructures were examined. Then, the origin of OER/ORR catalytic activity was analyzed. Among all the Fe, Co, Ni, Cu, Pd and Pt-doped ZnS and ZnSe nanostructures, Ni and Pt doping obviously improves the OER/ORR catalytic activity of ZnS and ZnSe nanostructures. Interestingly, Ni-doped ZnS and ZnSe nanostructures possess the lowest OER/ORR overpotentials (0.28/0.30 and 0.31/0.31 V), suggesting their promising bifunctional catalytic performance. This work not only provides helpful theoretical guidance for developing inexpensive and effective OER/ORR electrocatalysts but also extends the possible application of ZnS- and ZnSe-based materials in water electrocatalysis fields.

2. Computational methods

All the geometry optimization and energy calculations were undertaken using the Vienna ab initio simulation package (VASP)33,34 based on the DFT. The electron exchange and correlation effects were handled with the Perdew–Burke–Ernzerhof (PBE) function of the generalized gradient approximation (GGA).35–37 Generally, the DFT + U method is more accurate than GGA density functionals to describe the partially filled d and f electronic states of transition metal elements.38 It has been reported that the effects of spin polarization on the free energies of each elementary step for OER/ORR can be ignored.39,40 However, our aim in this work is to screen the ideal OER/ORR catalysts from various ZnS and ZnSe candidates, and we are particularly interested in the catalytic tendency of these candidates for OER/ORR. Therefore, the spin polarization method was not adopted in this study. The van der Waals interaction was considered using Grimme's scheme (DFT-D3).41 The cutoff energy for the plane-wave basis was set as 500 eV. The energy and force convergence criteria were set to 10−6 eV and 0.01 eV Å−1, respectively. The vacuum space of 18 Å was selected to avoid the interaction between the periodical images. We performed our calculations using a 3 × 3 × 1 supercell of TM (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe nanostructures. Brillouin zone sampling is employed with Monkhorst–Pack special k-point meshes,42 and the 6 × 6 × 1 k-grid was chosen for the structural relaxation and frequency calculation.

To evaluate the stability of the TM (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe nanostructures, the binding energy (Eb) was calculated using the following equation:1 Eb = ETM–ZnX − ETM − Evacancy-ZnX

where ETM–ZnX, ETM and Evacancy-ZnX are the calculated DFT energies of TM doped substrate, isolated TM atom and Zn-vacancy ZnX (X = S and Se) nanostructure, respectively.

Free energies of OER/ORR intermediates in electrochemical reaction pathways were evaluated according to the computational hydrogen electrode (CHE) model developed by Nørskov et al.43,44 using the reversible hydrogen electrode (RHE) as a reference whose chemical potential (H+ + e−) is equivalent to that of 1/2H2 at 0 V and at all pH values. The Gibbs free energy change (ΔG) of each elementary step during the OER/ORR is defined as follows:2 ΔG = ΔE + ΔEZPE − TΔS + ΔGU + ΔGpH

in which ΔE is the total energy obtained by DFT calculations. ΔEZPE and ΔS are the zero-point energy corrections and entropy contributions at 298.15 K, respectively, deriving from the computed vibrational frequencies for the adsorbed *OH, *O and *OOH intermediates. ΔGU and ΔGpH are the applied electrode potential and the free energy correction of the pH (pH = 0 in this work), respectively. To avoid the inaccurate calculation of the free energy of O2 gas, the experimental reaction energy of 2H2O → O2 + 2H2 (4.92 eV) is involved. The vibrational frequencies of *OH, *O and *OOH intermediates were analyzed for OER/ORR using VASPKIT software, which is a pre and post-processing program.45

3. Results and discussion

3.1 Stabilities and catalytic active sites

The stability of ZnS and ZnSe surface structures was considered by the phonon spectra before their calculation of the electrocatalytic activity. To verify the dynamical stability of the ZnS and ZnSe surface structures as well as whether their structures are optimized, the phonon spectra were calculated based on the same DFT level and are shown in Fig. S1.† Note that there is no imaginary frequency in both the ZnS and ZnSe systems, which confirms their kinetic stability and well-optimized structures. To tune the catalytic performance of ZnS and ZnSe nanostructures, the substitutional doping method was employed. As shown in Fig. 1a, TM (TM = Fe, Co, Ni, Cu, Pd and Pt) atom substitutes the Zn atom in the first surface of the ZnS nanostructure and the similar substitution for ZnSe nanostructures (Fig. S2a†). The stable combination of TM atom as the active site on substrate materials is the premise for catalysts to maintain high catalytic activity for a long time.46 However, owing to the accumulation of a large number of positive charges, the TM atom is usually considered as the active site on substrate materials for OER/ORR,47 which is beneficial to form *OH, *O and *OOH intermediates. Therefore, Fe, Co, Ni, Cu, Pd and Pt sites were considered as the catalytic active site on the surface of the doped ZnS and ZnSe nanostructures during the OER/ORR process. Moreover, the binding energy (Eb) of the TM-doped ZnS and ZnSe nanostructures were calculated to investigate the stabilities of these materials (shown in Fig. 1b and S2b†). Results show that the Eb of all the TM atoms to the ZnS and ZnSe nanostructures is negative, thus demonstrating the better stability of these nanostructures, which can be considered as the substrate materials for catalysts.

Fig. 1 (a) Top and side views of the TM (TM = Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS nanostructures, and the pink, gray and yellow atoms represent TM, Zn and S. (b) The binding energy of Fe, Co, Ni, Cu, Pd and Pt-doped ZnS nanostructures.

3.2 OER/ORR catalytic activity

After confirming the stabilities and catalytic active sites of TM-doped ZnS and ZnSe nanostructures, the electrocatalytic activity of these materials for OER/ORR was systematically investigated. It is reported that the OER Gibbs free energy changes along the four-electron path for the previously screened systems.48,49 The computational results23–25,50 suggested that the four-electron step occurs on the surface of ZnS, ZnSe and CdS materials during the OER process, and TM atoms can improve their catalytic performance. Thus, the four elementary electrons for OER/ORR were adopted on the surface ZnS and ZnSe nanostructures in this work, as shown in Fig. 2, and the important *OH, *O, and *OOH intermediates are involved. Accordingly, the detailed reaction pathways, free energy equations and overpotentials for OER and ORR under an acidic condition are presented in Table 1. Moreover, free energy changes of each elementary step and overpotentials of ZnX (X = S and Se) nanostructures for OER and ORR at U = 0 V and pH = 0 are shown in Tables S1 and S2.†

Fig. 2 Schematic of the mechanism of OER/ORR on the surface of TM-doped ZnS and ZnSe nanostructures. The brown, gray, yellow, red and white atoms are TM, Zn, S, O and H.

The reaction pathways, free energies equations and overpotentials for OER and ORR processesa

Reaction	Reaction pathway	ΔG/eV	η/V	
OER	H2O(1) + * → *OH + H+ + e−	ΔG1 = ΔG*OH		
*OH → *O + H+ + e−	ΔG2 = ΔG*O − ΔG*OH	
H2O(1) + *O → *OOH + H+ + e−	ΔG3 = ΔG*OOH − ΔG*O	
*OOH → * + O2 + H+ + e−	ΔG4 = 4.92 − ΔG*OOH	
ORR	* + O2 + H+ + e− → *OOH	ΔG1 = ΔG*OOH − 4.92		
*OOH + H+ + e− → H2O(1) + *O	ΔG2 = ΔG*O − ΔG*OOH	
*O + H+ + e− → *OH	ΔG3 = ΔG*OH − ΔG*O	
*OH + H+ + e− → H2O(1) + *	ΔG4 = –ΔG*OH	
a Where ΔG*OH, ΔG*O, and ΔG*OOH are the adsorption energies of *OH, *O, and *OOH intermediates. These adsorption energies are defined as follows: ΔG*OH = G*OH + 0.5GH2 − G* − GH2O, ΔG*O = G*O + GH2 − G* − GH2O, ΔG*OOH = G*OOH + 1.5GH2 − G* − 2GH2O, where GH2O and GH2 are the free energies of H2O and H2 molecules, respectively.

Furthermore, the OER/ORR catalytic activity can be directly determined by the adsorption Gibbs free energies of *OH, *O and *OOH intermediates (ΔG*OH, ΔG*O and ΔG*OOH), and too-strong or too-weak bonding between intermediates and reaction sites negatively impact the proceeding of reactions. Accordingly, ΔG*OH, ΔG*O and ΔG*OOH of pristine and TM-doped ZnS and ZnSe nanostructures for OER/ORR were calculated using the computational hydrogen electrode model proposed by Nørskov43 and presented in Table 2, which are indispensable for evaluating the ORR/OER activity. Moreover, the zero point energy (ZPE) corrections and entropy contributions (TS) of *OH, *O, and *OOH on the surface of ZnX (X = S and Se) nanostructures at 298.15 K by DFT calculations are presented in Tables S3 and S4.† As expected, *OH, *O, and *OOH intermediates are preferably adsorbed on the TM site of ZnS and ZnSe nanostructures owing to the moderate adsorption of Gibbs free energy (shown in Table 2), and the corresponding configurations are displayed in Fig. S3 and S4.† Moreover, the electrochemical interface between ZnX (X = S and Se) and water was investigated to determine how the inclusion of an explicit solvent (undissociated water molecule) influences the stability of adsorbed intermediates (*OH, *O and *OOH) in the OER and ORR. The results indicated that the inclusion of an explicit solvent can considerably stabilize the structure of adsorbed intermediates (*OH, *O and *OOH), which are capable of donating hydrogen bonds (as shown in Fig. S5†).

The adsorption-free energies of *OH, *O, and *OOH intermediates on the surface of ZnS and ZnSe nanostructures

Systems	ZnS	ZnSe	
ΔG*OH/eV	ΔG*O/eV	ΔG*OOH/eV	ΔG*OH/eV	ΔG*O/eV	ΔG*OOH/eV	
Pristine	1.89	4.18	4.73	1.71	2.27	4.67	
Fe-doped	0.01	1.43	3.29	0.02	1.43	3.36	
Co-doped	0.50	1.76	3.69	0.46	1.61	3.67	
Ni-doped	1.02	2.49	3.99	0.92	2.43	3.98	
Cu-doped	1.72	3.74	4.73	1.82	2.47	4.75	
Pd-doped	1.80	2.94	4.16	1.14	2.91	4.25	
Pt-doped	1.24	2.56	4.06	0.85	2.39	4.12	

To clarify OER/ORR catalytic activity and pathways, free-energy diagrams of pristine and TM-doped ZnS and ZnSe nanostructures for OER/ORR at U = 0 and 1.23 V are plotted and listed in Fig. 3–5. As mentioned above, the OER process involves four steps, including the formation of H2O → *OH → *O → *OOH → O2, while the ORR process is reversed. For pristine ZnS and ZnSe materials during the OER process, the potential-determining step (PDS) is the formation of *OH to *O and *O to *OOH, corresponding to the overpotential (ηOER) of 1.05 and 1.17 V, respectively. The steps of *OH to *O and *O to *OOH are uphill and endothermic at an equilibrium potential of 1.23 V, suggesting that additional applied electrode potential is required to ensure all the elementary steps are spontaneous. However, the ORR PDS is O2 → *OOH both for pristine ZnS and ZnSe nanostructures, and this step is downhill and endothermic, corresponding to the overpotentials (ηORR) of 1.04 and 0.98 V, respectively. The high ηOER (1.05 and 1.17 V) and ηORR (1.04 and 0.98 V) values suggest that the pristine ZnS and ZnSe nanostructures are not suitable for efficient oxygen electrocatalysts. Moreover, the associative and dissociative mechanism is involved in the ORR and thus the interaction between O2 molecule and ZnX (X = S and Se) nanostructures was considered. Our calculation suggests that, for both ZnS and ZnSe nanostructures, the *O2 dissociative pathway (*O2 → *O–*O) is endothermic, while the associative pathway (*O2 → *OOH) is exothermic, as shown in Fig. S6 in the ESI.† Therefore, the associative mechanism will be concerned with Fe, Co, Ni, Cu, Pd and Pt-doped ZnX (X = S and Se) nanostructures during the ORR process.

Fig. 3 Free energy diagram of (a) ZnS and (b) ZnSe nanostructures for OER and ORR at U = 0 and 1.23 V. The red and blue lines are the potential-determining step (PDS) for OER and ORR, respectively.

Fig. 4 Free energy diagram of (a) Fe, (b) Co, (c) Ni, (d) Cu, (e) Pd and (f) Pt-doped ZnS nanostructures for OER and ORR at U = 0 and 1.23 V, respectively. The red and blue lines are the potential-determining steps (PDS) for OER and ORR.

Fig. 5 Free energy diagram of (a) Fe, (b) Co, (c) Ni, (d) Cu, (e) Pd and (f) Pt-doped ZnSe nanostructures for OER and ORR at U = 0 and 1.23 V. The red and blue lines are the potential-determining steps (PDS) for OER and ORR, respectively.

As shown in Fig. 4 and 5, the four steps of the OER on the surface of Fe, Co, Ni, Cu, Pd and Pt-doped ZnS and ZnSe nanostructures were uphill at U of 0 V, confirming that the whole process was endothermic. When the electrode potential is applied to the equilibrium potential of 1.23 V, there are energy barriers for OER (Fig. 4 and 5), disclosing that additional applied electrode potential is required to make all the elementary steps spontaneous. Interestingly, the formation of *O to *OOH is the PDS for all TM-doped ZnS and ZnSe nanostructures except for Cu/ZnS, Pd/ZnS and Pd/ZnSe during the OER process. For the ORR process, we found that the PDS of Fe and Co-doped ZnS occurred in the process of *OH → H2O, while that of Ni, Cu, Pd and Pt-doped ZnS is the formation of O2 to *OOH (Fig. 4). As can be seen in Fig. 5, during the process of ORR, the PDS of Fe-, Co- and Ni-doped ZnSe is the formation of *OH → H2O, while that of Cu- Pd- and Pt-doped ZnSe is the formation of O2 → *OOH. In general, the adsorption strength of the reaction intermediates plays a crucial role in the OER/ORR electrocatalytic activity. According to the adsorption Gibbs free energies of *OH, *O and *OOH calculated using the equations in Table 1, the computed ηOER and ηORR on Fe, Co, Ni, Cu, Pd and Pt-doped ZnS and ZnSe nanostructures are labeled in Fig. 4 and 5, respectively. Note that the OER/ORR overpotentials of Ni-doped ZnS and ZnSe are 0.30/0.28 and 0.31/0.31 V, respectively, which are lower than those of IrO2 (ηOER = 0.56 V) and Pt (ηORR = 0.45 V).51 This phenomenon shows that Ni/ZnS and Ni/ZnSe can be practically used as promising bifunctional oxygen electrocatalysts. Moreover, Pd and Pt substitutional doping obviously decreases the overpotential of ZnS and ZnSe nanostructures for OER and ORR, indicating their potential bifunctional catalytic performance.

3.3 Origin of OER/ORR catalytic activity

It has been well accepted that the scaling relationships among the adsorption-free energies of *OH, *O, and *OOH intermediates can provide insight into understanding the origin of catalytic activity.52 According to the above discussions, all the *OH, *O, and *OOH intermediates preferentially adsorb on the TM atoms site with TM–O binding (Fig. 1a and S1a†). Therefore, the scaling relationships of ΔG*OH with ΔG*O and ΔG*OOH were plotted and are displayed in Fig. 6 and 7. Remarkably, the significant correlation linear relationships among ΔG*OH, ΔG*O and ΔG*OOH can be obtained on ZnS nanostructures, namely, ΔG*O = 1.29ΔG*OH + 1.22 eV (R2 = 0.83) and ΔG*OOH = 0.68ΔG*OH + 3.30 eV (R2 = 0.83), as shown in Fig. 6a and b. However, the linear relationship between ΔG*OH and ΔG*O is poor with a weak determination coefficient (R2 = 0.39) for ZnSe nanostructures, while that between ΔG*OH and ΔG*OOH is good with a strong determination coefficient (R2 = 0.98). According to the universally recognized Sabatier principle, the adsorption strength of the *OH, *O, and *OOH intermediates plays a crucial role in the catalyst's OER/ORR electrochemical activity.53–56 Too-strong or too-weak adsorption of intermediates (*OH, *O and *OOH) on the surface of the catalyst will negatively impact the proceeding of OER and ORR. This indicates that an ideal catalyst should have a moderate ΔGi value to ensure suitable adsorption of intermediates during the reaction process. As shown in Table 2, the calculated adsorption Gibbs free energies of *OH, *O, and *OOH intermediates for Ni/ZnS (ΔG*OH = 1.02 eV, ΔG*O = 2.49 eV, ΔG*OOH = 3.99 eV) and Ni/ZnSe (ΔG*OH = 0.92 eV, ΔG*O = 2.43 eV, ΔG*OOH = 3.98 eV) are approximately equivalent to the optimal values for ideal OER/ORR catalysts.

Fig. 6 Scaling relationship between the adsorption-free energies of (a) ΔG*OHvs. ΔG*O and (b) ΔG*OHvs. ΔG*OOH on the surface of pristine and TM-doped ZnS nanostructures. The volcano plot of (c) ηOERvs. ΔG*O − ΔG*OH and (d) ηORRvs. ΔG*OH for OER and ORR on the surface of pristine and TM-doped ZnS nanostructures.

Fig. 7 Scaling relationship between the adsorption-free energies of (a) ΔG*OHvs. ΔG*O and (b) ΔG*OHvs. ΔG*OOH on the surface of pristine and TM-doped ZnSe nanostructures. The volcano plot of (c) ηOERvs. ΔG*O − ΔG*OH and (d) ηORRvs. ΔG*OH for OER and ORR on the surface of pristine and TM-doped ZnSe nanostructures.

Moreover, the volcano plots are considered a powerful guide to select and develop efficient electrocatalysts. Thus, the overpotentials volcano plots are displayed in Fig. 6c, d, 7c and d, which are established by ΔG*OH and ΔG*O − ΔG*OH as the descriptors of the OER/ORR catalytic activity of different ZnS and ZnSe nanostructures. Usually, the best-performance catalysts with the most moderate adsorption energy are located at the top of the volcano plot. The best OER catalyst has the ΔG*O − ΔG*OH value of ∼1.50 eV, and the ΔG*OH value of ∼1.00 eV is required for the best ORR catalyst. Obviously, the Ni-doped ZnS and ZnSe are at the top of the OER and ORR volcano plots, elucidating their superior catalytic performance for OER/ORR and implying the promising potential as bifunctional oxygen electrocatalysts. Moreover, Pt/ZnS (ηOER/ηORR = 0.27/0.37 V) and Pt/ZnSe (ηOER/ηORR = 0.43/0.50 V) are close to the peak of volcano plots with lower overpotentials, which is considerably better than other ZnS and ZnSe nanostructures for OER and ORR. As discussed above, the lower overpotential is favorable for OER and ORR, which can be used as an indicator of the catalytic performance of the electrocatalysts. Remarkably, Ni-doped ZnS and ZnSe nanostructures possess the lowest OER and ORR overpotential, suggesting that they are both selective for OER and ORR. Our results demonstrate that Ni/ZnS and Ni/ZnSe nanostructures can act as excellent bifunctional catalysts for OER and ORR, which can be beneficial for synthesizing efficient performance electrocatalysts in the future.

The insight into the origin of OER/ORR catalytic activity can provide a theoretical guide to improve the performance of catalysts. In general, the adsorption of *OH, *O, and *OOH intermediates on catalytic active sites is the key factor influencing the catalytic activity of OER/ORR. However, the d-band center of transition metal atoms (εd, the average position of the TM-d orbital) has been widely used to describe the bond strength between the adsorbed intermediates (*OH, *O and *OOH) and the surface of catalysts.57 Therefore, to understand the influence of different TM atoms as catalytic active sites on the catalyst performance, the εd of Fe, Co, Ni, Cu, Pd, and Pt in our studied systems were calculated and are shown in Fig. S7.† It can be found that different TM (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe can result in the different adsorption strengthen of intermediates. In particular, the Ni atom in the Ni/ZnS and Ni/ZnSe nanostructures with the moderate d-band center position displays appropriate adsorption strength, which leads to the lower overpotential (ηOER/ηORR = 0.28/0.30 V), suggesting the promising catalytic performance for OER/ORR.

4. Conclusions

In conclusion, the catalytic performance of TM (Fe, Co, Ni, Cu, Pd and Pt)-doped ZnS and ZnSe nanostructures for OER/ORR has been systematically evaluated by high-throughput theoretical screening. Our calculated results suggest that Fe, Co, Ni, Cu, Pd and Pt can be treated as the catalytic active sites on the surface of ZnS and ZnSe nanostructures during the OER/ORR process, which is available for the formation of *OH, *O and *OOH intermediates. Moreover, the free-energy diagrams show that the pristine ZnS and ZnSe nanostructures own high OER/ORR overpotentials, implying inefficient catalytic performance. In contrast, the OER/ORR overpotentials of Co, Ni, Pd and Pt-doped ZnS and ZnSe nanostructures obviously decrease. Particularly, the Ni/ZnS and Ni/ZnSe nanostructures own the lower OER/ORR overpotentials of 0.28/0.30 V and 0.31/0.31 V, respectively, demonstrating their promising bifunctional catalytic activity. Our findings not only provide a theoretical strategy to design the bifunctional oxygen electrocatalysts but also shed insight into the electrocatalytic application of ZnS and ZnSe-based materials in water electrolysis.

Data availability

The authors confirm that the data supporting the findings of this study are available within its ESI.†

Conflicts of interest

There are no conflicts to declare.

Supplementary Material

RA-014-D4RA04011B-s001

This work was supported by the National Natural Science Foundation of China (No. 21703087).
==== Refs
Notes and references

Chu S. Cui Y. Liu N. Nat. Mater. 2017 16 16 22
Panwar N. L. Kaushik S. C. Kothari S. Renewable Sustainable Energy Rev. 2011 15 1513 1524
Liu G. S. You S. J. Tan Y. Ren N. Q. Environ. Sci. Technol. 2017 51 2339 2346 28117980
Shu Z. Chen H. F. Liu X. Jia H. X. Yan H. J. Cai Y. Q. Adv. Funct. Mater. 2023 35 2301493
Shu Z. Yan H. J. Chen H. F. Cai Y. Q. J. Mater. Chem. A 2022 10 5470 5478
Rai V. Lee K. P. Safanama D. Adams S. Blackwood D. J. ACS Appl. Energy Mater. 2020 3 9417 9427
Tahir M. Pan L. Idrees F. Zhang X. Wang L. Zou J.-J. Wang Z. L. Nano Energy 2017 37 136 157
Stacy J. Regmi Y. N. Leonard B. Fan M. Renewable Sustainable Energy Rev. 2017 69 401 414
Sun H. N. Jung W. C. J. Mater. Chem. A 2021 9 15506 15521
Ying J. Chen J. B. Xiao Y. X. Cordoba de Torresi S. I. Ozoemena K. I. Yang X. Y. J. Mater. Chem. A 2023 11 1634 1650
Dang Q. Lin H. P. Fan Z. L. Ma L. Shao Q. Ji Y. J. Zheng F. F. Geng S. Z. Yang S. Z. Kong N. N. Zhu W. X. Li Y. Y. Liao F. Huang X. Q. Shao M. W. Nat. Commun. 2021 12 6007 34650084
Zhang C. L. Shen X. C. Pan Y. B. Peng Z. M. Front. Energy 2017 11 268 285
Zhang J. W. Yuan Y. L. Gao L. Zeng G. M. Li M. F. Huang H. W. Adv. Mater. 2021 33 2006494
Majhi K. C. Yadav M. ACS Eng. Au 2023 3 278 284
Ingsel T. and Gupta R. K. , Nanomaterials for Electrocatalysis (Micro and Nano Technologies) . 2022, pp. 83–111
Liu Y. J. Guo Y. J. Liu Y. R. Wei Z. L. Wang K. Shi Z. Q. Energy Fuels 2023 37 2608 2630
Azam A. Yang J. Li W. X. Huang J.-K. Li S. Prog. Mater. Sci. 2022 132 101042
Zhang P. Xiang H. Y. Tao L. Dong H. J. Zhou Y. G. Hu T. S. Chen X. L. Liu S. Wang S. Y. Garaj S. Nano Energy 2019 57 535 541
Zheng X. Han X. Liu H. Chen J. Fu D. Wang J. Zhong C. Deng Y. Hu W. ACS Appl. Mater. Interfaces 2018 10 13675 13684 29616794
Jo S. Lee K. B. Sohn J. I. ACS Sustainable Chem. Eng. 2021 9 14911 14917
Wang Z. Y. Liu S. L. Duan W. Xing Y. C. L Hu Y. Ma Y. J. Int. J. Hydrogen Energy 2024 60 1414 1432
Peng X. Yan Y. J. Jin X. Huang C. Jin W. H. Gao B. Chu P. K. Nano Energy 2020 78 105234
Chen M. Yang Y. J. Xiong B. Huang Y. H. Liu J. Wang G. B. Int. J. Hydrogen Energy 2024 51B 540 550
Cui L. F. Xiang K. Kang X. M. Zhi K. K. Wang L. Zhang J. J. Fu X. Z. Luo J. L. J. Colloid Interface Sci. 2022 609 868 877 34839920
Yang L. F. Yao X. Y. Du C. L. Han Z. L. Jin M. W. Peng S. C. Ma X. L. Zhu Y. Q. Zou M. S. Cao C. B. Chem. Eng. J. 2024 481 148598
Xia F. F. Shu L. Wen Y. P. Yang F. L. Zheng C. Z. Mater. Adv. 2022 3 5772 5777
Hu L. Zhong P. Zhang X. Xiang Y. Balogun M.-S. Tong Y. X. Yang H. Appl. Surf. Sci. 2023 623 157040
Hussain S. Yang X. T. Yang J. J. Li Q. Y. Mater. Today Sustain. 2024 25 100686
Ede S. R. Luo Z. P. J. Mater. Chem. A 2021 9 20131 20163
Du X. Q. Su H. Zhang X. H. ACS Sustainable Chem. Eng. 2019 7 16917 16926
Hohenberg P. Kohn W. Phys. Rev. B 1964 136 B864 B871
Kohn W. Sham L. Phys. Rev. 1965 140 A1133 A1138
Kresse G. Joubert D. Phys. Rev. B: Condens. Matter Mater. Phys. 1999 59 1758
Kresse G. Furthmüller J. Phys. Rev. B: Condens. Matter Mater. Phys. 1996 54 11169
Perdew J. P. Burke K. Ernzerhof M. Phys. Rev. Lett. 1996 77 3865 3868 10062328
Perdew J. Wang Y. Phys. Rev. B: Condens. Matter Mater. Phys. 1992 45 13244 13249
Hammer B. Jacobsen K. W. Nørskov J. K. Phys. Rev. Lett. 1993 70 3971 3974 10054012
Cococcioni M. de Gironcoli S. Phys. Rev. B: Condens. Matter Mater. Phys. 2005 71 035105
Qin Z. Wang Z. X. Zhao J. X. Nanoscale 2022 14 6902 6911 35446333
Wang Y. W. Tian W. Wan J. Xiong G. Q. Zhang H. J. Wang Y. Phys. Chem. Chem. Phys. 2022 24 10325 10333 35438086
Grimme S. J. Comput. Chem. 2006 27 1787 1799 16955487
Monkhorst H. J. Pack J. D. Phys. Rev. B: Solid State 1976 13 5188 5192
Nørskov J. K. Bligaard T. Logadottir A. Kitchin J. R. Chen J. G. Pandelov S. Stimming U. J. Electrochem. Soc. 2005 152 J23
Nørskov J. K. Rossmeisl J. Logadottir A. Lindqvist L. Kitchin J. R. Bligaard T. Jonsson H. J. Phys. Chem. B 2004 108 17886 17892
Wang V. Xu N. Liu J. C. Tang G. Geng W. T. Comput. Phys. Commun. 2021 267 108033
Wang Y. Mao J. Meng X. G. Yu L. Deng D. Bao X. H. Chem. Rev. 2019 119 1806 1854 30575386
Wang D. Y. Wang K. X. Liu Q. T. Liu M. Y. Zhang G. Q. Feng K. Wang K. Ding X. W. Zhu H. M. Yang S. Liu Y. H. Li T. H. Gong P. Wang M. L. Wang P. G. Jin H. Z. Zhao W. Yu F. Adv. Sci. 2020 7 1901382
Feng C. Faheem M. B. Fu J. Xiao Y. Q. Li C. L. Li Y. B. ACS Catal. 2020 10 4019 4047
Spöri C. Briois P. Nong H. N. Reier T. Billard A. Kühl S. Teschner D. Strasser P. ACS Catal. 2019 9 6653 6663
Garg P. Nair A. S. Rawat K. S. Pathak B. J. Phys. Chem. C 2019 123 13419 13427
Rossmeisl J. Qu Z. W. Zhu H. Kroes G. J. Nørskov J. K. J. Electroanal. Chem. 2007 607 83 89
Qin Z. Zhao J. J. Colloid Interface Sci. 2022 605 155 162 34311310
Kari J. Olsen J. P. Jensen K. Badino S. F. Krogh K. B. R. M. Borch K. Westh P. ACS Catal. 2018 8 11966 11972
Bajdich M. García-Mota M. Vojvodic A. Nørskov J. K. Bell A. T. J. Am. Chem. Soc. 2013 135 13521 13530 23944254
Man I. C. Su H. Y. Calle-Vallejo F. Hansen H. A. Martínez J. I. Inoglu N. G. Kitchin J. Jaramillo T. F. Nørskov J. K. Rossmeisl J. ChemCatChem 2011 3 1159 1165
Zhang T. Zhang B. Peng Q. Zhou J. Sun Z. J. Mater. Chem. A 2021 9 433 441
Niu H. Wan X. H. Wang X. T. Shao C. Robertson J. Zhang Z. F. Guo Y. Z. ACS Sustain. Chem. Eng. 2021 9 3590 3599
