alphaFeedback

WorksMathematics

Chaotic CLT using Poisson Fourth Moment Theorem and Graph Theory

AI contributions
Not recorded for this paper
View details

Abstract

In this article, we consider the parabolic Anderson model and the hyperbolic Anderson model in dimension 1, driven by a time-independent Lévy noise, which is either white, or colored with integrable coloration kernel κκ. Using a Poisson version of the Fourth Moment Theorem due to Döbler, Vidotto and Zheng (2018), we prove that the spatial integral of the solution satisfies a central limit theorem. We prove that this result also holds for the hyperbolic Anderson model in dimension d=2d=2, if in addition κκ is bounded. To estimate the fourth cumulant of the projection of the spatial integral on a fixed chaos, we develop a novel technique using graph theory, which allows us to handle some intricate integrals involving the modified contractions.

Introduction

In this article, we consider the parabolic Anderson model in dimension d=1d=1:

{∂u∂t(t,x)=12∂2u∂x2(t,x)+u(t,x)X˙(x),t>0, x∈R,u(0,x)=1,x∈R.(1)\begin{cases} \frac{\partial u}{\partial t}(t,x)=\frac{1}{2}\frac{\partial^2u}{\partial x^2}(t,x)+u(t,x)\dot{X}(x), & t>0,\ x\in\mathbb{R},\\ u(0,x)=1, & x\in\mathbb{R}. \tag*{(1)} \end{cases}

and the hyperbolic Anderson model in dimension d∈{1,2}d\in\{1,2\}:

{∂2u∂t2(t,x)=Δu(t,x)+u(t,x)X˙(x),t>0, x∈Rd,u(0,x)=1,∂u∂t(0,x)=0,x∈Rd.(2)\begin{cases} \frac{\partial^2u}{\partial t^2}(t,x)=\Delta u(t,x)+u(t,x)\dot{X}(x), & t>0,\ x\in\mathbb{R}^d,\\ u(0,x)=1,\quad\frac{\partial u}{\partial t}(0,x)=0, & x\in\mathbb{R}^d. \tag*{(2)} \end{cases}

driven by the (time-independent) Lévy colored noise introduced in the recent preprint [4]. The goal is to prove a central limit theorem (as R→∞R\to\infty) for the spatial average of the solution:

FR(t)=∫BR(u(t,x)−1) dx,(3)F_R(t)=\int_{B_R}\bigl(u(t,x)-1\bigr)\,dx, \tag*{(3)}

where BR={x∈Rd;∣x∣≤R}B_R=\{x\in\mathbb{R}^d;|x|\le R\} and ∣⋅∣|\cdot| is the Euclidean norm. We let σR2(t)=Var⁡(FR(t))\sigma_R^2(t)=\operatorname{Var}(F_R(t)).

This problem has received a lot of attention in the recent years, especially in the case of the Gaussian noise, starting with the seminal article [14]. We refer the reader to [9, 10, 15, 17, 18, 19] for a sample of relevant references. These references rely on a combination of Malliavin calculus with Stein’s method, which yields the precise rate of convergence in the central limit theorem, a result which is now called the Quantitative Central Limit Theorem (QCLT). The major tool is a key estimate for the pp-th moment of the Malliavin derivative of the solution. In the Gaussian case, the time-independent noise can be viewed as a “black” noise, being the opposite of the white noise: if the noise is a fractional Brownian motion in time with index H∈(0,1)H\in(0,1), the case H=1/2H=1/2 corresponds formally to the white noise, while H=1H=1 corresponds to the time-independent noise.

Recently, the QCLT problem has been studied also for the spatial average of the solution of a stochastic partial differential equation (SPDE) driven by a finite-variance Lévy noise, using Malliavin calculus with respect to the underlying Poisson random measure of the noise. References [7, 8] studied the case of the Lévy white noise, while [3] treated the Lévy colored noise. In the present article, we continue these investigations in the case of the time-independent Lévy noise, which may be white or colored in space. This case introduces new challenges which we explain below.

First, defining an Itô integral with respect to the time-independent noise (Gaussian or Lévy) is a delicate issue. In the recent article [1], in dimension d=1d=1, we introduced an Itô integral with respect to a (finite-variance) time-independent Lévy noise LL, using a spatially-indexed filtration (Fx)x∈R(\mathcal{F}_x)_{x\in\mathbb{R}}. But this integral cannot be used for defining the mild solution of the SPDE. Consider for instance the equation Lu(t,x)=u(t,x)L(x)\mathcal{L}u(t,x)=u(t,x)\mathcal{L}(x) with u(0,x)=1u(0,x)=1, where L\mathcal{L} is a differential operator. One cannot define recursively the Picard iteration sequence

un+1(t,x)=1+∫0t∫RGt−s(x−y)un(s,y)L(dy)ds(u0=1),u_{n+1}(t,x)=1+\int_0^t\int_{\mathbb{R}}G_{t-s}(x-y)u_n(s,y)L(dy)ds\quad(u_0=1),

since the integrand is not Fy\mathcal{F}_y-measurable. For instance, in the case of the wave equation, u1(t,x)=1+12L((x−t+s,x+t−s))u_1(t,x)=1+\frac{1}{2}L((x-t+s,x+t-s)) is Fx+t−s\mathcal{F}_{x+t-s}-measurable. Therefore, the only integral that we can consider for defining the solution is the Skorohod integral. This is reviewed in Section 2 below.

Secondly, the traditional methods for proving a QCLT in the Gaussian case do not work in the case of the time-independent Lévy noise. The major problem is that, in the absence of hypercontractivity and Burkholder-Davis-Gundy inequality, there are no techniques for handling pp-th moments of stochastic integrals. The recent article [20] gives a powerful version of the second-order Poincaré inequality for the Poisson case, which has been applied in [7] for proving the QCLT for the wave equation driven by Lévy white noise. But applying this inequality requires an estimate for the pp-th moments of the Malliavin derivative of the solution, which is not available in the time-independent case.

We continue now by describing our framework. We start with the definition of the noise. We let S(Rd)\mathcal{S}(\mathbb{R}^d) be the space of rapidly decreasing functions on Rd\mathbb{R}^d, and S′(Rd)\mathcal{S}'(\mathbb{R}^d) be its dual, the space of tempered distributions. The (time-independent) Lévy colored noise is the collection X={X(φ);φ∈S(Rd)}X=\{X(\varphi);\varphi\in\mathcal{S}(\mathbb{R}^d)\} of random variables, given by:

X(φ)=L(φ∗κ),X(\varphi)=L(\varphi\ast\kappa),

where κ\kappa is a “coloration” kernel, and LL is a (time-independent) Lévy white noise given by:

L(ϕ)=∫Rd×R0ϕ(x)zN^(dx,dz)for all ϕ∈L2(Rd).(4)L(\phi)=\int_{\mathbb{R}^d\times\mathbb{R}_0}\phi(x)z\widehat{N}(dx,dz)\quad\text{for all }\phi\in L^2(\mathbb{R}^d). \tag*{(4)}

Here NN is a Poisson random measure on Rd×R0\mathbb{R}^d\times\mathbb{R}_0 of intensity m=Leb⁡×ν\mathbf{m}=\operatorname{Leb}\times\nu, defined on a complete probability space (Ω,F,P)(\Omega,\mathcal{F},\mathbb{P}), and N^(F)=N(F)−m(F)\widehat{N}(F)=N(F)-\mathbf{m}(F) is the compensated version of NN. The space R0:=R\{0}\mathbb{R}_0:=\mathbb{R}\backslash\{0\} is endowed with the distance d(x,y)=∣x−1−y−1∣d(x,y)=|x^{-1}-y^{-1}|.

In this article, we assume that κ∈L1(Rd)\kappa\in L^1(\mathbb{R}^d), or d=1d=1 and κ=δ0\kappa=\delta_0 (i.e. X=LX=L). Moreover, the measure ν\nu satisfies the condition:

m2:=∫R0∣z∣2ν(dz)<∞.(5)m_2:=\int_{\mathbb{R}_0}|z|^2\nu(dz)<\infty. \tag*{(5)}

If κ∈L1(Rd)\kappa\in L^1(\mathbb{R}^d), X(φ)X(\varphi) is well-defined for any φ∈L2(Rd)\varphi\in L^2(\mathbb{R}^d), since φ∗κ∈L2(Rd)\varphi\ast\kappa\in L^2(\mathbb{R}^d). Letting f=κ∗κ~f=\kappa\ast\widetilde{\kappa}, where κ~(x)=κ(−x)\widetilde{\kappa}(x)=\kappa(-x) for all x∈Rdx\in\mathbb{R}^d, we have: for any φ,ψ∈S(Rd)\varphi,\psi\in\mathcal{S}(\mathbb{R}^d),

⟨φ,ψ⟩0:=∫Rd∫Rdφ(x)ψ(y)f(x−y) dx dy=∫RdFφ(ξ)Fψ(ξ) μ(dξ),(6)\langle\varphi,\psi\rangle_0:=\int_{\mathbb{R}^d}\int_{\mathbb{R}^d}\varphi(x)\psi(y)f(x-y)\,dx\,dy=\int_{\mathbb{R}^d}\mathcal{F}\varphi(\xi)\mathcal{F}\psi(\xi)\,\mu(d\xi), \tag*{(6)}

where μ(dξ)=1(2π)d∣Fκ(ξ)∣2dξ\mu(d\xi)=\frac{1}{(2\pi)^d}|\mathcal{F}\kappa(\xi)|^2d\xi. Here Fφ(x)=∫Rde−iξ⋅xφ(x) dx\mathcal{F}\varphi(x)=\int_{\mathbb{R}^d}e^{-i\xi\cdot x}\varphi(x)\,dx is the Fourier transform of function φ∈L1(Rd)\varphi\in L^1(\mathbb{R}^d), and ξ⋅x\xi\cdot x is the inner product in Rd\mathbb{R}^d. For any φ,ψ∈S(Rd)\varphi,\psi\in\mathcal{S}(\mathbb{R}^d),

E[X(φ)X(ψ)]=m2⟨φ∗κ,ψ∗κ⟩L2(Rd)=m2⟨φ,ψ⟩0.\mathbb{E}[X(\varphi)X(\psi)]=m_2\langle\varphi\ast\kappa,\psi\ast\kappa\rangle_{L^2(\mathbb{R}^d)}=m_2\langle\varphi,\psi\rangle_0.

Intuitively, the solution of (1) should satisfy:

u(t,x)=1+∫0t∫RdGt−s(x−y)u(s,y)X(dy) ds,(7)u(t,x)=1+\int_0^t\int_{\mathbb{R}^d}G_{t-s}(x-y)u(s,y)X(dy)\,ds, \tag*{(7)}

where GtG_t is the fundamental solution of the heat equation:

Gt(x)=1(2πt)d/2exp⁡(−∣x∣22t),G_t(x) = \frac{1}{(2\pi t)^{d/2}} \exp\left(-\frac{|x|^2}{2t}\right),

and ∣⋅∣|\cdot| is the Euclidean norm. Similarly, the solution of (2) should satisfy (7), where GtG_t is the fundamental solution of the wave equation:

Gt(x)=121{∣x∣<t}if d=1,G_t(x) = \frac{1}{2}\mathbf{1}_{\{|x|<t\}} \quad\text{if } d=1,
Gt(x)=12π1t2−∣x∣21{∣x∣<t}if d=2.G_t(x) = \frac{1}{2\pi}\frac{1}{\sqrt{t^2-|x|^2}}\mathbf{1}_{\{|x|<t\}} \quad\text{if } d=2.

The rigorous definition of the solution is given Section 2 below. From [4], it is known that the solutions of (1) and (2) exist, provided that Dalang’s condition holds, i.e.

∫Rd11+∣ξ∣2 μ(dξ)<∞.(8)\int_{\mathbb{R}^d} \frac{1}{1+|\xi|^2}\,\mu(\mathrm{d}\xi)<\infty. \tag*{(8)}

Note that condition (8) always holds in dimension d=1d=1.

We are now ready to state our main results.

References

  1. [1]Balan, R.M. and Garza, J. (2026). An Itô integral for a two-sided Lévy process. J. Stoch. Anal. 7, no.2.DOI
  2. [2]Balan, R.M. and Jiménez, J.J. (2026). Moment estimates for solutions of SPDEs with Lévy colored noise. Stoch. PDEs: Anal. & Comp. 14, 1910-1969.arxiv.org/abs/2504.21648
  3. [3]Balan, R.M. and Stephenson, W. (2026). Gaussian fluctuations for hyperbolic Anderson model with Lévy colored noise. Preprint arxiv:2602.23137.arxiv.org/abs/2602.23137
  4. [4]Balan, R.M. and Wang, J. (2026). SPDEs with time-independent Lévy colored noise. Preprint arXiv:2604.24914.arxiv.org/abs/2604.24914
  5. [5]Balan, R.M. and Yuan, W. (2022). Spatial integral of the solution to hyperbolic Anderson model with time-independent noise. Stoch. Proc. Their Appl. 152, 177-207.arxiv.org/abs/2201.02319
  6. [6]Balan, R.M. and Yuan, W. (2023). Central limit theorems for heat equation with time-independent noise: The regular and rough cases. Inf. Dim. Anal. Quantum Probab. Rel. Topics 26, no. 2, 2250029.arxiv.org/abs/2205.13105
  7. [7]Balan, R.M. and Zheng, G. (2024). Hyperbolic Anderson model with Lévy white noise: spatial ergodicity and fluctuation. Trans. AMS 377, 4171-4221.DOI
  8. [8]Balan, R.M. and Zheng, G. (2026). Central limit theorem for stochastic nonlinear wave equation with pure-jump Lévy white noise. Electr. J. Probab. 31, article no. 70, 1-46.DOI
  9. [9]Bolaños-Guerrero, R., Nualart, D. and Zheng, G. (2021). Averaging 2d stochastic wave equation. Electr. J. Probab. 26, paper no. 102, 1-32.DOI
  10. [10]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.
  11. [11]Döbler, C. and Peccati, G. (2018). The fourth moment theorem on the Poisson space. Ann. Probab. 46, 1878-1916.arxiv.org/abs/1701.03120
  12. [12]Döbler, C. and Peccati, G. (2018). Fourth moment theorems on the Poisson space: analytic statements via product formulae. Electr. Comm. Probab. 23, no. 91, 1-12.arxiv.org/abs/1808.01836
  13. [13]Döbler, C., Vidotto, A. and Zheng, G. (2018). Fourth moment theorems on the Poisson space in any dimension. Electr. J. Probab. 23, no. 36, 1-27.DOI
  14. [14]Huang J., Nualart D. and Viitasaari L. (2020). A central limit theorem for the stochastic heat equation. Stoch. Proc. Their Appl. 130, 7170-7184.arxiv.org/abs/1810.09492
  15. [15]Huang, J., Nualart, D., Viitasaari L. and Zheng, G. (2020). Gaussian fluctuations for the stochastic heat equation with colored noise. Stoch. PDE: Anal. Comp. 8, 402-421.
  16. [16]Nourdin, I. and Peccati, G. (2012). Normal Approximations with Malliavin Calculus. Cambridge University Press, Cambridge.DOI
  17. [17]Nualart, D., Song, X. and Zheng, G. (2021). Spatial averages for the Parabolic Anderson model driven by rough noise. ALEA, Latin Amer. J. Probab. Math. Stat. 18, 907-943.DOI
  18. [18]Nualart, D., Xia, P. and Zheng, G. (2022). Quantitative central limit theorems for the parabolic Anderson model driven by colored noises. Electr. J. Probab. 27, paper No. 120, 1-43.DOI
  19. [19]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
  20. [20]Trauthwein, T. (2025). Quantitative CLTs on the Poisson space via Skorohod estimates and p-Poincaré inequalities. Ann. Appl. Probab. 35, 1716-1754.arxiv.org/abs/2212.03782

Paper details

vlnc:W_e9mmkzt6
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