This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

A note on the Fredrickson-Andersen one spin facilitated model in stationarity

Abstract.

This note discusses three problems related to the Fredrickson-Andersen one spin facilitated model in stationarity. The first, considered in 2008 in a paper of Cancrini, Martinelli, Roberto and Toninelli, is the spectral gap of the model’s infinitesimal generator. They study the decay of this spectral gap when the density is large, but in dimensions 33 and higher, they do not find the exact exponent. They also show that the persistence function of the model has exponential tail, but the typical decay time is not analyzed. We will see that the correct exponent for the decay of the spectral gap in dimension 33 and higher is 22, and discover how the time over which the persistence function decays diverges in high densities. We also discuss the scaling of the spectral gap in finite graphs.

The author acknowledges the support of the ERC Starting Grant 680275 MALIG

1. Introduction and results

The purpose of this note is to present three small results on the Fredrickson-Andersen one spin facilitated model (FA1f), following [2], addressing two problems that have not been solved there and one tightly related problem in a slightly different setting studied in [4, 5, 6]. Since it is, in a sense, an extension of [2], the reader is referred to [2] for the relevant background, references, and complete introduction of the model and notation.

We will only briefly remind here that sites in d\mathbb{Z}^{d} could be either occupied or empty, with equilibrium probabilities 1q1-q and qq respectively (where qq is thought of as small). When a site has at least one empty neighbor, it is being resampled from equilibrium with rate 11, and otherwise its occupation cannot change. The process is reversible with respect to the invariant measure μ\mu, given by an independent product of Bernoulli random variables with parameter 1q1-q. Probabilities and expected values with respect to the stochastic process are denoted μ\mathbb{P}_{\mu} and 𝔼μ\mathbb{E}_{\mu}, where the subscript μ\mu indicates that the initial configuration is drawn from equilibrium.

The first result here completes Theorem 6.4 of [2], which bounds the spectral gap of the FA1f model. It is shown there that the gap decays polynomially as the parameter qq tends to 0, and for dimensions 11 and 22 the exact exponent is identified, up to a logarithmic correction in dimension 22. For dimension d3d\geq 3, however, the exponent is bounded between 1+2/d1+2/d and 22, and its exact value is not determined. The following theorem shows that the correct scaling is q2q^{2}

Theorem 1.1.

Consider the setting of [2, Theorem 6.4], in dimension d3d\geq 3. Then there exists a positive constant CC (possibly depending on dd) such that

gap()Cq2.\operatorname{gap}(\mathcal{L})\leq C\,q^{2}.

The second result presented here concerns with the persistence function. Recall that

F0(t)\displaystyle F_{0}(t) =μ(τ0>t),\displaystyle=\mathbb{P}_{\mu}(\tau_{0}>t),
τ0\displaystyle\tau_{0} =inf{t:origin is empty at time t}.\displaystyle=\inf\{t:\text{origin is empty at time }t\}.

In general, when the spectral gap is positive, F0(t)et/τ¯F_{0}(t)\leq e^{-t/\overline{\tau}}, where τ¯\overline{\tau} could be chosen to be equal 1gapq\frac{1}{\text{gap}\cdot q}. We will see that for the FA11f model on d\mathbb{Z}^{d} this choice is not optimal, and that the typical time to empty the origin scales (perhaps with lower order corrections) like the inverse of the spectral gap, without the additional factor of qq.

Theorem 1.2.

Consider the FA11f model. Then there exists a positive constant CC such that

F0(t)\displaystyle F_{0}(t) eCq3t\displaystyle\leq e^{-Cq^{3}\,t}\qquad d=1,\displaystyle d=1,
F0(t)\displaystyle F_{0}(t) eCq2log(1/q)t\displaystyle\leq e^{-C\frac{q^{2}}{\log(1/q)}\,t}\qquad d=2,\displaystyle d=2,
F0(t)\displaystyle F_{0}(t) eCq2t\displaystyle\leq e^{-Cq^{2}\,t}\qquad d3.\displaystyle d\geq 3.

Moreover, for a different positive constant CC,

𝔼μ(τ0)\displaystyle\mathbb{E}_{\mu}(\tau_{0}) Cq3\displaystyle\geq Cq^{-3} d=1,\displaystyle d=1,
𝔼μ(τ0)\displaystyle\mathbb{E}_{\mu}(\tau_{0}) Cq2log(1/q)\displaystyle\geq Cq^{-2}\log(1/q)\qquad d=2,\displaystyle d=2,
𝔼μ(τ0)\displaystyle\mathbb{E}_{\mu}(\tau_{0}) Cq2\displaystyle\geq Cq^{-2} d=3.\displaystyle d=3.

The last theorem that will be presented is in a slightly different setting – the Fredrickson-Andersen model on a finite graph GG. A particularly interesting case, studied in [5, 6] and more recently in [4], is when |V(G)|=cq1\left|V(G)\right|=cq^{-1} for some positive constant cc, where V(G)V(G) denotes the set of vertices of GG. A lower bound on the spectral gap is given in [5, 6] and later on in [4] by suggesting a relaxation mechanism in which vacancies travel as random walkers on GG. The next theorem will bound the spectral gap of this model from above, showing that this mechanism has a leading contribution. Consider two independent continuous time random walks on GG, namely, each of the two walkers moves to each neighboring site with rate 11. For two vertices x,yV(G)x,y\in V(G), we define τmeet(x,y)\tau_{\operatorname{meet}}(x,y) to be the expected time that it takes for two such random walkers starting at xx and yy to reach distance at most 11. Let τ¯meet\overline{\tau}_{\operatorname{meet}} be its expected value, starting at two random positions, i.e., τ¯meet=|V(G)|2x,yV(G)τmeet(x,y)\overline{\tau}_{\operatorname{meet}}=\left|V(G)\right|^{-2}\sum_{x,y\in V(G)}\tau_{\operatorname{meet}}(x,y).

Theorem 1.3.

Consider the FA11f model on a finite graph GG, and assume that |V(G)|=c/q\left|V(G)\right|=c/q for some positive constant cc. Let gap(G)\operatorname{gap}(\mathcal{L}_{G}) denote its spectral gap with respect to the product measure conditioned on having at least one vacancy. Then

gap(G)Cqτ¯meet\operatorname{gap}(\mathcal{L}_{G})\leq\frac{Cq}{\overline{\tau}_{\operatorname{meet}}}

for some positive constant CC.

Remark 1.4.

[4] give a lower bound on the spectral gap, which in various graphs is of the same order as the upper bound given in Theorem 1.3. In particular, on the two dimensional torus 𝕋2=2/2\mathbb{T}^{2}=\mathbb{Z}^{2}/\ell\mathbb{Z}^{2}, =cq1/2\ell=cq^{-1/2}, both the upper and lower bounds scale like q2/log(1/q)q^{2}/\log(1/q). In view of this scaling, and the relaxation mechanism reflected in its proof, it seems that the correct scaling of the spectral gap is q2/log(1/q)q^{2}/\log(1/q) also in 2\mathbb{Z}^{2}, coinciding with the lower bound of [2]. Unfortunately, the ideas in the proof of Theorem 1.3 do not seem to be easily adapted for the model on 2\mathbb{Z}^{2}, and the problem remains open.

2. Notation

We will now recall some of the notation in [2] that will be used in this note.

  • For >0\ell>0, Λ={0,,1}d\Lambda_{\ell}=\{0,\dots,\ell-1\}^{d} and d()=d\mathbb{Z}^{d}(\ell)=\ell\mathbb{Z}^{d}.

  • For >0\ell>0 and xd()x\in\mathbb{Z}^{d}(\ell) we denote Λx=x+Λ\Lambda_{x}=x+\Lambda_{\ell} (where the \ell-dependence of Λx\Lambda_{x} is implicit).

  • The configuration space is Ω={0,1}d\Omega=\{0,1\}^{\mathbb{Z}^{d}}, and the measure μ\mu on this space is a product measure of Bernoulli random variables with parameter 1q1-q.

  • For ηΩ\eta\in\Omega and xdx\in\mathbb{Z}^{d}, the configuration which equals η\eta outside xx and different from η\eta at xx is denoted ηx\eta^{x}.

  • The FA11f critical length is denoted q=(log(1q0)log(1q))1/dCq1/d\ell_{q}=\left(\frac{\log(1-q_{0})}{\log(1-q)}\right)^{1/d}\approx Cq^{-1/d}, where q0(0,1)q_{0}\in(0,1) is given in [2, Theorem 4.1], and does not depend on qq.

  • The constraint of the FA11f dynamics, for a configuration ηΩ={0,1}d\eta\in\Omega=\{0,1\}^{\mathbb{Z}^{d}} and xdx\in\mathbb{Z}^{d}, is

    cx(η)={1y such that yx1=1 and η(y)=0,0otherwise.c_{x}(\eta)=\begin{cases}1&\exists y\text{ such that }\left\|y-x\right\|_{1}=1\text{ and }\eta(y)=0,\\ 0&\text{otherwise}.\end{cases}
  • The Dirichlet form of FA11f operating on a local function f:Ωf:\Omega\rightarrow\mathbb{R} is given by

    𝒟(f)\displaystyle\mathcal{D}(f) =xdμ(cxVarx(f))=q(1q)xdμ(cx(f(η)f(ηx))2).\displaystyle=\sum_{x\in\mathbb{Z}^{d}}\mu\left(c_{x}\operatorname{Var}_{x}(f)\right)=q(1-q)\sum_{x\in\mathbb{Z}^{d}}\mu\left(c_{x}(f(\eta)-f(\eta^{x}))^{2}\right).

For the FA11f model on a finite graph GG with vertex set V(G)V(G) denote ΩG={0,1}V(G)\Omega_{G}=\{0,1\}^{V(G)}; and μG\mu_{G} the product measure of Bernoulli random variables with parameter 1q1-q, conditioned on having at least one empty site. The constraint is defined in the same manner as in d\mathbb{Z}^{d}, except that yx1\left\|y-x\right\|_{1} should be understood as the graph distance between xx and yy, that we denote d(x,y)d(x,y). The Dirichlet form operating on f:ΩGf:\Omega_{G}\rightarrow\mathbb{R} is given by

𝒟G(f)=q(1q)xV(G)μG(cx(f(η)f(ηx))2).\mathcal{D}_{G}(f)=q(1-q)\sum_{x\in V(G)}\mu_{G}\left(c_{x}(f(\eta)-f(\eta^{x}))^{2}\right).

Throughout the proof CC will denote a generic positive constant, and qq is assumed to be small enough.

3. Proof of Theorem 1.1

In order to bound the spectral gap from above, we need to find an appropriate test function ff, such that

𝒟(f)Cq2Var(f).\mathcal{D}(f)\leq Cq^{2}\,\operatorname{Var}(f).

Consider the box Λ=Λ\Lambda=\Lambda_{\ell} for =1/q\ell=\left\lfloor 1/q\right\rfloor. For a configuration η\eta and a site xΛx\in\Lambda, the connected cluster of xx, denoted 𝒞x(η)\mathcal{C}_{x}(\eta), is defined as the set of sites yΛy\in\Lambda that are connected to xx via a path of empty sites in Λ\Lambda. If η(x)=1\eta(x)=1, its connected cluster is the empty set. This way, the set of empty sites in Λ\Lambda is partitioned in connected clusters, and we define:

f(η)=#connected clusters in Λ.f(\eta)=\#\text{connected clusters in }\Lambda. (3.1)
Proposition 3.1.

For the test function ff defined in equation (3.1),

Var(f)Cqd.\operatorname{Var}(f)\geq C\,q\,\ell^{d}. (3.2)
Proof.

This result is shown in [3] for the case of Bernoulli bond percolation. We will repeat their argument applied to our case for completeness.

First, note that we may write

f(η)=xΛ1η(x)|𝒞x(η)|,f(\eta)=\sum_{x\in\Lambda}\frac{1-\eta(x)}{|\mathcal{C}_{x}(\eta)|},

where, when η(x)=1\eta(x)=1 (and therefore 𝒞x(η)=\mathcal{C}_{x}(\eta)=\emptyset), we define 1η(x)|𝒞x(η)|=0\frac{1-\eta(x)}{|\mathcal{C}_{x}(\eta)|}=0.

Let G=3dΛG=3\mathbb{Z}^{d}\cap\Lambda, and for AGA\subseteq G, define χA(η)\chi_{A}(\eta) to be the indicator of the event, that the set {xΛ:η(y)=1y such that yx1=1}\{x\in\Lambda:\eta(y)=1\,\,\forall y\text{ such that }\left\|y-x\right\|_{1}=1\} is equal AA. Note that μ(χA)=(1q)2d|A|(1(1q)2d)|G||A|\mu(\chi_{A})=(1-q)^{2d|A|}\cdot\left(1-(1-q)^{2d}\right)^{|G|-|A|}. For such a set AA, let D(A)={yΛ:yx1=1 for some xA}D(A)=\{y\in\Lambda:\left\|y-x\right\|_{1}=1\text{ for some }x\in A\}, and define

fA(η)=xΛD(A)1η(x)|𝒞x(η)|.f_{A}(\eta)=\sum_{x\in\Lambda\setminus D(A)}\frac{1-\eta(x)}{\left|\mathcal{C}_{x}(\eta)\right|}.

When η\eta is such that χA(η)=1\chi_{A}(\eta)=1,

f(η)=fA(η)+xA(1η(x))fA(η)+nA(η).f(\eta)=f_{A}(\eta)+\sum_{x\in A}(1-\eta(x))\eqqcolon f_{A}(\eta)+n_{A}(\eta).

In order to use this identity, we split the variance over the different choices of AA:

Var(f)=AGμ((fμ(f))2χA).\operatorname{Var}(f)=\sum_{A\subseteq G}\mu\left((f-\mu(f))^{2}\chi_{A}\right).

Consider one of the summands in the above expression –

μ((fμ(f))2χA)=μ((fA(μ(f)+μ(nA))+nAμ(nA))2χA)=μ((fA(μ(f)+μ(nA))2χA)+μ((nAμ(nA))2χA)+μ((fA(μ(f)+μ(nA))(nAμ(nA))χA).\mu\left((f-\mu(f))^{2}\chi_{A}\right)=\mu\left((f_{A}-(\mu(f)+\mu(n_{A}))+n_{A}-\mu(n_{A}))^{2}\chi_{A}\right)\\ =\mu\left((f_{A}-(\mu(f)+\mu(n_{A}))^{2}\chi_{A}\right)+\mu\left((n_{A}-\mu(n_{A}))^{2}\chi_{A}\right)\\ +\mu\left((f_{A}-(\mu(f)+\mu(n_{A}))(n_{A}-\mu(n_{A}))\chi_{A}\right).

The first term is positive, and we will simply bound it by 0. In order to find the second term, we note that the variables nAn_{A} and χA\chi_{A} are independent, and therefore

μ((nAμ(nA))2χA)=μ(χA)Var(nA)=(1q)2d|A|(1(1q)2d)|G||A||A|q(1q).\mu\left((n_{A}-\mu(n_{A}))^{2}\chi_{A}\right)=\mu(\chi_{A})\operatorname{Var}(n_{A})=(1-q)^{2d|A|}\cdot\left(1-(1-q)^{2d}\right)^{|G|-|A|}\cdot\left|A\right|q(1-q).

Finally, since under the event {χA=1}\{\chi_{A}=1\} the variables fAf_{A} and nAn_{A} are independent, the third term vanishes. Therefore,

Var(f)AG(1q)2d|A|(1(1q)2d)|G||A||A|q(1q)=q(1q)2d+1|G|.\operatorname{Var}(f)\geq\sum_{A\subseteq G}(1-q)^{2d|A|}\cdot\left(1-(1-q)^{2d}\right)^{|G|-|A|}\cdot\left|A\right|q(1-q)=q(1-q)^{2d+1}\left|G\right|.

This establishes inequality (3.2). ∎

Proposition 3.2.

For the test function ff defined in equation (3.1),

𝒟(f)Cq3d.\mathcal{D}(f)\leq Cq^{3-d}. (3.3)
Proof.

Recall first that

𝒟(f)=q(1q)xdμ(cx(η)(f(ηx)f(η))2).\mathcal{D}(f)=q(1-q)\sum_{x\in\mathbb{Z}^{d}}\mu\left(c_{x}(\eta)\left(f(\eta^{x})-f(\eta)\right)^{2}\right).

Consider a single term in that sum. First, note that by flipping a single site ff could change by at most 2d2d. If xx is outside Λ\Lambda, flipping it could not change the number of clusters in Λ\Lambda and its contribution would be 0. If xx is on the boundary of Λ\Lambda (i.e., it is in Λ\Lambda and has a neighbor outside Λ\Lambda), then

μ(cx(η)(f(ηx)f(η))2)2dμ(cx(η))Cq.\mu\left(c_{x}(\eta)\left(f(\eta^{x})-f(\eta)\right)^{2}\right)\leq 2d\mu\left(c_{x}(\eta)\right)\leq C\,q.

Finally, if xx is in Λ\Lambda but has no neighbors outside Λ\Lambda, the number of open clusters could only change if it has at least two empty neighbors –

μ(cx(η)(f(ηx)f(η))2)2dμ(𝟙x has at least 2 empty neighbors)Cq2.\mu\left(c_{x}(\eta)\left(f(\eta^{x})-f(\eta)\right)^{2}\right)\leq 2d\,\mu(\mathds{1}_{x\text{ has at least }2\text{ empty neighbors}})\leq C\,q^{2}.

The proof is now concluded by summing these options –

xdμ(cx(η)(f(ηx)f(η))2)Cd1q+Cdq2=Cqd+2.\sum_{x\in\mathbb{Z}^{d}}\mu\left(c_{x}(\eta)\left(f(\eta^{x})-f(\eta)\right)^{2}\right)\leq C\ell^{d-1}q+C\ell^{d}q^{2}=Cq^{-d+2}.\qed

Theorem 1.1 follows from equations (3.2) and (3.3), together with the variational characterization of the spectral gap. ∎

4. Proof of Theorem 1.2

4.1. Upper bound

The basic tool for the proof of the upper bounds on F0(t)F_{0}(t) is the following result of [1] (see also [8, Section 4]):

Lemma 4.1.

Assume that, for some τ¯>0\overline{\tau}>0 and any local function ff which vanishes on the event {η0=0}\{\eta_{0}=0\},

μ(f2)τ¯𝒟(f).\mu(f^{2})\leq\overline{\tau}\mathcal{D}(f). (4.1)

Then F0(t)et/τ¯F_{0}(t)\leq e^{-t/\overline{\tau}}.

We will use a path argument, similar to [2, proof of Proposition 6.6], proving inequality (4.1) with the appropriate τ¯\overline{\tau}.

We start by defining a canonical geometric path, which is a discrete approximation of a straight segment. More precisely, for any zdz\in\mathbb{Z}^{d}, we will construct a nearest neighbor path γ(z)=(γ0(z),,γn(z))\gamma(z)=(\gamma_{0}(z),\dots,\gamma_{n}(z)) with γ0(z)=0\gamma_{0}(z)=0 and γn(z)=z\gamma_{n}(z)=z whose distance from the line segment [0,z]d[0,z]\in\mathbb{R}^{d} is small. The exact definition is rather cumbersome, and a reader who accepts that such a path could be constructed satisfying Observation 4.3 and Claim 4.6 (see Definition 4.5) may skip the technicalities involved in their proofs.

Definition 4.2.

Fix z=(z1,,zd)dz=(z_{1},\dots,z_{d})\in\mathbb{Z}^{d} with z1=n\left\|z\right\|_{1}=n. The canonical geometric path connecting zz to the origin is the path γ(z)=(γ0(z),,γn(z))\gamma(z)=(\gamma_{0}(z),\dots,\gamma_{n}(z)) constructed as follows – consider the set S(0,1]×{1,,d}S\subseteq(0,1]\times\{1,\dots,d\} defined as

S={(s,α):szα}.S=\left\{(s,\alpha):sz_{\alpha}\in\mathbb{Z}\right\}.

For each α\alpha there are zαz_{\alpha} values of ss for which szαsz_{\alpha}\in\mathbb{Z}, hence |S|=n|S|=n. We will order SS according to the lexicographic order, (s1,α1)<<(sn,αn)(s_{1},\alpha_{1})<\dots<(s_{n},\alpha_{n}), so that sisi+1s_{i}\leq s_{i+1} for all ii, and in case of equality αi<αi+1\alpha_{i}<\alpha_{i+1}. Then

γ0(z)\displaystyle\gamma_{0}(z) =0,\displaystyle=0,
γi(z)\displaystyle\gamma_{i}(z) =γi1(z)+eαii1.\displaystyle=\gamma_{i-1}(z)+\vec{e}_{\alpha_{i}}\quad i\geq 1.
Observation 4.3.

Fix z=(z1,,zd)dz=(z_{1},\dots,z_{d})\in\mathbb{Z}^{d} and 0iz10\leq i\leq\left\|z\right\|_{1}. Then γi(z)1=i\left\|\gamma_{i}(z)\right\|_{1}=i, i.e., the sites of the path are indexed by their norm.

Claim 4.4.

Fix z=(z1,,zd)dz=(z_{1},\dots,z_{d})\in\mathbb{Z}^{d} with z1=n\left\|z\right\|_{1}=n, and let (s1,α1)<<(sn,αn)(s_{1},\alpha_{1})<\dots<(s_{n},\alpha_{n}) be as in Definition 4.2. Then for all 1in1\leq i\leq n,

γi(z)=siz(αi),\gamma_{i}(z)=\left\lfloor s_{i}z\right\rfloor^{{\scriptscriptstyle(\alpha_{i})}},

where y(α)\left\lfloor y\right\rfloor^{{\scriptscriptstyle(\alpha)}}, for y=(y1,,yd)dy=(y_{1},\dots,y_{d})\in\mathbb{R}^{d}, is defined as

y(α)=(y1,,yα,yα+11,yd1).\left\lfloor y\right\rfloor^{{\scriptscriptstyle(\alpha)}}=(\left\lfloor y_{1}\right\rfloor,\dots,\left\lfloor y_{\alpha}\right\rfloor,\left\lceil y_{\alpha+1}\right\rceil-1,\dots\left\lceil y_{d}\right\rceil-1).
Proof.

We show this by induction. Start with i=1i=1, and consider the vector s1zs_{1}z. By the construction of s1s_{1},

0<s1zα<1\displaystyle 0<s_{1}z_{\alpha}<1 for α<α1,\displaystyle\qquad\text{for }\alpha<\alpha_{1},
s1zα=1\displaystyle s_{1}z_{\alpha}=1 for α=α1,\displaystyle\qquad\text{for }\alpha=\alpha_{1},
0<s1zα1\displaystyle 0<s_{1}z_{\alpha}\leq 1 for α>α1;\displaystyle\qquad\text{for }\alpha>\alpha_{1};

and indeed siz(α)=eα\left\lfloor s_{i}z\right\rfloor^{{\scriptscriptstyle(\alpha)}}=e_{\alpha}.

For i>1i>1, there are two options – either si=si1s_{i}=s_{i-1} and αi>αi1\alpha_{i}>\alpha_{i-1}, or si>si1s_{i}>s_{i-1}. In the first case, by induction and letting y=si1zy=s_{i-1}z,

γi(z)\displaystyle\gamma_{i}(z) =γi1(z)+eαi\displaystyle=\gamma_{i-1}(z)+e_{\alpha_{i}}
=(y1,,yαi1,yαi1+11,,yαi,,yd1).\displaystyle=(\left\lfloor y_{1}\right\rfloor,\dots,\left\lfloor y_{\alpha_{i-1}}\right\rfloor,\left\lceil y_{\alpha_{i-1}+1}\right\rceil-1,\dots,\left\lceil y_{\alpha_{i}}\right\rceil,\dots,\left\lceil y_{d}\right\rceil-1).

Since the coordinates between αi1\alpha_{i-1} and αi\alpha_{i} are not integer, we can replace 1\left\lceil\cdot\right\rceil-1 by \left\lfloor\cdot\right\rfloor, and since yαiy_{\alpha_{i}} is integer we may replace yαi\left\lceil y_{\alpha_{i}}\right\rceil by yαi\left\lfloor y_{\alpha_{i}}\right\rfloor. That is, γi(z)=y(αi)=siz(α)\gamma_{i}(z)=\left\lfloor y\right\rfloor^{{\scriptscriptstyle(\alpha_{i})}}=\left\lfloor s_{i}z\right\rfloor^{{\scriptscriptstyle(\alpha)}}.

Let us now consider the second case, where si>si1s_{i}>s_{i-1}. First, by induction, noting that the coordinates after αi1\alpha_{i-1} of si1zs_{i-1}z are not integer,

γi1(z)=(si1z1,,si1zd).\gamma_{i-1}(z)=(\left\lfloor s_{i-1}z_{1}\right\rfloor,\dots,\left\lfloor s_{i-1}z_{d}\right\rfloor).

On the other hand, we know that

si1zα<sizα<si1zα+1\displaystyle s_{i-1}z_{\alpha}<s_{i}z_{\alpha}<\left\lfloor s_{i-1}z_{\alpha}\right\rfloor+1 for α<α1,\displaystyle\qquad\text{for }\alpha<\alpha_{1},
sizα=si1zα+1\displaystyle s_{i}z_{\alpha}=\left\lfloor s_{i-1}z_{\alpha}\right\rfloor+1 for α=α1,\displaystyle\qquad\text{for }\alpha=\alpha_{1},
si1zα<sizαsi1zα+1\displaystyle s_{i-1}z_{\alpha}<s_{i}z_{\alpha}\leq\left\lfloor s_{i-1}z_{\alpha}\right\rfloor+1 for α>α1;\displaystyle\qquad\text{for }\alpha>\alpha_{1};

so siz(α)=γi1(z)+eα\left\lfloor s_{i}z\right\rfloor^{{\scriptscriptstyle(\alpha)}}=\gamma_{i-1}(z)+\vec{e}_{\alpha} and the proof is complete. ∎

Definition 4.5.

Fix >0\ell>0 and ydy\in\mathbb{Z}^{d} with y1=m\left\|y\right\|_{1}=m\leq\ell. The the \ell-cone of yy is the set

Cy()={zd:m<z1,γm(z)=y}.C_{y}^{(\ell)}=\{z\in\mathbb{Z}^{d}:m<\left\|z\right\|_{1}\leq\ell,\gamma_{m}(z)=y\}.
Claim 4.6.

Fix >0\ell>0 and yΛy\in\Lambda_{\ell} such that y1\left\|y\right\|_{1}\leq\ell. Then |Cy()|dy1d1+1|C_{y}^{(\ell)}|\leq\frac{\ell^{d}}{\left\|y\right\|_{1}^{d-1}+1}.

Proof.

First, since for y=0y=0 the cone Cy()C_{y}^{(\ell)} consists of the points in Λ\Lambda_{\ell} of norm smaller than \ell, its size is smaller than d\ell^{d}, so in what follows we may assume y0y\neq 0; and in this case we will show that the stronger inequality |Cy|dy1d1|C_{y}|\leq\frac{\ell^{d}}{\left\|y\right\|_{1}^{d-1}} holds.

Let zCyz\in C_{y}, i.e., γm(z)=y\gamma_{m}(z)=y, so by Claim 4.4 there exist ss and α\alpha such that szα=yαsz_{\alpha}=y_{\alpha} and

y=szα.y=\left\lfloor sz\right\rfloor^{\alpha}.

Assume first α=1\alpha=1, so in particular y10y_{1}\neq 0 by the construction of the geometric path. ss must be contained in {y1k+y1}k\{\frac{y_{1}}{k+y_{1}}\}_{k\in\mathbb{N}}; so we fix kk and let s=y1k+y1s=\frac{y_{1}}{k+y_{1}}, such that

z=k+y1y1(y+δ)z=\frac{k+y_{1}}{y_{1}}(y+\delta)

for some δ{0}×[0,1]d1\delta\in\{0\}\times[0,1]^{d-1}. That is, for all α>1\alpha>1,

zα(ky1+1)yα+[0,ky1+1],z_{\alpha}\in(\frac{k}{y_{1}}+1)\,y_{\alpha}+[0,\frac{k}{y_{1}}+1],

allowing at most (ky1+1)d1(\frac{k}{y_{1}}+1)^{d-1} integer choices of zz. Finally, since z1\left\|z\right\|_{1}\leq\ell, necessarily kkmax=(y11)y1k\leq k_{\text{max}}=(\frac{\ell}{\left\|y\right\|_{1}}-1)y_{1}, so, still for α=1\alpha=1, the number of possibilities for zz is bounded by

(y11)y1(kmaxy1+1)d1y1(y1)d.(\frac{\ell}{\left\|y\right\|_{1}}-1)y_{1}\cdot(\frac{k_{\text{max}}}{y_{1}}+1)^{d-1}\leq y_{1}(\frac{\ell}{\left\|y\right\|_{1}})^{d}.

Finally, summing over all possible values of α\alpha,

|Cy|dy1d(α=1dyα)=dy1d1.|C_{y}|\leq\frac{\ell^{d}}{\left\|y\right\|_{1}^{d}}\left(\sum_{\alpha=1}^{d}y_{\alpha}\right)=\frac{\ell^{d}}{\left\|y\right\|_{1}^{d-1}}.\qed

For any >0\ell>0, let

χ(η)=𝟙{xΛ,η(x)=0}.\chi_{\ell}(\eta)=\mathds{1}_{\{\exists x\in\Lambda_{\ell},\,\eta(x)=0\}}.
Claim 4.7.

Fix ηΩ\eta\in\Omega and >0\ell>0 such that χ=1\chi_{\ell}=1. Then there exists a path of configurations η(0),,η(j)\eta^{(0)},\dots,\eta^{(j)} and a sequence of sites x0,,xj1x_{0},\dots,x_{j-1} such that:

  1. (1)

    η(0)=η\eta^{(0)}=\eta and η0(j)=0\eta_{0}^{(j)}=0.

  2. (2)

    For any ii, η(i+1)=(η(i))xi\eta^{(i+1)}=\left(\eta^{(i)}\right)^{x_{i}}, and cxi(η(i))=1.c_{x_{i}}(\eta^{(i)})=1.

  3. (3)

    The sites x0,,xj1x_{0},\dots,x_{j-1} all belong to the geometric path γ(x1)\gamma(x_{1}). Moreover, each site of γ(x1)\gamma(x_{1}) appears at most twice in the sequence x0,,xj1x_{0},\dots,x_{j-1}, and in particular j2j\leq 2\ell.

  4. (4)

    For all ii, the number of sites in Λ{xi}\Lambda_{\ell}\setminus\{x_{i}\} which are empty for η(i)\eta^{(i)} is at most the number of sites in Λ{xi}\Lambda_{\ell}\setminus\{x_{i}\} which are empty for η\eta.

  5. (5)

    Fix zz, η\eta^{\prime}, xx^{\prime}. Then there exist at most one configuration η\eta and one index ii such that z=x1z=x_{1}, η=η(i)\eta^{\prime}=\eta^{(i)} and x=xix^{\prime}=x_{i}. We write (η,i)(η,x,z)(\eta,i)\sim(\eta^{\prime},x^{\prime},z).

Proof.

The path is constructed in the same manner as [2, proof of Proposition 6.6] – let zz be an empty site in Λ\Lambda_{\ell} with minimal 11-norm, and denote z1=n\left\|z\right\|_{1}=n. Then set, for i{0,,2n2}i\in\{0,\dots,2n-2\},

xi={γni21(z)i even,γni12(z)i odd.x_{i}=\begin{cases}\gamma_{n-\frac{i}{2}-1}(z)&i\text{ even},\\ \gamma_{n-\frac{i-1}{2}}(z)&i\text{ odd}.\end{cases}

This sequence defines a path η(0),,η(j)\eta^{(0)},\dots,\eta^{(j)} that indeed satisfied the conditions of the claim. ∎

Proposition 4.8.

For any local function ff that vanishes on the event {η0=0}\{\eta_{0}=0\} and for any >0\ell>0,

μ(χf2)\displaystyle\mu(\chi_{\ell}f^{2}) τ𝒟(f),\displaystyle\leq\tau_{\ell}\mathcal{D}(f),
τ\displaystyle\tau_{\ell} ={C2q1d=1,Clog2q1d=2,Cdq1d=3.\displaystyle=\begin{cases}C\ell^{2}q^{-1}&d=1,\\ C\log\ell\,\ell^{2}q^{-1}&d=2,\\ C\ell^{d}q^{-1}&d=3.\end{cases}
Proof.

First, since ff is local, we may restrict ourselves to proving the inequality for FA11f on a large finite set VdV\subset\mathbb{Z}^{d}, so the configuration space is ΩV={0,1}V\Omega_{V}=\{0,1\}^{V}. This allows us to write the Dirichlet form as

𝒟(f)\displaystyle\mathcal{D}(f) =12ηΩVxVR(ηx,η)(f(ηx)f(η))2,\displaystyle=\frac{1}{2}\sum_{\eta\in\Omega_{V}}\sum_{x\in V}R(\eta^{x},\eta)\left(f(\eta^{x})-f(\eta)\right)^{2},
R(ηx,η)\displaystyle R(\eta^{x},\eta) =R(η,ηx)=cx(η)q(1q)(μ(η)+μ(ηx)).\displaystyle=R(\eta,\eta^{x})=c_{x}(\eta)q(1-q)(\mu(\eta)+\mu(\eta^{x})).

Consider, for any η\eta, the path constructed in Claim 4.7, and for i{0,,j1}i\in\{0,\dots,j-1\} set

wi=w(xi1)=(xi1+1)(d1)/2.w_{i}=w(\left\|x_{i}\right\|_{1})=(\left\|x_{i}\right\|_{1}+1)^{(d-1)/2}.

Note that we can bound, uniformly in jj,

i=1jwi2\displaystyle\sum_{i=1}^{j}w_{i}^{-2} 2k=0w(k)W,\displaystyle\leq 2\sum_{k=0}^{\ell}w(k)\leq W,
W\displaystyle W ={2d=1,Clogd=2,Cd3;\displaystyle=\begin{cases}2\ell&d=1,\\ C\log\ell&d=2,\\ C&d\geq 3;\end{cases}

and by Claim 4.6, for every yΛy\in\Lambda_{\ell},

|Cy(d)|w(y)2\displaystyle|C_{y}^{(d\ell)}|\,w(y)^{2} Cd.\displaystyle\leq C\ell^{d}.

By the Cauchy-Schwarz inequality and the properties of the path,

μ(χf2)\displaystyle\mu\left(\chi_{\ell}f^{2}\right) =μ(χ(i=1j1wiwi(f(η(i))f(η(i1))))2)\displaystyle=\mu\left(\chi_{\ell}\left(\sum_{i=1}^{j}\frac{1}{w_{i}}\,w_{i}(f(\eta^{(i)})-f(\eta^{(i-1)}))\right)^{2}\right) (4.2)
Wiμ(wi2cxi(η(i))(f(η(i))f(η(i1)))2)\displaystyle\leq W\,\sum_{i}\mu\left(w_{i}^{2}c_{x_{i}}(\eta^{(i)})(f(\eta^{(i)})-f(\eta^{(i-1)}))^{2}\right)
Wi=12ηΩVμ(η)ηΩVxVzCx(d)𝟙η=ηi,xi=x,z=x1w(x)2cx(η)(f(ηx)f(η))2\displaystyle\leq W\sum_{i=1}^{2\ell}\sum_{\eta\in\Omega_{V}}\mu(\eta)\sum_{\eta^{\prime}\in\Omega_{V}}\sum_{x^{\prime}\in V}\sum_{z\in C_{x^{\prime}}^{(d\ell)}}\mathds{1}_{\eta^{\prime}=\eta_{i},\,x_{i}=x^{\prime},\,z=x_{1}}w(x^{\prime})^{2}c_{x^{\prime}}(\eta^{\prime})(f(\eta^{\prime x^{\prime}})-f(\eta^{\prime}))^{2}
=WηxR(ηx,η)zCx(h)iημ(η)R(ηx,η) 1(η,i)(η,x,z)w(x)2cx(η)(f(ηx)f(η))2.\displaystyle=W\sum_{\eta^{\prime}}\sum_{x^{\prime}}R(\eta^{\prime x^{\prime}},\eta^{\prime})\sum_{z\in C_{x^{\prime}}^{(h\ell)}}\sum_{i}\sum_{\eta}\frac{\mu(\eta)}{R(\eta^{\prime x^{\prime}},\eta^{\prime})}\,\mathds{1}_{(\eta,i)\sim(\eta^{\prime},x^{\prime},z)}w(x^{\prime})^{2}c_{x^{\prime}}(\eta^{\prime})(f(\eta^{\prime x^{\prime}})-f(\eta^{\prime}))^{2}.

Note that we are allowed to divide by R(ηx,η)R(\eta^{\prime x^{\prime}},\eta^{\prime}) since cx(η)=1c_{x^{\prime}}(\eta^{\prime})=1, and hence it is non-zero. We can estimate R(ηx,η)R(\eta^{\prime x^{\prime}},\eta^{\prime}) more precisely:

R(ηx,η)=q(1q)yx((1q)η(y)+q(1η(y))),R(\eta^{\prime x^{\prime}},\eta^{\prime})=q(1-q)\prod_{y\neq x^{\prime}}\left((1-q)\eta^{\prime}(y)+q(1-\eta^{\prime}(y))\right),

so

μ(η)R(ηx,η)=(1q)η(x)+q(1η(x))q(1q)yx(1q)η(y)+q(1η(y))(1q)η(y)+q(1η(y)).\frac{\mu(\eta)}{R(\eta^{\prime x^{\prime}},\eta^{\prime})}=\frac{(1-q)\eta(x^{\prime})+q(1-\eta(x^{\prime}))}{q(1-q)}\,\prod_{y\neq x^{\prime}}\frac{(1-q)\eta(y)+q(1-\eta(y))}{(1-q)\eta^{\prime}(y)+q(1-\eta^{\prime}(y))}.

By property 4 of the path we obtain

μ(η)R(ηx,η)q1.\frac{\mu(\eta)}{R(\eta^{\prime x^{\prime}},\eta^{\prime})}\leq q^{-1}.

We now conclude by continuing the estimate (4.2) –

μ(χf2)\displaystyle\mu\left(\chi_{\ell}f^{2}\right) CdWq1ηxR(ηx,η)cx(η)(f(ηx)f(η))2=Cdq1W𝒟(f).\displaystyle\leq C\ell^{d}Wq^{-1}\sum_{\eta^{\prime}}\sum_{x^{\prime}}R(\eta^{\prime x^{\prime}},\eta^{\prime})c_{x^{\prime}}(\eta^{\prime})(f(\eta^{\prime x^{\prime}})-f(\eta^{\prime}))^{2}=C\ell^{d}q^{-1}W\,\mathcal{D}(f).\qed
Remark 4.9.

If, rather than μ(f2)\mu(f^{2}) in inequality (4.1), we would like to bound μ(f)2\mu(f)^{2}, we could use Proposition 4.8 directly. By the Cauchy-Schwarz inequality

μ(f)2\displaystyle\mu(f)^{2} 2μ(χf)2+2μ((1χ)f)22μ(χf2)+2μ(1χ)μ(f2).\displaystyle\leq 2\mu\left(\chi_{\ell}f\right)^{2}+2\mu\left((1-\chi_{\ell})f\right)^{2}\leq 2\mu\left(\chi_{\ell}f^{2}\right)+2\mu(1-\chi_{\ell})\mu(f^{2}).

Choosing =C/q\ell=C/q with CC small enough, μ(1χ)\mu(1-\chi_{\ell}) is bounded below 1/41/4, so

μ(f)24μ(χf2)+Var(f)Cτ𝒟(f).\mu(f)^{2}\leq 4\mu\left(\chi_{\ell}f^{2}\right)+\operatorname{Var}(f)\leq C\tau_{\ell}\mathcal{D}(f).

This inequality is not entirely worthless, and it does bound 𝔼μ(τ0)\mathbb{E}_{\mu}(\tau_{0}) from above by CτC\tau_{\ell} (see [8, Section 4] and equation (4.3) in the following section). However, in order to obtain the exponential tail in Theorem 1.2 a more sophisticated approach is required.

From now on we set =q\ell=\ell_{q}. The lower bound on the spectral gap of [2, Theorem 6.4] is proven by introducing an auxiliary dynamics with large spectral gap and then comparing it to the FA11f dynamics. For that objective they define the constraints {c~x}xd()\{\tilde{c}_{x}\}_{x\in\mathbb{Z}^{d}(\ell)}, stating that none of the boxes Λy\Lambda_{y} is entirely occupied for yx+{e1,,ed}y\in x+\{\ell e_{1},\dots,\ell e_{d}\}; with the associated Dirichlet form

𝒟~(f)=xd()μ(c~xVarΛx(f)).\tilde{\mathcal{D}}(f)=\sum_{x\in\mathbb{Z}^{d}(\ell)}\mu\left(\tilde{c}_{x}\operatorname{Var}_{\Lambda_{x}}(f)\right).

The following Lemma is given in [2, equation (5.1)]:

Lemma 4.10.

The spectral gap associated with 𝒟~\tilde{\mathcal{D}} is at least 1/21/2.

Then, [2, Theorem 6.4] is proved by showing:

Proposition 4.11.

For any local function ff,

𝒟~(f)Cτ𝒟(f),\tilde{\mathcal{D}}(f)\leq C\tau_{\ell}\mathcal{D}(f),

where τ\tau_{\ell} is defined in Proposition 4.8.

We will use Lemma 4.10 in order to prove the following claim:

Claim 4.12.

Assume that qq is small enough, and let gg be a function vanishing on the event {χ=1}\{\chi_{\ell}=1\}. Then

μ(g2)C𝒟~(g).\mu(g^{2})\leq C\tilde{\mathcal{D}}(g).
Proof.

By Lemma 4.10, the spectral gap of the dynamics described by the Dirichlet form 𝒟~\tilde{\mathcal{D}} is at least 1/21/2. A simple application of Chebyshev’s inequality (see [8, Claim 4.11]) then yields

μ(g2)1+μ(χ)μ(χ) 2𝒟~(g),\mu(g^{2})\leq\frac{1+\mu(\chi_{\ell})}{\mu(\chi_{\ell})}\,2\mathcal{\tilde{D}}(g),

and the result follows since =q\ell=\ell_{q}, hence μ(χ)\mu(\chi_{\ell}) is bounded away from 0. ∎

We are now ready to prove the upper bound on the persistence function – consider ff which vanishes on {η0=0}\{\eta_{0}=0\}. Then, by Claim 4.12 and the fact that χ\chi_{\ell} does not depend on the occupation in Λx\Lambda_{x} for x0x\neq 0,

μ((1χ)f2)\displaystyle\mu\left((1-\chi_{\ell})f^{2}\right) C𝒟~((1χ)f)=Cxdμ(c~xVarΛx((1χ)f))\displaystyle\leq C\tilde{\mathcal{D}}\left((1-\chi_{\ell})f\right)=C\sum_{x\in\mathbb{Z}^{d}}\mu\left(\tilde{c}_{x}\operatorname{Var}_{\Lambda_{x}}\left((1-\chi_{\ell})f\right)\right)
=Cxd{0}μ(c~xVarΛx((1χ)f))+Cμ(c~0VarΛ0((1χ)f))\displaystyle=C\sum_{x\in\mathbb{Z}^{d}\setminus\{0\}}\mu\left(\tilde{c}_{x}\operatorname{Var}_{\Lambda_{x}}\left((1-\chi_{\ell})f\right)\right)+C\mu\left(\tilde{c}_{0}\operatorname{Var}_{\Lambda_{0}}\left((1-\chi_{\ell})f\right)\right)
Cxd{0}μ(c~xVarΛx(f))+CxΛ0μ(c~0Varx((1χ)f))\displaystyle\leq C\sum_{x\in\mathbb{Z}^{d}\setminus\{0\}}\mu\left(\tilde{c}_{x}\operatorname{Var}_{\Lambda_{x}}(f)\right)+C\sum_{x\in\Lambda_{0}}\mu\left(\tilde{c}_{0}\operatorname{Var}_{x}\left((1-\chi_{\ell})f\right)\right)
C𝒟~(f)+CxΛ0μ(c~0Varx((1χ)f)).\displaystyle\leq C\tilde{\mathcal{D}}(f)+C\sum_{x\in\Lambda_{0}}\mu\left(\tilde{c}_{0}\operatorname{Var}_{x}\left((1-\chi_{\ell})f\right)\right).

The first term could be bounded using Proposition 4.11. In order to bound the second term, we note that when χ(ηx)χ(η)\chi_{\ell}(\eta^{x})\neq\chi_{\ell}(\eta), necessarily χ(η1x)=0\chi_{\ell}(\eta^{1\leftarrow x})=0, where η1x\eta^{1\leftarrow x} is the configuration that equals η\eta outside xx and 11 at xx. Therefore,

Varx((1χ)f)\displaystyle\operatorname{Var}_{x}\left((1-\chi_{\ell})f\right) q(1q)((1χ(ηx))f(ηx)(1χ(η))f(η))2\displaystyle\leq q(1-q)\left((1-\chi_{\ell}(\eta^{x}))f(\eta^{x})-(1-\chi_{\ell}(\eta))f(\eta)\right)^{2}
=q(1q)𝟙χ(ηx)χ(η)(f(η1x))2\displaystyle=q(1-q)\mathds{1}_{\chi_{\ell}(\eta^{x})\neq\chi_{\ell}(\eta)}\left(f(\eta^{1\leftarrow x})\right)^{2}
q((1q)(f(η1x))2+q(f(η0x))2)\displaystyle\leq q\left((1-q)\left(f(\eta^{1\leftarrow x})\right)^{2}+q\left(f(\eta^{0\leftarrow x})\right)^{2}\right)
=qμx(f2).\displaystyle=q\mu_{x}(f^{2}).

Since c~0\tilde{c}_{0} does not depend on the occupation in Λ0\Lambda_{0},

xΛ0μ(c~0Varx((1χ)f))xΛ0qμ(c~0f2)xΛ0qμ(χ2f2)Cμ(χ2f2),\sum_{x\in\Lambda_{0}}\mu\left(\tilde{c}_{0}\operatorname{Var}_{x}\left((1-\chi_{\ell})f\right)\right)\leq\sum_{x\in\Lambda_{0}}q\mu\left(\tilde{c}_{0}f^{2}\right)\leq\sum_{x\in\Lambda_{0}}q\mu\left(\chi_{2\ell}f^{2}\right)\leq C\mu\left(\chi_{2\ell}f^{2}\right),

which could be bounded using Proposition 4.8.

We have so far shown that

μ((1χ)f2)Cτ𝒟(f)+Cτ2𝒟(f)Cτ𝒟(f).\mu\left((1-\chi_{\ell})f^{2}\right)\leq C\tau_{\ell}\mathcal{D}(f)+C\tau_{2\ell}\mathcal{D}(f)\leq C\tau_{\ell}\mathcal{D}(f).

Using Proposition 4.8 again μ(χf2)τ𝒟(f)\mu\left(\chi_{\ell}f^{2}\right)\leq\tau_{\ell}\mathcal{D}(f), and therefore

μ(f2)=μ((1χ)f2)+μ(χf2)Cτ𝒟(f),\mu(f^{2})=\mu\left((1-\chi_{\ell})f^{2}\right)+\mu\left(\chi_{\ell}f^{2}\right)\leq C\tau_{\ell}\mathcal{D}(f),

i.e., inequality (4.1) holds with τ¯=Cτ\overline{\tau}=C\,\tau_{\ell}; and Lemma 4.1 concludes the proof of the upper bound.∎

4.2. Lower bound

In order to bound the expected value of τ0\tau_{0} from below, we will use the following variational principle (see [8, Proposition 4.7]):

Lemma 4.13.

Let V0V_{0} be the space of local functions that vanish on the event {η(0)=0}\{\eta(0)=0\}. Then

𝔼μ(τ0)=supfV0(2μ(f)𝒟(f)).\mathbb{E}_{\mu}(\tau_{0})=\sup_{f\in V_{0}}(2\mu(f)-\mathcal{D}(f)).
Remark 4.14.

It would be more convenient to use a homogeneous version of that variational principle –

𝔼μ(τ0)\displaystyle\mathbb{E}_{\mu}(\tau_{0}) =supfV0supλ(2μ(λf)𝒟(λf))=supfV0supλ(2λμ(f)λ2𝒟(f)),\displaystyle=\sup_{f\in V_{0}}\sup_{\lambda\in\mathbb{R}}(2\mu(\lambda f)-\mathcal{D}(\lambda f))=\sup_{f\in V_{0}}\sup_{\lambda\in\mathbb{R}}(2\lambda\mu(f)-\lambda^{2}\mathcal{D}(f)),

and since for fixed ff the expression is maximized for λ=μ(f)𝒟(f)\lambda=\frac{\mu(f)}{\mathcal{D}(f)},

𝔼μ(τ0)=supfV0μ(f)2𝒟(f).\mathbb{E}_{\mu}(\tau_{0})=\sup_{f\in V_{0}}\frac{\mu(f)^{2}}{\mathcal{D}(f)}. (4.3)

We will now treat separately three different cases: d3d\geq 3, d=1d=1, and d=2d=2.

4.2.1. d3d\geq 3

For high dimensions, we will use the test function f(η)=η0.f(\eta)=\eta_{0}. Its expected value is 1q1-q, and its Dirichlet form is given by 𝒟(f)=q(1q)μ(c0)2q2.\mathcal{D}(f)=q(1-q)\mu(c_{0})\leq 2q^{2}. Equation (4.3) now concludes the proof of this case. ∎

4.2.2. d=1d=1

In the one dimensional case we will use a test function similar to [2, proof of Theorem 6.4]. Let =1/q\ell=\left\lceil 1/q\right\rceil,

ξ(η)=inf{|x|:ηx=0},\xi(\eta)=\inf\{|x|:\eta_{x}=0\},

and

f(η)=ξ𝟙ξ<+(2ξ)𝟙ξ<2.f(\eta)=\xi\mathds{1}_{\xi<\ell}+(2\ell-\xi)\mathds{1}_{\ell\leq\xi<2\ell}. (4.4)
Proposition 4.15.

Consider ff defined in equation (4.4). Then

μ(f)C.\mu(f)\geq C\ell.
Proof.

First, note that ξ\xi is a geometric random variable with parameter 1(1q)21-(1-q)^{2}, so we can calculate explicitly

μ(f>/2)=μ(2<ξ<32)=(1q)2/2(1(1q)2)>C,\mu(f>\ell/2)=\mu(\frac{\ell}{2}<\xi<\frac{3\ell}{2})=(1-q)^{2\,\ell/2}(1-(1-q)^{2\,\ell})>C,

and since ff is positive μ(f)C\mu(f)\geq C\ell. ∎

Proposition 4.16.

Consider ff defined in equation (4.4). Then

𝒟(f)4q.\mathcal{D}(f)\leq 4q.
Proof.

In order to bound 𝒟(f)\mathcal{D}(f) we make the following observations –

  1. (1)

    For fixed η\eta, if f(η)f(ηx)f(\eta)\neq f(\eta^{x}) then either ξ(η)=|x|\xi(\eta)=|x|, or ξ(η)=|x|+1\xi(\eta)=|x|+1.

  2. (2)

    For fixed η\eta, if f(η)f(ηx)f(\eta)\neq f(\eta^{x}) then (f(η)f(ηx))2=1\left(f(\eta)-f(\eta^{x})\right)^{2}=1.

With these observations in mind,

𝒟(f)\displaystyle\mathcal{D}(f) =q(1q)xμ(cx(f(ηx)f(η))2)\displaystyle=q(1-q)\sum_{x}\mu\left(c_{x}(f(\eta^{x})-f(\eta))^{2}\right)
q(1q)xμ(cx(𝟙ξ=|x|+𝟙ξ=|x|+1))\displaystyle\leq q(1-q)\sum_{x}\mu\left(c_{x}(\mathds{1}_{\xi=|x|}+\mathds{1}_{\xi=|x|+1})\right)
4q(1q)\displaystyle\leq 4q(1-q)\qed

Using these two propositions and equation (4.3), the case d=1d=1 is concluded.∎

4.2.3. d=2d=2

Let =q\ell=\ell_{q}, recalling q=Cq1/2\ell_{q}=Cq^{-1/2}, and Λ={x2:x1}\Lambda=\{x\in\mathbb{Z}^{2}:\left\|x\right\|_{1}\leq\ell\}. The test function we will use is

f(η)\displaystyle f(\eta) =infx2η(x)=0log(1+x1).\displaystyle=\inf_{\begin{subarray}{c}x\in\mathbb{Z}^{2}\\ \eta(x)=0\end{subarray}}\log(1+\left\|x\right\|_{1}\wedge\ell). (4.5)

Note that it vanishes on the event {η(0)=0}\{\eta(0)=0\}, and that it depends only on the occupation in Λ\Lambda.

Remark 4.17.

The function log(1+x1)\log(1+\left\|x\right\|_{1}) is used in [7] in a different context, in order to bound the relaxation time of the simple random walk on a certain graph that consists of two copies of Λn\Lambda_{n} (for some nn\in\mathbb{N}). Though presented differently, the proof there is based on the fact that this function serves as a test function for the hitting time at 0 of the random walk on Λn\Lambda_{n}; and that for the dynamics to relax a random walk in one of the two copies of Λn\Lambda_{n} must first hit 0. Indeed, the bound Cq2log(1/q)Cq^{-2}\log(1/q) obtained scales as the expected hitting time at the origin for a random walk in Λ\Lambda with jump rate qq.

Proposition 4.18.

Consider ff defined in equation (4.5). Then

μ(f)Clog(1/q).\mu(f)\geq C\log(1/q).
Proof.

The proof is based on the fact, that the probability that Λ\Lambda is entirely occupied, given by (1q)|Λ|(1-q)^{\left|\Lambda\right|}, is bounded away from 0 uniformly in qq (thanks to the choice =Cq1/2\ell=Cq^{-1/2}). In this case, ff equals log(1+)\log(1+\ell), which is greater than 12log(1/q)\frac{1}{2}\log(1/q). ∎

Proposition 4.19.

Consider ff defined in equation (4.5). Then

𝒟(f)Cq2log(1/q).\mathcal{D}(f)\leq Cq^{2}\log(1/q).
Proof.

The proof is based on the following observation:

Observation 4.20.

Fix ηΩ\eta\in\Omega and xΛx\in\Lambda such that cx(η)=1c_{x}(\eta)=1, η(x)=0\eta(x)=0, and f(η)f(ηx)f(\eta)\neq f(\eta^{x}). Then f(η)=log(1+x1)f(\eta)=\log(1+\left\|x\right\|_{1}) and f(ηx)=log(2+x1)f(\eta^{x})=\log(2+\left\|x\right\|_{1}).

Proof.

Since f(η)f(ηx)f(\eta)\neq f(\eta^{x}), there can be no vertex yy with η(y)=0\eta(y)=0 and y1x1\left\|y\right\|_{1}\leq\left\|x\right\|_{1}, so in particular f(η)=log(1+x1)f(\eta)=\log(1+\left\|x\right\|_{1}). Moreover, since cx=1c_{x}=1, it must have an empty neighbor zz, and since this neighbor has norm greater than x1\left\|x\right\|_{1}, necessarily z1=x1+1\left\|z\right\|_{1}=\left\|x\right\|_{1}+1. Since, in addition, no empty site for ηx\eta^{x} has norm strictly smaller than x1+1\left\|x\right\|_{1}+1, we conclude that f(ηx)=log(2+x1)f(\eta^{x})=\log(2+\left\|x\right\|_{1}). ∎

This observation implies in particular that for all ηΩ\eta\in\Omega and xΛx\in\Lambda such that cx(η)=1c_{x}(\eta)=1

(f(ηx)f(η))2(1+x1)2.\left(f(\eta^{x})-f(\eta)\right)^{2}\leq\left(1+\left\|x\right\|_{1}\right)^{-2}.

Using this estimate,

𝒟(f)\displaystyle\mathcal{D}(f) =q(1q)xμ(cx(f(ηx)f(η))2)\displaystyle=q(1-q)\sum_{x}\mu\left(c_{x}(f(\eta^{x})-f(\eta))^{2}\right)
q(1q)xΛμ(cx)(1+x1)2\displaystyle\leq q(1-q)\sum_{x\in\Lambda}\mu\left(c_{x}\right)\left(1+\left\|x\right\|_{1}\right)^{-2}
Cq2log()=Cq2log(1/q)\displaystyle\leq Cq^{2}\log(\ell)=Cq^{2}\log(1/q)\qed

The proof of Theorem 1.2 is then concluded by the last two propositions and equation (4.3).∎

5. Proof of Theorem 1.3

As in the proof of Theorem 1.1, we look for f:ΩGf:\Omega_{G}\rightarrow\mathbb{R} such that

𝒟G(f)Cdmaxqτ¯meet1Var(f),\mathcal{D}_{G}(f)\leq Cd_{\max}q\overline{\tau}_{\operatorname{meet}}^{-1}\,\operatorname{Var}(f),

where the variance is understood with respect to the measure μG\mu_{G}.

The test function that we use is

f(η)=maxx,yV(G)η(x)=η(y)=0τmeet(x,y).f(\eta)=\max_{\begin{subarray}{c}x,y\in V(G)\\ \eta(x)=\eta(y)=0\end{subarray}}\tau_{\text{meet}}(x,y). (5.1)

Before analyzing the function ff, note that τmeet\tau_{\text{meet}} solves the following Poisson problem:

RW(τmeet(x,y))\displaystyle-\mathcal{L}_{\text{RW}}\left(\tau_{\text{meet}}(x,y)\right) =1,d(x,y)>1,\displaystyle=1,\qquad d(x,y)>1,
τmeet(x,y)\displaystyle\tau_{\text{meet}}(x,y) =0,d(x,y)1;\displaystyle=0,\qquad d(x,y)\leq 1;

where RW\mathcal{L}_{\text{RW}} is the infinitesimal generator of two independent random walks on GG. Multiplying both sides by τmeet(x,y)\tau_{\text{meet}}(x,y) and averaging over xx and yy we obtain

τ¯meet=𝒟RW(τmeet),\overline{\tau}_{\text{meet}}=\mathcal{D}_{\text{RW}}\left(\tau_{\text{meet}}\right), (5.2)

where the Dirichlet form is given for every g:V(G)×V(G)g:V(G)\times V(G)\rightarrow\mathbb{R} by

𝒟RW(g)=12|V(G)|2xy[xx(g(x,y)g(x,y))2+yy(g(x,y)g(x,y))2].\mathcal{D}_{\text{RW}}\left(g\right)=\frac{1}{2\left|V(G)\right|^{2}}\sum_{x}\sum_{y}\left[\sum_{x^{\prime}\sim x}\left(g(x^{\prime},y)-g(x,y)\right)^{2}+\sum_{y^{\prime}\sim y}\left(g(x,y^{\prime})-g(x,y)\right)^{2}\right].
Proposition 5.1.

For ff defined in equation (5.1),

Var(f)Cτ¯meet2.\operatorname{Var}(f)\geq C\overline{\tau}_{\operatorname{meet}}^{2}.
Proof.

Under μG,\mu_{G}, the probability that exactly one site is empty is of order 11 (i.e., bounded away from 0 uniformly in qq). When this happens, f=0f=0, so

Var(f)\displaystyle\operatorname{Var}(f) =μG[(fμG(f))2]μG[(fμG(f))2𝟙x(1η(x))=1]\displaystyle=\mu_{G}\left[\left(f-\mu_{G}(f)\right)^{2}\right]\geq\mu_{G}\left[\left(f-\mu_{G}(f)\right)^{2}\mathds{1}_{\sum_{x}(1-\eta(x))=1}\right]
=μG(f)2μG[x(1η(x))=1]CμG(f)2.\displaystyle=\mu_{G}(f)^{2}\,\mu_{G}\left[\sum_{x}(1-\eta(x))=1\right]\geq C\mu_{G}(f)^{2}.

Moreover, the probability that there are exactly two vacancies is also of order 11. Under this event,

μG[f(η)|𝟙x(1η(x))=2]=1|V(G)|(|V(G)|1)x,yV(G)xyτmeet(x,y)τ¯meet.\mu_{G}\left[f(\eta)\middle|\mathds{1}_{\sum_{x}(1-\eta(x))=2}\right]=\frac{1}{\left|V(G)\right|(\left|V(G)\right|-1)}\sum_{\begin{subarray}{c}x,y\in V(G)\\ x\neq y\end{subarray}}\tau_{\text{meet}}(x,y)\geq\overline{\tau}_{\text{meet}}.\qed
Proposition 5.2.

For ff defined in equation (5.1),

𝒟(f)Cqτ¯meet.\mathcal{D}(f)\leq Cq\overline{\tau}_{\operatorname{meet}}.
Proof.

We start the proof with an observation:

Observation 5.3.

Let ηΩG\eta\in\Omega_{G} and xV(G)x\in V(G) such that cx(η)=1c_{x}(\eta)=1, η(x)=0\eta(x)=0, and f(η)f(ηx)f(\eta)\neq f(\eta^{x}). Then there exist x,yV(G)x^{\prime},y\in V(G) such that xxx^{\prime}\sim x, d(x,y)>1d(x,y)>1, η(y)=η(x)=0\eta(y)=\eta(x^{\prime})=0, and

τmeet(x,y)f(ηx)<f(η)=τmeet(x,y).\tau_{\operatorname{meet}}(x^{\prime},y)\leq f(\eta^{x})<f(\eta)=\tau_{\operatorname{meet}}(x,y).
Proof.

First, recalling equation (5.1), when filling an empty site ff could only decrease, and since f(η)f(ηx)f(\eta)\neq f(\eta^{x}) necessarily f(η)>f(ηx)f(\eta)>f(\eta^{x}). Moreover, ff could only change if the maximum is attained at the pair x,yx,y for some yV(G)y\in V(G), i.e., f(η)=τmeet(x,y)f(\eta)=\tau_{\text{meet}}(x,y). Note that f(η)f(\eta) is non-zero, hence d(x,y)>1d(x,y)>1. Finally, cx(η)=1c_{x}(\eta)=1 means that xx has an empty neighbor xx^{\prime}; and since in the configuration ηx\eta^{x} both xx^{\prime} and yy are empty f(ηx)τmeet(x,y)f(\eta^{x})\geq\tau_{\text{meet}}(x^{\prime},y). ∎

As a consequence of this observation, for all ηΩG\eta\in\Omega_{G} and xx such that cx(η)=1c_{x}(\eta)=1,

(f(ηx)f(η))2yV(G)d(x,y)>1xV(G)xx(1η(y))(1η(x))(τmeet(x,y)τmeet(x,y))2.\left(f(\eta^{x})-f(\eta)\right)^{2}\leq\sum_{\begin{subarray}{c}y\in V(G)\\ d(x,y)>1\end{subarray}}\sum_{\begin{subarray}{c}x^{\prime}\in V(G)\\ x^{\prime}\sim x\end{subarray}}(1-\eta(y))(1-\eta(x^{\prime}))\,\left(\tau_{\text{meet}}(x,y)-\tau_{\text{meet}}(x^{\prime},y)\right)^{2}.

We can now use this estimate and calculate the Dirichlet form:

𝒟(f)\displaystyle\mathcal{D}(f) =q(1q)xμG[cx(f(ηx)f(η))2]\displaystyle=q(1-q)\sum_{x}\mu_{G}\left[c_{x}\left(f(\eta^{x})-f(\eta)\right)^{2}\right]
q(1q)xμG[yV(G)d(x,y)>1xx(1η(y))(1η(x))(τmeet(x,y)τmeet(x,y))2]\displaystyle\leq q(1-q)\sum_{x}\mu_{G}\left[\sum_{\begin{subarray}{c}y\in V(G)\\ d(x,y)>1\end{subarray}}\sum_{x^{\prime}\sim x}(1-\eta(y))(1-\eta(x^{\prime}))\,\left(\tau_{\text{meet}}(x,y)-\tau_{\text{meet}}(x^{\prime},y)\right)^{2}\right]
q3(1q)xyV(G)d(x,y)>1xx(τmeet(x,y)τmeet(x,y))2\displaystyle\leq q^{3}(1-q)\sum_{x}\sum_{\begin{subarray}{c}y\in V(G)\\ d(x,y)>1\end{subarray}}\sum_{x^{\prime}\sim x}\left(\tau_{\text{meet}}(x,y)-\tau_{\text{meet}}(x^{\prime},y)\right)^{2}
Cq𝒟RW(τmeet),\displaystyle\leq Cq\,\mathcal{D}_{\text{RW}}(\tau_{\text{meet}}),

and the proposition follows form equation (5.2). ∎

Theorem 1.3 is a consequence of Propositions 5.1 and 5.2, using the variational characterization of the spectral gap. ∎

Remark 5.4.

For the case of the two dimensional torus discussed in Remark 1.4, it is possible to bound τ¯meet\overline{\tau}_{\text{meet}} from below using the equivalent of equation (4.3):

τ¯meet=supg(1|V(G)|2x,yg(x,y))2𝒟RW(g),\overline{\tau}_{\operatorname{meet}}=\sup_{g}\frac{\left(\frac{1}{\left|V(G)\right|^{2}}\sum_{x,y}g(x,y)\right)^{2}}{\mathcal{D}_{\text{RW}}(g)},

where the supremum is taken over functions g:V(G)×V(G)g:V(G)\times V(G)\rightarrow\mathbb{R}, for which g(x,y)=0g(x,y)=0 whenever d(x,y)1d(x,y)\leq 1. Taking the test function g(x,y)=log(d(x,y)1)g(x,y)=\log(d(x,y)\lor 1) (cf. Remark 4.17) yields the bound q2/log(1/q)q^{2}/\log(1/q) on the spectral gap.

Acknowledgments

I wish to thank Ivailo Hartarsky and Fabio Martinelli for the discussions and comments.

References

  • [1] Amine Asselah and Paolo Dai Pra, Quasi-stationary measures for conservative dynamics in the infinite lattice, Ann. Probab. 29 (2001), no. 4, 1733–1754. MR 1880240
  • [2] Nicoletta Cancrini, Fabio Martinelli, Cyril Roberto, and Cristina Toninelli, Kinetically constrained spin models, Probab. Theory Related Fields 140 (2008), no. 3-4, 459–504. MR 2365481
  • [3] J. Theodore Cox and Geoffrey Grimmett, Central limit theorems for associated random variables and the percolation model, Ann. Probab. 12 (1984), no. 2, 514–528. MR 735851
  • [4] Ivailo Hartarsky, Fabio Martinelli, and Cristina Toninelli, Coalescing and branching simple symmetric exclusion process, arXiv preprint arXiv:2006.01426 (2020).
  • [5] Natesh S. Pillai and Aaron Smith, Mixing times for a constrained Ising process on the torus at low density, Ann. Probab. 45 (2017), no. 2, 1003–1070.
  • [6] by same author, Mixing times for a constrained Ising process on the two-dimensional torus at low density, Ann. Inst. H. Poincaré Probab. Statist. 55 (2019), no. 3, 1649–1678.
  • [7] Laurent Saloff-Coste, Lectures on finite markov chains, Lecture Notes in Mathematics, vol. 1665, pp. 301–413, Springer, 1996.
  • [8] Assaf Shapira, Kinetically constrained models with random constraints, Ann. Appl. Probab. 30 (2020), no. 2, 987–1006.