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

Tightness of exponential metrics for log-correlated Gaussian fields in arbitrary dimension

Jian Ding  Ewain Gwynne  Zijie Zhuang Peking UniversityUniversity of ChicagoUniversity of Pennsylvania
(March 2024)
Abstract

We prove the tightness of a natural approximation scheme for an analog of the Liouville quantum gravity metric on d\mathbb{R}^{d} for arbitrary d2d\geq 2. More precisely, let {hn}n1\{h_{n}\}_{n\geq 1} be a suitable sequence of Gaussian random functions which approximates a log-correlated Gaussian field on d\mathbb{R}^{d}. Consider the family of random metrics on d\mathbb{R}^{d} obtained by weighting the lengths of paths by eξhne^{\xi h_{n}}, where ξ>0\xi>0 is a parameter. We prove that if ξ\xi belongs to the subcritical phase (which is defined by the condition that the distance exponent Q(ξ)Q(\xi) is greater than 2d\sqrt{2d}), then after appropriate re-scaling, these metrics are tight and that every subsequential limit is a metric on d\mathbb{R}^{d} which induces the Euclidean topology. We include a substantial list of open problems.


Acknowledgments. We thank Karim Adiprasito, Timothy Budd, Hugo Falconet, Josh Pfeffer, Scott Sheffield, and Xin Sun for helpful discussions. J.D. is partially supported by NSFC Key Program Project No. 12231002. E.G. was partially supported by a Clay research fellowship and by NSF grant DMS-2245832. Z.Z. was partially supported by NSF grant DMS-1953848.

1 Introduction

There has been an enormous amount of research in the past several decades concerning random geometry in two dimensions. Some of the major topics in this subject include Schramm-Loewner evolution, conformal field theory, statistical mechanics models on planar lattices, random planar maps, Liouville quantum gravity, and random geometries related to the KPZ universality class. We will not attempt to survey this vast literature here, but see, e.g., [sheffield-icm, gwynne-ams-survey, bp-lqg-notes, bn-sle-notes, pw-gff-notes, ghs-mating-survey, vargas-dozz-notes, legall-sphere-survey, ganguly-dl-survey] for some recent expository articles. However, most of the results in this area have not been extended to higher dimensions. One reason for this is that conformal invariance (or covariance) plays a central role in many of the results in two dimensions, and there are no non-trivial conformal maps in higher dimensions. Another reason is that many of the arguments in the two-dimensional case rely on topological properties which are not true in higher dimensions, e.g., the Jordan curve theorem.

In this paper, we consider the problem of constructing an analog of the Liouville quantum gravity (LQG) metric on d\mathbb{R}^{d}, for arbitrary d2d\geq 2. Heuristically speaking, LQG is the random geometry described by the random Riemannian metric tensor

eγh(dx2+dy2)e^{\gamma h}(dx^{2}+dy^{2}) (1.1)

where γ(0,2]\gamma\in(0,2] is a parameter, dx2+dy2dx^{2}+dy^{2} is the Euclidean metric tensor, and hh is a variant of the Gaussian free field (GFF) on 2\mathbb{R}^{2} (or more generally on a Riemann surface). See, e.g., [shef-gff, bp-lqg-notes, pw-gff-notes] for an introduction to the GFF. The definition (1.1) does not make literal sense since hh is a generalized function (distribution) instead of a true function, so its exponential cannot be defined pointwise. Nevertheless, one can define various objects associated with (1.1) by replacing hh with a sequence of continuous functions which approximate hh, then taking an appropriate limit.

Perhaps the easiest object to construct in this way is the LQG area measure, which is a limit of regularized versions of eγhdxdye^{\gamma h}\,dx\,dy (where dxdydx\,dy denotes Lebesgue measure). The construction of this measure is a special case of the theory of Gaussian multiplicative chaos (GMC), which allows one to make sense of random measures of the form eαh(x)dσ(x)e^{\alpha h(x)}\,d\sigma(x) for α>0\alpha>0, whenever hh is a log-correlated Gaussian field on a domain UdU\subset\mathbb{R}^{d} (for arbitrary d1d\geq 1) and σ\sigma is an appropriate deterministic base measure on UU. See [shef-kpz, rhodes-vargas-review, bp-lqg-notes] for more on Gaussian multiplicative chaos and the LQG area measure.

Recent works have also constructed the Riemannian distance function associated with (1.1), i.e., the LQG metric. This is a random metric DhD_{h} on 2\mathbb{R}^{2} constructed as follows. For ε>0\varepsilon>0, let hεh_{\varepsilon} be the convolution of the Gaussian free field with the heat kernel pε2/2(z)=1πε2e|z|2/ε2p_{\varepsilon^{2}/2}(z)=\frac{1}{\pi\varepsilon^{2}}e^{-|z|^{2}/\varepsilon^{2}}. Also let ξ=ξ(γ)=γ/dγ\xi=\xi(\gamma)=\gamma/d_{\gamma}, where dγd_{\gamma} is the so-called LQG dimension exponent [dg-lqg-dim]. Then, let

Dhε(z,w):=infP:zw01eξhε(P(t))|P(t)|dt,z,w2,D_{h}^{\varepsilon}(z,w):=\inf_{P:z\to w}\int_{0}^{1}e^{\xi h_{\varepsilon}(P(t))}|P^{\prime}(t)|\,dt,\quad\forall z,w\in\mathbb{R}^{2}, (1.2)

where the infimum is over all piecewise continuously differentiable paths P:[0,1]2P:[0,1]\to\mathbbm{R}^{2} from zz to ww. The papers [dddf-lfpp, gm-uniqueness] prove that there exist normalizing constants {𝔞ε}ε>0\{\mathfrak{a}_{\varepsilon}\}_{\varepsilon>0} such that 𝔞ε1Dhε\mathfrak{a}_{\varepsilon}^{-1}D_{h}^{\varepsilon} converges in probability to a limiting metric with respect to the topology of uniform convergence on compact subsets of 2×2\mathbb{R}^{2}\times\mathbb{R}^{2} (the convergence in probability was recently improved to a.s. convergence in [devlin-lfpp-as]). In particular, it was shown in [dddf-lfpp] that the approximating metrics are tight, and in [gm-uniqueness] (building on [local-metrics, gm-confluence, lqg-metric-estimates]) that the subsequential limit is unique. The proofs in these papers are much more difficult than the proofs in the construction of the LQG area measure. Intuitively, this is because the minimizing path in (1.2) depends on ε\varepsilon. See [ddg-metric-survey] for a survey of known results about the LQG metric.

In light of the theory of Gaussian multiplicative chaos, it is natural to wonder whether there is an analogous theory of exponential metrics associated with log-correlated Gaussian fields on d\mathbb{R}^{d} for arbitrary111Note that when d=1d=1, the metric induced by eξhe^{\xi h} is simply given by the one-dimensional GMC measure, as any path in \mathbb{R} is an interval. d2d\geq 2, which generalizes the LQG metric. The construction of such a theory is listed as Problem 7.19 in [gm-uniqueness].

This paper carries out the first major step toward such a theory: namely, we prove the tightness of a natural approximation scheme similar to (1.2) for log-correlated Gaussian fields on d\mathbb{R}^{d} (in the full subcritical phase of ξ\xi values). That is, we carry out the higher-dimensional analog of [dddf-lfpp]. See Theorem 1.2 below for a precise statement. We expect that it will be challenging, but possible to prove that the subsequential limit is unique (and characterized by a list of axioms similar to the ones which characterize the LQG metric in dimension two [gm-uniqueness]) by adapting the arguments in the two-dimensional case [gm-uniqueness, dg-uniqueness]. Indeed, these arguments do not use two-dimensionality in as fundamental a way as the proof of tightness in [dddf-lfpp]. See Problem LABEL:prob:uniqueness for further discussion.

More speculatively, our limiting metric might have connections to other higher-dimensional extensions of objects related to LQG, e.g., Liouville conformal field theory in even dimensions [cercle-higher-dimension, dhks-even-dim], the higher-dimensional analogs of the Brownian map considered in [ml-iterated-folding], uniform samples from various classes of triangulations of higher-dimensional spheres (see, e.g., [bz-locally-constructible, dj-3-manifolds]), higher-dimensional analogs of random planar maps constructed from trees [BC23, budd-lionni-3-spheres], and random graphs in d\mathbb{R}^{d} arising from sphere packings (see, e.g., [bc-sphere-packing]). See Subsection 1.2 for more details.

The problem of constructing natural random Riemannian metrics in dimension d3d\geq 3 is also of substantial interest in theoretical physics in the context of quantum gravity (see, e.g., the books [gh-quantum-gravity, adj-quantum-geometry, rovelli-quantum-geometry]). We refer to the introductions of [BC23, budd-lionni-3-spheres] for additional relevant discussion and references.

The proofs in this paper are by necessity substantially different than those in the two-dimensional case [dddf-lfpp]. In particular, we do not have an a priori Russo-Seymour-Welsh (RSW) type estimate (which in the two-dimensional case comes from a conformal invariance argument), and various path-joining arguments in [dddf-lfpp] do not work in higher dimensions. For these reasons, we use a fundamentally novel approach to proving tightness which bypasses any direct proofs of RSW estimates as well as the use of the Efron-Stein inequality. See Subsection 1.3 for details.

The results of this paper open up a number of interesting questions about random metrics on d\mathbb{R}^{d}. See Section LABEL:sec:open-problem for a discussion of some open problems.

1.1 Definitions and main result

We now introduce some notation and state the main result of this paper. We consider the space d\mathbb{R}^{d} with d2d\geq 2 and define the box

Br(x):=x+(r,r)d,xd,r>0.B_{r}(x):=x+(-r,r)^{d},\quad\forall x\in\mathbb{R}^{d},\quad\forall r>0. (1.3)

Fix a smooth function K:d[0,){\hyperref@@ii[K-condition1]{\mathfrak{K}}}:\mathbb{R}^{d}\rightarrow[0,\infty) and r0>0{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}>0 such that

  1. 1.

    K is radially symmetric, meaning that K(x)=K(y){\hyperref@@ii[K-condition1]{\mathfrak{K}}}(x)={\hyperref@@ii[K-condition1]{\mathfrak{K}}}(y) for any x,ydx,y\in\mathbb{R}^{d} with the same Euclidean norm.

  2. 2.

    K is supported222 We expect that our arguments can be adapted to the case when K is not compactly supported but has sufficiently rapid decay at \infty. This would require some added technicalities similar to ones encountered in [dddf-lfpp]. However, the choice of K in this paper is in some sense unimportant since, regardless of the choice of K, the fields we consider are closely related to the canonical log-correlated Gaussian field on d\mathbbm{R}^{d} considered in [lgf-survey, fgf-survey] (see Remark 1.4). in the box Br0(0)B_{{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}}(0).

  3. 3.

    K is normalized such that dK(x)2dx=1\int_{\mathbb{R}^{d}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}(x)^{2}dx=1.

We also let WW be a space-time white noise on d\mathbbm{R}^{d}. That is, WW is the Gaussian random generalized function on d×(0,)\mathbbm{R}^{d}\times(0,\infty) such that for any fL2(d×(0,))f\in L^{2}(\mathbbm{R}^{d}\times(0,\infty)), the formal integral d0f(y,t)W(dy,dt)\int_{\mathbbm{R}^{d}}\int_{0}^{\infty}f(y,t)W(dy,dt) is centered Gaussian with variance fL22\|f\|_{L^{2}}^{2}.

We consider a log-correlated Gaussian field hh and its approximation hnh_{n}, defined as follows:

h(x)\displaystyle h(x) =d01K(yxt)td+12W(dy,dt)and\displaystyle=\int_{\mathbb{R}^{d}}\int_{0}^{1}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}\big{(}\frac{y-x}{t}\big{)}t^{-\frac{d+1}{2}}W(dy,dt)\quad\mbox{and} (1.4)
hn(x)\displaystyle h_{n}(x) =d2n1K(yxt)td+12W(dy,dt)\displaystyle=\int_{\mathbb{R}^{d}}\int_{2^{-n}}^{1}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}\big{(}\frac{y-x}{t}\big{)}t^{-\frac{d+1}{2}}W(dy,dt)

for xdx\in\mathbb{R}^{d} and integer n1n\geq 1. From the definition of WW, we see that hh and hnh_{n} are centered Gaussian processes with covariance kernels

Cov(h(x1),h(x2))\displaystyle\mathrm{Cov}(h(x_{1}),h(x_{2})) =011t(KK)(x1x2t)dt\displaystyle=\int_{0}^{1}\frac{1}{t}({\hyperref@@ii[K-condition1]{\mathfrak{K}}}*{\hyperref@@ii[K-condition1]{\mathfrak{K}}})\left(\frac{x_{1}-x_{2}}{t}\right)\,dt
Cov(hn(x1),hn(x2))\displaystyle\mathrm{Cov}(h_{n}(x_{1}),h_{n}(x_{2})) =2n11t(KK)(x1x2t)dt,\displaystyle=\int_{2^{-n}}^{1}\frac{1}{t}({\hyperref@@ii[K-condition1]{\mathfrak{K}}}*{\hyperref@@ii[K-condition1]{\mathfrak{K}}})\left(\frac{x_{1}-x_{2}}{t}\right)\,dt\,, (1.5)

where KK{\hyperref@@ii[K-condition1]{\mathfrak{K}}}*{\hyperref@@ii[K-condition1]{\mathfrak{K}}} denotes the convolution. Using the representation (1.4) and the fact that WW is a random tempered distribution (see e.g. Section 2.3 of [fgf-survey]), one can verify that each hnh_{n} has a modification which is a smooth function (see also Proposition 2.1 of [df-lqg-metric]). We henceforth assume that each hnh_{n} has been replaced by such a modification. Furthermore, from (1.1) we get Varhn(x)=nlog2\operatorname{Var}h_{n}(x)=n\log 2 for each xdx\in\mathbb{R}^{d}. The process hh is interpreted as a generalized function, and is closely related to the log-correlated Gaussian field on d\mathbb{R}^{d} considered in [lgf-survey, fgf-survey] (see Remark 1.4).

Analogously333As explained in [dddf-lfpp] (see also [cg-support-thm, Section 3.1]), in the two-dimensional case, the convolution of the planar Gaussian free field with the heat kernel (at an appropriate nn-dependent time) has the same law as the field hnh_{n} of (1.4) with K(x)=2πe|x|2{\hyperref@@ii[K-condition1]{\mathfrak{K}}}(x)=\sqrt{\frac{2}{\pi}}e^{-|x|^{2}}, up to adding a random continuous function. Hence (1.6) is directly analogous to (1.2). To avoid unnecessary technical work, in this paper we require that K is compactly supported, but we expect that our results can be fairly easily extended to the case where K is not compactly supported but has sufficiently fast decay at \infty. to (1.2), for a parameter ξ>0\xi>0, we define the exponential metric associated with hnh_{n} as follows:

Dn(z,w):=infP:zw01eξhn(P(t))|P(t)|dt,z,wd,D_{n}(z,w):=\inf_{P:z\to w}\int_{0}^{1}e^{\xi h_{n}(P(t))}|P^{\prime}(t)|dt\,,\quad\forall z,w\in\mathbb{R}^{d}, (1.6)

where the infimum is taken over all piecewise continuously differentiable paths P:[0,1]dP:[0,1]\rightarrow\mathbb{R}^{d} joining z,wz,w. This can be interpreted as an approximation of the random metric formally given by reweighting the Euclidean lengths of paths by eξhe^{\xi h}. We will be interested in (subsequential) limits of the renormalized metrics λn1Dn\lambda_{n}^{-1}D_{n}, where the normalizing constant λn\lambda_{n}444For technical reasons, we first work with this particular choice of normalizing constant. However, in the end, we can choose any reasonable normalizing constant, such as the median of Dn(0,e1)D_{n}(0,e_{1}) or Dn(B1(0),B2(0))D_{n}(\partial B_{1}(0),\partial B_{2}(0)). is defined as:

λn:=median of Dn(0,e1;B2(0)),\lambda_{n}:=\mbox{median of }D_{n}(0,e_{1};B_{2}(0))\,, (1.7)

where Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)) denotes the minimal DnD_{n}-length of a path joining 0 and e1:=(1,0,,0)e_{1}:=(1,0,\ldots,0) inside the box B2(0)B_{2}(0).

In Section 3, we will prove the following.

Proposition 1.1.

For each ξ>0\xi>0, there exists Q=Q(ξ)Q=Q(\xi)\in\mathbbm{R} such that

λn=2(1ξQ)n+o(n)as n.\lambda_{n}=2^{-(1-\xi Q)n+o(n)}\quad\mbox{as }n\rightarrow\infty\,. (1.8)

Furthermore, ξQ(ξ)\xi\mapsto Q(\xi) is a continuous, non-increasing function and we have

1ξ2dQ(ξ)1ξ+2,ξ>0.\frac{1}{\xi}-\sqrt{2d}\leq Q(\xi)\leq\frac{1}{\xi}+\sqrt{2}\,,\quad\forall\xi>0. (1.9)

The proof of Proposition 1.1 is via a subadditivity argument. Just like in the two-dimensional case, we do not know the value of Q(ξ)Q(\xi) explicitly (see Problems LABEL:prob:positive-Q and LABEL:prob:special). Analogously to the two-dimensional case (see [dg-supercritical-lfpp, Equation (1.4)]), we define the critical value

ξcrit:=sup{ξ>0:Q(ξ)>2d}.\xi_{\mathrm{crit}}:=\sup\left\{\xi>0:Q(\xi)>\sqrt{2d}\right\}. (1.10)

See Remark 1.3 for some discussion of why this value is critical. We note that ξ<ξcrit\xi<\xi_{\mathrm{crit}} if and only if Q(ξ)>2dQ(\xi)>\sqrt{2d}. The lower bound Q(ξ)1ξ2dQ(\xi)\geq\frac{1}{\xi}-\sqrt{2d} from (1.9) implies that ξcrit122d\xi_{\mathrm{crit}}\geq\frac{1}{2\sqrt{2d}}, and the upper bound Q(ξ)1ξ+2Q(\xi)\leq\frac{1}{\xi}+\sqrt{2} implies that ξcrit12d2<\xi_{\rm crit}\leq\frac{1}{\sqrt{2d}-\sqrt{2}}<\infty. The main result of this paper is the tightness of our approximating metrics in the full subcritical phase.

Theorem 1.2.

When ξ<ξcrit\xi<\xi_{\mathrm{crit}}, equivalently Q(ξ)>2dQ(\xi)>\sqrt{2d}, the sequence of metrics {λn1Dn(,)}n1\{\lambda_{n}^{-1}D_{n}(\cdot,\cdot)\}_{n\geq 1} is tight with respect to the topology of uniform convergence on compact subsets of d×d\mathbb{R}^{d}\times\mathbb{R}^{d}. Furthermore, each possible subsequential limit (in distribution) is a metric on d\mathbb{R}^{d} which induces the Euclidean topology.

Remark 1.3.

When Q(ξ)<2dQ(\xi)<\sqrt{2d}, we expect that the metrics λn1Dn\lambda_{n}^{-1}D_{n} are not tight with respect to the topology of uniform convergence on compact subsets of d×d\mathbb{R}^{d}\times\mathbb{R}^{d}. So, our result should be optimal modulo the critical case when Q(ξ)=2dQ(\xi)=\sqrt{2d}. Indeed, the maximum of hnh_{n} on a fixed bounded open set UdU\subset\mathbb{R}^{d} should grow like (2d+o(1))(log2)n(\sqrt{2d}+o(1))(\log 2)n as nn\to\infty, see e.g. [Mad15]. From this and the continuity properties of hnh_{n} (Claim (2) of Lemma 2.3), if ε>0\varepsilon>0 is fixed, then when nn is large with high probability there exists zUz\in U such that hn(w)(2dε)(log2)nh_{n}(w)\geq(\sqrt{2d}-\varepsilon)(\log 2)n for each ww in the box B2n(z)B_{2^{-n}}(z). For this choice of zz, the definition of DnD_{n} shows that the DnD_{n}-distance from zz to B2n(z)\partial B_{2^{-n}}(z) is at least 2[(2dε)ξ1]n2^{[(\sqrt{2d}-\varepsilon)\xi-1]n}. By (1.8),

λn1Dn(z,B2n(z))2(2dQε+o(1))ξn.\lambda_{n}^{-1}D_{n}(z,\partial B_{2^{-n}}(z))\geq 2^{(\sqrt{2d}-Q-\varepsilon+o(1))\xi n}.

If Q(ξ)<2dQ(\xi)<\sqrt{2d}, then for a small enough choice of ε\varepsilon, this goes to \infty as nn\to\infty, which means that λn1Dn\lambda_{n}^{-1}D_{n} cannot be tight with respect to the local uniform topology.

In the two-dimensional case, it was shown in [dg-supercritical-lfpp, dg-uniqueness] that the re-scaled approximating metrics converge with respect to the topology on lower semicontinuous functions for all ξ>0\xi>0 (including when Q(ξ)2Q(\xi)\leq 2). However, when Q(ξ)<2Q(\xi)<2, the limiting metric does not induce the Euclidean topology on 2\mathbb{R}^{2}. Rather, there are uncountably many “singular points” which lie at infinite distance from every other point. It is plausible that similar statements are true for general d2d\geq 2, but we do not address this in the present paper. See Problem LABEL:prob:supercritical.

Remark 1.4.

The field hh of (1.4) is closely related to the log-correlated Gaussian field on d\mathbb{R}^{d} considered in [lgf-survey, fgf-survey]. Indeed, define the random generalized function hh^{\infty} in the same manner as hh in (1.4), but with tt integrated over (0,)(0,\infty) instead of over (0,1)(0,1). Then, a short computation shows that for any choice of the kernel K above, one can make sense of hh^{\infty} as a random generalized function viewed modulo additive constant555That is, dh(x)g(x)dx\int_{\mathbb{R}^{d}}h^{\infty}(x)g(x)\,dx makes sense whenever gg is smooth and compactly supported with dg(x)dx=0\int_{\mathbb{R}^{d}}g(x)\,dx=0. and that hh^{\infty} agrees in law, modulo additive constant, with the log-correlated Gaussian field from [lgf-survey, fgf-survey]. Furthermore, hhh^{\infty}-h has a modification which is a continuous function, viewed modulo additive constant. This was discussed in [lgf-survey, Section 4.1.1] and explained in detail in the two-dimensional case in [afs-metric-ball, Appendix B] (the same proof works for any dimension). Due to the continuity of hhh^{\infty}-h, one can easily deduce from Theorem 1.2 that a natural approximation scheme for the exponential metric associated with hh^{\infty} is also tight.

1.2 Related models

Since the construction of the LQG metric in [dddf-lfpp, gm-uniqueness], there have been several additional works which prove tightness and/or uniqueness for various random fractal metrics. Examples include the supercritical LQG metric [dg-supercritical-lfpp, dg-uniqueness] (as mentioned in Remark 1.3), the conformal loop ensemble chemical distance [miller-cle-metric], and the limit of long-range percolation on d\mathbb{Z}^{d} [baumler-long-range-perc, dfh-long-range-perc]. We also mention the directed landscape, a random directed metric on 2\mathbb{R}^{2} related to the KPZ universality class [dov-dl].

An important feature of LQG is its relation with two-dimensional Liouville conformal field theory (LCFT) rigorously constructed in [dkrv-lqg-sphere] and follow-up works. The framework of LCFT can produce exact solvability results for the area and length measures associated with LQG surfaces when the underlying field is well chosen. Recently, two-dimensional LCFT has been extended to even dimensions d4d\geq 4 in the papers [cercle-higher-dimension, dhks-even-dim]. Both of these works construct a log-correlated Gaussian field on a dd-dimensional manifold whose law is re-weighted according to the so-called Liouville action. In other words, these works carry out analogs of [dkrv-lqg-sphere] on certain dd-manifolds. It should be possible to use the results of the present paper to associate a random metric with the fields considered in [cercle-higher-dimension, dhks-even-dim], at least as a subsequential limit. As in the two-dimensional case, the exponent QQ of Proposition 1.1 should correspond to the background charge in [cercle-higher-dimension] (which is also called QQ).

In two dimensions, LQG is conjectured to describe the scaling limit of random planar maps. In particular, the LQG metric is believed to describe the scaling limit of the random planar maps equipped with their graph distance in, e.g., the Gromov-Hausdorff sense. This convergence has been rigorously established for uniform random planar maps toward LQG with γ=8/3\gamma=\sqrt{8/3} (ξ=1/6\xi=1/\sqrt{6}), but is open for other values of γ\gamma. More precisely, it was shown in [legall-uniqueness, miermont-brownian-map] that uniform random planar maps converge to a random metric space called the Brownian map, and in [lqg-tbm1, lqg-tbm2] that the Brownian map is equivalent to 8/3\sqrt{8/3}-LQG, as a metric space. See also [hs-cardy-embedding] for a stronger topology of convergence and Section 2.4 of [ddg-metric-survey] for further discussions.

It would be extremely interesting to find a natural discrete random geometry in dimension d3d\geq 3 whose scaling limit is described by one of the exponential random metrics considered in this paper (or some minor variant thereof).

In analogy with the case of uniform triangulations in two-dimensions (which converge to 8/3\sqrt{8/3}-LQG), a natural discrete model to consider is uniform triangulations of the dd-dimensional sphere, with nn\in\mathbbm{N} total dd-simplices. Such triangulations appear to be very difficult to analyze. For example, it is a well-known open problem to determine whether the number of triangulations of the three-sphere with nn total tetrahedra grows exponentially or superexponentially [dj-3-manifolds, gromov-spaces]. Moreover, simulations suggest that uniform triangulations of the three-sphere may not have interesting scaling limits when viewed as metric spaces, see, e.g., [bk-3d-simplicial, av-3d-simplicial, abkv-3d-vacuum, ckr-3d-entropy, hty-3d-phases, hin-3d-simulation]. We refer to the introductions of [dj-3-manifolds, bz-locally-constructible, budd-lionni-3-spheres] and the references therein for further discussion.

On the other hand, there are natural restricted classes of triangulations of dd-spheres which appear to be more tractable, and whose cardinality can be shown to grow exponentially in nn. Examples include locally constructible, constructible, shellable, and vertex-decomposable triangulations [dj-3-manifolds, bz-locally-constructible]. One could ask whether a uniform sample from any of these restricted classes converges in the Gromov-Hausdorff sense to the exponential metric associated with a log-correlated Gaussian field (or a field which locally looks like a log-correlated Gaussian field).

Another interesting class of dd-dimensional triangulations which one could consider are those which can be represented as the tangency graph of a sphere packing in d\mathbbm{R}^{d} (or the dd-sphere). Unlike for d=2d=2, there is not a simple criterion for when a graph can be represented by a sphere packing in d\mathbbm{R}^{d} for d3d\geq 3. In fact, for several values of dd, it is known that the problem of determining whether a given graph admits a sphere packing representation in d\mathbbm{R}^{d} is NP hard [hk-disks-and-balls]. In dimension two, circle packings and their links to random conformal geometry are fairly well-understood (see, e.g., the survey [nachmias-circle-packing]). In higher dimensions the theory is much less well-developed and likely to be much more difficult. But, a few results can be found, e.g., in [cr-rigidity, bc-sphere-packing, lee-conformal-growth]. One could look for a natural model of random sphere packings in d\mathbbm{R}^{d} whose scaling limit is described by the exponential of a log-correlated Gaussian field.

In another direction, connections between γ\gamma-LQG and random planar maps for general γ(0,2)\gamma\in(0,2) have been obtained using the framework of mating-of-trees theory [wedges], see the survey [ghs-mating-survey]. The recent work [BC23] presents an analog of mating-of-trees constructions in three dimensions. In a similar vein, the paper [budd-lionni-3-spheres] introduces a model of random triangulations of the three-sphere, decorated by a pair of trees, which is combinatorially tractable and has interesting geometric features. It is natural to wonder if either of these models are related to the exponential metrics for log-correlated fields in dimension three.

Recently, an analog of the Brownian map in dimension d3d\geq 3 was proposed in [ml-iterated-folding]. It is also natural to wonder whether this random metric space has any relation to the exponential metrics of log-correlated Gaussian fields, analogous to the aforementioned relationship between the Brownian map and 8/3\sqrt{8/3}-LQG.

1.3 Outline

Here, we outline the proof strategy of Theorem 1.2 and describe the content of each subsequent section.

1.3.1 Comparison to the two-dimensional case

First, let us highlight the main differences between the method in this paper and the methods used in the earlier works [ding-dunlap-lqg-fpp, ding-dunlap-lgd, df-lqg-metric, dddf-lfpp, dg-supercritical-lfpp] to establish the tightness of approximations of exponential metrics for log-correlated fields in two-dimensions. All the results in two dimensions rely crucially on RSW estimates, which give up-to-constants comparisons between quantiles of DnD_{n}-crossing lengths of rectangles in the “easy direction” and the “hard direction”, see e.g. Section 3 of [dddf-lfpp]. The arguments to prove these estimates are based on either approximate conformal invariance or on forcing paths to cross each other, neither of which works in higher dimensions. For this reason, we will use a fundamentally different approach to prove tightness which bypasses any direct proof of RSW estimates.

The first difference in our approach as compared to the two-dimensional case is that we initially use the median of the point-to-point distance, namely λn\lambda_{n} from (1.7), as the normalizing constant. In contrast, previous works use the median of the left-right crossing distance within a box as their normalizing constant (although these two medians are eventually proved to be equivalent up to a constant). The point-to-point distance is typically larger than the left-right crossing distance, which makes it easier to upper-bound other types of distances in terms of λn\lambda_{n}. We choose to work with the internal point-to-point distance inside a box to ensure that we have long-range independence, which allows us to employ percolation arguments. To use the percolation argument, we will actually work with the qq-quantile of Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)) for qq close to one, but not depending on nn, in most parts of the proof.

The second difference arises from the lack of concentration results for dimension 3\geq 3. In previous works, the authors have derived upper tail estimates for the left-to-right crossing distance, see Sections 4 and 5 of [dddf-lfpp] and Section 3 of [dg-supercritical-lfpp] for these types of results. These estimates are based on the RSW estimates, percolation arguments, and the Efron-Stein inequality. Despite the fact that the RSW argument is not applicable in our case, we can still hope to use the percolation argument to achieve an upper tail estimate in our setting. Simplistically, if we can compare λn\lambda_{n} and λnk\lambda_{n-k} (which should differ by at most a constant if kk is fixed), then we can divide a box into 2k2^{k} pieces. By using the scaling property of hnh_{n} (Lemma 2.2), the definition (1.7) of λnk\lambda_{n-k} (actually we will use a large quantile instead of the median) and percolation arguments, we can deduce an upper tail estimate for the DnD_{n}-distance across a hypercubic shell666A hypercubic shell is the domain between two concentric boxes, which is the dd-dimensional analog of a square annulus., e.g. Dn(B1(0),B2(0))D_{n}(\partial B_{1}(0),\partial B_{2}(0)), in terms of λn\lambda_{n}. It turns out that a specific comparison bound between λn\lambda_{n} and λnk\lambda_{n-k} for all integers 1kn1\leq k\leq n, as detailed in Proposition LABEL:prop:compare, is sufficient to achieve an upper tail estimate for the diameter of a box. This in turn ensures the tightness of the metric. Deriving this comparison is the most technical part of this paper and is detailed in Section LABEL:sec:compare. We will actually derive a comparison between the metrics DnD_{n} and DnkD_{n-k} in that section, which may also be of independent interest.

The third difference also arises from the lack of concentration results for dimension 3\geq 3 and our choice of the normalizing constant. In previous works, a lower tail estimate for the left-right crossing distance follows from the RSW argument, see Section 4 of [dddf-lfpp], and this implies that each subsequential limit is a metric. Here we will use a different approach to demonstrate this. Note that, a prior, the point-to-point distance can be much larger than the left-to-right crossing distance of a box. Our strategy begins by showing that the distance across a hypercubic shell is positive with non-zero probability. Combining with a zero-one law argument (Lemma LABEL:lem:zero-one), we can increase this probability to one, thereby establishing that each subsequential limit is a metric. This will be detailed in Subsection LABEL:subsec:non-degenerate.

1.3.2 Detailed outline

Next, we describe our strategy in more detail and outline the content of each section. More comprehensive overviews can be found at the beginning of each respective section and subsection.

In Section 2, we provide preliminaries and fix some notation. Let W(dx,dt)W(dx,dt) be the space-time white noise. Throughout this paper, we will work with the approximation of the log-correlated Gaussian field

hm,n(x)=d2n2mK(xyt)td+12dydth_{m,n}(x)=\int_{\mathbb{R}^{d}}\int_{2^{-n}}^{2^{-m}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}(\frac{x-y}{t})t^{-\frac{d+1}{2}}\,dy\,dt

for integers n>m0n>m\geq 0 (note that h0,n=hnh_{0,n}=h_{n}, as defined in (1.4)). Some basic properties and estimates of hm,nh_{m,n} are provided in Subsection 2.2. In Subsection 2.3, we define Dm,nD_{m,n} as the exponential metrics associated with hm,nh_{m,n} and establish some basic properties of these metrics, including a Gaussian concentration bound (Lemma 2.7). Subsection 2.4 collects basic arguments about percolation with finite range of dependence, which will play a crucial role in Sections 3 and LABEL:sec:bound-distance.

In Section 3, we will prove Proposition 3.1, which establishes the existence of an exponent QQ satisfying (1.8). This follows from the subadditivity inequality: λneCn2/3λmλnm\lambda_{n}\leq e^{Cn^{2/3}}\lambda_{m}\lambda_{n-m} for integers n>m1n>m\geq 1, and the proof is similar to Proposition 2.5 of [dg-supercritical-lfpp]. This inequality requires constructing a path connecting 0 and e1e_{1} of typical DnD_{n}-length. The construction essentially follows two steps. The first step is to construct a path on 2md2^{-m}\mathbb{Z}^{d} such that its DmD_{m}-length can be upper-bounded. The second step is to modify the path locally so that its DnD_{n}-length can be controlled using a percolation argument on a refined lattice. To use the percolation argument, we actually consider a large quantile of Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)). In Lemma LABEL:lem:Q-lower, we will establish basic properties of Q(ξ)Q(\xi), based on estimates for hnh_{n}. This will conclude the proof of Proposition 1.1.

In Section LABEL:sec:bound-distance, we will establish a chaining argument similar to the ones in Section 6.3 of [ding-dunlap-lqg-fpp] and Section 6.1 of [df-lqg-metric], and derive bounds for different types of distances. First, we present the chaining argument in Subsection LABEL:subsec:chaining. We use paths of typical DnD_{n}-length at different scales to connect any two points in a box, and thus establish an upper-bound for the DnD_{n}-diameter of a box in terms of the large quantiles of Dnm(0,e1;B2(0))D_{n-m}(0,e_{1};B_{2}(0)) for 0mn0\leq m\leq n. This result will be used subsequently in two places. First, we will use it in Subsections LABEL:subsec:bound-diameter and LABEL:subsec:cross to show that the medians of the diameter of a box or distance across a hypercubic shell all satisfy the relation in (1.8). Secondly, this result will be used to prove the tightness of the metric in Subsection LABEL:subsec:tightness after a comparison between quantiles of Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)) for different values of nn is achieved. In Subsection LABEL:subsec:super-exponential, we will establish super-exponential concentration bounds for distances across and around hypercubic shells, which will be used in Section LABEL:sec:compare.

In Section LABEL:sec:compare, we will compare the metrics DnD_{n} and Dn+kD_{n+k} for integers nk1n\geq k\geq 1. We will briefly describe the strategy here, and refer to Subsection LABEL:subsec:sec5-strategy for a more detailed outline of this section. The comparison is achieved by controlling the behavior of the field hn,n+kh_{n,n+k} (note that Dn+kD_{n+k} is obtained from DnD_{n} by adding hn,n+kh_{n,n+k} to the field). In most parts of the space, hn,n+kh_{n,n+k} behaves well, and Dn+kD_{n+k} and DnD_{n} satisfy the desired bound given in Proposition LABEL:prop:compare. However, there are places where hn,n+kh_{n,n+k} does not behave well, and a priori, it is possible that a DnD_{n}- or Dn+kD_{n+k}-geodesic spends most of its time in these problematic regions. Our main effort is to address these regions. The proof essentially involves two steps. In the first step, we use a coarse-graining argument to show that, with high probability, we can find boxes at different scales to cover the problematic regions. Importantly, all these boxes satisfy the condition that the DnD_{n}-distance around the hypercubic shell enclosing the box can be upper-bounded by the DnD_{n}-distance across a larger hypercubic shell. In the second step, we use this condition to show that the ill-behaved field within these boxes has a minor impact on the metric DnD_{n}. Specifically, paths can be modified to avoid these boxes, and their DnD_{n}-length will increase by no more than a constant factor. Moreover, for paths entirely contained within the domain where hn,n+kh_{n,n+k} behaves well, by adjusting the paths, their DnD_{n}-length and Dn+kD_{n+k}-length satisfy the desired bound. This leads to a comparison between DnD_{n} and Dn+kD_{n+k}.

In Section LABEL:sec:final-proof, we prove Theorem 1.2. The proof consists of two parts. In Subsection LABEL:subsec:tightness, we combine results from the chaining argument in Subsection LABEL:subsec:chaining and the comparison of quantiles from Section LABEL:sec:compare to demonstrate the tightness of DnD_{n} when normalized by the qq-quantile of Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)), where qq is close to one but independent from nn. In Subsection LABEL:subsec:non-degenerate, we will establish that each possible subsequential limit is a metric. From the definition of quantiles and the positive association (FKG) property for positively correlated Gaussian processes, we first show that the distance across a hypercubic shell is bounded away from zero with positive probability. Another crucial input is a zero-one law (Lemma LABEL:lem:zero-one), which is derived from the locality property of the metric. By using this argument, we can increase the probability to one. Applying this to countably many hypercubic shells shows that the subsequential limit is a metric, which in turn implies an up-to-constants comparison between the median λn\lambda_{n} and the qq-quantile of Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)). This gives tightness when we normalize by λn\lambda_{n} instead of by the qq-quantile.

In Section LABEL:sec:open-problem, we list some open problems related to the metric we constructed. Appendix LABEL:appendix:index includes a list of notation that we will use in this paper.

2 Preliminaries

2.1 Basic notation

Numbers

We write ={1,2,}\mathbb{N}=\{1,2,\ldots\}. Without specific mention, the logarithm in this paper will be taken with respect to the base ee. For aa\in\mathbb{R}, we use a\lfloor a\rfloor to represent the largest integer not greater than aa. For a random variable XX, we will use Med(X){\rm Med}(X) to represent its median.

Metrics

Let (X,D)(X,D) be a metric space. For a curve P:[a,b]XP:[a,b]\rightarrow X, the DD-length of PP is defined as

len(P;D):=supTi=1nD(P(ti),P(ti1)){\rm len}(P;D):=\sup_{T}\sum_{i=1}^{n}D(P(t_{i}),P(t_{i-1}))

where the supremum is taken over all partitions T:a=t0<t1<<tn=bT:a=t_{0}<t_{1}<\ldots<t_{n}=b of [a,b][a,b]. The DD-length of a curve may be infinite.

For a curve P:[a,b]XP:[a,b]\rightarrow X and a set YXY\subset X, consider the pre-image P1(Y)[a,b]P^{-1}(Y)\subset[a,b]. Write the interior of P1(Y)P^{-1}(Y) as the disjoint union of countably many open intervals {(ai,bi)}i1\{(a_{i},b_{i})\}_{i\geq 1}. We define the restriction of PP to YY as P|Y:=i1P[ai,bi]P|_{Y}:=\cup_{i\geq 1}P[a_{i},b_{i}], which is the union of a family of curves, and its length is defined as

len(P|Y;D):=i1len(P[ai,bi];D).{\rm len}(P|_{Y};D):=\sum_{i\geq 1}{\rm len}(P[a_{i},b_{i}];D)\,. (2.1)

Note that P|YP|_{Y} and P(P1(Y))P(P^{-1}(Y)) are the same up to the inclusion of end points (of intervals in P1(Y)P^{-1}(Y)) or single points (i.e., each interval containing them is not a subset of P1(Y)P^{-1}(Y)). For the sets YY that we will consider in this paper, the lengths of P|YP|_{Y} and P(P1(Y))P(P^{-1}(Y)) will be the same.

For YXY\subset X, the internal metric of DD on YY is defined as

D(x,y;Y):=infPYlen(P;D),x,yYD(x,y;Y):=\inf_{P\subset Y}{\rm len}(P;D),\quad\forall x,y\in Y (2.2)

where the infimum is taken over all paths PP in YY from xx to yy. Then D(,;Y)D(\cdot,\cdot;Y) is a metric on YY, allowing the distance between two points to be infinite.

We say DD is a length metric if for all x,yXx,y\in X and δ>0\delta>0, there exists a curve with DD-length at most D(x,y)+δD(x,y)+\delta connecting xx and yy. We say DD is a geodesic metric if for each x,yXx,y\in X, there exists a curve with DD-length precisely D(x,y)D(x,y) connecting xx and yy.

Subsets of Euclidean space

In this paper, we consider the space d\mathbb{R}^{d} where d2d\geq 2 is a fixed dimension. For zdz\in\mathbb{R}^{d}, we write z=(z1,,zd)z=(z_{1},\ldots,z_{d}) for its coordinates. We use the notation ||1|\cdot|_{1}, ||2|\cdot|_{2}, and |||\cdot|_{\infty} to represent the l1l^{1}-, l2l^{2}-, and ll^{\infty}-norms, respectively. We use 𝔡1\mathfrak{d}_{1}, 𝔡2\mathfrak{d}_{2}, and 𝔡\mathfrak{d}_{\infty} to denote the distances associated with these norms. Without specific mention, the distance that we use is the ll^{\infty}-distance. For a set AdA\subset\mathbb{R}^{d} and r>0r>0, we define the ll^{\infty}-neighborhood

Br(A):={zd:𝔡(z,A)<r}.B_{r}(A):=\{z\in\mathbb{R}^{d}:\mathfrak{d}_{\infty}(z,A)<r\}\,.

As in (1.3), for zdz\in\mathbb{R}^{d}, we write Br(z)=Br({z})=z+(r,r)dB_{r}(z)=B_{r}(\{z\})=z+(-r,r)^{d} for the open box centered at zz with side-length 2r2r. We call a domain AdA\subset\mathbb{R}^{d} a hypercubic shell if A=Br1(x)\Br2(x)A=B_{r_{1}}(x)\backslash B_{r_{2}}(x) for some xdx\in\mathbb{R}^{d} and r1>r2>0r_{1}>r_{2}>0.

We extend the notation of |||\cdot|_{\infty}, ||1|\cdot|_{1}, 𝔡\mathfrak{d}_{\infty}, and 𝔡1\mathfrak{d}_{1} to the integer lattice d\mathbb{Z}^{d}. For xdx\in\mathbb{Z}^{d} and an integer n0n\geq 0, we define Bn(x)B_{n}(x) as the box centered at xx with side-length 2n2n. Namely,

Bn(x):={zd:|zx|n}.B_{n}(x):=\{z\in\mathbb{Z}^{d}:|z-x|_{\infty}\leq n\}\,.

We will clarify in the context whether we are considering xx as a point in d\mathbb{R}^{d} or as a vertex in d\mathbb{Z}^{d}. For an integer n0n\geq 0, define the set

Ln:=2ndB2(0).{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{n}:=2^{-n}\mathbb{Z}^{d}\cap B_{2}(0)\,. (2.3)

Typically, we consider Ln{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{n} as a subset of d\mathbb{R}^{d}. However, when analyzing (*-)paths or (*-)clusters on the rescaled lattice 2nd2^{-n}\mathbb{Z}^{d}, as defined in Subsection 2.4, we view Ln{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{n} as a subset of 2nd2^{-n}\mathbb{Z}^{d}. In this paper, we will also consider the graph distance on the rescaled lattice εd\varepsilon\mathbb{Z}^{d} which is defined as 1/ε1/\varepsilon times the ll^{\infty}-distance when considering εd\varepsilon\mathbb{Z}^{d} as a subset of d\mathbb{R}^{d}.

Convention about constants

Constants like c,c,C,Cc,c^{\prime},C,C^{\prime} may change from place to place, while constants with subscripts like c1,C1c_{1},C_{1} remain fixed throughout the article. All constants may implicitly rely on the dimension dd, the kernel K, r0, and ξ\xi. The dependence on additional variables will be indicated at the first occurrence of each constant.

2.2 Approximation of a log-correlated Gaussian field

In this subsection, we establish some basic properties of the Gaussian random functions hnh_{n} introduced in Subsection 1.1. Let us fix a convolution kernel K:d[0,){\hyperref@@ii[K-condition1]{\mathfrak{K}}}:\mathbb{R}^{d}\rightarrow[0,\infty) and a constant r0>0{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}>0 which satisfy the conditions 1, 2, 3 in Subsection 1.1. Let WW be a white noise on d×(0,)\mathbb{R}^{d}\times(0,\infty) and we define hnh_{n} and hh as in (1.4). We will also have occasion to consider the following additional functions.

Definition 2.1.

For integers nm0n\geq m\geq 0 and xdx\in\mathbb{R}^{d}, we define

hm,n(x):=hn(x)hm(x)=d2n2mK(yxt)td+12W(dy,dt).h_{m,n}(x):=h_{n}(x)-h_{m}(x)=\int_{\mathbb{R}^{d}}\int_{2^{-n}}^{2^{-m}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}\big{(}\frac{y-x}{t}\big{)}t^{-\frac{d+1}{2}}W(dy,dt)\,.

Note that h0,n(x)=hn(x)h_{0,n}(x)=h_{n}(x).

The following properties of hm,nh_{m,n} follow directly from its definition and the conditions on K. We omit the proof here.

Lemma 2.2.

For integers n>m0n>m\geq 0, we have

  1. 1.

    hm,nh_{m,n} is smooth.

  2. 2.

    The law of hm,nh_{m,n} is invariant under translation and rotation of d\mathbb{R}^{d}.

  3. 3.

    For any U,VdU,V\subset\mathbb{R}^{d} with 𝔡(U,V)2r02m\mathfrak{d}_{\infty}(U,V)\geq 2{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}\cdot 2^{-m}, the fields hm,n|Uh_{m,n}|_{U} and hm,n|Vh_{m,n}|_{V}, which are obtained by restricting hm,nh_{m,n} to the domains UU and VV, are independent.

  4. 4.

    The fields satisfy the scaling property: (hm,n(x))xd=d(h0,nm(x2m))xd(h_{m,n}(x))_{x\in\mathbb{R}^{d}}\overset{d}{=}(h_{0,n-m}(x2^{m}))_{x\in\mathbb{R}^{d}}.

We collect some basic estimates about the field in the following lemma:

Lemma 2.3.
  1. 1.

    For any integers n>m0n>m\geq 0 and xdx\in\mathbb{R}^{d}, Var(hm,n(x))=(nm)log2{\rm Var}(h_{m,n}(x))=(n-m)\log 2.

  2. 2.

    There exists some constant C>0C>0 such that for all n1n\geq 1 and t>0t>0:

    [supx2nB1(0)|hn(x)|2nt]Cet2C.\mathbb{P}\Big{[}\sup_{x\in 2^{-n}B_{1}(0)}|\nabla h_{n}(x)|_{\infty}\geq 2^{n}t\Big{]}\leq Ce^{-\frac{t^{2}}{C}}. (2.4)
  3. 3.

    There exists some constant C>0C>0 such that for all n1n\geq 1, we have

    𝔼[supxB1(0)hn(x)]n2dlog2+Cn.\mathbb{E}\Big{[}\sup_{x\in B_{1}(0)}h_{n}(x)\Big{]}\leq n\sqrt{2d}\log 2+C\sqrt{n}\,.
  4. 4.

    (Borell-TIS inequality) For all u>0u>0 and integer n1n\geq 1, we have

    [supxB1(0)hn(x)𝔼supxB1(0)hn(x)+u]exp(u22log2n).\mathbb{P}\Big{[}\sup_{x\in B_{1}(0)}h_{n}(x)\geq\mathbb{E}\sup_{x\in B_{1}(0)}h_{n}(x)+u\Big{]}\leq\exp\big{(}-\frac{u^{2}}{2\log 2\cdot n}\big{)}\,. (2.5)
  5. 5.

    There exists some constant C>0C>0 such that for all u>0u>0 and integer n1n\geq 1, we have

    [supx2nB1(0)hn(x)>u]Cexp((uu2/3)22log2n)+Cexp(u4/3C).\mathbb{P}\Big{[}\sup_{x\in 2^{-n}B_{1}(0)}h_{n}(x)>u\Big{]}\leq C\exp\big{(}-\frac{(u-u^{2/3})^{2}}{2\log 2\cdot n}\big{)}+C\exp\big{(}-\frac{u^{4/3}}{C}\big{)}\,.
Proof.

We first prove Claim (1). By using the property of white noise and the identity dK(x)2dx=1\int_{\mathbb{R}^{d}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}(x)^{2}dx=1 from condition 3, we obtain:

Var(hm,n(x))\displaystyle\quad{\rm Var}(h_{m,n}(x))
=𝔼[d2n2mK(yxt)td+12W(dy,dt)\displaystyle=\mathbb{E}\bigg{[}\int_{\mathbb{R}^{d}}\int_{2^{-n}}^{2^{-m}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}\big{(}\frac{y-x}{t}\big{)}t^{-\frac{d+1}{2}}W(dy,dt)
×d2n2mK(yxt)td+12W(dy,dt)]\displaystyle\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times\int_{\mathbb{R}^{d}}\int_{2^{-n}}^{2^{-m}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}\big{(}\frac{y^{\prime}-x}{t^{\prime}}\big{)}t^{\prime-\frac{d+1}{2}}W(dy^{\prime},dt^{\prime})\bigg{]}
=d2n2mK(yxt)2td1dydt=2n2mt1dt=(nm)log2.\displaystyle=\int_{\mathbb{R}^{d}}\int_{2^{-n}}^{2^{-m}}{\hyperref@@ii[K-condition1]{\mathfrak{K}}}\big{(}\frac{y-x}{t}\big{)}^{2}t^{-d-1}dydt=\int_{2^{-n}}^{2^{-m}}t^{-1}dt=(n-m)\log 2\,.

We now prove Claim (2). Using the smoothness of K and Fernique’s theorem (see e.g. [fernique-criterion]), we have a tail estimate for h1h_{1}. That is, there exists a constant A>0A>0 such that for all t>0t>0:

[supxB1(0)|h1(x)|t]Aet2A.\mathbb{P}\Big{[}\sup_{x\in B_{1}(0)}|\nabla h_{1}(x)|_{\infty}\geq t\Big{]}\leq Ae^{-\frac{t^{2}}{A}}. (2.6)

By Claim (4) in Lemma 2.2, for any integer kk, we have (hk,k+1(x))x2kB1(0)=d(h0,1(x2k))x2kB1(0)(h_{k,k+1}(x))_{x\in 2^{-k}B_{1}(0)}\overset{d}{=}(h_{0,1}(x2^{k}))_{x\in 2^{-k}B_{1}(0)}. Therefore,

supx2kB1(0)|hk,k+1(x)|=d2ksupyB1(0)|h0,1(y)|.\sup_{x\in 2^{-k}B_{1}(0)}|\nabla h_{k,k+1}(x)|_{\infty}\overset{d}{=}2^{k}\sup_{y\in B_{1}(0)}|\nabla h_{0,1}(y)|_{\infty}\,.

Combining this with (2.6), we obtain that for all integer k0k\geq 0 and t>0t>0:

[supx2kB1(0)|hk,k+1(x)|2kt]Aet2A.\mathbb{P}\Big{[}\sup_{x\in 2^{-k}B_{1}(0)}|\nabla h_{k,k+1}(x)|_{\infty}\geq 2^{k}t\Big{]}\leq Ae^{-\frac{t^{2}}{A}}. (2.7)

Let us first prove (2.4) in the case where t1t\geq 1. Using the facts

supx2nB1(0)|hn(x)|k=0n1supx2kB1(0)|hk,k+1(x)|\sup_{x\in 2^{-n}B_{1}(0)}|\nabla h_{n}(x)|_{\infty}\leq\sum_{k=0}^{n-1}\sup_{x\in 2^{-k}B_{1}(0)}|\nabla h_{k,k+1}(x)|_{\infty}

and k=0n12k2nk242n\sum_{k=0}^{n-1}2^{k}\cdot 2^{\frac{n-k}{2}}\leq 4\cdot 2^{n}, we obtain that for all integer n1n\geq 1:

[supx2nB1(0)|hn(x)|2nt]\displaystyle\mathbb{P}\Big{[}\sup_{x\in 2^{-n}B_{1}(0)}|\nabla h_{n}(x)|_{\infty}\geq 2^{n}t\Big{]}
[k=0n1supx2kB1(0)|hk,k+1(x)|k=0n12k2nk2(t/4)]\displaystyle\qquad\leq\mathbb{P}\Big{[}\sum_{k=0}^{n-1}\sup_{x\in 2^{-k}B_{1}(0)}|\nabla h_{k,k+1}(x)|_{\infty}\geq\sum_{k=0}^{n-1}2^{k}\cdot 2^{\frac{n-k}{2}}(t/4)\Big{]}
k=0n1[supx2kB1(0)|hk,k+1(x)|2k2nk2(t/4)].\displaystyle\qquad\leq\sum_{k=0}^{n-1}\mathbb{P}\Big{[}\sup_{x\in 2^{-k}B_{1}(0)}|\nabla h_{k,k+1}(x)|_{\infty}\geq 2^{k}\cdot 2^{\frac{n-k}{2}}(t/4)\Big{]}.

Using (2.7) and the fact that t1t\geq 1, we can choose a constant C>0C>0 depending only on AA such that:

[supx2nB1(0)|hn(x)|2nt]k=0n1Ae2nk(t/4)2ACet2C.\displaystyle\mathbb{P}\Big{[}\sup_{x\in 2^{-n}B_{1}(0)}|\nabla h_{n}(x)|_{\infty}\geq 2^{n}t\Big{]}\leq\sum_{k=0}^{n-1}Ae^{-\frac{2^{n-k}(t/4)^{2}}{A}}\leq Ce^{-\frac{t^{2}}{C}}.

This result can be extended to all t>0t>0 by enlarging the value of CC, thereby proving Claim (2).

Next, we prove Claim (3). Using the fact

supxB1(0)hn(x)supxB1(0)2ndhn(x)+2ndsupyB2n(x)xB1(0)2nd|hn(y)|,\sup_{x\in B_{1}(0)}h_{n}(x)\leq\sup_{x\in B_{1}(0)\cap 2^{-n}\mathbb{Z}^{d}}h_{n}(x)+2^{-n}d\sup_{\begin{subarray}{c}y\in B_{2^{-n}}(x)\\ x\in B_{1}(0)\cap 2^{-n}\mathbb{Z}^{d}\end{subarray}}|\nabla h_{n}(y)|_{\infty}\,,

we obtain that for all integer n1n\geq 1 and s>0s>0

[supxB1(0)hn(x)n2dlog2+sn]\displaystyle\mathbb{P}\Big{[}\sup_{x\in B_{1}(0)}h_{n}(x)\geq n\sqrt{2d}\log 2+s\sqrt{n}\Big{]}
[supxB1(0)2ndhn(x)n2dlog2+sn/2]\displaystyle\qquad\leq\mathbb{P}\Big{[}\sup_{x\in B_{1}(0)\cap 2^{-n}\mathbb{Z}^{d}}h_{n}(x)\geq n\sqrt{2d}\log 2+s\sqrt{n}/2\Big{]}
+[2ndsupyB2n(x)xB1(0)2nd|hn(y)|sn/2].\displaystyle\qquad\qquad+\mathbb{P}\Big{[}2^{-n}d\sup_{\begin{subarray}{c}y\in B_{2^{-n}}(x)\\ x\in B_{1}(0)\cap 2^{-n}\mathbb{Z}^{d}\end{subarray}}|\nabla h_{n}(y)|_{\infty}\geq s\sqrt{n}/2\Big{]}.

Using Claims (1) and (2), translation invariance of hnh_{n}, and the fact that |B1(0)2nd|C2dn|B_{1}(0)\cap 2^{-n}\mathbb{Z}^{d}|\leq C2^{dn}, we have

[supxB1(0)hn(x)n2dlog2+sn]\displaystyle\mathbb{P}\Big{[}\sup_{x\in B_{1}(0)}h_{n}(x)\geq n\sqrt{2d}\log 2+s\sqrt{n}\Big{]} (2.8)
C2dnexp((n2dlog2+sn/2)22log2n)+C2dnexp(s2nC),\displaystyle\qquad\leq C2^{dn}\exp\Big{(}-\frac{(n\sqrt{2d}\log 2+s\sqrt{n}/2)^{2}}{2\log 2\cdot n}\Big{)}+C2^{dn}\exp\Big{(}-\frac{s^{2}n}{C}\Big{)}\,,

where we enlarged the value of CC. When ss is large enough (independent of nn), the right-hand side is smaller than Ces2/CCe^{-s^{2}/C}. By integrating (2.8) with respect to ss, we obtain Claim (3).

Claim (4) follows from the Borell-TIS inequality (see [borell-tis1, TIS76], and also [adler-taylor-fields, Theorem 2.1.1]) and the fact that Var(hn(x))=nlog2{\rm Var}(h_{n}(x))=n\log 2 as stated in Claim (1).

Finally, we prove Claim (5). Using the fact

supx2nB1(0)hn(x)hn(0)+d2nsupx2nB1(0)|hn(x)|,\sup_{x\in 2^{-n}B_{1}(0)}h_{n}(x)\leq h_{n}(0)+d2^{-n}\sup_{x\in 2^{-n}B_{1}(0)}|\nabla h_{n}(x)|_{\infty}\,,

we obtain

[supx2nB1(0)hn(x)>u]\displaystyle\mathbb{P}\Big{[}\sup_{x\in 2^{-n}B_{1}(0)}h_{n}(x)>u\Big{]}
[hn(0)>uu2/3]+[supx2nB1(0)2n|hn(x)|>u2/3/d].\displaystyle\qquad\leq\mathbb{P}\big{[}h_{n}(0)>u-u^{2/3}\big{]}+\mathbb{P}\Big{[}\sup_{x\in 2^{-n}B_{1}(0)}2^{-n}|\nabla h_{n}(x)|_{\infty}>u^{2/3}/d\Big{]}\,.

Applying Claims (1) and (2) gives the desired result. ∎

2.3 Definition of the exponential metric

In this subsection, we introduce the exponential metric associated with hm,nh_{m,n}, which is the main focus of this paper. We also establish some of its basic properties.

Definition 2.4.

Fix ξ>0\xi>0. For integers nm0n\geq m\geq 0, we define the exponential metric associated with the field hm,nh_{m,n} from Definition 2.1 as follows:

Dm,n(z,w):=infP:zw01eξhm,n(P(t))|P(t)|dt,D_{m,n}(z,w):=\inf_{P:z\to w}\int_{0}^{1}e^{\xi h_{m,n}(P(t))}|P^{\prime}(t)|dt\,,

where the infimum is taken over all piecewise continuously differentiable paths P:[0,1]dP:[0,1]\rightarrow\mathbb{R}^{d} joining z,wz,w. For an open set UdU\subset\mathbb{R}^{d}, we define the internal metric Dm,n(,;U)D_{m,n}(\cdot,\cdot;U) as described in (2.2). When m=0m=0, the metric D0,nD_{0,n} is the same as the metric DnD_{n} introduced in (1.6). When m=nm=n, Dm,nD_{m,n} is equivalent to the Euclidean metric.

The following lemma is a direct consequence of Claims (2) and (4) in Lemma 2.2. We omit the proof here.

Lemma 2.5.

For integers nm0n\geq m\geq 0 and any open set UdU\subset\mathbb{R}^{d} (including U=dU=\mathbb{R}^{d}), we have

  1. 1.

    The law of Dm,n(,;U)D_{m,n}(\cdot,\cdot;U) is invariant under translation and rotation of d\mathbb{R}^{d}.

  2. 2.

    The law of Dm,n(,;U)D_{m,n}(\cdot,\cdot;U) satisfies the scaling property:

    (Dm,n(x,y;U))x,yU=d2m(Dnm(2mx,2my,2mU))x,yU.(D_{m,n}(x,y;U))_{x,y\in U}\overset{d}{=}2^{-m}(D_{n-m}(2^{m}x,2^{m}y,2^{m}U))_{x,y\in U}\,.

As a corollary of Claim (3) in Lemma 2.2, we have that the internal metrics of Dm,nD_{m,n} are independent within two domains located far from each other.

Lemma 2.6.

For integers n>m0n>m\geq 0 and any open sets U,VdU,V\subset\mathbb{R}^{d} with 𝔡(U,V)2r02m\mathfrak{d}_{\infty}(U,V)\geq 2{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}\cdot 2^{-m}, the internal metrics Dm,n(,;U)D_{m,n}(\cdot,\cdot;U) and Dm,n(,;V)D_{m,n}(\cdot,\cdot;V) are independent.

Proof.

The internal metric Dm,n(,;U)D_{m,n}(\cdot,\cdot;U) is determined by hm,n|Uh_{m,n}|_{U}, and the internal metric Dm,n(,;V)D_{m,n}(\cdot,\cdot;V) is determined by hm,n|Vh_{m,n}|_{V}. By Claim (3) in Lemma 2.2, we obtain the result. ∎

We prove a concentration bound for the exponential metric. The proof is similar to that of [dddf-lfpp, Lemma 23].

Lemma 2.7.

For all open subset UdU\subset\mathbb{R}^{d} or U=dU=\mathbb{R}^{d}, compact subsets K1,K2UK_{1},K_{2}\subset U that are connected in UU and K1K2=K_{1}\cap K_{2}=\emptyset, integers n>m0n>m\geq 0, and t>0t>0, the following concentration bound holds:

[|logDm,n(K1,K2;U)𝔼logDm,n(K1,K2;U)|>t]2et22ξ2log2(nm).\mathbb{P}\big{[}|\log D_{m,n}(K_{1},K_{2};U)-\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)|>t\big{]}\leq 2e^{-\frac{t^{2}}{2\xi^{2}\log 2\cdot(n-m)}}. (2.9)
Proof.

We first show that |𝔼logDm,n(K1,K2;U)|<|\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)|<\infty. Let us begin with the upper bound. By the assumption, there exists a large constant NN such that K1K_{1} and K2K_{2} are connected by a path of Euclidean length at most NN in UBN(0)U\cap B_{N}(0). Therefore,

𝔼logDm,n(K1,K2;U)𝔼log(NeξsupxUBN(0)hm,n(x))<.\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)\leq\mathbb{E}\log\big{(}Ne^{\xi\sup_{x\in U\cap B_{N}(0)}h_{m,n}(x)}\big{)}<\infty\,.

The last inequality follows from the Gaussian tail of supxUBN(0)hm,n(x)\sup_{x\in U\cap B_{N}(0)}h_{m,n}(x), as indicated by Claim (2) in Lemma 2.3. Furthermore, there exists a large constant MM such that any path connecting K1K_{1} and K2K_{2} must have a Euclidean length of at least 1M\frac{1}{M} within UBM(0)U\cap B_{M}(0). Therefore,

𝔼logDm,n(K1,K2;U)𝔼log(1MeξinfxUBM(0)hm,n(x))>.\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)\geq\mathbb{E}\log\big{(}\frac{1}{M}e^{\xi\inf_{x\in U\cap B_{M}(0)}h_{m,n}(x)}\big{)}>-\infty\,.

Combining the above two inequalities yields |𝔼logDm,n(K1,K2;U)|<|\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)|<\infty.

We now prove (2.9) first for a bounded open set UU. For integer k1k\geq 1, let Dm,n(k)D_{m,n}^{(k)} be the exponential metric associated with hm,n(k)h_{m,n}^{(k)}, where hm,n(k)h_{m,n}^{(k)} is piecewise constant and takes the value hm,n(x)h_{m,n}(x) on each dyadic box B2k(x)B_{2^{-k}}(x) for xd2kdx\in\mathbb{R}^{d}\cap 2^{-k}\mathbb{Z}^{d}. Then, supxU|hm,n(x)hm,n(k)(x)|d2ksupxB1(U)|hm,n(x)|\sup_{x\in U}|h_{m,n}(x)-h_{m,n}^{(k)}(x)|\leq d2^{-k}\sup_{x\in B_{1}(U)}|\nabla h_{m,n}(x)|_{\infty}. This, combined with Definition 2.4, implies that

eξd2ksupxB1(U)|hm,n(x)|Dm,n(K1,K2;U)Dm,n(k)(K1,K2;U)eξd2ksupxB1(U)|hm,n(x)|.e^{-\xi d2^{-k}\sup_{x\in B_{1}(U)}|\nabla h_{m,n}(x)|_{\infty}}\leq\frac{D_{m,n}(K_{1},K_{2};U)}{D_{m,n}^{(k)}(K_{1},K_{2};U)}\leq e^{\xi d2^{-k}\sup_{x\in B_{1}(U)}|\nabla h_{m,n}(x)|_{\infty}}.

Together with the fact that 𝔼[supxB1(U)|hm,n(x)|]<\mathbb{E}[\sup_{x\in B_{1}(U)}|\nabla h_{m,n}(x)|_{\infty}]<\infty (because it has a Gaussian tail, as indicated by Claim (2) in Lemma 2.3), we obtain:

limkDm,n(k)(K1,K2;U)=Dm,n(K1,K2;U)and\displaystyle\lim_{k\rightarrow\infty}D_{m,n}^{(k)}(K_{1},K_{2};U)=D_{m,n}(K_{1},K_{2};U)\quad\mbox{and} (2.10)
limk𝔼logDm,n(k)(K1,K2;U)=𝔼logDm,n(K1,K2;U).\displaystyle\lim_{k\rightarrow\infty}\mathbb{E}\log D_{m,n}^{(k)}(K_{1},K_{2};U)=\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)\,.

By definition, logDm,n(k)(K1,K2;U)\log D_{m,n}^{(k)}(K_{1},K_{2};U) is ξ\xi-Lipschitz as a function of

(Y1,,Yp):=(hm,n(k)(x))xB1(U)2kd(Y_{1},\ldots,Y_{p}):=(h_{m,n}^{(k)}(x))_{x\in B_{1}(U)\cap 2^{-k}\mathbb{Z}^{d}}

in terms of the ll^{\infty}-norm. In addition, there exists a p×pp\times p matrix AA such that (Y1,,Yp)=dA(X1,,Xp)(Y_{1},\ldots,Y_{p})^{\intercal}\overset{d}{=}A(X_{1},\ldots,X_{p})^{\intercal}, where X1,,XpX_{1},\ldots,X_{p} are i.i.d. standard Gaussian random variables. By Claim (1) in Lemma 2.3, the l2l^{2}-norm of each row of AA equals to Varhm,n(k)(x)=(nm)log2\sqrt{\operatorname{Var}h_{m,n}^{(k)}(x)}=\sqrt{(n-m)\log 2}. Therefore, logDm,n(k)(K1,K2;U)\log D_{m,n}^{(k)}(K_{1},K_{2};U), as a function of (X1,,Xp)(X_{1},\ldots,X_{p}), is ξ(nm)log2\xi\sqrt{(n-m)\log 2}-Lipschitz in terms of the l2l^{2}-norm. By the Gaussian concentration inequality (see [borell-tis1, TIS76], and also [adler-taylor-fields, Lemma 2.1.6]), we have for all t>0t>0:

[|logDm,n(k)(K1,K2;U)𝔼logDm,n(k)(K1,K2;U)|>t]2et22ξ2log2(nm).\mathbb{P}\big{[}|\log D_{m,n}^{(k)}(K_{1},K_{2};U)-\mathbb{E}\log D_{m,n}^{(k)}(K_{1},K_{2};U)|>t\big{]}\leq 2e^{-\frac{t^{2}}{2\xi^{2}\log 2\cdot(n-m)}}. (2.11)

By sending kk to infinity and combining with (2.10), we obtain the desired lemma in the case where UU is bounded.

We can extend the result to arbitrary UU by considering the truncation UBN(0)U\cap B_{N}(0) for integers N1N\geq 1. Since Dm,n(K1,K2;UBN(0))D_{m,n}(K_{1},K_{2};U\cap B_{N}(0)) decreases to Dm,n(K1,K2;U)D_{m,n}(K_{1},K_{2};U) as NN goes to infinity, we have:

limNDm,n(K1,K2;UBN(0))=Dm,n(K1,K2;U),and\lim_{N\rightarrow\infty}D_{m,n}(K_{1},K_{2};U\cap B_{N}(0))=D_{m,n}(K_{1},K_{2};U)\,,\quad\mbox{and}
limN𝔼logDm,n(K1,K2;UBN(0))=𝔼logDm,n(K1,K2;U).\lim_{N\rightarrow\infty}\mathbb{E}\log D_{m,n}(K_{1},K_{2};U\cap B_{N}(0))=\mathbb{E}\log D_{m,n}(K_{1},K_{2};U)\,.

Note that the inequality (2.11) holds for Dm,n(K1,K2;UBN(0))D_{m,n}(K_{1},K_{2};U\cap B_{N}(0)) as long as K1K_{1} and K2K_{2} are connected in UBN(0)U\cap B_{N}(0) which holds for all sufficiently large NN. Therefore, applying (2.11) with UBN(0)U\cap B_{N}(0) instead of UU and then sending NN to infinity gives the desired lemma. ∎

2.4 Percolation with finite range of dependence

In this subsection, we consider the integer lattice d\mathbb{Z}^{d} with d2d\geq 2 and establish some results about percolation with finite range of dependence. The definitions and results from this subsection can be naturally adapted to the rescaled lattice 2nd2^{-n}\mathbb{Z}^{d} for any integer n1n\geq 1. These results will play an important role in Sections 3 and LABEL:sec:bound-distance.

Let M1M\geq 1 be an integer, and consider a probability measure μ\mu on the configuration ω{0,1}d\omega\in\{0,1\}^{\mathbb{Z}^{d}}. We say that μ\mu is MM-dependent if for any two subsets U,VdU,V\subset\mathbb{Z}^{d} with 𝔡(U,V)>M\mathfrak{d}_{\infty}(U,V)>M, the restrictions ω|U\omega|_{U} and ω|V\omega|_{V} are independent. A vertex xx is called open if ω(x)=1\omega(x)=1, and closed if ω(x)=0\omega(x)=0. A path (resp. *-path) is a sequence of vertices x1,,xnx_{1},\ldots,x_{n} such that |xixi+1|1=1|x_{i}-x_{i+1}|_{1}=1 (resp. |xixi+1|=1|x_{i}-x_{i+1}|_{\infty}=1) for any 1in11\leq i\leq{n-1}. A path is called open if all the vertices contained in it are open, and closed if all the vertices contained in it are closed. Similarly, we can define an open *-path and a closed *-path. For a subset UdU\subset\mathbb{Z}^{d}, we use U:={xU:yd\U such that xyE}\partial U:=\{x\in U:\exists y\in\mathbb{Z}^{d}\backslash U\mbox{ such that }xy\in E\} to denote its interior boundary, where EE is the edge set of d\mathbb{Z}^{d}.

We begin with a lemma about the exponential decay of the probability of long closed *-paths when MM is fixed and all the vertices have a probability close to one of being open. The proof follows from an elementary path-counting argument.

Lemma 2.8.

Fix an integer M1M\geq 1. There exist two constants c1(0,1)c_{1}\in(0,1) and C>0C>0 depending only on MM such that for any MM-dependent measure μ\mu satisfying infxdμ[w(x)=1]>c1\inf_{x\in\mathbb{Z}^{d}}\mu[w(x)=1]>c_{1}, we have

μ[There exists a closed -path connecting 0 and BN(0)]CeN/C,N1.\mu\big{[}\mbox{There exists a closed }{\rm*}\mbox{-path connecting }0\mbox{ and }\partial B_{N}(0)\big{]}\leq Ce^{-N/C},\forall N\geq 1.
Proof.

Let p(0,1)p\in(0,1) be a constant to be chosen. Assume that μ\mu is an MM-dependent measure with infxdμ[ω(x)=1]>p\inf_{x\in\mathbb{Z}^{d}}\mu[\omega(x)=1]>p. Let x1,x2,,xnx_{1},x_{2},\ldots,x_{n} be any *-path connecting 0 and BN(0)\partial B_{N}(0). Then, we have

x1=0,xnBN(0),and|xixi+1|=11in1.x_{1}=0\,,\quad x_{n}\in\partial B_{N}(0)\,,\quad\mbox{and}\quad|x_{i}-x_{i+1}|_{\infty}=1\quad\forall 1\leq i\leq n-1\,.

We consider a subset of this path defined inductively as follows: first, take i1=1i_{1}=1, and for j2j\geq 2 define

ij:=max{ij1kn:|xkxij1|M}+1.i_{j}:=\max\{i_{j-1}\leq k\leq n:|x_{k}-x_{i_{j-1}}|_{\infty}\leq M\}+1\,. (2.12)

We stop the induction when ij=n+1i_{j}=n+1. Let us consider the obtained sequence (y1,,ym):=(xi1,,xim)(y_{1},\ldots,y_{m}):=(x_{i_{1}},\ldots,x_{i_{m}}). Then, y1=xi1=0y_{1}=x_{i_{1}}=0. We further have:

min1i<jm|yiyj|>M,max1jm1|yjyj+1|M+1,\displaystyle\min_{1\leq i<j\leq m}|y_{i}-y_{j}|_{\infty}>M\,,\quad\max_{1\leq j\leq m-1}|y_{j}-y_{j+1}|_{\infty}\leq M+1\,, (2.13)
and𝔡(ym,BN(0))M.\displaystyle\quad\mbox{and}\quad\mathfrak{d}_{\infty}(y_{m},\partial B_{N}(0))\leq M\,.

The first property follows directly from (2.12). The second property is because, by (2.12), |yjyj+1|=|xijxij+1||xijxij+11|+1M+1|y_{j}-y_{j+1}|_{\infty}=|x_{i_{j}}-x_{i_{j+1}}|_{\infty}\leq|x_{i_{j}}-x_{i_{j+1}-1}|_{\infty}+1\leq M+1. The last property is because when the iteration stops, we have |ymxn|M|y_{m}-x_{n}|_{\infty}\leq M.

We now upper bound the probability that there exists a closed sequence satisfying (2.13). First, we have mNM+1m\geq\frac{N}{M+1}, which follows from the following inequality:

N1=𝔡(0,BN(0))\displaystyle N-1=\mathfrak{d}_{\infty}(0,\partial B_{N}(0)) i=1m1|yiyi+1|+𝔡(ym,BN(0))\displaystyle\leq\sum_{i=1}^{m-1}|y_{i}-y_{i+1}|_{\infty}+\mathfrak{d}_{\infty}(y_{m},\partial B_{N}(0))
(M+1)(m1)+M.\displaystyle\leq(M+1)(m-1)+M.

For fixed mm, we know that the number of sequences satisfying (2.13) is at most (2M+3)d(m1)(2M+3)^{d(m-1)} since y1=0y_{1}=0 and |yiyi+1|M+1|y_{i}-y_{i+1}|_{\infty}\leq M+1. Furthermore, for a fixed choice of the sequence, the probability that all the vertices contained in it are closed is at most (1p)m(1-p)^{m} by the MM-dependent property of μ\mu and the fact that supxdμ[ω(x)=0]<1p\sup_{x\in\mathbb{Z}^{d}}\mu[\omega(x)=0]<1-p. Therefore, when pp is close enough to one, we have

μ[There exists a closed -path connecting 0 and BN(0)]\displaystyle\quad\mu\big{[}\mbox{There exists a closed }{\rm*}\mbox{-path connecting }0\mbox{ and }\partial B_{N}(0)\big{]}
<mN/(M+1)(2M+3)d(m1)×(1p)mCeN/CN1.\displaystyle<\sum_{m\geq N/(M+1)}(2M+3)^{d(m-1)}\times(1-p)^{m}\leq Ce^{-N/C}\quad\forall N\geq 1\,.\qed

Next, we prove two corollaries from the above lemma, which will be used later. An open (resp. closed) cluster is a connected component of open (resp. closed) vertices. Similarly, we define the open (resp. closed) *-cluster which is a connected component of open (resp. closed) vertices where two vertices x,yx,y are considered to be neighboring each other if |xy|=1|x-y|_{\infty}=1. We define the diameter of a cluster or *-cluster with respect to the ll^{\infty}-distance on d\mathbb{Z}^{d}.

Lemma 2.9.

For an integer M1M\geq 1 and c1=c1(M)c_{1}=c_{1}(M) as defined in Lemma 2.8, let μ\mu be an MM-dependent measure that satisfies infxdμ[w(x)=1]>c1\inf_{x\in\mathbb{Z}^{d}}\mu[w(x)=1]>c_{1}. Then, for each K,N1K,N\geq 1,

μ[Each closed -cluster in BN(0) has diameter at most K]1CNdeK/C.\mu\big{[}\mbox{Each closed }{\rm*}\mbox{-cluster in }B_{N}(0)\mbox{ has diameter at most K}\big{]}\geq 1-CN^{d}e^{-K/C}.

Here, the constant CC may depend on MM, but is independent of μ\mu.

Proof.

If there exists a closed *-cluster in BN(0)B_{N}(0) with diameter at least K+1K+1, then we can find a vertex xBN(0)x\in B_{N}(0) such that xx is connected to BK+1(x)\partial B_{K+1}(x) with a closed *-path. Summing over all the possible choices of xx and applying Lemma 2.8 using K+1K+1 instead of NN, we obtain the desired result. ∎

Lemma 2.10.

For any integer M1M\geq 1 and ε>0\varepsilon>0, there exists a constant c2=c2(M,ε)(0,1)c_{2}=c_{2}(M,\varepsilon)\in(0,1) such that for any MM-dependent measure μ\mu satisfying infxdμ[w(x)=1]>c2\inf_{x\in\mathbb{Z}^{d}}\mu[w(x)=1]>c_{2}, we have

μ[There exists an infinite open cluster containing 0]1ε.\mu\big{[}\mbox{There exists an infinite open cluster containing }0\big{]}\geq 1-\varepsilon\,.
Proof.

Recall from Lemma 2.8 the constant c1c_{1}, which depends on MM. Let p(c1,1)p\in(c_{1},1) be a constant to be chosen. Let μ\mu be an MM-dependent measure with infxdμ[ω(x)=1]>p\inf_{x\in\mathbb{Z}^{d}}\mu[\omega(x)=1]>p. Let NN be a large integer to be chosen. Define the events

𝒦1\displaystyle\mathcal{K}_{1} :={All vertices in BN(0) are open},\displaystyle:=\{\mbox{All vertices in }B_{N}(0)\mbox{ are open}\}\,,
𝒦2\displaystyle\mathcal{K}_{2} :={There exists a closed -cluster enclosing BN(0)}.\displaystyle:=\{\mbox{There exists a closed }{\rm*}\mbox{-cluster enclosing }B_{N}(0)\}\,.

By duality, we know that on the event 𝒦1\𝒦2\mathcal{K}_{1}\backslash\mathcal{K}_{2}, all vertices in BN(0)B_{N}(0) are open and are connected to infinity by an open path. Hence, we only need to show that

μ[𝒦1\𝒦2]1ε.\mu[\mathcal{K}_{1}\backslash\mathcal{K}_{2}]\geq 1-\varepsilon. (2.14)

First, we prove a lower bound for the probability of 𝒦1\mathcal{K}_{1}. Using the assumption that infxdμ[ω(x)=1]>p\inf_{x\in\mathbb{Z}^{d}}\mu[\omega(x)=1]>p, we obtain

μ[𝒦1]1xBN(0)μ[ω(x)=0]1(2N+1)d(1p).\mu[\mathcal{K}_{1}]\geq 1-\sum_{x\in B_{N}(0)}\mu[\omega(x)=0]\geq 1-(2N+1)^{d}(1-p)\,. (2.15)

Next, we establish an upper bound for the probability of 𝒦2\mathcal{K}_{2}. If the event 𝒦2\mathcal{K}_{2} happens, then the closed *-cluster must intersect the set {xd:x1N+1,x2=x3==xd=0}\{x\in\mathbb{Z}^{d}:x_{1}\geq N+1,x_{2}=x_{3}=\ldots=x_{d}=0\}. Let x=(m,0,,0)x=(m,0,\ldots,0) be an intersection point where mN+1m\geq N+1. Then there exists a closed *-path from xx to Bm(x)\partial B_{m}(x). Thus, by applying Lemma 2.8 using mm instead of NN, we obtain:

μ[𝒦2]\displaystyle\mu[\mathcal{K}_{2}] mN+1μ[There exists a closed -path from x to Bm(x)]\displaystyle\leq\sum_{m\geq N+1}\mu\big{[}\mbox{There exists a closed }{\rm*}\mbox{-path from }x\mbox{ to }\partial B_{m}(x)\big{]} (2.16)
mN+1Cem/CCeN/C.\displaystyle\leq\sum_{m\geq N+1}Ce^{-m/C}\leq Ce^{-N/C}.

Combining (2.15) and (2.16), and first taking NN to be large and then taking pp close to 11, yields (2.14). In particular, the choice of pp depends only on MM and ε\varepsilon. This concludes the lemma. ∎

3 Existence of an exponent

In this section, we will first prove the existence of an exponent Q=Q(ξ)Q=Q(\xi)\in\mathbb{R} such that (1.8) holds (Proposition 3.1). This exponent governs the internal DnD_{n}-distance between two points in a box as nn grows. Furthermore, Lemma LABEL:lem:distance-any-point extends this result to any pair of points, and Lemma LABEL:lem:Q-lower establishes some basic properties about Q(ξ)Q(\xi). Combining these results gives Proposition 1.1.

We first introduce some notations. For each 1id1\leq i\leq d, let

ei:=the i-th standard basis vector in d.e_{i}:=\mbox{the }i\mbox{-th standard basis vector in }\mathbb{R}^{d}. (3.1)

That is, eie_{i} is a {0,1}\{0,1\} valued vector in d\mathbb{R}^{d} where only the ii-th coordinate is equal to 11. For integer n1n\geq 1 and p(0,1)p\in(0,1), let an(p)a_{n}^{(p)} represent the pp-th quantile of the internal distance Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)), defined as

an(p):=inf{l>0:[Dn(0,e1;B2(0))l]>p}.a_{n}^{(p)}:=\inf\{l>0:\mathbb{P}[D_{n}(0,e_{1};B_{2}(0))\leq l]>p\}\,. (3.2)

We have [Dn(0,e1;B2(0))l]=p\mathbb{P}[D_{n}(0,e_{1};B_{2}(0))\leq l]=p since Dn(0,e1;B2(0))D_{n}(0,e_{1};B_{2}(0)) is a continuous random variable. When p=1/2p=1/2, the number λn\lambda_{n} from (1.7) satisfies λn=an(1/2)=Med(Dn(0,e1;B2(0)))\lambda_{n}=a_{n}^{(1/2)}={\rm Med}(D_{n}(0,e_{1};B_{2}(0))).

Proposition 3.1.

There exists an exponent Q=Q(ξ)Q=Q(\xi)\in\mathbb{R} such that

λn=2(1ξQ)n+o(n)as n.\lambda_{n}=2^{-(1-\xi Q)n+o(n)}\quad\mbox{as }n\rightarrow\infty\,. (3.3)

The proof of Proposition 3.1 is via a subadditivity argument. We will use Lemmas 3.2 and 3.3 below. The former directly follows from the concentration bound in Lemma 2.7. The latter employs a percolation argument from Subsection 2.4 and follows an approach similar to that of [dg-supercritical-lfpp, Lemma 2.9].

Lemma 3.2.

For fixed p1,p2(0,1)p_{1},p_{2}\in(0,1), there exists a constant C>0C>0 depending only on p1p_{1} and p2p_{2} such that for all integer n1n\geq 1, we have

eCnan(p1)an(p2)eCnan(p1).e^{-C\sqrt{n}}a_{n}^{(p_{1})}\leq a_{n}^{(p_{2})}\leq e^{C\sqrt{n}}a_{n}^{(p_{1})}.
Proof.

Applying the concentration bound from Lemma 2.7 with K1={0},K2={e1}K_{1}=\{0\},K_{2}=\{e_{1}\}, and U=B2(0)U=B_{2}(0), yields that

[|logDn(0,e1;B2(0))𝔼logDn(0,e1;B2(0))|t]Cet2Cnt>0.\mathbb{P}\big{[}|\log D_{n}(0,e_{1};B_{2}(0))-\mathbb{E}\log D_{n}(0,e_{1};B_{2}(0))|\geq t\big{]}\leq Ce^{-\frac{t^{2}}{Cn}}\quad\forall t>0\,.

Hence, for any fixed p(0,1)p\in(0,1), the following inequality holds:

|logan(p)𝔼logDn(0,e1;B2(0))|Cn,|\log a_{n}^{(p)}-\mathbb{E}\log D_{n}(0,e_{1};B_{2}(0))|\leq C\sqrt{n}\,,

where the constant CC depends on pp, but is independent of nn. This implies the lemma. ∎

We now present a key lemma. It will imply Proposition 3.1 when combined with Lemma 3.2.

Lemma 3.3.

There exist c3(0,1)c_{3}\in(0,1) and a constant C>0C>0 such that for all integers n>m1n>m\geq 1:

λneCn2/3am(c3)anm(c3).\lambda_{n}\leq e^{Cn^{2/3}}a_{m}^{(c_{3})}a_{n-m}^{(c_{3})}\,. (3.4)

We will first use a subadditivity argument to prove Proposition 3.1 based on this lemma, and then provide the proof of Lemma 3.3.

Proof of Proposition 3.1.

Combining Lemmas 3.3 and 3.2, we obtain that for all n>m1n>m\geq 1:

λneCn2/3am(c3)anm(c3)eCn2/3λmλnm.\lambda_{n}\leq e^{Cn^{2/3}}a_{m}^{(c_{3})}a_{n-m}^{(c_{3})}\leq e^{Cn^{2/3}}\lambda_{m}\lambda_{n-m}\,.

Combining this inequality with Lemma 6.4.10 in [dembo-ld], applied to logλn\log\lambda_{n}, implies the existence of aa\in\mathbb{R} such that:

λn=ean+o(n)as n.\lambda_{n}=e^{an+o(n)}\quad\mbox{as }n\rightarrow\infty\,.

Taking QQ\in\mathbb{R} such that ea=2(1ξQ)e^{a}=2^{-(1-\xi Q)} yields the desired result. ∎

Next, we proceed to the proof of Lemma 3.3. First, we present two auxiliary results. In Lemma 3.4, we provide estimates for the field hnm,nh_{n-m,n}. Subsequently, we use these estimates in Lemma 3.5 to compare an(p)a_{n}^{(p)} and anm(p)a_{n-m}^{(p)}.

Lemma 3.4.

There exist constants C1>0C_{1}>0 and C>0C>0 such that for all integers n>m1n>m\geq 1:

[supxB2(0)hnm,n(x)C1mn]Cen/C.\mathbb{P}\Big{[}\sup_{x\in B_{2}(0)}h_{n-m,n}(x)\geq C_{1}\sqrt{mn}\Big{]}\leq Ce^{-n/C}.
Proof.

By Claim (4) in Lemma 2.2, we have that supyB2mn(0)hnm,n(y)=dsupxB1(0)h0,m(x)\sup_{y\in B_{2^{m-n}}(0)}h_{n-m,n}(y)\overset{d}{=}\sup_{x\in B_{1}(0)}h_{0,m}(x). Therefore,

[supyB2mn(0)hnm,n(y)s]=[supxB1(0)hm(x)s]s>0.\mathbb{P}\big{[}\sup_{y\in B_{2^{m-n}}(0)}h_{n-m,n}(y)\geq s\big{]}=\mathbb{P}\big{[}\sup_{x\in B_{1}(0)}h_{m}(x)\geq s\big{]}\quad\forall s>0\,.

Using Claims (3) and (4) from Lemma 2.3, we get that for all s(1+2dlog2)ms\geq(1+\sqrt{2d}\log 2)m:

[supyB2mn(0)hnm,n(y)s]Cexp(s2Cm).\mathbb{P}\big{[}\sup_{y\in B_{2^{m-n}}(0)}h_{n-m,n}(y)\geq s\big{]}\leq C\exp\big{(}-\frac{s^{2}}{Cm})\,. (3.5)

Hence, for all t>1+2dlog2t>1+\sqrt{2d}\log 2:

[supxB2(0)hnm,n(x)tmn]\displaystyle\quad\mathbb{P}\Big{[}\sup_{x\in B_{2}(0)}h_{n-m,n}(x)\geq t\sqrt{mn}\Big{]}
=[supxB2(0)2mndsupyB2mn(x)hnm,n(y)tmn]\displaystyle=\mathbb{P}\Big{[}\sup_{x\in B_{2}(0)\cap 2^{m-n}\mathbb{Z}^{d}}\sup_{y\in B_{2^{m-n}}(x)}h_{n-m,n}(y)\geq t\sqrt{mn}\Big{]}
xB2(0)2mnd[supyB2mn(0)hnm,n(y)tmn]C2nd×exp(t2mnCm).\displaystyle\leq\sum_{x\in B_{2}(0)\cap 2^{m-n}\mathbb{Z}^{d}}\mathbb{P}\Big{[}\sup_{y\in B_{2^{m-n}}(0)}h_{n-m,n}(y)\geq t\sqrt{mn}\Big{]}\leq C2^{nd}\times\exp(-\frac{t^{2}mn}{Cm})\,.

In the last inequality, we used (3.5), as well as the facts that n>mn>m and |B2(0)2mnd|C2nd|B_{2}(0)\cap 2^{m-n}\mathbb{Z}^{d}|\leq C2^{nd}. By choosing a sufficiently large tt, we obtain the desired result. ∎

We now provide a comparison between an(p)a_{n}^{(p)} and anm(p)a_{n-m}^{(p)} based on the above lemma.

Lemma 3.5.

For a fixed p(0,1)p\in(0,1), there exists a constant C=C(p)>0C=C(p)>0 such that for all integers n>m1n>m\geq 1:

eCmnanm(p)an(p)eCmnanm(p).e^{-C\sqrt{mn}}a_{n-m}^{(p)}\leq a_{n}^{(p)}\leq e^{C\sqrt{mn}}a_{n-m}^{(p)}\,.
Proof.

Based on the definition of an(p)a_{n}^{(p)} given by (3.2), we obtain:

[Dnm(0,e1;B2(0))anm(p/2)]=1p/2.\mathbb{P}\big{[}D_{n-m}(0,e_{1};B_{2}(0))\geq a_{n-m}^{(p/2)}\big{]}=1-p/2\,. (3.6)

By using Lemma 3.4 and the symmetry of hm,nh_{m,n}, there exists a constant A>0A>0 such that for all n>m1n>m\geq 1:

[infxB2(0)hnm,n(x)Amn]>1p/2.\mathbb{P}\big{[}\inf_{x\in B_{2}(0)}h_{n-m,n}(x)\geq-A\sqrt{mn}\big{]}>1-p/2\,. (3.7)

Since hn=hnm+hnm,nh_{n}=h_{n-m}+h_{n-m,n}, we have

Dn(0,e1;B2(0))Dnm(0,e1;B2(0))eξinfxB2(0)hnm,n(x).D_{n}(0,e_{1};B_{2}(0))\geq D_{n-m}(0,e_{1};B_{2}(0))e^{\xi\inf_{x\in B_{2}(0)}h_{n-m,n}(x)}.

Therefore, for all s>0s>0:

[Dn(0,e1;B2(0))esmnanm(p/2)]\displaystyle\quad\mathbb{P}\big{[}D_{n}(0,e_{1};B_{2}(0))\geq e^{-s\sqrt{mn}}a_{n-m}^{(p/2)}\big{]}
[{Dnm(0,e1;B2(0))anm(p/2)}{infxB2(0)hnm,n(x)smn/ξ}].\displaystyle\geq\mathbb{P}\Big{[}\big{\{}D_{n-m}(0,e_{1};B_{2}(0))\geq a_{n-m}^{(p/2)}\big{\}}\cap\big{\{}\inf_{x\in B_{2}(0)}h_{n-m,n}(x)\geq-s\sqrt{mn}/\xi\big{\}}\Big{]}\,.

Combining this with (3.6) and (3.7), we obtain that for all s>Aξs>A\xi, with AA being the constant from (3.7),

[Dn(0,e1;B2(0))esmnanm(p/2)]>1p/2p/2=1p.\displaystyle\mathbb{P}\big{[}D_{n}(0,e_{1};B_{2}(0))\geq e^{-s\sqrt{mn}}a_{n-m}^{(p/2)}\big{]}>1-p/2-p/2=1-p\,.

Combining with the definition of an(p)a_{n}^{(p)} and Lemma 3.2 yields that

an(p)e(Aξ+1)mnanm(p/2)eCmnanm(p).a_{n}^{(p)}\geq e^{-(A\xi+1)\sqrt{mn}}a_{n-m}^{(p/2)}\geq e^{-C\sqrt{mn}}a_{n-m}^{(p)}. (3.8)

Similarly, for sufficiently large s>0s>0, we can show that

[Dn(0,e1;B2(0))esmnanm(p+12)]\displaystyle\mathbb{P}\big{[}D_{n}(0,e_{1};B_{2}(0))\leq e^{s\sqrt{mn}}a_{n-m}^{(\frac{p+1}{2})}\big{]}
[{Dnm(0,e1;B2(0))anm(p+12)}{supxB2(0)hnm,n(x)smn/ξ}]\displaystyle\qquad\geq\mathbb{P}\Big{[}\big{\{}D_{n-m}(0,e_{1};B_{2}(0))\leq a_{n-m}^{(\frac{p+1}{2})}\big{\}}\cap\big{\{}\sup_{x\in B_{2}(0)}h_{n-m,n}(x)\leq s\sqrt{mn}/\xi\big{\}}\Big{]}
>p+121p2=p.\displaystyle\qquad>\frac{p+1}{2}-\frac{1-p}{2}=p\,.

This, together with the definition of an(p)a_{n}^{(p)} and Lemma 3.2, implies that

an(p)eCmnanm(p+12)eCmnanm(p).a_{n}^{(p)}\leq e^{C\sqrt{mn}}a_{n-m}^{(\frac{p+1}{2})}\leq e^{C\sqrt{mn}}a_{n-m}^{(p)}. (3.9)

Combining (3.8) and (3.9) yields the desired result. ∎

We now turn to the proof of Lemma 3.3. The proof follows a similar approach to that of [dg-supercritical-lfpp, Lemma 2.9]. Our goal is to construct a path that connects 0 and e1e_{1} within the box B2(0)B_{2}(0), such that the DnD_{n}-length of this path can be upper-bounded by am(p)a_{m}^{(p)} and anm(p)a_{n-m}^{(p)} with high probability provided that pp is sufficiently large. (We will actually use anmk(p)a_{n-m-k}^{(p)}, with k=(logm)2k=\lfloor(\log m)^{2}\rfloor, instead of anm(p)a_{n-m}^{(p)}. However, by Lemma 3.5, they do not differ much.) The construction will consist of four steps. In Step 1, we introduce some regularity events for the field, which all happen with high probability. In Step 2, we construct a discrete path on Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} (recall its definition from (2.3)) whose DmD_{m}-length can be upper-bounded by am(p)a_{m}^{(p)}. Step 3 involves local modifications to the discrete path so that its Dm+k,nD_{m+k,n}-length can be upper-bounded. We will use a percolation argument for the rescaled lattice 2mkd2^{-m-k}\mathbb{Z}^{d} to achieve this. The introduction of the auxiliary scale k=(logm)2k=\lfloor(\log m)^{2}\rfloor is mainly for this step. In Step 4, we control the DnD_{n}-length of the resulting path using the regularity events.

Proof of Lemma 3.3.

Let p(0,1)p\in(0,1) be a constant to be chosen. Define the integer

k:=(logm)2.k:=\lfloor(\log m)^{2}\rfloor\,. (3.10)

We assume that

m>100andn>m+k.m>100\quad\mbox{and}\quad n>m+k\,.

Otherwise, Equation (3.4) can be deduced from Lemmas 3.2 and 3.5 by choosing a sufficiently large CC. This is because, for a fixed pp, by Lemmas 3.2 and 3.5, we have

λneCnan(p)eCnanm(p)eCnam(p)anm(p)1m100,\displaystyle\lambda_{n}\leq e^{C\sqrt{n}}a_{n}^{(p)}\leq e^{C\sqrt{n}}a_{n-m}^{(p)}\leq e^{C\sqrt{n}}a_{m}^{(p)}a_{n-m}^{(p)}\quad\forall 1\leq m\leq 100\,, (3.11)
λneCnan(p)a1(p)eCnkam(p)anm(p)m<nm+k,\displaystyle\lambda_{n}\leq e^{C\sqrt{n}}a_{n}^{(p)}a_{1}^{(p)}\leq e^{C\sqrt{nk}}a_{m}^{(p)}a_{n-m}^{(p)}\quad\forall m<n\leq m+k\,,

and Cn,CnkCn2/3C\sqrt{n},C\sqrt{nk}\leq Cn^{2/3}.

Next, we will construct a path connecting 0 and e1e_{1} within B2(0)B_{2}(0). When pp is sufficiently close to one (not depending on nn), the DnD_{n}-length of this path will be at most eCn2/3am(p)anmk(p)e^{Cn^{2/3}}a_{m}^{(p)}a_{n-m-k}^{(p)} with probability at least 1/21/2. Therefore,

λneCn2/3am(p)anmk(p).\lambda_{n}\leq e^{Cn^{2/3}}a_{m}^{(p)}a_{n-m-k}^{(p)}\,. (3.12)

Combining this with Lemma 3.5, we obtain Lemma 3.3.

As announced earlier, the construction consists of four steps:

Step 1: Regularity event for hmh_{m} and hm,m+kh_{m,m+k}. Define the event

1:={2msupxB2(0)|hm(x)|n2/3}{supxB2(0)hm,m+k(x)C1k(m+k)},\mathcal{E}_{1}:=\big{\{}2^{-m}\sup_{x\in B_{2}(0)}|\nabla h_{m}(x)|_{\infty}\leq n^{2/3}\big{\}}\cap\big{\{}\sup_{x\in B_{2}(0)}h_{m,m+k}(x)\leq C_{1}\sqrt{k(m+k)}\big{\}}\,, (3.13)

where C1C_{1} is the constant defined in Lemma 3.4. Using the fact that |2mdB2(0)|C2md|2^{-m}\mathbb{Z}^{d}\cap B_{2}(0)|\leq C2^{md} and Claim (2) in Lemma 2.3, we obtain

[2msupxB2(0)|hm(x)|n2/3]\displaystyle\mathbb{P}\Big{[}2^{-m}\sup_{x\in B_{2}(0)}|\nabla h_{m}(x)|_{\infty}\leq n^{2/3}\Big{]} (3.14)
1x2mdB2(0)[supyB2m(x)|hm(y)|>2mn2/3]\displaystyle\qquad\geq 1-\sum_{x\in 2^{-m}\mathbb{Z}^{d}\cap B_{2}(0)}\mathbb{P}\Big{[}\sup_{y\in B_{2^{-m}}(x)}|\nabla h_{m}(y)|_{\infty}>2^{m}n^{2/3}\Big{]}
1C2md×Cen4/3/C1Cen4/3/C.\displaystyle\qquad\geq 1-C2^{md}\times Ce^{-n^{4/3}/C}\geq 1-Ce^{-n^{4/3}/C}.

Combining this with Lemma 3.4, applied with (m+k,k)(m+k,k) instead of (n,m)(n,m), yields that

[1]1Cen4/3/CCem/C1Cem/C.\mathbb{P}[\mathcal{E}_{1}]\geq 1-Ce^{-n^{4/3}/C}-Ce^{-m/C}\geq 1-Ce^{-m/C}. (3.15)

Step 2: Discretize the DmD_{m}-geodesic between 0 and e1e_{1} on Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}. Define the event

2:={Dm(0,e1;B2(0))am(p)}.\mathcal{E}_{2}:=\{D_{m}(0,e_{1};B_{2}(0))\leq a_{m}^{(p)}\}\,. (3.16)

By (3.2), we have

[2]=p.\mathbb{P}[\mathcal{E}_{2}]=p\,. (3.17)

On the event 2\mathcal{E}_{2}, there exists a piecewise continuously differentiable path P:[0,1]B2(0)P:[0,1]\rightarrow B_{2}(0) from 0 to e1e_{1} such that

len(P;Dm)=01eξhm(P(t))|P(t)|dt2am(p).{\rm len}(P;D_{m})=\int_{0}^{1}e^{\xi h_{m}(P(t))}|P^{\prime}(t)|dt\leq 2a_{m}^{(p)}. (3.18)

Recall from (2.3) that Lm=2mdB2(0){\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}=2^{-m}\mathbb{Z}^{d}\cap B_{2}(0). Then, we have 0,e1Lm0,e_{1}\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}. For an illustration, we refer to Figure 1. We consider Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} as a subset of d\mathbb{R}^{d}. Sometimes, we will consider (*-)paths or (*-)clusters on the rescaled lattice 2md2^{-m}\mathbb{Z}^{d}, as defined in Subsection 2.4, and only in these cases, we regard Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} as a subset of 2md2^{-m}\mathbb{Z}^{d}. We now construct, on the event 12\mathcal{E}_{1}\cap\mathcal{E}_{2}, a self-avoiding path on Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} as a discrete approximation of the path PP. See Figure 1 for an illustration.777For illustrative purposes, we depict planar graphs, but all these arguments hold for dimensions greater than two.

Refer to caption
Figure 1: Illustration of the sets Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} and Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime}, and the path PP and its discrete approximation: (x1,,xJ)(x_{1},\ldots,x_{J}). The dotted red lines represent the edges between neighboring vertices in Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime}. The path (x1,,xJ)(x_{1},\ldots,x_{J}), as illustrated by the red curve, is a self-avoiding path on Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime} connecting 0 and e1e_{1}.

Let Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime} be a subset of Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} defined as follows:

Lm:={xLm:PB2m(x)¯},{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime}:=\{x\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}:P\cap\overline{B_{2^{-m}}(x)}\neq\emptyset\}\,,

where B2m(x)¯\overline{B_{2^{-m}}(x)} represents the closure of B2m(x)B_{2^{-m}}(x). It follows that 0,e1Lm0,e_{1}\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime}, and there exists a discrete path in Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime} connecting them. This is because for any xLmx\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime} considering the first exit time of PP from the box B2m(x)¯\overline{B_{2^{-m}}(x)}, we can find a vertex yLmy\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} such that |xy|1=2m|x-y|_{1}=2^{-m}, and PP also enters the box B2m(y)¯\overline{B_{2^{-m}}(y)}. By doing this procedure iteratively, we obtain a discrete path in Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime} that connects 0 and e1e_{1}. Taking any path in Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime} connecting 0 and e1e_{1}, and applying the loop erasure procedure similar to (2.12), yields a self-avoiding path connecting 0 and e1e_{1} in Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}^{\prime}. That is, there exists a self-avoiding path 0=x1,,xJ=e10=x_{1},\ldots,x_{J}=e_{1} satisfying the properties that

xiLmandPB2m(xi)¯1iJ,and\displaystyle x_{i}\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}\quad\mbox{and}\quad P\cap\overline{B_{2^{-m}}(x_{i})}\neq\emptyset\quad\forall 1\leq i\leq J\,,\quad\mbox{and} (3.19)
|xixi+1|1=2m1iJ1.\displaystyle|x_{i}-x_{i+1}|_{1}=2^{-m}\quad\forall 1\leq i\leq J-1\,.

We now show that on the event 12\mathcal{E}_{1}\cap\mathcal{E}_{2}, we have

j=1J2meξh0,m(xj)am(p)eCn2/3.\sum_{j=1}^{J}2^{-m}e^{\xi h_{0,m}(x_{j})}\leq a_{m}^{(p)}e^{Cn^{2/3}}. (3.20)

This is because for each 1iJ1\leq i\leq J, the second property PB2m(xi)¯P\cap\overline{B_{2^{-m}}(x_{i})}\neq\emptyset in (3.19) ensures that the path PP must cross the hypercubic shell B2m+1(xi)\B2m(xi)B_{2^{-m+1}}(x_{i})\backslash B_{2^{-m}}(x_{i}). This segment has Euclidean length of at least 2m2^{-m}. By the event 1\mathcal{E}_{1}, for some C>0C>0, we have:

infzB2m+1(xi)h0,m(z)h0,m(xi)Cn2/3.\inf_{z\in B_{2^{-m+1}}(x_{i})}h_{0,m}(z)\geq h_{0,m}(x_{i})-Cn^{2/3}.

Therefore, this segment has a DmD_{m}-length of at least

2meξh0,m(xi)Cn2/3.2^{-m}e^{\xi h_{0,m}(x_{i})-Cn^{2/3}}. (3.21)

Furthermore, each point on PP is contained in at most 5d5^{d} such hypercubic shells. Combining this fact with (3.21) and (3.18), we obtain (3.20).

Step 3: Modify the path on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}. Recall from (2.3) that Lm+k=2mkdB2(0){\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}=2^{-m-k}\mathbb{Z}^{d}\cap B_{2}(0). It follows that LmLm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}\subset{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}. We now construct a path on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} that closely follows the path (x1,,xJ)(x_{1},\ldots,x_{J}) and has typical Dm+k,nD_{m+k,n}-length. We call a vertex xLm+kx\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} open if for all σ{1,1}\sigma\in\{1,-1\} and 1id1\leq i\leq d

Dm+k,n(x,x+σei2mk;B2mk+1(x))2mkanmk(p),D_{m+k,n}(x,x+\sigma e_{i}2^{-m-k};B_{2^{-m-k+1}}(x))\leq 2^{-m-k}a_{n-m-k}^{(p)}, (3.22)

and closed otherwise. We assume that all the vertices in 2mkd\Lm+k2^{-m-k}\mathbb{Z}^{d}\backslash{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} are open. Using the translation and rotational invariance and the scaling property from Lemma 2.5, we have

Dm+k,n(x,x+σei2mk;B2mk+1(x))=d2mkDnmk(0,e1;B2(0)).D_{m+k,n}(x,x+\sigma e_{i}2^{-m-k};B_{2^{-m-k+1}}(x))\overset{d}{=}2^{-m-k}D_{n-m-k}(0,e_{1};B_{2}(0))\,.

Combining this with the definition of anmk(p)a_{n-m-k}^{(p)} from (3.2), we obtain that for all xLm+kx\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}

[x is open]\displaystyle\quad\mathbb{P}[x\mbox{ is open}] (3.23)
1σ=±1,1id[Dm+k,n(x,x+σei2mk;B2mk+1(x))>2mkanmk(p)]\displaystyle\geq 1-\sum_{\sigma=\pm 1,1\leq i\leq d}\mathbb{P}\big{[}D_{m+k,n}(x,x+\sigma e_{i}2^{-m-k};B_{2^{-m-k+1}}(x))>2^{-m-k}a_{n-m-k}^{(p)}\big{]}
=12d[D0,nmk(0,e1;B2(0))>anmk(p)]=12d(1p).\displaystyle=1-2d\cdot\mathbb{P}\big{[}D_{0,n-m-k}(0,e_{1};B_{2}(0))>a_{n-m-k}^{(p)}\big{]}=1-2d(1-p)\,.

In particular, as pp approaches one, this probability also tends to one. Recalling the notation in Subsection 2.4, we similarly define open (or closed) (*-)paths and (*-)clusters on the rescaled lattice 2mkd2^{-m-k}\mathbb{Z}^{d}. Define the event

3\displaystyle\mathcal{E}_{3} :={Both 0 and e1 are contained in infinite open clusters on 2mkd,\displaystyle:=\{\mbox{Both }0\mbox{ and }e_{1}\mbox{ are contained in infinite open clusters on }2^{-m-k}\mathbb{Z}^{d}, (3.24)
 and each closed -cluster has diameter at most 2k2}.\displaystyle\qquad\mbox{ and each closed }{\rm*}\mbox{-cluster has diameter at most }2^{k-2}\}\,.

Here, the diameter is associated with the graph distance on the rescaled lattice 2mkd2^{-m-k}\mathbb{Z}^{d}.

By the definition in (3.22), whether a vertex xx is open is determined by the field hm+k,nh_{m+k,n} restricted to the domain B2mk+1(x)B_{2^{-m-k+1}}(x). So, according to Lemma 2.6, for two subsets U,V2mkdU,V\subset 2^{-m-k}\mathbb{Z}^{d} with graph distance at least 2r0+42{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}+4, the statuses of the vertices in UU being open or closed are independent of the statuses of those within VV. Therefore, \mathbb{P} induces an MM-dependent measure on {0,1}2mkd\{0,1\}^{2^{-m-k}\mathbb{Z}^{d}} (where 0 represents closed and 11 represents open) with M=2r0+4+1M=\lfloor 2{\hyperref@@ii[K-condition2]{\mathfrak{r}_{0}}}+4\rfloor+1. As a result, we can apply the percolation result in Subsection 2.4. By using (3.23), Lemma 2.9, and Lemma 2.10 (with ε=0.01\varepsilon=0.01), we can show the existence of c2(0,1)c_{2}^{\prime}\in(0,1) such that when pc2p\geq c_{2}^{\prime}, the following inequality holds:

[3]12×0.01C2d(m+k)e2k2/C10.02Ce2k/C.\mathbb{P}[\mathcal{E}_{3}]\geq 1-2\times 0.01-C2^{d(m+k)}e^{-2^{k-2}/C}\geq 1-0.02-Ce^{-2^{k}/C}. (3.25)

The last inequality is due to the fact that k(logm)21k\geq(\log m)^{2}-1. From now on, we take

p=max{c2,0.99}.p=\max\{c_{2}^{\prime},0.99\}. (3.26)

On the event 3\mathcal{E}_{3}, for each xLmLm+kx\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m}\subset{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}, there is no closed *-cluster on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} that crosses the hypercubic shell B2m(x)\B2m1(x)B_{2^{-m}}(x)\backslash B_{2^{-m-1}}(x)888For xLmx\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} with 𝔡(x,B2(0))=2m\mathfrak{d}_{\infty}(x,\partial B_{2}(0))=2^{-m}, we consider the hypercubic shell (B2m(x)\B2m1(x)){y:𝔡(y,B2(0))>2mk}(B_{2^{-m}}(x)\backslash B_{2^{-m-1}}(x))\cap\{y:\mathfrak{d}_{\infty}(y,\partial B_{2}(0))>2^{-m-k}\} instead. This ensures that for any yLm+ky\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} under consideration, we have B2mk+1(y)B2(0)B_{2^{-m-k+1}}(y)\subset B_{2}(0). or encloses B2m1(x)B_{2^{-m-1}}(x). Therefore, by duality, there exists a unique open cluster on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} that encloses B2m1(x)B_{2^{-m-1}}(x) within the hypercubic shell B2m(x)\B2m1(x)B_{2^{-m}}(x)\backslash B_{2^{-m-1}}(x). Furthermore, the open clusters corresponding to neighboring vertices on Lm{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m} intersect, as illustrated in Figure 2. Since both 0 and e1e_{1} are contained in infinite open clusters, we can find open paths on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}  see the brown curves in Figure 2. These paths connect 0 and e1e_{1} to their corresponding open clusters that enclose B2m1(0)B_{2^{-m-1}}(0) or B2m1(e1)B_{2^{-m-1}}(e_{1}), respectively.

By joining these open paths and clusters together, and applying the loop erasure procedure, we can construct a self-avoiding open path on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} that connects 0 and e1e_{1}, closely following the sequence (x1,,xJ)(x_{1},\ldots,x_{J}). Let us denote the resulting path as 0=y1,y2,,yK=e10=y_{1},y_{2},\ldots,y_{K}=e_{1}. It satisfies the condition that for each 1iK1\leq i\leq K:

yiLm+k is open,andmin1lJ|yixl|2m.y_{i}\in{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k}\mbox{ is open}\,,\quad\mbox{and}\quad\min_{1\leq l\leq J}|y_{i}-x_{l}|_{\infty}\leq 2^{-m}. (3.27)
Refer to caption
Figure 2: The red path corresponds to (x1,x2,,xJ)(x_{1},x_{2},\ldots,x_{J}). The open clusters on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} that enclose B2m1(xi)B_{2^{-m-1}}(x_{i}) are depicted in purple, and the two brown curves represent the open paths that connect 0 and e1e_{1} to their corresponding open clusters. By joining these clusters, we can construct a path on Lm+k{\hyperref@@ii[eq:def-rescaled-lattice]{\mathscr{L}}}_{m+k} connecting 0 and e1e_{1} that closely follows the red curve and has typical Dm+k,nD_{m+k,n}-length.

Step 4: Concatenate the geodesic and upper bound the DnD_{n}-length. In the final step, we join the geodesics between yiy_{i} and yi+1y_{i+1} for 1iK11\leq i\leq K-1 and upper bound its DnD_{n}-length. Assume that

123 happens.\mathcal{E}_{1}\cap\mathcal{E}_{2}\cap\mathcal{E}_{3}\mbox{ happens}\,.

By using (3.27) and the definition of open vertices from (3.22), for each 1iK11\leq i\leq K-1, there exists a piecewise continuously differentiable path Pi:[0,1]B2mk+1(yi)P_{i}:[0,1]\rightarrow B_{2^{-m-k+1}}(y_{i}) that connects yiy_{i} and yi+1y_{i+1} and satisfies:

len(Pi;Dm+k,n)=01eξhm+k,n(Pi(t))|Pi(t)|dt21mkanmk(p).{\rm len}(P_{i};D_{m+k,n})=\int_{0}^{1}e^{\xi h_{m+k,n}(P_{i}(t))}|P_{i}^{\prime}(t)|dt\leq 2^{1-m-k}a_{n-m-k}^{(p)}. (3.28)

By concatenating the paths P1,P2,,PK1P_{1},P_{2},\ldots,P_{K-1}, we obtain a path P~\widetilde{P} that connects 0 and e1e_{1} within B2(0)B_{2}(0).

We now upper bound the DnD_{n}-length of P~\widetilde{P} on the event 123\mathcal{E}_{1}\cap\mathcal{E}_{2}\cap\mathcal{E}_{3}. For each 1iK1\leq i\leq K, by (3.27), we can choose 1kiJ1\leq k_{i}\leq J such that

|yixki|2m.|y_{i}-x_{k_{i}}|_{\infty}\leq 2^{-m}. (3.29)

Since h0,n=h0,m+hm,m+k+hm+k,nh_{0,n}=h_{0,m}+h_{m,m+k}+h_{m+k,n}, we have

len(P~;Dn)\displaystyle{\rm len}(\widetilde{P};D_{n}) =i=1K101eξh0,n(Pi(t))|Pi(t)|dt\displaystyle=\sum_{i=1}^{K-1}\int_{0}^{1}e^{\xi h_{0,n}(P_{i}(t))}|P_{i}^{\prime}(t)|dt (3.30)
=i=1K101eξh0,m(Pi(t))eξhm,m+k(Pi(t))eξhm+k,n(Pi(t))|Pi(t)|dt.\displaystyle=\sum_{i=1}^{K-1}\int_{0}^{1}e^{\xi h_{0,m}(P_{i}(t))}e^{\xi h_{m,m+k}(P_{i}(t))}e^{\xi h_{m+k,n}(P_{i}(t))}|P_{i}^{\prime}(t)|dt\,.

By (3.29) and the event 1\mathcal{E}_{1} defined in (3.13), we obtain that there exists a constant C>0C>0 (not depending on n,mn,m) such that for all 1iK1\leq i\leq K and 0t10\leq t\leq 1,

|h0,m(Pi(t))h0,m(xki)|Cn2/3andhm,m+k(Pi(t))Ck(m+k).|h_{0,m}(P_{i}(t))-h_{0,m}(x_{k_{i}})|\leq Cn^{2/3}\quad\mbox{and}\quad h_{m,m+k}(P_{i}(t))\leq C\sqrt{k(m+k)}\,.

Combining this with (3.30) yields that

len(P~;Dn)\displaystyle{\rm len}(\widetilde{P};D_{n}) i=1K1eCn2/3+Ck(m+k)eξh0,m(xki)01eξhm+k,n(Pi(t))|Pi(t)|dt.\displaystyle\leq\sum_{i=1}^{K-1}e^{Cn^{2/3}+C\sqrt{k(m+k)}}e^{\xi h_{0,m}(x_{k_{i}})}\int_{0}^{1}e^{\xi h_{m+k,n}(P_{i}(t))}|P_{i}^{\prime}(t)|dt\,.

Combining this with (3.28) and (3.10), we further have

len(P~;Dn)eCn2/3anmk(p)i=1K12mkeξh0,m(xki).{\rm len}(\widetilde{P};D_{n})\leq e^{Cn^{2/3}}a_{n-m-k}^{(p)}\sum_{i=1}^{K-1}2^{-m-k}e^{\xi h_{0,m}(x_{k_{i}})}.

For each x2.3