Efficient Simulation of Sparse, Non-Local Fermion Models

本文提出了一种利用辅助费米子消除稀疏非局域费米子模型中约旦-维格纳弦的编码方案，从而将长时 Trotter 化时间演化的电路深度开销从乘性 $O(\log N)$ 因子降低为加性项，并在量子比特硬件上实现了渐近最优性能。

要点

想象一下，你正试图用一台标准的量子计算机来模拟一种称为费米子（物质的基本构成单元，如电子）的复杂粒子舞蹈。这些计算机使用的语言与费米子不同；它们使用“量子比特”（可以是 0、1 或两者同时存在的比特）。

为了让计算机理解费米子，科学家必须将费米子的规则翻译成量子比特的规则。问题在于，费米子有一条非常具体且棘手的规则：如果你交换其中两个，整个系统的符号就会翻转。在标准的翻译方法（称为Jordan-Wigner 变换）中，这条规则迫使计算机检查两个粒子之间的每一个量子比特，以确保符号正确。

问题：“长链”

把这想象成在一个巨大的体育场里玩“传话”游戏。如果玩家 A（在一端）想与玩家 B（在另一端）交谈，他们必须让消息通过站在他们之间的每一个人低声传递。在量子术语中，这是一串“长链”操作。

如果粒子相距甚远，这条“链”就会变得极其漫长。在量子计算机上，长链意味着模拟需要很长时间并消耗大量资源。这对于稀疏模型尤其糟糕，在这些模型中，粒子可能只与少数特定的邻居相互作用，但这些邻居可能位于系统的任何位置。

解决方案：添加“助手”

本文的作者 Reinis Irmejs 和 J. Ignacio Cirac 想出了一个巧妙的技巧来缩短这些长链。

1. 设置：添加“辅助”邻居
想象你系统中的每个粒子都有一支由辅助粒子（称为辅助费米子）组成的小团队，它们就住在粒子旁边。这些助手不会改变系统的物理性质；它们只是在那里帮助翻译。

2. 魔法技巧：稳定子
作者创建了一组特殊的规则，称为稳定子。可以将这些视为助手之间的“握手”协议。

• 在模拟开始之前，他们将所有助手准备在一个非常具体且同步的状态中，使它们都同意握手规则。
• 一旦设定好这种状态，助手们就充当桥梁。它们允许远处的粒子通过其本地助手直接通信，从而无需通过整个体育场低声传递。

3. 结果：切断链条
由于这种设置，“长链”操作消失了。计算机不再需要检查两个粒子之间的每一个量子比特，而只需检查固定数量的量子比特（本地粒子及其直接助手）。

代价：一次性费用

这里有一个陷阱，但这是一笔公平的交换。

• 设置成本：准备那些同步的助手需要一些时间和努力，而且是在最开始的时候。这就像在戏剧开始前搭建一个复杂的舞台。随着系统变大，这个初始设置需要的时间会稍长一些（具体来说，它与系统大小的对数成正比，即 $O(\log N)$）。
• 回报：一旦舞台搭建完毕，助手们将永远保持在那个完美的状态中。它们不需要在模拟的每一步都重置或重新准备。

为什么这很重要

过去，在量子比特计算机上模拟这些稀疏系统，比在理论上的“理想”费米子计算机上模拟要慢，其速度差距随着系统大小的增长而扩大（存在一个乘性的 $O(\log N)$ 惩罚）。

使用这种新方法：

• 初始设置是唯一具有该惩罚的部分。
• 对于长时间模拟（长时间运行舞蹈），每一步的成本变为常数。
• 现在，在标准量子比特计算机上运行模拟的总时间，除了一个小常数因子外，已与理想费米子计算机的性能相匹配。

结论

这篇论文证明，你不需要一台专门的“仅费米子”计算机来获得最佳结果。通过添加少量辅助粒子并进行一次性设置，你可以使标准量子比特计算机几乎像理论上的理想硬件一样高效地模拟稀疏费米子系统。它将一个“缓慢且随规模增长”的问题，转变为针对长时间模拟的“快速且恒定”的问题。

技术摘要

以下是 Reinis Irmejs 和 J. Ignacio Cirac 所著论文《稀疏非局域费米子模型的高效模拟》的详细技术总结。

1. 问题陈述

模拟相互作用的费米子系统是近期量子计算机的主要应用场景之一。然而，将费米子算符映射到量子比特硬件存在一个主要瓶颈。

• Jordan-Wigner (JW) 开销： 标准编码（如 JW）利用非局域的“弦”算符（Z 门的乘积）将费米子的反对易关系映射到量子比特算符。对于站点 $i$ 和 $j$ 之间的非局域相互作用，泡利权重按 $O(|i-j|)$ 缩放。
• 稀疏非局域模型： 许多物理相关的模型（例如 SYK 模型的稀疏实例、随机图上的费米 - 哈伯德模型）涉及远距离站点之间的相互作用，但它们是“稀疏”的（每个模式仅与常数 $d$ 个其他模式相互作用）。
• 当前局限性： 在全连接量子比特硬件上，模拟这些稀疏非局域模型通常会在每个 Trotter 步产生 $O(d \log N)$ 的电路深度开销，这是由 JW 弦引起的。与理想费米子硬件（可实现 $O(d)$ 深度）相比，这是一种乘法开销。此前试图消除这种开销的尝试需要大量的辅助量子比特，或者导致初始化成本过高，无法对长时间演化进行有利扩展。

2. 方法论

作者提出了一种新颖的编码方案，引入辅助费米子模式以消除 JW 弦，将问题转化为具有恒定泡利权重项的问题。

A. 系统设置

• 物理模式： $N$ 个物理费米子模式 $a_i$。
• 辅助模式： 每个物理站点 $i$ 增加 $\nu$ 个辅助费米子模式 $b^{(\ell)}_i$（其中 $1 \le \ell \le \nu$）。
• 编码： 对物理模式和辅助模式的组合集应用 JW 排序。

B. 稳定子构建

核心创新在于构建源自辅助马约拉纳算符（$c^{(\ell)}_i, d^{(\ell)}_i$）的稳定子算符 $P^{(\ell)}_{ij}$。

• 定义： $P^{(\ell)}_{ij} = i c^{(\ell)}_i d^{(\ell)}_j$。
• 变换： 原始哈密顿量项 $H_{ij}$ 变换为 $H_{ij} \to H_{ij} P^{(\ell)}_{ij}$。
• 效果： 由于 $P^{(\ell)}_{ij}$ 具有与相互作用项 $H_{ij}$ 相同的 JW 弦结构，将它们相乘会抵消扩展的弦。所得算符仅在 $i, j$ 处的物理量子比特以及涉及的特定辅助量子比特上局域作用，实现了与距离 $|i-j|$ 无关的恒定泡利权重。

C. 态制备与对易性

为了确保物理性质保持不变，系统必须初始化为所有稳定子的 $+1$ 本征空间（$P^{(\ell)}_{ij} |\Psi_{aux}\rangle = |\Psi_{aux}\rangle$）。

• 图着色： 为了确保所有稳定子对易（这是同时制备的要求），作者将相互作用图 $G$ 映射为边着色问题。
• 色指数（$\chi$）： 相互作用图的边被着色，使得任意顶点处的两条边不共享颜色。所需的最小颜色数即为色指数 $\chi$。
• 辅助计数： 每个站点的辅助寄存器数量由 $\nu = \lceil \chi / 2 \rceil$ 确定。两种颜色（互易边集）可以打包到单个辅助层 $\ell$ 中。
• 初始化协议：
  1. 置换： 使用广义费米子置换（深度 $O(\log N)$）重新排序辅助模式，使相互作用的配对相邻。
  2. 有序制备： 使用深度为 $O(\log \nu)$ 的电路为有序配对制备状态。
  3. 总初始化深度： 制备辅助状态的总深度为 $O(\chi \log N + \chi \log \chi)$。关键在于，该状态在变换后哈密顿量的时间演化下是不变的。

3. 主要贡献

1. 加法与乘法开销： 该论文证明，对于稀疏非局域模型，与 JW 弦相关的 $O(\log N)$ 开销可以从每个 Trotter 步的乘法因子转变为加法因子（一次性初始化成本）。
2. 渐近最优性： 对于 Trotter 步数 $M = \Omega(\log N)$ 的长时间模拟，量子比特硬件上的总电路深度变为 $O(M \cdot d \log d)$。这与理想费米子硬件的渐近缩放（$O(M \cdot d)$）相匹配，仅相差常数因子和度数 $d$ 中的对数项。
3. 通用适用性： 该方法适用于一般的稀疏哈密顿量，包括 $d$-正则图上的费米 - 哈伯德模型和稀疏 SYK 模型。

4. 结果与资源缩放

作者提供了详细的资源分析（总结于论文表 I）：

• 辅助量子比特： 需要 $(2\nu + 1)N \approx O(dN)$ 个辅助量子比特（与系统规模呈线性关系）。
• 初始化成本：
  • 深度：$O(\chi \log N + \chi \log \chi)$。
  • 这是一次性成本。
• 时间演化（每 Trotter 步）：
  • 深度：$O(\chi \log \chi)$。
  • 由于对于正则图 $\chi \approx d$，步深度为 $O(d \log d)$。
  • 步深度中没有 $\log N$ 因子。
• $M$ 步的总成本：
  $$\text{Depth}_{total} = O(\chi \log N + M \cdot \chi \log \chi)$$
  在 $M \gg \log N$ 的机制下，成本主要由 $M \cdot \chi \log \chi$ 主导，与费米子硬件的性能相匹配。

5. 意义

• 缩小差距： 这项工作确立了基于量子比特的量子计算机可以模拟稀疏费米子模型，其渐近效率可与理想费米子硬件相媲美，前提是愿意使用线性数量的辅助量子比特。
• 对近期设备的实用性： 通过从重复的时间演化步骤中移除 $O(\log N)$ 依赖，该方法显著减少了长时间动力学模拟所需的电路深度，这对于变分量子本征求解器（VQE）算法或在噪声中等规模量子（NISQ）设备上的实时动力学至关重要。
• 权衡： 该方法用一次性初始化开销和量子比特数量的线性增加（$O(dN)$ 个辅助比特）换取了模拟时间缩放的急剧减少，使其成为模拟复杂非局域量子材料的可行途径。

总之，Irmejs 和 Cirac 提供了一个严格的框架，用于在量子比特硬件上“去弦化”费米子相互作用，证明了在拥有足够的辅助资源的情况下，量子比特编码对于稀疏非局域模型的基本限制可以在渐近意义上被克服。