[論文レビュー] Inapproximability After Uniqueness Phase Transition in Two-Spin Systems

Jin‐Yi Cai, Xi Chen+2 more|arXiv (Cornell University)|May 14, 2012

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

ひとこと要約

本稿は、一意性位相遷移の閾値を超えて、2スピン系の分配関数の計算的近似不能性を確立する。E2LIN2問題からの還元とグラフ集合の確率的解析を用いて、d-正則グラフにおいてd = Ω(Δ(β,γ))のとき、相対誤差ϵ = 10⁻⁴で分配関数を近似できるランダム化多項式時間アルゴリズムが存在しない（NP ≠ RPでない限り）ことを証明し、統計力学の位相遷移と計算複雑性の間の長年の予想を確認する。

ABSTRACT

A two-state spin system is specified by a 2 x 2 matrix A = {A_{0,0} A_{0,1}, A_{1,0} A_{1,1}} = {β1, 1 γ} where β, γ\ge 0. Given an input graph G=(V,E), the partition function Z_A(G) of a system is defined as Z_A(G) = \sum_{σ: V -> {0,1}} \prod_{(u,v) \in E} A_{σ(u), σ(v)} We prove inapproximability results for the partition function in the region specified by the non-uniqueness condition from phase transition for the Gibbs measure. More specifically, assuming NP e RP, for any fixed β, γin the unit square, there is no randomized polynomial-time algorithm that approximates Z_A(G) for d-regular graphs G with relative error ε= 10^{-4}, if d = Ω(Δ(β,γ)), where Δ(β,γ) > 1/(1-βγ) is the uniqueness threshold. Up to a constant factor, this hardness result confirms the conjecture that the uniqueness phase transition coincides with the transition from computational tractability to intractability for Z_A(G). We also show a matching inapproximability result for a region of parameters β, γoutside the unit square, and all our results generalize to partition functions with an external field.

研究の動機と目的

• 2スピン系における一意性位相遷移が、分配関数近似の計算的非効率性の閾値と一致するという予想を解決すること。
• 特に有界次数のグラフにおいて、反強磁性2スピン系の非一意的領域における強い近似不能性結果を確立すること。
• 単位正方形（β, γ ∈ [0,1]）を超えて、この領域外のパラメータへの近似不能性結果を拡張すること。
• 外部場を含む系へと硬度結果を一般化し、重み付きグラフ同型写像やイジング型モデルへの適用範囲を広げること。
• 硬度が正確に一意性閾値で始まることを証明することで、統計力学的位相遷移と計算複雑性の間にきわめて明確な接続を確立すること。

提案手法

• E2LIN2問題をランダム化されたグラフ構築により2スピン系に還元し、E2LIN2の難易度をd-正則グラフの分配関数に埋め込む。
• 制御された次数と局所的構造（Ui,k, Vi,k）を持つランダムグラフ集合H(N, ∆)を構築し、ギブス測度をシミュレートし期待重みを分析する。
• モーメントの方法と指数的モーメントの上限を用いて、分配関数ZA(G)の期待値を推定し、エントロピー関数H(x) = −x ln x − (1−x) ln(1−x)を用いて指数部の主要項に注目する。
• チェルノフ不等式と集中不等式を適用し、局所的近傍における0割り当て頂点数がその平均のまわりにきわめて集中していることを示す。
• マーカフの不等式と和集合の上限を用いて、まれな構成からの分配関数寄与の尾部挙動を制御する。
• 期待分配関数の指数部に対するタイトな上界を導出し、それが1.22·dim未満であることを示し、典型的な構成とまれな構成との間で指数的分離が生じることを示す。

実験結果

リサーチクエスチョン

• RQ1 2スピン系における一意性位相遷移が、分配関数近似の計算的容易性と困難性の正確な境界を示しているか？
• RQ2 d-正則グラフにおいて、dが一意性閾値Δ(β,γ)を超えるときでさえ、反強磁性系においても近似不能性が証明可能か？
• RQ3 硬度結果は単位正方形（β, γ ∈ [0,1]）を超えて、特に強磁性または混合領域においても頑健か？
• RQ4 外部場を含む系へと近似不能性結果を拡張可能か？また、これにより計算の閾値はどのように変化するか？
• RQ5 グラフの構造（例：正則性、局所的木的性質）は、分配関数近似の難易度にどの程度影響を与えるか？

主な発見

• β, γ ∈ [0,1] かつ (β, γ) ≠ (0,0), (1,1) である任意の組に対して、d = Ω(1/(1−βγ)) であるd-正則グラフにおいて、相対誤差ϵ = 10⁻⁴でZA(G)を近似できるランダム化多項式時間アルゴリズムは存在しない（NP = RPでない限り）。
• 硬度閾値d = Ω(Δ(β,γ))は、一意性閾値Δ(β,γ)と定数倍の違いを除いて一致しており、位相遷移と複雑性遷移の一致という予想を確認する。
• 近似不能性結果はd-正則グラフに制限されても成立し、硬度が次数が一意性閾値を超えたことによるものであり、無限大の次数によるものではないことを示している。
• 本稿は、単位正方形外のパラメータに対しても一致する近似不能性結果を確立し、位相遷移と複雑性の一致が反強磁性領域を超えて成立することを示している。
• 典型的な構成では分配関数の期待値がexp(1.22·dim)で抑えられ、まれな構成からの寄与は指数的に小さくなることが示され、近似フレームワークにおける強い分離が可能になる。
• 一意性条件を満たす反強磁性2スピン系に対するアルゴリズム的FPTASが最適であることが分析で確認され、一意性条件が満たされない（すなわちd ≥ Δ(β,γ)）と直ちに近似不能性が生じることを示している。