[論文レビュー] The computational hardness of counting in two-spin models on d-regular graphs

Allan Sly, Nike Sun|arXiv (Cornell University)|Mar 12, 2012

Markov Chains and Monte Carlo Methods参考文献 17被引用数 60

ひとこと要約

本稿では、d-正則グラフ上で、非一意性がd-正則木上で成立する場合、2スピン系（例：ハードコアモデルや反強磁性イジング模型）の分配関数の近似がNP困難であることを確立している。著者らは、局所的に木構造を持つ拡張グラフを用い、自由エネルギー密度がベーテ予測に収束することを示し、従来の2次モーメント法の手法を回避することで、最大カットへの新たな確率的還元を介してこの結果を証明した。

ABSTRACT

The class of two-spin systems contains several important models, including random independent sets and the Ising model of statistical physics. We show that for both the hard-core (independent set) model and the anti-ferromagnetic Ising model with arbitrary external field, it is NP-hard to approximate the partition function or approximately sample from the model on d-regular graphs when the model has non-uniqueness on the d-regular tree. Together with results of Jerrum--Sinclair, Weitz, and Sinclair--Srivastava--Thurley giving FPRAS's for all other two-spin systems except at the uniqueness threshold, this gives an almost complete classification of the computational complexity of two-spin systems on bounded-degree graphs. Our proof establishes that the normalized log-partition function of any two-spin system on bipartite locally tree-like graphs converges to a limiting "free energy density" which coincides with the (non-rigorous) Bethe prediction of statistical physics. We use this result to characterize the local structure of two-spin systems on locally tree-like bipartite expander graphs, which then become the basic gadgets in a randomized reduction to approximate MAX-CUT. Our approach is novel in that it makes no use of the second moment method employed in previous works on these questions.

研究の動機と目的

• 2スピン系の分配関数の近似の計算複雑性を、有界次数のグラフ上で分類すること。
• d-正則木上の一意性閾値が、2スピン系の計算的遷移点を正確に示すかという未解決の問題を解明すること。
• 従来の結果で用いられた2次モーメント法に依存しない、新たな証明技術を提供すること。
• 外部場が任意の値をとる反強磁性イジング模型に対しても、硬度結果を拡張すること。
• 同質の2スピン系がd-正則グラフ上で示す計算複雑性の分類を完了させること（一意性閾値を除く）。

提案手法

• 硬度還元のための道具として、双方向で局所的に木構造を持つ拡張グラフを構築する。
• 2スピン系の正規化された対数分配関数が、統計物理学におけるベーテ予測に一致する極限自由エネルギー密度に収束することを証明する。
• 自由エネルギーの収束を用いて、局所的に木構造を持つ双方向拡張グラフ上での2スピン系の局所的構造を特徴付ける。
• 2スピン系の分配関数の近似への最大カットの近似への確率的還元を設計する。
• 3-正則グラフ $ H $ の各頂点を基本グラフ $ G $ のコピーに置き換え、対応する境界集合間に $ 2k $ 条の辺を接続することで、グラフ構成 $ H^G $ を定義する。
• 2スピン系から導かれるパrameter $ \Gamma $ と $ \Theta $ を用いて、$ Z_{H^G} $ と $ Z_{\widehat{H}^G} $ の比に上限を設定し、分配関数と $ H $ の最大カットを関連付ける。

実験結果

リサーチクエスチョン

• RQ1 d-正則木上の一意性閾値は、d-正則グラフ上での2スピン系の近似の計算的容易性と困難性の正確な境界を示すか？
• RQ2 2次モーメント法の計算に依存せずに、分配関数の近似の計算的硬度を確立できるか？
• RQ3 局所的に木構造を持つグラフ上での2スピン系の漸近的挙動を、ベーテ予測による自由エネルギー密度が正確に記述できるか？
• RQ4 拡張グラフ上での2スピン系を用いて、最大カットへの確率的還元を構築し、近似の硬度を証明できるか？
• RQ5 外部場が任意の値をとる反強磁性イジング模型の計算的遷移点も、d-正則木上の一意性閾値に一致するか？

主な発見

• d \geq 3 かつ \lambda > \lambda_c(d) = \frac{(d-1)^{d-1}}{(d-2)^d} のとき、d-正則グラフ上のハードコアモデルに対して、FPRASが存在するなら NP = RP でなければならない。
• d \geq 3、B \in \mathbb{R}、かつ \beta < \beta_{c,\text{af}}(B,d) のとき、外部場 B を持つ反強磁性イジング模型に対して、d-正則グラフ上にFPRASが存在するなら NP = RP でなければならない。
• 双方向で局所的に木構造を持つグラフ上での任意の2スピン系について、正規化された対数分配関数は、ベーテ予測に一致する極限自由エネルギー密度に収束する。
• 著者らは、$ H^G $ 上での2スピンモデルの分配関数が、$ H $ の最大カットと密接に関係しており、比が $ (1 \pm \epsilon)^m $ の範囲に収束することを確立した。これにより、確率的還元が可能となった。
• 証明は2次モーメント法を回避しており、従来の細密な2次モーメント計算に依存する手法に比べ、より概念的かつ一般的な枠組みを提供している。
• 結果として、同質の2スピン系がd-正則グラフ上で示す計算複雑性の分類が完了し、非一意性領域では硬度が示され、一意性領域ではFPTASが既知である。