局域势调控下石墨烯/h-BN异质结的能带结构和量子相变

Band Structure and Quantum Phase Transition of Graphene/h-BN Heterojunction under Local Potential Control

摘要: 采用紧束缚方法研究了局域势调控下石墨烯/六方氮化硼异质结的能带结构及量子相变。结果表明，在给定石墨烯层内禀自旋轨道耦合强度时，石墨烯层处于量子自旋霍尔态而六方氮化硼层处于绝缘态，增大石墨烯层局域势，系统会从量子自旋霍尔态变到半导体态。调节六方氮化硼层局域势可以使异质结产生新的无间隙边缘态，从而使系统中石墨烯层的量子自旋霍尔态变成由层内和层间的边缘态共同构成的量子自旋霍尔态。

Abstract: The band structure and quantum phase transition of graphene/hexagonal boron nitride heterojunction (h-BN) under local potentials are studied by the tight-binding method. The result shows that the graphene layer is in the quantum spin Hall state and the h-BN layer is in the insulating state when the intrinsic spin-orbit coupling strength of the graphene layer is given. As the local potential of the graphene layer increases, the system will change from a quantum spin Hall state to a semiconductor state. New gapless edge states can be generated by tuning the local potential of the h-BN layer, so that the quantum spin Hall state of the graphene layer becomes the quantum spin Hall state composed of the edge states within and between layers.

1. 介绍

石墨烯是一种令人着迷的材料，有着较高的载流子迁移率，有希望成为新一代电子元件 [1] [2] 。但是石墨烯是零带隙，在制作元件时会导致开/关电流比很小，所以打开石墨烯带隙变得至关重要。通过在石墨烯底部施加衬底，破坏石墨烯空间反演对称性从而可以打开带隙 [3] [4] [5] 。研究人员已经通过理论计算和实验证明了六方氮化硼(h-BN)作为石墨烯衬底非常合适。h-BN是有着较大带隙的绝缘体，与石墨烯形成范德瓦尔斯异质结不仅可以破坏石墨烯空间反演对称性打开能隙，还可以保持石墨烯高迁移率等特点 [6] [7] [8] 。对于石墨烯/h-BN异质结在外场调控下的能带结构以及电子和输运已经有了广泛的研究 [9] [10] ，但是与其自旋相关的性质还了解的比较少，有待进一步深入研究。

拓扑绝缘体是一种新的拓扑状态，它的内部是绝缘状态，而边缘是导电状态，电荷只能通过边缘传播。2005年，Kane和Mele在单层石墨烯中施加内禀自旋轨道耦合(ISOC)，首次预言了一种非平庸拓扑状态，量子自旋霍尔(QSH)状态 [11] [12] 。QSH状态是一种二维拓扑绝缘态，以自旋和速度方向锁定的螺旋形边缘态为特征。QSH效应受到时间反演对称性保护，其Z₂拓扑不变量ν = 1。当自旋守恒时，还可以引入自旋Chern数(Cs)来表征此状态，此时Cs与Z₂等价。但与Z₂不变量不同的是，即使时间反演对称性被破坏，Cs依然有效 [13] [14] [15] 。在之前的工作中，我们已经知道二维电子系统考虑到ISOC时可以实现QSH效应，产生非零的量子化自旋霍尔电导，而在此基础上考虑强Rashba SOC或强无序会破坏QSH效应的稳定性，从而产生非量子化的自旋霍尔电导，形成自旋霍尔效应 [14] [16] [17] 。另外，外场的调控会使处于QSH效应或weak-QSH效应的单双层石墨烯发生拓扑相变，形成新的拓扑状态 [18] [19] [20] 。这一系列的研究都推动了自旋电子学的进步，有着重要的理论和应用价值。

基于石墨烯的这些特殊性质，研究者对石墨烯/h-BN异质结进行了大量的研究，发现了很多有趣的性质，例如次级狄拉克锥的出现，拓扑电流 [21] [22] [23] 等。但是石墨烯/h-BN纳米带(GBNNR)的许多性质仍不清楚，特别是受局域势调控下的石墨烯/h-BN异质结的能带结构和拓扑量子相变性质。在本文中，我们在晶格匹配的条件下采用紧束缚方法研究Bernal堆叠zigzag边界GBNNR在局域势调控下的能带结构及量子相变。

2. 模型与方法

单层石墨烯和单层h-BN在局域势作用下的紧束缚哈密顿量分别为：

$H_{MLG}=-t_1{\displaystyle \underset{\langle i,j\rangle ,\alpha }{\sum }{C}_{i\alpha }^{†}}{C}_{j\alpha }+i\lambda {\displaystyle \underset{\langle \langle i,j\rangle \rangle ,\alpha \beta }{\sum }{\upsilon }_{ij}{C}_{i\alpha }^{†}{\left(S_z\right)}_{\alpha \beta }{C}_{j\beta }+{\displaystyle \underset{i,\alpha }{\sum }{V}_i}}{C}_{i\alpha }^{†}{C}_{i\alpha }$ (1)

$H_{h\text{-}BN}=-t_2{\displaystyle \underset{\langle i,j\rangle \alpha }{\sum }{C}_{i\alpha }^{†}{C}_{j\alpha }+{\displaystyle \underset{i,\alpha }{\sum }{E}_i{C}_{i\alpha }^{†}{C}_{i\alpha }}}+{\displaystyle \underset{i,\alpha }{\sum }{V}_i{C}_{i\alpha }^{†}{C}_{i\alpha }}$ (2)

其中 $H_{MLB}$ 和 $H_{h\text{-}BN}$ 分别表示单层石墨烯和单层h-BN的紧束缚哈密顿量。(1)和(2)中第一项描述最近邻格点之间的跃迁，石墨烯层和h-BN层最近邻跃迁能分别为t₁和t₂， $C_{i\alpha }^{†}\left(C_{j\alpha }\right)$ 表示自旋为 $\alpha$ 电子的产生(湮灭)算符；第三项中V_i表示石墨烯层和h-BN层原子在i格点的局域势。(1)中第二项为ISOC项， $\lambda$ 为石墨烯层ISOC强度， $\langle \langle i,j\rangle \rangle$ 表示次近邻原子跃迁 ${\upsilon }_{ij}=\pm 1$ ，表示次近邻原子跃迁沿着逆时针(顺时针)为+1 (−1)， $\alpha \left(\beta \right)=\pm$ 是自旋指标。(2)中第二项E_i表示h-BN层B原子和N原子在i格点的在位能，分别用E_B和E_N表示。

Bernal堆叠石墨烯/h-BN异质结的哈密顿量为：

$H=H_{MLB}+H_{h\text{-}BN}+t_{\perp }{\displaystyle \underset{i\in T,j\in B,\alpha }{\sum }\left(C_{i\alpha }^{†}{C}_{j\alpha }+C_{j\alpha }^{†}{C}_{i\alpha }\right)}$ (3)

第一、二项分别为(1)和(2)中单层石墨烯和单层h-BN的哈密顿量，最后一项为层间相互作用项，B和T分别表示为h-BN层(底层)和石墨烯层(顶层)。

在研究过程中，我们通过下面公式计算Chern数 [24] [25] [26] ：

$C_{\alpha \beta }=\frac{1}{2\pi }{\displaystyle \underset{n}{\sum }{\displaystyle {\int }_{BZ}\text{d}{k}_x}}\text{d}{k}_y{\left({\Omega }_{xy}\right)}_{\alpha \beta }$ (4)

其中， $C_{\alpha \beta }$ 是自旋为 $\alpha \beta$ 部分的Chern数， $\alpha \left(\beta \right)=\pm$ 是自旋指标， ${\Omega }_{xy}$ 是第n个能带在动量空间的Berry曲率 [19] [27] ，表达式如下：

${\Omega }_{xy}=-{\displaystyle \underset{n^{\prime }\ne n}{\sum }\frac{2\mathrm{Im}\langle {\Psi }_{nk}|v_x|{\Psi }_{n^{\prime }k}\rangle \langle {\Psi }_{n^{\prime }k}|v_y|{\Psi }_{nk}\rangle }{{\omega }_n-{\omega }_{n^{\prime }}}}$ (5)

求和是对费米能以下所有被占据能带进行求和， $v_{x,y}=\frac{\partial H}{\partial k_{x,y}}$ 是速度算符。当自旋守恒时，Cs可以写成被占据空间自旋向上(C₊)和自旋向下(C₋)之差的一半 $Cs=\left(C_{+}-C_{-}\right)/2$ 。

Figure 1. (a) Schematic diagram of Bernal stacking zigzag-edged GBNNR. The bottom layer is h-BN layer and the top layer is graphene layer. y is the periodic direction and the x direction contains N = 80 sites. Green, yellow and blue solid circles represent N, B and C atoms, respectively; (b) Schematic diagram of graphene/h-BN heterojunction under local potentials.V₁ and V₂ are applied in the upper and lower half regions of the graphene layer respectively, and V₃ and V₄ are applied in the upper and lower half regions of the h-BN layer respectively

图1. (a) Bernal堆叠zigzag边界GBNNR示意图。底层为h-BN层，顶层为石墨烯层。y为周期性方向，x方向包含N = 80个格点。绿色、黄色和蓝色实心圆分别表示N原子，B原子和C原子；(b) 局域势作用于石墨烯/h-BN异质结示意图。V₁和V₂分别被施加在石墨烯层的上半部分区域和下半部分区域，V₃和V₄分别被施加在h-BN层的上半部分区域和下半部分区域

3. 结果与讨论

论文研究的石墨烯/h-BN异质结是由单层石墨烯和单层h-BN通过范德瓦尔斯力结合而成。异质结的顶层为石墨烯层，底层为h-BN层，这种堆叠最稳定的构型为Bernal堆叠，即一个碳原子位于硼原子顶部，而另一个碳原子位于h-BN环中心上方如图1(a)。GBNNR中y为周期性方向，x为开边界方向包含40条锯齿链(即N = 80个格点)。局域势V₁和V₂分别施加在石墨烯层的上半部分和下半部分，V₃和V₄分别施加在h-BN层的上半部分和下半部分如图1(b)。

石墨烯和h-BN层相应的面内最近邻跃迁能分别为 $t_{\text{1}}=2.64\text{eV}$ 和 $t_2=2.79\text{eV}$ ，最近邻层间跃迁能为 $t_{\perp }=0.43\text{eV}$ ，在位能 $E_B=3.34\text{eV}$ 和 $E_N=-1.40\text{eV}$ ，如表1所示。

Table 1. Tight-binding parameters of graphene/h-BN heterojunctions calculated by density functional theory [28] (Unit of Energy: eV)

表1. 密度泛函理论计算出的石墨烯/h-BN异质结紧束缚参数 [28] (能量单位：eV)

基于以上参数我们采用紧束缚模型计算了Bernal堆叠GBNNR能带结构，如图2(a)所示。当局域势和ISOC均为零时，能带图中有一个有微小带隙，大小约为 ${\Delta }_E=20.5\text{meV}$ ，这表明GBNNR是半导体。改变石墨烯层局域势V₁，使其从零开始增大，发现带隙 ${\Delta }_E$ 先增大后减小如图2(b)中红色曲线，系统始终处于半导体相。

Figure 2. (a) Band structure of GBNNR with zero local potential and ISOC; (b) Band gap of GBNNR Δ_E as a function of V₁

图2. (a) 局域势和ISOC均为零时GBNNR能带结构；(b) GBNNR带隙Δ_E随局域势V₁变化图

在GBNNR中，考虑到石墨烯层ISOC (自旋轨道耦合强度 $\lambda =0.02t_1$ )，则系统由半导体相变为QSH相。通过改变局域势V₁，可以调控GBNNR的能带结构，如图3所示。单层石墨烯中ISOC打开的体间隙为 ${\Delta }_g=2{\Delta }_{SO}=0.55\text{eV}$ ，其中 ${\Delta }_{SO}=3\sqrt{3}\lambda$ ，而在GBNNR 中，当V₁=0时，由于异质结受到层间相互作用的影响，ISOC打开的体间隙实际上为 ${\Delta }_g=2{\Delta }_{SO}-\delta$ ，大小约为0.49eV。图3(a)中，当费米能 $\left|E_f\right|<{\Delta }_g$ 时，系统处于QSH态，此时自旋陈数为非零值，Cs = 1。

如图3(c)，当 $V_{\text{1}}=0.2\text{eV}$ 时，系统依然具有较大的体间隙 ${\Delta }_g=0.44\text{eV}$ ，局域势将体间隙中自旋向上和自旋向下的能带左右分开。费米能 $\left|E_f\right|<{\Delta }_g$ 范围内与能带有四个交点，用A，B，C，D来表示，我们计算了这四个点对应的状态概率分布如图3(d)所示，A和B表示自旋向上的状态，分别分布在石墨烯层的两个边缘处，即N = 80和1附近，C和D表示自旋向下的状态，分别分布在石墨烯层的两个边缘处，即N = 1和80附近。根据 ${\nu }_y=\frac{\partial E_k}{\partial k_y}$ 得到这四个边界态在纳米带中的流向示意如图3(f)，A和D位于石墨烯层相同边界上，具有相反的自旋方向和传播方向，B和C与之类似，从而形成了螺旋性边缘态。这表明此时石墨烯层依然处于QSH态，自旋陈数Cs = 1。

Figure 3. Left column is the band structure of GBNNR. The red and black curves represent the spin-up and spin-down energy bands respectively. The parameters are: λ = 0.02t₁, (a) V₁ = 0 eV; (c) V₁ = 0.2 eV; (e) V₁ = 1.2 eV. (b) The blue curve is the bandgap of the heterojunction under periodic boundary conditions, i.e. the bulk gap; the red curve is the GBNNR bandgap. (d) and (f) are the schematic diagrams of the probability distributions of edge states and the edge state propagation in (c), respectively

图3. 左列为GBNNR能带图。红色和黑色曲线分别表示自旋向上和自旋向下的能带。参数为：λ = 0.02t₁，(a) V₁ = 0 eV；(c) V₁ = 0.2 eV；(e) V₁ = 1.2 eV。(b) 蓝色曲线为异质结在周期性边界条件下的带隙，即体间隙；红色曲线为GBNNR带隙。(d)和(f)分别是(c)中边缘态的波函数概率分布与流向示意图

继续增大局域势V₁，体间隙 ${\Delta }_g$ 不断减小，当 $V_1=0.3~0.9\text{eV}$ 时，体间隙与局域势几乎成线性关系： ${\Delta }_g\approx 2{\Delta }_{SO}-\delta -\frac{V_1}{4}$ ，当 $V_1\ge 1.0\text{eV}$ 时， ${\Delta }_g=0$ ，如图3(b)中蓝色曲线。而当 $V_1=0~0.5\text{eV}$ 时GBNNR带隙 ${\Delta }_E=0$ ，当 $V_1>0.5\text{eV}$ 时系统打开一个较小的带隙，如图3(b)中红色曲线。这说明在 $V_1=0~0.5\text{eV}$ 时，体带隙中会有无间隙的边缘态出现，此时GBNNR处于QSH相，当继续增大V₁时，无间隙的边缘态消失，系统从QSH相变成半导体相，如图3(e)。

Figure 4. (a) Band structure of GBNNR with λ = 0.02 t₁, V₄ = −3.665 eV. The red and black curves represent the spin-up and spin-down energy bands respectively. (b) and (c) are schematic diagrams of the probability distributions of edge states and the edge state propagation in (a), respectively

图4. (a) GBNNR能带图，红色和黑色曲线分别表示自旋向上和自旋向下的能带，λ = 0.02t₁，V₄ = −3.665 eV。(b)和(c)分别是(a)中边缘态的波函数概率分布与流向示意图

在GBNNR的能带中，低能区域的能带主要来源于石墨烯层，h-BN层的能带则主要分布在高能区域，接下来我们通过调节局域势V₄，能使h-BN层的能带向低能区域移动，从而发现了有趣的能带演化现象。如图4(a)所示，当局域势V₄ = −3.665eV时，能带结构中出现了由ISOC打开的带隙，在此带隙中间有无间隙的边缘态存在。通过计算费米能E_f与能带相交的A，B，C，D四点的波函数概率分布可以知道这些态均为边缘态，如图4(b)所示。A和D分别表示自旋向上和自旋向下的状态，它们的传播方向相反，但是都分布在石墨烯层的同一个边缘处，即N=1附近。B和C分别表示自旋向上和自旋向下的状态，它们的传播方向也相反，但是它们同时分布在石墨烯层和h-BN层的同一个边缘处，即N = 80和160附近。根据边缘态的分布情况，我们可以给出边缘态在异质结中的流向图，如图4(c)所示。图中自旋向上的电流(红色虚线)和自旋向下的电流(黑色虚线)都各自形成两个闭合回路，其中一个回路在石墨烯层，另一个回路在横跨石墨烯层与h-BN层。每一个回路中，边缘态都是螺旋型的。从波函数的概率分布值上可以发现，边缘态主要分布于石墨烯层与h-BN层的两个不同边缘处，即N = 1和160处。将两个回路合在一起，系统的边缘态依然是量子化的，对应的状态为QSH态，自旋陈数Cs = 1。但是此QSH态，与单层的QSH态不同，它是由层内和层间的边缘态共同构成的。

4. 结论

本文在晶格匹配的条件下，采用紧束缚模型研究了Bernal堆叠zigzag边界GBNNR在局域势调控下的能带结构及量子相变。我们发现当石墨烯层内禀自旋轨道耦合强度给定时，石墨烯层处于量子自旋霍尔态而六方氮化硼层处于绝缘态，此时增大石墨烯层局域势强度，系统会从量子自旋霍尔相变到半导体相。而调节h-BN层局域势可以使GBNNR产生新的无间隙边缘态，从而使系统中石墨烯层的量子自旋霍尔态变成由层内和层间的边缘态共同构成的量子自旋霍尔态。

致谢

本文受到国家自然科学基金资助No. 11864039。

NOTES

$$

^*通讯作者。