Local structure and lithium diffusion pathways in Li4Mn2O5 high

7 hours ago - ... structure can be achieved by a combination of Near-Edge X-ray Absorption Spectroscopy and Pair Distribution Function analysis of neu...
0 downloads 8 Views 2MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Local structure and lithium diffusion pathways in Li4Mn2O5 high capacity cathode probed by total scattering and XANES Maria Diaz-Lopez, Melanie Freire, Yves Joly, Claire V. Colin, Henry E Fischer, Nils Blanc, Nathalie Boudet, Valerie Pralong, and Pierre Bordet Chem. Mater., Just Accepted Manuscript • DOI: 10.1021/acs.chemmater.8b00827 • Publication Date (Web): 19 Apr 2018 Downloaded from http://pubs.acs.org on April 19, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Local structure and lithium diffusion pathways in Li4Mn2O5 high capacity cathode probed by total scattering and XANES Maria Diaz-Lopez1*, Melanie Freire2, Yves Joly1, Claire V. Colin1, Henry E. Fischer3, Nils Blanc1, Nathalie Boudet1, Valerie Pralong2 and Pierre Bordet1. 1 Université Grenoble Alpes, Institut Néel, F-38000 Grenoble, France CNRS, Institut Néel, F-38000 Grenoble, France. 2 Laboratoire de Cristallographie et Sciences des Matériaux CRISMAT, ENSICAEN, Normandie Université, CNRS, Université de Caen Normandie, 6 Bd Maréchal Juin, F-14050 Caen, France. 3 Institut Laue-Langevin, 71 avenue des Martyrs, CS 20156, 38042 Grenoble cedex 9, France. *e-mail : [email protected]

Abstract In the constant race for more efficient Li-ion batteries, extensive research has focused on the design of new, more competitive cathode materials, currently limiting the battery performance. The improvement of cathode materials demands the detailed understanding of the complex structural mechanisms at play during battery operation, i.e. when Li+ ions are inserted and extracted from the cathode. Moreover, new cathode designs involve more and more disordered/nano-sized materials for enhanced Li+ cation diffusion and larger specific surfaces. This trend poses new challenges for the structural investigation methods employed, which mostly rely on the periodic and long range ordered nature of the compounds under study. This is specially the case of the recently discovered nanostructured Li4Mn2O5 high capacity cathode material, which shows record reversible capacities superior to the state-of-the-art Li-Mn-O electrodes and displays a strongly disordered rock salt-type structure. This last feature, mainly due to its synthetic route involving high energy milling, prevented from reaching a full understanding of the lithium exchange mechanism of 1 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

particular interest in this 3D framework compound. Here, we demonstrate that a thorough description of such a disordered structure can be achieved by a combination of Near-Edge X-ray Absorption Spectroscopy and Pair Distribution Function analysis of neutron and X-ray total scattering data, which ultimately lead to the elucidation of the Li cation diffusion pathways.

Text The raise of lithium-ion rechargeable battery technology as a promising candidate to meet the current demands for clean and renewable energy sources1 has led to an extensive research on the design of new, more competitive materials with improved performance2. Amongst the several cathode materials explored, most efforts were focused on transition metal oxides intercalation compounds with 2D structures having acceptable capacities (100-250 mAhg-1) and operating voltages (2.5-5 V) 3. Recently, the focus has changed to layered Li-excess oxides with even higher capacities (250-300 mAhg-1)4–7. Cells with cathode materials based on spinel 3D structures also have adequate performances8–10 (e.g. LixMn2O4 has a reversible capacity of 120 mAhg-1 over extensive cycling)11,12, which proves that Li-conduction Is not exclusive to layered materials. The recent discovery of Li1.25Nb0.25M0.5O2 with M = V13 and Mn14, Li1.9Mn0.95O2.05F0.9515 and Li4Mn2O516 rock-salt-type compounds with high capacities constitute new additions to the less studied family of cathode materials with 3D frameworks. Although reported to be efficient, the migration of Li in the non-stoichiometric Li4Mn2O5 (Li4) rock-salt structure is not yet understood, partly due to the lack of a detailed knowledge of this disordered structure at the atomic level. In this article, we present a comprehensive study of Li4’s structure to unravel the mechanism by which lithium is able to diffuse in this new and intriguing compound. Previous Transmission Electron Microscopy (TEM) and combined neutron and X-ray powder diffraction studies of Li4-xMn2O5 (x=0 and 4)16 revealed this material shows on average an MnO-type rock-salt structure, albeit highly distorted at the local scale due to the substitution of 2/3 manganese for lithium and the presence of 1/6 oxygen vacancies. In fact, the degree of disorder is such that the PDF data up to 12 Å could no longer be described with the rock-salt model17. The full structural description of Li4 is greatly challenged by its high degree of complexity involving the presence of oxygen vacancies, substitutional disorder on the cation site and nanostructuration. In this work, we approach the structural study of Li4 by the combination of Near-Edge X-ray Absorption Spectroscopy (XANES) at the Mn K-edge and combined neutron and X-ray total scattering measurements. First, the coordination environment of manganese was extracted from XANES, and this information was used to build the starting model for the Pair Distribution Function (PDF) analysis of the total scattering data. Then, 2 ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

the PDF analysis was carried out using the RMC modelling method, capable of generating a high degree of local disorder from highly ordered starting configurations18–21. Several refinements were performed on a ‘large-box’ with dimensions similar to the structural coherence length of Li417, simultaneously fitting both neutron and X-ray reduced pair distribution G(r) and total scattering S(Q) functions. The RMC-refined model was finally validated by the comparison of simulated and experimental XANES spectra. Additionally, a quantitative analysis of the RMC outputs was performed to understand the distribution and local environments of the Mn-framework and the mobile Li ions in the validated model. Finally, we introduce the Bond Valence Site Energy (BVSE) analysis of the RMC-refined model, where accessible Li+ sites were identified by the calculation of the mismatch of the Bond Valence Sum22 (BVS) with respect to Li+ formal valence23,24 and converted to energy units25. The BVSE map of Li revealed a plausible 3D Li diffusion pathway connecting large volumes of 5-coordinated lithium clustered around oxygen vacancies.

Results The starting model In order to propose a sensible initial configuration for RMC refinement, all the available structural information on the material was gathered and thoroughly analysed. This included real and reciprocal space Rietveld refinements17 and XANES data at the Mn K-edge. Following the subtraction of the Li2O impurity from the diffraction pattern, the composition of the main phase was also corrected. To this end, the ICP results referring to the overall composition including the main and impurity phases were inputted into the Rietveld refinement. The composition of the main phase was then determined to be Li3.6(2)Mn2.4(2)O5.4(2). The initial model was built based on this corrected composition, although we continue to address Li3.6Mn2.4O5.4 as ‘Li4’, as opposed to ‘Li3.6’, for clarity. The Li3.6-xMn2.4O5.4 x-PDF data showed how Li removal is accompanied by an isotropic contraction (breathing) of the cubic Mn-framework. Surprisingly, Li extraction can be carried out without the occurrence of major structural rearrangements as evidenced by the constant value of the refined coherence length. However, no detailed local structural information, particularly on the atom coordination environments could be extracted from x-PDF data alone due to the overlapping of all cation-cation, oxygen-oxygen contributions in this highly disordered rock-salt structure and the limited sensitivity to Li scatterers. Thus, we included neutron total scattering data, with a higher sensitivity to presence of the lighter Li and O atoms, in the PDF analysis of Li4. Additionally, we have probed the coordination environment of Mn by the analysis of XANES data performed at the Mn K-edge on Li4 and delithiated Li017. 3 ACS Paragon Plus Environment

Chemistry of Materials

The mere visual inspection of the XANES data (see Figure 1) reveals a shift to higher energies of the Li0 spectra due to an increase in the oxidation state from Mn3+ in Li4 to Mn~4+ 16 (see Figure S1). In both Li4 and Li0 spectra very similar pre-edge and overall spectra shapes are observed, albeit the expansion of the Li0 spectra in the energy scale (i.e. the separation between Li4 and Li0 peak maxima increases with the increasing energy). These observations suggest there is a high resemblance between the local structures of the two compositions and a contraction of the lattice for Li0 according to Natoli’s rule26, which is in good agreement with the smaller cell parameters of Li0 determined by real and reciprocal space Rietveld refinements17.

Normalized absorption (a.u.)

(a)

1.6

Li0 Li4 1.2

0.8

0.08

0.4

0.04

0.0

(b) 1.6 Normalized absorption (a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

6540

6540 6560 6580

6544 6600 6620

6640

MnO 1.2

0.8 0.08

0.4 0.04 6540

0.0 6550

6600

6544

6650

6700

6750

6800

E(eV) Figure 1. (a) Mn K-edge XANES of Li4 (green) and Li0 (purple). (b) Data expanded in the energy scale showing the flattening of the spectral features of Li4 and Li0 with respect to highly ordered MnO (black) denoting the presence of structural disorder. The insets in both figures are the magnified pre-edge regions. 4 ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Since the average structure of Li4-xMn2O5 can be described with the MnO rock-salt structure, its spectra is also introduced in Figure 1 for comparison. The significant differences between the spectra denote that the local structure is different to that of a rock-salt, in good agreement with the high mismatch observed between the PDF data and the average model up to ~12 Å17. Moreover, the high degree of local disorder in this cathode material is further supported by the flattening of the spectral features in Li4 and Li0 with respect to MnO (see Figure 1.b). Although small in amplitude, the pre-edge signal carries valuable information regarding the coordination environment of the absorber atom, i.e. Mn. For the K-edge of the first row transition metals with an open 3d shell, the 1s to 3d transition is responsible for the pre-edge signal (see Figure S2). Since this transition is forbidden by dipole selection rules its intensity is small for centrosymmetric environments, and strongly increases in intensity as progressively moving to non-centrosymmetric environments due to the hybridization between the 3d and 4p orbitals. Therefore, the intensity of the pre-edge for an octahedron is smaller than for a square pyramid and becomes rather large for a tetrahedral environment. See for example, the difference in intensities of the pre-edge in Mn2+O and in Ba3Mn5+2O8 with a relative intensity of 7 and 50 % respectively in Figure S3. The small intensity of the pre-edge signals in Li4 and Li0 confirmed the absence of tetrahedrally coordinated Mn. The presence of Mn in a tetrahedral coordination would entail rearrangements of the structure which would be detrimental for electrochemical cyclability and were neither observed on the short or long r-range PDF data, nor were they supported by the similar coherence lengths of Li4 and Li0. Alternatively, we considered the possibility of having Mn in octahedral, square pyramid and square planar environments, which may be expected in the presence of oxygen vacancies. Several models were built starting from a 2×2×2 supercell of the rock-salt type MnO cell, which is the smallest cell allowing to accommodate the Li3.6Mn2.4O5.4 stoichiometry of Li4, and their XANES spectra were simulated (see Figure 2). These model supercells include Li/Mn substitutional disorder on the cation sites and oxygen vacancy disorder on the anion sites without any atomic displacement from the average MnO structure. The best agreement with the data was obtained with the model for which Mn was found exclusively in an octahedral environment. The distance between the two pre-edge peaks, ~2 eV apart, is also in good agreement with the typical energy gap for octahedral configurations. The broadening of the preedge peaks was attributed to the distribution in energy of the t2g and eg states given by the 8 non-equivalent Mn sites in the initial model.

5 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2 (a) Pre-edge region of several simulations with Mn in different configurations: square pyramid (pink), square planar (blue), a 1:1 ratio mixture of octahedral and square pyramid (orange) and purely octahedral (black) compared to the Li4 data (green). (b) Partial Density of States (DOS) of the 8 nonequivalent Mn sites in the octahedral model showing a wide distribution of t2g and eg energy levels responsible for the broad pre-edge peaks. In both figures the energy of the Fermi level (EF) are denoted by vertical dashed lines. As a consequence of the octahedral coordination of manganese, lithium is clustered around the oxygen vacancies where it is at most 5-coordinated. In the starting model, the oxygen vacancies were distributed on the anion network so as to avoid the formation of large clusters of lithium. Such configuration is probable to allow the reversible lithium insertion/extraction for several cycles, as experimentally observed. The resulting initial model with the composition Li19Mn13O29, a = 8.332 Å with a Pm space group was expanded to form a large box for subsequent RMC refinement. This box with a = b = c = 58.324 Å (close to the coherence length

6 ACS Paragon Plus Environment

Page 6 of 28

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

of ~ 6 nm) contained a total of 20 923 atoms. It should be noted that the lattice parameters used here were determined by the combined Rietveld refinement of neutron and x-ray data with a very high accuracy.

The RMC Refinement In order to lead the RMC modelling algorithm into outputting physically and chemically sensible refined models minimum distance cut off constraints and soft BVS constraints were included in the refinement. The choice of cation-oxygen cut-off constraints was performed on the basis of the shortest distance values extracted from the G(r) function (i.e. 1.7 Å). For other atom pairs with overlapped signals, the choice of cutoff distances was performed on distances observed for related Li-Mn-O systems, subtracted by ~0.2 Å. BVS constraints were used as implemented in RMCProfile27 with R0 and b bond valence parameters for Mn3+, Li+ and O2- taken from 28. The refinement was further constrained by the simultaneous fitting of all the available data i.e. neutron and x-ray G(r) and S(Q), allowing for the short, medium and long range structures to be simultaneously modelled. Several tests were carried out before the final refinement to determine the correct weighting of each data sets. Once the refinement parameters were optimized and after ~40 ·106 atom movements were tested the overall goodness of fit of the refinement χ2 of ~8.5 did not decrease any further. Thereby refined models produced excellent fits to all data sets (see Figure 3).

7 ACS Paragon Plus Environment

Chemistry of Materials

1.5 1.0

Experiment RMC fit

Neutrons 0.5

1.0

0.0 0.5 -0.5 0.0

-1.0

S(Q)

-2

G(r) (Å )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 28

6

2.0 1.5

X-rays 3

1.0

0

0.5

-3

0.0

-6

0

5

10

15

20

-0.5

25

5

r(Å)

10

-1

Q(Å )

15

20

Figure 3. RMC fits of Li4 neutron (top) and X-ray (bottom) G(r) (left) and S(Q) (right) at 297 K.

Quantitative evaluation of the local structure In this section we introduce the quantitative analysis of the geometry of the constituent polyhedra of the outputted RMC refined models. The reported values account for the study of a total of ten refined models providing partial correlation functions with good statistical quality. It should be noted that each refinement produced consistent and reproducible results (to verify this, the interested reader is referred to Figure S4). We first evaluate the reliability of the refined models by computing the BVS of each atom (see Figure 4). A narrow distribution (~±0.5 v.u.) around the expected BVS values was observed for oxygen (1.9 v.u.), 6coordinated manganese (2.8 v.u.) and 5-coordinated lithium (1.0 v.u.). The BVS of the minority 6coordinated lithium (~5 %) of 1.1 v.u. was further away from its ideal value, yet reasonable. The calculated BVS plots as a function of the coordination number indicated the quasi-absence of 4-coordinated manganese and lithium, as expected by the absence of a peak in the PDF at these shorter cation-oxygen distances. The appearance of 4-coordinated sites and 5-coordinated manganese in the refined model, which are small in number and have broad BVS distributions, could be related to the thermal vibrations of the atoms around their ideal positions. Although the formation of such sites with a lower coordination number is captured in the RMC model, its presence is neither statistically meaningful nor relevant for the interpretation of the refined model. Therefore, a quantitative evaluation of only 6-coordinated manganese and 6 and 5coordinated lithium was carried out. 8 ACS Paragon Plus Environment

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Figure 4. BVS plots for (a) oxygen in grey and for (b) manganese (c) lithium in grey, orange and blue for 6, 5 and 4-coordinated sites respectively. The vertical red lines indicate the initial BVS values in the starting model. For the interpretation of the RMC outputted model the partial atomic pair distribution function gij(r)s is given in Figure 5.a, where we observe the contribution of the different pairs of atoms to the PDF. Here, we can clearly identify the sharpness of the Mn contributing peaks with respect to Li and O, denoting a higher degree of ordering within the cubic Mn-framework, corroborating our expectations17. This observation is in good agreement with the smaller shifts between the refined atomic positions and the initial positions within the rock-salt structure calculated for Mn with respect to the larger shifts observed for Li and O (see Figure 6.b).

9 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. (a) Partial atomic g(r)s of the Li4 RMC refined model at 297 K. (b) Lognormal distribution of shifts between the atomic positions in the refined and starting models. The calculated modes are 0.48, 0.41 and 0.26 Å and the mean values 0.67, 0.65 and 0.57 Å for O, Li and Mn respectively. The first two peaks in Figure 5.a centred at ~2 Å account to the nearest Li-O and Mn-O neighbours’ distances (see also Figure 6.a and b), where once more the sharpness of the Mn-O peak is in clear contrast with the broader Li-O peak. For Mn-O, a narrow distribution centred at 2.01 Å and a smaller peak at a longer 2.31 Å is observed. For Li-O, a broader and asymmetric peak with a maximum at 1.87 and a shoulder at 2.15 Å is observed for 5 and 6-coordinated sites (see Figure 6.b).

10 ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28

(a) 30000 Number of Mn

25000 20000 15000

4000 Coordination number 6 5 4

3000

4000 3000

5000

1000

0

0

(b) 30000

2000 1000 0

5000

25000 Number of Li

5000

2000

10000

800

4000

20000

600 3000

15000

400

2000

10000

200

1000

5000 0

0 0.0

0.2

Å

0 1.8 2.0 2.2 2.4 Averaged M-O bond distance ( )

Å

1.8 2.0 2.2 2.4 2.6 2.8 M-O bond distance ( ) Å

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

0.4

0.6

σ M-O ( )

Figure 6. Histograms of all Mn-O (a) and Li-O (b) distances in the refined model. From left to right: distribution of M-O distances, M-O distances averaged per M site and corresponding standard deviations (σMO)

given as a function of the coordination number.

We now proceed to evaluate the Mn/Li-O distribution of distances as a function of the coordination number. To quantify the degree of distortion of the octahedrally coordinated Mn and Li sites, we calculated the Minimum Bounding Ellipsoid (MBE)29 of every octahedron as follows. First, the ellipsoid principal axis lengths radii were calculated and ordered as R1≥R2≥R3. Then, the S parameter, where S = R3/R2 – R2/R130, was used to address the shape of the octahedra in the refined model as: axially compressed or prolate (S0) or regular (S=0). For a representation of the ellipsoidal fits of distorted octahedral the interested reader is referred to Figure S5. In Figure 7, the plot of S for Mn shows a narrow and symmetric distribution around S = 0, which indicates the longer Mn-O distances of ~2.31 Å are randomly distributed and do not give rise to prolate or oblates distortions in Li4. This observation was also supported by the smooth distribution of average Mn-O distances in Figure 6.a, and BVS in Figure 4.b. For Li, as well as for Mn, a narrow and symmetric distribution of the S value around zero is observed, which suggest that Ohcoordinated Li sites also tend to have spherical shapes.

11 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Histograms of calculated SMn (a) and SLi(b). The red vertical line indicates the initial value of S = 0 in the starting model The majority of Li (95 %) was 5-coordinated with square pyramidal geometry. To better understand the distortions in the 5-coordinated Li sites in the refined model, we evaluated the relaxation of Li and O atoms surrounding the oxygen vacancies (□), whose positions were fixed during the refinement. The distribution of refined □-O and □-Li distances is introduced in Figure 9, which shows that while the average refined values of 2.95 and 2.11 Å are identical/close to the initial values of 2.95 and 2.08 Å for □-O and □-Li respectively, the distribution of □-O distances is not symmetrical. Thus, from the inspection of the distribution of □-O in Figure 9, we can conclude that oxygen within the pyramids base are drawn in towards the vacancy as given by the shift of the maxima from 2.95 Å in the starting configuration to 2.81 Å. On the other hand, the refined □-Li distances displayed an almost symmetrical distribution around the initial value of 2.08 Å.

12 ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28

Number of -Li

(a)1500

1000

500

0 0.5

(b)

1.0

1.5

2.0

2.5

3.0

1500

Number of -O

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

1000 500 0 1.0

1.5

2.0

2.5

3.0

3.5

distance (Å) Figure 8. Histograms of □-Li (a) and □-O (b) distances. The red vertical lines indicate the initial values on the starting model. The small deviations of the □-Li distances from their initial values does not explain the splitting of Li-O distances in the refined model. To find the origin for the distribution of short and long Li-O distances centred at 1.87 and 2.15 Å in Figures 5.a and 6.b, we computed the projection of the Li shift along and perpendicular to the □-Li direction (See Figures 9.a and b). The projection of the refined lithium position along the □-Li direction shows a symmetrical distribution centred at zero, which is in contrast with the larger log-normally distributed shifts (mode = 0.25 Å, mean = 0.62 Å) perpendicular to this direction. In order to compensate its BVS, the 5-coordinatd Li was found to move towards some of the oxygen within the square pyramid’s base and away from others, resulting in the formation of largely distorted square pyramids with a wide distribution of Li-O distances. We quantified the direction of the observed shift perpendicular to the □-Li direction by projecting the refined lithium position in the pyramid’s base and computing the angle of the shift in Figure 9.c, where 0° accounts for displacement of Li along one oxygen atom and 45° is perpendicular to the pyramid’s base edge. Figure 9 shows a homogeneous distribution of Li over the surface of the pyramid’s base which points to a random displacement of Li in every direction.

13 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 9. The displacement of lithium was calculated with respect to its initial configuration at the centre of the base of the square-based pyramid of oxygen as illustrated in the top right figure. Red, shaded and full green spheres denote O, Liinitial and Lirefined respectively.

The displacement vector from Liinitial to Lirefined (in

black) was projected along and perpendicular to the vacancy-Liinit direction. (a) Projection of the refined lithium position along the □-Li direction, where positive values denote the displacement of lithium towards the oxygen vacancy, and (b) perpendicular to □-Li. (c) Histogram of calculated angles for the projection of the refined Li position in the pyramid’s base (where 0° accounts for a displacement along one oxygen atom and 45° is perpendicular to the pyramid’s base edge.

Confirmation of the validity of the RMC refined model In the XANES region the photoelectron has a low kinetic energy and is able to travel a long mean free path interacting with a large number of atoms. For this reason, the XANES data at the Mn K-edge are not only

14 ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

useful to determine the local symmetry and change in oxidation state of Mn from Li4, but can also serve us as a whole-structure probe. We have used approximately a quarter of the volume of the box for a single configuration probing a total of a thousand manganese sites to reproduce the experimental XANES spectra introduced in Figure 10 (see model A). The good agreement between the simulated and experimental XANES and the significant improvement of the simulation using the refined model with respect to the initial one (see Figure S6) are a robust proof of the validity of the RMC refined model.

Figure 10. (a) Experimental and simulated XANES spectra using the refinement strategies A-D. (b) 3 × difference calculated for the simulated spectra with respect to the experiment, where the vertical dashed line indicates the position of the white line. Note that the changes in the simulated spectra are more significant to the left of this line corresponding to the pre-edge region. The creation of a sensible starting configuration for RMC modelling was key to the success of the refinement. Although in principle, the RMC method is capable of fitting a given set of data from any starting configuration19; in practice, this approach could either be too time consuming to be feasible or, as in this case, it could be that the refinement of the diffraction data alone is not able to discriminate the correct 15 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 28

coordination environment of manganese. We introduce here the two refinements of the initial models in Figure 2.a with 50 % and 100 % of 5-coordinated manganese (namely B and C in Figure 10 and Table 1) using the same refinement approach as in model A; and a third refinement where we have used model A as the starting configuration, but where the swapping between lithium and manganese was allowed during the refinement (refinement approach D). From the results shown in table 1 we conclude that in this particular case, the RMC modelling approach can successfully generate disorder within the model to fit the data regardless of the coordination environment of manganese in the starting model. When the swapping between cations was allowed in refinement D, the refinement tends to homogeneously distribute Mn and Li amongst the 5-coordinated sites, also providing good fits to the data. Even though all the refinements in A-D produced good fits (notice the very similar overall χ2 in Table 1), the refined model produced in the refinement strategy A was the only approach which successfully described the XANES pre-edge and generated chemically plausible models (see refined BVS values in Table 1). Table 1 Goodness of fit and calculated BVS for the RMC refinement using different initial configurations.

Refinement Initial % of

Cations

ID

swapping

5-

χ2

Initial/Refined BVS Mn

Li

O

coordinated (Y/N) Mn

A

0

N

8.5(1) 2.5062/2.719(3) 0.9353/1.041(3) 1.7482/1.901(1)

B

50

N

8.6(1) 2.2813/2.609(3) 1.0230/1.078(2) 1.6929/1.875(1)

C

100

N

8.5(2) 2.1207/2.529(4) 1.0726/1.110(2) 1.6534/1.861(1)

D

0

Y

8.1(7) 2.5062/2.644(9) 0.9353/1.062(4) 1.7482/1.881(3)

3D Diffusion pathway of lithium The clustering of lithium in a square pyramid coordination around the oxygen vacancies together with the systematic displacement of these atoms away from the centre of the pyramid’s base lead to the formation of highly distorted □Li6 environments (i.e. the octahedron formed by the 6 Li atoms surrounding a vacant site). Moreover, the direction of the lithium shifts for given □Li6 environments seems random, as illustrated in Figure 9.c. Therefore, the better suited BVSE map calculation approach introduced in 2.4 was undertaken for the prediction the most plausible lithium diffusion pathways in Li4. This approach exploits the analogy between the squared BV mismatch and the Morse-type potential31, thereby transforming valence units into energy units. As for BV mismatch maps32, the graphical

16 ACS Paragon Plus Environment

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

representation of isosurfaces with constant BVSE can provide a clear image of the lithium diffusion pathways in the material studied33–35.

17 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

18 ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Figure 11. Graphical representations of isosurfaces with BVSE of 2.5 eV along a, b and c axes (from top to bottom respectively) in a Li4 4ax4ax4a supercell section of the RMC refined ‘big box’. Note that in the graphical representation of BVSE maps the energy threshold was lowered from 3 eV used in the calculations in Table 1 to illustrate the percolation more clearly. Green, purple and red spheres denote Li, Mn and O atoms respectively. The BVSE map in Figure 12 was performed in a 4ax4ax4a (where a is the lattice parameter of the ‘MnO’ cell) section of the RMC refined model introduced in 3.3 containing 488 atoms. The yellow isosurfaces represent regions in which Li+ could diffuse for a 2.5 eV energy threshold. This choice of energy threshold revealed an extensive 3D pathway network for Li+ transport where all the 5-coordinated lithium, constituting 95 % of the total lithium content, becomes accessible. The accessible volume elements belonging to the same ‘pathway cluster’ consisted of ‘…-Li-□-Li-Li-□-Li-…’ intercrossing chains, which percolate along the three axes. The pathway clusters shared common faces and edges with vacancy-lithium octahedra ‘□Li6’. The diffusion mechanism in Li4 can be better understood when looking at □Li6 environments individually. In Figure 13.b it can be clearly noticed how the Li+ isosurfaces describe a sphere around the oxygen atoms bridging the different Li in □Li6. The perspective side view of □Li6 in Figure 13.d shows how such spheres intersect in the empty tetrahedral interstices. Thus, the transport of Li+ between three lithium sites within the □Li6’s faces could occur via the momentary occupation of the neighbouring vacant tetrahedral interstices. The Li+ diffusion mechanism in Li4 bears a resemblance to that of other rock-salt type high capacity cathodes lacking the presence of oxygen vacancies (e.g. α-LiFeO2 or spinel-like LiCoO236,37). The crystal structure of these layered and spinel transitional metal oxides consist of a close-packed oxygen sublattice with octahedral and tetrahedral interstitial sites occupied by TM /Li with an ABC stacking sequence. In these structures , the lithium migration is reported to occur between two octahedral sites via a tetrahedral activated state38, mediated by the presence of cation vacancy clusters (divacancies for layered materials and di- and triple vacancies for spinels39). However, due to presence of oxygen vacancies in Li4, lithium migration rather occurs between two square-based pyramidal sites via a trigonal planar activated state, and contrarily to the aforementioned crystal structures, the presence of cation vacancies does not play a role in the lithium diffusion mechanism for this particular compound. The existence of pathway clusters alone is not sufficient to sustain Li+ migration throughout the material. Such clusters need to be arranged in a percolating network40,41 to enable long-range lithium transport. To fully address the lithium migration in Li4 we calculated the percolation energy, as the minimum energy above the ground state connecting the BVSE isosurfaces along a particular direction in the unit cell, and the accessible volume fraction (F), as the fraction of the unit cell that belongs to the infinite pathway33 for a 3 eV energy threshold. The promising percolation thresholds and F-values for Li4, comparable to those for 19 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 28

representative state-of-the-art cathode materials introduced in table 1, point out the high potential of Li4 warranting further investigations into this system.

Figure 12. (a, c) □Li6 connectivity and (b, d) superimposed BVSE map, Li-O connectivity and tetrahedral interstices (black spheres). The black dashed lines in d denote the intersection of three spheres centred at the neighbouring O. Table 2. Calculated volume fraction for ion mobility, minimum and percolation energies for several representative lithium transition metal oxides from models reported in the literature (column Ref.). The calculated BVSE maps are given in Figure S7.

Percolation E (eV) along

Ref.

Material

F (%)

Minimum E (eV)

a

b

c

Li4

14(3)

-5.3(4)

2.0(6)

2.0(6)

2.0(6)

This work

LiMn2O4

13.28

-4.4175

0.54

0.54

0.54

42

20 ACS Paragon Plus Environment

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

LiNi0.5Mn0.5O2

17.13

-3.1193

0.30

0.30

0.30

43

α-MnO2

19.44

-3.0518

No

No

0.18

44

LiMnBO3

4.5

-8.9629

No

No

1.81

45

LiCoO2

11.17

-5.1456

1.96

1.96

No

46

LiFePO4 olivine

9.86

-4.586

No

1.04

No

47

V2O5

10.48

-4.349

No

1.38

2.26

48

Li2FeSiO4

19.44

-3.8476

0.96

No

0.96

49

LiV3O8

11.67

-5.0514

No

1.34

2.46

50

Conclusion The complete understanding of the local disorder in nanostructured Li4 was reached via the RMC modelling of combined neutron and X-ray total scattering data. The success of the refinement relied on building a sensible starting model based on the knowledge of the average structure from diffraction experiments, the local environment of Mn from XANES and the implementation of an intelligent set of geometrical and chemical constraints. The analysis of the RMC refined model indicated that Mn is only present in an octahedral coordination that while distorted, is on average spherically shaped. Li is clustered around the oxygen vacancies and it is predominantly 5-coordinated with an unusual and highly distorted square pyramidal geometry. The 5coordinated Li is displaced away from the centre of the pyramid’s base leading to a wide distribution of Li-O distances. BVSE landscapes revealed large pathways for Li-diffusion surrounding all 5-coordinated Li sites which amounted to 95 % of the Li content constituting ~60 % of the cations in Li4. Moreover, these large isosurfaces percolated along the three axes providing Li4 with a unique 3D Li-ion diffusion network only ever reported for spinel electrodes.

Methods Synthesis. 5 g of nanostructured Li4 were synthesized by a two-step route, the first of which involved the synthesis of HT-LiMnO2 by the solid-state reaction method. Stoichiometric amounts of MnO, MnO2 and LiOH (with 5 wt% excess to compensate for evaporation at high temperatures) were mixed and heated at 1000 °C under an argon flow for 8 h. The second step involved the mechanical grinding of stoichiometric amounts of 21 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

HT-LiMnO2 and Li2O in a 20 ml capacity WC bowl with four WC grinding balls (with a 10 mm diameter) using a Fritsch Planetary Micro Mill PULVERISETTE 7 premium at 700 rpm for 10 h. Li4 was electrochemically oxidized versus lithium in a Swagelok type cell. After the first electrochemical charge to 4.8 V all lithium was subtracted from Li4. The battery was then paused and 20 mg of delithiated material (Li0) powder were recuperated and washed with acetonitrile for XAS experiment. The Li0 sample studied here corresponds to that in 17 , the Li:Mn ratio in ‘Li0’ is equal to 0.2±0.05:1 as determined by ICPOES. Near-edge X-ray Absorption Spectroscopy and FDMNES modelling. XANES data at the Mn K-edge was acquired at the French CRG-FAME beamline (BM30) at the European Synchrotron Radiation Facility (ESRF). The beam energy was selected using a cryo-cooled double crystal monochromator whose energy value was calibrated with a Mn-foil standard. A total of seven samples including five Mn standards in different oxidation states (Mn2+O, Mn3+2O3, Mn4+O2 and Ba3Mn5+O8), oxidized Li0 and pristine Li4 were diluted (~5 wt%) with boron nitride and measured in quartz capillaries. The capillaries for the air sensitive Li4 and Li0 were carefully prepared under inert conditions and epoxy sealed. The spectra were collected at room temperature using a 30-element Canberra Ge solid state detector. All the spectra herein introduced correspond to at least six scans that were averaged, background and self-absorption corrected and normalized with the ATHENA software51. Our data processing protocol is supported by the good agreement between our spectra of the standards (see Figure S3) and others found in the literature52. Experimental XANES spectra were simulated with the FDMNES software53 using the multiple-scattering theory under the muffin-tin approximation on the potential shape, i.e. the potential is spherically averaged in the atoms and constant between them. The code runs in direct space and the calculations were limited to a cluster embedded in a sphere centred on the absorbing atom. When several non-equivalent atoms were present in the unit cell, the corresponding number of simulations were performed and summed up to get the total spectra. We have first checked the convergence using increasing radius and found that for cluster size of 6.1 Å around the absorber atom corresponding to ~90 atoms, the spectra of the Mn3+2O3 standard was successfully reproduced (see Figure S8). Simultaneous RMC modelling of n-PDF and x-PDF. High-Energy Synchrotron X-ray powder diffraction (HESXRPD) patterns with high spatial resolution and counting statistics were acquired at the French CRG-D2AM beamline (BM02-ESRF) using λ = 0.4959 Å radiation for subsequent Rietveld and PDF analysis. The samples were prepared under inert conditions in epoxy-sealed borosilicate glass capillaries with an outside diameter of 0.7 mm. The data were collected at room temperature in the 0.8-128.0º 2θº range with 2 h of data

22 ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

collection time using a 2D pixel detector XPAD354. In order to collect this wide angular range, 120 diffraction images were taken at each degree 2θº. They were assembled, averaged and rebinned using the PyFAI55 software to produce a single 1D pattern. In addition, the diffraction pattern of an empty capillary was also measured. The data were calibrated using a LaB6 standard collected for 4 h. A Ni standard sample was also used to calibrate the experimental effects on the PDF data. The HE-SXRPD data were corrected and Fourier Transformed using PDFgetX3 software56 to yield PDFs with Qmax = 21.4 Å-1. Neutron powder diffraction patterns57 were acquired at the D4c diffractometer58 at the Institut LaueLangevin (ILL) with a monochromatized wavelength of 0.4990 Å. Approximately 1.5 g of powders were prepared in cylindrical vanadium cans of 6 mm of diameter sealed with Helicoflex in an inert Ar atmosphere. The data were collected at room temperature in the 1.7-135.0º 2θº range with 3 h of data collection time per sample. The data were corrected for background scattering, vanadium can and self-absorption, and normalized with the Correct software59. The corrected data was Fourier Transformed using the STOG software embedded in the RMCProfile suite 60 using Qmin = 0.4 Å-1, Qmax = 24 Å-1. The combined refinements of the neutron and x-ray G(r) and S(Q) data sets were carried out with the RMCProfile software60. To account for the finite size of the simulation box, S(Q) was convoluted with a box function. Note that although it is desirable to include the Bragg data in the refinement to further constrain the long-range average structure, this was not possible in the present study due to the current version of the software lacking the profile shape of constant wavelength neutrons and anisotropic size broadening. Due to the impossibility to perform a multiphase refinement with this piece of software, the contribution of a minor crystalline Li2O impurity phase was removed from the diffraction patterns as follows: the diffraction signals from the impurity phase were first well described by Rietveld refinements performed with the Fullprof software61 and then the simulated patterns using the refined parameters were subtracted from the full diffraction patterns (see Figure S9). For the interpretation of the RMC outputted models, a quantitative analysis of the refined atomic positions inside the large box was carried out with a specifically written code using MatLab62. Li+ diffusion in Li4. BVSE pathway models (also known as energy-scaled bond-valence mismatch landscapes or energy landscapes)34 were derived from the RMC refined static structure model for the elucidation of the Li migration pathways in Li4. BVSE maps were constructed by summing the BVSE contributions to a hypothetical Li+ ion at each point with a 7·10-4 Å3 volume up to a distance of 8 Å. SoftBV parameters24 were used to account for the influence of higher coordination shells (i.e. to avoid abrupt changes in the BVSE by cut-off effects at the boundaries of coordination shells), which could lead to the formation of artefacts in the BVSE maps.

23 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Additionally, the available volume fraction of a percolating pathway with a 3 eV threshold, and the percolation thresholds energy values along each axes were calculated from Monte-Carlo simulations using the BondStr software embedded in Fullprof-Suite61. The figures illustrating the local environments and diffusion pathways in this publication were produced with the VESTA software63.

Acknowledgements The authors acknowledge Dr Olivier Proux and Dr Jean-Louis Hazemann for their assistance during the experiment at FAME (CRG-ESRF, France). The refinements introduced in this publication were performed with the RMC-Profile software modified by Dr Wojciech Slawinski to allow for the simultaneous refinement of two S(Q) data sets, and his contribution is gratefully acknowledged. This work was supported by the ANR grant ANR15-CE05-0006-01 DAME.

Author contribution M.D., P.B. and C.V.C. were the experimentalists in the neutron and X-ray diffraction and XANES experiments, M.D. and P.B. performed the RMC refinement and model interpretation, M.D. and Y.J. performed the XANES simulations, M.F. and V.P. contributed to the synthesis of the material, H.F. lead the experiment at D4 and contributed to the data processing and correction for total scattering analysis, N.B. and N.B. designed the experimental set-up at D2AM, M.D and N.Boudet contributed to the data processing of X-ray data. M.D, C.V.C and P.B. contributed to writing the paper. P.B. and V.P. conceived and designed the project.

Competing interest The authors declare no competing financial interests.

References 1. Larcher, D.; Tarascon, J.-M. Towards greener and more sustainable batteries for electrical energy storage. Nat. Chem. 2014, 7, 19–29.

24 ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

2.

Armand, M.; Tarascon, J.-M. Building better batteries. Nature. 2008, 451, 652–657.

3. Ohzuku, T.; Makimura, Y. Layered Lithium Insertion Material of LiCo1/3Ni1/3Mn1/3O2 for Lithium-Ion Batteries. Chem. Lett. 2001, 642-643. 4. Armstrong, A. R.; Holzapfel, M.; Novák, P.; Johnson, C. S.; Kang, S.-H.; Thackeray, M. M.; Bruce, P. G. Demonstrating oxygen loss and associated structural reorganization in the lithium battery cathode Li[Ni0.2Li0.2Mn 0.6]O2. J. Am. Chem. Soc. 2006, 128, 8694–8698. 5. Thackeray, M. M.; Kang, S.-H.; Johnson, C. S.; Vaughey, J. T.; Benedek, R.; Hackney, S. A. Li2MnO3stabilized LiMO2 (M = Mn, Ni, Co) electrodes for lithium-ion batteries. J. Mater. Chem. 2007, 17, 3112-3125. 6. Johnson, C. S.; Li, N.; Lefief, C.; Vaughey, J. T.; Thackeray, M. M. Synthesis, Characterization and Electrochemistry of Lithium. Chem. Mater. 2008, 20, 6095–6106. 7. Sathiya, M.; Rousse, G.; Ramesha, K.; Laisa, C. P.; Vezin, H. ; Sougrati, M. T. ; Doublet, M.-L. ; Foix, D. ; Gonbeau, D. ; Walker, W. ; Prakash, A. S.; Ben Hassine, M. ; Dupont, L. ; Tarascon, J.-M. Reversible anionic redox chemistry in high-capacity layered-oxide electrodes. Nat. Mater. 2013, 12, 827–835. 8. Kawai, H.; Nagata, M.; Tabuchi, M.; Tukamoto, H. & West, A. R. Novel 5 V Spinel Cathode Li2FeMn3 O8 for Lithium Ion Batteries. Chem. Mater. 1998, 10, 3266–3268. 9. Fey, G. T.-K.; Wu, L.; Dhan, R. LiNiVO4: A 4.8 Volt Electrode Material for Lithium Cells. J. Electrochem. Soc. 1994, 141, 2279-2282. 10. Sigala, C.; Guyomard, D.; Verbaere, A.; Piffard, Y.; Tournoux, M. Positive Electrode Materials with High Operating Voltage for Lithium Batteries: LiCryMn2-yO4. Solid State Ionics. 1995, 81, 167–170. 11. Thackeray, M. M.; David, W. I. F.; Bruce, P. G.; Goodenough, J. B. Lithium insertion into manganese spinels. Mater. Res. Bull. 1983, 18, 461–472. 12. Amatucci, G.; Tarascon, J.-M. Optimization of Insertion Compounds Such as LiMn2O4 for Li-Ion Batteries. J. Electrochem. Soc. 2002, 149, K31–K46. 13. Nakajima, M.; Yabuuchi, N. Lithium-Excess Cation-Disordered Rocksalt-Type Oxide with Nanoscale Phase Segregation: Li1.25Nb0.25V0.5O2. Chem. Mater. 2017, 29, 6927–6935. 14. Wang, R.; Li, X.; Liu, L.; Lee, J.; Seo, D.-H.; Bo, S.-H.; Urban, A.; Ceder, G. A disordered rock-salt Liexcess cathode material with high capacity and substantial oxygen redox activity: Li1.25Nb0.25Mn0.5O2. Electrochem. commun. 2015, 60, 70–73. 15. House, R. A.; Jin, L.; Maitra, U.; Tsuruta, K.; Somerville, J.; Förstermann, D.; Massel, F.; Duda, L.; Roberts, M. R.; Bruce, P. G. Lithium manganese oxyfluoride as a new cathode material exhibiting oxygen redox. Energy Environ. Sci. 2018, 1-7. 16. Freire, M.; Kosova, N. V.; Jordy, C. ; Chateigner, D.; Lebedev, O. I.; Maignan, A.; Pralong, V. A new active Li–Mn–O compound for high energy density Li-ion batteries. Nat. Mater. 2015, 15, 1–6. 17. Freire, M.; Diaz-Lopez, M.; Bordet, P.; Colin, C. V.; Lebedev, O. I.; Kosova, N. V.; Jordy, C.; Chateigner, D.; Chuvilin, A. L.; Maignan, A.; Pralong, V. Investigation of the exceptional charge performance of the 0.93Li 4−x Mn 2 O 5 –0.07Li 2 O composite cathode for Li-ion batteries. J. Mater. Chem. A. 2018, 6, 5156–5165. 18. McGreevy, R. L.; Pusztai, L. Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures. Molecular Simulation. 1988, 1, 359–367.

25 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

19.

McGreevy, R. L. Reverse Monte Carlo modelling. J. Phys. Condens. Matter. 2001, 13, R877–R913.

20. Mcgreevy, R. L.; Zetterstrom, P. To RMC or not to RMC ? The use of reverse Monte Carlo modelling. Solid State Mater. Sci. 2003, 7, 41–47. 21. Keen, D. A.; Tucker, M. G.; Dove, M. T. Reverse Monte Carlo modelling of crystalline disorder. J. Phys. Condens. Matter. 2005, 17, S15–S22. 22. Brown, I. D. The chemical bond in inorganic chemistry: the bond valence model. Oxford University Press, Oxford. 2002. 23. Adams, S. Relationship between bond valence and bond softness of alkali halides and chalcogenides. Acta Cryst. 2001, 57, 278–287. 24.

Adams, S. softBVS web page, 2003. Available at: http://kristall.uni-mki.gwdg.de/softBV/.

25. Adams, S. From bond valence maps to energy landscapes for mobile ions in ion-conducting solids. Solid State Ion. 2006, 177, 1625–1630. 26. Bianconi, A.; Dell’Ariccia, M.; Gargano, A.; Natoli, C. R. Bond length determination using XANES. EXAFS and Near Edge Structure, Springer Series in Chemical Physics, Berlin. 1983, 27, 57–61. 27. Norberg, S. T.; Tucker, G.; Hull, S. Bond valence sum: a new soft chemical constraint for RMCProfile research papers. J. Appl. Crystallogr. 2009, 8898, 179–184. 28.

Brown, I. D. bvparm2016.cif. (2016).

29. Todd, M. J.; Yıldırım, E. A. On Khachiyan’s Algorithm for the Computation of Minimum Volume Enclosing Ellipsoids. Oper. Res. 2005, 1–18. 30.

Cumby, J.; Attfield, J. P. Ellipsoidal analysis of coordination polyhedra. Nat. Commun. 2017, 8, 1-8.

31. Adams, S. Practical Considerations in Determining Bond Valence Parameters. Brown, I.; Poeppelmeier, K. (eds). Bond Valences. Structure and Bonding, vol 158. Springer, Berlin, Heidelberg. 2013. 32. Adams, S. Bond valence analysis of structure – property relationships in solid electrolytes. J. Power Sources. 2006, 159, 200–204. 33. Adams, S.; Rao, P. Understanding Ionic Conduction and Energy Storage Materials with Bond-Valence Based Methods. Struct. Bond. 2006, 119, 193–223. 34. Adams, S.; Rao, R. P. Transport pathways for mobile ions in disordered solids from the analysis of energy-scaled bond-valence mismatch landscapes. Phys. Chem. Chem. Phys. 2009, 11, 3010-3216. 35. Adams, S.; Rao, R. P. High power lithium ion battery materials by computational design. Phys. Status Solidi. 2011, 208, 1746–1753. 36. Van Der Ven, A.; Ceder, G. Lithium Diffusion in Layered LixCoO2. Electrochem. Solid-State Lett. 2000, 3, 301–304. 37. Van der Ven, A.; Ceder, G. Lithium diffusion mechanisms in layered intercalation compounds. J. Power Sources. 2001, 97–98, 529–531. 38. Seo, D.-H. Lee, J.; Urban, A.; Malik, R.; Kang, S.; Ceder, G. The structural and chemical origin of the oxygen redox activity in layered and cation-disordered Li-excess cathode materials. Nat. Chem. 2016, 8, 692– 697. 26 ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

39. Urban, A.; Lee, J.; Ceder, G. The configurational space of rocksalt-type oxides for high-capacity lithium battery electrodes. Adv. Energy Mater. 2014, 4, 1–9. 40. Isichenko, M. B. Percolation, statistical topography, and transport in random media. Rev. Mod. Phys. 1992, 64, 961–1043. 41. Black, R. Introduction to Percolation theory. Davenas, J. (eds). Percolation Theory. Martinus Nijhoff Publishers, Dordrecht. 1986. 42. Oikawa, K.; Kamiyama, T.; Izumi, F.; Nakazato, D.; Ikuta, H.; Wakihara, M. Neutron and X-Ray Powder Diffraction Studies of LiMn2−yCryO4. J. Solid State Chem. 1999, 146, 322–328. 43. Lu, Z.; Chen, Z.; Dahn, J. R. Lack of Cation Clustering in Li[Ni xLi1/3-2x/3Mn2/3-x/3O2 (0 < x ≤ 1/2) and Li[CrxLi(1-x)/3Mn(2-2x)/3]O2 (0 < x < 1). Chem. Mater. 2003, 15, 3214–3220. 44. Kondrashev, Y. D.; Zaslavskii, A. I. The structure of the modifications of manganese(IV) oxide. zv. AN SSSR, Ser. Fiz. 1951, 15, 179–186. 45. Tao, L.; Neilson, J. R.; Melot, B. C.; McQueen, T. M.; Masquelier, C.; Rousse, G. Magnetic structures of LiMBO3 (M = Mn, Fe, Co) lithiated transition metal borates. Inorg. Chem. 2013, 52, 11966–11974. 46. Akimoto, J.; Gotoh, Y.; Oosawa, Y. Synthesis and structure refinement of LiCoO2 single crystals. J. Solid State Chem. 1998, 141, 298–302. 47. Osorio-Guillén, J. M.; Holm, B.; Ahuja, R.; Johansson, B. A theoretical study of olivine LiMPO4 cathodes. Solid State Ionics. 2004, 167, 221–227. 48. Shklover, V.; Haibach, T. Crystal Structure of the Product of Mg2+ Insertion into V2O5 Single Crystals. J. Solid State Chem. 1996, 323, 317–323. 49. Arroyo-deDompablo, M. E.; Dominko, R.; Gallardo-Amores, J. M.; Dupont, L.; Mali, G.; Ehrenberg, H.; Jamnik, J.; Morán, E. On the Energetic Stability and Electrochemistry of Li2MnSiO4 Polymorphs. Chem. Mater. 2008, 20, 5574–5584. 50. Wadsley, A. D. Crystal chemistry of non-stoichiometric pentavalent vandadium oxides: crystal structure of Li1+xV3O8. Acta Crystallogr. 1957, 10, 261–267. 51. Ravel, B.; Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: Data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat. 2005, 12, 537–541. 52. Xu, L.; Liu, W.; Li, X.; Rashid, S.; Shen, C.; Wen, Y. Fabrication of MnOx heterogeneous catalysts from wet sludge for degradation of azo dyes by activated peroxymonosulfate. RSC Adv. 2015, 5, 12248–12256. 53. Bunau, O.; Joly, Y. Self-consistent aspects of x-ray absorption calculations. J. Phys. Condens. Matter. 2009, 21, 345501. 54. Pangaud, P.; Basolo, S.; Boudet, N.; Berar, J.-F. ; Benoit, C.; Clemens, J.-C.; Delpierre, P.; Medjoubi, K.; Hustache, S.; Menouni, M.; Morel, C. XPAD3-S : A fast hybrid pixel readout chip for X-ray synchrotron facilities. 2008, 591, 159–162. 55. Ashiotis, G.; Deschildre, A.; Nawaz, Z.; Wright, J. P.; Karkoulis, D.; Picca, F. E.; Kieffer, J. The fast azimuthal integration Python library: PyFAI. J. Appl. Crystallogr. 2015, 48, 510–519. 56. Davis, T., Farrow, C. L.; Billinge, S. J. L. PDFgetX3 : a rapid and highly automatable program for processing powder diffraction data into total scattering pair distribution functions. 2013, 46, 560–566. 27 ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

57. Bordet, P.; Colin, C.; Diaz-Lopez, M.; Fischer, H. E.; Freire, M.; Maignan, A.; Pralong, V. Local structure of the new high capacity cathode material Li4-xMn2O5 by neutron PDF analysis. Inst. Laue-Langevin (2017). doi:10.5291/ILL-DATA.6-06-473 58. Fischer, H. E.; Cuello, G. J.; Palleau, P.; Feltin, D.; Barnes, A. C.; Badyal, Y. S.; Simonson, J. M. D4c: A very high precision diffractometer for disordered materials. 162, 160–162 (2008). 59.

Howe, M. A.; Mcgreevy, R. L. CORRECT : A correction program for neutron diffraction data. 1998.

60. Tucker, M. G.; Keen, D. A.; Dove, M. T.; Goodwin, A. L.; Hui, Q. RMCProfile: reverse Monte Carlo for polycrystalline materials. J. Phys. Condens. Matter. 2007, 19, 335218. 61. Rodríguez-Carvajal, J. Recent advances in magnetic structure determination by neutron powder diffraction. Phys. B Condens. Matter. 1993, 192, 55–69. 62.

The MathWorks, I. Matlab and Statistics Toolbox release R2016b. (2016).

63. Momma, K.; Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 2011, 44, 1272–1276.

TOC graphic:

Supplementary information: Figure S1: Energy tresholds and voltage composition profile Figure S2: Quadrupole and dipole transitions Figure S3: Mn K-edge of standards Figure S4: Quantitative analysis of a single refinement Figure S5: Ellipsoidal fits of distorted octahedra Figure S6: XANES simulation of initial and refined models Figure S7: BVSE maps of several lithium transition metal oxides Figure S8: XANES simulation of Mn2O3 Figure S9: Data correction 28 ACS Paragon Plus Environment

Page 28 of 28