来源论文: https://arxiv.org/abs/2604.16293v1 生成时间: Apr 20, 2026 18:05

0. 执行摘要

高温超导体的物理起源一直是凝聚态物理领域最核心的难题之一。铜氧化物（Cuprates）展现出的复杂相图——包括 $d$ 波超导（dSC）、伪能隙（Pseudogap）、电荷密度波（CDW）以及奇异金属行为——向传统的费米液体理论提出了巨大挑战。其中，对密度波（Pair Density Wave, PDW）作为一种空间调制的超导序，被认为是连接这些奇特物相的关键环节。

本研究利用 state-of-the-art 的热态张量网络（Thermal Tensor Network, ThermoTN）模拟，结合切空间张量重整化群（tanTRG）技术，对存在次近邻跳符 $t'$ 的单带 Hubbard 模型进行了大尺度有限温度数值模拟。研究发现，在空穴掺杂侧（hole-doped side），系统并未表现出稳健的 $d$ 波超导，而是在伪能隙区域的下半部分涌现出了强烈的涨落 PDW 态。这种 PDW 态具有动量 $\mathbf{Q}_{PDW} \approx (0, \pi)$ 的跨弧配对特征，与电子掺杂侧（electron-doped side）主导的零动量 $d$ 波超导形成了鲜明对比。这一发现不仅协调了长期以来数值模拟与实验观察之间的不一致性，也为理解高温超导中的配对不稳定性提供了全新的有限温度视角。

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

核心科学问题：粒子-空穴不对称性与配对竞争

在铜氧化物超导体的研究中，Hubbard 模型及其变体（如 $t-J$ 模型）被视为捕获其基本物理性质的最小模型。然而，一个长期存在的争议是：简单的单带 Hubbard 模型能否同时描述电子掺杂和空穴掺杂下的物理行为？

之前的数值研究（如 CP-AFQMC, DMRG, iPEPS）在空穴掺杂侧是否产生稳健的 $d$ 波超导上给出了相互矛盾的结果。特别是，实验上空穴掺杂的铜氧化物具有更高的超导转变温度，但在数值模拟中，其 $d$ 波超导关联往往弱于电子掺杂侧，甚至被电荷序（CDW）或条纹相（Stripes）掩盖。本项研究的核心目标是探索有限温度下的物理演化，查明在基态 CDW 之上是否存在某种前驱配对态，并理清配对关联与费米面拓扑结构之间的内在联系。

理论基础：$t-t'$ Hubbard 模型

模型哈密顿量如式 (1) 所示：

$$ H = - \sum_{ij\sigma} t_{ij} c^\dagger_{i\sigma} c_{j\sigma} + U \sum_{i} n_{i\uparrow} n_{i\downarrow} - \mu \sum_{i} n_i $$

其中，$t_{ij}$ 包含了最近邻跳符 $t=1$ 和次近邻跳符 $t' = -0.2$。$U=8$ 是典型的中等强度排斥。$t'$ 的引入至关重要，它打破了粒子-空穴对称性，使得费米面发生形变。研究采用了准二维的圆柱几何（Cylindrical Geometry），宽度为 $W=6, 8$，长度为 $L=18$。

技术难点：量子蒙特卡洛的“符号问题”与张量网络的熵增长

  1. 符号问题（Sign Problem）： 对于掺杂的 Hubbard 模型，传统的决定论量子蒙特卡洛（DQMC）在低温下会遭遇严重的指数级符号问题，导致计算无法收敛。
  2. 熵增长： 在有限温度张量网络模拟中，随着温度降低，热态密度矩阵 $\rho(\beta) = e^{-\beta H}$ 的纠缠熵（特别是算符空间纠缠熵）会迅速增加，这对张量收缩的精度和键维（Bond Dimension）提出了极高要求。

方法细节：tanTRG 与 ThermoTN

本研究采用了切空间张量重整化群（tanTRG），这是近年来张量网络领域的重要突破。其核心思想是将热态 $\rho(\beta)$ 表示为矩阵乘积算子（MPO），并利用时间演化变分原理（TDVP）在切空间中进行虚时间演化。相比于传统的线性 TRG，tanTRG 能够更有效地保留低能态信息，并支持极大的键维（本研究中达到了 $D=32768$）。

为了克服有限尺寸效应并提升动量分辨率，研究人员采用了周期性边界条件（PBC）与反周期性边界条件（APBC）的平均方案。通过虚时间代理函数（Imaginary-time proxy）来推断物理量：

  • 单粒子谱权重代理： $\beta G(\mathbf{k}, \beta/2)$ 对应于 $\omega=0$ 处的谱函数 $A(\mathbf{k}, \omega=0)$。
  • 配对结构因子代理： $\beta \Phi_{\alpha\alpha}(\mathbf{q}, \beta/2)$ 反映了零频率配对易受性，能够有效过滤高频噪声。

2. 关键 benchmark 体系，计算所得数据，性能数据

关键 Benchmark 体系：粒子-空穴不对称性的相图

研究构建了完整的温度-掺杂（$T-\delta$）相图。在电子掺杂侧（$\delta < 0$），在 $\delta \approx -1/9$ 附近发现了明显的 $d$ 波超导“穹顶”（dSC Dome）。而在空穴掺杂侧（$\delta > 0$），基态被电荷密度波（CDW）占据，但在 CDW 之上（$T/t \sim 0.05$ 到 $0.1$ 之间），发现了一个广阔的 PDW 涨落区。

计算所得关键数据：

  1. 谱权重分布（Node-Antinode Dichotomy）：

    • 电子掺杂： 谱权重主要集中在反节点区域（Antinodal, $(0, \pi)$），这有利于形成传统的零动量 $d$ 波超导。
    • 空穴掺杂： 反节点区域的谱权重被严重压制，形成了“费米弧”（Fermi Arcs），谱权重集中在节点区域（Nodal, $(\pi/2, \pi/2)$）。

  2. 配对关联强度：

    • 在电子掺杂侧，$\mathbf{q}=0$ 的 $d$ 波配对关联随温度降低而发散，符合超导不稳定性。
    • 在空穴掺杂侧，$\mathbf{q}=0$ 处的关联微弱，但在 $\mathbf{Q}_{PDW} \approx (0, \pi)$ 处观察到了显著的配对峰。数据表明这种 PDW 关联是“跨弧配对”（inter-arc pairing）的结果，即两个费米弧上的节点电子相互配对，产生净动量。

  3. 实空间配对模式：

    • 电子掺杂侧展现出均匀的 $d$ 波配对，相邻键上的配对相位相反。
    • 空穴掺杂侧展示了 $y$ 方向上的 2 周期调制，即配对振幅在实空间交替变换正负号，这是 PDW 的典型特征。

性能数据：

  • 键维（$D$）： 最大键维达到 32,768。这是目前同类模型 MPO 模拟的顶峰。
  • 时空体积： 模拟覆盖了 $W \times L \times \beta = 8 \times 18 \times 32 \approx 4600$ 的有效体积，远超传统 DQMC 的能力（通常限制在 $8 \times 8 \times 5$）。
  • 收敛性： 研究通过外推 $1/D \to 0$ 验证了 PDW 峰的稳健性，证明有限键维效应并未定性改变结论。

代码实现架构

该研究主要基于 Julia 语言开发的高性能张量网络库。Julia 凭借其“接近 C 的速度，接近 Python 的易用性”，成为现代计算物理的首选。

核心软件包：FiniteMPS.jl

主要计算逻辑实现在 FiniteMPS.jl 库中。这是一个专门为一维和准二维量子多体系统设计的张量网络库，支持 MPS（矩阵乘积态）和 MPO 的各种操作。

  • 主要特性：
    • 实现了基于 TDVP 的 tanTRG 算法。
    • 完备的对称性支持：显式利用了 $U(1)_{charge} \times SU(2)_{spin}$ 对称性。这是处理 $D$ 超过 30000 的关键，对称性分区块显著减少了内存占用和收缩成本。
    • 变分算符展开技术，用于控制演化过程中的键维增长。

复现指南：

  1. 环境准备： 安装 Julia 1.10+。克隆项目仓库：git clone https://github.com/Qiaoyi-Li/FiniteMPS.jl。
  2. 配置模型： 在输入脚本中定义 Hubbard 模型的参数 $t=1.0, t'=-0.2, U=8.0$。设置圆柱宽度（如 $W=6$）和边界条件。
  3. 执行演化： 调用 tanTRG 函数。建议先从小键维（$D=1024$）开始预热，然后逐步增加 $D$。对于 $D=32768$ 的生产计算，需要至少 512GB 内存及多核并行环境。
  4. 数据后处理： 利用公式 (2) 和 (3) 计算 Green 函数和配对关联函数。通过 Fourier 变换获取动量空间谱图。

开源链接：

4. 关键引用文献，以及你对这项工作局限性的评论

关键引用文献：

  1. [17, 18]： 代表了之前的基态 DMRG 和 AFQMC 研究，本工作通过有限温度模拟回应了其中的配对强度之争。
  2. [56, 57]： 本研究方法的基础，即 tanTRG 和热态张量网络演化协议。
  3. [58, 60]： 关于 PDW 理论的经典综述，提供了将实验观察到的伪能隙与 PDW 涨落联系起来的理论框架。
  4. [61, 62]： 为模型选择了与实验一致的 $U$ 和 $t'$ 参数。

局限性评论：

  1. 单带模型的局限： 尽管单带 $t-t'$ Hubbard 模型捕捉了粒子-空穴不对称性，但真实的铜氧化物是多带体系（氧的 $p$ 轨道起重要作用）。特别是氧轨道上的配对和电荷转移能级可能对 PDW 的形成有进一步修正。
  2. 圆柱几何的离散动量： 圆柱宽度 $W$ 限制了 $k_y$ 动量的取值。虽然通过 PBC/APBC 平均有所改善，但在研究跨弧配对的精细动量结构时，仍然受限于动量点的稀疏性。
  3. 准二维 vs 真二维： 尽管数值上已经达到很大规模，但系统本质上仍是长条形。真实二维极限下的超导转变（BKT 转变）在目前框架下难以直接探测，只能通过关联函数的增长趋势进行推断。
  4. 配对场的引入： 为了计算静态响应，研究引入了局部配对场。虽然这有助于直接观察相位，但可能在某种程度上诱导了不稳定性，尽管作者通过改变场强验证了线性响应的可靠性。

5. 其他必要补充：物理机制探讨

为什么是 PDW 而不是 $d$ 轮超导？

这是本研究最深刻的物理启示。在空穴掺杂侧，由于范霍夫奇点（Van Hove Singularity）的位置以及 $t'$ 的影响，费米面在反节点区域发生了严重的拓扑形变或被伪能隙压制。传统的 $d$ 波超导需要反节点处的准粒子进行 $(k, -k)$ 配对，因为那里的态密度最高且配对势最强。

当反节点“失效”后，系统被迫在节点区域寻找配对机会。然而，节点区域的费米子具有非零的群速度和复杂的相位分布。本工作通过计算 anomalous Green’s function $F(k; Q)$ 揭示了，空穴掺杂侧的配对发生在两个费米弧（Fermi Arcs）之间。这种“跨弧配对”必然带有净动量，因为它无法在布里渊区中心找到完全对称的补偿。这直接导致了 PDW 态的出现。

与实验的关联

该工作预测的 PDW 涨落位于伪能隙区域的低能端。这与近年来的 STM 实验以及扫频 X 射线散射实验观察到的“超导涨落但在动量空间不均匀”的现象高度吻合。它暗示了 PDW 可能不是一个独立的竞争相，而是高温超导发生过程中的一个内禀阶段，是由费米面拓扑结构决定的必然结果。

总结与展望

这项工作展示了热态张量网络在处理强关联电子系统中的巨大潜力。随着计算能力的提升，未来我们可以期待更复杂的多带模型（如 Emery Model）在有限温度下的精准模拟，从而最终彻底揭开高温超导的神秘面纱。