[论文解读] Optimal Mixing Time for the Ising Model in the Uniqueness Regime

Xiaoyu Chen, Weiming Feng+2 more|arXiv (Cornell University)|Nov 4, 2021

Markov Chains and Monte Carlo Methods参考文献 25被引用 4

一句话总结

本文在一般图的伊辛模型唯一性区域中，建立了格罗伯动力学的最优 $O(n\log n)$ 混合时间，即使最大度数 $\Delta$ 无界。通过改进利用熵独立性和新型修正对数 Sobolev（MLS）常数定理的增强技术，作者获得了与理论下界匹配的紧致界，解决了统计物理与 MCMC 采样领域长期存在的开放问题。

ABSTRACT

We prove an optimal $O(n \log n)$ mixing time of the Glauber dynamics for the Ising models with edge activity $β\in \left(\frac{Δ-2}Δ, \fracΔ{Δ-2} ight)$. This mixing time bound holds even if the maximum degree $Δ$ is unbounded. We refine the boosting technique developed in [CFYZ21], and prove a new boosting theorem by utilizing the entropic independence defined in [AJK+21]. The theorem relates the modified log-Sobolev (MLS) constant of the Glauber dynamics for a near-critical Ising model to that for an Ising model in a sub-critical regime.

研究动机与目标

• 解决已知多项式混合时间界与伊辛模型格罗伯动力学在唯一性区域中理论 $\Omega(n\log n)$ 下界之间差距的开放问题。
• 建立一个适用于最大度数 $\Delta$ 无界的任意图上铁磁与反铁磁伊辛模型的最优 $O(n\log n)$ 混合时间界。
• 通过利用熵独立性和新的 MLS 常数定理，改进先前工作的增强技术，以获得更紧致的界。
• 提供统一的分析方法，同时处理常数局部场及彼此间相差常数因子的场，优于以往的 $n^{O(1/\delta)}$ 界。
• 将谱独立性和基于 MLS 的方法的应用扩展至唯一性区域内部，包括临界附近参数。

提出的方法

• 提出一个新的增强定理，通过熵独立性将近临界伊辛模型的修正对数 Sobolev（MLS）常数与子临界模型的 MLS 常数联系起来。
• 利用 $\alpha$-谱独立性与多勃鲁什金影响矩阵，通过边际下界 $\alpha$ 和谱范数 $\|A\|_2$ 来界定 MLS 常数。
• 引入翻转分布 $\nu = \text{flip}(\mu, \chi_I)$ 以稳定边际界，并确保所有顶点上的边际概率具有统一的下界。
• 应用引理 8.5，利用 $\alpha \geq C / (2 \times 10^4)$ 和 $\|A\|_2 \leq 3/5$ 来下界 MLS 常数，从而得到 $\Omega(1/n)$ 的 MLS 常数。
• 利用关系式 $\rho_{\text{GD}}(\pi_\tau) = \frac{m}{n} \rho_{\text{GD}}(\pi_\tau^\Lambda)$，将子图的 MLS 常数扩展至整个系统。
• 通过 $T_{\text{mix}}(\varepsilon) \leq \frac{1}{\rho_{\text{GD}}(\mu)} \left( \log \log \frac{1}{\mu_{\min}} + \log \frac{1}{2\varepsilon^2} \right)$ 建立紧致的混合时间界，其中 $\mu_{\min}$ 以 $\lambda_{\min}/\lambda_{\max}$ 的形式有下界。

实验结果

研究问题

• RQ1 在最大度数 $\Delta$ 无界的情况下，伊辛模型格罗伯动力学在唯一性区域中的最优混合时间是什么？
• RQ2 对于 $\beta \in \left[ \frac{\Delta-2+\delta}{\Delta-\delta}, \frac{\Delta-\delta}{\Delta-2+\delta} \right]$，能否将 $n^{O(1/\delta)}$ 的混合时间界改进为 $O(n\log n)$？
• RQ3 如何利用熵独立性和修正对数 Sobolev 常数，实现对近临界伊辛模型的紧致增强定理？
• RQ4 混合时间界是否依赖于比值 $\lambda_{\max}/\lambda_{\min}$？若依赖，如何在分析中控制该依赖？
• RQ5 该分析能否扩展至获得一个 $O(n\log^2 n)$ 的采样算法，且常数依赖于 $\delta$？

主要发现

• 本文证明了在唯一性区域中，伊辛模型格罗伯动力学的混合时间达到最优 $O(n\log n)$，与理论下界完全匹配。
• 对于任意 $\delta \in (0,1)$，混合时间满足 $T_{\text{mix}}(\varepsilon) \leq C_\delta \cdot \frac{\lambda_{\max}}{\lambda_{\min}} \cdot n \left( \log \frac{n}{\varepsilon} + \log \log \frac{2\lambda_{\max}}{\lambda_{\min}} \right)$，其中 $C_\delta = \exp(O(1/\delta))$。
• 该界对最大度数 $\Delta$ 无界的任意图上的铁磁与反铁磁伊辛模型均成立。
• 证明了修正对数 Sobolev（MLS）常数至少为 $\Omega(1/n)$，这直接通过标准界推出 $O(n\log n)$ 的混合时间。
• 通过利用熵独立性的改进增强定理，实现了从子临界模型到近临界模型的 MLS 界转移，克服了以往的局限性。
• 通过结合当前结果与场动力学，展示了可实现一个 $O_\delta(n\log^2 n)$ 的采样算法，且常数仅依赖于 $\delta$。