Article Cite This: Chem. Mater. 2018, 30, 3060−3070
pubs.acs.org/cm
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† †
Université Grenoble Alpes, Institut Néel, F-38000 Grenoble, France CNRS, Institut Néel, F-38000 Grenoble, France 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 § Institut Laue-Langevin, 71 avenue des Martyrs, CS 20156, 38042 Grenoble cedex 9, France Downloaded via KAOHSIUNG MEDICAL UNIV on June 30, 2018 at 19:01:38 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.
‡
S Supporting Information *
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, that is, when Li+ ions are inserted and extracted from the cathode. Moreover, new cathode designs involve more and more disordered/nanosized 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 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.
T
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 model.17 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
he 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 performance.2 Among 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 Mn,14 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 nonstoichiometric Li4Mn2O5 (Li4) rock-salt structure is not yet understood, partly due to the lack of a detailed knowledge of © 2018 American Chemical Society
Received: February 26, 2018 Revised: April 19, 2018 Published: April 19, 2018 3060
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials (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, the PDF analysis was carried out using the RMC modeling method, capable of generating a high degree of local disorder from highly ordered starting configurations.18−21 Several refinements were performed on a “large-box” with dimensions similar to the structural coherence length of Li4,17 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 Mnframework 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 units.25 The BVSE map of Li revealed a plausible 3D Li diffusion pathway connecting large volumes of 5-coordinated lithium clustered around oxygen vacancies.
■
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.
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 analyzed. 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 Li0.17 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 Supporting Information (SI) 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 rule,26 which is in good agreement with the smaller cell parameters of Li0 determined by real and reciprocal space Rietveld refinements.17 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 1b). Although small in amplitude, the pre-edge signal carries valuable information regarding the coordination environment of the absorber atom, that is, 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 SI 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 noncentrosymmetric 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 SI 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, 3061
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
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 pre-edge peaks was attributed to the distribution in energy of the t2g and eg states given by the eight nonequivalent Mn sites in the initial model. 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 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 modeling 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 cutoff 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-
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
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 eight 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.
Figure 3. RMC fits of Li4 neutron (top) and X-ray (bottom) G(r) (left) and S(Q) (right) at 297 K. 3062
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
evaluation of only 6-coordinated manganese and 6 and 5coordinated lithium was carried out. For the interpretation of the RMC outputted model the partial atomic pair distribution function gij(r)s is given in Figure 5a, where we observe the contribution of the different pairs of
range structures to be simultaneously modeled. 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). 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 10 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 SI Figure S4). We first evaluate the reliability of the refined models by computing the BVS of each atom (see Figure 4). A narrow
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.
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 Mnframework, corroborating our expectations.17 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 6b). The first two peaks in Figure 5a centered at ∼ 2 Å account to the nearest Li-O and Mn-O neighbors’ distances (see also Figure 6a 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 centered 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 6b). 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/R1,30 was used to address the shape of the octahedra in the refined model as: axially compressed or prolate (S < 0), axially stretched or oblate (S > 0) or regular (S = 0). For a representation of the ellipsoidal fits of distorted octahedral the interested reader is referred to SI Figure S5. In Figure 7, the plot of S for Mn shows a narrow and symmetric distribution
Figure 4. BVS plots for (a) oxygen in gray and for (b) manganese (c) lithium in gray, orange and blue for 6-, 5-, and 4-coordinated sites respectively. The vertical red lines indicate the initial BVS values in the starting model.
distribution (∼±0.5 v.u.) around the expected BVS values was observed for oxygen (1.9 v.u.), 6-coordinated manganese (2.8 v.u.) and 5-coordinated lithium (1.0 v.u.). The BVS of the minority 6-coordinated 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 3063
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
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 (σM‑O) given as a function of the coordination number.
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
Figure 8. Histograms of □-Li (a) and □-O (b) distances. The red vertical lines indicate the initial values on the starting model.
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 6a, and BVS in Figure 4b. For Li, as well as for Mn, a narrow and symmetric distribution of the S value around zero is observed, which suggest that Oh-coordinated Li sites also tend to have spherical shapes. 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 8, 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 8, we can conclude that oxygen within the pyramids base are drawn in toward 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 Å. The small deviations of the □-Li distances from their initial values do 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 centered at 1.87 and 2.15 Å in Figures 5a and 6b, we computed the projection of the Li shift along and perpendicular to the □-Li direction (see Figure 9a and b). The projection of the refined lithium position along the □-Li direction shows a symmetrical distribution centered 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 toward 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 3064
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
Figure 9. Displacement of lithium was calculated with respect to its initial configuration at the center 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-Liinitial direction. (a) Projection of the refined lithium position along the □-Li direction, where positive values denote the displacement of lithium toward 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.
of the shift in Figure 9c, 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. 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 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 SI Figure S6) are a robust proof of the validity of the RMC refined model. The creation of a sensible starting configuration for RMC modeling 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 configuration;19 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 coordination environment of manganese. We introduce here the two refinements of the initial models in Figure 2a 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
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.
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 modeling approach can successfully 3065
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials Table 1. Goodness of Fit and Calculated BVS for the RMC Refinement Using Different Initial Configurations initial/refined BVS refinement ID
initial % of 5-coordinated Mn
cations swapping (Y/N)
χ2
Mn
Li
O
A B C D
0 50 100 0
N N N Y
8.5(1) 8.6(1) 8.5(2) 8.1(7)
2.5062/2.719(3) 2.2813/2.609(3) 2.1207/2.529(4) 2.5062/2.644(9)
0.9353/1.041(3) 1.0230/1.078(2) 1.0726/1.110(2) 0.9353/1.062(4)
1.7482/1.901(1) 1.6929/1.875(1) 1.6534/1.861(1) 1.7482/1.881(3)
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 among 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). 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 center 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 9c. 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 potential,31 thereby transforming valence units into energy units. As for BV mismatch maps,32 the graphical representation of isosurfaces with constant BVSE can provide a clear image of the lithium diffusion pathways in the material studied.33−35 The BVSE map in Figure 11 was performed in a 4ax4ax4a (where a is the lattice parameter of the “MnO” cell) section of the RMC refined model 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 12b 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 12d 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 neighboring 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 state,38 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 representative state-of-the-art cathode materials introduced in Table 1, point out the high potential of Li4 warranting further investigations into this system. Table 2
■
CONCLUSION The complete understanding of the local disorder in nanostructured Li4 was reached via the RMC modeling 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 5-coordinated Li is displaced away from the center 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. Five g of nanostructured Li4 were synthesized by a twostep route, the first of which involved the synthesis of HT-LiMnO2 by the solid-state reaction method. Stoichiometric amounts of MnO, 3066
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
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 centered at the neighboring O. 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 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 ICP-OES. Near-Edge X-ray Absorption Spectroscopy and FDMNES Modeling. 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 software.51 Our data processing protocol is supported by the good agreement between our spectra of the standards (see SI Figure S3) and others found in the literature.52 Experimental XANES spectra were simulated with the FDMNES software53 using the multiple-scattering theory under the muffin-tin approximation on the potential shape, that is, the potential is spherically averaged in the atoms and constant between them. The
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. 3067
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
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.)a percolation E (eV) along
a
material
F (%)
minimum E (eV)
a
b
c
ref
Li4 LiMn2O4 LiNi0.5Mn0.5O2 α-MnO2 LiMnBO3 LiCoO2 LiFePO4 olivine V2O5 Li2FeSiO4 LiV3O8
14(3) 13.28 17.13 19.44 4.5 11.17 9.86 10.48 19.44 11.67
−5.3(4) −4.4175 −3.1193 −3.0518 −8.9629 −5.1456 −4.586 −4.349 −3.8476 −5.0514
2.0(6) 0.54 0.30 No No 1.96 No No 0.96 No
2.0(6) 0.54 0.30 No No 1.96 1.04 1.38 No 1.34
2.0(6) 0.54 0.30 0.18 1.81 No No 2.26 0.96 2.46
this work 42 43 44 45 46 47 48 49 50
The calculated BVSE maps are given in SI Figure S7.
code runs in direct space and the calculations were limited to a cluster embedded in a sphere centered on the absorbing atom. When several nonequivalent 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 a cluster size of 6.1 Å around the absorber atom corresponding to ∼ 90 atoms, the spectra of the Mn3+2O3 standard was successfully reproduced (see SI Figure S8). Simultaneous RMC Modeling of n-PDF and x-PDF. Highenergy synchrotron X-ray powder diffraction (HE-SXRPD) 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 collection time using a 2D pixel detector XPAD3.54 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 diffractometer 58 at the Institut Laue-Langevin (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 software.59 The corrected data was Fourier transformed using the STOG software embedded in the RMCProfile suite60 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 software.60 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 SI 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 MatLab.62 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 cutoff effects at the boundaries of coordination shells), which could lead to the formation of artifacts in the BVSE maps. 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-Suite.61 The figures illustrating the local environments and diffusion pathways in this publication were produced with the VESTA software.63
■
ASSOCIATED CONTENT
S Supporting Information *
The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.chemmater.8b00827. 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 (PDF)
■
AUTHOR INFORMATION
Corresponding Author
*E-mail:
[email protected]. ORCID
Maria Diaz-Lopez: 0000-0001-5670-1859 Yves Joly: 0000-0002-1872-6112 Claire V. Colin: 0000-0003-1332-7929 Valerie Pralong: 0000-0003-4644-8006 Pierre Bordet: 0000-0002-1488-2257 3068
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials Author Contributions
(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 Li-excess 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, 11, 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. 2016, 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. Mol. Simul. 1988, 1, 359−367. (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. Curr. Opin. 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 Crystallogr., Sect. B: Struct. Sci. 2001, 57, 278−287. (24) Adams, S. softBVS web page, 2003; 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 Ionics 2006, 177, 1625−1630. (26) Bianconi, A.; Dell’Ariccia, M.; Gargano, A.; Natoli, C. R. Bond length determination using XANES. In EXAFS and Near Edge Structure; Springer Series in Chemical Physics: Berlin, 1983; Vol. 27, pp 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, 155, 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. In Bond Valences. Structure and Bonding; Brown, I.; Poeppelmeier, K., Eds.; Springer: Berlin, Heidelberg, 2013; Vol 158. (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. Bonding (Berlin, Ger.) 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 A 2011, 208, 1746−1753. (36) Van Der Ven, A.; Ceder, G. Lithium Diffusion in Layered LixCoO2. Electrochem. Solid-State Lett. 2000, 3, 301−304.
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 setup at D2AM, M.D., and N.B. 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. Notes
The authors declare no competing financial interest.
■
ACKNOWLEDGMENTS
■
REFERENCES
We 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.
(1) Larcher, D.; Tarascon, J.-M. Towards greener and more sustainable batteries for electrical energy storage. Nat. Chem. 2015, 7, 19−29. (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, 30, 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. Li2MnO3-stabilized 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 V 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. 3069
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070
Article
Chemistry of Materials
(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 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.
(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. (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. In Percolation Theory; Davenas, J.; Ed.; 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. Nucl. Instrum. Methods Phys. Res., Sect. A 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. J. Appl. Crystallogr. 2013, 46, 560−566. (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. LaueLangevin 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. Appl. Phys. A: Mater. Sci. Process. 2008, 162, 160−162. 3070
DOI: 10.1021/acs.chemmater.8b00827 Chem. Mater. 2018, 30, 3060−3070