[论文解读] Simulations of the adiabatic quantum optimization for the Set Partition Problem

Vadim Smelyanskiy, U. von Toussaint+1 more|ArXiv.org|Dec 22, 2001

Quantum Computing Algorithms and Architecture参考文献 14被引用 10

一句话总结

本文通过求解含时与定态薛定谔方程的数值模拟，研究了对NP完全问题集合划分问题（SPP）的绝热量子优化。引入对数尺度的代价函数以减少不同能量值的数量，并发现该算法的复杂度随问题规模呈指数增长，与先前对其他NP完全问题的研究结果相矛盾。

ABSTRACT

We analyze the complexity of the quantum optimization algorithm based on adiabatic evolution for the set partition problem. We introduce a cost function defined on a logarithmic scale of the partition residues so that the total number of values of the cost function is of the order of the problem size. We simulate the behavior of the algorithm by numerical solution of the time-dependent Schroedinger equation as well as the stationary equation for the adiabatic eigenvalues. The numerical results for the time-dependent quantum evolution indicate that the complexity of the algorithm scales exponentially with the problem size.This result appears to contradict the recent numerical results for complexity of quantum adiabatic algorithm applied to a different NP-complete problem (Farhi et al, Science 292, p.472 (2001)).

研究动机与目标

• 分析绝热量子优化在NP完全集合划分问题（SPP）中的复杂度。
• 研究绝热量子算法是否能像某些NP完全问题所假设的那样，实现SPP的多项式时间性能。
• 检验能量景观结构（尤其是回避交叉的分布）决定算法复杂度的假设。
• 探讨对划分余数进行粗粒化处理如何影响绝热演化过程的动力学与复杂度。
• 比较含时与定态薛定谔方程方法在分析SPP绝热演化动力学时的表现。

提出的方法

• 使用问题哈密顿量 $ H_P $ 将SPP表述为一个量子优化问题，该哈密顿量编码了划分余数的能量函数。
• 引入对数尺度的代价函数，将不同能量值的数量限制在问题规模量级，从而降低有效希尔伯特空间的复杂度。
• 对含时薛定谔方程进行数值模拟，以追踪绝热演化过程中量子态的演化。
• 求解定态薛定谔方程，计算绝热本征值，并分析小规模 $ n \leq 10 $ 时的回避能级交叉。
• 利用微扰理论，通过最小能隙 $ \Delta g_{\text{min}} $ 和矩阵元 $ V_{0k} $ 估算非绝热跃迁概率。
• 分析概率分布 $ |\langle z|\psi(T)\rangle|^2 $，以评估在不同代价函数值下最终态的局域化与结构特征。

实验结果

研究问题

• RQ1 集合划分问题的绝热量子优化在问题规模增大时，运行时间是呈现多项式还是指数增长？
• RQ2 选择对数尺度代价函数如何影响能量景观结构及算法的复杂度？
• RQ3 回避交叉与最小能隙在决定绝热算法运行时间复杂度方面起什么作用？
• RQ4 尽管最小能隙较小，为何算法仍无法实现高效性能？这是否与波函数中缺乏结构有关？
• RQ5 含时与定态薛定谔方程方法在捕捉SPP绝热演化动力学方面，其结果有何异同？

主要发现

• 集合划分问题的绝热量子优化算法复杂度随问题规模 $ n $ 呈指数增长，该结论由含时薛定谔方程的数值模拟结果支持。
• 对于 $ n $ 从11到17的中位运行时间复杂度，其拟合结果符合指数函数，最小二乘法拟合显示出明显的指数增长趋势。
• 概率分布 $ |\langle z|\psi(T)\rangle|^2 $ 形成与不同代价函数值相对应的‘台阶’结构，且每个台阶内部无精细结构。
• 在相同代价函数值下，波函数中缺乏精细结构，被确定为导致指数复杂度的主要原因。
• 对于 $ n=10 $，定态薛定谔方程分析揭示了涉及多个激发态的复杂回避交叉结构，偏离了标准的两能级模型。
• 基态与相关激发态之间的最小能隙显著大于其他激发态之间的间距，表明在回避交叉区域存在非平凡的多体效应。