来源论文: https://arxiv.org/abs/2604.08937v1 生成时间: Apr 13, 2026 06:57

0. 执行摘要

在强关联电子系统研究中，寻找和理解掺杂 Mott 绝缘体中的非常规超导性（SC）是凝聚态物理的核心挑战。本文聚焦于蜂巢晶格（Honeycomb Lattice）上的扩展 $t$-$J$ 模型，通过大尺度密度矩阵重整化群（DMRG）模拟和从属玻色子平均场理论（SBMFT）的结合，系统性地研究了次近邻（NNN）电子跃迁 $t'$ 和相应的超交换相互作用 $J'$ 对系统基态的影响。研究发现，在特定的 YC 型柱状几何结构下，系统表现出稳健的准长程 $d$ 波超导序，并伴随着亚领先的扶手椅型条纹（a-stripe）电荷序。特别地，超导 Luttinger 指数 $K_{sc}$ 随 $t'$ 呈现非单调的“谷状”依赖关系，在 $t' \approx 0.4$ 时达到最优超导增强。此外，通过对比 XC 型柱状结构，揭示了边界几何在稳定竞争相中的关键作用。SBMFT 的结果进一步证实了 $a$-stripe 在二维极限下的稳定性，并在大 $t'$ 极限下预言了向均匀对等 $d$ 波超导相的演化。这项工作为理解莫尔（Moiré）系统（如扭曲双层石墨烯或 TMDs）中的超导机制提供了重要的理论支撑。

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

1.1 核心科学问题

蜂巢晶格作为石墨烯及多种类石墨烯材料的结构基础，其物理性质在掺杂 Mott 绝缘体背景下显得尤为复杂。虽然在弱耦合极限下（如随机相位近似、重整化群方法）普遍预测会出现 $d+id$ 波超导序，但在强耦合 $t$-$J$ 模型中，电荷密度波（CDW）序、条纹序（Stripe）与超导序之间的竞争关系一直存在巨大争议。本文的核心科学问题在于：次近邻跃迁 $t'$ 如何重塑蜂巢晶格上的相图？它在促进或抑制超导性方面扮演了什么角色？此外，计算模拟中常用的有限尺寸柱状几何（Cylinders）如何影响这些竞争相的稳定性？

1.2 理论基础：扩展 t-J 模型

研究采用的哈密顿量为扩展 $t$-$J$ 模型：

$$H = \mathcal{P} H_t \mathcal{P} + H_s$$

其中，$H_t = -\sum_{ij\sigma} t_{ij} (c_{i\sigma}^\dagger c_{j\sigma} + h.c.)$ 描述电子跃迁，$\mathcal{P}$ 为 Gutzwiller 投影算符，禁止电子双占据。$H_s = \sum_{ij} J_{ij} (\mathbf{S}_i \cdot \mathbf{S}_j - \frac{1}{4} n_i n_j)$ 为自旋相互作用。由于引入了次近邻效应，$t_{ij}$ 在最近邻（NN）为 $t$，在次近邻（NNN）为 $t'$。基于超交换关系，次近邻超交换作用设为 $J' = (t'/t)^2 J$。实验和理论背景设定 $t=2, J=1$，这对应于有效 Hubbard 模型中的 $U/t=8$，属于典型的强关联区。

1.3 技术难点：有限尺寸效应与纠缠熵

DMRG 虽然是处理一维和准一维系统的利器，但在处理二维晶格模型时面临严重的指数级计算复杂度提升。对于蜂巢晶格，不同的柱状几何切割方式（如 YC 和 XC 型）会对能谱和电荷分布产生显著的边界诱导偏置。此外，超导序与 CDW 序的能量尺度极其接近，需要极高的截断精度（Truncation Error）和庞大的键维（Bond Dimension $m$）才能分辨真实的基态。在本文中，为了克服这些困难，作者使用了高达 $m=24000$ 的态数，并进行了详尽的截断误差外推。

1.4 方法细节：双重手段验证

  1. DMRG 模拟：在 YC4-0、XC8-0 和 YC4-4 三种几何结构上进行。YC4-0 具有沿 $y$ 方向的周期性边界和沿 $x$ 方向的开放边界。通过测量 pair-pair 相关函数 $\Phi_{\alpha\beta}(r)$ 和局部电荷密度 $n(x,y)$，提取 Luttinger 指数 $K_{sc}$ 和 $K_c$。
  2. SBMFT (Slave-Boson Mean-Field Theory)：为了探究二维极限（$L_y \to \infty$），引入从属玻色子算符 $c_{i\sigma} = b_i f_{i\sigma}^\dagger$。在平均场近似下，将哈密顿量解耦为粒子-空穴通道（自旋子跃迁 $\chi_{ij}$）和粒子-粒子通道（对配对 $\Delta_{ij}$）。通过自洽迭代寻找能量最低的解，涵盖了 DMRG 中观察到的所有条纹图案。

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

2.1 YC4-0 柱状结构中的超导相

在 YC4-0 柱状几何上，作者发现了稳健的超导信号。关键数据点如下：

  • 超导关联函数：$\Phi_{aa}(r)$ 表现出清晰的幂律衰减 $\sim r^{-K_{sc}}$。对于掺杂浓度 $\delta = 1/12$，$K_{sc}$ 在 $t'=0$ 时为 0.88(2)，在 $t'=0.4$ 时下降至最小值 0.82(2)，随后在 $t'=0.9$ 时回升至 0.93(4)。这表明 $t'=0.4$ 附近存在一个“最佳超导区”。
  • 对对称性：测量发现 $\Phi_{aa}(r) \approx -2\Phi_{ab}(r)$，且符号发生改变，这强烈暗示了对配对具有向列型（Nematic）$d$ 波对称性，而非各向同性的 $s$ 波。
  • 电荷条纹序：系统形成了所谓的 $a$-stripe，即沿扶手椅方向均匀、沿锯齿方向呈现周期性波动的条纹。条纹波长 $\lambda$ 与掺杂浓度 $\delta$ 满足 $\lambda = \sqrt{3}/(4\delta)$。这种条纹序与超导序共存，但 $K_{sc} < 1 < K_c$ 的关系表明超导性在长程尺度上占据主导。

2.2 XC8-0 中的相竞争与相分离

相比之下，XC8-0 几何的表现完全不同，凸显了几何结构的影响：

  • 在 $t'=0$ 时，系统表现出双向 CDW 序。
  • 在 $0.2 < t' < 0.5$ 的区间内，系统进入了明显的**相分离（Phase Separation）**区，电荷分布极度不均。
  • 当 $t' > 0.5$ 时，出现了一种长程锯齿形条纹（z-stripe），但其超导关联函数表现为指数衰减 $\Phi(r) \sim e^{-r/\xi}$，相关长度 $\xi$ 仅为 0.4 个单位胞，说明超导在该几何下被完全抑制。

2.3 SBMFT 的 2D 外推性能

SBMFT 的相图（Fig. 5）显示：

  • 对于 $\delta = 1/12$，$a$-stripe 相在整个研究的 $t'$ 范围内都是能量最低的，支持了 YC4-0 结果的可靠性。
  • 对于 $\delta = 1/8$，在 $t' < 0.5$ 时 $a$-stripe 稳定，但在 $t' > 0.5$ 时系统转变为均匀的向列型 d 波超导相。这一转变点与 DMRG 中观察到的 CDW 受到强烈抑制的趋势相吻合。

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

3.1 核心算法：DMRG 与 MPS

虽然论文未直接给出开源仓库链接，但此类研究通常基于以下成熟框架：

  • ITensor (C++ / Julia)：目前处理 $t$-$J$ 模型最为高效的张量网络库。其内置的 $U(1)$ 电荷守恒和 $SU(2)$ 自旋对称性可以极大减小希尔伯特空间。
  • TeNPy (Python)：适合快速构建复杂晶格模型并进行有限键维模拟。

3.2 复现指南

若要复现 YC4-0 柱状模拟，需注意以下配置：

  1. 晶格定义：蜂巢晶格需转化为具有最近邻和次近邻连接的矩阵形式。YC4 意味着 $y$ 方向有 4 个单位胞（8 个位点），$x$ 方向取 $L_x = 48$。
  2. 算符处理：必须实施 Gutzwiller 投影。在张量网络中，这通常通过定义局部物理位点空间为 $d=3$（空穴、电子自旋向上、电子自旋向下）来实现，从而自动排除了双占据态。
  3. 收敛参数：
    • m (Bond Dimension): 初始可取 2000，逐步增加至 10000+。对于 $t$-$J$ 模型，10000 是观察超导序的门槛。
    • sweeps: 至少进行 80 次扫描，且每次扫描需配合微小的噪声（Noise）项以避免陷入局部极小值。
    • truncation_error: 目标值应小于 $10^{-6}$。

3.3 SBMFT 计算指南

SBMFT 涉及自洽非线性方程组的求解。建议使用 Python 的 Scipy 优化库：

  1. 在扩大的单位胞（基于 $a$-stripe 的周期）内定义初值 $\chi_{ij}$ 和 $\Delta_{ij}$。
  2. 构建平均场哈密顿矩阵并对角化，计算费米能级。
  3. 更新序参数并重复直至能量收敛。

4. 关键引用文献与评论

4.1 关键参考文献

  1. White (1992) [Ref 75]: DMRG 算法的奠基之作，定义了整个数值强关联物理的范式。
  2. Jiang & Devereaux (2019) [Ref 11]: 在方晶格 $t$-$J$ 模型中系统性地证明了 $t'$ 增强超导的先驱工作，本文可视为其在蜂巢晶格上的重要扩展。
  3. Cao et al. (2018) [Ref 21, 22]: 扭曲双层石墨烯（TBG）中关联绝缘相和超导性的实验发现，是本文研究的最直接物理动机。
  4. Gu et al. (2013) [Ref 47]: 早期关于蜂巢晶格 $t$-$J$ 模型的研究，本文在此基础上通过引入 $t'$ 完善了相图。

4.2 局限性评论

尽管本工作在数值和理论结合上非常出色，但仍存在以下局限：

  • 有限尺寸的诱导序：YC4 和 XC8 的结果大相径庭，说明有限的宽度对电荷分布有极强的锁定效应。虽然 SBMFT 尝试弥补这一不足，但 SBMFT 本身忽略了自旋和电荷的动力学涨落，可能低估了某些相互竞争的无序相的能量。
  • J’ 的处理：作者假设 $J' = (t'/t)^2 J$，这在 Hubbard 模型推导中是合理的，但在真实的莫尔系统中，$J'$ 的具体数值可能独立于跃迁能。如果 $J'$ 与 $t'$ 的关系偏离该假设，超导相的范围可能会大幅收缩或扩张。
  • 三维耦合缺失：真实的层状材料具有层间耦合，这对于稳定超导长程序至关重要，而本文仅讨论了二维平面。

5. 补充：物理机制深度解析

5.1 $t'$ 增强超导的物理图像

在蜂巢晶格中，$t'$ 的引入实际上改变了费米面的几何形状。当 $t' > 0$ 时，空穴的动能项得到了有效的重整化。由于 $J'$ 的引入增强了自旋子在次近邻位点的配对趋势，这种效应与最近邻的 $J$ 配对产生了协同作用。特别是在 AFM（反铁磁）背景下，适度的 $t'$ 增加了磁涨落，而这些涨落恰恰是形成 Singlet Pair 的“粘合剂”。

5.2 几何结构的敏感性解析

为什么 YC4-0 稳定了超导而 XC8-0 没有？这涉及到“条纹几何适应性”。$a$-stripe 的波长和方向在 YC 几何下可以自然地与周期性边界条件匹配，从而最小化由于边界切割导致的能量损失。而在 XC 几何中，扶手椅方向的切口迫使条纹必须以不自然的角速度发生扭曲，导致系统宁愿选择相分离或不具备超导性的 $z$-stripe。这一发现警示所有数值模拟研究者：单一种类的晶格几何所得结论可能具有误导性，必须进行多几何校验。

5.3 对未来莫尔系统实验的启示

扭曲双层石墨烯（TBG）中的有效物理可以被映射到类似的蜂巢模型上。本文预言的 $t' \approx 0.4t$ 附近存在最佳超导区，这为实验中通过栅极电压调节能带倾斜度（从而调节 $t'$）来寻找更高 $T_c$ 提供了明确的方向。此外，向列型 $d$ 波的预测也可以通过量子振荡或 STM 准粒子干扰成像进行验证。