来源论文: https://arxiv.org/abs/2603.12309v1 生成时间: Mar 22, 2026 09:16

相互作用驱动的扩展哈伯德模型中的亚铁磁条纹态：深度解析

0. 执行摘要

在强关联电子体系的研究中，二维哈伯德模型（Hubbard Model）一直是理解铜氧化物高温超导体等材料中涌现现象的核心框架。传统研究多聚焦于现场库仑排斥力 $U$ 诱导的反铁磁（AFM）条纹序。然而，真实材料中的长程库仑相互作用——尤其是最近邻排斥力 $V$——对基态特性的定性影响长期以来由于计算复杂性而悬而未决。

近期由 Chunhan Feng、Miguel A. Morales 和 Shiwei Zhang 发表的研究成果表明，在二维平方晶格的扩展哈伯德模型中，当最近邻排斥与现场排斥的比值 $V/U$ 超过临界点（约 0.25）时，系统会发生剧烈的基态重构。研究通过结合高精度的辅助场量子蒙特卡洛（AFQMC）和密度矩阵重整化群（DMRG）方法，发现了一种新型的调制亚铁磁条纹（Modulated Ferrimagnetic Stripes），该序与棋盘格电荷密度波（CDW）相互交织。这一发现不仅挑战了传统的条纹相认知，也为在冷原子模拟和莫尔超晶格材料中寻找新型磁性织构提供了理论指引。

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

1.1 核心科学问题

哈伯德模型在掺杂情况下展现出的磁、荷条纹序被认为是理解伪能隙（Pseudogap）和非常规超导的关键。传统的 $t-U$ 模型预测了波长为 $\lambda = 2/\delta$ 的反铁磁条纹。然而，最近邻排斥 $V$ 的引入会产生竞争性倾向：$U$ 倾向于 AFM 关联，而 $V$ 倾向于棋盘格 CDW。核心问题在于：当这两种倾向在掺杂系统中正面交锋时，原有的条纹序会塌陷、增强还是演化成全新的物态？

1.2 理论基础：扩展哈伯德模型

该研究探讨的哈密顿量如下：

$$ H = -t \sum_{\langle i,j \rangle, \sigma} (\hat{c}_{i\sigma}^\dagger \hat{c}_{j\sigma} + \text{h.c.}) + U \sum_i \hat{n}_{i\uparrow} \hat{n}_{i\downarrow} + V \sum_{\langle i,j \rangle} \hat{n}_i \hat{n}_j $$

其中 $V$ 项惩罚了占据相邻格点的电子对。在半满（Half-filling）情况下，$V \gtrsim U/4$ 会稳定 CDW。但在掺杂（Away from half-filling）情况下，系统的自由度显著增加，电荷流体、自旋激发与非局部相互作用的耦合变得极其复杂。

1.3 技术难点：量子多体计算的挑战

  1. 负符号问题（Sign Problem）：在引入 $V$ 项后，辅助场量子蒙特卡洛（AFQMC）的辅助场数量增加，导致费米子符号问题显著恶化。常规方法在低温和低掺杂区难以收敛。
  2. 相互作用范围与纠缠：DMRG 在处理一维或准一维（圆柱形几何）时非常强大，但随着相互作用范围的扩大（如 $V$ 项），算符的纠缠熵迅速增加，要求极大的键维（Bond Dimension）才能保证精度。
  3. 多相竞争：条纹序与超导、电荷波等相的能量差极小（通常在 $10^{-3}t$ 数量级），需要极高的数值分辨率来确定基态。

1.4 方法细节：CP-AFQMC 与 DMRG 的协同效应

研究者采用了一种“双管齐下”的策略：

  • 约束路径 AFQMC (CP-AFQMC)：通过引入试验波函数 $\Psi_T$ 来消除符号问题。为了提高精度，研究采用了自洽 CP-AFQMC 方案。其核心在于利用 AFQMC 计算得到的密度矩阵去更新 Hartree-Fock (UHF/GHF) 计算中的有效参数 $U_{\text{eff}}$ 和 $V_{\text{eff}}$，从而迭代优化试验波函数，减小约束路径带来的系统误差。
  • DMRG 辅助验证：在较小的圆柱形格点（如 $16 \times 4$）上，DMRG 提供了精确的参考基准，用于验证 AFQMC 的自洽解是否陷入局部极小值。一旦在小尺寸上取得一致，AFQMC 即可扩展到更大的二维格点（如 $24 \times 8$ 或 $36 \times 12$）以排除尺寸效应。

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

2.1 Benchmark 体系设置

  • 晶格几何：主要采用方形晶格的圆柱形几何（Lx 沿 x 轴，Ly 沿 y 轴）。典型尺寸包括 $16 \times 4$、$24 \times 8$ 以及 $36 \times 12$。
  • 参数空间：固定 $U/t = 8$（强关联区），扫描 $V/U$ 从 0 到 0.5，掺杂率 $\delta$ 涵盖 $1/12$ 到 $1/4$。
  • 边界条件：x 方向为开边界（OBC），y 方向为周期性边界（PBC）。同时施加微弱的边界钉扎场（$v_p = 0.1$）以辅助可视化条纹图样，但不影响体相性质。

2.2 核心计算数据

  1. 相图重构：当 $V/U < 0.25$ 时，系统维持传统的 AFM 条纹。当 $V/U \gtrsim 0.25$ 时，相图进入“亚铁磁条纹区”。此时，电荷密度不再是简单的正弦调制，而是展现出强烈的棋盘格 CDW 构型。
  2. 自旋密度特征：在自旋平衡（$N_\uparrow = N_\downarrow$）的扇区中，自旋密度在“正-近零”和“负-近零”畴之间交替。在单个畴内，自旋排列呈现亚铁磁特征而非纯粹的反铁磁性。当自旋约束解除时，系统选择具有净磁化强度的均匀亚铁磁态。
  3. 结构因子 $S(k)$ 与 $C(k)$：
    • 在 $V=0$ 时，$S(k)$ 在 $((1-\delta)\pi, \pi)$ 处有峰值。
    • 在 $V/U = 0.4$ 时，出现新峰值 $(\delta\pi, 0)$，标志着畴间的周期性交替，同时 $(\pi, \pi)$ 处的强度显著增强。

2.3 计算性能与收敛性

  • 自洽收敛：图 2 展示了 AFQMC 经过 3-4 次自洽迭代后即可与 DMRG 结果完美重合。这种稳定性证明了该物态不是试验波函数的赝像，而是关联驱动的稳态。
  • 尺寸稳定性：在 $Lx = 36$ 的超长圆柱上，条纹周期保持稳定，证明了结论在热力学极限下的可靠性。

3.1 核心算法实现路径

复现该研究需要两套核心工具：一套支持扩展哈伯德相互作用的 DMRG 代码，以及一套高性能的 CP-AFQMC 代码。

3.1.1 DMRG 部分

推荐使用 ITensor 库（基于 C++ 或 Julia）。

  • 实现细节：
    • 需要手动构建包含 $V \sum \hat{n}_i \hat{n}_j$ 项的 MPO（Matrix Product Operator）。对于二维圆柱，跳跃项和 $V$ 项会涉及跨行的长程相互作用，需优化 MPO 的排布。
    • 键维设置：对于 $Ly=4$ 的系统，$M \approx 2000$ 即可收敛；对于 $Ly=6$ 或更高，需要 $M > 8000$ 并结合并行化。

3.1.2 AFQMC 部分

该研究使用的是 Flatiron Institute 内部开发的 AFQMC 软件包。对于开源界，最接近的替代品是 QMCPACK 或 pyscf-qmc。

  • 复现关键点：
    • Trial Wavefunction：必须实现广义 Hartree-Fock (GHF) 以允许非共线的磁性配置。
    • Local Energy Evaluation：由于引入了非局部相互作用 $V$，局部能量的计算复杂度从 $O(N)$ 提升至 $O(N^2)$，需要针对四费米子项进行张量分解优化。

3.2 复现指南步骤

  1. 初始化：在给定的 $\delta, U, V$ 下运行 UHF 程序，施加钉扎场获得初始试验波函数。
  2. 第一轮 AFQMC：使用该 UHF 波函数作为约束路径。记录输出的 $\langle n_{i\sigma} \rangle$。
  3. 自洽更新：调整 HF 模型中的 $U_{\text{eff}}$ 和 $V_{\text{eff}}$，使得 HF 密度分布与上一步 AFQMC 输出最接近。
  4. 循环：重复上述过程直至能量误差低于 $10^{-4}t$。

3.3 相关资源链接

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

4.1 关键引用文献

  1. Hubbard (1963) [Ref 1]: 定义了基础模型。
  2. Kivelson et al. (2003) [Ref 2]: 确立了条纹序在高温超导研究中的地位。
  3. White & Scalapino (1992, 1993) [Ref 27, 28]: 开创了 DMRG 在哈伯德模型中的应用。
  4. Zhang et al. (1995, 1997) [Ref 25, 26]: 约束路径 AFQMC 的奠基性工作。
  5. Qin et al. (2020) [Ref 5]: 确定了单纯 $t-U$ 模型中不存在超导的长程关联，为引入 $V$ 项的研究奠定了背景。

4.2 工作局限性评论

尽管该研究在数值上具有极高说服力，但仍存在以下局限性：

  1. $V/U$ 的物理范围：研究中定义的临界点 $V/U \sim 0.25$ 在真实铜氧化物中是否可达？真实材料的有效排斥力通常受到复杂屏蔽效应影响，$V$ 的取值仍有争议。
  2. 超导关联的缺失：虽然揭示了新型磁性条纹，但文章未详细讨论该亚铁磁序对超导配对（Pairing）的影响。通常 CDW 会抑制超导，而亚铁磁性的引入可能进一步破坏库珀对的单态配对。
  3. 自洽路径的局限：自洽 CP-AFQMC 虽然减小了误差，但其基准仍依赖于 HF 型试验波函数。对于极其强关联的系统，是否漏掉了某些非 HF 能描述的非平凡拓扑序或液态特征？

5. 其他必要补充：物理直觉与未来展望

5.1 物理直觉：为什么会出现亚铁磁性？

我们可以通过电荷排布来直观理解：$V$ 项强烈排斥相邻电子，迫使电子形成棋盘格排列（即子格 A 占据多，子格 B 占据少）。在反铁磁背景下，这种棋盘格电荷调制打破了子格对称性。由于两个子格上的平均占据数不等，导致原本相互抵消的自旋分量出现了净余，从而演化成亚铁磁性。这是一种典型的“几何挫折”与“电荷关联”共同作用的结果。

5.2 实验可观测性

  1. 冷原子平台：利用里德堡原子（Rydberg Atoms）可以实现长程排斥力，该研究预测的 $V/U$ 范围在目前的冷原子实验（如 Harvard 或 Max Planck 组的实验装置）中是完全可调控的。
  2. 莫尔超晶格：在扭曲双层石墨烯或过渡金属硫族化合物（TMDs）中，通过栅极电压调节载流子浓度和屏蔽环境，可以直接观测到这种交织序。其特征磁化强度可通过自旋极化扫描隧道显微镜（SP-STM）探测。

5.3 结论与启示

这项工作标志着哈伯德模型研究从“单纯的现场关联”迈向“非局部相互作用主导”的新阶段。它告诉我们，电子系统中微小的参数扰动（如 $V/U$ 的微调）足以诱导出全新的对称性破缺态。对于致力于量子材料设计的人员而言，这意味着除了调节 $U$ 和 $t$，控制非局部库仑相互作用将是操纵复杂自旋织构的又一核心手段。

通过对 $t'$（次近邻跳跃）的额外分析，研究进一步证明了这种亚铁磁条纹在不同空穴/电子掺杂下的鲁棒性。这为未来构建更完整的强关联电子相图提供了至关重要的拼图。