Renormalizing the Kardar-Parisi-Zhang equation in $d\geq 3$ in weak disorder

TL;DR

This paper investigates the behavior of the KPZ equation in three or more dimensions under weak noise, identifying its limiting distribution and establishing tail bounds and moment existence.

Contribution

It characterizes the limiting distribution of the KPZ equation in higher dimensions with weak noise and provides strong approximations and tail bounds.

Findings

Limit law identified for the KPZ equation in $d\,\geq\,3$

Limit has sub-Gaussian lower tails

All moments exist for the limit distribution

Abstract

We study Kardar-Parisi-Zhang equation in spatial dimension 3 or larger driven by a Gaussian space-time white noise with a small convolution in space. When the noise intensity is small, it is known that the solutions converge to a random limit as the smoothing parameter is turned off. We identify this limit, in the case of general initial conditions ranging from flat to droplet. We provide strong approximations of the solution which obey exactly the limit law. We prove that this limit has sub-Gaussian lower tails, implying existence of all negative (and positive) moments.

Equations

$$\frac{\partial}{\partial t}h=\frac{1}{2}\Delta h+\bigg{[}\frac{1}{2}|\nabla h|^{2}-\infty\bigg{]}+\xi\qquad u\colon{\mathbb{R}}^{d}\times\mathbb{R}_{+}\to\mathbb{R}$$

$$\xi (\phi )={\int }_0^{\infty }{\int }_{{\mathbb{R}}^d}\text{d}t\text{d}x\xi (t,x)\phi (t,x)$$

$$\mathbb{E}\big{[}\xi(\varphi_{1})\,\,\xi(\varphi_{2})\big{]}=\int_{0}^{\infty}\int_{{\mathbb{R}}^{d}}\varphi_{1}(t,x)\varphi_{2}(t,x)\,\text{\rm d}t\,\text{\rm d}x.$$

$$\frac{\partial}{\partial t}h_{\varepsilon}=\frac{1}{2}\Delta h_{\varepsilon}+\bigg{[}\frac{1}{2}|\nabla h_{\varepsilon}|^{2}-C_{\varepsilon}\bigg{]}+\beta\varepsilon^{\frac{d-2}{2}}\xi_{\varepsilon}\;,\qquad\,\,h_{\varepsilon}(0,x)=0,$$

$${\xi }_{\epsilon }(t,x)=(\xi (t,\cdot )\star {\varphi }_{\epsilon })(x)=\int {\varphi }_{\epsilon }(x-y)\xi (t,y)\text{d}y.$$

$${{\mathbb{E}[\xi_{\varepsilon}(t,x)\xi_{\varepsilon}(s,y)]=\delta({t-s})\,\varepsilon^{-d}V\big{(}(x-y)/\varepsilon\big{)},}}$$

$$\xi^{{\scriptscriptstyle{({\varepsilon,t,x}})}}(\varphi)\stackrel{{\scriptstyle\mathrm{(def)}}}{{=}}\,\,\varepsilon^{-\frac{d+2}{2}}\int_{\mathbb{R}}\int_{{\mathbb{R}}^{d}}\varphi\big{(}\varepsilon^{-2}(t-s),\varepsilon^{-1}(y-x)\big{)}\xi(s,y)ds\,dy.$$

$$\xi^{{\scriptscriptstyle{({\varepsilon,t,x}})}}(s,y)=\varepsilon^{\frac{d+2}{2}}\xi\bigg{(}\varepsilon^{2}\Big{(}\frac{t}{\varepsilon^{2}}-s\Big{)},\varepsilon\Big{(}y-\frac{x}{\varepsilon}\Big{)}\bigg{)}$$

$${\xi }^{(\epsilon ,t)}={\xi }^{(\epsilon ,t,0)}.$$

$$u_{\epsilon }=\exp h_{\epsilon }$$

$$\frac{\partial }{\partial t}{u}_{\epsilon }=\frac{1}{2}\Delta u_{\epsilon }+\beta {\epsilon }^{\frac{d-2}{2}}{u}_{\epsilon }{\xi }_{\epsilon },u_{\epsilon }(0,x)=1,$$

$$C_{\epsilon }={\beta }^2(\varphi \star \varphi )(0){\epsilon }^{-2}/2={\beta }^{2}V(0){\epsilon }^{-2}/2$$

$$u_{\varepsilon}(t,x)=E_{x}\bigg{[}\exp\bigg{\{}\beta\varepsilon^{\frac{d-2}{2}}\,\int_{0}^{t}\int_{{\mathbb{R}}^{d}}\,\phi_{\varepsilon}(W_{{t-s}}-y)\xi(s,y)\,\text{\rm d}s\,\text{\rm d}y-\frac{\beta^{2}\,t\,\varepsilon^{-2}}{2}\,\,V(0)\bigg{\}}\bigg{]}\;.$$

$$u_{\epsilon }(t,x)={\mathcal{Z}}_{\frac{t}{{\epsilon }^2}}\left({\xi }^{(\epsilon ,t)};\frac{x}{\epsilon }\right)$$

$$\mathscr{Z}_{T}(x)=\mathscr{Z}_{T}(\xi;x)=E_{x}\bigg{[}\exp\bigg{\{}\beta\,\int_{0}^{T}\int_{{\mathbb{R}}^{d}}\,\phi(W_{{s}}-y)\xi(s,y)\,\,\text{\rm d}s\,\text{\rm d}y-\frac{\beta^{2}\,T}{2}\,\,V(0)\bigg{\}}\bigg{]},$$

$${\mathcal{Z}}_T(x)\to \{\begin{array}{ll}{\mathcal{Z}}_{\infty }(x)& \text{mbox}if\beta \in (0,{\beta }_c),\\ 0& \text{mbox}if\beta \in ({\beta }_c,\infty ).\end{array}$$

$$\mathfrak{u}={\mathfrak{u}}_{\beta ,\varphi }:{\mathcal{C}}^{\alpha }(\mathbb{R}\times {\mathbb{R}}^d)\to (0,\infty ),$$

$$\mathfrak{h}=\log \mathfrak{u}.$$

$$\beta_{L^{2}}=\sup\left\{\beta>0:E_{0}\bigg{[}e^{\beta^{2}\int_{0}^{\infty}V(\sqrt{2}W_{s})\,\,\text{\rm d}s}\bigg{]}<\infty\right\}$$

$$h_{\varepsilon}(t,x)-{\mathfrak{h}}\big{(}\xi^{{\scriptscriptstyle{({\varepsilon,t,x}})}}\big{)}\stackrel{{\scriptstyle\mathbb{P}}}{{\longrightarrow}}0\;.$$

$$h_{\epsilon }(t,x)-\mathfrak{h}({\xi }^{(\epsilon ,t,x)})-\log \overline{u}(t,x)\stackrel{\mathbb{P}}{⟶}0,$$

$${\partial }_t\overline{u}=\frac{1}{2}\Delta \overline{u},\overline{u}(0,x)=\exp h_0(x).$$

$$\underset{t\searrow 0}{\lim }\exp h_{\epsilon }(t,\cdot )={\delta }_{x_0}(\cdot )$$

$$h_{\epsilon }(t,x)-\mathfrak{h}({\xi }^{(\epsilon ,t,x)})-\mathfrak{h}({\xi }_{(\epsilon ,x_0)})-\log \rho (t,x-x_0)\stackrel{\mathbb{P}}{⟶}0,$$

$${\xi }_{(\epsilon ,x_0)}(s,x)={\epsilon }^{\frac{d+2}{2}}\xi ({\epsilon }^{2}s,x_0+\epsilon x)$$

$$\underset{\beta \to 0}{\lim }\underset{\epsilon \to 0}{\lim }{h}_{\epsilon }=\underset{\epsilon \to 0}{\lim }\underset{\beta \to 0}{\lim }{h}_{\epsilon }.$$

$$\underset{e^{h_0}\to {\delta }_{x_0}}{\lim }\underset{\epsilon \to 0}{\lim }{h}_{\epsilon }^{(h_0)}\ne \underset{\epsilon \to 0}{\lim }\underset{e^{h_0}\to {\delta }_{x_0}}{\lim }{h}_{\epsilon }^{(h_0)}$$

$$\mathbb{P}[\mathfrak{h}\le -\theta ]\le C{\mathrm{e}}^{-{\theta }^2/2}.$$

$$\underset{e^{h_0}\to {\delta }_{x_0}}{\lim }\underset{\epsilon \to 0}{\lim }\mathbb{E}{h}_{\epsilon }^{(h_0)}-\underset{\epsilon \to 0}{\lim }\underset{e^{h_0}\to {\delta }_{x_0}}{\lim }\mathbb{E}{h}_{\epsilon }^{(h_0)}=-\mathbb{E}\mathfrak{h}>0.$$

$${\sigma }^2(\beta )=C_d{\int }_{{\mathbb{R}}^d}\text{d}yV(\sqrt{2}y)E_y[\exp \{{\beta }^2{\int }_0^{\infty }V(\sqrt{2}{W}_s)\text{d}s\}]$$

Full text

${}^{~{}1}$Université Paris Diderot

Laboratoire de Probabilités, Statistique et Modélisation

LPSM (UMR 8001 CNRS, SU, UPD)

Bâtiment Sophie Germain, 8 place Aurélie Nemours, 75013 Paris

comets@lpsm.paris, ccosco@lpsm.paris

${}^{~{}2}$University of Münster

Fachbereich Mathematik und Informatik

Einsteinstraße 62, Münster, D-48149

chiranjib.mukherjee@uni-muenster.de

${}^{~{}3}$ NYU-ECNU Institute of Mathematical Sciences at NYU Shanghai

**Keywords: SPDE, Kardar-Parisi-Zhang equation, stochastic heat equation, directed polymers, random environment, weak disorder, Edwards-Wilkinson limit

AMS 2010 subject classifications: Primary 60K35. Secondary 35R60, 35Q82, 60H15, 82D60**

1 Introduction and main results.

1.1 KPZ equation and its regularization.

We consider the Kardar-Parisi-Zhang (KPZ) equation written informally as

[TABLE]

and driven by a totally uncorrelated Gaussian space-time white noise $\xi$. More precisely, $\xi$ on ${\mathbb{R}}_{+}\times{\mathbb{R}}^{d}$ is a family $\{\xi(\varphi)\}_{\varphi\in\mathcal{S}({\mathbb{R}}_{+}\times{\mathbb{R}}^{d})}$ of Gaussian random variables

[TABLE]

with mean 0 and covariance

[TABLE]

The equation (1.1) describes the evolution of a growing interface in $d+1$ dimension [17, 22] and also appears as the scaling limit for $d=1$ of front propagation of the certain exclusion processes ([3, 9]) as well as that of the free energy of the discrete directed polymer ([1]). It should be noted that, on a rigorous level, only distribution-valued solutions are expected for (1.1), and thus it is already ill-posed in $d=1$ stemming from the inherent non-linearity of the equation and the fundamental problem of squaring or multiplying random distributions. For $d=1$, studies related to the above equation have enjoyed a huge resurgence of interest in the last decade starting with the important work [13] which gave an intrinsic precise notion of a solution to (1.1).

We now fix a spatial dimension $d\geq 3$. As remarked earlier, since (1.1) is a-priori ill-posed, we will study its regularized version

[TABLE]

which is driven by the spatially mollified noise

[TABLE]

with $\phi_{\varepsilon}(\cdot)=\varepsilon^{-d}\phi(\cdot/\varepsilon)$ being a suitable approximation of the Dirac measure $\delta_{0}$ and $C_{\varepsilon}$ being a suitable divergent (renormalization) constant. We will work with a fixed mollifier $\phi:{\mathbb{R}}^{d}\to{\mathbb{R}}_{+}$ which is smooth and spherically symmetric, with $\mathrm{supp}(\phi)\subset B(0,\frac{1}{2})$ and $\int_{{\mathbb{R}}^{d}}\phi(x)\,\text{\rm d}x=1$. Then, $\{\xi_{\varepsilon}(t,x)\}$ is a centered Gaussian field with covariance

[TABLE]

where $V=\phi\star\phi$ is a smooth function supported in $B(0,1)$. We also remark that in (1.2), the multiplicative parameter $\beta$ can be taken to be positive without loss of generality, while by rescaling, no multiplicative parameter is needed in (1.1), see [21]. Also in spatial dimensions $d\geq 3$, the factor $\varepsilon^{\frac{d-2}{2}}$ is the correct scaling – a small enough $\beta>0$ guarantees a non-trivial random limit of $h_{\varepsilon}$ as $\varepsilon\to 0$, see the discussion in Section 1.3.

The goal of the present article is to consider general solutions of (1.2), namely the solutions of (1.2) with various initial conditions and prove that as the mollification parameter $\varepsilon$ is turned off, the renormalized solution of (1.2) converges to a meaningful random limit as long as $\beta$ remains small enough. We use Feynman-Kac representation of the solution of stochastic heat equation and results from directed polymers. Not only do we identify the distributional limit of $h_{\varepsilon}$, but we also provide a sequence (indexed by $\varepsilon$) of functions of the noise such that

• •

it is a strong approximation of $h_{\varepsilon}$, i.e. the difference tends to 0 in norm,

• •

all terms in the sequence obey the limit law.

The above functions for the flat initial condition are defined from the martingale limit of a random polymer model taken at some rescaled, shifted and time-reversed version of the noise. The similar functions for other intial conditions can be derived from the martingale limit taken at various version of the noise and the heat kernel. We finally show that it has sub-Gaussian lower tails in this regime, which implies existence of all negative and positive moments of this object. Besides new contributions, we gather and reformulate results which are atomized in the literature, often stated in a primitive form and hidden by necessary technicalities. We end up the introduction with a rather complete account on the state-of-the-art. We now turn to a more precise description of our main results.

1.2 Main results.

In order to state our main results, we will introduce the following notation which will be consistently used throughout the sequel. Recall the definition of the space-time white noise $\xi\in\mathcal{S}^{\prime}(\mathbb{R}\times{\mathbb{R}}^{d})$ which is a random tempered distribution (defined in all times, including negative ones), and for any $\varphi\in\mathcal{S}(\mathbb{R}\times{\mathbb{R}}^{d})$, $\varepsilon>0$, $t\in\mathbb{R}$ and $x\in{\mathbb{R}}^{d}$,

[TABLE]

Equivalently,

[TABLE]

so that by invariance under space-time diffusive rescaling, time-reversal and spatially translation, $\xi^{{\scriptscriptstyle{({\varepsilon,t,x}})}}$ is itself a Gaussian white noise and possesses the same law as $\xi$. This is also the reason why we define the noise above also for negative times. To abbreviate notation, we will also write

[TABLE]

We also need specify the definition(s) of the critical disorder parameter. Note that (1.2) is inherently non-linear. The Hopf-Cole transformation suggests that

[TABLE]

solves the linear multiplicative noise stochastic heat equation (SHE)

[TABLE]

provided that the stochastic integral in (1.6) is interpreted in the classical Itô-Skorohod sense and that we choose

[TABLE]

equal to the Itô correction below. Then, the generalized Feynman-Kac formula ([16, Theorem 6.2.5]) provides a solution to (1.6)

[TABLE]

with $E_{x}$ denoting expectation with respect to the law $P_{x}$ of a $d$-dimensional Brownian path $W=(W_{s})_{s\geq 0}$ starting at $x\in{\mathbb{R}}^{d}$, which is independent of the noise $\xi$. Plugging (1.3) in the previous formula, using Brownian scaling and time-reversal, we get the a.s. equality

[TABLE]

where

[TABLE]

is the normalized partition function of the Brownian directed polymer is a white noise environment $\xi$, or equivalently, the total-mass of a Gaussian multiplicative chaos in the Wiener space ([20, Section 4]).

It follows that there exists $\beta_{c}\in(0,\infty)$ and a strictly positive non-degenerate random variable $\mathscr{Z}_{\infty}(x)$ so that, a.s. as $T\to\infty$,

[TABLE]

See [20], or [7] for a general reference. Moreover, $(\mathscr{Z}_{T})_{T\geq 0}$ is uniformly integrable for $\beta<\beta_{c}$, which we will always assume from now on. Now, let $\mathcal{C}^{\alpha}(\mathbb{R}\times{\mathbb{R}}^{d})$ denote the path-space of the white noise (see Appendix for a precise definition) and

[TABLE]

be any arbitrary representative of the random limit $\mathscr{Z}_{\infty}=\mathscr{Z}_{\infty}(0)$; in particular $\mathfrak{u}(\xi)=\mathscr{Z}_{\infty}$. Then, $\mathbb{E}[\mathfrak{u}]=1$, and throughout the sequel we will write (recall (1.5) and (1.8))

[TABLE]

Since $\mathfrak{u}$ is non constant with $\mathbb{E}\mathfrak{u}=1$, we have $\mathbb{E}\mathfrak{h}<0$.

Finally, we also define another critical disorder parameter:

[TABLE]

which corresponds to the $L^{2}$-region of the polymer model (see (1.14)). In $d\geq 3$, it is easy to see that for $\beta$ small enough, $\sup_{x\in{\mathbb{R}}^{d}}E_{x}[\beta\int_{0}^{\infty}V(W_{s})\,\,\text{\rm d}s]<1$, so that by Khas’minskii’s lemma, E_{0}\big{[}\exp\big{\{}\beta\int_{0}^{\infty}V(W_{s})\,\,\text{\rm d}s\big{\}}\big{]}<\infty, so this implies that $\beta_{L^{2}}>0$. Furthermore, for $\beta<\beta_{L^{2}}$, convergence (1.10) becomes an $L^{2}$-convergence, hence $0<\beta_{L^{2}}<\beta_{c}<\infty$.

We are now ready to state our main results.

Theorem 1.1**.**

Assume $d\geq 3$ and recall that $\mathfrak{h}$ is defined in (1.11).

• •

(Flat initial condition.) Fix $\beta\in(0,\beta_{c})$ and consider the solution $h_{\varepsilon}$ to (1.2) with $h_{\varepsilon}(0,\cdot)=0$. Then, for all $t>0,x\in{\mathbb{R}}^{d}$, we have as $\varepsilon\to 0$,

[TABLE]

• •

(General initial condition.) Fix $\beta\in(0,\beta_{L^{2}})$ and consider the solution $h_{\varepsilon}$ to (1.2) with $h_{\varepsilon}(0,\cdot)=h_{0}(\cdot)$ for some $h_{0}:{\mathbb{R}}^{d}\to{\mathbb{R}}$ which is continuous and bounded from above. Then, for all $t>0,x\in{\mathbb{R}}^{d}$, we have as $\varepsilon\to 0$,

[TABLE]

where

[TABLE]

• •

(Droplet or narrow-wedge initial condition.) Fix $\beta\in(0,\beta_{L^{2}})$ and consider the solution $h_{\varepsilon}$ to (1.2) such that

[TABLE]

for some $x_{0}\in{\mathbb{R}}^{d}$. Then, for all $t>0,x\in{\mathbb{R}}^{d}$, we have as $\varepsilon\to 0$,

[TABLE]

where $\rho$ is the $d$-dimensional Gaussian kernel, and

[TABLE]

is a space-time Gaussian white noise.

The deterministic terms in (1.12) and (1.13) are logarithms of solutions to heat equation without noise, and one can see that ${\mathfrak{h}}\to 0$ in $L^{2}$ as $\beta\to 0$. This implies that

[TABLE]

We obtain an immediate corollary to Theorem 1.1.

Corollary 1.2**.**

Fix $\beta\in(0,\beta_{L^{2}})$ and denote by $h_{\varepsilon}^{{\scriptscriptstyle{({h_{0}}})}}$ the solution of (1.2) with initial condition $h_{\varepsilon}(0,\cdot)=h_{0}(\cdot)$, where $h_{0}$ is continuous and bounded from above. Then for any $x_{0}\in{\mathbb{R}}^{d}$,

[TABLE]

Our next main result is the following which provides a sub-Gaussian upper tail estimate on the limit $\mathfrak{h}$ defined in (1.11).

Theorem 1.3**.**

Let $d\geq 3$ and $\beta\in(0,\beta_{L^{2}})$. Then for any $\theta>0$, there exists a constant $C\in(0,\infty)$ so that

[TABLE]

In particular, $\mathfrak{h}\in L^{p}(\mathbb{P})$ for any $p\in\mathbb{R}$.

From Theorems 1.1 and 1.3, we derive

Corollary 1.4**.**

In the hypothesis of Corollary 1.2, we have for any $x_{0}\in{\mathbb{R}}^{d}$,

[TABLE]

1.3 Literature remarks and discussion.

In the present set up, by finding a non-trivial limit when letting the regularization parameter vanish we have obtained a non-trivial renormalization of KPZ equation (1.1). Let us stress the main specificity of Theorem 1.1. The approximating sequence $({\mathfrak{h}}(\xi^{(\varepsilon,t,x)});\varepsilon>0)$ in the case of flat initial condition combines two interesting properties:

• •

it is constant in law for all $(\varepsilon,t,x)$, with law given by the one of $\log\mathscr{Z}_{\infty}$;

• •

it approximates $h_{\varepsilon}(t,x)$ in probability.

(Similar properties hold for the other initial conditions). Since it depends on $\varepsilon$, it is not a (strong) limit, but it can be used similarly. In particular, fluctuations can be studied as shown in [8]. This is quite different from using a deterministic centering, e.g., $\tilde{h}_{\varepsilon}(t,x)=h_{\varepsilon}(t,x)-\mathbb{E}h_{\varepsilon}(t,x)$. As mentioned in [10], $\tilde{h}_{\varepsilon}$ does not converge to 0 pointwise, but it does as a distribution. Integrating $\tilde{h}_{\varepsilon}$ in space against test functions cause oscillations to cancel. On the contrary, in our result $h_{\varepsilon}(t,x)-{\mathfrak{h}}(\xi^{(\varepsilon,t,x)})\to 0$ pointwise, and we do not need any averaging in space.

We also emphasize that our results concern studying the asymptotic behavior of the solution to the non-linear equation (1.2), and are not restricted to the linear multiplicative noise stochastic heat equation (see (1.6)). Furthermore, the statements of the results concern the solution itself, without need of integrating spatially against test functions. However, note that the limit obtained in Theorem 1.1 does depend on the smoothing procedure $\phi$ as well as on the disorder parameter $\beta$ and it is not universal (in particular, for $\beta<\beta_{L^{2}}$, the variance of $\exp(\mathfrak{h})$ can be computed from the RHS of (1.14) for $x=0$, and it depends on the mollification). Thus the present scenario lies in total contrast with the $1$-dimensional spatial case where the limit can be defined by a chaos expansion [1, 2] (with the parameter $\beta$ absorbed by scaling) or via the theory of regularity structures ([14]) which also produces a renormalized limit which does not depend on the mollification scheme.

In [8] we have also investigated the rate of the convergence of $h_{\varepsilon}\to\mathfrak{h}$ for small enough $\beta$ and showed that $\varepsilon^{\frac{d-2}{2}}[h_{\varepsilon}(t,x)-\mathfrak{h}(\xi^{{\scriptscriptstyle{({\varepsilon,t,x}})}})]\stackrel{{\scriptstyle\rm{law}}}{{\longrightarrow}}N(0,\sigma^{2}(\beta))$ for each fixed $x\in{\mathbb{R}}^{d}$ and $t>0$. For larger $\beta$, the so-called KPZ regime is expected to take place with different limits, different scaling exponents and non-Gaussian limiting distributions. In particular, the variance in the above Gaussian distribution is given by

[TABLE]

which already diverges for $\beta>\beta_{L^{2}}$ indicating that the amplitude of the fluctuations, or at least their distributional nature, changes at this point. However the KPZ regime is not expected before the critical value $\beta_{c}$. Hence this region $\beta\in(\beta_{L^{2}},\beta_{c})$ remains mysterious.

Finally, we remark on the correlation structure of the limit $\mathfrak{u}$ which were computed in [8]. It was shown that, for $\beta$ small enough,

[TABLE]

with ${\mathfrak{C}}_{1}=E_{{\bf e_{1}}/\sqrt{2}}[\exp\{\beta^{2}\int_{0}^{\infty}V(\sqrt{2}W_{s})ds\}-1]$. The above correlation structure also underlines that solution $u_{\varepsilon}(t,x)$ and $u_{\varepsilon}(t,y)$ becomes asymptotically independent so that the so that the spatial averages $\int_{{\mathbb{R}}^{d}}f(x)\,u_{\varepsilon}(t,x)\,\text{\rm d}x\to\int f(x)\overline{u}(t,x)\,\text{\rm d}x$ become deterministic and $\overline{u}$ solves the unperturbed heat equation $\partial_{t}\overline{u}=\frac{1}{2}\Delta\overline{u}$. As remarked earlier, the spatially averaged fluctuations $\varepsilon^{1-\frac{d}{2}}\int_{{\mathbb{R}}^{d}}f(x)[u_{\varepsilon}(t,x)-\overline{u}(t,x)]\,\text{\rm d}x$ were shown to converge ([18, 12, 10]) to the averages of the heat equation with additive space-time white noise with variance given by (a constant multiple of) $\sigma^{2}(\beta)$, which also underlines the Edwards-Wilkinson regime in weak disorder. For averaged fluctuations of similar nature in $d=2$ we refer to [5, 4, 11].

2 Proof of Theorem 1.1

We now consider the regularized KPZ equation (1.6) as before, but with different initial data and identify the limit of the solution up to leading order. For notational brevity, we will write

[TABLE]

where $V=\phi\star\phi$ so that $\mathscr{Z}_{T}(x)=\mathscr{Z}_{T}(\xi;x)=E_{x}[\Phi_{T}]$ and $\mathbb{E}[\mathscr{Z}_{T}]=1$.

We also remind the reader that $u_{\varepsilon}$ solves (1.6) and $h_{\varepsilon}=\exp[u_{\varepsilon}]$ solves (1.2) with $C_{\varepsilon}=\frac{\beta^{2}\varepsilon^{-2}}{2}V(0)$. Finally, recall that u_{\varepsilon}(t,x)=\mathscr{Z}_{\frac{t}{\varepsilon^{2}}}\big{(}\xi^{{\scriptscriptstyle{({\varepsilon,t}})}};\frac{x}{\varepsilon}\big{)} with $\xi^{{\scriptscriptstyle{({\varepsilon,t}})}}$ given by (1.3) and (1.4).

2.1 General initial condition: Proof of (1.12).

Fix continuous functions $u_{0}:{\mathbb{R}}^{d}\to(0,+\infty)$ and $h_{0}:{\mathbb{R}}^{d}\to{\mathbb{R}}$ which are bounded from above, consider the solution of SHE

[TABLE]

or, equivalently by the relations $u_{\varepsilon}=\exp h_{\varepsilon}$ and $u_{0}=\exp h_{0}$, the solution of KPZ

[TABLE]

As before, we have the Feynman-Kac representation

[TABLE]

with $\xi^{(\varepsilon,t)}$ as above.

Lemma 2.1**.**

For $\beta\in(0,\beta_{L^{2}})$,

[TABLE]

where $\theta_{x}$ denotes the canonical spatial translation in the path space $\mathscr{C}^{\alpha}$ of the white noise and $\overline{u}$ solves $\partial_{t}\overline{u}=\frac{1}{2}\Delta\overline{u}$ with $\overline{u}(0,\cdot)=u_{0}(\cdot)$.

Proof.

Note that

[TABLE]

Then

[TABLE]

Furthermore,

[TABLE]

where the right hand side is a triplet of three independent random variables, with $Z_{t}^{(1)}$ and $Z_{t}^{(2)}$ distributed as $W_{t}$. Hence, expectation (2.5) vanishes as $\varepsilon\to 0$, provided that $u_{0}$ is bounded and continuous, and because of uniform integrability which is implied by

[TABLE]

for $\beta<\beta_{L^{2}}$ and $\delta>0$ small enough. The proof is concluded by the observation that $E_{x/\varepsilon}\left[\Phi_{\varepsilon^{-2}t}(\xi;W)\right]-{\mathfrak{u}}(\xi\circ\theta_{x/\varepsilon})\stackrel{{\scriptstyle L^{2}}}{{\to}}0$.

We now end the

**Proof of (1.12): ** For $\beta<\beta_{L^{2}}$, for all $t,x$, as $\varepsilon\to 0$, we first show

[TABLE]

Note that (2.7) follows directly from Lemma 2.1 and (2.4). Then, since $\mathfrak{u}>0$, taking logarithm we deduce the convergence in probability (1.12).

Remark 2.2**.**

Recall that in [20] it was shown that, for any smooth function $f$ with compact support, $\int_{{\mathbb{R}}^{d}}u_{\varepsilon}(t,x)f(x)\,\,\text{\rm d}x\to\int_{{\mathbb{R}}^{d}}\overline{u}(t,x)f(x)\,\,\text{\rm d}x$. Note that unlike the latter statement, no smoothing in space is needed in the present context. In fact, we can recover the spatially averaged statement from above. Indeed, fast decorrelation in space of $\xi^{(\varepsilon,t,x)}$ as $\varepsilon\to 0$, ergodicity and smoothness justify the equivalence below:

[TABLE]

2.2 Narrow-wedge initial condition: Proof of (1.13).

Fix $x_{0}\in{\mathbb{R}}^{d}$, and consider the solution of SHE

[TABLE]

or, equivalently by the relation $u_{\varepsilon}=\exp h_{\varepsilon}$, the solution of KPZ

[TABLE]

By Feynman-Kac formula, the solution of SHE now admits a Brownian bridge representation:

[TABLE]

where $E_{0,x}^{t,y}$ denotes expectation with respect to a Brownian bridge starting at $x$ and conditioned to be found at $y$ at time $t$, $\rho$ is the $d$-dimensional Gaussian kernel and

[TABLE]

so that we again have $\xi_{(\varepsilon,x_{0})}\stackrel{{\scriptstyle\rm{law}}}{{=}}\xi$.

The following Lemma follows the approach for proving the local limit theorem as in [23, 25].

Lemma 2.3**.**

For $\beta\in(0,\beta_{L^{2}})$, for any $A>0$,

[TABLE]

Proof.

We will write $X=\varepsilon^{-1}x$, $T=\varepsilon^{-2}t$ and let $m=m_{\varepsilon}$ be a time parameter, such that $m_{\varepsilon}\to\infty$ and $m_{\varepsilon}=o({T})$, as $\varepsilon\to 0$. We use the notation:

[TABLE]

Step 1:

We first want to approximate $E_{0,0}^{T,X}[\Phi_{T}]$ by $E_{0,0}^{T,X}\left[\Phi_{m}\Phi_{T-m,T}\right]$ in $L^{2}$-norm, so we compute the difference:

[TABLE]

To show that the right hand side which goes to [math] as $\varepsilon\to 0$, it suffices to observe that, for all $a>0$,

[TABLE]

To prove this, we use Hölder’s inequality similarly to (2.6), and apply [8, Lemma 3.5] (alternatively [25, Corollary 3.8]) and transience of Brownian motion for $d\geq 3$, which implies, since $m\to\infty$, that

[TABLE]

Step 2:

We wish to use Markov property and symmetry of the Brownian bridge to show that $E_{0,0}^{T,X}\left[\Phi_{m}\Phi_{T-m,T}\right]$ factorizes asymptotically into the product $E_{0}\left[\Phi_{m}\right]E_{0}[\Phi_{m}(\xi^{(1,T,X)})]$, which satisfies:

[TABLE]

as $\varepsilon\to 0$ for $\beta<\beta_{L^{2}}$, by Cauchy-Schwarz inequality and invariance in law of the white noise with a shift by $X$. Hence, we compute

[TABLE]

After change of variable by setting $Y=\sqrt{T}y$ in the above integral and since $\rho(Ts,\sqrt{T}z)=T^{-d/2}\rho(s,z)$, observe that by Jensen’s inequality and dominated convergence, we can prove that

[TABLE]

if we can show that for all fixed $y\in\mathbb{R}^{d}$,

[TABLE]

To prove this, we use the density of the Brownian bridge, at truncated time horizon, with respect to the Brownian motion ([8, Lemma 3.4]), to get that:

[TABLE]

After rescaling, the difference inside the parenthesis goes almost surely to [math], for $y$ fixed and uniformly in $|x|\leq A$; we conclude the proof of the lemma using Hölder’s inequality.

We can now conclude the

Proof of (1.13): Let $h_{\varepsilon}$ be the narrow-wedge height function solution of of (2.9). We need to show that for $\beta<\beta_{L^{2}}$, for all $t,x$, as $\varepsilon\to 0$,

[TABLE]

with $\xi_{(\varepsilon,x_{0})}$ in (2.12). We use the representation (2.11) and the property $\xi_{(\varepsilon,x_{0})}\stackrel{{\scriptstyle\rm{law}}}{{=}}\xi$, so that we can exchange $\xi$ with $\xi_{(\varepsilon,x_{0})}$, in convergence (2.13) taken with endpoint $\varepsilon^{-1}(x-x_{0})$. This leads to the above convergence in probability for the logarithm, proving (1.13).

3 Proof of Theorem 1.3.

We focus on showing that for $\beta<\beta_{L^{2}}$, $\log{\mathscr{Z}}_{\infty}$ admits a sub-Gaussian lower tails estimate, that is, for some $C\in(0,\infty)$ and any $\theta>0$,

[TABLE]

We invoke a second moment method combined with the Talagrand’s concentration inequality as in [6] (see also [19, Section 2.2]). As the Cauchy-Schwarz inequality, which is a central tool in this proof is not directly available in the continuous setting, we choose to introduce a discretization of the white noise to recover it.

Consider $\mathcal{R}_{n}$ a tiling of $[0,2^{n}]\times[-2^{n},2^{n}]^{d}$, composed of cubes of length $2^{-n}$, such that every cube of $\mathcal{R}_{n}$ can be divided in $2^{d+1}$ cubes of $\mathcal{R}_{n+1}$. We define a discrete version of ${\mathscr{Z}}_{T}$ through:

[TABLE]

where,

• •

$\phi_{W}(s,y):=\phi(W_{s}-y)$,

• •

$\phi_{W}^{{\scriptscriptstyle{({n}})}}(s,y)=\inf_{R}\phi_{W}$ if $(s,y)$ is in a cube $R$ of $\mathcal{R}_{n}$, and [math] otherwise.

We stress that $\phi_{W}^{{\scriptscriptstyle{({n}})}}$ is non-decreasing with $n$ and converges almost surely to $\phi_{W}$.

Using the Gaussian covariance structure, we have that

[TABLE]

For $\beta<\beta_{L^{2}}$ and since $\phi_{W}^{{\scriptscriptstyle{({n}})}}\leq\phi_{W}$, we immediatly obtain from the monotone convergence theorem that the right-hand side goes to zero in the limit $T\to\infty$ followed by $n\to\infty$. By Doob’s $L^{2}$ inequality applied to the martingale ${\mathscr{Z}}-{\mathscr{Z}}^{{\scriptscriptstyle{({n}})}}$, this implies in particular that

[TABLE]

Since $\phi_{W}^{{\scriptscriptstyle{({n}})}}$ is set to be [math] outside of $[0,2^{n}]\times[-2^{n},2^{n}]^{d}$, the set $\mathcal{C}_{n}$ containing the centers of the cubes of $\mathcal{R}_{n}$ is finite. Hence, as each cube has volume $2^{-(d+1)n}$, we can write that

[TABLE]

where the $\xi_{n}(i,x)$ are independent centered Gaussian random variables of variance $1$. Then, we define the polymer measure $\widehat{P}_{\beta,\xi_{n}}$ of renormalized partition function ${\mathscr{Z}}_{\infty}^{{\scriptscriptstyle{({n}})}}$:

[TABLE]

where we have set, using (3.3),

[TABLE]

Finally, we let $\widehat{E}_{\beta,\xi_{n}}$ denote expectation corresponding to $\widehat{P}_{\beta,\xi_{n}}$.

Now, we can compare the free energies of two realizations of the noise $\xi_{n}(i,x)$ and $\xi^{\prime}_{n}(i,x)$:

[TABLE]

where $\mathrm{d}(\cdot,\cdot)$ denotes the euclidean distance on $\mathbb{R}^{\mathrm{Card}(\mathcal{C}_{n})}$, and where we used Jensen’s and Cauchy-Schwarz inequalities for respectively the first and second lower bounds.

Let $m$ and $C$ be two positive constants, and consider the set:

[TABLE]

For $\xi^{\prime}_{n}\in\mathcal{A}_{n}$, the above computation implies that

[TABLE]

therefore, assuming

[TABLE]

property (3.1) results from (3.4), (3.2) and the following Gaussian concentration inequality (Lemma 4.1 in [6], extracted from [24]):

[TABLE]

We now prove (3.5). By convergence (3.2), and since ${\mathscr{Z}}_{\infty}>0$ a.s., we can find $m>0$, such that for $n$ large enough,

[TABLE]

Then, for a large enough $C$,

[TABLE]

In the above display, the first lower bound follows from the definition of $\mathcal{A}$, while we used $\mathbb{P}[A\cap B]\geq\mathbb{P}[A]-\mathbb{P}[B^{c}]$ in the second lower bound. The third inequality comes from Markov’s inequality and the upper-bound $\phi_{W}^{(n)}\leq\phi_{W}$. Positivity of the left hand-side of the third line is assured for $C$ large enough, provided that $\beta<\beta_{L^{2}}$.

This entails that the KPZ limit $|\mathfrak{h}|$ (recall (1.11)) has all positive and negative moments for all $\beta<\beta_{L^{2}}$. Indeed, letting $\log_{-}=\log\wedge\,0$, the sub-Gaussian decay of the left tail of $\log{\mathscr{Z}}_{\infty}$ (3.1) gives $\mathbb{E}[\exp\{\nu\log_{-}{\mathscr{Z}}_{\infty}\}]<\infty$, for all $\nu\in\mathbb{R}$. Moreover, by definition of the $L^{2}$ region, we have $\mathbb{E}[\exp\{2\log{\mathscr{Z}}_{\infty}\}]<\infty$. Hence, $\log{\mathscr{Z}}_{\infty}$ admits all positive and negative moments.

Remark 3.1**.**

After finishing the writing of the present article, we learned that another proof of the negative moments of $\mathscr{Z}_{\infty}$ has been recently proposed in [15, 10] using a continuous approximation of the white noise. A proof of the corresponding result for the KPZ equation in dimension 2, which relies on the convexity of the free energy and the Malliavin derivative, can be found in [4].

Appendix A Appendix.

We will include some elementary facts regarding the regularity properties of space-time white noise $\xi$. For any $z,z^{\prime}\in{\mathbb{R}}\times{\mathbb{R}}^{d}$, we will denote by $\|\cdot\|$ the parabolic distance given by $\|z-z^{\prime}\|=|t-t^{\prime}|^{1/2}+\sum_{i=1}^{d}|x_{j}-x_{j}^{\prime}|$, where $z=(t,x)$ and $z^{\prime}=(t^{\prime},x^{\prime})$. Recall that the Hölder space of positive exponent $\alpha\in(0,1)$ consists of all functions $u:{\mathbb{R}}\times{\mathbb{R}}^{d}\to{\mathbb{R}}$ such that for any compact set $K\subset{\mathbb{R}}\times{\mathbb{R}}^{d}$,

[TABLE]

The corresponding Hölder (Besov) space of negative regularity is defined as follows. First for any $k\in\mathbb{N}$, let $B_{k}$ denote the space of all smooth functions $\varphi\colon{\mathbb{R}}\times{\mathbb{R}}^{d}\to{\mathbb{R}}$ which are supported on the unit ball in $({\mathbb{R}}\times{\mathbb{R}}^{d},\|\cdot\|)$ such that

[TABLE]

Then for any fixed $\alpha<0$, we define the space $\mathscr{C}^{\alpha}$ to be the space of all tempered distributions $\eta\in\mathcal{S}^{\prime}({\mathbb{R}}\times{\mathbb{R}}^{d})$ such that for any compact set $K\subset{\mathbb{R}}\times{\mathbb{R}}^{d}$,

[TABLE]

where $k=\lceil-\alpha\rceil$ and

[TABLE]

A crucial estimate on the Besov norm $\|\cdot\|_{\mathscr{C}^{\alpha}(K)}$ is given by

[TABLE]

where $\widetilde{K}$ is also a compact set slightly larger than $K$ and $u$ is a single, well-chosen test function (which can be constructed by wavelets, see [14]). Recall that if $\xi$ is space-time white noise (i.e, $\mathbb{E}[\langle\xi,\varphi_{1}\rangle\langle\xi,\varphi_{2}\rangle]=\int\varphi_{1}(t,x)\varphi_{2}(t,x)\,\,\text{\rm d}t\,\text{\rm d}x$), then with

[TABLE]

The following result, which is a consequence of Kolmogorov’s lemma and (A.1), then implies the desired regularity property of $\xi$.

Lemma A.1**.**

Fix $\alpha<0$ and $p\geq 1$ and let $\eta$ be a linear map from $\mathcal{S}({\mathbb{R}}\times{\mathbb{R}}^{d})$ to the space of random variables. Suppose there exists $C\in(0,\infty)$ such that for all $z\in{\mathbb{R}}\times{\mathbb{R}}^{d}$ and all $u\in\mathcal{S}({\mathbb{R}}\times{\mathbb{R}}^{d})$ with compact support in ${\mathbb{R}}\times{\mathbb{R}}^{d}$ with $\sup_{w}|u(w)|\leq 1$ one has

[TABLE]

Then there exists a random distribution $\widetilde{\eta}$ in $\mathcal{S}({\mathbb{R}}\times{\mathbb{R}}^{d})$ such that for all $\alpha^{\prime}<\alpha-\frac{(d+2)}{p}$ and compact set $K$,

[TABLE]

Since for any $p\geq 1$, $\mathbb{E}[\langle\xi_{\lambda},\varphi\rangle^{p}]\leq C_{p}\mathbb{E}[\langle\xi_{\lambda},\varphi\rangle^{2}]^{p/2}$, then Lemma A.1 and (A.2) imply that $\xi$ has regularity $\mathscr{C}^{-\frac{d}{2}-1-\delta}$ for any $\delta>0$.

Compliance with Ethical Standards: The authors declare that they have no conflict of interest.