Tuning of lattice oxygen reactivity and scaling relation to construct better oxygen evolution electrocatalyst

Developing efficient and low-cost electrocatalysts for oxygen evolution reaction is crucial in realizing practical energy systems for sustainable fuel production and energy storage from renewable energy sources. However, the inherent linear scaling relation for most catalytic materials imposes a theoretical overpotential ceiling, limiting the development of efficient electrocatalysts. Herein, using modeled NaxMn3O7 materials, we report an effective strategy to construct better oxygen evolution electrocatalyst through tuning both lattice oxygen reactivity and scaling relation via alkali metal ion mediation. Specifically, the number of Na+ is linked with lattice oxygen reactivity, which is determined by the number of oxygen hole in oxygen lone-pair states formed by native Mn vacancies, governing the barrier symmetry between O–H bond cleavage and O–O bond formation. On the other hand, the presence of Na+ could have specific noncovalent interaction with pendant oxygen in *OOH to overcome the limitation from linear scaling relation, reducing the overpotential ceiling. Combining in situ spectroscopy-based characterization with first-principles calculations, we demonstrate that an intermediate level of Na+ mediation (NaMn3O7) exhibits the optimum oxygen evolution activity. This work provides a new rational recipe to develop highly efficient catalyst towards water oxidation or other oxidative reactions through tuning lattice oxygen reactivity and scaling relation.

T he oxygen evolution reaction (OER) is a key reaction and constitutes the bottleneck in many energy conversion and storage systems such as water electrolyzers, rechargeable metal-air batteries and regenerative fuel cells [1][2][3] , due to its intrinsically sluggish kinetics 2,4-6 . Considering the origin of O 2 product, there are two widely accepted OER mechanisms including adsorbate evolution mechanism and lattice oxygen oxidation mechanism 5,7 . Regardless of which OER mechanism is applicable on a catalyst surface, it has been reported that O-O bond formation can follow two different pathways, i.e., acid-base nucleophilic attack and O-O direct coupling 7,8 . For the former, there is an inherent linear scaling relation (LSR) between the adsorption energy of *OOH and *OH intermediates. One implication of the above LSR is that the key steps of O-H bond cleavage and *OOH formation are mutually competing, rendering a minimum theoretical overpotential of~0.4 eV even for the best possible material 5,6 . For the latter, it is not subject to such LSR constraint, but the specific catalytic structural motif to trigger O-O direct coupling is difficult to realize for most catalytic materials [7][8][9] . Therefore, current research efforts are mainly directed to optimize the pathway of acid-base nucleophilic attack and overcome the limitation from such LSR for developing practical electrocatalysts.
Activating lattice oxygen to generate spin-characteristic ligand holes can tune the lattice oxygen reactivity that links to energy barrier symmetry between O-H bond cleavage and *OOH formation. Further reducing overpotential ceiling requires selective stabilization of *OOH over *OH to overcome such LSR 5,10,11 . Prior studies have indicated that the introduction of hydrated alkali metal ions (A + ) as promoter can stabilize the key intermediates or transition states via the noncovalent interaction [12][13][14] . Inspired by this, A x Mn 3 O 7 (0 < x ≤ 2) materials, as one type of alkali metal-incorporated metal oxides, can be a good platform for unveiling how to rationally design better OER electrocatalysts through tuning lattice oxygen reactivity and scaling relation mediated by alkali metal ion. On one hand, the native Mn vacancies in MnO 2 layers generate oxygen lone-pair states (|O 2p ), which provides a necessary condition for activating lattice oxygen in view of structural stability 15,16 . On the other hand, the alkali metal ions are directly incorporated into the MnO 2 interlayers, offering the noncovalent interaction between alkali metal and *OOH, and a reduced theoretical overpotential ceiling can be expected. Additionally, the Mn migration would be suppressed if the ionic radii of A and Mn differ largely 17 . In this regard, Na + (102 pm) shows larger contrast in ionic radii with high-valence Mn 4+ (54 pm) and its use in the compound looks promising.
Herein, we use Na x Mn 3 O 7 with tunable number of Na + as model to unlock the specific coordination configuration that can regulate the barrier symmetry between O-H bond cleavage and *OOH formation on the basis of overcoming the LSR between *OOH and *OH. Combining theoretical and experimental approaches, we reveal that the number of Na + is critical to the overall activity improvement. In terms of electronic effect, the O-O bond formation is promoted as the number of Na + reduces, because of the increased number of oxygen holes in |O 2p upon activating lattice oxygen. Correspondingly, the relative barrier between O-H bond cleavage and O-O bond formation is regulated. Contrarily, in terms of geometric effect, the overpotential ceiling increases as the number of Na + reduces, because of the weakening of Na + -specific stabilizing effect on pendant oxygen in *OOH. As a result of the above two opposite effects, an intermediate level of Na + mediation, in this case, NaMn 3 O 7 , exhibits the optimum OER activity. This work provides a guideline for the development of better catalysts towards water oxidation or other oxidative reactions through tuning both lattice oxygen reactivity and scaling relation.

Result
Pathway competition for O-O bond formation. Using density functional theory (DFT) calculations, we firstly probe the regions of space where oxygen lone-pair states locate and then unraveled how spin-characteristic ligand holes is generated upon activating lattice oxygen ( Fig. 1 and Supplementary Figs. 1-3). For modeled Na x Mn 3 O 7 (Na 2x/7 (Mn 6/7 □ 1/7 )O 2 , □ represents Mn vacancy, x = 2, 1.5, 1, and 0.5) slabs, there are two kinds of oxygen ions ( Fig. 1a and Supplementary Fig. 3), in which O1 is coordinated with three Mn ions and O2 is coordinated with two Mn ions, respectively. As such, one of the O(2p) orbitals pointing toward Mn vacancy in O2 coordination environment is non-bonded. According to the amplitude of charge transfer energy (Δ) and d-d Coulomb interaction (U) 18 , Na 2 Mn 3 O 7 (U > Δ) is located at charge-transfer regime, showing an empty metallic band lying above the fully filled |O 2p band (Fig. 1a), which is evidenced by the projected DOS of the Mn(3d) and O(2p) orbitals (shaded region around E F with the dominant oxygen character in Fig. 1b). As confirmed from the projected density of states and partial charge density near E F (Fig. 1b, c and Supplementary Table 1), more oxygen holes from |O 2p are generated upon activating lattice oxygen as the number of Na + reduces. Moreover, such generated oxygen holes are stabilized through π-type interaction between the occupied O 2p and the occupied Mn-t 2g in Mn oct O 6 where Mn would be in their maximum achievable oxidation states of 4+ 19,20 . More specifically, we investigated six representative O2 coordination environments (from S1 to S6 based on the number in neighboring Na + , Fig. 1c) for Na x Mn 3 O 7 to quantify the lattice oxygen reactivity using the number of oxygen holes (h O ). As seen, the magnetization moments of these oxygen ions increased from 0.11 μ B (S1) to 0.60 μ B (S6), which is further confirmed by Bader charge analysis as the charges of such oxygen ions increase from −1.03 e (S1) to −0.66 e (S6). More details about calculating the number of oxygen hole via crystal field theory, magnetization moment, and bader charge can be found in Supplementary Fig. 1.
Considering the dual possible roles of oxygen holes in O-O bond formation via either acid-base nucleophilic attack or O-O direct coupling, we further unravel the pathway competition from the perspectives of both thermodynamics and kinetics. In detail, we systematically correlate both pathways with the O2 coordination environments by comparing the relative stabilities and activation free energy barriers (ΔG ‡ ) between the corresponding isomeric intermediates of A3 and R4 ( Fig. 1d (Fig. 1d). For S6, A3 is 0.53 eV/intermediate higher in energy than R4 due to the increased numbers of oxygen hole, however, the ΔG ‡ for direct O-O coupling is as high as 1.20 eV ( Fig. 1e and Supplementary Fig. 7). All these calculations suggest that the acid-base nucleophilic attack prevails owing to the high energy penalty for breaking the directional Mn-O bond for the O-O direct coupling. Moreover, because the ligand holes provide the spin-controlled electron transfer channels between catalyst and reactant 23 , the ΔG *OOH for *OOH formation is lowered gradually with the greater number of oxygen holes (Fig. 1f).
Tuning of lattice oxygen reactivity and scaling relation via alkali metal mediation. The number of Na + provides the lever to optimize the overall activity through tuning the lattice oxygen reactivity and scaling relation ( Fig. 2a and Supplementary Fig. 8).
To correlate the change of barrier symmetry resulted from oxygen hole 24 , the electronic parameter (P1) is defined as the variation of ΔG *O -ΔG *OH with reference to S1. As seen, the greater number of oxygen holes, the higher value of P1 ( Fig. 2b and Supplementary Fig. 9a and Tables 2, 3). Due to the strong LSR of ΔG *OOH = ΔG *OH + 3.24 eV, the overpotential ceiling of~0.39 eV is ultimately approached when the magnetization moment (charge) of oxygen ions increases above 0.23 µB (−0.94 e). To demonstrate Na + -specific noncovalent interaction with *OOH 24 , the geometric parameter (P2) is then defined as the variation of ΔG *OOH with reference to S6. As seen, the greater number of Na + around lattice oxygen, the lower value of P2 ( Fig. 2c and Supplementary Tables 2, 4). Owing to the regulation of the LSR, the overpotential ceiling is gradually reduced from 0.39 to 0.19 V. The decreased intersection angles between O-O and H in *OOH further confirm the enhanced electrostatic interaction between *OOH and Na + (Supplementary Fig. 9b).
We further build a dynamic volcano plot to decipher the optimal coordination environment for the lowest OER overpotential ( Fig. 2d and Supplementary Fig. 10 and Tables 2-5) 11,25,26 . As shown in Supplementary Tables 3 and 4, the individual P1 or P2 optimization ([0, P2] or [P1, 0]) cannot guarantee to achieve low overpotential. This is because what defines an ideal OER catalyst is not that the ΔG *OOH -ΔG *OH be 2.46 eV or ΔG *O -ΔG *OH be 1.23 eV but rather that all of free energy change of OER steps are numerically equal to the equilibrium potential 10 . Various overpotentials of 0.37, 0.45, 0.33, and 0.41 eV can be achieved in the coordination environments from S2 to S5 (Fig. 2d) with the variation of both P1 and P2. All  these suggest that tuning of lattice oxygen reactivity and scaling relation via alkali metal mediation can construct better oxygen evolution electrocatalyst. For a better illustration of the overall activity trend, a binary descriptor of P1-P2 that takes the barrier symmetry between the O-H bond cleavage and *OOH formation and overcoming the limitation from LSR into consideration, is further proposed. As shown in Fig. 2e, the overpotentials on different oxygen coordination environments in Na x Mn 3 O 7 are plotted as a volcano-like function of the calculated P1-P2. More specifically, S4 is closest to the apex and thus considered as the most reactive oxygen site for OER. On the left branch, the OER activity is constrained by lattice oxygen reactivity with the ratelimiting step of *OOH formation; on the right branch, the OER activity is constrained by the scaling relation with the ratelimiting step of O-H bond cleavage. All these demonstrate that by tuning the lattice oxygen reactivity and scaling relation, NaMn 3 O 7 is predicted to be the optimal electrocatalyst in the modeled Na x Mn 3 O 7 materials.
Synthesis and characterization of Na x Mn 3 O 7 . Na x Mn 3 O 7 (x = 2.0, 1.5, 1.0, and 0.7) materials with different geometric and electronic environments of lattice oxygen are obtained using a solid-state reaction (see "Methods" and Supplementary Information). XRD pattern of Na 2 Mn 3 O 7 (Fig. 3a) demonstrates that the as-synthesized material is triclinic P1-phase Na 2 Mn 3 O 7 without the presence of impurity phase 17,27 . XRD patterns of other materials ( Supplementary Fig. 11) are similar to Na 2 Mn 3 O 7 with no impurity peak, which is ascribed to the structural flexibility due to the presence of native Mn vacancy in MnO 2 layer 27 . Furthermore, the diffraction peaks (16.3°, 32.5°, and 38.8°) exhibit a little shift to a higher angle, suggesting the gradual lattice contraction with the decreased number of Na + . The Na/ Mn ratios obtained by inductively coupled plasma optical emission spectrometry (ICP-MS) are roughly close to the starting materials (Supplementary Table 6). SEM and transmission electron microscopy (TEM) images combined with energy dispersive spectrometer (EDS) mapping further demonstrate the layered structure of the samples with a homogenous distribution of Na, Mn, and O elements . We further investigated the electronic and coordination structures of Na x Mn 3 O 7 from bulk to surface using hard and soft X-ray absorption spectroscopy (XAS). As shown in X-ray absorption near edge structure (XANES) spectra, the energy of Mn K-edge for Na 2 Mn 3 O 7 is very close to that of the referenced MnO 2 (Fig. 3b), indicating it predominantly consists of Mn 4+ ions. As the number of Na + decreases, the Mn K-edge shifts to a lower energy region and the intensity of white line peak is decreased. This indicates that Mn oxidation state in Na x Mn 3 O 7 is lowered with the presence of oxygen vacancy as the number of Na + decreases 28 . In addition, the almost overlapped Mn K-edge oscillation curves (Fig. 3c) indicate their similar geometric structures, in accordance with XRD patterns ( Supplementary  Fig. 11). Fourier-transformed magnitude plot of the Mn K-edge extended X-ray absorption fine structure (EXAFS) spectra are shown in Fig. 3d, where the two shells at 1.2 and 2.5 Å represent the Mn-O scattering path and Mn-Mn scattering path 29 ,   Table 7), further suggesting the presence of oxygen vacancy and Mn 3+ ions in the bulk 30 . Meanwhile, the negligible changes of Mn-Mn peak imply unchanged Mn-Mn interatomic distance and coordination number, reflecting the strong structural flexibility to accommodate distortion on transition metals from the native Mn vacancy in MnO 2 layer 17,27 . Soft X-ray absorption spectroscopy that is more surface sensitive, were further performed on Na 2 Mn 3 O 7 and NaMn 3 O 7 . From the Mn L-edge spectra (Fig. 3e), the two main peaks at the L III and L II edges, which arises from the transition of a 2p electron to the partially filled 3d shell 31 , are significantly changed in shape (less sharp e g /t 2g feature) and position (toward lower energy direction), suggesting the increased delocalized feature of electronic state and decreased oxidation state of Mn ions for NaMn 3 O 7 compared with Na 2 Mn 3 O 7 . From the O K preedge spectra (Fig. 3f), the characteristic peaks between 528 and 534 eV represent the spectroscopic excitations to the hybridized state of O-2p and Mn-3d, which are split by the crystal field of the local Mn-O coordination geometry 32 . A clear loss of intensity on Mn 3d-O 2p hybridization feature is observed when the number of Na + decreases, implying a decrease of Mn oxidation state on the surface, in line with the results from the Mn K-edge XANES spectra (Fig. 3b). The hybridization parameters (defined as absorbance/(e g holes + 1/4t 2g holes)) of Mn-O bonds of Na 2 Mn 3 O 7 and NaMn 3 O 7 are calculated to be about 0.45 and 0.51 33 . Moreover, the evident shift of the O K pre-edge to a higher energy region further confirms the decrease of Mn oxidation state owing to the increase of oxygen vacancy when the number of Na + decreases, in agreement with the Mn L-edge XAS spectra (Fig. 3e). As mentioned later, the OH -(aq.) tends to spontaneously fill the oxygen vacancy sites of NaMn 3 O 7 under electrochemical OER conditions. On further electrochemical deprotonation, the lattice oxygen ions coordinated with two Mn ions would begin to be oxidized, producing oxygen hole states in|O 2p .
Electrocatalytic measurement for water oxidation. To demonstrate the benefits of alkali metal-mediation, the as-synthesized materials were evaluated for electrocatalytic water oxidation by the rotating disk electrode technique in O 2 -saturated 1 M KOH electrolytes (the experimental details are shown in Supplementary  Fig. 15 and "Methods"). Figure 4a shows the resistance-(iR)corrected polarization curves of Na x Mn 3 O 7 , where the currents are normalized by Brunauer-Emmett-Teller (BET) surface areas to reflect the intrinsic activity ( Supplementary Fig. 16 and Table 6) 34 . The overpotentials for reaching a specific current density of 0.25 mA cm −2 ox are used for activity comparison. Na 2 Mn 3 O 7 shows the lowest activity with an overpotential of 370 mV. As the value of x decreases, the overpotentials are significantly decreased to 300 and 280 mV for Na 1.5 Mn 3 O 7 and NaMn 3 O 7 , respectively. However, further decreasing x significantly lowers the activity, with the high overpotential of 340 mV for Na 0.7 Mn 3 O 7 . Figure 4b further compares the specific current densities at a constant overpotential of 320 mV. As expected, NaMn 3 O 7 gives the highest specific current density of 1.08 mA cm −2 ox , which is 36.0, 2.16, and 12.0 times higher than Na 2 Mn 3 O 7 (0.03 mA cm −2 ox ), Na 1.5 Mn 3 O 7 (0.45 mA cm −2 ox ) and Na 0.7 Mn 3 O 7 (0.09 mA cm −2 ox ), respectively. Considering the intrinsic activity of the catalyst is controlled by both geometric and electronic properties of active sites, the descriptor of P1-P2 is averaged on different oxygen sites in Na x Mn 3 O 7 (Fig. 4b). As seen, P1-P2 scales with the number of Na + in a linear manner, which can explain the optimal performance of NaMn 3 O 7 , consistent with the predictions from DFT calculations (Fig. 2).
Tafel plots (Fig. 4c) are further derived from the polarization curves of Na x Mn 3 O 7 with the previously benchmarked IrO 2 as reference 35 . As seen, the electrocatalytic activity of NaMn 3   overpotentials to achieve the same kinetic OER current density. In addition, the Faradaic efficiency of 97% is measured for the best-performing catalyst of NaMn 3 O 7 , indicating that the measured current is primarily originated from water oxidation. The stability measurements were further carried out for Na x Mn 3 O 7 at a constant overpotential of 320 mV for 10 h. For NaMn 3 O 7 , it maintains 95% of its initial specific activity (Fig. 4d). Similar observations can be made for Na 1.5 Mn 3 O 7 and Na 0.7 Mn 3 O 7 (Supplementary Fig. 17). Besides, no evident change in pseudocapacitive and OER currents for NaMn 3 O 7 during cyclic voltammetry (CV) scans (Fig. 4d, inset) were observed as a good indicator of structural stability of the catalysts 28,36 . Neither evident surface amorphization in HRTEM images nor peak variation in XRD patterns was observed for the post-cycled catalyst after OER measurement ( Supplementary Fig. 18). The absence of signal from the XPS spectra of K 2p ( Supplementary  Fig. 19) on the cycled NaMn 3 O 7 after OER measurement shows that the intercalation of K + is below the detection limit and negligible. ICP-MS test on the electrolyte further demonstrates that no evident leaching of Na and Mn cations as the OER proceeds on Na x Mn 3 O 7 (x = 1.5, 1, and 0.7). Actually, Mn ions are unlikely to be able to migrate to the interlayer sites for NaMn 3 O 7 with the calculated energy barrier as high as 2.80 eV due to the size mismatch between the Mn and Na ions 27,37,38 . Moreover, the charge disproportionation reaction that typically leads to the dissolution of Mn ions is unfavorable in alkaline media 39,40 . As such, we attribute the OER durability and structural stability to the ordered native vacancies in NaMn 3 O 7 that can self-regulate its deformation and electrochemical reversibility 17,27 . On the contrary, Na 2 Mn 3 O 7 shows evident enhanced activity during OER measurement (Supplementary Fig. 17). ICP-MS test on the electrolyte collected after OER measurement on the Na 2 Mn 3 O 7 , demonstrates the evident leaching of Na + after OER measurement, whereas the negligible leaching of Mn cations can be found (Supplementary Table 8). The diffraction peak of the post-cycled Na 2 Mn 3 O 7 exhibits a little shift to a higher angle, confirming the interlayer Na + in the lattice is predominantly leached (Supplementary Fig. 20a). Partial surface amorphization with the thickness of 3-5 nm was also observed from HRTEM image for the post-cycled Na 2 Mn 3 O 7 (Supplementary Fig. 20b). We deduce that the high oxidative OER potential drives the Na + leaching 38 . As a result, the activated lattice oxygen from the Na + leaching contributes to the enhanced activity for Na 2 Mn 3 O 7 .
Verification of oxygen evolution mechanism and active site. Resolving the near-surface structures under electrochemical condition of the catalyst in its highest metastable catalytic state is a prerequisite for the understanding of the OER mechanism and related active site 41 . As such, the in situ X-ray photoelectron spectroscopy (XPS) measurements were performed on NaMn 3 O 7 . From Mn 2p XPS spectra (Fig. 5a), the binding energy shifts to a higher energy at an applied potential of 1.25 V compared with that collected at open circuit, indicating the oxidation state of surface Mn increases. With further increase of potential to 1.55 V, no evident variation of Mn 2p spectra indicates that the Mn ions of the catalyst are structurally and electronically similar to that of 1.25 V. From O 1s XPS spectra, the characteristic peak at 531.2 eV corresponding to oxygen vacancy diminishes at an applied potential of 1.25 V compared with that collected at open circuit. Moreover, the characteristic peak at 529.3 eV corresponding to lattice oxygen shifts to higher energy upon the increase of potential, indicating the oxidation of lattice oxygen. All these indicate that the refilling of oxygen vacancy with OH -(aq.) and the subsequent deprotonation occur before the electrochemical OER process, in agreement with the DFT calculations and in situ XAS measurement (Fig. 5c and Supplementary Figs. 21, 22). As shown in Fig. 5c, OH − (aq.) tends to spontaneously refill the oxygen vacancy sites of Na x Mn 3 O 7 under 1.23 V vs. RHE. From Pourbaix diagram ( Supplementary Fig. 21), the surface termination of Na x Mn 3 O 7 shows the stoichiometric composition under electrochemical OER conditions. Accompanied by the decreased interatomic Mn-O distance from OCV to 1.55 V ( Supplementary  Fig. 21), the more electrophilic oxygen species with oxygen hole can be formed on the highly covalent oxides after the deprotonation and involve in the subsequent O-O bond formation 8,28 . Combining all these results, we demonstrate that activating lattice oxygen leads to the enhanced OER activity, as the Fermi level enters the |O 2p states for Na x Mn 3 O 7 (x < 2) due to the charge compensation and redistribution, creating the reactive oxygen radicals on the surface which behave as electrophilic centers prone to nucleophilic attack from the oxygen lone pairs of OH -. The study of the kinetic isotope effect and pH dependence on OER activity can provide further insight into the reaction mechanism and related intermediates on NaMn 3 O 7 6,42,43 . Because proton mobility in deuterated water solutions can be 1.6-5.0 times slower than that in various protonated water electrolytes, the use of D 2 O can effectively slow down the protontransfer kinetics. As shown in Supplementary Fig. 23, the deuterium kinetic isotope effect (at the potential of 1.55 V) is 2.91 for NaMn 3 O 7 in KOH solution with the increased Tafel slope from 36.4 to 65.1 mV dec −1 , validating the proposed ratelimiting step involves cleavage of the O-H bonds. Such effect is absent in the case of Na 2 Mn 3 O 7 . From pH-dependent OER measurements on RHE scale, NaMn 3 O 7 shows the enhanced activity with the increase of pH from 12.5 to 14. The strong pH dependence indicates the chemical deprotonation step is rate-limiting (Fig. 5d). Correspondingly, one possible OER pathway on NaMn 3 O 7 is proposed in which the deprotonation of *OH involves only proton transfer (M-OH + OH − → M-O − + H 2 O) and is decoupled from the subsequent electron transfer during the release of oxygen (Fig. 5e). Similar phenomenon is also reported for other highly covalent oxides, in which the oxygen redox generates the negatively charged oxygenated species and therefore results in weak OHaffinity 42,44 . The partial charge of such species is not directly measurable, but DFT calculations have indeed shown that these intermediates bind to the catalyst surface as a peroxo-or supero-like anion 6,8,42,45 (Supplementary Fig. 24). To track these charged intermediates on NaMn 3 O 7 during OER, tetramethylammonium cation (TMA + ) as a chemical probe is introduced to the solution because of its specific electrostatic interaction with negative oxygenated intermediates 8,46 . As expected from the Raman spectra (Fig. 5f), there are three new peaks appear at 451, 753, and 951 cm −1 , coinciding with the characteristic peaks of TMA + , when the NaMn 3 O 7 electrode was operated at a constant potential of 1.50 V versus RHE in 1 M tetramethylammonium hydroxide (TMAOH) electrolyte. We further compare the OER activities of NaMn 3 O 7 in 1 M KOH and TMAOH solutions ( Supplementary Fig. 25). A drop in OER activity with the change of Tafel slope from 48.3 to 50.1 mV dec −1 can be observed in the case of TMAOH because of the partial inhibition of the OER, resulted from strong electrostatic interaction between TMA + and negative oxygenated intermediates.  In summary, through a combination of in-situ spectroscopybased characterization and first-principles calculations, we have employed modeled Na x Mn 3 O 7 materials to illustrate the crucial roles of alkali metal mediation to tune lattice oxygen reactivity and scaling relation for the rational design of better OER electrocatalysts. More specifically, an intermediate level of Na + mediation (NaMn 3 O 7 ) manifests the optimal activity, which is due to the regulation of the barrier symmetry between O-H bond cleavage and *OOH formation on the basis of overcoming the preexisting scaling relation. In addition, the pH-dependent experiment and Raman spectra further demonstrate it works in a decoupled proton/electron route with the presence of negatively charged oxidized oxygen species. This work provides a guideline for the rational design of better catalysts towards electrocatalytic water oxidation or other oxidative reactions through tuning lattice oxygen reactivity and scaling relation.

Methods
Systhesis of Na x Mn 3 O 7 . Na x Mn 3 O 7 (x = 2.0, 1.5, 1.0 and 0.7) materials were synthesized by a modified method of solid-state reaction 17,27,47 . In detail, the starting materials of NaNO 3 and MnCO 3 with the desired Na/Mn ratios were thoroughly mixed in an agate mortar and pressed into pellets under pressure of 10 MPa. Then, the Na x Mn 3 O 7 materials can be obtained by calcining such pellets at 500-650°C in tube furnace with O 2 atmosphere for 5-10 h. Before use, the Na x Mn 3 O 7 materials were stored in a glovebox with Ar atmosphere.
Composition and structure characterization. X-ray diffraction (XRD) patterns were collected with a Bruker D8 FOCUS equipped with nickel-filtered Cu Kα radiation (λ = 1.541 Å). Field emission SEM characterization was performed with a Hitachi S-4800 SEM. TEM and EDS characterizations were performed with a JEM-2100F transmission electron microscope. Elemental composition was analyzed using the techniques of Vista-MPX EL02115765 Inductively coupled plasma spectrometry and PerkinElmer NexION 350× Inductively coupled plasma-mass spectrometry (ICP-MS). For ICP-MS measurement, the standard curve is linearly fitted in the range of 0.1, 1.0, 10.0, 100, 500, and 1000 ppb with the internal standard of Rh (10 μg L −1 ). The BET surface areas were obtained from N 2 sorption isotherm measurements on Micrometrics TriStar 3000 equipment.
Mn K-edge X-ray absorption spectra (XAS) were collected at the X-ray Absorption Fine structure for catalysis (XAFCA) beamline at the Singapore Synchrotron Light Source (SSLS) using the transmission mode 48 . The photon energy is ranged from 1.2 to 12.8 keV achieved by two sets of monochromator crystals of Si (111) and KTiOPO 4 crystal. In-situ XAS measurements were conducted with a home-made electrochemical cell in fluorescence mode. The fluorescence yields were collected with silicon drift detector (Bruker Xflash 6|100). The catalysts were coated on carbon paper with mass loading of 2 mg cm −2 as the working electrode. Acquired EXAFS data were analyzed using ATHENA module implemented in the IFEFFIT software packages 49 . The k-and R-ranges to fit the EXAFS data were set as 2-12 Å −1 and 1.0-3.5 Å, respectively. Mn L-edge and O Kedge XAS were collected at Soft X-ray-ultraviolet (SUV) beamline at SSLS. X-ray photoelectron spectroscopy (XPS) characterizations were performed on a Thermo ESCALAB 250Xi X-ray photoelectron spectroscope using the home-made X-ray cell. At different applied potentials, the working electrodes were first stabilized to reach a steady state, then the corresponding XPS signals were collected and analyzed.
Electrochemical measurements. Electrochemical measurements were conducted in a three-electrode setup with graphite rod and Hg/HgO (1 M KOH) as the counter and reference electrode. The working electrode was prepared by coating catalyst ink on a glassy carbon with mass loading of 0.204 mg ox cm −2 . Typically, 4 mg catalyst and 0.8 mg acetylene black were suspended in 2 mL mixture solution of isopropanol, water and Nafion and violently ultrasounded for 3 h to form a homogeneous ink. To decrease the influence of capacitive current and gas bubbles, the linear sweep voltammetry measurements were performed with the scan rate of 2 mV s −1 and rotational speed of 1600 r.p.m. The Tafel plots were derived from the polarization curves as the function between overpotential and the log current (η = blog[J] + a), where b represents as Tafel slope. CV measurements were conducted at a scan rate of 10 mV s −1 for NaMn 3 O 7 to investigate the pseudocapacitive charge preceding the OER region. All the used potentials were calibrated based on the RHE. The Faradic efficiency was measured using in-line gas chromatograph, which is defined as the ratio between the amount of experimentally measured O 2 and the amount of theoretically produced O 2 from the reaction.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request. Source data are provided with this paper.