磁场中拓扑物态的量子输运*

强晓斌 卢海舟†

1) (南方科技大学物理系, 深圳量子科学与工程研究院, 深圳 518055)

2) (深圳市量子科学与工程重点实验室, 深圳 518055)

1 引 言

拓扑半金属的导带和价带在节点(节线)处接触[1,2], 节点处低能激发的准粒子(外尔/狄拉克费米子)具有线性的能带色散关系, 非常类似于三维(3D)的石墨烯模型.大量的理论[3−13] 和实验观测[14−25]均给出了相关材料存在的证据.

拓扑物态由于具有独特的能带结构和表面态行为导致了许多新颖的输运现象[26−68].本文基于在Front.Phys.发表的两篇论文[69,70], 以拓扑半金属为主回顾了有关拓扑物态中量子输运的最新工作[71−87], 相关的综述文章也可以参考[88−93].依据磁场强度的大小本文的结构安排如下: 第2 节, 介绍用于描述外尔、狄拉克、节线半金属及拓扑绝缘体的相关模型; 第3 节, 介绍非线性霍尔效应的最新进展, 并给出半导体-超导体纳米线中马约拉纳振荡的理论描述; 第4 节, 总结拓扑半金属和拓扑绝缘体中的弱局域化和弱反局域化理论; 第5 节,给出贝里曲率与负磁阻之间的关系, 介绍平面霍尔效应的最新工作; 第6 节, 讨论拓扑半金属中量子振荡的相移规则, 给出两种实现3D 量子霍尔效应的新机制; 第7 节, 总结量子极限下的输运理论,包括背散射禁止导致的电阻下降, 外尔费米子的湮灭, 外尔半金属中的不饱和磁化现象, 以及反常热电行为; 第8 节给出评论和展望.

2 有效模型

2.1 外尔半金属

外尔半金属的最小模型为[94]

其中 σx,y,z是泡利矩阵, k =(kx,ky,kz) 是波矢, A ,M , kw是模型参数.两个能带在 (0,0,±kw) 处相交(图1).

在两个节点 (0,0,±kw) 附近, H 约化为两个独立的局部模型

图1 拓扑半金属的能带结构和贝里曲率 (a)拓扑半金属的能谱示意图, (kx,ky,kz) 为 波矢, ; (b)贝里曲率矢量场, 拓扑半金属的导带和价带在外尔点处接触,且在该处存在一对单极子.转载自文献 [56]Fig.1.The band structure and Berry curvature of the topological semimetal: (a) The energy spectrum of a topological semimetal, (kx,ky,kz) is the wave vector, (b) the vector field of the Berry curvature.The conduction and valence bands of a topological semimetal touch at the Weyl nodes, and there is a pair of monopoles.Reproduced with permission from Ref.[56].

贝里曲率[95]Ω (k)=∇k×A(k) 可以给出体系的拓扑性质, 其中 A (k)=i〈u(k)|∇k|u(k)〉 为贝里联络.对于外尔半金属的+能带|u(k)〉=[cos(Θ/2),sin(Θ/2)eiφ]T(其中 cos Θ ≡Mk/E+, tan φ ≡ky/kx).对单个外尔点做曲面积分, 结果表明在 ± kw处产生了符号相反的拓扑荷 ∓ sgn(M) , 对应动量空间中的一对“磁单极子”.

在离散情况下, 陈数定义为在任意可定向闭曲面上正方格子的贝里曲率通量之和[96].这里取具有轮胎面拓扑结构的连续参数空间, 原因是某些物理参数空间实际上确实具有这种拓扑结构, 且对应的陈数具有物理意义.对于2D 晶格的布里渊区,其中晶格动量 (kx,ky) , (kx+2π,ky) 及(kx,ky+2π)是等价的.自然地, 连续情况下对贝里曲率通量的求和被整个参数空间 P 上对贝里曲率的面积分所代替

由于其可以解释为离散陈数的连续极限, 所以(3)式具有和离散陈数相同的性质, 即连续陈数是一个规范不变的整数.

对于外尔半金属, 若给定 kz, 陈数可以表征kx-ky平面上的拓扑性质[97],

当 − kw＜kz＜kw时, 陈数 nc(kz)=−sgn(M) , 否则nc(kz)=0[4].根据体边对应关系, 非零陈数对应于与 kz相关的一组边界态(费米弧)[4].

为了得到费米弧表面态的解, 在 y =0 处应用开放边界条件, 得到表面态的色散[76,94]

对应的波函数[98,99]

外尔半金属的单节点有效模型为

其中得到贝里联络为 A =(0,0,cos2(θ/2)/(k sin θ)).由此得到贝里曲率为, 求得单极子荷 N =−1.在另一个手性相反的外尔点处, 有效哈密顿量 H =−vk·σ , 单极子荷为1.因此, 对于双节点模型总单极子荷为零[26].

在 z 方向的磁场中, 外尔半金属的能谱被量子化为一组沿 kz色散的一维(1D)朗道能带(图2(a)).考虑沿 z 方向施加磁场 B =(0,0,B) , 选择朗道规范 A =(−yB,0,0).得到本征能量为[74,100−102]

2.2 狄拉克半金属

对于狄拉克半金属, 由于时间反演和反演对称性的保护, 狄拉克点是多重简并的[8].其模型可以用方程(1)中的 H (k) 及其时间反演 H∗(−k) 描述[60,61,103−105], 其中星号表示复共轭.一个简单的模型如下[76]:

其中狄拉克矩阵为 αx=σx⊗σx, αy=σx⊗σy,β =σz⊗σ0, σ0是 2 ×2 的单位矩阵.方向的表面态由两个具有相反自旋和相反有效速度的分支组成.通过幺正变换将模型写成对角形式

狄拉克点是四重简并的[94], 对于特定的 kz两个狄拉克点之间存在拓扑不变量, 根据体边对应关系,每个陈数对应一个边界态, 两个节点之间存在关于kz的一组边界态.这些边界态的色散构成了在费米面附近连接狄拉克点的费米弧, 这是拓扑狄拉克半金属的关键特征之一.

2.3 双外尔半金属

外尔半金属中的每个外尔点都具有1 或–1 的单极子荷.在双外尔半金属中, 单极子荷为2 或–2[6,106,107].对于外尔半金属或双外尔半金属的单一能谷, 最小模型为

其中 k±=kx±iky, χ =±1 是能谷指标, k 为外尔点处的动量.N =1,2 分别对应单和双外尔半金属.色散关系为在能谷 χ =+1 处导带的本征态为

其中 cos θ ≡vzkz/Ek, tan φ ≡ky/kx.通过在包围外尔点的任意费米球 Σ 上对贝里曲率积分得到单极子荷

等式右边 ± 对应 ± 能谷.

图2 在沿 z 方向的磁场 B 下, 外尔和狄拉克半金属的最小模型中沿 kz 色散的朗道能带.转载自文献[74]Fig.2.The Landau energy bands along the kz dispersion in the minimum model of Weyl and Dirac semimetals under the magnetic field B along the z direction.Reproduced with permission from Ref.[74].

2.4 拓扑节线半金属

节线半金属[108−111]的导带和价带在动量空间交叉形成一个闭合环[112−115], 并形成类似鼓形的表面态[116,117].当对称性破缺时, 节线半金属可能变为狄拉克半金属, 拓扑绝缘体或表面陈绝缘体[118].理论预测了不少存在这种奇异拓扑性质的材料[119−139], 部分已被实验观测所证实[140−146].

节线半金属可以用一个双能带有效模型来描述[147],

其 中 k =(kx,ky,kz) 是波矢, τx,y对应于双带赝自旋空间的泡利矩阵, 两个带在以及kz=0 时相交.如果费米能满足通过变量代换kx=(k0+κ cos φ)cosθ,ky=(k0+κcosφ)×sinθ, kz=κsinφ/α (其中 v0=2λk0, α=v/v0, v 和v0分别为沿着 z 方向和x-y 平面内的有效速度), 哈密顿量(15)式可以线性化为更为简单的形式:

对于更一般的四带模型, 节线半金属的哈密顿量可以进一步写为[124,143]

其中 λ , m , 和 u 是模型参数[124,143].哈密顿量的本征值为当 u 为正时, 两个能带在零能处相交, 其半径为时(对应于载流子密度较低的情况[140,141]), 费米面形成一个轮胎面, 当 EF＞u 时,费米面演化成鼓状结构.模型(17)式具有镜面反射对称性[124,143], 除此之外, 节线也可由其他对称性保护[110,113,121,123,140,141,148].

2.5 拓扑绝缘体

3D 拓扑绝缘体可以用 k·p 哈密顿量来描述[94,149,150]:

对于本征磁性拓扑绝缘体, 例如MnBi2Te3,具有面外反铁磁序, 需要额外引入[156]

其中 m0是层内反铁磁序的强度, d 是MnBi2Te3样品7 元层的厚度, τ0是 2 ×2 的单位矩阵, 这一模型可以解析求解, 为理解其他的本征磁性拓扑物态提供理论基础.

3 零磁场

3.1 非线性霍尔效应

霍尔效应由于本质上与对称和拓扑的深刻联系[95,157]已然成为凝聚态物理的基本研究范式.所有已知的可测霍尔效应都需要磁场或磁性掺杂来破缺时间反演对称性[157−159].最近, 通过测量霍尔电导的零频或倍频分量对低频驱动电场的响应(图3)发现了一种新的霍尔效应[160], 即非线性霍尔效应, 其特别之处在于不需要破缺时间反演对称性, 而是破缺反演对称性.这一新现象的发现为探索新奇物态的对称性和拓扑性质开辟了新的途径.在文献[82]中, 我们研究了非线性霍尔效应和二维(2D)系统能带特征之间的关系.由于非线性霍尔响应与贝里偶极子相关[160], 通过分析倾斜有质量狄拉克费米子的一般模型, 我们发现倾斜能带在反交叉和能带反转附近存在非常强的贝里偶极子,可以给出明显的非线性霍尔信号.

图3 测量非线性霍尔效应的示意图.转载自文献[82]Fig.3.Schematic of how to measure the nonlinear Hall effect.Reproduced with permission from Ref.[82].

其中 e 为电子电量, τ 为动量弛豫时间, 在一般实验条件下 ω τ ≪1 可忽略.贝里偶极子定义为

其中 ε abc 为全反对称张量, a ,b,c,d ∈{x,y,z} , feq为费米分布函数,为贝里曲率.注意在时间反演下以及,即(21)式在时间反演对称的体系中不为0.

由(21)式知, 为了给出不为0 的贝里偶极子,需要各向异性的能带结构, 通过在2D 有质量狄拉克模型[94]中引入倾斜项来实现这一点,

其 中 (kx,ky) 为波矢, σx,y,z为泡利矩阵, η =±1 ,m是能隙, t 使狄拉克锥沿 x 方向倾斜即破缺了反演对称性.如果没有倾斜(t =0), 贝里曲率将对称地集中在能带边缘, 因此没有贝里偶极子.当狄拉克锥倾斜时, 贝里曲率不再对称分布, 贝里偶极子随着倾斜的增加而增加, 并在狄拉克锥完全倾斜(t ＞v) 的 临界点 (t=v) 附近达到最大值.利用这一模型我们计算了双层 WTe2中的非线性霍尔响应[82], 给出的结果与最近实验[161,162]中的非线性霍尔响应及其角度依赖性高度一致, 这将启发更多关于非线性霍尔效应的相关实验和理论.

除了由倾斜能带引起的贝里偶极子导致的本征非线性霍尔效应以外, 无序在其中也起着极其重要的作用.在文献[83]中我们系统地计算了由无序导致的非线性霍尔效应, 并给出了霍尔电导的表达式.不同于一般的线性霍尔效应, 由于非线性霍尔效应的机制需要费米能量切过能带, 主要的物理发生在费米面上, 这就导致即使在领头阶无序也扮演着非常重要的角色.更重要的是, 我们提出了非线性霍尔效应的标度律, 这有助于识别不同机制的贡献, 解释未来实验中温度和厚度的依赖性.需要指出的是, 无序是如何以一种特定的形式对非线性霍尔信号做出贡献的, 目前尚不清楚, 这也是未来研究的重点.

3.2 马约拉纳振荡

在半导体-超导体纳米线中的马约拉纳零能模[163−168]总是成对出现, 并位于纳米线的两端.由于在真实纳米线中的有限尺度效应, 马约拉纳束缚态的能量是杂化的.杂化能量 E0因为是塞曼能量、化学势以及纳米线线长的函数而出现振荡[169−171],称为马约拉纳振荡.然而实验中观察到 E0作为磁场的函数[172−177]其振荡与理论预测的截然相反[170].

我们通过假设沿纳米线的自旋轨道耦合具有阶梯状分布[84](图4(a))解释了实验中的振荡模式[172−177](图4(b)—图(d)).这种自旋轨道耦合的阶梯状分布的假设是合理的, 因为栅极施加的是非均匀的静电势, 而自旋轨道耦合又取决于垂直于纳米线的静电场[178−180].此外, 由于超导体和半导体之间的屏蔽效应以及功函数的失配, 超导体的存在可以极大地改变纳米线中的静电场[181].因此, 纳米线从覆盖有超导体的部分到没有超导体的部分(隧道势垒区), 自旋轨道耦合的非均匀性是很好的理论假设, 但这一点被大多数理论所忽略, 导致了对马约拉纳振荡的理论理解困难.

图4 (a)半导体-超导体纳米线结构示意图[172−177], 两端可能存在一对马约拉纳束缚态; (b)−(d)杂化能随着磁场变化的振荡曲线.红色和黑色曲线为实验数据[172], 蓝色为理论曲线.转载自文献[84]Fig.4.(a) Schematic of the semiconductor-superconductor nanowire structure[172−177], its two ends may host a pair of Majorana bound states; (b)−(d) oscillation curves of hybridization energy vary with magnetic field.The red and black curves are experimental data adapted from Ref.[172].The blue curves are the theoritical results.Reproduced with permission from Ref.[84].

为了验证其中的物理图像, 使用阶梯状自旋轨道耦合进行计算.纳米线的哈密顿量为

参照GB-T14452-93标准A型试样，将钎焊后试样线切割为长约34 mm、宽5±0.25 mm、高5±0.25 mm的条状小试样进行三点弯曲试验，测定钎焊接头的抗弯强度。

其中 A ,α0,xL/R,λL/R为描述自旋轨道耦合随空间分布变化的参数.在实际实验中, 当改变栅极电压时, 参数会相互影响[181−183], 超导体的存在也会引起参数的变化[184,185].通过对角化晶格上的 H , 得到能谱和波函数, 其中最低的能量即束缚态杂化能量为 E0.

图4(b)—图4(d)中的蓝色曲线分别给出了使用三组模型参数的数值结果.数值结果与实验保持了很高的一致性, 不仅在衰减振幅方面, 而且还包括最低能量交叉(图4(b)和图4(c))、反交叉(图4(d))以及磁场中振荡周期的增加(图4(c)).注意到我们的结果是一般性的, 不依赖于具体的参数.

文献[84]进一步分析了这种阶梯状的自旋轨道耦合也能导致安德列夫束缚态光谱的衰减振荡.对于由马约拉纳束缚态影响的库仑阻塞[186−190], 预测峰间距通过 π /2 的相移与峰高相关联, 这在最近的实验中是不明确的, 但依然可以用阶梯状自旋轨道耦合来解释.期望这一结果将激发更多的工作来重新检验这种通常存在于真实实验装置中的非均匀自旋轨道耦合效应.

4 近零磁场

4.1 外尔半金属中的弱(反)局域化

弱(反)局域化是存在于无序金属[191]中由量子干涉导致的输运现象.如果量子干涉的修正为正, 即产生弱反局域化现象, 如果为负则产生弱局域化现象.无论是弱局域化还是弱反局域化本质上均取决于体系的对称性(表1).外尔费米子单一能谷的低能描述为, 只具有时间反演对称性.拓扑半金属中已广泛观察到弱反局域化现象[41−43,46,55,56,60,61].这一节介绍处理这一问题的基本方法[73]和主要结论, 最近的相关研究[192,193]利用相似的方法系统地研究了更一般的3D 狄拉克材料中的弱(反)局域化行为, 给出了与实验非常符合的理论结果.

表1 对称类(正交、辛和幺正)[194]与弱局域化(WL)和弱反局域化(WAL)之间的关系[195].转载自文献 [88]Table 1.The relation between the symmetry classes (orthogonal, symplectic and unitary) [194] and weak localization (WL) and anti-localization (WAL) [195].Reproduced with permission from Ref.[88].

图5 总结了用于研究量子干涉和相互作用引起的弱局域化和反局域化的费曼图[73].在计算中,对电导率有贡献的主要有三项, 包括半经典的电导率(图5(a))、量子干涉修正(图5(b))以及电子-电子相互作用修正(图5(c)).

图5 在无序(虚线)和电子-电子相互作用(波浪线)下,计算3D 外尔半金属电导率的费曼图[71,72,191,196−199], 有向直线代表格林函数.转载自文献 [73]Fig.5.In the disorder (dashed lines) and electron-electron interaction (wavy lines), the Feynman diagram[71,72,191,196−199]of the conductivity of 3D Weyl semimetal, and the directed line represents the Green"s function.Reproduced with permission from Ref [73].

图6 不同条件下的磁导 δ σqi(B) 对参数的依赖关系 (a) η I =η∗=0 时不同的相干长度 lϕ ; (b) η ∗=0 时不同的 ηI ; (c)有限η∗ 时不同的 ηI ; (d) ηI 和 η∗ 之间的差异, 其中 ηI 与能谷间散射相关, 而 η∗ 与能谷内散射相关.虚线表示两个散射过程的相关性,ν =±是能谷指标.转载自文献[73]Fig.6.The dependence of magnetoconductivity δ σqi(B) on parameters under different conditions: (a) Different coherence length lϕ at η I =η∗=0 ; (b) different ηI at η ∗=0 ; (c) different ηI at finite η∗ ; (d) the difference between ηI and η∗ , where ηI is correlated with intervalley scattering and η∗ is correlated with intravalley scattering.The dashed lines indicate the correlation between the two scattering processes.ν =± is the valley index.Reproduced with permission from Ref.[73].

图6 给出了理论计算得到的主要结果.可以看出, 在没有能谷间散射的情况下, 对于较大的 lϕ(相位相干长度), 磁导 δ σqi(B) 为负且与成正比, 显示了3D 外尔费米子的弱反局域化特征.这种的相关性与实验非常符合[41,42].随着 lϕ变短,从到 − B2的变化是很明显的, 最终 δqi(B) 在lϕ和平均自由程 le相等时消失, 因为该系统不再处于量子干涉状态而进入半经典扩散状态.

4.2 双外尔半金属中的弱局域化

对于单一能谷连接费米球上从 k 到 − k 的背散射中间态的每条路径(在图7(a)中标记为 P), 均在原点处包含了单极子荷, 并且存在一个对应的时间反演 P′.量子干涉是由两个时间反演路径 P 和P′之间的相位差决定的, 该相位差为沿 P 和≡−P′形成的环路累积的贝里相位[95,196,200−202], 即从−k 到 k 的对应路径:

值得注意的是, 该贝里相位仅取决于单极子荷, 而不是特定的环路形状[75]对于双外尔半金属, 单极子荷 N =2 , 贝里相位为 2π.在 2π 的贝里相位下,时间反演散射态之间产生相长干涉, 从而导致弱局域化效应.但是, 对于单外尔半金属, 单极子荷为N =1 且贝里相位为 π , 导致弱反局域化效应.这里对贝里相位的分析与前文提到的对称性分类得到的结果是一致的[203].由于贝里相位是单极子荷产生的贝里曲率场的直接结果, 因此我们在弱(反)局域化效应与单极子荷 N 的奇偶性之间建立了直接的联系.

通过计算最大交叉图得到对双外尔半金属电导率的修正.最大交叉图的核心计算可以表达为粒子-粒子关联, 称为Cooperon, 对双外尔半金属其为[75]

其中 q =k1+k2k1k2是Cooperon 波矢, 和分别为初态和末态的波矢, φ1和 φ2为相应波矢的方位角, D1=8τEFv///(3π) 和是扩散系数,NF是态密度, τ 是动量弛豫时间.

图7 3D 拓扑半金属动量空间中的费米球, 其中位于原点的点表示单极子荷 N (a) P 表示从波矢 k 到标记为(k1,k2,··· ,kn)的中间态的背散射, P ′ 表示 P 的时间反演;(b) P 和 P ′ 之间的相位差等效于在环路 周围累计的贝里相位.转载自文献 [75]Fig.7.The Fermi sphere in 3D topological semimetal momentum space, where the dot at the origin represents monopole charge N : (a) P is the backscattering from the wave vector k to − k via intermediate states labeled as (k1,k2,··· ,kn), P ′ represents the time-reversal of P ; (b) the phase difference between P and P ′ is equivalent to the Berry phase accumulated around loop .Reproduced with permission from Ref.[75].

当 q →0 , 即 k1=−k2时, Cooperon 发散, 成为背散射的最主要贡献.在此极限下, φ2=φ1+π (其中k cos θ , 并通过设置 φ →φ+π 和 θ →π−θ 得到 − k):

注意(27)式前的正号, 其对应于弱局域化效应, 即产生正的磁导, 这是双外尔半金属弱局域化的另一个特征.在低温极限下 lϕ≫lB≫lz, 磁导(B)∝, 在极限 lB≫lϕ以及 lB≫lz时,

根据文献[73,75]的理论结果, 我们提出了一个公式来拟合3D 弱(反)局域化引起的磁导:

存在相互作用的情况下, 对于外尔费米子的一个能谷, 电导率随温度的变化可以概括为

其中 cee和 cqi都是正的参数.(29)式描述了由相互作用引起的弱局域化和由干涉引起的弱反局域化之间的竞争(图8).临界温度 Tc=(cee/p·cqi)2/(p−1)因为 ceecqi＞0 , 这意味着只要 p ＞1 , 总有一个临界温度低于这个温度, 电导率就会随着温度的降低而下降.对于3D 中已知的退相干机制, p 总是大于1[191].通过一组典型参数的拟合, 发现 Tc≈0.4— 1 06K[73].

图8 电导率 Δ σ 随温度 T 变化的示意图.选择 c ee =cqi ,Tc是电导率随温度降低而下降的临界温度.转载自文献 [73]Fig.8.The schematic diagram of conductivity Δ σ changes with temperature T.We choose c ee =cqi , Tc is the critical temperature at which the conductivity drops with temperature.Reproduced with permission from Ref.[73].

能谷间的散射和关联也会导致弱局域化, 如图6(b)所示, 随着 ηI的增加, 负的 δ σqi被抑制.此外, 图6(c)显示, 当 ηI+η∗=3/2 时, 磁导可以变为正, 即局域化趋势.因此, 强的能谷间散射和关联的结合将增强无序外尔半金属的局域化倾向并可能导致金属-绝缘体相变[39].

4.3 节线半金属中的弱(反)局域化

模型(15)具有两个基本的对称性[147], 即=H(−k) 和=H(−δk) , 其中T1=K , T2=iτyK 为时间反演, k 是复共轭算符,δk=κ(cos φ,sin φ/α).由表1 可知, T1和 T2分别对应于正交和辛对称类, 分别导致弱局域化和弱反局域化(图9(b)—图9(c)).简单起见, 考虑低能激发即 EF相对较小, 费米面的主半径满足k0≫1/rsc≫κ , rsc为散射势的范围.此时不同的散射态k之间可以用一个单一的参数 θ 来描述, 即uk=u(θ)=u0fΔ(θ) , 其中 fΔ(x)=Θ(x+Δ)Θ(−x+Δ) ,Θ(x) 为阶跃函数.当 Δ =π 时为短程散射极限,Δ →0时为长程散射极限.

与第4.1 节类似, 在文献[147]中我们利用费曼图方法计算了电导率的量子干涉修正, 计算中主要包括了三个领头阶的图[195,196].对于零温电导率,总的量子修正为

图9 (a)节线半金属的轮胎状费米面, 小半径 κ , 主半径k0, 极向角 φ , 环面角 θ ; (b)对于短程杂质势导致在环形方向上产生从 k 到 − k 的相干背散射; (c)在长程杂质势作用下, 沿极向的 δ k 到 − δk 的散射, 此过程积累一个大小为π的贝里相位.转载自文献[147]Fig.9.(a) Torus-shaped Fermi surface of nodal-line semimetals, with minor radius κ , major radius k0 , poloidal angle φ, and toroidal angle θ ; (b) a coherent backscattering from wave vector k to − k around the toroidal direction for shortranged impurity potentials; (c) backscattering from wave vector δ k to − δk along the poloidal direction under long-ranged impurity potentials.The process contributs a π Berry phase.Reproduced with permission from Ref.[147].

其中 s =2 为自旋简并度, ηz=2 是对速 度的 顶点修正.推迟(R)和超前(A)格林函数利用一阶波恩近似求解得到, 即是弹性散射的弛豫时间, 其中=是最大交叉图的Cooperon[204].在短程散射极限下 Δ =π , 最主要的散射路径为(k1,k2,··· ,kn) 及其时间反演(−kn,··· ,−k2,−k1)(图9(b)), 需要指出的是, 这些路径没有包围节线,即不会积累贝里相位.直流极限下(ω →0)[147]

长程极限时 Δ →0 , 量子干涉修正主要由路径(δk1,δk2,··· ,δkn) 和 (−δkn,··· ,−δk2,−δk1) 所贡献(图9(c)), 此路径由于包围节线会积累一个大小为π的贝里相位.类似地可以得到[147]

引入磁场后, 退相干机制抑制了无序诱导的电导率量子修正[204].在短程极限下得到磁导率为

类似地在长程极限下

其中 Ψ (x) 为双伽马函数.图10 给出了(33)式和(34)式对磁导率的理论预测.短程极限下弱局域化效应导致正的磁导(图10(a)), 而在长程极限下的弱反局域化效应导致了负的磁导(图10(b)).实验中可以通过对温度的调节来控制相位相干长度lϕ的大小以匹配不同的曲线, 在低温及 B ∼0.1 —1.0 T时, 由于 lϕ∼100 nm—1 µm, lB∼10 nm, 即 lϕ≫lB.通常, 与3D 扩散相比, 2D 扩散会导致更大的弱局域化或弱反局域化效应.此外, 节线半金属的3D 扩散行为可以视为大量的2D 有效子系统, 导致磁导率的显著提高[205].实际上, 在图10 中可以观察到,大3 个数量级, 这表明弱反局域化效应在节线半金属的长程散射极限下是个非常强的特征.

图10 不同相位相干长度 lϕ 下的磁导率, 短程极限(a)和长程极限(b)分别对应(33)式和(34)式.转载自文献[147]Fig.10.The magnetoconductivity in the (a) short range limit Eq.(33) and (b) long range limit Eq.(34) for different phase coherence lengths lϕ.Reproduced with permission from Ref.[147].

4.4 拓扑绝缘体中的弱(反)局域化

基于费曼图方法我们计算了拓扑绝缘体磁掺杂表面态的磁导[71], 磁掺杂的引入以两种方式影响了体系的性质, 首先磁性杂质导致的平均场在狄拉克点处打开能隙, 其次平均场的局部涨落以随机方式散射表面态电子.经过计算发现, 除了无能隙狄拉克费米子产生的弱反局域化效应, 能隙的打开导致了额外的弱局域化项, 随着能隙的进一步增加将使体系展现出完全的弱局域化.

零温磁导公式为

其中Ψ是双伽马函数,

对非磁性掺杂的拓扑绝缘体的电子输运实验展现出非常大的矛盾[206−209], 磁导率随磁场的增加而下降, 即展现出反局域化的行为(图11(a)红色曲线).然而, 当温度降低时, 测得的电导率呈对数下降, 显示出普通无序金属中弱局域化的典型特征(图11(b)蓝色曲线).为了引入温度的效应, 我们考虑电子-电子相互作用对电导率的修正[72]:

图11 (a)磁导率 δ σ ≡σ(B)−σ(0) 与磁场 B 的关系; (b) 电导率 σ 与温度 T 的关系.转载自文献[72]Fig.11.The magnetoconductivity δ σ ≡σ(B)−σ(0) ;(b) conductivity σ vs temperature T.Reproduced with permission from Ref.[72].

ε0是真空介电常数, εr是考虑晶格离子和价电子影响的相对介电常数,.电导率公式是温度T和磁场 B 的函数, 并且依赖于狄拉克模型参数Δ/(2EF) 和 γ , 以及与样品相关的参数 le,lϕ,F.(38)式明确阐明了由相互作用主导的电导率的温度依赖性, 而磁导率主要由量子干涉贡献, 解决了拓扑绝缘体输运实验上的矛盾.

5 弱磁场

5.1 拓扑半金属中的负磁阻和平面霍尔效应

在拓扑半金属中非平庸的贝里曲率可以将外部磁场与电子速度耦合, 产生随着磁场的增加而增长的额外导电性[28,29].由于它与成对外尔点之间的手征电荷转移相关, 负磁阻也被认为是手征反常的一个特征[26,210,211].许多拓扑半金属材料中都观察到了负磁阻现象[41−43,46,55−57,60,61,63,66,68].除了在纵向磁场下产生随磁场增大的手征电流以外, 最近的理论发现, 外尔半金属中可以存在平面霍尔效应[212,213], 即磁场平行于电场时横向有电荷积累即产生霍尔电导率 σyx.不同于经典的霍尔效应, 平面霍尔效应不是由洛伦兹力导致的, 因此吸引了很多理论和实验的关注[214−216].文献[217]研究了普通外尔半金属和倾斜外尔半金属中的平面霍尔效应, 相比于前者, 倾斜外尔半金属的纵向电导率σxx和霍尔电导率 σyx均给出了不同的磁场 B 依赖关系.

我们从电子的半经典运动方程开始[95,218−220]:

其中 Dk=1+(e/ℏ)B·Ωk是对相空间体积的修正项, 包含 Ω 的各项分别导致了反常霍尔效应[221−223]、手征磁效应[36]、以及负磁阻[28,29].

由于电导率是流-流关联(图5(a)), 速度与B的线性相关性导致电导率的 B2相关性.在2.1 节,已经证明贝里曲率与 1 /k2成正比.考虑到电导率公式的3D 积分中存在 Ω2以及 k2, 则反常电导率部分应与费米波矢成反比, 与 B2成正比, 即

(41)式与文献[28,29]中得到的公式是一致的.电导率随着 B2的增加而增加, 产生负磁阻.由于非平庸的贝里曲率在外尔点处发散, 电导率随费米波矢和载流子密度的减小而增大.3D 时, 载流子密度n与成正比, 即

因此, 为了验证非平庸贝里曲率的负磁阻特性, 有必要对三个性质进行检验, 即角度依赖性、B2的磁场依赖性以及 n−2/3的载流子密度依赖性.在实验[61]中, 通过比较不同样品的结果检验了磁阻对载流子密度的依赖性.

在半经典区域, 电流密度矢量为[95]

括号内第一项表示电场 E 沿着 v 的分量引起的载流子分布偏离, 第二项表示手征化学势导致的偏离.

对于由最小有效模型(7)描述的外尔半金属,由于手征化学势的影响, 电子通过电场和磁场的共同作用从一个外尔锥转移到另一个手性相反的锥.如(44)式所示, 偏离平衡分布的这部分与 B ·E 成正比, 即当 B 和 E 不垂直时, 与磁场的大小成正比.在没有倾斜的情况下, 外尔锥在空间上是均匀的,并且速度 v 的大小在费米面上没有变化, 因此 v 在单个外尔锥的贡献为零.而且由于电子分布的偏离, 与磁场 B 成正比的反常速度项 (Ωk·v)B 对两个外尔锥的影响不能相互抵消.手征化学势和反常速度的共同作用导致了 B2依赖的电流密度矢量(43)式, 其纵向分量导致了磁场平方依赖的纵向磁导率 σxx, 即手征反常.它的横向分量产生了磁场平方依赖的霍尔电导率 σyx, 即平面霍尔效应.

对于倾斜的外尔半金属, 其中一个外尔点上的最小有效模型可以写为[217]

其中参数 t 为描述倾斜程度的矢量, 当|t|＜v 时该模型描述Ⅰ型外尔半金属, 以下讨论主要基于此类体系.

引入倾斜项后两个具有相反手性的外尔锥沿着相反的方向发生倾斜.假设 E 沿着倾斜的方向,(44)式括号中第一项显示了由电场 E 导致的分布偏离.因此, 两个外尔锥中占据态的反常速度的贡献不会抵消.但这一项不依赖于磁场 B , 只有反常速度项 (Ωk·v)B 线性依赖于磁场 B , 由此导致的纵向和横向电流分别产生了线性依赖于磁场的σxx和 σyx.考虑到(44)式括号中第二项手征化学势的影响后, 不同于普通外尔半金属, 此时 v 的大小在费米面上是变化的, 即 v 在单个外尔锥上的贡献不会相互抵消.只要 B 和 E 不垂直, 则手征化学势导致的载流子分布偏离正比于磁场的大小.导致除了由反常速度引起的2 次依赖于磁场的项外进一步产生了线性依赖于磁场的 σxx和 σyx.

文献[217]进一步分析了磁场的角度和相空间体积的修正对电子输运的影响, 并指出该半经典理论可以推广到Ⅱ型外尔半金属, 为未来可能的实验验证提供指导.

5.2 拓扑绝缘体中的负磁阻

拓扑半金属中的手征反常[26,210,211]被广泛认为是产生负磁阻的原因[41−43,46,55,56,60,61,63,65,66,68,224].然而, 在无法定义手征的拓扑绝缘体中, 也可以观察到负磁阻现象, 这导致了对负磁阻进行理论解释的巨大混乱[152−155,225−228].在文献[80]中, 我们使用具有贝里曲率和轨道磁矩修正后的半经典玻尔兹曼方程来解释拓扑绝缘体中的负磁阻, 并与实验结果进行了定量的比较(图12).

实验中, 负磁阻存在于 T =100 K[154]的温度以上, 因此可以排除量子干涉机制.此外, 由于拓扑绝缘体Bi2Te3和Bi2Se3[229]的迁移率较差, 当磁场达到6 T 时, 朗道能级不能很好地形成.半经典运动方程(39)的速度项为

图12 理论计算的负磁阻与实验[152−154]的比较.转载自文献[80]Fig.12.The comparison between the theoretical negative magnetoresistance and the experiments[152−154].Reproduced with permission from Ref.[80].

在玻尔兹曼理论中, 纵向电导率 σµµ由穿过费米能级的所有能带决定, 对于能带 n[223],

其中 feq是平衡态的费米分布, 弛豫时间 τ 假定为常数[29].

图12 表明了负磁阻的实验和数值计算结果之间的一致性.在实验[152]中, 温度为1.8 K, 因此原始数据在零场附近具有正磁阻, 这是由弱反局域化导致的[206,208,230−233].在实验[153,154]中负磁阻随温度变化不大这与理论保持一致, 显示了负磁阻的半经典性质.在大多数使用半经典输运理论研究磁阻的文献中, 轨道磁矩均被忽略[28,223].如文献[234,235],轨道磁矩对于外尔半金属中磁阻的各向异性是必不可少的.文献[80]表明, m ·B 对 ε (k) 的修正可以有效地增加负磁阻.

其他的一些工作工作[236,237]发展了一般3D 非磁性金属的具有二阶精度的半经典理论.这些工作给出了几个令人惊讶的结果.首先, 由于 δ σint/σ0与弛豫时间无关, 存在本征磁导 δ σint.第二, 显著的δσint项可能导致对科勒规则的偏离.第三, δ σint可能导致正的纵向磁导.在仅占据最低朗道能带的量子极限下, 磁阻微妙地依赖于散射机制[74,76,238], 而不是贝里曲率和轨道磁矩.最近的理论工作[227,239]为理解负磁阻提供了新的观点.值得注意的是文献[240]利用格林函数方法系统地计算了有质量狄拉克费米子在磁场下的输运行为, 该理论给出了负的且正比于 B2的纵向磁阻, 这一结论有望澄清长久以来对拓扑物态中负磁阻解释的困难.

6 强磁场

6.1 外尔和狄拉克半金属中的量子振荡

由2.1 节可知, 在 z 方向的磁场中能谱演化为一系列1D 的朗道能带[74,76](图2), 这导致了电阻的SdH 振荡.一般用栗弗席兹-科塞维奇(Lifshitz–Kosevich)公式[241]来描述:

式中, ϕ 为相移, F 为振荡频率, B 为磁场大小.每个频率分量的相移为 − 1/2+ϕB/2π+ϕ3D, 其中ϕ3D=∓1/8 为仅针对3D 的修正, ϕB为贝里相位[95,242].费米面沿磁场方向的曲率决定了 ϕ3D的符号[243,244].拓扑物态为探索非平庸贝里相位提供了一个新的平台[50−52,54,67,245−255].拓扑半金属中具有额外大小为 π 的贝里相位[122,242], 因此在2D 和3D 体系中相移分别为 0[256]和 ± 1/8[257].在文献[77]中, 我们发现在栗弗席兹点附近, 量子振荡的相移可以超过一般的 ± 1/8 或 ± 5/8 , 并且非单调地向更宽范围 ± 7/8 和 ± 9/8 之间移动.到目前为止, 在HfSiS[258−264], ZrSi(Se/Te)[265]和ZrGe(S/Se/Te)[266]中对量子振荡都进行了实验研究, 但它们的相移并不相同.表2 总结了具有线性和抛物线色散的 2 维和 3 维能带的相移.

表2 对于具有不同色散和维度的系统, (48)式中的相移 ϕ.B z 和 B // 是节线平面内外的磁场.α , β , γ , δ 对应于图14 中费米面的截面.转载自文献[79]Table 2.For systems with different dispersion and dimensions, the phase shift φ in Eq.(48).B z and B // are magnetic fields outside and inside the nodal-line plane.α , β ,γ , and δ correspond to the cross sections of Fermi surface in Fig.14.Reproduced with permission from Ref.[79].

实验中, 在 B 轴上的峰值位置或谷值位置给出整数朗道指数 n , 进而从 n 和 1 /B 的图中得到相移和频率 F.尽管如此, 对于是峰值[50−52,66,67,246,257]还是谷值[54,248,251,267]应该被赋予朗道指数尚不明晰.我们的计算结果表明, 在朗道能带边缘出现了与整数朗道指数对应的电阻率 ρzz和 ρxx的峰值.利用线性响应理论[268−271]估算电阻率分量[272,273].对于纵向情况, 电阻率 ρzz= 1 /σzz.在带边附近, 电导率 σzz呈现谷值, 因此 ρzz呈现峰值.对于横向情况,, 纵向电导率和场致霍尔电导率分别为

式中 δ ≪1 为振荡部分, σ0为零场电导率.通过详细分析[74,269,274]可知, 峰的位置遵循关系式ρzz∼, 因 此 ρzz和 ρxx在朗道能带边缘出现的峰及其相移是相同的.

由图13 可见, 数值结果与理论分析预测高度一致.定义 EA=Akw和对于随着拍频模式的出现, ϕ -EF曲线开始出现断裂.在图13(c)中, 对于 EA＜EM相移降到 − 5/8 以下, 而不是在栗弗席兹点周围从–1/8 单调地移到–5/8(即EF=EM), 这是因为没有简单的依赖关系[77].在栗弗席兹点, 可以解析地证明相移为 − 9/8 , 这与图13(c)中的相移一致.它相当于 − 1/8 , 通常认为它起源于大小为 π 的贝里相位.

由(11)式描述的狄拉克半金属具有时间反演对称性[105,275]].对于狄拉克半金属, 总相移可以取两个值, α ∈[0,1/4] 和 [3/4,1] 时为 − 1/8 , α ∈[1/4,3/4] 为 − 5/8.在栗弗席兹点附近, 总相移可能在两个值之间变化.

一般情况下电子载流子应该产生负相移, 空穴载流子应该产生正相移[257].然而, 在狄拉克半金属Cd3As2的实验中, 电子载流子的相移为正值[50,51,54].一种可能的解释是, 对于实验中的相移1/8—3/8,由于 2π 的周期性它们的实际值应为–7/8—–5/8.表3 列出了相移实验值的对应项.在文献[77]中,除了轨道量子干涉[276]、塞曼劈裂[52,247,277]和费米面嵌套[51]等机制, 我们还证明了由于能带反转而出现的拍频模式.

表3 从Cd3As2 的实验中得到的相移 ϕ exp.转载自文献[77]Table 3.The phase shift ϕ exp obtained from the experiment of Cd3As2.Reproduced with permission from Ref.[77].

6.2 节线半金属中的量子振荡

对于节线半金属无论是面内还是面外, 轮胎状费米面都具有最大截面和最小截面.其次, 贝里相位沿与节线平行的圆为0, 沿包围节线的圆为 π.根据昂萨格关系得到, 其中 A 是费米面上垂直于磁场的极值截面面积.当节线平面垂直于磁场, 即磁场取 Bz时, kz=0 平面上有两个极值截面(图14(b)).通过理论分析发现, 外圆的高频项为, 内圆的低频项为Fβ=m(u−EF)/(ℏe).这两个频率可能导致拍频模式的出现[278].在文献[79]中, 通过分析电阻率的计算结果得到了节线半金属的相移和频率.需要指出的是, 当保护节线的对称性破缺时, 出现大小为Δ的能隙将导带和价带分开, 此时贝里相位为表4 进一步总结了随机情况下相移的一般规则, 这些一般的规则可以帮助分析一些其他的材料, 例如ZrSiS 和Cu3PdN[259−264].

图13 对于(1)式描述的外尔半金属, 数值(散点)和解析(实线)得到的频率 F 的曲线 (a)固定 M 对应不同的 A ; (b)固定A对应不同的 M.(c)固定 E M 不同的 E A 对应的相移 ϕ 的曲线.曲线断裂是因为在拍频模式出现时, F 和 ϕ 无法拟合.垂直虚线表示栗弗席兹点.转载自文献[77]Fig.13.For the Weyl semimetal described in Eq.(1), the frequency F obtained by numerical (scatters) and analytical (solid curves): (a) Fixed M corresponds to different A ; (b) fixed A corresponds to different M.(c) Fixed E M , for different E A corresponds to the curve of phase shift ϕ.The curve breaks because F and ϕ can not fit when beating patterns occur.The vertical dashed lines represents the Lifshitz point.Reproduced with permission from Ref.[77].

对于节线半金属, 大多数量子振荡实验都是针对ZrSiS 族材料[258−266], 其中费米能处既有电子型也有空穴型口袋[140,141,259,264].文献[79]进一步分析了Cu3PdN[120,121]的相移.

图14 (a) (17)式中节线半金属的节线(虚线环), 轮胎状和鼓形费米面, E F 是费米能, u 是模型参数; (b)轮胎状费米面在节线平面内的最大(α)和最小(β)截面; (c)轮胎状费米面在节线平面外的最大(γ)和最小(δ)截面.转载自文献[79]Fig.14.(a) In the model of nodal-line semimetal Eq.(17),the nodal line (dashed ring), torus and drum Fermi surface,EF is Fermi energy, u is model parameter; (b) the maximum(α) and minimum (β) cross sections of the torus Fermi surface; (c) the maximum (γ) and minimum (δ) cross sections of the Fermi surface outside the nodal-line plane.Reproduced with permission from Ref.[79].

表4 节线半金属的相移 ϕ.α ,β,γ,δ 是图14 中的极值截面.转载自文献[79]Table 4.The phase shift ϕ of the nodal-line semimetal.α ,β,γ,δ are the extremal cross sections in Fig.14.Reproduced with permission from Ref.[79].

6.3 3D 量子霍尔效应

2D 量子霍尔效应的发现打开了拓扑物态领域的大门[157,279].在3D 电子气中, 由于沿磁场方向的波矢 k 是好量子数, 量子霍尔效应通常只能在2D系统中观察到[157,256,280−282].在文献[78,283]中, 我们描述了拓扑半金属中的3D 量子霍尔效应.拓扑半金属在平行于外尔点方向的表面处有拓扑保护的表面态, 即费米弧[284,285−288](图15(a)和图15(b)).

图15 (a)外尔半金属中的费米弧和体态的色散, k // 表示 (kx,ky) ; (b)在 k z-kx 平面上 y =L/2 , E F =Ew 处的费米弧; (c)宽度为 W , 厚度为 L 的外尔半金属板; (d)在EF = 处的费米弧(实线); 波函数在 处沿 y 轴的分布; (h) 3D 量子霍尔效应中的朗道能级和边缘态;(i)单一表面的电子无法被y 方向的磁场 B 驱动完成一个完整的回旋运动.转载自文献[78,283]Fig.15.(a) The energy dispersions of the Fermi arc and bulk states in a Weyl semimetal, k // stands for (kx,ky) ;(b) the Fermi arc at y =L/2 and E F =Ew in the kz-kx plane; (c) a Weyl semimetal slab with width W and thickness L; (d) Fermi arc (solid) at the distribution of wave function along y -axis at k z =0 ; (h) the Landau levels and edge states in the 3D quantum Hall effect; (i) an electron in single surface could not be driven in a y-direction magnetic field B to perform a complete cyclotron motion.Reproduced with permission from Refs.[78,283].

单一表面上, 费米弧无法形成一个闭合的费米环, 而这一点对量子霍尔效应至关重要.然而, 在外尔半金属中, 来自相反表面的费米弧(图15(c))可以形成所需的闭合费米环(图15(d)).电子通过外尔点在相对表面的费米弧之间实现隧穿(图15(e)—图15(g)), 通过这一机制电子可以完成完整的回旋运动, 从而实现3D 量子霍尔效应.在实际材料中,隧穿距离受平均自由程的限制, 在高迁移率的拓扑半金属[60,103]中, 平均自由程大约为 1 00 nm , 在一些材料中甚至可以达到 1 µm[289], 在我们的计算中块体厚度为100 nm.为使量子霍尔效应仅来自费米弧的贡献需将费米能级调到外尔点上来耗尽体态载流子[290].外尔半金属TaAs 家族[291−295]和狄拉克半金属Cd3As2和Na3Bi 具有量子霍尔效应所要求的极高迁移率[46,48,50,51,54]、低载流子密度[224]等特性.通过费米弧的3D 量子霍尔效应预计可以存在于Cd3As2[224,296−298], Na3Bi 以及TaAs家族的材料中[61,290].

根据久保(Kubo)公式, 我们计算了霍尔电导[78,299−303].当费米能远离外尔点时, 霍尔电导服从通常的 1 /B 依赖关系.当费米能向外尔点移动时, 开始出现的量子化平台.

从图15(c)可见, 费米弧的边界态具有独特的3D 空间分布.具体而言, 上边界态向左传播(绿色箭头), 下边界态向右传播(橙色箭头).这种费米弧边缘态的独特3D 分布可以通过扫描隧道显微镜[304]或微波阻抗显微镜[305]来探测.与拓扑绝缘体[281,282]不同, 基于费米弧的量子霍尔效应需要两个表面的共同作用.最近的工作[306]系统地探讨了外尔半金属中3D 量子霍尔效应的边缘态图像.外尔半金属中不同手性外尔点处的速度相反, 强磁场使体态进入手征朗道能带, 这一点显著地影响了边缘态的载流子输运.除了上下表面的费米弧通过隧穿形成一个闭环[78]外, 在磁场作用下拓扑半金属中还观察到其他奇异的量子霍尔现象[296,307,308].在垂直磁场 B 中, 基于半经典运动方程从理论上给出了边缘态的闭合3D 轨道.在 x 方向, 由于手征朗道能带的影响, 上(下)表面费米弧态的半圆轨道形成跳跃边缘态.在 z 方向, 手征朗道能带和侧面的费米弧态共同构成闭合轨道.此外, 由于手征朗道能带沿磁场方向色散, 边缘态的量子通道可以由倾斜的磁场控制, 导致边缘态的分布和霍尔电导的行为发生了显著的变化.当调整 Bz并保持 By不变时, 边缘态的位置会从一侧变化到另一侧, 霍尔电导也随之改变符号.文献[306]利用数值方法进一步分析了局域态密度分布, 以及霍尔电导随磁场大小和方向的变化.

3D 量子霍尔效应也可以由电荷密度波(CDW)机制[309]实现, 最近在ZrTe5[310]的实验中观察到了这一机制, 并看到了清晰的量子化平台.在文献[87]中我们详细地探讨了这一机制, 并给出了一个完备的理论可以解释实验的主要特征.我们发现电荷密度波可以在1D 的朗道能带上产生, 而且强烈依赖于磁场的大小.磁场诱导了3D 电子气向电荷密度波相的二阶相变, 导致朗道能带打开一个能隙从而使体态绝缘, 即3D 电子气变为一系列的2D 量子霍尔态, 从而实现3D 量子霍尔效应（图16）.

通过沿 z 方向的电子-电子或电子-声子相互作用, 可使 kF和 − kF附近的电子之间产生耦合从而打开电荷密度波能隙(由序参量 Δ 描述).电子-电子相互作用为[311,312]

图16 左图: 2D 电子气在磁场中形成量子霍尔态.中间图: 3D 时朗道能级变为一系列2D 的朗道能带.右图:电荷密度波使朗道能带打开能隙, 使体态绝缘, 可以观察到3D 量子霍尔效应.转载自文献[87]Fig.16.Left: the quantum Hall state in 2D electron gas under magnetic field.Center: in 3D, the Landau levels turn to one dimensional Landau bands.Right: the charge density wave gap the Landau band, so that the bulk is insulating and the 3D quantum Hall effect can be observed.Reproduced with permission from Ref.[87].

在实验中, 霍尔电阻率的平台覆盖了1.7—2.1 T 的较宽范围, 以 (e2/h) 为单位的霍尔电导率由电荷密度波的层数决定.电荷密度波波长 λcdw与费米波长相关:

在1.7—2.1 T 之间观察到的 ρxx平台意味着存在公度的电荷密度波, 即电荷密度波波长与晶格常数a的整数倍相称.通过比较公度和非公度的电荷密度波在2.1 T 附近的基态能量, 发现对于公度和非公度的电荷密度波在磁场分别在 [1.7,2.1] T 以及[2.1,3.0] T 的区间时具有较低的能量, 因此在两种电荷密度波之间存在交叉.在 B ∈[1.7,2.1] T 范围内固定的 λcdw表示费米能不变, 即系统属于巨正则系综载流子的数量可以改变.相比之下, 非公度电荷密度波中的载流子数量不能改变.因此, 电子能量的变化导致 B ∈[1.7,2.1] T 范围内公度的电荷密度波基态能量降低.进一步将磁场增加到2.1 T 以上, 磁场将把费米能推得更低(最终到达带底), 进一步说明由磁场诱导的电荷密度波存在从公度到非公度的交叉.

7 极强磁场

7.1 量子极限下的负磁阻

选择随机高斯势来探讨量子极限下无序散射对磁阻的影响

式中, ui表示随机分布杂质在 Ri处的散射强度, d是确定散射势范围的参数.杂质势范围 d 和磁长度lB定义了两个区域, 长程势区域 d ≫lB和短程势极限 d ≪lB.在强磁场极限下, 如果 B ＞10 T, 磁长度 lB小于 1 0 nm.

通过考虑散射时间对磁场的依赖性, 我们发现在强场极限(B →∞)下, 即对于第0 个朗道能带,费米面上的态间散射的弛豫时间为得到强场极限下的电导率

当磁场为零时, 磁长度发散, d /lB→0 , 给出最小电导率[76]

注意外尔点处态密度在零磁场时趋于0, 类似的结果在没有朗道能级的情况下也曾被发现过[316].通过进一步的计算得到磁导[76]

图17 给出了不同极限下的纵向和横向磁导.特别地, 图17(f)中长程势极限下的 σxx∝1/B.垂直于 x -y 平面的场中, 还存在霍尔电导率σyx=sgn(M)(kw/π)e2/h+en0/B, 其中第一项是反常霍尔电导率, 第二项是经典电导率.在弱场中, 经典霍尔效应占主导地位, σxx和 σyx都与 1 /B 成 正比,在磁场 B 中电阻是线性的.注意, 这里垂直场中的线性磁阻与以前的机制并不相同[274,317].

7.2 量子极限下的背散射禁止

拓扑绝缘体作为拓扑物态中最先被关注的体系, 为探索奇异的拓扑物态[318−325]提供了诸多线索.在强磁场中, 2D 拓扑绝缘体的最低朗道能级相互交叉, 该特征可以作为量子自旋霍尔相的标志[326,327].同样在2D 情况下, 干涉效应[328]也可以用来探测量子自旋霍尔相.然而在3D 时, 如何利用最低朗道能带来识别拓扑绝缘体却很少被讨论.在文献[81]中, 我们研究了3D 拓扑绝缘体在强场量子极限下的磁阻(图18(b)).在临界磁场中, 背散射在量子极限下被完全抑制, 这一效应可以用来识别拓扑绝缘体.不过这种禁止的背散射在拓扑半金属中是不存在的[74,76,238].这一理论与最近的实验展现出惊人的一致性(图18(c)和图18(d)).此外, 这一机制可以很好地应用于那些能隙较小的材料, 例如ZrTe5[329,330]和Ag2Te[331].

在沿 z 方向的强磁场 B 中, 能谱量子化为一系列1D 的朗道能带(图18(a)和图18(b)).最低的两个朗道能带分别表示为 0 + 和 0 − , 其能量为E0±=质量项m=当费米能仅与 0 + 朗道能带相交, 其本征态为

图17 不同的势范围下, 外尔半金属在 方向磁场B 中的纵向电导率 σ zz 和横向电导率 σ xx.转载自文献[76]Fig.17.The longitudinal conductivity σ zz and transverse conductivity σ xx of the Weyl semimetal in the -direction magnetic field B under the different potential ranges.Reproduced with permission from Ref.[76].

图18 (a) 3D 拓扑绝缘体的零场能谱(k x =ky =0); (b)在强磁场中, 费米能只穿过 0 + 朗道能带; (c)实验测得Pb1–xSnxSe 的磁阻[226]; (d)理论计算出的磁阻.转载自文献[81]Fig.18.(a) The zero field energy spectrum of 3D topological insulator (k x =ky =0); (b) in a strong magnetic field, fermi energy EF can only crosses the 0 + Landau energy band; (c) the magnetoresistance of Pb1–xSnxSe in experiment[226]; (d) the theoretical calculated magnetoresistance.Reproduced with permission from Ref.[81].

7.3 外尔费米子湮灭

当外界强磁场与外尔费米子耦合时其手性消失并获得质量, 而这是外尔费米子最明显的两个特征, 在这个意义上, 外尔费米子被湮灭.在外尔半金属TaP 中我们实现了这种磁耦合[64].

2.1 节给出了沿 z 方向磁场中的朗道能带.现在把这种情况推广到垂直于 y 方向B =B(sin ϕ,0,cos ϕ) 的 任意场, 其中 ϕ 是 z 方向和磁场之间的夹角.朗道规范可选为 A =(−Bzy,0,Bxy) ,在皮尔斯(Pierls)代换下

哈密顿量变为

图19 给出了外尔半金属的朗道能带.当磁场沿 z 方向(θ =0)时, 最低的朗道能带(红色)穿过费米能(虚线).随着磁场从 z 方向旋转到 x 方向(θ =π/2), 最低朗道能带发生了移动.当磁场沿x方向时, 朗道能带的能谱是粒子-空穴对称的, 并且由于外尔费米子之间的耦合存在一个能隙.这就是为什么在外尔半金属TaP[64]的强场量子极限下,霍尔电阻出现明显的符号反转.由于能隙的存在,外尔费米子获得质量并失去其手性.由于具有手性和无质量是外尔费米子的两个主要特征, 因此霍尔信号表明外尔费米子被湮灭.

7.4 外尔半金属中的不饱和磁化

强磁场中外尔费米子的第0 个朗道带是纯手征的[27,210,211].然而, 量子极限下的的磁阻复杂地依赖于杂质散射的性质, 不能确定地给出准粒子的能谱信息[74,76,77,238].相比之下, 电子的磁响应要简单得多, 因为它们不会与杂质散射相互作用.事实上,它们仅由电子的热力学势 Ω 对于磁场 H 的导数决定, 而且已被证明可以用来探测费米面的性质[241,332,333].在文献[85]中我们提出了一种在量子极限下利用磁响应探测相对论性准粒子的方法.我们证明外尔半金属TaAs 的有效横向磁化强度MT和有效平行磁化强度 M//是超越量子极限的准线性场依赖的, 这些非饱和磁响应明显不同于非相对论性准粒子.该方法作为一种热力学探测手段, 可用于识别拓扑材料中的相对论性准粒子.

图19 外尔半金属在垂直 y 方向的磁场 B 中的朗道能带,其中 k // ≡kx sin ϕ+kz cos ϕ 为平行于 B 的波矢, tan ϕ=Bx/Bz.红色曲线是第0 个朗道能带, 虚线是费米能.转载自文献[70]Fig.19.The Landau energy band of Weyl semimetal in the magnetic field B perpendicular to the y direction, where k// ≡kx sin ϕ+kz cos ϕ is the wave vector parallel to B ,tan ϕ=Bx/Bz.The red curve is the 0 th Landau energy band, and the dashed line represents the Fermi energy.Reproduced with permission from Ref.[70].

磁化强度来自于导带和价带的共同贡献.当磁场增加时, 朗道能带相继离开费米面, 在只有最低的朗道能带被占据之前 M//在零附近振荡 MT随着磁场的增加而下降.在磁场较弱时有很多朗道能带被占据, 费米能 EF几乎保持不变, 因此对于固定费米能 EF和固定载流子浓度 Nc两种情况弱场的计算结果几乎相同.然而, 当系统进入量子极限时,EF在强磁场中的变化将不再可以忽略.如果 EF固定, 零点能将使 EF之上的第0 个朗道能带继续升高以致超过临界磁场, 导致 M//和 MT消失.如果Nc固定, 则无论磁场增加多少, EF都将被固定在第0 个朗道能带的边缘, 导致 M//出现饱和, MT在一个常值上保持不变.这种强场行为可以归结为价带和导带的贡献在强场下饱和, 换句话说, 由于非相对论性的能带色散, M//和 MT在量子极限下是不变的.对于铋、硫掺杂的Bi2Te3和InSb 这种3D 有能隙的体系, 实验观察到了 M//和 MT在量子极限下明显的饱和行为[333,334].

相比之下, 由于能带反转和相对论性能谱, 无能隙的外尔费米子以一种不同的方式对磁化强度产生贡献.在弱场中, M//和 MT都围绕一个平均值振荡.当除了第0 个以外的所有朗道能带都离开EF时, 两种手征模式都对磁响应有贡献, 导致强磁场中不饱和的 M//和 MT.

先前的理论[335−337]只针对不同类型的能带结构在全磁场范围内计算了 M//, 这里为任意角度磁场下的 M//, MT建立一个更一般的磁化理论.详细的计算参考文献[85], 这里给出强场极限下的主要结论.由方程(2) 描述的单节点3D 无质量外尔费米子的哈密顿量可以写为

以前的研究主要集中于弱磁场下的磁化, 如文献[338]所指出的外尔半金属的贝里顺磁性.此外, 最近对晶格模型的理论计算表明, 轨道对磁化强度大小和符号的贡献取决于外尔/狄拉克半金属中费米能的精确值[339].对于模型(59)量子极限下的平行磁化强度为

横向磁化强度为

图20 给出了在强磁场中对TaAs 样品进行测量得到的 M//和 MT.由图可见实验数据可以很好地用(60)式和(61)式进行拟合.我们对3D 外尔费米子的理论计算再现了弱场中的德哈斯-范阿尔芬(dHvA)振荡以及量子极限下 M//显著增强的行为.需要强调的是, 不同于非相对论性电子的饱和磁化, 对于固定 EF和固定 Nc两种情况 M//和 MT在强磁场中均给出了类似的趋势.

图20 实验与理论的比较 (a) M // (黑线); (b) M T (黑线).转载自文献[85]Fig.20.Comparison between experiments and theory for the M // and M T of TaAs: (a) M // (black line); (b) MT(black line).Reproduced with permission from Ref.[85].

7.5 反常热电效应

热电效应是对材料拓扑性质更敏感的一种探测手段.非平庸电子的存在导致拓扑材料呈现出反常的热电行为.在较高的磁场下, 能斯特信号在热电势接近零的临界磁场处出现符号反转.在强场量子极限下, 热电势和能斯特信号都呈现出奇异的峰值.我们认为这些反常行为可以归因于拓扑材料中第0 个朗道能带的能隙闭合和能带反转[86].

过渡金属五碲化物(ZrTe5, HfTe5等)作为非常接近拓扑相变临界的拓扑材料引起了相当大的关注[43−45,340−347].近40 年来, 作为热电材料的ZrTe5以其巨大的热电势而闻名[348].然而, 关于ZrTe5在磁场中的热电性质却很少有人研究.在存在垂直磁场和纵向热梯度的情况下, 载流子的扩散可以产生纵向电场 Ex=−Sxx·|∇T|(热电势)和横向电场Ey=Sxy·|∇T|(能斯特效应)[349].由于这些热电效应与电导率的导数成正比, 它们对异常的载流子输运尤为敏感, 并已被用于研究各种半金属和拓扑材料[350−359].

我们详细探讨了ZrTe5在磁场下的热电效应[86].图21 给出了热电势 − Sxx的温度依赖性.高温时,负的 − Sxx表明主要载流子是空穴.随着温度的降低, − Sxx在温度大约为 T∗=132 K 时变为正, 并且电阻率达到峰值, 表明ZrTe5从 p 型半导体转变为 n 型半金属, 这一点与以前的研究相一致[346,360].由于不同迁移率的电子和空穴载流子的存在, ZrTe5在室温下表现出巨大的磁热电势.其热电势在9 T 的磁场中可增加至370 µV/K, 比大多数材料的热电势都大, 但比外尔半金属的热电势要小[361].

图21 零磁场下ZrTe5 的电阻率 ρ (T) (黑色)和塞贝克系数 - Sxx(T) (蓝色)的温度依赖曲线.插图:实验测量示意图,B 是磁场, ∇ T 是温度梯度, a , b 和 c 是晶轴.转载自文献[86]Fig.21.Temperature dependence of the electrical resistivity ρ (T) (black) and Seebeck coefficient - Sxx(T) (blue)of ZrTe 5 at zero magnetic field.Inset: the measurement setup.B is the magnetic field and ∇ T is the temperature gradient.a ,b and c are crystallographic axes.Reproduced with permission from Ref.[86].

在更低温下, − Sxx甚至随着磁场的增加而准线性增长, 直至达到量子极限.这种线性热电势成为ZrTe5区别于其他材料的主要特征[353,354,356,357].在低温下实验样品中的平庸电子和非平庸电子均对热输运有贡献[330,360].为了区分常规和反常的能斯特信号, 使用经验公式来拟合弱场数据[357]

其中 µ 是载流子迁移率, B0是饱和磁场.和分别表示常规和反常能斯特信号的幅度.在弱场下, 这种经验表达式可以很好地拟合实验数据.根据我们的分析, 反常的能斯特信号出现在30 K以下, 并且随着温度的降低,急剧增加.在最低温度下, 反常信号和常规能斯特信号的幅度相当.所有结果表明, 在30 K 以下, 非平庸电子开始主要支配磁热电势的特性, 导致了 - Sxx(B) 和Sxy(B)的反常行为.

最近的研究表明, ZrTe5的整数朗道指数可以利用电阻率的峰值位置来确定[310,330,347,362].当磁场达到5.2 T 以上时, 所有电子占据最低朗道能带,即系统处于量子极限.此时随着磁场的增加, 狄拉克半金属或外尔半金属的热电势预期会非饱和地线性增长[363].

为了进一步阐明量子极限下不寻常的热电响应, 继续增加磁场至33 T (图22(a)), 此时 - Sxx开始增加.相应地, 在系统进入量子极限后, 能斯特信号出现一个峰值.另一个反常特征是能斯特信号Sxy在 B∗=14 T 附近符号发生了反转, 同时不同温度下的热电势 - Sxx收敛到零.之前, 在外尔半金属TaP 中观察到场诱导的 ρyx符号发生变化, 其中最低朗道能带在极强磁场下移动到化学势之上, 电子载流子数量急剧减少[64].载流子类型的改变可能是导致能斯特或霍尔信号的符号出现反转的原因.然而, 对于我们的样品, 低温下的电荷载流子只是电子[360].此外, 霍尔电阻 ρyx在 B∗附近变化平稳.因此, 载流子类型的变化不可能解释 Sxy的反常符号反转.

图22 (a)不同温度下对 - Sxx 的强场测量; (b)不同温度下 对 S xy 的强场测量.临界场 B ∗ 附近 S xy 改变符号, -Sxx收敛到零.转载自文献[86]Fig.22.(a) High-field measurements of - Sxx at several temperatures; (b) high-field measurements of S xy at several temperatures.Near the critical field B ∗ , S xy changes its sign and - Sxx converges to zero.Reproduced with permission from Ref.[86].

使用3D 有质量狄拉克模型探讨量子极限下反常 - Sxx和 Sxy的潜在机制[44].可以得到第0 个(ν =0)朗道能带为[86]

其中 s =±1 代表电子和空穴带, λ =±1 , v 是费米速度, m 是描述 γ 点附近小能隙的狄拉克质量[329].带 底时

因此, 对于有限质量 m , 第0 个电子和空穴的朗道能带之间的能隙随磁场先减小后增大.到达临界磁场 B∗=2m/(gµB) 时体系处于无能隙状态.相反, 在普通材料中, 第0 个电子和空穴的朗道能带之间的能隙随磁场单调增加.由于ZrTe5的狄拉克质量小, g 因子大, 在中强磁场下很容易观察到这种朗道能隙的闭合.此外, 对于这种低密度系统在强磁场下载流子密度是固定的,

其中 f (x) 是费米-狄拉克分布.由于第0 个朗道能带能隙闭合, 这个固定的 n 导致费米能级显著下降到量子极限之外.能斯特信号中的两项Sxy=ρxxαxy−ρyxαxx有相反的符号,当朗道能带无能隙时, 它们相互竞争并具有相同的绝对值, 导致 Sxy的符号反转.当磁场大于 B∗后最低朗道能带的能隙增大, 导致了热电势的增长.

其他狄拉克材料也观察到了类似的行为[226,354].我们对ZrTe5中反常热电性质的理解为探索拓扑物态中的新奇物理开辟了一条新途径.从这个意义上说, 磁热电研究为我们提供了一个强大、简单且有效的实验工具来探索拓扑材料中的奇异行为.

8 评论和展望

量子极限中的大多数理论都采用玻恩近似, 当磁长度比无序势的范围短得多时, 电子可能被同一杂质多次散射.2D 和3D 时玻恩近似可能不再是好的处理手段[317,364,365], 超越玻恩近似将是3D 系统在极强磁场下的一个具有挑战性的课题.

有关量子振荡的理论需要进一步完善[366−369],相移规则能否推广到Ⅱ型外尔半金属[367,370,371]以及其他拓扑物态体系[372,373]将是一个值得探讨的问题.

研究几何相位时, 参数空间通常是实空间或动量空间[95].但费米弧和外尔点形成的外尔轨道是一种新的物理现象, 几何相位在实空间和动量空间均有积累.特别是, 当电子在上表面和下表面之间隧穿时, 几何相位沿路径的积累具有厚度依赖性[291,374].最近, 一项新的实验利用这种厚度相关的相移来证明外尔轨道对观测到的量子化霍尔电阻的贡献[298], 未来有望会有更多的工作来验证这些理论.

在量子极限下, 我们的理论表明, 如果体系是3D 强拓扑绝缘体, 则可能出现最多两次电阻下降[81].令人惊讶的是, 最近的实验表明在量子极限下最多可以有5 次振荡[375,376].振荡作为磁场的函数遵循标度不变性定律, 与冷原子体系中的Efimov 束缚态非常相似, 类Efimov 束缚态可以用来理解偶然振荡[377,378].实验中电阻的直接拟合也证明了这一点[375,376,378].然而, 其他可能导致标度不变性振荡的机制将是人们广泛感兴趣的话题.

固体中由空间群描述的丰富对称性保护的能量节点, 可能导致在高能物理中没有对应粒子的无质量费米子[379−382].文献[383,384]提出了具有三重简并的费米子.作为介于狄拉克费米子和外尔费米子的中间费米子, 其对应的半金属体系具有新奇的输运特征[385−387].

更多的课题也将引起广泛的兴趣, 包括磁性外尔半金属的量子输运[388−397], 笼目(Kagome)铁磁体Fe3Sn2[398], 双外尔半金属[6,107], Ⅱ型外尔半金属[294,399−408], Hopflink 节线半金属[409−412], 超越量子极限的量子振荡[413], 拓扑半金属中的普适电导涨落[414], 狄拉克半金属中的AB 效应[415], 外尔半金属的表面准粒子干涉[416−419], Fano 效应[420].