Na 2.4 Al 0.4 Mn 2.6 O 7 anionic redox cathode material for sodium-ion batteries – a combined experimental and theoretical approach to elucidate its charge storage mechanism †

Here we report the synthesis via ceramic methods of the high-performance Mn-rich Na 2


Introduction
Sodium-ion batteries (SIBs) are a promising candidate for nextgeneration batteries, as sodium is an abundant element with an easy-access source, which translates into a low-cost and sustainable alternative energy storage technology.Moreover, the chemical similarities of Li and Na coupled with the already well-established knowledge in lithium-ion battery (LIB) technology allow for rapid implementation of SIBs in applications such as the grid or low-cost transportation. 1,2Among the cathode materials researched in the eld of SIBs, layered Na-Mn-O ternary compounds have been particularly relevant due to their high capacity, abundance, non-toxicity and low cost.
The Na x MnO 2 (x > 0.5) composition has been the most studied, showing a maximum theoretical capacity of 244 mA h g À1 based on the Mn 3+/4+ redox-active couple.Nonetheless, these materials generally present low capacity retention, e.g., b-NaMnO 2 shows 53% capacity retention at C/20 aer 100 cycles, with an initial capacity of 175 mA h g À1 (when cycled between 4.2 V and 2.0 V); 3 and a-NaMnO 2 shows 70% capacity retention at C/30 aer 20 cycles, with an initial capacity of 185 mA h g À1 (when cycled between 3.8 V and 2.0 V). 4 The poor cycling stability of Mn-based compounds is generally attributed to structural instabilities associated with the Jahn-Teller (JT) effect of the Mn 3+ ions formed during cycling, particularly under a low working voltage range (1.5-2.1 V). 5 Na 2 Mn 3 O 7 (or Na 0.57 Mn 0.86 [ ] 0.14 O 2 ) (P 1 space group) has been extensively studied in the past years due to the large anionic redox contribution (75 mA h g À1 above 4.0 V vs. Na + /Na) to the overall charge stored in the material (theoretical capacity ¼ 155 mA h g À1 ), which results in high energy densities. 6,7][10][11] Despite its large anionic redox capacity, Na 2 Mn 3 O 7 does not fully use its redoxactive oxygen ions.Earlier studies on Na 2 Mn 3 O 7 revealed that three non-bonding oxygen atoms per formula unit are redoxactive, where theoretically, a maximum of only two of them can be oxidized/reduced during cycling. 6,8This shows that, once all non-bonding oxygen atoms of Na 2 Mn 3 O 7 become active, its anionic redox capacity can be increased by z50%.One effective way to further use the redox-active oxygen atoms is by aliovalent substitution, using elements with oxidation states lower than Mn 4+ , such as Mg 2+ , B 3+ , Al 3+ , etc. Thus, to satisfy charge neutrality, the aliovalent-doped Na 2+a M y Mn 3Ày O 7 (M ¼ elements with oxidation states lower than 4+) should have a greater number of Na + ions or, equivalently, higher anionic redox capacity.Another approach to enhance anionic redox reactions is by substitution of TMs with more electronegative elements. 12,13This concept stems from the fact that highly electronegative dopants can weaken the ionicity of Mn-O bonds and destabilize the energy of non-bonding O 2p bands, enabling anionic redox reactions.
In this regard, Al 3+ ions, with low oxidation state and high electronegativity compared to Mn 4+ ions, are ideal substituents to improve anionic redox reactions occurring in Na 2 Mn 3 O 7 .9][20][21] This causes a lack of understanding of the effect of Al 3+ ion substitution on the electrochemical performance of cathode materials with anionic redox reactions.][24][25] In this work, we aim to investigate the role of Al substitution in the charge storage mechanism of the Na 2 Mn 3 O 7 cathode.To this end, using Na 2.4 Al 0.4 Mn 2.6 O 7 as a testbed material, we study the effect of Al 3+ ion substitution on structural changes, such as the interlayer spacing and Na + ion insertion sequences.Using electrochemical, structural and hybrid density functional theory (DFT) calculations, we analyse the changes in Mn/oxygen redox reactions and associated O 2 release during battery cycling.Based on these ndings, we discuss the role of Al on the electrochemical behaviour of Mn-based cathodes and its potential application for developing high-performance cathodes for SIBs.
Powders were pressed into 13 mm pellets at 5 tones per cm 2 and calcined under an oxygen ow at 700 C for 12 h followed by 1000 C for 6 h, with an intermediate ball-milling step at 400 rpm for 1 h for the Na 2.4 Al 0.4 Mn 2.6 O 7 synthesis, and under an oxygen ow at 600 C for 4 h for the Na 2 Mn 3 O 7 synthesis.The heating and cooling rate used during the annealing process was 10 C min À1 .
Powder X-ray diffraction (PXRD) data for Na 2.4 Al 0.4 Mn 2.6 O 7 and Na 2 Mn 3 O 7 pristine materials, and cast Na 2.4 Al 0.4 Mn 2.6 O 7 electrodes at different states of charge were collected at ambient temperature from 10 to 90 2q degrees in a q/2q Bragg-Brentano geometry using an X-ray diffractometer (Miniex, Rigaku) with a Cu X-ray source (0.6 kW) (l K a1 ¼ 1.5406 Å, l K a2 ¼ 1.5444 Å) equipped with a D/teX-ULTRA 250 high-speed position sensitive detector system.The X-ray diffractometer was enclosed in an N 2 -lled MBraun glovebox.For the ex situ measurements, electrodes were galvanostatically cycled at 10 mA g À1 to selected voltages, extracted from the coin cells, rinsed three times with diethyl carbonate (DEC), and dried at 80 C under vacuum in the antechamber of an argon-lled glovebox (MBraun) ([H 2 O] and [O 2 ] # 0.1 ppm).To disregard any preferred orientation in the asprepared Na 2.4 Al 0.4 Mn 2.6 O 7 material, XRD data were also acquired in transmission mode, using a spinning borosilicate capillary (0.5 mm diameter) placed in a 4-circle Kappa precision goniometer in an Agilent Supernova diffractometer with a Mo Xray source and an Atlas S2 135 mm CCD area detector.XRD data were t with the Le Bail method, using the FullProf Suite soware. 26VESTA 3.5.7 soware was used to visualise the crystal structures. 27ample microstructure was analysed using eld emission scanning electron microscopy (FE-SEM) with a JEOL JSM-7800F operating at 5.0 kV.Elemental analysis was performed at 10 kV using X-ray Energy Dispersive Spectrometer (EDS) (X-Max50, large area 50 mm 2 Silicon Dri Detector (SDD) from Oxford Instruments).Sample preparation consisted of uniformly coating the samples onto carbon tabs (G3348N, Agar Scientic).
X-ray uorescence (XRF) spectrometry was carried out using an EDXRF Shimadzu EDX 8000 spectrometer with an SSD detector and an Rh X-ray source.The sample was pelletized and manipulated under an argon atmosphere without air exposure.Quantication of elemental Na was not possible due to the overlapping signal with Rh escape peaks originating from the Xray source.

Electrochemical characterisation
The electrochemical properties of Na 2.4 Al 0.4 Mn 2.6 O 7 and Na 2 Mn 3 O 7 were assessed in sodium half-cells using stainless steel 2032-type coin-cells (Tob New Energy).Electrode preparation and coin-cell assembly were carried out in an argon-lled glovebox (MBraun) ([H 2 O] and [O 2 ] # 0.1 ppm).Electrode slurries were rst prepared by mixing a homogenous powdered mixture (500 mg) consisting of 80 wt% active material (Na 2.4 -Al 0.4 Mn 2.6 O 7 ), 10 wt% carbon black (>99.0%,Alfa Aesar) and 10 wt% polyvinylidene uoride (PVDF) (>99.0%,Alfa Aesar) with a mortar and pestle.Then, the powdered mixture was dissolved in ca.1.2 mL N-methyl-2-pyrrolidone (NMP) (anhydrous, >99.0%,Alfa Aesar) and the homogeneous slurry mixture was cast onto an aluminium foil using a doctor blade and then dried under vacuum at 80 C for 12 h in the glovebox antechamber.Once dried, the 12 mm diameter electrodes were cut and pressed under z174 MPa, with a typical active material mass loading of 2.4-3.0 mg cm À2 .Sodium metal disks (15 mm in diameter) (Alfa Aesar Merck) were used as reference electrodes, and 1 M NaPF 6 (99%, Alfa Aesar) in the anhydrous organic solution of ethylene carbonate/diethylene carbonate (EC : DEC 1 : 1 w/w%, Gotion) as electrolyte.The electrolyte was dried under molecular sieves for a week before use.
Galvanostatic measurements were carried out on a battery tester (Neware) in the voltage range of 1.5-4.7 V vs. Na + /Na at selected current rates (C/20 to 5/C, 1C ¼ 160 and 155 mA g À1 for Na 2.4 Al 0.4 Mn 2.6 O 7 and Na 2 Mn 3 O 7 , respectively).
Cyclic voltammetry (CV) experiments were conducted using an IviumStat potentiostat (Alvatek), at a scan rate of 0.1 mV s À1 in the voltage range of 1.5-4.7 V vs. Na + /Na.
Galvanostatic intermittent titration technique (GITT) measurements during charge/discharge were carried out by applying a pulse of AE1.7 mA for 1 h followed by an open circuit relaxation of 2 h.The procedure was carried out in the voltage range of 1.5-4.7 V vs. Na + /Na.The diffusion coefficient for Na 2.4 Al 0.4 Mn 2.6 O 7 , D, was calculated using eqn (S1), † where we assumed a constant molar mass (320.81 g mol À1 ) and cell volume (169.81 cm 3 mol À1 ) during the cycling process.The surface area was calculated using two different approaches: (1) using the average particle size value obtained from SEM images (0.77 mm diameter particles), whereby we assumed that particles were spherical and had a smooth surface and an intimate contact with the electrolyte; and (2) considering a at surface contact area (electrode area).
Electrochemical impedance spectroscopy (EIS) studies were carried out on Na 2.4 Al 0.4 Mn 2.6 O 2 and Na 2 Mn 3 O 7 with two different electrochemical cell set-ups.EIS data were collected in the 100 kHz to 10 mHz frequency range with an amplitude of 10 mV and aer a 1 h voltage pulse using a Biologic VMP-300 potentiostat and an Ivium Octostat potentiostat.EIS data on the rst cycle were collected using a three-electrode cell (EL-cell), where the active electrode materials were cycled using metallic Na disks as counter and reference electrodes.The referece is a pre-assembled and commercial Na ring (El-cell) enclosed in the cell around the glass bre separator.On the other hand, EIS data on the 100 th cycle were collected using two-electrode coin cells.
All the electrochemical testing was conducted at room temperature.
Nuclear magnetic resonance (NMR) 23 Na magic-angle spinning (MAS) NMR spectra of pristine  23 Na MAS NMR spectra were acquired using a single pulse or Hahn echo pulse sequence, with a recycle delay of 0.5 s.Spectra referencing was done using a secondary solid reference of NaCl (d iso ¼ 7.2 ppm).Due to the nature of the samples and the experimental parameters used, the relative intensities may not match contributions in the materials and therefore spectra are not quantitative.Additionally, the distorted baseline of the single pulse experiments can affect the observed position of peaks compared to the atter baseline of the Hahn echo experiments.Ex situ 23 Na MAS NMR spectra were obtained at 9.4 T, using a Hahn echo pulse sequence and 30 kHz MAS. 1 kHz line broadening was used in spectra processing.Ex situ samples were run in powder form in coin cells (active material : carbon black weight ratio of 80 : 20) using around 8 mg of active material per cell.Na metal was used as a counter and reference electrode and 1 M NaPF 6 in EC : DEC (1 : 1 vol%) as the electrolyte.Aer cycling the cells to a dened voltage, the cells were opened, and powder was recovered.Powders were washed three times with DEC inside an argon-lled glovebox (MBraun), followed by the removal of excess solvent and drying under vacuum at room temperature in the glovebox antechamber before the NMR measurement.
Ex situ synchrotron X-ray absorption near-edge spectroscopy (XANES) X-ray absorption near edge structure (XANES) data were collected at the B18 beamline at Diamond Light Source (Harwell, UK).For the measurement, cast electrodes were cycled at different states of charge, and then, extracted from the coin cells and rinsed three times with DEC in an argon-lled glovebox (MBraun) ([H 2 O] and [O 2 ] # 0.1 ppm).Aer this, electrodes were dried at 80 C under vacuum in the glovebox antechamber and sealed in individual aluminium laminated pouches under vacuum to minimise air exposure during data acquisition.Measurements were performed at the Mn K-edge energies above and below the absorption edges (ca.6540 eV) and electrode data were collected in uorescence mode, using a 36-element Ge detector.Mn metal foil was used as an internal standard and reference data for MnO, Mn 3 O 2 and MnO 2 powders were collected on transmission mode using gas-lled ionisation chambers to measure both incident and transmitted intensities.Three repetitions of each sample were measured, totalling a 15 min measurement.Athena soware in the Demeter package was used to perform energy calibration, background subtraction and normalization. 28,29Linear combination tting (LCF) was performed using the Athena soware within an energy range of À20 eV below to +30 eV above the edge.Each sample spectra were tted for every combination of the standard spectra.The combination of standards that resulted in the best t (smallest R-factor and reduced c 2 ) was chosen as the most likely representation of the sample.

Ab initio calculations
All the calculations performed in this work employed DFT as implemented in the Vienna Ab initio Simulation Package code. 30,31Interactions between core and valence electrons were described using the projector augmented wave (PAW) method. 32a (2p 6 3s 1 ), Mn (3p 6 3d 5 4s 2 ), Al (3s 2 3p 1 ), and O (2s 2 2p 4 ) electron congurations were treated as the valence electrons.Convergence to plane wave energy was checked, with a cut-off of 500 eV found to be sufficient to converge the total energy to within 0.01 eV atom À1 .Brillouin zones for all compounds were sampled such that the k-points were converged with an accuracy of the total energy in 0.001 eV per atom, which corresponds to k-grids of 3 Â 3 Â 3 and 2 Â 2 Â 2 for primitive (24 atoms) and supercell (120 atoms) structures.
To analyse the Na + ion insertion sequences and theoretical voltages of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 upon cycling, we predicted the optimized geometries of partially and fully sodiated Na x Mn 3 O 7 and Na x Al 0.4 Mn 2.6 O 7 using the PBE sol +U method. 33This method provides computationally cost-effective corrections on highly correlated Mn d orbitals, which allows rapid modeling of materials with TMs.The detailed procedures are described in two sections of the ESI, † i.e., Computational details for the construction of model structures (Fig. S1 and S2 †) and Computational details for the theoretical voltage calculations (Fig. S3 †).
While the PBE sol +U method has been shown to provide an accurate geometry for a wide variety of solid systems with a relatively low computational cost, it can result in quantitatively incorrect descriptions of electronic structures, especially when predicting the defect state of highly localized electrons. 34,35Therefore, to better simulate O loss behaviour during desodiation of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 , we used the screened hybrid functional (HSE06), 36 with 25% of exact nonlocal Fock exchange added to the PBE 37 functional, that consistently produces more accurate electronic structures than standard DFT approaches. 38The systems before O 2 gas release were rst modelled by selecting four different structures of Na 2 Mn 3 O 7 , Na 2.4 Al 0.4 Mn 2.6 O 7 , NaMn 3 O 7 and Na 1.4 Al 0.4 Mn 2.6 O 7 from PBE sol +U optimized supercells, and relaxing their structures again with HSE06 functional.The desodiation and associated O loss were then simulated by creating Na and O vacancies in each supercell, followed by the structural relaxations of defective structures.Calculations for defective structures were also done with HSE06 functional to make total energies consistent with host supercells.All Na and O sites of host supercells were classied based on their local bonding environment (Tables 1 and 2).The energy required for desodiation, E f (D), was estimated by calculating defect formation energies for Na vacancies for individual Na sites according to: where E tot (bulk) and E tot (D) are the total energies of the supercells before and aer defect formation, respectively.E tot (Na) corresponds to the total energy of a Na atom in Na metal, whereas n Na is the number of Na atoms removed from host supercells.
To calculate the minimum energy path and corresponding activation barrier for O 2 dimer formation aer desodiation, climbing image nudged elastic band calculations were carried out for partially desodiated NaMn 3 O 7 using HSE06 functional.For this purpose, we rst optimized the geometry of NaMn 3 O 7 before and aer O 2 dimer formation.Using the optimized geometries as initial and nal states, ve intermediate images were then generated using linear interpolation.Lastly, the minimum energy path was searched by minimizing the spring forces parallel to the transition pathways and real forces projected perpendicular to the strings.For all HSE06 calculations, the Brillouin zone was sampled using only the G point.

Crystal structure and microstructure
The TM-decient layered compound Na 2.4 Mn 2.6 Al 0.4 O 7 (Na 0.68 -Mn 0.74 Al 0.11 [ ] 0.15 O 2 in conventional Na x MO 2 notation) was synthesised by solid-state reaction methods, as detailed in the Experimental section.Powder X-ray diffraction (XRD) data of asprepared Na 2.4 Mn 2.6 Al 0.4 O 7 show a sharp and well-dened ( 101) reection at 2q ¼ 15.9 and broader and less intense reections at higher 2q values (Fig. 1a).The latter implies a high level of long-range layer stacking disorder combined with poor crystallinity, hindering an accurate structural renement using the Rietveld method.Neutron and X-ray diffraction studies on the Na 2 Mn 3 O 7 parent structure have shown that the [Mn 3 O 7 ] 2À layers are not perfectly aligned throughout the structure, losing the stacking sequence across the c direction. 39XRD data were t with the Le Bail method, using published crystallographic data on Na 2 Mn 3 O 7 . 6As a result, XRD data were indexed in the triclinic P 1 space group, with cell parameters a ¼ 6.587(2) Å, b ¼ 6.984(2) Å and c ¼ 7.540(3) Å, and a ¼ 106.07(2) , b ¼ 107.32(2) , and g ¼ 110.39 (3) .R wp and R p agreement factors are 3.70 and 2.02%, respectively, indicating a good quality t.No secondary phases were observed in the X-ray diffractogram, suggesting high sample purity.Na 2.4 Mn 2.6 Al 0.4 O 7 shows a slightly larger cell volume of 280.94(9)1][42] The presence of Al 3+ ions does not affect the volume of the crystallographic unit cell since both Al 3+ and Mn 4+ ions have a very similar ionic radius when they are octahedrally coordinated with oxygen (z0.53 Å). 43 Subsequently, we attribute these differences to the extra Na + ions per formula unit (0.4 Na + ) in Na 2.4 Mn 2.6 Al 0.4 O 7 .These additional Na + ions are required to charge-balance the material, due to the aliovalent substitution of Mn 4+ ions by Al 3+ ions.For comparison purposes, only a 0.3% increase in volume per formula unit was observed in a similar layered material, LiNi 0.7 Co 0.3 O 2 , when substituting Co 3+ ions (z0.54 Å) with 10% Al 3+ ions. 20XRD data were also collected in a spinning capillary to discard any preferred orientation in the sample attributed to the at sample holder used (Fig. S5 †).Similar relative X-ray diffraction peak intensities were observed in both datasets, suggesting limited preferred orientation.The FE-SEM image of pristine Na 2.4 Mn 2.6 Al 0.4 O 7 shows a mixture of lamellar and irregular primary particles with an average diameter size of 0.77 mm that aggregate into larger irregular-shaped secondary particles with a broad size distribution between 4.2 and 18.1 mm (average of 9.1 mm) (Fig. 1b).EDS elemental mapping on selected particles shows a homogeneous distribution of sodium, aluminium, manganese, and oxygen atoms (Fig. 1c and S6 †), which conrms the presence of a single-phase product, as observed in the X-ray diffraction data (Fig. 1a).A Mn : Al ratio of 2.58 : 0.42 was determined from XRF spectrometry in the sample, which is close to the expected nominal stoichiometry. 23Na solid-state NMR spectra were obtained to examine local Na environments of Na 2.4 Al 0.4 Mn 2.6 O 7 and Na 2 Mn 3 O 7 , for comparison. 23Na MAS NMR spectra for Na-Mn-M-O materials (where M ¼ e.g., Mg and Al) typically show peaks that can be assigned to details such as layering arrangements, defects, and Mn oxidation states.S7. † The spectrum for Na 2 Mn 3 O 7 (Fig. S7a †) shows an isotropic peak at $1500 ppm.This is consistent with a large Fermi contact interaction that is expected for Na within a predominately Mn 4+ setting, as expected for this material. 3,44Two sharp peaks close to 0 ppm surrounded by a spinning sideband manifold are also observed.
Resonances close to 0 ppm suggest the presence of Na environments within a diamagnetic impurity phase. 3Although no secondary phases were observed in the XRD data, one explanation could be a low-level amorphous phase.The 23  setting.In summary, the comparison of the 23 Na spectra for both materials shows the extra local environments for Na + ions and lower shis corresponding to Mn 3+ upon Al 3+ ion substitution, where the former and latter results corroborate our DFT data (Fig. 1e, f and 4) and XANES data (Fig. S10 †), respectively.However, the spectral broadening prevents more detailed analysis.
To gain further insight into the crystal structure of Na 2.4 -Mn 2.6 Al 0.4 O 7 , DFT calculations were used to determine the atomic positions of additional Na + ions within the structure, assuming a disordered distribution of Al 3+ and Mn 4+ ions within the 2i sites (0.07, 0.07, 0.22) in the vacancy-ordered TM layer. 8,42The insertion sequences of Na  with ten formula units.In the supercell, four different Na sites were observed: three different 2i sites in the Na layer (Na1 (0.25, 0.70 and 0.50), Na2 (0.10, 0.6, 0.02) and Na4 (0.45, 0.72 and 0.87)), which may accommodate up to ten Na + ions per site; and one 1 h site in the TM layer, Na3 (0.5, 0.5, 0.5), which may uptake ve Na + ions.These sites are identical to those observed in Na 2 Mn 3 O 7 , 6,39 except for the Na2 site which is found in a slightly different position i.e., Na2 (0.33, 0.83, 0.07) due to excess Na in sites Na3 and Na4, which may repel some preexisting Na + ions in the Na2 site.
Further computational studies show that Na1 and Na2 sites of Na x Al 0.4 Mn 2.6 O 7 are occupied sequentially until x ¼ 2, with a slight difference in Na active sites near x ¼ 1 compared to Na x Mn 3 O 7 (Fig. 1e and f).Further Na + ion insertion occurs into the vacancy sites in the TM layer (Na3 site) until Na content reaches x ¼ 2.2.Na + ion insertion into vacancies in the TM layer may increase the TM-Na repulsion within layers and Na-Na repulsion between adjacent layers, increasing the cell volume, as observed in the XRD data (Fig. 1a).Last, Na + ion insertion occurs into site Na4 up to x ¼ 2.4.The presence of Na in sites Na3 and Na4 causes slight displacements of Na sites such that the Td (site Na1) and Oh (site Na2) sites become trigonal prismatic (Table 1).The structure of Na 2.4 Al 0.4 Mn 2.6 O 7 , with Na + ions (ca.0.2 Na) located in the TM layer resembles that of other cathode materials with demonstrated oxygen redox activity such as Li-rich layered oxides [45][46][47] and other Na layered compounds. 48,49

Charge carrier pathways and structural changes
To provide a better insight into the changes in diffusion pathways aer Al 3+ ion substitution, we represented the potential Na + ion diffusion pathways of both Na x Mn 3 O 7 and Na x Al 0.4 -Mn 2.6 O 7 .This was done by classifying the active Na sites in terms of their shape and displaying them on the TM layer  (Fig. 2).In Na x Mn 3 O 7 , active Na sites at low sodiation levels (0 < x # 1.0) are characterized by prismatic Na1 sites above unoccupied Mn sites (Fig. 2a and b).Based on the distribution of active Na sites, two diffusion pathways can be suggested: inplane 1D diffusion along with the prismatic Na1 sites, and interplanar diffusion through unoccupied Mn sites.As Na + ions continue to intercalate into Na x Mn 3 O 7 (1.0 < x # 2.0), repulsion between Na + ions increases and pushes Na1 sites toward unoccupied Mn sites, changing their coordination to tetrahedral.The so-formed tetrahedral Na1 sites are likely to be xed under x $ 1.0 (Fig. 1e) and block the diffusion pathways, which slows down the diffusion kinetics.Na + ion diffusion in highly sodiated Na x Mn 3 O 7 , therefore, tends to follow 1D diffusion pathways through octahedral Na2 sites, which is expected to be more sluggish than prismatic-prismatic Na + ion diffusion at low sodiation content (0 < x # 1.0). 50n the case of Na x Al 0.4 Mn 2.6 O 7 , Na + ion diffusion at high voltages (0 < x # 0.7) follows two diffusion pathways, similar to Na x Mn 3 O 7 under 0 < x # 1.0 (Fig. 2c).However, with an additional supply of Na + ions (0.7 < x # 1.2), both prismatic Na1 and Na2 sites become active, unlike in Na x Mn 3 O 7 (Fig. 1f).This is speculated to be due to the higher electronegativity of Al compared to Mn, which destabilizes O 2p orbitals of Al-O bonds and, in turn, makes Na + ions stable on site 2 above AlO 6 octahedra.These active Na2 sites bridge 1D diffusion pathways between Na1 sites and allow a broader range of diffusion along the in-planar direction, resulting in enhanced rate performance of Na x Al 0.4 -Mn 2.6 O 7 .Under highly sodiated conditions (1.2 < x # 2.0), similar to Na x Mn 3 O 7 , Na1 sites are xed as tetrahedral sites near unoccupied Mn sites, causing Na + ion diffusion to occur via octahedral Na2 sites.Once all Na1 and Na2 sites are occupied, excess Na + ions intercalate into the unoccupied Mn site (i.e., Na3) and nearby Na2 site (i.e., Na4) (Fig. 2d).This process is speculated to occur with the cooperative motion of pre-existing Na1 or Na2 as most free spaces are already occupied by Na1 and Na2.
The above discussion suggests that Na + ion diffusion pathways are widened due to the high electronegativity of Al 3+ ions.Moreover, subsequent structural analysis on the Al-substituted Na 2 Mn 3 O 7 cathode revealed that Al shortens the lattice along the in-planar direction of the TM layers, while increasing the interlayer spacing of these TM layers (see Discussion on the structural changes aer Al substitution section and Fig. S11 and S12 in ESI †).Since most Na + ion diffusion pathways of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 cathodes are aligned along the in-planar direction, these structural changes will shorten the diffusion pathways while widening the diffusion channels, further improving the Na + ion diffusion kinetics.This may explain the enhanced rate performance observed for Na 2.4 Al 0.4 Mn 2.6 O 7 Fig. 3  during rate capability tests, as it will be later discussed (Fig. 3d), and agrees with experimental D Na + values obtained from GITT data (Fig. S13 and S14 †), which are among the largest reported for sodium layered oxides (e.g., D z 10 À11 cm 2 s À1 for P2-Na 7/9 Cu 2/ 9 Fe 1/9 Mn 2/3 O 2 and D z 10 À12 to 10 À11 cm 2 s À1 for Na 0.7 Mn 0.7 -Ni 0.3Àx Co x O 2 , x ¼ 0, 0.1, 0.3). 51,52ectrochemical performance Fig. 3a shows cyclic voltammetry data of Na 2.4 Mn 2.6 Al 0.4 O 7 in the 1 st , 2 nd , 5 th , and 10 th cycles in the voltage range of 1.5-4.7 V vs. Na + /Na at a scan rate of 0.1 mV s À1 .Several redox peaks are observed in the CV data, indicating a multistep Na + insertion/ extraction process.The rst anodic sweep shows three distinct regions with oxidation peaks at potentials of 2.8 V (A ox ), 3.2 V and 3.6 V (B ox region), and above 4 V (C ox region).By direct comparison with other manganese-based systems such as P2-Na 2/3 -Mn 1Àx Al x O 2 , 16 Na 0.67 [Fe 0.5 Mn 0.5 ]O 2 53 and P2-Na 2/3 Mn 0.9 Mg 0.1 O 2 , 44 the rst anodic peak (A ox ) can be related to the oxidation reaction of Mn 3+ /Mn 4+ ions in the structure.Data are in accordance with NMR and XANES data, which show Mn 3+ /Mn 4+ cation mixing in the pristine sample (Fig. S9 and S10c †).In subsequent cycles, A ox shis to higher potentials upon cycling (by z0.15 V) and its maximum peak current intensity decreases upon cycling from 75 mA g À1 in the 2 nd cycle to 58 mA g À1 in the 10 th cycle.On the cathodic sweep, the corresponding A red reduction peak appears at 1.9 V and shis to higher potentials upon cycling, occurring at 2.2 V at the 10 th cycle.The increased polarisation could be attributed to structural changes occurring upon cycling. 54Anodic and cathodic peaks in the 2.5-3.7 V voltage range (region B in Fig. 3a) did not change in potential over the rst 10 cycles.GITT data in this voltage range showed a decrease in D Na + values by two orders of magnitude, which might be attributed to sodium/vacancy ordering structural arrangements (Fig. S14 †).The observed anodic peaks in region C (C ox ) are similar to those attributed to an oxygen redox reaction arising from the presence of non-bonding O 2p orbitals in Na 2 Mn 3 O 7 (Fig. S15a †). 6,8The fast loss of peak current intensity upon cycling together with an increase in voltage polarization reects the partial irreversibility of the process due to phenomena such as oxygen gas formation from the particle surface. 55Nevertheless, earlier works on cathode materials with oxygen-redox activity, e.g., [56][57][58] have shown oxygen redox activity, even though the associated anionic redox-related plateaus (in galvanostatic data) or current peaks (in CV data) are not discernible anymore.Therefore, we cannot rule out oxygen redox activity in Na 2.4 Al 0.4 Mn 2.6 O 7 aer the rst cycle by the absence of the highvoltage redox peaks in the CV data.
Galvanostatic charge/discharge proles of the rst two cycles of Na 2.4 Mn 2.6 Al 0.4 O 7 in the voltage range 1.5-4.7 V vs. Na + /Na at a constant current of 8 mA g À1 (C/20 rate) are shown in Fig. 3b.The material was initially charged from its open-circuit value (OCV) (z2.65 V) to 4.7 V, showing a capacity of 163 mA h g À1 which corresponds to the extraction of two Na + ions per formula unit to yield Na 0.4 Mn 2.6 Al 0.4 O 7 .The voltage prole starts with a moderate slope from 2.8 V to 4 V followed by a plateau that terminates at 4.7 V.The voltage proles are signicantly different to our experimental (Fig. S15b †) and literature data on Na 2 Mn 3 O 7 , which shows two high-voltage plateaus during charge. 6,8alvanostatic data for Na 2.4 Mn 2.6 Al 0.4 O 7 match the cyclic voltammetry data, which shows the oxidation of Mn 3+ ions to Mn 4+ starting at low potentials (z2.8 V), followed by the oxidation of oxygen anions at high-potentials (V > 4.0 V) (Fig. 3a).The sloppy load prole of Na 2.4 Mn 2.6 Al 0.4 O 7 compared to reported load curves for Na 2 Mn 3 O 7 indicates that Al 3+ ion substitution for Mn 4+ ions in the TM layer alleviates the structural transitions occurring upon cycling in Na 2 Mn 3 O 7 . 8,59On the subsequent discharge cycle, the same processes are observed: the reduction of oxygen anions above 4.0 V, reected as a short plateau, and the reduction of Mn 4+ to Mn 3+ ions below this voltage.A higher discharge capacity (i.e., 215 mA h g À1 ) compared to the rst charge capacity, which results in a coulombic efficiency (CE) of 130%, was observed (see Fig. 3c).This can be explained by the surplus of Mn 3+ ions formed and the ability of uptaking a higher sodium content compared to the as-synthesised material, due to the low cut-off potential choice.The charge capacity above 4 V is reduced from the rst cycle (55 mA h g À1 ) to the second cycle (26 mA h g À1 ), indicating that the high-voltage redox processes are not fully reversible, most likely reecting the irreversibility of the oxygen redox process and/or oxygen evolution reaction.The reversibility of the reaction will be further discussed with the help of ab initio calculations.Overall, a capacity retention of 90% with respect to the rst discharge cycle and a CE of ca.100% was achieved aer 40 cycles (Fig. 3c).For comparison, as-synthesised Na 2 Mn 3 O 7 was galvanostatically cycled using an identical protocol (i.e., 40 cycles at C/20 in the voltage range 1.5-4.7 V), showing a capacity retention of 59% with a CE of 96% (Fig. S15b †).These data agree with previous reports on Na 2 Mn 3 O 6.85 , where a capacity retention of 50% with a CE of 90% aer 50 cycles in the 1.5-4.5 V voltage range was observed. 59ate capability tests were carried out at different current densities from C/20 to 5C (i.e., 8 mA g À1 to 800 mA g À1 ) in the 1.5-4.7 V vs. Na + /Na potential window, acquiring data for ve consecutive cycles at each rate.Fig. 3d shows discharge capacities at different applied rates (solid squares), together with the CE (open circles).The rst discharge cycle shows a capacity of 205 mA h g À1 when applying a current of C/20 to the electrode.Increasing the applied current to C/10, C/5, C/2, 1C, and 5C, and back to C/20 leads to specic capacities (at the 5 th cycle) of 210, 200, 182, 171, 130 and 187 mA h g À1 , respectively.The discharge capacity of 130 mA h g À1 observed at 5C corresponds to 64% of the rst discharge capacity obtained at C/20.By contrast, at 5C the capacity for Na 2 Mn 3 O 7 was 51% of that obtained in the rst cycle at C/20 (Fig. S15d †).When the current returns to C/20, a reversible discharge capacity of 187 mA h g À1 was achieved.This corresponds to a capacity retention of 92% with respect to the initial discharge cycle.As the charge-discharge rate increased from C/5 to 5C, the voltage prole curves changed slightly, shortening the charge/discharge plateau at 4.0 V (Fig. S16 †).Therefore, the kinetically limited process seems to be related to the high voltage oxygen redox activity, which has been described as a rate-determining step in the charge/ discharge process. 60Nyquist plots for Na 2.4 Al 0.4 Mn 2.6 O 2 and Na 2 Mn 3 O 7 showed larger charge-transfer resistance values for the former, becoming even more evident aer 100 cycles (Fig. S17 and S18 †).EIS data for other Al-substituted compounds, such as P2-Na 0.5 Mn 0.5Àx Al x Co 0.5 O 2 15 and Aldoped LiNi 0.5 Mn 1.5 O 4 61 have also shown an increase in charge transfer resistance with Al 3+ ion substitution, which has been attributed to the formation of a thicker cathode-electrolyte interface (CEI) layer at high states of charge and/or a more unstable surface. 15,61,62arge compensation mechanism Ex situ XRD data were collected on cycled electrodes at different states of charge during the rst cycle, to examine the structural arrangements occurring in Na 2.4 Al 0.4 Mn 2.6 O 2 (Fig. S19 †).XRD data on the electrode charged to 4.7 V (red data) shows broadened peaks and changes in the relative intensities, suggesting an increase in disorder and atomic occupations with respect to the pristine electrode.Furthermore, an increase in the cell volume parameters attributed to an increase of the interlayer distance was observed.4][65] No extra phases were observed, indicating a solidsolution behaviour.On the other hand, XRD data on the electrode that was charged and discharged (4.7-1.5 V voltage range) shows two extra peaks at z38.5 and 42.8 2q (indicated with a diamond symbol in Fig. S19 †).These are in good agreement with the orthorhombic Na 0.44 MnO 2 phase (Pbma space group), previously observed during discharge in Na 2 Mn 3 O 7 (Fig. S20 †). 66able S2 † summarises the crystallographic data obtained by indexing the XRD data shown in Fig. S19 † using the Le Bail method.
Ex situ 23 Na solid-state NMR spectra for Na 2.4 Al 0.4 Mn 2.6 O 7 electrodes (pristine, charged and discharged) are shown in Fig. S21.† The spectrum for the pristine cathode is similar to that of the Na 2.4 Al 0.4 Mn 2.6 O 7 powder with distinct features at $500 ppm and $730 ppm, and with an added intensity around 0 ppm due to residual electrolyte.The resonance at 330 ppm is hidden by a spinning sideband.In the charged electrode material, the 500 and 730 ppm features are not observed, and the spectrum instead shows intensity at a higher shi of 910 ppm.Although this shi is lower than the 1500 ppm observed for the undoped material, the increase in shi compared to the pristine doped material is consistent with the removal of Na + ions from the structure, resulting in a more Mn 4+ ion rich local environment. 3,44The spectrum for the discharged material, whilst broadened, is similar to the pristine electrode material.For comparison, ex situ 23 Na SSNMR spectra were also obtained for Na 2 Mn 3 O 7 electrode materials (Fig. S22 †).As for the Na 2.4 Al 0.4 Mn 2.6 O 7 materials, there is an intense peak close to 0 ppm in the spectra, due to residual electrolyte and the diamagnetic component discussed previously.The pristine electrode material shows a main peak at $1500 ppm, consistent with the Na 2 Mn 3 O 7 powder.Spectra obtained at different MAS rates show that the shi of this peak may be temperature-dependent.The spectrum of the charged electrode material shows small added features at 930 ppm and 650 ppm (with one peak a spinning sideband of the other).The spectrum for the discharged electrode is different with the peak at $1500 ppm not observed and the spinning sideband pattern of the broader component of the spectrum shied to lower ppm and over a wider shi range.This is consistent with the discharged material having an increased sodium content, and therefore different local Na environments when compared to the pristine material.
Although new Na environments can be detected with NMR, it has not been possible to decide their nature.Therefore, to gain further understanding of the charge/discharge behaviour in Na 2.4 Al 0.4 Mn 2.6 O 7 and its Na 2 Mn 3 O 7 parent material, we examined via computational methods the insertion sequences of Na atoms and relate the results to the voltage proles of the cathode materials upon cycling.
For this purpose, the Na sites of Na 2 Mn 3 O 7 and Na 2 Al 0.4 -Mn 2.6 O 7 occupied during charging/discharging were predicted by optimizing the geometries of partially sodiated Na x Mn 3 O 7 and Na x Al 0.4 Mn 2.6 O 7 supercells.Upon cycling of the sodium layered electrodes, the local bonding between Na atoms and TM changes largely, whereas TM-O bonds mostly maintain their octahedral coordination.For such systems, the electrochemical potentials are primarily determined by the bond dissociation energy between active Na atoms and TM layers.From this perspective, the difference in Na + ion insertion sequences between Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 and their differing Na bonding environment during sodiation can result in differing voltage proles.To verify this relationship, we rst summarized all local bonding environments of Na sites in Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 structures (Table 1) and assessed how strongly Na + ions are bounded in each coordination environment by comparing their Na vacancy formation energies (Fig. 4a).Comparison in Na vacancy formation energies reveals that, of two Na sites in Na 2 Mn 3 O 7 , Na + ions in site Na1 are more strongly bonded to TM layers, requiring higher potentials to remove them compared to those in site Na2.This matches the two voltage plateaus observed, where Na + ions from site Na2 are rst removed at low voltage (Na 2 Mn 3 O 7 / NaMn 3 O 7 ), followed by the removal of Na + ions from site Na1 at high voltage (NaMn 3 O 7 / Mn 3 O 7 ).
For Na 2.4 Al 0.4 Mn 2.6 O 7 , the presence of Al 3+ ions and Na + ion excess changes the structure such that Na + ions in Na1 and Na2 sites have three and two different variations of local bonding environment, respectively (represented by prime symbols in Table 1).Consequently, Na vacancy formation energies may vary even for the same Na site: 4.74-5.17eV for Na1 site and 4.31-4.81eV for Na2 site.This shows that Na + ion extraction from these sites during charging occurs throughout various voltage values, resulting in sloping voltage proles of Na x Al 0.4 Mn 2.6 O 7 in the range of 0.4 < x < 2.0 (Fig. 3b).The relatively high Na vacancy formation energies from Na1 and Na2 indicate that the removal of these Na sites will require high charge voltages (4.2-5 V), which corresponds to the high voltage plateau in Fig. 3b and redox peaks in region C in the CV data (Fig. 3a).On the other hand, the Na vacancy formation energies of sites Na3 and Na4 are lower than those of sites Na1 and Na2, which is due to electrostatic repulsions with adjacent Na atoms in Na1 and Na2 sites.Such weakly bonded Na atoms are likely to be removed rst during charging, which is represented as the low voltage anodic peaks in region B of the cyclic voltammetry data (Fig. 3a).Additional theoretical voltage calculations reveal that the model structures resemble well the rst charge curves of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 (Fig. 4b and c): (1) predicted voltage proles agree with the overall shape of the charging curves, including the two-stepwise proles of Na 2 Mn 3 O 7 and the sloping curves in the 0.4 # x # 2.0 region and the distinct low voltage plateau in the 2.0 # x # 2.4 region of Na x Al 0.4 Mn 2.6 O 7 .
(2) When experimental voltage values are averaged over the x range of each voltage plateau, the averaged voltage values are close to the predicted voltages by less than 0.5 V.However, the above predictions are no longer consistent with the experimental data as charging and discharging processes proceed, due to structural changes which are not considered in these calculations.
Furthermore, we simulated the rst desodiation process of Na 2 Mn 3 O 7 and Na  S3 †).Considering that the energy to generate a positive hole is lower for lower Madelung energy, 38,67 electron holes tend to be localized on 2p states of O near the unoccupied Mn sites.Such hole localizations can also be predicted from the partial density of states of O atoms (Fig. S24 68 which is in good agreement with previous reports on Na 2 Mn 3 O 7 . 8,59pon further desodiation, electron holes continuously accumulate in non-bonding O atoms and, in turn, induce some structural distortion in the lattice of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 -Mn 2.6 O 7 .To understand this coupling between excess holes and structural changes, we rst analysed the structural changes near the non-bonding O atoms.2), the threecoordinated O2 (O coordinated by two Mn and one Na) is the most preferential hole site because of its low Madelung   2 and Fig. 5e), where each oxygen atom has differing coordination environment; two Mn and Na atoms for O2 0 ; Mn, Al and Na atoms for O3; and Mn, Al, and two Na atoms for O5 0 .Madelung potentials of Al-coordinated O (17.37 and 20.85 V for O3 and O5, respectively) are lower than that of two-Mncoordinated O2 0 (20.86 V), indicating that, upon further charging, excess holes are likely to be localized in O atoms adjacent to Al.Additional calculations on the O vacancy formation energies show that among all O atoms nearby the unoccupied Mn sites (O1, O2, O2 0 , O3, O5 0 and O5), partially oxidized O atoms, especially those coordinated by Al (O3 and O5 0 ), are more readily removed upon charging.This suggests that owing to the preferential hole formation excess in Alcoordinated O atoms during charging, Na 2.4 Al 0.4 Mn 2.6 O 7 can also suffer from the formation of O-O dimers and associated O loss upon charging.
Calculations on the Na and O vacancies of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 enabled us to reveal the early stage of the anionic redox reactions during charge.Na + ion extraction during charging partially oxidizes non-bonding oxygen atoms adjacent to the unoccupied Mn sites, causing the oxygen atoms to have localized excess holes (Fig. 5a-d).The distribution of excess holes can be modied by external dopants, such as Al, where holes are preferentially located on non-bonding oxygen atoms coordinated by Al.Next, the so-formed electron holes in oxygen atoms shorten the O-O bond distance up to 2.4 Å (Fig. 5e).Compared to other oxygen atoms, those comprising the shortened O-O bonds can be readily removed from their original sites (Table 2), which suggest the formation of oxygen dimers with shorter bond lengths, e.g., peroxo-like O-O dimers, and O 2 gas.Song et al. hypothesised on the possibility of O 2 gas evolution from mostly surface but also bulk Na 2 Mn 3 O 7 particles together with a decrease of crystallinity and increase of particle stress, although without providing any further experimental evidence. 6he above ndings of contracted O-O bonds and their weak binding to the TM layer were observed in both Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 , suggesting that these are susceptible to O loss upon the rst charge.Such O loss can also be supported from experiments: both cathodes display shortened high voltage plateaus during the rst discharge process and subsequent charge-discharge cycles (Fig. 4b and c), which indirectly indicates a decrease in the amount of redox-active O. Furthermore, additional calculations on the electronic structures of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 show that the removal of nonbonding O atoms leads to the formation of localized excess electrons on two adjacent Mn 4+ ions (Fig. 6).The so-formed excess electrons reduce both Mn 4+ ions to Mn 3+ , where excess electrons positioned in Mn with lower Madelung potentials exhibit higher energy (Table S4 and Fig. S29-S31 †).This suggests that O loss at high voltages causes the reduction of Mn 4+ to Mn 3+ to maintain charge neutrality, which can explain the mixed oxidation states of Mn at high voltages observed in XANES spectra.Fig. S10a and b † illustrates XANES spectra of Mn K-edge at various states of charge, together with the MnO (Mn 2+ ), Mn 2 O 3 (Mn 3+ ) and MnO 2 (Mn 4+ ) references.Linear combination t analyses (LCF) were used to semi-quantify the oxidation state of Mn at the different points of charge and discharge and the results are shown in Fig. S10c.† A report of the goodness of t parameters (R-factor and reduced c 2 ) along with the per cent that each model contributes to the t is summarised in Table S5.† A similar reduction of the TM at highly charged states due to the oxygen redox activity has also been observed for Ni 4+/3+ in Li 1.2 Ni 0.4 Ru 0.4 O 2 . 71Considering that the Mn 3+ O 6 octahedra are Jahn-Teller active, 72 such reduction of Mn 4+ ions might have detrimental effects on the cycling stability of the studied cathodes.
Although the above analyses have focused on O loss and subsequent effect on Mn reduction, electron holes formed in O atoms upon desodiation can cause other structural changes.In addition to O 2 gas release, these excess holes can shorten the O-O bond lengths to form various O 2 dimer species, such as peroxo-like O 2 nÀ or peroxide O 2 2À -species. 10Previous studies on Li-ion layered cathodes have reported that the shortened O-O bond can trigger TM migration to octahedral sites in the Li layer and subsequent formation of spinel-like phases. 73Other studies have also suggested the possibility of cation dopant migration to the tetrahedral site of the Na layer in layered oxides, which prevents further TM migration and stabilizes cyclability of cathodes. 74

Conclusions
In summary, we have studied the effects of Al 3+ ion substitution and Na + ion excess in the Mn-rich, Na 2 Mn 3 O 7 cathode material with triclinic P 1 space group.We synthesised Na 2.4 Al 0.4 Mn 2.6 O 7 (Na 0.68 Mn 0.74 Al 0.11 [ ] 0.15 O 2 in conventional Na x MO 2 notation) as a reference compound by a ceramic method and its crystal structure was characterised with X-ray diffraction.Given the highly disordered nature of the cathode material, we used DFT calculations to provide an accurate crystallographic model of the structure.We then assessed the electrochemical properties of the Na 2.4 Al 0.4 Mn 2.6 O 7 cathode in Na half-cells using different electrochemical methods, including galvanostatic cycling, cyclic voltammetry, electrochemical impedance spectroscopy, and galvanostatic intermittent titration.Cyclic voltammetric data showed the presence of both Mn 3+ /Mn 4+ and O 2À /O nÀ redox-active species in the material during Na cycling, which were further conrmed using XANES spectroscopy and computational studies.On the other hand, galvanostatic measurements showed an initial specic capacity of 215 mA h g À1 in the 1.5-4.7 voltage range at C/20 and a capacity retention of 90% aer 40 cycles.Experimental load curves were compared to computationally generated ones to gain insight into the Na + ion insertion/extraction sequences generated upon cycling.Subsequent examinations on the modelled structures revealed that owing to the high electronegativity of Al 3+ compared to the Mn 4+ counterpart, Na 2.4 Al 0.4 Mn 2.6 O 7 allows Na + ions to diffuse throughout various directions, parallel to the TM layers.The shorter Al-O bond lengths and lower oxidation state of Al 3+ ions distort the Na + ion migration pathways such that these are shorter and have wider diffusion channels, which further improve Na + ion diffusion kinetics.These results suggest selection criteria of dopants for high-rate performance layered oxides cathodes, where these should have: (1) higher electronegativity, (2) smaller ionic radius, and (3) lower oxidation states than TM elements.Therefore, we expect that this concept can be extended for developing future fast charging layered oxide cathodes.
To understand the underlying chemistry under oxygen loss and capacity fading, we performed partial charge density calculations on electron holes of desodiated Na  10 and TM migration into the alkali metal layer, 73 these studies may guide future studies on developing cathode materials that are resistive to oxygen loss and structural degradation.In conjunction with future calculations on structural changes during anionic redox (e.g., O, cation, and dopant migration), these ndings may enable researchers to uncover the effect of dopants on anionic redox reactions.

Na 2 . 4
Mn 2.6 Al 0.4 O 7 and Na 2 Mn 3 O 7 materials were obtained at 9.4 T on a Bruker Avance III 400 MHz spectrometer working at a Larmor frequency of 105.8 MHz.Powdered samples were packed into 2.5 mm rotors, in an argon-lled MBraun glovebox ([H 2 O] and [O 2 ] # 0.1 ppm) and kept static or spun up to 30 kHz.
Na MAS NMR spectrum for Na 2.4 Al 0.4 Mn 2.6 O 7 (Fig. S7b †) shows a wide spectral envelope, consisting of broadened intensity centred around 1000 ppm, and spinning sideband manifolds centred around 500 and 0 ppm, the latter again showing the presence of a diamagnetic impurity.Comparison of spectra recorded at variable MAS rate (Fig. S8 †) shows that the broad sideband manifold has isotropic peaks at 730 ppm, 500 ppm and 330 ppm.Comparison of both static NMR spectra (Fig. S9 †) conrms the lower shis for Na 2.4 Al 0.4 Mn 2.6 O 7 and the presence of two main groups at $1000 and $500 ppm.The lower shis observed for Na 2.4 Al 0.4 Mn 2.6 O 7 indicate a weaker Fermi contact interaction and the isotropic peaks at 500 and 330 ppm are consistent with the shi range for Na + ions within a local Mn 3+ 2.4 Al 0.4 Mn 2.6 O 7 upon discharge were then predicted by modelling partially desodiated Na x Al 0.4 Mn 2.6 O 7 (see Computational details for the construction of model structures, Fig. S1 and S2 in ESI †).Fig. 1d displays the symmetrically inequivalent Na sites predicted from Na 2.4 Al 0.4 Mn 2.6 O 7 , represented in the unit cell and supercell

Fig. 1
Fig. 1 (a) XRD data of layered Na 2.4 Mn 2.6 Al 0.4 O 7 and corresponding Le Bail fit using P1space group.Experimental data are presented in red, calculated data are shown in black and the difference between these two is shown in blue.Tick marks indicate the position of the calculated Bragg reflections.Miller indices are shown for the most intense reflections.(b) FE-SEM image of Na 2.4 Mn 2.6 Al 0.4 O 7 .(c) EDS mapping images showing O, Na, Mn and Al elements.(d) Schematics of predicted Na sites in Na x Mn 3 O 7 and Na x Al 0.4 Mn 2.6 O 7 .Top views of the Na and TM layers of supercells used in this study are shown.Number of occupied Na sites of the optimized (e) Na x Mn 3 O 7 and (f) Na x Al 0.4 Mn 2.6 O 7 supercell structures, showing Na + ion insertion sequences during sodiation.

Fig. 2
Fig. 2 (a) Changes in number and type of Na sites of Na x Al 0.4 Mn 2.6 O 7 under differing Na content.Results are plotted only for Na x Al 0.4 Mn 2.6 O 7 as Na x Mn 2.6 O 7 shows the same trends.Planar view of (b) Na x Mn 2.6 O 7 and (c) Na x Al 0.4 Mn 2.6 O 7 with respect to Na content, showing the active Na sites and associated Na + ion diffusion pathways.(d) Schematics of suggested Na + ion diffusion pathways of Na x Al 0.4 Mn 2.6 O 7 when excess Na + ions are intercalated into Na 2 Al 0.4 Mn 2.6 O 7 .
Fig. 3 (a) CV curves of Na 2.4 Mn 2.6 Al 0.4 O 7 at a scan rate of 0.1 mV s À1 in the voltage range 1.5-4.7 V vs. Na + /Na.(b) Galvanostatic chargedischarge voltage profiles in the voltage range 1.5-4.7 V vs. Na + /Na at C/20 (8 mA g À1 ) for cycles 1 and 2. (c) Discharge capacity (solid squares) and coulombic efficiency (open circles) as a function of cycle number in the voltage range 1.5-4.7 V vs. Na + /Na at a current density of C/20 (8 mA g À1 ).(d) Rate capability tests in the voltage range of 1.5-4.7 V vs. Na + /Na using current densities ranging from C/20 to 5C.

2 . 4
Al 0.4 Mn 2.6 O 7 and analysed its effect on electronic and atomic structures.Before desodiation, the valence band maxima (VBM) of both materials are predominantly formed by O 2p states and the conduction band minima are occupied by p-d hybridized states of O and Mn atoms (Fig. S23 †).Density of states calculations show that, aer removing Na atoms, electron holes are likely to be centred on the 2p states of oxygen.Additional analysis of the local Madelung potentials of O revealed that O atoms from Mn-O-Mn or Mn-O-Al near the unoccupied Mn sites in the TM layers (also referred to as nonbonding oxygen atoms 8 ) have lower Madelung potentials (23.93-25.88V) than those of other O atoms (27.39-29.86V) (Table †), where high-energy electrons are concentrated on the non-bonding O atoms.As a result, the desodiation of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 causes the localized holes at 2p states of non-bonding O atoms, which can be observed aer removing a Na atom from Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 (Fig. 5a-d).Such hole localization behaviours occur for all Na vacancies in Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 , regardless of the type of Na sites (Fig. S25 and S26 †).The above results indicate that the desodiation of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 proceeds by redox reactions from 2p states of non-bonding O atoms, Fig. 5e shows the changes in bond length between non-bonding O atoms in Na 2 Mn 3 O 7 and Na 2.4 -Al 0.4 Mn 2.6 O 7 before and aer desodiation.In the case of Na 2 Mn 3 O 7 , all O-O bonds near the unoccupied Mn sites have

Fig. 4
Fig. 4 (a) Local bonding environment of Na sites and their predicted Na vacancy formation energies in Na x Mn 3 O 7 and Na x Al 0.4 Mn 2.6 O 7 .Theoretical and experimental voltage profiles of (b) Na x Mn 3 O 7 and (c) Na x Al 0.4 Mn 2.6 O 7 .Electrochemically active Na sites for each Na content range are highlighted in differing colours.

Fig. 5
Fig. 5 Projected density of states calculated for (a) Na 1.9 Mn 3 O 7 and (b) Na 2.3 Al 0.4 Mn 2.6 O 7 .Insets show the magnified density of states near hole levels.Partial charge densities (green isosurfaces) of (c) Na 1.9 Mn 3 O 7 and (d) Na 2.3 Al 0.4 Mn 2.6 O 7 calculated for the hole bands located z1.5 eV above the VBM.(e) Distributions of O-O bond lengths near the unoccupied Mn sites before and after the desodiation of Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 .Distribution of hole polarons after desodiation are shown as insets, with labels for the corresponding types of O listed in Table 2.In the insets, the lengths of contracted O-O bonds are indicated by numbers.

Fig. 6
Fig. 6 Projected density of states of Mn atoms before (Mn host ) and after (Mn1 and Mn2) the formation of O vacancies in (a) Na 2 Mn 3 O 7 and (b) Na 2.4 Al 0.4 Mn 2.6 O 7 .Insets show O vacancy positions, reduced Mn atoms, and magnified density of states near excess electrons.

2 . 4
Al 0.4 Mn 2.6 O 7 , and conrmed oxygen redox reactions to stem from 2p states of non-bonding O atoms, analogously to Na 2 Mn 3 O 7 .In conjunction with O vacancy formation energies, energy barrier calculations suggested that the O atoms with electron holes during desodiation form "peroxo-like" O-O dimer species by O-O bond shortening.Furthermore, we showed that the preferential positions for O-O dimer formation can be adjusted by Al 3+ ion substitution in the TM layer.Such behaviour arises from the change in electrostatic energy (Madelung potential) of O atoms neighbouring Al atoms.This suggests that one can adjust the O-O dimer distribution of desodiated cathodes by careful dopant selection, i.e., by using dopants with low oxidation states and large ionic radii that lower the Madelung potential of adjacent oxygen atoms, one can force O-O dimer species to form nearby dopants, and vice versa.Considering that O-O dimers are a major cause of O 2 gas release

Table 1
List of Na sites in Na 2 Mn 3 O 7 and Na 2.4 Al 0.4 Mn 2.6 O 7 , showing the neighbouring structures and coordination of each Na site.The formation energies of Na vacancies are also shown for reference.All structures and vacancy formation energies were calculated using HSE06 functional

Table 2
70st of O sites in NaMn 3 O 7 , Na 2 Mn 3 O 7 , Na 1.4 Al 0.4 Mn 2.6 O 7 , and Na 2.4 Al 0.4 Mn 2.6 O 7 , showing neighbouring atoms and oxidation states of O sites.Relative formation energies of O vacancies are shown for reference TablesS3 and S4†) and its proximity to Na vacancies.Such selective partial oxidization of O 2À to O nÀ causes the contraction of bond distance to occur locally near O2 atoms (Fig.5eand S27 †).Additional calculations showed that the partially oxidized O nÀ acted as main sites for the formation of O-O dimer species: (1) the vacancy formation energies of O sites showed that the partially oxidized O nÀ required the lowest energy to be removed compared to other O sites (Table2); and (2) nudged elastic band calculations showed that the formation of the O 2 dimer from the non-bonding O nÀ species was a thermodynamically driven process (with reaction enthalpy of À0.29 eV), with a low activation barrier of 0.78 eV (Fig.S28 †).Together with Mn ion migration to the Na layer, the kinetic barrier for O-O dimerization can be further lowered.69Theabovendings,i.e., shortened bond distances of O2-O2 bonds, their low vacancy formation energies, and low activation barriers for O 2 dimer formation, suggest the formation of "peroxo-like" O-O dimers (and associated O 2 gas release) aer the rst charge of Na 2 Mn 3 O 7 .The formation of peroxo-like dimers has been demonstrated with resonant inelastic X-ray scattering (RIXS) in other layered oxides such as P2-Na x Li y Mn 1Ày O 2 (y ¼ 0.2, 0.25) with non-bonding O atoms.70The partial oxidation of non-bonding O atoms and associated O-O bond shortening can also be observed in Na 2.4 Al 0.4 -Mn 2.6 O 7 (Fig.5e).Overall, O-O bonds near the unoccupied Mn sites are slightly distorted compared to Na 2 Mn 3 O 7 such that the O-O bond lengths vary from 2.78 to 2.84 Å. Aer desodiation, the O-O bonds of Na 1.4 Al 0.4 Mn 2.6 O 7 are shortened to the lowest length of 2.42 Å.Such contraction in O-O bonds arise mainly from their partially oxidized O sites, i.e., O2 0 , O3, and O5 0 (Table