**1. Introduction**

In this paper we study in detail the solution of the time-fractional equation

$$\frac{\partial^{\nu}u}{\partial t^{\nu}} = c^2 \sum\_{j=1}^{d} \frac{\partial^2 u}{\partial x\_j^2} \tag{1}$$

for 1 < *ν* < 2 under the initial conditions

$$\begin{cases} u(\mathbf{x},0) = \delta(\mathbf{x}) \\ u\_t(\mathbf{x},0) = 0. \end{cases} \tag{2}$$

The time-fractional derivative is hereafter understood in the Caputo sense:

$$\frac{\partial^{\nu}u}{\partial t^{\nu}} = \frac{1}{\Gamma(m-\nu)} \int\_{0}^{t} \frac{\partial^{m}}{\partial t^{m}} u(\mathbf{x}, \mathbf{s}) \left(t-s\right)^{m-\nu-1} \, \text{ds} \qquad m-1 < \nu < m. \tag{3}$$

We first prove that the Fourier transform of the solution of the Cauchy problem (1) and (2) is

$$\mathcal{U}(\gamma\_1, \dots, \gamma\_d, t) = \mathcal{U}(\gamma, t) = E\_{\nu, 1}(-c^2 ||\gamma||^2 t^{\nu}) \tag{4}$$

where

$$E\_{\nu,1}(\mathbf{x}) = \sum\_{k=0}^{\infty} \frac{\mathbf{x}^k}{\Gamma(\nu k + 1)} \qquad \mathbf{x} \in \mathbb{R} \tag{5}$$

is the one-parameter Mittag-Leffler function, first introduced in [1]. The representation of (4) as a contour integral on the Hankel path *Ha*

$$E\_{\nu,1}(-c^2||\gamma||^2t^{\nu}) = \frac{1}{2\pi i} \int\_{H\_0} \frac{e^w w^{\nu-1}}{w^{\nu} + t^{\nu}c^2 ||\gamma||^2} \,\mathrm{d}w \tag{6}$$

permits us to obtain a representation of (6) as

$$\begin{split} E\_{\nu,1}(-\varepsilon^{2}||\gamma||^{2}t^{\nu}) &= \frac{\sin\pi\nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1}\varepsilon^{-tz\varepsilon^{2/\nu}} \|\gamma\|^{2/\nu}}{z^{2\nu} + 2z^{\nu}\cos\pi\nu + 1} dz \\ &+ \frac{1}{\nu} \left[ \varepsilon^{\varepsilon^{2/\nu} \|\gamma\|^{2/\nu} t e^{i\pi/\nu}} + \varepsilon^{\varepsilon^{2/\nu} \|\gamma\|^{2/\nu} t e^{-i\pi/\nu}} \right]. \end{split} \tag{7}$$

Some details about the representation (6) and the Hankel path can be found in [2]. For *d* = 1, 1 < *ν* < 2 the inversion of (7) is presented in [3] with the conclusion that the solution of (1) is the distribution of a stable symmetric process of order 2/*<sup>ν</sup>*.

We here show that for *d* > 1, 1 < *ν* < 2 the solution can be expressed in terms of the law of a *d*− dimensional stable process *<sup>S</sup>α*(*t*) with a suitable choice of the measure Γ appearing in

$$\mathbb{E}e^{i\gamma \cdot S\_{\mathbb{R}}(t)} = e^{-t\int\_{\mathbb{R}^{d-1}} \|\gamma \cdot s\|^{a} \left(1 - i \operatorname{sign}(\gamma \cdot s) \tan\frac{\pi a}{\mathbb{Z}}\right) \Gamma(ds)}.\tag{8}$$

In particular, for Γ uniform on the upper and lower hemispheres of S*d*−<sup>1</sup> = {*s* ∈ R*<sup>d</sup>* : *s* = <sup>1</sup>}, we prove that (8) yields the characteristic functions in square brackets of Formula (7). We give also the explicit forms of *<sup>u</sup>*(*<sup>x</sup>*, *t*) of the solution of (1) in terms of Bessel functions *J d*2 −<sup>1</sup>(*ρx*), which for *d* = 1 can be reduced to Fujita's result. Some results concerning wave equations of fractional type can be found, e.g., in [4].

#### **2. The Fractional Wave Equation**

In this note we present some relationship between stable processes (and their inverses) with fractional equations. Stable processes are studied in depth in the monograph [5]. Some simple and well known results state that a symmetric stable process *<sup>S</sup>α*(*t*), 0 < *α* ≤ 2 with characteristic function

$$\mathbb{E}\mathcal{E}^{i\gamma\mathbb{S}\_{\mathfrak{a}}(t)} = \mathfrak{e}^{-|\gamma|^{\mathfrak{a}}t} \qquad \gamma \in \mathbb{R}, t > 0 \tag{9}$$

has distribution *pα*(*<sup>x</sup>*, *t*), *x* ∈ R, *t* > 0, satisfying the fractional equation

$$\frac{\partial p}{\partial t} = \frac{\partial^{\kappa}}{\partial |\mathbf{x}|^{\alpha}} p \tag{10}$$

where *∂α ∂*|*x*|*<sup>α</sup>* is the Riesz fractional derivative usually defined as

$$\frac{\partial^a}{\partial |\mathbf{x}|^a} f(\mathbf{x}) = \frac{1}{2\Gamma(m-a)\cos(\pi a/2)} \frac{d^m}{d\mathbf{x}^m} \int\_{-\infty}^{+\infty} \frac{f(y)}{|\mathbf{x}-y|^{a+1-m}} dy \qquad m-1 < a < m \tag{11}$$

with Fourier transform

$$\int\_{-\infty}^{+\infty} e^{i\gamma x} \frac{\partial^a}{\partial |x|^a} f(x) \mathrm{d}x = -|\gamma|^a \int\_{-\infty}^{+\infty} e^{i\gamma x} f(x) \mathrm{d}x. \tag{12}$$

For the *d*-dimensional isotropic stable process *<sup>S</sup>α*(*t*)=(*S*1*α*(*t*), ... , *<sup>S</sup>dα*(*t*)) with characteristic function,

$$\mathbb{E}\mathcal{e}^{i\gamma S\_{\mathfrak{a}}^d(t)} = \mathfrak{e}^{-||\gamma||^a t} \qquad \gamma \in \mathbb{R}^d, t > 0 \tag{13}$$

The corresponding probability law *pα*(*<sup>x</sup>*1,..., *xd*, *t*) = *pα*(*<sup>x</sup>*, *t*) satisfies the equation

$$\frac{\partial p}{\partial t} = -(-\Delta)^a p \tag{14}$$

where −(−<sup>Δ</sup>)*<sup>α</sup>* is the fractional Laplacian defined as the operator such that

$$-\left(-\Delta\right)^{a}f(\mathbf{x}) = \frac{1}{(2\pi)^{d}}\int\_{\mathbb{R}^{d}} e^{-i\left(\gamma,\mathbf{x}\right)} \|\gamma\|^{a}\widehat{f}(\gamma)\mathbf{d}\gamma \; , \qquad f \in \text{Dom}\left((-\Delta)^{a}\right) \tag{15}$$

where ˆ *f*(*γ*) is the Fourier transform of a function *f*(*x*), *x* ∈ R*<sup>d</sup>* and the domain of the operator *Dom*((−<sup>Δ</sup>)*α*) is

$$Dom((-\Delta)^{a}) = \left\{ f \in L^{1}\_{loc}(\mathbb{R}^{n}) : \int\_{\mathbb{R}^{n}} |\widehat{f}(\gamma)|^{2} (1 + \|\gamma\|^{2a}) \, \mathrm{d}\gamma < \infty \right\}$$

(on this point see for example [6]). The connection between fractional operators and stochastic processes is explored, e.g., in [7]. A detailed comparison of the several possible definitions of the fractional Laplacian can be found in [8]. For the time-fractional equation (see [9]),

$$\begin{cases} \frac{\partial^p p}{\partial \mathbf{i}^p} = c^2 \frac{\partial^2 p}{\partial \mathbf{x}^2} & 0 < \boldsymbol{\nu} \le \mathbf{2}, \, \mathbf{x} \in \mathbb{R}, \, t > 0 \\ u(\mathbf{x}, 0) = \boldsymbol{\delta}(\mathbf{x}) \\ u\_t(\mathbf{x}, 0) = 0, \end{cases} \tag{16}$$

we have that the solution of the Cauchy problem is explicitly given by

$$u(\mathbf{x},t) = \frac{1}{2ct^{\nu/2}} \mathcal{W}\_{-\nu/2, 1-\nu/2} \left(-\frac{|\mathbf{x}|}{ct^{\nu}}\right) \tag{17}$$

where

$$\mathcal{W}\_{\mathfrak{a},\mathfrak{\beta}}(\mathbf{x}) = \sum\_{k=0}^{\infty} \frac{\mathfrak{x}^k}{\Gamma(\alpha k + \beta)} \qquad \mathfrak{a} > -1, b > 0, \mathfrak{x} \in \mathbb{R}$$

is the Wright function. The *d*−dimensional counterpart of (16) is

$$\begin{cases} \frac{\partial^{\nu} p}{\partial t^{\nu}} = c^2 \Delta u & 0 < \nu \le 2, \, \mathbf{x} \in \mathbb{R}^d, \, t > 0 \\ u(\mathbf{x}, 0) = \delta(\mathbf{x}) \\ u\_t(\mathbf{x}, 0) = 0. \end{cases} \tag{18}$$

Some details about time-fractional derivatives can be found in [10]. For 0 < *ν* < 1 the solution of (18) corresponds to the distribution of the vector process

$$B\left(L\_{\mathcal{V}}(t)\right), \quad t \ge 0 \tag{19}$$

where *<sup>B</sup>*(*t*)=(*<sup>B</sup>*1(*t*), ... , *Bd*(*t*)) is the *d*−dimensional Brownian motion and *<sup>L</sup>ν*(*t*) is the inverse of the stable subordinator *<sup>H</sup>ν*(*t*) (see [11]).

In the more general case

$$\begin{cases} \frac{\partial^r p}{\partial t^r} = -\varepsilon^2 (-\Delta)^a u & 0 < \nu \le 2, \ 0 < \mu \le 1, \ x \in \mathbb{R}^d, \ t > 0 \\ u(\mathbf{x}, 0) = \delta(\mathbf{x}) \\ u\_t(\mathbf{x}, 0) = 0. \end{cases} \tag{20}$$

The solution of the Cauchy problem (20) is the probability density of the process

$$\mathcal{S}\_{\mathfrak{a}}(L\_{\nu}(t)), \ t \ge 0 \tag{21}$$

where *<sup>S</sup><sup>α</sup>*(*t*)=(*S*1*α*(*t*), ... , *<sup>S</sup>dα*(*t*)) is an isotropic stable process (see [11]). We here consider the case where in (18) and (20) the order of the fractional derivative is 1 < *ν* ≤ 2.We start first with (18) and observe that the Laplace–Fourier transform of the solution *<sup>u</sup>α*(*<sup>x</sup>*, *t*) is

$$\int\_0^\infty e^{-\lambda t} \int\_{\mathbb{R}^d} e^{i\gamma \cdot x} u(x, t) dx = \frac{\lambda^{\nu - 1}}{\lambda^\nu + c^2 ||\gamma||^2} \tag{22}$$

and the Fourier transform reads

$$\int\_{\mathbb{R}^d} e^{i\gamma \cdot x} u(x, t) \, \mathrm{d}x = E\_{\nu, 1}(-c^2 \|\gamma\|^2 t^{\nu}). \tag{23}$$

The Mittag-Leffler function *<sup>E</sup><sup>ν</sup>*,<sup>1</sup>(−*c*<sup>2</sup>*γ*<sup>2</sup>*t<sup>ν</sup>*) can be represented as a contour integral on the Hankel path as

$$E\_{\nu,1}(-c^2||\gamma||^2t^{\nu}) = \frac{1}{2\pi i} \int\_{H\_a} \frac{e^w w^{\nu-1}}{w^{\nu} + t^{\nu}c^2||\gamma||^2} \,\mathrm{d}w \tag{24}$$

where *Ha* is the contour in the complex plane represented in Figure 1.

**Figure 1.** Hankel path in the complex plane.

The representation (24) is a consequence of the integral representation of the inverse of the Gamma function

$$\frac{1}{\Gamma(\nu)} = \frac{1}{2\pi i} \int\_{Ha} e^w w^{-\nu} \mathrm{d}w.$$

The integral in (24) can be developed by inserting a ring of radius < *R*.

The contour *C* is composed by the circumferences *CR* and *C* with two segments joining *Re*−*i<sup>π</sup>* with *e*<sup>−</sup>*i<sup>π</sup>* and *Reiπ* with *ei<sup>π</sup>*, and is run counterclockwise. See Figure 2.

In order to evaluate

$$\frac{1}{2\pi i} \int\_{\mathcal{C}} \frac{c^{\nu} w^{\nu - 1}}{w^{\nu} + t^{\nu} c^{2} ||\gamma||^{2}} \,\mathrm{d}w \tag{25}$$

we perform the transformation *w<sup>ν</sup>* = *z*2*<sup>m</sup>* for 2*m* − 1 < *ν* < 2*m*.

The contour of Figure 2 after the transformation *w<sup>ν</sup>* = *z*2*<sup>m</sup>* takes the form shown in Figure 3.

**Figure 2.** Representation of the contour *C*.

**Figure 3.** Representation of the contour *C*, with *R* = *R ν*2*m* and = *ν*2*m* .

Therefore, the horizontal segments of Figure 2 are rotated by an angle of amplitude ±*πν*/2*m* and the radii are subject to contraction or dilation according to the value of *ν*. The integral on *C* thus obtained from (25) is

$$\frac{1}{2\pi i} \int\_{\mathcal{C}} \frac{\varepsilon^{\text{av}} w^{\nu - 1}}{w^{\nu} + t^{\nu} c^2 ||\gamma||^2} \, \text{d}w = \frac{2m}{2\pi \text{v} i} \int\_{\mathcal{C}'} \frac{\varepsilon^{2m/\nu} z^{2m - 1}}{z^{2m} + t^{\nu} c^2 ||\gamma||^2} \, \text{d}z. \tag{26}$$

The integral on the right side of (26) can be evaluated by means of the Cauchy residue theorem. The function

$$f(z) = \frac{e^{z^{2m/\nu}} z^{2m-1}}{z^{2m} + t^{\nu} c^2 ||\gamma||^2} \qquad z \in \mathbb{C} \tag{27}$$

has 2*m* poles at points *zk* = *ei<sup>π</sup>* (2*k*+<sup>1</sup>) 2*m* (*c*<sup>2</sup>*γ*<sup>2</sup>*t<sup>ν</sup>*) 12*m* for 0 ≤ *k* ≤ 2*m* − 1. It is easy to show that the residues of (27) at the poles *zk* are given by

$$\lim\_{z \to z\_k} (z - z\_k) \frac{e^{z^{2m/\nu}} z^{2m - 1}}{z^{2m} + t^{\nu} c^2 ||\gamma||^2} = \frac{e^{z\_k^{2m/\nu}}}{2m}. \tag{28}$$

Thus the integral (26) can be written as

$$\frac{1}{2\pi i} \int\_{C'} f(z) \,\mathrm{d}z = \frac{1}{2m} \sum\_{k=0}^{2m-1} \varepsilon^{z\_k^{2m/\nu}} \qquad \text{where} \quad z\_k = \varepsilon^{i\pi \frac{2k+1}{2m}} (c^2 ||\gamma||^2 t^\nu)^{\frac{1}{2\pi}} \tag{29}$$

By adding the contribution of the segments (*Re*−*i<sup>π</sup>*, *e*<sup>−</sup>*iπ*) and (*Rei<sup>π</sup>*, *eiπ*) for *R* → ∞ and → 0 we obtain

$$E\_{\nu,1}(-c^2 \|\gamma\|^2 t^{\nu}) = \frac{\sin \pi \nu}{\pi} \int\_0^\infty \frac{z^{\nu-1} e^{-tzx^{2/\nu}} \|\gamma\|^{2/\nu}}{z^{2\nu} + 2z^\nu \cos \pi \nu + 1} dz + \frac{2m}{\nu} \left[ \frac{1}{2m} \sum\_{k=0}^{2m-1} e^{z\_k^{2m/\nu}} \right]. \tag{30}$$

For *m* = 1 we must distinguish the cases 0 < *ν* < 1 where

$$E\_{\nu,1}(-c^2 \|\gamma\|^2 t^{\nu}) = \frac{\sin \pi \nu}{\pi} \int\_0^\infty \frac{z^{\nu-1} e^{-tz^{2/\nu}} \|\gamma\|^{2/\nu}}{z^{2\nu} + 2z^\nu \cos \pi \nu + 1} dz \tag{31}$$

and 1 < *ν* < 2, where

$$\begin{split} E\_{\nu,1}(-\varepsilon^{2}||\gamma||^{2}t^{\nu}) &= \frac{\sin\pi\nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1}e^{-tz^{2/\nu}} \|\gamma\|^{2/\nu}}{z^{2\nu} + 2z^{\nu}\cos\pi\nu + 1} dz \\ &+ \frac{1}{\nu} \left[ \varepsilon^{2/\nu} \|\gamma\|^{2/\nu} t^{\sin\gamma} + \varepsilon^{2/\nu} \|\gamma\|^{2/\nu} t e^{-i\pi/\nu} \right]. \end{split} \tag{32}$$

In order to simplify the formulas involved in the analysis we take *c* = 1.

For *m* = 2 we have the subcases 2 < *ν* < 3 and 3 < *ν* < 4. In the first case the contour integral of Figure 3 involves two poles and thus yields two additional terms in the representation of the Mittag-Leffler function (32). In the second case we have the contribution of four poles in the contour integral of Figure 3, so that for 3 < *ν* < 4

$$\begin{split} E\_{\nu,1}(-\varepsilon^{2}||\gamma||^{2}t^{\nu}) &= \frac{\sin\pi\nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1}e^{-tz^{2/\nu}} \|\gamma\|^{2/\nu}}{z^{2\nu} + 2z^{\nu}\cos\pi\nu + 1} dz \\ &+ \frac{1}{\nu} \left[ \epsilon^{\|\gamma\|^{2/\gamma}t\varepsilon^{\text{ir}/4\nu}} + \epsilon^{\|\gamma\|^{2/\gamma}t\varepsilon^{3\text{i}\pi/4\nu}} + \epsilon^{\|\gamma\|^{2/\gamma}t\varepsilon^{\text{i}\pi/4\nu}} + \epsilon^{\|\gamma\|^{2/\gamma}t\varepsilon^{\text{i}\pi/4\nu}} \right]. \end{split} \tag{33}$$

The contours for 2 < *ν* < 3 and 3 < *ν* < 4 are depicted below (Figure 4).

> < *ν* <

**Figure 4.** Representation of the contour *C* for *m* = 2. The dots indicate the poles of *f*(*z*).

The substantial difference between the cases 1 < *ν* < 2 and *ν* > 2 is that in the first case we have that %(*e*±*iπ*/*ν*) is negative and the contribution of the poles correspond to the characteristic function of stable processes, whereas for 2 < *ν* < 3, %(*eiπ*/*ν*) and %(*e*7*iπ*/*ν*) are positive and thus are not characteristic functions of random variables. Let us now concentrate our attention on the integrals in Equations (31) and (32) (which is also true in the general case for *ν* > 2). If we write

$$\begin{split} &\frac{\sin\pi\nu}{\pi}\int\_{0}^{\infty}\frac{z^{\nu-1}e^{-tz}|\gamma|^{2/\nu}}{z^{2\nu}+2z^{\nu}\cos\pi\nu+1}\mathrm{d}z \\ &=\frac{\sin\pi\nu}{\pi\nu}\int\_{0}^{\infty}\frac{e^{-tw^{1/\nu}}|\gamma|^{2/\nu}}{w^{2}+2w\cos\pi\nu+1}\mathrm{d}w \\ &=\frac{\sin\pi\nu}{\pi\nu}\int\_{0}^{\infty}\frac{e^{-tw^{1/\nu}}|\gamma|^{2/\nu}}{(w+\cos\pi\nu)^{2}+\sin^{2}\pi\nu}\mathrm{d}w \\ &=\mathbb{E}e^{-t||\gamma||^{2/\nu}W^{1/\nu}} \end{split} \tag{34}$$

 < *ν* <

where *W* is a non-negative r.v. with density

$$f(w) = \frac{\sin \pi \nu}{\pi \nu} \frac{\mathrm{d}w}{(w + \cos \pi \nu)^2 + \sin^2 \pi \nu} \qquad w > 0, 0 < \nu < 1$$

Note that for 1 < *ν* < 2 the function (2) is negative on (0, <sup>∞</sup>). We note also that the r.v. *W<sup>ν</sup>*, 0 < *ν* < 1 with density

$$P(\mathcal{W}\_{\boldsymbol{V}} \in \mathrm{d}w)/\mathrm{d}w = \frac{\sin \pi \boldsymbol{v}}{\pi} \frac{w^{\nu - 1}}{1 + w^{2\nu} + 2w^{\nu} \cos \pi \boldsymbol{v}} \qquad w > 0 \tag{35}$$

appearing in (35) has the same distribution as the ratio of two independent stable subordinators of degree 0 < *ν* < 1.

We now give the inverse Fourier transform of (34) for 0 < *ν* < 1, *d* = 1.

$$p\_{\nu}(\mathbf{x},t) = \frac{1}{2\pi} \int\_{-\infty}^{\infty} e^{-i\gamma x} E\_{\nu,1}(-\gamma^2 t^{\nu}) \mathrm{d}\gamma \tag{36}$$

$$= \frac{1}{2\pi} \int\_{-\infty}^{\infty} e^{-i\gamma x} \frac{\sin \pi \nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1} e^{-tz\gamma^{2/\nu}}}{z^{2\nu} + 2z^{\nu} \cos \pi \nu + 1} \mathrm{d}z \, \mathrm{d}\gamma = (t\boldsymbol{z}\gamma^{\frac{2}{\nu}} = w)$$

$$= \frac{1}{2\pi} \int\_{-\infty}^{\infty} e^{-i\gamma x} \frac{\sin \pi \nu}{\pi} \int\_{0}^{\infty} \frac{e^{-w} w^{\nu-1} t^{\nu} \gamma^{2}}{w^{2\nu} + t^{2\nu} \gamma^{4} + 2t^{\nu} \gamma^{2} w^{\nu} \cos \pi \nu} \mathrm{d}w \, \mathrm{d}\gamma.$$

We start by evaluating the following integral

$$\begin{split} &\frac{1}{2\pi} \int\_{-\infty}^{\infty} \varepsilon^{-i\gamma x} \frac{t^{\nu} \gamma^{2}}{w^{2\upsilon} + t^{2\upsilon} \gamma^{4} + 2t^{\upsilon} \gamma^{2} w^{\upsilon} \cos \pi \upsilon} \mathrm{d}\gamma^{\upsilon} \\ &= \frac{1}{2\pi t^{\upsilon/2}} \int\_{-\infty}^{\infty} \frac{\varepsilon^{-i x \gamma' t^{-\upsilon/2}} \gamma^{\upsilon}^{2}}{\gamma^{\upsilon 4} + 2 \gamma^{\upsilon 2} w^{\upsilon} \cos \pi \upsilon + w^{2\upsilon}} \mathrm{d}\gamma^{\upsilon} \\ &= \frac{1}{2\pi (wt)^{\upsilon/2}} \int\_{-\infty}^{\infty} \frac{\varepsilon^{-i x \gamma (w/t)^{\upsilon/2}} \gamma^{2}}{1 + \gamma^{4} + 2 \gamma^{2} \cos \pi \upsilon} \mathrm{d}\gamma \\ &= \frac{1}{2\pi (wt)^{\upsilon/2}} \int\_{-\infty}^{\infty} \frac{\varepsilon^{-i \gamma A} \gamma^{2}}{1 + \gamma^{4} + 2 \gamma^{2} \cos \pi \upsilon} \mathrm{d}\gamma \end{split}$$

where *A* = (*w*/*t*)*ν*/2*x*.

> We must now evaluate the integral of

$$f(z) = \frac{z^2 e^{-izA}}{z^4 + 2z^2 \cos \pi \nu + 1} \qquad z \in \mathbb{C}$$

on a suitable contour *CR*. The four roots of *z*4 + 2*z*<sup>2</sup> cos *πν* + 1 = 0 are

$$\begin{cases} z\_1 = e^{\frac{i\pi\nu}{2} - i\frac{\pi}{2}} = \sin\frac{\pi\nu}{2} - i\cos\frac{\pi\nu}{2} \\ z\_2 = e^{i\frac{\pi\nu}{2} + i\frac{\pi}{2}} = -\sin\frac{\pi\nu}{2} + i\cos\frac{\pi\nu}{2} \\ z\_3 = e^{-i\frac{\pi\nu}{2} - i\frac{\pi}{2}} = -\sin\frac{\pi\nu}{2} - i\cos\frac{\pi\nu}{2} \\ z\_4 = e^{-i\frac{\pi\nu}{2} + i\frac{\pi}{2}} = \sin\frac{\pi\nu}{2} + i\cos\frac{\pi\nu}{2} \end{cases}$$

and are located in *C* as in Figure 5, because 1 < *ν* < 2.

**Figure 5.** Integration contour for *x* > 0.

We observe that *e*<sup>−</sup>*izA* = *e*<sup>−</sup>*<sup>i</sup>*(*u*+*iv*)*x*(*w*/*t*)*ν*/2 = *evAe*−*iuA* for *x* > 0 and *v* < 0, the curvilinear integral on the half-circle *Reiθ*, 0 ≤ *θ* ≤ *π* tends to zero as *R* → ∞. By the residue theorem we thus have

$$\int\_{\mathcal{C}\_{\rm R}} f(z) \, \mathrm{d}z = -2\pi i (R\_{z\_2} + R\_{z\_4}) \tag{37}$$

The minus sign is due to the fact that the contour in Figure 5 is run clockwise. The residues *Rz*2 and *Rz*4 have the following values

$$R\_{z\_2} = -\frac{e^{i\frac{\pi\nu}{2}}e^{A e^{i\pi\nu/2}}}{4\sin\pi\nu} \qquad R\_{z\_4} = \frac{e^{-i\frac{\pi\nu}{2}}e^{A e^{-i\pi\nu/2}}}{4\sin\pi\nu} \tag{38}$$

and thus

$$\begin{split} \int\_{\mathbb{C}\_{\mathbb{R}}} f(z) \, \mathrm{d}z &= -\frac{2\pi i}{2^{2} \sin \pi \upsilon} \left( e^{-i\frac{\pi \upsilon}{2}} e^{A e^{-i\upsilon \upsilon/2}} - e^{i\frac{\pi \upsilon}{2}} e^{A e^{i\upsilon \upsilon/2}} \right) \\ &= -\frac{\pi}{\sin \pi \upsilon} e^{A \cos \frac{\pi \upsilon}{2}} \sin \left( \frac{\pi \upsilon}{2} + A \sin \frac{\pi \upsilon}{2} \right) \end{split} \tag{39}$$

For *x* < 0, the integration of *f*(*z*) must be performed on the contour of Figure 6 and

$$\int\_{\mathbb{C}\_{\mathbb{R}}} f(z) \,\mathrm{d}z = -2\pi i (\mathcal{R}\_{z\_1} + \mathcal{R}\_{z\_3})\,\tag{40}$$

the sign being in this case positive because the path is run counterclockwise.

**Figure 6.** Integration contour for *x* < 0.

The residues in this case are

$$R\_{z\_1} = \frac{e^{i\pi\nu/2}e^{-Ae^{i\pi\nu/2}}}{4\sin\pi\nu} \qquad R\_{z\_3} = -\frac{e^{-i\pi\nu/2}e^{-Ae^{-i\pi\nu/2}}}{4\sin\pi\nu}$$

The integral (40), therefore, takes the form

$$\int\_{\mathbb{C}\_{\mathbb{R}}} f(z) \, \mathrm{d}z = \frac{2\pi i}{4 \sin \pi \nu} \left( e^{i\pi \nu/2} e^{-A e^{i\pi \nu/2}} - e^{-i\pi \nu/2} e^{-A e^{-i\pi \nu/2}} \right) \tag{41}$$

$$= -\frac{\pi}{\sin \pi \nu} e^{-A \cos \frac{\pi \nu}{2}} \sin \left( \frac{\pi \nu}{2} - A \sin \frac{\pi \nu}{2} \right)$$

In conclusion we have that

$$\int\_{-\infty}^{\infty} \frac{e^{-i\mathbf{x}\gamma(w/t)^{v/2}}\gamma^2}{1+\gamma^4+2\gamma^2\cos\pi\nu} \mathrm{d}\gamma \tag{42}$$

$$\mathbf{1} = -\frac{\pi}{\sin\pi\nu} e^{|\mathbf{x}|\cos\frac{\pi\nu}{2}(w/t)^{v/2}} \sin\left(\frac{\pi\nu}{2}+|\mathbf{x}|\left(\frac{w}{t}\right)^{\frac{v}{2}}\sin\frac{\pi\nu}{2}\right).$$

We now consider the integration with respect to *w* in (36). This leads to the evaluation of the following integral

$$\begin{split} &\frac{1}{2\pi} \int\_0^\infty e^{-w} w^{\nu-1} \frac{\sin \pi \nu}{\pi (wt)^{\nu/2}} \left[ -\frac{\pi}{\sin \pi \nu} e^{|x| \cos \frac{\pi \nu}{2} (w/t)^{\nu/2}} \sin \left( \frac{\pi \nu}{2} + |x| \left( \frac{w}{t} \right)^{\frac{\nu}{2}} \sin \frac{\pi \nu}{2} \right) \right] dw &= \\ &= -\frac{1}{2\pi t^{\frac{\nu}{2}}} \int\_0^\infty e^{-w} w^{\frac{\nu}{2} - 1} e^{|x| \cos \frac{\pi \nu}{2} (|\frac{\pi}{2}|^{\nu/2} \sin \left( \frac{\pi \nu}{2} + |x| \left( \frac{w}{t} \right)^{\frac{\nu}{2}} \sin \frac{\pi \nu}{2} \right)} dw \\ &= -\frac{1}{\pi \nu t^{\frac{\nu}{2}}} \int\_0^\infty e^{-tw^{\frac{\nu}{2}}} (w^{\frac{\nu}{2}} t)^{\frac{\nu}{2} - 1} \frac{2}{\nu} t w^{\frac{\nu}{2} - 1} e^{w |x| \cos \frac{\pi \nu}{2}} \sin \left( \frac{\pi \nu}{2} + w|x| \sin \frac{\pi \nu}{2} \right) dw \\ &= -\frac{1}{\pi \nu} \int\_0^\infty e^{-tw^{\frac{\nu}{2}}} e^{w |x| \cos \frac{\pi \nu}{2}} \sin \left( \frac{\pi \nu}{2} + w|x| \sin \frac{\pi \nu}{2} \right) dw \end{split}$$

The last step of (43) can be developed as follows

$$\begin{split} & -\frac{1}{\pi\upsilon} \int\_{0}^{\infty} e^{-tw^{2/\upsilon}} e^{w|\mathbf{x}|\cos\frac{\pi\upsilon}{2}} \left[ e^{i\left(\frac{\pi\upsilon}{2} + w|\mathbf{x}|\sin\frac{\pi\upsilon}{2}\right)} - e^{-i\left(\frac{\pi\upsilon}{2} + w|\mathbf{x}|\sin\frac{\pi\upsilon}{2}\right)} \right] \frac{\mathrm{d}w}{2i} \\ & = -\frac{e^{i\pi\upsilon/2}}{2i\pi\upsilon} \int\_{0}^{\infty} e^{-tw^{2/\upsilon}} e^{w|\mathbf{x}|e^{i\pi\upsilon/2}} \mathrm{d}w + \frac{e^{-i\pi\upsilon/2}}{2i\pi\upsilon} \int\_{0}^{\infty} e^{-tw^{2/\upsilon}} e^{w|\mathbf{x}|e^{-i\pi\upsilon/2}} \mathrm{d}w \end{split} \tag{44}$$

We evaluate the first integral in (44) by taking the contour integral of the function

$$f(z) = e^{-tz^{2/\nu} + |x|ze^{in\nu^2}} \qquad z \in \mathbb{C}$$

along the path *CR* depicted in Figure 7.

**Figure 7.** Path *CR* corresponding to the change of variables *z* = *ei π*2 (<sup>1</sup>−*<sup>ν</sup>*)*z* in the first integral of (44).

By the Cauchy theorem we have that

$$\int\_{C\_{\mathbb{R}}} f(z) \, \mathrm{d}z = \int\_{0}^{R} f(w) \, \mathrm{d}w + \int\_{C} f(R e^{i\theta}) \, \mathrm{d}\theta + \int\_{R}^{0} f(z e^{i\frac{\pi}{2} - i\frac{\pi v}{2}}) e^{i\frac{\pi}{2} - i\frac{\pi v}{2}} \, \mathrm{d}z = 0\tag{45}$$

The integral on the arc *C* tends to zero because

$$\left| \int\_{\mathbb{C}} f(R\varepsilon^{i\theta}) \, \mathrm{d}\theta \right| \leq \int\_{\mathbb{C}} |f(R\varepsilon^{i\theta})| \, \mathrm{d}\theta = \int\_{\frac{(1-\nu)\pi}{2}}^{0} e^{-tR\frac{\hat{\theta}}{\varepsilon} \cos\frac{2\theta}{\nu}} \, e^{|\mathbf{x}|R\cos\left(\theta + \frac{\pi\nu}{2}\right)} \, \mathrm{d}\theta \tag{46}$$

Since *π* > *θ* + *πν*2 > *π*2 as *θ* ∈ ((<sup>1</sup>−*<sup>ν</sup>*)*<sup>π</sup>* 2 , <sup>0</sup>), the exponent of the second factor of (46) is negative as well as the first one because cos(2*θ*/*ν*) for *θ* ranging in the same interval.

We thus conclude that (46) converges to zero as *R* → ∞ and thus (45) yields

$$\int\_0^\infty e^{-tw^{2/\nu}} e^{w|\mathbf{x}|e^{i\frac{\pi\nu}{2}}} \,\mathrm{d}w = \int\_0^\infty e^{i\frac{\pi}{2} - i\frac{\pi\nu}{2}} e^{-t\left(ze^{i\frac{\pi}{2} - i\frac{\pi\nu}{2}}\right)^{2/\nu} + i|\mathbf{x}|z} \,\mathrm{d}z \tag{47}$$
 
$$= i e^{-i\frac{\pi\nu}{2}} \int\_0^\infty e^{i|\mathbf{x}|z} e^{tz^{2/\nu}e^{i\pi/\nu}} \,\mathrm{d}z$$

In order to evaluate the second integral of (44) we integrate

$$f(z) = e^{-tz^{2/\nu} + |x|z e^{-in\nu 2}} \qquad z \in \mathbb{C}$$

along the contour of Figure 8.

**Figure 8.** Path *CR* corresponding to the change of variables *z* = *ei π*2 (*<sup>ν</sup>*−<sup>1</sup>)*z* in the second integral of (44).

By performing the same steps as above we obtain

$$\int\_0^\infty e^{-tw^{2/\nu}} e^{w|\mathbf{x}|e^{-i\frac{\pi\mathbf{w}}{2}}} \,\mathrm{d}w = \int\_0^\infty e^{i\frac{\pi\nu}{2} - i\frac{\pi}{2}} e^{-t\left(ze^{i\frac{\pi\nu}{2} - i\frac{\pi}{2}}\right)^{2/\nu} + |\mathbf{x}|z e^{i\frac{\pi\nu}{2} - i\frac{\pi}{2}} e^{-i\frac{\pi\nu}{2}}} \,\mathrm{d}z \tag{48}$$
 
$$= -ie^{i\frac{\pi\nu}{2}} \int\_0^\infty e^{-i|\mathbf{x}|z} e^{tz^{2/\nu}e^{-i\pi/\nu}} \,\mathrm{d}z$$

in view of (47) and (48) the integral (44) becomes

$$\begin{split} & -\frac{1}{2\pi\nu} \int\_0^\infty e^{i|x|z} e^{tz^{2/\nu}e^{i\pi/\nu}} \, \mathrm{d}z - \frac{1}{2\pi\nu} \int\_0^\infty e^{-i|x|z} e^{tz^{2/\nu}e^{-i\pi/\nu}} \, \mathrm{d}z \\ & = -\frac{1}{2\pi\nu} \int\_0^\infty e^{i|x|z} e^{tz^{2/\nu}e^{i\pi/\nu}} \, \mathrm{d}z - \frac{1}{2\pi\nu} \int\_{-\infty}^0 e^{i|x|z} e^{t(-z)^{2/\nu}e^{-i\pi/\nu}} \, \mathrm{d}z \\ & = -\frac{1}{2\pi\nu} \int\_{-\infty}^{+\infty} e^{i|x|z} e^{t|z|^{2/\nu}e^{i\tilde{\mathbb{P}}\mathrm{sign}z}} \, \mathrm{d}z \\ & = -\frac{1}{2\pi\nu} \int\_{-\infty}^{+\infty} e^{-i|x|z} e^{t|z|^{\tilde{\mathbb{P}}}e^{-i\tilde{\mathbb{P}}\mathrm{sign}z}} \, \mathrm{d}z \end{split}$$

From (32), *d* = 1, we conclude that

$$\begin{split} p\_{\boldsymbol{r}}(\boldsymbol{x},t) &= \frac{1}{2\pi} \int\_{-\infty}^{+\infty} \boldsymbol{\varepsilon}^{-i\gamma\boldsymbol{x}} E\_{\boldsymbol{r},1} (-\gamma^{2} \boldsymbol{t}^{\boldsymbol{r}}) \, d\boldsymbol{\gamma} \\ &= \frac{1}{2\pi\nu} \left[ -\int\_{-\infty}^{+\infty} \boldsymbol{\varepsilon}^{-|\boldsymbol{x}|\gamma} \boldsymbol{\varepsilon}^{\boldsymbol{t}|\gamma|^{2\boldsymbol{\beta}\boldsymbol{\gamma}} - \overleftrightarrow{\boldsymbol{\nabla}}^{\boldsymbol{\theta}\boldsymbol{q}\boldsymbol{\eta}}} \, d\boldsymbol{\gamma} + \int\_{-\infty}^{+\infty} \boldsymbol{\varepsilon}^{-i\boldsymbol{x}\gamma} \boldsymbol{\varepsilon}^{\boldsymbol{t}|\gamma|^{2\boldsymbol{\beta}\boldsymbol{\gamma}} \boldsymbol{t}^{\boldsymbol{\beta}\boldsymbol{q}\boldsymbol{\eta}\boldsymbol{\gamma}}} \, d\boldsymbol{\gamma} + \int\_{-\infty}^{+\infty} \boldsymbol{\varepsilon}^{-i\boldsymbol{x}\gamma} \boldsymbol{J}^{\boldsymbol{t}|\gamma|^{2\boldsymbol{\beta}\boldsymbol{\gamma}} \boldsymbol{t}^{-\boldsymbol{\beta}\boldsymbol{\gamma}}} \, d\boldsymbol{\gamma} \right] \\ &= \frac{1}{2\pi\nu} \int\_{-\infty}^{+\infty} \boldsymbol{\varepsilon}^{\boldsymbol{t}|\boldsymbol{x}\gamma} \boldsymbol{\varepsilon}^{\boldsymbol{t}|\gamma|^{2\boldsymbol{\beta}\boldsymbol{\gamma}} \boldsymbol{t}^{-\boldsymbol{\beta}\boldsymbol{\gamma}} \boldsymbol{t}^{-\boldsymbol{\beta}\boldsymbol{\gamma}} \, d\boldsymbol{\gamma} \end{split} \tag{49}$$

**Remark 1.** *The function h*(*<sup>γ</sup>*, 2*ν* ) = *e<sup>t</sup>*|*γ*|2/*<sup>ν</sup>e*<sup>−</sup> *iπν signγ* = *e<sup>t</sup>*|*γ*| 2*ν cos πν* (<sup>1</sup>−*isign<sup>γ</sup>* tan *πν* ) *is the characteristic function of a stable random variable of order* 1 < 2/*ν* < 2 *with symmetry parameter β* = 1*. Many details about the properties of such densities can be found in [12]. The function p* 2 *ν* (*<sup>x</sup>*, *t*) = 12*π* +∞ −∞ *<sup>e</sup>*<sup>−</sup>*iγxh*(*<sup>γ</sup>*, 2/*ν*)d*γ is unimodal with a positive maximal point and is such that* ∞0 *p* 2*ν* (*<sup>x</sup>*, *t*) = d*x* = *<sup>ν</sup>*/2*. Analogously, the function h*(*<sup>γ</sup>*, −2/*ν*) *is the characteristic function of a negatively skewed random variable. This implies that the function pν can be seen as the superposition of the densities of stable random variables with index β* = ±1*, conditional to be respectively positive or negative.*

#### **3. The Multidimensional Case for** 1 < *ν* < 2

The Cauchy problem

$$\begin{cases} \frac{\partial^r u}{\partial t^j} = \sum\_{j=1}^d \frac{\partial^2 u}{\partial x\_j^2} & \mathbf{x} \in \mathbb{R}^d, \ t > 0 \\ u(\mathbf{x}, 0) = \boldsymbol{\delta}(\mathbf{x}) \\ u\_l(\mathbf{x}, 0) = 0. \end{cases} \tag{50}$$

has solution with Fourier transform

$$\int\_{\mathbb{R}^d} \boldsymbol{\varepsilon}^{j\gamma \cdot x} u(\mathbf{x}, t) \, \mathbf{d}x = E\_{\nu, 1}(-\left\|\gamma\right\|^2 \boldsymbol{t}^{\nu}) \tag{51}$$

$$= \frac{\sin \pi \nu}{\pi} \int\_0^\infty \frac{z^{\nu - 1} e^{-tz} \|\gamma\|^{2/\nu}}{z^{2\nu} + 2z^\nu \cos \pi \nu + 1} \, \mathbf{d}z + \frac{1}{\nu} \left[ \boldsymbol{\varepsilon}^{\|\gamma\|^{2/\nu}} t^{\bar{\boldsymbol{t}} \boldsymbol{t} \cdot \boldsymbol{\tau}} + \boldsymbol{\varepsilon}^{\|\gamma\|^{2/\nu} t \boldsymbol{t} - i \boldsymbol{\tau} \cdot \boldsymbol{\tau}} \right].$$

as shown in the analysis presented above. Thus, the solution to the Cauchy problem (50) reads

$$\begin{split} u(\mathbf{x},t) &= \frac{1}{(2\pi)^{d}} \int\_{\mathbb{R}^{d}} \mathbf{e}^{-i\gamma \cdot \mathbf{x}} \\ &\times \left\{ \frac{\sin \pi \nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1} e^{-tz \|\gamma\|^{2/\nu}}}{z^{2\nu} + 2z^{\nu} \cos \pi \nu + 1} \mathrm{d}z + \frac{1}{\nu} \left[ \boldsymbol{\varepsilon}^{\|\gamma\|^{2/\nu} t e^{ix/\nu}} + \boldsymbol{\varepsilon}^{\|\gamma\|^{2/\nu} t e^{-i\pi/\nu}} \right] \right\} \mathrm{d}\gamma \end{split} \tag{52}$$

We must therefore evaluate the following three *d*−dimensional integrals, the first one being a function of *z*.

$$\int\_{\mathbb{R}^d} e^{-i\boldsymbol{\gamma}\cdot\boldsymbol{x}} e^{-t||\boldsymbol{\gamma}||^{2/\nu}z} \, \mathrm{d}\boldsymbol{\gamma} \, \qquad \int\_{\mathbb{R}^d} e^{-i\boldsymbol{\gamma}\cdot\boldsymbol{x}} e^{||\boldsymbol{\gamma}||^{2/\nu}t e^{i\boldsymbol{x}/\nu}} \, \mathrm{d}\boldsymbol{\gamma} \, \qquad \int\_{\mathbb{R}^d} e^{-i\boldsymbol{\gamma}\cdot\boldsymbol{x}} e^{||\boldsymbol{\gamma}||^{2/\nu}t e^{-i\boldsymbol{x}/\nu}} \, \mathrm{d}\boldsymbol{\gamma} \tag{53}$$

Since the three integrals (53) are substantially similar, we restrict ourselves to the evaluation of the first one. In spherical coordinates we have that

$$\begin{split} & \quad \int\_{\mathbb{R}^d} e^{-i\boldsymbol{\rho}\cdot\mathbf{x}} e^{-\|\boldsymbol{\rho}\|^{\frac{2}{2}}\tau\_d} \mathbf{d}\boldsymbol{\rho} \\ &= \int\_0^\infty \rho^{d-1} \mathbf{d}\rho \int\_0^\pi \mathbf{d}\theta\_1 \cdots \int\_0^\pi \mathbf{d}\theta\_{d-2} \int\_0^{2\pi} \mathbf{d}\phi \sin\theta\_1^{d-2} \cdots \sin\theta\_{d-2} \\ & \quad \times e^{-i\boldsymbol{\rho}\left(\boldsymbol{x}\_d\sin\theta\_1\cdots\sin\theta\_{d-2}\sin\phi+\mathbf{x}\_{d-1}\sin\theta\_1\cdots\sin\theta\_{d-2}\cos\phi+\cdots+\mathbf{x}\_2\sin\theta\_1\cos\theta\_2+\mathbf{x}\_1\cos\theta\_1\right)} \mathbf{e}^{-t\boldsymbol{\rho}^{2/2}\boldsymbol{z}} \\ &= \int\_0^\infty \rho^{d-1} e^{-t\boldsymbol{\rho}^{2/2}\boldsymbol{z}} \frac{(2\pi)^{\frac{d}{2}}J\_{\frac{d}{2}-1}(\boldsymbol{\rho}\|\boldsymbol{\varepsilon}\|)}{\left(\rho\|\boldsymbol{\varepsilon}\|\|\right)^{\frac{d}{2}-1}} \mathbf{d}\boldsymbol{\rho} \end{split}$$

The last step is the hyperspherical integral

$$\int\_{\{\gamma\_1,\dots,\gamma\_d \le \sum\_{j=1}^d \gamma\_j^2 = \rho^2\}} \epsilon^{-i \sum\_{j=1}^d x\_j \gamma\_j} \mathrm{d}\gamma\_1 \dots \mathrm{d}\gamma\_d \tag{55}$$

$$\begin{split} \epsilon^{\pi} &= \int\_0^\pi \mathrm{d}\theta\_1 \cdots \cdot \int\_0^\pi \mathrm{d}\theta\_{d-2} \int\_0^{2\pi} \mathrm{d}\phi \sin\theta\_1^{d-2} \cdots \cdot \sin\theta\_{d-2} \\ &\times \epsilon^{-i\rho} (x\_d \sin\theta\_1 \cdots \sin\theta\_{d-2} \sin\phi + x\_{d-1} \sin\theta\_1 \cdots \sin\theta\_{d-2} \cos\phi + \cdots + x\_2 \sin\theta\_1 \cos\theta\_2 + x\_1 \cos\theta\_1) \end{split}$$

$$\begin{split} \epsilon^{\pi} &= \frac{(2\pi)^{\frac{d}{2}} I\_{\frac{d}{2}-1} (\rho \|\|\mathbf{x}\|)}{(\rho \|\|\mathbf{x}\|)^{\frac{d}{2}-1}} \end{split}$$

which is proven in detail in [13], Formula (2.151).

By inserting (54) into (52) we have that

$$\begin{split} u(\mathbf{x},t) &= \frac{1}{(2\pi)^{d}} \int\_{0}^{\infty} \rho^{d-1} \frac{(2\pi)^{\frac{d}{2}} I\_{\frac{d}{2}-1}(\rho \|\mathbf{x}\|)}{(\rho \|\mathbf{x}\|)^{\frac{d}{2}-1}} d\rho \\ &\times \left\{ \frac{\sin \pi\nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1} e^{-t\rho^{2/\nu}z}}{z^{2\nu} + 2z^{\nu}\cos \pi\nu + 1} d\mathbf{z} + \frac{1}{\nu} \left[ e^{\rho^{2/\nu}t\mathbf{z}^{i\pi/\nu}} + e^{\rho^{2/\nu}t\mathbf{z}^{-i\pi/\nu}} \right] \right\} \\ &= \frac{1}{(2\pi)^{\frac{d}{2}} \|\mathbf{x}\|^{\frac{d}{2}-1}} \int\_{0}^{\infty} \rho^{\frac{d}{2}} I\_{\frac{d}{2}-1}(\rho \|\mathbf{x}\|) \mathbf{d}\rho \\ &\times \left\{ \frac{\sin \pi\nu}{\pi} \int\_{0}^{\infty} \frac{z^{\nu-1} e^{-t\rho^{2/\nu}z}}{z^{2\nu} + 2z^{\nu}\cos \pi\nu + 1} d\mathbf{z} + \frac{1}{\nu} \left[ e^{\rho^{2/\nu}t\mathbf{z}^{i\pi/\nu}} + e^{\rho^{2/\nu}t\mathbf{z}^{-i\pi/\nu}} \right] \right\} \end{split} \tag{36}$$

Note that the integral in *z* after the change of variable *z<sup>ν</sup>* = *z* becomes

$$\frac{\sin \pi \nu}{\pi \nu} \int\_0^\infty \frac{e^{-t\rho^{2/\nu} z^{1/\nu}}}{z^2 + 2z \cos \pi \nu + 1} d\mathbf{z} = \frac{1}{\nu} \int\_0^\infty e^{-t\rho^{2/\nu} z^{1/\nu}} \left[ \frac{1}{\pi} \frac{\sin \pi \nu}{(z - \cos \pi \nu)^2 + \sin^2 \pi \nu} \right] d\mathbf{z} \tag{57}$$

since 1 < *ν* < 2 the function

$$f(z) = \frac{1}{\pi} \frac{\sin \pi \nu}{(z - \cos \pi \nu)^2 + \sin^2 \pi \nu} \tag{58}$$

has the form shown in Figure 9.

**Figure 9.** Plot of the function *f* given in (58).

In conclusion

$$\frac{1}{\nu} \int\_0^\infty f(z) \, \mathrm{d}z = \begin{cases} -1 & \frac{3}{2} < \nu < 2 \\ -1 + \frac{1}{\nu} & 1 < \nu < \frac{3}{2} \end{cases} \tag{59}$$

We recall now the definition of an *<sup>α</sup>*−stable *d*−dimensional process *<sup>S</sup><sup>α</sup>*(*t*)=(*S<sup>α</sup>*1 (*t*), ... , *Sαd* (*t*)), 0 < *α* < 2.

Its characteristic function has the following form [14]

$$\mathbb{E}\mathbf{z}^{i\gamma\cdot S^{a}(t)} = \begin{cases} \mathcal{e}^{-t} \int\_{\mathbb{S}^{d-1}} |\gamma \cdot s|^{a} (1 - i \operatorname{sign}(\gamma \cdot s) \tan\left(\frac{\mu}{\mu}\right) \Gamma(\mathrm{ds}) + i\gamma \cdot \mu & \mathfrak{a} \neq \mathbf{1} \\\\ \mathcal{e}^{-t} \int\_{\mathbb{S}^{d-1}} |\gamma \cdot s| (1 - i \frac{\mu}{\pi} \operatorname{sign}(\gamma \cdot s) \log(\gamma \cdot s)) \Gamma(\mathrm{ds}) + i\gamma \cdot \mu & \mathfrak{a} = \mathbf{1} \end{cases} \tag{60}$$

where *μ* ∈ <sup>R</sup>*d*, Γ is a finite measure on the sphere S*d*−<sup>1</sup> = {*s* ∈ R*<sup>d</sup>* : *s* = <sup>1</sup>}.

Since *γ* · *s* = *γs* cos *θ* = *γ* cos *θ*, where *s* ∈ S*d*−<sup>1</sup> so that *s* = 1. Furthermore sign (*γ* · *s*) = sign cos *θ*. We can assume *γ* oriented through the north pole of S*d*−<sup>1</sup> and thus *θ* can be viewed as the latitude of vector *s*, as shown in Figure 10.

**Figure 10.** Domain of integration in (60), upon suitable rotation of the axis of the sphere.

We take the case 1 < *α* < 2 and rewrite the characteristic function as

$$\begin{split} \mathbb{E}e^{i\gamma \cdot S^{a}(t)} &= e^{-t\left\|\gamma\right\|^{a} \int\_{\mathbb{S}^{d-1}} \cos^{a}\theta \left(\cos\frac{\pi a}{2} - i \operatorname{sign}(\cos\theta)\sin\frac{\pi a}{2}\right) \frac{\Gamma(d\nu)}{\cos\frac{\pi b}{2}} + i\gamma \cdot \mu} \\ &= e^{-t\left\|\gamma\right\|^{a} \sigma^{a} \int\_{\mathbb{S}^{d-1}} \cos^{a}\theta \epsilon^{-i\frac{\pi b}{2}\mathfrak{a}\sin(\cos\theta)} \Gamma(d\nu) + i\gamma \cdot \mu} \end{split} \tag{61}$$

where *σ<sup>α</sup>* = 1 cos *πα*2

.

 For simplicity we assume *d* = 3 and suppose that Γ is a uniformly distributed measure on the upper hemisphere of the unit sphere. Thus

$$
\Gamma(\text{ds}) = \frac{\sin \theta}{2\pi} \text{d}\theta \text{d}\phi \qquad 0 < \theta < \frac{\pi}{2}, 0 < \phi < 2\pi. \tag{62}
$$

with sign(cos *θ*) = 1. For *μ* = 0 the integral in (61) becomes

$$\frac{e^{-i\frac{\pi\mu}{2}}}{2\pi} \int\_{\mathbb{S}^2} \cos^a \sin \theta \,\mathrm{d}\phi \mathrm{d}\theta = \frac{e^{-i\frac{\pi\mu}{2}}}{2\pi} \int\_0^{\frac{\pi}{2}} \sin \theta \cos^a \theta \int\_0^{2\pi} \mathrm{d}\phi = \frac{e^{-i\frac{\pi\mu}{2}}}{a+1} \tag{63}$$

The characteristic function (61) turns out to be

$$\mathbb{E}e^{i\gamma \cdot S^{a}(t)} = e^{-t\left\|\gamma\right\|^{a} \frac{1}{\operatorname{csm}\frac{\mathsf{W}\mathsf{F}}{2}} \frac{e^{-i\frac{\mathsf{F}\mathsf{A}}{2}}}{a+1}} = e^{t\left\|\gamma\right\|^{\frac{1}{\mathsf{F}}} \frac{1}{\left\|\operatorname{csm}\frac{\mathsf{F}}{\mathsf{F}}\right\|^{a} + 1} = e^{t\left\|\gamma\right\|^{\frac{1}{\mathsf{F}}} \frac{1}{\left\|\operatorname{csm}\frac{\mathsf{F}}{\mathsf{F}}\right\|^{a}} e^{-i\frac{\mathsf{F}}{\mathsf{F}}} \sigma^{\mathsf{F}}}$$

Since *α* = 2/*ν* and 1 < *ν* < 2, we have *π* > *π*/*ν* > *π*/2 so cos *απ*2 = cos *πν* is negative. If Γ is distributed on the lower hemisphere sign(cos *θ*) = −1, and in the same way we have that

$$\mathbb{E}\boldsymbol{\varepsilon}^{i\gamma\cdot S^{a}(t)} = \boldsymbol{\varepsilon}^{-t||\gamma||^{a}\frac{1}{\cos\frac{\pi\theta}{2}}\frac{\boldsymbol{\varepsilon}^{\frac{\pi\pi a}{2}}}{a+1}} = \boldsymbol{\varepsilon}^{t||\gamma||\frac{\tilde{\boldsymbol{\theta}}}{|\cos\frac{\pi\theta}{2}|}\frac{1}{a+1}} = \boldsymbol{\varepsilon}^{t||\gamma||\frac{\tilde{\boldsymbol{\theta}}}{|\cos\frac{\pi\theta}{2}|}\boldsymbol{\varepsilon}^{i\frac{\pi}{\theta}}\boldsymbol{\sigma}^{a}}$$

The situation in the space <sup>S</sup>*d*−1, *d* > 3 is quite similar with the integral in (61) evaluated in hyperspherical coordinates.

**Author Contributions:** Conceptualization, F.I. and E.O.; Writing—original draft, F.I. and E.O.; Writing—review & editing, F.I. and E.O. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Conflicts of Interest:** The authors declare no conflict of interest.
