蔡潇潇 罗国语 李志强 贺言

(四川大学物理学院,成都 610065)

理论研究了转角双层石墨烯在施加不同单轴应变下的能带结构和光电导率,用连续模型分别计算了转角为1.05°和1.47°的转角双层石墨烯在应变下的能带、态密度以及光电导率,发现这些量随应变的变化是连续且显著的.通过对能带的分析以及光电导率的测量能够获得应变对平带产生的实际影响,这为今后实验对应变与平带的研究打下基础;此外样品往往受到具有空间不均匀性的应变作用,测量其局域的光电导率便能够估计应变的空间分布大小;同时应变对能带的调制为原位调控转角双层石墨烯的强关联、拓扑以及量子效应提供了思路.

1 引 言

转角双层石墨烯(twisted bilayer grapllene,TBG)是由两层石墨烯以一定转角堆叠在一起的二维材料,其转角带来的晶格失配使得表面形成莫尔条纹[1],这种莫尔周期性结构对TBG的电子能带结构有很大的影响,其中最为重要的是在零能点附近形成两条莫尔能带.当转角慢慢接近某些特定的角度(魔角)附近,这两条能带逐渐变得平缓形成平带,而在这样具有高态密度的平带中,电子之间的相互作用变得不可忽略,这使TBG 成为研究电子强关联态和超导态的重要体系[2-4].不仅如此,将公度堆叠的魔角TBG 放在六方氮化硼(hBN)基底上,由于彼此的晶格失配将破坏TBG 面内的C2对称性,使狄拉克点打开能隙,石墨烯每个谷被分开的两条平带将具有相反的陈数(±1)[5,6];实验发现在平带3/4 填充时,由于自旋极化效应TBG将表现出铁磁性[7],随后便在这一基础上发现了该体系中的量子反常霍尔效应[8].

无论是对于TBG 还是TBG-hBN 体系,平带对于这些丰富物理现象的产生都是至关重要的.除了用转角来调控平带,应变也是一种重要的方法,在样品的制备过程中,应变的产生几乎是不可避免的,比如基底带来的应变作用:当对两层石墨烯施加不同的应变时,其莫尔能带将会受到调制,原本非魔角下的TBG 在此应变下也能出现平带[9,10];同时,应变不仅能够使材料能带拓扑数发生改变[10],也能够使晶格产生重构,引发许多新奇的物理现象如孤子[11]、光子晶体[12]等;应变相对于转角更易调控,只需要具有压电性质的基底即可,这使得原位调控范德瓦耳斯材料中的强关联、拓扑以及量子效应成为可能.由此可见,研究应变对于TBG的影响是在理论和实验上都非常有意义的.

本文研究应变对TBG 这一重要体系的能带以及光电导率[13]的影响:光电导率作为实验研究上的重要物理量,不仅能够表征材料光学性质,也能够在一定程度上来揭示应变对于材料的电子能带结构的变化.之前已有理论工作研究了TBG 在应变作用下的能带结构[10]以及电子能谱[9,14],但是TBG 在应变作用下的光电导率还有待研究.有鉴于此,本文利用连续模型[10,15-17]理论计算了转角为1.05°(魔角)和1.47°的TBG 在其两层施加相反单轴应变时的能带结构、态密度以及低温下的光电导率:在应变的作用下,TBG的导带与价带的间隔增大,同时伴随着范霍夫奇点的位置的移动;另外,导带和价带仍然由两个狄拉克点相连,但是两个狄拉克点的位置发生了偏移,在能量上不再简并,这些变化都会对TBG的电子性质以及光学性质产生一定的影响.本文首先简单介绍TBG 在施加应变后的连续模型,然后对TBG 在应变下的能带、态密度以及低温下的光电导率进行讨论.

2 理论方法

图1 展示了转角为5°的TBG 晶格结构,其莫尔条纹周期清晰可见,随着TBG的转角减小,所形成的莫尔条纹周期将会远远大于晶格常数,导致其不同谷之间的相互作用可以忽略,这时可以用连续模型来计算TBG的能带结构.由连续模型给出的TBG的哈密顿量形式为[10,15-17]

图1 TBG 在转角为5°时的结构示意图,其莫尔周期结构清晰可见Fig.1.Schematic of TBG structure at θ=5°,the resulted moiré pattern can be clearly seen.

式中,h1(k),h2(k) 分别为两层单层石墨烯自身的哈密顿量;T(r) 代表由莫尔超晶格引起的具有空间分布的层间隧穿[16],其表达式为

由于应变和转角的作用,单层石墨烯的哈密顿量hl(k)的表达式变为[10]

式中,l=1,2 代表第一、第二层石墨烯;费米速度vF=106m/s,由此得到ℏvF/a=2.68eV;I为2×2的单位矩阵;ε代表了应变与转角的作用项,其表达式为

Dl,ξ为在应变和转角下的狄拉克费米子的位置[10],其形式为

其中Kξ为无转角和应变下的石墨烯两个谷的位置,β=3.14 为跳跃因子.

这里考虑对两层石墨烯分别施加大小相同方向相反的单轴应变与转角,即ε2=−ε1=,单轴应变S由幅度λ和方向φ决定,其形式为[10]

其中v=0.16 为石墨烯的泊松比.

利用连续模型得到的本征波函数再结合库伯公式可以计算出TBG的光电导率,其形式为[18]

其中,积分取遍所有本征态,S代表积分区域面积,f(ε)为费米分布函数,εα(εβ)和|α〉(|β〉) 分别代表本征值和本征态,vx=−(i/ℏ)[x,H] 为速度算符,η取3 meV 代表展宽.本文计算的所有光电导率值都是以单层石墨烯能带线性区域所贡献的光电导率值σmono为单位:

其中gs和gv分别为自旋与谷自由度,其大小为2.

有实验[14]表明TBG的单轴应变的幅度范围大小0%—7%,所以本文分别计算了当应变幅度λ为0%,0.3%,0.6%,应变方向φ为30°时,转角为1.05°和1.47° TBG的能带以及光电导率.

3 理论结果与讨论

图2(a)展示了TBG 在未考虑应变时的第一莫尔布里渊区示意图,由于施加的应变幅度较小,在应变下的第一莫尔布里渊区形状相比之前差别并不大.我们沿着“K1-Γ-K2-M2-Γ-M1-K1”的路径绘制了转角为1.05°和1.47°时TBG 在不同应变下的能带(图3).可见在转角为1.05°且未施加应变时,零能附近形成了两条平带(图3(a)),带宽约为17 meV,在施加应变后(图3(b)和图3(c)),两条平带间隔增大,变化量与施加的应变大小成正比;其他位置处能带之间的间距也随着施加应变的增加而增大.同样对于转角为1.47°的能带(图3(d)—(f)),除了本身未形成平带之外,应变带来的影响与1.05°时几乎相同.值得一提的是,对于魔角1.05°的TBG,其莫尔能带在应变的作用下依然保持平带的特征,这为平带的调控提供了额外的思路.另外,单轴应变的施加并未打破TBG的二重旋转与时间反演(C2Γ)对称性,这使连接导带和价带的两个狄拉克点依然存在,但TBG的面内镜面对称性受到破坏[19,20],狄拉克点的位置将发生偏移,偏移的大小与施加的应变有关,图2(b)展示了转角为1.05°、应变大小为6%时狄拉克点附近的能带,所沿路径为图2(a)中的虚线,此时狄拉克点位置由虚线路径中的黑点标注.

图2 (a) TBG 莫尔布里渊区示意图:两个大正六边形代表上下两层单层石墨烯的第一布里渊区,小正六边形为转角形成的莫尔布里渊区;(b) 偏移后的狄拉克点附近能带:红色和蓝色曲线分别代表图(a)中莫尔布里渊区同色虚线路径的能带,虚线路径中的黑点代表狄拉克点的位置Fig.2.(a) Schematic of TBG moiré Brillouin zone:The two large regular hexagons represent the first Brillouin zone of the upper and lower graphene layers,the small regular hexagons refer to the moiré Brillouin zone.(b) The band structures near the shifted Dirac points:Red and blue curve lines represent the band structures follow the same colored dashed lines path in the panel(a) respectively,the Dirac points are marked by the black dots in dash lines.

图3 TBG的能带 (a)—(c) TBG 在转角为1.05°,应变大小分别为0%,3%和6%时的能带;(d)—(f) TBG 在转角为1.47°,应变大小分别为0%,3%和6%时的能带Fig.3.Band structures of TBG:(a)-(c) The band structures with 0%,0.3%,0.6% uniaxial heterostrain at twisted angle θ=1.05°,respectively;(d)-(f) the band structures with 0%,0.3%,0.6% uniaxial heterostrain of at twisted angle θ=1.47°,respectively.

图4 展示了转角为1.05°的TBG 在施加0.6%应变后的能带(图4(a))、态密度(图4(b))以及低温下的光电导率(图4(c)),目的是为了分析能带、态密度与光电导率的对应关系.态密度是单位能量中态的个数,即能带越平,态密度越大,如图4(a)能带中的两个莫尔平带贡献了态密度(图4(b))中最大的两个范霍夫奇点,其他能带的鞍点也对应于态密度的小峰值.本文所计算的光电导率源于费米面附近的带内跃迁以及导带与价带的带间跃迁(直接跃迁),带内跃迁只主导零能附近的光电导率,因此对光电导率其他非零处吸收峰的贡献均来自具有相同动量的两个态之间的带间跃迁,也就是说影响非零能处光电导率的因素有两个:1) 相同动量处基态与激发态的态密度;2) 跃迁概率.图4(c)用绿色、红色与蓝色箭头分别表示光电导率的3 个主特征峰,对应于图4(a)中同色箭头所代表的带间跃迁,其他未标注的吸收峰同样也来自不同能带之间的带间跃迁:图4(c)中红色箭头表示的特征峰对应于莫尔价带对第一非莫尔导带之间的跃迁,贡献了光电导率中最大的峰值;蓝色箭头则对应跃迁能量在200 meV 之上最显著的吸收峰;绿色箭头表示的特征峰以及其附近的小峰均源于两个莫尔能带之间的跃迁,而这两个能带在相同动量处的态密度极大,却并未贡献出最大的吸收峰,这是因为两个莫尔能带之间的跃迁概率是很小的.特别地,当哈密顿量((1)式)具有电子空穴对称性时,价带与相对应的导带之间的带间跃迁概率为零[18],比如两个莫尔能带的M点之间并不能发生带间跃迁.

图4 TBG 在转角为1.05°、施加0.6%大小的应变时的能带(a)、态密度(b)以及光电导(c),图(c)中绿色、红色与蓝色箭头对应的吸收峰分别对应于图(a)中的同色箭头代表的带间跃迁Fig.4.(a) Band structure,(b) density of states and(c) corresponding optical conductivity of TBG with 0.6% uniaxial heterostrain at 1.05°.The green,red and blue arrows in panel(c) correspond to the interband transition marked with arrows of the same color in panel(a).

为了探究应变对于TBG 光学性质的影响,比较了TBG 在不同应变下的态密度(图5)以及光电导率实部(图6),两个图中的曲线从下到上分别代表施加的应变幅度为0%,3%,6%.在两个转角下,态密度极值峰的半高全宽随着应变的增大而逐渐增大(图5),对应于能带带宽的增大,同时也伴随着态密度峰值大小的减小—这些变化将使光学吸收峰变得扁平(图6).从图6 可见,TBG 光电导率的峰值集中在跃迁能量600 meV 以内,而在大于600 meV的区域,光电导值约为2 倍σmono.转角为1.05°时的主峰位于86 meV 能量处,1.47°时的主峰位于188 meV 能量处,在应变作用下,这两个峰的能量位置变化不大,峰型变得扁平.其他位置的光电导率特征峰在施加应变后一部分发生合并,一部分发生劈裂,如图6(a)的①区域中原本处于能量259 和302 meV 处的两个峰在施加应变后合并成了在270 meV 处的单峰(在1.47°时也有类似的现象),这些变化均源于能带在应变下的改变,而从态密度中可以更加直观地展现这些变化:峰的合并源于态密度在施加应变后峰型变得扁平,如图5(a)的①区域,该区域所代表的能带与第一非莫尔价带跃迁贡献出图6(a)中的特征峰值;峰的劈裂在1.47°时尤其突出,如图5(b)中的②区域,其态密度自下而上从未施加应变时的“尖锐边缘型”劈裂为施加应变后的“阶梯型”,而由此区域能带与莫尔价带之间的跃迁贡献了图6(b)中的光电导率特征峰,所以光电导率的特征峰也由“边缘型”转变为“阶梯型”.对于1.05°的TBG,原本未施加应变时零能处的单个范霍夫奇点在应变下产生劈裂(图5(a)),莫尔能带在应变下分离(图3),这种结果将导致电子的关联效应减弱,不利于实验对电子强关联性质的研究;而伴随着单个范霍夫奇点的劈裂,1.05°下零能处的Drude峰(图6(a)中黑色曲线的零点峰)也一并消失,此时光学吸收谱只包含带间跃迁的贡献,这一效应在之前扫描隧道显微镜的实验[14]中得到证实.

图5 TBG的态密度 (a) TBG 在转角为1.05°,应变大小分别为0%(黑色),0.3%(蓝色),0.6%(红色)时的态密度;(b) TBG 在转角为1.47°,应变大小分别为0%(黑色),0.3%(蓝色),0.6%(红色)时的态密度;蓝色虚线表示正文中所讨论的①,②区域的边界Fig.5.(a) Density of states(DOS) of TBG with 0%(black curve),0.3%(blue curve) and 0.6%(red curve) uniaxial heterostrain at 1.05°;(b) DOS of TBG with 0%(black curve),0.3%(blue curve) and 0.6%(red curve) uniaxial heterostrain at 1.47°.Blue dash lines in panel(a) and(b) represent the boundary of ①,② region,respectively.

图6 TBG的光电导 (a) TBG 在转角为1.05°,应变大小分别为0%(黑色),0.3%(蓝色),0.6%(红色)时的光电导;(b) TBG 在转角为1.47°,应变大小分别为0%(黑色),0.3%(蓝色),0.6%(红色)时的光电导;蓝色虚线表示正文中所讨论的①,②区域的边界Fig.6.(a) Optical conductivity of TBG with 0%(black curve),0.3%(blue curve) and 0.6%(red curve) uniaxial heterostrain at 1.05°;(b) optical conductivity of TBG with 0%(black curve),0.3%(blue curve) and 0.6%(red curve) uniaxial heterostrain at 1.47°.Blue dash lines in panel(a) and(b) represent the boundary of ①,② region,respectively.

由此可见,施加一个比较大的应变对TBG的能带和光电导率均能产生较大的影响,所以通过改变施加应变的大小便能够调控TBG的电子性质.在实验上,TBG 样品往往受到具有空间不均匀性的应变作用,而本文的结果表明测量TBG 局域的光电导率就可以估计其局域的应变大小,从而得到应变的空间分布.

4 总 结

本文用连续模型从理论上研究了转角为1.05°和1.47°的TBG 光电导率在施加不同大小的单轴应变下所产生的变化.可见在两个转角下电导率谱主要由600 meV 以下的几个特征峰表征,这些特征峰与态密度中的一些范霍夫奇点能量位置之差相对应.随着应变的逐渐施加,能带与电导率谱也发生几乎连续的变化,而这些变化是比较明显的.也就是说,在实验上若不考虑应变对于TBG的影响将会使理论结果与实际偏离.对于魔角TBG,应变使得平带带宽增大,这不利于研究强关联效应,所以实验在对于魔角TBG 系统中强关联态进行研究时应尽量减小应变;而对于非魔角处的石墨烯,应变的施加也可以使其形成平带[9,10].这表明通过原位改变应变大小来调控TBG 电子的强关联效应是可行的.本文的理论结果也为用局域光电导率确定TBG的空间应变大小分布提供了参考.