Worldline-Susceptibility Scheduling for Quantum Annealing Beyond Local-Adiabatic Evolution

技术摘要：用于量子退火的世界线磁化率调度方案

问题陈述
量子退火（QA）的性能高度依赖于退火调度，该调度决定了总运行时间如何在演化路径上进行分配。虽然 Roland–Cerf 局部绝热调度在理论上是最优的——即将时间分配与瞬时能隙（$\Delta(s)^{-2}$）成反比——但在处理大规模优化问题时，这在实际中是不可行的。构建此类调度需要对瞬时低能谱的完整了解，这需要对指数级规模的哈密顿量进行重复的对角化运算。因此，迫切需要一种计算成本低廉的替代方案，能够在无需显式谱信息的情况下识别退火过程中的关键区域。

方法论
作者提出了一种基于模拟量子退火（SQA）过程中测量的世界线磁化率（$\chi_m$）的替代调度框架。该方法流程如下：

1. SQA 模拟： 量子退火过程通过 Suzuki–Trotter 分解被映射为一个等效的经典统计力学系统。这使得横场伊辛模型转化为一个具有 $M$ 个耦合副本（Trotter 切片）的 $(d+1)$ 维经典系统。
2. 观测量测量： 使用经典蒙特卡洛技术对平衡构型进行采样。从中可以计算出世界线磁化强度 $m$，并通过涨落-耗散定理导出其磁化率：$\chi_m = NM(\langle m^2 \rangle - \langle |m| \rangle^2)$。
3. 替代方案构建： 测得的磁化率 $\chi_m(s)$ 作为能隙平方倒数的代理指标。定义时间分配权重 $w(s) = \chi_m(s) + \chi_0$（其中 $\chi_0$ 为正则化常数）。通过对该权重进行积分构建累积时间分配函数 $\tau(s)$，进而通过反函数获得最终的调度方案 $s(t)$。
4. 验证： 该框架已在开源软件 Qanneal 中实现，通过针对 Sherrington–Kirkler 随机自旋玻璃实例的精确对角化（ED）结果来验证其替代调度方案。通过求解含时薛定谔方程，以线性调度、精确的 Roland–Cerf 调度以及所提出的替代调度方案为基准，通过最终基态概率（$P_{GS}$）来量化性能。

关键结果
研究表明，对于很大一部分 SK 实例，世界线磁化率调度方案的表现始终优于传统的线性退火，甚至优于直觉上的精确 Roland–Cerf 调度。

• 关键区域识别： 世界线磁化率的峰值能够可靠地识别出能隙极小值所在的宽泛区域，尽管它并不能定量地重现详细的能隙剖面。
• 卓越的有限时间性能： 对于代表性实例（$n=10$ 和 $n=12$），在各种总退火时间下，替代调度方案均能获得比精确 Roland–Cerf 调度更高的基态概率。
• 失效模式发现： 作者识别出了精确局部绝热调度的两种不同有限时间失效模式，这解释了其表现不佳的原因：
  1. 边界能隙陷阱（Boundary-Gap Trap）： 当最小能隙出现在退火末端（$s \approx 1$）时，Roland–Cerfeld 调度会将大部分运行时间分配给横场驱动器已基本消失的区域。这使得增加的时间在动力学上是无效的，因为此时量子涨落已被抑制。基于磁化率的调度通过将时间集中在量子涨落依然活跃的更宽阔的交叉区域，从而避免了这一问题。
  2. 振荡不稳定性（Oscillatory Instability）： 当最小能隙位于中间位置时，精确调度会将运行时间过于集中地聚集在能隙附近，从而诱发相干多级干涉。这导致成功概率随总退火时间呈现非单调的振荡依赖关系，而较平滑的磁化率调度方案则避免了这一现象。
• 统计鲁棒性： 在涵盖不同无序实现和系统规模（$n=10$ 至 $n=20$）的大规模数值实验中，研究证实了替代调度方案的优越性是一个普遍特征，而非特定实例的特例。随着系统规模的增加，边界能隙实例出现的频率上升，导致精确 Roland–Cerf 调度的性能相对于替代调度方案进一步下降。

意义与主张
本文认为，在现实的有限时间条件下，调度的鲁棒性比精确重现瞬时能隙剖面更为重要。作者主张：

• 低成本的平衡态观测量（特别是来自 SQA 的世界线磁化率）可以作为基于能隙策略的有效替代方案。
• 精确的局部绝热条件虽然在渐近极限下是理论最优的，但在有限时间机制下，由于“边界能隙陷阱”和“振荡不稳定性”的存在，可能会产生适得其反的效果。
• 所提出的方法论提供了一条可扩展的调度优化路径，绕过了精确谱计算在计算上的不可行性，使其能够应用于超出精确对角化能力的更大规模问题。

该工作已在开源 Qanneal 框架中实现，该框架为 SQA、磁化率测量和调度构建提供了一个统一的平台，从而促进了量子退火动力学领域的可重复研究。