alphaFeedback

WorksMathematics

A functional central limit theorem for the KPZ equation with narrow wedge initial data

AI contributions
Not recorded for this paper
View details

Abstract

In this article, we consider the Kardar-Parisi-Zhang (KPZ) equation with narrow wedge initial data. We prove a functional central limit theorem for the spatial average of the solution, which solves an open problem from Gu and Pu (2025). Our main tool is Clark-Ocone formula from Malliavin calculus, combined with an in-depth analysis of the increments of the solution. The main step is a moment bound for time increments of the normalized spatial integral. For this, we introduce a new technique, called a ``common-ray'' change of variables. Using the fact that Malliavin derivative of the solution satisfies a Chapman-Kolmogorov relation, we reduce the kernel differences at two times to spatial increments at one common time.

Introduction

Spatial averages of stochastic partial differential equations (SPDEs) are known to exhibit macroscopic Gaussian fluctuations, i.e. they satisfy a central limit theorem. For the Kardar-Parisi-Zhang (KPZ) equation with narrow wedge initial data, finite-dimensional convergence is known from [21], but temporal increments require separate control to obtain convergence of paths. The corresponding functional limit theorem has been left open in [21]. The goal of this article is to prove the needed tightness on compact time intervals bounded away from zero.

A starting point for central limit theorems for spatial averages is the work of [23] on the stochastic heat equation (SHE) with flat initial condition u(0,x)=1u(0, x) = 1 and a Lipschitz noise coefficient σ\sigma. In the linear case σ(u)=u\sigma(u)=u, this is the parabolic Anderson model (PAM):

∂tu(t,x)=12∂x2u(t,x)+u(t,x)ξ(t,x),t>0,x∈R,u(0,⋅)=1.(1)\partial_t u(t,x)=\frac{1}{2}\partial_x^2u(t,x)+u(t,x)\xi(t,x),\qquad t>0,x\in\mathbb{R},\qquad u(0,\cdot)=1. \tag*{(1)}

The mild solution of (1) is a jointly measurable, (Ft)(\mathcal{F}_t)-adapted random field satisfying

u(t,x)=1+∫0t∫Rpt−s(x−y)u(s,y)ξ(ds,dy),u(t,x)=1+\int_0^t\int_{\mathbb{R}}p_{t-s}(x-y)u(s,y)\xi(\mathrm{d}s,\mathrm{d}y),

where the integral is interpreted in the Itô sense, (Ft)t≥0(\mathcal{F}_t)_{t\geq0} is the filtration induced by ξ\xi (given by (9) below), and pt(x)p_t(x) is the heat kernel:

pt(x)=12πtexp⁡(−∣x∣22t).p_t(x)=\frac{1}{\sqrt{2\pi t}}\exp\left(-\frac{|x|^2}{2t}\right).

The authors of [23] studied the asymptotic behaviour as N→∞N\to\infty of the spatial average:

FN(t)=∫0N(u(t,x)−1)dx.F_N(t)=\int_0^N\bigl(u(t,x)-1\bigr)\mathrm{d}x.

Their study was motivated by the fact that, after proper normalization, FN(t)F_N(t) is asymptotically equivalent to FN∗(t)F_N^*(t), the spatial average of the solution uN∗(t,x)u_N^*(t,x) of a localized version of the original equation, which has the property that {uN∗(t,xi)}i≥1\{u_N^*(t,x_i)\}_{i\geq1} are independent if the xix_i’s are far from each other. The precise meaning of this localization was given in [16]. By writing the interval [0,N][0,N] as a union of alternating big and small blocks, we obtain enough spacing between the big blocks to prove a central limit theorem (CLT). The small blocks become asymptotically negligible. The normalization is given by the variance of FN(t)F_N(t). Since {u(t,x)}x∈R\{u(t,x)\}_{x\in\mathbb{R}} is strictly stationary,

ρtu(x):=Cov⁡(u(t,x),u(t,0))=∫0tp2(t−s)(x)E∣u(s,0)∣2ds\rho_t^u(x):=\operatorname{Cov}(u(t,x),u(t,0))=\int_0^t p_{2(t-s)}(x)\mathbb{E}|u(s,0)|^2\mathrm{d}s (by Itô’s isometry), and

Var⁡(FN(t))=∫0N∫0Nρtu(x−y)dxdy=2∫0Nρtu(x)(N−x)dx∼NKt,\operatorname{Var}(F_N(t))=\int_0^N\int_0^N\rho_t^u(x-y)\mathrm{d}x\mathrm{d}y=2\int_0^N\rho_t^u(x)(N-x)\mathrm{d}x\sim NK_t,

as N→∞N\to\infty, using the fact that Kt:=∫Rρtu(x)dx=∫0tE∣u(s,0)∣2ds<∞K_t:=\int_{\mathbb{R}}\rho_t^u(x)\mathrm{d}x=\int_0^t\mathbb{E}|u(s,0)|^2\mathrm{d}s<\infty since E∣u(t,0)∣2=2et/4Φ(t/2)\mathbb{E}|u(t,0)|^2=2e^{t/4}\Phi(\sqrt{t}/2), where Φ\Phi is the standard normal c.d.f. The novelty is that this CLT is proved in [23] using a combination between Stein method and Malliavin calculus, which provides a rate for the convergence in the total variation distance, called a quantitative CLT. Moreover, the authors of [23] proved a corresponding functional CLT:

1NFN(⋅)→dG(⋅)in C[0,T],\frac{1}{\sqrt{N}}F_N(\cdot)\xrightarrow{d}G(\cdot)\quad\text{in }C[0,T],

where GG is a centered Gaussian process with covariance E[G(t)G(s)]=∫0t∧sE∣u(r,0)∣2dr\mathbb{E}[G(t)G(s)]=\int_0^{t\wedge s}\mathbb{E}|u(r,0)|^2\mathrm{d}r. Here C[a,b]C[a,b] is the space of continuous functions f:[a,b]→Rf:[a,b]\to\mathbb{R} equipped with the topology of uniform convergence.

Article [23] sparked the development of a large area dedicated to similar QCLTs and FCLTs for a variety of equations and noises. We refer the reader to [3, 4, 5, 6, 7, 8, 17, 18, 19, 20, 27, 29, 30, 31, 32] for a sample of references. The intuition behind these results is the same: the noise may be correlated in space, but as long as the correlation between the big blocks decays fast enough, the resulting variables remain weakly dependent. Moreover, the same intuition applies to a function gg of the solution. In [12], the authors proved that

1Nd/2∫[0,N]d(g(u(t,x))−E[g(u(t,x))]) dx→dN(0,ct)\frac{1}{N^{d/2}}\int_{[0,N]^d}\left(g(u(t,x))-\mathbb{E}[g(u(t,x))]\right)\,dx \xrightarrow{d} N(0,c_t)

where uu is the solution of (SHE) with flat initial condition, driven by a spatially homogeneous Gaussian noise (with integrable covariance kernel), and gg is a Lipschitz function. This result was extended in [14] to locally Lipschitz functions gg, such as g(u)=log⁡ug(u)=\log u. This case is particularly important since, by the Cole-Hopf transform, hflat=log⁡uh^{\mathrm{flat}}=\log u is the solution of the Kardar-Parisi-Zhang (KPZ) equation introduced in [26].

The KPZ equation is formally written as:

∂th(t,x)=12∂x2h(t,x)+12(∂xh(t,x))2+ξ(t,x),t>0,x∈R,(2)\partial_t h(t,x)=\frac{1}{2}\partial_x^2h(t,x)+\frac{1}{2}\left(\partial_xh(t,x)\right)^2+\xi(t,x),\quad t>0,x\in\mathbb{R}, \tag*{(2)}

and models the growth of a randomly evolving surface, for instance a surface formed by particles landing and accumulating at random. Its solution h(t,x)h(t,x) represents the height of the surface at time tt and position xx. The diffusion smooths the surface, the nonlinear term makes the growth depend on the local slope, and the noise introduces random irregularities. Much of the interest in the KPZ equation comes from universality: many systems with very different microscopic rules are expected to show the same large-scale fluctuation behaviour. KPZ also reveals connections between random growth, interacting particle systems, directed polymers, and random matrix theory. But defining the solution rigorously is challenging: the solution is too rough for (∂xh)2(\partial_xh)^2 to make sense, so this term requires re-normalization; see [22].

In the present article, we consider the Cole-Hopf transform

h(t,x)=log⁡Z(t,x)(3)h(t,x)=\log Z(t,x) \tag*{(3)}

of the solution of (PAM) with Dirac delta initial condition:

∂tZ(t,x)=12∂x2Z(t,x)+Z(t,x)ξ(t,x),t>0, x∈RZ(0,⋅)=δ0.(4)\partial_tZ(t,x)=\frac{1}{2}\partial_x^2Z(t,x)+Z(t,x)\xi(t,x),\quad t>0,\ x\in\mathbb{R}\qquad Z(0,\cdot)=\delta_0. \tag*{(4)}

We say that h=hnwh=h^{\mathrm{nw}} is the solution of the KPZ (1.2) with narrow wedge initial data, interpreted by approximation (since we cannot take the logarithm of δ0\delta_0): if Zε(0,x)=12εexp⁡(−∣x∣ε)Z_\varepsilon(0,x)=\frac{1}{2\varepsilon}\exp\left(-\frac{|x|}{\varepsilon}\right), then h(0,⋅)=lim⁡ε→0+hε(0,⋅)h(0,\cdot)=\lim_{\varepsilon\to0+}h_\varepsilon(0,\cdot), where

hε(0,x)=log⁡Zε(0,x)=−log⁡(2ε)−∣x∣ε.h_\varepsilon(0,x)=\log Z_\varepsilon(0,x)=-\log(2\varepsilon)-\frac{|x|}{\varepsilon}.

Figure 1 shows some profiles of hεh_\varepsilon which justify the name: as ε\varepsilon decreases, the wedge becomes narrowed and its peak rises.

Narrow wedge profiles

Figure 1. Narrow wedge profiles hε(0,x)=−log⁡(2ε)−∣x∣/εh_\varepsilon(0,x)=-\log(2\varepsilon)-|x|/\varepsilon for ε=1\varepsilon=1, 0.50.5, 0.20.2.

Following Walsh [33], we interpret the solution of (1.4) in the mild sense, i.e. we assume that {Z(t,x);t≥0,x∈R}\{Z(t,x);t\geq0,x\in\mathbb{R}\} is a measurable and (Ft)(\mathcal{F}_t)-adapted random field which satisfies:

Z(t,x)=pt(x)+∫0t∫Rpt−s(x−y)Z(s,y)ξ(ds,dy).(5)Z(t,x)=p_t(x)+\int_0^t\int_{\mathbb{R}}p_{t-s}(x-y)Z(s,y)\xi(\mathrm{d}s,\mathrm{d}y). \tag*{(5)}

where the stochastic integral is interpreted in the Itô sense. The normalization below removes the deterministic spatial profile and allows us to use stationarity. Indeed, the process

Z‾(t,x)=Z(t,x)pt(x),x∈R,\overline{Z}(t,x)=\frac{Z(t,x)}{p_t(x)},\qquad x\in\mathbb{R},

is strictly stationary (see Proposition 1.4 of [2]). Moreover, Z‾\overline{Z} satisfies:

Z‾(t,x)=1+∫0t∫Rps(t−s)t(y−stx)Z‾(s,y)ξ(ds,dy),\overline{Z}(t,x)=1+\int_0^t\int_{\mathbb{R}}p_{\frac{s(t-s)}{t}}\left(y-\frac{s}{t}x\right)\overline{Z}(s,y)\xi(\mathrm{d}s,\mathrm{d}y),

using property (3.3) below of the heat kernel. In [13], the authors proved that its spatial average S‾N(t)=∫0N(Z‾(t,x)−1)dx\overline{S}_N(t)=\int_0^N(\overline{Z}(t,x)-1)\mathrm{d}x satisfies the following FCLT:

1Nlog⁡NS‾N(⋅)→d2B(⋅)in C[0,T],\frac{1}{\sqrt{N\log N}}\overline{S}_N(\cdot)\xrightarrow{d}\sqrt{2}B(\cdot)\quad\text{in }C[0,T],

where BB is a standard Brownian motion. The normalization comes from the fact that Var⁡(S‾N(t))∼2tNlog⁡N\operatorname{Var}(\overline{S}_N(t))\sim2tN\log N, according to Proposition 4.1 of [13]. Alternatively, this follows from stationarity and the fact that

References

  1. [1]Alberts, T., Janjigian, C., Rassoul-Agha, F. and Seppäläinen, T. (2026). The Green’s function of the parabolic Anderson model and the continuum directed polymer. Stoch. Partial Differ. Equ. Anal. Comput. Published online 20 August 2026. doi:10.1007/s40072-026-00440-0.DOI
  2. [2]Amir, G., Corwin, I. and Quastel, J. (2011). Probability distribution of the free energy of the continuum directed random polymer in 1 + 1 dimensions. Comm. Pure Appl. Math. 64, 466–537.
  3. [3]Balan, R. M., Huang, J., Wang, X., Xia, P. and Yuan, W. (2025). Gaussian fluctuations for the wave equation under rough random perturbations. Stochastic Process. Appl. 182, Paper no. 104569.arxiv.org/abs/2307.00103
  4. [4]Balan, R. M., Nualart, D., Quer-Sardanyons, L. and Zheng, G. (2022). The hyperbolic Anderson model: moment estimates of the Malliavin derivatives and applications. Stoch. Partial Differ. Equ. Anal. Comput. 10, 757–827.DOI
  5. [5]Balan, R. M. and Salins, M. (2026). Gaussian fluctuations for the nonlinear stochastic heat equation with drift. Electron. J. Probab. 31, Paper No. 122, 1–31.arxiv.org/abs/2512.12119
  6. [6]Balan, R. M. and Yuan, W. (2024). Hyperbolic Anderson model with time-independent rough noise: Gaussian fluctuations. Electron. J. Probab. 29, Paper No. 177, 1–52.arxiv.org/abs/2305.05043
  7. [7]Balan, R. M. and Zheng, G. (2024). Hyperbolic Anderson model with Lévy white noise: spatial ergodicity and fluctuations. Trans. Amer. Math. Soc. 377, 4171–4221.DOI
  8. [8]Bolaños-Guerrero, R., Nualart, D. and Zheng, G. (2021). Averaging 2d stochastic wave equation. Electron. J. Probab. 26, Paper No. 102, 1–32.DOI
  9. [9]Chen, L. and Dalang, R. C. (2015). Moments and growth indices for the nonlinear stochastic heat equation with rough initial conditions. Ann. Probab. 43, 3006–3051.DOI
  10. [10]Chen, L. and Huang, J. (2019). Comparison principle for stochastic heat equation on ℝ^d. Ann. Probab. 47, 989–1035.arxiv.org/abs/1607.03998
  11. [11]Chen, L. and Jiménez, J. J. (2026). Spatial covariance of KPZ from flat initial profile. Bernoulli, to appear.arxiv.org/abs/2603.14174
  12. [12]Chen, L., Khoshnevisan, D., Nualart, D. and Pu, F. (2022). Central limit theorems for parabolic stochastic partial differential equations. Ann. Inst. Henri Poincaré Probab. Stat. 58, 1052–1077.arxiv.org/abs/1912.01482
  13. [13]Chen, L., Khoshnevisan, D., Nualart, D. and Pu, F. (2022). Spatial ergodicity and central limit theorems for parabolic Anderson model with delta initial condition. J. Funct. Anal. 282, Paper No. 109290.arxiv.org/abs/2005.10417
  14. [14]Chen, L., Khoshnevisan, D., Nualart, D. and Pu, F. (2023). Central limit theorems for spatial averages of the stochastic heat equation via Malliavin-Stein’s method. Stoch. Partial Differ. Equ. Anal. Comput. 11, 122–176.
  15. [15]Chen, L., Ouyang, C., Tindel, S. and Xia, P. (2026). A class of d-dimensional directed polymers in a Gaussian environment. Preprint arXiv:2603.06574.arxiv.org/abs/2603.06574
  16. [16]Conus, D., Joseph, M. and Khoshnevisan, D. (2013). On the chaotic character of the stochastic heat equation, before the onset of intermittency. Ann. Probab. 41, 2225–2260.
  17. [17]Delgado-Vences, F., Nualart, D. and Zheng, G. (2020). A central limit theorem for the stochastic wave equation with fractional noise. Ann. Inst. Henri Poincaré Probab. Stat. 56, 3020–3042.
  18. [18]Ebina, M. (2024). Central limit theorems for nonlinear stochastic wave equations in dimension three. Stoch. Partial Differ. Equ. Anal. Comput. 12, 1141–1200.arxiv.org/abs/2206.12957
  19. [19]Ebina, M. (2025). Central limit theorems for stochastic wave equations in high dimensions. Electron. J. Probab. 30, Paper No. 60, 1–56.arxiv.org/abs/2308.05716
  20. [20]Gu, Y. and Li, J. (2020). Fluctuations of nonlinear stochastic heat equation in dimensions three and higher. SIAM J. Math. Anal. 52, 5422–5440.arxiv.org/abs/1911.01416
  21. [21]Gu, Y. and Pu, F. (2025). Spatial decorrelation of KPZ from narrow wedge. Preprint arXiv:2506.23065.arxiv.org/abs/2506.23065
  22. [22]Hairer, M. (2013). Solving the KPZ equation. Ann. of Math. (2) 178, 559–664.
  23. [23]Huang, J., Nualart, D., and Viitasaari, L. (2020). A central limit theorem for the stochastic heat equation. Stochastic Process. Appl. 130, 7170–7184.arxiv.org/abs/1810.09492
  24. [25]Kallenberg, O. (2002). Foundations of Modern Probability. Second edition. Springer, New York.
  25. [26]Kardar, M., Parisi, G. and Zhang, Y.-C. (1986). Dynamic scaling of growing interfaces. Phys. Rev. Lett. 56, 889–892.pubmed.ncbi.nlm.nih.gov/10033312
  26. [27]Li, J. and Zhang, Y. (2023). The law of the iterated logarithm for spatial averages of the stochastic heat equation. Acta Math. Sci. 43B, 907–918.DOI
  27. [28]Nualart, D. (2006). The Malliavin Calculus and Related Topics. Second edition. Springer-Verlag, Berlin Heidelberg.
  28. [29]Nualart, D., Song, X. and Zheng, G. (2021). Spatial averages for the Parabolic Anderson model driven by rough noise. ALEA Lat. Am. J. Probab. Math. Stat. 18, 907–943.DOI
  29. [30]Nualart, D., Xia, P. and Zheng, G. (2022). Quantitative central limit theorems for the parabolic Anderson model driven by colored noises. Electron. J. Probab. 27, Paper No. 120, 1–43.DOI
  30. [31]Nualart, D. and Zheng, G. (2020). Averaging Gaussian functionals. Electron. J. Probab. 25, Paper No. 48, 1–54.DOI
  31. [32]Nualart, D. and Zheng, G. (2022). Central limit theorems for stochastic wave equations in dimensions one and two. Stoch. Partial Differ. Equ. Anal. Comput. 10, 392–418.DOI
  32. [33]Walsh, J. B. (1986). An introduction to stochastic partial differential equations. École d’Été de Probabilités de Saint-Flour XIV. Lecture Notes in Math. 1180, 265–439. Springer-Verlag, Berlin.

Paper details

vlnc:W_k8wny2qm
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents