[論文レビュー] 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

ひとこと要約

本稿では、最小の節-スピン比 α ≈ 1 における一意の満たし可能な割り当て（USA）を有する K-充足可能性（KSAT）問題の最も計算的に難しい例を特定するため、マルチカノニカル集合の技術を用いたモンテカルロ探索手法を導入する。この手法により、エネルギー E=1 における状態密度 Ω₁ が Ω₁ ∝ exp[rN] の形で指数関数的に発散することが判明し、K=2 のとき r=0.348 から K=6 のとき r=0.680 に増加するレート定数 r を示す。これは 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）における状態密度 Ω₁ = g(E=1) を特徴づけ、系のサイズ N に伴うその指数的発散を定量化すること。
• 量子な断続的計算および量子アニーリングの効率をテストするための困難な問題集合を生成すること。
• 既知の複雑性スケーリングを持つ、物理的に妥当で計算的に困難なテストベッドを、量子探索アルゴリズムのためのものとして提供すること。

提案手法

• 統計物理学におけるマルチカノニカル集合モンテカルロ法を用い、一意の満たし可能な割り当て（USA）を有するレアな KSAT 実現をサンプリングする。
• スピン写像 sᵢ = 2Xᵢ − 1 を用いてブール変数をイジングスピンに変換し、 ground state が E=0 に位置するイジングハミルトニアン H_KSAT を構築する。
• 制約条件を課す：(1) 一意の基底状態（g(E=0)=1）、(2) 最小の α = M/N > 1、(3) 第一励起状態における位相空間体積（Ω₁）を最大化。
• 小規模系（2 ≤ N ≤ 18）に対して数値的シミュレーションを実施し、状態密度 g(E) を計算して Ω₁ = g(E=1) を抽出する。
• Ω₁ の N に伴う指数的スケーリング Ω₁ ∝ exp[rN] を分析し、K=2 から 6 までのレート定数 r(K) を決定する。
• 断続的量子計算フレームワークにおける量子分配関数 Z_Q を用いて、量子ギャップ-相関長 Δξ_GAP 及びそのスケーリングを推定する。

実験結果

リサーチクエスチョン

• RQ1 ランダム KSAT 問題において、一意の満たし可能な割り当て（USA）が存在する最小の節-スピン比 α = M/N は何か？ また、これは系のサイズ N にどのように依存するか？
• RQ2 一意の満たし可能な割り当てを有する KSAT 問題において、第一励起エネルギー準位（E=1）における状態密度 Ω₁ = g(E=1) は、系のサイズ N にどのように依存するか？
• RQ3 K=2 から 6 までの指数的スケーリング Ω₁ ∝ exp[rN] のレート定数 r(K) は何か？ また、K が増加するにつれて ln 2 に近づく様子はいかなるか？
• RQ4 困難な KSAT エンsemble における量子ギャップ-相関長 Δξ_GAP は、系のサイズ N にどのように依存するか？ これは量子アニーリングの実行時間にどのような意味を持つか？
• RQ5 E=1 において最大の Ω₁ を有する困難な KSAT 問題の集合は、量子探索アルゴリズムのための堅牢なベンチマークとして機能できるか？

主な発見

• エネルギー E=1 における状態密度 Ω₁ = g(E=1) は、Ω₁ ∝ exp[rN] の形で系のサイズ N に伴い指数関数的に発散し、K=2 のとき r=0.348 から K=6 のとき r=0.680 に増加する。
• K=2 から 6 までのレート定数 r(K) は、ln 2 ≈ 0.693 に近づく。これは、K が大きい場合には、唯一の基底状態（E=0）を除き、ほとんどすべてのスピン配置が E=1 に位置することを示唆する。
• K=3 の場合、量子ギャップ-相関長 Δξ_GAP は Δξ_GAP ∝ exp[r_Q N] とスケーリングし、r_Q ≈ 0.60(1) となる。これは古典的レート r(K=3)=0.569(1) に近く、強い量子的複雑性スケーリングを示唆する。
• 最小の α ≈ 1 において一意の満たし可能な割り当てを有する困難な KSAT 問題の集合は、ランダム KSAT エンsemble において指数的にまれであるため、量子探索アルゴリズムのテストに最適である。
• Ω₁ の指数的特異性は、古典的探索時間 τ_classical ∝ exp[rN] に比例し、量子アニーリングの実行時間は Δξ_GAP の二乗に比例する可能性があるため、指数的スローダウンを引き起こす。
• 高次スピン相互作用を含むこの困難な問題集合でさえ、2スピン相互作用および磁場項のみを用いた最大独立集合（MIS）形式に変換可能であり、その際も難易度は保たれる。