交错磁体Alpha-MnTe的第一性原理磁振子谱研究

PHYS. REV. B 113, 144301 (2026)

交错磁体Alpha-MnTe的第一性原理磁振子谱研究

First-Principles Study of the Magnon Spectrum of Altermagnetic MnTe

导读 导读：交错磁体（Altermagnet）的电子能带劈裂已被广泛研究，但其集体自旋激发--磁振子--的DFT研究一直受困于方法学挑战。本文通过自主开发的DFPT代码，首次从第一性原理计算了Alpha-MnTe的完整磁振子谱。核心发现：手性磁振子在四个节面上简并，在Brillouin区边界等能异强，在其他区域完全劈裂，磁振子谱的等能面呈现g波空间模式。对称性分析揭示，磁振子手性劈裂源于连接相反自旋子晶格的旋转操作（而非反演/平移），这是AM区别于传统反铁磁体的本质特征。

一、前言背景

交错磁体中的磁振子：从电子到集体激发

交错磁体（Altermagnet, AM）自2022年概念提出以来，大量研究集中在电子结构（能带劈裂、自旋-动量锁定）和输运性质（反常霍尔效应、自旋霍尔效应）。然而，磁振子（magnon）--自旋波的量子化准粒子--作为磁性材料中最重要的集体激发，在AM中的研究才刚刚起步。

Alpha-MnTe是AM的典型代表材料。它具有六角NiAs结构（P63/mmc空间群），Neel温度约310 K。实验上，ARPES已确认其g波形式的电子能带劈裂。但磁振子谱的DFT研究一直受限于方法学困难--传统DFT难以直接计算磁振子。

本文作者利用其自主开发的DFPT代码（基于VASP 5.4.4），首次从第一性原理计算了Alpha-MnTe的完整磁振子谱，揭示了手性磁振子的劈裂行为和g波空间模式，并与非弹性中子散射实验结果高度吻合。

DFPT磁振子计算：方法学突破

核心问题：在AM中，手性磁振子（左旋和右旋）是否劈裂？劈裂的对称性根源是什么？磁振子谱是否也呈现与电子能带相同的g波模式？

方法体系：VASP（LDA+U基态计算） -> 自研DFPT代码（横向自旋磁化率计算） -> 对称性分析（自旋空间群[C2||A]约束） -> 实验验证（非弹性中子散射对比）。

关键发现：(1) 手性磁振子谱在四个节面（qz=0平面）上简并，在Brillouin区边界处能量相同但强度不同，在其他区域完全劈裂；(2) 磁振子谱的等能面呈现g波空间模式（六次符号翻转）；(3) 对称性分析表明，磁振子劈裂源于连接相反自旋子晶格的操作是旋转而非反演或半格平移。

Alpha-MnTe磁振子谱DFPT研究流程。VASP LDA+U（基态电子结构、对称性分析） -> DFPT横向自旋磁化率计算（Sternheimer方程） -> Lorentz拟合+三次样条插值 -> 手性磁振子劈裂分析 + g波模式验证。核心发现：磁振子手性劈裂源于反演/半格平移对称性的破缺。

二、研究方法

DFPT计算横向自旋磁化率：原理与实现

DFPT（密度泛函微扰理论）是计算材料对外场线性响应的标准方法。在磁振子问题中，需要计算横向自旋磁化率Im[Chi+-(omega,q)]和Im[Chi-+(omega,q)]，它们分别对应右旋和左旋手性磁振子。

计算设置：LDA泛函 + PAW赝势, ENCUT=500 eV, 9x9x4 Gamma-centered k点, 实验晶格常数a=4.134 A, c=6.652 A。Mn-3d轨道加Hubbard Ueff=4.0 eV。频率范围-40 meV到40 meV，间隔5 meV，在谱峰附近细化到1 meV。

关键技术点：Sternheimer方程中展宽参数eta=1 meV（防止发散），自洽求解一阶波函数和诱导密度，然后通过线性响应理论计算自旋磁化率。

Sternheimer方程：一阶波函数delta u_nk在外场下的响应，eta为展宽参数

横向自旋磁化率Im[Chi+-]的Lorentz拟合形式：omega_R/L为手性磁振子能量，R/L为振幅

对称性约束：存在[C2||A]对称性时，Im[Chi+-(omega,q)] = Im[Chi-+(omega,Aq)]

Kohn-Sham横向自旋磁化率：只有自旋翻转的态到态跃迁有贡献

三、核心结果

图 1：(a) Alpha-MnTe晶体结构（NiAs型，P63/mmc）。(b) 第一Brillouin区和四个节面（黄色和绿色）。(c,d) 沿高对称路径Mbar-Gamma-M和N' -Gamma' -N的电子能带结构。(e,f) E=-1eV时的自旋等能面，展示g波空间模式。

电子能带结构：g波交错磁性的DFT确认

沿高对称路径Mbar-Gamma-M，自旋保持简并（这是节面方向）。沿非高对称路径N'-Gamma'-N，自旋交替极化清晰可见。这确认了MnTe的g波交错磁性特征。

在kz=0的节面上，自旋完全简并。在kz=0.3的非节面上，自旋极化六次翻转，呈现g波空间对称性。这与ARPES实验结果一致，验证了MnTe作为g波AM的经典地位。

图 2：沿不同路径的磁振子谱。Im[Chi+-]（红色）和Im[Chi-+]（蓝色），颜色强度由|Im[Chi+-]-Im[Chi-+]|决定。(a) Mbar-Gamma-M: 手性简并。(b) L'-A-L: 强度略有差异但能量相等。(c) N'-Gamma'-N: 交替手性劈裂。(d) Lbar-Gamma-L: 对称劈裂。

磁振子谱的四种行为模式

模式一（节面简并）：沿Mbar-Gamma-M，Im[Chi+-]=Im[Chi-+]，手性磁振子完全简并。Gamma点磁振子能量为零（Goldstone定理），M点能量为31.5 meV，与中子散射实验一致。

模式二（BZ边界等能异强）：沿L'-A-L，手性磁振子能量相等但强度不同。这是因为omega_R=omega_L（BZ周期性），但自旋磁化率缺乏周期性。

模式三（交替手性劈裂）：沿N'-Gamma'-N，手性磁振子在能量和强度上都劈裂，且劈裂方向沿qz方向交替。

模式四（对称劈裂）：沿Lbar-Gamma-L，手性磁振子对称劈裂，因为MnTe保留反演对称性P：Im[Chi+-(omega,q)]=Im[Chi+-(omega,Pq)]。

图 3：一维恒定Q切割的磁振子谱。M=(0.5,0,0), N=(0.25,0,0.5), N'=(0.25,0,-0.5), O=(0.25,0,0.25), O'=(0.25,0,-0.25)。蓝色和红色分别为Im[Chi+-]和Im[Chi-+]。圆圈/星号为计算数据，实线/虚线为Lorentz拟合。

对称性分析的严格推导

从Kohn-Sham磁化率出发，利用AM的对称性[C2||A]（C2是垂直于自旋的180度旋转，A是仅作用于实空间的变换），严格推导了Im[Chi+-(omega,q)]=Im[Chi-+(omega,Aq)]，即omega_R(q)=omega_L(Aq)。

对于MnTe，A可以是四个镜面操作之一：{m001|0 0 1/2}, {m1-10|0 0 1/2}, {m120|0 0 1/2}, {m210|0 0 1/2}。以{m001|0 0 1/2}为例，Aq=(qx,qy,-qz)，因此qz=0平面上的磁振子手性简并。

关键的物理洞察：在传统反铁磁体中，A是反演或半格平移，使得Aq=q对所有q成立，因此磁振子永远手性简并。在AM中，A只能是旋转操作，Aq不等于q对所有q成立，因此手性劈裂出现。

图 4：qz=-1.33非节面上的等能面切片。(a1-a3) omega=28 meV, (b1-b3) omega=32 meV, (c1-c3) omega=34 meV。(a1-c1) 总磁振子谱Im[Chi+-+Chi-+]。(a2-c2) 手性差异谱Im[Chi+--Chi-+]。(a3-c3) 手性磁振子能带等能面。

g波磁振子模式的实验验证

在omega=32 meV和34 meV时，总磁振子谱Im[Chi+-+Chi-+]出现六重高强度节点，与手性磁振子能带的简并点位置一致。手性差异谱Im[Chi+--Chi-+]呈现六次符号翻转，确认为g波模式。

橙色虚线（节面处）在所有能量下保持不变，是手性磁振子的对称性保护简并线。绿色实线随能量变化，是手性磁振子的偶然简并线。

本文结果与中子散射实验（Liu et al., PRL 133, 156702, 2024）高度吻合，验证了DFPT方法在磁振子计算中的可靠性。

图 5：本文计算结果与实验数据（Ref.[32]）的直接对比。热图为实验数据，标记点为本文计算数据。

DFT Tips

【DFT Tip 1】DFPT磁振子计算：VASP内置功能 vs 自研代码

标准VASP不直接支持DFPT磁振子计算。本文使用作者自研的DFPT代码（基于VASP 5.4.4），实现了横向自旋磁化率的自洽计算。如果你需要计算磁振子谱，有三个选择：

(1) Frozen Magnon方法：构建不同自旋螺旋构型的超胞，计算总能量差，精度高但。(2) 四态力方法（four-state force method）：从DFT力常数提取交换参数，再解自旋波方程。(3) 本文的DFPT方法：直接计算自旋磁化率，精度最高但需要自研代码。

如果只需要磁振子色散关系（而非手性劈裂），Frozen Magnon方法（VASP I_CONSTRAINED_M=1/2）是最实用的选择。

【DFT Tip 2】LDA vs GGA在MnTe中的选择

本文使用LDA而非GGA（PBE），这是关键选择。MnTe是窄带隙半导体，PBE通常低估带隙但高估晶格常数。LDA虽然低估晶格常数，但对MnTe的电子结构描述更准确。

作者使用了实验晶格常数（a=4.134 A, c=6.652 A），而非DFT优化的晶格常数。这是常见的做法--当DFT优化的晶格常数与实验偏差较大时，用实验晶格常数做电子结构计算更可靠。

常见错误：用PBE优化晶格常数后直接做PBE+U计算，晶格常数误差可能传播到电子结构和磁性质。建议始终对比实验晶格常数和DFT优化晶格常数下的结果。

【DFT Tip 3】Hubbard U对MnTe磁振子谱的影响

本文使用Ueff=4.0 eV（Mn-3d），这是MnTe文献中最常用的值（Krempasky et al., Nature 2024; Gonzalez Betancourt et al., PRL 2023）。U值直接影响d轨道能级位置和磁矩大小。

注意：U值的选择会影响磁振子能量。U增大通常导致磁矩增大、交换耦合增强、磁振子硬化。在论文中，作者应与实验磁振子能量（如中子散射31.5 meV at M点）进行对比，验证U值的合理性。

如果你要计算不同AM材料的磁振子谱，建议先用线性响应方法确定U值，或者至少扫描U=2-6 eV范围，比较磁振子能量的变化趋势。

【DFT Tip 4】Sternheimer方程中的展宽参数eta

DFPT中Sternheimer方程的展宽参数eta=1 meV是平衡精度和数值稳定性的关键。eta太小会导致发散（分母趋近零），eta太大导致谱峰展宽过度，丢失精细结构。

对于绝缘体/半导体，eta可以设置得较小（0.1-1 meV）。对于金属，eta需要更大（10-50 meV）以模拟电子散射。MnTe是半导体，eta=1 meV是合理选择。

在磁振子谱峰附近，频率间隔从5 meV细化到1 meV，以确保精确确定磁振子能量。这是数值精确性的关键细节，在你的计算中也应遵循。

【DFT Tip 5】k点采样对磁振子DFPT计算的影响

本文使用9x9x4 Gamma-centered k点网格，对于MnTe的六角晶胞来说不算特别密集。但DFPT的计算量随k点数量线性增长，9x9x4=324个k点已经足够捕获磁振子色散的主要特征。

如果你的体系晶胞更小（如立方钙钛矿），k点可以更密。如果晶胞更大（如超胞），k点可以适当减少。关键原则是：磁振子色散的最高能量点（如BZ边界M点）应该收敛到1 meV以内。

常见陷阱：用Monkhorst-Pack网格（而非Gamma-centered）做DFPT计算，导致Gamma点附近磁振子色散不准确。因为磁振子在Gamma点趋于零（Goldstone模式），Gamma-centered网格更精确。

【DFT Tip 6】自旋磁化率的Lorentz拟合技巧

DFPT直接输出的是Im[Chi+-(omega,q)]的离散数据点。为获得高分辨率谱，作者用Lorentz函数拟合，提取磁振子能量omega_R/L和振幅R/L，然后做三次样条插值重构谱。

Lorentz拟合假设每个磁振子峰是孤立且对称的。当两个磁振子峰重叠时，需要双Lorentz（或更多）拟合。拟合质量可以通过残差分析验证。

在你的计算中，如果磁振子谱出现不对称峰形（如Fano线形），说明磁振子与电子-空穴对连续谱耦合，此时单Lorentz拟合不适用。

【DFT Tip 7】DFPT磁振子计算中的频率范围选择

本文频率范围-40 meV到40 meV（约-320到320 cm^-1），覆盖了MnTe磁振子的整个色散范围（最大能量约31.5 meV at M点）。负频率对应反斯托克斯过程（磁振子吸收）。

频率步长5 meV是粗扫描，1 meV是精扫描。如果你的体系磁振子能量更高（如氧化物反铁磁体可达100 meV），需要相应扩大频率范围。

注意：DFPT磁振子计算的CPU时间与频率点数成正比，但内存与k点数成正比。先做粗扫描确定磁振子能量范围，再精细扫描。

【DFT Tip 8】AM磁振子手性劈裂的对称性判断

判断一个AM材料是否会出现磁振子手性劈裂，只需要检查连接相反自旋子晶格的对称操作A。如果A是反演或半格平移（传统AFM），则磁振子永远手性简并。如果A是旋转操作（AM），则磁振子在Aq不等于q的q点手性劈裂。

对于MnTe，A={m001|0 0 1/2}镜面操作，Aq=(qx,qy,-qz)。在qz=0平面上Aq=q，手性简并。在其他q点，Aq不等于q，手性劈裂。

你可以用这个简单规则快速判断任何AM材料：找到磁空间群中连接相反自旋子晶格的操作，检查它是否只包含旋转/镜面（不含反演/平移），然后判断哪些q点满足Aq=q。

知识扩展

【知识扩展 1】DFPT方法的发展历程与应用

【理论解释】DFPT由Baroni、Giannozzi和Testa于1987年提出（PRL 58, 1861），最初用于计算声子谱。其核心思想是：不直接计算体系对外场的响应函数，而是通过Sternheimer方程自洽求解一阶波函数，然后通过线性响应理论得到任意响应函数。

【方法比较】DFPT vs 有限差分法：DFPT的计算量与k点数成正比，与扰动类型无关；有限差分法需要构建超胞，计算量与超胞大小的立方成正比。对于磁振子（长程交换耦合），有限差分法需要非常大的超胞，DFPT是更高效的选择。

【经典参考】Baroni et al., Rev. Mod. Phys. 73, 515 (2001) - DFPT综述；Liu et al., PRB 108, 094405 (2023) - 本文作者之前的DFPT代码论文；Cao et al., Comput. Phys. Commun. 286, 108680 (2023) - HP代码（DFPT+U声子）。

【迁移能力】DFPT方法不仅适用于磁振子，还适用于声子、介电函数、EPC、NMR化学位移等。掌握了DFPT，你就掌握了DFT线性响应计算的通用框架。

【知识扩展 2】磁振子的手性：从经典到量子

【理论解释】磁振子（magnon）是自旋波的量子化准粒子。手性磁振子（chiral magnon）指自旋进动方向不同的磁振子：右旋（顺时针）和左旋（逆时针）。在传统铁磁体中，只有一种手性的磁振子（与磁化方向一致）；在反铁磁体中，两种手性的磁振子简并。

【AM中的手性劈裂】在AM中，由于连接相反自旋子晶格的对称操作是旋转而非反演/平移，手性磁振子在能量上可以劈裂。这是AM磁振子区别于传统AFM磁振子的最显著特征。

【实验探测】手性磁振子可以通过极化中子散射（需极化分析）直接探测。本文与非极化中子散射的比较只能验证能量，手性的直接验证需要极化中子实验。

【经典参考】Smejkal et al., PRL 131, 256703 (2023) - RuO2手性磁振子理论；Liu et al., PRL 133, 156702 (2024) - MnTe磁振子实验；Sandratskii et al., PRB 111, 184436 (2025) - Frozen Magnon方法计算AM磁振子。

科研经验

【科研经验 1】DFPT磁振子计算中的收敛问题

问题：DFPT磁振子计算中，Sternheimer方程自洽迭代不收敛，尤其是在磁振子能量接近零的低频区域。

原因：Gamma点附近磁振子能量趋于零（Goldstone模式），Sternheimer方程的分母趋近零，导致数值不稳定。此外，展宽参数eta太小也会导致发散。

解决方案：(1) 适当增大eta（本文用1 meV）；(2) 在Gamma点附近增加频率采样密度；(3) 对于Gamma点附近的点，可以用更小的k点网格单独计算；(4) 检查自洽场收敛标准（EDIFF），建议设置EDIFF=1e-8或更严格。

建议：如果DFPT计算在Gamma点附近仍然不收敛，可以先跳过Gamma点，用外推法得到Gamma点的磁振子能量（理论值为零）。

【科研经验 2】AM磁振子研究的实验与理论对比

问题：DFPT计算的磁振子能量与中子散射实验对比时，如何判断一致性？

原因：DFT计算存在系统误差（带隙低估、U值不确定性），中子散射实验存在分辨率限制和背底噪声。两者直接比较时，2-3 meV的偏差在合理范围内。

解决方案：(1) 比较磁振子色散的形状和趋势，而非绝对值；(2) 比较不同q点的相对能量差；(3) 如果可能，用HSE06或GW校正带隙后再做DFPT；(4) 同时报告计算和实验的不确定度。

建议：DFPT磁振子计算目前仍处于发展阶段，将计算结果视为定性/半定量预测更合理。与实验比较时，关注色散形状、手性劈裂的存在与否、以及相对趋势，而非绝对值。

如果是我，我还会继续算

【继续算 1】Frozen Magnon方法验证DFPT结果

为什么值得算：DFPT是相对较新的磁振子计算方法，用更成熟的Frozen Magnon方法验证可以增强结果的可信度。两种方法相互独立，一致的结果是对磁振子色散的强有力确认。

能回答的问题：DFPT和Frozen Magnon给出的磁振子能量是否一致？手性劈裂是否在两种方法中都出现？

适合体系：所有磁性材料。输入：VASP I_CONSTRAINED_M=1/2，需构建不同自旋螺旋构型的超胞。

【继续算 2】磁振子线宽与寿命：从DFPT到多体微扰

为什么值得算：本文计算的磁振子谱是无限寿命的（eta=1 meV是人为展宽）。真实的磁振子有有限寿命，由磁振子-磁振子散射和磁振子-电子散射决定。计算磁振子线宽可以预测磁振子的输运性质。

能回答的问题：MnTe中磁振子的本征寿命是多少？不同q点的磁振子寿命有何差异？手性磁振子的寿命是否劈裂？

适合体系：所有磁性材料。输入：需要多体微扰理论（如GW近似下的Bethe-Salpeter方程）。

【继续算 3】磁振子热导率：DFPT到Boltzmann输运

为什么值得算：磁振子可以携带热流，磁振子热导率是自旋热电子学（spin caloritronics）的核心物理量。AM中手性磁振子的热输运可能具有独特的各向异性。

能回答的问题：MnTe的磁振子热导率有多大？手性磁振子对热导率的贡献是否不同？热导率是否有各向异性？

适合体系：所有磁性绝缘体/半导体。输入：从DFPT磁振子谱提取群速度和弛豫时间，输入Boltzmann输运方程。

【继续算 4】磁振子-声子耦合：DFPT的扩展

为什么值得算：在有限温度下，磁振子和声子可以通过非谐耦合交换能量，影响热输运和自旋弛豫。DFPT框架可以自然扩展到自旋-晶格耦合的计算。

能回答的问题：MnTe中哪些声子模式与磁振子耦合最强？磁振子-声子耦合是否影响磁振子手性劈裂？

适合体系：所有磁性材料。输入：需要计算自旋-声子耦合矩阵元（DFPT的扩展）。

【继续算 5】应变调控磁振子谱：AM的机械调控

为什么值得算：AM的磁振子手性劈裂由对称性决定，应变可以改变晶格对称性，从而调控磁振子谱。应变调控是AM器件应用的重要方向。

能回答的问题：单轴/双轴应变如何影响磁振子手性劈裂？是否存在临界应变使手性劈裂消失或反转？

适合体系：所有2D和层状AM材料。输入：不同应变下的DFPT磁振子计算。

【继续算 6】自旋Nernst效应：手性磁振子的输运特征

为什么值得算：手性磁振子的劈裂必然导致自旋Nernst效应（温度梯度驱动的横向自旋流）。计算自旋Nernst效应可以建立AM磁振子与自旋电子学器件之间的直接联系。

能回答的问题：MnTe的自旋Nernst系数有多大？手性磁振子对自旋Nernst效应的贡献是否可区分？

适合体系：所有AM材料。输入：磁振子Boltzmann输运 + 自旋Berry曲率。

Li, Liu, Qiao | Phys. Rev. B 113, 144301 (2026) | 交错磁体 磁振子 DFPT 手性劈裂 g波