Conversion of Catalytically Inert 2D Bismuth Oxide Nanosheets for Effective Electrochemical Hydrogen Evolution Reaction Catalysis via Oxygen Vacancy Concentration Modulation

Ziyang Wu, Ting Liao,2✉, Sen Wang, Janith Adikaram Mudiyanselage,Aaron S. Micallef,5, Wei Li, Anthony P. O’Mullane2,, Jianping Yang, Wei Luo,Kostya Ostrikov2,, Yuantong Gu,2, Ziqi Sun2,✉

ABSTRACT Oxygen vacancies(Vo) in electrocatalysts are closely correlated with the hydrogen evolution reaction (HER) activity. The role of vacancy defects and the effect of their concentration, however, yet remains unclear. Herein,Bi2O3, an unfavorable electrocatalyst for the HER due to a less than ideal hydrogen adsorption Gibbs free energy (ΔGH*), is utilized as a perfect model to explore the function of Vo on HER performance.Through a facile plasma irradiation strategy, Bi2O3 nanosheets with different Vo concentrations are fabricated to evaluate the influence of defects on the HER process.Unexpectedly, while the generated oxygen vacancies contribute to the enhanced HER performance, higher Vo concentrations beyond a saturation value result in a significant drop in HER activity. By tunning the Vo concentration in the Bi2O3 nanosheets via adjusting the treatment time, the Bi2O3 catalyst with an optimized oxygen vacancy concentration and detectable charge carrier concentration of 1.52 × 1024 cm−3 demonstrates enhanced HER performance with an overpotential of 174.2 mV to reach 10 mA cm−2, a Tafel slope of 80 mV dec−1, and an exchange current density of 316 mA cm−2 in an alkaline solution, which approaches the top-tier activity among Bi-based HER electrocatalysts. Density-functional theory calculations confirm the preferred adsorption of H* onto Bi2O3 as a function of oxygen chemical potential (∆μO) and oxygen partial potential (PO2) and reveal that high Vo concentrations result in excessive stability of adsorbed hydrogen and hence the inferior HER activity. This study reveals the oxygen vacancy concentration-HER catalytic activity relationship and provides insights into activating catalytically inert materials into highly efficient electrocatalysts.

KEYWORDS Alkaline hydrogen evolution reaction; Bismuth oxide; Plasma irradiation; 2D materials; Oxygen vacancy

1 Introduction

Green hydrogen has been identified as a promising energy source for sustainable development and an ideal candidate to replace traditional fossil fuels to address global environmental concerns [1]. Among various hydrogen production technologies, electrochemical water splitting encompassing the electrochemical hydrogen evolution reaction (HER) has been demonstrated to be a simple yet efficient route for the conversion of water into H2[2—4]. Compared with HER in an acidic electrolyte, where the intermediate hydrogen is generated from proton reduction (H++ e−→ H*), the alkalinebased reaction needs to overcome a higher energy barrier for the additional water dissociation step (H2O + e− → H* + OH−)[5—7]. However, hydrogen production in alkaline media is desired for industrial production, due to the fact that often milder reaction conditions and higher anode catalyst activity are achievable for the catalysts towards the oxygen evolution reaction (OER) in alkaline electrolytes [8]. Noble metal electrocatalysts have so far attracted great attention for their high catalytic activity, but scarcity and their precious nature inhibit their large-scale application in electrocatalysis [9—12]. As a result, rational design and development of low-cost and highperformance electrocatalysts, e.g., transition metal oxide or carbon-based catalysts, for alkaline electrocatalysis, will play a vital role in the hydrogen economy [13—19].

Apart from exploring efficient earth-abundant metal oxidebased catalysts as potential Pt-free catalysts [20—22], various strategies have also been developed to enhance the activity of some existing electrocatalysts for alkaline conditions, such as multiple active sites tailoring [23, 24], electronic structure engineering [25—27], wettability control [28—30], etc., which can either optimize the interactions between the catalyst’s surface and the reaction intermediates or reduce the reaction energy barriers [31]. It should be noted that the oxygen or oxygen defective sites of metal oxide catalysts are strongly associated with electrocatalytic performance, especially in the water dissociation step of the HER. It has been well studied that the generation of oxygen vacancies (Vo) can not only improve electronic conductivity but also change the interactions between the metal 3dand the O-2pband of the catalysts, which plays a vital role in the lattice oxygen involved in the OER [32]. Oxygen vacancies have also been reported as anchoring sites for constructing active single-atoms. In a recent work, Pt single atoms were anchored onto theVosites of MoO2to achieve superior HER activity, where theVonot only stabilizes the single atoms but also modulates the electronic structure of the catalysts [33]. Thus, creatingVoin metal oxides is an efficient strategy for modulating the surface electronic structure and tuning the intrinsic catalytic activity of electrocatalysts [34, 35]. Typical methods for creatingVoin the catalysts include thermal treatment, reduction reaction, doping, high-energy particle irradiation, etc.[36]. Although oxygen defects have been widely reported in enhancing the electrocatalysis performance, the role of oxygen vacancies in functionalizing the catalysts and the underlying mechanism are far from being fully understood.Thus, it is essential to create an effective model to accurately identify the role ofVofor the HER under alkaline conditions.

Bismuth (Bi), aP-block post-transition metal and a remarkably harmless element among the toxic heavy metals, has sparked interest in areas varying from medical to industrial chemistry [37]. Specifically, it has been demonstrated to be an outstanding electrocatalyst for carbon dioxide [38, 39] and nitrogen reduction [40—42], due to its inert HER activity and unfavorable free energy of hydrogen adsorption (∆GH*) which inhibits this competing reaction. Nørskov et al. have reported that the Bi catalyst has ultrahigh binding energy towardH*intermediates among various metal surfaces [43]. In addition,Bi2O3is easily reduced into the metallic state and can accept electrons or loose oxygen ions to form a metal-like conductive surface [44, 45]. Compared with some other transition metals which usually need a high temperature reduction process or the use of a reducing agent [46], the property of Bi2O3enables the formation of reductive defects under mild conditions. For this reason, defects can be conveniently introduced into bismuth oxide via facile methods [47]. Therefore, the inherently sluggish HER activity and the simplicity of constructingVoof 2D Bi2O3nanosheets endow this material as an ideal candidate to reveal the relationship between the concentration ofVoand HER performance.

Herein, based on 2D Bi2O3nanosheets, we proposed a low-temperature plasma enabled approach to tune the oxygen vacancies concentration in the catalysts by plasma irradiation and investigated the corresponding alkaline HER performance to elucidate how the oxygen vacancies improve the alkaline activity of the catalyst, and tried to provide a principle to tune the inactive materials into active catalysts while maintaining their inherent chemical durability. This work discovered that the oxygen vacancy concentration is closely related with the alkaline HER catalytic activity of Bi2O3, but there is a saturation value ofVoto achieve the highest activity, over which the activity starts to drop. We then studied the relationship of the adsorption of the hydrogen intermediate onto the various active sites as a variation of partial oxygen pressure and the oxygen chemical potential by applying density functional theory (DFT) calculations. The theoretical investigation reveals that the single H* on O site (H*O)is the preferred adsorption model at a lowVoconcentration whereas the adsorption of single H* on twin oxygen vacancies site (2Vo− H*O) is dominant at a highVoconcentration. However, theVoconcentration cannot be too high, otherwise the desorption energy of H* will be too high to desorb for the following step. As a result, the 2D Bi2O3nanosheets with a mediumVoconcentration and a detectable charge carrier concentration of 1.52 × 1024cm−3presented a superior HER performance with a small overpotential (174.2 mV at 10 mA cm−2) and a low Tafel slope(80 mV dec−1) but a high exchange current density of 320 mA cm−2in 1 M KOH achieved in this HER inert catalyst.

2 Experimental Section

2.1 Materials

The chemicals used for the synthesis of 2D Bi2O3nanosheets include bismuth nitrate pentahydrate(Bi(NO3)3·5H2O), polyethylene oxide-polypropylene oxide-polyethylene oxide (PEO20-PPO70-PEO20, Pluronic P123), ethylene glycol (C2H6O2, EG), Nafion solution(117 solution, Aldrich), and deionized (DI) water.

2.2 Synthesis of Bi2O2CO3@Ni Foam (RT), Bi2O3@Ni Foam (Pl‑0)

First, P123 (0.3 g) was dissolved in a mixed solution of ethanol(6 mL) and water (4.25 mL) (Solution A). Then, 13 mL EG was used to dissolve Bi(NO3)3·5H2O (0.2 g) to form homogeneous Solution B, which was mixed with solution A after 5 min ultrasonication. The mixed solution was kept at room temperature with an aging time of 24 h. Ni foam substrates were subsequently cleaned with 1 M HCl solution, acetone,ethanol, and DI water in an ultrasonic bath for 10 min. The cleaned Ni foam (1 × 1 cm2) was placed into the autoclave with the precursor solution, and the time of hydrothermal reaction maintained at 170 °C was set as 3 h for the uniform growth of 2D nanosheets. The sample was subsequently collected after complete cooling of the reactor, washed with ethanol twice and dried in an oven at 80 °C for 24 h (the as-synthesized sample was denoted as RT). The dried powders were then heat-treated in a furnace at 400 °C for 2 h in air with a ramping rate of 5 °C min−1to get the bismuth oxide sample on Ni foam (Pl-0).

2.3 Plasma Irradiation Treatment

The plasma treatment was conducted in an atmosphericpressure quartz dielectric barrier discharge plasma reactor,where the quartz reaction chamber was sealed by a quartz lid cover. A CTP 2000 K power supply was employed to provide an amplitude of 26 kV and frequency of 9.1 kHz. A Rigol DS6104 oscilloscope was used to record the voltage and current changes. The feeding gas used in this study was N2.The Bi2O3covered Ni foam samples at a size around 1 × 1 cm2were placed at the center of the reactor and exposed with the desired time of 15, 30, and 60 min at a set power of 300 W. The samples after treatment were denoted as Pl-15, Pl-30, and Pl-60. Catalyst loadings after plasma treatment were weighted, which were about 4.8, 4.4, 4.4, 4.1,and 3.9 mg cm−2, respectively, for the samples of RT, Pl-0,Pl-15, Pl-30, and Pl-60. Andor SR-500 spectrometer was utilized to record the optical emission spectrum (OES) of the plasma. The Newton CCD camera with fiber optic cable were equipped on the device for recording and safety issue.

2.4 Materials Characterizations

A field emission scanning electron microscope (FESEM, JSM-7001F, JEOL) was used to characterize the surface structure.A transmission electron microscope (TEM, JEOL 2100) was employed to characterize the morphology and crystal information. The sample surface element valence and chemistry property were measured by X-ray photoelectron spectroscopy(XPS, Kratos AXIS Supra photoelectron with Al-Kαradiation (hν = 1486.6 eV)). The calibration was performed by the containment carbon with a value of 284.8 eV. Grazing incidence diffraction (GID) geometry-based X-ray diffraction(XRD) was employed for the identifying of Ni foam loaded Bi2O3nanosheets, and the spectrum was collected by a Rigaku Smartlab diffractometer (Cu radiation, 40 kV and 40 mA). The XRD patterns were collected at a scanning speed of 2° 2θ min−1with the range from 5 to 90° 2θ maintained at a 0.02° step size.The phases of the samples were identified by JADE (V2010,Materials Data Inc.) and EVA (V5, Bruker) with PDF4 + 2021 reference database. Rietveld refinement was performed with a TOPAS package (V6, Bruker) for the quantitative analysis.Renishaw Raman microscope was used for the collecting of Raman spectra. The sample surface area was revealed by a Micromeritics Tristar II 3020 Surface Area and Porosity Analyzer. Brunauer—Emmett—Teller method was employed for the related surface area calculation. NT-MDT Solver Pro atomic force microscope (AFM) was employed to characterize the Bi2O3nanosheets thickness before and after plasma treatment.Samples were prepared for electron paramagnetic resonance(EPR) spectroscopy by scratching Bi2O3powder from the surface of the paramagnetic Ni foams. EPR spectrometer (Magnettech MiniScope MS400) with a benchtop setting was used to collect the EPR spectra. The X-band (ν) was 9.4 GHz, and the measurements were conducted at room temperature. For better comparation of the EPR results, the samples were tested with same weight (13 mg), filling almost 1.5 cm-height at the bottom of the sample tube (diameter = 3 mm).

2.5 Electrochemical Measurements

A CHI760E workstation (CH Instruments, Inc.) was used to perform the electrochemical tests in a three-electrodesetup cell with a graphite counter electrode, a saturated calomel electrode (SCE) reference electrode, a 1 × 1 cm2Bi2O3loaded Ni foam working electrode. Conversion equation(ERHE=ESCE+ 0.244 + 0.059 × pH) was used to calculate the potentials relative to reversible hydrogen electrode (RHE).Pt/C (20 wt%) at a mass load of 4 mg cm−2was used as the benchmark to evaluate the HER activity. 1.0 M KOH solution was used as electrolyte, and the scan rate of the polarization measurement was 5 mV s−1. iR-correction (95%) was applied during electrocatalytic measurements. Electrochemical impedance spectroscopy (EIS) measurements were performed over 0.1 to 100,000 Hz at an overpotential of 10 mV.The Mott-Schottky measurements were performed between− 1 and 0.5 V (vs. SCE) at a frequency of 1000 Hz. Cyclic voltammetry curves (CVs) were collected within a non-Faradaic potential range under the open-circuit potential of the system for estimating the electrical double-layer capacitance(Cdl) and the electrochemical surface area (ECSA).

2.6 Theoretical Calculations

A Quantum-Espresso package with a spin-polarized DFT framework was employed for the theoretical calculations proposed in this work [48]. Ultrasoft pseudopotentials was used to describe the electron—ion interactions, and the generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) functional was applied for the exchangecorrelation interactions [49, 50]. Plane waves (PWs) basis sets were employed for the calculations of the Kohn—Sham(K—S) orbitals and the charge density. Bi2O3(010) supercells in either defect free state or with one or two oxygen vacancies involved were built to study their hydrogen adsorption stability, where a vacuum space of 15 Å along the c direction of the supercell was designed to avoid the mirror interactions. The Brillouin zone with a Γ and a 5 × 5 × 1k-point mesh was used for geometric optimization and electronic structure analysis [51]. 1 × 10−7eV for the energy and 1 × 10−4eV Å−1for the force were set as the convergence criteria in structural optimization.

The stability of the hydrogen adsorption on either the defect free or the oxygen vacancies involving the Bi2O3surface was calculated by using atomistic thermodynamics methods [52,53]. The formation energies were defined as follows:

In whichρ0is the reference state pressure. The chemical potentialsμo(T,ρ0) at the desired temperature and the reference pressureρ0were obtained from the thermochemical reference tables [54].

3 Results and Discussion

3.1 Characterization of Plasma Irradiated Catalysts

3.1.1 SEM andTEM Characterizations

The fabrication of 2D bismuth oxide nanosheets on Ni foam was carried out by a hydrothermal molecular self-assembly method [55], followed by high temperature calcination for conversion into Bi2O3at 400 °C in air for 2 h (Fig. S1). Subsequently, plasma processing was conducted under N2plasma with different irradiation durations (15, 30, and 60 min) to achieve different oxygen vacancy concentrations (Fig. 1).Gas-phase plasma usually interacts with the substrate surface via reactive intermediates derived from the gas. To characterize the activated species during the N2plasma, an OES was recorded, and the typical band system at 300—400 nm (316,337, 357, and 380 nm) evidently proved the second positive system from N2during the plasma activation process (Fig. S2)[56]. Plasma has been reported to be an effective approach to create oxygen vacancies in metal oxides [57—59].

To reveal the surface morphology evolution during the plasma treatment, SEM characterization was conducted.The as-synthesized ultrathin 2D Bi2O3nanosheets on Ni foam demonstrated a uniform graphene-like structure (Fig.S3a, b), and no collapse occurred after calcination at 400 °C(Fig. 2a). The plasma irradiated Bi2O3nanosheets had no significant change in their 2D form but the formation of through-holes with increasing numbers and pore size was observed with prolonged treatment time (Figs. 2 and S3c, d).

Low-magnification TEM images confirmed the structural changes with the plasma treatment (Fig. 2b, e, and h), where the 2D morphology of the samples were maintained with the formation of nanoholes. High-resolution TEM (HRTEM)detected phase transformation with the plasma treatment. For the as-synthesized Bi2O3nanosheets without treatment (Pl-0,Fig. 2c), the lattice fringes with a spacing of 0.39 nm can be assigned to the (200) planes ofβ-Bi2O3[45]. After plasma processing of 15 min (Pl-15, Fig. 2f) and 30 min (Pl-30, Fig. 2i),the (002) plane of theβ-Bi2O3phase with a spacing of 0.35 nm was identified. The co-existence ofα-Bi2O3andβ-Bi2O3phases implies that high-power plasma irradiation contributes to the atomic rearrangement of Bi2O3[60]. The SAED patterns collected on the Pl-30 sample also confirmed the formation of theα-Bi2O3phase (Fig. 2j), where the (120) planes forα-Bi2O3were identified together with the lattice planes forβ-Bi2O3[61].

The thickness and the pore size of the 2D Bi2O3nanosheets with plasma treatment at different durations were examined by AFM. As shown in Fig. 2k, increasing the irradiation time from 0 to 60 min resulted in the formation of more and larger nanopores in the nanosheets and a slight increase of the thickness of the nanosheets from 1.2 to 1.5 nm (Fig. S4a). The dark-field TEM observation with corresponding EDS mapping on Pl-30 indicate that the chemical distribution has not been disturbed by the plasma irradiation and uniform distributions of Bi, O, and C elements were clearly identified (Fig. S4b).

3.1.2 XRD andRaman Characterizations

To further evaluate the phase transformation during the plasma irradiation process, XRD data were collected for qualitative and semi-quantitative analysis of the phases. As shown in Fig. 3a, in the Pl-0 sample, theβ-Bi2O3(PDF#04-008-7003, P-421c,a= 0.7755 nm,c= 0.5659 nm) phase dominated in the Bi2O3nanosheets co-existing with traces of NiBi (PDF#04-007-2591) and the Ni metal substrate(PDF#04-004-6807). After 15 min of plasma treatment,only a very small amount of the α-Bi2O3phase was detected(PDF#04-017-2112, P21/c,a= 0.5854 nm,b= 0.8165 nm,c= 0.7508 nm,β= 112.82°) (Fig. S5). However, when the treatment time increased to 30 min and 60 min, the ratios ofα-Bi2O3toβ-Bi2O3were 1.00 and 1.32 (Fig. S5), respectively, which could be the result of the continuous escape and rearrangement of O atoms during the plasma treatment[62]. It is interesting to note that, with a substantial increase of plasma etching time, an obvious shift of the (201) plane ofβ-Bi2O3and the appearance of new peaks for the (002),(111), and (120) planes ofα-Bi2O3were identified, which were the result of phase transformation and lattice distortion(Fig. 3b). This result implies the expansion of the lattice along specific planes, due to oxygen vacancies generation activated by the plasma process [63]. Based on the TEM characterizations, the nanosheets presented a polycrystalline feature. The average crystalline size obtained from the XRD patterns by Rietveld refinement are of 41.2 and 25.7 nm forα-Bi2O3andβ-Bi2O3, respectively, which are coincided with the TEM observations.

Raman characterization was also conducted to identify the surface changes introduced by the plasma (Fig. 3c). For sample Pl-0, the Bi—O stretching modes at 90 and 120 cm−1fromβ-Bi2O3were clearly observed. It should be noted that the vibrational bands at 90, 120, and 316 cm−1are shared by bothβ-Bi2O3andα-Bi2O3, while the newly formed peaks at 182 and 210 cm−1were as a result of the phase transformation into the monoclinic phaseα-Bi2O3[64]. Compared with Pl-0, the plasma processed samples demonstrated a red shift and peak broadening, as shown in the inset in Fig. 3c, which are usually associated with structural softening induced by heating, defects, etc., in the crystal structure [65, 66]. Therefore, both XRD and Raman characterizations confirmed the formation of oxygen vacancies in the Bi2O3lattice. Owing to the creation of a large number of pores in the nanosheets,the specific surface area increased stepwise from 25.8 to 64.5 m2g−1upon prolonging the treatment time to 60 min(Fig. 3d).

3.1.3 XPS, EPR andMott‑Schottky Measurements onOxygen Vacancy Concentrations

The chemical compositions and the defects of the 2D Bi2O3nanosheets after plasma processing were identified by XPS measurements. The element survey spectrum indicated the presence of Ni, Bi, and O on the surface (Fig. S6). Even though N2-plasma has been used in this treatment, no N-doping was detected in all treated samples. In the Bi 4fspectrum (Fig. 4a), two peaks at 158.7 and 164 eV (∆E= 5.3 eV)assigned to the spin-orbit doublet of Bi 4f7/2and Bi 4f5/2for Bi3+were observed on Pl-0. Upon plasma treatment, these two peaks shifted to lower binding energy. For the Pl-60, the peaks were located at 158.2 and 163.5 eV, corresponding to a 0.5 eV shift. This type of down shifting is due to lowering the coordination of Bi atoms by the formation of plasma induced oxygen vacancies [58, 67].

Fig. 1 Schematic illustration of plasma irradiation on the formation of oxygen vacancies in 2D Bi2O3 nanosheets for electrocatalytic HER

Fig. 2 Morphology evolution of 2D Bi2O3 nanosheets upon plasma irradiations. SEM images (a, d and g), low-magnification TEM images (b, e and h) and high-resolution TEM images (c, f and i) of Bi2O3 before irradiation (Pl-0), irradiated for 15 min (Pl-15), and 30 min (Pl-30), respectively. SAED patterns (j), AFM images on the changes of thickness and pore sizes with different plasma irradiation durations (k)

The O 1s core levels can provide more information on the formation of oxygen vacancies (Fig. 4b), where the O 1s state located at around 530.8 eV is associated withVoor more accurately the change of chemical states of the lattice oxygen resulted by the formation of oxygen vacancies, and the one at 529.6 eV corresponds to lattice oxygen [62, 68]. Although the assignment ofVohas been disturbed by attached the hydroxyls, water molecules, and organic contaminants [69], the same fabrication, storage,and treatment conditions of the samples ensure that the increase of the high-energy core level of O 1s after plasma irradiation were ascribed by surface oxygen vacancies formation [58, 70, 71]. Therefore, we did not deconvolute the oxygen peak at 530—533 eV but exclude the contributions originated from the hydroxyls and the adsorbed water molecules by taking the non-treated sample Pl-0 as the baseline to calculate the generated oxygen vacancy concentrations resulted by the plasma irradiations. It is clear that the intensity and the proportion of the deconvolution associated withVoincreased significantly with the prolonged plasma treatment, corresponding to the increase of oxygen vacancy concentrations to 29.3% after exposing under plasma up to 60 min.

X-band EPR spectroscopy was also utilized to assess the formation of oxygen vacancies (Fig. 4c) [72]. The characteristic signal observed at g ≈ 2.001 is consistent with the formation ofVodefects upon plasma etching. The observed increase of the EPR signal intensity with plasma treatment time suggested that the number of oxygen vacancies is proportional to the treatment time. Table S2 displays the detailed values of the EPR variations.

The Mott-Schottky (M-S) test can provide quantitative information about the effective charge carriers generated within the semiconductors, and the increase of the carriers generated by the plasma irradiation could be considered as the generation of oxygen vacancies. Figure 4d presents the carrier density (N) as a function of the applied bias voltages recorded in the M-S tests, andNcan be evaluated from the slope of M-S curves based on Eq. (4) [62].

Fig. 3 Phase evolution of 2D Bi2O3 nanosheets upon plasma irradiations. XRD spectrum (a), locally enlarged image of XRD patterns (b),Raman spectrum with the inset showing the peak shifting and broadening (c), and N2 adsorption and desorption spectrum (d) of samples Pl-0,Pl-15, Pl-30, and Pl-60

Here,εis the dielectric constant of the material, which is 18.2 for Bi2O3[73, 74];ε0is the vacuum permittivity(8.85 × 10−12F m−1);eis the electron charge (1.6 × 10−19C); andVis the applied bias. Figure 4d shows the plot of 1/C2vsVof different samples measured in 1 M KOH electrolyte. The negative slopes indicate thep-type semiconductor behavior of Bi2O3[75]. The calculated charge carrier densityNfor the 30 min treated Bi2O3(Pl-30) is 1.52 × 1024cm−3, which is nearly 1 order of magnitude higher than that of Pl-0 (2.53 × 1023cm−3), confirming the significant increase ofVoupon prolonging the plasma treatment. A further increase in irradiation time to 60 min, resulted in the carrier density decreasing slightly to 1.34 × 1024cm−3, although a higherVoconcentration was indicated in both the XPS and EPR characterizations in Pl-60, which may be a result from over irradiation induced phase transformation and structure deformation,which may hinder electron transfer as observed previously in Xiao’s work [76]. Figure 4e summarizes the relationship between the plasma inducedVoconcentration and the irradiation time evaluated by these different approaches.Clearly, theVoconcentration calculated from XPS and EPR confirmed a nearly linear relationship with plasma irradiation time, demonstrating the effectiveness of the plasma treatment on generatingVoin Bi2O3nanosheets.

3.2 Electrocatalysis Performance of the Catalysts

The effect of tailoredVoon the electrocatalytic performance of the HER-inert Bi2O3was evaluated in a 1 M KOH electrolyte with a standard three electrodes configuration. Figure 5a exhibits the linear sweep voltammetry(LSV) curves of the samples with the different plasma irradiation time, namely, with the differentVocontents. It is interesting that the Bi2O3treated for 15—30 min presented significantly improved HER performance, and particularly, the Pl-30 sample showed the lowest onset potential after iR-compensation, demonstrating its lowest energy barrier to trigger the HER. To reach a current density of 10 mA cm−2, the requested overpotentials were 286.1,262.1, 174.2, and 181.2 mV, respectively, for Pl-0, Pl-15,Pl-30, and Pl-60, respectively. The Pl-30 sample exhibited the lowest overpotential to reach 10 mA cm−2, which is more than a 110 mV improvement compared with the as-synthesized catalyst. While the HER activity of Bi2O3after plasma treatment is still lower than that of a commercial Pt/C electrocatalyst (Fig. S7), the tuned Bi2O3material exhibited promising potential for high current density catalysis, and only an overpotential of 310 mV is required to reach a current density of 300 mA cm−2. Most critically,this investigation reveals the success of the strategy of tuning the HER-inert catalysts into active electrocatalysts via oxygen vacancy modulation. Compared with previously very few reported Bi-based electrocatalysts for alkaline HER (Fig. S8), the plasma inducedVoactivated 2D Bi2O3presented the best performance.

For a better understanding on howVoconcentration contributes to the HER activity, we plotted the overpotential at 10 mA cm−2as a variation of theVoconcentration measured by using XPS, EPR, and M-S techniques (Fig. 5b). A near linear relationship between the overpotential and theVoin the Pl-0, Pl-25, and Pl-30 catalysts was identified for the data acquired by these methods, and Pl-30 presented the best performance. However, the Pl-60 which possess the highest absoluteVoidentified by XPS and EPR exhibited decreased activity toward the HER. From this result, we can understand that there is a saturation level for theVogenerated in these materials [77]. If we look at the charge carrier density measured by the M-S test, a very close linear relationship exists between the overpotential and the charge carrier density for all samples, which reveals that only the effective carrier that can move freely can contribute to the HER catalysis, and therefore, the M-S measurement could be a more reliable method to identify the available activeVosites for electrocatalysis. Therefore, Pl-60 possessed the highestVoconcentration, but someVoare trapped and combined to become unmovable “dead” defects which cannot contribute to the HER.

Fig. 4 Chemical states and oxygen vacancy evolution of 2D Bi2O3 nanosheets upon plasma irradiations. High resolution XPS spectra of Bi 4f(a) and O 1s (b) bands, electron paramagnetic resonance (EPR) spectra (c), the Mott-Schottky plots (d) of Pl-0, Pl-15, Pl-30, and Pl-60 samples.The variation of oxygen vacancy concentration of 2D Bi2O3 with plasma irradiation time calculated based on the XPS, EPR, and Mott-Schottky measurements (e)

The HER kinetics as a function ofVoconcentration can also be determined from the Tafel slopes and the exchange current density (Figs. 5c and S9). The Pl-30 catalyst presented the lowest Tafel slope of 80 mV dec−1, which is comparable to Pl-60 (90 mV dec−1) but much smaller than that of Pl-0 (107 mV dec−1) and Pl-15 (103 mV dec−1), and the highest exchange current density (316 mA cm−2), a descriptor of catalytic activity, among all the samples. The Tafel slope can be used to address the elementary steps and the rate determining steps during the catalytic reactions [78].For the HER electrocatalysis in alkaline condition involving both a water dissociation step and the adsorption and desorption process of hydrogen intermediates, water dissociation could be the rate determining step when its value is around 120 mV dec−1[78, 79]. Here, the relatively lower Tafel slope of 80 mV dec−1suggests that the water dissociation step became more favorable in theVoconcentration tuned electrocatalysts. As observed in the LSV measurements, a properVoconcentration is critical for alkaline HER. Pl-30 with the most suitableVoconcentration presented the most favorable HER kinetics and the highest catalytic activity to deliver the most optimal HER reactivity.

The electrochemically active surface area (ECSA) provides a connection between the coverage of active sites and the catalytic activity, which can be estimated from the double-layer capacitance (Cdl) through CVs measured within a non-Faradaic region at different scan rates. Figure 5d displays theCdlcalculated based on a series of CV measurements on the 2D Bi2O3catalysts (Fig. S10), in which a higherCdlmeans the higher available electrochemically active surfaces during HER due to the positive correlation between theCdland the ECSA. In this case, the Pl-30 showed the highestCdlvalue of 68.5 mF cm−2, which is superior to Pl-60 (60.2 mF cm−2), Pl-15 (51.2 mF cm−2), and PL-0 (47.7 mF cm−2). This result indicates that the plasma processing can increase the number of electrochemically active sites, but excessive irradiation decreases the effectiveness of the catalysts even with an increased physical surface area. In addition, the superhydrophilic surfaces of the samples (Fig. S11) associated with the formation of oxygen vacancies also contributed to more active reaction sites during the HER process.

As an important indicator for the electrode kinetics of the HER, EIS provides detailed information on charge transfer kinetics. Figure 5e shows the Nyquist plot with an equivalent circuit and the simulated ohmic resistance (Rs), charge transfer resistance (Rct) as well as the constant phase element (RCPE) to illustrate the reaction rates. It can be deduced from theRctvalues of Pl-0 (4.1 Ω), Pl-15 (3.9 Ω), Pl-30(2.9 Ω), and Pl-60 (3.0 Ω) that the 30 min plasma irradiated sample displayed the highest charge transfer kinetics as a result of the plasma inducedVo. The stability of the oxygen vacancy tuned electrocatalysts was also monitored by a chronoamperometry test (Fig. 5f). Continuous running of the HER test for 50 h of Pl-30 (− 1.3 V vs SCE) was performed,and the overpotential variations at 10 mA cm−2demonstrated the good durability of the catalyst, in which only a 2.3% increasing of the overpotential from 174 to 178 mV were recorded before and after the long-time durability test.Previous reports informed that it is difficult to maintainVovacancies promoted HER performance for long-term, as theVointends to be consumed during the HER process. Here,the plasma induced highVoconcentration was retained for achieving excellent catalysis stability. We also examined the morphology and surface composition variation of the catalysts after the stability test. All the microstructural and compositional characterizations, including XRD (Fig. S12),SEM (Fig. S13), and XPS (Fig. S14), indicated almost no changes occurred on Pl-30, except for partial reduction of Bi2O3into metallic bismuth on the surface.

3.3 DFT Calculations Study

Based on the above experimental results, it is interesting that a properVomodulation in the 2D Bi2O3catalysts is very critical to performance. DFT calculations were, therefore,performed to understand the mechanism of the enhanced electrocatalytic HER performance in Bi2O3with the variation of the oxygen vacancy concentration. The active hydrogen adsorption sites related with both Bi and O atoms were used to find the most stable adsorption sites by evaluating the surface formation free energies (Fig. 6).

The considered hydrogen intermediate adsorption types,as the representative models shown in Fig. 6a, include single H* on O site (H*O), single H* on single oxygen vacancy site (Vo− H*Vo), single H* on twin oxygen vacancies site(2Vo− H*Vo), single H* on Bi site (H*Bi), single H* on the bismuth site beside the single oxygen vacancy (Vo− H*Bi),single H* on the Bi site beside the twin oxygen vacancies site (2Vo− H*Bi), two H* on O site (2H*O), two H* on single oxygen vacancy site (Vo− 2H*Vo), two H* on twin oxygen vacancies site (2Vo− 2H*Vo), two H* on Bi site (2H*Bi), two H* on the bismuth site beside the single oxygen vacancy(Vo− 2H*Bi), and two H* on the Bi site beside the twin oxygen vacancies site (2Vo− 2H*Bi).

Fig. 5 Electrocatalysis performance of 2D Bi2O3 nanosheets in 1.0 M KOH. Polarization curves collected at a scan rate of 5 mV s−1 (a). The relationship of the overpotentials needed to reach j = 10 mA cm−2 and the relative oxygen vacancy concentrations obtained from EPR, XPS and EPR measurements (b). Tafel slopes and exchange current densities for 2D Bi2O3 nanosheets with different oxygen vacancy concentrations (c).Calculated Cdl plots of 2D Bi2O3 nanosheets at 10 mA cm−2 (d). EIS measurements on different 2D Bi2O3 nanosheets. (f) Stability test on the Pl-30 up to 50 h (e)

Fig. 6 DFT calculations on the formation and stability of Vo within Bi2O3. From left to right, hydrogen atom adsorbed Bi2O3 (010) surface models in defect-free state and with one or two oxygen vacancy incorporated (a), the calculated relative surface energies as functions of the chemical potential of oxygen (∆μO) (b), and the spin-polarized electronic density of states of Bi2O3 surfaces with different oxygen vacancies coverage (c)

As shown in Fig. 6b, the calculated surface formation free energies of all types of possible hydrogen intermediate adsorption on the active sites of the Bi2O3surface as the variation of the oxygen chemical potential (ΔμO) or the oxygen partial pressure (PO2) were presented. In this plot,the oxygen chemical potential ΔμO(the bottomx-axis) is converted into the dependence on the oxygen partial pressurePO2(the topx-axis) according to the ideal gas laws at a temperature of 800 K. Higher ΔμOandPO2correspond to lower oxygen vacancy concentration in the catalysts. The energetically most stable hydrogen adsorbed models in the allowed chemical potential zone were also highlighted with thick yellow lines. In the high oxygen chemical potential range (ΔμO> − 1.5 eV orPO2> 102Pa), where the catalyst has very low concentration of oxygen vacancies, the hydrogen intermediate is preferred to attach onto the active O site rather than the Bi sites, and the most possible adsorption model is the single H* on O site (H*O). At a low oxygen chemical potential range (ΔμO< − 1.5 eV orPO2< 102Pa), where a high oxygen vacancies concentration exists, the adsorption of single H* on twin oxygen vacancies site (2Vo− H*Vo) yields the lowest formation energy. While the high oxygen vacancy concentration leads to favorable and high number of H* adsorption, the stability of the 2Vo− H*Vosignificantly increases as the dramatic drop of the formation energy with the decrease of ∆μOIt is confirmed that theVoconcentration cannot be too high, otherwise the desorption of H* becomes unfavorable and thus deactivate the HER catalysis. This result well explains the existence of the saturation ofVo-induced activity in the 2D Bi2O3nanosheets for triggering HER.

Figure 6c shows the calculated density of states (DOS)of Bi2O3with different oxygen concentrations. A clear reduced band gap of the Bi2O3resulted from the appearance of the defect states of oxygen vacancies at the top of valence band, which contributes to improved charge transfer required for the HER.

4 Conclusion

In conclusion, 2D Bi2O3nanosheets with modulated oxygen vacancies via the customized atmosphericpressure plasma irradiation were studied for activated electrocatalytic HER in alkaline media, which cannot only balance the hydrogen intermediates adsorption but also improved the charge transfer kinetics. A linear relationship between the plasma-inducedVoconcentration and the HER performance was revealed, ascribed by the optimized H* adsorption energy with improved charge transfer, before reaching a saturation concentration. After that, the increasedVolevel resulted in trapped defects,which lowered the charge carrier density, and induced highly stable H* adsorption on the active sites, and thus decreased the catalytic efficiency. By tuning the oxygen vacancy concentration, the HER-inert Bi2O3was activated into an active HER catalyst with a low overpotential of 174.2 mV to reach 10 mA cm−2, a low Tafel slope of 80 mV dec−1, a high exchange current density of 316 mA cm−2, and excellent stability. This work, thus,paves a way to activating inherently inert metal oxide materials into high performance electrocatalysts and provides some insights into the engineering ofVofor energy conversion and storage.

AcknowledgementsThis work was financially supported by the Australian Research Council (ARC) through Future Fellowship grants (FT180100387 and FT160100281) and Discovery Projects(DP200103568, DP210100472, and DP200102546). WL thanks the support of the Science and Technology Commission of Shanghai Municipality (19520713200). The authors acknowledge the support from Central Analytical Research Facility (CARF) in QUT for material characterizations and the generous grants of CPU time from the Australian National Computational Infrastructure Facility and the High-performance Computing Centre of QUT.

Authors’ Contibution ZW designed and conducted this study;TL and ZS supervised the whole project, conceived the concept,and revised the manuscript; KO supervised the plasma processing experiments, contributed to the results interpretation, and revised the manuscript; APO provided support to the electrocatalysis test and revised the manuscript; SW, JAM, AM, WL, JY, WL, and YG provided their support to the experiments and discussed the results.

FundingOpen access funding provided by Shanghai Jiao Tong University.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons. org/licen ses/ by/4. 0/.

Supplementary InformationThe online version contains supplementary material available at https:// doi. org/ 10. 1007/s40820- 022- 00832-6.