摘要: 基于非平衡格林函数方法, 理论研究了
Z 轴方向局域交换场和电场对锗烯纳米带中电子自旋极化输运性质的影响. 结果表明对锗烯的边缘区域施加强度大于其2倍有效自旋轨道耦合强度的交换场, 可使自旋向上和向下电子的能带在不同的能量区间产生带隙, 从而实现对不同自旋取向电子的100%过滤. 提出了一种利用电场辅助降低自旋过滤效应所需阈值交换场强度的方法. 研究表明, 同时对中心器件区域施加局域交换场和电场, 可以在小于有效自旋轨道耦合的弱交换场强度下, 在较大的能量窗口区域过滤自旋向上或向下的电子. 增大局域交换场强度, 自旋过滤效应所对应的能量区间显著增大.
关键词: 锗烯 /
局域交换场 /
自旋过滤效应 English Abstract Spin filter effect of germanene nanoribbon controlled by local exchange field and electric field Xiang Yang 1 ,Zheng Jun 1 ,Li Chun-Lei 2 ,Guo Yong 3,4 1.College of Mathematics and Physics, Bohai University, Jinzhou 121013, China 2.College of Elementary Education, Capital Normal University, Beijing 100048, China 3.Department of Physics, Tsinghua University, Beijing 100084, China 4.Collaborative Innovation Center of Quantum Matter, Beijing 100084, China Fund Project: Project supported by the National Natural Science Foundation of China (Grant Nos. 11604021, 11574173), the Open Project of State Key Laboratory of Low-Dimensional Quantum Physics, China (Grant No. KF201806), and the Science-Technology Foundation from Education Commission of Beijing, China (Grant No. KM20180028022). Received Date: 27 May 2019Accepted Date: 17 June 2019Available Online: 01 September 2019Published Online: 20 September 2019 Abstract: Germanene, which has been synthesized recently, is a single-layered material composed of germanium atoms. Almost all the striking properties of graphene can be transferred to germanene, because both of them have the same honeycomb lattice structure. In contrast with graphene, germanene has a sizable band gap and spin dependent helical edge states, which make it attractive candidate for spintronic applications. By using the nonequilibrium Green’s function method, the effects of Z -axis local exchange field and electric field on spin-polarized transport properties in germanene nanoribbon are studied theoretically. The results reveal that by exerting an exchange field with a strength greater than twice the effective spin-orbit coupling to the edge region of germanene, the spin resolved band gaps can be generated in different energy ranges, and thus 100% filtering of spin-up or spin-down electrons can be achieved. We theoretically propose a method to reduce the threshold exchange field strength for the spin filter effect by using the electric field. The results show that the spin-filter effect can be achieved in a larger energy range under a weaker exchange field when the electric field are applied to the entire central device region. With the increase of the local exchange field intensity, the energy range corresponding to the spin filter effect will increase significantly. Keywords: germanene /local exchange field /spin-filter effect 全文HTML --> --> --> 1.引 言 随着对二维石墨烯研究的不断深入[1 ,2 ] , 研究人员把目光转向了与碳同族的硅、锗、锡元素的单层结构. 锗烯具有与石墨烯类似的六角蜂窝状晶格结构, 使得石墨烯中基于狄拉克电子的新奇物理性质几乎都可以在锗烯中实现. 另一方面, 由于锗烯具有非共面的翘曲结构[3 ,4 ] 和较强的本征自旋轨道耦合[5 ,6 ] , 锗烯拥有石墨烯所不具备的一些优异性质. 锗烯具有比石墨烯和硅烯更大的体带隙以及受时间反演对称性保护的自旋动量锁定的边缘态[7 ] . 边缘态完全依赖于材料的体拓扑结构, 其电子输运无耗散且不受杂质和缺陷的影响. 2014年以后多个研究组相继在实验中成功制备了硅烯、锗烯和锡烯[8 -16 ] . 近些年, IV族拓扑绝缘体纳米结构的电子自旋输运特性受到越来越多的理论和实验研究关注[6 ,17 -28 ] . 电子具有电荷和自旋两种属性, 传统的半导体器件主要利用了电子荷载电荷的特性. 自旋电子学强调自旋自由度在信息存储、传输等方面的应用[29 ] . 与传统电子元件相比, 以自旋自由度为信息载体的器件具有信息不易丢失、处理速度快、能耗小等优点, 有助于实现高速度、低能耗的量子器件. 如何有效地实现自旋的输运和操作已经成为自旋电子学领域中的重要研究课题[30 ,31 ] . 自旋过滤效应是指在一定的能量范围仅让拥有某种自旋取向的电子导通, 而阻止其他自旋取向的电子通过. 自旋过滤能得到高度自旋极化的载流子, 从而实现逻辑运算电路中的“0”或“1”. 本文主要研究局域交换场和电场对锗烯纳米带自旋相关电导和能带特性的影响. 如图1 所示, 将锗烯纳米带分成左电极、中间器件及右电极三个区域, 在中间器件区域的下半部分(即Y 轴方向$ [1– N_y/2]$ 原子区域)施加Z 轴方向的交换场, 对整个器件区域通过背电极施加电场. 将重点讨论自旋相关带隙的形成与自旋轨道耦合及交换场强度的关系, 以及如何通过电场辅助在较大能量窗口内实现自旋向上或自旋向下电子的完全过滤. 图 1 (a)单层锗烯自旋过滤器件结构示意图(左、右两侧蓝色区域为左、右电极区域; 中间红色实线包围区域内为中心器件区; 灰色和黄色分别表示Z 轴方向交换场和电场的施加区域); (b)单层锗烯自旋过滤器件侧视图(红色和蓝色小球分别代表锗烯翘曲层顶部和底部的锗原子) Figure1. (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 ] 其中$ H_{\rm L} $ ,$ H_{\rm R} $ ,$ H_{\rm C} $ 分别表示左端导线、右端导线和锗烯中间器件区域的哈密顿量; $ c_{i\sigma}^\dagger $ ($ c_{i\sigma} $ )表示在第i 个晶格格点产生(湮灭)一个自旋方向为σ 在位能为$ \varepsilon_i $ 的电子; 指数σ 表示自旋方向, $ \sigma = \pm 1 $ 分别对应自旋向上$ \uparrow $ 和自旋向下$ \downarrow $ 的电子; t 表示最近邻原子间的跃迁能; $ \left\langle i,j \right\rangle $ 表示考虑所有最近邻格点间的相互作用; $ \lambda_{\rm so} $ 是有效自旋轨道耦合强度; $\langle\langle i,j \rangle\rangle$ 表示只考虑次近邻的格点间的相互作用; 当次近邻电子之间的跃迁沿顺时针(逆时针)方向时$ \upsilon_{ij} = +1 $ ($ \upsilon_{ij} = -1 $ ). $ H_{\rm C} $ 式中第二项表示在中间器件区域的下半部分($ i\in[1\sim N_y\times0.5]$ )施加大小为M , 方向沿Z 轴的交换场, $ { \sigma}^z $ 为Z 轴方向的泡利矩阵. 第三项表示在整体器件区域$ ( i\in[1\sim N_y] )$ 施加强度为$ \lambda_E= \ell E_z$ 的Z 轴方向电场, $ \mu_i = \pm 1 $ 分别对应翘曲层顶部(图1(b) 中红色)和底部(蓝色)格点, 顶部和底部原子间距$\ell $ = 0.33 ?. 流经电极的电流可以定义为单位时间内电极中粒子数的变化, 利用稳态时的电流守恒, 流入左电极的电流可表示为Landauer-Büttiker公式的形式[34 ] 其中, $ f_{{\rm L}({\rm R})\sigma}(E) = \dfrac{1}{{\rm exp}[(E-E_{{\rm FL}({\rm R})\sigma})/k_{\rm B}T]+1} $ 是左(右)电极中电子的费米分布函数, $T_{{\rm LR}\sigma} = $ $ {\rm Tr}[\varGamma_{{\rm L}\sigma}G_\sigma^r\varGamma_{{\rm R}\sigma}G_\sigma^a] $ 是自旋为σ 的电子从左(L)电极输运到右(R)电极的隧穿函数. 中心器件区域的超前(延迟)格林函数 线宽函数为$ \varGamma_{{\rm L}({\rm R})\sigma} = {\rm i}[\varSigma_{{\rm L}({\rm R})\sigma}^r-\varSigma_{{\rm L}({\rm R})\sigma}^a] $ , 自能项$ \varSigma_{{\rm L}({\rm R})\sigma}^r = H_{{\rm L}({\rm R}){\rm C}}^\dagger g_{{\rm L}({\rm R})\sigma}^rH_{{\rm L}({\rm R}){\rm C}} $ . 表面格林函数$ g_{{\rm L}({\rm R})\sigma}^r $ 可以利用格林函数或传递矩阵方法求解[35 ,36 ] . 假设左右电极之间存在微小电势差和温度差, 将费米函数围绕费米能量$ E_{\rm F} $ 和温度T 做线性展开, 令$f_{0} = $ $ \{{\rm exp}[(E-E_{\rm F})/k_{\rm B}T]+1\}^{-1} $ 为温度梯度和偏压等于零时的费米分布函数, 可将自旋相关电导表示为3.结果与讨论 在数值计算中将锗烯最近邻跃迁能t = 1.3 eV作为能量单位, 则有效自旋轨道耦合强度$ \lambda_{\rm so} =$ 4.3 meV$ \approx0.033t $ [26 ] . 结构示意图1(a) 所示中心器件区域(黄色区域)X 轴和Y 轴方向的锗原子数分别为$ N_x = 6 $ 和$ N_y = 16 $ , 在计算中令$ N_x = N_y = $ 80, 并固定系统的平衡温度为液氦温度T = 4.2 K. 首先讨论仅对中心器件的下边缘区域施加Z 轴方向的交换场时, 交换场强度对电子自旋输运性质的影响. 图2 给出了不同交换场强度M = 0.01t , 0.03t , 0.05t , 0.07t (图2(a) )和0.09t (图2(b) )条件下, 自旋相关电导$ G_\sigma $ 随费米能量$ E_{\rm F} $ 的变化关系. 当交换场强度较弱时没有自旋极化现象出现, 如图2(a) 内插图所示, 自旋向上与自旋向下的电导值相等, 此时不同自旋取向的能带简并, 自旋极化率为0. 增大交换场强度, 使其大于2倍的自旋轨道耦合强度(M > 2$\lambda_{\rm so} $ )时, 如图2(a) 和图2(b) 所示, 自旋向上(下)的电导与自旋向下(上)的电导分别以$ E_{\rm F} = M/2 $ 和$ E_{\rm F} = -M/2 $ 为中心出现一段电导为0 (1)的能量窗口区域. 因此仅通过施加Z 轴方向的交换场, 可以实现同时在不同能量区间分别过滤自旋向上或自旋向下的电子, 并且随着交换场强度的增加, 过滤效应对应的能量窗口也随之变宽, 电子的自旋过滤效应也越明显. 图 2 施加不同强度的局域交换场M 时, 电导G 与电子费米能量$E_{\rm F}$ 之间的关系(红色实线表示自旋向上的电子, 蓝色虚线表示自旋向下的电子) (a) M = 0.01t , 0.03t , 0.05t , 0.07t 时的自旋相关电导; (b) M = 0.09t 时的自旋相关电导 Figure2. Relationship between the conductance G and electron Fermi energy $E_{\rm F}$ 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 spin-down electrons. 为了进一步理解自旋过滤效应产生的原因, 计算了交换场强度$ M = 0 $ 和$ M = 0.09t $ 时对应的能带和电导图, 如图3 所示. 图3(a) 给出的是无外场影响下锗烯电极的能带, 由于$ M = 0 $ 自旋向上(下)的上边缘态与自旋向下(上)的下边缘态简并, 且都相交于$ k_x(a) = \text{π} $ 处. 此外, 在边缘态中电子移动的方向与自旋相关: 上边缘态支持自旋向上(下)的电子从左至右(从右至左)移动, 下边缘态与之相反; 而体能带中电子输运方向不受限制. 施加Z 轴方向的局域交换场, 会破坏系统的时间反演对称性, 从而使得简并的边缘态发生劈裂. 当交换场的强度$ M>2\lambda_{\rm 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, 但是由于边缘态的自旋动量锁定特性, 自旋向上的电子是通过上边缘态正向输运, 自旋向下的电子是利用体能带传导. 图 3 (a)无外场时锗烯的电子能带结构; (b)交换场强度为$M=0.09t$ 时的电子能带结构; 其中实心(空心)三角表示上(下)边缘态, 圆点表示体能带, 红色(蓝色)对应自旋向上(向下)的电子; (c)交换场强度为$M=0.09t$ 时自旋相关电导随着电子能量的变化 Figure3. (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$ . 接下来讨论在电场辅助调控下, 自旋过滤效应产生的情况. 这里令电场和交换场大小相等, 且同为Z 轴方向. 图4(a) 给出了交换场强度$M = \lambda_E = $ $ 0.01t $ 和$ M = \lambda_E = 0.03t $ 时自旋相关电导$ G_{\sigma} $ 随费米能量$ E_{\rm F} $ 的变化. 如图4(a) 所示, 当施加很小的外场($ M = \lambda_E = 0.01t $ )时, 自旋向上与自旋向下的电导平台值均为$ G_\sigma = e^2/h $ . 但是当M = $ \lambda_E\approx \lambda_{\rm so} $ 时, 自旋向下的电导在能量$ E\in [-M/2,0] $ 附近区域出现$ G_\sigma = 0 $ 的零电导平台, 而此时自旋向上的电导值恒为1. 这是因为在交换场和电场共同作用下, 下边缘态能带发生如图4(c) 所示的劈裂. 对于自旋向上电子, 从图3(a) 和图4(c) 可以看出, 在各个能量区域都有供电子进行正向输运的边缘态或体能带, 考虑到左右电极的能带如图3(a) 所讨论过的情况, 自旋向上的电子可以从左电极隧穿到右电极. 但是由于电极与中心区域能带在部分能量区间失配, 电导平台上形成一系列的共振隧穿峰. 随着交换场与电场强度的增加(图4(b) ), 自旋相关的带隙宽度随之增大, 即使在$ M = \lambda_E<2\lambda_{\rm so} $ 条件下, 也能在较大能量窗口内实现自旋向上或自旋向下电子的100%过滤. 改变交换场的方向(从Z 轴方向调整为–Z 轴方向)可以在一定能量区间实现对自旋向下电子的过滤. 需要强调的是, 虽然在Z 轴方向交换场和电场的共同作用下, 只能过滤一种自旋方向的电子, 但是此时自旋过滤效应产生需要的交换场条件明显弱于仅有局域交换场存在的情况. 图 4 (a)局域交换场和电场强度为$M=\lambda_E=0.01t$ 和$M=\lambda_E=0.03t$ 时自旋相关电导随着电子费米能量$E_{\rm F}$ 的变化; (b)交换场和电场强度为$M=\lambda_E=0.05t$ 时自旋相关电导随着电子费米能量$E_{\rm F}$ 的变化; (c)交换场和电场强度为$M=\lambda_E=0.03t$ 时的电子能带结构; 其中实心(空心)三角表示上(下)边缘态, 圆点表示体能带, 红色(蓝色)对应自旋向上(向下)的电子 Figure4. (a) The spin-dependent conductance as a function of electron Fermi energy $E_{\rm F}$ when $M=\lambda_E=0.01t$ and $M=\lambda_E=0.03t$ ; (b) the spin-dependent conductance vs electron Fermi energy $E_{\rm F}$ when $M=\lambda_E=0.05t$ ; (c) energy-band diagram of germanene with local exchange field and electric field $M=\lambda_E=0.03t$ . 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. 4.结 论 理论研究了Z 轴方向局域交换场和电场对锗烯纳米带中自旋极化输运性质的影响, 提出了两种实现自旋过滤效应的方法. 研究发现, 仅对锗烯纳米带中心器件边缘区域施加Z 轴方向交换场, 当场强大于2倍本征自旋轨道耦合强度($ M>2\lambda_{\rm so} $ )时, 自旋向上和向下电子的能带分别在$ E_{\rm F} = M/2 $ 和$ E_{\rm F} = -M/2 $ 附近产生自旋相关的带隙, 因此可以在不同能量区间实现不同自旋方向电子的过滤效应. 此外研究表明, 电场可以有效降低自旋过滤效应所需的强交换场条件要求, 对整个中心区域施加电场后, 可以在较弱交换场条件下实现自旋向上或自旋向下电子的100%过滤. 增大局域交换场强度, 两种方法引起的自旋过滤效应都会显著增强.