基于对称性指标预测拓扑材料 | 众妙之门

我们世界中的材料千变万化，丰富多彩。要想系统地理解这千千万万各种各样的材料，我们需要一些指导性原理、来开发描述这些材料的理论。以前的一个重要指导原理、就是朗道的对称性破缺理论。他指出材料的很多重要性质、取决于材料中的对称性。近二三十年来对量子材料的研究取得了很多进展。我们发现朗道的对称性破缺理论完全不够用。因为量子材料中有很多量子纠缠，而材料的量子性质，主要起源于这些量子纠缠。这样我们就需要一套全新的理论，如拓扑序理论和张量范畴理论，来描写这些量子纠缠的构型。这些理论都很有挑战性，而且起点也很高（见头条叶鹏的文章）。

可自然界中有很大一类材料，其中电子的相互作用并不重要，甚至可以忽略。这些材料中量子纠缠的构型，可以通过无相互作用电子能带，特别是能带的各种各样的拓扑性质来决定。我们知道能带的结构，是由电子在材料中所满足的各种各样的波动方程来决定的。对能带拓扑结构的研究，实际上就是对电子波动方程的拓扑结构的研究。相对于强关联电子的问题，这是一个简单得多的问题。最近几年，物理学家在这方面取得了很大的进展，使我们几乎可以完全系统的理解和计算、各种无相互作用电子材料的各种各样的拓扑性质。这些进展使我们能够发现和预言成千上万个有不同性质的拓扑材料，并展现其丰富多彩的拓扑性质。这篇文章，成文于这些进展的推动者，其详细系统地介绍了这方面最新的进展和成果。对于想了解和进入这一研究领域的学生和研究人员应该会很有帮助。

——文小刚

作者 | 唐峰1, 2 万贤纲1, 2

（1 南京大学物理学院固体微结构物理国家重点实验室）

（2 南京大学人工微结构科学与技术协同创新中心）

来源：中国物理学会期刊网

1 引 言

拓扑原本是一个数学概念，表示一类由系统整体决定而随细节不敏感的性质。这个概念融入到凝聚态物理归因于量子霍尔效应的发现。1980年Klitzing等人发现了二维电子气在强磁场中霍尔电导展现出不依赖于样品细节的量子化平台[1]。人们发现其中的相变不涉及任何对称性的破缺，脱离了朗道相变理论的范畴。1982 年Thouless 等人推导了二维电子气中的霍尔电导公式，并且将它与拓扑不变量——第一类陈数(∈ Z )联系起来[2]，成功地解释了量子霍尔电导的鲁棒性。这一工作开创性地将数学中的拓扑概念应用到固体晶格动量空间中波函数性质的研究，是数学与物理相结合的一大典范[2]。自此凝聚态中的拓扑性质受到广泛关注。6 年后，Haldane 发现在六角晶格中加入交错磁通也可实现量子霍尔效应，表明在凝聚态晶格体系中存在无朗道能级的量子霍尔效应[3]。但是实现量子霍尔效应需要施加巨大的外磁场[1]，这大大地限制了它在实际器件中的应用。这个方向近年来的一个重大突破是“量子自旋霍尔绝缘体(即二维拓扑绝缘体)”的发现[4， 5]。这一里程碑性的工作由美国宾夕法尼亚大学的Kane 和Mele[6，7]以及美国斯坦福大学的Bernevig，Hughes和Zhang[8]分别在石墨烯以及HgTe/CdTe 量子阱中完成，推动了拓扑量子物理成为凝聚态物理的前沿和主流。与量子霍尔效应不同，量子自旋霍尔效应无需外磁场，受到时间反演对称的保护，对应的拓扑不变量是Z2拓扑不变量 (只能取0或1) [7]。此后人们进一步将二维拓扑绝缘体推广到三维，得到了三维拓扑绝缘体[9—11]。三维拓扑绝缘体的拓扑不变量是4 个

(ν0；ν1，ν2，ν3) ，其中第一个不变量表征材料是否为强拓扑绝缘体，后面三个表征材料是否为弱拓扑绝缘体[4，5]。强拓扑绝缘体在任何面上均具有奇数个Dirac 点，表面态不受无序影响，而弱拓扑绝缘体只有在特定表面才具有拓扑表面态[4，5]。

进一步的研究表明，对称性以及空间维度对物理体系的拓扑性质有重要的影响。通过分析是否具有时间反演对称、粒子空穴对称以及手性对称性，人们对不同维度 (1D，2D，3D) 体系可能存在的拓扑相进行了系统的分类[12—15]。值得提及的是，除了上述的内禀对称性外，凝聚态体系中还存在大量的晶格对称性。2011 年，傅亮发现晶格对称对拓扑态有重要影响，从而开辟了拓扑晶体绝缘体这一新的研究方向[16]。目前包括镜面陈绝缘体[17]、Hourglass 拓扑晶体绝缘体[18]、Wallpaper拓扑晶体绝缘体[19]、高阶拓扑绝缘体[20—24]等大量拓扑晶体绝缘相相继被提出来[25， 26]。在上述各类拓扑的绝缘体中，体相是绝缘体，但是在材料的表面或者边界会出现拓扑保护的无能隙电子态[4，5，26]。

上述拓扑绝缘相的一个重要特性是，保持对称性的微扰对其拓扑性质不会有影响，具体地，只要这个微扰不关闭绝缘体的能隙，绝缘体的拓扑性质就不会发生改变。因此人们普遍认为拓扑物理只存在于绝缘体中。2011 年，通过研究烧绿石结构5d 过渡金属Ir 氧化物，人们发现随着电子关联强度U在合理范围内的变化，材料占据和非占据电子态的宇称将发生翻转[27]，从而发现体系有Weyl 点（即两个非简并线性色散能带的交点）。这时其低能电子行为可以用Weyl 方程描述，因而这类材料可命名为Weyl 半金属。Weyl 点是Berry曲率的源点，围绕一个Weyl点的闭合曲面上Berry曲率的积分绝对值为2π，该积分符号为Weyl 点所对应的“手性”符号[27]。由于Berry 曲率沿整个布里渊区边界面的面积分为零，所以对应的总的“手性”必须为零，在Weyl 半金属中，Weyl 点成对出现，外界微扰只能使Weyl 点移动，只有当两个手性相反的Weyl 点重合，Weyl 点才会消失。也就是说Weyl 点受拓扑保护稳定，和体系的对称性无关[27]。此外Weyl 半金属有着奇特的表面态Fermi 弧：Weyl 半金属表面态的费米面不是闭合的，而是一个开放的线段，这个开放的线段连接两个手性相反的Weyl 点在这个表面上的投影点[27]。Weyl 半金属只需要体系具有平移不变性，不一定需要其他的对称性。事实上它必须要求PT对称的破缺（P：空间反演，T：时间反演）。

除了Weyl 半金属外，大量其他类型的拓扑半金属相也被发现。通过分析布里渊区(BZ)中高对称点[28]或者高对称线上[29，30]的四度简并的能带交点，人们提出了低能激发满足Dirac 方程的Dirac半金属。考虑到能带接触点的类型，人们也提出了nodal line 半金属[31， 32]。进一步考虑若干nodalline 在BZ的位置以及它们之间的相对关系，人们又提出了nodal chain[33]、nodal link 和Hopf link 半金属[34—36]。在非简单群对称体系中，非简单操作要求能带形成Hourglass 型连接，人们据此提出了体Hourglass 半金属[37]。考虑到晶格对称保护的多重 (3，6，8 重) 简并点，人们也提出了多重简并费米子[38]。值得提及的是，除了Weyl 半金属外，上述这些拓扑 (半) 金属均形成需要晶格对称的保护才能稳定存在。

和普通材料相比，拓扑材料具有大量独特的物性。拓扑绝缘体特有的受拓扑保护的表面电子态具有自旋—动量绑定，背散射被禁闭[4，5，39]，因此拓扑绝缘体能够实现无能耗电子输运；拓扑绝缘体和超导的界面可以实现Majorana 费米子，这有望应用到量子计算中[40]；拓扑绝缘体表面Dirac电子态也具有反常的磁电耦合[41]等等，这一系列新奇的性质，具有广阔的应用前景。拓扑晶体绝缘体的无能隙表面态只在具有特定对称性的表面出现。由于晶体对称性更容易被外电场、应力所改变，拓扑晶体绝缘体很容易应用于场效应器件和压力感应器件[26]。Weyl 半金属也具有大量独有的性质：其具有和手性反常相关的高度各向异性磁电阻[42]；它也具有新颖的量子振荡现象[43]；通过上下两个表面Fermi 弧在Weyl 点投影处的隧穿行为，人们也理论设计[44]并且实验实现了三维量子霍尔效应[45]；其霍尔电导只和Weyl 点的符号和位置有关[46]；研究发现对Weyl 半金属掺杂可能实现拓扑超导体[47]，另外其光学性质也是人们关注的目标[48]。在拓扑半金属中实现Klein 隧穿型的输运现象[49]也引来了人们的广泛关注。Nodal line 半金属的拓扑表面态在费米面上有很大的态密度，所以它可能表现较高临界温度的超导行为[50]。这些新颖的物性可以预期在下一代电子、信息及至能源领域有很大应用潜力。事实上，这些对材料的缺陷、杂质细节不敏感的拓扑量子材料越来越成为物理基础及工程应用领域的前沿课题与追求。过去十年拓扑量子领域在世界上已经取得了巨大的研究进展，理论预言了大量拓扑量子材料[27—38，51—68]，很多预言以及一些关键性的拓扑量子效应已被实验证实[69—80]。

与常规的材料体系不同，拓扑体系的拓扑行为和能带细节无关，可用不同的拓扑不变量来描述。因此确定物理体系的拓扑不变量是判定其拓扑性质的核心。但是，实验判定拓扑不变量在技术上面临很大挑战，所以这一领域的一大特点就是通过理论计算拓扑不变量来预言新型拓扑材料，然后角分辨光发射谱、输运等实验证实[4，5，26，66，79，80]。这昭示出拓扑不变量的计算在找寻新型拓扑材料中起着至关重要的作用。但是，通过直接计算拓扑不变量来判断一个体系拓扑性质的方法效率较低、有很大的局限性：(1)计算拓扑不变量通常需要进行基于波函数的积分，这一计算工作量极大[81]；(2)对于一个材料，需要对所有的拓扑不变量进行计算，才能确定拓扑性质。这一问题也带来了很大的工作量。所以和人们已经成功生长出来的材料数目相比较，已知的拓扑材料数量还较为有限。此外，虽然人们已经发现了许多拓扑材料[4，5，26，66，79，80]，但是这些材料都不适合实际应用。比如拓扑绝缘体Bi2Se3 体系中容易形成缺陷，导致费米能级处于体能带的价带或者导带内[82]，拓扑晶体绝缘体SnTe 是p 型半导体，呈现了较好的体的导电性[83]。因此发展出新的理论方法，从而找到理想的、有实用价值的拓扑材料体系有着重要的科学价值和广阔的应用前景。

为了克服上述理论方法的困难，人们做了大量的工作。一个重要的进展是，傅亮等人发现在有中心反演对称的材料中，我们只需要计算空间反演不变点处价带的宇称，就可以得到Z2拓扑不变量，从而快速地判定体系是否为拓扑绝缘体[84]。除了中心反演对称性外，凝聚态体系还有很多其他空间对称性操作。基于材料电子能带在这些空间对称性操作的本征值，方辰等人发展了计算若干拓扑不变量的方案[85，93]。这些理论方法不需要计算耗时的波函数积分，可大大加快体系拓扑性质的判定过程。

另外，2017 年基于材料电子能带价带的实空间和倒空间图像的比较，人们发展了一套全面的拓扑分类方案：普林斯顿大学Benervig 研究组基于图论方法研究了具有实空间图像的能带在其布里渊区高对称点小群不可约表示之间的连接方式，提供了10403 个这样的基本连接方式，称之为“基本能带表示”[86]。他们进一步提出可以通过比较材料的价带和10403 个“基本能带表示”来判定拓扑性质。同时哈佛大学Vishwanath 研究组提出了能带指标理论，该理论将绝缘体能带用一组整数表示，并且把其拓展成一个群。通过这个群和它具有实空间图像的子群求商，从而对所有的空间群可能的拓扑行为进行了分类[87]。该理论方法给出了描述能带结构的商群分类，但是不能直接应用到拓扑材料搜索中[86]。Benervig 等人的理论方案确实可以用于筛选拓扑材料，但是“基本能带表示”的过多数目一定程度上限制了大规模应用。下面我们先介绍对称性指标理论，然后给出基于这一理论，我们发展出来从原子绝缘体基组出发实现高效判断材料拓扑的新方法[88]以及找到的典型拓扑材料。

2 对称性指标理论

凝聚态体系的电子特性可以用其电子的能带结构（能带色散和能带波函数）来描述，而电子的拓扑性质源于布洛赫本征态的特性。根据空间群对称，人们能够得到布洛赫本征态的变换性质，即本征态从属的不可约表示。在凝聚态体系布里渊区（BZ）里面有高对称点、高对称线、高对称面以及没有任何对称操作的普遍点。高对称点上能带可以用其对应小群的不可约表示来标记（简称为高对称点不可约表示），高对称点上能带不可约表示确定后，根据相容性关系能带在高对称线以及高对称面的不可约表示也可确定。而当费米面上下两个能带的不可约表示发生变化时，其拓扑性质发生变化。所以非常自然，一个绝缘体材料可以用其价带在所有高对称点的不可约表示情况来标记。当这个标记发生变化时预示着拓扑性质的转变。

以硅的能带为例子，硅的晶格属于第227 号空间群 (记为SG227) 。这个空间群的布里渊区有Γ，X，L，W这4 个高对称点。第一性原理计算得到的能带如图1(a)所示，其中费米能级取为0，费米能级以下有8 条能带。如表1 所示，Γ点具有6 个不可约表示，它们的维度分别为2，2，4，2，2，4，在图1(a)中我们标记了第一性原理计算得到的Γ点价带能级从属的不可约表示，其中Γ1表示相应能级对应于Γ点第1 个不可约表示 (即表1中抽象群G968 的第4 个不可约表示，参见表1) ，类似地，对于其他高对称点我们也标记了不可约表示的信息。根据前面的讨论，在分析材料拓扑性质时，可以不关心价带在高对称点的能级顺序以及高对称点间能带的具体色散这些细节，而只关心价带在高对称点每个不可约表示出现的次数或称为占据数nikj：这里， kj 标记第j 个高对称点，i 标记该高对称点第i 个不可约表示。当这个数字发生变化时预示着可能有拓扑性质的变化。对于硅，根据图1(a)，在Γ点，

表1 SG227 四个高对称点的坐标，以及它们小群不可约表示。这里高对称点坐标的基矢取为惯用元胞对应的倒空间基矢。根据文献[89]，高对称点小群不可约表示(Irrep)可以从某一个抽象群(记为Gnm )的不可约表示得到，高对称点小群的不可约表示只能取该抽象群的若干不可约表示[89]，我们也给出了这些不可约表示的维度(第三列括号中的数字)，比如，Γ点第1 个不可约表示为抽象群[89]G896 的第4 个不可约表示，其维度为2，Γ点第2 个不可约表示为抽象群[89]G896的第5个不可约表示，其维度为2

硅的8 个价带可被图1(a)中红色虚线分割成上下两部分，每一部分的对称性标记分别为：

显而易见，

即增加能带数目，只需对其对应的整数集做加法即可。

3 基于原子绝缘体基组的拓扑判断算法

若一个体系可以在满足对称性的情况下绝热地变化到原子极限，则这个体系等价于原子绝缘体，不表现出任何拓扑性质。设想有一个体系原胞里面只有一个原子，而这个原子只有一个s 轨道。这时其能带有如图2 的几种情况：(a)没有能带交点；(b)在高对称线上有交点；(c)在高对称点有能带交点。若费米面在图2 的

处，即2 个能带均被占据，则上述的3 种情况均可绝热变化到原子绝缘体，但是如果费米面在EF 处，体系不可能在不改变对称性的情况下变到原子绝缘体。

图2 某空间群1a Wyckoff 位置上有一个s 轨道的能带情况。考虑电子全部占据(费米面在图中紫色实线位置)，那么图中3 种情形都绝热等价于原子绝缘体。然而考虑半占据(费米面在图中紫色虚线位置)，图中3 种情况都无法退化到原子极限，分别对应于：(a)拓扑(晶体)绝缘体；(b)高对称线/面上出现能带交点；(c)高对称点上出现能带交点

对应非磁性材料，一共只有230 个空间群。我们基于对称性指标理论，高效判断拓扑材料算法的核心是构造出原子绝缘体基组[88]。值得指出的是，AI 基组与具体材料无关，对应230 个空间群，存在230 套AI 基组。相比基于图论的10403个基本能带表示 (这些基本能带表示都可以在230套AI 基组上展开) [86]，我们只需得到230 个原子绝缘体基组，然后通过实际材料不可约表示标记在AI基组上的展开系数，来判定其是否为拓扑材料[88]，具有可操作性。

4 构造原子绝缘体基组

原子绝缘体基组可以从下面3个步骤来构造：

第一步，根据空间群的Wyckoff 位置列表 (可参考网站http://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-table?from=getwp) 得到对应位群的不可约表示。例如对于SG227 Wyckoff 位置，它可以取：8a，8b，16c 到192i 等9 个值。每一个Wyckoff 位置的位群必然同构于32 个点群之一，其不可约表示也可以通过查阅群论书得到[89]。

第二步，对每一个Wyckoff 位置分别一一放入其位群不可约表示的函数，组成能带 (如果存在时间反演对称T，需要考虑T 的要求：同一个不可约表示出现两次或者两个不同不可约表示通过T联系[89]) 。通过上述实空间基组构建布洛赫和，得到BZ中所有高对称k 点小群不可约表示的占据数，从而得到

即得到该Wyckoff 位置被占据某位群不可约表示时对应原子绝缘体的不可约表示标记。对应SG227，BZ 具有4 个高对称点，这些高对称点小群以及不可约表示信息可以查阅文献得到[89]，表1 中列出了这些高对称点的坐标和不可约表示。

(a)任何一个AI 均可表示为这组AI 基矢的线性组合，且组合系数为整数，在物理上它表示不同原子极限的叠加；

(b)其基矢具有最大可能的公因子和最多的大于1的公因子数；

图4 算法流程图[88]：对于任何一材料，使用能带方法可以求出其价带在BZ高对称点不可约表示出现的情况n，然后在该空间群的AI 基组(与具体材料无关)上做线性展开，得到3 种情形，其中第一种等价于一个原子绝缘体(即文中的情况(i))，第二种对应于拓扑(晶体)绝缘体(即文中的情况(ii))，第三种对应于拓扑半金属(即文中的情况(iii))

我们判断拓扑性质的流程图如图4所示。相关的程序源代码参见文献[91]。

主要步骤如下：对于任何一个材料，基于第一性原理计算自洽的结果，可以求出价带在BZ 中各个高对称点上不可约表示的占据数。值得提及的是，这里不需要事先分析这组价带是不是和导带分开，即是否具有直接带隙。使用常规计算不可约表示出现次数的方法[91]，可以求得如(1)式形式的n，如果所有分量都是整数，说明所有高对称点上都具有有限的直接带隙，相反，如果不是整数说明这个非整数对应的高对称点具有能带交叉 (比如图2(c)) 。需要注意的是，即使所有n 的分量都是整数，n 也可能并不对应一个绝缘体 (比如图2(b)) 。我们把第一性原理计算得到材料的n 用该材料所属空间群对应的原子绝缘体基组 ( 即{ai， i=1，2，…，dAI}) 做展开：

1) n 的某个分量为非整数，说明非整数分量对应的高对称点存在能带交叉；

2) 所有n的分量都是整数，说明BZ所有高对称点都存在直接带隙，但是在高对称线/面上存在能带交叉。

其中，情形(i)中仍可能具有非平庸拓扑，需要进一步的分析，但是情形(ii)和(iii)必然具有非平庸拓扑，因此我们主要关心情形(ii)和(iii)对应的拓扑材料。

5 大规模拓扑材料搜索

众所周知，一个有中心反演对称性的材料，其原子绝缘体在演化的过程中如果其价带和导带的宇称在奇数个时间反演不变点处发生了翻转，该材料变成强拓扑绝缘体[84]，可对应于

这一方法非常高效，因此我们也用其对ICSD晶体库[94]中所有的材料进行了拓扑性质的判断和筛选[90，91]。ICSD晶体库搜集了人们已经成功生长出来的大约有18 万个晶体材料，然而其中大量的材料原子占据不是整数，基于密度泛函理论的第一性原理计算方法处理这些体系有一定困难，所以在研究中排除掉了这些体系。LDA和GGA势函数在处理电子关联强的体系有很大的局限性，我们也排除掉了含有3d、4f 以及部分4d/5d 元素 (Tc，Ru，Rh，Os，Ir) 的材料体系。5f 锕系核材料体系不利于实际应用，也被排除掉了。最终通过上述初筛，得到19143 个晶体材料，然后基于我们的拓扑材料搜索算法，对这些材料进行快速的拓扑分类，最终得到10897个有趣的材料，其中包含1123 个拓扑绝缘体，709 个拓扑晶体绝缘体和9065 个在费米面附近存在能带交点的体系 (材料的结构和能带信息见ccmp.nju.edu.cn) 。如图5 所示，大规模搜索结果表明拓扑材料并非人们以前认为的那么罕见[90]。

图5 ICSD晶体库拓扑材料大规模搜索结果统计[90]

6 典型材料分析

下面将分析我们找到的几个典型的拓扑材料来进一步展示如何利用对称性指标来判断各种各样的拓扑材料，以及在得到拓扑归类以后所需的进一步分析。

6.1 拓扑晶体绝缘体

为了验证上述高阶绝缘体特有的“体—边”对应，我们构建了同时具有空间反演和螺旋对称的棱柱 (图6(c)) ，并且计算了棱柱的电子结构，如图6(d)所示，根据波函数实空间的分布区分出来铰链态，在图6(d)中以红线表示 (两个helical 的铰链态在棱柱中的具体位置见图6(c)) 。值得提及的是，这里的结果被Wang等人的工作证实[101]。

6.2 拓扑绝缘体

在文献[90]中我们也列出一些潜在的适合应用的拓扑绝缘体材料，它们的能隙大于室温对应的热激发能 (~25 meV) ，如表4所示。

表4 我们预言的一些具有较大能隙(≥25 meV)的拓扑绝缘体材料，在“材料”列中，每一个化学式后面给出了所在的空间群

6.3 Dirac半金属

图10 Dirac半金属OPd(SG131)的电子能带图

7 总结和展望

本文详细地介绍了我们基于对称性指标理论，发展的使用原子绝缘体基组高效判断材料拓扑性质的方法，以及使用该算法对非磁性材料库进行拓扑筛选的结果 (见ccmp.nju.edu.cn) 。值得提及的是，中国科学院物理研究所方辰、翁红明、方忠团队[105]和普林斯顿Benervig，Zhijun Wang 团队[106]分别利用对称性指标[93，107]和图论方法[86]进行了晶体库的拓扑材料搜索，并且分别在http://materiae.iphy.ac.cn/[105] 和https://www.topologicalquantumchemistry.com/[106]中列出了所有发现的拓扑材料。在这些材料中，找寻性能优越具有潜在应用前景的拓扑材料将是一件非常有意思的课题。

文中介绍的方案也可以应用到二维材料，也就是需要考虑80 个layer 群[108]。此外，考虑磁性材料，该方案还可以被应用到1651 个磁性空间群[109]。除了稳定的拓扑，情形(i)中还可能存在不稳定的，所谓“易碎的”拓扑[110]，寻找具有这类拓扑的材料实现也非常有意思，而相关算法也可以很容易基于原子绝缘体基组来实现。此外，我们的方法对于构建具有期望拓扑性质的有效模型 (全BZ) 也非常有用。除了上述单电子图像下的拓扑性质，考虑多电子相互作用的拓扑材料的寻找也值得作为未来进一步的课题。

另外， 除了电子拓扑材料， 近年来， 光子、声子、磁振子等玻色系统的拓扑性质也引来了很多研究兴趣，而光子的空间群理论[111，112]也已经被提出来。基于对称性指标寻找实现拓扑超导的材料[113]，也是未来非常值得研究的课题。

计算方法说明

我们的第一性电子结构计算基于WIEN2k 软件[114]，并且考虑自旋轨道耦合，使用了Perdew—Burke—Ernzerhof (PBE) 采取广义梯度近似 (GGA) 的交换关联势[115]。在文献[90]中对挑选的好的拓扑 (晶体) 绝缘体进行了进一步mBJ (modifiedBecke—Johnson) [116]计算。

致谢 感谢Ashvin Vishwanath 和Hoi ChunPo在系列工作[88，90，91]中的合作。

参考文献

[1] v Klitzing K，Dorda G，Pepper M. Phys. Rev. Lett.，1980，45：494

[2] Thouless D J，Kohmoto M，Nightingale M P et al. Phys. Rev.Lett.，1982，49：405

[3] Haldane F D M. Phys. Rev. Lett.，1988，61：2015

[4] Hasan M Z，Kane C L. Rev. Mod. Phys.，2010，82：3045

[5] Qi X L，Zhang S C. Rev. Mod. Phys.，2011，83：1057

[6] Kane C L，Mele E J. Phys. Rev. Lett.，2005，95：226801

[7] Kane C L，Mele E J. Phys. Rev. Lett.，2005，95：146802

[8] Bernevig B A，Hughes T，Zhang S C. Science，2006，314：1757

[9] Fu L，Kane C L，Kane E J. Phy. Rev. Lett.，2007，98：106803

[10] Moore J E，Balents L. Phys. Rev. B，2007，75：121306 (R)

[11] Roy R. Phys. Rev. B，2009，79：195322

[12] Schnyder A P，Ryu S，Furusaki A et al. Phys. Rev. B，2008，78：195125

[13] Ryu S，Schnyder A P，Furusaki A et al. New Journal of Physics，2010，12：065010

[14] Kitaev A. AIP Conference Proceedings，2009，1134：22

[15] Chiu C K，Teo J C Y，Schnyder A P et al. Rev. Mod. Phys.，2016，88：035005

[16] Fu L. Phys. Rev. Lett.，2011，106：106802

[17] Hsieh T H，Lin H，Liu J et al. Nature Commun.，2012，3：982

[18] Wang Z，Alexandradinata A，Cava R J et al. Nature，2016，532：189

[19]Wieder B J et al. Science，2018，361：246

[20] Schindler F et al. Science Advances，2018，4：eaat0346

[21] Benalcazar W A，Bernevig B A，Hughes T L. Science，2017，357：61

[22] Benalcazar W A，Bernevig B A，Hughes T L. Phys. Rev. B，2017，96：245115

[23] Song Z，Fang Z，Fang C. Phys. Rev. Lett.，2017，119：246402

[24] Langbehn J，Peng Y，Trifunovic L et al. Phys. Rev. Lett.，2017，119：246401

[25] Slager R J，Mesarosm A，Juričic V et al. Nature Physics，2013，9：98

[26] Ando Y，Fu L. Annu. Rev. Condens. Matter Phys.，2015，6：361

[27] Wan X，Turner AM，Vishwanath A et al. Phys. Rev. B，2011，83：205101；万贤纲. 物理，2015，44：427

[28] Young S M，Zaheer S，Teo J C Y et al. Phys. Rev. Lett.，2012，108：140405

[29]Wang Z et al. Phys. Rev. B，2012，85：195320

[30]Wang Z，Weng H，Wu Q et al. Phys. Rev. B，2013，88：125427

[31] Burkov A A，Hook M D，Balents L. Phys. Rev. B，2011，84：235126

[32] Weng H，Liang Y，Xu Q et al. Phys. Rev. B，2015，92：045108；Fang C，Chen Y，Kee H Y et al. Phys. Rev. B，2015，92：081201

[33] Bzdusek T，Wu Q S，Ruegg A et al. Nature，2016，538：75

[34] Chen W，Lu H Z，Hou J M. Phys. Rev. B，2017，96：041102(R)；Sun X Q，Lian B，Zhang S C. Phys. Rev. Lett.，2017，119：147001；Yan Z，Bi R，Shen H et al. Phys. Rev. B，2017，96：041103(R)

[35] Chang PY，Yee C H. Phys. Rev. B，2017，96：081114(R)

[36] Ezawa M. Phys. Rev. B，2017，96：041202(R)

[37]Wang L，Jian S K，Yao H. Phys. Rev. B，2017，96：075110

[38] Bradlyn B et al. Science，2016，353：aaf5037；Weng H，Fang C，Fang Z et al. Phys. Rev. B，2016，93：241202(R)；Weng H，FangC，Fang Z et al. Phys. Rev. B，2016，94：165201

[39] Zhang T et al. Phys. Rev. Lett.，2009，103：266803

[40] Fu L，Kane C L. Phys. Rev. Lett.，2008，100：096407；AkhmerovA R，Nilsson J，Beenakker C W J. Phys. Rev. Lett.，2009，102：216404

[41] Qi X L，Hughes T L，Zhang S C. Phys. Rev. B，2008，78：195424；Essin A，Moore J E，Vanderbilt D. Phys. Rev. Lett.，2009，102：146805

[42] Son D T，Spivak B Z. Phys. Rev. B，2013，88：104412

[43] Potter A C，Kimchi I，Vishwanath A. Nature Commun.，2014，5：5161'· 354 ··48卷(2019 年) 6 期

[44] Wang C M，Sun H P，Lu H Z et al. Phys. Rev. Lett.，2017，119：136806

[45] Zhang C et al. Nature，2019，565：331

[46] Yang K Y，Lu YM，Ran Y. Phys. Rev. B，2011，84：075129

[47] Hosur P，Dai X，Fang Z et al. Phys. Rev. B，2014，90：045130

[48]Wu L et al. Nature Physics，2017，13：350

[49] Katsnelson M I，Novoselov K S，Geim A K. Nature Physics，2006，2：620

[50] Kopnin N B，Heikkila T T，Volovik G E. Phys. Rev. B，2011，83：220503(R)

[51] Zhang H，Liu C X，Qi X L et al. Nature Physics，2009，5：438

[52] Weng H，Fang C，Fang Z et al. Phys. Rev. X，2015，5：011029；Huang S M et al. Nature Commun.，2015，6：7373

[53] Xu G，Weng H，Wang Z et al. Phys. Rev. Lett.，2011，107：186806

[54] Yu R，Weng H，Fang Z et al. Phys. Rev. Lett.，2015，115：036807；Kim Y，Wieder B J，Kane C L et al. Phys. Rev. Lett.，2015，115：036806

[55] Yu R，ZhangW，Zhang H J et al. Science，2010，329：61

[56] Du Y，Wan B，Wang D et al. Scientific Reports，2015，5：14423；Du Y et al. npj Quantum Materials，2017，2：3

[57] Liu C C，FengW，Yao Y. Phys. Rev. Lett.，2011，107：076802

[58] Zhang H J et al. Phys. Rev. Lett.，2011，106：156402

[59] Qian X，Liu J，Fu L et al. Science，2014，346：1344

[60] ZhangWet al. Phys. Rev. Lett.，2011，106：156808

[61] Neupane M et al. Phys. Rev. B，2012，85：235406

[62] Kuroda K et al. Phys. Rev. Lett.，2012，108：206803

[63] Yan B et al. Europhys. Lett.，2010，90：37002

[64] Sato T et al. Phys. Rev. Lett.，2010，105：136802

[65] Xu Y et al. Phys. Rev. Lett.，2013，111：136804

[66] Bansil A，Lin H，Das T. Rev. Mod. Phys.，2016，88：021004

[67] Tang P，Zhou Q，Zhang S C. Phys. Rev. Lett.，2017，119：206402

[68] Schindler F et al. Nature Physics，2018，14：918

[69] Hsieh D et al. Nature，2009，460：1101；Chen Y et al. Science，2009，325：178；Xia Y et al. Nature Physics，2009，5：398

[70] Liu Z K et al. Science，2014，343：864；Liu Z K et al. Nature Materials，2014，13：677

[71] Lv B Q et al. Phys. Rev. X，2015，5：031013；Xu S Y et al. Science，2015，349：613；Lv B Q et al. Nature，2017，546：627

[72] Zhang C L et al. Nature Commun.，2016，7：10735；Huang Xet al. Phys. Rev. X，2015，5：031023

[73] Fei Z et al. Nature Physics，2017，13：677；Tang S et al. NaturePhysics，2017，13：683

[74] Zhu F F et al. Nature Materials，2015，14：1020

[75] Chang C Z et al. Science，2013，340：167

[76] Takane D et al. Phys. Rev. Lett.，2019，122：076402

[77] Rao Z et al. Nature，2019，567：496

[78] Sanchez D S et al. Nature，2019，567：500

[79] Ando Y. J. Phys. Soc. Jpn.，2013，82：102001

[80] Armitage N P，Mele E J，Vishwanath A. Rev. Mod. Phys.，2018，90：015001

[81] Feng W，Wen J，Zhou J et al. Computer Physics Communications，2012，183：1849

[82] Scanlon D O et al. Adv. Mater.，2012，24：2154

[83]Wang N et al. Phys. Rev. B，2014，89：045142

[84] Fu L，Kane C L. Phys. Rev. B，2007，76：045302

[85] Fang C，Gilbert M J，Bernevig B A. Phys. Rev. B，2012，86：115112

[86] Bradlyn B，Elcoro L，Cano J et al. Nature，2017，547：298

[87] Po H C，Vishwanath A，Watanabe H. Nature Commun.，2017，8：50

[88] Tang F，Po H C，Vishwanath A et al. Nature Physics，2019，15：470

[89] Bradley C J，Cracknell A P. The Mathematical Theory of Symmetryin Solids：Representation Theory for Point Groups andSpace Groups. Oxford University Press，1972

[90] Tang F，Po H C，Vishwanath A et al. Nature，2019，566：486

[91] Tang F，Po H C，Vishwanath A et al. Science Advances，2019，5：eaau8725

[92] Khalaf E，Po H C，Vishwanath A et al. Phys. Rev. X，2018，8：031070

[93] Song Z，Zhang T，Fang Z et al. Nature Commun.，2018，9：3530

[94] Hellenbrandt M. Crystallogr. Rev.，2004，10：17

[95] Brixner L H. J. Inorg. Nucl. Chem.，1962，24：257

[96] Brown B E. Acta Crystallogr.，1966，20：268

[97] Clarke R，Marseglia E，Hughes H P. Philos. Mag. B，1978，38：121

[98] Deng K et al. Nature Physics，2016，12：1105

[99] Tamai A et al. Phys. Rev. X，2016，6：031021

[100] Huang L et al. Nat. Mater.，2016，15：1155；Jiang J et al. Nat.Commun.，2017，8：13973

[101]Wang Z，Wieder B J，Li J et al. arXiv 1806.11116 (2018)

[102] Fang C，Fu L. arXiv 1709.01929 (2017)

[103] Hsu C H et al. 2D Mater.，2019，6：031004· 355 ·评述·48卷(2019 年) 6 期

[104] Tang F，Wan X. Frontiers of Physics，2019，14：4

[105] Zhang T，Jiang Y，Song Z et al. Nature，2019，566：475

[106] Vergniory M G，Elcoro L，Felser C et al. Nature，2019，566：480

[107] Song Z，Zhang T，Fang C. Phys. Rev. X，2018，8：031069

[108]Wang D et al. arXiv 1906. 01283 (2019)

[109] Watanabe H，Po H C，Vishwanath A. Science Advances，2018，4：eaat8685

[110] Po H C，Watanabe H，Vishwanath A. Phys. Rev. Lett.，2018，121：126402

[111]Watanabe H，Lu L. Phys. Rev. Lett.，2018，121：263903

[112] de Paz M B，Vergniory M G，Bercioux D et al. arXiv 1903.02562 (2019)

[113] Ono S，Yanase Y，Watanabe H. arXiv 1811.08712 (2018)

[114] Blaha P，Schwarz K，Madsen G K H et al. WIEN2k：An AugmentedPlane Wave + Local Orbitals Program for CalculatingCrystal Properties. Karlheinz Schwarz，Technische Univ. Wien，Autria，2001

[115] Perdew J P，Burke K，Ernzerhof M. Phys. Rev. Lett.，1996，77：3865

[116] Koller D，Tran F，Blaha P. Phys. Rev. B，2001，83：195134

本文经授权转载自微信公众号「中国物理学会期刊网」