PRB 交错磁体MnTe的多极矩分析与磁自旋霍尔效应PRB 113, 155104 (2026)交错磁体MnTe的多极矩分析与磁自旋霍尔效应Multipole Analysis of Spin Currents in Altermagnetic MnTe导读交错磁体中不同的磁构型产生不同秩的序参数这直接决定了自旋霍尔效应的存在性和各向异性。本文通过多极矩框架结合Kubo公式第一性原理计算系统分析了alpha-MnTe中N||y和N||x两种磁构型的电荷输运、反常霍尔效应和磁自旋霍尔效应。核心发现N||y构型磁偶极Mz同时具有AHE和磁SHE磁自旋霍尔角高达16%超过Pt的2倍N||x构型磁八极Mxyz只有磁SHE而无AHE。磁SHE的各向异性为实验判别序参数类型提供了实用方法。一、前言背景交错磁体MnTe从d波自旋劈裂到磁多极矩交错磁体Altermagnet是2022年以来凝聚态物理最受关注的发现之一。与传统反铁磁体由平移或反演连接反平行自旋不同交错磁体中反平行自旋由旋转或旋转-平移联合操作连接导致动量空间中产生非相对论自旋劈裂。alpha-MnTe是典型的g波交错磁体Neel温度~310 K具有六方NiAs型结构。ARPES实验已明确观测到其Kramers自旋简并的解除。本文考虑两种Neel矢量取向N||y实验可实现的磁构型和N||x尚未实现但理论预测的磁构型。本文的核心创新在于将多极矩multipole框架引入交错磁体的自旋输运分析建立了磁序参数与自旋霍尔效应SHE之间的直接对应关系并通过PAOFLOW Kubo公式进行了第一性原理验证。多极矩框架从对称性到输运响应的桥梁核心问题不同的磁构型N||y vs N||x如何影响自旋电流的产生磁自旋霍尔效应的各向异性如何反映序参数的类型方法体系VASPDFTUSOC基态电子结构 - PAOFLOW紧束缚投影构建Wannier-like Hamiltonian - Kubo公式线性响应电导率/自旋霍尔电导率计算 - 多极矩对称性分析。关键发现(1) N||y构型具有磁偶极矩Mz作为序参数同时出现AHE和磁SHE(2) N||x构型具有磁八极矩Mxyz作为序参数只有磁SHE而无AHE(3) 磁自旋霍尔角SHA高达16%超过Pt5-10%与beta-Ta12-15%相当。MnTe磁自旋霍尔效应研究流程。VASP DFTUSOC电子结构 - PAOFLOW紧束缚投影 - Kubo公式线性响应电导率、AHE、磁SHE、内禀SHE - 多极矩对称性分析。核心发现N||y构型同时具有AHE和磁SHESHA16%N||x构型只有磁SHESHA8.2%磁SHE各向异性可用于判别序参数类型。二、研究方法VASP DFTUSOC基态电子结构计算计算设置PBE泛函 PAW赝势, ENCUT350 eV, Gamma-centered 20x20x12 k点网格, 实验晶格参数a4.134 A, c6.652 A。SOC在所有计算中自洽包含。Hubbard U处理Mn-3d轨道采用DFTUDudarev方案U4.0 eV, J0.97 eV。该值在MnTe的常用范围内与GW计算和ARPES实验结果一致。注意U值的选择对MnTe的价带结构有显著影响但定性趋势对U值不敏感。两种磁构型(1) N||y磁点群mmm磁空间群Cmcm实验可实现(2) N||x磁点群mmm磁空间群Cmcm尚未实验实现但可通过应力-磁场组合诱导。PAOFLOW Kubo从DFT到输运响应PAOFLOW方法将VASP波函数投影到赝原子轨道PAO基组上构建紧束缚哈密顿量。Mn基组[1s,2s,2p,3s,3p,3d,4s,4p,4d,5s]Te基组[1s,2s,2p,3s,3p,3d,4s,4p,4d,5s,5p,5d,6s,6p]。通过Fourier插值零填充将哈密顿量插值到200x200x124的稠密k点网格。Kubo公式在恒定弛豫时间近似下计算六种线性响应耗散电导率T-even、AHET-odd、耗散REET-even、磁REET-odd、内禀SHET-even、磁SHET-odd。散射率eta0.01 eV由实验电导率确定。关键理解磁SHE是T-odd的耗散依赖响应由时间反演对称性破缺和交换相互作用驱动其物理起源与内禀SHET-evenSOC驱动完全不同。Kubo公式耗散贡献eta为散射率A和B为与物理量和外场对应的算符Kubo公式非耗散贡献在eta-0极限下独立于eta对应电场驱动的响应磁自旋霍尔角SHA无量纲量表征电荷-自旋转换效率beta_s(J),s_ij磁八极矩M_xyz最低阶交错磁序参数由电四极矩Q与磁偶极矩M_z的耦合产生三、核心结果图 1alpha-MnTe的晶体和磁结构。(a) 晶体结构NiAs型六方。(b) N||y磁构型磁矩沿y/-y方向交替排列。(c) N||x磁构型磁矩沿x/-x方向交替排列。图中标注了不可约对称操作。磁构型与多极矩序参数N||y构型磁点群mmm磁偶极矩Mz为最低阶序参数同时存在磁八极矩M_alpha_z和M_beta_z。这意味着该构型可以共享铁磁体沿z轴磁化的对称性无需额外对称性破缺。N||x构型磁点群mmm无磁偶极矩最低阶交错磁序参数为磁八极矩Mxyz。该构型中T-even和T-odd多极矩无法区分因为时间反演不是对称操作。多极矩分析的意义不同秩的序参数直接决定了响应张量的非零元素和对称性形式。例如磁偶极Mz允许AHEsigma_xy而磁八极Mxyz禁止AHE但允许特定分量的磁SHE。图 2自旋极化能带结构。(a) N||y构型的自旋极化能带颜色表示sz分量。(b) N||x构型的自旋极化能带。(c) 第一布里渊区和对称点。图 3自旋-动量锁定。(a) N||y构型在k_x-k_y平面上的sz(k)分布呈现sd波对称性Q0Qv主导。(b) N||x构型的sz(k)分布呈现d波对称性Qxy主导。自旋动量锁定的对称性特征N||ysz分量呈现sd波对称性。s波分量Q0对称性来自Mz但幅度很小d波分量Qv对称性来自M_beta_z是主导贡献。这种sd波混合是N||y构型的独特特征。N||xsz分量呈现纯d波对称性Qxy对称性源自Mxyz。sx和sy分量也各具有不同的宇称Qyz和Qzx。实验意义自旋动量锁定的不同宇称直接反映在磁SHE的各向异性上。通过测量磁SHE的非零分量可以反推序参数类型和自旋动量锁定的对称性。图 4(a) 耗散电导率sigma_xx^(J)对散射率eta的依赖性。(b,c) N||y和N||x构型的反常霍尔电导率sigma_ij^(E)随化学势E-EF的变化。AHE偶极vs八极的明确判据N||yAHC仅在xy平面非零sigma_xy^(E) -sigma_yx^(E)峰值约-227 S/cmE-EF-1.0 eV处。这直接对应磁偶极Mz作为序参数。在带隙内0-0.6 eVAHC为零。N||xAHC恒为零。这是因为磁八极Mxyz作为序参数不具有净磁化因此不产生AHE。这是区分N||y和N||x的最直接实验判据。物理本质AHE与净磁化共享相同的对称性要求。在交错磁体中AHE的存在与否直接反映了磁偶极是否作为序参数。图 5磁自旋霍尔电导率sigma_ij^(s(J),s)。(a-c) N||y构型sx,y,z。(d-f) N||x构型sx,y,z。磁SHE各向异性揭示序参数秩N||y磁SHE具有丰富的非零分量。sigma_yz^(s(J),y)在E-EF-0.8 eV处达到-1.4x10^3 S/cm。sigma_xx^(s(J),z)和sigma_yy^(s(J),z)也非零但各向异性显著。N||x磁SHE非零分量集不同。sigma_yx^(s(J),z)在E-EF-1.0 eV处达到1.5x10^3 S/cm。与N||y相比N||x的磁SHE分量与内禀SHE分量重叠因为T-even和T-odd多极矩无法区分。关键趋势自旋极化沿Neel矢量方向的分量更大的磁SHE因为费米面附近自旋极化沿Neel矢量方向更强。空穴掺杂比电子掺杂产生更大的响应。图 6磁自旋霍尔角beta_ij^(s(J),s)。(a) N||y构型最大SHA16%。(b) N||x构型最大SHA8.2%。16%自旋霍尔角超越Pt的效率N||y的最大SHA16%E-EF-0.04 eV是Pt5-10%的2倍以上与beta-Ta12-15%和AuW合金10%相当或更优。虽然不及beta-W的30%但MnTe的优势在于零净磁矩和丰富的可调控性。N||x的最大SHA8.2%E-EF-0.09 eV与Pt相当。两种构型的SHA差异再次反映了序参数的不同秩对自旋输运的影响。自旋极化方向依赖性沿Neel矢量方向的SHA分量更大。N||y中y极化分量的SHA大于x极化分量N||x中x极化分量的SHA大于y极化分量。图 7耗散电导率sigma_ij^(J)。(a) N||y构型。(b) N||x构型。两者相似因为电导率是T-even量不受磁序参数类型影响。图 8内禀自旋霍尔电导率sigma_ij^(s(E),s)。(a-c) N||y构型。(d-f) N||x构型。两种构型的内禀SHE相同因为除去时间反演操作外两者对称性相同。DFT Tips【DFT Tip 1】PAOFLOW紧束缚投影从VASP到输运计算的关键步骤PAOFLOW是连接DFT电子结构和输运性质计算的桥梁。它将VASP波函数投影到赝原子轨道PAO基组上构建实空间紧束缚哈密顿量然后通过Fourier插值获得任意稠密k点的能带和波函数。关键参数(1) PAO基组必须足够大以准确再现DFT能带。MnTe使用了Mn(10轨道)和Te(14轨道)这远大于最小基组确保带隙和能带色散的准确性。(2) 插值k点密度200x200x124需要足够高以保证输运系数的收敛。常见陷阱PAO基组不够大会导致紧束缚能带与DFT能带偏差大特别是在高能区域。建议在PAOFLOW中先对比投影能带与DFT能带确保关键能量窗口-2到2 eV内误差小于10 meV。【DFT Tip 2】Kubo公式中散射率eta的选择eta是Kubo公式中最重要的参数之一。它与实验弛豫时间tau的关系为etahbar/tau。本文中eta0.01 eV对应tau约66 fs是通过匹配实验电导率0.2-20 S/cm确定的。eta依赖性的物理(1) 耗散电导率sigma_xx^(J)在小eta下与1/eta成正比金属行为在大eta下偏离线性(2) MnTe在EF处表现为绝缘体行为sigma随eta增大而增大因为EF位于价带顶小eta时贡献态少。常见错误盲目使用eta0.01 eV而不验证。不同材料的eta应通过实验电导率标定。对于高电阻率材料eta可能需要更小对于金属eta可以更大。【DFT Tip 3】MnTe的DFTU设置U值的影响与验证本文使用U4.0 eV, J0.97 eVDudarev方案UeffU-J3.03 eV。该值在MnTe的常用范围内Ueff3-5 eV与GW计算和ARPES实验结果一致。U值的影响(1) 增大U会导致Mn-3d能级更局域带隙增大磁矩增大(2) 自旋劈裂幅度随U增大而增大(3) 输运系数AHE、SHE的定量值依赖于U但定性趋势如N||y vs N||x的差异对U值不敏感。建议对不同U值Ueff2,3,4,5 eV进行系统性测试验证输运系数的定性结论是否稳健。如果U值的改变导致序参数类型判断反转则需要更谨慎地选择U值。【DFT Tip 4】SOC自洽计算在交错磁体中的重要性本文所有计算都自洽包含SOC而非后处理添加。这是因为SOC不仅影响能带结构自旋劈裂、带隙还影响自旋极化方向和磁矩的大小。非自洽SOC的局限在非自洽SOC计算中电荷密度来自非SOC自洽计算SOC对电荷密度的反馈被忽略。对于MnTe这样的体系SOC对Mn-3d轨道的占据数和磁矩大小有不可忽略的影响。建议对于包含重元素Te的体系以及研究自旋相关输运性质时始终使用自洽SOC。对于轻元素体系非自洽SOC可能足够。【DFT Tip 5】磁构型验证MAGMOM与对称性一致性在DFT中设置交错磁序需要手动指定MAGMOM标签。对于MnTe N||yMAGMOM 0 4 0 0 -4 0每个Mn原子指定磁矩方向和大小。对于N||xMAGMOM 4 0 0 -4 0 0。验证步骤(1) 检查自洽后OSZICAR中的磁矩是否与预期一致方向和大小(2) 检查自旋极化能带是否展示预期的自旋劈裂模式(3) 验证磁矩方向在自洽过程中保持稳定不旋转。常见错误MAGMOM初始值设置过小可能导致收敛到非磁态或错误磁态。对于MnTe建议初始磁矩设置为4-5 muB/Mn大于自洽收敛值以确保收敛到正确磁态。【DFT Tip 6】线性响应计算中的T-even和T-odd分离Kubo公式中的响应函数可以根据时间反演对称性分为T-even和T-odd两类。T-even响应耗散电导率、内禀SHE在时间反演下不变T-odd响应AHE、磁SHE在时间反演下变号。在交错磁体中T-even和T-odd混合是一个关键特征。对于N||x构型T-even和T-odd多极矩无法区分导致磁SHE和内禀SHE在相同张量分量中出现。这是交错磁体区别于普通反铁磁体的重要特征。建议在分析输运数据时始终进行T-even和T-odd分解以理解不同物理机制的贡献。PAOFLOW中可以通过设置不同的算符组合来实现这一分解。【DFT Tip 7】磁自旋霍尔角SHA的实验验证SHA是连接理论计算和实验测量的关键量。理论上SHA (e/hbar) * sigma_sHE / sigma_xx其中sigma_sHE和sigma_xx都是能量的函数。实验中的SHA通过自旋扭矩铁磁共振ST-FMR或自旋泵浦测量。理论SHA与实验比较时需注意(1) 理论SHA是零温单能量值实验SHA是有限温度的能量积分(2) 理论SHA对散射率eta敏感eta的标定直接影响SHA值。建议在论文中同时报告SHA的能量依赖曲线而非仅报告最大值并讨论eta不确定性对SHA的影响。【DFT Tip 8】多极矩分析的实用工具与流程多极矩分析是理解交错磁体输运性质的有力工具。实用流程(1) 确定磁点群/磁空间群使用Bilbao Crystallographic Server或FINDSYM(2) 找出属于全对称表示的多极矩使用MULTIPOLAR或自写脚本(3) 根据多极矩-响应张量对应表推导允许的非零响应分量。关键参考Yatsushiro et al., PRB 104, 054412 (2021) - 122个磁点群的多极矩分类Hayami et al., JPSJ 93, 072001 (2024) - 多极矩统一描述综述。注意多极矩分析给出的是对称性允许的响应分量不给出数值大小。数值大小需要DFT线性响应计算。两者结合才能给出完整的物理图像。知识扩展【知识扩展 1】Kubo线性响应理论从微观到宏观【理论解释】Kubo公式是量子统计力学中计算线性响应的标准框架。核心思想外场如电场微扰下物理量如电流的期望值变化与未微扰系统的关联函数成正比。Kubo公式的第一性原理实现需要DFT波函数和能带作为输入。【关键近似】恒定弛豫时间近似CRTA假设所有态的弛豫时间相同用一个常数etahbar/tau代替复杂的散射机制。CRTA的合理性取决于散射机制的细节在杂质散射主导的低温下CRTA适用在电子-声子散射主导的高温下需要更复杂的处理。【方法比较】Kubo公式 vs Boltzmann输运方程Kubo公式自动包含带间跃迁interband和带内跃迁intraband贡献而Boltzmann方程通常只处理带内输运。对于拓扑材料如交错磁体带间贡献Berry曲率相关至关重要因此Kubo公式更合适。【经典参考】Kubo, JPSJ 12, 570 (1957) - 原始论文Freimuth et al., PRB 90, 174423 (2014) - 第一性原理Kubo公式的现代实现。【知识扩展 2】自旋霍尔效应从内禀到磁性【理论解释】自旋霍尔效应SHE是指电荷流在SOC作用下产生横向自旋流的现象。内禀SHE由SOC和能带Berry曲率驱动是T-even的时间反演下不变。磁SHE由时间反演对称性破缺和交换相互作用驱动是T-odd的。【物理机制对比】内禀SHESOC导致自旋依赖的横向偏转类似AHE但自旋流替代电荷流。磁SHE费米面在外电场下发生位移由于自旋依赖的费米速度产生净自旋流。磁SHE不需要SOC理论上可以存在于无SOC的体系。【材料对比】传统SHE材料PtSHA~5-10%、beta-W~30%、beta-Ta12-15%。交错磁体SHEMnTe16%、RuO2理论预测、CrSb理论预测。交错磁体的优势在于零净磁矩和丰富的对称性调控。【经典参考】Sinova et al., RMP 87, 1213 (2015) - SHE综述Gonzalez-Hernandez et al., PRL 126, 127701 (2021) - 交错磁体磁SHE。科研经验【科研经验 1】交错磁体研究中磁构型的确定问题如何确定交错磁体的磁基态和可能的亚稳态磁构型原因交错磁体的磁性由交换耦合、SOC和晶体场共同决定。不同的磁构型可能具有相近的能量DFT计算需要仔细收敛。解决方案(1) 从非磁态出发弛豫晶格常数(2) 根据对称性分析列出所有可能的共线磁构型使用ISODISTORT或人工枚举(3) 对每种磁构型进行DFTUSOC计算比较总能量(4) 验证能量最低构型的动力学稳定性声子谱无虚频。建议对于MnTeN||y是实验已知的磁基态能量最低。N||x是更高能量的亚稳态但可以通过应力-磁场组合诱导。在计算N||x构型时确保MAGMOM设置正确且自洽收敛到目标磁态。【科研经验 2】自旋霍尔电导率的收敛性检验问题Kubo公式计算的自旋霍尔电导率对k点密度和eta参数敏感如何确保收敛原因SHC涉及带间跃迁的矩阵元在能带简并或近简并处对k点采样和eta展宽特别敏感。收敛性不足会导致SHC的虚假振荡或错误符号。解决方案(1) 逐步增加k点密度100x100x62 - 200x200x124 - 300x300x186观察SHC是否收敛(2) 对不同的eta值0.005, 0.01, 0.02 eV计算SHC验证定性趋势是否一致(3) 检查SHC的能量依赖曲线是否平滑异常尖峰可能表明收敛不足。建议对于金属体系SHC的收敛需要比带结构更密的k点。至少需要200x200x124本文设置或更密。如果计算资源有限优先保证SHC峰值的收敛而不是全能量范围的收敛。如果是我我还会继续算【继续算 1】Wannier90紧束缚模型验证PAOFLOW结果为什么值得算PAOFLOW和Wannier90都是构建紧束缚模型的方法但投影基组不同PAO vs MLWF。交叉验证可以确认输运结果的可靠性排除方法依赖的伪影。能回答的问题PAOFLOW和Wannier90的SHC是否一致Mn-d和Te-p轨道对SHC的贡献有何不同适合体系所有需要输运计算的体系。输入Wannier90Mn-d, Te-p轨道。【继续算 2】Berry曲率分布AHE和SHE的k空间起源为什么值得算AHE和SHC是Berry曲率在Brillouin区的积分。计算k分辨的Berry曲率分布可以揭示哪些k点和哪些能带对输运贡献最大。能回答的问题AHE的k空间热点在哪里哪些能带交叉点贡献最大的Berry曲率N||y和N||x的Berry曲率分布有何不同适合体系所有拓扑材料。输入Wannier90 WannierTools/WannierBerri。【继续算 3】有限温度效应费米-狄拉克展宽为什么值得算本文计算的是零温响应。常温下费米-狄拉克分布展宽会平滑能量依赖的尖锐特征。实验SHA是温度依赖的需要理论预测温度效应。能回答的问题SHA在300 K时是否仍然显著温度展宽是否改变SHA的符号或峰值位置适合体系所有需要与实验比较的体系。输入在Kubo公式中引入费米-狄拉克展宽。【继续算 4】应力调控诱导N||x构型并调控磁SHE为什么值得算N||x构型尚未实验实现但理论预测可以通过剪切应力磁场组合诱导。DFT可以预测实现N||x所需的应力大小和方向。能回答的问题需要多大的应力才能稳定N||x构型应力如何影响磁SHE的大小和各向异性适合体系所有交错磁体。输入不同应力下的DFTUSOC计算。【继续算 5】其他交错磁体的多极矩分析从MnTe到RuO2、CrSb为什么值得算本文的多极矩框架是通用的可以迁移到其他交错磁体。对RuO2d波、CrSbg波、Mn5Si3等体系进行类似分析建立交错磁体输运的系统性认识。能回答的问题不同交错磁体的磁SHE各向异性有何规律是否存在SHA超过30%的交错磁体候选材料适合体系所有交错磁体候选材料。输入高通量DFT PAOFLOW但可并行。【继续算 6】自旋累积与非局域输运器件模拟为什么值得算磁SHE的最终目标是器件应用。将DFT计算的SHC作为输入通过漂移-扩散方程模拟MnTe/重金属异质结中的自旋累积和自旋扭矩。能回答的问题MnTe基自旋霍尔器件中的自旋扩散长度是多少实现磁化翻转所需的临界电流密度是多少适合体系所有自旋电子学材料。输入SHC张量自旋扩散方程。Hirakida et al., PRB 113, 155104 (2026) | 交错磁体 多极矩 磁自旋霍尔效应 Kubo公式 PAOFLOW PRB 交错磁体MnTe的多极矩分析与磁自旋霍尔效应
PRB 交错磁体MnTe的多极矩分析与磁自旋霍尔效应PRB 113, 155104 (2026)交错磁体MnTe的多极矩分析与磁自旋霍尔效应Multipole Analysis of Spin Currents in Altermagnetic MnTe导读交错磁体中不同的磁构型产生不同秩的序参数这直接决定了自旋霍尔效应的存在性和各向异性。本文通过多极矩框架结合Kubo公式第一性原理计算系统分析了alpha-MnTe中N||y和N||x两种磁构型的电荷输运、反常霍尔效应和磁自旋霍尔效应。核心发现N||y构型磁偶极Mz同时具有AHE和磁SHE磁自旋霍尔角高达16%超过Pt的2倍N||x构型磁八极Mxyz只有磁SHE而无AHE。磁SHE的各向异性为实验判别序参数类型提供了实用方法。一、前言背景交错磁体MnTe从d波自旋劈裂到磁多极矩交错磁体Altermagnet是2022年以来凝聚态物理最受关注的发现之一。与传统反铁磁体由平移或反演连接反平行自旋不同交错磁体中反平行自旋由旋转或旋转-平移联合操作连接导致动量空间中产生非相对论自旋劈裂。alpha-MnTe是典型的g波交错磁体Neel温度~310 K具有六方NiAs型结构。ARPES实验已明确观测到其Kramers自旋简并的解除。本文考虑两种Neel矢量取向N||y实验可实现的磁构型和N||x尚未实现但理论预测的磁构型。本文的核心创新在于将多极矩multipole框架引入交错磁体的自旋输运分析建立了磁序参数与自旋霍尔效应SHE之间的直接对应关系并通过PAOFLOW Kubo公式进行了第一性原理验证。多极矩框架从对称性到输运响应的桥梁核心问题不同的磁构型N||y vs N||x如何影响自旋电流的产生磁自旋霍尔效应的各向异性如何反映序参数的类型方法体系VASPDFTUSOC基态电子结构 - PAOFLOW紧束缚投影构建Wannier-like Hamiltonian - Kubo公式线性响应电导率/自旋霍尔电导率计算 - 多极矩对称性分析。关键发现(1) N||y构型具有磁偶极矩Mz作为序参数同时出现AHE和磁SHE(2) N||x构型具有磁八极矩Mxyz作为序参数只有磁SHE而无AHE(3) 磁自旋霍尔角SHA高达16%超过Pt5-10%与beta-Ta12-15%相当。MnTe磁自旋霍尔效应研究流程。VASP DFTUSOC电子结构 - PAOFLOW紧束缚投影 - Kubo公式线性响应电导率、AHE、磁SHE、内禀SHE - 多极矩对称性分析。核心发现N||y构型同时具有AHE和磁SHESHA16%N||x构型只有磁SHESHA8.2%磁SHE各向异性可用于判别序参数类型。二、研究方法VASP DFTUSOC基态电子结构计算计算设置PBE泛函 PAW赝势, ENCUT350 eV, Gamma-centered 20x20x12 k点网格, 实验晶格参数a4.134 A, c6.652 A。SOC在所有计算中自洽包含。Hubbard U处理Mn-3d轨道采用DFTUDudarev方案U4.0 eV, J0.97 eV。该值在MnTe的常用范围内与GW计算和ARPES实验结果一致。注意U值的选择对MnTe的价带结构有显著影响但定性趋势对U值不敏感。两种磁构型(1) N||y磁点群mmm磁空间群Cmcm实验可实现(2) N||x磁点群mmm磁空间群Cmcm尚未实验实现但可通过应力-磁场组合诱导。PAOFLOW Kubo从DFT到输运响应PAOFLOW方法将VASP波函数投影到赝原子轨道PAO基组上构建紧束缚哈密顿量。Mn基组[1s,2s,2p,3s,3p,3d,4s,4p,4d,5s]Te基组[1s,2s,2p,3s,3p,3d,4s,4p,4d,5s,5p,5d,6s,6p]。通过Fourier插值零填充将哈密顿量插值到200x200x124的稠密k点网格。Kubo公式在恒定弛豫时间近似下计算六种线性响应耗散电导率T-even、AHET-odd、耗散REET-even、磁REET-odd、内禀SHET-even、磁SHET-odd。散射率eta0.01 eV由实验电导率确定。关键理解磁SHE是T-odd的耗散依赖响应由时间反演对称性破缺和交换相互作用驱动其物理起源与内禀SHET-evenSOC驱动完全不同。Kubo公式耗散贡献eta为散射率A和B为与物理量和外场对应的算符Kubo公式非耗散贡献在eta-0极限下独立于eta对应电场驱动的响应磁自旋霍尔角SHA无量纲量表征电荷-自旋转换效率beta_s(J),s_ij磁八极矩M_xyz最低阶交错磁序参数由电四极矩Q与磁偶极矩M_z的耦合产生三、核心结果图 1alpha-MnTe的晶体和磁结构。(a) 晶体结构NiAs型六方。(b) N||y磁构型磁矩沿y/-y方向交替排列。(c) N||x磁构型磁矩沿x/-x方向交替排列。图中标注了不可约对称操作。磁构型与多极矩序参数N||y构型磁点群mmm磁偶极矩Mz为最低阶序参数同时存在磁八极矩M_alpha_z和M_beta_z。这意味着该构型可以共享铁磁体沿z轴磁化的对称性无需额外对称性破缺。N||x构型磁点群mmm无磁偶极矩最低阶交错磁序参数为磁八极矩Mxyz。该构型中T-even和T-odd多极矩无法区分因为时间反演不是对称操作。多极矩分析的意义不同秩的序参数直接决定了响应张量的非零元素和对称性形式。例如磁偶极Mz允许AHEsigma_xy而磁八极Mxyz禁止AHE但允许特定分量的磁SHE。图 2自旋极化能带结构。(a) N||y构型的自旋极化能带颜色表示sz分量。(b) N||x构型的自旋极化能带。(c) 第一布里渊区和对称点。图 3自旋-动量锁定。(a) N||y构型在k_x-k_y平面上的sz(k)分布呈现sd波对称性Q0Qv主导。(b) N||x构型的sz(k)分布呈现d波对称性Qxy主导。自旋动量锁定的对称性特征N||ysz分量呈现sd波对称性。s波分量Q0对称性来自Mz但幅度很小d波分量Qv对称性来自M_beta_z是主导贡献。这种sd波混合是N||y构型的独特特征。N||xsz分量呈现纯d波对称性Qxy对称性源自Mxyz。sx和sy分量也各具有不同的宇称Qyz和Qzx。实验意义自旋动量锁定的不同宇称直接反映在磁SHE的各向异性上。通过测量磁SHE的非零分量可以反推序参数类型和自旋动量锁定的对称性。图 4(a) 耗散电导率sigma_xx^(J)对散射率eta的依赖性。(b,c) N||y和N||x构型的反常霍尔电导率sigma_ij^(E)随化学势E-EF的变化。AHE偶极vs八极的明确判据N||yAHC仅在xy平面非零sigma_xy^(E) -sigma_yx^(E)峰值约-227 S/cmE-EF-1.0 eV处。这直接对应磁偶极Mz作为序参数。在带隙内0-0.6 eVAHC为零。N||xAHC恒为零。这是因为磁八极Mxyz作为序参数不具有净磁化因此不产生AHE。这是区分N||y和N||x的最直接实验判据。物理本质AHE与净磁化共享相同的对称性要求。在交错磁体中AHE的存在与否直接反映了磁偶极是否作为序参数。图 5磁自旋霍尔电导率sigma_ij^(s(J),s)。(a-c) N||y构型sx,y,z。(d-f) N||x构型sx,y,z。磁SHE各向异性揭示序参数秩N||y磁SHE具有丰富的非零分量。sigma_yz^(s(J),y)在E-EF-0.8 eV处达到-1.4x10^3 S/cm。sigma_xx^(s(J),z)和sigma_yy^(s(J),z)也非零但各向异性显著。N||x磁SHE非零分量集不同。sigma_yx^(s(J),z)在E-EF-1.0 eV处达到1.5x10^3 S/cm。与N||y相比N||x的磁SHE分量与内禀SHE分量重叠因为T-even和T-odd多极矩无法区分。关键趋势自旋极化沿Neel矢量方向的分量更大的磁SHE因为费米面附近自旋极化沿Neel矢量方向更强。空穴掺杂比电子掺杂产生更大的响应。图 6磁自旋霍尔角beta_ij^(s(J),s)。(a) N||y构型最大SHA16%。(b) N||x构型最大SHA8.2%。16%自旋霍尔角超越Pt的效率N||y的最大SHA16%E-EF-0.04 eV是Pt5-10%的2倍以上与beta-Ta12-15%和AuW合金10%相当或更优。虽然不及beta-W的30%但MnTe的优势在于零净磁矩和丰富的可调控性。N||x的最大SHA8.2%E-EF-0.09 eV与Pt相当。两种构型的SHA差异再次反映了序参数的不同秩对自旋输运的影响。自旋极化方向依赖性沿Neel矢量方向的SHA分量更大。N||y中y极化分量的SHA大于x极化分量N||x中x极化分量的SHA大于y极化分量。图 7耗散电导率sigma_ij^(J)。(a) N||y构型。(b) N||x构型。两者相似因为电导率是T-even量不受磁序参数类型影响。图 8内禀自旋霍尔电导率sigma_ij^(s(E),s)。(a-c) N||y构型。(d-f) N||x构型。两种构型的内禀SHE相同因为除去时间反演操作外两者对称性相同。DFT Tips【DFT Tip 1】PAOFLOW紧束缚投影从VASP到输运计算的关键步骤PAOFLOW是连接DFT电子结构和输运性质计算的桥梁。它将VASP波函数投影到赝原子轨道PAO基组上构建实空间紧束缚哈密顿量然后通过Fourier插值获得任意稠密k点的能带和波函数。关键参数(1) PAO基组必须足够大以准确再现DFT能带。MnTe使用了Mn(10轨道)和Te(14轨道)这远大于最小基组确保带隙和能带色散的准确性。(2) 插值k点密度200x200x124需要足够高以保证输运系数的收敛。常见陷阱PAO基组不够大会导致紧束缚能带与DFT能带偏差大特别是在高能区域。建议在PAOFLOW中先对比投影能带与DFT能带确保关键能量窗口-2到2 eV内误差小于10 meV。【DFT Tip 2】Kubo公式中散射率eta的选择eta是Kubo公式中最重要的参数之一。它与实验弛豫时间tau的关系为etahbar/tau。本文中eta0.01 eV对应tau约66 fs是通过匹配实验电导率0.2-20 S/cm确定的。eta依赖性的物理(1) 耗散电导率sigma_xx^(J)在小eta下与1/eta成正比金属行为在大eta下偏离线性(2) MnTe在EF处表现为绝缘体行为sigma随eta增大而增大因为EF位于价带顶小eta时贡献态少。常见错误盲目使用eta0.01 eV而不验证。不同材料的eta应通过实验电导率标定。对于高电阻率材料eta可能需要更小对于金属eta可以更大。【DFT Tip 3】MnTe的DFTU设置U值的影响与验证本文使用U4.0 eV, J0.97 eVDudarev方案UeffU-J3.03 eV。该值在MnTe的常用范围内Ueff3-5 eV与GW计算和ARPES实验结果一致。U值的影响(1) 增大U会导致Mn-3d能级更局域带隙增大磁矩增大(2) 自旋劈裂幅度随U增大而增大(3) 输运系数AHE、SHE的定量值依赖于U但定性趋势如N||y vs N||x的差异对U值不敏感。建议对不同U值Ueff2,3,4,5 eV进行系统性测试验证输运系数的定性结论是否稳健。如果U值的改变导致序参数类型判断反转则需要更谨慎地选择U值。【DFT Tip 4】SOC自洽计算在交错磁体中的重要性本文所有计算都自洽包含SOC而非后处理添加。这是因为SOC不仅影响能带结构自旋劈裂、带隙还影响自旋极化方向和磁矩的大小。非自洽SOC的局限在非自洽SOC计算中电荷密度来自非SOC自洽计算SOC对电荷密度的反馈被忽略。对于MnTe这样的体系SOC对Mn-3d轨道的占据数和磁矩大小有不可忽略的影响。建议对于包含重元素Te的体系以及研究自旋相关输运性质时始终使用自洽SOC。对于轻元素体系非自洽SOC可能足够。【DFT Tip 5】磁构型验证MAGMOM与对称性一致性在DFT中设置交错磁序需要手动指定MAGMOM标签。对于MnTe N||yMAGMOM 0 4 0 0 -4 0每个Mn原子指定磁矩方向和大小。对于N||xMAGMOM 4 0 0 -4 0 0。验证步骤(1) 检查自洽后OSZICAR中的磁矩是否与预期一致方向和大小(2) 检查自旋极化能带是否展示预期的自旋劈裂模式(3) 验证磁矩方向在自洽过程中保持稳定不旋转。常见错误MAGMOM初始值设置过小可能导致收敛到非磁态或错误磁态。对于MnTe建议初始磁矩设置为4-5 muB/Mn大于自洽收敛值以确保收敛到正确磁态。【DFT Tip 6】线性响应计算中的T-even和T-odd分离Kubo公式中的响应函数可以根据时间反演对称性分为T-even和T-odd两类。T-even响应耗散电导率、内禀SHE在时间反演下不变T-odd响应AHE、磁SHE在时间反演下变号。在交错磁体中T-even和T-odd混合是一个关键特征。对于N||x构型T-even和T-odd多极矩无法区分导致磁SHE和内禀SHE在相同张量分量中出现。这是交错磁体区别于普通反铁磁体的重要特征。建议在分析输运数据时始终进行T-even和T-odd分解以理解不同物理机制的贡献。PAOFLOW中可以通过设置不同的算符组合来实现这一分解。【DFT Tip 7】磁自旋霍尔角SHA的实验验证SHA是连接理论计算和实验测量的关键量。理论上SHA (e/hbar) * sigma_sHE / sigma_xx其中sigma_sHE和sigma_xx都是能量的函数。实验中的SHA通过自旋扭矩铁磁共振ST-FMR或自旋泵浦测量。理论SHA与实验比较时需注意(1) 理论SHA是零温单能量值实验SHA是有限温度的能量积分(2) 理论SHA对散射率eta敏感eta的标定直接影响SHA值。建议在论文中同时报告SHA的能量依赖曲线而非仅报告最大值并讨论eta不确定性对SHA的影响。【DFT Tip 8】多极矩分析的实用工具与流程多极矩分析是理解交错磁体输运性质的有力工具。实用流程(1) 确定磁点群/磁空间群使用Bilbao Crystallographic Server或FINDSYM(2) 找出属于全对称表示的多极矩使用MULTIPOLAR或自写脚本(3) 根据多极矩-响应张量对应表推导允许的非零响应分量。关键参考Yatsushiro et al., PRB 104, 054412 (2021) - 122个磁点群的多极矩分类Hayami et al., JPSJ 93, 072001 (2024) - 多极矩统一描述综述。注意多极矩分析给出的是对称性允许的响应分量不给出数值大小。数值大小需要DFT线性响应计算。两者结合才能给出完整的物理图像。知识扩展【知识扩展 1】Kubo线性响应理论从微观到宏观【理论解释】Kubo公式是量子统计力学中计算线性响应的标准框架。核心思想外场如电场微扰下物理量如电流的期望值变化与未微扰系统的关联函数成正比。Kubo公式的第一性原理实现需要DFT波函数和能带作为输入。【关键近似】恒定弛豫时间近似CRTA假设所有态的弛豫时间相同用一个常数etahbar/tau代替复杂的散射机制。CRTA的合理性取决于散射机制的细节在杂质散射主导的低温下CRTA适用在电子-声子散射主导的高温下需要更复杂的处理。【方法比较】Kubo公式 vs Boltzmann输运方程Kubo公式自动包含带间跃迁interband和带内跃迁intraband贡献而Boltzmann方程通常只处理带内输运。对于拓扑材料如交错磁体带间贡献Berry曲率相关至关重要因此Kubo公式更合适。【经典参考】Kubo, JPSJ 12, 570 (1957) - 原始论文Freimuth et al., PRB 90, 174423 (2014) - 第一性原理Kubo公式的现代实现。【知识扩展 2】自旋霍尔效应从内禀到磁性【理论解释】自旋霍尔效应SHE是指电荷流在SOC作用下产生横向自旋流的现象。内禀SHE由SOC和能带Berry曲率驱动是T-even的时间反演下不变。磁SHE由时间反演对称性破缺和交换相互作用驱动是T-odd的。【物理机制对比】内禀SHESOC导致自旋依赖的横向偏转类似AHE但自旋流替代电荷流。磁SHE费米面在外电场下发生位移由于自旋依赖的费米速度产生净自旋流。磁SHE不需要SOC理论上可以存在于无SOC的体系。【材料对比】传统SHE材料PtSHA~5-10%、beta-W~30%、beta-Ta12-15%。交错磁体SHEMnTe16%、RuO2理论预测、CrSb理论预测。交错磁体的优势在于零净磁矩和丰富的对称性调控。【经典参考】Sinova et al., RMP 87, 1213 (2015) - SHE综述Gonzalez-Hernandez et al., PRL 126, 127701 (2021) - 交错磁体磁SHE。科研经验【科研经验 1】交错磁体研究中磁构型的确定问题如何确定交错磁体的磁基态和可能的亚稳态磁构型原因交错磁体的磁性由交换耦合、SOC和晶体场共同决定。不同的磁构型可能具有相近的能量DFT计算需要仔细收敛。解决方案(1) 从非磁态出发弛豫晶格常数(2) 根据对称性分析列出所有可能的共线磁构型使用ISODISTORT或人工枚举(3) 对每种磁构型进行DFTUSOC计算比较总能量(4) 验证能量最低构型的动力学稳定性声子谱无虚频。建议对于MnTeN||y是实验已知的磁基态能量最低。N||x是更高能量的亚稳态但可以通过应力-磁场组合诱导。在计算N||x构型时确保MAGMOM设置正确且自洽收敛到目标磁态。【科研经验 2】自旋霍尔电导率的收敛性检验问题Kubo公式计算的自旋霍尔电导率对k点密度和eta参数敏感如何确保收敛原因SHC涉及带间跃迁的矩阵元在能带简并或近简并处对k点采样和eta展宽特别敏感。收敛性不足会导致SHC的虚假振荡或错误符号。解决方案(1) 逐步增加k点密度100x100x62 - 200x200x124 - 300x300x186观察SHC是否收敛(2) 对不同的eta值0.005, 0.01, 0.02 eV计算SHC验证定性趋势是否一致(3) 检查SHC的能量依赖曲线是否平滑异常尖峰可能表明收敛不足。建议对于金属体系SHC的收敛需要比带结构更密的k点。至少需要200x200x124本文设置或更密。如果计算资源有限优先保证SHC峰值的收敛而不是全能量范围的收敛。如果是我我还会继续算【继续算 1】Wannier90紧束缚模型验证PAOFLOW结果为什么值得算PAOFLOW和Wannier90都是构建紧束缚模型的方法但投影基组不同PAO vs MLWF。交叉验证可以确认输运结果的可靠性排除方法依赖的伪影。能回答的问题PAOFLOW和Wannier90的SHC是否一致Mn-d和Te-p轨道对SHC的贡献有何不同适合体系所有需要输运计算的体系。输入Wannier90Mn-d, Te-p轨道。【继续算 2】Berry曲率分布AHE和SHE的k空间起源为什么值得算AHE和SHC是Berry曲率在Brillouin区的积分。计算k分辨的Berry曲率分布可以揭示哪些k点和哪些能带对输运贡献最大。能回答的问题AHE的k空间热点在哪里哪些能带交叉点贡献最大的Berry曲率N||y和N||x的Berry曲率分布有何不同适合体系所有拓扑材料。输入Wannier90 WannierTools/WannierBerri。【继续算 3】有限温度效应费米-狄拉克展宽为什么值得算本文计算的是零温响应。常温下费米-狄拉克分布展宽会平滑能量依赖的尖锐特征。实验SHA是温度依赖的需要理论预测温度效应。能回答的问题SHA在300 K时是否仍然显著温度展宽是否改变SHA的符号或峰值位置适合体系所有需要与实验比较的体系。输入在Kubo公式中引入费米-狄拉克展宽。【继续算 4】应力调控诱导N||x构型并调控磁SHE为什么值得算N||x构型尚未实验实现但理论预测可以通过剪切应力磁场组合诱导。DFT可以预测实现N||x所需的应力大小和方向。能回答的问题需要多大的应力才能稳定N||x构型应力如何影响磁SHE的大小和各向异性适合体系所有交错磁体。输入不同应力下的DFTUSOC计算。【继续算 5】其他交错磁体的多极矩分析从MnTe到RuO2、CrSb为什么值得算本文的多极矩框架是通用的可以迁移到其他交错磁体。对RuO2d波、CrSbg波、Mn5Si3等体系进行类似分析建立交错磁体输运的系统性认识。能回答的问题不同交错磁体的磁SHE各向异性有何规律是否存在SHA超过30%的交错磁体候选材料适合体系所有交错磁体候选材料。输入高通量DFT PAOFLOW但可并行。【继续算 6】自旋累积与非局域输运器件模拟为什么值得算磁SHE的最终目标是器件应用。将DFT计算的SHC作为输入通过漂移-扩散方程模拟MnTe/重金属异质结中的自旋累积和自旋扭矩。能回答的问题MnTe基自旋霍尔器件中的自旋扩散长度是多少实现磁化翻转所需的临界电流密度是多少适合体系所有自旋电子学材料。输入SHC张量自旋扩散方程。Hirakida et al., PRB 113, 155104 (2026) | 交错磁体 多极矩 磁自旋霍尔效应 Kubo公式 PAOFLOW