A Cluster Expansion and the Decay of Correlations of the 1d Long-Range Ising Model at low temperatures


Abstract

In this work, a convergent low-temperature cluster expansion of the one-dimensional long-range ferromagnetic Ising model with polynomial decay \(\alpha\in (1,2]\) is developed; that is, \(J(r)=r^{-\alpha}\). As an application, the \(n\)-point correlations are studied and the two-point correlation is shown to be algebraic with rate of decay exactly \(\alpha\).

1 Introduction↩︎

1.1 A brief history of the one-dimensional long-range Ising model↩︎

In 1967, Gallavotti and Miracle-Sole [1] showed the existence of the thermodynamic limit and infinite-volume correlation functions for a large class of lattice models with spin space \(\{0,1\}\) and translation-invariant potentials. In particular, they showed the existence of the thermodynamic limit and infinite-volume correlation functions for regular long-range ferromagnetic Ising models. A year later, Ruelle [2] showed the absence of spontaneous magnetization if the long-range interaction decayed too fast, that is, if \(\alpha>2\). In 1969, employing a hierarchical model, Dyson [3] gave the first proof of the existence of phase transition for a ferromagnet with long-range polynomial decay \(\alpha\in (1,2)\). In the following years, progressively stronger criteria for the lack of long-range order were proved in [4], [5].

In 1980, Bricmont, Lebowitz and Pfister [6] showed that, under Ruelle’s condition, any equilibrium state is translation invariant. In 1982, Fannes, Vanheuverzwijn and Verbeure [7] established a criterion for translation invariance of equilibrium states in terms of correlation inequalities, which show, in particular, that any equilibrium state of a regular one-dimensional long-range ferromagnetic Ising model is translation invariant. A similar result (translation invariance of equilibrium states in the one-dimensional case) was obtained independently by Burkov and Sinai in [8].

The matter of phase-transition for critical decay remained open until [9]. In 1982, Fröhlich and Spencer successfully employed a new contour notion, which had already proved very fruitful in [10], where the first rigorous proof of the BKT transition for the planar XY model was given, to show the existence of long-range order at small temperature for \(\alpha=2\). In [11], building on the work of Fröhlich and Spencer, Imbrie provided a converging cluster expansion for the critical exponent and, as an application, proved the algebraic decay of the two-point function.

In 1988, employing techniques of percolation theory and the hypothesis that the interaction between first neighbors be arbitrarily large (which was absent in [9], [11]), Aizenman, Chayes, Chayes and Newman [12] proved the non-continuity of the magnetization for the critical polynomial decay (Thouless effect). In the same year and building on the previous work, Imbrie and Newman [13] proved the existence of an intermediary phase, where correlations decay algebraically with an exponent smaller than \(2\); that is, a BKT-like transition.

Nearly twenty years later, interest in the one-dimensional long-range ferromagnetic Ising model was renewed, when Cassandro, Ferrari, Merola and Presutti [14] gave a proof via contours of phase transition in case \(\alpha\in (3/2,2]\) (the same proof works for \(\alpha\in(3-\log_23,2]\)). Similarly to the works cited in the previous paragraph, they assumed that the interaction between first neighbors is sufficiently large. Under this hypothesis and employing the same contours, Cassandro, Orlandi and Picco [15] proved the stability of phase-transition under the influence of an external random field, provided \(\alpha\in(3-\log_23,3/2)\) and the variance \(\varepsilon\) be sufficiently small. Recently, this result was extended to \(\alpha\in (1,3/2)\) by Ding, Maia and Huang in [16] without the introduction of the first-neighbor interaction as a perturbative factor.

In [17], Cassandro, Merola, Picco and Rozikov provided a convergent cluster expansion, and as an application studied the phase separation point in the one-dimensional long-range Ising model (we recall the absence of non-translation invariant equilibrium states in 1D). They also enunciate a bound for the \(n\)-site correlation function (see Theorem 2.5 in [17]). In [18], the authors dealt with the question of phase separation for \(\alpha\in (3-\log_23,2)\), once again employing cluster expansion techniques.

A first attempt to eliminate the hypothesis that the interaction between first-neighbors is sufficiently large can be found in [19]. Bissacot, Endo, van Enter, Kimura and Ruszel were partially successful: employing the same contours as [14], [15], [17], [18], the authors proved the existence of phase-transition provided \(\alpha\in (\alpha^*,2]\), where \(\zeta(\alpha^*)=2\). We note that \(3-\log_23<\alpha^*\). Another work dealing with the same contours is [20], Littin and Picco show that the energy estimates break down for \(\alpha<3-\log_23\).

Motivated by the very fruitful adaptation of the Fröhlich-Spencer contours to the multidimensional case [21][26], Affonso, Bissacot, Corsini and Welsch [27] obtained, with appropriate modifications to energy estimates, a proof of phase transition via (Fröhlich-Spencer) contours for the region \(\alpha\in (1,2]\), without the introduction of the nearest neighbor interaction as a perturbative factor. This fact motivates revisiting, and complementing, whenever possible, many results where, while a contour approach was successful, the introduction of a perturbative first-neighbors interaction was necessary.

In this work, we provide a converging cluster expansion and, as an application, show that at low temperature the two-point correlation decays polynomially, with the same decay rate as \(J(r)=r^{-\alpha}\); the lower-bound was already proved in the seventies by Iagolnitzer and Souillard [28]. In order to study the \(n\)-point correlation function, we also revisit (abstractly) the technique of cluster expansion in the very particular scenario, where compatibility is (graph)-transitive.

On the one hand, summing over leaves of trees of polymers is very standard (see, for example, [29][34]). The convergence of the cluster expansion present in this work can be proved with the traditional criteria (Kotecký-Preiss [32], Dobrushin [29], [30], Fernández-Procacci [31]). For a comparison between traditional methods, see [35]. For recent technical improvements, see [36][38]. On the other hand, given the form of the activity of this particular gas of polymers, we sum over trees of contours; that is, trees where compatibility conditions are global. Hence, it is reasonable (and useful) to introduce more general operations which are represented graphically by removing vertices of an arbitrary degree.

1.2 Description of the results↩︎

In this work, we develop a convergent low-temperature cluster expansion of the one-dimensional long-range ferromagnetic Ising model with polynomial decay \(\alpha\in (1,2]\) in terms of polymer of Fröhlich-Spencer contours. In Section 2, we revisit the definition of the contours first introduced in [9] and list the crucial estimates that ensure the convergence of the cluster expansion.

In Section 3, we derive a representation of the one-dimensional Ising model in terms of a hard-core gas of polymers, which are collections of contours satisfying a "positive" compatibility relation. In this section, we also prove an upper bound to the activity of the polymer gas which is more convenient for computations.

In Section 4, we prove the usual "one-vertex" bound employing a convergence criterion similar to Kotecký-Preiss [32]; however, we note that our computation is closer to that of Fernández-Procacci [31]. We also prove (many-vertex) estimates for a special type of expansion in terms of trees where compatibility is global. This is useful because of the particular form of the upper bound to the activity of the polymer gas. In Section 5, we apply the abstract bounds present in Section 4 to our particular case. Transforming trees of polymers into trees of contours, we prove the convergence of the representation of the pressure given by the cluster expansion. That is, we prove the following theorem:

Theorem 1. If \(\alpha\in(1,2]\), then there exist \(\beta_0=\beta_0(\alpha)>0\) such that if \(\beta\geq\beta_0\), then \[\begin{align} \sup_{\Lambda\Subset\mathbb{Z}}\frac{1}{\#\Lambda}\sum_{n=1}^\infty\frac{1}{n!}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n_{\Lambda}}\lvert{\phi_n(\mathbf{\Gamma})}\rvert\prod_{i=1}^nz_{\beta}(\Gamma_i)<\infty \end{align}\] If \(A\subset\Lambda\) and \(\beta\geq\beta_0\), there exists \(r_0\mathrel{\vcenter{:}}= r_0(\beta,\lvert A\rvert)>0\) such that \[\begin{align} h\mapsto p_\beta(h)\mathrel{\vcenter{:}}=\lim_{\Lambda\uparrow\mathbb{Z}}\frac{1}{\#\Lambda}\sum_{n=1}^\infty\frac{1}{n!}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n_\Lambda}\phi_n(\mathbf{\Gamma})\prod_{i=1}^nz^+_{\beta,h}(\Gamma_i) \end{align}\] is analytic in the polydisk \(\mathbb{D}_{r_0}^A\mathrel{\vcenter{:}}=\{h\in\mathbb{C}^A:\lVert{h}\rVert_\infty\leq r_0\}\).

In Section 6, we establish some estimates in terms of the lattice sites in the spirit of [11]. Furthermore, we prove a lemma that changes estimates in terms of contours into estimates in terms of sites. In Section 7, we apply Theorem 1 and the contour-site estimates to obtain an upper bound to the decay rate of the correlation functions of the one-dimensional long-range ferromagnetic Ising model with polynomial decay \(\alpha\in (1,2]\). Thus, we obtain the following results:

Theorem 2. There exist \(c_3>0\) and \(\beta_0>0\) such that if \(\beta\geq\beta_0\), \(x\neq y\in\mathbb{Z}\) and \(\{x,y\}\subset\Lambda\), then \[\label{eq:32two-point-correlations} \rvert\langle\sigma_x;\sigma_y\rangle^+_{\Lambda,\beta}\rvert\leq 2e^{-c_3\beta}\lvert x-y\rvert^{-\alpha}.\tag{1}\]

Theorem 3. For all \(N\geq 3\), there exist \(\beta_0\mathrel{\vcenter{:}}=\beta_0(N)\) and \(c_4\mathrel{\vcenter{:}}= c_4(N)>0\) such that if \(\beta\geq\beta_0\), \(A\subset \Lambda\Subset \mathbb{Z}\) and \(\lvert A\rvert=N\), then \[\label{eq:32many-point-correlations} \lvert\langle:\sigma_A:\rangle^+_{\Lambda,\beta}\rvert\leq 2e^{-c_4\beta}\sum_{T\in\mathbbold{T}(A)}\prod_{\{a_1,a_2\}\in E(T)}\lvert a_1-a_2\rvert^{-\alpha}.\tag{2}\]

We note that combining Theorem 2 with [28], we have thus proven that, for \(\beta\) large, the two-point correlation decay rate is exactly \(\alpha\). We recall the existence of an intermediary phase in the case of \(\alpha=2\) [13].

2 Preliminaries↩︎

2.1 The long-range one-dimensional ferromagnetic Ising model↩︎

The configuration space of the one-dimensional Ising model is given by \(\Omega=\{-1,+1\}^\mathbb{Z}.\) Given \(\omega\in\Omega\) and \(\Lambda\Subset\mathbb{Z}\), we denote the set of configurations satisfying the \(\omega\)-boundary condition outside of \(\Lambda\) by \[\Omega^\omega_\Lambda\mathrel{\vcenter{:}}=\{\sigma\in\Omega:(x\notin\Lambda\Rightarrow\sigma_x=\omega_x)\}.\] We recall that \(\Lambda\Subset\mathbb{Z}\) denotes that \(\Lambda\) is a finite subset of \(\mathbb{Z}\). In a similar fashion, we denote the set of all configurations satisfying the \(\omega\)-boundary condition by \[\Omega^\omega\mathrel{\vcenter{:}}=\bigcup_{\Lambda\Subset\mathbb{Z}}\Omega^\omega_\Lambda.\]If \(\omega_x=+1\) for all \(x\in\mathbb{Z}\), we write \(\Omega^+\) instead of \(\Omega^\omega\) and say that \(\sigma\in\Omega^+\) satisfies the plus boundary condition. If \(\omega_x=-1\) for all \(x\in\mathbb{Z}\), we proceed in an analogous way. Henceforth, we restrict ourselves to dealing with homogeneous boundary conditions, that is, either the plus boundary condition or the minus boundary condition.

Let \(J:\mathbb{Z}\times\mathbb{Z}\rightarrow\mathbb{R}\) be a coupling function and \(\omega\in\Omega\) be a boundary condition, we define the Hamiltonian \(H^\omega_J:\Omega^\omega\rightarrow\mathbb{R}\) by \[\label{eq:32def-hamiltonian} H^\omega_J(\sigma)\mathrel{\vcenter{:}}=\frac{1}{2}\sum_{(x,y)}J_{xy}(\omega_x\omega_y-\sigma_x\sigma_y).\tag{3}\] In this work, we assume that\[\label{eq:32polynomial-decay-coupling} J_{xy}=\begin{cases} 0 &\textrm{if }x=y,\\ \lvert x-y\rvert^{-\alpha}&\textrm{otherwise,} \end{cases}\tag{4}\] where \(1<\alpha\leq2\). That is, we assume that the model is translation invariant, ferromagnetic, with a two-body regular essentially long-range interaction. Henceforth, we omit \(J\) from the notation for better writing.

Let \(h:\mathbb{Z}\rightarrow\mathbb{R}\) be an external field and \(\omega\in\Omega\) be a boundary condition, we define the Hamiltonian \(H^\omega_h:\Omega^\omega\rightarrow\mathbb{R}\) by \[\label{eq:32def-hamiltonian-ext-field} H^\omega_h(\sigma)\mathrel{\vcenter{:}}=\frac{1}{2}\sum_{(x,y)}J_{xy}(\omega_x\omega_y-\sigma_x\sigma_y)+\sum_{x}h_x(\omega_x-\sigma_x)\eqqcolon H^\omega(\sigma)+E^\omega_h(\sigma).\tag{5}\] Given \(\Lambda\Subset\mathbb{Z}\) and \(\beta>0\), we denote the partition function at the inverse temperature \(\beta\) by \[\label{eq:32def-partition-function} Z^\omega_{\Lambda,\beta}\mathrel{\vcenter{:}}=\sum_{\sigma\in\Omega^\omega_\Lambda}e^{-\beta H^\omega(\sigma)}.\tag{6}\] By the global spin-flip map \(\sigma\mapsto-\sigma\), it is easy to see that \(Z^-_{\Lambda,\beta}=Z^+_{\Lambda,\beta}\). If the external field \(h\) is not identically zero, we define \[Z^\omega_{\Lambda,\beta,h}\mathrel{\vcenter{:}}=\sum_{\sigma\in\Omega^\omega_\Lambda}e^{-\beta H_h^\omega(\sigma)}.\]We note that if \(h\not\equiv0\), then \(Z^-_{\Lambda,\beta,h}\neq Z^+_{\Lambda,\beta,h}\).

2.2 Configurations as collections of spin flips↩︎

In this work, following [9], we think of configurations of the one-dimensional Ising model as collections of spin flips. We associate a configuration \(\sigma\in\Omega\) to \[\partial\sigma\mathrel{\vcenter{:}}=\{e=\{x,y\}\in\mathbb{Z}^*:\sigma_x\neq\sigma_y\},\] where \(\mathbb{Z}^*\simeq\mathbb{Z}+1/2\) denotes the dual graph of \(\mathbb{Z}\). This application is two-to-one and invariant under the global spin-flip map \(\sigma\mapsto -\sigma\).

Since we shall be dealing with configurations with homogeneous boundary conditions, we define our configuration space as \[\label{eq:32configuration-space} \Omega^*\mathrel{\vcenter{:}}=\{\sigma^*\in\mathrm{P}(\mathbb{Z}^*):\lvert\sigma^*\rvert\textrm{ is even}\},\tag{7}\] where \(\mathrm{P}(X)\) denotes the set of all finite parts of \(X\). We define the canonical identification \(\sigma^+:\Omega^*\rightarrow\Omega^+\) by \[\label{eq:32identification-config-spaces} \sigma^+(\partial\sigma)=\sigma,\tag{8}\] for all \(\sigma\in\Omega^+\), that is, if \(\sigma\) satisfies the plus boundary condition. For the rest of this work, unless noted otherwise, we identify collections of spin flips with configurations satisfying the plus boundary condition.

Given \(\sigma^*\in\Omega^*\), we define its \(-\)-interior by \[\label{eq:32minus-interior} \mathrm{I}_-(\sigma^*)\mathrel{\vcenter{:}}=\{x\in\mathbb{Z}:\sigma^+(\sigma^*)=-1\}.\tag{9}\] We note that for any coupling function \[\label{eq:32dual-hamiltonian} H(\sigma^*)\mathrel{\vcenter{:}}=\sum_{\substack{x\in \mathrm{I}_-(\sigma^*)\\y\in\mathrm{I}_-(\sigma^*)^c}}2J_{xy}=H^+(\sigma^+),\tag{10}\] where \(\sigma^+\equiv\sigma^+(\sigma^*)\). In particular, this holds for the case 4 .

Remark 1. This geometric description (independent of the finite box considered) of the Hamiltonian justifies the definition present in Equation 3 . In any case, without this normalization, which is not detected by Gibbs’ measures, there is no hope that a cluster expansion would converge.

If \(\sigma^*_1,\sigma^*_2\in\Omega^*\) satisfy \(\mathrm{I}_-(\sigma_1^*)\cap\mathrm{I}_-(\sigma^*_2)=\emptyset\), it is also useful to define \[\label{eq:32definition-phi} \Phi(\sigma^*_1,\sigma^*_2)\mathrel{\vcenter{:}}=\sum_{\substack{x\in \mathrm{I}_-(\sigma^*_1)\\y\in\mathrm{I}_-(\sigma^*_2)}}4J_{xy}.\tag{11}\] Finally, if \(h\not\equiv0\) is an external field, we define \[\label{eq:32external-fiels-component-definition} E_h(\sigma^*)\mathrel{\vcenter{:}}=\sum_{x\in\mathrm{I}_-(\sigma^*)}2h_x=E_h^+(\sigma^+)=-E_h^-(\sigma^-)\tag{12}\] and \(H_h(\sigma^*)=H(\sigma^*)+E_h(\sigma^*)\).

2.3 Long-range contours↩︎

In [9], Fröhlich and Spencer showed the existence of a first-order phase transition for the critical exponent \(\alpha=2\) using the notion of a geometric contour that need not be connected, which they first introduced a year earlier in [10] to give the first rigorous proof of the BKT transition to the first-neighbor two-dimensional XY model. Such contours are defined by a fixed distancing exponent \(a=3/2\) and a linear parameter \(M>1\) whose value is chosen during the proof.

In [22], the authors managed to prove the existence of a first-order phase transition in the region \(\alpha>d\geq 2\), where \(d\) denotes the dimension, employing a modification to the definition originally given by Fröhlich and Spencer in which the choice of the distancing exponent itself depended on \((d,\alpha)\). It turns out that, at least in one dimension, this does not need to be the case. In [27] it was shown that the same contours employed by Fröhlich and Spencer are sufficient to prove the existence of a first-order phase transition in the whole region of polynomial decay where such a phenomenon occurs. In fact, any distancing exponent \(a\in(1,2)\) is sufficient.

We recall the slight modification to the definition of an irreducible contour introduced by Imbrie in [11], where, by employing a cluster expansion method, he established sharp estimates of the correlations in the case \(\alpha=2\).

Notation 1. Let \(n\in\mathbb{N}\), we write \([n]\mathrel{\vcenter{:}}=\{1,...,n\}\).

Notation 1. Let \(V\) be a finite set. We denote by \(\mathbbold{G}(V)\) the set of unordered connected graphs \(G\) such that \(V(G)=V\). Similarly, we denote by \(\mathbbold{T}(V)\) the set of all unordered connected trees \(T\) such that \(V(T)=V\). Finally, we denote the set of unordered partitions of \(V\) by \(\mathbbold{P}(V)\). If \(V=[n]\), we write \(\mathbbold{G}_n\), \(\mathbbold{T}_n\) and \(\mathbbold{P}_n\) instead.

Convention 1. If \(V\) is a finite set and \(P\in\mathbbold{P}(V)\), then \(\emptyset\notin P\). If \(\lvert P\rvert=1\), we say that \(P\) is the trivial partition of \(V\).

Definition 1 (\(M\)-contour). Fix \(M>1\). We say that \(\emptyset\neq\gamma\in\Omega^*\) is an \(M\)-irreducible contour if for any \(P\in\mathbbold{P}(\gamma)\cap\mathrm{P}(\Omega^*)\) not trivial, there exists \(\gamma_1\neq\gamma_2\in P\) such that\[\mathrm{dist}(\gamma_1,\gamma_2)\leq M(\min\{\mathrm{diam}(\gamma_1),\mathrm{diam}(\gamma_2)\})^{3/2}.\]We denote the set of all \(M\)-irreducible contours by \(\mathscr{C}_M\).

If \(M>1\) is clear from context, we omit \(M\) from the notation and say that \(\gamma\in\mathscr{C}\) is a contour. If \(\gamma_1,\gamma_2\in\mathscr{C}\) are not empty and \[\mathrm{dist}(\gamma_1,\gamma_2)>M(\min\{\mathrm{diam}(\gamma_1),\mathrm{diam}(\gamma_2)\})^{3/2},\]we say that \(\gamma_1\) and \(\gamma_2\) are compatible, denoted by \(\gamma_1\sim\gamma_2\). Note that \(\gamma\not\sim\gamma\).

Definition 2 (\(M\)-partition). Let \(\emptyset\neq\sigma^*\in\Omega^*\). We say that \(\Gamma\in \mathbbold{P}(\sigma^*)\cap\mathrm{P}(\mathscr{C})\) is an \(M\)-partition of \(\sigma^*\) if \[(\gamma_1,\gamma_2\in\Gamma)\wedge(\gamma_1\neq\gamma_2)\Rightarrow\gamma_1\sim\gamma_2.\]We denote the set of all \(M\)-partitions by \(\mathscr{E}\).

Following 10 , we define \(H_h:\mathscr{E}\rightarrow\mathbb{R}\) by \[\label{eq:32hamiltonian-partition} H_h(\Gamma)\mathrel{\vcenter{:}}= H_h\left(\sqcup_{\gamma\in\Gamma}\gamma\right).\tag{13}\]

Lemma 1. Let \(\emptyset\neq \sigma^*\in\Omega^*\), then there is a unique \(M\)-partition of \(\sigma^*\).

Proof. We refer the reader to Proposition 2.1 in [11]. ◻

Corollary 1. Let \(\Gamma_1,\Gamma_2\in\mathscr{E}\) such that \((\gamma_1\in\Gamma_1)\wedge(\gamma_2\in\Gamma_2)\Rightarrow\gamma_1\sim\gamma_2\), then \(\Gamma_1\sqcup\Gamma_2\in\mathscr{E}\).

Convention 2. We establish \(\emptyset\notin\mathscr{C}\) and \(\emptyset\notin\mathscr{E}\). The reason for this shall be clear when dealing with collections of compatible polymers.

Definition 3. If \(\gamma_1,\gamma_2\in\mathscr{C}\) satisfy \(\gamma_1\sim\gamma_2\) and \(\mathrm{I}_-(\gamma_1)\cap\mathrm{I}_-(\gamma_2)=\emptyset,\) we say that they are positively compatible, which we denote by \(\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_2.\)

If \(\Gamma_1,\Gamma_2\in\mathscr{E}\) are such that \((\gamma_1\in\Gamma_1)\wedge(\gamma_2\wedge\Gamma_2)\Rightarrow\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_2\), we define \[\label{eq:32Phi-definition-for-collections-of-contours} \Phi(\Gamma_1,\Gamma_2)\mathrel{\vcenter{:}}=\Phi(\sqcup_{\gamma_1\in\Gamma_1}\gamma_1,\sqcup_{\gamma_2\in\Gamma_2}\gamma_2).\tag{14}\] We say that \(\Gamma\in\mathscr{E}\) is positive if \[\gamma_1\neq\gamma_2\in\Gamma\Rightarrow\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_2.\]We denote the set of all positive \(M\)-partitions by \(\mathscr{P}\).

Figure 1: All three contours are compatible, but \gamma_1 is not positively compatible with \gamma_2. On the other hand, \gamma_3 is positively compatible with both \gamma_1 and \gamma_2. Hence, \Gamma=\{\gamma_1,\gamma_2,\gamma_3\} is not a positive collection.

Remark 1. For those who are colorblind, let us give an alternative description of Figure 1. Ordering spin flips so that \(i<j\) implies that \(b_i\) lies to the left of \(b_j\), we have \[\begin{align} &\gamma_1=\{b_1,b_4,b_7,b_8\}, &\gamma_2=\{b_2,b_3\},\\ &\gamma_3=\{b_5,b_6\}. \end{align}\]

Lemma 2. Let \(\Gamma_1,\Gamma_2\in\mathscr{P}\) such that \((\gamma_1\in\Gamma_1)\wedge(\gamma_2\in\Gamma_2)\Rightarrow\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_2\), then \(\Gamma_1\sqcup\Gamma_2\in\mathscr{P}\). In this case, we write \(\Gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\Gamma_2\).

Proof. This follows directly from the fact that compatibility is checked in pairs. That is, \(\Gamma\in\mathscr{E}\) if and only if \(\gamma_1\neq\gamma_2\in\Gamma\Rightarrow\gamma_1\sim\gamma_2.\) ◻

Notation 1. If \(\Lambda\Subset\mathbb{Z}\), we write \[\begin{align} &\mathscr{C}_\Lambda\mathrel{\vcenter{:}}=\{\gamma\in\mathscr{C}:\mathrm{I}_-(\gamma)\subset\Lambda\},\\ &\mathscr{E}_\Lambda\mathrel{\vcenter{:}}=\{\Gamma\in\mathscr{E}:\gamma\in\Gamma\Rightarrow\gamma\in\mathscr{C}_\Lambda\},\\ &\mathscr{P}_\Lambda\mathrel{\vcenter{:}}=\mathscr{P}\cap\mathscr{E}_\Lambda. \end{align}\]

Definition 4. Let \(\gamma_1,\gamma_2\in\mathscr{C}\). If \(\mathrm{I}_-(\gamma_1)\subset\mathrm{I}_-(\gamma_2)\), we write \(\gamma_1\trianglelefteq\gamma_2\). This defines a partial ordering in \(\mathscr{C}\) such that \[\gamma_1\sim\gamma_2\Rightarrow(\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_2)\vee(\gamma_1\trianglelefteq\gamma_2)\vee(\gamma_2\trianglelefteq\gamma_1).\]

The partial ordering \(\trianglelefteq\) defines a natural projection \(\partial:\mathscr{E}\rightarrow\mathscr{P}\) in the following way: \[\label{def:32positive-contours-of-Gamma} \partial\Gamma\mathrel{\vcenter{:}}=\{\gamma\in\Gamma:(\gamma'\in\Gamma)\wedge(\gamma\neq\gamma')\Rightarrow(\gamma'\trianglelefteq\gamma)\vee(\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma')\}.\tag{15}\] Note that \[\mathrm{I}_-(\Gamma)\mathrel{\vcenter{:}}=\mathrm{I}_-\left(\sqcup_{\gamma\in\Gamma}\gamma\right)\subset\bigsqcup_{\gamma\in\partial\Gamma}\mathrm{I}_-(\gamma).\]Furthermore, \(\partial\Gamma\) is the smallest subset of \(\Gamma\) with this property. Equality (\(\Gamma=\partial\Gamma\)) holds if and only if \(\Gamma\in\mathscr{P}\).

Figure 2: Take \Gamma as in Figure 1, then \partial\Gamma=\{\gamma_1,\gamma_3\}.

Remark 1. Recalling Definition 2.7 in [25], \(\partial\Gamma\) denotes the set of external contours of \(\Gamma\).

Definition 5. We say that \(\Gamma\in\mathscr{E}\) is internal to \(\gamma_0\in\mathscr{C}\), if \(\gamma\in\Gamma\Rightarrow(\gamma\sim\gamma_0)\wedge(\gamma\trianglelefteq\gamma_0)\), which we denote by \(\Gamma\trianglelefteq\gamma_0\). We denote the collection of all \(M\)-partitions internal to \(\gamma_0\) by \[\label{eq:32internal-collection} \mathscr{I}(\gamma_0)\mathrel{\vcenter{:}}=\{\Gamma\in\mathscr{E}:\Gamma\trianglelefteq\gamma_0\}.\tag{16}\] Since \(\gamma_0 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1\Rightarrow\mathscr{I}(\gamma_0)\cap\mathscr{I}(\gamma_1)=\emptyset\), this notion extends to \(\Gamma\in\mathscr{P}\) by \[\label{eq:32internal-collection-2} \mathscr{I}(\Gamma)\mathrel{\vcenter{:}}=\bigsqcup_{\gamma\in\Gamma}\mathscr{I}(\gamma).\tag{17}\]

Remark 1. If \(\Gamma\in\mathscr{P}\) and \(\Gamma_1\in\mathscr{I}(\Gamma)\), then there is a unique \(\gamma\in\Gamma\) such that \(\Gamma_1\in\mathscr{I}(\gamma)\).

Given \(\gamma_0\in\mathscr{C}\), we define \(\iota^*_{\gamma_0}:\mathscr{E}\rightarrow\mathscr{I}(\gamma_0)\sqcup\{\emptyset\}\) and \(\iota_{\gamma_0}:\mathscr{E}\rightarrow\mathscr{E}\) by\[\begin{align} &\iota^*_{\gamma_0}(\Gamma)\mathrel{\vcenter{:}}=\{\gamma\in\Gamma:(\gamma\sim\gamma_0)\wedge(\gamma\trianglelefteq\gamma_0)\}, &\iota_{\gamma_0}(\Gamma)\mathrel{\vcenter{:}}=\{\gamma_0\}\sqcup\iota_{\gamma_0}^*(\Gamma). \end{align}\] These definitions naturally extend to positive collections \(\Gamma_0\in\mathscr{P}\). We write \[\begin{align} &\iota^*_{\Gamma_0}(\Gamma)\mathrel{\vcenter{:}}=\bigsqcup_{\gamma_0\in\Gamma_0}\iota^*_{\gamma_0}(\Gamma), &\iota_{\Gamma_0}(\Gamma)\mathrel{\vcenter{:}}=\bigsqcup_{\gamma_0\in\Gamma_0}\iota_{\gamma_0}(\Gamma). \end{align}\]Finally, if \(\Gamma\in\mathscr{E}\), then \[\label{eq:32relation-minus-interior-Gamma-interior-gamma} \mathrm{I}_-(\Gamma)=\bigsqcup_{\gamma\in\partial\Gamma}\mathrm{I}_-(\iota_\gamma(\Gamma))=\mathrm{I}_-(\partial\Gamma)\cap \bigcap_{\gamma\in\partial\Gamma}\mathrm{I}_-(\iota^*_\gamma(\Gamma))^c.\tag{18}\]

Notation 1. Henceforth, we write (recall Equation 11 ) \[\begin{align} &H_h(\iota_{\gamma_0}(\Gamma))\equiv H_{h,\gamma_0}(\Gamma),\\ &H_h(\iota^*_{\gamma_0}(\Gamma))\equiv H^*_{h,\gamma_0}(\Gamma),\\ &\Phi\left(\iota_{\gamma_1}(\Gamma),\iota_{\gamma_2}(\Gamma)\right)\equiv\Phi_{\gamma_1,\gamma_2}(\Gamma). \end{align}\] If \(h\equiv0\), we omit it from the notation.

Lemma 3. If \(\Gamma\in\mathscr{E}\), then \[\label{eq:32energy95positive-partition} H_h(\Gamma)=\sum_{\gamma\in\partial\Gamma}H_{\gamma,h}(\Gamma)-\sum_{\substack{\gamma_1,\gamma_2\in\partial\Gamma\\\gamma_1\neq\gamma_2}}\frac{1}{2}\Phi_{\gamma_1,\gamma_2}(\Gamma).\tag{19}\] In particular, if \(\Gamma\in\mathscr{P}\), then \[\label{eq:32energy-positive-partition-2} H_h(\Gamma)=\sum_{\gamma\in\Gamma}H_h(\gamma)-\sum_{\substack{\gamma_1,\gamma_2\in\Gamma\\\gamma_1\neq\gamma_2}}\frac{1}{2}\Phi(\gamma_1,\gamma_2).\tag{20}\]

Proof. Let \(\{A_1,...,A_n\}\) be a finite collection of non-empty disjoint sets. Then, \[\begin{align} \left(\sqcup_{i=1}^nA_i\right)^c&=\cap_{i=1}^nA_i^c=A_i^c\setminus\left(\sqcup_{j\neq i}A_j\right), \end{align}\] for all \(i\in[n]\). Hence, \[\label{eq:32set-relation} \begin{align} \left(\sqcup_{i=1}^nA_i\right)\times\left[\left(\sqcup_{i=1}^nA_i\right)^c\right]&=\sqcup_{i=1}^n\left[A_i\times\left(A_i^c\setminus\sqcup_{j\neq i}A_j\right)\right]\\ &=\left[\sqcup_{i=1}^n(A_i\times A_i^c)\right]\setminus\left[\sqcup_{i=1}^n\sqcup_{j\neq i}A_i\times A_j\right]. \end{align}\tag{21}\] Let \(A_\gamma=\mathrm{I}_-(\iota_\gamma(\Gamma))\), where \(\gamma\in\partial\Gamma\). Combining Equations 18 and 21 , we obtain \[\begin{align} H_h(\Gamma)=\sum_{\substack{x\in\mathrm{I}_-(\Gamma)\\y\in\mathrm{I}_-(\Gamma)^c}}2J_{xy}=\sum_{\gamma\in\partial\Gamma}\sum_{\substack{x\in\mathrm{I}_-(\iota_\gamma(\Gamma))\\y\in\mathrm{I}_-(\iota_\gamma(\Gamma))^c}}2J_{xy}-\sum_{\substack{\substack{\gamma,\gamma'\in\partial\Gamma\\\gamma\neq\gamma'}}}\sum_{\substack{x\in\mathrm{I}_-(\iota_\gamma(\Gamma))\\y\in\mathrm{I}_-(\iota_{\gamma'}(\Gamma))}}2J_{xy}, \end{align}\]The lemma has thus been proved. ◻

Finally, note that \(\mathscr{C}\) is invariant under integer translations. If \(\gamma_0\in\mathscr{C}\), we define \[\begin{align} \label{eq:32translation-invariance} &[\gamma_0]\mathrel{\vcenter{:}}=\{\gamma_0+k:k\in\mathbb{Z}\}, &\overline{\mathscr{C}}\mathrel{\vcenter{:}}=\{[\gamma_0]:\gamma_0\in\mathscr{C}\}. \end{align}\tag{22}\]

2.4 Crucial estimates↩︎

More important than the geometric objects themselves are the relevant estimates that they satisfy. We list the three that are key to proving the convergence of the cluster expansion.

Hypothesis 1. There exists \(M=M(\alpha)>1\) such that \[H(\Gamma\setminus\gamma)+(7/8)H(\gamma)\leq H(\Gamma)\leq H(\Gamma\setminus\gamma)+H(\gamma),\] for any \(\Gamma\in\mathscr{E}\) and \(\gamma\in\partial\Gamma\). Furthermore, there exists \(\varepsilon>0\) such that \(H(\gamma)\geq \varepsilon\) for all \(\gamma\in\mathscr{C}\). In our case, we may take \(\varepsilon=2\).

Hypothesis 2. There exist \(c_0,c_1>0\), depending on \((M,\alpha)\), and a function \(N:\mathscr{C}\rightarrow\mathbb{N}\), invariant under integer translations, such that \[\label{eq:32cover-size} 1+\log_2\mathrm{diam}(\gamma)\leq N(\gamma)\leq c_0(\log 2)^{-1}H(\gamma)\tag{23}\] and, for all \(n\in\mathbb{N}\), \[\lvert\{\gamma\in\mathscr{C}:(0\in\mathrm{I}_-(\gamma))\wedge(N(\gamma)=n)\}\rvert\leq 2^{c_1n}.\]We note that \[\label{eq:32diameter-hypothesis} \gamma\in\mathscr{C}\Rightarrow \mathrm{diam}(\gamma)\leq e^{c_0 H(\gamma)}.\tag{24}\]

Remark 1. In the one-dimensional case, one can choose the cover size introduced by Fröhlich and Spencer in [9] as the function \(N\).

Hypothesis 3. There exist \(\beta_0=\beta_0(M,\alpha)>0\) and \(c_2=c_2(M,\alpha)>0\) such that \[\label{eq:32pressure-upper-bound} \beta\geq\beta_0\Rightarrow\sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{0\in\mathrm{I}_-(\gamma)}e^{-\beta H(\gamma)}\leq e^{-c_2\beta}.\tag{25}\]

These hypotheses are the crucial ingredients needed for the energy-entropy argument used to prove the existence of a first-order phase transition via contours. If \(\alpha=2\), all were shown to hold in [9]. They were sufficient to obtain a convergent cluster expansion in this case; see [11]. For the region \(\alpha\in(1,2]\), Hypotheses 1, 2 and 3 were shown to hold in [27]. We note that a Peierls’ argument is possible with different types of energy estimates, i.e, modifications to Hypothesis 1, see for instance [14], [19]. Unfortunately, in that case, a restriction on \(\alpha\) appeared, see [20].

Not only are these hypotheses sufficient to prove phase transition, they guaranty the absolute convergence of the cluster expansion. We show that this is the case following a similar approach to that of [17], [18], which deals with a smaller region of polynomial decay \(\alpha\). The same approach was successfully extended to the multidimensional case in [25].

Although Hypotheses 1 and 2 imply Hypothesis 3, it is convenient to have it as a separate hypothesis. As a final observation before describing the appropriate polymer decomposition, we note that \[\begin{align} \label{eq:32hypothesis-3-modified} \sum_{[\gamma]\in\overline{\mathscr{C}}}e^{-\beta H([\gamma])}\leq \sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{0\in\mathrm{I}_-(\gamma)}e^{-\beta H(\gamma)}=\sum_{[\gamma]\in\overline{\mathscr{C}}}\lvert\mathrm{I}_-(\gamma)\rvert e^{-\beta H([\gamma])}. \end{align}\tag{26}\] We will often use this alternative form of Hypothesis 3.

3 Polymers↩︎

In this Section, we describe a polymer construction that is appropriate to the model at hand. We follow very closely the approach developed in [17], [18], [25]: a polymer is a positive collection of contours, as in Definition 3. For this reason, we refer henceforth to positive collections of contours as polymers. Nevertheless, a different type of compatibility must be introduced, for reasons that become clear in the proof of the Lemma 5.

Definition 6. Let \(\Gamma_1,\Gamma_2\in\mathscr{P}\). They are polymerally compatible, which we denote by \(\Gamma_1\sim\Gamma_2\), if one of the following mutually exclusive conditions is satisfied:

  1. \(\Gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\Gamma_2\);

  2. \(\Gamma_2\in\mathscr{I}(\Gamma_1)\), i.e., there exists \(\gamma_1\in\Gamma_1\textrm{ such that }\Gamma_2\in\mathscr{I}(\gamma_1)\);

  3. \(\Gamma_1\in\mathscr{I}(\Gamma_2)\).

We say that \(X\in\mathrm{P}(\mathscr{P})\) is a collection of compatible polymers if \[(\Gamma_1,\Gamma_2\in X)\wedge(\Gamma_1\neq\Gamma_2)\Rightarrow\Gamma_1\sim\Gamma_2.\] We denote the set of all collections of compatible polymers by \(\mathscr{X}\). Similarly to Notation 1, if \(\Lambda\Subset\mathbb{Z}\), we write \[\mathscr{X}_\Lambda\mathrel{\vcenter{:}}=\{X\in\mathscr{X}:\Gamma\in X\Rightarrow\Gamma\in\mathscr{P}_\Lambda\}.\]

Remark 1. We note that if \(\Gamma_1=\{\gamma_1\}\) and \(\Gamma_2=\{\gamma_2\}\), then \(\Gamma_1\sim\Gamma_2\) if and only if \(\gamma_1\sim\gamma_2\).

Remark 1. For each \(\Gamma\in\mathscr{E}\), there exists at least one decomposition \(X\) of \(\Gamma\) into compatible polymers. Indeed, just consider \(X=\{\{\gamma\}:\gamma\in\Gamma\}\).

Convention 3. We establish \(\emptyset\in\mathscr{X}_\Lambda\), for all \(\Lambda\Subset\mathbb{Z}\). In particular, \(\emptyset\in\mathscr{X}\).

Remark 1. Let \(\Gamma_1,\Gamma_2\in\mathscr{P}\). Then \[\label{eq:32characteristic-function-polymer} \mathbbold{1}_{\Gamma_1\sim\Gamma_2}=\left(\mathbbold{1}_{\Gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\Gamma_2}+\sum_{\gamma'_1\in\Gamma_1}\mathbbold{1}_{\Gamma_2\in\mathscr{I}(\gamma'_1)}+\sum_{\gamma'_2\in\Gamma_2}\mathbbold{1}_{\Gamma_1\in\mathscr{I}(\gamma'_2)}\right)\prod_{\substack{\gamma_1\in\Gamma_1\\\gamma_2\in\Gamma_2}}\mathbbold{1}_{\gamma_1\sim\gamma_2}.\qquad{(1)}\] Since all indicative functions inside parentheses are mutually exclusive, we have \[\label{eq:32characteristic-function-not-polymer} \mathbbold{1}_{\Gamma_1\not\sim\Gamma_2}=\mathbbold{1}_{\left\{\substack{\exists\gamma_1\in\Gamma_1,\gamma_2\in\Gamma_2\\\textrm{s.t. }\gamma_1\not\sim\gamma_2}\right\}}+\prod_{\substack{\gamma_1,\gamma'_1\in\Gamma_1\\\gamma_2,\gamma'_2\in\Gamma_2}}\mathbbold{1}_{\gamma_1\sim\gamma_2}\mathbbold{1}_{\Gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr \hss/\hss\cr } } \Gamma_2}\mathbbold{1}_{\Gamma_2\not\in\mathscr{I}(\gamma'_1)}\mathbbold{1}_{\Gamma_1\not\in\mathscr{I}(\gamma'_2)}.\qquad{(2)}\]

Since we may see collections of compatible polymers as partitions of elements of \(\mathscr{E}\), the notion of compatible polymers implies that \[(X\in\mathscr{X})\wedge(Y\preceq X)\Rightarrow Y\in\mathscr{X}.\] That is, if \(X\) is a collection of compatible polymers and \(Y\) is finer than \(X\) (if \(A\in Y\), then there is \(B\in X\) such that \(A\subset B\)), then \(Y\) is also a collection of compatible polymers.

There is also a subjacent natural order to collections of compatible polymers. More precisely, if \(X\in\mathscr{X}\), it is reasonable to define \[\partial X\mathrel{\vcenter{:}}=\{\Gamma_1\in X:[(\Gamma_2\in X)\wedge(\mathrm{I}_-(\Gamma_1)\subset\mathrm{I}_-(\Gamma_2))]\Rightarrow\Gamma_1=\Gamma_2\}.\] In this case, we define recursively \(X_1=\partial X\) and \(X_{k+1}=\partial(X\setminus\sqcup_{i=1}^k X_i)\). It is easy to see that \[X=\bigsqcup_{k=1}^\infty X_k.\] This allows us to see \(X\) as a (graph) forest whose set of vertices is \(X\), and whose set of edges is \[\begin{align} E(X)\mathrel{\vcenter{:}}=\bigsqcup_{k=1}^\infty\left\{\{\Gamma_k,\Gamma_{k+1}\}:(\exists\gamma_k\in\Gamma_k\in X_k:\Gamma_{k+1}\in\mathscr{I}(\gamma_k)\cap X_{k+1})\right\}. \end{align}\]In fact, we can attribute the "color" \(\gamma_k\) to the edge \(\{\Gamma_k,\Gamma_{k+1}\}\), if \(\gamma_k\in\Gamma_k\) and \(\Gamma_{k+1}\in\mathscr{I}(\gamma_k)\).

This very same representation hints at the existence of a single coarsest collection of compatible polymers \(X_\Gamma\) representing the same configuration \(\Gamma\in\mathscr{E}\) as \(X\). If the edges \(\{\Gamma_k,\Gamma_{k+1,1}\}\neq\{\Gamma_k,\Gamma_{k+1,2}\}\) both have the color \(\gamma_k\in\Gamma_k\), it must be true that \(\Gamma_{k+1,1} \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\Gamma_{k+1,2}\), that is, \(\Gamma_{k+1,1}\sqcup\Gamma_{k+1,2}\) is also a polymer. It is easy to see that \[X'=X\sqcup\{\Gamma_{k+1,1}\sqcup\Gamma_{k+1,2}\}\setminus\{\Gamma_{k+1,1},\Gamma_{k+1,2}\}\] is also a collection of compatible polymers. Furthermore, this change does not affect the "grading" described earlier.

One may repeat the procedure described in the previous paragraph as long as there remain two edges with the same color. At this point, since we have not yet made any changes to \(\partial X\), we may not yet have obtained \(X_\Gamma\). The last step is simply to combine all polymers in \(\partial X\) into a single one. Thus, we obtain the coarsest collection of compatible polymers \(X_\Gamma\) representing \(\Gamma\in\mathscr{E}\). Note that the forest representation of \(X_\Gamma\) is a tree.

As a matter of concreteness, let us give two examples.

Figure 3: An example of the procedure whereof one obtains the coarsest decomposition of \Gamma\in\mathscr{E} into compatible polymers.

Let us explain the content of Figure 3. Firstly, for those who are colorblind, let us give an alternative description of the contours. Ordering spin flips so that \(i<j\) implies that \(b_i\) lies to the left of \(b_j\), we have\[\begin{align} &\gamma_1=\{b_1,b_8,b_{11},b_{16}\} &\gamma_2=\{b_2,b_5,b_6,b_7\}\\ &\gamma_3=\{b_3,b_4\} &\gamma_4=\{b_9,b_{10}\}\\ &\gamma_5=\{b_{12},b_{13}\} &\gamma_6=\{b_{14},b_{15}\}\\ &\gamma_7=\{b_{17},b_{18}\}. \end{align}\] Let \(\Gamma=\{\gamma_1,\gamma_2,\gamma_3,\gamma_4,\gamma_5,\gamma_6,\gamma_7\}\) and \(X=\{\Gamma_1,\Gamma_2,\Gamma_3,\Gamma_4,\Gamma_5\}\), where \[\begin{align} &\Gamma_1=\{\gamma_1,\gamma_7\}, &\Gamma_2=\{\gamma_2,\gamma_5\},\\ &\Gamma_3=\{\gamma_3\}, &\Gamma_4=\{\gamma_4\},\\ &\Gamma_5=\{\gamma_6\}. \end{align}\]Then \(X\) is a decomposition of \(\Gamma\) into compatible polymers. We may then fuse \(\Gamma_2=\{\gamma_2,\gamma_5\}\) and \(\Gamma_5=\{\gamma_6\}\), and thus obtain \(X'\), which is a decomposition of \(\Gamma\) into compatible polymers, coarser than \(X\). Finally, we may fuse \(\Gamma_1=\{\gamma_1,\gamma_7\}\) and \(\Gamma_4=\{\gamma_4\}\), and thus obtain \(X_\Gamma\), the coarsest decomposition of \(\Gamma\) into compatible polymers.

Figure 4: Another example of the procedure whereof one obtains the coarsest decomposition of \Gamma\in\mathscr{E} into compatible polymers.

Once again, let us explain the content of Figure 4. For those who are colorblind, let us give an alternative description of the contours. Ordering spin flips so that \(i<j\) implies that \(b_i\) lies to the left of \(b_j\), we have \[\begin{align} &\gamma_1 = \{b_1,b_{12}\},&\gamma_2=\{b_2,b_7\},\\ &\gamma_3=\{b_3,b_4\},&\gamma_4=\{b_5,b_6\},\\ &\gamma_5=\{b_8,b_{11}\},&\gamma_6=\{b_9,b_{10}\}. \end{align}\]Let \(\Gamma=\{\gamma_1,\gamma_2,\gamma_3,\gamma_4,\gamma_5,\gamma_6\}\) and \(X=\{\Gamma_i:i\in[6]\}\), where \(\Gamma_i=\{\gamma_i\}\). We may fuse \(\Gamma_3\) and \(\Gamma_4\); similarly, we may fuse \(\Gamma_2\) and \(\Gamma_5\). Note that \(\Gamma_6\) cannot be glued to \(\Gamma'\mathrel{\vcenter{:}}=\Gamma_3\sqcup\Gamma_4\), for that would break the property of compatibility of the collection of polymers (recall Definition 6). Thus, we obtain the coarsest decomposition of \(\Gamma\).

Lemma 4. Let \(\Gamma\in\mathscr{E}\). Define \[\mathscr{X}_\Gamma\mathrel{\vcenter{:}}=\{X\in\mathbbold{P}(\Gamma):X\preceq X_\Gamma\},\] where \(X_\Gamma\) denotes the coarsest decomposition of \(\Gamma\) into compatible polymers. Then, \[\mathscr{X}=\{\emptyset\}\sqcup\bigsqcup_{\Gamma\in\mathscr{E}}\mathscr{X}_\Gamma.\]Similarly, for any \(\Lambda\Subset\mathbb{Z}\), \[\mathscr{X}_\Lambda=\{\emptyset\}\sqcup\bigsqcup_{\Gamma\in\mathscr{E}_\Lambda}\mathscr{X}_\Gamma.\]

3.1 A polymer gas associated to the long-range Ising model↩︎

In this subsection, we show that, for all \(\Lambda\Subset\mathbb{Z}\), the partition function of the long-range ferromagnetic Ising model at inverse temperature \(\beta>0\) and external field \(h\) is identical to that of a hard-core polymer gas. This decomposition is similar to the one employed in [17], [18], [25].

Let \(\Gamma\in\mathscr{P}\), we define the activity of \(\Gamma\) at temperature \(\beta\) by \[\label{eq:32z-activity} z^+_{\beta,h}(\Gamma)\mathrel{\vcenter{:}}=\frac{Z_{\beta,h}(\Gamma)}{Z^*_{\beta,h}(\Gamma)},\tag{27}\] where (recall Notation 1) \[\begin{align} \label{eq:32bigZ-defin} Z_{\beta,h}(\Gamma)\mathrel{\vcenter{:}}= \sum_{G\in\mathbbold{G}(\Gamma)}\sum_{\substack{\overline{\Gamma}\in\mathscr{E}\\\partial\overline{\Gamma}=\Gamma}}\prod_{\gamma\in\Gamma}e^{-\beta H_{h,\gamma}(\overline{\Gamma})}\prod_{e\in E(G)}\left[e^{\beta\Phi_e(\overline{\Gamma})}-1\right] \end{align}\tag{28}\] and \[\begin{align} \label{eq:32bigZ-star-defin} Z^*_{\beta,h}(\Gamma)\mathrel{\vcenter{:}}= \prod_{\gamma\in\Gamma}Z^*_{\beta,h}(\gamma)\mathrel{\vcenter{:}}=\prod_{\gamma\in\Gamma}\left[\sum_{\overline{\Gamma}\in\mathscr{E}}\mathbbold{1}_{\partial\overline{\Gamma}=\{\gamma\}}e^{-\beta H_{h,\gamma}^*(\overline{\Gamma})}\right]. \end{align}\tag{29}\]

Lemma 5. Let \(\Lambda\Subset\mathbb{Z}\). Then \[\label{eq:32polymer-gas-partition-function} Z^+_{\beta,h,\Lambda}=\sum_{X\in\mathscr{X}_\Lambda}\prod_{\Gamma\in X}z^+_{\beta,h}(\Gamma)\eqqcolon\sum_{X\in\mathrm{P}(\mathscr{P})}z^+_{\beta,h}(X)\mathbbold{1}_{X\in\mathscr{X}_\Lambda}.\tag{30}\]

Proof. Given any finite set \(A\), we write \[\mathrm{P}_2(A)\mathrel{\vcenter{:}}=\{e\in\mathrm{P}(A):\lvert e\rvert=2\}.\]That is, we denote by \(\mathrm{P}_2(A)\) the collection of all subsets of \(A\) that contain exactly two elements. Let \(\{x_e:e\in \mathrm{P}_2(A)\}\subset\mathbb{R}\), then (recall Notation 1)\[\begin{align} \label{eq:32mayer-trick} \prod_{e\in\mathrm{P}_2(A)}(1+x_e)=\sum_{P\in\mathbbold{P}(A)}\prod_{Q\in P}\left[\sum_{G\in\mathbbold{G}(Q)}\prod_{e\in E(G)}x_e\right]. \end{align}\tag{31}\] Let \(\Lambda\Subset\mathbb{Z}\). Then, \[\begin{align} Z^+_{\Lambda,\beta,h}=\sum_{\sigma^+\in\Omega^+_\Lambda}e^{-\beta H_h^+(\sigma^+)}=1+\sum_{\overline{\Gamma}\in\mathscr{E}_\Lambda}e^{-\beta H_h(\overline{\Gamma})}=1+\sum_{\Gamma'_1\in\mathscr{P}_\Lambda}\sum_{\substack{\overline{\Gamma}_1\in\mathscr{E}\\\partial\overline{\Gamma}_1=\Gamma'_1}}e^{-\beta H_h(\overline{\Gamma}_1)}. \end{align}\] By Lemma 3 and Mayer’s trick (Equation 31 ), we have \[e^{-\beta H_h(\overline{\Gamma})}=\prod_{\gamma\in\partial\overline{\Gamma}}e^{-\beta H_{h,\gamma}(\overline{\Gamma})}\sum_{P\in\mathbbold{P}(\partial\overline{\Gamma})}\prod_{Q\in P}\sum_{G\in\mathbbold{G}(Q)}\prod_{e\in E(G)}\left(e^{\beta\Phi_e(\overline{\Gamma})}-1\right),\] for every \(\Gamma\in\mathscr{E}\). Hence,\[\begin{align} \sum_{\substack{\overline{\Gamma}_1\in\mathscr{E}\\\partial\overline{\Gamma}_1=\Gamma'_1}}e^{-\beta H_h(\overline{\Gamma}_1)}=\sum_{\substack{\overline{\Gamma}_1\in\mathscr{E}\\\partial\overline{\Gamma}_1=\Gamma'_1}}\prod_{\gamma_1\in\Gamma'_1}e^{-\beta H_{h,\gamma_1}(\overline{\Gamma}_1)}\sum_{X_1\in\mathbbold{P}(\Gamma'_1)}\prod_{\Gamma_1\in X_1}\sum_{G_1\in\mathbbold{G}(\Gamma_1)}\prod_{e_1\in E(G_1)}\left(e^{\beta\Phi_{e_1}(\overline{\Gamma}_1)}-1\right). \end{align}\]

Let us comment on the chosen notation. We write \(X_1\) in the hope that the reader identifies it with the same notation used earlier, where a graphical description of compatible collections of polymers is provided.

Back to the decomposition algorithm: note that, for each \(X_1\in\mathbbold{P}(\Gamma'_1)\), there is a natural bijection between \(\{\overline{\Gamma}_1\in\mathscr{E}:\partial\overline{\Gamma}_1=\Gamma'_1\}\) and \(\times_{\Gamma_1\in X_1}\{\overline{\Gamma}(\Gamma_1)\in\mathscr{E}:\partial(\overline{\Gamma}(\Gamma_1))=\Gamma_1\}\); the canonical identification being \[(\overline{\Gamma}(\Gamma_1):\Gamma_1\in X_1)\mapsto\overline{\Gamma}_1=\sqcup_{\Gamma_1\in X_1}\overline{\Gamma}(\Gamma_1),\] whose inverse is given by \[\overline{\Gamma}_1\mapsto\left(\iota_{\Gamma_1}(\overline{\Gamma}_1):\Gamma_1\in X_1\right).\] This bijection leaves the values of \(H_{h,\gamma_1}(\overline{\Gamma}_1)\) and \(\Phi_{e_1}(\overline{\Gamma}_1)\) invariant. That is, if \(\gamma_1\in\Gamma_1\), then \[\begin{align} H_{h,\gamma_1}(\overline{\Gamma}_1)=H_{h,\gamma_1}(\overline{\Gamma}(\Gamma_1)). \end{align}\]Similarly, if \(e_1\in\mathrm{P}_2(\Gamma_1)\), then \[\begin{align} \Phi_{e_1}(\overline{\Gamma}_1)=\Phi_{e_1}(\overline{\Gamma}(\Gamma_1)). \end{align}\] In this case, we arrive at \[\begin{align} \label{eq:32iteration-polymer-factorization} \sum_{\substack{\overline{\Gamma}_1\in\mathscr{E}\\\partial\overline{\Gamma}_1=\Gamma'_1}}e^{-\beta H_h(\overline{\Gamma}_1)}=\sum_{X_1\in\mathbbold{P}(\Gamma'_1)}\prod_{\Gamma_1\in X_1}Z_{\beta,h}(\Gamma_1)=\sum_{X_1\in\mathbbold{P}(\Gamma'_1)}\prod_{\Gamma_1\in X_1}z^+_{\beta,h}(\Gamma_1 )Z^*_{\beta,h}(\Gamma_1). \end{align}\tag{32}\] On the other hand, for each \(\gamma_1\in\Gamma_1\), we have \[\begin{align} Z^*_{\beta,h}(\gamma_1)&=\sum_{\substack{\overline{\Gamma}_1\in\mathscr{E}\\\partial\overline{\Gamma}_1=\{\gamma_1\}}}e^{-\beta H^*_{h,\gamma_1}(\overline{\Gamma}_1)}=1+\sum_{\substack{\Gamma'_2\in\mathscr{P}\\\Gamma'_2\in\mathscr{I}(\gamma_1)}}\sum_{\substack{\overline{\Gamma}_2\in\mathscr{E}\\\partial\overline{\Gamma}_2=\Gamma'_2}}e^{-\beta H_{h}(\overline{\Gamma}_2)}. \end{align}\]Since \(\mathscr{I}(\gamma_1)\) might be empty, the sum can be \(0\). This happens if, for example, \(\mathrm{diam}(\gamma_1)\leq 2M\). If not, we proceed with the decomposition by applying Equation 32 . We have thus provided an algorithm that writes the initial partition function in terms of a gas of polymers. ◻

Remark 1. Firstly, note that\[\gamma_{k+1}\in\Gamma'_{k+1}\in\mathscr{I}(\gamma_k)\Rightarrow\mathrm{diam(\gamma_{k+1})}\leq \frac{1}{2M}\mathrm{diam}(\gamma_k)<\frac{1}{2}\mathrm{diam}(\gamma_k).\] This means that this algorithm ends after at most \(\log_2\mathrm{diam}(\Lambda)\) steps. In this case, what does* such an algorithm produce? Certainly, it cannot produce \(z^+_{\beta,h}(X)\), lest \(X\in\mathscr{X}_\Lambda\). That is, it does not produce a collection of polymers that are not compatible. Hence, we obtain an expression of the type \[Z^+_{\Lambda,\beta,h}=1+\sum_{X\in\mathscr{X}_\Lambda}c_\Lambda(X)\prod_{\Gamma\in X}z^+_{\beta,h}(\Gamma).\]Since \(c_\Lambda(\emptyset)=0\), we may naturally identify "\(1\)" outside the summation with \(z_{\beta,h}(\emptyset)\).*

Consider now \(\emptyset\neq X\in\mathscr{X}_\Lambda\). Is \(z_{\beta,h}(X)\) produced by the algorithm described earlier? How many times? The subjacent order to any collection of compatible polymers described after Remark 1 guaranties that should \(z_{\beta,h}(X)\) appear, it does so only once. The fact that it does* appear is a consequence of the term "\(+1\)" in \(Z_{\beta,h}^*(\gamma)\), which gives us an arbitrary stopping step as we run the decomposition algorithm described in the proof.*

3.2 The formal representation of the pressure↩︎

The representation of the model in terms of a gas of polymers is useful as long as the series computing the pressure of the system is absolutely convergent. That is, if the following formal series converges absolutely: \[\label{eq:32pressure-polymer-gas} P_{\beta,h}(\Lambda)\mathrel{\vcenter{:}}=\frac{1}{\lvert\Lambda\rvert}\log Z^+_{\Lambda,\beta,h}=\frac{1}{\lvert\Lambda\rvert}\sum_{n=1}^\infty\frac{1}{n!}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n_\Lambda}\phi_n(\mathbf{\Gamma})\prod_{i=1}^nz^+_{\beta, h}(\Gamma_i),\tag{33}\] where \(\phi_n:\mathscr{P}_\Lambda^n\rightarrow\mathbb{R}\) denotes a Ursell’s function\[\label{eq:32ursell95function} \phi_n(\mathbf{\Gamma})\mathrel{\vcenter{:}}=\sum_{G\in\mathbbold{G}_n}\prod_{\{i,j\}\in E(G)}\left(-\mathbbold{1}_{\Gamma_i\not\sim\Gamma_j}\right)=\sum_{\substack{G\in\mathbbold{G}_n\\E(G)\subset E(\mathbf{\Gamma})}}(-1)^{\lvert E(G)\rvert},\tag{34}\] and \[\label{eq:32edges-vector-Gamma} E(\mathbf{\Gamma})=\{\{i,j\}\in\mathrm{P}_2([n]):\Gamma_i\not\sim\Gamma_j\}.\tag{35}\] We note that \(\phi_n\) is invariant under the action of the symmetric group \(\mathrm{Sym}(n)\), and that \(\phi_n(\mathbf{\Gamma})=0\) if the graph with vertex set \([n]\) and edge set \(E(\mathbf{\Gamma})\) is not connected. One may consult [34], for example, for a derivation of Equation 33 ; see Section 3.1.2.

3.3 Upper bounds to the activity↩︎

Although the decomposition derived in Section 3.1 is exact, it is much more reasonable to prove the absolute convergence of the cluster expansion by studying the convergence of the pressure of another hard-core gas of polymers with a simpler activity \(\bar{z}_{\beta,h}(\cdot)\) that bounds \(z^+_{\beta,h}(\cdot)\) uniformly.

Lemma 6. Let \(\Gamma\in\mathscr{P}\). If \(\beta\geq 8\log 2\), then \[\label{eq:32upper-activity-no-field} z^+_\beta(\Gamma)\leq\prod_{\gamma\in\Gamma}e^{-\frac{1}{2}\beta H(\gamma)}\sum_{T\in\mathbbold{T}(\Gamma)}\prod_{e\in E(T)}\Phi(e)\eqqcolon\bar{z}_\beta(\Gamma).\tag{36}\] Furthermore, if the external field is not trivial, then \[\begin{align} \label{eq:32upper-activity-with-field} z^+_{\beta,h}(\Gamma)\leq e^{\beta E_{\lvert h\rvert}(\Gamma)}\bar{z}_\beta(\Gamma)\eqqcolon\bar{z}_{\beta,h}(\Gamma), \end{align}\tag{37}\] where (recall Equation 12 ) \[\begin{align} \label{eq:32external-field-component-Gamma} E_{\lvert h\rvert}(\Gamma)\mathrel{\vcenter{:}}=\sum_{\gamma\in\Gamma}E_{\lvert h\rvert}(\gamma)=\sum_{\gamma\in\Gamma}\sum_{x\in\mathrm{I}_-(\gamma)}2\lvert h_x\rvert. \end{align}\tag{38}\]

Proof. Let \(\overline{\Gamma}\in\mathscr{E}\) such that \(\partial\overline{\Gamma}=\Gamma\). We begin by establishing a uniform (on \(\overline{\Gamma}\)) bound on \[\label{eq:32g-overline-Gamma-definition} g_\beta(\overline{\Gamma})\mathrel{\vcenter{:}}=\sum_{G\in\mathbbold{G}(\Gamma)}\prod_{e\in E(G)}\left(e^{\beta \Phi_e(\overline{\Gamma})}-1\right).\tag{39}\] Since \(\Phi_e(\overline{\Gamma})\leq \Phi(e)\), for any \(e\in\mathrm{P}_2(\Gamma)\), and every connected graph \(G\in\mathbbold{G}(\Gamma)\) has at least one generating tree \(T\in\mathbbold{T}(\Gamma)\), we have \[g_\beta(\overline{\Gamma})\leq\sum_{T\in\mathbbold{T}(\Gamma)}\prod_{e\in E(T)}\left(e^{\beta\Phi(e)}-1\right)\sum_{\substack{G\in\mathbbold{G}(\Gamma)\\E(T)\subset E(G)}}\prod_{e\in E(G)\setminus E(T)}\left(e^{\beta \Phi(e)}-1\right).\]The non-negativity of \(\Phi(e)\) implies that\[\sum_{\substack{G\in\mathbbold{G}(\Gamma)\\E(T)\subset E(G)}}\prod_{e\in E(G)\setminus E(T)}(e^{\beta\Phi(e)}-1)\leq\sum_{\substack{E\subset\mathrm{P}_2(\Gamma)\\E\cap E(T)=\emptyset}}\prod_{e\in E}(e^{\beta\Phi(e)}-1)=\prod_{e\in\mathrm{P_2}(\Gamma)\setminus E(T)}e^{\beta\Phi(e)}.\]On the other hand, since \(x\geq 0\Rightarrow e^x-1\leq xe^x\), we have \[\prod_{e\in E(T)}\left(e^{\beta\Phi(e)}-1\right)\leq \left(\prod_{e\in E(T)}\beta\Phi(e)\right)\exp\left(\sum_{e\in E(T)}\beta\Phi(e)\right).\]Hence \[g_\beta(\overline{\Gamma})\leq \sum_{T\in\mathbbold{T}(\Gamma)}\prod_{e\in E(T)}\beta\Phi(e)\prod_{e\in P_2(\Gamma)}e^{\beta\Phi(e)}.\]By Lemma 3 and Hypothesis 1, \[\begin{align} \sum_{e\in\mathrm{P}_2(\Gamma)}\Phi(e)=-H(\Gamma)+\sum_{\gamma\in\Gamma}H(\gamma)\leq\frac{1}{8}\sum_{\gamma\in\Gamma}H(\gamma). \end{align}\]Since \(\inf_{\gamma\in\mathscr{C}}H(\gamma)\geq 2\), \(\lvert E(T)\rvert=\lvert \Gamma\rvert-1\) and \(x\geq 0\Rightarrow x\leq e^x\), we have \[\beta\geq1\Rightarrow\prod_{e\in E(T)}\beta\leq\prod_{\gamma\in\Gamma}4(\beta/4)\leq \prod_{\gamma\in\Gamma}4e^{\frac{1}{4}\beta}\leq \prod_{\gamma\in\Gamma}e^{\log4+\frac{1}{8}\beta H(\gamma)}.\] Thus, we obtain \[\label{eq:32g-beta-bound-upper-activity} g_\beta(\overline{\Gamma})\leq \left(\prod_{\gamma\in\Gamma}e^{\log 4+\frac{1}{4}\beta H(\gamma)}\right)\sum_{T\in\mathbbold{T}(\Gamma)}\prod_{e\in E(T)}\Phi(e).\tag{40}\] Let \(\gamma\in\Gamma\). By Hypothesis 1, we have \[\label{eq:32hyp-M-to-choice-upper-activity} -H_\gamma(\overline{\Gamma})\leq -H^*_\gamma(\overline{\Gamma})-(7/8)H(\gamma).\tag{41}\] Combining Equations 40 and 41 with Equation 27 , we obtain\[z^+_\beta(\Gamma)\leq \left(\prod_{\gamma\in\Gamma}e^{\log 4-\frac{5}{8}\beta H(\gamma)}\right)\sum_{T\in\mathbbold{T}(\Gamma)}\prod_{e\in E(T)}\Phi(e).\]Choosing \(\beta\geq 8\log2\), we arrive at Equation 36 .

By Equation 18 , we have \[\label{eq:32field-activity-dealing-with-the-term} H^*_{h,\gamma}(\overline{\Gamma})-H_{h,\gamma}(\overline{\Gamma})=H^*_\gamma(\overline{\Gamma})-H_\gamma(\overline{\Gamma})+\sum_{x\in\mathrm{I}_-(\iota^*_\gamma(\overline{\Gamma}))}2h_x-\sum_{x\in\mathrm{I}_-(\iota_\gamma(\overline{\Gamma}))}2h_x\leq H^*_\gamma(\overline{\Gamma})-H_\gamma(\overline{\Gamma})+\sum_{x\in\mathrm{I}_-(\gamma)}2\lvert h_x\rvert.\tag{42}\] Combining Equation 42 with Equation 36 , we obtain Equation 37 . This concludes the proof of the lemma. ◻

4 Sums over trees↩︎

This Section should be understood as a systematization and extension of the methods employed in [17], [18], [25] to prove the absolute convergence of the cluster expansion. Although it is the traditional method of summing over leaves of trees, in the application to the long-range Ising model, we transform trees of polymers into trees of contours. This implies a change from local, as in for each edge of the tree, (in)compatibility conditions to global compatibility conditions. We can, under this hypothesis and another appropriate hypothesis (see Hypothesis 7), remove vertices of arbitrary degree and prove a rather useful result dealing with trees of contours (see Theorem 4). We deal with trees whose vertices and edges are assigned weights.

Let \(\mathscr{V}\) be a countable set. Consider a family of functions \(v_\beta:\mathscr{V}\rightarrow\mathbb{R}_{>0}\) and another (single) symmetric function \(e:\mathscr{V}^2\rightarrow\mathbb{R}_{\geq 0}\), where \(\beta>1\). We call \(\mathfrak{V}\mathrel{\vcenter{:}}=\{v_\beta:\beta> 1\}\) the vertex function family and \(e(\cdot,\cdot)\) the edge weight function. We define the set of compatible pairs of vertices by \[\mathscr{V}_2\mathrel{\vcenter{:}}=\{(\vartheta_1,\vartheta_2)\in\mathscr{V}^2:e(\vartheta_1,\vartheta_2)\neq 0\}.\]It is a symmetric set. We demand that \(\pi_1\mathscr{V}_2=\mathscr{V}\), where \(\pi_1\) denotes the first projection.

Since the function \(\vartheta\mapsto \log v_{\beta/2}(\vartheta)-\log v_\beta(\vartheta)\) will appear often, let us give it a name. We define the modified vertex function family \(\mathfrak{V}'\mathrel{\vcenter{:}}=\{v'_\beta:\beta>1\}\), where \(v'_\beta:\mathscr{V}\rightarrow\mathbb{R}\) is given by \[\label{eq:32modified-vertex-function} v'_\beta(\vartheta)\mathrel{\vcenter{:}}=\log v_{\beta/2}(\vartheta)-\log v_\beta(\vartheta).\tag{43}\] For \(k\geq 3,\) we also define \[\label{eq:32global-vertex-compatibility} \mathscr{V}_k\mathrel{\vcenter{:}}=\{\mathbf{v}=(\vartheta_1,...,\vartheta_k)\in\mathscr{V}^k:i\neq j\in[k]\Rightarrow(\vartheta_i,\vartheta_j)\in \mathscr{V}_2\}.\tag{44}\] We demand that the (modified) vertex function family \(\mathfrak{V}'\) satisfies the following hypothesis:

Hypothesis 4. The function \(\beta\mapsto v'_\beta(\vartheta)\) is an increasing function for all \(\vartheta\in\mathscr{V}\). Furthermore, \[\lim_{\beta\rightarrow\infty}\inf_{\vartheta\in\mathscr{V}}v'_\beta(\vartheta)=\infty.\]

Remark 1. In general, \(\mathscr{V}_n\subsetneq\mathscr{V}^n\). Nevertheless, it is natural to write \(\mathscr{V}=\mathscr{V}^1=\mathscr{V}_1\). If \(\mathcal{I}\) is a finite set of indices, we may naturally define \(\mathscr{V}^\mathcal{I}\) and \(\mathscr{V}_\mathcal{I}\), accordingly.

Remark 1. In our concrete applications, the modified vertex function \(v'_\beta(\cdot)\) is proportional to \(\beta H(\cdot)\), see Equations 77 and 80 .

Before proceeding, let us fix some useful notation regarding trees.

Notation 1. Let \(V\) be a finite set and \(T\in\mathbbold{T}(V)\). Let \(v\in V\), we define \[\begin{align} &E_T(v)\mathrel{\vcenter{:}}=\{e\in E(T):v\in e\},\\ &N_T(v)\mathrel{\vcenter{:}}=\{w\in V:\{v,w\}\in E(T)\}. \end{align}\]

Definition 7. Let \(T\in\mathbbold{T}(V)\) and \(v\in V\), we define \(\delta_T:V\rightarrow\mathbb{N}\) by \[\label{eq:32degree-definition} \delta_T(v)\mathrel{\vcenter{:}}=\lvert E_T(v)\rvert=\lvert N_T(v)\rvert.\tag{45}\] That is, we denote the function that maps a vertex \(v\) to its degree (in the tree \(T\)) by \(\delta_T\).

Notation 1. Let \(V\) be a finite set. We denote the set of leaves of \(T\in\mathbbold{T}(V)\) by \[L(T)\mathrel{\vcenter{:}}=\{v\in V:\delta_T(v)=1\}.\]Let \(A\subset B\Subset\mathbb{N}\), we define \[\mathbbold{T}(B,A)\mathrel{\vcenter{:}}=\{T\in\mathbbold{T}(B):L(T)= A\}\]If \(B=[n]\), we write \(\mathbbold{T}_n(A)\) instead of \(\mathbbold{T}([n],A)\). We also define \[\overline{\mathbbold{T}}(B,A)\mathrel{\vcenter{:}}= \{T\in\mathbbold{T}(B):L(T)\subseteq A\}.\]If \(B=[n]\), we write \(\overline{\mathbbold{T}}_n(A)\) instead of \(\overline{\mathbbold{T}}([n],A)\). If \(B=[n]\) and \(A=[m]\), we write \(\overline{\mathbbold{T}}_{n,m}\) instead of \(\overline{\mathbbold{T}}_n([m])\). Finally, we define \[\underline{\mathbbold{T}}_{n,m}\mathrel{\vcenter{:}}= \bigsqcup_{A\subset [m]\subset B\subset[n]}\mathbbold{T}(B,A)\]and \[\underline{\mathbbold{T}}_m\mathrel{\vcenter{:}}=\bigcup_{n:n\geq m}\underline{\mathbbold{T}}_{n,m}.\]Note that \(m\leq n_1\leq n_2\Rightarrow\underline{\mathbbold{T}}_{n_1,m}\subset\underline{\mathbbold{T}}_{n_2,m}\) and \(m_1\leq m_2\Rightarrow\underline{\mathbbold{T}}_{m_1}\subset\underline{\mathbbold{T}}_{m_2}\).

Figure 5: A concrete example of Notation 1 where B=\{1,2,4\}.

Let \(A\subset [n]\) and \(T\in\mathbbold{T}_n\). We denote by \(T[A]\) the smallest sub-tree of \(T\) whose set of vertices contains \(A\). If \(A=[m]\), we write \(T[m]\) instead of \(T[[m]]\). Note that \[\label{eq:32leaves-of-smallest-subtree} L(T[A])\subseteq A.\tag{46}\] Let \(V\Subset\mathbb{N}\), \(T\in\mathbbold{T}(V)\) and \(\mathbf{v}\in\mathscr{V}^V\), we define the density of \((T,\mathbf{v})\) by \[w_\beta(T,\mathbf{v})\mathrel{\vcenter{:}}=\prod_{\vartheta\in V}v_\beta(\vartheta)\prod_{\{i,j\}\in E(T)}e(\vartheta_i,\vartheta_j).\]Note that if \(\{i,j\}\in E(T)\) and \(\mathbbold{1}_{\mathscr{V}_2}(\vartheta_i,\vartheta_j)=0\), then \(w_\beta(T,\mathbf{v})=0\).

Let \(\mathbf{w}=(\varOmega_1,...,\varOmega_m)\in\mathscr{V}^m\). Let \(T\in\mathbbold{T}_n\) such that \(m\leq n\), we define the weight of \((T,\mathbf{w})\) by\[\begin{align} \overline{w}_\beta[\mathbf{w}](T)\mathrel{\vcenter{:}}=\sum_{\varphi\in \mathrm{Inj}(m,n)}\frac{1}{n!}\sum_{\mathbf{v}\in\mathscr{V}^n}\prod_{l=1}^m\mathbbold{1}_{\vartheta_{\varphi(l)}=\varOmega_l}w_\beta(T,\mathbf{v}), \end{align}\] where \(\mathrm{Inj}(m,n)\) denotes the set of all injections from \([m]\) to \([n]\).

We have just defined the usual weight seen in a cluster expansion. Let us now define another one, which takes into account the global compatibility condition announced earlier, which holds for the trees of contours, our ultimate interest. We define the modified weight of \((T,\mathbf{w})\) by \[w_\beta[\mathbf{w}](T)\mathrel{\vcenter{:}}=\sum_{\varphi\in \mathrm{Inj}(m,n)}\frac{1}{n!}\sum_{\mathbf{v}\in\mathscr{V}_n}\prod_{l=1}^m\mathbbold{1}_{\vartheta_{\varphi(l)}=\varOmega_l}w_\beta(T,\mathbf{v}).\]Note the change in the sum from \(\mathbf{v}\in\mathscr{V}^n\) to \(\mathbf{v}\in\mathscr{V}_n\).

Figure 6: Visual description of the injection of vertices on a tree.

The term \((n!)^{-1}\) comes from the labeling introduced in the set of trees, and without which there is no hope of computing sums of trees. In which case, fixing the value of specific vertices can be described by injections. For a graphical description, see Figure 6.

Lemma 7. Let \(\mathbf{w}\in\mathscr{V}^m\) and \(m\leq n\). Then \[\label{eq:32bar-w-normal-w-mute-indices} \begin{align} \overline{\mathrm{W}}_{\beta,n}(\mathbf{w})\mathrel{\vcenter{:}}=\sum_{T\in\mathbbold{T}_n}\overline{w}_\beta[\mathbf{w}](T)&=\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}w_\beta(T,\mathbf{v}),\\ \mathrm{W}_{\beta,n}(\mathbf{w})\mathrel{\vcenter{:}}=\sum_{T\in\mathbbold{T}_n}w_\beta[\mathbf{w}](T)&=\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}w_\beta(T,\mathbf{v}). \end{align}\tag{47}\]

Proof. Let \(\varphi\in\mathrm{Inj}(m,n)\) be an injection, we need to define a canonical permutation \(\tilde{\varphi}\mathrel{\vcenter{:}}=\tilde{\varphi}(\varphi)\in \mathrm{Sym}(n)\) such that \((\tilde{\varphi}\circ\varphi)(i)=i\), for all \(i\in [m]\). If \(j\in\mathrm{Im}(\varphi)\), we define \(\tilde{\varphi}(j)=\varphi^{-1}(j)\in[m]\). If \(j\in [n]\setminus\mathrm{Im}(\varphi)\), we simply extend \(\tilde{\varphi}\) such that if \(j,k\in[n]\setminus\mathrm{Im}(\varphi)\) and \(j<k\), then \(\tilde{\varphi}(j)<\tilde{\varphi}(k)\). We have thus constructed a canonical injection \(\mathrm{Inj}(m,n)\ni\varphi\rightarrow\tilde{\varphi}\in \mathrm{Sym}(n)\).

Denoting the usual left-action of \(\mathrm{Sym}(n)\) on \(\mathscr{V}^n\) by \(F_{\cdot}\), the permutation \(\tilde{\varphi}\) induces a bijection \(F_{\tilde{\varphi}}:\mathscr{V}^n\rightarrow\mathscr{V}^n\) given by\[(F_{\tilde{\varphi}}\mathbf{v})_i=\vartheta_{\tilde{\varphi}(i)}.\]Introducing the variable change \(\mathbf{v}'=F_{\tilde{\varphi}^{-1}}\mathbf{v}\), we obtain \[\overline{\mathrm{W}}_{\beta,n}(\mathbf{w})=\sum_{\varphi\in\mathrm{Inj}(m,n)}\frac{1}{n!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{v}'\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}'|_{[m]}=\mathbf{w}}w_\beta (T,F_{\tilde{\varphi}}\mathbf{v}').\]Analogously, consider the left-action \(F^\star_{\cdot}\) of \(\mathrm{Sym}(n)\) on \(\mathbbold{T}_n\) given by \[E(F^\star_{\tilde{\varphi}}T)\mathrel{\vcenter{:}}=\left\{\{\tilde{\varphi}(j),\tilde{\varphi}(k)\}:\{j,k\}\in E(T)\right\}.\] In this case, we have \(w_\beta(T,F_{\tilde{\varphi}}\mathbf{v}')=w_\beta(F^\star_{\tilde{\varphi}}T,\mathbf{v'})\).

Figure 7: We have \tilde{\varphi}(1)=3, \tilde{\varphi}(2)=4,\tilde{\varphi}(3)=2,\tilde{\varphi}(4)=1, \tilde{\varphi}(5)=5,\tilde{\varphi}(6)=6. Note that T\neq F^\star_{\tilde{\varphi}} T.

Since \(F^\star_{\tilde{\varphi}}\) is a bijection for any \(\tilde{\varphi}\in\mathrm{Sym}(n)\), introducing the variable change \(T'=F^\star_{\tilde{\varphi}}T\), we have \[\label{eq:32W-bar-n-definition} \overline{\mathrm{W}}_{\beta,n}(\mathbf{w})=\sum_{\varphi\in\mathrm{Inj}(m,n)}\frac{1}{n!}\sum_{T'\in\mathbbold{T}_n}\sum_{\mathbf{v'\in\mathscr{V}^n}}\mathbbold{1}_{\mathbf{v}'|_{[m]}=\mathbf{w}}w_\beta(T',\mathbf{v}').\tag{48}\] A simple counting argument and the relabeling of mute indices gives us the first part of Equation 47 . The second equality follows from the observation that \(\mathscr{V}_n\) is invariant under permutations of indices. The lemmas has thus been proved. ◻

Given \(\mathbf{w}\in\mathscr{V}^m\), our objective is to estimate\[\begin{align} \overline{\mathrm{W}}_\beta(\mathbf{w})\mathrel{\vcenter{:}}=\sum_{n=m}^\infty\overline{\mathrm{W}}_{\beta,n}(\mathbf{w})\textrm{ and }\mathrm{W}_\beta(\mathbf{w})\mathrel{\vcenter{:}}=\sum_{n=m}^\infty\mathrm{W}_{\beta,n}(\mathbf{w}). \end{align}\]That is, by Lemma 7, the goal is to estimate \[\label{eq:32mute-indices-cluster-formula} \begin{align} \overline{\mathrm{W}}_\beta(\mathbf{w})&=\sum_{n=m}^\infty\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}w_\beta(T,\mathbf{v})=\sum_{n=m}^\infty\overline{\mathrm{W}}_{\beta,n}(\mathbf{w}),\\ \mathrm{W}_\beta(\mathbf{w})&=\sum_{n=m}^\infty\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}w_\beta(T,\mathbf{v})=\sum_{n=m}^\infty\mathrm{W}_\beta(\mathbf{w}). \end{align}\tag{49}\]

Definition 8. We call \(\mathfrak{\overline{W}}\mathrel{\vcenter{:}}=\{\overline{\mathrm{W}}_\beta(\mathbf{w}):\beta>1, m\geq 1, \mathbf{w}\in\mathscr{V}^m\}\) the family of locally compatible sums associated to \((\mathscr{V},\mathfrak{V},e)\).

Similarly, we call \(\mathfrak{W}\mathrel{\vcenter{:}}=\{\mathrm{W}_\beta(\mathbf{w}):\beta>1, m\geq 1, \mathbf{w}\in\mathscr{V}_m\}\) the family of globally compatible sums associated to \((\mathscr{V},\mathfrak{V},e)\).

4.1 One-vertex estimates↩︎

Let us begin by tackling the case \(m=1\). Let \(\mathbf{w}=\varOmega\in\mathscr{V}\). Since we can see any \(T\in\mathbbold{T}_n\) as a tree rooted in \(1\), there is a one-to-one correspondence between \(T\in\mathbbold{T}_n\) and

  1. a number \(t\in [n-1]\), that is, the generation of the rooted tree \(T\);

  2. an ordered partition \((P_1,...,P_t)\) of \([2,n]\);

  3. a family of functions \(\varphi_p:P_p\rightarrow P_{p-1}\), with the convention that \(P_0\equiv\{1\}\). That is, the edges of the tree \(T\) are given by \(\{k_p,\varphi_p(k_p)\}\), where \(p\in [t]\) and \(k_p\in P_p\).

Hence, by separating the sum over \(T\in\mathbbold{T}_1\) and making the change of indices \(n\mapsto n-1\) (this implies the convention \(P_0\equiv\{0\}\) and \(\vartheta_0\equiv\varOmega\)), we obtain \[\begin{align} \frac{\overline{\mathrm{W}}_\beta(\varOmega)}{v_\beta(\varOmega)}=1+\sum_{n=1}^\infty\frac{1}{n!}\sum_{t=1}^n\sum_{P_1,...,P_t}\sum_{\varphi_1,...,\varphi_t}\sum_{\mathbf{v}\in\mathscr{V}^n}\prod_{p=1}^t\prod_{k_p\in P_p}v_\beta(\vartheta_{k_p})e(\vartheta_{k_p},\vartheta_{\varphi_p(k_p)}). \end{align}\] Since every index belonging to \([n]\) is mute, this sum is much easier to compute than it might look at first glance. One must simply count the number of ordered partitions of \([n]\) and note that \[\sum_{\varphi\in B^A}\prod_{a\in A}x_{a,\varphi(a)}=\prod_{a\in A}\sum_{b\in B} x_{a,b}.\] Thus, the series is equal to\[\sum_{n=1}^\infty\sum_{t=1}^n\sum_{\substack{n_1+...+n_t=n\\n_i\geq 1}}\prod_{p=1}^t\frac{1}{n_p!}\prod_{k_p=1}^{n_p}\sum_{l_{p-1}=1}^{n_{p-1}}\sum_{\vartheta_{p,k_p}}v_\beta(\vartheta_{p,k_p})e(\vartheta_{p,k_p},\vartheta_{p-1,l_{p-1}}).\] The first product should be understood as compute everything involving the index \(p=t\), then everything involving the index \(p=t-1\), etc. Exchanging the order of summation and summing over the superfluous index \(n\), we obtain \[\begin{align} \label{eq:32bar-S-definition-leaf-op} \sum_{t=1}^\infty\prod_{p=1}^t\sum_{n_p=1}^\infty\frac{1}{n_p!}\prod_{k_p=1}^{n_p}\sum_{l_{p-1}=1}^{n_{p-1}}\sum_{\vartheta_{p,k_p}}\ell_\beta(\vartheta_{p,k_p},\vartheta_{p-1,l_{p-1}})\eqqcolon\sum^\infty_{t=1}\mathrm{S}_{\beta,t}(\varOmega), \end{align}\tag{50}\] with the convention \(n_0\equiv1\), and \(\vartheta_{0,1}\equiv\varOmega\), and \[\ell_\beta(\vartheta_1,\vartheta_2)\mathrel{\vcenter{:}}= v_\beta(\vartheta_1)e(\vartheta_1,\vartheta_2).\] By Hypothesis 4, if \(\beta>2\), then, for all \((\vartheta_1,\vartheta_2)\in\mathscr{V}^2\), have\[\label{eq:32ell-beta-ell-beta-half} \ell_{\beta}(\vartheta_1,\vartheta_2)\leq\ell_{\beta/2}(\vartheta_1,\vartheta_2).\tag{51}\]

Define the leaf pruning function family \(\mathfrak{L}=\{\mathcal{L}_\beta:\beta>1\}\), where \(\mathcal{L}_\beta:\mathscr{V}\rightarrow\mathbb{R}_{\geq 0}\) is given by\[\label{eq:32leaf-pruning-operator} \mathcal{L}_\beta (\vartheta_1)\mathrel{\vcenter{:}}=\sum_{\vartheta\in\mathscr{V}}\ell_\beta(\vartheta,\vartheta_1).\tag{52}\] Suppose that this function satisfies the following hypothesis.

Hypothesis 5. The family function \(\mathfrak{V}'\) dominates the family function \(\mathfrak{L}\). That is,\[\lim_{\beta\rightarrow\infty}\sup_{\vartheta\in\mathscr{V}}\left\lvert\frac{\mathcal{L}_\beta(\vartheta)}{v'_\beta(\vartheta)}\right\rvert=0,\]where we recall the definition of the modified vertex function present in Equation 43 .

Remark 1. This is very similar to the convergence criterion given by Kotecký-Preiss in [32]. On the other hand, our computation is much closer to that of Fernandéz-Procacci in [31].

Lemma 8. For all \(\eta>0\), there exists \(\beta_0\mathrel{\vcenter{:}}=\beta_0(\eta)>0\) such that if \(\beta\geq\beta_0\) and \(\varOmega\in\mathscr{V}\), then \[\label{eq:32reducing-sums-over-trees} \overline{\mathrm{W}}_\beta(\varOmega)\leq \eta v_{\beta/2}(\varOmega).\tag{53}\]

Remark 1. The natural graphical representation of Lemma 8 is the contraction of trees into a single vertex (with fixed value \(\vartheta=\varOmega\)).

Proof. Recalling Equation 51 , we have \[\begin{align} \mathrm{S}_{\beta,t}(\varOmega)\leq&\prod_{p=1}^{t-1}\sum_{n_p=1}^\infty\frac{1}{n_p!}\prod_{k_p=1}^{n_p}\left(\sum_{l_{p-1}=1}^{n_{p-1}}\sum_{\vartheta_{p,k_p}}\ell_\beta(\vartheta_{p,k_p},\vartheta_{p-1,l_{p-1}})\right)\\&\times \sum_{n_t=1}^\infty\frac{1}{n_t!}\prod_{k_t=1}^{n_t}\sum_{l_{t-1}=1}^{n_{t-1}}\sum_{\vartheta_{t,k_t}}\ell_{\beta/2}(\vartheta_{t,k_t},\vartheta_{t-1,l_{t-1}})\eqqcolon\bar{\mathrm{S}}_{\beta,t}(\varOmega). \end{align}\] We prove the validity of Equation 53 by showing that there exists \(\beta_0>0\) such that if \(\beta\geq \beta_0\), then\[\begin{align} \label{eq:32S-bar-recurrence} \bar{\mathrm{S}}_{\beta,t+1}(\varOmega)\leq 2^{-1}\bar{\mathrm{S}}_{\beta,t}(\varOmega), \end{align}\tag{54}\] for all \(t\in\mathbb{N}\) and all \(\varOmega\in\mathscr{V}\). In this case, \[\begin{align} \frac{\overline{\mathrm{W}}_\beta(\varOmega)}{v_\beta(\varOmega)}\leq 1+\sum_{t=1}^\infty \bar{\mathrm{S}}_{\beta,t}(\varOmega)\leq 1+\sum_{t=1}^\infty2^{-(t-1)}\bar{\mathrm{S}}_{\beta,1}(\varOmega)\leq 2(1+\bar{\mathrm{S}}_{\beta,1}(\varOmega)). \end{align}\]By the definition of \(\bar{\mathrm{S}}_{\beta,1}(\varOmega)\) and Equation 52 , we have \[\begin{align} 1+\bar{\mathrm{S}}_{\beta,1}(\varOmega)&=1+\sum_{n=1}^\infty\frac{1}{n!}\prod_{k=1}^n\sum_{\vartheta_k}\ell_{\beta/2}(\vartheta_k,\varOmega)=\sum_{n=0}^\infty\frac{1}{n!}(\mathcal{L}_{\beta/2}(\varOmega))^n=\exp\mathcal{L}_{\beta/2}(\varOmega). \end{align}\]By Hypothesis 5, if \(\beta\geq\beta_0\) is sufficiently large, then, for every \(\varOmega\in\mathscr{V}\), we have \[\mathcal{L}_{\beta/2}(\varOmega)\leq\frac{1}{2}v'_{\beta/2}(\varOmega)\leq\frac{1}{2} v'_\beta(\varOmega),\]where the second inequality is a consequence of Hypothesis 4. In this case, \(\beta_0>0\) must be chosen so that \[2v_\beta(\varOmega)\exp(2^{-1} v'_\beta(\varOmega))\leq\eta v_{\beta/2}(\varOmega)\iff 4\eta^{-2}\leq\exp(v'_\beta(\varOmega))\]holds for \(\beta\geq\beta_0\) and all \(\varOmega\in\mathscr{V}\). This can certainly be done by Hypothesis 4.

Let us now prove that if \(\beta_0\) is sufficiently large, then \(\bar{\mathrm{S}}_{\beta,t+1}(\varOmega)\leq 2^{-1}\bar{\mathrm{S}}_{\beta,t}(\varOmega)\) for all \(t\geq 1\) and all \(\varOmega\in\mathscr{V}\). By the definition of \(\mathcal{L}_{\beta/2}(\cdot)\), we have \[\sum_{n_{t+1}=1}^{\infty}\frac{1}{n_{t+1}!}\prod_{k_{t+1}=1}^{n_{t+1}}\sum_{l_t=1}^{n_t}\sum_{\vartheta_{t+1,k_{t+1}}}\ell_{\beta/2}(\vartheta_{t+1,k_{t+1}},\vartheta_{t,l_t})= e^{\sum_{l_t=1}^{n_t}\mathcal{L}_{\beta/2}(\vartheta_{t,l_t})}-1.\] Recalling that \(x>0\Rightarrow(e^x-1<xe^x)\wedge(x<e^x)\), we arrive at \[\label{eq:32generation-eating-operation} \sum_{n_{t+1}=1}^{\infty}\frac{1}{n_{t+1}!}\prod_{k_{t+1}=1}^{n_{t+1}}\sum_{l_t=1}^{n_t}\sum_{\vartheta_{t+1,k_{t+1}}}\ell_{\beta/2}(\vartheta_{t+1,k_{t+1}},\vartheta_{t,l_t})\leq \frac{1}{2}\prod_{l_t=1}^{n_t}e^{3\mathcal{L}_{\beta/2}(\vartheta_{t,l_t})}.\tag{55}\] By Hypotheses 4 and 5, there exists \(\beta_0>0\) such that if \(\beta\geq\beta_0\), then \[\begin{align} e^{3\mathcal{L}_{\beta/2}(\vartheta_1)}\ell_\beta(\vartheta_1,\vartheta_2)\leq\ell_{\beta/2}(\vartheta_1,\vartheta_2), \end{align}\]for all \(\vartheta_1,\vartheta_2\in\mathscr{V}\). Combining this inequality with Equation 55 , we obtain Equation 54 . Thus, the lemma has been proved. ◻

4.2 Two-vertex estimates↩︎

We now tackle the case \(m=2\). Recalling Notation 1 and Equation 46 , if \(2\leq m\leq n\), we may define a canonical application \(\Psi_{n,m}:\mathbbold{T}_n\mapsto\underline{\mathbbold{T}}_{n,m}\) by \[\begin{align} \Psi_{n,m}(T)\mathrel{\vcenter{:}}= T[m]. \end{align}\] In particular, if \(m=2\), then this application maps \(\mathbbold{T}_n\) to linear trees with leaves \(\{1,2\}\). The matter is understanding how such an application interacts with the sums over trees.

Figure 8: Note that there might not exist n'\in\mathbb{N} such that \Psi_{n,2}(T)\in\mathbbold{T}_n.

We define the total edge function family \(\mathfrak{\overline{E}}=\{\overline{\mathrm{E}}_\beta:\beta>1\}\), where \(\overline{\mathrm{E}}_\beta:\mathscr{V}^2\rightarrow\mathbb{R}_{\geq0}\) is given by\[\label{eq:32total-edge-function} \overline{\mathrm{E}}_\beta(\varOmega_1,\varOmega_2)\mathrel{\vcenter{:}}=\sum_{n=2}^\infty\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_n=\varOmega_2}\prod_{i=2}^{n-1}\overline{\mathrm{W}}_\beta(\vartheta_i)\prod_{j=1}^{n-1}e(\vartheta_j,\vartheta_{j+1}).\tag{56}\] It is also useful to define \[\overline{\mathrm{E}}_{\beta,n}(\varOmega_1,\varOmega_2)\mathrel{\vcenter{:}}=\sum_{\mathbf{v}\in\mathscr{V}^{n+2}}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_{n+2}=\varOmega_2}\prod_{i=2}^{n+1}\overline{\mathrm{W}}_\beta(\vartheta_i)\prod_{j=1}^{n+1}e(\vartheta_j,\vartheta_{j+1}),\]in which case, we have \[\overline{\mathrm{E}}_{\beta}(\varOmega_1,\varOmega_2)\mathrel{\vcenter{:}}= e(\varOmega_1,\varOmega_2)+\sum_{n=1}^\infty\overline{\mathrm{E}}_{\beta,n}(\varOmega_1,\varOmega_2).\]Note that \(e(\varOmega_1,\varOmega_2)=0\) does not imply that \(\overline{\mathrm{E}}_\beta(\varOmega_1,\varOmega_2)=0\).

We also define the modified total edge function family \(\mathfrak{E}\mathrel{\vcenter{:}}=\{\mathrm{E}_\beta:\beta>1\}\), where \(\mathrm{E}_\beta:\mathscr{V}^2\rightarrow\mathbb{R}_{\geq0}\) is given by\[\mathrm{E}_\beta(\varOmega_1,\varOmega_2)\mathrel{\vcenter{:}}=\sum_{n=2}^\infty\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_n=\varOmega_2}\prod_{i=2}^{n-1}\mathrm{W}_\beta(\vartheta_i)\prod_{j=1}^{n-1}e(\vartheta_j,\vartheta_{j+1}).\]Once again, it is useful to define \[\mathrm{E}_{\beta,n}(\varOmega_1,\varOmega_2)\mathrel{\vcenter{:}}=\sum_{\mathbf{v}\in\mathscr{V}_{n+2}}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_{n+2}=\varOmega_2}\prod_{i=2}^{n+1}\mathrm{W}_\beta(\vartheta_i)\prod_{j=1}^{n+1}e(\vartheta_j,\vartheta_{j+1}),\]in which case, we have \[\mathrm{E}_{\beta}(\varOmega_1,\varOmega_2)\mathrel{\vcenter{:}}= e(\varOmega_1,\varOmega_2)+\sum_{n=1}^\infty\mathrm{E}_{\beta,n}(\varOmega_1,\varOmega_2).\]Note that, in this case, \(e(\varOmega_1,\varOmega_2)=0\) does imply that \(\mathrm{E}_\beta(\varOmega_1,\varOmega_2)=0\).

Lemma 9. Let \(\mathbf{w}=(\varOmega_1,\varOmega_2)\in\mathscr{V}^2\). Then, \[\begin{align} \overline{\mathrm{W}}_\beta(\mathbf{w})&=\overline{\mathrm{W}}_\beta(\varOmega_1)\overline{\mathrm{E}}_\beta(\varOmega_1,\varOmega_2)\overline{\mathrm{W}}_\beta(\varOmega_2),\\ \mathrm{W}_\beta(\mathbf{w})&\leq\mathrm{W}_\beta(\varOmega_1)\mathrm{E}_\beta(\varOmega_1,\varOmega_2)\mathrm{W}_\beta(\varOmega_2). \end{align}\]

Remark 1. The inequality in the case of global compatibility conditions is natural: it is not necessary that gluing globally compatible trees generates another globally compatible tree.

Proof. By Equation 49 , we have \[\overline{\mathrm{W}}_\beta(\mathbf{w})=\sum_{n=2}^\infty\frac{1}{(n-2)!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[2]}=\mathbf{w}}w_\beta(T,\mathbf{v}).\]Introducing the application \(\Psi_{n,2}\), we obtain \[\overline{\mathrm{W}}_\beta(\mathbf{w})=\sum_{n=2}^\infty\frac{1}{(n-2)!}\sum_{T\in\mathbbold{T}_n}\sum_{\underline{T}\in\underline{\mathbbold{T}}_{n,2}}\mathbbold{1}_{\Psi_{n,2}(T)=\underline{T}}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[2]}=\mathbf{w}}w_\beta(T,\mathbf{v}).\]Inverting the order of the sums and recalling that \(\Psi_{n,2}(T)=T[2]\), we have \[\overline{\mathrm{W}}_\beta(\mathbf{w})=\sum_{\underline{T}\in\underline{\mathbbold{T}}_2}\sum_{n=2}^\infty\frac{1}{(n-2)!}\sum_{T\in\mathbbold{T}_n}\mathbbold{1}_{T[2]=\underline{T}}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[2]}=\mathbf{w}}w_\beta(T,\mathbf{v}).\]We note that \(\underline{\mathbbold{T}}_2\) is in one-to-one correspondence with pairs \((\underline{n},\underline{\varphi})\), where \(\underline{n}\geq 2\) is a natural number and \(\underline{\varphi}\) is an injection from \([\underline{n}]\) to \(\mathbb{N}\) such that \(\underline{\varphi}(1)=1\) and \(\underline{\varphi}(\underline{n})=2\). This bijection is given by

  1. \(V(\underline{T})=\underline{\varphi}([\underline{n}])\),

  2. \(E(\underline{T})=\{\{\underline{\varphi}(i),\underline{\varphi}(i+1)\}:i\in [\underline{n}-1]\}\).

We also note that if \(T[2]=\underline{T}\), then \(\textrm{Im}(\underline{\varphi})\subset [n]\). That is, if \(\underline{N}\mathrel{\vcenter{:}}= \max\textrm{Im}(\underline{\varphi})\leq n\).

Figure 9: An example of the correspondences between \underline{T}\in\underline{\mathbbold{T}}_2 and (\underline{n},\underline{\varphi}), and between T\in\Psi_{n,2}^{-1}(\underline{T}) and (\mathbf{P,\mathbf{T}}).

For each \(\underline{T}\in\underline{\mathbbold{T}}_2\), provided \(n\geq\underline{N}\), the set \(\Psi_{n,2}^{-1}(\underline{T})\) is in a one-to-one correspondence with

  1. an ordered partition \(\mathbf{P}=(P_1,...,P_{\underline{n}})\) of \([n]\) such that \(i\in[\underline{n}]\Rightarrow\underline{\varphi}(i)\in P_i\);

  2. a vector of trees \(\mathbf{T}=(T_1,...,T_{\underline{n}})\), such that \(T_i\in\mathbbold{T}(P_i)\).

The correspondence is given by \[E(T)=E(\underline{T})\sqcup\bigsqcup_{i=1}^{\underline{n}}E(T_i).\]In this case, if \(T\in\Psi_{n,2}^{-1}(\underline{T})\), then there exists \((\mathbf{P},\mathbf{T})\) such that \[\label{eq:32weight-usual-cluster-2} w_\beta(T,\mathbf{v})=\prod_{i=1}^{\underline{n}}w_\beta(T_i,\mathbf{v}|_{P_i})\prod_{j=1}^{\underline{n}-1}e(\vartheta_{\underline{\varphi}(j)},\vartheta_{\underline{\varphi}(j+1)}).\tag{57}\] By Equation 57 , we have \[\overline{\mathrm{W}}_\beta(\mathbf{w})=\sum_{\underline{n}=2}^\infty\sum_{\underline{\varphi}}\sum_{n=\underline{N}}^\infty\frac{1}{(n-2)!}\sum_{\mathbf{P}}\sum_{\mathbf{T}}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[2]}=\mathbf{w}}\prod_{i=1}^{\underline{n}}w_\beta(T_i,\mathbf{v}|_{P_i})\prod_{j=1}^{\underline{n}-1}e(\vartheta_{\underline{\varphi}(j)},\vartheta_{\underline{\varphi}(j+1)}),\]and an analogous formula for \(\mathrm{W}_\beta(\mathbf{w})\).

It is useful to separate the sum term indexed by \(\underline{n}=2\). Note that in this case there is a single injection \(\underline{\varphi}\) and \(\underline{N}=2\). Furthermore, the term involving edge functions simplifies to \(e(\varOmega_1,\varOmega_2)\). Since the indices are mute, the matter becomes counting the set of partitions \(\mathbf{P}=(P_1,P_2)\) of \([n]\) such that \(1\in P_1\) and \(2\in P_2\). One readily finds that the term indexed by \(\underline{n}=2\) is equal to \[\overline{\mathrm{W}}_\beta(\varOmega_1)e(\varOmega_1,\varOmega_2)\overline{\mathrm{W}}_\beta(\varOmega_2).\]Since \(\mathscr{V}_{n_1+n_2}\subset\mathscr{V}_{n_1}\times\mathscr{V}_{n_2}\), where equality is not necessary, in the case of \(\mathrm{W}\), the same argument yields an upper bound: the term indexed by \(\underline{n}=2\) is bounded above by \[\mathrm{W}_\beta(\varOmega_1)e(\varOmega_1,\varOmega_2)\mathrm{W}_\beta(\varOmega_2).\]

Henceforth, we suppose that \(\underline{n}\geq 3\). Let us begin by eliminating the sum over \(\underline{\varphi}\) and transform it into a sum over \(\underline{N}\geq\underline{n}\). For each \(\underline{\varphi}\), this will be done by relabeling \(i\in[n]\) (not \(\underline{n}\) or \(\underline{N}\)) so that for \(i\in[\underline{n}]\), \(\underline{\varphi}(i)=i\). This operation will introduce a combinatorial factor which we compute in the following. Given \(\underline{N}\), one must simply count the amount of injections from \([\underline{n}]\) to \(\mathbb{N}\) such that \(\underline{\varphi}(1)=1\), \(\underline{\varphi}(\underline{n})=2\) and \(\max\textrm{Im}(\underline{\varphi})=\underline{N}\). There are \((\underline{n}-2)\) possible indices such that \(\underline{\varphi}(i)=\underline{N}\), in which case one must simply distribute \((\underline{N}-3)\) values over the remaining \((\underline{n}-3)\) indices. That is, one obtains the following combinatorial factor \[(\underline{n}-2)\times(\underline{N}-3)...(\underline{N}-\underline{n}+1)=(\underline{n}-2)\frac{(\underline{N}-3)!}{(\underline{N}-\underline{n})!}.\] The sum over \(\underline{n}\geq3\) becomes \[\sum_{\underline{n}=3}^\infty\sum_{\underline{N}=\underline{n}}^\infty\sum_{n=\underline{N}}^\infty\frac{(\underline{N}-3)!(\underline{n}-2)}{(\underline{N}-\underline{n})!(n-2)!}\sum_{\mathbf{P}}\sum_{\mathbf{T}}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_{\underline{n}}=\varOmega_2}\prod_{i=1}^{\underline{n}}w_\beta(T_i,\mathbf{v}|_{P_i})\prod_{j=1}^{\underline{n}-1}e(\vartheta_{j},\vartheta_{j+1}).\]We note that the ordered partition now satisfies \(i\in P_i\), for \(i\in [\underline{n}]\).

The next step is to relabel the ordered partitions. Consider a vector \(\mathbf{n}=(n_1,...,n_{\underline{n}})\) such that \(n_i\geq 1\) and \(n_1+...+n_{\underline{n}}=n\). Since \(i\in P_i\), there are \[\frac{(n-\underline{n})!}{\prod_{i=1}^{\underline{n}}(n_i-1)!}\]ordered partitions \(\mathbf{P}\) such that \(\lvert P_i\rvert=n_i\). In this case, relabeling them in terms of \(\mathbf{n}\) and recalling Equation 49 , the sum becomes \[\sum_{\underline{n}=3}^\infty\sum_{\underline{N}=\underline{n}}^\infty\sum_{n=\underline{N}}^\infty C(\underline{n},\underline{N},n)\sum_{\mathbf{n}}\sum_{\mathbf{v}\in\mathscr{V}^{\underline{n}}}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_{\underline{n}}=\varOmega_2}\prod_{i=1}^{\underline{n}}\overline{\mathrm{W}}_{\beta,n_i}(\vartheta_i)\prod_{j=1}^{\underline{n}-1}e(\vartheta_j,\vartheta_{j+1}),\]where \[C(\underline{n},\underline{N},n)\mathrel{\vcenter{:}}= \frac{(\underline{n}-2)(\underline{N}-3)!(n-\underline{n})!}{(\underline{N}-\underline{n})!(n-2)!}=\binom{\underline{N}-3}{\underline{n}-3}\binom{n-2}{\underline{n}-2}^{-1}.\]Once again, in the case of \(\mathrm{W}_\beta(\cdot)\), the analogous formula is an upper bound.

The final step consists of showing that, for all \(\mathbf{v}\in\mathscr{V}^{\underline{n}}\), \[\sum_{\underline{N}=\underline{n}}^\infty\sum_{n=\underline{N}}^\infty C(\underline{n},\underline{N},n)\sum_{\mathbf{n}}\prod_{i=1}^{\underline{n}}\overline{\mathrm{W}}_{\beta,n_i}(\vartheta_i)=\prod_{i=1}^{\underline{n}}\sum_{n_i=1}^\infty\overline{\mathrm{W}}_{\beta, n_i}(\vartheta_i)=\prod_{i=1}^{\underline{n}}\overline{\mathrm{W}}_\beta(\vartheta_i),\]where we recall Equation 49 . Noting that \[\sum_{\underline{N}=\underline{n}}^\infty\sum_{n=\underline{N}}^\infty\sum_{\mathbf{n}}=\sum_{n_1=...=n_{\underline{n}}=1}^\infty\sum_{n=\underline{n}}^\infty\sum_{\underline{N}=\underline{n}}^n\mathbbold{1}_{n=n_1+...+n_{\underline{n}}},\]it is sufficient to show that \[\label{eq:32computation-mute-indices} \sum_{n=\underline{n}}^\infty\sum_{\underline{N}=\underline{n}}^n\mathbbold{1}_{n=n_1+...+n_{\underline{n}}}C(\underline{n},\underline{N},n)=1.\tag{58}\] We recall that if \(r,s\geq 0\), then \[\sum_{t=0}^r\binom{s+t}{t}=\binom{r+s+1}{r+1}.\]Hence, \[\sum_{\underline{N}=\underline{n}}^n\binom{\underline{N}-3}{\underline{n}-3}=\sum_{\underline{N}=0}^{n-\underline{n}}\binom{\underline{n}-3+\underline{N}}{\underline{N}}=\binom{n-2}{\underline{n}-2}.\]We have thus shown that Equation 58 holds. The Lemma has thus been proved. ◻

By Lemma 9, it is clear that we must understand the (modified) total edge function. In general, we want to prove a bound of the type \[\label{eq:32total-edge-function-bound} \overline{\mathrm{E}}_\beta(\varOmega_1,\varOmega_2)\leq \eta v_{\beta_1}(\varOmega_1)^{-1}\tilde{e}(\varOmega_1,\varOmega_2)v_{\beta_1}(\varOmega_2)^{-1},\tag{59}\] where \(\beta_1\in\mathbb{R}\) is fixed, \(\tilde{e}:\mathscr{V}^2\rightarrow\mathbb{R}_+\) is a good edge function, and \(\eta>0\) can be made arbitrarily small, provided \(\beta\geq\beta_0=\beta_0(\eta)\) is sufficiently large. Although we have not found a reasonable abstract assumption that proves this bound, we are able to prove it in our case. On the other hand, in the case of the modified total edge function, we can prove, under reasonable hypotheses, that \[\label{eq:32modified-tota-edge-bound} \mathrm{E}_\beta(\varOmega_1,\varOmega_2)\leq \eta v'_{\beta/2}(\varOmega_1)^2e(\varOmega_1,\varOmega_2)v'_{\beta/2}(\varOmega_2)^2.\tag{60}\]

Choosing \(\eta=1\) in Lemma 8, if \(\beta_0\) is sufficiently large, then \(\mathrm{W}_\beta(\vartheta)\leq v_{\beta/2}(\vartheta)\), for all \(\beta\geq\beta_0\) and \(\vartheta\in\mathscr{V}\). In this case, if \(\beta\geq\beta_0\), then \[\label{eq:32E-E-prime-bound} \mathrm{E}_{\beta,n}(\varOmega_1,\varOmega_2)\leq \sum_{\mathbf{v}\in\mathscr{V}_{n+2}}\mathbbold{1}_{\vartheta_1=\varOmega_1}\mathbbold{1}_{\vartheta_{n+2}=\varOmega_2}\prod_{i=2}^{n+1}v_{\beta/2}(\vartheta_i)\prod_{j=1}^{n+1}e(\vartheta_j,\vartheta_{j+1})\eqqcolon \mathrm{E}'_{\beta/2,n}(\varOmega_1,\varOmega_2).\tag{61}\] This motivates us to introduce the two-vertex contracting function family \(\mathfrak{R}^2\mathrel{\vcenter{:}}=\{\mathcal{R}_\beta^2:\beta>1\}\), where \(\mathcal{R}_{\beta}^2:\mathscr{V}^2\rightarrow \mathbb{R}_{\geq 0}\) is given by \[\label{eq:32reducing-operator-definition} \mathcal{R}_{\beta}^2 (\vartheta_1,\vartheta_2)\mathrel{\vcenter{:}}=\sum_{\vartheta\in\mathscr{V}}e(\vartheta_1,\vartheta)v_{\beta}(\vartheta)e(\vartheta,\vartheta_2).\tag{62}\] By Hypothesis 4, the function \(\beta\mapsto \mathcal{R}^2_ \beta(\vartheta_1,\vartheta_2)\) is increasing for any \((\vartheta_1,\vartheta_2)\in\mathscr{V}^2\).

Hypothesis 6. The family function \(\mathfrak{V}'_2\) dominates the family function \(\mathfrak{R}^2\). That is, \[\lim_{\beta\rightarrow\infty}\sup_{(\vartheta_1,\vartheta_2)\in\mathscr{V}_2}\left\lvert\frac{\mathcal{R}_{\beta}^2(\vartheta_1,\vartheta_2)}{e(\vartheta_1,\vartheta_2)[v'_\beta(\vartheta_1)+v'_{\beta}(\vartheta_2)]}\right\rvert=0.\]

Remark 1. We will employ this object only in the case of contours. By Notation 1 and Equation 80 , Hypothesis 6 means that for all \(\eta>0\), there exists \(\beta_0\) such that if \(\beta\geq \beta_0\) and \((\gamma_1,\gamma_2)\in\mathscr{C}_2\), then \[\mathcal{R}^2_\beta(\gamma_1,\gamma_2)\leq \eta[\beta H(\gamma_1)+\beta H(\gamma_2)] \Phi(\gamma_1,\gamma_2).\]

Lemma 10. Let \(\eta\in (0,1)\). There exists \(\beta_0=\beta_0(\eta)\) such that if \(\beta\geq\beta_0\), \(n\in\mathbb{N}\) and \((\varOmega_1,\varOmega_2)\in\mathscr{V}_2\), then \[\mathrm{E}_{\beta,n}(\varOmega_1,\varOmega_2)\leq \eta^nv'_{\beta/2}(\varOmega_1)^2e(\varOmega_1,\varOmega_2)v'_{\beta/2}(\varOmega_2)^2.\]

Proof. First, employing an idea very similar to that in the proof of Proposition 3.7 in [25], we obtain a bound for \(\mathrm{E}'_{\beta,n}(\varOmega_1,\varOmega_2)\). Then, we conclude by Equation 61 . We provide the algorithm operated on \(\mathrm{E}'_{\beta,n}(\varOmega_1,\varOmega_2)\) below. Its crucial feature is the fact that no vertex is visited more than twice.

Figure 10: A graphical description of the algorithm in case n=2. The number to the left is the current step (0 denotes the starting position). A dotted circle represents a weight of 1, a blank circle represents a weight of v_\beta(\cdot) and a gray circle represents a weight of v_{\beta/2}(\cdot). A blank lozenge or square represents a weight of v'_\beta(\cdot). Superposition of figures denotes product of weights.

There is a state \(s\) that takes values in \(\{L,R\}\), which indicates in which direction the pointer will move next: \(L\) is for left and \(R\) is for right. The starting position of the pointer is \(n_0=\lfloor n/2\rfloor+1\) (recall that there are \(n+2\) vertices). If \(n\) is odd, the starting state is \(s_0=L\). If \(n\) is even, the starting state is \(s_0=R\).

After setting up the starting position, we run the following monotone algorithm. More precisely, we run an algorithm that can only increase the value of the estimated sum.

In the first step, we check whether the pointer is over an extremal vertex. If it is, then the algorithm ends. Otherwise, we proceed to the second step.

In the second step, we update the weight of the vertex indicated by the pointer. More precisely, its current weight, which is equal to \(v_\beta(\cdot)v'_\beta(\cdot)^{\xi}\), where \(\xi\in\{0,1,2\}\), is updated to \(2v_{\beta/2}(\cdot)\). This is possible due to Hypothesis 4 and \(x^2\leq 2e^x\), if \(x>0\).

In the third step, the pointer moves in the direction indicated by the state, then it updates (\(L\rightarrow R\) or \(R\rightarrow L\)).

In the fourth step, if \(s=L\) (\(s=R\)), then we remove the vertex to the left (right) of the pointer employing Hypothesis 6. More precisely, if \(\beta\geq \beta_0(\eta)\), then \[\mathcal{R}^2_{\beta/2}(\vartheta_1,\vartheta_2)\leq (\eta/4)[v'_{\beta/2}(\vartheta_1)+v'_{\beta/2}(\vartheta_2)]e(\vartheta_1,\vartheta_2)\leq (\eta/4)[v'_\beta(\vartheta_1)+v'_\beta(\vartheta_2)]e(\vartheta_1,\vartheta_2),\]where the second inequality follows from Hypothesis 4. In this way, a single sum is split into two: in the first, the weight of the vertex indicated by the pointer has been multiplied by \(v'_\beta(\cdot)\); in the second, the vertex to the direction indicated by the state of the pointer has been multiplied by \(v'_\beta(\cdot)\). We move to the first step.

For a detailed graphical description of the algorithm, see Figure 10. By the end, we obtain \[\begin{align} \mathrm{E}'_{\beta,n}(\varOmega_1,\varOmega_2)&\leq 2^n2^n(\eta/4)^n\max\left\{v'_\beta(\varOmega_1)^{\xi_1}v'_\beta(\varOmega_2)^{\xi_2}:\xi_1,\xi_2\in\{0,1,2\}\right\}e(\Omega_1,\Omega_2)\\ &=\eta^nv'_\beta(\varOmega_1)^2e(\varOmega_1,\varOmega_2)v'_\beta(\varOmega_2)^2, \end{align}\]where the equality is a consequence of Hypothesis 4, as long as \(\beta\geq\beta_0\) is sufficiently large. The first term \(2^n\) comes from the number of sums generated as the algorithm runs. The second term \(2^n\) comes from the second step. Recalling Equation 61 , this concludes the proof of the lemma. ◻

Corollary 2. There exists \(\beta_0>0\) such that if \(\beta\geq \beta_0\) and \((\varOmega_1,\varOmega_2)\in\mathscr{V}_2\), then \[\label{eq:32reduction-operator-sums} \mathrm{W}_\beta (\varOmega_{1},\varOmega_2)\leq v_{\beta/4}(\varOmega_{1})e(\varOmega_{1},\varOmega_2)v_{\beta/4}(\varOmega_2).\tag{63}\]

Proof. The proof is a direct application of Lemmas 8 and 10. We have\[\begin{align} \mathrm{W}_\beta(\varOmega_1,\varOmega_2)&\leq v_{\beta/2}(\varOmega_1)e(\varOmega_1,\varOmega_2)v_{\beta/2}(\varOmega_2)+\sum_{n=1}^\infty v_{\beta/2}(\varOmega_1)\mathrm{E}'_{\beta/2,n}(\varOmega_1,\varOmega_2)v_{\beta/2}(\varOmega_2). \end{align}\]By Hypothesis 4, we can bound the first term by \[\begin{align} v_{\beta/2}(\varOmega_1)e(\varOmega_1,\varOmega_2)v_{\beta/2}(\varOmega_2)\leq 2^{-1}v_{\beta/4}(\varOmega_1)e(\varOmega_1,\varOmega_2)v_{\beta/4}(\varOmega_2). \end{align}\]Let us now bound the second term. Choose \(\eta=9^{-1}\) in Lemma 10. Recalling that \[v'_{\beta/2}(\vartheta)^2v_{\beta/2}(\vartheta)\leq 2!v_{\beta/4}(\vartheta)=2v_{\beta/4}(\vartheta),\] we obtain \[\sum_{n=1}^\infty v_{\beta/2}(\varOmega_1)\mathrm{E}'_{\beta/2,n}(\varOmega_1,\varOmega_2)v_{\beta/2}(\varOmega_2)\leq 2^{-1}v_{\beta/4}(\varOmega_1)e(\varOmega_1,\varOmega_2)v_{\beta/4}(\varOmega_2).\]The corollary has thus been proved. ◻

4.3 Many-vertex estimates↩︎

We now tackle the case \(m\geq 3\). It is simpler than it seems, we just need a generalization of Lemma 9 and a generalization of the contraction function to an arbitrary (finite) number of vertices. Let us first prove a simple result on tree sizes that is very useful in accomplishing the first task.

Lemma 11. Let \(m,n\in\mathbb{N}\) such that \(2\leq m\leq n\). Let \(T\in\mathbbold{T}_n\) such that

  1. \(L(T)\subset [m]\);

  2. \(v\in [m+1,n]\Rightarrow \delta_T(v)\geq 3\).

Then, \(n\leq 2m-2\).

Proof. On the one hand, since \(T\) is a tree, we have \(\lvert E(T)\rvert=n-1\). On the other hand, we have \[\begin{align} \lvert E(T)\rvert&=\frac{1}{2}\sum_{v\in [n]}\delta_T(v)=\frac{1}{2}\sum_{v\in [m]}\delta_T(v)+\frac{1}{2}\sum_{v=m+1}^n\delta_T(v)\geq \frac{1}{2}m+\frac{3}{2}(n-m)=\frac{1}{2}(3n-2m). \end{align}\]We thus obtain \(n\leq 2m-2\). ◻

Remark 1. It is not difficult to construct a tree that satisfies the hypotheses of Lemma 11 such that \(n=2m-2\). That is, \(n=2m-2\) is the optimal bound.

Notation 1. Let \(m\leq n\leq 2m-2\), we denote by \(\mathbbold{T}^*_{n,m}\) the set of elements of \(\overline{\mathbbold{T}}_{n,m}\) (recall Notation 1) satisfying the conditions of Lemma 11.

Let \([m]\subset B\) such that \(\lvert B\rvert\leq 2m-2\), we denote by \(\mathbbold{T}^*_m(B)\) the set of elements \(T\in \overline{\mathbbold{T}}(B,[m])\) such that \(b\in B\setminus[m]\Rightarrow \delta_T(b)\geq 3\). Finally, we also write \[\mathbbold{T}^*_m\mathrel{\vcenter{:}}=\bigsqcup_{n=m}^{2m-2}\mathbbold{T}^*_{n,m}\textrm{ and }\overline{\mathbbold{T}}^*_m\mathrel{\vcenter{:}}=\bigsqcup_{\substack{[m]\subset B\\\lvert B\rvert\leq 2m-2}}\mathbbold{T}^*_m(B).\]For a graphical representation of trees in \(\mathbbold{T}^*_m\) for \(m\) small, see Figure 11 below.

Figure 11: All examples of trees in \mathbbold{T}^*_{n,m} for 2\leq m\leq 5. White vertices correspond to those labeled by [m] and black vertices correspond to those labeled by [m+1,n].

Lemma 12. Let \(\mathbf{w}\in(\varOmega_1,...,\varOmega_m)\in\mathscr{V}^m\), where \(m\geq 2\). Then, \[\label{eq:32many-point-arbitrary-sums} \begin{align} \overline{\mathrm{W}}_\beta(\mathbf{w})&=\sum_{n=m}^{2m-2}\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}^*_{n,m}}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}\prod_{i=1}^n\overline{\mathrm{W}}_\beta(\vartheta_i)\prod_{\{j,k\}\in E(T)}\overline{\mathrm{E}}_\beta(\vartheta_j,\vartheta_k),\\ \mathrm{W}_\beta(\mathbf{w})&\leq\sum_{n=m}^{2m-2}\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}^*_{n,m}}\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}\prod_{i=1}^n\mathrm{W}_\beta(\vartheta_i)\prod_{\{j,k\}\in E(T)}\mathrm{E}_\beta(\vartheta_j,\vartheta_k). \end{align}\tag{64}\]

Remark 1. Note that, due to Lemma 11, the combinatorial factor in Equation 64 is immaterial, when it comes to perturbative results.

Proof. Since the combinatorial part of the argument is identical to that of the proof of Lemma 9, we just describe the argument in terms of operations of trees.

In the first step, we contract every tree to ensure that the set of leaves is a subset of \([m]\). Then, we relabel them to ensure that their vertex set is given by \([n]\), for some \(n\geq m\). That is, recalling Notation 1, we obtain \[\label{eq:32m-points-abstract-sum-step-one} \overline{\mathrm{W}}_\beta(\mathbf{w})=\sum_{n=m}^\infty\frac{1}{(n-m)!}\sum_{T\in\overline{\mathbbold{T}}_{n,m}}\sum_{\mathbf{v}\in\mathscr{V}^n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}\prod_{i=1}^n\overline{\mathrm{W}}_{\beta}(\vartheta_i)\prod_{\{j,k\}\in E(T)}e(\vartheta_j,\vartheta_k).\tag{65}\]

In the second step, we contract lines to guaranty that if a vertex is not labeled by an element of \([m]\), then it has degree at least \(3\). Let \(\Psi^*_m:\overline{\mathbbold{T}}_m\rightarrow\mathbbold{T}^*_m\) be the application defined by the composition of the following two maps. First, let \(T\in\overline{\mathbbold{T}}_m\), define \[V(T^*)\mathrel{\vcenter{:}}=[m]\sqcup\{i\in V(T): (i>m)\wedge (\delta_T(i)\geq 3)\}.\]Let \(i,j\in V(T^*)\), we denote the (unordered) path in \(T\) connecting \(i\) and \(j\) by \(\mathfrak{p}_T(i,j)\). In this case, we also define \[E(T^*)=\{\{i,j\}\in\mathrm{P}_2(V(T^*)):V(\mathfrak{p}_T(i,j))\cap V(T^*)=\{i,j\}\}.\]Since \(T^*\mathrel{\vcenter{:}}= (V(T^*),E(T^*))\in\overline{\mathbbold{T}}^*_m\), we have just constructed an application \(\overline{\Psi}^*_m:\overline{\mathbbold{T}}_m\rightarrow\overline{\mathbbold{T}}^*_m\).

Figure 12: A graphic description of the procedure where m=4.

Let \([m]\subset B\) such that \(\lvert B\rvert\leq 2m-2\). We denote by \(\varphi_B:B\rightarrow[\lvert B\rvert]\) the canonical bijection satisfying \[b_1<b_2\in B\Rightarrow\varphi_B(b_1)<\varphi_B(b_2).\]The family \(\{\varphi_B:([m]\subset B\in\mathrm{P}(\mathbb{N}))\wedge(\lvert B\rvert\leq 2m-2)\}\) induces a canonical relabeling application \(\overline{\varphi}_m:\overline{\mathbbold{T}}^*_m\rightarrow\mathbbold{T}^*_m\). Finally, we define \(\Psi^*_m\mathrel{\vcenter{:}}=\overline{\varphi}_m\circ\overline{\Psi}^*_m\). For a visual description, see Figure 12.

Proceeding as in the proof of Lemma 9, we transform Equation 65 into Equation 64 . Once again, we note that in the case of the global compatibility condition, we obtain upper bounds. ◻

Fix \(m\geq 3\). By Hypothesis 4, Lemma [lem:32reducing-sums-over-trees] and Equation 60 (consequence of Lemma 10), if \(\beta\geq\beta_0\) is sufficiently large, we may suppose that \[\begin{align} &\mathrm{W}_\beta(\vartheta)\leq [(8m-12)!]^{-1/m}v_{\beta/2}(\vartheta),\\ &\mathrm{E}_\beta(\vartheta_1,\vartheta_2)\leq v'_{\beta/2}(\vartheta_1)^2e(\vartheta_1,\vartheta_2)v'_{\beta/2}(\vartheta_2)^2, \end{align}\]for all \(\vartheta,\vartheta_1,\vartheta_2\in\mathscr{V}\). By Lemma 12, for all \(\mathbf{w}\in\mathscr{V}_m\), we obtain \[\mathrm{W}_\beta(\mathbf{w})\leq\sum_{n=m}^{2m-2}\frac{1}{(8m-12)!}\sum_{T^*\in\mathbbold{T}^*_{n,m}}\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}\prod_{i=1}^n v'_{\beta/2}(\vartheta_i)^{2\delta_{T^*}(i)}v_{\beta/2}(\vartheta_i)\prod_{\{j,k\}\in E(T^*)}e(\vartheta_j,\vartheta_k).\]On the other hand, since \(x>0\Rightarrow x^n\leq n!e^x\), we have \[v'_{\beta/2}(\vartheta_i)^{2\delta_{T^*}(i)}v_{\beta/2}(\vartheta_i)\leq (2\delta_{T^*}(i))!v_{\beta/4}(\vartheta_i),\]and \[\sum_{i=1}^n2\delta_{T^*}(i)=4(n-1)\leq 8m-12\Rightarrow\prod_{i=1}^n(2\delta_{T^*}(i))!\leq(8m-12)!.\]Therefore, \[\label{eq:32W39-definition} \mathrm{W}_\beta(\mathbf{w})\leq\sum_{n=m}^{2m-2}\sum_{T^*\in\mathbbold{T}^*_{n,m}}\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}\prod_{i=1}^nv_{\beta/4}(\vartheta_i)\prod_{\{j,k\}\in E(T^*)}e(\vartheta_j,\vartheta_k)\eqqcolon\sum_{n=m}^{2m-2}\mathrm{W}^*_{\beta/4,n}(\mathbf{w}).\tag{66}\] This motivates us to define the following generalization of the contracting function introduced in Equation 62 . Let \(k\geq 1\), we define the \(k\)-th contracting function family \(\mathfrak{R}^k=\{\mathcal{R}^k_\beta:\beta>1\}\), where \(\mathcal{R}^k_{\beta}:\mathscr{V}^k\rightarrow \mathbb{R}_{\geq 0}\) is given by \[\label{eq:32kth-contracting-function-definition} \mathcal{R}^k_{\beta}(\vartheta_1,...,\vartheta_k)\mathrel{\vcenter{:}}=\sum_{\vartheta\in\mathscr{V}}v_\beta(\vartheta)\prod_{i=1}^ke(\vartheta,\vartheta_i).\tag{67}\] Note that \(\mathfrak{L}\equiv\mathfrak{R}^1\). We suppose that it satisfies the following generalization of Hypothesis 6.

Hypothesis 7. For all \(k\geq1\), The family function \(\mathfrak{V}'_k\) dominates the family function \(\mathfrak{R}^k\). That is, \[\label{eq:32kth-contracting-function-hypothesis} \lim_{\beta\rightarrow\infty}\sup_{\mathbf{v}\in\mathscr{V}_k}\left\lvert\frac{\mathcal{R}_{\beta}^k(\mathbf{v})}{\sum_{\varphi\in\mathrm{Sym}(k)}v'_\beta(\vartheta_{\varphi(k)})\prod_{i=1}^{k-1}e(\vartheta_{\varphi(i)},\vartheta_{\varphi(i+1)})}\right\rvert=0,\tag{68}\] where we recall that \(S_k\) denotes the \(k\)-th symmetric group.

We note that Remark 1 extends to \(k\)-th contracting functions, where \(k\geq3\). Nevertheless, the graphical interpretation is more subtle in this case. If \(m\leq n\) and \(T\in\mathbbold{T}_n\), then we write \[\mathrm{W}_\beta(\mathbf{w},T)\mathrel{\vcenter{:}}=\sum_{\mathbf{v}\in\mathscr{V}_n}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}w_\beta(T,\mathbf{v}).\]Consider the case \(m<n\) and \(\delta_T(n)\geq3\). Recalling Notation 1, we write \[\delta\equiv\delta_T(n),\textrm{ }E\equiv E_T(n)\textrm{ and }N_T(n)=\{l_1,...,l_\delta\},\] where \(l_1<...<l_\delta\). Thus, \[\mathrm{W}_\beta(\mathbf{w},T)=\sum_{\mathbf{v}\in\mathscr{V}_{n-1}}\mathbbold{1}_{\mathbf{v}|_{[m]}=\mathbf{w}}\prod_{i=1}^{n-1}v_\beta(\mathbf{\vartheta}_i)\prod_{\{j,k\}\in E(T)\setminus E}e(\vartheta_j,\vartheta_k)\times \mathcal{R}^\delta_\beta(\vartheta_{l_1},....,\vartheta_{l_\delta}).\]We begin by noting that removing the vertex \(n\) induces an ordered partition \(\mathbf{P}=(Q_1,...,Q_\delta)\) of \([n-1]\) such that if \(p\in[\delta]\), then \(l_p\in Q_p\) and \(T|_{Q_p}\in\mathbbold{T}(Q_p)\), where \(E(T|_{Q_p})\mathrel{\vcenter{:}}= E(T)\cap\mathbbold{P}(Q_p)\). Hence, \[\prod_{i=1}^{n-1}v_\beta(\vartheta_i)\prod_{\{j,k\}\in E(T)\setminus E}e(\vartheta_j,\vartheta_k)=\prod_{p=1}^\delta w_\beta(\mathbf{v}|_{Q_p},T|_{Q_p}).\]If \(\eta_\delta>0\), then, by Hypothesis 7, there exists \(\beta_0=\beta_0(\delta,\eta_\delta)\) such that \(\beta\geq\beta_0\) implies \[\mathcal{R}^\delta_\beta(\vartheta_{i_1},...,\vartheta_{i_\delta})\leq \eta_\delta\sum_{\varphi\in\mathrm{Sym}(\delta)}v'_\beta(\vartheta_{l_{\varphi(\delta)}})\prod_{p=1}^{\delta-1}e(\vartheta_{l_{\varphi(p)}},\vartheta_{l_{\varphi(p+1)}}).\] Combining this fact with Hypothesis 4 and \(v_\beta(\cdot)v'_\beta(\cdot)\leq v_{\beta/2}(\cdot)\), we obtain \[\mathrm{W}_\beta(\mathbf{w},T)\leq \eta_\delta\sum_{\varphi\in\mathrm{Sym}(\delta)}\sum_{\mathbf{v}\in\mathscr{V}_{n-1}}\prod_{p=1}^\delta w_{\beta/2}(\mathbf{v}|_{Q_p},T|_{Q_p})\prod_{q=1}^{\delta-1}e(\vartheta_{l_{\varphi(q)}},\vartheta_{l_{\varphi(q+1)}}),\] which is analogous to Equation 57 . That is, for every \(\varphi\in\mathrm{Sym}(\delta)\), there exists a (unique) tree \(T_\varphi\in\mathbbold{T}_{n-1}\) such that if \(\beta>0\) and \(\mathbf{v}\in\mathscr{V}_{n-1}\), then\[w_\beta(\mathbf{w},T_\varphi)=w_{\beta}(\mathbf{v}|_{Q_p},T|_{Q_p})\prod_{q=1}^{\delta-1}e(\vartheta_{l_{\varphi(q)}},\vartheta_{l_{\varphi(q+1)}}).\]and \(E(T_\varphi)=(E(T)\setminus E)\cup \{\{l_{\varphi(p)},l_{\varphi(p+1)}\}:p\in[\delta-1]\}\).

If \(\varphi_1,\varphi_2\in\mathrm{Sym}(\delta)\), then \(T_{\varphi_1}=T_{\varphi_2}\) if, and only if, \[\{\{{\varphi_1(p)},{\varphi_1(p+1)}\}:p\in[\delta-1]\}=\{\{{\varphi_2(p)},{\varphi_2(p+1)}\}:p\in[\delta-1]\}.\]It is not hard to see that either \(\varphi_1=\varphi_2\) or \(\varphi_1=\varphi_2\circ \varphi_{\delta,r}\), where \(\varphi_{\delta,r}\in\mathrm{Sym}(\delta)\) denotes the reversal permutation; that is, \(\varphi_{\delta,r}(i)=\delta-i\) for \(i\in[\delta]\). There is a final, and crucial, observation to be made: if \(T\in\mathbbold{T}_n\) and \(\varphi\in\mathrm{Sym}(\delta)\), then \(\delta_T(i)\leq \delta_{T_\varphi}(i)\) for all \(i\in[n-1]\). That is, if \(m<n\leq 2m-2\) and \(T^*\in\mathbbold{T}^*_{n,m}\), then \(T^*_\varphi\in\mathbbold{T}^*_{n-1,m}\) for all \(\varphi\in\mathrm{Sym}(\delta)\). We conclude that if \(m<n\leq 2m-2\), \(T^*\in\mathbbold{T}^*_{n,m}\) and \(\beta\geq\beta_0(\delta,\eta)\), then\[\label{eq:32W-beta-T-star} \mathrm{W}_\beta(\mathbf{w},T^*)\leq \eta_\delta\sum_{\varphi\in\mathrm{Sym}(\delta)}\mathrm{W}_{\beta/2}(\mathbf{w},T^*_\varphi)\leq 2\eta_\delta\sum_{\tilde{T}^*\in\mathbbold{T}^*_{n-1,m}}\mathrm{W}_{\beta/2}(\mathbf{w},\tilde{T}^*)=2\eta_\delta\mathrm{W}_{\beta/2,n-1}^*(\mathbf{w}).\tag{69}\] Figure 13 provides us with a useful visual description of the operation on trees associated with the \(k\)-th contracting function, and how it relates to Hypothesis 7. Although these graph operations are very complex (resolving \(\mathcal{R}_{\beta}^k\) generates \(k!/2\) distinct trees), this does not pose us any problems. That is because of Lemma 11.

Figure 13: Graphs arising after removing the black vertex. Squares indicate the vertices carrying the extra v'_\beta(\cdot) term. Note that this graph operation preserves the set of leaves. Note that the degree of every remaining vertex does not decrease.

Theorem 4. Let \(m\geq 1\). There exists \(\beta_0\mathrel{\vcenter{:}}=\beta_0(m)>0\) such that if \(\beta\geq\beta_0\) and \(\mathbf{w}\in\mathscr{V}_m\), then \[\label{eq:32m-point-bound} \mathrm{W}_\beta(\mathbf{w})\leq \sum_{T\in\mathbbold{T}_m}\prod_{i=1}^mv_{2^{-m}\beta}(\varOmega_i)\prod_{\{j,k\}\in E(T)}e(\varOmega_j,\varOmega_k).\tag{70}\]

Proof. The cases \(m=1\) and \(m=2\) have already been dealt with, and hence we assume that \(m\geq 3\). The starting point is Equation 66 . Our strategy is to show that if \(m\leq n\leq n+1\leq 2m-2\), then \[\label{eq:32iterative-step-m-point-bound} \mathrm{W}^*_{\beta,n+1}(\mathbf{w})\leq \frac{1}{2}\mathrm{W}^*_{\beta/2,n}(\mathbf{w}).\tag{71}\] This is indeed sufficient: assume that Equation 71 holds. An easy induction argument proves that if \(0\leq k\leq m-2\), then \(\mathrm{W}^*_{\beta,m+k}(\mathbf{w})\leq 2^{-k}\mathrm{W}^*_{2^{-k}\beta,m}(\mathbf{w})\). Therefore, \[\label{eq:32consequence-of-iterative-step-m-point-bound} \mathrm{W}_\beta(\mathbf{w})\leq\sum_{n=m}^{2m-2}\mathrm{W}^*_{\beta/4,n}(\mathbf{w})= \sum_{k=0}^{m-2}\mathrm{W}^*_{\beta/4,m+k}(\mathbf{w})\leq\sum_{k=0}^{m-2}\frac{1}{2^k}\mathrm{W}^*_{2^{-2-k}\beta,m}(\mathbf{w}).\tag{72}\] Since \(\mathbbold{T}^*_{m,m}=\mathbbold{T}_m\), we have \[\mathrm{W}^*_{2^{-2-k}\beta,m}(\mathbf{w})=\sum_{T\in\mathbbold{T}_m}\prod_{i=1}^mv_{2^{-2-k}\beta}(\varOmega_i)\prod_{\{j,k\}\in E(T)}e(\varOmega_j,\varOmega_k).\]Thus, if \(0\leq k<m-2\), then, by Hypothesis 4, we may suppose that \[\label{eq:32W39-same-index-decreasing-beta-bound} \mathrm{W}^*_{2^{-2-k}\beta,m}(\mathbf{w})\leq\frac{1}{2}\sum_{T\in\mathbbold{T}_m}\prod_{i=1}^mv_{2^{-m}\beta}(\varOmega_i)\prod_{\{j,k\}\in E(T)}e(\varOmega_j,\varOmega_k)=\frac{1}{2}\mathrm{W}^*_{2^{-m}\beta,m}(\mathbf{w}).\tag{73}\] Inserting Equation 73 into Equation 72 , we obtain \[\mathrm{W}_\beta(\mathbf{w})\leq\left(\sum_{k=0}^{m-3}\frac{1}{2^{k+1}}+\frac{1}{2^{m-2}}\right)\mathrm{W}^*_{2^{-m}\beta,m}(\mathbf{w})=\sum_{T\in\mathbbold{T}_m}\prod_{i=1}^mv_{2^{-m}\beta}(\varOmega_i)\prod_{\{j,k\}\in E(T)}e(\varOmega_j,\varOmega_k).\]To finish the proof, let us show that Equation 71 holds.

By Hypothesis 7, there exists \(\beta_0\) such that if \(\beta\geq\beta_0\), then \[\mathcal{R}^k_\beta(\mathbf{v})\leq (4\lvert\mathbbold{T}_{2m-2}\rvert)^{-1}\sum_{\varphi\in\mathrm{Sym}(k)}v'_\beta(\vartheta_{\varphi(k)})\prod_{i=1}^{k-1}e(\vartheta_{\varphi(i)},\vartheta_{\varphi(i+1)}),\]for \(3\leq k\leq 2m-3\) and all \(\mathbf{v}\in\mathscr{V}_k\). Hence, recalling Equation 69 , we obtain \[\mathrm{W}^*_{\beta,n+1}(\mathbf{w})=\sum_{T^*\in\mathbbold{T}^*_{n+1,m}}\mathrm{W}_\beta(\mathbf{w},T^*)\leq \sum_{T^*\in\mathbbold{T}^*_{n+1,m}}(2\lvert\mathbbold{T}_{2m-2}\rvert)^{-1}\mathrm{W}^*_{\beta/2,n}(\mathbf{w})\leq\frac{1}{2}\mathrm{W}^*_{\beta/2,n}(\mathbf{w}),\]where the final inequality holds since \(T^*\in\mathbbold{T}^*_{n+1,m}\) implies \(\delta_{T^*}(n+1)\in [3,2m-3]\) and \(m\leq n\leq 2m-2\) implies \(\lvert\mathbbold{T}^*_{n,m}\rvert\leq\lvert\mathbbold{T}_{2m-2}\rvert\). The theorem has thus been proved. ◻

5 Convergence of the polymer gas expansion↩︎

In this Section, we prove that the formal series in Equation 33 converges absolutely provided \(\beta\) is sufficiently large. Employing a graph-tree bound, we transform the sum over graphs into a sum over trees, which allows us to prove absolute convergence by applying twice the framework developed in Section 4.

5.1 Implementing the program↩︎

First, we recall the following graph-tree bound (see [39], or Proposition 5 in [31], or the generalist approach in [40]): \[\label{eq:32graph-tree-bound} \lvert\phi_n(\mathbf{\Gamma})\rvert\leq \sum_{T\in\mathbbold{T}_n}\mathbbold{1}_{E(T)\subset E(\mathbf{\Gamma})}.\tag{74}\] Therefore, recalling Equation 33 and Lemma 6, to establish the absolute convergence of the pressure, it is sufficient to show that \[\label{eq:32upper-pressure} \overline{P}_{\beta,h}(\Lambda)\mathrel{\vcenter{:}}=\frac{1}{\lvert \Lambda\rvert}\sum_{n=1}^\infty\frac{1}{n!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n_\Lambda}\prod_{i=1}^n\bar{z}_{\beta,h}(\Gamma_i)\prod_{\{j,k\}\in E(T)}\mathbbold{1}_{\Gamma_j\not\sim\Gamma_k}\tag{75}\] is uniformly bounded with respect to \(\Lambda\in\mathrm{P}(\mathbb{Z})\). Since we are interested in applying the cluster expansion to compute the correlations of the systems, we limit ourselves to dealing with the case where \(h\) is non-trivial in a finite (fixed) set of lattice sites. For this reason, we may suppose (see the proof of Theorem 1) that \(h\equiv 0\), which we shall do for the rest of this section.

We are interested in two types of sums of trees: trees of polymers and trees of contours. Recalling Definition 8, in the case of polymers, we are interested in the family of locally compatible sums \(\mathfrak{\overline{W}}\) associated to the triple \[(\mathscr{P},\{\bar{z}_\beta(\cdot)\},\mathbbold{1}_{\cdot\not\sim\cdot}).\] In the case of contours, we are interested in the family of globally compatible sums \(\mathfrak{W}\) associated to the triple (recall Definition 3) \[(\mathscr{C},\{e^{-\beta H(\cdot)}\},\Phi(\cdot,\cdot)\mathbbold{1}_{\cdot \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\cdot}).\] Since each sum is of a different type, we may establish the following convention:

Convention 4. Whenever we write \(\overline{\mathrm{W}}_\beta(\cdot)\), we mean a (locally compatible) sum over trees of polymers, and whenever we write \(\mathrm{W}_\beta(\cdot)\), we mean a (globally compatible) sum over trees of contours.

One should readily recognize that \[\label{eq:32pressure-tree-function} \overline{P}_\beta(\Lambda)\leq\frac{1}{\lvert\Lambda\rvert}\sum_{\Gamma\in\mathscr{P}_\Lambda}\overline{\mathrm{W}}_{\beta}(\Gamma).\tag{76}\] In order to apply Lemma 8, which shall be sufficient to prove the absolute convergence of the cluster expansion, we must show that Hypothesis 4 holds for \(\{\bar{z}_\beta(\cdot)\}\) and the appropriate leaf pruning function satisfies Hypothesis 5.

Convention 5. In order to distinguish between the leaf pruning function for polymers and for contours, we will denote the first by \(\overline{\mathcal{L}}_\beta\) and the second by \(\mathcal{L}_\beta\).

The first condition is very easy to check. For any \(\Gamma\in\mathscr{P}\),\[\begin{align} \label{eq:32v-prime-function-polymers} \beta\mapsto\log\frac{\bar{z}_{\beta/2}(\Gamma)}{\bar{z}_{\beta}(\Gamma)}=\sum_{\gamma\in\Gamma}\frac{1}{4}\beta H(\gamma)\textrm{ is increasing and}\inf_{\Gamma\in\mathscr{P}}\log\frac{\bar{z}_{\beta/2}(\Gamma)}{\bar{z}_{\beta}(\Gamma)}\geq\frac{1}{2}\beta\xrightarrow{\beta\rightarrow\infty}\infty. \end{align}\tag{77}\] The second condition is the content of Lemma 15. It is reasonable that it should be proved later, for its proof rests on the application of the machinery of Section 4 to trees of contours: the content of the next subsection.

5.2 Sums of trees of contours↩︎

Assume that \(\overline{\mathcal{L}}_\beta\) satisfies Hypothesis 5. In this case, choosing \(\eta=1\) in Lemma 8, for \(\beta\geq \beta_0\) sufficiently large, we have \[\begin{align} \label{eq:32pressure-intermediary-bound} \overline{P}_\beta(\Lambda)\leq \frac{1}{\lvert\Lambda\rvert}\sum_{\Gamma\in\mathscr{P}_\Lambda}\bar{z}_{\beta/2}(\Gamma). \end{align}\tag{78}\]

Remark 1. Although Equation 78 is only valid if \(\overline{\mathcal{L}}_\beta\) satisfies Hypothesis [hyp:32leaf-pruning-operator], every estimate involving sums over trees of contours shown in this section is true regardless. In fact, as mentioned earlier, the proof of Lemma 15 is based on such estimates.

Notation 1. We write \[\mathscr{C}_k=\{(\gamma_1,...,\gamma_k)\in\mathscr{C}^k:i\neq j\Rightarrow\gamma_i \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_j\}.\]

Recalling Equation 36 , if \(m\in\mathbb{N}\) and \(\mathbf{w}=(\gamma_1,...,\gamma_n)\in\mathscr{C}_m\), then\[\label{eq:32bar-zeta-compared-to-W} \sum_{\Gamma\in\mathscr{P}}\bar{z}_\beta(\Gamma)\prod_{i=1}^m\mathbbold{1}_{\gamma_i\in\Gamma}= \mathrm{W}_{\beta/2}(\mathbf{w}),\tag{79}\] which justifies the choice of the triple \((\mathscr{C},\{e^{-\beta H(\cdot)}\},\Phi(\cdot,\cdot)\mathbbold{1}_{\cdot \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\cdot})\). Similarly to the case of polymers, we must show that \(\{e^{-\beta H(\cdot)}\}\) satisfies Hypothesis 4. If \(\gamma\in\mathscr{C}\), then \[\label{eq:32modified-vertex-function-contours} \log\frac{e^{-\frac{1}{2}\beta H(\gamma)}}{e^{-\beta H(\gamma)}}=\frac{1}{2}\beta H(\gamma)\textrm{ is increasing in \beta and}\inf_{\gamma\in\mathscr{C}}\log\frac{e^{-\frac{1}{2}\beta H(\gamma)}}{e^{-\beta H(\gamma)}}\geq\beta\xrightarrow{\beta\rightarrow\infty}\infty.\tag{80}\] In order to apply Lemma 8, we must now show that the appropriate leaf pruning function \(\mathcal{L}_\beta:\mathscr{C}\rightarrow \mathbb{R}\) satisfies Hypothesis 5. This is the content of the following lemma:

Lemma 13. There exists \(\beta_0>0\) such that, if \(\beta\geq\beta_0\) and \(\gamma_1\in\mathscr{C}\), then \[\mathcal{L}_\beta(\gamma_1)\leq 2e^{-c_2\beta/2} H(\gamma_1).\]

Proof. By definition, we have \[\mathcal{L}_\beta(\gamma_1)=\sum_{\gamma\in\mathscr{C}}e^{-\beta H(\gamma)}\Phi(\gamma,\gamma_1)\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}=\sum_{[\gamma]\in\bar{\mathscr{C}}}e^{-\beta H([\gamma])}\sum_{\gamma\in[\gamma]}\Phi(\gamma,\gamma_1)\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1},\]where we introduce the equivalence classes of contours under integer translations.

Let us estimate the sum over an equivalence class \([\gamma]\). We have \[\sum_{\substack{\gamma\in[\gamma]}}\Phi(\gamma,\gamma_1)\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}=\sum_{\substack{x\in\mathrm{I}_-(\gamma_1)\\y\in\mathrm{I}_-(\gamma_1)^c}}4J_{xy}\sum_{\gamma\in[\gamma]}\mathbbold{1}_{y\in\mathrm{I}_-(\gamma)}\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}.\] Since \(\mathbbold{1}_{\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma}\leq 1\) and at most \(\lvert\mathrm{I}_-(\gamma)\rvert\leq \mathrm{diam}(\gamma)\) elements \(\gamma\in[\gamma]\) satisfy \(y\in\mathrm{I}_-(\gamma)\), we obtain \[\sum_{\substack{\gamma\in[\gamma]}}\Phi(\gamma,\gamma_1)\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}\leq 2\mathrm{diam}(\gamma)H(\gamma_1).\] Finally, by Hypotheses 2 and 3, we obtain \[\begin{align} \mathcal{L}_\beta(\gamma)\leq 2\sum_{[\gamma_1]\in\overline{\mathscr{C}}}e^{-(\beta-c_0)H([\gamma_1])} H(\gamma)\leq 2e^{-c_2\beta/2} H(\gamma), \end{align}\] where the last inequality follows by choosing \(\beta_0\geq 2c_0\) (recall Hypothesis 3). The lemma has thus been proved. ◻

Corollary 3. There exists \(\beta_0>0\) such that, if \(\beta\geq\beta_0\) and \(\gamma_1\in\mathscr{C}\), then \[\mathrm{W}_\beta(\gamma_1)\leq e^{-\frac{1}{2}\beta H(\gamma_1)}.\]

Proof. This is a direct consequence of Lemmas 8 and 13. ◻

Remark 1. With Corollary 3 proven, the only missing piece needed to prove the absolute convergence of pressure is Lemma 15.

Let us now estimate \(\mathrm{W}_\beta(\mathbf{w})\), where \(\mathbf{w}=(\gamma_1,...,\gamma_m)\in\mathscr{C}_m\) and \(m\geq 2\). That is, let us tackle the problem of multiple contour estimates. We shall need these estimates in the proof of Lemma 15. They are also crucial to computing the decay of correlations.

Definition 9. Let \(A,B\Subset\mathbb{Z}\) be disjointed and not empty. Define \[\Phi(A,B)\mathrel{\vcenter{:}}= \sum_{\substack{x\in A\\y\in B}}4J(\lvert x-y\rvert)=\sum_{\substack{x\in A\\y\in B}}\frac{4}{\lvert x-y\rvert^\alpha}.\]

The following Lemma is a generalization of Lemma 3.5 in [25].

Lemma 14. Let \(k\geq 2\) and \(\{A_i:i\in[k]\}\) be a family of pairwise nonempty disjoint finite subsets of \(\mathbb{Z}\). Let \(B\) be another nonempty finite subset of \(\mathbb{Z}\) such that \(A_i\cap B=\emptyset\), for all \(i\in [k]\). Then \[\label{eq:32phi95abc95bound} \prod_{i=1}^k\Phi(A_i,B)\leq f_k(B)\sum_{\varphi\in \mathrm{Sym}(k)}\frac{1}{\mathrm{dist}( A_{\varphi(k)},B)^\alpha}\prod_{i=1}^{k-1}\Phi(A_{\varphi(i)},A_{\varphi(i+1)}),\tag{81}\] where \(\mathrm{Sym}(k)\) denotes the \(k\)-th symmetric group and \[\begin{align} f_k(B)=16^{\alpha k}\lvert B\rvert^k(\mathrm{diam}(B)+1)^{\alpha k}. \end{align}\]

Proof. Let \(x,y,z\in\mathbb{R}\), then \[\begin{align} \label{eq:32reverse-holder} \lvert x-y\rvert^{\alpha}&\leq (\lvert x-z\rvert+\lvert z-x\rvert)^\alpha\leq 2^{\alpha-1}\lvert x-z\rvert^\alpha+2^{\alpha-1}\lvert z-y\rvert^\alpha, \end{align}\tag{82}\] where the second inequality is a reverse Hölder inequality. For each permutation \(\varphi\in \mathrm{Sym}(k)\), define \[\begin{align} E_\varphi=\{i\leq j\in [k]\Rightarrow\lvert a_{\varphi(i)}-b_{\varphi(i)}\rvert\geq\lvert a_{\varphi(j)}-b_{\varphi(j)}\rvert\}\subset\bigtimes_{i=1}^k(A_i\times B). \end{align}\] Let \((a_1,b_1,...,a_k,b_k)\in E_\varphi\) and \(i\in[k-1]\), then, by Equation 82 , \[\lvert a_{\varphi(i)}-a_{\varphi(i+1)}\rvert^\alpha\leq \left(2\lvert a_{\varphi(i)}-b_{\varphi(i)}\rvert+\mathrm{diam}(B)\right)^\alpha\leq 2^{2\alpha-1}\lvert a_{\varphi(i)}-b_{\varphi(i)}\rvert^\alpha+2^{\alpha-1}[1+\mathrm{diam}(B)]^\alpha.\]On the other hand, since \(x,y\geq2\Rightarrow x+y\leq xy\), we obtain \[\lvert a_{\varphi(i)}-a_{\varphi(i+1)}\rvert^\alpha\leq 4^{\alpha}(1+\mathrm{diam}(B))^\alpha\lvert a_{\varphi(i)}-b_{\varphi(i)}\rvert^\alpha.\]Hence, \[\begin{align} \prod_{i=1}^{k-1}\lvert a_{\varphi(i)}-a_{\varphi(i+1)}\rvert^\alpha\leq 4^{\alpha k}(1+\mathrm{diam}(B))^{\alpha k}\prod_{i=1}^{k-1}\lvert a_{\varphi(i)}-b_{\varphi(i)}\rvert^\alpha. \end{align}\] In this case, for any \((a_1,b_1,...,a_k,b_k)\in E_\varphi\), we have \[\begin{align} \frac{\prod_{i=1}^{k-1}\lvert a_{\varphi(i)}-a_{\varphi(i+1)}\rvert^{\alpha}}{\prod_{i=1}^k\lvert a_i-b_i\rvert^{\alpha}}\leq \frac{\left(4^{\alpha}(1+\mathrm{diam}(B))^{\alpha}\right)^{k}}{\lvert a_{\varphi(k)}-b_{\varphi(k)}\rvert^\alpha}. \end{align}\]On the other hand, since \(1\leq \sum_{\varphi\in S_k} \mathbbold{1}_{E_\varphi}\), we arrive at \[\begin{align} \prod_{i=1}^k\Phi(A_i,B)\leq \left(16^{\alpha}\lvert B\rvert(1+\mathrm{diam}(B))^\alpha\right)^{k}\sum_{\varphi\in \mathrm{Sym}(k)}\sum_{\substack{a_i\in A_i\\b\in B\\i\in[k]}}\frac{1}{\lvert a_{\varphi(k)}-b\rvert^{\alpha}}\prod_{j=1}^{k-1}\frac{1}{\lvert a_{\varphi(j)}-a_{\varphi(j+1)}\lvert^\alpha}. \end{align}\]Finally, since \(\mathrm{dist}(A_{\varphi(k)},B)\leq\lvert a_{\varphi(k)}-b\rvert\) and \[\begin{align} \sum_{\substack{a_i\in A_i\\i\in[k]}}\prod_{j=1}^{k-1}\frac{1}{\lvert a_{\varphi(j)}-a_{\varphi(j+1)}\lvert^\alpha}\leq \prod_{j=1}^{k-1}\sum_{\substack{a_{\varphi(j)}\in A_{\varphi(j)}\\a_{\varphi(j+1)}\in A_{\varphi(j+1)}}}\frac{1}{\lvert a_{\varphi(j)}-a_{\varphi(j+1)}\rvert^{\alpha}}, \end{align}\]we arrive at Equation 81 . The lemma has thus been proved. ◻

Corollary 4. For all \(k\geq2\), there exists \(\beta_0=\beta_0(\alpha,k)>0\) such that if \(\beta\geq\beta_0\) and \(\mathbf{w}\in\mathscr{C}_k\), then \[\mathcal{R}^k_{\beta}(\mathbf{w})\leq \frac{16^{\alpha k}\alpha}{\alpha-1}e^{-\frac{1}{4}c_2\beta}\sum_{\varphi\in \mathrm{Sym}(k)} H\left(\gamma_{\varphi(k)}\right)\prod_{i=1}^{k-1}\Phi\left(\gamma_{\varphi(i)},\gamma_{\varphi(i+1)}\right),\]where we recall the definition of the contracting functions \(\mathcal{R}^k_\beta\) in Equation 67 .

Proof. Introducing equivalence classes as in the proof of Lemma 13, we have \[\mathcal{R}^k_ {\beta}(\gamma_1,...,\gamma_k)=\sum_{[\gamma]\in\overline{\mathscr{C}}}e^{-\beta H([\gamma])}\sum_{\gamma\in [\gamma]}\prod_{i=1}^k\Phi(\gamma_i,\gamma)\mathbbold{1}_{\gamma_i \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma}.\] On the other hand, by Lemma 14, we obtain \[\prod_{i=1}^k\Phi(\gamma_i,\gamma)\mathbbold{1}_{\gamma_i \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma}\leq f_k([\gamma])\sum_{\varphi\in \mathrm{Sym}(k)}\frac{\mathbbold{1}_{\gamma_{\varphi(k)} \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma}}{\mathrm{dist}\left(\mathrm{I}_-(\gamma_{\varphi(k)}),\mathrm{I}_-(\gamma)\right)^\alpha}\prod_{j=1}^{k-1}\Phi\left(\gamma_{\varphi(j)},\gamma_{\varphi(j+1)}\right),\] where, by Hypothesis 2, \[\begin{align} f_k([\gamma])&=16^{\alpha k}\mathrm{diam}([\gamma])^{2\alpha k}\leq 16^{\alpha k}e^{2\alpha c_0kH([\gamma])} \end{align}\]For any \(l\in[k]\) and \(n\in\mathbb{N}\), there are at most \(2\lvert \gamma_l\rvert\) elements of \([\gamma]\) that satisfy \(\mathrm{dist}(\mathrm{I}_-(\gamma_l),\mathrm{I}_-(\gamma))=n\). Hence, \[\sum_{\gamma\in[\gamma]}\frac{\mathbbold{1}_{\gamma_l \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma}}{\mathrm{dist}(\mathrm{I}_-(\gamma_l),\mathrm{I}_-(\gamma))^\alpha}\leq 2\lvert\gamma_l\rvert\sum_{n=1}^\infty\frac{1}{n^\alpha}\leq\frac{\alpha}{\alpha-1}H(\gamma_l).\]Finally, we arrive at \[\begin{align} \mathcal{R}_{\beta}^k(\gamma_1,...,\gamma_k)&\leq\frac{16^{\alpha k}\alpha}{\alpha-1}\sum_{\varphi\in \mathrm{Sym}(k)}\prod_{j=1}^{k-1}\Phi(\gamma_{\varphi(j)},\gamma_{\varphi(j+1)})\sum_{[\gamma]\in\overline{\mathscr{C}}}e^{-(\beta-2\alpha c_0k)H([\gamma])} H\left(\gamma_{\varphi(k)}\right)\\ &\leq \frac{16^{\alpha k}\alpha}{\alpha-1}e^{-\frac{1}{2}c_2\beta}\sum_{\varphi\in \mathrm{Sym}(k)} H\left(\gamma_{\varphi(k)}\right)\prod_{j=1}^{k-1}\Phi\left(\gamma_{\varphi(j)},\gamma_{\varphi(j+1)}\right), \end{align}\]where the last inequality is a consequence of Hypothesis 3 and \(\beta\geq 4\alpha c_0k\). The corollary has thus been proved. ◻

Corollary 5. There exists \(\beta_0>0\) such that if \(\beta\geq\beta_0\) and \((\gamma_1,\gamma_2)\in\mathscr{C}_2\), then \[\mathrm{W}_\beta(\gamma_1,\gamma_2)\leq\Phi(\gamma_1,\gamma_2)e^{-\frac{1}{4}\beta [H(\gamma_1)+H(\gamma_2)]}.\]

Proof. This is a direct consequence of Corollaries 2 and 4. ◻

Corollary 6. Let \(m\geq 3\). Then, there exists \(\beta_0=\beta_0(\alpha,m)>0\) such that if \(\beta\geq\beta_0\) and \(\mathbf{w}\in\mathscr{C}_m\), then \[\mathrm{W}_\beta(\mathbf{w})\leq \sum_{T\in\mathbbold{T}_m}\prod_{i=1}^me^{-\beta H(\gamma_i)/2^{m}}\prod_{\{j,k\}\in E(T)}\Phi(\gamma_j,\gamma_k)\]

Proof. This is a direct consequence of Theorem 4 and Corollary 4. ◻

5.3 Convergence of the pressure↩︎

In this subsection, we prove the last lemma necessary to obtain the absolute convergence of the pressure. It is equivalent to Proposition 5.3 in [17] and Proposition 2.20 in [25]. We conclude with the proof of Theorem 1.

Lemma 15. There exists \(\beta_0>0\) such that if \(\beta\geq\beta_0\) and \(\Gamma_1\in\mathscr{P}\), then \[\label{eq:32notjoin-leaf-operator-satisfies-hypothesis} \overline{\mathcal{L}}_\beta(\Gamma_1)=\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\Gamma_1\not\sim\Gamma}\bar{z}_\beta(\Gamma)\leq 4e^{-c_2\beta/8}\sum_{\gamma_1\in\Gamma_1} H(\gamma_1).\tag{83}\]

Proof. We have \[\begin{align} \sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\Gamma_1\not\sim\Gamma}\bar{z}_\beta(\Gamma)&=\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\left\{\substack{\exists\gamma_1\in\Gamma_1,\gamma\in\Gamma\\\textrm{s.t. }\gamma_1\not\sim\gamma}\right\}}\bar{z}_\beta(\Gamma)+\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\left\{\substack{\forall\gamma_1\in\Gamma_1,\gamma\in\Gamma\\\Rightarrow\gamma_1\sim\gamma}\right\}}\mathbbold{1}_{\Gamma_1\not\sim\Gamma}\bar{z}_\beta(\Gamma). \end{align}\] Let us begin by dealing with the first sum. Recalling Equation 79 , Corollary 3 and introducing equivalence classes of contours under integer translations, we obtain \[\mathrm{S}_0\mathrel{\vcenter{:}}=\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\left\{\substack{\exists\gamma_1\in\Gamma_1,\gamma\in\Gamma\\\textrm{s.t. }\gamma_1\not\sim\gamma}\right\}}\bar{z}_\beta(\Gamma)\leq\sum_{\gamma_1\in\Gamma_1}\sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{\gamma\not\sim\gamma_1}\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\gamma\in\Gamma}\bar{z}_\beta(\Gamma)\leq \sum_{\gamma_1\in\Gamma_1}\sum_{[\gamma]\in\overline{\mathscr{C}}}e^{-\frac{1}{4}\beta H([\gamma])}\sum_{\gamma\in[\gamma]}\mathbbold{1}_{\gamma\not\sim\gamma_1}.\]By Definition 1, we have \[\sum_{\gamma\in[\gamma]}\mathbbold{1}_{\gamma\not\sim\gamma_1}\leq\sum_{x^*_1\in\gamma_1}\sum_{\gamma\in[\gamma]}\mathbbold{1}_{\mathrm{dist}(x^*_1,\gamma)\leq M\mathrm{diam}(\gamma)^{3/2}}\leq 3M\lvert\gamma_1\rvert\mathrm{diam}([\gamma])^{3/2}\leq 2M\mathrm{diam}([\gamma])^{3/2}H(\gamma_1).\]Combining Hypotheses 2 and 3, choosing \(\beta\geq 12c_0+4\log(2M)\), we obtain\[\mathrm{S}_0\leq e^{-c_2\beta/8}\sum_{\gamma_1\in\Gamma_1} H(\gamma_1).\]For the second sum, we again divide it into two parts. If \(\Gamma_1,\Gamma\in\mathscr{P}\) are such that \(\Gamma_1\cap\Gamma=\emptyset\), \(\Gamma_1\sqcup\Gamma\in\mathscr{E}\) and \(\Gamma_1\not\sim\Gamma\), then one of the following assertions holds:

  1. there exist \(\gamma\in\Gamma\) and \(\gamma_1,\gamma'_1\in\Gamma_1\) such that \(\gamma_1\trianglelefteq\gamma\) and \(\gamma_1' \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma\);

  2. or, there exists \(\gamma,\gamma'\in\Gamma\) and \(\gamma_1\in\Gamma_1\) such that \(\gamma\trianglelefteq\gamma_1\) and \(\gamma' \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1\).

That is, \(\Gamma\in A_1(\Gamma_1)\cup A_2(\Gamma_1)\), where \[A_1(\Gamma_1)\mathrel{\vcenter{:}}=\{\Gamma\in\mathscr{P}:\exists\gamma\in\Gamma,\exists\gamma_1,\gamma_1'\in\Gamma_1\textrm{ s.t. }(\gamma_1\trianglelefteq\gamma)\wedge(\gamma'_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma)\},\] and \[A_2(\Gamma_1)\mathrel{\vcenter{:}}=\{\Gamma\in\mathscr{P}:\exists\gamma,\gamma'\in\Gamma,\exists\gamma_1\in\Gamma_1\textrm{ s.t. }(\gamma\trianglelefteq\gamma_1)\wedge(\gamma' \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1)\}.\] We have, by Corollary 3, \[\mathrm{S}_1\mathrel{\vcenter{:}}=\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\Gamma\in A_1(\Gamma_1)}\bar{z}_\beta(\Gamma)\leq\sum_{\gamma_1\in\Gamma_1}\sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{\gamma_1\trianglelefteq\gamma}\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\gamma\in\Gamma}\bar{z}_\beta(\Gamma)\leq\sum_{\gamma_1\in\Gamma_1}\sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{\gamma_1\trianglelefteq\gamma}e^{-\frac{1}{4}\beta H(\gamma)}.\] Introducing equivalence classes of contours under integer translations, combining Hypotheses 2 and 3 and choosing \(\beta\geq8c_0\), we obtain \[\begin{align} \mathrm{S}_1\leq \sum_{\gamma_1\in\Gamma_1}\sum_{[\gamma]\in\overline{\mathscr{C}}}e^{-\frac{1}{4}\beta H([\gamma])}\sum_{\gamma\in[\gamma]}\mathbbold{1}_{\gamma_1\trianglelefteq\gamma}\leq\sum_{\gamma_1\in\Gamma_1}\sum_{[\gamma]\in\overline{\mathscr{C}}}\mathrm{diam}([\gamma])e^{-\frac{1}{4}\beta H([\gamma])}\leq e^{-c_2\beta/8}\sum_{\gamma_1\in\Gamma_1} H(\gamma_1). \end{align}\] On the other hand, we have \[\mathrm{S}_2\coloneq\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\Gamma\in A_2(\Gamma_1)}\bar{z}_\beta(\Gamma)\leq \sum_{\gamma_1\in\Gamma_1}\sum_{\gamma,\gamma'\in\mathscr{C}}\mathbbold{1}_{\gamma\trianglelefteq\gamma_1}\mathbbold{1}_{\gamma' \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma'}\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{\gamma,\gamma'\in\Gamma}\bar{z}_\beta(\Gamma).\]By Corollary 5, we obtain \[\mathrm{S}_2\leq \sum_{\gamma_1\in\Gamma_1}\sum_{\gamma,\gamma'\in\mathscr{C}}\mathbbold{1}_{\gamma\trianglelefteq\gamma_1}\mathbbold{1}_{\gamma' \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}\mathbbold{1}_{\gamma \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma'}\Phi(\gamma,\gamma')e^{-\frac{1}{8}\beta[H(\gamma)+H(\gamma')]}.\]Introducing equivalence classes of contours under integer translations and bounding \(\Phi(\gamma,\gamma')\) under the constraint \(\mathbbold{1}_{\gamma\trianglelefteq\gamma_1}\mathbbold{1}_{\gamma' \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_1}\), we obtain \[\mathrm{S}_2\leq\sum_{\gamma_1\in\Gamma_1}\sum_{[\gamma],[\gamma']\in\overline{\mathscr{C}}}e^{-\frac{1}{8}\beta[H([\gamma])+H([\gamma'])]}\sum_{\substack{\gamma\in[\gamma]\\\gamma'\in[\gamma']}}\sum_{\substack{x\in\mathrm{I}_-(\gamma_1)\\y\in\mathrm{I}_-(\gamma_1)^c}}\frac{4\mathbbold{1}_{x\in\mathrm{I}_-(\gamma)}\mathbbold{1}_{y\in\mathrm{I}_-(\gamma')}}{\lvert x-y\rvert^\alpha}.\]By Hypotheses 2 and 3, choosing \(\beta\geq 16c_0\), we obtain\[\mathrm{S}_2\leq \sum_{\gamma_1\in\Gamma_1} 2H(\gamma_1)\left(\sum_{[\gamma]\in\overline{\mathscr{C}}}\mathrm{diam}([\gamma])e^{-\frac{1}{8}\beta H([\gamma])}\right)^2\leq e^{-c_2\beta/8}\sum_{\gamma_1\in\Gamma_1} 2H(\gamma_1).\]Combining the bounds for \(\mathrm{S}_i\), \(i=0,1,2\), we obtain Equation 83 and Lemma 15 has thus been proved. ◻

Corollary 7. For all \(\eta>0\), there exists \(\beta_0=\beta_0(\eta)>0\) such that, if \(\beta\geq\beta_0\) and \(\varGamma\in\mathscr{P}\), then \[\overline{\mathrm{W}}_\beta(\varGamma)\leq \eta\bar{z}_{\beta/2}(\varGamma).\]

Proof. This is a direct consequence of Lemmas 8 and 15. ◻

Proof of Theorem 1. First, let us deal with the matter of convergence in the absence of an external field. By Equation 76 , it is sufficient to apply the machinery developed in Section 4. More precisely, with Lemma 15, we transform the sum over trees of polymers into a sum over trees of contours. Once again, we apply the machinery developed in Section 4; for all the details, see Section 5.2. That is, \[\begin{align} \lvert \Lambda\rvert P_\beta(\Lambda)\leq\sum_{\Gamma\in\mathscr{P}_\Lambda}\overline{\mathrm{W}}_\beta(\Gamma)\leq\sum_{\Gamma\in\mathscr{P}_\Lambda}\bar{z}_{\beta/2}(\Gamma)\leq\sum_{\gamma\in\mathscr{C}_\Lambda}\mathrm{W}_{\beta/4}(\gamma)\leq \sum_{\gamma\in\mathscr{C}_\Lambda}e^{-\frac{1}{8}\beta H(\gamma)}. \end{align}\] The first inequality is a consequence of Equation 76 . The second inequality is a consequence of Corollary 7 (\(\eta=1\)). The third inequality is a consequence of Equation 79 . Finally, the fourth inequality is a consequence of Corollary 3.

Let \([\gamma]\) be an equivalence class of contours associated to the action of \(\mathbb{Z}\) by translation. There are at most \(\lvert\Lambda\rvert\) elements of \([\gamma]\) in \(\mathscr{C}_\Lambda\). Hence, \[\begin{align} P_\beta(\Lambda)\leq \sum_{[\gamma]\in\overline{\mathscr{C}}}e^{-\frac{1}{8}\beta H([\gamma])}\leq e^{-\frac{1}{8}c_2\beta}, \end{align}\]where the second inequality is a consequence of Hypothesis 3 (recall Equation 26 ). The first part of Theorem 1 has thus been proved.

For the second part of Theorem 1, it is a direct application of Lemma 6 and \(\lvert e^z\rvert\leq e^{\lvert z\rvert}\). Take \(r_0\mathrel{\vcenter{:}}= (8\lvert A\rvert)^{-1}\). Since \(H(\gamma)\geq 2\), for all \(\gamma\in\mathscr{C}\), we have \[\beta E_{\lvert h\rvert}(\gamma)+\log2\leq 2\beta \lvert h\rvert \lvert A\rvert+\log2\leq \frac{1}{4}\beta H(\gamma),\]where the last inequality holds if \(\beta\geq 4\log2\). That is, if \(h\in\mathbb{D}_{A,r_0}\), then \[\lvert\bar{z}_{\beta,h}(\Gamma)\rvert\leq \bar{z}_{\beta/2}(\Gamma),\] for all \(\Gamma\in\mathscr{P}\). We thus conclude the proof of Theorem 1. ◻

6 Lattice site estimates↩︎

In this section, we establish some estimates involving lattice sites necessary to compute correlations of \(n\)-points. This is done because we transform polymers into contours, and then contours into lattice sites. The approach developed in this section can be seen as an extension of the one employed by Imbrie in [11] to estimate the two-point correlations for the exponent \(\alpha=2\). We begin with an auxiliary lemma.

Lemma 16. Let \(\mathbb{L}\) be a set and \(C:\mathbb{L}^2\rightarrow\mathbb{R}_{>0}\) be a function such that \(C(a,a)=1\) for all \(a\in\mathbb{L}\). If \(k\in\mathbb{N}\), \(\mathbf{x}=(x_1,...,x_k)\in\mathbb{L}^k\) and \(y\in\mathbb{L}\), then \[\sum_{\varphi\in\mathrm{Sym}(l)}C(\bar{x}_{\varphi(l)},y)\prod_{i=1}^{l-1}C(\bar{x}_{\varphi(i)},\bar{x}_{\varphi(i+1)})\leq \sum_{\varphi\in\mathrm{Sym}(k)}C(x_{\varphi(k)},y)\prod_{j=1}^{k-1}C(x_{\varphi(j)},x_{\varphi(j+1)}),\]where all the entries in \(\bar{\mathbf{x}}=(\bar{x}_1,...,\bar{x}_l)\) are distinct and \(\{\bar{x}_i:i\in[l]\}=\{x_j:j\in [k]\}\).

Proof. Let \(\mathbf{x}\) and \(\bar{\mathbf{x}}\) satisfy the hypotheses of the lemma. We shall construct an injection from \(\mathrm{Sym}(l)\) to \(\mathrm{Sym}(k)\), denoted by \(\varphi\mapsto\tilde{\varphi}\), such that \(\bar{x}_{\varphi(l)}=x_{\tilde{\varphi}(k)}\) and \[\label{eq:32C-x-x-C-bar-x-bar-x} \prod_{i=1}^{l-1}C(\bar{x}_{\varphi(i)},\bar{x}_{\varphi(i+1)})=\prod_{j=1}^{k-1}C(x_{\tilde{\varphi}(j)},x_{\tilde{\varphi}(j+1)}).\tag{84}\] Since \(C\) is positive, the existence of such an injection proves the lemma.

For \(i\in[l]\), define \(K_i\mathrel{\vcenter{:}}=\{j\in[k]:x_j=\bar{x}_i\}\).

Denote by \(\phi_i:[\lvert K_i\rvert]\rightarrow K_i\) the bijection that orders \(K_i\); that is, \(j<j'\implies \phi_i(j)<\phi_i(j')\).

Let \(\varphi\in\mathrm{Sym}(l)\). It is convenient to write \[k[\varphi](i)\mathrel{\vcenter{:}}=\sum_{p=1}^i\lvert K_{\varphi(p)}\rvert,\] with the convention \(k[\varphi](0)=0\). Note that \(k[\varphi](l)=k\) and that \(i\mapsto k[\varphi](i)\) is strictly increasing. In order to define \(\tilde{\varphi}\), we proceed in blocks: if \(j\in[k]\), then there exists a unique \(i(j)\in[l]\) such that \(k[\varphi](i(j)-1)<j\leq k[\varphi](i(j))\). Define \[\tilde{\varphi}(j)\mathrel{\vcenter{:}}= \phi_{\varphi(i(j))}\left(j-k[\varphi](i(j))\right).\]While this may seem rather involved, the graphical description is rather straightforward, see Figure 14.

Figure 14: Here, \mathbf{x}=(a,b,c,a,a,b,b,c), \bar{\mathbf{x}}=(c,a,b) and \varphi\in\mathrm{Sym(3)} denotes the trivial permutation.

Let us first show that \(\varphi\mapsto\tilde{\varphi}\) is an injection. If \(\varphi_1\neq\varphi_2\in\mathrm{Sym}(l)\), then there exists a smallest \(i\in[l]\) such that \(\varphi_1(i)\neq\varphi_2(i)\). Thus, if \(j=k[\varphi_1](i-1)+1=k[\varphi_2](i-1)+1\), then \(\tilde{\varphi}_1(j)\in K_{\varphi_1(i)}\) and \(\tilde{\varphi}_2(j)\in K_{\varphi_2(i)}\). Since \(K_{\varphi_1(i)}\cap K_{\varphi_2(i)}=\emptyset\), we conclude that \(\varphi\mapsto\tilde{\varphi}\) is an injection.

Finally, let us now show that Equation 84 holds. If \(i\in[l]\), then \[\prod_{j=k[\varphi](i-1)+1}^{k[\varphi](i)-1}C(x_{\tilde{\varphi}(j)},x_{\tilde{\varphi}(j+1)})=\prod_{j=k[\varphi](i-1)+1}^{k[\varphi](i)-1}C(\bar{x}_{\varphi(i)},\bar{x}_{\varphi(i)})=1.\]Therefore, we obtain \[\prod_{j=1}^{k-1}C(x_{\tilde{\varphi}(j)},x_{\tilde{\varphi}(j+1)})=\prod_{i=1}^{l-1}C(x_{\tilde{\varphi}(k[\varphi](i))},x_{\tilde{\varphi}(k[\varphi](i)+1)})=\prod_{i=1}^{l-1}C(\bar{x}_{\varphi(i)},\bar{x}_{\varphi(i+1)}).\]The reader may readily check that \(\bar{x}_{\varphi(l)}=x_{\tilde{\varphi}(k)}\). The lemma has thus been proved. ◻

Of course, the interest is not an arbitrary function \(C\). Following Imbrie in [11], we define\[C(x,y)=\begin{cases} \lvert x-y\rvert^{-\alpha}, &\textrm{if }x\neq y\\ 1, &\textrm{if }x=y. \end{cases}\]In general, the natural definition is \(C(x,y)=\mathbbold{1}_{x=y}+\mathbbold{1}_{x\neq y}J_{x,y}\).

Let us begin by proving an analogous result of Lemma 14 for lattice sites.

Lemma 17. Let \(k\geq 2\). Let \(\mathbf{x}=(x_1,...,x_k)\in\mathbb{Z}^k\) and \(y\in\mathbb{Z}\), then \[\label{eq:32C-xyz-bound} \prod_{i=1}^kC(x_i,y)\leq 16^{\alpha k}\sum_{\varphi\in \mathrm{Sym}(k)}C(x_{\varphi(k)},y)\prod_{i=1}^{k-1}C(x_{\varphi(i)},x_{\varphi(i+1)}).\tag{85}\] In particular, we have \[\label{eq:32triple-c-estimate} \sum_{y}\prod_{i=1}^kC(x_i,y)\leq 16^{\alpha k}(1+2\zeta(\alpha))\sum_{\varphi\in \mathrm{Sym}(k)}\prod_{i=1}^{k-1}C(x_{\varphi(i)},x_{\varphi(i+1)}).\tag{86}\]

Remark 1. If \(\{x_i\}\) and \(y\) are all distinct, this is simply Lemma 14.

Proof. Let \(\bar{\mathbf{x}}=(\bar{x}_1,...,\bar{x}_l)\) such that \(i\neq j\Rightarrow\bar{x}_i\neq\bar{x}_j\) and \(\{x_i:i\in[k]\}=\{\bar{x}_i:i\in[l]\}\). Then \[\begin{align} \prod_{i=1}^kC(x_i,y)=\prod_{j=1}^lC(\bar{x}_j,y)^{k_j}, \end{align}\]where \(k_j=\lvert\{i\in[k]:x_i=\bar{x}_j\}\rvert\). Hence, \[\prod_{i=1}^kC(x_i,y)\leq\prod_{j=1}^lC(\bar{x}_j,y)=\prod_{j=1}^l\frac{\mathbbold{1}_{\bar{x}_j\neq y}}{\lvert \bar{x}_j-y\rvert^\alpha}+\sum_{i=1}^l\mathbbold{1}_{\bar{x}_i=y}\prod_{j\neq i}\frac{\mathbbold{1}_{\bar{x}_j\neq y }}{\lvert y-\bar{x}_j\rvert^\alpha}.\]By Lemma 14, we bound the first term by \[\prod_{j=1}^l\frac{\mathbbold{1}_{\bar{x}_j\neq y}}{\lvert \bar{x}_j-y\rvert^\alpha}\leq 16^{\alpha l}\prod_{i=1}^l\mathbbold{1}_{\bar{x}_i\neq y}\sum_{\varphi\in \mathrm{Sym}(l)}C(\bar{x}_{\varphi(l)},y)\prod_{j=1}^{l-1}C(\bar{x}_{\varphi(j)},\bar{x}_{\varphi(j+1)}).\]Let us now bound the second term. For \(i\in[l]\), define \(\mathbf{z}^i=(\bar{x}_1,...,\bar{x}_{i-1},\bar{x}_{i+1},...,\bar{x}_l)\). Thus, the same Lemma 14 implies that \[\mathbbold{1}_{\bar{x}_i=y}\prod_{j=1}^{l-1}\frac{\mathbbold{1}_{z^i_j\neq y}}{\lvert z^i_j-y\rvert^{\alpha}}\leq 16^{\alpha(l-1)}\mathbbold{1}_{\bar{x}_i=y}\prod_{j\neq i}\mathbbold{1}_{\bar{x}_j\neq y}\sum_{\varphi\in\mathrm{Sym}(l-1)}\prod_{j=1}^{l-2}C(z^{i}_{\varphi(j)},z^i_{\varphi(j+1)})\cdot C(z^i_{\varphi(l-1)},\bar{x}_i).\]There is a canonical bijection \(\mathrm{Sym}(l-1)\ni\varphi\leftrightarrow\tilde{\varphi}\in\mathrm{Sym}(l)\) given by \[\tilde{\varphi}(j)\mathrel{\vcenter{:}}=\begin{cases} i, &\textrm{if }j=l;\\ \varphi(j), &\textrm{if }\varphi(j)<i;\\ \varphi(j)+1, &\textrm{if }\varphi(j)\geq i. \end{cases}\]Recalling that \(\mathbbold{1}_{\bar{x}_i=y}=\mathbbold{1}_{\bar{x}_i=y}C(\bar{x}_i,y)\), the right-hand side can be written as\[16^{\alpha (l-1)}\mathbbold{1}_{\bar{x}_i=y}\prod_{j\neq i}\mathbbold{1}_{\bar{x}_i\neq\bar{x}_j}\sum_{i=1}^l\sum_{\substack{\varphi\in\mathrm{Sym}(l)\\\varphi(l)=i}}C(\bar{x}_{\varphi(l)},y)\prod_{j=1}^{l-1}C(\bar{x}_{\varphi(j)},\bar{x}_{\varphi(j+1)}).\]Finally, since \[\prod_{i=1}^l\mathbbold{1}_{\bar{x}_i\neq y}+\sum_{i=1}^l\mathbbold{1}_{\bar{x}_i=y}\prod_{j\neq i}\mathbbold{1}_{\bar{x}_i\neq \bar{x}_j}\leq 1,\]we obtain that\[\prod_{i=1}^kC(x_i,y)\leq 16^{\alpha l}\sum_{\varphi\in\mathrm{Sym}(l)}C(\bar{x}_{\varphi(l)},y)\prod_{i=1}^{l-1} C(\bar{x}_{\varphi(i)},\bar{x}_{\varphi(i+1)}).\]By Lemma 16, the proof is finished. ◻

Let \(m<n\) be natural numbers, \(T\in\overline{\mathbbold{T}}_{n,m}\) (recall Notation 1) and \(\mathbf{x}=(x_1,...,x_m)\in\mathbb{Z}^m\). We want to estimate \[w(\mathbf{x},T)\mathrel{\vcenter{:}}=\sum_{\mathbf{y}\in\mathbb{Z}^n}\mathbbold{1}_{\mathbf{y}|_{[m]}=\mathbf{x}}\prod_{\{p,q\}\in E(T)}C(y_p,y_q).\]The algorithm is rather straightforward: we first sum over \(y_n\), then \(y_{n-1}\), until \(y_{m+1}\).

We proceed as in Section 4.3. Write \(\delta_n\mathrel{\vcenter{:}}=\delta_T(n)\) and \(N_T(n)=\{l_1,...,l_{\delta_n}\}\), where \(l_1<...<l_{\delta_n}\). By Lemma 17, we have \[w(\mathbf{x},T)\leq 16^{\alpha \delta_n}(1+2\zeta(\alpha))\sum_{\substack{\mathbf{y}\in\mathbb{Z}^{n-1}\\\mathbf{y}|_{[m]}=\mathbf{x}}}\prod_{\{p,q\}\in E(T)\setminus E_T(n)}C(y_p,y_q)\sum_{\varphi_n\in\mathrm{Sym}(\delta_n)}\prod_{r=1}^{\delta_n-1}C(y_{l_{\varphi(r)}},y_{l_{\varphi(r+1)}}).\]Recall that \(\varphi\mapsto T_\varphi\in\overline{\mathbbold{T}}_{n-1,m}\) is two-to-one, where \[E(T_\varphi)\mathrel{\vcenter{:}}=(E(T)\setminus E_T(n))\cup\{\{l_{\varphi(i)},l_{\varphi(i+1)}\}:i\in[\delta_n-1]\},\] and that \(\delta_n\leq n-1\). Therefore, \[\label{eq:32iterative-step-contracting-tree-of-points} w(\mathbf{x},T)\leq 16^{\alpha(n-1)}2(1+2\zeta(\alpha))\sum_{T'\in\overline{\mathbbold{T}}_{n-1,m}}w(\mathbf{x},T')\leq 16^{\alpha n}(1+2\zeta(\alpha))\sum_{T'\in\overline{\mathbbold{T}}_{n-1,m}}w(\mathbf{x},T').\tag{87}\]

Lemma 18. Let \(m<n\) be natural numbers, \(T\in\overline{\mathbbold{T}}_{n,m}\) and \(\mathbf{x}=(x_1,...,x_m)\in\mathbb{Z}^m\). Then, \[\label{eq:32contracting-tree-of-points} w(\mathbf{x},T)\leq \left[\prod_{l=m+1}^{n-1}\lvert\overline{\mathbbold{T}}_{l,m}\rvert\right] 8^{\alpha n^2}(1+2\zeta(\alpha))^n\sum_{T'\in\mathbbold{T}_m}\prod_{\{p,q\}\in E(T')}C(x_p,x_q).\tag{88}\]

Proof. Equation 88 follows from a direct iteration of Equation 87 . ◻

Let us finish this section with the description of a computationally more convenient way to detect incompatibility between polymers in terms of contours and relate it to estimates involving the function \(C(\cdot,\cdot)\).

Definition 10. We define \[\label{eq:32incompatibility-distance-contours-polymers} \mathscr{A}\mathrel{\vcenter{:}}=\left\{(\gamma_1,\gamma_2)\in\mathscr{C}^2:\mathrm{dist}(\mathrm{I}_-(\gamma_1),\mathrm{I}_-(\gamma_2))\leq 2M\min\{\mathrm{diam}(\gamma_1),\mathrm{diam}(\gamma_2)\}^{3/2}\right\}.\tag{89}\]

Lemma 19. Let \(\Gamma_1,\Gamma_2\in\mathscr{P}\). If \(\Gamma_1\not\sim\Gamma_2\), then there exists \((\gamma_1,\gamma_2)\in(\Gamma_1\times\Gamma_2)\cap\mathscr{A}\).

Proof. By Remark 1, there are two cases to consider. First, if there exist \(\gamma_1\in\Gamma_1\) and \(\gamma_2\in\Gamma_2\) such that \(\gamma_1\not\sim\gamma_2\), then \[\begin{align} \mathrm{dist}(\mathrm{I}_-(\gamma_1),\mathrm{I}_-(\gamma_2))&\leq 1+\mathrm{dist}(\gamma_1,\gamma_2)\leq 1+M\min\{\mathrm{diam}(\gamma_1),\mathrm{diam(\gamma_2)}\}^{3/2}\\&\leq 2M\min\{\mathrm{diam}(\gamma_1),\mathrm{diam(\gamma_2)}\}^{3/2}. \end{align}\]That is, \((\gamma_1,\gamma_2)\in\mathscr{A}\). If \(\Gamma_1\not\sim\Gamma_2\) and \((\gamma_1,\gamma_2)\in\Gamma_1\times\Gamma_2\Rightarrow \gamma_1\sim\gamma_2\), then there exist \(\gamma_1\in\Gamma_1\) and \(\gamma_2\in\Gamma_2\) such that \(\mathrm{I}_-(\gamma_1)\cap\mathrm{I}_-(\gamma_2)\neq\emptyset\). That is, \((\gamma_1,\gamma_2)\in\mathscr{A}\). The lemma has thus been proved. ◻

Let us now prove the estimates that allow us to translate estimates in terms of contours into estimates in terms of lattice sites.

Lemma 20. Let \(x,y\in\mathbb{Z}\).

  1. If \(\gamma\in\mathscr{C}\) and \(x,y\in\mathrm{I}_-(\gamma)\), then \(1\leq e^{\alpha c_0H(\gamma)} C(x,y)\).

  2. If \((\gamma_1,\gamma_2)\in\mathscr{C}_2\), \(x\in\mathrm{I}_-(\gamma_1)\) and \(y\in\mathrm{I}_-(\gamma_2)\), then \(\Phi(\gamma_1,\gamma_2)\leq 4^{\alpha}e^{2\alpha c_0[H(\gamma_1)+H(\gamma_2)]}C(x,y)\).

  3. If \((\gamma_1,\gamma_2)\in\mathscr{A}\), \(x\in\mathrm{I}_-(\gamma_1)\) and \(y\in\mathrm{I}_-(\gamma_2)\), then \(1\leq (4M)^\alpha e^{\alpha c_0[H(\gamma_1)+H(\gamma_2)]}C(x,y).\)

Proof. The first bound follows directly from Hypothesis 2.

Let us prove the second one. In general, we have \[\label{eq:32trivial-phi-bound} \Phi(\gamma_1,\gamma_2)\leq\frac{\mathrm{diam}(\gamma_1)\mathrm{diam}(\gamma_2)}{\mathrm{dist}(\gamma_1,\gamma_2)^\alpha}.\tag{90}\] On the other hand, by the triangular inequality, we have \[\lvert x-y\rvert\leq \mathrm{diam}(\gamma_1)+\mathrm{dist}(\gamma_1,\gamma_2)+\mathrm{diam}(\gamma_2).\]Dividing and multiplying by \(2\) and using \(x,y\geq 2\Rightarrow x+y\leq xy\), we arrive at \[\label{eq:32crude-x-y-bound} \lvert x-y\rvert^\alpha\leq 2^{2\alpha}\mathrm{dist}(\gamma_1,\gamma_2)^\alpha\mathrm{diam}(\gamma_1)^\alpha\mathrm{diam}(\gamma_2)^\alpha.\tag{91}\] Combining Hypothesis 2 with Equations 90 and 91 , we obtain the second bound.

Finally, let us prove the third bound. Let \((\gamma_1,\gamma_2)\in\mathscr{A}\), \(x_1\in\mathrm{I}_-(\gamma_1)\) and \(x_2\in\mathrm{I}_-(\gamma_2)\). Then \[\lvert x-y\rvert\leq \mathrm{diam}(\gamma_1)+\mathrm{dist}(\mathrm{I}_-(\gamma_1),\mathrm{I}_-(\gamma_2))+\mathrm{diam}(\gamma_2).\]By the definition of \(\mathcal{A}\), \[\mathrm{dist}(\mathrm{I}_-(\gamma_1),\mathrm{I}_-(\gamma_2))\leq 2M\mathrm{diam}(\gamma_1)^{3/4}\mathrm{diam}(\gamma_2)^{3/4}.\] A clever application of \(x,y\geq 2\Rightarrow x+y\leq xy\) gives us \[M\mathrm{diam}(\gamma_1)^{3/4}\mathrm{diam}(\gamma_2)^{3/4}+\mathrm{diam}(\gamma_1)\leq 2M\mathrm{diam}(\gamma_1)\mathrm{diam}(\gamma_2)^{3/4}.\] Thus, \[\lvert x-y\rvert\leq M\mathrm{diam}(\gamma_1)^{3/4}\mathrm{diam}(\gamma_2)^{3/4}(2\mathrm{diam}(\gamma_1)^{1/4}+2\mathrm{diam}(\gamma_2)^{1/4}).\]Finally, combining another application of \(x,y\geq 2\Rightarrow x+y\leq xy\) with Hypothesis 2, we obtain the third bound. The lemma has thus been proved. ◻

Corollary 8. Let \(\varGamma_1,\varGamma_2\in\mathscr{P}\), then \[\label{eq:32polymer-incompatibility-point-estimate} \mathbbold{1}_{\varGamma_1\not\sim\varGamma_2}\leq (4M)^\alpha\sum_{\substack{\gamma_1\in\Gamma_1\\\gamma_2\in\Gamma_2}}e^{\alpha c_0 H(\gamma_1)}e^{\alpha c_0 H(\gamma_2)}\sum_{\substack{x_1\in\mathrm{I}_-(\gamma_1)\\x_2\in\mathrm{I}_-(\gamma_2)}}C(x_1,x_2)\eqqcolon(4M)^\alpha\tilde{\mathrm{E}}_\alpha(\varGamma_1,\varGamma_2).\tag{92}\]

Proof. The corollary follows from Lemmas 19 and 20 (item 3): \[\begin{align} \mathbbold{1}_{\varGamma_1\not\sim\varGamma_2}&\leq \sum_{\substack{\gamma_1\in\varGamma_1\\\gamma_2\in\varGamma_2}}\mathbbold{1}_{\mathscr{A}}(\gamma_1,\gamma_2)\leq\sum_{\substack{\gamma_1\in\varGamma_1\\\gamma_2\in\varGamma_2}}\sum_{\substack{x_1\in\mathrm{I}_-(\gamma_1)\\x_2\in\mathrm{I}_-(\gamma_2)}}\mathbbold{1}_{\mathscr{A}}(\gamma_1,\gamma_2)\mathbbold{1}_{x_1\in\mathrm{I}_-(\gamma_1)}\mathbbold{1}_{x_2\in\mathrm{I}_-(\gamma_2)}\\ &\leq (4M)^\alpha\sum_{\substack{\gamma_1\in\varGamma_1\\\gamma_2\in\varGamma_2}}e^{\alpha c_0 H(\gamma_1)}e^{\alpha c_0H(\gamma_2)}\sum_{\substack{x_1\in\gamma_1\\x_2\in\gamma_2}}C(x_1,x_2). \end{align}\] ◻

All the ingredients needed to give an upper bound to the decay rate of the \(n\)-point correlation functions have now been described.

7 Correlations of n-points↩︎

In this section, we apply the cluster expansion developed in Sections 3 and 5 and the site estimates proved in Section 6 to obtain an upper bound to the \(n\)-point correlations of the long-range one-dimensional ferromagnetic Ising model with a homogeneous boundary condition. That is, we are interested in estimating \[\label{eq:32n-point-correlations} \langle :\sigma_A:\rangle^+_{\Lambda,\beta}\mathrel{\vcenter{:}}=\left\langle\prod_{a\in A}\left(\sigma_a-\langle\sigma_a\rangle^+_{\Lambda,\beta}\right)\right\rangle^+_{\Lambda,\beta},\tag{93}\] where \(A\) is an arbitrary subset of \(\Lambda\). Since \(A\) is finite, by Theorem 1, the function \[\mathbb{D}^A_{r_0}\ni (h_a)_{a\in A}\mapsto \sum_{n=1}^\infty\frac{1}{n!}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n_\Lambda}\phi_n(\mathbf{\Gamma})\prod_{i=1}^nz^+_{\beta,h}(\Gamma_i)\]is analytic. Therefore, for any \(a\in A\), we have \[\langle\sigma_a\rangle^+_{\Lambda,\beta}=\frac{1}{(-2\beta)}\partial_a|_{h=0}\log Z^+_{\Lambda,\beta, h},\]where \(\partial_a\equiv\partial/\partial h_a\). A standard argument (see Chapter 2 of [33]) gives us \[\label{eq:32A-point-wick-product} \langle :\sigma_A:\rangle^+_{\Lambda,\beta}=\frac{1}{(-2\beta)^{\lvert A\rvert}}\partial_A|_{h=0}\log Z^+_{\Lambda,\beta,h},\tag{94}\] where \(\partial_A\equiv \prod_{a\in A}\partial_a\).

Let us now define the measures induced by polymers that are useful in describing the \(n\)-point correlation functions. Let \(\Gamma\in\mathscr{P}\) be a polymer. We define \[\label{eq:32union-of-interiors} \mathscr{U}(\Gamma)\mathrel{\vcenter{:}}=\left\{\overline{\Gamma}\in\mathscr{E}:\overline{\Gamma}=\sqcup_{\gamma\in\Gamma}\iota^*_\gamma(\overline{\Gamma})\right\}.\tag{95}\] Recalling Notation 1 and Equations 28 and 29 , we define two families of atomic probability measures on \(\mathscr{U}(\Gamma)\) indexed by \(\beta\) in the following way: \[\label{eq:32not-star-polymer-measure} \begin{align} &\mu_{\Gamma,\beta}(\overline{\Gamma})\mathrel{\vcenter{:}}=\frac{1}{Z_\beta(\Gamma)}\sum_{G\in\mathbbold{G}(\Gamma)}\prod_{\gamma\in\Gamma}e^{-\beta H_\gamma(\overline{\Gamma})}\prod_{e\in E(G)}\left[e^{\beta\Phi_e(\overline{\Gamma})}-1\right],\\ &\mu^*_{\Gamma,\beta}(\overline{\Gamma})\mathrel{\vcenter{:}}=\frac{1}{Z^*_{\beta}(\Gamma)}\prod_{\gamma\in\Gamma}e^{-\beta H^*_\gamma(\overline{\Gamma})}. \end{align}\tag{96}\] Let \(A\in\mathrm{P}(\mathbb{Z})\). An easy exercise of calculus gives us \[\begin{align} &\partial_A|_{h=0}Z_{\beta,h}(\Gamma)=\mathbbold{1}_{A\subset\mathrm{I}_-(\Gamma)}(-2\beta)^{\lvert A\rvert}\mu_{\Gamma,\beta}\{A\cap\mathrm{I}_-(\overline{\Gamma})=\emptyset\}Z_\beta(\Gamma),\\ &\partial_A|_{h=0}Z^*_{\beta,h}(\Gamma)=\mathbbold{1}_{A\subset\mathrm{I}_-(\Gamma)}(-2\beta)^{\lvert A\rvert}\mu^*_{\Gamma,\beta}\{A\subset\mathrm{I}_-(\overline{\Gamma})\}Z^*_{\beta}(\Gamma). \end{align}\] Therefore, \[\partial_A|_{h=0}z^+_{\beta,h}(\Gamma)=\mathbbold{1}_{A\subset\mathrm{I}_-(\Gamma)} (-2\beta)^{\lvert A\rvert}\left[\mu_{\Gamma,\beta}\{A\cap\mathrm{I}_-(\overline{\Gamma})=\emptyset\}-\mu^*_{\Gamma,\beta }\{A\subset\mathrm{I}_-(\overline{\Gamma})\}\right]z^+_{\beta}(\Gamma).\] Let \(\mathbf{\Gamma}=(\Gamma_1,...,\Gamma_n)\in\mathscr{P}^n\), we obtain \[\label{eq:32derivative-of-the-activity} \left\lvert \frac{\partial}{\partial h_A}\left(\prod_{i=1}^nz^+_{\beta,h}(\Gamma_i)\right)_{h=0}\right\rvert\leq (2\beta)^{\lvert A\rvert}\sum_{\mathbf{A}}\prod_{i=1}^n\mathbbold{1}_{A_i\subset\mathrm{I}_-(\Gamma_i)}z^+_\beta(\Gamma_i),\tag{97}\] where \(\mathbf{A}\in\{(A_1,...,A_n)\in(\mathrm{P}(\mathbb{Z}))^n:(i\neq j\Rightarrow A_i\cap A_j=\emptyset)\wedge(\sqcup_iA_i=A)\}\).

For every partition \(P\in\mathbbold{P}(A)\), define \[\label{eq:32W-of-a-partition-definition} \overline{\mathrm{W}}_\beta(P)\mathrel{\vcenter{:}}=\sum_{n=m}^\infty\frac{1}{(n-m)!}\sum_{T\in\mathbbold{T}_n}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n}\prod_{i=1}^m\mathbbold{1}_{A_i\subset\mathrm{I}_-(\Gamma_i)}\prod_{j=1}^n\bar{z}_\beta(\Gamma_i)\prod_{\{p,q\}\in E(T)}\mathbbold{1}_{\Gamma_p\not\sim\Gamma_q},\tag{98}\] where \((A_1,...,A_{\lvert P\rvert})\) is the canonical ordering of \(P\) given by \(i<j\) implies \(\min A_i<\min A_j\).

By Lemma 6 and Equations 74 and 94 , we obtain \[\label{eq:32wick-product-bounded-by-sums-over-trees} \lvert\langle:\sigma_A:\rangle^+_{\Lambda,\beta}\rvert\leq\sum_{P\in\mathbbold{P}(A)}\overline{\mathrm{W}}_\beta(P),\tag{99}\] which we naturally identified with a family of sums of trees. Our objective is to describe it in terms of contours, which we can control very well (see Section 5.2).

If \(\lvert P\rvert\geq 2\), by a similar argument that proved Lemma 12, we have \[\label{eq:32W-overline-partition-star-trees} \overline{\mathrm{W}}_\beta(P)\leq\sum_{n=m}^{2m-2}\sum_{T\in\mathbbold{T}^*_{n,m}}\sum_{\mathbf{\Gamma}\in\mathscr{P}^n}\prod_{i=1}^m\mathbbold{1}_{A_i\subset\mathrm{I}_-(\Gamma_i)}\prod_{j=1}^n\overline{\mathrm{W}}_\beta(\Gamma_i)\prod_{\{p,q\}\in E(T)}\overline{\mathrm{E}}_\beta(\Gamma_p,\Gamma_q),\tag{100}\] where, for all \((\varGamma_1,\varGamma_2)\in\mathscr{P}^2\), we recall Equation 56 : \[\label{eq:32concrete-total-edge-function} \overline{\mathrm{E}}_\beta(\varGamma_1,\varGamma_2)\mathrel{\vcenter{:}}=\mathbbold{1}_{\varGamma_1\not\sim\varGamma_2}+\sum_{n=1}^\infty\sum_{\mathbf{\Gamma}\in\mathscr{P}^{n+2}}\mathbbold{1}_{\Gamma_1=\varGamma_1}\mathbbold{1}_{\Gamma_{n+2}=\varGamma_2}\prod_{i=2}^{n+1}\overline{\mathrm{W}}_\beta(\Gamma_i)\prod_{j=1}^{n+1}\mathbbold{1}_{\Gamma_j\not\sim\Gamma_{j+1}}.\tag{101}\] Evidently, the first task at hand is to obtain a decent bound for Equation 101 . Once again, we define for \(n\in\mathbb{N}\) \[\label{eq:32n-concrete-total-edge-function} \overline{\mathrm{E}}_{\beta,n}(\varGamma_1,\varGamma_2)\mathrel{\vcenter{:}}=\sum_{\mathbf{\Gamma}\in\mathscr{P}^{n+2}}\mathbbold{1}_{\Gamma_1=\varGamma_1}\mathbbold{1}_{\Gamma_{n+2}=\varGamma_2}\prod_{i=2}^{n+1}\overline{\mathrm{W}}_\beta(\Gamma_i)\prod_{j=1}^{n+1}\mathbbold{1}_{\Gamma_j\not\sim\Gamma_{j+1}}.\tag{102}\]

7.1 Total edge function bounds↩︎

Let us now bound the total edge function for polymer trees (see Equation 56 ). Let \(\Gamma_0,\Gamma_1,\Gamma_2\) be polymers such that \(\Gamma_0\not\sim\Gamma_1\) and \(\Gamma_1\not\sim\Gamma_2\). By Lemma 19, there exist \((\gamma_{0,1},\gamma_{1,0})\in(\Gamma_0\times\Gamma_1)\cap\mathscr{A}\) and \((\gamma_{1,2},\gamma_{2,1})\in(\Gamma_1\times\Gamma_2)\cap\mathscr{A}\). There are two cases to consider: \(\gamma_{1,0}=\gamma_{1,2}\) or \(\gamma_{1,0} \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_{1,2}\), that is, \((\gamma_{1,0},\gamma_{1,2})\in\mathscr{C}_2\). For a graphical description, see Figure 15.

Figure 15: Graphical description of the translation of polymer incompatibility in terms of contours.

In this case, we obtain\[\mathbbold{1}_{\Gamma_0\not\sim\Gamma_1}\mathbbold{1}_{\Gamma_1\not\sim\Gamma_2}\leq\sum_{\Upsilon}\mathbbold{1}_\mathscr{A}(\gamma_{0,1},\gamma_{1,0})\tilde{\mathbbold{1}}(\gamma_{1,0},\gamma_{1,2})\mathbbold{1}_{\mathscr{A}}(\gamma_{1,2},\gamma_{2,1})\mathbbold{1}_{\gamma_{0,1}\in\Gamma_0}\mathbbold{1}_{\gamma_{1,0},\gamma_{1,2}\in\Gamma_1}\mathbbold{1}_{\gamma_{2,1}\in\Gamma_2},\]where \(\Upsilon=(\gamma_{0,1},\gamma_{1,0},\gamma_{1,2},\gamma_{2,1})\in\mathscr{C}^4\) and \(\tilde{\mathbbold{1}}:\mathscr{C}^2\rightarrow\{0,1\}\) is given by \[\label{eq:32tilde-indicator-definition} \tilde{\mathbbold{1}}(\gamma_1,\gamma_2)\mathrel{\vcenter{:}}=\mathbbold{1}_{\gamma_1=\gamma_2}+\mathbbold{1}_{\gamma_1 \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_2}.\tag{103}\] This generalizes rather easily to linear trees: let \((\Gamma_0,\Gamma_1,...,\Gamma_{n+1})\in\mathscr{P}^{n+2}\), then \[\label{eq:32first-try-linear-incompatibilities} \prod_{i=0}^{n}\mathbbold{1}_{\Gamma_i\not\sim\Gamma_{i+1}}\leq\sum_{\Upsilon}\prod_{i=0}^{n}\mathbbold{1}_{\mathscr{A}}(\gamma_{i,i+1},\gamma_{i+1,i})\prod_{j=1}^{n}\tilde{\mathbbold{1}}(\gamma_{j,j-1},\gamma_{j,j+1})\prod_{k=0}^{n}\mathbbold{1}_{\gamma_{k,k+1}\in\Gamma_k}\mathbbold{1}_{\gamma_{k+1,k}\in\Gamma_{k+1}},\tag{104}\] where \(\Upsilon=(\gamma_{0,1},\gamma_{1,2},...,\gamma_{n,n+1},\gamma_{n+1,n})\in\mathscr{C}^{2n+2}\). Writing \[\label{eq:32J-tilde-indicator-function} \prod_{j=1}^{n}\tilde{\mathbbold{1}}(\gamma_{j,j-1},\gamma_{j,j+1})=\sum_{J\subset [n]}\prod_{j_1\in J}\mathbbold{1}_{\gamma_{j_1-1,j_1} \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_{j_1,j_1+1}}\prod_{j_2\in[n]\setminus J}\mathbbold{1}_{\gamma_{j_2,j_2-1}=\gamma_{j_2,j_2+1}}\eqqcolon\sum_{J\subset[n]}\tilde{\mathbbold{1}}_J(\Upsilon),\tag{105}\] we obtain \[\label{eq:32linear-incompatibilities-final-formula} \prod_{i=0}^{n}\mathbbold{1}_{\Gamma_i\not\sim\Gamma_{i+1}}\leq \sum_{J\subset [n]}\sum_{\Upsilon}\tilde{\mathbbold{1}}_J(\Upsilon)\tilde{\mathbbold{1}}_{\mathscr{A}}(\Upsilon)\prod_{k=0}^{n-1}\mathbbold{1}_{\gamma_{k,k+1}\in\Gamma_k}\mathbbold{1}_{\gamma_{k+1,k}\in\Gamma_{k+1}},\tag{106}\] where \[\tilde{\mathbbold{1}}_\mathscr{A}(\Upsilon)\mathrel{\vcenter{:}}=\prod_{i=0}^{n}\mathbbold{1}_{\mathscr{A}}(\gamma_{i,i+1},\gamma_{i+1,i}).\]

Lemma 21. There exists \(\beta_0>0\) such that if \(\beta\geq\beta_0\) and \(\varGamma_1,\varGamma_2\in\mathscr{P}\), then \[\label{eq:32concrete-overline-E-beta-bound} \overline{\mathrm{E}}_\beta(\varGamma_1,\varGamma_2)\leq (8M)^\alpha \tilde{\mathrm{E}}_\alpha(\varGamma_1,\varGamma_2),\tag{107}\] where we recall Equation 92 .

Remark 1. This is very similar to the type of bound announced in Equation 59 .

Proof. By Corollary 7 and Equations 102 and 106 , we have \[\overline{\mathrm{E}}_{\beta,n}(\varGamma_1,\varGamma_2)\leq\sum_{J\subset [n]}\overline{\mathrm{E}}_{\beta,J}(\varGamma_1,\varGamma_2),\]where \(\Upsilon=(\gamma_{0,1},\gamma_{1,0},...,\gamma_{n,n+1},\gamma_{n+1,n})\in\varGamma_1\times \mathscr{C}^{2n}\times\varGamma_2\), \(\mathbf{\Gamma}=(\Gamma_1,...,\Gamma_n)\in\mathscr{P}^n\) and \[\overline{\mathrm{E}}_{\beta,J}(\varGamma_1,\varGamma_2)\mathrel{\vcenter{:}}=\sum_{\Upsilon}\tilde{\mathbbold{1}}_J(\Upsilon)\tilde{\mathbbold{1}}_\mathscr{A}(\Upsilon)\sum_{\mathbf{\Gamma}}\bar{z}_{\beta/2} (\Gamma_k)\prod_{k=1}^n\mathbbold{1}_{\gamma_{k,k-1},\gamma_{k,k+1}\in\Gamma_k}.\]By Corollaries 3 and 5, we obtain \[\overline{\mathrm{E}}_{\beta, J}(\varGamma_1,\varGamma_2)\leq\sum_{\Upsilon}\overline{\mathbbold{1}}_J(\Upsilon)\tilde{\mathbbold{1}}_{\mathscr{A}}(\Upsilon)\prod_{i=1}^ne^{-\frac{1}{16}\beta[H(\gamma_{i,i-1})+H(\gamma_{i,i+1})]}\prod_{j\in J}\Phi(\gamma_{j,j-1},\gamma_{j,j+1}).\]We have thus bounded \(\overline{\mathrm{E}}_{\beta,n}\) in terms of contours.

Let us transform this bound in terms of contours into a bound in terms of lattice sites. Since \(\gamma\in\mathscr{C}\Rightarrow\lvert\mathrm{I}_-(\gamma)\rvert\geq 1\), we have \[\label{eq:32contours-to-points} 1\leq \sum_{\mathbf{x}}\prod_{k=0}^n\mathbbold{1}_{x_{k,k+1}\in\mathrm{I}_-(\gamma_{k,k+1})}\mathbbold{1}_{x_{k+1,k}\in\mathrm{I}_-(\gamma_{k+1,k})},\tag{108}\] where \(\mathbf{x}=(x_{0,1},x_{1,0},...,x_{n,n+1},x_{n+1,n})\in\mathbb{Z}^{2n+2}\).

By Lemma 20, if \(j\in J\), \(x_{j,j-1}\in\mathrm{I}_-(\gamma_{j,j-1})\) and \(x_{j,j+1}\in\mathrm{I}_-(\gamma_{j,j+1})\), then \[\label{eq:32in-J-contours-to-points} \Phi(\gamma_{j,j-1},\gamma_{j,j+1})\leq 4^{\alpha}e^{2\alpha c_0\left[H(\gamma_{j,j-1})+H(\gamma_{j,j+1})\right]}C(x_{j,j-1},x_{j,j+1}).\tag{109}\] If \(j\in [n]\setminus J\), \(x_{j,j-1}\in\mathrm{I}_-(\gamma_{j,j-1})\) and \(x_{j,j+1}\in\mathrm{I}_-(\gamma_{j,j+1})\), then \[\label{eq:32out-J-contours-to-points} 1\leq e^{\alpha c_0[H(\gamma_{j,j-1})+H(\gamma_{j,j+1})]}C(x_{j,j-1},x_{j,j+1}).\tag{110}\] If \((\gamma_{i,i+1},\gamma_{i+1,i})\in\mathscr{A}\), \(x_{i,i+1}\in\mathrm{I}_-(\gamma_{i,i+1})\) and \(x_{i+1,i}\in\mathrm{I}_-(\gamma_{i+1,i})\), then \[\label{eq:32in-scrA-contours-to-points} 1\leq (4M)^\alpha e^{\alpha c_0[H(\gamma_{i,i+1})+H(\gamma_{i+1,i})]}C(x_{i,i+1},x_{i+1,i}).\tag{111}\]

Figure 16: Condensed notation employed in this proof.

Condensing notation as in Figure 16, combining Equations 109 , 110 and 111 , we obtain \[\overline{\mathrm{E}}_{\beta,J}(\varGamma_1,\varGamma_2)\leq (8M)^{2\alpha n}\sum_{\mathbf{x}}\sum_{\Upsilon}\overline{\mathbbold{1}}_J(\Upsilon,\mathbf{x})C(\mathbf{x})\epsilon_\alpha(\Upsilon)\epsilon_\beta(\Upsilon).\] Let us now sum over \(\Upsilon'=(\gamma_{1,0},\gamma_{1,2},...,\gamma_{n,n-1},\gamma_{n,n+1})\in\mathscr{C}^{2n}\). As often done with estimates with contours, let us introduce equivalence classes of contours under integer translations. Noting that \(x_{i,j}\) belongs to the minus interior of \(\lvert\mathrm{I}_-(\gamma_{i,j})\rvert\leq e^{\alpha c_0H(\gamma_{i,j})}\) elements of \([\gamma_{i,j}]\), by Hypothesis 3, we obtain \[\label{eq:32summing-over-Upsilon-prime} \sum_{\Upsilon'}\overline{\mathbbold{1}}_J(\Upsilon,\mathbf{x})\epsilon_\beta(\Upsilon)\leq \mathbbold{1}_{x_{0,1}\in\mathrm{I}_-(\gamma_{0,1})}\mathbbold{1}_{x_{n+1,n}\in\mathrm{I}_-(\gamma_{n+1,n})}e^{-n\frac{c_2}{8}(\beta-64\alpha c_0)}.\tag{112}\] We must now sum over \(\mathbf{x}'=(x_{1,0},x_{1,2}...,x_{n,n-1},x_{n,n+1})\in\mathbb{Z}^{2n}\). Iterating Lemma 17, we obtain \[\label{eq:32summing-over-x-prime} \sum_{\mathbf{x}'}C(\mathbf{x})\leq 2^{8\alpha n}(1+2\zeta(\alpha))^{2n}C(x_{0,1},x_{n+1,n}).\tag{113}\] Choosing \(\beta\geq 128\alpha c_0+(32\alpha/c_2)\log[2^7M(1+2\zeta(\alpha))]\), and combining Equations 92 , 112 and 113 , we obtain \[\overline{\mathrm{E}}_{\beta,J}(\varGamma_1,\varGamma_2)\leq e^{-nc_2\beta/16}\tilde{\mathrm{E}}_\alpha(\varGamma_1,\varGamma_2).\]Since \(\lvert \mathrm{P}([n])\rvert=2^n\), if \(\beta\geq (16/c_2)\log 2\), then \[\label{eq:32summing-over-n-E-bar-n} \sum_{n=1}^\infty\overline{\mathrm{E}}_{\beta,n}(\varGamma_1,\varGamma_2)\leq \tilde{\mathrm{E}}_\alpha(\varGamma_1,\varGamma_2).\tag{114}\] Combining Equations 92 and@eq:eq:32summing-over-n-E-bar-n , we obtain Equation 107 . The lemma has thus been proved. ◻

7.2 One-polymer bounds↩︎

Consider a finite, totally ordered set \(A\) with cardinality \(N\mathscr{}in\mathbb{N}\). Given \(P\in\mathbbold{P}(A)\), it is useful to define \(T_P:\mathbbold{T}(P)\rightarrow\mathbbold{T}(A)\) in the following way. First, order each element of \(Q\in P\); that is, \(Q=\{a_{Q,1},...,a_{Q,\lvert Q\rvert}\}\), where \(i<j\) implies \(a_{Q,i}<a_{Q,j}\). Given \(T\in\mathbbold{T}(P)\), define \[E(T_P(T))\mathrel{\vcenter{:}}=\{\{a_{Q,1},a_{R,1}\}:\{Q,R\}\in E(T)\}\sqcup\bigsqcup_{S\in P}\{\{a_{S,i},a_{S,{i+1}}\}:i\in[\lvert S\rvert-1]\}.\]For a graphical description, see Figure 17.

Figure 17: Example of T_P, where A=[8] and P=\{\{1,5\},\{2,3,6\},\{4\},\{7,8\}\}.

We are now in a position to tackle the case of a single set partition in Equation 98 .

Lemma 22. There exists \(\beta_0>0\) such that if \(\beta\geq\beta_0\) and \(A\in\mathrm{P}(\mathbb{Z})\), then \[\label{eq:32A-in-a-single-contour} \sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{A\subset\mathrm{I}_-(\gamma)}e^{-\beta H(\gamma)}\leq [2(1+\mathrm{diam}(A))]^{-\frac{c_0\beta}{2}}.\tag{115}\]

Proof. If \(A\subset\mathrm{I}_-(\gamma)\), then \(\mathrm{diam}(A)+1\leq \mathrm{diam}(\gamma)\). Therefore, recalling Hypothesis 2, we have\[\sum_{\gamma\in\mathscr{C}}\mathbbold{1}_{A\subset\mathrm{I}_-(\gamma)}e^{-\beta H(\gamma)}\leq \sum_{n=\lceil 1+\log_2(1+\mathrm{diam}(A))\rceil}2^{c_1n}2^{-c_0\beta n}\leq [2(1+\mathrm{diam}(A))]^{-\frac{c_0\beta}{2}},\]where the second inequality holds if \(\beta\geq 2(1+c_1)\). The lemma has thus been proved. ◻

Lemma 23. Let \(A\in\mathrm{P}(\mathbb{Z})\), then there exists \(\beta_0=\beta_0(\lvert A\rvert)>0\) such that \[\label{eq:32trivial-partition-correlation-bound} \overline{\mathrm{W}}_\beta(\{A\})\leq2^{-\frac{c_0\beta}{2^{\lvert A\rvert+5}}}\sum_{T\in\mathbbold{T}(A)}\prod_{\{a_1,a_2\}\in E(T)}\lvert a_1-a_2\rvert^{-\alpha}.\tag{116}\]

Proof. By Corollary 7, we have \[\overline{\mathrm{W}}_\beta(\{A\})\leq\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{A\subset\mathrm{I}_-(\Gamma)}\overline{\mathrm{W}}_\beta(\Gamma)\leq \sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{A\subset\mathrm{I}_-(\Gamma)}\bar{z}_{\beta/2}(\Gamma).\] By the definition of polymer, we may write \[\sum_{\Gamma\in\mathscr{P}}\mathbbold{1}_{A\subset\mathrm{I}_-(\Gamma)}\bar{z}_{\beta/2}(\Gamma)=\sum_{P\in\mathbbold{P}(A)}\sum_{\Upsilon}\tilde{\mathbbold{1}}_P(\{\gamma_a\})\sum_{\Gamma\in\mathscr{P}}\bar{z}_{\beta/2}(\Gamma)\prod_{a\in A}\mathbbold{1}_{a\in\mathrm{I}_-(\gamma_a)}\mathbbold{1}_{\gamma_a\in\Gamma},\]where \(\Upsilon=(\gamma_a:a\in A)\in\mathscr{C}^A\) and \[\tilde{\mathbbold{1}}_P(\Upsilon)\mathrel{\vcenter{:}}=\prod_{Q\in P}\prod_{q_1,q_2\in Q}\mathbbold{1}_{\gamma_{q_1}=\gamma_{q_2}}\prod_{\substack{R,S\in P\\R\neq S}}\prod_{\substack{r\in R\\s\in S}}\mathbbold{1}_{\gamma_r \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_s}.\]Note that if \(\tilde{\mathbbold{1}}_P(\Upsilon)=1\), then there are just \(\lvert P\rvert\) diverse contours. Thus, we write\[\label{eq:32Upsilon-partition} \Upsilon_P=(\gamma_Q:Q\in P)\in\mathscr{C}_P\mathrel{\vcenter{:}}=\{\Upsilon\in\mathscr{C}^P:(Q\neq R\in P)\Rightarrow\gamma_Q \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_R\}.\tag{117}\]

By Corollaries 3, 5 and 6, we obtain \[\overline{\mathrm{W}}_\beta(\{A\})\leq\sum_{P\in\mathbbold{P}(A)}\sum_{T\in\mathbbold{T}(P)}\sum_{\Upsilon_P}\prod_{Q\in P}\mathbbold{1}_{Q\subset\mathrm{I}_-(\gamma_Q)}e^{-\frac{\beta H(\gamma_Q)}{2^{\lvert P\rvert+2}}}\prod_{\{R,S\}\in E(T)}\Phi(\gamma_R,\gamma_S).\] Consider \(\{R,S\}\in E(T)\). Since \(R\subset \mathrm{I}_-(\gamma_R)\) and \(S\subset\mathrm{I}_-(\gamma_S)\), there exist \(x_S\in\mathrm{I}_-(\gamma_S)\) and \(x_R\in\mathrm{I}_-(\gamma_R)\) such that \(\lvert x_R-x_S\rvert=\mathrm{diam}(R\sqcup S)\). By Lemma 20, we obtain \[\Phi(\gamma_R,\gamma_S)\leq 2^{2\alpha}\mathrm{diam}(\gamma_1)^{1+\alpha}\mathrm{diam}(\gamma_2)^{1+\alpha}\mathrm{diam}(R\sqcup S)^{-\alpha}.\]Fix \(Q\in P\). By Hypothesis 2 and Lemma 22, if \(\beta\geq 2^{\lvert A\rvert+2}c_0(1+\alpha)\lvert A\rvert\), then \[\sum_{\gamma_Q\in\mathscr{C}}\mathbbold{1}_{Q\subset\mathrm{I}_-(\gamma_i)}\mathrm{diam}(\gamma_Q)^{(1+\alpha)\delta_T(Q)}e^{-\frac{\beta H(\gamma_Q)}{2^{\lvert P\rvert+2}}}\leq [2(1+\mathrm{diam}(Q))]^{-\frac{c_0\beta}{2^{\lvert P\rvert+4}}}.\] Combining all these estimates, we obtain \[\label{eq:32W-bar-A-in-term-of-trees} \overline{\mathrm{W}}_\beta(\{A\})\leq 2^{3\alpha\lvert A\rvert}\sum_{P\in\mathbbold{P}(A)}\sum_{T\in\mathbbold{T}(P)}\prod_{Q\in P}[2(1+\mathrm{diam}(Q))]^{-\frac{c_0\beta}{2^{\lvert P\rvert+4}}}\prod_{\{R,S\}\in E(T)}\mathrm{diam}(R\sqcup S)^{-\alpha}.\tag{118}\]

We will conclude by employing the canonical application \(T_P:\mathbbold{T}(P)\rightarrow\mathbbold{T}(A)\) introduced at the beginning of this subsection. Endow \(A\) with the usual ordering of \(\mathbb{Z}\). Let \(P\in\mathbbold{P}(A)\) and \(Q\in P\), then \[\prod_{q=1}^{\lvert Q\rvert-1}\lvert a_{Q,q}-a_{Q,q+1}\rvert\leq \left(1+\sum_{q=1}^{\lvert Q\rvert-1}\lvert a_{Q,q}-a_{Q,q+1}\rvert\right)^{\lvert Q\rvert}\]Hence, if \(\beta\geq 2^{\lvert A\rvert+4}\lvert A\rvert \alpha c_0^{-1}\), then\[(1+\mathrm{diam}(Q))^{-\frac{c_0\beta}{2^{\lvert P\vert+4}}}\leq \prod_{q=1}^{\lvert Q\rvert-1}\lvert a_{Q,q}-a_{Q,q+1}\rvert^{-\alpha}.\]Furthermore, since \(\lvert a_{R,1}-a_{S,1}\rvert\leq\mathrm{diam}(R\sqcup S)\), for every \(P\in\mathbbold{P}(A)\), we have \[\prod_{Q\in P}(1+\mathrm{diam}(Q))^{-\frac{c_0\beta}{2^{\lvert P\rvert+4}}}\prod_{\{R,S\}\in E(T)}\mathrm{diam}(R\sqcup S)^{-\alpha}\leq\prod_{\{a_1,a_2\}\in E(T_P)}\lvert a_1-a_2\rvert^{-\alpha}.\] Finally, since the application \(T_P:\mathbbold{T}(P)\rightarrow\mathbbold{T}(A)\) is injective, we obtain \[\overline{\mathrm{W}}_\beta(\{A\})\leq\sum_{T\in\mathbbold{T}(A)}\left[\sum_{P\in\mathbbold{P}(A)}2^{3\alpha\lvert A\rvert-2^{-\lvert A\rvert-4}c_0\beta}\right]\prod_{\{a_1,a_2\}\in E(T)}\lvert a_1-a_2\rvert^{-\alpha}.\]Choosing \(\beta\geq 2^{\lvert A\rvert+5}\left[3\alpha\lvert A\rvert+\log_2\lvert\mathbbold{P}(A)\rvert\right]c_0^{-1}\), we obtain Equation 116 . The lemma has thus been proved. ◻

7.3 Two-point correlations↩︎

In this subsection, we employ Lemmas 21 and 23 to obtain an upper bound to the two-point correlation function.

Proof of Theorem 2. By Equation 99 , we have \[\lvert\langle\sigma_x;\sigma_y\rangle^+_{\Lambda,\beta}\rvert\leq\overline{\mathrm{W}}_\beta(\{x,y\})+\overline{\mathrm{W}}_\beta(\{x\},\{y\}),\]where we write \(\overline{\mathrm{W}}_\beta(A_1,...,A_n)\) instead of \(\overline{\mathrm{W}}_\beta(\{A_1,...,A_n\})\).

We have already dealt with \(\overline{\mathrm{W}}_\beta(\{x,y\})\) in Lemma 23. Let us bound the term \(\overline{\mathrm{W}}_\beta(\{x\},\{y\})\). By Equation 100 , we have \[\overline{\mathrm{W}}_\beta(\{x\},\{y\})=\sum_{\Gamma_1,\Gamma_2\in\mathscr{P}}\mathbbold{1}_{x\in\mathrm{I}_-(\Gamma_1)}\mathbbold{1}_{y\in\mathrm{I}_-(\Gamma_2)}\overline{\mathrm{W}}_\beta(\Gamma_1)\overline{\mathrm{E}}_\beta(\Gamma_1,\Gamma_2)\overline{\mathrm{W}}_\beta(\Gamma_2).\]Let us write \[\begin{align} &\mathbbold{1}(\mathbf{\Gamma},\Upsilon,\mathbf{x})\mathrel{\vcenter{:}}=\mathbbold{1}_{x\in\mathrm{I}_-(\Gamma_1)}\mathbbold{1}_{y\in\mathrm{I}_-(\Gamma_2)}\mathbbold{1}_{\gamma_1\in\Gamma_1}\mathbbold{1}_{\gamma_2\in\Gamma_2}\mathbbold{1}_{x_1\in\mathrm{I}_-(\gamma_1)}\mathbbold{1}_{x_2\in\gamma_2},\\ &\epsilon_\alpha(\Upsilon)\mathrel{\vcenter{:}}= e^{\alpha c_0 H(\gamma_1)}e^{\alpha c_0 H(\gamma_2)},\\ &\xi_\beta(\mathbf{\Gamma})\mathrel{\vcenter{:}}= \bar{z}_{\beta/2}(\Gamma_1)\bar{z}_{\beta/2}(\Gamma_2), \end{align}\]where \(\mathbf{\Gamma}=(\Gamma_1,\Gamma_2)\in\mathscr{P}^2\), \(\Upsilon=(\gamma_1,\gamma_2)\in\mathscr{C}^2\) and \(\mathbf{x}=(x_1,x_2)\in\mathbb{Z}^2\).

By Corollary 7 and Lemma 21, we obtain \[\overline{\mathrm{W}}_\beta(\{x\},\{y\})\leq (8M)^\alpha\sum_{\mathbf{\Gamma}}\sum_{\Upsilon}\sum_{\mathbf{x}}\mathbbold{1}(\mathbf{\Gamma},\Upsilon,\mathbf{x})\xi_\beta(\mathbf{\Gamma})\epsilon_\alpha(\Upsilon)C(x_1,x_2).\]Let us sum over the polymers. Writing \(\Upsilon'=(\gamma_1',\gamma_2')\in\mathscr{C}^2\), recalling Equation 103 , and introducing a sum over \(J\subset [2]\) (as in the proof of Lemma 21), we have \[\mathbbold{1}(\mathbf{\Gamma},\Upsilon,\mathbf{x})=\sum_J\sum_{\Upsilon'}\mathbbold{1}_J(\Upsilon,\Upsilon',\mathbf{x})\mathbbold{1}_{\gamma_1,\gamma'_1\in\Gamma_1}\mathbbold{1}_{\gamma_2,\gamma'_2\in\Gamma_2},\]where \[\mathbbold{1}_J(\Upsilon,\Upsilon',\mathbf{x})\mathrel{\vcenter{:}}=\mathbbold{1}_{x\in\mathrm{I}_-(\gamma'_1)}\mathbbold{1}_{y\in\mathrm{I}_-(\gamma'_2)}\mathbbold{1}_{x_1\in\mathrm{I}_-(\gamma_1)}\mathbbold{1}_{x_2\in\gamma_2}\prod_{j\in J}\mathbbold{1}_{\gamma_j \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma'_j}\prod_{j\in[2]\setminus J}\mathbbold{1}_{\gamma_j=\gamma'_j}.\]Thus, by Corollaries 3 and 5, we obtain \[\label{eq:32two-point-correlation-summed-over-polymers} \overline{\mathrm{W}}_\beta(\{x\},\{y\})\leq (8M)^\alpha\sum_J\sum_{\Upsilon,\Upsilon'}\sum_{\mathbf{x}}\mathbbold{1}_J(\Upsilon,\Upsilon',\mathbf{x})\epsilon_\alpha(\Upsilon)\epsilon_\beta(\Upsilon,\Upsilon')C(x_1,x_2)\prod_{j\in J}\Phi(\gamma_j,\gamma'_j),\tag{119}\] where \[\epsilon_\beta(\Upsilon,\Upsilon')\mathrel{\vcenter{:}}= e^{-\frac{1}{8}\beta H(\gamma_1)}e^{-\frac{1}{8}\beta H(\gamma'_1)}e^{-\frac{1}{8}\beta H(\gamma_2)}e^{-\frac{1}{8}\beta H(\gamma'_2)}.\]By Lemma 20, we have \[\label{eq:32two-point-correlation-contours-to-points} \mathbbold{1}_{J}(\Upsilon,\Upsilon',\mathbf{x})C(x_1,x_2)\leq 16^\alpha \epsilon'_\alpha(\Upsilon,\Upsilon')C(x,x_1)C(x_1,x_2)C(x_2,y)\tag{120}\] where \[\epsilon'_\alpha(\Upsilon,\Upsilon')\mathrel{\vcenter{:}}= e^{2\alpha c_0H(\gamma_1)}e^{2\alpha c_0H(\gamma'_1)}e^{2\alpha c_0H(\gamma_2)}e^{2\alpha c_0H(\gamma'_2)}.\]Combining Hypotheses 2 and 3 with Equations 119 and 120 , and choosing \(\beta\geq 48\alpha c_0\), we obtain \[\overline{\mathrm{W}}_\beta(\{x\},\{y\})\leq (2^7M)^\alpha e^{-\frac{1}{4}c_2\beta}\sum_J\sum_{\mathbf{x}}C(x,x_1)C(x_1,x_2)C(x_2,y).\] Finally, by Lemma 17 and \(\lvert \mathrm{P}([2])\rvert=4<2^{2\alpha}\), we obtain \[\label{eq:32two-point-correlation-W-bar-bound-x-y} \overline{\mathrm{W}}_\beta(\{x\},\{y\})\leq (2^{17}M)^\alpha(1+2\zeta(\alpha))^2e^{-\frac{1}{4}c_2\beta}C(x,y).\tag{121}\]

Choosing \(\beta\geq(8\alpha/c_2)\log[2^{17}M(1+2\zeta(\alpha))]\) and \(c_3=\min\{2^{-7}(\log2)c_0,2^{-3}c_2\}\), by Lemma 23 (applied to \(A=\{x,y\}\)) and Equation 121 , we obtain Equation 1 . The theorem has thus been proved. ◻

7.4 Many-point correlations↩︎

We close this work by tackling the matter of many-point correlation functions. The starting point is applying Corollary 7 and Lemma 21 to Equation 100 . Thus, if \(A\in\mathrm{P}(\mathbb{Z})\) and \(P\in\mathbbold{P}(A)\) ordered canonically (see Equation 98 ) such that \(\lvert A\rvert=N\geq 3\) and \(\lvert P\rvert=m\geq 2\), then \[\overline{\mathrm{W}}_\beta(P)\leq\sum_{n=m}^{2m-2}\sum_{T\in\mathbbold{T}^*_{n,m}}\overline{\mathrm{W}}'_{\beta}(P,T)\mathrel{\vcenter{:}}=\sum_{n=m}^{2m-2}\sum_{T\in\mathbbold{T}^*_{n,m}}\sum_{\substack{\mathbf{\Gamma}\in\mathscr{P}^n}}\prod_{i=1}^m\mathbbold{1}_{A_i\subset\mathrm{I}_-(\Gamma_i)}\prod_{j=1}^n\overline{z}_{\beta/2}(\Gamma_j)\prod_{\{p,q\}\in E(T)}\tilde{\mathrm{E}}_\alpha(\Gamma_p,\Gamma_q).\]

Lemma 24. Let \(A\in\mathrm{P}(\mathbb{Z})\), \(P\in\mathbbold{P}(A)\) and \(T\in\mathbbold{T}^*_{n,m}\) such that \(\lvert A\rvert=N\geq 3\), \(\lvert P\rvert=m\geq 2\) and \(m\leq n\leq 2m-2\). Then \[\label{eq:32overline-W-P-T-bound} \overline{\mathrm{W}}'_\beta(P,T)\leq F(N)e^{-\frac{\beta c_2}{2^{3N+3}}}\sum_{T'\in\mathbbold{T}(A)}\prod_{\{a_1,a_2\}\in E(T')}\lvert a_1-a_2\rvert^{-\alpha},\tag{122}\] where \[F(N)=2^{29\alpha N^2}(1+2\zeta(\alpha))^{3N}\left[\lvert\mathbbold{T}_{3N}\rvert\right]^{3N}\lvert \mathbbold{P}_{3N}\rvert.\]

Proof. Recalling Notation 1, we write \(\Upsilon=(\Upsilon_1,...,\Upsilon_n)\) and \(\mathbf{x}=(\mathbf{x}_1,...,\mathbf{x}_n)\), where \(\Upsilon_i\in\mathscr{C}^{E_T(i)}\) and \(\mathbf{x}_i\in\mathbb{Z}^{E_T(i)}\). Then \[\prod_{\{p,q\}\in E(T)}\tilde{\mathrm{E}}_\alpha(\Gamma_p,\Gamma_q)=\sum_{\Upsilon}\sum_{\mathbf{x}}\epsilon_\alpha(\Upsilon)C(\mathbf{x})\mathbbold{1}(\mathbf{\Gamma},\Upsilon)\mathbbold{1}(\Upsilon,\mathbf{x}),\]where \[\begin{align} &\epsilon_\alpha(\Upsilon)\mathrel{\vcenter{:}}=\prod_{i=1}^n\prod_{e\in E_T(i)}e^{\alpha c_0 H(\gamma_{i,e})},&\mathbbold{1}(\mathbf{\Gamma},\Upsilon)\mathrel{\vcenter{:}}=\prod_{i=1}^n\prod_{e\in E_T(i)}\mathbbold{1}_{\gamma_{i,e}\in\Gamma_i}\eqqcolon\prod_{i=1}^n\mathbbold{1}(\Gamma_i,\Upsilon_i),\\ &C(\mathbf{x})\mathrel{\vcenter{:}}=\prod_{\{p,q\}=e\in E(T)}C(x_{p,e},x_{q,e}), &\mathbbold{1}(\Upsilon,\mathbf{x})\mathrel{\vcenter{:}}=\prod_{i=1}^n\prod_{e\in E_T(i)}\mathbbold{1}_{x_{i,e}\in\mathrm{I}_-(\gamma_{i,e})}\eqqcolon\prod_{i=1}^n\mathbbold{1}(\Upsilon_i,\mathbf{x}_i). \end{align}\]On the other hand, we have \[\prod_{i=1}^m\mathbbold{1}_{A_i\in\mathrm{I}_-(\Gamma_i)}=\sum_{\Upsilon'}\sum_{\mathbf{x}'}\mathbbold{1}(\mathbf{x'},\mathbf{A})\mathbbold{1}(\Upsilon',\mathbf{x'})\mathbbold{1}(\mathbf{\Gamma},\Upsilon'),\]where \(\Upsilon'=(\Upsilon'_1,...,\Upsilon'_m)\) and \(\Upsilon'_i=(\gamma_{i,a_i}:a_i\in A_i)\in\mathscr{C}^{A_i}\), \(\mathbf{x}'=(\mathbf{x}'_1,...,\mathbf{x}'_m)\) and \(\mathbf{x}'_i=(x_{i,a_i}:a_i\in A_i)\in\mathbb{Z}^{A_i}\), and \[\begin{align} &\mathbbold{1}(\mathbf{x}',\mathbf{A})\mathrel{\vcenter{:}}= \prod_{i=1}^m\prod_{a_i\in A_i}\mathbbold{1}_{x_{i,a_i}=a_i}\eqqcolon\prod_{i=1}^m\mathbbold{1}(\mathbf{x}'_i,A_i),\\ &\mathbbold{1}(\Upsilon',\mathbf{x}')\mathrel{\vcenter{:}}=\prod_{i=1}^m\prod_{a_i\in A_i}\mathbbold{1}_{x_{i,a_i}\in\mathrm{I}_-(\gamma_{i,a_i})}\eqqcolon\prod_{i=1}^m\mathbbold{1}(\Upsilon'_i,\mathbf{x}'_i),\\ &\mathbbold{1}(\mathbf{\Gamma},\Upsilon')\mathrel{\vcenter{:}}=\prod_{i=1}^m\prod_{a_i\in A_i}\mathbbold{1}_{\gamma_{i,a_i}\in\Gamma_i}. \end{align}\]

We must now write \(\mathbbold{1}(\Gamma_i,\Upsilon_i)\mathbbold{1}(\Gamma_i,\Upsilon'_i)\) in a way that allows us to transform polymers into contours. There are two cases to consider: \(i\in [m]\) and \(m<i\leq n\). We begin with \(i\in[m]\). Let \[\mathbbold{1}(\Gamma_i,\Upsilon_i)\mathbbold{1}(\Gamma_i,\Upsilon'_i)=\mathbbold{1}(\Gamma_i,\Upsilon_i,\Upsilon'_i)=\prod_{e\in E_T(i)}\mathbbold{1}_{\gamma_{i,e}\in\Gamma_i}\prod_{a_i\in A_i}\mathbbold{1}_{\gamma_{i,a_i}\in\Gamma_i}.\]Introducing partitions \(P_i\in\mathbbold{P}(A_i\sqcup E_T(i))\), we have \[\mathbbold{1}(\Gamma_i,\Upsilon_i,\Upsilon'_i)=\sum_{P_i}\mathbbold{1}_{P_i}(\Upsilon_i,\Upsilon'_i)\mathbbold{1}(\Gamma_i,\Upsilon_i,\Upsilon'_i),\]where\[\mathbbold{1}_{P_i}(\Upsilon_i,\Upsilon'_i)\mathrel{\vcenter{:}}=\prod_{Q_i\in P_i}\prod_{q_i,q'_i\in Q_i}\mathbbold{1}_{\gamma_{i,q_i}=\gamma_{i,q'_i}}\prod_{\substack{R_i,S_i\in P_i\\R_i\neq S_i}}\prod_{\substack{r_i\in R_i\\s_i\in S_i}}\mathbbold{1}_{\gamma_{i,r_i} \mathrel{ \ooalign{ \hss\raisebox{-1pt}{\sim}\hss\cr \hss\raisebox{3.5pt}{\scalebox{0.4}{\triangle}}\hss\cr\\ } }\gamma_{i,s_i}}.\] Note that if \(\mathbbold{1}_{P_i}(\Upsilon_i,\Upsilon'_i)=1\), then there are \(\lvert P_i\rvert\) distinct contours. In this case, we have \[\mathbbold{1}_{P_i}(\Upsilon_i,\Upsilon'_i)\mathbbold{1}(\Gamma_i,\Upsilon_i,\Upsilon'_i)=\sum_{\Upsilon_{i,P_i}}\mathbbold{1}(\Upsilon_i,\Upsilon'_i,\Upsilon_{i,P_i})\mathbbold{1}(\Gamma_i,\Upsilon_{i,P_i}),\]where \(\Upsilon_{i,P_i}\in\mathscr{C}_{P_i}\) (see Equation 117 ), and \[\begin{align} &\mathbbold{1}(\Upsilon_i,\Upsilon'_i,\Upsilon_{i,P_i})=\prod_{Q_i\in P_i}\prod_{e\in E_T(i)}\mathbbold{1}_{\gamma_{i,e}=\gamma_{i,Q_i}}\prod_{a_i\in Q_i\cap A_i}\mathbbold{1}_{\gamma_{i,a_i}=\gamma_{i,Q_i}},\\&\mathbbold{1}(\Gamma_i,\Upsilon_{i,P_i})=\prod_{Q_i\in P_i}\mathbbold{1}_{\gamma_{i,Q_i}\in\Gamma_i}. \end{align}\] Thus, by Equation 79 and Corollaries 3, 5 and 6, we obtain \[\sum_{\Gamma_i}\bar{z}_{\beta/2}(\Gamma_i)\mathbbold{1}(\Gamma_i,\Upsilon_{i,P_i})\leq\sum_{T_i\in\mathbbold{T}(P_i)}\prod_{Q_i\in P_i}e^{-\frac{\beta H(\gamma_{i,Q_i})}{2^{\lvert P_i\rvert+2}}}\prod_{\{R_i,S_i\}\in E(T_i)}\Phi(\gamma_{i,R_i},\gamma_{i,S_i}).\]The case \(m<i\leq n\) is a simplification: we may formally consider \(A_i\equiv \emptyset\) and proceed exactly as we have done for the case \(i\in[m]\).

We must now employ Lemma 20 to sum over contours. More precisely, we must bound \[\label{eq:32many-point-correlations-sum-to-eliminate-contours} F(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T})\mathrel{\vcenter{:}}= C(\mathbf{x})\sum_{\Upsilon}\sum_{\Upsilon'}\sum_{\Upsilon_{\mathbf{P}}}\mathbbold{1}(\Upsilon',\mathbf{x}')\mathbbold{1}(\Upsilon,\mathbf{x})\mathbbold{1}(\Upsilon,\Upsilon',\Upsilon_{\mathbf{P}})\epsilon_\beta(\Upsilon_\mathbf{P})\varPhi(\Upsilon_\mathbf{P},\mathbf{T}),\tag{123}\] where \(\Upsilon_{\mathbf{P}}=(\Upsilon_{1,P_1},...,\Upsilon_{n,P_n})\), \(\mathbf{T}=(T_1,...,T_n)\) and \[\begin{align} \mathbbold{1}(\Upsilon,\Upsilon',\Upsilon_{\mathbf{P}})&\mathrel{\vcenter{:}}=\prod_{i=1}^n\mathbbold{1}(\Upsilon_i,\Upsilon'_i,\Upsilon_{i,P_i}),\\ \epsilon_\beta(\Upsilon_\mathbf{P})&\mathrel{\vcenter{:}}=\prod_{i=1}^n\prod_{Q_i\in P_i}e^{-\frac{\beta H(\gamma_{i,Q_i})}{2^{\lvert P_i\rvert+2}}},\\ \varPhi(\Upsilon_{\mathbf{P}},\mathbf{T})&\mathrel{\vcenter{:}}=\prod_{i=1}^n\prod_{\{R_i,S_i\}\in E(T_i)}\Phi(\gamma_{i,R_i},\gamma_{i,S_i}). \end{align}\]

First, for each \(i\in[n]\) and each \(P_i\in\mathbbold{P}(A_i\sqcup E_T(i))\), we order every \(Q_i\in P_i\) canonically as \(Q_i=\{q_{Q_i,1},...,q_{Q_{i},\lvert Q_i\rvert}\}\) such that

  1. if \(\lvert Q_i\cap E_T(i)\rvert\mathrel{\vcenter{:}}= n(Q_i,T)>0\), then \(Q_i\cap E_T(i)=\{y_{Q_i,q_i}:q_i\in[n(Q_i,T)]\}\). Furthermore, the edges are ordered in ascending order: if \(y_{Q_i,q_i}=\{i,j_{Q_i,q_i}\}\) and \(y_{Q_i,q'_i}=\{i,j_{Q_i,q'_i}\}\), then \(q_i<q'_i\leq q_i(T)\) if and only if \(j_{Q_i,q_i}<j_{Q_i,q'_i}\);

  2. the set \(\{y_{Q_i,q_i}:n(Q_i,T)<q_i\leq \lvert Q_i\rvert\}\subset A_i\) is also ordered in ascending order: if \(n(Q_i,T)<q_i<q'_i\), then \(y_{Q_i,q_i}<y_{Q_i,q'_i}\).

For reasons that shall be clear shortly, arguing as in Lemma 20, we obtain the following "phantom" bound\[\label{eq:32phantom-one-inside-partition-bound} 1\leq e^{(\lvert Q_i\rvert-1)\alpha c_0H(\gamma_{i,Q_i})}\prod_{q_i=1}^{\lvert Q_i\rvert-1}C\left(x_{i,y_{Q_i,q_i}},x_{i,y_{Q_i,q_i+1}}\right).\tag{124}\] Let us now bound \(\varPhi(\Upsilon_\mathbf{P},\mathbf{T})\) employing the canonical ordering described above.

Let \(\{R_i,S_i\}\in E(T_i)\), we write \(R_i=\{y_{R_i,r_i}:r_i\in [\lvert R_i\rvert]\}\) and \(S_i=\{y_{S_i,s_i}:s_i\in[\lvert S_i]\}\), in accordance with the canonical ordering. By Lemma 20, we have \[\label{eq:32T-i-edge-bounds} \Phi(\gamma_{i,R_i},\gamma_{i,S_i})\leq 4^\alpha e^{2\alpha c_0H(\gamma_{i,R_i})}e^{2\alpha c_0H(\gamma_{i,S_i})}C\left(x_{i,y_{R_i,1}},x_{i, y_{S_i,1}}\right).\tag{125}\] Since \(\lvert Q_i\rvert+\delta_{T_i}(Q_i)\leq 2N\), combining Equation 124 and 125 , we have \[\label{eq:32C-varPhi-bound} C(\mathbf{x})\varPhi(\Upsilon_{\mathbf{P}},\mathbf{T})\leq 4^{N\alpha}\bar\epsilon_{\alpha}(\Upsilon_{\mathbf{P}})C(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T}),\tag{126}\] where \[\begin{align} &\bar\epsilon_{\alpha}(\Upsilon_{\mathbf{P}})\mathrel{\vcenter{:}}=\prod_{i=1}^n\prod_{Q_i\in P_i}e^{4N\alpha c_0H(\gamma_{i,Q_i})},\\ &\frac{C(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T})}{C(\mathbf{x})}=\prod_{i=1}^n\prod_{Q_i\in P_i}\prod_{q_i=1}^{\lvert Q_i\rvert-1}C(x_{i,y_{Q_i,q_i}},x_{i,y_{Q_i,q_i+1}})\prod_{\{R_i,S_i\}\in E(T)}C(x_{i,y_{R_i,1}},x_{i,y_{S_i,1}}). \end{align}\]This might seem very complicated, but concretely we have \[C(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T})=\prod_{\{t,u\}\in E(T(\mathbf{P},\mathbf{T}))}C(x_t,x_u),\]where \(T(\mathbf{P},\mathbf{T})\) is the tree with vertex set \[V(T(\mathbf{P},\mathbf{T}))= \bigsqcup_{i=1}^n\{i\}\times (E_T(i)\sqcup A_i)\] and with the following set of edges\[\begin{align} E(T(\mathbf{P},\mathbf{T}))=&\left\{\{(i,e),(j,e)\}:\{i,j\}=e\in E(T)\right\}\\ &\sqcup\{\{(i,y_{Q_i,q_i}),(i,y_{Q_i,q_i+1})\}:i\in [n], Q_i\in P_i, q_i\in[\lvert Q_i\rvert-1]\}\\ &\sqcup\{\{(i,y_{R_i,1}),(i,y_{S_i,1})\}:i\in [n], \{R_i,S_i\}\in E(T_i)\}. \end{align}\]Note that \[\label{eq:32cardinality-special-tree-many-point-correlation} \lvert V(T(\mathbf{P},\mathbf{T}))\rvert=2(n-1)+N\tag{127}\] and\[\label{eq:32leaves-special-tree-many-point-correlation} \mathrm{L}(T(\mathbf{P},\mathbf{T}))\subset\sqcup_{i=1}^m\{i\}\times A_i\tag{128}\]

We are now in a position to sum over contours and obtain a decent upper bound for Equation 123 . By Equation 126 , we have \[\frac{F(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T})}{C(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T})}\leq 4^{N\alpha}\sum_{\Upsilon}\sum_{\Upsilon'}\sum_{\Upsilon_{\mathbf{P}}}\mathbbold{1}(\Upsilon',\mathbf{x}')\mathbbold{1}(\Upsilon,\mathbf{x})\mathbbold{1}(\Upsilon,\Upsilon',\Upsilon_{\mathbf{P}})\epsilon_\beta(\Upsilon_\mathbf{P})\bar\epsilon_\alpha(\Upsilon_\mathbf{P})\eqqcolon 4^{N\alpha}F'(\mathbf{P},\mathbf{T}).\]In fact, \(F'(\mathbf{P},\mathbf{T})\) is rather easy to bound. First, note that \[\mathbbold{1}(\Upsilon',\mathbf{x}')\mathbbold{1}(\Upsilon,\mathbf{x})\mathbbold{1}(\Upsilon,\Upsilon',\Upsilon_\mathbf{P})=\mathbbold{1}(\Upsilon,\Upsilon',\Upsilon_\mathbf{P})\prod_{i=1}^n\prod_{Q_i\in P_i}\prod_{q_i=1}^{\lvert Q_i\rvert}\mathbbold{1}_{x_{i,y_{Q_i,q_i}}\in \mathrm{I}_-(\gamma_{i,Q_i})}.\]Thus, \[F'(\mathbf{P},\mathbf{T})\leq \prod_{i=1}^n\prod_{Q_i\in P_i}\sum_{\gamma_{i,Q_i}}\mathbbold{1}_{x_{i,y_{Q_i,1}}\in\mathrm{I}_-(\gamma_{i,Q_i})}e^{2N\alpha c_0H(\gamma_{i,Q_i})}e^{-\frac{\beta H(\gamma_{i,Q_i})}{2^{\lvert P_i\rvert+2}}}.\]Introducing equivalence classes of contours under integer translation (note that \(\lvert\mathrm{I}_-(\gamma_{i,Q_i})\rvert\) elements of \([\gamma_{i,Q_i}]\) satisfy \(x_{i,y_{Q_i,1}}\in\mathrm{I}_-(\gamma_{i,Q_i})\)), recalling Hypotheses 2 and 3, and, finally, choosing \(\beta\geq 2^{N+5}N\alpha c_0\), we obtain\[F'(\mathbf{P},\mathbf{T})\leq \prod_{i=1}^n\prod_{Q_i\in P_i}e^{-\frac{\beta c_2}{2^{\lvert P_i\rvert+3}}}=\prod_{i=1}^ne^{-\frac{\lvert P_i\rvert\beta c_2}{2^{\lvert P_i\rvert+3}}}.\]Having eliminated the sums over contours, we have \[\label{eq:32overline-W-P-T-intermediary-bound} \overline{\mathrm{W}}_\beta(P,T)\leq 4^{N\alpha}\sum_{\mathbf{P}}\prod_{i=1}^ne^{-\frac{\lvert P_i\rvert\beta c_2}{2^{\lvert P_i\rvert+3}}}\sum_{\mathbf{T}}\sum_{\mathbf{x}'}\sum_{\mathbf{x}}\mathbbold{1}(\mathbf{x}',\mathbf{A})C(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T}).\tag{129}\] The final step is simply to apply Lemma 18, which can be done because of Equation 128 . Since \(n\leq 2m-2\leq 2N\), recalling Equation 127 , we obtain \[\label{eq:32C-x-x-prime-P-T-bound} \sum_{\mathbf{x}}C(\mathbf{x},\mathbf{x}',\mathbf{P},\mathbf{T})\leq \left[\lvert\mathbbold{T}_{3N}\rvert\right]^{2N}2^{27\alpha N^2}(1+2\zeta(\alpha))^{3N}\sum_{T'\in\mathbbold{T}(A_P)}\prod_{\{t,u\}\in E(T')}C(x_t,x_u),\tag{130}\] where \(A_P\mathrel{\vcenter{:}}=\sqcup_{i=1}^m\{i\}\times A_i\). On the other hand, \[\label{eq:32crudest-bound-over-partitions-trees} \sum_{\mathbf{P}}\sum_{\mathbf{T}}\prod_{i=1}^ne^{-\frac{\lvert P_i\rvert\beta c_2}{2^{\lvert P_i\rvert+3}}}\leq \left[\lvert \mathbbold{T}_{3N}\rvert\right]^N\lvert \mathbbold{P}_{3N}\rvert e^{-\frac{\beta c_2}{2^{3N+3}}}.\tag{131}\] Finally, we have \[\label{eq:32many-point-correlations-indicator-to-sum-of-trees} \sum_{\mathbf{x}'}\mathbbold{1}(\mathbf{x}',\mathbf{A})\sum_{T'\in\mathbbold{T}(A_P)}\prod_{\{t,u\}\in E(T')}C(x_t,x_u)=\sum_{T'\in\mathbbold{T}(A)}\prod_{\{a_1,a_2\}\in E(T')}\lvert a_1-a_2\rvert^{-\alpha}.\tag{132}\]

Combining Equations 129 , 130 , 131 and 132 , we obtain Equation 122 . The lemma has thus been proved. ◻

Proof of Theorem 3. The proof of the theorem follows from a rather direct application of Lemma 24 into Equation 99 . We have \[\lvert\langle:\sigma_A:\rangle^+_\Lambda\rvert\leq \overline{\mathrm{W}}_\beta(\{A\})+\sum_{m=2}^{N}\sum_{\substack{P\in\mathbbold{P}(A)\\\lvert P\rvert =m}}\sum_{n=m}^{2m-2}\sum_{T\in\mathbbold{T}^*_{n,m}}\frac{1}{(n-m)!}\overline{\mathrm{W}}'_\beta(P,T)\eqqcolon\overline{\mathrm{W}}_\beta(\{A\})+F_\beta(A).\]Bounding the first term is a direct application of Lemma 23. On the other hand, we have the following trivial bounds: \(\lvert\mathbbold{T}^*_{n,m}\rvert\leq \lvert \mathbbold{T}_{2N}\rvert\), \(m\leq N\) and \(\lvert\mathbbold{P}(A)\rvert\leq\lvert \mathbbold{P}_{N}\rvert\). Thus, by Lemma 24, \[F_\beta(A)\leq N\lvert\mathbbold{P}_N\rvert\lvert \mathbbold{T}_{2N}\rvert F(N)e^{-\frac{\beta c_2}{2^{3N+3}}}\sum_{T'\in\mathbbold{T}(A)}\prod_{\{a_1,a_2\}\in E(T')}\lvert a_1-a_2\rvert^{-\alpha}.\]Choosing \(\beta\geq 2^{3N+4}c_2^{-1}\log[N\lvert\mathbbold{P}_N\rvert\lvert\mathbbold{T}_{2N}\rvert F(N)]\) and \(c_4=\min\{2^{-3N-4}c_2,2^{-N-5}c_0\log2\}\), we obtain Equation 2 . The theorem has thus been proved. ◻

8 Concluding Remarks↩︎

This work successfully continues the program of obtaining results for the entire region of phase transition \(1<\alpha\leq 2\), whenever they are true, and erasing the hypothesis that the interaction between first-neighbors should be arbitrarily large.

The matter of phase-separation remains open in the case \(\alpha\in (1,3-3\log2]\sqcup\{2\}\): the case \(1<\alpha <3-3\log2\) should be rather easy given [18] and this work. We point out that another relatively natural avenue of further research is the study of the decay of correlations in the presence of a random field when \(\alpha\in (1,3/2)\), especially given recent work by Ding, Huang, and Maia [16].

Acknowledgements↩︎

RB was supported by CNPq grants 311658/2025-3 and 407527/2025-7. HC was supported by CAPES grant 160295/2024-6. We would like to thank Lucas Affonso, João Maia, Pierre Picco and Kelvyn Welsch for many discussions about the long-range Ising model, the contour and cluster expansion arguments employed in studying it. Without their extensive prior work, this paper would not have been possible.

References↩︎

[1]
G. Gallavotti and S. Miracle-Sole, “Statistical mechanics of lattice systems,” Communications in Mathemtatical Physics, vol. 5, pp. 317–323, 1967.
[2]
D. Ruelle, “Statistical mechanics of a one-dimensional lattice gas,” Communications in Mathematical Physics, vol. 9, pp. 267–278, 1968.
[3]
F. J. Dyson, “Existence of a phase-transition in a one-dimensional ising ferromagnet,” Communications in Mathematical Physics, vol. 12, pp. 91–107, 1969.
[4]
F. J. Dyson, “Non-existence of spontaneous magnetization in a one-dimensional ising ferromagnet,” Communications in Mathematical Physics, vol. 12, pp. 212–215, 1969.
[5]
J. B. Rogers and C. J. Thompson, “Absence of long-range order in one-dimensional spin systems,” Journal of Statistical Physics, vol. 25, no. 4, pp. 669–678, 1981, doi: 10.1007/BF01022361.
[6]
J. Bricmont, J. Lebowitz, and C.-E. Pfister, “On the surface tension of lattice systems,” Annals of the New York Academy of Sciences, vol. 337, p. 214, 1980.
[7]
M. Fannes, P. Vanheuverzwijn, and A. Verbeure, “Energy-entropy inequalities for classical lattice systems,” Journal of Statistical Physics, vol. 29, pp. 547–560, 1982.
[8]
S. E. Burkov and Y. G. Sinai, “Phase diagrams of one-dimensional lattice models with long-range anti-ferromagnetic interaction,” Russian Mathematical Surveys, vol. 38, 1983.
[9]
J. Fröhlich and T. Spencer, “The phase transition in the one-dimensional ising model with \(1/r^2\) interaction energy,” Communications in Mathematical Physics, vol. 84, pp. 87–101, 1982.
[10]
J. Fröhlich and T. Spencer, “The kosterlitz-thouless transition in two-dimensional abelian spin systems and the coulomb gas,” Communications in Mathematical Physics, vol. 81, pp. 527–602, 1981.
[11]
J. Z. Imbrie, “Decay of correlations in the one dimensional ising model with \({J_{ij}=|i-j|^{-2}}\),” Communications in Mathematical Physics, vol. 85, pp. 491–515, 1982.
[12]
M. Aizenman, J. T. Chayes, L. Chayes, and C. M. Newman, Discontinuity of the magnetization in one-dimensional \(1/|x-y|^2\) Ising and Potts models,” Journal of Statistical Physics, vol. 50, no. 1–2, pp. 1–40, 1988, doi: 10.1007/bf01022985.
[13]
J. Z. Imbrie and C. M. Newman, “An intermediate phase with slow decay of correlations in one dimensional \({1/|x- y|^2}\) percolation, ising and potts models,” Communications in Mathematical Physics, vol. 118, pp. 303–336, 1988.
[14]
M. Cassandro, P. A. Ferrari, I. Merola, and E. Presutti, “Geometry of contours and peierls estimates in \(d=1\) ising models with long range interactions,” Journal of Mathematical Physics, vol. 46, p. 053305, 2005.
[15]
M. Cassandro, E. Orlandi, and P. Picco, “Phase transition in the 1d random field ising model with long range interaction,” Communications in Mathematical Physics, vol. 288, no. 2, pp. 731–744, 2009.
[16]
J. Ding, F. Huang, and J. Maia, “Phase transitions in low-dimensional long-range random field ising models.” arXiv, 2024, doi: 10.48550/arXiv.2412.19281.
[17]
M. Cassandro, I. Merola, P. Picco, and U. Rozikov, “One-dimensional ising models with long range interactions: Cluster expansion, phase-separating point,” Communications in Mathematical Physics, vol. 327, pp. 951–991, 2014.
[18]
M. Cassandro, I. Merola, and P. Picco, “Phase separation for the long range one-dimensional ising model,” Journal of Statistical Physics, vol. 167, pp. 351–382, 2017.
[19]
R. Bissacot, E. O. Endo, A. C. D. van Enter, B. Kimura, and W. M. Ruszel, “Contour methods for long-range ising models: Weakening nearest-neighbor interactions and adding decaying fields,” Annales Henri Poincaré, vol. 19, no. 8, pp. 2557–2574, 2018.
[20]
J. Littin and P. Picco, “Quasi-additive estimates on the hamiltonian for the one-dimensional long range ising model,” Journal of Mathematical Physics, vol. 58, 2017, doi: 10.1063/1.4994034.
[21]
L. Affonso, “Multidimensional contours à la fröhlich-spencer and boundary conditions for quantum spin systems. Ph.d. Thesis - university of são paulo.” 2023, [Online]. Available: https://arxiv.org/abs/2310.07946.
[22]
L. Affonso, R. Bissacot, E. O. Endo, and S. Handa, “Long-range ising models: Contours, phase transitions and decaying fields,” Journal of the European Mathematical Society, vol. 27, no. 4, pp. 1679–1714, 2025, doi: 10.4171/jems/1529.
[23]
L. Affonso, R. Bissacot, G. Faria, and K. Welsch, “Phase transition in long-range ferromagnetic \(q-\)state models: Contours, random and decaying fields.” 2024.
[24]
L. Affonso, R. Bissacot, and J. Maia, “Phase transitions in multidimensional long-range random field ising models.” arXiv, 2024, doi: 10.48550/arXiv.2307.14150.
[25]
L. Affonso, R. Bissacot, J. Maia, J. Rodrigues, and K. Welsch, “Cluster expansion and decay of correlations for multidimensional long-range ising models,” Preprint. 2025, [Online]. Available: https://arxiv.org/abs/2508.15666.
[26]
J. Maia, “Phase transitions in ising models: The semi-infinite with decaying field and the random field long-range. Ph.d. Thesis - University of São Paulo.” 2024, [Online]. Available: https://arxiv.org/abs/2403.04921.
[27]
L. Affonso, R. Bissacot, H. Corsini, and K. Welsch, “Phase transitions on 1d long-range ising models with decaying fields: A direct proof via contours.” 2024, [Online]. Available: https://arxiv.org/abs/2412.07098.
[28]
D. Iagolnitzer and B. Souillard, Decay of correlations for slowly decreasing potentials,” Physical Review A, vol. 16, pp. 1700–4, 1977.
[29]
R. L. Dobrushin, “Estimates of semi-invariants for the ising model at low temperatures,” in Topics in statistical and theoretical physics, vol. 177, American Mathematical Society, 1996, pp. 59–82.
[30]
R. L. Dobrushin, “Perturbation methods of the theory of gibbsian fields,” in Ecole d’Été de probabilité de saint-flour XXIV - 1994, vol. 1648, Springer-Verlag, 1996, pp. 1–66.
[31]
R. Fernández and A. Procacci, “Cluster Expansion for Abstract Polymer Models. New Bounds from an Old Approach,” Communications in Mathematical Physics, vol. 274, no. 1, pp. 123–140, 2007, doi: 10.1007/s00220-007-0279-2.
[32]
R. Kotecký and D. Preiss, “Cluster expansion for abstract polymer models,” Communications in Mathematical Physics, vol. 103, no. 3, pp. 491–498, 1986, doi: 10.1007/bf01211762.
[33]
V. A. Malyshev and R. A. Minlos, Gibbs random fields, 1991st ed. Dordrecht, Netherlands: Springer, 1991.
[34]
A. Procacci, “Cluster expansion methods in rigorous statistical mechanics.” 2023, [Online]. Available: https://arxiv.org/abs/2308.06380.
[35]
R. Bissacot, R. Fernández, and A. Procacci, “On the convergence of cluster expansions for polymer gases,” Journal of Statistical Physics, vol. 139, pp. 598–617, 2010.
[36]
S. Jansen and L. Kolesnikov, “Cluster Expansions: Necessary and Sufficient Convergence Conditions,” Journal of Statistical Physics, vol. 189, no. 3, 2022, doi: 10.1007/s10955-022-02992-6.
[37]
C. Temmel, “Properties and applications of bernoulli random fields with strong dependency graph. Ph.d thesis - aix-marseille université and graz university of technology.” 2012.
[38]
C. Temmel, “Sufficient Conditions for Uniform Bounds in Abstract Polymer Systems and Explorative Partition Schemes,” Journal of Statistical Physics, vol. 157, no. 6, pp. 1225–1254, 2014, doi: 10.1007/s10955-014-1108-6.
[39]
O. Penrose, “Convergence of fugacity expansions for fluids and lattice gases,” Journal of Mathematical Physics, vol. 4, no. 10, pp. 1312–1320, 1963, doi: 10.1063/1.1703906.
[40]
A. D. Scott and A. D. Sokal, “The Repulsive Lattice Gas, the Independent-Set Polynomial, and the Lovász Local Lemma,” Journal of Statistical Physics, vol. 118, no. 5–6, pp. 1151–1261, 2005, doi: 10.1007/s10955-004-2055-4.