Rapid Mixing of Glauber Dynamics up to Uniqueness via Contraction

Abstract

For general antiferromagnetic $2$-spin systems, including the hardcore model on weighted independent sets and the antiferromagnetic Ising model, there is an $\mathrm{𝖥𝖯𝖳𝖠𝖲}$ for the partition function on graphs of maximum degree $\mathrm{\Delta }$ 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 $\mathrm{𝖥𝖯𝖱𝖠𝖲}$ to estimate the partition function unless $\mathrm{𝖭𝖯}=\mathrm{𝖱𝖯}$. 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 $\mathrm{\Delta }$. 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 of other vertices 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 so-called 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 the correlation decay or polynomial interpolation approach immediately imply rapid mixing of the Glauber dynamics. Our proof utilizes that the graph partition function is a divisor of the partition function for Weitz’s self-avoiding walk tree. This fact leads to new tools for the analysis of the influence of vertices, and may be of independent interest for the study of complex zeros.

1 Introduction

A remarkable connection has been established between the computational complexity of approximate counting problems in general graphs of maximum degree $\mathrm{\Delta }$ and the statistical physics phase transition on infinite, regular trees of degree $\mathrm{\Delta }$ (or up to $\mathrm{\Delta }$ in the more general case). This connection holds for 2-state antiferromagnetic spin systems – the hardcore model on independent sets and the Ising model are the most interesting examples of such systems.

Given an $n$-vertex graph $G=(V,E)$, configurations of the 2-spin model are the $2^n$ assignments of spins $0,1$ to the vertices. A 2-spin system is defined by three parameters: edge weights $\beta ,\gamma >0$ and a vertex weight $\lambda >0$. Edge parameter $\beta$ controls the (relative) strength of interaction between neighboring $1$-spins, $\gamma$ corresponds to neighboring $0$-spins, and $\lambda$ is the external field applied to vertices with $1$-spins.

Every spin configuration $\sigma \in {\{0,1\}}^V$ is assigned a weight

$$w_G(\sigma )={\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\lambda }^{n_1(\sigma )},$$

where, for spin $s\in \{0,1\}$, $m_s(\sigma )=\mathrm{\#}\{uv\in E:{\sigma }_u={\sigma }_v=s\}$ is the number of monochromatic edges with spin $s$, and $n_1(\sigma )=\mathrm{\#}\{v\in V:{\sigma }_v=1\}$ is the number of vertices with spin $1$ (as is standard, the parameters are normalized so we can avoid two additional parameters). The Gibbs distribution over spin configurations is given by ${\mu }_G(\sigma )=\frac{w_G(\sigma )}{Z_G(\beta ,\gamma ,\lambda )},$ where $Z_G(\beta ,\gamma ,\lambda )={\sum }_{\sigma \in {\{0,1\}}^V}{\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\lambda }^{n_1(\sigma )}$ is the partition function.

There are two examples of particular interest: the hardcore model and the Ising model. When $\beta =0$ and $\gamma =1$ then the only configurations with non-zero weight are independent sets of $G$ and the weight of an independent set $\sigma$ is $w(\sigma )={\lambda }^{|\sigma |}$; this example is known as the hardcore model where the parameter $\lambda$ corresponds to the fugacity.

In the case $\beta =\gamma$ then the important quantity is the total number of monochromatic edges $m(\sigma )=m_0(\sigma )+m_1(\sigma )$ and the weight of a configuration $\sigma$ is $w(\sigma )={\beta }^{m(\sigma )}{\lambda }^{n_1(\sigma )}$; this is the classical Ising model where the parameter $\beta$ corresponds to the inverse temperature and $\lambda$ is the external field ($\lambda =1$ means no external field). Note, when $\beta >1$ then the model is ferromagnetic as neighboring vertices prefer to have the same spin, and is the antiferromagnetic Ising model. In the general $2$-spin system, the model is ferromagnetic when $\beta \gamma >1$ and antiferromagnetic when . (When $\beta \gamma =1$ we get a trivial product distribution.)

The fundamental algorithmic tasks are to sample from the Gibbs distribution and to estimate the partition function. For the approximate sampling problem we are given a graph $G$ and an $ϵ>0$ and our goal is to generate a sample from a distribution $\pi$ which is within total variation distance $\le ϵ$ of the Gibbs distribution ${\mu }_G$ in time $\mathrm{poly}(n,\log (1/ϵ))$. An efficient approximate sampling algorithm implies an $\mathrm{𝖥𝖯𝖱𝖠𝖲}$ (fully-polynomial randomized approximation scheme) for the approximate counting problem [JVV86, ŠVV09]. Recall, given an $n$-vertex graph $G$, and $ϵ,\delta >0$, an $\mathrm{𝖥𝖯𝖱𝖠𝖲}$ outputs a $(1\pm ϵ)$-approximation of $Z_G$ with probability $\ge 1-\delta$ in time $\mathrm{poly}(n,1/ϵ,\log (1/\delta ))$, whereas an $\mathrm{𝖥𝖯𝖳𝖠𝖲}$ is the deterministic analog (i.e., $\delta =0$).

A standard approach to the approximate sampling problem is the Markov Chain Monte Carlo (MCMC) method; in fact there is a simple Markov chain known as the Glauber dynamics. The Glauber dynamics works as follows: from a configuration $X_t$ at time $t$, choose a random vertex $v$, we then set $X_{t+1}(w)=X_t(w)$ for all $w\ne v$, and finally we choose $X_{t+1}(v)$ from the conditional distribution of $\mu ({\sigma }_v|{\sigma }_w=X_{t+1}(w)\text{for all}w\ne v)$. For the case of the hardcore model, then $X_{t+1}(v)$ is set to occupied (i.e., spin $1$) with probability $\lambda /(1+\lambda )$ if no neighbors are currently occupied, and otherwise it is set to unoccupied.

It is straightforward to verify that the Glauber dynamics is ergodic with the Gibbs distribution as the unique stationary distribution. The mixing time is the minimum number of steps to guarantee, from the worst initial state $X_0$, that the distribution of $X_t$ is within total variation distance $\le 1/4$ of the Gibbs distribution. The goal is to prove that the mixing time is polynomial in $n$, in which case the chain is said to be rapidly mixing.

For the case of the ferromagnetic Ising model (with or without an external field), a classical result of Jerrum and Sinclair [JS93] gives an $\mathrm{𝖥𝖯𝖱𝖠𝖲}$ for all graphs via the MCMC method. This is the only case with an efficient algorithm for general graphs. For antiferromagnetic 2-spin models the picture is closely tied to statistical physics phase transitions on the regular tree.

The uniqueness/non-uniqueness phase transition is nicely illustrated for the case of the hardcore model. Consider the infinite $\mathrm{\Delta }$-regular tree $T$ rooted at $r$, and let $T_h$ denote the tree truncated at the first $h$ levels. This phase transition captures whether the configuration at the leaves of $T_h$ “influences” the root, in the limit $h\to \mathrm{\infty }$. For the hardcore model we can consider even height trees (corresponding to the all even boundary condition) versus odd height trees. Let $p_h$ denote the marginal probability that the root is occupied in the Gibbs distribution ${\mu }_{T_h}$. Let $p_{\mathrm{𝖾𝗏𝖾𝗇}}={lim}_{h\to \mathrm{\infty }}{p}_{2h}$ and $p_{\mathrm{𝗈𝖽𝖽}}={lim}_{h\to \mathrm{\infty }}{p}_{2h+1}$. We say that tree uniqueness holds if $p_{\mathrm{𝖾𝗏𝖾𝗇}}=p_{\mathrm{𝗈𝖽𝖽}}$ and tree non-uniqueness holds if they are not equal. For all $\mathrm{\Delta }\ge 3$ there exists a critical fugacity ${\lambda }_c(\mathrm{\Delta })={(\mathrm{\Delta }-1)}^{\mathrm{\Delta }-1}/{(\mathrm{\Delta }-2)}^{\mathrm{\Delta }})$ [Kel85], where tree uniqueness holds iff $\lambda \le {\lambda }_c(\mathrm{\Delta })$.

The remarkable connection is that an algorithmic phase transition for general graphs of maximum degree $\mathrm{\Delta }$ occurs at this same tree critical point. For all constant $\mathrm{\Delta }$, all $\delta >0$, all , all graphs of maximum degree $\mathrm{\Delta }$, [Wei06] presented an $\mathrm{𝖥𝖯𝖳𝖠𝖲}$ for approximating the partition function. On the other side, for all $\delta >0$, all $\lambda >(1+\delta ){\lambda }_c(\mathrm{\Delta })$, [Sly10, SS14, GŠV16] proved that, unless $\mathrm{𝖭𝖯}=\mathrm{𝖱𝖯}$, there is no $\mathrm{𝖥𝖯𝖱𝖠𝖲}$ for estimating the partition function.

One important caveat is that the running time of Weitz’s algorithm is ${(n/ϵ)}^{C\log \mathrm{\Delta }}$ where the approximation factor is $(1\pm ϵ)$ and the constant $C$ depends polynomially on the gap $\delta$ (recall, ). Weitz’s correlation decay algorithm was extended to the antiferromagnetic Ising model in the tree uniqueness region by Sinclair et al. [SST14], and to all antiferromagnetic 2-spin systems in the corresponding tree uniqueness region (as we detail below) by Li, Lu, and Yin [LLY13].

An intriguing new algorithmic approach was presented by Barvinok [Bar16] and refined by Patel and Regts [PR17], utilizing the absence of zeros of the partition function in the complex plane to efficiently approximate a suitable transformation of the logarithm of the partition function using Taylor approximation. This polynomial interpolation approach was shown to be efficient in the same tree uniqueness region as for Weitz’s result by Peters and Regts [PR19], although the exponent in the running time depends exponentially on $\mathrm{\Delta }$.

It was long conjectured that the simple Glauber dynamics is rapidly mixing in the tree uniqueness region. This was recently proved by Anari, Liu, and Oveis Gharan [ALO20]; they proved, for all $\delta >0$, the mixing time is $n^{O(\exp (1/\delta ))}$ whenever . We improve this result. First, we improve the mixing time from $n^{O(\exp (1/\delta ))}$ to $n^{O(1/\delta )}$ as detailed in the following theorem.

Theorem 1 (Hardcore model).

Let $\mathrm{\Delta }\ge 3$ be an integer and $\delta \in (0,1)$. For every $n$-vertex graph $G$ of maximum degree $\mathrm{\Delta }$ and every , the mixing time of the Glauber dynamics for the hardcore model on $G$ with fugacity $\lambda$ is $O(n^{2+32/\delta })$.

This bound is optimal barring further improvements in the local-to-global arguments from [AL20]. Our improved result follows from a simpler, cleaner proof approach which enables us to extend our result to a wide variety of 2-spin models, matching the key results for the correlation decay algorithm with vastly improved running times.

Our proof approach unifies the three major algorithmic tools for approximate counting: correlation decay, polynomial interpolation, and MCMC. Most known results for both correlation decay and polynomial interpolation approach are proved by showing contraction of a suitably defined potential function on the so-called tree recursions; the tree recursions arise as a result of Weitz’s self-avoiding walk tree that we will describe in more detail later in this paper. A recent work of Shao and Sun [SS20] unifies these two approaches by showing that the contraction which is normally used to prove efficiency of the correlation decay algorithm, also implies (under some additional analytic conditions) that the polynomial interpolation approach is efficient.

Here we prove that this same contraction of a potential function also implies rapid mixing of the Glauber dynamics, with our improved running time that is independent of $\mathrm{\Delta }$; see 4 and 5 for a detailed statement. Our proof utilizes several new tools concerning Weitz’s self-avoiding walk tree, which are detailed in Section 3. In particular, we show that the partition function of a graph $G$ divides the partition function of Weitz’s self-avoiding walk tree; see 8. This result is potentially of independent interest for establishing absence of zeros for the partition function with complex parameters, as it enables one to consider the self-avoiding walk tree. This result also yields a new, useful equivalence for bounding the influence in a graph in terms of the self-avoiding tree, which strengthens the previously known connection by Weitz [Wei06]; see 8 for details.

As an easy consequence we obtain rapid mixing for the Glauber dynamics for the antiferromagnetic Ising model in the tree uniqueness region. In terms of the edge activity, the two critical points for the Ising model on the $\mathrm{\Delta }$-regular tree are at ${\beta }_c(\mathrm{\Delta })=\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and ${\overline{\beta }}_c(\mathrm{\Delta })=\frac{1}{{\beta }_c(\mathrm{\Delta })}=\frac{\mathrm{\Delta }}{\mathrm{\Delta }-2}$; the first lies in the antiferromagnetic regime, while the second lies in the ferromagnetic regime. If , then uniqueness holds for all external field $\lambda$ on the $\mathrm{\Delta }$-regular tree.

As mentioned earlier, for the ferromagnetic Ising model, an $\mathrm{𝖥𝖯𝖱𝖠𝖲}$ was known for general graphs [JS93]. Furthermore, Mossel and Sly [MS13] proved $O(n\log n)$ mixing time of the Glauber dynamics for the ferromagnetic Ising model when . However, rapid mixing for the antiferromagnetic Ising model in the tree uniqueness region was not known.

We provide the following mixing result for the case $\beta >{\beta }_c(\mathrm{\Delta })$. Note, when $\beta \le {\beta }_c$ there is an additional uniqueness region for certain values of the external field $\lambda$; this region is covered by 3.

Theorem 2 (Antiferromagnetic Ising Model).

Let $\mathrm{\Delta }\ge 3$ be an integer and $\delta \in (0,1)$. Assume that $1>\beta \ge {\beta }_c(\mathrm{\Delta })+\delta (1-{\beta }_c(\mathrm{\Delta }))$ and $\lambda >0$. Then for every $n$-vertex graph $G$ of maximum degree $\mathrm{\Delta }$, the mixing time of the Glauber dynamics for the Ising model on $G$ with edge weight $\beta$ and external field $\lambda$ is $O(n^{2+1.5/\delta })$.

Our results for the hardcore and Ising models fit within a larger framework of general antiferromagnetic 2-spin systems. Recall that the antiferromagnetic case is when .

For general 2-spin systems the appropriate tree phase transition is more complicated as there are models where the tree uniqueness threshold is not monotone in $\mathrm{\Delta }$. Hence the appropriate notion is “up-to-$\mathrm{\Delta }$ uniqueness” as considered by [LLY13]. Roughly speaking, we say uniqueness with gap $\delta \in (0,1)$ holds on the $d$-regular tree if for every integer $\mathrm{\ell }\ge 1$, all vertices at distance $\mathrm{\ell }$ from the root have total “influence” $\lesssim {(1-\delta )}^{\mathrm{\ell }}$ on the marginal of the root. We say up-to-$\mathrm{\Delta }$ uniqueness with gap $\delta$ holds if uniqueness with gap $\delta$ holds on the $d$-regular tree for all $1\le d\le \mathrm{\Delta }$; see Section 2 for the precise definition.

Both 1 and 2 are corollaries of the following general rapid mixing result which holds for general antiferromagnetic $2$-spin systems in the entire tree uniqueness region.

Theorem 3 (General antiferromagnetic $2$-spin system).

Let $\mathrm{\Delta }\ge 3$ be an integer and $\delta \in (0,1)$. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. Assume that the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta$. Then for every $n$-vertex graph $G$ of maximum degree $\mathrm{\Delta }$, the mixing time of the Glauber dynamics for the antiferromagnetic $2$-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $O(n^{2+72/\delta })$.

We also match existing correlation decay results [GL18, SS20] for ferromagnetic 2-spin models; see Section 8 for results, and Appendix F for proofs.

1.1 Mixing by the potential method

The tree recursion is very useful in the study of approximating counting. Consider a tree rooted at $r$. Suppose that $r$ has $d$ children, denoted by $v_1,\mathrm{\dots },v_d$. For $1\le i\le {\mathrm{\Delta }}_i$ we define $T_{v_i}$ to be the subtree of $T$ rooted at $v_i$ that contains all descendant of $v_i$. Let $$Rr=μT​(σr​ = ​1)/μT​(σr​ = ​0) denote the marginal ratio of the root, and $$Rvi=μTvi​(σvi​ = ​1)/μTvi​(σvi​ = ​0) for each subtree. The tree recursion is a formula that computes $R_r$ given $R_{v_1},\mathrm{\dots },R_{v_d}$, due to the independence of $T_{v_i}$’s. More specifically, we can write $R_r=F_d(R_{v_1},\mathrm{\dots },R_{v_d})$ where $F_d:{[0,+\mathrm{\infty }]}^d\to [0,+\mathrm{\infty }]$ is a multivariate function such that for $(x_1,\mathrm{\dots },x_d)\in {[0,+\mathrm{\infty }]}^d$,

$$F_d(x_1,\mathrm{\dots },x_d)=\lambda \prod _{i=1}^d\frac{\beta x_i+1}{x_i+\gamma }.$$

In this paper, however, we pay particular interest in the log of marginal ratios. The reason is that we will carefully study the pairwise influence matrix ${ℐ}_G$ of the Gibbs distribution ${\mu }_G$, introduced in [ALO20] and defined as for every $r,v\in V$

$${ℐ}_G(r\to v)={\mu }_G({\sigma }_v=1\mid {\sigma }_r=1)-{\mu }_G({\sigma }_v=1\mid {\sigma }_r=0).$$

In [ALO20], the authors show that if the maximum eigenvalue of ${ℐ}_G$ is bounded appropriately, then the Glauber dynamics is rapid mixing. One crucial observation we make in this paper is that the influence $$ℐG​(r​ → ​v) of $r$ on $v$ can be viewed as the derivative of $\log R_r$ with respect to the log external field at $v$ (see 12). Thus, it is more convenient for us to work with the log ratios. To this end, we rewrite the tree recursion as $\log R_v=H_d(\log R_{v_1},\mathrm{\dots },\log R_{v_d})$ where $H_d:{[-\mathrm{\infty },+\mathrm{\infty }]}^d\to [-\mathrm{\infty },+\mathrm{\infty }]$ is a function such that for $(y_1,\mathrm{\dots },y_d)\in {[-\mathrm{\infty },+\mathrm{\infty }]}^d$,

$$H_d(y_1,\mathrm{\dots },y_d)=\log \lambda +\sum _{i=1}^d\log \left(\frac{\beta e^{y_i}+1}{e^{y_i}+\gamma }\right).$$

Observe that $H=\log \circ F\circ \exp$. Moreover, we define

$$h(y)=-\frac{(1-\beta \gamma )e^y}{(\beta e^y+1)(e^y+\gamma )}$$

for $y\in [-\mathrm{\infty },+\mathrm{\infty }]$, so that $\frac{\partial }{\partial y_i}{H}_d(y_1,\mathrm{\dots },y_d)=h(y_i)$ for each $i$.

To prove our main results, we use the potential method, which has been widely used to establish the decay of correlation. By choosing a suitable potential function for the log ratios, we show that the total influence from a given vertex decays exponentially with the distance, and thus establish rapid mixing of the Glauber dynamics. Let us first specify our requirements on the potential. For every integer $d\ge 0$, we define a bounded interval $J_d$ which contains all log ratios at a vertex of degree $d$. More specifically, we let $J_d=[\log (\lambda {\beta }^d),\log (\lambda /{\gamma }^d)]$ when , and $J_d=[\log (\lambda /{\gamma }^d),\log (\lambda {\beta }^d)]$ when $\beta \gamma >1$. Furthermore, define $J={\bigcup }_{d=0}^{\mathrm{\Delta }-1}{J}_d$ to be the interval containing all log ratios with degree less than $\mathrm{\Delta }$.

Definition 4 ($(\alpha ,c)$-Potential function).

Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$ and $\lambda >0$. Let $\mathrm{\Psi }:[-\mathrm{\infty },+\mathrm{\infty }]\to (-\mathrm{\infty },+\mathrm{\infty })$ be a differentiable and increasing function with image $S=\mathrm{\Psi }[-\mathrm{\infty },+\mathrm{\infty }]$ and derivative $\psi ={\mathrm{\Psi }}^{\prime }$. For any $\alpha \in (0,1)$ and $c>0$, we say $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential function with respect to $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$ if it satisfies the following conditions:

1. (Contraction) For every integer $d$ such that and every $({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\in S^d$, we have

   $${\Vert \nabla H_d^{\mathrm{\Psi }}({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\Vert }_1=\sum _{i=1}^d\frac{\psi (y)}{\psi (y_i)}\cdot |h(y_i)|\le 1-\alpha$$

   where $H_d^{\mathrm{\Psi }}=\mathrm{\Psi }\circ H_d\circ {\mathrm{\Psi }}^{-1}$, $y_i={\mathrm{\Psi }}^{-1}({\tilde{y}}_i)$ for $1\le i\le d$, and $y=H_d(y_1,\mathrm{\dots },y_d)$.

2. (Boundedness) For every $y_1,y_2\in J$, we have

   $$\frac{\psi (y_2)}{\psi (y_1)}\cdot \left|h(y_1)\right|\le \frac{c}{\mathrm{\Delta }}.$$

In the definition of $(\alpha ,c)$-potential, one should think of $y$ as the log marginal ratio at a vertex and the potential function is of $\log R$. The following theorem establishes rapid mixing of the Glauber dynamics given an $(\alpha ,c)$-potential function.

Theorem 5.

Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$ and $\lambda >0$. Suppose that there is an $(\alpha ,c)$-potential with respect to $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$ for some $\alpha \in (0,1)$ and $c>0$. Then for every $n$-vertex graph $G$ of maximum degree $\mathrm{\Delta }$, the mixing time of the Glauber dynamics for the $2$-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $O(n^{2+c/\alpha })$.

We outline our proofs in Section 3. Note that in both 4 and 5, the constant $c$ is allowed to depend on the maximum degree $\mathrm{\Delta }$ and parameters $(\beta ,\gamma ,\lambda )$ in general. For example, a straightforward black-box application of the potential in [LLY13] would give $c=\mathrm{\Theta }(\mathrm{\Delta })$ for the Boundedness condition, resulting in $n^{\mathrm{\Theta }(\mathrm{\Delta })}$ mixing. However, this is undesirable for graphs with potentially unbounded degrees. One of our contributions is that we show the Boundedness condition holds for a universal constant $c$ independent of $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$. Thus, our mixing time is $O(n^{2+c/\delta })$ with no parameters in the exponent except for $1/\delta$.

In Section 7, we give a slightly more general definition of $(\alpha ,c)$-potentials, which relaxes the Boundedness condition, and is necessary for our analysis of antiferromagnetic $2$-spin systems with . 5 still holds for this larger class of potentials.

We remark that in all previous works of the potential method, results and proofs are always presented in terms of $F_d$, the tree recursion of $R$, and $\mathrm{\Phi }$, a potential function of $R$. In fact, our results can also be translated into the language of $(F_d,\mathrm{\Phi })$. To see this, since $H_d=\log \circ F_d\circ \exp$, it is straightforward to check that $H_d^{\mathrm{\Psi }}=\mathrm{\Psi }\circ H_d\circ {\mathrm{\Psi }}^{-1}=\mathrm{\Phi }\circ F_d\circ {\mathrm{\Phi }}^{-1}=F_d^{\mathrm{\Phi }}$ if we pick $\mathrm{\Phi }=\mathrm{\Psi }\circ \log$, and thereby $\nabla H_d^{\mathrm{\Psi }}=\nabla F_d^{\mathrm{\Phi }}$. This implies that the Contraction condition in 4 holds for $(H_d,\mathrm{\Psi })$ if and only if the corresponding contraction condition holds for $(F_d,\mathrm{\Phi })$. The Boundedness condition can also be stated equivalently for $(F_d,\mathrm{\Phi })$. Nevertheless, in this paper we choose to work with $(H_d,\mathrm{\Psi })$ for the following two reasons. First, as mentioned earlier, the fact that $$ℐG​(r​ → ​v) is a derivative of $\log R_r$ makes it natural to consider the tree recursion for the log ratios. Indeed, it is easier and cleaner to present our results and proofs using $(H_d,\mathrm{\Psi })$ directly rather than switching to $(F_d,\mathrm{\Phi })$. Second, the potential function $\mathrm{\Psi }$ we will use is obtained from the exact potential $\mathrm{\Phi }$ in [LLY13], by the transformation $\mathrm{\Psi }=\mathrm{\Phi }\circ \exp$. It is intriguing to notice that the derivative of this potential is simply $\psi =\sqrt{|h|}$. Then the Contraction condition has a nice form: ${\sum }_{i=1}^d\sqrt{h(y)h(y_i)}\le 1-\alpha$; and the Boundedness condition only involves an upper bound on $h(y)$. This seems to shed some light on the mysterious potential function $\mathrm{\Phi }$ from [LLY13], and also indicates that $H_d$ is a meaningful variant of the tree recursion to consider. To add one more evidence, for a lot of cases (e.g., ) where the potential $\mathrm{\Phi }=\log$ is picked, that just means we can pick $\mathrm{\Psi }$ to be the identity function and $H_d$ itself is contracting without any nontrivial potential.

Revision in July 2021.

After the publication of this paper in FOCS 2020, a small error was found in [LLY13] regarding descriptions of the uniqueness region for antiferromagnetic 2-spin systems. The error was fixed in the latest version of [LLY13]. In this revision, we update corresponding results and proofs in Section 7 and Appendix E that rely on the changes in [LLY13]; in particular, 36 is adjusted in accordance with the current description of uniqueness regions. We remark that these changes are purely technical and do not affect the validity of our main results like 5.

Acknowledgments.

We would like to thank Shayan Oveis Gharan and Nima Anari for stimulating discussions. We also thank the anonymous referees for helpful comments and suggestions. We are grateful to Yitong Yin for communicating with us about the latest update of [LLY13] and for providing helpful instructions on modifying statements and proofs of results in Appendix E, particularly 36.

2 Preliminaries

Mixing time and spectral gap

Let $P$ be the transition matrix of an ergodic (i.e., irreducible and aperiodic) Markov chain on a finite state space $\mathrm{\Omega }$ with stationary distribution $\mu$. Let $P^t(x_0,\cdot )$ denote the distribution of the chain after $t$ steps starting from $x_0\in \mathrm{\Omega }$. The mixing time of $P$ is defined as

$$T_{\mathrm{mix}}(P)=\underset{x_0\in \mathrm{\Omega }}{\max }\min \left\{t\ge 0:{\Vert P^t(x_0,\cdot )-\mu (\cdot )\Vert }_{\mathrm{TV}}\le \frac{1}{4}\right\}.$$

We say $P$ is reversible if $\mu (x)P(x,y)=\mu (y)P(y,x)$ for all $x,y\in \mathrm{\Omega }$. If $P$ is reversible, then $P$ has only real eigenvalues which can be denoted by $1={\lambda }_1\ge \mathrm{\cdots }\ge {\lambda }_{|\mathrm{\Omega }|}\ge -1$. The spectral gap of $P$ is defined to be $1-{\lambda }_2$ and the absolute spectral gap of $P$ is defined as ${\lambda }^{\ast }(P)=1-\max \{|{\lambda }_2|,|{\lambda }_{|\mathrm{\Omega }|}|\}$. If $P$ is also positive semidefinite with respect to the inner product ${⟨\cdot ,\cdot ⟩}_{\mu }$, then all eigenvalues of $P$ are nonnegative and thus ${\lambda }^{\ast }(P)=1-{\lambda }_2$. Finally, the mixing time and the absolute spectral gap are related by

$$T_{\mathrm{mix}}(P)\le \frac{1}{{\lambda }^{\ast }(P)}\log \left(\frac{4}{{\min }_{x\in \mathrm{\Omega }}\mu (x)}\right).$$ (1)

Uniqueness

Let $\mathrm{\Delta }\ge 3$ be an integer or $\mathrm{\Delta }=\mathrm{\infty }$. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. For , define

$$f_d(R)=\lambda {\left(\frac{\beta R+1}{R+\gamma }\right)}^d$$

and denote the unique fixed point of $f_d$ by $R_d^{\ast }$. For $\delta \in (0,1)$, we say the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta$ if for all .

Ratio and influence

Consider the $2$-spin system on a graph $G=(V,E)$. Let $\mathrm{\Lambda }\subseteq V$ and ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$. For all $v\in V\backslash \mathrm{\Lambda }$, we define the marginal ratio at $v$ to be

$$R_G^{{\sigma }_{\mathrm{\Lambda }}}(v)=\frac{{\mu }_G({\sigma }_v=1\mid {\sigma }_{\mathrm{\Lambda }})}{{\mu }_G({\sigma }_v=0\mid {\sigma }_{\mathrm{\Lambda }})}.$$

For all $u,v\in V\backslash \mathrm{\Lambda }$, we define the (pairwise) influence of $u$ on $v$ by

$${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(u\to v)={\mu }_G({\sigma }_v=1\mid {\sigma }_u=1,{\sigma }_{\mathrm{\Lambda }})-{\mu }_G({\sigma }_v=1\mid {\sigma }_u=0,{\sigma }_{\mathrm{\Lambda }}).$$

Write ${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}$ for the (pairwise) influence matrix whose entries are given by $$ℐGσΛ​(u​ → ​v).

Weitz’s self-avoiding walk tree

Let $G=(V,E)$ be a connected graph and $r\in V$ be a vertex of $G$. The self-avoiding walk (SAW) tree is defined as follows. Suppose that there is a total ordering of the vertex set $V$. A self-avoiding walk from $r$ is a path $r=v_0-v_1-\mathrm{\cdots }-v_{\mathrm{\ell }}$ such that $v_i\ne v_j$ for all . The SAW tree $T_{\text{saw}}(G,r)$ is a tree rooted at $r$, consisting of all self-avoiding walks $r=v_0-v_1-\mathrm{\cdots }-v_{\mathrm{\ell }}$ with $\mathrm{deg}(v_{\mathrm{\ell }})=1$, and those appended with one more vertex that closes the cycle (i.e., $r=v_0-v_1-\mathrm{\cdots }-v_{\mathrm{\ell }}-v_i$ for some $0\le i\le \mathrm{\ell }-2$ such that $\{v_{\mathrm{\ell }},v_i\}\in E$). Note that a vertex of $G$ might have many copies in the SAW tree, and the degrees of vertices are preserved except for leaves. See Fig. 1 for an example.

We can define a $2$-spin system on $T_{\text{saw}}(G,r)$ with the same parameters $(\beta ,\gamma ,\lambda )$, in which some of the leaves are fixed to a particular spin. More specifically, for a self-avoiding walk $r=v_0-v_1-\mathrm{\cdots }-v_{\mathrm{\ell }}$ appended with $v_i$, we fix $v_i$ to be spin $1$ if with respect to the total ordering on $V$, and spin $0$ if $v_{i+1}>v_{\mathrm{\ell }}$. For each $v\in V$ we denote the set of all free (unfixed) copies of $v$ in $T_{\text{saw}}(G,r)$ by ${𝒞}_v$. For $\mathrm{\Lambda }\subseteq V$ and a partial configuration ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, we define the SAW tree with conditioning ${\sigma }_{\mathrm{\Lambda }}$ by assigning the spin ${\sigma }_v$ to every copy $\widehat{v}$ of $v$ from ${𝒞}_v$ and removing all descendants of $\widehat{v}$, for each $v\in \mathrm{\Lambda }$. Note that in general, different copies of $v$ from ${𝒞}_v$ can receive different spin assignments. Finally, in the case that every vertex $v$ has a distinct field ${\lambda }_v$, all copies of $v$ from ${𝒞}_v$ will have the same field ${\lambda }_v$ in the SAW tree.

Figure 1: A graph $G$ and the self-avoiding walk tree $T_{\text{saw}}(G,r)$ rooted at $r$. Vertices with the same label in $T_{\text{saw}}(G,r)$ are copies of the same vertex from $G$. ($\mathrm{●}$/$\mathrm{○}$: fixed to spin $1$/$0$.)

3 Proof outline for main results

Step 1 ([ALO20]): Spectral Independence implies rapid mixing.

Our proof builds on [ALO20] who showed that the Glauber dynamics for sampling from the hardcore distribution on graphs of maximum degree at most $\mathrm{\Delta }$ mixes in $O(n^{\exp (O(1/\delta ))})$ steps whenever $\lambda \le (1-\delta ){\lambda }_c(\mathrm{\Delta })$. One of the key ingredients of their proof is a notion they call spectral independence. [ALO20] shows that the spectral independence property implies rapid mixing. Note that the diagonal entries of ${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}$ are $1$, as opposed to $0$ in the original definition in [ALO20].

Definition 6 (Spectral Independence [ALO20]).

We say the Gibbs distribution ${\mu }_G$ on an $n$-vertex graph $G$ is $({\eta }_0,\mathrm{\dots },{\eta }_{n-2})$-spectrally independent, if for every $0\le k\le n-2$, $\mathrm{\Lambda }\subseteq V$ of size $k$ and ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, one has ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})-1\le {\eta }_k$.

Theorem 7 ([ALO20]).

If $\mu$ is an $({\eta }_0,\mathrm{\dots },{\eta }_{n-2})$-spectrally independent distribution, then the Glauber dynamics for sampling from $\mu$ has spectral gap at least

$$\frac{1}{n}\prod _{i=0}^{n-2}\left(1-\frac{{\eta }_i}{n-i-1}\right).$$

Our primary goal now is to bound the maximum eigenvalue of ${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}$.

Step 2: Self-avoiding walk trees preserve influences.

From standard linear algebra, we know that the maximum eigenvalue of ${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}$ is upper bounded by both the $1$-norm $$‖ℐGσΛ‖1=maxr∈V​∑v∈V|ℐGσΛ​(v​ → ​r)|, which corresponds to total influences on a vertex $r$, and the infinity-norm $$‖ℐGσΛ‖∞=maxr∈V​∑v∈V|ℐGσΛ​(r​ → ​v)|, corresponding to total influences of $r$. In [ALO20] the authors use ${\Vert {ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}\Vert }_1$ as an upper bound on ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})$. Roughly speaking, they show that the sum of absolute influences on a fixed vertex $r$, is upper bounded by the maximum absolute influences on $r$ in the self-avoiding walk tree rooted at $r$, over all boundary conditions. Here in this paper, we will use ${\Vert {ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}\Vert }_{\mathrm{\infty }}$ to upper bound ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})$ instead. In fact, much more is true if we look at the influences from $r$ in the self-avoiding tree. We show that for every vertex $v\in V$, the influence $$ℐGσΛ​(r​ → ​v) in $G$ is preserved in the self-avoiding walk tree $T=T_{\text{saw}}(G,r)$ rooted at $r$, in the form of sum of influences $$ℐTσΛ​(r​ → ​v^) over all copies $\widehat{v}$ of $v$.

The way we establish this fact is by viewing the partition function as a polynomial in $\lambda$. In fact, it will be useful to consider the more general case with an arbitrary external field ${\lambda }_v$ for every $v\in V$. Let $𝝀=\{{\lambda }_v:v\in V\}$ denote the fields. For $\mathrm{\Lambda }\subseteq V$ and ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, the weight of $\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}$ conditional on ${\sigma }_{\mathrm{\Lambda }}$ is defined to be $w_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})={\beta }^{m_1(\sigma \mid {\sigma }_{\mathrm{\Lambda }})}{\gamma }^{m_0(\sigma \mid {\sigma }_{\mathrm{\Lambda }})}{\prod }_{v\in V\backslash \mathrm{\Lambda }}{\lambda }_v^{{\sigma }_v}$ where $m_i(\cdot \mid {\sigma }_{\mathrm{\Lambda }})$ is the number of $i$-$i$ edges with at least one endpoint in $V\backslash \mathrm{\Lambda }$ for $i=0,1$. Furthermore, $Z_G^{{\sigma }_{\mathrm{\Lambda }}}={\sum }_{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}{w}_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})$ is the partition function conditioned on ${\sigma }_{\mathrm{\Lambda }}$. We shall view $\beta$ and $\gamma$ as some fixed constants and think of $𝝀$ as $n=|V|$ variables. In this sense, we regard the weights $w_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})$ as monomials in $𝝀$ and the partition function $Z_G^{{\sigma }_{\mathrm{\Lambda }}}$ as a polynomial in $𝝀$. Moreover, the marginal ratios $R_G^{{\sigma }_{\mathrm{\Lambda }}}(v)$ and the influences $$ℐGσΛ​(r​ → ​v) for $r,v\in V$ are all functions in $𝝀$. Our main result is that the partition function of $G$ divides that of $T_{\text{saw}}(G,r)$ for each $r\in V$. From that, we show that the SAW tree preserves influences of the root, as well as re-establishing Weitz’s celebrated result [Wei06], see 13.

Lemma 8.

Let $G=(V,E)$ be a connected graph, $r\in V$ be a vertex and $\mathrm{\Lambda }\subseteq V\backslash \{r\}$ such that $G\backslash \mathrm{\Lambda }$ is connected. Let $T=T_{\text{saw}}(G,r)$ be the self-avoiding walk tree of $G$ rooted at $r$. Then for every ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, $Z_G^{{\sigma }_{\mathrm{\Lambda }}}$ divides $Z_T^{{\sigma }_{\mathrm{\Lambda }}}$. More precisely, there exists a polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}(𝛌)$ independent of ${\lambda }_r$ such that

$Z_T^{{\sigma }_{\mathrm{\Lambda }}}=Z_G^{{\sigma }_{\mathrm{\Lambda }}}\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}.$ (2)

As a corollary, for each vertex $v\in V$,

${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)={\displaystyle \sum _{\widehat{v}\in {𝒞}_v}}{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to \widehat{v}),$ (3)

where ${𝒞}_v$ is the set of all free (unfixed) copies of $v$ in $T$.

Remark 1.

We emphasize that for the purposes of bounding the total influence of a vertex in $G$, only Eq. 3 of 8 is needed, which can be proved in a purely combinatorial fashion. However, we believe the divisibility property Eq. 2 of the multivariate partition function of $G$ and its self-avoiding walk tree may be of independent interest.

We note that a univariate version of the divisibility statement Eq. 2 has already appeared in [Ben18] for the hardcore model and [LSS19] for the zero-field Ising model in the study of complex roots of the partition function. From 8, we can get ${\sum }_{v\in V}|$ℐGσΛ​(r​ → ​v)|≤∑v∈VT|ℐTσΛ​(r​ → ​v)| for any fixed $r$. That means, we only need to upper bound the sum of all influences for trees, in order to get an upper bound on ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})$.

Step 3: Decay of influences given a good potential.

The tree recursion provides us a great tool for computing the (log) ratios of vertices recursively for trees. As we show in 12, the influence $$ℐGσΛ​(r​ → ​v) is in fact a version of derivative of the log marginal ratio at $r$. Thus, the tree recursion can be used naturally to relate these influences. We then apply the potential method, which has been widely used in literature to establish the decay of correlations (strong spatial mixing). The following lemma shows that the sum of absolute influences to distance $k$ has exponential decay with $k$, which can be thought of as the decay of pairwise influences.

Lemma 9.

If there exists an $(\alpha ,c)$-potential function $\mathrm{\Psi }$ with respect to $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$ where $\alpha \in (0,1)$ and $c>0$, then for every $\mathrm{\Lambda }\subseteq V_T\backslash \{r\}$, ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ and all integers $k\ge 1$,

$$\sum _{v\in L_r(k)}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|\le c\cdot {(1-\alpha )}^{k-1}$$

where $L_r(k)$ denote the set of all free vertices at distance $k$ away from $r$.

5 is then proved by combining 7, 8 and 9. We leave its proof to Appendix A.

Step 4: Find a good potential.

As our final step, we need to find an $(\alpha ,c)$-potential function as defined in 4. The potential $\mathrm{\Psi }$ we choose is exactly the one from [LLY13], adapted to the log marginal ratios and the tree recursion $H$ (see Section 6 for more details). We show that if the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$ and either $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ or $\gamma \le 1$, then $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential.

Lemma 10.

Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. Assume that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$. Define the function $\mathrm{\Psi }$ implicitly by

$${\mathrm{\Psi }}^{\prime }(y)=\psi (y)=\sqrt{\frac{(1-\beta \gamma )e^y}{(\beta e^y+1)(e^y+\gamma )}}=\sqrt{\left|h(y)\right|},\mathrm{\Psi }(0)=0.$$ (4)

If $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, then $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential function with $\alpha \ge \delta /2$ and $c\le 1.5$. If $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma \le 1$, then $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential with $\alpha \ge \delta /2$ and $c\le 18$; we can further take $c\le 4$ if $\beta =0$.

We deduce 3 for the case $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ or $\gamma \le 1$ from 5 and 10. The proof of it can be found in Appendix A. The case that $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$ is trickier. As discussed in Section 5 of [LLY13], when $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$, for some $\lambda >0$ the spin system lies in the uniqueness region for arbitrary graphs, even with unbounded degrees (i.e., up-to-$\mathrm{\infty }$ unique). Thus, in this case the total influences of a vertex can be as large as $\mathrm{\Theta }(\mathrm{\Delta }/\delta )$, resulting in $n^{\mathrm{\Theta }(\mathrm{\Delta }/\delta )}$ mixing time. To deal with this, we consider a suitably weighted sum of absolute influences of a fixed vertex, which also upper bounds the maximum eigenvalue of the influence matrix. 4 and 5 are then modified to a slightly stronger version. The statements and proofs for this case are presented in Section 7 and Appendix D.

The rest of the paper is organized as follows. In Section 4 we prove 8 about properties of the SAW tree. In Section 5 we establish 9 regarding the decay of influences by the potential method. We verify the Contraction condition in Section 6 for our choice of potential. Section 7 is devoted to the case that $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$, where a more general version of 4 and 5 is required; missing proofs can be found in Appendix D. In Appendix E we verify the Boundedness condition and its generalization for our potential in all cases. We consider ferromagnetic spin systems in Section 8 and the proofs are left to Appendix F. We prove all of our main results in Appendix A.

4 Preservation of influences for self-avoiding walk trees

In this section we show that the self-avoiding walk (SAW) tree, introduced in [Wei06] (see also [SS05]), maintains all the influence of the root, and thus establishes 8. To do this, we show that the partition function of $G$, viewed as a polynomial of the external fields $𝝀$, divides that of the SAW tree. From there we prove that the influence of the root vertex $r$ on another vertex $v$ in $G$, is exactly equal to that on all copies of $v$ in the SAW tree. Using our proof approach, we show that the marginal of the root is maintained in the SAW tree, re-establishing Weitz’s celebrated result [Wei06], and also all pairwise covariances concerned with $v$ are preserved.

Theorem 11.

Let $G=(V,E)$ be a connected graph, $r\in V$ be a vertex and $\mathrm{\Lambda }\subseteq V\backslash \{r\}$ such that $G\backslash \mathrm{\Lambda }$ is connected. Let $T=T_{\text{saw}}(G,r)$ be the self-avoiding walk tree of $G$ rooted at $r$. Then for every ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, $Z_G^{{\sigma }_{\mathrm{\Lambda }}}$ divides $Z_T^{{\sigma }_{\mathrm{\Lambda }}}$. More precisely, there exists a polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}(𝛌)$ such that

$$Z_T^{{\sigma }_{\mathrm{\Lambda }}}=Z_G^{{\sigma }_{\mathrm{\Lambda }}}\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}.$$

Moreover, the polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}$ is independent of ${\lambda }_r$.

Remark 2.

The proof of 11 can be adapted to give a purely combinatorial proof of Eq. 3 in 8. Like in the proof of [Wei06, Theorem 3.1], one can proceed via vertex splitting and telescoping, where instead of telescoping a product of marginal ratios, one instead telescopes a sum of single-vertex influences.

We remark that [Ben18] proved a univariate version of 11 for the hardcore model, and [LSS19] showed a similar result for the zero-field Ising model with a uniform edge weight. Our result holds for all $2$-spin systems and arbitrary fields for each vertex. We can also generalize it to arbitrary edge weights for each edge in a straightforward fashion. It is crucial that the quotient polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}$ is independent of the field ${\lambda }_r$ at the root, from which we can deduce the preservation of marginal and influences of the root immediately.

Before proving 11, we first give a few consequences of it. For all $u,v\in V\backslash \mathrm{\Lambda }$, we define the marginal at $v$ as $$MGσΛ​(v)=μG​(v​ = ​1∣σΛ) (henceforth we write $v=i$ for the event ${\sigma }_v=i$ for convenience), and the covariance of $u$ and $v$ as

$$K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,v)={\mu }_G(u=v=1\mid {\sigma }_{\mathrm{\Lambda }})-{\mu }_G(u=1\mid {\sigma }_{\mathrm{\Lambda }}){\mu }_G(v=1\mid {\sigma }_{\mathrm{\Lambda }}).$$

The following lemma relates the quantities we are interested in with appropriate derivatives of the (log) partition function. Parts 1 and 2 of the lemma are folklore.

Lemma 12.

For every graph $G=(V,E)$, $\mathrm{\Lambda }\subseteq V$ and ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, the following holds:

1. For all $v\in V$,

   $$\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)\log Z_G^{{\sigma }_{\mathrm{\Lambda }}}=M_G^{{\sigma }_{\mathrm{\Lambda }}}(v);$$

2. For all $u,v\in V$,

   $$\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)\left({\lambda }_u\frac{\partial }{\partial {\lambda }_u}\right)\log Z_G^{{\sigma }_{\mathrm{\Lambda }}}=\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)=K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,v);$$

3. For all $u,v\in V$,

   $$\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)\log R_G^{{\sigma }_{\mathrm{\Lambda }}}(u)={ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(u\to v).$$

Proof.

The first two parts are standard. We include the proofs of these two facts in Appendix B for completeness. For Part 3, we deduce from Part 2 that

$\left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\log R_G^{{\sigma }_{\mathrm{\Lambda }}}(u)=\left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\log \left({\displaystyle \frac{M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)}{1-M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)}}\right)={\displaystyle \frac{\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)}{M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)\left(1-M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)\right)}}={\displaystyle \frac{K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,v)}{K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,u)}}.$

It remains to show that

$${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(u\to v)=\frac{K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,v)}{K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,u)},$$

which actually holds for any two binary random variables. To see this, we first compute $$KGσΛ​(u,u)⋅ℐGσΛ​(u​ → ​v) by definition:

$K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,u)\cdot {ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(u\to v)$

$=$ ${\mu }_G(u=1\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0\mid {\sigma }_{\mathrm{\Lambda }})\cdot \left[{\mu }_G(v=1\mid u=1,{\sigma }_{\mathrm{\Lambda }})-{\mu }_G(v=1\mid u=0,{\sigma }_{\mathrm{\Lambda }})\right]$

$=$ ${\mu }_G(u=1,v=1\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0\mid {\sigma }_{\mathrm{\Lambda }})-{\mu }_G(u=1\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0,v=1\mid {\sigma }_{\mathrm{\Lambda }})$

$=$ ${\mu }_G(u=1,v=1\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0,v=0\mid {\sigma }_{\mathrm{\Lambda }})-{\mu }_G(u=1,v=0\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0,v=1\mid {\sigma }_{\mathrm{\Lambda }}).$

Meanwhile, the covariance can be written as

$K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,v)$ $={\mu }_G(u=1,v=1\mid {\sigma }_{\mathrm{\Lambda }})-{\mu }_G(u=1\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(v=1\mid {\sigma }_{\mathrm{\Lambda }})$

$={\mu }_G(u=1,v=1\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0,v=0\mid {\sigma }_{\mathrm{\Lambda }})-{\mu }_G(u=1,v=0\mid {\sigma }_{\mathrm{\Lambda }})\cdot {\mu }_G(u=0,v=1\mid {\sigma }_{\mathrm{\Lambda }}).$

This shows that $$ℐGσΛ​(u​ → ​v)=KGσΛ​(u,v)/KGσΛ​(u,u) and thus establishes Part 3. ∎

We deduce 8 from 11 and the second item of the following lemma. The proof of 11 is presented in Section 4.1.

Lemma 13.

Let $G=(V,E)$ be a connected graph, $r\in V$ be a vertex and $\mathrm{\Lambda }\subseteq V\backslash \{r\}$ such that $G\backslash \mathrm{\Lambda }$ is connected. Let $T=T_{\text{saw}}(G,r)$ be the self-avoiding walk tree of $G$ rooted at $r$. Then for every ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ we have:

1. ([Wei06, Theorem 3.1]) Preservation of marginal of the root $r$:

   $$M_G^{{\sigma }_{\mathrm{\Lambda }}}(r)=M_T^{{\sigma }_{\mathrm{\Lambda }}}(r)\text{and}{R}_G^{{\sigma }_{\mathrm{\Lambda }}}(r)=R_T^{{\sigma }_{\mathrm{\Lambda }}}(r);$$

2. Preservation of covariances and influences of $r$: for every $v\in V$,

   $$K_G^{{\sigma }_{\mathrm{\Lambda }}}(r,v)=\sum _{\widehat{v}\in {𝒞}_v}{K}_T^{{\sigma }_{\mathrm{\Lambda }}}(r,\widehat{v})\text{and}{ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)=\sum _{\widehat{v}\in {𝒞}_v}{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to \widehat{v}).$$

   where ${𝒞}_v$ is the set of all free (unfixed) copies of $v$ in $T$.

Proof.

By 11, there exists a polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}(𝝀)$ such that $Z_T^{{\sigma }_{\mathrm{\Lambda }}}=Z_G^{{\sigma }_{\mathrm{\Lambda }}}\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}$ and $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}$ is independent of ${\lambda }_r$. Then it follows from 12 that

$$M_T^{{\sigma }_{\mathrm{\Lambda }}}(r)=\left({\lambda }_r\frac{\partial }{\partial {\lambda }_r}\right)\log Z_T^{{\sigma }_{\mathrm{\Lambda }}}=\left({\lambda }_r\frac{\partial }{\partial {\lambda }_r}\right)\left(\log Z_G^{{\sigma }_{\mathrm{\Lambda }}}+\log P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}\right)=\left({\lambda }_r\frac{\partial }{\partial {\lambda }_r}\right)\log Z_G^{{\sigma }_{\mathrm{\Lambda }}}=M_G^{{\sigma }_{\mathrm{\Lambda }}}(r),$$

and therefore $R_T^{{\sigma }_{\mathrm{\Lambda }}}(r)=R_G^{{\sigma }_{\mathrm{\Lambda }}}(r)$. For the second item, again from 12 we get

$$K_G^{{\sigma }_{\mathrm{\Lambda }}}(r,v)=\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)M_G^{{\sigma }_{\mathrm{\Lambda }}}(r)=\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)M_T^{{\sigma }_{\mathrm{\Lambda }}}(r).$$

Recall that for the spin system on the SAW tree $T$, every free copy $\widehat{v}$ of $v$ from ${𝒞}_v$ has the same external field ${\lambda }_{\widehat{v}}={\lambda }_v$. Then, by the chain rule of derivatives and 12, we deduce that

$$K_G^{{\sigma }_{\mathrm{\Lambda }}}(r,v)=\sum _{\widehat{v}\in {𝒞}_v}\left({\lambda }_{\widehat{v}}\frac{\partial }{\partial {\lambda }_{\widehat{v}}}\right)M_T^{{\sigma }_{\mathrm{\Lambda }}}(r)\cdot \frac{\partial {\lambda }_{\widehat{v}}}{\partial {\lambda }_v}\cdot \frac{{\lambda }_v}{{\lambda }_{\widehat{v}}}=\sum _{\widehat{v}\in {𝒞}_v}{K}_T^{{\sigma }_{\mathrm{\Lambda }}}(r,\widehat{v}).$$

Finally, we have

$${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)=\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)\log R_G^{{\sigma }_{\mathrm{\Lambda }}}(r)=\left({\lambda }_v\frac{\partial }{\partial {\lambda }_v}\right)\log R_T^{{\sigma }_{\mathrm{\Lambda }}}(r)=\sum _{\widehat{v}\in {𝒞}_v}{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to \widehat{v}),$$

where the last equality follows as above. ∎

4.1 Proof of 11

Before presenting our proof, let us first review the notations and definitions introduced earlier. Denote the set of fields at all vertices by $𝝀=\{{\lambda }_v:v\in V\}$. For $\mathrm{\Lambda }\subseteq V$ and ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$, the weight of $\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}$ conditional on ${\sigma }_{\mathrm{\Lambda }}$ is given by

$$w_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})={\beta }^{m_1(\sigma \mid {\sigma }_{\mathrm{\Lambda }})}{\gamma }^{m_0(\sigma \mid {\sigma }_{\mathrm{\Lambda }})}\prod _{v\in V\backslash \mathrm{\Lambda }}{\lambda }_v^{{\sigma }_v},$$

where for $i=0,1$, $m_i(\cdot \mid {\sigma }_{\mathrm{\Lambda }})$ denotes the number of edges such that both endpoints receive the spin $i$ and at least one of them is in $V\backslash \mathrm{\Lambda }$. The partition function conditional on ${\sigma }_{\mathrm{\Lambda }}$ is defined as $Z_G^{{\sigma }_{\mathrm{\Lambda }}}={\sum }_{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}{w}_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})$. For the SAW tree, we define the conditional weights and partition function in the same way. In particular, recall that when we fix a conditioning ${\sigma }_{\mathrm{\Lambda }}$ on the SAW tree, we also remove all descendants of $\widehat{v}\in {𝒞}_v$ for each $v\in \mathrm{\Lambda }$.

For every $v\in V\backslash \mathrm{\Lambda }$ and $i\in \{0,1\}$, we shall write $v=i$ to represent the set of configurations such that ${\sigma }_v=i$ (i.e., $\{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}:{\sigma }_v=i\}$) and let $$ZGσΛ​(v​ = ​i) be sum of weights of all configurations with $$v​ = ​i. We further extend this notation and write $$ZGσΛ​(U​ = ​σU) for every $U\subseteq V\backslash \mathrm{\Lambda }$ and ${\sigma }_U\in {\{0,1\}}^U$. For the SAW tree we adopt the same notations as well.

Proof of 11.

We will show that there exists a polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}(𝝀)$, independent of ${\lambda }_r$, such that

$$Z_T^{{\sigma }_{\mathrm{\Lambda }}}(r=1)=Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}\text{and}{Z}_T^{{\sigma }_{\mathrm{\Lambda }}}(r=0)=Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=0)\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}.$$ (5)

The high-level proof idea of Eq. 5 is similar to the corresponding result in [Wei06, Theorem 3.1]. Let $m$ be the number of edges with at least one endpoint in $V\backslash \mathrm{\Lambda }$. We use induction on $m$. When $m=0$ the statement is trivial since $T=G$. Assume that Eq. 5 holds for all graphs and all conditioning with less than $m$ edges. Suppose that the root $r$ has $d$ neighbors $v_1,\mathrm{\dots },v_d$. Define $G^{\prime }$ to be the graph obtained by replacing the vertex $r$ with $d$ vertices $r_1,\mathrm{\dots },r_d$ and then connecting $\{r_i,d_i\}$ for $1\le i\le d$.

Consider first the case where $(G\backslash \{r\})\backslash \mathrm{\Lambda }$ is still connected. For each $i$, let $G_i=G^{\prime }-r_i$. Define the $2$-spin system on $G_i$ with the same parameters $(\beta ,\gamma ,𝝀)$, plus an additional conditioning that the vertices $r_1,\mathrm{\dots },r_{i-1}$ are fixed to spin $0$ while $r_{i+1},\mathrm{\dots },r_d$ are fixed to spin $1$; we denote this conditioning by ${\sigma }_{U_i}$ with $U_i=\{v_1,\mathrm{\dots },v_d\}\backslash \{v_i\}$. Then, $T=T_{\text{saw}}(G,r)$ can be generated by the following recursive procedure. Also see Fig. 2 for an illustration.

Figure 2: A recursive construction of the self-avoiding walk (SAW) tree. Here $T_i$ is the SAW tree of $G_i$ rooted at $v_i$ for $i=1,2,3$. ($\mathrm{●}$/$\mathrm{○}$: fixed to spin $1$/$0$.)

Algorithm: $T_{\text{saw}}(G,r)$

1. For each $i$, let $T_i=T_{\text{saw}}(G_i,v_i)$ plus the conditioning ${\sigma }_{U_i}$;

2. Let $T=T_{\text{saw}}(G,r)$ be the union of $r$ and $T_1,\mathrm{\dots },T_d$ by connecting $\{r,v_i\}$ for $1\le i\le d$; output $T$.

For the purpose of proof, we also consider the $2$-spin system on $G^{\prime }$ with the same parameters $(\beta ,\gamma ,𝝀)$, with an exception that we let the vertices $r_1,\mathrm{\dots },r_d$ have no fields (i.e., setting ${\lambda }_{r_i}=1$ for $1\le i\le d$ instead of ${\lambda }_r$). We then observe that

$$Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)={\lambda }_r\cdot Z_{G^{\prime }}^{{\sigma }_{\mathrm{\Lambda }}}(r_1=1,\mathrm{\dots },r_d=1),$$

and the same holds with spin $1$ replaced by $0$. For $1\le i\le d$, let ${\sigma }_{{\mathrm{\Lambda }}_i}$ denote the union of the conditioning ${\sigma }_{\mathrm{\Lambda }}$ and ${\sigma }_{U_i}$, where ${\mathrm{\Lambda }}_i=\mathrm{\Lambda }\cup U_i$. Then for every $1\le i\le d$ we have

$$Z_{G^{\prime }}^{{\sigma }_{\mathrm{\Lambda }}}(r_1=0,\mathrm{\dots },r_{i-1}=0,r_i=1,\mathrm{\dots },r_d=1)=\beta \cdot Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=1)+Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=0).$$

Notice that both sides are independent of the field ${\lambda }_r$: for the left side, all $r_i$’s do not have a field for the spin system on $G^{\prime }$; for the right side, recall that we do not count the weight of fixed vertices for the conditional partition function for each $G_i$. Now define $Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}(𝝀)$ by

$$Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=\prod _{i=2}^d{Z}_{G^{\prime }}^{{\sigma }_{\mathrm{\Lambda }}}(r_1=0,\mathrm{\dots },r_{i-1}=0,r_i=1,\mathrm{\dots },r_d=1),$$

which is independent of ${\lambda }_r$. Then we get

$Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)\cdot Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}$ $={\lambda }_r\cdot {\displaystyle \prod _{i=1}^d}{Z}_{G^{\prime }}^{{\sigma }_{\mathrm{\Lambda }}}(r_1=0,\mathrm{\dots },r_{i-1}=0,r_i=1,\mathrm{\dots },r_d=1)$

$={\lambda }_r\cdot {\displaystyle \prod _{i=1}^d}\left(\beta \cdot Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=1)+Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=0)\right).$

Using a similar argument, we also have

$Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=0)\cdot Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}$ $={\displaystyle \prod _{i=1}^d}{Z}_{G^{\prime }}^{{\sigma }_{\mathrm{\Lambda }}}(r_1=0,\mathrm{\dots },r_i=0,r_{i+1}=1,\mathrm{\dots },r_d=1)$

$={\displaystyle \prod _{i=1}^d}\left(Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=1)+\gamma \cdot Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=0)\right).$

Since we assume that $(G\backslash \{r\})\backslash \mathrm{\Lambda }$ is connected, the graph $G_i\backslash \mathrm{\Lambda }$ is also connected for each $i$. Then, by the induction hypothesis, for each $i$ there exists a polynomial $P_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}=P_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(𝝀)$ such that

$$Z_{T_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(r=1)=Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(r=1)\cdot P_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}\text{and}{Z}_{T_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(r=0)=Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(r=0)\cdot P_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}};$$

these polynomials are independent of ${\lambda }_r$ since the conditional partition functions for $G_i$’s do not involve ${\lambda }_r$. Now if we let

$$P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}\cdot \prod _{i=1}^d{P}_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}},$$

then it follows from the tree recursion that

$Z_T^{{\sigma }_{\mathrm{\Lambda }}}(r=1)$ $={\lambda }_r\cdot {\displaystyle \prod _{i=1}^d}\left(\beta \cdot Z_{T_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=1)+Z_{T_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=0)\right)$

$={\lambda }_r\cdot {\displaystyle \prod _{i=1}^d}\left(\beta \cdot Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=1)\cdot P_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}+Z_{G_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}(v_i=0)\cdot P_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}\right)$

$=Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)\cdot Q_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}\cdot {\displaystyle \prod _{i=1}^d}{P}_{G_i,v_i}^{{\sigma }_{{\mathrm{\Lambda }}_i}}$

$=Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}.$

The other equality $$ZTσΛ​(r​ = ​0)=ZGσΛ​(r​ = ​0)⋅PG,rσΛ is established in the same way. This completes the proof for the case that $(G\backslash \{r\})\backslash \mathrm{\Lambda }$ is connected.

If $(G\backslash \{r\})\backslash \mathrm{\Lambda }$ has two or more connected components, then we can construct $T_{\text{saw}}(G,r)$ by the SAW tree of each component. Recall that $G^{\prime }$ is defined by splitting the vertex $r$ into $d$ copies in the graph $G$. Suppose that $G^{\prime }\backslash \mathrm{\Lambda }$ has $k$ connected component for an integer $k\ge 2$. Let $G_{(1)}^{\prime },\mathrm{\dots },G_{(k)}^{\prime }$ be the subgraphs induced by each component, along with vertices from $\mathrm{\Lambda }$ that are adjacent to it. For each $j$, let $G_{(j)}$ be the graph obtained from $G_{(j)}^{\prime }$ by contracting all copies of $r$ into one vertex $r_{(j)}$, and let $T_{(j)}=T_{\text{saw}}(G_{(j)}^{\prime },r_{(j)})$. Observe that once we contract the roots $r_{(1)},\mathrm{\dots },r_{(k)}$ of $T_{(1)},\mathrm{\dots },T_{(k)}$, the resulting tree is $T_{\text{saw}}(G,r)$.

We define the $2$-spin system on each $G_{(j)}$ with the same parameters $(\beta ,\gamma ,𝝀)$, except that the vertex $r_{(j)}$ does not have a field (i.e., ${\lambda }_{r_{(j)}}=1$ instead of ${\lambda }_r$). For $1\le j\le k$, let ${\mathrm{\Lambda }}_{(j)}=\mathrm{\Lambda }\cap V(G_{(j)})$ and ${\sigma }_{{\mathrm{\Lambda }}_{(j)}}$ be the configuration ${\sigma }_{\mathrm{\Lambda }}$ restricted on ${\mathrm{\Lambda }}_{(j)}$. Then $G_{(j)}\backslash {\mathrm{\Lambda }}_{(j)}$ is connected for every $j$ and, since $k\ge 2$, each $G_{(j)}$ with conditioning ${\sigma }_{{\mathrm{\Lambda }}_{(j)}}$ has fewer than $m$ edges. Thus, we can apply the induction hypothesis; namely, for $1\le j\le k$ there exists a polynomial $P_{G_{(i)},r_{(i)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}=P_{G_{(i)},r_{(i)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(𝝀)$, which is independent of ${\lambda }_r$, such that

$$Z_{T_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(r_{(j)}=1)=Z_{G_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(r_{(j)}=1)\cdot P_{G_{(j)},r_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}\text{and}{Z}_{T_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(r_{(j)}=0)=Z_{G_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(r_{(j)}=0)\cdot P_{G_{(j)},r_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}.$$

We define the polynomial $P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}(𝝀)$ to be

$$P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}}=\prod _{j=1}^k{P}_{G_{(j)},r_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}.$$

It is then easy to check that

$Z_T^{{\sigma }_{\mathrm{\Lambda }}}(r=1)$ $={\lambda }_r\cdot {\displaystyle \prod _{j=1}^k}{Z}_{T_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(r_{(j)}=1)={\lambda }_r\cdot {\displaystyle \prod _{j=1}^k}\left(Z_{G_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}(r_{(j)}=1)\cdot P_{G_{(j)},r_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}\right)$

$=Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)\cdot {\displaystyle \prod _{j=1}^k}{P}_{G_{(j)},r_{(j)}}^{{\sigma }_{{\mathrm{\Lambda }}_{(j)}}}=Z_G^{{\sigma }_{\mathrm{\Lambda }}}(r=1)\cdot P_{G,r}^{{\sigma }_{\mathrm{\Lambda }}},$

and similarly $$ZTσΛ​(r​ = ​0)=ZGσΛ​(r​ = ​0)⋅PG,rσΛ. The theorem then follows. ∎

5 Influence bound for trees

In this section, we study the influences of the root on other vertices in a tree. We give an upper bound on the total influences of the root on all vertices at a fixed distance away. To do this, we apply the potential method, which has been used to establish the correlation decay property (see, e.g., [LLY12, LLY13, GL18]). Given an arbitrary potential function $\mathrm{\Psi }$, our upper bound is in terms of properties of $\mathrm{\Psi }$, involving bounds on ${\Vert \nabla H_d^{\mathrm{\Psi }}\Vert }_1$ and $|\psi |$ where $\psi ={\mathrm{\Psi }}^{\prime }$. We then deduce 9 in the case that $\mathrm{\Psi }$ an $(\alpha ,c)$-potential.

Assume that $T=(V_T,E_T)$ is a tree rooted at $r$ of maximum degree at most $\mathrm{\Delta }$. Let $\mathrm{\Lambda }\subseteq V_T\backslash \{r\}$ and ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ be arbitrary and fixed. Consider the $2$-spin system on $T$ with parameters $(\beta ,\gamma ,\lambda )$, conditioned on ${\sigma }_{\mathrm{\Lambda }}$. We need to bound the influence $$ℐTσΛ​(r​ → ​v) from the root $r$ to another vertex $v\in V_T$. Notice that if $v$ is disconnected from $r$ when $\mathrm{\Lambda }$ is removed, then $$ℐTσΛ​(r​ → ​v)=0 by the Markov property of spin systems. Therefore, we may assume that, by removing all such vertices, $\mathrm{\Lambda }$ contains only leaves of $T$.

For a vertex $v\in V_T$, let $T_v=(V_{T_v},E_{T_v})$ be the subtree of $T$ rooted at $v$ that contains all descendant of $v$; note that $T_r=T$. We will write $L_v(k)\subseteq V_T\backslash \mathrm{\Lambda }$ for the set of all free vertices at distance $k$ away from $v$ in $T_v$. We pay particular interest in the marginal ratio at $v$ in the subtree $T_v$, and write $R_v=R_{T_v}^{{\sigma }_{\mathrm{\Lambda }}}(v)$ for simplicity. The $\log R_v$’s are related by the tree recursion $H$. If a vertex $v$ has $d$ children, denoted by $u_1,\mathrm{\dots },u_d$, then the tree recursion is given by

$$\log R_v=H_d(\log R_{u_1},\mathrm{\dots },\log R_{u_d}),$$

where for $1\le d\le \mathrm{\Delta }$ and $(y_1,\mathrm{\dots },y_d)\in {[-\mathrm{\infty },+\mathrm{\infty }]}^d$,

$$H_d(y_1,\mathrm{\dots },y_d)=\log \lambda +\sum _{i=1}^d\log \left(\frac{\beta e^{y_i}+1}{e^{y_i}+\gamma }\right).$$

Also recall that for $y\in [-\mathrm{\infty },+\mathrm{\infty }]$, we define

$$h(y)=-\frac{(1-\beta \gamma )e^y}{(\beta e^y+1)(e^y+\gamma )}$$

and $\frac{\partial }{\partial y_i}{H}_d(y_1,\mathrm{\dots },y_d)=h(y_i)$ for all $1\le i\le d\le \mathrm{\Delta }$.

The following lemma allows us to bound the sum of all influences from the root to distance $k$, using an arbitrary potential function.

Lemma 14.

Let $\mathrm{\Psi }:[-\mathrm{\infty },+\mathrm{\infty }]\to (-\mathrm{\infty },+\mathrm{\infty })$ be a differentiable and increasing (potential) function with image $S=\mathrm{\Psi }[-\mathrm{\infty },+\mathrm{\infty }]$ and derivative $\psi ={\mathrm{\Psi }}^{\prime }$. Denote the degree of the root $r$ by ${\mathrm{\Delta }}_r$. Then for every integer $k\ge 1$,

where

$$A_{\mathrm{\Psi }}=\underset{u\in L_r(1)}{\max }\left\{\frac{|h(\log R_u)|}{\psi (\log R_u)}\right\}\text{and}{B}_{\mathrm{\Psi }}=\underset{v\in L_r(k)}{\max }\left\{\psi (\log R_v)\right\}.$$

Before proving 14, we first present two useful properties of the influences on trees. Firstly, it was shown in [ALO20] that the influences satisfy the following form of chain rule on trees.

Lemma 15 ([ALO20, Lemma B.2]).

Suppose that $u,v,w\in V_T$ are three distinct vertices such that $u$ is on the unique path from $v$ to $w$. Then

$${ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(v\to w)={ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(v\to u)\cdot {ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u\to w).$$

Secondly, for two adjacent vertices on a tree, the influence from one to the other is given by the function $h$.

Lemma 16.

Let $v\in V_T$ and $u$ be a child of $v$ in the subtree $T_v$. Then

$${ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(v\to u)=h(\log R_u).$$

Proof.

The lemma can be proved through an explicit computation of the influence. Here we present a more delicate proof utilizing 12, which gives some insights into the relation between the influence and the function $h$. We assume that $v$ has $d$ children in the subtree $T_v$, denoted by $u_1=u$ and $u_2,\mathrm{\dots },u_d$ respectively. We also assume, as a more general setting than uniform fields, that each vertex $w$ is attached to a field ${\lambda }_w$ of its own. Then 12 and the tree recursion imply that

${ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(v\to u)$ $={ℐ}_{T_v}^{{\sigma }_{\mathrm{\Lambda }}}(v\to u)=\left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)\log R_v$

$=\left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)H_d(\log R_{u_1},\mathrm{\dots },\log R_{u_d})$

$={\displaystyle \sum _{i=1}^d}{\displaystyle \frac{\partial }{\partial \log R_{u_i}}}{H}_d(\log R_{u_1},\mathrm{\dots },\log R_{u_d})\cdot \left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)\log R_{u_i}$

$={\displaystyle \sum _{i=1}^d}h(\log R_{u_i})\cdot {ℐ}_{T_{u_i}}^{{\sigma }_{\mathrm{\Lambda }}}(u_i\to u)=h(\log R_u),$

where the last equality is because $$ℐTuiσΛ​(ui​ → ​u)=0 for $u_i\ne u$ and $$ℐTuσΛ​(u​ → ​u)=1. ∎

We are now ready to prove 14.

Proof of 14.

For a vertex $v\in V_T$, denote the number of its children by $d_v$; note that $d_r={\mathrm{\Delta }}_r$. Let $u_1,\mathrm{\dots },u_{{\mathrm{\Delta }}_r}$ be the children of the root $r$. We may assume that all these children of $r$ are free, since if $u_i$ is fixed then $$ℐTσΛ​(r​ → ​ui)=0 by definition. Then by 15 and 16, we get

${\displaystyle \sum _{v\in L_r(k)}}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|$ $={\displaystyle \sum _{i=1}^{{\mathrm{\Delta }}_r}}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to u_i)\right|{\displaystyle \sum _{v\in L_{u_i}(k-1)}}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u_i\to v)\right|$

$={\displaystyle \sum _{i=1}^{{\mathrm{\Delta }}_r}}\left|h(\log R_{u_i})\right|{\displaystyle \sum _{v\in L_{u_i}(k-1)}}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u_i\to v)\right|$

$={\displaystyle \sum _{i=1}^{{\mathrm{\Delta }}_r}}{\displaystyle \frac{\left|h(\log R_{u_i})\right|}{\psi (\log R_{u_i})}}{\displaystyle \sum _{v\in L_{u_i}(k-1)}}\psi (\log R_{u_i})\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u_i\to v)\right|.$

Hence, we obtain that

$$\sum _{v\in L_r(k)}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|\le {\mathrm{\Delta }}_r\cdot \underset{1\le i\le {\mathrm{\Delta }}_r}{\max }\left\{\frac{|h(\log R_{u_i})|}{\psi (\log R_{u_i})}\right\}\cdot \underset{1\le i\le {\mathrm{\Delta }}_r}{\max }\left\{\sum _{v\in L_{u_i}(k-1)}\psi (\log R_{u_i})\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u_i\to v)\right|\right\}.$$ (6)

Next, we show by induction that for every vertex $u\in V_T\backslash \{r\}$ and every integer $k\ge 0$ we have

$$\sum _{v\in L_u(k)}\psi (\log R_u)\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u\to v)\right|\le \underset{v\in L_u(k)}{\max }\left\{\psi (\log R_v)\right\}\cdot {\left(\underset{w\in V_{T_u}}{\max }\underset{\widetilde{𝒚}\in S^{d_w}}{sup}{\Vert \nabla H_{d_w}^{\mathrm{\Psi }}(\widetilde{𝒚})\Vert }_1\right)}^k.$$ (7)

Observe that once we establish Eq. 7, the lemma follows immediately by plugging Eq. 7 into Eq. 6. We will use induction on $k$ to prove Eq. 7. When $k=0$, if $u\in \mathrm{\Lambda }$ is fixed then $L_u(0)=\mathrm{\varnothing }$ and there is nothing to show; otherwise, Eq. 7 becomes

$$\psi (\log R_u)\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u\to u)\right|\le \psi (\log R_u),$$

which holds with equality since $$ℐTσΛ​(u​ → ​u)=1. Now suppose that Eq. 7 holds for some integer $k-1\ge 0$ (and for every vertex $u\in V_T\backslash \{r\}$). Let $u\in V_T\backslash \{r\}$ be arbitrary and denote the children of $u$ by $w_1,\mathrm{\dots },w_d$, where (if $d=0$ then $L_u(k)=\mathrm{\varnothing }$ and Eq. 7 holds trivially). Again by 15 and 16 we have

${\displaystyle \sum _{v\in L_u(k)}}\psi (\log R_u)\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u\to v)\right|$ $={\displaystyle \sum _{i=1}^d}\psi (\log R_u)\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u\to w_i)\right|{\displaystyle \sum _{v\in L_{w_i}(k-1)}}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(w_i\to v)\right|$

$={\displaystyle \sum _{i=1}^d}{\displaystyle \frac{\psi (\log R_u)}{\psi (\log R_{w_i})}}\left|h(\log R_{w_i})\right|{\displaystyle \sum _{v\in L_{w_i}(k-1)}}\psi (\log R_{w_i})\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(w_i\to v)\right|.$

Using the induction hypothesis, we get

${\displaystyle \sum _{v\in L_u(k)}}\psi (\log R_u)\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(u\to v)\right|$

$\le$ ${\displaystyle \sum _{i=1}^d}{\displaystyle \frac{\psi (\log R_u)}{\psi (\log R_{w_i})}}\left|h(\log R_{w_i})\right|\cdot \underset{v\in L_{w_i}(k-1)}{\max }\left\{\psi (\log R_v)\right\}\cdot {\left(\underset{w\in V_{T_{w_i}}}{\max }\underset{\widetilde{𝒚}\in S^{d_w}}{sup}{\Vert \nabla H_{d_w}^{\mathrm{\Psi }}(\widetilde{𝒚})\Vert }_1\right)}^{k-1}$

$\le$ $\underset{v\in L_u(k)}{\max }\left\{\psi (\log R_v)\right\}\cdot {\left(\underset{w\in V_{T_u}\backslash \{u\}}{\max }\underset{\widetilde{𝒚}\in S^{d_w}}{sup}{\Vert \nabla H_{d_w}^{\mathrm{\Psi }}(\widetilde{𝒚})\Vert }_1\right)}^{k-1}\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{\psi (\log R_u)}{\psi (\log R_{w_i})}}\left|h(\log R_{w_i})\right|$

$\le$ $\underset{v\in L_u(k)}{\max }\left\{\psi (\log R_v)\right\}\cdot {\left(\underset{w\in V_{T_u}}{\max }\underset{\widetilde{𝒚}\in S^{d_w}}{sup}{\Vert \nabla H_{d_w}^{\mathrm{\Psi }}(\widetilde{𝒚})\Vert }_1\right)}^k,$

where the last inequality follows from that

${\displaystyle \sum _{i=1}^d}{\displaystyle \frac{\psi (\log R_u)}{\psi (\log R_{w_i})}}\left|h(\log R_{w_i})\right|$ $={\displaystyle \sum _{i=1}^d}\left|{\displaystyle \frac{\partial }{\partial \mathrm{\Psi }(\log R_{w_i})}}{H}_d^{\mathrm{\Psi }}(\mathrm{\Psi }(\log R_{w_1}),\mathrm{\dots },\mathrm{\Psi }(\log R_{w_d}))\right|$

$={\Vert \nabla H_d^{\mathrm{\Psi }}(\mathrm{\Psi }(\log R_{w_1}),\mathrm{\dots },\mathrm{\Psi }(\log R_{w_d}))\Vert }_1.$

This establishes Eq. 7, and thus completes the proof of the lemma. ∎

We then derive 9 as a corollary.

Proof of 9.

Since $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential, the Contraction condition implies that

Meanwhile, since the degree of a vertex $v\in V_T\backslash \{r\}$ in the subtree $T_v$ is less than $\mathrm{\Delta }$, we have $\log R_v\in J$. Then the Boundedness condition implies that for all $u\in L_r(1)$ and $v\in L_r(k)$,

$$\frac{\psi (\log R_v)}{\psi (\log R_u)}\cdot |h(\log R_u)|\le \frac{c}{\mathrm{\Delta }}.$$

Therefore, we get

$${\mathrm{\Delta }}_r{A}_{\mathrm{\Psi }}{B}_{\mathrm{\Psi }}={\mathrm{\Delta }}_r\cdot \underset{u\in L_r(1)}{\max }\left\{\frac{|h(\log R_u)|}{\psi (\log R_u)}\right\}\cdot \underset{v\in L_r(k)}{\max }\left\{\psi (\log R_v)\right\}\le c.$$

The lemma then follows immediately from 14. ∎

6 Verifying a good potential: Contraction

In this section, we make a first step for proving 10. Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. Recall that define our potential function $\mathrm{\Psi }:[-\mathrm{\infty },+\mathrm{\infty }]\to (-\mathrm{\infty },+\mathrm{\infty })$ through its derivative by

$${\mathrm{\Psi }}^{\prime }(y)=\psi (y)=\sqrt{\frac{(1-\beta \gamma )e^y}{(\beta e^y+1)(e^y+\gamma )}},\mathrm{\Psi }(0)=0.$$ (1)

We include a short proof in Appendix C to show that $\mathrm{\Psi }$ is well-defined. If $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$, then we show that $\mathrm{\Psi }$ satisfies the Contraction condition for $\alpha =\delta /2$. This holds for all parameters $(\beta ,\gamma ,\lambda )$ in the uniqueness region, without requiring that $\gamma \le 1$. Later in Appendix E, we establish the Boundedness condition for $\mathrm{\Psi }$ when $\gamma \le 1$, completing the proof of 10. The case of $\gamma >1$ is more complicated and is left to Section 7.

Before giving our proof, we first point out that the potential function $\mathrm{\Psi }$ is essentially the same potential function $\mathrm{\Phi }$ used in [LLY13] (notice that [LLY13] uses $\phi$ as the notation of the potential function and $\mathrm{\Phi }={\phi }^{\prime }$ for its derivative). Recall that the tree recursion for the marginal ratios is given by the function $F_d:{[0,+\mathrm{\infty }]}^d\to [0,+\mathrm{\infty }]$ where $1\le d\le \mathrm{\Delta }$ such that for all $(x_1,\mathrm{\dots },x_d)\in {[0,+\mathrm{\infty }]}^d$,

$$F_d(x_1,\mathrm{\dots },x_d)=\lambda \prod _{i=1}^d\frac{\beta x_i+1}{x_i+\gamma }.$$

The potential function $\mathrm{\Phi }:[0,+\mathrm{\infty }]\to (-\mathrm{\infty },+\mathrm{\infty })$ from [LLY13] is defined implicitly via its derivative as

$${\mathrm{\Phi }}^{\prime }(x)=\phi (x)=\frac{1}{\sqrt{x(\beta x+1)(x+\gamma )}},\mathrm{\Phi }(1)=0.$$

The follows lemma explains how we obtain our potential $\mathrm{\Psi }$ from $\mathrm{\Phi }$.

Lemma 17.

We have $\mathrm{\Psi }=\sqrt{1-\beta \gamma }\cdot (\mathrm{\Phi }\circ \exp )$; namely, $\mathrm{\Psi }(y)=\sqrt{1-\beta \gamma }\cdot \mathrm{\Phi }(e^y)$ for all $y\in [-\mathrm{\infty },+\mathrm{\infty }]$.

Proof.

It is straightforward to check that

$$\psi (y)=\sqrt{\frac{(1-\beta \gamma )e^y}{(\beta e^y+1)(e^y+\gamma )}}=\sqrt{1-\beta \gamma }\cdot e^y\cdot \sqrt{\frac{1}{e^y(\beta e^y+1)(e^y+\gamma )}}=\sqrt{1-\beta \gamma }\cdot e^y\phi (e^y).$$

Therefore,

$$\mathrm{\Psi }(y)={\int }_0^y\psi (t)dt=\sqrt{1-\beta \gamma }\cdot {\int }_0^y{e}^t\phi (e^t)dt=\sqrt{1-\beta \gamma }\cdot {\int }_1^{e^y}\phi (s)ds=\sqrt{1-\beta \gamma }\cdot \mathrm{\Phi }(e^y).\mathit{∎}$$

Combining the results of Lemmas 12, 13 and 14 from [LLY13], we get that the potential function $\mathrm{\Phi }$ satisfies the following gradient bound when $(\beta ,\gamma ,\lambda )$ is in the uniqueness region. Note that this can be regarded as the Contraction condition but for $\mathrm{\Phi }$ and $F_d$.

Theorem 18 ([LLY13]).

Let $S_{\mathrm{\Phi }}=\mathrm{\Phi }[0,+\mathrm{\infty }]$ be the image of $\mathrm{\Phi }$. If the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$, then for every integer $d$ such that and every $({\tilde{x}}_1,\mathrm{\dots },{\tilde{x}}_d)\in S_{\mathrm{\Phi }}^d$,

$${\Vert \nabla F_d^{\mathrm{\Phi }}({\tilde{x}}_1,\mathrm{\dots },{\tilde{x}}_d)\Vert }_1\le \sqrt{1-\delta }$$

where $F_d^{\mathrm{\Phi }}=\mathrm{\Phi }\circ F_d\circ {\mathrm{\Phi }}^{-1}$.

Recall our definition from Section 1.1. The tree recursion, in terms of the log marginal ratios, is described by the function $H_d:{[-\mathrm{\infty },+\mathrm{\infty }]}^d\to [-\mathrm{\infty },+\mathrm{\infty }]$ where $1\le d\le \mathrm{\Delta }$ such that for every $(y_1,\mathrm{\dots },y_d)\in {[-\mathrm{\infty },+\mathrm{\infty }]}^d$,

$$H_d(y_1,\mathrm{\dots },y_d)=\log \lambda +\sum _{i=1}^d\log \left(\frac{\beta e^{y_i}+1}{e^{y_i}+\gamma }\right).$$

Observe that $H_d=\log \circ F_d\circ \exp$, since we move from ratios to log ratios. We are now ready to establish the Contraction condition for $\mathrm{\Psi }$.

Lemma 19.

Let $S_{\mathrm{\Psi }}=\mathrm{\Psi }[-\mathrm{\infty },+\mathrm{\infty }]$ be the image of $\mathrm{\Psi }$. If the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$, then for every integer $d$ such that and every $({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\in S_{\mathrm{\Psi }}^d$,

$${\Vert \nabla H_d^{\mathrm{\Psi }}({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\Vert }_1\le \sqrt{1-\delta }$$

where $H_d^{\mathrm{\Psi }}=\mathrm{\Psi }\circ H_d\circ {\mathrm{\Psi }}^{-1}$.

Proof.

Define the linear function $a:ℝ\to ℝ$ to be $a(x)=\sqrt{1-\beta \gamma }\cdot x$ for $x\in ℝ$. Then 17 gives $\mathrm{\Psi }=a\circ \mathrm{\Phi }\circ \exp$, and thereby $\mathrm{\Psi }\circ \log =a\circ \mathrm{\Phi }$. It follows that for every ,

$$H_d^{\mathrm{\Psi }}=\mathrm{\Psi }\circ H_d\circ {\mathrm{\Psi }}^{-1}=\mathrm{\Psi }\circ \log \circ F_d\circ \exp \circ {\mathrm{\Psi }}^{-1}=a\circ \mathrm{\Phi }\circ F_d\circ {\mathrm{\Phi }}^{-1}\circ a^{-1}=a\circ F_d^{\mathrm{\Phi }}\circ a^{-1}.$$

That means, for every $({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\in S_{\mathrm{\Psi }}^d$ we have

$$H_d^{\mathrm{\Psi }}({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)=\sqrt{1-\beta \gamma }\cdot F_d^{\mathrm{\Phi }}({\tilde{x}}_1,\mathrm{\dots },{\tilde{x}}_d)$$

where ${\tilde{x}}_i={\tilde{y}}_i/\sqrt{1-\beta \gamma }$ for $1\le i\le d$. Then, for each $i$,

$$\frac{\partial }{\partial {\tilde{y}}_i}{H}_d^{\mathrm{\Psi }}({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)=\sqrt{1-\beta \gamma }\cdot \frac{\partial }{\partial {\tilde{x}}_i}{F}_d^{\mathrm{\Phi }}({\tilde{x}}_1,\mathrm{\dots },{\tilde{x}}_d)\cdot \frac{\mathrm{d}{\tilde{x}}_i}{\mathrm{d}{\tilde{y}}_i}=\frac{\partial }{\partial {\tilde{x}}_i}{F}_d^{\mathrm{\Phi }}({\tilde{x}}_1,\mathrm{\dots },{\tilde{x}}_d).$$

This implies that $\nabla H_d^{\mathrm{\Psi }}({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)=\nabla F_d^{\mathrm{\Phi }}({\tilde{x}}_1,\mathrm{\dots },{\tilde{x}}_d)$ for all $({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\in S_{\mathrm{\Psi }}^d$, and the lemma then follows from 18. ∎

7 Remaining antiferromagnetic cases: $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$

In this section, we discuss the case where $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$. As studied in [LLY13], in this case the uniqueness region is more complicated. For example, there exists a critical ${\lambda }_c^{\ast }>0$ such that the $2$-spin system with is in the uniqueness region for arbitrary graphs; namely, $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\infty }$ unique. To deal with large degrees, we need to relax the Boundedness condition in 4 and define a more general version of $(\alpha ,c)$-potentials. We shall see that 5 still holds for this general $(\alpha ,c)$-potential. The reason behind it is that in order to bound the maximum eigenvalue of the influence matrix, it suffices to consider a vertex-weighted sum of absolute influences of a vertex with large degree.

Remark 20.

We give more background on the uniqueness region in Section E.1. Note that in a recent revision of [LLY13], the authors updated the descriptions of the uniqueness region for the case $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$, fixing a small error in the previous version. Statements and proofs in this section and Appendix E of this paper are also adjusted accordingly based on the new version of [LLY13].

Recall that our goal is to bound the maximum eigenvalue of the matrix ${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}$. We can do this by upper bounding the absolute row sum ${\sum }_{v\in V\backslash \mathrm{\Lambda }}|$ℐGσΛ​(r​ → ​v)| for fixed $r$, thereby giving us a valid upper bound on ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})$. However, this approach does not work when $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$. In this case, the potential $\mathrm{\Psi }$ fails to be an $(\alpha ,c)$-potential for a universal constant $c$ independent of $\mathrm{\Delta }$. In fact, no such $(\alpha ,c)$-potentials exist as the absolute row sum ${\sum }_{v\in V\backslash \mathrm{\Lambda }}|$ℐGσΛ​(r​ → ​v)| can be as large as $\mathrm{\Theta }(\mathrm{\Delta })$. Especially, if the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\infty }$ unique, which means the spin system has uniqueness for arbitrary graphs, then the absolute row sum ${\sum }_{v\in V\backslash \mathrm{\Lambda }}|$ℐGσΛ​(r​ → ​v)| can be $\mathrm{\Theta }(n)$ where $n=|V|$. We give a specific example where this is the case.

Example 21.

Consider the antiferromagnetic 2-spin system specified by parameters $\beta =0$, $\gamma >1$ and $\lambda >0$ on the star graph centered at $r$ with $\mathrm{\Delta }$ leaves. A simple calculation reveals that $|$ℐG​(r​ → ​v)|=λλ+γ for any leaf vertex $v\ne r$. Hence, ${\sum }_{v\ne r}|$ℐG​(r​ → ​v)|=Δ⋅λλ+γ. Now, since $\gamma >1$, we have

forcing ${\sum }_{v\ne r}|$ℐG​(r​ → ​v)|=Θγ​(Δ) even when lies in the uniqueness region. However, we still have ${\lambda }_{\max }({ℐ}_G)=O(1)$ since ${\sum }_{v\ne r}|$ℐG​(v​ → ​r)|=O​(1).

To solve this issue, one might want to consider the absolute column sum, involving the sum of absolute influences on a fixed vertex. However, this will not allow us to use the beautiful connection between graphs and SAW trees as showed in 8. Instead, we consider here a vertex-weighted version of the absolute row sum of ${ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}$, which also upper bounds the maximum eigenvalue.

Lemma 22.

Let $\rho :V\to {ℝ}^{+}$ be a positive weight function of vertices. If there is a constant $\xi >0$ such that for every $r\in V$ we have

$$\sum _{v\in V\backslash \mathrm{\Lambda }}{\rho }_v\cdot \left|{ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|\le \xi \cdot {\rho }_r,$$ (8)

then ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})\le \xi$.

Proof.

Let $𝒫=\mathrm{diag}\{{\rho }_v:v\in V\backslash \mathrm{\Lambda }\}$. Then the assumption is equivalent to ${\Vert {𝒫}^{-1}{ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}𝒫\Vert }_{\mathrm{\infty }}\le \xi$. It follows that ${\lambda }_{\max }({ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}})={\lambda }_{\max }({𝒫}^{-1}{ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}𝒫)\le \xi$. ∎

We then modify our definition of $(\alpha ,c)$-potentials from 4 which allows a weaker Boundedness condition. We remark that the only two differences between 23 and 4 is that: we allow $\mathrm{\Delta }=\mathrm{\infty }$; and the Boundedness condition is relaxed to what we call General Boundedness. Recall that for every , we let $J_d=[\log (\lambda {\beta }^d),\log (\lambda /{\gamma }^d)]$ when , and $J_d=[\log (\lambda /{\gamma }^d),\log (\lambda {\beta }^d)]$ when $\beta \gamma >1$.

Definition 23 (General $(\alpha ,c)$-potential function).

Let $\mathrm{\Delta }\ge 3$ be an integer or $\mathrm{\Delta }=\mathrm{\infty }$. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$ and $\lambda >0$. Let $\mathrm{\Psi }:[-\mathrm{\infty },+\mathrm{\infty }]\to (-\mathrm{\infty },+\mathrm{\infty })$ be a differentiable and increasing function with image $S=\mathrm{\Psi }[-\mathrm{\infty },+\mathrm{\infty }]$ and derivative $\psi ={\mathrm{\Psi }}^{\prime }$. For any $\alpha \in (0,1)$ and $c>0$, we say $\mathrm{\Psi }$ is a general $(\alpha ,c)$-potential function with respect to $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$ if it satisfies the following conditions:

1. (Contraction) For every integer $d$ such that and every $({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\in S^d$, we have

   $${\Vert \nabla H_d^{\mathrm{\Psi }}({\tilde{y}}_1,\mathrm{\dots },{\tilde{y}}_d)\Vert }_1=\sum _{i=1}^d\frac{\psi (y)}{\psi (y_i)}\cdot |h(y_i)|\le 1-\alpha$$

   where $H_d^{\mathrm{\Psi }}=\mathrm{\Psi }\circ H_d\circ {\mathrm{\Psi }}^{-1}$, $y_i={\mathrm{\Psi }}^{-1}({\tilde{y}}_i)$ for $1\le i\le d$, and $y=H_d(y_1,\mathrm{\dots },y_d)$.

2. (General Boundedness) For all integers $d_1,d_2$ such that , and all reals $y_1\in J_{d_1},y_2\in J_{d_2}$, we have

   $$\frac{\psi (y_2)}{\psi (y_1)}\cdot |h(y_1)|\le \frac{2c}{d_1+d_2+2}.$$

Notice that General Boundedness is a weaker condition than Boundedness. To see this, if a potential function $\mathrm{\Psi }$ satisfies Boundedness with parameter $c$, then for every and every $y_i\in J_{d_i}$ where $i=1,2$ we have

$$\frac{\psi (y_2)}{\psi (y_1)}\cdot |h(y_1)|\le \frac{c}{\mathrm{\Delta }}\le \frac{2c}{d_1+d_2+2}.$$

The following theorem generalizes 5 and shows that a general $(\alpha ,c)$-potential function is sufficient to establish rapid mixing of the Glauber dynamics.

Theorem 24.

Let $\mathrm{\Delta }\ge 3$ be an integer or $\mathrm{\Delta }=+\mathrm{\infty }$. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$ and $\lambda >0$. Suppose that there is a general $(\alpha ,c)$-potential with respect to $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$ for some $\alpha \in (0,1)$ and $c>0$. Then for every $n$-vertex graph $G$ of maximum degree $\mathrm{\Delta }$, the mixing time of the Glauber dynamics for the $2$-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $O(n^{2+2c/\alpha })$.

We then give a counterpart of 10, showing that $\mathrm{\Psi }$ is a general $(\alpha ,c)$-potential when $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$. 3 for this case is then obtained from 24 and 25.

Lemma 25.

Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that and $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$. Assume that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$. Then the function $\mathrm{\Psi }$ defined implicitly by Eq. 4 is a general $(\alpha ,c)$-potential function with $\alpha \ge \delta /2$ and $c\le 18$; we can further take $c\le 4$ if $\beta =0$.

The proof of 24 can be found in Appendix D. For 25, the Contraction condition of $\mathrm{\Psi }$ follows from 19, and General Boundedness is proved in Appendix E together with all other cases.

8 Ferromagnetic cases

In the ferromagnetic case, the best known correlation decay results are given in [GL18, SS20]. Using the potential functions in [GL18] and [SS20], we show the following two results, which match the known correlation decay results. In fact, the potential function from [SS20] turns out to be an $(\alpha ,c)$-potential function for constants $\alpha =\mathrm{\Theta }(\delta )$ and $c\le O(1)$.

Theorem 26.

Fix an integer $\mathrm{\Delta }\ge 3$, positive real numbers $\beta ,\gamma ,\lambda$ and , and assume $(\beta ,\gamma ,\lambda )$ satisfies one of the following three conditions:

1. $\frac{\mathrm{\Delta }-2+\delta }{\mathrm{\Delta }-\delta }\le \sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-\delta }{\mathrm{\Delta }-2+\delta }$, and $\lambda$ is arbitrary;

2. $\sqrt{\beta \gamma }\ge \frac{\mathrm{\Delta }}{\mathrm{\Delta }-2}$ and $\lambda \le (1-\delta )\frac{\gamma }{\max \{1,{\beta }^{\mathrm{\Delta }-1}\}\cdot ((\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta })}$;

3. $\sqrt{\beta \gamma }\ge \frac{\mathrm{\Delta }}{\mathrm{\Delta }-2}$ and $\lambda \ge \frac{1}{1-\delta }\cdot \frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta \cdot \min \{1,1/{\gamma }^{\mathrm{\Delta }-1}\}}$.

Then the identity function $\mathrm{\Psi }(y)=y$ (based on the potential given in [SS20]) is an $(\alpha ,c)$-potential function for $\alpha =\mathrm{\Theta }(\delta )$ and $c\le O(1)$. Furthermore, for every $n$-vertex graph $G$ of maximum degree at most $\mathrm{\Delta }$, the mixing time of the Glauber dynamics for the 2-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $O(n^{2+c/\delta })$, for a universal constant $c>0$.

Remark 3.

Condition 1 includes both the ferromagnetic case and the antiferromagnetic case . Note that in both cases $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap $\delta$. For the antiferromagnetic case, the identity function $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential with $c\le 1.5$ and a better contraction rate $\alpha \ge \delta$, compared with the bound $\alpha \ge \delta /2$ of the potential $\mathrm{\Psi }$ given by Eq. 4 in 10. For the ferromagnetic case with $\beta =\gamma >1$ (Ising model), a stronger result by [MS13] was known, which gives $O(n\log n)$ mixing.

The potential function from [GL18] is indeed an $(\alpha ,c)$-potential, but $c$ must, unfortunately, depend on $\mathrm{\Delta }$. We have the following result, which is weaker than the correlation decay algorithm in [GL18] for unbounded degree graphs.

Theorem 27.

Fix an integer $\mathrm{\Delta }\ge 3$, and nonnegative real numbers $\beta ,\gamma ,\lambda$ satisfying $\beta \le 1\le \gamma$, $\sqrt{\beta \gamma }\ge \frac{\mathrm{\Delta }}{\mathrm{\Delta }-2}$, and . Then for every $n$-vertex graph $G$ with maximum degree at most $\mathrm{\Delta }$, the mixing time of the Glauber dynamics for the ferromagnetic 2-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $O(n^C)$, for a constant $C$ depending only on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$, but not $n$.

Proofs of these theorems are provided in Appendix F.

Appendix A Proof of main results

In this section we give the proofs of 1, 2, 3 and 5.

Proof of 5.

Note that since the transition matrix $P$ for the Glauber dynamics has all nonnegative eigenvalues, we have that ${\lambda }^{\ast }(P)=1-{\lambda }_2(P)$ and so in order to deduce mixing, it suffices to lower bound $1-{\lambda }_2(P)$. We do this by employing 7. It suffices to show $({\eta }_0,\mathrm{\dots },{\eta }_{n-2})$-spectrally independence for sufficiently small ${\eta }_i$.

To bound ${\eta }_i$, it suffices to bound ${\sum }_{v\in V\backslash \{r\}}|$ℐGσΛ​(r​ → ​v)| for all graphs $G=(V,E)$ with $n=|V|$ vertices and all boundary conditions ${\sigma }_{\mathrm{\Lambda }}$ on a subset $\mathrm{\Lambda }$ of $i$ vertices. We claim the following:

$$\sum _{v\in V\backslash \{r\}}\left|{ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|\le \min \{\frac{c}{\alpha },C(n-i-1)\}$$ (9)

where $C\in (0,1)$ is a constant depending only on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$. The first upper bound $\frac{c}{\delta }$ is deduced by

${\displaystyle \sum _{v\in V\backslash \{r\}}}\left|{ℐ}_G^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|$ $\le {\displaystyle \sum _{v\in V_T\backslash \{r\}}}\left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|$ (8; $T=T_{\text{saw}}(G,r)$)

$={\displaystyle \sum _{k=1}^{\mathrm{\infty }}}{\displaystyle \sum _{v\in L_r(k)}}\left|{ℐ}_T^{{\sigma }_{{\mathrm{\Lambda }}_{\mathrm{\Lambda }}}}(r\to v)\right|$ (split the sum by levels)

$\le c{\displaystyle \sum _{k=1}^{\mathrm{\infty }}}{(1-\alpha )}^{k-1}$ (9)

$={\displaystyle \frac{c}{\alpha }}.$

The second upper bound $C(n-i-1)$ is more trivial. Intuitively, it means each absolute pairwise influence $|$ℐGσΛ​(r​ → ​v)| is at most some constant $C$ and hence the sum of absolute influences is upper bounded by $C(n-i-1)$. The following two claims, whose proofs are provided in Section A.2, give a more precise statement.

Claim 28 (Antiferromagnetic Case).

Fix an integer $\mathrm{\Delta }\ge 3$ and real numbers $\beta ,\gamma ,\lambda$, and assume $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. Then for every $n$-vertex graph $G$ of maximum degree at most $\mathrm{\Delta }$, the antiferromagnetic $2$-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $Cn$-spectrally independent, for a constant depending only on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$. Furthermore, if $(\beta ,\gamma ,\mathrm{\Delta })$ is up-to-$\mathrm{\Delta }$ unique, then we can drop the dependence on $\mathrm{\Delta }$.

Claim 29 (Ferromagnetic Case).

Fix an integer $\mathrm{\Delta }\ge 3$ and positive real numbers $\beta ,\gamma ,\lambda$, and assume $\beta \le \gamma$ and $\beta \gamma >1$. Then for every $n$-vertex graph $G$ of maximum degree at most $\mathrm{\Delta }$, the ferromagnetic $2$-spin system on $G$ with parameters $(\beta ,\gamma ,\lambda )$ is $Cn$-spectrally independent, for a constant depending only on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$.

With Eq. 9 in hand, we immediately see that by 7,

$$1-{\lambda }_2(P)\ge \frac{1}{n}\prod _{i=0}^{n-2}\left(1-\frac{{\eta }_i}{n-i-1}\right)\ge \frac{1}{n}\cdot {(1-C)}^{2\lceil c/\alpha \rceil -1}\cdot \prod _{i=0}^{n-2\lceil c/\alpha \rceil -1}\left(1-\frac{c}{\alpha }\cdot \frac{1}{n-i-1}\right).$$

Using the fact that $1-x\ge \exp (-x-x^2)$ for all $0\le x\le \frac{1}{2}$ (which can be proved straightforwardly by calculus), we get

${\displaystyle \prod _{i=0}^{n-2\lceil c/\alpha \rceil -1}}\left(1-{\displaystyle \frac{c}{\alpha }}\cdot {\displaystyle \frac{1}{n-i-1}}\right)={\displaystyle \prod _{j=2\lceil c/\alpha \rceil }^{n-1}}\left(1-{\displaystyle \frac{c}{\alpha }}\cdot {\displaystyle \frac{1}{j}}\right)\ge \exp \left(-{\displaystyle \frac{c}{\alpha }}{\displaystyle \sum _{j=2\lceil c/\alpha \rceil }^{n-1}}{\displaystyle \frac{1}{j}}-{\displaystyle \frac{c^2}{{\alpha }^2}}{\displaystyle \sum _{j=2\lceil c/\alpha \rceil }^{n-1}}{\displaystyle \frac{1}{j^2}}\right).$

Now since

$$\sum _{j=2\lceil c/\alpha \rceil }^{n-1}\frac{1}{j}\le \sum _{j=2}^n\frac{1}{j}\le {\int }_1^n\frac{dx}{x}=\log n$$

and

$$\sum _{j=2\lceil c/\alpha \rceil }^{n-1}\frac{1}{j^2}\le \sum _{j=2}^{\mathrm{\infty }}\frac{1}{j(j-1)}=1,$$

we deduce that

$$1-{\lambda }_2(P)\ge {(1-C)}^{2\lceil c/\alpha \rceil -1}\cdot e^{-{(c/\alpha )}^2}\cdot n^{-(1+c/\alpha )}.$$

The theorem then follows from Eq. 1. ∎

Proof of 3.

We leverage 5 and 24, which shows $O(n^{2+\frac{c}{\alpha }})$ mixing as long as there is an $(\alpha ,c)$-potential, or $O(n^{2+\frac{2c}{\alpha }})$ mixing if there is a general $(\alpha ,c)$-potential. We use the potential given by Eq. 4, which is an adaptation of the potential function in [LLY13] to the log marginal ratios. When $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$, it is an $(\alpha ,c)$-potential or a general $(\alpha ,c)$-potential by 10 and 25, with $\alpha \ge \delta /2$ and $c$ a universal constant specified by the range of parameters. The theorem then follows. ∎

Proof of 1.

By 30 later in Section A.1, $\lambda \le (1-\delta ){\lambda }_c(\mathrm{\Delta })$ implies up-to-$\mathrm{\Delta }$ uniqueness with gap $\ge \delta /4$. Since $\gamma \le 1$, we can again appeal to 10 to obtain an $(\alpha ,c)$-potential with $\alpha \ge \delta /8$ and $c\le 4$. 1 then follows by 5 with $O(n^{2+32/\delta })$ mixing. ∎

Proof of 2.

By 31 later in Section A.1, $\beta \ge {\beta }_c(\mathrm{\Delta })+\delta (1-{\beta }_c(\mathrm{\Delta }))$ implies up-to-$\mathrm{\Delta }$ uniqueness with gap $\delta$. Again, appealing to 10, we obtain an $(\alpha ,c)$-potential with $\alpha \ge \delta /2$ and $c\le 1.5$. 2 then follows by 5 with $O(n^{2+3/\delta })$ mixing.

Though we technically get $O(n^{2+3/\delta })$ by using the [LLY13] potential, we can improve it to $O(n^{2+1.5/\delta })$ mixing by using the trivial identity function as the potential. See the first case of 26 (proved in Section F.1) and Remark 3. ∎

A.1 Uniqueness gaps in terms of parameter paps

In this section we state and prove 30 and 31, which relate the parameter gaps with the uniqueness gaps.

Claim 30 (Hardcore Model; Lemma C.1 from [ALO20]).

Fix an integer $\mathrm{\Delta }\ge 3$, , and $\beta =0,\gamma >0$. If $\lambda \le (1-\delta ){\lambda }_c(\gamma ,\mathrm{\Delta })$, then $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap $\delta /4$.

Claim 31 (Large $\sqrt{\beta \gamma }$).

Fix an integer $\mathrm{\Delta }\ge 3$, and . If $\sqrt{\beta \gamma }\ge \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}+\delta \left(1-\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}\right)=\frac{\mathrm{\Delta }-2(1-\delta )}{\mathrm{\Delta }}$, then $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap for all $\lambda$. Note if $\beta =\gamma$, this is precisely the condition $\beta \ge {\beta }_c(\mathrm{\Delta })+\delta (1-{\beta }_c(\mathrm{\Delta }))$.

Proof.

Consider the univariate recursion for the marginal ratios with children $f_d(R)=\lambda {\left(\frac{\beta R+1}{R+\gamma }\right)}^d$. Differentiating, we have

$f_d^{\prime }(R)$ $=d\lambda {\left({\displaystyle \frac{\beta R+1}{R+\gamma }}\right)}^{d-1}\cdot \left({\displaystyle \frac{\beta }{R+\gamma }}-{\displaystyle \frac{\beta R+1}{{(R+\gamma )}^2}}\right)=-d(1-\beta \gamma )\lambda {\left({\displaystyle \frac{\beta R+1}{R+\gamma }}\right)}^d\cdot {\displaystyle \frac{1}{(\beta R+1)(R+\gamma )}}$

$=-d(1-\beta \gamma )\cdot {\displaystyle \frac{f_d(R)}{(\beta R+1)(R+\gamma )}}.$

At the unique fixed point $R_d^{\ast }$, we have $f_d(R_d^{\ast })=R_d^{\ast }$ so

$\left|f_d^{\prime }(R_d^{\ast })\right|=d(1-\beta \gamma ){\displaystyle \frac{R_d^{\ast }}{(\beta R_d^{\ast }+1)(R_d^{\ast }+\gamma )}}.$

By 37, we have the upper bound

$\left|f_d^{\prime }(R_d^{\ast })\right|\le d\cdot {\displaystyle \frac{1-\beta \gamma }{{(1+\sqrt{\beta \gamma })}^2}}=d\cdot {\displaystyle \frac{1-\sqrt{\beta \gamma }}{1+\sqrt{\beta \gamma }}}.$

Since we assumed $\sqrt{\beta \gamma }\ge \frac{\mathrm{\Delta }-2(1-\delta )}{\mathrm{\Delta }}$, we obtain

$d\cdot {\displaystyle \frac{1-\sqrt{\beta \gamma }}{1+\sqrt{\beta \gamma }}}\le d\cdot {\displaystyle \frac{\mathrm{\Delta }-(\mathrm{\Delta }-2(1-\delta ))}{\mathrm{\Delta }+(\mathrm{\Delta }-2(1-\delta ))}}=d\cdot {\displaystyle \frac{1-\delta }{\mathrm{\Delta }-1+\delta }}\le (1-\delta ){\displaystyle \frac{d}{\mathrm{\Delta }-1}}.$

As this is at most $1-\delta$ for all , we have up-to-$\mathrm{\Delta }$ uniqueness with gap $\delta$. ∎

A.2 Spectral independence bounds for constant-size graphs

In this section, we prove spectral independence bounds for graphs with fewer than $O(c/\alpha )$-many vertices, since for graphs with such few vertices, our bounds based on contraction of the tree recursions become trivial.

Proof of 28.

If $R_v$ denotes the marginal ratio of a vertex $v\in G$, then $R_v\ge \lambda {\beta }^{\mathrm{\Delta }}$. In the case $\gamma \le 1$, we have $R_v\le \lambda /{\gamma }^{\mathrm{\Delta }}$ as well; if $\gamma >1$, we have $R_v\le \lambda$. It follows that we immediately have the bounds

for all $u,v\in G$. Note that these constants are less than $1$, and only depend on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$, yielding the first claim.

Now, we proceed to remove the dependence on $\mathrm{\Delta }$ when up-to-$\mathrm{\Delta }$ uniqueness holds. We have the following cases:

1. If $\gamma >1$, we immediately obtain a bound of $\frac{\lambda }{1+\lambda }$ which is independent of $\mathrm{\Delta }$.

2. If $\beta =0$ and $\gamma \le 1$, then $\frac{\lambda (1-{\beta }^{\mathrm{\Delta }}{\gamma }^{\mathrm{\Delta }})}{(\lambda +{\gamma }^{\mathrm{\Delta }})(1+\lambda {\beta }^{\mathrm{\Delta }})}=\frac{\lambda }{\lambda +{\gamma }^{\mathrm{\Delta }}}\le \frac{\lambda }{{\gamma }^{\mathrm{\Delta }}}$. Since $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, we must have . It follows that $\frac{\lambda }{{\gamma }^{\mathrm{\Delta }}}\le O(1/\mathrm{\Delta })$.

3. If $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma \le 1$, then

   ${\displaystyle \frac{\lambda (1-{\beta }^{\mathrm{\Delta }}{\gamma }^{\mathrm{\Delta }})}{(\lambda +{\gamma }^{\mathrm{\Delta }})(1+\lambda {\beta }^{\mathrm{\Delta }})}}\le 1-{\beta }^{\mathrm{\Delta }}{\gamma }^{\mathrm{\Delta }}\approx 1-e^{-2}.$

4. If $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, then let ${\mathrm{\Delta }}_0$ be the maximal such that $\sqrt{\beta \gamma }>\frac{d-2}{d}$. If $\lambda \le {\lambda }_c(\beta ,\gamma ,\mathrm{\Delta })$, then by 35, we have

   ${\displaystyle \frac{\lambda (1-{\beta }^{\mathrm{\Delta }}{\gamma }^{\mathrm{\Delta }})}{(\lambda +{\gamma }^{\mathrm{\Delta }})(1+\lambda {\beta }^{\mathrm{\Delta }})}}\le {\displaystyle \frac{\lambda }{{\gamma }^{\mathrm{\Delta }}}}\le O({\mathrm{\Delta }}_0/\mathrm{\Delta }).$

   If $\lambda \ge {\overline{\lambda }}_c(\beta ,\gamma ,\mathrm{\Delta })$, then again by 35, we have

   $$\frac{\lambda (1-{\beta }^{\mathrm{\Delta }}{\gamma }^{\mathrm{\Delta }})}{(\lambda +{\gamma }^{\mathrm{\Delta }})(1+\lambda {\beta }^{\mathrm{\Delta }})}\le \frac{1}{\lambda {\beta }^{\mathrm{\Delta }}}\le O({\mathrm{\Delta }}_0/\mathrm{\Delta }).\mathit{∎}$$

Proof of 29.

The proof is identical to the antiferromagnetic case and we omit it here. ∎

Appendix B Proof of 12 (Parts 1 and 2)

Proof of 12 (Parts 1 and 2).

To see the first equality, we compute directly and get

$\left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\log Z_G^{{\sigma }_{\mathrm{\Lambda }}}$ $={\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}\cdot \left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)Z_G^{{\sigma }_{\mathrm{\Lambda }}}$

$={\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}{\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}\left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\left({\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\displaystyle \prod _{w\in V}}{\lambda }_w^{{\sigma }_w}\right)$

$={\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}{\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}{\sigma }_v\left({\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\displaystyle \prod _{w\in V}}{\lambda }_w^{{\sigma }_w}\right)$

$={\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}{\sigma }_v\cdot {\mu }_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})=M_G^{{\sigma }_{\mathrm{\Lambda }}}(v).$

For Part 2, using the result above, we can also get

$\left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)\log Z_G^{{\sigma }_{\mathrm{\Lambda }}}$

$=$ $\left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\left({\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}\cdot \left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)Z_G^{{\sigma }_{\mathrm{\Lambda }}}\right)$

$=$ ${\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}\cdot \left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)Z_G^{{\sigma }_{\mathrm{\Lambda }}}-{\displaystyle \frac{1}{{(Z_G^{{\sigma }_{\mathrm{\Lambda }}})}^2}}\cdot \left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)Z_G^{{\sigma }_{\mathrm{\Lambda }}}\cdot \left({\lambda }_u{\displaystyle \frac{\partial }{\partial {\lambda }_u}}\right)Z_G^{{\sigma }_{\mathrm{\Lambda }}}$

$=$ ${\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}\cdot \left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\left({\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}{\sigma }_u\left({\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\displaystyle \prod _{w\in V}}{\lambda }_w^{{\sigma }_w}\right)\right)-M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)\cdot M_G^{{\sigma }_{\mathrm{\Lambda }}}(v)$

$=$ ${\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}{\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}{\sigma }_u\cdot \left({\lambda }_v{\displaystyle \frac{\partial }{\partial {\lambda }_v}}\right)\left({\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\displaystyle \prod _{w\in V}}{\lambda }_w^{{\sigma }_w}\right)-M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)\cdot M_G^{{\sigma }_{\mathrm{\Lambda }}}(v)$

$=$ ${\displaystyle \frac{1}{Z_G^{{\sigma }_{\mathrm{\Lambda }}}}}{\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}{\sigma }_u\cdot {\sigma }_v\left({\beta }^{m_1(\sigma )}{\gamma }^{m_0(\sigma )}{\displaystyle \prod _{w\in V}}{\lambda }_w^{{\sigma }_w}\right)-M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)\cdot M_G^{{\sigma }_{\mathrm{\Lambda }}}(v)$

$=$ ${\displaystyle \sum _{\sigma \in {\{0,1\}}^{V\backslash \mathrm{\Lambda }}}}{\sigma }_u\cdot {\sigma }_v\cdot {\mu }_G(\sigma \mid {\sigma }_{\mathrm{\Lambda }})-M_G^{{\sigma }_{\mathrm{\Lambda }}}(u)\cdot M_G^{{\sigma }_{\mathrm{\Lambda }}}(v)$

$=$ $K_G^{{\sigma }_{\mathrm{\Lambda }}}(u,v).\mathit{∎}$

Appendix C A technical lemma for $\mathrm{\Psi }$

The following lemma implies that the potential $\mathrm{\Psi }$ given by Eq. 4 is well-defined.

Lemma 32.

For all $\beta ,\gamma >0$ such that , we have

Proof.

For the $+\mathrm{\infty }$ side we have

Similarly, for the $-\mathrm{\infty }$ side we have

Appendix D Mixing by the potential method: Proof of 24

In this section, we prove 24 in the same way of 5, as outlined in Section 3. The major difference here is that we consider a weighted sum of absolute influences ${\sum }_{v\in V\backslash \mathrm{\Lambda }}$ρv⋅|ℐGσΛ​(r​ → ​v)| where $\rho :V\to {ℝ}^{+}$ is a weight function. This is sufficient for us to bound the eigenvalue of the influence matrix, as indicated by 22. We will choose the weight of a vertex $v$ to be ${\rho }_v={\mathrm{\Delta }}_v$, the degree of $v$. The following lemma provides us an upper bound on the weighted sum of absolute influences to distance $k$, given a general $(\alpha ,c)$-potential. In particular, it generalizes 9.

Lemma 33.

If there exists a general $(\alpha ,c)$-potential function $\mathrm{\Psi }$ with respect to $\mathrm{\Delta }$ and $(\beta ,\gamma ,\lambda )$ where $\alpha \in (0,1)$ and $c>0$, then for every $\mathrm{\Lambda }\subseteq V_T\backslash \{r\}$, ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ and all integers $k\ge 1$,

$$\sum _{v\in L_r(k)}{\mathrm{\Delta }}_v\cdot \left|{ℐ}_T^{{\sigma }_{\mathrm{\Lambda }}}(r\to v)\right|\le 2c\cdot {(1-\alpha )}^{k-1}\cdot {\mathrm{\Delta }}_r$$

where $L_r(k)$ denote the set of all free vertices at distance $k$ away from $r$.

To prove 33, we first state the following generalization of 14 for any weight function $\rho$. The proof of 34 is identical to 14 and we omit here.

Lemma 34.

Let $\mathrm{\Psi }:[-\mathrm{\infty },+\mathrm{\infty }]\to (-\mathrm{\infty },+\mathrm{\infty })$ be a differentiable and increasing (potential) function with image $S=\mathrm{\Psi }[-\mathrm{\infty },+\mathrm{\infty }]$ and derivative $\psi ={\mathrm{\Psi }}^{\prime }$. Denote the degree of the root $r$ by ${\mathrm{\Delta }}_r$. Then for every integer $k\ge 1$,

where

$$A_{\mathrm{\Psi }}=\underset{u\in L_r(1)}{\max }\left\{\frac{|h(\log R_u)|}{\psi (\log R_u)}\right\}\text{and}{B}_{\mathrm{\Psi }}^{\rho }=\underset{v\in L_r(k)}{\max }\left\{{\rho }_v\cdot \psi (\log R_v)\right\}.$$

We then prove 33 and 24.

Proof of 33.

Denote the degree of a vertex $v\in V_T\backslash \{r\}$ by ${\mathrm{\Delta }}_v$, and the degree of $v$ in the subtree $T_v$ by $d_v={\mathrm{\Delta }}_v-1$. Pick the weights of vertices to be ${\rho }_v={\mathrm{\Delta }}_v$ for all $v\in V_T$. Since $\mathrm{\Psi }$ is a general $(\alpha ,c)$-potential, the Contraction condition implies that

Since $\log R_v\in J_{d_v}$ by the definition of $J_d$, the General Boundedness condition implies that for all $u\in L_r(1)$ and $v\in L_r(k)$,

$$\frac{\psi (\log R_v)}{\psi (\log R_u)}\cdot |h(\log R_u)|\le \frac{2c}{{\mathrm{\Delta }}_u+{\mathrm{\Delta }}_v}.$$

Therefore, we get

$${\mathrm{\Delta }}_r{A}_{\mathrm{\Psi }}{B}_{\mathrm{\Psi }}^{\rho }={\mathrm{\Delta }}_r\cdot \underset{u\in L_r(1)}{\max }\left\{\frac{|h(\log R_u)|}{\psi (\log R_u)}\right\}\cdot \underset{v\in L_r(k)}{\max }\left\{{\mathrm{\Delta }}_v\cdot \psi (\log R_v)\right\}\le 2c\cdot {\mathrm{\Delta }}_r.$$

The lemma then follows immediately from 34. ∎

Proof of 24.

The proof of 24 is almost identical to 5. We point out that the only difference here is that we consider the weighted sum of absolute influences of a given vertex. Since the SAW tree preserve degrees of vertices, we can still apply 8. Then, combining 7, 22, 8 and 33, we complete the proof of the theorem. ∎

Appendix E Verifying a good potential: Boundedness

In this subsection, we show the Boundedness or General Boundedness condition for our potential function $\mathrm{\Psi }$ defined by Eq. 4 in different ranges of parameters. Combining 19, we complete the proofs of 10 and 25.

In Section E.1 we give background on the uniqueness region of the parameters $(\beta ,\gamma ,\lambda )$, based on the work of [LLY13]. We then show Boundedness and General Boundedness in Section E.2. Proofs of technical lemmas are left to Section E.3.

E.1 Preliminaries of the uniqueness region

In this section we give a brief description of the uniqueness region of parameters $(\beta ,\gamma ,\lambda )$. All the results here, and also their proofs, can be found in Lemma 21 from the latest version of [LLY13].

Let $\mathrm{\Delta }\ge 3$ be an integer and $\beta ,\gamma ,\lambda$ be reals. We assume that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. For $1\le d\le \mathrm{\Delta }$ define

$$f_d(R)=\lambda {\left(\frac{\beta R+1}{R+\gamma }\right)}^d$$

and denote the unique fixed point of $f_d$ by $R_d^{\ast }$. Recall that the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$ if for all .

When $\beta =0$, the spin system is called a hard-constraint model. In this case, there exists a critical threshold for the external field defined as

such that the parameters $(0,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique if and only if . In particular, when $\gamma \le 1$ the critical field is given by

$${\lambda }_c={\lambda }_c(\gamma ,\mathrm{\Delta })=\frac{{\gamma }^{\mathrm{\Delta }}{(\mathrm{\Delta }-1)}^{\mathrm{\Delta }-1}}{{(\mathrm{\Delta }-2)}^{\mathrm{\Delta }}}.$$

When $\beta >0$, the spin system is called a soft-constraint model. If $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, then $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique for all $\lambda >0$. If $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ the uniqueness region is more complicated which we now describe. Let

$$\overline{\mathrm{\Delta }}=\frac{1+\sqrt{\beta \gamma }}{1-\sqrt{\beta \gamma }},$$

so that for every we have , and for every $d\ge \overline{\mathrm{\Delta }}$ we have $d\cdot \frac{1-\sqrt{\beta \gamma }}{1+\sqrt{\beta \gamma }}\ge 1$. For every , we define $x_1(d)\le x_2(d)$ to be the two positive roots of the quadratic equation

$$\frac{d(1-\beta \gamma )x}{(\beta x+1)(x+\gamma )}=1.$$

More specifically, $x_1(d)$ and $x_2(d)$ are given by

$$x_1(d)=\frac{\theta (d)-\sqrt{\theta {(d)}^2-4\beta \gamma }}{2\beta }\text{and}{x}_2(d)=\frac{\theta (d)+\sqrt{\theta {(d)}^2-4\beta \gamma }}{2\beta }$$

where

$$\theta (d)=d(1-\beta \gamma )-(1+\beta \gamma ).$$

Notice that $\theta (d)\ge 2\sqrt{\beta \gamma }$ for all $d\ge \overline{\mathrm{\Delta }}$. For $i=1,2$ we let

$${\lambda }_i(d)=x_i(d){\left(\frac{x_i(d)+\gamma }{\beta x_i(d)+1}\right)}^d.$$

Then, the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique if and only if $\lambda$ belongs to the following regime

(10)

In particular, when $\gamma \le 1$ there are two critical thresholds such that the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique if and only if or $\lambda >{\overline{\lambda }}_c$ (i.e., $𝒜=(0,{\lambda }_c)\cup ({\overline{\lambda }}_c,\mathrm{\infty })$), where

The following bounds on the critical fields are helpful for our proofs later.

Lemma 35.

1. If $\beta =0$, then for every integer $d$ such that we have

   $${\lambda }_c\le \frac{4{\gamma }^{d+1}}{d-1}.$$

2. If $\beta >0$ and $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, then for every integer $d$ such that we have

   $${\lambda }_1(d)\le \frac{18{\gamma }^{d+1}}{\theta (d)}\text{and}{\lambda }_2(d)\ge \frac{\theta (d)}{18{\beta }^{d+1}}$$

   where $\theta (d)=d(1-\beta \gamma )-(1+\beta \gamma )$.

The proof of 35 is postponed to Section E.3.

E.2 Proofs of boundedness

In this section we complete the proofs of 10 and 25 by establishing Boundedness and General Boundedness in the corresponding range of parameters.

Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. Recall that the potential function $\mathrm{\Psi }$ is defined by

$${\mathrm{\Psi }}^{\prime }(y)=\psi (y)=\sqrt{\frac{(1-\beta \gamma )e^y}{(\beta e^y+1)(e^y+\gamma )}}=\sqrt{\left|h(y)\right|},\mathrm{\Psi }(0)=0.$$ (1)

It is surprising to find out that $\psi =\sqrt{|h|}$, as the potential $\mathrm{\Psi }$ is exactly the one from [LLY13] as indicated by 17. This seems not to be a coincidence, and it provides some intuition why the potential from [LLY13] works. More importantly, the fact that $\psi =\sqrt{|h|}$ is helpful in our proof of Boundedness and General Boundedness. Recall that for and we let $J_d=[\log (\lambda {\beta }^d),\log (\lambda /{\gamma }^d)]$ to be the range of log marginal ratios of a vertex with $d$ children. Then for every and $y_i\in J_{d_i}$ where $i=1,2$, we have

$$\frac{\psi (y_2)}{\psi (y_1)}\cdot |h(y_1)|=\sqrt{|h(y_1)|\cdot |h(y_2)|}.$$ (11)

The following lemma gives upper bounds on $\sqrt{|h(y_1)|\cdot |h(y_2)|}$, from which and Eq. 11 we deduce Boundedness and General Boundedness immediately. The brackets in the lemma indicate which lemma the bound is applied to.

Lemma 36.

Let $\mathrm{\Delta }\ge 3$ be an integer. Let $\beta ,\gamma ,\lambda$ be reals such that $0\le \beta \le \gamma$, $\gamma >0$, and $\lambda >0$. Assume that the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique with gap $\delta \in (0,1)$. Then for all integers $d_1,d_2$ such that , and all reals $y_i\in J_{d_i}$ where $i=1,2$, the following holds:

1. Hard-constraint models: $\beta =0$ and .

   1. (10) If $\gamma \le 1$, then

      $$|h(y_1)|\le \frac{4}{\mathrm{\Delta }}.$$

   2. (25) If $\gamma >1$, then

      $$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{8}{d_1+d_2+2}.$$

2. Soft-constraint models: $\beta >0$ and $\lambda \in 𝒜$.

   1. (10) If $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, then

      $$|h(y_1)|\le \frac{1.5}{\mathrm{\Delta }}.$$

   2. (10) If $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma \le 1$, then

      $$|h(y_1)|\le \frac{18}{\mathrm{\Delta }}.$$

   3. (25) If $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$, then

      $$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{36}{d_1+d_2+2}.$$

The following lemma, whose proof can be found in Section E.3, is helpful.

Lemma 37.

For every $y\in [-\mathrm{\infty },+\mathrm{\infty }]$ we have

$$|h(y)|=\frac{|1-\beta \gamma |e^y}{(\beta e^y+1)(e^y+\gamma )}\le \frac{|1-\sqrt{\beta \gamma }|}{1+\sqrt{\beta \gamma }}.$$

We present here the proof of 36.

Proof of 36.

We use notations and results from Section E.1.

H. Hard-constraint models: $\beta =0$ and .

H.1. $\gamma \le 1$.

For every $y_1\in J_{d_1}$ we deduce from 35 that

$$e^{y_1}\le \frac{\lambda }{{\gamma }^{d_1}}\le \frac{{\lambda }_c}{{\gamma }^{\mathrm{\Delta }-1}}\le \frac{4\gamma }{\mathrm{\Delta }-2}.$$

Hence,

$$|h(y_1)|=\frac{e^{y_1}}{e^{y_1}+\gamma }\le \frac{\frac{4\gamma }{\mathrm{\Delta }-2}}{\frac{4\gamma }{\mathrm{\Delta }-2}+\gamma }=\frac{4}{\mathrm{\Delta }+2}\le \frac{4}{\mathrm{\Delta }}.$$

H.2. $\gamma >1$.

Let $\overline{y}=\frac{y_1+y_2}{2}$ and $\overline{d}=\frac{d_1+d_2}{2}$. Then we get

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}=\sqrt{\frac{e^{y_1}}{e^{y_1}+\gamma }}\cdot \sqrt{\frac{e^{y_2}}{e^{y_2}+\gamma }}=\frac{1}{\sqrt{(1+\gamma e^{-y_1})(1+\gamma e^{-y_2})}}\le \frac{1}{1+\gamma e^{-\overline{y}}},$$

where the last inequality follows from the AM–GM inequality by

$$(1+\gamma e^{-y_1})(1+\gamma e^{-y_2})=1+\gamma (e^{-y_1}+e^{-y_2})+{\gamma }^2{e}^{-2\overline{y}}\ge 1+2\gamma e^{-\overline{y}}+{\gamma }^2{e}^{-2\overline{y}}={(1+\gamma e^{-\overline{y}})}^2.$$

Since $y_i\in J_{d_i}$ for $i=1,2$, we have

$$e^{\overline{y}}=\sqrt{e^{y_1}\cdot e^{y_2}}\le \sqrt{\frac{\lambda }{{\gamma }^{d_1}}\cdot \frac{\lambda }{{\gamma }^{d_2}}}=\frac{\lambda }{{\gamma }^{\overline{d}}}.$$

If $\overline{d}\ge 2$, then we deduce from 35 and $\gamma >1$ that

$$e^{\overline{y}}\le \frac{{\lambda }_c}{{\gamma }^{\lfloor \overline{d}\rfloor }}\le \frac{4\gamma }{\lfloor \overline{d}\rfloor -1}.$$

It follows that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1}{1+\gamma e^{-\overline{y}}}\le \frac{1}{1+\frac{\lfloor \overline{d}\rfloor -1}{4}}=\frac{4}{\lfloor \overline{d}\rfloor +3}\le \frac{8}{d_1+d_2+2}.$$

If , then it is easy to see that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le 1\le \frac{8}{d_1+d_2+2}.$$

S. Soft-constraint models: $\beta >0$ and $\lambda \in 𝒜$.

S.1. $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$.

For every $y_1\in J$ we deduce from 37 that

$$|h(y_1)|\le \frac{1-\sqrt{\beta \gamma }}{1+\sqrt{\beta \gamma }}\le \frac{1}{\mathrm{\Delta }-1}\le \frac{1.5}{\mathrm{\Delta }}.$$

S.2. $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma \le 1$.

In this case, we have either or $\lambda >{\overline{\lambda }}_c$ where ${\lambda }_c,{\overline{\lambda }}_c$ are the two critical fields. Consider first $\lambda >{\overline{\lambda }}_c$. For every $y_1\in J_{d_1}$ we deduce from 35 and that

$$e^{y_1}\ge \lambda {\beta }^{d_1}\ge {\overline{\lambda }}_c{\beta }^{\mathrm{\Delta }-1}\ge \frac{\theta (\mathrm{\Delta }-1)}{18\beta }$$

where $\theta (d)=d(1-\beta \gamma )-(1+\beta \gamma )$. Hence,

$|h(y_1)|={\displaystyle \frac{(1-\beta \gamma )e^{y_1}}{(\beta e^{y_1}+1)(e^{y_1}+\gamma )}}$ $={\displaystyle \frac{1-\beta \gamma }{\beta e^{y_1}+\gamma e^{-y_1}+(1+\beta \gamma )}}$

$\le {\displaystyle \frac{1-\beta \gamma }{\frac{\theta (\mathrm{\Delta }-1)}{18}+(1+\beta \gamma )}}={\displaystyle \frac{18(1-\beta \gamma )}{(\mathrm{\Delta }-1)(1-\beta \gamma )+17(1+\beta \gamma )}}\le {\displaystyle \frac{18}{\mathrm{\Delta }}}.$

Next we consider . For every $y_1\in J_{d_1}$ we deduce from 35 and $\gamma \le 1$ that

$$e^{y_1}\le \frac{\lambda }{{\gamma }^{d_1}}\le \frac{{\lambda }_c}{{\gamma }^{\mathrm{\Delta }-1}}\le \frac{18\gamma }{\theta (\mathrm{\Delta }-1)}.$$

Hence,

$$|h(y_1)|=\frac{1-\beta \gamma }{\beta e^{y_1}+\gamma e^{-y_1}+(1+\beta \gamma )}\le \frac{1-\beta \gamma }{\frac{\theta (\mathrm{\Delta }-1)}{18}+(1+\beta \gamma )}\le \frac{18}{\mathrm{\Delta }}.$$

S.3. $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ and $\gamma >1$.

Let $\overline{y}=\frac{y_1+y_2}{2}$, $\overline{d}=\frac{d_1+d_2}{2}$, $d_L=\lfloor \overline{d}\rfloor$, and $d_R=\lceil \overline{d}\rceil$. We first consider some trivial cases. If $\overline{d}\le 2$ then it is easy to see that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le 1\le \frac{6}{d_1+d_2+2}.$$

If $\overline{d}>2$ and $d_L\le \overline{\mathrm{\Delta }}$, then we deduce from 37 that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\sqrt{\beta \gamma }}{1+\sqrt{\beta \gamma }}=\frac{1}{\overline{\mathrm{\Delta }}}\le \frac{2}{d_1+d_2-2}\le \frac{6}{d_1+d_2+2}.$$

Hence, in the following we may assume that $\overline{d}>2$ and $d_L>\overline{\mathrm{\Delta }}$.

Since the parameters $(\beta ,\gamma ,\lambda )$ are up-to-$\mathrm{\Delta }$ unique, we have $\lambda \in 𝒜$ where the regime $𝒜$ is given by Eq. 10. Observe that

$$𝒜\subseteq (0,{\lambda }_1(d_L))\cup ({\lambda }_2(d_R),\mathrm{\infty })\cup ({\lambda }_2(d_L),{\lambda }_1(d_R))$$

where the last interval is nonempty only when . This means that $\lambda$ is contained in at least one of the three intervals. We establish the bound by considering these three cases separately.

Case 1: . By the Cauchy-Schwarz inequality, we have

$\sqrt{|h(y_1)|\cdot |h(y_2)|}$ $=\sqrt{{\displaystyle \frac{1-\beta \gamma }{\beta e^{y_1}+\gamma e^{-y_1}+(1+\beta \gamma )}}}\cdot \sqrt{{\displaystyle \frac{1-\beta \gamma }{\beta e^{y_2}+\gamma e^{-y_2}+(1+\beta \gamma )}}}$

$\le {\displaystyle \frac{1-\beta \gamma }{\sqrt{(\beta e^{y_1}+\gamma e^{-y_1})(\beta e^{y_2}+\gamma e^{-y_2})}+(1+\beta \gamma )}}.$ (12)

Therefore, we get

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\beta \gamma }{\gamma e^{-\overline{y}}+(1+\beta \gamma )}.$$

Since $y_i\in J_{d_i}$ for $i=1,2$ and $\gamma >1$, we deduce from 35 that

$$e^{\overline{y}}\le \frac{\lambda }{{\gamma }^{\overline{d}}}\le \frac{{\lambda }_1(d_L)}{{\gamma }^{d_L}}\le \frac{18\gamma }{\theta (d_L)},$$

where $\theta (d_L)=d_L(1-\beta \gamma )-(1+\beta \gamma )$. It follows that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\beta \gamma }{\gamma e^{-\overline{y}}+(1+\beta \gamma )}\le \frac{1-\beta \gamma }{\frac{\theta (d_L)}{18}+(1+\beta \gamma )}\le \frac{36}{d_1+d_2+2}.$$

Case 2: $\lambda >{\lambda }_2(d_R)$. Similarly, we obtain from Eq. 12 that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\beta \gamma }{\beta e^{\overline{y}}+(1+\beta \gamma )}.$$

Since $y_i\in J_{d_i}$ for $i=1,2$ and , we deduce from 35 that

$$e^{\overline{y}}\ge \lambda {\beta }^{\overline{d}}\ge {\lambda }_2(d_R){\beta }^{d_R}\ge \frac{\theta (d_R)}{18\beta },$$

where $\theta (d_R)=d_R(1-\beta \gamma )-(1+\beta \gamma )$. It follows that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\beta \gamma }{\beta e^{\overline{y}}+(1+\beta \gamma )}\le \frac{1-\beta \gamma }{\frac{\theta (d_R)}{18}+(1+\beta \gamma )}\le \frac{36}{d_1+d_2+2}.$$

Case 3: . We may assume that $d_1\ge d_2$. By Eq. 12, we obtain

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\beta \gamma }{\sqrt{\beta \gamma }{e}^{\frac{y_2-y_1}{2}}+(1+\beta \gamma )}.$$

Since $y_i\in J_{d_i}$ for $i=1,2$ and , we have

$$e^{y_2-y_1}\ge {\beta }^{d_2}{\gamma }^{d_1}\ge {\beta }^{d_L}{\gamma }^{d_R}.$$

Meanwhile, we deduce from 35 that

which implies

$$\sqrt{\beta \gamma }{e}^{\frac{y_2-y_1}{2}}\ge \sqrt{{\beta }^{d_L+1}{\gamma }^{d_R+1}}\ge \frac{\sqrt{\theta (d_L)\theta (d_R)}}{18}\ge \frac{\theta (d_L)}{18}.$$

It follows that

$$\sqrt{|h(y_1)|\cdot |h(y_2)|}\le \frac{1-\beta \gamma }{\sqrt{\beta \gamma }{e}^{\frac{y_2-y_1}{2}}+(1+\beta \gamma )}\le \frac{1-\beta \gamma }{\frac{\theta (d_L)}{18}+(1+\beta \gamma )}\le \frac{36}{d_1+d_2+2}.\mathit{∎}$$

E.3 Proofs of technical lemmas

Proof of 35.

1. For every we have

$${\lambda }_c\le \frac{{\gamma }^{d+1}{d}^d}{{(d-1)}^{d+1}}=\frac{{\gamma }^{d+1}}{d-1}{\left(\frac{d}{d-1}\right)}^d\le \frac{4{\gamma }^{d+1}}{d-1},$$

where the last inequality follows from that ${(\frac{d}{d-1})}^d\le 4$ for all integer $d>1$.

2. For every we have

$$x_1(d)=\frac{2\gamma }{\theta (d)+\sqrt{\theta {(d)}^2-4\beta \gamma }}\le \frac{2\gamma }{\theta (d)}.$$

Observe that the function $\frac{x+\gamma }{\beta x+1}$ is monotone increasing in $x$ when , and thus we deduce that

$$\frac{x_1(d)+\gamma }{\beta x_1(d)+1}\le \frac{\frac{2\gamma }{\theta (d)}+\gamma }{\frac{2\beta \gamma }{\theta (d)}+1}=\gamma \cdot \frac{2+d(1-\beta \gamma )-(1+\beta \gamma )}{2\beta \gamma +d(1-\beta \gamma )-(1+\beta \gamma )}=\gamma \cdot \frac{d+1}{d-1}.$$

Therefore,

$${\lambda }_1(d)=x_1(d){\left(\frac{x_1(d)+\gamma }{\beta x_1(d)+1}\right)}^d\le \frac{2\gamma }{\theta (d)}\cdot {\gamma }^d\cdot {\left(\frac{d+1}{d-1}\right)}^d\le \frac{18{\gamma }^{d+1}}{\theta (d)}$$

where the last inequality follows from that ${(\frac{d+1}{d-1})}^d\le 9$ for all integer $d>1$.

The second part can be proved similarly. For every we have

$$x_2(d)=\frac{\theta (d)+\sqrt{\theta {(d)}^2-4\beta \gamma }}{2\beta }\ge \frac{\theta (d)}{2\beta },$$

and hence,

$$\frac{x_2(d)+\gamma }{\beta x_2(d)+1}\ge \frac{\frac{\theta (d)}{2\beta }+\gamma }{\frac{\theta (d)}{2}+1}=\frac{1}{\beta }\cdot \frac{d(1-\beta \gamma )-(1+\beta \gamma )+2\beta \gamma }{d(1-\beta \gamma )-(1+\beta \gamma )+2}=\frac{1}{\beta }\cdot \frac{d-1}{d+1}.$$

We then conclude that

$${\lambda }_2(d)=x_2(d){\left(\frac{x_2(d)+\gamma }{\beta x_2(d)+1}\right)}^d\ge \frac{\theta (d)}{2\beta }\cdot \frac{1}{{\beta }^d}\cdot {\left(\frac{d-1}{d+1}\right)}^d\ge \frac{\theta (d)}{18{\beta }^{d+1}},$$

where the last inequality again follows from that ${(\frac{d+1}{d-1})}^d\le 9$ for all integer $d>1$. ∎

Proof of 37.

We deduce from the AM–GM inequality that

$$|h(y)|=\frac{|1-\beta \gamma |}{\beta e^y+\gamma e^{-y}+1+\beta }\le \frac{|1-\beta \gamma |}{2\sqrt{\beta \gamma }+1+\beta }=\frac{|1-\sqrt{\beta \gamma }|}{1+\sqrt{\beta \gamma }}.\mathit{∎}$$

Appendix F Proofs for ferromagnetic cases

F.1 Proof of 26

Proof of 26.

Throughout, we use the “trivial potential” function $\mathrm{\Psi }(y)=y$. Note that then, $\psi (y)=1$ is a constant function. Now, we prove Contraction and Boundedness. We split into the three cases.

1. We first prove the Contraction part. By 37, for all $y\in [-\mathrm{\infty },+\mathrm{\infty }]$ we have

   $$\left|h(y)\right|\le \frac{|1-\sqrt{\beta \gamma }|}{1+\sqrt{\beta \gamma }}\le \frac{1-\delta }{\mathrm{\Delta }-1}.$$

   Now let us prove the Boundedness condition. From the above inequality we have

   $$\left|h(y)\right|\le \frac{1}{\mathrm{\Delta }-1}\le \frac{1.5}{\mathrm{\Delta }}$$

   for $\mathrm{\Delta }\ge 3$.

2. For the Contraction part, since $\log (\lambda \max \{1,1/{\gamma }^{\mathrm{\Delta }-1}\})\le y_i\le \log (\lambda \max \{1,{\beta }^{\mathrm{\Delta }-1}\})$, we have

   $\left|{\displaystyle \frac{\partial H_d(𝒚)}{\partial y_i}}\right|$ $=\left|h(y_i)\right|={\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\gamma e^{-y_i}+\beta e^{y_i}}}\le {\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\gamma e^{-y_i}}}$

   $\le {\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\frac{\gamma }{\lambda \max \{1,{\beta }^{\mathrm{\Delta }-1}\}}}}.$

   Since we assumed $\lambda \le (1-\delta )\frac{\gamma }{\max \{1,{\beta }^{\mathrm{\Delta }-1}\}\cdot ((\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta })}$, it follows that we have the upper bound

   ${\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{1-\delta }}}$ $=(1-\delta ){\displaystyle \frac{\beta \gamma -1}{(\mathrm{\Delta }-1-\delta )\beta \gamma -(\mathrm{\Delta }-1+\delta )}}$

   $=(1-\delta ){\displaystyle \frac{\beta \gamma -1}{(\mathrm{\Delta }-1-\delta )(\beta \gamma -1)+2\delta }}$

   $\le {\displaystyle \frac{1-\delta }{\mathrm{\Delta }-1-\delta }}\le (1-\mathrm{\Theta }(\delta )){\displaystyle \frac{1}{\mathrm{\Delta }-1}}.$

   Now, we prove the Boundedness condition. Note that since $\lambda \le \frac{\gamma }{\max \{1,{\beta }^{\mathrm{\Delta }-1}\}\cdot ((\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta })}$, it follows that $y\le \log (\lambda \max \{1,{\beta }^{\mathrm{\Delta }-1}\})\le \log \left(\frac{\gamma }{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}\right)$. A simple calculation reveals that $\frac{\gamma }{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}\le \sqrt{\frac{\gamma }{\beta }}$ and so by 37, we have

   $\left|h(y)\right|$ $\le \left|h\left(\log \left({\displaystyle \frac{\gamma }{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}}\right)\right)\right|\le {\displaystyle \frac{(\beta \gamma -1)e^{\log \left(\frac{\gamma }{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}\right)}}{e^{\log \left(\frac{\gamma }{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}\right)}+\gamma }}$

   $=(\beta \gamma -1){\displaystyle \frac{1}{1+(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}}={\displaystyle \frac{\beta \gamma -1}{(\mathrm{\Delta }-2)(\beta \gamma -1)-1}}\le O(1/\mathrm{\Delta }).$

3. For the Contraction part, since $\log (\lambda \max \{1,1/{\gamma }^{\mathrm{\Delta }-1}\})\le y_i\le \log (\lambda \max \{1,{\beta }^{\mathrm{\Delta }-1}\})$, we have

   $\left|{\displaystyle \frac{\partial H_d(𝒚)}{\partial y_i}}\right|$ $=\left|h(y_i)\right|={\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\gamma e^{-y_i}+\beta e^{y_i}}}\le {\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\beta e^{y_i}}}$

   $\le {\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\beta \lambda \max \{1,1/{\gamma }^{\mathrm{\Delta }-1}\}}}.$

   Since we assumed $\lambda \ge \frac{1}{1-\delta }\cdot \frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta \cdot \min \{1,1/{\gamma }^{\mathrm{\Delta }-1}\}}$, it follows that we have the upper bound

   ${\displaystyle \frac{\beta \gamma -1}{1+\beta \gamma +\frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{1-\delta }}}$

   which is again is upper bounded by $(1-\mathrm{\Theta }(\delta ))\frac{1}{\mathrm{\Delta }-1}$ as we calculated in case 2 above.

   Now, we prove the Boundedness condition. Note that since $\lambda \ge \frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta \min \{1,1/{\gamma }^{\mathrm{\Delta }-2}}$, it follows that $y\ge \log (\lambda \min \{1,1/{\gamma }^{\mathrm{\Delta }-1}\}\ge \log \left(\frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta }\right)$. A simple calculation reveals that $\frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta }\ge \sqrt{\frac{\gamma }{\beta }}$ and so by 37, we have

   $\left|h(y)\right|$ $\le \left|h\left(\log \left({\displaystyle \frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta }}\right)\right)\right|\le (\beta \gamma -1){\displaystyle \frac{1}{\beta \cdot \frac{(\mathrm{\Delta }-2)\beta \gamma -\mathrm{\Delta }}{\beta }+1}}$

   $={\displaystyle \frac{\beta \gamma -1}{(\mathrm{\Delta }-2)(\beta \gamma -1)-1}}\le O(1/\mathrm{\Delta }).\mathit{∎}$

F.2 Proof of 27

In this subsection, we use results from [GL18] to prove 27. Their potential function is implicitly defined by its derivative for the marginal ratios as

${\mathrm{\Phi }}^{\prime }(R)=\varphi (R)=\min \{{\displaystyle \frac{\beta \gamma -1}{\alpha \gamma \log \frac{\lambda +\gamma }{\beta \lambda +1}}},{\displaystyle \frac{1}{R\log \frac{\lambda }{R}}}\}$

for a constant $0\le \alpha \le 1$ depending only on $\beta ,\gamma ,\lambda$ (see [GL18] for a precise definition). In our context, the corresponding potential for the log ratios is

${\mathrm{\Psi }}^{\prime }(y)=\psi (y)=e^y\varphi (e^y)=\min \{{\displaystyle \frac{\beta \gamma -1}{\alpha \gamma \log \frac{\lambda +\gamma }{\beta \lambda +1}}}{e}^y,{\displaystyle \frac{1}{\log \frac{\lambda }{e^y}}}\}$

and is bounded by constants depending on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$ for $\log (\lambda /{\gamma }^{\mathrm{\Delta }-1})\le y\le \log \lambda$.

One of the main technical results in [GL18] is showing that the tree recursion is contracting with the potential function $\mathrm{\Phi }$, and the derivative $\varphi$ is bounded in the sense that there exist positive constants $C_1,C_2$ depending only on $\beta ,\gamma ,\lambda$ such that $C_1\le \varphi (R)\le C_2$ for all $0\le R\le \lambda$. [GL18] refer to such a function as a universal potential function.

In our context, we get that $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential function which satisfies 4, but with a constant $c$ that depends on $\gamma ,\mathrm{\Delta }$. Indeed, worst case, we have

$\underset{y_1,y_2}{\max }{\displaystyle \frac{\psi (y_2)}{\psi (y_1)}}\ge {\displaystyle \frac{\psi (\log \lambda )}{\psi (\log (\lambda /{\gamma }^{\mathrm{\Delta }-1}))}}={\displaystyle \frac{\lambda \frac{\beta \gamma -1}{\alpha \gamma \log \frac{\lambda +\gamma }{\beta \lambda +1}}}{\frac{\beta \gamma -1}{\alpha \log \frac{\lambda +\gamma }{\beta \lambda +1}}\cdot \frac{\lambda }{{\gamma }^{\mathrm{\Delta }}}}}={\gamma }^{\mathrm{\Delta }-1}.$

More precisely, we have the following result from [GL18], stated in terms of the log marginal ratios.

Theorem 38.

Assume $\beta ,\gamma ,\lambda$ are nonnegative real numbers satisfying $\beta \le 1\le \gamma$, $\sqrt{\beta \gamma }\ge 1$, and . Then the function $\mathrm{\Psi }$ is an $(\alpha ,c)$-potential function for a constant depending on $\beta ,\gamma ,\lambda$, and a constant $c>0$ depending on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$.

Combined with 5, this gives $O(n^C)$ mixing with a constant $C$ depending only on $\beta ,\gamma ,\lambda ,\mathrm{\Delta }$. We note this is weaker than the correlation decay result in [GL18], since there, $C$ does not depend on $\mathrm{\Delta }$, and hence is efficient for arbitrary graphs.

Appendix G Slightly faster mixing

In this section, we slightly optimize our mixing time results for certain antiferromagnetic 2-spin systems by more carefully taking into account the tradeoff between the (nontrivial) spectral independence bound we prove based on contraction, and the (trivial) spectral independence bound we obtained in Section A.2 for handling constant-sized graphs.

Proposition 39.

Suppose a distribution $\mu$ on subsets of $[n]$ is $({\eta }_0,\mathrm{\dots },{\eta }_{n-2})$-spectrally independent for ${\eta }_i\le \min \{a,(n-i-1)b\}$, for some $a\ge 0$ and $0\le b\le 1$. Then the Glauber dynamics for sampling from $\mu$ has spectral gap at least $\frac{1}{n}\cdot \mathrm{\Omega }{\left(\frac{a}{bn}\right)}^a$

Proof.

Suppose we have already conditioned on $c$-fraction of elements to be “in/out. The resulting distribution is both $b(1-c)n$-spectrally independent and $a$-spectrally independent. The exact threshold $c$ for which the bound $b(1-c)n$ is better than $a$ is given by

$c=1-{\displaystyle \frac{a}{bn}}$

We note such a $c$ only makes sense when $0\le 1-\frac{a}{bn}\le 1$, or equivalently, $bn\ge a$. Now, we apply the $a$-spectral independence bound for all conditional distributions based on fixing at most $c$-fraction of vertices. We apply the $(n-i-1)b$-spectral independence otherwise. We obtain a final spectral gap lower bound of

${\displaystyle \frac{1}{n}}\cdot {(1-b)}^{(1-c)n}\cdot {\displaystyle \prod _{k=0}^{cn}}\left(1-{\displaystyle \frac{a}{n-k-1}}\right)$

Observe that

${(1-b)}^{(1-c)n}={(1-b)}^{\frac{a}{b}}\gtrsim \exp (-a)$

We also have

${\displaystyle \prod _{k=0}^{cn}}\left(1-{\displaystyle \frac{a}{n-k-1}}\right)$ $\gtrsim \exp \left(-a{\displaystyle \sum _{k=0}^{cn}}{\displaystyle \frac{1}{n-k-1}}\right)$

$\gtrsim \exp \left(-a\left(\underset{\approx \log n}{\underbrace{{\displaystyle \sum _{k=0}^{n-2}}{\displaystyle \frac{1}{n-k-1}}}}-\underset{\approx \log (1-c)n}{\underbrace{{\displaystyle \sum _{k=cn+1}^{n-2}}{\displaystyle \frac{1}{n-k-1}}}}\right)\right)$

$\gtrsim \exp \left(-a\cdot \log {\displaystyle \frac{1}{1-c}}\right)$

$\gtrsim \exp \left(-a\log {\displaystyle \frac{bn}{a}}\right)$

$\gtrsim {\left({\displaystyle \frac{a}{bn}}\right)}^a$

Putting these together, we obtain the desired lower bound. ∎

With this result, we can apply it to the antiferromagnetic models with $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }},\gamma \le 1$ and $\beta =0,\gamma \le 1$, since looking in the proof of 28, we have such systems are $Cn$-spectrally independent roughly with $C\le O(1/\mathrm{\Delta })$.

Corollary 40 (Soft Constraints).

Fix integers $\mathrm{\Delta }\ge 3$, . Let $\beta ,\gamma ,\lambda \ge 0$ be nonnegative real numbers satisfying $\frac{\overline{\mathrm{\Delta }}-2}{\overline{\mathrm{\Delta }}}\le \sqrt{\beta \gamma }\le \frac{\overline{\mathrm{\Delta }}-1}{\overline{\mathrm{\Delta }}+1}$ and $\gamma \le 1$. Assume further that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique with gap . Then for every $n$-vertex graph $G$ with maximum degree at most $\mathrm{\Delta }$, the Glauber dynamics for sampling from the antiferromagnetic 2-spin system with parameters $(\beta ,\gamma ,\lambda )$ mixes in $O{\left(\frac{\overline{\mathrm{\Delta }}\cdot n}{\mathrm{\Delta }}\right)}^{O(1/\delta )}$ steps.

Corollary 41 (Hard Constraints).

Fix an integer $\mathrm{\Delta }\ge 3$, fix $\beta =0$, and let $0\le \gamma \le 1,\lambda \ge 0$ be up-to-$\mathrm{\Delta }$ unique with gap . Then for every $n$-vertex graph $G$ with maximum degree at most $\mathrm{\Delta }$, the Glauber dynamics for sampling from the antiferromagnetic 2-spin system with parameters $(\beta ,\gamma ,\lambda )$-mixes in $O{\left(\frac{n}{\mathrm{\Delta }}\right)}^{O(1/\delta )}$ steps.