The scaling limit of the KPZ equation in space dimension 3 and higher
Jacques Magnen and Jérémie Unterberger
Centre de Physique Théorique,^{1}^{1}1Laboratoire associé au CNRS UMR 7644 Ecole Polytechnique,
91128 Palaiseau Cedex, France
jacques.magnen@cpht.polytechnique.fr
Institut Elie Cartan,^{2}^{2}2Laboratoire associé au CNRS UMR 7502. J. Unterberger acknowledges the support of the ANR, via the ANR project ANR16CE40002001. Université de Lorraine,
B.P. 239, F – 54506 VandœuvrelèsNancy Cedex, France
jeremie.unterberger@univlorraine.fr
We study in the present article the KardarParisiZhang (KPZ) equation
in dimensions in the perturbative regime, i.e. for small enough and a smooth, bounded, integrable initial condition . The forcing term in the righthand side is a regularized spacetime white noise. The exponential of – its socalled ColeHopf transform – is known to satisfy a linear PDE with multiplicative noise. We prove a largescale diffusive limit for the solution, in particular a timeintegrated heatkernel behavior for the covariance in a parabolic scaling.
The proof is based on a rigorous implementation of K. Wilson’s renormalization group scheme. A double cluster/momentumdecoupling expansion allows for perturbative estimates of the bare resolvent of the ColeHopf linear PDE in the smallfield region where the noise is not too large, following the broad lines of IagolnitzerMagnen [42]. Standard large deviation estimates for make it possible to extend the above estimates to the largefield region. Finally, we show, by resumming all the byproducts of the expansion, that the solution may be written in the largescale limit (after a suitable Galilei transformation) as a small perturbation of the solution of the underlying linear EdwardsWilkinson model () with renormalized coefficients .
Keywords: KPZ equation, ColeHopf transformation, directed polymer, constructive field theory, renormalization, cluster expansion, resolvent, large deviation estimates.
Mathematics Subject Classification (2010): 35B50, 35B51, 35D40, 35K55, 35R60, 35Q82, 60H15, 81T08, 81T16, 81T18, 82C41.
Contents
0 Introduction
The KPZ equation [50] is a stochastic partial differential equation describing the growth by normal deposition of an interface in space dimensions, see e.g. [7, 19]. By definition the time evolution of the height , , is given by
(0.1) 
where is a regularized white noise, and are constant. Three terms contribute to eq. (0.1): a viscous term proportional to the viscosity , leading to a smoothening of the interface; a growth by normal deposition with rate , called deposition rate, and playing the rôle of a coupling constant; and a random rise or lowering of the interface modelling molecular diffusivity, with coefficient called noise strength. In a related context, also represents the free energy of directed polymers in a random environment [43, 20, 26]. It makes sense to consider more general nonlinearities of the form with , say, positive and convex, instead of , which is in any case an approximation of , assuming that the gradient (the slope of the interface) remains throughout small enough so that the evolution makes physically sense, precluding e.g overhangs.
The interest is here in the largescale limit of this equation, for and/or large. A wellknown naive rescaling argument gives some ideas about the dependence on the dimension of this limit. Namely, the linearized equation, a stochastic heat (or infinitedimensional OrnsteinUhlenbeck [62]) equation called EdwardsWilkinson model [7] in the physics literature,
(0.2) 
– where requires no regularization – is invariant under the rescaling ; we used here the equality in distribution, . Assuming that is a solution of the KPZ equation instead yields after rescaling
(0.3) 
where (up to change of regularization) . For , vanishes in the limit ; in other terms, the KPZ equation is infrared superrenormalizable, hence (powerlike) asymptotically free at large scales in dimensions, i.e. expected to behave, in a small coupling (also called small disorder) regime where , like the corresponding linearized equation up to a redefinition (called renormalization) of the diffusion constant and of the noise strength .
Let us emphasize the striking difference with the onedimensional equation. For this equation, scaling behaviors, see (0.3), are reversed with respect to , in other words, KPZ is (powerlike) asymptotically free at small scales (i.e. in the ultraviolet), or equivalently (in the PDE analysts’ terminology) subcritical. A large part of the interest for this equation comes from the fact that the largescale strongly coupled theory [3, 26] is understood by comparison with integrable discrete statistical physics models [27, 64, 68, 71] relating to weakly asymmetric exclusion process [11] or the TracyWidom distribution of the largest eigenvalue of random matrices connected with Bethe Ansatz [71], free fermions and determinantal processes [44],… Note that is believed by perturbative QFT arguments to be strongly coupled at large scales [7, 19] and its largescale limit is not at all understood.
We prove the diffusive limit of dimensional KPZ with small coupling in the present work, thus establishing on firm mathematical ground old predictions of physicists, see e.g. Cardy [19]. The space dimension does not really matter as long as . In the smallcoupling regime, contrary to the case, we fall into the EdwardsWilkinson universality class.
**************************************************************************
In comparison with the achievements made in the study of strongly coupled largescale , this problem looks at first sight of lesser importance and difficulty. We believe that the interest of our result lies in the precision of our asymptotics, and in the potential wide scope of applicability of our methods.
Namely, the KPZ model is one particular instance of a large variety of dynamical problems in statistical physics, modelized as interacting particle systems, or as parabolic SPDEs heuristically derived by some mesoscopic limit, which have been turned into a functional integral form analogous to the Gibbs measure of equilibrium statistical mechanics, , using the socalled response field (RF), or MartinSiggiaRose (MSR) formalism and studied by using standard perturbative expansions originated from quantum field theory (QFT); for reviews see e.g. [19] or [2]. Despite the lack of mathematical rigor, this formalism yields a correct description of the qualitative behaviour of such dynamical problems in the large scale limit.
The Feynman perturbative approach, see e.g. [55], consists in expanding into a series in and making a clever resummation of some truncation of it into socalled counterterms, represented in terms of a sum of diagrams; as such, it is nonrigorous, since it yields point functions in terms of an asymptotic expansion in the coupling parameter which is divergent in all interesting cases (at least for bosonic theories). A few years ago, however, Gubinelli, M. Hairer, H. Weber,… [38, 39, 40, 6, 17, 18, 22, 23, 41, 59, 21], drawing sometimes on a dynamical approach to the construction of equilibrium measures advocated by Nelson [61], ParisiWu [63], and JonaLasinio, Mitter and Sénéor [45, 46, 47], have started developing this philosophy in a systematic way to solve subcritical parabolic SPDEs rigorously, i.e. beyond perturbation theory. Such SPDEs have only a finite number of counterterms, each counterterm being the sum of a finite number of terms (that can be interpreted in terms of Feynman diagrams), which makes the task considerably easier, but still far from trivial.
**************************************************************************
Constructive approaches developed in the context of statistical physics by mathematical physicists from the mid60es, see e.g. [28, 32, 33, 30, 31, 36, 42, 56, 57, 60] and surveys [35, 58, 66, 67, 72], have developed sophisticated, systematic truncation methods making it possible to control the error terms. The partial resummations are interpreted in the manner of K. Wilson [75, 76] as a scalebyscale, finite renormalization of the parameters of the Lagrangian . In many instances it has proved possible to subtract scale counterterms explicitly by hand and prove that the remainder is finite, yielding some description of the effective, largescale theory, see e.g. works in diverse contexts – random walks in random environment, KAM theory, etc. – by Bricmont, Gawedzki, Kupiainen and coauthors [15, 14, 16], and recent extensions to the study of subcritical parabolic PDEs [51, 52], as an alternative to the ”global counterterm” strategy mentioned in the last paragraph. However, the implementation of a fullfledged, multiscale constructive scheme is for the moment limited to equilibrium statistical physics models.
The present work is, to the best of our knowledge, the first attempt to use such a scheme in the context of nonequilibrium statistical mechanics, here for a parabolic SPDE. Instead of using the MSR formalism, we develop (as all previously mentioned mathematically rigorous approaches do) a more straightforward approach, starting directly from the equation and cutting the propagator into scales. We actually work on the following model.
**************************************************************************
The model. Let . We consider the following equation on ,
(0.4) 
where is a white noise regularized in time and in space; is a smooth, bounded, integrable initial condition, i.e. are ; is small enough; and is a constant, average interface velocity which we shall fix later on.
The precise choice of regularization for the white noise is unimportant; one should just keep in mind that local (in time and space) solvability of (0.1) in a strong sense requires that, for every compact set (equivalently, for any as in Definition 2.1 (iii)), is locally integrable. For simplicity of exposition, we define to be a smooth, stationary Gaussian noise with shortrange covariance. To be definite:
we fix a smooth, isotropic (i.e. invariant under space rotations) function with support and norm , and let
(0.5) 
Our main result is the following. Gaussian expectation with respect to is denoted either by , or or also if one wants to emphasize the dependence on the parameters ; the result also depends obviously on the initial condition . By convention, refers to the expectation with respect to the measure of the EdwardsWilkinson equation with zero initial condition, where is a standard (unregularized) spacetime white noise; for this equation we implicitly set . By definition, is a centered Gaussian process.
Theorem 0.1
(Main Theorem).
Let . Fix and a smooth, bounded, integrable initial condition . Let be small enough, . Then there exist three coefficients , and , all independent of
the initial condition ,
such that the solution of the KPZ equation (0.4) satisfies the following asymptotic properties:

for all with ,
(0.6) 
for all , with , and , letting ,
(0.7)
Since is a Gaussian measure, 2. may be rephrased as follows. Let
(0.8) 
(, ). Then
(0.9) 
and
(0.10) 
where the sum ranges over all pairings of the indices .
In other words, up to a Galilei transformation , the point functions of the KPZ equation behave asymptotically in the largescale limit as the point functions of the solution of the EdwardsWilkinson equation with renormalized coefficients ,
(0.11) 
where requires no regularization. Generally speaking, main corrections to the above asymptotic behaviour (0.6,0.10) are smaller by as proved in §5.3 D. Effective coefficients have a (diverging) asymptotic expansion in terms of ; lowestorder corrections in are computed in (4.29) and (LABEL:eq:leadingDeff). The term in (0.6) is a contribution due to the initial condition; further contributions of the initial condition to point functions come with an extra multiplicative factor in , which is the scaling of the vertex. Corrections to Gaussianity of point functions, of order , are examined in (2) a few pages below. Furthermore, our multiscale scheme actually involves an effective propagator differing slightly from the effective EdwardsWilkinson propagator , see section 7; this implies a correction w.r. to the r.h.s. of (0.10) with a small extra prefactor, which is proved to be a but could easily be improved to with arbitrary large.
Remark. A more common choice of regularization for is to take a discretized ”kick force”, namely, we pave by unit size intervals , , and let , be independent, centered Gaussian fields on which are constant in time and have smooth, spacetranslation invariant covariance kernel with finite range, for instance. This does not change the conclusion of Theorem 0.1, except that, the law of being now only periodic in time, is now a periodic function instead of the constant . This regularization has several advantages (see section 1); it allows in particular an explicit representation of in probabilistic terms. The scheme of proof extends without any significant modification if the covariance kernel decreases heatkernellike in space, e.g. if where is a standard space white noise, and is some constant.
Furthermore, it follows from the proof (see section 5) that the value of may be obtained by equating it to the constant such that independently of , in coherence with the value obtained in CarmonaHu [20] in a discrete setting for a random directed polymer measure (see section 2.1), where is the ColeHopf transform of (see below). Let us note that the equality between and points out to the fact that we are in a weak disorder regime in which the annealed and quenched free energies coincide. However, our proof is independent of that of Carmona and Hu (see [20], Theorem 1.5), based on Gaussian concentration inequalities.
**************************************************************************
The proof follows closely the article by IagolnitzerMagnen [42] on weakly selfavoiding polymers in four dimensions, which is the main reference for the present work. Namely, up to the change of function (called ColeHopf transform) and of coupling constant, , the KPZ equation is equivalent to the linear equation , solved as , where is a random resolvent. Formally then, our problem is a parabolic counterpart to the largescale analysis of polymers in a weak random potential solved in [42] by studying the equilibrium resolvent , where the ””coefficient is the Edwards model representation of the selfavoiding condition (the model is solved for but the selfavoiding condition is recovered for ). Though the two models are physically unrelated, one must analyze similar mathematical objects. As is often the case, the model with a time evolution (i.e. the parabolic one) turns out to be easier than the equilibrium model (i.e. the elliptic one), because of the causality constraint.
The general scheme of proof, following, as mentioned above, the philosophy of constructive field theory, is to introduce a multiscale expansion and define a renormalization mapping, , or equivalently (later on interpreted as the flow of the coupling constant through the ColeHopf transform), ensuring the convergence of the expansion at each scale and allowing to control error terms. The average interface velocity is fixed by requiring that the asymptotic velocity vanishes. The original parameters , called bare parameters, describe the theory at scale , while the EdwardsWilkinson model with scale parameters and drift velocity give a good approximation of the theory at time distances of order , which becomes asymptotically exact in the infrared limit, when . This goal is achieved in general by using a phasespace expansion, i.e. a horizontal cluster expansion casting into the form of a series the interactions at a given energymomentum level between the degrees of freedom, and a vertical cluster or momentumdecoupling expansion separating the different energymomentum levels. Energy, resp. momentum, are the Fourier conjugate variables of time and space; here a given energymomentum level is adequately defined by considering heatkernel propagators
with . Then the above series (roughly speaking, a truncated power series in the coupling constants with a bounded integral, Taylorlike remainder) converge if the bare coupling constant is small enough.
With our choice of covariance function for , however, the flow of the parameters is actually trivial starting from , i.e. for , and the noise strength , defined by resumming connected diagrams with four external legs, though scaledependent, requires no renormalization at all, because the equation is infrared superrenormalizable, and the total correction (obtained by summing over scales) is finite. This, and also the causality condition preventing the socalled lowmomentum field accumulation problem [36, 30, 72], leads to a much simplified framework, from which the phase space analysis has almost disappeared. Only scale 0, twopoint diagrams need to be renormalized, with a contribution at near zero momentum
leaving a remainder of parabolic order three in the momenta, i.e. or . Scale 0 diagrams are connected by ”lowmomentum” heatkernel propagators with , . A crucial point in the proof is that, thanks to the , remainders integrated over spacetime cost a factor , namely (see (2.19) and (5.21))
or, simply said, . What is left of the cluster expansions is adequately resummed as in [42] into the random resolvent in the form of localized ”vertex insertions” (see section 5), thereby suppressing combinatorial factors which make the series divergent. Then the contribution of all vertex insertions is bounded by some contour integral of a modified resolvent through the use of Cauchy’s formula.
An extra complication comes however from the inverse ColeHopf transform. Applying cluster expansions – which is done in practice by differentiation with respect to some additional parameters – to leads to rational expressions of the form , where the ’s are differential operators, acting on ”replicas” of . Then the scale 0 diagrams requiring renormalization can be factorized, hence averaged with respect to the measure . Remaining terms are shown to yield a convergent series in the form of a sum over ”polymers” for small enough.
The and prefactors contained in Theorem 0.1 may be guessed from the following guiding principles, put into light by the cluster expansion.
(1) First, the twopoint function of the renormalized EdwardsWilkinson equation,
(0.13) 
(), scales like , as can be seen by simply rescaling variables in the integral. There are two regimes: the equilibrium regime (), in which is essentially the equilibrium Green function of the Laplacian; the dynamical regime (), in which .
(2) The connected quantities (also called truncated point functions) are . Namely, Gaussian pairwise contractions yield the expected scaling in , i.e. per link, as expected from (1); whereas the connected expectation requires supplementary links and twice as much vertices (since these are not present in the linear theory) in the expansion, contributing an extra small prefactor. The cluster expansion makes it possible to develop those links explicitly.
The plan of the article is as follows. We start by recalling the ColeHopf transform in section 1, and make the bridge to previous results on the subject stated in terms of the associated directed polymer measure. We then introduce in section 2 a multiscale expansion for the propagators, together with multiscale estimates (also called ”powercounting”), which are the building blocks of our approach. Sections 3, 4, and 5 are the heart of the article. The dressed equation, and the cluster expansion thereof, is presented in section 3. Section 4 is dedicated to renormalization; the scale 0 counterterms obtained by factorizing twopoint functions through a supplementary Mayer expansion are bounded. Then we show in section 5 how to bound the sum of all terms produced by the expansion, and obtain final bounds for point functions, proving thus our main result, Theorem 0.1. Finally, there are two appendices. In the first one, we provide detailed combinatorial formulas for the horizontal and Mayer cluster expansions. The second one is merely dedicated to a technical result. Pictures are provided, which are there to help the reader visualize the outcome of the various expansions.
Notations.

(parabolic distance) Let . Similarly, for , , (Hausdorff distance). Then is the space projection of the distance , i.e. , etc.

Let be two functions on some set . We write if there exists some inessential constant (possibly depending on the parameters and on the space dimension ), uniform in for small enough, such that . Then, by definition, . If and , we write .

In many situations, one obtains dependent functions such that decays Gaussianlike, for some positive constant bounded away from . We then write without further specifying the value of , which may change from line to line. For instance, if is the heat kernel, then we may write , leaving out the dependence in the parameter as explained in 2. Note however that, if , , whereas the inequality does not hold uniformly in because the space decay of is slower than that of .
Acknowledgements. We wish to thank H. Spohn, F. Toninelli and the referee for numerous discussions, suggestions and corrections, which have hopefully contributed in particular to the readability of the paper.
1 ColeHopf transform
It is wellknown that is a solution of the linear equation with multiplicative noise,
(1.1) 
where
(1.2) 
plays the rôle of a bare coupling constant, from which (representing the solution as a Wiener integral by FeynmanKac’s formula)
(1.3) 
where the expectation is relative to the Wiener measure on dimensional Brownian paths issued from with normalization, i.e. , . Thus may be interpreted as the partition function of a directed polymer, see e.g. [20] and references within, but we shall not need this interpretation in the article. Note that , where is now a standard Brownian motion, from which – forgetting about the regularization and using the variable instead of –
Thus may be expanded in a series in the parameter . Similarly, , or conversely , from which
Without using the general theory developed in [73, 74], eq. (1.3) and (LABEL:eq:4.3) show that a.s. , exist and are for , say, and compactly supported. The ColeHopf solution coincides with the solution defined for more general HamiltonJacobi equations in [73, 74].
For the rest of the subsection only, we assume that is a discretized ”kick force”, i.e. are independent and constant in time, in order to compare with the existing literature. Since are independent fields, letting , where
(1.5) 
leads to for any and if , whence more generally
(1.6) 
Expanding the exponential in (1.5) and using
(1.7) 
one gets: , whence .
Let us state an easy preliminary result, adapted from Carmona and Hu [20].
Lemma 1.1
There exists some positive constant such that the solution of the KPZ equation with zero bare velocity,
(1.8) 
verifies
(1.9) 
Furthermore, .
Proof (see [20], Lemma 3.1) Let, for general forcing term,
(1.10) 
and
(1.11) 
Conditioning with respect to the terminal condition, , means that we average with respect to the law of the Brownian bridge from to (see e.g. [49]). Then, for ,
(1.12)  
where . By construction, . Hence (by concavity of the log)
(1.13) 
Taking the expectation with respect to the noise and using independence of from , together with space translation invariance, one gets the superadditive inequality,
(1.14) 
On the other hand, by convexity of exp, . Fekete’s superadditive lemma allows us to conclude to the existence of some constant verifying (1.9). This is the constant whose existence is asserted in Main Theorem (see (0.6)). Furthermore, by Jensen’s inequality, , as observed already in [20], Prop. 1.4.
As mentioned in the Introduction, Carmona and Hu [20] actually prove the existence of a limit random variable a.s.lim for the solution of the KPZ equation with velocity