Simulating fermionic fractional Chern insulators with infinite projected entangled-pair states

该研究通过变分优化具有$U(1)$对称性的费米子无限投影纠缠对态（iPEPS），成功模拟了分数陈绝缘体（FCI）相，揭示了表征该相所需的临界键维，并引入压缩方案有效计算了大单元胞的纠缠谱以验证其拓扑特性。

要点

想象一下，你正试图理解一群隐形粒子所进行的复杂舞蹈。这些被称为**费米子（fermions）的粒子有一个非常严格的规则：它们讨厌待在彼此相同的地点，并且其运动方式会创造出一种旋转的、手性的（chiral）模式。当这些粒子在网格上以特定方式排列时，它们会形成一种被称为分数量子陈绝缘体（Fractional Chern Insulator, FCI）**的物质态。这种状态非常奇特，因为它具有“分数”特性（例如携带电子的分数电荷），并且表现得像一种即使没有磁场也不会被卡住的流体。

问题在于，在计算机上模拟这场舞蹈极其困难。传统的方法就像是通过只看一个 3x3 的微小方块来观察这场舞蹈。你会错过大局，并被边缘效应所迷惑。

新方法：无限的地板
本文作者开发了一种新方法，使用一种名为 iPEPS（无限投影纠缠对态）的工具来模拟这场舞蹈。请将 iPEPS 想象成不是在看一个房间的快照，而是在绘制一张无限地板的蓝图。因为它专为无限网格而设计，所以它不会受到那些会让其他方法感到困惑的“边缘效应”的影响。它能让研究者看到粒子真实的、无限的舞蹈。

挑战：“不可行性”之墙
物理学中存在一个已知的规则（“不可行性定理”），即你无法使用简单的有限蓝图来完美描述这些旋转的、手性的舞蹈。这就像试图用直线画出一个完美的圆；你可以接近它，但总会留下一些锯齿状的边缘。

为了绕过这一点，团队使用了一个巧妙的技巧，即数学压缩。他们构建蓝图时使用了一个“键维数”（我们可以称之为细节水平或 D）：

• 低细节（D=4 到 6）： 蓝图太模糊了。粒子看起来像是形成了奇怪的、团块状的模式，这与真实的物理现象不符。
• 高细节（D=7 及以上）： 一旦他们将细节水平提高到 7，蓝图突然变得清晰起来。粒子的行为开始完全符合 FCI 应有的表现。作者发现，D=7 是“临界阈值”；低于此值，模拟是错误的；高于此值，模拟则能忠实于现实。

他们如何检验工作
为了确保他们的无限蓝图是正确的，他们观察了这场舞蹈的三个特定“特征”：

1. 格林函数（“回声”）： 他们检查了一个粒子的影响是如何扩散的。在健康的 FCI 中，这种影响应该迅速消散（就像安静房间里的喊叫声），但会留下一个微小的、模糊的“薄纱”般的尾迹。这个尾迹实际上是他们的模拟正在撞击“不可行性之墙”的迹象，但它非常微小，不足以破坏整体图像。
2. 配相关（“个人空间”）： 他们测量了两个粒子靠近彼此的可能性。在这种状态下，粒子会保持一定的距离（一个“相关空洞”），就像聚会上那些本能地避免站得太近的人一样。他们的模拟与一个著名的理论态——“劳格林态（Laughlin state）”的数学模型几乎完美契合。
3. 纠缠谱（“边缘指纹”）： 这是最重要的测试。尽管他们模拟的是无限地板，但他们可以通过数学手段进行“切割”来观察边缘。这些边缘的能量级就像是一个指纹。他们发现数据中出现了特定的数字序列 (1, 1, 2, 3, 5)。这个特定的序列是分数量子陈绝缘体的独特签名，证明他们成功模拟了这种奇异相。

秘诀：压缩
模拟这些无限地板需要巨大的计算能力。为了处理“边缘指纹”的数学运算，作者发明了一种压缩方案。想象一下尝试解决一个拥有百万块拼图碎片的拼图。他们并没有同时观察所有碎片，而是找到了一种方法，将碎片组合成更小、更易处理的块，同时又不丢失整体图像。这使得他们能够在足够宽的晶格“圆柱体”上运行模拟，从而观察到真实的模式，而这是以前的方法难以轻易做到的。

底线
这篇文章之所以是一项突破，是因为它首次成功使用了一种强大的新模拟工具（iPEPS）来建模一种非常困难的费米子分数量子陈绝缘体。他们证明了，只要给予模拟足够的“细节”（键维数至少为 7），它就能准确地重现这些粒子复杂的旋转行为，并与理论预测及其他数值方法相吻合。这为科学家们以更高的精度研究这些奇异材料打开了大门。

技术摘要

技术摘要：利用无限投影纠缠对态（iPEPS）模拟费米子分数量子陈绝缘体

问题陈述
分数量子陈绝缘体（FCIs）是分数填充窄陈能带中出现的分数量子霍尔（FQH）态的晶格模拟物。虽然近期二维莫尔材料中的实验实现重新激发了对这些相的研究兴趣，但模拟产生 FCIs 的微观哈密顿量仍然是一个重大的数值挑战。现有方法，如精确对角化（ED）和无限密度矩阵重整化群（iDMRG），受限于小簇或有限宽度的圆柱体，容易受到有限尺寸效应的影响。此外，尽管无限投影纠缠对态（iPEPS）已成功模拟了玻色子手征拓扑序（例如手征自旋液体），但此前尚未被应用于费米子拓扑序。一个具体的理论障碍是“禁止定理”（no-go theorem），该定理禁止存在具有精确有限相关长度的 iPEPS 形式的手征相，这往往导致标准算法的收敛缓慢。

方法论
作者通过集成多种算法进展，扩展了 iPEPS 框架以处理费米子系统，旨在优化针对蜂窝晶格上具有密度-密度相互作用的无自旋费米子 Haldane 模型的 $U(1)$ 对称费米子 iPEPS。

1. 模型与晶格构建： 研究重点关注一个包含最近邻（NN）、次近邻（2nd-NN）和第三近邻（3rd-NN）跳跃项的哈密顿量，其中包含诱导陈能带的复相位。蜂握晶格被粗粒化为一个布拉维晶格，其局部物理希尔伯特空间维度为 $d=4$（代表空态、A 位点占据、B 位点占据以及双占据态）。
2. 费米子编码与对称性： 为了高效处理费米子统计特性，作者利用了 $U(1)$ 对称张量。这种方法强制执行每个单元格内精确的粒子数守恒，避免了化学势调节。对于目标填充因子 $\nu=1/3$，构建了一个由三个不等价张量（$A_1, A_2, A_3$）组成的 $3 \times 3$ 单元格，并通过在 $A_1$ 上添加一个辅助腿来注入一个费米子，从而确保每个单元格恰好有三个费米子。
3. 通过固定点自动微分（AD）进行优化： 标准的虚时演化被基于自动微分（AD）的梯度优化所取代。一项关键创新是使用了固定点 AD 方法。由于 CTMRG 在处理手征态时需要大量的迭代次数（$O(10^3)$）才能收敛（因为长相关长度），标准的通过角转移矩阵重整化群（CTMRG）算法的反向传播在内存上是难以承受的。固定点方法直接对收敛后的固定点进行微分，消除了存储中间 CTMRG 步骤的需求，并绕过了内存瓶颈。
4. 纠缠谱（ES）压缩： 为了计算大单元格的动量分辨边缘纠缠谱，作者引入了一种压缩方案。矩阵乘积算符（MPO）表示的约化密度矩阵从大单元格（维度为 $D^M$）粗粒化为单个位点（截断维度为 $D_c$），使用的是基于高阶奇异值分解（HOSVD）的等距映射。这使得计算在更宽圆柱体（$W=5$）上的 ES 成为可能，而这在之前是计算上不可行的。

核心贡献与结果
作者成功优化了键维数高达 $D=9$ 的 iPEPS 态，并确定了能够忠实表征 FCI 相所需的临界键维数。

• 临界键维数 ($D_{min} = 7$)： 通过分析基态能量密度和电荷不均匀性，作者发现当键维数 $D < 7$ 时，无法捕捉到 FCI 相，而是倾向于具有电荷调制（charge modulations）的态。只有当 $D \ge 7$ 时，能量才会降至第一激发能隙之下（与 ED 结果相比），且电荷分布变得均匀。
• 体观测量：
  • 格林函数： 等时单粒子格林函数表现出快速的指数衰减（体相关长度 $\xi_{bulk}/a \approx 0.63$），随后是一个长程的“薄纱”（gossamer）尾部，这与手征态的禁止定理一致。
  • 对相关函数： 对相关函数 $g(x)$ 在短距离处显示出一个相关空洞，并迅速收敛至 1，这与 $\nu=1/3$ Laughlin 态的定性特征相匹配。