スヴェンセン・ワンのアルゴリズム

スヴェンセン・ワンのアルゴリズム (Swendsen–Wang algorithm) は、イジング模型モンテカルロシミュレーション用のアルゴリズムの一つで、サイト全体を同一スピンのクラスターに分割し、クラスターそれぞれに新たなスピンを割り当てるアルゴリズムである。ウォルフのアルゴリズムも参照。

このアルゴリズムは、一回の状態変化において広域にわたる変化を実行するアルゴリズムとしては初期の物である。もともとのアルゴリズムはイジング模型ポッツ模型用であったが、後に一般化され、ウォルフのアルゴリズムによるXY模型流体粒子など、他の系にも用いられるようになった。この手法のカギは、Fortuin[訳語疑問点]と Kasteleyn[訳語疑問点] に帰せられるイジング模型やポッツ模型を結合のパーコレーションでモデル化する方法である。結合されたサイトはクラスターを形成する。等しいスピンを持つサイト同士は次に示す確率で結合される。

P = 1 − exp(−2J/(kBT))

ここで、 J強磁性イジング模型におけるカップリング定数、 T温度kBボルツマン定数である。クラスターは同じ確率で「反転」される。この手法は二次相転移点付近において臨界減速を克服でき、最も効率がよい手法である。

Barbu & Zhu 2005で一般化がなされ、任意のサンプリング確率分布の場合にも提案されたモンテカルロ遷移の採択確率を計算することができるようになり、メトロポリス・ヘイスティングス法を実行できるようになった。

出典

外部リンク