噪声环境下基态制备量子算法的基准测试

论文信息

标题: Benchmark of quantum algorithms for ground state preparation in the presence of noise

作者: Daniel Molpeceres, Sirui Lu, J. Ignacio Cirac, et al.

发布日期: 2026-06-18

arXiv ID: 2606.20551v1

PDF 链接: 下载 PDF

3 分钟速览

研究问题：在噪声环境下，冷却（cooling）、绝热演化（adiabatic）和变分优化（QAOA）三类量子算法，哪一种能以最低的相对能量制备基态？
核心方法：以一类严格可解的二次费米子模型为试验平台，解析推导相对能量 $e$ 随噪声率 $\kappa$ 的标度律，并用数值模拟在相同噪声模型下公平比较三种算法。
关键结果：最优策略依赖于量子相：在平凡相绝热算法更优，在拓扑相多频率冷却算法更具竞争力或更优（见表 I、图 2）。
主要局限：模型始终是自由费米子，未涉及相互作用或超越可积性的系统；数值优化使用的 BFGS 算法未必达到全局最优。
适合读者：从事量子模拟、量子算法设计与含噪量子计算研究，希望理解算法选择与物理相之间关系的研究者。

论文背景和研究动机

基态制备是量子模拟与量子计算的核心任务之一。现有的制备策略可大致归为三类：绝热演化、以 QAOA 为代表的变分算法，以及通过耦合冷辅助浴实现能量抽取的冷却算法或耗散态制备方案。这些算法在无噪声的理想模型中各有优势，但在当前实验条件下，退相干与参数不完美等噪声因素不可避免，直接限制了可达到的态制备品质。

论文试图回答一个实用而根本的问题：在给定噪声模型和硬件约束下，不同算法的最终表现如何随噪声强度变化？哪种策略在特定物理条件下最具优势？已有工作多聚焦于单一算法的噪声容忍度，缺乏不同算法在统一噪声框架下的横向比较。该文由此切入，建立了一套可扩展的基准测试方案，在保持足够解析可控性的同时，揭示算法性能与系统量子相之间的深刻关联。

核心方法和技术细节

论文基于一个平移不变的二次费米子哈密顿量 $H_S(\theta)$ ，参数 $\theta \in [0, \pi/2]$ 控制模型从拓扑相到平凡相的过渡，并在 $\theta = \pi/4$ 处经历量子相变。通过傅里叶变换，多体动力学分解为独立的 $(k, -k)$ 模式对，每一对进一步限定在偶宇称子空间中，形成一组等效的自旋 -1/2 系统。这一特性使解析推导和中等规模数值模拟同时可行。

噪声统一采用作用在每个模式上的去极化信道，噪声率为 $\kappa$ 。对于冷却算法，解析工作集中在弱耦合极限 $(gt)^2 \ll 1$ 、冷却极限 $\epsilon_m t \gg 1$ 和低噪声极限 $\kappa t \ll 1$ 下完成。通过戴森展开和随机化循环时间 $\tau$ 与浴频率 $\Delta$ ，文章推导出有效冷却率和加热率，并在此基础上建立含噪稳态相对能量的解析表达式。绝热算法则分别针对无相变跨越的平凡相（利用绝热微扰论得到 $T^{-2}$ 修正）和穿过临界点的拓扑相（借助朗道-齐纳理论得到 $T^{-1/2}$ 修正）加以分析。两种情形下，非绝热误差与噪声引起的指数衰减 $e^{-2\kappa T}$ 竞争，形成最优演化时间。

变分算法 QAOA 采用 $p$ 层含参量子电路，每层由两个非对易的二次哈密顿量生成的幺正演化构成。电路参数通过 BFGS 算法优化，且总参数量与冷却算法的自由参数数目保持一致，以保证公平比较。所有算法均采用冷却算法在 $\theta = \pi/2$ 处产生的混合态 $\rho_{\mathrm{cool}}(\pi/2)$ 作为初态，排除了理想初态带来的不公平优势。

创新点和贡献

论文的突出贡献在于建立了一个可关联解析推导与数值验证的比较框架。该框架的主要创新体现在三个方面。首先，文章从冷却算法中解析地导出了相对能量对噪声率的最优标度律 $e \propto \kappa^{1/3}$ ，并与绝热算法在平凡相和拓扑相中分别获得的 $\kappa^{2/3}$ 和 $\kappa^{1/3}$ 标度进行对比，这在同类研究中尚属首次。表格 I 清晰总结了这些标度，并指出冷却算法的标度与量子相无关，而绝热算法的性能强烈依赖于是否穿过量子相变。

其次，文章系统地分析了算法性能的相依赖性。在平凡相，绝热算法凭借较小的前置因子和 $\kappa^{2/3}$ 的较优标度，在较低的噪声水平下占据优势；在拓扑相，冷却算法不仅在最优参数下给出更低的相对能量，而且对参数偏差表现出更强的鲁棒性（见图 2）。这种相敏感的行为与绝热路径是否跨越最小能隙闭合点直接相关，揭示了基态制备策略选择的深层物理依据。

第三，论文揭示并量化了冷却算法的鲁棒性优势。在扰动幅度 $\delta = 0.05$ 时，冷却算法相对于 QAOA 的领先区域显著扩大（图 2, c 和 f）。这一性质源于冷却算法的固定点结构：稳态是信道的吸引子，小参数扰动仅轻微移动吸引点，而 QAOA 通过单一幺正操作生成目标态，能量对参数扰动更为敏感。

实验结果分析

解析推导与数值模拟高度吻合。图 4 展示了多频率冷却与理想及现实条件下的绝热算法在若干目标点 $\theta_f$ 处的相对能量最优值。 $\theta_f = \pi/3$ （平凡相）时，绝热算法的能量始终低于冷却算法；当 $\theta_f$ 深入拓扑相，两条曲线逐渐接近并最终变为冷却占优。在极低噪声区域，由于 $k = N/2$ 模式在临界点能隙完全闭合而不可避免地激发，绝热曲线出现与系统尺寸相关的地板效应，该效应在热力学极限下消失（图 5）。

图 6 比较了优化冷却与 QAOA 在参数数量同为 20 时的整体表现。冷却在 $\theta < \pi/4$ 的拓扑相中普遍得到更低的相对能量，尤其在 $\theta \to 0$ （基塔耶夫链极限）时优势显著，这与 QAOA 电路深度固定时难以制备对称性保护拓扑态的已知约束一致。论文没有给出两种算法在单次优化运行中的迭代次数或收敛时间，也没有提供不同随机初始点下能量分布的方差或散度数据。

实践建议

对于当前量子硬件实现基态制备的工程师和实验物理学家，本研究提供了几项值得关注的指导性结论。第一，算法选择应根据目标哈密顿量的相图特性来决策：如果系统位于或接近拓扑相，优先考虑冷却类协议可能获得更高的态制备精度；如果目标处于充分有能隙的平凡相，简单的线性绝热演化即可提供有竞争力的性能。第二，在实验参数标定存在一定波动或系统漂移的现实条件下，冷却算法因其固定点结构的吸引性，对参数不完美的容忍能力明显优于 QAOA，这对提高实验可重复性和降低校准复杂度是有利信号。第三，在工程实现中，冷却协议所需的浴频率数目 $R$ 应根据噪声水平权衡：在噪声较强时，可适当减少频率数目以降低总运行时间中噪声的累积影响，论文在附录 A 开头 4 小节中讨论了这种多频到少频的过渡及其对能量标度律的修正。这些结论虽然来自一个特解的可积模型，但所揭示的相依赖性和鲁棒性趋势对更一般的量子多体系统具有一定的启发性。