2D arrays of hollow carbon nanoboxes: outward contraction-induced hollowing mechanism in Fe–N–C catalysts

Maximizing the utilization efficiency of monatomic Fe sites in Fe–N–C catalysts poses a significant challenge for their commercial applications. Herein, a structural and electronic dual-modulation is achieved on a Fe–N–C catalyst to substantially enhance its catalytic performance. We develop a facile multi-component ice-templating co-assembly (MIC) strategy to construct two-dimensional (2D) arrays of monatomic Fe-anchored hollow carbon nanoboxes (Fe-HCBA) via a novel dual-outward interfacial contraction hollowing mechanism. The pore engineering not only enlarges the physical surface area and pore volume but also doubles the electrochemically active specific surface area. Additionally, the unique 2D carbon array structure reduces interfacial resistance and promotes electron/mass transfer. Consequently, the Fe-HCBA catalysts exhibit superior oxygen reduction performance with a six-fold enhancement in both mass activity (1.84 A cm−2) and turnover frequency (0.048 e− site−1 s−1), compared to microporous Fe–N–C catalysts. Moreover, the incorporation of phosphorus further enhances the total electrocatalytic performance by three times by regulating the electron structure of Fe–N4 sites. Benefitting from these outstanding characteristics, the optimal 2D P/Fe-HCBA catalyst exhibits great applicability in rechargeable liquid- and solid-state zinc–air batteries with peak power densities of 186 and 44.5 mW cm−2, respectively.


Introduction
Electrochemical oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) assume a pivotal role in the eld of sustainable energy conversion and storage devices, such as zinc-air batteries (ZABs). 1,2To date, platinum group metal (PGM) catalysts are used as the most efficient catalysts for ORR and OER. 3,4However, suffering from their scarcity and unsustainability, it is of great signicance to develop non-PGM ORR catalysts with high activity and durability.[10][11] Additionally, the overall electrocatalytic performance is determined not only by the intrinsic activity of active sites but also by the number and utilization efficiency of accessible active sites.1][22] Therefore, it is imperative to design hierarchically porous structures incorporating both hollow cavities and/or highly open nanoarchitectures for Fe-N-C catalysts.
Hollow carbon is typically prepared using the hardtemplating method, which involves employing templates such as silica spheres (SiO 2 ), 23 polystyrene (PS) spheres, 24 ZnO spheres, 25,26 or CdS spheres. 27This synthesis strategy allows for the preparation of either monodisperse zero-dimensional (0D) hollow carbon nanoparticles or three-dimensional (3D) macroporous carbon materials.For example, using SiO 2 as a template and ZIF as a carbon source, 0D hollow carbon nanoparticles or 3D interconnected carbon frameworks can be successfully synthesized. 28,29Furthermore, 3D-ordered macroporous carbon nanocrystals can be prepared through the PS colloidal crystal templating method. 30However, these hard-templating methods oen involve complex synthesis processes and environmentally unfriendly template removal steps.Achieving two-dimensional (2D) monolayered hollow carbon materials, with a thickness equivalent to that of a single nanoparticle, remains a signicant challenge due to the difficulty in the self-assembly of template nanoparticles into 2D monolayered assemblies. 31lternatively, ZIF nanocrystals can serve as sacricial templates to prepare hollow carbon nanoarchitectures, such as 0D hollow carbon polyhedra, 32 one-dimensional (1D) carbon nanobers, 33 and 3D carbon aerogels, 34 through an outward contraction stress-induced synthesis approach.In this process, another carbon source (e.g., polyvinyl pyrrolidone (PVP), polyacrylamide (PAN), polydopamine (PDA)) is coated on the surface of ZIF nanocrystals using various coating methods to selectively form thin carbon shells.These shells induce the decomposition of ZIF nanocrystals, resulting in the formation of large cavities. 35However, the fabrication of ZIF-templated 2D monolayered hollow carbon nanoarchitecture has been rarely reported and is mainly hindered by the difficulties in assembling ZIF nanocrystals and additional carbon sources into 2D monolayers. 36everaging a comprehensive understanding of the mechanism underlying hollowing in ZIF-derived hollow carbon and capitalizing on the advantages of ice-templating self-assembly strategy for large-scale fabrication of 2D metal-organic framework monolayer superstructures, we meticulously engineer 2D arrays of hollow carbon nanoboxes anchored with monatomic Fe (Fe-HCBA).Through structural and electronic dualmodulation, we signicantly enhance its catalytic performance to meet the requirements for application in ZABs.A facile and efficient multi-component ice-templating coassembly (MIC) approach is utilized to prepare 2D Fe 3+ /PVP/ ZIF-8 hybrid monolayers, which are subsequently transferred into 2D Fe-HCBA via a dual-outward contraction stress-induced hollowing mechanism (Fig. 1a).During the pyrolysis process, ultrathin caron walls are preferentially formed around ZIF-8 polyhedra in the monolayer, inducing outward contraction stress (F outward 1).Simultaneously, as the carbon wall develops, the initially isolated ZIF-8 units approach each other closely, resulting in a second outward contraction stress (F outward 2) originating from the interconnected facets.By harnessing the synergy between F outward 1 and F outward 2 forces, robust 2D Fe-HCBA composed of well-arranged hollow carbon nanoboxes with ultrathin carbon walls and highly porous carbon frames in the cavities are constructed.Furthermore, the incorporation of Fe 3+ facilitates the generation of abundant Fe-N 4 active sites on both carbon walls and frames.The outward contraction hollowing mechanism-based pore engineering not only enlarges the physical surface area and pore volume to maximize the exposure of active sites but also increases the electrochemically accessible surface area by 2.7 times to promote utilization efficiency.The unique 2D monolayered carbon array architecture effectively promotes electron transfer, reduces interfacial resistance, and facilitates mass transfer in the three-phase reaction, thereby exhibiting excellent catalytic performance.Moreover, the doping of phosphorus regulates the electron structure of Fe-N 4 active sites, further enhancing the intrinsic activity of the P/Fe-HCBA catalyst and leading to satisfactory performance in rechargeable ZABs with both liquid and solidstate electrolytes.

Fabrication and structural characterization
The MIC of ZIF-8 colloids, PVP, and Fe 3+ with varying mass ratios results in the formation of a series of hybrid Fe 3+ /PVP/ZIF-8 monolayers, designated as Fe 3+ /PVP/ZIF-8-x (for details, see ESI †).Low magnication scanning electron microscope (LM-SEM) images demonstrate the large-scale 2D thin nanosheet morphology of Fe 3+ /PVP/ZIF-8-2 (Fig. 1b and S1a †).Scanning electron microscope (SEM) and transmission electron microscope (TEM) images show that truncated rhombic dodecahedral ZIF-8 nanocrystals are well-arranged and coated with PVP in a monolayer of Fe 3+ /PVP/ZIF-8-2 (Fig. 1c and S1b-g †).Elemental mapping images show that the Fe element is uniformly distributed on the surface and in the gap region among ZIF-8 polyhedra (Fig. S1h †).Subsequently, pyrolysis of Fe 3+ /PVP/ZIF-8-2 was conducted in a N 2 atmosphere at 900 C to obtain Fe-HCBA-2, and its 2D large nanosheet-like morphology is well maintained (Fig. S2a †).SEM images exhibit the carbon monolayer of well-arranged carbon nanoboxes (Fig. 1d, e and S2b-d †).TEM images clearly show that hollow carbon nanoboxes with distinct carbon walls are tightly interconnected in the monolayered architectures (Fig. 1f, g and S2e-g †).Scanning TEM (STEM) image provides further evidence for the ultrathin carbon walls and highly porous carbon frames within the Fe-HCBA-2 (Fig. 1h).Elemental mapping conrms the homogeneous dispersion of Fe, C, and N elements in Fe-HCBA-2 (Fig. 1i).Notably, atomic-resolution HAADF-STEM image of Fe-HCBA-2 (Fig. 1j) illustrates that the bright dots on the carbon matrix can be identied as atomically dispersed Fe single atoms (SAs).
To understand the evolutionary process of hollow architecture, Fe-HCBA-1 and Fe-HCBA-3 were prepared from Fe 3+ /PVP/ ZIF-8-1 and Fe 3+ /PVP/ZIF-8-3, respectively.The successful coassembly of ZIF-8 with increasing amounts of PVP is evidenced by the gradual reduction in peak intensity observed in powder X-ray diffraction (PXRD) patterns and decreased Brunauer-Emmett-Teller (BET) surface areas calculated from N 2 sorption measurements of Fe 3+ /PVP/ZIF-8-x (Fig. S3 and S4 †).SEM and TEM images for Fe-HCBA-1 and Fe-HCBA-3 are presented in Fig. S5 and S6.† A dual-outward contraction stressinduced hollowing mechanism is illustrated in Fig. 2a.At low PVP loading amounts, the in situ formed carbon layer may exhibit insufficient thickness to provide a robust interfacial interaction with ZIF-8 polyhedron, resulting in a weak outward contraction stress F outward 1, leading to an underdeveloped hollow structure of Fe-HCBA-1 (Fig. 2b).In this case, F outward 2 predominantly governs the hollowing process of Fe-HCBA-1.
With increasing PVP content, the thickened carbon walls offer enhanced interaction and signicantly amplify F outward 1, which synergistically contributes along with F outward 2 to welldeveloped hollow architectures with distinct carbon walls and inner carbon frames (Fig. 2c).Upon further increase in PVP content, the overall outward contraction stress is further augmented while maintaining a constant thickness of the carbon wall; however, the inner carbon frames are almost absence due to the full decomposition of ZIF-8.This results in the formation of 2D arrays of fully hollow carbon nanoboxes (Fig. 2d).This observation strongly suggests that the evolution of the Fe-HCBA-x is governed by the synergy between F outward 1 and F outward 2. The temperature-dependent hollowing process was also investigated.With an increase in pyrolysis temperature from 500 to 900 C, the decomposition of ZIF-8 polyhedra gradually occurs starting from the center point, resulting in the formation of highly porous carbon frames within the preferentially formed carbon nanoboxes (Fig. S7 †).This observation further conrms the dual-outward contraction stress-induced hollowing mechanism.
The sample of ZIF-derived microporous N-doped carbon anchored with Fe SAs (Fe-NC) was prepared according to the detailed procedure provided in the ESI (Fig. S8).† The PXRD patterns (Fig. S9 †) demonstrate that the broad peak at approximately 24 in all Fe-HCBA-x and Fe-NC samples corresponds to the (002) plane of graphitic carbon, while no peaks attributed to metallic Fe nanoparticles are detected.This observation indicates that the Fe species exist as atomically dispersed Fe SAs, consistent with the atomic-resolution HAADF-STEM characterization.The pore characteristics and BET surface areas of Fe-HCBA-x and Fe-NC were analyzed through N 2 sorption measurements, with detailed information provided in Table S1.† Fe-NC displays a typical type-I isotherm, indicating its micropore-dominated pore nature.In contrast, the isotherms for all Fe-HCBA-x samples show large hysteresis loops (Fig. 3a) and indicate the coexistence of micropores, mesopores, and macropores.The BET surface area of Fe-HCBA-2 is signicantly higher at 1570 m 2 g À1 , surpassing that of Fe-NC (670 m 2 g À1 ) by 2.3 times and exceeding those of Fe-HCBA-1 (1530 m 2 g À1 ) and Fe-HCBA-3 (1082 m 2 g À1 ).This variation in BET surface area is not only correlated with the micropore contents but also with the highly porous carbon in the hollow cavities.Additionally, the pore size distribution analysis based on nonlocal density functional theory (NLDFT) reveals that Fe-HCBA-x samples exhibit a signicant abundance of mesopores and macropores in the range of 10-100 nm, while Fe-NC only shows micropores (Fig. 3b).The Fe content in the Fe-HCBA-x and Fe-NC samples was determined using coupled plasma optical emission spectrometry (ICP-OES) (Table S2 †).Given a constant Fe content in Fe 3+ /PVP/ZIF-8-x precursors, the notable increase in Fe from 0.55 wt%, 0.72 wt%, to 1.06 wt% for Fe-HCBA-1, -2 and -3 is attributed to the difference in carbon mass aer pyrolysis (Fig. S10 †), providing further evidence of the degree of hollowness exhibited by these three Fe-HCBA-x samples.
X-ray photoelectron spectroscopy (XPS) measurements were conducted to further investigate the elemental composition and chemical state of the Fe-HCBA-x and Fe-NC catalysts.The highresolution Fe 2p spectra (Fig. S11 †) indicate that all catalysts contain two pairs of peaks for Fe 2+ (711.0 and 723.6 eV) and Fe 3+ (715.5 and 728.4 eV), suggesting the electron-decient characteristics of these cationic Fe species.Meanwhile, the highresolution N 1s spectra of all catalysts can be divided into ve peaks at 398.6, 399.5, 400.5, 401.3 and 403.2 eV, which originate from pyridinic-N, Fe-N, pyrrolic-N, graphitic-N and oxidized-N, respectively (Fig. 3c).The relative atomic percentages of all catalysts (Fig. S12a and Table S3 †) show a positive correlation between the nitrogen content and the enhanced degree of hollowness in Fe-HCBA-x samples, which can be attributed to the extent of ZIF-8's decomposition.In addition, Fe-HCBA-2 has the highest content of pyridinic-N atoms (23.60%) (Fig. S12b †), which not only serve as the active sites but also coordinate with Fe atoms to form Fe-N x active sites, thereby improving the ORR activity of the catalyst. 37o precisely identify the local atomic and electronic structures of Fe species in Fe-HCBA-2, synchrotron X-ray absorption near-edge structure (XANES) and extended X-ray absorption ne structure (EXAFS) spectroscopy were performed.As revealed in the Fe K-edge XANES spectra (Fig. 3d), the absorption edge of Fe-HCBA-2 lies around that of FePc and between those of the standard Fe foil and Fe 2 O 3 , indicating that the positive charge of Fe species in Fe-HCBA-2 is between 0 and +3.The Fouriertransformed EXAFS R-space plot in Fig. 3e displays that only a predominant peak of Fe-HCBA-2 is located at 1.53 Å, which is in line with the Fe-N rst coordination sphere.Meanwhile, there are no peaks attributed to Fe-Fe and Fe-O-Fe in Fe-HCBA-2, and the FT-EXAFS tting results (Fig. S13 and Table S4 †) reveal a Fe-N 4 coordination environment of Fe-HCBA-2.Moreover, the EXAFS wavelet transform (WT) spectra in Fig. 3f illustrate that a maximum k value of 5.5 ÅÀ1 in Fe-HCBA-2 is identied as Fe-N(O) path, which can be distinguished from the Fe-Fe (Fe foil, 8.1 ÅÀ1 ) and Fe-O-Fe (Fe 2 O 3 , 7.8 ÅÀ1 ) paths.
Collectively, these ndings provide compelling evidence supporting the existence of atomically dispersed Fe-N 4 sites in Fe-HCBA-2.

Electrocatalytic performance
Inspired by the desirable Fe-N 4 active sites and structural advantages offered by the 2D array architecture, the electrocatalytic activity of Fe-HCBA-x and Fe-NC catalysts was explored.The ORR activity was rst evaluated by linear sweep voltammetry (LSV) curves using a rotating disk electrode (RDE) in 0.1 M KOH at 1600 rpm.Fig. 4a indicates that Fe-HCBA-2 exhibits the most positive onset potential (E onset , 0.978 V) compared to Fe-HCBA-1, Fe-HCBA-3, Fe-NC and Pt/C catalysts.Notably, Fe-HCBA-2 shows a superior half-wave potential (E 1/2 ) of 0.894 V, outperforming those of Fe-NC (0.830 V), Fe-HCBA-1 (0.849 V), Fe-HCBA-3 (0.878 V), and Pt/C (0.876 V).As shown in Fig. S14a, † the kinetic current density (J k ) was extracted to evaluate the activity and mass transfer performance.As expected, the Fe-HCBA-2 exhibits an excellent J k of 10.6 mA cm À2 at 0.88 V, which is 7.6, 5.9, 2.0, and 2.1 times higher than those of Fe-NC, Fe-HCBA-1, Fe-HCBA-3, and Pt/C, respectively.Furthermore, the enhanced reaction kinetics of Fe-HCBA-2 was veried by its Tafel slope (51 mV dec À1 ) in Fig. 4b, which is much lower than that of Fe-NC (76 mV dec À1 ).Considering that Fe-NC and Fe-HCBA-1 possess similar Fe content, it is noteworthy that despite their comparable J k values, the Tafel slope of Fe-NC surpasses that of Fe-HCBA-1 signicantly.This observation suggests a substantial hindrance in the utilization efficacy of Fe SAs active sites within the microporous structure of Fe-NC, consequently leading to sluggish kinetics during the ORR reaction.In addition, despite having a higher content of Fe-N 4 sites and an almost entirely hollow structure, the ORR performance of Fe-HCBA-3 is inferior to that of Fe-HCBA-2.This observation highlights the crucial role played by the highly porous carbon frames within the hollow nanoboxes, which can signicantly enlarge the BET surface area and pore volume, thereby facilitating the exposure of Fe-N 4 active sites to reactants.The electrochemical active specic surface area (ECSA) serves as a credible metric for evaluating the surface properties of catalysts in an electrolyte, which can be deduced through the computation of double-layer capacitance (C dl ) values (Fig. S15 †).As depicted in Fig. 4c, the calculated C dl of Fe-HCBA-2 (34.7 mF cm À2 ) is 2.7 times higher than that of Fe-NC (12.9 mF cm À2 ), aligning well with the observed enlargement of physical BET surface area by a factor of 2.3.This result highlights the superior ECSA and enhanced accessibility of catalytic active sites for Fe-HCBA-2.
Electrochemical impedance spectroscopy (EIS) tests were employed to investigate electron transfer and interfacial kinetics properties (Fig. S14b †).Impressively, Fe-HCBA-2 shows a signicantly lower charge transfer resistance (R ct , 67.8 U) than that of Fe-NC (128.1 U), indicating that the facilitated electron transport efficiency in 2D array architecture can be attributed to its well-arranged particle conguration.In constract, the random stacking conguration of Fe-NC results in higher impedance for electron transfer. 38Furthermore, the boosted macroscopic reactant diffusion performance is evidenced by the quantied diffusion coefficients of hydroxide ion (D OH À) in the plot of Z 0 against u À1/2 . 36,38As shown in Fig. 4d, the calculated D OH À of Fe-HCBA-2 (1.67 Â 10 À14 cm 2 S À1 ) is ten times higher than that of Fe-NC (1.21 Â 10 À15 cm 2 S À1 ), while exhibiting similar values to other 2D array catalysts such as Fe-HCBA-1 (2.06 Â 10 À14 cm 2 S À1 ) and Fe-HCBA-3 (1.25 Â 10 À14 cm 2 S À1 ).This signicant enhancement in hydroxide ion transfer efficiency highlights the indispensable role played by the 2D array architectures in improving electrocatalytic performance.To garner a more profound understanding of the utilization of active sites, the mass activity (MA) and turnover frequency (TOF) are derived by normalizing the kinetic current density at 0.88 V with respect to the corresponding Fe quantity (Fig. S14c †).As summarized in Fig. 4e, Fe-HCBA-2 presents a signicantly enhanced value in MA of 1.84 A cm À2 and TOF of 0.048 e À site À1 s À1 , which are approximately 6-fold higher than those of Fe-NC, respectively.This highlights the effectiveness of the meticulously designed 2D monolayered interconnected hierarchically porous carbon nanobox architecture in maximizing the utilization of exposed active sites, facilitating electron and mass transfer, and enhancing O 2 diffusion and transportation, thereby demonstrating the signicant impact of structural modulation effect.
Considering that phosphorus (P) doping is one of the most effective strategies for enhancing the intrinsic activity of Fe-N 4 moieties by modulating their electronic structure toward both ORR and OER, the P element was further incorporated into Fe-HCBA-2 to obtain P/Fe-HCBA-2. 39The characterization results of P/Fe-HCBA-2 are shown in Fig. S16-S18.† The 2D array architecture of hollow carbon nanoboxes remains well-preserved aer phosphating treatment, demonstrating excellent structural integrity.The uniform distribution of the P element is observed throughout the carbon skeleton, indicating successful P doping.Additionally, the high-resolution P 2p spectrum of P/ Fe-HCBA-2 displays three distinct peaks corresponding to the P-C (132.3 eV), P-N (133.2 eV), and P-O (134.1 eV) groups, respectively. 9Furthermore, the analysis of multiple N species in P/Fe-HCBA-2 reveals that apart from iron-nitrogen species (16.06%), pyridinic-N (22.79%) and graphitic-N (30.20%) are the predominant species, which are conducive to the electronic conduction and structural stability. 40he LSV curves in Fig. 4a show that the E 1/2 and limit current density (J L ) of P/Fe-HCBA-2 are 0.910 V and 6.3 mA cm À2 , respectively, surpassing those of Fe-HCBA-2 (E 1/2 , 0.894 V; J L , 5.5 mA cm À2 ) and other advanced Fe-N-C catalysts (Table S5 †).The radar map in Fig. 4f provides a comparative analysis of ORR performance indices for P/Fe-HCBA-2 and Fe-HCBA-2.Specically, the J k and Tafel slope of P/Fe-HCBA-2 are superior to those of Fe-HCBA-2.More importantly, the MA and TOF metrics of P/ Fe-HCBA-2 have respectively escalated to be 4.94 A cm À2 and 0.128 e À site À1 s À1 , which are approximately three-fold higher than those of Fe-HCBA-2.These ndings suggest that the incorporation of P atoms exerts a pivotal inuence in modulating the coordination environment of the Fe-N 4 active sites, and thus enhancing the intrinsic activity of the catalyst.
DFT calculations were performed to further explore the inuence of electronic modulation in the process of ORR.According to experimental measurements, the reasonable structure of a P-doped catalyst (Fe-N 4 -P) should consist of 4 N atoms coordinating with 1 Fe atoms on the graphene sheet, and the P atom is coordinated with N and C atoms in the second shell (Fig. S19a †).For comparison, the Fe-N 4 model was also considered (Fig. S19b †).To unveil the origin of the observed enhancement in catalytic activity and accelerated reaction kinetics of ORR for Fe-N 4 -P, Gibbs free energy diagrams of ORR were investigated.As presented in Fig. 4g, at U ¼ 1.23 V, the desorption process of *OH is recognized to be the ratedetermining step. 41Compared to that on Fe-N 4 , the overpotential of ORR on Fe-N 4 -P was lower by 0.06 V. To further probe the inuence of coordination microenvironment of Fe-N 4 -P on the ORR process, the electronic structures of comparative models were investigated by calculating the difference charge densities and Bader charges of the metal active centers (Fig. 4h).Charges of 1.02 e À are transferred from the Fe atom to graphene over Fe-N 4 -P, lower than that of Fe-N 4 (1.08 e À ), which veries that the P atoms can promote electron accumulation in Fe sites.Furthermore, it is noted that the P doping is favorable to the electron donation from P atoms to the Fe sites, promoting the electron density in the Fe centers and thus enhancing intermediate absorption over the Fe sites.The DFT calculation results, in conjunction with the experimental ndings, provide compelling evidence that P doping signicantly enhances the catalytic performance of monatomic Fe.
To gain deeper insights into the perfect retention of the monatomic Fe-N 4 active sites in P/Fe-HCBA-2, a KSCN poisoning experiment was performed.Upon the injection of KSCN solution, P/Fe-HCBA-2 displays a remarkable negative shi of 43 mV in E 1/2 and a decrease in J L to 3.7 mA cm À2 , proving compelling evidence for the primary contribution of these monatomic Fe-N 4 sites (Fig. S20 †).Then, the rotating ring disk electrode (RRDE) measurements were carried out to determine the ORR pathways.Particularly, the peroxide yield of P/Fe-HCBA-2 is below 1% over the potential range of 0-0.8 V (Fig. S14d †), and its electron-transfer number (n) can be calculated as 3.99, indicating an excellent 4e À selectivity during the ORR process (Fig. S21 †).Fig. S22 † presents the current retention situation derived from the chronoamperometric measurements to evaluate the electrocatalytic stability.Compared to Pt/ C (84.7%), P/Fe-HCBA-2 exhibits remarkable stability (94.7%), highlighting its potential as a highly efficient and cost-effective alternative to traditional platinum-based catalysts.Moreover, the catalyst exhibits unparalleled robustness, as evidenced by the absence of discernible morphological or structural changes following the rigorous durability test (Fig. S23 †).The remarkable methanol tolerance ability of P/Fe-HCBA-2 was further substantiated by the absence of a signicant decline in current density aer methanol addition (Fig. S24 †).
Subsequently, the electrocatalytic OER activities of the catalysts were investigated in 1.0 M KOH electrolyte solutions.The LSV curves indicate that P/Fe-HCBA-2 exhibits an overpotential of 360 mV at 10 mA cm À2 (E j ¼ 10), which is comparable to the benchmark RuO 2 (322 mV) and much smaller than that of Fe-HCBA-2 (421 mV) (Fig. S25a †).Further exploration of OER kinetics in Fig. S25b † shows that the Tafel slope of P/Fe-HCBA-2 is determined to be 110 mV dec À1 , which is smaller than those of RuO 2 (130 mV dec À1 ) and Fe-HCBA-2 (163 mV dec À1 ), respectively.This highlights the signicant enhancement in the sluggish OER kinetics of Fe-HCBA-2 resulting from the incorporation of phosphorus atoms.Generally, the energy gap (DE) between the E 1/2 for ORR and the E j ¼ 10 for OER is adopted to estimate the bifunctional catalytic activity. 42As illustrated in Fig. S25c †, the P/ Fe-HCBA-2 shows an exceptional DE of 0.68 V, surpassing other reported bifunctional catalysts summarized in Table S6, † thereby indicating its prominent reversible electrocatalytic performance.The comprehensive electrochemical measurements unequivocally validate that the synthesis strategy, involving the construction of a 2D array architecture coupled with pore engineering and heteroatom doping, can simultaneously optimize intrinsic activity, active site utilization efficiency, and electron/mass transfer to achieve high catalytic performance.

Rechargeable liquid and solid ZABs
To evaluate the practical applicability of P/Fe-HCBA-2 in sustainable energy storage devices, we rst assembled the rechargeable liquid ZABs.Fig. 5a depicts the structural schematic diagram of the liquid ZABs, featuring P/Fe-HCBA-2, Fe-HCBA-2, and Pt/C + RuO 2 benchmark catalysts as the air cathode for comparative analysis.Consistent with the electrochemical test results, the performance of P/Fe-HCBA-2-based liquid ZABs surpasses that of Fe-HCBA-2 and Pt/C + RuO 2 benchmark catalysts across various metrics.As illustrated in Fig. 5b-e, liquid ZABs equipped with P/Fe-HCBA-2 cathode shows the highest open-circuit potential (OCP) of 1.497 V, the lowest chargedischarge voltage gap (x) especially at relatively high current density, the largest peak power density of 186 mW cm À2 , and optimal rate capability at different discharging current densities.In addition, the specic capacity of the rechargeable liquid ZABs was calculated by normalizing the Zn consumption at various discharging current densities.Encouragingly, the P/Fe-HCBA-2based liquid ZABs exhibit a maximum specic capacity of 800.2 mA h g À1 and a corresponding highest energy density of 1024.3W h kg À1 at 10 mA cm À2 (Fig. 5f), surpassing the performance of Pt/C + RuO 2 cathode even at 25 mA cm À2 (Fig. S26 †).Additionally, durability tests were conducted to assess battery output performance through charge-discharge cycles at constant current (Fig. 5g-i).Remarkably, the P/Fe-HCBA-2-based ZABs show a lower x of 0.92 V and higher energy utilization efficiency (3) of 54.1% compared to the benchmark Pt/C + RuO 2 catalyst for the 401 st cycle; furthermore, a negligible decrease of x and 3 is observed even at the 1000 th cycle.Continuous cycling tests at varying current densities were also conducted (Fig. S27 †), revealing excellent durability and reversibility of this rechargeable liquid ZABs.Motivated by the demands of miniaturized and portable electronic devices, rechargeable exible solid-state ZABs were fabricated using carbon cloth coated with P/Fe-HCBA-2 as the air cathode and the polyacrylic acid/KOH/Zn(CH 3 COO) 2 hydrogel polymer as the electrolyte (Fig. 6a).Flexible ZABs equipped with Pt/C + RuO 2 cathode were also assembled.As depicted in Fig. 6b, the P/Fe-HCBA-2-based exible ZABs deliver a higher OCV of 1.46 V compared to the Pt/C + RuO 2 -based one (1.41V).Of particular interest, at bending angles of 0 , 45 , 90 , 135 , 180 and back to 0 , its OCV values remain consistent at levels of 1.463 V, 1.475 V, 1.474 V, 1.462 V, 1.462 V, and nally reaching again at a value of approximately 1.466 V (Fig. S28 †).Additionally, the practicality of the P/Fe-HCBA-2-based exible ZABs is further validated by their integration into a series circuit.The two-battery assembly of the exible ZABs maintains a high OCV ranging from 2.88 to 2.90 V at various bending angles (Fig. 6b and S29 †), which is approximately twice that of a single exible battery, demonstrating its capability to power a series of LEDs (Fig. S30 †) and highlighting its exceptional exibility and versatility.As shown in Fig. 6c and d, P/Fe-HCBA-2-based exible ZABs display a peak power density (44.5 mW cm À2 ) and specic capacity (777.1 mA h g À1 ), surpassing those achieved with Pt/C + RuO 2 cathode (21.7 mW cm À2 and 654.8 mA h g À1 , respectively), as well as the most recently reported ZABs (Table S6 †).To comprehensively evaluate the exibility of the P/Fe-HCBA-2-based exible ZABs under external stress, we conducted in-depth investigations into their cycling stability under various bending angles.The x and 3 of the exible ZABs are well maintained aer 125 cycles at a current density of 1 mA cm À2 compared to the Pt/C + RuO 2 benchmark (Fig. 6e and S31 †).Notably, this exible ZABs shows consistent charge-discharge processes for up to 150 cycles at a current density of 2 mA cm À2 (Fig. S32 †), highlighting its immense prospects for portable electronic devices application.These results unequivocally demonstrate that the P/Fe-HCBA-2 catalyst possesses remarkable potential for the development of nextgeneration energy conversion devices.

Conclusions
In summary, we have fabricated 2D Fe 3+ /PVP/ZIF-8 hybrid monolayers with well-arranged ZIF nanopolyhedra coated with PVP and Fe 3+ species using a facial multi-component icetemplating co-assembly approach.These monolayers serve as precursors for the construction of 2D monolayered arrays of hollow carbon nanoboxes via a dual-outward contractioninduced hollowing mechanism.The validity of this hollowing mechanism is demonstrated by harnessing the synergistic interfacial interaction forces.The structurally and electronically dually-modulated P/Fe-HCBA-2 catalyst exhibits remarkable catalytic performance towards both ORR and OER, which can be attributed to several fundamental strengths: (i) hollow carbon nanoboxes with ultrathin carbon walls and highly porous carbon frames provide large physical and electrochemical active specic surface area that maximizes the exposure and utilization efficiency of Fe-N 4 active sites; (ii) the unique 2D monolayered array structure promotes electron/mass transfer efficiency; (iii) P-doping can signicantly enhance the intrinsic activity of Fe-N-C catalyst by modulating electronic structure of Fe-N 4 sites, thus enhancing the reaction kinetics.This work demonstrates a universal approach to rationally design advanced 2D monolayered hollow monatomic metal-nitrogen-carbon catalysts (M-N-C), while also emphasizing the signicance of optimizing utilization efficiency and intrinsic activity of active sites, as well as facilitating electron/mass transfer of catalysts to achieve superior overall electrocatalytic performance.

Fig. 6
Fig. 6 Performance evaluation of rechargeable flexible ZABs.(a) Schematic illustration of the rechargeable flexible ZABs.(b) Open-circuit voltage (inset: the stable voltage provided by two P/Fe-HCBA-2-based flexible ZABs in series).(c) Discharge polarization curves and corresponding power density plots.(d) Galvanostatic discharge curves at 1 mA cm À2 .(e) The discharge/charge cycling stability of flexible ZABs with P/ Fe-HCBA-2 cathode at 1 mA cm À2 under different bending states.