来源论文: https://arxiv.org/abs/2604.03531v1 生成时间: Apr 07, 2026 15:35

Kagome 格点上的原发性对密度波（PDW）：双轨道 Hubbard 模型与奇异模式重整化群的深度解析

0. 执行摘要

对密度波（Pair Density Wave, PDW）是凝聚态物理中一种极具异域色彩的超导态，其核心特征是 Cooper 对在零磁场下具有非零的总中心动量（CMM）。与传统的超导态（动量为零）或由外场诱导的 FFLO 态不同，原发性 PDW 态作为基态的实现极其罕见。本文基于最新的研究成果，解析了在 Kagome 格点双轨道 Hubbard 模型中，如何通过奇异模式函数重整化群（SM-FRG）方法发现稳健的原发性 PDW 相。

研究表明，Kagome 格点独特的几何结构所带来的强子格点极化（Sublattice Polarization）与多轨道电子结构的轨道极化（Orbital Polarization）相互耦合，有效地抑制了零动量的 Cooper 对形成（受限于原位库仑排斥 $U$），转而支持不同子格点间在 Brillouin 区边界 $M$ 点的配对。这种 PDW 态不仅是基态，还展现了丰富的手性拓扑特性，并与自旋密度波（SDW）和电荷密度波（CDW）呈现出复杂的交织关系。这一发现为 $CsCr_3Sb_5$ 等新型 Kagome 材料的实验观测提供了坚实的理论支撑。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：什么是“真正的”原发性 PDW？

在超导研究中，PDW 通常分为两种：

  1. 次级 PDW（Secondary PDW）：由于系统中已经存在电荷密度波（CDW）等空间调制背景，超导序参数被迫在其上产生波纹。这种情况下，PDW 不是物理驱动力，而是被动的响应。
  2. 原发性 PDW（Primary PDW）：系统自发打破平移对称性，Cooper 对直接以非零动量 $Q_P$ 凝聚。这是电子关联效应直接导致的奇异量子态。

本研究的核心问题在于：能否在一个具体且现实的微观模型中，在不依赖外场或预设密度波背景的前提下，获得作为基态的原发性 PDW？

1.2 理论基础：Kagome 格点与双轨道 Hubbard 模型

Kagome 格点（由交错三角形组成的六角格点）因其几何受挫特性和特殊的能带结构（包括平带、Dirac 点和范霍夫奇异点 vHS）而备受关注。本工作构建了一个双轨道 Hubbard 模型，其动能部分由 Slater-Koster 规则定义：

$$H_0 = \sum_{\langle ij \rangle abs} (t_{ij}^{ab} c_{ias}^\dagger c_{jbs} + h.c.) + \sum_{ias} (E_a - \mu) n_{ias}$$

其中 $a, b$ 代表局部轨道 $x, y$。相互作用部分包括标准的多轨道 Hubbard 项：

$$H_I = U \sum_{ia} n_{ia\uparrow} n_{ia\downarrow} + U' \sum_{i, a满足对称性要求 $U = U' + 2J_H$。该模型的核心物理在于：电子在 Fermi 面上的分布具有高度的轨道和子格点选择性（Sublattice selectivity）。

1.3 技术难点：为什么 PDW 在理论上难以获得？

实现 PDW 面临两大严峻挑战：

  1. 动量简并挑战：在标准 BCS 理论中，时间反演对称的两个状态 $|k, \uparrow\rangle$ 和 $|-k, \downarrow\rangle$ 具有完全相同的能量，这导致零动量配对在无穷小吸引作用下就会发生红外发散。而非零动量配对（$k$ 与 $-k+Q$）通常不占据简并能级，其极化率（Susceptibility）不发散，因此在竞争中天生处于劣势。
  2. 相互作用竞争挑战：粒子-空穴（PH）通道的交互作用通常倾向于产生 CDW 或 SDW，而这些序往往会先于 PDW 出现并占据主导地位。

1.4 方法细节：奇异模式函数重整化群（SM-FRG）

为了处理这种强关联竞争，作者采用了 SM-FRG。其核心思想是监控有效相互作用顶点 $\Gamma$ 随能量尺度 $\Lambda$（红外截止）降低时的演化过程。顶点在配对（P）、交叉（C）和直接（D）通道中进行分解：

  • 顶点演化：通过求解 $\partial \Gamma / \partial \Lambda$ 的演化方程，考虑粒子-粒子（PP）和粒子-空穴（PH）通道的二阶贡献（Loop diagrams）。
  • 奇异值分解（SVD）：在动量空间对散射矩阵进行分解：$V_{mn}(q) = \sum_{\alpha} \phi_{\alpha m}(q) S_{\alpha}(q) \phi_{\alpha n}^*(q)$。其中 $S_{\alpha}$ 是特征值，表示关联强度；$\phi_{\alpha}$ 是特征模式，决定了序参数的结构。
  • 自洽关联：SM-FRG 能够无偏见地处理所有通道的竞争，其发散点 $\Lambda_c$ 标志着系统进入了特定的有序态。

2. 关键 Benchmark 体系、计算所得数据与性能数据

2.1 体系参数设置

研究选取了能够模拟 $CsCr_3Sb_5$ 费米面拓扑的参数：

  • 斯莱特-科斯特参数：$t_{pp\sigma} = 5/8, t_{pp\pi} = -1/8$。
  • 晶体场能：$E_x = -0.055, E_y = 1.182$。
  • 化学势 $\mu = 1.36$，使得 Fermi 面避开了范霍夫奇异点（vHS），从而避免了平凡的 CDW 占优。

2.2 Fermi 面特征与子格点极化

  • $\\Gamma$ 洞穴袋：主要由变形的平带形成，子格点分量在动量空间中分布较为分散。
  • $M$ 电子口袋：位于 Brillouin 区边界，展现出极强的子格点纯度。例如，在 $M_1$ 附近的电子几乎完全分布在子格点 1 上。
  • 物理推论：由于 $|k \uparrow\rangle$ 和 $|-k \downarrow\rangle$ 倾向于分布在同一个子格点上，原位排斥 $U$ 会剧烈抑制零动量的传统超导。而对于动量为 $Q=M$ 的配对，它可以涉及不同子格点之间的电子，从而规避原位 $U$ 的惩罚。

2.3 关键数据分析

  1. 极化率分析（Fig. 2a）：

    • $\chi_{ph}$ 在 $M$ 点附近有峰值，但由于子格点选择性被显著削弱。
    • $\chi_{pp}$ 在 $\Gamma$ 点最强（BCS 不稳定性），但在考虑了原位相互作用后，该模式被抑制。而在 $M$ 点，由于涉及 NN（最近邻）键上的不同子格点配对，展现出了宽广的增强区域。

  2. FRG 流演化（Fig. 2b）：

    • 设置 $U=3.0, J_H=0.35$。随着 $\Lambda$ 降低，SDW 序最先增长，但超导（SC）序在低能尺度突然发生交叉并最快发散。
    • 发散特征值的动量分布在 $M_1, M_2, M_3$ 处呈现尖锐峰值，直接证明了 PDW 是主导序。

  3. 能隙（Gap）函数与单粒子谱（Fig. 4）：

    • 在 PDW 相中，费米面上的准粒子能隙展现出高度各向异性。$M$ 口袋处的能隙较大，而 $\Gamma$ 口袋处的能隙较小。
    • 态密度（DOS）分析显示，PDW 态仅部分抑制了费米面处的低能激发，留下了类似“残留费米面”的特征，这与 $KV_3Sb_5$ 等材料的实验观测一致。

3. 代码实现细节与复现指南

3.1 算法流程复现

复现该研究的核心是编写一个多轨道的 SM-FRG 求解器：

  1. 格点初始化：构建 Kagome 格点的实空间坐标（3 子格点单元），定义局部轨道基组。
  2. 能带对角化：通过 Slater-Koster 矩阵计算 Bloch 波函数和能带能量 $\epsilon_{nk}$。建立从轨道基组到能带基组的变换矩阵 $U_{o,n}(k)$。
  3. 顶点离散化：将动量空间划分为足够精细的网格（通常为 $24 \times 24$ 或更高）。由于是多轨道体系，顶点 $\Gamma$ 是一个张量。
  4. 求解演化方程：使用多级 Runge-Kutta 算法迭代 $\Lambda$。在每一步，计算 PP 通道极化率 $\chi_{pp}$ 和 PH 通道极化率 $\chi_{ph}$。
  5. 奇异值截断：为了提升性能，只保留散射矩阵中绝对值最大的前若干个奇异值模式（即“奇异模式”名称的来源）。对于本模型，通常保留前 100-200 个模式即可捕获主要物理。

3.2 软件包建议

目前学术界常用的开源或半开源 FRG 框架包括：

  • Pymatgen/Quantumespresso：用于前期的能带结构参数提取。
  • TPRIV (Functional Renormalization Group for Interacting Vertex)：一些研究组开发的私有 C++/Fortran 库。推荐参考南京大学或维尔茨堡大学公开的 FRG 实现思路。
  • Julia 实现：考虑到多轨道的大型矩阵运算，Julia 语言在处理 SM-FRG 的高维张量收缩时具有极佳的性能优势。

3.3 计算成本

  • 内存要求：由于需要存储 $\Gamma(q, k_1, k_2)$，且每个动量点都有轨道/子格点索引，内存消耗随网格密度和轨道数呈 $O(N_k^2 N_{orb}^4)$ 增长。本研究的双轨道模型在标准服务器（128GB-256GB RAM）上可运行。
  • 并行化：动量点的 Loop 是高度可并行的，通常使用 MPI 或 OpenMP 优化。

4. 关键引用文献与局限性评论

4.1 关键引用文献

  1. [1, 2] FFLO 原型：Fulde, Ferrell, Larkin, Ovchinnikov 的早期工作，确立了非零动量配对的概念。
  2. [12, 13] TMDCs 中的 PDW：在过渡金属硫族化合物中观测到二次 PDW 的实验文献。
  3. [26, 27] Kagome 实验证据：近期在 Kagome 超导体中观测到 PDW 信号的 Nature 文献，是本理论工作的直接诱因。
  4. [53] $CsCr_3Sb_5$：作为本文模型最直接的现实候选材料。
  5. [54, 55] SM-FRG 理论框架：汪强华教授团队关于奇异模式 FRG 方法论的基础文献。

4.2 局限性评论

  1. 重整化群截断：SM-FRG 忽略了三粒子交互作用（6 费米子项）以及更高阶的自能修正。在极强关联区域，这些项可能影响结果的定量准确性。
  2. 频率依赖性：本文采用了静态顶点（Static vertex）近似。虽然这在处理费米面不稳定性时是标准做法，但对于动态关联效应较强的体系，频率相关的 FRG 会更完备。
  3. 平均场后续处理：在 FRG 顶点发散后，为了计算 Gap 和谱函数，作者使用了平均场理论。这意味着 PDW 态本身的涨落效应未被完全包含在最终的谱分析中。

5. 补充：PDW 的拓扑性与交织序（Intertwined Orders）

5.1 手性 PDW 与拓扑

由于 Kagome 格点具有 $C_6$ 对称性，$M_1, M_2, M_3$ 三个动量点是简并的。研究发现，这些 PDW 组件可以线性组合形成一种手性 PDW（Chiral PDW）：

$$\Delta(r) = \Delta_{M1} e^{i M_1 \cdot r} + \Delta_{M2} e^{i M_2 \cdot r + i\theta} + \Delta_{M3} e^{i M_3 \cdot r + i2\theta}$$

这种态自发打破了时间反演对称性，并可能具有非零的 Chern 数，展现出拓扑超导的特征。这为实验上观测自发环流（Loop currents）提供了理论依据。

5.2 PDW 与 SDW/CDW 的量子交织

一个迷人的物理细节是：虽然 PDW 是主导序，但它会诱导出关联的密度波：

  • 诱导 CDW：由 PDW 序参数的复合项产生 $\rho_{2M} \sim \Delta_M^* \Delta_{-M}$。在本文中，由于原位 $U$ 对电荷密度的抑制，诱导出的 CDW 以**键流（Bond Current）**的形式存在，而非单纯的点电荷密度调制。
  • SDW 的角色：SDW 在高能尺度为 PDW 的形成提供了自旋涨落媒介。这种“交织”意味着 PDW 材料在进入超导态之前，往往已经展现出极强的自旋/电荷涨落特征。

5.3 对冷原子实验的启示

除了固体材料，该模型完全可以在光学格点中的多组分超冷原子系统中实现。通过调节 Feshbach 共振改变相互作用强度，或者通过格点辅助隧道效应调节跳符参数，可以精确复现本文预测的相图边界。

结语：这项工作标志着我们在理解非零动量超导这一复杂量子现象上迈出了关键一步。通过挖掘格点几何与轨道自由度的内在关联，证明了 PDW 不仅仅是实验室中的孤例，而是强关联物理中一种自然且稳健的基态选择。