交错磁体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波模式方法体系VASPLDAU基态计算 - 自研DFPT代码横向自旋磁化率计算 - 对称性分析自旋空间群[C2||A]约束 - 实验验证非弹性中子散射对比。关键发现(1) 手性磁振子谱在四个节面qz0平面上简并在Brillouin区边界处能量相同但强度不同在其他区域完全劈裂(2) 磁振子谱的等能面呈现g波空间模式六次符号翻转(3) 对称性分析表明磁振子劈裂源于连接相反自旋子晶格的操作是旋转而非反演或半格平移。Alpha-MnTe磁振子谱DFPT研究流程。VASP LDAU基态电子结构、对称性分析 - DFPT横向自旋磁化率计算Sternheimer方程 - Lorentz拟合三次样条插值 - 手性磁振子劈裂分析 g波模式验证。核心发现磁振子手性劈裂源于反演/半格平移对称性的破缺。二、研究方法DFPT计算横向自旋磁化率原理与实现DFPT密度泛函微扰理论是计算材料对外场线性响应的标准方法。在磁振子问题中需要计算横向自旋磁化率Im[Chi-(omega,q)]和Im[Chi-(omega,q)]它们分别对应右旋和左旋手性磁振子。计算设置LDA泛函 PAW赝势, ENCUT500 eV, 9x9x4 Gamma-centered k点, 实验晶格常数a4.134 A, c6.652 A。Mn-3d轨道加Hubbard Ueff4.0 eV。频率范围-40 meV到40 meV间隔5 meV在谱峰附近细化到1 meV。关键技术点Sternheimer方程中展宽参数eta1 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波交错磁性特征。在kz0的节面上自旋完全简并。在kz0.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-MIm[Chi-]Im[Chi-]手性磁振子完全简并。Gamma点磁振子能量为零Goldstone定理M点能量为31.5 meV与中子散射实验一致。模式二BZ边界等能异强沿L-A-L手性磁振子能量相等但强度不同。这是因为omega_Romega_LBZ周期性但自旋磁化率缺乏周期性。模式三交替手性劈裂沿N-Gamma-N手性磁振子在能量和强度上都劈裂且劈裂方向沿qz方向交替。模式四对称劈裂沿Lbar-Gamma-L手性磁振子对称劈裂因为MnTe保留反演对称性PIm[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)。对于MnTeA可以是四个镜面操作之一{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)因此qz0平面上的磁振子手性简并。关键的物理洞察在传统反铁磁体中A是反演或半格平移使得Aqq对所有q成立因此磁振子永远手性简并。在AM中A只能是旋转操作Aq不等于q对所有q成立因此手性劈裂出现。图 4qz-1.33非节面上的等能面切片。(a1-a3) omega28 meV, (b1-b3) omega32 meV, (c1-c3) omega34 meV。(a1-c1) 总磁振子谱Im[Chi-Chi-]。(a2-c2) 手性差异谱Im[Chi--Chi-]。(a3-c3) 手性磁振子能带等能面。g波磁振子模式的实验验证在omega32 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_M1/2是最实用的选择。【DFT Tip 2】LDA vs GGA在MnTe中的选择本文使用LDA而非GGAPBE这是关键选择。MnTe是窄带隙半导体PBE通常低估带隙但高估晶格常数。LDA虽然低估晶格常数但对MnTe的电子结构描述更准确。作者使用了实验晶格常数a4.134 A, c6.652 A而非DFT优化的晶格常数。这是常见的做法--当DFT优化的晶格常数与实验偏差较大时用实验晶格常数做电子结构计算更可靠。常见错误用PBE优化晶格常数后直接做PBEU计算晶格常数误差可能传播到电子结构和磁性质。建议始终对比实验晶格常数和DFT优化晶格常数下的结果。【DFT Tip 3】Hubbard U对MnTe磁振子谱的影响本文使用Ueff4.0 eVMn-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值或者至少扫描U2-6 eV范围比较磁振子能量的变化趋势。【DFT Tip 4】Sternheimer方程中的展宽参数etaDFPT中Sternheimer方程的展宽参数eta1 meV是平衡精度和数值稳定性的关键。eta太小会导致发散分母趋近零eta太大导致谱峰展宽过度丢失精细结构。对于绝缘体/半导体eta可以设置得较小0.1-1 meV。对于金属eta需要更大10-50 meV以模拟电子散射。MnTe是半导体eta1 meV是合理选择。在磁振子谱峰附近频率间隔从5 meV细化到1 meV以确保精确确定磁振子能量。这是数值精确性的关键细节在你的计算中也应遵循。【DFT Tip 5】k点采样对磁振子DFPT计算的影响本文使用9x9x4 Gamma-centered k点网格对于MnTe的六角晶胞来说不算特别密集。但DFPT的计算量随k点数量线性增长9x9x4324个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点手性劈裂。对于MnTeA{m001|0 0 1/2}镜面操作Aq(qx,qy,-qz)。在qz0平面上Aqq手性简并。在其他q点Aq不等于q手性劈裂。你可以用这个简单规则快速判断任何AM材料找到磁空间群中连接相反自旋子晶格的操作检查它是否只包含旋转/镜面不含反演/平移然后判断哪些q点满足Aqq。知识扩展【知识扩展 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代码DFPTU声子。【迁移能力】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建议设置EDIFF1e-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_M1/2需构建不同自旋螺旋构型的超胞。【继续算 2】磁振子线宽与寿命从DFPT到多体微扰为什么值得算本文计算的磁振子谱是无限寿命的eta1 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波