自旋轨道耦合系统中的整数量子霍尔效应

梁滔 李铭

(华南师范大学物理与电信工程学院物理系,广东省量子调控工程与材料重点实验室,广州 510006)

1 引 言

整数量子霍尔效应是由德国物理学家Klitzing等[1]于1980年在研究二维电子气时发现的.在低温和强磁场下二维电子气的霍尔电阻呈现跳跃式变化,在某些磁场区域出现(平台,) 而且在平台位置纵向电阻为零.霍尔电阻以h/e2为单位量化,即后来大量的实验表明,霍尔平台的电阻是高度精确的,与具体材料无关,以至有关国际组织把h/e2作为新的电阻标准.继 Klitzing 的这个发现之后,Tsui等[2]又在更强磁场下(30 T左右)发现了分数量子霍尔效应.霍尔电阻在ν为1/3等分数值也出现霍尔平台[2−4].2013年清华大学Xue等[5]首次观测到反常量子霍尔效应,他们在极低温(30 mK)和无外磁场的条件下,通过改变所施加的栅极电压调节费米面所在的位置,在磁性掺杂的拓扑绝缘体Cr0.15(Bi0.1Sb0.9)1.85Te3薄膜中观测到了整数为1的霍尔电阻平台.现在人们还在努力寻找具有整数为2以上的霍尔电阻平台的量子反常霍尔效应.

目前人们对整数量子霍尔效应提出了多种理论解释,如拓扑不变性理论、安德森局域化理论、量子输运理论等.1982年索利斯和合作者发现整数量子霍尔平台实际上是一种拓扑量子态,对应的拓扑不变量是第一类拓扑陈数(称为TKNN数)[6,7],这种霍尔电导和陈数的关系为TKNN数)只适用于无穷大二维晶体,而且难以解释纵向电阻为0的事实.

整数量子霍尔效应的安德森局域化解释是被广泛接受的一种理论解释,写进了很多教材[8].这个解释假设材料含有缺陷(杂质),缺陷使朗道能级展宽成有限宽度的能带[8−12].能带中央是扩展态,边沿是局域态.随着磁场强度增大,朗道能级(能带)一个接一个扫过费米能级.当费米能级位于局域态中,局域电子对电导无贡献,霍尔电阻不变从而呈现霍尔平台.当费米能级位于扩展态中,扩展态电子对电流有贡献,产生电阻,从而霍尔电阻从一个平台上升到另一个平台[8−12].这种局域化解释依赖杂质的存在.下面将会看到,整数量子霍尔效应实际上并不需要杂质的存在,也不需要安德森局域化的能带模型.

人们用介观系统的电子输运原理对整数量子霍尔效应现象做了解释,见Datta的著作[13].在强磁场作用下的一个有限宽度薄导体中,具有相反大波矢的两电子状态在空间上是分离的,即从左到右传导电流的电子局域在导体一个边缘上,而从右到左传导的电子局域在另一条边缘上.这样的电流不能发生背散射,因而不能产生纵向电阻.电子输运给出的霍尔平台电导为是在费米(能级以)下的朗道能级数目.这些霍尔平台以2e2/h为单位量化,其中的2 是因为自旋在每个朗道能级上简并.Datta在其著作中指出,自旋态在磁场作用下发生塞曼分裂,简并解除,从而获得通常的量子化平台但没有做更细节的讨论[13].

但是,简单的计算表明,自旋态的塞曼分裂正好跟朗道能级间距相同,分裂后的朗道能级仍然是简并的,霍尔电导只出现奇数平台.因此,本文在介观系统中量子输运的理论基础之上,尝试加入Rashba自旋轨道耦合[14−16],在强磁场和塞曼效应的共同作用下,解释二维无限长条形样品的整数量子霍尔效应.

2 电子输运

考虑一个二维无限长条形样品.在与样品平面垂直方向施加均匀磁场B=.设样品的x方向是无限长的,y方向具有宽度W.矢势选择为朗道规范A=(−By,0).系统哈密顿量写为

其中m∗是样品中电子的有效质量,第二项是电子的Rashba自旋轨道耦合项,第三项是外磁场下自旋的塞曼项,第四项是霍尔势,其中Vh是霍尔电压,最后一项是样品边界产生的边界势Vc(y).边界势选为硬墙形式,即

由于系统的哈密顿量中不显含x坐标,电子的定态波函数具有以下形式:

把这个波函数代入薛定谔方程可得η(y),ξ(y) 满足的本征方程.为了方便计算,对薛定谔方程做无量纲化处理,即用无量纲的长度y代替y/a,用无量纲的波矢k代替ka,长度单位a约为 6 mm.于是得到无量纲的本征方程如下:

首先不考虑自旋轨道耦合 (λh=0).得到的能谱如图1(a)所示.容易看出,由于样品两边缘的限制,朗道能级在大波矢区间快速上升,在中间部分由于霍尔势的作用略微倾斜,形成了能带.后面将看到,正是这些上升部分的能级提供了传导电流.朗道能级在磁场作用下发生自旋分裂.最低的朗道能级只含有下自旋态,但从第二个朗道能级开始,能级仍然是上下自旋简并的.这是因为朗道能级的间距恰好等于自旋的塞曼分裂,即其中ωc=eB/m∗是电子回旋频率.

图1 (a)自旋简并的 3 个朗道能级在磁场中分裂,但仍有能级简并; (b)第一朗道能级波函数的模,参数 α=0.06,λh=0,W=50,µ1=0.49,µ2=0.46Fig.1.(a) The lowest three spin degenerate Landau levels are split apart in the strong magnetic field,but the degeneracy is not fully resolved; (b) the modulus of the wave functions of the first Landau level for α=0.06,λh=0,W=50,µ1=0.49,µ2=0.46.

图1(b)是最低朗道能级波函数的模的分布图.每个波矢k对应的波函数在y方向几乎都是局域的.尤其大波矢区间波函数都局域在样品的两个边缘上.所以,传导电流主要发生在样品的边缘上.具有相反大波矢的电子在空间上是完全分离的.它们产生相反方向的电流.正是这个原因,大波矢电子不能在晶格边沿上发生背散射从而不能产生电阻.而朗道能级的中间略微倾斜的部分对应的波函数局域在样品内部.这部分的波函数重叠非常小.其他朗道能级的波函数也是同样的情况.

处于能带右边的斜率为正的上升大波矢部分的电子局域在样品上边缘且向x轴正向(右)移动,而处于左边斜率为负的大波矢部分的电子局域在下边缘且向左移动,如图2(a)所示.在电子输运过程中,电子可看做准经典波包,波包的群速度由能带的斜率确定:

所以,处于能带中间略微倾斜部分的电子群速度很小,流动性很小,传导电流主要发生在样品的两条边沿上.x方向的总电流由下式给出:

图2 (a)在均匀磁场中通电的长条形,箭头表示电子的输运方向; (b)自旋简并下霍尔电导的奇数霍尔平台.参数kBT=0.002,µ1=0.49,µ2=0.46Fig.2.(a) An long slab with current flows in a strong magnetic field,arrows indicate the direction of electron transport; (b) odd plateaus of Hall conductance due to the spin degeneracy.kBT=0.002,µ1=0.49,µ2=0.46.

其中En0是能带n的最低能量值.

上下两边缘的电子不能背散射,只能单向流动,因而只能分别来自左右两端电极.它们不产生电阻,因而分别具有左右电极中电子的化学势.于是,样品两电极的化学势差正好给出霍尔电压Vh,即−eVh=µ1− µ2.无量纲化的化学势给出无量纲参数Fh=−(µ1− µ2)/W.最后(7)式的计算给出霍尔电导σ=I/Vh,如图2(b)所示.

现在在上面的计算中加入电子的Rashba自旋轨道耦合.自旋轨道耦合对系统能量有显著影响[17].这个耦合产生的附加能量彻底解除了能级简并,如图3(a)所示.这正是理解整数量子霍尔效应的关键所在.这时,这些能级的波函数都是上下自旋混合的,边界态都包含上下自旋.取自旋轨道耦合系数λh=0.7,重新数值计算(7)式,得到霍尔电导如图3(b)所示.其中霍尔平台准确给出仍然是费米能级所切割的能带数目.这样就得到了完整的整数量子霍尔效应.

从(7)式可以看出,霍尔平台上电流的大小取决于费米能级切割的能带数目,而与能带的斜率变化快慢无关.所以,无论取硬墙形式的边界势,还是取其他光滑的边界势,霍尔平台都会出现在同样的磁场范围.虽然Rashba自旋轨道耦合使朗道能级进一步劈裂的过程会导致朗道能级的间隙不相等,但这对整数量子霍尔效应的霍尔平台的位置没有影响,只会影响霍尔平台的宽度.只要费米能级位于两个朗道能级之间,电导就会出现霍尔平台.而且随着磁场的增大,朗道能级的间距整体上是变大的,霍尔平台的长度随着磁场增强而变长.可见,整数量子霍尔效应不过是个无散射的边界效应.这个解释不需要引入杂质,也不需要能级展宽.

(7)式也很清楚地表明,温度对霍尔平台有影响,计算结果如图4(a)所示.当样品温度超过一个临界值Tc(在本文参数下kBTc约为 0.003),霍尔平台的霍尔电导开始显著偏离整数.根据前面的能量单位,对于一般材料(设电子有效质量m∗=0.15m)这个临界温度大约为Tc=0.5 K.量子霍尔效应作为国际电阻标准,样品超过这个临界温度会带来显著的误差.在此临界温度以下,霍尔平台的值也因温度效应有一定误差.在温度远低于临界温度时,霍尔电导近似为

图3 (a)自旋轨道耦合作用下最低的 6个朗道能级; (b)自旋轨道耦合作用下的霍尔电导(蓝线)和纵向电阻(红线).参数α=0.06,λh=0.7,µ1=−0.02,µ2=−0.04,kBT=0.002Fig.3.(a) The lowest six Landau levels under spin-orbit coupling; (b) corresponding Hall plateaus of the conductance (blue) and the longitudinal resistance (red).α=0.06,λh=0.7,µ1=−0.02,µ2=−0.04,kBT=0.002.

图4 (a) 温度对霍尔平台的影响; (b) 第二朗道能级上几个波矢 k 对应的波函数的模 (上自旋分量).参数 λh=0.7,µ1=−0.02,µ2= −0.04Fig.4.(a) The effect of temperature on Hall plateaus; (b) modulus of the wave functions of a few k vectors on the second Landau level.λh=0.7,µ1=−0.02,µ2=−0.04.

3 纵向电阻和背散射效应

量子霍尔效应的纵向电阻表现出明显的振荡行为[18],即在霍尔电阻平台对应的区域纵向电阻为零,而在霍尔平台之间会出现峰值.从朗道能级的形成来看,当费米能级从一个朗道能级移动到下一个朗道能级的中心时,纵向电阻出现一次0电阻低谷.那么,在以上的Rashba自旋轨道耦合系统中,当费米能级非常靠近朗道能级的底部时,具有相反波矢的波函数在空间上显著重叠,从而发生背散射,产生了纵向电阻.

在半导体材料中,一般主要散射是杂质散射和晶格散射[19].然而,在如上所述的Rashba自旋轨道耦合作用下的二维无限长条形样品中,整数量子霍尔效应的计算和解释并不需要引入杂质.为了方便计算纵向电阻,假设该样品是没有杂质的,我们只考虑晶格散射对纵向电阻的影响.晶格振动对电子散射的跃迁几率为

其中 ℏω是声子能量,是散射势,ϑnkn′k′(b=1),Γnkn′k′(b=−1),分别为发射和吸收声子时的跃迁几率.

根据玻尔兹曼方程,电导可由费米面上电子碰撞的弛豫时间τF决定[20,21],即σ=ne2τF/m,而弛豫时间由下式确定[22,23]:

其中θkk′是散射角.由于电子输运是一维的,θkk′=0,π.只有当θkk′= π,即k与k′异号时,(10)式的积分才能对背散射有贡献,从而产生纵向电阻.为了简化,计算中散射势取为一个常数,声子能量最大取到 2 0kBT.对(10)式进行数值计算得到纵向电阻随磁场的变化如图3(b)所示.纵向电阻在霍尔平台降为 0,而在霍尔平台跳跃的时候出现,与实验结果一致.

第二朗道能级上几个波矢对应的波函数(上自旋部分)的模显示在图4(b)中.容易看出,具有相反大波矢的两波函数都局域在样品两边沿,没有重叠,因而近 似为 0 ; 具有相反小波矢的两波函数局域在在样品中心区域,有显著重叠因而对纵向电阻有贡献.只有当费米能级靠近朗道能级的底部时(对应图3(b)中霍尔平台跳跃区),才能发生声子的发射和吸收,从而产生纵向电阻.

4 结 论

研究了Rashba自旋轨道耦合作用下二维无限长条形样品中的电子输运,计算了样品的霍尔电导和纵向电阻,得到了完整的整数量子霍尔效应.在一定强磁场范围内,由于样品两边缘的限制,能级在大波矢范围快速上升,在小波矢范围形成平坦的朗道能级.强磁场下自旋轨道耦合完全解除了自旋简并.这是形成整数量子霍尔效应的关键.位于朗道能级上升和下降区域的电子局域在样品边沿,形成了分离的双向电流.计算结果表明,霍尔电导呈现台阶型,平台出现在e2/h的整数倍位置,形成霍尔平台.温度对霍尔平台的电导有一定影响.在某临界温度以下,霍尔平台电导可以达到10−9以上的精度.本文最后分析了声子发射和吸收产生整数量子霍尔效应的纵向电阻的机制,近似计算了弛豫时间,得到了纵向电阻.结果表明,纵向电阻在霍尔平台区域为零,而在霍尔平台之间出现峰值.