Enhanced Redox Electrocatalysis in High-Entropy Perovskite Fluorides by Tailoring d-p Hybridization

Xudong Li, Zhuomin Qiang, Guokang Han, Shuyun Guan, Yang Zhao,Shuaifeng Lou, Yongming Zhu✉

ABSTRACT High-entropy catalysts featuring exceptional properties are, in no doubt, playing an increasingly significant role in aprotic lithium-oxygen batteries.Despite extensive effort devoted to tracing the origin of their unparalleled performance, the relationships between multiple active sites and reaction intermediates are still obscure.Here,enlightened by theoretical screening, we tailor a high-entropy perovskite fluoride (KCoMnNiMgZnF3-HEC) with various active sites to overcome the limitations of conventional catalysts in redox process.The entropy effect modulates the d-band center and d orbital occupancy of active centers, which optimizes the d-p hybridization between catalytic sites and key intermediates, enabling a moderate adsorption of LiO2 and thus reinforcing the reaction kinetics.As a result, the Li-O2 battery with KCoMnNiMgZnF3-HEC catalyst delivers a minimal discharge/charge polarization and long-term cycle stability, preceding majority of traditional catalysts reported.These encouraging results provide inspiring insights into the electron manipulation and d orbital structure optimization for advanced electrocatalyst.

KEYWORDS Lithium-oxygen batteries; KCoMnNiMgZnF3-HEC perovskite fluoride; Entropy effect; Catalytic kinetics; d-p orbital hybridization

1 Introduction

?Rechargeable lithium-oxygen battery (LOB) has drawn intriguing attraction owing to its ultrahigh theoretical energy density (3500 Wh kg-1), which can be expected to become ultimate candidate to replace current state-ofthe-art lithium-ion battery [1].Nevertheless, LOB suffers from low round-trip efficiency and poor cyclability,accompanied with large polarization.Those intractable issues mainly lie in the sluggish oxygen reduction reaction (ORR) and oxygen evolution reaction (OER), closely related to the complex O2/Li2O2redox chemistry in the air cathodes [2].The reversibility of Li2O2deposition/decomposition directly determines the performance of LOB, highlighting the necessity of ameliorating reaction kinetics through constructing well-tailored catalysts.It is widely accepted that catalyst activity depends substantially on the adsorption/desorption of intermediates on surfaceactive sites, which can be modulated via tuning orbital occupancy of active centers [3-5].Various elaborate catalysts withd-block transition metals, as LiO2-affinity regulators, have been rapidly explored in recent years to boosting catalytic conversion of O2/Li2O2in LOBs [6-10].Despite the considerable optimization in catalytic kinetics,the functionality and availability of such catalysts still fall short of expectations, mainly due to the limited number of accessible active sites.

High entropy compounds (HECs), a new frontier in catalysis with broad chemical space tunability, have shown particularly high catalytic activities compared to conventional mono- and bimetallic nanocrystals [11, 12].Infinite elemental combinations enrich the diversities of active sites and local electronic structures, which provides more possibilities for tailoring of functional properties, and thus expands a wide design platform for catalysts with desirable activity, durability, and effectivity.Theoretically, synergistic effect and interaction of dissimilar species can enhance site-to-site electron transfer, allowing simultaneous stabilization of reaction intermediates with moderate binding energies in multi-steps/electron/phase redox conversions of LOB electrochemistry.A majority of works have been dedicated to the exploration of catalytic performance for HECs, and however, the understanding of activity origins and correlations between intrinsic electronic structure and reaction intermediates in HECs is greatly neglected[13-15].Furthermore, the influence of multiple active sites on the nucleation/growth kinetics of Li2O2in redox processes remains ambiguous, hindering the selection and rational design of HECs catalysts for LOBs.

Herein, a sophisticated KCoMnNiMgZnF3, and HECperovskite fluoride, was first introduced into LOBs to overcome the limitations of conventional catalysts in redox process.We systematically investigated the local electronic environments of cation sites and identifiedd-band-dependent preferences for the adsorption of reaction intermediates during catalysis.The density function theory (DFT)calculations revealed that thed-band center anddorbital occupancy of active centers can be regulated by entropy effect, which engenders a favorabled-porbital hybridization between multiple active cations in HEC and LiO2.This moderate interaction grants KCoMnNiMgZnF3-HEC the enhanced catalytic kinetics, effectively lowering the redox energy barrier.Experimental analysis combined with COMSOL multiphysics demonstrated that the multiple active sites in KCoMnNiMgZnF3-HEC moiety served as nucleation seeds synergistically facilitate the homogeneous formation of Li2O2, which in favor of the mass transfer, leading to an impressive electrochemical performance of LOB.The LOB with KCoMnNiMgZnF3-HEC catalyst achieves a very low overall discharge/charge overpotential (0.7 V), extremely high discharge capacity(22,104 mAh g-1at a current density of 200 mA g-1) and outstanding long-term cyclability (over 500 cycles at a current density of 1000 mA g-1with a capacity limitation of 1000 mAh g-1), which is superior to that of the current traditional catalysts.Our work suggests the possibility of manipulating the reaction kinetics and nucleation/growth mechanisms of Li2O2by the unique entropy effect of multiple sites in HEC, providing guidance to the rational design of efficient active sites for LOBs.

2 Experimental Section

2.1 Materials

Cobalt chloride hexahydrate (CoCl2·6H2O, AR),nickel(II) chloride hexahydrate (NiCl2·6H2O, AR), manganese chloride tetrahydrate (MnCl2·4H2O, AR), zinc chloride (ZnCl2, AR), Potassium fluoride dihydrate(KF·2H2O, AR), polyvinyl pyrrolidone (PVP-K29-32, AR) and ethylene glycol (C2H6O2, AR) were supplied from Shanghai Aladdin Biochemical Technology Co., Ltd.(China).The CR2032 cell, glass fiber separator (GF/B), lithium bis(trifluoromethanesulphonyl)imide(LiN(CF3SO2)2), ≥ 99.9 wt%), tetraethylene glycol dimethyl ether (C10H22O5, ≥ 99.5%(GC)), Ketjen Black (KB EC-600JD), lithium foil and carbon paper (H060) were purchased by Aisim (Shenzhen) Technology Co., Ltd.(China).

2.2 Preparation of GO, HCPA and GO/HCPA Nanocomposite Papers

2.2.1 Preparations of the KCoF3, KCoMnNiF3,KCoMnNiMgF3 and KCoMnNiZnF3 Samples

Perovskite fluorides were synthesized accordant with reference as reported, but with some modifications here.Taking the KCoMnNiZnF3-HEC as a representative, 0.8 mmol CoCl2·6H2O, 0.8 mmol NiCl2·6H2O, 0.8 mmol MnCl2·4H2O,0.8 mmol MgCl2, 0.8 mmol ZnCl2, 10 mmol KF·2H2O and 0.4 g PVP-K29-32 were added into in 72 mL ethylene glycol solvent.After dissolving completely, the mixture was treated solvent thermally in a teflon-lined stainless-steel autoclave at 180 °C for 12 h in a homogeneous reactor, and allowed to cool-down naturally.Subsequently, the sample was washed with absolute alcohol, repeated several times until the solution clean up.Finally, the compound was dried in vacuum at 80 °C, obtaining the KCoMnNiZnF3-HEC material.As for other samples, the synthetic process complied with the above procedure, according to their metallic element ratios with the total molar ratio of 4 mmol.

2.2.2 Materials Characterizations

The X-ray diffraction (XRD) was performed on a Rigaku D/max-rb diffractometer equipped with Cu Kα radiation(λ= 1.5406 Å, 40 kV, 20 mA).Transmission electron microscopy (TEM) and high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images were recorded on JEOL JEM2100 microscope operated at 200 kV.Energy dispersive spectroscopy (EDS) was conducted on NORAN system 7 equipped with TEM.PerkinElmer Optima 5300DV ICP-OES system was used to determine the mass fraction of elements.The X-ray photoelectron spectroscopy(XPS) was carried out on Physical Electronics PHI model 5700 instrument with Al Kα radiation.The XAFS spectra were measured on Super Photon ring-8 of the Japan Synchrotron Radiation Research Institute.The electron paramagnetic resonance (EPR) spectroscopy was determined with JEOL JES-X-320 operating at X-band frequency of 9.79 GHz.The Co K-edge XANES data were accounted in transmission mode with Co foil, CoO and Co3O4as references.The obtained EXAFS data was processed based on the standard routine utilizing the Athena software packages.The EXAFS spectra were obtained by deducing the post-edge background from the overall absorption and then normalizing in respect of the edgejump step.Raman spectroscopy (Raman) was were collected on Invia-Reflex of Renishaw with wavelength of 532 nm.The scanning electron microscopy (SEM) was conducted on Helios Nanolab 600i at 20 kV.The in-situ differential electrochemical mass spectrometry was tested at HPR20-EGA.

2.2.3 Electrochemical Measurements

Li-O2 cell assembly and testingthe perovskite fluorides catalysts (30%) and Ketjen Black EC-600JD (KB, 50%) were uniform mixed with polyvinyl pyrrolidone (PVP, 20%) binder in mixture (isopropanol/water = 1:3) to form a homogeneous ink,which was sprayed on a piece of carbon paper (D= 14 mm).After vacuum drying at 80 °C for 12 h, the O2-electrodes were obtained.The electrolyte adopted 1 M lithium bis(trifluoromethane)sulfonimide (LiTFSI) in triethylene glycol dimethyl ether (TEGDME).The anode and separator used lithium foils and Whatman GF/D glass fibers, respectively.The 2032-tyoe coin cells were assembled in an argon-filled glovebox with oxygen and moisture contents below 0.1 ppm.The galvanostatic charge-discharge performances were conducted on a battery test system (LAND, China) at different current densities in the voltage scope of 2.0-4.5 V versus Li/Li+, or in a capacity limitation.The cyclic stability measures were conducted under a limited capacity.All the evaluation tests mentioned above were carried out at room temperature.

2.2.4 Computational Calculations and Simulations

The density functional theory(DFT)calculationsthe DFT were conducted on the Vienna ab initio simulation package (VASP).The models were constructed based on KCoF3(110) (17.7 Å × 8.35 Å × 11.25 Å), and the vacuum depth was set to 20 Å.The 3d/4selectrons of Co, Ni, Mn and Zn atoms, 2s/2pelectrons for F, 2s/2p/3selectrons for Mg, and 4s/3p/4selectrons for K were treated as valence electrons.The generalized gradient approximation (GGA) were used as the gradient correction function (Perdew-Burke-Ernzerhof (PBE).The interaction between valence electrons and core electrons was described by the plane-mode conserved pseudopotential (PAW).The dispersion interactions were described by Grimme’s DFT-D3 methodology.The 2 × 3 × 1 grid points were performed using the Monkhorst-Pack method in the Brillouin region with plane wave truncation energy of 400 eV.In the optimization process, the energy convergence accuracy was set to 1.0 × 10-5eV atom-1, and the structure relaxation was conducted until the maximum force on each atom below 0.02 eV Å-1.

The energy of adsorption for reaction intermediates was calculated as follows:

whereEtis the total energy of the adsorbed system,EaandEbrepresent the total energy of free species and bare surface,respectively.

The phase diagram of the discharge products is built based on the adsorption energy calculation.The SSLOB ORR/OER reaction free energy were determined by Norskov’s calculation method.

Computational simulationsthe simulations were conducted on the COMSOL Multiphysics 6.0 software, and the operational details were shown in support information.

3 Results and Discussion

3.1 Activity Origins of HECs

To reveal the influence of local electronic structure in HECs on the catalytic conversion of O2/Li2O2, a series of unitary KMF3(M = Co, Mn, Ni, Mg, Zn) and quinary KCoMnNiMgZnF3-HEC models were established for DFT calculations (Fig.1a, b).Electron localization functions (Fig.1ch,) and Bader charge (Table S1) elucidate the electron density redistribution of interfacial metal atoms in HEC, which modulates thed-band center, as evidenced by the projected density of states (PDOS, Fig.1i).Figures 1j-s and S1 illustrate the diagrams of charge density difference for the LiO2adsorbed on different sites and the corresponding adsorption energy (ΔEa), in that the blue area and yellow area presents the lost electron and the gained electron, respectively.Clearly,abundant electron transfers occur on metal sites, manifesting the strong electronic interaction between reaction sites and LiO2.It’s worth noting that the Co*, Mn*, Ni*, Mg*, Zn* site in KCoMnNiMgZnF3-HEC achieves more appealing ΔEbof- 2.59, - 3.11, - 1.40, - 2.45, - 2.83 eV, respectively, compared to those in KMF3candidates, enabling a moderate binding strength of LiO2.In fact, the O2/Li2O2conversion in LOB system involves thed-porbital hybridization between active site and key intermediate LiO2.This interaction is related to the occupancy state of antibonding orbitals, which depends on the energy level ofd-band, according to thed-band theory [8,16, 17].Since the antibonding states are always above thedstates in terms of energy, thed-band center model would be a significant descriptor of the interaction between the catalyst and intermediate.Taking Co* adsorbed configuration as a representative,d-band center in HEC shows a slight upshift trend(Fig.1i), implying the less filling of the antibonding states,leading to an elevated coupling of LiO2(Fig.1j, o).Therefore, it is believed that the entropy effect manipulates the local electron configuration and thedorbital electron filling, which endows multiple metal active sites in KCoMnNiMgZnF3-HEC with a favorable adsorbed energy of LiO2.

To profoundly understand the ORR/OER thermodynamics on multiple sites in HEC, we established the phase diagram of possible intermediates (LixOy,x= 1,y= 2;x= 3,y= 4) and discharge products (Fig.1t-x).The adsorption energy calculation is defined as ΔE= ΔE0+neU, wheren,e, andUrepresents the charge transfer number, electron charge, and electrode potential, respectively [18].The formation of LiO2is thermodynamically advantageous on Co*, Mn*, Ni*, Mg*, Zn* site when the discharge voltage is below 2.56, 3.08, 0.83, 2.43, 1.28 V, respectively.Notably, the LiO2can maintain stability on Mn* site at potentials greater than 1.80 V, while Li3O4dominates at potentials greater than 1.92, 1.33, 1.14 V on Co*, Mg*,Zn* site, respectively.These results substantiate the ability of multiple sites in HEC for stabilizing Li2-xO2intermediates.The difference (ΔU) between minimum OER potential and maximum ORR potential (Fig.1y), as a theoretical descriptor of catalytic activity, is predicted on the proposed models based on the Gibbs free energy of each endothermic step (Tables S2, S3).Apparently, Co*, Mn*,Ni*, Mg*, Zn* site in KCoMnNiMgZnF3-HEC shows a lower ΔUcompared with that in KMF3, suggesting the enhanced reaction kinetics.Therefore, the appropriate electronic environment of the KCoMnNiMgZnF3-HEC triggered by entropy effect can realize the multi-active sites for moderate adsorption of key intermediates during the electrocatalysis, which maximizes the utilization of surface electroactivity, dramatically optimizing the reaction pathway.

Fig.1 The optimized structural models of a KMF3 and b KCoMnNiMgZnF3-HEC.Cross-section of the charge density difference diagrams of c KCoF3, d KMnF3, e KNiF3, f KMgF3, g KZnF3, h KCoMnNiMgZnF3-HEC.i PDOS and corresponding calculated d-band centers of Co 3d in different models.Optimized structures and the corresponding binding energy of intermediate LiO2 on the sites of j Co in KCoF3, k Ni in KNiF3,l Ni in KNiF3, m Mg in KMgF3, n Zn in KZnF3, and o Co*, p Mn*, q Ni*, r Mg* and s Zn* in KCoMnNiMgZnF3-HEC.The phase diagram of possible intermediates and discharge products on the sites of t Co*, u Mn*, v Ni*, w Mg* and x Zn* in KCoMnNiMgZnF3-HEC.y The contours of predicted catalytic activity for different sites in various models

3.2 Synthesis and Structural Characterization

Inspired by the highlighted catalytic activity of KCoMnNiMgZnF3-HEC, we would like to shed light on the essentials of synergistic regulation mechanisms in HEC.Based on the DFT prediction, the model electrocatalysts,KCoF3, KCoMnNiF3, KCoMnNiMgF3, KCoMnNiZnF3,KCoMnNiMgZnF3-HEC were experimentally prepared through a one-pot solvothermal route, as illustrated in Fig.2a-d.In XRD patterns with detailed Rietveld refinements (Fig.2e), the KCoMnNiMgZnF3-HEC shares the same characteristics to that of KCoF3, suggesting the cubic perovskite structures.Obviously, the diffraction peak originating from the (110) plane slightly shift to lower angles after changing in the numbers of elements, demonstrating a strong high-entropy effect.The SEM and TEM images(Figs.2f, i and S2) exhibit a uniform cubic shape with an average size of 20 nm.The lattice fringe spacing of 0.288 and 0.292 nm can be assigned to (110) planes of KCoF3(Figs.2g and S3a) and KCoMnNiMgZnF3-HEC (Figs.2j and S3b), respectively, in good consistent with the XRD results.The inductively coupled plasma-optical emission spectrometer (ICP-OES, Table S4) determines the chemical formula of K0.95Co0.24Mn0.20Ni0.22Mg0.18Zn0.16F3.05.The corresponding fast Fourier transform (FFT) patterns (Fig.2h,k) further confirm the single crystal structures of the products.Energy dispersive X-ray spectroscopy (EDS) mapping images (Fig.2l, m) present the highly homogeneous elemental distribution, validating the complete formation of singlephase KCoMnNiMgZnF3-HEC.

The XPS was carried out to investigate the chemical states and electronic effects of KCoMnNiMgZnF3-HEC.Figure S4 demonstrates the presence of these Co, Mn, Ni, Mg,and Zn elements in KCoMnNiMgZnF3-HEC.From the Co 2pspectrum (Figs.3a and S5), the main peaks at 780.0,796.4, 782.2, 798.6, 786.3 and 802.7 eV correspond to Co3+2p3/2, Co3+2p1/2, Co2+2p3/2, Co2+2p1/2, Co Sat.2p3/2,Co Sat.2p1/2, respectively [19].Compared with KCoF3,Co* in KCoMnNiF3, KCoMnNiMgF3, KCoMnNiZnF3,KCoMnNiMgZnF3-HEC show negative shifts, indicating the electron donating nature.In comparison, the electron interaction in HEC is most moderate as evidenced by a very minor variation in binding energy.To determine the spin polarization regulation of the as-prepared KCoMnNiMgZnF3-HEC,the electron paramagnetic resonance (EPR) spectroscopy was carried out.The major characteristic signals with ag-factor of 2.037 (Eq.S1) correspond to the spin Co, Mn,Ni, Mg, Zn species were detected (Figs.3b and S6), manifesting the presence of unpaired electrons.Meanwhile, the X-ray absorption spectroscopy (XAS) was performed to understand the critical role of high-entropy effect in manipulating the electronic structure of KCoMnNiMgZnF3-HEC.X-ray absorption near edge structure (XANES) spectra of Co K-edge (Fig.3c) reveal a positive oxidation state of Co2+in KCoF3and KCoMnNiMgZnF3-HEC.In contrast with KCoF3, the peak of derivative XANES of KCoMnNiMgZnF3-HEC (Fig.3d) was shifted to high energy, implying the higher valence state of Co, consistent with the XPS analysis results.In the Fourier-transformed EXAFS spectra (Fig.3e), KCoMnNiMgZnF3-HEC presents only one main peak at ca.1.65 Å, corresponding to the combined Co-F coordination [20].This coordination path can be better visualized by the maxima in the wavelet transform(WT) of the EXAFS (Fig.3f), which is quite different compared to that of Co foil, CoO and Co3O4.We therefore conclude that the high-entropy effect can modulate the chemical environment and local structure, facilitating the 3delectron redistribution and further optimizesd-porbital hybridization during electrocatalysis process.

3.3 Evaluation of Electrochemical Performance

To make comprehensive analysis of the structure-performance relationship for multi-sites electrocatalyst, the electrocatalytic properties of the as-prepared electrocatalysts were evaluated in a LOB system (Fig.4a).Figure 4b depicts the first discharge-charge voltage curves at a current density of 500 mA g-1under an upperlimit capacity of 1000 mAh g-1.It is evident that the KCoMnNiMgZnF3-HEC cathode displays the smallest discharge-charge polarization (0.24/0.46 V), far lower than that of other candidates, illustrating the synergetic action among multiple sites on reaction kinetics of O2/Li2O2conversion.The corresponding details of ORR process can be visualized in Fig.4c.In deep discharge-charge process(Fig.4d), the KCoMnNiMgZnF3-HEC cathode delivers an exceptionally high discharge specific capacity of 22,104 mAh g-1at current density of 200 mA g-1, much superior to that of KCoF3(6710 mAh g-1), KCoMnNiF3(7574 mAh g-1), KCoMnNiMgF3(12,423 mAh g-1) and KCoMnNiZnF3(17,482 mAh g-1).Significantly, the KCoMnNiMgZnF3-HEC cathode can maintain its higher capacity retentions and coulombic efficiencies than the other counterparts even at a current density of 400, 800 and 1000 mA g-1(Figs.4e and S7-S11, Table S5), which corroborates prominent advantage of KCoMnNiMgZnF3-HEC in boosting electrochemical performance.Electrochemical impedance spectra (EIS) were performed to survey the reversibility of LOB.As depicted in Fig.S12 (Table S6),the KCoMnNiMgZnF3-HEC-based LOB exhibits a low charge transfer resistance (Rct) at full discharging, indicating a rapid charge transfer process at the three-phase interface of cathode-electrolyte-Li2O2.Upon recharging, the KCoMnNiMgZnF3-HEC-based LOB almost fully recovers to its initial state, which stands in stark contrast to that of KCoF3and highlights its exceptional reversibility.

Fig.3 a High-resolution Co 2p XPS spectra of KCoF3, KCoMnNiF3, KCoMnNiMgF3, KCoMnNiZnF3 and KCoMnNiMgZnF3-HEC.b EPR spectra of KCoF3 and KCoMnNiMgZnF3-HEC, c XANES spectra, d first derivatives of XANES and e Fourier transformed EXAFS for Co foil,Co3O4, KCoF3, KCoMnNiMgZnF3-HEC, and f the corresponding WT contour plots

Fig.4 a Schematic of a LOB and ORR/OER reaction mechanism.b The first discharge-charge voltage curves at a current density of 500 mA g-1 under an upper-limit capacity of 1000 mAh g-1, and c the corresponding details of the selected region in (b), d the first galvanostatic discharge/charge curves at a current density of 200 mA g-1 between 2.0 and 4.5 V, e discharge capacities at various current densities,f the cycle stability and g-k the corresponding galvanostatic discharge/charge curves extracted from whole cycling of the Li-O2 cells with KCoF3, KCoMnNiF3, KCoMnNiMgF3, KCoMnNiZnF3, and KCoMnNiMgZnF3-HEC.l Performance comparisons (capacity and cycle stability)of KCoMnNiMgZnF3-HEC and other representative published efforts

The cycle life as another critical parameter is measured by long-term discharge/charge test.As revealed in Fig.4f, the KCoMnNiMgZnF3-HEC based LOB achieves a long cycle lifetime of 1000 h (500 cycles) at a high current density of 1000 mA g-1with a limited capacity of 1000 mAh g-1, 1.3,2.2, 2.6, 8.3 times higher than that of KCoF3, KCoMnNiF3,KCoMnNiMgF3, and KCoMnNiZnF3, respectively.Figure 4g-k presents the corresponding discharge-charge curves of those candidates during long cycling.It is apparent that the KCoMnNiMgZnF3-HEC cathode remains a relatively lower overpotential over long periods, demonstrating extraordinary reversibility.As a result, the notorious parasitic reactions from electrolyte degradation and catalyst poisoning can be effectively suppressed, which contributes to the enhancement of round-trip efficiency and reversibility for LOB.Combining previous theoretical simulations (Fig.1)and experimental determinations (Fig.3), the excellent electrochemical performance of KCoMnNiMgZnF3-HEC can be attributed to the synergetic action among multiple sites, in which the entropy effect engenders an increased population of unpaired electrons, leading to the electron rearrangement with an elevated spin polarization, dramatically optimizing thed-porbital hybridization between cation sites and key intermediate.To our knowledge, those performance parameters of the LOB with KCoMnNiMgZnF3-HEC precede majority of traditional catalysts reported in previous studies(Fig.4l and Table S7) [3, 4, 21-27].

3.4 Evolution of Cathodes Products during Discharge and Charge Process

Inspired by the impressive electrochemical performance of KCoMnNiMgZnF3-HEC, the morphological evolution and chemical characteristics of cathodes products at various discharge stages was further investigated.Through observation on the SEM images (Fig.5a-e), it is found that the grainshaped products appear on those cathodes after discharged to 1000 mAh g-1.As revealed in the date visualization(Fig.5f-j), the distribution of products becomes more and more uniform as the number of metallic elements increment in catalyst, suggesting that the multiple active sites in HEC are capable of promoting the homogeneous nucleation of products.After full discharging, the products on the KCoF3(Fig.5k, p) and KCoMnNiF3(Fig.5l, q) cathodes evolve into a closely packed amorphous morphology, which can be identified as Li2O2and Li2CO3by high-resolution XPS analysis(Fig.5u, v).In stark contrast, toroidal Li2O2accommodations are deposited evenly on the KCoMnNiMgF3(Fig.5m,r), KCoMnNiZnF3(Fig.5n, s) and KCoMnNiMgZnF3-HEC(Fig.5o, t) cathode with the absence of obvious signal of Li2CO3(Fig.5w-y).It is worth mentioning that the order degree of cathode products correlated strongly with the number of atomic species in catalyst, verifying the modulating effect of the exposed multiple active sites in HEC on the growth pathway of Li2O2(Schematic diagram in Fig.5z).Essentially speaking, the Li2O2morphology and structure are the most important indexes in dominating LOB performance, including overpotential, capacity, cyclability,etc., while the Li2O2formation mechanism is influenced by the adsorption behavior of key intermediates, which is closely relevant to thed-porbital hybridization between catalyst and LiO2.In combination with DFT calculations(Fig.1i-s), it is reasonable to deduce that mixed elements in KCoMnNiMgZnF3-HEC engender significant 3dcharge redistribution and create multiple active sites with optimized energy barriers for stabilizing LiO2intermediates, which dramatically improves the electrochemical performance of LOB.

To thoroughly elucidate the relationship between the Li2O2growth model and the electrochemistry toward ORR kinetics, the Li2O2deposition behavior on O2-electrode was simulated using COMSOL Multiphysics (Fig.S13,Table S8) [28].Based on Kolmogrov’s phase transformation theory [29] and Sampson and Lynden models [30] combined with experimental results, two 2D transient matrixes were constructed (Figs.S14 and S15).Figure 6a displays the contours of O2concentration distributions and Li2O2porosity on O2-electrode at different discharge stages, where the heterogeneous nucleation of Li2O2on KCoF3cathode donates as model I and homogeneous nucleation of Li2O2KCoMnNiMgZnF3-HEC cathode donates as model II, taking the previously obtained SEM characterizations (Figs.5a,k, e, o) as evidence.It is clear that the O2transfer kinetics slows down with the increase of depth of discharge (DOD),which is attributed to the clogging of porous O2-electrode due to the accumulation of Li2O2.We therefore speculate that the formation of O2“dead zone” in cathode was the most common cause of LOB failure.In fact, the accumulation of discharge products on O2-electrode cuts offthe O2-electrode-eletrolyte interaction, which greatly discounts the efficiency of electron and mass transfer, resulting in a constantly-increasing polarization.Videos S1 and S2 exhibited the animate simulations of model I and model II,respectively.From the comparison, the O2diffusion depth,O2concentration (Fig.6b, c), cathode porosity (Fig.6d, e)and Li2O2porosity (Fig.6f, g) on model II were significantly higher than that of model I, suggesting that the homogeneous deposition of Li2O2triggered by multiple active sites in KCoMnNiMgZnF3-HEC can enhance the mass transfer at the three-phase interface among gas-Li2O2-electrolyte, and thereby alleviating the passivation of O2-electrode.Under such circumstance, the LOB enables a fast reaction kinetics,which holds attractive electrochemical performance, well consistent with simulated results (Fig.S16).

Fig.5 SEM images of cathodes after discharged to 1000 mA h g-1 at a current density of 500 mA g.-1 with a KCoF3, b KCoMnNiF3, c KCoMnNiMgF3, d KCoMnNiZnF3, e KCoMnNiMgZnF3-HEC, and f-j The corresponding date visualization.SEM images of cathodes after full discharge with k KCoF3, l KCoMnNiF3, m KCoMnNiMgF3, n KCoMnNiZnF3, o KCoMnNiMgZnF3-HEC, and p-t The corresponding date visualization.High-resolution Li 1s XPS spectra of cathodes after full discharge with u KCoF3, v KCoMnNiF3, w KCoMnNiMgF3, x KCoMnNiZnF3, y KCoMnNiMgZnF3-HEC.z Schematic illustration of discharge reaction mechanisms on KCoF3 (lift) and KCoMnNiMgZnF3-HEC (right)

Fig.6 COMSOL simulations of Li2O2 growth kinetics on cathodes with KCoF3 (model I) and KCoMnNiMgZnF3-HEC (model II) for a O2 concentration distributions and Li2O2 porosity at different discharge stages (0% DOD, 25% DOD, 75% DOD, 100% DOD), b and c the corresponding curves of O2 concentration distributions, d and e cathode porosity and f and g Li2O2 porosity

3.5 Function of Entropy Effect on Enhanced Cyclic Performance

Further investigations were centered on the reversibility of those kinds of cathode.To identify the cathode surface condition, the cycled LOBs were disassembled and analyzed.XRD (Fig.7a) and Raman patterns (Fig.7b) reveal the formation of Li2CO3, LiOH, LiAc, and LiRCO3as main side-products on KCoF3, KCoMnNiF3, KCoMnNiMgF3and KCoMnNiZnF3cathode [31, 32].In striking contrast,no signal peak assigned to those byproducts emerges on KCoMnNiMgZnF3-HEC cathode, even after 500th cycles,unveiling the suppression effect of multiple active sites in HEC on side reactions.It is clear from Fig.7c, d that accumulative products completely cover the KCoF3cathode after cycling.Since the attachment of those indissoluble products on cathode surface during cycling, catalytic active centers were gradually blocked, leading to a restricted electrochemistry and aggravated side reactions.In comparison,the KCoMnNiMgZnF3-HEC cathode (Fig.7e, f) basically maintains a fresh interface condition throughout the entire cycles, demonstrating the excellent anti-passivation power.Moreover, the cycled KCoMnNiMgZnF3-HEC exhibits a superior preservation of its pristine structure compared to KCoF3, as evidenced by Fig.S17, suggesting the excellent stability.Those results provide obvious proof of significant improvement on the long-term durability of LOB with KCoMnNiMgZnF3-HEC catalyst.Accordingly, a schematic illustration of the cycling electrochemistry occurring on the KCoF3and KCoMnNiMgZnF3-HEC cathode is demonstrated in Fig.7g.

To determine the reaction pathway and conversion kinetics, high-resolution ex-situ XPS techniques were conducted on the KCoMnNiMgZnF3-HEC cathode at selected pivotal states (Fig.7h).In the Li 1sspectra (Fig.7i), the peaks around 55.9 and 54.8 eV can be assigned to Li2-xO2and Li2O2, respectively.At the beginning of discharge (stage I), the signals of Li2-xO2and Li2O2appeared in the Li 1sspectrum, certifying a multi-step process of ORR.As the discharge proceeded, the ratio of Li2-xO2/ Li2O2decreased significantly (stage II), indicating the transformation of Li2-xO2intermediate into Li2O2product.Therefore, the possible ORR process on KCoMnNiMgZnF3-HEC cathode can be described by the following Eqs.(2-4):

Fig.7 a XRD, and b Raman spectra of cathodes with KCoF3, KCoMnNiF3, KCoMnNiMgF3, KCoMnNiZnF3, KCoMnNiMgZnF3-HEC after cycling at a current density of 1000 mA g-1 under an upper-limit capacity of 1000 mAh g-1.c SEM images and d EDS mapping of cathodes after cycling with KCoF3.e SEM images and f EDS mapping of cathodes after cycling with KCoMnNiMgZnF3-HEC.g Schematic illustration of cycling electrochemistry on KCoF3 (lift) and KCoMnNiMgZnF3-HEC (right).h Galvanostatic discharge-charge curves of the KCoMnNiMgZnF3-HEC electrode at a current of 500 mA g-1 with an upper-limit capacity of 1000 mAh g-1, and i high-resolution XPS spectra of Li 1s at selected pivotal states as indexed in (h).In-situ DEMS curves of LOB with KCoMnNiMgZnF3-HEC during charging process for j 1st cycling, k 100th cycling.l Schematic illustration of ORR/OER reaction mechanisms on KCoMnNiMgZnF3-HEC cathode

When recharged to state III, the intermediate Li2-xO2reappears on KCoMnNiMgZnF3-HEC cathode, implying the delithium process of Li2O2.At the final stage of recharging, the absence of Li 1speak presents the complete decomposition of Li2O2.The results verify the ability to stabilize key intermediates on multiple active sites in KCoMnNiMgZnF3-HEC cathode.

Moreover, the underlying reaction kinetics was analyzed with the assistance of in-situ differential electrochemical mass spectrometry (DEMS).The ratio of 2e-/O2on KCoMnNiMgZnF3-HEC cathode was displayed in Fig.S18(Eq.S2), supporting the discharging mechanism involving the reduction of O2along the two-electron pathway(O2+ 2Li++ 2e-→ Li2O2) [33].Notably, during charging process, the oxidation processes were predominantly governed by oxygen releasement, as evidenced by the ratio of 2e-/O2in Fig.7j, suggesting the electrochemical decomposition of Li2O2( Li2O2*-2e-→2 Li++O2), without involvement of other side reactions, aligning well with the XRD (Fig.7a) and Raman (Fig.7b) results.Due to the stabilization of intermediates on KCoMnNiMgZnF3-HEC cathode, two-electron transfer can be achieved most effectively during ORR/OER process, which is of vital importance in boosting reaction kinetics.Encouragingly, the gas analysis plots of O2evolution maintained well in the 100th cycles(Fig.7k), where the ratio ofv(e-):v(O2) is determined as 2.02:1, manifesting considerable reversibility of Li2O2formation/decomposition on KCoMnNiMgZnF3-HEC cathode during long-term cycling.We speculate that the improvement of the durability of LOB is closely correlated to the modulation effect of multiple active sites on the reaction pathway, as illustrated in Fig.7l.

3.6 Modulation Mechanism of Entropy Effect on Catalytic Kinetics

To in-depth shed light on the modulation mechanism of entropy effect on catalytic kinetics, we further carried out the DFT simulation and analysis.The local electron configurations of Co* site in KCoF3, KCoMnNiF3, KCoMnNiMgF3KCoMnNiZnF3and KCoMnNiMgZnF3-HEC were first examined through the PDOS calculations.The Co* 3dPDOS (Fig.8a) reveals the rearrangement of 3dorbital(dxz/xy/yz,dx2-y2,dz2) electron in both spin channels after introducing more metallic elements, signifying the strong interactions among elements.Especially, thedx2-y2anddz2in KCoMnNiMgZnF3-HEC leap into Fermi energy level compared to that of KCoF3, demonstrating the electron transition from low-energy orbitals to high-energy orbitals [34-36].Combined with the calculation ofg-factor(Fig.3b), it can be concluded that the local charge redistribution increases the unpaired electron density indxz/xy/yzorbitals, enabling the alteration in the spin polarization.The difference charge density maps (the insets of Fig.8a) present more electrons transfer of Co* in KCoMnNiMgZnF3than that of KCoF3, which points to a higher positive valence of Co* in HEC, well consistent with XANES results (Fig.3c,d).Based on the calculated electron spin states (Tables S9-S13), the possible electron arrangements of Co* in all those candidates were identified, as illustrated in Fig.8b.Considering the moderate spin polarization of active sites in KCoMnNiMgZnF3-HEC, it is believed that the adsorption of key intermediates in O2/Li2O2conversion can be greatly optimized (Figs.1j, o and S19), thus bringing about an enhanced intrinsic activity.Base on hybrid orbital theory, the orbital interactions between Co* centers and LiO2intermediate were further discussed.As revealed by Fig.8c,thedxz/xy/yzorbitals in KCoMnNiMgZnF3-HEC match better with the 2porbital of adsorbed LiO2in contrast with that of other counterparts, unveiling the optimald-porbital hybridization between HEC and intermediate.Due to the 3dorbital coupling of Co* in KCoMnNiMgZnF3-HEC, the electron delocalization decreased thedxz/xy/yzorbitals energy level,which enables easier electron injection from the Co*dxz/xy/yzorbitals into O 2porbital of LiO2, thereby facilitating the conversion kinetics.On this basis, thed-porbital hybridization during ORR process is visualized in Fig.8d.In detail,the antibonding orbitals induced by the coupling between O 2pand Co* 3dlocated below the Fermi level.Owing to the entropy effect in KCoMnNiMgZnF3-HEC, the Co*site denotes a higherd-band center and the electron donating nature, which engenders a favorable binding interaction between key intermediate LiO2and Co* site, reinforcing the intrinsic reaction activity.

Figures S19a-d and 8e picturizes the reaction pathways toward the formation of (Li2O2)2clusters at different overpotentials for the LOBs on Co* site in those various catalysts built on the calculated free energy (Table S14).The discharge and charge overpotentials were referred to evaluate the catalytic activity, defined as ΔUdc=Udc-Ueq,ΔUc=Uc-Ueq, whereUdc,UeqandUcrepresents the highest discharge potential, equilibrium potential and the lowest charge potential, respectively.It is evident from Fig.8f-g that the overpotential of KCoMnNiMgZnF3-HEC(ΔUdc/ΔUc= - 0.57/0.40 V) is relatively lower than that of KCoF3(ΔUdc/ΔUc= - 1.33/1.21 V), KCoMnNiF3(ΔUdc/ΔUc= - 0.70/1.14 V), KCoMnNiMgF3(ΔUdc/ΔUc= - 0.99/0.61 V) and KCoMnNiZnF3(ΔUdc/ΔUc= - 0.60/0.66 V), marking the highest catalytic activity (Fig.8h).Those DFT results elucidate the modulating effect of the exposed multiple active sites in KCoMnNiMgZnF3-HEC on enhancing ORR/OER kinetics, which strongly corroborated the obtained experimental determinations.

4 Conclusions

In summary, guided by the DFT screening, we tailored an ideal KCoMnNiMgZnF3-HEC catalyst to support the ORR and OER reactions for the development of high-performance LOBs.It is revealed that the entropy effect of multiple sites in KCoMnNiMgZnF3-HEC triggers the appropriate regulation of 3dorbital structure, leading to a moderate hybridization with theporbital of key intermediates.As a result,KCoMnNiMgZnF3-HEC-based LOB exerts an ultrahigh discharge capacity and outstanding long-term cyclability with a lowered overpotentials.Accordingly, we may draw a logical conclusion that, the KCoMnNiMgZnF3-HEC catalyst is, momentous to positively manipulate discharge/charge behavior of LOBs.This study demonstrates that thedorbital occupancy unidirectionally affectsd-porbital hybridization between intermediates and catalysts and determines the catalytic kinetics.Tuning the active centers with desirable electron distribution is able to improve the electrochemical performance of LOBs.These findings provide an in-depth understanding of the correlation between multiple active sites in HEC and reaction intermediates, bringing new horizon to the electron modulation andd-band center optimization for electrocatalyst.

Although our proposed high-entropy perovskite fluoride catalyst has shown promising results in LOB catalysis, further research is required to fully comprehend the relationship between structural engineering and the optimization of catalytic behaviors, specifically in identifying site function.Furthermore, there is a need to bridge the gap between experimental design and theoretical prediction to achieve seamless connections in the catalytic applications of HEC catalysts.From the long-term point of view, advanced synthesis strategy, high-throughput screening technology and nondestructive characterization are essential for the functionalization of HEC catalysts, both for catalytic purposes and beyond.Through these collaborative efforts, we believe that the materials science and catalysis communities can come together to explore new frontiers in the discovery and design of high-efficient HEC catalysts.

AcknowledgementsP.G.acknowledges the financial support from the Youth Foundation of Shandong Natural Science Foundation (No.ZR2023OB230), National Natural Science Foundation(No.22309035) and Double First-class Discipline Construction Fund Project of Harbin Institute of Technology at Weihai (No.2023SYLHY11).

Declarations

Conflict of interestThe authors declare no interest conflict.They have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

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

Supplementary InformationThe online version contains supplementary material available at https:// doi.org/ 10.1007/s40820- 023- 01275-3.