RuO2的交错磁性

 

PHYSICAL REVIEW B (2026)

RuO2的交错磁性:从争议到共识

Altermagnetism in RuO2: A Comprehensive Review

导读:RuO2是交错磁性(Altermagnetism)研究的"氢原子"--最简单的AM候选材料,承载着AM领域最核心的争议与共识。2024年,Nature同期发表两篇论文验证RuO2的AM特征,但随后中子散射实验发现其磁矩仅~0.05 muB(远小于DFT预测的~1 muB),引发激烈讨论。本文系统梳理RuO2的AM研究现状:从晶体结构与对称性分析,到磁基态争议,到自旋劈裂的d波特征,再到拓扑性质(Dirac节点线、Weyl点)和反常输运(AHE、SHE、晶体Hall效应),最后展望AM自旋电子学器件应用。

图片

一、前言背景

RuO2:交错磁性的"氢原子"

RuO2(二氧化钌)是金红石型氧化物,长期以来被视为最简单的反铁磁体。但2022年以来,随着交错磁性(Altermagnetism)概念的提出,RuO2被重新审视并成为AM研究的"氢原子"--它是最简单的AM候选材料,具有最清晰的对称性分析。

RuO2的金红石结构(空间群P4_2/mnm)包含两个不等价的Ru子晶格,由非 symmorphic 的四度螺旋轴4_2连接。在Neel温度以上为顺磁金属,Neel温度以下为共线反铁磁体。磁矩沿[001]方向排列,净磁矩为零。

但RuO2的AM身份存在争议:2024年多篇实验论文声称观察到d波自旋劈裂和非相对论AM特征,但后续中子散射实验发现RuO2的磁矩远小于预期(~0.05 muB而非~1 muB),挑战了长程共线AFM序的描述。

核心科学问题:RuO2是AM、AFM还是其他?

本文是一篇综述性论文,系统梳理了RuO2的AM研究现状和争议。核心问题包括:(1) RuO2的基态磁序到底是什么?(2) 自旋劈裂的起源是非相对论AM效应还是SOC+磁空间群效应?(3) 实验观测到的反常输运(AHE、自旋Hall效应)是否来自AM序?

关键争议点:(1) 磁矩大小--DFT预测~1 muB,但中子散射测量~0.05 muB;(2) 自旋劈裂的大小和节点结构--不同实验给出的d波劈裂幅值差异显著;(3) 缺陷效应--Ru空位可能改变磁基态。

综述涵盖:晶体结构与对称性、磁基态的理论与实验对比、电子结构(自旋劈裂)、拓扑性质(Dirac节点线)、输运性质(AHE、自旋Hall效应、晶体Hall效应)、超导邻近效应、以及RuO2在自旋电子学器件中的应用前景。

figure

RuO2交错磁性研究综述。晶体结构(金红石P4_2/mnm)-> 磁基态(Neel AFM)-> 自旋劈裂(d波)-> 拓扑性质(Dirac节点线)-> 输运性质(AHE/SHE)-> 超导邻近效应(RuO2/TiO2)-> 器件应用。核心争议:磁矩大小(DFT vs 实验)和自旋劈裂起源(非相对论AM vs SOC+磁空间群)。

二、核心内容

图片

图 1:RuO2的晶体结构和磁构型。金红石结构,两个Ru子晶格由4_2螺旋轴连接,形成共线Neel AFM序。

晶体结构与对称性:为什么RuO2是AM的"氢原子"

RuO2的金红石结构(P4_2/mnm, No. 136)是AM对称性分析的最简单范例。关键对称性操作:两个Ru子晶格由四度螺旋轴4_2 = {C4z|t(0,0,1/2)}连接--即绕z轴旋转90度后沿z方向平移半个晶格常数。

AM判据:两个Ru子晶格的自旋方向相反(共线AFM),但由晶体旋转操作(而非PT或纯平移+反演)连接。这满足AM的定义:补偿磁序(净磁矩为零)+ 非相对论自旋劈裂。

自旋劈裂的对称性:在RuO2中,4_2操作将Ru1(^)映射到Ru2(v),但同时在倒空间中产生d波形式的自旋劈裂。自旋劈裂的节点沿[110]和[1-10]方向,劈裂幅值在[100]和[010]方向最大。

图片

图 2:RuO2的电子结构。(a) 无SOC能带结构。(b) 含SOC能带结构。(c) 自旋分辨能带,展示d波自旋劈裂。(d) DOS和轨道投影。

电子结构:自旋劈裂的大小与节点

RuO2的能带结构在费米面附近主要由Ru-4d t2g轨道(dxy, dyz, dxz)贡献。无SOC时,能带存在非相对论自旋劈裂--这是AM的标志性特征。SOC进一步打开某些能带交叉处的能隙,改变自旋劈裂的细节但非本质。

自旋劈裂幅值:不同DFT计算给出的自旋劈裂值在100-500 meV之间,取决于泛函(PBE vs HSE vs GW)和U值。实验上,自旋分辨ARPES测得的自旋劈裂小于理论预测,引发争议。

节点结构:d波自旋劈裂有4个节点,沿[110]和[1-10]方向。在节点处,自旋向上和向下的能带简并。节点的存在是AM区别于FM(无节点)和传统AFM+SOC(无节点或节点结构不同)的关键特征。

图片

图 3:磁基态争议。(a) DFT计算的不同磁构型能量。(b) 中子散射实验。(c) 磁矩的温度依赖性。

磁基态争议:DFT vs 中子散射

DFT计算(PBE/PBE+U/HSE)一致预测RuO2的Neel AFM基态,Ru磁矩约0.8-1.2 muB。但中子散射实验测得的磁矩仅为~0.05 muB,两者相差约20倍。

可能的解释:(1) 强量子涨落压低了磁矩--RuO2接近金属-绝缘体转变,量子涨落效应强烈;(2) 缺陷(Ru空位、O空位)改变了磁基态;(3) 磁有序可能不是长程共线的,而是短程或非共线的;(4) 中子散射可能探测到了不同的磁关联长度。

这一争议是AM领域的核心问题之一。如果RuO2的磁矩确实如此之小,其作为AM候选材料的地位将受到挑战:AM效应(自旋劈裂、AHE)的强度与磁矩大小成正比。

图片

图 4:拓扑性质。(a) Dirac节点线在k空间的位置。(b) Weyl点的分布。(c) Berry曲率。

拓扑性质:Dirac节点线、Weyl点和Berry曲率

RuO2在SOC下展现出丰富的拓扑能带结构。在金红石结构的高对称面上,能带交叉形成Dirac节点线。这些节点线受非symmorphic对称性保护,对微扰鲁棒。

当磁矩沿[001]方向时,某些Dirac节点线在SOC下打开能隙,形成Weyl点。Weyl点的手性电荷产生Berry曲率,驱动反常输运效应。

注意:RuO2的Dirac节点线同时存在于无SOC和含SOC的能带中,但它们的起源不同。无SOC的节点线由能带对称性(轨道自由度)保护,含SOC的节点线由磁空间群对称性保护。

图片

图 5:反常Hall效应和自旋Hall效应。(a) AHE电导率。(b) SHE电导率。(c) 晶体Hall效应。

反常输运:AHE、SHE和晶体Hall效应

反常Hall效应(AHE):在RuO2中,由于SOC产生的Berry曲率,施加面内电场会产生面外Hall电压。AHE的大小与磁矩方向相关,可用于电学读取AM的磁态。

自旋Hall效应(SHE):RuO2中自旋轨道耦合产生自旋Hall效应,自旋向上和向下的电子向相反方向偏转,产生纯自旋流。AM的自旋劈裂可能增强SHE。

晶体Hall效应(crystal Hall effect):这是AM特有的输运效应--在零净磁矩下,由于自旋劈裂的节点结构,不同自旋通道的Hall电导不完全抵消,产生非零的净Hall响应。这是区分AM和传统AFM的"smoking gun"信号。

图片

图 6:RuO2/TiO2超导邻近效应。(a) 异质结结构。(b) 超导能隙。(c) Andreev束缚态。

超导邻近效应:AM/超导体界面的新物理

将RuO2与超导体(如TiO2-x)接触,可以在界面处诱导超导邻近效应。AM的自旋劈裂对超导配对对称性产生独特影响:不同自旋通道的Andreev反射不对称,导致自旋极化的超导电流。

这一效应为AM自旋电子学器件提供了新思路:利用AM/超导体界面实现自旋极化超导电流的注入和检测。但RuO2/TiO2界面的质量(缺陷、互扩散)是实现这一目标的关键挑战。

与FM/超导体界面的对比:FM/超导体界面中,交换场抑制s波超导配对,产生Fulde-Ferrell-Larkin-Ovchinnikov(FFLO)态。AM/超导体界面中,由于自旋劈裂的节点结构,FFLO态的稳定性可能不同。

figure

图 7:AM自旋电子学器件。(a) 自旋劈裂的电场调控。(b) AM隧道结。(c) 自旋注入/检测方案。

RuO2自旋电子学器件:AM的优势与挑战

AM器件相比传统FM器件的优势:(1) 零杂散场--不对相邻器件产生磁干扰;(2) THz级自旋动力学--反铁磁共振频率远高于FM的GHz级;(3) 自旋劈裂可被电场/应力调控--提供了额外的操控自由度。

RuO2的具体器件方案:(1) AM隧道结--利用自旋劈裂实现自旋过滤;(2) 自旋劈裂FET--通过栅压调控费米面位置,选择性地开关自旋极化电流;(3) AM/重金属双层--利用自旋Hall效应实现自旋轨道转矩翻转。

主要挑战:RuO2薄膜的制备质量(金红石相 vs 其他相)、磁矩太小导致的信号微弱、以及器件工作温度(需要低于Neel温度)。

figure

图 8:RuO2的AM研究时间线。(a) 关键发现的时间线。(b) 理论预测与实验验证的对比。(c) 未来展望。

DFT Tips

【DFT Tip 1】RuO2的DFT+U:U值选择与磁矩敏感性

RuO2的DFT+U计算中,U值的选择对磁矩大小有显著影响。PBE无U时磁矩约0.5 muB,PBE+U(J=0, U=2-4 eV)时磁矩约0.8-1.2 muB。但实验磁矩~0.05 muB远小于DFT预测。

建议:对RuO2做系统U值扫描(U=0,1,2,3,4,5 eV),分析磁矩、能带结构和自旋劈裂的演化。U值越大磁矩越大,但这是DFT+U方法的已知偏倚--它倾向于过分局域化电子。

更精确的方法:使用HSE06杂化泛函或GW方法,它们不依赖经验的U参数。但HSE06和GW对RuO2可能预测比PBE+U更小的磁矩,更接近实验。

【DFT Tip 2】非symmorphic对称性与Dirac节点线

RuO2的P4_2/mnm空间群包含非symmorphic对称性操作(4_2螺旋轴和n-glide镜面)。这些操作在倒空间中产生能带简并,形成Dirac节点线。

在DFT中识别节点线:(1) 做密集k点能带计算,寻找能带交叉点;(2) 使用WannierTools或irrep分析能带的对称性标记;(3) 计算节点线附近的Berry相(pi或0)来验证拓扑非平庸。

常见陷阱:非symmorphic对称性在slab模型中可能被破坏,导致节点线在表面计算中消失。对于薄膜/表面计算,必须使用包含完整对称性的体材料模型。

【DFT Tip 3】AM自旋劈裂的定量表征:劈裂幅值与节点

AM自旋劈裂的定量分析需要:(1) 计算自旋分辨能带,提取劈裂能量Delta E(k) = E_up(k) - E_down(k);(2) 在k空间绘制Delta E(k)的分布图,识别节点位置和劈裂的波对称性(d波/g波/i波);(3) 计算劈裂的k空间平均值和最大值。

VASP中自旋分辨能带计算:使用LNONCOLLINEAR=.TRUE. 和LSORBIT=.TRUE.,在能带计算中输出自旋分量。后处理时提取每个k点的自旋向上和向下能带能量。

注意:SOC下的自旋劈裂包含非相对论AM贡献和相对论SOC贡献的混合。分离两者的方法:分别计算无SOC和含SOC的能带,差值即为纯SOC贡献。

【DFT Tip 4】AHE和SHE的DFT计算:Berry曲率与输运系数

AHE和SHE的DFT计算流程:(1) 构建Wannier紧束缚模型(含SOC);(2) 在密集k网格上计算Berry曲率Omega_n(k);(3) 积分得到AHE电导率sigma_xy = -(e^2/h) * sum_n integral Omega_n(k) f(epsilon_n(k)) d^2k。

SHE的计算更复杂,需要自旋分辨的Berry曲率。Wannier90和WannierTools可以同时计算AHE和SHE。对于AM材料,特别注意晶体Hall效应的贡献--它来自不同自旋通道的Hall电导的不完全抵消。

k点收敛:AHE对k点采样极其敏感,Berry曲率在能带交叉点附近发散。需要极密的k网格(如1000x1000x1000)才能收敛,这只有通过Wannier插值才能实现。

【DFT Tip 5】RuO2中U值的争议:PBE vs HSE vs GW vs DMFT

计算层级对比:(1) PBE:计算成本最低,但倾向于低估磁矩和带隙(RuO2在PBE下为金属);(2) PBE+U:成本低,磁矩可调,但U值选择主观;(3) HSE06:成本中等,包含部分精确交换,可能给出更合理的磁矩;(4) GW:成本高,能带和磁矩最准确但计算量极大;(5) DFT+DMFT:成本极高,同时处理局域和巡游关联,是解决RuO2磁矩争议的最可靠方法。

当前RuO2的DFT+DMFT结果:预测磁矩约0.3-0.5 muB,介于PBE+U(~1 muB)和实验(~0.05 muB)之间。DMFT中的动态自旋涨落压低了平均磁矩。

建议:对于RuO2的计算,至少使用PBE+U和HSE06两种方法,比较结果的一致性。如果资源允许,DFT+DMFT是解决争议的最有力工具。

【DFT Tip 6】缺陷对RuO2磁性的DFT建模

实验中的RuO2薄膜和晶体通常含有Ru空位和O空位。缺陷对磁性的影响至关重要:(1) Ru空位可能破坏补偿磁序,产生净磁矩;(2) O空位引入额外电子,可能改变Ru的价态和磁矩。

DFT缺陷建模:(1) 构建超胞(至少2x2x2),引入单个缺陷;(2) 弛豫缺陷周围的原子位置;(3) 计算缺陷形成能和磁矩的变化。注意:缺陷浓度在DFT超胞中通常远高于实验,可能高估缺陷效应。

常见陷阱:带电缺陷的DFT计算需要背景电荷补偿(jellium background),这可能引入人工的静电效应。对于金属性RuO2,带电缺陷的处理相对简单(电荷被巡游电子屏蔽)。

【DFT Tip 7】AM能带结构的对称性分析工具

AM材料的对称性分析需要自旋空间群(spin group)而非传统的磁空间群(magnetic group)。关键工具:(1) spglib(支持自旋空间群的最新版本);(2) irrep(能带不可约表示分析);(3) Bilbao Crystallographic Server(磁对称性分析)。

自旋空间群与传统磁空间群的区别:自旋空间群允许自旋空间和实空间的独立操作,而磁空间群只允许耦合操作。AM的对称性只能在自旋空间群中完整描述。

实际操作:在DFT计算中,首先通过MAGMOM设置目标磁构型,然后使用对称性分析工具确认该磁构型属于哪个自旋空间群。

【DFT Tip 8】金红石结构氧化物的DFT计算规范

金红石结构(TiO2, RuO2, CrO2, MnO2等)的DFT计算有两个关键点:(1) 四方晶胞的c/a比的准确弛豫--c/a比影响能带结构和磁交换耦合;(2) 氧八面体的旋转/畸变模式--金红石结构中可以出现软模。

对于RuO2,建议使用以下计算设置作为起点:ENCUT=520 eV, k点8x8x12(体材料),ISIF=3(完全弛豫),EDIFF=1E-6, EDIFFG=-0.001。

实验晶格常数:a=4.49 A, c=3.11 A, c/a=0.69, u=0.305(氧位置参数)。DFT弛豫应接近这些值,偏差>2%需要检查计算设置。

知识扩展

【知识扩展 1】自旋空间群 vs 磁空间群:AM的对称性语言

【理论解释】磁空间群(magnetic space group)是描述磁性材料对称性的传统框架,它将自旋反转(时间反演T)与晶格操作耦合。但磁空间群只允许自旋翻转(spin flip),不允许自旋旋转(spin rotation)。自旋空间群(spin space group)允许自旋空间和实空间的独立操作,更完整地描述了磁性材料的对称性。

【AM的对称性本质】在磁空间群中,AM的对称性被不完整描述:自旋劈裂看似"破缺"了某些对称性。但在自旋空间群中,AM保留了C4T等联合操作,自旋劈裂是自然的对称性结果,不破缺任何对称性。

【经典参考】Smejkal et al., PRX 12, 031042 (2022) -- AM的对称性理论;Chen et al., Nature 640, 349 (2025) -- 自旋空间群实验验证;Jiang et al., PRX 14, 031039 (2024) -- 自旋空间群枚举。

【迁移能力】自旋空间群的概念适用于所有AM材料,也可用于理解非共线磁结构(如skyrmion晶格)的对称性。

【知识扩展 2】晶体Hall效应:AM的"smoking gun"信号

【理论解释】晶体Hall效应(crystal Hall effect)是AM特有的反常输运效应。在传统FM中,AHE正比于净磁矩。在传统AFM中,AHE为零(两个子晶格贡献抵消)。在AM中,尽管净磁矩为零,但由于自旋劈裂的节点结构,不同k点的Hall电导贡献不完全抵消,产生非零净Hall响应。

【实验意义】晶体Hall效应是区分AM和传统AFM的最直接电学信号。在AFM中,任何Hall信号都必须在SOC下产生,且与磁矩方向无关。在AM中,Hall信号与Neel矢量方向耦合,且可以在无SOC极限下存活。

【经典参考】Smejkal et al., PRX 10, 011035 (2020) -- 晶体Hall效应理论;Feng et al., Nat. Electron. 5, 735 (2022) -- RuO2中的晶体Hall效应实验。

【迁移能力】晶体Hall效应是AM材料的通用输运特征,可用于筛选和验证AM候选材料。

科研经验

【科研经验 1】AM材料验证的"黄金标准":理论-实验的多重对应

问题:仅凭自旋分辨ARPES看到自旋劈裂,是否足以声称一个材料是AM?

原因:传统AFM在SOC下也能产生自旋劈裂,但这是相对论效应(~10 meV),而非AM的非相对论效应(~100 meV)。仅凭自旋劈裂的存在无法区分AM和AFM+SOC。

解决方案:AM材料的验证需要多重量子:(1) 无SOC的DFT能带展示显著自旋劈裂;(2) 自旋劈裂的节点结构与AM对称性一致(d波/g波/i波);(3) 晶体Hall效应的电学测量;(4) 自旋分辨ARPES测量自旋劈裂的k空间分布;(5) 中子散射确认补偿磁序。任何一个单独证据都不足以确认AM。

建议:RuO2的争议正是由于不同实验给出了相互矛盾的结果。在AM研究领域,理论预测和实验验证之间的闭环至关重要。

【科研经验 2】DFT磁矩与实验磁矩不一致的应对策略

问题:DFT预测RuO2磁矩~1 muB,但中子散射测量~0.05 muB。这种量级差异在DFT研究中如何处理?

原因:DFT+PBE+U倾向于高估磁矩,因为它将电子过度局域化。实验上,RuO2接近金属-绝缘体转变,量子涨落和巡游效应强烈,局域磁矩图像可能不适用。

解决方案:(1) 使用更高级的方法(HSE06, GW, DMFT)校核磁矩;(2) 计算磁矩的温度依赖性(Monte Carlo + 交换耦合);(3) 在论文中明确讨论DFT预测与实验差异的可能原因,而非回避。

建议:不与实验比较绝对值,而是比较趋势和排序。例如,DFT预测的磁矩的各向异性(不同方向磁矩的差异)可能与实验更一致。

 

RuO2 Altermagnetism Review (2026) | 金红石 磁空间群 自旋劈裂 晶体Hall效应