删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Li(Na)AuS体系拓扑绝缘体材料的能带结构

本站小编 Free考研考试/2021-12-29

摘要:在拓扑领域中发现可以通过大数据搜索拓扑绝缘体, 使得此领域对材料的探索转变为对材料性质的研究. 半Heusler合金体系是非平庸拓扑绝缘体材料的重要载体. 通过全势线性缀加平面波方法计算Li(Na)AuS体系拓扑绝缘体材料的能带结构. 采用各种关联泛函计算LiAuS的平衡晶格常数, 发现得到的能带图均为具有反带结构的拓扑绝缘体, 而且打开了自然带隙. 较小的单轴应力破坏立方结构后也破坏了此类拓扑绝缘体的自然带隙, 通过施加单轴拉应力直到四方结构的平衡位置时, 系统带隙值约为0.2 eV, 这与立方结构平衡位置得到的带隙结果一致. 运用同族元素替代的手段, 实现了在保证材料拓扑绝缘体性质的同时, 不改变立方结构, 在体系的平衡晶格常数下使得材料的带隙打开, 从而提高了实验合成拓扑绝缘体材料的可行性.
关键词: 拓扑绝缘体/
Heusler合金/
第一性原理/
能带结构

English Abstract


--> --> -->
拓扑绝缘体HgTe的发现引发了世界的轰动, 尤其是最近我国在拓扑绝缘体领域取得了重大发现, 预测了上千种可能的新型拓扑绝缘体材料, 这就使得对这一领域的研究方向由材料的探索转变为材料性质的研究[1,2]. 半Heusler化合物与闪锌矿结构的相似性注定了其非凡的使命, 大量的研究致力于在Heusler材料家族中寻找非平庸拓扑绝缘体材料[3-5]. 具有18个电子的三元半Heusler合金可以通过合适的应力作用打开带隙, 由此开启了一个三维拓扑绝缘体的新系列[6-8]. 由于第一原理计算已证实此结构体系拓扑Z2不变量没有反对称性的量子本质, 加上这个系列材料如磁性和超导性等其他有趣的物理性质, 于是基于密度泛函理论的第一性原理计算方法将三元半Heusler化合物转变为三维拓扑绝缘体成为研究热点[9,10]. 对于三维拓扑绝缘体半Heusler化合物的研究, 结合这些材料的磁性和超导性质, 在实验方面也成为了一个有巨大价值的研究平台[11,12].
目前通过第一原理计算, 大量的半Heusler三元化合物已经被发现具有反带结构, 例如LuPtSb, ScPtBi, YPdBi, ThPtPb, 都可以通过单轴应力成为三维拓扑绝缘体的备选材料[13]. 通常, 对于具有小带隙的半Heusler化合物, 可以通过施加静压应力改变带序, 以及利用单轴应力打开能隙的综合方法实现拓扑绝缘体相[3]. 姚裕贵研究组[10]发现: LaPtBi在单轴应力下可以实现拓扑绝缘体相, 未形变的的LaPtBi的电子结构作为原型系统显示出明显的能带翻转特征. 通过施加[001]方向的单轴应变可以实现三维拓扑绝缘体, 在保留翻转带序的同时可以打开带隙, 并且通过计算不具有反对称性的系统中的拓扑不变量Z2证明LaPtBi是三维拓扑绝缘体. 张小明等[14]又模拟了LaPtBi在等体积下施加不同单轴拉压应力来破坏立方对称性, 并实现了能带的调控, 还对未形变体积下的单轴应力也进行了详细的研究. 关于半Heusler型化合物的能带翻转机制在参考文献[15]中有较为详细的论述, 而且给出了一些设计的规则, 另外文中将HgTe一类反转带序称为β型填充四面体结构的拓扑半金属, 这一类结构已经有人进行了大量的理论和实验研究; 还有具有另一类反转带序的Heusler体系拓扑绝缘体被称为α型填充四面体结构, 此类结构的材料预测目前研究较少, 实验研究几乎还没有开展.
Heusler化合物研究始于由4个面心立方子晶格组成的L21结构, 基于相似的面心晶格, 4个Wyckoff坐标分别为4a (0, 0, 0), 4b (1/2, 1/2, 1/2), 4c (1/4, 1/4, 1/4) 和4d (3/4, 3/4, 3/4). 4个子晶格之中有一个没有被填充, 就是C1b结构, 这种结构通常也称为半Heusler结构. 半Heusler化合物的化学式一般写为XYZ, 空间群为216 ($ F\bar43m $). 如果用A, B, CD表示4个坐标位置, 可以把4个面心立方子晶格看成4个等效晶位, 用晶胞对角线连接起来, 半Heusler结构C位没有填充. 由于目前的Heusler化合物的组合方法一般是XY元素属于在化学元素周期表中的过渡金属元素, 但是, 随着材料模拟领域的发展以及实验技术的进步, 其他元素如碱金属等也被用于Heusler化合物的研究. 于是这就大大的扩展了Heusler化合物的种类, 也极大地增大了各种Heusler化合物的研究潜能. 2015年林诗源等[16]采用赝势高通量方法计算了2295种半Heusler化合物, 并统计了 779 种拓扑绝缘体材料的能带反转强度和体结构带隙, 发现LiAuS 和 NaAuS 体系是非常优异的拓扑绝缘体材料, 本部分工作主要参考采用赝势的VASP软件高通量计算结果, 根据扩展的元素种类用全势方法研究α型填充四面体结构的拓扑绝缘体材料, 计算Li(Na)AuS体系拓扑绝缘体材料的能带结构, 由于Heusler化合物研究很少有人利用S放在主族元素所占的Z位置上, 用S作为主族元素, 采用贵金属Au占据A位, 为保证18个价电子, X元素用比较活泼的碱金属进行替代, 研究此类拓扑绝缘体的特点, 对拓扑绝缘体的理论和实验研究具有参考价值.
本工作中能带结构的计算是利用全势线性缀加平面波方法的WIEN2k来进行的. 交换关联泛函主要是采用了传统的局域自旋密度近似(LSDA)和修正的Becke-Johnson 交换势 (MBJ). 在计算中均对重金属元素考虑了自旋轨道耦合的效果. 软件利用四面体网格法对于波矢进行积分, 选取5000个k 点. 截断参数选取Rmt × Kmax = 9. 为了验证计算结果的可靠性, 采用LSDA和MBJ验算了拓扑绝缘体材料LiAuS的能带结构, 对比文献中的采用广义梯度近似(GGA)和杂化泛函(HSE)的计算结果.
2
3.1.不同关联泛函LiAuS能带结构
-->通过计算发现LiAuS和NaAuS具有反带结构, 这个结论和文献[15]中采用GW近似方法得到的结果一致; 另外, 对于反带结构和反转强度的讨论, 文献[16]中进行了细致的讨论, 不再作为本文讨论的重点. 利用GGA, LSDA和MBJ三种关联泛函计算了LiAuS. 各个化合物的平衡晶格常数都是通过能量优化得到的. 通过比较能量优化得到的晶格常数, 发现采用全势GGA泛函得到的平衡晶格常数为0.602 nm, 稍大于VASP赝势软件GGA计算值的0.599 nm, 而采用LSDA泛函优化得到的平衡晶格常数为0.586 nm, 可见局域密度泛函得到的平衡晶格常数为最小. 通过各自平衡晶格常数计算得到的能带图结构均为具有反带结构的拓扑绝缘体结构, 且均打开了自然带隙, 结果和文献一致, 但是在带隙数值上MBJ并没有显现出优势. 对于LaPtBi材料一般在平衡状态为反带结构, 但是并没有带隙, 普遍认为具有此类反带结构的立方半Heusler化合物要想打开体能隙, 常用的方法就是施加等体积单轴应力, 施加拉应力可以增大带隙, 从而得到真正意义上的拓扑绝缘体材料, 这种方法的实质是通过破坏立方对称性的保护而打开带隙[17,18]. 图1表示对LiAuS在保证体积不变的基础上对c轴施加单轴应力得到的能量优化曲线, 由于施加拉应力可以增大带隙, 所以应力施加是以拉应力为主, 图中横坐标的–10表示施加10%的单轴压应力, 而横坐标的正值均表示拉应力. 由图中可知, 当单轴拉力施加到约41%处能量最低, 此时破坏了立方对称性之后成为四方结构后空间群为$119 (I\bar4m2)$, 且平衡晶格常数为a = 0.52477 nm, c/a = 1.41时系统又趋于平衡状态.
图 1 对LiAuS在保证体积不变的基础上对c轴施加单轴应力得到的能量优化曲线, c/a = 1.41时得到平衡晶格常数为a = 0.52477 nm
Figure1. Calculated total energies as functions of the uniaxial strain along [001] direction with constant volume for LiAuS, the equilibrium lattice constant is a = 0.52477 nm and c/a = 1.41.

研究认为有必要对碱金属系列半Heusler化合物的拓扑绝缘体进行计算和分析, 另外由于全势的WIEN2k软件包的MBJ是原生的, 因此采用LSDA和MBJ两种关联泛函计算此系列的能带进行比较. 图2给出了利用WIEN2k软件分别采用LSDA和MBJ两种关联泛函计算得到的对LiAuS施加1%的单轴拉力得到的能带结构图, 可以发现较小的单轴应力破坏立方结构后也破坏了此类拓扑绝缘体的自然带隙, 通过计算得到施加较小的拉力和压力均会使系统由拓扑绝缘体转变为拓扑金属.
图 2 施加1%拉应力后破坏立方结构的能带结构, 左图为LSDA 计算能带图, 右图为MBJ计算得到的能带图
Figure2. Band structure of the LiAuS compound with 1% uniaxial tensile stress, on the left with LSDA, and on the right with MBJ.

图3给出了单轴应力下平衡状态下采用LSDA和MBJ两种关联泛函计算得到的LiAuS的能带结构图. 从图3可以看出, 通过施加单轴拉应力直到四方结构的平衡位置时, 系统带隙的值为0.2 eV左右, 这与立方结构平衡位置得到的带隙结果一致. 一般情况下, 计算半导体和绝缘体材料, 利用GGA方法计算得到的能带带隙要比MBJ的要小, 也就是通常所说的GGA会低估带隙, 发现在此系列的拓扑绝缘体带隙计算中, 每种关联泛函得到的带隙是一致的, 所以证实这是材料本身的属性. 但是, 在结果中通过MBJ泛函计算得到的能带图中, 费米能级明显穿过了价带. 通过计算发现, 在此体系的半Heusler化合物中, MBJ泛函在拓扑绝缘体的计算中并没有明显优势. 而在LSDA计算得到的能带图中, 在拉力作用下, Au的部分s轨道明显能量下降, 在局域密度近似当中尤其明显, 拉力为84%严重破坏立方对称性时, 带隙值仍然约为0.2 eV.
图 3 四方结构平衡晶格常数(a = 0.52477 nm, c/a = 1.41)下的能带结构 (a) LSDA 计算能带图; (b) MBJ计算得到的能带图
Figure3. Band structure of the tetragonal structure LiAuS compound with the equilibrium lattice constant (a = 0.52477 nm and c/a = 1.41): (a) LSDA; (b) MBJ.

2
3.2.Li(Na)AuS体系的化学替代
-->从前面的计算中可以看出, 等体积的单轴拉应力作用下破坏了立方对称性从而打开了拓扑绝缘体的能隙, 但是在实际材料的制备过程中, 等体积的单轴拉应力的实现存在困难. 对于Heusler化合物体系, 在实验上通过对某个原子同族的元素替换进行掺杂是比较常用的手段, 在保证价电子浓度不变的前提下, 通过掺杂可以改变晶格常数、电负性, 甚至实现半金属性, 还可以找到具有相变的新材料[19]. 实验上采用化学掺杂在Heusler化合物中获得成功的例子数不胜数, 由于化学掺杂也是对体系立方对称性的破坏, 因此运用同族元素替代的手段, 期望在保证材料拓扑绝缘体性质的同时, 不改变立方结构, 在体系的平衡晶格常数下使得材料的带隙打开, 从而提高实验合成拓扑绝缘体材料的可行性[20-22].
选取上文论述中两种本身都具有反带结构的材料, 观察Li和Na相互取代过程中是否能够打开带隙, 依次采用LSDA计算了相互替代比例为12.5%, 25.0%, 37.5%, 50.0%, 62.5%, 75.0%和87.5%的超胞结构的能带. 首先进行扩胞形成2 × 2 × 1超胞结构, 原子替换依据原有的套嵌结构由外侧立方顶角位置开始依次更换, 对于不同替代比例相应的超胞结构首先进行能量优化获得平衡晶格常数, 然后利用各自的平衡体积进行能带计算, 并且考虑自旋轨道耦合, 其他计算参数与前面的计算一致. 相互替代时没有考虑无序状态, 只研究由于掺杂引起体系立方对称性的破坏对能隙稳定性的影响. 如图4所示, 结果发现Li0.125Na0.875AuS, Li0.375Na0.625AuS, Li0.5Na0.5AuS, Li0.625Na0.375AuS, Li0.875Na0.125AuS都可以在平衡位置打开能隙, 而且在费米面位置处导带底和价带顶的能带十分稳定, 均为直接带隙, 带隙值约为0.2 eV左右. 进行替换后, 化合物材料的空间群由原来的216($ F\bar43m $)变为111($ P\bar42m $), 对称性稍有下降; 而替换50%时生成结构的空间群为115($ P\bar4m2 $), 立方对称性遭到破坏的程度较为强烈, 不同比例下由于进行元素替代导致的对称性变化引起了能带图沿着高对称线能带的简并度的不同, 另外由于样品对应的平衡晶格常数不同, 也使能带图出现微小变化. 本文得到的大量研究结果可以对半Heusler合金结构拓扑绝缘体的进一步理论研究和实验研究提供依据.
图 4 掺杂结构在各自平衡晶格常数下的能带结构 (a) Li0.125Na0.875AuS; (b) Li0.375Na0.625AuS; (c) Li0.5Na0.5AuS; (d) Li0.875Na0.125AuS
Figure4. Band structure of the LixNa1–xAuS compound with the equilibrium lattice constant: (a) Li0.125Na0.875AuS; (b) Li0.375Na0.625AuS; (c) Li0.5Na0.5AuS; (d) Li0.875Na0.125AuS.

通过全势线性缀加平面波方法计算Li(Na)AuS体系拓扑绝缘体材料的能带结构, 通过比较能量优化得到的晶格常数, 发现采用全势GGA泛函得到的平衡晶格常数大于赝势计算值, 局域自旋密度LSDA泛函得到的平衡晶格常数为最小. 通过各自平衡晶格常数计算得到的能带图结构均为具有反带结构的拓扑绝缘体结构, 而且均打开了自然带隙, 但是在带隙数值上MBJ并没有显现出优势. 较小的单轴应力破坏立方结构后也破坏了此类拓扑绝缘体的自然带隙, 通过施加单轴拉应力直到四方结构的平衡位置时, 系统带隙值约为0.2 eV, 这与立方结构平衡位置得到的带隙结果一致. 由于化学掺杂也是对体系立方对称性的破坏, 运用同族元素替代的手段, 实现了在保证材料拓扑绝缘体性质的同时, 不改变立方结构, 在体系的平衡晶格常数下使得材料的带隙打开.
感谢河北工业大学材料科学与工程学院张小明老师的指导和帮助.
相关话题/结构 计算 材料 系统 优化

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 吸收波长和吸收效率可控的超材料吸收器
    摘要:为了使超材料完美吸收器(metamaterialperfectabsorber,MPA)能够同时实现吸收效率和吸收波长的控制,本文提出利用二氧化钒(VO2)和石墨烯作为MPA的材料,通过对MPA的结构设计,在红外波段实现了高吸收,吸收效率最高可达99%.研究发现通过改变VO2的温度和石墨烯的化 ...
    本站小编 Free考研考试 2021-12-29
  • 基于二维六方氮化硼材料的光子晶体非对称传输异质结构设计
    摘要:二维六方氮化硼(hexagonalboronnitride,hBN)材料在产生光学稳定的超亮量子单光子光源领域有着潜在应用,有望用于量子计算和信息处理平台,已成为研究热点.而光学非对称传输设备是集成量子计算芯片中的关键器件之一.本文从理论上提出了一种基于hBN材料光子晶体异质结构的纳米光子学非 ...
    本站小编 Free考研考试 2021-12-29
  • 高温超导体电子结构和超导机理的角分辨光电子能谱研究
    摘要:超导是一种奇异的宏观量子现象.100多年来,已发现的超导体主要分为两类:以金属或者合金为代表的常规超导体以及以铜氧化物和铁基高温超导体为代表的非常规超导体.常规超导体的超导机理能被BCS超导理论完美解释,但高温超导体的超导机理至今仍未达成共识,已经成为凝聚态物理领域中长期争论且充满挑战的重大科 ...
    本站小编 Free考研考试 2021-12-29
  • 重费米子超导理论和材料研究进展
    摘要:重费米子超导体是一类典型的强关联和非常规超导系统,超导的产生与量子临界涨落有着紧密的关系.在实际材料中,不同结构体系的重费米子超导体往往表现出非常不同的竞争序和超导性质,表明f电子的行为对材料的结构特征具有敏感依赖性.特别是最近几年的超导实验研究,表明具体材料的实际电子结构对重费米子超导的性质 ...
    本站小编 Free考研考试 2021-12-29
  • 非对称自旋-轨道耦合系统的多体量子相干含时演化
    摘要:本文以具有非对称性自旋相互作用的三体自旋系统为研究对象,重点研究了三体量子相干含时演化规律.采用精确量子对角化和基于量子主方程的数值模拟方法,讨论了三体量子系统中多种量子相干组分及其退相干.研究发现,量子相干组分的含时演化与整个系统的初态量子特性紧密相关.当初态为可分离纯态时,在较短时间内,非 ...
    本站小编 Free考研考试 2021-12-29
  • 超级质子-质子对撞机中束流热屏的热-结构耦合模拟分析
    摘要:束流热屏(beamscreen)是新一代高能粒子对撞机中的重要部件,用于将束流在管道中运行时产生的热量转移到冷却系统中,同时通过束流热屏上的排气孔将残余气体输送至冷管壁上,维持良好的真空度.然而,在转移热负载的过程中,温度变化产生的形变会影响束流热屏的结构稳定性.如何在保证束流热屏良好传热性能 ...
    本站小编 Free考研考试 2021-12-29
  • 基于夹层结构的偏振无关1×2定向耦合型解复用器的设计
    摘要:提出一种基于夹层结构的偏振无关1×2定向耦合型解复用器,用于分离1310nm和1550nm两个波长.通过合理选择夹层结构芯区的折射率及波导间隙,可以调节同一波长两个正交偏振模的耦合长度相等,实现偏振无关;通过合理选择夹层结构波导宽度,可以使两个波长分别从不同输出波导端口输出,实现解复用功能.运 ...
    本站小编 Free考研考试 2021-12-29
  • 具有非广延分布电子的碰撞等离子体磁鞘的结构
    摘要:很多关于等离子体鞘层的研究工作都是基于电子满足经典的麦克斯韦速度分布函数,而等离子体中的粒子具有长程电磁相互作用,使用Tsallis提出的非广延分布来描述电子更为恰当.本文建立一个具有非广延分布电子的碰撞等离子体磁鞘模型,理论推导出受非广延参数q影响的玻姆判据,离子马赫数的下限数值会随着参数q ...
    本站小编 Free考研考试 2021-12-29
  • 基于2-2型压电复合材料的新型宽频带径向振动超声换能器
    摘要:本文提出了一种基于2-2型压电复合材料的新型宽频带径向振动超声换能器,它主要由内金属圆环和外压电陶瓷复合材料圆环组成.首先利用Newnham串并联理论和均匀场理论推导了2-2型压电复合材料的等效参数;其次利用解析法得到了金属圆环和径向极化压电复合陶瓷圆环径向振动的机电等效电路;最后得到了换能器 ...
    本站小编 Free考研考试 2021-12-29
  • 基于深度学习压缩感知与复合混沌系统的通用图像加密算法
    摘要:提出一种适用于灰度图像与RGB格式彩色图像的通用图像加密算法.利用双线性插值Bilinear与卷积神经网络对图像进行压缩,再使用二维云模型与Logistic组成的复合混沌系统对压缩图像加解密(滑动置乱与矢量分解),最后对解密图像进行重构.重构网络中,由卷积神经网络与双线性插值Bilinear主 ...
    本站小编 Free考研考试 2021-12-29