Desiderata for Fractional Derivatives and Integrals

: The purpose of this brief article is to initiate discussions in this special issue by proposing desiderata for calling an operator a fractional derivative or a fractional integral. Our desiderata are neither axioms nor do they deﬁne fractional derivatives or integrals uniquely. Instead they intend to stimulate the ﬁeld by providing guidelines based on a small number of time honoured and well established criteria.

A list of six desiderata 1 is proposed that in our opinion would justify calling an operator D α (or I α ) a fractional derivative (or a fractional integral) of non-integer order α / ∈ N. Derivatives and integrals of fractional order have a long history and, up until the recent proliferation of novel fractional derivatives, most definitions and interpretations of fractional operators seem to implicitly assume the desiderata of an operational calculus as formulated in this article.
Mathematical terms used in the formulation of our desiderata are defined in Appendix A. A family {D α , I α } of operators with α ∈ Q, R or C is proposed to be called a family of fractional derivatives D α and integrals I α of order α (with Re α ≥ 0) 2 if and only if it satisfies the following six desiderata: holds true for Re α ≥ 0 and Re β ≥ 0, where D(I α ) denotes the domain of I α , and • denotes composition of operators. (c) Restricted to a suitable subset G (c) ⊂ D(I α ) of the domain of I α the fractional derivatives D α of order α operate as left inverses for all α with Re α ≥ 0, where 1 G (c) is the identity on G (c) . 1 properties to be desired. 2 It is common to use only one of the symbols I or D in the sense that either D α = I −α or I α = D −α . In this paper we keep the distinction between I and D by assuming Re α ≥ 0 unless otherwise specified. This entails discussing the case Re α = 0 separately whenever necessary.
(d) There is a subset G (d) ⊂ D(D α ) of the domain of D α such that the limits exist in some sense and define linear maps D 1 : This means it is possible to define a multiplication · : holds for all f , g ∈ G (d) .
If the semigroup law (1) can be extended to all α ∈ R or α ∈ C, we propose to speak of fractional calculus. Our desiderata are obviously inspired by operational calculus. Recall that an operational calculus is a continuous one-to-one mapping between an algebra of functions and an algebra of operators such that the neutral elements match and algebraic relations are preserved. Extending the algebra from polynomial functions to convergent power series suffices for an operational calculus. More singular functions, namely non-analytic power functions, are required for fractional calculus.
Desiderata differing substantially from those above have been formulated in [1] (p. 5) and [2] (p. 5). Envisaging exclusively analytic functions the criteria given in [1] (p. 5) are extremely restrictive. In theory and applications it is nowadays imperative to include more general functions, measures and also distributions into the purview.
Given the extreme restrictions in [1] (p. 5) a more recent proposal [2] went to opposite extremes. Little or no attention is given to a domain of definition for the fractional derivatives in [2]. Our desiderata for fractional derivatives again differ substantially from those in [2]. Rather than requiring some form of the generalized Leibniz rule 2 P5 in [2] for all α we desire the Leibniz rule only for α = 1, and that differs not only from 2 P5, but also from 2 P3. In addition the identity rule 2 P2 in [2] does not restrict the admissible operators at all. As long as there is no continuity in α or a well defined limit, the identity rule can always be fulfilled, simply by setting D 0 = 1. More generally, Ref. [2] seems to neglect parameters other than α, or the topological and operator-theoretic implications of the limit in Equation (2) [2].
Our desiderata do not include non-locality of fractional derivatives. Fractional derivatives, that are local operators, were introduced in [3,4] and are discussed further in [5] and Section 7 of [6]. Contrary to 2 P3 in [2] we do not constrain the limits α → n with n ∈ N for n ≥ 2, because we wish to allow more generality.
To illustrate our desiderata in a simple case consider fractional operators of Riemann-Liouville type for complex valued functions f : when the fractional order α is real and restricted to 0 ≤ α ≤ 1. The (right-sided) Riemann-Liouville fractional integrals I α a+ of order α ≥ 0 with lower limit a are defined by setting I 0 for α > 0 and x ≤ b. The (right-sided) Riemann-Liouville-type fractional derivatives D α a+ of order 0 ≤ α ≤ 1 are defined as [5] (p. 434) where the number 0 ≤ β ≤ 1 parametrizes different types of fractional derivatives. The classical Riemann-Liouville derivative is of type β = 0, while the popular Liouville-Caputo derivative has type β = 1 [7] (p. 10). Both operator families are linear and desideratum (a) can be fulfilled on numerous linear spaces due to the compactness of [a, b]. Examples are the Lebesgue spaces The Riemann-Liouville fractional integrals (extended from 0 ≤ α ≤ 1 to α ≥ 0) are a strongly continuous semigroup of operators with respect to the parameter α ≥ 0, and obey the index law (1) ([a, b]) or suitable subspaces. The desiderata (d), (e) and (f) can then be derived with the help of the semigroup property. For the Riemann-Liouville operators they hold e.g. for smooth (infinitely often differentiable) functions f ∈ G For infinite intervals or for generalized functions the problem of domains may become more involved and our desiderata may become more restrictive. As an example consider the family of symmetric Riesz operators on the real line. In this case the limiting operator D 1 = lim α→−1 I α is again well defined, but does not fulfill desideratum (f). Instead, it fulfills Leibniz' formula for D n ( f g) with n = 2, i.e., We propose to call such operators obeying Leibniz formula for D n ( f g) with n ≥ 2 pseudofractional derivatives or fractional pseudoderivatives.
Of course, our desiderata do not define fractional derivatives and integrals in a unique way. Still, they considerably restrict the set of admissible operators as seen above. In our opinion the above desiderata formulate crucial constraints for the development of a meaningful mathematical theory of fractional calculus and its reasonable applications.
Much work has been done on mathematical interpretations of fractional derivatives and integrals. The results are documented in numerous texts and treatises (see [6,8] for recent reviews). It seems however, that the connection (or not) of classical and recent fractional calculi with historical and contemporary forms of operational and functional calculi such as Heaviside-Mikusinski calculus, Dunford-Schwarz calculus, or Hille-Phillips calculus is a rich source of numerous open problems whose speedy solution would seem pertinent to advance and ultimately consolidate the field. We hope that the desiderata above are sufficiently restrictive to initiate a discussion of these pressing problems, and thereby stimulate readers and contributors to address some of these open problems in their areas of expertise and interest.

Appendix A
For the convenience of readers from non-mathematical disciplines we recall some definitions: A real (or complex) linear space (or vector space) over the field R (or C) of real (or complex) numbers is a non-empty set X with two operations called addition and scalar multiplication fulfilling the usual rules of vector addition and multiplication of vectors with numbers 5 . 5 (a) for all f , g ∈ X also f + g ∈ X, (b) f + g = g + f , (c) f + (g + h) = ( f + g) + g, (d) there exists an element 0 ∈ X (called origin ) such that f + 0 = f for all f ∈ X, (e) for all f ∈ X there is an element − f ∈ X such that f + (− f ) = 0 (f) for all a ∈ R (or a ∈ C) and f ∈ X an element a f ∈ X is defined, (g) for all a ∈ R (or a ∈ C) and f , g ∈ X one has a( f + g) = a f + ag, (h) for all a, b ∈ R (or a, b ∈ C) and f ∈ X one has (a + b) f = a f + b f , (i) for all a, b ∈ R (or a, b ∈ C) and f ∈ X one has a(b f ) = (ab) f , and (j) 1 f = f for all f ∈ X.
Let X, Y, Z be linear spaces (vector spaces). A linear operator A : X → Y is a linear subspace of the direct sum X ⊕ Y, where is the linear space of pairs ( f , g) with f ∈ X, g ∈ Y and addition defined as ( f 1 , g 1 ) + ( f 2 , g 2 ) = ( f 1 + f 2 , g 1 + g 2 ) for all f i ∈ X and g i ∈ Y. The identity operator 1 X : X → X is defined as The domain D(A) and range R(A) of a linear operator A : X → Y are The inverse A −1 of A is defined as and is a Banach space for the norm For α = k + γ > 1 with k = 0, 1, 2, ... and 0 < γ ≤ 1 it is defined as where f (k) is the k-th derivative of f . A family {T(t)}, t ≥ 0 of bounded linear operators T(t) : X → X on a Banach space X is called a strongly continuous one-parameter semigroup if it satisfies: