搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

富锂锰基三元材料Li1.208Ni0.333Co0.042Mn0.417O2的电子结构和缺陷性质

黄文军 王亚平 曹昕睿 吴顺情 朱梓忠

引用本文:
Citation:

富锂锰基三元材料Li1.208Ni0.333Co0.042Mn0.417O2的电子结构和缺陷性质

黄文军, 王亚平, 曹昕睿, 吴顺情, 朱梓忠

Electronic structures and defect properties of lithium-rich manganese-based ternary material Li1.208Ni0.333Co0.042Mn0.417O2

Huang Wen-Jun, Wang Ya-Ping, Cao Xin-Rui, Wu Shun-Qing, Zhu Zi-Zhong
PDF
HTML
导出引用
  • 采用基于密度泛函理论的第一性原理方法, 计算了$ R\bar{3}m $相的富锂锰基三元正极材料Li1.208Ni0.333Co0.042Mn0.417O2的晶体结构、电子结构以及缺陷性质. 结果表明, Li1.208Ni0.333Co0.042Mn0.417O2是一种具有直接带隙的磁性半导体材料, 其价带顶是O的px, py, pz轨道与Ni的dxy, dyz, dxz原子轨道之间的杂化, 导带底除了有价带顶的特征外, 还有部分的Ni-$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $和Mn-$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $, Mn-$ {3\mathrm{d}}_{yz} $ 轨道参与杂化. 差分电荷密度图显示, 金属原子与O原子之间的成键方式是共价键和离子键的混合. 本文还计算了脱去单个过渡金属原子的空位形成能. 脱去Mn原子时的缺陷体系的体积变化最大, 而脱去Co原子的缺陷体系的体积则几乎不变. 得到的金属空位形成能的大小依次为Ef (Mn-空位) > Ef (Co-空位) > Ef (Ni-空位). 差分电荷密度图还显示, 空位的产生仅强烈影响了空位附近的部分O原子, 体现了金属空位影响的局域性特征.
    Lithium-rich manganese-based ternary cathode material for lithium-ion batteries, Li1.208Ni0.333Co0.042Mn0.417O2, has excellent structural stability and electrochemical stability due to its high Ni content. In order to understand the physical properties of this lithium-rich material, its crystal structure, electronic structure and defect properties are calculated by employing the first-principles method based on the density functional theory. The obtained electronic structure shows that Li1.208Ni0.333Co0.042Mn0.417O2 is a magnetic semiconductor with a direct band gap of 0.47 eV. The analysis of the electronic state suggests that the electronic state at the valence band maximum (VBM) is the hybridization of px, py, pz orbitals of oxygen and the dxy, dyz, dxz orbitals of Ni-atom. The electronic state at the conduction band minimum has similar characteristics to those at the VBM, however, part of Ni-${3{\rm{d}}}_{{x}^{2}-{y}^{2}}$ and Mn-${3{\rm{d}}}_{{x}^{2}-{y}^{2}}$, and Mn-${3{\rm{d}}}_{yz}$ also contribute to the electronic hybridizations. The charge density difference calculations indicate that the bonding between O and transition metal atoms are through the mixture of covalent bond with ionic bond. The vacancy formation of a single metal atom is also calculated. The results show that the volumes of the defect systems containing metal vacancies are all reduced in comparison with the volume of perfect lattice. The volume change is the largest for the formation of Mn-vacancy, while the volume is almost unchanged with Co atoms extracted. The vacancy formation energies of the metals are Ef (Mn) > Ef (Co) > Ef (Ni), and the vacancy formation energy of Mn is significantly higher than those of Ni and Co, indicating that the presence of Mn provides a major structural stability for the material. The calculated charge density differences also show that the formation of metal vacancies influences only the charge distribution of the oxygen atoms around the vacancy, showing the local character of the vacancy effect. Since the formation of metal vacancy breaks the bonding between the metal and the surrounding oxygen atoms, the O-2p states near the Fermi surface are significantly increased as shown in the calculated electronic density of states. Such a picture suggests that the electrons on oxygen atoms in vicinity of the metal vacancies become freer.
      通信作者: 朱梓忠, zzhu@xmu.edu.cn
    • 基金项目: 国家重点研发计划 (批准号: 2016YFA0202601, 2016YFB0901502)资助的课题
      Corresponding author: Zhu Zi-Zhong, zzhu@xmu.edu.cn
    • Funds: Project supported by the National Key R&D Program of China (Grant Nos. 2016YFA0202601, 2016YFB0901502)
    [1]

    Whittingham M S 2004 Chem. Rev. 104 4271Google Scholar

    [2]

    Hong J, Gwon H, Jung S K, Ku K, Kang K 2015 J. Electrochem. Soc. 162 A2447Google Scholar

    [3]

    He P, Yu H J, Li D, Zhou H S 2012 J. Mater. Chem. 22 3680Google Scholar

    [4]

    Lu J, Chen Z H, Ma Z F, Pan F, Curtiss L A, Amine K 2016 Nat. Nanotechnol. 11 1031Google Scholar

    [5]

    Deng Y W, Feng C L, Jiaqiang E, Wei K X, Zhang B, Zhang Z Q, Han D D, Zhao X H, Xu W W 2019 Energy 183 869Google Scholar

    [6]

    Cano Z P, Banham D, Ye S Y, Hintennach A, Lu J, Fowler M, Chen Z W 2018 Nat. Energy 3 279Google Scholar

    [7]

    Zeng X Q, Li M, Abd El-Hady D, Alshitari W, Al-Bogami A S, Lu J, Amine K 2019 Adv. Energy Mater. 9 1900161Google Scholar

    [8]

    Goodenough J B, Park K S 2013 J. Am. Chem. Soc. 135 1167Google Scholar

    [9]

    Hou P Y, Zhang H Z, Zi Z Y, Zhang L Q, Xu X J 2017 J. Mater. Chem. A 5 4254Google Scholar

    [10]

    Hou P Y, Yin J M, Ding M, Huang J Z, Xu X J 2017 Small 13 1701802Google Scholar

    [11]

    Kim H, Lee E J, Sun Y K 2014 Mater. Today 17 285Google Scholar

    [12]

    Lee J K, Oh C, Kim N, Hwang J Y, Sun Y K 2016 J. Mater. Chem. A 4 5366Google Scholar

    [13]

    Lee J H, Yoon C S, Hwang J Y, Kim S J, Maglia F, Lamp P, Myung S T, Sun Y K 2016 Energy Environ. Sci. 9 2152Google Scholar

    [14]

    Lu Z H, MacNeil D D, Dahn J R 2001 Electrochem. Solid-State Lett. 4 A191Google Scholar

    [15]

    Etacheri V, Marom R, Elazari R, Salitra G, Aurbach D 2011 Energy Environ. Sci. 4 3243Google Scholar

    [16]

    Lu Z H, MacNeil D D, Dahn J R 2001 Electrochem. Solid-State Lett. 4 A200Google Scholar

    [17]

    Thackeray M M, David W I F, Bruce P G, Goodenough J B 1983 Mater. Res. Bull. 18 461Google Scholar

    [18]

    Manthiram A, Chemelewski K, Lee E S 2014 Energy Environ. Sci. 7 1339Google Scholar

    [19]

    Padhi A K, Nanjundaswamy K S, Goodenough J B 1997 J. Electrochem. Soc. 144 1188Google Scholar

    [20]

    Yuan L X, Wang Z H, Zhang W X, Hu X L, Chen J T, Huang Y H, Goodenough J B 2011 Energy Environ. Sci. 4 269Google Scholar

    [21]

    Nyten A, Kamali S, Haggstrom L, Gustafsson T, Thomas J O 2006 J. Mater. Chem. 16 2266Google Scholar

    [22]

    Ohzuku T, Nagayama M, Tsuji K, Ariyoshi K 2011 J. Mater. Chem. 21 10179Google Scholar

    [23]

    Lee J, Zhang Q, Kim J, Dupre N, Avdeev M, Jeong M, Yoon W S, Gu L, Kang B 2020 Adv. Energy Mater. 10 1902231Google Scholar

    [24]

    Wang J, He X, Paillard E, Laszczynski N, Li J, Passerini S 2016 Adv. Energy Mater. 6 1600906Google Scholar

    [25]

    Zheng H, Hu Z, Liu P, Xu W, Xie Q, He W, Luo Q, Wang L, Gu D, Qu B 2020 Energy Storage Mater. 25 76Google Scholar

    [26]

    Liu W, Oh P, Liu X, Lee M J, Cho W, Chae S, Kim Y, Cho J 2015 Angew Chem. Int. Ed. 54 4440Google Scholar

    [27]

    Gu M, Belharouak I, Genc A, Wang Z G, Wang D P, Amine K, Gao F, Zhou G W, Thevuthasan S, Baer D R, Zhang J G, Browning N D, Liu J, Wang C M 2012 Nano Lett. 12 5186Google Scholar

    [28]

    Kresse G, Furthmuller J 1996 Comput. Mater. Sci. 6 15Google Scholar

    [29]

    Kresse G, Furthmuller J 1996 Phys. Rev. B 54 11169Google Scholar

    [30]

    Blochl P E 1994 Phys. Rev. B 50 17953Google Scholar

    [31]

    Kresse G, Joubert D 1999 Phys. Rev. B 59 1758

    [32]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [33]

    Perdew J P, Ruzsinszky A, Csonka G I, Vydrov O A, Scuseria G E, Constantin L A, Zhou X, Burke K 2008 Phys. Rev. Lett. 100 136406Google Scholar

    [34]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [35]

    Dudarev S L, Botton G A, Savrasov S Y, Humphreys C J, Sutton A P 1998 Phys. Rev. B 57 1505Google Scholar

    [36]

    Jain A, Hautier G, Ong S P, Moore C J, Fischer C C, Persson K A, Ceder G 2011 Phys. Rev. B 84 045115Google Scholar

    [37]

    Chong S K, Liu Y N, Yan W W, Chen Y Z 2016 RSC. Adv. 6 53662Google Scholar

    [38]

    Sun G, Yu F D, Que L F, Deng L, Wang M J, Jiang Y S, Shao G, Wang Z B 2019 Nano Energy 66 104102Google Scholar

    [39]

    Wu S Q, Cai N L, Zhu Z Z, Yang Y 2008 Electrochim. Acta 53 7915Google Scholar

    [40]

    Hu C H, Yang Y, Zhu Z Z 2010 Solid State Commun. 150 669Google Scholar

    [41]

    Turner D E, Zhu Z Z, Chan C T, Mo K M 1997 Phys. Rev. B 55 13842Google Scholar

  • 图 1  Li1.208Ni0.333Co0.042Mn0.417O2的晶体结构, 其中LiO6, NiO6, CoO6和MnO6八面体分别标识为绿色、灰色、蓝色和紫色

    Fig. 1.  Crystal structure of Li1.208Ni0.333Co0.042Mn0.417O2, where LiO6, NiO6, CoO6 and MnO6 octahedra are marked by green, gray, blue and purple, respectively.

    图 2  Li1.208Ni0.333Co0.042Mn0.417O2的能带结构 (a) 自旋向上的能带; (b) 自旋向下的能带. 红色虚线表示费米能级

    Fig. 2.  Band structures of Li1.208Ni0.333Co0.042Mn0.417O2: (a) Spin-up bands; (b) spin-down bands. Fermi level represented by the red dotted line.

    图 3  (a) Li1.208Ni0.333Co0.042Mn0.417O2材料的总态密度和各原子分态密度; (b)—(d) 金属空位形成前后材料中空位周围的6个氧原子的2p电子态密度和, 分别用黑色和红色实线表示

    Fig. 3.  (a) Total and atomic-decomposed partial density of states for Li1.208Ni0.333Co0.042Mn0.417O2; (b)–(d) sum of the density of states of 2p electronic states of the six oxygen atoms around the M-vacancy before and after the formation of M vacancy, respectively. Black and red solid lines represent DOSs before and after the vacancy formation, respectively.

    图 4  材料中三种过渡金属差分电荷密度对比 (a)—(c) 空位形成前; (d)—(f) 空位形成后. 所画的平面是金属附近4个O所在的平面. 图中实线和虚线分别表示电荷聚集$ (\Delta \rho > 0) $与电荷移出$ (\Delta \rho < 0) $

    Fig. 4.  Contour plots of the charge density differences: (a)–(c) Before the M-vacancy formation; (d)–(f) after the M-vacancy formation. The solid lines and the dot lines represent the charge accumulation (∆ρ > 0) and the charge depletion (∆ρ < 0), respectively.

    表 1  空位形成前后Li1.208Ni0.333Co0.042Mn0.417O2的晶格常数, 原胞体积和空位形成能的计算值

    Table 1.  Computational lattice constants, unit cell volume and vacancy formation energies of Li1.208Ni0.333Co0.042Mn0.417O2 before and after the defect formations.

    材料体系晶格常数体积/Å3体积变化量(VvacV0)/V0空位形成能$ /\mathrm{e}\mathrm{V} $
    abc
    无空位时5.78611.55214.302822.285
    含Ni空位5.77811.49214.268815.131–0.87%5.31
    含Co空位5.79611.43414.343821.290–0.12%6.52
    含Mn空位5.73211.52314.231811.067–1.36%8.27
    下载: 导出CSV

    表 2  各原子的电子轨道对价带顶和导带底上电子态的贡献

    Table 2.  Contribution of electronic orbital of various types of atoms for the electronic states at VBM and CBM.

    $ {\mathrm{p}}_{y} $$ {\mathrm{p}}_{z} $$ {\mathrm{p}}_{x} $$ {\mathrm{d}}_{xy} $$ {\mathrm{d}}_{yz} $$ {\mathrm{d}}_{{z}^{2}} $$ {\mathrm{d}}_{xz} $$ {\mathrm{d}}_{{x}^{2}-{y}^{2}} $
    VBM-O0.1380.2060.26800000
    CBM-O0.1340.1640.09700000
    VBM-Ni0000.0480.0290.0040.0950.01
    CBM-Ni0000.0280.14100.0680.102
    CBM-Mn0000.0030.0150.0070.0030.034
    下载: 导出CSV
  • [1]

    Whittingham M S 2004 Chem. Rev. 104 4271Google Scholar

    [2]

    Hong J, Gwon H, Jung S K, Ku K, Kang K 2015 J. Electrochem. Soc. 162 A2447Google Scholar

    [3]

    He P, Yu H J, Li D, Zhou H S 2012 J. Mater. Chem. 22 3680Google Scholar

    [4]

    Lu J, Chen Z H, Ma Z F, Pan F, Curtiss L A, Amine K 2016 Nat. Nanotechnol. 11 1031Google Scholar

    [5]

    Deng Y W, Feng C L, Jiaqiang E, Wei K X, Zhang B, Zhang Z Q, Han D D, Zhao X H, Xu W W 2019 Energy 183 869Google Scholar

    [6]

    Cano Z P, Banham D, Ye S Y, Hintennach A, Lu J, Fowler M, Chen Z W 2018 Nat. Energy 3 279Google Scholar

    [7]

    Zeng X Q, Li M, Abd El-Hady D, Alshitari W, Al-Bogami A S, Lu J, Amine K 2019 Adv. Energy Mater. 9 1900161Google Scholar

    [8]

    Goodenough J B, Park K S 2013 J. Am. Chem. Soc. 135 1167Google Scholar

    [9]

    Hou P Y, Zhang H Z, Zi Z Y, Zhang L Q, Xu X J 2017 J. Mater. Chem. A 5 4254Google Scholar

    [10]

    Hou P Y, Yin J M, Ding M, Huang J Z, Xu X J 2017 Small 13 1701802Google Scholar

    [11]

    Kim H, Lee E J, Sun Y K 2014 Mater. Today 17 285Google Scholar

    [12]

    Lee J K, Oh C, Kim N, Hwang J Y, Sun Y K 2016 J. Mater. Chem. A 4 5366Google Scholar

    [13]

    Lee J H, Yoon C S, Hwang J Y, Kim S J, Maglia F, Lamp P, Myung S T, Sun Y K 2016 Energy Environ. Sci. 9 2152Google Scholar

    [14]

    Lu Z H, MacNeil D D, Dahn J R 2001 Electrochem. Solid-State Lett. 4 A191Google Scholar

    [15]

    Etacheri V, Marom R, Elazari R, Salitra G, Aurbach D 2011 Energy Environ. Sci. 4 3243Google Scholar

    [16]

    Lu Z H, MacNeil D D, Dahn J R 2001 Electrochem. Solid-State Lett. 4 A200Google Scholar

    [17]

    Thackeray M M, David W I F, Bruce P G, Goodenough J B 1983 Mater. Res. Bull. 18 461Google Scholar

    [18]

    Manthiram A, Chemelewski K, Lee E S 2014 Energy Environ. Sci. 7 1339Google Scholar

    [19]

    Padhi A K, Nanjundaswamy K S, Goodenough J B 1997 J. Electrochem. Soc. 144 1188Google Scholar

    [20]

    Yuan L X, Wang Z H, Zhang W X, Hu X L, Chen J T, Huang Y H, Goodenough J B 2011 Energy Environ. Sci. 4 269Google Scholar

    [21]

    Nyten A, Kamali S, Haggstrom L, Gustafsson T, Thomas J O 2006 J. Mater. Chem. 16 2266Google Scholar

    [22]

    Ohzuku T, Nagayama M, Tsuji K, Ariyoshi K 2011 J. Mater. Chem. 21 10179Google Scholar

    [23]

    Lee J, Zhang Q, Kim J, Dupre N, Avdeev M, Jeong M, Yoon W S, Gu L, Kang B 2020 Adv. Energy Mater. 10 1902231Google Scholar

    [24]

    Wang J, He X, Paillard E, Laszczynski N, Li J, Passerini S 2016 Adv. Energy Mater. 6 1600906Google Scholar

    [25]

    Zheng H, Hu Z, Liu P, Xu W, Xie Q, He W, Luo Q, Wang L, Gu D, Qu B 2020 Energy Storage Mater. 25 76Google Scholar

    [26]

    Liu W, Oh P, Liu X, Lee M J, Cho W, Chae S, Kim Y, Cho J 2015 Angew Chem. Int. Ed. 54 4440Google Scholar

    [27]

    Gu M, Belharouak I, Genc A, Wang Z G, Wang D P, Amine K, Gao F, Zhou G W, Thevuthasan S, Baer D R, Zhang J G, Browning N D, Liu J, Wang C M 2012 Nano Lett. 12 5186Google Scholar

    [28]

    Kresse G, Furthmuller J 1996 Comput. Mater. Sci. 6 15Google Scholar

    [29]

    Kresse G, Furthmuller J 1996 Phys. Rev. B 54 11169Google Scholar

    [30]

    Blochl P E 1994 Phys. Rev. B 50 17953Google Scholar

    [31]

    Kresse G, Joubert D 1999 Phys. Rev. B 59 1758

    [32]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [33]

    Perdew J P, Ruzsinszky A, Csonka G I, Vydrov O A, Scuseria G E, Constantin L A, Zhou X, Burke K 2008 Phys. Rev. Lett. 100 136406Google Scholar

    [34]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [35]

    Dudarev S L, Botton G A, Savrasov S Y, Humphreys C J, Sutton A P 1998 Phys. Rev. B 57 1505Google Scholar

    [36]

    Jain A, Hautier G, Ong S P, Moore C J, Fischer C C, Persson K A, Ceder G 2011 Phys. Rev. B 84 045115Google Scholar

    [37]

    Chong S K, Liu Y N, Yan W W, Chen Y Z 2016 RSC. Adv. 6 53662Google Scholar

    [38]

    Sun G, Yu F D, Que L F, Deng L, Wang M J, Jiang Y S, Shao G, Wang Z B 2019 Nano Energy 66 104102Google Scholar

    [39]

    Wu S Q, Cai N L, Zhu Z Z, Yang Y 2008 Electrochim. Acta 53 7915Google Scholar

    [40]

    Hu C H, Yang Y, Zhu Z Z 2010 Solid State Commun. 150 669Google Scholar

    [41]

    Turner D E, Zhu Z Z, Chan C T, Mo K M 1997 Phys. Rev. B 55 13842Google Scholar

  • [1] 吴洪芬, 冯盼君, 张烁, 刘大鹏, 高淼, 闫循旺. 铁原子吸附联苯烯单层电子结构的第一性原理. 物理学报, 2022, 71(3): 036801. doi: 10.7498/aps.71.20211631
    [2] 史晓红, 陈京金, 曹昕睿, 吴顺情, 朱梓忠. 富锂锰基三元材料Li1.167Ni0.167Co0.167Mn0.5O2中的氧空位形成:第一性原理计算. 物理学报, 2022, 0(0): 0-0. doi: 10.7498/aps.71.20220274
    [3] 吴洪芬, 冯盼君, 张烁, 刘大鹏, 高淼, 闫循旺. 铁原子吸附联苯烯单层电子结构的第一性原理研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211631
    [4] 黄炳铨, 周铁戈, 吴道雄, 张召富, 李百奎. 空位及氮掺杂二维ZnO单层材料性质:第一性原理计算与分子轨道分析. 物理学报, 2019, 68(24): 246301. doi: 10.7498/aps.68.20191258
    [5] 闫小童, 侯育花, 郑寿红, 黄有林, 陶小马. Ga, Ge, As掺杂对锂离子电池正极材料Li2CoSiO4的电化学特性和电子结构影响的第一性原理研究. 物理学报, 2019, 68(18): 187101. doi: 10.7498/aps.68.20190503
    [6] 刘琪, 管鹏飞. La65X35(X=Ni,Al)非晶合金原子结构的第一性原理研究. 物理学报, 2018, 67(17): 178101. doi: 10.7498/aps.67.20180992
    [7] 彭琼, 何朝宇, 李金, 钟建新. MoSi2薄膜电子性质的第一性原理研究. 物理学报, 2015, 64(4): 047102. doi: 10.7498/aps.64.047102
    [8] 王平, 郭立新, 杨银堂, 张志勇. 铝氮共掺杂氧化锌纳米管电子结构的第一性原理研究. 物理学报, 2013, 62(5): 056105. doi: 10.7498/aps.62.056105
    [9] 焦照勇, 郭永亮, 牛毅君, 张现周. 缺陷黄铜矿结构Xga2S4 (X=Zn, Cd, Hg)晶体电子结构和光学性质的第一性原理研究. 物理学报, 2013, 62(7): 073101. doi: 10.7498/aps.62.073101
    [10] 邓娇娇, 刘波, 顾牡, 刘小林, 黄世明, 倪晨. 伽马CuX(X=Cl,Br,I)的电子结构和光学性质的第一性原理计算. 物理学报, 2012, 61(3): 036105. doi: 10.7498/aps.61.036105
    [11] 焦照勇, 杨继飞, 张现周, 马淑红, 郭永亮. 闪锌矿GaN弹性性质、电子结构和光学性质外压力效应的理论研究. 物理学报, 2011, 60(11): 117103. doi: 10.7498/aps.60.117103
    [12] 吴红丽, 赵新青, 宫声凯. Nb掺杂影响NiTi金属间化合物电子结构的第一性原理计算. 物理学报, 2010, 59(1): 515-520. doi: 10.7498/aps.59.515
    [13] 杨天兴, 成强, 许红斌, 王渊旭. 几种三元过渡金属碳化物弹性及电子结构的第一性原理研究. 物理学报, 2010, 59(7): 4919-4924. doi: 10.7498/aps.59.4919
    [14] 张学军, 高攀, 柳清菊. 氮铁共掺锐钛矿相TiO2电子结构和光学性质的第一性原理研究. 物理学报, 2010, 59(7): 4930-4938. doi: 10.7498/aps.59.4930
    [15] 李沛娟, 周薇薇, 唐元昊, 张华, 施思齐. CeO2的电子结构,光学和晶格动力学性质:第一性原理研究. 物理学报, 2010, 59(5): 3426-3431. doi: 10.7498/aps.59.3426
    [16] 汪志刚, 张杨, 文玉华, 朱梓忠. ZnO原子链的结构稳定性和电子性质的第一性原理研究. 物理学报, 2010, 59(3): 2051-2056. doi: 10.7498/aps.59.2051
    [17] 胡方, 明星, 范厚刚, 陈岗, 王春忠, 魏英进, 黄祖飞. 梯形化合物NaV2O4F电子结构的第一性原理研究. 物理学报, 2009, 58(2): 1173-1178. doi: 10.7498/aps.58.1173
    [18] 宋庆功, 王延峰, 宋庆龙, 康建海, 褚 勇. 插层化合物Ag1/4TiSe2电子结构的第一性原理研究. 物理学报, 2008, 57(12): 7827-7832. doi: 10.7498/aps.57.7827
    [19] 黄 丹, 邵元智, 陈弟虎, 郭 进, 黎光旭. 纤锌矿结构Zn1-xMgxO电子结构及吸收光谱的第一性原理研究. 物理学报, 2008, 57(2): 1078-1083. doi: 10.7498/aps.57.1078
    [20] 吴红丽, 赵新青, 宫声凯. Nb掺杂对TiO2/NiTi界面电子结构影响的第一性原理计算. 物理学报, 2008, 57(12): 7794-7799. doi: 10.7498/aps.57.7794
计量
  • 文章访问数:  991
  • PDF下载量:  25
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-02
  • 修回日期:  2021-06-03
  • 上网日期:  2021-08-15
  • 刊出日期:  2021-10-20

富锂锰基三元材料Li1.208Ni0.333Co0.042Mn0.417O2的电子结构和缺陷性质

  • 厦门大学物理科学与技术学院, 厦门 361005
  • 通信作者: 朱梓忠, zzhu@xmu.edu.cn
    基金项目: 国家重点研发计划 (批准号: 2016YFA0202601, 2016YFB0901502)资助的课题

摘要: 采用基于密度泛函理论的第一性原理方法, 计算了$ R\bar{3}m $相的富锂锰基三元正极材料Li1.208Ni0.333Co0.042Mn0.417O2的晶体结构、电子结构以及缺陷性质. 结果表明, Li1.208Ni0.333Co0.042Mn0.417O2是一种具有直接带隙的磁性半导体材料, 其价带顶是O的px, py, pz轨道与Ni的dxy, dyz, dxz原子轨道之间的杂化, 导带底除了有价带顶的特征外, 还有部分的Ni-$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $和Mn-$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $, Mn-$ {3\mathrm{d}}_{yz} $ 轨道参与杂化. 差分电荷密度图显示, 金属原子与O原子之间的成键方式是共价键和离子键的混合. 本文还计算了脱去单个过渡金属原子的空位形成能. 脱去Mn原子时的缺陷体系的体积变化最大, 而脱去Co原子的缺陷体系的体积则几乎不变. 得到的金属空位形成能的大小依次为Ef (Mn-空位) > Ef (Co-空位) > Ef (Ni-空位). 差分电荷密度图还显示, 空位的产生仅强烈影响了空位附近的部分O原子, 体现了金属空位影响的局域性特征.

English Abstract

    • 锂离子电池自商业化以来, 因具有相对较高的比能量、良好的倍率性能和有效的循环稳定性而得以迅速发展, 并广泛应用于便携式电子产品、电动汽车和航空航天等领域[1-5]. 随着大规模储能应用和技术的迅速发展, 当前的锂离子电池已经不能满足日益增长的对能量密度的要求, 研发更高能量密度的下一代锂离子电池至关重要[6,7]. 电极材料在电池的电化学行为中起着关键的作用, 其中第一个商业化的 LiCoO2正极和碳负极仍然主导着锂离子电池领域[8-10]. 但是, 与已开发的可逆容量超过1000 mA/g且循环稳定性和低成本的硅/碳复合负极材料相比[11,12], 层状LiCoO2正极实际上只能提供约150 mA/g (其理论比容量可以达到275 mA/g)[9,13]. 此外, 地球上钴源的稀缺性导致LiCoO2材料的高成本, 也限制了它在电动汽车中的应用.

      在过去的几十年中, 开发高容量、低成本且具有持久的循环寿命以及优异的充放速率和安全性能的正极材料一直是锂离子电池研究的重要方向. 在众多的正极候选材料中, 层状过渡金属氧化物, 包括Li(NixCoyMn1–xy)O2 (即NCM三元材料) 和富锂的Li1+z(Mn1–xMx)1–zO2 (M = Ni, Co, Ru, Sn等)[14-16], 尖晶石型的LiM2O4 (M = Ni, Mn)[17,18] 和橄榄石型的LiMXO4 (M = Fe, Mn, Co; X = P, Si)[19-21]正极材料得到了特别的重视. 富锂三元层状氧化物正极材料, 由于其高于300 mA/g的比容量, 大于1000 W·h/kg的能量密度以及丰富锰源带来的低成本, 被认为是下一代锂离子电池中最具有应用前景的正极材料之一[22-24]. 本文的研究对象就是一种富锂锰基的三元层状氧化物Li1.208Ni0.333Co0.042Mn0.417O2, 该材料有着高的放电电压、优异的结构和循环稳定性、高的电压/容量保持率[25]. 当它用于包覆其他的富锂三元材料时, 通过表面富Ni的设计能够增加整个体系的结构稳定性, 抑制其容量减小以及电压衰减[25]. 钴元素不仅是活性材料, 而且可以有效抑制阳离子的混排 (Li/Ni) 并稳定材料的结构, 有助于改善材料的深度放电特性并提高倍率性能; 锰元素对材料有支撑作用, 也可以抑制在大量脱锂时阳离子的混排, 但是锰含量的上升会部分降低材料的容量; 此外, 镍元素含量的上升能够提高材料的工作电压, 但由于Ni在Li层中较小的扩散势垒更容易在材料表面选择性偏析, 降低材料的结构稳定性[3,26,27]. 通过第一性原理的方法, 本文对Li1.208Ni0.333Co0.042Mn0.417O2材料的结构和电子结构性质进行了仔细的计算, 表明了该组分的层状三元材料是一个具有直接带隙的磁性半导体. 通过计算各个过渡金属原子的空位形成能, 探讨了Ni, Co和Mn元素在材料中的键合强度, 结果表明Ni元素最容易从材料中脱出, 而Mn元素最难从材料中脱出. 这些研究能够帮助理解这个具有复杂组分的三元材料的基本性质.

    • 本文采用的是基于密度泛函理论 (DFT) 的第一性原理计算方法, 通过从头计算的VASP程序包[28,29]来实现. 该程序包基于平面波展开和投影缀加波 (projector augmented-wave, PAW)[30,31]方法. 计算中, 交换关联泛函使用的是广义梯度近似 (GGA)[32]下的PBE (Perdew-Burke-Ernzerhof)[33]表述, 波函数则按平面波基展开, 展开时的平面波截断动能为550 eV. 布里渊区 (BZ) 上的积分替换为倒格空间特殊k点集上的离散求和, 采用了Monkhorst-pack 的5 × 2 × 2k 点网格[34]. 采用的超晶胞原胞的晶格常数和原胞内的原子坐标都进行了充分的弛豫, 直到原胞内所有原子上的Hellmann-Feynman力均小于0.2 eV/nm为止. 此外, 计算还采用了GGA + U [35]的方法来处理具有局域性特点的Ni, Co和Mn的3d电子的强关联性质, 其中有效相互作用参数Ueff值参考了Jain等[36]报道的数值, 即对Ni, Co 和Mn原子分别使用了Ueff = 6.0, 3.4和3.9 eV. 由于过渡金属d电子的局域性强关联特点, 使用GGA + U方法是必要的. 整个计算是在自旋极化的框架下进行的, 因为材料的电子结构会受到Ni, Co, Mn原子磁性的明显影响.

    • Li1.208Ni0.333Co0.042Mn0.417O2具有与LiCoO2类似的α-NaFeO2型层状结构, 属于六方晶系, 空间群为$ R\bar{3}m $. 每个原胞含有24个单位分子式 (单位分子式可以从LiCoO2来看, 当体系富锂加上三元时即为Li1.208Ni0.333Co0.042Mn0.417O2), 共有96个原子/原胞, 其中含有29个Li原子, 8个Ni原子, 1个Co原子, 10个Mn原子以及48个O原子, 晶体结构如图1所示. 可以看出, 所有Li离子和过渡金属 (TM = Ni, Co, Mn) 离子都处在密堆积氧晶格八面体的中心位置而形成LiO6八面体和TMO6八面体. 作为富锂的层状结构, 如果将Li1.208Ni0.333Co0.042Mn0.417O2重新写为LiTMO2 (其中TM = Li0.208Ni0.333Co0.042Mn0.417) 的结构式, 则晶格可看作是紧密堆积的Li-O-TM-O层的交替排列, 其中的TM层由金属元素和富余的Li有序镶嵌形成. 整个晶体可以看作是由共享边的LiO6八面体层和TMO6八面体层在c方向上的交替堆垛而成, 这种结构非常适合于锂离子的嵌入与脱出.

      图  1  Li1.208Ni0.333Co0.042Mn0.417O2的晶体结构, 其中LiO6, NiO6, CoO6和MnO6八面体分别标识为绿色、灰色、蓝色和紫色

      Figure 1.  Crystal structure of Li1.208Ni0.333Co0.042Mn0.417O2, where LiO6, NiO6, CoO6 and MnO6 octahedra are marked by green, gray, blue and purple, respectively.

      表1中给出了Li1.208Ni0.333Co0.042Mn0.417O2以及带有过渡金属单空位缺陷材料的晶格常数、原胞体积、体积变化率和空位形成能. 空位的计算时, 直接在计算完整晶体 (空位形成前) 的原胞中拿去一个相应的金属原子而形成空位. 所以, 空位计算的超原胞是95原子/原胞. 从表1中的结果可以看出, 脱去一个Mn原子, 体积变小最多; 脱去一个Ni原子体积变小次之; 而脱去单个Co原子后, 体积几乎没有变化. 再来看一下单个金属M的空位形成能, 其公式如下:

      材料体系晶格常数体积/Å3体积变化量(VvacV0)/V0空位形成能$ /\mathrm{e}\mathrm{V} $
      abc
      无空位时5.78611.55214.302822.285
      含Ni空位5.77811.49214.268815.131–0.87%5.31
      含Co空位5.79611.43414.343821.290–0.12%6.52
      含Mn空位5.73211.52314.231811.067–1.36%8.27

      表 1  空位形成前后Li1.208Ni0.333Co0.042Mn0.417O2的晶格常数, 原胞体积和空位形成能的计算值

      Table 1.  Computational lattice constants, unit cell volume and vacancy formation energies of Li1.208Ni0.333Co0.042Mn0.417O2 before and after the defect formations.

      $ \Delta E=E\left(\mathrm{d}\mathrm{e}\mathrm{f}/M\right)-E\left(\mathrm{p}\mathrm{e}\mathrm{f}\right)+E\left(M\right), $

      其中$ E\left(\mathrm{p}\mathrm{e}\mathrm{f}\right) $ 表示完整晶体的总结合能, $ E\left(\mathrm{d}\mathrm{e}\mathrm{f}/M\right) $ 表示脱去单个金属原子后缺陷晶体的总结合能, $ E\left(M\right) $M的体相中单个过渡金属原子的结合能 (M = Ni, Co, Mn). 计算结果表明, Mn空位的形成能最大, 意味着Mn空位最难形成, 这与体系是锰基材料相一致. Mn对材料的结构稳定性贡献最大. 此外, Ni空位的形成能最小, 这意味着Ni元素最容易脱离材料, 也是三个金属元素中最活跃的. Ni2+ 离子的半径约为0.69 Å, 这与Li+ 的半径 0.73 Å相近. 实验结果表明, Ni2+容易跳出TM层而与Li层的Li+ 发生交换, 形成阳离子混合[27,37], 这与计算的Ni空位形成能值较低相一致. 而且前面还提到, Mn含量的上升有利于提高结构稳定性, 这也与理论计算的Mn缺陷形成能较大相一致.

      Li1.208Ni0.333Co0.042Mn0.417O2材料的电子结构计算是在自旋极化的密度泛函理论框架下进行的, 这是因为材料中含有的Ni, Co和Mn是3d过渡金属元素, 这些元素使所构成的材料有强烈的自旋极化倾向. 图2给出了自旋向上和自旋向下的Li1.208Ni0.333Co0.042Mn0.417O2材料的能带结构图. 从图2可以看出, 自旋向上和自旋向下的能带表现出非常不同的电子行为, 表明该材料是一个磁性半导体. 能带结构的价带顶 (VBM) 和导带底 (CBM) 都位于Γ点位置, 所以Li1.208Ni0.333Co0.042Mn0.417O2是一个直接带隙的磁性半导体, 且在Γ点位置的带隙为0.47 eV.

      图  2  Li1.208Ni0.333Co0.042Mn0.417O2的能带结构 (a) 自旋向上的能带; (b) 自旋向下的能带. 红色虚线表示费米能级

      Figure 2.  Band structures of Li1.208Ni0.333Co0.042Mn0.417O2: (a) Spin-up bands; (b) spin-down bands. Fermi level represented by the red dotted line.

      这里重点分析价带顶和导带底的电子波函数的电子轨道贡献 (即来自各原子各电子轨道的贡献情况), 具体表达式为

      $ {\varPsi }_{nk}\left({\boldsymbol{r}}\right)=\sum _{\mu =1}^{N}\sum _{lm}{C}_{lm}(n,{\boldsymbol{k}},{{\boldsymbol{R}}}_{\mu }){\varphi }_{lm}^{\mathrm{a}\mathrm{t}\mathrm{o}\mathrm{m}\mathrm{i}\mathrm{c}}({\boldsymbol{r}}-{{\boldsymbol{R}}}_{\mu }), $

      式中, 左边是 (n, k) 态的电子波函数, 右边是该波函数对各原子的各电子轨道 $ {\varphi }_{lm}^{\mathrm{a}\mathrm{t}\mathrm{o}\mathrm{m}\mathrm{i}\mathrm{c}}({\boldsymbol{r}}-{{\boldsymbol{R}}}_{\mu }) $的展开, 展开系数$ {C}_{lm}(n, {\boldsymbol{k}}, {{\boldsymbol{R}}}_{\mu }) $ 表达了某个电子轨道对某波函数的贡献程度. 表2列出了价带顶和导带底的电子轨道的贡献. 由表中结果可见, 价带顶的电子态主要是O的2px, 2py, 2pz轨道与Ni的3dxy, 3dyz, 3dxz原子轨道之间的杂化, Co和Mn原子的电子轨道基本上没有参与; 导带底电子态的原子轨道贡献除了有价带顶的特征外, 还有部分的Ni-$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $和Mn-$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $, Mn-$ {3\mathrm{d}}_{yz} $ 轨道参与杂化. 由过渡金属中贡献最多的是Ni的3d轨道电子可见, Ni在Li1.208Ni0.333Co0.042Mn0.417O2材料的过渡金属中电子性质最为活泼. 这与各过渡金属原子在三元材料中不同价态的能级排列特点[38] 是一致的 (即 Ni2+, Co3+, Mn4+费米面下的能级顺序分别是越来越远离费米面的).

      $ {\mathrm{p}}_{y} $$ {\mathrm{p}}_{z} $$ {\mathrm{p}}_{x} $$ {\mathrm{d}}_{xy} $$ {\mathrm{d}}_{yz} $$ {\mathrm{d}}_{{z}^{2}} $$ {\mathrm{d}}_{xz} $$ {\mathrm{d}}_{{x}^{2}-{y}^{2}} $
      VBM-O0.1380.2060.26800000
      CBM-O0.1340.1640.09700000
      VBM-Ni0000.0480.0290.0040.0950.01
      CBM-Ni0000.0280.14100.0680.102
      CBM-Mn0000.0030.0150.0070.0030.034

      表 2  各原子的电子轨道对价带顶和导带底上电子态的贡献

      Table 2.  Contribution of electronic orbital of various types of atoms for the electronic states at VBM and CBM.

      计算得到的Li1.208Ni0.333Co0.042Mn0.417O2材料的总电子态密度 (TDOS) 和各原子分态密度 (PDOS) 如图3(a)所示. 从图3(a)可以清楚地看到, 该材料是一个半导体, 而且是一个自旋极化的磁性半导体. 在图3(b)图3(d)中, 分别给出了Li1.208Ni0.333Co0.042Mn0.417O2中空位周围的6个O原子的2p电子分态密度和在过渡金属空位 (即Ni, Co和Mn空位) 形成之前和形成之后的对比. 可以看到, 当金属空位形成后, 费米面附近的O-2p占据态 (DOS峰) 明显增加. 这是因为, 金属空位的形成打断了金属与周围O原子的成键, 导致周围这些氧原子上的电子更加自由, 能量上就更加接近于费米面. 带有Co, Ni 空位的体系的态密度特征类似, 即费米面附近有了很高的态密度值 (较自由的电子较多). 而Mn空位体系的DOS稍有不同, 其高态密度的峰离费米面稍远 (电子被束缚的较紧).

      图  3  (a) Li1.208Ni0.333Co0.042Mn0.417O2材料的总态密度和各原子分态密度; (b)—(d) 金属空位形成前后材料中空位周围的6个氧原子的2p电子态密度和, 分别用黑色和红色实线表示

      Figure 3.  (a) Total and atomic-decomposed partial density of states for Li1.208Ni0.333Co0.042Mn0.417O2; (b)–(d) sum of the density of states of 2p electronic states of the six oxygen atoms around the M-vacancy before and after the formation of M vacancy, respectively. Black and red solid lines represent DOSs before and after the vacancy formation, respectively.

      根据差分电荷密度图可以很好地了解Li1.208Ni0.333Co0.042Mn0.417O2材料中原子之间的成键特征以及缺陷对原子间键合性质的影响. 图4(a)图4(f)分别给出了Li1.208Ni0.333Co0.042Mn0.417O2材料中在三种过渡金属空位形成之前和形成之后的差分电荷密度对比, 所画的平面是金属附近4个O所在的平面. 差分电荷密度定义为总自洽电荷密度与独立原子电荷密度的叠加之差, 即 $\Delta \rho ({\boldsymbol{r}}) = $$ \rho \left({\boldsymbol{r}}\right)-\displaystyle\sum\nolimits _{\mu }{\rho }_{\mathrm{a}\mathrm{t}\mathrm{o}\mathrm{m}}\left({\boldsymbol{r}}-{{\boldsymbol{R}}}_{\mu }\right)$, 其中$ \rho \left({\boldsymbol{r}}\right) $表示体系的自洽总电荷密度, $\displaystyle\sum\nolimits _{\mu }{\rho }_{\mathrm{a}\mathrm{t}\mathrm{o}\mathrm{m}}\left({\boldsymbol{r}}-{{\boldsymbol{R}}}_{\mu }\right)$表示各原子的电荷密度的叠加, $ {{\boldsymbol{R}}}_{\mu } $ 是原子位置[39,40]. 为了有效地比较各个面内电荷分布的差异, 图4中各图的差分电荷密度等高线的取值范围均相同 (即: –1.0—0.7 e3) 且等高线的间隔值 (0.05 e3) 也相同. 图4中红色实线表示该区域的电荷聚集 (相对于独立原子的电荷密度之和而言), 而蓝色虚线的区域表示着电荷的缺失. 从图4可见, 金属元素M (M = Ni, Co, Mn) 周围有明显的电荷移出, 而电荷较明显地聚集在O原子周围. 还可以清楚地看到, M-O原子之间的电荷分布还表现出显著的方向性 (共价性), 所以M (M = Ni, Co, Mn) 与O之间的键合特征是很强的共价键与离子键的混合. 当形成金属元素M的空位后, 从图4可以看到, 金属空位所产生的影响在这个材料中是非常局域的 (即仅影响空位附近的电荷密度分布), 这种情况与简单金属中金属空位的影响[41]是非常不同的 (金属铝中铝原子空位的影响是相当长程的). 对比不同过渡金属空位的情况 (图4) 可以看到, Co空位的产生对其周围O原子电荷密度的影响最为明显, 而Ni空位的出现对其周边的O原子的影响最小. Mn, Co空位的出现都对其周围不等价氧原子有着明显的不同影响. 三种金属空位对邻近的过渡金属虽然有一定的影响, 但都比较小 (距离较远).

      图  4  材料中三种过渡金属差分电荷密度对比 (a)—(c) 空位形成前; (d)—(f) 空位形成后. 所画的平面是金属附近4个O所在的平面. 图中实线和虚线分别表示电荷聚集$ (\Delta \rho > 0) $与电荷移出$ (\Delta \rho < 0) $

      Figure 4.  Contour plots of the charge density differences: (a)–(c) Before the M-vacancy formation; (d)–(f) after the M-vacancy formation. The solid lines and the dot lines represent the charge accumulation (∆ρ > 0) and the charge depletion (∆ρ < 0), respectively.

    • 本文基于密度泛函理论的第一性原理方法, 计算了$ R\bar{3}m $相的富锂锰基三元材料Li1.208Ni0.333Co0.042Mn0.417O2及其含过渡金属原子缺陷体系的晶体结构、电子结构和缺陷性质. 计算的能带结构表明, 该材料为直接带隙的磁性半导体, 在Γ点的直接带隙为0.47 eV. 价带顶上的电子轨道贡献主要来自O-2p和Ni-3d的轨道电子, 而Co和Mn原子的3d电子基本上没有贡献; 导带底上的电子轨道贡献也主要来自O-2p和Ni-3d的轨道电子, 但此时 Mn的 $ {3\mathrm{d}}_{yz} $$ {3\mathrm{d}}_{{x}^{2}-{y}^{2}} $ 的轨道电子也有部分参与. 计算结果还表明, 锰空位形成时的体积变化率最大且空位的形成能也最大, 意味着Mn空位最难形成, 而且对材料的结构稳定性贡献最大. Ni空位的形成能最小, 这意味着Ni元素最容易脱离材料的. 差分电荷密度的计算显示, 金属原子与O原子之间的成键方式主要是共价键和离子键的混合. 金属空位的产生仅强烈影响空位附近O原子周围的电荷密度分布, 这体现了金属空位影响的局域性特征. 本研究可以为理解富锂锰基三元正极材料的组分、结构、电子结构和缺陷性质提供参考.

参考文献 (41)

目录

    /

    返回文章
    返回