[论文解读] Optimal Cooperation and Submodularity for Computing Potts' Partition Functions with a Large Number of State

J-Ch. Anglès d’Auriac, Ferenc Iglói+2 more|arXiv (Cornell University)|Apr 2, 2002

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

一句话总结

本文提出了一种多项式时间的组合优化算法——最优协作算法，用于在大q极限下计算具有随机铁磁耦合的q-state Potts模型的配分函数的主导贡献。通过将问题简化为在图上最小化一个子模函数，该方法能够以整数因子为界，精确且高效地计算配分函数，计算时间在L×L晶格上近似按L⁴增长。

ABSTRACT

The partition function of the q-state Potts model with random ferromagnetic couplings in the large-q limit is generally dominated by the contribution of a single diagram of the high temperature expansion. Computing this dominant diagram amounts to minimizing a particular submodular function. We provide a combinatorial optimization algorithm, the optimal cooperation algorithm, which works in polynomial time for any lattice. Practical implementation and the speed of the method is also discussed.

研究动机与目标

• 在大q极限下计算具有随机铁磁耦合的q-state Potts模型的配分函数。
• 解决当q较大且项数呈指数增长时配分函数计算的不可行性问题。
• 为最多包含200,000个自旋位的有限晶格开发一种既精确又高效的计算方法。
• 为Potts模型高温展开的特定结构提供一个组合优化框架。
• 通过精确计算主导图贡献，实现对无序自旋系统中相变与临界现象的研究。

提出的方法

• 将问题简化为从Potts模型配分函数的高温展开中导出的子模函数的最小化。
• 最优协作算法通过逐个添加顶点到不断增长的子图中来逐步构建解，并在每一步保持最优性。
• 在每一步中，通过使用Goldberg-Tarjan算法求解构造网络中的最大流问题。
• 通过在连通分量的生成树中逐步增加边容量来构建网络，并通过显式收缩分量来提升性能。
• 利用图子模函数的结构特性，避免了一般子模最小化方法的开销。
• 该方法使用最大流求解器实现，并通过分量收缩进行优化，相较于通用方法实现了显著加速。

实验结果

研究问题

• RQ1 对于具有随机耦合的有限晶格，能否精确且高效地计算大q极限下Potts模型配分函数的主导贡献？
• RQ2 在高温展开中出现的子模函数最小化问题，是否可在任意晶格上以多项式时间求解？
• RQ3 最优集合的结构（例如在临界点处）如何反映无序系统中相变的底层物理机制？
• RQ4 该算法在晶格尺寸增大时的计算复杂度如何？能否处理包含数十万自旋位的晶格？
• RQ5 该算法是否可被调整以提取有限尺寸标度数据，用于临界指数和相边界分析？

主要发现

• 在大q极限下，q-state Potts模型的配分函数可通过在晶格图上最小化一个子模函数，精确确定到一个整数乘法因子。
• 所提出的最优协作算法可在多项式时间内计算该最小值，从而实现对最多包含200,000个自旋位的晶格的精确计算。
• 对于L×L周期性正方晶格，计算时间近似按L⁴增长，512×512晶格在Pentium III 800 MHz处理器上约需两周时间。
• 与通用子模最小化方法相比，该算法显著更快，其中显式收缩连通分量是最高效的实现策略。
• 在临界点（w₁ + w₂ = 1）处，计算时间约为有序相（w₁ + w₂ > 1）的十倍，以及无序相（w₁ + w₂ < 1）的三倍。
• 该方法可在计算过程中提取所有中间子晶格（如L−1尺寸）的最优集合，从而便于有限尺寸标度分析。