空头 18 腿 t-J 圆柱体上的大孤立条纹

Large Isolated Stripes on Short 18-leg $t$-$J$ Cylinders

arXiv: 2512.15714v1

论文信息

标题: Large Isolated Stripes on Short 18-leg tt-JJ Cylinders

作者: Tizian Blatz, Sebastian Paeckel, Ulrich Schollwöck, et al.

发布日期: 2025-12-17

arXiv ID: 2512.15714v1

PDF 链接: 下载 PDF

3 分钟速览

  • 研究问题：这篇论文要解决铜氧化物高温超导体理论模型中「条纹相」的填充分数范围及其微观形成机制问题，尤其是如何克服传统数值研究中严重的有限尺寸效应。
  • 核心方法：作者采用密度矩阵重整化群算法，在一种非标准几何结构——18 腿圆柱体（宽度远大于长度）上计算基态，并通过定制的矩阵乘积态映射将计算瓶颈从周期方向长度转移到短边方向。
  • 关键结果：在空穴掺杂侧（t′<0t'<0），单个条纹可以在填充分数 ν≈1/3\nu\approx 1/3 到 ν=1\nu=1 的宽范围内稳定形成；而在电子掺杂侧（t′>0t'>0），条纹仅在 ν=1\nu=1 附近狭窄区间出现（图 1）。
  • 主要局限：系统的 xx 方向长度仅为 4 个格点，虽然 yy 方向达到 18 个格点，但极端宽高比的几何结构可能影响向真实二维热力学极限的外推。
  • 适合读者：对强关联电子体系、铜氧化物超导机理、张量网络数值方法感兴趣的理论凝聚态物理研究者和研究生。

论文背景和研究动机

铜氧化物高温超导体在低温下除超导序外，最显著的序参量之一是自旋-电荷条纹序。这种条纹相的特征是在反铁磁背景中，电荷密度峰将不同方向的反铁磁畴分隔开来。理解条纹相的形成机制及其与超导的共存关系，是揭示非常规超导机理的关键线索。

然而，数值上研究条纹相面临一个严重挑战：即使最先进的多体计算方法所能达到的系统尺寸也相当有限，而条纹序对有限尺寸效应极其敏感。在这些小系统中，条纹通常很短，难以将配对物理与条纹物理清晰分离。不同填充分数条纹之间的基态能量差可能极小，使得基态收敛困难，并强烈依赖于掺杂数目、平行于条纹波矢方向的系统尺寸等额外的有限尺寸参量。

面对这些困难，近年来的大量工作要么通过精细的有限尺寸外推和多方法对比来逼近热力学极限，要么利用变分方法解锁前所未有的系统尺寸。但现有的主流研究均采用低宽高比（Ly/Lx≲1L_y / L_x \lesssim 1）的圆柱几何，条纹长度受限，且可实现的填充分数被粗粒化。

本文作者提出一个互补视角：与其在有限掺杂下研究多条短条纹，不如研究一条长而孤立的条纹如何形成。通过分析单个条纹在反铁磁背景上的稳定填充区间，可以从更基本的层面理解条纹的微观起源，同时规避多条条纹耦合带来的复杂性。

核心方法和技术细节

模型哈密顿量

论文采用 tt-t′t'-JJ 模型作为描述铜氧化物物理的微观模型。哈密顿量写为：

H^t−t′−J=−t∑⟨i,j⟩,σ(c~i,σ†c~j,σ+H.c.)−t′∑⟨⟨i,j⟩⟩,σ(c~i,σ†c~j,σ+H.c.)+J∑⟨i,j⟩(S^i⋅S^j−n~in~j4)\hat{H}_{t-t'-J} = -t \sum_{\langle i,j\rangle,\sigma}(\tilde{c}^\dagger_{i,\sigma}\tilde{c}_{j,\sigma}+ \text{H.c.}) - t' \sum_{\langle\langle i,j\rangle\rangle,\sigma}(\tilde{c}^\dagger_{i,\sigma}\tilde{c}_{j,\sigma}+ \text{H.c.}) + J \sum_{\langle i,j\rangle}\left(\hat{\mathbf{S}}_i\cdot\hat{\mathbf{S}}_j - \frac{\tilde{n}_i\tilde{n}_j}{4}\right)

其中 c~i,σ†=c^i,σ†(1−n^i,−σ)\tilde{c}^\dagger_{i,\sigma} = \hat{c}^\dagger_{i,\sigma}(1-\hat{n}_{i,-\sigma}) 是禁止双占的投影费米子产生算符。参数固定为 t=1t=1（能量单位）、J=1/2J=1/2、t′∈[0;±0.2]t'\in[0; \pm0.2]。改变 t′t' 的符号等价于通过粒子-空穴变换交换电子和空穴掺杂，因此可以同时探索掺杂相图的两侧。

非标准几何与 MPS 映射

该项工作的关键技术突破在于采用非标准的圆柱几何：Lx×Ly=4×18L_x \times L_y = 4 \times 18，即宽高比 Ly/Lx=4.5L_y/L_x = 4.5，远大于传统 DMRG 研究中的 ≲1\lesssim 1。这种结构使得条纹在周期方向上可长达 18 个格点，明显长于文献中通常的条纹。

然而，宽圆柱意味着周期方向长，传统的一维矩阵乘积态路径映射会使最近邻相互作用的 MPS 范围达到 LyL_y 量级，导致计算代价呈指数增长。为解决这一瓶颈，作者设计了一种定制的 MPS 映射方案，其核心思想是将 yy 方向（周期方向）按类似一维周期系统的方式折叠。如此一来，yy 方向最近邻键的 MPS 范围变为 2Lx2L_x，而 xx 方向的最近邻键范围仅为 1。由于 Lx=4L_x=4，计算复杂度由指数依赖于 LyL_y 转变为指数依赖于 2Lx=82L_x=8。这使得研究 Ly=18L_y=18 的圆柱在计算代价上等同于传统映射下 Ly=8L_y=8 的圆柱。

条纹的识别

作者通过交错跨晶格自旋关联函数来识别条纹的形成：

CS^z(Lx)=1Ly(−1)Lx−1∑y=1Ly⟨S^(1,y)zS^(Lx,y)z⟩C_{\hat{S}^z}^{(L_x)} = \frac{1}{L_y}(-1)^{L_x-1} \sum_{y=1}^{L_y} \langle \hat{S}^z_{(1,y)} \hat{S}^z_{(L_x,y)} \rangle

正值对应未被扰动的反铁磁序，负值标志着反铁磁自旋畴壁的出现，即条纹的形成。这一简单而有效的判据构成了后续全部相图分析的基础。

创新点和贡献

1. 提供互补视角并验证文献结果

本文从单条纹形成这一完全不同视角逼近条纹相问题，所得的填充分数范围与大量已有文献高度一致（图 1）：在空穴掺杂侧可稳定 1/3≲ν≲11/3 \lesssim \nu \lesssim 1，电子掺杂侧仅在 ν=1\nu=1 附近。这一非平凡的一致性表明，条纹相的整体物理由单个条纹的稳定性主导，条纹间耦合起次要作用。这为理解早期提出的邻近效应机制和近期报道的条纹周期高压缩性提供了微观基础。

2. 揭示两种不同的条纹形成机制

通过分析自旋关联随填充分数的变化，并结合高概率直积态的快照，论文明确区分出两个条纹形成区间：

  • 高填充区间（ν≈1\nu\approx 1）：可由简单的「挤压空间」图像理解——在反铁磁背景中插入一行空穴，自旋关联被整体位移一个格点，自然产生畴壁。该机制几乎不依赖 t′t'，体现为不同 t′t' 下自旋关联的曲线高度重合。
  • 低填充区间（0<ν<10<\nu<1）：条纹形成源自配对物理。在 t′≤0t'\le 0 时，扩展的距离-3 对可促成条纹形成，出现的阈值可低至 ν=1/3\nu=1/3；在 t′>0t'>0 时，掺杂空穴形成紧束缚对且不破坏反铁磁序，CS^z(Lx)C_{\hat{S}^z}^{(L_x)} 保持正值，条纹不能形成。这种明显的 t′t' 依赖性直接对应到单对掺杂体在不同 t′t' 下的实空间结构差异（图 2）。

3. 建立从微观配对到宏观条纹的桥梁

论文的一个重要成就是将条纹相的复杂多体行为追溯至其基本构成单元——单对掺杂体。前期研究已发现，电子掺杂侧倾向于形成紧束缚对，空穴掺杂侧则对应扩展对并伴随自旋畴壁的形成。本文在更宽的系统上证实，多个这样的对要么形成集体条纹不稳定性（t′≤0t'\le0），要么和平共存于未被扰动的反铁磁背景中（t′>0t'>0）。由此，非整数对条纹（nIPS）可被理解为松散束缚对的集体不稳定性表现。

实验结果分析

条纹形成相图

图 1(b)展示了系统的中心结果。每个点表示一次独立的基态计算，颜色代表 CS^z(Lx)C_{\hat{S}^z}^{(L_x)} 的值。冷色调（负值）对应条纹态，暖色调（正值）对应无条纹态。在 t′=−0.2t'=-0.2（空穴掺杂侧），条纹填充分数从约 ν=1/3\nu=1/3 一直延伸到 ν=1\nu=1，涵盖文献中报道的几乎所有填充分数。在 t′=0t'=0，条纹仅出现在 ν≈1/2\nu\approx 1/2 至 11 区间。在 t′=0.2t'=0.2（电子掺杂侧），仅 ν=1\nu=1 附近出现负关联。

微观结构的高概率直积态分析

作者采用高概率直积态方法探测基态的电荷与自旋构型。图 2 展示了固定 ν=1/2\nu=1/2 时三种 t′t' 下的典型构型：

  • t′=−0.2t'=-0.2：掺杂空穴基本分离，形成清晰的自旋畴壁
  • t′=0.0t'=0.0：掺杂空穴形成紧束缚的次近邻对，但仍能产生畴壁
  • t′=0.2t'=0.2：掺杂空穴结合成紧束缚对，但反铁磁背景保持完整，无畴壁

该分析直接揭示了从单对到多对的过渡如何产生或抑制条纹序。

挤压空间模型的定量验证

图 3 对比了真实的 CS^z(Lx)C_{\hat{S}^z}^{(L_x)} 与仅基于电荷结构和未掺杂反铁磁关联的挤压空间预测。预测模型捕捉到了 ν=1\nu=1 附近的畴壁形成及其对 t′t' 的不敏感性（不同 t′t' 曲线在此重合），但在低填充区间完全失效——真实的关联呈现强 t′t' 依赖性。这种预测与实际的偏差精确地标记出配对物理起作用的分界线。

局限与待解决问题

本文的系统几何 Lx=4L_x=4，虽然 LyL_y 达到 18，可视为宽度方向上的强约束。向真正的二维热力学极限外推时，LxL_x 的增加是否会改变条纹填充区间的定量边界，目前尚无答案。论文通过对比更小系统的计算结果说明 4×18 结果已具有较好的定性鲁棒性，但仍无法完全排除 LxL_x 增大时新物理出现的可能性。

挤压空间模型成功描述 ν=1\nu=1 附近的物理，但在低填充区间失效，说明需要一个能同时包含配对关联与自旋关联的统一理论框架。目前论文仅通过高概率直积态快照给出物理解释，尚未提出可定量预测低填充条纹形成阈值的有效模型。

从量子模拟实验角度看，论文指出低填充条纹的形成依赖于配对能尺度（约 0.1t−0.2t0.1t-0.2t），当前冷原子实验的温标刚进入此区间。如何在实际实验中区分挤压空间机制主导的填充条纹和配对机制主导的部分填充条纹，以及它们在输运和谱学响应上的可观测差异，是有待进一步探讨的问题。

本文的工作仅针对铜氧化物超导体正方的 tt-t′t'-JJ 模型。对于更现实的材料描述，如多层结构、额外跃迁项、电子-声子耦合等的影响，均未被考虑。此外，条纹与超导之间的动态共存、条纹序的熔化温度等有限温度性质也超出本工作的基态研究范畴。