插层TMD中的交错磁性、Kagome平带与Weyl费米子

插层TMD中的交错磁性、Kagome平带与Weyl费米子

NPJ QUANTUM MATER. 11, 63 (2026)

插层TMD中的交错磁性、Kagome平带与Weyl费米子态

Altermagnetism, Kagome Flat Band, and Weyl Fermion States in Intercalated TMDs

导读 导读：交错磁体（Altermagnet）兼具铁磁自旋劈裂和反铁磁零净磁矩，是下一代自旋电子学的关键材料。本文在插层TMD家族XY4Z8（X=Mn,Fe,Co,Ni,Cr,V; Y=Nb,Ta; Z=Se,S）中系统搜索AM候选材料，发现FeNb4S8、CoNb4Se8、CoTa4Se8、FeNb4Se8为AM化合物。u/v比（层间/层内间距）作为结构描述符可预测磁基态。更引人注目的是，这些材料同时实现了AM、Kagome平带和Weyl费米子的三重量子态共存。

一、前言背景

交错磁性：从概念到材料实现

交错磁体（Altermagnet, AM）是2021年以来凝聚态物理最重要的发现之一。它突破了传统铁磁/反铁磁的二分法：实空间中自旋反平行（净磁矩为零），但倒空间中能带出现非相对论自旋劈裂。这种独特的d波/g波/i波形式的自旋劈裂源于自旋空间群（Spin Space Group）的对称性。

2022年Smejkal、Sinova、Jungwirth在PRX正式提出AM的分类框架后，全球掀起了AM材料搜索热潮。至今已发现MnTe、RuO2、CrSb、Mn5Si3、V2Te2O等数百种候选材料。2024年Nature连续发表MnTe和RuO2的实验验证，标志着AM研究进入实验时代。

本文在插层TMD（过渡金属二硫族化物）家族XY4Z8中系统搜索AM候选材料，发现了一个全新的AM材料设计平台，同时揭示了AM与kagome平带、Weyl费米子的三重共存。

三重量子态共存：AM + Kagome平带 + Weyl费米子

核心创新：在同一材料家族中同时实现交错磁性、kagome平带和Weyl费米子态。这三种量子态通常需要在不同材料中分别实现，XY4Z8将它们统一在一个平台中。

结构-磁性关联：u/v（层间/层内间距比）作为结构描述符，u/v<0.91对应A型AFM（AM候选），u/v>=0.91对应FM。这一简单的几何判据为高通量筛选AM材料提供了实用工具。

三个物理层次：(1) 非相对论层次：AM自旋劈裂（g波对称性）；(2) SOC层次：AM节点线转化为Weyl点+表面费米弧；(3) 轨道层次：插层诱导有效kagome晶格->平带+Dirac锥共存。

XY4Z8交错磁性研究流程。VASP PBE+U（晶体结构、u/v比分析、磁基态扫描）-> AM自旋劈裂分析（节点面/非节点面、Fermi面自旋纹理、自旋空间群对称性）-> Weyl费米子态（SOC转化节点线为Weyl点+表面费米弧）-> Kagome平带（插层诱导有效kagome晶格）。

二、研究方法

VASP PBE+U：高通量磁基态扫描

计算设置：PBE泛函+PAW赝势, ENCUT=450 eV, 9x9x4 Monkhorst-Pack k点。2x2x1超胞包含两个插层原子。弛豫时同时采样FM和AFM自旋排列以确定能量最低的磁基态。

Hubbard U校正：对插层过渡金属d轨道使用Dudarev的GGA+U方法，Ueff=0-3 eV系统扫描。结果表明磁基态对U值不敏感，AM候选材料的定性结论在U值范围内稳健。

u/v比作为结构描述符：u为层间间距（相邻vdW间隙中插层原子的垂直距离），v为层内间距（同一vdW间隙中插层原子的面内距离）。u/v<0.91->A型AFM, u/v>=0.91->FM。NiNb4Se8是唯一例外（u/v<0.91但FM），因为Ni的近满d壳层抑制了超交换。

u/v比与磁基态的经验关系：u/v<0.91对应A型AFM（AM候选），u/v>=0.91对应FM。来源于X-Nb-X键角控制的超交换与RKKY竞争

自旋空间群（SSG）理论：AM的对称性根源

SSG框架：在无SOC条件下，对称性操作为[R_s||R_r]，其中R_r为实空间对称性，R_sinSO(3)为独立自旋旋转。对于A型AFM结构，SSG为P-1_63/-1m1m-1cinfinitym1。

关键对称性：[C2||M_z]和[C2||C_6zt]连接相反自旋子晶格，保持磁构型不变。在k_z=0节点面上，[C2||M_z]留下k_z=0平面不变，产生Kramers-like二重简并，保护自旋简并。在k_z=pi/2c非节点面上，M_z不再保持面内动量不变，允许自旋劈裂。

SOC极限：自旋旋转不再独立于晶格旋转，SSG退化为磁空间群（MSG）P6'_3/m'm'c。

自旋空间群操作：[R_s||R_r]，R_s为独立自旋旋转，R_r为实空间对称性。AM的g波自旋劈裂来源于SSG的特定对称性约束

紧束缚模型：Kagome平带的起源

从六方到Kagome：插层原子位于部分Nb位正上方，将其On-site势能提高（epsilon_Nb1 >> epsilon_Nb2,3,4），有效"淬灭"中心Nb位。剩余三个Nb位形成有效kagome晶格。

TB哈密顿量：4x4矩阵描述四个Nb位之间的跃迁，t为kagome位之间的跃迁，t'为中心位与kagome位之间的跃迁。当t'->0或epsilon_Nb1->infinity时，模型退化为理想kagome晶格，出现完全平带。

DFT验证：轨道投影能带结构显示，中心Nb位的d_xy轨道对Dirac锥和平带的贡献被Co插层淬灭，kagome Nb位的d_xy轨道贡献主导。TB模型与DFT结果定性一致。

紧束缚哈密顿量：描述四个Nb位之间的跃迁，t为kagome位间跃迁，t0为中心位与kagome位间跃迁。当t0->0时退化为理想kagome模型

三、核心结果

图 1：XY4Z8的晶体结构和磁构型。(a) 侧视图，展示插层原子在vdW间隙中的位置，u和v定义层间/层内间距。(b) 俯视图，展示插层原子在面内的排列。(c) AFM/AM和FM磁构型示意图。(d) 六方晶格第一Brillouin区和高对称路径。

u/v比：从结构到磁性的预测性描述符

u/v<0.91体系：FeNb4S8(0.877), CoNb4Se8(0.893), CoTa4Se8(0.906), FeNb4Se8(0.909)。X-Nb-X键角较小，X-Nb-X-Se-Nb-X超交换路径增强，层间AFM耦合占优，形成A型AFM->AM候选。

u/v>=0.91体系：CrNb4Se8(0.911), VNb4Se8(0.919), MnNb4Se8(0.928), MnTa4S8(0.935)。层间间距增大，垂直交换减弱，面内RKKY占优，FM基态。

NiNb4Se8例外：u/v=0.897<0.91但FM基态。Ni的d轨道近乎填满，Pauli阻塞抑制了超交换路径，RKKY型交换占主导。FM和AFM能量差极小，表明Ni化合物处于磁相边界。

图 2：AM候选材料的自旋分辨能带结构（含SOC）。(a) CoNb4Se8, (b) FeNb4Se8, (c) CoTa4Se8, (d) FeNb4S8。左：k_z=0节点面（Gamma-M-K-Gamma），自旋简并。右：k_z=pi/2c非节点面（Gamma'-M'-K'-Gamma'），明显的自旋劈裂。

g波对称性自旋劈裂：AM的电子指纹

节点面k_z=0：镜面对称性[C2||M_z]保护整个平面的自旋简并，Gamma-M-K-Gamma路径上spin-up和spin-down能带完全重合。这是AM区别于FM的关键特征--FM在所有k点都有自旋劈裂。

非节点面k_z=pi/2c：M_z不再保持面内动量不变，自旋简并被解除。沿Gamma'-M'路径出现显著自旋劈裂，且劈裂模式呈现g波对称性（4个节点）。

CoTa4Se8的自旋劈裂最大：在费米能级处约110 meV，在-0.25 eV处约154 meV。其余AM候选材料（FeNb4Se8, CoNb4Se8, FeNb4S8）的劈裂幅度较小但定性相似。

图 3：AM材料的Fermi面自旋纹理。(a,b) CoNb4Se8和FeNb4Se8在无SOC时的自旋分辨Fermi面。(c,d) 对应含SOC的Fermi面。k_z=0面完全自旋简并，k_z=pi/2c面明显自旋极化。

Fermi面自旋纹理：AM vs FM的本质区别

AM的Fermi面：在k_z=0面完全自旋简并（红蓝重叠），在k_z=pi/2c面出现强自旋极化。SOC下虽然s_z不再是好量子数，但残余的自旋纹理与无SOC时高度相似。

FM的Fermi面：Zeeman-like的全Brillouin区均匀自旋劈裂，无节点面。所有k点自旋简并均被解除。这是AM区别于FM的Fermi面指纹。

实验验证展望：ARPES可以直接测量Fermi面的自旋纹理。k_z=0面的自旋简并和k_z=pi/2c面的自旋劈裂构成AM的明确ARPES证据。

图 4：CoNb4Se8中的Weyl点和拓扑表面态。(a) 体Brillouin区中Weyl点分布及(001)面投影。(b) 表面谱和Weyl锥色散。(c) 表面Fermi弧（含双费米弧）。(d) 杂化Wannier电荷中心演化，确认chi=+/-1拓扑电荷。

SOC驱动AM节点线->Weyl点的拓扑转化

无SOC时：AM的自旋劈裂导致相反自旋的能带在动量空间中形成延伸的节点线（nodal lines）。这些节点线由自旋通道的退耦合保护。

SOC引入后：SOC耦合自旋通道，打开大部分节点线简并。但反幺正对称性操作（空间对称性+时间反演）保护了某些孤立动量点的简并。连续节点线坍缩为离散的Weyl点。

CoNb4Se8中：12对Weyl点（chi=+1和chi=-1）分布在费米能级+/-0.2 eV范围内。表面态显示双费米弧（两个相同手性Weyl点的投影重叠），确认了体-边界对应。

图 5：CoNb4Se8中的Kagome平带和晶格结构。(a) 原子投影能带。(b,c) 中心Nb位和Kagome Nb位的轨道分辨能带。(d) 六方->Kagome晶格映射。(e) 紧束缚模型谱：六方->Kagome-like->理想Kagome的演化。

Kagome平带：插层诱导的轨道选择性电子结构

中心Nb位淬灭：位于Co插层正下方的Nb原子（Site 1）与Co d_z2轨道强烈杂化，其d_xy轨道对Dirac锥和平带（顶）的谱权重被显著抑制。

Kagome Nb位主导：其余三个Nb原子（Sites 2,3,4）的d_xy轨道贡献了Dirac锥和平带（顶）的大部分谱权重。部分平带（底）则由中心Nb位贡献。

TB模型验证：当中心Nb位On-site势能以epsilon_Nb1/t~-2.2抑制且t'/t~-0.4时，TB模型精确重现了DFT的Kagome-like能带特征。在理想Kagome极限（t'=0），平带变得完全平坦。

DFT Tips

【DFT Tip 1】高通量磁基态扫描：2x2x1超胞策略

对于插层TMD体系XY4Z8，使用2x2x1超胞包含两个插层原子，可以同时采样FM和AFM（A型、C型、G型）磁构型。超胞大小必须足够大以容纳至少两个磁性子晶格的原子。

弛豫策略：(1) 先固定晶格常数弛豫原子位置（ISIF=2）；(2) 再比较不同磁构型的总能量。不要先确定磁态再弛豫--磁态和结构是耦合的，不同磁态可能弛豫到不同的平衡结构。

常见错误：在1x1x1原胞中计算磁基态（原胞中只有一个插层原子），无法描述层间AFM耦合。必须使用包含至少两个插层原子的超胞。

【DFT Tip 2】AM自旋劈裂的DFT验证：节点面与非节点面

验证AM的关键步骤：(1) 计算k_z=0和k_z=pi/2c两个平面的能带结构；(2) k_z=0面应自旋简并，k_z=pi/2c面应出现自旋劈裂；(3) 如果两个面都自旋简并->常规AFM；如果两个面都自旋劈裂->FM。

能带计算路径：k_z=0面使用Gamma-M-K-Gamma；k_z=pi/2c面使用Gamma'-M'-K'-Gamma'（Gamma'=(0,0,pi/2c), M'=(1/2,0,pi/2c), K'=(1/3,1/3,pi/2c)）。需要在KPOINTS中手动指定这些k点。

注意：含SOC计算中s_z不再是好量子数，但VASP输出的投影仍可用于定性判断自旋纹理。真正严格的AM验证需要自旋空间群分析。

【DFT Tip 3】GGA+U中Ueff的系统扫描策略

本文对插层过渡金属使用Ueff=0-3 eV系统扫描，验证了磁基态对U值不敏感。这是高通量AM搜索的最佳实践--先确定U值范围内的定性结论，再选择代表性U值进行详细分析。

注意：对于4d/5d过渡金属（如本文的Nb、Ta），U值通常小于3d金属。Nb的Ueff通常为0-2 eV，Ta的Ueff为0-1 eV。对于插层3d金属（Mn, Fe, Co, Ni, Cr, V），Ueff=2-4 eV是合理范围。

常见陷阱：只使用单一U值计算磁基态，忽略了U值可能的敏感性。特别是对于接近磁相边界的化合物（如NiNb4Se8），U值可能改变磁基态。

【DFT Tip 4】Weyl点搜索：从节点线到SOC的拓扑转化

Weyl点搜索策略：(1) 先计算无SOC的能带结构，识别节点线（自旋简并的能带交叉）；(2) 开启SOC后，追踪节点线坍缩后的残余能带交叉点；(3) 用Wilson loop或杂化Wannier电荷中心计算手性电荷chi。

关键参数：Weyl点搜索需要极密的k点网格（至少100x100x100），最好使用Wannier90紧束缚模型插值。直接在DFT k点网格上搜索Weyl点几乎不可能准确定位。

常见错误：将SOC打开后的普通能带交叉误认为Weyl点。真正的Weyl点必须：(1) 位于一般k点（非高对称点）；(2) 线性色散；(3) 非零手性电荷chi=+/-1。

【DFT Tip 5】Kagome平带的DFT识别：轨道投影与位置分辨

Kagome平带的特征：(1) 在K点出现Dirac锥；(2) Dirac锥上方或下方出现弱色散的平带；(3) 平带由特定子晶格（kagome位）的轨道主导。

位置分辨PDOS：在POSCAR中区分中心Nb位和kagome Nb位，在INCAR中设置LORBIT=11输出原子投影DOS。通过比较不同位置Nb原子的PDOS，可以确认平带的kagome起源。

TB模型验证：构建简化的紧束缚模型，通过调整On-site势能和跃迁参数，重现DFT的能带特征。TB模型提供了对平带起源的物理直觉，是DFT计算的重要补充。

【DFT Tip 6】vdW层状材料的DFT弛豫：vdW修正的选择

插层TMD是典型的vdW层状材料，层间由弱vdW力结合。PBE泛函不包含vdW相互作用，可能低估层间距。建议使用DFT-D3或optB88-vdW等vdW修正。

本文未明确提及vdW修正，但PBE优化的结构通常能给出合理的u/v比趋势。对于需要精确层间距的研究（如压力效应、层间电荷转移），vdW修正是必须的。

注意：vdW修正可能影响层间距u，进而影响u/v比和磁基态预测。如果使用不同的vdW修正方案，需要重新验证u/v比与磁基态的关联。

【DFT Tip 7】自旋空间群分析工具：FindSpinGroup

FindSpinGroup（findspingroup.com）是在线自旋空间群识别工具，输入晶体结构和磁构型后自动输出SSG符号。对于AM候选材料的验证，SSG分析是理论上的"金标准"。

使用步骤：(1) 准备POSCAR格式的晶体结构（含磁矩信息）；(2) 上传到FindSpinGroup；(3) 检查输出的SSG是否允许动量依赖的自旋劈裂。

常见陷阱：SOC计算中SSG退化为MSG，自旋劈裂可能被SOC增强或改变。AM的"非相对论"本质应在无SOC计算中验证。

【DFT Tip 8】插层TMD的Fermi面计算：k_z依赖的截面

Fermi面计算需要3D k点网格（至少21x21x21），使用VASPKIT或FermiSurfer可视化。对于AM体系，必须分别绘制k_z=0和k_z=pi/2c两个截面的Fermi面。

自旋分辨Fermi面：在无SOC计算中，分开绘制spin-up和spin-down的Fermi面。在k_z=0面两者应完全重合，在k_z=pi/2c面应明显分离。这是AM最直接的Fermi面指纹。

注意：SOC下自旋不是好量子数，Fermi面不再有spin-up/spin-down之分。但VASP可以输出投影，通过颜色编码显示残余的自旋纹理。

知识扩展

【知识扩展 1】自旋空间群（SSG）与AM分类

【理论解释】传统磁空间群（MSG）将自旋旋转与晶格旋转绑定，不足以描述AM的非相对论自旋劈裂。自旋空间群（SSG）允许自旋空间和实空间的独立操作，是AM的对称性理论基础。

【方法比较】MSG->SSG->自旋点群->自旋Laue群，逐级递进。FindSpinGroup和Bilbao Crystallographic Server提供了AM分类的在线工具。

【经典参考】Litvin, Acta Cryst. A33, 279 (1977)--SSG奠基工作；Smejkal et al., PRX 12, 031042 (2022)--AM的SSG分类；Chen et al., PRX 14, 031038 (2024)--SSG的枚举与表示理论；Cheong & Huang, npj Quantum Mater. 10, 38 (2025)--AM分类综述。

【迁移能力】SSG框架适用于所有AM候选材料。对于任何声称的AM材料，SSG分析是理论验证的必要步骤。

【知识扩展 2】Kagome晶格中的平带物理

【理论解释】Kagome晶格由共享顶点的三角形组成，具有几何阻挫特征。在紧束缚模型中，Kagome晶格天然支持平带、Dirac锥和鞍点（van Hove奇点）的共存。平带来源于晶格中波函数的破坏性干涉导致的局域化。

【发现历程】Kagome平带最早在Mielke (1991)和Tasaki (1992)的紧束缚模型中被发现。近年来在CoSn、FeSn、AV3Sb5、YMn6Sn6等真实材料中实验观测到Kagome平带。平带的高态密度为关联电子态（磁性、超导、电荷序）提供了理想平台。

【经典参考】Mielke, J. Phys. A 24, L73 (1991)；Kang et al., Nat. Commun. 11, 4004 (2020)--CoSn中的Kagome平带；Li et al., Nat. Commun. 12, 3129 (2021)--YMn6Sn6中的Kagome物理。

【迁移能力】本文的"插层->有效Kagome"策略可推广到其他层状材料体系，通过插层或掺杂工程在原本非Kagome的晶格中诱导Kagome物理。

Sah, Lim, Conner, Chakraborty, Vignale, Chang, Sukhachov, Bian | npj Quantum Mater. 11, 63 (2026) | 交错磁性 Kagome平带 Weyl费米子 插层TMD