Better Estimate of Four Moments Theorems and Its Application to Polynomials

1. Introduction

The aim of this paper is to find the better estiamte of the four moments theorems of a random variable belonging to Markov chaos studied by Bourguin et al. in the paper [3]. The first study in this field is the central limit theorem, called the fourth moment theorem, in [18] studied by Nualart and Peccati. These authors found a necessary and sufficient condition such that a sequence of random variables, belonging to a fixed Wiener chaos, converges in distribution to a Gaussian random variable. More precisely, let $(X=\{X\left(h\right),h\in \mathfrak{H}\left\}\right)$ be an isonormal Gaussian process defined on a probability space $(\Omega ,\mathfrak{F},\mathbb{P})$, where $\mathfrak{H}$ is a real separable Hilbert space.

Theorem 1. 

[Fourth moment theorem]Fix an integer $q\ge 2$, and let $\{F_n,n\ge 1\}$ be a sequence of random variables belonging to the qth Wiener chaos with $E\left[F_n^2\right]=1$ for all $n\ge 1$. Then $F_n\stackrel{\mathcal{L}}{⟶}Z$ if and only if $E\left[F_n^4\right]\to 3$, where Z is a standard Gaussian random variable and the notation $\stackrel{\mathcal{L}}{⟶}$ denotes the convergence in distribution.

Such a result gives a dramatic simplication of the method of moments from the point of view of convergence in distribution. The above fourth moment theorem is expressed in terms of Malliavin derivative in [17]. However, the results given in [17,18] do not provide any information about the rate of convergenc, whereas, in the paper [10], the authors prove that Theorem 1 can be recovered from the estimate of the Kolmogorov (or total variation, Wasserstein) distance obtained by using the techniques based on the combination between Malliavin calculus (see, e.g., [13,15,16]) and Stein’s method for normal approximation (see, e.g., [4,20,21]). For more explanation of these techniques, we refer to the papers [6,9,10,11,12,13,14].

One of the remarkable achievements of Nourdin-Peccati approach (see Theorem 3.1 in [10]) is the quantification of fourth moment theorem for functionals of Gaussain fields. In the particular case where F is an element in the qth Wiener chaos of X with $\mathbb{E}\left[F^2\right]=1$, the upper bound of Kolmogorov distance is given by

$$d_{Kol}(F,Z)\le \sqrt{\frac{q-1}{3q}}\sqrt{\mathbb{E}\left[F^4\right]-3}.$$

(1)

Here $\mathbb{E}\left[F^4\right]-3$ is just the fourth cumulant ${\kappa }_4\left(F\right)$ of F.

Recently, the author in [8] proves that the fourth moment theorem also holds in the general framework of Markov diffusion generators. More precisely, under a certain spectral condition on Markov diffusion generator, a sequence of eigenfunctions of such a generator satisfies the bound given in (1). In particular, this new method may avoid the use of complicated product formula of multiple integrals. After this work, the authors in [1] introduce a Markov choas of eigenfunctions being less restrictive than Markov chaos defined in [8]. Using this Markov chaos, they derive the quantitative four moments theorem for convergence of the eigenfuctions towards Gaussian, Gamma, Beta distributions. Furthermore, the authors in [3] that the convergence of the elements of a Markov chaos to a Pearson distribution can be still bounded with just the first four moments by using the new concept of chaos grade.

For the purposes of this paper, we will start by referring to the estimate given in Theorem 3.9 obtained by Bourguin et al. in [3]. Pearson diffusions are Itô diffusion given by the following stochastic differential equation(sde)

$$d{X}_t=a\left(X_t\right)dt+\sqrt{2\theta b\left(X_t\right)}d{B}_t,$$

(2)

where $a\left(x\right)=-\theta (x-m)$ and $b\left(x\right)=b_2{x}^2+b_{1}x+b_0$. Given the generator L defined on $L^2(E,\mu )$ by

$$Lf\left(x\right)=-\theta (x-m)f^{\prime }\left(x\right)+\theta b\left(x\right)f^{\prime \prime }\left(x\right),$$

(3)

its invariant measure $\mu$ is a Pearson distribution and the set of eigenvalue of L is given by

$$\Lambda =\left\{-n(1-(n-1)b_2)\theta :n\in {\mathbb{N}}_0,b_2<\frac{1}{2n-1}\right\}.$$

(4)

Theorem 2. 

(Bourguin et al. (2019)) Let ν be a Pearson distribution associated to the diffusion given by sde (2). Let F be a chaotic eigenfunction of generator L with eigenvalue $-\lambda$, chaos grade and moments up to 4. Set $G=F+m$ and $\xi =\mathfrak{u}-2(1-b_2)$ . Then, it holds

$$\begin{array}{ccc}& & \mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-b\left(G\right)\right)}^2\right]\hfill \\ & \le & 2\left(1-b_2-\frac{\mathfrak{u}}{4}\right)\mathbb{E}\left[U\left(G\right)\right]+\frac{{\left(\xi \right)}_{+}(1-b_2)}{2}\mathbb{E}\left[Q^2\left(G\right)\right],\hfill \end{array}$$

(5)

where ${\left(\xi \right)}_{+}=\xi$ for $\xi >0$ and 0 for $\xi \le 0$, and the polynomials Q and U are given by

$$\begin{array}{ccc}\hfill Q\left(x\right)& =& x^2+\frac{2(b_1+m)}{2b_2-1}x+\frac{1}{b_2-1}\left(b_0+\frac{m(b_1+m)}{2b_2-1}\right),\hfill \end{array}$$

(6)

$$\begin{array}{ccc}\hfill U\left(x\right)& =& (1-b_2)Q^2\left(x\right)-\frac{1}{12}{\left(Q^{\prime }\left(x\right)\right)}^3(x-m).\hfill \end{array}$$

(7)

The notations $\Gamma$ and $L^{-1}$ in the above theorem, related to Markov generator, are explained in Section 2.

In this paper, we improve the estimate given in Theorem 2 by introducing the notion of the lower chaos grade in the set of eigenvalues of generator L. For example, if the target distribution $\nu$ in Theorem 2 is a standard Gaussian measure, then the diffusion coefficients are given as $b_2=b_1=0$ and $b_0=1$. Since a chaotic random variable $F=I_q\left(f\right)$, $f\in {\mathfrak{H}}^{\otimes q}$ with $\mathbb{E}\left[F^2\right]=1$, has the chaos grade $\mathfrak{u}=2$, the second term in the bound (5) is vanished and the bound is given as follows:

$$\mathbb{E}\left[{\left(\Gamma (F,-L^{-1}F)-1\right)}^2\right]\le \frac{1}{3}\left(\mathbb{E}\left[F^4\right]-3\right).$$

(8)

Note that $d_{Kol}(F,Z)\le \sqrt{\mathbb{E}\left[{(\Gamma (F,-L^{-1}F)-1)}^2\right]}$. Hence the bound in (1) provides a better estimate for four moments theorem in comparison with bound of (8) in the case of $F=I_q\left(f\right)$. In this paper, we will develop a new technique that provides more improved bounds as above.

Also we give two bounds, called the four moments theorem and fourth moment theorem respectively, for the normal approximation of the case where a random variable F comes from eigenfunctions of a Jacobi generator. One of the bounds is from our main result, Theorem 3 below, and the other bound, obtained using the result in [7]. shows that the fourth moment theorem holds even if the upper chaos grade is greater than two.

The rest of the paper is organized as follows: Section 2 reviews some basic notations and results of Markov diffusion generator. Our main result, in particular the bound in Theorem 3, is presented in Section 3, Finally, as an application of our main results, in Section 4, we consider the case where a random variable G in Theorem 2 comes from an eiegnfunction of a generator associated to a Pearson distribution.

2. Preliminaries

In this section, we recall some basic facts about Markov diffusion generator. The reader is referred to [2] for a more detailed explanation. We begin by the definition of Markov triple $(E,\mathfrak{F},\mu )$ in the sense of [2]. For the infinitesimal generator L of a Markov semigroup $\mathbf{P}={\left(P_t\right)}_{t\ge 0}$ with ${\mathbb{L}}^2\left(\mu \right)$-domain $\mathcal{D}\left(L\right)$, we associated a bilinear form $\Gamma$. Assume that we are given a vector space ${\mathcal{A}}_0$ of $\mathcal{D}\left(L\right)$ such that for every $(F,G)$ of random variables defined on a probability space $(E,\mathfrak{F},\mu )$, the product $FG$ is in $\mathcal{D}\left(L\right)$ (${\mathcal{A}}_0$ is an algebra). On this algebra ${\mathcal{A}}_0$, the bilinear map (carré du champ operator) $\Gamma$ is defined

$$\Gamma (F,G)=\frac{1}{2}(L\left(FG\right)-FLG-GLF).$$

for every $(F,G)\in {\mathcal{A}}_0\times {\mathcal{A}}_0$. As the carré du champ operator $\Gamma$ and the measure $\mu$ completely determine the symmetric Markov generator L, we will work throughout this paper with Markov triple $(E,\mathfrak{F},\mu )$ equipped with a probability measure $\mu$ on a state space $(E,\mathfrak{F})$ and a symmetric bilinear map $\Gamma :{\mathcal{A}}_0\times {\mathcal{A}}_0$ such that $\Gamma (F,F)\ge 0$.

Next, we construct domain $\mathcal{D}\left(\mathcal{E}\right)$ of the Dirichlet form $\mathcal{E}$ by completion of ${\mathcal{A}}_0$, and then obtain, from this Dirchlet domain, domain $\mathcal{D}\left(L\right)$ of L. Recall the Dirchlet form $\mathcal{E}$ as

$$\mathcal{E}(F,G)=\mathbb{E}\left[\Gamma (F,G)\right]for(F,G)\in {\mathcal{A}}_0\times {\mathcal{A}}_0.$$

(9)

If ${\mathcal{A}}_0$ is endowed with the norm

$${\parallel F\parallel }_{\mathcal{E}}={\left[{\parallel F\parallel }_2^2+\mathcal{E}(F,F)\right]}^{1/2},$$

(10)

the completion of ${\mathcal{A}}_0$ with respect to this norm turns it into a Hilbert space embedded in ${\mathbb{L}}^2\left(\mu \right)$. Once the Dirchlet domian $\mathcal{D}\left(\mathcal{E}\right)$ is contructed, the domaion $\mathcal{D}\left(L\right)\subseteq \mathcal{D}\left(\mathcal{E}\right)$ is defined as all elements $F\in \mathcal{D}\left(\mathcal{E}\right)$ such that

$$\left|\mathcal{E}(F,G)\right|\le c_F\mathbb{E}\left[G^2\right]$$

for all $G\in \mathcal{D}\left(\mathcal{E}\right)$, where $c_F$ is a finite constant only depending on F. On these domains, a relation of L and $\Gamma$ holds, namely the integration by parts formula

$$\mathbb{E}\left[\Gamma \right(F,G\left)\right]=-\mathbb{E}\left[FLG\right]=-\mathbb{E}\left[GLF\right].$$

(11)

By the integration by parts Formula (11) and $\Gamma (F,F)\ge 0$, the operator $-L$ is nonnegative and symmetric, and therefore the spectrum of $-L$ is contained $\mathcal{S}\subseteq [0,\infty )$ We assume that $-L$ has discrete spectrum $\mathcal{S}=\{{\lambda }_k,k\ge 0\}$. Obviously, the zero is always an eigenfunction such that $-L\left(1\right)=0$.

A Full Markov triple is a Standard Markov triple for which there is an extended algebra ${\mathcal{A}}_0\subset \mathcal{A}$, with no requirement of integrability for elements of $\mathcal{A}$, satisfying the requirements given in Section 3.4.3 of [2]. In particular, the diffusion property holds: for any ${\mathcal{C}}^{\infty }$ function $\Psi ;{\mathbb{R}}^k\to \mathbb{R}$, and $F_1,\dots ,F_k,G\in \mathcal{A}$,

$$\Gamma (\Psi (F_1,\dots F_k),G)=\sum _{i=1}^k{\partial }_i\Psi (F_1,\dots F_k)\Gamma (F_i,G),$$

(12)

and

$$\begin{array}{ccc}\hfill L(\Psi (F_1,\dots F_k)& =& \sum _{i=1}^k{\partial }_i\Psi (F_1,\dots F_k)L{F}_i\hfill \\ & & +\sum _{i,j=1}^k{\partial }_{ij}\Psi (F_1,\dots F_k)\Gamma (F_i,F_j).\hfill \end{array}$$

(13)

We also define the operator $L^{-1}$, called the pseudo-inverse of L, satisfying for any $F\in \mathcal{D}\left(L\right)$,

$$L{L}^{-1}F=L^{-1}LF=F-\mathbb{E}\left[F\right].$$

(14)

Obviously, this pseudo-inverse $L^{-1}$ is naturally constructed and defined on $\mathcal{D}\left(L\right)$ by a self-adjointness of the operator L.

3. Main Results

We denote the set of eigenvalues of the generator L by $\Lambda \subseteq (-\infty ,0]$. Then chaotic random variables are defined as follows:

Definition 1. 

An eigenfunction F with respect to an eigenvalue λ of the generator $-L$ is calledchaoticif there exists $\mathfrak{u}>1$ and $\mathfrak{e}\le 1$ such that $-\mathfrak{u}\lambda$ and $-\mathfrak{e}\lambda$ are eigenvalues of L, and

$$F^2\in \underset{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda \setminus \left\{0\right\}}{\mathfrak{e}\lambda \le \kappa \le \mathfrak{u}\lambda }}{⨁}\mathrm{Ker}(L+\kappa Id)⨁Ker\left(L\right).$$

(15)

The smallest $\mathfrak{u}$ satisfying (15) is called the upper chaos grade of F, and the largest $\mathfrak{e}$ satisfying (15) is called the lower chaos grade of F.

Now we improve the estimate given in Theorem 2 described in the introduction.

Theorem 3. 

Let ν be a Pearson distribution associated to the diffusion given by sde (2). Let F be a chaotic eigenfunction of generator L with eigenvalue $-\lambda$, chaos grade $\mathfrak{u}$ and moments up to 4. Set $G=F+m$. Then, we have

$$\begin{array}{ccc}& & \mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-b\left(G\right)\right)}^2\right]\hfill \\ & \le & 2\left(1-b_2-\frac{\mathfrak{u}+\mathfrak{e}}{4}\right)\left(\frac{1}{3}-b_2\right)\mathbb{E}\left[\tilde{U}\left(G\right)\right]\hfill \\ & & +\left(\frac{1-b_2}{2}-\frac{\mathfrak{e}}{4}\right)\left(\mathfrak{u}-2(1-b_2)\right)\mathbb{E}\left[Q^2\left(G\right)\right],\hfill \end{array}$$

(16)

where Q is given by (6), and

$$\begin{array}{ccc}\hfill \tilde{U}\left(x\right)& =& x^4+\frac{3(1-b_2)}{1-3b_2}(Q^2\left(x\right)-x^4)\hfill \\ & & -\frac{1}{4(1-3b_2)}\left({\left(Q^{\prime }\left(x\right)\right)}^3)(x-m)-8x^4\right).\hfill \end{array}$$

Proof: From the proof of Theorem 3.9 in [3], we write

$$\Gamma (G,-L^{-1}G)-b\left(G\right)=\frac{1}{2\lambda }\left(L+2(1-b_2)\lambda Id\right)\left(Q\left(G\right)\right),$$

(17)

where $Q\left(x\right)$ is a quadratic polynomial given by (6). By the assumption,

$$Q\left(G\right)=\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}{J}_{\kappa }Q\left(G\right).$$

(18)

Direct computations yield, together with (18), that

$$\begin{array}{ccc}& & \mathbb{E}\left[{(\Gamma (G,-L^{-1}G)-b\left(G\right))}^2\right]\hfill \\ & =& \frac{1}{4{\lambda }^2}\mathbb{E}\left[\left(LQ\left(G\right)+2(1-b_2)\lambda Q\left(G\right)\right)\left(L+2(1-b_2)\lambda \right)Q\left(G\right)\right]\hfill \\ & =& \frac{1}{4{\lambda }^2}\{\mathbb{E}\left[LQ\left(G\right)\left(L+2(1-b_2)\lambda \right)\left(Q\left(G\right)\right)\right]\hfill \\ & & +2(1-b_2)\lambda \mathbb{E}\left[Q\left(G\right)\left(L+2(1-b_2)\lambda \right)Q\left(G\right)\right]\}\hfill \\ & =& \frac{1}{4{\lambda }^2}\{\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(-\kappa )(2(1-b_2)\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \end{array}$$

(19)

$$\begin{array}{ccc}& & +2(1-b_2)\lambda \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(2(1-b_2)\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\}\hfill \\ & =& \frac{1}{4{\lambda }^2}\{\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(-\kappa )(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & +(2(1-b_2)-\mathfrak{u})\lambda \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(-\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & +2(1-b_2)\lambda \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & +2(1-b_2)(2(1-b_2)-\mathfrak{u}){\lambda }^2\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\}\hfill \\ & =& \frac{1}{4{\lambda }^2}\{\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(-\kappa )(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & +(4(1-b_2)-\mathfrak{u})\lambda \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & -(2(1-b_2)-\mathfrak{u})\mathfrak{u}{\lambda }^2\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\}\hfill \end{array}$$

(20)

$$\begin{array}{ccc}& & +2(1-b_2)(2(1-b_2)-\mathfrak{u}){\lambda }^2\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\}.\hfill \end{array}$$

(21)

Since $-\kappa \le -\mathfrak{e}\lambda$ for all $-\kappa \in \Lambda \setminus \left\{0\right\}$, we have that

$$\begin{array}{ccc}& & \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(-\kappa )(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & =& \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda \setminus \left\{0\right\}}{\kappa \le \mathfrak{u}\lambda }}(-\kappa )(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & \le & -\mathfrak{e}\lambda \sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right].\hfill \end{array}$$

(22)

Using (21) yields that

$$\begin{array}{ccc}& & \mathbb{E}\left[{(\Gamma (G,-L^{-1}G)-b\left(G\right))}^2\right]\hfill \\ & \le & \frac{1}{4{\lambda }^2}\{4\lambda \left((1-b_2)-\frac{\mathfrak{u}+\mathfrak{e}}{4}\right)\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}(\mathfrak{u}\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & -(2(1-b_2)-\mathfrak{u})\mathfrak{u}{\lambda }^2\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\hfill \\ & & +2(1-b_2)(2(1-b_2)-\mathfrak{u}){\lambda }^2\sum _{\genfrac{}{}{0pt}{}{-\kappa \in \Lambda }{\kappa \le \mathfrak{u}\lambda }}\mathbb{E}\left[J_{\kappa }{\left(Q\left(G\right)\right)}^2\right]\}\hfill \\ & =& \frac{1}{4{\lambda }^2}\{4\lambda \left((1-b_2)-\frac{\mathfrak{u}+\mathfrak{e}}{4}\right)\mathbb{E}\left[Q\left(G\right)(L+\mathfrak{u}\lambda Id)Q\left(G\right)\right]\hfill \\ & & +(\mathfrak{u}-2(1-b_2))\mathfrak{u}{\lambda }^2\mathbb{E}\left[Q^2\left(G\right)\right]\hfill \\ & & -2(1-b_2)(\mathfrak{u}-2(1-b_2)){\lambda }^2\mathbb{E}\left[Q^2\left(G\right)\right]\}\hfill \\ & =& \frac{1}{4{\lambda }^2}\{4\lambda \left((1-b_2)-\frac{\mathfrak{u}+\mathfrak{e}}{4}\right)\mathbb{E}\left[Q\left(G\right)(L+2(1-b_2)\lambda Id)Q\left(G\right)\right]\hfill \\ & & +4\left(\frac{1-b_2}{2}-\frac{\mathfrak{e}}{4}\right)(\mathfrak{u}-2(1-b_2)){\lambda }^2\mathbb{E}\left[Q^2\left(G\right)\right]\}\hfill \\ & =& 2\left(1-b_2-\frac{\mathfrak{u}+\mathfrak{e}}{4}\right)\left(\frac{1}{3}-b_2\right)\mathbb{E}\left[\tilde{U}\left(G\right)\right]\hfill \\ & & +\left(\frac{1-b_2}{2}-\frac{\mathfrak{e}}{4}\right)(\mathfrak{u}-2(1-b_2))\mathbb{E}\left[Q^2\left(G\right)\right].\hfill \end{array}$$

(23)

Here, for the last equality in (22), we use the following equality obtained from the proof of Theorem 3.9 in [3],

$$\begin{array}{ccc}\hfill \mathbb{E}\left[Q\right(G\left)\right(L+\mathfrak{u}\lambda Id\left)Q\right(G\left)\right]& =& \mathbb{E}\left[Q\left(G\right)(L+2(1-b_2)\lambda Id)Q\left(G\right)\right]\hfill \\ & & +(\mathfrak{u}-2(1-b_2))\lambda \mathbb{E}\left[Q^2\left(G\right)\right]\hfill \\ & =& 2\lambda \mathbb{E}\left[U\left(G\right)\right]+(\mathfrak{u}-2(1-b_2))\lambda \mathbb{E}\left[Q^2\left(G\right)\right].\hfill \end{array}$$

4. Application to Three Polynomials

In this section, three examples will be given in order to illustrate the estimate (16) with the explict expression. For this, we consider the case where a random variable F in Theorem 3 comes from eigenfunctions of a generator associated to a Pearson distribution. For simplicity, we only consider one-dimensional case, analogus results in finite or infinite dimensional case can be extended in a similar way.

4.1. Ornstein-Uhlenbeck Generator

We consider the one dimensional Ornstein-Uhlenbeck generator L, defined for any test function f by

$$Lf\left(x\right)=f^{\prime \prime }\left(x\right)-x{f}^{\prime }\left(x\right)forx\in R,$$

action on $L^2(\mathbb{R},\mu )$, where

$$\mu \left(dx\right)=\frac{1}{\sqrt{2\pi }}{e}^{-\frac{x^2}{2}}dx.$$

Let us set $F=H_q\left(x\right)$ where $H_q$ denotes the Hermite polynomial of order q. Then we have that $F\in Ker(L+qId)$.

Corollary 1. 

Let ν be a Gaussain distribution associated with the diffusion given by (2) with mean m and $b_0={\sigma }^2$. If $F=H_q\left(x\right)$, $q\ge 2$, and $G=F+m$, then we have

$$\mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-{\sigma }^2\right)}^2\right]\le \left(\frac{q-1}{3q}\right)\left(\mathbb{E}\left[F^4\right]-6{\sigma }^2\mathbb{E}\left[F^2\right]+3{\sigma }^4\right),$$

(24)

Proof. 

By the well-known product formula, the square of F can be expressed as a linear combination of Hermite polynomials up to order $2q$ such as

$$H_q^2\left(x\right)=\sum _{r=0}^{q}r!{\left(\genfrac{}{}{0pt}{}{q}{r}\right)}^2{H}_{2(q-r)}\left(x\right).$$

(25)

This product formula (24) gives that the upper chaos grade and lower chaos grade of $H_q$ are $\mathfrak{u}=2$ and $\mathfrak{e}=\frac{2}{q}\le 1$ for $q\ge 2$. Hence Theorem 3 yields that

$$\begin{array}{ccc}\hfill \mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-{\sigma }^2\right)}^2\right]& \le & \frac{2}{3}\left(1-\frac{2+\frac{2}{q}}{4}\right)\mathbb{E}\left[\tilde{U}\left(G\right)\right]\hfill \\ & \le & \left(\frac{q-1}{3q}\right)\mathbb{E}\left[\tilde{U}\left(G\right)\right].\hfill \end{array}$$

(26)

When $b_2=b_1=0$ and $b_0={\sigma }^2$, a directed computation yields that

$$U\left(x\right)=\frac{1}{3}\left[{(x-m)}^4-6{\sigma }^2{(x-m)}^2+3{\sigma }^4\right],$$

so that

$$\mathbb{E}\left[\tilde{U}\left(G\right)\right]=\mathbb{E}\left[F^4\right]-6{\sigma }^2\mathbb{E}\left[F^2\right]+3{\sigma }^4,$$

(27)

From (25) and (26), the proof of the result (23) is completed. □

Remark 1. 

When L is the infinite dimensional Ornstein-Uhlenbeck generator, then $L{I}_q\left(f\right)=-q{I}_q\left(f\right)$, $q=0,1,\dots ,$. Hence the spectrum of L consists of zero and the negative integers with the eigenfunctions being represented by mutiple stochastic integrals. The product formula of the multiple stochastic integrals gives that

$$I_q{\left(f\right)}^2=\sum _{r=0}^{q}r!{\left(\genfrac{}{}{0pt}{}{q}{r}\right)}^2{I}_{2q-2r}\left(f{\otimes }_{r}f\right).$$

This formula shows that theupper chaos gradeandlower chaos gradeof $I_q\left(f\right)$ are still given by $\mathfrak{u}=2$ and $\mathfrak{e}=\frac{2}{q}$ as the one-dimensional case. The upper bound in (1) can be obtained from Theorem 3. □

4.2. Jacobi Generator

We consider the one-dimensional Jacobi generator $L_{\alpha ,\beta }$ defined on $L^2([0,1],{\mu }_{\alpha ,\beta })$ by

$$L_{\alpha ,\beta }f\left(x\right)=(\alpha -(\alpha +\beta )x)f^{\prime }\left(x\right)+x(1-x)f^{\prime \prime }\left(x\right),$$

(28)

where

$${\mu }_{\alpha ,\beta }\left(dx\right)=\frac{\Gamma (\alpha +\beta )}{\Gamma \left(\alpha \right)\Gamma \left(\beta \right)}{x}^{\alpha -1}{(1-x)}^{\beta -1}{\mathbf{1}}_{[0,1]}dx.$$

Its spectrum $\Lambda$ is of the form

$$\Lambda =\left\{-n(n+\alpha +\beta -1):n\in {\mathbb{N}}_0\right\}.$$

(29)

Set ${\lambda }_n=n(n+\alpha +\beta -1)$, $n=0,1,\dots$. Then, we have that

$$L^2([0,1],{\mu }_{\alpha ,\beta })=\underset{n=0}{\overset{\infty }{⨁}}\mathrm{Ker}(L_{\alpha ,\beta }+{\lambda }_{n}Id),$$

and the kernels are given by

$$\mathrm{Ker}(L_{\alpha ,\beta }+{\lambda }_{n}Id)=\left\{a{P}_n^{(\alpha -1,\beta -1)}(1-2x);a\in \mathbb{R}\right\},$$

(30)

where $P_n^{(\alpha ,\beta )}\left(x\right)$ denotes the nth Jacobi polynomials

$$P_n^{(\alpha ,\beta )}\left(x\right)=\frac{{(-1)}^n}{2^{n}n!}{(1-x)}^{-\alpha }{(1+x)}^{-\beta }\frac{d^n}{d{x}^n}\left({(1-x)}^{\alpha +n}{(1+x)}^{\beta +n}\right).$$

(31)

Recall that ${}_p{F}_q$ denotes the generalized hypergeometric function with p numerator and q denominator, given by

$${}_p{F}_q\left(\genfrac{}{}{0pt}{}{\left(a_p\right)}{\left(b_q\right)}|x\right)=\sum _{k=0}^{\infty }\frac{{\left(a_1\right)}_k{\left(a_2\right)}_k\cdots {\left(a_p\right)}_k}{{\left(b_1\right)}_k{\left(b_2\right)}_k\cdots {\left(b_p\right)}_k}\frac{x^k}{k!},$$

where the notation $\left(a_p\right)$ denotes the array of p parameter $a_1,\dots ,a_p$ and

$${\left(\alpha \right)}_n=\frac{\Gamma (\alpha +n)}{\Gamma \left(\alpha \right)}.$$

Then Jacobi polynomials are given by

$$P_n^{(\alpha ,\beta )}\left(x\right)=\frac{{(\alpha +1)}_n}{n!}{}_2{F}_1\left(\genfrac{}{}{0pt}{}{-n,\alpha +\beta +n+1}{\alpha +1}|\frac{1-x}{2}\right).$$

4.2.1. Beta Approximation

In this section, we consider the case when the target distribution $\nu$ is a Beta distribution.

Corollary 2. 

Let ν be the Beta distribution associated to the the diffusion given by (2) with mean ,

$$m=\frac{\alpha }{\alpha +\beta },b_2=-\frac{1}{\alpha +\beta },b_1=\frac{1}{\alpha +\beta }and{b}_0=0.$$

(32)

Let $F=P_n^{(\alpha -1,\beta -1)}(1-2x)$, $n\ge 2$, and set

$$G=F+\frac{\alpha }{\alpha +\beta }for\alpha ,\beta >0.$$

Then we have

$$\begin{array}{ccc}& & \mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-\frac{1}{\alpha +\beta }G(1-G)\right)}^2\right]\hfill \\ & \le & \frac{2}{3}\left(1+\frac{1}{\alpha +\beta }-\frac{2n(2n+\alpha +\beta -1)+\alpha +\beta }{4n(n+\alpha +\beta -1)}\right)\left(\frac{\alpha +\beta +3}{\alpha +\beta }\right)\mathbb{E}\left[\tilde{U}\left(G\right)\right]\hfill \\ & & +\left(\frac{\alpha +\beta +1}{2(\alpha +\beta )}-\frac{\alpha +\beta }{4n(n+\alpha +\beta -1)}\right)\hfill \\ & & \times \left(\frac{2(2n+\alpha +\beta -1)}{(n+\alpha +\beta -1)}-\frac{2(\alpha +\beta +1)}{\alpha +\beta }\right)\mathbb{E}\left[Q^2\left(G\right)\right],\hfill \end{array}$$

(33)

where the constants $b_2$, $b_1$, $b_0$ and m in $\tilde{U}\left(x\right)$ and $Q\left(x\right)$ are given by (31).

Proof. 

The square of a Jacobi polynomial $P_n^{(\alpha ,\beta )}\left(x\right)$ can be expressed as a linear combination of Jacobi polynomials up to order $2n$ as follows:

$$P_n^{(\alpha -1,\beta -1)}{(1-2x)}^2=\sum _{k=0}^{2n}{c}_{n,k}{P}_k^{(\alpha -1,\beta -1)}(1-2x),$$

(34)

where the linearization coefficients $c_{n,k}$ are explicitly given in the paper [5]. This product Formula (33) shows that the upper chaos grade $\mathfrak{u}$ and the lower chaos grade$\mathfrak{e}$ of $P_n^{\alpha -1,\beta -1}$ are given by

$$\begin{array}{ccc}\hfill \mathfrak{u}& =& \frac{{\lambda }_{2n}}{{\lambda }_n}=\frac{2(2n+\alpha +\beta -1)}{(n+\alpha +\beta -1)},\hfill \end{array}$$

(35)

$$\begin{array}{ccc}\hfill \mathfrak{e}& =& \frac{{\lambda }_1}{{\lambda }_n}=\frac{\alpha +\beta }{n(n+\alpha +\beta -1)}.\hfill \end{array}$$

(36)

Hence from (34) and (35) together with $b_2=-\frac{1}{\alpha +\beta }$, the upper bound (32) follows. □

4.2.2. Normal Approximation

In this section, we consider the case when the target distribution $\nu$ is a standard Gaussian measure. Then the diffusion coefficients are given as $b_2=b_1=0$ and $b_0=1$. For simplicity, we will deal with the second Jacobi polynomials $P_2^{\alpha -1,\alpha -1}(1-2x)$ for $\alpha >0$, defined on $L^2([0,1],{\mu }_{\alpha ,\alpha })$, for the case n=2 and $\alpha =\beta$ in (30). Let us set

$$F=\frac{P_2^{\alpha -1,\alpha -1}(1-2x)}{\parallel P_2^{\alpha -1,\alpha -1}{(1-2·)\parallel }_{L^2([0,1],{\mu }_{\alpha ,\alpha })}}.$$

(37)

Then it is obvious that F has $\mathbb{E}\left[F\right]=0$ and $\mathbb{E}\left[F^2\right]=1$. From (34) and (35), it follows that

$$\begin{array}{ccc}\hfill \mathfrak{u}& =& \frac{{\lambda }_4}{{\lambda }_2}=\frac{2(2\alpha +3)}{(2\alpha +1)},\hfill \end{array}$$

(38)

$$\begin{array}{ccc}\hfill \mathfrak{e}& =& \frac{{\lambda }_1}{{\lambda }_2}=\frac{\alpha }{(2\alpha +1)}.\hfill \end{array}$$

(39)

This implies that the upper chaos grade has $\mathfrak{u}>2$ and the lower chaos grade $\mathfrak{e}<1$. By Theorem 3, the bound is given as follows:

$$\begin{array}{ccc}& & \mathbb{E}\left[{\left(\Gamma (F,-L^{-1}F)-1\right)}^2\right]\hfill \\ & \le & \frac{2}{3}\left(1-\frac{\mathfrak{u}+\mathfrak{e}}{4}\right)\left(\mathbb{E}\left[F^4\right]-3\right)+\frac{(2-\mathfrak{e})(\mathfrak{u}-2)}{4}\mathbb{V}ar\left(F^2\right).\hfill \end{array}$$

(40)

Even when the fourth cumulant of F in the first term of (39) is 0, we may not be able to guarantee that F has a standard Gaussian distribution because of the second term in (39). This shows that the fourth moment theorem of Theorem 1 may not hold.

To overcome this problem, a new techique, in [7], has been developed to show that the fourth moment theorem (Theorem 4 below) holds even though the chaos grade is greater than two. Let F be a chaotic eigenfunction of $-L$ with respect to $\lambda$ with $\mathbb{E}\left[F\right]=0$ and $\mathbb{E}\left[F^2\right]=1$. We define a linear function $\varphi \left(x\right)=mx+b$, where

$$\begin{array}{ccc}\hfill m& =& \frac{1}{4\lambda }\sum _{\mathfrak{e}\lambda \le \kappa \le \mathfrak{u}\lambda }(2\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(F^2\right)}^2\right],\hfill \\ \hfill b& =& -\frac{1}{4{\lambda }^2}\sum _{\mathfrak{e}\lambda \le \kappa \le \mathfrak{u}\lambda }\kappa (2\lambda -\kappa )\mathbb{E}\left[J_{\kappa }{\left(F^2\right)}^2\right].\hfill \end{array}$$

Here $J_{\kappa }\left(F^2\right)$ denotes the projection of $F^2$ on $\mathrm{Ker}(L+\kappa Id)$.

Theorem 4. 

If $m\ne 0$, then we have

$$\mathbb{E}\left[{(\Gamma (F,-L^{-1}F)-1)}^2\right]=\frac{2-c_{m,b}}{6}\left(\mathbb{E}\left[F^4\right]-3\right),$$

(41)

where ${\mathfrak{c}}_{m,b}$ is a constsnt such that $\varphi \left({\mathfrak{c}}_{m,b}\right)=0$.

Proof. 

Using the argument in the proof of Theorem 3 in [7] shows that , for any $x\in \mathbb{R}$,

$$\mathbb{E}\left[{(\Gamma (F,-L^{-1}F)-1)}^2\right]=2m-xm+\varphi \left(x\right).$$

(42)

Since $m\ne 0$, there exists a constant $c_{m,b}$, depending on m and b, such that $\varphi \left(c_{m,b}\right)=0$. Also the proof of Theorem 3 in [7] shows that

$$m=\frac{1}{6}\left(\mathbb{E}\left[F^4\right]-3{\left(\mathbb{E}\left[F^2\right]\right)}^2\right).$$

(43)

Plugging $c_{m,b}$ into x in (41) yields, together with (42), that (4) holds. □

We will use Theorem 4 to find that, given F as (36), under which conditions the fourth moment theorem is working by removing the second term in (39). Define a linear function $\varphi \left(x\right)=mx+b$, where the slope m and the intercept b are

$$\begin{array}{ccc}\hfill m& =& 16{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^2\{(\alpha +1){(2\alpha -1)}^2\hfill \end{array}$$

$$\begin{array}{ccc}& & -6\times {\left(24\right)}^2{(\alpha +2)}^5{(\alpha +3)}^7(2\alpha +1)(2\alpha +3)(2\alpha +5)\}.\hfill \\ \hfill b& =& 16{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^2\{12\times {\left(24\right)}^2{(\alpha +2)}^5{(\alpha +3)}^7{(2\alpha +3)}^2(2\alpha +5)\hfill \end{array}$$

(44)

$$\begin{array}{ccc}& & -(\alpha +1){(2\alpha -1)}^2(2\alpha +1)\}.\hfill \end{array}$$

(45)

Theorem 5. 

Let F be a chaotic random variable given by (36). If $\alpha >0$, one has that,

$$\mathbb{E}\left[{(\Gamma (F,-L^{-1}F)-1)}^2\right]=\frac{{\mathfrak{c}}_{\alpha }}{6}\left(\mathbb{E}\left[F^4\right]-3\right).$$

(46)

Here ${\mathfrak{c}}_{\alpha }$ is a positive constant given by

$$\begin{array}{ccc}\hfill c_{\alpha }& =& -\frac{4\vartheta \left(\alpha \right)-{(2\alpha -1)}^3(\alpha +1)}{\vartheta \left(\alpha \right)(2\alpha +1)-{(2\alpha -1)}^2(\alpha +1)},\hfill \end{array}$$

(47)

where

$$\vartheta \left(\alpha \right)=6\times {\left(24\right)}^2{(\alpha +2)}^5{(\alpha +3)}^7(2\alpha +3)(2\alpha +5).$$

Proof. 

When $n=2$ in (33), the linearization coefficients $c_{n,k}$ are given by

$$\begin{array}{ccc}\hfill c_{2,0}& =& 2{(2\alpha +1)}_2{\left[{\left(\alpha \right)}_2\right]}^2,\hfill \end{array}$$

(48)

$$\begin{array}{ccc}\hfill c_{2,2}& =& 8{\left[{(\alpha +1)}_1\right]}^3{\left(\alpha \right)}_3{(2\alpha -1)}_2{(2\alpha +1)}_1,\hfill \end{array}$$

(49)

$$\begin{array}{ccc}\hfill c_{2,4}& =& 24{(\alpha +2)}_2{\left(\alpha \right)}_2{\left(\alpha \right)}_4{\left[{(2\alpha +1)}_2\right]}^2,\hfill \end{array}$$

(50)

and $c_{n,k}=0$ if k is odd, where ${\left(x\right)}_n=x(x+1)\cdots (x+n-1)$. Note that the general form $c_{n,0}$ of (47) is also given by

$$c_{n,0}=n!{(n+2\alpha -1)}_n{\left[{\left(\alpha \right)}_n\right]}^2.$$

(51)

Since $P_0^{(\alpha -1,\alpha -1)}(1-2x)=1$, we have, from (33), that

$${\int }_0^1{P}_2^{(\alpha -1,\alpha -1)}{(1-2x)}^2\nu \left(dx\right)=c_{2,0}.$$

(52)

By orthogonality, we have that

$$\begin{array}{ccc}\hfill {\int }_0^1{P}_4^{(\alpha -1,\alpha -1)}{(1-2x)}^2\nu \left(dx\right)& =& \sum _{k=0}^8{c}_{4,k}{\int }_0^1{P}_k^{(\alpha -1,\alpha -1)}(1-2x)\nu \left(dx\right)\hfill \\ & =& c_{4,0}.\hfill \end{array}$$

(53)

Since $c_{n,k}=0$ for $k=1,3$ and

$$F=\frac{P_2^{\alpha -1,\alpha -1}(1-2x)}{\sqrt{c_{2,0}}},$$

the intercept of a linear function $\varphi$ can be written, using (51) and (52), as

$$\begin{array}{ccc}\hfill 4{\lambda }_2^{2}b& =& -{\lambda }_2(2{\lambda }_2-{\lambda }_2)\mathbb{E}\left[J_2{\left(F^2\right)}^2\right]-{\lambda }_4(2{\lambda }_2-{\lambda }_4)\mathbb{E}\left[J_4{\left(F^2\right)}^2\right]\hfill \\ & =& -{\lambda }_2(2{\lambda }_2-{\lambda }_2){\left(\frac{c_{2,2}}{c_{2,0}}\right)}^2{\int }_0^1{P}_2^{(\alpha -1,\alpha -1)}{(1-2x)}^2\nu \left(dx\right)\hfill \\ & & -{\lambda }_4(2{\lambda }_2-{\lambda }_4){\left(\frac{c_{2,4}}{c_{2,0}}\right)}^2{\int }_0^1{P}_4^{(\alpha -1,\alpha -1)}{(1-2x)}^2\nu \left(dx\right)\hfill \\ & =& -{\lambda }_2(2{\lambda }_2-{\lambda }_2)\frac{c_{2,2}^2}{c_{2,0}}-{\lambda }_4(2{\lambda }_2-{\lambda }_4){\left(\frac{c_{2,4}}{c_{2,0}}\right)}^2{c}_{4,0}.\hfill \end{array}$$

(54)

Using (47)–(50), the right-hand side of (53) can be computed as

$$\begin{array}{ccc}\hfill 4{\lambda }_2^{2}b& =& -64{\lambda }_2^2{\alpha }^2{(\alpha +1)}^5{(\alpha +2)}^2{(2\alpha -1)}^2(2\alpha +1)\hfill \\ & & +128\times {\left(24\right)}^3{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^7{(\alpha +3)}^7{(2\alpha +1)}^2\hfill \\ & & \times {(2\alpha +3)}^2(2\alpha +5).\hfill \end{array}$$

(55)

Hence we have

$$\begin{array}{ccc}\hfill b& =& -16{\alpha }^2{(\alpha +1)}^5{(\alpha +2)}^2{(2\alpha -1)}^2(2\alpha +1)\hfill \\ & & +8\times {\left(24\right)}^3{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^7{(\alpha +3)}^7{(2\alpha +3)}^2(2\alpha +5)\hfill \\ & =& 16{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^2\{12\times {\left(24\right)}^2{(\alpha +2)}^5{(\alpha +3)}^7{(2\alpha +3)}^2(2\alpha +5)\hfill \\ & & -(\alpha +1){(2\alpha -1)}^2(2\alpha +1)\}.\hfill \end{array}$$

(56)

From (56), we see that $b>0$ for $\alpha >0$. Using the same arguments as for the case of b shows that

$$\begin{array}{ccc}\hfill 4{\lambda }_{2}m& =& {\lambda }_2\frac{c_{2,2}^2}{c_{2,0}}+(2{\lambda }_2-{\lambda }_4){\left(\frac{c_{2,4}}{c_{2,0}}\right)}^2{c}_{4,0}\hfill \\ & =& 64{\lambda }_2{\alpha }^2{(\alpha +1)}^5{(\alpha +2)}^2{(2\alpha -1)}^2(2\alpha +1)\hfill \\ & & -32\times {\left(24\right)}^3{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^7{(\alpha +3)}^7{(2\alpha +1)}^2\hfill \\ & & \times (2\alpha +3)(2\alpha +5).\hfill \end{array}$$

(57)

So

$$\begin{array}{ccc}\hfill m& =& 16{\alpha }^2{(\alpha +1)}^5{(\alpha +2)}^2{(2\alpha -1)}^2\hfill \\ & & -4\times {\left(24\right)}^3{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^7{(\alpha +3)}^7(2\alpha +1)\hfill \\ & & \times (2\alpha +3)(2\alpha +5)\hfill \\ & =& 16{\alpha }^2{(\alpha +1)}^4{(\alpha +2)}^2\{(\alpha +1){(2\alpha -1)}^2\hfill \\ & & -6\times {\left(24\right)}^2{(\alpha +2)}^5{(\alpha +3)}^7(2\alpha +1)(2\alpha +3)(2\alpha +5)\}.\hfill \end{array}$$

(58)

Obviously, the right-hand side of (57) shows that $m<0$ for $\alpha >0$. Now, we will find a point x such that $\varphi \left(x\right)=0$. From (55) and (57), the solution of $\varphi \left(x\right)=0$ is given by

$$\begin{array}{ccc}\hfill x& =& -\frac{b}{m}\hfill \\ & =& \frac{2\vartheta \left(\alpha \right)(2\alpha +3)-{(2\alpha -1)}^2(\alpha +1)(2\alpha +1)}{\vartheta \left(\alpha \right)(2\alpha +1)-{(2\alpha -1)}^2(\alpha +1)}.\hfill \end{array}$$

(59)

Hence, it follows from (58) that

$$\begin{array}{ccc}\hfill c_{\alpha }& =& 2+\frac{b}{m}\hfill \\ & =& -\frac{4\vartheta \left(\alpha \right)-{(2\alpha -1)}^3(\alpha +1)}{\vartheta \left(\alpha \right)(2\alpha +1)-{(2\alpha -1)}^2(\alpha +1)}.\hfill \end{array}$$

(60)

□

Remark 2. 

In Theorem 5, we assume that $\alpha =\beta$ (ultraspherical case). This assumption shows that the factor ${}_9{F}_8$ being the generalized hypergeometric function with 9 numerators and 8 denominator parameters, given in the paper [19], vanishes, so that Rahman’s formula is considerably simplified. This assumption allows us to find a point x satisfying $\varphi \left(x\right)\le 0$ quickly and explicitly. □

4.3. Romanovski-Routh Generator

We consider the one-dimensional generator $L_{\alpha ,\beta }$, action on $L^2(\mathbb{R},{\mu }_{p,q})$, where

$${\mu }_{p,q}\left(dx\right)=\frac{\Gamma \left(p+\frac{iq}{2}\right)\Gamma \left(p-\frac{iq}{2}\right)}{2^{2(1-p)}\Gamma (2p-1)}{(x^2+1)}^{-p}exp\left(qarctan\left(x\right)\right)dx,$$

defined by

$$L_{p,q}f\left(x\right)=-\frac{x^2}{2(1-p)}{f}^{\prime \prime }\left(x\right)-\left(x-\frac{q}{2(p-1)}\right)f^{\prime }\left(x\right).$$

(61)

Its spectrum $\Lambda$ is of the form

$$\Lambda =\left\{-n\left(1+\frac{n-1}{2(1-p)}\right):n\in {\mathbb{N}}_0\right\}.$$

(62)

Corollary 3. 

Let ν be the skew t-distribution with mean $\frac{q}{2(p-1)}$ and diffusion coefficients given by

$$b_2=\frac{1}{2(p-1)}and{b}_0=\frac{1}{2(p-1)}.$$

Let $F=R_n^{p,q}\left(x\right)$, $p,q\in \mathbb{R}$ and $n\ge 2$ and set

$$G=F+\frac{q}{2(p-1)}forp,q\in \mathbb{R}.$$

If $p>2n+\frac{1}{2}$, then we have

$$\begin{array}{ccc}& & \mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-\frac{1}{2(p-1)}{G}^2-\frac{1}{2(p-1)}\right)}^2\right]\hfill \\ & \le & 2\left(\frac{2p-3}{2(p-1)}-\frac{4n-4p+3}{4(n-2p+1)}\right)\left(\frac{2p-5}{6(p-1)}\right)\hfill \\ & & \times \mathbb{E}\left[\tilde{U}\left(G\right)\right]\hfill \\ & & +\left(\frac{2p-3}{4(p-1)}-\frac{1}{4(n-2p+1)}\right)\hfill \\ & & \times \left(\frac{2(2n-2p+1)}{(n-2p+1)}-\frac{2p-3}{p-1}\right)\mathbb{E}\left[Q^2\left(G\right)\right],\hfill \end{array}$$

(63)

where the constants $b_2$, $b_1$, $b_0$ and m in $\tilde{U}\left(x\right)$ and $Q\left(x\right)$ are given by

$$b_2=\frac{1}{2(p-1)},b_1=0,b_0=\frac{1}{2(p-1)}andm=\frac{q}{2(p-1)}.$$

Proof. 

First note that the Romanovski-Routh polynomials $R_n^{(p,q)}\left(x\right)$ can be represented by complexified Jacobi polynomials:

$$R_n^{(p,q)}\left(x\right)=n!{\left(2i\right)}^n{P}_n^{(-p+\frac{qi}{2},-p-\frac{qi}{2})}\left(ix\right),$$

(64)

where $P_n^{(p,q)}\left(x\right)$ are the well-known Jacobi polynomials and $i=\sqrt{-1}$. By using (33) and (63), the square of a Ramanovski-Routh polynomial $R_n^{(p.q)}\left(x\right)$ can be expressed as a linear combination of Jacobi polynomials up to order $2n$ as follows:

$$\begin{array}{ccc}\hfill R_n^{(p,q)}{\left(x\right)}^2& =& {(n!)}^2{\left(2i\right)}^{2n}{P}_n^{(-p+\frac{qi}{2},-p-\frac{qi}{2})}{\left(ix\right)}^2\hfill \\ & =& {(n!)}^2{(-4)}^n\sum _{k=0}^{2n}{c}_{n,k}{P}_k^{(-p+\frac{qi}{2},-p-\frac{qi}{2})}\left(ix\right)\hfill \\ & =& {(n!)}^2{(-4)}^n\sum _{k=0}^{2n}\frac{c_{n,k}}{k!{\left(2i\right)}^k}{R}_k^{(p,q)}\left(x\right).\hfill \end{array}$$

(65)

where the linearization coefficients $c_{n,k}$ are explicitly given in the paper [5]. By Proposition 4.2 in [3], the random variable F is chaotic. This product formula (64) shows that the upper chaos grade $\mathfrak{u}$ and the lower chaos grade$\mathfrak{e}$ of $R_n^{(p,q)}$ are

$$\begin{array}{ccc}\hfill \mathfrak{u}& =& \frac{{\lambda }_{2n}}{{\lambda }_n}=\frac{-2n\left(1+\frac{2n-1}{2(1-p)}\right)}{-n\left(1+\frac{n-1}{2(1-p)}\right)}=\frac{2(2n-2p+1)}{n-2p+1},\hfill \end{array}$$

(66)

$$\begin{array}{ccc}\hfill \mathfrak{e}& =& \frac{{\lambda }_1}{{\lambda }_n}=\frac{-1}{-n\left(1+\frac{n-1}{2(1-p)}\right)}=\frac{1}{n-2p+1}.\hfill \end{array}$$

(67)

Hence it follows, from (65) and (66) together with $b_2=\frac{1}{2(p-1)}$, that

$$\begin{array}{ccc}& & \mathbb{E}\left[{\left(\Gamma (G,-L^{-1}G)-\frac{1}{2(p-1)}{G}^2-\frac{1}{2(p-1)}\right)}^2\right]\hfill \\ & \le & 2\left(\frac{2p-3}{2(p-1)}-\frac{4n-4p+3}{4(n-2p+1)}\right)\left(\frac{2p-5}{6(p-1)}\right)\hfill \\ & & \times \mathbb{E}\left[\tilde{U}\left(G\right)\right]\hfill \\ & & +\left(\frac{2p-3}{4(p-1)}-\frac{1}{4(n-2p+1)}\right)\hfill \\ & & \times \left(\frac{2(2n-2p+1)}{(n-2p+1)}-\frac{2p-3}{p-1}\right)\mathbb{E}\left[Q^2\left(G\right)\right],\hfill \end{array}$$

(68)

where the constants $b_2$, $b_1$, $b_0$ and m in $\tilde{U}\left(x\right)$ and $Q\left(x\right)$ are given by

$$b_2=\frac{1}{2(p-1)},b_1=0,b_0=\frac{1}{2(p-1)}andm=\frac{q}{2(p-1)}.$$

□

5. Conclusions and Future Works

The motivation of this study is that the bound in (1) provides a better estimate for four moments theorem in comparison with bound of (8) in the case of $F=I_q\left(f\right)$. We need to develope a new method for obtaining a more improved bound than the bound given in [3]. For this, we find the largest number except zero in the set of eigenvalues corresponding to its eigenfunction in the case where the square of a random variable F, coming from a Markov triple structure, can be expressed as a sum of eigenfunctions,

Future works will be carried out in two directions: (1) We will develop a new technique that can show that the fourth moment theorem like Theorem 4 holds even when the target distribution is not Gaussian. (2) We will study how the second term of the bound (16) in Theorem 3 can be removed even though the chaos grade is greater than two