MoTe2-x Kagome单层中的Kagome能带与磁性PHYS. REV. B 113, 155405 (2026)MoTe2-x Kagome单层中的Kagome能带与磁性Kagome Bands and Magnetism in MoTe2-x Kagome Monolayers导读 导读Kagome晶格因其平带、Dirac锥和Van Hove奇点而成为凝聚态物理的圣杯体系。但体材料中Kagome能带常与体态杂化难以独立调控。本文提出创新方案利用MoTe2单层中的镜面孪晶界MTB环形结构构建13种Kagome超晶格通过调控MTB环的尺寸N和排列T实现Kagome能带的可编程设计。核心发现4种构型展现非零Z2拓扑不变量量子自旋Hall绝缘体候选2种构型出现Stoner铁磁性为2D Kagome材料的设计开辟了新范式。【一、前言背景】Kagome晶格从几何阻挫到拓扑量子态Kagome晶格六角星形晶格因其独特的几何结构而成为凝聚态物理的圣杯体系。它天然支持三种标志性电子特征平带flat band、Dirac锥和Van Hove奇点。这些特征为磁有序、超导、电荷密度波CDW等量子相提供了理想舞台。然而体材料Kagome体系面临一个根本挑战Kagome特征能带特别是平带经常与体态能带杂化导致Kagome电子特征被抑制或偏离费米面。Co3Sn2S2磁性Weyl半金属和AV3Sb5超导CDW是著名的体Kagome材料但它们的Kagome特征能带都受到capping层的干扰。二维Kagome材料是突破这一困境的自然方案。但传统的2D Kagome材料如Nb3Cl8难以从体材料中剥离非vdW层间相互作用而moire超晶格虽然能产生电子Kagome晶格却面临结构重构和稳定性问题。MTB-loop策略在MoTe2中构建Kagome超晶格核心创新利用MoTe2单层中的镜面孪晶界MTB环形结构MTB loops构建Kagome超晶格。MTB是MoTe2中常见的线缺陷通过移除部分Te原子形成。不同大小N1,3,6,10,15,21和排列方式T0,T1,T3,T6的MTB三角形环构成13种MoTe2-x单层。方法体系VASPPBEvdW-DF/optB86b, ENCUT500 eV, 5x5x1/3x3x1 k点- Wannier90Mo-4dTe-5p MLWFs, 含SOC- WannierToolsZ2拓扑不变量- STM实验4.3K, 分子束外延生长。关键发现(1) 在MoTe2本征带隙内识别出两组Kagome能带KBS1和KBS2分别来自MTB边缘和顶点的Mo-d轨道(2) 4种构型展现非零Z2拓扑不变量暗示拓扑绝缘体相(3) 2种构型为金属性并出现Stoner铁磁性。MoTe2-x Kagome超晶格研究流程。MTB-loop构建N尺寸排列方式T值- 13种构型 - VASP DFTPBEvdW-DF- Kagome能带识别KBS1KBS2- Wannier90SOCZ2拓扑不变量- 形成能分析。核心发现4种拓扑相Z2非零 2种磁性相Stoner铁磁。【二、研究方法】VASP DFT计算设置vdW修正与Kagome超胞计算设置PBE泛函 PAW赝势, ENCUT500 eV弛豫和电子结构, 力收敛0.01 eV/A。vdW修正使用vdW-DF级别的optB86b泛函这对于MoTe2层状材料至关重要。真空层15 A。k点N1/N3用5x5x1N6/N10/N15/N21用3x3x1。MTB-loop分类N值三角形数序列{1,3,6,10,15,21}表示三角环内Te2二聚体数量。T值{0,1,3,6}表示由三个MTB三角顶点包围的Te2-V二聚体数量反映排列方式。T不能超过N/2除N1T1外。空间群分析N1T1、N3T6、N6T3属于P-62mNo.189, D3h点群具有镜面sigmah对称性可能形成Kagome晶格。其余属于P-6No.174, C3h点群破缺面外镜面对称性出现呼吸畸变。Wannier90 WannierTools拓扑不变量计算Wannier90紧束缚模型基于Mo-4d和Te-5p轨道构建MLWFs包含SOC。SOC是计算Z2拓扑不变量的必要条件因为Z2拓扑由SOC打开的能隙保护。Z2不变量计算通过WannierTools计算Wilson loop的演化确定Z2拓扑分类。对于时间反演不变的2D绝缘体Z20为普通绝缘体Z21为量子自旋Hall绝缘体。Stoner铁磁性判据对于金属性构型计算D(EF)*I 1Stoner判据判断是否出现自旋极化基态。D(EF)为费米面态密度I为Stoner交换参数。Stoner铁磁性判据D(E_F)为费米面态密度I为Stoner交换参数。乘积1时非磁态不稳定体系自发磁化Z2拓扑不变量通过Wilson loop沿半个布里渊区的绝热演化计算模2后为拓扑不变量【三、核心结果】图 1MTB三角形环的尺寸和排列。(a) N1,3,6,10的MTB环原子模型。(b-d) N6T0、N6T1、N6T3排列方式的原子模型。MTB-loop的几何分类从原子缺陷到Kagome超晶格MTB-loop的核心几何参数N环内Te2二聚体数量决定环的大小T三环排列中顶点包围的Te2-V数量决定排列方式。N和T的序列都遵循三角形数{1,3,6,10,15,21,...}。超晶格尺寸从12.66 A到30.49 A覆盖了从近原子尺度到介观尺度的范围。不同的N和T组合产生不同的空间群P-62m vs P-6决定了Kagome晶格是否存在以及呼吸畸变的程度。N1T1Mo5Te8是CTColoring Triangular单层具有最高的对称性D3h点群是唯一具有完整镜面sigmah对称性的构型。它展现出两组Kagome能带KBS1和KBS2中间有~0.1 eV的带隙。图 2不同MTB-loop尺寸和排列的MoTe2-x单层能带结构。(a) N1T1的能带结构B1-B6标记六条特征能带标注Z2拓扑不变量。(b,c) B2/B3和B6在Gamma点的波函数模方。(d-g) N3,6,10,15的NxT1 NM能带。(h-j) N1,3,6,10的NxT0 NM能带。Kagome能带的实空间起源边缘Mo vs 顶点MoKBS1红色能带组主要来自MTB环边缘的Mo-dz2轨道与dxy/dx2-y2杂化。这些边缘Mo原子形成CT晶格Coloring Triangular lattice贡献了第一组Kagome能带。KBS2蓝色能带组主要来自MTB环顶点的Mo面外d轨道。这些顶点原子形成独立的Kagome子晶格贡献了第二组Kagome能带。两组Kagome能带在能量上分离实现双Kagome能带结构。这一发现的意义Kagome能带不是来自整个MoTe2晶格而是来自MTB缺陷结构。这意味着通过设计MTB环的尺寸和排列可以精确调控Kagome能带的能量位置、带宽和拓扑性质。图 3形成能和STM图像。(a) 10种单层从Te-deficient到Te-rich条件的形成能。(b) (a)中虚线框区域的放大。(c-f) N1T1、N6T1、N10T1、N15T1的STM形貌图偏压1.0V。形成能与化学势稳定性窗口不同MTB构型的稳定性取决于Te化学势。在Te-deficient条件下N较大的构型更多Te空位更稳定在Te-rich条件下N较小的构型更稳定。4种构型在特定化学势窗口内具有最优稳定性。STM实验验证N1T1、N6T1、N10T1、N15T1的STM形貌图与DFT模拟的原子结构高度一致确认了MTB-loop构型的可实现性。STM图像中的亮斑对应MTB环的顶点和边缘。实验生长条件Te/Mo通量比30衬底温度~200degC随后350-500degC退火1小时。通过调控退火温度可以获得不同有序度的MTB超结构。图 4N6T1和N10T1的Stoner铁磁性。(a,b) N6T1和N10T1的自旋极化能带结构展示自旋劈裂。(c) N6T1的自旋密度分布。(d) N10T1的自旋密度分布。(e) Stoner参数D(EF)*I。Stoner铁磁性Kagome平带驱动的磁性N6T1和N10T1是金属性构型出现Stoner铁磁性。物理根源Kagome平带产生高费米面态密度D(EF)满足Stoner判据D(EF)*I1体系自发磁化。自旋密度分布磁性主要局域在MTB环边缘的Mo原子上形成铁磁有序。这与Kagome晶格中平带的实空间局域化特征一致平带波函数在实空间高度局域导致强的在位库仑相互作用。与体材料Kagome磁体的比较体材料中Kagome磁性与巡游电子耦合更复杂而2D MTB-loop体系中磁性更局域、更可控。通过改变N和T可以调节D(EF)从而开关磁性。【DFT Tips】【DFT Tip 1】vdW修正的选择optB86b-vdW vs DFT-D3本文使用vdW-DF/optB86b这是vdW-DF家族中的交换泛函。对于MoTe2这类层状TMD材料vdW修正至关重要。optB86b在描述层状材料的层间距和结合能方面表现优异。对比DFT-D3Grimme是半经验色散修正但对金属体系可能不准确。vdW-DF是非局域密度泛函更物理但2-3倍。对于TMD材料推荐vdW-DF/optB86b或rev-vdW-DF2。常见陷阱vdW修正对晶格常数影响可达1-2%对应变计算和电子结构有显著影响。必须在所有计算中一致使用vdW修正包括弛豫、电子结构和声子计算。【DFT Tip 2】MTB缺陷建模的超胞大小选择MTB-loop超胞尺寸从12.66 A到30.49 A对应的k点采样从5x5x1到3x3x1。超胞越大k点可以越稀疏但需要确保总能收敛。缺陷超胞的周期性镜像相互作用MTB环之间的弹性相互作用通过周期性边界条件传播。如果超胞不够大相邻MTB环的应力场会重叠导致能量的系统误差。对于N21构型30.49 A的超胞尺寸应足够。验证方法对于最小构型N1T1用更大的超胞如2x2验证总能和能带结构是否收敛。如果对N1T1收敛对大构型也应收敛。【DFT Tip 3】Wannier90中SOC处理与Z2拓扑计算Z2拓扑不变量的计算必须在包含SOC的Wannier紧束缚模型上进行。VASP中SOC计算需要非共线磁结构LNONCOLLINEAR.TRUE. LSORBIT.TRUE.。Wannier90设置对于Mo-4d和Te-5p轨道构建MLWFs初始投影需覆盖所有自旋分量spinor投影。disentanglement窗口需包含费米面附近的所有Kagome能带。常见陷阱SOC下Wannier拟合的spread通常比无SOC时大因为自旋混合增加了波函数的复杂度。如果拟合质量差可以增加inner窗口的范围或增加初始投影数量。拟合质量通过对比DFT和Wannier插值的能带结构来评估。【DFT Tip 4】Stoner判据的DFT计算实践Stoner判据D(EF)*I1中D(EF)从非磁态DOS获得但Stoner参数I不能直接从DFT获取。常见做法通过比较NM和FM总能量差来反推有效I或使用文献值。注意DFT中的Stoner铁磁性是平均场图像忽略了自旋涨落。在2D体系中Mermin-Wagner定理禁止有限温度下的长程磁有序因此DFT预测的FM基态可能在实际温度下被自旋涨落破坏。建议对于DFT预测的2D铁磁体务必补充Monte Carlo模拟基于DFT交换耦合来估计磁转变温度评估磁有序在实际温度下是否可行。【DFT Tip 5】Kagome能带的识别平带、Dirac锥、Van Hove奇点在DFT能带中识别Kagome特征(1) 平带在较大k空间范围内色散极小的能带0.1 eV通常出现在费米面附近(2) Dirac锥能带在K点线性交叉无带隙(3) Van Hove奇点能带在马鞍点处的态密度对数发散。Kagome能带的折叠分析MTB超晶格导致能带折叠折叠后的能带可能来自原始MoTe2能带的不同区域。需要通过轨道投影和波函数可视化来区分真正的Kagome能带和折叠伪影。常见陷阱并非所有平带都是Kagome平带。缺陷态、表面态、悬挂键态都可以产生平带。Kagome平带的特征是波函数在实空间局域在Kagome晶格的六边形上能带在动量空间具有Kagome特有的紧束缚色散关系。【DFT Tip 6】形成能计算中的化学势参考形成能E_form E(MoTe2-x) - E(MoTe2) x*mu_Te其中mu_Te是Te的化学势。化学势的参考态选择至关重要(1) Te-rich极限mu_Te E(Te_bulk)/2Te晶体(2) Te-poor极限mu_Te E(MoTe2) - E(Mo_bulk) - 2*mu_Te由MoTe2的稳定性决定。在不同化学势下不同构型的相对稳定性发生变化。本文的相图覆盖了从Te-deficient到Te-rich的整个范围为实验上控制Te化学势来选择特定MTB构型提供了指导。常见错误使用元素基态能量作为化学势如Te二聚体在气相中的能量这对应极端Te-rich条件在固体生长中不现实。【DFT Tip 7】2D Kagome材料的SOC效应拓扑能隙打开在Kagome晶格中SOC可以打开Dirac锥处的能隙产生拓扑非平庸的能带结构。SOC导致的能隙大小取决于组成元素的原子序数TeZ52的SOC较强MoZ42中等。对于MoTe2-xTe的SOC约为0.5-0.8 eV原子能级但在固体中由于晶体场效应会减小。SOC打开的能隙通常在10-100 meV量级取决于具体构型。SOC计算中的k点收敛SOC能隙对k点采样敏感建议使用比无SOC时更密的k点网格如将3x3x1提升到5x5x1。【DFT Tip 8】STM图像模拟Tersoff-Hamann近似本文比较了实验STM图像和DFT模拟的原子结构。STM图像模拟通常使用Tersoff-Hamann近似隧穿电流正比于针尖费米能级处样品的局域态密度积分。VASP中进行STM模拟的步骤(1) 自洽计算得到CHGCAR(2) 使用p4vasp或VASPKIT后处理在特定偏压窗口内积分态密度得到STM图像。注意偏压窗口的选择对应实验的偏压值本文1.0V。常见陷阱STM模拟对针尖高度和偏压极其敏感。不同的针尖-样品距离会产生截然不同的图像。建议对多个高度进行模拟选择与实验最匹配的结果。【知识扩展】【知识扩展 1】Kagome晶格的三重电子特征平带、Dirac、vHs【理论解释】Kagome晶格是二维晶格中具有最小带宽平带的晶格之一。三个特征来自紧束缚模型的解析解(1) 平带色散完全平坦群速度为零电子高度局域(2) Dirac锥能带在K点线性交叉电子表现为无质量Dirac费米子(3) Van Hove奇点态密度在M点呈现对数发散。【平带的物理后果】平带导致态密度发散增强电子关联效应。当平带在费米面附近时可以驱动磁有序、超导、Wigner晶化、分数量子Hall效应等强关联量子相。【经典参考】Mielke, JPA 24, 3311 (1991) -- Kagome平带紧束缚模型Yin et al., Nature 583, 533 (2020) -- CoSn KagomeOrtiz et al., PRL 125, 247002 (2020) -- AV3Sb5超导。【迁移能力】Kagome晶格的概念可推广至Lieb晶格、pyrochlore晶格、hyperkagome晶格等几何阻挫体系。【知识扩展 2】MTB镜面孪晶界在TMD中的物理【理论解释】MTB是TMD材料中最常见的线缺陷之一。在H相MoTe2中MTB将Te原子从三角棱柱配位转变为八面体配位局部形成1T相。MTB沿镜面方向延伸因此得名。【MTB的电子特征】MTB在TMD带隙中引入一维金属态形成量子线。这些MTB态与体态在能量上分离几乎不受衬底影响。正是这一特性使得MTB成为构建Kagome超晶格的理想构件。【经典参考】Lin et al., Nat. Nanotechnol. 9, 436 (2014) -- MoSe2 MTB实验Ma et al., ACS Nano 11, 513 (2017) -- MoTe2 MTBBarja et al., Nat. Phys. 12, 751 (2016) -- MTB金属态。【迁移能力】MTB-loop策略可推广至其他TMD材料MoS2, MoSe2, WS2, WSe2, WTe2为构建2D Kagome体系提供通用方案。【科研经验】【科研经验 1】从DFT能带中识别Kagome特征的挑战问题DFT能带结构复杂如何区分真正的Kagome特征和偶然的能带折叠/杂化效应原因MTB超晶格导致能带严重折叠折叠后的能带中可能出现伪平带或伪Dirac锥。仅凭能带色散无法可靠识别Kagome特征。解决方案(1) 做轨道投影能带fat band确认Kagome能带的轨道成分来自MTB边缘/顶点(2) 计算波函数模方partial charge density在实空间确认电子密度局域在Kagome晶格上(3) 构建紧束缚模型验证Kagome能带的跃迁参数是否与Kagome晶格的紧束缚解一致。建议本文的波函数可视化图2b,c是识别Kagome能带的黄金标准。在论文中务必展示实空间波函数分布而不仅是能带图。【科研经验 2】2D缺陷超晶格的热力学稳定性评估问题DFT预测的MTB-loop构型是否真的能在实验中稳定存在仅凭总能比较是否足够原因MTB-loop的形成涉及大量Te空位在生长过程中Te原子的扩散动力学决定了是否能形成有序的MTB超晶格。总能最低的构型可能在动力学上不可及。解决方案(1) 计算形成能作为Te化学势的函数确定稳定化学势窗口(2) 计算MTB-loop形成的能垒NEB方法评估动力学可行性(3) 做AIMD模拟验证构型在生长温度下的热稳定性。建议本文的STM实验验证是理论预测的最强有力支持。纯理论预测的2D缺陷超晶格如果没有实验验证需要格外谨慎。【如果是我我还会继续算】【继续算 1】HSE06杂化泛函验证Kagome能带为什么值得算PBE对平带和能隙的描述可能不准确HSE06可以给出更可靠的Kagome能带结构特别是确认Z2拓扑的能隙是否在HSE06中仍然存在。能回答的问题HSE06下的Z2拓扑分类是否与PBE一致Kagome平带的能量位置是否改变适合体系所有半导体Kagome材料。输入HSE06计算减少k点【继续算 2】Berry曲率与反常Hall效应为什么值得算对于磁性Kagome构型N6T1, N10T1计算Berry曲率可以揭示反常Hall效应AHE和拓扑Hall效应。Kagome晶格是AHE的天然平台。能回答的问题磁性构型的AHE有多大Berry曲率在动量空间如何分布是否存在量子反常Hall效应适合体系所有磁性Kagome材料。输入Wannier90 WannierToolsBerry曲率AHE。【继续算 3】EPC与超导可能性为什么值得算Kagome平带产生高态密度如果同时存在较强的EPC可能驱动超导。AV3Sb5的超导性就与Kagome能带相关。能回答的问题MoTe2-x Kagome构型是否存在电子-声子耦合驱动的超导EPC强度和Tc估计是多少适合体系所有Kagome金属。输入QE PHonon EPW。【继续算 4】边缘态与量子自旋Hall效应为什么值得算对于Z21的构型计算纳米带的边缘态可以直接验证量子自旋Hall效应。边缘态是拓扑非平庸的最直接证据。能回答的问题Z2非零构型的边缘态是否无带隙边缘态的穿透深度和自旋纹理如何适合体系所有2D拓扑绝缘体。输入构建纳米带模型DFT计算或Wannier紧束缚模型到高。【继续算 5】外场调控电场和应力对Kagome能带的影响为什么值得算垂直电场可以打破面外对称性调控Kagome能带的SOC能隙。双轴应力可以改变跃迁参数调控平带宽度和Dirac锥位置。能回答的问题电场能否实现拓扑相变应力能否将半导体的Kagome构型转变为金属或磁性适合体系所有2D材料。输入DFT电场EFIELD标签或应力计算到中等。【继续算 6】光学性质Kagome能带的激子效应为什么值得算N1T1的CT构型中VBM和CBM分别位于两个平带形成~0.1 eV的带隙。这种能带排列可能支持激子绝缘体相。计算光吸收谱和激子结合能可以验证。能回答的问题CT构型是否支持激子绝缘体激子结合能多大光吸收谱的特征峰对应哪些Kagome能带间的跃迁适合体系所有半导体Kagome材料。输入DFTBSEGWBSE或使用简化的模型激子计算。Dai, Zhang, Pan, Wang, Zhang, Cheng, Ji | Phys. Rev. B 113, 155405 (2026) | Kagome MTB MoTe2 拓扑