← 返回 PaperDaily 大模型与智能体

磁体第三家族迎实锤:g波多极首次“现形”

如果只说“铁磁性”和“反铁磁性”,那对磁有序世界的理解可是缺了一大块。铁磁体像一群“死磕同一个方向”的小磁针,一上外场就齐齐转向,磁化强度很大;反铁磁体则是“相邻小磁针反着排”,宏观上互相抵消,净磁化近似为零。

磁体第三家族迎实锤:g波多极首次“现形”
原论文信息如下:
论文标题:
Observation of g-wave altermagnetic multipole
发表日期:
2026年09月
发表单位:
The University of Tokyo, Hokkaido University, Tohoku University, Institut Laue-Langevin, Gdańsk University of Technology, Karlsruhe Institute of Technology, Nanyang Technological University, RIKEN Center for Emergent Matter Science
原文链接:
https://arxiv.org/pdf/2609.01969v1.pdf

如果只说“铁磁性”和“反铁磁性”,那对磁有序世界的理解可是缺了一大块。铁磁体像一群“死磕同一个方向”的小磁针,一上外场就齐齐转向,磁化强度很大;反铁磁体则是“相邻小磁针反着排”,宏观上互相抵消,净磁化近似为零。然而近年来凝聚态物理里蹿红的交替磁体(altermagnet),却把人们的直觉往奇怪的方向拉了一把:它明明也是共线的反平行排列、净磁化几乎为零,却能在能带中劈裂出可观的、与自旋方向强相关的电子态,甚至劈裂能量能到电子伏特量级。这种既“反铁磁”又“铁磁般好用”的特性,让它一下成了自旋电子学圈子的新宠。
问题来了:一个没有净磁化的体系,凭什么打破时间反演对称性?传统教科书会说“反铁磁体时间反演后经过一个晶格平移就能回到原样”,但交替磁体做不到。真正的“幕后黑手”,并不是宏观磁矩,而是原子周围那种非常规的多极矩(multipole)——可以粗略理解成电子密度围绕原子核形成的各种非球形“花纹”。铁磁体的自旋和电荷分布紧密耦合成各向同性的s波(l=0),而交替磁体则可能携带d波(l=2)、g波(l=4)甚至i波(l=6)形状的花样。这些高阶花样平时藏在价电子云里,看不见摸不着,所以也被称为“隐藏序”(hidden order)。
东京大学领衔的团队这次干了一件相当硬核的事。他们把目光锁定在CrSb(铬锑)这种室温交替磁体上,利用高能同步辐射X射线衍射和极化中子衍射,第一次在实空间里直接“看见”并定量拆解了g波交替磁体的多极矩。底下这张图把“抽象概念”翻译成了比较直观的画面:g波多极以yz(3x²-y²)这种四下伸展、带着凹凸的瓣状各向异性,在晶格中一层正一层负地交替排列,自旋则在上面“见缝插针”。时间反演让自旋翻转,却不让多极形状跟着翻转,于是两条原本应该成对出现的电子态就不再等价,自旋劈裂就这么诞生了。
图1:g波交替磁体的多极矩描述。(a)沿z轴交替排列的g波多极与自旋,自旋在时间反演下翻转,多极不翻转,从而打破时间反演与平移组合对称性。(b)g波交替磁体的典型六角晶体结构。(c)动量空间中自旋劈裂的费米面示意图,其各向异性与实空间g波多极一致。
图1:g波交替磁体的多极矩描述。(a)g波多极与自旋沿z轴交替排列的示意;时间反演让自旋反转但多极花纹不变,于是两个自旋配对态不再等价,组合平移对称性也被打破。(b)g波交替磁体的典型六角晶体结构示意图。(c)动量空间中费米面的自旋劈裂,其各向异性与实空间g波多极形状一致。
换句话说,这篇论文不只是拿一台显微镜拍了个“亮晶晶”的照片,而是把“隐藏序”从数学预言推到了可测、可量、可证的位置。它解决的是整个交替磁体领域最关键也最悬空的一个问题:自旋劈裂的对称性来源到底长什么样、有多大?下面咱们一层层拆开看。

交替磁体中的“隐藏秩序”:多极矩的直接可视化

先解释一个容易绕晕的词:多极矩(multipole)是什么意思。
一个原子孤立在真空中时,周围电子云的角分布是球对称的,像一颗均匀的毛线球。可一旦放进晶体,周围配位原子形成的电场就会把简并的轨道能量劈开,电子云也跟着变得不是球了。对任意角分布做“拆解”——就像把一组复杂波形拆成不同阶次的正弦波——可以得到一整套基函数,也就是球谐函数Y_lm。l=0是各向同性的单极项,俗称s波;l=1是像哑铃一样的偶极项,叫p波;l=2是四瓣状的四极项,叫d波;l=3称八极,叫f波;l=4是更高阶的十六极(hexadecapole),光谱学记号正好轮到g,这就是“g波”的由来。
这类高阶项在f电子体系里被研究了很久,经常作为“隐藏电子相”的序参量,但一直很难直接观测。交替磁体理论则指出,每个Cr位点上的电子云可以携带Q43∝yz(3x²-y²)这样的g波多极;两个Cr子晶格上的Q43符号相反,形成反铁性(antiferroic)排列,也就是“一个凸、一个凹、互相顶着”。传统反铁磁体在时间反演加半平移下能完整映射回自身,而这里多极花样不跟着自旋翻转,导致整体对称性被实实在在地破坏。这样看来,探测g波多极就等于给交替磁体的“隐藏秩序”拍一张实空间证件照。
选CrSb当“模特”也很讲究:它是被反复验证过的室温g波交替磁体,早已有角分辨光电子能谱实验发现其动量空间存在巨大的g波自旋劈裂。但此前的实验都停留在“倒空间”的能带照片,从没有人把实空间里那个本应存在的多极直接抠出来。这篇论文最大的看点,就是完成了这一跳。

实验突破:如何“看见”价电子密度的各向异性

想要直接看多极,就需要测量原子周围的价电子密度(valence electron density,简称VED),即参与成键和磁性的外层电子在实空间里的分布。X射线衍射在这里是“王牌探针”:X射线的散射强度正比于电子密度,理论上只要测到足够多衍射点,再用傅里叶合成就能还原出晶体里的三维电子密度,精度甚至可以到原子尺度以下。
麻烦的是,X射线看到的是全部电子,其中贡献最大的通常是紧贴原子核的内层芯电子,而真正携带各向异性信息的价电子只占很小部分。团队采用的方法是芯电子差分傅里叶合成(core differential Fourier synthesis,简称CDFS):先把总结构因子里的芯电子贡献精确扣除,再做逆傅里叶变换,从而只留下价电子的密度分布。这个方法的精妙之处在于不预设价态模型,尽量让实验数据自己“说话”。
图2:同步辐射X射线衍射实现价电子密度(VED)实空间可视化。(a)CrSb的六角P63/mmc晶体结构。(b)同步辐射X射线衍射实验示意。(c)CDFS方法示意:扣除芯电子贡献后重建价电子密度。(d)100 K下高质量结构精修结果,R因子低于1%。
图2:同步辐射X射线衍射实现价电子密度实空间可视化。(a)CrSb的六角P6₃/mmc晶体结构,Cr形成两类三角对称子晶格。(b)高能同步辐射X射线衍射实验示意。(c)CDFS方法示意:扣除芯电子密度后从总结构因子中重建价电子密度。(d)100 K下衍射精修结果,观测与计算强度的一致性因子R低于1%。
从数学上看,CDFS重建价电子密度的过程可以写成下面的形式:
公式:价电子密度重建公式,即对扣除芯电子贡献后的结构因子做逆傅里叶变换。
其中F_obs是实验观测的总结构因子,F_calc^core是根据芯电子波函数算出的芯电子结构因子,两者相减后取逆傅里叶变换,再加上电子总数归一化项n/V,就得到实空间的价电子密度ρ(r)。结构因子的精度直接决定最后密度图的质量。
关键前提是把晶体结构精修到足够好。团队在100 K下做完整结构精修,观测与计算的衍射强度一致性因子只有0.97%,相当漂亮。只有基底足够干净,后面扣芯电子、画价电子等值面才不会被噪声带偏。

g波多极的实空间指纹:从X射线衍射到多极分解

当实验价电子密度图真正画出来时,Cr1和Cr2两类位点上的等值面呈现出明显的、沿c轴方向交替变化的各向异性:每个Cr周围有六个低密度“空洞”,正好指向最近邻的六个Sb配体。乍一看,这符合“离子键视角”的预期——价电子在指向配体的方向上被挖走,而不会像强共价图像那样在Cr-Sb连线上堆积。理论计算的价电子密度在侧视和俯视两个方向上都复现出了这种独特的花样。
图3:g波交替磁体中各向异性价电子密度的直接观测。(a, b)实验与理论计算得到的Cr位点VED侧视图。(c, d)Cr1层VED俯视图。(e)三角畸变D3d环境下Cr d轨道能级分裂。(f)eg与eg‘轨道混合形成的e_g^1-4本征态能级。(g, h)实验与拟合得到的不同角度下VED角分布对比。
图3:g波交替磁体中各向异性价电子密度的直接观测。(a, b)实验与理论计算的VED侧视图,两个Cr子晶格的价电子云沿c轴交替变化。(c, d)实验与理论给出的Cr1层VED俯视图,保留三重旋转特征。(e)三角畸变D3d晶场下Cr 3d轨道劈裂为a₁g、e_g、e_g′等能级。(f)e_g与e_g′轨道混合后形成的e_g¹⁻⁴本征态能级图。(g, h)以Cr1为球心、固定半径处的VED角分布实验图与拟合图,后面分析正是对这类角分布做多极展开。
这种看似“绕来绕去”的密度花样,本质上源于Cr d轨道在三角畸变晶体场下的轨道混合。CrSb中,Cr的3d轨道被六个Sb组成的三角畸变八面体包围,能级先分成a₁g、e_g和e_g′等几组。其中两组二重简并的e_g和e_g′能量差极小(约15 meV),实际本征态会互相“混血”,形成两组新的二重态e_g¹、e_g²和e_g³、e_g⁴。对于Cr1子晶格,第一组新本征态可以写成:
公式1:e_g^1 = α c1 |d_x²-y²⟩ + c2 |d_yz⟩。
公式1中|d_x²-y²⟩代表d轨道波函数,c1和c2表示两套轨道之间的混合权重,α=±1则是子晶格自由度:两个Cr子晶格靠一个垂直于c轴的镜面关联,所以α在这里会反号。
公式2:e_g^2 = α c2 |d_zx⟩ + c1 |d_xy⟩。
公式2中的e_g²由另一组d_zx和d_xy轨道混合而成,系数与公式1相互对称。重点来了:把e_g¹和e_g²的电子密度相加后,交叉干涉项会给出一个形状非常明确的贡献:
公式3:α{yz·(x²-y²) + 2zx·xy} = α yz(3x²-y²),即g波多极Q43的角分布。
公式3的右侧正是yz(3x²-y²),也就是论文反复强调的g波多极Q43的角分布。这里有一个容易被忽略的奥妙:如果e_g和e_g′两组轨道不混合(即c2=0),那电子云只会有d轨道级别(l=2)的普通四极形变;一旦混合,干涉项就会直接催生l=4的十六极分量。公式前的α表示它在Cr1和Cr2两个子晶格上符号相反,对应着反铁性g波多极序
为了进一步确认,团队做了严格的多极分解。对称分析表明,在不考虑自旋轨道耦合的顺磁d轨道流形中,只有四个电多极是“活跃”的:单极Q00、四极Q20和两个十六极Q40、Q43。由于两个Cr子晶格由镜面关联,其中只有Q43会在镜面操作下变号,因此也只有它能形成反铁性排列;其余三个多极在对称性上只能“整整齐齐”地铁性排列。换句话说,只要发现VED中有反铁性交替的高阶成分,几乎就是Q43在“实名出场”。
图4:g波交替磁体VED的多极分解。(a, b)两个Cr位点VED实验中黄色等值面及截面图。(c-f)CrSb中活跃的电多极Q00、Q20、Q40与Q43,只有Q43可以反铁性有序。(g)基于对称性紧束缚Wannier建模得到的q_lm及期望值。(h, i)包含所有Q_lm后重建的VED与Q43正瓣对比。(j, k)去掉Q43后重建结果与实验不一致。
图4:g波交替磁体VED的多极分解。(a, b)实验测得的两个Cr位点VED等值面与截面。(c-f)对称性允许的电多极Q00、Q20、Q40与Q43,其中Q43∝yz(3x²-y²)是唯一能形成反铁性有序的十六极。(g)基于对称性紧束缚Wannier建模得到的各多极系数q_lm与期望值。(h, i)包含所有多极后重建的VED,两个密度极大值正好对应Q43的正瓣。(j, k)将Q43去掉后重建结果与实验偏差明显,证明g波多极不可或缺。
定量一环则借助了对称性紧束缚Wannier建模(symmetry-adapted closest Wannier modeling,简称SCW)。它将第一性原理能带信息投影到完整多极基组上,从而把晶体场哈密顿量写成多极“算符”的线性叠加:
公式4:H(α) = q00 Q00 + q20 Q20 + q40 Q40 + q43 Q43(α),晶体场哈密顿量按多极算符展开。
公式4中的Q̂_lm(α)代表各大电多极对应的晶体场算符,q_lm是它们的系数。SCW给出的结果相当惊艳:g波多极系数q43约为0.14 eV,不仅和Q40同量级,甚至能跟低阶的Q20掰手腕。如果重建时强行去掉Q43,VED的“洞”会偏离Sb配体方向,两个密度极大值的位置也会错位——这就像拼图少了一块,整体花纹立刻对不上。

自旋密度的铁性有序:极化中子衍射的独特视角

价电子密度反映的是全部电子电荷的空间分布,可科学家真正想拍的,还有另一半:自旋密度(spin density,简称SD),也就是“自旋向上的电子比自旋向下的电子多出来的那部分空间分布”。X射线几乎不区分自旋,这时就得请出带磁矩的中子。极化中子衍射通过测量中子自旋向上与自旋向下两种散射强度的比值——也就是翻转比(flipping ratio,简称FR)——来提取非常微弱的磁性散射信号。
一个容易绕进去的反转点在于:电荷密度里的Q43是反铁性排列的,到了自旋密度里,Q43反而变成铁性排列。原因是时间反演操作对磁多极额外贡献一个负号,把“反铁性”和“铁性”的阵营整体互换了。物理上更直观的理解是:两个Cr位点的磁偶极本来就是反平行(反铁性)的,要让整体磁多极在相加时不抵消、出现铁性的高阶磁多极,电荷那部分Q43就必须反铁性地“配合”,两者锁定在一起,才同时给出非零的铁性磁多极。自旋密度ρ_s(α)可以近似写成α·ρ(α),α本身就在两个子晶格之间交替。
图5:极化中子衍射探测磁多极的铁性有序。(a)300 K下用各向同性自旋密度模型拟合翻转比FR,偶l反射计算值为1,实验数据明显偏离。(b)加入铁性磁多极Q43的分量后,模型与l=even反射数据吻合更好,R_FR从3.48%降到2.13%。
图5:极化中子衍射探测磁多极的铁性有序。(a)300 K下仅用各向同性磁偶极模型拟合翻转比FR时,偶指标l=even反射的预测值严格为FR=1,但实验数据明显偏离,说明高阶磁多极确实存在。(b)纳入Q43所对应的铁性磁多极分量后,l=even反射的拟合显著改善,整体拟合优度R_FR从3.48%下降到2.13%。
实验数据在300 K下呈现一个清晰判据:若CrSb只是普通反铁磁体,那么磁偶极只影响l为奇数的反射,l为偶数的反射预期翻转比严格等于1。实际测量中,l=even反射的FR显著偏离1——这是铁性高阶磁多极“反常规”存在的直接证据。把Q43相关分量加进模型后,拟合优度从3.48%改善到2.13%,且Q43的拟合值为0.21(6),与理论预测的0.25吻合良好。两条独立实验线在此汇合:电多极反铁性有序 + 磁偶极反铁性有序 → 磁多极铁性有序,这正是交替磁体的核心指纹。

从实空间到动量空间:多极与自旋劈裂的深层联系

为什么一定要分清“哪个多极在反铁性排列”?因为实空间的密度花纹和动量空间能带的自旋劈裂,很可能是同一个物理根源的两张脸。CrSb早前被角分辨光电子能谱测出巨大的g波自旋劈裂,形状为k_y k_z(3k_x²-k_y²)。这种“倒空间风车”般的劈裂,常见解释聚焦于相邻Cr位点之间各向异性的hopping(跳跃积分)项。但本文指出,同一原子内的多极晶体场与子晶格间跳跃配合,同样可以产生g波劈裂。
微观模型计算显示,最低阶的劈裂项正比于:
公式5:g(k) ∝ J_z q43 t² k_z k_y (3k_x² - k_y²),自旋劈裂正比于多极系数q43。
公式5中J_z是交换相互作用,q43是刚才从SCW分析中得到的原子多极强度,t则是子晶格之间的最近邻跳跃振幅。这个式子把“看到的高阶多极”和“已知的自旋劈裂”直接连成一个乘积:q43越大,劈裂越大。此前不少理论把g波劈裂完全归因于跳跃的g波对称性,而论文展示了一条互补路径——就算hopping本身不含g波各向异性,只要原子位点上有足够强的Q43晶体场,再叠加上层间耦合t,动量空间的g波劈裂照样会被“点亮”。这为理解CrSb实验中电子伏级别的巨劈裂提供了新的微观抓手。

展望:多极成像策略的广泛应用前景

这篇论文更大的价值,是给出了一套可以“搬到别处”的方法论:高质量的同步辐射X射线衍射负责捕捉价电子密度,CDFS负责把芯电子背景抹掉,Wannier多极建模负责把密度拆成序参量,极化中子衍射负责给磁性那一半做交叉验证。这套组合拳不挑食,原则上可以推广到d波交替磁体、超胞交替磁体,甚至能推向那些翻转对称性破缺的奇宇称磁体,以及f电子体系里困扰人们已久的各种隐藏多极序。
对应用端来说,交替磁体之所以受宠,是因为它结合了反铁磁体“没有杂散场、抗磁场干扰、读写快”的优点,又拥有铁磁体那种自旋劈裂的操控性。如果能通过晶体对称性工程进一步操控g波多极排列的方向,CrSb这类室温材料在下一代自旋存储和太赫兹自旋电子学中的想象空间会相当可观。当然,这次工作还只是“看清了秩序”,距离像操纵铁磁畴一样操纵多极畴,中间还有不少路要走。

龙迷三问

下面是龙哥对于大家可能的一些问题的解答:
这篇论文到底在解决什么问题?研究用同步辐射X射线实测CrSb价电子密度,首次在实空间直接看到g波多极(十六极)在两个Cr位点上的反铁性交替;再结合极化中子证明铁性磁多极。该技术可为量子材料“隐藏序”探测提供通用路径。
这篇工作最值得看的点是什么?论文通过实验直接观察到了CrSb中g波交替磁多极的价电子密度各向异性,并通过极化中子衍射证实了铁性磁多极的存在,实验结果与理论计算高度一致。
这篇工作的边界或风险在哪里?优点:首次直接可视化并量化了g波交替磁多极,建立了交替磁性的实空间指纹;结合X射线衍射、中子衍射和理论建模,方法全面且系统。缺点:VED分析在Sb位点因5p轨道高度离域而无法重构;实验需要大型同步辐射和中子源设施,技术门槛高。
如果你还有哪些想要了解的,欢迎在评论区留言或者讨论~

龙哥点评

论文创新性分数:★★★★☆

通过高能同步辐射X射线衍射结合核心差分傅里叶合成方法直接可视化CrSb中Cr位点的价电子密度分布,并利用多极展开和对称性适配Wannier建模量化g波交替磁多极,同时结合极化中子衍射探测铁性磁多极。

实验合理度:★★★☆☆

现有材料未完整覆盖数据划分、基线公平性和统计显著性,因此按中性评价处理。

学术研究价值:★★★★☆

通过高能同步辐射X射线衍射结合核心差分傅里叶合成方法直接可视化CrSb中Cr位点的价电子密度分布,并利用多极展开和对称性适配Wannier建模量化g波交替磁多极,同时结合极化中子衍射探测铁性磁多极;更关键的是问题定义是否可复用到同类任务。

稳定性:★★★☆☆

现有材料未提供充分的极端条件、重复运行或扰动测试,稳定性暂按中性评价。

适应性以及泛化能力:★★★☆☆

现有材料未完整展示跨数据集、跨场景或分布外实验,泛化能力仍需进一步验证。

硬件需求及成本:★★★☆☆

现有材料缺少完整训练资源、参数量、显存和推理时延信息,成本暂按中性评价。

复现难度:★★★☆☆

现有材料未确认完整代码、配置、数据处理脚本和权重是否齐备,复现难度暂按中性评价。

产品化成熟度:★★★☆☆

论文验证以研究实验为主,真实部署中的时延、成本、维护和异常场景仍需补充验证。

可能的问题:VED分析在Sb位点因5p轨道高度离域而无法重构;实验需要大型同步辐射和中子源设施,技术门槛高。

主要参考文献

[1] R. Misawa, R. Oiwa, S. Kitou, et al. Observation of g-wave altermagnetic multipole, arXiv:2609.01969v1 (2026).
[2] L. Šmejkal, J. Sinova, T. Jungwirth. Emerging research landscape of altermagnetism, Phys. Rev. X 12, 040501 (2022).
[3] X. H. Verbeek, D. Voderholzer, S. Schären, et al. Nonrelativistic ferromagnetotriakontadipolar order and spin splitting in hematite, Phys. Rev. Res. 6, 043157 (2024).
[4] S. Reimers, L. Odenbreit, L. Šmejkal, et al. Direct observation of altermagnetic band splitting in CrSb thin films, Nat. Commun. 15, 2116 (2024).


end
欢迎加入龙哥读论文粉丝群,扫描下方二维码或者添加龙哥助手微信号加群:kangjinlonghelper。一定要备注:研究方向+地点+学校/

*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击"阅读原文",查看更多原论文细节哦!       

转发文章 微博 X LinkedIn Facebook
龙哥读论文 · PaperDaily

本文基于龙哥读论文 PaperDaily 数据库整理,结合论文原文与工程视角进行解读。