全文HTML
--> --> --> -->2.1.极性半导体III-V族氮化物极化场起源
第三代半导体材料III-V族氮化物在白光照明和大功率电子器件中有重要而广泛的应用. 不同于一般的III-V族化合物半导体, 纤锌矿III-V族氮化物具有极强的自发极化与压电极化.其中自发极化源于平衡态纤锌矿晶格结构本身的成键对称性破缺, 以GaN为例, Ga原子与最近邻的4个N原子构成四面体结构, Ga原子作为正电中心, 其实际位置与N原子的有效负电中心并不重合, 正电中心与负电中心空间上的分离就造成了一个天然的电偶极矩, 从而形成沿[0001]晶向的宏观有效极化场(见图1(a)), 成为纤锌矿III-V族自发极化的微观起源. 而压电极化则源于晶格内部电偶极矩对晶格失配或外部机械应力的响应. 相比于自发极化, 压电极化的存在更为广泛, 仅要求晶格中心反演对称性破缺, 因此在III-V族氮化物材料中, 自发极化与压电极化常常是同时存在的. 值得注意的是, 计算表明[15], III-V族氮化物纤锌矿半导体的压电极化强度要比一般的III-V族或II-VI族半导体材料高出一个数量级. 界面晶格匹配关系可以设计压应力或张应力, 产生与自发极化或同向或反向的压电极化, 从而增强或减弱系统宏观极化场, 构成了纤锌矿极性半导体极性界面调控的物理基础(见图1(b)).
图 1 纤锌矿III-V族氮化物中自发极化(a)与压电极化(b)
Figure1. Schematics of spontaneous polarization (a) and piezoelectric polarization (b). The orientations of piezoelectric polarization can be tuned by external strain.
2
2.2.GaN/InN/GaN量子阱能隙调控
长期以来, 纤锌矿III族氮化物, 包括InN、GaN与AlN等材料, 由于在白光照明、高频大功率晶体管等方面的广泛应用引起人们广泛地关注. 这里我们首先考虑GaN/InN/GaN量子阱结构, 其中GaN与InN构成了I型异质结(参见图2). GaN为带隙约3.44 eV的直接带隙半导体, 其晶格常数为a = b = 3.182?; c = 5.175?; 而InN晶格常数为a = b = 3.548?; c = 5.751?, 为直接带隙0.64 eV的半导体. GaN/InN/GaN量子阱结构如图2左侧所示, 量子阱生长方向为与自发极化平行的[0001]方向. 由于InN与GaN平面内晶格失配约为12%, 极性界面应力引致的压电极化不容忽视, 与自发极化以及量子限域效应竞争, 共同决定了GaN/InN/GaN量子阱的电子结构.图 2 GaN/InN/GaN[0001]方向生长量子阱示意图(左), GaN、InN带阶关系(右)
Figure2. Schematics of GaN/InN/GaN QW (left panel) and the relative band alignment (right panel).
图3中(a)、(b)所示分别为InN层厚为2层和4层时GaN/InN/GaN量子阱的能隙, 二者直接带隙分别为0.8 eV和0 eV, 单层情形下能隙约为2 eV, 此时量子限域效应占优势. 随着InN层数的增加, 量子限域效应衰减, 极化场显著降低系统能隙直至关闭能隙, 甚至能带翻转. 图3(c)清楚地表明极性面和非极性界面对能隙的影响. 对比沿自发极性方向[0001]和非极性方向[1010]生长的GaN/InN/GaN量子阱带隙随InN层厚的改变, 非极性方向生长的量子阱, 其带隙总是高于极性方向生长的量子阱, 系统能隙无法被关闭. 而极性界面情形, 能隙在三层InN厚度时, 能隙接近零. 图3(d)展示了GaN/InN/GaN量子阱內建极化场随InN层厚改变的变化情况. 量子结构中內建极化场是自发极化与压电极化叠加的结果. 对于GaN/InN/GaN量子阱, GaN区域的压电极化与自发极化同向, 而InN区域的压电极化与自发极化反向, 因此对于沿[0001]方向生长的量子阱, 內建极化场强度随着InN层厚增加下降较[1010]方向量子阱更快, 但內建电场强度均为10 MV/cm量级, 远远超过了半导体栅极技术所能达到的强度.
图 3 InN厚度为2层(a)、4层(b)情况下GaN/InN/GaN量子阱Γ点直接带隙. (c) 极性方向[0001]与非极性方向[1010]生长GaN/InN/GaN量子阱带隙随InN层厚变化关系, 内插[1010]方向生长量子阱晶格示意图. (d) 內建极化场随量子阱中InN厚度变化关系, 蓝色与橙色点线分别表示极性方向[0001]和非极性方向[1010]生长量子阱情况. GaN、InN区域极化场方向由内插图例示意
Figure3. (a) and (b) Band structure of a GaN/InN/GaN QW around the Γ point for 2 and 4 ML of InN, based on first-principles DFT-HSE calculations. The green lines represent electron states, red lines light-hole states, and blue lines heavy-hole states. (c) calculated energy gaps as a function of the thickness of InN layers for polar ([0001]) (blue squares) and nonpolar ([1010]) (orange diamonds) GaN/InN/GaN QWs. The thicknesses were scaled to the thickness of 1 ML of InN in a polar QW and therefore correspond to the number of InN layers in the polar case. The inset shows a ([1010]) QW with two InN layers. (d) polarization field as a function of number of inserted InN layers calculated by DFT-HSE (blue squares) and based on theoretical polarization constants (orange circles). The blue and red block arrows in the inset show the polarization direction of GaN and InN regions.
2
2.3.GaN/InN/GaN量子阱中的拓扑相: k.p理论模型
基于群论不变量理论发展的多带k.p理论模型具有物理图像清晰、计算量相对较小、可考虑尺寸较大的量子体系、外场效应(如: 电场、磁场和应变场等)易于计入等特点, 已经广泛成功应用于半导体低维量子结构[16-18]的研究. 并且晶体的对称性要求使得利用k.p模型中少许待定参量就可以将量子体系在布里渊区高对称点的能谱性质完全刻画, 近年来我们采用30带的kp模型可以刻画整个布里渊区的电子能谱.我们发展了考虑应变效应的纤锌矿半导体结构八带 k.p哈密顿量. 在GaN/InN/GaN量子阱系统[19]中, InN的轻空穴和重空穴之间的劈裂很小, 这样相对于HgTe的四带BHZ模型[20], 我们构建一个有效六带k.p模型研究InN量子阱在极化电场下的电子结构.
图4(a) 表述的是量子阱能级随阱宽的变化关系, 从八带k.p模型计算结果看出, GaN/InN/GaN量子阱系统, 在阱宽大于15.5?的时候, 带隙反转. 这一临界转变厚度, 大约为5个InN单层厚度. k.p模型使得我们能够对压变, 极化电场和自旋轨道耦合作用相互之间的竞争和调制进行详细研究, 同时也方便于我们揭示拓扑相变背后蕴藏的物理图像. 图4(b)为GaN/InN/GaN量子阱的反转带隙能谱图, 从图可以看出具有Γ6对称性的电子态E1(红颜色标记)在具有Γ8对称性的空穴态(重空穴态HH1绿线标记, 轻空穴态LH1蓝线标记)下面, 能带发生反转. 内插图为E1, LH1和HH1在k = 0 点处的波函数分布. 自旋轨道耦合在拓扑相变中起着十分关键的作用. InN和GaN的本征自旋轨道耦合非常小, 数量级为大概几个meV. 我们发现强烈的界面极化电场诱导了一个可观的Rashba自旋轨道耦合. 量子阱中的Rashba自旋轨道耦合吸引了人们的广泛注意, 因为Rashba自旋轨道耦合可以通过门电压[21,22]和能带工程[23,24]调控. 在图4(d)中, 我们从八带Kane模型出发, 估算Rashba自旋轨道耦合的强度为1—2 meV数量级, 这与InAs和HgTe量子阱中的由于外加电场引起的Rashba自旋劈裂相当[25-27]. 通常这种较强的Rashba自旋轨道耦合只能发生在重原子系统中.
图 4 (a)GaN/InN/ GaN量子阱能级随阱宽的变化关系; (b)从八带Kane哈密顿量得到的16? GaN/InN/ GaN QW的能带结构, 红线, 蓝线和绿线分别代表电子态E1, 轻空穴态LH1和重空穴态HH1. 内插图为E1, LH1和HH1在波矢k = 0处的波函数分布; (c)通过求解有效六带模型得到的霍尔条带的能带结构. 左右插图显示了一对Kramers边缘态的密度分布: 左边是k|| = –0.1?–1处的自旋态, 右边是k|| = 0.1?–1处的自旋向下态. 中间插图为无限长霍尔条带的示意图, 其宽度(沿y)为1000埃. 沿[0001]生长方向的厚度(标记为z)包括InN量子阱区、加上两侧的两个200 ?厚的GaN势垒层. 黄线表示螺旋边缘态, 绿色箭头表示自旋方向. 短黑箭头表示极化电场方向. (d)电子(绿色), HH(蓝色)和LH(红色)子带的Rashba自旋劈裂(RSS)
Figure4. (a) The energy level varies with GaN/InN/GaN quantum wells width ; (b) the energy band structure of 16? GaN/InN/GaN QW obtained from the eight-band Kane Hamiltonian, the red line, the blue line and the green line represent the electronic state E1, the light hole state LH1 and the heavy hole state HH1, respectively. The inner illustrations are the wave function distribution of E1, LH1 and HH1 at the wave vector k = 0; (c) band structure of the Hall bar obtained by solving the effective six-band model. The left and right insets show the density distributions of one Kramers pair of edge states: on the left the spin-up state at k|| = –0.1?–1, on the right the spin-down state at k|| = –0.1?–1; The middle inset is schematic of an infinite long spin Hall bar with a width (along y) of 1000 ?. The thickness along the [0001] growth direction (labeled as z) comprises the InN QW plus two 200-? -thick GaN barrier layers on either side. The yellow lines show the helical edge states, and the green arrows show the spin orientation. The short black arrows indicate the polarization induced electric field; (d) rashba spin splitting (RSS) of electron (green), HH (blue) and LH (red) subbands
更为直接的证据拓扑相变是在拓扑绝缘体中, 在边界处总会存在奇数对Kramers简并对[5,28]. 边缘态能谱可以通过解有限宽度的自旋霍尔条带结构的有效六带模型而得到. 我们构造与生长方向垂直的自旋霍尔条带结构系统(见图4(c)内插图), 并计算了该系统的边缘态. 能带结构在图4(c)中呈现, 图4(c)左右内插图还给出了边缘态的Kramers对的密度分布. 结果明显的表明: 在电子态E1和轻空穴态LH1之间的带隙内出现了边缘态, 边缘态的分布在实空间是局域在自旋霍尔条带结构的边界. 该系统中拓扑非平庸的带隙大概为15 meV, 这比HgTe/CdTe系统中相当[29], 但InN量子阱中无需重原子元素存在, 并可以被集成于InN基的晶体管中, 原则上可以被广泛应用于高频和高功率器件中[30].
GaN/InN/GaN量子阱结构目前已经被实验成功制备[31], 其中InN的厚度可以少至1—2个原子层, 并且GaN、InN极性界面具有可原子级分辨的平整度, 我们所提出的GaN/InN/GaN量子阱中强烈的极化场及其能隙的大范围调控[19]为实验所验证[32]. 这一量子相变是由材料本身具有的强的极化电场驱动的, 这是首次在弱本征自旋轨道耦合和常规轻原子半导体系统中提出利用本征极化电场实现拓扑绝缘体相.
-->
3.1.IV族非极性半导体极化场起源
我们考虑GaAs-Ge极性界面系统, GaAs与Ge晶格常数均为5.65?, 二者晶格完美匹配, 同时, Ge本身具有中心反演对称性, 压电极化在该系统中的贡献微乎其微. 我们设计的GaAs/Ge/GaAs量子阱沿[111]方向生长, 这样的设计基于以下考虑: (1) 处于GaAs中间的Ge层, 其上下界面成键方式不同, 如图5(a)和5(c)所示, 上界面处Ge与As原子成极性键, 界面处积累电子, 下界面处Ge与Ga原子成极性键, 界面处积累空穴, 在Ge的量子结构产生内建极化场; (2) 极性界面使GaAs与Ge的布里渊区同时发生折叠, 其结果是倒空间L点的能带被折叠到Γ点, 如图5(b)所示, 由于体材料的Ge是导带底在L点的间接能隙半导体, 因此这样的界面设计可以使Ge的量子结构变为直接能隙.图 5 (a) GaAs/Ge/GaAs[111]量子阱示意图, (b)布里渊区折叠示意图, (c) GaAs/Ge/GaAs[111]量子阱上下极性界面电荷积累示意图.
Figure5. (a) Schematic of the structure of an ultrathin Ge layer sandwiched by thick GaAs layers (the upper left-hand panel). The upper right-hand panel amplifies the atomic configuration of the GaAs/Ge/GaAs quantum well containing four bilayer Ge. Notice that the Ga and As atoms locate at the opposite interfaces which leads to a charge accumulation schematically shown in the left-hand panel. (b) the BZ of bulk Ge and the folded BZ of GaAs = Ge = GaAs QW along the [111] crystallographic direction. (c) the charge accumulation at two opposite interfaces obtained from the first-principles calculation. The red and green isosurfaces describe the positive and negative charge accumulations at opposite interfaces.
另外需要指出的是, GaAs/Ge界面晶格常数的完美匹配使这一方案相较大晶格失配的GaN/InN/GaN量子阱更易生长.
2
3.2.GaAs/Ge/GaAs极性界面电子结构调控
尽管GaAs/Ge/GaAs量子阱中极化场的来源与GaN/InN/GaN量子阱截然不同, 极化场的强度仍可比拟, 约为10 MV/cm量级. GaAs/Ge/GaAs量子阱的带隙由量子限域效应和极化场竞争决定, 随着Ge层厚的增加, 带隙由1.6 eV单调下降至 0 eV(见图6(a)、(b)). 此外, GaAs/Ge/GaAs量子阱破坏了Ge材料本身的中心反演对称性, 体系也具有较强的压电效应. 由于GaAs与Ge晶格的完美匹配, 我们可以通过掺杂、机械应力等方式对量子阱施加平面内应力, 平面内的张应力会增强极化场, 使能隙降低, 而压应力则会削弱极化场, 是能隙增大. 对于1%范围内的形变, 带隙的变化是线性的, 且变化范围在0.13 eV左右(见图6(c)), 这使得我们可能通过Ge的层厚以及应力手段, 连续覆盖0–1.6 eV波段, 准确设计满足特定波段光电应用需求的GaAs/Ge/GaAs量子阱.图 6 不同Ge层厚GaAs/Ge/GaAs[111]量子阱的带隙, (a)从左至右分别对应Ge层厚为2、4个双原子层以及4个双原子层量子阱并考虑3%平面内张应力情况. 左下(b)为量子阱能隙以及內建电场大小随Ge层厚增加的变化关系. (d)平面内应力对能隙的改变.
Figure6. (a) Band structures of GaAs/Ge/GaAs sandwiched structures with different Ge portions obtained from the first-principles HSE calculations. From left to right, two Ge bilayers, four Ge bilayers, and four Ge bilayers with 3% in-plane tensile strain. (b) the band gap (purple line with diamonds) and the inner polarization field strength (blue line with squares) as functions of the number of Ge bilayers. (c) the variation of band gap Eg = Estrain g-Eg as a function of in-plane strain.
2
3.3.GaAs/Ge/GaAs量子阱的拓扑相: k.p理论模型
类比GaN/InN/GaN量子阱系统, 我们也尝试建立GaAs/Ge/GaAs量子阱k.p理论模型来刻画量子阱中的拓扑相变, 与InN不同的是Ge为间接带隙半导体材料. 为了完整刻画材料体系整个布里渊区的能带信息, 我们成功的发展了研究间接带隙半导体的三十带k.p理论模型[34]. 三十带k.p模型可以有效在整个布里渊区中较准确的描述半导体的能带结构.图7(a)为三十带k.p理论量子阱能带结构. 计算发现带隙附近最重要的态是自旋向上和向下的电子态(
图 7 (a)从30带k.p模型获得的具有18? 的Ge层厚度的GaAs/Ge/GaAs QW结构的能带结构. 插图显示了电子, 重空穴(HH)和轻空穴(LH)的Rashba自旋分裂(RSS). (b)通过求解有效四带模型得到的量子线的能带结构. 无能隙的边缘态由图中的红线表示. 中心插图显示了量子线的示意图和螺旋边缘态. 右(红色)和左峰(蓝色)分别描述了在k|| = -0.01?-1处的自旋向上和自旋向下边缘状态的密度分布.
Figure7. (a) Band structures of a GaAs/Ge/GaAs QW structure with 18 ? Ge layer thickness obtained from the 30-band k.p model. The inset shows the Rashba spin splitting (RSS) of electron, heavy hole (HH), and light hole (LH). (b) Band structure of the quantum wire obtained by solving the effective four-band model. The gapless edge states are shown by the red line. The central inset shows the schematic of the quantum wire and the helical edge states. The right (red) and left peaks (blue) describe the density distribution of the spin-up and spin-down edge states at k|| = -0.01?-1, respectively.
在GaAs/Ge/GaAs量子阱中, 由于具有Γ5对称性的轻空穴带远离对拓扑相变起主要作用的电子带和重空穴带, 而我们只关注体带隙(10 meV)内的边缘态, 因此我们将三十带哈密顿量约化到一个类似于BHZ的四带k.p模型, 以研究自旋霍尔条带中的边缘态性质. 二维拓扑绝缘体的基本特征是在二维拓扑绝缘体边界附近存在螺旋边缘态. 我们考虑宽度为1000?的量子线结构. 沿(111)生长方向上, 是由厚度为18?的Ge层夹在两个200?的GaAs层之间组成. 图7(b)为上述量子线的能带结构以及一对Kramers边缘态的密度分布. 如图所示, 新的能量分支出现并扫过整个体带隙, 这些态局域在量子线的边缘附近. 具有相同的面内动量k||, 自旋向上和向下的边缘态分布在量子线在相对边缘处. 这些螺旋边缘态的存在, 清楚地证明了这种二维GaAs/Ge/GaAs量子阱系统已被驱动之拓扑绝缘体相.