Correlation Decay up to Uniqueness in Spin Systems

Abstract

We give a complete characterization of the two-state anti-ferromagnetic spin systems which are of strong spatial mixing on general graphs. We show that a two-state anti-ferromagnetic spin system is of strong spatial mixing on all graphs of maximum degree at most $\mathrm{\Delta }$ if and only if the system has a unique Gibbs measure on infinite regular trees of degree up to $\mathrm{\Delta }$, where $\mathrm{\Delta }$ can be either bounded or unbounded. As a consequence, there exists an FPTAS for the partition function of a two-state anti-ferromagnetic spin system on graphs of maximum degree at most $\mathrm{\Delta }$ when the uniqueness condition is satisfied on infinite regular trees of degree up to $\mathrm{\Delta }$. In particular, an FPTAS exists for arbitrary graphs if the uniqueness is satisfied on all infinite regular trees. This covers as special cases all previous algorithmic results for two-state anti-ferromagnetic systems on general-structure graphs.

Combining with the FPRAS for two-state ferromagnetic spin systems of Jerrum-Sinclair and Goldberg-Jerrum-Paterson, and the very recent hardness results of Sly-Sun and independently of Galanis-S̆tefankovic̆-Vigoda, this gives a complete classification, except at the phase transition boundary, of the approximability of all two-state spin systems, on either degree-bounded families of graphs or family of all graphs.

1 Introduction

Spin systems are well studied in the areas of Statistical Physics, Applied Probability and Computer Science as a general framework to capture the essence of how local interactions and constrains affect the macroscopic properties of particle systems. A system is usually described by a graph, with each vertex in one of a fixed number of states called spins, and edges specifying the neighborhood relation of the system.

Let $G(V,E)$ be a graph and $q$ be the number of spin states. A configuration of the system is one of the $q^{|V|}$ possible assignments $\sigma :V\to [q]$. Each configuration $\sigma$ has an energy $H(\sigma )$ as a sum over all edges and vertices, such that the contribution of an edge $(u,v)\in E$ is determined by a symmetric function of the spin states $\sigma (u)$ and $\sigma (v)$, and the contribution of a vertex $v\in V$ is determined by a function of its spin state $\sigma (v)$. The weight of a configuration $\sigma$ is $w(\sigma )=\exp (-\frac{H(\sigma )}{T})$, where $T$ is the temperature. We focus on the two-state spin systems. Up to normalization, a two-state spin system is fully captured by three parameters $(\beta ,\gamma ,\lambda )$, where $\beta$ and $\gamma$ determine the symmetric function for edge contribution and $\lambda$, also known as the external field, determines the function for vertex contribution. The Gibbs measure is a natural probability distribution over all configurations such that the probability of a configuration $\sigma$ is $\frac{w(\sigma )}{Z}$, where the normalizing factor $Z=\sum w(\sigma )$ is called the partition function. The partition function encodes rich information regarding the macroscopic behavior of the spin system. However, for almost all nontrivial settings it is #P-hard to compute the precise value of partition functions.

One of the most important properties of spin systems is the correlation decay, which says that the correlation between the marginal Gibbs distributions of two vertices decays rapidly with respect to the distance between the two vertices. This property is also called weak spatial mixing [26]. Of greater algorithmic significance is the strong spatial mixing, which says that the correlations decay in the presence of arbitrary fixed spins at other vertices. For two-state spin systems, the strong spatial mixing may imply efficient approximation algorithms for the partition function. It is then an important question to characterize the systems which exhibit strong spatial mixing on arbitrary instances of graphs by the parameters of the systems.

A two-state spin system is ferromagnetic if adjacent vertices favor agreement of spins, and is anti-ferromagnetic if otherwise. For all two-state ferromagnetic spin systems, the partition functions can be efficiently approximated [12, 10]. In the anti-ferromagnetic region, the correlation decay plays a central role in the approximability of partition functions. It is believed (see [19]) that for such models the approximability of the partition function is characterized by the uniqueness of Gibbs measure on infinite regular trees, which is equivalent to the weak spatial mixing on the infinite regular trees. This condition is called the uniqueness condition.

In a seminal work [26], Weitz shows that for the hardcore model strong spatial mixing is characterized by the uniqueness condition, and in a recent work of Sinclair, Srivastava, and Thurley [21] the same characterization is proved for the anti-ferromagnetic Ising model. Both models are important special two-state anti-ferromagnetic spin systems. On the hardness side, it is proved for the hardcore model by Sly [23] and Galanis et al. [6], and very recently for the general two-state anti-ferromagnetic spin systems by Sly and Sun [24] and independently by Galanis, S̆tefankovic̆, and Vigoda [7] that violating the uniqueness condition implies the inapproximability of partition functions. Two questions remain open for our complete understanding of the correlation decay and computation in two-state spin systems: the characterizations of the strong spatial mixing and approximability of general two-state spin systems.

1.1 Our results

We characterize the two-state anti-ferromagnetic spin systems which exhibit strong spatial mixing on general degree-bounded graphs or arbitrary graphs by the uniqueness of Gibbs measure on infinite regular trees.

Theorem 1.

For any finite $\mathrm{\Delta }\ge 2$ or $\mathrm{\Delta }=\mathrm{\infty }$, a two-state anti-ferromagnetic spin system is of strong spatial mixing on graphs of maximum degree at most $\mathrm{\Delta }$ if and only if the system exhibits uniqueness on infinite $d$-regular trees for all $d\le \mathrm{\Delta }$.

Due to Weitz’s self-avoiding walk tree construction [26], the strong spatial mixing of a two-state spin system on degree-bounded graphs immediately implies an FPTAS for the partition function. Indeed, we show an even stronger notion of correlation decay introduced in a previous work [15], namely the computationally efficient correlation decay, which gives FPTAS not only for the degree-bounded graphs but also for arbitrary graphs with unbounded degrees, when the corresponding uniqueness condition is satisfied.

Theorem 2.

For any finite $\mathrm{\Delta }\ge 3$ or $\mathrm{\Delta }=\mathrm{\infty }$, there exists an FPTAS for the partition function of the two-state anti-ferromagnetic spin system on graphs of maximum degree at most $\mathrm{\Delta }$ if for all $d\le \mathrm{\Delta }$ the system parameters lie in the interior of the uniqueness region of infinite $d$-regular tree.

In the above two theorems, the $\mathrm{\Delta }=\mathrm{\infty }$ case represents the graphs of unbounded degree.

Due to a very recent hardness results of Sly and Sun [24] for general two-state spin systems and an independent result of Galanis, S̆tefankovic̆ and Vigoda [7] for a less general setting, violating the uniqueness condition implies inapproximability of the partition function.

Theorem 3 (Due to [24] and [7]).

For any finite $\mathrm{\Delta }\ge 3$ or $\mathrm{\Delta }=\mathrm{\infty }$, unless $NP=RP$ there does not exist an FPRAS for the partition function of the two-state anti-ferromagnetic spin system on graphs of maximum degree at most $\mathrm{\Delta }$ if for some $d\le \mathrm{\Delta }$ the system parameters lie in the interior of the non-uniqueness region of infinite $d$-regular tree.

The original theorem in [24] holds for $d$-regular graphs with fixed $d$, which immediately implies the hardness for degree-bounded graphs or arbitrary graphs. And the hardness condition in both [24] and [7] requires the non-uniqueness on a $d$-regular tree with $d\ge 3$. But the uniqueness on the infinite 2-regular tree (i.e. the infinite path) always holds for any two-state anti-ferromagnetic spin system, thus the condition in Theorem 3 suffices.

For graphs of maximum degree 2 or less, the partition function can be computed exactly in polynomial time. Thus Theorem 2 and Theorem 3 together with the FPRAS for two-state ferromagnetic spin system [12, 10] give an almost complete (except at the phase transition boundary) classification of the approximability of the partition function of all two-state spin systems, on either degree-bounded families of graphs or family of all graphs.

1.1.1 Regularity and monotonicity

In Statistical Physics, the correlation decay is usually studied on regular graphs or even structurally symmetric graphs (e.g. Bethe lattice). Although for hardness results considering only regular graphs will strengthen the lower bounds, from algorithmic perspective it is more general to consider spin systems on general graphs. The approximation algorithms in [26] for the hardcore model and [21] for the anti-ferromagnetic Ising model are both for general graphs with bounded maximum degrees. As discussed in [21, 24], up to translation of parameters the hardcore and Ising models together are complete for general two-state spin systems on $d$-regular graphs with fixed $d$. However, this does not cover the most general case, namely the general two-state spin system on general graphs, of either bounded or unbounded degree. A fundamental reason for this gap is the non-monotonicity of general spin systems.

The well studied hardcore and anti-ferromagnetic Ising models (along with all two-state spin systems with $\beta ,\gamma \le 1$) are both monotone spin systems, in the sense that the uniqueness on infinite $d$-regular tree implies the uniqueness on all infinite regular trees with smaller degrees. This monotonicity does not necessarily hold in general two-state spin systems.

In [26], Weitz established the following implication in the hardcore model:

Claim (Theorem 2.3 in [26]).

Strong spatial mixing on a $d$-regular tree implies the strong spatial mixing on all graphs of maximum degree at most $d$.

In [26], Weitz also remarked without proof that this implication holds for all two-state spin systems (indeed this is rigorously proved for anti-ferromagnetic Ising model in [21]). With this to be true, devising approximation algorithms for two-state spin systems on degree-bounded graphs is reduced to verifying the strong spatial mixing on $d$-regular trees. This has been accepted as a fact about the two-state spin systems and has become a building block for approximation algorithms for such systems (see Theorem 2.8 in [21] and the discussions in [22, 23]). It was also raised as a conjecture in [22] whether the claim holds for general multi-state spin systems.

As a byproduct of our analysis (see Section 5), we find that this well-believed implication between the strong spatial mixing on $d$-regular tree and on graphs of maximum degree at most $d$ holds only for monotone spin systems (including the hardcore and anti-ferromagnetic Ising models). For general two-state spin systems the worst case for uniqueness as well as strong spatial mixing among all degree-bounded graphs is indeed a regular tree, but may no longer be the one of the highest degree. A bright side of this complication is that higher degrees may yield much faster correlation decay, making possible the FPTAS for graphs with unbounded degrees.

These new phenomena suggest that the general two-state spin systems have much richer structure than the well-studied monotone spin systems such as the hardcore and anti-ferromagnetic Ising models. The former approach via the strong spatial mixing on $d$-regular tree which succeeds in monotone spin systems on general graphs and general spin systems on regular graphs, meets a barrier when dealing with general spin systems on general graphs. We give a unified approach to the correlation decay in general two-state spin systems, through the strong spatial mixing on arbitrary trees instead of $d$-regular trees. This approach was initiated in our previous work [15] dealing with graphs of unbounded degrees. In this paper, we devise a unified potential-based analysis which adapts to both the irregularity of the arbitrary tree and the non-monotonicity of general two-state spin system and give tight correlation decay results for all two-state anti-ferromagnetic spin systems on degree-bounded families of graphs and family of all graphs.

1.1.2 Implications of the main results

By solving the uniqueness condition we can restate our main results in various threshold forms. Theorem 2 covers as special cases all previous algorithmic results for two-state anti-ferromagnetic spin systems on general-structure graphs as well as clears up previously uncovered cases.

In terms of interactions.

We can fix the external field $\lambda$ and discuss the tractable region of $(\beta ,\gamma )$. Since the roles of $\beta$ and $\gamma$ are symmetric, we can further fix one of them and discuss the tractable range of the other. This formulation was used in [10, 15].

Our main result can be restated as follows: for any $\mathrm{\Delta }$, there is a critical threshold ${\gamma }_c(\beta ,\lambda ,\mathrm{\Delta })$ for the uniqueness on infinite regular trees up to degree $\mathrm{\Delta }$ such that if there is an FPTAS for graphs of maximum degree at most $\mathrm{\Delta }$; and in particular, ${\gamma }_c={\gamma }_c(\beta ,\lambda ,\mathrm{\infty })>1$ is an absolute constant such that if $\gamma \in ({\gamma }_c,\frac{1}{\beta })$ there is an FPTAS for arbitrary graphs.

This covers as special cases all algorithmic results in [10] regarding the anti-ferromagnetic systems and all results in [15], extends the tractable regions in these previous works, and considers the degree-bounded graphs as well.

In terms of external field.

Motivated by the studies of hardcore and anti-ferromagnetic Ising models, we can fix $(\beta ,\gamma )$ and discuss the tractable range of external field.

Due to the symmetric role of $\beta$ and $\gamma$, we may assume that $\beta \le \gamma$. Our main results can be restated in specific settings as follows:

• Hard constraints (when $\beta =0$): For any $\mathrm{\Delta }$, is a critical threshold for the uniqueness on infinite regular trees up to degree $\mathrm{\Delta }$ such that if there exists an FPTAS for graphs of maximum degree at most $\mathrm{\Delta }$.

  For $\gamma \le 1$, the critical threshold equals ${\lambda }_c(\gamma ,\mathrm{\Delta })=\frac{{\gamma }^{\mathrm{\Delta }}{(\mathrm{\Delta }-1)}^{\mathrm{\Delta }-1}}{{(\mathrm{\Delta }-2)}^{\mathrm{\Delta }}}$. There is no external field $\lambda >0$ satisfying uniqueness on infinite regular trees of unbounded degrees. This is consistent with the hardness result for the hardcore model without degree bound [4, 23]. One particularly interesting case is when $\gamma =1$, in which case the model is the hardcore model with fugacity $\lambda$, and ${\lambda }_c(1,\mathrm{\Delta })=\frac{{(\mathrm{\Delta }-1)}^{\mathrm{\Delta }-1}}{{(\mathrm{\Delta }-2)}^{\mathrm{\Delta }}}$ is the critical threshold. This covers the result of [26].

  For $\gamma >1$, in addition to the results for degree-bounded graphs, there exists an absolute positive constant ${\lambda }_c(\gamma )={\min }_{d>1}\frac{{\gamma }^{d+1}{d}^d}{{(d-1)}^{d+1}}$ which lower bounds ${\lambda }_c(\gamma ,\mathrm{\Delta })$ for all $\mathrm{\Delta }$, such that if there exists an FPTAS for arbitrary graphs.

• Soft constraints (when $\beta >0$): If $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$ the uniqueness condition holds for any external field $\lambda >0$ and there always exists an FPTAS for graphs of maximum degree at most $\mathrm{\Delta }$.

  Now suppose $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$. Let $\overline{\mathrm{\Delta }}$ be the smallest $d$ satisfying $\sqrt{\beta \gamma }\le \frac{d-1}{d+1}$. For every integer , there exist two thresholds ${\lambda }_1(\beta ,\gamma ,d)$ and ${\lambda }_2(\beta ,\gamma ,d)$ for the uniqueness on the infinite $(d+1)$-regular tree, such that if or $\lambda >{\lambda }_2(\beta ,\gamma ,d)$ for all integers , there is an FPTAS for graphs of maximum degree at most $\mathrm{\Delta }$.

  Furthermore, if $\gamma \le 1$, the above uniqueness condition can be simplified as that or $\lambda >{\overline{\lambda }}_c$, where and . In particular, for $\beta =\gamma$, the anti-ferromagnetic Ising model, we have ${\overline{\lambda }}_c=1/{\lambda }_c$, and the uniqueness holds when $|\log \lambda |>\log {\overline{\lambda }}_c$. This covers the result of anti-ferro-Ising model in [21].

  For unbounded maximum degree $\mathrm{\Delta }$, if $\gamma \le 1$, there is no $\lambda >0$ satisfying the uniqueness on infinite regular trees of all degrees, which is consistent with the hardness results in [10]; and if $\gamma >1$, the uniqueness for all infinite regular trees holds when or $\lambda >{\lambda }_2(\beta ,\gamma ,d)$ for all $d\ge \overline{\mathrm{\Delta }}$, under which condition there exists an FPTAS for arbitrary graphs.

1.2 Related works

The approximation of partition functions of spin systems has been extensively studied[12, 13, 9, 3, 11, 25, 4, 5, 17]. In a seminal work [12], Jerrum and Sinclair devised a fully polynomial-time randomized approximation scheme (FPRAS) for the ferromagnetic Ising model. Later in [10], the FPRAS was extended to all two-state ferromagnetic spin systems by translating the parameters to the ferromagnetic Ising model. Also in [10], Goldberg, Jerrum, and Paterson gave an FPRAS and inapproximability results for two-state anti-ferromagnetic spin systems on arbitrary graphs. A gadget based on random regular bipartite graphs was proposed by Dyer, Frieze, and Jerrum in [4] and was also analyzed by Mossel, Weitz, and Wormald in [19] to study the inapproximability on degree-bounded graphs. It is widely believed that the transition of approximability of anti-ferromagnetic spin systems is captured by the phase transition of uniqueness on infinite trees. This was raised openly as a conjecture in [19]. The conjecture was proved by Sly in [23] for the hardcore model. This result was improved by Galanis et al. in [6] to a wide range of parameters. Very recently, Sly and Sun [24] proved the hardness of all two-state anti-ferromagnetic spin systems of non-uniqueness on infinite regular trees. A same result for anti-ferromagnetic Ising model without external field was independently proved by Galanis, S̆tefankovic̆, and Vigoda in [7].

The correlation decay technique developed independently by Weitz [26] and Bandyopadhyay and Gamarnik [1] is a powerful tool for devising deterministic fully polynomial-time approximation schemes (FPTAS) for partition functions (other important examples include [8, 2]). In [26], Weitz introduced the concept of strong spatial mixing and used it to devise FPTAS for the hardcore model up to the uniqueness threshold. The other most important two-state anti-ferromagnetic spin system, the anti-ferromagnetic Ising model, was studied recently by Sinclair, Srivastava, and Thurley in [21], where a more powerful message-decay method was introduced to analyze the strong spatial mixing and give FPTAS up to uniqueness threshold. A powerful technique was developed by Restrepo et al. in [20] which makes use of the specific structure of graphs for strong spatial mixing. A broader tractable region than the region of uniqueness is achieved on grid lattice by exploiting the structure of the graph. In a previous work [15], we gave an FPTAS for two-state anti-ferromagnetic spin systems without external field on arbitrary graphs with unbounded degrees, up to a continuous relaxation of the uniqueness threshold. The approach used in the current paper was initiated in [15], however the analysis in [15] cannot separate the uniqueness up to certain degree, thus fails in dealing with degree-bounded families of graphs.

2 Definitions and preliminaries

A two-state spin system is described by a graph $G=(V,E)$. A configuration of the system is one of the $2^{|V|}$ possible assignments $\sigma :V\to \{0,1\}$ of states to vertices. We also use two colors blue and green to denote these two states. The weight of a configuration can be described as a product of contributions of individual edges and vertices. Let and $b=(b_0,b_1)$. The weight of a configuration $\sigma :V\to \{0,1\}$ is given by

$w(\sigma )={\displaystyle \prod _{(u,v)\in E}}{A}_{\sigma (u),\sigma (v)}{\displaystyle \prod _{v\in V}}{b}_{\sigma (v)}.$

The Gibbs measure is a probability distribution over all configurations defined by $\rho (\sigma )=\frac{w(\sigma )}{Z(G)}$. The normalization factor $Z(G)={\sum }_{\sigma }w(\sigma )$ is called the partition function.

We can normalize the contributions of a $\{\mathrm{𝚋𝚕𝚞𝚎},\mathrm{𝚐𝚛𝚎𝚎𝚗}\}$ edge and of a green vertex to be 1. So that for some $\beta ,\gamma \ge 0$, and $b=(b_0,b_1)=(\lambda ,1)$ for some $\lambda >0$. Since the roles of blue and green are symmetric, we can assume that $\beta \le \gamma$ without loss of generality. The three parameters $(\beta ,\gamma ,\lambda )$ with $0\le \beta \le \gamma$ and $\lambda >0$ completely specify a two-state spin system. A two-state spin system with $\beta =\gamma$ is an Ising model and a two-state spin system with $\beta =0,\gamma =1$ or symmetrically $\beta =1,\gamma =0$ is a hardcore model.

A two-state spin system is called anti-ferromagnetic if adjacent vertices favor disagreeing spins, i.e. if . Without loss of generality, we focus on the cases that $\beta \le \gamma$.

Definition 4.

$(\beta ,\gamma ,\lambda )$ is anti-ferromagnetic if $0\le \beta \le \gamma$, $\gamma >0$, , and $\lambda >0$.

By the symmetry of $\beta$ and $\gamma$ and the triviality of the case $\beta =\gamma =0$, this definition is complete for all nontrivial two-state anti-ferromagnetic systems

2.1 Correlation decay

The Gibbs measure defines a marginal distribution of state for each vertex. Let $p_v$ denote the probability of vertex $v$ to be colored blue. Let ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ be a configuration of $\mathrm{\Lambda }\subset V$. We call vertices $v\in \mathrm{\Lambda }$ fixed vertices, and $v\notin \mathrm{\Lambda }$ free vertices. We use $p_v^{{\sigma }_{\mathrm{\Lambda }}}$ to denote the marginal probability of $v$ to be colored blue conditioning on the configuration of $\mathrm{\Lambda }$ being fixed as ${\sigma }_{\mathrm{\Lambda }}$.

Definition 5.

A spin system on a family of graphs is said to be of strong spatial mixing if for any graph $G=(V,E)$ in the family, any $v\in V,\mathrm{\Lambda }\subset V$ and ${\sigma }_{\mathrm{\Lambda }},{\tau }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$,

$|p_v^{{\sigma }_{\mathrm{\Lambda }}}-p_v^{{\tau }_{\mathrm{\Lambda }}}|\le \exp (-\mathrm{\Omega }(\mathrm{dist}(v,S))),$

where $S\subset \mathrm{\Lambda }$ is the subset on which ${\sigma }_{\mathrm{\Lambda }}$ and ${\tau }_{\mathrm{\Lambda }}$ differ, and $\mathrm{dist}(v,S)$ is the shortest distance from $v$ to any vertex in $S$.

The weak spatial mixing can be defined by measuring the decay with respect to $\mathrm{dist}(v,\mathrm{\Lambda })$ instead of $\mathrm{dist}(v,S)$. The spatial mixing property is also called correlation decay in Statistical Physics.

Let $T$ be a tree rooted by $v$. Define $R_T^{{\sigma }_{\mathrm{\Lambda }}}=p_v^{{\sigma }_{\mathrm{\Lambda }}}/(1-p_v^{{\sigma }_{\mathrm{\Lambda }}})$ to be the ratio between the probabilities that the root $v$ is blue and green, when imposing the condition ${\sigma }_{\mathrm{\Lambda }}$ (with the convention that $R_T^{{\sigma }_{\mathrm{\Lambda }}}=\mathrm{\infty }$ when $p_v^{{\sigma }_{\mathrm{\Lambda }}}=1$). Suppose that $v$ has $d$ children and $T_i$ is the subtree rooted by the $i$-th child. Due to the independence of subtrees, we have an easy recursion for calculating $R_T^{{\sigma }_{\mathrm{\Lambda }}}$:

$R_T^{{\sigma }_{\mathrm{\Lambda }}}$ $=\lambda {\displaystyle \prod _{i=1}^d}{\displaystyle \frac{\beta R_{T_i}^{{\sigma }_{\mathrm{\Lambda }}}+1}{R_{T_i}^{{\sigma }_{\mathrm{\Lambda }}}+\gamma }}.$ (1)

Let $G(V,E)$ be a graph. Similarly define that $R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}=p_v^{{\sigma }_{\mathrm{\Lambda }}}/(1-p_v^{{\sigma }_{\mathrm{\Lambda }}})$. In contrast to the case of tree, there is no easy recursion for calculating $R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}$ for a general graph $G$ because of the dependencies caused by cycles. In [26], a construction called the self-avoiding walk (SAW) tree was introduced which reduces the computing of marginal distribution in a general graph to that in a tree. Specifically, given a graph $G(V,E)$ and a vertex $v\in V$. The SAW tree $T_{\text{SAW}}(G,v)$ is a tree rooted at $v$ with a new vertex set $V_{\text{SAW}}$ (which effectively enumerates all paths originating from $v$ in $G$ and may include fixed leaves). Moreover, any vertex sets $S\subset \mathrm{\Lambda }\subset V$ are mapped to respective $S_{\text{SAW}}\subset {\mathrm{\Lambda }}_{\text{SAW}}\subset V_{\text{SAW}}$ and any configuration ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ is mapped to a ${\sigma }_{{\mathrm{\Lambda }}_{\text{SAW}}}\in {\{0,1\}}^{{\mathrm{\Lambda }}_{\text{SAW}}}$. We abuse the notation and write $S=S_{\text{SAW}}$ and ${\sigma }_{\mathrm{\Lambda }}={\sigma }_{{\mathrm{\Lambda }}_{\text{SAW}}}$ if no ambiguity is caused. Given a graph $G(V,E)$, $v\in V$ and $S\subset V$, let ${\mathrm{dist}}_G(v,S)$ be the shortest distance in $G$ from $v$ to any vertex in $S$.

Theorem 6 (Theorem 3.1 of Weitz [26]).

Let $G(V,E)$ be a graph, $v\in V$, ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ a configuration on $\mathrm{\Lambda }\subset V$, and $S\subset V$. Let $T=T_{\mathrm{SAW}}(G,v)$. It holds that the maximum degree of $T$ equals the maximum degree of $G$, ${\mathrm{dist}}_G(v,S)={\mathrm{dist}}_T(v,S)$, and $R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}=R_T^{{\sigma }_{\mathrm{\Lambda }}}$. Moreover, any neighborhood of $v$ in $T$ can be constructed in time proportional to the size of the neighborhood.

2.2 The uniqueness condition

We consider the uniqueness of Gibbs measure on the infinite $(d+1)$-regular trees ${\widehat{𝕋}}^d$, in which the recursion is given by $f_d(x)\triangleq \lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d$ due to the symmetric structure of ${\widehat{𝕋}}^d$.

Let ${\widehat{x}}_d$ be the positive fixed point of $f_d(x)$, that is, ${\widehat{x}}_d=f({\widehat{x}}_d)$. It is known [14, 18] that the two-state anti-ferromagnetic spin system on ${\widehat{𝕋}}^d$ undergoes a phase transition at $|f_d^{\prime }({\widehat{x}}_d)|=1$ with uniqueness when $|f_d^{\prime }({\widehat{x}}_d)|\le 1$. This motivates the following definition.

Definition 7.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. Let ${\widehat{x}}_d$ be the positive fixed point of function $f_d(x)=\lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d$. We say that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, if for all integers ,

In particular, $(\beta ,\gamma ,\lambda )$ is universally unique if it is up-to-$\mathrm{\infty }$ unique.

Being up-to-$\mathrm{\Delta }$ unique is equivalent to that the system is of weak spatial mixing on infinite regular trees up to degree $\mathrm{\Delta }$. The uniqueness condition can be described in various threshold forms, which are given in Appendix A.

The uniqueness is defined by requiring that . The following lemma states that $|f_d^{\prime }({\widehat{x}}_d)|$ is bounded by an absolute constant as long as the uniqueness condition holds.

Lemma 8.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. If $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique then there exists an absolute constant which depends only on $\beta$, $\gamma$, $\lambda$ and $\mathrm{\Delta }$, such that $|f_d^{\prime }({\widehat{x}}_d)|\le c$ for all .

Proof.

The lemma holds trivially for finite $\mathrm{\Delta }$. It then remains to show that in case of universal uniqueness, $|f_d^{\prime }({\widehat{x}}_d)|$ cannot be arbitrarily close to 1 as $d$ grows to infinity. If $(\beta ,\gamma ,\lambda )$ is universally unique, due to Lemma 21.2, we must have $\gamma >1$. For anti-ferromagnetic $(\beta ,\gamma ,\lambda )$, $\beta \le \frac{1}{\gamma }$, thus the fixed point ${\widehat{x}}_d=\lambda {\left(\frac{\beta {\widehat{x}}_d+1}{{\widehat{x}}_d+\gamma }\right)}^d\le \frac{\lambda }{{\gamma }^d}$, therefore $|f_d^{\prime }({\widehat{x}}_d)|=\frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}\le \frac{d\lambda }{{\gamma }^d}$. The lemma follows. ∎

3 The strong spatial mixing on general graphs

In this section we prove Theorem 1. The necessity of the uniqueness condition is trivial since strong spatial mixing on general graphs implies weak spatial mixing on regular trees. It then remains to prove the following theorem.

Theorem 9.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. For any finite $\mathrm{\Delta }\ge 2$ or $\mathrm{\Delta }=\mathrm{\infty }$, if $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, then the two-state spin system of parameters $(\beta ,\gamma ,\lambda )$ is of strong spatial mixing on graphs of maximum degree at most $\mathrm{\Delta }$.

Our approach is to prove the strong spatial mixing on arbitrary trees of maximum degree at most $\mathrm{\Delta }$, which by the self-avoiding walk tree construction implies the theorem. Note that unlike in [26] and [21], we analyze the decay on arbitrary tree instead of regular tree. This is because for general two-state anti-ferromagnetic spin systems the worst case for strong spatial mixing among all graphs of maximum degree at most $d$ may no longer be the $d$-regular tree. We will explain this in detail in Section 5.

Given any graph $G(V,E)$ of maximum degree at most $\mathrm{\Delta }$, any configuration ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ on $\mathrm{\Lambda }\subset V$ and any $S\subset \mathrm{\Lambda }$, fixing an arbitrary vertex $v\in V$, by Theorem 6, a self-avoiding walk tree $T=T_{\text{SAW}}(G,v)$ can be constructed such that the maximum degree of $T$ is bounded by $\mathrm{\Delta }$, ${\mathrm{dist}}_G(v,S)={\mathrm{dist}}_T(v,S)$ and $R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}=R_T^{{\sigma }_{\mathrm{\Lambda }}}$. Recall that $R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}=p_v^{{\sigma }_{\mathrm{\Lambda }}}/(1-p_v^{{\sigma }_{\mathrm{\Lambda }}})$ thus $p_v^{{\sigma }_{\mathrm{\Lambda }}}=R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}/(1+R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}})$. For any ${\sigma }_{\mathrm{\Lambda }}$ and ${\tau }_{\mathrm{\Lambda }}$,

$\left|p_v^{{\sigma }_{\mathrm{\Lambda }}}-p_v^{{\tau }_{\mathrm{\Lambda }}}\right|=\left|{\displaystyle \frac{R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}}{1+R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}}}-{\displaystyle \frac{R_{G,v}^{{\tau }_{\mathrm{\Lambda }}}}{1+R_{G,v}^{{\tau }_{\mathrm{\Lambda }}}}}\right|\le \left|R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}-R_{G,v}^{{\tau }_{\mathrm{\Lambda }}}\right|=\left|R_T^{{\sigma }_{\mathrm{\Lambda }}}-R_T^{{\tau }_{\mathrm{\Lambda }}}\right|.$

This motivates the following definition.

Definition 10.

Let $T$ be a tree rooted by vertex $v$, ${\tau }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ be a configuration on $\mathrm{\Lambda }\subset V$ and $S\subseteq \mathrm{\Lambda }$ be a vertex set. Define $R_v$ and ${\delta }_v$ as that $R_v\le R_T^{{\tau }_{\mathrm{\Lambda }}}\le R_v+{\delta }_v$ for all ${\sigma }_{\mathrm{\Lambda }}\in {\{0,1\}}^{\mathrm{\Lambda }}$ which differ from ${\tau }_{\mathrm{\Lambda }}$ only on $S$.

It is then sufficient to prove Theorem 9 by constructing such $R_v$ and ${\delta }_v$ for $T=T_{\text{SAW}}(G,v)$ and showing that ${\delta }_v\le \exp (-\mathrm{\Omega }(\mathrm{dist}(v,S)))$.

Let $T$ be a tree rooted by $v$, who has $d$ children $v_1,\mathrm{\dots },v_d$, and $T_i$ be the subtree rooted by $v_i$. It holds that

$R_T^{{\sigma }_{\mathrm{\Lambda }}}$ $=f(R_{T_1}^{{\sigma }_{\mathrm{\Lambda }}},\mathrm{\dots },R_{T_d}^{{\sigma }_{\mathrm{\Lambda }}})\triangleq \lambda {\displaystyle \prod _{i=1}^d}{\displaystyle \frac{\beta R_{T_i}^{{\sigma }_{\mathrm{\Lambda }}}+1}{R_{T_i}^{{\sigma }_{\mathrm{\Lambda }}}+\gamma }}.$

The lower and upper bounds $R_v$ and $R_v+{\delta }_v$ can be recursively constructed as follows. The base cases are: (1) $v\in S$, in which case $R_v=0$ and ${\delta }_v=\mathrm{\infty }$; and (2) $v\in \mathrm{\Lambda }\setminus S$, i.e. $v$ is fixed to be the same color in all ${\sigma }_{\mathrm{\Lambda }}$, in which case ${\delta }_v=0$ and $R_v=\mathrm{\infty }$ (or $R_v=0$) if $v$ is fixed to be blue (or green). For $v\notin \mathrm{\Lambda }$, since $f(R_1,\mathrm{\dots },R_d)$ is monotonically decreasing for anti-ferromagnetic $(\beta ,\gamma ,\lambda )$,

$R_v$ $=f(R_{v_1}+{\delta }_{v_1},\mathrm{\dots },R_{v_d}+{\delta }_{v_d}),$

$R_v+{\delta }_v$ $=f(R_{v_1},\mathrm{\dots },R_{v_d}),$

where $R_{v_i}$ and $R_{v_i}+{\delta }_{v_i}$ are the corresponding lower and upper bounds for $R_{T_i}^{{\sigma }_{\mathrm{\Lambda }}}$, $1\le i\le d$. In particular, when $d=0$ the empty product equals 1 by convention, thus $R_v=R_v+{\delta }_v=\lambda$, which is consistent with the case that $v$ is a free vertex having no children.

By the monotonicity of $f(R_1,\mathrm{\dots },R_d)$, it is easy to check that the $R_v$ and $R_v+{\delta }_v$ constructed above satisfy the requirement of Definition 10. Our goal is to show that ${\delta }_v$ decays exponentially in depth of recursion when the uniqueness holds. A straightforward approach is trying to prove that $\delta$ contracts at a constant rate in each recursion step. But this does not have to be true to guarantee the exponential decay. Indeed there are cases that the error does not decay in single steps but decay in a long run. To overcome this, we use a potential $\mathrm{\Phi }$ to amortize the contraction and show that $\delta \cdot \mathrm{\Phi }$ contracts at a constant rate. We choose the potential function to be

$\mathrm{\Phi }(R)={\displaystyle \frac{1}{\sqrt{R(\beta R+1)(R+\gamma )}}}.$

We are analyzing the decay on an arbitrary tree with irregular degrees. In order to adapt this irregularity, the potential function cannot have $d$ as an input, but only caries the information regarding the distribution at the current vertex, yet it has to be able to provide correct compensation to the step-wise decay at any state of $R$ and for all spin systems satisfying sufficient uniqueness. A heuristic procedure which leads us to this good potential function is discussed in Appendix C.

Let $\phi (R)$ be a monotone function satisfying that ${\phi }^{\prime }(R)=\mathrm{\Phi }(R)$. We define that

${ϵ}_v\triangleq \phi (R_v+{\delta }_v)-\phi (R_v),$

and accordingly, ${ϵ}_{v_i}\triangleq \phi (R_{v_i}+{\delta }_{v_i})-\phi (R_{v_i})$, $1\le i\le d$.

We define a function $\alpha (d;x_1,\mathrm{\dots },x_d)$ as follows:

$\alpha (d;x_1,\mathrm{\dots },x_d)$ $\triangleq {\displaystyle \frac{(1-\beta \gamma ){\left(\lambda {\prod }_{i=1}^d\frac{\beta x_i+1}{x_i+\gamma }\right)}^{\frac{1}{2}}}{{\left(\beta \lambda {\prod }_{i=1}^d\frac{\beta x_i+1}{x_i+\gamma }+1\right)}^{\frac{1}{2}}{\left(\lambda {\prod }_{i=1}^d\frac{\beta x_i+1}{x_i+\gamma }+\gamma \right)}^{\frac{1}{2}}}}\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{x_i^{\frac{1}{2}}}{{(\beta x_i+1)}^{\frac{1}{2}}{(x_i+\gamma )}^{\frac{1}{2}}}}.$

The following lemma is obtained from applying the Mean Value Theorem. Similar routines were previously used in [15, 20]. The proof of the lemma is in Appendix B.

Lemma 11.

The followings hold for ${ϵ}_v$.

1. (relation to ${\delta }_v$) ${ϵ}_v={\delta }_v\cdot \mathrm{\Phi }(\tilde{R})$ for some $\tilde{R}\in [R_v,R_v+{\delta }_v]$.

2. (absolute bound) Assuming that $\gamma >1$ or the maximum degree of $T$ is bounded by a constant, if $v\notin \mathrm{\Lambda }$ then $R_v+{\delta }_v=O(1)$ and if $v_i\notin S$ for all $1\le i\le d$ then ${ϵ}_v=O(1)$.

3. (stepwise contraction) There exist $\widetilde{R_i}\in [R_{v_i},R_{v_i}+{\delta }_{v_i}]$, $1\le i\le d$, such that

   ${ϵ}_v\le \alpha (d;{\tilde{R}}_1,\mathrm{\dots },{\tilde{R}}_d)\cdot \underset{1\le i\le d}{\max }\{{ϵ}_{v_i}\}.$

To prove the strong spatial mixing, we first relate ${\delta }_v$ to ${ϵ}_v$ by Item 1 of Lemma 11, and then apply induction on the depth in $T$, with Item 2 of Lemma 11 as basis and Item 3 of Lemma 11 as induction step. We then need to bound the contraction rate $\alpha (d;x_1,\mathrm{\dots },x_d)$. Note that $\frac{z}{(\beta z+1)(z+\gamma )}\le {\left(1+\sqrt{\beta \gamma }\right)}^{-2}\le 1$ for $z\in [0,\mathrm{\infty })$, thus it holds unconditionally for all $x_i\in [0,\mathrm{\infty })$, $1\le i\le d$, that

$\alpha (d;x_1,\mathrm{\dots },x_d)$ $\le d,$ (2)

$\alpha (d;x_1,\mathrm{\dots },x_d)$ $\le d\cdot \sqrt{{\displaystyle \frac{\lambda }{{\gamma }^{d+1}}}}.$ (3)

With the uniqueness, the following much tighter contraction bound can be proved.

Lemma 12.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. If $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, then there exists a constant such that for any integer and any $x_1,\mathrm{\dots },x_d\in [0,+\mathrm{\infty })$, $1\le i\le d$, it holds that $\alpha (d;x_1,\mathrm{\dots },x_d)\le \alpha$.

This lemma is the technical core of our analysis. It crucially relies on the choice of potential function. Before delving into the formal proof of Lemma 12, we note that Theorem 9 can be implied by this lemma.

Proof of Theorem 9.

Let $T=T_{\text{SAW}}(G,v)$ for a $G$ whose maximum degree is at most $\mathrm{\Delta }$. Then the maximum degree of $T$ is at most $\mathrm{\Delta }$, thus the root $v$ has at most $\mathrm{\Delta }$ children in $T$, and every other vertex in $T$ has less than $\mathrm{\Delta }$ children. We recursively construct $R_u$, ${\delta }_u$ and ${ϵ}_u$ for every subtree in $T$.

Let $t=\mathrm{dist}(v,S)$. By repeatedly applying Item 3 of Lemma 11, without loss of generality, we have a path $u_1{u}_2\mathrm{\cdots }{u}_{t-2}$ in $T$ with $u_1=v$ such that ${ϵ}_{u_j}\le \alpha (d_j;x_1,\mathrm{\dots },x_{d_j})\cdot {ϵ}_{u_{j+1}}$ for $j=1,2,\mathrm{\dots },t-3$, where $d_j$ is the number of children of $u_j$ and $x_i\in [0,\mathrm{\infty })$, $1\le i\le d_j$.

Note that $d_1\le \mathrm{\Delta }$, and for all other $j$. Assume that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique. If $\mathrm{\Delta }$ is bounded, then by Lemma 12 there exists a constant , such that ${ϵ}_{u_j}\le \alpha \cdot {ϵ}_{u_{j+1}}$ for $2\le j\le t-3$, and ${ϵ}_v\le d_1\cdot {ϵ}_{u_2}\le \mathrm{\Delta }\cdot {ϵ}_2$ due to (2); and if $\mathrm{\Delta }=\mathrm{\infty }$, then by Lemma 12, ${ϵ}_{u_j}\le \alpha \cdot {ϵ}_{u_{j+1}}$ for all $1\le j\le t-3$. In both cases we have ${ϵ}_v=O({\alpha }^t\cdot {ϵ}_{u_{t-2}})$.

Due to Item 1 of Lemma 11, ${\delta }_v=\frac{{ϵ}_v}{\mathrm{\Phi }(\tilde{R})}=O\left(\frac{1}{\mathrm{\Phi }(\tilde{R})}\cdot {\alpha }^t{ϵ}_{u_{t-2}}\right)$ for some $\tilde{R}\in [R_v,R_v+{\delta }_v]$. We then bound $\mathrm{\Phi }(\tilde{R})$ and ${ϵ}_{u_{t-2}}$. Due to Item 2 of Lemma 21, the fact that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique implies that either $\mathrm{\Delta }$ is bounded or $\gamma >1$. Note that $v$ must be free or the theorem is trivial to prove, and none of $u_{t-2}$’s children is in $S$ because $\mathrm{dist}(v,S)=t$. Thus by Item 2 of Lemma 11, ${ϵ}_{u_{t-2}}=O(1)$ and $\tilde{R}\le R_v+{\delta }_v=O(1)$, which implies that $\mathrm{\Phi }(\tilde{R})=\frac{1}{\sqrt{\tilde{R}(\beta \tilde{R}+1)(\tilde{R}+\gamma )}}=\mathrm{\Omega }(1)$.

In conclusion, if $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, there exists a constant , such that ${\delta }_v=O\left({\alpha }^t\right)$. As discussed in the beginning of this section, this proves Theorem 9.

Proof of Lemma 12.

The rest of this section is devoted to the proof of Lemma 12. Given that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, there exists an absolute constant such that $\alpha (d;x_1,\mathrm{\dots },x_d)\le \alpha$ for any and $x_1,\mathrm{\dots },x_d\ge 0$.

We define the symmetric version of $\alpha (d;x_1,\mathrm{\dots },x_d)$:

${\alpha }_d(x)$ $\triangleq \alpha (d;\underset{d}{\underbrace{x,\mathrm{\dots },x}})={\displaystyle \frac{d(1-\beta \gamma ){\left(x\cdot \lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d\right)}^{\frac{1}{2}}}{{(\beta x+1)}^{\frac{1}{2}}{(x+\gamma )}^{\frac{1}{2}}{\left(\beta \lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d+1\right)}^{\frac{1}{2}}{\left(\lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d+\gamma \right)}^{\frac{1}{2}}}}.$

The following lemma shows that the symmetric case dominates the maximum of $\alpha (d;x_1,\mathrm{\dots },x_d)$ by using the inequalities of Cauchy-Schwarz and arithmetic and geometric means.

Lemma 13.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. Then for any integer $d$ and any $x_1,\mathrm{\dots },x_d\in [0,+\mathrm{\infty })$, there exists an $\overline{x}\in [0,+\mathrm{\infty })$ such that $\alpha (d;x_1,\mathrm{\dots },x_d)\le \alpha (d,\overline{x})$.

Proof.

Let $z_i=\frac{\beta x_i+1}{x_i+\gamma }$. Then $z_i\in (\beta ,\frac{1}{\gamma }]$ and $x_i=\frac{1-\gamma z_i}{z_i-\beta }$. Express $\alpha (d;x_1,\mathrm{\dots },x_d)$ in terms of $z_i$:

$\alpha (d;x_1,\mathrm{\dots },x_d)$ $={\displaystyle \frac{{\left(\lambda {\prod }_{i=1}^d{z}_i\right)}^{\frac{1}{2}}}{{\left(\beta \lambda {\prod }_{i=1}^d{z}_i+1\right)}^{\frac{1}{2}}{\left(\lambda {\prod }_{i=1}^d{z}_i+\gamma \right)}^{\frac{1}{2}}}}\cdot {\displaystyle \sum _{i=1}^d}{(z_i^{-1}-\gamma )}^{\frac{1}{2}}{(z_i-\beta )}^{\frac{1}{2}}.$

Due to Cauchy-Schwarz inequality,

${\displaystyle \sum _{i=1}^d}{(z_i^{-1}-\gamma )}^{\frac{1}{2}}{(z_i-\beta )}^{\frac{1}{2}}\le d{\left({\displaystyle \frac{1}{d}}{\displaystyle \sum _{i=1}^d}(z_i^{-1}-\gamma )(z_i-\beta )\right)}^{\frac{1}{2}}=d{\left(1+\beta \gamma -{\displaystyle \frac{1}{d}}{\displaystyle \sum _{i=1}^d}(z_i\gamma +\beta z_i^{-1})\right)}^{\frac{1}{2}}.$

Due to the inequality of arithmetic and geometric means,

$d{\left(1+\beta \gamma -{\displaystyle \frac{1}{d}}{\displaystyle \sum _{i=1}^d}(z_i\gamma +\beta z_i^{-1})\right)}^{\frac{1}{2}}$ $\le d{\left(1+\beta \gamma -\gamma {\left({\displaystyle \prod _{i=1}^d}{z}_i\right)}^{\frac{1}{d}}-\beta {\left({\displaystyle \prod _{i=1}^d}{z}_i\right)}^{-\frac{1}{d}}\right)}^{\frac{1}{2}}.$

Let $\overline{z}={\left({\prod }_{i=1}^d{z}_i\right)}^{\frac{1}{d}}$. Then combining the above calculations,

$\alpha (d;x_1,\mathrm{\dots },x_d)$ $\le {\displaystyle \frac{{(\lambda {\overline{z}}^d)}^{\frac{1}{2}}\cdot d{(1+\beta \gamma -\gamma \overline{z}-\beta {\overline{z}}^{-1})}^{\frac{1}{2}}}{{(\beta \lambda {\overline{z}}^d+1)}^{\frac{1}{2}}{(\lambda {\overline{z}}^d+\gamma )}^{\frac{1}{2}}}}=d\cdot \sqrt{{\displaystyle \frac{\lambda {\overline{z}}^d({\overline{z}}^{-1}-\gamma )(\overline{z}-\beta )}{(\beta \lambda {\overline{z}}^d+1)(\lambda {\overline{z}}^d+\gamma )}}}.$

Let $\overline{x}$ be such that $\frac{\beta \overline{x}+1}{\overline{x}+\gamma }=\overline{z}$. Then $\overline{x}\in [0,+\mathrm{\infty })$ and by substituting $\frac{\beta \overline{x}+1}{\overline{x}+\gamma }$ for $\overline{z}$, we have

$\alpha (d;x_1,\mathrm{\dots },x_d)$ $\le {\displaystyle \frac{d(1-\beta \gamma ){\left(\overline{x}\cdot \lambda {\left(\frac{\beta \overline{x}+1}{\overline{x}+\gamma }\right)}^d\right)}^{\frac{1}{2}}}{{(\beta \overline{x}+1)}^{\frac{1}{2}}{(\overline{x}+\gamma )}^{\frac{1}{2}}{\left(\beta \lambda {\left(\frac{\beta \overline{x}+1}{\overline{x}+\gamma }\right)}^d+1\right)}^{\frac{1}{2}}{\left(\lambda {\left(\frac{\beta \overline{x}+1}{\overline{x}+\gamma }\right)}^d+\gamma \right)}^{\frac{1}{2}}}}={\alpha }_d(\overline{x}).$

∎

Lemma 14.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. If $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, then there exists a constant such that for any integer , it holds that ${\alpha }_d(x)\le \alpha$ for all $x\ge 0$.

Proof.

Fix $d$ to be any positive integer. We characterize the value of $x$ at which ${\alpha }_d(x)$ achieves its maximum. We denote that $f_d(x)=\lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d$. Taking derivative of ${\alpha }_d(x)$ with respect to $x$, we get that

${\alpha }_d^{\prime }(x)=d(1-\beta \gamma )\cdot {\displaystyle \frac{G^{\prime }(x)}{2\sqrt{G(x)}}},$

where $G(x)=\frac{x{f}_d(x)}{(\beta f_d(x)+1)(f_d(x)+\gamma )(\beta x+1)(x+\gamma )}$, whose derivative is

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

As $x$ ranges over $[0,\mathrm{\infty })$, the function $\frac{\gamma -\beta x^2}{d(1-\beta \gamma )x}$ is strictly decreasing in $x$ and ranges from $+\mathrm{\infty }$ to $-\mathrm{\infty }$, and the function $\frac{\gamma -\beta f_d{(x)}^2}{\left(\beta f_d(x)+1\right)\left(f_d(x)+\gamma \right)}$ is strictly increasing in $x$ and has a bounded range. Thus, the equation

${\displaystyle \frac{\gamma -\beta x^2}{d(1-\beta \gamma )x}}$ $={\displaystyle \frac{\gamma -\beta f_d{(x)}^2}{\left(\beta f_d(x)+1\right)\left(f_d(x)+\gamma \right)}}.$ (4)

has unique solution in $(0,\mathrm{\infty })$, denoted by $x_d$. Moreover, it holds that

(5)

The same also holds for ${\alpha }_d^{\prime }(x)$. Thus, for any fixed $d$, ${\alpha }_d(x)$ achieves its maximum when $x=x_d$.

Therefore, for all $x\ge 0$,

${\alpha }_d(x)\le {\alpha }_d(x_d)$ $=d(1-\beta \gamma ){\left({\displaystyle \frac{x_d{f}_d(x_d)}{(\beta x_d+1)(x_d+\gamma )\left(\beta f_d(x_d)+1\right)\left(f_d(x_d)+\gamma \right)}}\right)}^{\frac{1}{2}}$

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

$\triangleq {\tilde{\alpha }}_d(x_d).$

Equation (6) is obtained by substituting $\left(\beta f_d(x_d)+1\right)\left(f_d(x_d)+\gamma \right)$ according to (4).

Let ${\widehat{x}}_d$ be the positive fixed point of $f_d(x)$, that is, ${\widehat{x}}_d=f_d({\widehat{x}}_d)$. We then claim that if $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, then ${\tilde{\alpha }}_d(x_d)\le {\tilde{\alpha }}_d({\widehat{x}}_d)$ for any integral , To see that this claim is sufficient to imply the lemma, note that after substituting ${\widehat{x}}_d=f_d({\widehat{x}}_d)$, we have ${\tilde{\alpha }}_d({\widehat{x}}_d)=\sqrt{\frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}}=\sqrt{\left|f_d^{\prime }({\widehat{x}}_d)\right|}$. And due to Lemma 8, if $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique then there exists a constant such that for all integer .

We then prove the claim. Assume that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique and . It is then sufficient to show that ${\tilde{\alpha }}_d(x)$ is decreasing if ${\widehat{x}}_d\le x_d$ and is increasing if ${\widehat{x}}_d>x_d$.

Case 1: ${\widehat{x}}_d\le x_d$. Due to (5), $G^{\prime }({\widehat{x}}_d)\ge 0$. Note that

$G^{\prime }({\widehat{x}}_d)={\displaystyle \frac{d(1-\beta \gamma )(\gamma -\beta {\widehat{x}}_d^2){\widehat{x}}_d^2}{{(\beta {\widehat{x}}_d+1)}^3{({\widehat{x}}_d+\gamma )}^3}}\cdot \left({\displaystyle \frac{1}{d(1-\beta \gamma ){\widehat{x}}_d}}-{\displaystyle \frac{1}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}}\right).$

Due to the uniqueness, , thus $\frac{1}{d(1-\beta \gamma ){\widehat{x}}_d}-\frac{1}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}>0$. Combining with that $G^{\prime }(x)\ge 0$, we have $\gamma -\beta {\widehat{x}}_d^2\ge 0$. Since the function $f_d(x)$ is monotonically decreasing and ${\widehat{x}}_d$ is its fixed point, $\gamma -\beta f_d{(x_d)}^2\ge \gamma -\beta f_d{({\widehat{x}}_d)}^2=\gamma -\beta {\widehat{x}}_d^2\ge 0$. Since $x_d$ satisfies (4), $\gamma -\beta x_d^2$ and $\gamma -\beta f_d{(x_d)}^2$ must be simultaneously positive or negative, thus it also holds that $\gamma -\beta x_d^2\ge 0$. Then both $\frac{(\gamma -\beta x^2)}{(\beta x+1)(x+\gamma )}$ and $\frac{f_d(x)}{\left(\gamma -\beta f_d{(x)}^2\right)}$ are positive and monotonically decreasing in $x\in [{\widehat{x}}_d,x_d]$. Therefore, ${\tilde{\alpha }}_d(x_d)\le {\tilde{\alpha }}_d({\widehat{x}}_d)$.

Case 2: ${\widehat{x}}_d>x_d$. By the same argument as above, it holds that , , and . Thus both $\frac{(\gamma -\beta x^2)}{(\beta x+1)(x+\gamma )}$ and $\frac{f_d(x)}{\left(\gamma -\beta f_d{(x)}^2\right)}$ are negative and monotonically decreasing in $x\in [x_d,{\widehat{x}}_d]$, hence their product is positive and increasing in $x\in [x_d,{\widehat{x}}_d]$. Therefore, ${\tilde{\alpha }}_d(x_d)\le {\tilde{\alpha }}_d({\widehat{x}}_d)$. ∎

Lemma 12 is proved by combining Lemma 13 and 14. This completes the proof of Theorem 9.

Strong spatial mixing on regular trees.

As a byproduct of our analysis, we prove a strong spatial mixing theorem for regular trees. When the graphs $G$ itself is a regular tree. All vertices (except the root) has the same arity. And all $d$s (excerpt the one of the root) that appear in the proof are the same and equal that arity. Then the condition that the uniqueness holds on all infinite regular trees of degree up to $\mathrm{\Delta }$ can be replaced by the uniqueness on infinite $\mathrm{\Delta }$-regular tree.

Theorem 15.

For two-state anti-ferromagnetic spin systems, on any infinite $\mathrm{\Delta }$-regular tree the uniqueness implies the strong spatial mixing.

The same result can be obtained by combining the same theorem for the hardcore model [26] and anti-ferromagnetic Ising model [21] and translating the parameters of general two-state anti-ferromagnetic spin systems to these models (as discussed in[24, 21]). However, unlike the hardcore and the anti-ferromagnetic Ising models, for general two-state anti-ferromagnetic spin systems Theorem 15 itself is not sufficient to imply the strong spatial mixing on graphs of maximum degree at most $\mathrm{\Delta }$. This is discussed in Section 5.

4 Algorithmic implications

In this section we prove Theorem 2. That is, if an anti-ferromagnetic $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique then there exists an FPTAS for the partition function $Z(G)$ for any graph $G$ of maximum degree at most $\mathrm{\Delta }$, and in particular the universal uniqueness implies an FPTAS for arbitrary graph $G$.

It is well-known that $Z(G)$ can be computed from $p_v^{{\sigma }_{\mathrm{\Lambda }}}$ by the following standard procedure. Let $v_1,\mathrm{\dots },v_n$ enumerate the vertices in $G$. For $0\le i\le n$, let ${\sigma }_i$ be the configuration fixing the first $i$ vertices $v_1,\mathrm{\dots },v_i$ as follows: ${\sigma }_i(v_j)={\sigma }_{i-1}(v_j)$ for $1\le j\le i-1$ and ${\sigma }_i(v_i)$ is fixed so that $p_i\triangleq \mathrm{Pr}[{\sigma }_i(v_i)\mid {\sigma }_{i-1}]\ge 1/3$. In particular, ${\sigma }_n\in {\{0,1\}}^V$ is a configuration of $V$. It holds for the Gibbs measure of ${\sigma }_n$ that $\rho ({\sigma }_n)=p_1{p}_2\mathrm{\cdots }{p}_n$ as well as that $\rho ({\sigma }_n)=\frac{w({\sigma }_n)}{Z(G)}$, thus $Z(G)=\frac{w({\sigma }_n)}{p_1{p}_2\mathrm{\cdots }{p}_n}$, where the weight $w({\sigma }_n)={\prod }_{(u,v)\in E}{A}_{{\sigma }_n(u),{\sigma }_n(v)}{\prod }_{v\in V}{b}_{{\sigma }_n(v)}$ can be computed precisely for any particular ${\sigma }_n$ in time polynomial in $n$. Note that $p_i$ equals to either $p_{v_i}^{{\sigma }_{i-1}}$ or $1-p_{v_i}^{{\sigma }_{i-1}}$. Therefore, if $p_v^{{\sigma }_{\mathrm{\Lambda }}}$ can be approximated within an additive error $ϵ$ in time polynomial in $n$ and $\frac{1}{ϵ}$, then the configurations ${\sigma }_i$ can be efficiently constructed such that all $p_i$ are bounded away from 0, thus the product $p_1{p}_2\mathrm{\cdots }{p}_n$ can be approximated within a factor of $(1\pm nϵ)$ in time polynomial in $n$ and $\frac{1}{ϵ}$, which implies an FPTAS for $Z(G)$.

Bounded degree graphs.

Let $G$ be a graph whose maximum degree is at most $\mathrm{\Delta }$ and $v$ be any vertex. A self-avoiding walk tree $T=T_{\text{SAW}}(G,v)$ can be constructed so that $R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}=R_T^{{\sigma }_{\mathrm{\Lambda }}}$. We can use the recursive procedure described in Section 3 to compute the upper and lower bounds of $R_T^{{\sigma }_{\mathrm{\Lambda }}}$, with the setting that for all the vertices more than $t$ steps away from the root $v$, the trivial bounds $0\le R_T^{{\sigma }_{\mathrm{\Lambda }}}\le \mathrm{\infty }$ is used. Then the proof of Theorem 9 shows that the recursive procedure returns $R_0$ and $R_1$ such that $R_0\le R_T^{{\sigma }_{\mathrm{\Lambda }}}\le R_1$, and $R_1-R_0=O({\alpha }^t)$ for some constant assuming that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique. Note that $R_T^{{\sigma }_{\mathrm{\Lambda }}}=R_{G,v}^{{\sigma }_{\mathrm{\Lambda }}}=\frac{p_v^{{\sigma }_{\mathrm{\Lambda }}}}{1-p_v^{{\sigma }_{\mathrm{\Lambda }}}}$. Let $p_0=\frac{R_0}{R_0+1}$ and $p_1=\frac{R_1}{R_1+1}$. Then $p_0\le p_v^{{\sigma }_{\mathrm{\Lambda }}}\le p_1$ and

$p_1-p_0={\displaystyle \frac{R_1}{R_1+1}}-{\displaystyle \frac{R_0}{R_0+1}}\le R_1-R_0=O({\alpha }^t).$ (7)

The recursive procedure runs in time $O({\mathrm{\Delta }}^t)$ since it only needs to construct the first $t$ levels of the self-avoiding walk tree. If $\mathrm{\Delta }$ is bounded, by setting $t={\log }_{\alpha }ϵ$, this gives an algorithm which approximates $p_v^{{\sigma }_{\mathrm{\Lambda }}}$ within an additive error $O(ϵ)$ in time polynomial in $n$ and $\frac{1}{ϵ}$, which implies an FPTAS for $Z(G)$.

Arbitrary graphs.

Let $G$ be an arbitrary graph and $v$ be any vertex. Let $T=T_{\text{SAW}}(G,v)$. We use the method of Computationally Efficient Correlation Decay introduced in [15] to deal with the unbounded degrees. Intuitively, using this method we observe correlation decay in a refined metric instead of graph distance such that in this new metric a neighborhood of moderate size is sufficient to guarantee desirable correlation decay.

Similarly, we use the recursive procedure described in Section 3 to compute the upper and lower bounds of $R_T^{{\sigma }_{\mathrm{\Lambda }}}$, but this time the termination condition relies on a new depth defined as follows.

Definition 16.

Let $T$ be a rooted tree and $M>1$ be a constant. For any vertex $v$ in $T$, define the $M$-based depth of $v$, denoted ${\mathrm{\ell }}_M(v)$, as such: ${\mathrm{\ell }}_M(v)=0$ if $v$ is the root, and ${\mathrm{\ell }}_M(v)={\mathrm{\ell }}_M(u)+\lceil {\log }_M(d+1)\rceil$ if $v$ is one of the $d$ children of $u$.

Let $M>1$ to be fixed. Denote by $B(\mathrm{\ell })$ the set of all vertices with $M$-based depth along with their children and grandchildren in $T$. It can be verified by induction that $|B(\mathrm{\ell })|\le n^2{M}^{\mathrm{\ell }}$. The recursion is applied to estimate the $R_T^{{\sigma }_{\mathrm{\Lambda }}}$ when the current $v\in B(\mathrm{\ell })$ until $v$ is no longer in $B(\mathrm{\ell })$ in which case the trivial bounds $0\le R_T^{{\sigma }_{\mathrm{\Lambda }}}\le \mathrm{\infty }$ is used.

Let ${ϵ}_v$ be defined as in Section 3. Repeatedly applying Item 3 of Lemma 11, without loss of generality, we have a path $u_1{u}_2\mathrm{\cdots }{u}_k$ in $T$ from the root $u_1=v$ to a $u_k$ with ${\mathrm{\ell }}_M(u_k)\ge \mathrm{\ell }$ and , such that ${ϵ}_{u_j}\le \alpha (d_j;x_1,\mathrm{\dots },x_{d_j})\cdot {ϵ}_{u_{j+1}}$ for $j=1,2,\mathrm{\dots },k$, where $d_j$ is the number of children of $u_j$ and $x_i\in [0,\mathrm{\infty })$, $1\le i\le d_j$. The key to overcome the explosion caused by unbounded degrees is to observe that the contraction $\alpha (d;x_1,\mathrm{\dots },x_d)$ decreases dramatically as the degree $d$ grows.

Lemma 17.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. If $(\beta ,\gamma ,\lambda )$ satisfies the universal uniqueness, then there exist constants and $M>1$ such that for any integer $d\ge 1$, and any $x_i\in [0,\mathrm{\infty })$, $1\le i\le d$, it holds that $\alpha (d;x_1,\mathrm{\dots },x_d)\le {\alpha }^{\lceil {\log }_M(d+1)\rceil }$.

Proof.

Assume the universal uniqueness of $(\beta ,\gamma ,\lambda )$. Due to Lemma 12, there exists a constant such that $\alpha (d;x_1,\mathrm{\dots },x_d)\le \alpha$. By Item 2 of Lemma 21, the universal uniqueness implies that $\gamma >1$, thus there exists a constant $M>1$ such that $d\cdot \sqrt{\frac{\lambda }{{\gamma }^{d+1}}}\le {\alpha }^{\lceil {\log }_M(d+1)\rceil }$ for all $d\ge M$. Due to (3), it holds that $\alpha (d;x_1,\mathrm{\dots },x_d)\le d\cdot \sqrt{\frac{\lambda }{{\gamma }^{d+1}}}\le {\alpha }^{\lceil {\log }_M(d+1)\rceil }$ for all $d\ge M$. Note that $\alpha (d;x_1,\mathrm{\dots },x_d)\le \alpha$ means that $\alpha (d;x_1,\mathrm{\dots },x_d)\le {\alpha }^{\lceil {\log }_M(d+1)\rceil }$ for . Therefore, $\alpha (d;x_1,\mathrm{\dots },x_d)\le {\alpha }^{\lceil {\log }_M(d+1)\rceil }$ for all $d$. ∎

By Lemma 17, there exist constants and $M>1$ such that

${ϵ}_v$ $\le {ϵ}_{u_k}\cdot {\displaystyle \prod _{j=1}^k}{\alpha }^{\lceil {\log }_M(d_j+1)\rceil }\le {ϵ}_{u_k}\cdot {\alpha }^{{\sum }_{j=1}^k\lceil {\log }_M(d_j+1)\rceil }={ϵ}_{u_k}\cdot {\alpha }^{{\mathrm{\ell }}_M(u_k)}\le {ϵ}_{u_k}\cdot {\alpha }^{\mathrm{\ell }}.$

With the notation used in Section 3, $S$ is the complement of $B(\mathrm{\ell })$. Note that all $u_k$’s children are in $B(\mathrm{\ell })$ thus none of them are in $S$, and by Item 2 of Lemma 21 the universal uniqueness implies that $\gamma >1$. Thus by Item 2 of Lemma 11 it holds that ${ϵ}_{u_k}=O(1)$. Therefore, ${ϵ}_v\le {ϵ}_{u_k}\cdot {\alpha }^{\mathrm{\ell }}=O({\alpha }^{\mathrm{\ell }})$.

Let ${\delta }_v=R_1-R_0$, where $R_0$ and $R_1$ are the bounds returned by the recursive procedure such that $R_0\le R_T^{{\sigma }_{\mathrm{\Lambda }}}\le R_1$. By the same analysis as in Section 3, ${\delta }_v=\frac{{ϵ}_v}{\mathrm{\Phi }(\tilde{R})}=O({ϵ}_v)=O({\alpha }^{\mathrm{\ell }})$. Then by (7), the marginal probability $p_v^{{\sigma }_{\mathrm{\Lambda }}}$ is approximated within an additive error of $O({\alpha }^{\mathrm{\ell }})$. The running time of the recursion is $O(nB(\mathrm{\ell }))=O(n^3{M}^{\mathrm{\ell }})$. By setting $\mathrm{\ell }={\log }_{\alpha }ϵ$, we have an algorithm which approximates $p_v^{{\sigma }_{\mathrm{\Lambda }}}$ within an additive error of $O(ϵ)$ in time polynomial in $n$ and $\frac{1}{ϵ}$, which implies an FPTAS for $Z(G)$ for arbitrary graph $G$.

Heterogeneous spin systems.

Our analysis in last and this sections actually holds for heterogeneous spin systems which allow that each vertex $v$ has a distinct constant external field ${\lambda }_v>0$.

Theorem 18.

For a two-state anti-ferromagnetic heterogeneous spin system with parameters $\beta$, $\gamma$, and external field ${\lambda }_v$ at each vertex $v$, for any finite $\mathrm{\Delta }\ge 2$ or $\mathrm{\Delta }=\mathrm{\infty }$, if for all $v$ the $(\beta ,\gamma ,{\lambda }_v)$ is up-to-$\mathrm{\Delta }$ unique then the spin system is of strong spatial mixing and has FPTAS on graphs of maximum degree at most $\mathrm{\Delta }$.

5 Non-monotonicity of general two-state spin system

In this section, we prove the following theorem.

Theorem 19.

There exist two-state anti-ferromagnetic spin systems which exhibit strong spatial mixing on infinite $d$-regular tree but does not exhibit strong spatial mixing on all graphs of maximum degree at most $d$.

In a seminal work [26], Weitz proved that for the hardcore model the strong spatial mixing on an infinite $d$-regular tree implies the strong spatial mixing on graphs of maximum degree at most $d$ (Theorem 2.3 in [26]). He further remarked that the same implication holds for all two-state spin systems. That is,

for any two-state spin system, strong spatial mixing on an infinite $d$-regular tree implies the strong spatial mixing on graphs of maximum degree at most $d$.

This claim played important roles in current understanding of correlation decay in two-state spin systems as well as devising FPTAS for such systems. An algorithmic form of this claim was cited in [21] as a theorem for all two-state spin systems (Theorem 2.8 in [21]) and was proved for the anti-ferromagnetic Ising model. It was raised as a conjecture in [22] whether this claim holds for multi-state spin systems.

Here we clarify that this claim holds only for the two-state spin systems under limited settings but does not hold for all general two-state spin systems. This disproves the conjecture in [22] and shows that the common belief that the $d$-regular tree represents the worst case for strong spatial mixing among all graphs of maximum degree at most $d$ cannot be generalized to general two-state spin systems. We first describe a region that the claim is true.

Lemma 20.

For $0\le \beta ,\gamma \le 1$, the strong spatial mixing on infinite $d$-regular tree implies the strong spatial mixing on trees of maximum degree at most $d$.

Proof.

Given a rooted tree of maximum degree at most $d$, for each vertex of $k$ children with , we can attach $d-1-k$ dummy children with fixed (distributions of) spin states. This is the method used in [26] and [21]: for the hardcore model the dummy children are fixed to be unoccupied and for the anti-ferromagnetic Ising model the dummy children are of uniform distributions over spin states. In both cases, the dummy children have no effect on their parent. In general, we fix the distribution to be $(p_0,p_1)$ at each dummy child satisfying

$p_0+p_1$ $=1,$

$\beta p_0+p_1$ $=p_0+\gamma p_1.$

When $0\le \beta ,\gamma \le 1$, this system has solutions in $0\le p_0,p_1\le 1$. With the ratio $R_i$ at the $i$-th child, $1\le i\le k$, and $R_i=p_0/p_1$ for the dummy children , the ratio at the parent is given by the recursion

$\lambda {\displaystyle \prod _{i=1}^{d-1}}{\displaystyle \frac{\beta R_i+1}{R_i+\gamma }}=\lambda {\displaystyle \prod _{i=1}^k}{\displaystyle \frac{\beta R_i+1}{R_i+\gamma }},$

which is identical to the original quantity. ∎

Due to the self-avoiding walk tree construction (Theorem 6), it holds that for $0\le \beta ,\gamma \le 1$, strong spatial mixing on infinite $d$-regular tree implies the strong spatial mixing and the FPTAS on graphs of maximum degree at most $d$.

For $0\le \beta ,\gamma \le 1$, the spin system shows the following monotone property: the uniqueness on infinite $d$-regular tree implies the uniqueness on all infinite regular trees of smaller degree. This can be verified by the following reasoning: due to Theorem 15, on infinite $d$-regular tree the uniqueness implies the strong spatial mixing, which for $0\le \beta ,\gamma \le 1$, implies the strong spatial mixing (including the uniqueness) on all infinite regular trees of smaller degree.

There exist two-state anti-ferromagnetic spin systems which are non-monotone. We can choose anti-ferromagnetic $(\beta ,\gamma ,\lambda )$ satisfying that $\gamma >1$ and $\lambda >{\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )$, where ${\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )$ is the upper threshold for universal uniqueness given in Item 8 of Lemma 21. Due to Item 3 of Lemma 21, for $\gamma >1$ the uniqueness holds on $d$-regular trees for all sufficiently large $d$. On the other hand, due to Item 8 of Lemma 21, $(\beta ,\gamma ,\lambda )$ cannot be universally unique when $\lambda >{\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )$, which means that there exists a finite $d^{\prime }$ such that the system is non-unique on $d^{\prime }$-regular tree.

For such non-monotone systems, due to Theorem 15, the uniqueness implies the strong spatial mixing on $d$-regular tree for sufficiently large $d$, but the strong spatial mixing does not hold on a regular tree of smaller degree (because of the non-uniqueness on the smaller tree). Therefore the implication between the strong spatial mixing on $d$-regular tree and on graphs of maximum degree at most $d$ does not hold for general two-state spin systems.

Acknowledgment.

We would like to thank Alistair Sinclair and Piyush Srivastava for the helpful comments to an early version of this paper.

Acknowledgment for the 2021 revision.

We would like to thank Xiaoyu Chen, Weiming Feng and Xinyuan Zhang for locating an error in Lemma 21 in the original (2013) version of the paper and suggesting its fix.

Appendix A The Uniqueness Thresholds

The following lemma translates the uniqueness condition into its various threshold forms. The lemma (except for item 2) is not used in the proofs of the main results but is used in the interpretation of the main results and comparisons to the previous results which are mostly stated in threshold forms.

Lemma 21.

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic.

1. $(\beta ,\gamma ,\lambda )$ is up-to-2 unique.

2. If $\gamma \le 1$, then the uniqueness does not hold on infinite $d$-regular tree for all sufficiently large $d$.

3. If $\gamma >1$, then the uniqueness holds on infinite $d$-regular tree for all sufficiently large $d$.

4. For any $\mathrm{\Delta }$ (including $\mathrm{\Delta }=\mathrm{\infty }$), there exists a critical threshold ${\gamma }_c={\gamma }_c(\beta ,\lambda ,\mathrm{\Delta })$ such that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\gamma \in ({\gamma }_c,\frac{1}{\beta })$. In particular, ${\gamma }_c(\beta ,\lambda ,\mathrm{\infty })>1$ and ${\gamma }_c(\beta ,\lambda ,\mathrm{\infty })={\gamma }_c(\beta ,\lambda ,\mathrm{\Delta })$ for some finite $\mathrm{\Delta }$.

5. If $\beta =0$, for any $\mathrm{\Delta }$ (including $\mathrm{\Delta }=\mathrm{\infty }$), there exists a critical threshold such that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\lambda \in (0,{\lambda }_c)$.

6. If $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, then for any external field $\lambda$, $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique.

7. If $\beta >0$, for any $\mathrm{\Delta }$ (including $\mathrm{\Delta }=\mathrm{\infty }$) that $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$, the regime for the $\lambda$ satisfying up-to-$\mathrm{\Delta }$ uniqueness is as follows.

   Let $\overline{\mathrm{\Delta }}\triangleq \frac{1+\sqrt{\beta \gamma }}{1-\sqrt{\beta \gamma }}$. For any integer , let $x_1(d)\le x_2(d)$ be the two positive roots of equation $\frac{d(1-\beta \gamma )x}{(\beta x+1)(x+\gamma )}=1$, more specifically,

   $x_1(d)={\displaystyle \frac{-1-\beta \gamma +d(1-\beta \gamma )-\sqrt{{(-1-\beta \gamma +d(1-\beta \gamma ))}^2-4\beta \gamma }}{2\beta }},$

   $x_2(d)={\displaystyle \frac{-1-\beta \gamma +d(1-\beta \gamma )+\sqrt{{(-1-\beta \gamma +d(1-\beta \gamma ))}^2-4\beta \gamma }}{2\beta }}.$

   Let ${\lambda }_i(d)\triangleq x_i(d){\left(\frac{x_i(d)+\gamma }{\beta x_i(d)+1}\right)}^d$, where $i=1,2$, be defined for integers .

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

   (8)

   And if furthermore $\gamma \le 1$, then $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\lambda \in (0,{\lambda }_c)\cup ({\overline{\lambda }}_c,\mathrm{\infty })$, where

   ${\lambda }_c$

   ${\overline{\lambda }}_c$

   In particular, when $\beta =\gamma$, it holds that ${\lambda }_c\cdot {\overline{\lambda }}_c=1$, and thus $(\beta ,\beta ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $|\log \lambda |>\log {\overline{\lambda }}_c$.

8. $(\beta ,\gamma ,\lambda )$ can be universally unique only when $\gamma >1$. If $\gamma >1$, there exist finite positive constants ${\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}={\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}(\beta ,\gamma )$ and ${\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}={\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )$ where ${\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}\le {\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}$, such that $(\beta ,\gamma ,\lambda )$ is universally unique if and $(\beta ,\gamma ,\lambda )$ is universally unique only if .

Remark (An error in the SODA’13 version).

In the earlier version of the paper, there was an error in Item 7 of Lemma 21 (which was Lemma 3.1 in the SODA’13 conference version [16]). In that version, the uniqueness regime (8) was unconditionally and wrongly simplified to the bi-interval simple form $(0,{\lambda }_c)\cup ({\overline{\lambda }}_c,\mathrm{\infty })$, which actually holds more restrictively when both $\beta ,\gamma \le 1$. However, when $\gamma >1$, there exist such that the uniqueness regime (8) may become more complicated than consisting of at most two intervals. For example, when $\beta =\frac{1}{122}$, $\gamma =100$, and $\mathrm{\Delta }=23$, the regime for the $\lambda$ such that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique, consists of three intervals $(0,102.664)\cup (104.339,106.967)\cup (109.444,\mathrm{\infty })$.

The current version corrects the error and changes the statements of Item 7 and Item 8 of Lemma 21. We further remark that the error does not affect the main result of the paper, namely the implication from the up-to-$\mathrm{\Delta }$ uniqueness to strong spatial mixing and FPTAS, because Lemma 21 is not used in the proofs of the main results (except for Item 2), but is only used for the interpretation of the uniqueness regime in various threshold forms and comparison with other known results.

Proof of Lemma 21.

Let $f_d(x)=\lambda {\left(\frac{\beta x_d+1}{x_d+\gamma }\right)}^d$ and ${\widehat{x}}_d=f_d({\widehat{x}}_d)$ be the positive fixed point of $f_d(x)$. Then

$|f_d({\widehat{x}}_d)|={\displaystyle \frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}}.$

Let $(\beta ,\gamma ,\lambda )$ be anti-ferromagnetic. That is, $0\le \beta \le \gamma$, $\gamma >0$, and , thus .

1. It is easy to verify that for , $(\beta x+1)(x+\gamma )-(1-\beta \gamma )x>0$ for any positive $x>0$. Therefore, when $d=1$, we have , which means that $(\beta ,\gamma ,\lambda )$ is up-to-2 unique.

2. For all sufficiently large $d$, it holds that

   $\lambda {\beta }^d\exp \left({\displaystyle \frac{d}{(1-\beta \gamma )d-3}}\right)$

   $\lambda \exp \left(-{\displaystyle \frac{\gamma d}{d(1-\beta \gamma )-3}}\right)$ $>{\displaystyle \frac{\gamma }{d(1-\beta \gamma )-3}}.$

   By contradiction, suppose that $|f_d({\widehat{x}}_d)|=\frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}\le 1$. Then,

   $1\ge {\displaystyle \frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}}={\displaystyle \frac{d(1-\beta \gamma )}{\beta {\widehat{x}}_d+\frac{\gamma }{{\widehat{x}}_d}+(1+\beta \gamma )}}\ge {\displaystyle \frac{d(1-\beta \gamma )}{\beta {\widehat{x}}_d+\frac{\gamma }{{\widehat{x}}_d}+2}}.$

   • ${\widehat{x}}_d\ge \gamma$. Then $\frac{\gamma }{{\widehat{x}}_d}\le 1$. Thus,

     $1\ge {\displaystyle \frac{d(1-\beta \gamma )}{\beta {\widehat{x}}_d+\frac{\gamma }{{\widehat{x}}_d}+2}}\ge {\displaystyle \frac{d(1-\beta \gamma )}{\beta {\widehat{x}}_d+3}},$

     which implies that ${\widehat{x}}_d\ge \frac{d(1-\beta \gamma )-3}{\beta }$. However, it holds that

     ${\widehat{x}}_d$ $=\lambda {\left({\displaystyle \frac{\beta {\widehat{x}}_d+1}{{\widehat{x}}_d+\gamma }}\right)}^d\le \lambda {\left({\displaystyle \frac{\beta {\widehat{x}}_d+1}{{\widehat{x}}_d}}\right)}^d\le \lambda {\left(\beta +{\displaystyle \frac{\beta }{d(1-\beta \gamma )-3}}\right)}^d$

     a contradiction.

   • . Then . Thus,

     $1\ge {\displaystyle \frac{d(1-\beta \gamma )}{\beta {\widehat{x}}_d+\frac{\gamma }{{\widehat{x}}_d}+2}}\ge {\displaystyle \frac{d(1-\beta \gamma )}{\frac{\gamma }{{\widehat{x}}_d}+3}},$

     which implies that ${\widehat{x}}_d\le \frac{\gamma }{d(1-\beta \gamma )-3}$. However, it holds that

     ${\widehat{x}}_d$ $=\lambda {\left({\displaystyle \frac{\beta {\widehat{x}}_d+1}{{\widehat{x}}_d+\gamma }}\right)}^d\ge {\displaystyle \frac{\lambda }{{({\widehat{x}}_d+1)}^d}}\ge \lambda {\left(1+{\displaystyle \frac{\gamma }{d(1-\beta \gamma )-3}}\right)}^{-d}$

     $\ge \lambda \cdot \exp \left(-{\displaystyle \frac{\gamma d}{d(1-\beta \gamma )-3}}\right)>{\displaystyle \frac{\gamma }{d(1-\beta \gamma )-3}},$

     a contradiction.

3. Let $\gamma >1$. The fixed point ${\widehat{x}}_d=\lambda {\left(\frac{\beta {\widehat{x}}_d+1}{{\widehat{x}}_d+\gamma }\right)}^d\le \frac{\lambda }{{\gamma }^d}$, and

   $|f_d^{\prime }({\widehat{x}}_d)|={\displaystyle \frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}}\le {\displaystyle \frac{d\lambda }{{\gamma }^d}},$

   which is strictly less than 1 for all sufficiently large $d$. Thus the uniqueness holds on infinite $d$-regular tree for all sufficiently large $d$.

4. We first show that there exists a critical threshold ${\gamma }_c={\gamma }_c(\beta ,\lambda ,\mathrm{\Delta })$ such that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\gamma \in ({\gamma }_c,\frac{1}{\beta })$. It is sufficient to show that if an anti-ferromagnetic $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique then $(\beta ,{\gamma }^{\prime },\lambda )$ is up-to-$\mathrm{\Delta }$ unique for any ${\gamma }^{\prime }>\gamma$ and .

   Recall that ${\widehat{x}}_d$ is the positive fixed point of $f_d(x)=\lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d$. Also let ${\widehat{x}}_d^{\prime }$ denote the positive solution to $x=\lambda {\left(\frac{\beta x+1}{x+{\gamma }^{\prime }}\right)}^d$.

   We first show that . By contradiction, assume that ${\widehat{x}}_d^{\prime }\ge {\widehat{x}}_d$. Since for anti-ferromagnetic $(\beta ,\gamma ,\lambda )$, $f_d(x)$ is monotonically decreasing, we have

   $${\widehat{x}}_d=\lambda {\left(\frac{\beta {\widehat{x}}_d+1}{{\widehat{x}}_d+\gamma }\right)}^d\ge \lambda {\left(\frac{\beta {\widehat{x}}_d^{\prime }+1}{{\widehat{x}}_d^{\prime }+\gamma }\right)}^d>\lambda {\left(\frac{\beta {\widehat{x}}_d^{\prime }+1}{{\widehat{x}}_d^{\prime }+{\gamma }^{\prime }}\right)}^d={\widehat{x}}_d^{\prime },$$

   a contradiction.

   Since

   we have

   For , it also holds that

   Multiplying the above two inequalities together, we have

   Note that these are the absolute derivatives at the respective fixed points when the parameters are $(\beta ,\gamma ,\lambda )$ and $(\beta ,{\gamma }^{\prime },\lambda )$. Therefore if $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique then $(\beta ,{\gamma }^{\prime },\lambda )$ is up-to-$\mathrm{\Delta }$ unique.

   Due to Part 2 of the lemma, if $\gamma \le 1$, $|f_d^{\prime }({\widehat{x}}_d)|>1$ for all sufficiently large $d$, thus ${\gamma }_c(\beta ,\lambda ,\mathrm{\infty })>1$. And due to Part 3 of the lemma, for any $\gamma \ge {\gamma }_c(\beta ,\lambda ,\mathrm{\infty })>1$, $|f_d^{\prime }({\widehat{x}}_d)|$ is arbitrarily close to 0 as $d$ grows to infinity, thus ${\gamma }_c(\beta ,\lambda ,\mathrm{\infty })={\gamma }_c(\beta ,\lambda ,\mathrm{\Delta })$ for a finite $\mathrm{\Delta }$.

5. When $\beta =0$, $|f_d^{\prime }({\widehat{x}}_d)|=\frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}=\frac{d{\widehat{x}}_d}{{\widehat{x}}_d+\gamma }$, the uniqueness condition is equivalent to that (here we assume $d>1$ since due to Part 1 of the lemma, for $d=1$ the uniqueness always holds). Recall that ${\widehat{x}}_d=\lambda {\left(\frac{1}{{\widehat{x}}_d+\gamma }\right)}^d$. Then if and only if

   Let . It holds that $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\lambda \in (0,{\lambda }_c)$.

6. We note that $|f_d^{\prime }({\widehat{x}}_d)|=\frac{d(1-\beta \gamma ){\widehat{x}}_d}{(\beta {\widehat{x}}_d+1)({\widehat{x}}_d+\gamma )}$ is not monotone in ${\widehat{x}}_d$. It achieves its maximum value at ${\widehat{x}}_d=\sqrt{\frac{\beta }{\gamma }}$. Therefore, if for any , for $x=\sqrt{\frac{\beta }{\gamma }}$, then $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique for any $\lambda$. This condition holds when $\sqrt{\beta \gamma }>\frac{d-1}{d+1}$ for all , i.e. when $\sqrt{\beta \gamma }>\frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$.

7. Let $\theta (d)\triangleq -1-\beta \gamma +d(1-\beta \gamma )$. It holds that $\theta (d)\ge 2\sqrt{\beta \gamma }$ for all $d\ge \overline{\mathrm{\Delta }}=\frac{1+\sqrt{\beta \gamma }}{1-\sqrt{\beta \gamma }}$. Thus, $x_1(d),x_2(d)$ are well-defined and can be expressed as:

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

   Recall that ${\lambda }_i(d)\triangleq x_i(d){\left(\frac{x_i(d)+\gamma }{\beta x_i(d)+1}\right)}^d$ for $i=1,2$.

   The following monotonicity and unboundedness of ${\lambda }_2(d)$ is easy to verify.

   Claim 22.

   ${\lambda }_2(d)$ is monotonically increasing in $d$ and goes to infinity as $d$ grows.

   Proof.

   Observe that for $0\le \beta \le \gamma$, and $d\ge \overline{\mathrm{\Delta }}$, it holds that and $\frac{x+\gamma }{\beta x+1}$ is increasing in $x$. Moreover,

   ${\displaystyle \frac{x_2(d)+\gamma }{\beta x_2(d)+1}}$ $=$ ${\displaystyle \frac{1}{\beta }}\cdot {\displaystyle \frac{(d-1)(1-\beta \gamma )+\sqrt{\theta {(d)}^2-4\beta \gamma }}{(d+1)(1-\beta \gamma )+\sqrt{\theta {(d)}^2-4\beta \gamma }}}$

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

   It is easy to verify that $x_2(d)$ is increasing in $d$ and unbounded as $d$ grows. Therefore, we can conclude that ${\lambda }_2(d)=x_2(d){\left(\frac{x_2(d)+\gamma }{\beta x_2(d)+1}\right)}^d$ is increasing in $d$ and is unbounded as $d$ grows. ∎

   We proceed to analyze the uniqueness regime. For , it holds that $\sqrt{\beta \gamma }>\frac{d-1}{d+1}$, in which case it always holds that due to Part 6 of this lemma.

   Recall that $\sqrt{\beta \gamma }\le \frac{\mathrm{\Delta }-2}{\mathrm{\Delta }}$. It remains to analyze the $|f_d^{\prime }({\widehat{x}}_d)|$ for those . For such $d$’s, we have $\sqrt{\beta \gamma }\le \frac{d-1}{d+1}$, and since $x_1(d)$ and $x_2(d)$ are the roots of equation $\frac{d(1-\beta \gamma )x}{(\beta x+1)(x+\gamma )}=1$, it holds that if and only if or ${\widehat{x}}_d>x_2(d)$. Note that $x{\left(\frac{x+\gamma }{\beta x+1}\right)}^d$ is monotonically increasing in $x$. Thus:

   ${\widehat{x}}_d>x_2(d)$ $⟺\lambda ={\widehat{x}}_d{\left({\displaystyle \frac{{\widehat{x}}_d+\gamma }{\beta {\widehat{x}}_d+1}}\right)}^d>x_2(d){\left({\displaystyle \frac{x_2(d)+\gamma }{\beta x_2(d)+1}}\right)}^d={\lambda }_2(d),$

   which means for any , if and only if $\lambda \in (0,{\lambda }_1(d))\cup ({\lambda }_2(d),\mathrm{\infty })$.

   Therefore, $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\lambda$ belongs to the regime

   which is exactly the uniqueness regime (8). Recall that

   where the last equation is due to the monotonicity of ${\lambda }_2(d)$.

   We then show that for $\gamma \le 1$, the uniqueness regime (8) becomes $(0,{\lambda }_c)\cup ({\overline{\lambda }}_c,\mathrm{\infty })$. We need the following claims, whose proofs are postponed to the end of this section.

   Claim 23.

   For , it always holds that

   furthermore, if ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all integers , then

   Claim 24.

   The followings hold:

   • if $\gamma \le 1$, then ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all integers $d\ge \overline{\mathrm{\Delta }}$;

   • if $\gamma >1$, there is a finite $d_0=d_0(\beta ,\gamma )\ge \overline{\mathrm{\Delta }}$ such that ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all $d\ge d_0$.

   Now assume that $\gamma \le 1$. By Claim 24, ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all integers . Then due to Claim 23, the uniqueness regime (8) becomes

   In particular, consider $\beta =\gamma$, the anti-ferromagnetic Ising model. It is easy to verify that $x_1(d)\cdot x_2(d)=\frac{\gamma }{\beta }=1$, and thus

   $${\lambda }_1(d)\cdot {\lambda }_2(d)=x_1(d)\cdot x_2(d){\left(\frac{x_1(d)+\gamma }{\beta x_1(d)+1}\cdot \frac{x_2(d)+\gamma }{\beta x_2(d)+1}\right)}^d=1.$$

   Therefore, ${\lambda }_1(d)=1/{\lambda }_2(d)$ is monotonically decreasing, due to the monotonicity of ${\lambda }_2(d)$ guaranteed by Claim 22. As a result, ${\lambda }_c={\lambda }_1(\mathrm{\Delta }-1)$, ${\overline{\lambda }}_c={\lambda }_2(\mathrm{\Delta }-1)$, and ${\lambda }_c\cdot {\overline{\lambda }}_c={\lambda }_1(\mathrm{\Delta }-1)\cdot {\lambda }_2(\mathrm{\Delta }-1)=1$. Then $(\beta ,\beta ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $|\log \lambda |>\log {\overline{\lambda }}_c$.

8. The case of $\beta =0$ is taken cared of in Part 5 of this lemma. In such case, due to Part 5 of the lemma, $(0,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\lambda \in (0,{\lambda }_c(\gamma ,\mathrm{\Delta }))$ where . In particular, ${\lambda }_c(\gamma ,\mathrm{\infty })={\min }_{d>1}\frac{{\gamma }^{d+1}{d}^d}{{(d-1)}^{d+1}}$ becomes 0 if $\gamma \le 1$ and becomes a positive and finite ${\lambda }_c(\gamma )>0$ if $\gamma >1$. Therefore, $(0,\gamma ,\lambda )$ is universally unique only if $\gamma >1$, and when $\gamma >1$, $(0,\gamma ,\lambda )$ is universally unique if and only if $\lambda \in (0,{\lambda }_c(\gamma ))$. We can let ${\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}(\gamma )={\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\gamma )={\lambda }_c(\gamma )$ and the lemma follows.

   For the case of $\beta >0$, due to Part 7 of the lemma, for all $\mathrm{\Delta }\ge \overline{\mathrm{\Delta }}\triangleq \frac{1+\sqrt{\beta \gamma }}{1-\sqrt{\beta \gamma }}$, the tuple $(\beta ,\gamma ,\lambda )$ is up-to-$\mathrm{\Delta }$ unique if and only if $\lambda$ belongs to the following regime:

   When $\gamma \le 1$, this regime becomes $(0,{\lambda }_c(\beta ,\gamma ,\mathrm{\Delta }))\cup (\overline{\lambda }(\beta ,\gamma ,\mathrm{\Delta }),\mathrm{\infty })$, where and .

   Claim 25.

   The followings hold:

   • if $\gamma \le 1$, then ${\lambda }_1(d)$ approaches 0 as $d$ grows to infinity;

   • if $\gamma >1$, then ${\lambda }_1(d)$ is unbounded as $d$ grows to infinity.

   Proof.

   For integers $d\ge \overline{\mathrm{\Delta }}$, let $\theta (d)\triangleq -1-\beta \gamma +d(1-\beta \gamma )$, and let

   ${\widehat{\lambda }}_1(d)\triangleq {\displaystyle \frac{2{\gamma }^{d+1}}{\theta (d)}}{\left({\displaystyle \frac{d+1}{d-1}}\right)}^d.$ (9)

   When $\gamma \le 1$, it is easy to verify that ${\widehat{\lambda }}_1(d)$ approaches 0 as $d$ grows to infinity.

   Moreover, ${\lambda }_1(d)$ is upper bounded by ${\widehat{\lambda }}_1(d)$. Specifically, for $d\ge \overline{\mathrm{\Delta }}$,

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

   Since $\frac{x+\gamma }{\beta x+1}$ is increasing in $x$ when , we have

   $$\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{d+1}{d-1}.$$

   Thus, for all $d\ge \overline{\mathrm{\Delta }}$,

   ${\lambda }_1(d)$ $=x_1(d){\left({\displaystyle \frac{x_1(d)+\gamma }{\beta x_1(d)+1}}\right)}^d\le {\displaystyle \frac{2{\gamma }^{d+1}}{\theta (d)}}{\left({\displaystyle \frac{d+1}{d-1}}\right)}^d={\widehat{\lambda }}_1(d).$ (10)

   Therefore, when $\gamma \le 1$, ${\lambda }_1(d)$ approaches 0 as $d$ grows to infinity.

   On the other hand, when $\gamma >1$, for $d\ge \overline{\mathrm{\Delta }}$,

   ${\lambda }_1(d)=x_1(d){\left({\displaystyle \frac{x_1(d)+\gamma }{\beta x_1(d)+1}}\right)}^d\ge x_1(d)\cdot {\gamma }^d={\displaystyle \frac{2{\gamma }^{d+1}}{\theta (d)+\sqrt{\theta {(d)}^2-4\beta \gamma }}},$

   which is obviously unbounded as $d$ grows. ∎

   When $\lambda \le 1$, due to above claim, ${\lambda }_1(d)$ approaches 0 as $d$ grows to infinity, which means ${\lambda }_c(\beta ,\gamma ,\mathrm{\infty })=0$; and by Claim 22, ${\lambda }_2(d)$ goes to infinity as $d$ grows to infinity, which means ${\overline{\lambda }}_c(\beta ,\gamma ,\mathrm{\infty })=\mathrm{\infty }$. Therefore, $(\beta ,\gamma ,\lambda )$ can be universally unique only when $\gamma >1$.

   Now assume $\gamma >1$. By Claim 25, ${\lambda }_1(d)$ is unbounded as $d$ grows. Thus, ${\lambda }_1(d)$ achieves its minimum value at a finite $d=d(\beta ,\gamma )\ge \overline{\mathrm{\Delta }}$. And it is obvious that ${\lambda }_1(d)$ is always finite and positive for any finite $d\ge \overline{\mathrm{\Delta }}$. Therefore, ${\lambda }_c(\beta ,\gamma ,\mathrm{\infty })={\lambda }_c(\beta ,\gamma ,d(\beta ,\gamma ))$ is finite and positive. Let ${\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}(\beta ,\gamma )={\lambda }_c(\beta ,\gamma ,\mathrm{\infty })>0$ denote this finite positive threshold.

   On the other hand, by Claim 24, for $\gamma >1$, there is a finite $d_0=d_0(\beta ,\gamma )\ge \overline{\mathrm{\Delta }}$ such that ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all $d\ge d_0$. Let ${\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )={\min }_{d\ge d_0}{\lambda }_1(d)$. Obviously, ${\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )\ge {\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}(\beta ,\gamma )>0$. And ${\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )\le {\lambda }_1(d_0(\beta ,\gamma ))$ is also finite.

   Due to Claim 23 and the monotonicity and unboundedness of ${\lambda }_2(d)$ by Claim 22, the regime for $\mathrm{\infty }$-uniqueness, namely ${\bigcap }_{d⩾\overline{\mathrm{\Delta }}}\left((0,{\lambda }_1(d))\cup ({\lambda }_2(d),\mathrm{\infty })\right)$, is sandwiched as follows:

   $(0,{\lambda }_c^{\mathrm{𝗅𝗈𝗐𝖾𝗋}}(\beta ,\gamma ))$ $\subseteq {\displaystyle \bigcap _{d⩾\overline{\mathrm{\Delta }}}}\left((0,{\lambda }_1(d))\cup ({\lambda }_2(d),\mathrm{\infty })\right)$

   $\subseteq {\displaystyle \bigcap _{d⩾d_0}}\left((0,{\lambda }_1(d))\cup ({\lambda }_2(d),\mathrm{\infty })\right)$

   $=(0,{\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma ))\cup ({\lambda }_2(\mathrm{\infty }),\mathrm{\infty })$

   $=(0,{\lambda }_c^{\mathrm{𝗎𝗉𝗉𝖾𝗋}}(\beta ,\gamma )).$

   Therefore, $(\beta ,\gamma ,\lambda )$ is universally unique, if and only if .

∎

It now remains to prove Claim 23 and Claim 24.

Proof of Claim 23.

First, it is easy to verify that for , it always holds that

since for all , it always holds that

Next, we prove that the two regimes are equal if ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all integers . This is proved by induction.

The equality in the claim is trivial to hold when , because in this case there is at most one integer $d$ satisfying . Then both sides of the equality becomes

$$\left((0,{\lambda }_1(d))\cup ({\lambda }_2(d),\mathrm{\infty })\right).$$

This establishes the induction basis.

For the induction hypothesis, assume that the equality holds for $d_1>d_0$:

(11)

where we denote

For the induction step, we want to establish the equality when $d_1$ is increased by 1. Thus in addition, we assume that ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for the unique integer $d\in [d_1-1,d_1)$, that is

${\lambda }_1(\lceil d_1\rceil )\le {\lambda }_2(\lceil d_1\rceil -1).$ (12)

To finish the induction, it is then sufficient to verify that

assuming the induction hypothesis (11). Therefore, it is sufficient to verify that

$=$ (13)

Note that $\lceil d_1\rceil$ is the unique integer in $[d_1,d_1+1)$. And due to the monotonicity of ${\lambda }_2(d)$, we have . Therefore, , and thus (13) holds as long as , which is guaranteed by the assumption (12). ∎

Proof of Claim 24.

Recall function ${\widehat{\lambda }}_1(d)$ defined in (9): For integers $d\ge \overline{\mathrm{\Delta }}$, let ${\widehat{\lambda }}_1(d)\triangleq \frac{2{\gamma }^{d+1}}{\theta (d)}{\left(\frac{d+1}{d-1}\right)}^d$, where $\theta (d)\triangleq -1-\beta \gamma +d(1-\beta \gamma )$. By (10), for all $d\ge \overline{\mathrm{\Delta }}$, we have ${\lambda }_1(d)\le {\widehat{\lambda }}_1(d)$.

It also holds that ${\widehat{\lambda }}_1(d)\le {\lambda }_2(d)$. For all $d\ge \overline{\mathrm{\Delta }}$, we have $x_2(d)=\frac{\theta (d)+\sqrt{\theta {(d)}^2-4\beta \gamma }}{2\beta }\ge \frac{\theta (d)}{2\beta }$. Since $\frac{x+\gamma }{\beta x+1}$ is monotonically increasing in $x$ when , it holds that $\frac{x_1(d)+\gamma }{\beta x_1(d)+1}\ge \frac{\frac{\theta (d)}{2\beta }+\gamma }{\frac{\theta (d)}{2}+1}=\frac{1}{\beta }\cdot \frac{d-1}{d+1}$. Meanwhile, for all $d\ge \overline{\mathrm{\Delta }}=\frac{1+\sqrt{\beta \gamma }}{1-\sqrt{\beta \gamma }}$, it holds that $\beta \gamma \le {\left(\frac{d-1}{d+1}\right)}^2$ and $\theta {(d)}^2\ge 4\beta \gamma$. Therefore,

${\lambda }_2(d)$ $=x_2(d){\left({\displaystyle \frac{x_2(d)+\gamma }{\beta x_2(d)+1}}\right)}^d\ge {\displaystyle \frac{\theta (d)}{2{\beta }^{d+1}}}{\left({\displaystyle \frac{d-1}{d+1}}\right)}^d\ge {\displaystyle \frac{2{\gamma }^{d+1}}{\theta (d)}}{\left({\displaystyle \frac{d+1}{d-1}}\right)}^d={\widehat{\lambda }}_1(d).$ (14)

Assume $\gamma \le 1$. It is also not hard to observe that ${\widehat{\lambda }}_1(d)$ is decreasing in $d\ge \overline{\mathrm{\Delta }}>1$. Indeed, if $\gamma \le 1$, then $\frac{2{\gamma }^{d+1}}{\theta (d)}>0$ is monotonically decreasing in $d\ge \overline{\mathrm{\Delta }}$; and we can verify the monotonicity of function $h(d)\triangleq {\left(\frac{d+1}{d-1}\right)}^d$ for $d\ge \overline{\mathrm{\Delta }}>1$ as follows:

${\left(\ln h(d)\right)}^{\prime }$ $=-{\displaystyle \frac{2d}{d^2-1}}+\ln \left({\displaystyle \frac{d+1}{d-1}}\right)\to 0\text{as}d\to \mathrm{\infty },$

${\left(\ln h(d)\right)}^{\prime \prime }$ $={\displaystyle \frac{4}{d^2-1}}>0,$

which guarantees that $h(d)$ is monotonically decreasing in $d$, and hence ${\widehat{\lambda }}_1(d)=\frac{2{\gamma }^{d+1}}{\theta (d)}\cdot h(d)$ is decreasing in $d$ for $d\ge \overline{\mathrm{\Delta }}$, since both $\frac{2{\gamma }^{d+1}}{\theta (d)}$ and $h(d)$ are positive and decreasing in $d$.

Therefore, when $\gamma \le 1$, due to (10), the monotonicity of ${\widehat{\lambda }}_1(d)$, and (14), for all $d\ge \overline{\mathrm{\Delta }}$,

$${\lambda }_1(d+1)\le {\widehat{\lambda }}_1(d+1)\le {\widehat{\lambda }}_1(d)\le {\lambda }_2(d).$$

This proves that ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all $d\ge \overline{\mathrm{\Delta }}$ when $\gamma \le 1$.

Finally, when $\gamma >1$, by contradiction suppose that ${\lambda }_1(d+1)>{\lambda }_2(d)$ for some $d\ge \overline{\mathrm{\Delta }}$. We then deduce a finite upper bound $d_0(\beta ,\gamma )$ on such bad $d$. By (10), we have ${\lambda }_1(d+1)\le {\widehat{\lambda }}_1(d+1)=\frac{2{\gamma }^{d+2}}{\theta (d+1)}\cdot h(d+1)\le \frac{2{\gamma }^{d+2}}{\theta (d)}\cdot h(d)$, where the last inequality is due to that $\theta (d)$ is increasing and $h(d)$ is decreasing; and by (14), we have ${\lambda }_2(d)\ge \frac{\theta (d)}{2{\beta }^{d+1}}\cdot \frac{1}{h(d)}$. Hence ${\lambda }_1(d+1)>{\lambda }_2(d)$ implies that

(15)

Recall that for $d\ge \overline{\mathrm{\Delta }}=\frac{1+\sqrt{\beta \gamma }}{1-\sqrt{\beta \gamma }}$, it always holds that $\theta {(d)}^2\ge 4\beta \gamma$. And without loss of generality, we only consider large enough $d\ge 2$, and thus $h(d)\le h(2)=9$ since $h(d)$ is decreasing. Therefore inequality (15) implies that

which gives us a necessary condition for ${\lambda }_1(d+1)>{\lambda }_2(d)$, among the $d$’s that $d\ge \overline{\mathrm{\Delta }}$ and $d\ge 2$. Let

$$d_0(\beta ,\gamma )\triangleq \max \{2,\overline{\mathrm{\Delta }},{\log }_{\frac{1}{\beta \gamma }}(81\gamma )\}.$$

It then must hold that ${\lambda }_1(d+1)\le {\lambda }_2(d)$ for all $d\ge d_0(\beta ,\gamma )$. ∎

Appendix B Proof of Lemma 11 (a Mean Value Theorem approach)

We prove Lemma 11. The notations are defined in Section 3.

1. Due to the Mean Value Theorem, there exists an $\tilde{R}\in [R_v,R_v+{\delta }_v]$ such that

   ${ϵ}_v$ $=\phi (R_v+{\delta }_v)-\phi (R_v)={\phi }^{\prime }(\tilde{R})\cdot {\delta }_v=\mathrm{\Phi }(\tilde{R})\cdot {\delta }_v.$

2. Suppose that each vertex has at most $k$ children. Then due to the assumption either $k$ is bounded or $\gamma >1$.

   If $v\notin \mathrm{\Lambda }$, then ${\delta }_v\le R_v+{\delta }_v=f(R_{v_1},\mathrm{\dots },R_{v_d})\le f(0,\mathrm{\dots },0)=\frac{\lambda }{{\gamma }^d}$, where $0\le d\le k$. If $\gamma >1$, then $R_v+{\delta }_v\le \lambda =O(1)$; and if $k$ is finite, then $R_v+{\delta }_v\le \max \{\frac{\lambda }{{\gamma }^k},\lambda \}=O(1)$.

   Due to the Mean Value Theorem, there exist $\widetilde{R_i}\in [R_{v_i},R_{v_i}+{\delta }_{v_i}]$, $1\le i\le d$, such that

   ${ϵ}_v$ $=\phi \left(f(R_{v_1},\mathrm{\dots },R_{v_d})\right)-\phi \left(f(R_{v_1}+{\delta }_{v_1},\mathrm{\dots },R_{v_d}+{\delta }_{v_d})\right)$

   $=-\nabla \phi \left(f(\widetilde{R_1},\mathrm{\dots },\widetilde{R_d})\right)\cdot ({\delta }_{v_1},\mathrm{\dots },{\delta }_{v_d})$

   $={\displaystyle \frac{(1-\beta \gamma )\cdot {\left(f(\widetilde{R_1},\mathrm{\dots },\widetilde{R_d})\right)}^{\frac{1}{2}}}{{\left(\beta f(\widetilde{R_1},\mathrm{\dots },\widetilde{R_d})+1\right)}^{\frac{1}{2}}{\left(f(\widetilde{R_1},\mathrm{\dots },\widetilde{R_d})+\gamma \right)}^{\frac{1}{2}}}}\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{{\delta }_{v_i}}{(\beta \widetilde{R_i}+1)(\widetilde{R_i}+\gamma )}}$

   $\le \sqrt{{\displaystyle \frac{f(0,\mathrm{\dots },0)}{\gamma }}}\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{{\delta }_{v_i}}{\gamma }}.$

   If $v_i\notin S$ for all $v_i$, then due to the above argument, ${\delta }_{v_i}\le \frac{\lambda }{{\gamma }^{d_i}}$ for free $v_i$, where $0\le d_i\le k$ is the number of children of $v_i$; and trivially ${\delta }_{v_i}=0$ for fixed $v_i$. Therefore, ${ϵ}_v\le \sqrt{\frac{\lambda }{{\gamma }^{d+1}}}{\sum }_{i=1}^d\frac{\lambda }{{\gamma }^{d_i+1}}\le {\lambda }^{\frac{3}{2}}d{\gamma }^{-\frac{d+3}{2}}\cdot \max \{\frac{1}{{\gamma }^k},1\}$. If $\gamma >1$, then it is easy to see that $d{\gamma }^{-\frac{d+3}{2}}$ is bounded by a constant, thus it holds that ${ϵ}_v=O(1)$; if $k$ is bounded, then $d\le k$ is also bounded, thus clearly ${ϵ}_v=O(1)$.

3. We then analyze the stepwise contraction of ${ϵ}_v$. Define that $y_v=\phi (R_v)$ and accordingly $y_{v_i}=\phi (R_{v_i})$, $1\le i\le d$. Then $y_v+{ϵ}_v=\phi (R_v+{\delta }_v)$ and $y_{v_i}+{ϵ}_{v_i}=\phi (R_{v_i}+{\delta }_{v_i})$, $1\le i\le d$. We have

   $y_v$ $=$ $\phi (f({\phi }^{-1}(y_{v_1}+{ϵ}_{v_1}),\mathrm{\dots },{\phi }^{-1}(y_{v_d}+{ϵ}_{v_d}))),$

   $y_v+{ϵ}_v$ $=$ $\phi (f({\phi }^{-1}(y_{v_1}),\mathrm{\dots },{\phi }^{-1}(y_{v_d}))).$

   Apply the Mean Value Theorem. There exist $\widetilde{y_i}\in [y_{v_i},y_{v_i}+{ϵ}_{v_i}]$ and corresponding $\widetilde{R_i}\in [R_{v_i},R_{v_i}+{\delta }_{v_i}]$ satisfying $\widetilde{y_i}=\phi (\widetilde{R_i}),1\le i\le d$, such that

   ${ϵ}_v$ $=\phi (f({\phi }^{-1}(y_{v_1}),\mathrm{\dots },{\phi }^{-1}(y_{v_d})))-\phi (f({\mathrm{\Phi }}^{-1}(y_{v_1}+{ϵ}_{v_1}),\mathrm{\dots },{\phi }^{-1}(y_{v_d}+{ϵ}_{v_d})))$

   $=-\nabla \phi (f({\phi }^{-1}(\widetilde{y_1}),\mathrm{\dots },{\phi }^{-1}(\widetilde{y_d})))\cdot ({ϵ}_{v_1},\mathrm{\dots },{ϵ}_{v_d})$

   $=-\mathrm{\Phi }(f(\widetilde{R_1},\mathrm{\dots },\widetilde{R_d}))\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{\partial f}{\partial R_i}}{\displaystyle \frac{1}{\mathrm{\Phi }(\widetilde{R_i})}}\cdot {ϵ}_i$

   $={\displaystyle \frac{(1-\beta \gamma ){\left(\lambda {\prod }_{i=1}^d\frac{\beta \widetilde{R_i}+1}{\widetilde{R_i}+\gamma }\right)}^{\frac{1}{2}}}{{\left(\beta \lambda {\prod }_{i=1}^d\frac{\beta \widetilde{R_i}+1}{\widetilde{R_i}+\gamma }+1\right)}^{\frac{1}{2}}{\left(\lambda {\prod }_{i=1}^d\frac{\beta \widetilde{R_i}+1}{\widetilde{R_i}+\gamma }+\gamma \right)}^{\frac{1}{2}}}}\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{{ϵ}_{v_i}{\widetilde{R_i}}^{\frac{1}{2}}}{{(\beta \widetilde{R_i}+1)}^{\frac{1}{2}}{(\widetilde{R_i}+\gamma )}^{\frac{1}{2}}}}$

   $\le \underset{1\le i\le d}{\max }\{{ϵ}_{v_i}\}\cdot {\displaystyle \frac{(1-\beta \gamma ){\left(\lambda {\prod }_{i=1}^d\frac{\beta \widetilde{R_i}+1}{\widetilde{R_i}+\gamma }\right)}^{\frac{1}{2}}}{{\left(\beta \lambda {\prod }_{i=1}^d\frac{\beta \widetilde{R_i}+1}{\widetilde{R_i}+\gamma }+1\right)}^{\frac{1}{2}}{\left(\lambda {\prod }_{i=1}^d\frac{\beta \widetilde{R_i}+1}{\widetilde{R_i}+\gamma }+\gamma \right)}^{\frac{1}{2}}}}\cdot {\displaystyle \sum _{i=1}^d}{\displaystyle \frac{{\widetilde{R_i}}^{\frac{1}{2}}}{{(\beta \widetilde{R_i}+1)}^{\frac{1}{2}}{(\widetilde{R_i}+\gamma )}^{\frac{1}{2}}}}$

   $=\alpha (d;{\tilde{R}}_1,\mathrm{\dots },{\tilde{R}}_d)\cdot \underset{1\le i\le d}{\max }\{{ϵ}_{v_i}\}.$

Appendix C Finding a good potential function heuristically

Perhaps the most mysterious step in our proof is the choice of the potential function $\mathrm{\Phi }(x)=\frac{1}{\sqrt{x(\beta x+1)(x+\gamma )}}$. As in many cases where potential analysis is applied, there is no standard routine for searching for a suitable potential function. On the other hand, it is quite unlikely that we can just guess such a fairly complicated formula without any hints. Here we present a heuristic approach which leads us to the discovery of a good potential function. This part is not rigorous and logically unnecessary for the soundness of our result. Nevertheless, this heuristic approach is general and interesting enough and may find its applications in other scenarios, thus deserves an exposition here.

The heuristics consists of three steps:

1. Find a necessary condition for the potential function, which is an equation related to the potential function at one point.

2. (heuristic step) Enhance the condition by assuming that the equation holds for the whole range of the variable, which gives us a differential equation.

3. Solve the differential equation and get a potential function.

We first assume that the system is at the boundary of uniqueness. This means that $d$ is the critical degree and $f^{\prime }(\widehat{x})=-1$, where $f(x)=\lambda {\left(\frac{\beta x+1}{x+\gamma }\right)}^d$ and $\widehat{x}=f(\widehat{x})$ is the positive fixed point.

We have the following identities:

$\widehat{x}=\lambda {\left({\displaystyle \frac{\beta \widehat{x}+1}{\widehat{x}+\gamma }}\right)}^d\text{and}1=\lambda d(1-\beta \gamma ){\displaystyle \frac{{(\beta \widehat{x}+1)}^{d-1}}{{(\widehat{x}+\gamma )}^{d+1}}},$

which together implies another identity

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

The goal is to find a potential function such that $\alpha (x)=|f^{\prime }(x)|\frac{\mathrm{\Phi }\left(f(x)\right)}{\mathrm{\Phi }(x)}\le 1$ for all $x$. On the other hand, we have $\alpha (\widehat{x})=|f^{\prime }(\widehat{x})|\frac{\mathrm{\Phi }(f(\widehat{x}))}{\mathrm{\Phi }(\widehat{x})}=1\cdot \frac{\mathrm{\Phi }(\widehat{x})}{\mathrm{\Phi }(\widehat{x})}=1$. So $\alpha (x)$ achieves its maximum when $x=\widehat{x}$. As a differentiable function, it holds that ${\alpha }^{\prime }(\widehat{x})=0$, that is

$${{\left[f^{\prime }(x)\cdot \frac{\mathrm{\Phi }\left(f(x)\right)}{\mathrm{\Phi }(x)}\right]}^{\prime }|}_{x=\widehat{x}}=0.$$

We have the following calculation:

${\left[f^{\prime }(x)\cdot {\displaystyle \frac{\mathrm{\Phi }(f(x))}{\mathrm{\Phi }(x)}}\right]}_{x=\widehat{x}}^{\prime }=0$

$\iff$ ${{\left[f^{\prime }(x)\mathrm{\Phi }(f(x))\right]}^{\prime }\mathrm{\Phi }(x)|}_{x=\widehat{x}}={f^{\prime }(x)\mathrm{\Phi }(f(x)){\mathrm{\Phi }}^{\prime }(x)|}_{x=\widehat{x}}$

$\iff$ $\left[f^{\prime \prime }(\widehat{x})\mathrm{\Phi }(f(\widehat{x}))+f^{\prime }(\widehat{x}){\mathrm{\Phi }}^{\prime }(f(\widehat{x}))f^{\prime }(\widehat{x})\right]\mathrm{\Phi }(\widehat{x})=f^{\prime }(\widehat{x})\mathrm{\Phi }(f(\widehat{x})){\mathrm{\Phi }}^{\prime }(\widehat{x})$

$\iff$ $f^{\prime \prime }(\widehat{x})\mathrm{\Phi }(\widehat{x})+{\mathrm{\Phi }}^{\prime }(\widehat{x})=-{\mathrm{\Phi }}^{\prime }(\widehat{x})$

$\iff$ $-{\displaystyle \frac{f^{\prime \prime }(\widehat{x})}{2}}={\displaystyle \frac{{\mathrm{\Phi }}^{\prime }(\widehat{x})}{\mathrm{\Phi }(\widehat{x})}}={\left(\ln (\mathrm{\Phi }(\widehat{x}))\right)}^{\prime },$

where we use the fact that $\widehat{x}=f(\widehat{x})$ and $f^{\prime }(\widehat{x})=-1$.

Taking the second derivative of $f(x)$, we have

$f^{\prime \prime }(x)=\lambda d(\beta \gamma -1){\displaystyle \frac{{(\beta x+1)}^{d-2}}{{(x+\gamma )}^{d+2}}}\left((d-1)\beta (x+\gamma )-(d+1)(\beta x+1)\right).$

This equation already gives an identity for $\mathrm{\Phi }$. But it is too complicated. We may further simplify the above formula at the point $x=\widehat{x}$.

$f^{\prime \prime }(\widehat{x})$ $=$ $\lambda d(\beta \gamma -1){\displaystyle \frac{{(\beta \widehat{x}+1)}^{d-2}}{{(\widehat{x}+\gamma )}^{d+2}}}\left((d-1)\beta (\widehat{x}+\gamma )-(d+1)(\beta \widehat{x}+1)\right)$

$=$ ${\displaystyle \frac{1}{(\beta \widehat{x}+1)(\widehat{x}+\gamma )}}\left((d+1)(\beta \widehat{x}+1)-(d-1)\beta (\widehat{x}+\gamma )\right)$

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

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

$=$ ${\displaystyle \frac{1}{\widehat{x}}}+{\displaystyle \frac{1}{\widehat{x}+\gamma }}+{\displaystyle \frac{\beta }{\beta \widehat{x}+1}}.$

So we have

${\left(\ln \left(\mathrm{\Phi }(\widehat{x})\right)\right)}^{\prime }=-{\displaystyle \frac{f^{\prime \prime }(\widehat{x})}{2}}=-{\displaystyle \frac{1}{2}}({\displaystyle \frac{1}{\widehat{x}}}+{\displaystyle \frac{1}{\widehat{x}+\gamma }}+{\displaystyle \frac{\beta }{\beta \widehat{x}+1}}).$ (17)

We apply the heuristics and assume that (17) simply holds for all $x$. This gives us a differential equation

$${\left(\ln (\mathrm{\Phi }(x))\right)}^{\prime }=-\frac{1}{2}(\frac{1}{x}+\frac{1}{x+\gamma }+\frac{\beta }{\beta x+1}).$$

The solution of this differential equation is

$$\ln (\mathrm{\Phi }(x))=-\frac{1}{2}\ln (x(x+\gamma )(\beta x+1))+C_1,$$

which gives that

$$\mathrm{\Phi }(x)=\frac{C_2}{\sqrt{x(\beta x+1)(x+\gamma )}},$$

where $C_1,C_2$ are some constants. This gives the potential function we used in the paper.

There are also other equations which hold for the fixed point $\widehat{x}$. Choosing which equation for $\widehat{x}$ to heuristically extend to all $x$ may affect the potential function we obtained. For example, (16) implies that

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

thus we can rewrite $f^{\prime \prime }(\widehat{x})$ as

$$f^{\prime \prime }(\widehat{x})=\frac{1}{\widehat{x}}+\frac{1}{\widehat{x}+\gamma }+\frac{\beta }{\beta \widehat{x}+1}=\frac{d+1}{d\widehat{x}}+\frac{2\beta }{\beta \widehat{x}+1}.$$

This gives that

$${\left(\ln (\mathrm{\Phi }(\widehat{x}))\right)}^{\prime }=-\frac{d+1}{2d\widehat{x}}-\frac{\beta }{\beta \widehat{x}+1}.$$

We can also treat this equation as a differential equation for variable $x$ and solving it gives us

$$\mathrm{\Phi }(x)=\frac{c}{x^{\frac{d+1}{2d}}(\beta x+1)},$$

which is the potential function used in [15].