[論文レビュー] Spontaneous magnetisation in the plane

Geoff K. Nicholls|arXiv (Cornell University)|Jul 19, 2000

Theoretical and Computational Physics参考文献 11被引用数 5

ひとこと要約

本稿では、平面に非交差色付き多角形を生成する可解な確率過程であるアラク過程から導かれる温度修正型連続的イジング模型を提案する。アラク過程をグランドカノニカル集合におけるギブス分布として解釈することで、非ゼロの臨界温度で相転移が発生することを証明した。モンテカルロシミュレーションにより、2次元イジング普遍性クラスと整合する臨界指数が得られた。

ABSTRACT

The Arak process is a solvable stochastic process which generates coloured patterns in the plane. Patterns are made up of a variable number of random non-intersecting polygons. We show that the distribution of Arak process states is the Gibbs distribution of its states in thermodynamic equilibrium in the grand canonical ensemble. The sequence of Gibbs distributions form a new model parameterised by temperature. We prove that there is a phase transition in this model, for some non-zero temperature. We illustrate this conclusion with simulation results. We measure the critical exponents of this off-lattice model and find they are consistent with those of the Ising model in two dimensions.

研究の動機と目的

• アラク過程の分布をグランドカノニカル集合におけるギブス測度として解釈することで、アラク過程と統計力学との間の接続を確立すること。
• 分配関数が解析的に扱えなくなるが、温度修正アラク模型において相転移が発生するかを調査すること。
• マルコフ連鎖モンテカルロ法を用いてモデルをシミュレートし、臨界指数を測定することで、2次元イジング模型との普遍性を検証すること。
• グラウンド状態の degeneracy と系の2次元性が、非格子的・多角形構造であるにもかかわらず、イジング的臨界行動を引き起こすかを検討すること。
• 臨界温度を推定し、磁化およびバインダーパラメータのデータコラプスと回帰分析を用いて臨界行動を検証すること。

提案手法

• アラク過程を、黒領域と白領域の間の全界面長に比例するエネルギーを持つギブス分布として再解釈する。
• モデルは温度でパrameter化され、元のアラク過程は T=1 に対応する。低温では閉形式の分配関数が失われるが、モデルの定義は保たれる。
• マルコフ連鎖モンテカルロ（MCMC）シミュレーションを用いて実現を生成し、頂点の生成・消滅、頂点の移動、辺の入れ替え、領域の再色分けなどの更新手順を用いる。
• サイズ d が異なる有限系に対して磁化およびバインダーパラメータを計算し、イジング模型の臨界指数を用いてデータコラプスを適用し、普遍性を検証する。
• 再スケーリングされた磁化およびバインダーパラメータのデータに対する回帰分析により臨界指数を推定し、適合度の評価に χ² 統計量を用いる。
• 角の分割とヤコビアンの制御を用いて多角形積分の理論的境界を導出し、分配関数の収束を確立する。

実験結果

リサーチクエスチョン

• RQ1 温度修正アラク模型は非ゼロの臨界温度で相転移を示すか？
• RQ2 温度修正アラク模型の臨界指数は2次元イジング模型のそれと整合的であり、普遍性を示すか？
• RQ3 本モデルの推定臨界温度 Tc は何か？ また、T=1 の元のアラク過程と比較するとどうなるか？
• RQ4 イジング臨界指数を用いて再スケーリングした場合、シミュレートされた磁化およびバインダーパラメータのデータはどの程度コラプスするか？
• RQ5 アラク過程はグランドカノニカル集合におけるギブス測度として解釈可能か？ また、この解釈は相共存の存在を支持するか？

主な発見

• 閉形式の分配関数が失われるにもかかわらず、厳密に証明されたように、温度修正アラク模型では非ゼロの臨界温度で相転移が発生する。
• 推定された臨界温度は Tc ≈ 0.6665(5) であり、シミュレーションでは T=1（元のアラク過程）では単一相領域であることが示された。
• イジング模型の指数を用いて再スケーリングした磁化およびバインダーパラメータの臨界指数は、統計的に妥当な適合度を示した（バインダーパラメータの適合では χ² = 38.5）。
• イジング指数を用いて再スケーリングした磁化データは適合度が悪い（χ² 統計量により適合が不十分であると示唆された）ため、有限サイズ効果や臨界領域におけるずれの可能性がある。
• モデルのグラウンド状態は二重に degenerate であり、系の2次元性がイジング普遍性の仮説を支持しており、観測された臨界指数と整合的である。
• T=1 におけるアラク過程は、既知の正規化定数を持つグランドカノニカルギブス測度に対応するが、T<1 の場合でもモデルは定義され、相共存を示す。