[論文レビュー] Optimal Mixing of Glauber Dynamics: Entropy Factorization via High-Dimensional Expansion

Zongchen Chen, Kuikui Liu+1 more|arXiv (Cornell University)|Nov 4, 2020

Markov Chains and Monte Carlo Methods参考文献 57被引用数 5

ひとこと要約

本稿は、有界次数のグラフ上の反鉄相反2スピン系および関連モデルにおけるGlauberダイナミクスの最適な $O(n\log n)$ 混合時間の確立を、高次元拡張に基づく新規なエントロピー因子化フレームワークを用いて行った。ブロック因子化によるエントロピーの近似的なテンソル化を証明し、局所的からグローバルなスペクトル拡張を活用することで、木の一意性領域におけるタイトな混合時間の上限を達成した。これは、臨界閾値未満のファージャシティにおけるハードコアモデルや、$q > \alpha\Delta$（$α \approx 1.763$）の三角形を含まないグラフの$q$-彩色に対しても成立する。

ABSTRACT

We prove an optimal mixing time bound on the single-site update Markov chain known as the Glauber dynamics or Gibbs sampling in a variety of settings. Our work presents an improved version of the spectral independence approach of Anari et al. (2020) and shows $O(n\log{n})$ mixing time on any $n$-vertex graph of bounded degree when the maximum eigenvalue of an associated influence matrix is bounded. As an application of our results, for the hard-core model on independent sets weighted by a fugacity $λ$, we establish $O(n\log{n})$ mixing time for the Glauber dynamics on any $n$-vertex graph of constant maximum degree $Δ$ when $λ<λ_c(Δ)$ where $λ_c(Δ)$ is the critical point for the uniqueness/non-uniqueness phase transition on the $Δ$-regular tree. More generally, for any antiferromagnetic 2-spin system we prove $O(n\log{n})$ mixing time of the Glauber dynamics on any bounded degree graph in the corresponding tree uniqueness region. Our results apply more broadly; for example, we also obtain $O(n\log{n})$ mixing for $q$-colorings of triangle-free graphs of maximum degree $Δ$ when the number of colors satisfies $q > αΔ$ where $α\approx 1.763$, and $O(m\log{n})$ mixing for generating random matchings of any graph with bounded degree and $m$ edges.

研究の動機と目的

• 有界次数のグラフ上のスピン系における木の一意性領域で、Glauberダイナミクスが $O(n\\log n)$ 時間で混合するという長年の予想を解決すること。
• ブロック因子化によるエントロピーの近似的なテンソル化を確立することで、最適混合時間の一般枠組みを構築すること。
• 従来の結果を超えて、ハードコアモデル、$q$-彩色、モノマー・ダイマー・モデルを含む広範なモデルへ、スペクトル独立性の技術を拡張すること。
• スペクトル独立性と高次元拡張子の先行研究を統合・一般化し、特にマルコフ連鎖の混合時間の文脈で統一すること。
• スペクトル独立性と周辺確率の有界性に依存する混合時間のタイトな定量的境界を提供し、従来の次数やパrameterに依存する指数的依存を改善すること。

提案手法

• 個々の頂点ではなく、固定された線形サイズのブロックへの因子化を通じて、エントロピーの近似的なテンソル化を証明する、エントロピー因子化の新手法を導入。
• AlevとLau（2020）の局所的からグローバルなスキームを、局所スペクトル拡張を満たす重み付き単体複体へ適応し、ブロックエントロピー因子化を確立。
• 高次元拡張の性質を用いて影響行列を制御し、最大固有値を評価することで、スペクトル独立性の議論を可能にする。
• スペクトル独立性と周辺確率の有界性を組み合わせ、定理1.11を用いて $O(n\log n)$ 混合時間の上限を導出。この定理は、これらの性質が $O(n\log n)$ 混合にどのように関連するかを示している。
• ハードコアモデル、三角形を含まないグラフの $q$-彩色、モノマー・ダイマー・モデルの各具体例に、フレームワークを適用。それぞれのケースでスペクトル独立性と周辺確率の有界性を検証。
• 既知の $q$-彩色およびモノマー・ダイマー・モデルのスペクトル独立性の結果を、新規なブロックベースのエントロピー因子化フレームワークに適応。

実験結果

リサーチクエスチョン

• RQ1 有界次数のグラフ上の反鉄相反2スピン系の木の一意性領域において、Glauberダイナミクスの $O(n\log n)$ 混合時間の上限が確立可能か。
• RQ2 修正された対数ソボレフ不等式の証明におけるエントロピー因子化ステップを、ブロックベースの因子化に置き換えることで、よりタイトな混合時間の上限が得られるか。
• RQ3 局所的からグローバルなスペクトル拡張フレームワークは、従来の設定を超えて一般の重み付き単体複体へ拡張可能か。
• RQ4 最大次数 $\Delta$ およびスペクトル独立性パrameter $\eta$ に依存する混合時間の依存関係を、指数的から多項式的へ改善可能か。
• RQ5 正則木上のモノマー・ダイマー・モデルの最適スペクトル独立性定数は何か。また、その値は影響行列ノルムから導かれる境界と一致するか。

主な発見

• 最大次数 $\Delta$ の任意の $n$ 頂点グラフにおけるハードコアモデルのGlauberダイナミクスは、$\lambda < \lambda_c(\Delta)$（$\Delta$-正則木における臨界ファージャシティ）のとき、$O(n\log n)$ 時間で混合する。
• 木の一意性領域に位置する反鉄相反2スピン系に対して、任意の有界次数のグラフ上でのGlauberダイナミクスは $O(n\log n)$ 時間で混合する。
• 三角形を含まないグラフの $q$-彩色で $q > \alpha\Delta$（$α \approx 1.763$）のとき、Glauberダイナミクスは $O(n\log n)$ 時間で混合する。
• $q \geq 3\Delta$ の場合、ドブリージンの一意性条件により、混合時間は $2n\log(n/\varepsilon)$ で有界である。
• $m$ 棄の辺を持つ有界次数のグラフ上のモノマー・ダイマー・モデルでは、Glauberダイナミクスは $O(m\log n)$ 時間で混合する。
• 無限大の $\Delta$-正則木上のモノマー・ダイマー・モデルのスペクトル独立性は、$\frac{2x}{1-x}$ に正確に等しく、$x = \frac{1}{\Delta-1}\left(1 - \frac{2}{\sqrt{1+4\lambda(\Delta-1)}+1}\right)$ である。$\Delta \geq 3$ のとき、これは $O(1/\Delta)$ の依存性を示している。