交变磁性(altermagnetism)作为第三种基本磁相,近年来在凝聚态物理和自旋电子学领域引起了广泛关注。本文解读的这项研究以 CrSb 为模型体系,通过对称性分析和第一性原理计算,揭示了对称性破缺如何影响交变磁性的形成,并首次发现了"碎片化节线曲线"(Fragmented Nodal Curves, FNCs)这一新奇量子现象,为交变磁体在量子器件中的应用开辟了新途径。

背景
交变磁性(Altermagnetism)是继铁磁性和反铁磁性之后的第三种基本磁相。与传统的反铁磁体不同,交变磁体在动量空间中展现出独特的自旋劈裂现象------即所谓的"交变磁自旋劈裂"(Altermagnetic Spin Splitting, AMSS)。这种劈裂类似于铁磁体中时间反演对称性破缺导致的能带劈裂,但在节面(nodal plane)两侧大小相等、方向相反,使得在整个布里渊区积分后净磁化为零,保持了反铁磁体的特征。
交变磁性是一种纯粹的非相对论量子现象,其动量空间自旋极化(Momentum Space Spin Polarization, MSSP)既不同于经典铁磁体,也不同于反铁磁体,为自旋电子学领域开辟了全新的量子输运途径。
CrSb 是研究最为广泛的交变磁体之一,属于 NiAs 型六方晶系(空间群 P63/mmc),具有以下显著特征:
◆ 巨大的交变磁自旋劈裂,量级可达约 1 eV
◆ 大的磁晶各向异性能(MAE),有利于信息存储
◆ 薄膜中 MAE 降低,便于操纵 Neel 矢量
◆ 可通过应变、电场、掺杂或电压重新定向 Neel 矢量
◆ 交变磁相、铁磁相和铁电相之间存在适中的能量势垒,适合多铁自旋电子学
◆ 最近 ARPES 实验证实了交变磁 Weyl 节点和拓扑费米弧的存在
这些特性使 CrSb 成为研究交变磁性形成机制以及通过对称性工程调控其性质的理想平台。
方法
第一性原理计算
计算框架:
▸ VASP(Vienna Ab-initio Simulation Package)
采用平面波基组的投影缀加波(PAW)方法,使用 PBE 广义梯度近似作为交换关联泛函。动能截断为 500 eV,布里渊区积分采用 Gamma 中心的 12x12x6 k 点网格。
▸ QUANTUM ESPRESSO + WANNIER90
用于计算反常霍尔电导率(AHC),采用模守恒赝势并包含自旋轨道耦合。通过最大局域化 Wannier 函数构建紧束缚哈密顿量,利用 WANNIERBERRI 代码计算 AHC。
▸ DFT+U 方法
对 RbMnPO4 的计算中,使用 Dudarev 的旋转不变泛函框架,有效 Hubbard 参数 Ueff = 3 eV。
模型结构设计
为系统研究对称性破缺对交变磁性的影响,从原始 CrSb(Cr2Sb2)出发,通过假设性的空位工程和间隙掺杂,设计了五种模型结构(Model Structures, MS-I 至 MS-V)

MS-I(Cr2Sb1)
移除一个 Sb 原子,空间群降至 P-6m2,自旋空间群 P-1-61m-12∞m1
。两个磁性子晶格通过 C2\|\|Mz 和 C2\|\|S6z 对称性连接。
MS-II(Cr2Sb3)
在六方空隙中心额外添加一个 Sb 原子,形成 Cr2Sb3。空间群和自旋空间群与 MS-I 相同,同样具有 S6z 旋转反演对称性。
MS-III(Cr2Sb4, 构型A)
在两个次近邻六方中心各添加一个 Sb。存在反演中心,导致体系为纯反铁磁体,无交变磁性特征。
MS-IV(Cr2Sb4, 构型B)
与 MS-III 化学式相同但原子排布不同。同样存在反演中心,为纯反铁磁体。
MS-V(Cr2Sb4, 构型C)
两个 Sb 占据最近邻六方中心位置,空间群 P21/m,自旋空间群 P-121/-1m∞m1
。仅有 C2z 和 Mz 对称性。

理论推导与对称性分析
在非相对论极限(Non-Relativistic Limit, NRL)下,时间反演对称性作用于动量空间,同时对自旋施加恒等操作,这一双重操作记为 E\|\|TR,L。其效果等同于空间反演算符 I:

因此,在 NRL 下,C6z 和 S6z(六重旋转反演)操作在动量空间中表现得完全相同。

对于具有 C6z 对称性的原始 CrSb,波函数在 C2\|\|C6z 和 E\|\|TR,L 联合操作下的变换为:

而对于具有 S6z 对称性的 MS-I 和 MS-II,波函数变换为:

两种情况下得到相同的能量本征值关系:

在决定交变磁性特征时,真正旋转(proper rotation)C6z 和非真旋转(improper rotation)S6z 在倒空间中产生完全等价的效果,都导致四个节面的出现------一个基面(kz=0)和三个对角节面。
当对称性从 C6z/S6z 降低到 C2z 时,情况发生了根本性变化。MS-V 具有 C2\|\|Mz 和 C2\|\|C2z 对称性,波函数变换如下:


结合 NRL 下的时间反演操作,得到支配 FNC 形成的核心能量关系:

这个关系式表明:对于任意给定的 kz,自旋相反的两个子晶格能带对必须在 kx-ky 平面内的某个特定点(节点)相交。当 kz 连续变化时,这些节点便形成了一条"碎片化节线曲线"(FNC)。与传统的节面不同,FNC 的位置是能带依赖的------每一对能带都有自己独特的 FNC 形状。这一现象在此前的文献中从未被报道过。

结果
原始 CrSb 的电子结构
原始 CrSb 的基态展现出典型的交变磁体能带结构,在对角节面两侧具有大小相等、方向相反的 AMSS。Cr 位的局域自旋矩约为 2.73 μB,表明高自旋构型。从化学键角度看,交变磁能带色散和强 AMSS 主要归因于 Sb-p 轨道之间第二近邻相互作用的决定性角色。
原始 CrSb 的节面示意图(a)及布里渊区高对称 k 路径(b)

原始 CrSb 在不同能量下的等能面,展示了四个节面的存在

模型结构的电子结构演化
随着从 MS-I 到 MS-V 对称性的逐步降低,电子结构发生了系统性变化:
▸ MS-I(Cr2Sb1)
移除一个 Sb 后,更多电子局域在 Cr 位,共价杂化减弱,带宽变窄。局域自旋矩增大至约 3.5 μB。AMSS 强度比原始体系小约一个数量级,归因于 Sb-Sb 第二近邻相互作用的破坏。
▸ MS-II(Cr2Sb3)
额外的 Sb 贡献更多价电子,能带更窄。局域自旋矩降至约 1 μB。与 MS-I 一样,S6z 和 Mz 对称性保证了四个节面的存在,但 AMSS 强度同样减弱。
▸ MS-III 和 MS-IV
由于反演中心的存在,自旋相反的子晶格能带完全简并,无交变磁性特征。
▸ MS-V
这是最重要的发现:仅有 C2z 和 Mz 对称性,三个对角节面消失,取而代之的是遍布整个布里渊区的 碎片化节线曲线(FNCs)。基面节面(kz=0)和沿 kz 轴的节线仍然保留。

MS-II 中 Cr 原子的自旋密度轮廓,展示了 π/6 旋转和 Mz 对称性

MS-V 中等能面的不同视图,清晰展示了对角节面的消失和 FNCs 的形成
FNCs 的形成并非仅限于假设模型结构------通过对原始 CrSb 施加沿 x 轴的单轴应变,可以在实际材料中实现这一现象。在 5% 应变下(a 从 4.11 Å 变为 4.32 Å),C6z 对称性被打破,体系保留了 C2z 和 Mz 对称性,等能面展现出与六重对称性的明显偏离和 FNCs 的形成。进一步的应变增大将导致节线曲线更加碎片化。

应变 CrSb 中的交变磁性和 FNCs 形成
为进一步验证 FNCs 的普适性,在对称性较低的反铁磁体,找到了实验合成的 RbMnPO4(空间群 P21,仅具有 C2y 对称性)。DFT+U 计算表明,RbMnPO4 的等能面同样展现出清晰的 FNCs 特征,证实了 FNCs 是 C2 旋转对称性下的普遍现象。

RbMnPO4 中的 FNC 驱动交变磁性
反常霍尔效应(AHE)的对称性调控
反常霍尔电导率(Anomalous Hall Conductivity, AHC)是交变磁体在自旋电子学应用中最重要的输运性质之一。通过 Kubo 公式系统地计算了不同对称性下的 AHC:

其中 Berry 曲率:

对称性与 AHC 的关系
AHC 的存在与否完全由体系的磁空间群(MSG)对称性决定。系统分析了原始 CrSb 和各模型结构在 Neel 矢量沿面外(c 轴)和面内(与 a 轴成 30°)两种取向下的 AHC 行为:
▸ 原始 CrSb(面外 Neel 矢量)
MSG 为 P6'3/m'm'c,对称性迫使 Berry 曲率在整个布里渊区抵消,AHC 为零。
▸ 原始 CrSb(面内 Neel 矢量)
MSG 降至 Cm'c'm,允许非零的 σzxy 分量。
▸ MS-I / MS-II(面外 Neel 矢量)
与原始情况类似,所有 AHC 分量消失。
▸ MS-I / MS-II(面内 Neel 矢量)
出现 σxyz 分量。
▸ MS-V / 应变 CrSb(面外 Neel 矢量)
同时出现 σxyz 和 σzxy 两个非零分量------这是对称性降低的最显著效应。
▸ MS-V / 应变 CrSb(面内 Neel 矢量)
出现 σyzx 分量。

最引人注目的发现是:MS-V 和应变 CrSb 在 Neel 矢量沿面外方向时仍能产生有限的 AHC,而原始 CrSb 在此取向下 AHC 为零。这意味着通过应变工程,可以在保持 Neel 矢量沿易轴(面外)的同时获得非零 AHC,这对于实际器件应用具有重大意义。

CrSb 体系在不同磁空间群和 Neel 矢量取向下的对称性分析与 AHC 张量分量
核心结论
- 非真旋转与真旋转的等价性
在非相对论极限下,S6z(六重旋转反演)与 C6z(六重旋转)在决定交变磁性特征时产生完全等价的效果,都导致四个节面(一个基面 + 三个对角节面)的出现。这一结论已通过严格的对称性分析得到数学证明。
- 碎片化节线曲线(FNCs)的首次发现
当对称性从 C6z/S6z 降低到 C2z 时,三个对角节面消失,取而代之的是遍布整个布里渊区的碎片化节线曲线(FNCs)。沿 FNC,自旋相反的子晶格能带对保持简并;在其他位置,它们以不同的强度发生劈裂。FNC 的位置是能带依赖的,每一对能带都有独特的 FNC 形状。这是此前文献中从未报道过的新奇量子现象。RbMnPO4 的计算结果进一步证实了 FNCs 是 C2 旋转对称性下的普遍现象。
- 应变工程的可行性验证
通过对原始 CrSb 施加沿基面晶轴的单轴应变(5%),可以实际实现对称性从 C6z 到 C2z 的降低,从而在真实材料中产生 FNCs。
- AHC 的灵活可调性
对称性降低使得 MS-V 和应变 CrSb 在 Neel 矢量沿面外和面内两种取向下都能产生非零 AHC。特别是面外取向下的非零 AHC 突破了原始 CrSb 的限制,为交变磁体在自旋电子学器件中的应用提供了新的自由度。应变 CrSb 在费米能级附近的 AHC 值甚至高于原始 CrSb。