罗开发 余睿

(武汉大学,物理科学与技术学院,武汉 430072)

利用凝聚态物理中紧束缚哈密顿量与集中参数电子线路中基尔霍夫方程的对应关系,可以在电子线路中设计出种类丰富的拓扑物态.本文详细介绍用电路实现一维SSH模型、三维结线半金属模型和外尔半金属模型的设计方案.在上述拓扑电路中可探测到端点态、表面鼓膜态、表面费米弧等体拓扑性质对应的界面态.由于电子线路对应的紧束缚哈密顿量中的跃迁项具有丰富的调控自由度,如强度、距离、维度等,容易推广到非厄密系统以及四维或更高维度的系统,使得人们能在电路中设计和验证传统凝聚态体系中难以实现或无法实现的新物态.此外,电子线路具备器件功能多样、制备工艺成熟可靠等优势,为探索新奇物态提供了一个便利的实验平台.

1 引 言

利用拓扑的概念对物理系统中的物态进行分类[1—3]在近十多年来取得了巨大成就.整数、分数量子霍尔效应、量子自旋霍尔效应、量子反常霍尔效应、狄拉克半金属、外尔半金属和拓扑结线半金属等拓扑态极大地促进了人们对电子材料物性的理解[4—15].由于拓扑态的出现主要取决于能带结构中的贝里曲率,于是拓扑能带论被迅速推广到冷原子体系[16—21]、光子[22—28]、声子[29—34]和机械系统[35—46]等经典体系.

近期,在经典线性电子线路中[47—60]研究人员设计出了多种拓扑态,包括:基于时间反演对称的Hofstadter模型[47,48]和Kane-Mele模型[61]的二维拓扑电路;一维和二维SSH拓扑电路[49,59];存在外尔点的三维拓扑电路[51,52];借助负电阻引入等效磁场后具有手性边界态的二维陈拓扑电路[62,63];模拟高阶拓扑角态[64,65];拓扑安德森绝缘态[66]和马约拉纳零能模[67]的电路.这些电路的主要构成单元是理想电容和电感,器件的连接方式和参数决定了系统的对称性和物理性质.上述工作采用的分析方法大多是电路理论中常用的电势运动方程[47,48]或拉普拉斯矩阵方程[51],与我们研究电子结构拓扑态时熟悉的哈密顿量有一定距离.为便于在哈密顿力学的框架下对电路系统的拓扑态进行设计和分析,我们介绍描述集中参数电子线路的基本方程-基尔霍夫方程与凝聚态物理中紧束缚哈密顿量之间的对应关系[68],从而将电子线路的问题转化成研究单粒子系统紧束缚哈密顿量的问题.本文主要介绍一维SSH电路系统、三维结线和外尔电路系统的设计方案[52].基于此,我们可以将凝聚态材料中已有的拓扑态推广到电路系统,也可以设计出传统体系中难以实现的一些新奇拓扑态,比如非厄密拓扑态,四维或更高维度的拓扑态等.

2 从基尔霍夫电流方程到紧束缚哈密顿量

集中参数电路中电流和电压的变化规律由基尔霍夫电流和电压定律刻画.前者源于载流子的电荷守恒定律,后者源于能量守恒定律.对电路网络中的节点进行编号α=1,2,···,N,任意两节点α和β间导纳记为yαβ.节点α上的电势记为vα(取地面为零势能参考点),流入该节点的净电流记为Iα.对节点α,从与之相连的节点α′流入的电流记为Iαα′(α′=α).欧姆定律给出Iαα′=(vα′-vα)yαα′,对所有支路求和得到yαα′=Iα.重复此过程,可得电路中所有节点上电压和流入电流之间的关系式LV=I,其中V=(v1,v2,···,vN)T为节点电压矢量,I=(I1,I2,···,IN)T为节点电流矢量.矩阵L包含了电路网络中的所有器件连接信息[69].

考虑一种简单情况:电路由一系列接地的LC谐振回路构成,LC谐振回路之间由电容器连接(见图1).后面我们会看到电路采取这种连接方式形成谐振频率色散的机制和晶体中电子形成能带的物理图像非常相似.在交流信号下电容C和电感L的导纳分别为yC=jωC,yL=1/(jωL),其中为虚数单位,ω为电路的谐振频率.如果线路不外接任何源和漏(实验上可以通过为电路输入脉冲信号而后撤去电源来实现),则电流矢量I=0,矩阵方程LV=0化简移项后可写成

上式与定态薛定谔方程H|ψ〉=E|ψ〉形式一致.这里谐振频率1/ω2L为本征值,节点电压矢量V为本征矢量,Y可看成哈密顿量矩阵.Y中的非对角元yαβ=-Cαβ可理解为准粒子在节点α和β之间的跃迁强度,正比于电容值;对角元yα为准粒子在轨道α上的在位能.由此,Y矩阵可类比为凝聚态物理中常见的紧束缚哈密顿量.求解此本征值方程,可得电路谐振频率构成的频谱.回想电子能带的形成过程:对于孤立原子,其中的电子具有离散能级.随着原子间距减小并逐渐形成晶体,电子可以在相邻原子间跃迁,孤立原子的离散能级展宽成固体中的电子能带.在电路中,孤立的LC谐振回路具有特定的谐振频率可类比于原子上的电子能级,当这些回路通过电容器C连接后,电容值的大小决定LC谐振回路间的耦合强度.当耦合足够强时,所有LC回路的谐振频率会形成由谐振频率构成的能带.因此,通过调节LC谐振回路的连接方式及强度,我们可以调控谐振频率色散,使其发生频带的反带或简并,实现拓扑绝缘态或者拓扑半金属态.这正是我们在电子线路中模拟SSH模型、结线和外尔态的物理基础.

图1 一维SSH电路.原胞(蓝色虚线框)内有A和B两个不等价节点,经并联的电感L和电容 C0 接 地.原胞内节点由电容C1相连,原胞间节点由电容 C2 相连Fig.1.The 1 D LC chain,in which a unit cell containing two inequivalent nodes A and B labeled by a dashed blue box.Each node A or B is grounded through a parallel connected inductor L and capacitor C0.All nodes are connected by C 1 and C 2 alternatively.

3 一维SSH电路中的拓扑态

3.1 一维SSH电路的模型哈密顿量

Su-Schrieffer-Heeger (SSH)模型是最简单的一维拓扑模型,在冷原子[70]、光子晶体[71]、声子晶体[72,73]和机械系统[36]中均有详细的理论和实验研究.在电子线路里,已有工作从系统阻抗频谱[51]的角度对一维LC链的拓扑端点态进行了考察.这里我们采用哈密顿矩阵的形式对这个典型的拓扑体系进行重新刻画和分析,展示这一方法的有效性和直观性.

图1为电路SSH链的示意图.原胞内包含两个不等价节点A和B,通过LC回路接地 (电感L,电容 C0).取原胞间距a为单位长度.原胞内的AB节点通过电容器 C1连接,原胞间的A-B节点通过电容器 C2连接.对于第n个原胞,根据基尔霍夫定律可写出流入原胞A,B节点的净电流分别为:

线路无源和漏接入时 IA(n)=IB(n)=0,两端同时除以 j ω 得到

其中 µ ≡C0+C1+C2.联立所有原胞内节点上的方程可得到形如 Y V=(ω2L)-1V 的矩阵方程,从Y中,可以提取出紧束缚模型中的跃迁项和在位能项:

其中R表示跃迁矢量.考虑到系统具有平移对称性,通过傅里叶变换(m,n=A,B)可得动量空间中的布洛赫哈密顿矩阵,

其中σ+=(σ1+iσ2)/2,σ0为单位阵,σ1,2,3为 泡利矩阵.这正是凝聚态物理中为人熟知的SSH哈密顿量,其色散关系为Ek=µ±|h+|=µ±

图2 (a) 上: 电容 C1 从 0 逐渐增加至超过 C2,频谱中两个端点态 (红线) 在能 C1 =C2 处消失,表明系统发生了拓扑相变.下:系统缠绕数从1到0的跃变与端点态的消失临界值一致.C1/C2=0.5 时 系统的等效极化矢量在 dx -dy 平面上随动量参数k从0连续变到2π时绕原点一周.红圈对应缠绕数为1;C1/C2=1.5 时 (dx,dy)绕原点0圈.黑色圈对应缠绕数为0.(b)取C1/C2=0.8,端点态(绿色和紫色)以及一个随机挑选的体态(灰色)对应的电势分布VFig.2.(a) Upper: increase the parameter of C1 from zero to exceed C 2,the end states (red) converge into the bulk states,indicating the topological transition.Bottom: the transition of winding number is consistent with the appearance and absence of end states.The effective polarization vector (dx,dy) winds the original a round when the momentum varies continuously from 0 to2π for C1=0.5 (left),while zero round for C1=1.5 (right).(b) The electric potential distributions of two end states (green and puple) and a randomly selected bulk state (grey).

3.2 系统的拓扑数

由于 HSSH中 µ σ0项不改变系统的拓扑性质,在下面的讨论中均将其忽略.系统具有手征对称性,其拓扑不变量由缠绕数来描述:

其中Θ(x)是Heaviside阶梯函数,当x＞0时Θ(x)=1,而x＜0时Θ (x)=0.上述结果可以通过赝自旋σ=(σx,σy)在赝磁场d(k)=(dx(k),dy(k))(dx=(h++h-)/2,dy=i(h+-h-)/2 )中的极化来理解((见图2)).赝磁场 d (k)随参数k的演化会导致赝自旋相应发生旋转,绕原点 d=0 的圈数对应缠绕数W.我们还可以用Zak相位来表征系统是否处于拓扑非平庸态,其 中是贝里联络,是一维SSH模型频谱中低频支对应的布洛赫电势分布函数.由此可得系统的Zak相位为

与缠绕数的计算一致.实验上,缠绕数和Zak相位可以通过观测拓扑泵浦[74]和表面阻抗[72]等物理量来测量.

3.3 受拓扑保护的端点态

从上面的讨论已经知道,|C1|＜|C2| 时系统处于缠绕数为1的拓扑非平庸态,存在受手征对称性保护的零能解.对于含 N原胞的 SSH链,H=记其零能模则方程H|0〉=0 给出如下系数间的关系式:

和边界条件

由 C1=0 可 推出 b1=aN=0,说明对有限长度的SSH链,(8)式只有平庸解.但系统趋于热力学极限 N → ∞ 时,可以给出a1,bN为非零有限值的两个非平庸解.通过测量A或B格点上边界态从端点往体内的电势分布衰减指数,可以与理论值 - C1/C2进行对比.

4 三维电路中的拓扑结线和外尔态

4.1 拓扑结线态

在电子系统中,结线态[75]可以视为其它拓扑态的母态.通过改变系统参数,如自旋轨道耦合强度和应力等,结线态可以演化成外尔态、狄拉克半金属态或拓扑绝缘态等.在电路中,我们首先考虑实现拓扑结线态,然后通过调控对称性,实现外尔态.

对于无自旋系统,一般有三种方案可以设计出拓扑结线态: 1) 时间反演加空间反演[76,77];2) 镜面对称性[78,79];3) 包含螺旋轴或滑移面的非简单空间对称群[80—82].由于无洛伦兹力和科里奥利力的系统中天然存在时间反演对称性,而且系统空间对称性越高相应地会对电路网络中的器件参数值的选取给出更严苛的要求,因此我们采用上述第一种方案.考虑无自旋的两能级系统,k参数空间中哈密顿量的一般形式为时间反演算符为T=K,即复共轭算符.空间反演算符I可表示为σ0或σ3,取决于两个基函数宇称相同或相反.当系统同时具有时间反演和空间反演对称性时,哈密顿量将受到如下约束:

时间反演要求d1,3(k)为偶函数,而d2(k)为奇函数.当I=σ0时,空间反演不产生约束,而I=σ3时,空间反演对称性要求d1,2(k)为奇函数,d3(k)为偶函数.综合上述约束方程得到d1(k)=0,d2(k)为奇函数,d3(k)为偶函数.结线态的出现要求d1(k)=d2(k)=d3(k)=0.在三维k空间中求解两个约束方程,可以得到无穷多个解,刚好是我们寻找的结线态.例如这两个函数同时为零的解是k空间中平面与二次曲面的交线,此交线即为结线型能谱简并点.下面我们按照这个思路在电路中设计结线态.

我们从二维六角蜂窝电路出发,图3(a)所示,a-b平面内所有A和B节点通过电容C1,2,3连接(这里取三种电容值相等).每个节点A (B)由并联的电感LA(LB)和电容CGA(CGB)接地.按照上一节介绍的方法,这些二维电路系统的基尔霍夫电流方程可写成与石墨烯哈密顿量完全一致的形式.类似于石墨烯的特征能带,二维六角蜂窝电路的共振频谱在第一布里渊区内存在两个闭合点.沿c方向堆叠,层间无耦合时二维布里渊区中的频谱交点将在三维布里渊区内形成图3(a)中小图所示的两条直线型结线.如果将最近邻层间的A和B节点用电容 C4连接(图3(b)),调节 C4的容值并不破坏系统的空间反演对称性,但直线状的结线可逐渐弯曲至闭合成为我们感兴趣的结环.如果用 CA(CB)连接层间的A-A(B-B)节点(图3(c)),选择参数CA=CB和CGA=CGB,此时空间反演对称性将被破缺,结线型交叉有可能退化成离散的外尔点.上述电路的等效紧束缚哈密顿矩阵为[52]:

图3 三维LC电路示意图.(a) 单层LC蜂窝电路沿c方向上无连接的堆叠起来,每个原胞内的两个不等价节点A和B在层内由 C1 ,2,3连接.每个节点A (B)都通过并联的 LA (LB) 和CGA(CGB) 回路接地.a,b和c表示格矢.小图: 单层LC电路频谱中包含两个简并点,沿着 kc 方向扩展将在布里渊区中形成两条直线状的结线(红线).(b) 电容 C4 连接最近邻层间的节点A和B,结线在给定合适的 C4时将弯曲成闭合环形.(c) 电容 CA (CB) 连 接最近邻层间的A-A(B-B)节点对,CA=CB 且 CGA=CGB 时空间反演对称破缺,环状结线可能退化成离散的外尔点.此外,LC网络可以变形成(d),简化实验装置的同时保证频谱不变Fig.3.Schematic setup of the 3 D LC circuit lattice.(a) LC honeycomb layers stacked along c-direction without any interlayer connection.Two inequivalent nodes A and B within a unit cell,linked by capacitors C 1,2,3.Each node A(B) is grounded through the parallel connected inductor LA (LB ) and capacitor C GA (C GB ).a,b,and c denote lattice vectors.Inset: spectra of a single layer LC lattice includes two band-crossing points,extended uniformly along kc-direction and form two straight nodal lines (red) in the BZ.(b) Connecting nodes A and B between the nearest neighbor-layers with C4.The straight lines could be curved to a closed ring given appropriate C4.(c) Connect node-pairs A-A(B-B) between the nearest neighbor-layers with CA (CB ),to break space inversion symmetry by setting CA=CB and CGA=CGB,and the nodal ring may be degenerated to discrete Weyl points.The LC network can be deformed into (d),which is convenient to construct circuit elements in experiments while spectrum invariant.

体边对应要求非平庸的结线对应着表面态的存在,沿(001)方向投影能带的结环内部会出现图4(b)中色散极小的鼓膜态.将三维布洛赫哈密顿矩阵参数化为等效的一维系统 H (ka,kb),沿kc方向进行积分即得到 ka-kb面上的贝里相位为频谱中低频支的布洛赫电势分布函数.积分路径位于结线内部时得到拓扑非平庸的贝里相位π,否则是平庸的相位0.

图4 (a) 结线(红色)及其在(001),(010)和(100)面上的投影(灰色).此时参数取为C1=1mF,C2=2mF,C3=1mF,C4=0.833mF,CGA=CGB=1mF 和LA=LB=1mH.(b) 沿 Γ-M-A-Y-B-M-Γ路径的频谱,其中A和B是(a)中结线与kc=0平面的交点.上: (001) 方向的体态 (灰色)及表面态 (红线),以及周期边界条件下的能带 (紫色).下: 积分路径在结线内部(外部)的贝里相位 θ k ‖ 等于拓扑非平庸的 π(平庸的0)Fig.4.(a) Nodal line (red) and its projections (grey) on the (001),(010),and (100) planes.The parameters are set as C1=1mF,C2=2mF,C3=1mF,C4=0.833mF,CGA=CGB=1mF,and LA=LB=1mH.(b) Bands along Γ-M-A-Y-B-MΓ,where A and B are two points with kc=0 on the nodal line as labeled in (a).Upper: the bulk bands(grey) with the drumheadlike surface states nestled inside the projection of the nodal ring (red) on the (001) surface and the two bands (purple) in periodic boundary condition.Bottom: the Berry phase θ k ‖ equals π(0) inside (outside) the nodal ring.

图5 (a)布里渊区中的四个外尔点及它们在(001),(010)和(100)方向的投影.取CA=0.2mF,CB=0.01mF 和CGA=0.77mF,其它参数均与图4中相同.外尔点是由 d1 ,2(k)=0 决 定的结线与 d3 (k)=0 决定的两平面的交点,它们的手性用蓝色五角星(χ=+1 ) 和红色圆点(χ=-1 )标记.(b) 左上: 在 (001) 表面上,费米弧连接了手性相反外尔点的投影点.右上:开边界时路径AB上(亮青色虚线)的无能隙表面态,绿色虚线标记外尔点所在频率.左下: 垂直于 k a 的 各平面上的陈数.沿ka方向移动,经过正(负)手性外尔点时平面上的陈数会增加1(减少1)Fig.5.(a) Four Weyl points in the Brillouin zone and their projections on (001),(010) and (100) direction.CA=0.2mF,CB=0.01mF,and CGA=0.77mF are used in the calculations.The other parameters are the same as Fig.4.The Weyl points are the intersection points between the nodal ring determined by d1 ,2(k)=0 and the two planes determined by d3 (k)=0.The chirality are indicated as blue stars (red points) for χ=+1 (- 1 ).(b) upper left: on the (001) surface,Fermi arcs connect the projections of the bulk Weyl nodes carrying opposite chiralities onto the surface.Upper right: the gapless surface band in theA-B path.Dashed green line denotes the frequency where Weyl points lie.Bottom left: the Chern numbers for planes perpendicular to ka.Moving along k a,the Chern number increases (decreases) when the plane passing through the Weyl points with+1(—1) chirality.

4.2 拓扑外尔态

取电容值CA=CB,CGA=CGB可以破缺空间反演对称性从而使结线态相变到外尔态.此时d3(k)不再恒为零,能隙闭合条件需要进一步要求coskc=1+(CGA-CGB)/2(CA-CB).选取电容参数使得等式右边小于1,d3(k)=0在三维布里渊区中给出两个垂直于kc的平面,当这两个平面与d1,2(k)=0给出的结线相交时便得到图5(a)中的4个交点,这正是我们寻找的外尔点.由于时间反演对称性会把一个位于k的外尔点映射到手性相同但位于-k的外尔点,而整个系统的外尔点手性之和恒为零.于是这四个外尔点的手性必然有2个为正,另2个为负,即这类具有时间反演对称性的系统中外尔点的最小数目为4.如图5(b)所示将ka作为参数,可得一系列由kb-kc张成的二维面,由这些二维面上陈数的变化可判断出外尔点的手性.在远离外尔点处二维面上的能带具有能隙,若存在非零陈数,则开边界时能隙中会出现非平庸的手征表面态.因此,和电子系统中的外尔半金属类似,电路系统中手性相反的两个外尔点在表面的投影也会被表面费米弧连接起来.

在上面的讨论中,我们考虑的是理想电容和电感.实际中的电子元件受到制造工艺等因素影响,其真实参数在标注值的基础上都存在一定范围的误差,需要进一步考虑参数误差对前面两类拓扑态的影响.由于结线态需要空间反演对称性保护,在大量器件都存在误差的前提下要保持此对称过于困难,因而一般是不稳定的.而外尔态只要求系统的平移对称性以保持动量空间的存在,不需要其他对称性,可以期待外尔点能抵抗一定程度的器件参数误差.

下面我们对上述电路进行器件误差分析.考虑一个 3 ×3×3 的超胞.其中电容参数在精确值基础上附加一定幅度的随机误差,即哈密顿矩阵中的化学势和跃迁项被附加随机数.误差幅度从0开始以精确参数值的5%为步长逐渐增加到100%,每个数据点随机选取均匀分布的误差重复100次并计算第 N/2 条 和第 N/2+1 条能带的频率差.如图6所示,误差强度小于30%时外尔点依然存在.

图6 3×3×3超胞的能隙随误差幅度的变化,其中N是能带条数.参数同图5理想情形.每种误差幅度随机重复计算100次.数值计算中误差幅度不超过30%时外尔点仍然存在.小图为电容误差幅度为 ± 20% 时两个外尔点所在路径的色散Fig.6.The band gap,as a function of the tolerance values for a 3 ×3×3 super cell,where N is the number of total bands.The parameters are the same as Fig.5 in ideal case.Each range of tolerance is calculated 100 times.Numerical results show that Weyl points survive for ranges less than the critical value 30%.Inset: the bands along the k-path crossing two Weyl points with 20% range of tolerance on the capacitors.

5 小结与展望

本文介绍了一类基尔霍夫方程可映射为紧束缚哈密顿量的LC网络,将经典电子线路与单粒子的量子问题联系起来.我们以SSH模型为例,展示了如何在电路中实现一维拓扑态.通过堆叠二维蜂窝LC网络,我们设计了一个在不同参数下分别具有拓扑结线态和外尔点的三维LC网络,调控电容能够实现两者之间的拓扑相变,在此电路网格中,可以观测到结线态对应的表面鼓膜态和外尔态对应的表面费米弧,

从上述例子可见,电路模型哈密顿量中跃迁项的设计自由度极大,跃迁强度、方向、距离等均能通过连接电容器件来调控,借以实现传统凝聚态体系中不易制备的新奇物态,例如最近得到广泛关注的非厄密拓扑态.非厄密物理的典型特性之一是参数空间中的非厄密能带闭合点[83],同时伴随着体费米弧[84]和反常边界态[85]等行为.非厄密系统的物理性质强烈依赖于边界条件[86—88],而在通常的量子体系和人工超材料体系中难以在周期边界条件和开放边界条件之间自由切换.但对于集中参数电子线路,可以非常方便的调控部分线路的开或断来改变边界条件.此外,在系统中添加运算放大器或电阻等有源/耗器件后,很容易为电路网络引入非厄密效应.通过引入不对称跃迁项[53,89],基于本文讨论的蜂窝电路,可以展示多种以非厄密简并点为端点或边界的体费米弧和体鼓膜态,该系统中体边对应的破缺和恢复也有望通过实验观测得到验证.

电子线路作为一个宏观稳定易操控的物理平台,其潜力远不止于用来模拟其它量子或经典系统中已被实验观测的物态.考虑到集中参数电路中元件的相互连接几乎不受空间维度和器件数量的限制,按照类似的思路不难设计出对应于准晶[90]和非晶[91]哈密顿量的电路网络.通过更复杂的连接,还可以设计四维拓扑绝缘态[56]和更多其他系统中难以研究的高维拓扑态.