Next Article in Journal
Impact of Wall Material-to-Active Ratio in the Stability of Spray-Dried Ascorbic Acid Using Maltodextrin and Gum Arabic
Previous Article in Journal
Enhancing Functional Properties and Protein Structure of Almond Protein Isolate Using High-Power Ultrasound Treatment
Previous Article in Special Issue
Computational Design of Ni6@Pt1M31 Clusters for Multifunctional Electrocatalysts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molybdenum and Vanadium-Codoped Cobalt Carbonate Nanosheets Deposited on Nickel Foam as a High-Efficient Bifunctional Catalyst for Overall Alkaline Water Splitting

1
State Key Laboratory of Analytical Chemistry for Life Science, Nanjing University, Nanjing 210023, China
2
School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
3
Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing 210023, China
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(15), 3591; https://doi.org/10.3390/molecules29153591
Submission received: 8 July 2024 / Revised: 25 July 2024 / Accepted: 28 July 2024 / Published: 30 July 2024

Abstract

:
To address issues of global energy sustainability, it is essential to develop highly efficient bifunctional transition metal-based electrocatalysts to accelerate the kinetics of both the hydrogen evolution reaction (HER) and the oxygen evolution reaction (OER). Herein, the heterogeneous molybdenum and vanadium codoped cobalt carbonate nanosheets loaded on nickel foam (VMoCoCOx@NF) are fabricated by facile hydrothermal deposition. Firstly, the mole ratio of V/Mo/Co in the composite is optimized by response surface methodology (RSM). When the optimized composite serves as a bifunctional catalyst, the water-splitting current density achieves 10 mA cm−2 and 100 mA cm−2 at cell voltages of 1.54 V and 1.61 V in a 1.0 M KOH electrolyte with robust stability. Furthermore, characterization is carried out using field emission scanning electron microscopy-energy dispersive spectroscopy (FESEM-EDS), high-resolution transmission electron microscopy (HRTEM), X-ray diffraction (XRD), and X-ray photoelectron spectroscopy (XPS). Density functional theory (DFT) calculations reveal that the fabricated VMoCoCOx@NF catalyst synergistically decreases the Gibbs free energy of hydrogen and oxygen-containing intermediates, thus accelerating OER/HER catalytic kinetics. Benefiting from the concerted advantages of porous NF substrates and clustered VMoCoCOx nanosheets, the fabricated catalyst exhibits superior electrocatalytic performance. This work presents a novel approach to developing transition metal catalysts for overall water splitting.

Graphical Abstract

1. Introduction

Recently, increasing populations and rapid consumption of fuel for industrial and agricultural production have boosted the energy crisis [1]. Electrochemical water electrolysis is an important strategy to convert water to hydrogen (H2) and oxygen (O2) in today’s need for clean and recyclable energy that has a high capability for storage [1]. H2 is one of the most important clean energy fuels, which can be used in fuel cells to produce electricity or as a fuel gas in heating devices. For large-scale H2 production, electrochemical water splitting offers the dual benefits of sustainable clean energy generation and an effective reduction in carbon discharge. Highly effective electrocatalysts are required in water electrolysis to accelerate the kinetics of hydrogen and oxygen evolution (HER and OER). Currently, noble metal-based materials are promising catalysts for use in electrochemical water splitting [2]. Platinum (Pt) is the best-ever high-efficiency HER catalyst, while IrO2 and RuO2 are deemed the best-known OER catalysts. However, due to their high cost and scarcity, these materials lack the feasibility of developing scalable electrocatalysts for new energy production [3,4]. To this end, it is urgent to develop durable, cheap, and efficient electrocatalysts from earth-abundant transition metals, which is vital for the development of the H2 economy [5].
Transition metal-based electrocatalysts, due to their Earth-abundant resources, cheapness, and comparable catalytic performance, have attracted enormous interest as alternatives to expensive and scarce noble metals [5]. However, their HER performances are not satisfactory under alkaline media, and the excessively high Gibbs energy barriers (ΔG) limit their HER activity. The OER suffers from more sluggish charge/electron transfer kinetics because the OER is a four-step charge/electron transfer, resulting in a high requirement for electrocatalysts. Therefore, developing cheap, stable, and active bifunctional electrocatalysts for the achievement of highly efficient electrochemical H2 production is challenging. For this purpose, transition metals/transition metallic compounds are introduced as potential materials for developing highly efficient electrocatalysts for water electrolysis. Transition metallic bifunctional electrocatalysts are advantageous owing to their special electronic structures of atoms. The electronic states of bimetallic/multimetallic catalysts can be tuned by heteroatom codoping [6,7] or oxygen vacancies [8]. Generally, heteroatom codoping is an effective approach for bifunctional non-noble metal catalysts to accelerate charge/electron transfer and expose more active sites.
Compared with expensive and scarce noble metals or noble metallic compounds, the lower band gaps and stronger electron-donating abilities of transition metals such as vanadium (V), molybdenum (Mo), and cobalt (Co) make them excellent dopant candidates to enhance catalytic activity [5]. Their electrocatalytic activities are mainly dependent on their electronic structures, active sites, charge/electron transfer, and surface energy/surface area [9]. The co-element has attracted great attention for its ability to accelerate charge/electron diffusion, resulting in high efficiency in overall water splitting [10]. Mo-based electrocatalytic materials have been fabricated as efficient HER and OER catalysts in alkaline media due to their multiple valence states and unique redox properties. V element with multiple valence states, rich earth abundance, and high electrochemical performance is a promising OER dopant to effectively tune the synergistic effect among different metal active sites [11]. Zhou et al. fabricated V-doped Ni3S2@NiFe LDH ultrathin nanosheets on porous Ni foam (NF) substrates as a highly effective electrocatalyst for water splitting. The presence of V guarantees a remarkable promotion of electrocatalytic kinetics and activity [12]. Yet, the intrinsic performance of HER was still inferior due to the weak Gibbs free energy change in hydrogen adsorption (ΔGH*) on the V atom [8]. The addition of transition metals can modify the electronic structure of electrocatalysts, thus exposing more active sites. Xiao et al. fabricated the MoCo(OH)2/CoP@NF, which exhibited enhanced catalytic activity toward the OER [13]. Yan et al. combined Co with Mo2C, which could decrease the ΔGH*, thus improving the catalytic activity [14]. Inspired by the above outcomes, multi-transition-metal composites can synergistically optimize the HER/OER catalytic activity. Doping high-valence V and Mo ions can optimize the electronic structure of transition metal materials owing to their electron-drawing capacity. The nanosheet structure of cobalt carbonate facilitates ion/electron diffusion and exposes more active sites. Moreover, porous NF is an intriguing substrate for enhancing electronic states and active sites due to its porous surface area [15]. To the best of our knowledge, the development of multimetallic codoped bifunctional electrocatalysts loaded on porous NF substrates for bifunctional water splitting has rarely been reported till now.
In this work, the V and Mo-codoped CoCOx nanosheets in situ loaded on a porous NF substrate (VMoCoCOx@NF) were fabricated via facile hydrothermal deposition (Scheme 1). Firstly, response surface methodology (RSM) was employed to achieve the optimum mole ratio of V/Mo/Co in the catalyst. More detailed illustrations of RSM, as well as related data graphs (Figures S1 and S2), were given in the Supporting Information (SI). The predicted RSM model responses of VMoCoCOx@NF were well-fitted for all the actual responses (Figure S2). The optimized catalyst showed remarkable HER/OER activities in a 1.0 mol L−1 (M) KOH electrolyte. Meanwhile, density functional theory (DFT) was used to fully illustrate the high HER/OER activities of the fabricated catalyst. Additionally, the excellent HER/OER activities for VMoCoCOx@NF enabled it to serve as both anode and cathode for the alkaline electrolyzer, a catalyst that required cell voltages of 1.54 V/1.61 V to reach current densities of 10 mA cm−2/100 mA cm−2 with excellent cycling stability.

2. Results and Discussion

2.1. Characterizations of Catalysts

The structure of the VMoCoCOx@NF was confirmed by X-ray diffraction (XRD) and high-resolution transmission electron microscopy (HRTEM). The crystallinity of the catalysts was characterized by the XRD technique within the 2θ range of 10–80°. The XRD patterns of VMoCoCOx@NF and VCoCOx@NF and the XRD positions for pure MoO3, CoCO3, and V2O5 references were depicted in Figure 1A. The strong diffraction lines at 44.2°, 51.8°, and 76.3° correspond to the standard pattern of Ni (PDF#04-0850) in NF substrate. For comparison, the XRD positions for pure MoO3, CoCO3, and V2O5 references were indicated according to the PDF# database. In Table S7, the diffraction lines at 19.9°, 26.5°, 30.8°, 33.7°, 47.1°, and 50.8° correspond to V2O5 (PDF#00-041-1426), the diffraction lines at 23.3°, 25.7°, and 27.3° are well indexed to MoO3 (PDF#05-0508), and the diffraction lines at 42.9°, 46.5°, 51.8°, and 53.5° are attributed to CoCO3 (PDF#11-0692). The correspondence with MoO3 marked on the picture was observed at 23.3°, which corresponds to MoO3 (PDF#05-0508). Importantly, in Figure 1A and Figure S3, the diffraction lines of Ni were removed during the XRD operation. The diffraction lines of VMoCoCOx@NF were observed at approximately 23.3°, 26.5°, 32.5°, 46.5°, and 53.5°. The diffraction lines of VCoCOx@NF were observed at approximately 26.5°, 33.7°, and 46.5°, indicating that the characteristic lines of MoO3, V2O5, and CoCO3 were visible at 23.3°, 26.5°, and 46.5° in the synthesized hybrid catalyst VMoCoCOx@NF. The difference was that the additional line towards the diffraction angle of 32.5° was for VMoCoCOx@NF, indicating the formation of a new phase (112) [16]. It suggests that, compared with VCoCOx@NF, an appropriate amount of Mo was assembled into the hybrid. The corresponding diffraction line at 32.5° might arise from the doping of Mo.
To further analyze the morphology of the sample, the TEM images in Figure 1B show crystalline nanosheets embedded in the NF substrate. The portions with agglomerates should be the nano-edge structure of NF. The HRTEM image of VMoCoCOx@NF reveals the well-resolved lattice fringes of 0.190 nm and 0.261 nm, which are in agreement with the (310) and (112) planes of VMoCoCOx@NF (Figure 1C and Figure S4A) [17], which is also in agreement with the XRD observation in Figure 1A. Notably, Figure S4B,C also clearly exhibits that the accurate measured values of 10 neighbor lattice fringes in the corresponding crystal plane are 1.90 nm and 2.61 nm, which are exactly ten times the single d-spacing values measured by HRTEM. The crystal structure of the VMoCoCOx@NF catalyst was constructed via strong coordination bond interaction among V, Mo, and Co; the lamellas were combined by an electrostatic force. Taking advantage of the synergistic effects of 3d and 4d vacancies and V/Mo codoping, the doped nanostructure was effectively modulated, and more active sites were exposed.
Moreover, the field emission scanning electron microscopy-energy dispersive spectroscopy (FESEM-EDS) elemental mapping images in Figure 1D–F validate the distribution of V, Mo, and Co elements over the surface of VMoCoCOx@NF. The EDS mappings indicate that the structure was mainly determined by CoCO3 since V2O5 and MoO3 were comparably more randomly distributed on the NF. Figure 1G–I characterize the morphology and nanostructure of the fabricated VMoCoCOx@NF catalyst by FESEM. The clustered hexagonal VMoCoCOx nanosheets were anchored on the surface of NF. It presented an interconnected porous architecture with abundant charge/electron transfer channels. From the FESEM image of VMoCoCOx@NF (Figure 1G), it can be seen that a small number of nanoparticles, that is, CoCO3 (Ksp = 1.4 × 10−13), are deposited onto the VMoCoCOx nanosheet surface, further being confirmed by the Fourier transform infrared spectrum (FTIR) of VMoCoCOx@NF (Figure S5). FTIR spectroscopy was utilized to further identify the functional groups on the surface of VMoCoCOx@NF, VCoCOx@NF, and V2O5@NF. The FTIR spectrum of VMoCoCOx@NF shows that the peaks at 1667.00 cm−1, 1078.15 cm−1, and 833 cm−1 are assigned to the stretching vibration of the C=O bond, the C-O bond, and the bending vibration of the CO32− functional group. The peak around 1100 cm−1 is distinctive for VMoCoCOx@NF. The peak at 478.7 cm−1 is attributed to the vibration of the VO43− functional group [18]. The FTIR measurement confirms the CO32− existing in the hybrid catalyst VMoCoCOx@NF. The FESEM image in Figure 1H demonstrates clustered hexagonal ultrathin nanosheet morphology with a thickness of approximately 100 nm. Notably, such a unique structure provides irregular charge/electron transfer channels and an abundant porous architecture, which is not only beneficial to the diffusion/penetration of electrolytes and the escape of H2 and O2 but also offers a large surface area with abundant exposed active sites for reaction intermediates.
The full survey X-ray photoelectron spectroscopy (XPS) spectra of V 2p, Co 2p, and Mo 3d for the fabricated VMoCoCOx@NF catalyst are shown in Figure 2A, confirming the main peaks and valence states of V, Mo, and Co elements. The Ni peak was removed during the XPS operation. In the analysis of the XRD results shown in Figure 1A, all the nickel peak intensities were very large. Considering that a large amount of Ni from the NF substrate could be easily introduced into the hydrothermal system, the Ni peak was omitted during the XPS scanning process due to the visualization effect of the target peaks. Figure 2B shows the V2p spectrum; the two peaks were fitted into V 2p3/2 and V 2p1/2 located at 516.9 eV and 524.4 eV with a spin–orbit splitting of 7.5 eV, corresponding to the occurrence of V+5 (VO43−) species. The V on the surface mainly existed in the V5+ species. The high-valence V5+ species is crucial to higher electric conductivity, which benefits the improvement of catalytic performance, which is consistent with previous works [19,20]. In the XPS diagram of the 2p orbital of Co (Figure 2C), 2p1/2 and 2p3/2 only deconvolve one peak. The major peaks of 781.7 and 782.1 eV are consistent with the deconvolution of Co 2p 3/2. Simultaneously, peaks of 796.4 and 797.09 eV conform to the spin–orbit characteristics of Co2p 1/2, and the two vibrating satellite peaks are observed at 786.23 eV and 802.61 eV for high-resolution Co 2p spectra [21]. Interestingly, the spin–orbit splitting of about 16 eV between Co 2p 3/2 and Co 2p 1/2 peaks is observed to prove the Co2+ species in VMoCoCOx@NF [10,17]. The Mo doping caused an upshift of Co 2p binding energy in VMoCoCOx@NF, indicating that the electronic structure of Co was greatly adjusted by the charge transfer from the electron cloud around Co to the doped atoms [4]. It is essential to verify the efficient doping of the elements and evaluate a possible synergistic effect that might influence HER/OER activity [22]. Figure 2D shows the Mo 3d region XPS spectrum for VMoCoCOx@NF, in which only two peaks of Mo 3d5/2 and Mo 3d3/2 at 232.5 eV and 235.6 eV were identified and assigned to Mo6+ (MoO42−). However, the peaks at ~228 eV and 232 eV, which correspond to lower bonding energies of Mo4+ (3d5/2 and 3d3/2), were not observed, indicating that Mo4+ did not exist in the catalyst [23]. In Figure 1A, VMoCoCOx@NF has a small diffraction peak at 23.3°. The XRD analysis shows the presence of MoO3 (PDF#05-0508). This result is also consistent with the result of the XPS analysis (Figure 2D), confirming the presence of only Mo6+ (MoO42−) species. The emergence of the MoO42− peaks suggested that Mo was doped into the catalyst. Doping high-valence Mo6+ and V5+ ions into CoCO3 altered the electronic structure of the 3d metal atom, owing to the capability of Mo6+ and V5+ to draw electrons. The XPS and FTIR results confirmed the MoO42−, VO43−, Co2+, and CO32− species presented in VMoCoCOx@NF. The synergistic effect should help elucidate the atomic interactions among V, Mo, and Co.

2.2. Electrocatalytic HER Performance

The electrocatalytic HER performance of the as-prepared electrocatalysts was measured using the standard three-electrode system in the N2-saturated 1.0 M KOH electrolyte. For comparison, the HER activities of Pt/C@NF, VCoCOx@NF, MoCO3@NF, CoCO3@NF, V2O5@NF, and NF were studied. Figure 3A displays the HER linear sweep voltammetry (LSV) curves of all samples recorded with a scan rate of 5 mV s−1 from 0.00 mV to −100 mV. As expected, the VMoCoCOx@NF catalyst displayed excellent alkaline HER activity, which agreed with the previous reports [17,22] (Table S8). The overpotential and Tafel slope are the most important parameters for clarifying the reaction kinetics of catalysts. Doping of two or more elements (salts) could efficiently facilitate the abundance of active sites, which resulted in high HER activity and significantly reduced the overpotential and Tafel slope [24]. To further obtain a clear comparison of the as-prepared catalysts, overpotentials and Tafel slopes required to reach a current density of 10 mA cm−2 were plotted in Figure 3B. Obviously, the overpotential of VMoCoCOx@NF was only 37 mV, which was much lower than that of VCoCOx@NF (76 mV), V2O5@NF (209 mV), CoCO3@NF (235 mV), and Mo(CO3)2@NF (189 mV). Nevertheless, the overpotential of commercial Pt/C@NF (26 mV) was found to be lower than that of VMoCoCOx@NF. This fact can only be explained by the highest HER activity of VMoCoCOx@NF among the as-prepared catalysts, except for Pt/C@NF. Additionally, the Tafel slopes for VMoCoCOx@NF, Pt/C@NF, VCoCOx@NF, V2O5@NF, CoCO3@NF, and Mo(CO3)2@NF were 45, 32, 82, 102, 112, and 107 mV dec−1, respectively. Figure S7 provides corresponding Tafel plots. We fitted the linear regions in the LSV curves of all of the samples to the Tafel equation to measure the corresponding Tafel slopes. As revealed in Figure S7A, the VMoCoCOx@NF exhibited a low Tafel slope of 45 mV dec−1 at the current density of 10 mA cm−2 during HER. The lower the Tafel slope, the faster the HER catalytic kinetics will be. The lowest Tafel slope of VMoCoCOx@NF among the catalysts also further proved that the grafting of V/Mo clustered nanosheets favorably promoted the HER catalytic kinetics. The electrochemical double layer capacitance (Cdl) test was carried out to further evaluate the HER activity [23]. The cyclic voltammetry (CV) curves of VMoCoCOx@NF, VCoCOx@NF, and V2O5@NF recorded at various scan rates (10, 20, 40, 60, 80, 100, and 120 mV s−1) are displayed in Figure S8A–C. For better interpretation and comparison of electrochemical surface areas (ECSAs), we measured the Cdl behavior of the catalysts at the same scan rates. At each scan rate, the VMoCoCOx@NF catalyst provided much higher anodic (ia) and cathodic (ic) current densities (CV areas) than those of the other catalysts, consistent with a much greater active surface area. From these linear plots in Figure 3C or Figure S8D, we obtained the Cdl values of VMoCoCOx@NF, VCoCOx@NF, and V2O5@NF, and the bare Ni foam of 15.4, 6.2, 2.7, and 0.249 mF cm−2, respectively. The considerably higher Cdl of the VMoCoCOx@NF suggested that it had the highest ECSA and, correspondingly, the highest surface roughness when compared with those of the other catalysts. Furthermore, an electrochemical impedance spectroscopy (EIS) plot is collected to assess the charge transfer features of materials between the electrode and electrolyte interface. As illustrated in Figure S6, the VMoCoCOx@NF electrode exhibits a smaller radius than those of VCoCOx@NF, CoCOx@NF, and V2O5@NF, implying the smallest charge transfer resistance (Rct) for the HER process, suggesting it has excellent electron transfer. Above all, the Mo modification accelerates charge transfer in the HER process, and the higher HER activity of the VMoCoCOx@NF electrode can be explained by the smallest Tafel slope and minimal Rct.
Stability is another key parameter for the catalyst in practical applications. The stability of VMoCoCOx@NF was measured by CV and chronoamperometry (i–t) methods for HER at a constant current density of 10 mA cm−2 [10]. As displayed in Figure 3D, the stability of VMoCoCOx@NF was evidenced for more than 100 h, with the current density at about 93.7% of its original value, indicating an outstanding stability of the use of VMoCoCOx@NF. The inset in Figure 3D presents the LSV curves of the VMoCoCOx@NF electrode initially and after 10,000 CV cycles, which basically overlap, proving the catalyst has good durability. In addition, to compare the catalytic performances of VMoCoCOx@NF with different mole ratios, the HER polarization curves of VMo0.5Co0.5COx@NF, VMo0.6Co0.4COx@NF, VMo0.4Co0.6COx@NF, and VMo0.3Co0.7COx@NF were displayed in Figure S9A. Among them, VMo0.5Co0.5COx@NF needed a small overpotential to reach 10 mA cm−2, indicating that its catalytic performance was greatly improved after the introduction of V, Mo, and Co with the optimized mole ratio 1:0.5:0.5. These results confirmed that the three elements incorporated with the optimized mole ratio of 1:0.5:0.5 effectively improved the HER electrocatalytic activity of the catalyst.

2.3. Electrocatalytic OER Measurement

Similarly, VMoCoCOx@NF exhibited better OER performance in an O2-saturated 1.0 M KOH electrolyte at a current density of 10 mA cm−2. To further ascertain the codoping contribution of trimetallic V, Mo, and Co loaded on the bare NF substrate, the OER performance was also assessed using the standard three-electrode system. As depicted in Figure 4A, VMoCoCOx@NF revealed the best OER performance with the lowest overpotential. In Figure 4B and Figure S7B, the overpotential of VMoCoCOx@NF only reached 264 mV, which was far less than those for VCoCO3@NF (298 mV), V2O5@NF (345 mV), CoCO3@NF (361 mV), and Mo(CO3)2@NF (354 mV), and close to RuO2@NF (252 mV). The Tafel slope value obtained for VMoCoCOx@NF (67 mV dec−2) was also much lower than those of VCoCO3@NF (85 mV dec−1), V2O5@NF (93 mV dec−1), CoCO3@NF (98 mV dec−1), and Mo(CO3)2@NF (101 mV dec−1). The lower value of the Tafel slope for VMoCoCOx@NF also reveals the better OER performance of VMoCoCOx@NF, in good agreement with the LSV data. It clearly revealed the superior activity of VMoCoCOx with Ni foam support [25]. The porous NF substrate plays an important role in improving electrocatalytic properties by enlarging the active surface area, accelerating charge/electron transfer, and providing excellent electrical conductivity. The active catalyst loaded on the NF substrate also improved the catalyst–substrate contact for effective charge/electron transfer and rapid release of H2/O2 during water electrolysis.
Furthermore, as shown in Figure 4C, VMoCoCOx@NF displayed an Rct of 0.66 Ω, which was smaller than those of VCoCOx@NF (19.8 Ω), V2O5@NF (27.8 Ω), and NF (>60 Ω). VMoCoCOx@NF featured the smallest semicircle diameter with a minimum Rct among the other catalysts, indicating that the dopants reduced the Rct, enhanced the electrocatalytic activity, and benefited the OER catalytic kinetics. The V, Mo, and Co elements acted as the exposed active sites for the OER reaction as well. Also, the active catalyst loaded on NF improved catalyst–substrate contact, enabling efficient charge/electron transfer during water electrolysis, thereby showing its better electronic conductivity and faster OER catalytic kinetics [15]. The Rct was also consistent with the Tafel slopes and overpotentials. Based on the above, in situ growth of VMoCoCOx@NF also contributed to a decrease in Rct. The trimetallic doping of VMoCoCOx@NF was critical to the enhancement of the OER catalytic kinetics.
The stability of VMoCoCOx@NF during OER was assessed using chronoamperometry (it) and CV tests. As shown in Figure 4D, its chronoamperometric curve displayed that the constant current density of 10 mA cm−2 was barely degenerated after the 100 h test. The inset in Figure 4D depicts that the LSV curves of the as-prepared electrode almost overlap with those of the initial and after 100 h tests, which proves the sample has outstanding stability. The aforementioned comparative study showed that the VMoCoCOx@NF had a smaller overpotential and Tafel slope than those of the compared catalysts, thus confirming its superior OER performance (Table S8). Moreover, as shown in Figure S9B, VMo0.5Co0.5COx@NF manifested better OER activity with a smaller overpotential, which was clearly less than those of VMo0.6Co0.4COx@NF, VMo0.4Co0.6COx@NF, and VMo0.3Co0.7COx@NF, implying a significantly enhanced electrocatalytic performance towards OER after the incorporation of trimetallic V, Mo, and Co with an optimized mole ratio of 1:0.5:0.5. The optimized experimental results were also consistent with the RSM simulation.

2.4. Electrocatalytic Performance for Overall Water Splitting

Electrocatalytic water splitting refers to the process by which water molecules are cleaved into H2 and O2 by catalyzing the separation of H2O through the catalyst electrodes in an electrolytic cell (2H2O → 2H2↑ + O2↑). Considering the excellent OER and HER activities of VMoCoCOx@NF, by using the as-prepared catalyst serving as both anode (4OH → 2H2O + O2↑ + 4e) and cathode (2H2O + 2e→ 2OH+ H2↑) for overall water splitting, a symmetric two-electrode water electrolyzer was constructed in a 1.0 M KOH electrolyte.
To investigate the influence of VMoCoCOx@NF with different V/Mo/Co mole ratios on the alkaline HER and OER, the mole ratios of metal precursors were optimized by RSM. Herein, single-factor experiments confirmed that the mole amounts of Mo and Co were sensitive factors for the catalyst formulation, and these sensitive factors were optimized by RSM. The VMo0.6Co0.4COx@NF, VMo0.5Co0.5COx@NF, VMo0.4Co0.6COx@NF, and VMo0.3Co0.7COx@NF were selected as target electrocatalysts (Table S1). When V was set to 1 mole, Mo mole (χ1) and Co mole (χ2) were selected as independent variables (Table S2). Mo/Co mole ratio, O2, and H2 production were chosen as responses. Figures S1 and S2 present the single and combining effects of independent variables on all responses. The 3D response surface contour plots (3D-RSCP) were generated accordingly. The adequacy of the model was rigorously evaluated by analysis of variance (ANOVA) with the corresponding results listed in Tables S4 and S5. The coefficients of determination R-squared (R2), adjusted R2 (adj. R2), and predicated R2 (pred. R2) were used to evaluate the fitting model. The difference between adj. R2 and pred. R2 was less than 0.2, which indicated that the model fitted all data (Figure S2 and Table S3 for details).
The generated volumes of H2 and O2 were the response values of optimization with help from RSM. In Table S1, the VMoCoCOx@NF fabricated with NH4VO3: (NH4)2MoO4:CoCl2 mole ratio of 1:0.5:0.5 showed maximum yields of H2 (96.3%) and O2 (92.5%), which exhibited excellent catalytic activity and should be a most promising bifunctional catalyst. Figure S2 illustrates a good match between predicted and actual generated volumes of O2 and H2. The predicted yields of H2 (94.0%) and O2 (90.2%) were similar to the actual yields of H2 (96.0%) and O2 (91.6%) (Table S6). The comparison points were evenly distributed on both sides of the fitted line. On the other hand, LSV measurement is one of the most common methods for rapid assessment of electrocatalyst activity. The LSV polarization curves of the optimized VMo0.5Co0.5COx@NF for both HER and OER showed lower overpotentials than the other catalysts, as depicted in Figure S9. Therefore, Equations (S1) and (S2) in the SI provide a good fit to the data and can be used to predict the uninvestigated experimental variables. The prediction model after the ANOVA on the responses predicted the optimal V:Mo:Co mole ratio of 1.0:0.538:0.483. The VMoCoCOx@NF with an optimal mole ratio of 1:0.5:0.5 revealed outstanding bifunctional electrocatalytic activity for overall water splitting by achieving 10 and 100 mA cm−2 at cell voltages of 1.54 V and 1.61 V in 1.0 M KOH electrolyte. Inspired by the superb OER and HER performances of VMoCoCOx@NF in alkaline media, the catalyst served as both anode and cathode by applying electric voltage to an alkaline electrolyzer. Numerous tiny bubbles of H2 and O2 were discharged from the electrode surface. Moreover, to investigate the electrocatalytic efficiency of the VMoCoCOx@NF||VMoCoCOx@NF electrolyzer, the generated H2 and O2 gases were monitored using gas chromatography (GC). The millimoles of O2 and H2 generated from the water splitting at 5 min intervals are intuitively observed in Figure 5A. As expected, after 30 min of water electrolysis, 0.298 millimoles of O2 and 0.694 millimoles of H2 were collected at the anode and the cathode, respectively. Concomitantly, the evolution millimolar ratio of the generated H2 and O2 bubbles was close to 2:1, which was consistent with the theoretically calculated value produced by water electrolysis. The agreement between experimental and theoretical values indicated the gas-generated efficiencies were 96.3% for H2 and 92.5% for O2, respectively, which were closer to 100%, proving that VMoCoCOx@NF was a promising bifunctional catalyst (Table S1). The contribution of the NF-supported VMoCoCOx@NF catalyst was also verified. Table S8 summarizes and compares the electrocatalytic activity (HER, OER, and overall water splitting) of the VMoCoCOx@NF catalyst with those of recently reported NF-supported transition metallic electrocatalysts in 0.1 M KOH electrolyte, from which it is found that the HER/OER results as well as the water-splitting results of VMoCoCOx@NF are better than those of the recently reported bifunctional electrocatalysts [17,22,26,27,28,29,30,31,32]. The overall water-splitting experiment was performed in a two-electrode electrolyzer consisting of a VMoCoCOx@NF catalyst as both the positive and negative electrodes. Notably, the water-splitting penitential of the VMoCoCOx@NF||VMoCoCOx@NF alkaline water electrolyzer was lower than that of the others, as displayed in Figure 5B. It needs a voltage of 1.54 V to afford a water-splitting current density of 10 mA cm−2 with tiny bubbles of H2 and O2 evolution on both electrodes. The developed catalyst had a promising application in place of the noble metal electrocatalyst, which might be related to the following factors [33,34]. First, the clustered VMoCoCOx nanosheet structure was verified by FESEM analysis. The catalyst could expose more active sites, thereby accelerating water electrolysis. Second, the multimetallic active sites facilitate the bifunctional electrocatalytic behavior due to heteroatom doping of V/Mo/Co. Third, discontinuous crystal fringes enlarged the surface area on VMoCoCOx@NF, favoring the adsorption, activation, and dissociation of water molecules. Fourth, porous NF structure also greatly facilitates the generation and diffusion of H2 and O2 bubbles and effectively enhances water electrolytic efficiency. Finally, the integrated catalytic material design reduces the charge/electron transfer resistance during the water electrolysis. Herein, a novel electrocatalyst with V and Mo codoping for water splitting was first developed. The excellent electrocatalytic performance of VMoCoCOx@NF was assigned to the synergistic effect of 2p, 3d, and 4d vacancies and multimetallic codoping, including increased valence states and enhanced catalytic kinetics. The electronic structure benefits from the assembled forest-like nanosheet architecture with more active sites, a short charge/electron transfer path, and easy infiltration of electrolytes [35].
Figure 6A displays the LSV curves of VMoCoCOx@NF||VMoCoCOx@NF, Pt/C@NF||RuO2@NF, and NF||NF electrolyzers. The VMoCoCOx@NF||VMoCoCOx@NF alkaline water electrolyzer required a cell voltage of 1.54 V to achieve a current density of 10 mA cm−2, which was lower than those of the commercial Pt/C@NF||RuO2@NF (1.56 V) and NF||NF (1.96 V) catalysts. The better performance outperformed most previously reported catalysts (Table S8) [5]. More significantly, the catalytic activity of VMoCoCOx@NF was increased when V and Mo were codoped into CoCOx, demonstrating that the synergetic effect of V and Mo was mainly responsible for improving the bifunctional performance. Figure 6B shows the VMoCoCOx@NF||VMoCoCOx@NF alkaline water electrolyzer operating at current densities of 10 and 100 mA cm−2, the voltages of 1.54 V and 1.61 V, respectively, which are superior to Pt/C@NF||RuO2@NF electrolyzer (voltages of 1.56 V and 1.67 V, respectively). Therefore, the VMoCoCOx@NF catalyst possesses promising applications for water electrolysis in reality.
The durability of the electrode is a key parameter to evaluate the performance of an electrocatalyst. The catalyst was subjected to continuous CV for durability testing, and the polarization curves were plotted after 10,000 CV cycles. Figure 6C shows almost similar polarization curves before and after the 120 h test, which is conducive to the long-term stable use of the catalyst. The inset in Figure 6C shows the overall water electrolytic cell device using VMoCoCOx@NF as both electrodes. In addition to activity, long-term durability is another critical aspect of catalysts, especially for nanostructural catalysts. Figure 3D and Figure 4D show the superior stabilities for HER and OER at a constant current density of 10 mA cm−2. Indeed, the stability measured at a constant current density is more reliable. Considering the large mass loading on the NF substrate, the cyclic stability of over-water splitting should be measured at a large current density. In Figure 6D, when constant current densities of 10 mA cm−2 and 100 mA cm−2 were applied for overall water splitting, the vt curves of the fabricated VMoCoCOx@NF required cell voltages of 1.54 V and 1.61 V. The potentials exhibited negligible changes over 120 h tests with robust stability. As displayed in Figure S10, there were no obvious changes in the structure (XRD), morphology (SEM), or chemical state (XPS) of the sample after stability testing, indicating that the catalyst had very good performance.
Interestingly, as shown in Figure 3B and Figure 4B, the VMoCoCOx@NF catalyst exhibited the highest HER and OER activities among the synthesized samples, except for the commercial Pt/C and RuO2. Nevertheless, the overall water splitting experiments were performed in a two-electrode configuration composed of VMoCoCOx@NF catalyst as anode and cathode. Figure 6A,B displays that the water-splitting potential of the VMoCoCOx@NF||VMoCoCOx@NF cell (1.54 V) is lower than that of the Pt/C@NF||RuO2@NF cell (1.59 V). This fact can be explained by the high over-water-splitting activity of VMoCoCOx@NF||VMoCoCOx@NF cell.

2.5. DFT Theoretical Calculation

The DFT theoretical calculation in this work was theoretically performed to elucidate the origin of the improved electrocatalytic performance of VMoCoCOx@NF in alkaline media. In order to understand its contribution to the electrocatalytic activity, in the computation hydrogenation electrode (CHE) modeling, a VMoCoCOx@NF surface slab with six-atomic hexagonal clustered nanosheets was first constructed to study the reaction intermediate steps involved in both HER and OER. As revealed in Figure 7A, the optimized geometry structures of OH*, H*, and H2O* intermediates on the active sites of VMoCoCOx were observed for HER. Fortunately, the H* spontaneously moved between the interfacial V and Mo, suggesting that the doping metallic atoms could provide active sites. Meanwhile, the structural features and the Gibbs free energy (ΔG) profiles of VMoCoCOx@NF were verified by constructing the correlative theoretical model. The alkaline catalytic process proceeds through either the Volmer–Heyrovrosky or Volmer–Tafel reactions (SI). The alkaline HER mechanism of VMoCoCOx@NF was proposed through the DFT calculation: VMoCoCOx@NF served to favor the adsorption, activation, and dissociation of H2O (* + H2O → *H2O → H* + *OH) to produce H* and *OH. The clustered VMoCoCOx nanosheets provided the active sites for H2 generation (*OH + e → OH; H* + H*→ H2↑). To examine the H* adsorption/H2 desorption step (the Heyrovsky step), the ΔGH* values for all possible active sites were calculated [36]. Both pathways involve the adsorption of H2O, electrochemical reduction into adsorbed H atoms and OH, and desorption of H2.
In an alkaline environment, OER has occurred on the anode. During the adsorption process, the velocity of the adsorbent decreases slowly; it finally stops on the surface slab of the catalyst, and part of the energy is released, which is the adsorption energy. The more negative the adsorption energy, the faster the charge/electron transfer from active sites to H2O molecules. The adsorption energy is the critical factor for validating these catalysts’ theoretical activity [20,22]. The V, Mo, and Co atoms of VMoCoCOx@NF are the exact active sites for OER. The adsorption energies at the active sites of V, Mo, and Co of VMoCoCOx@NF are −6.88 eV, −7.39 eV, and −5.43 eV, respectively. In this case, the geometric structures of OER intermediate adsorption on VMoCoCOx@NF are depicted in Figure 7B. The OER catalytic activity can follow four steps: (1)–(4).
M* + OH→M-OH* + e  [OH first adsorption step (1)]
M-OH* + OH →M-O* + H2O + e  [O-H bond breaking step (2)]
M-O* + OH →M-OOH* + e  [OH second adsorption step (3)]
M-OOH* + OH →M* + O2↑ + H2O + e  [oxygen forming step (4)]
where M is the active site of V, Mo, and Co elements. First, OH in the electrolyte is adsorbed by the catalytically excited active sites (M*) to produce M-OH*. Second, the produced M-OH* reacts with OH to form M-O*, and the produced M-O* further reacts with OH to produce intermediate M-OOH*. Finally, the intermediate M-OOH* reacts with OH to generate O2 [37].
For the HER, Figure 7C displays the reaction intermediate steps and corresponding ΔGH* values involved in the HER process on the catalyst, which unveils the origination of electrocatalytic HER activity. The RDS for H2 generation was the conversion of H* to H2, which was the Tafel step. These profiles compared the four computed ΔGH* values for different possible adsorption sites of CoCOx@NF, VMoCOx@NF, VMoCoCOx@NF, and VCoCOx@NF, respectively. According to the ΔGH* values of HER on the as-prepared catalysts, RDS was the step involved in the formation of H*. The ΔGH* of RDS were −0.04 eV, 0.07, 0.23, and 0.37 eV for the HER activities of VMoCoCOx@NF, VCoCOx@NF, VMoCO3@NF, and CoCO3@NF, respectively. Similarly, V doping of CoCO3@NF significantly reduced the ΔGH* from 0.37 eV (CoCO3@NF) to 0.07 eV (VCoCOx@NF). The lower the ΔGH* is, the higher the HER activity. The dissociated H* were adsorbed by the V and Mo atoms on the surface of the VMoCoCOx@NF electrode. The ΔGH* of VMoCoCOx@NF for water dissociation to generate H2 was −0.04 eV, which was close to the optimum value ΔGH* = 0 for an ideal HER electrocatalyst, and the H2O could be effectively activated on its surface. The above results suggested that the ΔGH* value of the VMoCoCOx@NF catalyst was lower than that of the VCoCOx@NF, VMoCO3@NF, and CoCO3@NF catalysts. The significantly decreased ΔGH* value demonstrated that the heterointerface effect between VMoCoCOx and NF also speeded up the H2 release on the developed catalyst, which further promoted HER activity in alkaline media. The DFT results revealed that the VMoCoCOx@NF exhibited excellent water-splitting capability. The structure of VMoCoCOx@NF was moderated by the dopants of V and Mo, reducing its adsorption of H2O and increasing the possibility of H2O dissociation. The theoretical calculation findings are basically consistent with the experimental data.
However, a theoretical study depicted the OER process in alkaline electrolytes. Water molecules underwent a complex four-step electron “cross-mountain” process (* → *OH → *O → *OOH → * + O2). The OER activity was related to the corresponding ΔG of *, *OH, *O, and *OOH intermediates. The total ΔG barrier of O2 generated from water electrolysis was 4.92 eV. In theoretical conditions, the ΔG value for each step of water electrolysis was 1.23 eV, but the actual OER process was usually complex, and the ΔG was different for each step. Actually, the alkaline OER process (4OH → O2↑ + 2H2O + 4e) can be described by the four-step water associative mechanism [4].
The OER activity of the catalyst was significantly related to the characteristics of the surface metal sites. For the DFT calculations, the adopted active sites were V, Mo, and Co on VMoCoCOx@NF. In Figure 7D, the affinity of VMoCoCOx@NF with *O in water was slightly higher; step (3) was the RDS. The discussion concerning the OER process showed the RDS was found to be deprotonation of absorbed *O, which was in good accordance with Ding’s finding [36]. The decreased ΔG of the RDS suggested a higher intrinsic OER activity.
In Figure 7D, the ΔG values of different steps for *, OH*, O*, and OOH* were calculated to reveal the origin of OER activity. The star * is the active site. The theoretical overpotential of OER on VMoCoCOx@NF was determined by the following equation:
ηOER = {max[ΔG1, ΔG2, ΔG3, ΔG4]/e} − 1.23 [V]
where ΔGn (n = 1–4) is the ΔG along the OER process (4.92 eV/4 = 1.23 eV). DFT calculation indicated that VMoCoCOx@NF processed a smaller ΔG than those of VMoCOx@NF and VCoCOx@NF for PDS in the four-step OER reactions, confirming the positive effect of VMoCoCOx@NF on catalyzing OER.
The ΔG diagrams of VMoCOx@NF, VMoCoCOx@NF, and VCoCOx@NF calculated for the OER are shown in Figure 7D. What is more, the value of ΔGO* − ΔGOH* is often chosen as the OER descriptor [4]. Herein, the RDS was the deprotonation step (HO* → O*+ H+ + e) with ΔG20 of 0.60 eV, 0.62 eV, and 0.64 eV for VMoCoCOx@NF (1.97 eV − 1.37 eV = 0.60 eV), VMoCOx@NF (2.27 − 1.65 = 0.62 eV), and VCoCOx@NF (1.82 − 1.18 = 0.64 eV), respectively. Also, as Figure 7D suggested, the *OH-M intermediate was thermodynamically slightly more stable than *OH-M because of its lower ΔG20. Therefore, the affinity of VMoCoCOx@NF with *OH in water was also slightly higher, implying higher surface wettability, which was further verified by the HRTEM test (0.190 nM, 0.285 nM). The codoped V and Mo almost directly took part in the OER procedure. According to DFT computations, the third step is the key process. The electronic modulation from co-doping heteroatoms into composites enhances the electrical conductivity owing to the synergistic effect among V2O5, MoO3, and/or CoCO3 that accelerates charge/electron transfer and provides more additional active sites.
The presence of V, Mo, and Co elements substantiated the trimetallic doping of VMoCoCOx@NF, which synergistically improved the active sites of the synthesized hybrid catalysts. For the OER and HER activities, the V, Mo, and Co ions in VMoCoCOx@NF were considered active sites among V2O5, MoO3, and/or CoCO3 for each reaction. The construction of multiple heterointerfaces in the VMoCoCOx@NF hybrid catalyst can synergistically improve the charge/electron transfer kinetics and expose more active sites for overall water splitting.
The factors leading to the good bifunctional HER/OER activity of the fabricated catalyst are to facilitate charge-transfer kinetics and enrich more electrically connected active sites. Herein, V, Mo, and Co dopants synergistically enrich more additional active sites on VMoCoCOx@NF, and NF acts as a conducting substrate. The advantage of this heterostructure is conveyed by its fabrication with enriched active sites acting as bridged sites from dopants, which facilitates charge-transfer behavior and ameliorates the catalytic activity of the catalyst for overall water splitting.

3. Experimental

Additional experimental details for “Chemicals and Materials” and “Density Functional Theory (DFT)” were provided in the SI.

3.1. Preparation of Electrocatalysts

Scheme 1 illustrates the fabrication procedures of VMoCoCOx@NF nanoarrays [17]. The Mo and V-codoped CoCOx nanosheets loaded on the NF substrate were designed to be fabricated as a highly efficient bifunctional catalyst by facile hydrothermal deposition. Prior to the preparation, a piece of NF with the size of 3 × 1 × 0.1 cm3 was treated with 1.0 M HCl to remove surface impurities by ultrasonication for 15 min, washing with deionized (DI) water and ethanol consecutively for 10 min, followed by drying in a vacuum oven for 8 h. Specifically, 15.0 mL DI water, 0.18 g NH4VO3, 0.15 g (NH4)2MoO4, 0.18 g CoCl2⋅6H2O, 0.90 g NaOH, and 1.00 g CH4N2O were successively added to a 100 mL beaker. Afterward, the treated NF was immersed in the beaker. Subsequently, the mixed solution was ultrasonicated for 5 min, followed by adding 1.00 g of (NH4)2S2O8. After stirring the solution at 25 °C for 12 h, the above solution with NF was transferred to a 25 mL Teflon-lined stainless-steel autoclave and maintained at 150 °C for 10 h. Finally, the as-prepared catalyst was washed with DI water and ethanol in sequence and dried in a vacuum oven at 60 °C for 10 h. Two control samples of VCoCOx@NF and VMoCOx@NF were also fabricated under the same conditions for comparison.

3.2. Characterization

The morphology, structure, and elemental information of the prepared catalysts were measured by FESEM-EDX (JEOL JSM-7600F, JEOL, Tokyo, Japan, operating at 5 kV). The crystalline/amorphous phase and lattice spacing were identified using HRTEM (JEM-2100, JEOL, Tokyo, Japan, operating at 200 kV). HRTEM simulation was performed using JEMS 2014 software. The crystallographic identification was acquired by XRD (Ultima IV, Rigaku, Tokyo, Japan, Cu Kα radiation, λ = 1.5406 Å). The test conditions were 40 kV, 40 mA, and a scan rate of 1 × 10−2 degree/second in the 2θ range from 10° to 80°. The surface elemental composition and valence states were investigated by XPS (PHI 5000 VersaProbe, ULVaC-Phi, Kanagawa, Japan, Al Kα radiation, = 1486.6 eV), and data were fitted by the XPSPEAK 41 software. FTIR analysis was conducted on an FTIR spectrometer (VERTEX 80V, Bruker, Ettlingen, Germany) from 450 to 4000 cm−1 at 0.5 cm−1 spectral resolution. An inductively coupled plasma–optical emission spectrometer (ICP-OES) (Avio 220, PerkinElmer, Waltham, MA, USA, operating under 15.0 L min−1 coolant Ar gas, 0.9 L min−1 plasma Ar gas, and 0.9 L min−1 nebulizer Ar gas) was used to detect the element molar ratios.

3.3. Electrochemical Measurements

Electrochemical measurements were conducted on an electrochemical workstation (CHI 760E, Chenhua Instrument Company, Shanghai, China) connected to a conventional three-electrode system in 1.0 M KOH electrolyte at ambient temperature (25 °C). The fabricated VMoCoCOx@NF catalyst was tailored as 1 × 1 cm2 for direct use as the working electrode. A graphite plate and normal calomel (Hg/HgO) were used as the counter electrode and reference electrode, respectively. The measured potentials (EHg/HgO) were calibrated to reversible hydrogen electrode potentials (ERHE) according to the Nernst equation: ERHE = EHg/HgO + 0.098 V + 0.059 × pH. All data presented were corrected for ohmic resistance. The HER/OER activities and Tafel slopes of the tested catalysts were compared, respectively. The HER and OER polarization tests were performed in N2-saturated and O2-saturated electrolytes, respectively. The LSV polarization curves of both OER and HER were recorded at a scan rate of 5 mV s−1. The EIS was measured from 100 kHz to 0.01 Hz with an amplitude of 5 mV. The cycling stability of VMoCoCOx@NF was measured at a scan rate of 50 mV s−1 by CV for 10,000 CV cycles. The acquisition of the Tafel equation, η = a + b log|j|, is based on the LSV calculation, where b is the Tafel slope and j is the current density. The ECSA of the electrode was evaluated using the electrochemical Cdl and the specific capacitance (Cs) (ECSA = Cdl/Cs). Cs is 0.04 mF cm−2 for transition metal-based materials in alkaline media [35]. Cdl was measured by a CV test under scan rates ranging from 10 to 120 mV s−1. The two-electrode configuration was used to perform overall water splitting. The permanency of the optimized VMoCoCOx@NF catalyst toward HER, OER, and overall water splitting was evaluated using the chronoamperometry (jt) or chronopotentiometry (vt) methods. The generated volumes of H2 and O2 gas were measured using GC (7890A, Agilent Technologies, Santa Clara, CA, USA) (GB/T 28124-2011) [38]. Besides VMoCoCOx@NF, VCoCOx, VMoCOx, CoCO3, V2O5, commercial Pt/C, and RuO2 loaded on NF were also used as working electrodes for comparison.

3.4. Density Functional Theory (DFT) Calculations

All spin-polarized DFT calculations were completed using the Vienna Ab initio software package (VASP5.4.4) [35]. The projector augmented wave (PAW) method was used to describe the electron–ion interaction, with 1s of H, 2s2p of C and O, 3d4s of V and Co, and 4d5s of Mo considered valence electrons. The electron–electron exchange-correlation functions were calculated by Perdew–Burke–Ernzerhof (PBE) in the generalized gradient approximation (GGA) for structure and energy optimization [39]. For an accurate description of V (3d), Co (3d), and Mo (4d) electrons, Hubbard effective terms Ueff (Co) = 3.32 eV, Ueff (Mo) = 4.38 eV, and Ueff (V) = 3.4 eV were added to the PBE functional via the rotationally invariant approach. The Ueff values were obtained from the Materials Project database [40]. For all periodic structures, the k point of the Brillouin zone based on monkhorst-pack grid mesh was 2 × 1 × 1, the truncation energy was 500 eV, and the convergence criteria of force and energy were 10−2 eV Å−1 and 10−5 eV, respectively [41].
However, the VMoCoCOx (112) plane reproduces the “distorted hexagonal clustered nanosheets” structure very well. The optimized bulk face-centered cubic V, Mo, and Co structure has lattice constants of 3.030 Å, 3.147 Å, and 2.507 Å, which are much closer to their experimental lattice constants (3.046 Å, 3.123 Å, and 2.531 Å). The lattice constants for optimized bulk hexagonal VMoCoCOx are a = 4.577 Å and c = 4.211 Å (experimental values: a = 4.562 Å and c = 4.127 Å). The optimized bulk structure was used to construct the surface slab model [42]. For VMoCoCOx@NF, a V- and Mo-terminated surface (112) slab of 3 layers was constructed. The VMoCoCOx hexagonal clustered nanosheets loaded on the top surface of the NF substrate (111), with V and Mo terminated, interact with the NF (111) surface to make an effective catalytic interface. Therefore, the VMoCoCOx nanosheets in the hybrid model experienced compressive strain from the NF substrate. The two sides of the clustered hexagonal nanosheets correspond to the V- and Mo-terminated surfaces of VMoCoCOx (112). In Figure S11A, the structural characteristics for constructing VMoCoCOx@NF were explored by employing the methods of model chemistry, including ab initio computation and DFT calculation. As illustrated in the VMoCoCOx@NF structure diagram, except for V, Mo, and Co, the CO32− anions composed of C and O should intersperse in the structure constructed by V, Mo, and Co.
In the H*Slab model, the adsorption energy of hydrogen (ΔEH*) was expressed as ΔEH* = EH*slab − (Eslab + 0.5EH2). Meanwhile, a catalyst with ΔG = 0, or 1.23 eV, is an excellent candidate for OER. The ΔG correlations were calculated by zero-point energy and entropic correction [43]. The climbing image nudged elastic band (CI-NEB) and the dimer approach were used to simulate the reaction pathway.

4. Conclusions

In this work, the molybdenum and vanadium codoped cobalt carbonate nanosheets loaded on nickel foam (VMoCoCOx@NF) were fabricated, and the mole ratio of V/Mo/Co was optimized by response surface methodology. The optimized VMoCoCOx@NF catalyst displays superior HER and OER performance with lower overpotentials and smaller Tafel slopes, along with robust stability in the 1.0 M KOH electrolyte, which is associated with the synergy of V2O5, MoO3, and CoCO3 and heterostructural features. Also, an alkaline water electrolyzer with VMoCoCOx@NF as both the cathode and anode enables current densities of 10 and 100 mA cm−2 with voltages of 1.54 V and 1.61 V, respectively. Through detailed investigation, the working mechanism of the yielded advanced HER and OER activity was clarified via the DFT simulations. More importantly, this work offers an excellent and highly cost-effective self-supported transition metal-based bifunctional catalyst material for practical use in clean energy H2 generation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29153591/s1. Optimize mole ratios of the fabricated electrocatalysts with response surface methodology (RSM); Detailed DFT calculations; Figure S1: 3-D response surface plots and contour plots for the as-prepared catalysts with different mole ratios of Mo and Co; Figure S2: Fit plots for predicted and actual values; Figure S3: XRD patterns of VMoCoCOx@NF, VCoCOx@NF, CoCOx@NF and V2O5@NF on the same coordinates; Figure S4: (A) HRTEM image of VMoCoCOx@NF. (B) The lattice spacing distributions corresponding to 0.261 nm for the (310) plane and (C) 0.190 nm for the (112) plane; Figure S5: The FTIR spectra of VMoCoCOx@NF, VCoCOx@NF and V2O5@NF; Figure S6: Nyquist plots of HER for the four catalytic electrodes; Figure S7: Comparison of Tafel slopes of the catalysts selected for HER (A) at a constant j = −10 mA cm−2 and OER (B) at a constant j = 10 mA cm−2; Figure S8: The original CV curves of VMoCoCOx@NF (A), VCoCOx@NF (B) and V2O5@NF (C) recorded at various scan rates in the non-faradaic region for the estimation of Cdl. (D) Estimation of Cdl values of VMoCoCOx@NF, VCoCOx@NF, CoCOx@NF and V2O5@NF by plotting the ratios of catalyst ΔJ to scan rate.; Figure S9: The polarization curves of VMo0.5Co0.5COx@NF, VMo0.6Co0.4COx@NF, VMo0.4Co0.6COx@NF, and VMo0.3Co0.7COx@NF for HER (A) and OER (B); Figure S10: The structure (XRD)(A), morphology (SEM)(B) and chemical state (XPS)(C) of the sample before and after stability testing; Figure S11: (A) The crystal plane of VMoCoCOx (112), (B) Crystal structure of NF (111) in top view. Table S1. Optimized the catalysts with different mole ratios of V, Mo and Co based on the experimental production of O2 and H2. Table S2. Levels of independent variables for the response surface design for VMoCoCOx@NF with different mole ratios of Mo and Co. Table S3. Response surface optimal (custom) design and the experimental and predicted production of O2 and H2 by VMoCoCOx@NF with different ratios of Mo and Co. Table S4. ANOVA for response surface quadratic model for O2 production of VMoCoCOx@NF with different ratios of Mo and Co. Table S5. ANOVA for response surface quadratic model for H2 production of VMoCoCOx@NF with different ratios of Mo and Co. Table S6. Optimum values for the ratio of Mo and Co for VMoCoCOx@NF and the production of O2 and H2. Table S7. XRD data of the target materials. Table S8. Comparison of the recently reported NF-supported transition metallic electrocatalysts for the electrocatalytic activities (HER, OER, and overall water splitting) in 0.1 M KOH electrolyte.

Author Contributions

Investigation, methodology, data curation, formal analysis, visualization, validation, and writing—original draft: W.W.; software and visualization: L.X.; investigation and data curation: R.Y. and P.Y.; supervision, validation, and funding acquisition: J.Z.; supervision, conceptualization, writing—review and editing, resources, validation, and funding acquisition: L.J.; supervision, conceptualization, methodology, validation, writing—review and editing, project administration, and funding acquisition: X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Nature Science Foundations of China (nos. 21775070 and 21673108).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created except for the data in the Supplementary Materials and the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Jamesh, M.I. Recent progress on earth abundant hydrogen evolution reaction and oxygen evolution reaction bifunctional electrocatalyst for overall water splitting in alkaline media. J. Power Sources 2016, 333, 213–236. [Google Scholar] [CrossRef]
  2. Sun, H.M.; Yan, Z.H.; Liu, F.M.; Xu, W.C.; Cheng, F.Y.; Chen, J. Self-supported transition-metal-based electrocatalysts for hydrogen and oxygen evolution. Adv. Mater. 2020, 32, 1806326. [Google Scholar] [CrossRef]
  3. He, T.O.; Wang, W.C.; Shi, F.L.; Yang, X.L.; Li, X.; Wu, J.B.; Yin, Y.D.; Jin, M.S. Mastering the surface strain of platinum catalysts for efficient electrocatalysis. Nature 2021, 598, 76–81. [Google Scholar] [CrossRef]
  4. Zeng, K.; Tian, M.; Chen, X.; Zhang, J.L.; Rummeli, M.H.; Strasserf, P.; Sun, J.Y.; Yang, R.Z. Strong electronic coupling between single Ru atoms and cobalt-vanadium layered double hydroxide harness efficient water splitting. Chem. Eng. J. 2023, 452, 139151. [Google Scholar] [CrossRef]
  5. Yu, Z.Y.; Duan, Y.; Feng, X.Y.; Yu, X.X.; Gao, M.R.; Yu, S.H. Clean and affordable hydrogen fuel from alkaline water splitting: Past, recent progress, and future prospects. Adv. Mater. 2021, 33, 2007100. [Google Scholar] [CrossRef]
  6. Hao, Y.W.; Du, G.H.; Fan, Y.; Jia, L.; Han, D.; Zhao, W.Q.; Su, Q.M.; Ding, S.K.; Xu, B.S. Mo/P Dual-doped Co/Oxygen-deficient Co3O4 core-shell nanorods supported on Ni foam for electrochemical overall water splitting. ACS Appl. Mater. Inter. 2021, 13, 55263–55271. [Google Scholar] [CrossRef]
  7. Rathore, D.; Ghosh, S.; Chowdhury, J.; Pande, S. Fe-doped NiCo2Se4 nanorod arrays as electrocatalysts for overall electrochemical water splitting. ACS Appl. Nano Mater. 2023, 6, 3095–3110. [Google Scholar] [CrossRef]
  8. Zhang, J.H.; Zhang, H.X.; Liu, M.M.; Xu, Q.C.; Jiang, H.; Li, C.Z. Cobalt-stabilized oxygen vacancy of V2O5 nanosheet arrays with delocalized valence electron for alkaline water splitting. Chem. Eng. Sci. 2020, 227, 115915. [Google Scholar] [CrossRef]
  9. Gong, T.Y.; Zhang, J.Y.; Liu, Y.; Hou, L.R.; Deng, J.L.; Yuan, C.Z. Construction of hetero-phase Mo2C-CoO@N-CNFs film as a self-supported bi-functional catalyst towards overall water splitting. Chem. Eng. J. 2023, 451, 139025. [Google Scholar] [CrossRef]
  10. Yu, S.H.; Shao, Y. A review on cobalt-based oxides electrocatalytic materials for electrochemical water splitting. Ionics 2023, 29, 1273–1284. [Google Scholar] [CrossRef]
  11. Kong, Y.X.; Xiong, D.K.; Lu, C.X.; Wang, J.; Liu, T.; Ying, S.L.; Ma, X.H.; Yi, F.Y. Vanadium-based trimetallic metal-organic-framework family as extremely high-performing and ultrastable electrocatalysts for water splitting. ACS Appl. Mater. Inter. 2022, 14, 37804–37813. [Google Scholar] [CrossRef] [PubMed]
  12. Zhou, J.; Yu, L.; Zhu, Q.; Huang, C.; Yu, Y. Defective and ultrathin NiFe LDH nanosheets decorated on V-doped Ni3S2 nanorod arrays: A 3D core–shell electrocatalyst for efficient water oxidation. J. Mater. Chem. A 2019, 7, 18118–18125. [Google Scholar] [CrossRef]
  13. Xiao, Y.; Chen, X.; Li, T.; Mao, Y.; Liu, C.; Chen, Y.; Wang, W. Mo-doped cobalt hydroxide nanosheets coupled with cobalt phosphide nanoarrays as bifunctional catalyst for efficient and high-stability overall water splitting. Int. J. Hydrogen Energ. 2022, 47, 9915–9924. [Google Scholar] [CrossRef]
  14. Yan, X.D.; Tian, L.H.; Atkins, S.; Liu, Y.; Murowchick, J.; Chen, X.B. Converting CoMoO4 into CoO/MoOx for overall water splitting by hydrogenation. ACS Sustain. Chem. Eng. 2016, 4, 3743–3749. [Google Scholar] [CrossRef]
  15. Xu, L.L.; Yang, P.; Ye, R.L.; Wu, X.C.; Tao, Y.R. NiFeLDH/Mo4/3B2−xTz/NF composite electrodes to enhance oxygen evolution performance. J. Mater. Chem. A 2024, 12, 9714–9722. [Google Scholar] [CrossRef]
  16. Ma, H.B.; Yang, X.J.; Wang, Z.L.; Jiang, Q. Engineering the interface of porous CoMoO3 nanosheets with Co3Mo nanoparticles for high-performance electrochemical overall water splitting. J. Mater. Sci. Technol. 2023, 137, 184–192. [Google Scholar] [CrossRef]
  17. Meena, A.; Thangavel, P.; Nissimagoudar, A.S.; Singh, A.N.; Jana, A.; Jeong, D.S.; Im, H.; Kim, K.S. Bifunctional oxovanadate doped cobalt carbonate for high-efficient overall water splitting in alkaline-anion-exchange-membrane water-electrolyzer. Chem. Eng. J. 2022, 430, 132623. [Google Scholar] [CrossRef]
  18. Wang, G.W.; Zhang, G.K.; Ke, X.X.; Chen, X.Y.; Chen, X.; Wang, Y.S.; Huang, G.Y.; Dong, J.C.; Chu, S.Q.; Sui, M.L. Direct synthesis of stable 1T-MoS2 doped with Ni single atoms for water splitting in alkaline media. Small 2022, 18, 2107238. [Google Scholar] [CrossRef]
  19. Bera, K.; Karmakar, A.; Kumaravel, S.; Sankar, S.S.; Madhu, R.; Dhandapani, H.N.; Nagappan, S.; Kundu, S. Vanadium-doped nickel cobalt layered double hydroxide: A high-performance oxygen evolution reaction electrocatalyst in alkaline medium. Inorg. Chem. 2022, 61, 4502–4512. [Google Scholar] [CrossRef]
  20. Chen, G.R.; Xu, J.Q.; Yu, H.J.; Guo, F.; Xie, J.Q.; Wang, Y.Q. Effect of the non-thermal plasma treatment on the structure and SCR activity of vanadium-based catalysts. Chem. Eng. J. 2020, 380, 122286. [Google Scholar] [CrossRef]
  21. Li, B.; Wang, K.F.; Ren, J.X.; Qu, P. NiOOH@cobalt copper carbonate hydroxide nanorods as bifunctional electrocatalysts for highly efficient water and hydrazine oxidation. New J. Chem. 2022, 46, 7615–7625. [Google Scholar] [CrossRef]
  22. Yang, H.; Hu, T.P.; Meng, R.Q.; Guo, L.J. Efficient Mo-Co(OH)2/Co3O4/Ni foam electrocatalyst for overall water splitting. J. Solid State Chem. 2023, 320, 123837. [Google Scholar] [CrossRef]
  23. Saad, A.; Gao, Y.; Owusu, K.A.; Liu, W.; Wu, Y.Y.; Ramiere, A.; Guo, H.C.; Tsiakaras, P.; Cai, X.K. Ternary Mo2NiB2 as a superior bifunctional electrocatalyst for overall water splitting. Small 2022, 18, 2104303. [Google Scholar] [CrossRef] [PubMed]
  24. Gong, Y.Q.; Zhi, Y.; Lin, Y.; Zhou, T.; Li, J.H.; Jiao, F.X.; Wang, W.F. Controlled synthesis of bifunctional particle-like Mo/Mn-NixSy/NF electrocatalyst for highly efficient overall water splitting. Dalton Trans. 2019, 20, 6718–6729. [Google Scholar] [CrossRef] [PubMed]
  25. Narasimman, R.; Jnanapriya, G.; Sujatha, S.; Ilangovan, S.A. Vertically aligned nickel-iron-phosphide nanosheets on three dimensional porous nickel foam for overall water splitting. J. Alloys Compd. 2023, 947, 169474. [Google Scholar] [CrossRef]
  26. Meena, A.; Ha, M.; Chandrasekaran, S.S.; Sultan, S.; Thangavel, P.; Harzandi, A.M.; Singh, B.; Tiwari, J.N.; Kim, K.S. Pt-like hydrogen evolution on a V2O5/Ni(OH)2 electrocatalyst. J. Mater. Chem. A 2019, 7, 15794. [Google Scholar] [CrossRef]
  27. Wang, Y.Q.; Zhao, L.; Sui, X.L.; Gu, D.M.; Wang, Z.B. Hierarchical CoP3/NiMoO4 heterostructures on Ni foam as an efficient bifunctional electrocatalyst for overall water splitting. Ceram. Int. 2019, 45, 17128–17136. [Google Scholar] [CrossRef]
  28. Wang, X.F.; Liu, Z.H.; Guo, Z.G.; Ge, L.; Liu, Z.F. NiO-CoFe2O4 electrocatalyst prepared on Ni foam by one-step hydrothermal method for efficient overall water splitting. J. Mater. Sci. 2021, 56, 8575–8587. [Google Scholar] [CrossRef]
  29. Xiong, S.S.; Weng, S.T.; Tang, Y.; Qian, L.; Xu, Y.Q.; Li, X.F.; Lin, H.J.; Xu, Y.C.; Jiao, Y.; Chen, J.R. Mo-doped Co3O4 ultrathin nanosheet arrays anchored on nickel foam as a bi-functional electrode for supercapacitor and overall water splitting. J. Colloid Interface Sci. 2021, 602, 355–366. [Google Scholar] [CrossRef]
  30. Yang, J.H.; Xu, X.H.; Chen, M.M.; Yang, D.; Lu, H.; Sun, Y.Z.; Shao, C.; Song, Q.Q.; Zhang, J.; Gao, L.; et al. Morphology-controllable nanocrystal β-Ni(OH)2/NF designed by hydrothermal etching method as high-efficiency electrocatalyst for overall water splitting. J. Electroanal. Chem. 2021, 882, 115035. [Google Scholar] [CrossRef]
  31. Shen, Y.J.; Chen, Y.; Fang, S.K.; Park, J.K.; Feng, K. Plasma-modified Ni foam-supported CuCo2S4 nanowires as bifunctional electrocatalysts for high-performance overall water splitting. Arab. J. Chem. 2023, 16, 104989. [Google Scholar] [CrossRef]
  32. Nagajyothi, P.C.; Pavani, K.; Ramaraghavulu, R.; Shim, J. Microwave synthesis of NiMn2O4/Ni-foam: Efficient bifunctional electrocatalysts for overall water splitting. Int. J. Hydrogen Energ. 2023, 54, 691–699. [Google Scholar] [CrossRef]
  33. Tian, J.Q.; Cheng, N.Y.; Liu, Q.; Sun, X.P.; He, Y.Q.; Asiri, A.M. Self-supported NiMo hollow nanorod array: An efficient 3D bifunctional catalytic electrode for overall water splitting. J. Mater. Chem. A 2015, 3, 20056–20059. [Google Scholar] [CrossRef]
  34. Jin, Y.S.; Wang, H.T.; Li, J.J.; Yue, X.; Han, Y.J.; Shen, P.K.; Cui, Y. Porous MoO2 nanosheets as non-noble bifunctional electrocatalysts for overall water splitting. Adv. Mater. 2016, 28, 3785–3790. [Google Scholar] [CrossRef] [PubMed]
  35. Al-Naggar, A.H.; Shinde, N.M.; Kim, J.S.; Mane, R.S. Water splitting performance of metal and non-metal-doped transition metal oxide electrocatalysts. Coordin. Chem. Rev. 2023, 474, 214864. [Google Scholar] [CrossRef]
  36. Ding, X.D.; Yu, J.; Huang, W.Q.; Chen, D.Y.; Lin, W.; Xie, Z.L. Modulation of the interfacial charge density on Fe2P-CoP by coupling CeO2 for accelerating alkaline electrocatalytic hydrogen evolution reaction and overall water splitting. Chem. Eng. J. 2023, 451, 138550. [Google Scholar] [CrossRef]
  37. Niyitanga, T.; Kim, H. Synergistic effect of trimetallic-based CuxMox/Co1−xO nanoparticles/reduced graphene oxide as high efficient electrocatalyst for oxygen evolution reaction. Chem. Eng. J. 2023, 451, 138649. [Google Scholar] [CrossRef]
  38. GB/T 28124-2011; Determination of Trace Hydrogen, Oxygen, Methane and Carbon Monoxide in the Inert Gases—Gas Chromatographic Method. National Standardization Technical Committee of the Gas: Beijing, China, 2011.
  39. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  40. Kresse, G.; Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  41. Ma, X.; Zhang, W.; Deng, Y.; Zhong, C.; Hu, W.B.; Han, X.P. Phase and composition controlled synthesis of cobalt sulfide hollow nanospheres for electrocatalytic water splitting. Nanoscale 2018, 10, 4816–4824. [Google Scholar] [CrossRef]
  42. Song, F.Z.; Li, W.; Yang, J.Q.; Han, G.Q.; Liao, P.L.; Sun, Y.J. Interfacing nickel nitride and nickel boosts both electrocatalytic hydrogen evolution and oxidation reactions. Nat. Commun. 2018, 9, 4531. [Google Scholar] [CrossRef] [PubMed]
  43. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef] [PubMed]
Scheme 1. Schematic diagram for the fabrication procedures of VMoCoCOx@NF.
Scheme 1. Schematic diagram for the fabrication procedures of VMoCoCOx@NF.
Molecules 29 03591 sch001
Figure 1. (A) XRD patterns of VCoCOx@NF (lake blue), VMoCoCOx@NF (navy blue), and references of pure MoO3, CoCO3, and V2O5. Notably, the diffraction lines of Ni were removed during the XRD operation. (B) TEM and (C) HRTEM images of the NF substrate and VMoCoCOx@NF, respectively. EDS elemental mapping images of V (D), Co (E), and Mo (F) for the surface of VMoCoCOx@NF. FESEM images of VMoCoCOx@NF scanned at 500 μm (G), 50 μm (H), and 10 μm scales (I).
Figure 1. (A) XRD patterns of VCoCOx@NF (lake blue), VMoCoCOx@NF (navy blue), and references of pure MoO3, CoCO3, and V2O5. Notably, the diffraction lines of Ni were removed during the XRD operation. (B) TEM and (C) HRTEM images of the NF substrate and VMoCoCOx@NF, respectively. EDS elemental mapping images of V (D), Co (E), and Mo (F) for the surface of VMoCoCOx@NF. FESEM images of VMoCoCOx@NF scanned at 500 μm (G), 50 μm (H), and 10 μm scales (I).
Molecules 29 03591 g001
Figure 2. (A) XPS survey for VMoCoCOx@NF. High-resolution XPS spectra of (B) V 2p, (C) Co 2p, and (D) Mo 3d for VMoCoCOx@NF.
Figure 2. (A) XPS survey for VMoCoCOx@NF. High-resolution XPS spectra of (B) V 2p, (C) Co 2p, and (D) Mo 3d for VMoCoCOx@NF.
Molecules 29 03591 g002
Figure 3. Electrocatalytic properties of the catalysts selected for HER. (A) LSV polarization curves. (B) Comparison of Tafel slopes and overpotentials of the catalysts at a constant j = −10 mA cm−2. (C) Estimation of Cdl by plotting the ratios of catalyst ΔJ to scan rate. (D) Durability test of VMoCoCOx@NF performed at a constant j = −10 mA cm−2 within 100 h. The inset showed the LSV curves initially and after 10,000 CV cycles.
Figure 3. Electrocatalytic properties of the catalysts selected for HER. (A) LSV polarization curves. (B) Comparison of Tafel slopes and overpotentials of the catalysts at a constant j = −10 mA cm−2. (C) Estimation of Cdl by plotting the ratios of catalyst ΔJ to scan rate. (D) Durability test of VMoCoCOx@NF performed at a constant j = −10 mA cm−2 within 100 h. The inset showed the LSV curves initially and after 10,000 CV cycles.
Molecules 29 03591 g003
Figure 4. Electrocatalytic properties of the catalysts selected for OER. (A) OER LSV polarization curves. (B) Comparison of Tafel slopes and overpotentials of the catalysts at a constant j = 10 mA cm−2. (C) EIS curves of the four catalysts. (D) Chronoamperometry stability test of VMoCoCOx@NF performed at a constant j = 10 mA cm−2. The inset showed the LSV curves initially and after 10,000 CV cycles.
Figure 4. Electrocatalytic properties of the catalysts selected for OER. (A) OER LSV polarization curves. (B) Comparison of Tafel slopes and overpotentials of the catalysts at a constant j = 10 mA cm−2. (C) EIS curves of the four catalysts. (D) Chronoamperometry stability test of VMoCoCOx@NF performed at a constant j = 10 mA cm−2. The inset showed the LSV curves initially and after 10,000 CV cycles.
Molecules 29 03591 g004
Figure 5. (A) The overall water splitting activity of VMoCoCOx@NF in the alkaline electrolyzer. The experimentally measured and theoretically calculated millimolars of H2 and O2 as a function of time were plotted at 5 min intervals during the water splitting at a constant j = 10 mA cm−2. (B) The comparison of potentials between the VMoCoCOx@NF electrolyzer and other reported NF-supported transition metal-based catalysts [17,22,26,27,28,29,30,31,32].
Figure 5. (A) The overall water splitting activity of VMoCoCOx@NF in the alkaline electrolyzer. The experimentally measured and theoretically calculated millimolars of H2 and O2 as a function of time were plotted at 5 min intervals during the water splitting at a constant j = 10 mA cm−2. (B) The comparison of potentials between the VMoCoCOx@NF electrolyzer and other reported NF-supported transition metal-based catalysts [17,22,26,27,28,29,30,31,32].
Molecules 29 03591 g005
Figure 6. (A) The polarization curves of VMoCoCOx@NF||VMoCoCOx@NF, Pt/C@NF||RuO2@NF and NF||NF-coupled catalysts in a two-electrode configuration. (B) Comparison of the overpotentials between VMoCoCOx@NF||VMoCoCOx@NF and RuO2@NF||Pt/C@NF at j = 10 mA cm−2 and 100 mA cm−2. (C) The polarization curves before and after the 120 h stability test. Inset (C) showed the VMoCoCOx@NF||VMoCoCOx@NF alkaline water electrolyzer. (D) Chronopotentiometry stability tests performed at constant j = 10 mA cm−2 and 100 mA cm−2 for 120 h.
Figure 6. (A) The polarization curves of VMoCoCOx@NF||VMoCoCOx@NF, Pt/C@NF||RuO2@NF and NF||NF-coupled catalysts in a two-electrode configuration. (B) Comparison of the overpotentials between VMoCoCOx@NF||VMoCoCOx@NF and RuO2@NF||Pt/C@NF at j = 10 mA cm−2 and 100 mA cm−2. (C) The polarization curves before and after the 120 h stability test. Inset (C) showed the VMoCoCOx@NF||VMoCoCOx@NF alkaline water electrolyzer. (D) Chronopotentiometry stability tests performed at constant j = 10 mA cm−2 and 100 mA cm−2 for 120 h.
Molecules 29 03591 g006
Figure 7. DFT calculation. (A) DFT-optimized structures of OH*, H*, and H2O* intermediates on the Co and Mo sites of VCMoCoCOx for HER. (B) DFT-optimized structures of *, OH*, O*, and OOH* intermediates on the Co and Mo sites of VCMoCoCOx for OER. Color code: The Co atom is green, the Mo atom is purple, the V atom is yellow, the O atom is red, and the H atom is black. (C) ΔGH* for the HER mechanism on the surface of CoCO3@NF, VMoCoCOx@NF, VMoCOx@NF, and VCoCOx@NF. (D) ΔG for the OER mechanism on the surfaces of VMoCOx@NF, VMoCoCOx@NF, and VCoCOx@NF. The dotted lines in (C,D) represent the ΔG jumps at each step.
Figure 7. DFT calculation. (A) DFT-optimized structures of OH*, H*, and H2O* intermediates on the Co and Mo sites of VCMoCoCOx for HER. (B) DFT-optimized structures of *, OH*, O*, and OOH* intermediates on the Co and Mo sites of VCMoCoCOx for OER. Color code: The Co atom is green, the Mo atom is purple, the V atom is yellow, the O atom is red, and the H atom is black. (C) ΔGH* for the HER mechanism on the surface of CoCO3@NF, VMoCoCOx@NF, VMoCOx@NF, and VCoCOx@NF. (D) ΔG for the OER mechanism on the surfaces of VMoCOx@NF, VMoCoCOx@NF, and VCoCOx@NF. The dotted lines in (C,D) represent the ΔG jumps at each step.
Molecules 29 03591 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, W.; Xu, L.; Ye, R.; Yang, P.; Zhu, J.; Jiang, L.; Wu, X. Molybdenum and Vanadium-Codoped Cobalt Carbonate Nanosheets Deposited on Nickel Foam as a High-Efficient Bifunctional Catalyst for Overall Alkaline Water Splitting. Molecules 2024, 29, 3591. https://doi.org/10.3390/molecules29153591

AMA Style

Wang W, Xu L, Ye R, Yang P, Zhu J, Jiang L, Wu X. Molybdenum and Vanadium-Codoped Cobalt Carbonate Nanosheets Deposited on Nickel Foam as a High-Efficient Bifunctional Catalyst for Overall Alkaline Water Splitting. Molecules. 2024; 29(15):3591. https://doi.org/10.3390/molecules29153591

Chicago/Turabian Style

Wang, Wenxin, Lulu Xu, Ruilong Ye, Peng Yang, Junjie Zhu, Liping Jiang, and Xingcai Wu. 2024. "Molybdenum and Vanadium-Codoped Cobalt Carbonate Nanosheets Deposited on Nickel Foam as a High-Efficient Bifunctional Catalyst for Overall Alkaline Water Splitting" Molecules 29, no. 15: 3591. https://doi.org/10.3390/molecules29153591

APA Style

Wang, W., Xu, L., Ye, R., Yang, P., Zhu, J., Jiang, L., & Wu, X. (2024). Molybdenum and Vanadium-Codoped Cobalt Carbonate Nanosheets Deposited on Nickel Foam as a High-Efficient Bifunctional Catalyst for Overall Alkaline Water Splitting. Molecules, 29(15), 3591. https://doi.org/10.3390/molecules29153591

Article Metrics

Back to TopTop