[论文解读] Monte Carlo Search for Very Hard KSAT Realizations for Use in Quantum Annealing

Neuhaus Thomas|arXiv (Cornell University)|Dec 17, 2014

Quantum Computing Algorithms and Architecture参考文献 3被引用 5

一句话总结

本文提出一种基于多Canonical系综技术的蒙特卡洛搜索方法，用于在最小的子句-自旋比（α ≈ 1）下识别具有唯一满足赋值（USA）的K-可满足性（KSAT）问题中最难的实例。研究发现，能量E=1处的态密度Ω₁随系统尺寸N呈指数发散，即Ω₁ ∝ exp[rN]，其中速率常数r从K=2时的0.348增至K=6时的0.680，趋近于ln 2，为测试量子退火效率提供了基准系综。

ABSTRACT

Using powerful Multicanonical Ensemble Monte Carlo methods from statistical physics we explore the realization space of random K satisfiability (KSAT) in search for computational hard problems, most likely the 'hardest problems'. We search for realizations with unique satisfying assignments (USA) at ratio of clause to spin number $α=M/N$ that is minimal. USA realizations are found for $α$-values that approach $α=1$ from above with increasing number of spins $N$. We consider small spin numbers in $2 \le N \le 18$. The ensemble mean exhibits very special properties. We find that the density of states of the first excited state with energy one $Ω_1=g(E=1)$ is consistent with an exponential divergence in $N$: $Ω_1 \propto { m exp} [+rN]$. The rate constants for $K=2,3,4,5$ and $K=6$ of KSAT with USA realizations at $α=1$ are determined numerically to be in the interval $r=0.348$ at $K=2$ and $r=0.680$ at $K=6$. These approach the unstructured search value ${ m ln}2$ with increasing $K$. Our ensemble of hard problems is expected to provide a test bed for studies of quantum searches with Hamiltonians that have the form of general Ising models.

研究动机与目标

• 识别随机K-可满足性（KSAT）问题中具有唯一满足赋值（USA）的最困难计算实例。
• 探索在子句-自旋比α = M/N从上方趋近于1时KSAT的实现空间。
• 表征第一激发能级（E=1）处的态密度Ω₁，并量化其随系统尺寸N的指数发散特性。
• 构建一个用于测试量子绝热计算与量子退火效率的困难问题系综。
• 为具有已知复杂度标度的量子搜索算法提供一个物理上有依据、计算上具有挑战性的测试平台。

提出的方法

• 采用统计物理中的多Canonical系综蒙特卡洛方法，对具有唯一满足赋值（USA）的罕见、困难KSAT实现进行采样。
• 通过自旋映射sᵢ = 2Xᵢ − 1将布尔变量转换为伊辛自旋，从而构建伊辛哈密顿量H_KSAT，其基态位于E=0。
• 施加约束条件：(1) 唯一基态（g(E=0)=1），(2) 最小α = M/N > 1，(3) E=1处相空间体积最大（Ω₁）。
• 对小系统（2 ≤ N ≤ 18）进行数值模拟，计算态密度g(E)，并提取Ω₁ = g(E=1)。
• 分析Ω₁随Ω₁ ∝ exp[rN]的指数标度行为，确定K=2至6的速率常数r(K)。
• 在绝热量子计算框架中使用量子配分函数Z_Q，估算量子能隙-相关长度Δξ_GAP及其标度行为。

实验结果

研究问题

• RQ1 在随机KSAT问题中，唯一满足赋值（USA）存在的最小子句-自旋比α = M/N是多少？其随系统尺寸N如何变化？
• RQ2 对于具有USA的KSAT问题，第一激发能级（E=1）处的态密度Ω₁ = g(E=1)如何随系统尺寸N变化？
• RQ3 在指数标度Ω₁ ∝ exp[rN]中，K=2至6的速率常数r(K)是多少？它们如何趋近于无结构搜索极限ln 2？
• RQ4 在困难KSAT系综中，量子能隙-相关长度Δξ_GAP如何随系统尺寸N变化？这对量子退火运行时间有何含义？
• RQ5 能否将E=1处Ω₁最大的困难KSAT问题系综作为量子搜索算法的稳健基准？

主要发现

• 能量E=1处的态密度Ω₁ = g(E=1)随系统尺寸N呈指数发散，即Ω₁ ∝ exp[rN]，其中速率常数r从K=2时的0.348增至K=6时的0.680。
• K=2至6的速率常数r(K)趋近于无结构搜索极限ln 2 ≈ 0.693，表明当K较大时，几乎所有自旋构型都处于E=1，仅有一个基态位于E=0。
• 对于K=3，量子能隙-相关长度Δξ_GAP的标度行为为Δξ_GAP ∝ exp[r_Q N]，其中r_Q ≈ 0.60(1)，接近经典速率r(K=3)=0.569(1)，表明存在强烈的量子复杂度标度行为。
• 在最小α ≈ 1处具有唯一满足赋值的困难KSAT问题系综在随机KSAT系综中呈指数稀少，因此特别适合作为测试量子搜索算法的基准。
• Ω₁的指数奇异性意味着经典搜索时间τ_classical ∝ exp[rN]，而量子退火运行时间可能与Δξ_GAP的平方成正比，导致指数级变慢。
• 尽管该困难问题系综由高阶自旋相互作用构建，但可转换为仅含两体自旋耦合与磁质场耦合的极大独立集（MIS）形式，且保持其困难特性。