A law of large numbers concerning the number of critical points of isotropic Gaussian functions

Abstract.

Let $\Phi$ be a smooth isotropic random Gaussian function $\Phi$ on ${\mathbb{R}}^{m}$. Denote by $Z_{N}(\Phi)$ the number of critical points of $\Phi$ inside the cube $[0,N]^{m}$. We prove that $N^{-m}Z_{N}(\Phi)$ converges $\mathrm{a.s.}$ and $L^{2}$ to a universal explicit constant $C_{m}(\Phi)$.

1. Introduction

Denote by $\operatorname{Meas}({\mathbb{R}}^{m})$ the space of finite Borel measures on ${\mathbb{R}}^{m}$. Suppose that $\mathfrak{a}:{\mathbb{R}}\to{\mathbb{R}}$ is an even Schwartz function such that $\mathfrak{a}(0)=1$. Consider the finite measure $\mu\in\operatorname{Meas}({\mathbb{R}}^{m})$

$$\mu\big{[}\,d\xi\,\big{]}=\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}=\frac{1}{(2\pi)^{m}}w_{\mathfrak{a},m}\big{(}\,\xi\,\big{)}{\boldsymbol{\lambda}}\big{[}\,d\xi\,\big{]},\;\;w_{\mathfrak{a},m}\big{(}\,\xi\,\big{)}=\mathfrak{a}\big{(}\,|\xi|\,\big{)}^{2}.$$

Its characteristic function is the nonnegative definite function

$${\boldsymbol{K}}\big{(}\,{\boldsymbol{x}}\,\big{)}={\boldsymbol{K}}_{\mathfrak{a}}\big{(}\,{\boldsymbol{x}}\,\big{)}=\int_{{\mathbb{R}}^{m}}e^{\boldsymbol{i}\langle\xi,{\boldsymbol{x}}\rangle}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}=\frac{1}{(2\pi)^{m}}\int_{{\mathbb{R}}^{m}}e^{\boldsymbol{i}\langle\xi,{\boldsymbol{x}}\rangle}\mathfrak{a}\big{(}\,|\xi|\,\big{)}^{2}{\boldsymbol{\lambda}}\big{[}\,d\xi\,\big{]}.$$ (1.1)

Clearly ${\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}})$ is an $O(n)$-invariant, real valued Schwartz function. Then ${\boldsymbol{K}}_{\mathfrak{a}}\big{(}\,{\boldsymbol{x}}-{\boldsymbol{y}}\,\big{)}$ is the covariance kernel of a real valued, smooth isotropic Gaussian function $\Phi=\Phi_{\mathfrak{a}}$ on ${\mathbb{R}}^{m}$ with spectral measure $\mu_{\mathfrak{a}}$.

A box is a subset of ${\mathbb{R}}^{m}$ of the form $[a_{1},b_{1}]\times\cdots\times[a_{m},b_{m}]$, $a_{i}<b_{i}$, $\forall i$. For any box $B$ we denote by $Z_{\mathfrak{a}}\big{[}\,B\,\big{]}$ the number of critical points of $\Phi_{\mathfrak{a}}$ in $B$. Then $\mathrm{a.s.}$ $\Phi_{\mathfrak{a}}$ has no critical points on $\partial B$ and all the critical points of $\Phi_{\mathfrak{a}}$ in $B$ are nondegenerate; see [1, 3].

There exists a universal explicit constant $C_{m}(\mathfrak{a})>0$ such that for any box $B$ we have

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B)\,\big{]}=C_{m}(\mathfrak{a})\operatorname{vol}\big{[}\,B\,\big{]}.$$

Fix numbers $a_{1},\dotsc,a_{m}>0$ and denote by $B$ the box $[0,a_{1}]\times\cdots[0,a_{m}]$. For ${\vec{\ell}}\in{{\mathbb{N}}}$ we denote by $B_{\vec{\ell}}$ the box

$$B_{\vec{\ell}}=\prod_{j=1}^{m}\big{[}\,(\ell_{j}-1)a_{j},\ell_{j}a_{j}\,\big{]}.$$

For $N\in{{\mathbb{N}}}$ we denote by $B_{N}$ the box

$$B_{N}=N\cdot B=\prod_{j=1}^{m}\big{[}\,0,Na_{j}\,\big{]}.$$

The main goal of this paper is to prove the following strong law of large numbers.

Theorem 1.1.

As $N\to\infty$, the random variable

$$\frac{1}{N^{m}}Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}$$

converges $\mathrm{a.s.}$ and $L^{2}$ to the (deterministic) constant

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B)\,\big{]}=C_{m}(\mathfrak{a})\operatorname{vol}\big{[}\,B\,\big{]}.$$

$\hbox to0.0pt{$\sqcap$\hss}\sqcup$

Let us offer a different perspective on this result. For $\hbar>0$ we set

$$\mathfrak{a}_{\hbar}(t):=\mathfrak{a}(\hbar t\,\big{)},\;\;\forall t\in{\mathbb{R}}.$$

Consider the finite Borel measure $\mu_{\mathfrak{a}}^{\hbar}\in\operatorname{Meas}({\mathbb{R}}^{m})$

$$\mu_{\mathfrak{a}}^{\hbar}\big{[}\,d\xi\,\big{]}=\frac{1}{(2\pi)^{m}}w_{{\boldsymbol{a}},m}\big{(}\,\hbar\xi\,\big{)}{\boldsymbol{\lambda}}\big{[}\,d\xi\,\big{]}=\frac{1}{(2\pi)^{m}}w_{{\boldsymbol{a}}_{\hbar},m}\big{(}\,\xi\,\big{)}{\boldsymbol{\lambda}}\big{[}\,d\xi\,\big{]}.$$

Its characteristic function is the nonnegative definite function

$${\boldsymbol{K}}^{\hbar}\big{(}\,{\boldsymbol{x}}\,\big{)}={\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}\big{(}\,{\boldsymbol{x}}\,\big{)}=\frac{1}{(2\pi)^{m}}\int_{{\mathbb{R}}^{m}}e^{\boldsymbol{i}\langle\xi,{\boldsymbol{x}}\rangle}\mathfrak{a}\big{(}\,|\hbar\xi|\,\big{)}^{2}d\xi.$$ (1.2)

We set $\omega:=\hbar\xi$ in (1.2) and we deduce

$${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}\big{(}\,{\boldsymbol{x}}\,\big{)}=\frac{1}{(2\pi\hbar)^{m}}\int_{{\mathbb{R}}^{m}}e^{\frac{\boldsymbol{i}}{\hbar}\langle\omega,{\boldsymbol{x}}\rangle}\mathfrak{a}\big{(}\,|\omega|\,\big{)}^{2}d\omega,$$

so that

$${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}({\boldsymbol{x}})=\hbar^{-m}{\boldsymbol{K}}_{\mathfrak{a}}\big{(}\,\hbar^{-1}{\boldsymbol{x}}\,\big{)}.$$

We deduce that ${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}\big{(}\,{\boldsymbol{x}}-{\boldsymbol{y}}\,\big{)}$ is the covariance kernel of the Gaussian function

$$\Phi_{\mathfrak{a}}^{\hbar}({\boldsymbol{x}}):=\hbar^{-m/2}\Phi_{\mathfrak{a}}\big{(}\,\hbar^{-1}{\boldsymbol{x}}\,\big{)}.$$

We want to investigate the behavior of ${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}({\boldsymbol{x}})$ as $\hbar\searrow 0$. Using the Fourier inversion formula we deduce

$$\int_{{\mathbb{R}}^{m}}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}})d{\boldsymbol{x}}=\mathfrak{a}(0)^{2}=1.$$

This implies that the measures ${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}\big{(}\,{\boldsymbol{x}}\,\big{)}dx$ converge weakly to the Dirac measure $\delta_{0}$. To use a terminology favored by physicists, we have

$${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}\big{(}\,{\boldsymbol{x}}\,\big{)}\to\delta({\boldsymbol{x}}),$$

where $\delta({\boldsymbol{x}})$ is Dirac’s mysterious Delta function. In particular,

$${\boldsymbol{K}}_{\mathfrak{a}}^{\hbar}\big{(}\,{\boldsymbol{x}}-{\boldsymbol{y}}\,\big{)}\to\delta({\boldsymbol{x}}-{\boldsymbol{y}}).$$

In other words, as $\hbar\searrow 0$, the Gaussian random function $\Phi_{\mathfrak{a}}^{\hbar}$ converges in some sense to a Gaussian random “function” $\Phi^{0}$ whose covariance kernel is ${\boldsymbol{K}}^{0}({\boldsymbol{x}}-{\boldsymbol{y}})=\delta({\boldsymbol{x}}-{\boldsymbol{y}})$. This is a Gaussian white noise.

Note that

$$Z_{\mathfrak{a}_{\hbar}}\big{[}\,B\,\big{]}=Z_{\mathfrak{a}}\big{[}\,\hbar^{-1}B\,\big{]}.$$

Consider the random measure

$$\nu_{\hbar}=\frac{\hbar^{m}}{C_{m}(\mathfrak{a})}\sum_{\nabla\Phi_{\mathfrak{a}_{\hbar}}({\boldsymbol{x}})=0}\delta_{\boldsymbol{x}}.$$

Then

$$\frac{\hbar^{m}}{C_{m}(\mathfrak{a})}Z_{\mathfrak{a}_{\hbar}}\big{[}\,B\,\big{]}=\nu_{\hbar}\big{[}\,B\,\big{]}.$$

Now let $\hbar=N^{-1}$. In this case $Z_{\mathfrak{a}_{\hbar}}\big{[}\,B\,\big{]}=Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}$. Theorem 1.1 shows that in the white noise limit ($N\to\infty$) the random measure $\nu_{N^{-1}}$ converges $\mathrm{a.s.}$ to the deterministic measure ${\boldsymbol{\lambda}}$, the Lebesgue measure on ${\mathbb{R}}^{m}$. In other words, the critical points of $\Phi_{\mathfrak{a}_{\hbar}}$ will $\mathrm{a.s.}$ equidistribute as $\hbar\searrow 0$.

In the earlier work [7] we proved that the random variable $Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}$ is highly concentrated around its mean. More precisely we showed that $N\to\infty$ we have

$$\operatorname{Var}\big{[}\,Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}\,\big{]}\sim c\cdot N^{m},$$

where $c$ is a positive constant and the renormalized random variables

$$N^{-m/2}\big{(}\,Z_{\mathfrak{a}}\big{[}\,B\,\big{]}-{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}\big{[}\,B\,\big{]}\,\big{]}\,\big{)}$$

converge in distribution to a normal random variable with nonzero standard deviation.

The strategy we use to prove Theorem 1.1 is very simple. Set ${\mathbb{I}}_{N}:=[1,N]\cap{{\mathbb{N}}}$.

$$\frac{1}{N^{m}}Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}=\frac{1}{N^{m}}\sum_{{\vec{\ell}}\in{\mathbb{I}}_{N}^{m}}Z_{\mathfrak{a}}\big{[}\,B_{\vec{\ell}}\,\big{]}.$$

Since the Gaussian function $\Phi_{\mathfrak{a}}$ is stationary, the family of random variables $\big{(}\,Z_{\mathfrak{a}}\big{[}\,B_{\vec{\ell}}\,\big{]}\,\big{)}_{{\vec{\ell}}\in{\mathbb{Z}}^{m}}$ is stationary. Using the multiparametric ergodic theorem [6, Chap.6, Thm.2.5] we conclude that $\frac{1}{N^{m}}Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}$ converge $\mathrm{a.s.}$ and $L^{1}$ to some random variable $Z_{\infty}$. To show that $Z_{\infty}$ is $\mathrm{a.s.}$ constant we prove that

$$\operatorname{Var}\big{[}\,N^{-m}Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}\,\big{]}\to 0\;\;\mbox{as $N\to\infty$}.$$

We achieve this by observing that

$$\operatorname{Var}\big{[}\,N^{-m}Z_{\mathfrak{a}}\big{[}\,N\cdot B\,\big{]}=N^{-2m}\sum_{{\vec{k}},{\vec{\ell}}\in{\mathbb{I}}_{N}^{m}}\underbrace{\operatorname{Cov}\big{[}\,Z_{\mathfrak{a}}\big{[}\,B_{\vec{k}}\,\big{]},Z_{\mathfrak{a}}\big{[}\,B_{\vec{\ell}}\,\big{]}}_{=:C({\vec{k}},{\vec{\ell}})}.$$

The above sum consists of $N^{2m}$ terms, but the terms that are far from the diagonal ${\vec{k}}={\vec{\ell}}$ contribute very little. More precisely, using the Kac-Rice formula we will prove that for any $p>0$ there exists $C=C(p)>0$ such that

$$\big{|}\,C({\vec{k}},{\vec{\ell}})\,\big{|}\leq\frac{C(p)}{\big{(}\,1+|{\vec{k}}-{\vec{\ell}}|_{1}\,\big{)}^{p}},\;\;\forall({\vec{k}},{\vec{\ell}})\in{{\mathbb{N}}}^{m}\times{{\mathbb{N}}}^{m},$$

where

$$\big{|}\,{\boldsymbol{x}}\,\big{|}_{1}=\sum_{i=1}^{m}\big{|}\,x_{i}\,\big{|},\;\;\forall{\boldsymbol{x}}=(x_{1},\dotsc,x_{m})\in{\mathbb{R}}^{m}.$$

According to Lemma 3.1, this implies that

$$\sum_{{\vec{k}},{\vec{\ell}}\in\mathcal{R}_{N}}C({\vec{k}},{\vec{\ell}})=O\big{(}\,N^{2m-1}\,\big{)}\;\;\mbox{as $N\to\infty$}.$$

Theorem 1.1 is valid for a wider family of homogeneous Gaussian functions on ${\mathbb{R}}^{m}$ with fewer constraints on their spectral functions. We chose to work in this more restrictive context for two reasons. First, relaxing the assumptions on the Gaussian function would have added a few more technical complications that would have blurred the main arguments.

The second reason is that the class of Gaussian functions discussed in this paper is exactly the class of function that arises when studying critical points of certain random Fourier series of eigenvalues on Riemann manifolds.

More precisely given a Riemann manifold $(M,g)$ with Laplacian $\Delta$ we denote by $\mathcal{K}^{g}_{\mathfrak{a}}$ the Schwartz kernel of the smoothing operator $\mathfrak{a}(\Delta_{g})^{2}$. This the the covariance kernel of a smooth random function on $M$ that can be represented as a random Fourier series. Upon rescaling the metric $g\to g_{\hbar}=\hbar^{-2}g$ we have $\Delta_{g_{\hbar}}=\hbar^{2}\Delta_{g}$ and, as $\hbar\searrow 0$, the kernel $\mathcal{K}_{\mathfrak{a}}^{g_{\hbar}}$ converges in a precise sense to the kernel ${\boldsymbol{K}}_{\mathfrak{a}}$ in (1.1).

The paper is organized as follows. In Section 2 we describe some basic properties of $\Phi_{\mathfrak{a}}$. Section 3 describes an abstract weak law of large numbers for families of mean zero random variables $(X_{\vec{\ell}})_{{\vec{\ell}}\in{\mathbb{Z}}^{m}}$ such that $X_{\vec{\ell}}$, $X_{\vec{k}}$ are weakly correlated if ${\vec{k}}$ and ${\vec{\ell}}$ are far apart. The main theorem is an immediate consequence of this result.

2. Some properties of the function $\Phi_{\mathfrak{a}}$

We collect in this section several useful properties of the random function $\Phi_{\mathfrak{a}}$ used throughout the paper. For $k\in{{\mathbb{N}}}$ we denote by $D^{k}\Phi_{\mathfrak{a}}$ the $k$-th order differential of $\Phi_{\mathfrak{a}}$

Proposition 2.1.

Let $N\in{{\mathbb{N}}}$. Then the following hold.

1. The function $\Phi_{\mathfrak{a}}$ is $N$-ample, i.e., for any are distinct points ${\boldsymbol{x}}_{1},\dotsc,{\boldsymbol{x}}_{N}\in{\mathbb{R}}^{m}$, the Gausian vector

   $$\big{(}\,\Phi_{\mathfrak{a}}({\boldsymbol{x}}_{1}),\dotsc,\Phi_{\mathfrak{a}}({\boldsymbol{x}}_{N})\,\big{)}.$$

   is nondegenerate.

2. The function $\Phi_{\mathfrak{a}}$ is $J_{N}$-ample, i.e., the Gaussian vector

   $$\Phi_{\mathfrak{a}}({\boldsymbol{x}})\oplus D\Phi_{\mathfrak{a}}({\boldsymbol{x}})\oplus\cdots\oplus D^{N}\Phi_{\mathfrak{a}}({\boldsymbol{x}})$$

   is nondegenerate for any ${\boldsymbol{x}}\in{\mathbb{R}}^{m}$.

Proof.

(i) Since $w_{\mathfrak{a},m}\big{(}\,\big{|}\,\xi\,\big{|}\,\big{)}=\mathfrak{a}\big{(}\,|\xi|\,\big{)}^{2}$ is positive on an open neighborhood of $0\in{\mathbb{R}}^{m}$ we deduce from [8, Thm. 6.8] that if ${\boldsymbol{x}}_{1},\dotsc,{\boldsymbol{x}}_{N}\in{\mathbb{R}}^{m}$ are distinct points, then the symmetric $N\times N$ matrix

$$\big{(}\,K({\boldsymbol{x}}_{i}-{\boldsymbol{x}}_{j}\,\big{)}_{1\leq i,j\leq N}$$

is positive definite. This matrix is the variance matrix of the Gaussian vector

$$\big{(}\,\Phi_{\mathfrak{a}}({\boldsymbol{x}}_{1}),\dotsc,\Phi_{\mathfrak{a}}({\boldsymbol{x}}_{N})\,\big{)}.$$

(ii) Observe that for any multi-indices $\alpha\in\big{(}\,{\mathbb{Z}}_{\geq 0}\,\big{)}^{m}$, we have

$${\mathbb{E}}\big{(}\,\partial^{\alpha}\Phi_{\mathfrak{a}}({\boldsymbol{x}})\partial^{\beta}\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{)}=\partial^{\alpha}_{x}\partial^{\beta}_{y}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}}-{\boldsymbol{y}}\,\big{)}\big{|}_{{\boldsymbol{x}}={\boldsymbol{y}}}$$

$$=\int_{{\mathbb{R}}^{m}}\xi^{\alpha}\xi^{\beta}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]},\;\;\xi^{\alpha}:=\xi_{1}^{\alpha_{1}}\cdots\xi_{m}^{\alpha_{m}}$$

This shows that for any $N\in{{\mathbb{N}}}$ and any ${\boldsymbol{x}}\in{\mathbb{R}}^{n}$ the variance the Gaussian vector $\big{(}\,\partial^{\alpha}\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{)}_{|\alpha|\leq N}$ is the Gramian matrix of the functions $\big{(}\,\xi^{\alpha}\,\big{)}_{|\alpha|\leq N}$ with respect to the inner product in $L^{2}\big{(}\,{\mathbb{R}}^{m},\mu_{\mathfrak{a}}\,\big{)}$. Since $\mathfrak{a}(0)=1$ we deduce that the functions $\xi^{\alpha}$ are linearly independent in $L^{2}\big{(}\,{\mathbb{R}}^{m},\mu_{\mathfrak{a}}\,\big{)}$ so the determinant of their Gramian matrix is nonzero. Hence the Gaussian vector

$$\Phi_{\mathfrak{a}}({\boldsymbol{x}})\oplus D\Phi_{\mathfrak{a}}({\boldsymbol{x}})\oplus\cdots\oplus D^{N}\Phi_{\mathfrak{a}}({\boldsymbol{x}})$$

is nondegenerate, for any $k\in{{\mathbb{N}}}$ and any ${\boldsymbol{x}}\in{\mathbb{R}}^{m}$. $\hbox to0.0pt{$\sqcap$\hss}\sqcup$

We deduce from the above proposition and the results of Ancona-Letendre [2] or Gass-Stecconi [4] that for any box $B$ the random variable $Z_{\mathfrak{a}}\big{[}\,B\,\big{]}$ has finite moments of any order. To compute the expectation of $Z_{\mathfrak{a}}(B)$ we rely on the Kac-Rice formula [1, 3]

Theorem 2.2.

Let ${\boldsymbol{V}}$ be a finite dimensional Euclidean space, $\mathcal{V}\subset{\boldsymbol{V}}$ an open set and $F:\mathcal{V}\to{\mathbb{R}}$ a Gaussian random function that is $\mathrm{a.s.}$ $C^{2}$ and such that the Gaussian vector $\nabla F({\boldsymbol{v}})$ is nondegenerate for any ${\boldsymbol{v}}\in\mathcal{V}$. We denote by $p_{\nabla F(v)}$ is probability density. For a subset $S\subset\mathcal{V}$ we denote by $Z(S,F)$ the number of critical points of $F$ in $S$. Then for any box $B\subset\mathcal{V}$, the function $F$ $\mathrm{a.s.}$ has no critical points on $\partial B$, all the critical points are nondegenerate and

$${\mathbb{E}}\big{[}\,Z(B,F)\,\big{]}=\int_{B}{\mathbb{E}}\big{[}\,|\det\operatorname{Hess}F({\boldsymbol{v}})|\,\,\big{|}\nabla F({\boldsymbol{v}})=0\,\big{]}p_{\nabla F(v)}(0)d{\boldsymbol{v}}.$$ (2.1)

$\hbox to0.0pt{$\sqcap$\hss}\sqcup$

We want to apply the above result to the function $\Phi_{\mathfrak{a}}$. Prooposition 2.1 shows that the Gaussian vector $\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})$ is nondegenerate for any ${\boldsymbol{x}}\in{\mathbb{R}}^{m}$.

For any multi-indices $\alpha,\beta\in\big{(}\,{\mathbb{Z}}_{\geq 0}\,\big{)}^{m}$ we have

$${\mathbb{E}}\big{[}\,\partial^{\alpha}\Phi_{\mathfrak{a}}({\boldsymbol{x}})\partial^{\beta}\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{]}_{{\boldsymbol{x}}={\boldsymbol{y}}}=\partial^{\alpha}_{\boldsymbol{x}}\partial^{\beta}_{\boldsymbol{y}}\mathcal{K}^{\mathfrak{a}}({\boldsymbol{x}},{\boldsymbol{y}})\big{|}\,_{{\boldsymbol{x}}={\boldsymbol{y}}}=\frac{(-1)^{|\beta||}\boldsymbol{i}^{|\alpha|+|\beta|}}{(2\pi)^{m}}\int_{{\mathbb{R}}^{m}}\xi^{\alpha+\beta}\mathfrak{a}(|\xi|)^{2}d\xi.$$ (2.2)

For any multi-index $\alpha\in\big{(}\,{\mathbb{Z}}_{\geq 0}\,\big{)}^{m}$ we set

$$M^{\mathfrak{a}}_{\gamma}:=\int_{{\mathbb{R}}^{m}}\xi^{\alpha}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}=\frac{1}{(2\pi)^{m}}\int_{{\mathbb{R}}^{m}}\xi^{\alpha}\mathfrak{a}(|\xi|)^{2}d\xi$$

We say that the multi-index $\alpha=(\alpha_{1},\dotsc,\alpha_{m})$ is even if $\alpha_{j}$ is even for any $j=1,\dotsc,m$. The multi-index $\alpha$ is called odd if it is not even. The radial symmetry of $\mathfrak{a}\big{(}\,|\xi|\,\big{)}$ implies that

$$M^{\mathfrak{a}}_{\alpha}=0\;\;\mbox{if $\alpha$ is odd}.$$ (2.3)

Using spherical coordinates on ${\mathbb{R}}^{m}$ we deduce that for any $\alpha$ we have

$$M^{\mathfrak{a}}_{\alpha}=\frac{1}{(2\pi)^{m}}\left(\int_{0}^{\infty}r^{m-1+|\alpha|}\mathfrak{a}(r)^{2}dr\right)\times\underbrace{\int_{S^{m-1}}\xi^{\alpha}\operatorname{vol}_{S^{m-1}}\big{[}\,d\xi\,\big{]}}_{=:\boldsymbol{m}_{\alpha}}.$$ (2.4)

Note that $\boldsymbol{m}_{\alpha}$ is independent of $\mathfrak{a}$. If we let $\mathfrak{a}_{0}:=(2\pi)^{m/2}e^{-\frac{t^{2}}{4}}$, then

$$M^{\mathfrak{a}_{0}}_{\alpha}=\int_{{\mathbb{R}}^{m}}\xi^{\gamma}e^{-|\xi|^{2}/2}d\xi=(2\pi)^{m/2}\prod_{j=1}^{m}\int_{\mathbb{R}}\xi^{\alpha_{j}}{\boldsymbol{\gamma}}_{1}\big{[}\,d\xi\,\big{]}$$

where ${\boldsymbol{\gamma}}_{1}$ denotes the Gaussian measure on ${\mathbb{R}}$ with mean zero and variance $1$. If $\alpha$ is even, $\alpha=2\kappa$, then

$$M^{\mathfrak{a}_{0}}_{2\kappa}=(2\pi)^{m/2}\prod_{j=1}^{m}(2\kappa_{j}-1)!!.$$

On the other hand, using (2.4) we deduce

$$M^{\mathfrak{a}_{0}}_{2\kappa}=\boldsymbol{m}_{2\kappa}\int_{0}^{\infty}r^{m+2|\kappa|-1}e^{-r^{2}/2}dr=\sqrt{\frac{\pi}{2}}\boldsymbol{m}_{2\kappa}\int_{\mathbb{R}}|x|^{m+2|\kappa|-1}{\boldsymbol{\gamma}}_{1}\big{[}\,dx\,\big{]}$$

$$=2^{|\kappa|+m/2-1}\boldsymbol{m}_{2\kappa}\Gamma\big{(}\,|\kappa|+m/2\,\big{)}.$$

Hence

$$\boldsymbol{m}_{2\kappa}=2^{|\kappa|+m/2-1}\frac{\prod_{j=1}^{m}(2\kappa_{j}-1)!!}{\Gamma\big{(}\,|\kappa|+m/2\,\big{)}}=\frac{2\prod_{j=1}\Gamma(\kappa_{j}+1/2)}{\Gamma\big{(}\,|\kappa|+m/2\,\big{)}}.$$ (2.5)

For every $k\in{\mathbb{Z}}_{\geq 0}$ we set

$$I_{k}(\mathfrak{a}):=\int_{0}^{\infty}r^{k}\mathfrak{a}(r)^{2}dr.$$

We deduce

$$(2\pi)^{m}M^{\mathfrak{a}}_{2\kappa}=2I_{m-1+2|\kappa|}(\mathfrak{a})\frac{2\prod_{j=1}\Gamma(\kappa_{j}+1/2)}{\Gamma\big{(}\,|\kappa|+m/2\,\big{)}}.$$ (2.6)

We set

$$s_{m}:=\int_{{\mathbb{R}}^{m}}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]},\;\;d_{m}=\int_{{\mathbb{R}}^{m}}\xi_{1}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]},\;\;h_{m}:=\int_{{\mathbb{R}}^{m}}\xi_{1}^{2}\xi_{2}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}.$$ (2.7)

Then

$$\int_{{\mathbb{R}}^{m}}\mathfrak{a}(|\xi|)^{2}d\xi=\frac{2\pi^{m/2}}{\Gamma(m/2)}I_{m-1}(\mathfrak{a})=(2\pi)^{m}s_{m},$$ (2.8)

$$\int_{{\mathbb{R}}^{m}}\xi_{j}^{2}\mathfrak{a}(|\xi|)^{2}d\xi=\frac{2\pi^{m/2}}{\Gamma(m/2+1)}I_{m+1}(\mathfrak{a})=(2\pi)^{m}d_{m},\;\;\forall j,$$ (2.9)

$$\int_{{\mathbb{R}}^{m}}\xi_{j}^{2}\xi_{k}^{2}\mathfrak{a}(|\xi|)^{2}d\xi=\frac{(2\pi)^{m/2}}{\Gamma(m/2+2)}I_{m+3}(\mathfrak{a})=(2\pi)^{m}h_{m},\;\;\forall j\neq k,$$ (2.10)

$$\int_{{\mathbb{R}}^{m}}\xi_{j}^{4}\mathfrak{a}(|\xi|)^{2}d\xi=\frac{6\pi^{m/2}}{\Gamma(m/2+1)}I_{m+3}(\mathfrak{a})=3(2\pi)^{m}h_{m},\;\;\forall j.$$ (2.11)

Using (2.2) and (2.3) we deduce that for any ${\boldsymbol{x}}\in{\mathbb{R}}^{m}$ the Gaussian vectors $\nabla\Phi({\boldsymbol{x}})$ and $\operatorname{Hess}_{\Phi}({\boldsymbol{x}})$ are independent.

Using (2.2)and (2.9) we deduce that the variance matrix of $\nabla\Phi({\boldsymbol{x}})$ is

$$\operatorname{Var}\big{[}\,\nabla\Phi({\boldsymbol{x}})\,\big{]}=d_{m}{\mathbbm{1}}_{m},\;\;\forall{\boldsymbol{x}}\in{\mathbb{R}}^{m},$$ (2.12)

where ${\mathbbm{1}}_{m}$ denotes the identity $m\times m$ matrix.

The space $\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})$ of real symmetric $m\times m$ matrices is equipped with an inner product $(A,B)=\operatorname{{\rm tr}}(A)$. Moreover the linear functions $\ell_{ij},\omega_{ij}:\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})\to{\mathbb{R}}$, $1\leq i\leq j\leq m$,

$$\ell_{ij}(A)=a_{ij},\;\;\omega_{ij}(A)=\begin{cases}a_{ii},&i=j,\\ \sqrt{2}a_{ij},&i<j\end{cases}$$ (2.13)

define coordinates on $\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})$ that are orthonormal with respect to the above inner product. We set

$$L_{ij}({\boldsymbol{x}}):=\ell_{ij}\big{(}\,\operatorname{Hess}_{\Phi}({\boldsymbol{x}})\,\big{)},\;\;\Omega_{ij}({\boldsymbol{x}}):=\omega_{ij}\big{(}\,\operatorname{Hess}_{\Phi}({\boldsymbol{x}})\,\big{)}.$$ (2.14)

Then

$${\mathbb{E}}\big{[}\,\partial^{2}_{x_{i}x_{j}}\Phi_{\mathfrak{a}}(x)\partial^{2}_{x_{k}x_{\ell}}\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}=\frac{1}{(2\pi)^{m}}\int_{{\mathbb{R}}^{m}}\xi_{i}\xi_{j}\xi_{k}\xi_{\ell}a\big{(}\,|\xi|^{2}\,\big{)}d\xi,\;\;i\leq j,\;\;k\leq\ell.$$

Note that if $i<j$, then the above integral is nonzero iff $(i,j)=(k,\ell)$ in which case

$${\mathbb{E}}\big{[}\,L_{ij}({\boldsymbol{x}})L_{ij}({\boldsymbol{x}})\,\big{]}={\mathbb{E}}\big{[}\,\big{(}\,\partial^{2}_{x_{i}x_{j}}\Phi_{\mathfrak{a}}(x)\,\big{)}^{2}\,\big{]}$$

$$=\frac{1}{(2\pi)^{m}}\int_{{\mathbb{R}}^{m}}\xi_{i}^{2}\xi_{j}^{2}\mathfrak{a}\big{(}\,|\xi|^{2}\,\big{)}d\xi\stackrel{{\scriptstyle(\ref{spec_mom_2})}}{{=}}h_{m}.$$

If $i=j$, then the above integral is nonzero iff $k=\ell$, in which case we deduce from (2.10) and (2.11)

$${\mathbb{E}}\big{[}\,\partial^{2}_{x_{i}}\Phi_{\mathfrak{a}}(x)\partial^{2}_{x_{k}}\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}=\begin{cases}h_{m}&i\neq k,\\ 3h_{m},&i=k.\end{cases}$$

The above equalities can be rewritten in the more compact form

$${\mathbb{E}}\big{[}\,L_{ij}({\boldsymbol{x}})L_{k\ell}({\boldsymbol{x}})\,\big{]}=h_{m}\big{(}\,\delta_{ij}\delta_{k\ell}+\delta_{ik}\delta_{j\ell}+\delta_{i\ell}\delta_{jk}\,\big{)},\;\;\forall i\leq j,\;k\leq\ell.$$ (2.15)

These equalities show that the off diagonal entries of $\operatorname{Hess}_{\Phi}$ are $\mathrm{i.i.d.}$, and also independent of the diagonal entries. The diagonal entries have identical distributions but are dependent. The parameter $h_{m}$ describes the various variances and covariances.

The Gaussian measure on $\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})$ determined by these covariance equalities is invariant with respect to the action of $O(m)$ by conjugation on $\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})$. For $v>0$ denote by $\mathcal{S}_{m}^{v}$ the space $\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})$ equipped with the $O(m)$-invariant Gaussian measure on $\operatorname{\mathbf{Sym}}({\mathbb{R}}^{m})$ determined by the covariances

$${\mathbb{E}}\big{[}\,L_{ij}({\boldsymbol{x}})L_{k\ell}({\boldsymbol{x}})\,\big{]}=v\big{(}\,\delta_{ij}\delta_{k\ell}+\delta_{ik}\delta_{j\ell}+\delta_{i\ell}\delta_{jk}\,\big{)},\;\;\forall i\leq j,\;k\leq\ell.$$

For any box $B\subset{\mathbb{R}}^{m}$ we denote by $Z_{\mathfrak{a}}(B)=Z_{\mathfrak{a},m}(B)$ the number of critical points of $\Phi_{\mathfrak{a}}$ in $B$. We deduce from the Kac-Rice formula (2.1) that

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B)\,\big{]}=\int_{B}{\mathbb{E}}_{\mathcal{S}_{m}^{h_{m}}}\big{[}\,|\det H|\,\big{]}p_{\nabla\Phi({\boldsymbol{x}})}(0){\boldsymbol{\lambda}}\big{[}\,d{\boldsymbol{x}}\,\big{]}$$

$$\stackrel{{\scriptstyle(\ref{grad_iso})}}{{=}}(2\pi d_{m})^{-m/2}{\mathbb{E}}_{\mathcal{S}_{m}^{h_{m}}}\big{[}\,|\det H|\,\big{]}\operatorname{vol}\big{[}\,B\,\big{]}$$

($X=(2h_{m})^{-1/2}H$)

$$=\underbrace{\left(\frac{h_{m}}{\pi d_{m}}\right)^{m/2}{\mathbb{E}}_{\mathcal{S}_{m}^{1/2}}\big{[}\,|\det X|\,\big{]}}_{=:C_{m}(\mathfrak{a})}\,\operatorname{vol}\big{[}\,B\,\big{]}.$$

Hence

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a},m}(B)\,\big{]}=C_{m}(\mathfrak{a})\operatorname{vol}\big{[}\,B\,\big{]}.$$ (2.16)

Remark 2.3.

One can prove that as $m\to\infty$

$$C_{m}(\mathfrak{a})\sim 2^{\frac{5}{2}}\Gamma\left(\frac{m+3}{2}\right)\left(\frac{h_{m}}{\pi d_{m}}\right)^{m/2}\left(\frac{1}{m+1}\right)^{1/2}.$$

Using (2.9) and (2.10) we deduce

$$\frac{h_{m}}{d_{m}}=\frac{\Gamma(1+m/2)}{\Gamma(2+m/2)}\times\frac{I_{m+3}(\mathfrak{a})}{I_{m+1}(\mathfrak{a})}=\frac{2I_{m+3}(\mathfrak{a})}{(m+2)I_{m+1}(\mathfrak{a})}.$$

Hence

$$\begin{split}C_{m}(\mathfrak{a})\sim 2^{5/2}\left(\frac{h_{m}(\mathfrak{a})}{d_{m}(\mathfrak{a})}\right)^{m/2}\Gamma\left(\frac{m+3}{2}\right)m^{-1/2}\hskip 72.26999pt\\ \sim 2^{\frac{m+5}{2}}\left(\frac{I_{m+3}(\mathfrak{a})}{(m+2)I_{m+1}(\mathfrak{a})}\right)^{m/2}\Gamma\left(\frac{m+3}{2}\right)m^{-1/2}\;\;\mbox{as $m\to\infty$}.\end{split}$$ (2.17)

The constant $C_{m}(\mathfrak{a})$ tends to grow very fast as $m\to\infty$, but its large $m$ behavior depends on the tail of the amplitude $\mathfrak{a}$. Roughly speaking, the slower the decay at $\infty$ of $\mathfrak{a}$ the faster the growth of $C_{m}(\mathfrak{a})$. $\hbox to0.0pt{$\sqcap$\hss}\sqcup$

3. An abstract law of large numbers

Fix $m\in{{\mathbb{N}}}$. For $N\in{{\mathbb{N}}}$, let ${\mathbb{I}}_{N}:=[1,N]\cap{{\mathbb{N}}}$. Suppose that we have an even continuous function $\rho:{\mathbb{R}}^{m}\to(0,\infty)$ that decays sufficiently fast to $0$ as $|x|\to\infty$. Then

$$\int_{NB\times NB}\rho({\boldsymbol{x}}-{\boldsymbol{y}})d{\boldsymbol{x}}d{\boldsymbol{y}}=N^{2m}\int_{B\times B}\rho_{N}\big{(}\,{\boldsymbol{u}}-{\boldsymbol{v}}\,\big{)}d{\boldsymbol{u}}d{\boldsymbol{v}}$$

where $\rho_{N}({\boldsymbol{x}})=\rho(N{\boldsymbol{x}})$. Observing that $\rho_{N}({\boldsymbol{x}})\to 0$ almost everywhere on $B$ we deduce from the dominated convergence theorem that

$$\int_{B\times B}\rho_{N}\big{(}\,{\boldsymbol{u}}-{\boldsymbol{v}}\,\big{)}d{\boldsymbol{u}}d{\boldsymbol{v}}\to 0$$

as $N\to\infty$. Hence

$$\int_{NB\times NB}\rho({\boldsymbol{x}}-{\boldsymbol{y}})d{\boldsymbol{x}}d{\boldsymbol{y}}=o\big{(}\,N^{2m}\,\big{)}\;\;\mbox{as $N\to\infty$}$$

The sum

$$\sum_{({\vec{k}},{\vec{\ell}})\in{\mathbb{I}}_{N}^{m}\times{\mathbb{I}}_{N}^{m}}\rho({\vec{k}}-{\vec{\ell}})$$

is a very rough Riemann sum approximation of the above integral when $B=[-1]^{N}$. The next result shows that under some conditions on the decay at infinity, this Riemann sum is also $o\big{(}\,N^{2m}\,\big{)}$ as $N\to\infty$.

Lemma 3.1.

Fix $\alpha,p>0$. For any $N\in{{\mathbb{N}}}$, $\alpha>0$ we set $\mathcal{R}_{N}:={\mathbb{I}}_{N}^{m}\times{\mathbb{I}}_{N}^{m}$. There exists a constant $K=K(\alpha,p,m)>0$ such that

$$\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}\leq K\begin{cases}N^{2m-\kappa(p)},\;\;\kappa(p)=\min(1,p),&p\neq m,\\ N^{2m-1}\log N,&p=m.\end{cases}$$

Proof.

For any $N\in{{\mathbb{N}}}$ define

$$D_{N}:=\big{\{}\,({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N};\;\;\exists j=1,\dotsc,m,\;\;k_{j}=\ell_{j}\,\big{\}},\;\;\mathcal{R}_{N}^{*}:=\mathcal{R}_{N}\setminus D_{N}.$$

Note that

$$D_{N}=\bigcup_{i=1}^{m}D_{N}^{i},\;\;D_{N}^{i}=\big{\{}\,({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N};\;\;k_{i}=\ell_{i}\,\big{\}}.$$

For $1\leq i_{1}<\cdots<i_{r}\leq m$ we have

$$\#\big{(}\,D_{N}^{i_{1}}\cap\cdots\cap D_{N}^{i_{p}}\,\big{)}=N^{2m-2r+1}$$

Using the Inclusion-Exclusion Principle we deduce that

$$\#D_{N}=\sum_{r=1}^{m}(-1)^{p-1}\binom{m}{p}N^{2m-2r+1}\leq 2^{m}N^{2m-1}.$$

We have

$$\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}=\underbrace{\sum_{({\vec{k}},{\vec{\ell}})\in D_{N}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}}_{=:Y_{N}}+\underbrace{\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N}^{*}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}}_{=:Z_{N}}.$$

Note that

$$Y_{N}\leq\#D_{N}\leq 2^{m}N^{2m-1}.$$

To estimate $Z_{N}$ we need to first analyze the structure of the region $\mathcal{R}_{N}^{*}$. Denote by $R_{i}$ the reflection

$$R_{i}:{\mathbb{R}}^{m}\times{\mathbb{R}}^{m}\to{\mathbb{R}}^{m}\times{\mathbb{R}}^{m},\;\;\left[\begin{array}[]{cc}x_{1}&y_{1}\\ \vdots&\vdots\\ \mathchoice{\leavevmode\hbox to13pt{\vbox to13pt{\pgfpicture\makeatletter\hbox{\hskip 6.49834pt\lower-6.49834pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{6.29834pt}{0.0pt}\pgfsys@curveto{6.29834pt}{3.47852pt}{3.47852pt}{6.29834pt}{0.0pt}{6.29834pt}\pgfsys@curveto{-3.47852pt}{6.29834pt}{-6.29834pt}{3.47852pt}{-6.29834pt}{0.0pt}\pgfsys@curveto{-6.29834pt}{-3.47852pt}{-3.47852pt}{-6.29834pt}{0.0pt}{-6.29834pt}\pgfsys@curveto{3.47852pt}{-6.29834pt}{6.29834pt}{-3.47852pt}{6.29834pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-3.82227pt}{-1.22945pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\displaystyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}{\leavevmode\hbox to13pt{\vbox to13pt{\pgfpicture\makeatletter\hbox{\hskip 6.49834pt\lower-6.49834pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{6.29834pt}{0.0pt}\pgfsys@curveto{6.29834pt}{3.47852pt}{3.47852pt}{6.29834pt}{0.0pt}{6.29834pt}\pgfsys@curveto{-3.47852pt}{6.29834pt}{-6.29834pt}{3.47852pt}{-6.29834pt}{0.0pt}\pgfsys@curveto{-6.29834pt}{-3.47852pt}{-3.47852pt}{-6.29834pt}{0.0pt}{-6.29834pt}\pgfsys@curveto{3.47852pt}{-6.29834pt}{6.29834pt}{-3.47852pt}{6.29834pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-3.82227pt}{-1.22945pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\textstyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}{\leavevmode\hbox to10.09pt{\vbox to10.09pt{\pgfpicture\makeatletter\hbox{\hskip 5.04326pt\lower-5.04326pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{4.84326pt}{0.0pt}\pgfsys@curveto{4.84326pt}{2.67488pt}{2.67488pt}{4.84326pt}{0.0pt}{4.84326pt}\pgfsys@curveto{-2.67488pt}{4.84326pt}{-4.84326pt}{2.67488pt}{-4.84326pt}{0.0pt}\pgfsys@curveto{-4.84326pt}{-2.67488pt}{-2.67488pt}{-4.84326pt}{0.0pt}{-4.84326pt}\pgfsys@curveto{2.67488pt}{-4.84326pt}{4.84326pt}{-2.67488pt}{4.84326pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-2.68936pt}{-0.84743pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\scriptstyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}{\leavevmode\hbox to8.67pt{\vbox to8.67pt{\pgfpicture\makeatletter\hbox{\hskip 4.3355pt\lower-4.3355pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{4.1355pt}{0.0pt}\pgfsys@curveto{4.1355pt}{2.284pt}{2.284pt}{4.1355pt}{0.0pt}{4.1355pt}\pgfsys@curveto{-2.284pt}{4.1355pt}{-4.1355pt}{2.284pt}{-4.1355pt}{0.0pt}\pgfsys@curveto{-4.1355pt}{-2.284pt}{-2.284pt}{-4.1355pt}{0.0pt}{-4.1355pt}\pgfsys@curveto{2.284pt}{-4.1355pt}{4.1355pt}{-2.284pt}{4.1355pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-2.11784pt}{-0.41687pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\scriptscriptstyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}&\boxed{y_{i}}\\ \vdots&\vdots\\ x_{m}&y_{m}\end{array}\right]\mapsto\left[\begin{array}[]{cc}x_{1}&y_{1}\\ \vdots&\vdots\\ \boxed{y_{i}}&\mathchoice{\leavevmode\hbox to13pt{\vbox to13pt{\pgfpicture\makeatletter\hbox{\hskip 6.49834pt\lower-6.49834pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{6.29834pt}{0.0pt}\pgfsys@curveto{6.29834pt}{3.47852pt}{3.47852pt}{6.29834pt}{0.0pt}{6.29834pt}\pgfsys@curveto{-3.47852pt}{6.29834pt}{-6.29834pt}{3.47852pt}{-6.29834pt}{0.0pt}\pgfsys@curveto{-6.29834pt}{-3.47852pt}{-3.47852pt}{-6.29834pt}{0.0pt}{-6.29834pt}\pgfsys@curveto{3.47852pt}{-6.29834pt}{6.29834pt}{-3.47852pt}{6.29834pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-3.82227pt}{-1.22945pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\displaystyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}{\leavevmode\hbox to13pt{\vbox to13pt{\pgfpicture\makeatletter\hbox{\hskip 6.49834pt\lower-6.49834pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{6.29834pt}{0.0pt}\pgfsys@curveto{6.29834pt}{3.47852pt}{3.47852pt}{6.29834pt}{0.0pt}{6.29834pt}\pgfsys@curveto{-3.47852pt}{6.29834pt}{-6.29834pt}{3.47852pt}{-6.29834pt}{0.0pt}\pgfsys@curveto{-6.29834pt}{-3.47852pt}{-3.47852pt}{-6.29834pt}{0.0pt}{-6.29834pt}\pgfsys@curveto{3.47852pt}{-6.29834pt}{6.29834pt}{-3.47852pt}{6.29834pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-3.82227pt}{-1.22945pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\textstyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}{\leavevmode\hbox to10.09pt{\vbox to10.09pt{\pgfpicture\makeatletter\hbox{\hskip 5.04326pt\lower-5.04326pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{4.84326pt}{0.0pt}\pgfsys@curveto{4.84326pt}{2.67488pt}{2.67488pt}{4.84326pt}{0.0pt}{4.84326pt}\pgfsys@curveto{-2.67488pt}{4.84326pt}{-4.84326pt}{2.67488pt}{-4.84326pt}{0.0pt}\pgfsys@curveto{-4.84326pt}{-2.67488pt}{-2.67488pt}{-4.84326pt}{0.0pt}{-4.84326pt}\pgfsys@curveto{2.67488pt}{-4.84326pt}{4.84326pt}{-2.67488pt}{4.84326pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-2.68936pt}{-0.84743pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\scriptstyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}{\leavevmode\hbox to8.67pt{\vbox to8.67pt{\pgfpicture\makeatletter\hbox{\hskip 4.3355pt\lower-4.3355pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}}\hbox{\hbox{{\pgfsys@beginscope\pgfsys@invoke{ }{{}{{{}}}{{}}{}{}{}{}{}{}{}{}{}{{}\pgfsys@moveto{4.1355pt}{0.0pt}\pgfsys@curveto{4.1355pt}{2.284pt}{2.284pt}{4.1355pt}{0.0pt}{4.1355pt}\pgfsys@curveto{-2.284pt}{4.1355pt}{-4.1355pt}{2.284pt}{-4.1355pt}{0.0pt}\pgfsys@curveto{-4.1355pt}{-2.284pt}{-2.284pt}{-4.1355pt}{0.0pt}{-4.1355pt}\pgfsys@curveto{2.284pt}{-4.1355pt}{4.1355pt}{-2.284pt}{4.1355pt}{0.0pt}\pgfsys@closepath\pgfsys@moveto{0.0pt}{0.0pt}\pgfsys@stroke\pgfsys@invoke{ } }{{{{}}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@transformcm{1.0}{0.0}{0.0}{1.0}{-2.11784pt}{-0.41687pt}\pgfsys@invoke{ }\hbox{{\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\hbox{{$\scriptscriptstyle x_{i}$}} }}\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope}}} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{{{}}}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\\ \vdots&\vdots\\ x_{m}&y_{m}\end{array}\right]. \textstyle x_{i} \textstyle x_{i} \textstyle x_{i} \textstyle x_{i}$$

Denote by $G_{m}$ the direct product of cyclic groups

$$G_{m}=\big{(}\,{\mathbb{Z}}/2{\mathbb{Z}}\,\big{)}^{m}=\big{\{}\,\vec{{\epsilon}}=({\epsilon}_{1},\dotsc,{\epsilon}_{m});\;\;{\epsilon}_{k}=0,1\,\big{\}}.$$

The group $G_{m}$ acts freely on $\mathcal{R}_{N}^{*}$

$$\vec{{\epsilon}}\cdot\big{(}\,{\vec{k}},{\vec{\ell}}\,\big{)}=R^{\vec{{\epsilon}}}\big{(}\,{\vec{k}},{\vec{\ell}}\,\big{)},\;\;R^{\vec{{\epsilon}}}=R_{1}^{{\epsilon}_{1}}\cdots R_{m}^{{\epsilon}_{m}}.$$

We denote by $\mathcal{C}_{N}^{+}$ the positive chamber of $\mathcal{R}_{N}^{*}$,

$$\mathcal{C}_{n}^{+}:=\big{\{}\,\big{(}\,({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N}^{*};\;\ell_{j}>k_{j},\;\;\forall 1\leq j\leq m\,\big{\}},$$

and we observe that

$$\mathcal{C}_{n}^{+}=\mathcal{T}_{N}^{m},\;\;\mathcal{T}_{N}:=\big{\{}\,(k,\ell)\in{\mathbb{I}}_{N}\times{\mathbb{I}}_{N};\;\;\ell>k\,\big{\}}.$$

We have

$$\mathcal{R}_{N}^{*}=\bigcup_{{\vec{{\epsilon}}}\in G_{m}}R^{\vec{{\epsilon}}}\mathcal{C}_{n}^{+}$$

The function $\rho$ is $G_{m}$-invariant so

$$\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N}^{*}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}=\sum_{{\vec{{\epsilon}}}\in G_{m}}\sum_{({\vec{k}},{\vec{\ell}})\in R^{\vec{{\epsilon}}}\mathcal{C}_{N}^{+}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}$$

$$=2^{m}\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{C}_{N}^{+}}\frac{1}{\big{(}\,1+\alpha|{\vec{\ell}}-{\vec{k}}|_{1}\,\big{)}^{p}}<\frac{2^{m}}{\alpha^{p}}\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{C}_{N}^{+}}\frac{1}{|{\vec{\ell}}-{\vec{k}}|^{p}_{1}}$$

(use the AM-GM inequality)

$$\leq\frac{2^{m}}{(m\alpha)^{p}}\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{C}_{N}^{+}}\Big{(}\,\prod_{j=1}^{m}(\ell_{j}-k_{j})\,\Big{)}^{-p/m}.$$

($\mathcal{C}_{n}^{+}=\mathcal{T}_{N}^{m}$)

$$=\frac{2^{m}}{(m\alpha)^{p}}\prod_{j=1}^{m}\sum_{(k_{j},\ell_{j})\in\mathcal{T}_{N}}\big{(}\,\ell_{j}-k_{j}\,\big{)}^{-p/m}=\frac{2^{m}}{(m\alpha)^{p}}\Big{(}\,\sum_{(k,\ell)\in\mathcal{T}_{N}}\big{(}\,\ell-k\,\big{)}^{-p/m}\,\Big{)}^{m}$$

. Now observe that

$$\sum_{(k,\ell)\in\mathcal{T}_{N}}\big{(}\,\ell-k\,\big{)}^{-p/m}=\sum_{k=1}^{N-1}\sum_{j=1}^{N-k}j^{-p/m}.$$

To proceed further we need to use the following result.

Sublemma 3.2.

Let $r\in{\mathbb{R}}$. Then for any $M\in{{\mathbb{N}}}$

$$S_{r}(M):=\sum_{j=1}^{M}j^{r}\leq u_{r}(M):=\begin{cases}\frac{1}{r+1}M^{r+1},&r\geq 0,\\ \frac{1}{r+1}M^{r+1}+1,&r\in(-1,0),\\ 1+\log M&r=-1\\ 2,&r<-1.\end{cases}$$

Proof.

Using approximations by Riemann sums for the integral

$$I_{r}(M):=\int_{1}^{M}x^{r}dx$$

we deduce

$$S_{r}(M)\leq\begin{cases}I_{r}(M),&r>0,\\ I_{r}(M)+1,&r<0\end{cases}=\begin{cases}\frac{1}{r+1}\big{(}\,M^{r+1}-1\,\big{)},&r\geq 0,\\ \frac{1}{r+1}\big{(}\,M^{r+1}-1\,\big{)}+1,&r\in(-1,0),\\ 1+\log M+1,&r=-1,\\ \frac{1}{|r+1|}\big{(}\,1-M^{r+1}\,\big{)}+1,&r<-1.\end{cases}$$

$$\leq\begin{cases}\frac{1}{r+1}M^{r+1},&r\geq 0,\\ \frac{1}{r+1}M^{r+1}+1,&r\in(-1,0),\\ 1+\log M,&r=-1,\\ 2,&r<-1.\end{cases}$$

$\hbox to0.0pt{$\sqcap$\hss}\sqcup$

Suppose that $p\neq m$. Using Sublemma 3.2 we deduce that

$$\sum_{j=1}^{N-k}j^{-p/m}\leq u_{-p/m}(N-k)=\begin{cases}\frac{1}{(1-p/m)}(N-k)^{1-p/m},&p/m<1,\\ 2&1<p/m.\end{cases}$$

Next, using the sublemma again we deduce

$$\sum_{k=1}^{N}u_{r}(N-k)\leq\begin{cases}\frac{1}{(1-p/m)(2-p/m)}N^{2-p/m},&p/m<1,\\ 2N&1<p/m.\end{cases}$$

Hence

$$\sum_{(k,\ell)\in\mathcal{T}_{N}}\big{(}\,\ell-k\,\big{)}^{-p/m}\leq\begin{cases}\frac{1}{(1-p/m)(2-p/m)}N^{2-p/m},&p/m<1,\\ 2N,&1<p/m.\end{cases}$$

and thus

$$Z_{N}=\sum_{({\vec{k}},{\vec{\ell}})\in\mathcal{R}_{N}^{*}}\frac{1}{(1+|{\vec{\ell}}-{\vec{k}}|_{1})^{p}}\leq\frac{2^{m}}{(m\alpha)^{p}}\leq C(m,\alpha,p)\begin{cases}N^{2m-p},&p<m,\\ N^{m}&p>m.\end{cases}$$

If $p=m$, then Sublemma 3.2 implies that

$$\sum_{j=1}^{N-k}j^{-1}\leq u_{-1}(N-k)=1+\log(N-k),$$

and

$$\sum_{k=1}^{N}\big{(}\,1+\log(N-k),\,\big{)}=N+\log N!.$$

The conclusion follows from Stirling’s formula which implies that

$$\log N!=O\big{(}\,N\log N\,\big{)}.$$

$\hbox to0.0pt{$\sqcap$\hss}\sqcup$

Corollary 3.3.

Consider a family of random variable $\big{(}\,X_{\vec{\ell}}\,\big{)}_{{\vec{\ell}}\in{{\mathbb{N}}}^{M}}$ defined on the same probability space $(\Omega,\mathcal{S},{{\mathbb{P}}})$ such that there exist constants $C,\alpha,p>0$ such that

$$\big{|}\,\operatorname{Cov}\big{[}\,X_{\vec{k}},X_{\vec{\ell}}\,\big{]}\,\big{|}\leq\frac{C}{\big{(}\,1+\alpha|{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}\,\big{)}^{p}},\;\;\forall{\vec{k}},{\vec{\ell}}\in{{\mathbb{N}}}^{m}.$$

Then, as $N\to\infty$

$$S_{N}(X):=\frac{1}{N^{m}}\sum_{{\vec{k}}\in{\mathbb{I}}_{N}^{m}}\big{(}\,X_{\vec{k}}-{\mathbb{E}}\big{[}\,X_{\vec{k}}\,\big{]}\,\big{)}\to 0$$

in $L^{2}$. In particular $S_{N}(X)$ converges in probability to zero as $N\to\infty$.

Proof.

We have

$${\mathbb{E}}\big{[}\,S_{N}(X)^{2}\,\big{]}=\frac{1}{N^{2m}}\sum_{({\vec{k}},{\vec{\ell}})\in{\mathbb{I}}_{N}^{m}\times{\mathbb{I}}_{N}}\operatorname{Cov}\big{[}\,X_{\vec{k}},X_{\vec{\ell}}\,\big{]}.$$

Lemma 3.1 implies that $S_{N}(X)\to 0$ in $L^{2}$. $\hbox to0.0pt{$\sqcap$\hss}\sqcup$

4. Proof of Theorem 1.1

We have

$$Z_{\mathfrak{a}}(N\cdot B)=\sum_{{\vec{\ell}}\in{\mathbb{I}}_{N}^{m}}Z_{\mathfrak{a}}(B_{\vec{\ell}})$$

since the probability that $\Phi_{\mathfrak{a}}$ has a critical point on the boundary of some $B_{\ell}$ is zero.

The Gaussian function $\Phi_{\mathfrak{a}}$ is a stationary is stationary, so the collection $\big{(}\,Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{)}_{{\vec{\ell}}\in{{\mathbb{N}}}^{m}}$ is stationary as well, i.e., for any $n\in{{\mathbb{N}}}$ and any ${\vec{\ell}},{\vec{\ell}}_{1},\dotsc,{\vec{\ell}}_{n}\in{{\mathbb{N}}}^{n}$ the random vectors

$$\big{(}\,Z_{\mathfrak{a}}(B_{{\vec{\ell}}_{1}}),\dotsc,Z_{\mathfrak{a}}(B_{{\vec{\ell}}_{n}})\,\big{)}\;\mbox{and}\;\big{(}\,Z_{\mathfrak{a}}(B_{{\vec{\ell}}_{1}+{\vec{\ell}}}),\dotsc,Z_{\mathfrak{a}}(B_{{\vec{\ell}}_{n}+{\vec{\ell}}})\,\big{)}$$

have the same distribution. We deduce from the ergodic theorem [6, Chap.6, Thm.2.5] that the averages

$$\frac{1}{\operatorname{vol}(NB)}Z_{\mathfrak{a}}(N\cdot B)=\frac{1}{N^{m}\operatorname{vol}(B)}\sum_{{\vec{\ell}}\in{\mathbb{I}}_{N}}Z_{\mathfrak{a}}(B_{\vec{\ell}})$$

converge $\mathrm{a.s.}$ and in $L^{1}$ to some random variable $Z_{\infty}$.

Set

$$X_{\vec{\ell}}:=Z_{\mathfrak{a}}(B_{\vec{\ell}})-{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{]}.$$

We need to prove that the limit $Z_{\infty}$ is the constant

$$\frac{1}{\operatorname{vol}(B)}{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B)\,\big{]}.$$

To achieve this

$$\lim_{N\to\infty}\frac{1}{N^{m}}\Bigg{\|}\sum_{{\vec{\ell}}\in{\mathbb{I}}_{N}^{m}}X_{\vec{\ell}}\Bigg{\|}_{L^{2}}\to 0.$$

We will deduce this from Corollary 3.3. We set

$$C({\vec{k}},{\vec{\ell}}):=\operatorname{Cov}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{k}}),Z_{\mathfrak{a}}(B_{\vec{\ell}})\big{[}\,={\mathbb{E}}\big{[}\,X_{\vec{k}}X_{\vec{\ell}}\,\big{]},\;\;\forall{\vec{k}},{\vec{\ell}}\in{\mathbb{I}}_{N}^{m}.$$

Notice that since $\Phi_{\mathfrak{a}}$ is a stationary Gaussian function we have

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\ell})^{p}\,\big{]}={\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B)^{p}\,\big{]}<\infty,\;\;\forall p\in[1,\infty),\;\;{\vec{\ell}}\in{{\mathbb{N}}}^{m}$$

and we deduce that

$$\exists K_{1}=K_{1}(\mathfrak{a},m)>0:\big{|}\,C({\vec{k}},{\vec{\ell}})\,\big{|}<K_{1},\;\;\forall{\vec{k}},{\vec{\ell}}\in{{\mathbb{N}}}^{m}.$$ (4.1)

To proceed further we define

$${\widehat{\Phi}}:{\mathbb{R}}^{m}\times{\mathbb{R}}^{m}\to{\mathbb{R}},\;\;{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})=\Phi_{\mathfrak{a}}({\boldsymbol{x}})+\Phi_{\mathfrak{a}}({\boldsymbol{y}})$$

and we set

$$\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}}):=\operatorname{Hess}_{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}}),\;\;H({\boldsymbol{x}}):=\operatorname{Hess}_{\Phi_{\mathfrak{a}}}({\boldsymbol{x}}).$$

Denote for any Borel subset $\widehat{B}\in{\mathbb{R}}^{m}\times{\mathbb{R}}^{m}$ by $\widehat{Z}\big{(}\,\widehat{B}\,\big{)}$ the number of critical points of ${\widehat{\Phi}}$ in $\widehat{B}$. Note that if $B_{\vec{k}}\cap B_{\vec{\ell}}=\emptyset$, then

$$\widehat{Z}(B_{\vec{k}}\times B_{\vec{\ell}})=Z_{\mathfrak{a}}(B_{\vec{k}})Z_{\mathfrak{a}}(B_{\vec{\ell}})$$

so

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{k}})Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{]}={\mathbb{E}}\big{[}\,\widehat{Z}(B_{\vec{k}}\times B_{\vec{\ell}})\,\big{]}.$$

We want to compute ${\mathbb{E}}\big{[}\,\widehat{Z}(B_{\vec{k}}\times B_{\vec{\ell}})\,\big{]}$ using the Kac-Rice formula. We first need to verify that $\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})$ is nondegenerate for any ${\boldsymbol{x}}\neq{\boldsymbol{y}}$.

Lemma 4.1.

For any ${\boldsymbol{x}},{\boldsymbol{y}}\in{\mathbb{R}}^{m}$, ${\boldsymbol{x}}\neq{\boldsymbol{y}}$, the Gaussian vector $\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})$ is nondegenerate.

Proof.

We have

$$\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=\left[\begin{array}[]{cc}\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}&\operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{]}\\ &\\ \operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}&\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{]}\end{array}\right].$$

As shown in(2.12), for any ${\boldsymbol{x}}\in{\mathbb{R}}^{n}$ we have

$$\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}=v_{m}{\mathbbm{1}}_{m},\;\;v_{m}=\int_{{\mathbb{R}}^{n}}\xi_{1}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}.$$

We have

$$\operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{]}=\big{(}\,\partial_{x_{j}}\partial_{y_{k}}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}}-{\boldsymbol{y}})\,\big{)}_{1\leq j,k\leq m}$$

and

$$\big{[}\,\partial_{x_{j}}\partial_{y_{k}}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}}-{\boldsymbol{y}})=\int_{{\mathbb{R}}^{m}}e^{-\boldsymbol{i}\langle\xi,{\boldsymbol{x}},{\boldsymbol{y}}\rangle}\xi_{j}\xi_{k}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}.$$ (4.2)

Since $\Phi_{\mathfrak{a}}$ is stationary it suffice to consider only the case ${\boldsymbol{x}}=0$. On the other hand, $\Phi_{\mathfrak{a}}$ is $O(m)$-invariant so, up to a rotation we can assume that ${\boldsymbol{x}}-{\boldsymbol{y}}=-t{\boldsymbol{e}}_{1}$, $t\neq 0$, where $\{{\boldsymbol{e}}_{1},\dotsc,{\boldsymbol{e}}_{m}\}$ is the canonical basis of ${\mathbb{R}}^{m}$. Hence

$$\partial_{x_{j}}\partial_{y_{k}}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}}-{\boldsymbol{y}})=\int_{{\mathbb{R}}^{m}}e^{\boldsymbol{i}t\xi_{1}}\xi_{j}\xi_{k}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}.$$

Let us observe that if $j\neq k$, then either $j\neq 1$, or $k\neq 1$. Suppose $j\neq 1$. The function $e^{\boldsymbol{i}t\xi_{1}}\xi_{j}\xi_{k}$ is odd with respect to the reflection $\xi_{j}\mapsto-\xi_{j}$ so

$$\partial_{x_{j}}\partial_{y_{k}}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}},{\boldsymbol{y}})=\int_{{\mathbb{R}}^{m}}e^{\boldsymbol{i}t\xi_{1}}\xi_{j}\xi_{k}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}=0,\;\;\forall j\neq k.$$

If $j=k$, then

$$v_{m}(j):=\partial_{x_{j}}\partial_{y_{j}}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{x}},{\boldsymbol{y}})=\int_{{\mathbb{R}}^{m}}e^{\boldsymbol{i}t\xi_{1}}\xi_{j}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}=\int_{{\mathbb{R}}^{m}}\cos(t\xi_{1})\xi_{j}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}$$

and we deduce

$$|v_{m}(j)|\leq\int_{{\mathbb{R}}^{m}}\big{|}\,\cos(t\xi_{1})\,\big{|}\xi_{j}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}<\int_{{\mathbb{R}}^{m}}\xi_{j}^{2}\mu_{\mathfrak{a}}\big{[}\,d\xi\,\big{]}=v_{m}.$$

After a reordering

$$\begin{array}[]{c}\big{(}\,\partial_{x_{1}}\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\dotsc,\partial_{x_{m}}\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\partial_{y_{1}}\Phi_{\mathfrak{a}}({\boldsymbol{y}}),\dotsc\partial_{y_{m}}\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{)}\\ \downarrow\\ \big{(}\,\partial_{x_{1}}\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\partial_{y_{1}}\Phi_{\mathfrak{a}}({\boldsymbol{y}}),\dotsc,\partial_{x_{m}}\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\partial_{y_{m}}\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{)}\end{array}$$

we see that

$$\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=\bigoplus_{j=1}^{m}\underbrace{\left[\begin{array}[]{cc}v_{m}&v_{m}(j)\\ v_{m}(j)&v_{m}\end{array}\right]}_{=:V_{j}}.$$

Note that, for each $j$, the symmetric matrix $V_{j}$ is positive definite since

$$\det V_{j}=v_{m}^{2}-v_{m}(j)^{2}>0.$$

$\hbox to0.0pt{$\sqcap$\hss}\sqcup$

Suppose that $B_{\vec{k}}\cap B_{\vec{\ell}}=\emptyset$. We deduce from Lemma 4.1 that $\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})$, is nondegenerate for any $({\boldsymbol{x}},{\boldsymbol{y}})\in B_{\vec{k}},\times B_{\vec{\ell}}$. We can the apply the Kac-Rice formula to deduce that

$$\begin{split}{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{k}})Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{]}={\mathbb{E}}\big{[}\,\widehat{Z}(B_{\vec{k}}\times B_{\vec{\ell}})\,\big{]}\hskip 142.26378pt\\ =\int_{B_{\vec{k}}\times B_{\vec{\ell}}}\underbrace{{\mathbb{E}}\big{[}\,|\det\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})|\,\big{|}\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})=0\,\big{]}p_{\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})}(0)}_{=\widehat{\rho}({\boldsymbol{x}},{\boldsymbol{y}})}{\boldsymbol{\lambda}}\big{[}\,d{\boldsymbol{x}}d{\boldsymbol{y}}\,\big{]},\\ \;\;\mbox{if $B_{\vec{k}}\cap B_{\vec{\ell}}=\emptyset$}\hskip 199.16928pt.\end{split}$$ (4.3)

For ${\boldsymbol{x}}\in{\mathbb{R}}^{m}$ we denote by $\big{|}\,{\boldsymbol{x}}\,\big{|}_{\infty}$ the sup-norm of ${\boldsymbol{x}}$

$$\big{|}\,{\boldsymbol{x}}\,\big{|}_{\infty}:=\max_{1\leq i\leq m}|x_{i}|.$$

Note that

$$\big{|}\,{\boldsymbol{x}}\,\big{|}_{1}\leq m\big{|}\,{\boldsymbol{x}}\,\big{|}_{\infty}.$$

and $B_{\vec{k}}\cap B_{\vec{\ell}}=\emptyset$ if $\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{\infty}>1$. Hence,

$$\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}>m\Rightarrow B_{\vec{k}}\cap B_{\vec{\ell}}=\emptyset.$$

Hence

$$\begin{split}{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{v}k)Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{]}=&\int_{B_{\vec{k}}\times B_{\vec{\ell}}}\underbrace{{\mathbb{E}}\big{[}\,|\det\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})|\,\big{|}\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})=0\,\big{]}p_{\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})}(0)}_{=\widehat{\rho}({\boldsymbol{x}},{\boldsymbol{y}})}{\boldsymbol{\lambda}}\big{[}\,d{\boldsymbol{x}}d{\boldsymbol{y}}\,\big{]},\\ \;\;\mbox{if $\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}>m$}.&\end{split}$$ (4.4)

Let us now express ${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{k}})\,\big{]}{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{]}$ as an integral over $B_{\vec{k}}\times B_{\vec{\ell}}$. Choose an independent copy $\Phi^{*}_{\mathfrak{a}}$ of $\Phi_{\mathfrak{a}}$. We set

$${\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}}):=\Phi_{\mathfrak{a}}({\boldsymbol{x}})+\Phi^{*}_{\mathfrak{a}}({\boldsymbol{y}}),\;\;\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}}):=\operatorname{Hess}_{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}}).$$

Then

$${\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{k}})\,\big{]}{\mathbb{E}}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{\ell}})\,\big{]}={\mathbb{E}}\big{[}\,Z_{{\widetilde{\Phi}}}(B_{\vec{k}}\times B_{\vec{\ell}})\,\big{]}$$

$$\begin{split}=&\int_{B_{\vec{k}}\times B_{\vec{\ell}}}\underbrace{{\mathbb{E}}\big{[}\,\big{|}\,\det\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{|}\,\big{\|}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})=0\,\big{]}p_{\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})}(0)}_{=\widetilde{\rho}({\boldsymbol{x}},{\boldsymbol{y}})}{\boldsymbol{\lambda}}\big{[}\,d{\boldsymbol{x}}d{\boldsymbol{y}}\,\big{]},\\ \;\;\mbox{$\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}>m$}.&\end{split}$$

Thus

$$\operatorname{Cov}\big{[}\,Z_{\mathfrak{a}}(B_{\vec{k}}),Z_{\mathfrak{a}}(B_{\vec{\ell}})]=\int_{B_{k}\times B_{\ell}}\big{(}\,\widehat{\rho}({\boldsymbol{x}},{\boldsymbol{y}})-\widetilde{\rho}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{)}{\boldsymbol{\lambda}}\big{[}\,{\boldsymbol{x}}d{\boldsymbol{y}}.\,\big{]}$$ (4.5)

To proceed further, we need to describe technical results that seems to be part of the mathematical folklore.

Digression 4.2.

Suppose that ${\boldsymbol{V}}$ is an $m$-dimensional real Euclidean space with inner product $(-,-)$. Denote by $S_{1}({\boldsymbol{V}})$ the unit sphere in ${\boldsymbol{V}}$ and by $\operatorname{\mathbf{Sym}}({\boldsymbol{V}})$ the space of symmetric operators ${\boldsymbol{V}}\to{\boldsymbol{V}}$ and by $\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})$ the cone of nonnegative ones. For $A\in\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})$ we denote by $\Gamma_{A}$ the centered Gaussian measure with variance $A$.

The space $\operatorname{\mathbf{Sym}}({\boldsymbol{V}})$ is equipped with an inner product

$$\big{(}\,A,B\,\big{)}_{\mathrm{op}}=\operatorname{{\rm tr}}(AB),\;\;\forall A,B\in\operatorname{\mathbf{Sym}}({\boldsymbol{V}}).$$

We denote by $\|-\|_{\mathrm{op}}$ the associated norm.

We have a natural map $\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})\to\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})$, $A\mapsto A^{1/2}$. We will need the following result, [5, Prop.2.1]

$$\mu\big{\|}\,A^{1/2}-B^{1/2}\,\big{\|}_{\mathrm{op}}\leq\big{\|}\,A-B\,\big{\|}^{1/2}_{\mathrm{op}},\;\;\forall A,B\in\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}}),\;\;A^{1/2}+B^{1/2}\geq\mu{\mathbbm{1}}$$ (4.6)

This ends the digression. $\hbox to0.0pt{$\sqcap$\hss}\sqcup$

Lemma 4.3.

Fix $\mu_{0}>0$ and $A_{0}\in\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})$ such that $A_{0}^{1/2}\geq\mu_{0}{\mathbbm{1}}$. Suppose that $f:{\boldsymbol{V}}\to{\mathbb{R}}$ is a locally Lipschitz function that is homogeneous of degree $k\geq 1$. For $A\in\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})$ we set

$${\mathcal{I}}_{A}(f):=\int_{\boldsymbol{V}}f({\boldsymbol{v}})\Gamma_{A}\big{[}\,d{\boldsymbol{v}}\,\big{]},$$

where $\Gamma_{A}$ denotes the centered Gaussian measure on ${\boldsymbol{V}}$ with variance $A$. Then for any $R\geq\|A_{0}\|_{\mathrm{op}}$ there exists a constant $C=C(f,R,\mu_{0})>0$ such that, $\forall A\in\operatorname{\mathbf{Sym}}_{\geq 0}({\boldsymbol{V}})$, $\|A\|_{\mathrm{op}}\leq R$

$$\big{|}\,{\mathcal{I}}_{A_{0}}(f)-{\mathcal{I}}_{A}(f)\,\big{|}\leq C\|A-A_{0}\|^{1/2}.$$ (4.7)

Proof.

The function $f$ is Lipschitz on the ball

$$B_{R}({\boldsymbol{V}}):=\big{\{}\,{\boldsymbol{v}}\in{\boldsymbol{V}};\|{\boldsymbol{v}}\|\leq R\,\big{\}},$$

so there exists $L=L(R)>0$ such that

$$\big{[}\,f({\boldsymbol{u}})-f({\boldsymbol{v}})\,\big{|}\leq L\|{\boldsymbol{u}}-{\boldsymbol{v}}\|,\;\;\forall{\boldsymbol{u}},{\boldsymbol{v}}\in B_{R}({\boldsymbol{V}}).$$ (4.8)

Note that

$${\mathcal{I}}_{A}(f)=\int_{\boldsymbol{V}}f\big{(}\,A^{1/2}{\boldsymbol{v}}\,\big{)}\Gamma_{{\mathbbm{1}}}\big{[}\,d{\boldsymbol{v}}\,\big{]},$$

so

$$\big{|}\,{\mathcal{I}}_{A_{0}}(f)-{\mathcal{I}}_{A}(f)\,\big{|}\leq\int_{\boldsymbol{V}}\big{|}\,f\big{(}\,A^{1/2}{\boldsymbol{v}}\,\big{)}-f\big{(}\,A_{0}^{1/2}{\boldsymbol{v}}\,\big{)}\,\big{|}\;\Gamma_{{\mathbbm{1}}}\big{[}\,d{\boldsymbol{v}}\,\big{]}$$

$$=\underbrace{\frac{1}{(2\pi)^{m/2}}\left(\int_{0}^{\infty}r^{n+k-1}e^{-r^{2}/2}dr\right)}_{C_{m,k}}\;\int_{S_{1}({\boldsymbol{V}})}\big{|}\,f\big{(}\,A^{1/2}{\boldsymbol{v}}\,\big{)}-f\big{(}\,A_{0}^{1/2}{\boldsymbol{v}}\,\big{)}\,\big{|}\operatorname{vol}_{S_{1}({\boldsymbol{V}})}\big{[}\,d{\boldsymbol{v}}\,\big{]}$$

$$\stackrel{{\scriptstyle(\ref{Lipschitz_gauss})}}{{\leq}}C_{m,k}L(R)\int_{S_{1}({\boldsymbol{V}})}\|A^{1/2}-A_{0}^{1/2}\|_{\mathrm{op}}\operatorname{vol}_{S^{1}({\boldsymbol{V}})}\big{[}\,d{\boldsymbol{v}}\,\big{]}\stackrel{{\scriptstyle(\ref{sq_holder})}}{{\leq}}C(k,R,\mu_{0})\|A-A_{0}\|_{\mathrm{op}}^{1/2}.$$

$\hbox to0.0pt{$\sqcap$\hss}\sqcup$

Note that

$$\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=\left[\begin{array}[]{cc}\operatorname{Var}\big{[}\,H({\boldsymbol{x}})\,\big{]}&0\\ 0&\operatorname{Var}\big{[}\,H({\boldsymbol{y}})\,\big{]}\end{array}\right].$$

Note that since $\Phi_{\mathfrak{a}}$ is stationary, $\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}$ is independent of ${\boldsymbol{x}}$ and ${\boldsymbol{y}}$. We have

$$\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=\left[\begin{array}[]{cc}\operatorname{Var}\big{[}\,H({\boldsymbol{x}})\,\big{]}&\operatorname{Cov}\big{[}\,H({\boldsymbol{x}}),H({\boldsymbol{y}})\,\big{]}\\ &\\ \operatorname{Cov}\big{[}\,H({\boldsymbol{y}}),H({\boldsymbol{x}})\,\big{]}&\operatorname{Var}\big{[}\,H({\boldsymbol{y}})\,\big{]}\end{array}\right]$$

$$=\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}+\underbrace{\left[\begin{array}[]{cc}0&\operatorname{Cov}\big{[}\,H({\boldsymbol{x}}),H({\boldsymbol{y}})\,\big{]}\\ &\\ \operatorname{Cov}\big{[}\,H({\boldsymbol{y}}),H({\boldsymbol{x}})\,\big{]}&0\end{array}\right]}_{=:R_{H}({\boldsymbol{x}},{\boldsymbol{y}})}.$$

For every ${\boldsymbol{z}}\in{\mathbb{R}}^{m}$ we set

$$T(z):=\sum_{|\alpha|\leq 4}\big{|}\,\partial^{\alpha}{\boldsymbol{K}}_{\mathfrak{a}}({\boldsymbol{z}})|.$$

Since ${\boldsymbol{K}}_{\mathfrak{a}}$ is a Schwartz function we deduce that

$$T({\boldsymbol{z}})=O\big{(}\,|{\boldsymbol{z}}|_{1}^{-\infty}\,\big{)}\;\;\mbox{as $|{\boldsymbol{z}}|_{1}\to\infty$}.$$

This means that

$$\forall N\in{{\mathbb{N}}},\;\;T({\boldsymbol{z}})=O\big{(}\,|{\boldsymbol{z}}|_{1}^{-N}\,\big{)}\;\;\mbox{as $|{\boldsymbol{z}}|_{1}\to\infty$}.$$

We deduce

$$\big{\|}\,\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}-\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\,\big{\|}_{\mathrm{op}}=\|R_{H}({\boldsymbol{x}},{\boldsymbol{y}})\|_{\mathrm{op}}=O\big{(}\,T({\boldsymbol{x}}-{\boldsymbol{y}})\,\big{)}.$$ (4.9)

Observe next that

$$\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=\left[\begin{array}[]{cc}\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}&0\\ 0&\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}\end{array}\right],$$

and

$$\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=\left[\begin{array}[]{cc}\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}&\operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{]}\\ &\\ \operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}&\operatorname{Var}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}\end{array}\right]$$

$$=\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}+\underbrace{\left[\begin{array}[]{cc}0&\operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}})\,\big{]}\\ &\\ \operatorname{Cov}\big{[}\,\nabla\Phi_{\mathfrak{a}}({\boldsymbol{y}}),\nabla\Phi_{\mathfrak{a}}({\boldsymbol{x}})\,\big{]}&0\end{array}\right]}_{=:R_{\nabla}({\boldsymbol{x}},{\boldsymbol{y}})}.$$

We deduce

$$\big{\|}\,\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}-\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\,\big{\|}_{\mathrm{op}}=\|R_{\nabla}({\boldsymbol{x}},{\boldsymbol{y}})\|_{\mathrm{op}}=O\big{(}\,T({\boldsymbol{x}}-{\boldsymbol{y}})\,\big{)}.$$ (4.10)

We denote by $\widetilde{H}({\boldsymbol{x}},y)^{\flat}$ the Gaussian random matrix obtained from $\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})$ by conditioning on $\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})=0$. Similarly, we denote by $\widehat{H}({\boldsymbol{x}},y)^{\flat}$ the Gaussian random matrix obtained from $\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})$ by conditioning on $\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})=0$. The distributions of $\widetilde{H}({\boldsymbol{x}},y)^{\flat}$ and $\widehat{H}({\boldsymbol{x}},y)^{\flat}$ are determined by the Gaussian regression formula.

Since $\widetilde{H}({\boldsymbol{x}},y)$ and $\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})$ are independent we deduce

$$\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})^{\flat}\,\big{]}=\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}.$$

From the Gaussian regression formula we deduce we deduce

$$\begin{split}\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})^{\flat}\,\big{]}=\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\hskip 85.35826pt&\\ -\operatorname{Cov}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}}),\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}\operatorname{Cov}\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}}),\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}&\end{split}$$

$$\begin{split}\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})^{\flat}\,\big{]}=\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}+R_{H}({\boldsymbol{x}},{\boldsymbol{y}})\hskip 28.45274pt&\\ -\operatorname{Cov}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}}),\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}\operatorname{Cov}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}}),\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}&\end{split}$$

Now observe that

$$\operatorname{Cov}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}}),\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=O\big{(}\,T({\boldsymbol{x}}-{\boldsymbol{y}})\,\big{)}$$

$$\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}=\Big{(}\,\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}+R_{\nabla}({\boldsymbol{x}},{\boldsymbol{y}})\,\Big{)}^{-1}$$

$$\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}\Big{(}\,{\mathbbm{1}}+\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}R_{\nabla}({\boldsymbol{x}},{\boldsymbol{y}})\,\Big{)}^{-1}$$

Since $\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}$ is independent of ${\boldsymbol{x}}$ and ${\boldsymbol{y}}$ we deduce

$$\operatorname{Cov}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}}),\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}\operatorname{Cov}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}}),\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}=O\big{(}\,T({\boldsymbol{x}},{\boldsymbol{y}})\,\big{)},$$

and

$$\Big{|}\,\det\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}-\det\operatorname{Var}\big{[}\,\nabla{\widetilde{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}^{-1}\,\Big{|}=O\big{(}\,T({\boldsymbol{x}}-{\boldsymbol{y}})\,\big{)}.$$ (4.11)

Since $\operatorname{Var}\big{[}\,\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}$ is independent of ${\boldsymbol{x}},{\boldsymbol{y}}$ we deduce that

$$\sup_{{\boldsymbol{x}}\neq{\boldsymbol{y}}}\|\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})^{\flat}\,\big{]}\|_{\mathrm{op}}<\infty,$$

and we deduce from (4.7) that

$$\Big{|}\,{\mathbb{E}}\big{[}\,|\det\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})|\,\big{]}-{\mathbb{E}}\big{[}\,|\det\widetilde{H}({\boldsymbol{x}},{\boldsymbol{y}})|\,\big{]}\,\Big{|}=O\big{(}\,T({\boldsymbol{x}}-{\boldsymbol{y}})^{1/2}\,\big{)}.$$ (4.12)

We can now estimate the right-hand-side of (4.5). Note that if $B_{{\vec{k}}}\cap B_{\vec{\ell}}=\emptyset$ we have

$$\sup_{({\boldsymbol{x}},{\boldsymbol{y}})\in B_{{\vec{k}}}\times B_{\vec{\ell}}}O\big{(}\,T({\vec{k}}-{\vec{\ell}})\,\big{)}=O\big{(}\,{\vec{k}}-{\vec{\ell}}|_{1}^{-\infty}\,\big{)}.$$

Since

$$\sup_{{\boldsymbol{x}}\neq{\boldsymbol{y}}}\big{(}\,\big{\|}\,\operatorname{Var}\big{[}\,\widehat{H}({\boldsymbol{x}},{\boldsymbol{y}})^{\flat}\,\big{]}\,\big{\|}_{\mathrm{op}}+\big{\|}\,\operatorname{Var}\big{[}\,\nabla{\widehat{\Phi}}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{]}\,\big{\|}_{\mathrm{op}}\,\big{)}<\infty$$

we deduce from (4.11) and (4.12) that

$$\sup_{{\boldsymbol{x}}\neq{\boldsymbol{y}}}\big{|}\,\widehat{\rho}({\boldsymbol{x}},{\boldsymbol{y}})-\widetilde{\rho}({\boldsymbol{x}},{\boldsymbol{y}})\,\big{|}=O\big{(}\,\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}^{-\infty}\,\big{)}.$$

Hence

$$C({\vec{k}},{\vec{\ell}})=O\big{(}\,\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}^{-\infty}\,\big{)}\;\;\mbox{as $\big{|}\,{\vec{k}}-{\vec{\ell}}\,\big{|}_{1}\to\infty$}.$$

The above estimate coupled with (4.1) implies that, for any $p>0$ ,there exists a constant $C=C(p)>0$ such that

$$\big{|}\,C({\vec{k}},{\vec{\ell}})\,\big{|}\leq\frac{C(p)}{\big{(}\,1+|{\vec{k}}-{\vec{\ell}}|_{1}\,\big{)}^{p}},\;\;\forall({\vec{k}},{\vec{\ell}})\in{{\mathbb{N}}}^{m}\times{{\mathbb{N}}}^{m}.$$

Theorem 1.1 now follows from the above estimate and Corollary 3.3.