[論文レビュー] Dynamic critical behavior of cluster algorithms for 2D Ashkin-Teller and Potts models

Jesús Salas|arXiv (Cornell University)|Apr 19, 2000

Theoretical and Computational Physics参考文献 27被引用数 3

ひとこと要約

本稿は、2次元アーシン＝トレラー模型におけるスウェンセン–ワングおよび関連するクラスターアルゴリズムの動的臨界行動を調査し、自己相関時間のリー–ソカラウの下界が鋭くはないが、ほぼ鋭いものであることが判明した。自己相関時間の統合値と比熱の比は、系サイズの対数的または小さな指数（0.05–0.12）に従って増大する。エネルギーの指数的および統合自己相関時間は、同じ動的臨界指数に従うことが示され、近似的に最適な動的スケーリング行動を示している。

ABSTRACT

We study the dynamic critical behavior of two algorithms: the Swendsen-Wang algorithm for the two-dimensional Potts model with q=2,3,4 and a Swendsen-Wang-type algorithm for the two-dimensional symmetric Ashkin-Teller model on the self-dual curve. We find that the Li--Sokal bound on the autocorrelation time τ_{{ m int},{\cal E}} \geq const imes C_H is almost, but not quite sharp. The ratio τ_{{ m int},{\cal E}}/C_H appears to tend to infinity either as a logarithm or as a small power (0.05 \ltapprox p \ltapprox 0.12). We also show that the exponential autocorrelation time τ_{{ m exp},{\cal E}} is proportional to the integrated autocorrelation time τ_{{ m int},{\cal E}}.

研究の動機と目的

• スウェンセン–ワングアルゴリズムが一部の模型ではうまく機能するが他の模型ではうまくいかない理由、特に臨界点付近での挙動を理解すること。
• 自己相関時間のリー–ソカラウ下界の鋭さを検証すること。この下界は、動的臨界指数と静的臨界指数を関連付ける。
• 2次元クラスターアルゴリズムにおいて、エネルギーの指数的および統合自己相関時間の動的臨界指数が同一であるかどうかを特定すること。
• 非局所モンテカルロアルゴリズムのダイナミクスにおける遅いモードの性質を調査すること。

提案手法

• 2次元q状態ポッツ模型（q=2,3,4）および自己双対曲線上の対称アーシン＝トレラー模型に対するスウェンセン–ワングおよび関連クラスターアルゴリズムを用いた数値モンテカルロシミュレーション。
• 長時間シミュレーションからの高精度データを用いて、エネルギー観測量の統合的および指数的自己相関時間を推定する。
• 系サイズLに伴う自己相関時間のスケーリングを用いて、有限サイズスケーリング解析を適用し、動的臨界指数を抽出する。
• 自己相関関数ρ_AA(t;L) ≈ |t|^{-p_A} h_A[t/τ_exp,A; ξ(L)/L] のスケーリングアンサッツを用いて、エネルギーの自己相関関数が純粋な指数関数的挙動を示すかどうかを検証する。
• τ_int,ℰ / C_H とリー–ソカラウ下界の比較を通じて、その鋭さを評価し、対数的およびべき乗則的フィッティング形式を用いる。
• ρ_ℰℰ(t) と t/τ_int,ℰ のデータコラプスが異なる格子サイズで成立するかを確認することで、アンサッツ τ_int,ℰ ≈ τ_exp,ℰ の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 2次元系におけるクラスターアルゴリズムに対して、自己相関時間のリー–ソカラウ下界は鋭いか？
• RQ2 スウェンセン–ワングアルゴリズムにおいて、エネルギーの指数的および統合自己相関時間は、同じ動的臨界指数に従うか？
• RQ3 比 τ_int,ℰ / C_H は系サイズLとともにどのようにスケーリングするか？ 対数的増大か、小さなべき乗（L^p）か？
• RQ4 2次元ポッツおよびアーシン＝トレラー模型のダイナミクスにおける最も遅いモードの性質は何か？
• RQ5 なぜリー–ソカラウ下界は2次元ではほぼ鋭いが、高次元では明確に鋭くないのか？

主な発見

• リー–ソカラウ下界 τ_int,ℰ ≥ const × C_H は鋭くないが、比 τ_int,ℰ / C_H は系サイズLとともに対数的または小さな指数（0.05 ≤ p ≤ 0.12）に従って増大する。
• 2次元ポッツ模型（q=2,3,4）およびアーシン＝トレラー模型におけるエネルギーの動的臨界指数 z_int,ℰ は z_int,ℰ ≈ z_exp,ℰ と一致しており、指数的および統合自己相関時間が同じ指数に従うことを示唆している。
• 2次元イジング模型（q=2）では、比 τ_int,ℰ / C_H が対数的に増大するため、アルゴリズムの動的スケーリングがほぼ最適であることが示唆される。
• スケーリング τ_int,ℰ / C_H ∼ L^p における指数pは、イジング模型の約0.05から4状態ポッツ模型の0.12へと増加し、モデル間で連続的な傾向を示している。
• L ≥ 64 において、ρ_ℰℰ(t) と t/τ_int,ℰ のデータコラプスが極めて良好であるため、エネルギー自己相関関数が純粋な指数関数的挙動を示すというアンサッツの妥当性が支持される。
• 本研究は、d=2ではエネルギーと結合占有確率が最も遅いモードに属するが、d>2ではそうではないことから、リー–ソカラウ下界の鋭さの違いが説明できると示唆している。