Inherent momentum-dependent gap structure of altermagnetic superconductors

基本信息

• 论文ID: 2509.03247
• 标题: Inherent momentum-dependent gap structure of altermagnetic superconductors
• 作者: Christian L. H. Rasmussen, Jannik Gondolf, Mats Barkman, Mercè Roig, Daniel F. Agterberg, Andreas Kreisel, Brian M. Andersen
• 机构: Niels Bohr Institute (Copenhagen), University of Wisconsin–Milwaukee
• 分类: cond-mat.supr-con (凝聚态物理-超导)
• 提交时间: 2025年11月14日
• 论文链接: https://arxiv.org/abs/2509.03247v2

摘要

本文研究了交替磁体金属(altermagnetic metals)中的超导态特性。交替磁体打破时间反演对称性，通过补偿的Néel有序共线磁矩产生自旋劈裂的费米面。研究基于包含亚晶格自由度的微观模型，发现亚晶格结构强烈影响超导能隙结构：对于动量无关的吸引相互作用，在布里渊区边界产生能隙节点；而扩展范围相互作用可稳定自旋单态和等自旋配对三重态。在大自旋劈裂极限下，等自旋配对三重态超导占主导，并表现出交替磁序导致的非幺正特性。

研究背景与动机

研究问题

交替磁体(altermagnets)作为继铁磁体和反铁磁体之外的第三类共线磁性材料，其独特之处在于：

1. 对称性特征：缺乏时间反演与反演(或平移)的联合对称性，但保持时间反演与旋转的联合对称性
2. 电子结构：即使无净磁化和自旋轨道耦合，也展现大的动量依赖自旋劈裂能带
3. 超导可能性：作为金属，低温下可能发生超导转变

问题重要性

1. 新物理平台：交替磁性超导体为研究磁性与超导的相互作用提供新平台
2. 自旋电子学应用：自旋劈裂费米面可用于自旋电子器件
3. 非常规超导：可能稳定FFLO态、拓扑超导、超导二极管效应等奇异现象

现有方法局限

现有理论研究主要采用简化单带模型：

• 在布里渊区中心展开的低能模型
• 施加动量依赖的Zeeman自旋劈裂场
• 忽略亚晶格自由度：但实际交替磁体材料必然包含多个磁性格点/单胞

实际交替磁体具有q=0的Néel序，不扩展晶胞，必须包含亚晶格结构。

研究动机

本文基于包含亚晶格自由度的最小模型，系统研究：

1. 亚晶格结构对超导配对的影响
2. 不同相互作用(现场吸引 vs 扩展相互作用)下的能隙结构
3. 等自旋配对三重态的非幺正性质

核心贡献

1. 建立完整理论框架：基于包含亚晶格自由度的最小模型，发展了交替磁性超导的自洽平均场理论
2. 发现固有能隙节点机制：对于动量无关的现场吸引相互作用(如声子介导)，亚晶格结构在布里渊区边界产生阻塞效应(blocking effect)，导致能隙节点
3. 揭示扩展相互作用下的丰富相图：
   • 可稳定自旋单态(d波、扩展s波)和等自旋配对三重态
   • 大自旋劈裂极限下，等自旋配对三重态占优
4. 发现非幺正三重态性质：等自旋配对态展现特征性非幺正分量 d×d* ≠ 0，直接反映交替磁序的自旋劈裂
5. 预言双相变：向完全能隙的等自旋配对态转变通过两个热力学相变实现

方法详解

任务定义

研究具有亚晶格结构的二维交替磁体金属的超导配对问题：

• 输入：交替磁体正常态哈密顿量(包含亚晶格A、B)、吸引相互作用类型
• 输出：自洽超导能隙结构Δ(k)、配对对称性、相图
• 约束：限制在q=0均匀超导态(主导相图)

模型架构

1. 交替磁体正常态模型

采用中心对称空间群的最小模型，磁性原子占据多重度为2的反演对称Wyckoff位置。二维方格子情况：

$H = \sum_k c_k^\dagger H_{MM} c_k$

其中 $c_k^\dagger = (c_{↑A,k}^\dagger, c_{↑B,k}^\dagger, c_{↓A,k}^\dagger, c_{↓B,k}^\dagger)$，哈密顿量矩阵：

$H_{MM} = \varepsilon_{0,k}\tau_0 + t_{x,k}\tau_x + t_{z,k}\tau_z + \tau_y\lambda_k \cdot \sigma + \tau_z N \cdot \sigma$

关键参数：

• $\tau_i$, $\sigma_i$：亚晶格和自旋Pauli矩阵
• $t_{x,k}$：亚晶格间跳跃
• $t_{z,k}$：亚晶格内各向异性跳跃(驱动自旋劈裂)
• $N$：Néel序参量
• $\lambda_k$：自旋轨道耦合(本文忽略)

2. 不同类型交替磁体

根据$t_{z,k}$的形式区分(见表I)：

本文主要研究 $t_{z,k} = 4t_4\sin k_x \sin k_y$ (dxy型)。

3. 亚晶格-能带基变换

通过规范变换消除相位，对角化得到能带本征态：

$E_{\sigma,k}^\lambda = \varepsilon_{0,k} \pm \sqrt{|t_{x,k}|^2 + (t_{z,k} + \sigma N)^2}$

关键的亚晶格权重函数：

$l_{\sigma,k} = \cos\frac{\theta_{\sigma,k}}{2}, \quad m_{\sigma,k} = \frac{t_x}{|t_x|}\sin\frac{\theta_{\sigma,k}}{2}$

其中： $\cos\frac{\theta_{\sigma,k}}{2} = \frac{1}{\sqrt{2}}\sqrt{1 + \frac{t_{z,k} + \sigma N}{\sqrt{|t_{x,k}|^2 + (t_{z,k} + \sigma N)^2}}}$

这些权重函数决定了超导配对的动量依赖性。

超导相互作用与能隙方程

1. 现场(OS)吸引

$H_{int}^{OS} = -\frac{V}{N}\sum_{s,k,k'} c_{↑,s,k}^\dagger c_{↓,s,-k}^\dagger c_{↓,s,-k'} c_{↑,s,k'}$

自洽能隙方程(α能带)：

$\Delta_{↑↓,k}^\alpha = -\frac{V}{N}\left[l_{↑,k}l_{↓,-k}\sum_{k'} l_{↑,k'}l_{↓,-k'}\frac{\Delta_{↑↓,k'}^{\alpha*}}{2E_{↑↓,k'}^{\alpha,0}}(1-2f(E_{↑↓,k'}^{\alpha,+})) + ...\right]$

关键特征：能隙通过 $l_{\sigma,k}l_{\sigma',-k}$ 和 $m_{\sigma,k}m_{\sigma',-k}$ 组合获得内禀动量依赖性。

2. 最近邻(NN)吸引

连接相反亚晶格，相互作用：

$V_{k,k',A,B}^{NN} = 4V^{NN}e^{-i(k_x+k_y)/2}e^{i(k_x'+k_y')/2}\cos\frac{k_x-k_x'}{2}\cos\frac{k_y-k_y'}{2}$

能隙方程包含不同的本征矢量组合 $l_{\sigma,k}m_{\sigma',-k}$。

3. 次近邻(NNN)吸引

连接相同亚晶格：

$V_{k,k',s,s}^{NNN} = 2V^{NNN}[\cos(k_x-k_x') + \cos(k_y-k_y')]$

技术创新点

1. 完整保留亚晶格自由度：不同于简化单带模型，保留了 $l_{\sigma,k}$ 和 $m_{\sigma,k}$ 的完整动量依赖
2. Fourier分解分析：将本征矢量组合分解为： $l_{↑,k}l_{↓,-k} = \sum_{\nu,\eta} c_{\nu,\eta}^l e^{-ik_x\nu}e^{-ik_y\eta}$ 揭示了配对的有效范围(图3)
3. 阻塞效应机制：在布里渊区边界，非公度对称性强制能带简并，导致：
   • 自旋和亚晶格权重纯化
   • 相同亚晶格相反自旋态重叠为零
   • 阻断现场吸引的自旋单态配对
4. 非幺正三重态参数化： $\Delta_{↑↑}^\pm \propto p_{x-y} \pm i\alpha p_{x+y}, \quad \Delta_{↓↓}^\pm \propto p_{x+y} \pm i\alpha p_{x-y}$ 非幺正分量：$d \times d^* \propto (1-\alpha^2)p_{x-y}p_{x+y}$

实验设置

模型参数

• 跳跃参数：$\{t_1, t_2, t_3, t_4\} = \{0.425, 0.05, -0.025, -0.075\}t$ (t为能量单位)
• Néel序参量：$N = 0.2t$ (贯穿全文)
• 化学势：研究三个代表性费米面
  • μ = 0.3t (小费米面)
  • μ = 0.6t (中等费米面，穿越布里渊区边界)
  • μ = 1.0t (大费米面)

相互作用强度

• 现场吸引：$V = 1.5t$
• 最近邻：$V^{NN} = 0.5-0.9t$
• 次近邻：$V^{NNN} = 1.0t$

计算方法

1. 自洽求解：迭代Bogoliubov-de Gennes方程直至收敛
2. 温度：零温极限，使用Fermi-Dirac分布
3. 动量网格：布里渊区精细网格
4. 初始条件：随机复数初值(用于探测简并态)

评价指标

• 超导能隙振幅 $|\Delta(k)|$ 及其动量分布
• 能隙节点位置
• d矢量分量 $(d_x, d_y)$ 和非幺正性 $d \times d^*$
• 配对对称性(s波、d波、p波)
• 相变温度(通过能隙幅度判断)

实验结果

主要结果1：现场吸引下的动量依赖能隙

能隙结构(图2)

对于动量无关的现场吸引，能隙展现显著动量依赖：

• 小费米面(μ=0.3t)：能隙相对均匀，$\Delta_{max} = 0.17t$
• 中等费米面(μ=0.6t)：穿越布里渊区边界，边界出现能隙节点，$\Delta_{max} = 0.24t$
• 大费米面(μ=1.0t)：在X、Y高对称点附近能隙显著抑制，$\Delta_{max} = 0.18t$

阻塞效应机制(图3)

沿Γ-X-M-Γ路径的分析揭示：

1. 布里渊区边界：
   • 自旋和亚晶格权重完全分离
   • $|d_{α,k,↑}|^2$ 在A亚晶格为1，B亚晶格为0
   • $|d_{α,k,↓}|^2$ 在B亚晶格为1，A亚晶格为0
2. Fourier分解：
   • 现场分量 $c_{00}^l = c_{00}^m = 0.414$ 最大
   • 但存在显著的扩展分量(最近邻0.05，次近邻0.02)
3. 物理机制：
   • 非公度对称性在边界强制 $t_{x,k}$ 和 $t_{z,k}$ 同时为零
   • 仅Néel序打破简并
   • 相同亚晶格相反自旋态无重叠 → 阻断配对

依赖性

• Néel序强度：N增大使阻塞区域向布里渊区中心扩展
• 交替磁体类型：
  • 非公度型(dxy)：整个边界阻塞
  • 公度型(dx²-y²)：仅M点最强阻塞

主要结果2：扩展相互作用下的配对对称性

最近邻吸引(图4)

在弱到中等自旋劈裂区，稳定自旋单态：

关键特征：

• 边界无节点：NN可配对不同亚晶格，绕过阻塞效应
• 动量依赖主要来自裸相互作用 $V_{k,k'}^{NN}$
• 费米面形状决定优势对称性

主要结果3：等自旋配对三重态(图5)

四重简并态

在大自旋劈裂极限($V^{NNN}=1t$, $N=0.2t$, $μ=0.6t$)，发现四个简并的等自旋配对态：

$\Delta_{↑↑}^+ \otimes \Delta_{↓↓}^-, \quad \Delta_{↑↑}^- \otimes \Delta_{↓↓}^+, \quad \Delta_{↑↑}^+ \otimes \Delta_{↓↓}^+, \quad \Delta_{↑↑}^- \otimes \Delta_{↓↓}^-$

分类：

• 前两个：手性态(保持正常态联合对称性)
• 后两个：螺旋态(破坏联合对称性)

非幺正性质

所有四态展现相同的非幺正分量结构(图5c)：

$i(d \times d^*)_z \propto (1-\alpha^2)p_{x-y}p_{x+y}$

特征：

• 振幅 $\sim 10^{-2}$ (相对于 $|d| \sim 0.2$)
• 空间分布反映自旋劈裂的C4对称破缺
• 符号由 $N_z t_z$ 乘积决定(磁序方向)

d矢量结构(图5b)

• 实部：主要在 $d_x$ 或 $d_y$，展现p波节点结构
• 虚部：相互正交，幅度由 $\alpha \neq 1$ 调制
• 四态具有相同 $|d|$ 和谱隙

双相变

• 高温：向第一个p波不可约表示转变
• 低温：向第二个不可约表示转变(若自旋扇区解耦)
• 耦合情况：低温转变变为交叉

消融分析

通过改变参数验证机制：

1. 去除亚晶格结构 ($t_z = 0$, $N = 0$)：
   • 现场吸引恢复均匀能隙
   • 证实动量依赖源于亚晶格
2. 改变Néel序强度：
   • $N = 0.1t$：阻塞区域缩小
   • $N = 0.4t$：阻塞区域扩大至更内部
   • 证实阻塞效应可调
3. 不同交替磁体类型：
   • dx²-y² 型：节点仅在M点最强
   • IV型(无自旋劈裂)：类似阻塞效应
   • 证实普适性

相关工作

交替磁性超导理论

1. FFLO态研究 15-24：
   • 自旋劈裂费米面不利于均匀自旋单态
   • 导致有限动量Cooper配对
   • Bogoliubov费米面和超导二极管效应
2. 非常规超导 18,29-36：
   • 等自旋配对p波三重态优势
   • 拓扑超导可能性
   • 本文贡献：揭示非幺正性质和双相变
3. 简化模型 18,36-38：
   • 单带模型 + 动量依赖Zeeman场
   • 忽略亚晶格 → 遗漏阻塞效应和能隙节点
   • 本文改进：完整亚晶格模型

亚晶格效应研究

4. 最近工作 23,49-51：
   • Ref.49：配对结构约束
   • Ref.23：FF相稳定
   • Ref.50：杂质束缚态
   • Ref.51：涨落介导配对
   • 本文互补：系统研究能隙结构

实验进展

5. 候选材料 57-63：
   • RuO2薄膜(应变下超导)
   • KV2Se2O, CrSb, MnTe, CoNb4Sb8
   • 本文预言可实验验证

结论与讨论

主要结论

1. 亚晶格驱动的能隙各向异性：
   • 现场吸引在布里渊区边界产生固有节点
   • 非公度交替磁体中效应最显著
   • 预期在实际材料中普遍存在(费米面常穿越边界)
2. 扩展相互作用的丰富相图：
   • 小自旋劈裂：自旋单态(s波、d波)
   • 大自旋劈裂：等自旋配对三重态
   • 费米面形状决定对称性选择
3. 非幺正三重态的交替磁指纹：
   • $d \times d^* \neq 0$ 直接反映自旋劈裂
   • 四重简并态(手性/螺旋)
   • 双相变机制
4. 普适性：结论适用于反铁磁金属超导体

局限性

1. 模型简化：
   • 限制在q=0均匀超导(忽略FFLO态)
   • 忽略自旋轨道耦合
   • 二维模型(实际材料三维)
2. 相互作用处理：
   • 平均场近似(忽略涨落)
   • 限制带内配对(弱耦合)
   • 唯象吸引势(未考虑微观机制)
3. 参数选择：
   • 固定Néel序强度N=0.2t
   • 特定跳跃参数组合
   • 零温计算(未研究有限温度)
4. 实验对应：
   • 未与具体材料定量比较
   • 杂质和无序效应未考虑

未来方向

1. 材料预言：
   • RuO2薄膜能隙结构的STM测量
   • CrSb等材料的超导配对对称性
   • 非幺正三重态的实验签名
2. 理论扩展：
   • 包含FFLO有限动量配对
   • 自旋轨道耦合效应
   • 三维模型和层间耦合
3. 拓扑性质：
   • 手性p波态的边缘态
   • Majorana零模
   • 拓扑相变
4. 动力学：
   • 超导涨落效应
   • 集体模式
   • 非平衡动力学

深度评价

优点

1. 方法创新性 ⭐⭐⭐⭐⭐

• 首次系统研究：完整包含亚晶格自由度的交替磁性超导理论
• 新物理机制：发现阻塞效应导致的固有能隙节点
• 技术严谨：从最小模型出发，自洽求解BdG方程
• 分析深入：Fourier分解揭示配对的有效范围

2. 物理洞察 ⭐⭐⭐⭐⭐

• 对称性分析：清晰阐明非公度对称性的关键作用
• 机制透明：图3直观展示边界处亚晶格/自旋权重分离
• 普适结论：适用于交替磁体和反铁磁金属
• 预言可检验：STM可直接观测能隙节点

3. 实验相关性 ⭐⭐⭐⭐

• 及时性：RuO2等材料刚实现超导
• 可操作性：
  • 能隙各向异性：角分辨光电子谱(ARPES)
  • 节点位置：热导率的角度依赖
  • 非幺正性：核磁共振Knight shift
• 材料指导：指出费米面位置的重要性

4. 理论完整性 ⭐⭐⭐⭐

• 多场景覆盖：现场/NN/NNN相互作用
• 不同区间：弱/强自旋劈裂极限
• 配对类型：自旋单态和三重态
• 数值充分：三个代表性费米面

不足

1. 模型局限 ⭐⭐⭐

• 维度限制：仅考虑二维，实际材料三维
• 相互作用简化：唯象吸引势，未考虑电子-声子/自旋涨落微观机制
• 参数固定：N=0.2t贯穿全文，缺乏N依赖的系统研究
• SOC忽略：实际材料中SOC可能重要

2. 定量比较缺失 ⭐⭐⭐

• 无材料对应：未与RuO2等具体材料参数定量比较
• 能隙尺度：$\Delta \sim 0.2t$ 对应实际温度未明确
• 相互作用强度：V的微观起源和合理范围未讨论
• 实验可观测量：缺乏具体数值预言(如Tc、节点角度)

3. 理论深度 ⭐⭐⭐

• 平均场近似：忽略配对涨落(可能在低维重要)
• 弱耦合假设：限制带内配对，强耦合效应未探讨
• 拓扑性质浅尝：仅提及手性态拓扑性，未深入分析
• 稳定性分析：缺乏对不同态的自由能比较

4. 呈现方式 ⭐⭐⭐⭐

• 图5信息密集：10个子图难以快速理解
• 数学细节：方程(16)-(28)冗长，可能放附录
• 物理图像：缺乏阻塞效应的实空间示意图
• 结论强度：部分陈述("expect on general grounds")缺乏定量支撑

影响力评估

学术影响 ⭐⭐⭐⭐⭐

• 开创性：建立交替磁性超导的标准理论框架
• 引用潜力：后续研究必引(亚晶格效应不可忽略)
• 方法示范：为其他磁性超导体研究提供模板
• 争议激发：非幺正三重态预言将引发实验竞赛

实用价值 ⭐⭐⭐⭐

• 材料设计：指导寻找具有特定能隙结构的超导体
• 器件应用：
  • 能隙节点 → Josephson结的相位敏感探测
  • 非幺正态 → 自旋流产生
  • 手性态 → 拓扑量子计算
• 实验指导：明确需要测量的物理量

可复现性 ⭐⭐⭐⭐

• 方法清晰：BdG方程和参数完全给出
• 代码可实现：基于标准数值方法
• 但缺乏：
  • 收敛判据未明确说明
  • 动量网格大小未给出
  • 初始条件选择策略未详述

适用场景

直接适用 ✅

1. 交替磁体超导体：RuO2, KV2Se2O, CrSb, MnTe
2. 反铁磁金属超导体：CeCoIn5, UPd2Al3等重费米子
3. 理论建模：任何具有亚晶格的磁性超导体

需要修正 ⚠️

1. 三维材料：需扩展到3D模型(层间耦合)
2. 强SOC系统：需包含λk项
3. 强关联材料：需超越平均场

不适用 ❌

1. 非中心对称超导体：对称性不同
2. 铁磁超导体：净磁化破坏配对
3. 一维系统：Luttinger液体物理主导

创新性评分

关键技术细节补充

阻塞效应的数学表述

在布里渊区边界(如X点 kx=π, ky=0)： $t_{x,k} = 0, \quad t_{z,k} = 0 \quad \Rightarrow \quad \theta_{\sigma,k} = \begin{cases} 0 & \sigma N > 0 \\ \pi & \sigma N < 0 \end{cases}$

因此： $l_{↑,k} = 1, m_{↑,k} = 0, \quad l_{↓,k} = 0, m_{↓,k} = 1$

导致： $l_{↑,k}l_{↓,-k} + m_{↑,k}m_{↓,-k} = 0 \quad \Rightarrow \quad \Delta_{OS}(k) = 0$

非幺正分量的物理意义

$\mathbf{d} \times \mathbf{d}^* = 2i\text{Im}(d_x d_y^*)\hat{z} \propto (1-\alpha^2)p_{x-y}p_{x+y}\hat{z}$

物理后果：

1. 自旋极化：$\langle S_z \rangle \propto d \times d^* \neq 0$
2. 破坏自旋旋转对称：即使总自旋S=1
3. 实验签名：NMR Knight shift不为零

相变温度估计

从BCS理论： $k_B T_c \sim 1.14\omega_D e^{-1/(N(0)V)}$

对于 $V=1.5t$，$N(0) \sim 0.1/t$： $T_c \sim 0.1t \sim 100\text{K} \quad (\text{若}\ t=1\text{eV})$

但实际Tc受自旋劈裂和配对竞争压制。

参考文献（精选）

6 L. Šmejkal et al., Phys. Rev. X 12, 040501 (2022) - 交替磁体综述
18 S. Hong et al., Phys. Rev. B 111, 054501 (2025) - 单带模型p波配对
47 M. Roig et al., Phys. Rev. B 110, 144412 (2024) - 最小模型
49 D. Chakraborty & A. M. Black-Schaffer, Phys. Rev. B 112, 014516 (2025) - 配对约束
57 M. Uchida et al., Phys. Rev. Lett. 125, 147001 (2020) - RuO2超导

总评：这是一篇高质量的理论物理论文，在交替磁性超导这一新兴领域做出了开创性贡献。通过包含亚晶格自由度的完整处理，发现了现有简化模型遗漏的重要物理(阻塞效应、非幺正三重态)。方法严谨，物理图像清晰，预言可实验检验。主要不足是缺乏与具体材料的定量对应和强耦合效应的探讨。预期将成为该领域的重要参考文献，并激发大量后续实验和理论工作。推荐发表于高水平期刊(如PRB或PRX)。