Green Measures for a Class of Non-markov Processes

1. Introduction

The goal of this paper (see Theorem 1 and Corollary 1 below) is to prove the existence of the Green measure for a class of non-Gaussian processes in ${\mathbb{R}}^d$, called generalized grey Brownian motion (ggBm for short). We denote this family of processes by $B_{\beta ,\alpha }$ with parameters $0<\beta \le 1$ and $0<\alpha \le 2$. More precisely, for a Borel function $f:{\mathbb{R}}^d⟶\mathbb{R}$, the potential of f (see [1,2] for details) is defined as

$$V_{\beta ,\alpha }(f,x)={\int }_0^{\infty }\mathbb{E}\left[f(x+B_{\beta ,\alpha }\left(t\right))\right]\mathrm{d}t,x\in {\mathbb{R}}^d.$$

(1)

We would like to investigate the class of functions f for which the potential of f has the representation

$$V_{\beta ,\alpha }(f,x)={\int }_{{\mathbb{R}}^d}f\left(y\right){\mathcal{G}}_{\beta ,\alpha }(x,\mathrm{d}y),$$

(2)

where $\mathcal{G}(x,·):={\mathcal{G}}_{\beta ,\alpha }(x,·)$ is a Radon measure on ${\mathbb{R}}^d$ called Green measure corresponding to the ggBm $B_{\beta ,\alpha }$, see Definition 2 below. First, we establish the existence of the perpetual integral (cf. Theorem 1)

$${\int }_0^{\infty }f(x+B_{\beta ,\alpha }\left(t\right))\mathrm{d}t$$

with probability 1. This leads to an explicit representation of the Green measure for ggBm, namely (cf. Corollary 1)

$${\mathcal{G}}_{\beta ,\alpha }(x,\mathrm{d}y)=\frac{D}{{|x-y|}^{d-2/\alpha }}\mathrm{d}y,d\alpha >2,1<\alpha \le 2,$$

where D is a constant that depends on $\beta ,\alpha$, and the dimension d; see Theorem 1 for the explicit expression. Note that as $d\alpha >2$ and $1<\alpha \le 2$, the Green measure ${\mathcal{G}}_{\beta ,\alpha }(x,·)$ exists for $d\ge 2$, since $d>2/\alpha \in [1,2)$. The Brownian case ($\alpha =1)$ is covered only for $d\ge 3$.

We emphasize that the existence of the Green measure for a given process X is not always guaranteed. As an example, the d-dimensional Brownian motion (Bm) starting at $x\in {\mathbb{R}}^d$ has a density given by $p_t(x,y)={\left(2\pi t\right)}^{-d/2}exp\left({-|x-y|}^2/\left(2t\right)\right)$, $y\in {\mathbb{R}}^d$. It is not difficult to see that ${\int }_0^{\infty }{p}_t(x,y)\mathrm{d}t$ does not exist for $d=1,2$. This implies the non-existence of the Green measure of Bm for $d=1,2$. On the other hand, for $d\ge 3$, the Green measure of Bm on ${\mathbb{R}}^d$ exists and is given by $\mathcal{G}(x,\mathrm{d}y)={C\left(d\right)|x-y|}^{2-d}\mathrm{d}y$, where $C\left(d\right)$ is a constant depending on the dimension d; see [3] and the references therein for more details. In a two-dimensional space, the Green measure of ggBm is determined by the parameter $\alpha$. The Green measure of ggBm for $d=1$ requires further analysis (for Bm see [4], Ch. 4) which we will postpone for a future paper.

The paper is organized as follows. In Section 2 we recall the definition and main properties of ggBm that will be needed later. In Section 3 we show the existence of the perpetual integral with probability 1, which leads to the explicit formula for the Green measure for ggBm. In Section 4, we discuss the results obtained, connect them with other topics, and draw conclusions.

2. Generalized Grey Brownian Motion

We recall the class of non-Gaussian processes, called the generalized grey Brownian motion, which we study below. This class of processes was first introduced by Schneider [5,6], and was generalized by Mura et al. (see [7,8]) as a stochastic model for slow/fast anomalous diffusion described by the time fractional diffusion equation.

2.1. Definition and Properties

For $0<\beta \le 1$ the Mittag-Leffler (entire) function $E_{\beta }$ is defined by the Taylor series

$$E_{\beta }\left(z\right):=\sum _{n=0}^{\infty }\frac{z^n}{\mathsf{\Gamma }(\beta n+1)},z\in \mathbb{C},$$

(3)

where

$$\mathsf{\Gamma }\left(z\right)={\int }_0^{\infty }{t}^{z-1}{e}^{-t}dt,z\in \mathbb{C},\mathrm{Re}\left(z\right)\ge 0$$

is the Euler gamma function.

The M-Wright function is a special case of the class of Wright functions $W_{\lambda ,\mu }$, $\lambda >-1$, $\mu \in \mathbb{C}$ via

$$M_{\beta }\left(z\right):=W_{-\beta ,1-\beta }(-z)=\sum _{n=0}^{\infty }\frac{{(-z)}^n}{n!\mathsf{\Gamma }(-\beta n+1-\beta )}.$$

The special choice $\beta =1/2$ yields the Gaussian density on $[0,\infty )$

$$M_{1/2}\left(z\right)=\frac{1}{\sqrt{\pi }}\exp \left(-\frac{z^2}{4}\right).$$

(4)

The Mittag-Leffler function $E_{\beta }$ is the Laplace transform of the M-Wright function, that is,

$$E_{\beta }(-s)={\int }_0^{\infty }{e}^{-s\tau }{M}_{\beta }\left(\tau \right)d\tau .$$

(5)

The generalized moments of the density $M_{\beta }$ of order $\delta >-1$ are finite and are given (see [7]) by

$${\int }_0^{\infty }{\tau }^{\delta }{M}_{\beta }\left(\tau \right)d\tau =\frac{\mathsf{\Gamma }(\delta +1)}{\mathsf{\Gamma }(\beta \delta +1)}.$$

(6)

Definition 1.

Let $0<\beta \le 1$ and $0<\alpha \le 2$ be given. A d-dimensional continuous stochastic process $B_{\beta ,\alpha }=\left\{B_{\beta ,\alpha }\left(t\right),t\ge 0\right\}$ starting at $0\in {\mathbb{R}}^d$ and defined on a complete probability space $(\mathsf{\Omega },\mathcal{F},\mathbb{P})$, is a ggBm in ${\mathbb{R}}^d$ (see [8] for $d=1$) if:

• $\mathbb{P}(B_{\beta ,\alpha }\left(0\right)=0)=1$, that is, $B_{\beta ,\alpha }$ starts at zero $\mathbb{P}$-almost surely ($\mathbb{P}$-a.s.).
• Any collection $\left\{B_{\beta ,\alpha }\left(t_1\right),\dots ,B_{\beta ,\alpha }\left(t_n\right)\right\}$ with $0\le t_1<t_2<\cdots <t_n<\infty$ has a characteristic function given, for any $\theta =({\theta }_1,\dots ,{\theta }_n)\in {\left({\mathbb{R}}^d\right)}^n$ with ${\theta }_k=({\theta }_{k,1},\cdots {\theta }_{k,d})$, $k=1,\cdots ,n$, by

  $$\mathbb{E}\left[\exp \left(\mathrm{i}\sum _{k=1}^n{({\theta }_k,B_{\beta ,\alpha }\left(t_k\right))}_{{\mathbb{R}}^d}\right)\right]=E_{\beta }\left[-\frac{1}{2}\sum _{j=1}^d{({\theta }_{.,j},{\gamma }_{\alpha }{\theta }_{.,j})}_{{\mathbb{R}}^n}\right],$$

  (7)
  where $\mathbb{E}$ denotes the expectation w.r.t. $\mathbb{P}$ and

  $${\gamma }_{\alpha }:={\gamma }_{\alpha ,n}:={\left(t_k^{\alpha }+t_j^{\alpha }-{|t_k-t_j|}^{\alpha }\right)}_{k,j=1}^n.$$
• The joint probability density function of $\left(B_{\beta ,\alpha }\left(t_1\right),\dots ,B_{\beta ,\alpha }\left(t_n\right)\right)$ is equal to

  $${\rho }_{\beta }(\theta ,{\gamma }_{\alpha })=\frac{{\left(2\pi \right)}^{-\frac{nd}{2}}}{{\left(\det {\gamma }_{\alpha }\right)}^{d/2}}{\int }_0^{\infty }{\tau }^{-\frac{nd}{2}}{\mathrm{e}}^{-\frac{1}{2\tau }{\sum }_{j=1}^d{({\theta }_{·,j},{\gamma }_{\alpha }^{-1}{\theta }_{·,j})}_{{\mathbb{R}}^n}}{M}_{\beta }\left(\tau \right)d\tau .$$

  (8)

The following are the most important key properties of ggBm:

(P1).

For each $t\ge 0$, the moments of any order of $B_{\beta ,\alpha }\left(t\right)$ are given by

$$\left\{\begin{array}{cc}\mathbb{E}\left[|B_{\beta ,\alpha }{\left(t\right)|}^{2n+1}\right]\hfill & =0,\hfill \\ \mathbb{E}\left[|B_{\beta ,\alpha }{\left(t\right)|}^{2n}\right]\hfill & =\frac{\left(2n\right)!}{2^n\mathsf{\Gamma }(\beta n+1)}{t}^{\alpha n}.\hfill \end{array}\right.$$

(P2).

The covariance function has the form

$$\mathbb{E}\left[\left(B_{\beta ,\alpha }\left(t\right),B_{\beta ,\alpha }\left(s\right)\right)\right]=\frac{d}{2\mathsf{\Gamma }(\beta +1)}\left(t^{\alpha }+s^{\alpha }-{|t-s|}^{\alpha }\right),t,s\ge 0.$$

(9)

(P3).

For each $t,s\ge 0$, the characteristic function of the increments is

$$\mathbb{E}\left[e^{i(k,B_{\beta ,\alpha }\left(t\right)-B_{\beta ,\alpha }\left(s\right))}\right]=E_{\beta }\left(-\frac{{\left|k\right|}^2}{2}{|t-s|}^{\alpha }\right),k\in {\mathbb{R}}^d.$$

(10)

(P4).

The process $B_{\beta ,\alpha }$ is non-Gaussian, $\alpha /2$-self-similar with stationary increments.

(P5).

The ggBm is not a semimartingale. Furthermore, $B_{\alpha ,\beta }$ cannot be of finite variation in $[0,1]$ and, by scaling and stationarity of the increment, on any interval in ${\mathbb{R}}^{+}$.

(P5).

For $n=1$, the density ${\rho }_{\beta }(x,t)$, $x\in {\mathbb{R}}^d$, $t>0$, is the fundamental solution of the following fractional differential equation (see [9])

$${\mathbb{D}}_t^{2\beta }{\rho }_{\beta }(x,t)={\Delta }_x{\rho }_{\beta }(x,t),$$

where ${\Delta }_x$ is the d-dimensional Laplacian in x and ${\mathbb{D}}_t^{2\beta }$ is the Caputo-Dzherbashian fractional derivative; see [10] for the definition and properties.

2.2. Representations of Generalized Grey Brownian Motion

The ggBm admits different representations in terms of well-known processes. It follows from (7) that ggBm has an elliptical distribution, see Section 2 in [11] or Section 3 in [12]. On the other hand, ggBm is also given as a product (see [7] for $d=1$) of two processes as follows

$$\left\{B_{\beta ,\alpha }\left(t\right),t\ge 0\right\}\stackrel{\mathcal{L}}{=}\left\{\sqrt{Y_{\beta }}{B}^{\alpha /2}\left(t\right),t\ge 0\right\}.$$

(11)

Here, $\stackrel{\mathcal{L}}{=}$ means equality in law, the nonnegative random variable $Y_{\beta }$ has density $M_{\beta }$ and $B^{\alpha /2}$ is a d-dimensional fBm with Hurst parameter $\alpha /2$ and independent of $Y_{\beta }$.

We give another representation of ggBm $B_{\beta ,\alpha }$ as a subordination of fBm (see Prop. 2.14 in [13] for $d=1$) which is used below. For completeness, we give the short proof.

Proposition 1.

The ggBm has the following representation

$$\left\{B_{\beta ,\alpha }\left(t\right),t\ge 0\right\}\stackrel{\mathcal{L}}{=}\left\{B^{\alpha /2}\left(t{Y}_{\beta }^{1/\alpha }\right),t\ge 0\right\}.$$

(12)

Proof. 

We must show that both representations (11) and (12) have the same finite-dimensional distribution. For every $\theta =({\theta }_1,\dots ,{\theta }_n)\in {\left({\mathbb{R}}^d\right)}^n$, we have

$$\begin{array}{ccc}\hfill \mathbb{E}\left[exp\left(\mathrm{i}\sum _{k=1}^n\left({\theta }_k,B^{\alpha /2}\left(t_k{Y}_{\beta }^{1/\alpha }\right)\right)\right)\right]& =& {\int }_0^{\infty }\mathbb{E}\left[exp\left(\mathrm{i}\sum _{k=1}^n\left({\theta }_k,B^{\alpha /2}\left(t_k{y}^{1/\alpha }\right)\right)\right)\right]M_{\beta }\left(y\right)\mathrm{d}y\hfill \\ & =& {\int }_0^{\infty }\mathbb{E}\left[exp\left(\mathrm{i}\sum _{k=1}^n\left({\theta }_k,y^{1/2}{B}^{\alpha /2}\left(t_k\right)\right)\right)\right]M_{\beta }\left(y\right)\mathrm{d}y\hfill \\ & =& \mathbb{E}\left[exp\left(\mathrm{i}\sum _{k=1}^n\left({\theta }_k,Y_{\beta }^{1/2}{B}^{\alpha /2}\left(t_k\right)\right)\right)\right].\hfill \end{array}$$

In the second equality, we used the $\alpha /2$-self-similarity of fBm. This completes the proof. □

3. The Green Measure for Generalized Grey Brownian Motion

In this section we show the existence of the Green measure for the ggBm, see (1) and (2). Let us begin by discussing the existence of the Green measure for a general stochastic process X.

Let $X=\left\{X\right(t),t\ge 0\}$ be a stochastic process in ${\mathbb{R}}^d$ starting from $x\in {\mathbb{R}}^d$. If $X\left(t\right)$, $t\ge 0$, has a probability distribution ${\rho }_{X\left(t\right)}(x,·)$, then Eq. (1) becomes

$$V_X(x,f)={\int }_0^{\infty }{\int }_{{\mathbb{R}}^d}f\left(y\right){\rho }_{X\left(t\right)}(x,\mathrm{d}y)\mathrm{d}t.$$

(13)

Then, applying the Fubini theorem, the Green measure ${\mathcal{G}}_X(x,·)$ of X is given by

$${\mathcal{G}}_X(x,\mathrm{d}y)={\int }_0^{\infty }{\rho }_{X\left(t\right)}(x,\mathrm{d}y)\mathrm{d}t,$$

assuming the existence of ${\mathcal{G}}_X(x,·)$ as a Radon measure on ${\mathbb{R}}^d$. That is, for every bounded Borel set $B\in {\mathcal{B}}_b\left({\mathbb{R}}^d\right)$ we have

$${\mathcal{G}}_X(x,B)={\int }_0^{\infty }{\rho }_{X\left(t\right)}(x,B)\mathrm{d}t<\infty .$$

If the probability distribution ${\rho }_{X\left(t\right)}(x,·)$ is also absolutely continuous with respect to the Lebesgue measure, say ${\rho }_{X\left(t\right)}(x,\mathrm{d}y)={\rho }_t(x,y)\mathrm{d}y$, then the function

$$g_X(x,y):={\int }_0^{\infty }{\rho }_t(x,y)\mathrm{d}t,\forall y\in {\mathbb{R}}^d,$$

(14)

is called the Green function of the stochastic process X. Moreover, the Green measure in this case is given by ${\mathcal{G}}_X(x,\mathrm{d}y)=g_X(x,y)\mathrm{d}y$.

This leads us to the following definition of the Green measure of a stochastic process $X.$

Definition 2.

Let $X=\left\{X\right(t),t\ge 0\}$ be a stochastic process on ${\mathbb{R}}^d$ starting from $x\in {\mathbb{R}}^d$ and ${\rho }_{X\left(t\right)}(x,·)$ be the probability distribution of $X\left(t\right)$, $t\ge 0$. The Green measure of X is defined as a Radon measure on ${\mathbb{R}}^d$ by

$${\mathcal{G}}_X(x,B):={\int }_0^{\infty }{\rho }_{X\left(t\right)}(x,B)\mathrm{d}t,B\in {\mathcal{B}}_b\left({\mathbb{R}}^d\right),$$

or

$${\int }_{{\mathbb{R}}^d}f\left(y\right){\mathcal{G}}_X(x,\mathrm{d}y)={\int }_{{\mathbb{R}}^d}f\left(y\right){\int }_0^{\infty }{\rho }_{X\left(t\right)}(x,dy)\mathrm{d}t,f\in C_0\left({\mathbb{R}}^d\right)$$

whenever these integrals exist.

In other words, ${\mathcal{G}}_X(x,B)$ is the expected length of time the process remains in B.

In order to state the main theorem which establishes the existence of the Green measure for ggBm, first, we introduce a proper Banach space of functions $f:{\mathbb{R}}^d⟶\mathbb{R}$ such that

$${\int }_0^{\infty }f\left(x+B_{\beta ,\alpha }\left(t\right)\right)\mathrm{d}t$$

is finite $\mathbb{P}$-a.s. Without loss of generality, we may assume that $f\ge 0$ above. We define the space $CL\left({\mathbb{R}}^d\right)$ of continuous real-valued on ${\mathbb{R}}^d$ by

$$CL\left({\mathbb{R}}^d\right):=\left\{f:{\mathbb{R}}^d\to \mathbb{R}\mid f\mathrm{is}\mathrm{continuous},\mathrm{bounded}\mathrm{and}f\in L^1\left({\mathbb{R}}^d\right)\right\}.$$

The space $CL\left({\mathbb{R}}^d\right)$ becomes a Banach space with the norm

$${\parallel f\parallel }_{CL}:={\parallel f\parallel }_{\infty }+{\parallel f\parallel }_1,\forall f\in CL\left({\mathbb{R}}^d\right),$$

where ${\parallel ·\parallel }_{\infty }$ denotes the sup-norm and ${\parallel ·\parallel }_1$ is the norm in $L^1({\mathbb{R}}^d$). The choice of $CL\left({\mathbb{R}}^d\right)$ allows us to show that the family of random variables (also known as perpetual integral functionals)

$${\int }_0^{\infty }f\left(x+B_{\beta ,\alpha }\left(t\right)\right)\mathrm{d}t,f\in CL\left({\mathbb{R}}^d\right)$$

have finite expectations $\mathbb{P}$-a.s.

Theorem 1.

Let $f\in CL\left({\mathbb{R}}^d\right)$ and $x\in {\mathbb{R}}^d$ be given and consider ggBm $B_{\beta ,\alpha }$ with $d\alpha >2$ and $1<\alpha \le 2$. Then, the perpetual integral functional ${\int }_0^{\infty }f(x+B\left(t\right))\mathrm{d}t$ is finite $\mathbb{P}$-a.s. and its expectation equals

$$\begin{array}{c}\mathbb{E}\left[{\int }_0^{\infty }f(x+B_{\beta ,\alpha }\left(t\right))\mathrm{d}t\right]=D{\int }_{{\mathbb{R}}^d}\frac{f(x+y)}{{\left|y\right|}^{d-2/\alpha }}\mathrm{d}y,\end{array}$$

(15)

where $D=D(\beta ,\alpha ,d)=\frac{1}{\alpha }{2}^{-1/\alpha }{\pi }^{-\frac{d}{2}}\mathsf{\Gamma }\left(\frac{d}{2}-\frac{1}{\alpha }\right)\frac{\mathsf{\Gamma }(1-\frac{1}{\alpha })}{\mathsf{\Gamma }(1-\frac{\beta }{\alpha })}$.

Proof. 

Given $x\in {\mathbb{R}}^d$ and $f\in CL\left({\mathbb{R}}^d\right)$ non-negative, let ${\rho }_{\beta }(·,t^{\alpha })$ denote the density of $B_{\beta ,\alpha }\left(t\right)$, $t\ge 0$, given by (see (8) with $n=1$)

$${\rho }_{\beta }(y,t^{\alpha })=\frac{1}{{\left(2\pi t^{\alpha }\right)}^{d/2}}{\int }_0^{\infty }{\tau }^{-d/2}{\mathrm{e}}^{-\frac{{\left|y\right|}^2}{2t^{\alpha }\tau }}{M}_{\beta }\left(\tau \right)\mathrm{d}\tau ,y\in {\mathbb{R}}^d.$$

First, we show the equality (15). It follows from the above considerations that

$$\begin{array}{c}\mathbb{E}\left[{\int }_0^{\infty }f(x+B_{\beta ,\alpha }\left(t\right))\mathrm{d}t\right]={\int }_0^{\infty }{\int }_{{\mathbb{R}}^d}f(x+y){\rho }_t^{\beta ,\alpha }\left(y\right)\mathrm{d}y\mathrm{d}t.\hfill \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}={\int }_0^{\infty }{\int }_{{\mathbb{R}}^d}f(x+y)\frac{1}{{\left(2\pi t^{\alpha }\right)}^{d/2}}{\int }_0^{\infty }{\tau }^{-d/2}{M}_{\beta }\left(\tau \right)e^{-\frac{{\left|y\right|}^2}{2t^{\alpha }\tau }}\mathrm{d}\tau \mathrm{d}y\mathrm{d}t.\hfill \end{array}$$

Using Fubini’s Theorem, we first compute the t-integral and use the assumption $d\alpha >2$. We obtain

$${\int }_0^{\infty }\frac{1}{{\left(2\pi t^{\alpha }\tau \right)}^{d/2}}{e}^{-\frac{{\left|y\right|}^2}{2t^{\alpha }\tau }}\mathrm{d}t=C(\alpha ,d)\frac{{\tau }^{-\frac{1}{\alpha }}}{{\left|y\right|}^{d-2/\alpha }},$$

where

$$C(\alpha ,d):=\frac{1}{\alpha }{2}^{-1/\alpha }{\pi }^{-\frac{d}{2}}\mathsf{\Gamma }\left(\frac{d}{2}-\frac{1}{\alpha }\right).$$

Next we compute the $\tau$-integral using (6) so that

$${\int }_0^{\infty }{\tau }^{-1/\alpha }{M}_{\beta }\left(\tau \right)\mathrm{d}\tau =\frac{\mathsf{\Gamma }(1-\frac{1}{\alpha })}{\mathsf{\Gamma }(1-\frac{\beta }{\alpha })},\alpha >1.$$

Combining gives

$$\mathbb{E}\left[{\int }_0^{\infty }f(x+B_{\beta ,\alpha }\left(t\right))\mathrm{d}t\right]=D{\int }_{{\mathbb{R}}^d}\frac{f(x+y)}{{\left|y\right|}^{d-2/\alpha }}\mathrm{d}y,$$

where $D=D(\beta ,\alpha ,d)=C(\alpha ,d)\frac{\mathsf{\Gamma }(1-\frac{1}{\alpha })}{\mathsf{\Gamma }(1-\frac{\beta }{\alpha })}$. Therefore, the equality (15) is shown.

Now we show that the right-hand side of (15) is finite for every non-negative $f\in CL\left({\mathbb{R}}^d\right)$. To see this, we may use the local integrability of ${\left|y\right|}^{d-2/\alpha }$ in y and obtain

$$\begin{array}{ccc}\hfill {\int }_{{\mathbb{R}}^d}\frac{f(x+y)}{{\left|y\right|}^{d-2/\alpha }}\mathrm{d}y& =& {\int }_{\left\{\right|y|\le 1\}}\frac{f(x+y)}{{\left|y\right|}^{d-2/\alpha }}\mathrm{d}y+{\int }_{\left\{\right|y|>1\}}\frac{f(x+y)}{{\left|y\right|}^{d-2/\alpha }}\mathrm{d}y\hfill \\ & \le & C_1{\parallel f\parallel }_{\infty }+C_2{\parallel f\parallel }_1\le C{\parallel f\parallel }_{CL}.\hfill \end{array}$$

Therefore, the integral in (15) is, in fact, well defined. In other words, the integral ${\int }_0^{\infty }f\left(x+B_{\beta ,\alpha }\left(t\right)\right)\mathrm{d}t$ exists with probability 1. This completes the proof. □

As a consequence of the above theorem, we immediately obtain the Green measure of ggBm $B_{\beta ,\alpha }$, that is, comparing (2) and (15).

Corollary 1.

The Green measure of ggBm $B_{\beta ,\alpha }$ for $d\alpha >2$ is given by

$${\mathcal{G}}_{\beta ,\alpha }(x,\mathrm{d}y)=\frac{D}{{|x-y|}^{d-2/\alpha }}\mathrm{d}y.$$

Remark 1.

• It is possible to show that given $f\ne 0$, the perpetual integral ${\int }_0^{\infty }f(x+B_{\beta ,\alpha }\left(t\right))\mathrm{d}t$ is a non-constant random variable. As a consequence, for $f\ge 0$ the variance of ${\int }_0^{\infty }f(x+B_{\beta ,\alpha }\left(t\right))\mathrm{d}t$ is strictly positive. The proof uses the notion of conditional full support of ggBm. We will not provide a detailed explanation of this result that closely follows the ideas of Theorem 2.2 in [3] to which we address the interested readers.
• Note also that the functional in (1)

  $$V_{\beta ,\alpha }(·,x):CL\left({\mathbb{R}}^d\right)⟶\mathbb{R}$$
  is continuous. In fact, from the proof of Theorem 1 for any $f\in CL\left({\mathbb{R}}^d\right)$ yields

  $$|V_{\beta ,\alpha }{(f,x)|\le K\parallel f\parallel }_{CL},$$
  where K is a constant depending on the parameters $\beta ,\alpha$, and d.

4. Discussion and Conclusions

We have derived the Green measure for the class of stochastic processes called the generalized grey Brownian motion in Euclidean space ${\mathbb{R}}^d$ for $d\ge 2$. This class includes, in particular, fractional Brownian motion and other non-Gaussian processes. To address the case where $d=1$, a renormalization process is needed. However, this will be postponed to future work. For $\beta =\alpha =1$ ggBm $B_{1,1}$ is nothing but a Brownian motion. In this case, the Green measure exists for $d\ge 3$.

The relationship between the Green measure and the local time of the ggBm can be described as follows. For any $T>0$ and a continuous function $f:{\mathbb{R}}^d⟶\mathbb{R}$, the integral functional

$${\int }_0^{T}f\left(B_{\beta ,\alpha }\left(t\right)\right)\mathrm{d}t$$

(16)

is well-defined. For $d=1$ the integral (16) with $f\in L^1\left(\mathbb{R}\right)$ is represented as

$${\int }_0^{T}f\left(B_{\beta ,\alpha }\left(t\right)\right)\mathrm{d}t={\int }_{\mathbb{R}}f\left(x\right)L_{\beta ,\alpha }(T,x)\mathrm{d}x,$$

where $L_{\beta ,\alpha }(T,x)$ is the local time of ggBm up to time T at the point x, see [11,12]. The Green measure corresponds to the asymptotic behaviour in T of the expectation of local time $L_{\beta ,\alpha }(T,x)$. The existence of this asymptotic depends on the dimension d and the transient or recurrent properties of the process.