Eigenfunctions and Fundamental Solutions of the Fractional Laplace and Dirac Operators: The Riemann-Liouville Case

In this paper we study eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator Δ+(α,β,γ):=Dx0+1+α+Dy0+1+β+Dz0+1+γ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta _+^{(\alpha , \beta , \gamma )}:= D_{x_0^+}^{1+\alpha } +D_{y_0^+}^{1+\beta } +D_{z_0^+}^{1+\gamma },$$\end{document} where (α,β,γ)∈]0,1]3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\alpha , \beta , \gamma ) \in \,]0,1]^3$$\end{document}, and the fractional derivatives Dx0+1+α,Dy0+1+β,Dz0+1+γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{x_0^+}^{1+\alpha }, D_{y_0^+}^{1+\beta }, D_{z_0^+}^{1+\gamma }$$\end{document} are in the Riemann–Liouville sense. Applying operational techniques via two-dimensional Laplace transform we describe a complete family of eigenfunctions and fundamental solutions of the operator Δ+(α,β,γ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta _+^{(\alpha ,\beta ,\gamma )}$$\end{document} in classes of functions admitting a summable fractional derivative. Making use of the Mittag–Leffler function, a symbolic operational form of the solutions is presented. From the obtained family of fundamental solutions we deduce a family of fundamental solutions of the fractional Dirac operator, which factorizes the fractional Laplace operator. We apply also the method of separation of variables to obtain eigenfunctions and fundamental solutions.


Introduction
In the last decades the interest in fractional calculus increased substantially. This fact is due to on the one hand different problems can be considered in the framework of fractional derivatives like, for example, in optics and quantum mechanics, and on the other hand fractional calculus gives us a new degree of freedom which can be used for more complete characterization of an object or as an additional encoding parameter. For more details about fractional partial differential equations, their applications and their numerical solutions see [8] and the references indicated there.
The problems with the fractional Laplace attracted in the last years a lot of attention, due especially to their large range of applications. The fractional Laplace appears e.g. in probabilistic framework as well as in mathematical finance as infinitesimal generators of the stable Lévy processes [1]. One can find problems involving the fractional Laplace in mechanics and in elastostatics, for example, a Signorini obstacle problem originating from linear elasticity [4]. The concerning fluid mechanics and hydrodynamics to the nonlocal fractional Laplace appears, for instance, in the quasigeostropic fractional Navier-Stokes equation [3] and in the hydrodynamic model of the flow in some porous media [16].
The connections between fractional calculus and physics are, in some sense, relatively new but, and more important for the community, a subject of strong interest. In [14] the author proposed a fractional Dirac equation of order 2/3 and established the relation between the corresponding γ μ α −matrix algebra and generalized Clifford algebras. This approach was generalized in [19], where the author found that relativistic covariant equations generated by taking the n−th root of the d'Alembert operator are fractional wave equations with an inherent SU(n) symmetry. The study of the fractional Dirac operator is important due to its physical and geometrical interpretations. Physically, this fractional differential operator is related with some aspects of fractional quantum mechanics such as the derivation of the fractal Schrödinger type wave equation, the resolution of the gauge hierarchy problem, and the study of super-symmetries. Geometrically, the fractional classical part of this operator may be identified with the scalar curvature in Riemannian geometry.
Clifford analysis is a generalization of classical complex analysis in the plane to the case of an arbitrary dimension. At the heart of the theory lies the Dirac operator D, a conformally invariant first-order differential operator which plays the same role as the Cauchy-Riemann operator in complex analysis. In [10,17] the authors studied the connections between Clifford analysis and fractional calculus, however, the fractional Dirac operator considered in these works do not coincide with the one used here.
The aim of this paper is to present an explicit expression for the family of eigenfunctions and fundamental solutions of the three-parameter fractional Laplace operator, as well as, a family of fundamental solutions of the fractional Dirac operator. For the sake of simplicity we restrict ourselves to the three dimensional case, however the results can be generalized for an arbitrary dimension. The two dimensional case was studied in [18] without considering the connections with Clifford analysis. The authors would like to point out that the fractional Laplace operator considered in this paper is different from the fractional Laplace operator defined via Fourier transform (see [8]). The deduction of the fundamental solution for the fractional Dirac operator defined via Riemann-Liouville derivatives is a completely new result in the context of fractional Clifford analysis. The fundamental solutions of the fractional Dirac operator obtained in this paper are the basis to develop an operator calculus theory in the context of fractional Clifford analysis.
The structure of the paper reads as follows: in the Preliminaries we recall some basic facts about fractional calculus, special functions and Clifford analysis, which are necessary for the development of this work. In Sect. 3 we use operational techniques for the two dimensional Laplace transformation and its extension to generalized functions to describe a complete family of eigenfunctions and fundamental solutions of the fractional Laplace operator. In the same section we compute the family of fundamental solutions for the fractional Dirac operator. In Sect. 4 we obtain the analogous of the results of Sect. 3, but via the method of separation of variables. Particular cases of solutions can be obtained using the obtained generic formulas and considering α = β = γ = 1.

Fractional Calculus and Special Functions
Let D α a + f (x) denote the fractional Riemann-Liouville derivative of order α > 0 (see [11]) where n = [α] + 1 and [α] means the integer part of α. When 0 < α < 1 then (2.1) takes the form (2. 2) The Riemann-Liouville fractional integral of order α > 0 is given by (see [11]) We denote by I α a + (L 1 ) the class of functions f represented by the fractional integral (2.3) of a summable function, that is f = I α a + ϕ, ϕ ∈ L 1 (a, b). A description of this class of functions was given in [15].   [15]) A function f ∈ L 1 (a, b) has a summable fractional deriva- If a function f admits a summable fractional derivative, then the composition of (2.1) and (2.3) can be written in the form (see [15,Thm. 2.4]) . Nevertheless we note that D α a + I α a + f = f in both cases. This is a particular case of a more general property (cf. [13, (2.114 It is important to remark that the semigroup property for the composition of fractional derivatives does not hold in general (see [13,Sect. 2.3.6]). In fact, the property holds whenever One important function used in this paper is the two-parameter Mittag-Leffler function E μ,ν (z) [7], which is defined in terms of the power series by In particular, the function E μ,ν (z) is entire of order ρ = 1 μ and type σ = 1. The exponential, trigonometric and hyperbolic functions are expressed through (2.8) as follows (see [7]): Two important fractional integral and differential formulae involving the twoparametric Mittag-Leffler function are the following (see [7, p. 61, p. 87]) The formal approach presented in Sects. 3 and 4 based on the Laplace transform leads to the solution of a linear Abel integral equation of the second kind.
has a unique solution

Clifford Analysis
Let {e 1 , . . . , e d } be the standard basis of the Euclidean vector space in R d . The associated Clifford algebra R 0,d is the free algebra generated by R d modulo x 2 = −||x|| 2 e 0 , where x ∈ R d and e 0 is the neutral element with respect to the multiplication operation in the Clifford algebra R 0,d . The defining relation induces the multiplication rules e i e j + e j e i = −2δ i j , (2.12) where δ i j denotes the Kronecker's delta. In particular, e 2 i = −1 for all i = 1, . . . , d. The standard basis vectors thus operate as imaginary units. A vector space basis for R 0,d is given by the set where e A = e l r e l r −1 . . . e l 1 , and e j = −e j for j = 1, . . . , d, e 0 = e 0 = 1. An important subspace of the real Clifford algebra R 0,d is the so-called space of paravectors R d 1 = R R d , being the sum of scalars and vectors. Each non-zero vector a ∈ R d 1 has a multiplicative inverse given by a ||a|| 2 . Clifford analysis can be regarded as a higher-dimensional generalization of complex function theory in the sense of the Riemann approach. An R 0,d −valued function f over ⊂ R d 1 has the representation f = A e A f A , with components f A : → R 0,d . Properties such as continuity or differentiability have to be understood componentwise. Next, we recall the Euclidean Dirac operator D = d j=1 e j ∂ x j , which factorizes the For more details about Clifford algebras and basic concepts of its associated function theory we refer the interested reader for example to [5,9].

Eigenfunctions and Fundamental Solution of the Fractional Laplace Operator
We consider the eigenfunction problem for the fractional Laplace operator Taking the integral operator I 1+α x + 0 from both sides of (3.1) and taking into account (2.4) we get Applying the integral operator I 1+β y + 0 to both sides of the previous expression and using Fubini's Theorem we get where we denote the Cauchy's fractional integral conditions by Finally, we apply I 1+γ z + 0 to both sides of equation (3.2) and we use again Fubini's Theorem to get which is equivalent to where we denote the Cauchy's fractional integral conditions by We observe that the fractional integrals in (3.4) are Laplace-transformable functions. Therefore, we may apply the two-dimensional Laplace transform to y and z: Taking into account its convolution and operational properties [6,11] we obtain the following relations: Proceeding in a similar way we obtain the Laplace transform of the remaining terms of (3.4). Combining all the resulting terms and multiplying by s we obtain the following second kind homogeneous integral equation of Volterra type: where Using (2.11) we have that the unique solution of (3.7) in the class of summable functions is which involves as the kernel the two-parameter Mittag-Leffler function (see (2.8)). Due to the convergence of the integrals and series that appear in (3.9), we can interchange them and rewrite (3.9) in the following way: In order to cancel the Laplace transform we need to take into account its distributional form in Zemanian's space (for more details about generalized integral transforms see [20]) and the following relations: where i = 0, 1, n ∈ N 0 , δ is Dirac's delta function, and the convergence is in D . Therefore, applying the multinomial theorem and after straightforward calculations we get the following family of eigenfunctions of (3.1) where the convergence of the series is in D . From the previous calculations we obtain the following theorem, where we describe the eigenfunctions in an operational form using the Mittag-Leffler function (2.8).
Proof We give a direct proof of the theorem. It is based on the fact that D 1+α We use also the fractional analogous formula for differentiation of integrals depending on a parameter where the upper limit also depends on the same parameter (see [ Rearranging the terms of the series we obtain We would like to remark that the eigenfunctions constructed in [18] for the two-parameter Laplace operator have a different structure and do not satisfy the eigenfunction equation, except when λ = 0. A family of generalized eigenfunctions for the fractional Laplace operator can be obtained considering λ = 0 as the next theorem states.

Theorem 3.2 The generalized fundamental solution of the fractional Laplace operator
(α,β,γ ) + is given by: (3.12) In a similar way, applying in (3.4) the two-dimensional Laplace transform with respect to x and y we obtain the following generalized fundamental solutions v 0 (x, y, z) = (y − y 0 ) β−1 E 1+β,β −(y − y 0 ) 1+β D 1+α Furthermore, a similar result can be obtained when we apply the two-dimensional Laplace transform with respect to x and z.

Remark 3.3
It is possible to obtain from (3.12) the fundamental solution for the Euclidean Laplace operator when α = β = γ = 1. Since the fundamental solution of the Euclidean Laplace operator in R 3 is given by which corresponds to the following convergent power series (y−y 0 ) 2 +(z−z 0 ) 2 < 1, we have to consider in (3.12) α = β = γ = 1, f 1 , h 1 , g 1 , h 0 , and g 0 the null functions, and f 0 such that The function f 0 (y, z) = (y − y 0 ) 2 + (z − z 0 ) 2 − 1 2 satisfies (3.15). To see this we recall that the p-dimensional Euclidean Laplace satisfies (see with r = x , x ∈ R p , n ∈ N. Therefore, for p = 2 and k = −1 we obtain which leads immediately to (3.15)

Fundamental Solution of the Fractional Dirac Operator
In this section we compute the fundamental solution for the three dimensional fractional left Dirac operator defined by This operator factorizes the fractional Laplace operator (α,β,γ ) + for Clifford valued functions f given by f (x, y, z) = A e A f A (x, y, z), where e A ∈ {1, e 1 , e 2 , e 3 , e 1 e 2 , e 1 e 3 , e 2 e 3 , e 1 e 2 e 3 }, f A ∈ AC 2 ( ) and In fact, for such functions we can apply the semigroup property (2.6) to obtain We remark that in the case of our fractional Dirac operator, condition (2.7) reduces If f ∈ AC 2 ( ) doesn't satisfy this last property it is always possible to define g(x, y, z) In a similar way we conclude that From (3.18), (3.19), (3.20), and the multiplication rules (2.12) of the Clifford algebra, we finally get i.e., the fractional Dirac operator factorizes the fractional Laplace operator. We remark that the factorization property (3.21) also holds for functions which are continuous on (x 0 , X 0 ] and integrable on any subinterval of [x 0 , X 0 ] since the semigroup property (2.6) remains valid in this case (see [12]). In order to get the fundamental solution of D (α,β,γ ) + we apply this operator to the fundamental solution (3.12). In the following theorem, to compute D 1+α 2 x + 0 u we make use of the derivation rule (2.10) and the fractional analogous formula for differentiation of integrals depending on a parameter where the upper limit also depends on the same parameter (see [

Eigenfunctions and Fundamental Solution of the Fractional Laplace Operator
Let us consider again equation (3.1) D 1+α and assume that u(x, y, z) = u 1 (x) u 2 (y) u 3 (z). Substituting in (3.4) and taking into account the initial conditions (3.3), (3.5), and (3.6) we obtain where a i = f i , b i = h i , c i = g i ∈ C, i = 1, 2, are constants defined by the initial conditions (3.3), (3.5), and (3.6). Supposing that I 1+α for (x, y, z) ∈ , we can divide (4.1) by this factor. Separating the variables we get the following three Abel's type second kind integral equations: (4.4) where λ, μ, ν ∈ C are constants. We observe that the equality I 1+α for at least one point (ξ, η, θ) agrees with (4.1), (4.2), (4.3), and (4.4). Solving the latter equations using (2.11) in Theorem 2.3 and after straightforward computations we obtain the following family of eigenfunctions u λ,μ,ν (x, y, z) = u 1 (x) u 2 (y) u 3 (z), Remark 4.1 In the special case of α = β = γ = 1 the functions u 1 , u 2 and u 3 take the form which are the components of the fundamental solution of the Laplace operator in R 3 obtained by the method of separation of variables.

Fundamental Solution of the Fractional Dirac Operator
Following the procedure presented in Sect. 3.2, the fundamental solution for the fractional Dirac operator D