Critical two-point function for long-range $O(n)$ models below the upper critical dimension

TL;DR

This paper rigorously proves that in long-range $O(n)$ models below the upper critical dimension, the critical two-point function decays with a mean-field exponent, confirming a prediction from physics literature using renormalisation group and cluster expansion techniques.

Contribution

It provides the first rigorous proof of the critical two-point function decay with mean-field exponent in long-range models below the upper critical dimension.

Findings

Critical two-point function decays as $r^{-(d-rac{1}{2}(d+ ext{epsilon}))}$

Validation of physics prediction from 1972 about critical exponents

Application of rigorous renormalisation group and cluster expansion methods

Abstract

We consider the $n$-component $|\varphi|^4$ lattice spin model ($n \ge 1$) and the weakly self-avoiding walk ($n=0$) on $\mathbb{Z}^d$, in dimensions $d=1,2,3$. We study long-range models based on the fractional Laplacian, with spin-spin interactions or walk step probabilities decaying with distance $r$ as $r^{-(d+\alpha)}$ with $\alpha \in (0,2)$. The upper critical dimension is $d_c=2\alpha$. For $\epsilon >0$, and $\alpha = \frac 12 (d+\epsilon)$, the dimension $d=d_c-\epsilon$ is below the upper critical dimension. For small $\epsilon$, weak coupling, and all integers $n \ge 0$, we prove that the two-point function at the critical point decays with distance as $r^{-(d-\alpha)}$. This "sticking" of the critical exponent at its mean-field value was first predicted in the physics literature in 1972. Our proof is based on a rigorous renormalisation group method. The treatment of observables differs from that used in recent work on the nearest-neighbour 4-dimensional case, via our use of a cluster expansion.

Equations

$$\lambda (k)=4\sum _{j=1}^d{\sin }^2(k_j/2)=2\sum _{j=1}^d(1-\cos k_j).$$

$$-{\Delta }_{x,y}=\frac{1}{(2\pi {)}^d}{\int }_{[-\pi ,\pi {]}^d}\lambda (k)e^{ik\cdot (x-y)}dk.$$

$$(-\Delta {)}_{x,y}^{\alpha /2}=\frac{1}{(2\pi {)}^d}{\int }_{[-\pi ,\pi {]}^d}\lambda (k{)}^{\alpha /2}{e}^{ik\cdot (x-y)}dk.$$

$$-(-\Delta {)}_{x,y}^{\alpha /2}\asymp |x-y{|}^{-(d+\alpha )}$$

$$((-\Delta {)}^{\alpha /2}+m^2{)}_{0,x}^{-1}\le c\frac{1}{|x{|}^{d-\alpha }}\frac{1}{1+m^4|x{|}^{2\alpha }},$$

$$((-\Delta {)}^{\alpha /2}{)}_{0,x}^{-1}\sim c_{d,\alpha }\frac{1}{|x{|}^{d-\alpha }}$$

$$(-{\Delta }_{{\Lambda }_N}{)}_{x,y}^{\alpha /2}=\sum _{z\in {\mathbb{Z}}^d}(-\Delta {)}_{x,y+z{L}^N}^{\alpha /2}(x,y\in {\Lambda }_N).$$

$$V(\varphi)=\sum_{x\in\Lambda}\big{(}\tfrac{1}{4}g|\varphi_{x}|^{4}+\textstyle{\frac{1}{2}}\nu|\varphi_{x}|^{2}+\textstyle{\frac{1}{2}}\varphi_{x}\cdot((-\Delta_{\Lambda})^{\alpha/2}\varphi)_{x}\big{)}.$$

$$Z_{g,\nu ,N}={\int }_{({\mathbb{R}}^n{)}^{\Lambda }}{e}^{-V(\phi )}d\phi ,$$

$$⟨F{⟩}_{g,\nu ,N}=\frac{1}{Z_{g,\nu ,N}}{\int }_{({\mathbb{R}}^n{)}^{\Lambda }}F(\phi )e^{-V(\phi )}d\phi .$$

$$G_{a,b,N}(g,\nu ;n)$$

$$G_{a,b}(g,\nu ;n)$$

$$\chi (g,\nu ;n)=\underset{N\to \infty }{\lim }\sum _{b\in {\Lambda }_N}{G}_{a,b,N}(g,\nu ;n)$$

$$P_a(X(t)=b)=E_a(\text{mathbbm}{1}_{X(t)=b})=(e^{tQ}{)}_{a,b}.$$

$$I_{T}=\sum_{x\in{\mathbb{Z}}^{d}}\big{(}L_{T}^{x}\big{)}^{2}=\int_{0}^{T}\!\!\int_{0}^{T}\mathbbm{1}_{X(t_{1})=X(t_{2})}\;dt_{1}\,dt_{2}.$$

$$G_{a,b}(g,\nu ;0)={\int }_0^{\infty }{E}_a\left(e^{-g{I}_T}\text{mathbbm}{1}_{X(T)=b}\right)e^{-\nu T}dT,$$

$$\chi (g,\nu ;0)=\sum _{b\in {\mathbb{Z}}^d}{G}_{a,b}(g,\nu ;0)={\int }_0^{\infty }{E}_a(e^{-g{I}_T})e^{-\nu T}dT.$$

$$C^{-1}{t}^{-(1+\frac{n+2}{n+8}\frac{ϵ}{\alpha }-C{ϵ}^2)}\le \chi (g,\nu ;n)\le C{t}^{-(1+\frac{n+2}{n+8}\frac{ϵ}{\alpha }+C{ϵ}^2)}.$$

$$\gamma =1+\frac{n+2}{n+8}\frac{ϵ}{\alpha }+O({ϵ}^2)(n\ge 0).$$

$$G_{a,b}(g,{\nu }_c;n)=(1+O(ϵ))((-\Delta {)}^{\alpha /2}{)}_{a,b}^{-1}\asymp \frac{1}{|a-b{|}^{d-\alpha }}.$$

$$U_0$$

$$G_{a,b,N}(g,{\nu }_0^c(m^2)+m^2;n)$$

$${\mathbb{E}}_{C_j}\theta \circ \cdots \circ {\mathbb{E}}_{C_1}\theta e^{-U_0}\approx e^{-U_j},$$

$$U_j$$

$$G_{a,b,N}(g,{\nu }_0^c+m^2;n)$$

$$(g_j,{\nu }_j,u_j,{\lambda }_{x,j},q_{x,j})\to (g_{j+1},{\nu }_{j+1},u_{j+1},{\lambda }_{x,j+1},q_{x,j+1}),$$

$$q_{x,j}\approx {\lambda }_{a,j_{ab}}{\lambda }_{b,j_{ab}}{w}_{j;a,b},$$

$${\lambda }_{a,j}=e^{u_j|\Lambda |}\bar{D}{D}_{{\sigma }_a}{e}^{-U_j}\approx e^{u_j|\Lambda |}\bar{D}{D}_{{\sigma }_a}{\mathbb{E}}_{w_j}\theta e^{-U_0}.$$

$$\bar{D}{D}_{{\sigma }_a}{\mathbb{E}}_{w_j}\theta e^{-U_0}$$

$${\lambda }_{a,j}\approx 1+w_j^{(1)}{\nu }_j.$$

$${\tau }_x=\frac{1}{2}|{\phi }_x{|}^2(n\ge 1).$$

Full text

1 Introduction and main result

Broadly speaking, the mathematical understanding of critical phenomena for spin systems has progressed in dimension $d=2$, where exact solutions and SLE are important tools; in dimensions $d>4$, where infrared bounds and the lace expansion are useful; and in dimension $d=4$, where renormalisation group (RG) methods have been applied. The physically most important case of $d=3$ is more difficult, and mathematical methods are scarce.

In the physics literature, the $\epsilon$-expansion was introduced to study non-integer dimensions slightly below $d=4$. An alternate approach is to consider long-range models, which change the upper critical dimension from $d_{c}=4$ to a lower value $d_{c}=2\alpha$ with $\alpha\in(0,2)$. By choosing $d=1,2,3$ and $\alpha=\frac{1}{2}(d+\epsilon)$ with small $\epsilon$, it is possible to study integer dimension $d$ which is slightly below the upper critical dimension $2\alpha=d+\epsilon$. In this paper, we consider $n$-component spins and the weakly self-avoiding walk in this long-range context, and prove that the critical two-point function has mean-field decay $r^{-(d-\alpha)}$ also below the upper critical dimension. Our method involves a RG analysis in the vicinity of a non-Gaussian fixed point.

1.1 Introduction

We consider long-range $O(n)$ models on ${{{\mathbb{Z}}}^{d}}$ for integers $n\geq 0$ and dimensions $d=1,2,3$. The case $n=0$ is the continuous-time weakly self-avoiding walk, and the case $n\geq 1$ is the $n$-component $|\varphi|^{4}$ lattice spin model. For $n=0$ the underlying random walk model takes steps of length $r$ with probabilities decaying as $r^{-(d+\alpha)}$ with $\alpha\in(0,2)$, and for $n\geq 1$ the spin-spin interaction in the Hamiltonian has that same decay. More precisely, the models are based on the fractional Laplacian $(-\Delta)^{\alpha/2}$, whose kernel decays at large distance as $r^{-(d+\alpha)}$.

The upper critical dimension is predicted to be $d_{c}=2\alpha$ for all $n\geq 0$. Thus, for $\alpha<\frac{d}{2}$, mean-field behaviour is predicted; this has been proved for self-avoiding walk, for the Ising model, for the 1-component $\varphi^{4}$ model, and for other models [3, 17, 22, 23]. In the physics literature, it is observed that below the upper critical dimension the critical two-point function continues to exhibit the mean-field decay $r^{-(d-\alpha)}$ for $\alpha\in(\frac{d}{2},2-\eta)$, and then crosses over to $r^{-(d-2+\eta)}$ decay for $\alpha\in(2-\eta,2)$. Here $\eta$ is the exponent for the nearest-neighbour model; for $n=1$ this is $\eta=\frac{1}{4}$ for $d=2$ [35], and a recent estimate for $d=3$ is $\eta=0.03631(3)$ [18]. The earliest paper to elucidate the critical behaviour of long-range models is [20], with [33] roughly contemporaneous and [29] providing further development. A very recent paper which analyses the crossover for the two-point function in detail for $n=1$ is [8]. At the crossover, when $\alpha=\alpha_{*}=2-\eta$, a logarithmic correction is predicted, with overall decay $\frac{1}{r^{d-\alpha_{*}}}\frac{1}{\log r}$ [11, 8]. The relationship with conformal invariance is explained in [28].

Let $n=0,1,2,\ldots$; $d=1,2,3$; and $\alpha=\frac{1}{2}(d+\epsilon)$. We use a rigorous RG argument to prove that for small $\epsilon>0$, the critical two-point function has decay $r^{-(d-\alpha)}$. This proves the “sticking” of the critical exponent at its mean-field value, for $\alpha$ slightly above $\frac{d}{2}$, or equivalently, for $d$ slightly below the upper critical dimension $d_{c}=2\alpha$. Our proof extends recent results and methods used to study the $\epsilon$-expansion for the critical exponents for the susceptibility and specific heat of the long-range models [31]. It also relies on results and techniques developed to study related problems for the 4-dimensional nearest-neighbour models [16, 5, 32]. However, our treatment of observables differs from that used in the 4-dimensional case, via our application of a cluster expansion.

Earlier mathematical work which applies RG methods to long-range models includes the construction of global RG trajectories for $n=0$ and $d=3$ [27], and for a continuum version of the $n=1$ model in [12, 1]. These references do not study critical exponents. The exponents for critical correlations in a certain hierarchical version of the model, for $d=3$ and $n=1$, are computed in [2]. For a closely related continuum model with $n=1$ in dimensions $d=2,3$, a proof of the “sticking” of the critical exponent for the critical two-point function was announced in a 2013 lecture [24].

1.2 Fractional Laplacian

The models we study are defined in terms of the fractional Laplacian. We now define the fractional Laplacian and list some of its properties. Further details can be found in [31, Sections LABEL:alpha-sec:fL–LABEL:alpha-sec:frd].

Let $d\geq 1$ and $\alpha\in(0,2)$. We write $|x|$ for the Euclidean norm of $x\in{\mathbb{Z}}^{d}$. Let $J$ be the ${{{\mathbb{Z}}}^{d}}\times{{{\mathbb{Z}}}^{d}}$ matrix with $J_{xy}=1$ if $|x-y|=1$, and otherwise $J_{xy}=0$. Let $I$ denote the identity matrix. The lattice Laplacian on ${{{\mathbb{Z}}}^{d}}$ is $\Delta=J-2dI$. For $k=(k_{1},\ldots,k_{d})\in[-\pi,\pi]^{d}$, let

[TABLE]

The matrix element $-\Delta_{x,y}$ can be written as the Fourier integral

[TABLE]

The fractional Laplacian is the matrix $(-\Delta)^{\alpha/2}$ defined by

[TABLE]

For $|x-y|\rightarrow\infty$, the fractional Laplacian decays as

[TABLE]

(see [31, Lemma LABEL:alpha-lem:fracLapdecay], or [10, Theorem 5.3] for a more precise and more general statement). Here, and in the following, we write $a\asymp b$ to denote the existence of $c>0$ such that $c^{-1}b\leq a\leq cb$. For $d\geq 1$, $\alpha\in(0,2\wedge d)$, $\bar{m}^{2}>0$, $m^{2}\in[0,\bar{m}^{2}]$, and $x\neq 0$, the resolvent obeys

[TABLE]

with $c$ depending on $d,\alpha,\bar{m}^{2}$ (see [31, Lemma LABEL:alpha-lem:covbd]). For $m^{2}=0$, an asymptotic formula

[TABLE]

is proven in [9, Theorem 2.4], with precise constant $c_{d,\alpha}$.

Given integers $L,N>1$, let $\Lambda=\Lambda_{N}={\mathbb{Z}}^{d}/L^{N}{\mathbb{Z}}^{d}$ denote the $d$-dimensional discrete torus of side length $L^{N}$. The torus fractional Laplacian is defined by

[TABLE]

The sum on the right-hand side of (1.7) converges, by (1.4).

1.3 The $|\varphi|^{4}$ model

We first define the model on the torus $\Lambda=\Lambda_{N}$, as usual for spin systems. Let $d\geq 1$ and $\alpha\in(0,2)$. Let $n\geq 1$. The spin field $\varphi$ is a function $\varphi:\Lambda\rightarrow{\mathbb{R}}^{n}$, denoted $x\mapsto\varphi_{x}$, which we may regard as an element $\varphi\in({\mathbb{R}}^{n})^{\Lambda}$. The Euclidean norm of $v=(v^{1},\ldots,v^{n})\in{\mathbb{R}}^{n}$ is $|v|=[\sum_{i=1}^{n}(v^{i})^{2}]^{1/2}$, with inner product $v\cdot w=\sum_{i=1}^{n}v^{i}w^{i}$. We extend the action of the fractional Laplacian to act on the spin field component-wise, namely $((-\Delta_{\Lambda})^{\alpha/2}\varphi)_{x}^{i}=\sum_{y\in\Lambda}(-\Delta_{\Lambda})^{\alpha/2}_{x,y}\varphi^{i}_{y}$.

Given $g>0$ and $\nu\in{\mathbb{R}}$, we define the interaction $V:({\mathbb{R}}^{n})^{\Lambda}\rightarrow{\mathbb{R}}$ by

[TABLE]

The partition function is defined by

[TABLE]

where $d\varphi$ is the Lebesgue measure on $({\mathbb{R}}^{n})^{\Lambda}$. The expectation of a random variable $F:({\mathbb{R}}^{n})^{\Lambda}\rightarrow{\mathbb{R}}$ is

[TABLE]

Given lattice points $a,b$, we define the finite- and infinite-volume two-point function by

[TABLE]

On the left-hand side of (1.12) we have $a,b\in{{{\mathbb{Z}}}^{d}}$, and on the right-hand side we identify these points with elements of $\Lambda_{N}$ for large $N$, by regarding the vertices of $\Lambda_{N}$ as a cube in ${\mathbb{Z}}^{d}$ (without boundaries identified) approximately centred at the origin. The susceptibility is defined by

[TABLE]

and can be used to identify the critical point of the model. By translation invariance, $\chi$ is independent of $a$. Existence of the infinite volume limits in (1.12)–(1.13), in our context, is discussed below.

1.4 Weakly self-avoiding walk

Let $d\geq 1$ and $\alpha\in(0,2)$. Let $X$ denote the continuous-time Markov chain with state space ${{{\mathbb{Z}}}^{d}}$ and infinitesimal generator $Q=-(-\Delta_{{{{\mathbb{Z}}}^{d}}})^{\alpha/2}$. Verification that $Q$ has the attributes required of a generator is given in [31, Lemma LABEL:alpha-lem:DJ]. Let $P$ be the probability measure associated with $X$, and $E$ the corresponding expectation; a subscript $a$ specifies $X(0)=a$. The transition probabilities are given by

[TABLE]

The local time of $X$ at $x$ up to time $T$ is the random variable $L_{T}^{x}=\int_{0}^{T}\mathbbm{1}_{X(t)=x}\;dt$. The self-intersection local time up to time $T$ is the random variable

[TABLE]

Given $g>0$, $\nu\in{\mathbb{R}}$, and $a,b\in{{{\mathbb{Z}}}^{d}}$, the continuous-time weakly self-avoiding walk two-point function is defined by the integral

[TABLE]

and the susceptibility is defined by

[TABLE]

The labels [math] on the left-hand sides of (1.16)–(1.17) reflect the fact that the weakly self-avoiding walk corresponds to the formal $n=0$ case of the $n$-component $|\varphi|^{4}$ model. As in earlier work on the 4-dimensional case, e.g., [32, 31], we treat both cases $n\geq 1$ (spins) and $n=0$ (self-avoiding walk) simultaneously and rigorously, via a supersymmetric spin representation for the weakly self-avoiding walk.

1.5 Susceptibility and critical point

Let $d=1,2,3$; $n\geq 0$; $L$ be sufficiently large; $\epsilon>0$ be sufficiently small; and $\alpha=\frac{1}{2}(d+\epsilon)$. Let $\tau^{(\alpha)}$ denote the diagonal element of the Green function, i.e., $\tau^{(\alpha)}=((-\Delta)^{\alpha/2})^{-1}_{00}$. One of the main results of [31] is that there exists $\bar{s}\asymp\epsilon$ such that, for $g\in[\frac{63}{64}\bar{s},\frac{65}{64}\bar{s}]$, there exist $\nu_{c}=\nu_{c}(g;n)=-(n+2)\tau^{(\alpha)}g(1+O(g))$ and $C>0$ such that for $\nu=\nu_{c}+t$ with $t\downarrow 0$,

[TABLE]

This is a statement that there is a critical point at $\nu=\nu_{c}$, and that the critical exponent $\gamma$ exists to order $\epsilon$, with

[TABLE]

It is part of the statement that for $n\geq 1$ the susceptibility is given by the infinite-volume limit (1.13), under the above hypotheses. The critical exponent for the specific heat is also computed to order $\epsilon$ in [31], for $n\geq 1$.

1.6 Main result

Our main result is the following theorem, which shows that just below the upper critical dimension, the exponent for the critical two-point function “sticks” at its mean-field value (see (1.6)), as predicted by [20]. The theorem applies for all $n\geq 0$, including the case $n=0$ of the weakly self-avoiding walk. The critical value $\nu_{c}=\nu_{c}(g;n)$ is the one mentioned in Section 1.5. As part of the proof of the theorem, it is shown that for $n\geq 1$ the infinite-volume limit (1.12) exists for $\nu=\nu_{c}$.

Theorem 1.1**.**

Let $d=1,2,3$; $n\geq 0$; $L$ be sufficiently large; $\epsilon>0$ be sufficiently small; and $\alpha=\frac{1}{2}(d+\epsilon)$. For $g\in[\frac{63}{64}\bar{s},\frac{65}{64}\bar{s}]$ the critical two-point function obeys, as $|a-b|\rightarrow\infty$,

[TABLE]

Note that Theorem 1.1 identifies the constant in the decay of the interacting two-point function only up to an error of order $\epsilon$. However, the error is uniformly bounded in $a,b$, so the power in the decay rate takes its mean-field value, and this is true to all orders in $\epsilon$.

1.7 Strategy of proof

The proof is based on a rigorous RG method developed in a series of papers by Bauerschmidt, Brydges and Slade, where the focus is on the nearest-neighbour models in dimension 4. The method is adapted to the long-range setting in [31].

Fix $g$ as in the statement of Theorem 1.1. In [31], given small $m^{2}>0$, a critical value $\nu_{0}^{c}(m^{2})$ is constructed, with the property that the critical point $\nu_{c}$ is given by $\nu_{c}=\lim_{m^{2}\downarrow 0}\nu_{0}^{c}(m^{2})$. Let

[TABLE]

For $x=a,b$, let $D_{\sigma_{x}}=\frac{\partial}{\partial\sigma_{x}}|_{\sigma_{a}=\sigma_{b}=0}$. For $n\geq 1$, the two-point function obeys

[TABLE]

where $\mathbb{E}_{C}$ denotes Gaussian expectation with covariance $C=((-\Delta_{\Lambda})^{\frac{\alpha}{2}}+m^{2})^{-1}$ ($m^{2}>0$ ensures existence of the inverse). Thus the two-point function is interpreted as a perturbation of a Gaussian expectation. A similar representation is valid for the weakly self-avoiding walk, using a Gaussian superexpectation.

Perturbation theory is performed inductively in a multi-scale fashion, using a finite-range decomposition $C=C_{1}+\cdots+C_{N}$, with $C_{j}$ of range $\sim L^{j}$. This is implemented via the Gaussian convolution identity $\mathbb{E}_{C}\theta=\mathbb{E}_{C_{N}}\theta\circ\cdots\circ\mathbb{E}_{C_{1}}\theta$, where $\mathbb{E}_{C}\theta$ denotes Gaussian convolution. At every step in the induction, we get a representation

[TABLE]

where the polynomial

[TABLE]

includes all Euclidean- and $O(n)$-invariant monomials that are relevant and marginal according to the RG philosophy. The error in this approximation is irrelevant in the RG sense and is controlled uniformly in the volume by parametrising it as a polymer gas. According to (1.23), after the final step of the induction has been performed, we obtain

[TABLE]

To control $q_{x,N}$ ($x=a,b$), we need to study the RG dynamical system

[TABLE]

and its non-perturbative corrections. The initial condition is $(g_{0},\nu_{0},u_{0},\lambda_{x,0},q_{x,0})=(g,\nu_{0}^{c},0,1,1,0,0)$. (In fact, the coupling constant $u_{j}$ does not play an important role for the two-point function.) For $d=4$, the dynamical system has a Gaussian fixed point. We use the adaptation of the RG method, as developed in [31], to the long-range setting below the upper critical dimension, where the fixed point is instead non-Gaussian. In [31] only the flow of $g_{j},\nu_{j},u_{j}$ was studied and $\lambda_{j},q_{j}$ did not appear, but the flow of $g_{j},\nu_{j},u_{j}$ remains identical when these additional coupling constants do appear. For the nearest-neighbour model on ${\mathbb{Z}}^{4}$, the RG method was applied in [5, 32] to prove $|a-b|^{-2}$ decay of the critical two-point function for all $n\geq 0$. We mainly follow the approach of [5, 32]. In particular, our treatment of the flow of $q_{x,j}$ remains the same and yields

[TABLE]

where $w_{j}=\sum_{k=1}^{j}C_{k}$, and where $j_{ab}=\lfloor\log_{L}(2|a-b|)\rfloor$ is the coalescence scale defined to ensure that $C_{k;a,b}=0$ when $k\leq j_{ab}$. By definition of $j_{ab}$, the right-hand side of (1.27) is zero for $j$ below the coalescence scale, and this remains true non-perturbatively as well: $q_{x,j}=0$ for scales $j\leq j_{ab}$.

The flow of $\lambda_{j}$ was analysed recursively for the Gaussian RG fixed point in [5, 32], but for the non-Gaussian fixed point in our current setting the recursive analysis cannot be applied due to the non-summability of remainder terms, and a different approach is needed. Let $\bar{D}=\sum_{x\in\Lambda}\frac{\partial}{\partial\varphi_{x}^{1}}|_{\varphi=0}$ and $\bar{D}^{2}=\sum_{x,y\in\Lambda}\frac{\partial^{2}}{\partial\varphi_{x}^{1}\partial\varphi_{y}^{1}}|_{\varphi=0}$. According to (1.23),

[TABLE]

Let $w_{j}^{(1)}=\sum_{x\in\Lambda}w_{j;a,x}$, which is independent of $a$. Using Gaussian integration by parts and translation invariance, we show in (4.7) that

[TABLE]

By using (1.23) to evaluate the two terms in the above right-hand side approximately, we thus obtain

[TABLE]

This relates $\lambda_{j}$ to the bulk coupling constants $g_{j},\nu_{j}$ whose flow is known from [31]. In particular, it is shown in [31] that $w_{j}^{(1)}\nu_{j}=O(\epsilon)$. All of the above is carried out uniformly in $m^{2}$, which permits the limit $m^{2}\downarrow 0$ to be taken after the infinite-volume limit. Since $\lim_{m^{2}\downarrow 0}\lim_{N\rightarrow\infty}w_{N:a,b}=((-\Delta_{{{{\mathbb{Z}}}^{d}}})^{\alpha/2})^{-1}_{a,b}\asymp|a-b|^{-(d-\alpha)}$, all this, together with the rigorous versions of (1.27) and (1.25), implies our main result (1.20). The non-perturbative corrections to (1.30) due to the irrelevant error coordinate are controlled using a cluster expansion. This is the main innovation in the proof of Theorem 1.1.

The remainder of the paper is organised as follows. In Section 2, we provide some background and definitions needed for the RG method. In Section 3, we formulate the RG map and state the main theorem which provides estimates on the RG map; this is an adaptation of the main result of [16] as applied to the long-range model in [31]. The main difference, compared to [31], is the inclusion of observables in the RG map. The flow of the observable coupling constant $\lambda_{j}$ is analysed in Section 4. The flow of the observable coupling constant $q_{j}$ is then analysed in Section 5.1, where the proof of Theorem 1.1 is completed.

2 Set-up for RG method

In this section, we summarise some notation and background for the RG method, needed for the proof of Theorem 1.1. Additional details can be found in [31].

2.1 Formula for two-point function

We begin with a formula for the two-point function that serves as our starting point.

2.1.1 The case $n\geq 1$

For $n\geq 1$, we define

[TABLE]

Given $g>0$, $\nu\in{\mathbb{R}}$, $m^{2}>0$, we set

[TABLE]

Given $a,b\in\Lambda$, we introduce observable fields $\sigma_{a},\sigma_{b}\in{\mathbb{R}}$, and define $V_{0}$ and $Z_{0}$ by

[TABLE]

with $V_{0}(\Lambda)=\sum_{x\in\Lambda}V_{0}(\varphi_{x})$.

Given a $\Lambda\times\Lambda$ covariance matrix $w$, let $\mathbb{E}_{w}$ denote the Gaussian expectation with covariance $w$. Let $C=((-\Delta_{\Lambda_{N}})^{\alpha/2}+m^{2})^{-1}$. By shifting part of the $|\varphi|^{2}$ term into the covariance, the expectation (1.10) can be rewritten as

[TABLE]

where $V^{\varnothing}_{0}(\Lambda)$ denotes the evaluation of $V_{0}(\Lambda)$ at $\sigma_{a}=\sigma_{b}=0$. When $F$ is a monomial, it is standard to write this ratio of expectations as a derivative of a logarithmic generating function. Let $D_{\sigma_{a}}$ denote the operator $\frac{\partial}{\partial\sigma_{a}}|_{\sigma_{a}=\sigma_{b}=0}$, and similarly for higher derivatives. Then the two-point function is given, for $n\geq 1$, by

[TABLE]

2.1.2 The case $n=0$

For $n=0$, as in several previous papers (e.g., [6, 5, 32]) we formulate the weakly self-avoiding walk model as the infinite-volume limit of a supersymmetric version of the $|\varphi|^{4}$ model. The supersymmetric model involves a complex boson field $(\phi_{x},\bar{\phi}_{x})_{x\in\Lambda}$ and a fermion field given by the 1-forms $\psi_{x}=\frac{1}{\sqrt{2\pi i}}d\phi_{x}$, $\bar{\psi}_{x}=\frac{1}{\sqrt{2\pi i}}d\bar{\phi}_{x}$. For $n=0$, in place of (2.1), we set

[TABLE]

and we replace $\varphi_{a}^{1},\varphi_{b}^{1}$ in (2.3) by $\bar{\phi}_{a},\phi_{b}$.

For $n=0$, a formula closely related to (2.5) is given, e.g., in [32, (LABEL:phi4-e:DaDbPNn0)], with $\mathbb{E}_{C}$ in (2.5) replaced by the Gaussian superexpectation. As in [32], our formalism applies to the supersymmetric model with only notational changes, with $n$ interpreted as $n=0$ in formulas such as (1.19), and with the Gaussian expectation replaced by a superexpectation. For notational simplicity, we concentrate throughout the paper on the case $n\geq 1$.

2.2 Progressive integration

In our version of the RG method, the expectation $\mathbb{E}_{C}e^{-V_{0}(\Lambda)}$ of (2.5) is evaluated in a multi-scale fashion, via a finite-range decomposition of the covariance $C$. We use the same finite-range decomposition $C=C_{1}+C_{2}+\cdots+C_{N-1}+C_{N,N}$ of the covariance $C=((-\Delta_{\Lambda_{N}})^{\alpha/2}+m^{2})^{-1}$ that is described and analysed in [31, Section LABEL:alpha-sec:frd]. A closely related decomposition was first introduced in [25] and subsequently developed in [26]. The covariances $C_{j}$ are translation invariant, and have the finite-range property

[TABLE]

Thus, we may regard $C_{j}$ either as a covariance on ${{{\mathbb{Z}}}^{d}}$ or on $\Lambda_{N}$, as long as $N>j$. Viewing the $C_{j}$ as covariances on ${{{\mathbb{Z}}}^{d}}$, we also have a decomposition of the infinite-volume covariance given by $((-\Delta)^{\alpha/2}+m^{2})^{-1}=\sum_{j=1}^{\infty}C_{j}$. We leave implicit the dependence of the covariance $C_{j}$ on $m^{2}$. According to [31, (LABEL:alpha-e:scaling-estimate)], for $m^{2}$ bounded, the covariances $C_{j}$ satisfy the estimates

[TABLE]

For $n\geq 1$, and for an integrable $F:({\mathbb{R}}^{n})^{\Lambda}\rightarrow{\mathbb{R}}$, we define the convolution $\mathbb{E}_{C}\theta F$ by

[TABLE]

where the expectation $\mathbb{E}_{C}$ on the right-hand side acts on $\zeta$ and leaves $\varphi$ fixed. A similar construction is used for $n=0$ (see, e.g., [6, Section LABEL:log-sec:ga]). By [13, Proposition LABEL:norm-prop:conv], the Gaussian convolution can be evaluated as

[TABLE]

with an abuse of notation where $C_{N}$ means $C_{N,N}$. To compute the expectation $\mathbb{E}_{C}e^{-V_{0}(\Lambda)}$ in (2.5), we use (2.10) to evaluate $\mathbb{E}_{C}\theta e^{-V_{0}(\Lambda)}$ progressively, as follows. We write $\mathbb{E}_{j}=\mathbb{E}_{C_{j}}$ and let

[TABLE]

with $Z_{0}=e^{-V_{0}(\Lambda)}$ as in (2.3). By (2.10), $\mathbb{E}_{C}Z_{0}$ is obtained by setting $\varphi=0$ in

[TABLE]

This leads us to study the recursion $Z_{j}\mapsto Z_{j+1}$.

2.3 Function space

The observable fields $\sigma_{a},\sigma_{b}$ are needed only for the purpose of evaluating the second derivative in (2.5). Therefore, dependence on the observable fields which is higher order than quadratic plays no role. We make use of this by defining the function space $\mathcal{N}$ as explained below. We also define the $T_{\varphi}$ seminorm on $\mathcal{N}$. These definitions are as in, e.g., [14, 32]. We focus on the case $n\geq 1$; the modifications needed for $n=0$ are as in, e.g., [15].

2.3.1 The space $\mathcal{N}$

Given $p_{\mathcal{N}}>0$, let

[TABLE]

As in [31, Section LABEL:alpha-sec:tfnorm], we fix any $p_{\mathcal{N}}\geq 10$. For $n=0$, $\mathcal{N}^{\varnothing}$ is instead a space of even differential forms with $p_{\mathcal{N}}$-times differentiable coefficients.

In order to treat functions of the observable fields $\sigma_{a},\sigma_{b}$, we define an extension $\mathcal{N}$ of $\mathcal{N}^{\varnothing}$ exactly as in [32, Section LABEL:phi4-sec:phi4observables_representation]. Namely, let $\mathcal{N}^{\prime}$ be the space of real-valued functions of $\varphi,\sigma_{a},\sigma_{b}$ which are $C^{p_{\mathcal{N}}}$ in $\varphi$ and $C^{\infty}$ in $\sigma_{a},\sigma_{b}$. An ideal $\mathcal{I}$ in $\mathcal{N}^{\prime}$ is formed by those elements of $\mathcal{N}^{\prime}$ whose formal power series expansion in the observable fields to order $1,\sigma_{a},\sigma_{b},\sigma_{a}\sigma_{b}$ is equal to zero. We define $\mathcal{N}$ as the quotient algebra $\mathcal{N}=\mathcal{N}^{\prime}/\mathcal{I}$. Then $\mathcal{N}$ has a direct sum decomposition

[TABLE]

where elements of $\mathcal{N}^{a},\mathcal{N}^{b},\mathcal{N}^{ab}$ are given by elements of $\mathcal{N}^{\varnothing}$ multiplied by $\sigma_{a}$, by $\sigma_{b}$, and by $\sigma_{a}\sigma_{b}$ respectively. Thus, elements of $\mathcal{N}$ can be identified with polynomials over $\mathcal{N}^{\varnothing}$ in the observable fields with terms only of order $1,\sigma_{a},\sigma_{b},\sigma_{a}\sigma_{b}$, i.e. $F\in\mathcal{N}$ can be written as

[TABLE]

with $F_{c}\in\mathcal{N}^{\varnothing}$ for each $c\in\{\varnothing,a,b,ab\}$. There are natural projections $\pi_{c}:\mathcal{N}\rightarrow\mathcal{N}^{c}$ defined for such $F$ by $\pi_{\varnothing}F=F_{\varnothing}$, $\pi_{a}F=\sigma_{a}F_{a}$, $\pi_{b}F=\sigma_{b}F_{b}$, and $\pi_{ab}F=\sigma_{a}\sigma_{b}F_{ab}$. We set $\pi_{*}=1-\pi_{\varnothing}$. The expectation $\mathbb{E}_{C}$ acts term-by-term on $F\in\mathcal{N}$, namely $\mathbb{E}_{C}F=\mathbb{E}_{C}F_{\varnothing}+\sigma_{a}\mathbb{E}_{C}F_{a}+\sigma_{b}\mathbb{E}_{C}F_{b}+\sigma_{a}\sigma_{b}\mathbb{E}_{C}F_{ab}$ for $F$ as in (2.15).

2.3.2 Seminorms

A family of seminorms is used to control the size of elements of $\mathcal{N}$. Let $\Lambda^{*}$ denote the set of sequences of any finite length (including length [math]), composed of elements of $\Lambda\times\{1,\ldots,n\}$. Let $\varphi\in({\mathbb{R}}^{n})^{\Lambda}$ be a field, and let $F\in\mathcal{N}^{\varnothing}$. Given $\vec{x}=((x_{1},i_{1}),\ldots,(x_{p},i_{p}))\in\Lambda^{*}$, we write $|\vec{x}|=p$ and let

[TABLE]

A test function $g$ is a mapping $g:\Lambda^{*}\rightarrow{\mathbb{R}}$, written $\vec{x}\mapsto g_{\vec{x}}$. We define the $\varphi$-pairing of $F$ with a test function $g$ by

[TABLE]

Given a parameter $\mathfrak{h}_{j}>0$, a scale-dependent norm $\|g\|_{\Phi}=\|g\|_{\Phi_{j}(\mathfrak{h}_{j})}$ is defined on test functions in [31, (LABEL:alpha-e:Phinorm)]. The $\Phi=\Phi_{j}(\mathfrak{h}_{j})$ norm controls the size of a test function and its discrete gradients up to order $p_{\Phi}=4$, but its precise definition is immaterial for the present discussion. With $B_{\Phi}(1)$ the unit ball in $\Phi$, we define the $T_{\varphi}=T_{\varphi,j}(\mathfrak{h}_{j})$ seminorm on $\mathcal{N}^{\varnothing}$ by

[TABLE]

Given an additional parameter $\mathfrak{h}_{\sigma}=\mathfrak{h}_{\sigma,j}$, we extend this definition to all of $\mathcal{N}$ exactly as in [14], i.e., the seminorm of $F$ of the form (2.15) is defined to be

[TABLE]

2.4 Blocks, polymers and scales

2.4.1 Blocks and polymers

The finite-range covariance decomposition is well-suited to a block decomposition of the torus $\Lambda_{N}$ of period $L^{N}$ into disjoint blocks of side $L^{j}$, for scales $0\leq j\leq N$. This decomposition is an important ingredient in our choice of the coordinates in which we represent the RG map. We now describe it in detail, along with a number of useful related definitions, as in [16].

We partition the torus $\Lambda_{N}$, which has period $L^{N}$, into disjoint $j$-blocks of side $L^{j}$ ($j\leq N$). Each $j$-block is a translate of the block $\{x\in\Lambda:0\leq x_{i}<L^{j},i=1,\ldots,d\}$. We denote the collection of $j$-blocks by $\mathcal{B}_{j}$.

A $j$-polymer is any (possibly empty) union of $j$-blocks, and $\mathcal{P}_{j}$ denotes the set of $j$-polymers. Given $X\in\mathcal{P}_{j}$, we denote by $\mathcal{B}_{j}(X)$ the set of $j$-blocks in $X$, and denote by $\mathcal{P}_{j}(X)$ the set of $j$-polymers in $X$. A nonempty polymer $X$ is connected if for any $x,x^{\prime}\in X$, there is a sequence $x=x_{0},\ldots,x_{n}=x^{\prime}\in X$ with $\|x_{i+1}-x_{i}\|_{\infty}=1$ for $i=0,\ldots,n-1$. Let $\mathcal{C}_{j}$ denote the set of connected $j$-polymers and, for any $X\in\mathcal{P}_{j}$, let ${\rm Comp}_{j}(X)\subset\mathcal{C}_{j}(X)$ be the set of connected components of $X$. The empty set $\varnothing$ is not in $\mathcal{C}_{j}$.

We say that two polymers $X,Y$ do not touch if $\min\{\|x-y\|_{\infty}:x\in X,y\in Y\}>1$. We call a connected polymer $X\in\mathcal{C}_{j}$ a small set if it consists of at most $2^{d}$ $j$-blocks, and write $\mathcal{S}_{j}$ for the collection of small sets in $\mathcal{C}_{j}$. The small-set neighbourhood of a polymer $X\in\mathcal{P}_{j}$ is $X^{\Box}=\cup_{Y\in\mathcal{S}_{j}(X):X\cap Y\neq\varnothing}Y$.

For $F_{1},F_{2}:\mathcal{P}_{j}\rightarrow\mathcal{N}$, we define the scale-$j$ circle product $F_{1}\circ F_{2}:\mathcal{P}_{j}\rightarrow\mathcal{N}$ by

[TABLE]

We only consider maps $F:\mathcal{P}_{j}\rightarrow\mathcal{N}$ with the property that $F(\varnothing)=1$. The identity element for the circle product is the map $\mathbbm{1}_{\varnothing}:\mathcal{P}_{j}\rightarrow\mathcal{N}$ defined by $\mathbbm{1}_{\varnothing}(\varnothing)=1$, $\mathbbm{1}_{\varnothing}(X)=0$ if $X\neq\varnothing$.

2.4.2 Mass and coalescence scales

Two scales play an important role for the nature of the RG recursion (2.11). We define the mass scale $j_{m}$ by

[TABLE]

By definition, $j_{m}$ is the smallest scale for which $m^{2}L^{\alpha(j_{m}-1)}\geq 1$. The mass scale is the scale beyond which the mass $m^{2}$ plays a significant helpful role in the decay of the covariance $C_{j}$. Indeed, by (2.8) and the elementary inequality (with notation $x_{+}=\max\{x,0\}$)

[TABLE]

we have

[TABLE]

We also define the coalescence scale $j_{ab}$ by

[TABLE]

By definition, $j_{ab}$ is the unique integer such that

[TABLE]

By (2.7), $C_{j;a,b}=0$ for all $j\leq j_{ab}$, and hence

[TABLE]

Ultimately, we take the limit $m^{2}\downarrow 0$ before considering large $|a-b|$, so we can and do assume that $j_{m}>j_{ab}$.

2.5 Localisation operator $\operatorname{Loc}$

We use the operator $\operatorname{Loc}$ defined and analysed in [14], to extract a local polynomial from an element $F\in\mathcal{N}$. For appropriate $X\subset\Lambda$, the local polynomial $\operatorname{Loc}_{X}F$ extracts the parts of $F$ that are relevant and marginal in the RG sense.

2.5.1 Local polynomials

The range of the operator $\operatorname{Loc}_{X}$ is a certain vector space $\mathcal{U}(X)$ of local polynomials in the field. We now define this vector space, taking into account that the elements $F\in\mathcal{N}$ to which $\operatorname{Loc}_{X}$ will be applied obey Euclidean covariance and $O(n)$ invariance on $\mathcal{N}^{\varnothing}$.

Given bulk coupling constants $g,\nu,u\in{\mathbb{R}}$; $a,b\in\Lambda$; observable fields $\sigma_{a},\sigma_{b}\in{\mathbb{R}}$; and the observable coupling constants $\lambda_{a},\lambda_{b},q_{a},q_{b}\in{\mathbb{R}}$; let

[TABLE]

The symbol $\varnothing$ denotes the bulk. (For $n=0$, we can take $u=0$ in $U^{\varnothing}$ due to supersymmetry; see [7].) For $U$ as in (2.28) and $X\subset\Lambda$, we write

[TABLE]

Let $\mathcal{U}(X)$ denote the space of polynomials of the form (2.29). Let $\mathcal{V}(X)\subset\mathcal{U}(X)$ be the subspace for which $u=q_{a}=q_{b}=0$. Note that $V_{0}$ of (2.3) obeys $V_{0}(X)\in\mathcal{V}(X)$ with $\lambda_{a}=\lambda_{b}=1$.

2.5.2 Definition of Loc

To define $\operatorname{Loc}$, we must first define a set of polynomial test functions, as in [14, Section LABEL:loc-sec:oploc]. Let $p>0$, let $a=(a_{1},\ldots,a_{p})$ with each $a_{r}\in{\mathbb{N}}_{0}^{d}$, and let $k=(k_{1},\ldots,k_{p})$ with each $k_{r}\in\{1,\ldots,n\}$. Let $\Lambda^{\prime}\subset\Lambda$ be a coordinate patch as defined in [14, Section LABEL:loc-sec:oploc] (e.g., $\Lambda^{\prime}$ can be any small set as defined in Section 2.4). Recall the set $\Lambda^{*}$ of sequences, defined in Section 2.3.2. We define a test function $q^{a,k}$, supported on sequences $\vec{x}=((x_{1},i_{1}),\ldots,(x_{p},i_{p}))\in\Lambda^{*}$ with each $x_{r}\in\Lambda^{\prime}$, by

[TABLE]

We include the case $p=0$ by interpreting (2.30) as the constant number $1$ in this case. The role of the coordinate patch, which cannot “wrap around” the torus, is to permit polynomial test functions such as (2.30) to be well-defined. We define the field dimension $[\varphi]=\frac{d-\alpha}{2}$. The dimension of the test function $q^{a,k}$ is defined to equal $p[\varphi]+\sum_{r=1}^{p}\sum_{l=1}^{d}a_{r,l}$. Given $d_{+}\geq 0$, we let $\Pi=\Pi^{d_{+}}[\Lambda^{\prime}]\subset\Phi$ denote the span of all test functions $q^{a,k}$ of dimension at most $d_{+}$.

Let $\mathcal{N}^{\varnothing}(\Lambda^{\prime})$ denote the space of functionals of the field that only depend on field values at points in $\Lambda^{\prime}$. By [14, Proposition LABEL:loc-prop:LTsymexists], for $X\subset\Lambda^{\prime}$, there is a unique operator $\operatorname{{\rm loc}}_{X}:\mathcal{N}^{\varnothing}(\Lambda^{\prime})\rightarrow\mathcal{V}(X)$ (independent of the choice of $\Lambda^{\prime}$) such that

[TABLE]

with the pairing given by (2.17) with $\varphi=0$. To extend $\operatorname{{\rm loc}}$ from $\mathcal{N}^{\varnothing}$ to $\mathcal{N}$, suppose we are given $d_{+}^{c}$ for $c=\varnothing,a,b,ab$. We let ${\rm loc}^{c}_{X}$ denote the operator ${\rm loc}_{X}$ with $d_{+}=d_{+}^{c}$ and, for $F$ as in (2.15), define

[TABLE]

It remains to specify the $d_{+}^{c}$.

In [31], for scales below the mass scale the range of the restriction of $\operatorname{Loc}$ to $\mathcal{N}^{\varnothing}$ is specified as the span of $\{1,\tau,\tau^{2}\}$. This corresponds to the choice $d_{+}^{\varnothing}=d$ (using symmetry considerations to disregard non-symmetric monomials). Above the mass scale, a different range for $\operatorname{Loc}$ is used, namely $d_{+}^{\varnothing}=d-\alpha$, due to the enhanced decay of the covariance decomposition (see [31, Section LABEL:alpha-sec:Loc]).

As in [7, Section LABEL:pt-sec:loc-specs], we set $d_{+}^{a}=d_{+}^{b}=[\varphi]$ when $\operatorname{Loc}$ acts at scales strictly less than $j_{ab}$, and set $d_{+}^{a}=d_{+}^{b}=0$ for larger scales. We always take $d_{+}^{ab}=0$.

The following elementary lemma will be useful. Let $\mathbbm{1}^{1}$ denote the constant test function supported on sequences of length $1$ and defined by $\mathbbm{1}^{1}_{(x,i)}=\delta_{i,1}$. Likewise, let $\mathbbm{1}^{2}$ denote the constant test function supported on sequences of length $2$ and defined by $\mathbbm{1}^{2}_{((x_{1},i_{1}),(x_{2},i_{2}))}=\delta_{i_{1},1}\delta_{i_{2},1}$. Note that $\mathbbm{1}^{1},\mathbbm{1}^{2}$ are each of the form (2.30), with respective dimensions $[\varphi]$ and $2[\varphi]$.

Lemma 2.1**.**

Given $x\in\Lambda$ and a coordinate patch $\Lambda^{\prime}\ni x$, suppose that $F\in\mathcal{N}^{\varnothing}(\Lambda^{\prime})$. For $j<N$ and $m=1,2$,

[TABLE]

Moreover, if $j<j_{ab}$, then for $c=a,b$,

[TABLE]

Proof.

The first statement is an immediate consequence of the definition of ${\rm Loc}$ together with the fact that the test function $\mathbbm{1}^{m}$ has dimension $m[\varphi]\leq d-\alpha\leq d_{+}^{\varnothing}$ for $m=1,2$. The second statement follows similarly using the fact that the dimension $[\varphi]$ of $\mathbbm{1}^{1}$ is equal to $d_{+}^{a}=d_{+}^{b}$ if $j<j_{ab}$.

3 RG map

In the absence of observables, i.e., with $\sigma_{a}=\sigma_{b}=0$, the RG map for the long-range models is constructed and bounded in [31] using the main theorem of [16]. The result is given in [31, Theorem LABEL:alpha-thm:step-mr]. The extension of this construction to the case of nonzero observable fields $\sigma_{a},\sigma_{b}$ follows a similar route as in the 4-dimensional nearest neighbour case in [5, 32], as we now explain. The coordinates for the RG map are discussed in Section 3.1, the domain of the RG map is discussed in Section 3.2, and the main estimates for the RG map are given in Theorem 3.3. These estimates, combined with a new estimate derived from a cluster expansion, are used in Sections 4–5 to control the flow generated by the RG map.

3.1 RG coordinates

The RG map will be defined so as to express the sequence $Z_{j}$ defined by (2.11) as

[TABLE]

for a real sequence $\zeta_{j}$ and sequences of maps $I_{j}:\mathcal{P}_{j}\rightarrow\mathcal{N}$ and $K_{j}:\mathcal{P}_{j}\rightarrow\mathcal{N}$. The perturbative coordinate $I_{j}$ is an explicit function of $V_{j}\in\mathcal{V}$, and

[TABLE]

The nonperturbative coordinate $K_{j}$ is discussed in detail below. By (2.3), (3.1) holds at scale $j=0$ with $u_{0}=q_{x,0}=0$, $I_{0}=e^{-V_{0}}$, and $K_{0}=\mathbbm{1}_{\varnothing}$. We sometimes write an element of $\mathcal{U}$ as $U=(\zeta,V)$ with $V\in\mathcal{V}$, where $\zeta$ encodes $u,q_{a},q_{b}$.

We express the map $Z_{j}\mapsto Z_{j+1}$ of (2.11) via a map $(V_{j},K_{j})\mapsto(\delta\zeta_{j+1},V_{j+1},K_{j+1})$, the renormalisation group (RG) map, in such a manner that

[TABLE]

with $I_{j+1}=I_{j+1}(V_{j+1})$ and $\delta\zeta_{j+1}=\zeta_{j+1}-\zeta_{j}$. This ensures that $Z_{j+1}$ has the form (3.1) with $\zeta_{j+1}=\zeta_{j}+\delta\zeta_{j+1}$.

3.1.1 Perturbative coordinate

The form of the perturbative coordinate $I_{j}$ is as follows. Given a $\Lambda\times\Lambda$ matrix $w$, we define the operator $\mathcal{L}_{w}=\frac{1}{2}\sum_{u,v\in\Lambda}w_{u,v}\sum_{i=1}^{n}\frac{\partial}{\partial\varphi_{u}^{i}}\frac{\partial}{\partial\varphi_{v}^{i}}$. Recall the projections defined in Section 2.3.1. Given $V^{\prime},V^{\prime\prime}\in\mathcal{V}$, we also define

[TABLE]

For $j\geq 0$, we write the partial sums of the covariance decomposition as

[TABLE]

As in [7, (LABEL:pt-e:WLTF)], for $B\in\mathcal{B}_{j}$ we define

[TABLE]

The polynomial $W_{j}(V,B)$ in the fields is thus an explicit quadratic function of $V$. In particular, $W_{j}^{\varnothing}$ is an even polynomial in the fields, and $W_{j}$ is quadratic in the coupling constants and is irrelevant in the RG sense. Finally, for $V\in\mathcal{V}$, we define $I_{j}=I_{j}(V,\cdot\,):\mathcal{P}_{j}\rightarrow\mathcal{N}$ by

[TABLE]

As in (2.5), we write $D_{\sigma_{a}}=\frac{\partial}{\partial\sigma_{a}}|_{\sigma_{a}=\sigma_{b}=0}$. We also write $\bar{D}=\sum_{x\in\Lambda}\frac{\partial}{\partial\varphi_{x}^{1}}|_{\varphi=0}$ and $\bar{D}^{2}=\sum_{x,y\in\Lambda}\frac{\partial^{2}}{\partial\varphi_{x}^{1}\partial\varphi_{y}^{1}}|_{\varphi=0}$. We will later make use of the following corollary of Lemma 2.1.

Corollary 3.1**.**

For $V\in\mathcal{V}$ and $x\in\Lambda$, and with $W_{j}=W_{j}(V,x)$,

[TABLE]

Moreover, if $j<j_{ab}$, then

[TABLE]

Proof.

The fact that $\bar{D}W^{\varnothing}_{j}=0$ is immediate since $W_{j}^{\varnothing}$ is even in the fields. Also, since $\bar{D}^{2}W_{j}^{\varnothing}=\langle W^{\varnothing}_{j},\mathbbm{1}^{2}\rangle_{0}$, it follows from Lemma 2.1 that $\bar{D}^{2}W_{j}^{\varnothing}=0$. To see that $\bar{D}D_{\sigma_{a}}W_{j}=0$, note that $D_{\sigma_{a}}W_{j}=\tfrac{1}{2}(1-{\rm loc}^{a}_{x})F_{a}$, with $F_{a}\in\mathcal{N}^{\varnothing}$ the coefficient of $\sigma_{a}$ in $F_{\pi,w_{j}}(V_{x},V(\Lambda))$. Thus, by definition (2.17) of the pairing, $\bar{D}D_{\sigma_{a}}W_{j}=\tfrac{1}{2}\langle(1-{\rm loc}^{a}_{x})F_{a},\mathbbm{1}^{1}\rangle_{0}$, which vanishes by Lemma 2.1 when $j<j_{ab}$. This completes the proof.

3.1.2 Nonperturbative coordinate

We now define the space $\mathcal{K}_{j}$ of maps $K:\mathcal{P}_{j}\rightarrow\mathcal{N}$ which contains the nonperturbative RG coordinate. With $\mathcal{N}$ replaced by $\mathcal{N}^{\varnothing}$, such a space is defined in [31, Definition LABEL:alpha-def:Kspace], and, as in [16], we extend it here to include observables. The symmetries (Euclidean covariance, gauge invariance, supersymmetry, and $O(n)$-invariance) used in Definition 3.2 are defined in [16, Section LABEL:step-sec:coordinates] and [4, Section LABEL:phi4-log-sec:margrel]. For $n\geq 1$, as a replacement for the gauge invariance which holds for $n=0$ we also introduce sign invariance, which is invariance under the map $(\sigma,\varphi)\mapsto(-\sigma,-\varphi)$. Note that $V_{0}$ of (2.3) is sign invariant. It can be verified that the property of sign invariance is preserved by the map $K\mapsto K_{+}$ of [16].

Definition 3.2**.**

For $j<N$, let $\mathcal{C}\mathcal{K}_{j}=\mathcal{C}\mathcal{K}_{j}(\Lambda)$ denote the real vector space of functions $K:\mathcal{C}_{j}\rightarrow\mathcal{N}$ with the following properties:

• •

Field Locality: For all $X\in\mathcal{P}_{j}(\Lambda)$, $K(X)\in\mathcal{N}(X^{\Box})$. Also, (i) $\pi_{a}K(X)=0$ unless $a\in X$, (ii) $\pi_{b}K(X)=0$ unless $b\in X$, and (iii) $\pi_{ab}K(X)=0$ unless $a\in X$ and $b\in X^{\Box}$ or vice versa, and $\pi_{ab}K(X)=0$ if $X\in\mathcal{S}_{j}$ and $j<j_{ab}$.

• •

Symmetry: (i) $\pi_{\varnothing}K$ is Euclidean covariant, (ii) if $n=0$, $\pi_{\varnothing}K$ is supersymmetric and $K$ is gauge invariant and has no constant part; if $n\geq 1$, $\pi_{\varnothing}K$ is $O(n)$-invariant and $K$ is sign invariant.

Let $\mathcal{K}_{j}=\mathcal{K}_{j}(\Lambda)$ be the real vector space of functions $K:\mathcal{P}_{j}\rightarrow\mathcal{N}$ which have the above field locality and symmetry properties, and, in addition:

• •

Component Factorisation: for all polymers $X$, $K(X)=\prod_{Y\in{\rm Comp}(X)}K(Y)$.

The nonperturbative coordinate $K_{j}$ appearing in (3.1) is an element of $\mathcal{K}_{j}$. An element of $\mathcal{K}_{j}$ determines an element of $\mathcal{C}\mathcal{K}_{j}$ by restriction to $X\in\mathcal{C}_{j}$. Also, an element of $\mathcal{C}\mathcal{K}_{j}$ determines an element of $\mathcal{K}_{j}$ by the factorisation property. The same symbol is used for both elements related by this correspondence. Since the empty set is not a connected set, $\mathbbm{1}_{\varnothing}\in\mathcal{K}_{j}$ becomes $0\in\mathcal{C}\mathcal{K}_{j}$ under this correspondence.

3.2 Norms and RG domain

We now specify the domain of the RG map, which requires specification of norms on the spaces $\mathcal{V}$ and $\mathcal{C}\mathcal{K}$. Without the observables fields, the norms are discussed in [31, Section LABEL:alpha-sec:nr]. For the nearest-neighbour 4-dimensional case, the adaptation of the norms to include observables is discussed, e.g., in [32, Section LABEL:phi4-sec:par]. For our current long-range setting, we need only adjust some norm parameters, compared to [32, Section LABEL:phi4-sec:par].

As in [31, (LABEL:alpha-e:sfix)], the small number $\bar{s}$ in Theorem 1.1 is given by

[TABLE]

with the constant $a$ specified in [31, Lemma LABEL:alpha-lem:beta-a0] (not to be confused with the point $a\in{{{\mathbb{Z}}}^{d}}$ used for the two-point function). Recall the mass scale $j_{m}$ defined in (2.21). Following [31, (LABEL:alpha-e:alphapdef)–(LABEL:alpha-e:elldefa), (LABEL:alpha-e:hdef)], we fix

[TABLE]

and define the bulk parameters

[TABLE]

Here $\ell_{0}$ can be chosen large (depending on $L$) and $k_{0}$ is a fixed (small) constant. We use $\mathfrak{h}_{j}$ to refer to either of the bulk parameters $\ell_{j},h_{j}$.

Now that observables are present, the pair of parameters $\mathfrak{h}_{j}$ is supplemented by the pair

[TABLE]

We only use $h_{\sigma,j}$ for $j\leq j_{m}$. Recall that we assume that the coalescence scale $j_{ab}$ is smaller than the mass scale $j_{m}$, since the limit $m^{2}\downarrow 0$ will be taken before considering arbitrarily large $|a-b|$.

For $U\in\mathcal{U}\simeq{\mathbb{R}}^{7}$, we define the scale-dependent norm

[TABLE]

We denote the restriction of $\|\cdot\|_{\mathcal{U}}$ to $\mathcal{V}$ by the same symbol. Given $C_{\mathcal{D}}>0$, we define the domain

[TABLE]

Note that $\mathcal{D}_{j}$ is a domain in $\mathcal{V}$, and as such, does not involve the coupling constants $u$ or $q$.

A sequence $\mathcal{W}^{\varnothing}_{j}$ of Banach spaces is defined in terms of the $T_{\varphi}(\mathfrak{h}_{j})$ seminorms in [31] (they are denoted $\mathcal{W}_{j}$ there). We extend $\mathcal{W}^{\varnothing}_{j}$ to a space $\mathcal{W}_{j}\subset\mathcal{C}\mathcal{K}_{j}$ whose definition is the same with the exception that the $T_{\varphi}$ seminorms are defined on the extended space $\mathcal{N}$. As in [31, Remark LABEL:alpha-rk:chiL], we define a sequence

[TABLE]

Given a parameter $t>0$, the domain of the RG map is defined by

[TABLE]

where $B_{\mathcal{W}_{j}}(r)$ is the open ball of radius $r$ in the Banach space $\mathcal{W}_{j}$.

3.3 Estimates on RG map

We now specify the RG map $(V_{j},K_{j})\mapsto(U_{j+1},K_{j+1})=(\delta\zeta_{j+1},V_{j+1},K_{j+1})$ and state our bounds on it. To shorten notation, we condense indices and write, e.g., $(V,K)$ for $(V_{j},K_{j})$ and $(U_{+},K_{+})$ for $(U_{j+1},K_{j+1})$. The definition of the maps $U_{+},K_{+}$ is described in a general setup in [7, 16], and is adapted to the long-range model with $\sigma_{a}=\sigma_{b}=0$ in [31]. The same definitions extend to include observables.

In particular, the map $(V,K)\mapsto U_{+}=(\delta\zeta_{+},V_{+})=\operatorname{PT}(V)+R_{+}(V,K)$ is explicit and consists of a perturbative part $\operatorname{PT}$, incorporating second-order perturbation theory, and a nonperturbative, third-order error $R_{+}$. The explicit map $\operatorname{PT}$ is the one defined in [7] for $n=0$, extended in [4] to $n\geq 1$, and used in [32] for general $n\geq 0$. Let $\lambda$ denote $\lambda_{a}$ or $\lambda_{b}$, and let $q$ denote $q_{a}$ or $q_{b}$. We denote the $\lambda,q$ components of the map $\operatorname{PT}$ by $\lambda_{\rm pt},\delta q_{\rm pt}$. For $j<N$, and with $w_{j}$ given by (3.6), let

[TABLE]

By [31, (LABEL:alpha-e:Greeknoprime) and Lemma LABEL:alpha-lem:wlims],

[TABLE]

Let

[TABLE]

Recall the definition of the coalescence scale $j_{ab}$ in (2.24). Then, as in [32, Proposition LABEL:phi4-prop:pt], for general $n\geq 0$ the observable part of the map $\operatorname{PT}$ is the map $V\mapsto(\lambda_{\rm pt},\delta q_{\rm pt})$ given by

[TABLE]

Note that $\lambda_{\rm pt}=\lambda$ for all scales $j\geq j_{ab}-1$, i.e., the flow of $\lambda$ stops evolving after scale $j_{ab}-1$. Conversely, since $C_{j+1;a,b}=0$ for $j+1\leq j_{ab}$, nonzero $\delta q_{{\rm pt}}$ can occur only at scales $j\geq j_{ab}$. The map $(V,K)\mapsto U_{+}$ is now defined by

[TABLE]

The localisation operator ${\rm Loc}_{Y,B}$ is defined in [14, Definition LABEL:loc-def:LTXYsym]. The higher-order correction $R_{+}:\mathcal{V}\rightarrow\mathcal{U}$ to the perturbative calculation is then defined by

[TABLE]

so that $U_{+}=\operatorname{PT}(V)+R_{+}(V,K)$. We do not need the explicit form of $R_{+}$ and only use the bounds of Theorem 3.3 below.

The map $(V,K)\rightarrow K_{+}$ is also given explicitly in [16], but it is complicated to write down. Like $R_{+}$, this nonperturbative part of the RG map is of order $O(\bar{s}^{3})$. It is part of the statement of Theorem 3.3 below that the formula for $K_{+}$ constructed in [16] is well-defined on the domain specified in Theorem 3.3. We do not need to know more here about $K_{+}$ than the estimates provided by Theorem 3.3.

The RG map depends on the mass $m^{2}$ through its dependence on the covariance $C_{+}$. We require continuity in the mass in the limit $m^{2}\downarrow 0$, which can only be taken after the infinite-volume limit $N\rightarrow\infty$. Given small $\delta>0$, we define the mass domain for the RG map by

[TABLE]

The special attention to $j=N$ is due to the fact that the final covariance $C_{N,N}$ is only defined for $m^{2}>0$, and it obeys good estimates for $m^{2}\in\mathbb{I}_{N}$.

The following extends [31, Theorem LABEL:alpha-thm:step-mr] to allow for the presence of observables. Its estimates appear identical to [31, Theorem LABEL:alpha-thm:step-mr], but it is in fact an extension since the domain and range of the RG map now include observables in $(V,K)\in\mathbb{D}$, as well as in $R_{+}$ and $K_{+}$. Note that the map $R_{+}$, which acts on $(V,K)$ with $V\in\mathcal{V}$, produces a polynomial in $\mathcal{U}$ which in particular contains the nonperturbative contributions to $\delta\zeta$. The bound (3.29) on $R_{+}$ controls these nonperturbative contributions to $\delta\zeta$. Note that the estimates (3.29) hold for $m^{2}\in\mathbb{I}_{+}$, but the continuity is in the smaller interval $m^{2}\in[0,L^{-\alpha j}]$. A restriction like this on the continuity interval is essential, because larger $m^{2}$ will put $j$ above the mass scale, at which point the spaces themselves become dependent on $m^{2}$ through their dependence on $\ell_{j}$ and a continuity statement becomes meaningless.

Theorem 3.3**.**

Let $d=1,2,3$; $n\geq 0$; $\alpha=\frac{1}{2}(d+\epsilon)$ and $j<N$. Let $C_{\mathcal{D}}$ and $L$ be sufficiently large, and let $\epsilon$ be sufficiently small. There exist $C_{\rm RG}>0$, $\delta>0$, such that, with the domain $\mathbb{D}$ defined using $t=4C_{\rm RG}$, the maps

[TABLE]

are analytic in $(V,K)$, provide a solution to (3.3), and satisfy the estimates

[TABLE]

The coordinate in $R_{+}$ corresponding to $\delta q_{a,j},\delta q_{b,j}$ is identically zero for $j\leq j_{ab}$, and the coordinate corresponding to $\lambda_{a,j},\lambda_{b,j}$ is identically zero for $j\geq j_{ab}$. In addition, $R_{+},K_{+}$ are jointly continuous in $m^{2}\in[0,L^{-\alpha j}],V,K$.

Proof.

The theorem is a consequence of the main result of [16], which focusses on the 4-dimensional nearest-neighbour case. For the long-range model, the appropriate modifications for the bulk part of the RG map are discussed in [31], and we assume familiarity with both the methodology and the modifications. In order to include observables, only minor further modifications are required, compared to [15, 16].

One requirement is to verify that, for $V\in\mathcal{D}$, the basic small parameters $\epsilon_{V}$ and $\bar{\epsilon}$ obey appropriate estimates when observables are present, as in [15, Sections LABEL:IE-sec:epVW–LABEL:IE-sec:epdV-app]. We verify this here; this verification validates our choice (3.15) for the norm parameters. (In fact, somewhat larger domains are used in [15, Sections LABEL:IE-sec:epVW–LABEL:IE-sec:epdV-app]; the main ideas are present for $V\in\mathcal{D}$, which we consider here, and the extension to the larger domains is a matter of bookkeeping.) A second requirement is to verify that the “crucial contraction” is maintained in the presence of observables, and we also verify this here.

Bound on $\epsilon_{V}$. Let $V\in\mathcal{D}$. For $\epsilon_{V}$, it suffices to observe that for $|\lambda_{a}|\leq C_{\mathcal{D}}\bar{s}\ell^{-1}\ell_{\sigma}^{-1}$,

[TABLE]

which implies stability on the domain $\mathcal{D}$ of (3.17), and complements the arguments of [16, 31].

Bound on $\bar{\epsilon}$. We must also verify the analogue of [15, Lemma LABEL:IE-lem:epdV]. To state the desired estimate, as in [31, (LABEL:alpha-e:ellhatdef)] we define the norm parameter

[TABLE]

and as in [31, (LABEL:alpha-e:epdVdef)] we define the small parameter

[TABLE]

We write $U_{\rm pt}=\operatorname{PT}(V)$ and let $\delta V=\theta V-U_{\rm pt}$. Our goal then is to show that, for $V\in\mathcal{D}$,

[TABLE]

where $C_{\delta V}$ is an $L$-dependent constant.

It is argued in [31, Section LABEL:alpha-sec:epdV] that (3.33) holds with $\delta V$ replaced by $\delta V^{\varnothing}=\pi_{\varnothing}\delta V$. Thus, it suffices to establish (3.33) with $\delta V$ replaced by $\delta V^{*}=\pi_{*}\delta V$. This can be done by writing

[TABLE]

and applying the triangle inequality to estimate each of the two terms on the right-hand side.

For instance, if $a\in B$, then the $\sigma_{a}$ term of $\theta V^{*}(B)-V^{*}(B)$ is $\lambda_{a}\sigma_{a}\zeta^{1}_{a}$. By definition of the norm, by (3.31), (3.15), (3.32), (3.18), and by the fact that $\alpha^{\prime}<\tfrac{1}{2}\alpha$,

[TABLE]

The $\sigma_{a}$ term of $V^{*}(B)-U_{\rm pt}^{*}(B)$ is zero above the coalescence scale, whereas if $j+1<j_{ab}$ then it is $\delta[\nu w^{(1)}]\lambda_{a}\sigma_{a}\varphi^{1}_{a}$, by (3.23). Thus, by (3.30), it is sufficient to show that

[TABLE]

By its definition in (3.22),

[TABLE]

By (3.17) and (2.23), and the finite-range property (2.7), the first term is bounded by

[TABLE]

and the second term is bounded by

[TABLE]

These bounds do better than what is required by (3.36).

For the $\sigma_{a}\sigma_{b}$ term, we can take $j\geq j_{ab}$. The $\sigma_{a}\sigma_{b}$ term of $\theta V^{*}-V^{*}$ is always [math] and the coefficient of $\sigma_{a}\sigma_{b}$ in $V^{*}(B)-U_{\rm pt}^{*}(B)$ is at most $|C_{+;a,b}\lambda_{a}\lambda_{b}|$. By (3.17), (2.23), and (3.13) and the fact that $\alpha^{\prime}<\frac{1}{2}\alpha$,

[TABLE]

When $\mathfrak{h}=h$ (hence $j<j_{m}$) this is $O(\bar{\epsilon}(h)^{2})$, and when $\mathfrak{h}=\ell$ it is $O(\bar{\epsilon}(\ell)^{2})$. This is better than what is required for (3.33).

Crucial contraction. The adaptation of the crucial contraction to the long-range model is provided for the bulk in [31, Section LABEL:alpha-sec:kappabms–LABEL:alpha-sec:kappapms]. We now extend the adaptation to include observables.

Below the mass scale, the least irrelevant of the sign invariant monomials involving the observable fields each have two additional spin fields compared to their marginal counterparts $\sigma_{a}\varphi_{a}^{1}$ and $\sigma_{a}\sigma_{b}$ (the latter occurs only above the coalescence scale), so have dimension which is larger by $2[\varphi]=d-\alpha$. Compared to [31, (LABEL:alpha-e:kappa)], this gives rise to $\gamma=L^{-(d-\alpha)}$, and there is no factor $L^{d}$ for observables, so the gain here is proportional to $L^{-(d-\alpha)}$. The worst $\gamma$ occurs for $d=1$, where we have $\gamma=L^{-\frac{d}{2}+\frac{\epsilon}{2}}=L^{-\frac{1}{2}+\frac{\epsilon}{2}}$. This is consistent with the values of $L^{d}\gamma$ reported for the bulk in [31, (LABEL:alpha-e:kappa)].

Above the mass scale, we extend the discussion in [31, Section LABEL:alpha-sec:kappapms], as follows. For the perturbative contribution to $K$, we have already verified that we can continue to use the $\bar{\epsilon}$ given by (3.32) when observables are present, and there is therefore no change to [31] concerning this issue. It remains to consider the crucial contraction.

We recall and invoke our assumption that $j_{ab}<j_{m}$. Now $d_{+}^{a}=0$, so the least irrelevant monomial in $\mathcal{N}^{a}$ is $\sigma_{a}\varphi$. This scales as

[TABLE]

A change from scale $j$ to scale $j+1$ in the above right-hand side gives rise to a factor $2L^{-\frac{1}{2}(d+\alpha^{\prime})}$. As in [31, (LABEL:alpha-e:kappajm)–(LABEL:alpha-e:kappa-above-jm-5)], the essential condition here is that the product of this factor with ${\vartheta}^{-3}=L^{\frac{3}{4}\alpha}$ should be bounded above by an inverse power of $L$. This condition is indeed satisfied, since

[TABLE]

Similarly, the least irrelevant monomial in $\mathcal{N}^{ab}$ that is sign invariant is of the form $\sigma_{a}\sigma_{b}\varphi\varphi$, and has scaling dimension twice that considered in the previous paragraph, so twice as good. Thus the crucial contraction is not harmed by the presence of observables.

Estimate for $R_{+}$ above the mass scale. Finally, we consider the extension of [31, Section LABEL:alpha-sec:Rbound] to include observables. The observable terms have the same $T_{0}$ and $\mathcal{U}$ norms: $\|\sigma_{a}\varphi_{a}^{1}\|_{T_{0}}=\ell_{\sigma}\ell=\|\sigma_{a}\varphi_{a}^{1}\|_{\mathcal{U}}$ and $\|\sigma_{a}\sigma_{b}\|_{T_{0}}=\ell_{\sigma}^{2}=\|\sigma_{a}\sigma_{b}\|_{\mathcal{U}}$. This leads to an extension to [31, Lemma LABEL:alpha-lem:monnormcomp], as follows. Let

[TABLE]

The estimates of [31, Lemma LABEL:alpha-lem:monnormcomp] now become

[TABLE]

i.e., the bound remains the same for $F_{1}$ but loses a helpful factor $L^{-\alpha^{\prime}(j-j_{m})_{+}}$ for $F_{2}$. The bound on $F_{1}$ then implies, as in [31, (LABEL:alpha-e:Rnormcomp)], that

[TABLE]

and the bound (3.29) follows from this as in [31, Section LABEL:alpha-sec:Rbound].

The introduction of observables does lead to a change in the bounds on $F,W,P$ in [31, Lemma LABEL:alpha-lem:W1], due to the weakened estimate for $F_{2}$ in (3.45). The change is to replace the factors $L^{-(\alpha+\alpha^{\prime})(j-j_{m})_{+}}$ and $(c/L)^{-(\alpha+\alpha^{\prime})(j-j_{m})_{+}}$ in the three upper bounds of [31, Lemma LABEL:alpha-lem:W1] by the worse factor ${\vartheta}^{2}=L^{-\frac{1}{2}\alpha(j-j_{m})_{+}}$. Since we seek an upper bound which includes the factor ${\vartheta}^{2}$ in $\bar{\epsilon}^{2}$, the weakened bounds remain more than good enough.

For general reasons, $\pi_{ab}W=0$ [15, Proposition LABEL:IE-prop:Wbounds], so there can be no such term in $W$. Thus, in the proof of [31, Lemma LABEL:alpha-lem:W1], only one factor $L^{-\alpha^{\prime}(j-j_{m})_{+}}$ can be lost by application of (3.45), not two. Also, by direct calculation, the relevant contribution to $F$ is $F_{C}(\lambda_{a}\sigma_{a}\varphi_{a}^{1},\lambda_{b}\sigma_{b}\varphi_{b}^{1})=\textstyle{\frac{1}{2}}\lambda_{a}\lambda_{b}C_{a,b}\sigma_{a}\sigma_{b}$, whose $T_{0}(\ell)$ norm is given as in the first inequality of (3.40) to be at most $L^{-(\alpha-\alpha^{\prime})(j-j_{m})}$, which is better than the required ${\vartheta}^{2}$. The bound on $P$ follows from the bounds on $F,W$ as in [15, Proposition 4.1].

In the absence of observables, Theorem 3.3 is used in [31] to construct a global RG flow $(g_{j},\nu_{j},K_{j}^{\varnothing})$ that remains in the RG domain for all $j$. This requires tuning the initial $\nu$ to a mass-dependent critical value $\nu_{0}^{c}(m^{2})$; this value converges to the critical point $\nu_{c}(g;n)$ as $m^{2}\downarrow 0$ (see [31, (LABEL:alpha-e:nuceta)–(LABEL:alpha-e:nuc)]). Throughout the remainder of the present paper, we always take $(g_{j},\nu_{j})$ to be this global flow of coupling constants. For general reasons this flow is the same in the presence of observables as in their absence: see [16, (LABEL:step-e:piVKplus)–(LABEL:step-e:plusindep)]. The main task for the proof of Theorem 1.1 is to apply the estimates of Theorem 3.3 to control, in addition, the flow of the observable coupling constants $\lambda$ and $\delta q$, and the observable part of the coordinate $K$. The flow of $\delta q$ and $K$ is analysed as in the 4-dimensional nearest-neighbour case [5, 32].

The flow of $\lambda$ is marginal, for the same reasons as in the 4-dimensional case. In [5, 32], the perturbative approximation (3.23) to the recursion for $\lambda$ is solved along the lines of the rough computation

[TABLE]

In [5, 32], the errors introduced by the map $R_{+}$ into (3.23) were summable over all scales because of the decay of the marginal coupling constant $g_{j}$ with the scale (Gaussian fixed point), and the above computation survives the introduction of these errors.

For the long-range model, the fixed point is non-Gaussian, and the corrections due to $R_{+}$ are not summable. Instead of trying to follow the route laid out in [5, 32], we derive an exact relation between $\lambda_{a,j}$ and the known bulk coupling constants, similar to (3.47), which gives better control of its flow than the recursion. This is done in Section 4.

4 Flow of $\lambda$

According to (3.23) and Theorem 3.3, the flow of $\lambda_{a,j}$ under the RG map is nontrivial only until scale $j_{ab}-1$, and stops beyond this scale. Conversely, $q_{a,j}=0$ for $j<j_{ab}$, and the flow of $q_{a,j}$ is nontrivial only for scales $j\geq j_{ab}$. Our goal now is to determine the form of the flow until scale $j_{ab}$. Since we later take the limit $m^{2}\downarrow 0$ before studying large $j_{ab}\sim\log_{L}|a-b|$, we can and do assume that $j_{ab}<j_{m}$. We will prove the following proposition.

Proposition 4.1**.**

Let $n\geq 0$, let $L$ be sufficiently large, let $\Lambda_{N}$ be the torus of period $L^{N}$, and let $\epsilon$ be sufficiently small. Let $g\in[\frac{63}{64}\bar{s},\frac{65}{64}\bar{s}]$, and let $m^{2}\in[L^{-\alpha(N-1)},\delta]$ with $\delta>0$ sufficiently small. Let $j_{ab}<j_{m}<N$. Let $g_{0}=g$, let $\nu_{0}$ be the critical value $\nu_{0}^{c}(m^{2})$ constructed in [31], and let $\lambda_{a,0}=\lambda_{b,0}=1$. Then the RG map can be iterated to scale $j_{ab}$, i.e., it produces a sequence $(V_{j},K_{j})\in\mathbb{D}_{j}$ with initial condition $(V_{0},0)$, such that (3.1) holds for all $j\leq j_{ab}$ with $I_{j}=I_{j}(V_{j})$ and $\zeta_{j}=\sum_{k=1}^{j}\delta\zeta_{j}$. Moreover, $q_{x,j}=0$, and for the component $\lambda_{x,j}$ of this flow we have the stronger statement

[TABLE]

The proof of Proposition 4.1 is given in Section 4.1 below. Its statement holds trivially at $j=0$, and will be established inductively for higher scales. The induction for the bulk quantities $U_{j}^{\varnothing},K_{j}^{\varnothing}$ is the result of [31], and is unaffected by the presence of observables.

The main additional ingredient for the induction of the observable parts is to establish the flow (4.1) of $\lambda_{a,j}$. To achieve this, in Lemma 4.2 we use integration by parts to obtain a relation between $\lambda_{a,j}$, quantities of the bulk flow, and the observable parts of the coordinate $K_{j}$. This is achieved by taking suitable derivatives of the identity $Z_{j}=\mathbb{E}_{w_{j}}\theta Z_{0}$. The contribution due to $K_{j}$ is bounded uniformly in the volume using a cluster expansion, in Section 4.2.

The formula (4.1) for $\lambda_{x,j}$ has a natural counterpart for the nearest-neighbour 4-dimensional case, with error term $O(g_{j}^{2})$ instead of $O(\bar{s}^{2})$. In that context, $-\nu_{j}w_{j}^{(1)}+O(g_{j}^{2})\rightarrow 0$ as $j\rightarrow\infty$. This provides insight into the fact that $\lim_{j\rightarrow\infty}\lambda_{x,j}=1$ in [5, Lemma LABEL:saw4-lem:lamlim] and [32, Corollary LABEL:phi4-cor:vx]. For the long-range model considered in Proposition 4.1, the non-Gaussian fixed point leads to a limit which is not equal to $1$.

4.1 Integration by parts

For notational convenience we restrict attention to $n\geq 1$; small modifications apply for $n=0$. Recall that $\bar{D}$ and $D_{\sigma_{a}}$ are defined above Corollary 3.1, and that $Z_{j}=\mathbb{E}_{w_{j}}\theta Z_{0}$. Let

[TABLE]

Then we have

[TABLE]

The existence of the logarithm $\mathcal{L}_{j}$ is discussed in Section 4.2, where it is constructed as an element of a Banach space $T_{0}(\ell_{j})$ which only examines derivatives at zero field, using a cluster expansion. Bounds on $\mathcal{L}_{j}$ and its derivatives at zero field are proved in Proposition 4.4 below.

Lemma 4.2**.**

The functions $I_{j}$ and $\mathcal{L}_{j}$ are related by the identity

[TABLE]

Proof.

By definition, followed by Gaussian integration by parts,

[TABLE]

On the right-hand side, $\frac{\partial}{\partial\zeta_{y}^{1}}$ can be replaced by $\frac{\partial}{\partial\varphi_{y}^{1}}$, and the latter commutes with the expectation. Then application of $\bar{D}=\sum_{x\in\Lambda}\frac{\partial}{\partial\varphi_{x}^{1}}|_{\varphi=0}$ gives

[TABLE]

which by translation invariance and by definition of $Z_{j}$ is the same as

[TABLE]

Now we divide both sides of (4.7) by $Z_{j}^{\varnothing}|_{\varphi=0}$ and use (4.3). Since $I_{j}|_{\varphi=0}=1$, and since $\bar{D}Z^{\varnothing}_{j}=\bar{D}I_{j}|_{\sigma_{a}=\sigma_{b}=0}=D_{\sigma_{a}}I_{j}|_{\varphi=0}=0$ by symmetry, the result is (4.4).

Note that the right-hand side of (4.4) involves only bulk quantities, while the left-hand side depends on $\lambda_{a,j}$ through $I_{j}(\Lambda)$ and $D_{\sigma_{a}}\mathcal{L}_{j}(\Lambda)$, and also on the observable part of the irrelevant coordinate $K_{j}$ (through $D_{\sigma_{a}}\mathcal{L}_{j}(\Lambda)$). For the explicit terms, we have the following identities.

Lemma 4.3**.**

For $j\leq N$ and $V\in\mathcal{V}$,

[TABLE]

and if $j<j_{ab}$ then

[TABLE]

Proof.

We differentiate the formula $I_{j}(\Lambda)=e^{-V_{j}(\Lambda)}\prod_{B\in\mathcal{B}_{j}(\Lambda)}(1+W_{j}(V_{j},B))$, which is (3.8). We apply the product rule, Corollary 3.1, and the facts that $\bar{D}V_{j}^{\varnothing}=0$ and $I_{j}|_{\varphi=0}=1$, to obtain

[TABLE]

Similarly, for $j<j_{ab}$, we also use $D_{\sigma_{a}}V_{j}(\Lambda)=-\lambda_{a,j}\varphi^{1}_{a}$ to obtain

[TABLE]

and the proof is complete.

We now state our bounds for the terms in (4.4) involving $\mathcal{L}_{j}$. The hypothesis $(V_{j},K_{j})\in\mathbb{D}_{j}$ of Proposition 4.4 will be established inductively.

Proposition 4.4**.**

Let $j\leq j_{ab}$, let $\mathcal{L}_{j}$ be defined as in (4.2), and assume that $Z_{j}=e^{\zeta_{j}}(I_{j}\circ K_{j})$ with $I_{j}=I_{j}(V_{j})$ and $(V_{j},K_{j})\in\mathbb{D}_{j}$. Then there is a constant $c_{1}>0$ such that

[TABLE]

We defer the proof of Proposition 4.4 to Section 4.2.

Proof of Proposition 4.1.

The proof is by induction on $j$. The statement of Proposition 4.1 for $j=0$ is trivial. Without loss of generality, we consider the case $x=a$. We assume that we have (3.1) for $Z_{k}$ with $(V_{k},K_{k})$ constructed inductively using the RG map for $k\leq j$, and we make the constant in the hypothesis (4.1) explicit by assuming that, with $c_{1}$ from (4.12),

[TABLE]

Then we have (3.1) with a pair of RG coordinates $(V_{j},K_{j})\in\mathbb{D}_{j}$, satisfying in addition (4.13). Theorem 3.3 guarantees the existence of RG coordinates $(U_{j+1},K_{j+1})=(\delta\zeta_{j+1}=\delta u_{j+1},V_{j+1},K_{j+1})$ at scale $j+1$ such that $Z_{j+1}$ obeys (3.1), with $U_{j+1}=\operatorname{PT}_{j}(V_{j})+R_{j+1}(V_{j},K_{j})$, and bounds on $R_{j+1}(V_{j},K_{j})$ and $K_{j+1}$ as in (3.29).

It has been proved in [31] that the bulk part of $V_{j+1}$ lies in $\mathbb{D}_{j+1}$. The second bound in (3.29) is sufficient to guarantee that $K_{j+1}$ also lies in $\mathbb{D}_{j+1}$. Therefore, to complete the proof that $(V_{j+1},K_{j+1})\in\mathbb{D}_{j+1}$, we only need to show that $|\lambda_{j+1}|<C_{\mathcal{D}}$, where $C_{\mathcal{D}}>1$ is the constant in (3.17). By (3.23), and by the first bound of (3.29) together with the definition of the norm in (3.16), we have $\lambda_{j+1}=(1+O(\bar{s}))\lambda_{j}+O(\bar{s}^{2})$. It now follows immediately from (4.13) that $0<\lambda_{j+1}=1+O(\bar{s})<C_{\mathcal{D}}$, since $\bar{s}$ can be chosen small enough. This proves that $(V_{j+1},K_{j+1})\in\mathbb{D}_{j+1}$.

To complete the induction, we must prove (4.13) with $j$ replaced by $j+1$. Since (4.1) is only required for scales below the coalescence scale, we may assume here that $j+1<j_{ab}$. The bounds of Proposition 4.4 at scale $j+1$ can be applied, since the hypothesis $(V_{j+1},K_{j+1})\in\mathbb{D}_{j+1}$ has now been verified. Also, the hypothesis $j+1<j_{ab}$ of Lemma 4.3 is satisfied. We use Lemma 4.3 in conjunction with (4.4), and apply Proposition 4.4. This gives

[TABLE]

by (3.21) and by taking $\bar{s}$ sufficiently small. This advances (4.13) to scale $j+1$, and completes the proof.

4.2 Cluster expansion

In this section, we use a cluster expansion to construct a formula for $\mathcal{L}_{j}=\log z_{j}$ and prove Proposition 4.4. Let $p(X)=K_{j}(X)/I_{j}(X)$. By (3.1), (4.2), and by definition of the circle product,

[TABLE]

where the term in the sum with $X=\varnothing$ is interpreted as $1$. In the sum, we decompose $X\in\mathcal{P}_{j}$ into its connected components $X_{1},\ldots,X_{n}\in\mathcal{C}_{j}$, which may be labelled in $n!$ different ways. For $X,X^{\prime}\in\mathcal{C}_{j}$, we set $g(X,X^{\prime})=-1$ if $X$ and $X^{\prime}$ touch, and otherwise set $g(X,X^{\prime})=0$. Using the component factorisation property of $K_{j}$, we obtain

[TABLE]

where the $n=0$ term is again interpreted as $1$. This has the form of the partition function of a polymer system, as defined, e.g., in [34, (1)]. It is a standard result, e.g., [34, 19], that $\log z_{j}$ can be written as a cluster expansion and accurately bounded, provided the polymer activities $p(X)$ obey suitable estimates. In the following proof, we discuss this in detail and invoke a convergence criterion from [34]; see also [21, 30] for pedagogical introductions to the cluster expansion. The verification of the criterion from [34] is an almost immediate consequence of the norm estimates in the definition of the domain $\mathbb{D}_{j}$.

Since we are interested only in the derivatives of $\mathcal{L}_{j}$ at zero external and observable fields, we do not construct $\mathcal{L}_{j}$ as a function of these fields (even though this would also be possible for suitably small fields), but rather as a Taylor polynomial (jet) of order $p_{\mathcal{N}}$ in the fields around zero. In other words, we work on the quotient of $\mathcal{N}$ by the ideal of elements of $F\in\mathcal{N}$ with $\|F\|_{T_{0}(\ell_{j})}=0$. On this quotient, the $T_{0}(\ell_{j})$ seminorm becomes a norm, and the quotient becomes a finite-dimensional Banach algebra. This is discussed in detail in [16, Section LABEL:step-sec:Knorms, Appendix LABEL:step-sec:Banach]. We adopt the point of view in the following that we work in this normed space, and write simply $\|\cdot\|$ for $\|\cdot\|_{T_{0}(\ell_{j})}$. Although the results of [34] are stated for complex-valued $p(X)$, the proofs hold verbatim for values in any Banach algebra. The completeness of the Banach algebra is important for the existence of $\mathcal{L}_{j}=\log z_{j}$, which is defined in terms of an infinite sum.

The estimates we use, for $(V_{j},K_{j})\in\mathbb{D}_{j}$ and $X\in\mathcal{C}_{j}$, are:

[TABLE]

where $a>0$ is small; here, $|X|=|X|_{j}$ denotes the number of $j$-blocks in $X$. The bound on $I_{j}^{-X}$ is a small adaptation of [15, Proposition LABEL:IE-prop:Istab] to our long-range setting, and the bound on $K_{j}$ follows from the definition of the $\mathcal{W}_{j}$ norm and (3.19). Absorbing the factor $2^{|X|}$ by replacing $M$ by $M^{\prime}>M$, replacing $a$ by $a^{\prime}\in(0,a)$, and using the fact that $\bar{s}$ is sufficiently small, we conclude that the polymer activity obeys the bound

[TABLE]

The following lemma uses this bound and will be employed to verify the hypothesis of [34, Theorem 1].

Lemma 4.5**.**

If $B\in\mathcal{B}_{j}$, then for $\bar{s}$ sufficiently small (depending only on $d$)

[TABLE]

where the constant depends on $d$.

Proof.

The number of connected polymers $Y\in\mathcal{C}_{j}$ that touch a block $Y$ and have size $|Y|=n$ is at most $A^{n}$ for some $d$-dependent constant $A$. Thus,

[TABLE]

We split the sum on the right-hand side into sums with $n\leq 2^{d}$ and $n>2^{d}$. The first of these is a constant that depends only on $d$. Taking $\bar{s}$ small so that $Ae\bar{s}^{a^{\prime}}<1$, the second sum is bounded as

[TABLE]

with a $d$-dependent constant, which suffices.

Proof of Proposition 4.4.

By Lemma 4.5, if $\bar{s}$ is sufficiently small, then

[TABLE]

which verifies the hypothesis [34, (3)] with $a(A)=|A|$. Also by Lemma 4.5,

[TABLE]

which verifies the other hypothesis of [34, Theorem 1].

Let $u(X_{1},\ldots,X_{n})$ denote the Ursell function, defined in [34, (2)] (with $g$ written as $\zeta$). We conclude from [34, Theorem 1] that $\mathcal{L}_{j}$ is given by the absolutely convergent sum

[TABLE]

and that for all $X_{1}\in\mathcal{C}_{j}$ we have

[TABLE]

with the $n=1$ term in (4.25) interpreted as $1$.

By (4.24)–(4.25) (the factor $n$ in (4.25) is not needed) and (4.23),

[TABLE]

Similarly, for $\|D_{\sigma_{a}}\mathcal{L}_{j}\|$, we use the product rule for differentiation, this time using the factor $n$ (due to the product rule) in (4.25). With the definition of the $T_{0}$ seminorm in (2.19) and of $\ell_{\sigma,j}$ in (3.15) for $j\leq j_{ab}$, we obtain

[TABLE]

For $F\in\mathcal{N}$, $\bar{D}F=\langle F,\mathbbm{1}^{1}\rangle_{0}$ and $\bar{D}^{2}F=\langle F,\mathbbm{1}^{2}\rangle_{0}$, with the test functions $\mathbbm{1}^{1},\mathbbm{1}^{2}$ of Lemma 2.1. These two test functions have $\Phi_{j}$-norms (as defined, e.g., in [31, (LABEL:alpha-e:Phinorm)])

[TABLE]

Therefore, for $m=1,2$,

[TABLE]

In particular, since $L^{Nd}=|\Lambda|$,

[TABLE]

and the proof is complete.

5 Full RG flow and proof of Theorem 1.1

In Proposition 4.1, the RG flow $(\zeta_{j},V_{j},K_{j})$ is constructed for scales $j\leq j_{ab}$. The sequence $\zeta_{j}$ of (3.2) contains in particular the coupling constants $q_{a,j},q_{b,j}$; recall that $q_{x,j}=0$ for $j\leq j_{ab}$. In Section 5.1, we apply Theorem 3.3 inductively to continue the RG flow $(\zeta_{j},V_{j},K_{j})$ to scales $j_{ab}<j\leq N$. Using the extended flow, we prove Theorem 1.1 in Section 5.2. The analysis proceeds as in [5, 32].

Once the RG flow has been extended to all scales, the combination of (2.12) and (3.1) gives, at the final scale $j=N$, the representation

[TABLE]

From this, we apply (2.5) to to calculate the two-point function as

[TABLE]

with

[TABLE]

5.1 Flow of $q$

The next proposition states that the RG flow exists for scales $j_{ab}\leq j\leq N$, and in particular analyses the flow of $q$ and establishes control on the terms of the right-hand side of (5.2), which is needed to prove Theorem 1.1.

Proposition 5.1**.**

Let $n\geq 0$, let $L$ be sufficiently large, let $\Lambda_{N}$ be the torus of period $L^{N}$, and let $\epsilon$ be sufficiently small. Let $g\in[\frac{63}{64}\bar{s},\frac{65}{64}\bar{s}]$, and let $m^{2}\in[L^{-\alpha(N-1)},\delta]$ with $\delta>0$ sufficiently small. Suppose that $j_{ab}<j_{m}$. Starting with $(V_{j_{ab}},K_{j_{ab}})$ produced by Proposition 4.1, the RG map can be iterated to scale $N$, i.e., it produces a sequence $(V_{j},K_{j})\in\mathbb{D}_{j}$ such that (3.1) holds for all $j\leq N$ with $I_{j}=I_{j}(V_{j})$ and $\zeta_{j}=\sum_{k=1}^{j}\delta\zeta_{j}$. The $q_{x,j}$ component of $\zeta_{j}$ is given by

[TABLE]

with

[TABLE]

Moreover,

[TABLE]

Proof.

For $j=j_{ab}$, we have $(V_{j_{ab}},K_{j_{ab}})\in\mathbb{D}_{j_{ab}}$ by Proposition 4.1. Also, (5.4)–(5.5) hold trivially, since $r_{x,j_{ab}}=0$ by Theorem 3.3 and hence $q_{x,j_{ab}}=\lambda_{a,j_{ab}}\lambda_{b,j_{ab}}w_{j;a,b}$ by (3.24).

We fix $j\geq j_{ab}$ and assume inductively that (3.1) holds with a pair of RG coordinates $(V_{j},K_{j})\in\mathbb{D}_{j}$ and that (5.4)–(5.5) hold. As in the proof of Proposition 4.1, Theorem 3.3 guarantees the existence of RG coordinates $(V_{j+1},K_{j+1})$ at scale $j+1$, with $V_{j+1}=\operatorname{PT}_{j+1}(V_{j})+R_{j+1}(V_{j},K_{j})$, and bounds on $R_{j+1}(V_{j},K_{j})$ and $K_{j+1}$ as in (3.29).

As before, it has been proved in [31] that the bulk part of $V_{j+1}$ lies in $\mathbb{D}_{j+1}$. The coordinate $\lambda_{a,j}=\lambda_{a,j_{ab}}$ remains constant for $j>j_{ab}$, and thus still lies in $\mathbb{D}_{j+1}$. As before, the second bound in (3.29) is sufficient to guarantee that $K_{j+1}$ also lies in $\mathbb{D}_{j+1}$. This shows that $(V_{j+1},K_{j+1})\in\mathbb{D}_{j+1}$.

We now show that $q_{a,j+1}$ satisfies (5.4) at scale $j+1$ and that (5.5) holds. Using (5.4) and (3.24), and denoting by $r_{a,j}$ the component of $R_{j+1}(U_{j},K_{j})$ corresponding to the component $q_{a}$, we see that

[TABLE]

verifying (5.4) at scale $j+1$. By definition of the norm in (3.16) and by our assumption that $j_{ab}<j_{m}$, Theorem 3.3 gives the bound

[TABLE]

which proves (5.5) since $L^{-j_{ab}(d-\alpha)}=O(|a-b|^{-(d-\alpha)})$ by (2.25).

Finally, we write $D^{k}_{\sigma}$ to mean no derivative for $k=0$, the derivative with respect to $\sigma_{a}$ or $\sigma_{b}$ for $k=1$, and the second derivative with respect to $\sigma_{a},\sigma_{b}$ for $k=2$. Since $(V_{N},K_{N})\in\mathbb{D}_{N}$, it follows from (3.29), with the fact that the $\mathcal{W}_{N}$ norm bounds the $T_{0}(\ell_{N})$ norm, that

[TABLE]

Since ${\vartheta}_{N}\rightarrow 0$ as $N\rightarrow\infty$, this implies (5.6) and completes the proof.

5.2 Proof of Theorem 1.1

With (5.2) and Proposition 5.1, it is now straightforward to complete the proof of our main result. In the proof, we write $C_{a,b}$ for the massless free two-point function $((-\Delta)^{\alpha/2})^{-1}_{a,b}$ on ${{{\mathbb{Z}}}^{d}}$. According to (1.6), $C_{a,b}\asymp|a-b|^{-(d-\alpha)}$. The proof uses the following lemma.

Lemma 5.2**.**

Let $\nu_{0}=\nu_{0}^{c}(m^{2})$. Then for any $j<N$ the map $m^{2}\mapsto(V_{j},K_{j})$ is continuous for $m^{2}\in[0,L^{-\alpha j}]$. Moreover, the sequence $V_{j}$ is independent of $N$. In particular, for any $j<\infty$, the maps $V_{j},K_{j},R_{+,j}$ depend continuously on $m^{2}$ at $m^{2}=0$.

Proof.

We show by induction that $(V_{j},K_{j})$ depends continuously on $m^{2}\in[0,L^{-\alpha j}]$. The case $j=0$ follows from [31, (LABEL:alpha-e:nu0)] and [31, Corollary LABEL:alpha-cor:mu0]. Now suppose the inductive hypothesis holds for some $j\geq 0$. Then the case $j+1$ follows from (3.23), [31, Lemma LABEL:alpha-lem:wlims], and Theorem 3.3. The fact that $V_{j}$ is independent of $N$ is [16, Proposition LABEL:step-prop:VZd].

Proof of Theorem 1.1.

We first take the limit $N\rightarrow\infty$, then take the limit $m^{2}\downarrow 0$, and finally consider large $|a-b|$. By (5.4) with $j=N$,

[TABLE]

By Proposition 5.1, the remainder term is bounded uniformly in $N$ and in $m^{2}\in[L^{-\alpha(N-1)},\delta]$ by

[TABLE]

By dominated convergence, and by the continuity of $r_{x,i}$ (a component of $R_{i+1}$) at $m^{2}=0$ guaranteed by Lemma 5.2, $\lim_{m^{2}\downarrow 0}\lim_{N\rightarrow\infty}\sum_{i=j_{ab}}^{N-1}r_{x,i}$ exists and is bounded by $O(\bar{s})C_{a,b}$. For the main term, since $\lambda_{j_{ab}}=1+O(\bar{s})$ by Proposition 4.1, it follows from the definition of $w_{N}$ in (3.6) (together with the fact that that the covariance appearing in (3.24) is always the infinite-volume one), that

[TABLE]

The existence of the above limit as $m^{2}\downarrow 0$ is a consequence of the fact that $w_{\infty,a,b}=((-\Lambda)^{\alpha/2}+m^{2})^{-1}_{ab}\rightarrow C_{a,b}$, together with the mass continuity of $\lambda_{j_{ab}}$, which follows from Lemma 5.2. We apply (5.6) in (5.2), and find that the critical two-point function obeys

[TABLE]

This completes the proof.

Acknowledgements

This work was supported in part by NSERC of Canada. We thank Slava Rychkov for helpful correspondence, and an anonymous referee for useful suggestions.