Coupling between oxygen redox and cation migration explains unusual electrochemistry in lithium-rich layered oxides

Lithium-rich layered transition metal oxide positive electrodes offer access to anion redox at high potentials, thereby promising high energy densities for lithium-ion batteries. However, anion redox is also associated with several unfavorable electrochemical properties, such as open-circuit voltage hysteresis. Here we reveal that in Li1.17–xNi0.21Co0.08Mn0.54O2, these properties arise from a strong coupling between anion redox and cation migration. We combine various X-ray spectroscopic, microscopic, and structural probes to show that partially reversible transition metal migration decreases the potential of the bulk oxygen redox couple by > 1 V, leading to a reordering in the anionic and cationic redox potentials during cycling. First principles calculations show that this is due to the drastic change in the local oxygen coordination environments associated with the transition metal migration. We propose that this mechanism is involved in stabilizing the oxygen redox couple, which we observe spectroscopically to persist for 500 charge/discharge cycles.

A more efficient and sustainable energy infrastructure requires a significant improvement in the cost and energy density of Li-ion batteries 1,2 . Lithium-rich layered oxide electrodes (Li 1+x M 1−x O 2 , 0 < x ≤ 1/3) have attracted interest as they offer access to substantially higher energy densities than conventional layered oxide electrodes 3,4 owing to the presence of a reversible anionic redox couple in addition to the usual transition metal (TM) redox couples [5][6][7][8][9][10][11][12][13][14][15] . Understanding the mechanism of anion redox in these systems is necessary to mitigate their unfavorable electrochemical properties 16,17 , which include significant charge/discharge voltage hysteresis (several hundred millivolts) and long-term voltage fade, even when cycled at vanishingly low rates [18][19][20] .
Establishing the oxygen redox mechanism in these electrode materials, however, has proved challenging. Tarascon and coworkers have suggested that oxidation of oxygen generally results in the pairing of O ions, resulting in an effective 2O 2-/O 2 nredox couple 5,[9][10][11][12][13][14]21 , which is stabilized against evolution as oxygen gas in the presence of 4d and 5d TMs due to the increased TM-O hybridization and improved band alignment over 3d TMs 22 . Ceder and coworkers, meanwhile, have predicted that O-O dimerization should only occur in the presence of d 0 and d 10 cations due to the rotational freedom of the O 2p orbitals, with a localized hole O 2-/Omechanism prevailing elsewhere 23 . Bruce and coworkers recently reported evidence for the localized O 2-/Omechanism in Li-and Mn-rich Ni/Mn/Co layered oxides (LMR-NMC) based on O K edge XAS 6,7 , yet similar observations have led others to conclude a 2O 2-/O 2 nredox couple in the same system 12,13 . Thus, a clear consensus on which mechanism prevails in which materials is lacking.
A key factor giving rise to the confusion is the difficulty of probing anionic redox species with conventional spectroscopic techniques, given the heterogeneous nature of many Li-rich materials. For example, Li-rich surface and near-surface regions often exhibit oxygen evolution and reconstruction during delithiation [24][25][26][27][28] , yet surface-sensitive X-ray photoelectron spectroscopy (XPS) is often used to infer the nature of bulk anion redox 5,9,11,29 . Likewise, correlations between spatially averaged Xray absorption spectra (XAS) are often used to assess hybridization 6,7,10,[12][13][14][30][31][32] , which assumes that redox chemistry occurs uniformly. The relatively subtle evolution in the spectra during anion redox, combined with their sensitivity to probing depth (a few nm for quantitative XPS and total-electronyield (TEY) XAS) and detection mode (e.g., self-absorption distortions in fluorescence-yield (FY) XAS) has contributed to the conflicting proposed mechanisms for oxygen redox in certain materials and, consequently, explanations for its role in determining capacity and electrochemical stability 15,21,22,32,33 .
Furthermore, these mechanisms are typically discussed without considering the evolution of local and average structures during delithiation 6,7,22,32 . Indeed, X-ray diffraction (XRD) [34][35][36][37][38][39][40] , transmission electron microscopy (TEM) 29,37,41 , neutron diffraction 42 , magnetic susceptibility 37,38,41,43 , and computational [44][45][46][47] studies have confirmed that the bulk structure in many Li-rich oxides evolves substantially during delithiation, primarily due to TM migration into Li sites between the TM-O layers ( Supplementary Fig. 1). Tarascon and coworkers have furthermore shown that the electrochemical signature of oxygen redox strongly suggests a coupling with structural evolution 16 . Nonetheless, the mechanistic link between structure dynamics and anion redox chemistry and stability remains largely unexplored and unclear, due to the aforementioned characterization challenges. Indeed, this connection is expected to be significant, as Ceder and coworkers have recently predicted that local coordination environment plays a crucial role in determining the oxygen redox potential 23 . Thus, clarifying the nature of anion redox and its effect on electrochemical stability requires an approach that simultaneously probes the spatial distribution of anion redox chemistry and the evolution of local structure.
In this work, we address this challenge by combining scanning transmission X-ray microscopy and nanoscale XAS (STXM-XAS) with resonant inelastic X-ray scattering (RIXS) and various structural probes to investigate the anion redox mechanism in LMR-NMC. The combination of STXM-XAS and RIXS provides the nanoscale distribution and emission signature of the features observed through conventional XAS, offering unprecedented insights into the nature of the oxidized oxygen species. We reveal that LMR-NMC exhibits two distinct reversible redox mechanisms during the first charge: (1) depopulation of hybridized TM 3d -O 2p bands whose redox potential is relatively constant ("TM-O redox") below 4.50 V; (2) depopulation of states with predominantly O 2p character ("O redox") above 4.50 V whose redox potential is dynamic and strongly coupled to TM migration and an associated change in the electronic structure of the material. We unambiguously confirm that both occur throughout the bulk of the primary particles, while oxygen evolution occurs only at high voltage in near-surface regions. We show that O redox persists for hundreds of cycles in uncoated LMR-NMC, which runs counter to recent predictions made for 3d systems under the assumption of a non-evolving structure 21,22,33 . Ab initio calculations reveal that the coupling between TM migration and O redox arises due to the dramatic change in local O coordination environment, which shifts the depopulated O 2p states to higher energy and lowers the O redox potential relative to that of the hybridized TM-O redox by > 1 V. Thus the O redox chemistry in LMR-NMC cannot be understood as a static O 2-→ O -+ eredox couple, but rather as a dynamic structure-redox coupled process described by {O 2-+ TM} → {O -+ TM mig } + e -. We suggest that this previously unconsidered structure-redox coupling plays an important role in stabilizing anion redox in LMR-NMC. Our results further suggest that anion redox chemistry can be tuned through control of the crystal structure and resulting TM migration pathways, providing an alternative route to improve Li-rich materials without altering TM-O bond covalency through substitution with heavier 4d and 5d TMs. More generally, we reveal a dynamic switching in cation and anion redox potentials, demonstrating the crucial role of simultaneous local and electronic structure rearrangement in determining electrochemical functionality in nominally topotactic intercalation materials.

Results
First cycle structural and electrochemical behavior. Uncoated secondary particles of Li 1.17 Ni 0.21 Co 0.08 Mn 0.54 O 2 were synthesized, as described in the Methods section. Rietveld refinement of the average structure in the C2/m space group to the synchrotron XRD pattern of the pristine powder yielded lattice parameters and Li-TM mixing (~2%) ( Supplementary Fig. 2), consistent with previous literature 35,36,42 . We note that a two-phase model consisting of a trigonal R3m phase and a monoclinic C2/m phase can also be used to describe the pristine LMR-NMC structure with similar accuracy 42 . However, similar Li-TM mixing is typically observed in these two-phase refinements, and so to improve the reliability of this refined parameter we employ the single-phase model which has fewer refined variables. Implications of the intermediate range in-plane cation ordering on the migration phenomenon investigated here are discussed later. To avoid Li transport limitations, electrochemical measurements were performed at a 4 mA g −1 current density (~C/68 where 1 C indicates the rate it takes to charge or discharge the electrode in 1 h), unless otherwise noted. Figure 1a shows the first and second cycle voltage and dQ dV −1 profiles. The sloping region during the first charge is entirely reversible if the material is cycled below 4.35 V, i.e., the peaks in dQ dV −1 profile are symmetric between charge and discharge. Upon charging to 4.60 V, the LMR-NMC electrode loses the honeycomb-like in-plane TM ordering in the TM layer, which manifests as the disappearance of the 9-11˚(20-23˚in Cu Kα) superstructure peaks in the XRD pattern (Fig. 1b) 48 . Simultaneously, structure refinement shows a substantial increase in the fraction of TMs in the Li layer from 2.8 % to 9.0 % (Fig. 1c, also see Supplementary Table 2), consistent with previous reports (see Methods) 35,36 . Upon subsequent discharge, the TM occupancy in the Li layer (4.7 %) does not return to that of the pristine material, the superstructure is not recovered, and the voltage profile is likewise permanently lowered in subsequent cycling. We note that this intra-cycle partially reversible bulk TM migration, which has been observed previously [34][35][36]40,43 , is distinct from the permanent TM migration that occurs at the surface during the first cycle, which is known to be the result of oxygen evolution and densification 27,[49][50][51] , and from the longer-term irreversible TM migration associated with the formation of a spinel-like structure and voltage fade 52 . Reversible TM migration is consistent with the TM site preference varying with electron count and/or Li content 53 , such that on charge certain electronic and structural conditions are met that promote TM migration, while on discharge these conditions are quenched by lithiation and the original TM sites become favored again. The partial irreversibility (i.e., the resulting in-plane disorder and residual TM occupancy in the Li layer) suggests some "trapping" of TM ions as well as hysteresis in the migration pathways during charge and discharge.
Oxygen oxidation during first charge voltage plateau. To understand the unusual redox chemistry of LMR-NMC we first investigate gross changes in the electrode electronic structure by following the average spectral evolution during the first cycle. Xray transparent samples were prepared ex situ through sonication and dispersion of electrodes harvested at various voltages ( Fig. 2a) in an Ar glovebox and transported to the microscope with minimal air exposure. The spatially averaged transmission XAS spectra obtained through STXM ( Fig. 2c-f), which represent the true depth-averaged absorption spectra, show that when charging below the plateau (pristine, 'P' → 4.35 V, '1') we observe primarily an inversion in the Ni L 3 edge peak ratio and a shift of the Co L 3 peak by + 0.4 eV, correlated with the simultaneous growth of a peak at 528.5 eV in the O K pre-edge. These changes, supported by the TM K edge XAS ( Supplementary Fig. 3), are wellunderstood to reflect the depopulation of hybridized O 2p -TM 3d (antibonding) bands [54][55][56][57] , constituting standard hybridized "TM-O redox".
By contrast, when charging through the voltage plateau (4.35 V, '1' → 4.60 V, '2'), the TM L 3 -edge spectra change minimally: the Co and Mn L 3 edges exhibit some broadening, and there is a small increase of the low energy peak at the Ni L 3 edge, suggesting a surprising reduction of the Ni during charge 12,30,31 . We later return to this important observation. By far the most significant change during the voltage plateau is the growth of a sharp and intense peak at 530.8 eV in the O K-edge spectrum, which indicates that predominantly O 2p character states are depopulated during the 4.50 V plateau 56,58,59 . The differential spectra in Fig. 2c clearly show that this feature is distinct from the peak at 528.5 eV, confirming that the two redox mechanisms are fundamentally different. We therefore refer to the 530.8 eV peak being the signature of "O redox", the nature of which we now investigate. Synchrotron XRD patterns throughout the first cycle showing the loss of inplane superstructure peaks (red dashed box) after charging through the 4.50 V plateau. The in-plane structure is not recovered during discharge, indicating an irreversible rearrangement in the cation ordering. c Fraction of the TMs in the Li layer (TM Li ) throughout the first cycle, obtained through refinement of the crystal structure to the XRD patterns, showing the substantial increase in out-of-plane disorder during the plateau. The TM migration is hysteretic and the material never fully recovers the pristine structure. Error bars representing the refinement error are smaller than the data symbols and are therefore omitted NATURE COMMUNICATIONS | DOI: 10.1038/s41467-017-02041-x ARTICLE We note that this is the first report of such a significant new peak emerging at the O K edge during the voltage plateau. Previous XAS results 6,7,10,12-14,30-32 mainly utilized electron and fluorescence based detection modes (Fig. 2b), which are distorted by surface contributions, self-absorption, and peak broadening effects and show only a weak and broad feature at the same energy. As a result, a definitive assessment of the plateau redox mechanism has been difficult. Figure 2b shows that we can reproduce these previously reported FY-and TEY-XAS results on the same electrodes, confirming that the discrepancy arises from the detection modes employed. By directly measuring the absorption coefficient, transmission XAS yields the true absorption spectrum.
O redox is a bulk phenomenon. To understand the chemical nature of O redox, it is important to know whether it is a bulk or surface phenomenon and whether it is spatially correlated to the subtle changes in the average TM spectra, which could be occurring in different regions of the electrode particles. We therefore map the nanoscale distribution of the Ni L 3 -and O Kedge spectroscopy at the single primary particle length scale during the first-charge voltage plateau. In Fig. 3 we apply principal component analysis (PCA) and non-negative matrix factorization (NMF) to a series of STXM spectro-images taken at six intermediate points along the 4.50 V plateau to identify the endmember spectra that describe the spectroscopic changes. At both the O and Ni edges, two end-member spectra were identified ( Fig. 3a): the O spectra are essentially with and without the peak at 530.8 eV, while the Ni spectra resemble Ni 4+ and Ni 2+ 54 . Using these end-member spectra, we obtain the average fraction ( Although all pixels measured with STXM do contain some surface signal, the surface contribution decreases with particle thickness, being lowest in the center of the particles' 2D projections. Since the ellipsoidal LMR-NMC particles are several hundreds of nanometers in thickness in the center, contribution from the surface (a few nanometers) is small, and the signal is dominated by the bulk. For the same reason, near the edge of the particles, surface signal dominates. Thus, from the nanoscale maps of the O end-member spectra (Fig. 3c, top), it is immediately evident that O oxidation occurs throughout the bulk of the primary particles during the voltage plateau, establishing unambiguously that O redox is a bulk process in LMR-NMC. On the other hand, the Ni chemical maps (Fig. 3c, bottom) reveal that bulk Ni ions remain in a 4 + oxidation state during the plateau. The variation observed in the average spectrum in Fig. 2d is in fact due to Ni reduction that is confined mostly to the primary particle surfaces. Spectral linescans in Fig. 3c show that the new peak at 530.8 eV at the O K edge is actually suppressed near the particle surfaces where Ni is reduced. In other words, although Ni reduction and O oxidation occur simultaneously during the voltage plateau, . The FY-XRM spectrum was collected on the same sample as the STXM spectrum, confirming sample preparation is not the cause of the discrepancy. c-f Spatially averaged transmission soft XAS throughout the first cycle at the c O K, d Ni L 3 , e Co L 3 , and f Mn L 3 edges, obtained by summing the spectra of all the pixels containing active material in each STXM spectro-image. Solid traces indicate electrodes harvested during charge, while dashed traces indicate electrodes harvested during discharge. The red and black traces below each plot indicate the differential spectra between the points indicated they are spatially separated. Its prevalence near the surface suggests that Ni reduction is more likely due to a steep gradient of oxygen non-stoichiometry due to oxygen evolution and material densification or reaction with the electrolyte 6,7,25,27,50,60 . This is further confirmed by the redox behavior during the first discharge and second charge (Fig. 4a), which reveals that discharging to 2.00 V reversibly reduces Mn near the surface with no change in the bulk. This is consistent with an electrochemically active but oxygen deficient near-surface structure that approaches the bulk composition within tens of nanometer.
O redox is non-rigid and stable in LMR-NMC. The spatially resolved nano-spectroscopy afforded by STXM reveals minimal pixel-wise correlation between the TM and O spectroscopy in the bulk during the voltage plateau, seen quantitatively in Fig. 4b and also in Supplementary Fig. 5  Furthermore, such a mechanism was previously predicted to be unstable against oxygen evolution in 3d layered oxides 21,33 . We therefore investigate the stability of oxygen redox over extended cycling using optimized 18650 cells (Samsung, see Methods), which exhibit 94% retention of their low-rate (intrinsic) capacity after 500 cycles at a 1C/2C charge/discharge rate. Figure 5b plots the O RIXS maps for the 501st cycle in the charged and discharged states and shows that even after 500 cycles the oxygen redox feature still reversibly appears. Thus, a significant fraction of the oxygen anions remains redox active after hundreds of cycles. These particles were not protected against oxygen evolution through coating. This is the first spectroscopic evidence of long-term stability of anion redox in 3d Li-rich materials and further suggests that the unstable rigid O 2-/Omechanism is unlikely.
Electronic reshuffling inverts redox sequence on discharge. The dynamic reshuffling of the O redox states has a profound effect on the LMR-NMC electrochemistry. Figure 2 reveals that the spectra on discharge (4.60 V, '2' → 3.65 V, '3' → 2.00 V, '4') do not follow a simple reversal of the changes observed during charge.    Fig. 2b and Supplementary Fig. 4a. The straight lines are due to pixels with 100% or 0% of a given oxidation state, shifted by the mean of the sample to which those pixels correspond Most notably, whereas O oxidation took place exclusively above 4.35 V during charge (after the Ni and Co had been fully oxidized to~4 + ), during discharge the Ni and Co are almost entirely reduced first in the 4.60-3.65 V range ('2 → '3'), with most of the O reduction taking place later in the 3.65 -2.00 V range ('3' → '4'). This massive > 1 V shift in the O redox voltage is illustrated by the color-coded regions in the dQ dV −1 plot in Fig. 2a. Thus the sequence of the redox couples (i.e., the order in which electronic states are (de)populated) is inverted after the 4.50 V plateau, and the inverted sequence (relative to the first charge) persists into the subsequent cycles ( Supplementary  Fig. 6), consistent with the altered structure and electrochemistry after the first cycle (Fig. 1). Spectro-imaging confirms that this is a bulk phenomenon ( Supplementary Fig. 7). This agrees well with our conclusion from XAS that some of the O 2p states depopulated during the plateau are shifted from an energy below the TM 3d -O 2p antibonding states to one above them, further confirming the electronic reshuffling associated with O redox.
O redox voltage modulated by TM migration. Figure 6a reveals a clear correlation between the hysteresis in the O redox and the fraction of TMs in the Li layer as calculated from Rietveld refinement of the synchrotron XRD patterns, suggesting that the electron state reshuffling arises from a strong coupling between O redox and TM migration. Supplementary Fig. 8 (Fig. 6b) to explore the effect of TM migration on the O projected density of states (pDOS). We did not include Co in our calculations due to its low concentration in our experimentally studied material. Ceder and coworkers have recently shown that anion redox chemistry is most strongly affected by the anion nearest-neighbor coordination environment 23 . We therefore designed the supercells to contain both possible oxygen coordination environments in LMR-NMC: O (1) , coordinated to two TMs and exhibiting the Li-O-Li geometry 23 , and O (2) , coordinated to three TMs. While there is a debate as to the intermediate range clustering of the TMs in LMR-NMC (i.e., whether it is a nano-composite of monoclinic Li 2 MnO 3 and rhombohedral NMC domains 3 , or a uniform solid solution monoclinic structure 67 ), we note that in both cases there are only the two types of oxygen environment mentioned, and thus the effect of TM migration on anion redox in either case can be understood by observing the effect on O (1) and O (2) sites. We first show in  Supplementary Fig. 9. Detailed discussion of all pathways investigated in this study is included in Supplementary Methods. While our treatment here is by no means exhaustive, we find that TM migration in all cases results in a shift to higher energy of the O pDOS regardless of the pathway examined and the identity of the migrating TM, suggesting that this observation is general and largely independent of the migration mechanism. We propose that this is due to 44 . An oxidized oxygen that was initially bonded to two TMs (O (1) in Fig. 6b) can become singly coordinated, with its depopulated O 2p states raised by several eV into the TM redox band. An oxygen that was initially unoxidized and bonded to three TMs (O (2) ) can become doubly coordinated and oxidized, transferring its electrons to another O or a TM. K-edge XAS ( Supplementary Fig. 3, cf. ref. 40 ) and extended X-ray absorption fine structure (EXAFS) measurements ( Supplementary Fig. 10) indicate that Ni migrates the most during the voltage plateau, suggesting an important structural role of Ni in controlling the O redox chemistry. Indeed, the shift of the O 2p states appears to be greater when the Ni migrates and the O (1) remains bonded to a Mn ion. The altered coordination environment of the shifted O 2p states due to TM migration is consistent with the distinct RIXS signature, while the shifting to higher energy is consistent with the inverted bulk redox sequence observed by STXM-XAS as well as the strong correlation between oxygen oxidation state and TM migration observed through XRD.
While we only observe a relatively small increase of TM in the Li layer (3-9%) during the plateau, each migration can decoordinate up to five oxygen ions (assuming the TMs move into neighboring octahedral sites), meaning that up to~20% of the oxygen sites can be directly affected. We also observe smaller indirect shifts in the oxygen redox states over longer ranges due to the distortion of the oxygen sub-lattice as a result of migration, and thus an even greater fraction of oxygen sites can be affected. We also note that in-plane TM rearrangement, which we observe as the loss of the superstructure peaks during the first cycle in Fig. 1b but cannot quantify due to their weak intensity and asymmetric shape, can have the same effect of reducing the overall O coordination number if, for example, the TMs move into tetrahedral sites. Thus additional in-plane TM migration beyond the 3-9 % out-of-plane migration may contribute to shifting a larger fraction of the O redox to lower voltage. Error bars indicating fitting residual and refinement error for the O oxidation state and TM migration fraction, respectively, are smaller than the data symbols and are therefore not shown. b pDOS for the TMs and the two-coordinate (O (1) , blue) and three-coordinate (O (2) , black) oxygen environments in the pristine delithiated state (top), and after Mn (middle) and Ni (bottom) migration into octahedral sites in the Li layer. The integrated charge density for the lowest unoccupied states in the pristine delithiated structure (blue shaded area of the pristine DOS) is shown in the top right inset. The total TM pDOS is normalized by the number of TMs in the supercell for comparison with the single-site oxygen pDOS. Schematics of each supercell used to generate the pDOS are shown to the left, with the plotted O (1) and O (2) oxygen environments circled. c Schematic of the reorganization of the electronic structure due to TM migration

Discussion
We therefore conclude that TM migration is intrinsically coupled to the bulk O redox potential, necessarily causing the reshuffling of the relative O and TM-O redox potentials. Thus, the O redox mechanism in LMR-NMC cannot be described as a static O 2-/Oredox couple, but rather as a dynamic {O 2-+ TM} → {O -+ TM mig } + eprocess, where TM mig indicates a migrated TM. In the specific case where a nominally 4 + Ni ion migrates into an octahedral Li vacancy in the Li layer, this can be written using Kröger-Vink notation 70 as: Since many different migration pathways for the various TMs are possible, we represent the substantial change and distribution in the oxygen coordination environments with {O -+ TM mig }, which modulates the oxygen redox potential and shifts a large fraction to lower voltage (<3.65 V) after the first charge. Figure 6c shows a schematic representation of this redox process from an electronic structure perspective. Importantly, if the O redox voltage is a primary factor affecting its stability, as has been suggested 22 , then the shifting of a substantial fraction of O redox to lower voltage coupled with TM migration may stabilize LMR-NMC against oxygen evolution, offering a possible explanation for the long-term anion redox stability observed here. Thus, while the long-term, irreversible TM migration associated with the formation of spinel-like structures is known to be associated with voltage fade 29 , the intra-cycle, partially reversible bulk TM migration may in fact be beneficial for the long-term cycling of LMR-NMC. While the role of the structure-redox coupling in voltage fade requires further study, we tentatively suggest that strategies aimed at increasing the latter while minimizing the former (e.g., by screening cation substitutions) may be effective in improving the cyclability of LMR-NMC. Further work will seek to understand the relative kinetics of O redox and TM migration, identify the intermediate states, and establish whether, for example, O 2p states are first depopulated, promoting TM migration through (de)stabilization of specific structural configurations, as has been predicted computationally 69 .
We note that contraction of the O-O distance to form O 2 nspecies has also been proposed as a non-rigid redox mechanism in certain Li-rich materials 5,[9][10][11][12][13][14]21,22 . Indeed, similar XAS features at 530.8 eV have been previously assigned to such species due to the comparable energy of the 1 s → σ* transition of Li 2 O 2 [12][13][14]71 . However, this transition in Li 2 O 2 is typically broad 71 , contrasting with the sharp resonance observed in LMR-NMC, and we furthermore show in Fig. 7 that the calculated XAS signature for O-O dimers only resides at this energy when they have bond lengths similar to that of Li 2 O 2 (~1.55 Å). At the larger O-O separations that we observe in our own calculations (~1.9 Å) and at the even larger separations observed in other Li-rich materials (≥2.2 Å 5,9 ) the 1 s → σ* transition lies at much lower energy (<529 eV) and therefore cannot explain our spectroscopic observations. Furthermore, Raman 6,7 and computational studies 23 have shown the formation of O-O dimers with true~1.55 Å bond lengths to be unlikely in LMR-NMC.
In summary, we have comprehensively established the bulk anion redox mechanism in LMR-NMC. By applying STXM-XAS and RIXS to resolve the spatial distribution and emission signatures of conventional XAS features, we overcome the longestablished challenges of material heterogeneity and spectroscopic ambiguity to demonstrate the first robust quantitative detection and analysis of bulk anion redox in Li-rich layered oxides. We reveal the inherent unreliability of quantifying anion redox through measurement of the overall O K pre-edge intensity alone 6,7 , which we show to be affected by both TM-O and O redox and has furthermore long been known to be heavily influenced by changes in TM-O hybridization strength 72,73 . We reveal that, unlike the lower voltage TM-O redox, the bulk 4.50 V plateau redox mechanism cannot be understood as a static redox couple. We show unambiguously that in this regime a reshuffling of the electron states and redox potentials occurs after depopulation of predominantly O 2p character states in LMR-NMC, giving rise to unique spectroscopic signatures that are distinct from previously reported rigid-band anion redox 56 . By correlating these observations with average and local structure probes, and in conjunction with ab initio calculations, we show that this dynamic character arises from a coupling between the O 2-/Oredox partner and TM migration, which reorganizes the electronic structure and shifts the O 2p states to higher energy relative to the TM-O states. This modulates the anion redox potential, inverting the TM-O/O redox sequence after the first charge and coupling the oxygen oxidation state to the fraction of migrated TMs. We therefore rewrite the mechanism as {O 2-+ TM} → {O -+ TM mig } + e -, which holistically describes the spectroscopic, electrochemical, and structural properties of LMR-NMC. We propose that this structure-redox coupling is involved in stabilizing the oxygen redox couple, which we observe to persist for 500 cycles in uncoated LMR-NMC despite predictions of intrinsic instability under the assumption of static structure. We therefore reveal that anion redox chemistry in many Li-rich materials cannot be fully understood by study of only their as-synthesized structures. Likewise, variations in structural behavior must be considered when comparing anion redox chemistry and stability between materials. Our results suggest that it may be possible to tune the stability and voltage of anion redox through control of the TM migration pathways. Thus, we suggest a new strategy for designing Li-rich layered oxides with improved cycling performance whereby the oxygen redox chemistry is tuned through structural modifications rather than the more common covalency modifications, which typically require substitution with rare 4d and 5d elements. More broadly, we demonstrate that atomic and electronic structural evolution during (de)intercalation need to be considered when assessing anion and cation redox chemistry even in nominally topotactic intercalation electrodes. The voltage curves at 4 mA g −1 were then measured in order to minimize the contribution of any kinetic overpotential to the measured voltage. All electrochemical measurements were performed in a temperature controlled chamber (Espec) set to 30°C. Details on the electrochemical measurements in Fig. 5b and general ex situ sample preparation are given in the Supplementary Methods.
X-ray diffraction. X-ray diffraction was performed at beamline 7-2 at the Stanford Synchrotron Radiation Lightsource (SSRL) at 14 keV (0.8856 Å) beam energy. The size of the X-ray beam was 100 × 100 μm 2 , and the distance between sample and the Pilatus-100K detector was 1000 mm. Electrodes were transferred to a Kapton adhesive film to remove the aluminum current collector. A pellet of the pristine powder was also measured, which showed negligible difference compared to the Kapton film sample. All of the samples were measured using the theta-2theta geometry. Rietveld Refinement details are given in the Supplementary Methods.
X-ray absorption measurements. XAS measurements in the hard X-ray regime were performed in transmission mode on the intact harvested electrodes (sealed under argon in a polymer pouch) at beam line 4-1 at the Stanford Synchrotron Radiation Lightsource at SLAC National Accelerator Laboratory using a Si (220) double crystal monochromator detuned to 50-60% of its original intensity to eliminate high order harmonics. The spectra of Ni, Co, and Mn reference foils were used to calibrate the photon energy by setting the first crossing of the second derivative of the absorbance spectrum to be 8333 eV, 7709 eV, and 6539 eV, respectively. Three ion chambers were used in series to simultaneously measure I 0 , I sample , and I ref . Spectrum normalization and alignment was performed using the Athena software package 74 .
Scanning transmission X-ray microscopy. Harvested electrodes were sonicated in DMC under argon for 2 h at 45°C to dislodge the LMR-NMC material from the electrode and break up the secondary particle structure. TEY XAS was performed on sonicated and un-sonicated materials and confirms that the procedure does not damage or alter the material, even in the~2 nm near the surface probed by TEY XAS (Supplementary Fig. 11). The primary particle suspension was drop-cast onto copper TEM grids with a carbon film (Ted Pella). The grids were then loaded onto a sample holder which was sealed under argon and brought to the beam line. STXM was performed at beam line 11.0.2 at the Advanced Light Source (ALS) in Lawrence Berkeley National Laboratory 75 , with preliminary measurements conducted at beam line 5.3.2.2 76 . Imaging was performed with either a 25 nm or 50 nm zone plate. To acquire an image, an interferometer-controlled stage raster-scans the position of the sample while a point detector measures the transmitted beam intensity 77 . While the step size was varied (35-80 nm) between edges and samples, the dwell time for each pixel was typically 1 ms. We confirmed that the ex situ imaging experiment did not damage the samples by measuring the Ni L 3 edge before and after imaging at the other three edges and observing no change to the average spectrum. STXM spectro-images were aligned in the aXis2000 software package. Subsequent data processing, including principle component analysis, was performed in the MATLAB 2016b software package and is detailed in the Supplementary Methods.
Resonant inelastic X-ray scattering. Soft X-ray RIXS maps (and conventional TEY and TFY XAS spectra) were collected in the newly commissioned ultra-high efficiency iRIXS endstation at Beamline 8.0.1 of the ALS 65 . All the samples were processed in glove box with high purity Ar environment, and the mounted samples were sealed and transferred into the experimental vacuum chamber through a specially designed sample transfer kit 78 . The kit was directly pumped down without stopping the pump while opening to the vacuum chamber to avoid any air exposure. Radiation damage was monitored by multiple scans of XAS and was considered negligible for these samples. Electrode sample surfaces were mounted at 45°to the incident X-ray beam. The outgoing photon direction (momentum transfer) along the RIXS spectrograph is 90°. RIXS resolving power and other technical details of the VLS spectrograph used for this work can be find in our previous report 79 . Details on how the RIXS maps are generated are given in the Supplementary Methods.
First-principles calculations. Lattice relaxations were performed using the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation 80 with the rotationally invariant Dudarev method to correct for self-interaction 81 as implemented in the Vienna Ab initio Simulation Package (VASP) 82 . A U values of 3.9 and 6.0 eV were used for Mn and Ni, respectively, according to previous literature 23 . Projector-augmented wave pseudopotentials were used 83,84 with a planewave basis set with a kinetic energy cut-off of 600 eV. Generally, all atomic positions and lattice parameters were allowed to relax until the forces on the atoms were less than 0.01 eV Å −1 . In certain cases, lattice parameters were fixed to prevent changes in stacking order during migration that otherwise proceed with no energetic barrier (discussed below). Since the structural composition is kept constant during migration, migration enthalpies are calculated as the difference in the cell energies before and after migration. Oxygen vacancy formation enthalpies are calculated with reference to molecular oxygen gas in a 10 Å × 8.5 Å × 10 Å cell. Ferromagnetic coupling was assumed. For electronic structure calculations, spin-polarized Density Funcitonal Theory (DFT) calculations employing HSE06 functionals were performed using the Quantum Espresso software package 85,86 . An exact exchange mixing parameter of 0.25 was used for all calculations. Fritz-Harber Institute pseudopotentials were used to replace the all-electron ion potentials, and the electronic wavefunctions were described using a plane-wave basis set with a kinetic energy cutoff of 80 Ry. A kpoint grid of 2 × 2 × 2 was used for typical supercell dimensions of~10 Å × 8.5 Å × 10 Å. The input structures were obtained from the PBE lattice relaxations, and the nuclear positions were allowed to further relax at the HSE level, keeping the lattice parameters fixed. Details on the methodology employed are discussed in Supplementary Note 2.
Data availability. The data supporting the findings of this study are available from the corresponding authors on request.