[論文レビュー] Convergence to equilibrium for finite Markov processes, with application to the Random Energy Model

Pierre Mathieu, Pierre Picco|arXiv (Cornell University)|Jul 10, 2003

Stochastic processes and statistical mechanics参考文献 11被引用数 7

ひとこと要約

この論文は、初期分布依存性を組み込み、スペクトルギャップとPoincaré不等式の手法を用いて、有限で可逆なマルコフ過程の平衡への収束に対する精密な上限を導出する。Derridaのランダムエネルギー模型（REM）のメトロポリスダイナミクスに適用することで、対数混合時間スケーリングのよりきつい上界 $eta^2$ を確立し、$eta \leq \beta_c$ の場合に先行研究より2倍の改善を達成する。また、確率的環境実現の影響を受けることなく、収束時間の鋭い推定値を提供する。

ABSTRACT

We estimate the distance in total variation between the law of a finite state Markov process at time t, starting from a given initial measure, and its unique invariant measure. We derive upper bounds for the time to reach the equilibrium. As an example of application we consider a special case of finite state Markov process in random environment: the Metropolis dynamics of the Random Energy Model. We also study the process of the environment as seen from the process.

研究の動機と目的

• 初期分布に明示的に依存する収束時間推定を、初期分布の法則に一様なものではなく、得ること。
• ランダムエネルギー模型（REM）のメトロポリスダイナミクスの混合時間に対する既存の上界を改善すること、特に低温領域（$\beta \leq \beta_c$）において。
• 環境実現 $\omega$ に依存しない収束時間推定 $T_N(\epsilon, c, \eta)$ を確立し、ほとんど確実な漸近的境界を可能にすること。
• スペクトルギャップ法を、状態空間上の経路の最適化を丁寧に行うことで組み合わせることで、分配関数に基づく推定よりも鋭い境界が得られることを示すこと。

提案手法

• 一般化されたPoincaré不等式とスペクトルギャップ推定を用いて、時間 $t$ における過程の分布と不変測度との全 Variation 距離を評価する。
• 環境実現の高確率で成り立つ決定的収束時間 $T_N(\epsilon, c, \eta)$ を導入し、ほとんど確実な漸近的境界を保証する。
• メトロポリス遷移行列によって誘導される重み付きグラフ上で経路に基づく解析を適用し、混合時間推定を最適化するための特定の経路を構築する。
• 指数マルコフ不等式と大偏差推定を用いて、REMにおけるエネルギー区間の訪問回数を制御し、制約付き分配関数の制御を可能にする。
• エネルギー空間を区間 $\Delta_k$ に分割し、集中不等式を用いて各区間の状態数の大きな偏差の確率を評価する。
• 制約付き分配関数 $Z_\alpha(\beta, \geq -d)$ の推定とスペクトルギャップを組み合わせ、最終的な混合時間の上限を導出する。

実験結果

リサーチクエスチョン

• RQ1 有限マルコフ連鎖の収束時間境界は、初期分布への依存性を組み込むことで改善可能か？
• RQ2 ランダムエネルギー模型のメトロポリスダイナミクスの混合時間の鋭い漸近的スケーリングは何か？
• RQ3 経路に基づく手法は、スピンガラス模型の混合時間推定において、分配関数に基づく手法と比べてどのように異なるか？
• RQ4 $β \leq \beta_c$ の場合に $φ\text{\textbackslash{}limsup}_{N\to\infty} \frac{1}{N}\log T_N(\epsilon, c, \eta) \leq \beta^2$ であるという境界は最適か？

主な発見

• この論文は、REMにおけるメトロポリスダイナミクスの収束時間に対して、新たな上界 $φ\text{\textbackslash{}limsup}_{N\to\infty} \frac{1}{N}\log T_N(\epsilon, c, \eta) \leq \beta^2$ を確立し、$β \leq \beta_c$ のすべての値に対して有効である。
• この境界は、$β \leq \beta_c$ の場合に、先行研究の $2\beta^2$ や $2\beta\beta_c$ の結果よりも指数部で2倍の改善を達成する。
• 境界 $β^2$ が鋭いと想定されており、$β \leq \beta_c$ の場合に $φ\text{\textbackslash{}lim}_{N\to\infty} \frac{1}{N}\log T_N(\epsilon, c, \eta) = \beta^2$ であると予想されている。
• 状態空間における特定の経路の構築を用いる手法は、分配関数に基づく推定手法よりも優れた推定をもたらすことが示されており、指数部の改善が得られている。
• $β > \beta_c$ の場合、境界 $β\beta_c$ は依然有効であるが、$β \leq \beta_c$ の場合に新しい $β^2$ の境界がよりきつい。
• 収束時間推定 $T_N(\epsilon, c, \eta)$ は決定的であり、環境実現の高確率で成り立つため、$T^\omega(\epsilon)$ とのほとんど確実な漸近的比較が可能になる。