[论文解读] 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

一句话总结

本文提出一种随机退火算法，用于研究自旋玻璃体系中基态能量涨落，重点关注Sherrington-Kirkpatrick（SK）和Edwards-Anderson（EA）模型。研究证实了能量涨落在不同维度下的幂律标度行为，其中EA模型（d=2）的指数ρ ≈ 0.501，d=2时的ω ≈ 1.25；并发现基态能量分布的极限形状在热力学极限下趋于高斯分布，与有限尺寸SK系统中观察到的Gumbel型行为形成对比。

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.

研究动机与目标

• 研究具有 quenched disorder 的自旋玻璃模型中基态能量涨落的有限尺寸标度行为。
• 验证一种随机退火算法在复杂自旋玻璃体系中计算基态的可行性。
• 确定标准偏差和能量密度均值的标度指数ρ和ω。
• 分析在热力学极限下基态能量概率分布p_N(e_N)的极限形状。

提出的方法

• 采用一种兼具贪婪与保守策略的随机退火算法，以在自旋玻璃哈密顿量中寻找全局能量最小值。
• 该算法被应用于平均场的Sherrington-Kirkpatrick（SK）模型和短程相互作用的Edwards-Anderson（EA）模型。
• 通过大量独立实现的quenched耦合J的 disorder 平均，计算能量密度涨落。
• 采用幂律拟合进行有限尺寸标度：ε_N = e_∞ + bN^(-ω) 和 σ_N = aN^(-ρ)，其中N为体系尺寸。
• 将分布p_N(e_N)重新标度为x_N = (e_N - ε_N)/σ_N，以研究其在大N极限下的极限形式。
• 计算统计矩（偏度与峰度）以评估其向高斯分布收敛的趋势。

实验结果

研究问题

• RQ1 在Edwards-Anderson模型中，基态能量的标准差σ_N和能量密度均值ε_N如何随体系尺寸N变化？
• RQ2 在热力学极限下，EA模型中基态能量密度p_N(e_N)的概率分布是否收敛于高斯分布？
• RQ3 EA模型（d=2和d=3）中，能量均值及其涨落的有限尺寸修正指数ω和ρ分别为多少？
• RQ4 随着体系尺寸增大，p_N(e_N)的形状如何演化？是否如极值理论所预测的，从Gumbel型行为过渡到高斯行为？
• RQ5 该随机退火算法在不同自旋玻璃模型中是否能有效且准确地计算基态能量？

主要发现

• 在二维Edwards-Anderson模型中，基态能量的标准差按σ_N ∝ N^(-0.501) 标度，最佳拟合指数为ρ = 0.501 ± 6×10⁻³。
• 在三维情况下，能量密度均值的有限尺寸修正满足ε_N = -1.698 + 2.077N^(-1)，与文献中估计的ω ≈ 0.967 一致。
• 基态能量分布的偏度与峰度随体系尺寸减小，支持其在热力学极限下向高斯分布收敛。
• 对于小尺寸体系，EA模型（d=2）的能量分布呈现Gumbel分布特征（m=6），但随着L增大，逐渐趋向高斯分布。
• SK模型的极限分布仍保持Gumbel型，而EA模型则表现出向高斯分布的转变，表明其能量贡献具有更弱的相关性。
• 随机退火算法的结果与先前研究一致，验证了其在自旋玻璃体系中计算基态能量的可靠性。