来源论文: https://arxiv.org/abs/2603.08650v2 生成时间: Mar 11, 2026 09:40

0. 执行摘要

三角晶格 S=1/2 反铁磁海森堡模型（TLHAF）是凝聚态物理中研究几何挫折与量子自旋液体（QSL）的范本体系。长期以来，在该模型中引入次近邻耦合 $J_2$ 被认为能在 $120^\circ$ 磁序与条纹序（Stripe order）之间诱导出非磁性的 QSL 相。然而，关于该 QSL 的本质——是带能隙的 $Z_2$ 自旋液体还是无能隙的 $U(1)$ 狄拉克自旋液体——学界一直存在激烈争论。此前大量的密度矩阵重整化群（DMRG）研究在有限宽度的柱状几何上发现两个能量极度接近的“拟简并”基态（通常称为偶数和奇数扇区），并多将其归因于 $Z_2$ 拓扑序的拓扑扇区。

本研究由慕尼黑大学（LMU）的 Oleksandra Kovalska 等人完成，发表于 2026 年（预印本 arXiv:2603.08650v2）。作者采用高精度的矩阵乘积态（MPS）模拟，配合 SU(2) 对称性保留及切空间 Krylov（TaSK）方法，对 $J_2/J_1 = 0.125$ 处的 YC6 柱状体系进行了系统对比。研究发现，这两个扇区的基态不仅在静态等时结构因子上表现出显著强度差异，其低能动力学激发谱更是呈现出定性的不同。这些证据表明，这两个基态不能简单地理解为同一个带隙 $Z_2$ 自旋液体的拓扑扇区。特别是奇数扇区展现出与 $120^\circ$ 磁序高度相似的特征，而偶数扇区则更符合 $U(1)$ 狄拉克自旋液体的预言。这一发现要求我们重新评估强挫折体系中“拟简并”基态的起源，并为实验探测 QSL 提供了新的理论判据。

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

1.1 核心科学问题

本项研究的核心在于：在强挫折三角晶格 J1-J2 模型中观察到的两个拟简并基态，究竟是拓扑序的体现，还是不同相竞争的结果？

在拓扑量子物态理论中，如果系统处于带能隙的 $Z_2$ 自旋液体相，那么在具有非平凡同伦群的流形（如圆柱面）上，应当存在多个简并的拓扑扇区。由于有限尺寸效应，这些扇区在数值模拟中表现为极其微小的能量差。然而，如果这两个态在局部可观测物理量（如键关联、动力学谱）上存在显著差异，那么“拓扑扇区”的假设将面临挑战。本工作正是要通过对比这两个扇区的静态和动态响应来验证这一假设。

1.2 理论基础：J1-J2 海森堡模型

研究对象为如下定义的 $J_1$-$J_2$ 海森堡哈密顿量：

$$H = J_1 \sum_{\langle i,j \rangle} \mathbf{S}_i \cdot \mathbf{S}_j + J_2 \sum_{\langle \langle i,j \rangle \rangle} \mathbf{S}_i \cdot \mathbf{S}_j$$

其中 $J_1$ 为最近邻（NN）交换相互作用，$J_2$ 为次近邻（NNN）交换相互作用。

  • 当 $J_2=0$ 时，基态表现为 coplanar $120^\circ$ 磁序。
  • 当 $J_2/J_1$ 足够大时，系统表现为 collinear stripe 磁序。
  • 在中间区域 $0.07 \lesssim J_2/J_1 \lesssim 0.15$，系统进入 QSL 候选区。

1.3 技术难点

  1. 能隙极小与高度简并：QSL 区域的能隙极小，传统的 DMRG 极易陷入局部极小值，或者在不同扇区之间跳变。
  2. SU(2) 对称性的保持：为了精确区分物理态并减少希尔伯特空间的维度，必须严格保持自旋旋转对称性。
  3. 动力学模拟的高分辨需求：传统的时间演化方法（如 tDMRG/TEBD）在处理长时关联和低能激发分辨力上存在局限，且计算量巨大。
  4. 扇区选择：在 SU(2) 保留的计算中，如何系统地“诱导”系统进入特定的扇区（尤其是能量略高的扇区）是一个技巧性难题。

1.4 方法细节：MPS 与 TaSK

作者采用了以下先进方案：

  • DMRG 与 SU(2) 对称性：使用 QSpace 张量库，利用非阿贝尔 SU(2) 对称性，这比单纯使用 $U(1)$ 对称性能更有效地压缩状态并确保物理态的纯正性。在扫频过程中引入了受控键扩展（controlled bond expansion）技术，防止状态过早收敛到错误的局域点。
  • 切空间 Krylov (TaSK) 方法：这是计算动力学结构因子（DSF）的核心技术。其基本思想是将哈密顿量投影到基态 MPS 的正交切空间 $\mathbb{V}^{\perp}$ 中。在切空间内构造 Krylov 基底： $$\mathcal{K}(\Psi_0) = \text{span} \{ |\Psi_0\rangle, H^{\perp}|\Psi_0\rangle, \dots, (H^{\perp})^{N_{kr}}|\Psi_0\rangle \}$$ 其中 $|\Psi_0\rangle = \mathbf{S}_k |\Psi_g\rangle$ 是初始激发态。通过 Lanczos 算法在此空间内对哈密顿量矩阵化，可以极高效率地获得低能频率响应。这种方法比基于时间演化的方法在频率分辨率上具有显著优势，尤其适合捕捉自旋液体中的连续谱特征。

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

2.1 Benchmark 体系配置

研究聚焦于 YC6 圆柱几何（宽度 $L_y=6$，长度 $L_x=36$）。圆柱绕 $a_2$ 方向闭合。参数选取为 $J_2/J_1 = 0.125$，这是经典模型的相变点，也是量子模型中 QSL 候选区的代表点。

2.2 能量与稳定性数据

作者在 DMRG 扫频过程中观察到了一个关键现象：

  • 当键维度 $D^* = 2048$ 时，能量密度出现突然下降，系统从 $|\psi_{\text{even}}\rangle$ 态跃迁到了 $|\psi_{\text{odd}}\rangle$ 态。
  • 能量差异：奇数扇区基态比偶数扇区低约 1.1%。
  • 正交性：两个态的保真度 $\mathcal{F} = |\langle \psi_{\text{even}} | \psi_{\text{odd}} \rangle| \sim 10^{-4}$，证明它们在热力学极限下是完全不同的物理态。

2.3 静态关联数据 (ETSF)

等时结构因子 $\chi(\mathbf{k})$ 的计算结果显示：

  • $120^\circ$ 序 ($J_2=0$)：在 Brillouin 区角落（K 和 K’ 点）有极尖锐的 Bragg 峰。
  • Stripe 序 ($J_2=0.2$)：在 M 和 M’ 点出现尖锐特征。
  • 偶数扇区 ($J_2=0.125$)：K 点和 M 点的强度分布相对均衡，表现出显著的软化，符合 $U(1)$ 狄拉克自旋液体的特征（具有宽广的关联分布）。
  • 奇数扇区 ($J_2=0.125$)：大部分光谱权重集中在 K 点，这与临近的 $120^\circ$ 磁序特征惊人地相似。最近邻键关联（Table I）显示，奇数扇区的键强更加各向同性（NN1, NN2, NN3 的差异极小），而偶数扇区则表现出明显的交替模式（Staggering pattern）。

2.4 动力学响应数据 (DSF)

动力学结构因子 $S(\mathbf{k}, \omega)$ 揭示了更深层次的差异：

  • 偶数扇区：在 K 点表现出最强的激发强度，但 M 点亦有显著贡献。这种分布与之前的研究一致，支持其作为 QSL 的身份。
  • 奇数扇区：展现出更类似于 $120^\circ$ 磁序的准粒子图像。虽然存在有限尺寸效应导致的能隙，但其色散曲线和权重分布高度暗示了其与磁序态的亲缘关系。

2.5 性能指标

  • 键维度 (Bond Dimension)：SU(2) 多重态数目达到 $D^* = 3000$，在动力学计算前压缩至 $D^* = 1000$。
  • Krylov 步数：$N_{kr} > 50$，使用连分数展开（CFE）进行谱平滑。

3.1 核心软件包：QSpace

本研究所依赖的核心张量库是 QSpace，由 Andreas Weichselbaum 开发。

  • 特点：支持任意非阿贝尔对称性（如 SU(N), Sp(N)），能够自动处理 Clebsch-Gordan 系数，极大地减少了强关联模拟中的内存开销和计算时间。
  • 获取途径：QSpace 部分算法集成在一些学术项目中，虽然其核心库可能不是完全开源的通用软件，但其算法逻辑在作者的文献 [42-44] 中有详尽描述。

3.2 动力学方法实现：TaSK

TaSK 方法是近年来 MPS 领域的重要进展，其复现关键在于切空间投影的实现。

  • 参考实现：可以参考 ITensor 社区中关于切空间算子（Tangent Space Operators）的实现思路，或者基于 Julia 语言的 MPSKit.jl 库。MPSKit.jl 提供了成熟的切空间和有限/无限 MPS 处理函数，非常适合进行类似的动力学模拟。
  • 开源 Repo 相关项：

3.3 复现指南（关键参数建议）

  1. 初始化：在 YC6 几何上，初态需要通过随机张量初始化，并确保在指定的对称性扇区（Even/Odd）内进行。
  2. 扫频策略：
    • 前几个 half-sweeps 使用较大的噪声项或混合参数（如 DMRG3S 逻辑 [58, 59]）来跳出局部极小值。
    • 键维度逐步增加：$512 \to 1024 \to 2048 \to 3072$。
  3. 动力学拓扑投影：执行 TaSK 时，必须先确保基态 MPS 已经完全收敛。在切空间内进行 Lanczos 迭代时，需保留至少 50-100 个 Krylov 向量以保证低能分辨。
  4. 后处理：使用连分数展开（Continued Fraction Expansion）处理 Lanczos 产生的 $a_n, b_n$ 系数，平滑参数建议设置为 $\sigma = 0.4$。

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

4.1 关键引用文献

  1. [1] P. Anderson (1973): RVB 理论的起点，提出 TLHAF 可能存在 QSL。
  2. [14] Z. Zhu and S. R. White (2015): 首次通过 DMRG 深入探讨 J1-J2 模型的 QSL 性质，并提到扇区问题。
  3. [33] M. Drescher et al. (2023): 讨论了对称性破缺与液体相的动力学特征，是本工作的直接理论背景。
  4. [47] O. Kovalska et al. (2025, arXiv): TaSK 方法的详细技术文档。
  5. [18] S. Hu et al. (2019): 关于狄拉克自旋液体的关键 DMRG 证据。

4.2 局限性评论

尽管本工作在技术上非常扎实，但仍存在以下局限：

  1. 有限宽度限制：所有的计算均在 YC6 圆柱（$L_y=6$）上完成。在三角晶格中，$L_y=6$ 仍然是一个相对窄的几何，有限尺寸诱导的各向异性可能极大地改变了能隙的大小和性质。正如作者在结论中提到的，当宽度增加到 YC8 或 YC10 时，偶数和奇数扇区的能量顺序是否会反转仍未可知。
  2. 真实基态的定义：奇数扇区能量更低，但这是否意味着它是真正的 2D 极限下的基态？或者是圆柱几何对 $120^\circ$ 序某种“残余”的偏好？目前的计算尚不能给出定论。
  3. U(1) vs. Z2 的最终判据：虽然动力学谱更支持偶数扇区作为狄拉克液体，但对于奇数扇区的“非拓扑性”解释更多是基于现象学的对比。目前仍然缺乏一个严格的理论框架来定量解释为什么奇数扇区会表现出如此强的准粒子特征，而它在物理上又为何能与偶数扇区如此接近。
  4. SU(2) 的限制：虽然保留 SU(2) 对称性非常优美，但在实验体系中，微小的各向异性（如 DM 相互作用）会破缺这一对称性。因此，该理论预言在实际材料（如 $NaYbSe_2$）中的适用性还需考虑对称性破缺的摄动。

5. 其他必要的补充

5.1 物理直觉：为什么会出现这种扇区差异？

在圆柱几何中，绕圆柱一周的路径是非平凡的。偶数和奇数扇区的区别本质上在于通过圆柱中心的磁通量（Flux）。在 $Z_2$ 自旋液体中，这对应于粘合 vison（$Z_2$ 旋涡）。然而，本工作揭示了一个更深刻的可能性：由于三角晶格的几何挫折，圆柱的周期性边界条件实际上对不同的磁序构型产生了不同的“张力”。奇数扇区可能更好地调和了这种张力，从而表现出类似磁序的特征；而偶数扇区则保留了更多的液体本性。这暗示了在数值模拟中，所谓的“拓扑扇区”有时可能是不同物理相在有限尺寸下的某种“投影”。

5.2 关于 TaSK 方法的先进性补充

TaSK 方法之所以优于传统的频率域方法，是因为它利用了 MPS 的流形几何性质。在 MPS 切空间内，自旋算子的作用被线性化了。这意味着我们可以直接处理哈密顿量的本征激发，而不需要像 FFT 那样依赖长时间的演化。对于研究自旋液体这种具有宽能谱连续区（Continuum）的体系，TaSK 能够清晰地辨别出光谱权重是集中在某些准粒子峰上，还是均匀散布在 Brillouin 区，这对于区分 $U(1)$ 液体至关重要。

5.3 对后续研究的启示

本工作对量子自旋液体搜索具有重要的指导意义：

  • 实验家：在分析中子散射数据时，不应只关注峰的位置，更应关注强度的各向同性程度以及色散的宽窄。如果观察到类似奇数扇区的特征，可能意味着系统更偏向于受挫磁序而非真正的 QSL。
  • 理论家：在进行 DMRG 计算时，必须极度小心能量极度接近的态。简单的能量判据不足以认定物相，必须考察动力学响应。同时，本工作展示了 SU(2) 对称性在处理这类敏感问题时的必要性，这应当成为未来此类研究的标准配置。

5.4 结论总结

Kovalska 等人的工作通过精密的数值实验，打破了 J1-J2 模型中拟简并态的“拓扑神话”。他们证明了在特定几何下，系统可以同时存在表现得像自旋液体的态和表现得像受挫磁序的态，且后者在能量上略胜一筹。这一发现不仅深化了我们对三角晶格强挫折物理的理解，也为未来在更大尺寸、更复杂模型下的量子模拟设立了新的技术标杆。