The SF, SFN, D-SFN (S_3−yFN_x, x = 0, 0.2; y = 0, 0.2), and BZCYYb materials were prepared by the sol-gel method [25]. For S_3−yFN_x, Sr(NO₃)₂ (AR, ≥ 99.5%), Fe(NO₃)₃·9H₂O (AR, ≥ 98.5%), C₁₀H₅NbO₂₀ (AR, ≥ 98%) as the chemical reagents were stoichiometrically dissolved in deionized water on a magnetic heating stirrer at 80 °C until gel formation occurred. The citric acid (AR, ≥ 99.5%) and EDTA (AR, ≥ 99.5%) were added as complexing agents, and ammonia water (25%-28%) was employed to make the solution neutral or slightly alkaline. Subsequently, the gel was heated in the oven at 180 °C for 10 h to obtain the precursor, which was then calcined at 1100 °C for 10 h to obtain the initial S_3-yFN_x powders.

An X-ray diffractometer (XRD, Rigaku SmartLab 9 kW) was employed to record diffraction patterns of the powder samples within the 10°-80° ranges using a step-scan mode with 0.02° intervals. X-ray photoelectron spectroscopy (XPS, Thermo Fisher Scientific Nexsa) was utilized to investigate the elemental valence states on the S_3−yFN_x surface. A field emission scanning electron microscope (SEM, Tescan MIRA) and scanning transmission electron microscopy (STEM, FEI Talos F200x) were used to observe the cross-sectional microstructure of the single cell. A high-resolution transmission electron microscope (HR-TEM, FEI Talos F200x) was employed to observe the lattice spacing of fresh powder samples, and energy-dispersive X-ray spectroscopy (EDX, FEI Talos F200x) mapping was used to obtain element distributions. Inductively coupled plasma optical emission spectroscopy (ICP-OES, Thermo Scientific iCAP 7600) was employed for precise compositional analysis. The oxygen vacancy content of S_3−yFN_x samples at room temperature was determined via iodometric titration. Detailed experimental information is available in the supporting information. Thermogravimetry analysis (TGA, TGA5500) was conducted from 40 to 700 °C. The average linear thermal expansion coefficient (TEC) of the S_3−yFN_x samples was measured using a dilatometer (DIL 402CL, Netzsch) in air, with a heating rate of 10 °C min⁻¹. The fourier transform infrared (FT-IR, Thermo Scientific Nicolet iS5) spectra were used to examine the hydration degree of the samples, with a scanning range from 1500 to 4000 cm⁻¹. Before FT-IR test, the samples were subjected to a 5-h drying process in a 100 °C oven to eliminate adsorbed water from the sample surfaces. The electrical conductivity of the S_3−yFN_x was examined by the DC four-probe method. The specific experimental details can be found in the supporting information.

Electrochemical impedance spectroscopy (EIS) test was done from 10⁵ to 10^-1 Hz by an electrochemical workstation (Solarton 1287 + 1260), with a signal amplitude of 30 mV. I-V-p curve and Ni-BZCYYb|BZCYYb|S_3−yFN_x single cell stability tests were measured by a digital source meter (Keithley 2440). The single cell was supplied with 50 sccm of H₂ as fuel and 100 sccm of synthetic air as oxidant. Water vapor was introduced via a water bath heating method. Furthermore, the distribution of relaxation time (DRT) method was employed to comprehensively understand the EIS results [26,27]. Details about DRT can be found in the supporting information.

The first-principles computations were conducted within the framework of DFT using VASP [28,29]. The exchange-correlation interaction was treated employing a generalized gradient approximation (GGA) as described by the Perdew-Burke-Ernzerhof (PBE) function [30]. The kinetic energy cutoff of the plane wave utilized to extend the Kohn-Sham electron wave function was set to 400 eV, with an iterative convergence of energy being 10^-5 eV. All atomic positions were allowed to relax until the Hellmann-Feynman force fell below 0.01 eV Å⁻¹. The Brillouin zone was sampled by a 2 × 2 × 2 k-point grid. Further specifics of the calculation are elucidated in the supporting information.

To illustrate the respective effects of Nb doping and A-site cation defects on the material properties, three perovskite precursors with different compositions were designed, namely, the initial Sr₃Fe₂O_7-δ (SF) material, the Nb-doped Sr₃Fe_1.8Nb_0.2O_7-δ (SFN), and the simultaneous Nb-doped and A-site deficient Sr_2.8Fe_1.8Nb_0.2O_7-δ (D-SFN). The successful synthesis and crystallinity of SF, SFN, and D-SFN compositions via the sol-gel method were verified using XRD, as shown in Fig. 2a. All patterns were indexed to the standard Ruddlesden-Popper A₃B₂O_7-δ perovskite structure (PDF#01-082-0427), with no detectable impurities. This confirms the incorporation of Nb into the SF lattice for SFN and D-SFN. In perovskites, the BO₆ octahedron structurally supports the framework, enabling A-site Sr vacancies in D-SFN without instability. Further Rietveld refinement quantified the lattice parameters (Fig. S1, Tables S1-S4). Excellent refinement fitting (R_P < 10%, R_WP < 15%) validated the reliability. The successful doping of Nb and the introduction of A-site defects are further supported by the partially expanded lattice after Nb doping (301.4 ų for SF compared to 305.5 ų for SFN) and the subsequent partially contracted unit cell of D-SFN (302.7 ų). STEM revealed uniform nanoparticle morphologies for D-SFN with diameters of ~ 200-300 nm (Fig. 2c). HR-TEM verified the crystallinity through measured d-spacings of 0.2641 and 0.2737 nm, matching the (112) and (110) planes (Fig. 2d, e). EDX confirmed the expected elemental composition (Fig. 2b). Given that EDX is a semi-quantitative analysis method, the elemental composition of fresh D-SFN samples was further determined using ICP-OES. The results indicate that the proportions of Sr, Fe, and Nb in the cationic total are 50.63%, 20.04%, and 4.00%, respectively, closely aligning with the nominal values of 51.49%, 21.10%, and 3.90%, with only minor deviations observed. EDX mapping further verified the successful homogeneous substitution of Nb into the parent SF lattice (Fig. 2f).

For the air electrode of RePCC, the excellent electronic conductivity helps to realize the electronic conducting path and provides the necessary conditions for the electrochemical reaction to proceed [31]. Herein, the electronic conductivity of S_3−yFN_x materials was investigated using the DC four-probe method within a temperature range of 250-750 °C (Fig. 3a). The electrical conductivity profile of SF material exhibits remarkable temperature-dependent behavior, initially increasing with temperature likely due to thermal excitation of charge carriers, followed by a noticeable decline around 400 °C. This non-linear behavior indicates a transition from semiconducting to metallic characteristics, which is due to the lattice oxygen loss at elevated temperatures according to the following defect reaction [32]:

As temperature rises, the reduced Fe facilitates the escape of lattice oxygen, forming oxygen vacancy (${V}_{O}^{\cdot \cdot }$) and releasing O₂. Concurrently, the electron holes (${h}^{\cdot }$) were consumed, causing the decreased conductivity. Incorporating Nb led to a substantial reduction in electrical conductivity compared to the pristine SF material. Over 250–750 °C, the conductivity of SFN was only 2.8–4.46 S cm⁻¹, much lower than the 23.04–39.06 S cm⁻¹ range exhibited by SF. This considerable decrease suggests Nb⁵⁺ ions inhibit electron mobility in the material [33].

However, combining Nb doping and A-site deficiency in D-SFN notably enhanced its electrical conductivity compared to SFN, doubling it to 4.79-9.73 S cm⁻¹ over 250-750 °C. This indicates that introducing Sr vacancies can compensate for the reduced conductivity induced by Nb doping. To comprehend the mechanisms contributing to the enhanced conductivity in D-SFN, an exploration of the electronic structures of SFN and D-SFN was undertaken utilizing density of states (DOS) plots. After calibration, 0 eV corresponds to the Fermi level (E_f). As depicted in Fig. S2, the conduction band near the Fermi level is predominantly shaped by Fe 3d orbitals, while the valence band primarily originates from O 2p orbitals. This underscores the vital role of Fe-O interactions in influencing the electronic conductivity of the material. The broadening of unoccupied states around the valence band in D-SFN implies an augmentation in charge transfer [34,35]. Introducing Sr defects leads to an increase in the Fe valence states, resulting in an amplified overlap of Fe 3d and O 2p orbitals near the Fermi level. This indicates a reinforcement of Fe-O covalency, favorably impacting electron conduction in the material [36,37]. Nonetheless, the conductivity remains inferior to that of pristine SF, highlighting the intricate interplay between doping and deficiency in influencing material properties. However, when the air electrode meets the electronic conductivity of 1 S cm⁻¹, it is enough to provide the required electron transport path for electrochemical reactions [38,39]. Therefore, the electronic conductivities of the three materials can meet the needs of high-performance air electrodes.

In addition to sufficient electronic conductivity, an ideal RePCC air electrode should also present high O²⁻/H⁺ conductivity and surface exchange performance. These properties are usually related to the B-site transition metal oxidation state, oxygen vacancy concentration, and material surface chemical properties [40]. Therefore, XPS was introduced to get information about surface chemistry and oxidation states in SF, SFN, and D-SFN. The Nb 3d peaks observed around 209 eV in the SFN and D-SFN survey scans (Fig. S3) confirm the successful doping of Nb into the parent SF material. Deconvoluted high-resolution XPS spectra of Fe 2p reveal a mixed Fe³⁺/Fe⁴⁺ valence state in all samples (Fig. 3b) [41,42]. Compared with SF, the Fe valence state of SFN is significantly reduced, while the Fe valence state of D-SFN is slightly increased. Relative Fe³⁺ and Fe⁴⁺ percentages were adopted to estimate δ of the surface, presented in Table S5. Nb doping induces a decrease in the δ from 0.48 in SF to 0.41 in SFN, indicating Nb incorporation reduces surface oxygen vacancies. However, A-site deficiency in D-SFN marginally increases δ to 0.49 and is even slightly higher than that of SF.

It is generally accepted that the air electrode of RePCC, in a humid atmosphere, undergoes hydration reactions by the combination of oxygen vacancies on the electrode surface with steam, resulting in the generation of proton defects [43,44], as indicated by the following equation:

Therefore, the surface oxygen vacancy content is generally indicative of the material's hydration performance partially. It is anticipated that D-SFN, possessing a higher concentration of oxygen vacancies due to Sr deficiency, would exhibit superior hydration performance compared to SFN. FT-IR testing was employed to ascertain the hydration characteristics of SFN and D-SFN samples. As shown in Fig. S4, both samples revealed infrared absorption peaks corresponding to ${{\text{OH}}}_{{\text{O}}}^{\cdot }$ in the range of 3200–3700 cm⁻¹ after complete hydration (3% H₂O-air for 150 h). Similar phenomena have been observed in other perovskite air electrode materials, such as BaCo_0.7Fe_0.2Zr_0.1O_3−δ [45] and BaCo_0.7Ce_0.24Y_0.06O_3−δ [46]. The ${{\text{OH}}}_{{\text{O}}}^{\cdot }$ absorption peak in D-SFN is more pronounced than in SFN, indicating enhanced hydration capabilities.

Considering that lattice oxygen also participates in the ORR/WOR reaction of RePCC, excellent bulk oxygen migration properties are also significant to the ORR/WOR activity of the material [47,48]. The bulk oxygen vacancy content of the materials was determined through the iodometric method. The average oxygen vacancy contents of SF, SFN, and D-SFN are marginally higher than the surface oxygen vacancy contents, measuring 0.51, 0.46, and 0.52, respectively. Discrepancies between XPS and iodometric results could stem from variations in the valence state of the Fe element between the bulk and surface regions of the materials [49]. Compared with Fe ions in the bulk phase, Fe ions on the surface are more fully in contact with air and are more easily oxidized [50].

In addition, the lattice oxygen activity of the materials at high temperatures was investigated based on TG. Shown in Fig. S5, the pronounced mass loss observed above ~ 400 °C in all samples has been attributed to the thermal reduction of Fe and subsequent release of lattice oxygen [51,52], which confirms the change of the conductivity in Fig. 3a.

Since RePCCs usually operate at intermediate temperatures of 400-700 °C, integrating iodometry and TGA results yielded the most representative oxygen vacancy content of the air electrode material in the operational state. The oxygen vacancy evolution of SF, SFN, and D-SFN in the corresponding operating temperature ranges was depicted in Fig. 3c. Compared to SF, SFN has a lower δ value at 400-700 °C. For example, at 550 °C, the δ value for SF is 0.81, whereas for SFN, it is only 0.77. Considering that the bond strength of the Nb-O bond is significantly higher than that of the Fe-O bond, the lattice oxygen combined with Nb is more stable than that combined with Fe [53,54]. Therefore, after doping Nb, the overall activity of lattice oxygen in the bulk phase is weakened, O²⁻ ions are less likely to escape from the lattice when the temperature increases and the oxygen vacancy concentration decreases. D-SFN has the highest δ value of 1.32 at 700 °C, which may be related to the increased lattice oxygen activity near the Sr defects.

The thermal compatibility between the air electrode and electrolyte material is crucial for the practical application of a single cell. Mismatched TEC between electrolyte and electrode can induce residual stresses, leading to interfacial delamination [55]. The thermal expansion behavior of the materials under an air atmosphere was examined using thermal expansion tests, as shown in Fig. S6. The average TEC values for SF, SFN, and D-SFN from 100 to 1000 °C were 18.3 × 10^-6, 16.5 × 10^-6, and 17.1 × 10^-6 K⁻¹, respectively, significantly lower than the TEC of some typical Co-based air electrodes [56,57]. Despite Sr defects causing unit cell contraction, partially aiding in the improvement of the TEC, the material's high-temperature expansion behavior is primarily attributed to the reduction of B-site elements and lattice oxygen loss [58,59]. At elevated temperatures, the TEC of D-SFN is slightly higher than SFN, possibly due to the greater reduction of Fe⁴⁺ in D-SFN at high temperatures, as evidenced by TG and XPS results. However, it should be noted that SFN and D-SFN exhibit very close TEC within the operational temperature range (500-750 °C), highlighting the beneficial impact of Nb incorporation. Considering the relatively low average TEC of BZCYYb (9.5 × 10^-6 K⁻¹) [60], the introduction of Nb can decrease the TEC of the SF-based air electrode, enhancing the long-term thermal compatibility with BZCYYb electrolyte during operation.

To further evaluate the oxygen surface exchange and bulk diffusion properties of the electrode material in a real environment, the chemical bulk diffusion and surface exchange kinetics of the material when the oxygen partial pressure changes were tested through ECR experiments. The experimental process has been given in detail in the experimental section. Table S6 summarizes the bulk oxygen diffusion coefficient (D_chem) and surface oxygen exchange coefficient (k_chem) of SF, SFN, and D-SFN between 500 and 650 °C quantified by the ECR curve (Fig. S7). The Arrhenius curves of k_chem, and D_chem of the three materials at 500-650 °C were plotted in Fig. 3d. Compared with SF, SFN exhibits lower D_chem and k_chem values between 500 and 650 °C, indicating that niobium doping may adversely affect bulk and surface oxygen transport [61]. The introduction of A-site defects in SFN can significantly improve the oxygen transport properties of the material. Compared to SFN, the D_chem, and k_chem of D-SFN increased by 1.9 and 1.7 times respectively, at 600 °C. The oxygen vacancies introduced by Sr defects and the improved lattice oxygen activity effectively improve the bulk diffusion and surface exchange kinetics of SFN [62]. These ECR results highlight the efficacy of A-site vacancies in enhancing the oxygen transport capability of niobium-doped S_3−yFN_x compositions.

Based on the above discussion, the mechanism by which niobium incorporation and strontium defects manipulate material properties can be represented by the schematic diagram in Fig. 3e and the following Kröger-Vink notation:

Nb doping:

A-site deficiency:

Nb⁵⁺ substitution on Fe sites consumes oxygen vacancies (${V}_{O}^{\cdot \cdot }$) for charge compensation (Eq. 3). Introducing Sr deficiency regenerates vacancies and partially elevates Fe valence via reactions (Eqs. 4, 5). Thereby, A-site vacancies counterbalance the Nb-induced vacancy reduction, modulating bulk transport and surface reactivity.

The compatibility of electrode candidate materials with electrolyte components is critical for long-term stable cell operation. To examine the chemical compatibility of SF-based materials and BZCYYb electrolytes, D-SFN and BZCYYb were mixed and co-fired at 1100 °C for 10 h. The powder XRD results (Fig. S8) confirmed the absence of any new phases, indicating excellent compatibility.

Different from the ORR in SOFC mode and OER in SOEC mode in the air electrode of traditional SOCs, due to the participation of water vapor, the air electrode of RePCC performs more complex ORR in PCFC mode and WOR in PCEC mode [51]. For illustration, Fig. 4a shows a schematic diagram of the ORR reaction occurring at the air electrode in PCFC mode, and the reverse process is WOR. The introduction of steam generates proton defects by promoting the hydration reaction between surface oxygen vacancies and water vapor, causing the three-phase boundary to expand from the original electrode-electrolyte interface under dry air (Fig. S9) to encompass the entire air electrode surface (Fig. 4a), playing a key role in enhancing ORR/WOR dynamics [63,64].

To reveal the electrochemical performance of S_3-yFN_x candidates as RePCC air electrodes, EIS tests were conducted on symmetrical cells under different air humidity (0, 3, and 6 vol%) and temperature (500-700 °C), and the Arrhenius curves of area specific resistance (ASR) are shown in Figs. 4b and S10. As the steam partial pressure in the air increases, the ASR and electrode reaction activation energy (E_a) of the three electrode materials are significantly optimized, mainly because hydration generates proton defects, which effectively improves the bulk proton conductivity of the materials.

To further distinguish the ORR/WOR reaction processes of the D-SFN air electrode, the EIS curve (Fig. 4c) was obtained under various vapor pressures. The DRT method was employed to discern specific electrochemical sub-processes, leveraging the characteristic frequency response (10⁵-10⁻¹ Hz) of the EIS curve [65,66], as illustrated in Fig. 4d. Typically, low frequency (LF, 10^-1-10¹ Hz) peaks signify gas diffusion or surface adsorption/desorption steps. Mid-frequency (MF, 10¹-10⁴ Hz) peaks correlate to surface oxygen exchange and/or bulk diffusion processes. High frequency (HF, 10⁴-10⁵ Hz) peaks represent charge transfer reactions. The peak area in a specific frequency range represents the polarization resistance of the corresponding electrochemical substep. It is noted that at different vapor partial pressures, the MF and LF peak areas constitute the main part of the entire DRT curve, indicating that gas diffusion, surface exchange, and bulk diffusion resistance are the main components of polarization resistance. After the introduction of water vapor, the MF and LF resistance dropped significantly, indicating that the introduction of water vapor significantly improved the diffusion of water, surface mass transfer, and bulk ion transport performance of the electrode [40].

In addition, the electrochemical performance of the three electrode materials was compared under the conventional operating atmosphere of the air electrode (3 vol% H₂O-air) to evaluate the impact of Nb doping and A-site defects on the ORR/WOR reaction of the electrode. As shown in Fig. 4e, at 500-700 °C, the electrode reaction activation energy of SFN and D-SFN decreases significantly, indicating that the Nb substitution is beneficial to low-temperature operation. However, the incorporation of Nb reduces the electrochemical activity of SF electrodes, as can be seen from the increased polarization resistance. The further introduction of A-site Sr defects enables the D-SFN electrode to exhibit optimal ORR/WOR performance in the entire temperature range. At 550 °C, the ASR values of SF, SFN, and D-SFN were 1.781, 2.885, and 1.209 Ω cm² respectively. To provide a more objective evaluation of the performance of the D-SFN electrode, Fig. 4f and Table S7 compare the ASR of D-SFN with that of high-performance air electrodes reported in the literature. The ASR of D-SFN is even significantly lower than many Co-based air electrodes, demonstrating outstanding electrocatalytic performance.

Besides, combined with DRT to further separate the electrode reaction process, as shown in Fig. 4g, the incorporation of Nb has adverse effects on both the surface mass transfer and the bulk diffusion process. The further introduction of Sr defects not only eliminates the adverse effects of Nb doping on the electrode reaction but also significantly reduces the resistance of all electrochemical sub-steps. The performance improvement can be attributed to the fact that the D-SFN material has more surface and bulk oxygen vacancies and presents the best surface exchange and bulk diffusion rates, which are key properties of the ORR/WOR reaction [8].

The ORR reaction process based on a four-electron mechanism further elucidates the ORR reaction mechanisms on the surfaces of three air electrodes [24,67]. The ORR reaction-free energies on the surfaces of SF, Nb-substituted SFN, and A-site deficient D-SFN are depicted in Fig. 4h. The ORR on the air electrode surface comprises four steps. According to the ORR reaction barrier, it is evident that the rate-determining step for SF, SFN, and D-SFN electrode is located at the H₂O formation step. For SF, the overpotential (η) is 1.95 V, for SFN, the overpotential is 2.12 V, and D-SFN exhibits a low η of 0.74 V, which is favorable for the rapid reaction kinetics of the air electrode.

Long-term polarization resistance tests were performed on symmetrical cells to further evaluate the stability of the D-SFN electrode material. To facilitate the explanation of the effects of Sr defects and Nb doping, the stability of SF and SFN materials was also checked. The formation energies of the three materials were initially calculated (Fig. 5a) and the calculation details can be found in the supporting information. Both SFN and D-SFN exhibited lower E_form than SF, indicating that the incorporation of Nb into the Fe site enhances the stability of the perovskite. The EIS tests were then conducted in 3% H₂O-air at 550 °C for 150 h, and the ASR evolution was shown in Fig. 5b. SF electrode exhibited continuously increasing ASR indicating performance degradation over time, while Nb-doped SFN and D-SFN demonstrated negligible ASR changes, highlighting remarkably stable operation. This result confirms the the variation in E_form depicted in Fig. 5a, signifying that the introduction of Nb enhances the stability of the material.

To probe the origin of the performance divergence during long-term tests, the phase structure of the three materials in the air was analyzed. First, S_3−yFN_x powders were treated at 550 °C for 150 h in dry air to explore the phase structure evolution. As shown in Fig. 5c, all materials maintain the initial phase structure, indicating that SF-based materials can maintain long-term phase structure stability in dry air. Phase stability was further explored under humid air for 150 h (Fig. 5d). SF material underwent significant phase decomposition, forming a large amount of Sr₃Fe₂(OH)₁₂ (PDF#01-072-2242), which could be attributed to the extensive hydration-induced decomposition that SF underwent [68,69]. The flexible SrO layer enables water incorporation, disrupting the parent lattice [70,71,72]. In contrast, only a small amount of the second phase appears in niobium-doped SFN and D-SFN, confirming their excellent structural stability.

The D-SFN powder, subjected to prolonged wet treatment, was collected for TEM analysis. As depicted in Fig. 5e, the surface of the particles reveals the (115) crystal plane spacing of the D-SFN phase and the (111) crystal plane spacing of the SFH phase. This suggests that the SFH phase forms on the surface of D-SFN after prolonged treatment, consistent with the results in Fig. 5d. Additionally, a comparison of the Fe 2p XPS spectra of D-SFN before and after wet treatment (Fig. 5f) demonstrates an increase in Fe³⁺ and a decrease in Fe⁴⁺, resulting in a decrease in the average oxidation state of Fe from + 3.56 to + 3.45. According to the literature, in the SFH phase, Fe ions coordinate directly with OH-, exhibiting a + 3-oxidation state [70]. The reduction in the average oxidation state of Fe after hydration also indirectly confirms the formation of the SFH phase.

EIS and DRT analysis (Fig. 5g, h) of the deteriorating SF electrode over time provides further evidence of continuously increasing LF and MF peak resistances. Since these revelations correlate to surface and bulk processes, excessive hydration likely hinders oxygen exchange and diffusion by occupying lattice oxygen sites and vacancies [73]. While initially enhancing performance, unrestrained hydration can thus degrade long-term stability. In contrast, Nb-doped D-SFN maintains high activity alongside exceptional steam stability by stabilizing the lattice oxygen to prevent undesirable hydroxide formation. Ultimately, elucidating the intricate effects of composition and hydration on long-term RePCC electrode performance provides invaluable guidance for designing highly stable and active air electrode materials.

To evaluate the promising D-SFN air electrode candidate material in a practical RePCC device, a Ni-BZCYYb fuel electrode-supported single cell was fabricated comprising a ~ 23 μm electrolyte and ~ 10 μm D-SFN air electrode (Fig. S11). A single cell with the same configuration using SF as the air electrode was also prepared, which has the same electrolyte thickness as the single cell with the D-SFN air electrode (Fig. S11). D-SFN based fuel cell (PCFC) i-V-P curves (Fig. 6a) exhibit exceptional peak power densities (PPD) of 596, 483, 361, 242, and 165 mW cm⁻² from 650 to 450 °C, respectively, using H₂-humidified air feeds. This significantly outperforms the same configuration of a single cell with the SF air electrode (Fig. S12), directly confirming the optimized ORR electrocatalytic performance of the D-SFN air electrode. The PPDs surpass numerous reported PCFCs employing advanced cathodes (Fig. 6b, Table S8), and the cell operates stably for approximately 142 h under 0.8 V (Fig. 6c), highlighting the potential of D-SFN for air electrodes in PCFCs.

Besides, RePCCs can not only work in PCFC mode, but also PCEC mode, so it is vital to investigate the performance of the D-SFN air electrode candidate in electrolysis mode (PCEC). Under electrolysis mode at 1.3 V (Fig. 6d), high current densities of − 0.17 to − 1.19 A cm⁻² were achieved from 450 to 650 °C, exhibiting much higher current densities than electrolysis cell with the SF air electrode (Fig. S12). The performance is even comparable to cells with top-tier air electrodes (Fig. 6e, Table S9). Excellent 109-h stability at 1.3 V (Fig. 6f) demonstrates negligible degradation. The i-V curve and electrolysis stability test prove the excellent OER activity and durability of the D-SFN air electrode.

In addition, considering the working scenario of RePCCs combined with renewable energy, it is critical to examine the reversibility of the D-SFN air electrode, because the cell needs to switch between PCFC and PCEC modes intermittently [74]. The cell underwent alternating cycles of operation in PCFC mode (at 0.8 V) and PCEC mode (at 1.3 V) for 4 h each, totaling 20 cycles and lasting 160 h (Fig. 6g). Despite repeated oxidation and reduction, only slight performance loss occurred in the cell, verifying the excellent stability of the electrode and cell reversibility.

By extensively evaluating the optimized D-SFN air electrode material in a single-cell configuration, this work substantiates its high activity for ORR/WOR electrocatalysis along with excellent stability for practical RePCC devices. The exceptional outputs highlight the efficacy of the co-substitution approach in developing optimized electrode materials.