Inducing Fe 3d Electron Delocalization and Spin-State Transition of FeN4 Species Boosts Oxygen Reduction Reaction for Wearable Zinc–Air Battery

Highlights The strong interaction between Ti3C2Sx and FeN4 species induces the central metal Fe(II) in FeN4 species with intermediate spin state transferred to high spin state, in which the latter is favorable to initiate the reduction of oxygen. This strong interaction induces a remarkable Fe 3d electron delocalization with d band center upshift, boosting oxygen-containing groups adsorption on FeN4 species and oxygen reduction reaction kinetics. The resulting FeN4–Ti3C2Sx with FeN4 moieties in high spin state exhibits high half-wave potential of 0.89 V vs. RHE and high limiting current density of 6.5 mA cm−2, enabling wearable zinc–air battery showing a good discharge performance with a maximum power density of 133.6 mW cm−2. Supplementary Information The online version contains supplementary material available at 10.1007/s40820-023-01014-8.

• The resulting FeN 4 -Ti 3 C 2 S x with FeN 4 moieties in high spin state exhibits high half-wave potential of 0.89 V vs. RHE and high lim- iting current density of 6.5 mA cm −2 , enabling wearable zinc-air battery showing a good discharge performance with a maximum power density of 133.6 mW cm −2 .
ABSTRACT Transition metal-nitrogen-carbon materials (M-N-Cs), particularly Fe-N-Cs, have been found to be electroactive for accelerating oxygen reduction reaction (ORR) kinetics.Although substantial efforts have been devoted to design Fe-N-Cs with increased active species content, surface area, and electronic conductivity, their performance is still far from satisfactory.Hitherto, there is limited research about regulation on the electronic spin states of Fe centers for Fe-N-Cs electrocatalysts to improve their catalytic performance.Here, we introduce Ti 3 C 2 MXene with sulfur terminals to regulate the electronic configuration of FeN 4 species and dramatically enhance catalytic activity toward ORR.The MXene with sulfur terminals induce the spin-state transition of FeN 4 species and Fe 3d electron delocalization with d band center upshift, enabling the Fe(II) ions to bind oxygen in the end-on Shengmei Chen, Xiongyi Liang, and Sixia Hu contribute equally to this work.
To date, massive efforts have been developed to incorporate FeN 4 species, improve their density, and design geometric structures for FeN 4 species exposed [24,25].Though previous works have proved that regulation on electronic structure of FeN 4 species such as manipulation of the metal center, atomic vacancy, and edge defect can improve their activity, the electron configuration of Fe(II) in such FeN [26,27].Some researchers suppose both the low spin and high spin states are favorable to catalyze ORR reaction because an empty 3d z2 orbital or occupied by a single electron enables Fe(II) ions to bind oxygen in the end-on adsorption mode which is more readily to absorb and desorb related reaction intermediates and thus facilitate ORR catalysis.Some researchers consider that the intermediate spin state with fully filled 3d z2 orbital prevents the end-on adsorption of oxygen on Fe(II) ions, which hinders the reduction of oxygen on this site [26][27][28][29].
However, the exact correlation between the electronic spin configuration of the active site and the ORR activity has remained inadequately understood, which impedes the rational design of high-performance ORR catalysts.Thus, directly regulating the electronic spin states of Fe centers for FeN 4 species' to study the ORR catalytic activity is necessary.To achieve this, it is promising to have a catalyst support with good electronic conductivity and mechanical stability for FeN 4 species dispersion and capable of inducing an interaction with the active species [30].
MXenes, as a new kind of 2D materials, are fabricated by selective extraction of A layer from M n+1 AX n phase (where M stands for Ti, V, Nb, Mo, etc.; A represents Al, Si, Ga, etc.; X is C or N; n = 1-3) [31][32][33].In light of electronegativity, surface hydrophilicity, good mechanical stability, and electronic conductivity, MXenes emerge as a fascinating catalyst support [31][32][33].In particular, the easy surface tunability of MXenes enables the easy tailoring ability of near-MXene surface environment and their supported catalyst [31].Thus, it is attractive to apply MXenes as FeN 4 species support.However, to modify the electronic structure of FeN 4 active species, further tailoring the MXenes structure is necessary.Incorporating sulfur terminal has been proved to adjust the surface polarities and electronic properties of carbon-based materials successfully [34][35][36].Therefore, we suppose that incorporating S in the MXenes terminal can also tailor the electronic properties of MXenes, enabling the modification of the electronic structure on the supported FeN 4 active species owing to the strong interaction between the active species and the MXenes support.
Herein, we realize a significant improvement for the intrinsic ORR activity of Fe-N-Cs electrocatalyst via introducing the sulfur-terminated Ti 3 C 2 MXene as the support to disperse the iron-nitrogen species.The resulting catalyst is denoted as FeN 4 -Ti 3 C 2 S x .We perform X-ray absorption fine spectroscopy (XAFS) and X-ray photoelectron spectroscopy ) confirmed via ultraviolet photoelectron spectroscopy (UPS), electron spin resonance (ESR) spectroscopy, temperature-dependent magnetic susceptibility (M − T) measurements, and density functional theory (DFT) calculations.A d z2 orbital occupied by a single electron enables their Fe(II) ions to bind oxygen in the end-on adsorption mode which is more readily to absorb and desorb related reaction intermediates and thus facilitates ORR catalysis.Besides, the remarkable Fe 3d electron delocalization with d band center upshift can optimize the orbital hybridization of Fe 3d with p orbital of oxygen-containing groups, boosting oxygen-containing groups adsorption on FeN 4 species and ORR kinetics.Our catalyst exhibits remarkable catalytic performance enhancement with positively shifted of 80 mV for half-wave potential compared with the one without sulfur terminals, named FeN 4 -Ti 3 C 2 .Furthermore, its catalytic activity is comparable to that of commercial 20% Pt-C for half-wave potential and limiting current density (half-wave potential: 0.89 V vs. RHE for FeN 4 -Ti 3 C 2 S x , 0.88 V vs. RHE for Pt-C; limiting current density: 6.5 mA cm −2 for FeN 4 -Ti 3 C 2 S x , 5.5 mA cm −2 for Pt-C) and its long-term stability is superior to that of commercial Pt-C.Besides, integrating this FeN 4 -Ti 3 C 2 S x catalyst into a wearable ZAB shows a good discharge performance with a maximum power density of 133.6 mW cm −2 and a high cycling stability with 110 h at 2 mA cm −2 , demonstrating the feasibility of FeN 4 -Ti 3 C 2 S x in ZAB applications.

Preparation of Ti 3 C 2 MXene and FeN 4 -Ti 3 C 2 S x 2.2.1 Synthesis of Ti 3 C 2 MXene
The Ti 3 C 2 MXene was obtained using a wet chemical method reported by Gogotsi [37].Typically, 1 g of LiF powder was dissolved in 10 mL of 9 M HCL via stirring for 30 min at room temperature.Then, 1 g of Ti 3 AlC 2 MAX powder was added slowly into the above solution under ice bath condition, stirred the sealing mixed solution continuously for 24 h at 35 °C, washed this mixture solution with DI water 10 times, sonicated the resulting sediment in DI for 20 min, and then centrifuged at 3500 rpm for 20 min, and the obtained black colloidal supernatant was the Ti 3 C 2 MXene solution.The Ti 3 C 2 MXene was obtained via freeze drying the black colloidal supernatant for 48 h.

Synthesis of FeN 4 -Ti 3 C 2 S x
For FeN 4 -Ti 3 C 2 S x synthesis, 1 mmol of FeCl 3 and 3 mmol of KSCN were dissolved in 25 mL of DI water via vigorous stirring for 30 min.Then, 50 mg of the as-prepared Ti 3 C 2 MXene dissolved in 25 mL enthanol was added into the above solution and stirred for another 30 min; then sonicated the mixed solution for 1 h; after that, added 2 mmol 2, 2'-bipyridine in the above mixture solution and continuously stirred for 20 h at 25 °C until getting a dark red slurry; and next, collected the slurry and dried at 80 °C in vacuum.Then, the dried precursor was ground in quartz mortar and annealed at 950 °C at a heating rate of 5 °C min −1 under Ar atmosphere for 2 h.Finally, the annealed product was soaked in 0.5 M H 2 SO 4 at 80 °C for 8 h to remove inactive iron species.The leached sample was washed to neutral with water and enthanol for three times and dried in vacuum at 60 °C for 24 h.The FeN 4 -Ti 3 C 2 reference sample was obtained via the similar processes without adding 3 mmol of KSCN.

Materials Characterization
The methodology and structure of the prepared materials were evaluated via using scanning electron microscopy (SEM) (Philips XL30 FEG) and transmission electron microscopy (TEM, FEI Tecnai G2 F30 with 300 kV of accelerating voltage).The crystal structure was analyzed by using X-ray diffraction (Bruker, D2 Phaser) with Cu Kα (λ = 1.5418Å) radiation and Raman spectroscopy (Renishaw in Via™ confocal Raman microscope with a excitation laser of 514 nm wavelength).The specific surface area was derived from the N 2 adsorption-desorption isotherms obtained with a Micrometric ASAP 2020 instrument.Besides, the Barrett-Joyner-Halenda (BJH) method was applied to obtain the pore size distribution.The chemical composition was evaluated by using X-ray photoelectron spectroscopy (XPS) (VG ESCALAB 220i-XL).The chemical coordinated information was obtained by X-ray absorption fine spectroscopy (XAFS) spectra conducted at the beamline 1W1B of Beijing Synchrotron Radiation Facility (BSRF) at Institute of High Energy Physics, Chinese Academy of Sciences.The storage rings of BSRF were conducted at 2.5 GeV with an average current of 250 mA.The data collection was conducted in transmission mode using ionization chamber when using Si (111) double-crystal monochromator.The data were processed and analyzed similar to previous procedures via using ATHENA and ARTEMIS for X-ray absorption near-edge structure (XANES) and the extended X-ray absorption fine structure (EXAFS) spectra, respectively [38].Ultraviolet photoemission spectroscopy (UPS) measurements were carried out on an ESCA LAB 250 Xi spectrometer with He I resonance lines (21.2 eV).Temperature-dependent magnetic susceptibility (M-T) measurements were conducted in the temperature range from 10 to 300 K with a physical property measurement system model 6000 (Quantum Design, USA).Electron spin resonance (ESR) spectra were obtained by an ER200-SRC-10/12 (Bruker, Germany) spectrometer at 300 K.

Electrochemical Measurements
All electrochemical measurements were conducted on a CHI 760E electrochemical workstation integrating a rotating ring disk electrode (RRDE) in a three electrodes system, in which a glassy carbon electrode (GCE) with diameter (3 mm) loaded with catalyst as working electrode, while Pt sheet and Ag/AgCl (3 M KCl) were used as counter electrode and reference electrode, respectively.The recorded potential was converted to reversible hydrogen electrode (RHE) potential according to the following equation: E RHE = E Ag/Ag/Cl + 0.059 × pH + 0.210.The loading mass was 280 ug cm −2 for all of the catalysts measured.

ORR Measurements
All electrochemical measurements were operated in the N 2 or O 2 saturated 0.1 M KOH electrolyte.The cyclic voltammetry (CV) measurements were recorded at a scan rate of 100 mv s −1 in N 2 or O 2 saturated electrolyte.The LSV measurements were performed at a scan rate of 5 mV s −1 in O 2 saturated electrolyte.Each catalyst repeated at least 3 times of each measurement to exclude possible incidental errors.The numbers (n) of electron transferred per O 2 molecule were calculated according to the following Koutecky-Levich (K-L) Eq. 1 [39]: where j is the measured electrode current density, j k is the kinetic current density, and ω is theelectrode rotating rate.B represents the Levich slope given by the following Eq.2: where n is the number of electrons transferred per oxygen molecule, F is the Faraday constant (F = 96,485 C mol −1 ), D O2 is the diffusion coefficient of O 2 in 0.1 M KOH, where D O2 is 1.9 × 10 -5 cm 2 s −1 , v is the kinetic viscosity (0.01 cm 2 s −1 ), and C O2 is the bulk concentration of O 2 , which is 1.2 × 10 -6 mol cm −3 .Constant 0.2 is used when the rotating speed is expressed in rpm.
The RRDE tests were performed using a Pt ring surrounded by 4 mm diameter GCE.The value was determined by the following Eq.3: The HO 2 − yield is decided by Eq. 4: where i d represents the disk current and i r represents the ring current.N represents the current collection efficiency of the Pt ring which is determined as 0.44.The long-term stability was analyzed by chronoamperometric tests at a fixed potential of 0.7 V vs. RHE and a rotation speed of 1600 rpm in O 2 saturated electrolyte.

Electrochemical Double Layer Capacitances (C dl )
The C dl was measured via a simple CV method. (1) The C dl of various materials can be determined from the CV measurement, which is expected to be linearly proportional to the electrochemical active surface areas.A potential range of −0.1-0.1 V vs. Ag/AgCl was selected for measuring capacitance because no obvious Faradaic currents were observed in this region.The absolute value capacitive currents of ΔJ@0 V/2 were plotted as a function of the CV scan rate of 20, 40, 60, 80, and 100 mV s −1 .The slopes of the fitting data line were the geometric C dl .

Electrochemical
Impedance Spectroscopy (EIS) The EIS measurements were carried out by applying an AC voltage with 5 mV amplitude in a frequency range from 100 to 100 mHz.

Aqueous Zinc-Air Battery (ZAB) Assembly
The air electrode was fabricated via spraying catalyst slurry on a clean carbon cloth with an active area of 1 cm 2 and then dried at room temperature for 24 h.The catalyst slurry was obtained via dispersing 8 mg of FeN 4 -Ti 3 C 2 S x or Pt-C catalysts into 1 ml of mixed solution containing 2-propanol, DI water, and Nafion solution (5 wt%) with a ratio of 10:40:3.The resulting loading mass was 1.0 mg cm −2 .The air electrode served as cathode, 6.0 M KOH with 0.2 M Zn(Ac) 2 addictive as electrolyte, and polished Zn plate electrode as anode to assemble an aqueous rechargeable ZAB.

Stretchable Solid-State Fiber-Shaped ZAB Assembly
The air electrode was assembled via spraying catalyst slurry on a carbon nanotube (CNT) paper with a loading mass of 1.0 mg cm −2 and then dried at room temperature for 24 h.The dual-network PANa-cellulose hydrogel was synthesized using our previous developed method and soaking with 6.0 M KOH with 0.2 M Zn(Ac) 2 addictive as stretchable solid-state electrolyte.The anode was a zinc spring.The stretchable solid-state fiber-shaped ZAB was assembled via the following process: (a) coated the relaxed Zn spring anode with dual-network PANa-cellulose hydrogel electrolyte; (b) stretched the spring-hydrogel system; and (c) coated the catalyst loading CNT paper on the stretched spring-hydrogel system, then released.Due to strong adhesion of hydrogel, the anode and cathode can be firmly adhered to hydrogel electrolyte.The galvanostatic tests were performed via using a Land 2001 A battery test system at room temperature.The charge-discharge polarization and A.C. impedance with 5 mV amplitude in a frequency range from 100 to 100 mHz were determined by using an electrochemical workstation (CHI 760e, Chenhua).

Computational Details
All the first principle calculations were conducted using spin-polarized DFT as implemented in Quantum Espresso [40,41].Generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) functional was chosen to represent the exchange-correlation interaction [42].Grimme's DFT-D3 method was used to determine van der Waals (vdW) interactions [43].A plane-wave cutoff of 55 Ry and a density cutoff of 550 Ry were applied based on standard solid-state pseudopotentials with projector augmented-wave (PAW) method [44,45].The DFT + U calculations were adopted to describe strong on-site Coulomb interaction of localized electrons.The value of Hubbard correction U for 3d orbitals of Ti and Fe was set to 3.0 and 5.0 eV, respectively.A 4 × 4 supercell of MXene monolayer with a sufficiently vacuum slab of 20 Å was built.For sampling the Brillouin zone, Monkhorst-Pack k-point was set as 2 × 2 × 1, and a larger 6 × 6 × 1 k-point was applied to study the electronic properties.All atoms were fully relaxed until the forces on each atom were less than 0.02 eV Å −1 .
The calculation of Gibbs free energy change (ΔG) for each elemental step was based on the computational hydrogen electrode (CHE) model [46], which could be expressed by ΔG = ΔE + ΔE ZPE -TΔS + eU + ΔG pH , where ∆E represents the electronic energy difference between the free standing and adsorption states of reaction intermediates; ∆E ZPE and ∆S represent the changes in zero point energies and entropy, respectively, which are obtained from the vibrational frequency calculations.T represents the temperature and here is set as 298.15K; e and U represent the number of electrons transferred and the electrode applied potential, respectively; ∆G pH represents the free energy correction of pH, which can be derived from: ∆G pH = K B T × pH × ln10.In this work, H 2 and H 2 O were used as the reference states; hence, a series of equivalent reactions for the ORR mechanism are applied to determine ΔG.The complete ORR catalytic process in alkaline condition includes the following five elementary steps [22]: where * indicates the adsorption site and steps ii − v represent the four-electron transfer processes.

Catalysts Fabrication and Structural Characterizations
The FeN 4 -Ti 3 C 2 S x sample was obtained by coating Ti 3 C 2 MXene with Fe salt, 2,2-bipyridine and potassium thiocyanate (KSCN), followed by pyrolysis at 900 °C in N 2 atmosphere and HF acid leaching process.FeN 4 -Ti 3 C 2 sample was obtained via similar method except that KSCN was not used.The SEM image of Ti 3 C 2 MXene, as shown in Fig. 1a, indicates that it is nanosheet structure with wrinkle, while after coating with Fe, N, and S salts before carbonization, it becomes thicker as indicated in Fig. S1.After carbonization, when the Fe, N, S dopant are introduced, the products become thicker and rough and few nanoparticles attached to the surface are clearly seen for FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x in Fig. 1b, c, respectively.The morphology of these two samples, without and with S dopant, is similar as also confirmed by the TEM images indicated in Fig. 1d for FeN 4 -Ti 3 C 2 , and in Fig. 1e for FeN 4 -Ti 3 C 2 S x .The XRD patterns in Fig. 1f show the peaks at 2θ = 6.7°, 17.8°, and 27.6°, which are characteristic of (002), (006), and (008) crystal planes of layered pristine Ti 3 C 2 MXene, respectively [47].Besides, the peaks of other two samples with Fe, N, and S dopants are similar to those of pristine Ti 3 C 2 , indicating that introducing dopant does not change the structure of pristine Ti 3 C 2 ; however, more anatase peaks appear indicating the doping process resulting in mild oxidation of Ti 3 C 2 flakes, which is also consistent with the SEM results that some nanoparticles appear in the rough surface of Ti 3 C 2 after introducing dopant.Besides, the broad peak at ~ 25° suggests the formation of amorphous carbon during the doping process owing to the mid oxidation of Ti 3 C 2 flakes.The Raman spectra in Fig. 1g indicate that the peaks of the samples with Fe, N, and S dopants are similar to those of pristine Ti 3 C 2 ; however, the peak intensity of the carbon skeleton (D/G bands) becomes more obvious in the spectra for samples with dopants compared with pristine Ti 3 C 2 .This is because more defects generated during dopant and carbonization process.Normally, the defects could affect the catalytic performance by influencing the electronic conductivity and active sites of samples [35,48].Figure 1h shows the N 2 adsorption and desorption isotherms of all three samples displaying a type IV isotherm with hysteresis, suggesting the existence of porosities in the structure.
The calculated BET surface areas of these samples follow the trend of Ti 3 C 2 > FeN 4 -Ti 3 C 2 > FeN 4 -Ti 3 C 2 S x , and mesoporous structures are dominant in all three samples as shown in Fig. 1i.The reduced BET surface areas and pore diameter as the increasing of the dopants' species may be caused by the dopants filling effect.The surface compositions of the prepared samples have been analyzed by XPS, as presented in Fig. 2. The survey XPS spectra in Fig. 2a show that the F 1s, O 1s, Ti 2p, and C 1s peaks exhibit in pristine Ti 3 C 2 , while F 1s is absent in samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x , suggesting that during dopants introduced process, the surface F atoms in pristine Ti 3 C 2 disappear.Instead, the N 1s and Fe 2p peaks exhibit in sample FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x , suggesting that the Fe and N atoms are introduced in the pristine Ti 3 C 2 MXene.Besides, we can see that the S 2p peak only exhibits in sample FeN 4 -Ti 3 C 2 S x , indicating that only sample FeN 4 -Ti 3 C 2 S x has S atom as expected.The corresponding elementary composition obtained from XPS survey spectra of sample FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x is provided in Table S1.The high-resolution XPS spectrum of N 1s for sample FeN 4 -Ti 3 C 2 S x shown in Fig. 2b shows a peak at around 396.5 eV that can be assigned to Ti-N species.The Ti-N peaks in high-resolution XPS spectra of N 1s (Fig. S2a) for sample FeN 4 -Ti 3 C 2 are similar to those in sample FeN 4 -Ti 3 C 2 S x .The peak at around 395.4 eV presented in sample FeN 4 -Ti 3 C 2 S x may relate to S coordinated with N atom species, denoted as S-N species here, since this peak is absent in sample FeN 4 -Ti 3 C 2 without sulfur dopant.Besides, the observation of Fe-N x peak centered at 399.1 eV in N 1s of these two samples suggests that the N atoms are coordinated with Fe atoms instead of S atoms coordinating with Fe atoms [22].The high-resolution XPS spectra of Fe 2p for samples FeN 4 -Ti 3 C 2 S x (Fig. 2c) and FeN 4 -Ti 3 C 2 (Fig. S2b) can be fitted to four peaks, with two spin-orbit doublets at around 710.5 and 724.3 eV, which corresponds to the Fe 2p 3/2 and Fe 2p 1/2 , respectively.Besides, the peak at around 714.3 eV is assigned to the Fe-N x configuration, which further confirms that the Fe atoms are coordinated with N atoms in these two samples [1,[49][50][51].The highesolution XPS spectrum of S 2p for sample FeN 4 -Ti 3 C 2 S x as 1 3 shown in Fig. 2d can be fitted to four peaks at binding energies of 162.6, 163.6, 165.2, and 168.8 eV, while the peak at 162.6 eV is attributed to Ti-S bond, indicating that the sulfur atoms substitute the terminal groups in the pristine Ti 3 C 2 MXene and form S-Ti-C bond in sample FeN 4 -Ti 3 C 2 S x [52,53].Since no peak related to S coordinated with Fe can be observed, the exotic S atoms tend to bond to surrounding N atoms instead directly to the Fe atoms.The peaks at 163.8 and 165.2 eV correspond to S 2p 3/2 and S 2p 1/2 of the C-S-C covalent bond of thiophene-S.The peak at 168.8 eV is assigned to oxidized sulfur [54].The high-resolution XPS spectra of Ti 2p for sample FeN 4 -Ti 3 C 2 (Fig. S2c) and sample FeN 4 -Ti 3 C 2 S x (Fig. S2c) can be fitted to four doublets, including two peaks assigned to Ti-C, two assigned to Ti 2+ , two assigned to Ti 3+ , and two assigned to Ti-O, in agreement with previous reports [55].
To further identify the coordinated status of Fe, N, and S atoms in these two samples, X-ray absorption near-edge structure (XANES) and Fourier transform extended X-ray absorption fine structure (EXAFS) spectrometry measurements are conducted.Figure 2e shows the Fourier transforms EXAFS at the Fe K-edge of samples FeN 4 -Ti 3 C 2 , FeN 4 -Ti 3 C 2 S x , commercial FePc (iron phthalocyanine), and Fe foil samples as references, which reveals that samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x exhibit a primary strong peak at ~ 1.5 Å, corresponding to Fe-N peak, originating x .f XRD patterns, and g Raman spectra.h N 2 adsorption-desorption isotherms, and corresponding i pore size distribution of the products from a nitrogen shell surrounding iron atoms in reference to that of FePc [27,56].The Fe-centered coordination number and bond distance of these samples derived from EXAFS fitting is shown in Fig. S3 and Table S2 and confirms the coordination number of FeN x for samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x is ~ 4 which is also similar to that of sample FePc.This confirms that the presence of FeN 4 moieties in samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x and FeN 4 species is usually considered as the main active species for catalysis ORR [57][58][59].The phase-corrected bond length of Fe-N shell is increased from 1.92 to 1.97 Å from sample FePc to FeN 4 -Ti 3 C 2 and then further increased to 1.98 Å for FeN 4 -Ti 3 C 2 S x but shorter than that of Fe-S (2.15 Å), suggesting only Fe-N in these two samples.This also indicates that the square-planar structure of Fe 2+ − N 4 moieties becomes distorted in sample FeN 4 -Ti 3 C 2 compared with that in commercial FePc, which is further distorted in sample FeN 4 -Ti 3 C 2 S x .Besides, the Fe-Fe peaks at ~ 2.2 Å are only presented in sample Fe foil, indicating no iron-based crystalline structures in both samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x [28,60].Figure 2f shows the XANES spectra at the Fe K-edge of samples FeN 4 -Ti 3 C 2 , FeN 4 -Ti 3 C 2 S x , and the commercial FePc, Fe foil samples as references; an enlarged version of this figure, with added y-axis offset, is included as Fig. S4 for clarity.These figures show that sample FePc presents a weak pre-edge peak at ~ 7113.3 eV, assigning to a fingerprint of square-planar Fe 2+ − N 4 moieties [22].This feature becomes less obvious but is still perceptible in sample FeN 4 -Ti 3 C 2 ; however, it is not perceptible in sample FeN 4 -Ti 3 C 2 S x suggesting that the square-planar structure of Fe 2+ − N 4 moieties becomes distorted in sample FeN 4 -Ti 3 C 2 , which is destructed in sample FeN 4 -Ti 3 C 2 S x [26].The distortion of the square-planar structure of Fe 2+ − N 4 moieties is caused by the distortion of the D 4h symmetry as a result of the central Fe away from the N 4 -plane due to the interaction of the Ti 3 C 2 support while incorporation of S terminals further enhances the interaction resulting in destruction of the D 4h symmetry and squareplanar structure [27].

ORR Performance
To investigate the ORR activity of the obtained materials, we conducted cyclic voltammetry (CV) and liner sweep voltammetry (LSV) measurements.The CV curves, recorded in O 2 saturated 0.1 M KOH, are presented in Fig. 3a and show that the oxygen reduction peak potential for pristine Ti 3 C 2 is 0.60 V vs. RHE (its enlarged version in Fig. S5c), indicating that pristine Ti 3 C 2 presents poor ORR activity.After the introduction of Fe and N dopants, the peak potential has been improved to 0.77 V vs. RHE for sample FeN 4 -Ti 3 C 2 .Interestingly, with additional S terminals, the peak potential has been further improved to 0.862 V vs. RHE for sample FeN 4 -Ti 3 C 2 S x , which is comparable to that of commercial Pt-C (0.860 V).However, the pristine Ti 3 C 2 with only S terminal (Ti 3 C 2 S x ) does not show obvious enhancement when compared with pristine Ti 3 C 2 (Fig. S5a).This suggests that the sample FeN 4 -Ti 3 C 2 S x processed the highest ORR activity ORR activity.This phenomenon is also proved in the LSV curves recorded at a rotation speed of 1600 rpm, as indicated in Fig. 3b.The half-wave potential for pristine Ti 3 C 2 is 0.64 V vs. RHE (its enlarged version in Fig. S5d), while for FeN 4 -Ti 3 C 2 it is improved to 0.81 V vs. RHE.A more positive half-wave potential of 0.89 V vs. RHE for sample FeN 4 -Ti 3 C 2 S x is observed, and the improvement is up to 80 mV of halfwave potential after introducing S terminals.Besides, the corresponding Tafel plots as shown in Fig. S5b also indicate that sample FeN 4 -Ti 3 C 2 S x possesses the smallest Tafel slope, suggesting the highest ORR kinetics.However, the half-wave potential and Tafel slope of Ti 3 C 2 S x are almost similar to those of pristine Ti 3 C 2 .We also conduct electrochemical double layer capacitance (C dl ) and electrochemical impedance spectroscopy (EIS) measurements for the above x is comparable to that of commercial Pt-C (0.88 V) and outperforms many other related electrocatalysts reported in recent literatures (Table S3).The rotating disk electrode (RDE) tests recorded at various rotation speeds are further performed to obtain more ORR kinetics information as indicated in Fig. S8.It is obvious that each sample FeN 4 -Ti 3 C 2 (Fig. S8b), FeN 4 -Ti 3 C 2 S x (Fig. S8e), and Pt-C (Fig. S8h) shows well-defined diffusion-limited platforms at rotation speed from 400 to 3600 rpm.The corresponding Koutecky-Levich plots (K-L plots) at various potentials exhibit good linearity.The calculated average electron transfer number (n) is 3.64 for sample FeN 4 -Ti 3 C 2 , 3.97 for sample FeN 4 -Ti 3 C 2 S x , and 3.99 for commercial Pt-C, as shown in Fig. S8c, f, and i, respectively.The electron transfer numbers of sample FeN 4 -Ti 3 C 2 S x and Pt-C are highly close to 4, indicating a four-electron reduction pathway.The parallel rotating ring disk electrode (RRDE) tests shown in Fig. 3c, d show that the electrons transfer numbers of sample FeN 4 -Ti 3 C 2 S x locate between 3.81 and 3.97 with a HO 2 − yield of 2.1-11.9%,which are close to those of commercial Pt-C (n: 3.81-3.99;HO 2 − yield: 0.7-9.1%).In addition, the long-term stability indicated in Fig. S9 indicates that sample FeN 4 -Ti 3 C 2 S x retains 84% of its initial current, which surpasses that of commercial Pt-C with 52% of its initial current retention in a continuous 24 h test, manifesting a better cycling stability for sample FeN 4 -Ti 3 C 2 S x .

Wearable Solid-state ZAB Performance
Since sample FeN 4 -Ti 3 C 2 S x presents high ORR activity, we further investigate the aqueous ZAB performance of this sample and commercial Pt-C-IrO 2 for comparison, as shown in Fig. S10.Notably, the energy efficiency, maximum power density, and cycling stability of FeN 4 -Ti 3 C 2 S x -based ZAB are better than those based on commercial Pt-C-IrO 2 electrodes emphasizing the promising potential of FeN 4 -Ti 3 C 2 S x catalyst for rechargeable ZAB applications.In addition, we also construct a stretchable and wearable fiber-shaped ZAB using our developed FeN 4 -Ti 3 C 2 S x catalyst and a recently developed alkaline tolerant dual-network PANa and cellulose hydrogel (PANa-cellulose) as stretchable solid-state electrolyte [13].
The stretchability of this PANa-cellulose hydrogel soaked with 6 M KOH + 0.2 M Zn(CH 3 COO) 2; Fig. 4a shows that it can be stretched over 1000% strain without any breakage and visible cracking demonstrating excellent stretchability performance.The structure of the fiber-shaped ZAB is depicted in Fig. 4b, using hydrogel electrolyte to coat the Zn spring electrode firstly and then stretch them and finally coat them by using the FeN 4 -Ti 3 C 2 S x loading carbon nanotube paper as air electrode.The charge-discharge profiles and corresponding power density of the fiber-shaped ZAB at initial and 800% stretched states are shown in Fig. 4c,  d, which indicates that the ZAB in stretched state exhibits increased energy efficiency and power density when compared with those in initial state.The increases in the energy efficiency and power density for the stretched state can be attributed to the increased contact areas between the hydrogel electrolyte and the active materials.The maximum power density for ZAB in initial state is 133.6 mW cm −2 and at 800 stretched state is 182.3 mW cm −2 , suggesting that the battery is stretchable and the electrochemical performance is good under stretched state.Besides, the battery shows excellent cycling stability with 110 h stable cycles at 2 mA cm −2 , as shown in Fig. 4e.To demonstrate its wearability, two fiber-shaped ZABs with a length of 10 cm and a diameter of 2 mm have been woven into a wristband and connected to a wearing glove, as shown in Fig. 4f, g, respectively.This wristband can power a set of LEDs on the wearing glove demonstrating the feasibility of such an efficient stretchable and wearable fiber-shaped ZABs based on the prepared FeN 4 -Ti 3 C 2 S x catalyst.

Origin of the ORR Electrocatalytic Activity
The remarkable ORR catalytic activity of our prepared   species and ORR kinetics [30].Besides, we have also performed ESR and M-T measurements to disclose the electron spin configuration of FeN 4 moieties in Ti 3 C 2 support with and without S terminals.The ESR spectra in Fig. S11 indicate an obvious shift and larger g factor for FeN 4 -Ti 3 C 2 S x compared with those of sample FeN 4 -Ti 3 C 2 .The shift of ESR spectra is due to stronger interaction between S terminated Ti 3 C 2 and unpaired electrons in FeN 4 active species for -Ti 3 C 2 S x compared to that of sample FeN 4 -Ti 3 C 2 .The larger g factor for sample FeN 4 -Ti 3 C 2 S x is due to more unpaired electrons occurrence caused by the stronger interaction between FeN 4 species and Ti 3 C 2 S x supporter.The M-T measurements are presented in Fig. S12.The corresponding effective magnetic moment (µ effect ) shown in Fig. 5c

Conclusions
Conventional strategies to design efficient Fe-N-Cs catalysts are based on incorporating FeN 4 species, improving their density, and designing geometric structures for exposing FeN 4 species; however, few works have focused on regulating the electronic structure of FeN 4 species, especially their electronic spin states for improved activity.We demonstrate that introducing sulfur-terminated Ti x catalyst into a wearable ZAB shows a good discharge performance and high cycling stability.This study endows a guideline for regulation on electronic structure of active species via coupling with their support, which is significant to enhance catalytic activity.

Fig. 1
Fig. 1 SEM images for a pristine Ti 3 C 2 , b FeN 4 -Ti 3 C 2 , and c FeN 4 -Ti 3 C 2 S x .TEM images for d FeN 4 -Ti 3 C 2 , and e FeN 4 -Ti 3 C 2 S x .f XRD patterns, and g Raman spectra.h N 2 adsorption-desorption isotherms, and corresponding i pore size distribution of the products

Fig. 2 a
Fig. 2 a XPS survey spectra of pristine Ti 3 C 2 , FeN 4 -Ti 3 C 2 , and FeN 4 -Ti 3 C 2 S x .High-resolution XPS spectra of b N 1s, c Fe 2p, and d S 2p for sample FeN 4 -Ti 3 C 2 S x .e Fourier transforms of Fe K-edge EXAFS spectroscopy oscillations of samples FeN 4 -Ti 3 C 2 , and FeN 4 -Ti 3 C 2 S x , with Fe foil and FePc as references.f Fe K-edge XANES spectra of samples FeN 4 -Ti 3 C 2 , and FeN 4 -Ti 3 C 2 S x , with Fe foil and FePc as references

Fig. 3 a
Fig. 3 a CV curves of pristine Ti 3 C 2 , FeN 4 -Ti 3 C 2 , FeN 4 -Ti 3 C 2 S x , and commercial Pt-C recorded at 100 mV s −1 in O 2 saturated 0.1 M KOH solution.b LSV curves of pristine Ti 3 C 2 , FeN 4 -Ti 3 C 2 , FeN 4 -Ti 3 C 2 S x , and commercial Pt-C at 1600 rpm rotation speeds.c RRDE measurements of samples FeN 4 -Ti 3 C 2 S x , and commercial Pt-C at 1600 rpm.d Electron transfer numbers (n) and HO 2 − yield derived from RRDE tests

FeN 4 - 3
Ti 3 C 2 S x should come from the possible interactions 1 between the FeN 4 species and the sulfur-terminated Ti 3 C 2 support.To verify this interaction, multiple spectroscopies are conducted to explore the electronic structure changes of FeN 4 species supported on pristine Ti 3 C 2 MXene between on sulfur terminated Ti 3 C 2 MXene.UPS is first conducted to study the band structure of samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x .As indicated in Fig. 5a, the cutoff energy (E cutoff ) of FeN 4 -Ti 3 C 2 is 17.1 while FeN 4 -Ti 3 C 2 S x is 17.23.The work function Φ can be calculated according to Φ = hv-|E cutoff -E F | [30], where hv represents the photon energy of the excitation light (21.22 eV) and E F here is 0 eV.Therefore, the work functions Φ are estimated as 4.12 eV for FeN 4 -Ti 3 C 2 and 3.99 eV for FeN 4 -Ti 3 C 2 S x .Furthermore, the valence band maximum (E V ) shifts to lower energy after the introduction of S into the Ti 3 C 2 support as shown in Fig. 5b where E V is estimated to be 2.56 and 2.74 eV for FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x , respectively.The decrease of Φ and E V shift to lower energy for FeN 4 -Ti 3 C 2 S x compared to FeN 4 -Ti 3 C 2 demonstrates that the electrons within FeN 4 moieties become more spatially stable after incorporation of S terminals in Ti 3 C 2 support and that the 3d band center

Fig. 4 aFig. 5
Fig.4a Tensile stress versus strain curve of the prepared PANa-cellulose hydrogel with 6 M KOH + 0.2 M Zn(CH 3 COO) 2 intake; the insets are optical photos of the initial and stretched states of this hydrogel electrolyte.b Schematic illustration of stretchable fiber-shaped ZAB.c Galvanodynamic charge-discharge profiles of fiber-shaped ZAB at initial and 800% stretched states.d Galvanodynamic discharge and corresponding power density curves of fiber-shaped ZAB at initial and 800% stretched states.e Cycling stability test of fiber-shaped ZAB at 2 mA cm −2 ; f Photographs of two fiber-shaped ZABs (length: 10 cm, diameter: 2 mm) woven into a wristband; g Photographs of this wristband connected to a glove.h Photographs of this wristband connected to a glove to power a set of LEDs d x2-y21 after introducing S terminals in sample FeN 4 -Ti 3 C 2 S x , which is consistent with ESR and M-T results.The introduction of S terminals affects the configuration way for Fe coordinating with N, resulting in different configuration of d orbitals for Fe as shown in Fig.5f, g, so their energy orders are different[24,51].The transfer of Fe(II) electron spin configuration from intermediate spin state in FeN 4 moieties for sample FeN 4 -Ti 3 C 2 to high spin state for sample FeN 4 -Ti 3 C 2 S x could enhance the adsorption of molecular oxygen on catalyst surface and thus improve the ORR performance[52][53][54].In addition, the calculated PDOS indicates the electron density of Fe atom reduces and the d band center rises from −2.49 to −2.15 eV.The Fe 3d electron delocalization and d band center upshift can optimize the orbital hybridization of Fe 3d with p orbital of oxygen-containing groups, boosting oxygen-containing groups adsorption on FeN 4 species and ORR kinetics[55][56][57].The free energy diagram of ORR processes for pristine Ti 3 C 2 , FeN 4 -Ti 3 C 2 , and FeN 4 -Ti 3 C 2 S x in alkaline condition (pH = 13) is shown in Fig.5h.It is noted that the formation of *OOH is the rate-determining step (RDS) for ORR of these three systems.After the introduction of S terminals, the adsorption of *OOH intermediate on FeN 4 -Ti 3 C 2 S x significantly increases, indicating that FeN 4 -Ti 3 C 2 S x system exhibits better catalytic activity for ORR, which agrees well with experimental results.

3 C 2
MXene to support FeN 4 species via fabricating FeN 4 -Ti 3 C 2 S x sample can regulate the electronic spin state of FeN 4 species and dramatically enhance catalytic activity toward ORR.Our experimental investigation and theoretical studies uncover that the sulfurterminated MXene induces the central metal Fe(II) in FeN 4 species with original intermediate spin state (d xy the latter the d z2 orbital occupied by a single electron enables their Fe(II) ions to bind oxygen in the end-on adsorption mode favorable to initiate the reduction of oxygen.Furthermore, it induces a remarkable Fe 3d electron delocalization with d band center upshift, optimizing the orbital hybridization of Fe 3d with p orbital of oxygen-containing groups, boosting oxygen-containing groups adsorption on FeN 4 species and ORR kinetics.The resulting FeN 4 -Ti 3 C 2 S x exhibits enhanced 80 mV of half-wave potential compared to that of the FeN 4 -Ti 3 C 2 and also comparable catalytic performance to those of commercial Pt-C.Besides, integrating this FeN 4 -Ti 3 C 2 S

1 3
(XPS) measurements on this catalyst to disclose that the iron is coordinated with nitrogen in the form of FeN 4 , while the sulfur in Ti 3 C 2 S x terminal is bond with N to manipulate the electronic structure of central Fe sites.Furthermore, such FeN 4 coupling withTi 3 C 2 S x leads to a remarkable Fe 3d electron delocalization with d band center upshift and the central metal Fe(II) in FeN 4 species from original intermediate spin state (d xy of Fe (II) changes.This means that the electron density of FeN 4 moieties distributes more thinly at Fe (II) center, while more densely at the nitrogen ligands, resulting in Fe 3d electron delocalization caused by the sulfur terminal with high electronegativity interacting strongly with FeN 4 moieties and stabilizing the valence band maximum by reducing the electron density of Fe(II).It is believed that a decreased electron density of Fe (II) center with strong delocalization can optimize the orbital overlap of Fe 3d with O 2 2p and thus favors the oxygen adsorption on FeN 4 1 3 indicates that the µ effect of sample FeN 4 -Ti 3 C 2 S x is larger than that of sample FeN 4 -Ti 3 C 2 .The larger µ effect demonstrates that the larger number of unpaired d electron of Fe(II) in sample FeN 4 -Ti 3 C 2 S x compared with that of in sample FeN 4 -Ti 3 C 2 .Both ESR and M-T measurements prove the spin state of central Fe(II) in FeN 4 moieties for sample FeN 4 -Ti 3 C 2 is higher than that of FeN 4 -Ti 3 C 2 , and thus, we speculate that the central Fe(II) in FeN 4 moieties for sample FeN 4 -Ti 3 C 2 may possess an intermediate spin state with electron configuration of d xy Ti 3 C 2 S x .To gain deeper insight into the electronic structure of FeN 4 moieties in samples FeN 4 -Ti 3 C 2 and FeN 4 -Ti 3 C 2 S x , DFT calculations are performed.The computational details are presented in Experimental Section.The FeN 4 moieties originally processes the square-planar structure with central Fe(II) in low spin state with electron configuration of d xy Fe ion slightly out of N 4 plane when the FeN 4 moieties anchored on the Ti 3 C 2 support, constrained by the structure of Ti 3 C 2 , as shown in Fig. 5d (left).The additional S atom adjacent to FeN 4 moiety will further pull the central Fe out of N 4 plane, as shown in Fig. 5d (right).The detail top view and side view of computational optimized atomic structure of pristine Ti 3 C 2 , FeN 4 -Ti 3 C 2 , and FeN 4 -Ti 3 C 2 S x are also indicated in Fig. S13.The calculated PDOS, shown in Fig. 5e (above), indicates that the FeN 4 moieties in sample FeN 4 -Ti 3 C 2 exhibit the magnetic moment due to the spin splitting of the partially occupied d xz and d yz orbitals of Fe(II).The d xy and d z xy orbital is also fully filled.The spin splitting of partially occupied d xz and d yz orbitals decrease, but still can be observed, as shown in Fig. 5e.This electron spin configuration in sample FeN 4 -Ti 3 C 2 S x is close to high spin states of central Fe (II) in FeN 4 moiety, as indicated in Fig. 5g.Thus, it suggests that the central Fe(II) in FeN 4 moieties for sample FeN 4 -Ti 3 C 2 processes an intermediate spin state with electron configuration of d xy