三聚化非厄密晶格中具有趋肤效应的拓扑边缘态*

许楠 张岩

(东北师范大学物理学院,长春 130024)

1 引 言

拓扑绝缘体是一种具有全新量子特性的物质态,在量子理论的发展过程中具有重要意义.从1980年Klitzing等[1,2]在二维电子气系统中发现了量子整数霍尔效应,到Haldane[3]发现了量子反常霍尔效应,以及之后的十年,拓扑性质相关的理论实验研究都得到了极大的发展[4-10].而作为其中最为经典的拓扑模型之一,一维的Aubry-André-Harper (AAH)模型在研究多聚化晶格方面得到了越来越多的关注与重视.该模型最初被用来研究类周期系统的局域性转变[11-15].AAH晶格的本征能谱由三条能带、两条能隙组成,两条能隙中分别存在可由参数调制的边缘态.因此,能谱可分为两种不同的相位：有边缘态的非平庸相和没有边缘态的平庸相.然而,在实际的拓扑实验中,由于总是与环境相互作用,所以需考虑耗散效应对系统的影响.其中,典型的例子便是针对非厄密系统(引入增益和损耗)的研究[16-27].尤其在近几年,均匀的增益或损耗被引入到一些非厄密的拓扑绝缘体理论或者实验中,如Su-Schrieffer-Heeger (SSH)模型[28,29]和AAH模型[30-33]的理论研究,以及一维耦合微腔[34,35]和一维耦合光波导的实验论证[36,37],都得到了很多厄密系统中没有的新奇现象.实验中通过设计一维耦合环腔,在虚数规范场中调制辅助腔的增益和损耗[38],使得二聚化模型系统中最近邻正反耦合系数出现差值,从而产生趋肤效应[39,40],该效应在近年获得极大关注.区别于同一耦合的正反向系数一致的拓扑系统具有一般性的体边对应关系,非厄密系统中通过改变晶格正反向耦合系数差,破坏系统的手性反演对称,从而展现了扩展的体边对应关系.主要表现为包含系统的拓扑保护边缘态和体态在内的所有本征态都局域于开放系统单侧边缘附近的趋肤效应.然而,三聚化非厄密系统的趋肤效应的研究还不充分.因此,本文旨在研究一维三聚化的非厄密晶格的具有趋肤效应的拓扑特性.将三聚化模型结构从一维链状拓展到类一维的锯齿形结构,通过适当调制胞内最近邻正反耦合系数差和次近邻耦合系数,观察三聚化模型中的拓扑特性变化.发现正反向系数差不为零的最近邻耦合对非厄密系统的能谱和拓扑保护边缘态都有显著影响,并引发新奇的趋肤效应；次近邻耦合则对系统能谱的上下能隙中的边缘态能带的宽度和系统本征态的局域效果有影响.最后,进行了总结并展望了拥有丰富结构模型的各种新型拓扑材料的发展前景.

2 理论模型与能谱分布

如图1所示,本文考虑的是类一维的非厄密三聚化晶格,并为实现可调的耦合,将晶胞结构设计为锯齿形.晶胞内最近邻正向耦合系数为g1+µ,反向系数为g1-µ,其中 2µ为最近邻正反向耦合系数差,次近邻耦合系数为ν,胞间耦合系数为g2.在该非厄密系统中,在每个晶胞内添加均匀的增益损耗——添加势能 iγ于增益格点A、零势能于中性格点B、势能 -iγ于损耗格点C.这些参数条件构成了本研究所需的扩展的一维非厄密的AAH模型.其中,令an,bn和cn表示第n个晶胞内格点的湮灭算符,总晶胞数为N.则该扩展的AAH模型的哈密顿量为

式内各项分别为AAH模型哈密顿量HAAH=次近邻项的哈密顿量和非厄密项的哈密顿量然后将哈密顿量代入薛定谔方程[41],便可得：

图1 扩展的三聚化非厄密晶格示意图Fig.1.Schematic of the generalized trimerized non-Hermitian lattice.

结合(1)式与(2)式,得到可描述该三聚化模型的本征能谱,如图2所示；并可对与最近邻耦合有关的趋肤效应和与次近邻耦合有关的拓扑保护边缘态进行深入的研究,如图3和图4所示.这里设相关参数为：N=30,g1=1,γ=0.2,g2≡g1-δcosθ(δ=0.5,θ∈{0,2π} ).

首先,通过调节胞内最近邻正反向耦合系数差,对系统本征能谱进行分析.图2(a),图2(c),图2(e)和图2(g)为系统本征能谱的实部,能带中包含上、下两个非平庸能隙,共四条边缘态.上(下)能隙中存在两条边缘态的本征值,分别用棕色虚线和绿色实线表示,且这两条边缘态简并.图2(b),图2(d),图2(f)和图2(h)为本征能谱的虚部,同样用绿色实线和棕色虚线表示对应边缘态的虚部.系统本征值的实部描述的是能谱中系统能带的位置；本征值的虚部描述的是能谱中能级的宽度.

此处,不考虑次近邻耦合(ν=0 ),体态的能谱具有对称特征.如图2(a)和图2(b)所示,当胞内最近邻正反向耦合相同(µ=0 )时,体态和边缘态的本征能谱均关于E=0 和θ=π 对称,能谱的三条能带在 0.5π 和 1.5π 附近区域结合,上下能隙中的四条边缘态具有相同长度.然而,如图2(c)—图2(f)所示,随着µ逐渐增大,虽然体态能谱仍保持原有对称性,但是用棕色虚线表示边缘态能谱的实部和虚部均会缩短,而用绿色表示的两条边缘态的实部和虚部均不变化.当最近邻正反向耦合差继续增大至µ＞0.2 时,如图2(g)和图2(h)所示,棕色虚线表示的两条边缘态消失,只存在绿色实线表示的两条边缘态.这表明边缘态的简并情况发生改变,边缘态数量由四条缩减到两条,该三聚化晶格产生了新的拓扑特性.这对应边缘态在晶格的分布情况就是晶格一侧的两条边缘态消失,如图3所示,此为边缘态的趋肤效应.

3 体态与边缘态的趋肤效应

在上一部分中,通过对动力学方程的求解,讨论了本征能谱在最近邻耦合影响下边缘态能带宽度以及数量发生的变化.本小节通过观察本征态的演变,继续讨论三聚化非厄密系统中的边缘态和体态的趋肤效应.图3为该三聚化模型包括体态(蓝点)和拓扑保护边缘态(绿线)在内的所有本征态的分布.如图3(a)所示,当µ=0 时,体态均匀分布,边缘态则对称地局域在晶格的两端.如图3(b)所示,随着µ增大,所有的本征态都开始向晶格的同一端附近局域,展现了趋肤效应.如图3(c)所示,当µ=0.2 时,本征态的局域性更好；体态继续向一端局域；边缘态完全消失,只存在于左端,数量也由四条变为两条.也就是说实现了边缘态数量的动态调节.如图3(d)所示,当最近邻正反耦合差继续增大,体态也完全局域在左端,也就是说三聚化非厄密晶格本征态的趋肤效应愈加显著.

图2 三聚化非厄密晶格的本征能谱 (a),(c),(e)和(g)分别为 µ=0,0.1,0.15,0.2 时本征能谱的实部；(b),(d),(f)和(h)分别为 µ=0,0.1,0.15,0.2 时本征能谱的虚部,边缘态能带用绿色实线和棕色虚线表示Fig.2.The eigen-energy spectrum of the trimerized non-Hermitian lattice.(a),(c),(e)and (g)show the real parts of the energy spectrum for µ=0,0.1,0.15,0.2,respectively；(b),(d),(f)and (h)show the imaginary parts of the energy spectrum forµ=0,0.1,0.15,0.2.Green solid lines and brown dash lines represent the energy band of the edge states.

图3 系统本征态的光子分布 (a)µ=0 ；(b)µ=0.1 ；(c)µ=0.2 ；(d)µ=0.3Fig.3.Photon distributions for eigenstates of the system：(a)µ=0 ；(b)µ=0.1 ；(c)µ=0.2 ；(d)µ=0.3 .

接下来,在趋肤效应下,继续观察次近邻耦合对本征能谱和边缘态的影响.在图4中,第一列描述的是受次近邻耦合影响的本征能谱的实部,第二列代表上能隙中边缘态的光子分布,第三列代表下能隙中边缘态的光子分布.当µ=0.5 时,次近邻耦合系数在适当增加,如图4(a),图4(d)和图4(g)所示,能谱的实部仅关于θ=π 对称,而不再关于E=0对称.值得注意的是,非平庸的上能隙开始变窄,导致非平庸拓扑区域缩小；非平庸的下能隙开始变宽,导致非平庸相位区域扩大.通过比较图4(b)和图4(c)发现,相对应的边缘态不同于体态的普通振荡,而是呈急速衰减.随着上下能隙中的边缘态能带数量的缩减,只在上下能隙中各存在一条.如图4(b),图4(e)和图4(h)所示,随着次近邻耦合的增强,由于趋肤效应而局域在左侧的上能隙边缘态的峰值在逐渐减小,局域性减弱；而由于趋肤效应而局域在左侧的下能隙边缘态变化正相反,峰值逐渐增大,局域效果逐渐增强.此外,如果次近邻耦合进一步增强,本征能谱上下能隙及其边缘态会继续呈现相反的变化趋势；最终上边缘态能带会缩短直至消失,下边缘态峰值逐渐增大,局域效果会明显增强.

4 结论与展望

本文探讨和比较了扩展的非厄密三聚化晶格的能谱和边缘态的演化和局域效果.首先,特殊调制晶胞内的最近邻耦合系数,即最近邻正反向耦合存在差异,会导致系统本征态发生趋肤效应.包括体态和边缘态在内的所有本征态都局域于系统的同一端附近；并且,边缘态的二重简并态消失,边缘态数量由四条减少为两条.其次,次近邻耦合对系统能谱的演化有着显著的影响.次近邻耦合会使系统非平庸能隙和相对应的具有趋肤效应的边缘态发生不对称变化,且会影响所有本征态的局域效果.未来,通过设计多样的多聚化非厄密晶格,并应用于拓扑量子材料的研究中,一定会发现更多有意义的、新奇的拓扑物理性质.

特别感谢东北师范大学物理学院吴金辉教授对文章写作的建议.

图4 次近邻耦合影响下系统的本征能谱和边缘态的光子分布 (a),(d)和(g)分别为 ν=0.1,0.15,0.2 时本征能谱的实部；(b),(e)和(h)分别为 ν=0.1,0.15,0.2 时上能隙中的边缘态；(c),(f)和(i)分别为 ν=0.1,0.15,0.2 时下能隙的边缘态Fig.4.The eigen-energy spectrum and photon distributions for edge states of the system under the influence of the next-nearestneighbor coupling：(a),(d)and (g)show the real parts of the energy spectrum for ν=0.1,0.15,0.2 ；(b),(e)and (h)show the edge states whose energy band are in the upper gaps for ν=0.1,0.15,0.2 ；(c),(f)and (i)show the edge states whose energy bands are in the lower gaps for ν=0.1,0.15,0.2 .