物理学におけるメトロポリス・モンテカルロ法

メトロポリスアルゴリズムは、統計物理学シミュレーションの主力であり、提案された移動をそのエネルギーコストに基づいて受理または拒否することにより、正しいボルツマン確率で配置をサンプリングするマルコフ連鎖を構築します。

Definition

メトロポリスアルゴリズムは、マルコフ連鎖モンテカルロ法の一種であり、局所的な変化を提案し、エネルギー変化のボルツマン因子によって設定された確率でそれらを受理することにより、その極限分布が正準アンサンブルである配置のシーケンスを生成します。

Scope

このトピックでは、物理システムに適用されるメトロポリスアルゴリズムとメトロポリス・ヘイスティングスアルゴリズムについて扱います。具体的には、受理規則、詳細釣り合いとエルゴード性、平衡化と自己相関、および熱平均とその統計誤差の推定についてです。これは、より広範なモンテカルロ領域の基礎となるサンプリング手法です。

Core questions

• 受理確率は、提案された移動のエネルギー変化にどのように依存しますか？
• なぜ詳細釣り合いは正しい定常分布を保証するのですか？
• 平衡化時間と自己相関時間はどのように診断され、考慮されますか？
• 相関のあるサンプルからモンテカルロ平均の統計誤差はどのように推定されますか？

Key theories

Detailed balance and stationarity

Choosing acceptance probabilities that satisfy detailed balance with respect to the Boltzmann distribution ensures that distribution is stationary under the Markov chain, so long-run averages converge to thermal expectation values.

Metropolis-Hastings generalization

Hastings generalized the acceptance rule to asymmetric proposal distributions, broadening the algorithm beyond symmetric local moves while preserving the target stationary distribution.

Autocorrelation and error estimation

Successive Metropolis samples are correlated, so the effective number of independent samples is reduced by the autocorrelation time, which must be measured to assign honest error bars to thermal averages.

Clinical relevance

メトロポリスサンプリングは、格子スピンモデル、流体、高分子の熱力学的量を計算し、相転移を特定し、モンテカルロ分子シミュレーションや多くの量子モンテカルロスキームの中核エンジンとして機能します。

History

1953年にロスアラモスのMANIACコンピュータで二次元ハードディスク流体の状態方程式を計算するために導入されたこのアルゴリズムは、1970年にヘイスティングスによって一般化され、統計物理学、そして後にベイズ統計学において最も広く使用されるシミュレーション手法となりました。

Key figures

• Nicholas Metropolis
• Arianna Rosenbluth
• W. Keith Hastings

Seminal works

• metropolis1953
• hastings1970

Frequently asked questions

なぜエネルギーを低下させる移動は常に受理されるのですか？

エネルギーを低下させる移動はボルツマン重みを増加させるため、それを受理することは常に連鎖をより確率の高い状態へと向かわせます。上り坂の移動は、エネルギー増加によって設定された確率で時々受理されるだけであり、これにより連鎖は下り坂にのみ転がるのではなく、完全な熱分布を探索することができます。

なぜ実行開始時にサンプルを破棄する必要があるのですか？

連鎖は、まだ平衡分布を代表していない任意の配置から開始されます。初期の平衡化期間またはバーンイン期間は破棄され、測定された平均が開始時のバイアスではなく、真の熱アンサンブルを反映するようにします。