相阳 郑军† 李春雷 郭永

1) (渤海大学数理学院,锦州 121013)

2) (首都师范大学初等教育学院,北京 100048)

3) (清华大学物理系,北京 100084)

4) (量子物质科学协同创新中心,北京 100084)

1 引 言

随着对二维石墨烯研究的不断深入[1,2],研究人员把目光转向了与碳同族的硅、锗、锡元素的单层结构.锗烯具有与石墨烯类似的六角蜂窝状晶格结构,使得石墨烯中基于狄拉克电子的新奇物理性质几乎都可以在锗烯中实现.另一方面,由于锗烯具有非共面的翘曲结构[3,4]和较强的本征自旋轨道耦合[5,6],锗烯拥有石墨烯所不具备的一些优异性质.锗烯具有比石墨烯和硅烯更大的体带隙以及受时间反演对称性保护的自旋动量锁定的边缘态[7].边缘态完全依赖于材料的体拓扑结构,其电子输运无耗散且不受杂质和缺陷的影响.2014年以后多个研究组相继在实验中成功制备了硅烯、锗烯和锡烯[8−16].近些年,IV族拓扑绝缘体纳米结构的电子自旋输运特性受到越来越多的理论和实验研究关注[6,17−28].

电子具有电荷和自旋两种属性,传统的半导体器件主要利用了电子荷载电荷的特性.自旋电子学强调自旋自由度在信息存储、传输等方面的应用[29].与传统电子元件相比,以自旋自由度为信息载体的器件具有信息不易丢失、处理速度快、能耗小等优点,有助于实现高速度、低能耗的量子器件.如何有效地实现自旋的输运和操作已经成为自旋电子学领域中的重要研究课题[30,31].自旋过滤效应是指在一定的能量范围仅让拥有某种自旋取向的电子导通,而阻止其他自旋取向的电子通过.自旋过滤能得到高度自旋极化的载流子,从而实现逻辑运算电路中的“0”或“1”.

本文主要研究局域交换场和电场对锗烯纳米带自旋相关电导和能带特性的影响.如图1所示,将锗烯纳米带分成左电极、中间器件及右电极三个区域,在中间器件区域的下半部分(即Y轴方向[1–Ny/2]原子区域)施加Z轴方向的交换场,对整个器件区域通过背电极施加电场.将重点讨论自旋相关带隙的形成与自旋轨道耦合及交换场强度的关系,以及如何通过电场辅助在较大能量窗口内实现自旋向上或自旋向下电子的完全过滤.

图1 (a)单层锗烯自旋过滤器件结构示意图(左、右两侧蓝色区域为左、右电极区域; 中间红色实线包围区域内为中心器件区; 灰色和黄色分别表示Z轴方向交换场和电场的施加区域); (b)单层锗烯自旋过滤器件侧视图(红色和蓝色小球分别代表锗烯翘曲层顶部和底部的锗原子)Fig.1.(a) Schematic structure of of germanene spin filter.The blue regions on the left and right side are left and right lead.The area surrounded by red solid line is the central device.The regions shown in gray and yellow are exerted to the Z-direction exchange field and electric field.(b) Side view of germanene spin filter.Red and yellow spheres represent atoms in the top and bottom buckled-layer of germanene.

2 研究方法

在紧束缚近似下,双端锗烯结构的哈密顿量可以表示为[32,33]

其中HL,HR,HC分别表示左端导线、右端导线和锗烯中间器件区域的哈密顿量表示在第i个晶格格点产生(湮灭)一个自旋方向为s在位能为εi的电子; 指数s表示自旋方向,σ=±1 分别对应自旋向上↑和自旋向下↓的电子;t表示最近邻原子间的跃迁能;表示考虑所有最近邻格点间的相互作用;λso是有效自旋轨道耦合强度;表示只考虑次近邻的格点间的相互作用; 当次近邻电子之间的跃迁沿顺时针(逆时针)方向时υij=+1 (υij=−1).HC式中第二项表示在中间器件区域的下半部分(i∈[1∼Ny×0.5])施加大小为M,方向沿Z轴的交换场,σz为Z轴方向的泡利矩阵.第三项表示在整体器件区域(i∈[1∼Ny])施加强度为的Z轴方向电场,µi=±1 分别对应翘曲层顶部(图1(b)中红色)和底部(蓝色)格点,顶部和底部原子间距ℓ=0.33 Å.

流经电极的电流可以定义为单位时间内电极中粒子数的变化,利用稳态时的电流守恒,流入左电极的电流可表示为Landauer-Büttiker公式的形式[34]

3 结果与讨论

在数值计算中将锗烯最近邻跃迁能t=1.3 eV作为能量单位,则有效自旋轨道耦合强度λso=4.3 meV ≈0.033t[26].结构示意图1(a)所示中心器件区域(黄色区域)X轴和Y轴方向的锗原子数分别为Nx=6和Ny=16 ,在计算中令Nx=Ny=80,并固定系统的平衡温度为液氦温度T=4.2 K.

首先讨论仅对中心器件的下边缘区域施加Z轴方向的交换场时,交换场强度对电子自旋输运性质的影响.图2给出了不同交换场强度M=0.01t,0.03t,0.05t,0.07t(图2(a))和0.09t(图2(b))条件下,自旋相关电导Gσ随费米能量EF的变化关系.当交换场强度较弱时没有自旋极化现象出现,如图2(a)内插图所示,自旋向上与自旋向下的电导值相等,此时不同自旋取向的能带简并,自旋极化率为0.增大交换场强度,使其大于2倍的自旋轨道耦合强度(M> 2λso)时,如图2(a)和图2(b)所示,自旋向上(下)的电导与自旋向下(上)的电导分别以EF=M/2和EF=−M/2 为中心出现一段电导为0 (1)的能量窗口区域.因此仅通过施加Z轴方向的交换场,可以实现同时在不同能量区间分别过滤自旋向上或自旋向下的电子,并且随着交换场强度的增加,过滤效应对应的能量窗口也随之变宽,电子的自旋过滤效应也越明显.

图2 施加不同强度的局域交换场M时,电导G与电子费米能量 EF 之间的关系(红色实线表示自旋向上的电子,蓝色虚线表示自旋向下的电子) (a) M=0.01t,0.03t,0.05t,0.07t时的自旋相关电导; (b) M=0.09t时的自旋相关电导Fig.2.Relationship between the conductance G and electron Fermi energy EF when applying different values of local exchange field M: (a) Spin-dependent conductance with M=0.01t,0.03t,0.05t,0.07t; (b) spin-dependent conductance with M=0.09t.The red solid line represents spin-up electrons and the blue dash line represents spindown electrons.

为了进一步理解自旋过滤效应产生的原因,计算了交换场强度M=0和M=0.09t时对应的能带和电导图,如图3所示.图3(a)给出的是无外场影响下锗烯电极的能带,由于M=0 自旋向上(下)的上边缘态与自旋向下(上)的下边缘态简并,且都相交于kx(a)=π 处.此外,在边缘态中电子移动的方向与自旋相关: 上边缘态支持自旋向上(下)的电子从左至右(从右至左)移动,下边缘态与之相反; 而体能带中电子输运方向不受限制.施加Z轴方向的局域交换场,会破坏系统的时间反演对称性,从而使得简并的边缘态发生劈裂.当交换场的强度M>2λso时,分别在能量E=M/2 和E=−M/2的附近产生自旋相关的带隙,并且带隙的宽度随交换场强度的增强而增大.如图3(b)所示,当施加强度为M=0.09t的交换场时,自旋向上与自旋向下的能带分别在E=0.045t与E=−0.045t为中心处产生带隙,此时左电极中能量为E=M/2(E=−M/2)的自旋向上(下)电子无法通过器件中心区域的带隙,但是自旋向下(上)的电子可以利用体能带进行输运,从而导致不同自旋取向的电子在不同能量窗口发生自旋过滤现象.值得注意的是,虽然费米面附近区域不同自旋方向的电导值均为1,但是由于边缘态的自旋动量锁定特性,自旋向上的电子是通过上边缘态正向输运,自旋向下的电子是利用体能带传导.

接下来讨论在电场辅助调控下,自旋过滤效应产生的情况.这里令电场和交换场大小相等,且同为Z轴方向.图4(a)给出了交换场强度M=λE=0.01t和M=λE=0.03t时自旋相关电导Gσ随费米能量EF的变化.如图4(a)所示,当施加很小的外场(M=λE=0.01t)时,自旋向上与自旋向下的电导平台值均为Gσ=e2/h.但是当M=λE≈λso时,自旋向下的电导在能量E∈[−M/2,0]附近区域出现Gσ=0 的零电导平台,而此时自旋向上的电导值恒为1.这是因为在交换场和电场共同作用下,下边缘态能带发生如图4(c)所示的劈裂.对于自旋向上电子,从图3(a)和图4(c)可以看出,在各个能量区域都有供电子进行正向输运的边缘态或体能带,考虑到左右电极的能带如图3(a)所讨论过的情况,自旋向上的电子可以从左电极隧穿到右电极.但是由于电极与中心区域能带在部分能量区间失配,电导平台上形成一系列的共振隧穿峰.随着交换场与电场强度的增加(图4(b)),自旋相关的带隙宽度随之增大,即使在M=λE<2λso条件下,也能在较大能量窗口内实现自旋向上或自旋向下电子的100%过滤.改变交换场的方向(从Z轴方向调整为–Z轴方向)可以在一定能量区间实现对自旋向下电子的过滤.需要强调的是,虽然在Z轴方向交换场和电场的共同作用下,只能过滤一种自旋方向的电子,但是此时自旋过滤效应产生需要的交换场条件明显弱于仅有局域交换场存在的情况.

图3 (a)无外场时锗烯的电子能带结构; (b)交换场强度为 M=0.09t 时的电子能带结构; 其中实心(空心)三角表示上(下)边缘态,圆点表示体能带,红色(蓝色)对应自旋向上(向下)的电子; (c)交换场强度为 M=0.09t 时自旋相关电导随着电子能量的变化Fig.3.(a) Energy-band diagram of germanene without external field; (b) energy-band diagram of germanene with local exchange field M=0.09t ; the solid (hollow) triangle represents the upper (lower) edge state,the dot represents the energy band,the red(blue) colour corresponds to the spin-up (spin-down) electron; (c) the spin-dependent conductance as a function of electronic energy when M=0.09t.

4 结 论

理论研究了Z轴方向局域交换场和电场对锗烯纳米带中自旋极化输运性质的影响,提出了两种实现自旋过滤效应的方法.研究发现,仅对锗烯纳米带中心器件边缘区域施加Z轴方向交换场,当场强大于2倍本征自旋轨道耦合强度(M>2λso)时,自旋向上和向下电子的能带分别在EF=M/2和EF=−M/2 附近产生自旋相关的带隙,因此可以在不同能量区间实现不同自旋方向电子的过滤效应.此外研究表明,电场可以有效降低自旋过滤效应所需的强交换场条件要求,对整个中心区域施加电场后,可以在较弱交换场条件下实现自旋向上或自旋向下电子的100%过滤.增大局域交换场强度,两种方法引起的自旋过滤效应都会显著增强.