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

许佳玲 贾利云† 刘超 吴佺 赵领军 马丽 侯登录

1) (河北建筑工程学院数理系, 张家口 075000)

2) (河北师范大学物理学院, 石家庄 050016)

1 引 言

拓扑绝缘体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如果用A, B, C 和D 表示4 个坐标位置, 可以把4 个面心立方子晶格看成4 个等效晶位, 用晶胞对角线连接起来, 半Heusler 结构C 位没有填充.由于目前的Heusler 化合物的组合方法一般是X 和Y 元素属于在化学元素周期表中的过渡金属元素,但是, 随着材料模拟领域的发展以及实验技术的进步, 其他元素如碱金属等也被用于Heusler 化合物的研究.于是这就大大的扩展了Heusler 化合物的种类, 也极大地增大了各种Heusler 化合物的研究潜能.2015 年林诗源等[16]采用赝势高通量方法计算了2295 种半Heusler 化合物, 并统计了 779 种拓扑绝缘体材料的能带反转强度和体结构带隙, 发现LiAuS 和 NaAuS 体系是非常优异的拓扑绝缘体材料, 本部分工作主要参考采用赝势的VASP软件高通量计算结果, 根据扩展的元素种类用全势方法研究α 型填充四面体结构的拓扑绝缘体材料,计算Li(Na)AuS 体系拓扑绝缘体材料的能带结构,由于Heusler 化合物研究很少有人利用S 放在主族元素所占的Z 位置上, 用S 作为主族元素, 采用贵金属Au 占据A 位, 为保证18 个价电子, X 元素用比较活泼的碱金属进行替代, 研究此类拓扑绝缘体的特点, 对拓扑绝缘体的理论和实验研究具有参考价值.

2 计算方法

本工作中能带结构的计算是利用全势线性缀加平面波方法的WIEN2k 来进行的.交换关联泛函主要是采用了传统的局域自旋密度近似(LSDA)和修正的Becke-Johnson 交换势 (MBJ).在计算中均对重金属元素考虑了自旋轨道耦合的效果.软件利用四面体网格法对于波矢进行积分, 选取5000 个k 点.截断参数选取Rmt× Kmax= 9.为了验证计算结果的可靠性, 采用LSDA 和MBJ 验算了拓扑绝缘体材料LiAuS 的能带结构, 对比文献中的采用广义梯度近似(GGA)和杂化泛函(HSE)的计算结果.

3 结果和讨论

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%处能量最低, 此时破坏了立方对称性之后成为四方结构后空间群为且平衡晶格常数为a = 0.52477 nm,c/a = 1.41 时系统又趋于平衡状态.

图1 对LiAuS 在保证体积不变的基础上对c 轴施加单轴应力得到的能量优化曲线, c/a = 1.41 时得到平衡晶格常数为a = 0.52477 nmFig.1.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%的单轴拉力得到的能带结构图, 可以发现较小的单轴应力破坏立方结构后也破坏了此类拓扑绝缘体的自然带隙, 通过计算得到施加较小的拉力和压力均会使系统由拓扑绝缘体转变为拓扑金属.

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

3.2 Li(Na)AuS 体系的化学替代

图2 施加1%拉应力后破坏立方结构的能带结构, 左图为LSDA 计算能带图, 右图为MBJ 计算得到的能带图Fig.2.Band structure of the LiAuS compound with 1% uniaxial tensile stress, on the left with LSDA, and on the right with MBJ.

图3 四方结构平衡晶格常数(a = 0.52477 nm, c/a = 1.41)下的能带结构 (a) LSDA 计算能带图; (b) MBJ 计算得到的能带图Fig.3.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.

图4 掺杂结构在各自平衡晶格常数下的能带结构 (a) Li0.125Na0.875AuS; (b) Li0.375Na0.625AuS; (c) Li0.5Na0.5AuS; (d) Li0.875Na0.125AuSFig.4.Band structure of the LixNa1 – xAuS compound with the equilibrium lattice constant: (a) Li0.125Na0.875AuS; (b) Li0.375 Na0.625AuS; (c) Li0.5Na0.5AuS; (d) Li0.875Na0.125AuS.

从前面的计算中可以看出, 等体积的单轴拉应力作用下破坏了立方对称性从而打开了拓扑绝缘体的能隙, 但是在实际材料的制备过程中, 等体积的单轴拉应力的实现存在困难.对于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变为对称性稍有下降; 而替换50%时生成结构的空间群为立方对称性遭到破坏的程度较为强烈, 不同比例下由于进行元素替代导致的对称性变化引起了能带图沿着高对称线能带的简并度的不同, 另外由于样品对应的平衡晶格常数不同, 也使能带图出现微小变化.本文得到的大量研究结果可以对半Heusler 合金结构拓扑绝缘体的进一步理论研究和实验研究提供依据.

4 结 论

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

感谢河北工业大学材料科学与工程学院张小明老师的指导和帮助.