反铁磁拓扑绝缘体与轴子绝缘体 ——MnBi2Te4系列磁性体系的研究进展
|作者：占国慧1 王怀强1 张海军1,2,†
(1 南京大学物理学院)
(2 南京大学固体微结构物理国家重点实验室)
本文选自《物理》2020年第12期
摘要 拓扑物质态是21世纪以来凝聚态物理领域最重要的前沿课题之一。它不仅深化了人们对宏观量子物质态的认识，同时也具有极大的潜在应用价值。目前，非磁性拓扑物质态的研究已经较为完善，而磁性拓扑物质态的研究仍处于初始阶段。近两年来，以MnBi2Te4系列体系为代表的本征磁性拓扑绝缘体的出现，迅速掀起了磁性拓扑绝缘体的研究热潮。文章从拓扑物质态的基本理论出发，介绍了近期反铁磁拓扑绝缘体方面的一些重要研究进展，着重阐述了MnBi2Te4系列的反铁磁拓扑绝缘体、静态轴子绝缘体以及动态轴子绝缘体，并对磁性拓扑绝缘体的下一步研究进行了展望。
关键词 拓扑物质态，反铁磁拓扑绝缘体，静态轴子绝缘体，动态轴子绝缘体
引 言
拓扑本是一个重要的数学几何概念，主要研究一个几何图形或者空间在进行连续改变后能够保持不变的某些性质。莫比乌斯纸带、多面体的欧拉定理、扭结理论等，都是大家广为熟知的拓扑学问题。在20世纪七八十年代，三位美国物理学家 David Thouless，Michael Kosterlitz 和 Duncan Haldane 创造性地运用拓扑理论来揭示凝聚态物理中的新奇现象，开创了全新的拓扑物理研究方向，并且由此获得了2016年的诺贝尔物理学奖。凝聚态物理中新发现的拓扑物质态不能由传统对称性破缺理论来描述，是一种全新的宏观量子态。当前，拓扑物质态和拓扑材料研究已经成为凝聚态物理领域最重要的研究内容之一。令人倍感骄傲的是，众多中国科学家站在了拓扑物理研究的科研最前沿，在国际上发挥了重要的引领作用。
近年来，人们发现了多种新奇的拓扑物质态，极大地丰富了凝聚态物理的研究内容[1—4]。2005年，量子自旋霍尔效应是由美国科学家 Kane 和 Mele 首先在理论上预言的拓扑物质态[5]，具有非零的拓扑不变量 Z2 = 1。从能谱上看，量子自旋霍尔态的体电子结构是绝缘态(具有非零体能隙)，在边界上有一对金属性的被时间反演对称性保护的一维螺旋(helical)边界态(图1(a))。在保证体系时间反演对称性和体能隙不闭合的条件下，量子自旋霍尔态不能“绝热”演化到拓扑平庸的绝缘态(Z2 = 0) ，体现了量子自旋霍尔效应的拓扑稳定性。2007年，量子自旋霍尔效应在斯坦福大学张首晟教授与合作者预言的HgTe/CdTe量子阱中被实验证实[6,7]。从此掀开了在凝聚态物理中探索拓扑物质态和拓扑材料的热潮。2010年，中国科学院物理研究所的方忠院士团队最早理论预言磁性原子掺杂的(BiSb)2Te3薄膜可以实现量子反常霍尔效应[8] (图1(b))，2013年清华大学薛其坤院士团队在该体系中成功观测到量子反常霍尔效应[9]。2009年，张首晟教授和方忠院士团队合作，理论上发现了迄今为止最具代表性的三维拓扑绝缘体Bi2Se3系列(图1(c))，引燃了国际上拓扑绝缘体研究的热潮[10]。2011年，南京大学万贤纲教授及其合作者，最早提出外尔(Weyl)半金属的概念(图1(d))，开辟了拓扑半金属的研究方向[11]。2015年，翁红明研究员、戴希教授和方忠院士团队发现TaAs系列Weyl半金属[12]，推动了国际上拓扑半金属的研究。除此之外，拓扑超导态、晶体对称性保护的拓扑绝缘体、磁性拓扑绝缘体、高阶拓扑绝缘体(图1(e))等全新拓扑物质态被相继发现[13—16]。
图1 丰富多彩的拓扑物质态 (a)量子自旋霍尔效应示意图，体系边界具有一对朝相反方向传输的螺旋边缘态；(b)量子霍尔效应/量子反常霍尔效应示意图，体系边界具有单向传输的手征边缘态；(c)典型的三维拓扑绝缘体Bi2Se3结构示意图；(d)Weyl半金属示意图，费米弧连接了一对具有相反手征的Weyl点；(e)高阶拓扑相示意图，体系在不同表面的连接处具有螺旋的棱态[16]
虽然拓扑物质态的探索取得了巨大成功，但是目前主要研究成果集中在非磁性拓扑物质态方面。而磁性拓扑绝缘体的研究尚处于初始阶段，需要人们进一步去探索和研究。主要原因是，一方面，磁性理论涉及强关联电子问题，一直是凝聚态物理中的难题之一，磁性拓扑物质态很难通过基于单电子近似的第一性原理计算准确地预言；另一方面，磁性自由度与体系空间群结合形成巨大的磁空间群，极大地丰富了磁性拓扑物质态。此外，因为实际磁性材料中经常伴随着各种磁畴和样品缺陷，合成纯净的磁性材料样品也是实验工作的巨大挑战。下面我们将介绍近期人们在磁性拓扑绝缘体研究方面取得的一些重要进展。限于篇幅，在保持内容完整性的基础上，以介绍本课题组的工作内容为主。
拓扑物质态理论
当前对于拓扑物质态的描述主要有两种理论：一种是基于单电子图像的拓扑能带理论[17,18]；另一种是基于低能长波条件下的有效拓扑场论[2,19]。拓扑能带理论主要适用于无相互作用体系，可以提供能带拓扑性质的简洁判据，定义各种拓扑不变量。拓扑场论可以刻画具有无序和相互作用的体系，并且可以用来描述与拓扑性质相关的各种物理响应行为。在无相互作用的极限下，两者对于系统拓扑性质的描述是等价的[20]。下面我们将分别介绍这两种拓扑理论。
2.1.1 量子(反常)霍尔效应与陈数
量子霍尔效应和量子反常霍尔效应都具有量子化的霍尔电导。它们的霍尔电导来源于无能隙的具有鲁棒性(robust)单向传输的一维手征(chiral)边缘态。通过体边对应关系，边缘态的数目可以通过计算体系的拓扑不变量(陈数)得到。陈数可以表示成贝利(Berry)曲率在整个布里渊区积分：
其中，n是能带指标，
分别表示定义在布洛赫波函数上的贝利曲率和贝利联络。对所有占据态的陈数进行求和就可以得到系统的总陈数，它和系统量子化霍尔电导的倍数相等。作为拓扑不变量，只要系统的体能隙不关闭，陈数就保持不变。需要注意的是，当体系具有时间反演对称性时，贝利曲率满足 Fn(k) = -Fn(-k)，从而使得陈数为0。因此要得到非零的陈数必须破坏系统的时间反演对称性，如通过外加磁场或者引入磁性长程序等。
2.1.2 拓扑绝缘体与拓扑不变量Z2
拓扑绝缘体的一个重要特征是布里渊区高对称性点附近的能带反转。比如，HgTe中s轨道主导的 Γ6 带和p轨道主导的 Γ8 带在Γ点的能带反转[6]。然而反带结构并不是体系具有非平庸拓扑性质的充要条件，需要进一步计算体系的拓扑不变量。这里我们主要介绍Kane等人基于能带理论定义的拓扑不变量Z2[18]。
对时间反演操作关联的k和-k两点的布洛赫态，可以定义如下的关联矩阵：
其中Θ代表时间反演算符。在时间反演不变点Λi处，可以证明wmn矩阵是反对称的，从而可以定义如下特征量：
这里 Pf 代表反对称矩阵的Pfaffian。因此，对于二维拓扑绝缘体，Z2不变量可以定义为所有4个时间反演不变点的δi乘积：
其中，ν = 1(ν = 0)代表拓扑非平庸(平庸)相。这里， ν = 1的二维拓扑绝缘态就是量子自旋霍尔态。
对于三维拓扑绝缘体，类似地，所有8个时间反演不变点的δi乘积可以定义出一个强拓扑不变量：
2.2.1 拓扑作用量
这里我们主要介绍三维拓扑绝缘体的低能有效拓扑场论[2,19]。通过对4(空间维度)+1(时间维度)维时间反演不变的拓扑绝缘体的 Chern—Simons作用量进行维度约化，便可以得到3+1维绝缘体的有效作用量：
其中 α 是精细结构常数，ϵμνρτ是四阶Levi—Civita张量，Fμν = ∂μAν - ∂νAμ 是电磁场张量，θ(r,t)是个赝标量场，它在场论中被称为轴子(axion)场。
轴子场θ有两个重要的性质。首先，对于具有周期性边界条件的体系，θ和θ+2π是等价的；其次，θ在时间反演或者空间反演操作下要反号变成 -θ 。因此，对于具有时间反演对称性的三维绝缘体，θ只能为0或者π，对应三维普通或拓扑绝缘体。这样拓扑场论中量子化的轴子场θ和拓扑能带理论定义的Z2拓扑不变量可以等价地描述三维拓扑绝缘体，有
2.2.2 轴子电动力学与动态轴子绝缘体
在高能物理领域，轴子是20世纪70年代为解决强相互作用中“空间宇称和电荷共轭的联合变换”守恒问题，而由Peccei—Quinn理论首先提出的一种假想粒子[21]。该理论预言，轴子具有极微小的质量(电子质量的10-11)，并且其电荷和自旋都为 0，通过引力和电磁力相互作用。轴子具有寿命长、温度低、无碰撞等类似于磁单极粒子的特性。近些年来，轴子被认为是宇宙模型中暗物质的可能构成粒子之一。可惜直到今天，轴子在自然界中是否真实存在仍然没有定论。1987年，理论物理学家Wilczek引入一个额外的麦克斯韦作用量 ΔL= κθE∙B 来研究轴子电动力学效应[22]，其中 κ 是耦合常数， θ 是轴子场， E 和 B 分别是电场和磁场。如果对比拓扑场论中三维拓扑绝缘体的拓扑作用量 Sθ 与 Wilczek 引入的轴子作用量ΔL，可以看到拓扑作用量本质上就是一个有效的轴子麦克斯韦作用量。轴子麦克斯韦作用量会导致新形式的麦克斯韦方程组：
其中 D 和 H 是普通的电位移矢量和磁场强度矢量，ρ 和 j 是通常的电荷密度和电流密度。从方程组(7)可明显地看出轴子场产生的新效应——额外的“电荷密度”和“电流密度”。磁场 B 通过空间变化的轴子场 ∇θ 诱导出额外的电荷密度，电场 E 和磁场 B 通过时空变化的轴子场( θ 和 ∇θ )诱导出额外的电流密度，反映了磁电耦合效应。我们还可以注意到，如果轴子场θ是常数，这些额外的“电荷密度”和“电流密度”必然为零，这时麦克斯韦方程组退化到通常的形式。具有随时空变化的轴子场 θ(r,t) 绝缘体被称为动态轴子绝缘体。通常情况下，传统磁性绝缘体中的轴子场非常微弱。2010年，张首晟教授与合作者理论研究表明，通过在三维拓扑绝缘体Bi2Se3中引入反铁磁序来设计强动态轴子绝缘体[23]，并且预言了轴子极化激元等新奇现象。
2.2.3 拓扑磁电效应与轴子绝缘体
三维拓扑绝缘体(θ = π)的拓扑作用量会导致很多有趣的物理效应[23—32]，包括表面的半整数量子霍尔效应[2]、拓扑磁电效应[25—28]以及拓扑磁光(Faraday和Kerr)效应[29—33](图2(f))。下面，我们介绍一下半整数量子霍尔效应。
半整数的量子化霍尔电导可以从拓扑绝缘体表面存在的奇数个无能隙的Dirac锥表面态去理解。通过在表面引入破坏时间反演的微扰项(如外磁场或者铁磁近邻效应)，可以打开原本无能隙表面态的带隙。每一个有能隙的Dirac准粒子会贡献1/2的量子化霍尔电导，从而使得表面霍尔电导为半整数。值得注意的是，表面半整数量子霍尔效应是不能够通过电子输运实验所观测的。因为一旦拓扑绝缘体的整个表面态被打开能隙，则意味着该二维表面没有开放的一维边界来承载霍尔电流。有意思的是，我们可以通过在拓扑绝缘体表面构造磁畴来实现表面态有效质量的正负转变(图2(c))，表现为整数量子霍尔效应。这也是实验上获得量子反常霍尔效应的一种可行手段。此外，实际上拓扑绝缘体总有上下两个表面，系统总的霍尔电导是这两个表面贡献之和。当上下表面磁性方向同向时，上下表面态的霍尔电导同为n/2，则此时系统的总霍尔电导为n，表现为量子反常霍尔效应(图2(a))；当上下表面磁性方向相反时，两个表面的分数霍尔电导符号相反，总的霍尔电导为0，表现为轴子绝缘态(图2(b))[19, 34]。
拓扑磁电效应是轴子绝缘体的最重要效应之一，即电场可以诱导磁极化，磁场也可诱导电极化。设计如图2(d)的圆柱体结构，当磁化强度指向圆柱体表面时，该表面态打开能隙，具有固定半整数的霍尔电导
为了实现轴子绝缘态和量子化拓扑磁电效应，需要满足三个基本要求[19, 28]：(1)体能带是绝缘的； (2)拓扑表面态要求也是绝缘的；(3)大的块材拓扑绝缘体，消除有限尺寸效应。其中，条件(1)比较容易实现，譬如通过精确调控拓扑绝缘体Bi2Se3就可以得到有能隙的体能带。条件(2)非常具有挑战性。目前实验上的做法是放宽条件(3)，通过尺寸效应和表面覆盖铁磁绝缘体来实现轴子绝缘态。但是受到尺寸效应的影响，严格来说，并不是实现了精确量子化的静态轴子绝缘体[35]。另一方面，由于在拓扑绝缘体表面覆盖铁磁绝缘体，一般情况下界面必然伴随着各种缺陷。最理想的情况是能够在本征磁性的体系中寻找轴子绝缘体，自然地满足条件(2)和(3)。下面我们将具体介绍该方面的最新研究进展 —— 反铁磁拓扑绝缘体MnBi2Te4。
03反铁磁轴子绝缘体MnBi2Te4
早在2010年，K. Mong等人指出，反铁磁绝缘体中的时间反演对称操作
首先是清华大学何珂教授与合作者率先利用分子束外延方法成功合成了MnBi2Te4薄膜[39]，他们发现交替生长Bi2Se3和MnTe层，体系会自发地形成全新的化合物MnBi2Te4，实验上通过角分辨光电子谱(ARPES)成功观测到了典型的Dirac表面态，表明该体系是一类新的磁性拓扑绝缘体，并结合理论计算指出了体系具有A型的反铁磁结构。南京大学张海军课题组和复旦大学王靖教授合作，聚焦在三维MnBi2Te4体系，理论上确认了该体系的A型反铁磁长程序，并且最早指出该体系是受S对称性保护的三维反铁磁轴子绝缘体[40]，其轴子场 θ = π 。与三维拓扑绝缘体的本质区别是，它只在特定的满足S对称的晶面上具有拓扑保护的无能隙表面态，其他晶面上表面态打开了能隙(图3(b))。如果样品表面避开S对称性的晶面，则表面态可以完全打开能隙，实现反铁磁轴子绝缘态和量子化的拓扑磁电效应。清华大学徐勇教授与合作者基于第一性原理计算[41]，深入研究了MnBi2Te4体系，独立地发现了该体系的重要磁性拓扑性质，特别是指出了该体系具有新奇的层厚奇偶性特征。国际上，西班牙Chulkov教授团队与国内同时期报道发现了MnBi2Te4的磁性拓扑绝缘体性质[42]。很快，国内外越来越多的研究组对MnBi2Te4体系展开系统的深入研究。
如图3(a)所示，MnBi2Te4是一种A-B-C型密堆叠层状磁性材料，(Te=Bi=Te=Mn=Te=Bi=Te)七原子层(septuple layer, SL)是结构的基本单元，Mn原子层位于每个SL单元的中心。SL单元内部，近邻原子之间形成强的共价键，而SL单元之间是弱的范德瓦耳斯相互作用，故而该体系易于实验剥离，是一类准二维体系。若忽略Mn原子的磁性，该体系空间对称性与三维拓扑绝缘体 Bi2Te3一致，属于空间群
这里，联合对称S“扮演着”时间反演算符
该体系还具有更为丰富的拓扑物态现象。虽然体系的磁性基态是A型反铁磁，但是由于磁矩间的交换作用不是太强，实验上完全可以通过外界磁场B来反转Mn原子的磁矩，实现铁磁结构(图4(d))。譬如，磁场B加在三度旋转轴方向时，基于第一性原理计算，我们发现在铁磁结构下，MnBi2Te4的能带结构中会在沿着磁场方向出现两个Weyl点，成为铁磁Weyl 半金属(图4(a—c))。非常有意思的是，此时体系中只有一对Weyl点，而且Weyl点精确地落在费米能级，是最理想的Weyl半金属态[40, 41]。此外， 当MnBi2Te4体系从三维降到二维薄膜时(图4(e))，具有非常有意思的奇偶特性[45]。薄膜具有奇数层SL单元时，因为相邻SL单元的磁性正好相反，所以体系会有剩余的磁矩，表现出亚铁磁性质。理论计算表明3SL及以上厚度的MnBi2Te4薄膜可以实现量子反常霍尔效应(图4(f))。薄膜具有偶数层SL单元时，体系没有剩余的磁矩，此时体系具有时间反演和空间反演的联合对称性——
图4 Weyl半金属态和陈绝缘体 (a)铁磁相MnBi2Te4体能带结构，Γ—Z(平行磁场方向)上具有一对处在同一能量处的理想Weyl点；(b)(011)方向表面的表面态，可以很明显看到kz方向上的一对Weyl点；(c)连接这对Weyl点的费米弧；(d)1SL单元结构示意图；(e)1—3SL的MnBi2Te4薄膜示意图；(f)3SL的霍尔电导，可以看到体系处在量子反常霍尔态，具有量子化的霍尔电导 σH = e2/h [40]
近期实验上取得了一系列重要进展[43,44,46]，展现了一系列新奇的磁性拓扑物态。复旦大学张远波教授和中国科学技术大学陈仙辉院士等人合作，实验上证实了Neel温度 TN = 25 K 的A型反铁磁结构[43]，率先观测到了奇数层MnBi2Te4薄膜的高温量子反常霍尔效应(1.4 K)(图5(a)和(b))。清华大学王亚愚教授和徐勇教授等人合作，发现在磁场下的MnBi2Te4薄膜(6SL)可实现轴子绝缘态与量子反常霍尔态之间的量子转变[44](图5(c)和(d))。北京大学王健教授课题组在基于MnBi2Te4的器件中实现了高陈数的量子反常霍尔态[46]。
图5 近期重要实验进展 (a, b)实验上观测到奇数层MnBi2Te4薄膜的高温量子反常霍尔效应的霍尔电阻(a)和纵向电阻(b)[43]；(c, d)实验上在6SL的MnBi2Te4薄膜中观测到的轴子绝缘体相(c)和陈绝缘体相(d)[44]
MnBi2Te4其实代表了一大类“Mn-Bi-Te”磁性体系，一般性的化学式表示为 (XB)m(A2B3)n ，其中 X =Mn/Eu， A =Sb/Bi， B =Se/Te， m 和 n 为整数。可以看出，这类化合物可以形成丰富的体系家族，譬如，已经被实验上展开深入研究的MnBi2Te4，MnBi4Te7，MnBi6Te10和MnBi8Te13[47—51](图6(b))，还有被预言为动态轴子绝缘体的Mn2Bi2Te5，Mn2Bi6Te11等[52, 53]。随着成分(A2B3)的增加，可以预期层间反铁磁作用将会逐步减弱，甚至会变为铁磁结构，更有利于外界磁场的调控。类似的材料体系如(GeTe)m(Sb2Te3)n和(GeTe)m(Bi2Te3)n已经作为重要的相变存储材料被广泛深入地研究[54]。人们可以看到该类体系的常见基本结构模块分别是5层单元A2B3，7层单元XA2B4和9层单元X2A2B5(图6(a))。基于这些基本结构模块，可以像搭建乐高积木一样，搭建出各种复杂结构。故而，Mn-Bi-Te磁性体系家族，非常值得理论和实验上继续深入探索和挖掘。
动态轴子绝缘体
我们前面提到，对于同时破坏空间反演和时间反演对称性的三维反铁磁绝缘体，轴子场不再被限制为量子化的0和π 。这类磁性体系中内禀存在的反铁磁序参量的涨落会导致动态轴子场[23,28,55,56]，表示为：θ(x,t) = θ0 + δθ(x,t)，其中θ0和δθ(x,t)分别是源于反铁磁基态的静态部分和源于反铁磁序涨落的动态部分。动态轴子场可以导致很多新奇有趣的物理效应，比如动态的手征磁效应和反常霍尔效应[22,52,57]，轴子和光线性耦合形成的轴子极化激元[23]，以及由轴子失稳诱导的电磁效应[58,59]。
目前，在动态轴子绝缘体领域有两个重要的科学问题依然悬而未决。其一，拓扑分类的问题。由于时间反演对称性被破坏，拓扑不变量Z2已经不再适用，能否从拓扑上区分以拓扑绝缘体为母体的动态轴子态和传统反铁磁绝缘体的动态轴子态？我们定义基于赝自旋陈数的隐含拓扑不变量来区分这两类动态轴子态[53]。拓扑非平庸的动态轴子绝缘体具有各向同性的强动态轴子场(θ0~ π)，而拓扑平庸的动态轴子绝缘体具有非常微弱的轴子场 (θ ~ 0) ，譬如传统反铁磁绝缘体氧化铬[60]。其二，尽管基于理论模型，人们预言了很多和动态轴子场有关的新奇物理效应，但是由于缺乏实际的材料，实验上一直没有能够直接观测和验证这些现象。我们理论预言了拓扑非平庸的动态轴子绝缘体Mn2Bi2Te5系列体系[52]以及MnBi2Te4/Bi2Te3异质结系列体系[53]，有望推动该方面的实验工作。
不失一般性，我们从三维拓扑绝缘体Bi2Te3的低能有效模型出发[10]，然后把反铁磁项看成系统微扰。在非磁态，体系具有空间反演对称性，能带中4个低能态来源于Bi原子p轨道的成键态
反铁磁序的存在同时破坏空间反演
通过相关理论计算可得三维绝缘体中的轴子场θ 的值。反铁磁质量项 m5 会导致对 θ的线性修正，即：
基于理论分析，我们指出获得强动态轴子场的一般性办法：首先，要求同时破坏时间和空间反演对称性、接近拓扑转变点、具有较小反带质量m的三维反铁磁绝缘体；其次，要求体系接近具有强磁涨落的反铁磁和顺磁转变点 (m5 ≈ 0) 。因此，我们认为具有弱反铁磁序并且接近拓扑相变点的本征反铁磁拓扑绝缘体具有机会实现强动态轴子绝缘体[52]。
由于反铁磁质量项m5的存在，当反带质量m从负(拓扑绝缘体区域)变到正(普通绝缘体区域) 时，体能隙一直不关闭。因此严格意义上没有发生电子态上的拓扑相变。但是我们发现，可以定义源于赝自旋陈数的隐含拓扑不变量来从拓扑上分类动态轴子绝缘态[53]。(赝)自旋陈数的概念最早主要是为了研究量子自旋霍尔效应在时间反演对称性被破坏时的拓扑稳定性[61]。我们把kz看成参数去计算二维kx—ky平面的赝自旋陈数。针对有效哈密顿量(8)式，我们得到kz依赖的赝自旋陈数表达式：
可以看到，赝自旋陈数不依赖于反铁磁质量项m5。由于B1和B2被设定为正数，如果反铁磁绝缘体质量项 m < 0 ，以kz为参量，必然存在赝自旋陈数 CS 的拓扑转变；而如果反铁磁绝缘体 m > 0 时，则不存在赝自旋陈数的拓扑转变。质量项m本质上反映了反铁磁绝缘体的能带反转特征。进一步我们可以定义一个类似Z2的隐含拓扑不变量ν：
ν = -1 (1)对应反铁磁绝缘体是拓扑非平庸(平庸)的动态轴子绝缘体。ν不是直接定义在电子态上的拓扑不变量，而是定义在电子态伴随的赝自旋态上的隐含拓扑不变量。有意思的是，ν = -1 可以直接与电子态中的能带反转相对应。意味着，通过在拓扑绝缘体中进行弱磁性掺杂，只要不改变体系的能带反转特征，必然得到具有非平庸的动态轴子绝缘体，譬如反铁磁掺杂的Bi2Se3模型[23]。
我们提出了两类有希望实现拓扑非平庸动态轴子态的磁性拓扑材料。一类是以Mn2Bi2Te5为代表的系列体系，它的基本结构是具有两层磁矩相反Mn原子的9层单元。另外一类是超晶格(MnBi2Te4)m(Bi2Se3)n系列体系，其中选取的 (m,n) 值必须使得体系破坏空间反演对称性，例如反铁磁结构超晶格(MnBi2Te4)2(Bi2Se3)[52, 53]。下面，我们将以此超晶格体系为例说明其动态轴子场性质。
(MnBi2Te4)2(Bi2Se3)的结构如图7(a)所示，它的基本单元是由一层Bi2Se3穿插在2SL的MnBi2Te4间形成的三明治结构。该结构的磁基态构型保持A型反铁磁，每SL层内部是铁磁序，磁矩平行于面外，相邻SL层间是反铁磁序。该体系同时破坏了时间反演
为了研究该体系中轴子场的性质，我们从低能有效模型出发，在图7(b)中给出了轴子场的静态部分θ0(蓝线)以及动态部分与反铁磁序涨落的比值1/g(红线)随自旋轨道耦合强度变化的曲线图。很明显，静态θ0偏离0和π，并且随着自旋轨道耦合强度的增加而增加，θ0从接近0的值逐渐增加到接近π的值。对于动态部分，1/|g| 在 λ > 0.68 的反带区域要比 λ < 0.68 的非反带区域大一个数量级。1/|g| 在从反带区域接近转变点的位置取到极大值，和我们的理论预期是一致的。另外，作为MnBi2Te4体系的姊妹体系Mn2Bi2Te5，也已经在实验上成功合成。因此，Mn2Bi2Te5和(MnBi2Te4)m(Bi2Se3)n系列磁性体系，有希望实现拓扑非平庸的强动态轴子态，为进一步实验探测和研究与轴子电动力学相关的新奇效应提供理想的平台。
总结与展望
本文简要回顾了磁性拓扑绝缘体的基本知识，重点介绍了反铁磁拓扑绝缘体的最新研究进展。 以MnBi2Te4为代表的系列磁性绝缘体为例，结合本课题组近期在理论方面的几个工作，着重讨论了反铁磁拓扑绝缘体、静态轴子绝缘体和动态轴子绝缘体。由于篇幅所限，未能涵盖磁性拓扑绝缘体全面相关的重要工作。当前反铁磁拓扑绝缘体领域的研究正在如火如荼地开展，重要研究成果不断涌现，但是其中也还有不少需要亟待解决的科学问题。譬如，尽管理论预言MnBi2Te4的表面态是有能隙的，但是最近几个课题组的角分辨光电子谱(ARPES)实验都观测到了无能隙的表面态[62—64]，这一现象需要进一步的理论解释和实验澄清，也许和复杂的表面磁结构有关。Mn-Bi-Te系列磁性体系的发现只是磁性拓扑物理研究的开端，预期更多新的磁性拓扑物理效应和理想磁性拓扑体系会被发现。最后，需要强调的是，动态轴子绝缘体的研究还处在理论研究的初期阶段，随着候选体系Mn2Bi2Te5和(MnBi2Te4)m(Bi2Te3)n的实验合成，必将开启人们在动态轴子态新奇电磁响应方面的实验研究，推动凝聚态物理与光物理的交叉研究，甚至可以跨领域推动宇宙暗物质的探测。
[2] Qi X L，Zhang S C. Rev. Mod. Phys.，2011，83：1057[3] Chiu C K，Teo J C Y, Schnyder A P et al. Rev. Mod. Phys.，2016，88：035005[4] Armitage N P，Mele E J，Vishwanath A. Rev. Mod. Phys.，2018，90：015001[5] Kane C L，Mele E J. Phys. Rev. Lett.，2005，95：226801[6] Bernevig B A，Hughes T L，Zhang S C. Science，2006，314：1757[7] König M，Wiedmann S，Brüne C et al. Science，2007，318：766[8] Yu R，Zhang W，Zhang H J et al.Science，2010，329：61[9] Chang C Z，Zhang J S，Feng X et al. Science，2013，340：167[10] Zhang H J，Liu C X，Qi X L et al. Nat.Phys.，2009，5：438[11] Wan X G，Turner A M，Vishwanath A et al. Phys.Rev. B，2011，83：205101[12] Weng H M，Fang C，Fang Z et al. Phys. Rev. X，2015，5：011029[13] Qi X L，Witten E，Zhang S C et al. Phys. Rev. B，2013，87：134519[14] Fu L. Phys. Rev. Lett.，2011，106：106802[15] Tokura Y，Yasuda K J，Tsukazaki A. Nat. Rev. Phys.，2019，1：126[16] Schindler F，Wang Z J，Vergniory M G et al. Nat. Phys.，2018，14：918[17] Fu L，Kane C L，Mele E J. Phys. Rev. Lett.，2007，98：106803[18] Fu L，Kane C L. Phys. Rev. B，2007，76：045302[19] Qi X L，Hughes T L，Zhang S C. Phys. Rev. B，2008，78：195424[20] Wang Z，Qi X L，Zhang S C. New. J. Phys.，2010，12：065007[21] Peccei R D，Quinn H R. Phys. Rev.Lett.，1977，38：440[22] Wilczek F. Phys. Rev. Lett.，1987，58：1799[23] Li R D，Wang J，Qi X L et al. Nat. Phys.，2010，6：284[24] Qi X L，Li R D，Zang J D et al. Science，2009，323：1184[25] Karch A. Phys. Rev. Lett.，2009，103：171601[26] Nomura K，Nagaosa N. Phys. Rev. Lett.，2011，106：166802[27] Maciejko J，Qi X L，Drew H D et al. Phys. Rev. Lett.，2010，105：166803[28] Wang J，Lian B，Qi X L et al. Phys. Rev. B，2015，92：081107[29] Tse W K，MacDonald A H. Phys. Rev. Lett.，2010，105：057401[30] Wu L，Salehi M，Koirala N et al. Science，2016，354：1124[31] Ochiai T. J. Phys. Soc. Jpn，2012, 81：4401[32] Mal'shukov A G，Skarsvag Hans，Brataas Arne. Phys. Rev. B，2013，88：245122[33] Okada K N，Takahashi Y，Mogi M et al. Nat. Comm.，2016，7：12245[34] Mogi M，Kawamura M，Yoshimi R et al. Nat. Mater.，2017，16：516[35] Xiao D，Jiang J，Shin J H et al. Phys. Rev. Lett.，2018，120：056801[36] Mong R S K，Essin A M，Moore J E. Phys. Rev. B，2010,81：245209[37] Lee D S，Kim T H，Park C H et al. Cryst. Eng. Comm，2013，15：5532[38] Otrokov M M，Menshchikova T V，Vergniory M G et al. 2D Mater.，2017，4：025082[39] Gong Y，Guo J W，Li J H et al. Chin. Phys. Lett.，2019，36：076801[40] Zhang D Q，Shi M J，Zhu T S et al. Phys. Rev. Lett.，2019，122：206401[41] Li J H，Li Y，Du S Q et al. Sci. Adv.，2019，5：eaaw5685[42] Otrokov M M，Klimovskikh I I，Bentmann H et al. Nature，2019，576：416[43] Shi Z M，Guo Z X，Xu Z H et al. Science，2020，367：895[44] Liu C，Wang Y C，Li H et al. Nat. Mater.，2020，19：522[45] Otrokov M M，Rusinov I P，Blanco- Rey M et al. Phys. Rev. Lett.，2019，122：107202[46] Ge J，Liu Y Z，Li J H et al. Nat.Sci. Rev.，2020，nwaa089[47] Aliev Z S，Amiraslanov I R，Nasonova D I et al. J. Alloy. Compd.，2019，789：443[48] Klimovskikh I I，Otrokov M M，Estyunin D et al. npj Quantum Mater.，2020，5：54[49] Shi M Z，Lei B，Zhu C S et al. Phys. Rev. B，2019，100：155144[50] Tian S J，Gao S，Nie S M et al. Phys. Rev. B, 2020，102：035144[51] Hu C W，Ding L，Gordon K N et al. Sci. Adv., 2020，22：eaba4275[52] Zhang J L，Wang D H，Shi M J et al. Chin. Phys. Lett.，2020，37(7)：077304[53] Wang H Q，Wang D H，Yang Z L et al. Phys. Rev. B，2020，101：081109[54] Wuttig M，Yamada N. Nat. Mater.，2007，6：824[55] Wang J，Lian B，Zhang S C. Phys. Rev. B，2016，93：045115[56] Wang J L，Lei C，MacDonald A H et al. 2019，arXiv：1901.08536. [57] Sekine A，Nomura K. Phys. Rev. Lett.，2016，116：096401[58] Ooguri H，Oshikawa M. Phys. Rev. Lett.，2012，108：161803[59] Imaeda T，Kawaguchi Y，Tanaka Y et al. J. Phys. Soc. Jpn.，2019，88：024402[60] Coh S，Vanderbilt D，Malashevich A et al. Phys. Rev. B，2011，83：085108[61] Sheng D N，Weng Z Y，Sheng L et al. Phys. Rev. Lett.，2006，97：036808[62] Hao Y J，Liu P f，Feng Y et al. Phys. Rev. X，2019，9：041038[63] Li H，Gao S Y，Duan S F et al. Phys. Rev. X，2019，9：041039[64] Chen Y J，Xu L X，Li J H et al. Phys. Rev. X，2019，9：041040
磁性拓扑材料专题：
磁性外尔半金属材料研究现状与展望
《物理》在淘宝店和微店上线，扫码即可购买过刊和现刊。