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

一句话总结

本文研究了聚类算法（特别是二维伊辛-塔勒模型的Swendsen–Wang及其相关算法）的动力临界行为，发现Li–Sokal对自相关时间的下界并不严格成立，但几乎成立，自相关时间积分与比热之比随系统尺寸对数增长或以较小幂次（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}}.

研究动机与目标

• 理解Swendsen–Wang算法为何在某些模型中表现良好，而在其他模型中表现不佳，特别是在临界点附近。
• 检验Li–Sokal对自相关时间的下界是否严格成立，该下界将动态临界指数与静态临界指数关联起来。
• 确定在二维聚类算法中，能量的指数自相关时间与积分自相关时间是否以相同的动态临界指数缩放。
• 研究非局部蒙特卡洛算法动力学中慢模式的性质。

提出的方法

• 在自 dual 曲线上对二维 q 状态 Potts 模型（q=2,3,4）和对称 Ashkin–Teller 模型进行 Swendsen–Wang 及其相关聚类算法的数值蒙特卡洛模拟。
• 利用长时间模拟的高精度数据，估计能量可观测量的积分自相关时间和指数自相关时间。
• 应用有限尺寸标度分析，通过自相关时间随系统尺寸 L 的缩放关系提取动态临界指数。
• 使用标度假设 ρ_AA(t;L) ≈ |t|^{-p_A} h_A[t/τ_exp,A; ξ(L)/L] 检验能量自相关函数是否表现为纯指数函数。
• 通过对数和幂律拟合形式，比较 τ_int,ℰ / C_H 与 Li–Sokal 下界的关系，以评估其严格性。
• 通过在不同晶格尺寸下检查 ρ_ℰℰ(t) 与 t/τ_int,ℰ 的数据坍缩，验证 τ_int,ℰ ≈ τ_exp,ℰ 的假设。

实验结果

研究问题

• RQ1 在二维系统中，聚类算法的 Li–Sokal 自相关时间下界是否严格成立？
• RQ2 在 Swendsen–Wang 算法中，能量的指数自相关时间与积分自相关时间是否以相同的动态临界指数缩放？
• RQ3 比值 τ_int,ℰ / C_H 如何随系统尺寸 L 缩放？其增长是否为对数或小幂次？
• RQ4 二维 Potts 模型与 Ashkin–Teller 模型动力学中最慢的模式具有何种性质？
• RQ5 为何 Li–Sokal 下界在二维中几乎严格成立，而在高维中明显不严格？

主要发现

• Li–Sokal 下界 τ_int,ℰ ≥ const × C_H 并不严格成立，但比值 τ_int,ℰ / C_H 随系统尺寸 L 增大而增长，表现为对数增长或较小幂次（0.05 ≤ p ≤ 0.12）。
• 在二维 Potts 模型（q=2,3,4）和 Ashkin–Teller 模型中，能量的动态临界指数 z_int,ℰ 与 z_exp,ℰ 一致，表明指数自相关时间与积分自相关时间以相同指数缩放。
• 对于二维伊辛模型（q=2），比值 τ_int,ℰ / C_H 呈对数增长，表明算法的动力临界标度近乎最优。
• 在 τ_int,ℰ / C_H ∼ L^p 的缩放关系中，幂次 p 从伊辛模型中的约 0.05 增加到四状态 Potts 模型中的 0.12，表明在不同模型间存在连续趋势。
• 对于 L ≥ 64，ρ_ℰℰ(t) 与 t/τ_int,ℰ 的数据坍缩效果极佳，支持能量自相关函数表现为纯指数函数的假设。
• 研究表明，在 d=2 时，能量和键占据数是响应最慢的模式，而在 d>2 时则不是，这解释了 Li–Sokal 下界在二维中近乎严格而高维中不严格的原因。