搜索

x

留言板

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

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

钛酸钡的光学性质及其体积效应

孙智征 荀威 张加永 刘传洋 仲嘉霖 吴银忠

引用本文:
Citation:

钛酸钡的光学性质及其体积效应

孙智征, 荀威, 张加永, 刘传洋, 仲嘉霖, 吴银忠

Optical properties of BaTiO3 and its volume effects

Sun Zhi-Zheng, Xun Wei, Zhang Jia-Yong, Liu Chuan-Yang, Zhong Jia-Lin, Wu Yin-Zhong
PDF
HTML
导出引用
  • 钛酸钡(BaTiO3, BTO)是铁电物理学和材料学领域具有代表性的研究对象. 本文基于准粒子相互作用的GW方法研究BTO材料的光学性质, 并研究了等应变情况下的体积效应. 第一性计算结果表明, 考虑了激子效应的GW (格林函数(G)-库仑势(W))方法得到的激发态性质更接近实验结果. 引入等应变调控, 发现体积膨胀会导致光学吸收峰红移, 体积压缩则光学吸收峰蓝移. 同时, 探究了体积变化对BTO块材的自发极化和电子结构的影响, 发现体积膨胀会使钛原子的d轨道和氧原子的p轨道杂化更显著, 进一步导致材料自发极化的增大, 而体积压缩对自发极化和dp杂化的影响正好相反. 通过比较研究, 还发现等应变的体积效应对极化的影响不如等应力体积效应明显.
    BaTiO3 (BTO) is a typical studying object both in ferroelectrics and in material science. By the GW method, optical property of BTO is investigated, and its volume effect under the case of iso-strain is also studied. It is found that the results of excited states are closer to the experimental results with the consideration of electron-hole interaction in the framework of GW method. Considering the volume effect, we obtain that the red shift of the peaks of optical absorption occurs under the expansion of volume, and the blue shift appears when the BTO is compressed. At the same time, the polarization and the hybridization between d orbital of Ti atom and p orbital of O atom are enhanced for the case of volume expansion, however, things will be opposite under the compression of volume. Furthermore, the volume effect in the iso-strain case is less dramatic than in the iso-stress case.
      通信作者: 吴银忠, yzwu@usts.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11274054)、江苏省十三五重点学科项目(批准号: 20168765)和苏州市低维光电材料与器件重点实验室(批准号: SZS201611)资助的课题.
      Corresponding author: Wu Yin-Zhong, yzwu@usts.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11274054), the Jiangsu Key Disciplines of the Thirteenth Five-Year Plan, China(Grant No. 20168765), and the Suzhou Key Laboratory for Low Dimensional Optoelectronic Materials and Devices, China (Grant No. SZS201611).
    [1]

    钟维烈 1998 铁电体物理学(北京: 科学出版社)

    Zhong W L 1998 Ferroelectric Physics (Beijing: Science Press) (in Chinese)

    [2]

    Cohen R E 1992 Nature 358 136Google Scholar

    [3]

    Ma C, Jin K J, Ge C, Yang G Z 2018 Phys. Rev. B 97 115103Google Scholar

    [4]

    Zhuravlev M Y, Sabirianov R F, Jaswal S S, Tsymbal E Y 2005 Phys. Rev. Lett. 94 246802Google Scholar

    [5]

    Yu Z L, Ma Q R, Zhao Y Q, Liu B, Cai M Q 2018 J. Phys. Chem. C 122 9275

    [6]

    Khenata R, Sahnoun M, Baltache H, Rerat M, Rashek A H, Illes N, Bouhafs B 2005 Solid State Commun. 136 120Google Scholar

    [7]

    Hosseini S M, Movlarooy T, Kompany A 2005 Euro. Phys. J B 46 463Google Scholar

    [8]

    Gupta G, Nautiyal T, Auluck S 2004 Phys. Rev. B 69 052101

    [9]

    Saha S, Sinha T P, Mookerjee A 2000 Phys. Rev. B 62 8828Google Scholar

    [10]

    Liu Q J, Zhang N C, Liu F S, Wang H Y, Liu Z T 2013 Opt. Mater. 35 2629Google Scholar

    [11]

    Wemple S H 1970 Phys. Rev. B 2 2679Google Scholar

    [12]

    Piskunov S, Heifets E, Eglitis R I, Borstel G 2004 Comp. Mater. Sci. 29 165Google Scholar

    [13]

    Hybertsen M S, Louie S G 1986 Phy. Rev. B 34 5390Google Scholar

    [14]

    Rohlfing M, Louie S G 2000 Phy. Rev. B 62 4927

    [15]

    Deslippe J, Samsonidze, Strubbe D A, Jain M, Cohen M L, Louie S G 2012 Comp. Mater. Sci. 183 1269

    [16]

    Sanna S, Thierfelder C, Wippermann S, Sinha T P, Schmidt W G 2011 Phys. Rev. B 83 054112Google Scholar

    [17]

    Kornev I A, Bellaiche L, Bouvier P, Janolin P E, Dkhil B, Kreisel J 2005 Phys. Rev. Lett. 95 196804Google Scholar

    [18]

    Moriwake H, Koyama Y, Matsunaga K, Hirayama T, Tanaka I 2008 J Phys.: Condens. Mat. 20 345207Google Scholar

    [19]

    Wang X X, Wang C L 2003 Chin. J. Chem. 21 1130

    [20]

    Nguyen M D, Dekkers M, Houwman E, Steenwelle R, Wan X, Roelofs A, Schmitz-Kempen T, Rijnders G 2011 Appl. Phys. Lett. 99 252904Google Scholar

    [21]

    Wen Z, Qiu X B, Li C, Zheng C Y, Ge X Y, Li A D, Wu D 2014 Appl. Phys. Lett. 104 042907Google Scholar

    [22]

    Li W, Weng G J 2001 J. Appl. Phys. 90 2484Google Scholar

    [23]

    Dieguez O, Tinte S, Antons A, Bungaro C, Neaton J B, Rabe K M, Vanderbilt D 2004 Phys. Rev. B 69 212101Google Scholar

    [24]

    Zhao R, Li W W, Lee J H, Choi E M, Liang Y, Zhang W, Tang R J, Wang H Y, Jia Q X, MacManus-Driscoll J L, Yang H 2014 Adv. Funct. Mater. 24 5240Google Scholar

    [25]

    Van Helvoort A T J, Da H L, Soleim B G, Holmestad R, Tybell T 2005 Appl. Phys. Lett. 86 092907Google Scholar

    [26]

    Chen A, Hu J M, Lu P, Yang T N, Zhang W R, Li L G, Ahmed T, Enriquez E, Weigand W, Su Q, Wang H Y, Zhu J Y, MacManus-Driscoll J L, Chen L Q, Yarotski D, Jia Q X 2016 Sci. Adv. 2 e1600245Google Scholar

    [27]

    Kohn W, Sham L J 1965 Phys. Rev. 140 A1133Google Scholar

    [28]

    Perdew J P, Yue W 1986 Phys. Rev. B 33 8820

    [29]

    Gao S Y, Yang L 2017 Phys. Rev. B 96 155410Google Scholar

    [30]

    Ehrenreich H, Cohen M H 1959 Phys. Rev. 115 1786

    [31]

    Bilc D I, Orlando R, Shaltaf R, Rignanese G M, Íñiguez J, Ghosez P 2008 Phys. Rev. B 77 165107Google Scholar

    [32]

    Wang F G, Grinberg I, Rappe A M 2014 Appl. Phys. Lett. 104 152903Google Scholar

    [33]

    Dong H F, Wu Z G, Wang S Y, Duan W H, Li J B 2013 Appl. Phys. Lett. 102 072905Google Scholar

    [34]

    Chernova E, Pacherova O, Chvostova D, Dejneka A, Kocourek T, Jelinek M, Tyunina M 2015 Appl. Phys. Lett. 106 192903Google Scholar

  • 图 1  BTO块材自发极化随体积的变化曲线

    Fig. 1.  Spontaneous polarization of bulk BTO versus its volume

    图 2  体积变化对Ti原子d轨道分波态密度的影响

    Fig. 2.  Projected DOS of ${\rm d}_{xy}$ and ${\rm d}_{xz}$ of Ti atom for different volume

    图 3  体积改变对BTO能带的影响

    Fig. 3.  Energy bands of BTO for selected volumes.

    图 4  不同体积BTO的介电常数虚部谱($\epsilon_2^{\bot}$为垂直于极化方向即xy平面内的介电常数虚部, $\epsilon_2^{^{_{/\!/}}}$z方向的介电常数虚部)

    Fig. 4.  Spectrum of imagine part of dielectric constant for selected BTO volumes, where $\epsilon_2^{\bot}$ corresponds the direction perpendicular to the polarization, and $\epsilon_2^{^{_{/\!/}}}$ stands for the direction parallel to the polarization

    图 5  不同体积下BTO块材的光学吸收谱($\alpha^{\bot}$为垂直于极化方向即xy平面内的光学吸收系数, $\alpha^{^{_{/\!/}}}$z方向的吸收系数)

    Fig. 5.  Spectrum of optical absorption of BTO for selected volumes, where $\alpha^{\bot}$ corresponds the direction perpendicular to the polarization, and $\alpha^{^{_{/\!/}}}$ stands for the direction parallel to the polarization

    表 1  不同体积下计算得到的BTO块材的能隙(单位:eV)

    Table 1.  Energy gaps (in eV) for selected BTO volumes and different theoretical method

    方法 V = 0.9V0 V = V0 V = 1.1V0
    PBE 1.89 1.76, 1.71[16] 1.66
    HSE 3.42 3.25, 3.12[32] 3.10
    GW 3.65 3.54, 3.90[16] 3.41
    实验 3.4[11]
    下载: 导出CSV
  • [1]

    钟维烈 1998 铁电体物理学(北京: 科学出版社)

    Zhong W L 1998 Ferroelectric Physics (Beijing: Science Press) (in Chinese)

    [2]

    Cohen R E 1992 Nature 358 136Google Scholar

    [3]

    Ma C, Jin K J, Ge C, Yang G Z 2018 Phys. Rev. B 97 115103Google Scholar

    [4]

    Zhuravlev M Y, Sabirianov R F, Jaswal S S, Tsymbal E Y 2005 Phys. Rev. Lett. 94 246802Google Scholar

    [5]

    Yu Z L, Ma Q R, Zhao Y Q, Liu B, Cai M Q 2018 J. Phys. Chem. C 122 9275

    [6]

    Khenata R, Sahnoun M, Baltache H, Rerat M, Rashek A H, Illes N, Bouhafs B 2005 Solid State Commun. 136 120Google Scholar

    [7]

    Hosseini S M, Movlarooy T, Kompany A 2005 Euro. Phys. J B 46 463Google Scholar

    [8]

    Gupta G, Nautiyal T, Auluck S 2004 Phys. Rev. B 69 052101

    [9]

    Saha S, Sinha T P, Mookerjee A 2000 Phys. Rev. B 62 8828Google Scholar

    [10]

    Liu Q J, Zhang N C, Liu F S, Wang H Y, Liu Z T 2013 Opt. Mater. 35 2629Google Scholar

    [11]

    Wemple S H 1970 Phys. Rev. B 2 2679Google Scholar

    [12]

    Piskunov S, Heifets E, Eglitis R I, Borstel G 2004 Comp. Mater. Sci. 29 165Google Scholar

    [13]

    Hybertsen M S, Louie S G 1986 Phy. Rev. B 34 5390Google Scholar

    [14]

    Rohlfing M, Louie S G 2000 Phy. Rev. B 62 4927

    [15]

    Deslippe J, Samsonidze, Strubbe D A, Jain M, Cohen M L, Louie S G 2012 Comp. Mater. Sci. 183 1269

    [16]

    Sanna S, Thierfelder C, Wippermann S, Sinha T P, Schmidt W G 2011 Phys. Rev. B 83 054112Google Scholar

    [17]

    Kornev I A, Bellaiche L, Bouvier P, Janolin P E, Dkhil B, Kreisel J 2005 Phys. Rev. Lett. 95 196804Google Scholar

    [18]

    Moriwake H, Koyama Y, Matsunaga K, Hirayama T, Tanaka I 2008 J Phys.: Condens. Mat. 20 345207Google Scholar

    [19]

    Wang X X, Wang C L 2003 Chin. J. Chem. 21 1130

    [20]

    Nguyen M D, Dekkers M, Houwman E, Steenwelle R, Wan X, Roelofs A, Schmitz-Kempen T, Rijnders G 2011 Appl. Phys. Lett. 99 252904Google Scholar

    [21]

    Wen Z, Qiu X B, Li C, Zheng C Y, Ge X Y, Li A D, Wu D 2014 Appl. Phys. Lett. 104 042907Google Scholar

    [22]

    Li W, Weng G J 2001 J. Appl. Phys. 90 2484Google Scholar

    [23]

    Dieguez O, Tinte S, Antons A, Bungaro C, Neaton J B, Rabe K M, Vanderbilt D 2004 Phys. Rev. B 69 212101Google Scholar

    [24]

    Zhao R, Li W W, Lee J H, Choi E M, Liang Y, Zhang W, Tang R J, Wang H Y, Jia Q X, MacManus-Driscoll J L, Yang H 2014 Adv. Funct. Mater. 24 5240Google Scholar

    [25]

    Van Helvoort A T J, Da H L, Soleim B G, Holmestad R, Tybell T 2005 Appl. Phys. Lett. 86 092907Google Scholar

    [26]

    Chen A, Hu J M, Lu P, Yang T N, Zhang W R, Li L G, Ahmed T, Enriquez E, Weigand W, Su Q, Wang H Y, Zhu J Y, MacManus-Driscoll J L, Chen L Q, Yarotski D, Jia Q X 2016 Sci. Adv. 2 e1600245Google Scholar

    [27]

    Kohn W, Sham L J 1965 Phys. Rev. 140 A1133Google Scholar

    [28]

    Perdew J P, Yue W 1986 Phys. Rev. B 33 8820

    [29]

    Gao S Y, Yang L 2017 Phys. Rev. B 96 155410Google Scholar

    [30]

    Ehrenreich H, Cohen M H 1959 Phys. Rev. 115 1786

    [31]

    Bilc D I, Orlando R, Shaltaf R, Rignanese G M, Íñiguez J, Ghosez P 2008 Phys. Rev. B 77 165107Google Scholar

    [32]

    Wang F G, Grinberg I, Rappe A M 2014 Appl. Phys. Lett. 104 152903Google Scholar

    [33]

    Dong H F, Wu Z G, Wang S Y, Duan W H, Li J B 2013 Appl. Phys. Lett. 102 072905Google Scholar

    [34]

    Chernova E, Pacherova O, Chvostova D, Dejneka A, Kocourek T, Jelinek M, Tyunina M 2015 Appl. Phys. Lett. 106 192903Google Scholar

  • [1] 梁帅博, 袁涛, 邱扬, 张震, 妙亚宁, 韩竞峰, 刘秀童, 姚春丽. 钛酸钡介电调控提升纸基摩擦纳米发电机输出性能. 物理学报, 2022, 71(7): 077701. doi: 10.7498/aps.71.20212022
    [2] 黄申洋, 张国伟, 汪凡洁, 雷雨晨, 晏湖根. 二维黑磷的光学性质. 物理学报, 2021, 70(2): 027802. doi: 10.7498/aps.70.20201497
    [3] 王娇, 刘少辉, 陈长青, 郝好山, 翟继卫. 钛酸钡基/聚偏氟乙烯复合介质材料的界面改性与储能性能. 物理学报, 2020, 69(21): 217702. doi: 10.7498/aps.69.20201031
    [4] 张丽丽, 夏桐, 刘桂安, 雷博程, 赵旭才, 王少霞, 黄以能. 第一性原理方法研究N-Pr共掺杂ZnO的电子结构和光学性质. 物理学报, 2019, 68(1): 017401. doi: 10.7498/aps.68.20181531
    [5] 景奇, 李晓娟. 多孔钛酸钡陶瓷制备及其增强的压电灵敏性. 物理学报, 2019, 68(5): 057701. doi: 10.7498/aps.68.20181790
    [6] 金芹, 董海明, 韩奎, 王雪峰. 石墨烯超快动态光学性质. 物理学报, 2015, 64(23): 237801. doi: 10.7498/aps.64.237801
    [7] 吴琼, 刘俊, 董前民, 刘阳, 梁培, 舒海波. 硫化锡电子结构和光学性质的量子尺寸效应. 物理学报, 2014, 63(6): 067101. doi: 10.7498/aps.63.067101
    [8] 冯现徉, 逯瑶, 蒋雷, 张国莲, 张昌文, 王培吉. In掺杂ZnO超晶格光学性质的研究. 物理学报, 2012, 61(5): 057101. doi: 10.7498/aps.61.057101
    [9] 高尚鹏, 祝桐. 基于自洽GW方法的碳化硅准粒子能带结构计算. 物理学报, 2012, 61(13): 137103. doi: 10.7498/aps.61.137103
    [10] 孙中华, 王红艳, 张志东, 张中月. 金纳米环结构的光学性质研究. 物理学报, 2011, 60(4): 047808. doi: 10.7498/aps.60.047808
    [11] 焦照勇, 杨继飞, 张现周, 马淑红, 郭永亮. 闪锌矿GaN弹性性质、电子结构和光学性质外压力效应的理论研究. 物理学报, 2011, 60(11): 117103. doi: 10.7498/aps.60.117103
    [12] 李旭珍, 谢泉, 陈茜, 赵凤娟, 崔冬萌. OsSi2电子结构和光学性质的研究. 物理学报, 2010, 59(3): 2016-2021. doi: 10.7498/aps.59.2016
    [13] 邢海英, 范广涵, 赵德刚, 何 苗, 章 勇, 周天明. Mn掺杂GaN电子结构和光学性质研究. 物理学报, 2008, 57(10): 6513-6519. doi: 10.7498/aps.57.6513
    [14] 段满益, 徐 明, 周海平, 陈青云, 胡志刚, 董成军. 碳掺杂ZnO的电子结构和光学性质. 物理学报, 2008, 57(10): 6520-6525. doi: 10.7498/aps.57.6520
    [15] 关 丽, 刘保亭, 李 旭, 赵庆勋, 王英龙, 郭建新, 王书彪. 萤石结构TiO2的电子结构和光学性质. 物理学报, 2008, 57(1): 482-487. doi: 10.7498/aps.57.482
    [16] 倪建刚, 刘 诺, 杨果来, 张 曦. 第一性原理研究BaTiO3(001)表面的电子结构. 物理学报, 2008, 57(7): 4434-4440. doi: 10.7498/aps.57.4434
    [17] 沈益斌, 周 勋, 徐 明, 丁迎春, 段满益, 令狐荣锋, 祝文军. 过渡金属掺杂ZnO的电子结构和光学性质. 物理学报, 2007, 56(6): 3440-3445. doi: 10.7498/aps.56.3440
    [18] 王培吉, 周忠祥, 苏 燕, 荣振宇, 赵 朋, 张奉军. 钽掺杂对钛酸钡导热性能影响的研究. 物理学报, 2006, 55(4): 1959-1964. doi: 10.7498/aps.55.1959
    [19] 伞海生, 李 斌, 冯博学, 何毓阳, 陈 冲. 由缺陷引起的Burstein-Moss和带隙收缩效应对CdIn2O4透明导电薄膜光带隙的影响. 物理学报, 2005, 54(2): 842-847. doi: 10.7498/aps.54.842
    [20] 李智强, 陆夏莲, 陈敏, 何山, 李景德. 钙钛矿结构中的简谐子软模. 物理学报, 2002, 51(7): 1581-1585. doi: 10.7498/aps.51.1581
计量
  • 文章访问数:  4249
  • PDF下载量:  45
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-24
  • 修回日期:  2019-02-04
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-04-20

钛酸钡的光学性质及其体积效应

  • 1. 苏州科技大学数理学院, 苏州 215009
  • 2. 苏州科技大学电子与信息工程学院, 苏州 215009
  • 通信作者: 吴银忠, yzwu@usts.edu.cn
    基金项目: 国家自然科学基金(批准号: 11274054)、江苏省十三五重点学科项目(批准号: 20168765)和苏州市低维光电材料与器件重点实验室(批准号: SZS201611)资助的课题.

摘要: 钛酸钡(BaTiO3, BTO)是铁电物理学和材料学领域具有代表性的研究对象. 本文基于准粒子相互作用的GW方法研究BTO材料的光学性质, 并研究了等应变情况下的体积效应. 第一性计算结果表明, 考虑了激子效应的GW (格林函数(G)-库仑势(W))方法得到的激发态性质更接近实验结果. 引入等应变调控, 发现体积膨胀会导致光学吸收峰红移, 体积压缩则光学吸收峰蓝移. 同时, 探究了体积变化对BTO块材的自发极化和电子结构的影响, 发现体积膨胀会使钛原子的d轨道和氧原子的p轨道杂化更显著, 进一步导致材料自发极化的增大, 而体积压缩对自发极化和dp杂化的影响正好相反. 通过比较研究, 还发现等应变的体积效应对极化的影响不如等应力体积效应明显.

English Abstract

    • 钛酸钡(BaTiO3, BTO)是典型的钙钛矿型铁电材料, 室温下结构为四方铁电相, 相变温度为393 K[1]. 由于其在电子学、电力转换、非线性光学和非挥发性存储器件等领域有广阔的应用前景, 因此实验和理论上BTO块材[2]和掺杂体系[3]以及由BTO垒组成的量子隧道结[4]一直是比较活跃的研究领域, 包括由此延伸出来的杂化钙钛矿材料的光学性质[5]. 本文主要研究四方相BTO体材料的光学性质并考察体积效应对材料自发极化和光学性质的影响. 之前理论上研究BTO材料的光学性质, 主要是利用密度泛函理论(DFT)计算出材料的介电常数和光学吸收系数[6-10]. 从太阳能应用角度来讲, BTO是宽禁带光伏材料, 其禁带宽度约为3.5 eV[11], DFT理论计算往往低估BTO材料的能隙, 而杂化泛函方法也不能得到与实验相吻合的晶格参数[12]. 考虑了准粒子相互作用的格林函数-库仑势 (GW)方法[13-15] 是把DFT的计算结果作为零阶近似, 把准粒子之间的相互作用作为一级近似, 可以得到更为准确的激发态信息和能隙, 能进一步通过解Bethe-Salpeter方程考虑激子效应. 因此, 在计算资源充足的情况下, GW方法是一个非常有效的理论方法, 对铁电材料计算得到的能隙也与实验结果很接近[16].

      铁电材料的铁电相的结构一般与高对称顺电相差别很小, 因此铁电材料的体积效应会非常显著. BTO材料的铁电性主要源于氧八面体中心的Ti 离子偏移中心对称位置, BTO被称为位移型铁电体. 从电子结构角度来说, BTO材料的自发极化主要源于Ti的3d 电子与O的2p 电子的杂化. 对铁电材料体积效应以及应变效应的报道很多, 大致可以分为以下三个方向的研究. 第一为通过加静水压的方法研究体积效应, 即考虑铁电晶体 在$ x $, $ y $, $ z $三个方向受到相等的应力, 研究静水压对铁电材料的结构[17]、自发极化[18]、电子结构[19]和光学性质[6]的影响. 第二类为仅通过基底的失配对铁电薄膜在$ xy $平面内形成晶格失配, 从而在薄膜的$ xy $平面内形成张应力或压应力[20,21]. 这个方向的实验研究比较容易实现. 理论上可以通过朗道唯象理论加以研究[22], 也有大量第一性原理计算的研究报道[23], 定性的结论为, $ xy $平面内压应力会增强$ z $方向的自发极化, 而张应力则会抑制$ z $方向的极化. 第三个方向为通过施加垂直应变来研究$ z $方向的应变对铁电材料极化的影响, 垂直应变可以通过压电响应力显微镜(PFM)探针直接施加于材料表面, 也可以通过垂直方向的晶格失配来调控[24-26]. 从实验技术的角度, 目前已经可以做到分别在$ xy $平面内和$ z $方向对铁电材料施加应变. 静水压实现的是三维等应力的效果, 本文则考虑在$ x $, $ y $, $ z $ 三个方向等应变的情况下, BTO体积变化对材料的极化、电子结构以及光学性质的影响. 同时, 本文采用的计算方法是考虑了准粒子相互作用的GW方法, 因此理论计算结果可以更好地预测材料的激发态性质.

    • GW方法是基于准粒子概念和格林函数的多体微扰理论[13-15], 与之对应的DFT理论(Kohn-Sham方法)只能得到比较准确的材料基态性质[27,28]. 格林函数(G)是描述准粒子的首选模型, 用屏蔽库仑势($ W $)描述体系准粒子之间的相互作用, 自能算符描述电子气中交换关联效应, 这就是GW方法名称的来源. GW方法可以得到更为准确的激发态性质, 已经成功应用于计算半导体材料以及新型功能二维材料的电子结构和能带[29]. 目前GW方法已经被纳入到第一性原理计算的VASP和ABINT软件包中, 其主要计算过程分为以下几个步 骤: 1)进行DFT计算, 得到体系的基态波函数和本征值; 2)基于上一步的结果, 利用无规相近似(RPA)计算屏蔽库伦势$ W $, 考虑单体格林函数, 获得准粒子的能量; 3)基于单体格林函数和准粒子能量, 考虑两体格林函数$ \rm G_2 $, 也就是求解Bethe-Salpeter 方程, 得到修正的电子空穴激发态.

    • VASP软件已经包含了与频率相关的介电常数的计算, $ \epsilon(\omega) = \epsilon_1(\omega)+{\rm i}\epsilon_2(\omega) $, 其中介电常数的虚部为[30]:

      $ \epsilon_2(\omega) = \frac{e^2\hbar}{{\text{π}} m^2\omega^2}\sum_{v, c}\int_{\rm BZ}|{M_{cv}(k)}|^2\delta[\omega_{cv}(k)-\omega]{\rm d}^3k. $

      积分遍及第一布里渊区, $ {M_{cv}(k)} = \left\langle u_{ck}|e\nabla|u_{vk}\right\rangle $是价带和导带之间的直接跃迁矩阵, $ \hbar\omega_{cv}(k) = E_{ck}- $$E_{vk}$是激发能量. 介电常数的实部$ \epsilon_1(\omega) $可以由Kramers-Kronig关系式给出[6]:

      $ \epsilon_1(\omega) = 1+\frac{{\text{π}}}{2}P\int_{0}^{\infty}\frac{\omega^\prime\epsilon_2(\omega^\prime)}{\omega^{\prime^2}-\omega^2}{\rm d}\omega^\prime. $

      基于介电常数的实部和虚部, 可以计算出材料的光学吸收谱[9]:

      $ \alpha(\omega) = \frac{\sqrt2\omega}{c}\sqrt{\epsilon_1(\omega)^2+\epsilon_2(\omega)^2-\epsilon_1( \omega)}. $

    • 本文以典型四方相BTO材料为例, 利用VASP软件包进行第一性计算, 势文件选取PBE (Perdew-Burke-Ernzerhof)势, 各原子参与自洽计算的电子如下: Ba 5s25p66s2; Ti 3s23p63d24s2; O 2s22p4. 平面波截断能量 $ E_{\rm cut} $ = 550 eV, 空带数$ {N_{\rm Bands}} $ = 384, 弛豫时原子受力标准Hellmann-Feynman力小于0.01 eV/Å. 布里渊区K点取样为$8\times 8\times 8$. 在GW和RPA计算中, 响应函数的存储和处理占据了主要的计算时间, 因此, 综合考虑计算的精确度和有限的计算资源, 选取能量截断值ENCUTGW = 150 eV, 总态数NBANDSGW = 30. 在BSE (Bethe-Salpeter equation)计算中, 由于计算时间与参与计算的总态数(占据态和空态)成三次方关系, 一般仅取靠近费米能级附近的量子态, 因此设置占据态数NBANDSO = 20 (文中参与计算的体系总电子数为40)和空态数NBANDSV = 10, 这样计入了全部占据态和10个空态来计算光学性质.

    • 在未施加应变的情况下, 弛豫得到的四方相BTO晶格常数为: $ a $ = $ b $ = 3.998 Å, $ c $ = 4.237 Å, $ V_0 = abc $, 这与文献[31]报道的理论结果一致. 利用Berry phase方法, 我们计算得到了BTO块材的自发极化$ P_0 $ $= 0.46\; {\rm C/m^2} $. 本文采用等应变的方式来研究体积效应, 即令晶格常数$ a $, $ b $, $ c $按同等比例压缩或者伸展(iso-strain). 图1中取了5个体积数据点来计算和分析, 分别是$ V/V_0 $ = 0.9, 0.95, 1.00, 1.05和1.10, 其中$ V_0 $为BTO未施加应变时的体积. 为了与等应力情形(iso-stress)进行比较, 同时选择了14, 5, –3和–4 GPa四种静水压力, 经过计算发现体系的体积变化正好为–10%, –5%, +5%, +10%. 由图1可以看出, 等应变或等应力的情形, BTO自发极化随着体积的增大而增大, 且等应变情况下的极化随体积的变化(实线)比等应力的情形(短划线)要缓慢些.

      图  1  BTO块材自发极化随体积的变化曲线

      Figure 1.  Spontaneous polarization of bulk BTO versus its volume

      BTO的极化主要来源于Ti 3d电子和O 2p电子的杂化, 对于四方相BTO, 2018年中国科学院物理研究所杨国桢院士课题组[3]用Ti的dxy和dyz, dxz轨道态密度的不对称性来表征掺杂BTO体系铁电性的强弱, 本文也采用这一表征方法观察体积变化对Ti原子的3d轨道的分波态密度dxy和dyz的影响(四方相的dxz和dyz态密度是对称的, 因此只要选择其中一个与dxy进行对比研究). 文献对BTO的体积效应的研究大多集中于等应力条件下[18, 19], 即用静水压方式实现. 如无特殊说明, 下文中涉及的体积效应均指等应变情形下的体积效应. 不同体积下四方相BTO中的Ti的3d电子分波态密度见图2. 图2(a), (b), (c)分别对应于$ V/V_0 $ = 0.9, $ V/V_0 $ = 1.0, $ V/V_0 $ = 1.1的情况, 不难发现Ti 3d 电子在费米面以下主要分布在–4 — –2 eV, 在费米面以上主要分布于2—4 eV. 其中实线对应于dxy轨道上的电子态密度, 短划线表示dxz轨道上的电子态密度. 比较图2(a), (b), (c)中的实线与短划线, 一方面随着体积的增加, 导带中的dxy和dxz峰值位置向低能区移动; 另一方面dxy和dxz的DOS (density of states) 分布错开得越来越明显, 也就是dxz分布相对于dxy分布位置蓝移得越来越多. 因为立方相BTO的dxy和dxz的态密度是重叠的, 因此根据文献[3]的观点, 随着体积的膨胀, BTO材料的铁电性增强, 这与图1中得到的自发极化增大的结论一致.

      图  2  体积变化对Ti原子d轨道分波态密度的影响

      Figure 2.  Projected DOS of ${\rm d}_{xy}$ and ${\rm d}_{xz}$ of Ti atom for different volume

    • 目前VASP软件包中的DFT计算可以直接给出体系的能带, 但VASP包中的GW计算要得到特定路径下的能带需要通过Wannier函数插值. 图3分别给出了不同体积下体系的GW能带图(实线)和DFT能带图(短划线). 由图3可以看出, 无论是DFT方法还是考虑了准粒子相互作用的GW方法, 计算得到的BTO能带均是间接带隙($A{\text{-}}\varGamma $). 首先分析无应变下体系的能带图(图3(b)), 其中DFT和GW方法得到的价带图像基本一致, 而对于导带位置, GW方法得到的导带能量蓝移了近1.8 eV ($ \varGamma{}$点). 在无应变情况下, GW方法得到的能隙为$ E_{\rm g} $ = 3.54 eV, 因此在DFT理论的基础上进一步考虑准粒子相互作用, 可以得到与实验结果(3.4 eV)[11]更加吻合的能隙, 且我们的理论计算结果与文献[16] 基本一致. 比较图3(a), (b), (c), 不难发现随着体积的压缩, 导带和价带的形状基本保持不变, 只是能隙变大, 这是由于随着体积的压缩, 导带底Ti的d轨道向高能区移动所导致的(亦可见图2), 具体的能隙变化整理于表1中. 表1中也同时给出了杂化泛函HSE的计算结果. 相对于PBE, 杂化泛函方法HSE 对能隙有了很好的修正, 但PBE和HSE均是基于DFT理论, GW则是基于多体格林函数方法, GW本质上应更加符合物理实际. 本文应用GW 方法研究了钙钛矿材料的等应变体积效应, 虽然定性的体积效应结论与DFT理论一致, 但由于GW 方法对激发态的计算更加准确, 因此可以给出更精确的光学性质. 下一节将会讨论BTO块材的光学性质及其体积效应.

      图  3  体积改变对BTO能带的影响

      Figure 3.  Energy bands of BTO for selected volumes.

      方法 V = 0.9V0 V = V0 V = 1.1V0
      PBE 1.89 1.76, 1.71[16] 1.66
      HSE 3.42 3.25, 3.12[32] 3.10
      GW 3.65 3.54, 3.90[16] 3.41
      实验 3.4[11]

      表 1  不同体积下计算得到的BTO块材的能隙(单位:eV)

      Table 1.  Energy gaps (in eV) for selected BTO volumes and different theoretical method

    • 应变调控可以有效提高光电材料的光学吸收效率[33], 对BTO材料的应变实验研究表明[34], 当BTO外延生长在SrTiO$ _3 $ 基底时, BTO在$ xy $平面内获取的应变大致在0.3%—1.3%. 实验发现, 由于应变的引入, 原体材料在5 eV处的吸收峰发生了0.2—0.4 eV 的蓝移, 而在8.5 eV附近的吸收峰几乎没有受到应变的影响. 在理论处理方法上, 针对立方相和四方相BTO光学性质的第一性计算研究报道见文献[6, 8−10, 16], 其中前四篇是基于DFT的计算, 最后一篇是利用GW + BSE方法分别给出了BTO的介电常数和光学吸收谱. 本文考虑在等应变情形下改变四方相BTO的体积, 研究体积对材料光学性质的影响. 图4图5分别给出了四方相BTO材料在V/V0 = 0.9, 1.0, 1.1情况下体系的介电常数虚部和光学吸收系数, 其中的蓝色实线是DFT-PBE的计算结果, 黑色实线代表的是考虑了准粒子相互作用的GW计算得到的数据结果, 红色实线表示的是在GW理论框架下进一步考虑电子-空穴相互作用(激子效应), 即求解BSE方程得到的四方相BTO的光学性质.

      图  4  不同体积BTO的介电常数虚部谱($\epsilon_2^{\bot}$为垂直于极化方向即xy平面内的介电常数虚部, $\epsilon_2^{^{_{/\!/}}}$z方向的介电常数虚部)

      Figure 4.  Spectrum of imagine part of dielectric constant for selected BTO volumes, where $\epsilon_2^{\bot}$ corresponds the direction perpendicular to the polarization, and $\epsilon_2^{^{_{/\!/}}}$ stands for the direction parallel to the polarization

      图  5  不同体积下BTO块材的光学吸收谱($\alpha^{\bot}$为垂直于极化方向即xy平面内的光学吸收系数, $\alpha^{^{_{/\!/}}}$z方向的吸收系数)

      Figure 5.  Spectrum of optical absorption of BTO for selected volumes, where $\alpha^{\bot}$ corresponds the direction perpendicular to the polarization, and $\alpha^{^{_{/\!/}}}$ stands for the direction parallel to the polarization

      通过对图4图5的分析, 得到的主要结果和结论如下. 1)无论是介电常数虚部还是光学吸收系数, 从起始峰角度来看, 当V = V$ _0 $ 时(无应变), 考虑了准粒子相互作用的GW 方法计算结果相对于DFT理论结果发生了2 eV的蓝移, 而激子效应引起起始峰相对于GW结果发生了1 eV左右的红移, 计算结果表明, 考虑了激子效应的理论计算结果与实验结果[11]更为接近. 2)由于四方相BTO原子结构和电子结构在平行于$ z $轴和垂直于$ z $轴呈现各向异性, 因此介电常数的虚部和光学吸收系数在两个方向上也不等价. 通过对照研究, 沿着$ z $方向的自发极化会使材料在$ z $方向产生内建电场, 从而导致介电常数虚部和光学吸收系数的起始峰在平行方向相对于垂直方向发生大约0.5 eV的蓝移. 3)体积压缩会使介电常数虚部和光学吸收系数发生蓝移(起始峰的蓝移现象最为显著), 而体积膨胀则导致介电常数虚部和光学吸收系数的红移, 我们的计算结果定性上与加静水压得到的结果一致[7].

    • 应用GW理论方法, 对BTO块材的电子结构和光学性质进行了详细研究, 同时考察了等应变情况下的体积效应. 研究发现: 体积膨胀会增大材料的自发极化, 增强Ti的d轨道和O的p轨道的杂化, 使材料的能隙降低, 光学吸收起始峰发生红移; 而体积压缩发生的现象与体积膨胀正好相反. 由于GW方法可以更好地描述材料的激发态性质, 原则上可以得到更加精确的光学性质及其体积效应. 在综合考虑计算资源的前提下, 如何进一步考虑声子对光学跃迁的贡献是下一步比较有意义的研究方向.

      感谢苏州大学雎胜教授对GW计算方法有益的交流和讨论.

参考文献 (34)

目录

    /

    返回文章
    返回