alphaFeedback

WorksMathematics

Local-to-Global Theorems Based on Cones with Applications to Graph Coloring

AI contributions
  • AI
View details

Abstract

Local-to-Global theorems are a cornerstone of many of recent developments in theoretical computer science from coding theory and PCPs all the way to analyzing mixing time of Markov chains. Building on trickledown theorems on simplicial complexes, several techniques have been established to obtain tighter local-to-global theorems to analyze the spectral gap of the Glauber dynamics, namely matrix trickledown, C{\mathcal C}-Lorentzian polynomials and more recently Bochner type inequalities. In this paper, we adopt non-local local-to-global/Bochner type inequalities on "cones", corresponding to the local neighborhood of a given vertex, to show that the Glauber dynamics mixes in polynomial time to sample proper colorings of maximum degree ΔΔ graphs as long as q≥1.7408Δq\geq 1.7408Δ.

Introduction

Let G=(V,E)G=(V,E) be a graph with nn vertices and maximum degree Δ>1\Delta>1 and a family of list of colors {Li}i∈G\{L_i\}_{i\in G}. A proper coloring is an assignment x∈∏iLi\mathbf{x}\in\prod_i L_i such that for any two adjacent vertices i∼ji\sim j, we have xi≠xjx_i\ne x_j. Let μ\mu be the uniform distribution over all proper colorings. One can run a simple Markov chain to generate a sample from μ\mu: Start with an arbitrary proper coloring x\mathbf{x}, each time choose a uniformly random vertex ii and re-sample its color given the colors of its neighbors. Let 1nP\frac{1}{n}\mathcal{P} be the transition probability matrix of the Glauber dynamics Markov chain (the factor 1/n1/n is for notational convenience as will be clear later; essentially it implies every vertex is chosen at a “unit rate”) then the mixing time is defined to be

max⁡xmin⁡{t:∥1xT(1nP)t−μ∥TV≤14}.\max_{\mathbf{x}}\min\left\{t:\left\|\mathbf{1}_{\mathbf{x}}^{\mathsf{T}}\left(\frac{1}{n}\mathcal{P}\right)^t-\mu\right\|_{\mathrm{TV}}\le\frac{1}{4}\right\}.

To bound the mixing time of 1nP\frac{1}{n}\mathcal{P} a method of choice is to bound its spectral gap.

Theorem 1.1. If the second largest eigenvalue of 1nP\frac{1}{n}\mathcal{P} is at most 1−ε1-\varepsilon, or equivalently the second smallest eigenvalue of the Laplacian matrix 1nD=I−1nP\frac{1}{n}\mathcal{D}=I-\frac{1}{n}\mathcal{P} is at least ε\varepsilon, then the mixing time of 1nP\frac{1}{n}\mathcal{P} is at most O(n⋅log⁡(1μmin⁡)/ε)O\left(n\cdot\log\left(\frac{1}{\mu_{\min}}\right)/\varepsilon\right) where μmin⁡=min⁡xμ(x)\mu_{\min}=\min_{\mathbf{x}}\mu(\mathbf{x}).

Suppose that the size of every list is at least qq. It is a long standing open problem that if q≥Δ+2q \ge\Delta+ 2, then the Glauber dynamics mentioned above mixes rapidly.

Conjecture 1.2. For any graph GG with maximum degree Δ\Delta, if q≥Δ+2q \ge\Delta+ 2, then the Glauber dynamics mixes in time polynomial in nn.

In the early 1990s, Jerrum [19] showed that q≥2Δ+1q \ge2\Delta+ 1 suffices via a simple coupling approach. Vigoda [25] improved this to 11Δ/611\Delta/6 using a more sophisticated coupling. Later in a sequence of works [CDM⁺19, JMV26, CV25] the bound was improved to (11/6−ε)Δ(11/6-\varepsilon)\Delta using various more involved coupling techniques.

Local to Global Theorems. With the introduction of the spectral independence technique [3] a new challenge started: To come up with new local-to-global arguments to analyze Markov chains. It is worth noting that even though classical proof techniques such as the correlation decay or path coupling methods could be translated to spectral independence [3, 12, 21 BCC⁺22] it seemed the main strength of the spectral independence machinery stems from the local-to-global technique known as the trickledown theorem [24] which was used to count bases of matroids [4].

Towards this goal, analyzing Glauber dynamics for graph coloring was considered as a benchmark. The first improvement over the classical trickledown technique came from the Matrix trickledown theorems [2, 1 CCC⁺25] which were used to obtain near optimal bounds for edge coloring problems [WZZ24], or to sample from the Sherrington-Kirkpatrick model [1]. Later, this was extended to the CC-Lorentzian polynomial framework which led to a “spectral” Dobrushin’s criterion [22, 23] that itself lead to improved analysis of the multi-state hardcore model or improved analysis of the Hard core model on planar graphs [16]. Later on, Göbel, Jenssen, Michelen, Pappik, Perkins and Schiller [GJM⁺26] observed that the matrix trickledown of [CCC⁺25] can be seen as a discrete Bochner inequality. This led to several improvement for various problems [11 BLO26, 10]. In this paper, we establish the discrete Bochner inequality over O(Δ)O(\Delta)-sized subgraphs and then use it to analyze the Glauber dynamics for sampling proper colorings of graphs.

Background on Bochner Inequality. The discrete Bochner type inequalities were recently recently used in [GJM⁺26] to re-prove the trickledown theorem of Chen, Chen, Chen, Yin, and Zhang [CCC⁺25]. For background on this technique see [6 Equation 1.16.6, Proposition 4.8.3] and references therein. Bochner-type methods have also been used to establish spectral-gap and mixing-time bounds in a variety of settings, including the Sherrington-Kirkpatrick model [Wan26, BLO26], matroid intersection [10], and several recent works by OpenAI.

Let PifP_i f be the function which averages ff on Li\mathcal{L}_i given the color of all other vertices, namely:

Pif(x)=Eμ[f(X1,…,Xn)∣Xj=xj,∀j≠i].P_i f(\mathbf{x}) = \mathbb{E}_{\mu}[f(X_1,\ldots,X_n)\mid X_j = x_j,\forall j \ne i].

Throughout the paper we use bold letters to denote vectors.

Then, we can write P=∑iPiP = \sum_i P_i. In other words the PiP_i operators give a decomposition of PP into projection matrices; and in most recent applications this has been perhaps the main starting point of Bochner type inequalities. Then, it follows that the transition probability matrix of the Glauber dynamics is 1nP\frac{1}{n}P. We let Di=I−PiD_i=I-P_i and D=∑iDiD=\sum_i D_i. It is not hard to see that the operators Pi,DiP_i,D_i are projections, i.e., for any function ff, Pi2f=PifP_i^2f=P_if. Define the inner-product

⟨f,g⟩μ=EX∼μ[f(X)g(X)].\langle f,g\rangle_\mu=\mathbb{E}_{X\sim\mu}[f(X)g(X)].

It is not hard to see that Pi,DiP_i,D_i are self-adjoint with respect to this inner-product, i.e., for two functions f,gf,g, ⟨f,Pig⟩=⟨Pif,g⟩\langle f,P_ig\rangle=\langle P_if,g\rangle, and since Di=I−PiD_i=I-P_i, the self-adjointness of DiD_i follows.

Fact 1.3. If for any function f:Ω→Rf:\Omega\to\mathbb{R},

⟨Df,Df⟩μ≥ε⋅⟨Df,f⟩μ=ε⋅∑i⟨Dif,f⟩μ,(1)\langle Df,Df\rangle_\mu\geq\varepsilon\cdot\langle Df,f\rangle_\mu=\varepsilon\cdot\sum_i\langle D_if,f\rangle_\mu, \tag*{(1)}

then the spectral gap of the Glauber dynamics 1nP\frac{1}{n}P is at least ε/n\varepsilon/n.

Proof. Given an arbitrary function gg, since D⪰0D\succeq0, we let f=D−1/2gf=D^{-1/2}g (to be precise, we need to take the pseudo-inverse of DD since DD has a 0 eigenvalue). Having this,

⟨Dg,g⟩=⟨Df,Df⟩≥ε⋅⟨Df,f⟩=ε⋅⟨g,g⟩,\langle Dg,g\rangle=\langle Df,Df\rangle\geq\varepsilon\cdot\langle Df,f\rangle=\varepsilon\cdot\langle g,g\rangle,

where all the inner-products are with respect to μ\mu. This concludes the proof.

Top Link Spectral Expansion. A natural approach to prove (1) is to lower-bound ⟨Dif,Djf⟩μ\langle D_if,D_jf\rangle_\mu for i,j∈Gi,j\in G. This is very analogous to analyzing links of co-dimension 2 in simplicial complexes, see e.g., [4]. See the following lemma and the proof.

Lemma 1.4. Given a graph HH with vertices 1 and 2. If HH has no edges then ⟨D1f,D2f⟩μ≥0\langle D_1f,D_2f\rangle_\mu\geq0. Otherwise, suppose s≥5s\geq5 and ∣L1∣,∣L2∣≥s+1|L_1|,|L_2|\geq s+1. Then for every f∈L2(μ)f\in L^2(\mu), we have

2⟨D1f,D2f⟩μ≥−1s(∥D1f∥μ2+∥D2f∥μ2).2\langle D_1f,D_2f\rangle_\mu\geq-\frac{1}{s}\left(\lVert D_1f\rVert_\mu^2+\lVert D_2f\rVert_\mu^2\right).

Proof. If HH has no edges then μ\mu is a product measure, so D1,D2D_1,D_2 are commuting orthogonal projections:

⟨D1f,D2f⟩μ=⟨f,D1D2f⟩μ=∥D1D2f∥μ2≥0.\langle D_1f,D_2f\rangle_\mu=\langle f,D_1D_2f\rangle_\mu=\lVert D_1D_2f\rVert_\mu^2\geq0.

Now suppose HH has an edge. We let PP be the simple random walk on the bipartite graph with classes {1}×L1\{1\}\times L_1 and {2}×L2\{2\}\times L_2, where (1,a)(1,a) and (2,b)(2,b) are adjacent exactly when a≠ba\ne b. By [11 Eq. (4)],

2⟨D1f,D2f⟩μ≥−λ2(P)(∥D1f∥μ2+∥D2f∥μ2),2\langle D_1f,D_2f\rangle_\mu\geq-\lambda_2(P)\left(\lVert D_1f\rVert_\mu^2+\lVert D_2f\rVert_\mu^2\right),

where λ2(P)\lambda_2(P) denotes the second largest eigenvalue. It follows from [2 proof of Proposition IV.4] that λ2(P)≤1(ℓ1−1)(ℓ2−1)≤1/s\lambda_2(P)\leq\frac{1}{\sqrt{(\ell_1-1)(\ell_2-1)}}\leq1/s.

Having this, we can obtain that the Glauber dynamics has a spectral gap at least 1/Δn1/\Delta n if q>2Δ+2q>2\Delta+2. In particular

⟨Df,Df⟩≥∑i⟨Dif,f⟩μ+∑{i,j}∈E(G)2⟨Dif,Djf⟩μ≥∑i(1−deg⁡(i)q−Δ)⋅⟨Dif,f⟩μ≥1Δ⟨Df,f⟩μ.\langle Df,Df\rangle\geq\sum_i\langle D_if,f\rangle_\mu+\sum_{\{i,j\}\in E(G)}2\langle D_if,D_jf\rangle_\mu\geq\sum_i\left(1-\frac{\deg(i)}{q-\Delta}\right)\cdot\langle D_if,f\rangle_\mu\geq\frac{1}{\Delta}\langle Df,f\rangle_\mu.

References

  1. [1][AKV24] Nima Anari, Frederic Koehler, and Thuy-Duong Vuong. Trickle-down in localization schemes and applications. In STOC’24—Proceedings of the 56th Annual ACM Symposium on Theory of Computing, pages 1094–1105. ACM, New York, 2024.
  2. [2][ALO22] Dorna Abdolazimi, Kuikui Liu, and Shayan Oveis Gharan. A matrix trickle-down theorem on simplicial complexes and applications to sampling colorings. In 2021 IEEE 62nd Annual Symposium on Foundations of Computer Science—FOCS 2021, pages 161–172. IEEE Computer Soc., Los Alamitos, CA, 2022.arxiv.org/abs/2106.03845
  3. [3]Nima Anari, Kuikui Liu, and Shayan Oveis Gharan. Spectral independence in high-dimensional expanders and applications to the hardcore model. SIAM J. Comput., 53(6):FOCS20–1–FOCS20–37, 2024.
  4. [4]Nima Anari, Kuikui Liu, Shayan Oveis Gharan, and Cynthia Vinzant. Log-concave polynomials II: High-dimensional walks and an FPRAS for counting bases of a matroid. Ann. of Math. (2), 199(1):259–299, 2024.
  5. [6]Dominique Bakry, Ivan Gentil, and Michel Ledoux. Analysis and geometry of Markov diffusion operators, volume 348 of Grundlehren der mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer, Cham, 2014.DOI
  6. [10]Xiaoyu Chen and Kuikui Liu. An FPRAS for Counting Common Bases of Two Matroids. arXiv preprint arXiv:2610.06724, 2026.arxiv.org/abs/2610.06724
  7. [11]Xiaoyu Chen and Kuikui Liu. Spectral gap of down-up walks via trickle-down: A simplified and sharpened analysis. arXiv preprint arXiv:2609.19514, 2026.arxiv.org/abs/2609.19514
  8. [12]Zongchen Chen, Kuikui Liu, and Eric Vigoda. Rapid mixing of Glauber dynamics up to uniqueness via contraction. SIAM J. Comput., 52(1):196–237, 2023.
  9. [14]Martin Dyer, Alan Frieze, and Mark Jerrum. On counting independent sets in sparse graphs. In 40th Annual Symposium on Foundations of Computer Science (New York, 1999), pages 210–217. IEEE Computer Soc., Los Alamitos, CA, 1999.
  10. [16]Heng Guo and Xinyuan Zhang. Approximating spin systems on planar graphs. arXiv preprint arXiv:2608.06172, 2026.arxiv.org/abs/2608.06172
  11. [17]Alireza Haqi. Sampling line-graph colorings with constant extra colors. arXiv preprint arXiv:2609.27440, 2026.arxiv.org/abs/2609.27440
  12. [18]Thomas P. Hayes and Eric Vigoda. Coupling with the stationary distribution and improved sampling for colorings and independent sets. Ann. Appl. Probab., 16(3):1297–1318, 2006.DOI
  13. [19]Mark Jerrum. A very simple algorithm for estimating the number of k-colorings of a low-degree graph. Random Structures Algorithms, 7(2):157–165, 1995.DOI
  14. [21]Kuikui Liu. From coupling to spectral independence and blackbox comparison with the down-up walk. In Approximation, randomization, and combinatorial optimization. Algorithms and techniques, volume 207 of LIPICs. Leibniz Int. Proc. Inform., pages Art. No. 32, 21. Schloss Dagstuhl. Leibniz-Zent. Inform., Wadern, 2021.
  15. [22]Jonathan Leake, Kasper Lindberg, and Shayan Oveis Gharan. Optimal trickle-down theorems for path complexes via C-Lorentzian polynomials with applications to sampling and log-concave sequences. In 2025 IEEE 66th Annual Symposium on Foundations of Computer Science—FOCS 2025, pages 1857–1870. IEEE Comput. Soc. Press, Los Alamitos, CA, 2025.arxiv.org/abs/2503.01005
  16. [23]Jonathan Leake and Shayan Oveis Gharan. Trickle-down theorems via C-Lorentzian polynomials II: Pairwise spectral influence and improved Dobrushin’s condition. arXiv preprint arXiv:2510.06549, 2025.arxiv.org/abs/2510.06549
  17. [24]Izhar Oppenheim. Local spectral expansion approach to high dimensional expanders Part I: Descent of spectral gaps. Discrete Comput. Geom., 59(2):293–330, 2018.
  18. [25]Eric Vigoda. Improved bounds for sampling colorings. J. Math. Phys., 41(3):1555–1569, 2000. Probabilistic techniques in equilibrium and nonequilibrium statistical physics.

Paper details

vlnc:W_qr7ewazh
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 1

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

AI · View statement →

Looking for related papers…
Contents