January 01, 1970
Consider the following plaquette model from statistical physics: a lamp lies at every vertex of the triangular lattice and a switch lies at every even vertex of the (bipartite) dual hexagonal lattice. Each switch toggles the three lamps on its face. The energy of a configuration is the number of ON lamps.
For the Glauber dynamics associated with the Gibbs measure defined by this Hamiltonian at any inverse temperature \(\beta>0\), we show that, in any dimension \(d\geqslant 2\), the infinite volume relaxation time satisfies \[e^{\beta^2/C}/C \leqslant T_{\mathrm{rel}}\leqslant Ce^{e^{C\beta}}\] for some \(C>0\). Our result entails that the Gibbs measure is unique. The \(e^{\beta^2}\) scaling was conjectured by Newman and Moore in 1999 and matches the behaviour of supercritical rooted kinetically constrained models such as the East model, thus recovering fragile glass phenomenology in the absence of kinetic constraints. More precisely, we show that, on a torus of side length \(2^k\), when \(\beta\to\infty\) and \(k/\beta\to0\), we have \(T_{\mathrm{rel}}=e^{2\beta k(1+o(1))}\). Quite surprisingly, however, we also prove that, on non-periodic finite domains of size \(n\leqslant e^{\beta/C}\) for large \(C>0\), we have the much larger asymptotics \(\ln T_{\mathrm{rel}}=\beta n^{\Theta(1)}\).
The main ingredients of the proofs are new results in extremal and enumerative combinatorics and rely on renormalisation ideas for the dynamics and its groundstates also known as the Ledrappier subshift. We note consequences of our results to geometric group theory (more precisely to the complexity of the word problem for the Baumslag finitely presented group) and to ergodic theory.
MSC2020: 60K35; 05A16; 05D99; 20F65; 37A25; 37B10; 60C05; 82C20
Keywords: triangular plaquette model; Newman-Moore model; relaxation time; Ledrappier subshift; Isodiametral function; Baumslag group
One of the great achievements of statistical physics is the development of discrete, finite models whose behaviour closely follows, and ideally explains, physical phenomena. Foremost among these, the Ising model Lenz20?, Ising25? in dimension \(d\geqslant 2\) exhibits phase transitions as manifested in magnetic bodies: there is a \(\pm1\) spin at every lattice vertex, and the system’s energy is computed from spin values and nearest-neighbour interactions. In a dynamical perspective, the spins evolve independently at random so as to minimise the energy, a temperature parameter dictating the likelihood of frustrated (\({+}{-}\)) bonds.
Several paradigms have been proposed to explain the more complex behaviour of glassy materials in condensed matter physics (see e.g. Arceri21? for an overview of this vast domain). One hallmark of such materials is a super-Arrhenius divergence of relaxation time scales at high inverse temperature \(\beta>0\). These time scales may be viewed as quantifying the exponential decay rate of spin correlations with time or the viscosity of the system. For instance, a quadratic scaling \(T_{\mathrm{rel}}\approx e^{\beta^2}\) can supply a good quantitative fit for the behaviour experimentally measured in so-called fragile glass materials. In any case, the emergence of energy barriers whose size diverges when approaching the possibly degenerate phase transition point is to be expected.
One such paradigm, introduced by Newman and Moore Newman99? in 1999 (a similar model goes back to Baxter and Wu Baxter73? from 1973) can be viewed as a natural generalisation of the one-dimensional Ising model to \(d\) dimensions with \((d+1)\)-body interactions. We concentrate on their triangular plaquette model, and present it informally (see Section 1.1 for a more formal definition). Consider \(\pm1\) spins arranged in a triangular lattice \(\Lambda\), in which the energy is determined by the number of frustrated upwards-pointing triangles; more precisely, the energy of the system with spins \(\sigma\in\{\pm1\}^\Lambda\) is written as a Hamiltonian \[H(\sigma)=-\sum_{\Delta\subset\Lambda}\prod_{x\in\Delta}\sigma_x\] (the sum over upwards pointing triangles \(\Delta\) is infinite, but only relative energies will be considered). The probability of a spin \(\sigma\) is given by the Gibbs measure \(\mu(\sigma)\propto\exp(-\beta H(\sigma))\). Spins may be made to evolve under Glauber dynamics: each site of \(\Lambda\) is equipped with a standard Poisson clock; when the clock at site \(x\) rings, \(\sigma\) is updated to \(\sigma'\) with \(\sigma'_x=-\sigma_x\) with probability \(1/(1+\exp(\beta(H(\sigma')-H(\sigma))))\). This leads to numerical estimates on the relaxation time, and experimental support for their glassy (super-Arrhenius) behaviour, for which we refer to Jack05?, Garrahan02?, Jack05a?, Newman99?, and Inack22? for a more recent study.
An equivalent description of the model, also going back to Newman99?, is as follows, see Figure 1. The spins are switches (\(\pm1\)) at all vertices of the triangular lattice \(\Lambda\); there are lamps at another triangular lattice, consisting of even positions of the dual lattice to \(\Lambda\), which are ON when the product of the neighbouring switches is \(-1\) and OFF otherwise; so toggling a switch toggles simultaneously three lamps. Each ON lamp consumes an amount \(\exp(2\beta)\) of energy. There are numerous popular games in which lamps are controlled by complicated configurations of switches, see e.g. Marcarini26? for a recent account. The triangular plaquette model studied in this paper corresponds to a directed (North-East) version of the ‘Lights Out’ game or its \(3\times 3\) predecessor, the ‘Magic square’ electronic game, from 1978. In general, the interest and difficulty of these games is related to the difficulty of finding canonical paths between given configurations.
Informally, our results are as follows (see Section 1.2 for formal statements). Firstly, if we work on a torus of side length \(2^k\) with \(k/\beta\to0\) and \(\beta\to\infty\), then \(\ln T_{\mathrm{rel}}\sim 2 k\beta\), proving the prediction of Newman99?, where precisely this setting was considered (namely, \(k\to\infty\) after \(\beta\to\infty\)). Secondly, if we work in infinite volume, for any dimension \(d\geqslant 2\), suitably large \(C>0\) and any \(\beta>0\) we have \[e^{\beta^2/C}/C\leqslant T_{\mathrm{rel}}\leqslant Ce^{e^{C\beta}}.\] While these bounds may seem distant, we show that, strikingly, in finite volume of size \(n\) without periodic boundary condition, the upper bound is close to the truth for \(\ln n/\beta\) small, while we expect the lower bound to be essentially sharp when \(\ln n/\beta\) is large. The latter would imply \(T_{\mathrm{rel}}=e^{\Theta(\beta^2)}\) in infinite volume, despite the much slower relaxation in finite volume.
Let \(d\) be a positive integer. We will work on \({\mathbb{Z}} ^d\) rather than the triangular lattice (and higher dimensional generalisations) for convenience of notation, but this is equivalent up to a linear transformation of the lattice. We denote the canonical basis of \({\mathbb{R}} ^d\) by \(e_1,\dots,e_d\). We introduce the (simplex) plaquette \[\label{eq:def:Td} T_d=\{0,e_1,e_2,\dots,e_d\}.\tag{1}\] For any finite set \(\Lambda\), let \(\Omega_\Lambda=\{1,-1\}^\Lambda\) and \(\Omega=\Omega_{{\mathbb{Z}} ^d}\). We denote by \(\mathbf{1}_\Lambda\in\Omega_\Lambda\) the constant configuration and omit the subscript \(\Lambda\) when it is clear from the context. For \(\sigma\in\Omega_\Lambda\) and \(x\in\Lambda\), we denote by \(\sigma_x\) the value of \(\sigma\) at site \(x\). For \(V\subseteq\Lambda\) and \(\sigma\in\Omega_\Lambda\), we write \(\sigma_V\in\Omega_V\) for the restriction of \(\sigma\) to \(V\). For \(\sigma\in\Omega_\Lambda\) and finite \(V\subseteq\Lambda\), we write \([\sigma]_V=\prod_{x\in V}\sigma_x\). For disjoint \(\Lambda,\Lambda'\subseteq {\mathbb{Z}} ^d\) and \(\sigma\in\Omega_\Lambda,\sigma'\in\Omega_{\Lambda'}\), we write \(\sigma\cdot\sigma'\in\Omega_{\Lambda\cup\Lambda'}\) for the configuration equal to \(\sigma\) on \(\Lambda\) and to \(\sigma'\) on \(\Lambda'\). For \(x\in\Lambda\) and \(\sigma\in\Omega_\Lambda\), we denote by \(\sigma^x\in\Omega_\Lambda\) the configuration such that \(\sigma^x_x=-\sigma_x\) and \(\sigma_{\Lambda\setminus\{x\}}=\sigma_{\Lambda\setminus\{x\}}^x\), that is, the result of flipping the state of \(x\). A function \(f\colon\Omega\to{\mathbb{R}}\) is local, if there exists a finite \(V\subset{\mathbb{Z}} ^d\) such that for any \(\sigma,\sigma'\in\Omega\) such that \(\sigma_V=\sigma'_V\), we have \(f(\sigma)=f(\sigma')\). We call the smallest set \(V\) as above the support of \(f\) and denote it by \(\mathop{\mathrm{supp}}f\). We denote the integral of a function \(f\) with respect to a measure \(\mu\) by \(\mu(f)\).
For finite \(\Lambda\) and \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus \Lambda}\) (viewed as a boundary condition) and \(\sigma\in\Omega_\Lambda\), we define the Hamiltonian, partition function and Boltzmann measure by \[\begin{align} \label{eq:def:H} H_\Lambda^\tau(\sigma)&{}=-\sum_{x\in\Lambda-T_d}[\sigma\cdot\tau]_{x+T_d},&Z_\Lambda^\tau&{}=\sum_{\sigma\in\Omega_\Lambda}e^{-\beta H_\Lambda^\tau(\sigma)},&\mu_{\Lambda}^{\tau}(\sigma)&{}=\frac{e^{-\beta H_\Lambda^\tau(\sigma)}}{Z_\Lambda^\tau}, \end{align}\tag{2}\] keeping the inverse temperature parameter \(\beta>0\) implicit. When \(\Lambda=\{x\}\) is a singleton, we write \(x\) instead of \(\Lambda\) in the above notation. The degenerate cases \(\beta=0\) and \(\beta=\infty\) (which we do not consider unless otherwise stated) correspond to \(\mu_\Lambda^\tau\) being the uniform measure on \(\Omega_\Lambda\) and on the groundstates \(\mathop{\mathrm{argmin}}H_\Lambda^\tau\), respectively. Notice that, for any local \(f\) and finite \(\Lambda\), \(\mu_\Lambda(f)\) is a function from \(\Omega_{{\mathbb{Z}} ^d\setminus\Lambda}\) to \({\mathbb{R}}\), which we identify with a function from \(\Omega\) to \({\mathbb{R}}\) by projection. We denote the variance associated to \(\mu_\Lambda^\tau\) by \(\mathop{\mathrm{Var}}_\Lambda^\tau:f\mapsto\mu_\Lambda^\tau(f^2)-(\mu_\Lambda^\tau(f))^2\).
We denote by \(\mu\) any infinite volume Gibbs measure, that is, any measure on \(\Omega\) such that, for all finite \(\Lambda\) and local function \(f\) with \(\mathop{\mathrm{supp}}f\subseteq\Lambda\), we have \[\label{eq:Gibbs:DLR} \int \mu_\Lambda^{\sigma_{{\mathbb{Z}} ^d\setminus\Lambda}}(f)\mathrm d\mu(\sigma)=\mu(f).\tag{3}\] It follows from our results that \(\mu\) is unique for any \(\beta>0\), but this is not clear a priori. In order to lighten notation, we set \[\mu_\Lambda(f)\mathrel{\vcenter{:}}=\mu_\Lambda^{\sigma_{{\mathbb{Z}} ^d\setminus\Lambda}}(f),\] whenever the configuration \(\sigma\) is clear from the context, e.g.@eq:eq:Gibbs:DLR reads \(\mu(\mu_\Lambda(f))=\mu(f)\).
The Glauber dynamics in a domain \(\Lambda\subseteq{\mathbb{Z}} ^d\) with boundary condition \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda}\) is the (continuous time) Markov process whose generator acts on local functions \(f\colon\Omega_\Lambda\to{\mathbb{R}}\) via \[\mathcal{L}_\Lambda^\tau f(\sigma)=\sum_{x\in\Lambda}\left(\mu_x^{\tau\cdot\sigma_{\Lambda\setminus\{x\}}}(f)-f\right)(\sigma).\] Equivalently, one may view this process via the following more intuitive graphical construction. Each site of \(\Lambda\) is equipped with a Poisson clock, which rings at independent intervals of time with exponential distribution of mean 1. When the clock at site \(x\) rings, the current configuration \(\sigma\) is updated to either \(\sigma\) or \(\sigma^x\) with probability such that \(\sigma_x\) has distribution \(\mu_x^{\tau\cdot\sigma_{\Lambda\setminus \{x\}}}\) after the update. It is classical that such finite range dynamics are well-defined in infinite volume (see e.g.Liggett05?) and the generator \(\mathcal{L}=\mathcal{L}_{{\mathbb{Z}} ^d}\) extends to a self-adjoint operator in \(L^2(\Omega,\mu)\).
The Dirichlet form associated to the generator is given by \[\label{eq:def:D} \mathcal{D}_\Lambda^\tau(f)=-\mu_\Lambda^\tau(f\mathcal{L}_\Lambda^\tau f)=\sum_{x\in\Lambda}\mu_\Lambda^\tau(\mathop{\mathrm{Var}}_x(f)),\tag{4}\] where, as before, \(\mathop{\mathrm{Var}}_x\) denotes the variance with respect to \(\mu_{\{x\}}^{\tau\cdot\sigma_{{\mathbb{Z}} ^d\setminus\Lambda}}\) with \(\sigma\) the configuration over which \(\mu_\Lambda^\tau\) runs. We may now define the relaxation time as \[\label{eq:def:trel} T_{\mathrm{rel}}^{\Lambda,\tau}=\left(\inf\left\{\frac{\mathcal{D}_\Lambda^\tau(f)}{\mathop{\mathrm{Var}}_\Lambda^\tau(f)}:f\text{ local},\mathop{\mathrm{supp}}f\subseteq\Lambda,\mathop{\mathrm{Var}}_\Lambda^\tau(f)\neq 0\right\}\right)^{-1}\in[0,\infty]\tag{5}\] and \(T_{\mathrm{rel}}=T_{\mathrm{rel}}^{{\mathbb{Z}} ^d}\). This fundamental quantity has many equivalent definitions, in particular its inverse is known as the spectral gap of \(\mathcal{L}_\Lambda^\tau\). It governs the exponential rate at which correlations decay in the system. We direct the reader to Martinelli99? for background on Glauber dynamics though we will not assume familiarity with this classical reference.
All the above quantities are also apply with periodic boundary condition as follows. For positive integer \(n\), let \[{\mathbb{T}} _n=({\mathbb{Z}} /n{\mathbb{Z}} )^d\] denote the torus of dimension \(d\) and (linear) size \(n\). In this case, no boundary condition needs to be specified, so we write \(\mu_{{\mathbb{T}} _n}\), \(\mathcal{L}_{{\mathbb{T}} _n}\) \(\mathcal{D}_{{\mathbb{T}} _n}\), \(T_{\mathrm{rel}}^{{\mathbb{T}} _n}\) for the corresponding Boltzmann measure, generator, Dirichlet form and relaxation time defined as above.
We denote by \(\Lambda_{d,n}\) the simplex of side length \(n\): \[\label{eq:def:Lambda:dn} \Lambda_{d,n}\mathrel{\vcenter{:}}=\left\{(x_1,\dots,x_d)\in \mathbb{N}^d \mid x_1+\dots+x_d\leqslant n\right\}\subset \mathbb{Z}^d.\tag{6}\] (throughout this text \(\mathbb{N}\) denotes the natural numbers and contains \(0\)) and we call simplex a translate \(x+\Lambda_{d,n}\) for some \(x\in{\mathbb{Z}} ^d\).
We shall consider at various moments switch or spin configurations, that take values in \(\{\pm1\}\), and lamp configurations, that take values in \(\{0,1\}\) or are represented as subsets via their characteristic function. Even though both are considered on domains \(\Lambda_{d,n}\), there should be no confusion.
We are now ready to state our main results. We begin with those pertaining to infinite-volume configurations.
Theorem 1 (Infinite volume relaxation time). For any \(d\geqslant 2\) there exists \(C=C(d)\geqslant 1\) such that, for all inverse temperature \(\beta>0\), we have \[e^{\beta^2/C}/C\leqslant T_{\mathrm{rel}}\leqslant Ce^{e^{C\beta}}\] and the Gibbs state \(\mu\) is unique.
From the physics perspective, the lower bound establishing slowness of the dynamics is the most important and matches the prediction of Newman and Moore Newman99?. The fact that \(T_{\mathrm{rel}}<\infty\) when \(d=2\) was established in Chleboun17?, while a quantitative upper bound of order \(e^{e^{C\beta}}\) follows immediately from that work by Martinelli99?*Theorem 3.8 and Proposition 4.10, since strong spatial mixing is known to imply finite relaxation time Martinelli94a? (see Ott25? for more recent work on mixing conditions in two dimensions). The approach of Chleboun17? does not generalise to higher dimension. The uniqueness of Gibbs state implies that the model does not undergo a phase transition.
Turning to finite volume, we start with (our formal interpretation of) the conjecture of Newman99?*Figure 3 and (13) concerning tori of size \(2^k\) with subcritical size.
Theorem 2 (Subcritical asymptotics for dyadic tori). In dimension \(d=2\), we have \[\lim_{\substack{\beta\to \infty\\k/\beta\to0}}\frac{\ln T_{\mathrm{rel}}^{{\mathbb{T}} _{2^k}}}{k\beta}=2.\]
We complement this result by a much larger lower bound essentially matching the upper bound in Theorem 1, for a domain \(\Lambda_{d,n}\) of subcritical size. It does not seem to have been predicted in the physics studies.
Theorem 3 (Stronger bottleneck). For any \(d\geqslant 2\) there exists \(C=C(d)\geqslant 1\) such that the following holds. If \(n\leqslant e^{\beta/C}\), then \[e^{\beta n^{1/C}}/C\leqslant T_{\mathrm{rel}}^{-\Lambda_{d,n},\mathbf{1}}\leqslant Ce^{Cn^{d}(\beta+1)}.\]
We note however that, because of the subcriticality assumption, despite 56 below, Theorem 3 does not imply that the same lower bound holds in infinite volume. Indeed, we rather expect a highly unusual speed-up to take place in larger volumes and the lower bound of Theorem 1 to be closer to the truth. The result of Theorem 3 is in fact not restricted to this particular shape of domain (e.g.\(\{1,\dots,n\}^d\) would work) or boundary condition, but is specific to the size specified in the statement.
Further consequences of our study in geometric group theory and ergodic theory will be presented directly in Sections 4 and 5, once the necessary background is recalled.
We mention another class of models that have been proposed to explain the behaviour of glassy materials, the kinetically constrained models. An example in one dimension, is the East model: a site is occupied (\(-1\)) or empty (\(+1\)); and a site may switch its occupation status as in Glauber dynamics, but only when its East (right) neighbour is occupied. Aldous and Diaconis Aldous02? showed, in their pioneering work (see Faggionato13? for an overview), that this model also follows a quadratic super-Arrhenius law \(T_{\mathrm{rel}}=e^{(C+o(1))\beta^2}\) as \(\beta\to\infty\) with \(C=1/(2\ln2)\), see Cancrini08?.
Within kinetically constrained models, an important role is played by those in the so-called ‘supercritical rooted’ universality class, including the East model. In two dimensions, Marêché, Martinelli, Morris and Toninelli showed that the same \(e^{\beta^2}\) behaviour holds Martinelli19a?, Mareche20combi?; we direct the reader to overviews from a physical Ritort03?, Garrahan11? and a mathematical HartarskybookKCM? perspective.
A major drawback of these kinetically constrained models is that their constraints are imposed ad hoc to reflect the intuition of dynamical facilitation and caging effects, without any microscopic justification for the emergence of such constraints. This was an important motivation to introduce a more physically plausible model in the form of plaquettes. While widely studied in physics (nowadays, quantum versions of these models are actively investigated in physics, see Sfairopoulos23?, Sfairopoulos25?) plaquette models remained largely intractable mathematically until the work of Chleboun, Faggionato, Martinelli and Toninelli Chleboun17?. The behaviour of such models is highly dependent on the shape of the plaquette. Some models, such as the square plaquette one, present behaviour similar to the ‘supercritical unrooted’ kinetically constrained models featuring an Arrhenius scaling \(T_{\mathrm{rel}}=e^{\beta}\) and are much more accessible Chleboun21?, Chleboun20?, Espriu04?, Mueller17?.
The triangular plaquette model is also closely related to important notions in quite different fields. In ergodic theory, the Ledrappier subshift, or \(3\)-dot subshift Ledrappier78? consists in spin configurations in which every triangle contains an odd number of \(+1\); equivalently, in minimal-energy configurations also known as groundstates. It is a prominent example of \(\mathbb{Z}^2\)-subshift which is mixing but not \(2\)-mixing. We refer to Section 5 for more background and consequences of our results to the description of configurations with small support.
Within group theory, Baumslag Baumslag72? gave an example of finitely-presented metabelian group \(\Gamma\) containing the wreath product \(\mathbb{Z}\wr\mathbb{Z}=\bigoplus_{\mathbb{Z}}\mathbb{Z}\rtimes\mathbb{Z}\). A variant contains the “lamplighter group”, namely the wreath product \((\mathbb{Z}/2\mathbb{Z})\wr\mathbb{Z}\). It may be given by the presentation \[\label{eq:Baumslag} \Gamma=\left\langle a,x,y\mid xy=yx,a^2=a\cdot x^{-1}a x\cdot y^{-1}a y=1\right\rangle\tag{7}\] or as \(\Gamma=S\rtimes\mathbb{Z}^2\) where \(S\) consists in equivalence classes of finitely supported lamp configurations (in other words, those that differ by an element of the Ledrappier subshift). We refer to Section 4 for more background and consequences or our results to the word problem of \(\Gamma\).
In Section 2, we present three combinatorial results, corresponding respectively to: the lower bounds of Theorems 1 and 2; the lower bound of Theorem 3; the upper bound of Theorem 1. The first (Section 2.1) may be stated as follows in terms of lamps and switches: converting the configuration with a lone ON lamp at the origin into a configuration with no ON lamp in a ball of radius \(n\) requires, at some intermediate step, to have at least \(\log_2 n\) ON lamps in that ball; so the energy of the system must rise significantly at some intermediate step. Since this is improbable, the relaxation time is bounded from below.
The second combinatorial result (Section 2.2) asserts that there exist configurations in a region of size \(n\) which can be reduced to the constant \(+1\) configuration, but require the creation of \(n^{\Theta(1)}\) additional ON lamps to do so; we call these entangled configurations.
The third combinatorial result (Section 2.3) considers cycles: switch configurations whose associated lamp configuration has all lamps OFF in a given region. We obtain estimates on a generating function counting cycles; namely, we show that, in a suitably large simplex, there are not too many cycles, when weighted by their number of \(-1\) switches. It follows that the partition functions associated with two different completions outside a large simplex differ very little.
In Section 3, we provide the probabilistic arguments that make use of the combinatorial results mentioned above. The proof of Theorems 1, 2 and 3 is completed there. Note however that the upper bounds of Theorems 2 and 3 are proved directly (by the bisection technique or canonical paths) without input from Section 2.
Finally, in Sections 4 and 5, we elaborate on the applications of the combinatorial results to group theory and subshifts.
In this section we focus on three purely deterministic statements of combinatorial nature which will be the key to the main results.
In this subsection, we only consider the lamp configurations; every lamp may be ON (\(1\)) or OFF (\(0\)). We call elements of \(\{0,1\}^{\mathbb{Z}^d}\) simply configurations, and we identify them with subsets of \(\mathbb{Z}^d\) by considering their support. We add configurations pointwise modulo \(2\); this corresponds to symmetric difference.
Recall the plaquette \(T_d=\{0,\, e_1,\, e_2,\dots,e_d\}\) from 1 . In one move we may choose a point \(x\in\mathbb{Z}^d\) and toggle the lamps at all points of the set \(x+T_d\) (i.e.add its indicator function modulo \(2\)). Any shift \(x+T_d\) of this set is called a plaquette.
Definition 1 (Chain). A sequence of configurations \((W_i)_{i=0}^n\) is called a chain if for every \(i\) the symmetric difference \[W_i\mathbin{\mathpalette\triangle\relax}W_{i+1}\] is a plaquette or \(\varnothing\).
We will use the closed \(\ell_{\infty}\)-ball in \(\mathbb{Z}^d\): \[\label{eq:ball} B(x,R) \mathrel{\vcenter{:}}= \left\{y\in\mathbb{Z}^d\mid \|x-y\|_{\infty} \leqslant R\right\}.\tag{8}\]
Theorem 4 (Combinatorial bottleneck). Denote the origin of \(\mathbb{Z}^d\) by \(O\). Consider \(k\in\mathbb{N}\), and set \(B_k = B(O,3d\,2^{k})\). Let \((W_0,\dots,W_n)\) be a chain such that \(W_0=\{O\}\) and \(W_n\cap B_k=\varnothing\). Then for some index \(i\) we have \[|W_i\cap B_k| > k.\]
We note that this result is sharp up to an additive constant in the last display when \(d=2\) and up to a multiplicative one for \(d>2\). This was already noted for \(d=2\) in Newman99?, and can be easily seen by induction, using the configuration of Figure 4 and 38 . Since we will not use this fact, we omit the details.
For \(x\in\mathbb{Z}^d\) define the inverted plaquette \[\label{eq:def:inverted} K_x \mathrel{\vcenter{:}}= 2x - T_d = \left\{2x,\, 2x-e_1,\, 2x-e_2, \dots, 2x-e_d\right\}.\tag{9}\]
Lemma 1 (Types of inverted plaquettes).
If \(x\neq y\) are in \(\mathbb{Z}^d\), then \(K_x\cap K_y=\varnothing\).
A plaquette intersects inverted plaquettes in exactly one of the following ways:
it intersects none;
it intersects exactly one inverted plaquette and does so in exactly two points;
it has the form \(2x+T_d\) for some \(x\in\mathbb{Z}^d\) and intersects exactly the \(d+1\) inverted plaquettes \(K_x,K_{x+e_1},K_{x+e_2},\dots, K_{x+e_d}\), each in one point.
Proof. [item:1]. If \(2x-t = 2y-s\) for some \(t,s\in T_d\), then \(2(x-y)=t-s\). But \(t-s\) has coordinates in \(\{-1,0,1\}\), hence cannot equal a nonzero even vector (that is, a vector in \(2\mathbb{Z}^d\)). Therefore \(t=s\) and \(x=y\).
[item:2]. Modulo \(2\) we have \(2y\equiv 0\) and \(2y-e_i\equiv e_i\), hence \[p\in\bigcup_y K_y \iff p \bmod 2 \in T_d.\]
For a plaquette \(B=x+T_d\) the \((d+1)\) residues are \((x\bmod 2)+T_d\), so the number of points of \(B\) lying in \(\bigcup_y K_y\) is \((d+1)\) iff \(x\) is even, \(2\) iff \(x\) has exactly one or two odd coordinates and otherwise \(0\).
Thus we get cases [item:intersection:c], [item:intersection:b], [item:intersection:a]. If \(x=2u{\in2\mathbb{Z}^d}\), then \[2u\in K_u,\qquad 2u+e_i\in K_{u+e_i}\;(i=1,\dots, d),\] so \(2u+T_d\) meets \(K_u,K_{u+e_1},\dots,K_{u+e_d}\). By [item:1]. of the proof, the intersections of \(2x+T_d\) with each of these inverted plaquettes are disjoint, so they consist of one point each.
If \(x = 2u-e_i\), then \(x+T_d\) has two common points \(\{2u,2u-e_i\}\) with \(K_{u}\). By [item:1]. of the proof, \(x+T_d\) has no intersection with any \(K_v\) for \(v\in\mathbb{Z}^d\setminus\{u\}\).
If \(x=2u-e_i-e_j\) for \(i\neq j\), then \(x+T_d\) has two common points \(\{2u-e_i,2u-e_j\}\) with \(K_u\). By [item:1]. of the proof, there \(x+T_d\) has no intersection with any \(K_v\) for \(v\in\mathbb{Z}^d\setminus\{u\}\).
In the other cases, \(x\) and \(x+e_i\) have at least two odd coordinates, so cannot belong to an inverted plaquette. ◻
We now define the renormalisation map \(\mathcal{R}\) from configurations to configurations as illustrated in Figure 2. Given a configuration \(W\), define \(\mathcal{R}(W)\) at a point \(x\in\mathbb{Z}^d\) as the parity of the number of ON lamps of \(W\) inside \(K_x\): \[\label{eq:renorm} \mathcal{R}(W)(x) \mathrel{\vcenter{:}}= \sum_{p\in K_x} W(p)\pmod 2.\tag{10}\] (this operator already appears, for \(d=2\), in the proof of Arenas-Carmona08?*Theorem 7.1). From Lemma 1 we obtain the following key fact.
Lemma 2 (Renormalisation). If \((W_0,\dots,W_n)\) is a chain, then \(\bigl(\mathcal{R}(W_0),\dots,\mathcal{R}(W_n)\bigr)\) is also a chain.0◻
Let us note that Lemma 2 applies without change to any plaquette shape instead of \(T_d\), although we will not need this fact.
In Lemma 1, we classified the possible intersection types between plaquettes and inverted plaquettes. We will now look at intersections with simplices, recalling 6 .
Lemma 3.
A plaquette is either entirely contained in a simplex, or it intersects the simplex in at most one point.
Two distinct plaquettes intersect in at most one point.
Proof. [item:intersect2:1]. It suffices to consider a simplex of the form \(\Lambda_{d,m}\). Suppose \(T_d+x\) has a nonempty intersection with \(\Lambda_{d,m}\) but is not contained in \(\Lambda_{d,m}\). There are two cases.
Assume \(x\in \Lambda_{d,m}\). Then \(x_1+\dots+x_d = m\) (otherwise \(x+e_i\in \Lambda_{d,m}\) for all \(i\)), and therefore \[(T_d+x)\cap\Lambda_{d,m} = \{x\}.\]
Assume \(x\notin \Lambda_{d,m}\) and \(x+e_i\in \Lambda_{d,m}\) for some \(i\). Then \(x_i = -1\), and \[(T_d+x)\cap \Lambda_{d,m} = \{x+e_i\}.\]
[item:intersect2:2]. This is immediate from [item:intersect2:1]., since the plaquette \(T_d+x\) is precisely \(x+\Lambda_{d,1}\). ◻
We construct recursively a sequence of nested simplices \(U_k\). We start with \(U_0\mathrel{\vcenter{:}}=\{O\}\). Assuming that the simplex \(U_k\) has been constructed, define \[U'_{k+1}\mathrel{\vcenter{:}}= \bigcup_{x\in U_{k}}K_x.\] Let \(U_{k+1}\) be the smallest simplex that contains, in full, every plaquette that intersects \(U'_{k+1}\).
Proposition 5 (Inductive statement of the combinatorial bottleneck). Let \((W_0,\dots,W_n)\) be a chain such that \(W_0=\{O\}\) and \(W_n\cap U_k=\varnothing\). Then for some index \(i\), we have \[|W_i\cap U_k| > k.\]
Before proving this result, let us deduce Theorem 4 from it.
Proof of Theorem 4. If \(U_k=x+\Lambda_{d,m}\) then \(U_{k+1}=x'+\Lambda_{d,m'}\) with \(m'=2m+3d\) and \(x'=2x-2(e_1+\dots+e_d)\). We obtain by induction \(U_k=\Lambda_{d,3d(2^k-1)}-2(2^k-1)(e_1+\dots+e_d)\subset B_k\), and we are done by Proposition 5. ◻
Definition 2 (Good chain). Given \(k\in\mathbb{N}\), a chain \(\mathbb{W} = (W_0,\dots,W_n)\) is good, if \(W_0=\{O\}\), \(W_n\cap U_k=\varnothing\) and \[|W_i\cap U_k|\leqslant k\] for all \(i=0,\dots,n\).
Definition 3 (Minimal good chain). Given \(k\in\mathbb{N}\), a good chain \((W_0,\dots,W_n)\) is minimal, if, for any other good chain \((W'_0,\dots,W'_{n'})\), one of the following holds: \[\begin{align} n&{}< n';\\ n&{}=n',&\sum_{i=0}^n|W_i\cap U_k|&{}<\sum_{i=0}^{n}|W'_i\cap U_k|;\\ n&{}=n',&\sum_{i=0}^n|W_i\cap U_k|&{}=\sum_{i=0}^{n}|W'_i\cap U_k|,&\sum_{i=0}^n|\mathcal{R}(W_i)\cap U_{k-1}|&{}\leqslant\sum_{i=0}^n|\mathcal{R}(W'_i)\cap U_{k-1}|. \end{align}\]
Definition 4 (Internal inverted plaquette). Given \(k\in\mathbb{N}\), if \(u\in U_{k-1}\), we call the inverted plaquette \(K_u\) internal. In what follows, we abbreviate “internal inverted plaquette” as iip.
Recall that the union of all iip is the set \(U'_k\subset U_k\).
Definition 5 (Critical configuration). Given \(k\in\mathbb{N}\), we call a configuration \(W\) critical if inside \(U'_k\) there are exactly \(k\) iip, each containing exactly one ON lamp, and there are no other ON lamps in \(U_k\).
Clearly, a configuration is critical if and only if \(|W\cap U_k|=|\mathcal{R}(W)\cap U_{k-1}| = k\).
Definition 6 (Spread plaquette). Given \(k\in\mathbb{N}\), we say that a plaquette is spread if it is contained in \(U'_k\).
Note that a spread plaquette is necessarily of type [item:intersection:c] in Lemma 1 (i.e.it intersects \((d+1)\) inverted plaquettes).
Lemma 4. Let \(k\in\mathbb{N}\) and let \(W\) be a critical configuration, let \(\mathcal{T}\) be a plaquette, and assume \(W\cap U_k\neq (W\mathbin{\mathpalette\triangle\relax}\mathcal{T})\cap U_k\) and \(|(W\mathbin{\mathpalette\triangle\relax}\mathcal{T})\cap U_k|\leqslant k\). Then \(\mathcal{T}\) is a spread plaquette.
Proof. \(\mathcal{T}\) intersects \(U_k\), so toggling \(\mathcal{T}\) must turn off at least one ON lamp of \(W\) in \(U_k\); otherwise we would only turn lamps on and hence would get \(|(W\mathbin{\mathpalette\triangle\relax}\mathcal{T})\cap U_k|>k\). Thus, \(\mathcal{T}\) intersects \(U'_{k}\), because \(W\subset U'_k\). Then, by the definition of \(U_k\), we get \(\mathcal{T}\subset U_k\).
Since \(d+1\geqslant 3\), applying \(\mathcal{T}\) must switch off at least two ON lamps of \(W\). Each inverted plaquette contains at most one ON lamp, so \(\mathcal{T}\) intersects at least two iip. This forces case [item:intersection:c] of Lemma 1: applying \(\mathcal{T}\) changes lamps in \((d+1)\) inverted plaquettes of the form \(K_x,K_{x+e_1},K_{x+e_2},\dots, K_{x+d}\).
Since \(U_{k-1}\) is a simplex, Lemma 3 implies \[|(x+T_d)\cap U_{k-1}| \in\{0, 1,d+1\}.\] The only case compatible with our assumptions is the third one. In this case, all \((d+1)\) inverted plaquettes intersecting \(\mathcal{T}\) are iip. ◻
Lemma 5 (Minimal good chains around a critical configuration). Let \(k\in\mathbb{N}\) and let \(\mathbb{W}=(W_0, W_1, \dots, W_n)\) be a minimal good chain. Let \(0 < i < n\) be an index such that \(|\mathcal{R}(W_i)\cap U_{k-1}| = k\). Set \(\mathcal{T}_1\mathrel{\vcenter{:}}= W_i\mathbin{\mathpalette\triangle\relax}W_{i-1}\) and \(\mathcal{T}_2\mathrel{\vcenter{:}}= W_{i+1}\mathbin{\mathpalette\triangle\relax}W_i\), and set \(W'\mathrel{\vcenter{:}}= W_{i-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_2\). Then
\(\mathcal{T}_1\) and \(\mathcal{T}_2\) are disjoint spread plaquettes;
\(|W_{i-1}\cap U_k| = |W_{i+1}\cap U_k| = k\);
the configurations \(W_{i\pm 1}\) contain no ON lamps in \(U_k\setminus U'_k\);
the configuration \(W'\) is also critical;
the sequence \((W_0, W_1, \dots, W_{i-1}, W', W_{i+1}, \dots, W_n)\) is also a minimal good chain.
Proof. Clearly, \(W_i\) is critical.
\(\mathcal{T}_1\) and \(\mathcal{T}_2\) both intersect \(U_k\); otherwise we could shorten the chain. By Lemma 4, it follows that \(\mathcal{T}_1\) and \(\mathcal{T}_2\) are spread plaquettes. In particular, they are contained in \(U'_k\), which proves 3.
If \(\mathcal{T}_1 = \mathcal{T}_2\), then \(W_{i-1} = W_{i+1}\) and we can again shorten the chain. Therefore, \(\mathcal{T}_1\) and \(\mathcal{T}_2\) are disjoint, proving 1. Indeed, distinct spread plaquettes are disjoint by [item:1]. of Lemma 1, which applies to spread plaquettes by symmetry.
We may apply \(\mathcal{T}_1\) and \(\mathcal{T}_2\) in the opposite order: consider the chain \[\mathbb{W}'\mathrel{\vcenter{:}}=(W_0, W_1, \dots, W_{i-1}, W', W_{i+1}, \dots, W_n).\] If \(|W'\cap U_k| < k\), then \(\mathbb{W}'\) is a good chain and this contradicts the minimality of \(\mathbb{W}\). Therefore, \(|W'\cap U_k| \geqslant k\). Since \(\mathcal{T}_1\) and \(\mathcal{T}_2\) are disjoint, \[|W_i\cap U_k| + |W'\cap U_k| = |W_{i-1}\cap U_k| + |W_{i+1}\cap U_k|.\] The left-hand side is at least \(2k\), while the right-hand side is at most \(2k\). Therefore, the only possibility is \[|W_{i-1}\cap U_k|=|W_{i+1}\cap U_k|=|W'\cap U_k|=k.\] This proves 2.
Next, \(|\mathcal{R}(W')\cap U_{k-1}| \leqslant|W'\cap U_{k}| = k\), so \(\mathbb{W}'\) is also a minimal good chain, proving 5. If \(|\mathcal{R}(W')\cap U_{k-1}| < |W'\cap U_{k}|\), we again contradict the minimality of \(\mathbb{W}\). Hence, \(|\mathcal{R}(W')\cap U_{k-1}| = k\), so \(W'\) is critical, proving 4. ◻
Proof of Proposition 5. We proceed by induction on \(k\). The base cases \(k=0,1\) are straightforward, using the assumption \(d\geqslant 2\). Fix \(k>1\). Assume the statement holds for \(k-1\), and let us prove it for \(k\).
Our goal is to show that no good chain exists. Suppose, for the sake of contradiction, that there exists a good chain \(\mathbb{W} = (W_0,\dots,W_n)\) that we may assume minimal.
Define \[V_i \mathrel{\vcenter{:}}= \mathcal{R}(W_i).\] Then \(V_0\) consists of the single point \(O\). Moreover, \(V_n\cap U_{k-1}=\varnothing\), because if \(x\in U_{k-1}\) then \(K_x\subseteq U'_k\subset U_k\), so \(W_n\) has no ON lamps in \(K_x\).
By Lemma 2, the sequence \((V_0,\dots,V_n)\) is a chain. By the induction hypothesis (applied to the chain \((V_0,\dots,V_n)\) and size parameter \(k-1\)), there exists an index \(i\) such that \[|V_i\cap U_{k-1}| \geqslant k.\] Then \(|V_i\cap U_{k-1}| = k\), as otherwise \(|W_i\cap U_k| > k\).
Since \(|V_0\cap U_{k-1}|=1\) and \(V_n\cap U_{k-1}=\varnothing\), we can choose an index \(t\) such that \[|V_{t-1}\cap U_{k-1}| < k \quad \text{and} \quad |V_t\cap U_{k-1}| = k.\] Since \(|W_t\cap U_k| \geqslant|V_t\cap U_{k-1}|\) and \(|W_t\cap U_k| \leqslant k\), the only possibility is \(|W_t\cap U_k| = k\).
Set \(\mathcal{T}_1\mathrel{\vcenter{:}}= W_t\mathbin{\mathpalette\triangle\relax}W_{t-1}\) and \(\mathcal{T}_2\mathrel{\vcenter{:}}= W_t\mathbin{\mathpalette\triangle\relax}W_{t+1}\). Then, by Lemma 5, \(\mathcal{T}_1\) and \(\mathcal{T}_2\) are spread plaquettes, and moreover \[|W_{t-1}\cap U_k|=|W_{t+1}\cap U_k|= k.\]
If the dimension \(d\) is even, we immediately obtain a contradiction with \(|W_t\cap U_k|=k\), since the plaquette toggles an odd number of lamps. For odd \(d\), we need a more delicate argument.
As in Lemma 5, set \(W'\mathrel{\vcenter{:}}= W_{t-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_2\). By Lemma 5, we have \(|W'\cap U_k| = k\). We will contradict the minimality of the original chain if we can show that \(|\mathcal{R}(W')\cap U_{k-1}| < |W'\cap U_k|\).
Since \(W_t\) is critical, all its ON lamps lie inside iip. By Lemma 4, \(\mathcal{T}_1\) contains no points outside iip. Therefore, all ON lamps of \(W_{t-1} = W_t\mathbin{\mathpalette\triangle\relax}\mathcal{T}_1\) also lie inside iip. At the same time, we know that \(|V_{t-1}\cap U_{k-1}| < |W_{t-1}\cap U_k|\). Hence, in some iip \(K_u\) the configuration \(W_{t-1}\) has at least two ON lamps. It is easy to see that there are exactly two, since applying \(\mathcal{T}_1\) changes the number of ON lamps in \(K_u\) by exactly \(1\), as \(\mathcal{T}_1\) is a spread plaquette touching \(K_u\).
For \(i=1,2\) set \(\mathcal{S}_i=\{x\in\mathbb{Z}^d:K_x\cap\mathcal{T}_i\neq\varnothing\}\), and note by Lemma 1 that each \(\mathcal{S}_i\) is a plaquette. Since \(W_t = W_{t-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_1\) is critical, we have \(u \in \mathcal{S}_1\); otherwise \(W_t\) would also contain two ON lamps inside \(K_u\). From Lemma 5 we know that \(W' = W_{t-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_{2}\) is also critical. Thus, by the same reasoning, \(u\in \mathcal{S}_2\).
Therefore, \(W_{t-1}\) contains two ON lamps inside \(K_u\); denote them by \(\ell_1\in\mathbb{Z}^d\) and \(\ell_2\in\mathbb{Z}^d\). Up to reordering, we may assume \(\ell_i\) lies in \(\mathcal{T}_i\) for \(i=1,2\), and there are no other intersections of \(\mathcal{T}_i\) with \(K_u\). Then \(W_t\) contains exactly one ON lamp inside \(K_u\), namely \(\ell_2\), while \(W'\cap K_u=\{\ell_1\}\), and \(W_{t+1}\) contains no ON lamps inside \(K_u\).
We now consider two cases.
Case 1. Assume that the configuration \(W_{t+1}\) is not critical. We apply to \(W_{t+1}\) the same reasoning as above for \(W_{t-1}\). The configuration \(W_{t+1}\) has exactly \(k\) ON lamps inside \(U'_k\), so in some iip \(K_v\) two lamps must be lit. As before, we obtain \(v\in \mathcal{S}_1\cap \mathcal{S}_2\). Moreover, \(v\neq u\), because \(W_{t+1}\) has no ON lamps inside \(K_u\). Therefore \(|\mathcal{S}_1\cap\mathcal{S}_2|\geqslant 2\), contradicting Lemma 3.
Case 2. Assume that the configuration \(W_{t+1}\) is critical. Set \[\mathcal{T}_3\mathrel{\vcenter{:}}= W_{t+2}\mathbin{\mathpalette\triangle\relax}W_{t+1}.\] Apply Lemma 5 at index \(i = t+1\). We obtain \(\mathcal{T}_2\neq \mathcal{T}_3\), and a minimal good chain \[\mathbb{W}''\mathrel{\vcenter{:}}=(W_0,\dots, W_{t-1}, W_t, W_t\mathbin{\mathpalette\triangle\relax}\mathcal{T}_3, W_{t+2}, \dots, W_n)\] Apply Lemma 5 again, this time to the sequence \(\mathbb{W}''\) at index \(t\). We obtain yet another minimal good chain \[\mathbb{W}'''\mathrel{\vcenter{:}}=(W_0,\dots, W_{t-1}, W_{t-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_3, W_t\mathbin{\mathpalette\triangle\relax}\mathcal{T}_3, W_{t+2}, \dots, W_n).\]
In addition, \(\mathcal{T}_1\neq \mathcal{T}_3\) (otherwise \(\mathbb{W}\) would not have been minimal), and the configuration \(W_{t-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_3\) is critical. Recall that \(W_{t-1}\) has ON lamps \(\ell_i\in K_u\cap\mathcal{T}_i\), for \(i=1,2\). Neither of these lamps lies in \(\mathcal{T}_3\), since distinct spread plaquettes are disjoint. Therefore, the configuration \(W_{t-1}\mathbin{\mathpalette\triangle\relax}\mathcal{T}_3\) cannot be critical, because it still contains at least two ON lamps inside \(K_u\).
This contradiction completes the inductive step and hence the proof of the proposition. ◻
Recall that a plaquette is any translate \(x+T_d\) of the set \(T_d=\{0,e_1,e_2,\dots,e_d\}\).
Definition 7 (Admissible configuration). A set \(W\subset \mathbb{Z}^d\) is called admissible if it is the sum modulo \(2\) of finitely many plaquettes.
The purpose of this subsection is to produce small, admissible configurations that require a large number of additional ON lamps to be turned completely off. We call these entangled configurations.
Admissible sets have a simple algebraic interpretation. Let \[R_d \mathrel{\vcenter{:}}= \mathbb{F}_2\left[x_1^{\pm 1},\ldots,x_d^{\pm 1}\right]\] be the ring of Laurent polynomials over \(\mathbb{F}_2\) in the variables \(x_1,\ldots,x_d\). Its elements are finite sums of the form \[\sum_{\alpha \in \mathbb{Z}^d} c_\alpha x^\alpha, \qquad c_\alpha \in \mathbb{F}_2, \qquad x^\alpha = x_1^{\alpha_1}\cdots x_d^{\alpha_d},\] where only finitely many coefficients \(c_\alpha\) are nonzero.
Proposition 6. A finite set \(W \subset \mathbb{Z}^d\) is admissible if and only if the Laurent polynomial \[\sum_{(t_1,\dots,t_d)\in W} \prod_{i=1}^d x_i^{t_i}\] is divisible in \(R_d\) by \(1+x_1+\cdots+x_d\).
Proof. This follows directly from the fact that a translate \(a+T_d\) corresponds to the Laurent monomial \(x^a\) multiplied by \(1+x_1+\cdots+x_d\), while summation of configurations modulo \(2\) corresponds to addition in \(R_d\). ◻
Thus, under this correspondence, finite configurations are represented by Laurent polynomials, and admissible configurations form the principal ideal generated by \(1+x_1+\cdots+x_d\). We shall not need any deeper algebraic properties of this representation.
Let \(k\geqslant 0\) be an integer. By the Frobenius endomorphism, \[(1+x_1+\dots+x_d)^{2^k} = 1+x_1^{2^k}+\dots+x_d^{2^k}.\] Therefore the set \[2^kT_d = \left\{0,2^ke_1,\dots, 2^ke_d\right\}\] is also admissible. We shall call these sets, and all their translates, large plaquettes.
For \(d=1\), an admissible set is simply any finite set of even cardinality. In higher dimensions, admissible sets are much more rigid. In particular, admissible sets of small cardinality can be completely described.
Lemma 6. Let \(d \geqslant 2\), and let \(X\subset \mathbb{Z}^d\) be an admissible set such that \[0 < |X| \leqslant d+1.\] Then \(X\) is a large plaquette.
Proof. We prove the statement by induction on \(d\).
Base case: \(d=2\). This case is proved, for instance, in Arenas-Carmona08?*Lemma 5.6. For the reader’s convenience, we recall one possible argument.
The Newton polytope of a Laurent polynomial is the convex hull of its support. Under multiplication of Laurent polynomials, Newton polytopes add by Minkowski summation. The Newton polytope of \(1+x_1+x_2\) is the triangle \(\operatorname{conv}(T_2)\). Hence Proposition 6 implies that every admissible set \(X\subset\mathbb{Z}^2\) has
an edge with outer normal \(-e_2\);
an edge with outer normal \(-e_1\);
an edge with outer normal \(e_1+e_2\).
It follows at once that if \(|X|\leqslant 3\), then \(X\) is homothetic to \(T_2\) with an integral scaling factor. After translating \(X\), we may assume that \[X = \{(0,0), (a,0), (0,a)\}\] for some positive integer \(a\). It remains to show that \(a\) is a power of \(2\). Set \[P(x_1,x_2)=1+x_1^a+x_2^a.\] Since \(P\) is divisible by \(1+x_1+x_2\), the polynomial \(P(x_1,1+x_1)\) must vanish identically. Over \(\mathbb{F}_2\), \[0=P(x_1,1+x_1) = 1+x_1^a+(1+x_1)^a = \sum_{i=1}^{a-1}\binom{a}{i}x_1^i.\] Thus all intermediate binomial coefficients \(\binom{a}{i}\), \(1\leqslant i\leqslant a-1\), are even. By Kummer’s theorem, equivalently by Lucas’ theorem, this is possible only when \(a\) is a power of \(2\). Hence \(X\) is a large plaquette.
Induction step. Assume that the lemma has already been proved in dimension \(d-1\), we prove it in dimension \(d\).
Let \(X\subset\mathbb{Z}^d\) be admissible. Translating \(X\), we may assume that the minimum first coordinate among points of \(X\) is \(0\). Consider the Laurent polynomial \[P(x_1,\dots,x_d) = \sum_{(t_1,\dots,t_d)\in X}\prod_{i=1}^d x_i^{t_i}.\] By Proposition 6, \[P(x_1,\dots,x_d) = (1+x_1+\dots+x_d)Q(x_1,\dots,x_d)\] for some Laurent polynomial \(Q\).
Let \(b\) be the minimum exponent of \(x_1\) among the nonzero monomials of \(Q\). Since the minimum exponent of \(x_1\) among the monomials of \(P\) is \(0\), we have \(b=0\). Write \[Q=Q_1+Q_2,\] where all monomials of \(Q_1\) have \(x_1\)-degree \(0\), and all monomials of \(Q_2\) have strictly positive \(x_1\)-degree. Then \[P_0\mathrel{\vcenter{:}}= Q_1(1+x_2+x_3+\dots+x_d)\] is precisely the sum of all monomials of \(P\) with \(x_1\)-degree \(0\). Therefore \(P_0\) has at most \(d+1\) monomials. If \(P_0\) had exactly \(d+1\) monomials, then, since \(|X|\leqslant d+1\), we would have \(P=P_0\). This is impossible: the Newton polytope of \(P_0\) has dimension at most \(d-1\), whereas every nonzero multiple of \(1+x_1+\dots+x_d\) has \(d\)-dimensional Newton polytope.
Thus \(P_0\) has at most \(d\) monomials. The support of \(P_0\) is an admissible set in \(\mathbb{Z}^{d-1}\), so by the induction hypothesis, \[P_0 = x_2^{t_2}\cdots x_d^{t_d} \left(1+x_2^{2^k}+\dots+x_d^{2^k}\right)\] for some integers \(t_2,\dots,t_d\) and some integer \(k\geqslant 0\). In particular, \(P_0\) has exactly \(d\) nonzero monomials, and \(P\) has exactly \(d+1\) nonzero monomials. At this point we have also proved that, in dimension \(d\), there are no nonempty admissible sets of cardinality strictly smaller than \(d+1\).
| Now consider $$P |
|---|
| x_2{t_2}x_d{t_d} |
| (1+x_1{2k}+x_2{2k}++x_d{2k}).$$ This Laurent |
| polynomial is divisible by \(1+x_1+\dots+x_d\). On the other hand, it has |
| at most two nonzero monomials. Since there are no nonempty admissible |
| sets of cardinality smaller than \(d+1\), this polynomial must be zero. |
| Hence $$P |
| = |
| x_2{t_2}x_d{t_d} |
| (1+x_1{2k}++x_d{2k}),$$ so \(X\) is a large plaquette. |
| The induction step is complete. ◻ |
Let \(d\in\mathbb{N}\). We call a vector long if it has one of the forms \[\pm 2^k e_i \qquad\text{or}\qquad \pm 2^k(e_i-e_j),\] where \(k\geqslant 0\) is an integer and \(1\leqslant i<j\leqslant d\). Equivalently, long vectors are precisely directed edges of large plaquettes.
Define a metric \(\rho\) on \(\mathbb{Z}^d\) as follows: \(\rho(x,y)\) is the minimum number of long vectors whose sum is \(x-y\).
Lemma 7 (Exponential decay in \(\rho\)). Let \[\lambda_d=\frac{d}{d+1}.\] Let \(X\subset\mathbb{Z}^d\) be a nonempty admissible set, and let \(\alpha\in X\). Then \[\label{eq:rho95sum} \sum_{x\in X\setminus\{\alpha\}} \lambda_d^{\rho(x,\alpha)-1} \geqslant 2.\tag{11}\]
Proof. We argue by induction on \(|X|\).
By Lemma 6, the smallest possible cardinality of a nonempty admissible set is \(d+1\). In that case \(X\) is a large plaquette. If \(\alpha\in X\), then all other \(d\) points of \(X\) are at \(\rho\)-distance \(1\) from \(\alpha\), and the left-hand side of 11 is exactly \(d\). This proves the base case.
Now assume that \(|X|>d+1\). Translating \(X\), we may suppose \(\alpha=0\in X\). We also assume that \(X\) contains a point with at least one odd coordinate. Indeed, if all coordinates of all points of \(X\) are even, then we may divide all coordinates by their largest common power of \(2\). The resulting set remains admissible, and the left-hand side of 11 does not increase.
Call the following \(d(d+1)\) points near: \[\alpha\pm e_i, \qquad \alpha\pm(e_i-e_j), \qquad 1\leqslant i<j\leqslant d.\] These points have \(\rho\)-distance \(1\) from \(\alpha\).
Case 1. Suppose that at least two near points belong to \(X\). Each of them contributes \(1\) to the sum in 11 , so the desired inequality follows immediately.
Case 2. Suppose that \(d=2\) and exactly one near point belongs to \(X\). Denote this point by \(\beta\). Choose a point \(\gamma\) such that \[\{\alpha,\beta,\gamma\}\] is a plaquette. Then \(\gamma\) is also near, and \(\gamma\notin X\).
Set \[X' = X \mathbin{\mathpalette\triangle\relax}\{\alpha,\beta,\gamma\},\] where \(\mathbin{\mathpalette\triangle\relax}\) denotes symmetric difference. Then \(X'\) is admissible, \(|X'|<|X|\), and \(\gamma\in X'\). By the induction hypothesis, \[\sum_{x\in X\setminus\{\alpha,\beta\}} \lambda_2^{\rho(x,\gamma)-1} = \sum_{y\in X'\setminus\{\gamma\}} \lambda_2^{\rho(y,\gamma)-1} \geqslant 2.\] Since \(\rho(\gamma,\alpha)=1\), we have \[\rho(x,\alpha)\leqslant\rho(x,\gamma)+1\] for every \(x\). Hence, \[\sum_{x\in X\setminus\{\alpha\}} \lambda_2^{\rho(x,\alpha)-1}\geqslant 1+ \sum_{x\in X\setminus\{\alpha,\beta\}} \lambda_2^{\rho(x,\alpha)-1} \geqslant 1+ \lambda_2 \sum_{x\in X\setminus\{\alpha,\beta\}} \lambda_2^{\rho(x,\gamma)-1} \geqslant 1+2\lambda_2 >2.\]
Case 3. Suppose either that \(d=2\) and no near point belongs to \(X\), or that \(d>2\) and at most one near point belongs to \(X\). If such a near point exists, denote it by \(\beta\).
Let \(\mathcal{I}=\{0,1\}^d\) be the set of all binary strings of length \(d\). Partition \(\mathbb{Z}^d\) into \(2^d\) parity classes, denoted by \(\mathbb{Z}_I^d\) for all \(I\in\mathcal{I}\). Put \[X_I \mathrel{\vcenter{:}}= (X\setminus\{\alpha=0\})\cap \mathbb{Z}_I^d.\] For \(I\in\mathcal{I}\), define \[S_I \mathrel{\vcenter{:}}= \sum_{x\in X_I} \lambda_d^{\rho(x,0)-1}.\] Define \(\mathbf{0}=00\dots 0\in\mathcal{I}\). Among the sums \(S_I\) with \(I\neq \mathbf{0}\), choose a maximal one and denote it by \(S_J\). If \(S_J=0\), then all points of \(X\) have all coordinates even, contrary to our reduction above. Therefore \(S_J>0\).
Claim 7. There exists \(\gamma\in T_d\) such that \(\gamma-T_d\) contains no point of \(\mathbb{Z}_J^d\) and \(\gamma-T_d\) contains no point of \(X\) other than possibly \(\alpha=0\).
Proof. There are \(d+1\) possible choices for \(\gamma\). For every such \(\gamma\), all points of \(\gamma-T_d\) are either \(0\) or near points. Among all the \(d(d+1)\) near points, at most two belong to the parity class \(\mathbb{Z}_J^d\). Moreover, by the assumptions of Case 3, at most one near point belongs to \(X\). Thus there are at most three bad near points when \(d>2\), and at most two when \(d=2\).
It remains to observe that each bad point rules out at most one choice of \(\gamma\in T_d\). Indeed, if the same bad point belonged to both \(\gamma_1-T_d\) and \(\gamma_2-T_d\), then the two distinct translates \(\gamma_1-T_d\) and \(\gamma_2-T_d\) would intersect in both that bad point and the origin. Equivalently, two distinct translates of \(T_d\) would intersect in two points, which is impossible because of Lemma 3. Hence at least one admissible choice of \(\gamma\) exists. ◻
Choose \(\gamma\) as in Claim 7. Let \(\mathcal{I}_\gamma\subset\mathcal{I}\) be the set of parity classes represented by the points of \(\gamma-T_d\). This set consists of the parity vector of \(\gamma\) and the \(d\) parity vectors obtained from it by flipping exactly one coordinate. In particular, \(|\mathcal{I}_\gamma|=d+1\), and \(\mathbf{0}\in\mathcal{I}_\gamma\). Also, by construction, \(J\notin\mathcal{I}_\gamma\).
Recalling the renormalisation map \(\mathcal{R}\) from 10 , set \[Y=\mathcal{R}(X-\gamma).\] By Lemma 2, \(Y\) is admissible. Since \[|X\cap(\gamma-T_d)|=1,\] we have \(0\in Y\).
For every \(y\in Y\), \[\gamma+2y-T_d\] contains an odd number of points of \(X\). Choose one of these points and denote it by \(\tau(y)\). The sets \(\gamma+2y-T_d\) are pairwise disjoint for distinct \(y\) by Lemma 1, so the points \(\tau(y)\) are pairwise distinct.
For \(I\in\mathcal{I}_\gamma\), define \(Y_I=\{y\in Y\setminus\{0\}\mid\tau(y)\in X_I\}\), so we have \[|Y_I|\leqslant|X_I|\] for every \(I\in\mathcal{I}\). Moreover, \(Y_I\) is empty unless \(I\in\mathcal{I}_\gamma\). In particular, \(Y_J\) is empty. Since \(S_J>0\), we have \(|X_J|>0\), and therefore \(|Y|<|X|\). Thus the induction hypothesis applies to \(Y\): \[\sum_{y\in Y\setminus\{0\}} \lambda_d^{\rho(y,0)-1} \geqslant 2.\]
We now compare the contributions of \(Y\) and \(X\). If \(\tau(y)\in X_{\mathbf{0}}\), then \[\rho(\tau(y),0)\leqslant\rho(y,0).\] If \(\tau(y)\notin X_{\mathbf{0}}\), then \[\rho(\tau(y),0)\leqslant\rho(y,0)+1.\] Consequently, \[S_{\mathbf{0}} \geqslant \sum_{y\in Y_{\mathbf{0}}} \lambda_d^{\rho(y,0)-1},\] and for every \(I\in\mathcal{I}_\gamma\setminus\{\mathbf{0}\}\), \[S_I \geqslant \lambda_d \sum_{y\in Y_I} \lambda_d^{\rho(y,0)-1}.\]
Since \(S_J\) is maximal among all \(S_I\) with \(I\neq\mathbf{0}\), and since \(J\notin\mathcal{I}_\gamma\), we get \[\begin{align} \sum_{x\in X\setminus\{\alpha\}} \lambda_d^{\rho(x,\alpha)-1} &\geqslant \sum_{I\in\mathcal{I}_\gamma}S_I+S_J \\ &\geqslant S_{\mathbf{0}} + \frac{d+1}{d} \sum_{I\in\mathcal{I}_\gamma\setminus\{\mathbf{0}\}}S_I \\ &\geqslant \sum_{y\in Y_{\mathbf{0}}} \lambda_d^{\rho(y,0)-1} + \frac{d+1}{d} \sum_{I\in\mathcal{I}_\gamma\setminus\{\mathbf{0}\}} \lambda_d \sum_{y\in Y_I} \lambda_d^{\rho(y,0)-1}. \end{align}\] By our choice \(\lambda_d=d/(d+1)\), the last expression equals \[\sum_{I\in\mathcal{I}_\gamma}\sum_{y\in Y_I}\lambda_d^{\rho(y,0)-1}=\sum_{y\in Y\setminus\{0\}} \lambda_d^{\rho(y,0)-1}.\] Therefore \[\sum_{x\in X\setminus\{\alpha\}} \lambda_d^{\rho(x,\alpha)-1} \geqslant \sum_{y\in Y\setminus\{0\}} \lambda_d^{\rho(y,0)-1} \geqslant 2.\] The lemma follows. ◻
Definition 8 (\(r\)-separated). Let \(r\in\mathbb{N}\). A subset \(X\subset \mathbb{Z}^d\) is called \(r\)-separated if \(\rho(x,y)\geqslant r\) for all distinct points \(x,y\in X\).
Lemma 8 (Large sets contain admissible ones). Let \(d,n\in\mathbb{N}\). If a subset \(X\subseteq \Lambda_{d,n}\) contains at least \[\binom{n+d-1}{d-1}+1\] points, then \(X\) contains a nonempty admissible subset.
Proof. Let \(S_0\) be the face of the simplex consisting of points whose first coordinate is zero. Then \[|S_0|=\binom{n+d-1}{d-1}.\]
Subsets of \(\Lambda_{d,n}\) form a vector space over \(\mathbb{F}_2\), and admissible sets form a linear subspace. It is enough to show that the codimension of this subspace is at most \(|S_0|\). For this, we show that any lamp configuration can be transformed, by symmetric differences with plaquettes, into a configuration supported inside \(S_0\).
For any point \(x\in\Lambda_{d,n}\setminus S_0\), the plaquette \(x-e_1+T_d\) is contained entirely in \(\Lambda_{d,n}\). Hence one can switch off any lamp with positive first coordinate without affecting lamps whose first coordinate is at least that of \(x\). Proceeding in decreasing order of the first coordinate transforms the configuration into one supported in \(S_0\). ◻
In the metric \(\rho\), balls are infinite. Nevertheless, we can estimate their density, that is, the size of their intersection with a simplex of a given size.
Lemma 9 (Truncated \(\rho\)-ball size). Let \(d,n\in\mathbb{N}\), and let \(0<c<1\). Let \(\alpha\in \Lambda_{d,n}\), and let \[A=\{x\in\mathbb{Z}^d\mid \rho(\alpha,x)\leqslant c\ln n\}.\] Then \[|A\cap \Lambda_{d,n}| \leqslant n^{3c\ln(4d)-c\ln c}.\]
Proof. Put \(r\mathrel{\vcenter{:}}=\lfloor c\ln n\rfloor\). We may assume that \(r\geqslant 1\). Choose \(k\) such that \[2^{k-1}\leqslant n<2^k.\] Consider the quotient map \[\pi\colon\mathbb{Z}^d\to \mathbb{Z}^d/2^k\mathbb{Z}^d.\] Since \(\pi\) does not identify distinct points of \(\Lambda_{d,n}\), it is enough to prove that \[|\pi(A)|\leqslant n^{3c\ln(4d)-c\ln c}.\] Under \(\pi\), all long vectors divisible by \(2^k\) map to zero. Consider the set of \[m\mathrel{\vcenter{:}}=1+kd(d+1)\] vectors \[V\mathrel{\vcenter{:}}=\{0\}\cup\left\{\pm2^a e_i\mid 1\leqslant i\leqslant d,\;0\leqslant a<k\right\} \cup\left\{\pm2^a(e_i-e_j)\mid 1\leqslant i<j\leqslant d,\;0\leqslant a<k\right\}.\] Clearly, \(|\pi(A)|\) is at most the number of ways to choose \(r\) elements from \(V\) with repetitions. Therefore \[\binom{m+r-1}{r} \leqslant\frac{(m+r)^r}{(r/e)^r} \leqslant\left(e\left(1+\frac{m}{r}\right)\right)^{c\ln n} = n^{c\left(1+\ln\left(1+\frac{m}{r}\right)\right)}.\] A routine estimate gives \[1+\frac{m}{r} \leqslant\frac{4(d+1)^2}{c}.\] After a further harmless simplification, we obtain \[|\pi(A)|\leqslant n^{3c\ln(4d)-c\ln c}.\qedhere\] ◻
Corollary 1 (Admissible separated sets exist). For any \(d,n\in\mathbb{N}\) there exists, inside \(\Lambda_{d,n}\), a nonempty admissible \((c_d\ln n)\)-separated set \(X\) of cardinality at most \[\binom{n+d-1}{d-1}+1.\] Here \(c_d>0\) is a constant depending only on \(d\).
Proof. Choose \(c_d>0\) sufficiently small so that, for all relevant \(n\), \[n^{3c_d\ln(4d)-c_d\ln c_d}<\frac{n+d}{d}.\] Equivalently, \[n^{3c_d\ln(4d)-c_d\ln c_d}\binom{n+d-1}{d-1}<|\Lambda_{d,n}|.\]
Put \(r=\lfloor c_d\ln n\rfloor\). We first construct an \(r\)-separated subset of \(\Lambda_{d,n}\) of cardinality \[\binom{n+d-1}{d-1}+1.\] Add points greedily. If no more than \(\binom{n+d-1}{d-1}\) points have already been chosen, Lemma 9 and the inequality above ensure that there is still a point of \(\Lambda_{d,n}\) outside all \(\rho\)-balls of radius \(r\) around the chosen points.
By Lemma 8, this \(r\)-separated set contains a nonempty admissible subset. Any subset of an \(r\)-separated set is again \(r\)-separated. ◻
Recall that chains were introduced in Definition 1, as sequences of configurations varying by a single switch, and that intermediate configurations need not be contained in \(\Lambda_{d,n}\).
Theorem 8 (Configurations with high energy barrier). For each \(d\geqslant 2\), there exist \(\alpha_d>0\) and \(n_d\in\mathbb{N}\) such that, for any \(n\geqslant n_d\), there exists an admissible set \(X\subseteq \Lambda_{d,n}\) with the following property. Any chain \((X_t)_{t=0}^T\) with \(X_0=X\), \(X_T=\varnothing\) satisfies \[|X_t\setminus X|\geqslant n^{\alpha_d}+3|X\setminus X_t|/2\] for some \(i\in\{0,\dots,T\}\). In particular, \(|X_t|\geqslant|X|+n^{\alpha_d}\).
Proof. By Corollary 1, there exists an admissible \((c_d\ln n)\)-separated subset \(X\subset \Lambda_{d,n}\) with \(|X|\leqslant C_d n^{d-1}\). Put \[r=c_d\ln n.\] Let \(|X|=m\), and write \[X=\{x_1,\dots,x_m\}.\] For each \(i\), define \[U_i\mathrel{\vcenter{:}}=\{y\in\mathbb{Z}^d\mid \rho(y,x_i)<r/3\}.\] The sets \(U_i\) are pairwise disjoint.
Let \[X=X_0,X_1,\dots,X_T=\varnothing\] be any sequence of configurations such that the symmetric difference of any two consecutive configurations is a plaquette. We do not require the intermediate configurations to be contained in \(\Lambda_{d,n}\).
Say that a configuration \(X_t\) behaves unusually inside \(U_j\) if \[|X_t\cap U_j|\leqslant 1 \qquad\text{and}\qquad X_t\cap U_j\neq \{x_j\}.\] Initially, the configuration behaves usually inside every \(U_j\), whereas the final configuration behaves unusually inside every \(U_j\). Hence, there is an index \(t\) such that \(X_{t-1}\) behaves usually inside all \(U_j\), while \(X_t\) behaves unusually inside at least one \(U_j\). This index \(j\) is unique, since a single plaquette cannot intersect two distinct sets \(U_j\) and \(U_{j'}\) when \(r\) is large enough. Without loss of generality, assume that \(j=1\).
Put \[X' = X_t \mathbin{\mathpalette\triangle\relax}X.\] Then \(X'\) is admissible, \(x_1\in X'\), and \[|X'\cap U_1|\leqslant 2.\] By Lemma 7, \[\sum_{y\in X'\setminus\{x_1\}} \lambda_d^{\rho(x_1,y)-1}\geqslant 2, \qquad \lambda_d=\frac{d}{d+1}.\] Inside \(U_1\), the set \(X'\setminus\{x_1\}\) contains at most one point, whose contribution to the sum is at most \(1\). Therefore, the points of \(X'\setminus U_1\) contribute at least \(1\) to the sum. Each such point contributes at most \(\lambda_d^{r/3-1}\), and hence \[|X'|\geqslant\lambda_d^{1-r/3} = C'_d\left(1+\frac{1}{d}\right)^{r/3} \geqslant n^{\alpha_d}\] for a suitable constant \(\alpha_d>0\) depending only on \(d\).
Call a point of \(X'\) old if it belongs to \(X\), and new otherwise. Suppose that \(x_j\) is an old point of \(X'\) with \(j\neq1\). Since \(X_t\) behaves usually inside \(U_j\), and \(x_j\notin X_t\), we must have \[|X_t\cap U_j|\geqslant 2.\] Thus, each old point other than possibly \(x_1\) forces at least two new points. Consequently, if \(O\) and \(N\) denote the numbers of old and new points in \(X'\), then \(N\geqslant 2(O-1)\), so \[N\geqslant 5(O-1)/3+N/6=3O/2+|X'|/6-5/3\geqslant 3O/2+n^{\alpha_d}/6-5/3\geqslant 3O/2+n^{\alpha_d/2},\] taking \(n\) large enough in the last inequality. ◻
Let \(A\) and \(B\) be two subsets of \(\mathbb{Z}^d\), which we think of, now, as a space of switches. A subset \(X\subseteq B\) is called an \((A,B)\)-cycle if every translate \(A+x\) that is entirely contained in \(B\) intersects \(X\) in an even number of points. We denote the set of all \((A,B)\)-cycles by \(\mathcal{C}_{(A,B)}\). Using symmetric difference \(\mathbin{\mathpalette\triangle\relax}\) for addition, \(\mathcal{C}_{(A,B)}\) is a vector space over the field \(\mathbb{F}_2\). In fact, set \(V=\mathbb{F}_2^B\) and \(W=\mathbb{F}_2^{\{x\mid A+x\subseteq B\}}\), with \(f\colon V\to W\) given by \(f(v)(x)=\sum_{y\in A+x}v(y)\); then \(\mathcal{C}_{(A,B)}=\ker f\). Here elements of \(V\) are switch configurations, whose image in \(W\) is the corresponding lamp configuration.
Recall the plaquette \(T_d\) from 1 and the size-\(n\) simplex \(\Lambda_{d,n}\) from 6 . We are interested in the set of \((T_d,\Lambda_{d,n})\)-cycles, which we denote by \(\mathcal{C}_{d,n}\). Define the generating polynomial \[\label{eq:def:fdn} G_{d,n}(t)=\sum_{X\in \mathcal{C}_{d,n}} t^{|X|}.\tag{12}\]
Our third main combinatorial result is the following bound:
Theorem 9 (Polynomial estimate). For every \(d\in \mathbb{N}\) there exists \(a_d\in \mathbb{N}\) such that the following holds. For every \(m\in \mathbb{N}\) and every \(s\geqslant a_d\), if \(t>0\) and \(n\in \mathbb{N}\) satisfy \[t \leqslant 1-\frac{1}{m} \qquad\text{and}\qquad n \geqslant m^s,\] then \[G_{d,n}(t) \leqslant 1+n^{-s}.\]
At a high level, the proof of Theorem 9 proceeds as follows (see Figure 3). We will prove a recursive bound (see Proposition 10) of the form \(G_{d,n}(t)\leqslant(G_{d,n/2}(t^{1.5}))^{2^{d-1}}\), provided \(G_{d-1,n}(t^{1/4})\) is sufficiently close to \(1\). To do so, we consider the restriction of a cycle to the even or odd hyperplanes perpendicular to a well chosen direction among \(e_1,\dots,e_d\). We notice that such sections form a collection of \(2^{d-1}\) independent cycles blown up by a factor \(2\) (see Lemma 12). Up to choosing the slicing direction and parity well, we are able to show that (see Lemma 13) the total contribution to the generating polynomial \(G_{d,n}(t)\) of cycles with identical restrictions is smaller than \(t^{1.5|Y|}\), where \(Y\) is the restriction. This relies on the bound on \(G_{d-1,n}(t^{1/4})\) and the fact that symmetric differences of such cycles boil down to \((d-1)\)-dimensional cycles living on a facet of the simplex. Once the recursive bound is established, Theorem 9 is follows by crude estimation carried out in Section 2.3.3. Morally, we just note that iterating the inequality yields \(G_{d,n}(t)\leqslant(1+t^{1.5^{\log_2 n}})^{2^{(d-1)\log_2n}}\to 1\) as \(n\to\infty\).
Lemma 10 (The first layer determines everything). Let \(d,n\in \mathbb{N}\), and let \(1\leqslant k\leqslant d\). Denote by \(S_0\) and \(S_1\) the sets of points in \(\Lambda_{d,n}\) whose \(k\)-th coordinate is equal to \(0\) and \(1\) respectively.
If \(X\in \mathcal{C}_{d,n}\) and \(X\cap S_0=\emptyset\), then \(X=\emptyset\).
If \(X\in \mathcal{C}_{d,n}\) and \(X\cap S_1=\emptyset\), then \(X\subseteq S_0\). Moreover, \(X\in \mathcal{C}_{(T',S_0)}\), where \(T'=T_d\setminus\{e_k\}\).
Proof. [item:cyc95layer:1]. Suppose that \(|X|>0\). Choose a point \(x\in X\) with minimal \(k\)-th coordinate. Since \(x\notin S_0\), the plaquette \(x-e_k+T_d\) is entirely contained in \(\Lambda_{d,n}\), and it contains exactly one point of \(X\), namely \(x\). This contradicts the cycle condition.
[item:cyc95layer:2]. Let \(S_+\mathrel{\vcenter{:}}=\Lambda_{d,n}\setminus S_0\), that is, the set of all points whose \(k\)-th coordinate is positive. Clearly, \(X\cap S_+\) is a \((T_d,S_+)\)-cycle. By part 1, already proved above, we have \(X\cap S_+=\emptyset\), hence \(X\subseteq S_0\).
For any translate \(T'+x\) contained entirely in \(S_0\), there is a point of \(S_1\) which completes it to a translate of \(T_d\). Since this point does not belong to \(X\), it follows that \(|(T'+x)\cap X|\) is even. Therefore \(X\in \mathcal{C}_{(T',S_0)}\). ◻
We now show that every cycle is either empty or has large cardinality.
Lemma 11 (Cycles are macroscopic). Let \(X\in \mathcal{C}_{d,n}\) and assume \(X\neq \emptyset\). Then \(|X|\geqslant n+1\).
Proof. We prove this by double induction: first on \(n\), and then on \(d\). The base case is \(n=0\) and arbitrary \(d\).
Let \(X\in \mathcal{C}_{d,n}\). Decompose \(\Lambda_{d,n}\) into two parts: \[\Lambda_{d,n}=S_0\sqcup S_+,\] where \(S_0\) consists of the points whose \(d\)-th coordinate is equal to \(0\), and \(S_+\) consists of all remaining points.
Observe that \(S_+\) is a translate of \(\Lambda_{d,n-1}\), while \(S_0\) is an embedding of \(\Lambda_{d-1,n}\) into \(\mathbb{Z}^d\). We assume that the induction hypothesis applies to these sets.
Case 1. Suppose that \(X\cap S_0=\emptyset\). This is impossible by part 1 of Lemma 10.
Case 2. Suppose that \(X\cap S_+=\emptyset\). Then by part 2 of Lemma 10, the set \(X\) is a \((T',S_0)\)-cycle. Applying the induction hypothesis for \(d-1\) and \(n\), we obtain \(|X|\geqslant n+1\).
Case 3. Suppose that both \(X\cap S_0\) and \(X\cap S_+\) are nonempty. Then \(X\cap S_+\) is a \((T_d,S_+)\)-cycle, and by the induction hypothesis we have \[|X\cap S_+|\geqslant n.\] Hence \(|X|\geqslant n+1\). ◻
For \(d=1\), the generating polynomial of cycles can be computed explicitly: the set \(\Lambda_{1,n}\) is an interval of length \(n\). In this case, \[\mathcal{C}_{1,n}=\{\Lambda_{1,n},\emptyset\},\] and therefore \[G_{1,n}(t)=1+t^{n+1}.\]
Proposition 10 (Recursive relation). Let \(d>1\). Assume that for some \(t\in (0,1)\) and \(n\in \mathbb{N}\) one has \[G_{d-1,n}\left(t^{1/4}\right)\leqslant 1+\frac{1}{4d}.\] Then \[G_{d,n}(t)\leqslant \left(G_{d,\left\lfloor \frac{n-\epsilon}{2} \right\rfloor}\left(t^{1.5}\right)\right)^{2^{d-1}}\] for some \(0\leqslant \epsilon\leqslant d\).
Proof. First, let us derive a simple consequence of the condition imposed on \(t\) and \(n\).
The simplex \(\Lambda_{d-1,n}\) contains the empty cycle. It also contains a cycle of size \(n+1\), namely the set of all points whose coordinates other than the first are equal to \(0\). Hence \[G_{d-1,n}\left(t^{1/4}\right)\geqslant 1+t^{(n+1)/4},\] and therefore \[t^{(n+1)/4}\leqslant \frac{1}{4d}. \label{ineq:t}\tag{13}\]
Now we examine the geometry of the simplex \(\Lambda_{d,n}\) more closely. Partition all points of \(\Lambda_{d,n}\) into \(2^d\) subsets according to the parities of their coordinates (see Figure 3). For \(r=(r_1,\dots,r_d)\in \{0,1\}^d\), define \[\label{eq:def:Lambda:dnr} \Lambda_{d,n}^r\mathrel{\vcenter{:}}=\left\{(x_1,\dots,x_d)\in \Lambda_{d,n}\mid x_i\equiv r_i \pmod 2 \text{ for all } 1\leqslant i\leqslant d\right\}.\tag{14}\]
Then, clearly, \[\Lambda_{d,n}=\bigsqcup_{r\in \{0,1\}^d}\Lambda_{d,n}^r.\]
For each \(r\in \{0,1\}^d\), the set \(\Lambda_{d,n}^r\) looks like a simplex scaled by a factor of \(2\), with side length approximately \(n/2\). More precisely, \[\Lambda_{d,n}^r=r+2\Lambda_{d,n_r}, \label{eq:rec}\tag{15}\] where \[n_r\mathrel{\vcenter{:}}=\left\lfloor \frac{n-(r_1+\dots+r_d)}{2}\right\rfloor.\]
Formula 15 defines a bijection \[\tau_r\colon\Lambda_{d,n}^r\to\Lambda_{d,n_r}.\]
If \(X\subseteq \Lambda_{d,n}\), define \[\tau_r(X)\mathrel{\vcenter{:}}=\tau_r\left(X\cap\Lambda_{d,n}^r\right).\]
Thus, we may view a subset \(X\subseteq \Lambda_{d,n}\) as a collection of \(2^d\) subsets of the form \(\tau_r(X)\).
Lemma 12 (Restrictions of cycles are renormalised cycles). If \(X\in \mathcal{C}_{d,n}\), then \(\tau_r(X)\in \mathcal{C}_{d,n_r}\) for all \(r\in \{0,1\}^d\).
Proof. The statement follows immediately from the identity \[T_d \mathbin{\mathpalette\triangle\relax}(e_1+T_d)\mathbin{\mathpalette\triangle\relax}\dots \mathbin{\mathpalette\triangle\relax}(e_d+T_d)=2T_d.\]
This identity is easy to verify directly: all points of the form \(e_i+e_j\), where \(i\neq j\), cancel out. It can also be interpreted as the Frobenius identity \[(1+x_1+\dots+x_d)^2=1+x_1^2+\dots+x_d^2\] over a field of characteristic \(2\). (This observation is not new; it is used, for instance, to construct a renormalisation for the Ledrappier subshift.) ◻
Consider \(1\leqslant k\leqslant d\) and \(\varepsilon\in \{0,1\}\), and define \[\label{eq:def:Lambda:dnke} \Lambda_{d,n}^{k,\varepsilon}\mathrel{\vcenter{:}}=\left\{(x_1,\dots,x_d)\in \Lambda_{d,n}\mid x_k\equiv \varepsilon \pmod 2\right\}\tag{16}\] (see Figure 3). We call subsets of \(\Lambda_{d,n}\) of this form binary layers. There are \(2d\) binary layers in total. Binary layers with the same \(k\) and different values of \(\varepsilon\) will be called complementary; the one with \(\varepsilon=0\) will be called large, and the one with \(\varepsilon=1\) small.
If in addition \(X\in \mathcal{C}_{d,n}\), then we call \(X\cap\Lambda_{d,n}^{k,\varepsilon}\) a partial cycle of type \((k,\varepsilon)\). The set of all such partial cycles is denoted by \(\mathcal{C}_{d,n}^{k,\varepsilon}\).
Lemma 13 (Cycles with fixed even or odd sections). Let \(1\leqslant k\leqslant d\), \(\varepsilon\in \{0,1\}\), and suppose that \(t\) and \(n\) satisfy the assumptions of Proposition 10. Let \(Y\) be a nonempty partial cycle of type \((k,\varepsilon)\), and let \(X_1,\dots,X_m\) be distinct cycles from \(\mathcal{C}_{d,n}\) such that \[X_i\cap\Lambda_{d,n}^{k,\varepsilon}=Y \quad\text{and}\quad |X_i|\geqslant 2|Y|\] for every \(i\). Then \[\sum_{i=1}^{m} t^{|X_i|}\leqslant \frac{t^{1.5|Y|}}{2d}.\]
Proof. Let \(S_0\) be the set of points in \(\Lambda_{d,n}\) whose \(k\)-th coordinate is equal to \(0\). We distinguish two cases depending on the value of \(\varepsilon\).
Case 1: \(\varepsilon=0\). We claim that in this case \(m\leqslant 1\). Indeed, suppose that \(X_1\cap\Lambda_{d,n}^{k,0}=X_2\cap\Lambda_{d,n}^{k,0}\). Consider the symmetric difference \[X'\mathrel{\vcenter{:}}= X_1\mathbin{\mathpalette\triangle\relax}X_2.\] Then \(X'\cap S_0=\emptyset\), hence \(X'=\emptyset\) by part 1 of Lemma 10. Therefore \(X_1=X_2\), so it remains to check that \[t^{|X_1|}\leqslant \frac{t^{1.5|Y|}}{2d}.\] We know that \(|X_1|\geqslant n+1\) by Lemma 11, and that \(|Y|\leqslant |X_1|/2\), hence, by 13 , \[t^{|X_1|-1.5|Y|} \leqslant t^{|X_1|/4} \leqslant t^{(n+1)/4} \leqslant \frac{1}{4d}.\]
Case 2: \(\varepsilon=1\). Without loss of generality, we may assume that \(|X_i|\geqslant |X_1|\) for all \(1<i\leqslant m\). As in the previous case, \[t^{|X_1|-1.5|Y|}\leqslant t^{(n+1)/4},\] so we only need to estimate the sum of the remaining terms.
Define \[L_i\mathrel{\vcenter{:}}= X_i\mathbin{\mathpalette\triangle\relax}X_1 \qquad (1\leqslant i\leqslant m),\] so in particular \(L_1=\emptyset\). Clearly, \(L_i\cap \Lambda_{d,n}^{k,1}=\emptyset\) for every \(i\), so part 2 of Lemma 10 implies that \[L_i\subseteq S_0 \quad\text{and}\quad L_i\in \mathcal{C}_{(T',S_0)},\] where \(T'=T_d\setminus\{e_k\}\). Thus each \(L_i\) can be viewed simply as a cycle in a simplex of smaller dimension. Therefore, for every \(\lambda\in (0,1)\), \[\sum_{i=2}^m \lambda^{|L_i|} \leqslant \sum_{\substack{X\in \mathcal{C}_{d-1,n}\\ X\neq \emptyset}} \lambda^{|X|} = G_{d-1,n}(\lambda)-1. \label{ineq:Li}\tag{17}\]
For \(1\leqslant i\leqslant m\), let \[Z_i\mathrel{\vcenter{:}}= X_i\cap\Lambda_{d,n}^{k,0},\] so \(Z_i=X_i\setminus Y\). For every \(i\) we have \(Z_i=Z_1\mathbin{\mathpalette\triangle\relax}L_i\), and we also have \[|Z_i|\geqslant |Z_1| \quad\text{and}\quad |Z_i|\geqslant |Y|.\]
We claim that \[|Z_i|\geqslant |L_i|/2 \qquad\text{for all } i.\] Indeed, if \(|L_i|\geqslant 2|Z_1|\), then \(|Z_i|\geqslant |L_i|-|Z_1|\geqslant |L_i|/2\), and if \(|L_i|\leqslant 2|Z_1|\) then \(|Z_i|\geqslant |Z_1|\geqslant |L_i|/2\). Combining this with \(|Z_i|\geqslant |Y|\), we have \[|Z_i|-|Y|/2 \geqslant |Z_i|/2 \geqslant |L_i|/4.\]
Using these inequalities, we obtain \[\sum_{i=2}^{m} t^{|X_i|-1.5|Y|} = \sum_{i=2}^{m} t^{|Z_i|-|Y|/2} \leqslant \sum_{i=2}^{m} t^{|L_i|/4} \leqslant G_{d-1,n}\left(t^{1/4}\right)-1.\]
Therefore, by 13 , \[\sum_{i=1}^{m} t^{|X_i|-1.5|Y|} \leqslant t^{(n+1)/4}+G_{d-1,n}\left(t^{1/4}\right)-1 \leqslant \frac{1}{2d}.\qedhere\] ◻
If \(1\leqslant k\leqslant d\) and \(\varepsilon\in \{0,1\}\), let \(R_{k,\varepsilon}\) denote the set of all binary strings of length \(d\) whose \(k\)-th coordinate is equal to \(\varepsilon\).
We will prove the following inequality: \[G_{d,n}(t)\leqslant \frac{1}{2d}\sum_{\substack{1\leqslant k\leqslant d\\ \varepsilon\in \{0,1\}}} \prod_{r\in R_{k,\varepsilon}} G_{d,n_r}\left(t^{1.5}\right). \label{eq}\tag{18}\]
Observe that this inequality immediately implies Proposition 10.
Lemma 14 (Slicing direction). Let \(X\in \mathcal{C}_{d,n}\) and assume that \(|X|>0\). Then there exists \(k\) such that \[X\cap\Lambda_{d,n}^{k,0}\neq\emptyset\text{ and }X\cap\Lambda_{d,n}^{k,1}\neq\emptyset.\]
Proof. Let \(x\in X\) and \(y\in \Lambda_{d,n-1}\) be such that \(x\in y+T_d\subseteq\Lambda_{d,n}\). Since \(X\in\mathcal{C}_{d,n}\), \(|(y+T_d)\cap X|\geqslant 2\), so there exists \(z\in(y+T_d)\setminus\{x\}\). Then any \(k\) such that \(\langle x-z,e_k\rangle\neq 0\) is as desired. ◻
Let \(X\in \mathcal{C}_{d,n}\) be a nonzero cycle. Choose, according to Lemma 14, values of \(k\) and \(\varepsilon\) such that \[0<\left|X\cap\Lambda_{d,n}^{k,\varepsilon}\right|\leqslant |X|/2.\] Let \[Y=X\cap\Lambda_{d,n}^{k,\varepsilon}.\] We say that the cycle \(X\) is of type \((k,\varepsilon,Y)\).
In this way, we define a type for every nonzero cycle.
We now group the nonconstant terms on the left-hand side of 18 by type, and estimate the sum within each type using Lemma 13: \[G_{d,n}(t) = \sum_{X\in \mathcal{C}_{d,n}} t^{|X|} = 1+\sum_{X \text{ is of type }(k,\varepsilon,Y)} t^{|X|} \leqslant 1+\frac{1}{2d}\sum_{\substack{1\leqslant k\leqslant d\\ \varepsilon\in \{0,1\}}} \sum_{\substack{Y\in \mathcal{C}_{d,n}^{k,\varepsilon}\\ Y\neq \emptyset}} t^{1.5|Y|}.\] Equivalently, \[G_{d,n}(t) \leqslant \frac{1}{2d}\sum_{\substack{1\leqslant k\leqslant d\\ \varepsilon\in \{0,1\}}} \sum_{Y\in \mathcal{C}_{d,n}^{k,\varepsilon}} t^{1.5|Y|}.\]
Now fix \(k\) and \(\varepsilon\). Consider a partial cycle \(Y\in \mathcal{C}_{d,n}^{k,\varepsilon}\). The support of \(Y\) lies inside \[\Lambda_{d,n}^{k,\varepsilon}=\bigsqcup_{r\in R_{k,\varepsilon}} \Lambda_{d,n}^r,\] and \[t^{1.5|Y|} = \prod_{r\in R_{k,\varepsilon}} \left(t^{1.5}\right)^{|\tau_r(Y)|}.\]
By Lemma 12, all sets \(\tau_r(Y)\) are cycles in the corresponding simplices. Therefore the monomial \(t^{1.5|Y|}\) appears in the expansion of \[\prod_{r\in R_{k,\varepsilon}} G_{d,n_r}\left(t^{1.5}\right) = \prod_{r\in R_{k,\varepsilon}} \sum_{Z\in \mathcal{C}_{d,n_r}} \left(t^{1.5}\right)^{|Z|},\] namely by choosing in each factor the term corresponding to \(\tau_r(Y)\).
Different \(Y\) correspond to different such choices. Therefore, for positive \(t\), we obtain \[\sum_{Y\in \mathcal{C}_{d,n}^{k,\varepsilon}} t^{1.5|Y|} \leqslant \prod_{r\in R_{k,\varepsilon}} G_{d,n_r}\left(t^{1.5}\right),\] which completes the proof of 18 , and hence also of Proposition 10. ◻
Lemma 15 (Crude bound). Let \(d,n\in \mathbb{N}\). Then the function \(G_{d,n}(t)\) is increasing on the interval \([0,1]\), and on this interval it satisfies \[G_{d,n}(t) \leqslant 1 + 2^{(n+1)^{d-1}} t^{n+1}.\]
Proof. Monotonicity is obvious. By Lemma 11, every nontrivial cycle has size at least \(n+1\), hence for \(0\leqslant t\leqslant 1\), \[G_{d,n}(t) \leqslant 1 + |\mathcal{C}_{d,n}|\, t^{n+1}.\]
It remains to estimate the total number of cycles. By part 1 of Lemma 10, a cycle is uniquely determined by its intersection with \(S_0\), and \[|S_0|\leqslant (n+1)^{d-1}.\] Therefore \[|\mathcal{C}_{d,n}| \leqslant 2^{(n+1)^{d-1}}.\qedhere\] ◻
We are now ready to prove Theorem 9 by induction.
Proof of Theorem 9. Base case: \(d=1\). We know that \[G_{1,n}(t)=1+t^{n+1}.\] Take \(a_1=6\). We must check that for \(s\geqslant 6\) and \(n\geqslant m^s\) one has \[\left(1-\frac{1}{m}\right)^{n+1} \leqslant n^{-s}.\] Using \((1-1/m)^m\leqslant e^{-1}\), it is enough to prove \[\frac{n+1}{m}\geqslant s\ln n.\] Since \(m\leqslant n^{1/s}\), it is enough to show that for \(n\geqslant m^s\geqslant 2^s\) we have \[n+1\geqslant s n^{1/s}\ln n;\] this holds for \(n=2^s\) since \(s\geqslant 6\), and also for larger \(n\) by convexity.
Induction step.
Fix \(a_{d-1}>d\) such that for all \(m\geqslant 2\), \(s\geqslant a_{d-1}\), and \(n>m^s\), one has \[G_{d-1,n}\!\left(1-\frac{1}{m}\right) < 1+n^{-s}.\]
Then, whenever \(t<1-\frac{1}{m}\) and \(n\geqslant m^{3a_{d-1}}\), we automatically have \[n \geqslant m^{3a_{d-1}} \geqslant (4m)^{a_{d-1}} \geqslant 8^d > 4d,\] and therefore \[G_{d-1,n}\left(\sqrt[4]{t}\right) \leqslant G_{d-1,n}\!\left(1-\frac{1}{4m}\right) < 1+n^{-a_{d-1}} < 1+\frac{1}{n} < 1+\frac{1}{4d}.\]
Hence the assumption of Proposition 10 is satisfied, so \[G_{d,n}(t) \leqslant \left( G_{d,\left\lfloor \frac{n-\epsilon}{2}\right\rfloor}\left(t^{1.5}\right) \right)^{2^{d-1}}\] for some \(0\leqslant \epsilon \leqslant d\). Since \(n>4d\), we have \[\left\lfloor \frac{n-\epsilon}{2}\right\rfloor \geqslant \frac{n}{4}.\]
Choose \(a_d\) large enough so that the following conditions hold: \[\begin{gather} \tag{19}a_d > 6a_{d-1};\\ \tag{20}x^{a_d\ln(1.5)/4} \geqslant 3x \cdot \ln 2 \cdot x^{3(d-1)a_{d-1}}\text{ for all }x\geqslant 2;\\ \tag{21}a_d > \frac{12}{\ln 1.5};\\ \tag{22}1.5^{x/24} \geqslant x\text{ for all }x\geqslant a_d;\\ \tag{23}x^{\ln1.5}/12 > 3\ln x\text{ for all }x\geqslant 2^{a_d};\\ \tag{24}1.5^{\frac{2}{3}x} \geqslant 3\ln 2\cdot dx\text{ for all }x\geqslant \dfrac{a_d}{8}. \end{gather}\]
These conditions will be used below. It is easy to check that each of them holds for all sufficiently large \(a_d\).
Assume now that \[t<1-\frac{1}{m},\qquad s\geqslant a_d,\qquad n\geqslant m^s.\]
We apply Proposition 10 repeatedly as long as possible. At each step, the parameter \(n\) decreases, and we may continue until \(n\) becomes smaller than \(m^{3a_{d-1}}\). The parameter \(t\) also decreases, but this does not interfere with the application of the lemma: \[\label{eq:fdn:recurrence} G_{d,n}(t) \leqslant \left(G_{d,n_1}\left(t^{1.5}\right)\right)^{2^{d-1}} \leqslant \left(G_{d,n_2}\left(t^{1.5^2}\right)\right)^{(2^{d-1})^2} \leqslant \dots \leqslant \left(G_{d,n_r}\left(t^{1.5^r}\right)\right)^{(2^{d-1})^r},\tag{25}\] where \((n_i)\) is a decreasing sequence of positive integers, with \[2 \leqslant \frac{n_i}{n_{i+1}} \leqslant 4, \qquad\text{and}\qquad n_r < m^{3a_{d-1}}.\]
This gives a lower bound on \(r\): \[r \geqslant \frac{\ln n - \ln m^{3a_{d-1}}}{\ln 4}.\]
By 19 , we have \(\ln m^{3a_{d-1}} < 0.5\ln n\), so \[r \geqslant \frac{\ln n}{4} \geqslant \frac{s\ln m}{4} \geqslant \frac{s}{8}. \label{neq:r}\tag{26}\]
By Lemma 15, for any \(0<\lambda<1\) we have \[G_{d,n_r}(\lambda) \leqslant 1 + 2^{(n_r+1)^{d-1}} \lambda^{n_r+1}\leqslant 1 + 2^{(m^{3a_{d-1}})^{d-1}}{\lambda},\] so 25 gives \[G_{d,n}(t) \leqslant \left( 1 + 2^{m^{3(d-1)a_{d-1}}} \, t^{1.5^r} \right)^{(2^{d-1})^r}.\]
We want to check that this is at most \(1+n^{-s}\).
If \(k\) is a positive integer and \(0 \leqslant 2kx<1\), then \[(1+x)^k \leqslant 1+2kx,\] which is easily proved by induction. So it is enough to verify \[2\cdot2^{m^{3(d-1)a_{d-1}}}\cdot t^{1.5^r}\cdot \left(2^{d-1}\right)^r < n^{-s}.\]
Taking logarithms, this becomes \[\ln 2 \cdot m^{3(d-1)a_{d-1}} + (\ln t)\cdot 1.5^r + \ln 2 \cdot (1+r(d-1)) < -s\ln n.\]
Since \(t<1-\frac{1}{m}\), we have \(\ln t < -\frac{1}{m}\). We also have \(r\geqslant \frac{s}{8} \geqslant 1\). Thus, it is enough to prove \[\label{eq:poly} \ln 2 \cdot m^{3(d-1)a_{d-1}} + s\ln n + \ln 2 \cdot dr < \frac{1.5^r}{m}.\tag{27}\]
To do this, we show that each of the three terms on the left is at most \(\frac{1.5^r}{3m}\).
By 26 and 20 , \[\label{eq:poly:1} 1.5^r \geqslant m^{s\ln1.5/4} \geqslant m^{a_d\ln1.5/4} \geqslant 3m\cdot \ln 2 \cdot m^{3(d-1)a_{d-1}}.\tag{28}\]
By 26 and 21 , \[\label{eq:2:a} \frac{r}{3} \geqslant \frac{s\ln m}{12}\geqslant\frac{a_d\ln m}{12}\geqslant\frac{\ln m}{\ln 1.5}.\tag{29}\] Again by 26 and 22 , we have \[\label{eq:2:b} 1.5^{r/3} \geqslant 1.5^{s/24} \geqslant s.\tag{30}\] By 26 and 23 , we also have \[\label{eq:2:c} \frac{r}{3} \geqslant \frac{\ln n}{12}\geqslant\frac{\ln(3\ln n)}{\ln 1.5}.\tag{31}\] Putting 29 , 30 and 31 together, we get \[\label{eq:poly:2} 3ms\ln n \leqslant 1.5^r.\tag{32}\]
By 29 , 26 and 24 , we have \[\label{eq:poly:3} 1.5^r \geqslant m\cdot1.5^{2r/3}\geqslant 3m\ln 2\cdot dr.\tag{33}\]
In this section we prove the lower bounds of Theorems 1 and 2. Given Theorem 4, the route we follow is relatively standard in the context of kinetically constrained models. The two proofs are very similar, but we start with the simpler one.
Proof of the lower bound of Theorem 2. The lower bound holds for general even \(d\geqslant 2\). For \(\varepsilon>0\) small enough (depending on \(d\)), consider \(\beta>1/\varepsilon^2\) and an integer \(k\) such that \(k\in(1/\varepsilon,\beta\varepsilon)\) (if it exists). We consider the torus \({\mathbb{T}} _n\) with \(n=2^k\). Set \(\Omega'={\mathbb{F}} _2^{{\mathbb{T}} _n}\), and define the mapping \[\label{eq:def:Phi} \Phi\colon\Omega_{{\mathbb{T}} _n}\to\Omega'\colon\sigma\mapsto\left(\frac{1-[\sigma]_{x+T_d}}{2}\right)_{x\in{\mathbb{T}} _n}.\tag{34}\] from the spin configuration to the lamp configuration. Then, recalling the definitions of \(T_d\) from 1 and \(H\) from 2 , we have \[H_{{\mathbb{T}} _n}(\sigma)=-|{\mathbb{T}} _n|+2\sum_{x\in{\mathbb{T}} _n}\Phi(\sigma)_x,\] while \(\Phi(\sigma^x)=\Phi(\sigma)+\mathbb{I}_{x-T_d}\); here and below addition takes place in \({\mathbb{F}} _2^{{\mathbb{T}} _n}\), and for a set or element \(S\) we denote by \(\mathbb{I}_S\) the characteristic function of \(S\). We then see that the image of the Glauber dynamics on \({\mathbb{T}} _n\) via \(\Phi\) is the Markov process with generator \[\label{eq:def:L:prime}\mathcal{L}'f(\omega)=\sum_{x\in{\mathbb{T}} _n}\frac{e^{-2\beta\sum_{y\in x-T_d}(1-\omega_y)}}{e^{-2\beta\sum_{y\in x-T_d}\omega_y}+e^{-2\beta\sum_{y\in x-T_d}(1-\omega_y)}}\left(f(\omega+\mathbb{I}_{x-T_d})-f(\omega)\right).\tag{35}\] We denote by \(\pi\) the product Bernoulli measure on \(\Omega'\) with parameter \[\label{eq:def:p} p=1/(1+e^{2\beta})\in\left[1/(2e^{2\beta}),e^{-2\beta}\right],\tag{36}\] which is clearly invariant for this dynamics.
Consider a central box \(B=\{-L,\dots,L\}^d\) with \(L=2^{k-2}\). Define the events \[\begin{align} \mathcal{E}_0&{}=\left\{\omega\in\Omega'\mid\forall x\in B,\omega_x=0\right\},& \mathcal{E}_1&{}=\left\{\omega\in\Omega'\mid\omega_B=\mathbb{I}_{O}\right\}. \end{align}\] Recalling Definition 1 and Theorem 4, we say that \(\omega\in\Omega'\) is good if there exists a chain \((\omega^{(i)})_{i=0}^m\) such that \(\omega^{(0)}=\omega\), \(\omega^{(m)}_B=\mathbb{I}_{O}\) and \(|\omega^{(i)}|\mathrel{\vcenter{:}}=\sum_{x\in B}\omega^{(i)}_x\leqslant k-d-4\) for all \(i\in\{0,\dots,m\}\). Let \(\mathcal{G}\) be the set of good configurations. By Theorem 4 applied to the chain \((\omega^{(m-i)}+\omega^{(m)}+\mathbb{I}_O)_{i=0}^m\), we obtain \(\mathcal{G}\cap\mathcal{E}_0=\varnothing\). Indeed, we may apply Theorem 4 to \({\mathbb{T}} _m\) rather than to \({\mathbb{Z}} ^d\) when \(m\geqslant 4+3d2^{k+1}\): consider a counterexample chain on the torus of minimal length; then all its configurations are supported in \(B_k+\{-1,0,1\}^d\), so it lifts to a chain on \({\mathbb{Z}} ^d\) that contradicts Theorem 4. Consequently, \[\mathop{\mathrm{Var}}'(\mathbb{I}_{\mathcal{G}})=\pi(\mathcal{G})(1-\pi(\mathcal{G}))\geqslant\pi(\mathcal{E}_1)\pi(\mathcal{E}_0)=p(1-p)^{2|B|-1}.\] Moreover, \[\begin{align} \mathcal{D}'\left(\mathbb{I}_{\mathcal{G}}\right)&{}=-\pi\left(\mathbb{I}_\mathcal{G}\mathcal{L}'\mathbb{I}_\mathcal{G}\right)\leqslant 2\sum_{x\in{\mathbb{T}} _n}\pi\left(\omega\in\mathcal{G}\not\ni\omega+\mathbb{I}_{x-T_d}\right)\\ &{}\leqslant 2|B+T_d|\pi\left(\sum_{x\in B}\omega_x>k-d-4-|T_d|\right) \leqslant 2|B+T_d|(p|B|)^{k-2d-4}. \end{align}\] Putting these inequalities together, we obtain \[\begin{align} \nonumber\frac{\mathop{\mathrm{Var}}'(\mathbb{I}_\mathcal{G})}{\mathcal{D}'(\mathbb{I}_\mathcal{G})}&{}\geqslant\frac{p(1-p)^{2|B|-1}}{2|B+T_d|(p|B|)^{k-2d-4}}\geqslant p^2\left(p2^{dk}\right)^{2d+4-k}\\ &{}\geqslant e^{(2-\varepsilon d\ln 2)\beta(k-2d-4)-4\beta}/4\geqslant e^{\beta k(2-\varepsilon d\ln 2)(1-\varepsilon(2d+9))}\geqslant e^{2\beta k (1-O(\varepsilon))}.\label{eq:trel:lower39} \end{align}\tag{37}\]
Lemma 16. Whenever \(n\) is a power of \(2\) and \(d\) is even, the map \(\Phi\) from 34 is bijective.
Proof. Note that \(|\Omega_{{\mathbb{T}} _n}|=|\Omega'|\), that \(\Phi\) is a group homomorphism, and that \(\Omega'_{{\mathbb{T}} _n}\) is a vector space with basis \((\mathbb{I}_x)_{x\in{\mathbb{T}} _n}\). Therefore, it is enough to show that, for any \(x\in{\mathbb{T}} _n\), there exists \(\sigma\in\Omega_{{\mathbb{T}} _n}\) such that \(\Phi(\sigma)=\mathbb{I}_{x}\). By symmetry we only show this for \(x=O\). Indeed, the Sierpinski gasket configuration (see Figure 4) \[\label{eq:Sierpinski} \sigma_y=(-1)^{\binom{-\sum_{i=1}^dy_i}{-y_1,\dots,-y_d}}\tag{38}\] for \(-y\in\{0,\dots,n-1\}^d\) verifies this property for \(d\) even and \(n\) power of \(2\). ◻
To prepare for the infinite volume case, we need the following substitute for Lemma 16 which fails in infinite volume.
Lemma 17. Consider \(\Omega'={\mathbb{F}} _2^{{\mathbb{Z}} ^d}\) and \[\Phi\colon\Omega\to\Omega'\colon\sigma\mapsto\left(\frac{1-[\sigma]_{x+T_d}}{2}\right)_{x\in{\mathbb{Z}} ^d}.\] Let \(\Phi(\mu)\) denote the push forward of \(\mu\) under \(\Phi\) and \(\pi\) denote the Bernoulli product measure on \({\mathbb{Z}} ^d\) with parameter \(p=1/(1+e^{2\beta})\). Then \(\Phi(\mu)=\pi\).
Proof. Fix \(\varepsilon>0\) and finite \(\Lambda\subset{\mathbb{Z}} ^d\). We will prove (independently) in Section 3.4 that the Gibbs measure \(\mu\) is unique. Consequently, we can find \(\Lambda'\supset\Lambda+\{-1,0,1\}^d=:\Lambda_+\) large enough so that \[\sup_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}}d_{\mathrm{TV}}\left(\mu^\tau_{\Lambda_+},\mu_{\Lambda_+}\right)<\varepsilon\] (\(d_{\mathrm{TV}}\) is the total variation distance, that is, half of the \(\ell^1\)-distance), where \(\mu^\tau_{\Lambda_+}\) is the marginal of \(\mu^\tau_{\Lambda'}\) on \(\Lambda_+\) and similarly for \(\mu\) instead of \(\mu^\tau_{\Lambda'}\). Note that \(\Omega\to\Omega'_\Lambda:\sigma\mapsto\Phi(\sigma)_{\Lambda}\) is local with support contained in \(\Lambda_+\) and surjective. Indeed, we may change the value at a single vertex in \(\Lambda\) by adding to \(\sigma\) a large enough Sierpiński gasket translated at that vertex (see 38 and Figure 4).
Consider the image under \(\Phi\) of the chain with generator \(\mathcal{L}_{\Lambda'}^\tau\). As above, it is a Markov chain, with invariant measure \(\pi^\tau_{\Lambda'}=\pi_{\Lambda'-T_d}(\cdot\mid\Omega'_{\Lambda',\tau})\), where \(\Omega'_{\Lambda',\tau}\subseteq\Omega'_{\Lambda'}\mathrel{\vcenter{:}}={\mathbb{F}} _2^{\Lambda'-T_d}\) is the irreducible component of \(\mathcal{L}\) induced by the boundary condition \(\tau\). In particular, \(\Omega'_{\Lambda',\tau}\) are either equal or disjoint for different \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}\). By the surjectivity above, \(\bigcup_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}}\Omega'_{\Lambda',\tau}=\Omega'_{\Lambda'}\), so we can find a probability measure \(\mathrm d\theta(\tau)\) on \(\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}\) such that \(\theta(\pi_{\Lambda'}^\tau)=\pi_{\Lambda'}\). Since, for any \(\tau,\tau'\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}\), \[d_{\mathrm{TV}}(\pi^\tau_{\Lambda},\pi^{\tau'}_{\Lambda})\leqslant d_{\mathrm{TV}}(\mu_{\Lambda+}^\tau,\mu_{\Lambda_+}^{\tau'})\leqslant 2\varepsilon,\] we deduce \[d_{\mathrm{TV}}\left(\pi_{\Lambda},\Phi(\mu)_\Lambda\right)\leqslant 2\varepsilon+\sup_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}}d_{\mathrm{TV}}\left(\pi_\Lambda^\tau,\Phi(\mu)_\Lambda\right)\leqslant 2\varepsilon+\sup_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}}d_{\mathrm{TV}}\left(\mu^\tau_{\Lambda_+},\mu_{\Lambda_+}\right)\leqslant 3\varepsilon.\] Therefore, \(\pi=\Phi(\mu)\). ◻
Proof of the lower bound of Theorem 1. The rest of the proof proceeds exactly as the lower bound of Theorem 2, so we only indicate the non-trivial changes.
Fix \(d\geqslant 2\) and \(\beta>0\). Let \(k=\lceil\beta/(d\ln 2)\rceil\), \(L=3d2^k\), \(B=\{-L,\dots,L\}^d\), so that \[e^{\beta}\leqslant L^d\leqslant|B|\leqslant(3L)^d\leqslant(18d)^de^{\beta}.\] Define \(\mathcal{G}\subset\Omega'\) as in the proof of Theorem 2 with chains featuring at most \(k\) lamps ON (instead of \(k-d-4\)). Then, we similarly get \[\label{eq:trel:lower} \frac{\mathop{\mathrm{Var}}'(\mathbb{I}_\mathcal{G})}{\mathcal{D}'(\mathbb{I}_\mathcal{G})}\geqslant\frac{p(1-p)^{2|B|-1}}{2|B+T_d|(p|B|)^{k-d}}\geqslant c e^{c\beta^2}\tag{39}\] for a small enough \(c=c(d)>0\). By Lemma 17 and 39 , plugging the local function \(\mathbb{I}_{\Phi^{-1}(\mathcal{G})}\) into 5 gives the desired lower bound on \(T_{\mathrm{rel}}\). ◻
In this section, we prove the lower bound of Theorem 3, using Theorem 8.
Proof of the lower bound of Theorem 3. Fix \(d\geqslant 2\), \(\alpha_d\leqslant 1\) and \(n_d\) as in Theorem 8. Let \(\beta>0\) and \(n\geqslant n_d\) satisfy \(n\leqslant e^{2\beta/(5d)}\). Let \(\Lambda=-\Lambda_{d,n}\). Define the mapping \[\Phi\colon\Omega_{\Lambda}\to{\mathbb{F}} _2^{{\mathbb{Z}} ^d}\colon\sigma\mapsto\left(\frac{1-[\sigma\cdot\mathbf{1}_{{\mathbb{Z}} ^d\setminus\Lambda}]_{x+T_d}}{2}\right)_{x\in{\mathbb{Z}} ^d}.\]
Lemma 18. Let \(\Omega'=\{X\subset{\mathbb{Z}} ^d\mid -X\subseteq\Lambda-T_d,X\text{ is admissible}\}\). Then the function \(\Phi\) is bijective onto its image \[\Phi(\Omega_\Lambda)=\left\{\mathbb{I}_{-X}\mid X\in\Omega'\right\},\] where we recall that \(\mathbb{I}_{-X}\) denotes the characteristic function of \(-X\).
Proof. We first claim that, for every \(\sigma\in\Omega_\Lambda\), \(\Phi(\sigma)=\mathbb{I}_{-X}\) for some \(X\in\Omega'\). Indeed, \[\Phi(\sigma)_x=\frac{1-[\sigma\cdot\mathbf{1}_{{\mathbb{Z}} ^d\setminus\Lambda}]_{x+T_d}}{2}=\frac{1-\prod_{y\in x+T_d}1}{2}=0\] for all \(x\in{\mathbb{Z}} ^d\setminus\Lambda-T_d\). To see that \(\Phi(\sigma)\) is admissible, it suffices to change one spin at a time starting with \(\mathbf{1}_{\Lambda}\) and take the image under \(\Phi\). Hence, \(\Phi(\Omega_\Lambda)\subseteq\{\mathbb{I}_{-X}\mid X\in\Omega'\}\).
Fix \(X\in\Omega'\). By Definition 7, \(\mathbb{I}_{-X}=\sum_{y\in Y}\mathbb{I}_{y-T_d}\) for a finite set \(Y\subset{\mathbb{Z}} ^d\). Assume for a contradiction that \(Y\not\subseteq\Lambda\). Assume that \(z\in Y\) is such that, for any \(z'\in Y\), if \(z_i'\geqslant z_i\) for all \(i\in\{1,\dots,d\}\), then \(z'=z\) and further assume that \(z_1> 0\) (we proceed analogously for other coordinates). Then \(\mathbb{I}_{-X}(z)=\sum_{y\in Y}\mathbb{I}_{y-T_d}(z)=\mathbb{I}_{z-T_d}(z)=1\), contradicting the assumption \(-X\subseteq\Lambda-T_d\). Assume that \(z\in Y\) is such that \(\langle z,(1,\dots,1)\rangle=\min_{y\in Y}\langle y,(1,\dots,1)\rangle<-n\) and \(z_1\leqslant z_1'\) for any \(z'\in Y\) with \(\langle z',(1,\dots,1)\rangle=\langle z,(1,\dots,1)\rangle\). Then \(\mathbb{I}_{-X}(z-e_1)=\sum_{y\in Y}\mathbb{I}_{y-T_d}(z-e_1)=\mathbb{I}_{z-T_d}(z-e_1)=1\), again contradicting \(-X\subseteq\Lambda-T_d\). This proves that \(Y\subseteq\Lambda\) (recall 6 ). Hence, \(\Phi((-\mathbf{1}_Y)\cdot\mathbf{1}_{\Lambda\setminus Y})=\mathbb{I}_{-X}\), proving surjectivity.
If there are two preimages of \(\mathbb{I}_{-X}\), taking their point-wise product, we obtain \(\sigma\in\Omega_\Lambda\) with \(\Phi(\sigma)=\mathbf{0}\). Proceeding as above, we get that the support of \(\sigma\) is contained in \(\Lambda\), but since \(\varnothing\) is admissible and contained in any translate of \(\Lambda\), this implies that \(\sigma=\mathbf{1}\), yielding the desired injectivity. ◻
By Theorem 8, let \(X\in\Omega'\) be such that any chain from \(X\) to \(\varnothing\) goes through \[\mathcal{X}\mathrel{\vcenter{:}}=\left\{X'\in\Omega'\middle|\;|X'\setminus X|\geqslant n^{\alpha_d}+3|X\setminus X'|/2\right\}.\] By Lemma 18, let \(\sigma=\Phi^{-1}(\mathbb{I}_{-X})\). Let \(\mathcal{G}\subset\Omega_\Lambda\) be the set of configurations which can be reached from \(\mathbf{1}_\Lambda\) by switching the state of one spin in \(\Lambda\) at a time without going through \(\partial\mathcal{G}\mathrel{\vcenter{:}}=\Phi^{-1}(\{\mathbb{I}_{-X'}\mid X'\in\mathcal{X}\})\). In particular, \(\sigma\not\in\mathcal{G}\ni\mathbf{1}_\Lambda\).
Thus, \[\label{eq:finite:lower:var} \mathop{\mathrm{Var}}_{\Lambda}^\mathbf{1}(\mathbb{I}_\mathcal{G})\geqslant\frac{1}{2}\min\left(\mu_\Lambda^\mathbf{1}(\mathbf{1}_\Lambda),\mu_{\Lambda}^\mathbf{1}(\sigma)\right)=\frac{\mu_\Lambda^\mathbf{1}(\sigma)}{2}.\tag{40}\] Moreover, by 2 , for \(X'\in\Omega'\), \[H_\Lambda^\mathbf{1}(\mathbf{1}_\Lambda)=-|\Lambda-T_d|=H_\Lambda^\mathbf{1}\left(\Phi^{-1}\left(\mathbb{I}_{-X'}\right)\right)-2|X'|,\] so \[\label{eq:finite:lower:mu:sigma} \mu^\mathbf{1}_\Lambda(\sigma)=\frac{e^{-2\beta|X|}}{\sum_{X'\in\Omega'}e^{-2\beta|X'|}}\geqslant\frac{e^{-2\beta|X|}}{(1+e^{-2\beta})^{|\Lambda-T_d|}}\geqslant\exp\left(-2\beta|X|-e^{-2\beta}|\Lambda-T_d|\right).\tag{41}\]
Similarly, recalling 4 , we get \[\begin{align} \nonumber\mathcal{D}_\Lambda^\mathbf{1}(\mathbb{I}_\mathcal{G})&{}\leqslant 2 \sum_{x\in\Lambda}\sum_{\sigma'\in\Omega_\Lambda}\mu_\Lambda^\mathbf{1}(\sigma')\mathbb{I}_{\sigma'\not\in\mathcal{G}\ni(\sigma')^x}\leqslant 2|\Lambda|\mu^\mathbf{1}_\Lambda(\partial\mathcal{G})=2|\Lambda|\sum_{X'\in\mathcal{X}}\mu_\Lambda^\mathbf{1}\left(\Phi^{-1}\left(\mathbb{I}_{-X'}\right)\right)\\ \nonumber&{}\leqslant 2|\Lambda|\sum_{X'\in\mathcal{X}}e^{-2\beta|X'|}\leqslant 2|\Lambda|e^{-2\beta|X|}\sum_{X'\in\mathcal{X}}e^{-2\beta(4n^{\alpha_d}/5+(|X'\setminus X|+|X\setminus X'|)/5)}\\ \nonumber&{}\leqslant 2|\Lambda|e^{-2\beta(|X|+4n^{\alpha_d}/5)}\sum_{Y\subseteq\Lambda-T_d}e^{-2\beta|Y|/5}=2|\Lambda|e^{-2\beta(|X|+4n^{\alpha_d}/5)}\left(1+e^{-2\beta/5}\right)^{|\Lambda-T_d|}\\ &{}\leqslant 2|\Lambda|\exp\left(-2\beta(|X|+4n^{\alpha_d}/5)+|\Lambda-T_d|e^{-2\beta/5}\right). \label{eq:finite:lower:D} \end{align}\tag{42}\]
Finally, recalling 5 and combining 40 , 41 and 42 , we obtain \[\begin{align} T_{\mathrm{rel}}^{\Lambda,\mathbf{1}}&{}\geqslant\frac{\exp(8\beta n^{\alpha_d}/5-|\Lambda-T_d|e^{-2\beta/5}-|\Lambda-T_d|e^{-2\beta})}{4|\Lambda|}\\&{}\geqslant\exp\left(8\beta n^{\alpha_d}/5-2^{d+1}-(d+2)\ln(2d)-2\beta/5\right)\geqslant\frac{e^{\beta n^{\alpha_d}}}{C}, \end{align}\] since \(|\Lambda|\leqslant|\Lambda-T_d|\leqslant(2n)^d\leqslant 2^de^{2\beta/5}\), taking \(C=(2d)^{d+2}e^{2^{d+1}}\). ◻
In this section we prove the upper bound of Theorem 2 via the bisection and canonical paths techniques.
Proof of the upper bound of Theorem 2. Fix \(d=2\) and \(\beta>0\). The result follows immediately, once we prove by induction that, for \(k\geqslant 0\) \[T_{\mathrm{rel}}^{{\mathbb{T}} _{2^k}}\leqslant\left(41472\cdot e^{2\beta}\right)^{k}.\] The base case \(k=0\) is clear and uses that \(d\) is even (the dynamics is not irreducible for \(d\) odd).
Assume \(k\geqslant 1\) and let \(n=2^k\). Recall \(\Phi\) and \(p\) from 34 and 36 . As in Section 3.1, let \(\pi\) denote the product Bernoulli measure with parameter \(p\) on \(\Omega'_{{\mathbb{T}} _{n}}={\mathbb{F}} _2^{{\mathbb{T}} _n}\) and \(\mathop{\mathrm{Var}}'\) be its variance. By Lemma 16, we may reason in terms of the lamp Markov chain with generator as in 35 . For \(r=(r_1,\dots,r_d)\in\{0,1\}^d\), let \[{\mathbb{T}} ^r_n=\left\{(x_1,\dots,x_d)\in{\mathbb{T}} _n\mid\forall i\in\{1,\dots,d\},x_i=r_i\pmod 2\right\}=r+2{\mathbb{T}} _n,\] taking into account that \(n\) is even. For \(A\subseteq{\mathbb{T}} _n\) and \(\omega\in\Omega'_{{\mathbb{T}} _n}\), we write \(\omega^A=\omega+\mathbb{I}_A\). Let \[\begin{align} \mathcal{L}^rf(\omega)&{}=\sum_{x\in{\mathbb{T}} _n^r}\frac{\pi(\omega^{x-2T_d})}{\pi(\omega)+\pi(\omega^{x-2T_d})}\left(f\left(\omega^{x-2T_d}\right)-f(\omega)\right),\\ \tag{43}\mathcal{D}^r(f)&{}=\sum_{\omega\in\Omega'_{{\mathbb{T}} _n}}\sum_{x\in{\mathbb{T}} _n^r}\frac{\pi(\omega)\pi(\omega^{x-2T_d})}{\pi(\omega)+\pi(\omega^{x-2T_d})}\left(f\left(\omega^{x-2T_d}\right)-f(\omega)\right)^2=-\pi\left(f\mathcal{L}^rf\right),\\ \tag{44}T_{\mathrm{rel}}^r&{}=\left(\inf\left\{\frac{\mathcal{D}^r(\pi_{{\mathbb{T}} _n\setminus{\mathbb{T}} _n^r}(f))}{\mathop{\mathrm{Var}}'\left(\pi_{{\mathbb{T}} _n\setminus{\mathbb{T}} _n^r}(f)\right)}\middle| f\colon\Omega'_{{\mathbb{T}} _n}\to{\mathbb{R}} ,\mathop{\mathrm{Var}}'(\pi_{{\mathbb{T}} _n\setminus{\mathbb{T}} _n^r}(f))\neq0\right\}\right)^{-1}. \end{align}\]
Fix \(f\colon\Omega'_{{\mathbb{T}} _n}\to{\mathbb{R}}\). By the Poincaré inequality for independent Glauber dynamics (the Efron–Stein inequality) and 44 , we have \[\label{eq:var39}\mathop{\mathrm{Var}}'(f)\leqslant\sum_{r\in\{0,1\}^d}\pi\left(\mathop{\mathrm{Var}}'_{{\mathbb{T}} ^r_{n}}(f)\right)\leqslant\sup_{r\in\{0,1\}^d}T_{\mathrm{rel}}^{r}\sum_{r\in\{0,1\}^d}\pi\left(\mathcal{D}^r(f)\right)=T_{\mathrm{rel}}^{{\mathbb{T}} _{n/2}}\sum_{r\in\{0,1\}^d}\pi\left(\mathcal{D}^r(f)\right).\tag{45}\]
We finally seek to relate each term in \(\pi(\mathcal{D}^r(f))\) to appropriate terms featuring in \(\mathcal{D}(f)\). Consider \(r\in\{0,1\}^d\), \(x\in{\mathbb{T}} ^r_{n}\) and \(\omega\in\Omega'_{{\mathbb{T}} _n}\). Set \(\Omega''=\{\omega\in\Omega'_{{\mathbb{T}} _n}\mid\sum_{y\in x-2T_d}\omega_y\in\{2,3\}\}\). In view of the symmetry between \(\omega\) and \(\omega^{x-2T_d}\) in 43 , we may assume \(\omega\in\Omega''\). Define the configurations \(\omega=\omega^{(0)},\omega^{(1)},\omega^{(2)},\omega^{(3)}=\omega^{x-2T_d}\) as indicated in Figure 5. The Cauchy–Schwarz inequality gives \[\left(f\left(\omega^{x-2T_d}\right)-f(\omega)\right)^2\leqslant 3\sum_{i=1}^3\left(f\left(\omega^{(i)}\right)-f\left(\omega^{(i-1)}\right)\right)^2.\] Summing over \(\omega\in\Omega''\), we get \[\begin{align} \label{eq:paths} &2\sum_{\omega\in\Omega''}\frac{\pi(\omega)\pi(\omega^{x-2T_d})}{\pi(\omega)+\pi(\omega^{x-2T_d})}\left(f\left(\omega^{x-2T_d}\right)-f(\omega)\right)^2\leqslant 2\sum_{\omega\in\Omega''}\pi(\omega)\left(f\left(\omega^{x-2T_d}\right)-f(\omega)\right)^2\\ \nonumber&{}\leqslant 6\sum_{i=1}^3\sum_{\omega\in\Omega''}\pi(\omega)\left(f\left(\omega^{(i)}\right)-f\left(\omega^{(i-1)}\right)\right)^2 \underbrace{\sum_{y\in x-T_d}\mathbb{I}_{\omega^{(i-1)}=(\omega^{(i)})^{y-T_d}}}_{=1}\underbrace{ \sum_{\omega'\in \Omega'_{{\mathbb{T}} _n}}\mathbb{I}_{\omega'=\omega^{(i)}}\frac{\pi(\omega')\pi((\omega')^{y-T_d})}{\pi(\omega')\pi((\omega')^{y-T_d})}}_{=1}\\ &{}= 6 \begin{multlined}[t][\dimexpr\textwidth - \lhswidth - 0.5em] \sum_{\omega'\in\Omega'_{{\mathbb{T}} _n}}\sum_{y\in x-T_d}\left(f\left(\omega'\right)-f\left(\left(\omega'\right)^{y-T_d}\right)\right)^2\frac{\pi(\omega')\pi((\omega')^{y-T_d})}{\pi(\omega')+\pi((\omega')^{y-T_d})}\\ \times\underbrace{\sum_{i=1}^3\sum_{\omega\in\Omega''}\mathbb{I}_{\omega^{(i)}=\omega',\omega^{(i-1)}=(\omega')^{y-T_d}}\frac{\pi(\omega)(\pi(\omega')+\pi((\omega')^{y-T_d}))}{\pi(\omega')\pi((\omega')^{y-T_d})}}_{=:K(\omega,\omega',i,y)}.\end{multlined} \nonumber \end{align}\tag{46}\]
A direct inspection of Figure 5 shows that, for any \(\omega\in\Omega''\) and any \(j\in\{0,\dots,3\}\) we have \(\|\omega^{(j)}\|_1-\|\omega\|_1\leqslant 1\) (this property is specific to \(d=2\)). Consequently, \[\begin{align} \sup_{\omega'\in\Omega'_{{\mathbb{T}} _n},y\in x-T_d}K\left(\omega,\omega',i,y\right)&{}\leqslant\frac{2}{p}\sup_{\omega'\in\Omega'_{{\mathbb{T}} _n}}\sum_{i=1}^3\sum_{\omega\in\Omega''}\mathbb{I}_{\omega^{(i)}=\omega'}\\ &{}\leqslant\frac{6}{p}\left|\left\{\omega\in\Omega''\middle|\omega_{{\mathbb{T}} _n\setminus (x-T_d-T_d)}=\omega'_{{\mathbb{T}} _n\setminus(x-T_d-T_d)}\right\}\right|=\frac{6\cdot 2^6}{2p}. \end{align}\]
Combining this with 5 , 4 , 45 , 46 , we obtain \[T_{\mathrm{rel}}^{{\mathbb{T}} _n}\leqslant T_{\mathrm{rel}}^{{\mathbb{T}} _{n/2}}\cdot \frac{1152}{p}\cdot6\cdot\sup_{y\in{\mathbb{T}} _n}\left|\left\{x\in{\mathbb{T}} _n\middle| y\in x-T_d\right\}\right|=T_{\mathrm{rel}}^{{\mathbb{T}} _{n/2}}\cdot\frac{20736}{p},\] completing the induction step in view of 36 . ◻
In this section we prove the upper bound of Theorem 1. The main ingredient is Theorem 9. As a warm-up, we prove the upper bound of Theorem 3, which is an immediate corollary of the following trivial bound (see Martinelli99?*Theorem 3.8 for a slightly better one which gives \(n^{d-1}\) instead of \(n^d\) in Theorem 3). We include its proof for completeness.
Lemma 19 (Finite volume upper bound). For any \(d\geqslant 2\), \(n\geqslant 1\), \(\beta>0\) and \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus(-\Lambda_{d,n})}\), we have \[T_{\mathrm{rel}}^{-\Lambda_{d,n},\tau}\leqslant e^{2(\beta+1)|\Lambda_{d,n+3}|}.\]
Proof. Recalling 5 , we consider \(f\colon\Omega_{-\Lambda_{d,n}}\to{\mathbb{R}}\) and may assume that there exist \(\sigma,\sigma'\in\Omega_{-\Lambda_{d,n}}\) with \(f(\sigma)=0\leqslant f(\sigma'')\leqslant f(\sigma')=1\) for all \(\sigma''\in\Omega_{-\Lambda_{d,n}}\). In particular, \(\mathop{\mathrm{Var}}_{-\Lambda_{d,n}}^\tau(f)\leqslant 1\). On the other hand, we may fix a sequence \(\sigma^{(0)},\dots,\sigma^{(m)}\) such that \(\sigma^{(0)}=\sigma\), \(\sigma^{(m)}=\sigma'\), \(\|\sigma^{(i)}-\sigma^{(i-1)}\|_1=1\) for \(i\in\{1,\dots,m\}\) (we change one spin at a time) and \(m\leqslant|\Lambda_{d,n}|\). Then there exists \(k\in\{1,\dots,m\}\) such that \(f(\sigma^{(k)})-f(\sigma^{(k-1)})\geqslant 1/m\). Write \(\sigma^{(k)}-\sigma^{(k-1)}=2\mathbb{I}_x\), with \(\mathbb{I}_x\) the characteristic function of \(x\in-\Lambda_{d,n}\). Then \[\begin{align} \mathcal{D}^{\tau}_{-\Lambda_{d,n}}(f)&{}\overset{\eqref{eq:def:D}}\geqslant\mu_{-\Lambda_{d,n}}^\tau\left(\sigma^{(k)}\right)\mathop{\mathrm{Var}}_x(f)\left(\sigma^{(k)}\right)\geqslant\mu_{-\Lambda_{d,n}}^\tau\left(\sigma^{(k)}\right)\cdot\frac{\min(\mu^\tau_{-\Lambda_{d,n}}(\sigma^{(k)}),\mu^\tau_{-\Lambda_{d,n}}(\sigma^{(k-1)}))}{2m^2(\mu^\tau_{-\Lambda_{d,n}}(\sigma^{(k)})+\mu^\tau_{-\Lambda_{d,n}}(\sigma^{(k-1)}))} \\ &{}\overset{\eqref{eq:def:H}}\geqslant\frac{\mu^\tau_{-\Lambda_{d,n}}(\sigma^{(k)})}{2m^2(1+e^{2(d+1)\beta})}\geqslant\frac{\min\{\mu^\tau_{-\Lambda_{d,n}}(\bar\sigma)\mid\bar \sigma\in\Omega_{-\Lambda_{d,n}}\}}{4|\Lambda_{d,n}|^2e^{2(d+1)\beta}} \\ &{} \overset{\eqref{eq:def:H}}\geqslant\frac{\min_{\bar\sigma}e^{-\beta H^\tau_{-\Lambda_{d,n}}(\bar\sigma)}}{|\Omega_{-\Lambda_{d,n}}|\max_{\bar\sigma} e^{-\beta H^\tau_{-\Lambda_{d,n}}(\bar\sigma)}}\cdot\frac{1}{4|\Lambda_{d,n}|^2e^{2(d+1)\beta}}\overset{\eqref{eq:def:H}}\geqslant\frac{e^{-2\beta|-\Lambda_{d,n}-T_d|}}{|\Omega_{-\Lambda_{d,n}}|\cdot |\Lambda_{d,n}|^2\cdot4e^{2(d+1)\beta}}\\ &{}\geqslant\frac{e^{-(2\beta+2)|\Lambda_{d,n+1}|}}{4e^{2(d+1)\beta}} \end{align}\] using \(n\geqslant 1,d\geqslant 2\) in the last inequality. Plugging this into 5 and using \(n\geqslant 1,d\geqslant 2\) again, we obtain the desired inequality. ◻
The upper bound of Theorem 1 will use Lemma 19 at a suitably chosen scale of order \(e^{C\beta}\) with large \(C\). Above this scale we will be able to apply Theorem 9 to deduce that boundary conditions have little influence on the partition function and so also on the Boltzmann measure. This will allow us to use the path coupling method in finite volume and then transfer fast convergence to equilibrium to infinite volume.
Proof of the upper bound of Theorem 1. The starting point is the same high temperature expansion as in Chleboun17?*Section 4.0.1, which we recall for completeness. Fix a positive integer \(d\), \(\beta>0\) and set \(\theta=\tanh(\beta)\). For any \(u\in\{-1,1\}\), we have \[e^{\beta u}=\cosh(\beta)(1+u\theta).\] Consequently, for any finite \(\Lambda\subset{\mathbb{Z}} ^d\) and \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda}\), we have \[\begin{align} Z^\tau_\Lambda&{}=\sum_{\sigma\in\Omega_\Lambda}\prod_{x\in\Lambda-T_d}e^{\beta[\sigma\cdot\tau]_{x+T_d}}=(\cosh(\beta))^{|\Lambda-T_d|}\sum_{X\subset\Lambda-T_d}\theta^{|X|}\sum_{\sigma\in\Omega_\Lambda}\prod_{x\in X}[\sigma\cdot\tau]_{x+T_d}\\ &{}=(\cosh(\beta))^{|\Lambda-T_d|}\sum_{X\in\mathcal{C}_{(-T_d,\Lambda-T_d)}}\theta^{|X|}\sum_{\sigma\in\Omega_\Lambda}\prod_{x\in X}[\sigma\cdot\tau]_{x+T_d}, \end{align}\]
recalling the set of \((A,B)\)-cycles \(\mathcal{C}_{A,B}\) from Section 2.3.1. We set \(\mathcal{Z}_\Lambda^\tau=Z_\Lambda^\tau2^{-|\Lambda|}(\cosh(\beta))^{-|\Lambda-T_d|}\) to obtain \[\max_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda}}\left|\mathcal{Z}_\Lambda^\tau-1\right|\leqslant\sum_{X\in\mathcal{C}_{(-T_d,\Lambda-T_d)}\setminus\{\varnothing\}}\theta^{|X|}.\]
Recalling 6 and 12 , we get \[\max_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus (-\Lambda_{d,n})}}\big|\mathcal{Z}^\tau_{-\Lambda_{d,n}}-1\big|\leqslant G_{d,n+1}(\theta)-1\] for any positive integer \(n\).
The rest of the proof diverges from Chleboun17?. Set \(m=\lceil1/(1-\theta)\rceil\geqslant 2\) and fix a positive integer \(a_d\geqslant d^2\) such that Theorem 9 holds and \[\label{eq:nNconditions} 2^{a_d-d^2}\geqslant 4(6d)^{d(d+1)}\tag{47}\] (the latter is clearly true for \(a_d\) large enough). Then Theorem 9 gives \[\label{eq:Z:close:to:1} n\geqslant m^{a_d}\Longrightarrow \max_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus(-\Lambda_{d,n})}}\left|\mathcal{Z}^\tau_{-\Lambda_{d,n}}-1\right|\leqslant n^{-a_d}.\tag{48}\] This result will combine with the following observation.
Lemma 20 (Boundary condition influence). Fix finite \(\Lambda'\subset\Lambda\subset{\mathbb{Z}} ^d\) and \(x\in{\mathbb{Z}} ^d\setminus\Lambda\) such that \(\min_{y\in\Lambda\setminus\Lambda'}\|x-y\|_1\geqslant 3\). For \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda}\), define the marginal of \(\mu^\tau_\Lambda\) on \(\Lambda\setminus\Lambda'\) by \[\mu^\tau_{\Lambda\setminus\Lambda'}\colon\Omega_{\Lambda\setminus\Lambda'}\to{\mathbb{R}} \colon\sigma\mapsto \sum_{\eta\in\Omega_{\Lambda'}}\mu^\tau_\Lambda(\eta\cdot\sigma).\] Then, for any \(\tau\in\Omega_{\Lambda\setminus\Lambda'}\), the marginal’s Radon–Nikodým derivative is uniformly bounded by \[\frac{\mathrm d\mu^\tau_{\Lambda\setminus\Lambda'}}{\mathrm d\mu^{\tau^x}_{\Lambda\setminus\Lambda'}}\geqslant\frac{\mathcal{Z}^{\tau^x}_\Lambda}{\mathcal{Z}^\tau_\Lambda}\min_{\sigma\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda'}}\frac{\mathcal{Z}^{\sigma^x}_{\Lambda'}}{\mathcal{Z}^{\sigma}_{\Lambda'}}.\]
Proof. Fix \(\tau\in\Omega_{{\mathbb{Z}} ^d\setminus\Lambda}\). For \(\sigma\in\Omega_{\Lambda\setminus\Lambda'}\) and \(\eta'\in\Omega_{\Lambda'}\), we have \[\begin{align} \mu^\tau_{\Lambda\setminus\Lambda'}(\sigma)&{}=\sum_{\eta\in\Omega_{\Lambda'}}\mu^\tau_\Lambda(\eta\cdot\sigma)=\sum_{\eta\in\Omega_{\Lambda'}}\frac{e^{-\beta H_\Lambda^\tau(\eta\cdot\sigma)}}{Z^\tau_\Lambda}= \frac{1}{Z^{\tau}_\Lambda}\sum_{\eta\in\Omega_{\Lambda'}}e^{-\beta H_{\Lambda'}^{\tau\cdot\sigma}(\eta)}e^{\beta\sum_{y\in(\Lambda-T_d)\setminus(\Lambda'-T_d)}[\tau\cdot\sigma\cdot\eta]_{y+T_d}}\\ &{}=\frac{Z_{\Lambda'}^{\tau\cdot\sigma}}{Z^{\tau}_\Lambda}e^{\beta\sum_{y\in(\Lambda-T_d)\setminus(\Lambda'-T_d)}[\tau^x\cdot\sigma\cdot \eta']_{y+T_d}}=\frac{Z^{\tau\cdot\sigma}_{\Lambda'}Z^{\tau^x}_\Lambda}{Z^\tau_\Lambda Z^{\tau^x\cdot\sigma}_{\Lambda'}}\mu^{\tau^x}_{\Lambda\setminus\Lambda'}(\sigma)=\frac{\mathcal{Z}^{\tau\cdot\sigma}_{\Lambda'}\mathcal{Z}^{\tau^x}_\Lambda}{\mathcal{Z}^{\tau^x\cdot\sigma}_{\Lambda'}\mathcal{Z}^\tau_\Lambda}\mu^{\tau^x}_{\Lambda\setminus\Lambda'}(\sigma), \end{align}\] where we used the fact that \(\sum_{y\in(\Lambda-T_d)\setminus(\Lambda'-T_d)}[\tau\cdot\sigma\cdot\eta]\) does not depend on \(\tau_x\) nor on \(\eta\) because of the range of \(y\) and the assumption on \(x\). ◻
Fix \[\begin{align} \label{eq:def:N} n&{}=m^{a_d},&N&{}=(3d)^{d+1}(n+1)^d, \end{align}\tag{49}\] so that \[\begin{align} \tag{50} N^d/(3d^d)&{}\geqslant 3^dN^{d-1}d(n+1)^d\geqslant(n+1)^d(d+1)(N+2)^{d-1}\\ \nonumber&{}\geqslant(N+2)^{2d-1}(d+1)N^{-d}\geqslant(N+2)^{2d-1}(d+1)(6d)^{-(d+1)d}n^{-d^2}\\ \tag{51}&{}\geqslant(N+1)^d(N+2)^{d-1}(d+1)4n^{-a_d}, \end{align}\] thanks to 47 in the last inequality. Let us consider the following block dynamics on a \(d\)-dimensional torus \({\mathbb{T}} ={\mathbb{T}} _L\) of side length \(L>N+1\) which we first describe via its graphical representation. Endow each \(x\in{\mathbb{T}}\) with a Poisson clock of unit rate which rings at times \(P_x\subset(0,\infty)\). When the clock at \(x\) rings, change the current configuration \(\sigma\) to \(\tau\cdot\sigma'\) with \(\tau=\sigma_{{\mathbb{T}} \setminus (x-\Lambda_{d,N})}\) and \(\sigma'\) chosen at random according to \(\mu_{x-\Lambda_{d,N}}^{\tau}\); that is, we fully resample the configuration in a simplex region conditionally on the configuration outside of it. In other words, this is the Markov process with generator and Dirichlet form given by \[\begin{align} \mathcal{L}_{\mathbb{T}} f&{}=\sum_{x\in{\mathbb{T}} }\left(\mu_{x-\Lambda_{d,N}}(f)-f\right),&\mathcal{D}_{\mathbb{T}} (f)&{}=\sum_{x\in{\mathbb{T}} }\mu_{\mathbb{T}} \left(\mathop{\mathrm{Var}}_{x-\Lambda_{d,N}}(f)\right). \end{align}\]
We next resort to the path coupling method (see Levin09?*Chapter 14). It is used to show that a Markov chain (the block dynamics we just defined) converges to its invariant measure exponentially fast. To be able to apply it, we need to couple two copies of the chain with different initial conditions in such a way that their \(\ell^1\)-distance (number of discrepancies) is contracted in expectation. By a well-known reduction, it suffices to treat initial conditions that only differ at a single site.
Consider a configuration \(\sigma\in\Omega_{\mathbb{T}}\) and \(x\in{\mathbb{T}}\). Recall that \(\sigma^x\) is the configuration that only differs from \(\sigma\) at \(x\). We construct a coupling between the Markov process \((\sigma(t))_{t\geqslant 0}\) with initial state \(\sigma\) and the one with initial state \(\sigma^x\) denoted by \((\sigma'(t))_{t\geqslant 0}\) as follows. Let \((\xi_{y,t})_{y\in{\mathbb{T}} ,t\in P_y}\) be a sequence of i.i.d.Bernoulli random variables with parameter \(1-4n^{-a_d}\). We use the same Poisson clocks \((P_y)_{y\in{\mathbb{T}} }\) for both processes. Both processes remain constant outside the clock ring times \(\bigcup_{y\in {\mathbb{T}} }P_y\). Assume the coupling is constructed until the clock ring time \(t\in P_y\) for some \(y\in{\mathbb{T}}\). Set \(\partial\Lambda=y+((-\Lambda_{d,N}+(-T_d)+T_d)\setminus(-\Lambda_{d,N}))\), \(\tau=\sigma_{{\mathbb{T}} \setminus(y-\Lambda_{d,N})}(t-)\), \(\tau'=\sigma_{{\mathbb{T}} \setminus(y-\Lambda_{d,N})}'(t-)\). If \(x\in\partial\Lambda\), fix \(z\in{\mathbb{T}}\) so that \(z-\Lambda_{d,n}\subset y-\Lambda_{d,N}\) and \(x\) is at \(\ell^1\)-distance at least 3 from \((y-\Lambda_{d,N})\setminus(z-\Lambda_{d,n})\), see Figure 6. We construct \(\sigma(t)\) and \(\sigma'(t)\) by the following rules.
If \(x\in {\mathbb{T}} \setminus(y-\Lambda_{d,N}+(-T_d)+T_d)\) and \(\tau^x=\tau'\), we define \(\sigma(t)=\tau\cdot\eta\) and \(\sigma'(t)=\tau'\cdot\eta\) with \(\eta\) distributed according to \(\mu^\tau_{y-\Lambda_{d,N}}=\mu^{\tau'}_{y-\Lambda_{d,N}}\). That is, if \(x\) is far from the update block, it does not influence the update.
If \(x\in y-\Lambda_{d,N}\) and \(\tau=\tau'\), we define \(\sigma(t)=\sigma'(t)=\tau\cdot\eta\) with \(\eta\) distributed according to \(\mu^\tau_{y-\Lambda_{d,N}}=\mu^{\tau'}_{y-\Lambda_{d,N}}\). That is, if \(x\) is inside the block and the boundary conditions coincide, the two processes coalesce.
If \(\tau'=\tau^x\), \(x\in\partial\Lambda\) and \(\xi_{y,t}=1\), we define \(\sigma(t)=\tau\cdot\kappa\cdot\eta\) and \(\sigma'(t)=\tau'\cdot\kappa\cdot\eta'\) with \(\kappa\) distributed according to \(\mu^\tau_{(y-\Lambda_{d,N})\setminus(z-\Lambda_{d,n})}\) and, conditionally on \(\kappa\), \((\eta,\eta')\) distributed according to \(\mu^{\tau\cdot\kappa}_{z-\Lambda_{d,n}}\otimes\mu^{\tau'\cdot\kappa}_{z-\Lambda_{d,n}}\). That is, if \(x\) is on the boundary of the block, the boundary conditions coincide except at \(x\), and we are lucky, we ensure that the two processes coincide outside the small simplex close to \(x\).
If \(\tau'=\tau^x\), \(x\in\partial\Lambda\) and \(\xi_{y,t}=0\), we define \(\sigma(t)=\tau\cdot\eta\) and \(\sigma'(t)=\tau'\cdot\eta'\) with \(\eta\) and \(\eta'\) independent with distributions such that the total marginals with the previous case are \(\mu^\tau_{y-\Lambda_{d,N}}\) and \(\mu^{\tau'}_{y-\Lambda_{d,N}}\). That is, if \(x\) is on the boundary of the block, the boundary conditions coincide except at \(x\), and we are unlucky, we perform the update somehow.
Otherwise, we define \(\sigma(t)=\tau\cdot\eta\) and \(\sigma'(t)=\tau'\cdot\eta'\) with \((\eta,\eta')\) with law \(\mu^\tau_{y-\Lambda_{d,N}}\otimes\mu^{\tau'}_{y-\Lambda_{d,N}}\). That is, in all other cases, the two processes perform the block update independently.
We note that this construction is possible, since Lemma 20 and 48 give \[\begin{align} \frac{\mathrm d\mu_{(y-\Lambda_{d,N})\setminus(z-\Lambda_{d,n})}^{\tau^x}}{\mathrm d\mu_{(y-\Lambda_{d,N})\setminus(z-\Lambda_{d,n})}^{\tau}}&{}\geqslant\frac{1-N^{-a_d}}{1+N^{-a_d}}\cdot\frac{1-n^{-a_d}}{1+n^{-a_d}}\overset{N\geqslant n}\geqslant\left(\frac{1-n^{-a_d}}{1+n^{-a_d}}\right)^2\\ &{}\geqslant(1-n^{-a_d})^4\geqslant 1-4n^{-a_d}{=\mathbb{P}(\xi_{y,t}=1)}. \end{align}\] We next claim that this coupling is contracting for the natural graph metric on \(\Omega_{\mathbb{T}}\).
Lemma 21 (Contraction). Under the above coupling, setting \(\phi\colon[0,\infty]\colon t\mapsto {\mathbb{E}} \|\sigma(t)-\sigma'(t)\|_1/2\), we have \[\phi'(0)\leqslant-\frac{N^d}{3d^d}\phi(0).\]
Proof. By construction, \(\phi(0)=1\), since \((\sigma(0))^x=\sigma'(0)\). We consider each of the first four cases in the coupling separately (the fifth one requires two updates to occur, so it does not feature in the derivative at 0). Notice that, in these cases, \(\phi\) increases by \(0\); \(-1\); at most \(|\Lambda_{d,n}|\); at most \(|\Lambda_{d,N}|\), respectively. Hence, taking into account the rate at which each updates of each type occur, we get \[\begin{align} &\phi'(0)\\ &{}\leqslant 0\cdot (|{\mathbb{T}} |-|\Lambda_{d,N}+(-T_d)+T_d|) - 1 \cdot|\Lambda_{d,N}|+|\Lambda_{d,n}|\cdot|\partial\Lambda|\cdot(1-4n^{-a_d})+|\Lambda_{d,N}|\cdot|\partial\Lambda|\cdot 4n^{-a_d}\\ &{}\leqslant-(N/d)^d+(n+1)^d\cdot (d+1)(N+2)^{d-1}+(N+1)^d\cdot (d+1)(N+2)^{d-1}\cdot 4n^{-a_d}\\ &{}\leqslant-N^d/(3d^d), \end{align}\] recalling 50 and 51 in the last inequality. ◻
By the standard technique of path coupling (see Levin09?*Theorem 14.6, which applies similarly to the continuous time setting), Lemma 21 implies that, for any probability measure \(\nu\) on \(\Omega_{\mathbb{T}}\) and \(t\geqslant 0\) \[\label{eq:dTV} d_{\mathrm{TV}}\left(\nu P^{\mathbb{T}} _t,\mu_{\mathbb{T}} \right)\leqslant W_1\left(\nu P_t^{\mathbb{T}} ,\mu_{\mathbb{T}} \right)\leqslant e^{-tN^d/(3d^d)}|{\mathbb{T}} |,\tag{52}\] where \(P_t^{\mathbb{T}}\) denotes the semi-group of \(\mathcal{L}_{\mathbb{T}}\) and \(2W_1\) is the \(\ell_1\)-Wasserstein distance and \(d_{\mathrm{TV}}\) is the total variation distance.
We next seek to compare the invariant measure \(\mu_{\mathbb{T}}\) on the torus to the semi-group \(P^{(N)}_t\) of the block dynamics on \({\mathbb{Z}} ^d\), whose Dirichlet form is given by \[\label{eq:def:DN} \mathcal{D}^{(N)}(f)=\sum_{x\in{\mathbb{Z}} ^d}\mu\left(\mathop{\mathrm{Var}}_{x-\Lambda_{d,N}}(f)\right)\tag{53}\] and whose relaxation time we denote by \(T_{\mathrm{rel}}^{(N)}\) (recall 5 ). This follows from the classical fact that, in a finite-range model, information cannot travel too fast, as we recall next.
Lemma 22 (Exponential ergodicity). Let \(a\) be a positive integer, \(f\) be a local function with \(\mathop{\mathrm{supp}}f \subset[-a,a)^d\) and \(\sup f-\inf f\leqslant 1\). Then, for any \(t\geqslant a\) large enough depending on \(d\), if the size \(L\) of \({\mathbb{T}} ={\mathbb{T}} _L\) satisfies \(L/t\in[(5N)^{5d},(6N)^{5d}]\), we have \[\left\|P_t^{(N)}f-\mu_{\mathbb{T}} (f)\right\|_\infty\leqslant e^{-tN^d/(4d^d)}.\]
Proof. We define \(f_{\mathbb{T}} \colon \Omega_{\mathbb{T}} \to {\mathbb{R}} :\sigma^{{\mathbb{T}} }\mapsto f(\sigma)\), where \(\sigma\) is any configuration in \(\Omega\) such that \(\sigma_x=\sigma^{\mathbb{T}} _x\) for \(x\in[-L/2,L/2)^d\cap{\mathbb{Z}} ^d\equiv{\mathbb{T}}\)(we assume \(L\) divisible by 2 for simplicity).
Consider a coupling of the block dynamics \((\sigma^{\mathbb{T}} (t))_{t\geqslant 0}\) on \({\mathbb{T}}\) and the one on \({\mathbb{Z}} ^d\) denoted \((\sigma(t))_{t\geqslant 0}\) with consistent initial conditions \(\sigma^{\mathbb{T}} (0)=\sigma(0)_{[-L/2,L/2)^d\cap{\mathbb{Z}} ^d}\), constructed as follows. Let \((P_x)_{x\in{\mathbb{Z}} ^d}\) be Poisson processes on \([0,\infty)\) with unit intensity. At times \(t\in P_x\) for \(x\in{\mathbb{T}} \equiv[-L/2,L/2)^d\cap{\mathbb{Z}} ^d\), we update both \(\sigma^{\mathbb{T}} (t-)\) and \(\sigma(t-)\) to the same state in the block \(x-\Lambda_{d,N}\), if the two processes have the same boundary condition on the boundary of this simplex (in the sense of Figure 6). If the boundary conditions do not match, the updates are performed independently. The updates at sites in \({\mathbb{Z}} ^d\setminus [-L/2,L/2)^d\) are performed only in \(\sigma(t)\), governed by independent Poisson processes.
Clearly, if \(f(\sigma(t))\neq f_{\mathbb{T}} (\sigma^{\mathbb{T}} (t))\) for some \(t\geqslant 0\), then there exists a sequence \((x_i)_{i=0}^m\) of points in \([-L/2,L/2)^d\cap{\mathbb{Z}} ^d\) and decreasing sequence of times \((t_i)_{i=0}^m\) such that: \(t_0\leqslant t\); \(\mathop{\mathrm{supp}}f\cap(x_0-\Lambda_{d,N})\neq \varnothing\); \(x_m-\Lambda_{d,N}+\Lambda_{d,N}-T_d+T_d\not \subset [-L/2,L/2)^d\); for all \(i\in\{1,\dots,m\}\) we have \((x_{i-1}-\Lambda_{d,N}-T_d+T_d)\cap(x_i-\Lambda_{d,N})\neq\varnothing\); for all \(i\in\{0,\dots,m\}\) we have \(t_i\in P_{x_i}\). Note that any such sequence satisfies \(m+2\geqslant(L/2-a)/(2N+1)\geqslant L/(5N)+1\).
The probability that clocks ring in order along a given sequence \((x_i)_{i=0}^k\) vertices before time \(t\) is exactly \({\mathbb{P}} (N_t\geqslant k+1)\), where \(N_t\) denotes a Poisson random variable of mean \(t\). In order to bound it, we recall the Bennett inequality: for \(\alpha\geqslant 1\), \[{\mathbb{P}} (N_t\geqslant t\alpha)\leqslant e^{-t(\alpha\log(\alpha)-\alpha+1)},\] which follows from Markov’s inequality applied to the moment generating function \({\mathbb{E}} [e^{xN_t}]=e^{(e^x-1)t}\). Hence, \[\begin{align} \nonumber\left\|P^{\mathbb{T}} _tf_{\mathbb{T}} -P^{(N)}_tf\right\|_\infty&{}\leqslant(2a)^d|\Lambda_{d,N}-\Lambda_{d,N}+T_d-T_d|^{L/(5N)}{\mathbb{P}} (N_t\geqslant L/(5N))\\ \nonumber&{}\leqslant(2a)^d(N+2)^{2dL/(5N)}e^{-L/(5N)(\ln (L/(5Nt))-1)}\\ \label{eq:PtTZd}&{}\leqslant(2a)^de^{-2L/(5N)}\leqslant e^{-L/(5N)}, \end{align}\tag{54}\] since \(L\geqslant(5N)^{5d}t\geqslant 5e^3tN(N+2)^{2d}\) and \(e^{L/(5N)}\geqslant e^{2dt}\geqslant e^{2da}\geqslant(2a)^d\). Then, \[\begin{align} \left\|P^{(N)}_tf-\mu_{\mathbb{T}} (f)\right\|_\infty&{}\overset{\eqref{eq:PtTZd}}\leqslant e^{-L/(5N)}+\left\|P_t^{\mathbb{T}} f-\mu_{\mathbb{T}} (f)\right\|_\infty\overset{\eqref{eq:dTV}}\leqslant e^{-t(5N)^{5d-1}}+e^{-tN^d/(3d^d)}L^d\\ &{}\leqslant e^{-tN^d}+e^{-t N^d/(3d^d)}(6N)^{5d^2}t^d\leqslant e^{-tN^d/(4d^d)}.\qedhere \end{align}\] ◻
Thanks to the range of possible values of \(L\) in Lemma 22, we conclude that \(\mu_{\mathbb{T}} (f)\) has a limit \(\mu(f)\) (as the size \(L\) of \({\mathbb{T}}\) diverges) and \(P_t^{(N)} f\) converges exponentially fast to it in \(\ell^\infty\) at rate \(N^d/(4d^d)\). In particular, the infinite volume process is ergodic, there is a unique invariant measure and therefore a unique Gibbs measure (see Liggett05?*Theorem IV.2.15). Moreover, this exponential ergodicity implies (see the remark at the end of Martinelli99?*Section 3.5) \[\label{eq:trelN:upper} T_{\mathrm{rel}}^{(N)}\leqslant 4d^d/N^d.\tag{55}\]
Furthermore, for any local function \(f\colon\Omega\to{\mathbb{R}}\) we have \[\begin{align} \label{eq:trel:small:scales}\mathop{\mathrm{Var}}(f)&{}\overset{\eqref{eq:def:trel}}\leqslant T_{\mathrm{rel}}^{(N)}\mathcal{D}^{(N)}(f)\overset{\eqref{eq:def:DN}}=T_{\mathrm{rel}}^{(N)}\sum_{x\in{\mathbb{Z}} ^d}\mu\left(\mathop{\mathrm{Var}}_{x-\Lambda_{d,N}}(f)\right)\\ \nonumber&{}\overset{\eqref{eq:def:trel}}\leqslant T_{\mathrm{rel}}^{(N)}\sum_{x\in{\mathbb{Z}} ^d}\sup_{\tau\in\Omega_{{\mathbb{Z}} ^d\setminus(-\Lambda_{d,N})}}T_{\mathrm{rel}}^{-\Lambda_{d,N},\tau}\mu\left(\mathcal{D}_{x-\Lambda_{d,N}}(f)\right)\\ \nonumber&{}\overset{\eqref{eq:def:D}}=T_{\mathrm{rel}}^{(N)}\sup_{\tau}T_{\mathrm{rel}}^{-\Lambda_{d,N},\tau}\sum_{x\in{\mathbb{Z}} ^d}\mu\left(\sum_{y\in x+\Lambda_{d,N}}\mu_{x-\Lambda_{d,N}}(\mathop{\mathrm{Var}}_y(f))\right)\\ &{}\overset{\eqref{eq:def:D}}=T_{\mathrm{rel}}^{(N)}\sup_{\tau}T_{\mathrm{rel}}^{-\Lambda_{d,N},\tau}|\Lambda_{d,N}|\mathcal{D}(f). \end{align}\tag{56}\]
Using 5 , 56 , 55 and Lemma 19, we get \[T_{\mathrm{rel}}\leqslant\frac{4d^d}{N^d}(N+1)^de^{2(\beta+1)(N+4)^d}\leqslant e^{(\beta+1)m^{C'}}\leqslant Ce^{e^{C\beta}},\] recalling 49 and \(m=\lceil 1/(1-\tanh(\beta))\rceil\geqslant 2\) and taking suitably large \(C',C>0\) depending only on \(d\) (and \(a_d\)). ◻
We saw in Section 1 how modifying switches to affect lamps can be understood as a game, whose goal is to turn off all lamps by an appropriate sequence of switches. Combinatorial group theory, insofar as it is concerned with computations in group presentations, generalises broadly this idea; see in particular the beautiful application by Conway and Lagarias to tiling problems Conway90?. We start by reviewing the basics.
Consider a group presentation \(G=\langle S\mid R\rangle\). This is a compact description of a group \(G\) with the property that every element may be represented as a word over \(S\cup S^{-1}\), and such that two such representations define the same group element if one word may be obtained from the other by repeated insertion of an word from \(R\cup R^{-1}\), and free cancellation \((s\cdot s^{-1}=s^{-1}\cdot s=1)\). For example, \(G=\langle x,y\mid x^{-1} y^{-1} x y\rangle\) defines the group \({\mathbb{Z}} ^2\); every element may be uniquely represented as a word of the form \(x^m y^n\) with \(m,n\in{\mathbb{Z}}\).
Consider a word \(w\) over \(S\cup S^{-1}\) that represents the identity in \(G\). The process of reducing \(w\) to \(1\) by repeated insertions of relators may be visualised in a van Kampen diagram: this is a planar graph all of whose edges are oriented and labeled by \(S\), all of whose faces read an element of \(R\cup R^{-1}\) on their boundary, if edges in reversed orientation are interpreted as \(S^{-1}\); and whose external boundary reads \(w\). The weight of \(w\) is the minimal number of faces of a van Kampen diagram for \(w\), or equivalently the minimal \(n\) such that \(w\) may be written in the free group \(F_S\) in the form \[\label{eq:Dehn} w=\prod_{i=1}^n u_i r_i u_i^{-1},\quad r_i\in R\cup R^{-1}.\tag{57}\] The Dehn function of \(G\) is the function \(\Delta(n)\) counting the maximal weight of all words of length at most \(n\) that represent the identity. For example, the Dehn function of \(\mathbb{Z}^2\), with its presentation given above, is \(\Delta(n)=\lfloor (n/4)^2\rfloor\); the weight of \(x^{-n}y^{-n}x^n y^n\) is \(n^2\) and maximises the weight among words of length at most \(4n\).
The Dehn function of a finitely presented group depends only mildly on the presentation, and is a crude measure of the complexity of the word problem (namely, the problem of testing whether two words define equal elements of the group). For example, some of the most computationally tractable groups, Gromov’s word hyperbolic groups, may be characterised by the fact that their Dehn function is bounded by a linear function. There exist finitely presented groups with unsolvable word problem, and therefore with uncomputable Dehn function.
There is nothing sacred about the Dehn function, and in fact many other invariants have been considered; see Gromov93?*§5 for such a collection. An interesting invariant defined there is the filling length, also considered by Gersten in Gersten95?. The filling length of a van Kampen diagram \(\mathcal{K}\) with basepoint \(*\) is the minimal \(\lambda\) such that there exists a sequence \(\gamma_0,\dots,\gamma_n\) of loops at \(*\) in \(\mathcal{K}\), such that
all \(\gamma_i\) have length at most \(\lambda\);
\(\gamma_0\) is the trivial loop at \(*\), and \(\gamma_n\) is the boundary loop of \(\mathcal{K}\);
\(\gamma_i\) and \(\gamma_{i+1}\) differ by exactly one cell, which is outside \(\gamma_i\) and inside \(\gamma_{i+1}\) with winding number \(1\)
(such a sequence of loops is called a combinatorial null-homotopy). As before, one defines the filling length of a word \(w\) representing the identity as the minimal filling length of van Kampen diagrams for \(w\), and the filling function \(\lambda(n)\) as the maximum of \(\lambda(w)\) over all words \(w\) of length at most \(n\) that represent the identity.
We may refine it as follows: firstly, for a generator \(s\), we let \(|w|_s\) denote the number of \(s^{\pm1}\) in a freely reduced word \(w\); then, for a van Kampen diagram \(\mathcal{K}\), we let \(\lambda_s(\mathcal{K})\) be the minimal number \(\lambda\) such that there exists a combinatorial null-homotopy \((\gamma_0,\dots,\gamma_n)\) with \(|\gamma_i|_s\leqslant\lambda\) for all \(i\), and as before \(\lambda_s(w)\) as the minimum of \(\lambda_s(\mathcal{K})\) over van Kampen diagrams for \(w\), and \(\lambda_s(n)\) as the maximum of \(\lambda_s(w)\) over all words \(w\) of length at most \(n\) that represent the identity.
Secondly, we denote by \(\gamma_i^\circ\) the part of \(\gamma_i\) that is not on the boundary of the van Kampen diagram (a “chord”). For example, van Kampen diagrams in groups admitting a free subgroup of finite index are thin, so they admit null homotopies of bounded length. We let \(\lambda^\circ_s(\mathcal{K})\) be the minimum, over combinatorial null-homotopies, of \(\max_i|\gamma_i^\circ|_s\), and define \(\lambda^\circ_s(w)\) and \(\lambda^\circ_s(n)\) as above.
We return to the finitely presented group \(\Gamma\) of 7 . Kassabov and Riley prove in Kassabov12? that its Dehn function satisfies \(\Delta(n)\precsim n^4\), and note that the actual value \(\Delta(n)\approx n^2\) follows from Bartholdi08?, Drutu04?.
To better understand van Kampen diagrams for \(\Gamma\), we note that every word \(w\) over \(\{a,x^{\pm1},y^{\pm1}\}\) that represents the identity has an associated lamp configuration: if \(w=w_1\dots w_n\), start at the origin and read the letters \(w_1,\dots,w_n\) in sequence, moving left/right/up/down when the letter is \(x^{-1},x,y,y^{-1}\) respectively, and toggling the lamp at the current position when the letter is \(a^{\pm1}\). After the last letter has been read, we have returned to the origin. The configuration diameter of \(w\) is the diameter of the associated lamp configuration, namely, the maximal distance between two ON lamps, in the \(\ell^\infty\) metric from 8 .
For example, the word \(a x^{-2^n} a x^{2^n}y^{-2^n} a y^{2^n}\) represents the identity in \(\Gamma\), and has configuration diameter \(2^n\) with three ON lamps. Figure 7 shows a van Kampen diagram for this word with \(n=3\).
We obtain the following consequence of Theorem 4:
Corollary 2. In the group \(\Gamma\), the \(a\)-filling length function \(\lambda^\circ_a(n)\) is bounded from below by \(\log_2(n)\).
Proof. Consider the word \(w=a x^{-2^n} a x^{2^n}y^{-2^n} a y^{2^n}\) mentioned above, which has size \(\mathcal{O}(2^n)\) and represents the identity in \(\Gamma\). Let \(\mathcal{K}\) be a van Kampen diagram for \(w_n\), and let \((\gamma_0,\dots,\gamma_\Delta)\) be a combinatorial null-homotopy: a sequence of loops in \(\mathcal{K}\), based at the start of \(w\), such that \(\gamma_0\) is the empty loop, \(\gamma_\Delta\) in the peripheral loop reading \(w\), and each \(\gamma_i\) differs from \(\gamma_{i-1}\) by enclosing a single extra cell \(C_i\). If \(u_i\) is the word read along \(\gamma_i\) from the basepoint to the beginning of the relator \(r_i\) around \(C_i\), we obtain a corresponding expression of \(w_n\) as a product of conjugates of relators \(r=a x a x^{-1} y a y^{-1}\) and \(s=[x,y]\), \[w_n=\prod_{i=1}^\Delta u_i r_i u_i^{-1},\quad r_i\in\{r^{\pm1},s^{\pm1}\}.\] Consider now a loop \(\gamma_i\); the word read along it also represents the identity in \(\Gamma\). Consider the associated lamp configuration: follow the \(x^{\pm1},y^{\pm1}\) letters of \(\gamma_i\), and flip the lamp at the current position when an \(a\) is read. This lamp configuration \(W_i\) is admissible, and the configurations \((W_0,\dots,W_\Delta)\) form a chain. The claim now follows from Theorem 4. ◻
In an entirely similar manner, we deduce from Theorem 8 the following result.
Corollary 3. There exist \(\alpha_2>0\) and \(n_2\in\mathbb{N}\) such that, for any \(n\geqslant n_2\), there exists \(w\in\{a,x^{\pm1},y^{\pm1}\}\) with the following properties:
it represents the identity in \(\Gamma\);
it has configuration diameter at most \(n\);
it contains \(L\) letters \(a\) or \(a^{-1}\);
its \(a\)-filling length function \(\lambda_a(w)\) is at least \(L+n^{\alpha_2}\).0◻
We note in passing that the homotopy in the van Kampen diagram contains strictly more information than the chain from Definition 1, since it specifies not only which lamps remain ON at an intermediate step of a turning-off, but also in which order and through which lamp switches they are visited. A tight reformulation of Theorem 4 in terms of group theory would be homological, rather than homotopical. In particular, a null-homotopy determines a chain, but not conversely.
We start with a brief introduction to the Ledrappier subshift from an ergodic theory perspective. Consider a measure-preserving action of a discrete group \(G\) on a probability space \((X,\mu)\), and recall that it is called \(r\)-mixing if for any measurable \(A_1,\dots,A_r\) one has \[\mu\left(\bigcap_{i=1}^r g_i\cdot A_i\right)\to\prod_{i=1}^r\mu(A_i)\text{ as }d(g_i,g_j)\to\infty\;\forall i\ne j;\] and “\(2\)-mixing” is simply called “mixing”. The question whether, for \({\mathbb{Z}}\)-actions, mixing necessarily implies \(r\)-mixing for all \(r\) is an old open problem in ergodic theory. Ledrappier gave in Ledrappier78? an example of mixing \({\mathbb{Z}} ^2\)-dynamical system which is not \(3\)-mixing — this is precisely the system of lamp configurations that are cycles, in our terminology; namely, \[X=\left\{x\colon{\mathbb{Z}} ^2\to{\mathbb{F}} _2\middle| x(m,n)+x(m+1,n)+x(m,n+1)=0\;\forall m,n\right\},\] with \(\mu\) its Haar measure. (To see that it is not \(3\)-mixing, consider \(A_1=A_2=A_3=[x(0,0)=1]\), and \(g_1=(0,0),g_2=(2^n,0),g_3=(0,2^n)\). Each of the events \(g_i\cdot A_i\) has measure \(1/2\), but their intersection is empty.)
In Arenas-Carmona08?, Arenas-Carmona, Berend and Bergelson show that the Ledrappier system is “almost \(r\)-mixing” for all \(r\geqslant 2\). To explain that notion, let \(d_{\mathrm{H}}(C,D)\) denote the Hausdorff distance between two elements or subsets \(C,D\) of \({\mathbb{Z}} ^2\), and let \(d_{\mathrm{H}}(C,\mathscr D)\) denote the minimal distance between \(C\subseteq{\mathbb{Z}} ^2\) and any element of a collection \(\mathscr D\) of subsets of \({\mathbb{Z}} ^2\). For \(r\geqslant 2\) and \(\mathscr M\subseteq G^r\), they define \(X\) to be \(r\)-mixing modulo \(\mathscr M\) if for any measurable \(A_1,\dots,A_r\) one has \[\mu\left(\bigcap_{i=1}^r g_i\cdot A_i\right)\to\prod_{i=1}^r\mu(A_i)\text{ as }d(g_i,g_j)\to\infty\;\forall i\ne j\text{ and }d_{\mathrm{H}}(\{g_1,\dots,g_r\},\mathscr M)\to\infty.\] Recall from Definition 7 that admissible configurations are finite sums of translates of the plaquette. In Arenas-Carmona08?’s terminology, an admissible configuration with \(r\) lamps is called a special \(r\)-gon. Their main result is that, for all \(r\geqslant 3\), the Ledrappier subshift is \(r\)-mixing modulo special \(r\)-gons. In particular, the Ledrappier shift is \(3\)-mixing modulo large plaquettes.
This led to a study of the structure of special \(r\)-gons, in Arenas-Carmona08?*§7. The results in Section 2.2.3 partly answer their questions, as we now show. We first recall:
Theorem 11 (Arenas-Carmona08?*Theorem 7.1). If \(X\subset \mathbb{Z}^{2}\) is an admissible set of cardinality \(n\), then \(X\) is a sum modulo \(2\) of at most \(n^3\) large plaquettes.
Let us concentrate on \(d=2\). Let \(h(n)\) denote the minimal number of large plaquettes sufficient to represent any admissible set of cardinality at most \(n\); so \(h(n)\leqslant n^3\). In Arenas-Carmona08?*Remark 7.2, the authors ask for lower bounds on \(h(n)\) and construct examples which, conjecturally, give \(h(n)\geqslant 3n-8\). Our construction proves a stronger lower bound:
Theorem 12 (The function \(h(n)\) grows superlinearly). For some \(\eta>0\) and any \(n\) large enough, we have \[\label{eq:h} h(n)\geqslant n^{1+\eta}.\tag{58}\]
Proof. By Corollary 1, there exists a nonempty admissible configuration \[X=\{x_1,x_2,\dots,x_m\}\] of at most \(n\) points such that any two distinct points have \(\rho\)-distance at least \(r\), where \(r\geqslant c'\ln(n-2)\geqslant 15+c\ln n\) for some positive constants \(c,c'\), as long as \(n\) is large enough. Up to repeatedly replacing \(X\) by \(X\sqcup(X+v)\) for some \(v\in\mathbb{Z}^2\) at \(\rho\)-distance at least \(r\) from \(X-X\), we may assume \(m\geqslant n/2\).
Since \(X\) is admissible, write it as the sum of at most \(h(n)\) large plaquettes. For every \(i=1,\dots,m\), let \(M_i\) be the set of those large plaquettes in the sum that have at least one vertex at \(\rho\)-distance \(<r/3\) from \(x_i\). If \(r\geqslant 3\), then the sets \(M_i\) are pairwise disjoint.
Let \(X_i\) be the sum modulo \(2\) of all large plaquettes from \(M_i\). All other points of \(X_i\) are at \(\rho\)-distance at least \(r/3\) from \(x_i\). Applying Lemma 7 to \(X_i\) and the point \(x_i\), we obtain \[|X_i|\geqslant\frac{3}{4}\left(\frac{3}{2}\right)^{r/3} \geqslant 4 n^{\eta}, \qquad \eta\mathrel{\vcenter{:}}=\frac{c}{3}\ln\left(\frac{3}{2}\right)>0.\] This implies \[|M_i|\geqslant 2 n^{\eta}.\] Summing over \(i=1,\dots,m\) and using the disjointness of the sets \(M_i\) gives \(h(n)\geqslant 2m\,n^\eta\geqslant n^{1+\eta}\), as claimed. ◻
L.B.and I.M.gratefully acknowledge support from the SNSF Advanced Grant TMAG-2_216487/1.
We thank Paul Chleboun, Fabio Martinelli and Marius Tiba for helpful discussions, as well as Guillaume Aubrun for sparking this collaboration. For open access purposes, the authors have applied a CC BY public copyright license to any author-accepted manuscript version arising from this submission.