[논문 리뷰] Rapid Mixing of Glauber Dynamics up to Uniqueness via Contraction

Zongchen Chen, Kuikui Liu외 1명|arXiv (Cornell University)|2020. 04. 20.

Markov Chains and Monte Carlo Methods참고 문헌 20인용 수 17

한 줄 요약

이 논문은 관련 감쇠와 다항식 보간에 사용되는 트리 재귀의 수축을 일반적인 2스핀 시스템에서 글라우버 동역학의 빠른 혼합성과 연결하는 통합 프레임워크를 수립한다. 잠재 함수가 트리 재귀에서 수축을 보장할 경우, 유일성 임계점 이하에서는 글라우버 동역학이 다항식 시간 내에 혼합됨을 증명하며, 이는 반자성 및 자성 모델 간의 알고리즘 접근법을 통합한다.

ABSTRACT

For general antiferromagnetic 2-spin systems, including the hardcore model and the antiferromagnetic Ising model, there is an $\mathsf{FPTAS}$ for the partition function on graphs of maximum degree $Δ$ when the infinite regular tree lies in the uniqueness region by Li et al. (2013). Moreover, in the tree non-uniqueness region, Sly (2010) showed that there is no $\mathsf{FPRAS}$ to estimate the partition function unless $\mathsf{NP}=\mathsf{RP}$. The algorithmic results follow from the correlation decay approach due to Weitz (2006) or the polynomial interpolation approach developed by Barvinok (2016). However the running time is only polynomial for constant $Δ$. For the hardcore model, recent work of Anari et al. (2020) establishes rapid mixing of the simple single-site Markov chain known as the Glauber dynamics in the tree uniqueness region. Our work simplifies their analysis of the Glauber dynamics by considering the total pairwise influence of a fixed vertex $v$ on other vertices, as opposed to the total influence on $v$, thereby extending their work to all 2-spin models and improving the mixing time. More importantly our proof ties together the three disparate algorithmic approaches: we show that contraction of the tree recursions with a suitable potential function, which is the primary technique for establishing efficiency of Weitz's correlation decay approach and Barvinok's polynomial interpolation approach, also establishes rapid mixing of the Glauber dynamics. We emphasize that this connection holds for all 2-spin models (both antiferromagnetic and ferromagnetic), and existing proofs for correlation decay or polynomial interpolation immediately imply rapid mixing of Glauber dynamics. Our proof utilizes that the graph partition function divides that of Weitz's self-avoiding walk trees, leading to new tools for analyzing influence of vertices.

연구 동기 및 목표

• 2스핀 시스템에 대한 세 가지 주요 알고리즘 접근법—관련 감쇠, 다항식 보간, 마르코프 체인 몬테카를로(MCMC)—의 분석을 통합하기 위해.
• 모든 2스핀 모델(반자성 및 자성 모두 포함)에 대해 트리 유일성 영역에서 글라우버 동역학의 빠른 혼합성을 확립하기 위해.
• 적절한 잠재 함수를 통해 트리 재귀의 수축이 글라우버 동역학의 다항식 혼합 시간을 암시함을 보여주기 위해.
• 정점별 총 상호영향을 사용한 影향 분석의 단순화를 통해 하드코어 모델에 대한 이전 결과를 일반 2스핀 시스템으로 확장하기 위해.
• 기존의 관련 감쇠 또는 다항식 보간 증명이 동일한 잠재 함수 프레임워크를 통해 글라우버 동역학의 빠른 혼합성을 즉시 암시함을 보여주기 위해.

제안 방법

• 2스핀 시스템에서 트리 재귀의 수축을 분석하기 위한 일반적인 $(\alpha, c)$-잠재 함수 프레임워크를 도입한다.
• 고정된 정점 $v$가 다른 정점에 미치는 총 상호영향을 정의하여, 이전의 $v$에 대한 영향 분석 방식을 대체함으로써 분석을 단순화한다.
• 위츠의 자가피避행렬 트리의 분할 함수가 그래프의 분할 함수를 나누므로, 영향 분석을 위한 새로운 도구를 유도한다.
• 잠재 함수 $\Psi$에 대한 수축 조건을 적용하여 스펙트럼 간격 분석을 통해 글라우버 동역학의 혼합 시간을 경계한다.
• 만약 $\Psi$가 수축 및 유계 조건을 만족하면, 글라우버 동역학은 $O(n^{2+2c/\alpha})$ 시간 내에 혼합됨을 증명한다.
• 특히 $\sqrt{\beta\gamma} \leq \frac{\Delta-2}{\Delta}$ 이고 $\gamma > 1$일 경우, 암시적 $\Psi$ 함수는 $\alpha \geq \delta/2$ 및 $c \leq 18$ (또는 $\beta = 0$일 경우 $c \leq 4$)를 만족하는 일반 $(\alpha, c)$-잠재 함수임을 확립한다.

실험 결과

연구 질문

• RQ1 관련 감쇠와 다항식 보간에 사용되는 트리 재귀의 수축이 글라우버 동역학의 빠른 혼합성과 직접적으로 연결될 수 있는가?
• RQ2 정점이 다른 정점에 미치는 총 상호영향이 글라우버 동역학 혼합 시간 분석을 위한 더 단순하고 일반적인 프레임워크를 제공하는가?
• RQ3 관련 감쇠나 보간 효율성을 보장하는 동일한 잠재 함수가 글라우버 동역학의 빠른 혼합성도 보장하는가?
• RQ4 일반 2스핀 시스템에서 트리 유일성 영역에서 글라우버 동역학의 혼합 시간은 무엇이며, 이는 유일성 갭 $\delta$에 어떻게 의존하는가?
• RQ5 결과들이 반자성 및 자성 2스핀 모델, 하드코어 모델 및 이징 모델을 포함하여 어떻게 확장되는가?

주요 결과

• 일반 $(\alpha, c)$-잠재 함수가 존재할 경우, 일반 2스핀 시스템의 글라우버 동역학은 $O(n^{2+2c/\alpha})$ 시간 내에 혼합된다.
• $\sqrt{\beta\gamma} \leq \frac{\Delta-2}{\Delta}$ 이고 $\gamma > 1$일 경우, 암시적 잠재 함수 $\Psi$는 $\alpha \geq \delta/2$ 및 $c \leq 18$를 만족하는 일반 $(\alpha, c)$-잠재 함수이며, 이는 이전의 경계를 향상시킨다.
• $\beta = 0$일 경우, 동일한 잠재 함수는 $c \leq 4$를 제공하여 하드코어 모델에 대해 더 탴지 혼합 시간 경계를 암시한다.
• 항등 함수 $\Psi(y) = y$는 특정 조건 하에서 자성 2스핀 시스템에 대해 $(\alpha, c)$-잠재 함수이며, $\alpha = \Theta(\delta)$ 및 $c \leq O(1)$를 만족하여 $O(n^{2+c/\delta})$ 혼합 시간을 이끈다.
• 자성 시스템에서 $\sqrt{\beta\gamma} \geq \frac{\Delta}{\Delta-2}$ 이고 $\lambda$가 유계일 경우, 혼합 시간은 $\beta, \gamma, \lambda, \Delta$에 따라 달라지는 상수 $C$에 대해 $O(n^C)$이다.
• 이 프레임워크는 기존의 관련 감쇠 및 보간 결과를 통합한다: 잠재 함수를 통한 효율적 관련 감쇠 또는 보간 증명은 즉시 동일한 프레임워크를 통해 글라우버 동역학의 빠른 혼합성을 암시한다.