[論文レビュー] A PTAS for the Classical Ising Spin Glass Problem on the Chimera Graph Structure

Rishi Saket|arXiv (Cornell University)|Jun 20, 2013

Markov Chains and Monte Carlo Methods参考文献 4被引用数 9

ひとこと要約

本稿では、平面グリッドにおけるBansalら（2009年）の手法を、量子アニーリングで用いられる非平面的構造であるチメラグラフにおけるイジングスピンガラス問題に適応することで、多項式時間近似スキーム（PTAS）を提示する。主な結果として、任意のε>0に対して時間O(n·2^(32/ε))で(1−ε)-近似が達成可能であり、自明な下界が最適解の定数倍（≤17.26）以内にあることが、グロテンジックの不等式およびK₄,₄部分グラフ解析を用いて証明された。

ABSTRACT

We present a polynomial time approximation scheme (PTAS) for the minimum value of the classical Ising Hamiltonian with linear terms on the Chimera graph structure as defined in the recent work of McGeoch and Wang. The result follows from a direct application of the techniques used by Bansal, Bravyi and Terhal who gave a PTAS for the same problem on planar and, in particular, grid graphs. We also show that on Chimera graphs, the trivial lower bound is within a constant factor of the optimum.

研究の動機と目的

• 量子アニーリングで用いられる非平面的チメラグラフ構造へ、平面グラフにおけるイジングスピンガラス問題のPTASを拡張すること。
• チメラグラフ上でのイジングハミルトニアンの最小化に対して、多項式時間近似スキーム（PTAS）の存在を確立すること。
• 地面状態エネルギーの自明な下界が真の最適解の定数倍以内にあることを証明し、近似保証を向上させること。
• Bansalら（2009年）の平面グリッドにおける分割および動的計画法の技術を、チメラグラフのレイヤー構造と多成分構造に適応すること。

提案手法

• E₀またはE₁からのわずかなエッジの削除によって、チメラグラフを定数幅のストリップに分解し、再帰的分割を可能にする。
• T=2/ε個のパーティションに対して確率的丸めを適用することで、最小ハミルトニアン値の(1−ε)-近似を達成する。
• 幅8T以下のストリップグラフ上で動的計画法を適用し、各パーティションの最小エネルギー割り当てを計算する。
• E₀₁の構造を互いに素なK₄,₄部分グラフとして利用し、グロテンジックの不等式と普遍定数Cを用いて双線形形式をバウンドする。
• E₀、E₁、E₀₁および線形項のバウンドを組み合わせることで、最小エネルギーの下界を確立する。
• ベクトル埋め込みと内積最大化を用いて、K₄,₄上での二次形式の最小値の下界を導出し、定数倍近似を可能にする。

実験結果

リサーチクエスチョン

• RQ1 非平面的チメラグラフ構造における古典的イジングスピンガラス問題に対して、PTASを構築可能か？
• RQ2 チメラグラフにおいて、地面状態エネルギーの自明な下界が真の最適解の定数倍以内にあるか？
• RQ3 Bansalら（2009年）の平面グリッドにおける技術を、チメラグラフのレイヤー構造と多成分構造に適応可能か？
• RQ4 最小エネルギーに対する定数倍下界を用いた場合、得られる最も良い近似要因は何か？

主な発見

• 任意のε>0に対して、チメラグラフ上でのイジングスピンガラス問題に対して(1−ε)-近似アルゴリズムが存在し、実行時間はO(n·2^(32/ε))である。
• 最小エネルギーの自明な下界は、真の最適解の(3C+4)/C ≤ 17.26倍以内であり、ここでC > ln(1+√2)/πである。
• ハミルトニアンH(S)の最小値は、−(C/(3C+4))(∑|c_uv| + ∑|d_u|)で上界から抑えられ、定数倍近似が確立された。
• O(1/ε)個のエッジを削除した後、定数幅ストリップへのグラフ分解によってPTASの存在が可能になった。
• 最小双線形形式のバウンドにはグロテンジックの不等式が用いられ、普遍定数Cが得られた。
• アルゴリズムはT=2/ε個のパーティション上でハミルトニアンを最小化し、最良の割り当てを選択することで近似を達成しており、測度集中の性質を満たしている。