Full picture discovery for mixed-fluorine anion effects on high-voltage spinel lithium nickel manganese oxide cathodes

Small amounts of fluorine substituting for oxygen deficiencies could reduce Mn dissolution, enhancing the cyclability in spinel-type lithium nickel manganese oxides (LiNi0.5Mn1.5O4). Fluorine anion incorporation simultaneously enhances the C-rate capability and specific capacity fading. We used experimental and theoretical approaches to obtain a full picture of the mixed-anion effects for LiNi0.5Mn1.5O4−xFx cathode materials. The fluorine anion reduced the activation barrier for lithiumion hopping along the most energetically preferable 8a-16c-8a route, enhancing the C-rate capability. Simultaneously, the coordination bond of the linear F−–Mn3+–F− (Mn@2F diagonal) arrangement increased the oxidation potential to 5.1 V (vs Li+/Li). This hampered full extraction of Li+ from the spinel lattice, which was triggered by the oxidation of Mn3+ below the cutoff voltage (3.5–4.8 V (vs Li+/Li)), leading to a capacity loss.


Introduction
Lithium (Li)-ion secondary batteries have improved our lifestyles over the past two decades.1 However, the technical requirements for Li-ion secondary batteries are entering a phase of major change with respect to the energy density, which will require higher-energy-density cathode materials.2 The cathode's energy density depends on the specific capacity, tap density, loading amount and average operating voltage of the active materials. High-voltage materials could reduce the number of serially connected cells and enhance the energy density.
Spinel-type LiNi0.5Mn1.5O4 (LNMO) is a 5-V level high-voltage active material because of the high redox potential of Ni2+/Ni4+ (vs Li+/Li).3, 4, 5, 6 LNMO can be crystallographically classified into two spinel structures: ordered P4332 space group and disordered Fd-3m space group (cubic symmetry). Formation of an oxygen vacancy at the 24e site is critical for determining the symmetry and accompanying the reduction of Mn4+ to Mn3+ because of the charge balance. Recently, we performed density function theory (DFT) calculations to evaluate the formation energy using oxygen vacancy and metal-excess models for ordered and disordered LiNi0.5Mn1.5O4−δ (LNMO4−δ).7 Oxygen vacancy formation reactions were unlikely, although interstitial cation occupation occurred at the octahedral vacancies in both P4332 and Fd-3m LNMO4−δ spinel compounds.
The formation of Mn3+ in the LNMO lattice typically decreases the operating voltage, making the energy density lesser than that of stoichiometric LNMO.8,9,10,11,12 An excess of Mn3+ could produce structural distortions and disproportionations, degrading the cycle capability during highvoltage operations.13, 14, 15 Alternatively, the ordered P4332 phase has a limited power density. Low electron conductivity (two orders of magnitude lower than that of the Fd-3m LNMO) caused capacity fading in the high C-rate charge-discharge reaction. 12,16,17 We showed that tuning the Mn3+ concentration is important for suitable electrochemical characteristics. Substitution of the transition metal cation with Ni2+ and/or Mn4+ is a plausible way to tune the valence state of the manganese ion, which distributes defects in the spinel structure. Partial substitution of a cation (Cr, Fe, Mg, Zn or Ru) into the Ni and Mn sites provides a better performance than that of the pure structure. 18,19,20,21,22,23,24,25,26 Mixed-anion compounds (e.g., oxyfluoride) could enhance the electrochemical performance. A small amount of fluorine substitution for oxygen could reduce Mn dissolution (from hydrogen fluoride attack) and enhance the cyclability.27, 28 The incorporation of fluorine anions (F−) simultaneously enhanced the C-rate capability and capacity fading.29, 30 However, the origin of the enhanced electrochemical performance is not fully understood in relation to the electrical and structural properties,31, 32, 33, 34 which are affected by the electronic and ionic conductivities. Because of the higher electronegativity, the F− could interact with the Li+, causing an increased activation barrier for Li+ hopping along the 8a-16c-8a path in the spinel framework.
Here, we address this contradictory issue for LiNi0.5Mn1.5O4−xFx (LNMOF) cathode materials using experimental and theoretical approaches. Many fundamental studies have been carried out on cation-substitution effects, but there are only a few reports on mixed-anion effects. Therefore, determining the full picture for mixed-anion effects of high-voltage spinel systems provides new methods for material design that overcome the current performance limitations.

Materials and methods
Flux growth and characterization of LNMOF crystals LNMO4−δ crystals were grown by LiCl-KCl flux growth. 35 We performed a two-step flux growth for the LNMOF crystals with different amounts of fluorine substitution, including the flux growth of LNMO4−δ crystals from the LiCl-KCl flux and subsequent fluorine anion incorporation in a molten KCl flux. For the incorporation of the fluorine anion, the LNMO4−δ crystal was reacted with LiF in a KCl flux (800 °C, 20 h) at different LiF concentrations. The powders were washed with warm water to remove the remaining KCl flux and annealed under an O2 atmosphere (700 °C, 10 h). The phases and structures of the crystals were identified using X-ray diffraction (XRD) analysis with a Cu-Kα radiation source. The X-ray diffractometer (Miniflex II; Rigaku, Tokyo, Japan) was operated at 30 kV and 20 mA, with 2θ=10-80° (scan step=0.02). Field emission scanning electron microscopy (JSM-7600F at 15 kV; JEOL, Tokyo, Japan) observations were conducted to examine the microstructural characteristics in the crystals. The valence state of the Mn ion and atomic percent of incorporated fluorine atoms were studied by X-ray photoelectron spectroscopy (XPS) with a monochromic Al-Kα source (JPS-9010MX at 15 kV, 15 mA; JEOL). All binding energies in the spectra were referenced to the C1s hydrocarbon peak at 284.5 eV. The etching thickness was controlled at 20 nm per Ar sputtering, which was calculated in relation to the silicon wafer. The orders of the Ni and Mn configurations were characterized by Raman spectroscopy with 532 nm excitation (LabRAM, Horiba, Japan). The electrochemical characteristics were evaluated using a coin-type cell (R2032-type). The cathode was prepared by a conventional pasting process and contained LNMOF crystals, acetylene black and polyvinylidene fluoride (90:5:5 wt%). The mixture was diluted with N-methyl-2-pyrrolidone to give a viscosity of ~5.1 Pa·s. The prepared paste was homogeneously coated onto a 20-μm-thick Al foil using a slit coater. The electrode density was adjusted to 3.0 g cm−3 using a press machine. All electrodes were dried at 120 °C in a vacuum oven before cell assembly. Li metal foil and a porous polypropylene film (#2500; Celgard, Charlotte, NC, USA) were used as the counter electrode and separator, respectively. A solution of 1 M LiPF6 in a mixture of ethylene carbonate and dimethyl carbonate (1:1 vol.%) was used as the electrolyte. The coin-type cells were assembled in an Ar-filled glovebox (MDB-2BL; Miwa Mfg, Osaka, Japan) with a controlled atmosphere containing <1 p.p.m. H2O and O2. The galvanostatic charge-discharge tests and electrochemical impedance spectroscopy (EIS) were performed using a potentio/galvanostat (HJ1020Msd8 (Hokuto Denko, Tokyo, Japan); VSP-300 (Bio-Logic, Seyssinet-Pariset, France)). The cutoff voltage range for the battery tests was controlled between 3.5 and 4.8 V vs Li/Li+. The frequency range for the EIS measurements was 1 MHz-0.001 Hz. All EIS measurements were performed at 23 °C in a constant-temperature chamber (SU-221; Espec, Hudsonville, MI, USA).
X-ray absorption fine-structure spectroscopy Ex situ Mn K edge X-ray absorption fine-structure spectroscopy (XAFS) spectra were measured to evaluate the changes in the electrochemical states of the LNMOF crystal-based cathodes with different states of charge (SOCs) using the BL15A1 beamline of the Photon Factory (High Energy Accelerator Research Organization (KEK) in Tsukuba, Japan). The XAFS spectra were recorded from 6300 to 6900 eV in transmission mode. Each sample was measured for ~90 min. The spectral energy was calibrated using that of Mn2O3. LNMOF electrodes were removed from a disassembled coin cell after three cycles and used as samples for the XAFS measurements. The SOCs were controlled at 0 and 75%, respectively, before disassembly. The capacities of these SOCs were carefully estimated from their first-cycle discharge capacities. We disassembled the cells and washed them with dimethyl carbonate in an Ar-filled glove box. The samples were kept in a dried Ar atmosphere until the XAFS measurements.

DFT simulation
The DFT calculations were performed using the generalized gradient approximation (GGA-PBEsol)+U and projector-augmented wave methods as implemented in the Vienna ab initio simulation package.36, 37, 38, 39 The U-values for the d-orbitals of Ni and Mn were set to 6.0 and 3.9 eV, respectively.40, 41, 42 The magnetic ground state of LNMO exhibits ferrimagnetic ordering (Ni↓Mn↑), as reported in previous experimental and computational studies.31, 43, 44, 45 Thus, the ferrimagnetic spin configurations were used for all the calculations in this study. An energy cutoff of 500 eV and a 3 × 3 × 3 k-point mesh were used for the superstructure of 56 atoms in a cubic spinel lattice of 8(LiNi0.5Mn1.5O4) with P4332 symmetry as a starting structure. To realize the atomic arrangement of 8(LiNi0.5Mn1.5O4−xFx) (0.125≤x≤0.375), we added 1-3 fluorine anions into the possible 24e and 8c O oxygen sites in the unit cell of LiNi0.5Mn1.5O4−xFx (0.125≤x≤0.375). The nudged elastic band method was used to investigate the minimum energy pathways of Li hopping from one lattice position to adjacent sites.46, 47 Relaxation of the crystal structure was allowed for all calculations, and the final energies of the optimized geometries were recalculated to correct for changes of the plane-wave basis during relaxation.

Surface energy of the bulk model and surface model
We manually arranged the atoms on the LNMOF surface while retaining their crystallographic symmetry and chemical stoichiometry for calculating Eslab (the energy of a slab with a surface facet). The crystallographic symmetry of the top and bottom slab surfaces is essential for achieving a rational computational prediction. Therefore, the same analogy for calculating Eslab cannot be applied to the structures for x=0.125 and 0.375 because of their asymmetric O/F arrangement. In this work, we evaluated the surface energies for only LNMO3.75F0.25 (x=0.25).
The surface energies, γ, are defined as Direct observation of an original atomic surface structure requires clean and electron transparent thin scanning transmission electron microscopy (STEM) specimens. We simply put the LNMOF powders in ethanol and dispersed them onto a perforated amorphous carbon grid, rather than using the conventional Ar-ion thinning method to avoid significant surface damages. The atomic-resolution annular bright-field (ABF) and annular dark-field (ADF) STEM images were acquired using an aberration-corrected JEM ARM300CF instrument (JEOL) installed at the University of Tokyo, equipped with a cold-field emission gun and electron energy-loss spectroscopy (EELS, Quantum; Gatan, Pleasanton, CA, USA). To suppress the beam damage to Li-containing materials, we operated our microscope at 80 kV at a relatively low-dose condition (~13 pA), and no significant structure changes were observed during image acquisition. The used convergence semiangle was 27 mrad and the collection angles were 13-27 and 70-200 mrad for ABF-and ADF-STEM imaging, respectively.

Results and discussion
Figures 1a-c shows the most stable O/F arrangements for (a) x=0.125, (b) x=0.25 and (c) x=0.375 in LNMOF. We considered all possible O/F arrangements to determine the most stable arrangements; 2, 13 and 9 different O/F arrangements were calculated for x=0.125, 0.25 and 0.375, respectively. The energy differences between the most stable and second most stable arrangements for x=0.125, 0.25 and 0.375 were 0.48, 0.29 and 0.07 eV, respectively. F− preferentially occupied the 24e O sites instead of the 8c sites, independent of the amount of F− substitution. The same tendency was noted for oxygen vacancy formation in the non-stoichiometric LNMO4−δ system (e.g., the formation energy for 24e O site occupation in LiNi0.5Mn1.5O3.875F0.125 (LNMOF0.1) was ~0.48 eV lower than that for 8c O site occupation). F− coordinated with two Mn3+ ions and one Ni2+ ion. Furthermore, Mn4+ in the LNMO4−δ lattice was reduced to Mn3+ as the same number of F− ions was introduced into the structure. The incorporated F− reduced the neighboring Mn4+ to Mn3+. A linear F−-Mn3+-F− cluster formed in LiNi0.5Mn1.5O3.75F0.25 (LNMOF0.2) along the b-axis in octahedral MnO6 (24e sites occupied by F− (Mn@2F diagonal)). Further, incorporated F− occupied other nearest-neighboring 24e site. Two different configurations were formed in the LiNi0.5Mn1.5O3.625F0.375 (LNMOF0.3) lattice (linear and bent for F−-Mn3+-F− (Mn@2F neighbor) clusters in octahedral MnO6). The bond length of Mn3+-F− in the Mn@2F diagonal arrangement was longer than that of the Mn@2F neighboring arrangement, indicating that Jahn-Teller distortion possibly elongated the direction of the Mn@2F diagonal arrangement (parallel to the b-axis). Experimental and computational studies on the fluorine anion distribution inside the spinel frame work and the crystals of LiNi0.5Mn1.5O4−xFx (LNMOF). All possible O/F arrangements were considered to determine the most stable 2, 13 and 9 different O/F arrangements, which were calculated for x=0.125, 0.25 and 0.375, respectively: (a) x=0.125, (b) x=0.25 and (c) x=0.375 in LiNi0.5Mn1.5O4−xFx. (d) X-ray diffraction (XRD) patterns of LNMOF crystals with various amounts of fluorine anion incorporation. XRD patterns taken from all of the prepared samples are in good agreement with the reference data of ICDD PDF 70-8650 (LiNi0.5Mn1.5O4) with a small shift from fluorine addition. (e) The lattice parameters of LNMOF determined from the XRD patterns and density function theory (DFT) calculations. Series of F1s X-ray photoelectron core-level spectra: (f) LNMOF0.1, (g) LNMOF0.2 and (h) LNMOF0.3 with varying Ar spattering etching thicknesses (black: as prepared; red: 20 nm; blue: 40 nm; and pink: 60 nm).
The powder XRD patterns for all flux-growth samples were in good agreement with the reference data (ICDD PDF 70-8650; Figure 1d). All F−-substituted LNMOF crystals were of single phase. Chemical composition analysis based on inductively coupled plasma optical emission spectrometry revealed that all the LNMOF crystals had almost stoichiometric compositions (Supplementary Table  S1). The difference in the F contents for data obtained from the experiments and that used in the calculations are negligible from a crystallographic point of view. We used a superstructure  3)) for our computation. The chemical formula can be described as LiNi0.5Mn1.5O3.875F0.125, where one fluorine anion is substituted with O in the unit cell. We believe that such a negligible difference will not exert any significant effects on these findings. Figure 1e represents the lattice parameters determined from the XRD patterns and DFT calculations. There was little difference for the LNMO4−δ crystals with no F− substitution, which shows the theoretical validity of our model. The XRD experiments indicated lattice parameter changes that occurred with increasing F− substitution. The relatively small increase of the lattice parameter in highly fluorinated LNMOF may be from passivation of the oxygen deficiency by the F− at the defect site. Our calculations predicted a crystal structure transition from a cubic to orthorhombic system with F− substitution, accompanied by an anisotropic lattice dilatation along the b-axis (caused by Jahn-Teller distortion).
We examined the XPS F1s spectral depth profile to elucidate the incompatibility between the experimental and computation results for the lattice parameters (Figures 1f-h). The changes of the F1s core-level spectra as a function of Ar sputtering time represent the intensity of the XPS F1s corelevel spectrum, which became extremely weak after the 20 nm etching, and no XPS F1s signals were detected after 40 nm of etching. The remarkable changes in the depth profile indicated that F− ions were disproportionally incorporated at the crystal surface (which is strongly supported by the DFT calculation). Supplementary Figure S1 represents the energetically most stable atomic arrangements for the surface level at the (a) (1 0 0), (b) (1 1 0) and (c) (1 1 1) plane and for the core level at the (d) (1 0 0), (e) (1 1 0) and (f) (1 1 1) plane in LiNi0.5Mn1.5O4−xFx (x=0.27 for (1 0 0), x=0.25 for (1 1 0) and (1 1 1)). Supplementary Table S2 shows the calculated surface energies from these models; the surface energy of the bulk model was higher than that of the surface model, independent of the crystal faces. F− preferentially passivated the oxygen vacancy formed at the surface to stabilize the crystal surface. Moreover, the antisite energy for Ni/Mn ordering decreased with respect to F− substitution (Supplementary Table S3). The highest F− substitution (x=0.375) made the antisite energy negative, indicating that Ni/Mn disordering preferentially occurred and promoted the phase transition from an ordered to disordered structure.
Mixed F− influenced the morphology of the LNMOF crystals. The field emission scanning electron microscopy observation clearly visualized the formation of a step-terrace structure at the flat {1 1 1} plane and truncation at the vertices and edges of the parental octahedral structure (Supplementary Figure S1). The truncated lattice faces were assigned to the {1 0 0} faces, bearing a higher surface energy than the {1 1 1} face. Thus, the incorporated F− potentially reduced the surface energy of the {1 0 0} faces compared with the standard LNMO phase. This is most evident in the atomic-resolution ADF and ABF-STEM images observed in the [1 1 0] crystallographic orientation of LNMOF0.2 ( Figures  2a and b). The surface is not flat but has two components of the {1 0 0} and {1 1 1} surface steps, which is consistent with the field emission scanning electron microscopy observations in Supplementary Figure S2. A close inspection of the topmost surface indicates that a rock-salt like (similar to nickel oxide) atomic structure was formed at a thickness of a few unit cells, as evident from the cell-averaged inset images in Figures 2a and b: a new atomic site appears on the diamond center, as marked by arrowheads. However, in the {1 1 1} topmost surface, there appears antisite defect spinel structures that seem to form at a thickness of a few unit cells; transition metals might be occupied at Li atomic sites, as marked by arrowheads in Figures 2c and d. Based on the formation of the rock-salt and antisite spinel structure, we conclude that the amount of Li ions is lacking within a few unit cells from the surface because of their elution into water during flux removal.  The effect of F− on the ordering of the Ni/Mn arrangement in the LNMO4−δ lattice was studied using Raman spectroscopy. Intense Raman signals at 498 and 638 cm−1 (assigned to ordered Ni-O and Mn-O stretching bands, respectively) in the spinel lattice were diminished by the incorporation of F− independent of the substitution (Figure 3a). The mixed F− may have highly randomized the Ni/Mn arrangement. The proportion of Mn3+/Mn4+ examined by XPS (Figure 3b) also strongly supports our consideration. Two peaks attributed to Mn3+ and Mn4+ were observed in the Mn 2p core-level spectra. The Mn3+/Mn4+ ratio was 0.72 for the LNMO4−δ crystals and increased with F− substitution, which agrees with our computational predictions (1.15 for LNMOF0.1, 1.36 for LNMOF0.2 and 1.85 for LNMOF0.3). Our DFT calculation also predicted that Mn4+ was reduced to Mn3+ in the LNMO4−δ lattice when the same number of F− was introduced. Impact of fluorine anion incorporation on the Ni/Mn arrangements. (a) Raman spectra and (b) X-ray photoelectron spectroscopy (XPS) Mn 2p spectra of LNMOF crystals. The broad peak can be deconvoluted into two peaks, which were assigned to Mn3+ (blue) and Mn4+ (red) species. (c) Partial densities of state (PDOSs) for the Mn-3d band in LiNi0.5Mn1.5O4−xFx (0≤x≤0.375) (black: Mn@0F; blue: Mn@1F; red: Mn@2F diagonal; and dark green: Mn@2F neighbor).
Full figure and legend (166K) Figure 3c shows the calculated partial densities of state (PDOSs) for the manganese ions in LiNi0.5Mn1.5O4−xFx (0≤x≤0.375). At x=0.125, one Mn-3d band appeared near the Fermi level, indicating that the electron was doped in Mn-3d through the formation of a Mn@1F bond, giving Mn3+. Mn@2F and Mn@1F formed at x=0.25. The Mn-3d bands appeared near the Fermi level. The Mn-3d bands originating from Mn@2F shifted from −0.2 to −1.2 eV, compared with that at x=0.125. Thus, the Mn3+ ions were stabilized and became inactive redox couples against the applied voltage of 4.0 V (vs Li+/Li). This stabilization effect is likely associated with the enhanced Jahn-Teller distortion of Mn3+ by the two linearly arranged F− ions coordinated with Mn3+ along the direction parallel to the distortion. At x=0.375, additional Mn-3d bands for Mn@2F appeared at −0.6 eV. Moreover, the antisite energy for Ni/Mn ordering decreased with respect to F− substitution (Supplementary Table S3). Finally, the highest F− substitution (x=0.375) made the antisite energy negative, indicating that Ni/Mn disordering preferentially occurred and promoted the phase transition from an ordered to disordered structure.
Galvanostatic charge-discharge tests at various C-rates were systematically performed to examine the effect of F− substitution on the battery performance, including the rate capability. Figure 4a shows the third charge-discharge cycle profiles at 30 mA g−1, corresponding to a 0.2 C-rate showed that the F− substitution affected the voltage slope characteristics. Discontinuous sharp changes of the flat voltage were observed in the LNMOF crystals during the charge-discharge reactions from 4.6 to 4.9 V in the discharge reaction, indicating that two redox couples (Ni2+/Ni3+ and Ni3+/Ni4+) separately appeared. An inflexion point in the voltage slope at ~4.0 V was detected, which was assigned to the redox response from Mn3+/Mn4+. Supplementary Figure S3 shows the cyclic voltammetry of the LNMO4−δ and LNMOF crystals. From 4.6 to 4.9 V, two peaks were separated, corresponding to the Ni2+/Ni3+ and Ni3+/Ni4+ redox pairs, which were the origin of the voltage step at the lower Li-content region in the voltage profiles (Figure 4a).31, 48 The Mn3+/Mn4+ redox peak increased as the mixed F− increased (caused by charge valence). Meanwhile, the fluorine anion narrowed the potential gap compared with that of LNMO4−δ (shown in the redox potential gap for Ni2+/Ni3+). The reduction of the redox potential gap implies a faster lithiation-delithiation reaction in the disordered Ni/Mn configuration by incorporating F−.49 As shown in Figure 4a, the capacities faded with respect to F− substitution. The LNMOF0.3 cathode/Li cell had the lowest discharge capacity, which was three-quarters of the discharge capacity in the LNMO4−δ cathode/Li cell. The DFT calculations provided a full picture of this intriguing phenomenon. The changes in the most stable atomic structure as a function of the Li composition were computed using the formation energies of LNMOF with different atomic arrangements (Supplementary Figure  S4). The plausible reaction routes during the charge-discharge reaction and the voltage slopes were computationally predicted (Figure 4b). Our computational studies duplicated the voltage profile characteristics of the charging process from 3.5 to 5.0 V. The LNMO electrode showed a flat voltage as delithiation at ~4.7 V, assigned to oxidation of Ni2+ to Ni4+, whereas the LNMOF electrodes displayed two distinct plateaus at ~4.7 V and an additional small plateau at ~4.0 V attributed to the oxidation of Mn3+ to Mn4+. There was no intermediate phase once the oxidation state of all manganese ions reached the tetravalent state. Thus, the two-phase region coexisted during the oxidation of Ni2+ to Ni4+ in the charge process.31, 50, 51 The computational and experimental results indicated that the delithiation energy of LNMO was almost constant for 0≤Li≤1. In contrast, the formation of an oxygen deficiency in the LNMO lattice changed the manganese oxidation state from Mn4+ to Mn3+; a small portion of the manganese remained as Mn3+ in the LNMO/Li cell. In fact, crystal growth in a molten salt tends to occur with the lack of oxygen partial pressure, leading to an oxygen defect spinel compared with the solid-state reaction in an oxygen atmosphere.
Our calculations provide new information about the mixed-F− effects on the electrochemical reaction, including capacity fading. Several voltage steps were observed in both the higher and lower Li-content regions, including a phase transformation triggered by the Mn3+/Mn4+, Ni2+/Ni3+ and Ni3+/Ni4+ oxidation reactions. Higher delithiation energies than that of fully-lithiated LNMO (Li=1) are required in the lower Li-content region (Li≤0.125) because of the increased Li vacancy formation energy. The Li ions stored in the LNMOF lattice were not completely electrochemically extracted, causing capacity fade.
The computationally predicted PDOS of LNMOF with different Li compositions showed that the Mn3+ ions were oxidized preferentially to Mn4+ during the delithiation reaction. For instance, manganese ions show little changes in PDOS when the Li composition is reduced, and, hence, their contribution to the oxidation reaction during delithiation is very small. Figure 5 shows the change in the PDOS for the Mn-3d bands at different Li compositions. The manganese ions formed the Mn@2F arrangement (red line in Figures 5c and d) and did not change their original PDOS profile, even in the delithiated phase, which indicates that they remain in the Mn3+ oxidation state. Thus, Mn3+ coordinated with two fluorine anions in Mn@2F provided a small contribution to the redox reaction in Li0.125Ni0.5Mn1.5O4−xFx (0.25<x≤0.375) through an electrochemical delithiation reaction within the cutoff voltage range of 3.5-4.9 V.  Figure S5. A sufficient tendency to explain the effects of fluorine substitution on the capacity fade was observed even though the peak positions with different amounts of F− substitution and the peak shift at different SOCs were very small. The positions of the absorption peaks at SOC=0% shifted to the lower energy side with increasing F−. The formation of the Mn-F coordination bond reduced Mn4+ to Mn3+. Furthermore, all peaks shifted to a higher energy site after charging; thus, Mn3+ was oxidized to Mn4+. Increasing the F− substitution decreased the amount of shift. These features strongly support our computational predictions and considerations for the capacity fade of LNMOF.
The C-rate capability was performed to examine the mixed-fluorine-anion effect on the kinetics of Li ion transport. The LNMO4−δ cathode showed the highest discharge capacity at a low current density. However, it showed a greater capacity fade with an increasing current density than that for LNMOF ( Figure 6a). While the series of LNMOF cathodes showed a lower discharge capacity at a low C-rate, their capacity dropped because the increased internal resistance with increasing current density markedly improved. The incorporation of F− affects the electrochemical Li intercalation kinetics, which are strongly associated with the bulk properties of the electrode materials (e.g., electron conductivity and Li ion conductivity), and the kinetic parameters at the electrode-electrolyte interface. Lithium ion transportational characteristics in the LiNi0.5Mn1.5O4−xFx (LNMOF) electrodes. (a) Crate capabilities and (b) Nyquist plot. The kinetic parameters were determined by electrochemical impedance spectroscopy (EIS) approaches. The Nyquist plots were fitted using the equivalent circuit model. The high-and low-frequency semicircles are attributed to surface film resistance (Rsf) and charge transfer resistance (Rct) at the electrode-electrolyte interface, respectively. The Li-ion diffusion coefficient from the relation between the real impedance and the angular frequency in the low-frequency region were further estimated to evaluate quantitatively the F−-substitution effects. (c) Lithium ion hoping energy as a function of F content. The calculated energy profiles for lithium-ion hopping along the 8c-4a/12d-8c route with the smallest activation barrier for ion hopping in LNMO and series of LNMOF crystals (12d-3 model in x=0.125, 4a model in x=0.25 and 4a model in x=0.375 in Supplementary Figure 6).

Full figure and legend (205K)
The kinetic parameters were determined by EIS approaches. Figure 6b shows the Nyquist plots of the LNMO4−δ and LNMOF electrodes. The impedance spectra were fitted using the equivalent circuit model. The high-and low-frequency semicircles are attributed to surface film resistance (Rsf) and charge transfer resistance (Rct), respectively, at the electrode-electrolyte interface. We further calculated the Li-ion diffusion coefficient from the relation between the real impedance and the angular frequency in the low-frequency region to evaluate quantitatively the F−-substitution effects. The kinetic parameters are summarized in Supplementary Table S4. Increasing the mixed-F− content effectively reduced Rsf and Rct, and they reached minimum values for LNMOF0.3. Further F− substitution increased the resistance, which might correspond to disturbed Li ion migration by excess amounts of F−.
The relationship between Li ion migration and fluorine anion substitution was considered using DFT calculations to predict Li-ion migration paths (Figure 6c). The Li ion in the spinel framework preferentially migrated along the 8c-4a/12d-8c route, independent of F− substitution. In the considerable migration models shown in Supplementary Figure S6, the activation barrier for ion hopping in the LNMO and series of LNMOF crystals was the smallest in the 12d-3 model at x=0.125 (0.37 eV), 4a model at x=0.25 (0.33 eV) and 4a model at x=0.375 (0.44 eV). In particular, the LNMOF0.2 cathode showed minimum activation energies and the highest Li ion conductivity, which agree with the C-rate capability test.
The migration paths and their corresponding activation energy profiles for Li+ hopping along 8c-4a/12d-8c in stoichiometric LNMO, non-stoichiometric LNMO4−δ and LNMOF0.1 are summarized in Supplementary Figure S7. These results indicate that the migration path and the activation energy for Li+ hopping strongly depend on the interaction between at the 8c-4a/12d-8c site and or at the 24e site. The Li+ ion in stoichiometric LNMO moves in a straight line along the migration path, as shown for spinel Li4Ti5O12.46 In contrast, Li+ in LNMO4−δ and LNMOF0.1 migrates tortuously at and , respectively. Replacing with mitigated the tortuous Li+ migration path and lowered the activation energy from 1 eV to ~0.5 eV, as required for Li+ hopping because the coordination number of the Li+ ion at 8c-4a/12d-8c in LNMOF0.1 remained three (two oxygen and one fluorine), even though that of LNMO was two (two oxygen).
Li-ion hopping can be interpreted as Li ion vacancy ( ) hopping. Generally, the valence state for is negatively charged. Therefore, slightly changed the charged state of the 4a/12d site ( ) and became more positive because of the different valence states. The activation barrier for Li-ion hopping along the 8c-4a/12d-8c path became larger because the Li ion path was more positively charged. In other words, is strongly trapped with because of defect aggregation. In contrast, the DFT calculations revealed that incorporating F− strongly affected Li ion conductivity in nonstoichiometric LNMO4−δ with . The charged state of the 4a/12d site in non-stoichiometric LNMO became much more positive after the formation of a at the 24e site compared with that of stoichiometric LNMO. Replacing at the 24e site with reduced the activation energy for Li-ion hopping because of the relaxed electrostatic repulsion at the 4a/12d site of LNMO4−δ, which was positively charged.
Fluorine anion substitution-dependent changes to the Li-ion hopping route were not seen; however, their activation barrier changed significantly with respect to the symmetry in the site potential of the two tetrahedral 8c sites adjacent to the octahedral 4a/12d sites in the Li-ion diffusion path. A symmetric site potential across 8c-4a-8c was formed in LNMOF0.2 by the substitution of a single F− in both tetrahedral 8c sites, leading to a site potential minimum (which is in contrast to those of asymmetric LNMOF0.1 and LNMOF0.3, whose activation barriers were higher than that of LNMOF0.2). Oxygen vacancies in the two tetrahedral 8c sites adjusted to the octahedral 4a/12d site and were inhomogeneously substituted with F− to form different Li-F coordination numbers in LNMOF0.1 (Li8cx@1F24e-Li8cx@O24ex) and LNMOF0.3 (Li8cx@1F24e-Li8cx@2F24e), which imbalanced the site potential and increased the activation barrier for ion hopping. These calculations agreed with the experimentally demonstrated C-rate capability and EIS results. Figure 7a shows the discharge capacity vs cycle number for LNMO4−δ and a series of LNMOF cathodes/Li cells cycled at 23 °C with a current corresponding to 0.2 °C. A significant capacity loss with cycle number was observed to 80% of the initial capacity for the LNMO4−δ/Li cells after 100 cycles; all of the fluorine-substituted LNMOF/Li cells showed an obvious improvement in the cyclability with capacity retentions >90%. The LNMOF0.2/Li cell showed almost no capacity loss during cycling. The coulomb efficiency of the cell was almost constant at ~0.990 during the cycles. Furthermore, EIS experiments for the cells after 100 cycles (Figure 7b) showed that the all-kinetic parameters (Rsf, Rct and DLi) degraded with cycling; however, the deterioration of the LNMOF0.2/Li cell was significantly smaller (Supplementary Table S5). Fluorine anion substitution likely suppressed the dissolution of the manganese ion into the electrolyte and metal-fluoride formation as a solid electrolyte interface by oxidative decomposition of the liquid electrolyte at high voltages. The computational studies implied that the stabilization of Mn3+ by neighboring F− (formation of coordination bonds and enhancement of Jahn-Teller distortion) primarily contributed to the significant cyclability improvement. Furthermore, we reported that surface modifications of the LNMO crystals with <2-nm-thick fluoroalkylsilane monolayers improved the cyclability by enhancing the lyophobic characteristics against liquid electrolytes.35 The mixed F− lowered the surface energy by inducing segregation near the surface. A similar surface modification was imparted by the fluoroalkylsilane monolayer and will contribute in improving the cycle characteristics. Based on the changes in the Mn3+/Mn4+ ratio of the XPS-Mn core-level spectra taken from the electrode surface after cycling, we showed that Mn3+ does not go through the notorious disproportionation reaction. The absence ratio of Mn3+ was sufficiently retained in the LNMO stabilized with the fluoroalkylsilane monolayer compared with that of bare LNMO. Furthermore, retention of dark pigments was observed on the polypropylene separator after 100 cycles in the LNMO/Li cell, unlike that in the LNMOF0.2/Li cells. The batteries were reusable after replacing an old separator with a new one. These results strongly suggest that the incorporation of F− into the LNMO lattice prevents the risk for higher capacity fading under high-voltage operation because of the occurrence of a short circuit caused by Mn3+ dissolution during cycling.

Figure 7.
Cycle performances of the LiNi0.5Mn1.5O4−xFx (LNMOF) electrodes. (a) Cyclabilities at a 0.2 C-rate, and the discharge capacity vs cycle number for LNMO4−δ and a series of LNMOF cathodes/Li cells cycled at 23 °C with a current corresponding to 0.2 °C. (b) Impedance spectra taken after the 100th cycling test at a 0.2 C-rate and 23 °C.

Full figure and legend (87K)
In summary, we determined the effect of mixed-fluorine anions in oxygen-deficient, spinel-type LNMO4−δ crystals on the structural, electrical and electrochemical characteristics by experimental and computational approaches. This is the first report that provides the full picture for mixedfluorine-anion effects on high-voltage spinel manganese cathodes. The anion-mixed effects are an exciting new direction for the development of high-voltage cathode materials.