Logarithmic Mixing of Random Walks on Dynamical Random Cluster Models


Abstract

We study random walks on dynamically evolving graphs, where the environment is given by a time-dependent subset of the edges of an underlying graph. Concretely, following the recently introduced framework of Lelli and Stauffer, we consider a random walk interacting with a dynamical random-cluster environment, in which edges are updated with rate \(\mu>0\) according to Glauber dynamics with parameters \(p\) and \(q\), and the walker moves at rate 1 but may only traverse edges that are present at the time of the move. This setting introduces strong dependencies between the walk and the environment, as edge-update probabilities depend on the global connectivity structure.

We focus on the case where the underlying graph is a random \(d\)-regular graph and the parameters lie in the subcritical regime \(p < p_\mathrm{u}(q, d)\) where it is known that the Glauber dynamics mixes quickly. Our main result is to show that for any \(\varepsilon >0\) and all \(q \ge 1\), for all \(p\) in the subcritical regime, the mixing time of the joint process is \(\Theta(\log n)\) (in continuous time) whenever \(\mu\geq \varepsilon \log n\). This matches the mixing time of the simple random walk on a static random regular graph, showing that in this regime the evolving environment does not slow down mixing. Our proof is based on a coupling argument that uses path-count techniques to overcome the dependencies in the edge dynamics by controlling the structure of the environment along typical trajectories.

1 Introduction↩︎

Random walks on graphs are a central object of study in computer science and probability theory, providing a fundamental tool for algorithms, sampling procedures, and stochastic processes on discrete structures. Classical theory, however, largely concentrates on static underlying graphs, whereas in many complex networks the underlying structure changes over time. This dynamic evolution introduces new challenges in characterising the mixing properties of random walks. In this paper, we study random walks on a dynamically evolving environment and analyse how changes in the underlying structure affect the walk; the “environment" for the random walk is an evolving edge set \((\eta_t)\) of an underlying graph \(G\) (later, we will focus on random regular graphs).

We consider a well-studied model in which the graph evolves as a random process on a fixed underlying graph \(G=(V,E)\). More precisely, the system is described by a continuous-time Markov process \((\eta_t,X_t)_{t\geq 0}\) where \(\eta_t\in \{0,1\}^E\) encodes the set of present edges at time \(t\) and \(X_t\in V\) denotes the location of the walker at time \(t\). The dynamics are driven by \(|E| + 1\) independent Poisson clocks – one clock \(\mathsf{C}_{\mathrm{w}}\) with rate \(1\) which updates the walker \(X_t\), and a clock \(\mathsf{C}_e\) with rate \(\mu>0\) for each edge \(e\in E\). When an edge clock \(\mathsf{C}_e\) rings, the environment is updated at edge \(e\). Let \((\eta_{t-},X_{t-})\) denote the state just before time \(t\). Whenever the walker clock \(\mathsf{C}_{\mathrm{w}}\) rings, say at time \(t\), the walker picks a neighbour \(u\) of \(X_{t-}\) in \(G\) u.a.r. If the edge \(\{X_{t-},u\}\) is present in \(\eta_{t-}\), then \(X_t=u\). Otherwise, the walker remains at its current location, i.e., \(X_t=X_{t-}\).

Key to the dynamic graph evolution are the edge updates, which are triggered whenever an edge clock \(\mathsf{C}_e\) rings. Previous work has mainly focused on the independent edge-update setting, based on the percolation model, where the edge \(e\) is included with probability \(p\) for some fixed \(p\in (0,1)\), independently from the rest of the configuration \(\eta_t\). Here, we focus on the recently introduced framework of Lelli and Stauffer in which the update of an edge is correlated with the current connected-component structure, according to the random cluster model. For a parameter \(q>1\), let \(p_{\min}\mathrel{\vcenter{:}}= \min\{p, \frac{p}{q(1-p) + p}\}\) and \(p_{\max}\mathrel{\vcenter{:}}= \max\{p, \frac{p}{q(1-p) + p}\}\). If \(e\) is a cut edge in \(\eta_{t-}\), we set \(\eta_t(e)=1\) w.p. \(p_{\min}\) and \(\eta_t(e)=0\) otherwise. If \(e\) is not a cut edge in \(\eta_{t-}\), we set \(\eta_t(e)=1\) w.p. \(p_{\max}\) and \(\eta_t(e)=0\) otherwise.

It can be shown that for \(d\)-regular graphs \(G=(V,E)\) and arbitrary \(q>0\) the process \((\eta_t,X_t)\) converges to the distribution \(\pi_{G,p,q}\times \pi_V\), where \(\pi_V\) is the uniform distribution over \(V\) and \(\pi_{G,p,q}\) is the random cluster measure where \(\kappa(\eta)\) is the number of connected components in \(\eta\) and \[\pi_{G,p,q}(\eta)\propto p^{|\eta|}(1-p)^{|E|-|\eta|}q^{\kappa(\eta)}\; for all \eta\in \{0,1\}^E.\] We are interested in the mixing time \(t_{\mathrm{mix}}^{(\mu, p, q)} (G):= \max_{\eta_0,X_0}\inf \big\{t:\;\|(\eta_t,X_t)-\pi_{G,p,q}\times \pi_V\|_{TV}\leq \tfrac{1}{4}\}\) to get close to stationary from an arbitrary starting state. A standard argument yields that \(t_{\mathrm{mix}}^{(\mu, p, q)} (G)=\Omega(\log n)\) for all graphs \(G\) with diameter \(\Omega(\log n)\), since the walker can move graph distance at most one per walker-clock ring; we will primarily focus on obtaining matching upper bounds in the dynamic setting.

A key feature of the model is the interplay between the walk and the edge dynamics. The process \(\eta_t\) corresponds to Glauber dynamics for the random-cluster model with update rate \(\mu\). When \(q=1\), this reduces to the independent-edge setting, where the environment mixes once each edge has been updated. For \(q>1\), however, the behaviour is substantially more complex: the mixing time of Glauber dynamics is known to be exponential when \(p>p_\mathrm{u}(q, d)\) and it remains open whether rapid mixing holds for all \(p\leq p_\mathrm{u}(q, d)\). For random \(d\)-regular graphs, Blanca and Gheissari [1] showed that the mixing time for \(p<p_\mathrm{u}(q, d)\) is \(O(\tfrac{1}{\mu}\log n)\). This motivates our focus on random regular graphs and the subcritical regime \(p<p_\mathrm{u}(q, d)\).

A key difficulty in analysing the dynamic random walk in the subcritical regime was already highlighted in the work of Lelli and Stauffer on the torus. In this regime, the configuration typically consists of small connected components that evolve over time. Although this prevents the walker from being trapped in moderately-large regions, the iterative reformation of components introduces subtle dependencies between the walk and the environment. In particular, the walker is confined to its current component, while the component structure itself influences the edge dynamics, making it challenging to control the mixing behaviour of the joint process.

Our main result establishes that the full joint process mixes in \(O(\log n)\) time for all \(p\) up to the critical threshold \(p_\mathrm{u}(q, d)\), matching the mixing time of the simple random walk on a static random regular graph. Specifically, for \(q\geq 1\), \(p<p_\mathrm{u}(q, d)\) and sufficiently large \(\mu\), we show that the mixing time of the joint process on a random \(d\)-regular graph is \(O(\log n)\) (with high probability over the choice of the random \(d\)-regular graph).

restatabletheoremmainThm Fix integer \(d\geq 3\) and reals \(q\geq 1\) and \(p,\varepsilon > 0\) such that \(p < p_\mathrm{u}(q, d)\). With probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for every edge update rate \(\mu \geq \varepsilon \log n\), we have \(t_{\mathrm{mix}}^{(\mu, p, q)} (G) = \Theta\left(\log n\right)\).

In particular, our result shows that, throughout the subcritical regime and for sufficiently fast edge updates, the evolving environment does not slow down mixing beyond the static case. Our assumption \(\mu = \Omega(\log n)\) for the edge updates ensures that the environment mixes in \(O(1)\) time, placing us in a regime where the dynamics of the walk, rather than the environment, govern the mixing behaviour. In this regime, we show that the walk mixes optimally, matching the static case. Extending these results to slower environments, in particular \(\mu = O(1)\), is an interesting open problem, since the environment may then become the dominant source of mixing bottlenecks.

Related Work. The independent edge-update setting was first considered by Peres, Stauffer, and Steif in [2], where they established the mixing time on the \(d\)-dimensional torus of side-length \(n\) for the subcritical regime. Subsequent work in this setting obtained sharp results on the mixing time (and other related quantities) for various families of graphs when \(q = 1\), notably the work of Peres, Sousi, and Steif [3] on the torus and Sousi and Thomas [4] on the complete graph in the supercritical regime. Hermon and Sousi [5] devised, for \(q=1\), a general comparison principle between the dynamic and static settings on arbitrary underlying graphs.

The analysis of dynamic random walks on evolving random cluster configurations for arbitrary \(q>0\) was initiated by Lelli and Stauffer [6], where they showed optimal mixing-time results on the torus for sufficiently small \(p>0\) and arbitrary \(\mu>0\) (resembling to a certain extent the subcritical regime we consider here). Their proof develops a non-Markovian coupling guided by a multi-scale space-time analysis of the environment, alternating between random-walk coupling steps in favourable regions and identity coupling near unfavourable regions. Our proof of Theorem [thm:main] is also based on a coupling approach, but adapts the path-count method of Lubetzky and Sly [7] to the dynamic setting; we discuss the main obstacles in the next section.

We finally remark that our notion of “dynamic graph" differs from the notion of dynamic graphs studied in other algorithmic settings where vertex/edge updates can be arbitrary. In the random walk context, such updates are typically not ergodic unless some extra assumptions are made. We refer the reader to the work of Sauerwald and Zanetti [8] for progress on this front (see also [9]) and related references. Various other properties beyond mixing time have been studied in this setting. For other dynamical MCMC variants and applications, see [10][15].

2 Proof Overview↩︎

In the rest of the paper we fix an integer \(d\geq 3\) and reals \(q\geq 1\) and \(p, \varepsilon > 0\) such that \(p < p_\mathrm{u}(q, d)\). Since \(q \geq 1\), \(p_{\min}= \tfrac{p}{q(1-p) + p}\) and \(p_{\max}= p\). We assume that \(G=(V,E)\) is sampled from \(\mathcal{G}(d, n)\), so it is \(d\)-regular and has \(n\) vertices. For vertices \(u, v \in V\), \(d(u, v)\) denotes the distance between \(u\) and \(v\) in \(G\). Given \(v \in V\) and a positive integer \(R\), the ball of radius \(R\) around \(v\), denoted by \(B_R(v)\), is the set of vertices with distance at most \(R\) from \(v\). The boundary \(\partial B_R(v)\) is the set of vertices with distance exactly \(R\) from \(v\). From now on, fix \(R = R(d,n) \mathrel{\vcenter{:}}= \frac{1}{5}\log_{d-1} n\). Where the base is not specified, logarithms are base \(e\).

High-Level Strategy. We prove Theorem [thm:main] via a coupling argument, inspired by the approach of Lubetzky and Sly [7] in the static setting. There, a main component in the analysis relies on the fact that random regular graphs are locally tree-like, allowing the walk to be approximated by a biased random walk on the infinite \(d\)-regular tree, whose behaviour can be tracked precisely.

In our dynamic setting, this approach breaks down. The local neighbourhood of the walker evolves over time and can be highly irregular; more importantly, edge updates are no longer independent and the probability that an edge opens depends on its cut status, which in turn depends more globally on the structure of the current configuration.

Our approach to bypass this difficulty is to establish precise lower bounds on the probability that the walker follows prescribed trajectories (paths). To this end, we decompose such trajectory events into realisations of the walker clock and transition choices, together with compatible environment updates. The main challenge is to control the dependence introduced by the random-cluster dynamics, where edge refresh probabilities differ between cut and non-cut edges. Thus, instead of the local tree analysis of the static setting, we need to devise a path-probability lower bound argument that handles more robustly the dependencies that arise in the dynamic setting.

We show that, with high probability, for a large set of “good" starting vertices, typical trajectories of length \(c\log n\) (for appropriate constant \(c\)) avoid small cycles and explore regions that remain tree-like over the relevant time. Environment-wise, we invoke sparsity results from  [1] to show that along such trajectories the edges encountered by the walker are, with high probability, cut edges at their last refresh time and hence are available to the walker with probability \(p_{\min}\). Consequently, most trajectories see cut edges and behave like a walk with bias parameter \(p_{\min}\), giving the key technical estimate that allows us to control transition probabilities despite the dependencies.

The proof is completed via a three-phase coupling. We couple two copies of the chain, one started from an arbitrary initial state and one from stationarity. In the first phase, the chains proceed independently and we show that with constant probability the walker reaches a “good" vertex where we can control the neighbourhood structure and presence of small cycles. In the second phase, we exploit the trajectory estimates above to ensure that the marginal distributions of the two walkers have substantial overlap. Finally, in a short third phase, we use the rapid mixing of the environment (which occurs in \(O(1)\) time in our regime) to couple the full system. Altogether, this yields a coupling with constant success probability in \(O(\log n)\) time.

2.0.0.1 Lower bounding walker trajectory probabilities.

To carry out the previous strategy, a key ingredient in our analysis is to obtain an accurate lower bound on the probability that the walker follows a given trajectory over a prescribed time interval.

For \(\alpha \in \mathbb{N}\), we let \(\overrightarrow{x_\alpha} \in V^{\alpha + 1}\) be a walk \((x_0, \dots, x_\alpha)\) on the underlying graph \(G\) (with repeated vertices allowed) and use \(\ell(\overrightarrow{x_\alpha})\) to denote the number of stationary transitions in the walk (i.e., the number of indices \(i\) such that \(x_i=x_{i+1}\)). For a time interval \((t, t')\), \(N_\mathrm{w}(t, t')\) denotes the number of rings of the walker clock \(\mathsf{C}_{\mathrm{w}}\).

Definition 1 (The walker trajectory event \(\widehat\Xi \left(\overrightarrow{x_\alpha}, t_0, T\right)\)). Let \(\alpha \in \mathbb{N}\) and consider a walk \(\overrightarrow{x_\alpha} \in V^{\alpha + 1}\). For all \(t_0 \geq 0\) and \(T > 0\), \(\widehat\Xi \left(\overrightarrow{x_\alpha}, t_0, T\right)\) is the event that \(N_\mathrm{w}(t_0, t_0 + T) = \alpha\) and the walker follows the walk \(\overrightarrow{x_\alpha}\) during the time interval \((t_0, t_0 + T)\).

When an edge refreshes in the random cluster configuration corresponding to the environment, it opens with probability \(p_{\min}\) if it is a cut edge and with probability \(p_{\max}\) otherwise. If \(q > 1\) then \(p_{\min}< p_{\max}\) so updates depend on the current configuration. Hence bounding the probability of a prescribed trajectory requires controlling whether or not examined edges were cut edges at their last refresh.

An edge fails to be a cut edge if it lies on a cycle (of the random cluster configuration) whose edges are all open at the time that the edge is refreshed. In the regime when \(p < p_\mathrm{u}(q, d)\), this is highly unlikely for long cycles, after a sufficiently long “burn-in" time of the environment dynamics. We take a different approach for short cycles. By restricting to walks that avoid vertices on small cycles, namely cycles of order \(O(\log \log n)\), we ensure that, with high probability, the examined edges are cut edges when they are refreshed. Consequently, along such trajectories, the walker effectively experiences the minimum opening probability \(p_{\min}\).

Definition 2 (\(r\), \(r\)-acyclic walk). Let \(r = {r(q, d, n)}\mathrel{\vcenter{:}}= \frac{3 \log_{d-1} \log n}{\log_{d-1} (2 / (1 + p_\mathrm{u}(q, d)))}\). A walk on a graph is said to be \(r\)-acyclic* if it never visits a cycle of length less than \(r\).*

In Section 3 we prove Lemma [lem:te95prob95lw95bdd], which gives a lower bound on the probability that the walker follows a prescribed \(r\)-acyclic trajectory over a moderate time interval, after the environment has undergone a sufficiently long burn-in. The proof exercises tight control over the probability that examined edges are cut edges, and this yields a probability bound in which the stationary transitions are weighted by \(1 - p_{\min}\) rather than the naive \(1 - p_{\max}\). This improved bound saves a polynomial factor2 in the coupling time (and hence in the mixing time bound), and is the main point of the cut edge analysis.

restatablelemmateProbLwBdd There is a constant \(C_0 \in (0, 1)\) such that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for every \(\mu \geq \varepsilon \log n\), every \(T > 0\) such that \(({1}/{(50 \;p_{\min})}) \log_{d - 1} n \leq T \leq ({4}/{p_{\min}}) \log_{d - 1} n\), every \(\alpha \in \mathbb{N}\) such that \(\alpha \leq ({4}/{p_{\min}}) \log_{d - 1} n\), every \(r\)-acyclic walk \(\overrightarrow{x_\alpha} \in V(G)^{\alpha + 1}\), and every \(\eta \in \{0, 1\}^{E(G)}\), \[\begin{align} \mathbb{P} [\widehat\Xi(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T) & \mid (\eta_0, X_0) = (\eta, x_0)] \geq C_0 \mathbb{P}[N_\mathrm{w}(C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) = \alpha] \left(\frac{p_{\min}^{\alpha - \ell(\overrightarrow{x_\alpha})} (1 - p_{\min})^{\ell(\overrightarrow{x_\alpha})}}{d^{\alpha - \ell(\overrightarrow{x_\alpha})}}\right). \end{align}\]

2.0.0.2 Counts of \(r\)-acyclic walks between “good" vertices.

Let \(G=(V,E)\). Our coupling hinges on the walker quickly hitting a “good" set of vertices \(V_{\mathrm{good}}\subseteq V\) of size \(n - o(n)\) and then becoming almost uniformly mixed over it in \(O(\log n)\) time. We will introduce two properties, namely \((h, i)\)-sparse and \(k\)-root vertices, which will be central to our definition of \(V_{\mathrm{good}}\). We first define some auxiliary notions needed for these properties.

It will be convenient to view a walk on a regular graph \(G\) as a walk on a non-backtracking walk tree (NBWT) of \(G\) (but allowing transitions from any node to its parent). Given \(G = (V, E)\) and \(u \in V\), the NBWT of \(G\) at \(u\), denoted by \(\mathcal{T}_u\), has nodes corresponding to finite non-backtracking walks starting at \(u\), each labelled by its terminal vertex. The root is labelled \(u\) since it corresponds to the trivial walk of length \(0\), and a node corresponding to \((u = x_0, x_1, \dots, x_k)\) has children given by the extensions \((u = x_0, x_1, \dots, x_k, x_{k + 1})\) such that \(\{x_k, x_{k+1}\} \in E\) and \(x_{k+1} \neq x_{k-1}\).

In the dynamic setting, a walk from \(u\) can then be naturally interpreted as a walk on \(\mathcal{T}_u\). We say that it is \((h,i)\)-constrained for non-negative integers \(h\) and \(i\) if it ends at depth \(h\) in \(\mathcal{T}_u\) and has \(h + 2i\) transitions that change the depth by plus or minus one.

Definition 3 (\((h, i)\)-constrained walk). Let \(h, i \in \mathbb{N}\). An \((h, i)\)-constrained walk* on a rooted tree is a walk of length \(h + 2i\) starting at the root of the tree, with exactly \(h + i\) transitions increasing the depth by \(+1\) and \(i\) transitions decreasing the depth by \(-1\). By \(\omega_{h, i}\) we denote the number of \((h,i)\)-constrained walks on the rooted infinite \(d\)-regular tree, while by \(\widetilde{\omega}_{h, i}\) we denote the number of \((h,i)\)-constrained walks on the rooted infinite \((d - 1)\)-ary tree.*

If for a suitably large range of values for \(h, i \in \mathbb{N}\) we have that almost all \((h, i)\)-constrained walks from a vertex \(u\) are \(r\)-acyclic, then Lemma [lem:te95prob95lw95bdd] can be applied for a large range of \(\alpha \in \mathbb{N}\) over almost all length-\(\alpha\) walks from \(u\). This motivates the following definition.

Definition 4 (\((h, i)\)-sparse vertex). A vertex \(u \in V\) is said to be \((h, i)\)-sparse* if at least \((1 - 1/(\log n)^3)\, \omega_{h,i}\) of the \((h, i)\)-constrained walks from \(u\) are \(r\)-acyclic.*

In Appendix 8.2 we will prove Lemma [lem:rrg95glbl95geom], which establishes that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), almost all vertices are \((h, i)\)-sparse for all \(h, i \leq ({4}/{p_{\min}}) \log_{d - 1} n\). Therefore, for every such vertex \(u\) and every \(\alpha \in \mathbb{N}\) such that \(\alpha \leq ({4}/{p_{\min}}) \log_{d - 1} n\), Lemma [lem:te95prob95lw95bdd] applies over almost all length-\(\alpha\) walks from \(u\).

restatablelemmarrgGlblGeom With probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for all \(h, i \in \mathbb{N}\) such that \(h, i \leq ({4}/{p_{\min}}) \log_{d - 1} n\), there are \(n - O\left((d-1)^r (\log n)^8\right)\) vertices which are \(\left(h, i\right)\)-sparse.

Definition 5 (\(k\), \(k\)-root). Let \(k = k(d, n) \mathrel{\vcenter{:}}= \left\lfloor\log_{d-1}\log n\right\rfloor\). A vertex \(v\) of a \(d\)-regular \(n\)-vertex graph \(G\) is a \(k\)-root* if, and only if, the induced subgraph \(G[B_k(v)]\) is a tree. \(V_{\mathrm{roots}}\) denotes the set of all \(k\)-roots in \(G\).*

From now on fix \(H_{\min}(d, n)\mathrel{\vcenter{:}}= \left\lfloor\log_{d-1} n\right\rfloor + 2 \left\lfloor\log_{d-1}\log n\right\rfloor\) and \(H_{\max}(d, n)\mathrel{\vcenter{:}}= \left\lfloor\log_{d-1} n\right\rfloor + \left\lfloor\frac{1}{10}\log_{d-1} n\right\rfloor - 1\). A result from Lubetzky and Sly ([7], Lemmas 3.2 and 3.5) shows that, with probability \(1 - o_n (1)\) over \(G \sim \mathcal{G}(d, n)\), almost all vertices are \(k\)-roots and for all \(u, v \in V_{\mathrm{roots}}\) and \(h, i \in \mathbb{N}\) satisfying \(d(u, v) > 2 k\) and \(h \in [H_{\min}(d, n), H_{\max}(d, n)]\), at least a \(\left(\frac{1 - o_n(1)}{n}\right)\) fraction of all \((h, i)\)-constrained walks from \(u\) end at \(v\). In other words, for a \(k\)-root \(u\) and \(h \in [H_{\min}(d, n), H_{\max}(d, n)]\), the terminal vertex of a uniformly chosen \((h, i)\)-constrained walk from \(u\) is close to uniform over \(n - o(n)\) vertices which are \(k\)-roots. This underpins our rationale for viewing walks of a given length from a vertex \(u\) as families of \((h, i)\)-constrained walks from \(u\). We are now in a position to define the notion of a “good" vertex, central to our analysis.

Definition 6 (Good vertex). A vertex \(u \in V\) is good* if it is a \(k\)-root and it is \((h, i)\)-sparse for all \(h, i \in \mathbb{N}\) such that \(h, i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\). \(V_{\mathrm{good}}\) is the set of all good vertices in the graph \(G\).*

2.0.0.3 Coupling phases.

We next describe our coupling in more detail. We will fix the times \[\begin{align} T_1 &\mathrel{\vcenter{:}}= 1 + C_{\mathrm{burn}}+ \left({d}/{(d - 2)}\right) \left({1}/{p_{\min}}\right) \left({1}/{40}\right) \log_{d - 1} n, \\ T_2 &\mathrel{\vcenter{:}}= T_1 + C_{\mathrm{burn}}+ \left({d}/{(d - 2)}\right) \left({1}/{p_{\min}}\right) \left(\left({81}/{80}\right) \log_{d- 1} n + \log_{d-1} \log n\right), \end{align}\] and an additional time \(T_3\) which we will define in the proof of Theorem [thm:main] — it will be \(T_2 + O(1)\).

Consider two copies of our full chain, \(M \mathrel{\vcenter{:}}= (\eta_t, X_t)\) started at a given state \((\eta, u)\) and \(M' \mathrel{\vcenter{:}}= (\eta'_t, X'_t)\) started from the stationary distribution \(\pi_{G, p, q} \times \pi_V\). We will consider a coupling of these two chains over three phases, of lengths \(T_1\), \(T_2 - T_1\) and \(T_3 - T_2\). We will first run the chains independently until time \(T_1\).

We will show that at time \(T_1\) the walker \(X\) in chain \(M\) is in \(V_{\mathrm{good}}\) with probability \(\Omega(1)\). Now, given \(v \in V_{\mathrm{good}}\) and \(\hat{\eta} \in \{0, 1\}^E\), we consider a maximal coupling between (1) the chain \(M\) under the conditional distribution where \(\left(\eta_{T_1}, X_{T_1}\right) = (\hat{\eta}, v)\) and (2) the chain \(M'\) under the stationary distribution. We will run this maximal coupling until time \(T_3\). We will show that by this point, the chain \(M\) is close to stationarity.

Lemma [lem:phase195walk95counts], which is proved in Section 8, establishes that from every vertex \(u\), among walks of length at most \({3R}/{5}\) which end at a distance \(\frac{R}{10} \leq h \leq \frac{R}{5}\) from \(u\), a constant fraction of them avoid cycles of length less than \(r = {r(q, d, n)}\) (except possibly at \(u\) if it is contained in such a small cycle) and end at a vertex in \(V_{\mathrm{good}}\).

restatablelemmaphaseOneWalkCounts With probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for all vertices \(u\), all \(h, i \in \mathbb{N}\) such that \({R}/{10} \leq h \leq {R}/{5}\), and all \(i \leq {R}/{5}\), there are at least \(d - 2\) neighbours \(w\) of \(u\) such that from each \(w\) there are \({\widetilde{\omega}_{h, i}} / {2}\) walks which are \((h, i)\)-constrained, \(r\)-acyclic and end at a vertex in \(V_{\mathrm{good}}\).

Consequently, from the trajectory probability lower bound in Lemma [lem:te95prob95lw95bdd], we obtain that, with probability \(\Omega(1)\), the walker \(X\) is at a “good" vertex at time \(T_1\). This is captured in Lemma [lem:phase195bdd], which is proved in Section [sec:phase195bdd].

restatablelemmaphaseOneBound There is a constant \(C_1 \in (0, 1)\) such that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for all \(\mu \geq \varepsilon \log n\) and all \((\eta, u) \in V(G) \times \{0, 1\}^{E(G)}\), \(\mathbb{P}\left[X_{T_1} \in V_{\mathrm{good}}\mid \left(\eta_0, X_0\right) = \left(\eta, u\right)\right] \geq C_1\).

We will show that, once the dynamics has reached a vertex \(v \in V_{\mathrm{good}}\) at time \(T_1\), there is a set \(S_v \subseteq V\) with \(|S_v| = n - o(n)\) such that for all \(x \in S_v\), for suitable values of \(h, i \in \mathbb{N}\), there are \(\Omega\left(\left({1}/{n}\right) \omega_{h, i}\right)\) walks from \(v\) to \(x\) which are \((h, i)\)-constrained and \(r\)-acyclic. This allows us to invoke Lemma [lem:te95prob95lw95bdd] over all admissible values of \(h\) and \(i\), giving us a lower bound on the success probability of phase two by aggregating over all such trajectories. This is Lemma [lem:sparse95good95set95size], which is proved in Section 8.4 using Lemmas [lem:rrg95glbl95geom] and 10.

restatablelemmasparseGoodSetSize With probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for any \(u \in V_{\mathrm{good}}\) there exists \(S_u \subseteq V\) such that \(\left|S_u\right| = n - o(n)\) and \(S_u\) has the following property. For all \(h\in \mathbb{N}\) such that \(H_{\min}(d, n)\leq h \leq H_{\max}(d, n)\), all \(i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\), and all \(v \in S_u\), there are at least \(({1}/{4n}) \widetilde{\omega}_{h, i}\) walks from \(u\) to \(v\) which are \((h, i)\)-constrained and \(r\)-acyclic.

Consequently, we obtain Lemma [lem:phase295bdd], which is proved in Section [sec:phase295bdd], and shows that the walker \(X\) is almost uniformly mixed over the set \(S_v\) by time \(T_2\).

restatablelemmaphaseTwoBound There is a constant \(C_2 \in (0, 1)\) such that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for all \(\mu \geq \varepsilon \log n\) and all \(v \in V_{\mathrm{good}}\), there is a set \(S_v \subseteq V(G)\) such that \(|S_v| = n - o(n)\) and, for all \(\left(\hat{\eta}, x\right) \in \{0, 1\}^{E(G)} \times S_v\), \(\mathbb{P}\left[X_{T_2} = x \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \geq {C_2}/{n}\).

2.1 Proof of Theorem [thm:main]↩︎

Proof. Let \(G=(V,E)\). By Lemma [lem:phase195bdd], there is a constant \(C_1 \in (0,1)\) such that, with probability \(1-o_n(1)\) over the choice of \(G\), for all \(u \in V\) and \(\eta \in \{0, 1\}^E\), \[\label{eq:main1} \mathbb{P}\left[X_{T_1} \in V_{\mathrm{good}}\mid \left(\eta_0, X_0\right) = \left(\eta, u\right)\right] \geq C_1.\tag{1}\] By Lemma [lem:phase295bdd] there is a constant \(C_2\in (0,1)\) such that, with probability \(1-o_n(1)\), for all \(v \in V_{\mathrm{good}}\) there is a set \(S_v \subseteq V\) such that \(|S_v| = n - o(n)\) and for all \(x \in S_v\) and \(\hat{\eta} \in \{0, 1\}^E\), \[\label{eq:main2} \mathbb{P}\left[X_{T_2} = x \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \geq {C_2}/{n}.\tag{2}\]

Blanca and Gheissari established in [1] that, with probability \(1-o_n(1)\) over the choice of \(G\), the mixing time of the discrete-time random cluster Glauber dynamics with parameters \((p, q)\) on \(G\) satisfies \(t_{\mathrm{mix}}^{\mathrm{GD}, (p, q)} (G) = \Theta(n \log n)\). This gives a \(\Theta(\log n)\) mixing time in continuous time, and a \(\Theta(1)\) mixing time in our setting with \(\mu \geq \epsilon \log n\). Consequently, there is a constant \(\tau > 0\) such that, with probability \(1 - o_n(1)\) over the choice of \(G\), for every \(\mu \geq \varepsilon \log n\) and every initial distribution \(\eta_0 \sim \rho\), the law of the environment dynamics \((\eta_t^{(\rho)})_{t \geq 0}\) at time \(\tau\) is within total variation distance \(1/8\) of \(\pi_{G, p, q}\).

Let \(T_3 \mathrel{\vcenter{:}}= T_2 + \tau\). Given \(v, x \in V\) and \(\hat{\eta} \in \{0, 1\}^E\), let \(\rho \mathrel{\vcenter{:}}= \rho\left(\hat{\eta}, v, x\right)\) be the law of \(\eta_{T_2}\) conditioned on \((\eta_{T_1}, X_{T_1}) = (\hat{\eta}, v)\) and \(X_{T_2} = x\). Let \[F_{\hat{\eta}, v, x} \! \mathrel{\vcenter{:}}= \! \left\{\omega \in \{0, 1\}^E \colon \mathbb{P}\left[\eta_{T_3} = \omega \mid \eta_{T_2} \sim \rho\left(\hat{\eta}, v, x\right)\right] \geq (1/2) \pi_{G, p, q}(\omega)\right\}.\] For all \(\omega\in \{0,1\}^E\setminus F_{\hat{\eta},v,x}\), \(\min\{\mathbb{P}\left[\eta_{T_3} = \omega \mid \eta_{T_2} \sim \rho(\hat{\eta}, v, x)\right], \pi_{G, p, q}(\omega)\}\) is at most \((1/2)\pi_{G,p,q}(\omega)\) and for \(\omega\in F_{\hat{\eta},v,x}\) this quantity is at most \(\pi_{G,p,q}(\omega)\). Let \(\mathcal{D}(\eta^{(\rho)}_{\tau})\) be the distribution of the environment dynamics, started from the initial distribution \(\rho\), at time \(\tau\). Since \(||\mathcal{D}(\eta^{(\rho)}_{\tau}) - \pi_{G, p, q}||_{\mathrm{TV}} \leq 1/8\) for our choice of \(\tau\), and since \(||\mathcal{D}(\eta^{(\rho)}_{\tau}) - \pi_{G, p, q}||_{\mathrm{TV}} \! = 1 - \!\!\!\!\!\! \sum\limits_{\omega \in \{0, 1\}^E} \!\!\!\!\!\! \min\left\{\mathbb{P}\left[\eta_{T_3} = \omega \mid \eta_{T_2} \sim \rho\right], \pi_{G, p, q}(\omega)\right\}\), we have that \(\pi_{G, p, q} \left(F_{\hat{\eta}, v, x}\right) \geq 3/4\). Since the environment evolves independently of the walker, if the walker clock does not ring between \(T_2\) and \(T_3\), the walker remains at the same vertex while the environment mixes. Thus, for all \(v, x \in V\), \(\hat{\eta} \in \{0, 1\}^E\) and \(\omega \in F_{\hat{\eta}, v, x}\), \[\begin{align} &\;\mathbb{P}\left[\left(\eta_{T_3}, X_{T_3}\right) = (\omega, x) \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \\ \geq& \;\mathbb{P}\left[\eta_{T_3} = \omega, N_{\mathrm{w}}(T_2, T_3) = 0 \mid X_{T_2} = x, \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \mathbb{P}\left[X_{T_2} = x \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \\ =& \;\mathbb{P}\left[N_{\mathrm{w}}(T_2, T_3) = 0\right]\mathbb{P}\left[\eta_{T_3} = \omega \mid X_{T_2} = x, \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \mathbb{P}\left[X_{T_2} = x \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \\ \geq& \;(1/2) e^{-\tau} \pi_{G, p, q}(\omega) \mathbb{P}\left[X_{T_2} = x \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \end{align}\] where the last inequality follows from the fact that \(\omega \in F_{\hat{\eta}, v,x}\) and the fact that the probability of a clock of rate \(1\) not ringing in \(\tau\) time is \(e^{-\tau}\). Combining this with (2 ), we get that for all \(v \in V_{\mathrm{good}}\) there is a set \(S_v \subseteq V\) satisfying \(|S_v| = n - o(n)\) such that, for all \(\left(\hat{\eta}, x\right) \in \{0, 1\}^E \times S_v\) and \(\omega \in F_{\hat{\eta}, v, x}\), \[\label{eq:main4} \mathbb{P}\left[\left(\eta_{T_3}, X_{T_3}\right) = (\omega, x) \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \geq \left(\frac{C_2 e^{-\tau}}{2 n}\right) \pi_{G, p, q}(\omega).\tag{3}\]

Under a maximal coupling from time \(T_1\) to time \(T_3\), from (3 ) and from the fact that \(\pi_V(x) = 1/n\), we have that \[\begin{align} &\;\mathbb{P}\left[\left(\eta_{T_3}, X_{T_3}\right) = \left(\eta'_{T_3}, X'_{T_3}\right) \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right]\\ =&\;\sum_{(\omega, x) \in \{0, 1\}^E \times V} \min\left\{\mathbb{P}\left[\left(\eta_{T_3}, X_{T_3}\right) = (\omega, x) \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right], \pi_{G, p, q} (\omega) \pi_V (x)\right\} \\ \geq& \sum_{x \in S_v} \sum_{\omega \in F_{\hat{\eta},v,x}} \min \{\frac{C_2 e^{-\tau}}{2 n} \pi_{G,p,q}(\omega), \pi_{G,p,q}(\omega) (1/n)\} \geq \;\frac{C_2 e^{-\tau}}{2 n} \sum_{x \in S_v} \pi_{G, p, q} \left(F_{\hat{\eta}, v, x}\right). \end{align}\] Since \(|S_v| = n - o(n)\) and \(\pi_{G, p, q} (F_{\hat{\eta}, v, x}) \geq \frac{3}{4}\) for all \(x\), this is at least \(3 C_2 e^{-\tau}({n - o(n)})/(8n) \geq {C_2 e^{-\tau}}/{8}\). Therefore, the probability of coupling at time \(T_3\) is \[\begin{align} & \;\mathbb{P}\left[\left(\eta_{T_3}, X_{T_3}\right) = \left(\eta'_{T_3}, X'_{T_3}\right) \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \\ \geq& \;\sum_{v \in V_{\mathrm{good}}} \sum_{\hat{\eta} \in \{0, 1\}^E} \mathbb{P}\left[\left(\eta_{T_3}, X_{T_3}\right) = \left(\eta'_{T_3}, X'_{T_3}\right), \left(\eta_{T_1}, X_{T_1}\right) = (\hat{\eta}, v) \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \\ \geq& \;\sum_{v \in V_{\mathrm{good}}} \sum_{\hat{\eta} \in \{0, 1\}^E} \frac{C_2 e^{-\tau}}{8} \;\mathbb{P} \left[\left(\eta_{T_1}, X_{T_1}\right) = (\hat{\eta}, v) \mid \left(\eta_0, X_0\right) = (\eta, u) \right] \\ =& \;\sum_{v \in V_{\mathrm{good}}} \frac{C_2 e^{-\tau}}{8} \;\mathbb{P} \left[X_{T_1} = v \mid \left(\eta_0, X_0\right) = (\eta, u) \right] \\ =& \;\frac{C_2 e^{-\tau}}{8} \;\mathbb{P} \left[X_{T_1} \in V_{\mathrm{good}}\mid \left(\eta_0, X_0\right) = (\eta, u) \right] \geq \frac{C_1 C_2 e^{-\tau}}{8} \end{align}\] where the last inequality follows from (1 ). Hence, with probability \(\Omega(1)\), the two chains couple at time \(T_3\) and therefore the mixing time is \(O(\log n)\). ◻

3 Lower bounding walker trajectory probabilities: Proof of Lemma [lem:te95prob95lw95bdd]↩︎

Consider a \(d\)-regular graph \(G = (V, E)\). Given an edge \(e \in E\) and \(t' > t \geq 0\), by \(\mathsf{C}_e(t, t')\) we denote the event that \(\mathsf{C}_e\) rings in \((t, t')\). Furthermore, conditioned on the event \(\mathsf{C}_e(t, t')\), let \(\tau_e (t, t')\) denote the last ring time of \(\mathsf{C}_e\) in \((t, t')\). Consider \(\alpha \in \mathbb{N}\) and let \(\overrightarrow{x_\alpha} = (x_0, x_1, \dots, x_{\alpha}) \in V^{\alpha + 1}\) be a walk on \(G\), with loops allowed. Suppose that the first \(\alpha\) walker clock rings occur at times \(t_1, \dots, t_\alpha\). Let \(\Delta_i \mathrel{\vcenter{:}}= (t_{i + 1} - t_i) \left(\frac{\varepsilon \log n}{\mu}\right)\). Note that since \(\mu \geq \varepsilon \log n\), then \((t_{i + 1} - \Delta_i, t_{i + 1}) \subseteq (t_i, t_{i + 1})\).

If the walker is at \(x_i\) at time \(t_i\) and \(x_{i + 1} \neq x_i\), for the walker to transition to \(x_{i + 1}\) at time \(t_{i + 1}\), it suffices that the following events hold: the edge \(\{x_i, x_{i + 1}\}\) refreshes during \((t_{i + 1} - \Delta_i, t_{i + 1})\) and it refreshes open during its last ring, and at time \(t_{i + 1}\) the walker examines the edge \(\{x_i, x_{i + 1}\}\), i.e., \(V_{x_i} (t_{i+1}) = x_{i + 1}\) (and hence the walker traverses the edge since it is open). This motivates the following event definition.

Definition 7 (\(\xi_{u \to u'} (t, t')\)). Let \(G = (V, E)\) be a \(d\)-regular graph. Let \(t' > t \geq 0\) and \(u, u' \in V\) such that \(e \mathrel{\vcenter{:}}= \{u, u'\} \in E\). The event \(\xi_{u \to u'} (t, t')\) is the event that all of the following hold: \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\), and \(e\) is refreshed open at time \(\tau_e(t, t')\).

The following lemma, proved in Section 7, gives a lower bound on the probability of the event \(\xi_{u \to u'} (t, t')\), in terms of the edge refresh rate \(\mu\) and the length of the time interval \((t, t')\).

restatablelemmatransitionLwBdd For every \(d\)-regular graph \(G = (V, E)\) and for all \(\mu > 0\), \(t' > t \geq 0\), \(\{u, u'\} \in E\) and \((\eta, v) \in \{0, 1\}^E \times V\), \(\mathbb{P} \left[\xi_{u \to u'} \left(t,t'\right) \mid (\eta_t, X_t) = (\eta, v)\right] \geq ({p_{\min}}/{d}) (1-e^{-\mu (t' - t)})\).

On the other hand, suppose that the walker is at \(x_i\) at time \(t_i\) and that \(x_{i + 1} = x_i\). For the walker to remain at \(x_i\) at time \(t_{i + 1}\), it suffices that for every neighbour \(u_{i + 1}\) of \(x_i\), if the walker examines the edge \(\{x_i, u_{i + 1}\}\) at time \(t_{i + 1}\) (i.e., \(V_{x_i} (t_{i+1}) = u_{i + 1}\)), then it holds that the edge \(\{x_i, u_{i + 1}\}\) was refreshed during \((t_{i + 1} - \Delta_i, t_{i + 1})\) and it was refreshed closed during its last ring (and hence the walker remains at \(x_i\) since the examined edge is closed). This motivates the following event definition.

Definition 8 (\(\xi_{u \not\to u'} (t, t')\)). Let \(G = (V, E)\) be a \(d\)-regular graph. Let \(t' > t \geq 0\) and \(u, u' \in V\) such that \(e \mathrel{\vcenter{:}}= \{u, u'\} \in E\). The event \(\xi_{u \not\to u'} (t, t')\) is the event that all of the following hold: \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\), and \(e\) is refreshed closed at time \(\tau_e(t, t')\).

To obtain good lower bounds on the probability of \(\xi_{u \not\to u'} (t, t')\), we need good lower bounds on the probability that \(e = \{u, u'\}\) refreshes closed at its last refresh time \(\tau_e(t, t')\) during \((t, t')\). Our analysis handles three separate cases, depending on whether \(e\) lies on a cycle, and on the length of that cycle. In particular, we will require the following event for handling long cycles (see Section 7.2 for further details of all the cases that arise).

Given a configuration \(\eta \in \{0, 1\}^E\) and a set of edges \(H \subseteq E\), let \(\eta^H\) be the configuration such that \(\eta^H (e) = 1\) if \(e \in H\), and \(\eta^H (e) = \eta(e)\) if \(e \in E \setminus H\). Given \(v \in V\), let \(E_v \mathrel{\vcenter{:}}= E\left(B_R (v)\right)\), and consider the continuous-time random cluster Glauber dynamics \(\left(\eta_{t}^{E_v}\right)_{t \geq 0}\) with parameters \((p, q)\), that is, the dynamics where all edges in \(E_v\) are wired open. We define \(S_{t} \left(B_R (v), K\right)\) as the event that, in the graph \(\left(V,\eta_{t}^{E_v} \setminus E_v\right)\), at most \(K\) vertices in \(\partial B_R (v)\) are in non-trivial components, and let \(S_{t} (R, K) \mathrel{\vcenter{:}}= \cap_{v \in V} S_{t} (B_R(v), K)\). By \(S_{[t', t'']} (R, K)\) we denote the event that \(S_t (R, K)\) holds for all \(t \in [t', t'']\). See Definition 9 in Section 6 for further details; we will eventually prove a slight adaptation of a result from Blanca and Gheissari ([1], Theorem 5), which establishes that, with high probability over \(G \sim \mathcal{G}(d, n)\), the event \(S_t (R, K)\) holds for a significantly long time interval after some burn-in time \(C_{\mathrm{burn}}\).

The following lemma, proved in Section 7, gives a lower bound on the probability of \(\xi_{u \not\to u'} (t, t')\), in terms of the edge refresh rate \(\mu\), the length of the time interval \((t, t')\), and the probability of the \(K\)-sparsity event \(S_{[t, t']} \left(R, K\right)\).

restatablelemmastationaryLwBdd Let \(n \in \mathbb{N}\) be sufficiently large and let \(G = (V, E)\) be a \(d\)-regular graph on \(n\) vertices. For all \(\mu > 0\), \(t' > t \geq 0\) such that \(\mu (t' - t) > 2 \log 2\), all \(e = \{u, u'\} \in E\) such that \(u\) does not lie on a cycle of length less than \(r\) and \(B_R(u)\) has at most one cycle, and all \((\eta, v) \in \{0, 1\}^E \times V\), \[\begin{align} &\;\mathbb{P} \left[\xi_{u \not\to u'} (t, t') \mid (\eta_t, X_t) = (\eta, v)\right] \\ \geq&\;\frac{1-p_{\min}}{d} \left(1 - \frac{1}{(\log n)^2} - 2(K + 1) e^{-\mu(t' - t)/2}\right) \mathbb{P}\left[S_{[t, t']} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right]. \end{align}\]

Let \(\mathcal{L}\left(\overrightarrow{x_\alpha}\right)\) be the set of indices such that \(i \in \mathcal{L}\left(\overrightarrow{x_\alpha}\right)\) if, and only if, \(x_i = x_{i+1}\). We call a sequence \(\overrightarrow{u_\alpha} = (u_1, \dots, u_\alpha) \in V^\alpha\) a realisation of the walk \(\overrightarrow{x_\alpha}\) if, and only if, \[u_{i+1} \in \begin{cases} \left\{u \in V\colon \{x_i, u\} \in E\right\}, & \mathrm{if} \;i \in \mathcal{L}\left(\overrightarrow{x_\alpha}\right) \\ \{x_{i+1}\}, & \mathrm{otherwise} \end{cases}.\]

We denote the set of all realisations by \(\mathcal{R}\left(\overrightarrow{x_\alpha}\right)\). Let \(\ell\left(\overrightarrow{x_\alpha}\right) = \left|\mathcal{L}\left(\overrightarrow{x_\alpha}\right)\right|\). When there is no room for ambiguity, we will simply write \(\mathcal{R}\), \(\mathcal{L}\) and \(\ell\). Observe that the number of realisations \(|\mathcal{R}|\) for a walk \(\overrightarrow{x_\alpha}\) is \(d^\ell\).

Thus, given that the walker is at \(x_0\) at time \(t_0\) (where \(t_0 < t_1\)), in order for the walker to follow the prescribed trajectory \(\overrightarrow{x_\alpha} = (x_0, \dots, x_\alpha)\) at the ring times \(t_1, \dots, t_\alpha\), it suffices that for some realisation \(\overrightarrow{u_\alpha} \in \mathcal{R}\), for each \(i \in \{0, \dots, \alpha - 1\}\), the appropriate transition event occurs between times \(t_i\) and \(t_{i + 1}\). Specifically, if \(i \notin \mathcal{L}\), then the event \(\xi_{x_i \to x_{i + 1}}(t_i, t_{i + 1})\) holds. On the other hand, if \(i \in \mathcal{L}\), then the event \(\xi_{x_i \not\to u_{i + 1}}(t_i, t_{i + 1})\) holds, where \(u_{i + 1}\) is the neighbour examined at time \(t_{i + 1}\) as specified by the realisation \(\overrightarrow{u_\alpha}\).

Given a realisation \(\overrightarrow{u_\alpha}\) of a walk \(\overrightarrow{x_\alpha}\), \(t_0 \geq 0\) and \(\alpha\) times \(\overrightarrow{t_\alpha} = (t_1, \dots, t_\alpha)\) such that \(t_0 < t_1 < \dots < t_\alpha\), let \(\Delta_i \mathrel{\vcenter{:}}= (t_{i + 1} - t_i) \left(\frac{\varepsilon \log n}{\mu}\right)\) for all \(0 \leq i < \alpha\) and define the trajectory realisation event \[\Xi \! \left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) \mathrel{\vcenter{:}}= \Biggl(\bigcap_{i \not\in \mathcal{L}} \xi_{x_i \to x_{i+1}} \left(t_{i + 1} - \Delta_i, t_{i+1}\right)\Biggr) \cap \Biggl(\bigcap_{i \in \mathcal{L}} \xi_{x_i \not\to u_{i+1}} \left(t_{i + 1} - \Delta_i, t_{i+1}\right)\Biggr).\]

From Lemmas [lem:transition95lw95bdd] and [lem:stationary95lw95bdd] we can deduce the following lower bound on the probability of \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right)\), which we prove in Section 3.1.

restatablelemmaLemmaFourA For all sufficiently large \(n \in \mathbb{N}\), every \(\mu \geq \varepsilon \log n\), every \(T_{\mathrm{max}}\geq T > 0\) and \(t_0 > 0\), every \(d\)-regular \(n\)-vertex graph \(G\) such that every radius-\(R\) ball has at most one cycle, every \(\alpha \in \mathbb{N}\), every \(r\)-acyclic walk \(\overrightarrow{x_\alpha} \in V(G)^{\alpha + 1}\) and realisation \(\overrightarrow{u_{\alpha}} \in \mathcal{R}\), every \(\overrightarrow{t_\alpha} = (t_1, t_2, \dots, t_\alpha)\) satisfying \(t_0 < t_1 < t_2 < \dots < t_\alpha < t_0 + T\) and \((t_{i + 1} - t_i) (\varepsilon \log n) > 2 \log 2\), and every \(\eta \in \{0, 1\}^{E(G)}\), \[\begin{align} &\;\mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \\ \geq& \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right) \! - \sum_{j \in \mathcal{L}} \mathbb{P}\left[\neg S_{[t_{j + 1} - \Delta_j, t_{j + 1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right] \end{align}\] where, for all \(0 \leq i \leq \alpha - 1\), \(T_i \mathrel{\vcenter{:}}= t_{i + 1} - t_i\) and \(\Delta_i \mathrel{\vcenter{:}}= T_i \left(\frac{\varepsilon \log n}{\mu}\right)\).

We define the trajectory event \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}\right) \mathrel{\vcenter{:}}= \bigcup_{\overrightarrow{u_\alpha} \in \mathcal{R}} \Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right)\). We note the following remark.

Remark 1. For any two distinct realisations \(\overrightarrow{u_\alpha}\) and \(\overrightarrow{w_\alpha}\) of a walk \(\overrightarrow{x_\alpha}\), the trajectory realisation events \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right)\) and \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{w_\alpha}\right)\) are disjoint, since there is at least one position \(i \in \mathcal{L}\) such that \(u_{i+1} \neq w_{i+1}\), and therefore the events \(V_{x_i}(t_{i+1}) = u_{i+1}\) and \(V_{x_i}(t_{i+1}) = w_{i+1}\) are disjoint.

Recall from Definition 1 that \(\widehat\Xi \left(\overrightarrow{x_\alpha}, t_0, T\right)\) is the event that \(\mathsf{C}_{\mathrm{w}}\) rings exactly \(\alpha\) times during \((t_0, t_0 + T)\) and the walker follows the walk \(\overrightarrow{x_\alpha}\). By our previous discussions, if the walker is at the initial vertex \(x_0\) of the walk \(\overrightarrow{x_\alpha}\) at time \(t_0\), then if \(\mathsf{C}_{\mathrm{w}}\) rings exactly \(\alpha\) times during \((t_0, t_0 + T)\) and for the \(\alpha\) ring times \(\overrightarrow{t_\alpha}\) in \((t_0, t_0 + T)\) of \(\mathsf{C}_{\mathrm{w}}\) the event \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}\right)\) holds, it follows that \(\widehat\Xi \left(\overrightarrow{x_\alpha}, t_0, T\right)\) holds and at time \(t_0 + T\) the walker is at the end vertex \(x_\alpha\) of the walk, i.e., \(X_{t_0 + T} = x_\alpha\). Consequently, Lemma [lem:lem4a] will serve as the basis of the proof of Lemma [lem:te95prob95lw95bdd].

To pass from fixed ring times of the walker clock to the event that the clock rings \(\alpha\) times in an interval of length \(T\), we must average the fixed-time trajectory bounds over the simplex of possible walker clock ring times. The following estimate, key to the proof of Lemma [lem:te95prob95lw95bdd], says that after excluding very short gaps between successive walker clock rings, the loss from requiring examined edges to refresh during each gap is only a constant factor. Lemma [lem:simplex95integral] is proved in Appendix 9.

restatablelemmasimplexIntegralLemma Fix \(C > 0\). There exists \(c > 0\) such that the following holds: For all \(c' > c\), all sufficiently large \(n \in \mathbb{N}\), all \(T \geq ({1}/{(50 \;p_{\min})}) \log_{d - 1} n\), and all \(\alpha \in \mathbb{N}\) such that \(\alpha \leq ({4}/{p_{\min}}) \log_{d - 1} n\), with \(\delta \mathrel{\vcenter{:}}= {c'}/{(\log n)}\), \[\int_{\substack{T_0, \dots, T_\alpha \geq\delta \\ \sum_{i = 0}^\alpha T_i = T}} \;\prod_{i=0}^\alpha \left(1 - C n^{-(\varepsilon / 2) T_i}\right) \;\mathrm{d} \overrightarrow{T_\alpha} \geq \frac{1}{4} \frac{T^{\alpha}}{\alpha!} \left(1 - \frac{(\alpha + 1) \delta}{T}\right)^\alpha.\]

Remark 2. Given \(t_0, T > 0\) and \(\alpha \in \mathbb{N}\), the simplex \(\{(t_1, \dots, t_\alpha) \colon t_0 < t_1 < \dots < t_\alpha < t_0 + T\}\) is naturally identified with the simplex \(\{T_0, \dots, T_\alpha > 0 \colon \sum\limits_{i = 0}^\alpha T_i = T\}\) via the change of variables \(T_i \mathrel{\vcenter{:}}= t_{i+1} - t_i\) for all \(i \in \{0, 1, \dots, \alpha - 1\}\) and \(T_\alpha = t_0 + T - t_\alpha\), which is linear with unit Jacobian and hence preserves \(\alpha\)-dimensional volume.

We are now in a position to prove our main result for this section.

Proof. Let \(\alpha_{\max}\mathrel{\vcenter{:}}= \lfloor ({4}/{p_{\min}}) \log_{d - 1} n\rfloor\) and \(T_{\mathrm{max}}\mathrel{\vcenter{:}}= ({4}/{p_{\min}}) \log_{d - 1} n\). By Lemma 2, with probability \(1 - O\left(n^{-1}\right)\) over \(G \sim \mathcal{G}(d, n)\), for all \(\mu \geq \varepsilon \log n\), for all \(t' > t \geq C_{\mathrm{burn}}\) such that \(t' - t \leq T_{\mathrm{max}}\) and for all \((\eta, v) \in \{0, 1\}^E \times V\), \[\label{eq:wtp1} \mathbb{P} \left[\neg S_{[t' - \Delta, t']} (R, K) \mid (\eta_0, X_0) = (\eta, v)\right] = O\left(\tfrac{1}{n} \left(\tfrac{p_{\min}}{4 d}\right)^{\alpha_{\max}}\right)\tag{4}\] where \(\Delta \mathrel{\vcenter{:}}= (t' - t) \left(\frac{\varepsilon \log n}{\mu}\right)\), noting that the environment dynamics evolve independently of the walker history. Furthermore, by Lemma 2.1 in [7], we have that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), every radius-\(R\) ball in \(G\) contains at most one cycle.

Hence, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), both of these properties hold. Fix such a graph \(G\). Fix \(\eta \in \{0, 1\}^E\), \(\mu \geq \varepsilon \log n\), \(T > 0\) such that \(({1}/{(50 \;p_{\min})}) \log_{d - 1} n \leq T \leq T_{\mathrm{max}}\), \(\alpha \in \mathbb{N}\) such that \(\alpha \leq \alpha_{\max}\), and let \(\overrightarrow{x_\alpha} = (x_0, \dots, x_\alpha) \in V^{\alpha + 1}\) be an \(r\)-acyclic walk.

Let \(\mathcal{N}_0\) and \(\mathcal{N}_\alpha\) denote the events \(N_\mathrm{w}[0, C_{\mathrm{burn}}] = 0\) and \(N_\mathrm{w}(C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) = \alpha\), respectively. These events are independent since they concern the ringing of the same Poisson clock over disjoint time intervals. Conditioned on the event \(\mathcal{N}_0\), the walker does not move during \([0, C_{\mathrm{burn}}]\), and hence \(X_{C_{\mathrm{burn}}} = x_0\). If, furthermore, the event \(\mathcal{N}_\alpha\) holds and the \(\alpha\) ordered ring times of \(\mathsf{C}_{\mathrm{w}}\) in \((C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T)\) are given by the vector \(\overrightarrow{\tau_\alpha}\), then the event \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{\tau_\alpha}\right)\) forces the walker to follow the trajectory \(\overrightarrow{x_\alpha}\) during \((C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T)\) such that \(X_{C_{\mathrm{burn}}+ T} = x_\alpha\). Therefore these events imply \(\widehat\Xi \left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right)\).

Conditioned on \(\mathcal{N}_\alpha\), the \(\alpha\) ordered ring times of \(\mathsf{C}_{\mathrm{w}}\) occurring during \((C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T)\) are uniformly distributed on the simplex \(\mathfrak{T}_\alpha \mathrel{\vcenter{:}}= \{(t_1, \dots, t_\alpha) \colon C_{\mathrm{burn}}< t_1 < \dots < t_\alpha < C_{\mathrm{burn}}+ T\}\) which has volume \(\mathrm{Vol}\left(\mathfrak{T}_\alpha\right) = {T^\alpha} / {\alpha!}\). We therefore have that \[\begin{align} &\;\mathbb{P}\left[\widehat\Xi \left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \nonumber \\ \geq&\;\mathbb{P}\left[\mathcal{N}_0, \mathcal{N}_\alpha \mid (\eta_0, X_0) = (\eta, x_0)\right] \frac{\alpha!}{T^{\alpha}} \!\! \int\limits_{\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha} \!\!\! \mathbb{P} \left[\Xi(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}) \mid \overrightarrow{\tau_\alpha} = \overrightarrow{t_\alpha}, \mathcal{N}_0, \mathcal{N}_\alpha, (\eta_0, X_0) = (\eta, x_0)\right] \mathrm{d} \overrightarrow{t_\alpha}. \end{align}\]

For fixed \(\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha\), the event \(\Xi(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha})\) is independent of the walker clock, and therefore we have that \(\mathbb{P} \left[\Xi(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}) \mid \overrightarrow{\tau_\alpha} = \overrightarrow{t_\alpha}, \mathcal{N}_0, \mathcal{N}_\alpha, (\eta_0, X_0) = (\eta, x_0)\right] = \mathbb{P} [\Xi(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}) \mid (\eta_0, X_0) = (\eta, x_0)]\). Furthermore, since \(C_{\mathrm{burn}}\) is constant, there exists some constant \(C_1 \in (0, 1)\) such that \(\mathbb{P}[\mathcal{N}_0] \geq C_1\) and hence, by independence of \(\mathcal{N}_0\) and \(\mathcal{N}_\alpha\), we have that \(\mathbb{P}\left[\mathcal{N}_0, \mathcal{N}_\alpha \mid (\eta_0, X_0) = (\eta, x_0)\right] \geq C_1 \mathbb{P}\left[\mathcal{N}_\alpha\right]\).

Combining everything together, it follows that \[\mathbb{P}\left[\widehat\Xi \left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \geq C_1 \mathbb{P}\left[\mathcal{N}_\alpha\right] \frac{\alpha!}{T^{\alpha}} \!\!\! \int\limits_{\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha} \!\!\! \mathbb{P} [\Xi(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}) \mid (\eta_0, X_0) = (\eta, x_0)] \mathrm{d} \overrightarrow{t_\alpha}. \label{eq:wtp2}\tag{5}\]

It will be convenient to define \(t_{\alpha + 1} \mathrel{\vcenter{:}}= C_{\mathrm{burn}}+ T\), and, given \(\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha\), to define, for all \(0 \leq i \leq \alpha\), \(T_i \mathrel{\vcenter{:}}= t_{i + 1} - t_i\) and \(\Delta_i \mathrel{\vcenter{:}}= T_i \left(\frac{\varepsilon \log n}{\mu}\right)\). By Lemma [lem:simplex95integral], there exists \(c > 0\) such that for all \(c' > c\), all \(n\) sufficiently large, with \(\delta \mathrel{\vcenter{:}}= \frac{c'}{\log n}\), \[\int_{\substack{T_0, \dots, T_\alpha \geq\delta \\ \sum_{i = 0}^\alpha T_i = T}} \;\prod_{i=0}^\alpha \left(1 - 4(K+1) n^{-(\varepsilon / 2) T_i}\right) \;\mathrm{d} \overrightarrow{T_\alpha} \geq \frac{1}{4} \frac{T^{\alpha}}{\alpha!} \left(1 - \frac{(\alpha + 1) \delta}{T}\right)^\alpha. \label{eq:wtp10}\tag{6}\]

In particular, we will choose \(c'\) sufficiently large such that \(\delta (\varepsilon \log n) > 2 \log (8(K + 1))\). Let \(\widetilde{\mathfrak{T}}_\alpha\) be the region of the simplex \(\mathfrak{T}_\alpha\) such that \((t_1, \dots, t_\alpha) \in \mathfrak{T}_\alpha\) is in \(\widetilde{\mathfrak{T}}_\alpha\) if, and only if, \(\min\limits_{0 \leq i \leq \alpha} T_i \geq \delta\).

Next fix a realisation \(\overrightarrow{u_\alpha}\) and \(\overrightarrow{t_\alpha} \in \widetilde{\mathfrak{T}}_\alpha\). For this choice of \(\overrightarrow{t_\alpha}\) it follows that, for all \(0 \leq i < \alpha\), since \(T_i (\varepsilon \log n) > 2 \log (8(K + 1))\), then \(T_i (\varepsilon \log n) > 2 \log 2\). By Lemma [lem:lem4a] we have that, \[\begin{align} &\;\mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \nonumber \\ \geq& \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right) \! - \sum_{j \in \mathcal{L}} \mathbb{P}\left[\neg S_{[t_{j + 1} - \Delta_j, t_{j + 1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right]. \label{eq:wtp6} \end{align}\tag{7}\]

From (4 ) we have, for all \(j \in \mathcal{L}\), \(\mathbb{P}\left[\neg S_{[t_{j + 1} - \Delta_j, t_{j + 1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right] = O\left(\frac{1}{n} \left(\frac{p_{\min}}{4 d}\right)^{\alpha_{\max}}\right)\). Furthermore, \(|\mathcal{L}| \leq \alpha = O(\log n)\). Hence, in conjunction with (7 ), \[\begin{align} & \;\mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \nonumber \\ \geq& \;\frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right) - \left(\frac{p_{\min}}{4d}\right)^{\alpha_{\max}} O\left(\tfrac{\log n}{n}\right). \label{eq:wtp7} \end{align}\tag{8}\]

For all \(0 \leq i < \alpha\), since \(T_i (\varepsilon \log n) > 2 \log (8(K + 1))\) then \[2(K + 1)n^{- (\varepsilon / 2) T_i} = 2 (K + 1) e^{- T_i (\varepsilon / 2) \log n} \leq 1/4\] and hence for all \(n\) sufficiently large we have that \(1 - (\log n)^{-2} - 2(K + 1)e^{-\frac{\mu T_i}{2}} \geq 3/4 - (\log n)^{-2} \geq 1/4\). Furthermore, since \(d \geq 3\) and \(p_{\min}< \frac{1}{d-1}\), it follows that \(1 - p_{\min}> p_{\min}\). Therefore, for all \(n\) sufficiently large, we have that \(p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell\prod_{i = 0}^{\alpha - 1} \left(1 - \frac{1}{(\log n)^2} - 2(K+1)e^{-\frac{\mu T_i}{2}}\right) \geq \left(p_{\min}/ 4\right)^{\alpha_{\max}}\). Hence from (8 ), for all \(n\) sufficiently large, \[\begin{align} \mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \geq& (1 - O(\tfrac{\log n}{n})) \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right) \nonumber \\ >& \;\frac{1}{2} \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right). \label{eq:wtp8} \end{align}\tag{9}\]

By the arbitrariness of the realisation \(\overrightarrow{u_\alpha}\) and the disjointness of trajectory realisation events from Remark 1, from (9 ) we obtain that for all \(n\) sufficiently large \[\begin{align} \mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \nonumber \geq& \! \sum_{\overrightarrow{u_\alpha} \in \mathcal{R}} \frac{1}{2} \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right) \nonumber \\ \geq& \;\frac{1}{4} \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^{\alpha - \ell}} \prod_{i = 0}^{\alpha - 1} \left(1 - \frac{4(K+1)}{n^{(\varepsilon / 2) T_i}}\right)\label{eq:wtp9} \end{align}\tag{10}\] where the last inequality follows from the facts that \(|\mathcal{R}| = d^\ell\), and that for all \(x > 0\) and all \(n\) sufficiently large the inequality \(1 - \frac{1}{(\log n)^2} - 2(K+1)n^{-x} \geq \left(1 - \frac{1}{(\log n)^2}\right) \left(1 - 4(K+1)n^{-x} \right)\) holds and since \(\alpha = O(\log n)\), then for all \(n\) sufficiently large we have that \(\left(1 - \frac{1}{(\log n)^2}\right)^\alpha > \frac{1}{2}\).

In conjunction with (6 ) and Remark 2, from (10 ) we have that \[\begin{align} \int\limits_{\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha} \mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \mathrm{d} \overrightarrow{t_\alpha} \geq&\;\frac{1}{4} \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^{\alpha - \ell}} \!\! \int\limits_{\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha} \!\! \prod_{i = 0}^{\alpha - 1} \left(1 - \frac{4(K+1)}{n^{(\varepsilon / 2) T_i}}\right) \mathrm{d} \overrightarrow{t_\alpha} \nonumber \\ \geq& \;\frac{1}{16} \frac{T^\alpha}{\alpha!} \left(1 - \frac{(\alpha + 1) \delta}{T}\right)^\alpha \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^{\alpha - \ell}}. \label{eq:wtp11} \end{align}\tag{11}\]

Since \(\alpha \leq (4/p_{\min}) \log_{d-1} n, T \geq (1/(50 \;p_{\min}) \log_{d - 1} n)\) and \(\delta = c' / \log n = (c' \log_{d-1} e) / \log_{d-1} n\), \(\left(1 - \frac{(\alpha + 1) \delta}{T}\right)^\alpha \geq \left(1 - \left(\frac{c' \log_{d-1} e}{(1/(50 \;p_{\min}))}\right) \left(\frac{1 + (4/p_{\min}) \log_{d-1} n}{(\log_{d-1} n)^2}\right)\right)^{\alpha_{\max}} \geq \left(1 - \frac{400 c' \log_{d-1} e}{\log_{d-1} n}\right)^{(4/p_{\min}) \log_{d-1} n}\) for \(n\) sufficiently large. Hence \(\left(1 - \frac{(\alpha + 1) \delta}{T}\right)^\alpha\) is lower bounded by some constant \(C_2 \in (0, 1)\) for \(n\) sufficiently large. Therefore from (11 ) it follows that for all \(n\) sufficiently large, \[\begin{align} \frac{\alpha!}{T^{\alpha}} \int\limits_{\overrightarrow{t_\alpha} \in \mathfrak{T}_\alpha} \mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \mathrm{d} \overrightarrow{t_\alpha} \geq (C_2 / 16) \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^{\alpha - \ell}}. \label{eq:wtp12} \end{align}\tag{12}\]

The result follows by substituting (12 ) in (5 ). ◻

3.1 Proof of Lemma [lem:lem4a]↩︎

Proof. We will introduce some compact notation. For all \(i \in \mathcal{L}\), let \(\xi_i \mathrel{\vcenter{:}}= \xi_{x_i \not\to u_{i+1}} \left(t_{i + 1} - \Delta_i, t_{i+1}\right)\) and \(\gamma_i \mathrel{\vcenter{:}}= \frac{1-p_{\min}}{d} \left(1 - (\log n)^{-2} - 2(K + 1)n^{-(\varepsilon / 2) T_i}\right)\). For all \(i \notin \mathcal{L}\), let \(\xi_i \mathrel{\vcenter{:}}= \xi_{x_i \to x_{i+1}} \left(t_{i + 1} - \Delta_i, t_{i+1}\right)\) and \(\gamma_i \mathrel{\vcenter{:}}= \frac{p_{\min}}{d} \left(1 - n^{-\varepsilon T_i}\right)\). For all \(0 \leq i < \alpha\), let \(\xi_{[i]} = \bigcap_{j = 0}^i \xi_j\). Observe that \(\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) = \xi_{[\alpha - 1]}\).

By the properties of the graph \(G\) and walk \(\overrightarrow{x_\alpha}\), the radius-\(R\) ball around every vertex in the walk has at most one cycle and does not lie on a cycle of length less than \(r\). Furthermore, for all \(0 \leq i \leq \alpha - 1\), since \(T_i (\varepsilon \log n) > 2 \log 2\) then \(\mu \Delta_i > 2 \log 2\). Therefore, by Lemma [lem:transition95lw95bdd] we have that if \(i \notin \mathcal{L}\) then \[\begin{align} \mathbb{P} \left[\xi_{[i]} \mid (\eta_0, X_0) = (\eta, x_0)\right] &= \mathbb{P} \left[\xi_i \mid \xi_{[i-1]}, (\eta_0, X_0) = (\eta, x_0)\right] \mathbb{P} \left[\xi_{[i-1]} \mid (\eta_0, X_0) = (\eta, x_0)\right] \\ &\geq\;\gamma_i \;\mathbb{P} \left[\xi_{[i-1]} \mid (\eta_0, X_0) = (\eta, x_0)\right]. \end{align}\]

Otherwise by Lemma [lem:stationary95lw95bdd], if \(i \in \mathcal{L}\) we have that \[\begin{align} &\;\mathbb{P} \left[\xi_{[i]} \mid (\eta_0, X_0) = (\eta, x_0)\right] \\ =&\;\mathbb{P} \left[\xi_i \mid \xi_{[i-1]}, (\eta_0, X_0) = (\eta, x_0)\right] \mathbb{P} \left[\xi_{[i-1]} \mid (\eta_0, X_0) = (\eta, x_0)\right] \\ \geq&\;\gamma_i \;\mathbb{P} \left[S_{[t_{i + 1} - \Delta_i, t_{i+1}]} (R, K), \xi_{[i-1]} \mid (\eta_0, X_0) = (\eta, x_0)\right] \\ \geq&\;\gamma_i \bigl(\mathbb{P} \left[\xi_{[i-1]} \mid (\eta_0, X_0) = (\eta, x_0)\right]\! - \mathbb{P}\left[\neg S_{[t_{i + 1} - \Delta_i, t_{i+1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right]\bigr). \end{align}\]

Therefore we inductively have that \[\mathbb{P} \left[\xi_{[\alpha - 1]} \mid (\eta_0, X_0) = (\eta, x_0)\right] \geq \prod_{i = 0}^{\alpha - 1} \gamma_i - \sum_{j \in \mathcal{L}} \left(\prod_{i = j}^{\alpha - 1} \gamma_i \right) \mathbb{P}\left[\neg S_{[t_{j + 1} - \Delta_j, t_{j + 1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right].\]

Since \(\gamma_i \in [0, 1]\) for all \(0 \leq i \leq \alpha - 1\), we obtain that \[\label{eq:l4a1} \mathbb{P} \left[\xi_{[\alpha - 1]} \mid (\eta_0, X_0) = (\eta, x_0)\right] \geq \prod_{i = 0}^{\alpha - 1} \gamma_i - \sum_{j \in \mathcal{L}} \mathbb{P}\left[\neg S_{[t_{j + 1} - \Delta_j, t_{j + 1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right].\tag{13}\]

Now, for all \(n\) sufficiently large, since for all \(0 \leq i \leq \alpha - 1\) we have that \(T_i \varepsilon \log n > 2 \log 2\), then \(1 - n^{-\varepsilon T_i} \geq 1 - n^{-(\varepsilon / 2) T_i} \geq 1 - (\log n)^{-2} - 2(K+1) n^{-(\varepsilon / 2) T_i}\). Then for all \(i \notin \mathcal{L}\) it follows that \(\gamma_i = \frac{p_{\min}}{d} \left(1 - n^{-\varepsilon T_i}\right) \geq \frac{p_{\min}}{d} \left(1 - (\log n)^{-2} - 2(K+1)n^{-(\varepsilon / 2) T_i}\right)\). By this lower bound on \(\gamma_i\) when \(i \notin \mathcal{L}\) and the definition of \(\gamma_i\) when \(i \in \mathcal{L}\), from (13 ) it follows that \[\begin{align} &\;\mathbb{P} \left[\Xi\left(\overrightarrow{x_\alpha}, \overrightarrow{t_\alpha} ; \overrightarrow{u_\alpha}\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \\ \geq& \frac{p_{\min}^{\alpha - \ell} (1 - p_{\min})^\ell}{d^\alpha} \prod_{i = 0}^{\alpha - 1} \left(1 - \tfrac{1}{(\log n)^2} - \frac{2(K+1)}{n^{(\varepsilon / 2) T_i}}\right) \! - \sum_{j \in \mathcal{L}} \mathbb{P}\left[\neg S_{[t_{j + 1} - \Delta_j, t_{j + 1}]} (R, K) \mid (\eta_0, X_0) = (\eta, x_0)\right] \end{align}\] as required. ◻

4 Coupling success probability lower bounds: Lemmas [lem:phase195bdd] and [lem:phase295bdd]↩︎

In this section we will prove the main technical results for our coupling, namely Lemmas [lem:phase195bdd] and [lem:phase295bdd]. We will first prove Lemma 1 — Lemmas [lem:phase195bdd] and Lemma [lem:phase295bdd] are applications of this lemma. Lemma 1 says that if a set \(S\) of vertices has a positive fraction \(\zeta \in (0, 1)\) of \(r\)-acyclic \((h, i)\)-constrained walks from a vertex \(u\), for suitable values of \(h\) and \(i\), then with probability proportional to \(\zeta\), after an appropriate amount of time, the walker transitions from \(u\) to \(S\).

In order to prove Lemma 1 (and hence Lemmas [lem:phase195bdd] and [lem:phase295bdd]), we first require the following two technical lemmas which are proved in Appendix 10. Lemma [lem:f95alpha95properties] will be used to show that the dominant contribution in the lower bound of Lemma 1 comes from walks of length \(\Theta(\log n)\) which are \((h, i)\)-constrained for “typical" values of \(h\) and \(i\), namely values in some window of order \(\Theta(\sqrt{\log n})\). Lemma [lem:poisson95concentration] ensures that the walker clock has constant probability of ringing the desired number of times in the desired time window.

restatablelemmafAlphaLem Let \(a, b\in (0, 1)\) such that \(a + b < {1}/{2}\) and \(\alpha \in \mathbb{N}\). For all \(x, y \in \mathbb{R}_+\) such that \(x + y < \alpha\), define: \[f_\alpha(x,y) \mathrel{\vcenter{:}}= x \log\left({a}/{x}\right) + y\log \left({b}/{y}\right) + (\alpha - x - y)\log\left(\frac{1-a-b}{\alpha - x - y}\right).\] Then,

(i) \(f_\alpha(x, y)\) attains a maximum at \((x, y) = \left(a \alpha, b \alpha\right)\).

(ii) For \(\alpha\) sufficiently large, for all \(0 \leq x, y \leq \sqrt{\alpha}\), \(f_\alpha\left(a \alpha, b \alpha\right) - f_\alpha\left(a (\alpha + x), b (\alpha + y)\right) \leq \frac{a+b}{2(1-a-b)}\).

restatablelemmapoissonConcentration Let \(C_{\max}> C_{\min}> 0\) and \(C_{\mathrm{mid}}> 0\) such that \(C_{\mathrm{mid}}\leq C_{\max}- C_{\min}\). Let \(n \in \mathbb{N}\) be sufficiently large and let \(k_{\min}, k_{\max}\in \mathbb{N}\) such that \(C_{\min}\log n \leq k_{\min}< k_{\max}\leq C_{\max}\log n\) and \(k_{\max}- k_{\min}\geq C_{\mathrm{mid}}\log n\). Let \(Y \sim \mathsf{Poisson}(1)\) and let \(N_Y (T)\) denote the number of rings of \(Y\) during an interval of length \(T > 0\). Then, \(\mathbb{P}\left[k_{\min}\leq N_Y \left(\frac{k_{\min}+ k_{\max}}{2}\right) \leq k_{\max}\right] \geq \frac{1}{2}\).

We are now in a position to prove Lemma 1.

Lemma 1. There exists a constant \(C \in (0, 1)\) such that for all \(n \in \mathbb{N}\) sufficiently large, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for all \(\mu \geq \varepsilon \log n\), and for every \(S \subseteq V\) and \(u \in V\) such that there exists \(\alpha_{\min}, \alpha_{\max}\in \mathbb{N}\) and \(\zeta \in (0, 1)\) satisfying the following:

(i) \(({1}/{(50 \;p_{\min})}) \log_{d - 1} n \leq \alpha_{\min}< \alpha_{\max}\leq ({4}/{p_{\min}}) \log_{d - 1} n\) and \(\alpha_{\max}- \alpha_{\min}\geq \frac{\log_{d-1} n}{100 \;p_{\min}}\)

(ii) for all \(\alpha, h, i \in \mathbb{N}\) such that \(\alpha \in [\alpha_{\min}, \alpha_{\max}]\), \(h \in \left[\left((d - 2) / d\right) p_{\min}\alpha, \left((d - 2) / d\right) p_{\min}\left(\alpha + \sqrt{\alpha}\right)\right]\) and \(i \in \left[(p_{\min}/ d) \alpha, (p_{\min}/ d) \left(\alpha + \sqrt{\alpha}\right)\right]\), there are at least \(\zeta \;\widetilde{\omega}_{h, i}\) walks from \(u\) to a vertex in \(S\) which are \(r\)-acyclic and are \((h, i)\)-constrained walks

then for all \(\eta \in \{0, 1\}^E\), \(\mathbb{P} \left[X_{C_{\mathrm{burn}}+ \frac{\alpha_{\max}+ \alpha_{\min}}{2}} \in S \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \geq C \zeta\).

Proof. By Lemma [lem:te95prob95lw95bdd], there exists \(C_0 \in (0, 1)\) such that for \(n \in \mathbb{N}\) sufficiently large, with probability \(1 - o_n(1)\) over \(G = (V, E) \sim \mathcal{G}(d, n)\), for every \(\mu \geq \varepsilon \log n\), \(T > 0\) such that \(({1}/{(50 \;p_{\min})}) \log_{d - 1} n \leq T \leq ({4}/{p_{\min}}) \log_{d - 1} n\) and \(\alpha \in \mathbb{N}\) such that \(\alpha \leq ({4}/{p_{\min}}) \log_{d - 1} n\), and for all \(r\)-acyclic walks \(\overrightarrow{x_\alpha} \in V^{\alpha + 1}\) and \(\eta \in \{0, 1\}^E\), \[\label{eq:dcw0} \mathbb{P} \left[\widehat\Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right) \mid (\eta_0, X_0) = (\eta, x_0)\right] \geq C_0 \mathbb{P}\left[N_\mathrm{w}(C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) = \alpha\right] \frac{p_{\min}^{\alpha - \ell}(1 - p_{\min})^{\ell}}{d^{\alpha - \ell}}.\tag{14}\]

Fix such a graph \(G\) and \(\mu \geq\varepsilon \log n\). Let \(S \subseteq V\) and \(u \in V\) such that there exists \(\alpha_{\min}, \alpha_{\max}\in \mathbb{N}\) and \(\zeta \in (0, 1)\) satisfying conditions (i) and (ii) in the theorem. For convenience, for all \(x > 0\) define the functions

i. \(h_{\min}(x) \mathrel{\vcenter{:}}= \left(\frac{d-2}{d}\right) p_{\min}x\) and \(h_{\max}(x) \mathrel{\vcenter{:}}= \left(\frac{d-2}{d}\right) p_{\min}\left(x + \sqrt{x}\right)\);

ii. \(i_{\min}(x) \mathrel{\vcenter{:}}= \frac{p_{\min}}{d} x\) and \(i_{\max}(x) \mathrel{\vcenter{:}}= \frac{p_{\min}}{d} \left(x + \sqrt{x}\right)\).

Fix \(T \mathrel{\vcenter{:}}= \frac{\alpha_{\min}+ \alpha_{\max}}{2}\), and for all \(\alpha \in \mathbb{N}\) such that \(\alpha \in [\alpha_{\min}, \alpha_{\max}]\), let \[\mathcal{I}\left( \alpha \right) \mathrel{\vcenter{:}}= \left\{(h, i, \alpha - h - 2i) \in \mathbb{N}^3 \colon h_{\min}(\alpha) \leq h \leq h_{\max}(\alpha), i_{\min}(\alpha) \leq i \leq i_{\max}(\alpha)\right\}.\]

Furthermore, for all \((h, i, \ell) \in \mathcal{I}\left( \alpha \right)\), let \(\mathcal{W}_{h, i, \ell}\) be the collection of walks of length \(\alpha = h + 2i + \ell\) from \(u\) to \(S\) which are \(r\)-acyclic, have exactly \(h + 2i\) transitions corresponding to an \((h, i)\)-constrained walk, and exactly \(\ell\) stationary transitions. Noting that there are \(\binom{h + 2i + \ell}{h + 2i}\) ways of inserting \(\ell\) stationary transitions along an \((h, i)\)-constrained walk, by (ii) it follows that \(\left|\mathcal{W}_{h, i, \ell}\right | \geq \zeta \binom{h + 2i + \ell}{h+2i} \widetilde{\omega}_{h, i}\).

In particular then, by (14 ), we have that for all \(\alpha \in \mathbb{N}\) such that \(\alpha_{\min}\leq \alpha \leq \alpha_{\max}\), for all \((h, i, \ell) \in \mathcal{I}\left( \alpha \right)\), for all \(\overrightarrow{x_\alpha} \in \mathcal{W}_{h, i, \ell}\) and for all \(\eta \in \{0, 1\}^E\): \[\label{eq:dcw1} \mathbb{P} \left[\widehat\Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right) \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \geq C_0 \mathbb{P} \left[N_\mathrm{w} \left(C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T\right) = \alpha\right] \left(\frac{p_{\min}}{d}\right)^{h + 2i} (1 - p_{\min})^\ell.\tag{15}\]

Next fix \(\eta \in \{0, 1\}^E\), \(\alpha \in \mathbb{N}\) and a walk \(\overrightarrow{x_\alpha}\) starting at \(u\) and ending at some vertex in \(S\). Recall that, by our discussion in Section 5, that if \(\widehat \Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T\right)\) holds then the walker follows the trajectory \(\overrightarrow{x_\alpha}\) during the time interval \((C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T)\) and is at the end-vertex of the walk, i.e. in \(S\), at time \(C_{\mathrm{burn}}+ T\). We therefore have that, \[\mathbb{P} \left[X_{C_{\mathrm{burn}}+ T} \in S \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \geq \mathbb{P} \biggl[ \;\bigcup_{\alpha = \alpha_{\min}}^{\alpha_{\max}} \bigcup_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \bigcup_{\overrightarrow{x_\alpha} \in \mathcal{W}_{h, i, \ell}} \!\!\!\! \widehat\Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T \right)\bigg| \left(\eta_0, X_0\right) = (\eta, u) \biggr]. \label{eq:dcw3}\tag{16}\]

Next note that for any \(\alpha, \alpha' \in \mathbb{N}\) such that \(\alpha \neq \alpha'\), the events \(N_\mathrm{w}(C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) = \alpha\) and \(N_\mathrm{w}(C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) = \alpha'\) are disjoint. Also, for any \((h, i, \ell), (h', i', l') \in \mathcal{I}\left( \alpha \right)\) such that \((h, i, \ell) \neq (h', i', l')\), given any \(\overrightarrow{x_\alpha} \in \mathcal{W}_{h, i, \ell}\) and \(\overrightarrow{y_\alpha} \in \mathcal{W}_{h', i', l'}\), the events \(\widehat\Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right)\) and \(\widehat\Xi\left(\overrightarrow{y_\alpha}, C_{\mathrm{burn}}, T\right)\) are disjoint. This follows from the fact that either \(h \neq h'\) and hence the number of non-stationary transitions is distinct, or else, if \(h' = h\) then \(l \neq l'\) and therefore the number of stationary transitions is distinct.

Lastly, we have that for any \(\overrightarrow{x_\alpha}, \overrightarrow{y_\alpha} \in \mathcal{W}_{h, i, \ell}\) such that \(\overrightarrow{x_\alpha} \neq \overrightarrow{y_\alpha}\), the events \(\widehat\Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right)\) and \(\widehat\Xi\left(\overrightarrow{y_\alpha}, C_{\mathrm{burn}}, T\right)\) are disjoint. This follows from the fact that either the non-stationary transitions of both walks differ in their indices, or else they match and hence the walks must differ in at least one vertex.

Therefore from (16 ) we have that, \[\mathbb{P} \left[X_{C_{\mathrm{burn}}+ T} \in S \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \geq \sum_{\alpha = \alpha_{\min}}^{\alpha_{\max}} \sum_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \sum_{\overrightarrow{x_\alpha} \in \mathcal{W}_{h, i, \ell}} \mathbb{P} \left[ \widehat\Xi\left(\overrightarrow{x_\alpha}, C_{\mathrm{burn}}, T\right) \mid \left(\eta_0, X_0\right) = (\eta, u) \right]. \label{eq:dcw5}\tag{17}\]

Together with (15 ) and the fact that for all \((h, i, \ell) \in \mathcal{I}\left( \alpha \right)\) we have that \(|\mathcal{W}_{h, i, \ell}| > \zeta \binom{h + 2i + \ell}{h + 2i} \widetilde{\omega}_{h, i}\), \[\begin{align} &\;\mathbb{P} \left[X_{C_{\mathrm{burn}}+ T} \in S \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \nonumber \\ \geq& \;C_0 \;\zeta \sum_{\alpha = \alpha_{\min}}^{\alpha_{\max}} \mathbb{P} \left[N_\mathrm{w} (C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) = \alpha\right] \sum_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \binom{h + 2i + \ell}{h+2i}\widetilde{\omega}_{h, i} \left(\frac{p_{\min}}{d}\right)^{h+2i} (1 - p_{\min})^\ell. \label{eq:dcw6} \end{align}\tag{18}\]

Let \(a \mathrel{\vcenter{:}}= \left(1-\frac{1}{d}\right) p_{\min}\) and \(b \mathrel{\vcenter{:}}= \frac{p_{\min}}{d}\). By the lower bound on \(\widetilde{\omega}_{h,i}\) in Lemma 8, our choice \(a\) and \(b\), Stirling’s approximation for the factorial, along with the definition of \(f_\alpha\) in Lemma [lem:f95alpha95properties], we have that for all \(n\) sufficiently large: \[\begin{align} &\sum_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \binom{h + 2i + \ell}{h+2i}\widetilde{\omega}_{h, i} \left(\frac{p_{\min}}{d}\right)^{h+2i} (1 - p_{\min})^\ell \\ \geq& \sum\limits_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \frac{h+1}{h+i+1} \frac{(h + 2i + \ell)!}{(h+i)! \;i! \;\ell!} \left(\frac{(d-1)p_{\min}}{d}\right)^{h+i} \left(\frac{p_{\min}}{d}\right)^{i} (1-p_{\min})^\ell \\ \geq& \;\frac{1}{2 \pi} \;\alpha^{\alpha} \! \sum\limits_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \frac{h+1}{h+i+1} \sqrt{\frac{\alpha}{(h+i) \;i \;\ell}} \left(\frac{a}{h+i}\right)^{h+i} \left(\frac{b}{i}\right)^{i} \left(\frac{1-a-b}{\alpha - h - 2i}\right)^{\alpha - h - 2i} \\ \geq& \;C_1 \;\alpha^{\alpha - 1} \! \sum\limits_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \exp\Big(f_\alpha (h+i, i)\Big) \end{align}\] where \(C_1 \mathrel{\vcenter{:}}= \frac{1}{3200 \pi}\), since \(\sqrt{\frac{\alpha}{(h+i) \;i \;\ell}} \geq \frac{1}{\alpha}\) and \(\frac{h+1}{h+i+1} \geq \frac{h_{\min}\left(({1}/{(50 \;p_{\min})}) \log n\right)+1}{h_{\max}\left(({4}/{p_{\min}}) \log_{d-1} n\right) + i_{\max}\left(({4}/{p_{\min}}) \log_{d-1} n\right) + 1} \geq 2 \pi C_1\).

Letting \(C_2 \mathrel{\vcenter{:}}= \left(\frac{d-2}{2}\right)\left(\frac{p_{\min}}{d}\right)^2\), for all \(n\) sufficiently large we have \(|\mathcal{I}\left( \alpha \right)| \geq C_2 \alpha\). Also, letting \(C_3 \mathrel{\vcenter{:}}= \exp\left(-\frac{a+b}{2(1-a-b)}\right)\), by Lemma [lem:f95alpha95properties] we have that for all \((h, i, \ell) \in \mathcal{I}\left( \alpha \right)\): \(\exp\left(f_\alpha (h + i, i)\right) \geq C_3 \exp\left(f_\alpha (a \alpha, b \alpha)\right) \geq C_3 \alpha^{-\alpha}\). Therefore, \[\sum\limits_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \exp\left(f_\alpha (h + i, i)\right) \geq C_2 C_3 \alpha \exp\left(f_\alpha (a \alpha, b \alpha)\right) \geq \frac{C_2 C_3}{\alpha^{\alpha-1}}.\]

Combining everything together, we obtain that: \[\label{eq:dcw7} \sum_{(h, i, \ell) \in \mathcal{I}\left( \alpha \right)} \binom{h + 2i + \ell}{h+2i}\widetilde{\omega}_{h, i} \left(\frac{p_{\min}}{d}\right)^{h+2i} (1 - p_{\min})^\ell \geq C_1 C_2 C_3.\tag{19}\]

From (18 ) and (19 ), defining \(C \mathrel{\vcenter{:}}= \frac{C_0 C_1 C_2 C_3}{2}\), we have that \[\begin{align} \mathbb{P} \left[X_{C_{\mathrm{burn}}+ T} \in S \mid \left(\eta_0, X_0\right) = (\eta, u)\right] &\geq 2 C \zeta \sum_{\alpha = \alpha_{\min}}^{\alpha_{\max}} \mathbb{P} \left[N_\mathrm{w} (C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) \alpha\right] \\ &= 2 C \zeta \cdot \mathbb{P} \left[N_\mathrm{w} (C_{\mathrm{burn}}, C_{\mathrm{burn}}+ T) \in \left[\alpha_{\min}, \alpha_{\max}\right]\right] \\ &\geq C \zeta \end{align}\] where the last inequality follows from Lemma [lem:poisson95concentration]. The result follows. ◻

We are now in a position to prove Lemmas [lem:phase195bdd] and [lem:phase295bdd]. Recall that we fixed the times \[\begin{align} T_1 &\mathrel{\vcenter{:}}= (1 + C_{\mathrm{burn}}) + \left({d}/{(d - 2)}\right) \left({1}/{p_{\min}}\right) \left({1}/{40}\right) \log_{d - 1} n \\ \mathrm{and} \quad T_2 &\mathrel{\vcenter{:}}= T_1 + C_{\mathrm{burn}}+ \left({d}/{(d - 2)}\right) \left({1}/{p_{\min}}\right) \left(\left({81}/{80}\right) \log_{d- 1} n + \log_{d-1} \log n\right). \end{align}\]

Proof. By Lemma [lem:phase195walk95counts], with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for any vertex \(u\) and \(h, i \in \mathbb{N}\) such that \({R}/{10} \leq h \leq {R}/{5}\) and \(i \leq {R}/{5}\), there are at least \(d - 2\) neighbours of \(u\) from which there are \(\frac{\widetilde{\omega}_{h, i}}{2}\) walks which are \((h, i)\)-constrained, \(r\)-acyclic and end at a vertex in \(V_{\mathrm{good}}\).

Fix such a graph \(G\), vertex \(u \in V\), and define \(\alpha_{\min}\mathrel{\vcenter{:}}= \left\lceil\frac{1}{50} \left({d}/{d - 2}\right) \left({1}/{p_{\min}}\right) \log_{d - 1} n\right\rceil\) and \(\alpha_{\max}\mathrel{\vcenter{:}}= \left\lfloor\frac{3}{100} \left({d}/{d - 2}\right) \left({1}/{p_{\min}}\right) \log_{d - 1} n \right\rfloor\). Observe that this choice of \(\alpha_{\min}\) and \(\alpha_{\max}\) satisfies condition (i) of Lemma 1. Furthermore, for all \(\alpha, h, i \in \mathbb{N}\) such that \(\alpha \in [\alpha_{\min}, \alpha_{\max}]\), \(i \in \left[\frac{p_{\min}}{d} \alpha, \frac{p_{\min}}{d} \left(\alpha + \sqrt{\alpha}\right)\right]\) and \(h \in \left[\left(\frac{d-2}{d}\right) p_{\min}\alpha, \left(\frac{d-2}{d}\right) p_{\min}\left(\alpha + \sqrt{\alpha}\right)\right]\), it follows that \(h\) and \(i\) satisfy \({R}/{10} \leq h \leq {R}/{5}\) and \(i \leq {R}/{5}\).

Therefore there are \(d - 2\) neighbours of \(u\), say \(u_1, \dots, u_{d-2}\), such that for all \(1 \leq j \leq d - 2\), there are \(\frac{\widetilde{\omega}_{h, i}}{2}\) walks from \(u_j\) to a vertex in \(V_{\mathrm{good}}\) which are \(r\)-acyclic and \((h, i)\)-constrained, and hence condition (ii) of Lemma 1 is satisfied for every \(u_j\) with \(S = V_{\mathrm{good}}\) and \(\zeta = 1/2\).

Hence by Lemma 1, with probability \(1 - o_n(1)\) over the choice of \(G\), noting that we invoke it at time \(t = 1\) and that \(T_1 = 1 + C_{\mathrm{burn}}+ \frac{\alpha_{\max}+ \alpha_{\min}}{2}\), there exists a constant \(C \in (0, 1)\) such that for all \(1 \leq j \leq d - 2\) we have that \[\label{eq:pob2} \mathbb{P} \left[X_{T_1} \in V_{\mathrm{good}}\mid X_1 = u_j \right] \geq {C}/{2}.\tag{20}\]

It remains to show that, with probability \(\Omega(1)\), the walker transitions during \((0, 1)\) from \(u\) to one of the neighbours \(u_1, \dots, u_{d - 2}\) and then stays there until time \(t = 1\). Recall the definition of the transition event \(\xi_{u \to u'} (t, t')\) as in Definition 7 and let \(\sigma_1\) denote the first ring time of the walker clock. If \(N_{\mathrm{w}}(0, 1) = 1\) and \(\xi_{u \to u_j} (0, \sigma_1)\) both hold for some \(1 \leq j \leq d - 2\), the walker clock rings exactly once in \((0, 1)\), the edge \(\{u, u_j\}\) refreshes open during \((0, \sigma_1)\), and at time \(\sigma_1\) the walker examines that edge. Hence the walker transitions from \(u\) to \(u_j\); since there are no further rings in \((\sigma_1, 1)\), it remains there until time \(t = 1\).

Conditioned on \(N_{\mathrm{w}}(0, 1) = 1\), the time \(\sigma_1\) is uniformly distributed on \((0, 1)\). Therefore, for all \(1 \leq j \leq d - 2\) and for all \(\eta \in \{0, 1\}^E\), \[\label{eq:pob3} \mathbb{P} \left[X_1 = u_j \mid \left(\eta_0, X_0\right) = (\eta, u) \right] \geq \int_0^1 \mathbb{P}\left[\xi_{u \to u_j} (0, x) \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \mathbb{P}\left[N_{\mathrm{w}}(0, 1) = 1\right] \mathrm{d} x.\tag{21}\]

By Lemma [lem:transition95lw95bdd] and the fact that \(\mathbb{P}\left[N_{\mathrm{w}}(0, 1) = 1\right] = {1}/{e}\), \[\label{eq:pob4} \left(\int_0^1 \mathbb{P}\left[\xi_{u \to u_j} (0, x) \mid \left(\eta_0, X_0\right) = (\eta, u)\right] \mathrm{d} x\right) \mathbb{P}\left[N_{\mathrm{w}}(0, 1) = 1\right] = \frac{p_{\min}}{e d} \left(1 - \frac{1 - e^{-\mu}}{\mu}\right).\tag{22}\]

Since \(\mu = \Omega(\log n)\), \(\left(1 - \frac{1 - e^{-\mu}}{\mu}\right) \geq {1}/{2}\) for \(n\) sufficiently large. Moreover, for each fixed \(x > 0\), the events \(\xi_{u \to u_1}(0, x), \dots, \xi_{u \to u_{d - 2}}(0, x)\) are pairwise disjoint. Therefore for all \(n\) sufficiently large, by (21 ) and (22 ), we have that \(\mathbb{P} \left[X_1 \in \{u_1, \dots, u_{d - 2}\} \mid \left(\eta_0, X_0\right) = (\eta, u) \right] \geq \left(\frac{d - 2}{d}\right) \frac{p_{\min}}{2 e}\). The result follows. ◻

Proof. By Lemma [lem:sparse95good95set95size], with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for any \(v \in V_{\mathrm{good}}\) there exists \(S_v \subseteq V\) such that \(\left|S_v\right| = n - o(n)\) and for all \(h, i \in \mathbb{N}\) such that \(H_{\min}(d, n)\leq h \leq H_{\max}(d, n)\) and \(i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\), the following holds: For every \(x \in S_v\) there are at least \(\frac{1}{4n} \widetilde{\omega}_{h, i}\) walks from \(v\) to \(x\) which are \((h, i)\)-constrained and \(r\)-acyclic.

Fix such a graph \(G\), vertex \(v \in V_{\mathrm{good}}\), and let \(\alpha_{\min}\mathrel{\vcenter{:}}= \left\lceil\frac{d}{d-2} \frac{1}{p_{\min}} \left(\log_{d-1} n + 2\log_{d-1}\log n \right)\right\rceil\) and \(\alpha_{\max}\mathrel{\vcenter{:}}= \left\lfloor\frac{d}{d-2} \frac{1}{p_{\min}}(1+\frac{1}{40})\log_{d-1} n\right\rfloor\). Observe that this choice of \(\alpha_{\min}\) and \(\alpha_{\max}\) satisfies condition (i) of Lemma 1. Furthermore, for all \(\alpha, i, h \in \mathbb{N}\) such that \(\alpha \in [\alpha_{\min}, \alpha_{\max}]\), \(i \in \left[\frac{p_{\min}}{d} \alpha, \frac{p_{\min}}{d} \left(\alpha + \sqrt{\alpha}\right)\right]\) and \(h \in \left[\left(\frac{d-2}{d}\right) p_{\min}\alpha, \left(\frac{d-2}{d}\right) p_{\min}\left(\alpha + \sqrt{\alpha}\right)\right]\), it follows that \(h\) and \(i\) satisfy \(H_{\min}(d, n)\leq h \leq H_{\max}(d, n)\) and \(i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\).

Therefore there exists \(S_v \subseteq V\) such that \(|S_v| = n - o(n)\) and for all \(x \in S_v\), there are \(\frac{1}{4n}\widetilde{\omega}_{h, i}\) walks from \(v\) to \(x\) which are \(r\)-acyclic and \((h, i)\)-constrained, and hence condition (ii) of Lemma 1 is satisfied with \(S = \{x\}\) and \(\zeta = {1}/{4n}\).

Hence by Lemma 1, with probability \(1 - o_n(1)\) over the choice of \(G\), noting that we invoke it at time \(t = T_1\) and that \(T_2 - T_1 = C_{\mathrm{burn}}+ \frac{\alpha_{\max}+ \alpha_{\min}}{2}\), there exists \(C \in (0, 1)\) such that for all \(x \in S_v\) and \(\hat{\eta} \in \{0, 1\}^E\), we have that \(\mathbb{P}\left[X_{T_2} = x \mid \left(\eta_{T_1}, X_{T_1}\right) = \left(\hat{\eta}, v\right)\right] \geq {C}/{4 n}\). The result follows. ◻

5 Detailed Chain Description↩︎

We next formalise the full Markov chain that we study. For all positive integers \(d\) and \(n\), \(\mathcal{G}(d, n)\) denotes the family of (simple) \(d\)-regular graphs with \(n\) vertices. Given a graph \(G = (V, E)\), an indicator vector \(\eta \in \{0, 1\}^E\) and an edge \(e \in E\), \(\eta+e\) denotes the indicator vector obtained from \(\eta\) by changing the entry of \(e\) to \(1\) (or leaving it alone if it is already \(1\)). Similarly, \(\eta-e\) is the configuration obtained from \(\eta\) by changing the entry of \(e\) to \(0\). Fix \(G = (V, E) \sim \mathcal{G}(d, n)\). We consider a continuous-time Markov chain \(\left(\eta_t, X_t\right)\) with state space \(\{0, 1\}^E \times V\). The indicator vector \(\eta_t\in \{0,1\}^E\) is the environment at time \(t\). It represents the configuration of edges in \(E\) which are open at time \(t\). The vertex \(X_t\) is the walker position at time \(t\).

The dynamics are driven by \(|E| + 1\) independent Poisson clocks – one clock \(\mathsf{C}_{\mathrm{w}}\) with rate \(1\) which updates the walker \(X_t\), and for each edge \(e\in E\), a clock \(\mathsf{C}_e\) with rate \(\mu\) which updates the entry of \(\eta_t\) corresponding to \(e\). Let \(p_{\min}\mathrel{\vcenter{:}}= \min\left\{p, \frac{p}{q(1-p) + p}\right\}\) and \(p_{\max}\mathrel{\vcenter{:}}= \max\left\{p, \frac{p}{q(1-p) + p}\right\}\). If the clock \(\mathsf{C}_e\) rings at time \(t\) then \(\eta_{t-}\) denotes the environment just before \(\mathsf{C}_e\) rings. The dynamics samples a uniform random variable \(U_e (t)\) on \((0, 1)\) to determine whether edge \(e\) is opened or closed in \(\eta_t\), depending on whether or not \(e\) is a cut edge in \(\eta_{t-}\), according to the procedure Environment_Update\((\eta_{t-},e)\).

\(U_e (t) \sim \text{Uniform}(0, 1)\) \(\eta_{t-} + e\) \(\eta_{t-} - e\) \(\eta_{t-} + e\) \(\eta_{t-} - e\)

Let \(X_{t-}\) denote the location of the walker just before \(\mathsf{C}_{\mathrm{w}}\) rings. Whenever the walker clock \(\mathsf{C}_{\mathrm{w}}\) rings, say at time \(t\), the chain samples \(V_{X_{t-}} (t)\) uniformly on the \(d\) neighbours of \(X_{t-}\) in \(G\).3 If the edge \(\{X_{t-}, V_{X_{t-}} (t)\}\) is open in \(\eta_{t-}\), then \(X_t = V_{X_{t-}} (t)\). Otherwise, the walker stays at the same location, i.e., \(X_t = X_{t-}\). The full Markov chain is thus described by the procedure full_chain\((\eta_0,X_0)\).

\(\forall u \in V \colon V_u (t) \sim \text{Uniform}\left(v \in V \colon \{u, v\} \in E\right)\) \(X_t \leftarrow V_{X_{t-}} (t)\) \(\eta_t \leftarrow\)

Observe that the environment dynamics \((\eta_t)\) is a Markov chain whose stationary distribution is the random cluster measure \(\pi_{G,p,q}\), but the walker position \((X_t)\) is not Markovian since its transitions depend on the current configuration. Nevertheless, the full system \((\eta_t, X_t)\) is a Markov chain. Let \(\pi^{\textrm{SRW}}_G\) denote the stationary distribution of the simple random walk on \(G\). Then the stationary distribution of the full system is \(\pi_{G, p, q} \times \pi^{\textrm{SRW}}_G\). Given a regular graph \(G\), \(\pi^{\textrm{SRW}}_G\) is simply the uniform distribution on the vertex set \(V\), which we denote by \(\pi_V\). We use \(t_{\mathrm{mix}}^{(\mu, p, q)} (G)\) to denote the mixing time of the full chain.

6 Sparse boundaries in burned-in environments↩︎

Given a graph \(G = (V, E)\), a configuration \(\eta \in \{0, 1\}^E\) and a set of edges \(H \subseteq E\), recall that by \(\eta^H\) we denote the configuration \(\eta^H (e) = 1\) if \(e \in H\), else \(\eta^H (e) = \eta(e)\) if \(e \in E \setminus H\).

Definition 9 (External \(K\)-sparsity events). Given \(v \in V\), let \(E_v \mathrel{\vcenter{:}}= E\left(B_R (v)\right)\) and consider the discrete-time random cluster Glauber dynamics \(\left(\eta_{\langle t \rangle}^{E_v}\right)_{t \in \mathbb{N}}\) with parameters \((p, q)\), that is, the dynamics where all edges in \(E_v\) are wired open.

For all \(F \subseteq \partial B_R (v)\), we define \(S^F_{\langle t \rangle} \left(B_R(v)\right)\) as the event that, in the graph \(\left(V,\eta_{\langle t \rangle}^{E_v} \setminus E_v\right)\), the vertices in \(\partial B_R (v)\) which are in non-trivial components are exactly \(F\).

We define \(S_{\langle t \rangle} \left(B_R(v), K\right)\) as the event that \(S^F_{\langle t \rangle} \left(B_R(v)\right)\) holds for some \(F \subseteq \partial B_R (u)\) such that \(|F| \leq K\). Furthermore, we define \(S_{\langle t \rangle} (R, K) \mathrel{\vcenter{:}}= \cap_{v \in V} S_{\langle t \rangle} (B_R(v), K)\). For \(t'' > t' > 0\), by \(S_{\langle[t', t'']\rangle} (R, K)\) we denote the event that \(S_{\langle t \rangle} (R, K)\) holds for all \(t \in [t', t'']\).

By \(S^F_t \left(B_R (v)\right)\), \(S_t \left(B_R (v), K\right)\), \(S_t (R, K)\) and \(S_{[t', t'']} (R, K)\) we will denote the continuous-time analogues of each event, respectively.

The following result, implicit from the proof of Proposition 2 and Theorem 5 in [1], asserts that after a sufficiently long “burn-in" for the discrete-time random cluster Glauber dynamics, for any fixed \(\delta, c_1, c_2 > 0\) there exists a constant \(K > 0\) such that with probability \(1 - O(n^{-c_1})\) over the choice of graph, for a given ball of radius \(R\) around a vertex, with probability \(1 - O(n^{-c_2})\) the number of vertices in non-trivial boundary components of this ball is at most \(K\).

restatabletheoremsparseBurnIn Let \(c_1, c_2 > 0\). There exist constants \(D_{\mathrm{burn}}, K > 0\) such that the following holds: For all \(n \in \mathbb{N}\) sufficiently large and for all \(t \in \mathbb{N}\) such that \(t \geq D_\mathrm{burn} n \log n\), with probability \(1 - O\left(n^{-c_1}\right)\) over \(G \sim \mathcal{G}(d, n)\), for all \(\eta' \in \{0, 1\}^E\), \[\mathbb{P} \left[\neg S_{\langle t \rangle} (R, K) \mid \eta_{\langle 0 \rangle} = \eta'\right] = O\left(n^{-c_2}\right).\]

Proof. Here we briefly remark why we can state the results of [1] in this slightly different form. For a boundary condition \(\eta\) on \(\partial B_R(v)\), by \(\mathcal{J}_{B_R (v), \eta}\) we will denote the set of vertices in \(B_R (v)\) which belong to non-singleton components of the partition \(\eta\left[\partial B_R (v)\right]\).

First, their proof relies on the discrete-time joint-revealing procedure on the ball \(B_R(v)\) and its boundary. In particular, their procedure does not reveal any information on the random variables associated with the edges inside \(B_R (v)\) up to time \(t\).

Instead, they consider the interior with all edges wired open, hence maximising connectivity through the ball, and only consider the random variables associated with the exterior edges for the dynamics. This gives a configuration \(\tilde{\eta}_{\langle t \rangle}\) which dominates the true configuration \(\eta_{\langle t \rangle}\), and hence sparsity on \(\tilde{\eta}_{\langle t \rangle}\) via the exterior implies sparsity on \(\eta_{\langle t \rangle}\) via the exterior.

Furthermore, in their proof they show that, jointly over the configuration model and \(\tilde{\eta}_{\langle t \rangle}\), for any \(k > 1\), the probability of \(\left\{\big|\mathcal{J}_{B_R (v), \eta_{\langle t \rangle}}\big| < k\right\}\) is dominated by the probability of the event \(S_{\langle t \rangle} \left(B_R (v), k\right)\), which in turn is bounded above by \(n^{-\delta k \wedge 4}\).

The constant \(4\) that appears in the exponent can be carefully controlled. We desire it to be \(c_1 + c_2 + 2\); let \(M \geq \left\lceil\frac{c_1 + c_2 + 2}{\delta}\right\rceil\). We can then set \(\Lambda\) sufficiently large in the proof of their Proposition 2 (dependent on \(M\), among other parameters), such that the contribution from the term \(\mathbb{P}\left[\sum_{k \leq k_\emptyset} |\mathcal{V}_k| \geq \Lambda \left|\mathcal{V}_0\right|\right]\) is bounded above by \(O(n^{-M})\) by their Lemma 6. The proof of Proposition 2 then proceeds as usual, with the final Markov inequality over the configuration model applied with \(n^{-(c_2 + 1)}\) instead of \(n^{-2}\).

This then gives that, with probability \(1 - O\left(n^{-(c_1 + 1)}\right)\) over \(G\) sampled from the configuration model, \(\sup\limits_{v \in V} \mathbb{P} \left[\neg S_{\langle t \rangle}(B_R(v), K) \mid \eta_{\langle 0 \rangle} = \eta'\right] \leq n^{-(c_2 + 1)}\). From this, one can deduce the analogue of their Theorem 5 by union bounding over the vertices and conditioning on simple graphs in the configuration model. ◻

The following result gives a continuous-time analogue of Theorem [thm:k95sparse95burn95in].

Lemma 2. There are constants \(C_{\mathrm{burn}}, K > 0\) such that the following holds: With probability \(1 - O(1/n)\) over \(G \sim \mathcal{G}(d, n)\), for all \(\mu \geq \varepsilon \log n\), for all \(t' > t \geq C_{\mathrm{burn}}\) such that \(t' - t \leq (4/p_{\min}) \log_{d-1} n\) and for all \(\eta \in \{0, 1\}^E\), \[\mathbb{P} [\neg S_{[t' - \Delta, t']} (R, K) \mid \eta_0 = \eta] = O \Big(\frac{1}{n}\Big(\frac{p_{\min}}{4 d}\Big)^{(4/{p_{\min}}) \log_{d - 1} n}\Big)\] where \(\Delta \mathrel{\vcenter{:}}= (t' - t) \left(\frac{\varepsilon \log n}{\mu}\right)\).

Proof. Recall that we will use \(S_t (R, K)\) and \(S_{\langle t \rangle} (R, K)\) to differentiate between the continuous and discrete-time analogues, respectively, of the external sparsity event in Definition 9.

Let \(\rho \mathrel{\vcenter{:}}= \left(\frac{4}{p_{\min}}\right) \frac{\log \left({4d}/{p_{\min}}\right)}{\log (d - 1)}\). By Theorem [thm:k95sparse95burn95in], there exist constants \(D_\mathrm{burn}, K > 0\) such that for all \(T \in \mathbb{N}\) satisfying \(T \geq D_\mathrm{burn} n \log n\): With probability \(1 - O\left(n^{-3}\right)\) over \(G \sim \mathcal{G}(d, n)\), \[\label{eq:ssp1} \mathbb{P} \left[\neg S_{\langle T \rangle} (R, K) \mid \eta_{\langle 0 \rangle} = \eta\right] = O\left(n^{-(\rho + 3)}\right), \qquad \forall \eta \in \{0, 1\}^E.\tag{23}\]

Let \(C > 0\) be a constant. Let the edge clocks refresh at rate \(\mu \geq\varepsilon \log n\). Let \(t' > t \geq C\) such that \(t' - t \leq (4/p_{\min}) \log_{d-1} n\). Define \(\Delta \mathrel{\vcenter{:}}= (t' - t) \left(\frac{\varepsilon \log n}{\mu}\right)\) and \(t_0 \mathrel{\vcenter{:}}= t' - \Delta - C \left(\frac{\varepsilon \log n}{\mu}\right)\); observe that \(t' - \Delta = t_0 + C \left(\frac{\varepsilon \log n}{\mu}\right)\) and that for all \(\mu \geq \varepsilon \log n\) we have that \(t_0 = t' - (t' - t + C) \left(\frac{\varepsilon \log n}{\mu}\right) \geq t' - (t' - t + C) \geq t - C \geq 0\) since \(t \geq C\).

By the Markov property, we have that, for all \(\eta \in \{0, 1\}^E\), \[\mathbb{P} [\neg S_{[t' - \Delta, t']} (R, K) \mid \eta_0 = \eta] = \sum_{\eta' \in \{0, 1\}^E} \mathbb{P} [\neg S_{[t' - \Delta, t']} (R, K) \mid \eta_{t_0} = \eta'] \mathbb{P}\left[\eta_{t_0} = \eta' \mid \eta_0 = \eta\right]\] and therefore it suffices to show that, for all \(\eta' \in \{0, 1\}^E\), \(\mathbb{P} \left[\neg S_{[t' - \Delta, t']} \left(R, K\right) \mid \eta_{t_0} = \eta'\right] = O\left(n^{-(\rho + 1)}\right)\).

During the interval \(\left[t_0, t' - \Delta\right)\), which has length \(C \left(\frac{\varepsilon \log n}{\mu}\right)\), the number of edge refreshes is \(N_1 (C) \sim \mathsf{Poisson}\left(\left({C d \varepsilon}/{2}\right) n \log n\right)\). Let \(T_\mathrm{burn} \mathrel{\vcenter{:}}= \left\lceil D_\mathrm{burn} n \log n\right\rceil\). Since \(\mathbb{E}\left[N_1 (C)\right] = \left({C d \varepsilon}/{2}\right) n \log n\), then by concentration there exists \(C_{\mathrm{burn}}> 0\) such that \(\mathbb{P}\left[N_1 (C_{\mathrm{burn}}) < T_\mathrm{burn} \right] \leq O\left(n^{-(\rho + 1)}\right)\). Also, there exists \(C_1 > 0\) such that \(\mathbb{P}\left[N_1(C_{\mathrm{burn}}) > C_1 n \log n\right] \leq O\left(n^{-(\rho + 1)}\right)\). Define \(T_{1, \mathrm{max}} \mathrel{\vcenter{:}}= \left\lfloor C_1 n \log n \right\rfloor\).

During the interval \(\left[t' - \Delta, t'\right]\), the number of edge refreshes is \(N_2 \sim \mathsf{Poisson}\left((d \varepsilon / 2) (t' - t) n \log n)\right)\). Since \(t' - t \leq (4 / p_{\min}) \log_{d-1} n\), we have that \(\mathbb{E}\left[N_2\right] \leq \left(\frac{2 d \varepsilon}{p_{\min}\log(d - 1)}\right) n (\log n)^2\), and by concentration there exists \(C_2 > 0\) such that \(\mathbb{P}\left[N_2 > C_2 n (\log n)^2 \right] \leq O\left(n^{-(\rho + 1)}\right)\). Let \(T_{2, \mathrm{max}} \mathrel{\vcenter{:}}= \left\lfloor C_2 n (\log n)^2\right\rfloor\).

Note that, since Poisson clocks are independent over disjoint time intervals, \(N_1 (C_{\mathrm{burn}})\) and \(N_2\) are independent. Given \(T_1, T_2 \in \mathbb{N}\) and \(\eta' \in \{0, 1\}^E\), conditioned on \(N_1 (C_{\mathrm{burn}}) = T_1\) and \(N_2 = T_2\), the continuous-time dynamics on \([t_0, t']\), started at \(\eta_{t_0} = \eta'\), is distributed as the discrete-time chain run for \(T_1 + T_2\) steps from the initial state \(\eta_{\langle 0 \rangle} = \eta'\), where at each step an edge is chosen independently and uniformly at random. Define \(T_\mathrm{max} \mathrel{\vcenter{:}}= T_{1, \mathrm{max}} + T_{2, \mathrm{max}}\). By a simple union bound, for all \(\eta' \in \{0, 1\}^E\), \[\begin{align} & \;\mathbb{P} \left[\neg S_{[t' - \Delta, t']} \left(R, K\right) \mid \eta_{t_0} = \eta'\right] \nonumber \\ =& \sum_{T_1 = 0}^\infty \sum_{T_2 = 0}^\infty \mathbb{P} \left[\neg S_{\langle[T_1, T_1 + T_2]\rangle} \left(R, K\right) \mid \eta_{\langle 0 \rangle} = \eta'\right] \mathbb{P}\left[N_1 (C_{\mathrm{burn}}) = T_1\right] \mathbb{P}\left[N_2 = T_2\right] \nonumber \\ \leq& \sum_{T_1 = 0}^\infty \sum_{T_2 = 0}^\infty \sum_{T = T_1}^{T_1 + T_2} \mathbb{P} \left[\neg S_{\langle T \rangle} \left(R, K\right) \mid \eta_{\langle 0 \rangle} = \eta'\right] \mathbb{P}\left[N_1 (C_{\mathrm{burn}}) = T_1\right] \mathbb{P}\left[N_2 = T_2\right] \nonumber \\ \leq& \;\mathbb{P}\left[N_1 (C_{\mathrm{burn}}) < T_\mathrm{burn}\right] + \mathbb{P}\left[N_1 (C_{\mathrm{burn}}) > T_{1, \mathrm{max}}\right] + \mathbb{P}\left[N_2 > T_{2, \mathrm{max}}\right] \nonumber \\ & \;+ \sum_{T_1 = T_\mathrm{burn}}^{T_{1, \mathrm{max}}} \sum_{T_2 = 0}^{T_{2, \mathrm{max}}} \sum_{T = T_1}^{T_1 + T_2} \mathbb{P} \left[\neg S_{\langle T \rangle} \left(R, K\right) \mid \eta_{\langle 0 \rangle} = \eta'\right] \mathbb{P}\left[N_1 (C_{\mathrm{burn}}) = T_1\right] \mathbb{P}\left[N_2 = T_2\right] \nonumber \\ \leq & \;\mathbb{P}\left[N_1 (C_{\mathrm{burn}}) \notin \left[T_\mathrm{burn}, T_{1, \mathrm{max}}\right]\right] + \mathbb{P}\left[N_2 > T_{2, \mathrm{max}}\right] + \sum_{T = T_\mathrm{burn}}^{T_\mathrm{max}} \mathbb{P} \left[\neg S_{\langle T \rangle} \left(R, K\right) \mid \eta_{\langle 0 \rangle} = \eta'\right]. \label{eq:ssp2} \end{align}\tag{24}\]

For all \(T \in \mathbb{N}\) such that \(T_\mathrm{burn} \leq T \leq T_\mathrm{max}\), let \(\mathcal{G}_{\langle T \rangle}\) be the family of graphs satisfying (23 ). Since \(T_\mathrm{max} = \Theta\left(n (\log n)^2\right)\), then by a simple union bound, with probability \(1 - O\left(n^{-1}\right)\) over \(G \sim \mathcal{G}(d, n)\), \(G\) is in \(\cap_{T = T_\mathrm{burn}}^{T_\mathrm{max}} \mathcal{G}_{\langle T \rangle}\), since for our choice of \(K\), the probability that \(G \notin \mathcal{G}_{\langle T \rangle}\) is \(O\left(n^{-3}\right)\). Now for every such graph \(G\) we have that for all \(T \in \mathbb{N}\) such that \(T_\mathrm{burn} \leq T \leq T_\mathrm{max}\) and for all \(\eta' \in \{0, 1\}^E\), \(\mathbb{P} \left[\neg S_{\langle T \rangle} (R, K) \mid \eta_{\langle 0 \rangle} = \eta'\right] = O\left(n^{-(\rho + 3)}\right)\). Furthermore, by our choices of \(C_{\mathrm{burn}}\), \(T_{1, \mathrm{max}}\) and \(T_{2, \mathrm{max}}\), we have that \(\mathbb{P}\left[N_1 (C_{\mathrm{burn}}) \notin \left[T_\mathrm{burn}, T_{1, \mathrm{max}}\right]\right] = O\left(n^{-(\rho + 1)}\right)\) and \(\mathbb{P}\left[N_2 > T_{2, \mathrm{max}}\right] = O\left(n^{-(\rho + 1)}\right)\).

Therefore, from these facts in conjunction with (24 ), \[\mathbb{P} \left[\neg S_{[t' - \Delta, t']} \left(R, K\right) \mid \eta_{t_0} = \eta'\right] = O\left(n^{-(\rho + 1)}\right) + \sum_{T = T_\mathrm{burn}}^{T_\mathrm{max}} O\left(n^{-(\rho + 3)}\right) = O\left(n^{-(\rho + 1)}\right)\] noting that \(T_\mathrm{max} = \Theta\left(n (\log n)^2\right)\). The result follows. ◻

7 Lower bounding walker transition events: Proofs of Lemmas [lem:transition95lw95bdd] and [lem:stationary95lw95bdd]↩︎

We next turn our attention back to our full chain, with both the walker and environment dynamics. Consider a \(d\)-regular graph \(G = (V, E)\) on \(n\) vertices. Recall that for an edge \(e \in E\) and \(t' > t \geq 0\), by \(\mathsf{C}_e(t, t')\) we denote the event that \(\mathsf{C}_e\) rings in \((t, t')\), and, conditioned on the event \(\mathsf{C}_e(t, t')\), let \(\tau_e (t, t')\) denote the last ring time of \(\mathsf{C}_e\) in \((t, t')\). Given \(e = \{u, u'\} \in E\) and \(t' > t \geq 0\), recall from Definition 7 that \(\xi_{u \to u'} (t, t')\) is the event that all of the following hold: \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\) and \(e\) is refreshed open at time \(\tau_e (t, t')\). Also, from Definition 8, recall that \(\xi_{u \not\to u'} (t, t')\) is the event that all of the following hold: \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\) and \(e\) is refreshed closed at time \(\tau_e (t, t')\).

We will prove Lemmas [lem:transition95lw95bdd] and [lem:stationary95lw95bdd], which establish lower bounds on the probability of the events \(\xi_{u \to u'} (t, t')\) and \(\xi_{u \not\to u'} (t, t')\), respectively.

Before proceeding further, we note the following useful fact about Poisson clocks.

Lemma 3. Let \(\mathsf{C}\) be a Poisson clock of rate \(\mu > 0\). Let \(T > 0\) and let \(\tau\) be the last time \(\mathsf{C}\) rings on \((0, T)\), conditioned on the event \(\mathsf{C}(0, T)\). Then \(\tau\) has cumulative distribution function \[\begin{align} \mathcal{F}_{\tau \vert \mathsf{C}(0, T)} (x) = \mathbb{P}\left[\tau \leq x \mid \mathsf{C}(0, T)\right] = \frac{e^{\mu x} - 1}{e^{\mu T} - 1},& &0 < x < T \end{align}\] and probability density function \(f_{\tau \vert \mathsf{C}(0, T)} (x) = \frac{\mu e^{\mu x}}{e^{\mu T} - 1}\), \(0 < x < T\).

Proof. First note that \(\mathcal{F}_{\tau \vert \mathsf{C}(0, T)} (x) = \mathbb{P}\left[\tau \leq x \mid \mathsf{C}(0, T)\right] = \frac{\mathbb{P}\left[\mathsf{C}(0, x]\right]\mathbb{P}\left[\neg\mathsf{C}(x, T)\right]}{\mathbb{P}\left[\mathsf{C}(0, T)\right]} = \frac{e^{\mu x} - 1}{e^{\mu T} - 1}\). Taking the derivative with respect to \(x\), we obtain that \(f_{\tau \vert \mathsf{C}(0, T)} (x) = \dfrac{\mu e^{\mu x}}{e^{\mu T} - 1}\). ◻

7.1 Proof of Lemma [lem:transition95lw95bdd]↩︎

Proof. Let \(e \mathrel{\vcenter{:}}= \{u, u'\}\). If the edge \(e\) refreshes during \((t, t')\), i.e. the event \(\mathsf{C}_e(t, t')\) holds, and during the last ring \(\tau_e (t, t')\) it holds that \(U_e \left(\tau_e (t, t')\right) < p_{\min}\), then \(e\) is refreshed open independent of whether it is a cut edge or not at time \(\tau_e (t, t')\). Therefore, the events \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\) and \(U_e \left(\tau_e (t, t')\right) < p_{\min}\) imply the event \(\xi_{u \to u'} (t, t')\). Hence, \[\mathbb{P} \left[\xi_{u \to u'} \left(t,t'\right) \mid (\eta_t, X_t) = (\eta, v)\right] \geq \mathbb{P} \left[V_{u}\left(t'\right) = u', \mathsf{C}_e\left(t,t'\right), U_e \left(\tau_e(t, t')\right) < p_{\min}\mid (\eta_t, X_t) = (\eta, v)\right].\]

The events \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\) and \(U_e \left(\tau_e (t, t')\right) < p_{\min}\) concern the ringing of Poisson clocks and sampling of uniforms during the time interval \((t, t')\), which are independent of the state at time \(t\). Furthermore, the random variable \(V_u\) is independent of the environment random variables, namely \(U_e\) and \(\mathsf{C}_e\).

We therefore have that, \[\begin{align} &\;\mathbb{P} \left[V_{u}\left(t'\right) = u', \mathsf{C}_e\left(t,t'\right), U_e \left(\tau_e(t, t')\right) < p_{\min}\mid (\eta_t, X_t) = (\eta, v)\right] \\ =&\;\mathbb{P}\left[V_{u}\left(t'\right) = u'\right] \mathbb{P}\left[U_e \left(\tau_e(t, t')\right) < p_{\min}\mid \mathsf{C}_e\left(t,t'\right)\right] \mathbb{P}\left[\mathsf{C}_e\left(t,t'\right)\right]. \end{align}\]

Since \(V_{u}\left(t'\right)\) is uniform on the \(d\) neighbours of \(u\), then \(\mathbb{P}\left[V_{u}\left(t'\right) = u'\right] = 1/d\). By the probability of a Poisson clock ringing in a given time interval, we have \(\mathbb{P}\left[\mathsf{C}_e\left(t,t'\right)\right] = 1 - e^{-\mu \left(t' - t\right)}\). Lastly, since \(U_e\) is uniform on \((0, 1)\), we have \(\mathbb{P}\left[U_e \left(\tau_e(t, t')\right) < p_{\min}\mid \mathsf{C}_e\left(t,t'\right)\right] = p_{\min}\). The result follows. ◻

7.2 Controlling cut edge event probabilities↩︎

In order to obtain a good lower bound on the probability of \(\xi_{u \not\to u'}(t, t')\), we need to carefully control the probability that \(e = \{u, u'\}\) is a cut edge in the configuration at its last refresh time \(\tau_e (t, t')\), conditioned on the event \(\mathsf{C}_e(t, t')\).

An edge \(e = \{u, v\}\) fails to be a cut edge in a configuration if there exists an open path in the configuration between \(u\) and \(v\) that does not use the edge \(e\). If, at the last refresh time of the edge \(e\), before it is examined by the walker, there were no such open paths between \(u\) and \(v\), then \(e\) was a cut edge at its last refresh time. It therefore suffices to control the probability of open paths between endpoints of edges incident to the walker vertex at the last time of refresh for the edge examined. We proceed to define a series of useful events which together guarantee that an edge is a cut edge in a given configuration. Given an edge \(e' \in E\) and \(t'' > t \geq 0\), let \(\mathcal{E}_{e'} (t, t'')\) be the event that \(\mathsf{C}_{e'} (t, t'')\) and \(U_{e'}(\tau_{e'}(t, t'')) \geq p_{\max}= p\) both hold. We therefore have that, \[\label{eq:ste1} \mathbb{P}\left[\neg \mathcal{E}_{e'}(t, t'')\right] = 1 - (1-p)(1-e^{-\mu (t'' - t)}) = p + (1-p)e^{-\mu (t'' - t)}.\tag{25}\]

In particular, if \(\mathcal{E}_{e'}(t, t'')\) holds then the edge \(e'\) is refreshed during \((t, t'')\) and at its last refresh it is closed independently of whether it was a cut edge or not, and remains closed until time \(t''\).

Remark 3. For any two edges \(e', e'' \in E\), the events \(\mathcal{E}_{e'} (t, t'')\) and \(\mathcal{E}_{e''} (t, t'')\) are independent, by the independence of the corresponding clocks and random variables. More generally, \(\mathcal{E}_{e'} (t, t'')\) is independent of any other event concerning the clocks and uniforms of edges distinct from \(e'\) during \((t, t'')\). Moreover, since \(\mathcal{E}_{e'} (t, t'')\) concerns the ringing of Poisson clocks and sampling of uniforms during the time interval \((t, t'')\), which are independent of the state at time \(t\), then \(\mathcal{E}_{e'} (t, t'')\) is independent of the state at time \(t\).

Given a vertex \(u\) and a cycle \(\mathbf{c}\) in the underlying graph \(G\), let \(E_{u, r} (\mathbf{c})\) be the set of edges on \(\mathbf{c}\) which are distance at most \(r\) away from \(u\) (i.e. for a given edge \(e\) on \(\mathbf{c}\), the shortest path from \(u\) to each end-vertex of \(e\) has length at most \(r\)). If a cycle \(\mathbf{c}\) passing through \(u\) has length greater than \(r\), then \(\left|E_{u, r} (\mathbf{c})\right| \geq r\). Also, given an edge \(e \in E\), let \(C_e (R)\) be the set of cycles in the underlying graph \(G\) with length strictly less than \(2R\) that contain \(e\). Then given a vertex \(u\) and edge \(e\) incident to \(u\), for all \(t'' > t \geq 0\) define the event \[\xi_{u, e}^{\mathrm{acyclic}}\left(t, t''\right) \mathrel{\vcenter{:}}= \bigcap_{\mathbf{c} \in C_e (R)} \bigcup_{\substack{e' \in E_{u, r} (\mathbf{c}) \\ e' \neq e}} \!\! \mathcal{E}_{e'} (t, t'').\] which is the event that during the time interval \((t, t'')\), for every cycle in \(C_e (R)\) at least one edge \(e' \neq e\) on the cycle, distance at most \(r\) away from \(u\), was refreshed and the last time this edge was refreshed in \((t, t'')\) it was closed (independently of whether it was a cut edge or not).

Remark 4. In particular, if \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t''\right)\) holds then for any cycle containing \(e\) which is entirely contained in the ball \(B_R(u)\) (i.e. has length less than \(2R\)), at least one edge, distinct from \(e\), on such a cycle must be closed at time \(t''\).

We next prove the following lower bound for the probability of the event \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right)\).

Lemma 4. For every \(d\)-regular graph \(G = (V, E)\) and for all \(\mu > 0\), \(t \geq 0\), \(x > 0\) and \(e = \{u, u'\} \in E\) such that \(u\) does not lie on a cycle of length less than \(r\) and \(B_R(u)\) has at most one cycle, then \(\mathbb{P} \left[\xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right) \right] \geq 1 - \left(p + (1-p)e^{-\mu x}\right)^{r - 1}\).

Proof. Since there is at most one cycle in \(B_R (u)\), then \(\left|C_e (R)\right| \leq 1\). Moreover, if such a cycle containing \(e\) exists in \(B_R (u)\), its length must be at least \(r\) (otherwise \(u\) lies on a cycle of length less than \(r\)) and hence \(\left|E_{u, r}(\mathbf{c})\right| \geq r\). By complementation and a simple union bound (noting the independence of the \(\mathcal{E}\) events for distinct edges in Remark 3), \[\mathbb{P} \left[\neg \xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right) \right] \leq \sum_{\mathbf{c} \in C_e(R)} \prod_{\substack{e' \in E_{u, r}(\mathbf{c}) \\ e' \neq e}} \mathbb{P} \left[\neg \mathcal{E}_{e'} (t, t + x)\right] \leq \left(p + (1-p)e^{-\mu x}\right)^{r - 1}\] where the last inequality follows using \(\left|C_{e} (R)\right| \leq 1\), the fact that any such cycle must have at least \(r - 1\) edges distinct from \(e\) that are a distance at most \(r\) from \(u\), and the upper bound on \(\mathbb{P}\left[\neg \mathcal{E}_{e'}(t, t + x)\right]\) given in (25 ). The result follows. ◻

Given a vertex \(u\) and a path \(\mathbf{p}\), let \(E_{u, (r, R)}(\mathbf{p})\) be the set of edges on \(\mathbf{p}\) which are at distance at least \(r\) but at most \(R\) away from \(u\) (i.e. for a given edge \(e\) on \(\mathbf{p}\), the shortest path from \(u\) to each end-vertex of \(e\) has length at least \(r\) but at most \(R\)). If a path \(\mathbf{p}\) has one end-vertex at \(u\) and the other in \(\partial B_R(u)\), then \(\left|E_{u, (r, R)}(\mathbf{p})\right| \geq R - r\). Furthermore, given \(F \subseteq \partial B_R (u)\), let \(P\left(B_R (u), F\right)\) be the set of paths in the graph from \(u\) to \(F\) which are contained entirely in \(B_R(u)\). Then, for any edge \(e\) incident to \(u\), any \(t'' > t \geq 0\), and any \(F \subseteq \partial B_R (u)\), define the event \[\mathcal{P}^F_{u, e} (t, t'') \mathrel{\vcenter{:}}= \bigcap_{\mathbf{p} \in P\left(B_R (u), F\right)} \;\bigcup_{\substack{e' \in E_{u, (r, R)}(\mathbf{p}) \\ e' \neq e}} \mathcal{E}_{e'} (t, t'')\] which is the event that during the interval \((t, t'')\), for every path in \(P\left(B_R (u), F\right)\), at least one edge \(e' \neq e\) on the path was refreshed and the last time this edge was refreshed in \((t, t'')\) it was closed (independently of whether it was a cut edge or not).

Next, recall the definition of the external sparsity event \(S^F_t \left(B_R (u)\right)\) from Definition 9. Given \(K > 0\) as in Lemma 2, define the event \[\xi_{u, e}^{\mathrm{path}}\left(t, t''\right) \mathrel{\vcenter{:}}= \bigcup_{\substack{F \subseteq \partial B_R (u) \\ |F| \leq K}} \mathcal{P}^{F}_{u, e} (t, t'') \cap S^F_{t''} \left(B_R (u)\right).\]

Remark 5. In particular, if \(\xi_{u, e}^{\mathrm{path}}\left(t, t''\right)\) holds then for any cycle in the underlying graph which contains \(e\) and is not entirely contained in the ball \(B_R(u)\) (i.e. has length greater than \(2R\)), there is at least one edge, distinct from \(e\), which is closed in the configuration at time \(t''\). Suppose that \(F \subseteq \partial B_R (u)\) is the set of non-trivial boundary vertices on \(B_R(u)\) at time \(t''\). Every such cycle must pass through the boundary of the ball \(B_R (u)\) through at least two distinct vertices. If the cycle does not pass through the boundary at a vertex in \(F\), then at least one edge (distinct from \(e\)) outside the ball on the cycle must be closed in the configuration at time \(t''\) (otherwise there is an open path outside of the ball between non-trivial boundary vertices). Otherwise, if the cycle passes through the boundary at a vertex in \(F\), then the segment of the cycle from \(u\) to this vertex which is contained in \(B_R(u)\) must have had an edge (distinct from \(e\)) refreshed closed during \((t, t'')\), since \(\mathcal{P}^{F}_{u, e} (t, t'')\) holds.

We will next prove the following lemma.

Lemma 5. Let \(n \in \mathbb{N}\) be sufficiently large and \(G = (V, E)\) a \(d\)-regular graph on \(n\) vertices. For all \(\mu > 0\), \(t \geq 0\), \(x > 0\), \(e = \{u, u'\} \in E\) such that \(B_R(u)\) has at most one cycle, and \((\eta, v) \in \{0, 1\}^E \times V\), then \[\begin{align} & \mathbb{P}\left[\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \\ \geq&\;\left(1 - 2K \left(p + (1-p)e^{-\mu x}\right)^{{R}/{2}}\right) \mathbb{P}\left[S_{t+x} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right]. \end{align}\]

Proof. First note that, by disjointness of events over \(F \subseteq \partial B_R (u)\), we have that \[\mathbb{P}\left[\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] = \sum_{\substack{F \subseteq \partial B_R (u) \\ |F| \leq K}} \mathbb{P}\left[\mathcal{P}^{F}_{u, e} (t, t + x) \cap S^F_{t + x} \left(B_R (u)\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right].\]

Now, by Remark 3, the event \(\mathcal{P}^{F}_{u, e} (t, t + x)\) is independent of the state at time \(t\). Furthermore, recall that the external sparsity event \(S^F_{t + x} \left(B_R (u)\right)\) is independent of the edge clocks and uniforms associated with edges inside \(B_R(u)\) up to time \(t + x\). On the other hand, \(\mathcal{P}^{F}_{u, e} (t, t + x)\) is an event concerning the clocks and uniforms associated with edges inside \(B_R(u)\) during \((t, t + x)\). Hence \(S^F_{t + x} \left(B_R (u)\right)\) and \(\mathcal{P}^{F}_{u, e} (t, t + x)\) are also independent. Consequently, we have that \[\mathbb{P}\left[\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right) \mid \left(\eta_t, X_t\right) = (\eta', v)\right] = \sum_{\substack{F \subseteq \partial B_R (u) \\ |F| \leq K}} \!\! \mathbb{P}\left[\mathcal{P}^{F}_{u, e} (t, t + x)\right] \mathbb{P}\left[S^F_{t + x} \left(B_R (u)\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right].\]

Since there is at most one cycle in \(B_R(u)\), then given \(F \subseteq \partial B_R (u)\) such that \(|F| \leq K\), the number of paths inside \(B_R (u)\) from \(u\) to a vertex in \(F\) is at most \(2K\) and hence \(\left|P\left(B_R(u), F\right)\right| \leq 2K\). Indeed, suppose that there is some vertex \(v\) in \(F\) such that there are three or more paths in \(B_R (u)\) from \(u\) to \(v\); then there are at least two cycles in \(B_R(u)\), which is a contradiction.

Next, by complementation and a simple union bound (noting the independence of the \(\mathcal{E}\) events for distinct edges in Remark 3), we have that \[\mathbb{P} \left[\neg \mathcal{P}^{F}_{u, e} (t, t + x) \right] \leq \sum_{\mathbf{p} \in P\left(B_R(u), F\right)} \prod_{\substack{e' \in E_{e, (r, R)}(\mathbf{p}) \\ e' \neq e}} \mathbb{P} \left[\neg \mathcal{E}_{e'} (t, t + x)\right] \text{d} x \leq 2K \left(p + (1-p)e^{-\mu x}\right)^{R - r - 1}\] where the last inequality follows from \(\left|P\left(B_R(u), F\right)\right| \leq 2K\), the fact that for any \(\mathbf{p} \in P\left(B_R(u), F\right)\) we have \(\left|E_{e, (r, R)} (\mathbf{p})\right| \geq R - r\), and the upper bound on \(\mathbb{P}\left[\neg \mathcal{E}_{e'}(t, t + x)\right]\) given in (25 ).

For \(n\) sufficiently large we have \(R - r - 1 \geq {R}/{2}\), since \(R = \Theta(\log n)\) and \(r = \Theta(\log \log n)\). Therefore, for \(n\) sufficiently large, \(\mathbb{P} \left[\neg \mathcal{P}^{F}_{u, e} (t, t + x) \right] \leq 2K \left(p + (1-p)e^{-\mu x}\right)^{{R}/{2}}\). From Definition 9, we have that \(S_{t+x} \left(B_R(u), K\right)\) is the disjoint union of the events \(S^F_{t+x} \left(B_R (u)\right)\) for all \(F \subseteq \partial B_R (u)\) such that \(|F| \leq K\). Combining everything together, we have that \[\begin{align} &\;\mathbb{P}\left[\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right) \mid \left(\eta_t, X_t\right) = (\eta', v)\right] \\ \geq& \left(1 - 2K \left(p + (1-p)e^{-\mu x}\right)^{{R}/{2}}\right) \sum_{\substack{F \subseteq \partial B_R (u) \\ |F| \leq K}} \mathbb{P}\left[S^F_{t + x} \left(B_R (u)\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \\ =& \left(1 - 2K \left(p + (1-p)e^{-\mu x}\right)^{{R}/{2}}\right) \mathbb{P}\left[S_{t+x} \left(B_R (u), K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \\ \geq& \left(1 - 2K \left(p + (1-p)e^{-\mu x}\right)^{{R}/{2}}\right) \mathbb{P}\left[S_{t+x} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \end{align}\] where the last inequality follows from the fact that \(S_{t+x} \left(R, K\right) \subseteq S_{t+x} \left(B_R (u), K\right)\). The result follows. ◻

If the events \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t''\right)\) and \(\xi_{u, e}^{\mathrm{path}}\left(t, t''\right)\) both hold, then by Remarks 4 and 5 all cycles passing through the edge \(e\) in the underlying graph have at least one edge, distinct from \(e\), closed at time \(t''\). Hence \(e\) is a cut edge in the configuration at that time. We therefore define the following so-called “cut edge event".

Definition 10 (\(\xi_{u, e}^{\mathrm{cut}} (t, t'')\)). Let \(G = (V, E)\) be a \(d\)-regular graph. Let \(t'' > t \geq 0\) and \(u, u' \in V\) such that \(e = \{u, u'\} \in E\). The event \(\xi_{u, e}^{\mathrm{cut}} (t, t'')\) is the event that \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t''\right)\) and \(\xi_{u, e}^{\mathrm{path}}\left(t, t''\right)\) both hold.

Before lower bounding the probability of our cut edge events, we will need the following lemma.

Lemma 6. For all \(n \in \mathbb{N}\) sufficiently large, \(\mu > 0\) and \(\Delta > 0\) such that \(\mu \Delta > 2 \log 2\), \[\int_{0}^{\Delta} \frac{\mu e^{\mu x}}{e^{\mu \Delta} - 1} \left(p + (1-p)e^{-\mu x}\right)^{r - 1} \text{d} x \leq e^{-{\mu \Delta}/{2}} + O\left((\log n)^{-3}\right).\]

Proof. Since \(e^{-y}\) is strictly decreasing in \(y\) and for all \(y \geq 0\), \(e^{-y} \leq 1\), we have that, in particular, \(p + (1-p)e^{-\mu x} \leq 1\) for \(x \in [0, {\Delta}/{2}]\), and \(p + (1-p)e^{-\mu x} \leq p + (1-p)e^{-{\mu \Delta}/{2}}\) for \(x \in [{\Delta}/{2}, \Delta]\).

Also, since \(\mu \Delta > 2 \log 2\), then \(p + (1-p) e^{-{\mu \Delta}/{2}} < {(1+p)}/{2} \leq {(1+p_\mathrm{u}(q, d))}/{2}\), which is also strictly less than \(1\). In particular, by our choice of \(r\), namely \(r \mathrel{\vcenter{:}}= \dfrac{3 \log_{d-1} \log n}{\log_{d-1} (2 / (1 + p_\mathrm{u}(q, d)))}\), there exists some \(C > 0\) such that for all \(n\) sufficiently large, \(\left({(1+p_\mathrm{u}(q, d))}/{2}\right)^{r - 1} < C(\log n)^{-3}\).

Therefore, splitting the integral over these two regions of \(x\), we have that, for all \(n\) sufficiently large, \[\begin{align} \int_{0}^{\Delta} \frac{\mu e^{\mu x}}{e^{\mu \Delta} - 1} \left(p + (1-p)e^{-\mu x}\right)^{r - 1} \text{d} x &\leq \frac{1}{e^{{\mu \Delta}/{2}} + 1} + \int_{{\Delta}/{2}}^{\Delta} \frac{\mu e^{\mu x}}{e^{\mu \Delta} - 1} \left(\frac{1+p_\mathrm{u}(q, d)}{2}\right)^{r - 1} \text{d} x \nonumber \\ &\leq e^{-{\mu \Delta}/{2}} + C(\log n)^{-3} \end{align}\] and hence the result follows. ◻

We are now in a position to prove the following result.

Lemma 7. Let \(n \in \mathbb{N}\) be sufficiently large and \(G = (V, E)\) a \(d\)-regular graph on \(n\) vertices. For all \(\mu > 0\), \(t \geq 0\), \(\Delta > 0\) such that \(\mu \Delta > 2 \log 2\), \(e = \{u, u'\} \in E\) such that \(u\) does not lie on a cycle of length less than \(r\) and \(B_R(u)\) has at most one cycle, and \((\eta, v) \in \{0, 1\}^E \times V\), \[\begin{align} &\;\mathbb{P}\left[\xi_{u, e}^{\mathrm{cut}} (t, \tau_e(t, t + \Delta))\mid \mathsf{C}_e\left(t, t + \Delta\right), (\eta_t, X_t) = (\eta, v)\right] \\ \geq&\;\mathbb{P}\left[S_{[t, t + \Delta]} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \left(1 - (2K + 1) e^{-{\mu \Delta}/{2}} - O\left((\log n)^{-3}\right)\right). \end{align}\]

Proof. First note that, for \(n\) sufficiently large, we have that \(R > r > 1\) and hence for all \(0 \leq x \leq \Delta\), the event \(\xi_{u, e}^{\mathrm{cut}} (t, t + x)\) is independent of the clock \(\mathsf{C}_e\) and uniform \(U_e\) associated with the edge \(e\). Indeed, \(\xi_{u, e}^{\mathrm{cut}} (t, t + x)\) is the event that \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right)\) and \(\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right)\) both hold, where \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right)\) is an event about clocks and uniforms during \((t, t + x)\) for edges in \(E\left(B_r (u)\right)\) distinct from \(e\) while \(\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right)\) is an event about the clocks and uniforms up to time \(t + x\) for edges in \(E\left(B_R (u)\right) \setminus E\left(B_r (u)\right)\).

Hence for all \(0 \leq x \leq \Delta\), we have that the events \(\xi_{u, e}^{\mathrm{cut}} (t, t + x)\) and \(\mathsf{C}_e(t, t + \Delta)\) are independent, the events \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right)\) and \(\xi_{u, e}^{\mathrm{path}}\left(t, t + x\right)\) are independent, and lastly the event \(\xi_{u, e}^{\mathrm{acyclic}}\left(t, t + x\right)\) is independent of the state at time \(t\). Therefore, in conjunction with Lemma 3, \[\begin{align} &\;\mathbb{P} \left[\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e(t, t + \Delta)\right) \mid \mathsf{C}_{e}\left(t, t + \Delta\right), (\eta_t, X_t) = (\eta, v)\right] \nonumber \\ =& \int_{0}^{\Delta} \!\! \frac{\mu e^{\mu x}}{e^{\mu \Delta} - 1} \;\mathbb{P}\left[\xi_{u, e}^{\mathrm{cut}} (t, t + x) \mid (\eta_t, X_t) = (\eta, v) \right] \text{d} x \nonumber \\ =& \int_{0}^{\Delta} \!\! \frac{\mu e^{\mu x}}{e^{\mu \Delta} - 1} \;\mathbb{P}\left[\xi_{u, e}^{\mathrm{acyclic}} (t, t+x) \right] \mathbb{P}\left[\xi_{u, e}^{\mathrm{path}}(t, t+x) \mid (\eta_t, X_t) = (\eta, v) \right] \text{d} x. \label{clp:eq1} \end{align}\tag{26}\]

Recall from Definition 9 that \(S_{[t, t + \Delta]} (R, K)\) is the event that \(S_{t+x} (R, K)\) holds for all \(0 \leq x \leq \Delta\). Then, \(\mathbb{P}\left[S_{t+x} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \geq \mathbb{P}\left[S_{[t, t + \Delta]} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right]\).

Together with Lemmas 4 and 5, we have that \[\begin{align} &\;\mathbb{P}\left[\xi_{u, e}^{\mathrm{acyclic}} (t, t+x) \right] \mathbb{P}\left[\xi_{u, e}^{\mathrm{path}}(t, t+x) \mid (\eta_t, X_t) = (\eta, v) \right] \nonumber \\ \geq&\;\mathbb{P}\left[S_{t+x} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \left(1 - 2K \left(p + (1-p)e^{-\mu x}\right)^{{R}/{2}}\right) \left(1 - \left(p + (1-p)e^{-\mu x}\right)^{r-1} \right). \nonumber \\ >&\;\mathbb{P}\left[S_{[t, t + \Delta]} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \left(1 - (2K + 1) \left(p + (1-p)e^{-\mu x}\right)^{r - 1}\right)\label{clp:eq2} \end{align}\tag{27}\] where the last inequality follows from the facts that \((1 - z)(1 - y) > 1 - z - y\) for \(z, y > 0\), and \(R/2 \geq r - 1\) for sufficiently large \(n\).

Hence, by (26 ) and (27 ), we obtain that \[\begin{align} &\;\mathbb{P} \left[\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e\right) \mid \mathsf{C}_{e}\left(t, t + \Delta\right), (\eta_t, X_t) = (\eta, v)\right] \\ \geq& \;\mathbb{P}\left[S_{[t, t + \Delta]} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \left(1 - (2K + 1) \int_{0}^{\Delta} \!\! \frac{\mu e^{\mu x}}{e^{\mu \Delta} - 1} \;\left(p + (1-p)e^{-\mu x}\right)^{r - 1} \;\mathrm{d} x \right). \end{align}\]

The result follows by Lemma 6. ◻

7.3 Proof of Lemma [lem:stationary95lw95bdd]↩︎

We are finally in a position to lower bound the probability of the event \(\xi_{u \not\to u'} (t, t')\).

Proof. Let \(e \mathrel{\vcenter{:}}= \{u, u'\}\). If the edge \(e\) refreshes during \((t, t')\), i.e. the event \(\mathsf{C}_e(t, t')\) holds, and during the last ring \(\tau_e (t, t')\) both \(U_e (\tau_e (t, t')) > p_{\min}\) and \(\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right)\) hold, then \(e\) is a cut edge in the configuration at time \(\tau_e (t, t')\), then \(e\) is refreshed closed at time \(\tau_e (t, t')\) and remains closed until time \(t'\).

Therefore, the events \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\), \(U_e (\tau_e (t, t')) > p_{\min}\) and \(\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right)\) imply the event \(\xi_{u \not\to u'} (t, t')\). Hence, \[\begin{align} &\;\mathbb{P}\left[\xi_{u \not\to u'} (t, t') \mid (\eta_t, X_t) = (\eta, v)\right] \\ \geq&\;\mathbb{P} \left[V_u (t') = u', \mathsf{C}_e(t, t'), U_e (\tau_e (t, t')) > p_{\min}, \xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right) \mid (\eta_t, X_t) = (\eta, v)\right] \end{align}\]

The events \(V_u (t') = u'\), \(\mathsf{C}_e(t, t')\) and \(U_e \left(\tau_e (t, t')\right) > p_{\min}\) concern the ringing of Poisson clocks and sampling of uniforms during \((t, t')\), which are independent of the state at time \(t\). Furthermore, the random variable \(V_u\) is independent of all environment random variables, namely \(U_e\) and \(\mathsf{C}_e\), and in particular is independent of the event \(\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right)\) which only concerns environment random variables. Also, the event \(U_e \left(\tau_e (t, t')\right) > p_{\min}\) is independent of \(\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right)\), since \(\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right)\) by construction is an event about the ringing of clocks and sampling of uniforms associated with edges distinct from \(e\) up to the last ring time \(\tau_e (t, t')\) of \(e\) in the time interval \((t, t')\). We therefore have that, \[\begin{align} &\;\mathbb{P} \left[V_u (t') = u', \mathsf{C}_e(t, t'), U_e (\tau_e (t, t')) > p_{\min}, \xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right) \mid (\eta_t, X_t) = (\eta, v)\right] \\ =&\;\mathbb{P}\left[V_{u}\left(t'\right) = u'\right] \mathbb{P}\left[U_e \left(\tau_e (t, t')\right) > p_{\min}\mid \mathsf{C}_e\left(t,t'\right)\right]\mathbb{P}\left[\mathsf{C}_e\left(t,t'\right)\right] \\ &\;\mathbb{P}\left[\xi_{u, e}^{\mathrm{cut}} \left(t, \tau_e (t, t')\right) \mid \mathsf{C}_e\left(t, t'\right), (\eta_t, X_t) = (\eta, v)\right]. \end{align}\]

Since \(V_{u}\left(t'\right)\) is uniform on the \(d\) neighbours of \(u\), then \(\mathbb{P}\left[V_{u}\left(t'\right) = u'\right] = 1/d\). By the probability of a Poisson clock ringing in a time interval, we have \(\mathbb{P}\left[\mathsf{C}_e\left(t,t'\right)\right] = 1 - e^{-\mu \left(t' - t\right)}\). Also, since \(U_e\) is uniform on \((0, 1)\), \(\mathbb{P}\left[U_e \left(\tau_e(t, t')\right) > p_{\min}\mid \mathsf{C}_e\left(t,t'\right)\right] = 1 - p_{\min}\). Lastly, by Lemma 7, noting that all necessary conditions are met, we have that for \(n\) sufficiently large, \[\begin{align} &\;\mathbb{P}\left[\xi_{u, e}^{\mathrm{cut}} (t, \tau_e)\mid \mathsf{C}_e\left(t, t'\right), (\eta_t, X_t) = (\eta, v)\right] \\ \geq&\;\mathbb{P}\left[S_{[t, t']} \left(R, K\right) \mid \left(\eta_t, X_t\right) = \left(\eta, v\right)\right] \left(1 - \frac{1}{(\log n)^2} - (2K + 1) e^{-\mu (t' - t)/2}\right). \end{align}\]

The result follows by combining all these lower bounds together, noting that \[(1 - e^{-a})(1 - c - (2b + 1) e^{-a/2}) \geq 1 - c - 2(b + 1)e^{-a/2}\] for all \(a, b > 0\) and \(0 < c < 1\). ◻

8 Global properties of random regular graphs: Lemmas [lem:phase195walk95counts] and [lem:sparse95good95set95size]↩︎

In this section we study several properties of random regular graphs, in order to carefully control the number of walks (under some constraints) which avoid small cycles in the graph, namely those of order \({r(q, d, n)}\). In particular, we will prove Lemmas [lem:phase195walk95counts] and [lem:sparse95good95set95size].

8.1 Count of \((h, i)\)-constrained walks in infinite \(d\)-regular and \((d-1)\)-ary trees↩︎

By \(\mathcal{T}\) we denote the infinite \(d\)-regular tree with a designated vertex \(\rho\) as the root. Also, by \(\widetilde{\mathcal{T}}\) we denote the infinite \((d-1)\)-ary tree with a designated vertex \(\widetilde{\rho}\) as the root.

We will study the counts of \((h, i)\)-constrained walks as in Definition 3. We first note the following observation.

Remark 6. By symmetry in the infinite \(d\)-regular tree \(\mathcal{T}\), given any vertex \(x \in V(\mathcal{T})\) the number of walks from \(x\) of length \(h+2i\) which end at a distance \(h\) away from \(x\) is equal to the number of such walks from the root \(\rho\), namely \(\omega_{h,i}\).

We have the following lemma on the number of \((h, i)\)-constrained walks on the \(d\)-regular and \((d-1)\)-ary infinite trees.

Lemma 8. For all \(h, i \in \mathbb{N}\), \(\omega_{h,i} > \widetilde{\omega}_{h,i} = \frac{h+1}{h+i+1} \binom{h + 2i}{h+i} (d-1)^{h+i}\).

Proof. First note that we can embed the infinite \((d-1)\)-ary tree into the infinite \(d\)-regular tree, with the roots of both trees identified. Hence every \((h, i)\)-constrained walk on \(\widetilde{\mathcal{T}}\) is uniquely mapped to an \((h, i)\)-constrained walk on \(\mathcal{T}\) under this embedding. Therefore, \(\omega_{h,i} > \widetilde{\omega}_{h,i}\).

We will next determine \(\widetilde{\omega}_{h,i}\). Starting at the root \(\widetilde{\rho}\) of \(\widetilde{\mathcal{T}}\), to reach depth \(h\) in a walk of length \(h +2i\), exactly \(h + i\) transitions must increase the depth by \(1\), and \(i\) transitions must decrease the depth by \(1\). The number of such walks is equivalent to the number of walks of length \(h + 2i\) on \(\mathbb{Z}\) starting at \(0\) and ending at \(h\) which never become negative. By Bertrand’s Ballot Theorem, the number of such walks is known to be \(\frac{h+1}{h+i+1} \binom{h + 2i}{h+i}\). For each of the \(h + i\) downward transitions we have \(d-1\) choices, since each vertex has exactly \(d-1\) children.

Therefore the total number of such walks is \(\frac{h+1}{h+i+1} \binom{h + 2i}{h+i} (d-1)^{h+i}\). The result follows. ◻

Remark 7. By symmetry in \(\mathcal{T}\) and the fact that in a tree there is a unique path of length \(h\) from the root to a vertex at depth \(h\), we note that given a set \(\mathcal{P}\) of paths from the root of \(\mathcal{T}\) to depth \(h\), the number of \((h, i)\)-constrained walks which follow a path in \(\mathcal{P}\) is \(\frac{|\mathcal{P}|}{d(d-1)^{h-1}} \;\omega_{h,i}\).

Analogous counts hold for \(\widetilde{\mathcal{T}}\).

8.2 Counts of \((h, i)\)-constrained walks avoiding small cycles in regular graph↩︎

Recall that we will view random walks on a regular graph starting at a vertex \(u\) as walks on the non-backtracking walk tree (NBWT) \(\mathcal{T}_u\) starting at the root, as introduced in Section 2. We will use NBWTs as a convenient tool to bound the number of \((h,i)\)-constrained walks from a given vertex in a regular graph that avoid a given subset of vertices, namely vertices which lie on small cycles of length at most \({r(q, d, n)}\).

For any two vertices \(u\) and \(v\) in a graph \(G\), let \(w_{h,i}(v, u)\) be the set of \((h, i)\)-constrained walks in the NBWT \(\mathcal{T}_v\) which pass through a vertex labelled \(u\). The following lemma upper bounds the total number of \((h,i)\)-constrained walks that pass through a vertex labelled \(u\) across all NBWTs of \(G\).

Lemma 9. Let \(h, i \in \mathbb{N}\) and \(u \in V\). Then, \[\sum_{v \in V} |w_{h,i}(v,u)| \leq \frac{1}{2}(h+i+1)(h+3i+2) \;\omega_{h,i}.\]

Proof. Fix a vertex \(v\) and a walk in \(\mathcal{T}_v\) of length \(h + 2i\) from the root to depth \(h\) that passes through a vertex \(\hat{\rho}\) labelled \(u\). Re-rooting the tree \(\mathcal{T}_v\) at \(\hat{\rho}\) results in another labelled tree isomorphic to \(\mathcal{T}_u\), with the image of the walk passing through the new root \(\hat{\rho}\).

In particular, the walks \(w_{h,i}(v,u)\) in \(\mathcal{T}_v\) correspond to walks in \(\mathcal{T}_u\) of length \(h + 2i\) which:

i. pass through the root \(\hat{\rho}\) of \(\mathcal{T}_u\),

ii. start at some vertex \(x \in V(\mathcal{T}_u)\) labelled \(v\),

iii. and end at a vertex distance \(h\) away from \(x\) in the tree \(\mathcal{T}_u\).

Note that a walk in \(w_{h,i}(v,u)\) on \(\mathcal{T}_v\) may correspond to multiple such walks on \(\mathcal{T}_u\), depending on how many times a vertex labelled \(u\) appears on the walk.

Therefore to upper bound \(\sum_{v \in V} |w_{h,i} (v,u)|\), it suffices to count the number of walks of length \(h+2i\) in \(\mathcal{T}_u\) which start at any vertex, pass through the root, and end at a vertex distance \(h\) away from the starting vertex. Let \(\mathcal{W}\) be the set of all such walks in \(\mathcal{T}_u\).

Given \(r \in \mathbb{N}\) and \(x \in V(\mathcal{T}_u)\), let \(\partial B(x,r)\) be the set of vertices in \(\mathcal{T}_u\) which are distance \(r\) away from \(x\). Then \(|\partial B(x, r)| = d(d-1)^{r-1}\) by regularity, and hence the cardinality is invariant of \(x\). Also, given \(t \in \mathbb{N}\) and \(x, y \in V(\mathcal{T}_u)\), let \(\mathcal{W}_{x, t} (y)\) be the set of all walks of length \(h+2i\) in \(\mathcal{T}_u\) starting at \(x\) which pass through \(y\) at step \(t\).

Now, any walk in \(\mathcal{W}\) must start at most a distance \(h+i\) away from the root \(\hat{\rho}\) of \(\mathcal{T}_u\), since they have to pass through the root and end a distance \(h\) from the start. Therefore, for every walk in \(\mathcal{W}\) there exists \(r \in \mathbb{N}\) such that \(0 \leq r \leq h + i\), the walk starts at a vertex \(x \in \partial B(\hat{\rho}, r)\), and the walk passes through the root \(\hat{\rho}\) at step \(t \in \mathbb{N}\) where \(r \leq t \leq h + 2i\). In particular we have that \[\label{eq:pcl1} \mathcal{W}= \bigcup_{r = 0}^{h+i} \bigcup_{t = r}^{h+2i} \bigcup_{x \in \partial B(\hat{\rho}, r)} \mathcal{W}_{x,t}(\hat{\rho}).\tag{28}\]

Let \(\mathcal{W}_{x,t}^r\) be the union of \(\mathcal{W}_{x,t} (y)\) for all \(y \in \partial B(x,r)\). Note that for any two distinct vertices \(y, z \in V(\mathcal{T}_u)\), the sets \(\mathcal{W}_{x,t} (y)\) and \(\mathcal{W}_{x,t} (z)\) are disjoint, since the walks differ at step \(t\). Therefore, \[\label{eq:pcl2} \left|\mathcal{W}_{x,t}^r\right| = \sum_{y \in \partial B(x,r)} \left|\mathcal{W}_{x,t} (y)\right|.\tag{29}\]

By symmetry in \(\mathcal{T}_u\), for any two vertices \(y, z \in \partial B(x,r)\), the number of walks of length \(h+2i\) which start at \(x\) and end at distance \(h\) from \(x\) while passing through \(y\) and \(z\) respectively at step \(t\) is equal, which is to say that \(|\mathcal{W}_{x,t} (y)| = |\mathcal{W}_{x,t} (z)|\). Therefore in conjunction with (29 ) and the fact that \(|\partial B(x, r)| = d(d-1)^{r-1}\), for all \(y \in \partial B(x,r)\) we have that \[\label{eq:pcl3} |\mathcal{W}_{x,t}^r| = d(d-1)^{r-1} |\mathcal{W}_{x,t} (y)|.\tag{30}\]

Furthermore, \(\mathcal{W}_{x,t}^r\) is a subset of all walks of length \(h+2i\) in \(\mathcal{T}_u\) starting at \(x\) which end a distance \(h\) away from \(x\), and hence by Remark 6 we have that \[\label{eq:pcl4} |\mathcal{W}_{x,t}^r| \leq \omega_{h,i}.\tag{31}\]

In particular, if the root \(\hat{\rho}\) is in \(\partial B (x,r)\) then by (30 ) and (31 ) we have that \[\label{eq:pcl5} |\mathcal{W}_{x,t} (\hat{\rho})| \leq \frac{\omega_{h,i}}{d(d-1)^{r-1}}.\tag{32}\]

Hence from (28 ) and (32 ) along with the fact that \(|\partial B(\hat{\rho}, r)| = d(d-1)^{r-1}\), we have that \[\begin{align} |\mathcal{W}| \leq \sum_{r = 0}^{h+i} \sum_{t = r}^{h+2i} \sum_{x \in \partial B(\hat{\rho}, r)} |\mathcal{W}_{x,t} (\hat{\rho})| \leq \sum_{r = 0}^{h+i} \sum_{t = r}^{h+2i} \sum_{x \in \partial B(\hat{\rho}, r)} \frac{\omega_{h,i}}{d(d-1)^{r-1}} = \frac{1}{2}(h+i+1)(h+3i+2) \;\omega_{h,i} \end{align}\] as required. ◻

Recall that a walk on a graph is said to be \(r\)-acyclic (Definition 2) if it never visits a cycle of length less than \({r(q, d, n)}\). We say that a vertex \(u\) is \((h, i)\)-sparse (Definition 4) if at least \(\left(1 - {1}/{(\log n)^3}\right) \omega_{h, i}\) of the \((h, i)\)-constrained walks from \(u\) are \(r\)-acyclic.

The previous lemma, in conjunction with an upper bound on the number of vertices on cycles of length less than \({r(q, d, n)}\), allows us to establish bounds on the number of \((h, i)\)-sparse vertices in a random regular graph, for appropriate values of \(h\) and \(i\). However, we first require upper bounds on the number of vertices contained in such small cycles. Given a random \(d\)-regular graph, by \(\mathcal{C}_k\) for \(k \geq 3\) we will denote the number of cycles of length \(k\). The following result establishes that for \(r, d \geq 3\) such that \(\sqrt{r}(d-1)^{3r/2 - 1} = o(n)\), the number of cycles \(\mathcal{C}_3, \dots, \mathcal{C}_r\) are asymptotically jointly distributed as independent Poisson random variables \(\mathcal{Z}_3, \dots, \mathcal{Z}_r\) with means \(\frac{(d-1)^k}{2k}\) for \(3 \leq k \leq r\).

Theorem 8 (Theorem 11, [16]). Let \(G\) be a random \(d\)-regular graph on \(n\) vertices with cycle counts \(\left(\mathcal{C}_k \colon k \geq 3\right)\). Let \(\left(\mathcal{Z}_k \colon k \geq 3\right)\) be independent Poisson random variables with \(\mathbb{E}\left[\mathcal{Z}_k\right] = \frac{(d-1)^k}{2k}\). For any \(n \geq 1\) and \(r, d \geq 3\), \(||\left(\mathcal{C}_3, \dots, \mathcal{C}_r\right) - \left(\mathcal{Z}_3, \dots, \mathcal{Z}_r\right)||_{\mathrm{TV}} = O\left(\frac{\sqrt{r} (d-1)^{3r/2 - 1}}{n}\right)\).

The following corollary follows immediately from Theorem 8.

Corollary 1. Let \(G\) be a random \(d\)-regular graph on \(n\) vertices with \(\mathcal{X}_{\leq k} (G)\) vertices on cycles of length at most \(k\). For any \(n \geq 1\) and \(r, d \geq 3\) such that \(\sqrt{r}(d-1)^{3r/2 - 1} = o(n)\), with probability \(1 - o(1)\) over \(G \sim \mathcal{G}(d, n)\), \(\mathcal{X}_{\leq r} (G) \leq (\log n) (d-1)^r\).

We are now in a position to prove the main result for this section.

Proof. Since \(r = O(\log \log n)\) then \(\sqrt{r} (d-1)^{3r/2 - 1} = o(n)\) and by Corollary 1, the number of vertices on cycles of length at most \(r\) is, with probability \(1 - o(1)\) over \(G \sim \mathcal{G}(d, n)\), at most \((\log n)(d-1)^r\). Let \(\mathcal{B}\) be the set of vertices contained in such cycles.

Fix \(h, i \in \mathbb{N}\) such that \(h, i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\). Given a vertex \(u \in \mathcal{B}\), for every vertex \(v \in V\) let \(w_{h,i}(v, u)\) be the set of \((h, i)\)-constrained walks from \(v\) which pass through \(u\).

For a vertex \(v\) to be \((h, i)\)-sparse, it suffices that \(\sum_{u \in \mathcal{B}} |w_{h,i}(v,u)| < ({1}/{(\log n)^3}) \omega_{h,i}\), which is to say that the \((h, i)\)-constrained walks from \(v\) which pass through a vertex in \(\mathcal{B}\) only make up at most a \({1}/{(\log n)^3}\) fraction of all possible \((h, i)\)-constrained walks from \(v\).

Define \(V_\textrm{bad} (h, i) \mathrel{\vcenter{:}}= \{v \in V \colon \sum_{u \in \mathcal{B}} |w_{h,i}(v,u)| \geq ({1}/{(\log n)^3}) \omega_{h,i}\}\). By the previous remark, every vertex in \(V \setminus V_\textrm{bad} (h, i)\) is \((h, i)\)-sparse. Hence it suffices to upper bound the size of \(V_\textrm{bad} (h, i)\). By Lemma 9 we have that \[\sum_{u \in \mathcal{B}} \sum_{v \in V} |w_{h,i}(v,u)| \leq |\mathcal{B}| (1/2)(h+i+1)(h+3i+2) \;\omega_{h,i} \leq \left(1 + {8}/{p_{\min}}\right)^2 |\mathcal{B}| (\log_{d - 1} n)^2 \omega_{h,i}\] where the last inequality follows from \(h, i \leq ({4}/{p_{\min}}) \log_{d - 1} n\). We also have the following lower bound, \[\sum_{u \in \mathcal{B}} \sum_{v \in V} |w_{h,i}(v,u)| \geq \sum_{v \in V_\mathrm{bad} (h, i)} \sum_{u \in \mathcal{B}} |w_{h,i}(v,u)| \geq |V_\textrm{bad} (h, i)| ({1}/{(\log n)^3}) \omega_{h,i}\] where the last inequality follows by the definition of \(V_\textrm{bad} (h, i)\).

Therefore \(|V_\textrm{bad} (h, i)| ({1}/{(\log n)^3}) \omega_{h,i} \leq \left(1 + {8}/{p_{\min}}\right)^2 |\mathcal{B}| (\log_{d - 1} n)^2 \omega_{h,i}\) and hence \[|V_\textrm{bad} (h, i)| \leq \left(\frac{1 + {8}/{p_{\min}}}{\log(d - 1)}\right)^2 (d - 1)^r (\log n)^6\] since \(\left|\mathcal{B}\right| \leq (\log n) (d-1)^r\). Let \(V_{\mathrm{bad}} = \bigcup\limits_{h, i \leq ({4}/{p_{\min}}) \log_{d-1} n} V_{\mathrm{bad}} (h, i)\). Then by the upper bound on each \(V_{\mathrm{bad}} (h, i)\), we have that \(\left|V_{\mathrm{bad}}\right| = O\left((d - 1)^r (\log n)^8\right)\). In particular, all the vertices in \(V \setminus V_{\mathrm{bad}}\) are \((h, i)\)-sparse for all \(h, i \leq ({4}/{p_{\min}}) \log_{d - 1} n\) and hence the result follows. ◻

8.3 Proof of Lemma [lem:phase195walk95counts]↩︎

From Lemma [lem:rrg95glbl95geom] we can prove the first of our two main results for this section.

Figure 1: Illustration of the paths on T_j reaching k-roots which are also in the set S. The ball radius \log_{d-1} \log n highlighted in blue illustrates that at suitable depths every vertex is a k-root. Furthermore, the vertices highlighted by \times are not in S.

Proof. By Lemma [lem:rrg95glbl95geom], with probability \(1 - o(1)\) over \(G \sim \mathcal{G}(d, n)\), there exist \(n - O\left(\mathrm{polylog}(n)\right)\) vertices which are \((h, i)\)-sparse for all \(h, i \in \mathbb{N}\) such that \(h, i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\).

Furthermore, by Lemma 2.1 in [7], we have that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), every ball of radius \(R\) in \(G\) contains at most one cycle. Hence, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), both of these properties hold. Fix such a graph \(G\).

Let \(S\) be the set of all vertices in \(G\) which are \((h, i)\)-sparse for all \(h, i \in \mathbb{N}\) such that \(h, i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\).

Fix a vertex \(u \in V\) and consider the subgraph \(B\) induced by the ball \(B_R (u)\). Since there is at most one cycle in this subgraph, deleting \(u\) from \(B\) results in at least \(d - 1\) components; in particular there must be \(d - 2\) rooted trees \(T_1, \dots, T_{d-2}\), whose roots are neighbours of \(u\), call them \(u_1, \dots, u_{d-2}\) respectively. Furthermore, the (at most) unique cycle in \(B\) is entirely contained in \(H \mathrel{\vcenter{:}}= B \setminus \cup_{j = 1}^{d-1} T_j\).

Fix \(1 \leq j \leq d-2\) and \(h, i \in \mathbb{N}\) such that \({R}/{10} \leq h \leq {R}/{5}\) and \(i \leq {R}/{5}\). The sub-tree \(T_j\) corresponding to \(u_j\) has degree \(d - 1\) at the root and degree \(d\) at all other non-leaf vertices. Hence \(T_j\) is a \((d-1)\)-ary tree of depth \(R - 1\). Since \(h + i \leq {2R}/{5} \leq R - 1\) and \((h, i)\)-constrained walks never exceed depth \(h + i\), the number of \((h, i)\)-constrained walks on \(T_j\) is equal to that on the infinite \((d-1)\)-ary tree, namely \(\widetilde{\omega}_{h, i}\).

Consider an \((h, i)\)-constrained walk on \(T_j\). Suppose a vertex in this walk is contained in a cycle of length less than \(r\) in \(G\). Since the walk reaches a depth of at most \(h + i \leq {2R}/{5}\) and the cycle has length at most \(r\) (which is less than \({R}/{10}\) for all \(n\) sufficiently large), then every vertex in the cycle is at distance at most \({R}/{2}\) away from \(u_j\) in \(G\).

Therefore, since \(u_j\) is a neighbour of \(u\), every vertex in the cycle is at distance at most \(({R}/{2}) + 1\) away from \(u\) in \(G\). Therefore the cycle is entirely contained in the subgraph \(B\) induced by the ball \(B_R (u)\). However, this cycle contains at least one vertex in \(T_j\), and the (at most) unique cycle in \(B\) must be entirely contained in \(H\), which is disjoint from \(T_j\). Therefore such a cycle cannot exist, and hence every such walk in \(T_j\) does not visit a vertex on a cycle of length less than \(r\).

Now for every vertex in \(T_j\) at depth greater than \(k = \log_{d-1} \log n\) but no more than \(R - k - 1\), the ball radius \(k\) around such a vertex is does not contain any cycle in \(G\), and hence every such vertex is a \(k\)-root. This is illustrated in Figure 1. Given \(n\) sufficiently large, it follows that for all \(h \in \mathbb{N}\) such that \({R}/{10} \leq h \leq {R}/{5}\), all the vertices in \(T_j\) at depth \(h\) are \(k\)-roots.

Furthermore, since \({R}/{10} \leq h\) and \(R = (1/5) \log_{d-1} n\), we have that \((d-1)^h = \Omega\left(n^{1/50}\right)\). Since \(|V \setminus S| = O\left(\mathrm{polylog}(n)\right)\), it follows that for \(n\) sufficiently large, \(|V \setminus S| < \frac{1}{2} (d-1)^h\). Therefore, at least half of the vertices at depth \(h\) in \(T_j\) are \(k\)-roots and in \(S\), i.e., in \(V_{\mathrm{good}}\).

Hence for all \(h, i \in \mathbb{N}\) satisfying \({R}/{10} \leq h \leq {R}/{5}\) and \(i \leq {R}/{5}\), at least half of the \((h, i)\)-constrained walks from \(u_j\) in \(T_j\) end at a vertex in \(V_{\mathrm{good}}\), and we have shown that such a walk on \(T_j\) never visits a cycle of length less than \(r\). Since there are at least \(d - 2\) such neighbours of \(u\), the result follows. ◻

8.4 Proof of Lemma [lem:sparse95good95set95size]↩︎

For any two vertices \(u, v \in V\) and \(h \in \mathbb{N}\), by \(\mathcal{P}_{h} (u, v)\) we denote the family of paths length \(h\) between \(u\) and \(v\). Recall that \(H_{\min}(d, n)\mathrel{\vcenter{:}}= \left\lfloor\log_{d-1} n\right\rfloor + 2 \left\lfloor\log_{d-1}\log n\right\rfloor\) and \(H_{\max}(d, n)\mathrel{\vcenter{:}}= \left\lfloor\log_{d-1} n\right\rfloor + \left\lfloor (1/10)\log_{d-1} n\right\rfloor - 1\). Before proving Lemma [lem:sparse95good95set95size], we require the following result.

Lemma 10 ([7], Lemma 3.5). With probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), for any \(u, v \in V_{\mathrm{roots}}\) and \(h \in \mathbb{N}\) satisfying \(d(u, v) > 2 k\) and \(h \in [H_{\min}(d, n), H_{\max}(d, n)]\): \(\left|\mathcal{P}_{h} (u, v)\right| \geq \left(\frac{1-o_n(1)}{n}\right) d(d-1)^{h - 1}\).

We also take note of the following useful observations about the counts of paths between \(k\)-roots and the number of \((h, i)\)-constrained walks between them.

Remark 9. Observe that every simple path in a \(d\)-regular graph \(G\) starting at a vertex \(u\) corresponds to a simple path in the non-backtracking walk tree \(\mathcal{T}_u\) starting at the root.

Consequently, for any two \(k\)-roots \(u\) and \(v\) as in Lemma 10 and \(h \in \mathbb{N}\) such that \(H_{\min}(d, n)\leq h \leq H_{\max}(d, n)\), each path in \(\mathcal{P}_{h} (u, v)\) maps to a distinct path in \(\mathcal{T}_u\) starting at the root. By Remark 7, for any \(i \in \mathbb{N}\), the number of \((h, i)\)-constrained walks which follow one of these paths in \(\mathcal{T}_u\) is at least \(\left(1-o(1)\right) (1/n) \omega_{h,i}\).

Lastly, we will require the following result which establishes that, with probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), most vertices are \(k\)-roots.

Theorem 10 ([7], Lemma 3.2). With probability \(1 - o_n(1)\) over \(G \sim \mathcal{G}(d, n)\), there are \(n - o(n)\) vertices in \(G\) which are \(k\)-roots.

We are now in a position to prove Lemma [lem:sparse95good95set95size], which gives an analogue of Lemma 10 for counts of \((h, i)\)-constrained walks avoiding small cycles of length less than \(r\).

Proof. By Lemmas 10 and 10, with probability \(1 - o(1)\) over \(G \sim \mathcal{G}(d, n)\), for a given \(k\)-root \(u\) there exists a set \(S'_u\) of \(k\)-roots with size \(n - o(n)\), such that for all \(h \in \mathbb{N}\) satisfying \(H_{\min}(d, n)\leq h \leq H_{\max}(d, n)\) and for all \(v \in S'_u\): \(\left|\mathcal{P}_{h} (u, v)\right| \geq \left(1-o(1)\right) \left({1}/{n}\right) d(d-1)^{h - 1}\). In particular, for all \(n\) sufficiently large, we have that \(\left|\mathcal{P}_{h} (u, v)\right| \geq ({1}/{2n}) d(d-1)^{h - 1}\).

Fix such a graph \(G\), \(u \in V_{\mathrm{good}}\) and \(v \in S'_u\). For all \(i \in \mathbb{N}\), by Remark 9 we have that the paths in \(\mathcal{P}_h(u,v)\) correspond to \(({1}/{2n}) \omega_{h, i}\) walks from \(u\) which are \((h, i)\)-constrained and end at \(v\).

Let \(\mathcal{W}_{h, i} (u, v)\) correspond to this collection of walks. By the previous remarks, \(\left|\mathcal{W}_{h, i} (u, v)\right| > ({1}/{2n}) \omega_{h, i}\). Clearly for any two distinct vertices \(v, w \in S'_u\), these collections are disjoint (namely since they end at distinct vertices). Let \(\mathcal{W}_{h, i} (u)\) be the collection of all \((h, i)\)-constrained walks from \(u\), where \(\left|\mathcal{W}_{h, i} (u)\right| = \omega_{h, i}\).

Next fix \(h, i \in \mathbb{N}\) such that \(H_{\min}(d, n)\leq h \leq H_{\max}(d, n)\) and \(i \leq \left({4}/{p_{\min}}\right) \log_{d - 1} n\). Since \(u \in V_{\mathrm{good}}\) then it is \((h, i)\)-sparse and hence at most \(\left({1}/{(\log n)^3}\right) \omega_{h, i}\) of the walks in \(\mathcal{W}_{h, i} (u)\) are not \(r\)-acyclic. Let this collection of “bad" walks be \(\mathcal{W}^\mathrm{bad}_{h, i} (u)\), and for all \(v \in S'_u\) let \(\mathcal{W}^\mathrm{bad}_{h, i} (u, v) \mathrel{\vcenter{:}}= \mathcal{W}^\mathrm{bad}_{h, i} (u) \cap \mathcal{W}_{h, i} (u, v)\). Note that all of these sets are pairwise disjoint and that their union is a subset of \(\mathcal{W}^\mathrm{bad}_{h, i} (u)\).

Let \(S^\mathrm{bad}_u (h, i)\) be the subset of \(S'_u\) such that \(v \in S^\mathrm{bad}_u (h, i)\) if, and only if, \(\left|\mathcal{W}^\mathrm{bad}_{h, i} (u, v)\right| > \frac{1}{2} \left|\mathcal{W}_{h, i} (u, v)\right|\). Then, \[\begin{align} \left|S^\mathrm{bad}_u (h, i)\right| \frac{\omega_{h, i}}{4n} \leq \sum_{v \in S^\mathrm{bad}_u (h, i)} \frac{\left|\mathcal{W}_{h, i} (u, v)\right|}{2} < \sum_{v \in S^\mathrm{bad}_u (h, i)} \left|\mathcal{W}^\mathrm{bad}_{h, i} (u, v)\right| \leq \left|\mathcal{W}^\mathrm{bad}_{h, i} (u)\right| < \frac{\omega_{h, i}}{(\log n)^3} \end{align}\] and hence \(\left|S^\mathrm{bad}_u (h, i)\right| = O\left({n}/{(\log n)^3}\right)\).

Define \(S_u \mathrel{\vcenter{:}}= S'_u \setminus \bigcup\limits_{h, i} S^\mathrm{bad}_u (h, i)\). Since \(\left|S^\mathrm{bad}_u (h, i)\right| = O\left({n}/{(\log n)^3}\right)\) for all \(H_{\min}(d, n)\leq h \leq h_{\max}\) and \(i \leq ({4}/{p_{\min}}) \log_{d - 1} n\), then \(|S_u| \geq |S'_u| - O\left({n}/{(\log n)}\right)\). Since \(|S'_u| = n - o(n)\) it follows that \(|S_u| = n - o(n)\).

We have therefore constructed a set \(S_u\) of size \(n - o(n)\), such that for all \(H_{\min}(d, n)\leq h \leq h_{\max}\) and \(i \leq ({4}/{p_{\min}}) \log_{d - 1} n\) and for all \(v \in S_u\), there are at least \(({1}/{4 n}) \omega_{h, i}\) walks from \(u\) to \(v\) which are \((h, i)\)-constrained and \(r\)-acyclic. The result follows. ◻

9 Proof of Lemma [lem:simplex95integral]↩︎

Proof. Fix \(c > 0\) and \(c' > c\). We will eventually choose an appropriate value of \(c\). Let \(\delta \mathrel{\vcenter{:}}= {c'}/{(\log n)}\). Let \(x_0, \dots, x_\alpha \in [0, 1]\) be the shifted and rescaled variables satisfying \(x_i = \frac{T_i - \delta}{T - (\alpha + 1)\delta}\) for all \(0 \leq i \leq \alpha\), and \(\sum_{i=0}^\alpha x_i = 1\). The Jacobian for this change of variables satisfies \(\mathrm{d} \overrightarrow{T_\alpha} = \left(T - (\alpha + 1)\delta\right)^\alpha \;\mathrm{d} \overrightarrow{x}\), and the simplex is bijectively mapped onto the unit simplex \[\mathcal{K}_{\alpha} \mathrel{\vcenter{:}}= \left\{x_0, \dots, x_\alpha \in [0, 1] \colon x_0 + \dots + x_\alpha = 1 \right\}\] with volume \(\mathrm{Vol}(\mathcal{K}_\alpha) = \frac{1}{\alpha!}\). Let \(\mu_0 \mathrel{\vcenter{:}}= \varepsilon \log n\) and \(\lambda \mathrel{\vcenter{:}}= \frac{1}{2} \mu_0 \left(T - (\alpha + 1) \delta\right)\). Then, for all \(0 \leq i \leq \alpha\), we have that \(1 - C n^{-(\varepsilon / 2) T_i} = 1 - C \exp\left(-(\mu_0 / 2) T_i\right) = 1 - C \exp\left(-(\mu_0 / 2) \delta - \lambda x_i\right)\), and therefore \[\label{eq:sil1} \int_{\substack{T_0, \dots, T_\alpha \geq \delta \\ \sum_{i = 0}^\alpha T_i = T}} \;\prod_{i=0}^\alpha \left(1 - C n^{-(\varepsilon / 2) T_i}\right) \;\mathrm{d} \overrightarrow{T_\alpha} = \left(T - (\alpha + 1)\delta\right)^\alpha \int_{\mathcal{K}_{\alpha}} \prod_{i = 0}^\alpha \left(1 - C e^{-(\mu_0 / 2) \delta - \lambda x_i}\right) \;\mathrm{d} \overrightarrow{x}.\tag{33}\]

Let \(\gamma(c) \mathrel{\vcenter{:}}= C e^{-c \left({\varepsilon}/{2}\right)}\). Since \(\mu_0 = \varepsilon \log n\) and \(\delta = \frac{c'}{\log n} > \frac{c}{\log n}\), then, for all \(0 \leq i \leq \alpha\), we have that \(1 - C e^{-(\mu_0 / 2) \delta - \lambda x_i} \geq 1 - \gamma(c) e^{-\lambda x_i}\). By inclusion-exclusion, we also have that \[\label{eq:sil2} \prod_{i = 0}^\alpha \left(1 - \gamma(c) e^{-\lambda x_i}\right) = \sum_{j = 0}^{\alpha+1} (-1)^j \left(\gamma(c)\right)^j \sum_{\substack{S \subseteq \{0, 1, \dots, \alpha\} \\ |S| = j}} e^{-\lambda \sum_{i\in S} x_i}\tag{34}\]

By symmetry over the simplex, we have that each subset of size \(j\) contributes equally to the summation above, and therefore we can rewrite (34 ) as \[\label{eq:sil3} \prod_{i = 0}^\alpha \left(1 - \gamma(c) e^{-\lambda x_i}\right) = \sum_{j = 0}^{\alpha+1} (-1)^j \left(\gamma(c)\right)^j \binom{\alpha + 1}{j} e^{-\lambda \sum_{i = 1}^j x_{i-1}}.\tag{35}\]

Letting \(M_j (\lambda) \mathrel{\vcenter{:}}= \left(\gamma(c)\right)^j \binom{\alpha + 1}{j} \int_{\mathcal{K}_{\alpha}} e^{-\lambda \sum_{i = 1}^j x_{i-1}} \;\mathrm{d} \overrightarrow{x}\) for \(0 \leq j \leq \alpha + 1\), we therefore have that \[\label{eq:sil4} \int_{\mathcal{K}_{\alpha}} \prod_{i = 0}^\alpha \left(1 - C e^{-(\mu_0 / 2) \delta - \lambda x_i}\right) \;\mathrm{d} \overrightarrow{x} \geq \sum_{j = 0}^{\alpha+1} (-1)^j M_j(\lambda).\tag{36}\]

We will proceed to show that for \(j \geq 1\), the differences \(M_j (\lambda) - M_{j+1} (\lambda)\) are non-negative. For \(1 \leq j \leq \alpha\), let \(X_j = x_0 + \dots + x_{j-1}\). Under the uniform distribution on \(\mathcal{K}_{\alpha}\), \((x_0, \dots, x_\alpha)\) follow a \(\mathrm{Dirichlet}(1, \dots, 1)\) distribution, and by the aggregation property for the Dirichlet distribution, \((X_j, 1 - X_j)\) follows a \(\mathrm{Dirichlet}(j, \alpha - j + 1)\) distribution. From the marginal properties of Dirichlet distributions, we therefore have that \(X_j\) follows a \(\mathrm{Beta}(j, \alpha - j + 1)\) distribution with probability density function \[\begin{align} f_{X_j}(x) = \frac{\alpha!}{(j-1)! \;(\alpha - j)!} x^{j - 1}(1 - x)^{\alpha - j}, & \forall x \in [0, 1]. \end{align}\]

We therefore have that \[M_j (\lambda) = \mathrm{Vol}(\mathcal{K}_{\alpha}) \cdot \left(\gamma(c)\right)^j \binom{\alpha + 1}{j} \cdot \frac{\alpha!}{(j-1)! \;(\alpha - j)!} \int_0^1 e^{-\lambda x} x^{j - 1} (1 - x)^{\alpha - j} \;\mathrm{d} x.\]

Next note that since \((1-x)^{\alpha-j} \leq 1\), then \[\int_0^1 e^{-\lambda x} x^{j - 1}(1 - x)^{\alpha - j} \;\mathrm{d} x \leq \int_0^1 e^{-\lambda x} x^{j - 1} \;\mathrm{d} x = \lambda^{-j} \bigl((j-1)! - \Gamma(j, \lambda)\bigr) \leq \lambda^{-j} (j-1)!\] where \(\Gamma(a, z)\) denotes the incomplete gamma function.

Hence \(M_j (\lambda)\) is upper bounded by \(\left({\gamma(c)}/{\lambda}\right)^j \binom{\alpha + 1}{j} \frac{1}{(\alpha-j)!}\). We next seek a lower bound on \(M_j (\lambda)\). First note that \(\lambda = \Omega\left((\log n)^2\right)\) and that \(j = O(\log n)\). Hence for all \(n\) sufficiently large, \(j \lambda^{-1} \in (0, 1)\) and therefore for \(x \in \left[0, j \lambda^{-1}\right]\) it follows that \(1 - x \geq 1 - j \lambda^{-1}\).

We therefore have that, for all \(n\) sufficiently large, \[\begin{align} \int_0^1 e^{-\lambda x} x^{j - 1}(1 - x)^{\alpha - j} \;\mathrm{d} x &\geq \int_0^{j \lambda^{-1}} e^{-\lambda x} x^{j - 1} \left(1 - j \lambda^{-1}\right)^{\alpha - j} \;\mathrm{d} x \\ &= \left(1 - j \lambda^{-1}\right)^{\alpha-j} \lambda^{-j} \bigl((j-1)! - \Gamma(j, j)\bigr). \end{align}\]

For all \(n\) sufficiently large, we have that \(\lambda > \frac{\varepsilon}{200 p_{\min}} (\log_{d-1} n) (\log n)\). Together with the fact that \(j \leq \alpha \leq ({4}/{p_{\min}}) \log_{d - 1} n\), we have that for sufficiently large \(n\), \[\left(1-j \lambda^{-1}\right)^{\alpha - j} \geq \left(1 - \frac{\left({800}/{\varepsilon}\right)}{\log n}\right)^{({4}/{p_{\min}}) \log_{d-1} n} \geq e^{-\frac{6400}{\varepsilon p_{\min}}}.\]

We also have that the term \((j-1)! - \Gamma(j, j)\) is lower bounded by \(\frac{(j-1)!}{2}\). Let \(\theta \mathrel{\vcenter{:}}= \frac{1}{2} e^{-\frac{6400}{\varepsilon p_{\min}}}\). Then \(M_j (\lambda)\) is lower bounded by \(\theta \left({\gamma(c)}/{\lambda}\right)^j\binom{\alpha + 1}{j}\frac{1}{(\alpha - j)!}\). Next note that, for all \(n\) sufficiently large, \[\begin{align} M_j (\lambda) - M_{j+1}(\lambda) &\geq \frac{\theta}{(\alpha-j)!} \left(\frac{\gamma(c)}{\lambda}\right)^j \binom{\alpha + 1}{j} - \left(\frac{\gamma(c)}{\lambda}\right)^{j+1} \frac{1}{(\alpha-j-1)!} \binom{\alpha + 1}{j+1} \nonumber \\ &= \frac{1}{(\alpha-j)!} \left(\frac{\gamma(c)}{\lambda}\right)^j \binom{\alpha + 1}{j} \left(\theta - \frac{\gamma(c) (\alpha-j)(\alpha-j+1)}{\lambda(j+1)}\right) \nonumber \\ &\geq \frac{1}{(\alpha-j)!} \left(\frac{\gamma(c)}{\lambda}\right)^j \binom{\alpha + 1}{j} \left(\theta - \gamma(c) \left(\frac{1600 ({4}/{p_{\min}} + {1}/{(\log n)})}{\varepsilon}\right)\right) \nonumber \\ &\geq \frac{1}{(\alpha-j)!} \left(\frac{\gamma(c)}{\lambda}\right)^j \binom{\alpha + 1}{j} \left(\theta - \gamma(c) \left(\frac{12800}{\varepsilon p_{\min}}\right)\right) \label{eq:sil5}. \end{align}\tag{37}\] and hence for \(M_j (\lambda) - M_{j+1}(\lambda) \geq 0\) we require that \(\theta \geq \gamma(c) \left(\frac{12800}{\varepsilon p_{\min}}\right)\).

Substituting for \(\theta\) and \(\gamma(c)\) then re-arranging, we require that \[\label{eq:sil6} e^{c \left({\varepsilon}/{2}\right)} \geq C \left(\frac{25600}{\varepsilon p_{\min}}\right)e^{\frac{6400}{\varepsilon p_{\min}}}.\tag{38}\]

In particular, there exists \(c > 0\) such that (38 ) is satisfied and hence \(M_j (\lambda) > M_{j+1}(\lambda)\) for all \(1 \leq j \leq \alpha - 1\) by (37 ). Therefore, for this choice of \(c\), together with (36 ) we have that \[\label{eq:sil7} \int_{\mathcal{K}_{\alpha}} \prod_{i = 0}^\alpha \left(1 - C e^{-(\mu_0 / 2) \delta - \lambda x_i}\right) \;\mathrm{d} \overrightarrow{x} \geq M_0 (\lambda) - M_1 (\lambda) - M_{\alpha + 1} (\lambda).\tag{39}\]

Now, \(M_0 (\lambda) = \frac{1}{\alpha!}\), \(M_1 (\lambda) \leq \frac{\gamma(c)}{\lambda} \frac{\alpha + 1}{(\alpha - 1)!} = \frac{1}{\alpha!} \left(\gamma(c)\lambda^{-1} \alpha (\alpha + 1)\right)\) and \(M_{\alpha + 1} (\lambda) = \frac{1}{\alpha!} \frac{\left(\gamma(c)\right)^{\alpha + 1}}{e^\lambda}\). Since \(\gamma(c)\) is a constant independent of \(n\), \(\alpha = \Theta(\log n)\) and \(\lambda = \Omega\left((\log n)^2\right)\), it follows that \(M_{\alpha + 1} (\lambda) = \frac{o(1)}{\alpha!}\). In particular then, for \(n\) sufficiently large, we have that \[\begin{align} M_0 (\lambda) - M_1 (\lambda) - M_{\alpha + 1} (\lambda) &\geq \frac{1}{\alpha!} \left(1 - \frac{\gamma(c) \alpha (\alpha + 1)}{\lambda} - o_n(1)\right) \nonumber \\ &\geq \frac{1}{\alpha!} \left(1 - o_n(1) - \gamma(c) \left(\frac{12800}{\varepsilon p_{\min}}\right)\right) \nonumber \\ &\geq \frac{1}{\alpha!} \left(\frac{1}{4} + \left(\theta - \gamma(c) \left(\frac{12800}{\varepsilon p_{\min}}\right)\right)\right) \nonumber \\ &\geq \frac{1}{4} \frac{1}{\alpha!} \label{eq:sil8} \end{align}\tag{40}\] where the first inequality follows from our bounds on each term, the second inequality follows from \(\alpha \leq ({4}/{p_{\min}}) \log_{d-1} n\) and \(\lambda > \frac{\varepsilon}{200 p_{\min}} (\log_{d-1} n)(\log n)\), the third inequality follows from the fact that \(\theta < \frac{1}{2} < \frac{3}{4} - o_n(1)\) for \(n\) sufficiently large, and the last inequality follows from our previous remark that our choice of \(c\) is such that \(\theta \geq \gamma(c) \left(\frac{12800}{\varepsilon p_{\min}}\right)\). The result follows from (33 ), (39 ) and (40 ). ◻

10 Proofs of Lemmas [lem:f95alpha95properties] and [lem:poisson95concentration]↩︎

Proof. First note that \(\frac{\partial}{\partial x} f_\alpha(x, y) = \log\left(\frac{a(\alpha - x - y)}{x(1-a-b)}\right)\) and \(\frac{\partial}{\partial y} f_\alpha(x, y) = \log\left(\frac{b(\alpha - x - y)}{y(1-a-b)}\right)\). Both of these are zero at \((x, y) = \left(a \alpha, b \alpha\right)\). Hence (i) follows.

Let \(c = 1 - a - b\) and define \(h_\alpha(x, y) \mathrel{\vcenter{:}}= f_\alpha\left(a \alpha, b \alpha\right) - f_\alpha\left(a (\alpha + x), b (\alpha + y)\right)\). We will show that, given \(\alpha\) sufficiently large, for all \(0 \leq x, y \leq \sqrt{\alpha}\) we have that \(h_\alpha (x, y) \leq h_\alpha (\sqrt{\alpha}, \sqrt{\alpha})\). It suffices to show that \(\left(\frac{\partial}{\partial x} h_\alpha\right) (x, y) \geq 0\) and \(\left(\frac{\partial}{\partial y} h_\alpha\right) (x, y) \geq 0\) over the region \(0 \leq x, y \leq \sqrt{\alpha}\). First, observe that, for all \(\alpha\) sufficiently large, we have that for all \(0 \leq x, y \leq \sqrt{\alpha}\):

i. \(\left(\frac{\partial^2}{\partial x \partial x} h_\alpha\right) (x, y) = \frac{a}{\alpha + x} + \frac{a^2}{c \alpha - a x - by} > 0\) and \(\left(\frac{\partial^2}{\partial y \partial y} h_\alpha\right) (x, y) = \frac{b}{\alpha + y} + \frac{b^2}{c \alpha - a x - by} > 0\);

ii. \(\left(\frac{\partial^2}{\partial x \partial y} h_\alpha\right) (x, y) = \frac{a b}{c \alpha - a x - by} > 0\).

Also, \[\begin{align} \left(\frac{\partial}{\partial x} h_\alpha\right) (x, y) &= a \left( \log\left(\frac{c}{c \alpha - a x - b y}\right) - \log\left(\frac{1}{\alpha + x}\right)\right) \\ \mathrm{and} \quad \left(\frac{\partial}{\partial y} h_\alpha\right) (x, y) &= b \left(\log\left(\frac{c}{c \alpha - a x - b y}\right) - \log\left(\frac{1}{\alpha + y}\right)\right) \end{align}\] which are both \(0\) at \((x, y) = (0, 0)\). Hence, in conjunction with the second-order derivatives, we have that for \(\alpha\) sufficiently large then \(\left(\frac{\partial}{\partial x} h_\alpha\right) (x, y) \geq 0\) and \(\left(\frac{\partial}{\partial y} h_\alpha\right) (x, y) \geq 0\) over the region \(0 \leq x, y \leq \sqrt{\alpha}\), and consequently \(h_\alpha (x, y) \leq h_\alpha (\sqrt{\alpha}, \sqrt{\alpha})\) over this region.

All that remains to prove (ii) is to show that for \(\alpha\) sufficiently large, \(h_\alpha \left(\sqrt{\alpha}, \sqrt{\alpha}\right) \leq \frac{a+b}{2(1-a-b)}\). Consider the function \[g(z) = \frac{(1-c)^2 \left((1-c + cz) \log \left(1+\frac{c z}{1-c}\right)+ c(1-z) \log (1-z)\right)}{c^2 z^2}\] where, after some simple algebraic manipulation, we have that \[\begin{align} h_\alpha (\sqrt{\alpha}, \sqrt{\alpha}) &= \sqrt{\alpha} \left(\left(1+\sqrt{\alpha}\right) (1-c) \log \left(1+\frac{1}{\sqrt{\alpha}}\right)+\left(c\sqrt{\alpha} + c - 1\right) \log \left(1-\frac{1-c}{c \sqrt{\alpha}}\right)\right)\\ &= g\left(\frac{1-c}{c \sqrt{\alpha}}\right). \end{align}\]

From the Taylor expansion of \(g(z)\), we have that \(g(z) = \frac{1-c}{2c} + O(z)\) and therefore \(\lim\limits_{\alpha \to \infty} g\left(\frac{1-c}{c \sqrt{\alpha}}\right) = \frac{1-c}{2c}\). We also have that the Taylor expansion of \(g'(z)\) is given by \[g'(z) = \frac{1-2c}{6c} + \frac{(1-c)^2}{c}\sum_{k = 1}^\infty \frac{k-1}{k(k+1)}\left(1-\frac{c^{k+2}}{(c-1)^{k+2}}\right)z^k\] which one can readily check is negative for \(z\) sufficiently small given \(\frac{1}{2} < c < 1\). In particular there exists \(\alpha_0 > 0\) such that for all \(\alpha > \alpha_0\) the function \(g\left(\frac{1-c}{c \sqrt{\alpha}}\right)\) is increasing. Therefore for all \(\alpha > \alpha_0\), \(g\left(\frac{1-c}{c \sqrt{\alpha}}\right) < \frac{1-c}{2c}\) and hence (ii) follows. ◻

Proof. For any \(Z \sim \mathsf{Poisson}(T)\) for some \(T > 0\), it is well known (see e.g. Proposition 2.10 in [17]) that for any \(x > 0\), \(\mathbb{P}\left[\left|Z - T\right| \geq x\right] \leq 2 \exp\left(-\frac{x^2}{2(T + x)}\right)\). Fix \(T = \frac{k_{\min}+ k_{\max}}{2}\) and \(x = \frac{k_{\max}- k_{\min}}{2}\); observe that \(C_{\min}\log n \leq T \leq C_{\max}\log n\) and \(\left({C_{\mathrm{mid}}}/{2}\right) \log n \leq x \leq \left({C_{\max}}/{2}\right) \log n\).

Letting \(C \mathrel{\vcenter{:}}= \frac{C_{\mathrm{mid}}^2}{12 C_{\max}}\), since \(Y\) has rate 1 then \(N_Y(T)\) is distributed as \(\mathsf{Poisson}(T)\) and hence \[\mathbb{P} \left[k_{\min}\leq N_Y (T) \leq k_{\max}\right] \geq 1 - \mathbb{P}\left[\left|N_Y (T) - T\right| \geq x\right] \geq 1 - 2n^{-C} \geq \frac{1}{2}\] where the last inequality follows holds for all \(n\) sufficiently large, as required. ◻

References↩︎

[1]
A. Blanca and R. Gheissari, Random-Cluster Dynamics on Random Regular Graphs in Tree Uniqueness,” Communications in Mathematical Physics, vol. 386, no. 2, pp. 1243–1287, 2021, [Online]. Available: https://link.springer.com/content/pdf/10.1007/s00220-021-04093-z.pdf.
[2]
Y. Peres, A. Stauffer, and J. E. Steif, Random walks on dynamical percolation: mixing times, mean squared displacement and hitting times,” Probability Theory and Related Fields, vol. 162, no. 3, pp. 487–530, Aug. 2015, doi: 10.1007/s00440-014-0578-4.
[3]
Y. Peres, P. Sousi, and J. E. Steif, Mixing time for random walk on supercritical dynamical percolation,” Probability Theory and Related Fields, vol. 176, no. 3, pp. 809–849, Apr. 2020, doi: 10.1007/s00440-019-00927-z.
[4]
P. Sousi and S. Thomas, Cutoff for random walk on dynamical Erdős–Rényi graph,” Annales de l’Institut Henri Poincaré, Probabilités et Statistiques, vol. 56, no. 4, Nov. 2020, doi: 10.1214/20-aihp1057.
[5]
J. Hermon and P. Sousi, A comparison principle for random walk on dynamical percolation,” The Annals of Probability, vol. 48, no. 6, Nov. 2020, doi: 10.1214/20-aop1441.
[6]
A. Lelli and A. Stauffer, “Mixing time of random walk on dynamical random cluster,” Probability Theory and Related Fields, vol. 189, no. 3, pp. 981–1043, Aug. 2024, doi: 10.1007/s00440-024-01262-8.
[7]
E. Lubetzky and A. Sly, Cutoff phenomena for random walks on random regular graphs,” Duke Mathematical Journal, vol. 153, no. 3, Jun. 2010, doi: 10.1215/00127094-2010-029.
[8]
T. Sauerwald and L. Zanetti, Random Walks on Dynamic Graphs: Mixing Times, Hitting Times, and Return Probabilities,” in 46th international colloquium on automata, languages, and programming (ICALP 2019), 2019, vol. 132, pp. 93:1–93:15, doi: 10.4230/LIPIcs.ICALP.2019.93.
[9]
L. Cai, T. Sauerwald, and L. Zanetti, “Random walks on randomly evolving graphs,” in Structural information and communication complexity, 2020, pp. 111–128.
[10]
S. Andres, N. Gantert, D. Schmid, and P. Sousi, Biased random walk on dynamical percolation,” The Annals of Probability, vol. 52, no. 6, Nov. 2024, doi: 10.1214/23-aop1679.
[11]
L. Avena, R. Baldasso, R. S. Hazra, F. den Hollander, and M. Quattropani, arXiv:2501.08703The voter model on random regular graphs with random rewiring.” 2025, [Online]. Available: https://arxiv.org/abs/2501.08703.
[12]
L. Avena, H. Güldaş, R. van der Hofstad, and F. den Hollander, “MIXING TIMES OF RANDOM WALKS ON DYNAMIC CONFIGURATION MODELS,” The Annals of Applied Probability, vol. 28, no. 4, pp. 1977–2002, 2018, Accessed: Sep. 17, 2025. [Online]. Available: https://www.jstor.org/stable/26542447.
[13]
L. Avena, H. Güldaş, R. van der Hofstad, and F. den Hollander, Random walks on dynamic configuration models: A trichotomy,” Stochastic Processes and their Applications, vol. 129, no. 9, pp. 3360–3375, 2019, doi: https://doi.org/10.1016/j.spa.2018.09.010.
[14]
L. Avena, R. van der Hofstad, F. den Hollander, and O. Nagy, “Mixing of fast random walks on dynamic random permutations,” Probability Theory and Related Fields, Apr. 2025, doi: 10.1007/s00440-025-01375-8.
[15]
M. Biskup and P.-F. Rodriguez, “Limit theory for random walks in degenerate time-dependent random environments,” Journal of Functional Analysis, vol. 274, no. 4, pp. 985–1046, 2018, doi: https://doi.org/10.1016/j.jfa.2017.12.002.
[16]
T. Johnson, Exchangeable Pairs, Switchings, and Random Regular Graphs,” The Electronic Journal of Combinatorics, vol. 22, no. 1, Feb. 2015, doi: 10.37236/4659.
[17]
M. Wainwright, High-dimensional statistics: a non-asymptotic viewpoint. Cambridge University Press, 2019.

  1. Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK↩︎

  2. If we omit the cut edge analysis and allow the effective random cluster parameters to differ so that some are \(p_{\min}\) and others are \(p_{\max}\), then the lower bound on the probability in Lemma 1 becomes smaller by a polynomial in \(n\), which leads to a corresponding factor in the coupling time.↩︎

  3. We will in fact assume that, whenever the walker clock \(\mathsf{C}_{\mathrm{w}}\) rings (say at time \(t\)), for all \(u\in V\) we sample \(V_u(t)\) to be a uniformly chosen neighbour of \(u\) in \(G\); the chain only uses \(V_{X_{t-}}(t)\) but the extra samples will be technically convenient later.↩︎