Open Access Article
Jiahui Luoab,
Baoxin Wuac,
Kejun Yanab,
Yongbiao Mu
ab,
Qing Zhangab,
Huanxiu Zoud,
Zhiqing Tangab,
Guobin Zhange,
Wanwisa Limphiratf,
Wenjia Li*ab and
Lin Zeng
*ab
aDepartment of Mechanical and Energy Engineering, Southern University of Science and Technology, Shenzhen 518055, China
bSUSTech Energy Institute for Carbon Neutrality, Southern University of Science and Technology, Shenzhen 518055, P. R. China
cDepartment of Mechanical Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong SAR, P. R. China
dKey Laboratory of Hunan Province for Advanced Carbon-Based Functional Materials, School of Chemistry and Chemical Engineering, Hunan Institute of Science and Technology, Yueyang 414006, P. R. China
eFuture Technology School, Shenzhen Technology University, Shenzhen, 518118, China
fSynchrotron Light Research Institute, 111 University Avenue, Muang District, Nakhon Ratchasima 30000, Thailand
First published on 22nd May 2026
Industrial anion exchange membrane water electrolysis (AEMWE) still lacks earth-abundant oxygen evolution reaction (OER) catalysts that simultaneously exhibit high intrinsic activity and long-term durability. This limitation mainly originates from the conventional adsorbate evolution mechanism (AEM), which is restricted by the linear scaling relationship between *OOH and *O intermediates. Herein, we propose a rapid diffusion strategy based on molten-salt-assisted approach to synthesize P/FeOx. The entire synthesis completes within minutes via high-temperature solid-state interdiffusion, enabling single-batch production at the kilogram scale with exceptional time and energy efficiency. Phosphorus doping successfully narrows the band gap between Fe 3d and O 2p, thus enhancing electronic conductivity. It also promotes the formation of Ov and induces a mechanistic transition from the AEM to the lattice-oxygen oxidation mechanism (LOM) pathway, simultaneously reducing the activation barrier for the *OH → *O step. The optimized catalyst achieves an overpotential of 256 mV at 100 mA cm−2 and maintains stable operation for over 2000 h at 100 mA cm−2 without noticeable degradation. When integrated into an AEMWE device, the membrane-electrode assembly requires merely 1.78 V to deliver 500 mA cm−2 and operates stably for 1000 h. In conclusion, this study presents a rational synthesis approach and practical design principles for engineering high-performance OER electrocatalysts tailored to industrial water electrolysis.
Extensive studies have demonstrated that the intrinsic oxygen evolution reaction (OER) activity of an electrocatalyst is primarily determined by the reaction pathway occurring at its active sites.31,32 In alkaline media, the OER process generally proceeds through two main pathways: the adsorbate evolution mechanism (AEM) and the lattice oxygen mechanism (LOM).19,33 In the conventional AEM process, the catalytic performance is governed by the binding strength of metal sites toward oxygenated intermediates, including *OH, *O, and *OOH. However, this pathway is inherently limited by the linear *OOH–*O scaling relationship, which constrains the theoretical overpotential to approximately 0.37 V. In contrast, the LOM bypasses the *OOH formation bottleneck by coupling *O with lattice oxygen to form O2, thereby achieving a lower theoretical overpotential.34 Consequently, regulating the reaction pathway through electronic structure modulation and surface engineering provides an effective strategy to lower the overall overpotential. A key characteristic of the LOM pathway is the absence of the *OOH intermediate; lattice oxygen participates directly in forming *OO species, circumventing the potential-determining O–O coupling step and minimizing both overpotential and energy loss.35 As a result, LOM-based electrocatalysts have emerged as a central research focus in water-splitting catalysis.36–38 For instance, Huang et al.39 employed 18O differential electrochemical mass spectrometry (DEMS) and density functional theory (DFT) analyses to reveal that abundant oxygen vacancies in Fe2O3–CeO2 heterojunctions facilitate a mechanistic transition from the AEM to LOM. Similarly, Tang et al.40 demonstrated that atomic-scale incorporation of rare-earth elements into NiCo2O4 enhances M–O covalency, thereby activating lattice oxygen. Both experimental and theoretical investigations consistently indicate that a vacancy-rich transition-metal (oxy)hydroxide matrix with strengthened M–O covalency is a prerequisite for efficient lattice-oxygen oxidation in transition-metal oxides.34,41–44
Molten-salt-assisted synthesis has recently emerged as an effective strategy for modulating the morphology and microstructure of catalysts.45–47 Although phosphorus incorporation can, in principle, be realized by conventional synthetic or post-treatment methods, such routes usually involve relatively slow diffusion processes and may lead to less controllable dopant distribution, limited defect regulation, or particle aggregation. In contrast, the molten-salt flash synthesis strategy provides a highly dynamic and non-equilibrium reaction environment, where rapid ion diffusion and efficient mass transport enable simultaneous surface-confined phosphorus incorporation, oxygen-vacancy generation, and morphology regulation within a short reaction time. Therefore, the optimized catalytic behavior of P/FeOx originates from the synergistic effect between phosphorus-induced electronic/defect modulation and the unique structural regulation enabled by the molten-salt flash synthesis route.48 In this work, we developed a rapid diffusion strategy based on a molten-salt-assisted approach to construct an iron oxide catalyst with in situ phosphorus doping. With high-temperature solid-state interdiffusion, we achieved simultaneous morphology regulation and oxygen-vacancy engineering. The molten-salt-mediated rapid annealing effectively relieves lattice strain and introduces multidimensional vacancies into the Fe2O3 matrix. Phosphorus incorporation downshifts the Fe 3d conduction band, narrows the bandgap, and thereby enhances electronic conductivity. Comprehensive characterization techniques, including electron paramagnetic resonance (EPR), pH-dependent activity measurements, TMAOH probing experiments, operando 18O-labeled DEMS, and DFT calculations, suggest that phosphorus incorporation promotes lattice-oxygen activation and enhances the contribution of catalyst-derived oxygen species during the OER. Rather than indicating a complete AEM-to-LOM transition, these results support a LOM-involved pathway with partial lattice-oxygen participation in P/FeOx. The optimized catalyst achieves overpotentials of only 208 and 256 mV at current densities of 10 and 100 mA cm−2, respectively, and demonstrates exceptional durability over 2000 h. When assembled into an alkaline AEMWE device, it delivers 500 mA cm−2 at 1.78 V (60 °C) and maintains stable operation for 1000 h without observable degradation. In summary, this work establishes a molten-salt-assisted protocol that merges facile operation with kilogram-scale throughput, leveraging the high-ion-conductivity melt to collapse diffusion barriers and accelerate reaction kinetics. It delivers practical design rules for the rapid deployment of high-activity, robust oxygen-evolution electrocatalysts in industrial water electrolysis.
Scanning electron microscopy (SEM) revealed that the phosphorus-rich oxide exhibited a number-weighted particle diameter of 50–100 nm, markedly smaller than that of its undoped counterpart (100–300 nm; Fig. S1 and S2, SI). Besides, laser particle size analysis (LPSA) further confirmed that phosphorus doping effectively reduced the average particle size and suppressed grain growth (Fig. S3, SI). Powder X-ray diffraction (XRD) was performed to examine the phase structure of P/FeOx. As shown in Fig. 1b, the diffraction peaks located at 2θ values of 23.7°, 32.7°, 35.2°, 49.1°, 53.7°, 62.1°, and 63.7° can be indexed to hematite Fe2O3 (JCPDS No. 88-2359). Compared with FeOx, P/FeOx exhibits slight peak broadening and a small shift toward lower 2θ angles, suggesting increased interplanar spacing and possible lattice expansion/strain after phosphorus incorporation. The aberration-corrected high-angle annular dark-field scanning TEM (AC-HAADF-STEM) analysis of P/FeOx revealed abundant crystalline–amorphous (c–a) interfaces (Fig. 1c, S4 and S5). Besides, further analysis clearly indicated the interlaced crystalline and amorphous structure (Fig. 1d), simply due to the flash heating and cooling process of the molten-salt reaction. To study the composition of the amorphous phase, electron paramagnetic resonance (EPR) spectroscopy was further employed to probe the concentration of oxygen vacancies (VO) in the catalysts (Fig. 1e). P/FeOx displayed a significantly intensified EPR signal compared with the undoped FeOx, confirming the enhanced formation of VO induced by phosphorus incorporation. These results collectively demonstrate the superior oxygen-vacancy-generating ability of phosphorus and its critical role in boosting OER performance. When focused on the crystal phase in the catalyst, high-resolution transmission electron microscopy (HRTEM) images showed that the FeOx particles were well faceted, while phosphorus incorporation disrupted the long-range order and produced an amorphous surface layer enriched with coordinatively unsaturated Fe sites (Fig. 1f, S6 and S7, SI).
In addition to the amorphous halo (Fig. 1f), lattice fringes corresponding to the (116) and (112) planes of P/FeOx were clearly identified. The lattice spacing of the (112) facet expanded from 0.377 nm in FeOx to 0.381 nm upon phosphorus incorporation (Fig. 1g and h), evidencing lattice strain associated with phosphorus enrichment. Furthermore, the phosphorus distribution in P/FeOx was investigated by electron energy-loss spectroscopy (EELS) line-scan across an individual particle. As shown in Fig. 1i, the P L2,3-edge signal is mainly detected within the first ∼2.0 nm from the surface and gradually decreases to the noise level at greater depths, whereas the Fe L2,3-edge and O K-edge signals remain relatively constant across the scanned region. These results indicate that phosphorus species are predominantly confined to the surface/near-surface region rather than being uniformly incorporated into the bulk lattice. Such a surface-confined distribution exposes P-modified Fe–O sites directly to the electrolyte, which is favorable for regulating the near-surface electronic structure, oxygen-vacancy formation, and interfacial OER kinetics. Elemental mapping by energy-dispersive X-ray spectroscopy (EDS) further confirmed the uniform distribution of Fe, P, and O within the nanoparticles (Fig. 1j and S8, SI), indicating a homogeneous single-phase structure. Quantitative analysis revealed approximate atomic contents of 3.29 at% P, 30.40 at% Fe, and 66.31 at% O (Table S1, SI). In addition, the Raman spectra of P/FeOx and FeOx exhibited seven characteristic phonon modes of Fe2O3 at 223, 238, 291, 408, 499, 609, and 657 cm−1 (Fig. S9, SI), corroborating the above XRD results. ICP-OES was further performed to quantify the actual phosphorus content in the final catalysts. For the optimized P/FeOx sample prepared with a nominal phosphorus feeding ratio of 2 wt%, the actual phosphorus content is 0.84 wt%, corresponding to a phosphorus retention efficiency of approximately 42%. The lower actual content than the nominal feeding amount suggests that part of the introduced phosphorus is removed during molten-salt synthesis and subsequent washing. These results provide quantitative information on the retained phosphorus content in P/FeOx and clarify the relationship between the precursor feeding amount and the actual catalyst composition.
X-ray photoelectron spectroscopy (XPS) was employed to probe the surface chemical states and near-surface electronic environment of the catalysts. Therefore, the Fe 2p and O 1s spectra are discussed mainly in terms of surface electronic modulation and surface oxygen-vacancy-related species. As shown in Fig. S10, SI, the peak located at approximately 133 eV corresponds to the P 2p signal, confirming the successful incorporation of phosphorus. The Fe 2p spectra of P/FeOx and FeOx are displayed in Fig. 2a. The characteristic peaks at 712.0/725.5 eV and 718.6/732.1 eV were assigned to Fe3+ (2p3/2/2p1/2) and their corresponding satellite peaks, respectively. Compared with FeOx, the Fe 2p peaks of P/FeOx shift slightly toward higher binding energy, indicating that phosphorus incorporation modifies the local electronic environment of Fe species. This shift reflects electronic redistribution within the Fe–O–P framework. As shown in Fig. 2b, the O 1s spectra were deconvoluted into three distinct components corresponding to Olattice, Ovacancy, and Odefect. The Olattice peak at ∼529.7 eV is associated with oxygen bound within the metal oxide lattice, while the Ovacancy peak centered at ∼531.2 eV arises from surface chemisorbed oxygen species (O22− or O−) linked to defect oxides or hydroxyl groups. The Odefect component is attributed to unavoidable surface contaminants such as residual water and carbon–oxygen species from adsorbed organics. The ratio of Ovacancy/Olattice serves as a qualitative indicator of surface oxygen-vacancy concentration, where a higher ratio corresponds to a greater vacancy density. As illustrated in Fig. 2b, FeOx exhibits a relatively low Ovacancy/Olattice ratio, while the ratio increases markedly in P/FeOx, indicating that P–O bonding effectively promotes the generation of surface oxygen vacancies.
To further investigate the average electronic structure and local coordination environment of Fe atoms, X-ray absorption spectroscopy (XAS) measurements were performed. The Fe K-edge XANES spectra provide information on the average oxidation state and electronic structure of Fe, while the EXAFS fitting results reveal the average Fe–O coordination environment, including the Fe–O coordination number and local structural disorder. The normalized Fe K-edge X-ray absorption near-edge structure (XANES) spectra of P/FeOx, and FeOx closely resemble that of Fe2O3 (Fig. 2c), confirming that Fe species predominantly exist in a positively charged state with an average oxidation state near +3. The oxidation states of P/FeOx (+3.05) and FeOx (+2.83) were calculated by linearly fitting the first derivative position (E0) (Fig. S12, SI). The increased Fe valence is mainly attributed to electron aggregation around phosphorus. Fourier-transformed extended X-ray absorption fine structure (FT-XANES) spectra at the Fe K-edge (Fig. 2d) show peaks at approximately 1.54 and 2.67 Å, corresponding to Fe–O and Fe–Fe coordination, respectively, consistent with Fe2O3. Upon phosphorus incorporation, the Fe–O bond length shifts to higher R values by ∼0.062 Å, reflecting lattice expansion due to substitutional incorporation of larger P atoms for O atoms. The EXAFS fitting results (Fig. S13 and Table S2, SI) indicate Fe–O coordination numbers of 6.00, 5.85, and 5.31 for Fe2O3, FeOx, and P/FeOx, respectively. Compared with FeOx, the further decrease in the Fe–O coordination number of P/FeOx suggests that phosphorus incorporation induces a more defective and disordered local Fe–O environment. Considering that both FeOx and P/FeOx were synthesized by the same molten-salt route, this additional reduction in coordination number indicates that phosphorus doping introduces extra perturbation to the Fe–O framework during oxide formation, leading to enhanced defect generation and local structural distortion in the catalyst. Besides, the corresponding k-space fitting curves (Fig. S14, SI) further validate the reliability of these results. Wavelet-transformed (WT) EXAFS spectra of P/FeOx (Fig. 2f) display contour intensity maxima at 4.23 Å−1 and 7.89 Å−1, closely matching those of FeOx and Fe2O3, corroborating the above XAS analysis and confirming the structural integrity of the doped catalyst.
The OER performance of P/FeOx and the reference sample was evaluated in 1 M KOH electrolyte using a standard three-electrode configuration. All catalysts were ultrasonically dispersed for 1 h and then uniformly deposited on 1 × 1 cm2 nickel foam substrates at a loading of 1 mg cm−2 (Fig. 3). As shown in Fig. 3a, the iR-corrected polarization curves demonstrate that P/FeOx achieves current densities of 10 and 100 mA cm−2 at overpotentials of 208 and 256 mV, respectively (Fig. S15, SI). These values are significantly lower than that of FeOx (258 and 334 mV) and bare nickel foam (293 and 427 mV), outperforming many recently reported OER catalysts (Fig. 3c, S29 and Table S4, SI). To determine the optimal phosphorus feeding ratio, a series of P/FeOx catalysts with different nominal phosphorus contents of 1, 2, 4, 6, and 8 wt% were synthesized under otherwise identical molten-salt flash synthesis conditions. Representative SEM images show that these samples generally maintain nanoparticle morphologies, suggesting that the overall morphology is mainly governed by the molten-salt reaction environment. Their OER activities were further evaluated in 1.0 M KOH. As shown in Fig. S16 and Table S5, the 2 wt% P/FeOx sample exhibits the best catalytic performance, with an overpotential of 256 mV at 100 mA cm−2. By contrast, both lower and higher phosphorus contents result in decreased activity. This result indicates that phosphorus incorporation has an optimal range: insufficient phosphorus limits electronic/defect modulation, whereas excessive phosphorus may disturb the optimized surface structure or partially block active sites, leading to inferior OER performance. Given the intrinsically sluggish OER kinetics, Tafel slopes (Δη/Δlog|j|) were calculated to assess reaction kinetics. As shown in Fig. 3b, P/FeOx exhibits an exceptionally low Tafel slope of ∼31.6 mV dec−1, substantially outperforming FeOx (∼73.9 mV dec−1), and nickel foam (∼50.3 mV dec−1), confirming its superior kinetic characteristics. Electrochemical impedance spectroscopy (EIS) was also employed to further evaluate the charge-transfer resistance (Rct). The Nyquist plots (Fig. S15b, SI) reveal that P/FeOx possesses the smallest Rct among the tested catalysts, indicating the most efficient interfacial electron-transfer dynamics.
The electrochemical double-layer capacitance (Cdl), determined from cyclic voltammetry (CV) tests at varying scan rates (Fig. S17, SI), was used as a proxy for the electrochemically active surface area. As shown in Fig. 3d, P/FeOx exhibits a markedly larger Cdl compared with the control sample, suggesting a higher density of accessible active sites. Catalyst durability was also assessed, as stability is a key performance criterion for OER catalysts. As shown in Fig. 3e, P/FeOx demonstrates outstanding long-term stability, retaining over 97% of its initial activity after 2000 h of continuous operation at a constant current density of 100 mA cm−2. Moreover, its linear sweep voltammetry (LSV) curve after 18
000 activation cycles nearly overlaps with the initial curve (Fig. 3e), indicating negligible degradation. In contrast, FeOx exhibits an overpotential decay of ∼40 mV at 100 mA cm−2 after cycling (Fig. S18, SI). To further verify the structural robustness of P/FeOx under prolonged electrochemical operation, post-cycling characterization studies were carried out after 18
000 LSV cycles. The XRD pattern of the cycled catalyst remains essentially unchanged compared with that of the pristine sample, with no obvious impurity peaks or detectable phase transformation. Furthermore, the post-cycling STEM and TEM images (Fig. S19, SI) show that the nanostructure is largely preserved and clear lattice fringes are still maintained, further confirming the excellent structural stability of P/FeOx during long-term OER cycling. In addition, XPS analysis after cycling shows that the Fe 2p, O 1s, and P 2p spectra exhibit no significant change compared with those of the fresh sample, and the Fe3+ peaks remain nearly unchanged, indicating good stability of the surface chemical states during prolonged electrochemical operation. When benchmarked against state-of-the-art OER catalysts reported in recent literature at 10 mA cm−2, P/FeOx maintains a clear performance advantage (Fig. S29, SI). Furthermore, faradaic efficiency measurements (Fig. S20, SI) confirm that both doped and undoped catalysts achieve nearly 100% selectivity toward oxygen evolution, indicating that phosphorus incorporation does not induce parasitic side reactions and that the measured current arises exclusively from the OER process.
The synthesized P/FeOx catalyst was further employed as the anode material in a customized anion exchange membrane water electrolysis (AEMWE) device (Fig. S21, SI). Commercial Pt/C served as the cathode to form the electrode pair (Pt/C‖P/FeOx), while Pt/C‖FeOx and Pt/C‖Ni foam configurations were used as a benchmark for comparison. The electrolysis performance of the assembled AEMWE systems is shown in Fig. 3f. Polarization measurements reveal that the Pt/C‖P/FeOx cell achieves cell voltages of 1.51 and 1.76 V at current densities of 100 and 500 mA cm−2, respectively, at 60 °C. These values are significantly lower than those of the FeOx-based electrolyzer (1.63 and 1.93 V at 100 and 500 mA cm−2) and the Ni foam-based one (1.69 and 2.05 V at 100 and 500 mA cm−2), clearly demonstrating the superior catalytic performance of P/FeOx. Moreover, this superior performance surpasses that of most recently reported AEMWE systems (Table S6, SI). Temperature-dependent linear sweep voltammetry (LSV) measurements further confirm that elevated operating temperatures were capable of decreasing the required cell voltage and enhancing catalytic kinetics (Fig. S22, SI). In addition to high activity, the P/FeOx-based AEMWE exhibits excellent durability, maintaining a nearly constant voltage of ∼1.78 V over 1000 h of continuous operation at 500 mA cm−2, with less than 5% performance degradation (Fig. 4g). These results collectively demonstrate the high efficiency and long-term stability of P/FeOx as a robust anode catalyst for AEMWE, underscoring its strong potential for sustainable, large-scale hydrogen production.
In situ Raman spectroscopy was employed to track the structural evolution of P/FeOx and FeOx during the anodic polarization process. At open-circuit potential, both samples exhibit the characteristic spectrum of α-Fe2O3 (218, 286, 400, 488 and 602 cm−1). For P/FeOx, these peaks begin to attenuate at 0.40 V vs. RHE and are completely extinguished by 0.60 V (Fig. 4a). The reversible disappearance of Raman intensity is attributed to potential-driven, concerted extraction of lattice oxygen and hole injection into Fe 3d states, which transiently breaks the translational symmetry required for Raman activity without provoking long-range amorphization or cation dissolution. Upon returning the potential to <0.40 V, the original phonon manifold is quantitatively restored, confirming that the Fe–O sublattice remains intact throughout the redox cycle. In contrast, the undoped FeOx retains its full Raman signature up to 0.50 V and only begins to weaken above 0.60 V (Fig. S23, SI), demonstrating the slow reaction potential of the undoped catalyst.
In our endeavor to delve deeper into the catalyst's active source, we undertook an in-depth investigation into the fundamental mechanisms that govern its activity enhancement and functional performance. The OER generally proceeds through two primary pathways: the adsorbate evolution mechanism (AEM) and the lattice oxygen mechanism (LOM). The key distinction between these two pathways lies in the nature of the reaction intermediates. In the AEM pathway, the OER progresses via the formation of the *OOH intermediate, which typically incurs a high overpotential due to the energy-intensive O–O bond formation step. In contrast, the LOM pathway circumvents the *OOH formation by directly coupling *O species with lattice oxygen to generate *OO, thereby reducing the overall energy barrier and facilitating oxygen evolution. To experimentally validate the mechanistic transition, the pH dependence of OER activity was systematically examined. LSV measurements were conducted across a pH range of 12.5–14 (Fig. 4b). As shown in Fig. 4c, P/FeOx exhibits a markedly stronger pH dependence compared to FeOx. Quantitative analysis based on the proton reaction order (ρRHE = ∂log
j/∂pH, referenced to the reversible hydrogen electrode) yields a ρRHE value of 0.59 for P/FeOx, substantially higher than those of FeOx (0.23). This elevated proton reaction order indicates that the OER kinetics of P/FeOx are strongly influenced by proton activity, consistent with a non-concerted proton–electron transfer process characteristic of the LOM, rather than the conventional AEM pathway.
Direct experimental evidence for LOM participation was further obtained through the detection of negatively charged oxygen intermediates (*O22−), which are characteristic of lattice oxygen involvement. Tetramethylammonium hydroxide (TMAOH) was employed as a mechanistic probe, as its cation (TMA+) can electrostatically bind to these negatively charged oxygen species, competitively blocking active sites and suppressing OER activity if the LOM pathway is operative (Fig. 4d). When the electrolyte was switched from 1 M KOH to 1 M TMAOH, P/FeOx exhibited a pronounced activity loss, requiring an additional 36 mV to reach 100 mA cm−2. In contrast, the overpotentials of FeOx decreased by only 15 mV. This sharp contrast strongly supports the predominance of the LOM pathway in P/FeOx. Kinetic behavior was further analyzed by examining the redox peaks in the OER pre-regime at scan rates ranging from 1 to 5 mV s−1 (Fig. S24, SI). The relationship between peak current (i) and scan rate (v) follows the equation i = avb, where the exponent b (obtained from the slope of log i versus log
v) reflects the rate-determining process: b ≈ 0.5 indicates diffusion-controlled behavior, while b ≈ 1.0 corresponds to surface-capacitive (adsorption-controlled) behavior. Deconvolution of the b values (Fig. S25, SI) reveals that approximately 67% of the OER process on P/FeOx is governed by surface-capacitive control, with only 33% limited by diffusion. In contrast, FeOx exhibits significantly lower surface-controlled contributions (45%), indicating that P–O bonding facilitates a transition from diffusion-limited to surface-dominated kinetics. This shift aligns well with the activation of the LOM pathway. To provide mechanistic insight into the involvement of catalyst-derived oxygen species, operando 18O-labeled differential electrochemical mass spectrometry (DEMS) measurements were performed using H218O in 1 M KOH. After 18O labeling, the catalyst was rinsed with H216O electrolyte five times and then subjected to OER measurements in 1 M H216O, while the evolved oxygen species were continuously monitored by online DEMS. During repeated electrochemical cycling, weakly adsorbed or loosely bound surface 18O species would be expected to be largely consumed at the early stage of the measurement. Therefore, although a minor contribution from surface-exchangeable oxygen species cannot be completely excluded, the stronger 16O18O signal observed for P/FeOx than for FeOx provides strong evidence for the participation of catalyst-derived oxygen species during the OER and is consistent with enhanced lattice-oxygen involvement. Combined with the pH-dependent behavior, TMAOH probing experiments, and 18O isotope labeling mass spectrometry, these observations support a LOM-involved OER process with increased lattice-oxygen participation in P/FeOx, rather than an exclusive LOM pathway.
The high density of oxygen vacancies, outstanding OER performance, and reduced Tafel slope collectively indicate modified reaction kinetics and the possible involvement of the LOM pathway. The doping approach has proven effective in promoting lattice oxygen release within iron oxide. Phosphorus doping significantly reduces the energy required for the formation of lattice oxygen vacancies. Computational analysis confirmed that the introduction of a single phosphorus atom can decrease the oxygen formation energy from 3.00 eV to 2.79 eV (Fig. 4g). The schematic diagram in Fig. 4g shows that the released lattice oxygen combines with adsorbed oxygen species to form oxygen molecules, thereby completing a full cycle of the LOM pathway. This mechanism not only enhances the catalytic activity of the iron oxide but also offers new perspectives for optimizing OER catalysts.
To further elucidate the origin of the exceptional intrinsic OER activity and durability of P/FeOx, density functional theory (DFT) calculations were conducted to comprehensively investigate its reaction mechanism. To reflect the experimentally observed surface-doped character of P/FeOx, the DFT calculations were performed on a defective α-Fe2O3 (012) slab in which the P atom was introduced only into the outermost surface layer. The corresponding model contains one surface oxygen vacancy and a surface P concentration of one P atom per 109.14 Å2. This model reasonably captures the surface-confined phosphorus modification and was used to evaluate its influence on the electronic structure and OER energetics (Fig. S27, SI). In the lattice oxygen mechanism, lattice oxygen directly participates in the OER as part of the reaction intermediates, and the direct O–O coupling effectively circumvents the linear scaling relationship that constrains the conventional adsorbate evolution mechanism. The optimized structural model and the corresponding reaction pathway are illustrated in Fig. 5a. A comparative analysis of the density of states (DOS) for P/FeOx and FeOx (Fig. 5b and c) reveals that FeOx possesses a bandgap of 1.24 eV between the Fe 3d and O 2p states. Upon the co-introduction of nonmetallic phosphorus and oxygen vacancies (Ov), the valence and conduction bands begin to overlap, effectively eliminating the bandgap. This overlap signifies enhanced electronic coupling between the Fe 3d and O 2p orbitals, thereby promoting charge delocalization and improving electrical conductivity. To gain deeper mechanistic insight into the LOM-driven OER process, Gibbs free energy (ΔG) profiles for each elementary step were computed for Fe active sites (Fig. 5d). For both P/FeOx and FeOx, the potential-determining step corresponds to the dehydrogenation of adsorbed *OH to *O. Notably, phosphorus doping and oxygen-vacancy formation substantially lower the ΔG for this step from 1.15 to 0.89 eV, underscoring the crucial role of phosphorus in facilitating *OH dehydrogenation.
To further probe the electronic origin of this enhancement, adsorption configurations of OH* on FeOx and P/FeOx were analyzed through differential charge density mapping (Fig. 5e). In FeOx–OH*, a charge of 0.32 eV is transferred from the catalyst to the OH* adsorbate, whereas this value increases to 0.35 eV in P/FeOx–OH*. Since the oxygen atom in OH* is directly bonded to the active Fe site, the increased electron density localized on this oxygen enhances electrostatic repulsion toward the positively charged H*, thereby promoting H* desorption. The strength of the O–H bond was further evaluated using projected crystal orbital Hamilton population (pCOHP) analysis. The integrated COHP value for the O–H bond increases from −7.86 to −7.77 upon phosphorus and oxygen-vacancy incorporation (Fig. 5f), indicating a weakened O–H bond and a reduced energy barrier for dehydrogenation. Moreover, d-band center calculations (Fig. 5g and S28, SI) reveal that phosphorus and Ov co-doping downshifts the O 2p band center from −2.70 to −3.26 eV. This downward shift was consistent with the electron-donating characteristics of phosphorus and oxygen vacancies. The enhanced O 2p antibonding state further weakens the generated O–H bond, promoting the cleavage of the O–H bond during the dehydrogenation process.
800 rpm for 3 min. The pellet was re-dispersed and washed twice more to eliminate residual NaNO3 and chloride. The final precipitate was dried overnight at 60 °C under dynamic vacuum (<10−2 mbar) to yield the red P/FeOx powder.
FeOx was prepared identically by omitting NaH2PO4 entirely, while all other parameters remained unchanged.
| This journal is © The Royal Society of Chemistry 2026 |