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

arXiv: 2606.20551v1

论文信息

标题: 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

量子基态制备算法的噪声基准测试：冷却、绝热与 QAOA 的全面对决

研究背景与核心问题

在量子计算和量子模拟领域，制备多体系统的基态是一项核心任务。无论是研究新材料的量子特性、优化复杂系统，还是实现量子化学计算，高效且准确的基态制备都是关键的第一步。目前，主流方法大致分为三类：绝热演化（通过缓慢变化哈密顿量从简单初态过渡到目标基态）、变分量子算法如量子近似优化算法（QAOA，通过优化量子电路参数最小化能量期望值），以及冷却或耗散方案（通过工程化系统与辅助环境的耦合来提取能量）。

然而，所有这些方法在当前的噪声中等规模量子（NISQ）设备上都面临严峻挑战。退相干和操作误差会降低制备态的质量，并限制可达到的能量精度。这自然引出了一个根本性问题：在给定噪声水平和硬件约束下，不同策略的最终表现如何？可达到的能量随噪声强度如何变化？哪种算法在什么条件下最优？

这篇由 Molpeceres 等人完成的论文，正是针对这一核心问题进行了系统性的基准测试。研究团队没有停留在单独分析某种算法的噪声影响，而是建立了一个统一框架，在完全相同的噪声环境和评价指标下，对冷却、绝热和 QAOA 三种代表性算法进行了公平的比较。这种跨算法、系统性的比较方法，填补了现有研究的空白。

模型体系与评价指标

为了能够进行解析推导和大规模数值模拟，研究者选择了一个巧妙且典型的模型系统：一维平移不变的二次费米子哈密顿量，由一个角度参数 θ\theta 控制。这个模型看似简单，却包含了丰富的物理内涵：

  • 当 θ=π/2\theta = \pi/2 时，系统处于 “平凡相”，基态是可以直积的真空态，制备相对容易
  • 当 θ=0\theta = 0 时，系统处于 “拓扑相”，对应 Kitaev 马约拉纳链的拓扑点，基态具有长程纠缠
  • 在 θ=π/4\theta = \pi/4 处，系统发生量子相变，能隙关闭

这种双相结构使得研究者能够测试算法在不同量子相下的表现差异——这是本文的核心洞察之一。噪声模型采用退极化通道，其优势在于与二次哈密顿量的酉演化对易，极大地简化了解析处理。

作为评价指标，研究者采用相对能量而非传统的全局保真度。相对能量定义为：

e=E−EGSEGSe = \frac{E - E_{GS}}{E_{GS}}

其中 EGSE_{GS} 是基态能量，EE 是制备态的能量期望值。这一选择的精妙之处在于：相对能量是内禀的、强度量，不随系统尺寸增大而失效，更符合量子模拟的实际需求。

核心方法：三种算法的噪声标度律

论文最具创新性的贡献之一是推导出了三种算法可达到的相对能量随噪声率 κ\kappa 的标度律。这一分析建立在弱耦合、低噪声等可控近似下，但结果揭示了深刻的物理机制。

冷却算法采用了研究团队此前提出的多频率、随机化时间方案。其基本思想是：在每轮冷却循环中，随机选择耦合时间 τ\tau 和辅助比特频率 Δ\Delta，使辅助比特与系统特定能量模式共振，从而选择性地提取该模式的能量。通过覆盖整个能谱的多个频率和充分长的循环时间，系统逐渐冷却至稳态。

在退极化噪声下，每个冷却子循环的噪声贡献与 κt\kappa t（tt 为平均循环时间）成正比，而冷却速率与 g2tg^2 t（gg 为耦合强度）成正比。在固定 gtgt 的约束下，存在一个最优循环时间平衡这两者。解析推导表明，无论处于平凡相还是拓扑相，冷却算法的相对能量均按 κ1/3\kappa^{1/3} 标度：

ecoolopt≈3C(A2η−2κ)1/3e_{\text{cool}}^{\text{opt}} \approx \frac{3}{C}(A^2\eta^{-2}\kappa)^{1/3}

其中 A,CA, C 是依赖于系统能谱的常数，η=gt\eta = gt。

绝热算法的表现则高度依赖于目标相。在平凡相中，由于不存在能级交叉，非绝热误差按 T−2T^{-2}（TT 为总演化时间）衰减；噪声效应按 κT\kappa T 增长。优化得到相对能量按 κ2/3\kappa^{2/3} 标度：

eprodopt≈3A~1/3κ2/3e_{\text{prod}}^{\text{opt}} \approx 3\tilde{A}^{1/3}\kappa^{2/3}

在拓扑相中，绝热路径必须穿过量子临界点，能隙关闭导致基态附近的模式极易被激发。借助朗道-齐纳理论，非绝热误差按 T−1/2T^{-1/2} 衰减，优化后相对能量按 κ1/3\kappa^{1/3} 标度，但前因子 B~\tilde{B} 随深入拓扑相而增大：

eentopt≈321/3B~2/3κ1/3e_{\text{ent}}^{\text{opt}} \approx \frac{3}{2^{1/3}}\tilde{B}^{2/3}\kappa^{1/3}

这一结果揭示了一个反直觉的结论：仅从标度指数看，拓扑相的绝热算法似乎更抗噪（κ1/3\kappa^{1/3} vs κ2/3\kappa^{2/3}），但由于前因子更大，实际能量反而更高。这提醒我们，在评估算法性能时，不能只看标度指数，前因子同样至关重要。

数值实验：QAOA 与优化冷却的对决

除了解析比较，研究者还进行了公平的数值对比：将 QAOA 与可优化冷却算法在相同数量的自由参数下进行调优。冷却协议允许两个冷却循环，每个循环拥有独立的辅助频率、循环时间和耦合参数；QAOA 则使用 pp 层线路，每层由两个非对易哈密顿量的演化组成。

关键发现：

  • 在平凡相，QAOA 与冷却表现相当，均优于基础绝热算法
  • 在拓扑相，冷却算法明显优于 QAOA。这与对称保护拓扑态的电路深度约束一致：从平凡直积态制备拓扑态需要深度随系统尺寸线性增长的 QAOA 线路，而本文考虑的固定深度线路难以满足
  • 更重要的是，当优化参数受到微小扰动时（模拟实验实现中的参数误差），冷却算法展现出显著更强的鲁棒性。图 2 显示，随着扰动强度从 0 增加到 5%，冷却优于 QAOA 的区域明显扩大

冷却算法的鲁棒性可归因于其吸引子机制：稳态是非平衡动力学的稳定固定点，参数微小变化只会略微移动这个固定点，而不会改变其存在性。QAOA 则缺乏这种稳定性保护。

创新贡献与实践启示

本文的创新点可归纳为以下四个层面：

统一基准框架：建立了在相同噪声模型和评价指标下比较不同基态制备算法的系统方法学，这在已有文献中尚属首次。

解析标度律推导：在弱耦合和低噪声极限下严格推导了冷却和绝热算法的能量−噪声标度关系，揭示了不同量子相下的机制差异。

相依赖性洞察：明确展示了最优算法依赖于目标量子相——平凡相绝热占优，拓扑相冷却占优。这一结论可能推广到更广泛的量子多体系统。

鲁棒性优势验证：通过引入参数扰动，定量证明了冷却算法在实际噪声环境中的工程优势，为实验实现提供了重要指导。

对于量子计算实践者，本文提供了几条有价值的指导：首先，在选择基态制备策略时，必须考虑目标系统的量子相特性；其次，在 NISQ 时代，冷却或耗散算法因其固定点稳定性和参数鲁棒性，可能是比变分算法更务实的选择；最后，在评估算法性能时，应综合考察标度律、前因子和鲁棒性，而非仅看单一指标。

未来展望与局限

本研究也留下了一些悬而未决的问题。最显著的是普适性检验：当前分析基于可二次化的、平移不变的费米子系统，对于更一般的、不可积的相互作用系统，比如哈伯德模型或自旋玻璃等，这些标度律是否仍然成立？量子相变的结构和能谱特征如何影响不同算法的相对优势？

另一个方向是算法改进：是否可以通过优化绝热路径（如采用反绝热驱动）或结合冷却与变分优化来进一步提升性能？机器学习和自动优化技术可能在冷却循环的在线优化中找到用武之地。

此外，本研究未考虑Trotter 化误差和门级噪声模型，这些是数字量子实现中的关键因素。连接连续时间分析与实际量子门操作的误差模型，将是推动这些算法走向实验实现的重要一步。

总的来说，Molpeceres 等人的工作不仅提供了具体的基准测试结果，更重要的是建立了一个可扩展的分析框架。这一框架可以应用于其他模型、噪声通道和算法变体，为量子基态制备的算法选择提供理性的、定量的指导。在量子计算从原理验证走向实用化的关键阶段，这种系统性的性能基准研究具有不可替代的价值。