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

Brownian loops on non-smooth surfaces
and the Polyakov-Alvarez formula

Minjae Park University of Chicago Joshua Pfeffer Columbia University Scott Sheffield MIT
Abstract

Let ρ\rho be compactly supported on D2D\subset\mathbb{R}^{2}. Endow 2\mathbb{R}^{2} with the metric eρ(dx12+dx22)e^{\rho}(dx_{1}^{2}+dx_{2}^{2}). As δ0\delta\to 0 the set of Brownian loops centered in DD with length at least δ\delta has measure

Vol(D)2πδ+148π(ρ,ρ)+o(1).\frac{\mathrm{Vol}(D)}{2\pi\delta}+\frac{1}{48\pi}(\rho,\rho)_{\nabla}+o(1).

When ρ\rho is smooth, this follows from the classical Polyakov-Alvarez formula. We show that the above also holds if ρ\rho is not smooth, e.g. if ρ\rho is only Lipschitz. This fact can alternatively be expressed in terms of heat kernel traces, eigenvalue asymptotics, or zeta regularized determinants. Variants of this statement apply to more general non-smooth manifolds on which one considers all loops (not only those centered in a domain DD).

We also show that the o(1)o(1) error is uniform for any family of ρ\rho satisfying certain conditions. This implies that if we weight a measure ν\nu on this family by the (δ\delta-truncated) Brownian loop soup partition function, and take the vague δ0\delta\to 0 limit, we obtain a measure whose Radon-Nikodym derivative with respect to ν\nu is exp(148π(ρ,ρ))\exp\bigl{(}\frac{1}{48\pi}(\rho,\rho)_{\nabla}\bigr{)}. When the measure is a certain regularized Liouville quantum gravity measure, a companion work [APPS20] shows that this weighting has the effect of changing the so-called central charge of the surface.


Acknowledgments. We thank Morris Ang, Ewain Gwynne, Camillo De Lellis, Sung-jin Oh, and Peter Sarnak for helpful comments. The authors were partially supported by NSF grants DMS 1712862 and DMS 2153742. J.P. was partially supported by a NSF Postdoctoral Research Fellowship under grant 2002159.

1 Introduction

Let us first recall a few standard definitions and observations. On a compact surface with boundary, the heat kernel trace can be written Z=Z(t)=spetΔ=etλnZ=Z(t)=\textrm{sp}\,e^{t\Delta}=\sum e^{t\lambda_{n}} where λn\lambda_{n} are the eigenvalues of the Laplace-Beltrami operator Δ\Delta. If I(s)I(s) is the number of λn-\lambda_{n} less than ss then

0etsI(s)𝑑s=0etsn=01s>λnds=n=0λnets𝑑s=1tetλn=Z(t)/t.\int_{0}^{\infty}e^{-ts}I(s)ds=\int_{0}^{\infty}e^{-ts}\sum_{n=0}^{\infty}1_{s>-\lambda_{n}}ds=\sum_{n=0}^{\infty}\int_{-\lambda_{n}}^{\infty}e^{-ts}ds=\frac{1}{t}\sum e^{t\lambda_{n}}=Z(t)/t.

In other words, Z/tZ/t is the Laplace transform of II. The asymptotics of ZZ (as t0t\to 0) are therefore closely related to the asymptotics of λn\lambda_{n} (as nn\to\infty). Weyl addressed the latter for bounded planar domains DD in 1911 [Wey11] (see discussion in [MS+67]) by showing λn2πnVol(D)-\lambda_{n}\sim\frac{2\pi n}{\mathrm{Vol}(D)} as nn\to\infty which is equivalent to

ZVol(D)4πtZ\sim\frac{\mathrm{Vol}(D)}{4\pi t} (1.1)

as t0t\to 0. In 1966 Kac gave higher order correction terms for ZZ on domains with piecewise linear boundaries (accounting for boundary length and corners) in his famously titled “Can you hear the shape of a drum?” which asks what features of the geometry of DD can be deduced from II (or equivalently from ZZ) [Kac66]. McKean and Singer extended these asymptotics from planar domains to smooth manifolds with non-zero curvature [MS+67] where the constant order correction term is a certain curvature integral. For two dimensional surfaces, with metric given by eρe^{\rho} times a flat metric, the integral δZ(t)t𝑑t\int_{\delta}^{\infty}\frac{Z(t)}{t}dt turns out to be a natural quantity whose small δ\delta asymptotics involve a constant order term that corresponds to the Dirichlet energy of ρ\rho (the so-called Polyakov-Alvarez formula, also known as the Polyakov-Ray-Singer or Weyl anomaly formula) [RS73, Pol81, Alv83, Sar87, OPS88]. This constant order Dirichlet energy term (which can also be formulated in terms of Brownian loop soups, see below) is the main concern of this paper.

Much of the literature assumes that ρ\rho is smooth and makes regular use of objects like curvature that are not well defined if ρ\rho is not C2C^{2}. But it is known [Hör68] that if ρ\rho is only C2C^{2} then Weyl’s law still holds, i.e. λn2πnVolρ(D)-\lambda_{n}\sim\frac{2\pi n}{\mathrm{Vol}_{\rho}(D)} [AHT18, Example 4.9] and Weyl’s law has been established in certain less smooth settings as well.111We remark that there is a general theory of metric measure spaces with the so-called “Riemannian curvature-dimension” (RCD) condition, not necessarily confined to conformal changes of flat metrics. They include Ricci limit spaces [Stu06, LV09], weighted Riemannian manifolds [Gri06], Alexandrov spaces [Pet10], and many others. A lower bound on Ricci curvature is a key ingredient of many useful estimates in geometric analysis, so Sturm, Lott and Villani [Stu06, LV09] initiated the study of a class of metric measure spaces with a generalized lower-Ricci-bound condition. This has been an active research topic for the last decade; see [Gig18] for an overview. The classical Weyl’s law and the short time asymptotics for heat kernels on these non-smooth metric measure spaces still hold [ZZ17, AHT18]. Many aspects of the theory are stable under the pointed measured Gromov-Hausdorff topology; for example eigenvalues, heat kernels, and Green’s function converge uniformly [Din02, ZZ17], Brownian motions converge weakly [Suz19], etc. Therefore, Weyl’s law holds for any RCD space with a measure that can be reasonably approximated. On the other hand, the short time expansion used to define the functional determinant does not exist in this non-smooth setting, so it is not clear if the zeta regularization procedure is also stable. The problem is somewhat different when the regularity is below C2C^{2}, since curvature is no longer well-defined everywhere and the relevant estimates no longer hold pointwise and instead hold in an average sense. The primary purpose of this note is to extend some of the basic results in this subject about the conformal anomaly (the Dirichlet energy of ρ\rho) to ρ\rho that are less regular—e.g., only Lipschitz—and to show that the rates of convergence can be made to hold uniformly across certain families of ρ\rho values.

This paper is motivated in part by another work by the authors [APPS20] in which similar results are formulated in terms of the so-called Brownian loop measures which were introduced in [LW04] and are related to heat kernel traces on planar domains in e.g. [Dub09, Wan18] as well as [APPS20]. The results here are useful in the context of [APPS20] because they strengthen the sense in which one can say that “decorating” regularized Liouville quantum gravity surfaces by Brownian loop soups has the effect of changing their central charge. We will formulate the results in this paper solely in terms of Brownian loop measures and their generalizations. (The relationship to heat kernels is explained in [APPS20].)

In addition to the weaker regularity assumptions and the use of generalized loop measures, there are several smaller differences between our presentation and the classical approach in e.g. [MS+67]: we work in the conformal gauge throughout and do all our calculations in terms of ρ\rho, we index loops by their Euclidean center rather than by a typical point on the loop (which would be more similar to the heat kernel approach), and we establish the Polyakov-Alvarez formula in terms of Dirichlet energy directly rather than first establishing an equivalent curvature integral.

Although we encounter some complexity due to the non-smoothness of ρ\rho, we also take advantage of the extra simplicity of the two-dimensional setting, where the manifold is completely determined by a conformal factor.

Finally, we note that there is a great deal of additional work in this area, and we cannot begin to survey it all. For example, reference texts such as [BGV03, Gil18] explore heat kernel traces in greater generality: dimensions other than 22, operators other than the Laplacian, etc. Other works extend the behavior known for compact smooth manifolds to specific non-smooth manifolds such as those with conical singularities or boundary corners (which both correspond to logarithmic singularities in ρ\rho) [Moo99, Kok13, She13, She15, AR18, Gre21, Kal21] or to non-compact surfaces [AAR13]. There are also many open problems in this subject, which spans probability, geometry, number theory, mathematical physics, and analysis. We present a few of these questions in Section 6. We hope that the techniques and perspectives presented here will facilitate progress on these problems and perhaps also find applications in other contexts where Weyl’s formula and the Polyakov-Alvarez term appear.

1.1 Main result

Let \mathcal{L} denote the set of zero-centered unit-length loops in 2\mathbb{R}^{2}. We define the Brownian loop measure in the plane by encoding each loop in the plane as an element of 2×(0,)×\mathbb{R}^{2}\times(0,\infty)\times\mathcal{L} and formulating the Brownian loop measure as a measure on this product space.

Definition 1.1.

We express every loop LL in 2\mathbb{R}^{2} by the triple (x,t,)(x,t,\ell), where

  • t=len(L)t=\mathrm{len}(L) is the length of LL, where we define the length of a path as the length of its parametrizing interval.

  • x=cen(L)x=\mathrm{cen}(L) is the center of LL, Euclidean center of mass of LL, which is equal to t10tL(s)𝑑st^{-1}\int_{0}^{t}L(s)ds.

  • \ell is the zero-centered unit-length loop st1/2(L(ts)x)s\mapsto t^{-1/2}(L(ts)-x) obtained from LL by translating the center to zero and rescaling time by t1t^{-1} and space by t1/2t^{-1/2}.

Definition 1.2.

We define the Brownian loop measure on 2\mathbb{R}^{2} as the measure on loops (x,t,)(x,t,\ell) in 2\mathbb{R}^{2} given by

12πt2dxdtd,\frac{1}{2\pi t^{2}}dx\,dt\,d\ell,

where dxdx denotes Lebesgue measure on 2\mathbb{R}^{2}, dtdt is Lebesgue measure on (0,)(0,\infty) and dd\ell is the probabilistic law of the random loop in \mathcal{L} obtained by first sampling a two-dimensional Brownian bridge on [0,1][0,1] and then subtracting its mean.222Equivalently dd\ell on the complex plane is the law of the complex-valued GFF indexed by the unit-length circle—with additive constant chosen to make the mean zero. In particular, dd\ell is invariant under rotations of that circle.

The mass of the set of Brownian loops centered in DD with size greater than δ\delta is given by

Dδ12πt2𝑑t𝑑x=Vol(D)2πδ.\int_{D}\int_{\delta}^{\infty}\frac{1}{2\pi t^{2}}\,dt\,dx=\frac{\mathrm{Vol}(D)}{2\pi\delta}. (1.2)

In particular, (1.2) implies that no matter how small δ\delta is, most loops with length t(δ,)t\in(\delta,\infty) have length of order δ\delta: half of them have t<2δt<2\delta, ninety-five percent have t<20δt<20\delta, and so forth. Also, the fact that (1.2) tends to 0 as δ\delta\to\infty informally means that there are very few large loops centered in DD.

Our main result describes how this mass of Brownian loops changes when we measure the length of loops with respect to a different metric eρ|dz|2e^{\rho}|dz|^{2} on the plane, for ρ\rho a Lipschitz function supported in DD. We begin by defining the length of a Brownian loop in the metric eρ|dz|2e^{\rho}|dz|^{2}, which we call its ρ\rho-length. If the loop were a smooth curve, we would compute its ρ\rho-length by integrating eρ/2e^{\rho/2} along the curve. Since Brownian loops have Hausdorff dimension 22, we instead define its ρ\rho-length by integrating eρe^{\rho} along the loop, so that it has the same scale factor as area.

Definition 1.3.

Let (M,g)(M,g) be a smooth two-dimensional Riemannian manifold, and let ρ\rho be a function on MM. We define the ρ\rho-length lenρ(L)\mathrm{len}_{\rho}(L) of a loop LL as 0len(L)eρ(L(s))𝑑s\int_{0}^{\mathrm{len}(L)}e^{\rho(L(s))}ds. We define the ρ\rho-volume form Volρ\mathrm{Vol}_{\rho} as the volume form associated to (M,eρg)(M,e^{\rho}g), and we write Vol:=Vol0\mathrm{Vol}:=\mathrm{Vol}_{0}.

Except in Theorem 1.9 and Section 4, we always take (M,g)(M,g) in Definition 1.3 to be the Euclidean plane. We first describe the space of functions in the scope of this section.

Theorem 1.4.

Let DD be a bounded open subset of 2\mathbb{R}^{2}, and Lip(D)\operatorname{Lip}(D) be the space of real-valued Lipschitz functions that vanish outside of DD. Suppose that Lip(D)\mathcal{B}\subset\operatorname{Lip}(D) is a collection of functions that (1) has uniformly bounded Lipschitz constants, and (2) is precompact in W1,1(D)W^{1,1}(D).

Then as δ0\delta\to 0, the μloop\mu^{\text{loop}}-mass of loops centered in DD with ρ\rho-length at least δ\delta, with respect to the Brownian loop measure, is given by

Volρ(D)2πδ+148π(ρ,ρ)+o(1)\frac{\mathrm{Vol}_{\rho}(D)}{2\pi\delta}+\frac{1}{48\pi}(\rho,\rho)_{\nabla}+o(1) (1.3)

with the convergence uniform over ρ\rho\in\mathcal{B}.

Remark 1.5 (Uniform boundedness).

The conditions (1) and (2) imply that the functions in \mathcal{B} are uniformly bounded.

Remark 1.6 (General pp).

Since we require for now that the Lipschitz constant is uniformly bounded and the domain is bounded, precompactness in W1,1W^{1,1} is equivalent to precompactness in W1,pW^{1,p} for any fixed p(1,)p\in(1,\infty). In particular, Theorem 1.4 could have been formulated using precompactness in W1,2W^{1,2} instead of W1,1W^{1,1}. Let us also remark that the space W1,2(D)W^{1,2}(D) is equivalent to the space H1(D)H^{1}(D) of ρ\rho for which (ρ,ρ)(\rho,\rho)_{\nabla} is finite.

Remark 1.7 (Precompactness and uniform equicontinuity).

Recall the Fréchet-Kolmogorov theorem (e.g., see [BB11]): let DnD\subset\mathbb{R}^{n} be a bounded domain, and 1p<1\leq p<\infty. A subset 𝒜Lp(D)\mathcal{A}\subset L^{p}(D) is precompact if and only if 𝒜\mathcal{A} is bounded in Lp(D)L^{p}(D) and

supu𝒜D|u(x+h)u(x)|p𝑑x0as h0,\sup_{u\in\mathcal{A}}\,\int_{D}\left|u(x+h)-u(x)\right|^{p}\,dx\to 0\quad\text{as }h\to 0, (1.4)

where uu is extended to the function on n\mathbb{R}^{n} whose value outside DD is zero.

We will use this equivalent characterization of precompactness in some of our proofs, usually referred as the uniform equicontinuity condition in LpL^{p}. In particular, we will apply this, in the case p=1p=1, to the set 𝒜\mathcal{A} of the gradients of the functions in the set \mathcal{B} from the statement of Theorem 1.4.

Remark 1.8 (The uniform equicontinuity condition is necessary).

As mentioned above, the precompactness hypothesis is equivalent to a type of uniform equicontinuity hypothesis. This hypothesis—or some similar condition on the functions ρ\rho\in\mathcal{B}—is necessary for the conclusion of Theorem 1.4 (or Theorem 1.12) to hold. Simply requiring all surfaces in \mathcal{B} to be C1C^{1} with a universal bound on |ρ||\nabla\rho| would not suffice. For example, in the Theorem 1.4 setting, \mathcal{B} could contain a sequence ρ1,ρ2,\rho_{1},\rho_{2},\ldots of C1C^{1} functions that converge uniformly to zero with (ρj,ρj)=1(\rho_{j},\rho_{j})_{\nabla}=1 for each jj. We can construct such a sequence of functions ρj\rho_{j} by arranging for ρj\nabla\rho_{j} to oscillate between fixed opposite values, with the oscillation rate becoming faster as jj\to\infty. (A simple example of such a family of functions on the torus [0,2π)2[0,2\pi)^{2} is given by a constant multiple of ρj((a,b)))=j1sin(ja)\rho_{j}\Bigl{(}(a,b)\Bigr{)})=j^{-1}\sin(ja); we can define ρj\rho_{j} similarly on the planar domain DD by tapering the sine functions to zero near the boundary of DD.) We can also perturb the functions to arrange that Volρj(D)=Volρ(D)\mathrm{Vol}_{\rho_{j}}(D)=\mathrm{Vol}_{\rho}(D) for all jj. This set of functions \mathcal{B} does not satisfy (1.3): for any fixed δ\delta, one can easily show that

limjμ{L:cen(L)D,lenρj(L)δ}Vol(D)δ=0,\lim_{j\to\infty}\mu\{L:\mathrm{cen}(L)\in D,\,\mathrm{len}_{\rho_{j}}(L)\geq\delta\}-\frac{\mathrm{Vol}(D)}{\delta}=0, (1.5)

even though for each fixed jj we have

limδ0(μ{L:cen(L)D,lenρj(L)δ}Vol(D)δ)=b/2.\lim_{\delta\to 0}\bigl{(}\mu\{L:\mathrm{cen}(L)\in D,\,\mathrm{len}_{\rho_{j}}(L)\geq\delta\}-\frac{\mathrm{Vol}(D)}{\delta}\bigr{)}=b/2. (1.6)

If the limit in (1.6) were uniform in jj, we could choose a δ\delta with μ{L:cen(L)D,lenρj(L)δ}>b/4\mu\{L:\mathrm{cen}(L)\in D,\,\mathrm{len}_{\rho_{j}}(L)\geq\delta\}>b/4 for all jj, and (1.5) would not hold for that δ\delta.333One might wonder whether it is enough have ρ\rho in the Hilbert space defined by the inner product (ρ,ρ)(\rho,\rho)_{\nabla}, i.e., the Sobolev space W1,2(D)=H1(D)W^{1,2}(D)=H^{1}(D), with (say) zero boundary conditions. Such a ρ\rho can be nowhere differentiable [Ser61], so one would also need to modify the condition on \mathcal{B}. See Question 6.1.

We extend this result to general surfaces. The statement of the theorem involves the notion of the zeta-regularized determinant of the Laplacian, as defined, e.g., in [Alv83, Sar87]. (However, it is not necessary to understand the definition of the zeta-regularized determinant to follow the proof of Theorem 1.9.)

Theorem 1.9.

Let (M,g)(M,g) be a fixed compact smooth two-dimensional Riemannian manifold, and we let μloop\mu^{\text{loop}} denote the Brownian loop measure on (M,g)(M,g). Let KK be the Gaussian curvature on MM, let Δ\Delta be the Laplacian associated to (M,g)(M,g), and let detζΔ\det_{\zeta}^{\prime}\Delta denote its zeta-regularized determinant. Let \mathcal{B} be a family of Lipschitz functions that (1) has uniformly bounded Lipshitz constants, and (2) is precompact in W1,1(M)W^{1,1}(M).

Then the μloop\mu^{\text{loop}}-mass of loops with ρ\rho-length between δ\delta and CC is given by

Volρ(M)2πδχ(M)6(logδ2+γ)+logC+γ+148πM(ρ2+2Kρ)Vol(dz)\displaystyle\frac{\mathrm{Vol}_{\rho}(M)}{2\pi\delta}-{\frac{\chi(M)}{6}}(\log\frac{\delta}{2}+\upgamma)+\log{C}+\upgamma+\frac{1}{48\pi}\int_{M}(\|\nabla\rho\|^{2}+2K\rho)\mathrm{Vol}(dz)
+logVol(M)logVolρ(M)logdetζΔ+oδ(1)+oC(C(1+ε)/2),\displaystyle\qquad+\log\mathrm{Vol}(M)-\log\mathrm{Vol}_{\rho}(M)-\log\det\nolimits_{\zeta}^{\prime}\Delta+o_{\delta}(1)+o_{C}(C^{(-1+\varepsilon)/2}), (1.7)

with the convergence as δ0\delta\to 0 and CC\to\infty uniform over ρ\rho\in\mathcal{B}, where γ0.5772\upgamma\approx 0.5772 is the Euler-Mascheroni constant.

For simplicity, we have addressed just the compact manifold case, but one could prove a similar result for manifolds with boundary, see Question 6.3 where we give a heuristic justification in the preceding paragraph; the resulting expression would include a boundary term which is of order δ1/2\delta^{-1/2}.

Observe that, for smooth ρ\rho, the expression in the second line (1.7) of the above expression is equal to logdetζΔρ-\log\det_{\zeta}^{\prime}\Delta_{\rho}, where Δρ\Delta_{\rho} is the Laplacian associated to (M,eρg)(M,e^{\rho}g). The expression (1.7) for logdetζΔρ-\log\det_{\zeta}^{\prime}\Delta_{\rho} is known as the Polyakov-Alvarez formula; see, e.g., [APPS20, Proposition 6.9]. Thus, for smooth ρ\rho, Theorem 1.9 reduces to a relation between the Brownian loop measure and the zeta-regularized Laplacian determinant, which was shown in [APPS20, Theorem 1.3].

In fact, we prove a slight generalization of Theorem 1.4 in which we consider a more general class of loop measures.

Definition 1.10.

The expected occupation measure of a random variable Z:[0,T]2Z:[0,T]\rightarrow\mathbb{R}^{2} is the function θ:2(0,)\theta:\mathbb{R}^{2}\rightarrow(0,\infty) such that, for each measurable set A2A\subset\mathbb{R}^{2}, the set {t[0,T]:Z(t)T}\{t\in[0,T]:Z(t)\in T\} has expected Lebesgue measure Aθ(x)𝑑x\int_{A}\theta(x)dx.

Definition 1.11.

We define a generalized loop measure μ\mu as a measure on loops (x,t,)(x,t,\ell) in 2\mathbb{R}^{2} given by

1t2dxdtd,\frac{1}{t^{2}}dx\,dt\,d\ell,

where dxdx denotes Lebesgue measure on 2\mathbb{R}^{2}, dtdt is Lebesgue measure on (0,)(0,\infty) and dd\ell is an arbitrary rotationally invariant measure on loops in \mathcal{L} whose expected occupation measure is a Schwartz distribution (but not necessarily Gaussian as for the Brownian loop measure). We denote by bb the second central moment of the first—or equivalently, second—coordinate of a random variable whose density is this expected occupation measure.

Definition 1.11 is the same as Definition 1.2, except that we no longer insist that dd\ell be the Brownian bridge (and we have removed the 2π2\pi factor as it is less natural for general dd\ell). The measure dd\ell can be supported on the space of circular loops, square-shaped loops, or outer boundaries of Brownian loops, etc. The μ\mu from Definition 1.11 need not have the same conformal symmetries as the Brownian loop measure. Even if dd\ell is supported on smooth loops (rather than Brownian loops) we parameterize the space of loops as in Definition 1.1, so that (0,t,)(0,t,\ell) represents the loop that traces t\sqrt{t}\ell over time duration tt. In the special case of the Brownian loop measure, b=1/12b=1/12. (See Proposition 3.9.) The Schwartz distribution assumption in Definition 1.11 does not seem necessary for Theorem 1.12 to hold, but we have included it to simplify the calculations in the proof of Proposition 3.1 below.

Theorem 1.12.

Let DD and \mathcal{B} be as in Theorem 1.4, and let μ\mu be a generalized loop measure in the sense of Definition 1.11. Then as δ0\delta\rightarrow 0, the μ\mu-mass of loops centered in DD with ρ\rho-length at least δ\delta is given by

Volρ(D)δ+b2(ρ,ρ)+o(1),\frac{\mathrm{Vol}_{\rho}(D)}{\delta}+\frac{b}{2}(\rho,\rho)_{\nabla}+o(1), (1.8)

with the convergence uniform over ρ\rho\in\mathcal{B}.

1.2 Proof outline

In this section, we let \mathcal{B} be a fixed collection of functions ρ\rho satisfying the hypotheses of Theorem 1.12. To prove Theorem 1.12, we compare lenρ\mathrm{len}_{\rho} to a simpler notion of the length of a loop with respect to eρ|dz|2e^{\rho}|dz|^{2}, in which we approximate eρe^{\rho} along the loop by its value at the center of the loop.

Definition 1.13.

We define the center ρ\rho-length clenρ(L)\mathrm{clen}_{\rho}(L) of a loop LL in 2\mathbb{R}^{2} centered at a point xx as 0teρ(x)𝑑s=eρ(x)t\int_{0}^{t}e^{\rho(x)}ds=e^{\rho(x)}t.

We observe that the cutoff clenρ(L)=δ\mathrm{clen}_{\rho}(L)=\delta corresponds to a unique value of len(L)\mathrm{len}(L):

Proposition 1.14.

Let δ>0\delta>0, xDx\in D and \ell\in\mathcal{L}, and set β:=eρ(x)δ\beta:=e^{-\rho(x)}\delta. The loop L=(x,t,)L=(x,t,\ell) satisfies clenρ(L)=δ\mathrm{clen}_{\rho}(L)=\delta iff t=βt=\beta, and clenρ(L)δ\mathrm{clen}_{\rho}(L)\geq\delta iff tβt\geq\beta.

Proof.

The result follows trivially from the definition of center ρ\rho-length. ∎

Proposition 1.14 immediately implies the following trivial analogue of Theorem 1.12 for center ρ\rho-length.

Proposition 1.15.

The mass of loops LL centered in DD with clenρδ\mathrm{clen}_{\rho}\geq\delta with respect to the Brownian loop measure is given by

Dβ1t2𝑑t𝑑𝑑x=Dβ1𝑑x=Deρ(x)δ𝑑x=Volρ(D)δ.\int_{D}\int_{\mathcal{L}}\int_{\beta}^{\infty}\frac{1}{t^{2}}\,dt\,d\ell\,dx=\int_{D}\beta^{-1}dx=\int_{D}\frac{e^{\rho(x)}}{\delta}dx=\frac{\mathrm{Vol}_{\rho}(D)}{\delta}. (1.9)
Proof.

The result is an immediate consequence of Proposition 1.14. ∎

We can therefore restate Theorem 1.12 as the assertion that if we change our notion of loop length from clenρ\mathrm{clen}_{\rho} to lenρ\mathrm{len}_{\rho}, the μ\mu-mass of loops with length δ\geq\delta increases by b2(ρ,ρ)\frac{b}{2}(\rho,\rho)_{\nabla}, up to an error that is o(1)o(1) as δ0\delta\to 0 uniformly in ρ\rho\in\mathcal{B}.

We divide the proof of Theorem 1.12 into two stages. First, in Section 2 we show that, up to a uniform o(1)o(1) error, replacing clenρ\mathrm{clen}_{\rho} with lenρ\mathrm{len}_{\rho} has the effect of subtracting δ1\delta^{-1} times the average discrepancy between the value of eρe^{\rho} along a Brownian loop and the value of eρe^{\rho} at its center.

Lemma 1.16.

Consider a loop sampled from μ\mu conditioned to have its center in DD and length β\beta. Let XX denote its center, and let ZZ denote a point on the loop sampled uniformly with respect to length. Then the μ\mu-mass of loops LL with center in DD and lenρ(L)δ\mathrm{len}_{\rho}(L)\geq\delta is equal to the μ\mu-mass of loops LL with center in DD and clenρ(L)δ\mathrm{clen}_{\rho}(L)\geq\delta, minus

1δ𝔼[eρ(Z)eρ(X)]+o(1),\frac{1}{\delta}\mathbb{E}[e^{\rho(Z)}-e^{\rho(X)}]+o(1), (1.10)

with the o(1)o(1) error tending to 0 as δ0\delta\to 0 at a rate that is uniform in ρ\rho\in\mathcal{B}. (In (1.10) the expectation is w.r.t. the overall law of XX and ZZ as described above.)

We then complete the proof of Theorem 1.12 in Section 3 by showing that the quantity (1.10) equals b2(ρ,ρ)\frac{b}{2}(\rho,\rho)_{\nabla} up to a uniform o(1)o(1) error.

Lemma 1.17.

The quantity (1.10) equals b2(ρ,ρ)\frac{b}{2}(\rho,\rho)_{\nabla} plus an error term that converges to 0 as δ0\delta\rightarrow 0 uniformly in ρ\rho\in\mathcal{B}.

2 Loop mass difference vs. expected length discrepancy

In this section, we prove Lemma 1.16 in three steps.

Step 1: Establishing a length threshold α>0\alpha>0 corresponding to ρ\rho-length δ\delta. We saw in Proposition 1.14 that, with β=eρ(x)δ\beta=e^{-\rho(x)}\delta, we have clenρ(L)=δ\mathrm{clen}_{\rho}(L)=\delta if and only if t=βt=\beta, and clen(L)δ\mathrm{clen}(L)\geq\delta if and only if tβt\geq\beta. To prove Lemma 1.16, we establish a similar result for ρ\rho-length. We will show in Proposition 2.4 below that, for xDx\in D and \ell\in\mathcal{L} with the diameter of (x,δ,)(x,\delta,\ell) sufficiently small, there exists a threshold α>0\alpha>0 such that lenρ((x,t,))=δ\mathrm{len}_{\rho}((x,t,\ell))=\delta if and only if t=αt=\alpha, and len((x,t,))δ\mathrm{len}((x,t,\ell))\geq\delta if and only if tαt\geq\alpha. We note that, unlike β\beta, the threshold α\alpha may depend on \ell as well as xx and δ\delta.

Step 2: Relating the difference in the masses of loops to the quantity α1\alpha^{-1}. We saw in Proposition 1.15 that the μ\mu-mass of loops with center ρ\rho-length δ\geq\delta can be expressed as an integral of β1\beta^{-1}. In Proposition 2.5 below, we similarly express the μ\mu-mass of loops with ρ\rho-length δ\geq\delta as an integral of α1\alpha^{-1}, plus a uniform o(1)o(1) error. This reduces the task of proving Lemma 1.16 to analyzing the difference of integrands α1β1\alpha^{-1}-\beta^{-1}.

Step 3: Expressing the difference α1β1\alpha^{-1}-\beta^{-1} in terms of a difference in lengths. We first express the difference αβ\alpha-\beta in terms of a difference between the ρ\rho-length and center ρ\rho-length of a loop (Proposition 2.7). We then apply this result in Proposition 2.8 to derive a similar expression for α1β1\alpha^{-1}-\beta^{-1}, which immediately yields Lemma 1.16.

Having described the main steps of the proof of Lemma 1.16, we now proceed with Step 1—proving the existence of the threshold α\alpha. As we observed in Proposition 1.14, the existence of the analogous threshold β\beta for center ρ\rho-length is trivial, since for fixed ρ,x,\rho,x, and \ell, the function tclenρ(L)t\mapsto\mathrm{clen}_{\rho}(L) is linear with slope eρ(x)e^{\rho(x)}. This is not the case for ρ\rho-length, so we proceed by showing its derivative as a function of tt is positive on a sufficiently large interval. We first observe that, since the functions ρ\rho\in\mathcal{B} are uniformly bounded above and below, we can crudely bound the function tlenρ(L)t\mapsto\mathrm{len}_{\rho}(L) between two linear functions uniformly in ρ,x,\rho,x, and \ell.

Proposition 2.1.

There exists a constant Λ>0\Lambda>0 such that, for each ρ\rho\in\mathcal{B}, xDx\in D and \ell\in\mathcal{L}, the loop L=(x,t,)L=(x,t,\ell) satisfies

Λ1lenρ(L)len(L)Λ.\Lambda^{-1}\leq\frac{\mathrm{len}_{\rho}(L)}{\mathrm{len}(L)}\leq\Lambda. (2.1)

In other words, len(L)\mathrm{len}(L) and lenρ(L)\mathrm{len}_{\rho}(L) length agree up to a universal constant factor.

Proof.

The lemma follows immediately from the fact that ρ\rho is bounded from above and below by a constant uniform in ρ\rho\in\mathcal{B}, as noted in Remark 1.5. ∎

In addition, the collection of functions {eρ}ρ\{e^{\rho}\}_{\rho\in\mathcal{B}} satisfies the same conditions of \mathcal{B}, possibly with different bounds.

Proposition 2.2.

The functions eρe^{\rho} for ρ\rho\in\mathcal{B} also have uniformly bounded Lipschitz constants and are precompact in W1,1(D)W^{1,1}(D).

Proof.

Since the functions ρ\rho\in\mathcal{B} are uniformly bounded, their images are contained in some finite closed interval. The exponential function is Lipschitz on any finite closed interval, so the composition expρ\exp\circ\rho is also Lipschitz. Other conditions are straightforward to check. ∎

We now apply the uniform boundedness of the Lipschitz constants of eρe^{\rho} for ρ\rho\in\mathcal{B} to show that, when diam(L)=diam()t\mathrm{diam}(L)=\mathrm{diam}(\ell)\sqrt{t} is not too large, the derivative of the function tlenρ(L)t\mapsto\mathrm{len}_{\rho}(L) is uniformly close to that of the linear function tclenρ(L)t\mapsto\mathrm{clen}_{\rho}(L).

Proposition 2.3.

For any ε>0\varepsilon>0, there exists d=d(ε)>0d=d(\varepsilon)>0 and a family of sets {Dρ,ε}ρ\{D^{\rho,\varepsilon}\}_{\rho\in\mathcal{B}} with Vol(Dρ,ε)ε\mathrm{Vol}(D^{\rho,\varepsilon})\leq\varepsilon such that the following is true. For each ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L}, the derivative

tlenρ((x,t,))\frac{\partial}{\partial t}\mathrm{len}_{\rho}((x,t,\ell))

exists and differs from

tclenρ((x,t,))=eρ(x)\frac{\partial}{\partial t}\mathrm{clen}_{\rho}((x,t,\ell))=e^{\rho(x)}

by at most εdiam()t\varepsilon\,\mathrm{diam}(\ell)\sqrt{t} for almost every t>0t>0 with diam()t<d\mathrm{diam}(\ell)\sqrt{t}<d. Furthermore, there exists a constant Λ~>0\widetilde{\Lambda}>0 such that the previous statement it true for arbitrary dd and Dρ,ε=D^{\rho,\varepsilon}=\emptyset when we choose ε=Λ~\varepsilon=\widetilde{\Lambda}.

Proof.

By Rademacher’s theorem, any Lipschitz function is differentiable almost everywhere. In particular, by Proposition 2.2, there exists some constant Λ~>0\widetilde{\Lambda}>0 that does not depend on ρ\rho such that |(eρ)|<23Λ~\left|\nabla(e^{\rho})\right|<\frac{2}{3}\widetilde{\Lambda} for almost every xDx\in D for all ρ\rho\in\mathcal{B}, where (eρ)\nabla(e^{\rho}) is a weak gradient of eρe^{\rho}. In addition, as \mathcal{B} is precompact in W1,1(D)W^{1,1}(D), it follow from (1.4) that there exists some d>0d>0 such that the set

Dρ,ε:={xD:esssup|h|d|(eρ)(x+h)eρ(x)|>23ε}D^{\rho,\varepsilon}:=\{x\in D:\operatorname*{ess\,sup}_{|h|\leq d}\left|(\nabla e^{\rho})(x+h)-\nabla e^{\rho}(x)\right|>\frac{2}{3}\varepsilon\} (2.2)

satisfies Vol(Dρ,ε)ε\mathrm{Vol}(D^{\rho,\varepsilon})\leq\varepsilon for each ρ\rho\in\mathcal{B}.

For fixed ρ\rho\in\mathcal{B}, xDx\in D, and \ell\in\mathcal{L}, we can write lenρ((x,t,))\mathrm{len}_{\rho}((x,t,\ell)) as A(t)tA(\sqrt{t})t, where A(r)=Aρ,x,(r):=01eρ(r(s)+x)𝑑sA(r)=A_{\rho,x,\ell}(r):=\int_{0}^{1}e^{\rho}(r\ell(s)+x)ds. We express a weak tt-derivative of lenρ((x,t,))\mathrm{len}_{\rho}((x,t,\ell)) in terms of AA as

tlenρ((x,t,))=t(A(t)t)=12t1/2A(t)t+A(t)=12A(t)t+A(t).\frac{\partial}{\partial t}\mathrm{len}_{\rho}((x,t,\ell))=\frac{\partial}{\partial t}(A(\sqrt{t})t)=\frac{1}{2}t^{-1/2}A^{\prime}(\sqrt{t})t+A(\sqrt{t})=\frac{1}{2}A^{\prime}(\sqrt{t})\sqrt{t}+A(\sqrt{t}). (2.3)

Since

A(r)\displaystyle A^{\prime}(r) =01(s)(eρ)(r(s)+x)𝑑s\displaystyle=\int_{0}^{1}\ell(s)\cdot(\nabla e^{\rho})(r\ell(s)+x)ds (2.4)
=01(s)((eρ)(r(s)+x)(eρ)(x))𝑑s,\displaystyle=\int_{0}^{1}\ell(s)\cdot((\nabla e^{\rho})(r\ell(s)+x)-(\nabla e^{\rho})(x))ds, (2.5)

we can bound |A(t)||A^{\prime}(\sqrt{t})| from above by

diam()suph2|(eρ)(x+h)|23Λ~diam()\mathrm{diam}(\ell)\sup_{h\in\mathbb{R}^{2}}\left|(\nabla e^{\rho})(x+h)\right|\leq\frac{2}{3}\widetilde{\Lambda}\,\mathrm{diam}(\ell)

for almost every t>0t>0, using (2.4). On the other hand, if xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, we use (2.2) and (2.5) to bound |A(t)||A^{\prime}(\sqrt{t})| from above by

diam()esssup|h|diam()t|(eρ)(x+h)eρ(x)|23εdiam()\mathrm{diam}(\ell)\operatorname*{ess\,sup}_{|h|\leq\mathrm{diam}(\ell)\sqrt{t}}\left|(\nabla e^{\rho})(x+h)-\nabla e^{\rho}(x)\right|\leq\frac{2}{3}\varepsilon\,\mathrm{diam}(\ell)

for almost every t>0t>0 with diam()t<d\mathrm{diam}(\ell)\sqrt{t}<d.

Note that (2.3) gives

|tlenρ((x,t,))A(t)|=12|A(t)|t,\left|\frac{\partial}{\partial t}\mathrm{len}_{\rho}((x,t,\ell))-A(\sqrt{t})\right|=\frac{1}{2}\left|A^{\prime}(\sqrt{t})\right|\sqrt{t},

and we also have

|A(t)eρ(x)|=|A(t)A(0)|0t|A(s)|𝑑ssup0st|A(s)|t.|A(\sqrt{t})-e^{\rho(x)}|=|A(\sqrt{t})-A(0)|\leq\int_{0}^{\sqrt{t}}|A^{\prime}(s)|ds\leq\sup_{0\leq s\leq\sqrt{t}}\left|A^{\prime}(s)\right|\sqrt{t}. (2.6)

Combining these two inequalities with the estimates for |A(t)|\left|A^{\prime}(\sqrt{t})\right| in two different cases gives the desired result. ∎

We use Proposition 2.3 to show that, just as for center ρ\rho-length, there is a unique value α\alpha of len(L)\mathrm{len}(L) corresponding to lenρ(L)=δ\mathrm{len}_{\rho}(L)=\delta.

Proposition 2.4.

We can choose d>0d_{*}>0 such that the following holds. For each ρ\rho\in\mathcal{B} and xDx\in D, if L=(x,δ,)L=(x,\delta,\ell) is a loop with diam()δ\mathrm{diam}(\ell)\sqrt{\delta} less than dd_{*}, then there is a unique positive α=α(x,δ,,ρ)\alpha=\alpha(x,\delta,\ell,\rho) such that lenρ((x,t,))=δ\mathrm{len}_{\rho}((x,t,\ell))=\delta iff t=αt=\alpha and lenρ((x,t,))>δ\mathrm{len}_{\rho}((x,t,\ell))>\delta iff t>αt>\alpha. (If diam()δ\mathrm{diam}(\ell)\sqrt{\delta} is d\geq d_{*}, then we arbitrarily set α=δ\alpha=\delta, so that α\alpha is defined for every loop L=(x,δ,)L=(x,\delta,\ell).)

Proof.

Let Λ\Lambda be the constant in Proposition 2.1. If tΛδt\geq\Lambda\delta, then lenρ((x,t,))Λ1tδ\mathrm{len}_{\rho}((x,t,\ell))\geq\Lambda^{-1}t\geq\delta by Proposition 2.1. Thus, it suffices to show that, for diam()δ\mathrm{diam}(\ell)\sqrt{\delta} sufficiently small, the function lenρ((x,t,))\mathrm{len}_{\rho}((x,t,\ell)) is strictly increasing in tt when t(0,Λδ)t\in(0,\Lambda\delta). We prove this fact by analyzing tlenρ((x,t,))\frac{\partial}{\partial t}\mathrm{len}_{\rho}((x,t,\ell)) using Proposition 2.3.

Let Λ~\widetilde{\Lambda} be the constant in Proposition 2.3 so that for all ρ\rho\in\mathcal{B}, xD,x\in D,\ell\in\mathcal{L}, and almost every t>0t>0,

|tlenρ((x,t,))eρ(x)|Λ~diam()t.\left|\frac{\partial}{\partial t}\mathrm{len}_{\rho}((x,t,\ell))-e^{\rho(x)}\right|\leq\widetilde{\Lambda}\,\mathrm{diam}(\ell)\sqrt{t}. (2.7)

We choose diam()δ\mathrm{diam}(\ell)\sqrt{\delta} sufficiently small less than d>0d_{*}>0 such that, for almost every tΛδt\leq\Lambda\delta, the bound Λ~diam()t\widetilde{\Lambda}\,\mathrm{diam}(\ell)\sqrt{t} in (2.7) is less than 12eρ(x)\frac{1}{2}e^{\rho(x)}. This means tlenρ((x,t,))\frac{\partial}{\partial t}\mathrm{len}_{\rho}((x,t,\ell)) is strictly positive for almost every tΛδt\leq\Lambda\delta. Since lenρ((x,t,))\mathrm{len}_{\rho}((x,t,\ell)) is a continous function in tt, we conclude that lenρ((x,t,))\mathrm{len}_{\rho}((x,t,\ell)) is strictly increasing on (0,Λδ)(0,\Lambda\delta). ∎

When xx and \ell are fixed, α\alpha gives the Euclidean tt value of len\mathrm{len} that corresponds to a δ\delta value of lenρ\mathrm{len}_{\rho}. We obtain the following analogue of Proposition 1.15 with lenρ\mathrm{len}_{\rho} in place of clenρ\mathrm{clen}_{\rho}.

Proposition 2.5.

The μ\mu-mass of the set of loops LL with center in DD and lenρ(L)δ\mathrm{len}_{\rho}(L)\geq\delta is equal to

Dα1t2𝑑t𝑑𝑑x=Dα1𝑑𝑑x\int_{D}\int_{\mathcal{L}}\int_{\alpha}^{\infty}\frac{1}{t^{2}}\,dt\,d\ell\,dx=\int_{D}\int_{\mathcal{L}}\alpha^{-1}\,d\ell\,dx (2.8)

plus an error term that tends to 0 as δ0\delta\rightarrow 0 at a rate that is uniform in ρ\rho\in\mathcal{B}.

The reason the mass of loops does not exactly equal (2.8) is that α\alpha is improperly defined when diam()δd\mathrm{diam}(\ell)\sqrt{\delta}\geq d_{*}, with dd_{*} as in Proposition 2.4. Proposition 2.5 asserts that the resulting error is negligible in the δ0\delta\to 0 limit.

To prove Proposition 2.5, we use the following bound on the dd\ell-measure of loops with large diameter.

Proposition 2.6.

The dd\ell-measure of loops \ell with diameter greater than KK tends to 0 as KK\rightarrow\infty faster than any negative power of KK.

Proof.

If \ell has diameter >K>K, then Proposition 2.1 implies that 0len()𝟏|(s)|>K/3𝑑t>cK\int_{0}^{\mathrm{len}(\ell)}\mathbf{1}_{|\ell(s)|>K/3}dt>cK for some constant c>0c>0. It follows from Markov’s inequality that the probability of this event is bounded from above by (cK)1𝔼(0len()𝟏|(s)|>K/3𝑑t)(cK)^{-1}\mathbb{E}(\int_{0}^{\mathrm{len}(\ell)}\mathbf{1}_{|\ell(s)|>K/3}dt). By the Schwartz condition in Definition 1.11, the latter decays as KK\rightarrow\infty faster than any negative power of KK. ∎

Proof of Proposition 2.5.

By definition of α\alpha, the integral in (2.8) gives the μ\mu-mass of loops LL with cen(L)D\mathrm{cen}(L)\in D and either

  • diam()δ<d\mathrm{diam}(\ell)\sqrt{\delta}<d_{*} and lenρ(L)δ\mathrm{len}_{\rho}(L)\geq\delta, or

  • diam()δd\mathrm{diam}(\ell)\sqrt{\delta}\geq d_{*} and tδt\geq\delta.

Thus, the error term—i.e., the difference between (2.8) and the mass of loops we consider in the proposition statement—is the μ\mu-mass of loops LL with diam()δd\mathrm{diam}(\ell)\sqrt{\delta}\geq d_{*} and for which either

  • tδt\geq\delta and lenρ(L)<δ\mathrm{len}_{\rho}(L)<\delta, or

  • t<δt<\delta and lenρ(L)δ\mathrm{len}_{\rho}(L)\geq\delta.

To bound this error, we recall from Proposition 2.1 that lenρ(L)\mathrm{len}_{\rho}(L) is >δ>\delta when t>Λδt>\Lambda\delta and <δ<\delta when t<Λ1δt<\Lambda^{-1}\delta. Thus, the error is at most the μ\mu-mass of loops LL with center in DD, diam()δd\mathrm{diam}(\ell)\sqrt{\delta}\geq d_{*} and Λ1δtΛδ\Lambda^{-1}\delta\leq t\leq\Lambda\delta. This mass is equal to the Lebesgue measure of DD times (ΛΛ1)/δ(\Lambda-\Lambda^{-1})/\delta times the dd\ell-measure of the set of loops \ell with diam()δ>d\mathrm{diam}(\ell)\sqrt{\delta}>d_{*}. From Proposition 2.6, we deduce that the error tends to 0 as δ0\delta\rightarrow 0 at a rate that is uniform in ρ\rho\in\mathcal{B}. ∎

Proposition 2.5 reduces the problem of proving Lemma 1.16 to the problem of analyzing the difference of integrands α1\alpha^{-1} and β1\beta^{-1} in (2.8) and (1.9). We first derive an expression for αβ\alpha-\beta in terms of the difference in the ρ\rho-length and center ρ\rho-length of the loop (x,β,)(x,\beta,\ell).

Proposition 2.7.

Let ε,K>0\varepsilon,K>0 be fixed. Then, for all \ell\in\mathcal{L} with diameter at most KK, as δ0\delta\to 0,

αβ=eρ(x)(lenρ(x,β,)clenρ(x,β,))+o(δ2),\alpha-\beta=e^{-\rho(x)}\Bigl{(}\mathrm{len}_{\rho}(x,\beta,\ell)-\mathrm{clen}_{\rho}(x,\beta,\ell)\Bigr{)}+o(\delta^{2}), (2.9)

with the error converging uniformly in the choice of ρ\rho\in\mathcal{B}, \ell with diameter at most KK, and xDDρ,εx\in D\setminus D^{\rho,\varepsilon} where {Dρ,ε}ρD\{D^{\rho,\varepsilon}\}_{\rho\in D} is defined as in Proposition 2.3. If we remove the restriction on diam()\mathrm{diam}(\ell), then the error is O(δ2diam()2)O(\delta^{2}\mathrm{diam}(\ell)^{2}) uniformly in ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L}..

Refer to caption
Figure 1: An illustration of the quantities we consider in Proposition 2.7. With ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L} fixed, the blue curve is the graph of tclenρ(L)t\mapsto\mathrm{clen}_{\rho}(L), and the red curve is the graph of tlenρ(L)t\mapsto\mathrm{len}_{\rho}(L). The blue curve is a line with slope eρ(x)e^{\rho(x)}, and the red curve is differentiable with the same derivative eρ(x)e^{\rho(x)} at the origin. Here, we consider δ>0\delta>0 for which diam()δ<d\mathrm{diam}(\ell)\sqrt{\delta}<d_{*}, with dd_{*} as in Proposition 2.4. By Proposition 2.4, the red curve intersects the horizontal (dotted) line of height δ\delta at the single red point (α,δ)(\alpha,\delta). The blue curve intersects the line of height δ\delta at the blue point (β,δ)(\beta,\delta). The green point is the point on the red curve directly above the blue point, and the purple line segment is the segment parallel to the blue line from the green point to the dotted line (purple point). Proposition 2.7 asserts that the distance between the blue and red points can be approximated by the distance between the blue and purple points. The latter distance is simply the length of the green segment divided by the slope eρ(x)e^{\rho(x)} of the blue line.
Proof.

Let d>0d_{*}>0 be as in Proposition 2.4. Observe that if we are restricting to \ell with diameter at most some constant KK, we automatically have diam()δ<d\mathrm{diam}(\ell)\sqrt{\delta}<d_{*} for uniformly small δ\delta. If we do not impose the restriction diam()K\mathrm{diam}(\ell)\leq K, then we could have diam()δ>d\mathrm{diam}(\ell)\sqrt{\delta}>d_{*} for arbitrarily small δ\delta. However, the proposition statement easily holds for this range of δ\delta and \ell: by Proposition 2.1, the error in the proposition statement must be bounded by a uniform constant times δ\delta, which is O(δ2diam()2)O(\delta^{2}\mathrm{diam}(\ell)^{2}) when diam()δ>d\mathrm{diam}(\ell)\sqrt{\delta}>d_{*} because δ2diam()2>d2δ\delta^{2}\mathrm{diam}(\ell)^{2}>d_{*}^{2}\delta. Thus, we may assume for the rest of the proof that diam()δ<d\mathrm{diam}(\ell)\sqrt{\delta}<d_{*}.

Throughout the proof, we refer to the graphs of the functions tclenρ(L)t\mapsto\mathrm{clen}_{\rho}(L) and tlenρ(L)t\mapsto\mathrm{len}_{\rho}(L) for fixed ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L} in Figure 1. See the caption of the figure for the definitions of the red, blue and green points. By Proposition 2.4, αβ\alpha-\beta is the distance between the red and blue points, and the quantity eρ(x)(lenρ(x,β,)clenρ(x,β,))e^{-\rho(x)}\left(\mathrm{len}_{\rho}(x,\beta,\ell)-\mathrm{clen}_{\rho}(x,\beta,\ell)\right) on the right-hand side of (2.9) is the distance between the blue and purple points. We can express these two distances in terms of the slopes of the two curves:

  • The distance between the blue and purple points is equal to the length of the green segment divided by the slope of the blue line (i.e., eρ(x)e^{\rho(x)}).

  • The distance between the blue and red points is equal to the length of the green segment divided by the average derivative of the red curve between the red and green points.

The error that we need to bound is the difference between these two distances—i.e., the distance between the red and purple points. By Proposition 2.3, the derivative of the red curve between the red and green points differs from the derivative of the blue line by at most cdiam()δc\,\mathrm{diam}(\ell)\sqrt{\delta}, where cc is a constant bounded uniformly in ρ\rho\in\mathcal{B} and xDDρ,εx\in D\setminus D^{\rho,\varepsilon} that can be made arbitrarily small for diam()δ\mathrm{diam}(\ell)\sqrt{\delta} sufficiently small. By Proposition 2.1, this implies that the inverses of these two derivatives differ by a uniform constant multiplied by cdiam()δc\,\mathrm{diam}(\ell)\sqrt{\delta}. Moreover, the length of the green segment is lenρ(x,β,)δ\mathrm{len}_{\rho}(x,\beta,\ell)-\delta, and by (2.6),

lenρ(x,β,)δ=β(A(β)eρ(x))cdiam()β3/2,\mathrm{len}_{\rho}(x,\beta,\ell)-\delta=\beta(A(\sqrt{\beta})-e^{\rho(x)})\leq c\,\mathrm{diam}(\ell)\beta^{3/2}, (2.10)

with cc as above. Thus, the error term—i.e., the distance between the red and purple points—is bounded by a uniform constant times cdiam()δcdiam()δ3/2=c2δ2diam()2c\,\mathrm{diam}(\ell)\sqrt{\delta}\cdot c\,\mathrm{diam}(\ell)\delta^{3/2}=c^{2}\delta^{2}\mathrm{diam}(\ell)^{2}. If we do not restrict to diam()K\mathrm{diam}(\ell)\leq K, the error is O(δ2diam()2)O(\delta^{2}\mathrm{diam}(\ell)^{2}) uniformly in ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L}. If we restrict to diam()K\mathrm{diam}(\ell)\leq K, then c0c\rightarrow 0 as δ0\delta\rightarrow 0 at a uniform rate, so the error is o(δ2)o(\delta^{2}) uniformly in ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell with diam()K\mathrm{diam}(\ell)\leq K. ∎

Proposition 2.7 immediately yields the following expression for α1β1\alpha^{-1}-\beta^{-1}.

Proposition 2.8.

Let ε,K>0\varepsilon,K>0 be fixed. Then, for all \ell with diameter at most KK, as δ0\delta\to 0,

α1β1=δ1β1(lenρ(x,β,)clenρ(x,β,))+o(1),\alpha^{-1}-\beta^{-1}=-\delta^{-1}\beta^{-1}\Bigl{(}\mathrm{len}_{\rho}(x,\beta,\ell)-\mathrm{clen}_{\rho}(x,\beta,\ell)\Bigr{)}+o(1), (2.11)

with the error converging uniformly in the choice of ρ\rho\in\mathcal{B}, \ell with diameter at most KK, and xDDρ,εx\in D\setminus D^{\rho,\varepsilon} where {Dρ,ε}ρD\{D^{\rho,\varepsilon}\}_{\rho\in D} is defined as in Proposition 2.3. If we remove the restriction on diam()\mathrm{diam}(\ell), then the error is O(diam()2)O(\mathrm{diam}(\ell)^{2}) + o(δ2)o(\delta^{2}) uniformly in ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L}.

Proof.

Throughout the following proof, each O()O(\cdot) and o()o(\cdot) error converges uniformly as δ0\delta\rightarrow 0 in the choice of ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell.

By the Taylor expansion of f(r)=1/rf(r)=1/r at β\beta, we have

α1β1=β2(αβ)+2γ3(αβ)2\alpha^{-1}-\beta^{-1}=-\beta^{-2}(\alpha-\beta)+2\gamma^{-3}(\alpha-\beta)^{2}

for some γ\gamma between α\alpha and β\beta. Proposition 2.7 gives

β2(αβ)=δ1β1(lenρ(x,β,)clenρ(x,β,))+o(1).-\beta^{-2}(\alpha-\beta)=-\delta^{-1}\beta^{-1}\Bigl{(}\mathrm{len}_{\rho}(x,\beta,\ell)-\mathrm{clen}_{\rho}(x,\beta,\ell)\Bigr{)}+o(1).

We now handle the 2γ3(αβ)22\gamma^{-3}(\alpha-\beta)^{2} term. From (2.1), we have Λ1α/δΛ\Lambda^{-1}\leq\alpha/\delta\leq\Lambda; therefore, γ3=O(δ3)\gamma^{-3}=O(\delta^{-3}). Next, by (2.10) and Proposition 2.7, (αβ)2=O(c2diam()2δ3)+o(δ4)(\alpha-\beta)^{2}=O(c^{2}\mathrm{diam}(\ell)^{2}\delta^{3})+o(\delta^{4}), where cc is a constant bounded uniformly in ρ\rho\in\mathcal{B} that can be made arbitrarily small for diam()δ\mathrm{diam}(\ell)\sqrt{\delta} sufficiently small. Hence, γ3(αβ)2=O(c2diam()2)+o(δ2)\gamma^{-3}(\alpha-\beta)^{2}=O(c^{2}\mathrm{diam}(\ell)^{2})+o(\delta^{2}). The latter is o(1)o(1) with the restriction on diam()\mathrm{diam}(\ell), and O(diam()2)+o(δ2)O(\mathrm{diam}(\ell)^{2})+o(\delta^{2}) otherwise. ∎

We also give a similar estimate for bookkeeping which might be useful for Question 6.5.

Proposition 2.9.

With the same assumption of Proposition 2.8 without restriction on diam()\mathrm{diam}(\ell), we have α2β2=O(cdiam()δ1/2)\alpha^{-2}-\beta^{-2}=O(c\,\mathrm{diam}(\ell)\delta^{-1/2}) uniformly in ρ\rho\in\mathcal{B}, xDDρ,εx\in D\setminus D^{\rho,\varepsilon}, and \ell\in\mathcal{L}, where cc is a constant bounded uniformly in ρ\rho\in\mathcal{B} that can be made arbitrarily small for diam()δ\mathrm{diam}(\ell)\sqrt{\delta} sufficiently small.

Proof.

By the Taylor expansion of f(r)=1/rf(r)=1/r at β\beta, we have α2β2=2γ3(αβ)\alpha^{-2}-\beta^{-2}=-2\gamma^{-3}(\alpha-\beta) for some γ\gamma between α\alpha and β\beta. From (2.1), we have Λ1α/δΛ\Lambda^{-1}\leq\alpha/\delta\leq\Lambda; therefore, γ3=O(δ3)\gamma^{-3}=O(\delta^{-3}). Next, by (2.10) and Proposition 2.7, αβ=O(cdiam()δ3/2)\alpha-\beta=O(c\,\mathrm{diam}(\ell)\delta^{3/2}), where cc is a constant bounded uniformly in ρ\rho\in\mathcal{B} that can be made arbitrarily small for diam()δ\mathrm{diam}(\ell)\sqrt{\delta} sufficiently small. Hence, γ3(αβ)=O(cdiam()δ1/2)\gamma^{-3}(\alpha-\beta)=O(c\,\mathrm{diam}(\ell)\delta^{-1/2}). ∎

Proof of Lemma 1.16.

Let ε>0\varepsilon>0 and {Dρ,ε}ρD\{D^{\rho,\varepsilon}\}_{\rho\in D} be defined as in Proposition 2.3 so that Vol(Dρ,ε)ε\mathrm{Vol}(D^{\rho,\varepsilon})\leq\varepsilon. By Propositions 1.15 and 2.5, integrating α1β1\alpha^{-1}-\beta^{-1} over xDDρ,εx\in D\setminus D^{\rho,\varepsilon} and \ell\in\mathcal{L} yields the μ\mu-mass of loops LL centered in DD with lenρ(L)δ\mathrm{len}_{\rho}(L)\geq\delta minus the μ\mu-mass of loops LL centered in DD with clenρ(L)δ\mathrm{clen}_{\rho}(L)\geq\delta, up to a uniform o(1)o(1) error. By Proposition 2.8, for each fixed xDDρ,εx\in D\setminus D^{\rho,\varepsilon} and \ell\in\mathcal{L}, the integrand α1β1\alpha^{-1}-\beta^{-1} is equal to

δ1(β1lenρ(x,eρ(x)δ,)eρ(x))-\delta^{-1}\Bigl{(}\beta^{-1}\mathrm{len}_{\rho}(x,e^{-\rho(x)}\delta,\ell)-e^{\rho(x)}\Bigr{)} (2.12)

plus an error term that has the following limiting behavior as δ0\delta\rightarrow 0:

  1. (a)

    The error is O(diam()2)+o(δ2)O(\mathrm{diam}(\ell)^{2})+o(\delta^{2}) uniformly in ρ\rho\in\mathcal{B} and (x,)(DDρ,ε)×(x,\ell)\in(D\setminus D^{\rho,\varepsilon})\times\mathcal{L}.

  2. (b)

    For each fixed K>0K>0, the error is o(1)o(1) uniformly in ρ\rho\in\mathcal{B} and (x,)(DDρ,ε)×(x,\ell)\in(D\setminus D^{\rho,\varepsilon})\times\mathcal{L} with diam()K\mathrm{diam}(\ell)\leq K.

We now integrate α1β1\alpha^{-1}-\beta^{-1} over xDDρ,εx\in D\setminus D^{\rho,\varepsilon} and \ell\in\mathcal{L} and take the δ0\delta\rightarrow 0 limit. The integral of (2.12) is exactly equal to the term 1δ𝔼[eρ(Z)eρ(X)]\frac{1}{\delta}\mathbb{E}[e^{\rho(Z)}-e^{\rho(X)}] in (1.10), so we just need to show that the integral of the error term in the expression for α1β1\alpha^{-1}-\beta^{-1} tends to 0 as δ0\delta\rightarrow 0 uniformly in ρ\rho\in\mathcal{B}.

To analyze the integral of this error term, we partition the domain of integration. For any function K(δ)K(\delta) of δ>0\delta>0, we can partition (DDρ,ε)×(D\setminus D^{\rho,\varepsilon})\times\mathcal{L} into two subdomains: the set of pairs (x,)(x,\ell) with diam()K(δ)\mathrm{diam}(\ell)\leq K(\delta), and the set of (x,)(x,\ell) with diam()>K(δ)\mathrm{diam}(\ell)>K(\delta). The bound (a) implies that the integral of the error over (x,)(x,\ell) with diam()>k\mathrm{diam}(\ell)>k equals a uniform constant times {diam()>k}(diam()2+δ2)𝑑\int_{\{\mathrm{diam}(\ell)>k\}}(\mathrm{diam}(\ell)^{2}+\delta^{2})d\ell, which tends to zero as kk\rightarrow\infty at a rate uniform in (small) δ>0\delta>0 and ρ\rho\in\mathcal{B}. Moreover, (b) implies that if k>0k>0 is fixed, the integral of the error over (x,)(x,\ell) with diam()k\mathrm{diam}(\ell)\leq k tends to 0 as δ0\delta\rightarrow 0 uniformly in ρ\rho\in\mathcal{B}. Hence, if we choose a function K(δ)K(\delta) that tends to infinity sufficiently slowly as δ0\delta\rightarrow 0, the integrals of the error term over both subdomains of (DDρ,ε)×(D\setminus D^{\rho,\varepsilon})\times\mathcal{L} tend to 0 as δ0\delta\rightarrow 0 uniformly in ρ\rho\in\mathcal{B}.

Finally, the integral of α1β1\alpha^{-1}-\beta^{-1} over xDρ,εx\in D^{\rho,\varepsilon} and \ell\in\mathcal{L} is bounded by a uniform constant times εδ\varepsilon\delta because of (2.1), which finishes the proof. ∎

3 Expected length discrepancy vs. Dirichlet energy

We now complete the proof of Theorem 1.12 by proving Lemma 1.17. The main step in proving this lemma is proving the following proposition.

Proposition 3.1.

With ZZ and XX defined as in Lemma 1.16, the quantity β1𝔼[ρ(Z)ρ(X)]\beta^{-1}\mathbb{E}[\rho(Z)-\rho(X)] converges to 0 as δ0\delta\rightarrow 0 uniformly in ρ\rho\in\mathcal{B}. The expectation here is w.r.t. to the overall law of XX and ZZ as defined in Lemma 1.16.

We prove Proposition 3.1 by a Fourier analysis argument. Throughout this section, we define the Fourier transform of a function ϕ\phi as ϕ^(ζ):=2eiζzϕ(z)𝑑z\widehat{\phi}(\zeta):=\int_{\mathbb{R}^{2}}e^{i\zeta\cdot z}\phi(z)\,dz (using the convention of characteristic functions). To prove Proposition 3.1, we express the conditional density of ZZ given X=xX=x in terms of the expected occupation measure of a loop sampled from μ\mu with given length and center. In the following proposition, we introduce some notation for this expected occupation measure and record some of its elementary properties.

Proposition 3.2.

For z2z\in\mathbb{R}^{2} and s>0s>0, let θ(z,s)\theta(z,s) be the expected occupation measure of a loop sampled from μ\mu and conditioned to have length ss and center zero. For fixed ss, the measure θ(z,s)\theta(z,s) is radially symmetric Schwartz, and a probability measure. Each coordinate has second central moment bsbs. Moreover, θ\theta satisfies the scaling relation θ(z/t,s)=tθ(z,ts)\theta(z/\sqrt{t},s)=t\theta(z,ts) for any t0t\geq 0. Finally, we have lims0θ^(z,s)=1\lim_{s\to 0}\widehat{\theta}(z,s)=1.

Proof.

Since θ(z,s)\theta(z,s) is the law of With \ell sampled from dd\ell, the loop ts(t/s)t\mapsto\sqrt{s}\ell(t/s) has law θ(z,s)\theta(z,s). It follows from Definition 1.11 that θ(z,s)\theta(z,s) is radially symmetric, Schwartz, and a probability measure with the desired second central moments. The scaling relation also follows immediately. Finally, the second central moments imply that the measures θ(z,s)\theta(z,s) converge weakly as s0s\rightarrow 0 to a point mass at 0, so their characteristic functions converge to 11. ∎

We first reduce the task of proving Proposition 3.1 to analyzing the Dirichlet energy of the inverse Laplacian of an appropriately chosen measure.

Proposition 3.3.

Suppose that ϕ:2\phi:\mathbb{R}^{2}\to\mathbb{R} is a function satisfying

Δϕ(y)=Dβ1(θ(yx,β)𝜹(yx))𝑑xandlim|y|ϕ(y)=0,\Delta\phi(y)=\int_{D}\beta^{-1}(\theta(y-x,\beta)-\bm{\delta}(y-x))\,dx\qquad\text{and}\qquad\lim_{\left|y\right|\to\infty}\phi(y)=0, (3.1)

where θ\theta is defined in Proposition 3.2 and 𝛅\bm{\delta} is the point mass at 0. Then

β1𝔼[ρ(Z)ρ(X)](ρ,ρ)2(ϕ,ϕ)2\beta^{-1}\mathbb{E}[\rho(Z)-\rho(X)]\leq(\rho,\rho)_{\nabla}^{2}(\phi,\phi)_{\nabla}^{2}
Proof.

The conditional density of ZZ given X=xX=x is θ(x,β)\theta(\cdot-x,\beta). Therefore,

β1𝔼[ρ(Z)ρ(X)]=D(ρ,β1(θ(x,β)𝜹(x)))dx=(ρ,Δϕ)\beta^{-1}\mathbb{E}[\rho(Z)-\rho(X)]=\int_{D}\Bigl{(}\rho,\beta^{-1}\bigl{(}\theta(\cdot-x,\beta)-\bm{\delta}(\cdot-x)\bigr{)}\Bigr{)}\,dx=(\rho,\Delta\phi)

by the definition of ϕ\phi. Integrating by parts (or applying Green’s identity) gives (ρ,Δϕ)=(ρ,ϕ)(\rho,\Delta\phi)=-(\rho,\phi)_{\nabla}, and

(ρ,ϕ)2(ρ,ρ)(ϕ,ϕ).(\rho,\phi)_{\nabla}^{2}\leq(\rho,\rho)_{\nabla}(\phi,\phi)_{\nabla}.

by applying the Cauchy-Schwarz inequality. ∎

Therefore, it is enough to show (ϕ,ϕ)0(\phi,\phi)_{\nabla}\to 0 uniformly for the proof. To express a function ϕ\phi satisfying (3.1) using Fourier integral444Alternatively, such ϕ\phi can be written in terms of the convolution with Green’s function., we first define a family of auxiliary functions that we label Θt\Theta_{t} indexed by t>0t>0.

Proposition 3.4.

With the notation θs^(ζ,s0):=.θ^s(ζ,s)|s=s0\widehat{\theta_{s}}(\zeta,s_{0}):=\bigl{.}\frac{\partial\widehat{\theta}}{\partial s}(\zeta,s)\bigr{|}_{s=s_{0}}, the function Θt\Theta_{t} for t>0t>0 defined as the inverse Fourier transform of |ζ|2θs^(ζ,t)\left|\zeta\right|^{-2}\widehat{\theta_{s}}(\zeta,t) is a radially symmetric Schwartz function. In addition, Θt(y)=Θ1(y/t)/t\Theta_{t}(y)=\Theta_{1}(y/\sqrt{t})/t and Θt^(ζ)\widehat{\Theta_{t}}(\zeta) is uniformly bounded for all t>0t>0 and ζ2\zeta\in\mathbb{R}^{2}.

Proof.

Fix t>0t>0. Recall that θ^(r,t)\widehat{\theta}(r,t) is Schwartz in r0r\geq 0, so we use the notation θr^(r,t):=.θ^|ζ|(|ζ|,t)||ζ|=r\widehat{\theta_{r}}(r,t):=\bigl{.}\frac{\partial\widehat{\theta}}{\partial\left|\zeta\right|}(\left|\zeta\right|,t)\Bigr{|}_{\left|\zeta\right|=r} in this proof. The scaling relation and the radial symmetry of θ\theta (Proposition 3.2) implies that Θt^(ζ)=|ζ|2θs^(ζ,t)=θs^(1,t|ζ|2)\widehat{\Theta_{t}}(\zeta)=\left|\zeta\right|^{-2}\widehat{\theta_{s}}(\zeta,t)=\widehat{\theta_{s}}(1,t\left|\zeta\right|^{2}). Hence, it is enough to show that Θt^(r)=θs^(1,tr2)\widehat{\Theta_{t}}(r)=\widehat{\theta_{s}}(1,tr^{2}), as a function of r0r\geq 0, is a Schwartz function. This will also imply that Θt^(ζ)=Θ1^(tζ)\widehat{\Theta_{t}}(\zeta)=\widehat{\Theta_{1}}(\sqrt{t}\zeta) is uniformly bounded by supt0|Θ1^(t)|\sup_{t\geq 0}\left|\widehat{\Theta_{1}}(t)\right|, and Θt(y)=Θ1(y/t)/t\Theta_{t}(y)=\Theta_{1}(y/\sqrt{t})/t from the Fourier inversion.

First we repeatedly differentiate both sides of the relation θ^(1,tr2)=θ^(r,t)\widehat{\theta}(1,tr^{2})=\widehat{\theta}(r,t) using the chain rule, and observe skθ^(1,tr2)r=0\partial_{s}^{k}\widehat{\theta}(1,tr^{2})\mid_{r=0} exists as a constant multiple of r2kθ^(r,t)r=0\partial_{r}^{2k}\widehat{\theta}(r,t)\mid_{r=0}. Also, the first differentiation implies

Θt^(r)=θs^(1,tr2)=θr^(r,t)2tr.\widehat{\Theta_{t}}(r)=\widehat{\theta_{s}}(1,tr^{2})=\frac{\widehat{\theta_{r}}(r,t)}{2tr}.

Let P(r)P(r) be any polynomial in rr, and kk be any nonnegative integer. On one hand, for a fixed ε>0\varepsilon>0, we have

supr>ε|P(r)dkdrkθr^(r,t)2tr|<.\sup_{r>\varepsilon}\left|P(r)\frac{d^{k}}{dr^{k}}\frac{\widehat{\theta_{r}}(r,t)}{2tr}\right|<\infty.

because θr^(r,t)\widehat{\theta_{r}}(r,t) and its derivatives with respect to rr are rapidly decreasing, On the other hand, note that |P(r)dkdrkθs^(1,tr2)|\left|P(r)\frac{d^{k}}{dr^{k}}\widehat{\theta_{s}}(1,tr^{2})\right| is continuous on (0,ε](0,\varepsilon] and bounded at 0. This proves that Θt\Theta_{t} is Schwartz. ∎

Proposition 3.5.

With Θt\Theta_{t} defined in Proposition 3.4, for x2x\in\mathbb{R}^{2}, let ϕx:2\phi^{x}:\mathbb{R}^{2}\to\mathbb{R} be defined as

ϕx(y):=1δ0δΘteρ(x)(yx)𝑑t.\phi^{x}(y):=-\frac{1}{\delta}\int_{0}^{\delta}\Theta_{te^{-\rho(x)}}(y-x)\,dt.

and ϕ:2\phi:\mathbb{R}^{2}\to\mathbb{R} defined as

ϕ(y):=Dϕx(y)𝑑x\phi(y):=\int_{D}\phi^{x}(y)\,dx (3.2)

for each y2y\in\mathbb{R}^{2}. Then ϕ\phi satisfies (3.1).

Proof.

By the Fourier transform, we have

ϕx^(ζ)\displaystyle\widehat{\phi^{x}}(\zeta) =1δ0δeiζxΘteρ(x)^(ζ)𝑑t\displaystyle=-\frac{1}{\delta}\int_{0}^{\delta}e^{i\zeta\cdot x}\widehat{\Theta_{te^{-\rho(x)}}}(\zeta)\,dt
=1δ0δeiζx|ζ|2θs^(ζ,teρ(x))𝑑t\displaystyle=-\frac{1}{\delta}\int_{0}^{\delta}e^{i\zeta\cdot x}\left|\zeta\right|^{-2}\widehat{\theta_{s}}\bigl{(}\zeta,te^{-\rho(x)}\bigr{)}\,dt
=eiζx|ζ|2β1(θ^(ζ,δeρ(x))1),\displaystyle=-e^{i\zeta\cdot x}\left|\zeta\right|^{-2}\beta^{-1}\bigl{(}\widehat{\theta}(\zeta,\delta e^{-\rho(x)})-1\bigr{)},

as we have defined β=δeρ(x)\beta=\delta e^{-\rho(x)}. Therefore, the Fourier inversion of Δϕx^=|ζ|2ϕx^\widehat{\Delta\phi^{x}}=-\left|\zeta\right|^{2}\widehat{\phi^{x}} gives the desired result. Furthermore, the first identity shows that ϕx^L1\widehat{\phi^{x}}\in L^{1} from Proposition 3.4. By the Riemann-Lebesgue lemma, we conclude that lim|y|ϕx(y)=0\lim_{\left|y\right|\to\infty}\phi^{x}(y)=0, and thus it follows from the dominated convergence theorem that lim|y|ϕ(y)=0\lim_{\left|y\right|\to\infty}\phi(y)=0. ∎

Proposition 3.6.

With ϕ\phi defined as (3.2) and Θt\Theta_{t} defined in Proposition 3.4, define

Gρ,ry(z):=e3ρ(yrz)/2Θ1(zeρ(yrz)/2)G_{\rho,r}^{y}(z):=e^{3\rho(y-rz)/2}\nabla\Theta_{1}(ze^{\rho(y-rz)/2}) (3.3)

for each ρ\rho\in\mathcal{B}, y2y\in\mathbb{R}^{2}, rr\in\mathbb{R}, and z2z\in\mathbb{R}^{2}. Then

yϕ(y)=1δ0δ1tyDtGρ,ty(z)𝑑z𝑑t.\nabla_{y}\phi(y)=\frac{1}{\delta}\int_{0}^{\delta}\frac{1}{\sqrt{t}}\int_{\frac{y-D}{\sqrt{t}}}G_{\rho,\sqrt{t}}^{y}(z)\,dz\,dt. (3.4)
Proof.

From Proposition 3.4,

ϕ(y)=1δ0δΘteρ(x)(yx)𝑑x𝑑t=1δ0δDΘ1((yx)(teρ(x))12)teρ(x)𝑑x𝑑t,\phi(y)=-\frac{1}{\delta}\int_{0}^{\delta}\Theta_{te^{-\rho(x)}}(y-x)\,dx\,dt=-\frac{1}{\delta}\int_{0}^{\delta}\int_{D}\frac{\Theta_{1}\bigl{(}(y-x)(te^{-\rho(x)})^{-\frac{1}{2}}\bigr{)}}{te^{-\rho(x)}}\,dx\,dt,

taking gradient and substituting yx=tzy-x=\sqrt{t}z gives

yϕ(y)=1δ0δ1tyDte3ρ(ytz)/2Θ1(zeρ(ytz)/2)𝑑z𝑑t,\nabla_{y}\phi(y)=\frac{1}{\delta}\int_{0}^{\delta}\frac{1}{\sqrt{t}}\int_{\frac{y-D}{\sqrt{t}}}e^{3\rho(y-\sqrt{t}z)/2}\nabla\Theta_{1}(ze^{\rho(y-\sqrt{t}z)/2})\,dz\,dt,

so the result follows. ∎

Proposition 3.7.

For each ρ\rho\in\mathcal{B}, rr\in\mathbb{R}, and y2y\in\mathbb{R}^{2}, let

Fρ,r(y):=2Gρ,ry(z)𝑑zF_{\rho,r}(y):=\int_{\mathbb{R}^{2}}G_{\rho,r}^{y}(z)\,dz (3.5)

where Gρ,ryG_{\rho,r}^{y} is defined as (3.3). Then Fρ,rF_{\rho,r} and rFρ,r\partial_{r}F_{\rho,r} both converges to 0 in L2(2)L^{2}(\mathbb{R}^{2}) as r0r\to 0 uniformly over ρ\rho\in\mathcal{B}.

Proof.

Since Θ1\nabla\Theta_{1} is rapidly decreasing, given any ε>0\varepsilon>0, there exists K=K(ϵ)K=K(\epsilon) such that

|z|K|Gρ,ry(z)|2𝑑z<ε\int_{\left|z\right|\geq K}\left|G_{\rho,r}^{y}(z)\right|^{2}\,dz<\varepsilon (3.6)

for all rr\in\mathbb{R}, y2y\in\mathbb{R}^{2}, and ρ\rho\in\mathcal{B} as ρ\rho is uniformly bounded. By radial symmetry, note that Fρ,0(y)=0F_{\rho,0}(y)=0 for all ρ\rho and yy because Θ1\nabla\Theta_{1} is an odd function. Now we divide the domain of integral (3.5) into two regions depending on the sign of z1z2z_{1}z_{2} when we write z=(z1,z2)2z=(z_{1},z_{2})\in\mathbb{R}^{2}, so that the integral on each domain becomes

±z1z2>0,|z|<KGρ,ry(z)𝑑z=z1>0,±z2>0,|z|<K(Gρ,ry(z)Gρ,ry(z))𝑑z.\int_{\pm z_{1}z_{2}>0,\left|z\right|<K}G_{\rho,r}^{y}(z)\,dz=\int_{z_{1}>0,\pm z_{2}>0,\left|z\right|<K}\bigl{(}G_{\rho,r}^{y}(z)-G_{\rho,-r}^{y}(z)\bigr{)}\,dz.

Since Θ1\nabla\Theta_{1} is Lipschitz as being Schwartz, the uniform equicontinuity of ρ\rho in L2(D)L^{2}(D) (in the sense of Remark 1.7) implies that Gρ,rG_{\rho,r}, as a function of rzrz, is also uniformly equicontinuous in L2(2)L^{2}(\mathbb{R}^{2}). Therefore, denoting A±={z2:z1>0,±z2>0,|z|<K}A_{\pm}=\{z\in\mathbb{R}^{2}:z_{1}>0,\pm z_{2}>0,\left|z\right|<K\}, there exists d=d(ε)>0d=d(\varepsilon)>0 such that

D||z|<KGρ,ry(z)|2𝑑y\displaystyle\int_{D}\left|\int_{\left|z\right|<K}G_{\rho,r}^{y}(z)\right|^{2}dy i{±}AiD|Gρ,ry(z)Gρ,ry(z)|2𝑑y𝑑z2ε\displaystyle\leq\sum_{i\in\{\pm\}}\int_{A_{i}}\int_{D}\left|G_{\rho,r}^{y}(z)-G_{\rho,-r}^{y}(z)\right|^{2}\,dy\,dz\leq 2\varepsilon

for all ρ\rho\in\mathcal{B} whenever |r|<d\left|r\right|<d. Combined with (3.6), we conclude that Fρ,rF_{\rho,r} converges to 0 in L2(D)L^{2}(D) uniformly over ρ\rho\in\mathcal{B} as r0r\to 0.

A similar argument applies to rFρ,r\partial_{r}F_{\rho,r}. In particular, as a function of ρ\rho and ρ\nabla\rho, note that Gρ\nabla G_{\rho} is Lipschitz. Hence the uniform equicontinuity of ρ\rho in W1,2(D)W^{1,2}(D) (in the sense of Remark 1.7) implies the uniform equicontinuity of Gρ,ry(z)\nabla G_{\rho,r}^{y}(z) in L2(2)L^{2}(\mathbb{R}^{2}) over ρ\rho\in\mathcal{B} as a function of rzrz. Since Θ1\Theta_{1} is rapidly decreasing, from integration by parts,

rFρ,r(y)|r=0=2zGρ,0y(z)z𝑑z=22Gρ,0y(z)𝑑z=0\left.\frac{\partial}{\partial r}F_{\rho,r}(y)\right|_{r=0}=-\int_{\mathbb{R}^{2}}\nabla_{z}G_{\rho,0}^{y}(z)\cdot z\,dz=\int_{\mathbb{R}^{2}}2G_{\rho,0}^{y}(z)\,dz=0

again by the oddity of Gρ,0yG_{\rho,0}^{y}. Repeating the previous argument, we conclude that rFρ,r\partial_{r}F_{\rho,r} converges to 0 in L2(2)L^{2}(\mathbb{R}^{2}) uniformly over ρ\rho\in\mathcal{B} as r0r\to 0. ∎

Proposition 3.8.

For each y2y\in\mathbb{R}^{2}, define

Φ(y):=1δ0δFρ,t(y)t𝑑t.\Phi(y):=\frac{1}{\delta}\int_{0}^{\delta}\frac{F_{\rho,\sqrt{t}}(y)}{\sqrt{t}}\,dt. (3.7)

Then ϕΦ\nabla\phi-\Phi converges to 0 in L2(2)L^{2}(\mathbb{R}^{2}) uniformly over ρ\rho\in\mathcal{B} as δ0\delta\to 0.

Proof.

Given ε>0\varepsilon>0, let K=K(ε)K=K(\varepsilon) be defined such that (3.6) holds. Choose d(ε)d(\varepsilon) small enough so that as long as 0<r<d0<r<d, we have |Gρ,ry(z)Gρ,0y(z))|<ε/πK2\left|G_{\rho,r}^{y}(z)-G_{\rho,0}^{y}(z))\right|<\varepsilon/\pi K^{2} and the diameter of (yD)/r(y-D)/r is greater than KK. Using the triangle inequality, we obtain the desired result from (3.4) and (3.7). ∎

Proof of Proposition 3.1.

By the mean value theorem, with Φ\Phi is defined in Proposition 3.8, we have

D|Φ(y)|2𝑑y1δ0δD|Fρ,t(y)t|2𝑑y𝑑t=D|rFρ,r(y)|2𝑑y\int_{D}\left|\Phi(y)\right|^{2}\,dy\leq\frac{1}{\delta}\int_{0}^{\delta}\int_{D}\left|\frac{F_{\rho,\sqrt{t}}(y)}{\sqrt{t}}\right|^{2}\,dy\,dt=\int_{D}\left|\partial_{r}{F_{\rho,r}(y)}\right|^{2}\,dy

for some r(0,δ)r\in(0,\sqrt{\delta}), that is Φ\Phi converges to 0 in L2(2)L^{2}(\mathbb{R}^{2}) as δ0\delta\to 0 uniformly over ρ\rho\in\mathcal{B} by Proposition 3.7. With Proposition 3.8, we conclude ϕ\nabla\phi also converges to 0 in L2(2)L^{2}(\mathbb{R}^{2}) as δ0\delta\to 0 uniformly over ρ\rho\in\mathcal{B}. This completes the proof. ∎

Proof of Lemma 1.17.

Let Y=ρ(Z)ρ(X)Y=\rho(Z)-\rho(X), so that

δ1𝔼[eρ(Z)eρ(X)]\displaystyle\delta^{-1}\mathbb{E}[e^{\rho(Z)}-e^{\rho(X)}] =δ1𝔼[eρ(X)(Y+Y22!+Y33!+)]\displaystyle=\delta^{-1}\mathbb{E}\left[e^{\rho(X)}\left(Y+\frac{Y^{2}}{2!}+\frac{Y^{3}}{3!}+\ldots\right)\right]
=δ1𝔼[eρ(X)Y]+δ1𝔼[eρ(X)(Y22!+O(Y3))])\displaystyle=\delta^{-1}\mathbb{E}[e^{\rho(X)}Y]+\delta^{-1}\mathbb{E}\left[e^{\rho(X)}\left(\frac{Y^{2}}{2!}+O(Y^{3})\right)\right]\Bigr{)} (3.8)

By Proposition 3.1, the term δ1𝔼[eρ(X)Y]\delta^{-1}\mathbb{E}[e^{\rho(X)}Y] is o(1)o(1) with the convergence uniform in ρ\rho\in\mathcal{B}. We now analyze the second term in (3.8). Since the gradients ρ\nabla\rho for ρ\rho\in\mathcal{B} are uniformly equicontinuous in L1L^{1} (in the sense of Remark 1.7), we can express YY as ρ(X)(ZX)+ε\nabla\rho(X)\cdot(Z-X)+\varepsilon almost surely, where ε\varepsilon is an error term uniformly o(|ZX|)o(|Z-X|) in expectation, as |ZX|0|Z-X|\to 0. Conditional on XX, each coordinate of ZXZ-X has mean XX and variance βb\beta b. Therefore, with UUniform[0,2π]U\sim\text{Uniform}[0,2\pi],

𝔼((ρ(X)(ZX))2|X)=|ρ(X)|2𝔼[|ZX|2|X]𝔼(cos2(U))=βb|ρ(X)|2.\mathbb{E}((\nabla\rho(X)\cdot(Z-X))^{2}|X)=|\nabla\rho(X)|^{2}\mathbb{E}\left[|Z-X|^{2}|X\right]\mathbb{E}(\cos^{2}(U))=\beta b|\nabla\rho(X)|^{2}.

Therefore,

𝔼[Y2|X]\displaystyle\mathbb{E}[Y^{2}|X] =𝔼[(ρ(X)(ZX)+ε)2|X]\displaystyle=\mathbb{E}\left[(\nabla\rho(X)\cdot(Z-X)+\varepsilon)^{2}|X\right]
=βb|ρ(X)|2+𝔼[2ερ(X)(ZX)+ε2|X]\displaystyle=\beta b|\nabla\rho(X)|^{2}+\mathbb{E}[2\varepsilon\nabla\rho(X)\cdot(Z-X)+\varepsilon^{2}|X]
=βb|ρ(X)|2+o(δ),\displaystyle=\beta b|\nabla\rho(X)|^{2}+o(\delta),

with the o(δ)o(\delta) error uniform in ρ\rho and XX. We can similarly show that 𝔼[Y3|X]\mathbb{E}[Y^{3}|X] is o(δ)o(\delta) uniformly in ρ\rho and XX. Multiplying both sides by δ1eρ(X)\delta^{-1}e^{\rho(X)} and taking the expectation, we deduce that the second term in (3.8) is equal to b2(ρ,ρ)\frac{b}{2}(\rho,\rho)_{\nabla} plus an o(1)o(1) error that converges uniformly in ρ\rho\in\mathcal{B}. ∎

Proof of Theorem 1.12.

The result immediately follows from combining Proposition 1.15, Lemma 1.16 and Lemma 1.17. ∎

Finally, to deduce Theorem 1.4 from Theorem 1.12, we explicitly characterize the expected occupation measure.

Proposition 3.9.

For the Brownian loop measure, the expected occupation measure of a loop sampled from (,d)(\mathcal{L},d\ell) is the density of a complex Gaussian with variance 1/121/12.

Proof.

The law of a loop sampled from (,d)(\mathcal{L},d\ell) is that of a Brownian bridge indexed by the circle minus its mean. The value of a Brownian bridge indexed by the circle at any given time minus the mean value is a complex mean-zero Gaussian random variable of variance 1/121/12; this calculation appears, for example, in [She07].555One way to see it is to consider a Gaussian free field hh indexed by the circle and observe that the Dirichlet energy on the circle parameterized by [1/2,1/2][-1/2,1/2] of the function f(x)=x2/2f(x)=x^{2}/2 is given by 1/121/12, so (f,f)=1/12(f,f)_{\nabla}=1/12. But for a function gg on the circle we have from integration by parts that (f(x),g(x))=g(0)1/21/2g(x)𝑑x(f(x),g(x))_{\nabla}=g(0)-\int_{-1/2}^{1/2}g(x)dx. Then using the above and the definition of the GFF we have Var((h,f))=Var(h(1/2)01h(x)𝑑x)=(f,f)=1/12\mathrm{Var}\bigl{(}(h,f)_{\nabla}\bigr{)}=\mathrm{Var}\bigl{(}h(1/2)-\int_{0}^{1}h(x)dx\bigr{)}=(f,f)_{\nabla}=1/12. By rotational symmetry, this holds if 0 is replaced by any other number in [1/2,1/2][-1/2,1/2]. The number 1/121/12 is also derived in [She07] by Fourier series, and we remark that comparing these two derivations is one way to prove n=1n2=π2/6\sum_{n=1}^{\infty}n^{-2}=\pi^{2}/6.

Proof of Theorem 1.12.

The result follows from Theorem 1.12 and Proposition 3.9. ∎

4 Brownian loops on surfaces

In this section, we prove Theorem 1.9. Throughout the section, we let (M,g)(M,g) be a fixed compact smooth two-dimensional Riemannian manifold, and we let μloop\mu^{\text{loop}} denote the Brownian loop measure on (M,g)(M,g). We also fix \mathcal{B} as a precompact set of Lipschitz functions in W1,1(M)W^{1,1}(M) with uniformly bounded Lipschitz constants.

To prove Theorem 1.9, we analyze the mass of “large” loops and the mass of “small” loops with respect to ρ\rho-length separately. We begin by analyzing the mass of large loops by proving a central limit theorem for ρ\rho-length along large loops:

Proposition 4.1.

We can choose a constant c>0c>0 such that, for every ε,t>0\varepsilon,t>0 and ρ\rho\in\mathcal{B} and every loop LL sampled from μ(z,z;t)\mu(z,z;t),

|0teρ(L(s))𝑑stVolρ(M)/Vol(M)|ct(1+ε)/2\left|\int_{0}^{t}e^{\rho(L(s))}ds-t\mathrm{Vol}_{\rho}(M)/\mathrm{Vol}(M)\right|\leq ct^{(1+\varepsilon)/2} (4.1)

with probability 1O(etε)1-O(e^{-t^{\varepsilon}}) as tt\to\infty, with the rate uniform in ρ\rho\in\mathcal{B}.

To prove Proposition 4.1, we apply the Markov central limit theorem to Brownian motion on (M,g)(M,g), and then compare Brownian motion to a loop sampled from μ(z,z;t)\mu(z,z;t) by using the following Radon-Nikodym estimate.

Proposition 4.2.

Let zMz\in M be fixed, and let LL be a loop sampled from μ(z,z;t)\mu(z,z;t). The Radon-Nikodym derivative of the law of L|[0,s]L|_{[0,s]} with respect to Brownian motion restricted to [0,s][0,s] is given by 1+O(eα(ts))1+O(e^{-\alpha(t-s)}) for some α=α(M,g)>0\alpha=\alpha(M,g)>0.

Proof.

By Proposition 4.3, μ(,z;ts)/μ(z,z;t)=1+O(eα(ts))\|\mu(\cdot,z;t-s)/\mu(z,z;t)\|_{\infty}=1+O(e^{-\alpha(t-s)}) for some α=α(M,g)>0\alpha=\alpha(M,g)>0, with the error uniform in zMz\in M. This implies the derivative bound. ∎

To apply the Markov central limit theorem to Brownian motion on (M,g)(M,g), we need the following convergence result for the Brownian transition kernel μ(,;t)\mu(\cdot,\cdot;t) as tt\to\infty.

Proposition 4.3.

We have

μ(,;t)(Vol(M))1L(M×M)=O(eαt).\|\mu(\cdot,\cdot;t)-(\mathrm{Vol}(M))^{-1}\|_{L^{\infty}(M\times M)}=O(e^{-\alpha t}).

for some constant α>0\alpha>0.

Proof.

Let 0=λ0λ1λ20=\lambda_{0}\leq\lambda_{1}\leq\lambda_{2}\leq\cdots denote the eigenvalues of Δ\Delta, and let {ϕn}\{\phi_{n}\} be a corresponding Hilbert basis of eigenfunctions. It follows from the heat equation that pp can be written as

p(x,y;t)=n=0etλnϕn(x)ϕn(y).p(x,y;t)=\sum_{n=0}^{\infty}e^{-t\lambda_{n}}\phi_{n}(x)\phi_{n}(y).

It is known that ϕncλn1/4\|\phi_{n}\|_{\infty}\leq c\lambda_{n}^{1/4} for some constant cc depending only on (M,g)(M,g) [Don01]. Hence, μ(,;t)(Vol(M))1L(M×M)\|\mu(\cdot,\cdot;t)-(\mathrm{Vol}(M))^{-1}\|_{L^{\infty}(M\times M)} is bounded from above by a function of tt that decays exponentially as tt\to\infty. ∎

Proof of Proposition 4.1.

(In the proof that follows, the rate of convergence of the O()O(\cdot) errors are uniform in the choice of ρ\rho\in\mathcal{B}.) Let BtB_{t} be a Brownian motion on (M,g)(M,g) started at a point sampled from the volume measure associated to (M,g)(M,g). Let Yn=B|[n,n+1]Y_{n}=B|_{[n,n+1]}, and let f(Yn)=nn+1eρ(Bt)𝑑tf(Y_{n})=\int_{n}^{n+1}e^{\rho(B_{t})}dt. For yMy\in M, let f(y)f_{*}(y) be the expected value of nn+1eρ(Bty)𝑑t\int_{n}^{n+1}e^{\rho(B^{y}_{t})}dt, where BtyB_{t}^{y} is a Brownian motion on (M,g)(M,g) started at yy. By Proposition 4.3, the conditional expectation of f(Yn)f(Y_{n}) given Y0Y_{0} is Mp(B1,y;n1)f(y)Vol(dy)=𝔼(f(Y0))+O(eαn)\int_{M}p(B_{1},y;n-1)f_{*}(y)\mathrm{Vol}(dy)=\mathbb{E}(f(Y_{0}))+O(e^{-\alpha n}) for some α=α(M,g)\alpha=\alpha(M,g). Hence, Cov(f(Y0),f(Yn))=O(eαn)\mathrm{Cov}(f(Y_{0}),f(Y_{n}))=O(e^{-\alpha n}) for some α=α(M,g)\alpha=\alpha(M,g). Thus, we may apply the Markov chain central limit theorem to deduce that

n(1nj=1nf(Yj)𝔼(f(Y0)))\sqrt{n}\left(\frac{1}{n}\sum_{j=1}^{n}f(Y_{j})-\mathbb{E}(f(Y_{0}))\right)

converges in the nn\to\infty limit to a centered Gaussian distribution whose variance is bounded uniformly in the choice of ρ\rho\in\mathcal{B}. Note that we have 𝔼(f(Y0))=1Vol(M)Meρ(x))Vol(dx)=Volρ(M)/Vol(M)\mathbb{E}(f(Y_{0}))=\frac{1}{\mathrm{Vol}(M)}\int_{M}e^{\rho(x)})\mathrm{Vol}(dx)=\mathrm{Vol}_{\rho}(M)/\mathrm{Vol}(M). Therefore, for each fixed ε,c>0\varepsilon,c>0,

|0neρ(Bt)𝑑tnVolρ(M)/Vol(M)|cn(1+ε)/2\left|\int_{0}^{n}e^{\rho(B_{t})}dt-n\mathrm{Vol}_{\rho}(M)/\mathrm{Vol}(M)\right|\leq cn^{(1+\varepsilon)/2}

with probability 1O(enε)1-O(e^{-n^{\varepsilon}}).

Combining this with Proposition 4.2, we deduce that, if nn is an integer with t2tnttt-2\sqrt{t}\leq n\leq t-\sqrt{t}, then

|0neρ(L(s))𝑑snVolρ(M)/Vol(M)|ct(1+ε)/2\left|\int_{0}^{n}e^{\rho(L(s))}ds-n\mathrm{Vol}_{\rho}(M)/\mathrm{Vol}(M)\right|\leq c\sqrt{t}^{(1+\varepsilon)/2}

with probability 1O(etε)1-O(e^{-t^{\varepsilon}}) as tt\to\infty. Since nteρ(L(s))𝑑s\int_{n}^{t}e^{\rho(L(s))}ds is bounded from above by t\sqrt{t} times a constant uniform in ρ\rho\in\mathcal{B}, this implies the proposition. ∎

We now apply Proposition 4.1 to analyze the mass of large loops.

Proposition 4.4.

For each ε>0\varepsilon>0, the symmetric difference between

  • the set of loops with length δ\geq\delta and length C\leq C , and

  • the set of loops with length δ\geq\delta and ρ\rho-length CVolρ(M)/Vol(M)\leq C\mathrm{Vol}_{\rho}(M)/\mathrm{Vol}(M),

has μloop\mu^{\text{loop}}-mass at most O(C(1+ε)/2)O(C^{(-1+\varepsilon)/2}), with the rate of convergence uniform in ρ\rho\in\mathcal{B}.

Proof.

(In the proof that follows, the rate of convergence of the O()O(\cdot) errors are uniform in the choice of ρ\rho\in\mathcal{B}.) Let 𝒮\mathcal{S} denote the symmetric difference between the two sets. By Proposition 2.1, each loop in 𝒮\mathcal{S} has length between Λ1C\Lambda^{-1}C and ΛC\Lambda C, with Λ\Lambda independent of the choice of ρ\rho\in\mathcal{B}. By Proposition 4.1, the subset 𝒮𝒮\mathcal{S}^{*}\subset\mathcal{S} of loops not satisfying (4.1) has μloop\mu^{\text{loop}}-mass at most Λ1CΛCO(etε/t)𝑑t=O(exp(ΛεCε)/C)\int_{\Lambda^{-1}C}^{\Lambda C}O(e^{-t^{\varepsilon}}/t)dt=O(\exp(-\Lambda^{-\varepsilon}C^{\varepsilon})/C). Moreover, with a±=CVolρ(M)/Vol(M)±c(C(1+ε)/2a_{\pm}=C\mathrm{Vol}_{\rho}(M)/\mathrm{Vol}(M)\pm c(C^{(1+\varepsilon)/2}, we can bound the μloop\mu^{\text{loop}}-mass of the set of loops in 𝒮\𝒮\mathcal{S}\backslash\mathcal{S}^{*} by

aa+Mt1μρ(z,z;t)Volρ(dz)𝑑tVolρ(M)aa+t1𝑑t=O(C(1+ε)/2).\int_{a_{-}}^{a_{+}}\int_{M}t^{-1}\mu_{\rho}(z,z;t)\mathrm{Vol}_{\rho}(dz)dt\leq\mathrm{Vol}_{\rho}(M)\int_{a_{-}}^{a_{+}}t^{-1}dt=O(C^{(-1+\varepsilon)/2}).\qed

Next, we analyze the mass of small loops.

Proposition 4.5.

For each ρ1,ρ2\rho_{1},\rho_{2}\in\mathcal{B}, the difference in the masses of the sets

{L:lenρjLδ,lenLC}j=1,2\{L:\mathrm{len}_{\rho_{j}}L\geq\delta,\mathrm{len}L\leq C\}\qquad j=1,2

under the Brownian loop measure in (M,g)(M,g) is equal to the difference in the expressions

Volρj(M)2πδ+148πM(ρj2+2Kρj)Vol(dz),j=1,2.\frac{\mathrm{Vol}_{\rho_{j}}(M)}{2\pi\delta}+\frac{1}{48\pi}\int_{M}(\|\nabla\rho_{j}\|^{2}+2K\rho_{j})\mathrm{Vol}(dz),\quad j=1,2. (4.2)

plus a term that tends to zero as δ0\delta\to 0 at a rate that is uniform in the choice of ρ1,ρ2\rho_{1},\rho_{2}\in\mathcal{B}.

To prove Proposition 4.5, we will apply the following pair of propositions.

Proposition 4.6.

Let (M,g)(M,g) be a smooth two-dimensional Riemannian manifold. Let UMU\subset M, and let 𝒮\mathcal{S} be a collection of loops in (M,g)(M,g) that intersects a closed set disjoint from U¯\overline{U}. Then the following holds for all δ>0\delta>0 sufficiently small. Let ρ1,ρ2\rho_{1},\rho_{2}\in\mathcal{B} be uniformly bounded functions that agree outside UU. Then, for each L𝒮L\in\mathcal{S}, we have lenρ1Lδ\mathrm{len}_{\rho_{1}}L\geq\delta iff lenρ2Lδ\mathrm{len}_{\rho_{2}}L\geq\delta.

Proof.

It is straightforward as \mathcal{B} is uniformly bounded and the length of LL outside UU is bounded below. ∎

Proposition 4.7.

Let ff be a compactly supported function on a region D2D\subset\mathbb{R}^{2} with finite Dirichlet energy, and let hh be a smooth compactly supported function on DD, and let h,νh,Kh\nabla_{h},\nu_{h},K_{h} denote the gradient, volume form and Gaussian curvature associated to (D,eh|dz|2)(D,e^{h}|dz|^{2}). Then

D(f+h)2𝑑zDh2𝑑z=D(hf2+2Khf)νh(dz).\int_{D}\|\nabla(f+h)\|^{2}dz-\int_{D}\|\nabla h\|^{2}dz=\int_{D}(\|\nabla_{h}f\|^{2}+2K_{h}f)\nu_{h}(dz).
Proof.

It follows from

D((f+h)2h2)𝑑z=D(fΔf+2fΔh)𝑑z\displaystyle\int_{D}(\|\nabla(f+h)\|^{2}-\|\nabla h\|^{2})dz=-\int_{D}(f\Delta f+2f\Delta h)dz
=\displaystyle= D(fehΔf+2fehΔh)νh(dz)=D(fΔhf2Khf)νh(dz)\displaystyle-\int_{D}(fe^{-h}\Delta f+2fe^{-h}\Delta h)\nu_{h}(dz)=-\int_{D}(f\Delta_{h}f-2K_{h}f)\nu_{h}(dz)
=\displaystyle= D(fehΔf2Khf)νh(dz)=D(hf2+2Khf)νh(dz).\displaystyle-\int_{D}(fe^{-h}\Delta f-2K_{h}f)\nu_{h}(dz)=\int_{D}(\|\nabla_{h}f\|^{2}+2K_{h}f)\nu_{h}(dz).\qed
Proof of Proposition 4.5.

In the proof that follows, we consider ρ\rho-lengths and ρ\rho-volume forms (as defined in Definition 1.3) for functions ρ\rho on both (M,g)(M,g) and a region of the Euclidean plane. To avoid confusion between the two settings, we use the notation lenρ\mathrm{len}_{\rho} and Volρ\mathrm{Vol}_{\rho} in the (M,g)(M,g) setting, and len~ρ\widetilde{\mathrm{len}}_{\rho} and Vol~ρ\widetilde{\mathrm{Vol}}_{\rho} in the Euclidean setting.

Let U,V,WU,V,W be open sets in MM with U¯V\overline{U}\subset V and V¯W\overline{V}\subset W, such that we can find a homeomorphism φ:WW~2\varphi:W\rightarrow\widetilde{W}\subset\mathbb{R}^{2}. By a partition of unity argument, it suffices to prove the proposition under the assumption that ρ1,ρ2\rho_{1},\rho_{2} agree on UU. So, we assume that this is the case. By Proposition 4.6, the difference in the masses of the sets

{L:lenρjLδ,lenLC}j=1,2\{L:\mathrm{len}_{\rho_{j}}L\geq\delta,\mathrm{len}L\leq C\}\qquad j=1,2 (4.3)

under the Brownian loop measure in (M,g)(M,g) is equal to the difference in masses with (4.3) replaced by

{L:lenρjLδ,lenLC,LV}j=1,2.\{L:\mathrm{len}_{\rho_{j}}L\geq\delta,\mathrm{len}L\leq C,L\subset V\}\qquad j=1,2. (4.4)

Since, by Proposition 2.1, lenL>C\mathrm{len}L>C automatically implies lenρjLδ\mathrm{len}_{\rho_{j}}L\geq\delta for sufficiently small δ\delta (in a manner that does not depend on the choice of ρ1,ρ2\rho_{1},\rho_{2}\in\mathcal{B}), we can replace (4.4) by the sets

{L:lenρjLδ,LV},j=1,2,\{L:\mathrm{len}_{\rho_{j}}L\geq\delta,L\subset V\},\qquad j=1,2, (4.5)

with the condition lenLC\mathrm{len}L\leq C in (4.4) removed. Let ρ~j\widetilde{\rho}_{j} be the pushforward of ρj\rho_{j} via φ\varphi, and let eσ|dz|2e^{\sigma}|dz|^{2} be the pushforward of the metric gg via φ\varphi. Also, set V~=φ(V)\widetilde{V}=\varphi(V). By conformal invariance of the Brownian loop measure [APPS20, Lemma 3.3], the masses of the sets (4.5) under the Brownian loop measure in (M,g)(M,g) equals the masses of the sets

{L:len~ρ~j+σLδ,LV~}j=1,2\{L:\widetilde{\mathrm{len}}_{\widetilde{\rho}_{j}+\sigma}L\geq\delta,L\subset\widetilde{V}\}\qquad j=1,2 (4.6)

under the Brownian loop measure in the Euclidean plane. Now, let θ\theta be a smooth function on 2\mathbb{R}^{2} that equals 11 on V~¯\overline{\widetilde{V}} and zero outside a compact subset of W~\widetilde{W}. Since θ1\theta\equiv 1 on V~\widetilde{V}, (4.6) is equal to the set

{L:len~θ(ρ~j+σ)Lδ,LV~}j=1,2\{L:\widetilde{\mathrm{len}}_{\theta(\widetilde{\rho}_{j}+\sigma)}L\geq\delta,L\subset\widetilde{V}\}\qquad j=1,2 (4.7)

By Proposition 4.6, the difference in these loop masses is unchanged if we replace the sets (4.7) by the sets

{L:len~θ(ρ~j+σ)Lδ,cen(L)W~},j=1,2.\{L:\widetilde{\mathrm{len}}_{\theta(\widetilde{\rho}_{j}+\sigma)}L\geq\delta,\mathrm{cen}(L)\in\widetilde{W}\},\qquad j=1,2.

By Theorem 1.4, this difference in loop masses is given by the difference in the quantities

Vol~θ(ρ~j+σ)(W~)2πδ+148π148πW~(θ(ρj+σ))2𝑑z,j=1,2.\frac{\widetilde{\mathrm{Vol}}_{\theta(\widetilde{\rho}_{j}+\sigma)}(\widetilde{W})}{2\pi\delta}+\frac{1}{48\pi}\frac{1}{48\pi}\int_{\widetilde{W}}\|\nabla(\theta(\rho_{j}+\sigma))\|^{2}dz,\qquad j=1,2. (4.8)

By Proposition 4.7 with f=θρ~jf=\theta\widetilde{\rho}_{j} and h=σh=\sigma, together with the fact that ρ~1ρ~2\widetilde{\rho}_{1}\equiv\widetilde{\rho}_{2} outside UU and θ1\theta\equiv 1 in VV, we can write the difference in the quantities (4.8) as the difference in the quantities

Vol~ρ~j+σ(W~)2πδ+148πW~(σρ~j2+2Kσρ~j)Vol~σ(dz),j=1,2,\frac{\widetilde{\mathrm{Vol}}_{\widetilde{\rho}_{j}+\sigma}(\widetilde{W})}{2\pi\delta}+\frac{1}{48\pi}\int_{\widetilde{W}}(\|\nabla_{\sigma}\widetilde{\rho}_{j}\|^{2}+2K_{\sigma}\widetilde{\rho}_{j})\widetilde{\mathrm{Vol}}_{\sigma}(dz),\quad j=1,2, (4.9)

where σ,Kσ\nabla_{\sigma},K_{\sigma} are the gradient and Gauss curvature associated to (W~,eσ|dz|2)(\widetilde{W},e^{\sigma}|dz|^{2}). Pulling back via φ\varphi, we can rewrite the expressions (4.9) as

Volρj(W)2πδ+148πW(ρj2+2Kρj)Vol(dz),j=1,2.\frac{\mathrm{Vol}_{\rho_{j}}(W)}{2\pi\delta}+\frac{1}{48\pi}\int_{W}(\|\nabla\rho_{j}\|^{2}+2K\rho_{j})\mathrm{Vol}(dz),\quad j=1,2.

We complete the proof by noting that, since ρ1ρ2\rho_{1}\equiv\rho_{2} outside WW, we can replace WW by MM in (4.2). ∎

Proof of Theorem 1.9.

By conformal invariance of the Brownian loop measure [APPS20, Lemma 3.3], the statement of the theorem is equivalent to the assertion that (1.7) holds up to scaling. The result holds for ρ0\rho\equiv 0 by [APPS20, Theorem 1.3].666We remark that the δ\delta in the statement of [APPS20, Theorem 1.3] represents the quadratic variation of Brownian loops, which is two times the time interval length we use in this paper. Thus, for our application, we need to substitute δ\delta there into δ/2\delta/2. We deduce the result for general ρ\rho\in\mathcal{B} by applying Propositions 4.4 and 4.5. ∎

5 Constructing square subdivision regularizations

Now what happens if ρ\rho is defined from a finite square subdivision as in [APPS20, Section 6], so that it has constant Laplacians on each square in a grid? These functions can be shown to be C1C^{1}, but along the edges of the squares they fail to be C2C^{2}. In this section, we explain enough about these function to make it clear that they fit into the framework of this paper, at least if one restricts attentions to those for which the square averages are restricted to a compact set.

One can construct functions with piecewise-constant Laplacian on squares somewhat explicitly. Consider the function on f(z)=Re[z2log(z2)]/π+|z|2/2f(z)=\mathrm{Re}[z^{2}\log(z^{2})]/\pi+|z|^{2}/2 restricted to the quadrant Q={z:arg(z)(π/4,π/4)}Q=\{z:\arg(z)\in(-\pi/4,\pi/4)\}. Note that on the boundary of the quadrant, z2z^{2} is purely imaginary, equal to |z2|i|z^{2}|i on the upper boundary ray and |z2|i-|z^{2}|i on the lower, while Imlog(z2)=arg(z2)\mathrm{Im}\log(z^{2})=\arg(z^{2}) is π/2\pi/2 on the upper boundary and π/2-\pi/2 on the lower. so that Re[z2log(z2)]/π=|z2|/2\mathrm{Re}[z^{2}\log(z^{2})]/\pi=-|z^{2}|/2 on both rays. Thus ff has constant Laplacian and equals zero on the boundary QQ. We can extend the definition of ff to the other three quadrants (iQiQ, iQ-iQ, and Q-Q) by imposing the relation f(iz)=f(z)f(iz)=-f(z). In a small neighborhood of the origin, the ff defined this way is negative on ±Q\pm Q and positive on ±iQ\pm iQ. The complex derivatives of g(z)=z2log(z2)=2z2log(z)g(z)=z^{2}\log(z^{2})=2z^{2}\log(z) are first g(z)=4zlog(z)+2zg^{\prime}(z)=4z\log(z)+2z, second g′′(z)=4log(z)+6g^{\prime\prime}(z)=4\log(z)+6, third g′′′(z)=4/zg^{\prime\prime\prime}(z)=4/z, etc. One can deduce from this that ff is differentiable as a real-valued function (with derivative 0 at 0) but that its second derivatives blow up slowly (logarithmically) near zero and also have discontinuities along the quadrant boundaries.

The function f(ωz)f(\omega z) (where ω\omega is a fixed eighth root of unity) is thus a C1C^{1} function that has piecewise constant Laplacian on the four standard quadrants, while being equal to zero on the boundaries of these quadrants.

By taking linear combinations of f(ωz)f(\omega z) and the four functions [max(Re(az),0)]2\bigl{[}\max\bigl{(}\mathrm{Re}(az),0\bigr{)}\bigr{]}^{2} (with a{±1,±i}a\in\{\pm 1,\pm i\}) we can get a differentiable function whose Laplacian matches any function that is constant on each of the four standard quadrants — in particular a function whose Laplacian is 11 on one quadrant and 0 on the other three. By taking differences of translates of this function, we can obtain a function whose Laplacian is 11 on a semi-strip or a rectangle (and 0 elsewhere). Linear combinations of these allow us to describe a function ϕ\phi whose Laplacian is any given function that is constant on the squares of the grid. Any other function with the same Laplacian has to then differ from ϕ\phi by a harmonic function. For example, by subtracting the harmonic extension of the values of ϕ\phi on D\partial D one obtains a function with the desired Laplacian whose boundary values on D\partial D are zero.

6 Open questions

In this section, we discuss some open questions. First, as already mentioned in Footnote 3, we expect that some of our main results can be extended to ρW1,2(D)\rho\in W^{1,2}(D). In addition to the fact that the Dirichlet energy is well-defined for such functions, here is another reason that the arguments from the Lipschitz case might carry over to this setting. Suppose that MρM-M\leq\rho\leq M. Then the function

f(x):={ex1 for x[M,M]0 otherwise.f(x):=\begin{cases}e^{x}-1&\text{ for }x\in[-M,M]\\ 0&\text{ otherwise.}\end{cases}

is a Lipschitz function on [M,M][-M,M] with f(0)=0f(0)=0. Therefore, the celebrated Stampacchia’s theorem [EG18, §4.2.2] asserts that ρW1,2(D)\rho\in W^{1,2}(D) implies eρ1e^{\rho}-1 is W1,2(D)W^{1,2}(D) and the chain rule of weak derivatives implies that

(eρ(x))=eρ(x)ρ(x)\nabla(e^{\rho(x)})=e^{\rho(x)}\nabla\rho(x)

holds for almost every xx. This already implies that many parts of our proofs carry over to the general case. As an example, Lemma 1.17 follows almost immediately with the condition of precompactness in W1,2(D)W^{1,2}(D), or equivalently the uniform equicontinuity in W1,2(D)W^{1,2}(D), in the sense of Remark 1.7. However, there are a few technical issues where we need to improve our estimates. For example, our proof of Lemma 1.16 is based on a pointwise estimate of ρ\nabla\rho, which would have to be replaced by an L2L^{2} estimate, and such an estimate does not immediately follow from our current arguments.

Question 6.1.

Prove that Theorem 1.4 holds for a larger class of functions. For example, one might consider functions in W1,2(D)W^{1,2}(D) and take \mathcal{B} to be any precompact subset of W1,2(D)W^{1,2}(D), possibly with some additional conditions.

Next, we note that Lemma 1.17 is proven for Schwartz functions θ\theta, which is not expected to be optimal in any sense other than making the proofs simple. In fact, we believe the result holds for a much more general class of functions. For instance, we apply the lemma to expected occupation measure, but it should also be true for a single occupation measure in some sense.

Question 6.2.

Prove that Lemma 1.17 holds for a more general class of functions. For example, prove a similar result for a random function describing the occupation measure of a Brownian loop.

Third, we may consider the case when the manifold MM has a boundary, say when M=M=\mathbb{H}. In principle, one could formulate the boundary problem in a style similar to that presented above, and try to weaken the required regularity along the boundary as well. To start, imagine the boundary line is the horizontal real axis. Let dbd\ell_{b} be the measure on unit length loops that hit the real axis and are centered at a point on the positive imaginary axis.

This measure can be obtained from by starting with dd\ell, then weighting by the gap between the minimum and maximum imaginary values obtained by the loop, then shifting the vertical height loop to a uniformly random height (within this range). The expected maximal height of a Brownian bridge is π/8\sqrt{\pi/8} (as can be proved using the reflection principle). The expected minimal height is thus minus that, and the expected difference 2π/82\sqrt{\pi/8}, which means that the expected vertical gap between the maximum height and the central height is again π/8\sqrt{\pi/8}.

So formally, instead of being a probability measure, dbd\ell_{b} is a measure with weight π/8\sqrt{\pi/8}. For real values xx, write bclen(x,t,b)=eρ(x)δ\mathrm{bclen}(x,t,\ell_{b})=e^{\rho(x)}\delta. Note that the set of loops of length greater than δ\delta is given by

δπ/8t12πt2𝑑t=2π/82πδ=122πδ\int_{\delta}^{\infty}\sqrt{\pi/8}\sqrt{t}\frac{1}{2\pi t^{2}}dt=\frac{2\sqrt{\pi/8}}{2\pi\sqrt{\delta}}=\frac{1}{2\sqrt{2\pi\delta}}

We can imagine ρ\rho is defined in a neighborhood of a real line segment, and then try to measure the mass of loops (of size greater than δ\delta) that hit the boundary line itself. The relevant quantities are the gradient in the parallel direction and the gradient in the normal direction (with the latter affecting whether the real axis is hit by more loops centered above or below).

Question 6.3.

Prove that Theorem 1.9 also holds for manifolds with boundaries. In other words, generalize the boundary case of [APPS20, Theorem 1.3]:777Similar to the case without boundary, we plug δ/2\delta/2 instead of δ\delta in the statement of [APPS20, Theorem 1.3], where δ\delta there represents the quadratic variation of Brownian loops, which is two times the time interval length we use in this paper.

Conjecture 6.4.

Let (M,g)(M,g) be a fixed compact smooth two-dimensional Riemannian manifold with smooth boundary, and we let μloop\mu^{\text{loop}} denote the Brownian loop measure on (M,g)(M,g). Let KK be the Gaussian curvature on MM, let Δ\Delta be the Laplacian associated to (M,g)(M,g), and let detζΔ\det_{\zeta}^{\prime}\Delta denote its zeta-regularized determinant. Let \mathcal{B} be a family of Lipschitz functions that (1) has uniformly bounded Lipshitz constants, and (2) is precompact in W1,1(M)W^{1,1}(M).

The μloop\mu^{\text{loop}}-mass of loops with ρ\rho-length greater than δ\delta is given by

Volρ(M)2πδLenρ(M)22πδχ(M)6(logδ2+γ)+148πM(ρ2+2Kρ)Vol(dz)\displaystyle\frac{\mathrm{Vol}_{\rho}(M)}{2\pi\delta}-\frac{\mathrm{Len}_{\rho}(\partial M)}{2\sqrt{2\pi\delta}}-{\frac{\chi(M)}{6}}(\log\frac{\delta}{2}+\upgamma)+\frac{1}{48\pi}\int_{M}(\|\nabla\rho\|^{2}+2K\rho)\,\mathrm{Vol}(dz)
+logVol(M)logVolρ(M)logdetζΔ+O(δ1/2),\displaystyle\qquad\qquad\qquad+\log\mathrm{Vol}(M)-\log\mathrm{Vol}_{\rho}(M)-\log\det\nolimits_{\zeta}^{\prime}\Delta+O(\delta^{1/2}),

with the convergence as δ0\delta\to 0 uniform over ρ\rho\in\mathcal{B}, where γ0.5772\upgamma\approx 0.5772 is the Euler-Mascheroni constant.

Finally, one naive approach to generalizing [APPS20, Proposition 6.9] for lower-regularity settings is to follow the zeta-regularization procedure verbatim, hoping each step works for the generalized settings in a similar way. The first obstacle in this direction is the lack of short-time expansions for the trace of heat kernels. For our application, what we need is if (M,g)(M,g) is a two-dimensional non-smooth manifold without boundary, and Δ\Delta is the associated Laplacian defined in terms of Brownian loop mass, then

tr(eδΔ/2)=Volg(M)/δ+χ(M)/6+o(1).\mathrm{tr}(e^{-\delta\Delta/2})=\mathrm{Vol}_{g}(M)/\delta+\chi(M)/6+o(1).
Question 6.5.

Prove the short time expansion for the (trace) heat kernel holds for lower regularity metrics.

References

  • [AAR13] P. Albin, C. L. Aldana, and F. Rochon. Ricci flow and the determinant of the Laplacian on non-compact surfaces. Communications in Partial Differential Equations, 38(4):711–749, 2013.
  • [AHT18] L. Ambrosio, S. Honda, and D. Tewodrose. Short-time behavior of the heat kernel and weyl’s law on RCD(k,n)\text{RCD}^{*}(k,n) spaces. Annals of Global Analysis and Geometry, 53(1):97–119, 2018.
  • [Alv83] O. Alvarez. Theory of strings with boundaries: Fluctuations, topology and quantum geometry. Nuclear Physics B, 216(1):125 – 184, 1983.
  • [APPS20] M. Ang, M. Park, J. Pfeffer, and S. Sheffield. Brownian loops and the central charge of a Liouville random surface. arXiv preprint arXiv:2005.11845, 2020.
  • [AR18] C. L. Aldana and J. Rowlett. A Polyakov formula for sectors. The Journal of Geometric Analysis, 28(2):1773–1839, 2018.
  • [BB11] H. Brezis and H. Brézis. Functional analysis, Sobolev spaces and partial differential equations, volume 2. Springer, 2011.
  • [BGV03] N. Berline, E. Getzler, and M. Vergne. Heat kernels and Dirac operators. Springer Science & Business Media, 2003.
  • [Din02] Y. Ding. Heat kernels and green’s functions on limit spaces. Communications in Analysis and Geometry, 10(3):475–514, 2002.
  • [Don01] H. Donnelly. Bounds for eigenfunctions of the laplacian on compact riemannian manifolds. Journal of Functional Analysis, 187(1):247–261, 2001.
  • [Dub09] J. Dubédat. SLE and the free field: partition functions and couplings. J. Amer. Math. Soc., 22(4):995–1054, 2009, 0712.3018. MR2525778 (2011d:60242)
  • [EG18] L. C. Evans and R. F. Garzepy. Measure theory and fine properties of functions. Routledge, 2018.
  • [Gig18] N. Gigli. Lecture notes on differential calculus on rcd spaces. Publications of the Research Institute for Mathematical Sciences, 54(4):855–918, 2018.
  • [Gil18] P. B. Gilkey. Invariance theory: the heat equation and the Atiyah-Singer index theorem, volume 16. CRC press, 2018.
  • [Gre21] R. L. Greenblatt. Discrete and zeta-regularized determinants of the Laplacian on polygonal domains with dirichlet boundary conditions. arXiv preprint arXiv:2102.04837, 2021.
  • [Gri06] A. Grigor’yan. Heat kernels on weighted manifolds and applications. Cont. Math, 398(2006):93–191, 2006.
  • [Hör68] L. Hörmander. The spectral function of an elliptic operator. Acta mathematica, 121(1):193–218, 1968.
  • [Kac66] M. Kac. Can one hear the shape of a drum? The american mathematical monthly, 73(4P2):1–23, 1966.
  • [Kal21] V. Kalvin. Polyakov-alvarez type comparison formulas for determinants of laplacians on riemann surfaces with conical singularities. Journal of Functional Analysis, 280(7):108866, 2021.
  • [Kok13] A. Kokotov. Polyhedral surfaces and determinant of Laplacian. Proceedings of the American Mathematical Society, 141(2):725–735, 2013.
  • [LV09] J. Lott and C. Villani. Ricci curvature for metric-measure spaces via optimal transport. Annals of Mathematics, pages 903–991, 2009.
  • [LW04] G. F. Lawler and W. Werner. The Brownian loop soup. Probability theory and related fields, 128(4):565–588, 2004.
  • [Moo99] E. A. Mooers. Heat kernel asymptotics on manifolds with conic singularities. Journal d’Analyse Mathematique, 78(1):1–36, 1999.
  • [MS+67] H. P. McKean, I. M. Singer, et al. Curvature and the eigenvalues of the Laplacian. J. Differential Geometry, 1(1):43–69, 1967.
  • [OPS88] B. Osgood, R. Phillips, and P. Sarnak. Extremals of determinants of Laplacians. Journal of Functional Analysis, 80(1):148 – 211, 1988.
  • [Pet10] A. Petrunin. Alexandrov meets lott–villani–sturm. arXiv preprint arXiv:1003.5948, 2010.
  • [Pol81] A. M. Polyakov. Quantum geometry of bosonic strings. Phys. Lett. B, 103(3):207–210, 1981. MR623209 (84h:81093a)
  • [RS73] D. Ray and I. Singer. Analytic torsion for complex manifolds. Annals Math., 98:154–177, 1973.
  • [Sar87] P. Sarnak. Determinants of Laplacians. Communications in mathematical physics, 110(1):113–120, 1987.
  • [Ser61] J. Serrin. On the differentiability of functions of several variables. Archive for Rational Mechanics and Analysis, 7(1):359–372, 1961.
  • [She07] S. Sheffield. Gaussian free fields for mathematicians. Probab. Theory Related Fields, 139(3-4):521–541, 2007, math/0312099. MR2322706 (2008d:60120)
  • [She13] D. Sher. The heat kernel on an asymptotically conic manifold. Analysis & PDE, 6(7):1755–1791, 2013.
  • [She15] D. A. Sher. Conic degeneration and the determinant of the Laplacian. Journal d’Analyse Mathématique, 126(1):175–226, 2015.
  • [Stu06] K.-T. Sturm. On the geometry of metric measure spaces. Acta mathematica, 196(1):65–131, 2006.
  • [Suz19] K. Suzuki. Convergence of brownian motions on metric measure spaces under riemannian curvature–dimension conditions. Electronic Journal of Probability, 24:1–36, 2019.
  • [Wan18] Y. Wang. Equivalent descriptions of the Loewner energy. Inventiones mathematicae, pages 1–49, 2018.
  • [Wey11] H. Weyl. Über die asymptotische verteilung der eigenwerte. Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1911:110–117, 1911.
  • [ZZ17] H.-C. Zhang and X.-P. Zhu. Weyl’s law on RCD(k,n)\text{RCD}^{*}(k,n) metric measure spaces. arXiv preprint arXiv:1701.01967, 2017.