ADV. SCI. 反铁磁体MnTe2中的Plaid-like自旋劈裂与手性磁振子

ADV. SCI. 反铁磁体MnTe2中的Plaid-like自旋劈裂与手性磁振子

ADV. SCI. (2026) e76555

反铁磁体MnTe2中的Plaid-like自旋劈裂与手性磁振子

Plaid-Like Spin Splitting and Chirality of Magnon Bands in Antiferromagnetic MnTe2

导读 导读：交错磁体（Altermagnet）是2022年提出的磁性材料第三极，兼具反铁磁体的补偿磁序和铁磁体的自旋劈裂。本文通过手性分辨Raman光谱+DFT+U+线性自旋波理论，首次在非共面反铁磁体MnTe2中观测到磁振子的手性劈裂--RL和LR圆偏振通道的Raman强度出现Stokes/anti-Stokes反转，角分布图案旋转90度。DFT+自旋波理论揭示了磁振子能带中的plaid-like（格子花纹）自旋劈裂纹理，与电子能带的自旋纹理完全一致，确证了MnTe2的非共线交错磁性。

一、前言背景

交错磁体（Altermagnet）：磁性材料的第三极

2022年，Smejkal、Sinova、Jungwirth在PRX上正式提出altermagnetism概念，打破了铁磁体和反铁磁体的传统二分法。交错磁体兼具反铁磁体的补偿磁序（零净磁矩）和铁磁体的自旋劈裂能带，其物理根源是非相对论性的晶体对称性（旋转/镜面操作连接相反自旋子晶格），而非SOC。

交错磁体的关键特征：(1) 动量空间中的d波、g波或i波自旋劈裂；(2) 非相对论性（不含SOC）；(3) 能带和磁振子均可能呈现自旋劈裂。这使得交错磁体在THz级自旋动力学、无杂散场器件和手性磁振子学方面具有独特优势。

MnTe2是交错磁体的重要候选材料：黄铁矿型结构（Pa-3），Mn2+(S=5/2)形成FCC网络，四个Mn子晶格的自旋沿不同111方向排列，构成非共面补偿反铁磁序。本文首次通过手性分辨Raman光谱揭示了MnTe2中磁振子的手性劈裂。

手性磁振子：交错磁体的"指纹"激发

磁振子（magnon）是磁性有序态中的自旋波量子激发。在传统反铁磁体中，左手和右手磁振子严格简并。但在交错磁体中，复合自旋-晶格对称性操作可以解除这种简并，导致omega(k) != omega(-k)--即非互易磁振子色散。

本文的核心实验方法：手性分辨Raman光谱（RL和LR圆偏振配置），直接探测磁振子的手性不对称性。关键发现：2.5 meV低能磁振子在Stokes和anti-Stokes侧出现RL/LR强度反转，这是手性磁振子的"指纹"信号。

理论方法链：VASP DFT+U（基态电子结构）-> OpenMX（全带Hamiltonian）-> Jx代码（磁力线性响应，提取交换参数）-> Holstein-Primakoff变换（自旋波理论）-> Colpa对角化（磁振子能带+自旋纹理）。

MnTe2交错磁体手性磁振子研究流程。VASP DFT+U（基态电子结构） -> OpenMX DFT+U（全带Hamiltonian）-> Jx代码（磁力线性响应理论，提取J/D/Gamma交换参数）-> Holstein-Primakoff变换 + Colpa对角化（线性自旋波理论，4支磁振子能带）-> plaid-like自旋劈裂纹理 + 手性磁振子omega(k)!=omega(-k) -> 手性分辨Raman光谱实验验证（RL/LR偏振，Stokes/anti-Stokes反转）。

二、研究方法

VASP DFT+U：非共线磁基态的第一性原理描述

计算设置：VASP PAW赝势, PBE泛函, Ueff=4.2 eV（Dudarev方案, Mn-3d电子）, ENCUT=500 eV, 10x10x10 Gamma-centered k点。实验晶格常数6.90 A固定，原子位置弛豫至力<1 meV/A。

MnTe2的磁结构：4个Mn子晶格（Wyckoff 4a位置），自旋分别沿111, 1-1-1, -11-1, -1-11方向。这是非共面（non-coplanar）反铁磁序，总磁矩为零。

U值选择：Ueff=4.2 eV是基于Mn2+（d5高自旋）的典型值。作者测试了Ueff=1.0-5.0 eV范围，发现|D1|/J1比值变化<20.3%，证明交换参数对U值不敏感。同时也用cFLL+U+J（U=5.0 eV, J=0.8 eV）进行了交叉验证。

磁力线性响应理论：从DFT到交换参数

Jx代码基于Liechtenstein磁力理论（1987），通过计算自旋旋转引起的能量变化来提取交换参数。与传统的能量映射法（构建不同磁构型超胞，计算总能量差）不同，磁力线性响应直接在DFT基态上计算，不需要构建多个超胞。

核心公式：J_ij = (1/4pi) * Im(TrDelta_i \* G_ij\^up \* Delta_j \* G_ji\^down)，其中Delta_i为在位交换劈裂，G_ij为Green函数。这种方法可以同时提取各向同性交换J、DM矢量D和对称各向异性交换Gamma张量。

本文的交换参数：J1=0.675 meV（最近邻AFM），J2=-0.029 meV（次近邻FM），|D1|=0.063 meV（DM矢量，稳定非共线序）。D2=0（次近邻存在反演中心，DM被禁止）。

有效自旋哈密顿量：包含各向同性交换J、DM矢量D和对称各向异性交换Gamma张量。最近邻(ij)和次近邻(ij)分别考虑。

线性自旋波理论：Holstein-Primakoff变换与Colpa对角化

Holstein-Primakoff（HP）变换将自旋算符映射为玻色子算符：S_i^+ = sqrt(2S)*a_i, S_i^- = sqrt(2S)*a_i^dagger, S_i^z = S - a_i^dagger*a_i。对于S=5/2的Mn2+，1/S量子修正约为20%，线性自旋波理论是合理的近似。

对于非共线磁序，需要在每个子晶格的局域自旋坐标系中定义HP变换，然后通过Colpa方法对角化玻色子Hamiltonian。Colpa方法（1978）是处理非对角玻色子Bogoliubov-de Gennes方程的标准方法。

结果：4支磁振子能带（对应4个Mn子晶格），Gamma点能量E1=E2=E3=3.54 meV（上三支简并），E4=2.40 meV（最低支），与实验观测的2.5 meV和3.6 meV一致。

Holstein-Primakoff变换：将自旋算符映射为玻色子产生/湮灭算符

手性磁振子判据：交错磁体中磁振子色散不满足omega(k)=omega(-k)

三、核心结果

图 1：MnTe2的晶体结构和Raman光谱。(a) 黄铁矿型晶体结构，绿色和橙色箭头示意RL和LR散射配置。(b) T=10K下RL(红)和LR(蓝)偏振的代表性光谱，插图为光学声子对称性指认。(c) RL和LR的强度差（RCD信号）。(d,e) RL和LR偏振下温度依赖Raman光谱的彩色等高线图。(f) 10K-305K范围内RCD信号的温度演化。

Raman圆二色性（RCD）：手性磁振子的直接探针

RCD = I_RL - I_LR，正比于磁振子的手性不对称性。本文的核心实验发现：2.5 meV低能磁振子在Stokes和anti-Stokes两侧呈现相反的RCD信号，这是手性磁振子的独有特征。

物理机制：RL和LR偏振分别耦合到相反手性的磁振子。在交错磁体中，复合自旋-晶格对称性使左手和右手磁振子不再简并，导致RL和LR通道的散射强度不同。Stokes/anti-Stokes反转来源于磁振子产生和湮灭过程的手性选择规则。

对比：3.6 meV高能磁振子和声子模式均未显示RCD，说明手性不对称性仅存在于特定磁振子支。这与线性自旋波理论预测的表面不一致（高能支预测有更大自旋极化），但Raman强度还受散射截面和Raman顶点结构影响。

图 2：偏振分辨Raman光谱的角依赖性和对称性分析。(a-d) 10K和305K下平行和交叉偏振的Raman数据。(e-h) Ag和Tg声子的角强度分布，符合Raman张量预期。(i,j) 3.6 meV磁振子的二重旋转对称性。(k,l) 2.5 meV Stokes侧磁振子的四重+二重卷积。(m,n) 2.5 meV anti-Stokes侧的90度旋转图案。

磁振子对称性破缺：从四重到二重

Tg对称性的磁振子在传统反铁磁体中应呈现四重旋转对称性。但在MnTe2中，2.5 meV磁振子出现明显的二重成分，3.6 meV磁振子甚至完全呈现二重对称性。这种对称性降低是交错磁序的直接结果。

d波对称性：Stokes和anti-Stokes侧的角分布图案旋转90度，对应d波序参量在动量空间中的符号反转。RL和LR偏振分别探测相反的自旋子晶格，在d波节点处产生强度最大值。

温度依赖性：二重各向异性随温度升高而减弱，在T=90K（略高于TN=87K）时恢复为近四重对称的准弹性自旋涨落。这确证了对称性破缺与磁有序的直接关联。

图 3：磁振子能带和自旋纹理。(a) Brillouin区高对称点。(b) (-X)-Gamma-X-M-Gamma-(-M)路径上的磁振子色散，按自旋z分量着色。(c) (-R)-Gamma-R路径上的磁振子色散，强自旋z纹理。(c1-4) 各支能带的自旋纹理分别展示。(d) kz*a=+/-0.1pi平面上E1-3磁振子能带的自旋纹理，展示plaid-like（格子花纹）劈裂结构。

Plaid-like自旋劈裂：磁振子的"交错磁性指纹"

Plaid-like（格子花纹）自旋劈裂是交错磁性磁振子的标志性特征：自旋纹理满足在k_i->-k_i下为偶函数，在k_j->-k_j（j!=i）下为奇函数。这与先前在MnTe2电子能带中观测到的plaid-like劈裂完全相同（Zhu et al., Nature 2024）。

对称性根源：复合操作2\|\|{2\|a/2(e_i+e_j)}--自旋绕y轴旋转pi + 晶格绕y轴旋转pi + 平移(0, 0.5a, 0.5a)--是MnTe2磁基态的对称性。这一操作直接导致动量空间中自旋纹理的奇偶交替。

物理意义：磁振子和电子能带共享相同的自旋纹理对称性，说明交错磁性是体系的固有属性，不依赖于具体激发类型。这为通过Raman光谱（探测磁振子）间接推断电子能带的自旋劈裂提供了新途径。

DFT Tips

【DFT Tip 1】DFT+U中Ueff的选择：Dudarev vs cFLL+U+J

本文使用Dudarev方案（Ueff=U-J=4.2 eV），仅依赖有效U参数。同时对比了cFLL+U+J方案（U=5.0 eV, J=0.8 eV），发现|D1|/J1比值变化约21%。

两种方案的区别：Dudarev方案简单（一个参数），适合定性趋势分析；cFLL+U+J方案分别控制U和J，可以更精确地描述轨道选择性，但需要两个参数。

选择建议：对于Mn2+（d5高自旋，轨道角动量淬灭），Dudarev方案通常足够。对于含轨道自由度的体系（如Co2+, Fe2+），cFLL+U+J更合适。无论哪种方案，必须测试U值在+/-1-2 eV范围内的鲁棒性。

【DFT Tip 2】磁力线性响应 vs 能量映射法：提取交换参数

能量映射法（energy mapping）：构建FM、AFM1、AFM2等不同磁构型，计算总能量差，解线性方程组拟合J。优点：直观，任何DFT代码可做。缺点：需要超胞，难以提取DM矢量和各向异性。

磁力线性响应（Jx代码）：直接在DFT基态上通过Green函数计算自旋旋转的线性响应。优点：不需要超胞，可同时提取J/D/Gamma；缺点：需要全带Hamiltonian（通常用OpenMX或Wannier90接口），对k点密度敏感。

本文的验证：J1=0.675 meV, J2=-0.029 meV与文献中的能量映射结果（J1~7.1K=0.61 meV, J2~-1.3K=-0.11 meV）定性一致，验证了磁力线性响应的可靠性。

【DFT Tip 3】非共线磁序的DFT设置：MAGMOM与SAXIS

MnTe2的4个Mn子晶格自旋沿不同111方向，在VASP中需要设置非共线计算（LNONCOLLINEAR=.TRUE.）和每个原子的初始磁矩方向（MAGMOM）。

INCAR关键设置：LNONCOLLINEAR=.TRUE., MAGMOM = 0 0 5 0 0 0 0 0 0 0 0 0（仅Mn原子有磁矩）。具体的自旋方向通过SAXIS或直接设置MAGMOM的x,y,z分量。

常见陷阱：(1) 忘记设置LNONCOLLINEAR，VASP默认共线计算，非共线磁序会被错误地投影到z轴；(2) 初始磁矩方向设置错误，可能收敛到亚稳态；(3) 非共线计算中LORBIT=11可以输出每个原子的轨道和自旋磁矩方向。

【DFT Tip 4】Holstein-Primakoff变换的适用条件与局限

HP变换将自旋算符近似为玻色子，适用于低温（T<>1/2）极限。对于S=5/2的Mn2+，1/S修正约为20%，线性自旋波理论是合理的。

局限：(1) 忽略了磁振子-磁振子相互作用（非线性效应），在高温或接近TN时失效；(2) 对于S=1/2体系（如Cu2+），量子涨落强，线性自旋波理论精度差；(3) 对于强阻挫体系，可能需要考虑超越线性近似的量子修正。

Colpa对角化：处理非对角玻色子Bogoliubov-de Gennes方程的标准方法。注意：只有正的本征值才是物理的磁振子能量，负本征值表示选择的磁基态不稳定。

【DFT Tip 5】DM相互作用的DFT提取：Moriya规则与对称性

DM矢量D_ij满足Moriya对称性规则：如果两个磁性离子之间存在反演中心，DM=0；如果存在镜面，DM垂直于镜面；如果存在二重旋转轴，DM垂直于旋转轴。

MnTe2中D1!=0（最近邻Mn-Mn对无反演中心），D2=0（次近邻Mn-Mn对存在反演中心）。这是Moriya规则的直接验证。

DFT中提取DM的方法：(1) 能量映射法：计算非共线磁构型的能量，包含SOC；(2) 磁力线性响应：直接计算自旋-轨道耦合对交换相互作用的修正。两种方法都需要包含SOC。注意：SOC计算中k点需要更密，因为DM矢量通常比J小一个数量级。

【DFT Tip 6】交错磁体磁振子的对称性分析：自旋空间群

传统磁空间群（magnetic space group）无法描述交错磁体的对称性，需要自旋空间群（spin space group）。自旋空间群允许自旋空间和实空间的独立操作，包括复合操作如2\|\|{2\|t}。

MnTe2的对称性操作2\|\|{2\|a/2(e_i+e_j)}：自旋绕y轴旋转pi，同时晶格绕y轴旋转pi并平移(0, 0.5a, 0.5a)。这一操作导致磁振子自旋纹理的奇偶交替--plaid-like pattern。

实践建议：在分析交错磁体时，先确定自旋空间群（可用FindSpinGroup或自定义脚本），然后推导自旋纹理的对称性约束。如果实验结果与对称性预测不一致，检查磁构型是否正确。

【DFT Tip 7】Raman光谱与DFT的互补性：从实验到理论

Raman光谱探测的是Gamma点（q=0）的磁振子，而DFT+自旋波理论可以计算整个Brillouin区的磁振子色散。两者的互补性：(1) Raman提供Gamma点磁振子能量的精确实验值；(2) DFT提供整个色散关系和自旋纹理。

Raman强度的DFT计算：需要计算Raman张量（磁振子的光耦合强度），这涉及自旋-声子耦合和磁光耦合的复杂计算，比能带计算困难得多。本文未进行Raman强度的第一性原理计算，而是专注于磁振子色散和自旋纹理。

交叉验证：Raman实验观测到2.5 meV和3.6 meV两个磁振子，DFT预测Gamma点E4=2.40 meV, E1=E2=E3=3.54 meV。定性一致，但定量差异约4-5%，可接受。

【DFT Tip 8】OpenMX在磁性体系中的优势：全带Hamiltonian

OpenMX使用赝原子轨道（PAO）基组，天然适合构建紧束缚型Hamiltonian。对于磁力线性响应计算，需要完整的Bloch波函数在实空间Wannier基中的表示，OpenMX的PAO基组可以直接提供。

VASP vs OpenMX的工作分工：VASP做结构弛豫和基态电子结构（PAW基组精度高），OpenMX做全带Hamiltonian构建和磁力计算（PAO基组便于后续处理）。这是磁性体系DFT计算的典型工作流。

基组选择：Mn s3p2d2f1, Te s3p3d2f1。f极化轨道对Mn（3d）的影响较小（<5%），但对Te（5p）有一定影响（~10%）。9x9x9 k点和1000 Ry截断是OpenMX的标准设置。

知识扩展

【知识扩展 1】交错磁体（Altermagnet）的分类与识别

【理论解释】交错磁性是自旋群（spin group）对称性分类的产物。根据自旋劈裂的角动量形式，分为d波、g波、i波交错磁体。d波（l=2）是最常见的，如MnTe（六方）和RuO2（四方）。MnTe2（立方）中的plaid-like劈裂实际上是一种更高阶的d波形式。

【识别方法】实验上识别交错磁体的方法：(1) ARPES（角分辨光电子能谱）直接观测电子能带自旋劈裂；(2) 中子散射观测磁振子手性劈裂；(3) Raman光谱观测磁振子圆二色性（本文方法）；(4) 输运测量（反常Hall效应、自旋Hall效应）。

【方法比较】ARPES最直接但需要大单晶和同步辐射；中子散射需要大样品和反应堆/散裂源；Raman光谱是桌面级方法，样品需求量小，但只能探测Gamma点磁振子。三种方法互补。

【经典参考】Smejkal et al., PRX 12, 031042 & 040501 (2022)；Krempasky et al., Nature 626, 517 (2024)；Zhu et al., Nature 626, 523 (2024)；Liu et al., PRX 12, 021016 (2022)。

【迁移能力】交错磁体概念适用于所有磁性体系。关键问题是：你的体系的自旋群是否允许非相对论自旋劈裂？可以使用在线工具或FindSpinGroup代码检查。

【知识扩展 2】线性自旋波理论：从海森堡模型到磁振子

【理论解释】自旋波是磁性有序态中的低能集体激发。在铁磁体中，自旋波色散为omega~k^2（无能隙）；在反铁磁体中，omega~k（线性色散，无能隙）。交错磁体中的自旋波继承了能带的自旋劈裂特征，出现手性劈裂。

【方法比较】计算磁振子色散的方法：(1) 线性自旋波理论（本文，低温，大S极限）；(2) 经典Monte Carlo + 动力学结构因子（任意温度，但）；(3) TDDFT（第一性原理，但，仅适用于小体系）；(4) 非弹性中子散射（实验）。

【Holstein-Primakoff vs Dyson-Maleev】HP变换保留玻色子对易关系但产生非对角项，Dyson-Maleev变换产生对角Hamiltonian但破坏厄米性。对于非共线磁序，HP+Colpa是标准选择。

【经典参考】Holstein & Primakoff, Phys. Rev. 58, 1098 (1940)；Colpa, Physica A 93, 327 (1978)；Toth & Lake, J. Phys.: Condens. Matter 27, 166002 (2015)--SpinW教程。

【迁移能力】线性自旋波理论适用于任何已知交换参数的磁性体系。对于S=1/2体系，建议用量子Monte Carlo或精确对角化验证。

科研经验

【科研经验 1】DFT提取交换参数的不确定性：三种方法的系统比较

问题：我用能量映射法、磁力线性响应和Green函数方法得到的J值各不相同，该信哪个？

原因：(1) 能量映射法依赖于超胞大小和磁构型选择，J3和更远的交换可能被错误地"折叠"到J1和J2中；(2) 磁力线性响应依赖于k点密度和能量窗口；(3) 不同方法对U值的敏感度不同。

解决方案：(1) 同时使用至少两种独立方法验证J值；(2) 对于能量映射法，使用多个超胞大小验证收敛性；(3) 对于磁力线性响应，测试k点从8x8x8到12x12x12的收敛性；(4) 报告J值的误差范围而非单一数值。

建议：本文的|D1|/J1比值在不同U值下变化<20.3%，不同U方案下变化<21%，这是合理的误差范围。如果J值变化超过50%，说明体系对U值敏感，需要更精确的交换关联泛函（如HSE06）。

【科研经验 2】Raman光谱中磁振子信号的识别：避免误判

问题：Raman光谱中低能区（<10 meV）的信号很多，如何区分磁振子和声子/杂质荧光？

原因：低能区可能存在：(1) 声学声子（但Gamma点声学声子能量为零）；(2) 杂质或缺陷的局域振动模式；(3) 准弹性散射（Rayleigh尾）；(4) 荧光背景。这些都可能被误判为磁振子。

解决方案：(1) 温度依赖性--磁振子强度在TN以下出现，在TN以上消失（或转变为顺磁振子）；(2) 偏振依赖性--声子遵循Raman张量，磁振子可能违反；(3) 磁场依赖性--外加磁场使磁振子Zeeman劈裂，声子不受影响；(4) 手性分辨--手性磁振子显示RL/LR不对称，声子通常不显示。

建议：本文的2.5 meV信号同时满足所有四条判据，确认为磁振子。如果你只看到其中一条，需要更多证据。

如果是我，我还会继续算

【继续算 1】非线性自旋波理论：磁振子-磁振子相互作用

为什么值得算：线性自旋波理论忽略了磁振子-磁振子相互作用。在高温下（接近TN），非线性效应变得重要。1/S展开（下一阶~1/S^2）可以给出磁振子寿命和能量重整化。

能回答的问题：非线性修正对Gamma点磁振子能量（2.5 meV和3.6 meV）有多大影响？磁振子寿命随温度如何变化？

适合体系：所有磁性体系。输入：交换参数+1/S展开代码（如SpinW的扩展）。

【继续算 2】Wannier90 + 电子能带自旋纹理：与磁振子对比

为什么值得算：本文展示了磁振子的plaid-like自旋纹理，但未直接计算电子能带的自旋纹理。Wannier90可以构建电子能带的紧束缚模型，直接计算电子自旋纹理，与磁振子自旋纹理进行定量对比。

能回答的问题：电子能带和磁振子能带的自旋纹理是否完全一致？是否存在定量差异（如劈裂幅度）？

适合体系：所有交错磁体。输入：Wannier90（Mn-d, Te-p轨道）。

【继续算 3】Monte Carlo模拟：磁转变温度与临界行为

为什么值得算：DFT不能直接预测TN=87K。使用提取的交换参数进行经典Monte Carlo模拟，可以预测TN和临界指数，与实验值直接对比。

能回答的问题：DFT+MC预测的TN与实验值（87K）是否一致？DM相互作用对TN的贡献有多大？

适合体系：所有磁性体系。输入：EspinS/UppASD（Heisenberg+DM模型）到中等。

【继续算 4】Berry曲率与反常Hall效应：交错磁体的拓扑输运

为什么值得算：交错磁体中的非零Berry曲率多极矩可以产生反常Hall效应。MnTe2的非共面磁序+plaid-like自旋纹理可能产生独特的拓扑输运性质。

能回答的问题：MnTe2的AHE有多大？AHE是否与磁振子手性劈裂有关联？

适合体系：所有交错磁体。输入：Wannier90+WannierTools到中等。

【继续算 5】应变工程：调控磁振子手性劈裂

为什么值得算：应变可以改变晶格对称性，从而改变自旋空间群的对称性。压缩或拉伸应变可能增强或抑制磁振子的手性劈裂，为器件应用提供调控手段。

能回答的问题：-2%到+2%应变下磁振子手性劈裂如何变化？是否存在对称性破缺导致的拓扑相变？

适合体系：所有2D和薄膜磁性体系。输入：应变下的DFT+U + 磁力线性响应 + 自旋波理论。

【继续算 6】自旋波输运：磁振子自旋流与非互易传播

为什么值得算：手性磁振子意味着omega(k)!=omega(-k)，这直接导致磁振子的非互易传播--自旋波沿+k和-k方向的传播速度不同。这对于磁振子逻辑器件至关重要。

能回答的问题：MnTe2中磁振子的群速度各向异性有多大？非互易比（v(+k)/v(-k)-1）是多少？

适合体系：所有手性磁振子体系。输入：自旋波色散+群速度计算。

Wulferding, An, Choi et al. | Adv. Sci. (2026) e76555 | 交错磁体 手性磁振子 Raman光谱 自旋波理论