搜索

x

留言板

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

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

In掺杂h-LuFeO3光吸收及极化性能的第一性原理计算

张小娅 宋佳讯 王鑫豪 王金斌 钟向丽

引用本文:
Citation:

In掺杂h-LuFeO3光吸收及极化性能的第一性原理计算

张小娅, 宋佳讯, 王鑫豪, 王金斌, 钟向丽

First principles calculation of optical absorption and polarization properties of In doped h-LuFeO3

Zhang Xiao-Ya, Song Jia-Xun, Wang Xin-Hao, Wang Jin-Bin, Zhong Xiang-Li
PDF
HTML
导出引用
  • h-LuFeO3是一种窄带隙铁电半导体材料, 已被证明在铁电光伏领域有较好的应用前景. 然而, 较低的极化强度使光生电子-空穴对复合率大, 限制了h-LuFeO3基铁电光伏电池效率的提高. 为改善h-LuFeO3的极化强度, 提高光吸收性质, 本文利用第一性原理计算方法研究了In原子在h-LuFeO3不同位置的掺杂形成能, 得到最稳定的掺杂位置, 比较了h-Lu1–xInxFeO3 (x = 0, 0.167, 0.333, 0.667)的带隙、光吸收性能及极化强度等性质. 计算结果表明, 随着In掺杂比例的增加, h-LuFeO3的晶格常数c/a比不断增大, 铁电极化强度得到提高. 当In∶Lu = 2∶1时, 材料杂质能级出现, Fe-O轨道杂化得到增强, 提高了h-LuFeO3的光吸收性能. 此工作证明了In掺杂是改善h-LuFeO3极化强度和光吸收系数的有效方法, 对铁电光伏性能的提高提供一种新途径.
    The h-LuFeO3 is a kind of narrow band gap hexagonal ferrite material, with a good application prospect in the field of ferroelectric photovoltaic. However, the low polarization intensity of h-LuFeO3 makes the recombination rate of photogenerated electrons and holes large, which is not conducive to the improvement of the efficiency of h-LuFeO3-based ferroelectric photovoltaic cells. In order to improve the ferroelectricity and optical absorption properties of h-LuFeO3, the first principles method is used to calculate the doping formation energy values of In atom at different positions of h-LuFeO3, and the most stable doping position is determined. The comparisons of band gap, optical absorption performance and polarization intensity among h-Lu1-xInxFeO3 (x = 0, 0.167, 0.333, 0.667) are made. With the increase of In doping, the cells of h-Lu1–xInxFeO3 stretch along the c-axis. The ratio of the lattice constant c/a increases from 1.94 at x = 0 to 2.04 at x = 0.667 when all the positions of In replace P1 position. Using the qualitative calculation of Berne effective charge, the results show that the ferroelectric polarization intensity of h-LuFeO3, h-Lu0.833In0.167FeO3, h-Lu0.667In0.333FeO3 and h-Lu0.333In0.667FeO3 along the c-axis are 3.93, 5.91, 7.92, and 11.02 μC·cm–2, respectively. Therefore, with the increase of the number of In atoms replacing Lu atoms, the lattice constant c/a ratio of h-Lu1–xInxFeO3 increases, which can improve the ferroelectric polarization strength of the material. By analyzing the density of states of h-LuFeO3 and h-Lu0.333In0.667FeO3, we can see that In doping enhances the Fe-O orbital hybridization in h-Lu0.333In0.667FeO3, and makes the optical absorption coefficient of h-Lu0.333In0.667FeO3 in the solar light range larger. In summary, In doped h-LuFeO3 is an effective method to improve its polarization intensity and optical absorption coefficient, which is of great significance for improving the performance of ferroelectric photovoltaic.
      通信作者: 钟向丽, xlzhong@xtu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 51872251, 11875229)和电子元器件可靠性物理及其应用技术重点实验室开放基金(批准号: ZHD201803)资助的课题
      Corresponding author: Zhong Xiang-Li, xlzhong@xtu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 51872251, 11875229) and the Opening Project of Science and Technology on Reliability Physics and Application Technology of Electronic Component Laboratory, China (Grant No. ZHD201803)
    [1]

    Ji Y, Gao T, Wang Z L, Yang Y 2019 Nano Energy 64 103909Google Scholar

    [2]

    Teh Y S, Bhattacharya K 2019 J. Appl. Phys. 125 064103Google Scholar

    [3]

    Pal S, Swain A B, Biswas P P, Murali D, Pal A, Nanda B R K, Murugavel P 2018 Sci. Rep. 8 17Google Scholar

    [4]

    Butler K T, Frost J M, Walsh A 2015 Energy Environ. Sci. 8 838848Google Scholar

    [5]

    蔡田怡, 雎胜 2018 物理学报 67 157801Google Scholar

    Cai T Y, Ju S 2018 Acta Phys. Sin. 67 157801Google Scholar

    [6]

    王婧, 吴霞, 邓朝勇, 朱孔军, 南策文 2014 无机材料学报 29 905911Google Scholar

    Wang J, Wu X, Deng C Y, Zhu K J, Nan C W 2014 J. Inorg. Mater. 29 905911Google Scholar

    [7]

    Chen Y, Chen J, Yang S, Li Y, Gao X, Zeng M, Fan Z, Gao X, Lu X, Liu J M 2018 Mater. Res. Bull. 107 456Google Scholar

    [8]

    Young S M, Zheng F, Rappe A M 2015 Phys. Rev. Appl. 4 054004Google Scholar

    [9]

    张兴良 2010 硕士学位论文 (武汉: 武汉理工大学)

    Zhang X L 2010 M. S. Thesis (Wuhan: Wuhan University Of Technology) (in Chinese)

    [10]

    Sinha K, Zhang Y, Jiang X, Wang H, Wang X, Zhang X, Ryan P J, Kim J W, Bowlan J, Yarotski D A, Li Y, DiChiara A D, Cheng X, Wu X, Xu X 2017 Phys. Rev. B 95 094110Google Scholar

    [11]

    Han H, Kim D, Chae S, Park J, Nam S Y, Choi M, Yong K, Kim H J, Son J, Jang H M 2018 Nanoscale 10 13261Google Scholar

    [12]

    Han H, Kim D, Chu K, Park J, Nam S Y, Heo S, Yang C, Jang H M 2018 ACS Appl. Mater. Interfaces 10 18461853Google Scholar

    [13]

    Akbashev A R, Semisalova A S, Perov N S, Kaul A R 2011 Appl. Phys. Lett. 99 122502Google Scholar

    [14]

    Wang W, Zhao J, Wang W, et al. 2013 Rev. Lett. 110 237601Google Scholar

    [15]

    Huang X, Paudel T R, Dong S, Tsymbal E Y 2015 Phys. Rev. B 92 125201Google Scholar

    [16]

    Lin L, Zhang H M, Liu M F, Shen S, Zhou S, Li D, Wang X, Yan Z B, Zhang Z D, Zhao J, Dong S, Liu J M 2016 Phys. Rev. B 93 075146Google Scholar

    [17]

    Liu J, Sun T L, Liu X Q, Tian H, Gao T T, Chen X M 2018 Adv. Funct. Mater. 28 1706062Google Scholar

    [18]

    Fu Z, Nair H S, Xiao Y, Senyshyn A, Pomjakushin V, Feng E, Pomjakushin V, Su Y, Jin W T, Bruckel T 2016 Phys. Rev. B 94 125150Google Scholar

    [19]

    Clark S J, Segall M D, Pickard C J, Hasnip P, Probert M I, Refson K, Payne M C 2005 Z. Kristallogr. 220 567

    [20]

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

    [21]

    Holinsworth B S, Mazumdar D, Brooks C M, Mundy J A, Das H, Cherian J G, McGill S A, Fennie C J, Schlom D G, Musfeldt J L 2015 Appl. Phys. Lett. 106 082902Google Scholar

    [22]

    Ridzwan M H, Yaakob M K, Taib M F M, Ali A M M, Hassan O H, Yahya M Z A 2017 Mater. Res. Express 4 044001Google Scholar

    [23]

    Disseler S M, Borchers J A, Brooks C M, Mundy J A, Moyer J A, Hillsberry D A, Thies E L, Tenne D A, Heron J, Holtz M, Clarkson J D, Stiehl G M, Schiffer P, Muller D A, Schlom D G, Ratcliff W D 2015 Phys. Rev. Lett. 114 217602Google Scholar

    [24]

    Wang W, Wang H, Xu X, Zhu L, He L, Wills E, Cheng X, Keavney D J, Shen J, Wu X, Xu X 2012 Appl. Phys. Lett. 101 241907Google Scholar

    [25]

    Das H, Wysocki A L, Geng Y, Wu W, Fennie C J 2014 Nat. Commun. 5 3998Google Scholar

    [26]

    Tu S, Zhang Y, Reshak A H, Auluck S, Ye L, Han X, Ma T, Huang H 2019 Nano Energy 56 840Google Scholar

    [27]

    Roy A, Mukherjee S, Gupta R, Auluck S, Prasad R, Garg A 2011 J. Phys. Condens. Matter 23 325902Google Scholar

  • 图 1  h-Lu1–xInxFeO3的球棍结构模型图 (a) P1, P2位置; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

    Fig. 1.  Model of the ball-and-stick structure of h-Lu1–xInxFeO3: (a) P1, P2 position; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

    图 2  h-Lu1–xFexO3的能带图 (a) h-Lu0.833In0.167FeO3; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

    Fig. 2.  Energy band diagrams of h-Lu1–xFexO3: (a) h-Lu0.833In0.167FeO3; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

    图 3  分布态密度图 (a) 未掺杂的h-LuFeO3; (b) h-Lu0.333In0.667FeO3

    Fig. 3.  Distribution density of states: (a) Undoped h-LuFeO3; (b) h-Lu0.333In0.667FeO3.

    图 4  In掺杂前后h-LuFeO3光学吸收系数随入射光子能量的变化

    Fig. 4.  Change of optical absorption coefficient of h-LuFeO3 with incident photon energy before and after In doping.

    图 5  不同In/Lu比的h-Lu1–xInxFeO3极化值(红色曲线)和晶格常数c/a比(蓝色曲线)

    Fig. 5.  Polarization values (red curve) and lattice constant c/a ratios (blue curve) of h-Lu1–xInxFeO3 with different In/Lu ratios.

    表 1  不同磁序下不同位置的In掺杂相对能量变化

    Table 1.  Relative energy changes of In doping at different positions under different magnetic orders.

    掺杂位置G型/eVC型/eVA型/eV铁磁型/eV
    未掺00.011.370.01
    P1 位置010.383.1412.51
    P2 位置00.030.251.41
    下载: 导出CSV

    表 2  In∶Lu为1∶2和2∶1时h-Lu1–xInxO3不同磁序的相对能量

    Table 2.  Relative energy of h-Lu1–xInxO3 with different magnetic sequencewhen In∶Lu is 1∶2 and 2∶1.

    材料G型/eVC型/eVA型/eV铁磁型/eV
    h-Lu2/3In1/3FeO300.01260.86040.8604
    h-Lu1/3In2/3FeO300.00960.89520.9487
    下载: 导出CSV

    表 3  h-Lu1–xInxO3的结构优化结果

    Table 3.  Structure optimization results of h-Lu1–xInxO3.

    材料晶格常数/Å体积/Å3轴角/(°)
    abcαβγ
    (h-LuFeO3)[15]5.9655.96511.702
    (h-LuFeO3)[23]5.9855.98511.770
    h-LuFeO36.0676.06711.756374.92690.00090.000119.993
    h-Lu0.833In0.167FeO36.0426.04211.880374.77090.04790.004120.214
    h-Lu0.667In0.333FeO36.0056.00711.935372.22190.01490.033120.166
    h-Lu0.333In0.667FeO35.9225.92312.119367.51689.99990.001120.173
    下载: 导出CSV

    表 4  本工作带隙计算结果与已发表结果对比

    Table 4.  Comparison of calculated band gap results with published results.

    CASTEP[22]WIEN2K[24]VASP[25]本工作
    交换关
    联泛函
    LDAGGA-
    PBE
    GGA-
    PBEsol
    GGA-
    PBE
    U34.54.614.5
    带隙/eV0.541.11.351.16
    下载: 导出CSV
  • [1]

    Ji Y, Gao T, Wang Z L, Yang Y 2019 Nano Energy 64 103909Google Scholar

    [2]

    Teh Y S, Bhattacharya K 2019 J. Appl. Phys. 125 064103Google Scholar

    [3]

    Pal S, Swain A B, Biswas P P, Murali D, Pal A, Nanda B R K, Murugavel P 2018 Sci. Rep. 8 17Google Scholar

    [4]

    Butler K T, Frost J M, Walsh A 2015 Energy Environ. Sci. 8 838848Google Scholar

    [5]

    蔡田怡, 雎胜 2018 物理学报 67 157801Google Scholar

    Cai T Y, Ju S 2018 Acta Phys. Sin. 67 157801Google Scholar

    [6]

    王婧, 吴霞, 邓朝勇, 朱孔军, 南策文 2014 无机材料学报 29 905911Google Scholar

    Wang J, Wu X, Deng C Y, Zhu K J, Nan C W 2014 J. Inorg. Mater. 29 905911Google Scholar

    [7]

    Chen Y, Chen J, Yang S, Li Y, Gao X, Zeng M, Fan Z, Gao X, Lu X, Liu J M 2018 Mater. Res. Bull. 107 456Google Scholar

    [8]

    Young S M, Zheng F, Rappe A M 2015 Phys. Rev. Appl. 4 054004Google Scholar

    [9]

    张兴良 2010 硕士学位论文 (武汉: 武汉理工大学)

    Zhang X L 2010 M. S. Thesis (Wuhan: Wuhan University Of Technology) (in Chinese)

    [10]

    Sinha K, Zhang Y, Jiang X, Wang H, Wang X, Zhang X, Ryan P J, Kim J W, Bowlan J, Yarotski D A, Li Y, DiChiara A D, Cheng X, Wu X, Xu X 2017 Phys. Rev. B 95 094110Google Scholar

    [11]

    Han H, Kim D, Chae S, Park J, Nam S Y, Choi M, Yong K, Kim H J, Son J, Jang H M 2018 Nanoscale 10 13261Google Scholar

    [12]

    Han H, Kim D, Chu K, Park J, Nam S Y, Heo S, Yang C, Jang H M 2018 ACS Appl. Mater. Interfaces 10 18461853Google Scholar

    [13]

    Akbashev A R, Semisalova A S, Perov N S, Kaul A R 2011 Appl. Phys. Lett. 99 122502Google Scholar

    [14]

    Wang W, Zhao J, Wang W, et al. 2013 Rev. Lett. 110 237601Google Scholar

    [15]

    Huang X, Paudel T R, Dong S, Tsymbal E Y 2015 Phys. Rev. B 92 125201Google Scholar

    [16]

    Lin L, Zhang H M, Liu M F, Shen S, Zhou S, Li D, Wang X, Yan Z B, Zhang Z D, Zhao J, Dong S, Liu J M 2016 Phys. Rev. B 93 075146Google Scholar

    [17]

    Liu J, Sun T L, Liu X Q, Tian H, Gao T T, Chen X M 2018 Adv. Funct. Mater. 28 1706062Google Scholar

    [18]

    Fu Z, Nair H S, Xiao Y, Senyshyn A, Pomjakushin V, Feng E, Pomjakushin V, Su Y, Jin W T, Bruckel T 2016 Phys. Rev. B 94 125150Google Scholar

    [19]

    Clark S J, Segall M D, Pickard C J, Hasnip P, Probert M I, Refson K, Payne M C 2005 Z. Kristallogr. 220 567

    [20]

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

    [21]

    Holinsworth B S, Mazumdar D, Brooks C M, Mundy J A, Das H, Cherian J G, McGill S A, Fennie C J, Schlom D G, Musfeldt J L 2015 Appl. Phys. Lett. 106 082902Google Scholar

    [22]

    Ridzwan M H, Yaakob M K, Taib M F M, Ali A M M, Hassan O H, Yahya M Z A 2017 Mater. Res. Express 4 044001Google Scholar

    [23]

    Disseler S M, Borchers J A, Brooks C M, Mundy J A, Moyer J A, Hillsberry D A, Thies E L, Tenne D A, Heron J, Holtz M, Clarkson J D, Stiehl G M, Schiffer P, Muller D A, Schlom D G, Ratcliff W D 2015 Phys. Rev. Lett. 114 217602Google Scholar

    [24]

    Wang W, Wang H, Xu X, Zhu L, He L, Wills E, Cheng X, Keavney D J, Shen J, Wu X, Xu X 2012 Appl. Phys. Lett. 101 241907Google Scholar

    [25]

    Das H, Wysocki A L, Geng Y, Wu W, Fennie C J 2014 Nat. Commun. 5 3998Google Scholar

    [26]

    Tu S, Zhang Y, Reshak A H, Auluck S, Ye L, Han X, Ma T, Huang H 2019 Nano Energy 56 840Google Scholar

    [27]

    Roy A, Mukherjee S, Gupta R, Auluck S, Prasad R, Garg A 2011 J. Phys. Condens. Matter 23 325902Google Scholar

计量
  • 文章访问数:  942
  • PDF下载量:  48
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-07
  • 修回日期:  2020-09-07
  • 上网日期:  2021-01-24
  • 刊出日期:  2021-02-05

In掺杂h-LuFeO3光吸收及极化性能的第一性原理计算

  • 湘潭大学材料科学与工程学院, 湘潭 411105
  • 通信作者: 钟向丽, xlzhong@xtu.edu.cn
    基金项目: 国家自然科学基金(批准号: 51872251, 11875229)和电子元器件可靠性物理及其应用技术重点实验室开放基金(批准号: ZHD201803)资助的课题

摘要: h-LuFeO3是一种窄带隙铁电半导体材料, 已被证明在铁电光伏领域有较好的应用前景. 然而, 较低的极化强度使光生电子-空穴对复合率大, 限制了h-LuFeO3基铁电光伏电池效率的提高. 为改善h-LuFeO3的极化强度, 提高光吸收性质, 本文利用第一性原理计算方法研究了In原子在h-LuFeO3不同位置的掺杂形成能, 得到最稳定的掺杂位置, 比较了h-Lu1–xInxFeO3 (x = 0, 0.167, 0.333, 0.667)的带隙、光吸收性能及极化强度等性质. 计算结果表明, 随着In掺杂比例的增加, h-LuFeO3的晶格常数c/a比不断增大, 铁电极化强度得到提高. 当In∶Lu = 2∶1时, 材料杂质能级出现, Fe-O轨道杂化得到增强, 提高了h-LuFeO3的光吸收性能. 此工作证明了In掺杂是改善h-LuFeO3极化强度和光吸收系数的有效方法, 对铁电光伏性能的提高提供一种新途径.

English Abstract

    • 铁电材料中的铁电极化可作为内建电场促进光生载流子的分离, 使铁电光伏电池能够实现光学带隙以上的开路电压[1-4]. 这种特殊的铁电光伏效应, 为突破传统p-n结型光伏电池的理论转换效率瓶颈提供了全新的方法, 在新能源领域表现出巨大的应用潜力[5]. 研究人员对众多铁电材料, 例如钛酸钡(BaTiO3)、锆钛酸铅(Pb(ZrxTi1–x)O3)、铌酸锂(LiNiO3)等的光伏性质进行了大量研究[6-8], 然而, 较大的带隙使得大部分太阳光能量无法被吸收(太阳光能量范围约为0.4—4 eV)[9]. h-LuFeO3是一种新型多铁材料, 光学带隙在2.0 eV左右[10], 在铁电光伏领域有较好的应用前景. 由于本征h-LuFeO3较低的极化强度(仅为4.7 μC·cm–2)[11], 使得h-LuFeO3基铁电光伏电池的光生载流子复合率较大, 不利于电池光电转换效率的提高. 目前实验上制备的h-LuFeO3基铁电光伏电池的光电转化效率在0.001%左右[12]. 先前的研究表明, 提高h-LuFeO3的极化强度能够有效降低光生载流子的复合率, 增加电池的开路电压[5], 改善其光吸收性质, 提高h-LuFeO3基铁电光伏电池的性能.

      室温条件下h-LuFeO3的铁电性不稳定, 可以通过应力使其稳定存在[13,14]. 然而采用外延应力调控h-LuFeO3铁电性能对实验条件的要求十分苛刻[15], 极大地限制了应力对h-LuFeO3的调控效果. 研究发现, 采用掺杂的方法不仅能对h-LuFeO3性能进行调控, 且能制备出铁电性在室温下稳定存在的h-LuFeO3[16-18]. 该方法对实验条件依赖小, 能够增强材料的铁电光伏相关性能. 在Lu位置置换离子半径较小的离子, 是获得稳定六角结构LuFeO3的一种重要方法. Lin等[16]通过第一性原理计算表明, Sc掺杂可以使LuFeO3更稳定, 而其多铁性不受影响. 但In离子半径比Lu和Sc都要小, 理论上InFeO3也是六角结构. Liu等[17]利用In掺杂h-LuFeO3能得到六角结构的h-Lu1–xInxFeO3 材料, X射线全谱拟合结果表明, 随着In含量的增加, 材料结构逐渐由极性的P 63cm (x = 0.4—0.6)转变为非极性的P 63/mmc (x = 0.75), 即顺电相. 而目前并未深入研究In掺杂对h-LuFeO3光学性能的影响. 本文基于密度泛函理论(density functional theory, DFT)的第一性原理计算方法对h-Lu1–xInxFeO3 (x = 0, 0.167, 0.333, 0.667)的磁结构、带隙、光吸收系数及自发极化等性质进行了研究, 揭示In掺杂对h-LuFeO3的极化性质及光学性质的影响, 为提高h-LuFeO3的铁电光伏性能提供可靠的理论依据.

    • 使用的第一性原理计算方法均基于Materials Studio软件的CASTEP模块[10], 采用广义梯度近似(generalized gradient approximation, GGA)和平面赝势波的方法处理电子与离子间的相互作用, 利用GGA下的Perdew-Burke-Ernzerhof (PBE)泛函作为交换关联泛函[19]. 赝势选取Ultrasolft, 截断能为572 eV, 计算中所有晶胞和原子都进行弛豫, 以达到最稳定状态. 电子自洽精度为 10–5 eV/atom, 原子受力小于0.05 eV/Å, 晶胞结构优化时第一布里渊区的k点采用4 × 4 × 2, 电子结构优化时k点采用6 × 6 × 3. 对于磁性元素Fe采用GGA+U的方法来处理, 其中U是位库仑势. 参考先前的报道, U值选择为4.5[20].

    • In在h-LuFeO3的掺杂位置有两种, 分别记为P1和P2[21], 如图1(a)所示. 在包含30个原子的晶胞模型中, 有4个P1位原子和2个P2位原子, 首先讨论h-Lu1–xFexO3最稳定的磁序. 采用A型、C型、G型三种反铁磁结构和铁磁结构对相同的晶胞模型进行优化. 为了得到最稳定的晶胞结构, 以G型反铁磁序h-Lu0.833In0.167FeO3能量为基准, 比较了不同掺杂位置的In原子替换单个Lu原子对h-LuFeO3磁序的影响, 并比较了相对总能量, 结果如表1所列. 结果表明, 四种磁序中非线性G型反铁磁序的h-Lu0.833In0.167FeO3能量最低, 这与先前报道的计算结论一致[22], 且无论是P1位掺杂还是P2位掺杂均没有改变磁序. 利用(1)式计算了不同掺杂位置的掺杂形成能:

      图  1  h-Lu1–xInxFeO3的球棍结构模型图 (a) P1, P2位置; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

      Figure 1.  Model of the ball-and-stick structure of h-Lu1–xInxFeO3: (a) P1, P2 position; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

      ${E_{\rm{f}}} = E{}_{{\rm{doped}}} - {E_{{\rm{pristine}}}} - {\mu _{{\rm{Lu}}}} + {\mu _{{\rm{In}}}},$

      其中Ef为In掺杂形成能, EdopedEpristine分别代表掺杂后和掺杂前的h-LuFeO3总能, μLuμIn分别为Lu和In的化学势. 形成能反映了掺杂的难易程度, 形成能越小越容易掺杂. 通过计算发现, P1位置的掺杂形成能为5.426 eV, P2位置的掺杂形成能为5.581 eV, 表明P1位置相对P2位置更容易掺杂.

      掺杂位置G型/eVC型/eVA型/eV铁磁型/eV
      未掺00.011.370.01
      P1 位置010.383.1412.51
      P2 位置00.030.251.41

      表 1  不同磁序下不同位置的In掺杂相对能量变化

      Table 1.  Relative energy changes of In doping at different positions under different magnetic orders.

      由于本征InFeO3稳定相是顺电相P 63/mmc, 计算时控制了In∶Lu最高为2∶1. h-LuFeO3的P1位置比P2位置更容易掺杂, 故均选择In原子替换P1位原子, In∶Lu分别为1∶2和2∶1时, 对应的晶胞优化结果如图1(b)图1(c)所示. 以非线性G型反铁磁序的h-Lu0.833In0.167FeO3能量为基准, 比较了不同磁序的晶胞结构优化后相对能量的变化, 结果如表2所列. 随着掺杂浓度的提高, h-Lu1–xInxFeO3 (x = 0.333, 0.667)磁序未发生变化, 非线性G型反铁磁序的h-Lu0.833In0.167FeO3能量最低.

      材料G型/eVC型/eVA型/eV铁磁型/eV
      h-Lu2/3In1/3FeO300.01260.86040.8604
      h-Lu1/3In2/3FeO300.00960.89520.9487

      表 2  In∶Lu为1∶2和2∶1时h-Lu1–xInxO3不同磁序的相对能量

      Table 2.  Relative energy of h-Lu1–xInxO3 with different magnetic sequencewhen In∶Lu is 1∶2 and 2∶1.

      确定了最稳定的磁序后, 计算P1位的In掺杂对h-LuFeO3晶格参数的影响, 与先前实验报道的h-LuFeO3晶格参数进行对比, 结果如表3所列, 表明非线性G型反铁磁序h-Lu0.833In0.167FeO3结构优化结果与实验值符合较好. 随着In原子含量的提高, h-Lu1–xInxFeO3晶胞沿c轴拉伸, c/a比从1.94增加到2.04, 这有利于材料自发极化率的增强. 而晶胞略微缩小, 我们猜测可能是由于In的离子半径比Lu离子半径小造成的.

      材料晶格常数/Å体积/Å3轴角/(°)
      abcαβγ
      (h-LuFeO3)[15]5.9655.96511.702
      (h-LuFeO3)[23]5.9855.98511.770
      h-LuFeO36.0676.06711.756374.92690.00090.000119.993
      h-Lu0.833In0.167FeO36.0426.04211.880374.77090.04790.004120.214
      h-Lu0.667In0.333FeO36.0056.00711.935372.22190.01490.033120.166
      h-Lu0.333In0.667FeO35.9225.92312.119367.51689.99990.001120.173

      表 3  h-Lu1–xInxO3的结构优化结果

      Table 3.  Structure optimization results of h-Lu1–xInxO3.

    • h-Lu0.833In0.167FeO3, h-Lu0.667In0.333FeO3以及h-Lu0.333In0.667FeO3的能带图见图2. 从图2(a)可知, h-Lu0.833In0.167FeO3的带隙大小为1.16 eV, 小于实验值(2.0 eV)[24], 与先前文献中利用了Local Density Approximation (LDA), Perdew-Burke-Ernzerhof for solid (PBEsol)等方法计算的h-LuFeO3带隙对比如表4所列. 带隙偏小是因为利用DFT方法会低估材料带隙, 然而掺杂计算的主要目的是为了预测相对变化趋势, 而非研究绝对值的大小, 因此选取DFT方法并不影响最终结论.

      图  2  h-Lu1–xFexO3的能带图 (a) h-Lu0.833In0.167FeO3; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

      Figure 2.  Energy band diagrams of h-Lu1–xFexO3: (a) h-Lu0.833In0.167FeO3; (b) h-Lu0.667In0.333FeO3; (c) h-Lu0.333In0.667FeO3

      CASTEP[22]WIEN2K[24]VASP[25]本工作
      交换关
      联泛函
      LDAGGA-
      PBE
      GGA-
      PBEsol
      GGA-
      PBE
      U34.54.614.5
      带隙/eV0.541.11.351.16

      表 4  本工作带隙计算结果与已发表结果对比

      Table 4.  Comparison of calculated band gap results with published results.

      三种计算模型的h-Lu1–xInxFeO3导带底和价带顶均在G点, 表明h-Lu1–xInxFeO3为直接带隙半导体材料. 如图2(b)所示, 当In∶Lu = 1∶2时, 即In占据晶胞中两个P1位置时, 带隙减小到1.05 eV, 表明In原子的替换可减小材料的带隙. 由图2(c)可知, 当In∶Lu = 2∶1时, In完全占据4个P1位置, 带隙的变化并不明显, 仅下降了0.01 eV; 但此时导带底(1—2 eV范围内)杂质能级的增多有利于提升材料的光吸收系数.

      为了深入研究In掺杂对h-LuFeO3带隙两端电子跃迁的影响, 分析了未掺杂h-LuFeO3h-Lu0.333In0.667FeO3在带隙附近, 即能量在–2— 2 eV间的态密度图, 如图3所示. 由图3(a)可知, 未掺杂h-LuFeO3的导带底(1—2 eV范围内)大部分由Fe 3d轨道占据, 小部分由O 2p轨道占据; 价带顶(–1—0 eV范围内)主要是O 2p轨道占据, 而Lu对价带顶和导带底的贡献有限, 从主要占据价带顶的O 2p轨道电子转变为主要占据导带底的Fe 3d轨道电子, 表现出Fe 3d轨道与O 2p轨道的高度杂化. 因此调控Fe 3d轨道与O 2p轨道杂化程度可以影响材料的光吸收性能.

      图  3  分布态密度图 (a) 未掺杂的h-LuFeO3; (b) h-Lu0.333In0.667FeO3

      Figure 3.  Distribution density of states: (a) Undoped h-LuFeO3; (b) h-Lu0.333In0.667FeO3.

      图3(b)h-Lu0.333In0.667FeO3的态密度图, 导带底(1—2 eV范围内)大部分由Fe 3d轨道占据, 小部分由O 2p轨道及In 4s轨道占据; 价带顶(–1—0 eV范围内)主要由O 2p轨道占据. 相对未掺杂的h-LuFeO3, h-Lu0.333In0.667FeO3杂质能级的出现及导带底Fe 3d轨道赝势峰强度的增大, 使材料带隙缩小, 接收价带电子的能力增强. 分析结果表明, In掺杂影响了材料Fe-O间的轨道杂化, 提高了h-Lu0.333In0.667FeO3的光吸收系数.

    • 图4为光沿c轴入射h-Lu1–xInxFeO3的光吸收系数变化情况, 发现在1.64—2.4 eV范围内, 掺杂后的光吸收系数比未掺杂的h-LuFeO3略微减小. 但在能量大于2.4 eV区域, 材料光吸收系数随着In/Lu比值的增大而增大, 表明In掺杂能够有效提高h-LuFeO3在太阳光能量范围内的光吸收效率.

      图  4  In掺杂前后h-LuFeO3光学吸收系数随入射光子能量的变化

      Figure 4.  Change of optical absorption coefficient of h-LuFeO3 with incident photon energy before and after In doping.

      铁电极化是影响铁电光伏性能的另一个重要性质. 铁电光伏材料以自身铁电极化为内建电场, 增大铁电极化能够减少载流子复合概率[26], 因此缩小带隙和增大铁电极化强度是提高铁电光伏材料性能的两大研究方向. 根据伯恩有效电荷定性计算自发极化强度计算公式为[27]

      $\Delta {P_\alpha } \cong \sum\limits_{j\beta } {\frac{{\partial {P_\alpha }}}{{\partial {u_{j\beta }}}}} \left( {{u_{j\beta }} - {u_{0j\beta }}} \right) = \frac{e}{\varOmega }\sum\limits_{j\beta } {Z_{_{j\alpha \beta }}^{^{\rm{*}}}} \Delta {u_{j\beta }},$

      其中e表示电子电荷量, $Z_{_{j\alpha \beta }}^{^{\rm{*}}}$表示伯恩有效电荷, Δu表示原子相对位移, Ω表示晶胞总体积. 将h-Lu1–xInxFeO3沿c轴的伯恩有效电荷代入(2)式, 计算得到h-LuFeO3, h-Lu0.833In0.167FeO3, h-Lu0.667In0.333FeO3h-Lu0.333In0.667FeO3沿c轴的铁电极化强度, 分别为3.93, 5.91, 7.92和11.02 μC·cm–2, 如图5中红色曲线所示. 图5中蓝色曲线表示随In/Lu比的增大, h-Lu1–xInxFeO3晶胞的晶格参数c/a比值的变化. 可见材料的极化值与c/a比有相同的变化趋势, 结果表明在保持铁电相的前提下, 提高材料的c/a比可提高材料的极化强度. 对铁电光伏材料而言, 极化强度的增大能够得到更高的光生电压及光生电流[11]. 此外, 改变c/a比相当于材料受力沿c轴产生拉应变, 而拉应变能够减小六角铁酸盐类多铁材料的带隙并增强光吸收性质[13].

      图  5  不同In/Lu比的h-Lu1–xInxFeO3极化值(红色曲线)和晶格常数c/a比(蓝色曲线)

      Figure 5.  Polarization values (red curve) and lattice constant c/a ratios (blue curve) of h-Lu1–xInxFeO3 with different In/Lu ratios.

    • 本文基于DFT的第一性原理计算方法, 研究了h-Lu1–xInxFeO3 (x = 0, 0.167, 0.333, 0.667)的带隙、极化强度及光吸收性质. 研究结果表明, In原子在h-LuFeO3晶胞中会优先替换P1位置, 随着In掺杂量的增加, h-Lu1–xInxFeO3的晶胞沿c轴拉伸. 当x达到0.667时, 晶格常数c/a比从1.94增大至2.04. 通过分析h-LuFeO3h-Lu0.333In0.667FeO3的分态密度可知, c/a比的增大能够增强h-Lu0.333In0.667FeO3层间的Fe-O轨道杂化程度, 提高在太阳光范围内的光吸收系数. 利用伯恩有效电荷定性地计算h-In1–xLuxFeO3的极化强度, 发现In掺杂能够提高材料沿c轴的极化强度, 且极化值的变化趋势与晶格常数c/a比的变化趋势相同. 因此, In掺杂能够有效提高h-LuFeO3的光吸收系数及铁电极化强度, 为提高铁电光伏性能提供重要理论指导.

参考文献 (27)

目录

    /

    返回文章
    返回