General Fractional Integrals and Derivatives of Arbitrary Order

: In this paper, we introduce the general fractional integrals and derivatives of arbitrary order and study some of their basic properties and particular cases. First, a suitable generalization of the Sonine condition is presented, and some important classes of the kernels that satisfy this condition are introduced. Whereas the kernels of the general fractional derivatives of arbitrary order possess integrable singularities at the point zero, the kernels of the general fractional integrals can—depending on their order—be both singular and continuous at the origin. For the general fractional integrals and derivatives of arbitrary order with the kernels introduced in this paper, two fundamental theorems of fractional calculus are formulated and proved. ﬁrst fundamental theorem of fractional calculus; second fundamental theorem of fractional calculus


Introduction
In his papers [1,2], Abel derived and studied a mathematical model for the tautochrone problem in the form of the following integral equation (with slightly different notations): In fact, he considered the even more general integral equation under an implicit restriction 0 < α < 1. It is easy to see that the right-hand side of (2) is the operator that is currently referred to as the Caputo fractional derivative * D α 0+ of the order α, 0 < α < 1. Abel's solution formula to Equation (2) is nothing else than the operator now called the Riemann-Liouville fractional integral I α 0+ of the order α > 0: In modern notation, Formulas (2) and (3) correspond to the second fundamental theorem of FC for the Caputo fractional derivative of a function that takes the value zero at the point zero: where the validity of the condition φ(0) = 0 follows from the construction of Abel's mathematical model for the tautochrone problem. For more details regarding Abel's results and derivations presented in [1,2], see the recent paper presented in [3]. To solve the integral Equation (2), in [2], published in 1826, Abel employed the relation where the operation * stands for the Laplace convolution, and {1} is the function that is identically equal to 1 for t ≥ 0.
In [4], published in 1884, Sonine recognized that the relation (5) is the most crucial ingredient of Abel's solution method that can be generalized and applied to the analytical treatment of a larger class of integral equations. In place of (5), Sonine considered a pair of functions κ, k (Sonine kernels) that satisfy the relation In what follows, we denote the set of the Sonine kernels by S. For a given Sonine kernel κ, the kernel k that satisfies the Sonine condition (7) is called its associate Sonine kernel. Following Abel's solution method, Sonine showed that the integral equation has a solution in the form provided that the kernels κ, k satisfy the Sonine condition (7). Indeed, we obtain which immediately leads to the Formula (9). Of course, any concrete realization of the Sonine schema requires a precise characterization of the Sonine kernels and the spaces of functions where the operators from the right-hand sides of (8) and (9) are well defined. In [4], Sonine introduced a large class of the Sonine kernels in the form where the functions κ 1 = κ 1 (t), k 1 = k 1 (t) are analytical on R and their coefficients are connected by the relations The most prominent pair of the Sonine kernels from this class is given by the formulas where are the Bessel and the modified Bessel functions, respectively. Later, the evolution equations with the integro-differential operators of the convolution type (compare them with the solution by Sonine in Formula (9)), were actively studied in the framework of the abstract Volterra integral equations on the Banach spaces (see [5] and references therein). For example, in [6], the case of operators with the completely positive kernels k ∈ L 1 (0, +∞) was considered. The kernels from this class satisfy the condition (compare it with the Sonine condition (7)) where a ≥ 0 and l ∈ L 1 (0, +∞) is a non-negative and non-increasing function. However, until recently, no interpretation of these general results in the framework of fractional calculus (FC) had been suggested. The situation changed with the publication of the paper presented in [7] (see also [8][9][10]). In [7], Kochubei introduced a class K of kernels that satisfy the following conditions: exists for all p > 0; (K2)k(p) is a Stieltjes function (see [11] for details regarding the Stieltjes functions); (K3)k(p) → 0 and pk(p) → +∞ as p → +∞; (K4)k(p) → +∞ and pk(p) → 0 as p → 0.
Using the technique of the complete Bernstein functions, Kochubei investigated the integro-differential operators in the form of (15) and their Caputo type modifications with the kernels from K. In [7], Kochubei showed the inclusion K ⊂ S, introduced the corresponding integral operator and proved the validity of the first fundamental theorem of FC; i.e., that the operators (15) and (18) are left-inverse to the integral operator (19) on the suitable spaces of functions. Moreover, Kochubei treated some basic ordinary and partial fractional differential equations with the time-derivative in the form of (18) and proved that the solution to the Cauchy problem for the relaxation equation with the operator (18) and a positive initial condition is completely monotonic and that the fundamental solution to the Cauchy problem for the fractional diffusion equation with the time-derivative in the form of (18) can be interpreted as a probability density function. These results justified calling the operators (15) and (18) the general fractional derivatives (GFDs) in the Riemann-Liouville and Caputo sense, respectively. The integral operator (19) was called the general fractional integral (GFI).
The GFDs (15) and (18) with the kernels k ∈ K ⊂ S possess a series of important properties. However, the conditions (K1)-(K4) are very strong (especially the condition (K2)), and thus in subsequent publications, the operators (15) and (18) with the Sonine kernels from some larger classes were considered from the viewpoint of FC and its applications. In [12], a class of the kernels was introduced that ensures the validity of a maximum principle for the general time-fractional diffusion equations with the operators of type (18). Another important class of the Sonine kernels was described in [8] in terms of the completely monotone functions. As shown in [8], any singular (unbounded in a neighborhood of the point zero) locally integrable completely monotone function κ is a Sonine kernel, and its associate kernel k is also a locally integrable completely monotone function.
In the recent publications presented in [9,13], the operators (15) and (18) with the Sonine kernels from the class S −1 ⊂ S that satisfy only some minimal restrictions were studied from the viewpoint of FC. The Sonine kernels κ, k ∈ S −1 are continuous on R + and possess the integrable singularities of the power function type at the point zero.
In particular, in [9], the first and second fundamental theorems of FC for the operators (15) and (18) with the kernels k ∈ S −1 were formulated and proved. In [13], an operational calculus of the Mikusiński type for the operators (18) with the Sonine kernels k ∈ S −1 was constructed and applied for the analytical treatment of some initial value problems for the fractional differential equations with these operators.
It is clear that weakening the Kochubei conditions (K1)-(K4) on the Sonine kernels from K leads to the abandonment of some properties that were derived in [7] for the GFDs (15) and (18). However, it was shown in [9,13] that the operators (15) and (18) with the Sonine kernels k ∈ S −1 and the corresponding integral operator (19) still satisfy the main properties that the fractional derivatives and integrals should fulfill (see [14] and the references therein). Thus, these operators can also be interpreted as the GFDs and GFIs.
Another important point concerns the "generalized order" of the GFDs (15) and (18) with the Sonine kernels from the classes mentioned above. While projecting these operators to the conventional Riemann-Liouville and Caputo fractional derivatives (the case of the kernel k(t) = h 1−α (t)), the derivatives' order is restricted only to the case of α ∈ (0, 1). The reason is that the Sonine condition (5) for the power functions h α and h 1−α holds true only in the case 0 < α < 1. Moreover, even in the definition of the Caputo type general fractional derivative (18), only one initial condition is contained, which again indicates that the "generalized order" of this operator does not exceed one.
Because the Riemann-Liouville fractional integral and the Riemann-Liouville and Caputo fractional derivatives are defined for arbitrary order α ≥ 0, an extension of the GFDs (15) and (18) to the case of arbitrary order is worthy of investigation.
In a recent paper [9], the n-fold GFIs and GFDs were introduced as an attempt to extend their order behind the interval (0, 1). For example, the two-fold general fractional derivative constructed for the operator (15) with the kernel κ(t) = h 1−α (t), 0 < α < 1 is the Riemann-Liouville fractional derivative of the order 2α: Thus, we cannot ensure that the order of this two-fold GFD is always greater than one. Depending on the values of α and n, the "generalized order" of the n-fold GFD can be any number in the interval (0, n).
The main objective of this paper is to introduce the GFIs and GFDs of an arbitrary order in analogy to the Riemann-Liouville fractional integral and the Riemann-Liouville and Caputo fractional derivatives. This is done by a suitable generalization of the Sonine condition (7) and by the corresponding adjustment of Formulas (15) and (18), which define the GFDs in the Riemann-Liouville and Caputo senses.
The rest of the paper is organized as follows. In Section 2, following [9,13], we provide some basic definitions and properties of the GFDs (15) and (18) with the Sonine kernels k ∈ S −1 . Section 3 presents our main results. First, a suitable generalization of the Sonine condition (7) is introduced and some examples of the kernels that satisfy this condition are discussed. Then, the GFDs of an arbitrary order with these kernels are defined and their properties are studied. The conventional Riemann-Liouville and Caputo fractional derivatives of arbitrary order are particular cases of these GFDs. Another important example is the integro-differential operators of convolution type with the Bessel and the modified Bessel functions in the kernels. The constructions introduced in this section allow the formulation of the fractional differential equations with the GFDs with a generalized order greater than one with several initial conditions.

General Fractional Integrals and Derivatives with the Sonine Kernels
In this section, we provide some basic definitions and results regarding the GFIs and GFDs with the Sonine kernels from the class S −1 introduced in [9]. For more details, other relevant results and the proofs, see [9,13].
For our aims, we also need another two-parameter family of sub-spaces of C α (0, +∞) that allows us to better control the behavior of the functions at the origin: In particular, the sub-space C −1,0 (0, +∞) contains the functions that are continuous on R + and possess the integrable singularities of the power function type at the origin.
As mentioned in [17] (see also [8]), any Sonine kernel has an integrable singularity at the point zero. On the other hand, the kernels of the fractional integrals and derivatives should be singular [18]. Thus, the fractional integrals and derivatives with the Sonine kernels are worthy of investigation. In what follows, we consider the GFI (19) and the GFDs (15) and (18) of the Riemann-Liouville and Caputo types, respectively, with the Sonine kernels κ and k that belong to the sub-space C −1,0 (0, +∞) of the space C −1 (0, +∞). Definition 1. Let κ, k ∈ C −1,0 (0, +∞) be a pair of the Sonine kernels; i.e., let the Sonine condition (7) be fulfilled. The set of such Sonine kernels is denoted by S −1 : Several important features of the GFI (19) on the space C −1 (0, +∞) follow from the well-known properties of the Laplace convolution. In particular, we mention the mapping property the commutativity law and the index law that are valid on the space C −1 (0, +∞). Let κ ∈ S −1 and k be its associate Sonine kernel. The GFDs of the Riemann-Liouville and the Caputo types associated to the GFI (19) are given by the Formulas (15) and (18), respectively. It is easy to see that the GFD (18) in the Caputo sense can be rewritten as a regularized GFD (15) in the Riemann-Liouville sense: For the functions from C 1 −1 (0, +∞), the Riemann-Liouville GFD (15) can be represented as which immediately leads to the useful representation of the Caputo type GFD (18) that is valid on the space C 1 −1 (0, +∞). In the rest of this section, we formulate the first and second fundamental theorems of FC for the GFDs in the Riemann-Liouville and Caputo senses.

Theorem 1 (First Fundamental Theorem for the GFD).
Let κ ∈ S −1 and k be its associate Sonine kernel.
As shown in [9], the space C −1,(k) (0, +∞) can be also characterized as follows: Now, we proceed with the second fundamental theorem of FC for the GFDs in the Riemann-Liouville and Caputo senses.
Theorem 2 (Second Fundamental Theorem for the GFD). Let κ ∈ S −1 and k be its associate Sonine kernel.
In [9,13], the n-fold GFIs and GFDs with the Sonine kernels from S −1 were introduced and studied. For more details, we refer interested readers to these publications.

General Fractional Integrals and Derivatives of Arbitrary Order
As already mentioned in the Introduction, the "generalized order" of the GFIs and GFDs introduced so far is restricted to the interval (0, 1). The order of the n-fold GFIs and GFDs recently introduced in [9] belongs to the interval (0, n). However, it is hardly possible to fix their order between two neighboring natural numbers as in the case of the conventional Riemann-Liouville and Caputo fractional derivatives and thus to study, for example, the fractional oscillator equations or the time-fractional diffusion-wave equations with the GFDs of the order from the interval (1,2).
In this section, we define the GFIs and GFDs of arbitrary order and study their basic properties. As in the case of the conventional Riemann-Liouville and Caputo fractional derivatives, for the GFDs, we also have to distinguish between two completely different cases; namely, between the case of the integer order and the case of non-integer order. In the first case, the conventional Riemann-Liouville and Caputo fractional derivatives are defined as the integer-order derivatives, while in the second case, they are non-local integro-differential operators. Because the conventional Riemann-Liouville and Caputo fractional derivatives are important particular cases of the GFDs, we have no other choice but to follow the same strategy; namely, to separately define the GFDs of integer order as the integer-order derivatives and the GFDs of non-integer order as some integro-differential operators. In what follows, we focus on the case of the GFDs of non-integer order (the integer-order GFDs are simply the integer-order derivatives).
To introduce the GFIs and the GFDs of arbitrary non-integer order, we first formulate a condition on their kernels that generalizes the Sonine condition (7): Evidently, the Sonine condition corresponds to the case n = 1 of the more general condition (35).
Another important ingredient of our definitions is a set of the kernels that satisfy the condition (35) and belong to the suitable spaces of functions.
The set of pairs (κ, k) of such kernels is denoted by L n .

Remark 1.
The set L 1 coincides with the set of the Sonine kernels S −1 discussed in the previous section (see Definition 1). Indeed, in this case, the kernel κ ∈ C −1 (0, +∞) is a Sonine kernel, and therefore it has an integrable singularity at the point zero. Thus, it belongs to the subspace C −1,0 (0, +∞) as required in Definition 1.

Remark 2.
For n > 1, Definition 2 is not symmetrical with respect to the kernels κ and k because of the non-symmetrical inclusions κ ∈ C −1 (0, +∞) and k ∈ C −1,0 (0, +∞) (in the case n = 1, Definition 1 is symmetrical and one can interchange the kernels κ and k). However, the same statement is valid for the kernel κ(t) = h α (t), α > 0 of the Riemann-Liouville integral I α 0+ and the kernel k(t) = h n−α (t) of the Riemann-Liouville and Caputo fractional derivatives of order α, n − 1 < α < n, n ∈ N, defined as follows: with I α 0+ being the Riemann-Liouville fractional integral of order α: The solution to defining the integer-order Riemann-Liouville and Caputo fractional derivatives consists of a separate definition of the Riemann-Liouville fractional integral of the order α = 0: Of course, the definition (39) is not arbitrary and is justified inter alia by the formula that is valid for f ∈ L 1 (0, T) in every Lebesgue point of f; i.e., almost everywhere on the interval (0, T), T > 0 (see, e.g., [19]).

Example 1.
The kernels κ(t) = h α (t), α > 0 and k(t) = h n−α (t), n − 1 < α < n, n ∈ N provide a first example of the kernels from L n . Please note that the power functions h α and h n−α build a pair of the Sonine kernels only in the case n = 1; i.e., only in the case when the fractional derivatives' order is less than one.
Because both the Sonine condition (7) and its generalization (35) contain the Laplace convolution of two kernels, it is very natural to transform them into the Laplace domain. Providing that the Laplace transformsκ,k of the functions κ and k exist, the convolution theorem for the Laplace transform leads to the relatioñ for the Laplace transforms of the Sonine kernels and to a more general relatioñ for the kernels from the set L n introduced in Definition 2.

Example 2.
Formula (42) along with the works in [20,21] for the direct and inverse Laplace transforms, respectively, can be used to deduce other nontrivial examples of the kernels from L n . For instance, we employ the Laplace transform formulas (see [20]) for the Bessel function J ν and the modified Bessel function I ν defined by the power series (14) to introduce the kernels These kernels satisfy the condition (42). Moreover, for n − 2 < ν < n − 1, n ∈ N, the inclusions κ ∈ C −1 (0, +∞) and k ∈ C −1,0 (0, +∞) hold true, and thus the pair of the kernels (κ, k) given by (43) is from L n . Now let us consider a pair of the Sonine kernels (κ, k) from L 1 (in [4,8,9,13,17] and other related publications, many pairs of such kernels were presented). There are at least two reasonable possibilities to construct a pair (κ n k n ) of the kernels from L n , n > 1 based on the Sonine kernels κ, k from L 1 .
The first strategy consists of building the kernels κ n = κ n and k n = k n . Evidently, the kernels κ n and k n satisfy the relation (35) because κ and k are the Sonine kernels: However, the pair (κ n , k n ) does not always belong to the set L n . This is the case only under an additional condition; namely, only when the inclusion k n ∈ C −1,0 (0, +∞) holds true (of course, κ n ∈ C −1 (0, +∞) for any n ∈ N). This is a very strong and restrictive condition. For example, in the case of the Riemann-Liouville fractional integral I α 0+ with the kernel κ(t) = h α (t), 0 < α < 1 and the Riemann-Liouville fractional derivative D α 0+ with the kernel k(t) = h 1−α , the kernel k n takes the form k n (t) = h n(1−α) (t). It belongs to the space C −1,0 (0, +∞) only under the condition 0 < n(1 − α) < 1; i.e., if 1 − 1 n < α < 1, which is very restrictive. Moreover, the example of the kernels (43) shows that not every pair of the kernels from L n can be represented in the form (κ n , k n ) with the kernels (κ, k) ∈ L 1 .
Another and even more general and important possibility for the construction of a pair (κ n , k n ) of the kernels from L n , n > 1 based on the Sonine kernels κ, k from L 1 is presented in the following theorem: Theorem 3. Let (κ, k) be a pair of the Sonine kernels from L 1 .
Then, the pair (κ n , k n ) of the kernels given by the formula belongs to the set L n .
Proof. First, we check that the kernels (45) satisfy the condition (35): It is worth mentioning that the Caputo type GFD (49) can be represented in a slightly different form: As regards the basic properties of the GFI (47) of an arbitrary order on C −1 (0, +∞), they follow from the well-known properties of the Laplace convolution (compare these to the properties of the GFI (19) of the order less than one): To justify this denotation of GFIs and GFDs, in the rest of this section, we formulate and prove the two fundamental theorems of FC for the GFDs (48) and (49) of the Riemann-Liouville and Caputo types.
Theorem 4 (First Fundamental Theorem for the GFD of an Arbitrary Order). Let (κ, k) be a pair of the kernels from L n .

Proof.
We start with a proof of the Formula (61): .
A function f ∈ C −1,(k) (0, +∞) can be represented in the form f (t) = (I (k) φ)(t), φ ∈ C −1 (0, +∞), and thus the following chain of equations is valid: The last relation implicates the inclusion I (κ) f ∈ C n −1 (0, +∞) and the relations To derive Formula (62), we employ the representation (57) of the GFD of the Caputo type, Formula (63) and the relation (61) that we already proved: Theorem 5 (Second Fundamental Theorem for the GFD of an Abitrary Order). Let (κ, k) be a pair of the kernels from L n .
Then, the relation Sonine kernels from L 1 , the pair (κ n , k n ) of the kernels given by the Formula (45) belongs to the set L n , n > 1. The GFI (47) with the kernel κ n = ({1} n−1 * κ)(t) takes the form whereas the GFDs of the Riemann-Liouville and Caputo types with the kernel k n = k can be represented as follows: As we see, these constructions are completely analogical to the definitions of the Riemann-Liouville fractional integral and the Riemann-Liouville and Caputo fractional derivatives of an arbitrary order.

Conclusions
Starting from the work presented in [7], the so-called GFDs of the Riemann-Liouville and Caputo types have become a topic of active research in FC. In particular, both the ordinary and the partial fractional differential equations with these derivatives have been considered (see [10] for a survey of some recent results). However, the GFDs introduced to date have been based on the classical Sonine condition, and thus their "generalized order" was restricted to the interval (0, 1). In particular, the initial value problems for the fractional differential equations with these derivatives permitted only one initial condition, and thus no models for the intermediate processes between diffusion and wave propagation could be formulated in terms of these GFDs.
The main contribution of this paper is an extension of the definitions of the GFIs and GFDs to the case of arbitrary order. To achieve this aim, a suitable generalization of the Sonine condition was introduced, and some important classes of the kernels that satisfy this generalized condition were described. The kernels of the GFDs of an arbitrary order possess integrable singularities at the point zero. However, the kernels of the GFIs can be both singular (in the case of an order less than one) and continuous (in the case of an order greater or equal to one) at the origin. The conventional Riemann-Liouville and Caputo fractional derivatives of arbitrary order are particular cases of these GFDs. Another important example is the integro-differential operators of the convolution type with the Bessel and the modified Bessel functions in the kernels.
To justify the denotation of GFIs and GFDs of arbitrary order, in this paper, two fundamental theorems of fractional calculus for these operators were formulated and proved. The constructions introduced in this paper allow the formulation of the initialvalue problems for the fractional differential equations with GFDs of a generalized order greater than one with several initial conditions. Thus, further research regarding the properties of the GFIs and GFDs of an arbitrary order introduced in this paper as well as applications of the fractional differential equations with the GFDs of arbitrary order to model, for instance, the processes intermediate between diffusion and wave propagation is needed.