[論文レビュー] Rapid Mixing of Glauber Dynamics up to Uniqueness via Contraction

Zongchen Chen, Kuikui Liu+1 more|arXiv (Cornell University)|Apr 20, 2020

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

ひとこと要約

本論文は、相関崩壊と多項式補間に用いられる木再帰の収縮と、一般の2スピン系におけるGlauberダイナミクスの高速混合を結びつける統一的な枠組みを確立する。収束性を保証するポテンシャル関数が存在する場合、Glauberダイナミクスは一意性閾値まで多項式時間で混合することを証明し、反強磁性および強磁性モデルの両方におけるアルゴリズム的手法を統一する。

ABSTRACT

For general antiferromagnetic 2-spin systems, including the hardcore model and the antiferromagnetic Ising model, there is an $\mathsf{FPTAS}$ for the partition function on graphs of maximum degree $Δ$ when the infinite regular tree lies in the uniqueness region by Li et al. (2013). Moreover, in the tree non-uniqueness region, Sly (2010) showed that there is no $\mathsf{FPRAS}$ to estimate the partition function unless $\mathsf{NP}=\mathsf{RP}$. The algorithmic results follow from the correlation decay approach due to Weitz (2006) or the polynomial interpolation approach developed by Barvinok (2016). However the running time is only polynomial for constant $Δ$. For the hardcore model, recent work of Anari et al. (2020) establishes rapid mixing of the simple single-site Markov chain known as the Glauber dynamics in the tree uniqueness region. Our work simplifies their analysis of the Glauber dynamics by considering the total pairwise influence of a fixed vertex $v$ on other vertices, as opposed to the total influence on $v$, thereby extending their work to all 2-spin models and improving the mixing time. More importantly our proof ties together the three disparate algorithmic approaches: we show that contraction of the tree recursions with a suitable potential function, which is the primary technique for establishing efficiency of Weitz's correlation decay approach and Barvinok's polynomial interpolation approach, also establishes rapid mixing of the Glauber dynamics. We emphasize that this connection holds for all 2-spin models (both antiferromagnetic and ferromagnetic), and existing proofs for correlation decay or polynomial interpolation immediately imply rapid mixing of Glauber dynamics. Our proof utilizes that the graph partition function divides that of Weitz's self-avoiding walk trees, leading to new tools for analyzing influence of vertices.

研究の動機と目的

• 2スピン系における3つの主要なアルゴリズム的手法—相関崩壊、多項式補間、およびマコフ連鎖モンテカルロ—の分析を統一すること。
• すべての2スピンモデル（反強磁性および強磁性を含む）において、木の唯一性領域におけるGlauberダイナミクスの高速混合を確立すること。
• 適切なポテンシャル関数による木再帰の収縮が、Glauberダイナミクスの多項式時間混合を示すこと。
• 頂点ごとの全ペアワイズ影響を用いた影響解析の簡素化により、ハードカーリングモデルに関する先行結果を一般の2スピン系に拡張すること。
• 相関崩壊または多項式補間の既存の証明が、同じポテンシャル関数枠組みを用いて即座にGlauberダイナミクスの高速混合を示すこと。

提案手法

• 2スピン系における木再帰の収縮を分析する一般の$(\alpha, c)$-ポテンシャル関数枠組みを導入する。
• 固定頂点$v$が他の頂点に与える全ペアワイズ影響を定義し、従来の$v$への影響分析を置き換えることで、解析を簡素化する。
• グラフの分割関数がワイツの自己回避ウォーク木の分割関数を割り切ることを用いて、影響解析のための新規な道具を導出する。
• ポテンシャル関数$\Psi$における収縮条件を適用し、スピンギャップ解析を用いてGlauberダイナミクスの混合時間をバウンドする。
• $\Psi$が収縮および有界性条件を満たす場合、Glauberダイナミクスは$O(n^{2+2c/\alpha})$時間で混合することを証明する。
• ${\sqrt{\beta\gamma} \leq \frac{\Delta-2}{\Delta}}$および$\gamma > 1$のとき、暗黙の$\Psi$関数は一般の$(\alpha, c)$-ポテンシャルであり、$\alpha \geq \delta/2$および$c \leq 18$（$\beta = 0$の場合は$c \leq 4$）を満たす。

実験結果

リサーチクエスチョン

• RQ1 相関崩壊や多項式補間に用いられる木再帰の収縮を、Glauberダイナミクスの高速混合に直接結びつけることは可能か？
• RQ2 頂点が他の頂点に与える全ペアワイズ影響は、Glauberダイナミクスの混合時間解析において、より単純かつ一般性の高い枠組みを提供するか？
• RQ3 相関崩壊や補間の効率性を保証する同じポテンシャル関数が、Glauberダイナミクスの高速混合を保証するか？
• RQ4 一般の2スピン系におけるGlauberダイナミクスの混合時間は、一意性ギャップ$\delta$にどのように依存するか？
• RQ5 結果は反強磁性および強磁性の両方の2スピンモデル、特にハードカーリングモデルやイジング模型にどのように拡張されるか？

主な発見

• 一般の$(\alpha, c)$-ポテンシャル関数が存在する場合、一般の2スピン系におけるGlauberダイナミクスは$O(n^{2+2c/\alpha})$時間で混合する。
• ${\sqrt{\beta\gamma} \leq \frac{\Delta-2}{\Delta}}$および$\gamma > 1$のとき、暗黙のポテンシャル$\Psi$は一般の$(\alpha, c)$-ポテンシャルであり、$\alpha \geq \delta/2$および$c \leq 18$を満たし、先行の境界を改善する。
• $\beta = 0$のとき、同じポテンシャルは$c \leq 4$をもたらし、ハードカーリングモデルにおけるよりタイトな混合時間境界を示す。
• 恒等関数$\Psi(y) = y$は、特定の条件下で強磁性2スピン系において$(\alpha, c)$-ポテンシャルであり、$\alpha = \Theta(\delta)$および$c \leq O(1)$を満たし、$O(n^{2+c/\delta})$の混合時間につながる。
• 強磁性系で${\sqrt{\beta\gamma} \geq \frac{\Delta}{\Delta-2}}$かつ$\lambda$が有界であるとき、混合時間は$O(n^C)$（$\beta, \gamma, \lambda, \Delta$に依存する定数$C$）である。
• この枠組みは、既存の相関崩壊および補間の結果を統合する：ポテンシャル関数を用いた相関崩壊または補間の有効性の証明は、即座にGlauberダイナミクスの高速混合を示す。