Subsurface Distribution of Antisite Defects in LiMnPO4: Direct

Oct 26, 2016 - We demonstrate a thermodynamically stable distribution of antisite exchange defects within a single LiMnPO4 crystal through a combinati...
0 downloads 4 Views 2MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Subsurface Distribution of Antisite Defects in LiMnPO: Direct Comparison with LiFePO 4

4

Sung-Yoon Chung, Si-Young Choi, Tae-Hwan Kim, and Seongsu Lee J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b09730 • Publication Date (Web): 26 Oct 2016 Downloaded from http://pubs.acs.org on October 31, 2016

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry C 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 17

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

The Journal of Physical Chemistry 1

Subsurface Distribution of Antisite Defects in LiMnPO4: Direct Comparison with LiFePO4 Sung-Yoon Chung,*,†, 4 Si-Young Choi,‡, 4 Tae-Hwan Kim,§ and Seongsu Lee§



Graduate School of EEWS, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea ‡ Korea Institute of Materials Science, Changwon 51508, Korea § Korea Atomic Energy Research Institute, Daejeon 34057, Korea

*Corresponding Author: [email protected] and [email protected] Theses authors contributed equally to this work.

4

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 2 of 17 2

Abstract We demonstrate a thermodynamically stable distribution of antisite exchange defects within a single LiMnPO4 crystal through a combination of atomic-resolution scanning transmission electron microcopy and ab initio calculations. In contrast to the strong segregation of exchange defects in subsurface layers with selective orientations in LiFePO4, the highly ordered arrangement between Li and Mn in the bulk is preserved near the surface of LiMnPO4 crystals, showing no preferential aggregation of defects irrespective of the location in a particle. Based on the identical crystal structure and the analogous compositional basis of the two olivine phosphates, the distinct distribution characteristics of the antisite exchange defects are notably unexpected findings, highlighting the value of direct visualization at an atomic scale.

ACS Paragon Plus Environment

Page 3 of 17

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

The Journal of Physical Chemistry 3

Introduction An ordered configuration of Li and transition metals (TM) among the oxygen interstitials is a significant structural feature in many Li-intercalation oxides, and results in rapid Li-ion mobility within a close-packed oxygen framework during the electrochemical Li (de)insertion reactions. In addition to the macroscopic information provided by the conventional powder diffraction and refinement analyses,1,2 direct identification of local variation in the cation ordering and segregation has been a central issue in structural chemistry dealing with Li-intercalation oxides over the past decade.3−9 In particular, unequivocal experimental evidence based on atomic columnby-column imaging in aberration-corrected scanning transmission electron microscopy (STEM) has provided an important contribution to a better understanding of the correlation between cation disordering and resulting electrochemical performance. Many notable findings observed via STEM at an atomic scale have been reported since 2008, visualizing the Li−TM antisite exchange defects,10−16 the local compositional segregation,8,17,18 and even structural transitions19−25 in various layered rocksalt-, olivine-, and spinel-type Li-intercalation oxides. The Li (de)intercalation reaction is initiated at the surface of particles during electrochemical cycling. As a result, the apparent Li diffusion and subsequent reaction kinetics can be greatly affected by the local Li−TM configuration perturbed from the complete ordering in the surface regions. A substantial number of studies have thus recently focused on elucidating the atomic configurational variations near the particle surface.5,8,9,15,16,21,22,25 Among many intriguing aspects, a strong orientation-dependent distribution of the Li−TM intermixing at the surface is one of the most noteworthy features.15,16,25 Furthermore, such anisotropy is not confined to a specific oxide system but rather commonly occurs in many layered oxides encompassing Li1+x(Mn,Ni)O2 and Li(Mn,Ni,Co)O2. Both the composition and the atomic structure at each surface in

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 4 of 17 4

ordered Li-intercalation oxides vary remarkably with the crystallographic orientation of the surface. Consequently, the facet-dependent anisotropic formation of Li−TM intermixing near the surface is understood in terms of an energetically favorable configuration for the thermodynamic stability of entire particles, as exemplified in Li1+x(Mn,Ni,Co)O2.16 Very recently, analogous behavior in the subsurface distribution of antisite exchange defects was demonstrated in olivine-type LiFePO4.25 For example, while LiFe−FeLi exchanges were hardly observed in the (100) surface region during the STEM analysis, a substantial amount of exchange defects were clearly imaged in the (001) surface regions, showing the highly anisotropic formation of defects. In this article, we report that antisite defects are not detected in significant amounts in the surface regions of LiMnPO4 particles, in stark contrast to LiFePO4. More importantly, such LiMn−MnLi exchanges are remarkably insensitive to the surface facet orientation. Recalling the compositional analogy and identical crystal structure of LiMnPO4 and LiFePO4, these two distinct subsurface defect formation characteristics appear to be unusual.

Experimental Methods Particle synthesis. LiMnPO4 and LiFePO4 crystalline particles were prepared via a solid-state reaction by using high-purity lithium carbonate (Li2CO3, Aldrich), manganese oxalate dihydrate (Mn(II)C2O4⋅2H2O, Alfa Aesar), iron oxalate dihydrate (Fe(II)C2O4⋅2H2O, Aldrich), and ammonium dihydrogenphosphate (NH4H2PO4, Aldrich). Stoichiometric powder mixtures of the starting materials were ball-milled in acetone for 24 h with zirconia milling media. Dried slurries were calcined at 350°C for 5 h under a flow of high-purity Ar (99.999%, 400 sccm). The calcined amorphous powder samples were further annealed at 500°C for LiMnPO4 and at 600C for LiFePO4 in the same Ar atmosphere for a sufficiently long time, 72 h, to exclude kinetic formation of exchange defects during synthesis26 and thus obtain crystalline particles

ACS Paragon Plus Environment

Page 5 of 17

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

The Journal of Physical Chemistry 5

having a thermodynamic equilibrium state of defect configurations. HAADF-STEM. Atomic-column-resolved HAADF-STEM images were acquired using a transmission electron microscope with a spherical-aberration corrector (CEOS GmbH) for an electron probe. The size of the probe in STEM mode was 0.96 Å. The collection semiangles of the HAADF detector were adjusted from 71 to 190 mrad in order to utilize incoherently scattered electrons at large angles. The obtained raw images were band-pass filtered to reduce background noise. DFT calculations. Ab initio calculations were carried out within the spin-polarized generalized-gradient approximation (GGA) along with the Perdew-Burke-Ernzerhof (PBE) functional for exchange correlation and the ultrasoft pseudopotentials for ionic cores, as implemented in the CASTEP code (Biovia). A sufficiently long slab along with a 10-Å vacuum layer was constructed as an optimum supercell for each calculation to make the relaxation layer of each slab more than 10 Å in thickness. As previously suggested, stoichiometric [LiMnPO4]x-type slabs were adopted as proper supercells for reasonable determination of the surface termination.25 To exclude any energy variation induced by change of the surface termination, a pair of LiMn−MnLi exchange defects was introduced into the second Li and Mn rows beneath the top-most surface plane, preserving the initial surface termination. The formation energy of a defect pair in the subsurfaece layer (ES) was obtained by ES = E(defects) − E(w/o defects) where E(defects) and E(w/o defects), respectively, are the total energies of a supercell with and without a defect pair after geometry optimization. The plane-wave basis set for the kinetic energy cutoff was 450 eV. Optimization of the internal coordinates for each case was performed using the BFGS algorithm with convergence tolerances of 0.1 eV/Å for the maximum ionic force, 5 × 10−5 eV/atom for the total energy, and 0.005 Å for the maximum ionic displacement.

Results and Discussion Figure 1 shows the atomic-scale HAADF-STEM images obtained near the (001) vicinal

surface

planes

of

LiFePO4

and

LiMnPO4

ACS Paragon Plus Environment

particles,

respectively.

The Journal of Physical Chemistry

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 6 of 17 6

Comprehensively different formation behavior of subsurface exchange defects between the two phosphates is readily recognized in this set of images. As consistently shown in the previous study,25 substantial column intensity in most of the Li sites (green arrows) is detected in the surface region of a LiFePO4 particle (Figure 1a), indicating the presence of a significant amount of antisite exchange defects. In contrast, Figure 1b reveals that the ordered image feature with invisible Li columns and bright Mn columns is maintained over the surface region of a LiMnPO4 particle. Very few Li columns showing sufficiently detectable intensity near this surface were found during the STEM observation. As demonstrated in Figure 2, as an example, merely several Li columns out of more than a hundred columns in a magnified image exhibited intensity that exceeds the background noise, verifying the distinct distribution of antisite defects in the subsurface region between LiMnPO4 and LiFePO4. For further comparison, intensity profiles of Li columns in LiMnPO4 and LiFePO4 are demonstrated in Figure S1 in Supporting Information. To examine whether any surface-orientation-dependent characteristics of the defect distribution are present in LiMnPO4, atomic-resolution HAADF-STEM images were collected near various surface facets of a single particle. The low-magnification STEM image in Figure 3 shows a LiMnPO4 particle in the [010] projection. A series of HAADF-STEM images were subsequently obtained from the three different surface regions denoted by orange rectangles with the numeric labels of 1, 2, and 3. The most noticeable aspect in this set of images is that the highly ordered configuration between Li and Mn observed in the bulk is preserved well up to the top-most surfaces, irrespective of the surface-facet orientation. Few Li columns with detectable intensity stemming from Li−Mn exchanges are identified in all three regions. The maintenance of a high degree of Li−Mn ordering demonstrated in Figure 3 implies

that

the

formation

of

exchange

defects

ACS Paragon Plus Environment

has

similar

likelihood

Page 7 of 17

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

The Journal of Physical Chemistry 7

thermodynamically regardless of the location in a particle. We thus carried out ab initio DFT calculations to theoretically compare the defect formation energies near surfaces of a different crystallographic orientation. As adopted in previous studies of olivine phosphates,25,27 (LiMnPO4)x-type stoichiometric slabs with a vacuum layer were utilized as proper supercells. In addition, an exchange-defect pair, LiMn−MnLi, was introduced into the subsurface layers, not at the top-most surface layers, in order to maintain the initial surface termination of each slab. Figure 4 illustrates the supercell models with seven different surface-facet orientations for DFT calculations. As denoted by a red ellipse in each of the cells, a pair of LiMn−MnLi was introduced beneath the top surface rather than the very first Li and Mn rows. Based on these two approaches for supercell construction, we can avoid unnecessary concern regarding the configurational change of surface termination. We can therefore exclude possible unknown influence induced by the termination change on the total energy variation, thereby allowing a reasonable comparison of the formation energies between subsurface regions. Three crystallographically distinct neighboring Mn atoms are placed around each Li in this ordered olivine-type lattice (see Figure S2 in Supporting Information). We thus calculated the LiMn−MnLi formation energies for 21 combinations in total, considering the three possible exchange geometries for each of the seven supercells. The lowest formation energy (Es) among the three values is listed in Table 1 in each case (see Table S1 in Supporting Information for the entire set of formation energies). For a direct comparison, the Es values of LiFePO4 obtained in our recent study in the same manner are also included in the right column of the table. One of the notable findings from the DFT calculations for LiMnPO4 is that the defect formation energy does not considerably vary with the surface orientation. While highly negative values of Es are obtained for the (001) and (111) surface planes of LiFePO4, the calculation results of LiMnPO4 in Table 1 show no negative formation

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 17 8

energies in any surface plane. Although the defect formation energies for the (010) and (001) surfaces in LiMnPO4 are comparatively lower than that in the bulk (0.75 eV), the overall Es does not deviate from a range of 0.39−0.76 eV, in striking contrast to the strong surface-orientation-dependent defect formation in LiFePO4. This relatively small variation in the defect formation energy in LiMnPO4 also agrees well with the HAADFSTEM observations shown in Figures 1 and 2, consistently supporting a distribution of antisite exchange defects insensitive to the location within a particle. Precise and direct probing of the initial atomic arrangements of Li and TM in pristine Li-intercalation oxides prior to any electrochemical cycling is of great significance as it is the first step to understand any plausible correlations between the thermodynamically stable configurations and the electrochemically induced structural changes. If two Li-intercalation oxides have an identical crystal structure and compositional analogy, it is usually accepted that their point defect structure and distribution will not completely differ in general. Contrary to this conventional wisdom, our STEM observations along with the DFT calculations in this study demonstrate an important discrimination in the distribution of antisite exchange defects between two identical olivine frameworks, LiMnPO4 and LiFePO4, highlighting the value of an atomic column-by-column analysis.

Conclusion Using HAADF-STEM, we have demonstrated the distinct distribution characteristics of antisite exchange defects between LiMnPO4 and LiFePO4 in the surface regions. While strong subsurface segregation of the defects was observable in some selective surface orientations in LiFePO4, no substantial differences in the highly ordered Li−Mn arrangement were identified near surface facets during the STEM analysis. Ab initio DFT calculations also showed that the formation energies of a pair of

ACS Paragon Plus Environment

Page 9 of 17

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

The Journal of Physical Chemistry 9

exchange defects do not vary significantly with the surface crystallographic orientation in LiMnPO4, consistently providing theoretical support to the experimental observations. This rapid communication revealing the unusual discrimination of defect distribution in two olivine phosphates emphasizes the necessity of direct atomic-scale probing for understanding the unexpected defect chemistry in Li-intercalation compounds.

Associated Content

Supporting Information The Supporting Information is available free of charge on the ACS Publication website at DOI: Unit cell structure of LiMnPO4; intensity profiles of Li columns; and additional table for defect formation energies (PDF)

Author Information

Corresponding Author *E-mail: [email protected] or [email protected] Author Contributions 4 S.-Y. Chung and S.-Y. Choi contributed equally to this work. Notes The authors declare no competing financial interest.

Acknowledgements

This work was supported by the National Research Foundation of Korea (NRF), grant numbers 2014R1A4A1003712 (BRL Program) and 2012M2A2A6002461 (Nuclear R&D Program) and also by the Ministry of Tarde, Industry, and Energy, grant number 10065691.

References 1. Lu, Z.; Beaulieu, L. Y.; Donaberger, R. A.; Thomas, C. L.; Dahn, J. R. Synthesis, Structure, and Electrochemical Behavior of Li[NixLi1/3−2x/3Mn2/3−x/3]O2. J. Electrochem. Soc. 2002, 149, A778−A791. 2. Shao-Horn, Y.; Levasseur, S.; Weill, F.; Delmas, C. Probing Lithium and Vacancy Ordering; An Electron Diffraction Study. J. Electrochem. Soc. 2003, 150, A366−A373. 3. Kang, K.; Meng, Y. S.; Bréger, J.; Grey, C. P.; Ceder, G. Electrodes with High Power and High Capacity for Rechargeable Lithium Batteries. Science 2006, 311, 977−980. 4. Chung, S.-Y.; Choi, S.-Y.; Yamamoto, T.; Ikuhara, Y. Orientation-Dependent

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 10 of 17 10

Arrangement of Antisite Defects in Lithium Iron(II) Phosphate Crystals. Angew. Chem. Int. Ed. 2009, 48, 543−546. 5. Xu, B.; Fell, C. R.; Chi, M.; Meng, Y. S. Identifying Surface Structure Changes in Layered Li-Excess Nickel Manganese Oxides in High Voltage Lithium Ion Batteries: A Joint Experimental and Theoretical Study. Energy Environ. Sci. 2011, 4, 2223−2233. 6. Javis, K. A.; Deng, Z.; Allard, L. F.; Manthiram, A.; Ferreira, P. J. Atomic Structure of a Lithium-Rich Layered Oxide Material for Lithium-Ion Batteries: Evidence of a Solid Solution. Chem. Mater. 2011, 23, 3614−3621. 7. Hoang, K.; Johannes, M. Tailoring Native Defects in LiFePO4: Insights from FirstPrinciples Calculations. Chem. Mater. 2011, 23, 3003−3013. 8. Gu, M.; Belharouak, I.; Genc, A.; Wang, Z.; Wang, D.; Amine, K.; Gao, F.; Zhou, G.; Thevuthasan, S.; Baer, D. R.; et al. Conflicting Roles of Nickel in Controlling Cathode Performance in Lithium Ion Batteries. Nano Lett. 2012, 12, 5186−5191. 9. Lee, S.; Yoon, G.; Jeong, M.; Lee, M.; Kang, K.; Cho, J. Hierarchical Surface Atomic Structure of a Manganese-Based Spinel Cathode for Lithium-Ion Batteries. Angew. Chem. Int. Ed. 2015, 54, 1153−1158. 10. Chung, S.-Y.; Choi, S.-Y.; Yamamoto, T.; Ikuhara, Y. Atomic-Scale Visualization of Antisite Defects in LiFePO4. Phys. Rev. Lett. 2008, 100, 125502. 11. Lee, J.; Zhou, W.; Idrobo, J. C.; Pennycook, S. J.; Patelides, S. T. Vacancy-Driven Anisotropic Defect Distribution in the Battery-Cathode Material LiFePO4. Phys. Rev. Lett. 2011, 107, 085507. 12. Chung, S.-Y.; Choi, S.-Y.; Lee, S.; Ikuhara, Y. Distinct Configurations of Antisite Defects in Ordered Metal Phosphates: Comparison between LiMnPO4 and LiFePO4. Phys. Rev. Lett. 2012, 108, 195501. 13. Guo, X.; Wang, M.; Huang, X.; Zhao, P.; Liu, X.; Che, R. Direct evidence of Antisite defects in LiFe0.5Mn0.5PO4 via Atomic-Level HAADF-EELS. J. Mater. Chem. A, 2013, 1, 8775−8781. 14. Truong, Q. D.; Devaraju, M. K.; Sasaki, Y.; Hyodo, H.; Tomai, T.; Honma, I. Relocation of Cobalt Ions in Electrochemically Delithiated LiCoPO4 Cathode Materials. Chem. Mater. 2014, 26, 2770−2773. 15. Lin, F.; Markus, I. M.; Nordlund, D.; Weng, T.-C.; Asta, M. D.; Xin, H. L.; Doeff, M. M. Surface Reconstruction and Chemical Evolution of Stoichiometric Layered Cathode Materials for Lithium-Ion Batteries. Nat. Commun. 2014, 5, 3529. 16. Dixit, H.; Zhou, W.; Idrobo, J.-C.; Nanda, J.; Cooper, V. R. Facet-Dependent Disorder in Pristine High-Voltage Lithium−Manganese-Rich Cathode Material. ACS Nano 2014, 8, 12710−12716.

ACS Paragon Plus Environment

Page 11 of 17

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

The Journal of Physical Chemistry 11

17. Yan, P.; Nie, A.; Zheng, J.; Zhou, Y.; Lu, D.; Zhang, X.; Xu, R.; Belharouak, I.; Zu, X.; Xiao, J.; et al. Evolution of Lattice Structure and Chemical Composition of the Surface Reconstruction Layer in Li1.2Ni0.2Mn0.6O2 Cathode Material for Lithium Ion Batteries. Nano Lett. 2015, 15, 514−522. 18. Yan, P.; Zheng, J.; Zheng, J.; Wang, Z.; Teng, G.; Kuppan, S.; Xiao, J.; Chen, G.; Pan, F.; Zhang, J.-G.; et al. Ni and Co Segregation on Selective Surface Facets and Rational Design of Layered Lithium Transition-Metal Oxide Cathodes. Adv. Energy Mater. 2016, 6, 1502455. 19. Gu. L.; Zhu, C.; Li, H.; Yu, Y.; Li, C.; Tsukimoto, S.; Maier, J.; Ikuhara, Y. Direct Observation of Lithium Staging in Partially Delithiated LiFePO4 at Atomic Resolution. J. Am. Chem. Soc. 2011, 133, 4661−4663. 20. Gu, M.; Belharouak, I.; Zheng, J.; Wu, H.; Xiao, J.; Genc, A.; Amine, K.; Thevuthasan, S.; Baer, D. R.; Zhang, J.-G.; et al. Formation of the Spinel Phase in the Layered Composite Used in Li-Ion Batteries. ACS Nano 2013, 7, 760−767 (2013). 21. Tang, D.; Sun, Y.; Yang, Z.; Ben, L.; Gu, L.; Huang, X.. Surface Structure Evolution of LiMn2O4 Cathode Material upon Charge/Discharge. Chem. Mater. 2014, 26, 3535−3543. 22. Lu, P.; Yan, P.; Romero, E.; Spoerke, E. D.; Zhang, J.-G.; Wang, C.-M. Observation of Electron-Beam-Induced Phase Evolution Mimicking the Effect of the Charge−Discharge Cycle in Li-Rich Layered Cathode Materials Used for Li Ion Batteries. Chem. Mater. 2015, 27, 1375−1380. 23. Ryoo, H.; Bae, H. B.; Kim, Y.-M.; Kim, J.-G.; Lee, S.; Chung, S.-Y. Frenkel-DefectMediated Chemical Ordering Transition in a Li−Mn−Ni Spinel Oxide. Angew. Chem. Int. Ed. 2015, 54, 7963−7967. 24. Lu, X.; Gu, L.; Hu, Y. S.; Chiu, H. C.; Li, H.; Demopoulos, G. P.; Chen, L. New Insight into the Atomic-Scale Bulk and Surface Structure Evolution of Li4Ti5O12 Anode. J. Am. Chem. Soc. 2015, 137, 1581−1586. 25. Chung, S.-Y.; Choi, S.-Y.; Kim, T.-H.; Lee, S. Surface-Orientation-Dependent Distribution of Subsurface Cation-Exchange Defects in Olivine-Phosphate Nanocrystals. ACS Nano 2015, 9, 850−859. 26. Chung, S.-Y.; Kim, Y.-M.; Lee, S.; Oh, S. H.; Kim, J.-G.; Choi, S.-Y.; Kim, Y.-J.; Kang, S.-J. L. Cation Disordering by Rapid Crystal Growth in Olivine-Phosphate Nanocrystals. Nano Lett. 2012, 12, 3068−3073. 27. Wang, L.; Zhou, F.; Meng, Y. S.; Ceder, G. First-Principles Study of Surface Properties of LiFePO4: Surface Energy, Structure, Wulff Shape, and Surface Redox Potential. Phys. Rev. B 2007, 76, 165435.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 12 of 17 12

Figure 1. HAADF-STEM images showing the (001) vicinal surface regions in LiFePO4 and LiMnPO4 along with their unit-cell structures. This set of two atomic-column images directly compares the distinct distribution of subsurface exchange defects. (a) Most Li columns (green arrows) in LiFePO4 reveal detectable bright intensity, indicating a large amount of Li−Fe exchange defects near the surface. (b) In contrast, few Li columns containing antisite Mn are identified in LiMnPO4. Schematic illustrations for atomic arrangement are superimposed on each of the images.

ACS Paragon Plus Environment

Page 13 of 17

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

The Journal of Physical Chemistry 13

Figure 2. HAADF-STEM image showing a Li column with detectable intensity (denoted by a yellow arrow) due to antisite Mn in the subsurface region of LiMnPO4. The intensity profile in the right panel consistently verifies the higher contrast in Column 4 (yellow arrow) among Columns 1−4.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 14 of 17 14

Figure 3. Series of HAADF-STEM images with different surface orientations in a LiMnPO4 particle. Three atomic-column images in the vicinal (001), (100) and (10 1 ) surfaces were obtained from Regions 1−3 denoted by orange rectangles in the lowmagnification image. The highly ordered feature between the bright Mn columns and the invisible Li columns is maintained well over all three images, irrespective of the surface orientation.

ACS Paragon Plus Environment

Page 15 of 17

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

The Journal of Physical Chemistry 15

Figure 4. Supercell models for DFT calculations of subsurface defect formation energies in LiMnPO4. A pair of antisite exchange defects, LiMn−MnLi, is denoted by a red ellipse in each of the illustrations. Note that this defect pair has been introduced in the subsurface layer rather than on the top-most surface so as not to disturb the initial atomic termination in each slab.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 17 16

Table 1. Formation energies of an antisite defect pair in each subsurface layer and differences with its formation energy in the bulk. Surface Plane

Defect Formation Energy (ES )*, eV LiMnPO4

LiFePO4

(100)

0.75

0.59

(010)

0.39

0.14

(001)

0.41

−0.56

(110)

0.76

0.15

(101)

0.64

0.38

(011)

0.43

0.25

(111)

0.50

−0.70

(*Defect Formation Energy of LiMnPO4 in Bulk, EB : 0.75 eV. The formation energy in red in LiFePO4 denotes a larger value than the formation energy in bulk (0.43 eV). The negative formation energies are also represented in blue.)

ACS Paragon Plus Environment

Page 17 of 17

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

The Journal of Physical Chemistry 17

TOC Graphic

ACS Paragon Plus Environment