Skip to main content

Research Article

Synergistic O2-/Li+ Dual Ion Transportation at Atomic Scale

F. Q. Meng1, Q. H. Zhang1,*, A. Gao2, X. Z. Liu1, J. N. Zhang1, S. Y. Peng3, X. Lu2,*, L. Gu1,4,5,*, and H. Li1,4,5

1Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
2School of Materials, Sun Yat-Sen University, Guangzhou 510275, China
3College of Materials Science and Engineering, Jiangxi University of Science and Technology, Jiangxi 341000, China
4Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China
5School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
*Corresponding author(s) email(s): and nc.ude.usys.liam@3aixul and


The ion migration during electrochemical process is a fundamental scientific issue for phase transition behavior and of technical importance for various functional devices, where cations or anions are active under electrical bias. Usually only one type of functional ion, O2- or Li+, is activated due to their different migration energy barriers, cooperated by the valence change of other immobile ions in the host lattice matrix, e.g., Co2+/Co3+ and Mn3+/Mn4+ redox couples, owing to the charge neutralization. Here we select spinel Li4Ti5O12 as anode and construct an all-solid-state battery under a transmission electron microscope; a synergistic transportation of O2- and Li+ driven by an electrical bias was directly observed at the atomic scale. A small amount of oxygen anions was extracted firstly as a result of its lowest vacancy formation energy under 2.2 V, leading to the vertical displacement of oxygen. Up to 2.7 V, an ordered phase with both Li- and O- deficiency formed. The Li+ and O2- ions are simultaneously extracted out from the [LiO4] tetrahedra due to the electroneutrality principle. The migration paths of O and Li have been proposed and verified by first-principles calculations. These results reveal a brand new synergistic ion migration manner and may provide up-to-date insights on the transportation process of lithium ion conductors.

1. Introduction

Electrochemical phase transition involves various migration processes of ions, which is the foundation of its functionality, especially in the operation of lithium-ion battery (LIB) materials where the Li ions shuffle between anode and cathode during charge-discharge process [13]. Usually only one type of functional ions responses to the external electric field, cooperated by the valence change of other immobile ions in the host lattice matrix due to the charge neutralization. For example, the extraction and insertion of lithium are accompanied by the oxidation and reduction of the transition metal ions, Co2+/Co3+ and Mn3+/Mn4+ redox couples, in typical anode materials LiCoO2 and LiMn2O4, respectively [4, 5]. However, the situation of lithium extraction without valence-changeable cations, where the charge compensation cannot be accomplished by the oxidation of transition metal ions, is a less investigated area but deserves deep understandings. At this point, we demonstrate the detailed migration process of Li+ and O2- when we further extracted the lithium from the stable spinel Li4Ti5O12 (LTO) with the chemically stabilized Ti4+ ions upon delithiation.

Crystallographically, LTO spinel structure can be expressed as [Li [Li1/3Ti5/3 [O4 as displayed in Figure 1(a), where one Li+ ion occupies the 8a sites with tetrahedral coordination, the remaining 1/3 Li+ and all Ti4+ ions being in the 16d sites with octahedral coordination and the O2- ions in the 32e sites. The [TiO6] octahedra share faces with the [LiO4] tetrahedra. There are also empty tetrahedral 8a and octahedral 16c sites interconnected by edge-sharing octahedral sites, providing the interstitial space for the ion migration. In fact, the Li4Ti5O12 is mostly used as an anode material in lithium ion batteries (LIBs) [6]. During lithiation, the Coulombic interaction between tetrahedral sites drives the migration of Li+ ions from the tetrahedral 8a site to the empty octahedral 16c sites with a small energy barrier ~ 0.3 eV, resulting in the Li7Ti5O12 phase with rock-salt structure in two-phase separation mechanism [79] and the partial reduction of Ti4+ to Ti3+ ions [7, 1012]. During the delithiation, the reverse ion migration process and phase transition from the rock-salt Li7Ti5O12 to the spinel Li4Ti5O12 occur, accompanying the oxidation of partial Ti3+ to tetravalent Ti ions. However, in practical operation of LIBs with the Li4Ti5O12 anode, the over-discharging is always unavoidable, which is directly correlating with the capacity fading upon cycling. The charge compensation for the further extraction of Li ions from the [Ti4+]5 seems significantly difficult from the Ti4+ ions due to its stabilized electron configuration, thus promising that from the lattice oxygen ions upon the Li extraction. However, it is unclear whether the oxygen ions would be oxidized to the O- state or escaped from the Li4Ti5O12 lattice as the O2 gas. Therefore, it is desirable to reveal the structure evolution under the electric extraction of Li ions to understand the degradation mechanism of Li4Ti5O12 spinel.

Figure 1

Alt text
Atomic structure of Li4Ti5O12 and its in-situ experiment setup. (a) Structural schematics of Li4Ti5O12 and the projection along the <110> direction. Ti columns with different projection density are labeled with Ti-a and Ti-b, the neighbored oxygen columns are also labeled with Oa and Ob, respectively. (b) In-situ setup for constructed batteries where Li4Ti5O12 serves as positive electrode, the LLZO acts as the electrolyte and Si as the counter-electrode. (c, d) HAADF and ABF images of Li4Ti5O12 viewed along the <110> direction, the Ti, O and Li columns have been identified at atomic scale. Scale bar: 1 nm.

2. Results

2.1. In Situ Setup for Ion Extraction

We constructed an all-solid-sate battery in the TEM by a newly developed chip-based in situ technique with focused ion beam (FIB) milling [13], which allows the atomic resolution imaging under electric bias, and investigated the structural evolution of Li4Ti5O12 during the lithium extraction. As shown in Figure 1(b), we used the LLZO as the electrolyte and Si as the counterelectrode and applied a positive bias on the Li4Ti5O12 electrode. Figure 1(c) demonstrates the high-angle dark-field (HAADF) image of Li4Ti5O12 viewed along the direction under this electrical configuration, where Ti columns are imaged as the strong white spots (Ti-a arrays) and weak white spots (Ti-b arrays) because the Ti content in Ti-a columns is roughly twice that in the Ti-b columns as indicated by red arrows, while the O and Li are hardly visualized due to the Z2 dependence of HAADF contrast (Z is elemental number). On the other hand, the angular bright-field (ABF) imaging geometry is extremely feasible to image light elements such as lithium, where the image contrast is of ~Z1/3 dependence [14, 15]. Consequently, the Li and O columns can be observed clearly in the ABF image as displayed in Figure 1(d). Here we define the oxygen in the Ti-a array as Oa and that in the Ti-b arrays as Ob. We can see that the vertical atomic displacement in the Oa arrays is more prominent than that in the Ob arrays within the Ti (001) plane, which is in good agreement with the pristine structure shown in Figure 1(a). It should be noted that although the Z of Li is smaller than that of O, the ABF contrast of Li shown here is comparable to that of O due to the thickness effect [16].

2.2. Formation of Oxygen Vacancies and Related Atomic Shift at 2.2 V

When an appropriate positive bias was applied on the top electrode, in LTO electrode the Li+ ions would move towards LLZO and O2- ions towards top electrode, as indicated by arrows in Figure 1(b). Then the voltage was increased from the 0 V to 2.2 V, where little change was detected from the Ti diamond of the spinel lattice as shown in the HAADF image (Figures 2(b) and 2(d)), while the vertical displacement of oxygen at the Oa sites becomes more obvious with respect to the vertical displacement in the Ob arrays decrease and lead to a nearly flat arrangement (Figures 2(a) and 2(c)) as depicted in the ABF images. We schematize this change in Figure 2(e), where the purple and dark-yellow arrows indicate the vertical displacement of the Oa and the Ob arrays under 2.2V, with reference to the pristine Li4Ti5O12 lattice. This atom transfer of oxygen may be induced by the presence of some oxygen vacancies—interstitial oxygen cluster as reported by Kanno et al. [17]. Meanwhile, the Li ions at 8a sites slightly shift away from their original positions due to the presence of O vacancies, similar to the situation in defected spinel Mn2O4-δ [18]. Thus, the atomic displacement of oxygen indicates the forever n of a certain amount of oxygen vacancies under a modest (2.2 V) electric bias.

Figure 2

Alt text
STEM images of Li4-x1Ti5O12-y1 under 2.2 V. (a, b) HAADF and ABF images under 2.2 V, respectively. (c, d) the enlarged HAADF and ABF images from the red dotted frame in (a, b). Blue arrows in ABF image (c) indicated the oxygen atomic shift. Schematics of oxygen atomic shift (e), brown and purple arrows are used to label the displacement direction in the Ti-a and Ti-b columns, respectively. Scale bar: 1 nm.

To find out which ion is more likely to be extracted firstly, we compared the vacancy formation energy at the Li- and O- sites in spinel Li4Ti5O12 by first-principles calculations (simulation details refer to our previous work [7, 9, 19]). The results are 4.88 eV, 5.10 eV, and 4.62 eV for the Li-Vacancy ( ) at 8a sites, at 16d sites, and O-Vacancy ( ) at 32e sites. Thus, the VO would be created firstly under electric bias due to its smallest formation energy, while the Li removal (8a site) will almost be accompanied with the O loss. This is in agreement with our observation of oxygen shift under 2.2 V due to the oxygen defects. Thus, the chemical formula of the defected spinel LTO can be depicted as Li4-xTi5O12-δ. Here the x indicates that some Li ions also should be extracted out from the lattice due to the charge compensation mechanism.

2.3. Li Extraction and Oxygen Vacancy Ordering at 2.7V

An obvious structural change emerges at 2.7V. As shown in Figure 3(a), a bright contrast in the inner of the Ti diamond appears in the ABF image, indicating lots of vacancies in the 8a-Li sites. From the enlarged ABF image in Figure 3(c), we find that the contrast of the Li ions at the upper 8a sites (outlined by a small green dotted circle) disappears, indicating the partial extraction of Li ions from the defected spinel under 2.7 V. The rest of the Li ions deviated from the lower 8a site and shifted towards one of the Oa sites (outlined by a small red dotted circle). Meanwhile, the oxygen at Ob sites also exhibits periodic change in ABF intensity as indicated by the yellow and small brown arrows in Figure 3(c), where an obvious decision decrease of oxygen occupation can be visualized in the Ob sites labeled by small brown arrows. Line profile of the ABF image contrast along the [Ti-b]-[Ob]-[Ti-b] chain (Figure 3(d)) from the positions indicated by the yellow arrows in Figure 3(a) also demonstrated an obvious periodic intensity modulation (yellow and small brown arrows) at the Ob sites, which reflects the occupancy variation of oxygen and formation of ordered oxygen vacancies. For clarity, O1 and O2 are used to distinguish these two types of Ob sites as labeled in Figure 3(d). It means that accompanying the Li extraction from one of the 8a sites, the oxygen in one of the Ob sites simultaneously migrated out to satisfy the charge compensation and resulted in an oxygen-vacancy ordered Li-deficient structure when the amount of the oxygen vacancy reached a threshold. Line scan EEL spectra have also been acquired along the [Ti-b]-[Ob]-[Ti-b] chain and compared with the pristine state of LTO at two distinct oxygen sites O1 and O2 can be identified. No obvious changes in the Ti L edges were observed, indicating the valence state invariance of Ti4+, as shown in Figure 3(e), while obvious decrease of peaks b, c, and d in the O K edges can be visualized in Figure 3(f) which manifests the presence of oxygen vacancies [20]. More importantly, the decrease of b, c, and d peak at the O1 sites is more prominent than that at the O2 sites, which is in good agreement with the modulation of oxygen content shown in ABF image (Figure 3(a)).

Figure 3

Alt text
STEM images and EEL spectra of Li4-x2Ti5O12-y2 under 2.7 V. (a) ABF and (b) HAADF images of the oxygen-vacancy ordered Li-deficient LTO. (c) Enlarged ABF image from the red frame in (a). (d) Line profile of inversed ABF contrast along the [Ti-b]-[Ob]-[Ti-b] chain. (e) Ti L edges and (f) O K edge of O1 sites and O2 sites, which are normalized by peak a. Scale bar: 1 nm.

3. Discussion

3.1. Discussion on Ion Migration Path

Then a 1×3×1 LTO supercell was built to simulate the Li and O transport features (details on the structural construction and Li-Ti 16d antisite optimization refer to our previous work [14]). Figure 4 shows the Li, O, and Li-O concerted diffusion inside the LTO supercell, respectively. Upon the formation of O vacancy, the simulation results demonstrated that one VO introduces obvious lattice distortion, leading to one five-coordinated [LiO5] polyhedra surrounded by three [TiO6] octahedra, by which the local structural distortion was generated to modulate the Li and/or O migration behavior [9]. From the STEM images as shown in Figures 2 and 3, the contrast decrease in O columns becomes increasingly obvious, indicating that the O ions inside the LTO spinel are gradually extracted out at a higher electric bias. From a thermodynamic point of view, the mobile O ions would probably interact with each other to release the O2 gas. At the same time, the delithiation process occurs to establish the charge balance. In this respect, the extraction of Li and O ions in LTO spinel could possibly happen in a synergetic way.

Figure 4

Alt text
Theoretical calculations on ion diffusion. The Li (a), O (b) Li-O (c) concerted diffusion inside a 1×3×1 LTO supercell with the corresponding activation energies (d).

Figures 4(a) and 4(b) manifest the single Li or O ions diffusion inside the LTO spinel with their activation energy shown in Figure 4(d). Upon a Li ion migration, the standard 8a-16c-8a 3D diffusion feature is derived from simulations. The activation energy of Li ion diffusion is estimated to be 0.48 eV, while the one step hopping of O ions needs to overcome an energy barrier of 1.33 eV, indicating the mobile possibility at such higher electric bias. Then if the synergetic migration of Li and O ions inside the LTO spinel is considered, the total energy barrier is approx. 1.48 eV, a little higher than the single O ion migration, which means the possible interaction between the Li and O ions during the dynamic transport process. This also accounts for the STEM observations in Figures 3 and 4. Further simulation results indicated that there is a higher activation energy for the Li-O dual ion migration process if the Li and O trajectories are very close with each other and the diffusion trajectory will be much more intricate as well.

Similar synergistic transport behavior of dual Li and O ions could be extended to other functional oxides without valence-changeable cations, where extra capacity can be expected due to the extra Li extraction. Contrary to O2 gas release process in other Li- and Mn-rich cathode materials [21], the oxygen loss here does not introduce obvious reduction of Ti and structural degradation, as shown in Figures 3(b) and 3(e). In addition, a shrinking of 1.25% along the short diagonal and an expansion of 4.67% along the long diagonal of the Ti diamond were obtained from the HAADF images under 2.7 V, respectively. This relative low lattice strain should be related to the simultaneous extraction of dual Li and O ions where the electrostatic repulsive force due to the Li extraction between O2- layers is cancelled to a certain extent, promising a good stability during the synergistic transport process.

3.2. Conclusion

The synergistic transport processes of O2- and Li+ ions under electric field are observed directly at the atomic scale, and an ordered oxygen deficient phase was identified subsequently at a higher voltage. A small amount of oxygen anions was extracted firstly and led to the vertical displacement of oxygen under 2.2 V; an ordered phase with both Li- and O- deficiency formed up to 2.7 V. The migration path of Li+  and O2- ions has been discussed by first-principle calculations. These results reveal that the Li+ and O2- ions are simultaneously extracted out from the [LiO4] tetrahedral, demonstrating a new synergistic migration manner of functional ions in lithium ion conductors.

4. Materials and Methods

4.1. STEM Experiments

High-resolution STEM investigation was carried out on a JEM-ARM200F microscope with double correctors for the condenser lens and objective lens. The available point resolution is better than 0.08 nm at an operating voltage of 200 kV. HAADF and ABF images were acquired at acceptance angles of 70-150 mrad and 11-22 mrad, respectively. All of the images presented here are Fourier-filtered to minimize the effect of the contrast noise. The filtering does not have any effect on the results of our measurements. In addition, we adjusted deliberately the brightness and contrast in order to better represent the atomic arrangements; this adjustment does not have any effect on the result.

4.2. DFT Calculations

Spin-polarized calculations were performed using the Vienna Ab Initio Simulation Package within the projector augmented-wave approach [2224]. Generalized gradient approximation in the parameterization of Perdew, Burke, and Ernzerh was used to describe the exchange-correlation potential [2527]. The cutoff of the kinetic energy was set to 600 eV for all calculations. Geometry optimization was performed using a 1a 3b 1c supercell and considered converged when the force on each atom was <1 meV Å−1. The Brillouin zone integration was performed with 3 1 3 and 1 1 1 Γ centered Monkhorst-Pack k-point meshes [28] in geometry optimization and climbing image nudged elastic band calculations, respectively.

Data Availability

All data needed to evaluate the conclusions in the paper are present in the paper. Additional data related to this paper may be requested from the authors (QH Zhang:

Conflicts of Interest

The authors declare no competing financial interests.

Authors’ Contributions

L. Gu and H. Li conceived this project. F. Q. Meng, X. Z. Liu, and Q. H. Zhang designed experimental procedures in detail. J. N. Zhang provided the LLZO sample. F. Q. Meng, S. Y. Peng, and Q. H. Zhang performed the characterization and data analysis. The first principle calculations were performed by A. Gao under the guidance of X. Lu. All of the authors participated in the preparation of the manuscript.


This work was supported by National Program on Key Basic Research Project (2014CB921002, 2015CB921700) and The Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDB07030200) and National Natural Science Foundation of China (51522212, 51421002, 51332001, and 11704019).

References and notes

  1. J. B. Goodenough and K.-S. Park, “The Li-ion rechargeable battery: a perspective,” Journal of the American Chemical Society, vol. 135, no. 4, pp. 1167–1176, 2013. View at Publisher · View at Scopus · View at Google Scholar
  2. J. B. Goodenough and Y. Kim, “Challenges for rechargeable Li batteries,” Chemistry of Materials, vol. 22, no. 3, pp. 587–603, 2010. View at Publisher · View at Scopus · View at Google Scholar
  3. T. Shang, Y. Wen, D. Xiao, L. Gu, Y. Hu, and H. Li, “Atomic-scale monitoring of electrode materials in lithium-ion batteries using in situ transmission electron microscopy,” Advanced Energy Materials, vol. 7, no. 23, p. 16, 2017. View at Publisher · View at Google Scholar
  4. T. Ohzuku and A. Ueda, “Solid-state redox reactions of LiCoO2 (R-3m) for 4 volt secondary lithium cells,” Journal of The Electrochemical Society, vol. 141, no. 11, pp. 2972–2977, 1994. View at Publisher · View at Scopus · View at Google Scholar
  5. A. R. Armstrong and P. G. Bruce, “Synthesis of layered LiMnO2 as an electrode for rechargeable lithium batteries,” Nature, vol. 381, no. 6582, pp. 499-500, 1996. View at Publisher · View at Scopus · View at Google Scholar
  6. L. Aldon, P. Kubiak, M. Womes et al., “Chemical and electrochemical Li-insertion into the Li4Ti5O12 spinel,” Chemistry of Materials, vol. 16, no. 26, pp. 5721–5725, 2004. View at Publisher · View at Scopus · View at Google Scholar
  7. X. Lu, L. Zhao, X. He et al., “Lithium storage in Li4Ti5O12 spinel: the full static picture from electron microscopy,” Advanced Materials, vol. 24, no. 24, pp. 3233–3238, 2012. View at Publisher · View at Scopus · View at Google Scholar
  8. X. Lu, L. Gu, Y.-S. Hu et al., “New insight into the atomic-scale bulk and surface structure evolution of Li4Ti5O12 anode,” Journal of the American Chemical Society, vol. 137, no. 4, pp. 1581–1586, 2015. View at Publisher · View at Scopus · View at Google Scholar
  9. Y. Sun, L. Zhao, H. Pan et al., “Direct atomic-scale confirmation of three-phase storage mechanism in Li4Ti5O12 anodes for room-temperature sodium-ion batteries,” Nature Communications, vol. 4, 2013. View at Scopus · View at Google Scholar
  10. H.-C. Chiu, X. Lu, J. Zhou et al., “Annealing-regulated elimination of residual strain-induced structural relaxation for stable high-power Li4Ti5O12 nanosheet anodes,” Nano Energy, vol. 32, pp. 533–541, 2017. View at Publisher · View at Scopus · View at Google Scholar
  11. H. Chiu, X. Lu, J. Zhou et al., “Capacity fade mechanism of Li4Ti5O12 nanosheet anode,” Advanced Energy Materials, vol. 7, no. 5, p. 1601825, 2017. View at Publisher · View at Google Scholar
  12. Y. Wen, X. Chen, X. Lu, and L. Gu, “Interface charges boosted ultrafast lithiation in Li4Ti5O12 revealed by in-situ electron holography,” Journal of Energy Chemistry, 2018. View at Scopus · View at Google Scholar
  13. Y. Gong, J. Zhang, L. Jiang et al., “In situ atomic-scale observation of electrochemical delithiation induced structure evolution of LiCoO2 cathode in a working all-solid-state battery,” Journal of the American Chemical Society, vol. 139, no. 12, pp. 4274–4277, 2017. View at Publisher · View at Scopus · View at Google Scholar
  14. S. D. Findlay, N. Shibata, H. Sawada et al., “Robust atomic resolution imaging of light elements using scanning transmission electron microscopy,” Applied Physics Letters, vol. 95, no. 19, p. 191913, 2009. View at Publisher · View at Google Scholar
  15. S. Kobayashi, S. D. Findlay, N. Shibata et al., “Simultaneous visualization of oxygen vacancies and the accompanying cation shifts in a perovskite oxide by combining annular imaging techniques,” Applied Physics Letters, vol. 100, no. 19, p. 193112, 2012. View at Publisher · View at Google Scholar
  16. X. He, L. Gu, C. Zhu et al., “Direct imaging of lithium ions using aberration-corrected annular-bright-field scanning transmission electron microscopy and Associated contrast Mechanisms,” Materials Express, vol. 1, no. 1, pp. 43–50, 2011. View at Publisher · View at Scopus · View at Google Scholar
  17. R. Kanno, A. Kondo, M. Yonemura et al., “The relationships between phases and structures of lithium manganese spinels,” Journal of Power Sources, vol. 81-82, pp. 542–546, 1999. View at Publisher · View at Scopus · View at Google Scholar
  18. R. Huang, Y. H. Ikuhara, T. Mizoguchi et al., “Oxygen-vacancy ordering at surfaces of lithium manganese(III,IV) oxide spinel nanoparticles,” Angewandte Chemie International Edition, vol. 50, no. 13, pp. 3053–3057, 2011. View at Publisher · View at Scopus · View at Google Scholar
  19. Q. Zhang, S. Zhang, F. Ning et al., “Calcium doping of lithium titanium oxide nanospheres: a combined first-principles and experimental study,” Energy Technology, vol. 5, no. 4, pp. 539–543, 2017. View at Publisher · View at Scopus · View at Google Scholar
  20. D. A. Muller, N. Nakagawa, A. Ohtomo, J. L. Grazul, and H. Y. Hwang, “Atomic-scale imaging of nanoengineered oxygen vacancy profiles in SrTiO3,” Nature, vol. 430, no. 7000, pp. 657–661, 2004. View at Publisher · View at Scopus · View at Google Scholar
  21. E. Hu, X. Yu, R. Lin et al., “Evolution of redox couples in Li- and Mn-rich cathode materials and mitigation of voltage fade by reducing oxygen release,” Nature Energy, vol. 3, no. 8, pp. 690–698, 2018. View at Publisher · View at Google Scholar
  22. G. Mills and H. Jónsson, “Quantum and thermal effects in H2 dissociative adsorption: Evaluation of free energy barriers in multidimensional quantum systems,” Physical Review Letters, vol. 72, no. 7, pp. 1124–1127, 1994. View at Publisher · View at Scopus · View at Google Scholar
  23. G. Kresse and J. Hafner, “Ab initio molecular dynamics for liquid metals,” Physical Review B: Condensed Matter and Materials Physics, vol. 47, no. 1, pp. 558–561, 1993. View at Publisher · View at Scopus · View at Google Scholar
  24. G. Kresse and J. Furthmüller, “Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set,” Physical Review B: Condensed Matter and Materials Physics, vol. 54, no. 16, pp. 11169–11186, 1996. View at Publisher · View at Scopus · View at Google Scholar
  25. J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996. View at Publisher · View at Scopus · View at Google Scholar
  26. P. E. Blöchl, “Projector augmented-wave method,” Physical Review B: Condensed Matter and Materials Physics, vol. 50, no. 24, pp. 17953–17979, 1994. View at Publisher · View at Scopus · View at Google Scholar
  27. G. Kresse and D. Joubert, “From ultrasoft pseudopotentials to the projector augmented-wave method,” Physical Review B: Condensed Matter and Materials Physics, vol. 59, no. 3, pp. 1758–1775, 1999. View at Publisher · View at Scopus · View at Google Scholar
  28. H. J. Monkhorst and J. D. Pack, “Special points for Brillouin-zone integrations,” Physical Review B: Condensed Matter and Materials Physics, vol. 13, no. 12, pp. 5188–5192, 1976. View at Publisher · View at MathSciNet · View at Scopus · View at Google Scholar

2019 © F. Q. Meng et al. Exclusive License Science and Technology Review Publishing House. Distributed under a Creative Commons Attribution License (CC BY 4.0).