Hostname: page-component-cd9895bd7-dzt6s Total loading time: 0 Render date: 2024-12-23T06:03:14.324Z Has data issue: false hasContentIssue false

Sampling from the random cluster model on random regular graphs at all temperatures via Glauber dynamics

Published online by Cambridge University Press:  09 December 2024

Andreas Galanis
Affiliation:
University of Oxford, Oxford, UK
Leslie Ann Goldberg
Affiliation:
University of Oxford, Oxford, UK
Paulina Smolarova*
Affiliation:
University of Oxford, Oxford, UK
*
Corresponding author: Paulina Smolarova; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

We consider the performance of Glauber dynamics for the random cluster model with real parameter $q\gt 1$ and temperature $\beta \gt 0$. Recent work by Helmuth, Jenssen, and Perkins detailed the ordered/disordered transition of the model on random $\Delta$-regular graphs for all sufficiently large $q$ and obtained an efficient sampling algorithm for all temperatures $\beta$ using cluster expansion methods. Despite this major progress, the performance of natural Markov chains, including Glauber dynamics, is not yet well understood on the random regular graph, partly because of the non-local nature of the model (especially at low temperatures) and partly because of severe bottleneck phenomena that emerge in a window around the ordered/disordered transition. Nevertheless, it is widely conjectured that the bottleneck phenomena that impede mixing from worst-case starting configurations can be avoided by initialising the chain more judiciously. Our main result establishes this conjecture for all sufficiently large $q$ (with respect to $\Delta$). Specifically, we consider the mixing time of Glauber dynamics initialised from the two extreme configurations, the all-in and all-out, and obtain a pair of fast mixing bounds which cover all temperatures $\beta$, including in particular the bottleneck window. Our result is inspired by the recent approach of Gheissari and Sinclair for the Ising model who obtained a similar flavoured mixing-time bound on the random regular graph for sufficiently low temperatures. To cover all temperatures in the RC model, we refine appropriately the structural results of Helmuth, Jenssen and Perkins about the ordered/disordered transition and show spatial mixing properties ‘within the phase’, which are then related to the evolution of the chain.

Type
Paper
Copyright
© The Author(s), 2024. Published by Cambridge University Press

1. Introduction

Glauber dynamics is a well-studied Markov chain that is widely used to sample from spin systems such as the Ising and Potts models. A particularly appealing feature of the chain is that it is usually very simple to implement; yet, establishing whether convergence to the equilibrium distribution is fast turns out to be significantly harder, and typically requires an in-depth understanding of the underlying distribution. Here, we focus on studying the performance of Glauber dynamics for the random cluster representation of the classical Potts model on the random regular graph, where the behaviour of the chain is underpinned by various phenomena that distinguish it from classical spin systems and pose new challenges in the analysis.

We begin with a few definitions. For real numbers $q, \beta \gt 0$ and a graph $G = (V,E)$ , the random cluster model on $G$ with parameters $q$ and $\beta$ is a probability distribution on the set $\Omega =\Omega _G$ of all assignments ${\mathcal{F}}\;:\;E\rightarrow \{0,1\}$ ; we typically refer to assignments in $\Omega$ as configurations. For a configuration $\mathcal{F}$ , we say that edges mapped to $1$ are in-edges, and edges mapped to $0$ are out-edges. We use ${\mathrm{In}}({\mathcal{F}})$ to denote the set of edges $e$ with ${\mathcal{F}}(e)=1$ , $\mathrm{Out}({\mathcal{F}})$ to denote the set of edges $e$ with ${\mathcal{F}}(e)=0$ , $|{\mathcal{F}}|$ for the cardinality of ${\mathrm{In}}({\mathcal{F}})$ and $c({\mathcal{F}})$ for the number of connected components in the graph $(V,{\mathrm{In}}({\mathcal{F}}))$ . Then, the weight of $\mathcal{F}$ in the RC model is given by $w_G({\mathcal{F}})= q^{c({\mathcal{F}})} (e^\beta -1)^{|{\mathcal{F}}|}$ .

For integer values of $q$ , the RC model is closely connected to the (ferromagnetic) Ising/Potts models; $q=2$ is the Ising model and $q\geq 3$ is the Potts model whose configurations are all possible assignments of $q$ colours to the vertices of the graph where an assignment $\sigma$ has weight proportional to $\mathrm{e}^{\beta m(\sigma )}$ with $m(\sigma )$ being the number of monochromatic edges under $\sigma$ . The RC model is an alternative edge representation of the models (for integer $q$ ) that has also been studied extensively in its own right due to its intricate behaviour (see, e.g., [Reference Grimmett18]).

We will be primarily interested in sampling from the so-called Gibbs distribution on $\Omega$ induced by these weights, denoted by $\pi _G(\cdot )$ , where for a configuration $\mathcal{F}$ , $\pi _G({\mathcal{F}}) =w_G({\mathcal{F}})/Z_G$ where the normalising factor $Z_G = \sum _{{\mathcal{F}}'\in \Omega _G}w_G({\mathcal{F}}')$ is the aggregate sum of weight of all configurations (known as the partition function). We focus on the Glauber dynamics which is a classical Markov chain for sampling from Gibbs distributions which is a particularly useful tool for developing approximate sampling algorithms. We will refer to Glauber dynamics for the RC model as the RC dynamics. Roughly, the RC dynamics is a Markov chain $(X_t)_{t\geq 0}$ initialised at some configuration $X_0$ which evolves by iteratively updating at each step $t\geq 1$ a randomly chosen edge based on whether its endpoints belong to the same component in the graph $(V,{\mathrm{In}}(X_t))$ . The mixing time of the chain is the number of steps to get within total variation distance $\leq 1/4$ from $\pi _G$ , see Section 2 for details.

Our goal is to obtain a fast algorithm for the RC model using Glauber dynamics on the random regular graph. There are two key obstacles that arise, especially at low temperatures (large $\beta$ ): (i) Glauber dynamics for the RC model has a non-local behaviour since its updates depend on the component structure of the running configuration, and (ii) there are severe bottleneck phenomena and worst-case graphs which prohibit a general fast-convergence result, and more generally an efficient algorithm. The random regular graph is a particularly interesting testbed in this front since it exhibits all the relevant phase transition phenomena and has also been used as the main gadget in hardness reductions [Reference Galanis, Štefankovic, Vigoda and Yang14].

To overview the phenomena that are most relevant for us, the following picture was detailed in a remarkable development by Jenssen, Helmuth, and Perkins [Reference Helmuth, Jenssen and Perkins21]: for $\Delta \geq 5$ and all sufficiently large $q$ , they established the ordered/disordered transition occurring at some $\beta _c$ satisfying $\beta _c=(1+o_q(1))\tfrac{2\log q}{\Delta }$ (see also [Reference Galanis, Štefankovic, Vigoda and Yang14] for integer $q\geq 3$ ).Footnote 1 Roughly, for $\beta \lt \beta _c$ a typical configuration of the model is disordered, whereas for $\beta \gt \beta _c$ it is ordered: disordered configurations resemble the all-out configuration (in that all components are of size $O(\log n)$ ) whereas ordered configurations resemble the all-in configuration (where there is a giant component with $\Omega (n)$ vertices). The two types of configurations coexist at $\beta =\beta _c$ , that is, each appears with some probability bounded away from zero. The methods in [Reference Helmuth, Jenssen and Perkins21] are based on cluster expansion techniques which also yielded an efficient sampling algorithm at all temperatures $\beta \gt 0$ . This is a surprising algorithmic result given that the coexistence causes multimodality in $\pi _G$ and severe bottleneck phenomena for Markov chains in a window around $\beta _c$ ; it was shown for instance in [Reference Helmuth, Jenssen and Perkins21] that the RC dynamics (and the related non-local Swendsen–Wang dynamics) have exponential mixing time, essentially because of the number of steps needed for the chain to move from ordered to disordered (and vice versa).

These results pose a rather bleak landscape for the RC dynamics; yet, on random regular graphs it is widely conjectured that the multimodality and the associated bottlenecks can be circumvented by initialising the chain more judiciously, in particular at either the all-out or the all-in configurations (depending on whether $\beta \leq \beta _c$ ). However the tools available for analysing Markov chains are typically insensitive to the initial configuration, and even more so when working at a critical range of the parameters.

Our main result establishes this conjecture for all $\Delta \geq 5$ and $q$ sufficiently large (conditions which we inherit from [Reference Helmuth, Jenssen and Perkins21]). For an integer $n$ such that $\Delta n$ is even, let ${\mathcal{G}}_{n,\Delta }$ denote the set of all $\Delta$ -regular graphs with $n$ vertices.Footnote 2 Throughout, we use $O(1)$ to denote a constant depending on $q,\beta, \Delta$ but independent of $n$ .

Theorem 1. Let $\Delta \geq 5$ be an integer. There exists $C=C(\Delta )\gt 0$ such that, for all sufficiently large $q$ , the following holds for any $\beta \gt 0$ , w.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ .

  1. 1. For $\beta \lt \beta _c$ , the mixing time of the RC dynamics starting from all-out is $O(n\log n).$

  2. 2. For $\beta \gt \beta _c$ , the mixing time of the RC dynamics starting from all-in is $O(n^{C})$ . For integer $q$ , the mixing time is in fact $O(n\log n)$ .

Note that Theorem1 implies a polynomial-time sampling algorithm from the Potts model for all $\beta \neq \beta _c$ (and all sufficiently large $q$ ). Intuitively, and as we will see later in more detail, Theorem1 asserts that the RC dynamics starting from all-in mixes quickly within the set of ordered configurations for $\beta \gt \beta _c$ , and similarly it mixes well within the disordered set of configurations starting from all out when $\beta \lt \beta _c$ . In fact, we show that both of these remain true even for $\beta =\beta _c$ and hence the RC dynamics can be used to sample even at criticality, starting from an appropriate mixture of the all-in and the all-out configurations. See Section 2 (and Section 2.4 in particular) for the exact statements.

Finally, note that the RC dynamics can be used analogously to Theorem1 to produce a sample within total variation distance $\varepsilon$ of $\pi _G$ for any $\varepsilon \geq \mathrm{e}^{-\Theta (n)}$ , by running it for a number of steps which is $\log (1/\varepsilon )$ times the corresponding mixing time bound.Footnote 3 The lower bound on the error comes from the total variation distance between $\pi _G$ and the conditional ‘ordered’ and ‘disordered’ configurations, see Lemma3.

1.1 Further related work

Our approach to proving Theorem1 is inspired from a recent paper by Gheissari and Sinclair [Reference Gheissari and Sinclair15] who established similar flavoured results for the Ising model ( $q=2$ ) on the random regular graph for large $\beta$ . To obtain our results for all $\beta$ , we adapt suitably their notion of ‘spatial mixing within the phase’, see Section 2.2 for details.

Among the results in [Reference Gheissari and Sinclair15], it was established that Glauber dynamics on the random regular graph, initialised appropriately, mixes in $O(n \log n)$ time when $\beta$ is sufficiently large.Footnote 4 More recently, Gheissari and Sinclair [Reference Gheissari and Sinclair16] obtained mixing-time bounds for the RC dynamics on the lattice $\mathbb{Z}^d$ under appropriate boundary conditions. They also analyse the mixing time starting from a mixture of the all-in/all-out initialisation. Note that the phase transition on grid lattices is qualitatively different than that of the random regular graph; there, instead of a window/interval of temperatures, the three points $\beta _u,\beta _u^*$ and $\beta _c$ all coincide into a single phase transition point. See also [Reference Borgs, Chayes, Helmuth, Perkins and Tetali6,Reference Helmuth, Perkins and Regts22] for related algorithmic results on $\mathbb{Z}^d$ using cluster expansion methods.

For the random regular graph, Blanca and Gheissari [Reference Blanca and Gheissari4] showed for all integer $\Delta \geq 3$ and real $q\geq 1$ that the mixing time is $O(n\log n)$ provided that $\beta \lt \beta _u(q,\Delta )$ where $\beta _u$ is the uniqueness threshold on the tree. A sampling algorithm (not based on MCMC) for $\beta \lt \beta _c(q,\Delta )$ and $q,\Delta \geq 3$ was designed by Efthymiou [Reference Efthymiou12] (see also [Reference Blanca, Galanis, Goldberg, Štefankovič, Vigoda and Yang3]), albeit achieving weaker approximation guarantees. Coja-Oghlan, Galanis, Goldberg, Raveloman, Štefankovič, and Vigoda [Reference Coja-Oghlan, Galanis, Goldberg, Ravelomanana, Štefankovič and Vigoda10] showed that, for all integer $q,\Delta \geq 3$ and $\beta \in (\beta _u,\beta _u')$ the mixing time is $\mathrm{e}^{\Omega (n)}$ where $\beta _u'=\log (1+\tfrac{q}{\Delta -2})\gt \beta _u$ is (conjectured to be) another uniqueness threshold on the tree (see [Reference Häggström20,Reference Jonasson24]). More generally, for integer $q\geq 3$ , the hardness results/techniques of [Reference Galanis, Štefankovic, Vigoda and Yang14,Reference Goldberg and Jerrum17] yield that for any $\beta \gt \beta _c$ , there are graphs $G$ where the mixing time of the RC dynamics is $\exp (n^{\Omega (1)})$ and the problem of appoximately sampling on graphs of max-degree $\Delta$ becomes #BIS-hard; on the other hand, for $\beta \leq (1-o_{q,\Delta }(1))\beta _c$ it has been shown in [Reference Carlson, Davies, Fraiman, Kolla, Potukuchi and Yap7,Reference Coulson, Davies, Kolla, Patel and Regts11] that the cluster expansion technique of [Reference Helmuth, Jenssen and Perkins21] yields a sampling algorithm on any max-degree $\Delta$ graph.

As a final note, another model of interest where analogous mixing results for Glauber dynamics (initialised appropriately) should be obtainable is for sampling independent sets on random bipartite regular graphs. However, in contrast to the RC/Potts models, the phase transition there is analogous to that of the Ising model, and hence, establishing the relevant spatial mixing properties close to the criticality threshold is likely to require different techniques, see, for example, [Reference Chen, Galanis, Štefankovič and Vigoda9] for more discussion.

In the next section, we outline the proof of Theorem1, explaining the main ingredients and showing how to combine them in order to conclude the proof. The rest of the paper is about establishing the main ingredients, see also Section 2.5 for a more detailed overview of the later parts.

1.2 Independent results of Blanca and gheissari

In an independent and simultaneous work, Blanca and Gheissari [Reference Blanca and Gheissari5] obtain related (but incomparable) results. For $\Delta \geq 3, q\geq 1$ and arbitrarily small $\tau \gt 0$ , they show for sufficiently large $\beta$ a mixing time bound of $O(n^{1+\tau })$ for the RC dynamics on the random regular graph starting from an arbitrary configuration (and obtain an analogous result for the grid and the Swendsen–Wang dynamics). Our result instead applies to all $\beta$ for the random regular graph (even the critical window) by taking into consideration the initial configuration; the two papers have different approaches to obtain the main ingredients.

2. Proof of Theorem1

We start with the formal description of the RC dynamics. Given a graph $G=(V,E)$ and an initial configuration $X_0:E\rightarrow \{0,1\}$ , the RC dynamics on $G$ is a Markov chain $(X_t)_{t\geq 0}$ on the set of configurations $\Omega _G$ . Let $p\;:\!=\;1-\mathrm{e}^{-\beta }$ and $\hat{p}\;:\!=\;\tfrac{p}{(1-p)q+p}$ (note that for $q\gt 1$ it holds that $\hat{p}\in (p/q,p)$ ). For $t\geq 0$ , to obtain $X_{t+1}$ from $X_t$ :

  1. 1. Choose u.a.r. an edge $e\in E$ . If $e$ is a cut-edge in the graph $(V,{\mathrm{In}}(X_t)\cup \{e\})$ , set $X_{t+1}(e)=1$ with probability $\hat{p}$ (and $X_{t+1}(e)=0$ otherwise). Else, set $X_{t+1}(e)=1$ with probability $p$ , and $X_{t+1}(e)=0$ otherwise.

  2. 2. Set $X_{t+1}(f)=X_t(f)$ for all $f\in E\backslash \{e\}$ .

It is a standard fact that the distribution of $X_t$ converges to the RC distribution $\pi _G$ . Let $T_{\mathrm{mix}}(G;\;X_0)=\min _{t\geq 0} \{t\mid \mathrm{dist}_{\mathrm{TV}}(X_t,\pi _G)\leq 1/4\}$ be the number of steps needed to get within total variation distance $\leq 1/4$ from $\pi _G$ starting from $X_0$ , and $T_{\mathrm{mix}}(G)=\max _{X_0}T_{\mathrm{mix}}(G;X_0)$ be the mixing time from the worst starting state.

2.1 The ordered and disordered phases on random regular graphs

We review in more detail the ordered/disordered transition, following [Reference Helmuth, Jenssen and Perkins21].

Definition 2. For $\Delta \geq 3$ , let $\eta =\eta (\Delta )\in (0,1/2)$ be a small constant (see Definition 19 ). For $G\in{\mathcal{G}}_{n,\Delta }$ , the ordered phase is the set of configurations $ \Omega ^{\mathrm{ord}}\;:\!=\; \{{\mathcal{F}}\in \Omega : |{\mathrm{In}}({\mathcal{F}})|\geq (1-\eta )|{E}|\}$ , whereas the disordered phase is the set $ \Omega ^{\mathrm{dis}}\;:\!=\; \{{\mathcal{F}}\in \Omega : |{\mathrm{In}}({\mathcal{F}})|\leq \eta |{E}| \}$ . For $q,\beta \gt 0$ , let $\pi ^{\mathrm{ord}}_G,\pi ^{\mathrm{dis}}_G$ be the conditional distributions of $\pi _G$ on $\Omega ^{\mathrm{ord}},\Omega ^{\mathrm{dis}}$ , respectively.

We will use the following result of Helmuth, Jenssen and Perkins [Reference Helmuth, Jenssen and Perkins21, Lemma9].

Lemma 3 ([Reference Helmuth, Jenssen and Perkins21, Theorem 1]). Let $\Delta \geq 5$ be an integer. Then, for all sufficiently large $q$ , there exists $\beta _c\gt 0$ satisfying $\beta _c=(1+o_q(1))\tfrac{2\log q}{\Delta }$ such that the following holds for any $\beta \gt 0$ w.h.p. for $G\sim{\mathcal{G}}_{n,\Delta }$ .

(1) \begin{equation} \mbox{if $\beta \lt \beta _c$, then $\left \| \pi _G-\pi _G^{\mathrm{dis}}\right \|_{\mathrm{TV}}=\mathrm{e}^{-\Omega (n)}$}; \qquad \mbox{if $\beta \gt \beta _c$, then $\left \| \pi _G-\pi _G^{\mathrm{ord}}\right \|_{\mathrm{TV}}=\mathrm{e}^{-\Omega (n)}$}. \end{equation}

Moreover, there exists $\zeta =\zeta (\Delta )\gt 0$ with $\zeta \lt \eta$ such that

(2) \begin{equation} \begin{aligned} \mbox{for $\beta \leq \beta _c$}, &\quad \pi _G^{\mathrm{dis}}\Big (|{\mathrm{In}}({\mathcal{F}})|\geq \zeta |E|\Big )=\mathrm{e}^{-\Omega (n)}, \quad \mbox{ and }\\[5pt] \mbox{for $\beta \geq \beta _c$}, &\quad \pi _G^{\mathrm{ord}}\Big (|{\mathrm{In}}({\mathcal{F}})|\leq (1-\zeta )|E|\Big )=\mathrm{e}^{-\Omega (n)}. \end{aligned} \end{equation}

Proof. The claims about the total variation distance are shown in [Reference Helmuth, Jenssen and Perkins21, Theorem1, Items (2), (3), (8)]. Equation (2) shows a bit of slack in the definitions of $\Omega ^{\mathrm{dis}}$ and $\Omega ^{{\mathrm{ord}}}$ that will be useful later; it follows essentially from the same theorem, we defer the details to Lemma24.

2.2 Main ingredient: weak spatial mixing within a phase

Let $G=(V,E)$ be a graph. For $v\in V$ and $r\geq 0$ , let $B_r(v)$ denote the set of all vertices in $V$ whose distance from $v$ is at most $r$ . Let $\pi =\pi _G$ be the RC distribution on $G$ and let $\pi _{{\mathcal{B}}_r^+(v)}$ be the conditional distribution of $\pi$ where all edges in $E\backslash E(B_r(v))$ are ‘in’. We define analogously $\pi _{{\mathcal{B}}_r^-(v)}$ by conditioning the edges in $E\backslash E(B_r(v))$ to be ‘out’.

Definition 4. Let $G$ be a graph with $m$ edges. Let $q,\beta \gt 0$ be reals and $r\geq 1$ be an integer. We say that the graph $G$ has WSM within the ordered phase at radius $r$ if for every $v\in V(G)$ and every edge $e$ incident to $v$ , $\Vert{\pi _{{\mathcal{B}}_r^+(v)}(e \mapsto \cdot ) - \pi ^{{\mathrm{ord}}}_G(e\mapsto \cdot ) }\Vert _{\mathrm{TV}} \leq \frac{1}{100m}$ . Analogously, we say that $G$ has WSM within the disordered phase at radius $r$ if for every $v\in V(G)$ and every edge $e$ incident to $v$ , $\Vert{\pi _{{\mathcal{B}}_r^-(v)}(e \mapsto \cdot ) - \pi ^{\mathrm{dis}}_G(e\mapsto \cdot ) }\Vert _{\mathrm{TV}} \leq \frac{1}{100m}$ .

The bulk of our arguments consists of showing the following two theorems.

Theorem 5. Let $\Delta \geq 5$ be an integer. There exists $M=M(\Delta )\gt 0$ such that for all $q$ sufficiently large, the following holds for any $\beta \geq \beta _c$ . W.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ , $G$ has WSM within the ordered phase at a radius $r$ which satisfies $r\leq \tfrac{M}{\beta }\log n$ .

The upper bound on the radius $r$ in terms of $1/\beta$ ensures that we can remove the dependence on $\beta$ of the mixing time in Theorem1 (caused by a loose bound on the mixing time on the tree, see Lemma9 below). For the disordered phase, we have

Theorem 6. For all integer $\Delta \geq 5$ , for all $q$ sufficiently large and any $\beta \leq \beta _c$ , w.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ , $G$ has WSM within the disordered phase at a radius $r$ which satisfies $r\leq \tfrac{1}{3}\log _{\Delta -1}\!n$ .

2.3 Second ingredient: local mixing on tree-like neighbourhoods

We first define a local version of RC dynamics where we perform only updates in a small ball around a vertex. Here, we need to consider the extreme boundary conditions that all vertices outside of the ball belong in distinct components (‘free boundary’) and where they belong to the same component (‘wired boundary’); we will refer to these two chains as the free and wired RC dynamics, respectively. For the random regular graph, these ‘local-mixing’ considerations are strongly connnected to the $\Delta$ -regular tree.

Formally, given a graph $G=(V,E)$ and a subset $U\subseteq V$ , let $G[U]$ be the induced subgraph of $G$ on $U$ . The tree excess of a connected graph $G$ is given by $|E|-|V|+1$ . For a vertex $v$ in $G$ and integer $r\geq 0$ , let $B_r(v)$ denote the set of vertices at distance at most $r$ from $v$ and $S_r(v)$ those at distance exactly $r$ from $v$ . For $K\gt 0$ , a max-degree $\Delta$ graph $G$ is locally $K$ -treelike if for every $v\in V$ and $r\leq \frac{1}{3}\log _{\Delta -1}|V|$ , the graph $G[B_r(v)]$ has tree excess $\leq K$ .

Lemma 7 (see, e.g., [Reference Gheissari and Sinclair15, Lemma 5.8], [Reference Blanca and Gheissari4, Fact 2.3]). For any integer $\Delta \geq 3$ , there is $K\gt 0$ such that w.h.p. $G\sim G_{n,\Delta }$ is locally $K$ -treelike.

For a graph $G$ , a vertex $\rho$ in $G$ and an integer $r\geq 1$ , the free RC dynamics on $B_r(\rho )$ is the RC dynamics where all edges outside of $B_r(\rho )$ are conditioned to be out and only edges of $G$ with both endpoints in $B_r(\rho )$ are updated.

Lemma 8 ([Reference Blanca and Gheissari4, Lemma 6.5]). Let $\Delta \geq 3$ be an integer, and $q,K\gt 1$ , $\beta \gt 0$ be reals. There exists $C\gt 0$ such that the following holds for any $\Delta$ -regular graph $G$ and integer $r\geq 1$ .

Suppose that $\rho \in V$ is such that $G[B_r(\rho )]$ is $K$ -treelike. Then, with $n=|B_r(\rho )|$ , the mixing time of the free RC dynamics on $B_r(\rho )$ is $\leq Cn\log n$ .

To define the wired RC dynamics, for a graph $G$ , a vertex $\rho$ in $G$ and an integer $r\geq 1$ , let $H$ be the graph obtained by removing all vertices and edges outside of $B_r(\rho )$ , and adding a new vertex $v_{\infty }$ connected to all vertices in $S_r(\rho )$ . The wired RC dynamics on $B_r(\rho )$ is the RC dynamics on $H$ where the edges adjacent to $v_\infty$ are conditioned to be in and only edges of $G$ with both endpoints in $B_r(\rho )$ are updated. Denote by $\hat{\pi }_{B_r(\rho )}$ the stationary distribution of the wired RC dynamics. Note that when the graph outside of $B_r(\rho )$ is connected, $\hat{\pi }_{B_r(\rho )}$ induces the same distribution as $\pi _{{\mathcal{B}}^+_r(\rho )}$ .Footnote 5

Lemma 9. Let $\Delta \geq 3$ be an integer, and $q,K\gt 1$ , $\beta \gt 0$ be reals. There exists $\hat{C}\gt 0$ such that the following holds for every $\Delta$ -regular graph $G=(V,E)$ and any integer $r\geq 1$ .

Suppose that $\rho \in V$ is such that $G[B_r(\rho )]$ is $K$ -treelike. Then, with $n=|B_r(\rho )|$ , the mixing time of the wired RC dynamics on $B_r(\rho )$ is $\leq \hat{C}n^3(q^{4}\mathrm{e}^{\beta })^{\Delta r}$ .

Lemma9 is proved in Appendix A.

Remark 10. For integer $q\gt 1$ , the mixing time bound in Lemma 9 can be improved to $O(n\log n)$ using results of [2], see Section 7.2 for details.

2.4 Proof of Theorem 1

Assuming for now the above ingredients, we will conclude the proof of Theorem1. For clarity, we break up the latter into the following pieces.

Theorem 11 (Convergence to ordered for real $q$ ). Let $\Delta \geq 5$ be an integer. Then, for all sufficiently large $q$ , there exists $C=C(\Delta )$ such that the following holds for $\beta \geq \beta _c$ , w.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ . The RC dynamics $(X_t)_{t\geq 0}$ on $G$ starting from all-in satisfies $\mathrm{dist}_{\mathrm{TV}}(X_T,\pi _G^{\mathrm{ord}})\leq 1/5$ for $T=O(n^{C})$ .

Theorem 12 (Convergence to disordered for real $q$ ). Let $\Delta \geq 5$ be an integer. Then, for all sufficiently large $q$ and $\beta \leq \beta _c$ , w.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ , the RC dynamics $(X_t)_{t\geq 0}$ on $G$ starting from all-in satisfies $\mathrm{dist}_{\mathrm{TV}}(X_T,\pi ^{\mathrm{dis}}_G)\leq 1/5$ for $T=O(n\log n)$ .

Theorem 13 (Faster convergence to ordered for integer $q$ ). Let $\Delta \geq 5$ be an integer. Then, for all sufficiently large integer $q$ and $\beta \geq \beta _c$ , w.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ , the RC dynamics $(X_t)_{t\geq 0}$ on $G$ starting from all-in satisfies $\mathrm{dist}_{\mathrm{TV}}(X_T,\pi ^{\mathrm{dis}}_G)\leq 1/5$ for $T=O(n\log n)$ .

Using these, the proof of Theorem1 can be concluded easily.

Proof of Theorem 1. Follows immediately by Theorems11-13, since by Lemma3 we have that $\left \| \pi _G-\pi _G^{\mathrm{ord}}\right \|_{\mathrm{TV}}=\mathrm{e}^{-\Omega (n)}$ when $\beta \gt \beta _c$ and $\left \| \pi _G-\pi _G^{\mathrm{dis}}\right \|_{\mathrm{TV}}=\mathrm{e}^{-\Omega (n)}$ .

We next show how to combine the ingredients of Sections 2.1-2.3 in order to conclude Theorem11. The proofs of Theorems12 and 13 are very similar, but involve one more technical tool in order to conclude the $O(n\log n)$ bounds. We defer their proof to Section 7.

Proof of Theorem 11. The argument resembles that of [Reference Gheissari and Sinclair15], a bit of care is required to combine the pieces. Consider $G=(V,E)\sim{\mathcal{G}}_{n,\Delta }$ with $n=|V|$ and $m=|E|$ . Let $q$ be sufficiently large so that both Lemma3 and Theorem5 apply; assume also that Lemma7 applies so that $G$ is locally $K$ -treelike. By Lemma3, for every $\beta \geq \beta _c$ , $\pi _G^{\mathrm{ord}}\big (|{\mathrm{In}}({\mathcal{F}})|\leq (1-\zeta )|E|\big )=\mathrm{e}^{-\Omega (n)}$ .

From now on consider arbitrary $\beta \geq \beta _c$ and set $\beta _0\;:\!=\;\log (q^{1.9/\Delta }+1)$ . Since $\beta _c=(1+o_q(1))\tfrac{2\log q}{\Delta }$ , we have that $\beta \geq \beta _0$ for all sufficiently large $q$ . Moreover, by Theorem5, $G$ has WSM within the ordered phase at radius $r$ for some $r\leq \tfrac{M}{\beta }\log n$ , where $M=M(\Delta )\gt 0$ is a constant independent of $\beta$ . Note that by taking $q$ large, we can ensure that $\beta _0$ and hence $\beta$ are at least $3M\log (\Delta -1)$ so that $r\leq \tfrac{1}{3}\log _{\Delta -1}n$ (and hence the radius- $r$ neighbourhood of an arbitrary vertex in $G$ is locally $K$ -treelike).

We will consider the RC dynamics $(X_t)_{t\geq 0}$ with $X_0$ being the all-in configuration on the edges. We will also consider the ‘ordered’ RC dynamics $\hat{X}_t$ with $\hat{X}_0\sim \pi ^{\mathrm{ord}}_G$ where we reject moves that lead to configurations outside of $\Omega ^{{\mathrm{ord}}}$ ; since $\hat{X}_0$ is stationary note that $\hat{X}_t\sim \pi ^{\mathrm{ord}}_G$ for all $t\geq 0$ . We will show how to couple these two dynamics in $O(n^C)$ steps. The theorem will follow by showing that ${\mathrm{dist}}_{\mathrm{TV}}(X_T,\hat{X}_T)\leq 1/5$ , where $T=O(n^{2+\log W})$ with $W=\Delta ^{2M/\beta _0} \mathrm{e}^{M \Delta (\Delta +1)}$ being independent of $\beta$ .

Consider the dynamics $(X_t)_{t\geq 0}$ and $(\hat{X}_t)_{t\geq 0}$ . For $t\geq 0$ , let $\mathcal{E}_t$ be the event that ${\mathrm{In}}(\hat{X}_t)\geq (1-\zeta )|E|$ and let $\mathcal{E}_{\lt t}\;:\!=\;\bigcap _{t'=0,\ldots, t-1} \mathcal{E}_{t'}$ . From Lemma3 we have that $\pi ^{{\mathrm{ord}}}_G(\mathcal{E}_t)\geq 1-\mathrm{e}^{-\Omega (n)}$ and hence by a union bound $\pi ^{{\mathrm{ord}}}_G(\mathcal{E}_{\lt t})\geq 1-t\mathrm{e}^{-\Omega (n)}$ as well.

We couple the evolution of $X_t$ and $\hat{X}_t$ using the monotone coupling, that is, at every step of the two chains choose the same edge $e_t$ to update and use the same uniform number $U_t\in [0,1]$ to decide whether to include $e_t$ in each of $X_{t+1},\hat{X}_{t+1}$ . Using the monotonicity of the model for $q\geq 1$ (and in particular that $p\gt \hat{p}$ ), under the monotone coupling, for all $t\geq 0$ such that $\mathcal{E}_{\lt t}$ holds (and hence no reject move has happened in $\hat{X}_t$ so far), we have that $\hat{X}_t\leq X_t$ (i.e., ${\mathrm{In}}(\hat{X}_t)\subseteq{\mathrm{In}}(X_t)$ ). To complete the proof, it therefore suffices to show that

(3) \begin{equation} {\mathbb{P}}(X_T\neq \hat{X}_T)\leq 1/4. \end{equation}

Consider an arbitrary time $t\geq 0$ . By a union bound, we have that

(4) \begin{equation} {\mathbb{P}}\big (X_{t}\neq \hat{X}_t\big )\leq \sum _{e}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\big )\leq m{\mathbb{P}}\big (\overline{ \mathcal{E}_{\lt t}}\big )+\sum _{e}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\mid \mathcal{E}_{\lt t}\big ). \end{equation}

Fix an arbitrary edge $e$ incident to some vertex $v$ , and let $(X^v_t)$ be the wired RC dynamics on $G[B_r(v)]$ . We couple the evolution of $(X_t^v)$ with that of $(X_t)$ and $(\hat{X}_t)$ using the monotone coupling analogously to above, where in $X_t^v$ we ignore updates of edges outside the ball $G[B_r(v)]$ ). We have $X_t^v\geq X_t$ for all $t\geq 0$ , and hence, conditioned on $\mathcal{E}_{\lt t}$ , we have that $X^v_t\geq X_t\geq \hat{X}_t$ . It follows that

\begin{align*}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\mid \mathcal{E}_{\lt t}\big )&={\mathbb{P}}(X_t(e)=1\mid \mathcal{E}_{\lt t})-{\mathbb{P}}(\hat{X}_t(e)=1\mid \mathcal{E}_{\lt t})\\[5pt] &\leq |{\mathbb{P}}(X^v_t(e)=1\mid \mathcal{E}_{\lt t})-{\mathbb{P}}(\hat{X}_t(e)=1\mid \mathcal{E}_{\lt t})|. \end{align*}

For any two events $A,B$ , we have $|{\mathbb{P}}(A)-{\mathbb{P}}(A\mid B)|\leq 2{\mathbb{P}}(\overline{B})$ , so using this for $B=\mathcal{E}_{\lt t}$ and $A$ the events $\{X^v_t(e)=1\},\{\hat{X}_t(e)=1\}$ , the triangle inequality gives

\begin{align*}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\mid \mathcal{E}_{\lt t}\big ) \leq 4{\mathbb{P}}\big (\overline{ \mathcal{E}_{\lt t}}\big )+|{\mathbb{P}}(X^v_t(e)=1)-{\mathbb{P}}(\hat{X}_t(e)=1)| \end{align*}

Note that ${\mathbb{P}}(\hat{X}_t(e)=1)=\pi ^{{\mathrm{ord}}}_G(e\mapsto 1)$ , so another application of triangle inequality gives

(5) \begin{equation} \begin{aligned}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\mid \mathcal{E}_{\lt t}\big )\leq 4{\mathbb{P}}\big (\overline{ \mathcal{E}_{\lt t}}\big )+\big |{\mathbb{P}}(X^v_t(e)&=1)-\pi _{{\mathcal{B}}^+_r(v)}(e\mapsto 1)\big |\\[5pt] &+\big |\pi _{{\mathcal{B}}^+_r(v)}(e\mapsto 1)-\pi ^{{\mathrm{ord}}}_G(e\mapsto 1)\big |. \end{aligned} \end{equation}

Since $G$ has WSM within the ordered phase at radius $r$ , we have that

(6) \begin{equation} \big |\pi _{{\mathcal{B}}^+_r(v)}(e\mapsto 1)-\pi ^{{\mathrm{ord}}}_G(e\mapsto 1)\big |\leq 1/(100 m). \end{equation}

Moreover, let $T_v$ be the mixing time of the wired RC dynamics on $G[B_r(v)]$ and let $N_v=|E(B_r(v))|\leq \Delta ^{r+1}$ . Since $r\leq \tfrac{1}{3}\log _{\Delta -1} n$ , $G[B_r(v)]$ is $K$ -treelike, so from Lemma9, with $\hat{C}=O(1)$ denoting the constant there (and absorbing a couple of factors of $\Delta$ into it),

\begin{equation*}T_v\leq \hat {C} (N_v)^3 (q^4 \mathrm {e}^\beta )^{\Delta r}\leq \hat {C} N_v\Delta ^{2r}q^{4\Delta r}\mathrm {e}^{\beta \Delta r}= \hat {C} N_v \Big (\Delta ^{2M/\beta } q^{4\Delta M/\beta }\mathrm {e}^{M \Delta }\Big )^{\log n}\leq \hat {C} N_v W^{\log n},\end{equation*}

where in the last inequality we used that $\beta \gt \beta _0$ , $\beta _0\gt \tfrac{1}{\Delta }\log q$ and $W=\Delta ^{2M/\beta _0} \mathrm{e}^{M \Delta (\Delta +1)}$ . For $T=\Theta (n^{2+\log W})$ , we have $T\geq 40T_v \tfrac{m}{N_v}\log m$ , so by Chernoff bounds, with probability $1-\exp (-n^{\Omega (1)})$ , we have at least $10T_v \log m$ updates inside $B_r(v)$ among $t=1,\ldots, T$ . For integer $k\geq 1$ the distance from stationarity after $kT_v$ steps is at most $(1/4)^k$ , we obtain

(7) \begin{equation} \big |{\mathbb{P}}(X^v_T(e)=1)-\pi _{{\mathcal{B}}^+_r(v)}(e\mapsto 1)\big |\leq \exp (-n^{\Omega (1)})+ \mathrm{e}^{-4\log m}\leq 1/m^3. \end{equation}

Plugging (6) and (7) into (5) for $t=T$ , and then back into (4), we obtain using ${\mathbb{P}}(\overline{\mathcal{E}_{\lt T}})\leq T\mathrm{e}^{-\Omega (n)}$ that ${\mathbb{P}}(X_T\neq \hat{X}_T)\leq 5m T \mathrm{e}^{-\Omega (n)}+m/m^3+1/100\leq 1/5$ , as needed.

Remark 14. The value $T$ in the proof of Theorem 11 gives the running time $O(n^{C(\Delta )})$ in the statement. In particular, $T = \Theta (n^{2 + \log W})$ where $W=\Delta ^{2M/\beta _0} \mathrm{e}^{M \Delta (\Delta +1)}$ , $\beta _0 \gt \log (q)/\Delta$ , and $M = \Theta (\Delta \log \Delta )$ (see the start of Section 5 ) so $C(\Delta ) = O(\Delta ^3 \log \Delta )$ .

To get the improved mixing time bounds of $O(n\log n)$ in Theorems12 and 13 the reasoning is very similar. The main difference from the above proof is that for any vertex $v$ the mixing time $T_v$ is bounded by $T_v=O( N_v \log N_v)$ , and therefore the above argument yields a mixing time upper bound of $O(n(\log n)^2)$ so a bit more care is needed to remove the extra $\log n$ factor using a log-Sobolev inequality, the details can be found in Section 7.3.

Finally, note that Theorems11, 12 and 13 all apply to the critical case $\beta =\beta _c$ as well. The only difference at criticality is that the two phases coexist [Reference Helmuth, Jenssen and Perkins21] (i.e., each appears with $\Omega (1)$ probability), so in order to obtain a sample from $\pi _G$ , one should output a sample for $\pi ^{{\mathrm{ord}}}_G$ with some probability $Q$ and otherwise a sample from $\pi ^{\mathrm{dis}}_G$ . The value of $Q$ can be computed in time $\tilde{O}(n^2)$ by approximating the corresponding partition functions, by using, for example, the algorithms in [Reference Chen, Galanis, Goldberg, Perkins, Stewart and Vigoda8,Reference Helmuth, Jenssen and Perkins21] (or even the RC dynamics itself). Precise results characterising the distribution of $Q$ can further be found in [Reference Helmuth, Jenssen and Perkins21, Theorems 2 & 3]; it is shown for example that $Q$ converges to $1/(q+1)$ as $q$ grows large.

2.5 Organisation of the rest of the paper

We first give a summary of how to obtain our WSM results for the ordered phase in the next section (Section 3); this is the most involved part of our arguments. Then, in Section 4, we revisit the polymer framework for the random cluster model and show some technical results that we will need to carry out the WSM proofs in Sections 5 and 6 (for the ordered and disordered regimes, respectively). We conclude with Section 7 which has some left-over technical pieces needed for the proof of Theorem1. Appendix A has, for completeness, a proof of the mixing-time bound on the tree stated in Lemma9 (uses relatively standard arguments) and Appenbdix B reviews some standard monotonicity properties of the random cluster model.

3. Proof outline of the WSM within the ordered phase

3.1 Locally tree-like expanders

Analogously to [Reference Helmuth, Jenssen and Perkins21], we work a bit more generally with $\Delta$ -regular expanders, which are also tree-like. The expansion profile of an $n$ -vertex graph $G=(V,E)$ for $\varepsilon \gt 0$ is given by

\begin{equation*}\phi _G(\varepsilon ) \;:\!=\; \min _{S\subseteq V;\ 0 \lt |S| \leq \varepsilon n} \frac {|E(S,V\backslash S)|}{\Delta |S|}.\end{equation*}

Then the classes $G_{\Delta, \delta }$ and ${\mathcal{G}}_{\Delta, \delta, K}$ are as follows.

Definition 15. Let $\Delta \geq 5$ be an integer, and $\delta \in (0,1/2), K\gt 0$ be reals. $G_{\Delta, \delta }$ is the class of $\Delta$ -regular graphs such that $\phi _G(1/2) \geq 1/10$ and $\phi _G(\delta ) \geq 5/9$ . ${\mathcal{G}}_{\Delta, \delta, K}$ is the class of all locally $K$ -treelike graphs $G\in{\mathcal{G}}_{\Delta, \delta }$ .

We use the following lemma.

Lemma 16 ([Reference Helmuth, Jenssen and Perkins21, Proposition 37]). Fix $\Delta \geq 5$ . There is a constant $\delta \in (0,1/2)$ such that w.h.p. a uniformly random $\Delta$ -regular graph belongs to ${\mathcal{G}}_{\Delta, \delta }$ .

Lemma16 and Lemma7 show that there is also a positive integer $K$ such that, w.h.p, $G\in \mathcal{G}_{\Delta, \delta, K}$ . Next we state an important property of expanders from [Reference Trevisan27].

Lemma 17 ([Reference Trevisan27, Lemma 2.3]). Let $G=(V,E)$ be a regular graph and consider $E'\subseteq E$ with $|E'|\leq \theta |E|$ for some $\theta \in (0,\phi _G(1/2))$ . Then $(V,E\backslash E')$ has a component of size at least $\big (1 - \tfrac{\theta }{2 \phi _G(1/2)}\big ) |V|$ .

We use Lemma17 to establish the existence of a giant component.

Definition 18. The size of a component of a graph is the number of vertices in the component. A giant component in an $n$ -vertex graph is a component whose size is greater than $n/2$ . Given a graph $G=(V,E)$ and a subset $F\subseteq E$ , $G[F]$ denotes the graph $(V,F)$ .

Definition 19. Fix $\Delta \geq 5$ . Fix $\delta \in (0,1/2)$ satisfying Lemma 16 . Let $\eta = \min (\delta /5, 1/100)$ .

Corollary 20. Fix integers $\Delta \geq 5$ and $K\geq 0$ and a real number $\delta \in (0,1/2)$ . Let $G$ be a graph in ${\mathcal{G}}_{\Delta, \delta, K}$ and let $\mathcal{F}$ be a configuration in $\Omega ^{\mathrm{ord}}$ or a partial configuration with $|{\mathrm{In}}({\mathcal{F}})|\geq (1-\eta )|E|$ . Then there is a giant component in $G[{\mathrm{In}}({\mathcal{F}})]$ whose size is at least $ (1-\delta )|V|$ .

Proof. Apply Lemma17 with $E' = \mathrm{Out}({\mathcal{F}})$ and $\theta = \eta = \min (\delta /5,1/100)$ . Note $|\mathrm{Out}({\mathcal{F}})|\leq \eta |E|$ and $\phi _G(1/2)\geq 1/10$ . Thus the lemma say that $G[{\mathrm{In}}({\mathcal{F}})]$ has a component of size at least $\left (1-\frac{\delta / 5}{2\cdot 1/10}\right )|V| = (1 - \delta )|V| \gt |V|/2$ .

3.2 Sketch of proof of Theorem5

Let $\Delta \geq 5$ be an integer. Consider any sufficiently large $q$ and any $\beta \geq \beta _c$ . For sufficiently large $n$ , choose a ‘radius’ $r \approx \tfrac{1}{\beta } \log n$ and let $G=(V,E)\sim{\mathcal{G}}_{n,\Delta }$ . Fix a vertex $v\in V$ and an edge $e$ incident to $v$ . We wish to show, with sufficiently high probability, that $\Vert{\pi _{{\mathcal{B}}_r^+(v)}(e \mapsto \cdot ) - \pi ^{{\mathrm{ord}}}_G(e\mapsto \cdot ) }\Vert _{\mathrm{TV}} \leq{1}/{(100 |E|)}$ .

Our goal is essentially to construct a coupling of ${\mathcal{F}}^+\sim \pi _{{\mathcal{B}}^+_{r}(v)}$ and ${\mathcal{F}}^{\mathrm{ord}}\sim \pi ^{\mathrm{ord}}$ , such that ${\mathbb{P}}({\mathcal{F}}^+(e)\neq{\mathcal{F}}^{\mathrm{ord}}(e))$ is sufficiently small. In order to construct the coupling, we take advantage of the fact that $G[B_r(v)]$ is locally tree-like. In fact, we identify a suitable subgraph of $G[B_r(v)]$ without cycles and restrict the coupling to this subgraph.

Consider a breadth-first search from $v$ in $G[B_r(v)]$ . Let $T_0$ be the rooted tree consisting of all forward edges in this breadth-first search. All other edges in $B_r(v)$ are called ‘excess edges’. W.h.p., since $G\sim{\mathcal{G}}_{n,\Delta }$ , there are at most $K$ excess edges in $B_r(v)$ for some absolute constant $K\gt 0$ . In particular, since $G$ is locally tree-like, we can identify integers $r_1$ and $r_2$ satisfying $r \geq r_1 \gt r_2 \geq 0$ such that $E(B_{r_1}(v))\setminus E(B_{r_2}(v))$ contains no excess edges and $r_1 - r_2 \geq r/(2K)=\Omega (r)$ . The fact that $r_1-r_2 = \Omega (r)$ ensures that $B_{r_1}(v) \setminus B_{r_2}(v)$ is a sufficiently large subgraph of $G$ , and the coupling focuses on this subgraph.

In order to describe the coupling process we need a small amount of notation. A partial configuration $\mathcal{F}$ is a map from the edges of $G$ to the set $\{0,1,*\}$ . In-edges and out-edges (those that are mapped to $1$ or to $0$ ) are ‘revealed’ and edges that are mapped to $*$ are ‘unrevealed’. A refinement of a partial configuration is obtained by revealing more edges. We use ${\mathcal{F}} \subseteq{\mathcal{F}}'$ to denote the fact that ${\mathcal{F}}'$ refines $\mathcal{F}$ .

In the coupling, we generate a sequence of edge subsets $F_0 \subseteq F_1 \subseteq \cdots \subseteq E$ such that, after iteration $i$ , the edges in $F_i$ are revealed. We also construct two sequences of partial configurations ${\mathcal{F}}_0^+\subseteq{\mathcal{F}}_1^+\subseteq \dots \subseteq{\mathcal{F}}^+$ and ${\mathcal{F}}_0^{\mathrm{ord}}\subseteq{\mathcal{F}}_1^{\mathrm{ord}}\subseteq \dots \subseteq{\mathcal{F}}^{\mathrm{ord}}$ , maintaining the invariant that the revealed edges in ${\mathcal{F}}_i^{\mathrm{ord}}$ and ${\mathcal{F}}_i^+$ are exactly the edges in $F_i$ . The coupling will have the crucial property that ${\mathcal{F}}^{\mathrm{ord}} \sim \pi ^{\mathrm{ord}}$ and ${\mathcal{F}}^+ \sim \pi _{{\mathcal{B}}_{r_1}^+(v)}$

  • The process starts with iteration $i=0$ . The initial set $F_0$ of revealed edges is all edges except those in $E(B_{r_1}(v))$ . In ${\mathcal{F}}_0^{\mathrm{ord}}$ these revealed edges are sampled from the projection $\pi ^{\mathrm{ord}}_{F_0}$ of $\pi ^{\mathrm{ord}}$ onto $F_0$ . It is likely that the configuration ${\mathcal{F}}_0^{\mathrm{ord}}$ has at least $(1-\eta )|E|$ in-edges. If not, then the coupling terminates (unsuccessfully), generating ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ from the right distributions. We will show that the probability of this unsuccessful termination is low. On the other hand, if ${\mathcal{F}}_0^{\mathrm{ord}}$ has a least $(1-\eta )|E|$ in-edges, then we are off to a good start. All configurations refining ${\mathcal{F}}_0^{\mathrm{ord}}$ are in $\Omega ^{\mathrm{ord}}$ , so the projection of $\pi$ and $\pi ^{\mathrm{ord}}$ onto subsequent edges that get revealed are the same (making it easier to continue the coupling). At this point ${\mathcal{F}}_0^+$ is taken to be the configuration with revealed edges $F_0$ where all revealed edges are in-edges.

  • After iteration $i=0$ , iterations continue with $i=1,2,\ldots$ until an edge is revealed whose distance from $v$ is at most $r_2$ or until the in-edges in ${\mathcal{F}}_i^{\mathrm{ord}}$ induce a giant component, and this giant component contains all vertices on the boundary of $F_i$ . We will show that it is very unlikely that an edge at distance at most $r_2$ from $v$ is reached. So it is likely the giant component in ${\mathcal{F}}_i^{\mathrm{ord}}$ contains all vertices on the boundary of $F_i$ . This is a good situation because the conditional distribution of $\pi$ , conditioned on refining ${\mathcal{F}}_i^{\mathrm{ord}}$ and the conditional distribution of ${\mathcal{F}}^+$ , conditioned on refining ${\mathcal{F}}_i^+$ induce the same distribution on edges incident to $v$ , which enables us to show that ${\mathbb{P}}({\mathcal{F}}^+(e)\neq{\mathcal{F}}^{\mathrm{ord}}(e))$ is sufficiently small.

  • The process at iteration $i+1$ is as follows. $W_i$ is taken to be the set of all vertices on the boundary of $F_i$ whose components (induced by the in-edges in ${\mathcal{F}}_i^{\mathrm{ord}}$ ) are all small. By ‘boundary’. we mean that vertices in $W_i$ are adjacent to revealed edges, and to unrevealed edges. If $W_i$ is empty, then the coupling finishes. Otherwise, a vertex $w_i\in W_i$ is chosen to be as far from $v$ as possible. The edges in the subtree of $T_0$ below the parent of $w_i$ are revealed in $F_{i+1}$ .

The main remaining ingredient in the proof is showing that the unsuccessful terminations of the coupling are unlikely. To do this, we use the polymer framework of [Reference Helmuth, Jenssen and Perkins21]. (Ordered) polymers are defined using an inductive definition. For a set of edges $A\subseteq E$ , let ${\mathcal{B}}_0(A) = A$ , and inductively for $j = 0,1,2,\dots$ define ${\mathcal{B}}_{j+1}(A)$ to be the set of all edges such that they are either in ${\mathcal{B}}_j(A)$ or edges that are incident to a vertex that has at least ${5\Delta }/{9}$ incident edges in ${\mathcal{B}}_j(A)$ . Let ${\mathcal{B}}_\infty (A) = \bigcup _{j\in \mathbb{N}}{\mathcal{B}}_j(A)$ . An ordered polymer of a configuration $\mathcal{F}$ is a connected component of $B_\infty (\mathrm{Out}({\mathcal{F}}))$ . The bulk of the work is to prove the following lemma, which is repeated in Section 5.2 (with more detail) as Lemma44.

Lemma 21. Fix $\Delta \geq 5$ and $K, M\gt 0$ . Suppose that $\beta \geq 3 M$ . Suppose that $n$ is sufficiently large so that $r\;:\!=\; \tfrac{M}\beta \log _{\Delta -1} n \gt K$ and $|B_{r}(v)| \leq 9 \Delta n/200$ . Define $r_1$ as above. Let ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ be generated by the process. Then at least one of the following conditions holds.

  1. 1. ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ agree on the edges that are incident to $v$ .

  2. 2. $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus E(B_{r_1}(v))| \lt (1-\eta ) |E|$ .

  3. 3. ${\mathcal{F}}^{\mathrm{ord}}$ contains a polymer of size at least $\frac{r}{400\Delta (1+K)}-1$ .

To complete the proof of Theorem5, we show that items 2 and 3 are unlikely. The proof that item 2 is unlikely, Lemma45, follows from the slack specified in Equation (2) of Lemma3. The proof that item 3 is unlikely, Lemma28, follows from an analysis on the size of polymers by adapting appropriately the cluster expansion techniques of [Reference Helmuth, Jenssen and Perkins21] (a bit of extra work is needed there to capture the $1/\beta$ dependence in the size of the polymer, see Lemma28).

The proof of WSM for the disordered phase (Theorem6) follows a similar strategy but the details are substantially simpler (due to a more straightforward notion of polymers), the argument can be found in Section 6.

4. Polymers for RC on expander graphs

All logarithms in this paper are with base $\mathrm{e}$ . Given a graph $G=(V,E)$ , a subset $U\subseteq V$ and a subset $F\subseteq E$ , we use $G[U]$ to denote the subgraph of $G$ induced by $U$ and we use $G[F]$ to denote the graph $(V,F)$ .

4.1 The RC model on expanders in the ordered regime

We use the following Lemma from [Reference Helmuth, Jenssen and Perkins21].

Lemma 22 ([Reference Helmuth, Jenssen and Perkins21, Lemma 8]). Let $\Delta \geq 5$ , $\delta \in (0, 1/2)$ , and $\zeta \in (0, 1/2)$ be reals. Suppose that $G=(V,E)\in{\mathcal{G}}_{\Delta, \delta }$ satisfies $|V|\geq 360/(\zeta \delta )$ and $A\subset E$ is such that $\zeta |E| \leq |A| \leq (1-\zeta )|E|$ . Then ${c(A)}/{|V|} +{|A|}/{|E|} \leq 1 - \zeta /40$ where $c(A)$ is the number of connected components in $G[A]$ .

Remark 23 (Reference Helmuth, Jenssen and Perkins21, Lemma 8). is proved for the special case of $\zeta = \min (1/100,\delta /5)$ , however the proof works for any positive $\zeta \lt 1/2$ .

Lemma 24 ([Reference Helmuth, Jenssen and Perkins21, Theorems 1 and 2]). Let $\Delta \geq 5$ be an integer. Fix $\delta \in (0,1/2)$ satisfying Lemma 16 . Let $\eta = \min \{\delta /5,1/100\}$ . Let $q_0 \geq e^{21 \Delta /\eta }$ be sufficiently large that Theorems1 and 2 of [21] hold. There is a positive number $\zeta \lt \eta$ such that w.h.p. for $G\sim{\mathcal{G}}_{n,\Delta }$ , it holds that

\begin{equation*} \begin {aligned} \mbox {for $\beta \geq \beta _c$}, &\quad \pi _G^{\mathrm {ord}}\Big (|{\mathrm {In}}({\mathcal {F}})|\leq (1-\zeta )|E|\Big )=\mathrm {e}^{-\Omega (n)},\\[5pt] \mbox {for $\beta \leq \beta _c$}, &\quad \pi _G^{\mathrm {dis}}\Big (|{\mathrm {In}}({\mathcal {F}})|\geq \zeta |E|\Big )=\mathrm {e}^{-\Omega (n)}. \end {aligned} \end{equation*}

Proof. We prove the first bound, the second is similar. Choose $\zeta = 20 \Delta /\log (q_0)$ . Note that $\zeta \lt \eta$ , as desired. Let $\Omega ' = \{{\mathcal{F}} \in \Omega ^{\mathrm{ord}}\;:\; |{\mathrm{In}}({\mathcal{F}})|\leq (1-\zeta )|E|\}$ and let

\begin{equation*}Z' = \sum _{{\mathcal {F}} \in \Omega '} w_G({\mathcal {F}}) = \sum _{{\mathcal {F}}\in \Omega '} q^{c({\mathcal {F}})} (\mathrm {e}^\beta - 1)^{|{\mathrm {In}}({\mathcal {F}})|}.\end{equation*}

Let $z = \max \{q,(\mathrm{e}^{\beta }-1)^{\Delta /2}\}$ . Since $Z \geq z^n$ every ${\mathcal{F}}\in \Omega$ has

\begin{equation*}w_G({\mathcal {F}})/Z \leq z^{c({\mathcal {F}}) + 2|{\mathrm {In}}({\mathcal {F}})|/\Delta - n} = z^{n(c({\mathcal {F}})/n + |{\mathrm {In}}({\mathcal {F}})|/|E| - 1)} .\end{equation*}

Since $\zeta \lt \eta \lt 1/2$ , for any ${\mathcal{F}}\in \Omega '$ , $\zeta |E| \leq (1-\eta )|E| \leq |{\mathrm{In}}({\mathcal{F}})| \leq (1-\zeta )|E|$ . By Lemma16, w.h.p., $G\in \mathcal{G}_{\Delta, \delta }$ . By Lemma22, for $n\geq 360/(\zeta \delta )$ , we have $c({\mathcal{F}})/n + |{\mathrm{In}}({\mathcal{F}})|/|E| \leq 1-\zeta /40$ so $Z'/Z \leq 2^{n\Delta /2} z^{n(-\zeta /40)}$ , which gives $Z'/Z = \mathrm{e}^{-\Omega (n)}$ since $q \gt \mathrm{e}^{20 \Delta /\zeta }$ . To finish, we will show that $Z/Z^{\mathrm{ord}} = O(1)$ , so that $\frac{Z'}{Z^{\mathrm{ord}}} = \frac{Z'}{Z} \times \frac{Z}{Z^{\mathrm{ord}}} = \mathrm{e}^{-\Omega (n)} O(1) = \mathrm{e}^{-\Omega (n)}$ . If $\beta \gt \beta _c$ , item 3 of [Reference Helmuth, Jenssen and Perkins21, Theorem1] says that there are positive numbers $n_0$ and $\xi$ such that, for $n\geq n_0$ , $(1/n) \log ((Z-Z^{\mathrm{ord}})/Z) \leq -\xi$ . So for $n\geq n_0$ , $Z^{\mathrm{ord}}/Z \geq 1 - e^{-\xi n}$ , which is at least $1/2$ for $n\geq 1/\xi$ . If $\beta = \beta _c$ , then items 1 and 3 of [Reference Helmuth, Jenssen and Perkins21, Theorem 2] show that $Z^{\mathrm{ord}}/Z$ converges in distribution to a random variable $Q/(Q+1)$ where (say) $Q\geq q/2$ for $q\geq q_0$ . So for sufficiently large $q$ , $Z^{\mathrm{ord}}/Z = O(1)$ .

4.2 Disordered polymers

We will use disordered polymers from [Reference Helmuth, Jenssen and Perkins21, Section 2.3]. Given a graph $G = (V,E)$ and a configuration ${\mathcal{F}}\in \Omega ^{\mathrm{dis}}$ , a disordered polymer of $\mathcal{F}$ is a connected component of $G[{\mathrm{In}}({\mathcal{F}})]$ . Since ${\mathcal{F}} \in \Omega ^{\mathrm{dis}}$ , each such connected component contains at most $\eta |E|$ edges.

Each disordered polymer $\gamma =(V(\gamma ),E(\gamma ))$ has an associated weight given by $ w_\gamma ^{\mathrm{dis}} \;:\!=\; q^{1-|V(\gamma )|} (e^{\beta }-1)^{|E(\gamma )|}$ . We observe the following relation between the weight of a disordered configuration and the weights of its polymers.

Observation 25. For ${\mathcal{F}}\in \Omega ^{\mathrm{ord}}$ , let $\Gamma ^{\mathrm{dis}}({\mathcal{F}})$ be the set of disordered polymers of $\mathcal{F}$ . Then

\begin{equation*}w_G({\mathcal {F}}) = q^n \prod _{\gamma \in \Gamma ^{\mathrm {dis}}({\mathcal {F}})} w_\gamma ^{\mathrm {dis}}.\end{equation*}

This observation follows from the fact each vertex of $G$ is contained in some polymer of $\mathcal{F}$ , each in-edge of $\mathcal{F}$ is an edge of some polymer of $\mathcal{F}$ , and each polymer of $\mathcal{F}$ corresponds to one component of $G[{\mathrm{In}}({\mathcal{F}})]$ .

Based on the work by Helmuth, Jenssen and Perkins [Reference Helmuth, Jenssen and Perkins21, Proposition 11], we can upper bound the probability that ${\mathcal{F}}\sim \pi ^{\mathrm{dis}}$ contains a large polymer.

Lemma 26. Let $\Delta \geq 5$ and $K\geq 0$ be integers, $\delta \in (0, 1/2)$ be real. Let $C$ be a constant. Then, for all sufficiently large $q$ , the following holds for all $\beta \leq \log (q^{2.1/\Delta }+1)$ . Let $n$ be sufficiently large. For $G\in{\mathcal{G}}_{\Delta, \delta, K}$ with $n=|V(G)|$ , for $s\;:\!=\;C \log n$ ,

\begin{equation*}\pi ^{\mathrm {dis}}_G\big (\mbox {${\mathcal {F}}$ contains a polymer with at least $s$ edges}\big )\leq n^{-3/2}/2.\end{equation*}

Proof. Let $q$ and $n$ be large enough so that Proposition 11 from [Reference Helmuth, Jenssen and Perkins21] applies, and so in particular [Reference Helmuth, Jenssen and Perkins21, Equation (13)] holds, and also that $s \gt 1$ and $q\geq e^{4\Delta (\frac{5}{2C}-1)}$ .

For a disordered configuration $\mathcal{F}$ , let $\Gamma ^{\mathrm{dis}}({\mathcal{F}})$ be the set of disordered polymers of $\mathcal{F}$ .

First we will bound, for a particular polymer $\gamma$ with $|E(\gamma )| \geq 1$ , ${\mathbb{P}}(\gamma \in \Gamma ^{\mathrm{dis}}(.))$ . Suppose ${\mathcal{F}}\in \Omega ^{\mathrm{dis}}$ has $\gamma \in \Gamma ^{\mathrm{dis}}({\mathcal{F}})$ . By Observation25, $w_G({\mathcal{F}}) = q^n \prod _{\gamma \in \Gamma ^{\mathrm{dis}}({\mathcal{F}})} w_\gamma ^{\mathrm{dis}}$ . Let ${\mathcal{F}}'$ be the (disordered) configuration defined by ${\mathrm{In}}({\mathcal{F}}') ={\mathrm{In}}({\mathcal{F}})\backslash E(\gamma )$ . Since $E(\gamma )\neq \emptyset$ , ${\mathcal{F}}\neq{\mathcal{F}}'$ . By construction, ${\mathcal{F}}'\in \Omega ^{\mathrm{dis}}$ . Since vertices in $\gamma$ are not incident to edges in ${\mathrm{In}}({\mathcal{F}}')$ , $\Gamma ^{\mathrm{dis}}({\mathcal{F}}') = (\Gamma ^{\mathrm{dis}}({\mathcal{F}})\setminus \{\gamma \})\cup _{v\in \gamma } \{v\}$ . A polymer $\gamma '$ consisting of a single vertex has weight $w_{\gamma '}^{\mathrm{dis}} = 1$ . Thus, $w_G({\mathcal{F}}') = w_G({\mathcal{F}})/w_\gamma ^{\mathrm{dis}}$ . Let $\Omega _\gamma$ be the set of all disordered configurations containing $\gamma$ as a polymer. Then $Z^{\mathrm{dis}} \geq \sum _{{\mathcal{F}}\in \Omega _\gamma } (w_G({\mathcal{F}}) + w_G({\mathcal{F}}')) = (1 + 1/w_\gamma ^{\mathrm{dis}})\sum _{{\mathcal{F}}\in \Omega _\gamma } w_G({\mathcal{F}})$ . Thus, the probability that ${\mathcal{F}}\sim \pi ^{\mathrm{dis}}$ contains the polymer $\gamma$ is at most $\frac{1}{1 + 1/w_\gamma ^{\mathrm{dis}}} \leq w_\gamma ^{\mathrm{dis}}$ .

Next, for $v\in V(G)$ let $p_v$ be the probability that ${\mathcal{F}}\sim \pi _G^{\mathrm{dis}}$ has a polymer with at least $s$ edges containing the vertex $v$ vertex $v$ . By a union bound, $p_v \leq P_v \;:\!=\; \sum _{\gamma \ni v: |E(\gamma )|\geq s} w_\gamma ^{\mathrm{dis}}$ . Let $\theta = \log (q)/(4 \Delta )$ . By [Reference Helmuth, Jenssen and Perkins21, Equation (13)],

\begin{equation*} \sum _{\gamma \ni v: |\gamma |\gt 1} e^{(1+\theta )|E(\gamma )|} w_\gamma ^{\mathrm {dis}} \leq \tfrac {1}{2}.\end{equation*}

Thus also

\begin{equation*}e^{(1+\theta )s} P_v = e^{(1+\theta )s}\sum _{\gamma \ni v: |E(\gamma )|\geq s} w_\gamma ^{\mathrm {dis}} \leq \sum _{\gamma \ni : |E(\gamma )|\geq s} e^{(1+ \theta )|E(\gamma )|} w_\gamma ^{\mathrm {dis}} \leq \tfrac 12.\end{equation*}

Hence $P_v \leq (1/2) e^{-(1+\theta )s}$ . By a union bound over all vertices, the probability in the statement of the lemma is at most $\sum _v P_v \leq (n/2) 2^{-(1+\theta )s}= (1/2) n^{1-C(1+\theta )}$ . The lemma follows by requiring $q$ to be sufficiently large with respect to $\Delta$ and $C$ that $C(1+\theta ) -1 \geq 3/2$ .

4.3 Ordered polymers

We use the definition of ordered polymers from [Reference Helmuth, Jenssen and Perkins21, Section 2.4.2]. Let $G=(V,E)$ be a graph of max degree $\Delta$ with $n=|V(G)|$ and $m=|E(G)|$ . For a set of edges $A\subseteq E$ , Let ${\mathcal{B}}_0(A) = A$ , and inductively for $i = 0,1,2,\dots$ define ${\mathcal{B}}_{i+1}(A)$ to be the set of all edges such that they are either in ${\mathcal{B}}_i(A)$ or edges that are incident to a vertex that has at least ${5\Delta }/{9}$ incident edges in ${\mathcal{B}}_i(A)$ . Let ${\mathcal{B}}_\infty (A) = \bigcup _{i\in \mathbb{N}}{\mathcal{B}}_i(A)$ . It is shown in [Reference Helmuth, Jenssen and Perkins21, Lemma 12] that $|{\mathcal{B}}_\infty (A)|\leq 10|{\mathcal{B}}_0(A)|$ .

An ordered polymer is a connected subgraph $\gamma =(V(\gamma ),E(\gamma ))$ of $G$ together with an edge configuration $\ell : E(\gamma )\rightarrow \{0,1\}$ subject to: (i) $\mathrm{Out}(\ell )\leq \eta m$ , and (ii) ${\mathcal{B}}_{\infty }(\mathrm{Out}(\ell ))=E(\gamma )$ . We let $E_{u}(\gamma )$ be the set of unoccupied edges $\mathrm{Out}(\ell )$ of the polymer and $c'(\gamma )$ be the number of components in the graph $G[E\backslash E_{u}(\gamma )]$ with fewer than $n/2$ vertices. The size of $\gamma$ is defined as $|E_\gamma |$ , whereas the weight of $\gamma$ is defined as $w_\gamma ^{\mathrm{ord}} \;:\!=\; q^{c'(\gamma )} (e^\beta -1)^{-E_u(\gamma )}$ .

Given a (partial) configuration $\mathcal{F}$ , the set of ordered polymers of $\mathcal{F}$ , denoted by $\Gamma ({\mathcal{F}})$ , consists of the connected components of $G[{\mathcal{B}}_\infty (\mathrm{Out}({\mathcal{F}}))]$ , each with the labelling on the edges induced by the corresponding assignment in $\mathcal{F}$ .

We will use a couple of polymer properties from [Reference Helmuth, Jenssen and Perkins21]. First, we note the following connection between the weight of a configuration $\mathcal{F}$ to the weight of its polymers $\Gamma ({\mathcal{F}})$ .

Lemma 27. Fix $\Delta \geq 5$ , $K\geq 0$ and $\delta \in (0, 1/2)$ . Let $G\in{\mathcal{G}}_{\Delta, \delta, K}$ and ${\mathcal{F}}\in \Omega ^{\mathrm{ord}}$ . Then $w_G({\mathcal{F}}) = q(\mathrm{e}^\beta - 1)^{|E|} \prod _{\gamma \in \Gamma ({\mathcal{F}})} w_\gamma ^{\mathrm{ord}}$ .

Proof. Note that $\mathrm{Out}({\mathcal{F}}) = \bigcup _{\gamma \in \Gamma ({\mathcal{F}})} E_u(\gamma )$ , so $(\mathrm{e}^\beta -1)^{|E| - \sum _{\gamma \in \Gamma ({\mathcal{F}})}|E_u(\gamma )|} = (\mathrm{e}^\beta -1)^{|{\mathrm{In}}({\mathcal{F}})|}.$ Hence, by [Reference Helmuth, Jenssen and Perkins21, Lemma21], $q^{\sum _{\gamma \in \Gamma ({\mathcal{F}})} c'(\gamma )} = q^{c(G[{\mathrm{In}}({\mathcal{F}})]) - 1}$ , so the result follows.

Lemma 28. Let $\Delta \geq 5$ and $K\geq 0$ be integers, $\delta \in (0, 1/2)$ be real. Then, for all sufficiently large $q$ , the following holds for all $\beta \geq \log (q^{1.9/\Delta }+1)$ . For $G\in{\mathcal{G}}_{\Delta, \delta, K}$ with $n=|V(G)|$ , for $s\;:\!=\;{2000}\log (n)/\beta$ ,

\begin{equation*} \pi ^{{\mathrm {ord}}}_G\big ( \mbox {${\mathcal {F}}$ contains a polymer with size at least $s$} \big )\leq 2n^{-3/2}. \end{equation*}

Proof. This proof closely follows the proof of Proposition 15 from [Reference Helmuth, Jenssen and Perkins21]. Suppose that $q$ is large enough so that $q^{1.9/\Delta }\geq (2\mathrm{e}\Delta )^{400}$ and fix arbitrary $\beta \geq \log (q^{1.9/\Delta }+1)$ .

Consider an arbitrary ordered polymer $\gamma$ . For any configuration ${\mathcal{F}}\in \Omega ^{{\mathrm{ord}}}$ with $\gamma \in \Gamma ({\mathcal{F}})$ , by Lemma27, it holds that $w_G({\mathcal{F}}) \leq q(\mathrm{e}^\beta - 1)^{|E|} w_\gamma ^{\mathrm{ord}}$ . Note that for any ${\mathcal{F}}\in \Omega ^{\mathrm{ord}}$ such that $\gamma \in \Gamma ({\mathcal{F}})$ , there is a configuration ${\mathcal{F}}'$ with $\Gamma ({\mathcal{F}}') = \Gamma ({\mathcal{F}})\setminus \{\gamma \}$ (obtained by setting ${\mathrm{In}}({\mathcal{F}}') ={\mathrm{In}}({\mathcal{F}})\cup E_u(\gamma )$ ), which therefore has weight $w_G({\mathcal{F}}') = w_G({\mathcal{F}})/w_\gamma ^{\mathrm{ord}}$ .

Let $\Omega _\gamma$ be the set of all ordered configurations with a polymer $\gamma$ . We can lower bound $Z^{\mathrm{ord}}$ by $\sum _{{\mathcal{F}}\in \Omega _\gamma } (w_G({\mathcal{F}}) + w_G({\mathcal{F}}')) = (1 + \frac{1}{w_\gamma ^{\mathrm{ord}}}) \sum _{{\mathcal{F}}\in \Omega _\gamma } w_G({\mathcal{F}})$ . Therefore,

\begin{equation*} \pi ^{\mathrm {ord}}(\gamma \in \Gamma (\cdot )) = \frac {\sum _{{\mathcal {F}}\in \Omega _\gamma } w_G({\mathcal {F}})}{Z^{\mathrm {ord}}} \leq \frac {1}{1 + \frac {1}{w_\gamma ^{\mathrm {ord}}}}\leq w_\gamma ^{\mathrm {ord}}, \end{equation*}

so we conclude that the probability that any fixed polymer $\gamma$ is a polymer of ${\mathcal{F}}\sim \pi ^{\mathrm{ord}}_G$ is at most $w_\gamma ^{\mathrm{ord}}$ .

Let $u$ be a vertex in $V(G)$ and denote by $\Gamma _u$ the set of polymers containing $u$ . Let $p_u$ be the probability that $u$ is contained in a polymer of size at least $s$ for ${\mathcal{F}}\sim \pi ^{\mathrm{ord}}_G$ . By a union bound over the polymers in $\Gamma _u$ , We have that $p_u\leq P_u$ , where

\begin{equation*} P_u\;:\!=\; \sum _{\gamma \in \Gamma _u;\ |E(\gamma )|\geq s} w_\gamma ^{\mathrm {ord}} \end{equation*}

By [Reference Helmuth, Jenssen and Perkins21, Lemma 14], $c'(\gamma )\leq 9|E_u(\gamma )|/(5 \Delta )$ . The bound on $\beta$ ensures that

\begin{equation*}w_\gamma ^{\mathrm {ord}} = q^{c'(\gamma )} (\mathrm {e}^\beta -1)^{-|E_u(\gamma )|} \leq (\mathrm {e}^\beta -1)^{(\frac {9}{5\Delta }\cdot \frac {\Delta }{1.9} - 1)|E_u(\gamma )|}\leq (\mathrm {e}^\beta -1)^{-|E_u(\gamma )|/20}.\end{equation*}

Hence we get

\begin{equation*} P_u \leq \sum _{\gamma \in \Gamma _u;\ |E(\gamma )|\geq s} (\mathrm {e}^\beta -1)^{-|E_u(\gamma )|/20}. \end{equation*}

By [Reference Helmuth, Jenssen and Perkins21, Lemma 12], for any ordered polymer $\gamma$ , $|E(\gamma )|\leq 10|E_u(\gamma )|$ . Thus

\begin{equation*} P_u \leq \sum _{k\geq s/10}\ \ \sum _{\gamma \in \Gamma _u;\ |E_u(\gamma )| = k} (\mathrm {e}^\beta -1)^{-|E_u(\gamma )|/20}. \end{equation*}

By [Reference Helmuth, Jenssen and Perkins21, Proof of Proposition 15], the number of polymers with $k$ unoccupied edges containing a particular vertex is at most $(2\mathrm{e}\Delta )^{10k}$ , hence

\begin{equation*} P_u\leq \sum _{k\geq s/10}\big ((2\mathrm {e}\Delta )^{10}(\mathrm {e}^{\beta }-1)^{-1/20}\big )^{k}.\end{equation*}

By the lower bounds on $\beta$ and $q$ , we have $\mathrm{e}^{\beta }-1\geq q^{1.9/\Delta }\geq (2\mathrm{e}\Delta )^{400}$ so

\begin{equation*} P_u\leq \sum _{k\geq s/10} (\mathrm {e}^{\beta }-1)^{-k/40} \leq \sum _{k\geq s/10} \mathrm {e}^{-\beta k/80}.\end{equation*}

The last sum is a geometric series with ratio $\mathrm{e}^{-\beta /80}\lt 1/2$ , so it is upper bounded by twice its largest term. Since $s=\tfrac{2000}{\beta }\log n$ , it follows that $P_u\leq 2/n^{5/2}$ . By a union bound over the vertices of $G$ , the probability that there is a polymer of ${\mathcal{F}}\sim \pi ^{{\mathrm{ord}}}_G$ which has size at least $s$ is therefore at most $2/n^{3/2}$ , as required.

Furthermore, we can prove the following lemma based on the proof of [Reference Helmuth, Jenssen and Perkins21, Lemma21].

Lemma 29. Fix $\Delta \geq 5$ , $K\geq 0$ and $\delta \in (0,1/2)$ . Suppose $n\geq 2/\delta$ . Let $G=(V,E)$ with $|V|=n$ be a graph in ${\mathcal{G}}_{\Delta, \delta, K}$ . Consider ${\mathcal{F}}\in \Omega ^{\mathrm{ord}}$ . Let $\kappa$ be the giant component of $G[{\mathrm{In}}({\mathcal{F}})]$ (which exists by Corollary 20 ).

  1. 1. If $u\notin \kappa$ and $e$ is incident to $u$ then $e$ is in a polymer of $\mathcal{F}$ .

  2. 2. Let $S$ be a non-empty set of edges with $|S| \lt{9\Delta n}/{200}$ . Let $\kappa '$ be a component of $G[{\mathrm{In}}({\mathcal{F}})\backslash S]$ with $|V(\kappa ')| \lt n/2$ . Then all but at most $45\Delta |S|$ vertices of $\kappa '$ are such that all their incident edges belong to a polymer of $\mathcal{F}$ .

  3. 3. For $S$ and $\kappa '$ as in Item 2 , there are at most $50\Delta ^2|S|$ components in $G[{\mathcal{B}}_\infty (\mathrm{Out}({\mathcal{F}}))]$ containing a vertex of $\kappa '$ , that is, there are at most $50\Delta ^2|S|$ polymers of $\mathcal{F}$ containing vertices of $\kappa '$ .

Proof. We start with the proof of Item 1. Consider a non-giant component $\kappa '$ of $G[{\mathrm{In}}({\mathcal{F}})]$ . Let $\tau$ be the set of vertices in $\kappa '$ which have an incident edge that is not in a polymer of $\mathcal{F}$ . For contradiction, assume that $\tau$ is non-empty. By Corollary20, $|\kappa '|\leq \delta |V|$ and thus $|\tau |\leq \delta |V|$ . The definition of $\tau$ and the fact that out-edges are in polymers imply that every edge in the cut $(\tau, V\backslash \tau )$ is in some polymer. By the expansion properties $|E(\tau, V \backslash \tau )|\geq 5\Delta |\tau |/9$ , hence there is a vertex $u\in \tau$ with at least ${5\Delta }/{9}$ incident edges in the cut, and thus, from the definition of polymers, all edges adjacent to $u$ are in polymers, contradicting that $u$ is in $\tau$ .

Before proving Item 2, we use it to prove Item 3. Let $\kappa '$ be a non-giant component in $G[{\mathrm{In}}({\mathcal{F}})\backslash S]$ . As before, let $\tau$ be the set of vertices in $\kappa '$ which have an incident edge that is not in a polymer of $\mathcal{F}$ . By Item 2, $|\tau | \leq 45\Delta |S|$ . Let $E_\tau$ be the set of edges that are incident to vertices in $\tau$ . Note that $G[V(\kappa '),E(\kappa ')]$ is connected, and removing an edge from a graph can increase the number of components by at most one. Thus, $(V(\kappa '),E(\kappa ')\setminus E_\tau )$ has at most $1+\Delta |\tau |$ components. Any edges in these components are in polymers of $\mathcal{F}$ . Thus there are at most $1 + \Delta |\tau | \leq 50\Delta ^2|S|$ components of $G[{\mathcal{B}}_\infty (\mathrm{Out}({\mathcal{F}}))]$ that contain vertices in $V(\kappa ')$ , proving Item 3.

To finish, we prove Item 2. First we bound the size of $\kappa '$ . Since ${\mathcal{F}} \in \Omega ^{\mathrm{ord}}$ , $|{\mathrm{In}}({\mathcal{F}})\backslash S| \geq (1 - \eta ) |E| - |S|$ . Since $\eta \leq 1/100$ , $\eta + \frac{|S|}{|E|} \lt 1/10 \leq \phi _G(1/2)$ , hence we can apply Lemma17 with $E' = \mathrm{Out}({\mathcal{F}}) \cup S$ and $\theta = |E'|/|E|$ . We find that $(V,{\mathrm{In}}({\mathcal{F}}) \backslash S)$ has a component of size at least $(1 - \frac{\theta }{2 \phi _G(1/2)}) n \geq (1 -{5 \theta } ) n$ so the number of vertices in $\kappa '$ is at most $5 \theta n \leq 5(\eta + |S|/|E|)n \leq \delta n + 10|S|/\Delta$ . Thus, $|\tau | \leq \delta n + 10 |S|/\Delta$ . Next let $\lambda = \lceil 10 |S|/\Delta \rceil$ and let $U$ be a subset of $\tau$ of size $|\tau | - \lambda \leq \delta n$ . Then

\begin{align*} |E(\tau, V\backslash \tau )| &\geq |E(U,V\setminus \tau )| \\[5pt] &\geq |E(U,V\setminus U)| - |E(U,\tau \setminus U)|\\[5pt] &\geq \tfrac 59 \Delta |U| - \lambda \Delta = \tfrac 59 \Delta |\tau | - 14 \lambda \Delta /9. \end{align*}

By the definition of $\tau$ , all edges in $E(\tau, V(\kappa ')\setminus \tau )$ are in polymers of $\mathcal{F}$ . Since $\kappa '$ is a component of $G[{\mathrm{In}}({\mathcal{F}})\setminus S]$ , there are at most $|S|$ edges in ${\mathrm{In}}({\mathcal{F}})$ between $\tau$ and $V\setminus V(\kappa ')$ . Thus the number of edges in $E(\tau, V\setminus \tau )$ that are in polymers is at least $|E(\tau, V\backslash \tau )| - |S|$ . This is at least $\tfrac 59 \Delta |\tau | - 14 \lambda \Delta /9 - |S| \geq \tfrac 59 \Delta |\tau | - 149 |S|/9 - 14\Delta /9 \geq \tfrac 59 \Delta |\tau | - 45 \Delta |S|/9$ , where the last inequality follows from $\Delta \geq 5$ , $|S|\geq 1$ .

Averaged over all vertices $w\in \tau$ , the number of edges in $E(\tau, V\setminus \tau )$ that are in polymers and are incident to $w$ is at least $\tfrac 59 \Delta - 45 \Delta |S|/(9\tau )$ . We may assume $\tau$ is non-empty, otherwise Item 2 directly follows. Thus, there is a vertex $w\in \tau$ such that the number of edges in $E(\tau, V\setminus \tau )$ that are in polymers and are incident to $w$ is at least $\tfrac 59 \Delta - 45 \Delta |S|/(9\tau )$ . Since this number is an integer, the number of edges in $E(\tau, V\setminus \tau )$ that are in polymers and are incident to $w$ is at least

\begin{equation*} W \;:\!=\; \Big \lceil \frac 59 \Delta - \frac {45\Delta |S|}{9|\tau |}\Big \rceil . \end{equation*}

Suppose for contradiction that $|\tau | \gt 45\Delta |S|$ . Then the term that is subtracted in the definition of $W$ is less than $1/9$ . On the other hand, if the first term, $5\Delta /9$ is not an integer, then its non-integer part is at least $1/9$ . Either way, $W = \lceil 5 \Delta /9\rceil \geq 5\Delta /9$ . As in the proof of Item 1, all incident edges of $w$ hence must be a in a polymer, contradicting that $w\in \tau$ .

5. Proving WSM within the ordered phase

In this section we will prove Theorem5. We will use the following notation.

Definition 30. Let $G=(V,E)$ be a graph. Given vertices $u,v\in V$ , let $d_G(v,u)$ be the distance in $G$ from $v$ to $u$ . For $U\subseteq V$ and $F\subseteq E$ , denote by $d_G(v,U)$ the min distance in $G$ from $v$ to any vertex in $U$ , and by $d_G(v,F)$ the min distance in $G$ to any endpoint of any edge in $F$ . Let $\partial _G(F)$ be the set of vertices in $V$ incident to both an edge in $F$ and an edge in $E\backslash F$ . For a rooted tree $T$ and a vertex $u\in V(T)$ , let $T_u$ denote the subtree of $T$ rooted at $u$ .

Fix positive integers $\Delta \geq 5$ and $K\geq 0$ . Let $M = M(K,\Delta )=2000 \times 400 \times (2+K) \Delta \log (\Delta -1)$ . Fix $q$ sufficiently large and $\beta \geq \beta _c$ . Fix a graph $G=(V,E)\in{\mathcal{G}}_{\Delta, \delta, K}$ with $n=|V|$ sufficiently large and fix a vertex $v\in V$ . Let $r \;:\!=\; \min (\tfrac{M}\beta, \tfrac 13) \log _{\Delta -1} n$ .

Definition 31. Consider a breadth-first search from $v$ in $G[B_r(v)]$ . Let $T_0$ be the rooted tree consisting of all forward edges in this breadth-first search. We refer to edges in $E(T_0)$ as the tree edges. We refer to the edges in $\mathcal{E} \;:\!=\; E(B_r(v)) \backslash E(T_0)$ as excess edges.

Observation 32. For every edge $\{u,u'\}$ in $\mathcal{E}$ , $d_{T_0}(v,u)$ and $d_{T_0}(v,u')$ differ by at most one. For any $u\in B_r(v)$ , $d_G(v,u) = d_{T_0}(v,u)$ .

Since ${\mathcal{G}}_{\Delta, \delta, K}$ , $|\mathcal{E}|\leq K$ so the ball $B_r(v)$ is close to being a tree. The next observation picks out a piece of $B_r(v)$ without excess edges.

Observation 33. Suppose $r \gt K+1$ . Then there are integers $r_1$ and $r_2$ satisfying $r \geq r_1 \gt r_2 \geq 0$ such that $E(B_{r_1}(v))\setminus E(B_{r_2}(v))$ contains no edges from $\mathcal{E}$ and $r_1 - r_2 \geq \frac{r}{1 + K} - 1$ .

Proof. For every edge $e$ of $G$ , let $d^+(v,e)$ denote the maximum distance in $G$ from $v$ to an endpoint of $e$ . Let $\mathcal{D} = \{ d^+_G(e) \mid e \in \mathcal{E}\}$ . Note that $|\mathcal{D}| \leq K$ . Thus the set $\{1,2,\ldots, \lfloor r\rfloor \} \backslash \mathcal{D}$ can be partitioned into at most $K + 1$ intervals of contiguous distances. Thus, there are integers $r_1$ and $r'_2$ with $r_1 \geq r_2'$ such that the set $\{r_2', r_2' + 1, \dots, r_1\}$ does not intersect $\mathcal{D}$ and $r_1 - r_2' + 1 \geq \frac{\lfloor r\rfloor - K}{K + 1} \geq \frac{r - 1 - K}{K + 1} = \frac{r}{K+1} - 1$ .

Let $r_2 = r_2' - 1$ . We claim that $(E(B_{r_1}(v))\backslash E(B_{r_2}(v)))\cap \mathcal{E} = \emptyset$ . To see this, consider any $ e \in \mathcal{E}$ . By the choice of $r_1$ and $ r_2$ , either $d^+(v,e)\gt r_1$ or $ d^+(v,e)\leq r_2$ . In the former case, an endpoint of $e$ is outside $V(B_{r_1}(v))$ , thus $e \notin B_{r_1}(v)$ . In the latter case, both endpoints of $e$ are in $B_{r_2}(v)$ , thus $ e\in B_{r_2}(v)$ . That concludes the proof.

5.1 Definitions/Notation

In addition to a configurations, we will work with partial configurations, where some edges are not yet revealed. We will use the symbol ‘ $*$ ’ to denote that an edge $e$ has not yet been revealed.

Definition 34. A partial configuration $\mathcal{A}$ is a function ${\mathcal{A}} : E \rightarrow \{0,1,*\}$ . We write $R({\mathcal{A}}) \;:\!=\;{\mathcal{A}}^{-1}(\{0,1\})$ and refer to it as the set of revealed edges in $\mathcal{A}$ . We denote by ${\mathrm{In}}({\mathcal{A}})$ the set of edges $e$ with ${\mathcal{A}}(e)=1$ (the ‘in’ edges) and by $\mathrm{Out}({\mathcal{A}})$ the set of edges $e$ with ${\mathcal{A}}(e)=0$ (the ‘out’ edges). We use $\Omega ^*$ to denote the set of all partial configurations.

Note that each configuration is also a partial configuration – these are precisely the partial configurations $\mathcal{A}$ with $R({\mathcal{A}})=E$ . Next we define notation for partial configurations.

Definition 35. Let ${\mathcal{A}}_1,{\mathcal{A}}_2$ be partial configurations. We say that ${\mathcal{A}}_2$ is a refinement of ${\mathcal{A}}_1$ , written ${\mathcal{A}}_1\subseteq{\mathcal{A}}_2$ , whenever ${\mathrm{In}}({\mathcal{A}}_1)\subseteq{\mathrm{In}}({\mathcal{A}}_2)$ and $\mathrm{Out}({\mathcal{A}}_1)\subseteq \mathrm{Out}({\mathcal{A}}_2)$ (and hence $R({\mathcal{A}}_1)\subseteq R({\mathcal{A}}_2)$ as well). If $R({\mathcal{A}}_1)$ and $R({\mathcal{A}}_2)$ are disjoint, the union of ${\mathcal{A}}_1$ and ${\mathcal{A}}_2$ , written ${\mathcal{A}}_1\cup{\mathcal{A}}_2$ , is the partial configuration with ${\mathrm{In}}({\mathcal{A}}_1\cup{\mathcal{A}}_2) ={\mathrm{In}}({\mathcal{A}}_1)\cup{\mathrm{In}}({\mathcal{A}}_2)$ and $\mathrm{Out}({\mathcal{A}}_1\cup{\mathcal{A}}_2)=\mathrm{Out}({\mathcal{A}}_1)\cup \mathrm{Out}({\mathcal{A}}_2)$ .

The following distributions on configurations and partial configurations will be useful.

Definition 36. Let $\mathcal{A}$ be a partial configuration in $\Omega ^*$ . Then $\Omega _{\mathcal{A}} \;:\!=\; \{{\mathcal{F}} \in \Omega \mid{\mathcal{A}} \subseteq{\mathcal{F}}\}$ will denote the set of all configurations that refine $\mathcal{A}$ and $\pi _{\mathcal{A}}$ is the conditional distribution of $\pi$ in $\Omega _{\mathcal{A}}$ , that is, for ${\mathcal{F}} \in \Omega _{\mathcal{A}}$ , $\pi _{\mathcal{A}}({\mathcal{F}}) = \pi ({\mathcal{F}}) /\sum _{{\mathcal{F}}' \in \Omega _{\mathcal{A}}} \pi ({\mathcal{F}}')$ . Similarly, if $\Omega _{\mathcal{A}} \cap \Omega ^{\mathrm{ord}}$ is non-empty, then $\pi ^{\mathrm{ord}}_{\mathcal{A}}$ is the conditional distribution of $\pi ^{\mathrm{ord}}$ in $\Omega _{\mathcal{A}}\cap \Omega ^{\mathrm{ord}}$ .

5.2 The coupling process

From now on, we assume that $n=|V(G)|$ is sufficiently large, so $r\gt K + 1$ . We define $r_1$ and $r_2$ as in Observation33. Note $r_1\geq \frac{r}{K+1}-1$ . Let $T = T_0[B_{r_1}(v)]$ be the first $r_1 + 1$ levels of the tree $T_0$ . To prove Theorem5, we construct a coupling of ${\mathcal{F}}^+\sim \pi _{{\mathcal{B}}^+_{r_1}(v)}$ and ${\mathcal{F}}^{\mathrm{ord}}\sim \pi ^{\mathrm{ord}}$ , such that ${\mathbb{P}}({\mathcal{F}}^+(e)\neq{\mathcal{F}}^{\mathrm{ord}}(e)) \leq 1/(100 |E|)$ .

To do this, we generate two sequences of coupled partial configurations, starting by revealing all edges outside the ball $B_{r_1}(v)$ and progressively revealing the edges of $T$ , starting from its leaves. We reveal edges one subtree at a time until either the boundary is contained in a giant component of in edges, or an edge incident to $B_{r_2}(v)$ is revealed. In the first case, we use Observation51 from Appendix B to couple edges incident to $v$ perfectly. In the other case, we will be able to show the existence of a large polymer (which is an unlikely event).

Formally, we will construct a sequence of edge subsets $F_0 \subseteq F_1 \subseteq \cdots \subseteq E$ , sequence of vertex subsets $V_0, V_1, \dots$ , and two sequences of partial configurations ${\mathcal{F}}_0^+\subseteq{\mathcal{F}}_1^+\subseteq \dots \subseteq{\mathcal{F}}^+$ and ${\mathcal{F}}_0^{\mathrm{ord}}\subseteq{\mathcal{F}}_1^{\mathrm{ord}}\subseteq \dots \subseteq{\mathcal{F}}^{\mathrm{ord}}$ , maintaining the following invariant for all integers $i\geq 0$ :

\begin{equation*} R({\mathcal {F}}_i^+) = R({\mathcal {F}}_i^{\mathrm {ord}})= F_i \mbox { and } {\mathrm {In}}({\mathcal {F}}_i^{\mathrm {ord}}) \subseteq {\mathrm {In}}({\mathcal {F}}_i^+) \end{equation*}

For each $i$ , all unrevealed edges of some subtree of $T$ get revealed. It will also hold, for each $i$ , that $E\setminus E(B_{r_1}(v)) \subseteq F_i$ . To facilitate this process of revealing the edges, the following definitions will be helpful.

Definition 37. Let $F$ be a subset of $E$ . Then $\Omega _{F} \;:\!=\; \{{\mathcal{A}} \in \Omega ^* \mid R({\mathcal{A}}) = F\}$ . The distribution $\pi _F$ on $\Omega _F$ is defined as follows. For ${\mathcal{A}}\in \Omega _F$ , $\pi _F({\mathcal{A}}) \;:\!=\; \pi (\Omega _{\mathcal{A}})$ . Similarly, $\pi ^{\mathrm{ord}}_F$ is the distribution on $\Omega _F$ so that $\pi _F^{\mathrm{ord}}({\mathcal{A}}) = \pi ^{\mathrm{ord}}(\Omega _{\mathcal{A}})$ .

Definition 38. Let ${\mathcal{A}}\in \Omega ^*$ be a partial configuration. Let $F$ be a subset of $E$ with $R({\mathcal{A}}) \subseteq F$ . Then $\Omega _{{\mathcal{A}},F} \;:\!=\; \{{\mathcal{A}}'\in \Omega _F \mid{\mathcal{A}}\subseteq{\mathcal{A}}', R({\mathcal{A}})=F\}$ . We define a distribution $\pi _{{\mathcal{A}},F}$ on $\Omega _{{\mathcal{A}},F}$ given by

\begin{equation*}\pi _{{\mathcal {A}},F}({\mathcal {A}}') \;:\!=\; \frac {\pi (\Omega _{{\mathcal {A}}'})}{\sum _{{\mathcal {A}}'' \in \Omega _{{\mathcal {A}},F}} \pi (\Omega _{{\mathcal {A}}''})} .\end{equation*}

We can now describe the process in detail.

  1. Step 1. Let $i \;:\!=\; 0$ . Let $F_0 \;:\!=\; E \setminus E(B_{r_1}(v))$ . Let ${\mathcal{F}}_0^{\mathrm{ord}}\sim \pi ^{\mathrm{ord}}_{F_0}$ . Let $V_0 \;:\!=\; \{ w\in V \mid d_G(v,w) = r_1\}$ .

    1. (a) If $|{\mathrm{In}}({\mathcal{F}}_0^{\mathrm{ord}})| \lt (1-\eta ) |E|$ : Generate ${\mathcal{F}}^{\mathrm{ord}} \sim \pi ^{\mathrm{ord}}_{{\mathcal{F}}_0^{\mathrm{ord}}}$ and generate ${\mathcal{F}}^+_0\sim \pi _{{\mathcal{B}}_{r+1}^+(v)}$ . Let ${\mathcal{F}}^+ \;:\!=\;{\mathcal{F}}^+_0$ . Terminate (unsuccessfully).

    2. (b) Otherwise: Define ${\mathcal{F}}_0^+$ by $R({\mathcal{F}}_0^+) ={\mathrm{In}}({\mathcal{F}}_0^+) = F_0$ .

  2. Step 2. Repeat until $d_G(F_i, v) \leq r_2$ :

    1. Let $W_i$ be the set of all vertices in $\partial _G(F_i)$ that are in a component of size less than $n/2$ in $G[{\mathrm{In}}({\mathcal{F}}_i^{\mathrm{ord}})]$ .

      1. (a) If $W_i$ is non-empty: Choose $w_i \in W_i$ to maximise $d_G(w_i,v)$ . Let $p_{i}$ be the parent of $w_i$ in $T$ . Let $F_{i+1} \;:\!=\; F_i \cup E(T_{p_i})$ . Generate, optimally coupled, ${\mathcal{F}}^{\mathrm{ord}}_{i+1}\sim \pi _{{\mathcal{F}}^{\mathrm{ord}}_i,F_{i+1}}$ and ${\mathcal{F}}_{i+1}^+ \sim \pi _{{\mathcal{F}}^+_i,F_{i+1}}$ . Let $V_{i+1} \;:\!=\; (V_{i}\backslash V(T_{p_i}))\cup \{p_i\}$ . Let $i \;:\!=\; i + 1$ .

      2. (b) Otherwise: Generate, optimally coupled, ${\mathcal{F}}^{\mathrm{ord}}\sim \pi _{{\mathcal{F}}_i^{\mathrm{ord}}}$ and ${\mathcal{F}}^+\sim \pi _{{\mathcal{F}}^+_i}$ and terminate (succesfully).

  3. Step 3. Generate, optimally coupled, ${\mathcal{F}}^{\mathrm{ord}}\sim \pi _{{\mathcal{F}}_i^{\mathrm{ord}}}$ and ${\mathcal{F}}^+\sim \pi _{{\mathcal{F}}^+_i}$ . Terminate (unsuccesfully).

We next derive some useful observations justifying the definition of the process. We will require the following definition.

Definition 39. Let $\mathcal{A}$ be a partial configuration. The boundary component set $\xi ({\mathcal{A}})$ of $\mathcal{A}$ is the set of equivalence classes corresponding to components of the boundary vertices, that is, $ \xi ({\mathcal{A}}) = \{ \kappa \cap \partial (R({\mathcal{A}})) \mid \kappa \text{ is a component in } G[{\mathrm{In}}({\mathcal{A}})]\}$ .

Observation 40. Suppose that the process does not terminate in Step 1a (with $i=0$ ). Then it maintains the following invariants for all $i\geq 0$ such that ${\mathcal{F}}_i^{\mathrm{ord}},{\mathcal{F}}_i^+, F_i$ and $V_i$ are defined.

  1. Inv1(i). $R({\mathcal{F}}_i^{\mathrm{ord}}) = R({\mathcal{F}}_i^+)=F_i$ , ${\mathrm{In}}({\mathcal{F}}_i^{\mathrm{ord}}) \subseteq{\mathrm{In}}({\mathcal{F}}_i^+)$ , and $E\setminus E(B_{r_1}(v))\subseteq F_i$ ,

  2. Inv2(i). For all distinct $u_1, u_2\in V_i$ , $V(T_{u_1})\cap V(T_{u_2}) = \emptyset$ ,

  3. Inv3(i). $\bigcup _{u\in V_i} V(T_u) = V(F_i)\cap V(T)$ and $\bigcup _{u\in V_i} E(T_u) \subseteq F_i$ ,

  4. Inv4(i). If $d_G(F_i, v) \gt r_2$ , then $\partial _G (F_i)\subseteq V_i$ ,

  5. Inv5(i). If $F_{i+1}$ is defined, then $F_i \subset F_{i+1}$ , and

  6. Inv6(i). If ${\mathcal{F}}_{i+1}^{\mathrm{ord}}$ and ${\mathcal{F}}_{i+1}^+$ are defined in the process, then ${\mathcal{F}}^{\mathrm{ord}}_i \subseteq{\mathcal{F}}^{\mathrm{ord}}_{i+1}$ and ${\mathcal{F}}^+_i \subseteq{\mathcal{F}}^+_{i+1}$

Proof. We do an induction on $i$ . For the base case let $i = 0$ . Inv1( $0$ ) holds by the definitions of $F_0$ , ${\mathcal{F}}_0^{\mathrm{ord}}$ and ${\mathcal{F}}_0^+$ . Inv2( $0$ ) follows from the fact that any $u\in V_0$ is a leaf of $T$ and thus $T_u$ contains a single vertex - $u$ .

The second part of Inv3( $0$ ) is trivial as neither of the subtrees have any edges. We show the first part of Inv3( $0$ ) by proving that for all $u\in V(T)$ , $u\in V_0$ if and only if $u\in V(F_0)$ . First, suppose that $u\in V(F_0)\cap V(T)$ . The definitions of $F_0$ and $T$ imply that $d_G(u,v) = r_1$ so $u\in V_0$ , as desired. For the other direction, consider $u\in V_0$ . We need to show that $G$ has an edge $\{u,u'\}$ with $d_G(u',v) \gt r_1$ . The vertex $u$ is a leaf of $T$ , hence $u$ has a single incident edge in $E(T)$ . Given that $G$ is $\Delta$ -regular with $\Delta \gt 1$ , there is some edge $\{u, u'\}$ , such that $u'$ is not $u$ ’s parent in $T$ . By Observation32, $d_G(u', v) \in \{r_1 - 1, r_1, r_1 + 1\}$ . By the choice of $r_1$ , $ r_2$ , and $u'$ , $\{u, u'\} \notin E(B_{r_1}(v))\setminus E(B_{r_2}(v))$ , so $d_G(u', v) = r_1 + 1 \gt r_1$ , and the result follows.

The invariant Inv4( $0$ ) follows from the fact that for any $u\in \partial _G (F_0)$ there are edges $\{u, u_1\} \in F_0$ and $\{u,u_2\}\notin F_0$ , so $d_G(u, v) \leq r_1$ and $d_G(u_2, v) \leq r_1$ , and thus also $d_G(u_1, v) \gt r_1$ , so by Observation32 it follows that $d_G(u, v) = r_1$ and thus $u\in V_0$ .

For the inductive step, suppose $F_{i+1}$ , ${\mathcal{F}}_{i+1}^{\mathrm{ord}}$ , $F_{i+1}^+$ and $V_{i+1}$ are all defined and suppose that Inv1( $i$ )-Inv4( $i$ ) hold. We now show Inv5( $i$ ), Inv6( $i$ ) and Inv1( $i+1$ )–Inv4( $i+1$ ).

For Inv5( $i$ ), $F_{i}\subseteq F_{i+1}$ follows by construction. For the strict inequality, it is enough to show that $\{w_i, p_i\}\in F_{i+1}\setminus F_i$ . Since $F_{i+1}$ is defined, $d_G(F_i, v) \gt r_2$ , otherwise the process would terminate before the $(i+1)$ th iteration of Step 3. Thus by Inv4( $i$ ) and by $w_i\in \partial _G (F_i)$ , $w_i\in V_i$ . So for any $u\in V_i$ , Inv2( $i$ ) implies that $p_i\notin V(T_u)$ . Hence Inv3( $i$ ) implies that $p_i\notin V(F_i)$ so $\{w_i, p_i\}\notin F_i$ . But $E(T_{p_i})\subseteq F_{i+1}$ , hence $\{w_i, p_i\}\in F_{i+1}$ , completing the proof of Inv5( $i$ ). Inv6( $i$ ) follows by construction.

The first and third parts of Inv1( $i+1$ ) follow by construction, the second part follows from Corollary53 (see Appendix B, using ${\mathrm{In}}({\mathcal{F}}_i^{\mathrm{ord}})\subseteq{\mathrm{In}}({\mathcal{F}}_i^+)$ and the fact that the coupling is optimal. Inv2( $i + 1$ ) follows by construction and Inv2( $i$ ). Also, the second part of Inv3( $i + 1$ ) follows by construction and Inv3( $i$ ).

For the first part of Inv3( $i+1$ ), by the definition of $F_{i+1}$ , if $u\in V(F_{i+1})\cap V(T)$ , then $u\in V(T_{p_i})$ or $u\in V(F_{i})\cap V(T)$ . In the first case, note $p_i\in V_{i+1}$ . Otherwise, by Inv3( $i$ ), there is a vertex $u'\in V_i$ with $u\in V(T_{u'})$ . If $u'$ is in $V(T_{p_i})$ , then so is $u$ , hence the first case applies. If not, then $u'\in V_{i+1}$ . For the converse, for $u\in V_{i+1}$ , there are two cases. If $u = p_i$ then $E(T_{p_i})\subseteq F_{i+1}$ so $V(T_{p_i}) \subseteq V(F_{i+1})\cap V(T)$ . Otherwise, $u\in V_i$ , so the result follows by Inv3( $i$ ) and Inv5( $i$ ).

Finally we prove Inv4( $i+1$ ). Suppose that $d_G(F_{i+1}, v) \gt r_2$ and let $u\in \partial _G F_{i+1}$ . Our goal is to show $u\in V_{i+1}$ Let $u_1, u_2\in V$ be such that $\{u, u_1\}\in F_{i+1}$ and $\{u, u_2\}\notin F_{i+1}$ . By Inv1( $i+1$ ), $E\setminus E(B_{r_1}(v))\subseteq F_{i+1}$ , and so $u, u_2 \in V(B_{r_1}(v)) = V(T)$ . Also, $d_G(u, v) \gt r_2$ implies $\{u, u_2\}\in E(B_{r_1}(v))\setminus E(B_{r_2}(v)$ . Since $\{u, u_1\}\in F_{i+1}$ , $u\in V(F_{i+1})$ , so by Inv3( $i+1$ ), there is a vertex $u'\in V_{i+1}$ such that $u\in V(T_{u'})$ . By the choice of $r_1$ and $r_2$ , $\{u, u_2\}$ is a tree edge, so if $u\neq u'$ , by the second part of Inv3( $i+1$ ), $\{u, u_2\}$ would be then in $F_{i+1}$ , contradicting the choice of $u_2$ . Thus $u = u' \in V_{i+1}$ .

Observation 41. The process eventually terminates, with ${\mathcal{F}}^{\mathrm{ord}} \sim \pi ^{\mathrm{ord}}$ and ${\mathcal{F}}^+\sim \pi _{{\mathcal{B}}_{r_1}^+(v)}$ .

Proof. By Inv5, the process always terminates. If the process terminates in Step 1a, then, by construction, the resulting configurations are drawn from the correct distributions. Otherwise, Step 2 is executed. Thus, $(1-\eta ) |E| \leq |{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}}_0)|$ . By Inv6, for every value of $i$ until termination, $(1-\eta ) |E| \leq |{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}}_0)| \leq |{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}}_i)|$ . This implies that every $\Omega _{{\mathcal{F}}_i^{\mathrm{ord}}} \subseteq \Omega ^{\mathrm{ord}}$ so $\pi _{{\mathcal{F}}_i^{\mathrm{ord}}}^{\mathrm{ord}} = \pi _{{\mathcal{F}}_i^{\mathrm{ord}}}$ . From this it follows that ${\mathcal{F}}^{\mathrm{ord}} \sim \pi ^{\mathrm{ord}}$ . Given how ${\mathcal{F}}_0^+$ is generated in Step 1b, it holds that, for any ${\mathcal{F}}\in \Omega$ , $\pi _{{\mathcal{B}}_{r_1}^+(v)}({\mathcal{F}}) \gt 0$ precisely when ${\mathcal{F}}_0^+\subseteq{\mathcal{F}}$ . Then ${\mathcal{F}}^+\sim \pi _{{\mathcal{B}}_{r_1}^+(v)}$ follows from the way that ${\mathcal{F}}^+_i$ is generated in the process.

Observation 42. If the process terminates in Step 2a, then there is a non-negative integer $i$ such that $G[{\mathrm{In}}({\mathcal{F}}_i^{\mathrm{ord}})]$ has a giant component containing all vertices in $\partial _G(F_i)$ . Also, ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ agree on the edges that are incident to $v$ .

Proof. The only way that the process can get to Step 2a, with some index $i$ , is if $d_G(F_i, v) \gt r_2$ and $W_i = \emptyset$ . In this case all vertices of the boundary $\partial _G(F_i)$ are contained in the giant component of $G[{\mathrm{In}}({\mathcal{F}}_i^{\mathrm{ord}})]$ .

Since $r_2 \geq 0$ , $d_G(F_i, v) \gt r_2 \geq 0$ , so the edges that are incident to $v$ are not in $F_i$ . Since $G$ is connected, the boundary $\partial _G (F_i)$ is non-empty.

Thus $\xi ({\mathcal{F}}_{i}^{\mathrm{ord}})$ consists of a single equivalence class containing all vertices of $\partial _G(F_i)$ . By Inv1( $i$ ), ${\mathrm{In}}({\mathcal{F}}_i^{\mathrm{ord}})\subseteq{\mathrm{In}}({\mathcal{F}}_i^+)$ and $R({\mathcal{F}}_i^{\mathrm{ord}}) = R({\mathcal{F}}_i^+)$ , so $G[{\mathrm{In}}({\mathcal{F}}_i^+)]$ has a giant component containing all vertices in $\partial _G(F_i)$ and $\xi ({\mathcal{F}}_i^+)$ contains a single equivalence class, containing all vetices in $\partial _G(F_i)$ . By Observation51, $\pi _{{\mathcal{F}}_i^{\mathrm{ord}}}$ and $\pi _{{\mathcal{F}}_i^+}$ have the same projection on $E\backslash F_i$ , meaning that for any partial configuration ${\mathcal{A}}\in \Omega _{E\backslash F_i}$ , $\pi _{{\mathcal{F}}_i^{\mathrm{ord}}}({\mathcal{F}}_i^{\mathrm{ord}} \cup{\mathcal{A}}) = \pi _{{\mathcal{F}}_i^+}({\mathcal{F}}_i^+ \cup{\mathcal{A}})$ . Since the coupling in Step 2a is optimal, ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ agree on the edges incident to $v$ (since they are not in $F_i$ ).

Observation 43. Suppose that Step 2a is executed with index $i$ and let $k = r_1 - d_G(w_i, v)$ . There are non-negative integers $j_0,\ldots, j_k$ with $j_k=i$ such that, for every $\ell \in \{0,\ldots, k-1\}$ , $j_{\ell } \lt j_{\ell +1}$ and $w_{j_{\ell +1}} = p_{j_\ell }$ .

Proof. By construction, $w_i\in \partial _G(F_i)$ and $d_G(F_i,v) \gt r_2$ . By Inv4( $i$ ), $\partial _G(F_i) \subseteq V_i\subseteq V(T)$ . Hence $d_G(w_i, v) \leq r_1$ , so $k \geq 0$ . The proof is by induction on $k$ .

For the base case, $k = 0$ , it suffices to define $j_0 \;:\!=\; i$ and the condition in the statement of the observation is vacuous.

For the induction step, fix $k\geq 0$ and assume that the observation holds whenever Step 2a is executed with any index $j$ such that $r_1 - d_G(w_j, v)\leq k$ . Suppose that Step 2a is executed with index $i$ such that $r_1 - d_G(w_i, v) = k+1$ . Define $j_{k+1} \;:\!=\; i$ . By Invariant Inv4( $i$ ), $w_i\in V_i$ . By the construction of $V_i$ , either $d_G(w_i, v) = r_1$ , or there is an index $i'\lt i$ such that $w_i = p_{i'}$ . The first case is ruled out by $d_G(w_i, v) = r_1 - (k+1) \lt r_1$ . In the second case, $w_i$ is the parent of $w_{i'}$ in the BFS tree, so $d_G(w_{i'}, v) = d_G(w_i, v) + 1$ and $k= r_1 - d_G(w_{i'},v)$ . We finish by applying the inductive hypothesis with index $i'$ .

Equipped with these results, we can now prove the following lemma.

Lemma 44. Fix $\Delta \geq 5$ and $K, M\gt 0$ . Suppose that $\beta \geq 3 M$ . Suppose that $n$ is sufficiently large so that $r\;:\!=\; \tfrac{M}\beta \log _{\Delta -1} n \gt K$ . and $|B_{r}(v)| \leq 9 \Delta n/200$ . Define $r_1$ as in Observation 33 so that $r_1\geq \frac{r}{K+1}-1$ . Let ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ be generated by the process. Then at least one of the following conditions holds.

  1. 1. ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ agree on the edges that are incident to $v$ .

  2. 2. $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus E(B_{r_1}(v))| \lt (1-\eta ) |E|$ .

  3. 3. ${\mathcal{F}}^{\mathrm{ord}}$ contains a polymer of size at least $\frac{r}{400\Delta (1+K)}-1$ .

Proof. By Observation41, the process eventually terminates. If it terminates in Step 1a, then $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus E(B_{r_1}(v))| \lt (1-\eta )|E|$ , thus Condition 2 is satisfied. If the process terminates in Step 2a, then by Observation42, Condition 1 holds.

For the rest of the proof we assume that the process terminates in Step 3. We will show that Condition 3 holds. Let $i'$ be the index that triggers the condition $d_G(F_{i'}, v) \leq r_2$ in Step 2. Let $i \;:\!=\; i' - 1$ so that step 2a is run for the last time with index $i$ . By the choice of $i'$ , $d_G(F_i, v) \gt r_2$ . Since $p_{i}$ is an endpoint of an edge in $F_{i'}$ , $d_G(p_i,v) \leq r_2$ . Since $w_i\in \partial _G(F_i)$ , $d_G(v,w_i) \gt r_2$ . Since $d_G(w_i, v) = d_G(p_i, v) + 1$ , we conclude that $d_G(w_i, v) = r_2 + 1$ . By Observation43 with index $i$ and $k = r_1 - (r_2+1)$ , there are non-negative integers $\{j_0,\ldots, j_k\}$ with $j_k=i$ such that, for every $\ell \in \{0,\ldots, k-1\}$ , $j_\ell \lt j_{\ell +1}$ and $p_{j_\ell } = w_{j_{\ell +1}}$ . Thus $w_{j_{\ell +1}} \in V_{\ell +1}$ . Note that $w_{j_0},\ldots, w_{j_k}$ is a path in $T$ and that $w_{j_k}$ is closest to the root $v$ in this path. By Observation33 and the definition of $k$ , $k \geq r/(K+1)-2$ . Let $r' = \lfloor k / 2\rfloor$ , $J' = \{j_{\ell } \mid r' \leq \ell \leq k\}$ and $W' = \{w_j \mid j\in J'\}$ . Since $V(T) \subseteq V(B_{r+1}(v))$ , we have established that for all $w\in W'$ $r_2+1 \leq d_G(v,w) \leq r_1$ . We now distinguish three cases.

Case 1. For all $j\in J'$ , $w_j$ is not in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ .

Let $F'$ be the set of edges with endpoints in $W'$ and note that $F'$ contains the edges in the path $w_{r'},\ldots, w_{k}$ . By Observation41, ${\mathcal{F}}^{\mathrm{ord}} \in \Omega ^{\mathrm{ord}}$ . Applying Item 1 of Lemma29 to each vertex $u\in W'$ , we find that every edge in $F'$ is contained in a polymer of ${\mathcal{F}}^{\mathrm{ord}}$ . Since the vertices in $W'$ form a path in $G$ , the edges in $F'$ are all in the same polymer of ${\mathcal{F}}^{\mathrm{ord}}$ . Thus, ${\mathcal{F}}^{\mathrm{ord}}$ has a polymer of size at least $k/2 \geq r/(2(K+1))-1$ .

Case 2. There is an index $j \in J'$ such that $w_j$ is in the giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ and the giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus \{w_j, p_j\}]$ .

We make the following claim.

Claim 2a: There is then a non-empty set $S$ of at most $9 \Delta |V|/200$ edges of $G$ such that the size of the component $\kappa '$ of $w_j$ in the graph $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\backslash S]$ is less than $n/2$ and furthermore $|V(\kappa ')| \geq |S|r'$ .

Before proving Claim 2a, we show that it implies Condition 3, completing the proof of Case 2. By Lemma29, item 2, all but at most $45\Delta |S|$ vertices of $\kappa '$ are such that all their incident edges belong to polymer. Thus, at least $|S|r' - 45\Delta |S|$ vertices of $\kappa '$ are such that all their incident edges are in a polymer. By Lemma29, item 3, there are at most $50\Delta ^2|S|$ polymers of ${\mathcal{F}}^{\mathrm{ord}}$ containing vertices of $\kappa '$ . We conclude that there is a polymer of ${\mathcal{F}}^{\mathrm{ord}}$ containing at least $\frac{|S|r' - 45\Delta |S|}{50\Delta ^2|S|}$ vertices of $\kappa '$ with the property that all of their incident edges are in the polymer.

Since the number of edges incident to a set of size $z$ is at least $\Delta z/2$ , there is a polymer of ${\mathcal{F}}^{\mathrm{ord}}$ containing a vertex of $B_{r_1}(v)$ with size at least

\begin{equation*}\frac {\Delta }{2} \times \frac {|S|r' - 45\Delta |S|}{50\Delta ^2|S|} \geq \frac {r' - 45\Delta }{100\Delta } \geq \frac {k/2 - 50\Delta }{100\Delta } \geq \frac {r}{400\Delta (1+K)} - 1.\end{equation*}

To conclude the proof of Case 2, we prove Claim 2a. Since the process ends in Step 3, Step 1a does not occur, so $|{\mathrm{In}}(F^{\mathrm{ord}}_0)|\geq (1-\eta )|E|$ . Since ${\mathcal{F}}^{{\mathrm{ord}}}_0\subseteq{\mathcal{F}}^{\mathrm{ord}}_j$ from Inv6, ${\mathrm{In}}(F^{{\mathrm{ord}}}_0)\subseteq{\mathrm{In}}(F^{\mathrm{ord}}_j)$ , so $|{\mathrm{In}}(F^{\mathrm{ord}}_j)|\geq (1-\eta )|E|$ . By Corollary20, $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}}_j)]$ has a (unique) giant component. Call this component $\kappa _j$ . Similarly, $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ has a giant component, call this component $\kappa$ . Note that $V(\kappa _j)\subseteq V(\kappa )$ . Since $w_j$ was chosen in Step 2a, $w_j\in V(\kappa )\backslash V(\kappa _j)$ .

Let $\Gamma$ be the set of all paths from $w_j$ to $V(\kappa _j)$ in the graph $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ . Every path $\gamma \in \Gamma$ contains an edge in $E\setminus F_j$ . Let $\{u_1^\gamma, u_2^\gamma \}$ be the first such edge in $\gamma$ , with $d_G(u_1^\gamma, w_j) \lt d_G(u_2^\gamma, w_j)$ . Let $S\;:\!=\; \{ \{u_1^\gamma, u_2^\gamma \} \mid \gamma \in \Gamma \}$ .

By the definition of $S$ , any path from $w_j$ to $\kappa _j$ goes through an edge in $S$ , thus $w_j$ is not connected to $\kappa _j$ in $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus S]$ . Let $\kappa '$ be the component of $w_j$ in $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus S]$ . The component $\kappa '$ contains no vertices of $\kappa _j$ . Since $|V(\kappa _j)|\gt n/2$ , we get $|V(\kappa ')| \lt n/2$ .

We conclude by we proving that $|V(\kappa ')\cap B_{r_1}(v)|\geq |S|r'$ . The proof follows from a sequence of claims.

Claim 2b: For any $\gamma \in \Gamma$ , $u_1^\gamma \in \partial _G(F_j)$ and $d_G(v,u_1^\gamma ) \gt r_2$ .

The proof of Claim 2b is as follows. If $u_1^\gamma = w_j$ , then $u_1^\gamma \in \partial _G(F_j)$ by the choice of $w_j$ . Otherwise, since $u_1^\gamma$ is the first edge on $\gamma$ that is not in $E\setminus F_j$ , the edges on the path from $w_j$ to $u_1^\gamma$ are in $F_j$ , so there is a at least one edge from $F_j$ incident to $u_1^\gamma$ , and the edge $\{u_1^\gamma, u_2^\gamma \}\notin F_j$ is also incident to $u_1^\gamma$ , so $u_1^\gamma \in \partial _G(F_j)$ . In either case, $d_G(v,u_1^\gamma ) \geq d_G(v,F_j)\gt r_2$ .

Claim 2c: For any $\gamma \in \Gamma$ , $d_G(v,u_1^\gamma ) \leq r_1$ , $d_G(v,u_2^\gamma ) \leq r_1$ , $u_1^\gamma \in V_j$ , and $u_2^\gamma$ is the parent of $u_1^\gamma$ in $V(T)$

To prove Claim 2c, consider $\gamma \in \Gamma$ . Since the edge $\{u_1^\gamma, u_2^\gamma \} \in E\setminus F_j$ , both $u_1^\gamma$ and $u_2^\gamma$ are in $V(B_{r_1}(v))=V(T)$ . Since Step 2 is executed with index $j$ , $d_G(F_j,v) \gt r_2$ . So by Claim 2b and Inv4( $j$ ), $u_1^\gamma \in V_j$ . By Invariant Inv3( $j$ ), $E(T_{u_1^\gamma })\subseteq F_j$ , so the only candidate for $u_2^\gamma$ is the parent of $u_1^\gamma$ in $T$ .

Claim 2d: $|S| = |\{u_1^\gamma \mid \gamma \in \Gamma \}|$ .

By Claim 2c, the mapping of a vertex $u \in \{u_1^\gamma \mid \gamma \in \Gamma \}$ to its parental edge in $T$ is a bijection between $ \{u_1^\gamma \mid \gamma \in \Gamma \}$ and $S= \{ \{u_1^\gamma, u_2^\gamma \} \mid \gamma \in \Gamma \}$ .

Claim 2e: For any $\gamma \in \Gamma$ , $|V(\kappa ')\cap V(T_{u_1^\gamma })| \geq r'$ .

To prove Claim 2e, consider $\gamma \in \Gamma$ . We consider two cases. First, suppose $u_1^\gamma \neq w_j$ . By construction $w_j \in V_j$ and from Claim 2c, $u_1^\gamma \in V_j$ . By Inv2( $j$ ), $V(T_{w_j})\cap V(T_{u_1^\gamma }) = \emptyset$ . At the beginning of the proof, we established $r_2+1 \leq d_G(v,w_j) \leq r_1$ . By Claims 2b and 2c, $r_2 + 1 \leq d_G(v,u_1^\gamma ) \leq r_1$ . Since the path $\gamma$ goes from $w_j$ to $u_1^\gamma$ without leaving $F_j$ , it goes from $w_j$ to a leaf of $T_{w_j}$ and it finishes by going from a leaf of $T_{u_1^\gamma }$ to $u_1^\gamma$ . All of these edges are in $\kappa '$ and there are at least $r_1 - d_G(v,u_1^\gamma )$ of them. Since, by Claim 2b, $u_1^\gamma \in \partial _G(F_j)$ but, by construction, it is not in $\kappa _j$ , by the definition of $W_j$ , $u_1^\gamma \in W_j$ . By the choice of $w_j$ , $d_G(u_1^\gamma, v) \leq d_G(w_j, v) \leq r_1 - r'$ . Thus $|V(\kappa ')\cap V(T_{u_1^\gamma })|\geq r_1 - (r_1 - r') = r'$ .

For the second case, suppose $u_1^\gamma = w_j$ . By the assumption of Case 2, $w_j$ is in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus \{w_j, p_j\}]$ , thus $S$ contains an edge other than the edge $\{u_1^{\gamma '},u_2^{\gamma '}\}$ with $\gamma ' \in \Gamma$ that is not equal to the edge $\{w_j,p_j\}$ . By Claim 2c, $u_1^{\gamma '} \neq w_j$ . Applying the argument from the first case to $\gamma '$ , the path $\gamma '$ starts by going from $w_j$ to a leaf of $T_{w_j}$ so it contains at least $r'$ edges of $T_{w_j}$ , all of which are in $\kappa _j$ .

We now use Claims 2d and 2e to finish the proof that $|V(\kappa ')| \geq |S| r'$ , which completes the proof of Claim 2a, and hence the proof of Case 2. Consider paths $\gamma$ and $\gamma '$ in $\Gamma$ with $u_1^\gamma \neq u_1^{\gamma '}$ . By Inv2( $j$ ), $V(T_{u_1^\gamma })\cap V(T_{u_1^{\gamma '}})=\emptyset$ , therefore $|V(\kappa ')| \geq \sum _{u\in \{u_1^\gamma \mid \gamma \in \Gamma \}} |V(\kappa ')\cap V(T_{u_1^\gamma })|$ . Claim 2e shows that each term in the sum is at least $r'$ . Claim 2d shows that $|S|$ is equal to the number of terms in the sum. Therefore we obtain $|V(\kappa ')| \geq |S| r'$ , as required.

Case 3. There is an index $j \in J'$ such that $w_j$ is in the giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ . However, for all $j\in J'$ , $w_j$ is not in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus \{w_j, p_j\}]$ .

Let $\kappa$ be the giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ . Let $j_\ell = \min \{j \in J' \mid w_j \in V(\kappa )\}$ .

Claim 3: For each $\ell '$ satisfying $\ell \leq \ell ' \lt k$ , $w_{j_{\ell '}}\in V(\kappa )$ and $\{w_{j_{\ell '}},w_{j_{\ell '+1}}\} \in{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})$ .

We prove Claim 3 by induction on $\ell '$ . The base case is $\ell ' = \ell$ . In this case, $w_{j_{\ell }}\in V(\kappa )$ is from the definition of $j_\ell$ . Then since $w_{j_\ell }$ is in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ but is not in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\backslash \{w_{j_\ell }, p_{j_{\ell }}\}]$ , the edge $\{w_{j_\ell }, p_{j_{\ell }}\} = \{w_{j_\ell }, w_{j_{\ell +1}}\} \in{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})$ , as required. For the induction step, fix $\ell '$ satisfying $\ell \leq \ell ' \lt k-1$ and assume Claim 3 for $\ell '$ . Since Claim 3 implies $w_{j_{\ell '}}\in \kappa$ and $\{w_{j_{\ell '}}, w_{j_{\ell '+1}}\}\in{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})$ , $w_{j_{\ell '}}$ and $w_{j_{\ell '+1}}$ are in the same component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ , namely, $\kappa$ . Then since $w_{j_{\ell '+1}}$ is in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})]$ but is not in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\backslash \{ w_{j_{\ell '+1}}, p_{j_{\ell '+1}}\}]$ , the edge $\{w_{j_{\ell '+1}}, p_{j_{\ell '+1}}\} = \{w_{j_{\ell '+1}}, w_{j_{\ell '+2}}\} \in{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})$ , as required.

Claim 3 implies that for every $j$ satisfying $j_\ell \leq j \leq k$ , $w_j\in V(\kappa )$ . Let $e = \{w_k,p_k\}$ . Since $w_k$ is not in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus \{e\}]$ and the edge $e$ is not on the path $w_{j_\ell },\ldots, w_k$ it follows that no vertices in $W'$ are in a giant component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus \{e\}]$ . We now consider two cases.

Case 3a. $k - \ell \geq r'/2$ .

Take $S = \{e\}$ . Suppose that $|V|$ is sufficiently large that $|S| \lt 9\Delta |V|/200$ . For any $j\in J'$ let $\kappa _j$ be the component containing $w_j$ in $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\setminus \{w_j,p_j\}]$ . Item 2 of Lemma29 implies that all but at most $45\Delta$ vertices of $\kappa _j$ are such that all their incident edges belong to a polymer of ${\mathcal{F}}^{\mathrm{ord}}$ . Item 3 implies that there are at most $50 \Delta ^2$ polymers containing vertices of $\kappa _j$ . Since $k-\ell \geq r'/2$ , we conclude that there are at least $ r'/2 - 45\Delta$ vertices of $W'\cap V(\kappa )$ such that all their incident edges are in a polymer, and at most $50\Delta ^2$ polymers containing them. Hence there a polymer of size at least $\frac{\Delta }{2} \times \frac{r'/2 - 45\Delta }{50\Delta ^2}\geq \frac{r' - 90\Delta }{200\Delta } \geq \frac{r}{400\Delta (1+K)}-1$ .

Case 3b. $k - \ell \lt r'/2$ .

We have already seen that $w_{r'},\ldots, w_{\ell -1}$ is a path in $G$ with length at least $\ell - r' \geq r'/2$ and its vertices are not in $V(\kappa )$ . By Item 1 of Lemma29, every edge incident to one of these vertices is contained in a polymer of ${\mathcal{F}}^{\mathrm{ord}}$ . Since the vertices form a path in $G$ , the edges adjacent to them are all in the same polymer of ${\mathcal{F}}^{\mathrm{ord}}$ . Thus, ${\mathcal{F}}^{\mathrm{ord}}$ has a polymer of size at least $(r'/2) \geq \frac{(k/2)-1}{2} \geq \frac{r}{4(K+1)}-1$ .

Next we will show that the termination condition from Step 1 of the process is unlikely to happen.

Lemma 45. Let $\Delta \geq 5$ and $K\geq 0$ be integers. Let $\delta \in (0,1/2)$ be a real number. Let $\eta = \min \{\delta /5,1/100\}$ There are positive real numbers $q_0\geq 1$ and $n_0$ such that the following holds for all $q\geq q_0$ and all $\beta \geq \beta _c$ . Let $G=(V,E)$ . Let $n=|V|$ and let $v$ be a vertex in $V$ . Let $r$ be a real number satisfying $r \leq \tfrac{1}{3}\log _{\Delta -1}(n)$ . Then $\pi _G^{\mathrm{ord}}(|{\mathrm{In}}({\mathcal{F}})\backslash E(B_r(v))| \lt (1-\eta )|E|) = e^{-\Omega (n)}$ .

Proof. Let $r_0(n) = \tfrac{1}{3} \log _{\Delta -1}(n)$ so $r(n) \leq r_0(n)$ . Let $n_0$ be sufficiently large that $\Delta ^{r(n_0)+1} \leq (\eta - \zeta ) n_0\Delta /2$ , where $\zeta$ is the constant from Lemma3. By Lemma3,

\begin{equation*}\pi _G^{\mathrm {ord}}\Big (|{\mathrm {In}}({\mathcal {F}})|\leq (1-\zeta )|E|\Big )=\mathrm {e}^{-\Omega (n)}.\end{equation*}

However, if $|{\mathrm{In}}({\mathcal{F}})| \gt (1-\zeta )|E|$ , then

\begin{equation*}|{\mathrm {In}}({\mathcal {F}})\backslash E(B_r(v))| \gt (1-\zeta )|E| - \Delta ^{r(n)+1} \geq (1-\zeta )|E| - \Delta ^{r(n_0)+1}\geq (1-\eta )|E|.\end{equation*}

5.3 Proof of WSM within the ordered phase

We can now prove Theorem5. We start with the following Lemma.

Lemma 46. Let $\Delta \geq 5$ and $K\geq 0$ be integers and let $\delta \in (0, 1/2)$ be a real. There is $M=M(\Delta, K)\gt 0$ such that the following holds for all sufficiently large $q$ and any $\beta \geq \beta _c$ . For all sufficiently large $n$ and any $n$ -vertex graph $G\in{\mathcal{G}}_{\Delta, \delta, K}$ , the RC model on $G$ with parameters $q$ and $\beta$ has WSM within the ordered phase at radius $r_1$ satisfying $r_1\leq \frac{{M}}{\beta }\log _{\Delta -1} (n)$ .

Proof. Recall that $\eta = \min \{\delta /5,1/100\}$ . Let $M = 2000 \times 400 \times (2+K) \Delta \log (\Delta -1)$ .

Let $r(n,\beta ) \;:\!=\; \frac{{M}}{\beta }\log _{\Delta -1}(n)$ . Let $\beta _0(q) \;:\!=\; \log (q^{1.9/\Delta }+1)$ . Let $q_0$ and $n_0$ be large enough that

  • $\beta _c(q_0) \geq \beta _0(q_0)$ .

  • Lemmas28, 45 and 44 apply

  • $\beta _0(q_0) \geq 3 M$ .

  • $2 n_0^{-3/2} + f(n_0) \leq 1/(50 \Delta n_0)$ where $f(n)$ is the $e^{-\Omega (n)}$ upper bound from Lemma45.

For every $\beta$ define $n_1(\beta )$ to be the smallest positive integer such that

  • $\frac{r(n_1(\beta ),\beta )}{400\Delta (1+K)}-1 \geq \frac{r(n_1(\beta ),\beta )}{400\Delta (2+K)}$

Now fix $q\geq q_0$ and $\beta \geq \beta _c(q)$ .

Consider any $n\geq \max \{n_0,n_1(\beta )\}$ . Consider an $n$ -vertex graph $G\in{\mathcal{G}}_{\Delta, \delta, K}$ . Let $r_1$ be as in Lemma44. Use the process to generate ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ . We wish to show that for every edge $e$ incident to $v$ ,

\begin{equation*} \Vert {\pi _{{\mathcal {B}}_{r_1}^+(v)}(e \mapsto \cdot ) - \pi ^{{\mathrm {ord}}}(e\mapsto \cdot ) }\Vert _{\mathrm {TV}} \leq 1/(100|E|) = 1/(50\Delta n) \end{equation*}

By Observation41 ${\mathcal{F}}^{\mathrm{ord}} \sim \pi ^{\mathrm{ord}}$ and ${\mathcal{F}}^+\sim \pi _{{\mathcal{B}}_{r+1}^+(v)}$ . By Lemma44, if ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ do not agree on an edge incident to $v$ , then ${\mathcal{F}}^{\mathrm{ord}}$ contains a polymer of size at least $\frac{r(n,\beta )}{400\Delta (1+K)}-1$ , or $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{ord}})\backslash E(B_{r_1}(v))| \lt (1-\eta )|E|$ . Thus the probability that ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ do not agree on an edge incident to $v$ is at most the sum of probabilities of those two events.

Applying Lemma28, noting that $\frac{r(n,\beta )}{400\Delta (1+K)}-1 \geq \frac{r(n,\beta )}{400\Delta (2+K)}$ and (given the lower bound on $M$ at the start of the proof) that this quantity is at least $2000\log n/\beta$ , the probability of the first event is at most $2n^{-3/2}$ .

Note that $r(n,\beta ) \leq \tfrac 13 \log _{\Delta -1}(n)$ . Applying Lemma45, the probability of the second event is $f(n) = e^{-\Omega (n)}$ . The result follows by summing these probabilities.

Finally, we use Lemma46 to prove Theorem5.

Proof. By Lemmas7 and 16, there exist $K=K(\Delta )\gt 0$ and $\delta =\delta (\Delta ) \in (0,1/2)$ such that, w.h.p., $G\sim{\mathcal{G}}_{\Delta, n}$ is in ${\mathcal{G}}_{\Delta, \delta, K}$ . Let $M = M(\Delta, K)$ be as in Lemma46, and note that $M$ is in fact only dependent on $\Delta$ . By Lemma46, for all $q$ sufficiently large and $\beta \geq \beta _c$ , for all $n$ sufficiently large and any $n$ -vertex $G\in{\mathcal{G}}_{\Delta, \delta, K}$ , $G$ has a WSM within the ordered phase at radius $r_1$ where $r_1 \leq \tfrac{M}\beta \log _{\Delta -1}(n)$ .

6. Proof of WSM within the disordered phase

We will prove that for $q$ large enough, all large enough graphs in ${\mathcal{G}}_{\Delta, \delta, K}$ have WSM within the disordered phase on all low enough temperatures. We use the same notion of partial configurations as in Section 5.1.

Definition 47. Let $\mathcal{A}$ be a partial configuration in $\Omega ^*$ . If $\Omega _{\mathcal{A}} \cap \Omega ^{\mathrm{dis}}$ is non-empty, then $\pi ^{\mathrm{dis}}_{\mathcal{A}}$ is the conditional distribution of $\pi ^{\mathrm{dis}}$ in $\Omega _{\mathcal{A}}\cap \Omega ^{\mathrm{dis}}$ .

For a set $F\subseteq E$ , define the distribution $\pi ^{\mathrm{dis}}_F$ to be the distribution on $\Omega _F$ , such that for ${\mathcal{A}}\in \Omega _F$ , $\pi _F^{\mathrm{dis}}({\mathcal{A}}) = \pi ^{\mathrm{dis}}(\Omega _{\mathcal{A}})$ .

First we prove the following lemma showing that a suitable coupling exists.

Lemma 48. Fix $\Delta \geq 5$ , $K\geq 0$ integers, and $\delta \in (0, 1/2)$ a real. Then, for all $q$ large enough and $\beta \leq \log (q^{2.1/\Delta }+1)$ , the following holds for all $G=(V,E)\in{\mathcal{G}}_{\Delta, \delta, K}$ with sufficiently many vertices. Let $v\in V$ and $r \;:\!=\; \frac{1}{3}\log _{\Delta -1} |V|$ . Then there is an integer $r_1$ satisfying $r \geq r_1 \geq \frac{r}{K+1}-1$ and there is a coupling $({\mathcal{F}}^-,{\mathcal{F}}^{\mathrm{dis}})$ such that ${\mathcal{F}}^-\sim \pi _{{\mathcal{B}}^-_{r_1}(v)}$ and ${\mathcal{F}}^{\mathrm{dis}}\sim \pi ^{\mathrm{dis}}$ , and moreover at least one of the following holds:

  1. 1. $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}}\backslash E(B_{r_1}(v))| \gt \eta |E| - |E(B_{r_1}(v))|$ .

  2. 2. For any edge $e$ incident to $v$ , ${\mathcal{F}}^-(e) ={\mathcal{F}}^{\mathrm{dis}}(e)$ .

  3. 3. $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}})]$ has a component with at least $\frac{r}{K+1}-2$ edges.

Proof. We use the definition of $T_0$ and $\mathcal{E}$ from Definition31. By Observation33, we get $r_1$ and $r_2$ satisfying $r\geq r_1 \gt r_2 \geq 0$ such that $r_1-r_2\geq \frac{r}{K+1}-1$ and $E(B_{r_1}(v))\setminus E(B_{r_2}(v))$ contains no edges from $\mathcal{E}$ . Let $T$ be the tree consisting of the first $r_1 + 1$ levels of $T_0$ , that is, $T_0[B_{r_1}(v)]$ .

Construct the coupling as follows: First, let $F_0\;:\!=\; E\setminus E(B_{r_1}(v))$ and ${\mathcal{F}}_0^{\mathrm{dis}}\sim \pi _{F_0}^{\mathrm{dis}}$ . If $|{\mathrm{In}}({\mathcal{F}}_0^{\mathrm{dis}})| \gt \eta |E| - |E(B_{r_1}(v))|$ , then let ${\mathcal{F}}^{\mathrm{dis}}\sim \pi ^{\mathrm{dis}}_{{\mathcal{F}}_0^{\mathrm{dis}}}$ and ${\mathcal{F}}^-\sim \pi ^{\mathrm{dis}}_{{\mathcal{B}}_{r_1}^-(v)}$ . Otherwise, let ${\mathcal{F}}_0^-$ be the partial configuration with $R({\mathcal{F}}_0^-) = \mathrm{Out}({\mathcal{F}}_0^-) = F_0$ , let $F_1 = E\setminus E(B_{r_2+1}(v))$ and generate, optimally coupled, ${\mathcal{F}}_1^{\mathrm{dis}}\sim \pi _{{\mathcal{F}}_0^{\mathrm{dis}},F_1}$ and ${\mathcal{F}}_1^-\sim \pi _{{\mathcal{F}}_0^-,F_1}$ . Finally, generate, optimally coupled ${\mathcal{F}}^-\sim \pi _{{\mathcal{F}}^-_1}$ and ${\mathcal{F}}^{\mathrm{dis}}\sim \pi _{{\mathcal{F}}^{\mathrm{dis}}_1}$ .

First note, that ${\mathcal{F}}^{\mathrm{dis}}\sim \pi ^{\mathrm{dis}}$ and ${\mathcal{F}}^-\sim \pi _{{\mathcal{B}}_{r_1}^-(v)}$ . In the case $|{\mathrm{In}}({\mathcal{F}}_0^{\mathrm{dis}})| \gt \eta |E| - |E(B_{r_1}(v))|$ it follows by construction. Otherwise, we may note that any refinement of ${\mathcal{F}}^{\mathrm{dis}}_0$ has at most $\eta |E|$ edges, as $|E\setminus R({\mathcal{F}}_0^{\mathrm{dis}})| = |E(B_{r_1}(v))|$ , thus $\pi _{{\mathcal{F}}_0^{\mathrm{dis}},F_1} = \pi _{{\mathcal{F}}_0^{\mathrm{dis}},F_1}^{\mathrm{dis}}$ . Similarly, $\pi _{{\mathcal{F}}^{\mathrm{dis}}_1} = \pi _{{\mathcal{F}}_1^{\mathrm{dis}}}^{\mathrm{dis}}$ . For ${\mathcal{F}}^-$ the result follows from the fact that a configuration $\mathcal{F}$ is a refinement of ${\mathcal{F}}_0^-$ if and only if $\pi _{{\mathcal{B}}^-_{r_1}(v)}({\mathcal{F}})\gt 0$ .

Now we proceed to show that the resulting configurations satisfy one of the conditions. We have three exhaustive cases, and we show that each corresponds to one of the conditions.

Case 1. $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}}_0)| \gt \eta |E| - |E(B_{r_1}(v))|$ .

Then its refinement ${\mathcal{F}}^{\mathrm{dis}}$ satisfies $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}})\backslash E(B_{r_1}(v))| \gt \eta |E| - |E(B_{r_1}(v))|$ , since $E(B_{r_1}(v))\cap R({\mathcal{F}}^{\mathrm{dis}}_0) = \emptyset$ , thus the first condition holds.

Case 2. $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}}_0)| \leq \eta |E| - |E(B_{r_1}(v))|$ , and $\xi ({\mathcal{F}}_1^{\mathrm{dis}})$ is a free boundary. That is, no two vertices in $\partial (F_1)$ are in the same component of $G[{\mathrm{In}}({\mathcal{F}}_1^{\mathrm{dis}})]$ .

Since $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}}_0)| \leq \eta |E| - |E(B_{r_1}(v))|$ , we have $R({\mathcal{F}}_0^{\mathrm{dis}}) = R({\mathcal{F}}_0^-) = F_0$ , and $R({\mathcal{F}}_1^{\mathrm{dis}}) = R({\mathcal{F}}_1^0) = F_1$ . Also, $\emptyset ={\mathrm{In}}({\mathcal{F}}_0^-)\subseteq{\mathrm{In}}({\mathcal{F}}_0^{\mathrm{dis}})$ , thus by Corollary53, and by the fact that ${\mathcal{F}}_1^{\mathrm{dis}}$ and ${\mathcal{F}}_1^-$ were optimally coupled, also ${\mathrm{In}}({\mathcal{F}}_1^-)\subseteq{\mathrm{In}}({\mathcal{F}}_1^{\mathrm{dis}})$ . So since no two vertices in $\partial (R({\mathcal{F}}_1^{\mathrm{dis}})) = \partial (F_1) = \partial (R({\mathcal{F}}_1^-))$ are in the same component of $G[{\mathrm{In}}({\mathcal{F}}_0^{\mathrm{dis}})]$ , it is also true that no two vertices in $\partial (F_1)$ are in the same component of $G[{\mathrm{In}}({\mathcal{F}}_0^-)]$ . So $\xi ({\mathcal{F}}^{\mathrm{dis}}_0) = \xi ({\mathcal{F}}^-_0)$ , thus by Observation51 and by the fact that they were optimally coupled, ${\mathcal{F}}^-$ and ${\mathcal{F}}^{\mathrm{dis}}$ agree on the edges of $E\setminus F_1 = E(B_{r_2+1}(v))$ . Note that since $r_2 + 1 \geq 1$ , all edges incident to $v$ are in $E(B_{r_2+1}(v))$ , thus ${\mathcal{F}}^{\mathrm{dis}}$ and ${\mathcal{F}}^-$ agree on them. This corresponds to condition two.

Case 3. $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}}_0)| \leq \eta |E| - |E(B_{r_1}(v))|$ and there are distinct vertices $u_1, u_2\in \partial (F_1)$ such that $u_1$ and $u_2$ are in the same component of $G[{\mathrm{In}}({\mathcal{F}}_1^{\mathrm{dis}})]$ .

Consider any path $\gamma$ between $u_1$ and $u_2$ in the $G[{\mathrm{In}}({\mathcal{F}}_1^{\mathrm{dis}})]$ . Since no edges in $B_{r_2}(v)$ are revealed in ${\mathcal{F}}_1^{\mathrm{dis}}$ , $\gamma$ is contained in $E\setminus E(B_{r_2}(v))$ . Note that since $T_{u_1}$ and $T_{u_2}$ - the subtrees of $T$ rooted in $u_1$ and $u_2$ , respectively, are disjoint, $V(\gamma )$ is not contained in $V(T_{u_1})$ . Let $e = \{w_1, w_2\}$ be the first edge on $\gamma$ such that $w_1\in V(T_{u_1})$ and $w_2\not \in V(T_{u_1})$ . Note that for any vertex $w$ with $r_1 \gt d_G(w, v) \gt r_2$ , all edges incident to $w$ are edges of $T$ , and in particular, if it also holds that if $w$ is a non-leaf vertex in $V(T_{u_1})$ , all incident edges to $w$ are in $E(T_{u_1})$ .

Then, by construction and the choice of $r_1$ and $ r_2$ , $\partial (F_1) = \{u\in V \mid d_G(u, v) = r_2 + 1\}$ . Also no edges incident to vertices with distance at most $r_2$ from $v$ were revealed, so it must be the case that $d_G(w_1, v) = r_1$ .

Hence there is a path containing at least $r_1 - (r_2 + 1) \geq \frac{r}{K+1} - 2$ edges in $G[{\mathrm{In}}({\mathcal{F}}_1^{\mathrm{dis}})]$ , and thus also in $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}})]$ . So there is a connected component of $G[{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}})]$ with size at least $\frac{r}{K + 1}-2$ , which corresponds to the third condition.

Next we prove Theorem6.

Theorem 6. For all integer $\Delta \geq 5$ , for all $q$ sufficiently large and any $\beta \leq \beta _c$ , w.h.p. over $G\sim{\mathcal{G}}_{n,\Delta }$ , $G$ has WSM within the disordered phase at a radius $r$ which satisfies $r\leq \tfrac{1}{3}\log _{\Delta -1}n$ .

Proof. By Lemmas7 and 16, there exist $K=K(\Delta )\gt 0$ and $\delta =\delta (\Delta ) \in (0,1/2)$ such that, w.h.p., $G\sim{\mathcal{G}}_{\Delta, n}$ is in ${\mathcal{G}}_{\Delta, \delta, K}$ . Let $\beta _1(q) \;:\!=\;\log (q^{2.1/\Delta }+1)$ . Let $\zeta$ be the constant from Lemma3 and let $C'$ be constant implicit in Equation 2 so that this equation guarantees $\pi _G^{\mathrm{dis}}\Big (|{\mathrm{In}}({\mathcal{F}})|\geq \zeta |E|\Big ) \leq e^{-C' n}$ .

Let $q_0$ and $n_0$ be large enough that

  • $\beta _c(q_0) \leq \beta _1(q_0)$ .

  • Lemmas3 and 48 apply. Lemma26 applies with $C = \frac{1}{3(2+K)\log (\Delta -1)}$ .

  • $(\eta - \zeta ) n_0 \geq \Delta n_0^{1/3}$ and $C'n_0 \geq \log (100\Delta n_0)$ and $n_0\geq (50\Delta )^2$ .

  • $\frac{\log _{\Delta -1} n_0}{3(K+1)}-2\geq \frac{\log _{\Delta -1} n_0}{3(K+2)}$ ,

Now fix $q\geq q_0$ and $\beta \leq \beta _c(q)$ .

Consider any $n\geq n_0$ . Consider an $n$ -vertex graph $G=(V,E)\in{\mathcal{G}}_{\Delta, \delta, K}$ . Let $r\;:\!=\; \frac 13 \log _{\Delta -1}(n)$ . Fix $v\in V$ . Let $r_1$ be as in Lemma48. Use the process to generate ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ . We wish to show that for every edge $e$ incident to $v$ ,

\begin{equation*} \Vert {\pi _{{\mathcal {B}}_{r_1}^-(v)}(e \mapsto \cdot ) - \pi ^{\mathrm {dis}}(e\mapsto \cdot ) }\Vert _{\mathrm {TV}} \leq 1/(100|E|) \end{equation*}

By Lemma48, if ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ do not agree on all edges incident to $v$ , then $|{\mathrm{In}}({\mathcal{F}}^{\mathrm{dis}})|\geq \eta |E| - |E(B_{r_1}(v))|$ or ${\mathcal{F}}^{\mathrm{ord}}$ contains a polymer of size at least $\frac{r}{K+1}-2\geq \frac{\log n}{3(K+2)\log (\Delta -1)} = C \log n$ . Thus the probability that ${\mathcal{F}}^{\mathrm{ord}}$ and ${\mathcal{F}}^+$ do not agree on all edges incident to $v$ is at most the sum of the probabilities of these two events.

Since $(\eta - \zeta ) n \geq \Delta n^{1/3}$ , $\zeta |E| \leq \eta |E| - |E(B_{r_1}(v))|$ , so Lemma3 guarantees that the probability of the first event is at most $e^{-C'n}\leq \frac{1}{100\Delta n} = \frac{1}{200|E|}$ .

By Lemma26, the probability of the second event is at most $n^{-3/2} / 2\leq \frac{1}{100\Delta n} = \frac{1}{200|E|}$ . The result follows by summing these probabilities.

7. Proof of Theorems12 and 13

In this section, we provide the last ingredients that were used in the proof of Theorem1 by showing the $O(n\log n)$ bounds in Theorems12 and 13. The overall argument is quite close to what was presented in Section 2.4 for the proof of Theorem11; the only extra argument required is a slightly more refined estimate from the distance to stationarity using log-Sobolev constants to save an $O(\log n)$ factor (a similar argument appears in [Reference Blanca and Gheissari4]). We present first the relevant tools in Sections 7.1 and 7.2, and then finish the proofs in Section 7.3.

7.1 The log-sobolev constant and mixing time inequality

Let $\mu$ be a distribution supported on a set $\Omega$ . For a function $g:\Omega \rightarrow \mathbb{R}_{\geq 0}$ , let

\begin{equation*}\mathrm {Ent}_{\mu }[g]\;:\!=\;\mathbb {E}_{\mu }[g \log g]-\mathbb {E}_{\mu }[g]\log \mathbb {E}_{\mu }[\log g],\end{equation*}

with the convention $0\log 0\;:\!=\;0$ . Moreover, for a reversible Markov chain $(X_t)_{t\geq 0}$ on $\Omega$ with transition matrix $P\in \mathbb{R}^{\Omega \times \Omega }_{\geq 0}$ , let $\mathcal{E}(g,g)\;:\!=\;\tfrac{1}{2}\sum _{\omega, \omega '}\mu (\omega )P(\omega, \omega ')(g(\omega )-g(\omega '))^2$ be the Dirichlet form for $P$ . The standard log-Sobolev constant of the chain is then defined as

\begin{equation*}\alpha (P)\;:\!=\;\min _{\substack {g:\Omega \rightarrow \mathbb {R}_{\geq 0};\\[5pt] \mathrm {Ent}_{\mu }[g]\neq 0}}\frac {\mathcal {E}(\sqrt {g},\sqrt {g})}{\mathrm {Ent}_\mu [g]}.\end{equation*}

The following well-known connection between the log-Sobolev constant and the distance from stationarity can be found, for example, in [Reference Blanca and Gheissari4, Fact 6.1]. Let $\mu _{\min }=\min _{\omega \in \Omega }. \mu (\omega )$ . Then, for any $\gamma \lt \alpha (P)$ , it holds that

(8) \begin{equation} \mbox{$\max _{X_0}$}\, \mathrm{dist}_{\mathrm{TV}}(X_t,\mu )\leq \mathrm{e}^{-\gamma t/2}\big (\log \tfrac{1}{\mu _{\min }}\big )^{1/2}. \end{equation}

7.2 Log-sobolev constant for free/wired tree-like neighbourhoods

Here we briefly discuss how to bound the log-Sobolev for the wired RC dynamics for tree-like neighbourhoods for integer $q\gt 1$ (cf. Remark10) using the results from [Reference Blanca, Chen, Štefankovič and Vigoda2] on the tree.

Lemma 49. Let $q\geq 2$ and $\Delta \geq 3$ be integers, and $K,\beta \gt 0$ be reals. There exists $\tilde{C}\gt 0$ such that the following holds for every $\Delta$ -regular graph $G=(V,E)$ and any integer $r\geq 1$ .

Suppose that $\rho \in V$ is such that $G[B_r(\rho )]$ is $K$ -treelike. Then, with $n=|B_r(\rho )|$ , the log-Sobolev constant for the wired RC dynamics on $B_r(\rho )$ is $\geq \tilde{C}/n$ .

Proof. Let $T=T_{\Delta }(\rho )$ denote the $\Delta$ -regular tree rooted at $\rho$ . Then, for integer $q\geq 2$ and $\beta \gt 0$ , it is shown in [Reference Blanca, Chen, Štefankovič and Vigoda2, Proof of Lemma 32] that there exists $C\gt 0$ such that the log-Sobolev constant for the wired RC dynamics on $T[B_r(\rho )]$ is $\geq C/N$ where $N=|B_r(\rho )|$ .

To translate this into a lower bound on the log-Sobolev constant of $G[B_r(\rho )]$ (which is $K$ -treelike) as in the statement of the lemma, we can just use the argument in [Reference Gheissari and Sinclair15, Proof of Lemma 4.4]. There, they show, for the Ising model with all plus boundary condition, it via a graph decomposition argument that there exists a constant $C'=C(q,\beta, K,\Delta )$ independent of $r$ such that $\alpha \big (P_{G[B_r(\rho )]}\big )\geq C\alpha \big (P_{T[B_r(\rho )]}\big )$ whenever $G[B_r(\rho )]$ is $K$ -treelike. The details of the graph decomposition do not depend on the Ising model, so the exact same strategy yields the analogue for the RC model.

The analogous result for the free boundary is available from [Reference Blanca and Gheissari4] for all $q\geq 1$ , we state here the following more precise version of Lemma9.

Lemma 50 ([Reference Blanca and Gheissari4, Lemma 6.5]). Let $\Delta \geq 3$ be an integer, and $q,K\gt 1$ , $\beta \gt 0$ be reals. There exists $C\gt 0$ such that the following holds for any $\Delta$ -regular graph $G$ and integer $r\geq 1$ .

Suppose that $\rho \in V$ is such that $G[B_r(\rho )]$ is $K$ -treelike. Then, with $n=|B_r(\rho )|$ , the log-Sobolev constant for the free RC dynamics on $B_r(\rho )$ is $\geq C/n$ .

7.3 Completing the proof of Theorems 12 and 13

Proof of Theorems 12 and 13. We first show Theorem12. Consider therefore $\Delta \geq 5$ , and let $q$ be sufficiently large so that both Lemma3 and Theorem6 apply, and suppose that $\beta \leq \beta _c$ .

Consider $G=(V,E)\sim{\mathcal{G}}_{n,\Delta }$ with $n=|V|$ and $m=|E|$ . By Lemma7, we can assume that $G$ is locally $K$ -treelike. By Lemma3, $\pi _G^{\mathrm{dis}}\big (|{\mathrm{In}}({\mathcal{F}})|\geq \zeta |E|\big )=\mathrm{e}^{-\Omega (n)}$ . By Theorem6, $G$ has WSM within the disordered phase at radius $r$ for some $r\leq \tfrac{1}{3}\log _{\Delta -1}n$ .

We will consider the RC dynamics $(X_t)_{t\geq 0}$ with $X_0$ being the all-out configuration on the edges. We will also consider the ‘disordered’ RC dynamics $(\hat{X}_t)_{t\geq 0}$ with $\hat{X}_0\sim \pi ^{\mathrm{dis}}_G$ where we reject moves of the chain that lead to configurations outside of $\Omega ^{\mathrm{dis}}$ . Note that $\hat{X}_t\sim \pi ^{\mathrm{dis}}_G$ for all $t\geq 0$ . We will show how to couple these two dynamics in $O(n \log n)$ time. The result will follow by showing that there is a coupling between $(X_t)_{t\geq 0}$ and $(\hat{X}_t)_{t\geq 0}$ such that, for $T=O(n\log n)$ , it holds that

(9) \begin{equation} {\mathbb{P}}(X_T\neq \hat{X}_T)\leq 1/4. \end{equation}

Analogously to the ordered case we use the monotone coupling, where at every step $t$ , the two chains choose the same edge $e_t$ to update and use the same uniform number $U_t\in [0,1]$ to decide whether to include $e_t$ in each of $X_{t+1},\hat{X}_{t+1}$ . For $t\geq 0$ , let $\mathcal{E}_t$ be the event that ${\mathrm{In}}(\hat{X}_t)\leq \zeta |E|$ and let $\mathcal{E}_{\lt t}\;:\!=\;\bigcap _{t'=0,\ldots, t-1} \mathcal{E}_{t'}$ . From Lemma3 we have that $\pi ^{\mathrm{dis}}_G(\mathcal{E}_{\lt t})\geq 1-t\mathrm{e}^{-\Omega (n)}$ . Using again the monotonicity of the model for $q\geq 1$ , under the monotone coupling, for all $t\geq 0$ such that $\mathcal{E}_{\lt t}$ holds (and hence no reject move has happened in $\hat{X}_t$ so far), we have that $X_t\leq \hat{X}_t$ (i.e., ${\mathrm{In}}(X_t)\subseteq{\mathrm{In}}(\hat{X}_t)$ ). We next proceed to bound the terms in the upper bound

(4) \begin{align}{\mathbb{P}}\big (X_{t}\neq \hat{X}_t\big )\leq \sum _{e}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\big )\leq m{\mathbb{P}}\big (\overline{ \mathcal{E}_{\lt t}}\big )+\sum _{e}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\mid \mathcal{E}_{\lt t}\big ). \end{align}

So, fix an arbitrary edge $e$ incident to some vertex $v$ , and let $(X^v_t)$ be the free RC dynamics on $G[B_r(v)]$ . We couple the evolution of $(X_t^v)$ with that of $(X_t)$ and $(\hat{X}_t)$ using the monotone coupling analogously to the ordered case, where in $X_t^v$ we ignore updates of edges outside the ball $G[B_r(v)]$ ). We have $X_t^v\leq X_t$ for all $t\geq 0$ , and hence, conditioned on $\mathcal{E}_{\lt t}$ , we have that $X^v_t\leq X_t\leq \hat{X}_t$ . Proceeding as in the proof of Theorem11, we therefore get the following analogue of (5):

(10) \begin{equation} \begin{aligned}{\mathbb{P}}\big (X_t(e)\neq \hat{X}_t(e)\mid \mathcal{E}_{\lt t}\big )\leq 4{\mathbb{P}}\big (\overline{ \mathcal{E}_{\lt t}}\big )+\big |{\mathbb{P}}(X^v_t(e)&=0)-\pi _{{\mathcal{B}}^-_r(v)}(e\mapsto 0)\big |\\[5pt] &+\big |\pi _{{\mathcal{B}}^-_r(v)}(e\mapsto 0)-\pi ^{\mathrm{dis}}_G(e\mapsto 0)\big |. \end{aligned} \end{equation}

Since $G$ has WSM within the disordered phase at radius $r$ , we have that

(11) \begin{equation} \big |\pi _{{\mathcal{B}}^-_r(v)}(e\mapsto 0)-\pi ^{\mathrm{dis}}_G(e\mapsto 0)\big |\leq 1/(100 m). \end{equation}

Let $C\gt 0$ be the constant in Lemma50 and set $N_v=|E(B_r(v))|\leq \Delta ^{r+1}$ . By Chernoff bounds, for $T=\Theta (n\log n)$ , we get at least $t_v\;:\!=\;\tfrac{40}{C}N_v\log n$ edge updates within the ball $B_r(v)$ with probability $1-\exp (-n^{\Omega (1)})$ . Since $r\leq \tfrac{1}{3}\log _{\Delta -1} n$ , $G[B_r(v)]$ is $K$ -treelike, so from Lemma50 and applying the log-Sobolev inequality (8) (with $\gamma =C/(2N_v)$ , $t=t_v$ and $\min _{{\mathcal{F}}}\pi _{{\mathcal{B}}^-_r(v)}({\mathcal{F}})\geq (2q \mathrm{e}^{\beta })^{-N_v}$ ), we have

(12) \begin{equation} \big |{\mathbb{P}}(X^v_T(e)=0)-\pi _{{\mathcal{B}}^-_r(v)}(e\mapsto 0)\big |\leq \exp (\!-n^{\Omega (1)})+ \mathrm{e}^{-C t_v/(4N_v) }\big (\log \tfrac{1}{(2q \mathrm{e}^{\beta })^{-N_v}} \big )^{1/2}\leq 1/m^3, \end{equation}

where the last inequality holds for all $m=\Delta n=\Omega (1)$ . Plugging (11) and (12) into (10) for $t=T$ , and then back into (4), we get (9), that is, ${\mathbb{P}}(X_T\neq \hat{X}_T)\leq 5m T \mathrm{e}^{-\Omega (n)}+m/m^3+1/100\leq 1/5$ . This finishes the proof of Theorem11.

For the proof of Theorem13 (integer $q$ and $\beta \geq \beta _c$ ), the argument is completely analogous to what was just presented for $\beta \leq \beta _c$ , using now the log-Sobolev bound of Lemma8 to get the analogue of (12) (see the proof in Section 2.4 and the inequality (7)).

Appendix A. Mixing on tree-like graphs for the wired RC dynamics

In this section, we prove Lemma9 which we restate here for convenience.

Lemma 9. Let $\Delta \geq 3$ be an integer, and $q,K\gt 1$ , $\beta \gt 0$ be reals. There exists $\hat{C}\gt 0$ such that the following holds for any $\Delta$ -regular graph $G$ and integer $r\geq 1$ .

Suppose that $\rho \in V$ is such that $G[B_r(\rho )]$ is $K$ -treelike. Then, with $n=|B_r(\rho )|$ , the mixing time of the free RC dynamics on $B_r(\rho )$ is $\leq \hat{C}n\log n$ .

Proof. We use a canonical paths argument, following largely [Reference Mossel and Sly25, Proof of Lemma8]. Let $G=(V,E)$ be a $\Delta$ -regular graph and $\rho \in V$ be such that $G[B_r(\rho )]$ is $K$ -treelike. Let $n=|B_r(\rho )|$ , $m=|E(B_r(\rho )|$ . Let $\hat{\pi }=\hat{\pi }_{B_r(\rho )}$ be the stationary distribution of the wired RC dynamics on $B_r(\rho )$ .

Consider a BFS tree $T$ for $G[B_r(\rho )]$ starting from $\rho$ , let $L$ be the set of the leaves of $T$ and let $\mathcal{E}=E(G[B_r(\rho )])\backslash E(T)$ be the set of excess edges. Note that the height of $T$ is at most $r$ , the set $S_r(\rho )$ of vertices at distance $r$ from $\rho$ satisfies $S_r(\rho )\subseteq L$ and there are at most $2K$ leaves of $T$ at depth less than $r$ (by the $\Delta$ -regularity of $G$ , each such leaf must be incident to an excess edge of $G$ ), so $|L\backslash S_r(\rho )|\leq 2K$ . For a vertex $u$ in $B_r(\rho )$ , we let $T_u$ be the subtree of $T$ rooted at $u$ , and $\mathrm{Anc}(u)$ be the ancestors of $u$ in $T$ (i.e., the vertices on the path from the root $\rho$ to $u$ , including $u$ ).

Consider next a depth-first-search traversal of the tree $T$ starting from $\rho$ , and let $v_1,\ldots, v_n$ be the order in which the vertices of $T$ were first visited. We write $v_i\lt v_j$ whenever $i\lt j$ ; moreover, for a subset of vertices $S$ and a vertex $u$ , we write $S\lt u$ to denote that for each $w\in S$ it holds that $w\lt u$ . Since the tree $T$ has depth $r$ , at any stage of the DFS traversal there are at most $r$ vertices which are visited but not fully explored (note that a leaf that is visited is automatically explored). Since the degree of any vertex is $\leq \Delta$ , it follows that

  1. (a) for $i=1,\ldots, n$ , there are at most $\Delta r$ tree edges with one endpoint in $\{v_1,\ldots, v_i\}$ and the other in $\{v_{i+1},\ldots, v_m\}$ (from the DFS traversal).

  2. (b) for an arbitrary vertex $u$ in $T$ and any vertex $w$ which is not an ancestor of $u$ (i.e., $w\notin \mathrm{Anc}(u)$ ), we have from the DFS traversal either that $V(T_w)\lt u$ or $V(T_w) \gt u$ .

We order the edges $e_1,\ldots, e_m$ in $G[{\mathcal{B}}_r(\rho )]$ (including excess edges) in lexicographic order in terms of the order of the vertices. So edges with an endpoint $v_1$ come first (ordered by the other endpoint) and so on.

For two configurations ${\mathcal{F}},{\mathcal{F}}':E\rightarrow \{0,1\}$ that differ on the edges $e_{i_1},\ldots, e_{i_k}$ (with $i_1\lt \cdots \lt i_k$ ), define the path $\mathrm{Path}({\mathcal{F}},{\mathcal{F}}')$ of configurations $\mathrm{Path}({\mathcal{F}},{\mathcal{F}}')={\mathcal{F}}_0\rightarrow{\mathcal{F}}_1\cdots \rightarrow{\mathcal{F}}_{k}$ by defining, for $j=0,\ldots, k$ , ${\mathcal{F}}_{j}$ to agree with $\mathcal{F}$ on the edges $\{e_{i_{j+1}},\ldots, e_{i_k}\}$ and to agree with ${\mathcal{F}}'$ on the edges $\{e_{i_1},\ldots, e_{i_j}\}$ . Note that ${\mathcal{F}}_0={\mathcal{F}},{\mathcal{F}}_k={\mathcal{F}}'$ and both ${\mathcal{F}},{\mathcal{F}}'$ and hence ${\mathcal{F}}_j$ as well agree on $E\backslash \{e_{i_1},\ldots, e_{i_k}\})$ . Then, the relaxation time of the chain, see for example [Reference Jerrum23, Chapter 5], is bounded by

\begin{equation*}\tau \leq m\max _{({\mathcal {F}},{\mathcal {F}}')}\sum _{\substack {{\mathcal {A}},{\mathcal {A}}':E\rightarrow \{0,1\};\\[5pt] ({\mathcal {F}}\rightarrow {\mathcal {F}}')\in \mathrm {Path}({\mathcal {A}},{\mathcal {A}}')}}\frac {\hat {\pi }({\mathcal {A}})\hat {\pi }({\mathcal {A}}')}{\hat {\pi }({\mathcal {F}}){\mathbb {P}}({\mathcal {F}},{\mathcal {F}}')},\end{equation*}

where the maximum is over all pairs of configurations ${\mathcal{F}},{\mathcal{F}}'$ that differ on a single edge and the summation is over all pairs of configurations ${\mathcal{A}},{\mathcal{A}}'$ whose corresponding $\mathrm{Path}({\mathcal{A}},{\mathcal{A}}')$ includes the transition ${\mathcal{F}}\rightarrow{\mathcal{F}}'$ as part of the path.

Consider an arbitrary pair of configurations $\kappa =({\mathcal{F}},{\mathcal{F}}')$ that differ at a single edge $f=\{v_j, v_{j'}\}$ for some $j\lt j'$ and suppose further that $f=e_\ell$ . For configurations ${\mathcal{A}},{\mathcal{A}}'$ such that $({\mathcal{F}}\rightarrow{\mathcal{F}}')\in \mathrm{Path}({\mathcal{A}},{\mathcal{A}}')$ define the configuration $g_{\kappa }({\mathcal{A}},{\mathcal{A}}')$ so that

(13) \begin{equation} \mbox{$g_{\kappa }({\mathcal{A}},{\mathcal{A}}')$ agrees with ${\mathcal{A}}$ on $\{e_1,\ldots, e_{\ell -1}\}$ and with ${\mathcal{A}}'$ on $\{e_\ell, \ldots, e_{m}\}$.} \end{equation}

From the fact that $({\mathcal{F}}\rightarrow{\mathcal{F}}')\in \mathrm{Path}({\mathcal{A}},{\mathcal{A}}')$ , we also have that:

(14) \begin{equation} \mbox{${\mathcal{A}}$ agrees with ${\mathcal{F}}$ on $\{e_{\ell },\ldots, e_{m}\}$},\ \ \mbox{${\mathcal{A}}'$ agrees with ${\mathcal{F}}'$ on $\{e_{1},\ldots, e_{\ell -1}\}$.} \end{equation}

It follows that the map $g_{\kappa }(\cdot, \cdot )$ is injective, that is, given its value ${\mathcal{F}}^*$ and the configurations ${\mathcal{F}},{\mathcal{F}}'$ , there is a unique pair $({\mathcal{A}},{\mathcal{A}}')$ such that ${\mathcal{F}}^*=g_{\kappa }({\mathcal{A}},{\mathcal{A}}')$ and $({\mathcal{F}}\rightarrow{\mathcal{F}}')\in \mathrm{Path}({\mathcal{A}},{\mathcal{A}}')$ . For a configuration $\mathcal{X}\;:\; E(B_r(\rho ))\rightarrow \{0,1\}$ , let $\hat{c}(\mathcal{X})$ denote the number of components in $(B_r(\rho ),{\mathrm{In}}(X))$ that do not include any of the vertices in $S_r(\rho )$ , cf. Footnote 5. Let $s\;:\!=\; \Delta r + K$ and $t \;:\!=\; 4 \Delta r + 12 K$ . The main step in the proof is to show that:

(15) \begin{equation} |{\mathrm{In}}({\mathcal{A}})|+|{\mathrm{In}}({\mathcal{A}}')|-|{\mathrm{In}}({\mathcal{F}})|-|{\mathrm{In}}(g_\kappa ({\mathcal{A}},{\mathcal{A}}'))|\leq s, \end{equation}
(16) \begin{equation} |\hat{c}({\mathcal{A}})+\hat{c}({\mathcal{A}}')-\hat{c}({\mathcal{F}})-\hat{c}(g_{\kappa }({\mathcal{A}},{\mathcal{A}}'))|\leq t. \end{equation}

Assuming these for the moment, we have that

\begin{equation*} \frac {\hat {\pi }({\mathcal {A}})\hat {\pi }({\mathcal {A}}')}{\hat {\pi }({\mathcal {F}})\hat {\pi }\big (g_{\kappa }({\mathcal {A}},{\mathcal {A}}')\big )}\leq q^{t}\mathrm {e}^{\beta s}.\end{equation*}

Recall that $p=1-\mathrm{e}^{-\beta }$ and that for $q\gt 1$ it holds that $\hat{p}\in (p/q,p)$ . Therefore, for the RC dynamics, we have that ${\mathbb{P}}({\mathcal{F}},{\mathcal{F}}')\geq \frac{\min \{p/q,(1-p)\}}{m} \geq (1-\mathrm{e}^{-\beta })/(q m \mathrm{e}^\beta )$ using that $\min \{p/q,(1-p)\}\geq \mathrm{e}^{-\beta }/q$ . Using this and the injectivity of $g$ , we can bound $\tau$ by

\begin{equation*}\tau \leq \frac {m^2 q^{t+1}\mathrm {e}^{\beta (s+1)}} {(1-\mathrm {e}^{-\beta })} \sum _{\substack {{\mathcal {A}},{\mathcal {A}}':E\rightarrow \{0,1\};\\[5pt] ({\mathcal {F}}\rightarrow {\mathcal {F}}')\in \mathrm {Path}({\mathcal {A}},{\mathcal {A}}')}} \hat {\pi }\big (g_{\kappa }({\mathcal {A}},{\mathcal {A}}')\big )\leq \frac {m^2 q^{t+1}\mathrm {e}^{\beta (s+1)}} {(1-\mathrm {e}^{-\beta })}.\end{equation*}

The mixing time is at most $ \tau (1 + \tfrac 12 \log (\min _{\mathcal{X}} \hat{\pi }(\mathcal{X})^{-1}))$ Since $\min _{\mathcal{X}}\hat{\pi }(\mathcal{X})\geq 2^{-n} q^{-n} e^{-\beta m}$ and $m\leq \Delta n$ , the mixing time is at most $n^3 \Delta ^2 q^{t+1}\mathrm{e}^{\beta (s+1)}\tfrac{1+ \log (q) + \beta \Delta }{1-\mathrm{e}^{-\beta }}$ so the statement for the mixing time follows by taking $ \hat{C}\;:\!=\; 2 \Delta ^3 (q \mathrm{e}^{\beta })^{12K+3}/{(1-\mathrm{e}^{-\beta })}$ .

It remains to show (15) and (16). For convenience let ${\mathcal{F}}^*=g_{\kappa }({\mathcal{A}},{\mathcal{A}}')$ . We will shortly prove the following facts for an arbitrary edge $e=\{u,w\}$ :

  1. 1. if ${\mathcal{A}}(e)={\mathcal{A}}'(e)$ , then ${\mathcal{A}}(e)={\mathcal{A}}'(e)={\mathcal{F}}(e)={\mathcal{F}}^*(e)$ .

  2. 2. if ${\mathcal{A}}(e)\neq{\mathcal{A}}'(e)$ , then $\big ({\mathcal{A}}(e),{\mathcal{A}}'(e)\big )= \big ({\mathcal{F}}^*(e),{\mathcal{F}}(e)\big )$ when $u,w\leq v_j$ and $\big ({\mathcal{A}}(e),{\mathcal{A}}'(e)\big )= \big ({\mathcal{F}}(e),{\mathcal{F}}^*(e)\big )$ when $u,w\gt v_j$ .

Indeed, for Item 1, consider $e$ with ${\mathcal{A}}(e) ={\mathcal{A}}'(e)$ and suppose first that ${\mathcal{A}}(e)={\mathcal{A}}'(e)=1$ . Then ${\mathcal{F}}(e)={\mathcal{F}}'(e)=1$ (since by the construction of $\mathrm{Path}({\mathcal{A}},{\mathcal{A}}')$ , all configurations in it agree on the edges where ${\mathcal{A}},{\mathcal{A}}'$ agree). Moreover, ${\mathcal{F}}^*(e)=1$ by construction (since ${\mathcal{F}}^*=g_{\kappa }({\mathcal{A}},{\mathcal{A}}')$ agrees with ${\mathcal{A}},{\mathcal{A}}'$ on the edges where the latter agree). The proof for the second case ${\mathcal{A}}(e)={\mathcal{A}}'(e)=0$ is analogous. Item 2 follows from (13) and (14) after observing that $u,w\leq v_j$ implies that $e=\{u,w\}\lt e_\ell$ , while $u,w\gt v_j$ gives that $e\gt e_\ell$ .

Let $W_j$ be the tree edges which are incident to a vertex in $\mathrm{Anc}(v_j)$ , so that $|W_j|\leq \Delta r$ . From Item (1), these are the only tree edges that can have one endpoint in $\{v_1,\ldots, v_j\}$ and the other among $\{v_{j+1},\ldots, v_n\}$ , along with any non-tree edges (i.e., excess edges $\mathcal{E}$ ). Therefore, from Items 1 and 2 above, it follows that for an edge $e\in E\backslash (W_j \cup \mathcal{E})$ it holds that

(17) \begin{equation} {\mathcal{A}}(e)+{\mathcal{A}}'(e)={\mathcal{F}}(e)+{\mathcal{F}}^*(e), \end{equation}

which establishes (15) since $|W_j \cup \mathcal{E}|\leq \Delta r+K$ .

To show (16), fix an arbitrary root-to-leaf path $P$ passing through $v_j$ , and denote by $E_P$ be the edges of the path. Now, consider the slightly ‘tweaked’ configurations $\hat{{\mathcal{A}}},\hat{{\mathcal{A}}}',\hat{{\mathcal{F}}},\hat{{\mathcal{F}}}^*$ obtained from ${\mathcal{A}},{\mathcal{A}}',{\mathcal{F}},{\mathcal{F}}^*$ , respectively, by setting, for each $\mathcal{X}\in \{{\mathcal{A}},{\mathcal{A}}',{\mathcal{F}},{\mathcal{F}}^*\}$ , $\hat{\mathcal{X}}(e)=1$ for $e\in E_P\cup W_j$ , $\hat{\mathcal{X}}(e)=0$ for $e\in \mathcal{E}$ and $\hat{\mathcal{X}}(e)=\mathcal{X}(e)$ for $e\notin E_P\cup W_j\cup \mathcal{E}$ . Note that, for each $\mathcal{X}\in \{{\mathcal{A}},{\mathcal{A}}',{\mathcal{F}},{\mathcal{F}}^*\}$ , we have $\big ||\hat{\mathcal{X}}|-|\mathcal{X}|\big |\leq \Delta r+K$ , so $\big ||\hat{c}(\hat{\mathcal{X}})|-|\hat{c}(\mathcal{X})|\big |\leq \Delta r+K$ as well. So, to prove (16), it suffices to show

(18) \begin{equation} |\hat{c}(\hat{\mathcal{A}})+\hat{c}(\hat{\mathcal{A}}')-\hat{c}(\hat{\mathcal{F}})-\hat{c}(\hat{\mathcal{F}}^*)|\leq 8K. \end{equation}

To show this, first note that $\hat{{\mathcal{A}}},\hat{{\mathcal{A}}}',\hat{{\mathcal{F}}},\hat{{\mathcal{F}}}^*$ still satisfy Items 1 and 2 above (replacing $\mathcal{A}$ with $\hat{\mathcal{A}}$ , and so on) since the only changes are for edges in $E_P\cup W_j\cup \mathcal{E}$ on which the assignments of the tweaked configurations are identical (and hence fall under Item 1). In fact, Items 1 and 2 now apply to all edges for the tweaked configurations (previously, edges in $W_j\cup \mathcal{E}$ were potentially not covered in Item 2, but now are covered since they fall under Item 1). So, the analogue of (17) for $\hat{{\mathcal{A}}},\hat{{\mathcal{A}}}',\hat{{\mathcal{F}}},\hat{{\mathcal{F}}}^*$ holds for every edge of $G[B_r(\rho )]$ , yielding that

(19) \begin{equation} |{\mathrm{In}}(\hat{\mathcal{A}})|+|{\mathrm{In}}(\hat{\mathcal{A}}')|=|{\mathrm{In}}(\hat{\mathcal{F}})|+|{\mathrm{In}}(\hat{\mathcal{F}}^*)|, \end{equation}

and, further, for any edge $e=\{u,w\}$ we have that

(20) \begin{equation} \big (\hat{\mathcal{A}}(e),\hat{\mathcal{A}}'(e)\big )=\begin{cases} \big (\hat{\mathcal{F}}^*(e),\hat{\mathcal{F}}(e)\big ) & \mbox{ if } u,w\leq v_j\\[5pt] \big (\hat{\mathcal{F}}(e),\hat{\mathcal{F}}^*(e)\big ),& \mbox{ otherwise. } \end{cases} \end{equation}

At this stage, it will be convenient to consider the graph $T^*$ obtained from $G[B_r(\rho )]$ by identifying all the leaves of $T$ into a single vertex $v^*$ (recall that the leaves of $T$ consist of the vertices $S_r(\rho )$ , which are at depth $r$ , together with at most $2K$ leaves that are at depth $\lt r$ ). For $\mathcal{X}\in \{\hat{\mathcal{A}},\hat{\mathcal{A}}',\hat{\mathcal{F}},\hat{\mathcal{F}}^*\}$ , let $c^*(\mathcal{X})$ be the number of components in the graph $T^*[{\mathrm{In}}(\mathcal{X})]$ that do not include $v^*$ and let ${\mathcal{C}}^*(\mathcal{X})=(V^*(\mathcal{X}),E^*(\mathcal{X}))$ be the connected component of $v^*$ . Note that any connected component of $(B_r(\rho ),{\mathrm{In}}(\mathcal{X}))$ that does not include a vertex from $S_r(\rho )$ is not connected to $v^*$ in $T^*[{\mathrm{In}}(\mathcal{X})]$ , unless it contains one of the $2K$ leaves at depth $\lt r$ . It follows that $|c^*(\mathcal{X})-\hat{c}(\mathcal{X})|\leq 2K$ , so to prove (18) it suffices to show

(21) \begin{equation} c^*(\hat{\mathcal{A}})+c^*(\hat{\mathcal{A}}')= c^*(\hat{\mathcal{F}})+ c^*(\hat{\mathcal{F}}^*). \end{equation}

For $\mathcal{X}\in \{\hat{\mathcal{A}},\hat{\mathcal{A}}',\hat{\mathcal{F}},\hat{\mathcal{F}}^*\}$ , note that a connected component $F$ of $\mathcal{X}$ that does not include $v^*$ is a subgraph of $T$ (using that $\mathcal{X}(e)=0$ for $e\in \mathcal{E}$ ) and hence a tree, so $1=|V(F)|-|E(F)|$ . Summing over all such components, we obtain that

\begin{equation*}c^*(\mathcal {X})=|V(T^*)|-|V^*(\mathcal {X})|-|{\mathrm {In}}(\mathcal {X})|+|E^*(\mathcal {X})|.\end{equation*}

Therefore, using (19), to finish the proof of (21) it suffices to show that for an arbitrary vertex $v$ and an arbitrary edge $e$ it holds that

(22) \begin{equation} \textbf{1}\{v\in V^*(\hat{\mathcal{A}})\}+\textbf{1}\{v\in V^*(\hat{\mathcal{A}}')\}=\textbf{1}\{v\in V^*(\hat{\mathcal{F}})\}+\textbf{1}\{v\in V^*(\hat{\mathcal{F}}^*)\}. \end{equation}
(23) \begin{equation} \textbf{1}\{e\in E^*(\hat{\mathcal{A}})\}+\textbf{1}\{e\in E^*(\hat{\mathcal{A}}')\}=\textbf{1}\{e\in E^*(\hat{\mathcal{F}})\}+\textbf{1}\{e\in E^*(\hat{\mathcal{F}}^*)\}. \end{equation}

Clearly (22) is true for any $v\in \mathrm{Anc}(v_j)$ since $v\in{\mathcal{C}}(\mathcal{X})$ for each $\mathcal{X}\in \{\hat{\mathcal{A}},\hat{\mathcal{A}}',\hat{\mathcal{F}},\hat{\mathcal{F}}^*\}$ by construction of the tweaked configurations. So, consider $v\notin \mathrm{Anc}(v_j)$ . First, suppose $v\lt v_j$ . Let $u$ be the ancestor of $v$ whose parent $w$ is an ancestor of $v_j$ . The key point is that for each $\mathcal{X}\in \{\hat{\mathcal{A}},\hat{\mathcal{A}}',\hat{\mathcal{F}},\hat{\mathcal{F}}^*\}$ we have $w\in V^*(\mathcal{X})$ (since $w$ is an ancestor of $v_j$ ) and the edge $e=\{u,w\}$ belongs to $W_j$ , so $\hat{\mathcal{X}}(e)=1$ . Therefore whether $v \in V^*(\mathcal{X})$ is determined by the set of those paths contained in the subtree $T_u$ which start from $v$ and end either at $u$ or a leaf of $T_u$ . Denote by $\mathcal{P}_{v}$ the set of all such paths. Since $v\lt v_j$ , we have that $u\lt v_j$ and that $u$ is not ancestor of $v_j$ , so from Item 2 it holds that $V(T_u)\lt v_j$ . It follows from (20) that for every edge $e\in E(T_u)$ we have $(\hat{\mathcal{A}}(e),\hat{\mathcal{A}}'(e))=({\mathcal{F}}^*(e),\hat{\mathcal{F}}^*(e))$ , so

\begin{equation*}v\in V^*(\hat {\mathcal {A}}) \Leftrightarrow \, \mbox {$\exists P'\in \mathcal {P}_{v}\;:\; E(P')\subseteq {\mathrm {In}}(\hat {\mathcal {A}})$} \Leftrightarrow \,\mbox {$\exists P'\in \mathcal {P}_{v}\:\; E(P')\subseteq {\mathrm {In}}(\hat {\mathcal {F}}^*)$} \Leftrightarrow v\in V^*(\hat {\mathcal {F}}^*),\end{equation*}

and similarly $ v\in V^*(\hat{\mathcal{A}}')\Leftrightarrow v\in V^*(\hat{\mathcal{F}})$ , proving (22) for $v\lt v_j$ . An analogous argument applies when $v\gt v_j$ , then $V(T_u)\gt v_j$ and, using (20) again, $v\in V^*(\hat{\mathcal{A}})\Leftrightarrow v\in V^*(\hat{\mathcal{F}})$ , $ v\in V^*(\hat{\mathcal{A}}')\Leftrightarrow v\in V^*(\hat{\mathcal{F}}^*)$ , finishing the proof of (22). The proof of (23) is very similar, after noting that for $\mathcal{X}\in \{\hat{\mathcal{A}},\hat{\mathcal{A}}',\hat{\mathcal{F}},\hat{\mathcal{F}}^*\}$ and an edge $e=\{u,w\}$ it holds that $e\in E^*(\mathcal{X})$ iff $ \mathcal{X}(e)=1$ and $u,w\in V^*(\mathcal{X})$ .

This completes the proof of (21), and hence the proof of (18) and (16) as well, concluding therefore the lemma.

Appendix B. Boundary conditions and monotonicity

We will be interested in vertices on the boundary of a partial configuration, and in the structure of connected components in the graph induced by in-edges. See Definition39.

Observation 51. Let ${\mathcal{F}}_1$ and ${\mathcal{F}}_2$ be two partial configurations such that $R({\mathcal{F}}_1) = R({\mathcal{F}}_2)$ and $\xi ({\mathcal{F}}_1)=\xi ({\mathcal{F}}_2)$ . Let $F = E\backslash R({\mathcal{F}}_1)$ . Then $\pi _{{\mathcal{F}}_1}$ and $\pi _{{\mathcal{F}}_2}$ have the same projection on $F$ , meaning that for any partial configuration ${\mathcal{A}}\in \Omega _F$ , $\pi _{{\mathcal{F}}_1}({\mathcal{F}}_1 \cup{\mathcal{A}}) = \pi _{{\mathcal{F}}_2}({\mathcal{F}}_2 \cup{\mathcal{A}})$ .

Proof. Let $\mathcal{A}$ be a partial configuration in $\Omega _F$ . Consider the components of $G[{\mathrm{In}}({\mathcal{F}}_1\cup{\mathcal{A}})]$ and $G[{\mathrm{In}}({\mathcal{F}}_2\cup{\mathcal{A}})]$ .

  • The components of $G[{\mathrm{In}}({\mathcal{F}}_1\cup{\mathcal{A}})]$ containing only vertices of $V_R({\mathcal{F}}_1)\backslash \partial{\mathcal{F}}_1$ do not depend on the partial configuration $\mathcal{A}$ . Let $C_1$ be the number of these components. Similarly, the components of $G[{\mathrm{In}}({\mathcal{F}}_2 \cup{\mathcal{A}})]$ containing only vertices of $V_R({\mathcal{F}}_2) \backslash \partial{\mathcal{F}}_2$ do not depend on $\mathcal{A}$ . Let $C_2$ be the number of these components.

  • The components of $G[{\mathrm{In}}({\mathcal{F}}_1\cup{\mathcal{A}})]$ containing only vertices of $V\backslash V_R({\mathcal{F}}_1)$ depend on $\mathcal{A}$ but not on $F$ , so the same components are in $G[{\mathrm{In}}({\mathcal{F}}_2\cup{\mathcal{A}})]$ .

  • Since $\xi ({\mathcal{F}}_1) = \xi ({\mathcal{F}}_2)$ , there is one-to-one correspondence between components of $G[{\mathrm{In}}({\mathcal{F}}_1 \cup{\mathcal{A}})]$ that contain vertices in $\partial{\mathcal{F}}_1$ and components of $G[{\mathrm{In}}({\mathcal{F}}_2 \cup{\mathcal{A}})]$ containing vertices in $\partial{\mathcal{F}}_2$ . Any two components that correspond to each other induce the same component on $F$ .

Thus, for any ${\mathcal{A}} \colon F \rightarrow \{0,1,*\}$ , the number of components in $G[{\mathrm{In}}({\mathcal{A}} \cup{\mathcal{F}}_1)]$ minus $C_1$ is equal to the number of components in $G[{\mathrm{In}}({\mathcal{A}} \cup{\mathcal{F}}_2)]$ minus $C_2$ . Hence $w_G({\mathcal{A}}\cup{\mathcal{F}}_1)$ and $w_G({\mathcal{A}}\cup{\mathcal{F}}_2)$ differ by a constant multiplicative factor, from which the result follows.

Cor 53 gives information about marginals in the case where $R({\mathcal{F}}_1) = R({\mathcal{F}}_2)$ and ${\mathrm{In}}({\mathcal{F}}_1) \subseteq{\mathrm{In}}({\mathcal{F}}_2)$ . The proof follows from Observation52.

Observation 52. Let $f$ be an edge in $E$ . Consider partial configurations ${\mathcal{F}}_1$ and ${\mathcal{F}}_2$ with $R({\mathcal{F}}_1) = R({\mathcal{F}}_2)= E\backslash \{f\}$ and ${\mathrm{In}}({\mathcal{F}}_1) \subseteq{\mathrm{In}}({\mathcal{F}}_2)$ . Let ${\mathcal{A}}_1 \sim \pi _{{\mathcal{F}}_1}$ and ${\mathcal{A}}_2 \sim \pi _{{\mathcal{F}}_2}$ . Then ${\mathbb{P}}( f\in{\mathrm{In}}({\mathcal{A}}_1)) \leq{\mathbb{P}}(f\in{\mathrm{In}}({\mathcal{A}}_2))$ .

Proof. First, observe that if $G[{\mathrm{In}}({\mathcal{F}}_1)]$ has two components connected by $f$ then the vertices in these components are connected in $G[{\mathrm{In}}({\mathcal{F}}_2)\cup \{f\}]$ .

For $j\in \{0,1\}$ , let ${\mathcal{F}}_{i,j}$ be ${\mathcal{F}}_i\cup \{f\mapsto j\}$ . Then, by the observation, $c({\mathcal{F}}_{2,0})-c({\mathcal{F}}_{2,1})\leq c({\mathcal{F}}_{1,0})-c({\mathcal{F}}_{1,1})$ . It follows that

\begin{align*}{\mathbb{P}}(f\in{\mathrm{In}}({\mathcal{A}}_1)) &= \frac{{\mathbb{P}}({\mathcal{F}}_{1,1})}{{\mathbb{P}}({\mathcal{F}}_{1,1})+{\mathbb{P}}({\mathcal{F}}_{1,0})} = \frac{(e^\beta -1)^{|{\mathcal{F}}_{1}^{-1}(1)|+1} q^{c({\mathcal{F}}_{1,1})}}{ (e^\beta -1)^{|{\mathcal{F}}_{1}^{-1}(1)|+1} q^{c({\mathcal{F}}_{1,1})} + (e^\beta -1)^{|{\mathcal{F}}_{1}^{-1}(1)|} q^{c({\mathcal{F}}_{1,0})}} \\[5pt] &= \frac{(e^\beta -1)q^{c({\mathcal{F}}_{1,1})-c({\mathcal{F}}_{1,0})}}{(e^\beta -1)q^{c({\mathcal{F}}_{1,1})-c({\mathcal{F}}_{1,0})} + 1} \leq \frac{(e^\beta -1)q^{c({\mathcal{F}}_{2,1})-c({\mathcal{F}}_{2,0})}}{(e^\beta -1)q^{c({\mathcal{F}}_{2,1})-c({\mathcal{F}}_{2,0})} + 1} \\[5pt] &={\mathbb{P}}(f\in{\mathrm{In}}({\mathcal{A}}_2)) \end{align*}

where the inequality follows from $c({\mathcal{F}}_{1,1})-c({\mathcal{F}}_{1,0})\leq c({\mathcal{F}}_{2,1})-c({\mathcal{F}}_{2,0})$ .

In the language of [Reference Grimmett18], Observation52 says that $\pi$ is $1$ -monotonic. It implies the following monotonicity result.

Corollary 53. Let $f$ be an edge of $E$ . Consider partial configurations ${\mathcal{F}}_1$ and ${\mathcal{F}}_2$ with $R({\mathcal{F}}_1) = R({\mathcal{F}}_2)$ and ${\mathrm{In}}({\mathcal{F}}_1) \subseteq{\mathrm{In}}({\mathcal{F}}_2)$ . Let ${\mathcal{A}}_1 \sim \pi _{{\mathcal{F}}_1}$ and ${\mathcal{A}}_2 \sim \pi _{{\mathcal{F}}_2}$ . Then ${\mathbb{P}}( f\in{\mathrm{In}}({\mathcal{A}}_1)) \leq{\mathbb{P}}(f\in{\mathrm{In}}({\mathcal{A}}_2))$ .

Proof. Use [Reference Grimmett18, Theorem 2.27]. It says that if a RCM measure is $1$ -monotonic, then the measure is monotonic, that is, $\pi _{{\mathcal{F}}_1}(A) \leq \pi _{{\mathcal{F}}_2}(A)$ for any increasing event – that is any $A\subseteq \Omega$ such that whenever ${\mathcal{F}}'_{\!\!1},{\mathcal{F}}'_{\!\!2}$ are configurations with ${\mathrm{In}}({\mathcal{F}}'_{\!\!1})\subseteq{\mathrm{In}}({\mathcal{F}}'_{\!\!2})$ and ${\mathcal{F}}'_{\!\!1}\in A$ , then also ${\mathcal{F}}'_{\!\!2}\in A$ .

Observation52 says that $\pi$ is $1$ -monotonic, thus monotonicity follows by the theorem. To conclude, take $A$ to be the increasing event ‘ $f$ is occupied’ – we can see that this is increasing as whenever, for configurations ${\mathcal{F}}'_{\!\!1},{\mathcal{F}}'_{\!\!2}$ , $f\in{\mathrm{In}}({\mathcal{F}}'_{\!\!1})$ , and ${\mathrm{In}}({\mathcal{F}}'_{\!\!1})\subseteq{\mathrm{In}}({\mathcal{F}}'_{\!\!2})$ , then also $f\in{\mathrm{In}}({\mathcal{F}}'_{\!\!2})$ .

Footnotes

1 Recent results of Bencs, Borbényi, and Csikvári [Reference Bencs, Borbényi and Csikvári1] yield the exact formula $\beta _c=\log \tfrac{q-2}{(q-1)^{1-2/\Delta }-1}$ for all $q\gt 2$ and $\Delta \geq 3$ , which was previously only known for integer $q$ [Reference Galanis, Štefankovic, Vigoda and Yang14].

2 We write $G\sim{\mathcal{G}}_{n,\Delta }$ to denote a graph in ${\mathcal{G}}_{n,\Delta }$ chosen uniformly at random, and we say that a property holds w.h.p. for $G\sim{\mathcal{G}}_{n,\Delta }$ as a shorthand for with probability $1-o_n(1)$ over a graph $G\in{\mathcal{G}}_{n,\Delta }$ chosen uniformly at random.

3 The standard submultiplicative argument to bootstrap the total variation distance goes through using the monotonicity of the RC model (to account for the constraint on the initial configuration), see also [Reference Peres and Winkler26].

4 Note that, for $q=2$ , an $O(n^{10})$ upper bound for the RC dynamics on any graph $G$ was previously known at all temperatures $\beta$ by Guo and Jerrum [Reference Guo and Jerrum19] (see also [Reference Feng, Guo and Wang13]).

5 More precisely, the weight of a configuration ${\mathcal{F}}\;:\;E(B_r(\rho ))\rightarrow \{0,1\}$ in $\hat{\pi }_{B_r(\rho )}$ is proportional to $q^{\hat{c}({\mathcal{F}})}(\mathrm{e}^{\beta }-1)^{|{\mathcal{F}}|}$ where $\hat{c}({\mathcal{F}})$ denotes the number of components in the graph $(B_r(\rho ),{\mathrm{In}}(F))$ that do not include any of the vertices in $S_r(\rho )$ (since all of these belong to the same component in the wired dynamics and hence contribute just a single extra factor of $q$ ).

References

Bencs, F., Borbényi, M. and Csikvári, P. (2022) Random cluster model on regular graphs. Commun. Math. Phys. 146.Google Scholar
Blanca, A., Chen, Z., Štefankovič, D. and Vigoda, E. (2021) The Swendsen-Wang dynamics on trees, In Approximation, Randomization, and Combinatorial Optimization. Algorithms and Techniques (APPROX/RANDOM 2021), vol. 207, pp. 115.Google Scholar
Blanca, A., Galanis, A., Goldberg, L. A., Štefankovič, D., Vigoda, E. and Yang, K. (2020) Sampling in uniqueness from the Potts and random-cluster models on random regular graphs. Siam. J. Discrete. Math. 34(1) 742793.CrossRefGoogle Scholar
Blanca, A. and Gheissari, R. (2021) Random-cluster dynamics on random regular graphs in tree uniqueness. Commun. Math. Phys. 386(2) 12431287.CrossRefGoogle Scholar
Blanca, A. and Gheissari, R. (2024) On the tractability of sampling from the Potts model at low temperatures via random-cluster dynamics. Probab. Theory Relat. Fields. https://doi.org/10.1007/s00440-024-01289-xCrossRefGoogle Scholar
Borgs, C., Chayes, J., Helmuth, T., Perkins, W. and Tetali, P. (2020) Efficient sampling and counting algorithms for the Potts model on ℤ $\mathbb{Z}$ d at all temperatures. In Proceedings of the 52nd Annual ACM SIGACT Symposium on Theory of Computing, (STOC ’20), pp. 738751.CrossRefGoogle Scholar
Carlson, C., Davies, E., Fraiman, N., Kolla, A., Potukuchi, A. and Yap, C. (2022) Algorithms for the ferromagnetic Potts model on expanders. In 63rd Annual Symposium on Foundations of Computer Science (FOCS 2022), pp. 344355.CrossRefGoogle Scholar
Chen, Z., Galanis, A., Goldberg, L. A., Perkins, W., Stewart, J. and Vigoda, E. (2021) Fast algorithms at low temperatures via markov chains†. Random Struct. Algor. 58(2) 294321.CrossRefGoogle Scholar
Chen, Zongchen, Galanis, Andreas, Štefankovič, Daniel and Vigoda, Eric (2022) Sampling colorings and independent sets of random regular bipartite graphs in the non-uniqueness region. In Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms, (SODA ’22), pp. 21982207.CrossRefGoogle Scholar
Coja-Oghlan, A., Galanis, A., Goldberg, L. A., Ravelomanana, J. B., Štefankovič, D. and Vigoda, E. (2023) Metastability of the Potts ferromagnet on random regular graphs. Commun. Math. Phys. 401 185225.CrossRefGoogle Scholar
Coulson, M., Davies, E., Kolla, A., Patel, V. and Regts, G. (2020) Statistical physics approaches to unique games. In 35th Computational Complexity Conference, (CCC 2020), volume 169 of LIPIcs, 127.Google Scholar
Efthymiou, C. (2022) On sampling symmetric Gibbs distributions on sparse random graphs and hypergraphs, 49th International Colloquium on Automata, Languages, and Programming (ICALP 2022), pp. 116.Google Scholar
Feng, W., Guo, H. and Wang, J. (2022) Sampling from the ferromagnetic Ising model with external fields. CoRR, abs/2205.01985, 2022. arXiv:2205.01985 Google Scholar
Galanis, A., Štefankovic, D., Vigoda, E. and Yang, L. (2016) Ferromagnetic Potts model: Refined #bis-hardness and related results. Siam. J. Comput. 45(6) 20042065.CrossRefGoogle Scholar
Gheissari, R. and Sinclair, A. (2022) Low-temperature Ising dynamics with random initializations. In 54th Annual ACM SIGACT Symposium on Theory of Computing, (STOC ’22), 14451458.CrossRefGoogle Scholar
Gheissari, R. and Sinclair, A. (2023) Spatial mixing and the random-cluster dynamics on lattices. In Proceedings of the 2023 Annual ACM-SIAM Symposium on Discrete Algorithms, (SODA ’23), 46064621.CrossRefGoogle Scholar
Goldberg, L. A. and Jerrum, M. (2012) Approximating the partition function of the ferromagnetic Potts model. J. ACM 59(5) 131.CrossRefGoogle Scholar
Grimmett, G. (2006) The random-cluster model. Springer, 333 CrossRefGoogle Scholar
Guo, H. and Jerrum, M. (2017) Random cluster dynamics for the Ising model is rapidly mixing. In Proceedings of the Twenty-Eighth Annual ACM-SIAM Symposium on Discrete Algorithms. (SODA 2017), pp. 18181827.CrossRefGoogle Scholar
Häggström, O. (1996) The random-cluster model on a homogeneous tree. Probab. Theory Rel. 104 231253.CrossRefGoogle Scholar
Helmuth, T., Jenssen, M. and Perkins, W. (2023) Finite-size scaling, phase coexistence, and algorithms for the random cluster model on random graphs. Annales de l’Institut Henri Poincare (B) Probabilites et statistiques 59(2) 817848.Google Scholar
Helmuth, T., Perkins, W. and Regts, G. (2019) Algorithmic Pirogov-Sinai theory. In Proceedings of the 51st Annual ACM SIGACT Symposium on Theory of Computing, (STOC ’19), pp. 10091020.CrossRefGoogle Scholar
Jerrum, M. (2003) Counting, sampling and integrating: algorithms and complexity. Basel: Springer Basel AG.CrossRefGoogle Scholar
Jonasson, J. (1999) The random cluster model on a general graph and a phase transition characterization of nonamenability. Stoch. Proc. Appl. 79(2) 335354.CrossRefGoogle Scholar
Mossel, E. and Sly, A. (2013) Exact thresholds for Ising Gibbs samplers on general graphs. Ann. Probab. 41(1) 294328.CrossRefGoogle Scholar
Peres, Y. and Winkler, P. (2013) Can extra updates delay mixing? Commun. Math. Phys. 323(3) 10071016.CrossRefGoogle Scholar
Trevisan, L. (2016) Lecture notes on graph partitioning, expanders and spectral methods. URL: https://lucatrevisan.github.io/books/expanders-2016.pdf.Google Scholar