搜索

x

留言板

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

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

新型层状Bi2Se3的第一性原理研究

郭宇 周思 赵纪军

新型层状Bi2Se3的第一性原理研究

郭宇, 周思, 赵纪军
PDF
HTML
导出引用
导出核心图
  • 近年来, 在石墨烯研究热潮的推动下, 众多种类丰富、性能各异的二维化合物材料相继被发现, 其中一些二维材料具有多种同素异构体, 进而呈现出更丰富的性质. 层状Bi2Se3由于其独特的物理性质, 受到人们广泛的关注, 而它的同素异构体尚未有人研究. 本文采用基于密度泛函理论的结构搜索方法, 预测了一个稳定的β-Bi2Se3新相, 它具有良好的动力学和热力学稳定性, 并在低Bi2Se3源化学势条件下容易形成. 单层β-Bi2Se3是一个直接带隙为2.44 eV的二维半导体, 其电子载流子有效质量低至0.52m0, 在可见光范围内具有高达105 cm–1的光吸收系数, 并且能带边缘位置适中, 可用于光催化水分解制氢气. 此外, 由于β-Bi2Se3在垂直层面方向的镜面对称性破缺, 能够产生面外极化强度, 具有0.58 pm/V的面外压电系数. 鉴于其新颖的电子特性, 二维β-Bi2Se3在未来的电子器件中可能发挥重要的作用.
      通信作者: 周思, sizhou@dlut.edu.cn
    • 基金项目: 中国博士后科学基金(批准号: BX20190052, 2020M670739)、国家自然科学基金(批准号: 11974068)和中央高校基础研究经费(批准号: DUT20LAB110)资助的课题
    [1]

    Zhang H, Liu C X, Qi X L, Dai X, Fang Z, Zhang S C 2009 Nat. Phys. 5 438

    [2]

    Kong D, Chen Y, Cha J J, Zhang Q, Analytis J G, Lai K, Liu Z, Hong S S, Koski K J, Mo S K 2011 Nat. Nanotechnol. 6 705

    [3]

    Brom J E, Ke Y, Du R, Won D, Weng X, Andre K, Gagnon J C, Mohney S E, Li Q, Chen K 2012 Appl. Phys. Lett. 100 162110

    [4]

    Alegria L D, Schroer M D, Chatterjee A, Poirier G R, Pretko M, Patel S K, Petta J R 2012 Nano Lett. 12 4711

    [5]

    Alegria L D, Petta J R 2012 Nanotechnology 23 435601

    [6]

    Le P H, Wu K H, Luo C W, Leu J 2013 Thin Solid Films 534 659

    [7]

    Hirahara T, Sakamoto Y, Takeichi Y, Miyazaki H, Kimura S, Matsuda I, Kakizaki A, Hasegawa S 2010 Phys. Rev. B 82 155309

    [8]

    Yu X, He L, Lang M, Jiang W, Xiu F, Liao Z, Wang Y, Kou X, Zhang P, Tang J 2012 Nanotechnology 24 015705

    [9]

    Li Y Y, Wang G, Zhu X G, Liu M H, Ye C, Chen X, Wang Y Y, He K, Wang L L, Ma X C 2010 Adv. Mater. 22 4002

    [10]

    Xia Y, Qian D, Hsieh D, Wray L, Pal A, Lin H, Bansil A, Grauer D, Hor Y S, Cava R J 2009 Nat. Phys. 5 398

    [11]

    Bansal N, Koirala N, Brahlek M, Han M G, Zhu Y, Cao Y, Waugh J, Dessau D S, Oh S 2014 Appl. Phys. Lett. 104 241606

    [12]

    Chen S, Zhao C, Li Y, Huang H, Lu S, Zhang H, Wen S 2014 Opt. Mater. Express 4 587

    [13]

    Sun Y, Cheng H, Gao S, Liu Q, Sun Z, Xiao C, Wu C, Wei S, Xie Y 2012 J. Am. Chem. Soc. 134 20294

    [14]

    Min Y, Park G, Kim B, Giri A, Zeng J, Roh J W, Kim S I, Lee K H, Jeong U 2015 ACS Nano 9 6843

    [15]

    Xu H, Chen G, Jin R, Chen D, Wang Y, Pei J, Zhang Y, Yan C, Qiu Z 2014 Crystengcomm 16 3965

    [16]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsov A A 2004 Science 306 666

    [17]

    Li Y, Xu L, Liu H, Li Y 2014 Chem. Soc. Rev. 43 2572

    [18]

    Li L, Yu Y, Ye G J, Ge Q, Ou X, Wu H, Feng D, Chen X H, Zhang Y 2014 Nat. Nanotechnol. 9 372

    [19]

    Qiao J, Kong X, Hu Z X, Yang F, Ji W 2014 Nat. Commun. 5 4475

    [20]

    Ghosh B, Nahas S, Bhowmick S, Agarwal A 2015 Phys. Rev. B 91 115433

    [21]

    Mogulkoc Y, Modarresi M, Mogulkoc A, Ciftci Y O 2016 Comput. Mater. Sci. 124 23

    [22]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169

    [23]

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

    [24]

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

    [25]

    Grimme S 2006 J. Comput. Chem. 27 1787

    [26]

    Grimme S, Antony J, Ehrlich S, Krieg H 2010 J. Chem. Phys. 132 154104

    [27]

    Togo A, Oba F, Tanaka I 2008 Phys. Rev. B 78 134106

    [28]

    Baroni S, De Gironcoli S, Dal Corso A, Giannozzi P 2001 Rev. Mod. Phys. 73 515

    [29]

    Barnett R N, Landman U 1993 Phys. Rev. B 48 2081

    [30]

    Martyna G J, Klein M L, Tuckerman M 1992 J. Chem. Phys. 97 2635

    [31]

    Wang Y, Lv J, Zhu L, Ma Y 2010 Phys. Rev. B 82 094116

    [32]

    Wang Y, Miao M, Lv J, Zhu L, Yin K, Liu H, Ma Y 2012 J. Chem. Phys. 137 224108

    [33]

    Han N, Liu H, Zhou S, Zhao J 2016 J. Phys. Chem. C 120 14699

    [34]

    Zhan L B, Yang C L, Wang M S, Ma X G 2020 Physica E 124 114272

    [35]

    Zhang Y, He K, Chang C Z, Song C L, Wang L L, Chen X, Jia J F, Fang Z, Dai X, Shan W Y, Shen S Q, Niu Q, Qi X L, Zhang S C, Ma X C, Xue Q K 2010 Nat. Phys. 6 584

    [36]

    Graziano G, Klimeš J, Fernandez Alonso F, Michaelides A 2012 J. Phys.-Condes. Matter 24 424216

    [37]

    Cai Y, Zhang G, Zhang Y W 2014 Sci. Rep. 4 6677

    [38]

    Chakrapani V, Angus J C, Anderson A B, Wolter S D, Stoner B R, Sumanasekera G U 2007 Science 318 1424

    [39]

    Zhuang H L, Hennig R G 2013 Chem. Mater. 25 3232

    [40]

    Ma Z, Zhuang J, Zhang X, Zhou Z 2018 Front. Phys. 13 138104

    [41]

    Zhang X, Zhang Z, Wu D, Zhang X, Zhao X, Zhou Z 2018 Small Methods 2 1700359

    [42]

    Beal A R, Hughes H P 1979 Solid State Phys. 12 881

    [43]

    Duerloo K N, Ong M T, Reed E J 2012 J. Phys. Chem. Lett. 3 2871

    [44]

    King Smith R D, Vanderbilt D 1993 Phys. Rev. B 47 1651

    [45]

    Hangleiter A, Hitzel F, Lahmann S, Rossow U 2003 Appl. Phys. Lett. 83 1169

    [46]

    Shimada K 2006 Jpn. J. Appl. Phys. 45 L358

    [47]

    Guo Y, Zhou S, Bai Y, Zhao J 2017 Appl. Phys. Lett. 110 163102

  • 图 1  (a) α-Bi2Se3的原子结构; (b)单层β-Bi2Se3结构的俯视图(上图)和侧视图(下图); (c)双层β-Bi2Se3结构的俯视图(上图)和侧视图(下图); (d)经过10 ps第一性原理分子动力学模拟, 得到了300 K时Bi2Se3单层的平衡结构; (e) β-Bi2Se3的声子谱; (f) β-Bi2Se3单层的电子局域函数

    Fig. 1.  (a) Atomic structure of α-Bi2Se3; (b) the top and side views of monolayer β-Bi2Se3; (c) the top and side views of bilayer β-Bi2Se3; (d) snapshots of the equilibrium structures of the β-Bi2Se3 monolayer at 300 K after 10 ps ab initio molecular dynamic simulation; (e) phonon dispersion of monolayer β-Bi2Se3; (f) electron localization function for monolayer β-Bi2Se3.

    图 A1  CALYPSO搜索得到的几个较低能量的Bi2Se3单层结构(a)及对应的声子谱(b), 其中Bi2Se3-1, Bi2Se3-2, Bi2Se3-3的形成能分别为–0.15, –0.12, –0.09 eV/atom

    Fig. A1.  Some typical low-energy structures (a) of freestanding Bi2Se3 monolayer predicted by the CALYPSO code and corresponding phonon dispersions (b). The formation energy of Bi2Se3-1, Bi2Se3-2, Bi2Se3-3 are –0.15, –0.12, –0.09 eV/atom respectively.

    图 A2  温度为300 K时β-Bi2Se3单层的能量-时间变化 曲线

    Fig. A2.  Variations of temperature and energy with the time of AIMD simulation for β-Bi2Se3 monolayer at 300 K.

    图 2  α-Bi2Se3β-Bi2Se3体系表面自由能的化学势相图

    Fig. 2.  Chemical potential phase diagram of surface free ener-gy for α-Bi2Se3 and β-Bi2Se3.

    图 3  (a)不考虑SOC和(b)考虑SOC时, 采用HSE06泛函计算得到的β-Bi2Se3的能带结构和LDOS

    Fig. 3.  The electronic band structures (left panel) and LDOS (right panel) (a) without and (b) with SOC effect for monolayer β-Bi2Se3 using HSE06 functional, respectively.

    图 4  (a)采用HSE06泛函并且考虑SOC效应的双层(左图)和块体(右图)β-Bi2Se3的能带结构; (b)单层β-Bi2Se3带隙随双轴应变的变化

    Fig. 4.  (a) The electronic band structures for bilayer (left panel) and bulk (right panel) β-Bi2Se3 based on HSE06 level with SOC effect; (d) effect of biaxial strain on band gap of monolayer β-Bi2Se3.

    图 A3  不同堆叠方式的双层β-Bi2Se3 (a)能量最低的β-Bi2Se3双层结构, 将它的能量设定为0 eV; (b)相对能量为0.32 eV;(c)相对能量为0.55 eV

    Fig. A3.  β-Bi2Se3 bilayer with different stacking types and their relative energies: (a) the atomic structure of β-Bi2Se3 bilayer with the lowest energy, and its energy is set to 0 eV; the bilayer structures with relative energies of 0.32 eV (b) and 0.55 eV (c), respectively.

    图 5  (a)单层β-Bi2Se3的VBM和CBM对比pH = 7和pH = 0的氧化还原电势; (b)单层β-Bi2Se3的光吸收系数, λ是波长, 虚线中间区域表示可见光区

    Fig. 5.  (a) The location of VBM and CBM relative to vacuum energy of monolayer β-Bi2Se3 at pH = 0 and 7; (b) optical absorption coefficient for monolayer β-Bi2Se3. λ is the wave length, and the area between the red and the purple represents the visible range

    表 1  单层、双层和块体β-Bi2Se3相对真空能级的价带顶VBM和导带底CBM, 空穴和电子沿着xy方向的有效质量(mxh, myh, mxe, mye). 载流子有效质量以自由电子的静止质量m0为单位

    Table 1.  The VBM and CBM related to vacuum level for monolayer, bilayer and bulk β-Bi2Se3, and the corresponding carrier effective mass. m0 is the electron rest mass.

    β-Bi2Se3VBM/eVCBM/eVmxhmxemyhmye
    Monolayer–5.82–3.437.880.705.690.66
    Bilayer–5.00–4.242.550.522.360.52
    Bulk0.650.630.650.67
    下载: 导出CSV
  • [1]

    Zhang H, Liu C X, Qi X L, Dai X, Fang Z, Zhang S C 2009 Nat. Phys. 5 438

    [2]

    Kong D, Chen Y, Cha J J, Zhang Q, Analytis J G, Lai K, Liu Z, Hong S S, Koski K J, Mo S K 2011 Nat. Nanotechnol. 6 705

    [3]

    Brom J E, Ke Y, Du R, Won D, Weng X, Andre K, Gagnon J C, Mohney S E, Li Q, Chen K 2012 Appl. Phys. Lett. 100 162110

    [4]

    Alegria L D, Schroer M D, Chatterjee A, Poirier G R, Pretko M, Patel S K, Petta J R 2012 Nano Lett. 12 4711

    [5]

    Alegria L D, Petta J R 2012 Nanotechnology 23 435601

    [6]

    Le P H, Wu K H, Luo C W, Leu J 2013 Thin Solid Films 534 659

    [7]

    Hirahara T, Sakamoto Y, Takeichi Y, Miyazaki H, Kimura S, Matsuda I, Kakizaki A, Hasegawa S 2010 Phys. Rev. B 82 155309

    [8]

    Yu X, He L, Lang M, Jiang W, Xiu F, Liao Z, Wang Y, Kou X, Zhang P, Tang J 2012 Nanotechnology 24 015705

    [9]

    Li Y Y, Wang G, Zhu X G, Liu M H, Ye C, Chen X, Wang Y Y, He K, Wang L L, Ma X C 2010 Adv. Mater. 22 4002

    [10]

    Xia Y, Qian D, Hsieh D, Wray L, Pal A, Lin H, Bansil A, Grauer D, Hor Y S, Cava R J 2009 Nat. Phys. 5 398

    [11]

    Bansal N, Koirala N, Brahlek M, Han M G, Zhu Y, Cao Y, Waugh J, Dessau D S, Oh S 2014 Appl. Phys. Lett. 104 241606

    [12]

    Chen S, Zhao C, Li Y, Huang H, Lu S, Zhang H, Wen S 2014 Opt. Mater. Express 4 587

    [13]

    Sun Y, Cheng H, Gao S, Liu Q, Sun Z, Xiao C, Wu C, Wei S, Xie Y 2012 J. Am. Chem. Soc. 134 20294

    [14]

    Min Y, Park G, Kim B, Giri A, Zeng J, Roh J W, Kim S I, Lee K H, Jeong U 2015 ACS Nano 9 6843

    [15]

    Xu H, Chen G, Jin R, Chen D, Wang Y, Pei J, Zhang Y, Yan C, Qiu Z 2014 Crystengcomm 16 3965

    [16]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsov A A 2004 Science 306 666

    [17]

    Li Y, Xu L, Liu H, Li Y 2014 Chem. Soc. Rev. 43 2572

    [18]

    Li L, Yu Y, Ye G J, Ge Q, Ou X, Wu H, Feng D, Chen X H, Zhang Y 2014 Nat. Nanotechnol. 9 372

    [19]

    Qiao J, Kong X, Hu Z X, Yang F, Ji W 2014 Nat. Commun. 5 4475

    [20]

    Ghosh B, Nahas S, Bhowmick S, Agarwal A 2015 Phys. Rev. B 91 115433

    [21]

    Mogulkoc Y, Modarresi M, Mogulkoc A, Ciftci Y O 2016 Comput. Mater. Sci. 124 23

    [22]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169

    [23]

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

    [24]

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

    [25]

    Grimme S 2006 J. Comput. Chem. 27 1787

    [26]

    Grimme S, Antony J, Ehrlich S, Krieg H 2010 J. Chem. Phys. 132 154104

    [27]

    Togo A, Oba F, Tanaka I 2008 Phys. Rev. B 78 134106

    [28]

    Baroni S, De Gironcoli S, Dal Corso A, Giannozzi P 2001 Rev. Mod. Phys. 73 515

    [29]

    Barnett R N, Landman U 1993 Phys. Rev. B 48 2081

    [30]

    Martyna G J, Klein M L, Tuckerman M 1992 J. Chem. Phys. 97 2635

    [31]

    Wang Y, Lv J, Zhu L, Ma Y 2010 Phys. Rev. B 82 094116

    [32]

    Wang Y, Miao M, Lv J, Zhu L, Yin K, Liu H, Ma Y 2012 J. Chem. Phys. 137 224108

    [33]

    Han N, Liu H, Zhou S, Zhao J 2016 J. Phys. Chem. C 120 14699

    [34]

    Zhan L B, Yang C L, Wang M S, Ma X G 2020 Physica E 124 114272

    [35]

    Zhang Y, He K, Chang C Z, Song C L, Wang L L, Chen X, Jia J F, Fang Z, Dai X, Shan W Y, Shen S Q, Niu Q, Qi X L, Zhang S C, Ma X C, Xue Q K 2010 Nat. Phys. 6 584

    [36]

    Graziano G, Klimeš J, Fernandez Alonso F, Michaelides A 2012 J. Phys.-Condes. Matter 24 424216

    [37]

    Cai Y, Zhang G, Zhang Y W 2014 Sci. Rep. 4 6677

    [38]

    Chakrapani V, Angus J C, Anderson A B, Wolter S D, Stoner B R, Sumanasekera G U 2007 Science 318 1424

    [39]

    Zhuang H L, Hennig R G 2013 Chem. Mater. 25 3232

    [40]

    Ma Z, Zhuang J, Zhang X, Zhou Z 2018 Front. Phys. 13 138104

    [41]

    Zhang X, Zhang Z, Wu D, Zhang X, Zhao X, Zhou Z 2018 Small Methods 2 1700359

    [42]

    Beal A R, Hughes H P 1979 Solid State Phys. 12 881

    [43]

    Duerloo K N, Ong M T, Reed E J 2012 J. Phys. Chem. Lett. 3 2871

    [44]

    King Smith R D, Vanderbilt D 1993 Phys. Rev. B 47 1651

    [45]

    Hangleiter A, Hitzel F, Lahmann S, Rossow U 2003 Appl. Phys. Lett. 83 1169

    [46]

    Shimada K 2006 Jpn. J. Appl. Phys. 45 L358

    [47]

    Guo Y, Zhou S, Bai Y, Zhao J 2017 Appl. Phys. Lett. 110 163102

  • 引用本文:
    Citation:
计量
  • 文章访问数:  701
  • PDF下载量:  91
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-31
  • 修回日期:  2020-09-19
  • 上网日期:  2021-01-09
  • 刊出日期:  2021-01-20

新型层状Bi2Se3的第一性原理研究

  • 大连理工大学物理学院, 三束材料改性教育部重点实验室, 大连 116024
  • 通信作者: 周思, sizhou@dlut.edu.cn
    基金项目: 中国博士后科学基金(批准号: BX20190052, 2020M670739)、国家自然科学基金(批准号: 11974068)和中央高校基础研究经费(批准号: DUT20LAB110)资助的课题

摘要: 近年来, 在石墨烯研究热潮的推动下, 众多种类丰富、性能各异的二维化合物材料相继被发现, 其中一些二维材料具有多种同素异构体, 进而呈现出更丰富的性质. 层状Bi2Se3由于其独特的物理性质, 受到人们广泛的关注, 而它的同素异构体尚未有人研究. 本文采用基于密度泛函理论的结构搜索方法, 预测了一个稳定的β-Bi2Se3新相, 它具有良好的动力学和热力学稳定性, 并在低Bi2Se3源化学势条件下容易形成. 单层β-Bi2Se3是一个直接带隙为2.44 eV的二维半导体, 其电子载流子有效质量低至0.52m0, 在可见光范围内具有高达105 cm–1的光吸收系数, 并且能带边缘位置适中, 可用于光催化水分解制氢气. 此外, 由于β-Bi2Se3在垂直层面方向的镜面对称性破缺, 能够产生面外极化强度, 具有0.58 pm/V的面外压电系数. 鉴于其新颖的电子特性, 二维β-Bi2Se3在未来的电子器件中可能发挥重要的作用.

English Abstract

    • Bi2Se3是由V-VI主族元素组成的化合物半导体, 属于六方晶系, 具有辉铋碲矿结构, 如图1(a)所示, 我们将其命名为α-Bi2Se3. 该固体具有层状结构, 每层由5个原子层(quintuple layer, QL)组成, 以Se-Bi-Se-Bi-Se方式排布, 层厚约0.96 nm, 由共价键和离子键结合而成, 而QL-QL之间则是通过范德瓦耳斯力结合, 因此可以通过机械剥离得到QL[1]. α-Bi2Se3可通过化学气相沉积[2]、物理化学混合气相沉积[3]、金属-有机化学气相沉积[4,5]、脉冲激光沉积[6]和分子束外延[7-9]等方法制备得到. 2009年, Zhang等[1]通过理论计算预言了α-Bi2Se3是拓扑绝缘体. 同年, 实验上通过角分辨光电子能谱仪第一次观察到其表面的狄拉克锥能带结构, 证明了α-Bi2Se3是一种三维拓扑绝缘体材料[10]. 作为表面态只有1个狄拉克点同时具有较大带隙的强拓扑绝缘体[10], α-Bi2Se3为拓扑绝缘体高性能光电器件的研究提供了可能性[11]. α-Bi2Se3对红外光和太赫兹波有高响应, 可用于红外和太赫兹探测和成像[12]. 此外, α-Bi2Se3薄膜在热电应用方面极具潜力, 其费米能级附近超高的态密度使得Seebeck系数显著增大[13-15].

      图  1  (a) α-Bi2Se3的原子结构; (b)单层β-Bi2Se3结构的俯视图(上图)和侧视图(下图); (c)双层β-Bi2Se3结构的俯视图(上图)和侧视图(下图); (d)经过10 ps第一性原理分子动力学模拟, 得到了300 K时Bi2Se3单层的平衡结构; (e) β-Bi2Se3的声子谱; (f) β-Bi2Se3单层的电子局域函数

      Figure 1.  (a) Atomic structure of α-Bi2Se3; (b) the top and side views of monolayer β-Bi2Se3; (c) the top and side views of bilayer β-Bi2Se3; (d) snapshots of the equilibrium structures of the β-Bi2Se3 monolayer at 300 K after 10 ps ab initio molecular dynamic simulation; (e) phonon dispersion of monolayer β-Bi2Se3; (f) electron localization function for monolayer β-Bi2Se3.

      对于具有特定化学计量比的二维材料, 由于原子杂化方式的不同, 可能存在多种同素异构体, 而异构体的不同原子结构赋予了它们独特的物理和化学性质. 例如石墨烯和石墨炔是二维碳材料的同素异构体: 石墨烯具有sp2杂化的蜂巢结构, 是零带隙半导体, 拥有线性色散关系的狄拉克锥[16]; 石墨炔是由碳的sp和sp2轨道杂化而成, 它具有不同方向变形的狄拉克锥[17]. 黑磷和蓝磷是二维磷烯的同素异构体: 黑磷具有约2 eV的直接带隙和高载流子迁移率等优点[18,19], 可用于可见光范围的光电器件; 蓝磷具有2.73 eV的间接带隙, 在紫外区域吸收较强[20,21], 是紫外探测器件的备选材料. 因此, 探寻新型二维材料可能存在的同素异构体对于研究新奇物理效应、开发具有特定功能的材料体系具有重要推动作用.

      本研究首次提出一种新型Bi2Se3同素异形体—β-Bi2Se3. 在密度泛函理论框架下, 采用基于粒子群算法的结构预测方法, 搜索了化学计量比为2∶3的Bi2Se3二维结构, 并进一步计算其稳定性和物理性质. 结果表明, 我们预测的β-Bi2Se3具有良好的动力学稳定性和热稳定性, 它是能隙为2.40 eV的直接半导体, 载流子有效质量低至0.52m0, 在可见光范围内具有高达105 cm–1的光吸收系数, 能带边缘位置适中, 可用于光催化水分解制氢气. 此外, 由于镜面对称性破缺, β-Bi2Se3具有垂直平面的面外压电系数(0.58 pm/V). 因此, β-Bi2Se3在微电子、光电子和压电器件, 以及光催化等领域应用上具有很好的潜力.

    • 本文的第一性原理计算采用VASP软件包来实现[22], 它是基于密度泛函理论的软件. 采用PAW势来描述电子-离子相互作用[23]. 交换关联泛函采用广义梯度近似下的PBE泛函[24]. 平面波基组的动能截断取为500 eV, 能量收敛判据为10–5 eV, 力收敛判据为0.01 eV/Å. 晶胞垂直层面的真空层为20 Å. 采用HSE06杂化泛函计算电子性质[24,25], k点在布里渊区的密度设为0.015/Å. 采用DFT-D3方法描述多层体系中的长程范德瓦耳斯作用[26].

      为了研究β-Bi2Se3的动力学稳定性, 我们采用基于密度泛函微扰理论(DFPT)和VASP相结合的Phonopy程序计算声子谱[27,28]. 采用以PAW方法和PBE泛函为基础的从头计算分子动力学(AIMD)方法预测Bi2Se3的热稳定性[29]. 在AIMD模拟中, 以3 × 3 × 1的Bi2Se3超胞为初始结构, 模拟了温度为300 K的NVT系综, 总时长为10 ps, 步长为1.0 fs, 温度采用Nosé-Hoover方法控制[30].

      采用CALYPSO程序中的粒子群优化(PSO)算法搜索二维(Bi2Se3)n (n = 1—4)单层的结构[31,32], 种群大小和代数都设置为30. 结构搜索的计算使用VASP5.4程序, 基于PBE泛函进行结构弛豫[22].

    • 在应用CALYPSO进行全局的结构搜索之后, 生成了数百个结构, 并比较了它们的能量(见附录图A1), 其中β-Bi2Se3是众多结构中最稳定的. 图1(b)给出了β-Bi2Se3单层结构, 它具有六角晶体结构, 并且沿z方向产生了褶皱, 单胞包含2个Bi原子和2个Se原子, 每个Bi原子与3个Se原子成键, 每个Se原子与2个Bi原子成键. 晶格常数a为7.17 Å, 厚度h为1.40 Å, Bi–Se键长为2.69 Å. 为了深入分析β-Bi2Se3的化学键性质, 我们计算了电子局域函数(ELF), 如图1(f)所示, 电子密度主要位于Se原子周围, 而在Bi和Se之间也布满了电子波函数, Bi—Se键的键级为0.4, 均表明β-Bi2Se3主要由极性共价键构成. β-Bi2Se3具有良好的动力学稳定性和热稳定性, 如图1(d)所示, β-Bi2Se3单层的声子谱无虚频, 且在300 K下运行10 ps的AIMD模拟后, 仍可保持原有的层状结构, 体系能量在10 ps内变化也很小(见附录图A2).

      图  A1  CALYPSO搜索得到的几个较低能量的Bi2Se3单层结构(a)及对应的声子谱(b), 其中Bi2Se3-1, Bi2Se3-2, Bi2Se3-3的形成能分别为–0.15, –0.12, –0.09 eV/atom

      Figure A1.  Some typical low-energy structures (a) of freestanding Bi2Se3 monolayer predicted by the CALYPSO code and corresponding phonon dispersions (b). The formation energy of Bi2Se3-1, Bi2Se3-2, Bi2Se3-3 are –0.15, –0.12, –0.09 eV/atom respectively.

      图  A2  温度为300 K时β-Bi2Se3单层的能量-时间变化 曲线

      Figure A2.  Variations of temperature and energy with the time of AIMD simulation for β-Bi2Se3 monolayer at 300 K.

      为了评估β-Bi2Se3的能量稳定性, 我们首先计算了它的形成能, 定义为

      $ {H_{\rm{f}}} = {{ }}({E_{{\rm{tot}}}} - {n_1} \times {E_{{\rm{Bi}}}} \times {n_2} \times {E_{{\rm{Se}}}})/n, $

      其中Etotβ-Bi2Se3的总能, EBiESe为Bi和Se固体中每个原子的能量, 因子n1n2表示β-Bi2Se3单胞中Bi和Se的原子数, 因子n表示单胞内的总原子数. 计算得到Hf = –0.20 eV/atom, 表明由Bi和Se固体形成β-Bi2Se3是放热的, 并且这个值接近实验上已获得的α-Bi2Se3单层的形成能(图1(a), Hf = –0.39 eV/atom). 为了评估α-Bi2Se3β-Bi2Se3在实验制备过程中的相互竞争, 我们计算了表面自由能(ΔG)随Bi2Se3单层的化学势(${\mu _{{\rm{B}}{{\rm{i}}_{\rm{2}}}{\rm{S}}{{\rm{e}}_3}}}$)的变化曲线[33]. 表面自由能定义为

      $ \Delta G = (G - {G_{\rm{\alpha }}} - N \times {\mu _{{\rm{B}}{{\rm{i}}_{\rm{2}}}{\rm{S}}{{\rm{e}}_3}}})/A, $

      其中GGα分别是预测体系和α-Bi2Se3单层的吉布斯自由能, N是超元胞内Bi2Se3分子式的数量, A是单层的表面积. 在给定化学势条件下, 在一定的衬底上制备Bi2Se3单层, 具有表面自由能的系统是最稳定的, 因此可以在相图中出现. 如图2所示, 在缺Bi2Se3的实验条件下, β-Bi2Se3容易在衬底上形成; 而在富Bi2Se3条件下, α-Bi2Se3容易形成. 两个异构体的ΔG曲线的交点给出了两种相的临界转变点(${\mu _{{\rm{B}}{{\rm{i}}_{\rm{2}}}{\rm{S}}{{\rm{e}}_3}}}$≈ –0.49 eV), 在这个交点处, 两种Bi2Se3同素异构单层具有相同的出现概率. 这些结果表明我们所预测的β-Bi2Se3在实验上很有可能被制备.

      图  2  α-Bi2Se3β-Bi2Se3体系表面自由能的化学势相图

      Figure 2.  Chemical potential phase diagram of surface free ener-gy for α-Bi2Se3 and β-Bi2Se3.

      接着计算β-Bi2Se3单层的电子能带结构和局域态密度(LDOS), 如图3所示. 采用HSE06杂化泛函, β-Bi2Se3Γ点具有2.66 eV的直接带隙; 考虑自旋轨道耦合效应(SOC)后, 带隙为2.40 eV, 比未考虑SOC减小了约10%, 说明β-Bi2Se3具有较强的SOC. 而考虑SOC情况下单层α-Bi2Se3仅具有0.88 eV的间接带隙, 这个结果与前人的理论计算结果相符合(大约为1.0 eV)[34], 同时实验上所制备的多层α-Bi2Se3带隙随着厚度的增加逐渐减小(0—0.25 eV)[35], 由此推断单层α-Bi2Se3将具有更大的带隙, 与本文计算值定性相符, 说明了我们计算结果的可靠性. LDOS分析表明, β-Bi2Se3价带的边缘主要由Se原子的4p轨道占据, 部分来自Bi原子的6p轨道, 而导带边缘是由Bi原子的6p轨道和Se原子的4p轨道共同贡献. 同时, 在费米能级附近的LDOS的大面积重叠表明: Bi和Se原子的轨道之间具有较强的杂化作用.

      图  3  (a)不考虑SOC和(b)考虑SOC时, 采用HSE06泛函计算得到的β-Bi2Se3的能带结构和LDOS

      Figure 3.  The electronic band structures (left panel) and LDOS (right panel) (a) without and (b) with SOC effect for monolayer β-Bi2Se3 using HSE06 functional, respectively.

      由于实验上合成的二维材料薄片可能为多层, 我们还考虑了双层β-Bi2Se3的几种高对称堆叠构型(见附录图A3), 其中图1(c)所示的堆叠方式是最稳定的双层结构, 晶格常数为7.01 Å, 层间距为3.14 Å, 大于Bi—Se键长(2.69 Å), 层间结合能为–0.065 eV/atom, 与石墨烯(–0.065 eV/atom)[36]、磷烯(–0.055 eV/atom)[37]、h-BN(–0.065 eV/atom)[36]的值相近, 说明层间作用为范德瓦耳斯力, 多层或块体β-Bi2Se3不会在层间形成共价键发生重构. 双层β-Bi2Se3的带隙为0.77 eV(见图4). 我们还考察了β-Bi2Se3块体的性质, 其晶格常数为7.00 Å, 层间距为2.37 Å, 考虑SOC计算得到的带隙为0.87 eV, 并且保持了直接带隙的特点. 由此可见, 二维β-Bi2Se3的层厚对电子结构有较大影响: 双层(0.77 eV)和块体(0.87 eV)的带隙很接近, 并且远小于单层的带隙(2.40 eV), 但是体系的厚度并不会影响直接带隙的特点. 电子和空穴的有效质量可以通过拟合导带最小值(CBM)和价带最大值(VBM)附近的二次函数得到. 计算β-Bi2Se3的单层、双层、块体的载流子有效质量范围分别为0.66m0—7.88m0, 0.52m0—2.55m0, 0.63m0—0.67m0(见表1). 值得注意的是, 沿y方向的载流子有效质量可低至0.52m0, 说明载流子在β-Bi2Se3中可能具有较高的迁移率. 而且对于单层和双层β-Bi2Se3, 空穴和电子载流子的有效质量具有很大差异, 将有利于对空穴和电子的筛选, 在太阳能转换方面具有优势. 为了考察应力对能带的影响, 我们进一步计算了单层β-Bi2Se3在双轴应变下的能带结构(见图4(b)), 发现在–3%到3%的应变下, 单层β-Bi2Se3的带隙从2.61 eV线性降低到2.22 eV, 并且依然保持直接带隙的特征.

      β-Bi2Se3VBM/eVCBM/eVmxhmxemyhmye
      Monolayer–5.82–3.437.880.705.690.66
      Bilayer–5.00–4.242.550.522.360.52
      Bulk0.650.630.650.67

      表 1  单层、双层和块体β-Bi2Se3相对真空能级的价带顶VBM和导带底CBM, 空穴和电子沿着xy方向的有效质量(mxh, myh, mxe, mye). 载流子有效质量以自由电子的静止质量m0为单位

      Table 1.  The VBM and CBM related to vacuum level for monolayer, bilayer and bulk β-Bi2Se3, and the corresponding carrier effective mass. m0 is the electron rest mass.

      图  4  (a)采用HSE06泛函并且考虑SOC效应的双层(左图)和块体(右图)β-Bi2Se3的能带结构; (b)单层β-Bi2Se3带隙随双轴应变的变化

      Figure 4.  (a) The electronic band structures for bilayer (left panel) and bulk (right panel) β-Bi2Se3 based on HSE06 level with SOC effect; (d) effect of biaxial strain on band gap of monolayer β-Bi2Se3.

      图  A3  不同堆叠方式的双层β-Bi2Se3 (a)能量最低的β-Bi2Se3双层结构, 将它的能量设定为0 eV; (b)相对能量为0.32 eV;(c)相对能量为0.55 eV

      Figure A3.  β-Bi2Se3 bilayer with different stacking types and their relative energies: (a) the atomic structure of β-Bi2Se3 bilayer with the lowest energy, and its energy is set to 0 eV; the bilayer structures with relative energies of 0.32 eV (b) and 0.55 eV (c), respectively.

      接着, 我们探索了β-Bi2Se3用于光催化水分解的可能, 其中一个重要条件是催化剂的带边必须跨越水的氧化还原电势. 对于水分解反应, 标准还原电势为${E_{{{\rm{H}}^{{ + }}}{{/}}{{\rm{H}}_2}}}$ =–4.44 eV + pH × 0.059 eV, 标准氧化电势为${E_{{{\rm{O}}_2}{{/}}{{\rm{H}}_2}{\rm{O}}}}$ = –5.67 eV + pH × 0.059 eV, 二者均取决于pH值[38-41]. 计算表明, 单层β-Bi2Se3的VBM和CBM(相对真空能级)分别为–5.82 eV和–3.43 eV (见表1), pH = 0和7时带边均跨过水的氧化还原电势(图5(a)), 这说明在酸性和中性环境下, 单层β-Bi2Se3均适用于光催化水分解. 由于双层β-Bi2Se3的带隙很小, 不足以提供水分解的氧化还原势, 因此我们接下来只讨论单层β-Bi2Se3的光催化水分解性质. 良好的光吸收性质是光催化水分解的另一个必要条件. 我们采用HSE06泛函计算了β-Bi2Se3单层的复介电函数, 进而得到光吸收系数(定义为光强在单位长度介质中扩散的衰减). 如图5(b)所示, 单层β-Bi2Se3在可见光范围内表现出很强的光学吸收能力, 可见光谱吸收系数可达105 cm–1, 高于α-Bi2Se3的光吸收系数(104—105 cm–1)[42], 因此可高效捕获大部分太阳光能量用于驱动水分解制氢气.

      图  5  (a)单层β-Bi2Se3的VBM和CBM对比pH = 7和pH = 0的氧化还原电势; (b)单层β-Bi2Se3的光吸收系数, λ是波长, 虚线中间区域表示可见光区

      Figure 5.  (a) The location of VBM and CBM relative to vacuum energy of monolayer β-Bi2Se3 at pH = 0 and 7; (b) optical absorption coefficient for monolayer β-Bi2Se3. λ is the wave length, and the area between the red and the purple represents the visible range

      最后, 我们研究了单层β-Bi2Se3的压电性质. 由于单层β-Bi2Se3属于C3V对称点群, 在z方向不具有对称性, 因此β-Bi2Se3具有垂直层面的压电效应. 为了计算β-Bi2Se3的压电系数, 我们采用Duerloo等[43]发展的方法估算不同单轴应变下的弹性常数和极性的变化. 面内的弹性常数可以通过下列公式得到:

      $\begin{split} &{C_{11}} = \frac{1}{{{A_0}}}\frac{{{\partial ^2}U}}{{\partial \varepsilon _{11}^2}},\;\;{C_{22}} = \frac{1}{{{A_0}}}\frac{{{\partial ^2}U}}{{\partial \varepsilon _{22}^2}},\\ &{C_{12}} = \frac{1}{{{A_0}}}\frac{{{\partial ^2}U}}{{\partial {\varepsilon _{11}}\partial {\varepsilon _{22}}}}, \end{split}$

      其中Uβ-Bi2Se3体系的总能; ε11, ε22是平面内沿着两个方向的应变; A0是平衡态下单胞的面积. ε11, ε22取–0.006 到0.006, 步长为0.002, 在每一个应变下, β-Bi2Se3的原子位置被充分优化. 结果显示, β-Bi2Se3单层的弹性常数C11C12分别为11 N/m和6 N/m(对于六角晶格, C11 = C22). 基于现代极化理论[44], 通过估算单轴应变下的极化强度的变化(eijk), 我们计算了β-Bi2Se3单层的压电系数(dijk), 根据以下公式得到:

      $ {e}_{ijk}=\frac{\partial {{P}}_{i}}{\partial { \varepsilon }_{jk}}{,} $

      $ {d}_{ijk}=\frac{\partial {{P}}_{i}}{\partial { \sigma }_{jk}}{. } $

      其中εjk为应变张量, σjk是应力张量, Pi 为本征极化张量. 由于单层β-Bi2Se3的对称性特点, 其压电系数仅有面外系数e31d31, 这些压电系数与弹性力学常数有关:

      $ {{e}_{31}=d}_{31}({C}_{11}+{C}_{12}){. } $

      根据以上公式, 通过P3ε11的线性关系得到e31 = 0.10 × 10–10 C/m, d31 = 0.58 pm/V, β-Bi2Se3的压电系数接近甚至优于一些常见压电材料的数值, 如块体BN(d31 = 0.33 pm/V)和GaN(d31 = 0.96 pm/V)[45,46], 也媲美三族硫属化物单层Janus结构(Ga2SSe, Ga2STe, Ga2SeTe, In2SSe, In2STe, In2SeTe, GaInS2, GaInSe2, GaInTe2)的值[47]. 相比二维三族硫属化物需构建Janus结构以打破镜面对称性、引发面外极化, 二维β-Bi2Se3具有本征的面外压电效应, 因而将更具应用前景.

    • 利用第一性原理结构搜索技术, 我们预测了Bi2Se3的一个新相—β-Bi2Se3, 通过化学势相图计算、声子谱分析、分子动力学模拟, 该结构具有良好的热稳定性和动力学稳定性. β-Bi2Se3单层具有2.40 eV的直接带隙和低至0.52m0的载流子有效质量. 并且在可见光范围内, 表现出超高(105 cm–1)的光吸收系数, 能带边缘位置适中, 适用于光催化水分解制氢气. 此外, 由于β-Bi2Se3在垂直层面方向的非对称性, 引起了面外极化, 面外压电系数可达0.58 pm/V. 这些优异的物理性质使得层状β-Bi2Se3有希望成为未来高速微电子器件、光电器件、能量转换装置的组成材料.

      感谢内布拉斯加大学林肯分校化学系Zeng Xiao Cheng教授的讨论.

      附录

参考文献 (47)

目录

    /

    返回文章
    返回