近似影子哈密顿模拟的高效算法

论文信息

标题: An efficient algorithm for approximate shadow Hamiltonian simulation

作者: Abhijit Chakraborty, Bharath Sambasivam, Karunya Shirali, et al.

发布日期: 2026-07-13

arXiv ID: 2607.11882v1

PDF 链接: 下载 PDF

3 分钟速览

研究问题：这篇论文要解决在相互作用多体系统中使用影子哈密顿量模拟方法时，所需算符代数呈指数增长、导致量子比特资源开销过大的瓶颈问题。原方法在无相互作用系统中效率很高，但一到有相互作用的物理系统就失去优势。
核心方法：作者提出三种「算符剪枝」策略——基于预定义基的图剪枝、基于算子 Lanczos 的 Krylov 基剪枝、以及将两者结合的混合方法。这三种方法系统地识别出对目标观测量最重要的算符子空间，从而构建大大缩小规模的近似影子哈密顿量。
关键结果：在一维混合场 Ising 模型中，仅用 10 个量子比特的影子寄存器就可跟踪 100 个物理比特系统的磁化动力学（见图 2）。二维系统中也观察到类似的资源节省（见表 1）。
主要局限：方法有效性依赖于哈密顿量中存在小参数（如较弱的横场）。当小参数接近 $\mathcal{O}(1)$ 量级时，剪枝效率急剧下降。纯 Krylov 基方法在 Pauli 基下需要存储指数增长的算符展开系数，限制了其在中等规模以上系统的可扩展性（见论文 Numerical results 倒数第二段）。
适合读者：从事量子算法设计、量子模拟或量子多体物理研究的科研人员。对量子计算在凝聚态物理、量子化学等领域的实际应用感兴趣的人也值得一读。

论文背景和研究动机

实时模拟多体哈密顿量的动力学是量子计算的核心应用方向之一，涉及高能物理、凝聚态物理、量子混沌和量子化学等多个领域。传统的量子态演化算法（如乘积公式、量子信号处理、变分方法等）通常需要用 $N$ 个物理量子比特直接编码 $N$ 量子比特系统的态，这在多数情况下是必要的，但也意味着量子资源开销与系统规模呈线性关系。

2025 年，Somma 等人在 Nature Communications 上提出的影子哈密顿量模拟（shadow Hamiltonian simulation）提供了一种新思路：与其跟踪完整的量子态，不如仅跟踪一组关心的可观测量的期望值演化。其核心思想是将观测量在哈密顿量对易作用下生成的算符代数的期望值向量——即「影子态」——编码到量子比特上，然后在影子哈密顿量的驱动下进行演化。

这一方法在无相互作用系统（如自由费米子模型）中表现出色，因为算符代数的维度随系统大小仅呈多项式增长，影子态编码所需的量子比特数仅为 $\mathcal{O}(\log_2 N)$ 。然而对于有相互作用的系统，算符代数的维度通常随 $N$ 呈指数增长，最坏情况下需要 $2N$ 个量子比特来精确编码影子态——这完全抵消了影子模拟潜在的优势。

本文正是在这一背景下，针对「影子模拟如何处理有相互作用系统」这一关键难题提供了系统性的解决方案。作者的核心洞察是：并非所有出现在算符代数中的算符都同等重要；通过有控制的近似，可以仅保留那些对目标观测量贡献最大的算符，从而大幅压缩所需代数维度。

核心方法和技术细节

论文的核心框架是对影子哈密顿量模拟进行「有控制的近似」。目标观测量 $\bar{O}$ 的海森堡演化可以写成嵌套对易子的级数展开：

\bar{O}(t) = \sum_{j=0}^{\infty} \frac{(it)^j}{j!} \mathcal{A}_H^j(\bar{O})

其中 $\mathcal{A}_H(\cdot) := [H, \cdot]$ 是对易超算符。精确的影子模拟要求包含所有出现在这些嵌套对易子中的算符，而近似方法则将其截断为最重要的一部分。

作者提出了三种剪枝方案：

预定义基方案（Predefined basis pruning）

该方法将算符代数的生成过程建模为一个带权重的有向层状图。图的根层节点为目标观测量在预定义基（本文中为 Pauli 基）下的展开分量。第 $i+1$ 层节点是 $\mathcal{A}_H$ 作用在第 $i$ 层节点上产生的新基算符。每条边的权重由对易子系数定义（即式 8），每个节点的累积权重定义为从根层到该节点的最大路径权重之积（式 9）。在生成图的过程中，凡累积权重低于预设阈值 $\epsilon$ 的节点及其连边都被剪除。

值得注意的是，对于摄动哈密顿量 $H = H_0 + \xi H_1$ ，若取 $\epsilon = \xi^M$ ，该剪枝方案等价于对 Dyson 级数做 $M$ 阶截断（见 Methods D）。作者由此推导出在时间主导区和精度主导区下 $M$ 的渐近标度（式 32），尽管该界限在实践中相当宽松（见图 8）。

Krylov 基方案（Krylov basis pruning）

该方法采用算子 Lanczos 迭代，从根观测量出发依次生成正交化的 Krylov 向量。每一步新算符先由 $\mathcal{A}_H$ 作用于上一个 Krylov 向量得到，再与前序向量做 Gram-Schmidt 正交化。与常规标准不同，作者采用的截断判据是乘性的： $w = \prod_{j=1}^k \|O_j^\perp\| \ge \epsilon$ （式 13），其中 $O_j^\perp$ 是每次正交化后的垂直分量。这一设计使得沿整条 Krylov 轨道的相对重要性得以累积追踪，比仅看单步范数的传统判据更为精细。

然而，该方案的一个致命弱点是：尽管 Krylov 基本身维度可能很小，但在 Pauli 基下展开每个 Krylov 向量时，可能包含多达 $4^N$ 项的不同 Pauli 串（见 Numerical results 讨论），导致经典存储成本和阴影态初始化成本都呈指数增长。

混合方案（Hybrid approach）

混合方案正是为了克服 Krylov 方案的上述局限而设计的。它先用预定义基方案的图剪枝确定一个规模受控的算符子空间 $\mathbf{S}_{\bar{O}}$ ，然后在该子空间内运行 Krylov 算法。这样做带来双重收益：(1) 每个 Krylov 向量最多只包含 $|\mathbf{S}_{\bar{O}}|$ 个 Pauli 串元素，彻底解决了存储爆炸问题；(2) 在相同模拟精度下，所需影子寄存器规模通常比纯预定义基方案更小（见表 1）。

创新点和贡献

首次将有相互作用的哈密顿量纳入影子模拟的可操作范畴：此前影子模拟被认为仅在无相互作用系统中有效，本文证明了通过系统的算符剪枝，甚至在一维、二维有相互作用的自旋系统中都能取得实质性的量子资源节省。
三种互补的剪枝策略的完整框架：预定义基方案提供了清晰的图表示和摄动论连接，Krylov 方案沿动力学轨道自适应选择最优基，混合方案则在可扩展性和精度间取得平衡。三种方法可根据实际问题的特点灵活选用或组合。
超越局域观测量的适用性：论文展示了剪枝方案不仅对单点磁化这种局域观测量有效，也能精确捕获两点自旋流关联函数（图 4）和四点 OTOC（图 5）的动力学。这在理论和应用层面都具有重要意义，因为自旋输运和量子混沌的诊断正是凝聚态物理和高能物理的前沿话题。
经典-量子协同的工作流程：剪枝过程完全在经典计算机上预处理完成，最终只向量子计算机提交一个维度大幅缩减的影子哈密顿量。这种分工在近期量子设备的噪声和比特数约束下尤为实用。

实验结果分析

一维混合场 Ising 模型

论文在一维 MFIM（ $J=1$ , $h_x=1$ , $h_z=0.3$ ）上系统测试了三种剪枝方案的性能。

对于预定义基方案（图 2a），在不同剪枝阈值 $\epsilon = h_z^m$ 下，影子寄存器规模均随物理系统大小 $N$ 增加而趋于饱和。例如，阈值 $\epsilon = h_z^3$ 时只需 10 个影子量子比特即可表示影子哈密顿量，而精确影子模拟则需要约 24 个比特（红色虚线）。图 2b 展示了 $N=12$ 系统上磁化的动力学演化，较低阈值 $\epsilon = h_z^4, h_z^5$ 的结果在 $t \le 20$ 范围内几乎与精确对角化一致。

Krylov 基方案（图 2c-d）在 $\epsilon = 10^{-12}$ 下仅需 7 个量子比特即可模拟 $N=8$ 系统，但如前所述，该方案因 Pauli 展开系数的爆炸式增长而在更大系统中无法直接使用（论文明确说明 $N>8$ 的所有案例均未采用纯 Krylov 方法）。

混合方案（图 2e-f）成为最具优势的选择：对于 $\epsilon = (h_z^3, 10^{-10})$ （括号中第一个阈值为预定义基剪枝，第二个为 Krylov 剪枝）， $N=12$ 系统仅需 7 个影子量子比特，且磁化精度优于同等比特数的纯预定义基方案。

二维 MFIM 与扩展测试

在二维正方格点 MFIM 上（图 3）， $h_z = 0.2$ 且阈值 $\epsilon = h_z^2$ 时，一个 $4\times4$ （16 比特）系统的平均磁化可由 14 个影子量子比特精确模拟，混合方案更可进一步降至 7 个量子比特（图 3b）。

对于非局域观测量，论文考察了两种重要场景：(1) 海森堡 XXZ 模型中的自旋流自关联函数（图 4），在弹道输运区 $J_z = 0.05$ 下用最多 10 个影子量子比特精确捕获了扩散系数 $\alpha=0.96$ （理论预期 $\alpha=1$ ）；(2) 五点 LFXXZ 模型的 OTOC（图 5），成功再现了早期指数的快速上升和后续的振荡，且混合方案中每个 Krylov 向量内的 Pauli 串数显著少于纯 Krylov 方案（内插图）。

局限与待解决问题

尽管剪枝策略在数值实验中表现优异，但论文坦诚指出了若干根本性局限。

首先，方法的有效性强烈依赖于哈密顿量中存在小参数。当 $h_z \to 1.0$ 时，完整的 Pauli 代数几乎都需要保留，近似影子模拟的优势荡然无存（见 Numerical results 第 4 段）。类似地，如果目标观测量是接近全满 Hamming 权重的 Pauli 串（如 $Z^{\otimes N}$ ），剪枝后集合的维度仍会以 $2^n$ （ $n > N$ ）的速率增长（见 Supplementary Information S.VI）。

其次，纯 Krylov 方法在 Pauli 基下的可扩展性存在本质障碍。每一轮 Krylov 迭代都可能将算符在 Pauli 基中的支持扩散到指数多分量上。虽然混合方案在本文的一维和二维例子中成功规避了这一问题，但对于更复杂的哈密顿量（如量子化学中的费米子哈密顿量），预定义基剪枝本身能否有效压縮代数空间仍有待研究。

第三，论文中所有数值模拟均使用精确对角化作为基准，且最大测试规模仅到 16 个物理比特（二维）或 100 个物理比特（一维，但此时无法做精确对角化验证）。对于影子哈密顿量在真实量子硬件上的实现——包括影子态的制备代价、影子哈密顿量的局部性（locality）和门复杂度分析，论文均悬而未决，留待未来工作。

最后，论文未讨论在剪枝后的影子哈密顿量维度不是 $2^n$ 形式时，如何高效地将其嵌入到物理量子比特上进行模拟。这在工程层面是一个不可回避的问题。