[論文レビュー] Numerical study of ground state energy fluctuations in spin glasses

Claudio Giberti, Cecilia Vernia|ArXiv.org|Jul 8, 2008

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

ひとこと要約

本稿では、スピンガラスにおける基底状態エネルギーの揺らぎを研究するための確率的アニーリングアルゴリズムを提示している。主に、シャーリングトン＝キルパトリック（SK）モデルとエドワーズ＝アンドリューズ（EA）モデルに焦点を当て、エネルギー揺らぎのべき乗則的スケーリングを確認した。その結果、d=2のEAモデルではρ ≈ 0.501、d=2のSKモデルではω ≈ 1.25の指数が得られ、熱力学的極限において基底状態エネルギー分布はガウス型に近づくことが示された。これは有限サイズのSK系で観察されるガムベル型の挙動とは対照的である。

ABSTRACT

Using a stochastic algorithm introduced in a previous paper, we study the finite size volume corrections and the fluctuations of the ground state energy in the Sherrington-Kirkpatrick and the Edwards-Anderson models at zero temperature. The algorithm is based on a suitable annealing procedure coupled with a balanced greedy-reluctant strategy that drives the systems towards the deepest minimum of the energy function.

研究の動機と目的

• クエンチド秩序を有するスピンガラスモデルにおける基底状態エネルギー揺らぎの有限サイズスケーリングを調査すること。
• 複雑なスピンガラス系における基底状態の計算に適した確率的アニーリングアルゴリズムの妥当性を検証すること。
• エネルギー密度の標準偏差と平均エネルギー密度のスケーリング指数ρおよびωを特定すること。
• 熱力学的極限における基底状態エネルギー確率分布p_N(e_N)の極限的形状を分析すること。

提案手法

• グリーディと不満な戦略をバランスさせた確率的アニーリングアルゴリズムを用い、スピンガラスハミルトニアンにおけるグローバルエネルギー最小値を探索した。
• アルゴリズムは、平均場型のシャーリングトン＝キルパトリック（SK）モデルと短距離相互作用型のエドワーズ＝アンドリューズ（EA）モデルの両方へ適用された。
• クエンチド結合Jの多数の独立した実現の平均を取ることで、エネルギー密度の揺らぎが計算された。
• 有限サイズスケーリングは、べき乗則フィットε_N = e_∞ + bN^(-ω) および σ_N = aN^(-ρ) を用い、Nは系のサイズを表す。
• 分布p_N(e_N)は、x_N = (e_N - ε_N)/σ_N としてスケーリングされ、大N極限における極限形の分析が行われた。
• ガウス型への収束を評価するために、歪度と尖度の統計的モーメントが計算された。

実験結果

リサーチクエスチョン

• RQ1 エドワーズ＝アンドリューズモデルにおいて、系のサイズNに伴い標準偏差σ_Nと平均エネルギー密度ε_Nはどのようにスケーリングするか？
• RQ2 EAモデルにおいて、熱力学的極限において基底状態エネルギー密度の確率分布p_N(e_N)はガウス分布に収束するか？
• RQ3 2次元および3次元のEAモデル（d=2, d=3）における平均エネルギーおよびその揺らぎの有限サイズ補正指数ωおよびρは何か？
• RQ4 p_N(e_N)の形状は系のサイズに伴いどのように変化するか？また、極値理論の予測通り、Gumbel型からガウス型への遷移が見られるか？
• RQ5 確率的アニーリングアルゴリズムは、異なるスピンガラスモデルにおいて基底状態エネルギーを正確に計算するのに有効か？

主な発見

• 2次元のエドワーズ＝アンドリューズモデルにおいて、基底状態エネルギーの標準偏差はσ_N ∝ N^(-0.501) とスケーリングし、最良のフィット指数はρ = 0.501 ± 6×10⁻³であった。
• 3次元の場合、平均エネルギー密度の有限サイズ補正はε_N = -1.698 + 2.077N^(-1) と記述され、文献のω ≈ 0.967の推定値と整合的であった。
• 基底状態エネルギー分布の歪度と尖度は、系のサイズが増大するに従い減少し、熱力学的極限においてガウス分布への収束が支持された。
• 小サイズ系では、EAモデル（d=2）のエネルギー分布はGumbel分布に類似していた（m=6）が、Lが増大するにつれてガウス型に近づいた。
• SKモデルの極限分布は依然としてガムベル型に類似しているが、EAモデルではガウス型への遷移が観察され、エネルギー寄与の相関が弱いことが示唆された。
• 確率的アニーリングアルゴリズムは、先行研究と整合する結果を生成し、スピンガラスにおける基底状態エネルギー計算に有効であることが妥当性を確認された。