[论文解读] The computational hardness of counting in two-spin models on d-regular graphs

Allan Sly, Nike Sun|arXiv (Cornell University)|Mar 12, 2012

Markov Chains and Monte Carlo Methods参考文献 17被引用 60

一句话总结

该论文证明了在 d-正则图上，当系统在 d-正则树上表现出非唯一性时，对自旋系统（如硬核模型和反铁磁伊辛模型）的划分函数进行近似是 NP-难的。作者通过一种新颖的随机化归约方法将问题归约为最大割问题，利用局部树状结构的扩展图，并证明自由能密度收敛于贝特定律预测值，从而绕过了传统的二阶矩方法技术。

ABSTRACT

The class of two-spin systems contains several important models, including random independent sets and the Ising model of statistical physics. We show that for both the hard-core (independent set) model and the anti-ferromagnetic Ising model with arbitrary external field, it is NP-hard to approximate the partition function or approximately sample from the model on d-regular graphs when the model has non-uniqueness on the d-regular tree. Together with results of Jerrum--Sinclair, Weitz, and Sinclair--Srivastava--Thurley giving FPRAS's for all other two-spin systems except at the uniqueness threshold, this gives an almost complete classification of the computational complexity of two-spin systems on bounded-degree graphs. Our proof establishes that the normalized log-partition function of any two-spin system on bipartite locally tree-like graphs converges to a limiting "free energy density" which coincides with the (non-rigorous) Bethe prediction of statistical physics. We use this result to characterize the local structure of two-spin systems on locally tree-like bipartite expander graphs, which then become the basic gadgets in a randomized reduction to approximate MAX-CUT. Our approach is novel in that it makes no use of the second moment method employed in previous works on these questions.

研究动机与目标

• 对有界度数图上两自旋系统的划分函数近似计算复杂性的分类。
• 解决关于 d-正则树上唯一性阈值是否为两自旋模型计算复杂性转变的确切边界的开放问题。
• 提出一种新证明技术，避免先前硬度结果中使用的二阶矩方法。
• 将硬度结果扩展到具有任意外场的反铁磁伊辛模型。
• 完成对 d-正则图上同质两自旋系统计算复杂性的分类，除唯一性阈值外均已完成。

提出的方法

• 构建作为硬度归约工具的双分图局部树状扩展图。
• 证明任意两自旋系统在双分图局部树状图上的归一化对数划分函数收敛于与统计物理中贝特定律预测一致的极限自由能密度。
• 利用自由能收敛性来表征局部树状双分图扩展图上两自旋系统的局部结构。
• 设计一种从近似最大割到近似两自旋模型划分函数的随机化归约。
• 通过将 3-正则图 $ H $ 的顶点替换为基图 $ G $ 的副本，并在对应边界集之间连接 $ 2k $ 条边，定义图构造 $ H^G $。
• 建立 $ Z_{H^G} $ 与 $ Z_{\widehat{H}^G} $ 之间的比值界限，该界限通过两自旋模型导出的参数 $ \Gamma $ 和 $ \Theta $ 将划分函数与 $ H $ 的最大割联系起来。

实验结果

研究问题

• RQ1 d-正则树上的唯一性阈值是否正是 d-正则图上两自旋模型可近似与不可近似之间的精确边界？
• RQ2 是否可以不依赖于二阶矩方法计算来建立划分函数近似计算的计算硬度？
• RQ3 贝特定律对自由能密度的预测是否准确描述了局部树状图上两自旋系统的渐近行为？
• RQ4 能否利用扩展图上的两自旋模型构造从最大割到近似的随机化归约，以证明近似硬度？
• RQ5 具有任意外场的反铁磁伊辛模型的计算转变点是否也位于 d-正则树上的唯一性阈值处？

主要发现

• 对于 $ d \geq 3 $ 且 $ \lambda > \lambda_c(d) = \frac{(d-1)^{d-1}}{(d-2)^d} $，除非 NP = RP，否则在 d-正则图上硬核模型不存在 FPRAS。
• 对于 $ d \geq 3 $，$ B \in \mathbb{R} $，且 $ \beta < \beta_{c,\text{af}}(B,d) $，除非 NP = RP，否则在 d-正则图上具有外场 $ B $ 的反铁磁伊辛模型不存在 FPRAS。
• 任意两自旋系统在双分图局部树状图上的归一化对数划分函数收敛于等于贝特定律预测的极限自由能密度。
• 作者证明了两自旋模型在 $ H^G $ 上的划分函数与 $ H $ 的最大割紧密相关，其比值被 $ (1 \pm \epsilon)^m $ 所限制，从而实现了随机化归约。
• 该证明避开了二阶矩方法，相比以往依赖精细二阶矩计算的工作，提供了更具概念性和普适性的方法。
• 结果完成了对 d-正则图上同质两自旋系统计算复杂性的分类，所有非唯一性区域的硬度均已确立，所有唯一性区域的 FPTAS 已知。