Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Robust Localization for Near- and Far-Field Signals with an Unknown Number of Sources
Next Article in Special Issue
A Rigorous Analysis of Integro-Differential Operators with Non-Singular Kernels
Previous Article in Journal
Existence Results for Caputo Tripled Fractional Differential Inclusions with Integral and Multi-Point Boundary Conditions
Previous Article in Special Issue
Stability and Controllability Study for Mixed Integral Fractional Delay Dynamic Systems Endowed with Impulsive Effects on Time Scales
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the 1st-Level General Fractional Derivatives of Arbitrary Order

Department of Mathematics, Physics and Chemistry, Berlin University of Applied Sciences and Technology, Luxemburger Str. 10, 13353 Berlin, Germany
Fractal Fract. 2023, 7(2), 183; https://doi.org/10.3390/fractalfract7020183
Submission received: 16 January 2023 / Revised: 2 February 2023 / Accepted: 10 February 2023 / Published: 12 February 2023
(This article belongs to the Special Issue New Trends on Generalized Fractional Calculus)

Abstract

:
In this paper, the 1st-level general fractional derivatives of arbitrary order are defined and investigated for the first time. We start with a generalization of the Sonin condition for the kernels of the general fractional integrals and derivatives and then specify a set of the kernels that satisfy this condition and possess an integrable singularity of the power law type at the origin. The 1st-level general fractional derivatives of arbitrary order are integro-differential operators of convolution type with the kernels from this set. They contain both the general fractional derivatives of arbitrary order of the Riemann–Liouville type and the regularized general fractional derivatives of arbitrary order considered in the literature so far. For the 1st-level general fractional derivatives of arbitrary order, some important properties, including the 1st and the 2nd fundamental theorems of fractional calculus, are formulated and proved.

1. Introduction

Within the last few years, a lot of attention in the fractional calculus (FC) literature has been devoted to the so-called generalized Riemann–Liouville fractional derivative, nowadays often referred to as the Hilfer fractional derivative that is defined as follows:
( D 0 + α , β f ) ( t ) : = ( I 0 + β ( n α ) d n d t n ( I 0 + ( 1 β ) ( n α ) f ) ) ( t ) , t > 0 , n 1 < α n , n N ,
where the Riemann–Liouville fractional integral I 0 + α of the order α , α > 0 is given by the formula
( I 0 + α f ) ( t ) : = τ α 1 Γ ( α ) f ( t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 f ( τ ) d τ , t > 0 .
Because the operator family I 0 + α tends to the identity operator, say, in the sense of the L p -norm as α 0 + , the Riemann–Liouville fractional integral of the order α = 0 is defined as the identity operator:
( I 0 + 0 f ) ( t ) : = f ( t ) , t > 0 .
In Ref. [1] (see also [2]), the operator (1) was treated in the case of the derivative order α ( 0 , 1 ] and type β [ 0 , 1 ] and in [3] in the case of an arbitrary non-negative order α ( n 1 < α n , n N ) and type β [ 0 , 1 ] .
The Hilfer fractional derivative with the order α ( 0 , 1 ] is a particular case of a more general operator introduced in the paper [4] by Djrbashian and Nersessian published in 1968. However, the operator (1) with the derivative order α > 1 seems to be not considered in the literature before publication of the paper [3].
It is worth mentioning that in [5], a different parametrization of the Djrbashian–Nersessian operators was introduced in the form
( D n L α , ( γ ) f ) ( t ) : = k = 1 n ( I γ k d d t ) ( I n α s n f ) ( t ) ,
where
s k = i = 1 k γ i , k = 1 , 2 , , n
and the parameters α ( 0 , 1 ] and γ 1 , γ 2 , , γ n R satisfy the conditions
0 γ k a n d α + s k k , k = 1 , 2 , , n .
Because the operators (4) are compositions of n first-order derivatives and n + 1 Riemann–Liouville fractional integrals, in [5], they were called the nth-level fractional derivatives of order α , 0 < α 1 and type γ = ( γ 1 , γ 2 , , γ n ) . It is easy to see that the Hilfer fractional derivative (1) with the order α ( 0 , 1 ] can be interpreted as a 1st level fractional derivative. In this paper, we introduce a generalization of the Hilfer fractional derivative of arbitrary order in the form of the integro-differential operators with some general kernels. By analogy with the 1st-level fractional derivative (the Hilfer fractional derivative), we call these operators the 1st-level general fractional derivatives.
The main advantage of the Hilfer fractional derivative (1) is that this operator contains both the Riemann–Liouville fractional derivative D 0 + α and the Caputo fractional derivative * D 0 + α as its particular cases ( n 1 < α n , n N ):
( D 0 + α , 0 f ) ( x ) = ( I 0 + 0 d n d t n I 0 + n α f ) ( t ) = d n d t n ( I 0 + n α f ) ( t ) = : ( D 0 + α f ) ( t ) ,
( D 0 + α , 1 f ) ( t ) = ( I 0 + n α d n d t n I 0 + 0 f ) ( t ) = ( I 0 + n α f ( n ) ) ( t ) = : ( * D 0 + α f ) ( t ) .
Moreover, the Hilfer fractional derivative (1) with any value of the parameter β [ 0 , 1 ] is a left-inverse operator to the Riemann–Liouville fractional integral (2) of the order α , α > 0 (see, for example, [5]) and thus, it can be interpreted as a fractional derivative of the order α . Because of the relations (5) and (6), any result derived for the Hilfer fractional derivative covers the analogous results for the Riemann–Liouville and for the Caputo fractional derivatives that are often obtained using the different methods and in the separate publications.
Another hot topic in the current FC literature is the so-called general fractional integrals and derivatives with the Sonin kernels that satisfy the Sonin condition ([6]):
( κ k ) ( t ) = 0 t κ ( t τ ) k ( τ ) d τ = { 1 } , t > 0 ,
where * denotes the Laplace convolution, and { 1 } stands for the function identically equal to one for t > 0 . For a certain kernel κ , the function k in the relation (7) is called a kernel associated to the kernel κ . A very important pair of the Sonin kernels is provided in terms of the power law functions
κ ( t ) = h α ( t ) , k ( t ) = h 1 α ( t ) , 0 < α < 1 ,
where
h α ( t ) : = t α 1 Γ ( α ) , α > 0 .
The relation (7) for the functions h α and h 1 α in the form
( h α h 1 α ) ( t ) = h 1 ( t ) = { 1 } , t > 0 , 0 < α < 1
was known already to Abel. In Ref. [7], Abel used this formula to derive a solution to the integro-differential equation (in slightly different notations)
f ( t ) = 1 Γ ( 1 α ) 0 t ( t τ ) α ϕ ( τ ) d τ , t > 0 , 0 < α < 1
with ϕ ( 0 ) = 0 in the form
ϕ ( t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 f ( τ ) d τ , t > 0 .
In Ref. [6], Sonin recognized that the same method can be applied to solve the integral equation
f ( t ) = ( κ ϕ ) ( t ) = 0 t κ ( t τ ) ϕ ( τ ) d τ , t > 0
with a kernel κ that satisfies the condition (7). Its solution can be represented in the form
ϕ ( t ) = d d t ( k f ) ( t ) = d d t 0 t k ( t τ ) f ( τ ) d τ , t > 0 ,
where k is the Sonin kernel associated to the kernel κ .
In Ref. [8] (see also [9,10,11]), Kochubei investigated the properties of the integro-differential operators of the convolution type with the kernels k that satisfy the conditions (K1)–(K4) below:
(K1)
The Laplace transform k ˜ of k
k ˜ ( p ) = ( L k ) ( p ) = 0 + k ( t ) e p t d t
exists for all p > 0 .
(K2)
k ˜ ( p ) is a Stieltjes function (see [12] for details regarding the Stieltjes functions).
(K3)
k ˜ ( p ) 0 and p k ˜ ( p ) + as p + .
(K4)
k ˜ ( p ) + and p k ˜ ( p ) 0 as p 0 .
It turns out that for any function k that satisfies the properties (K1)–(K4), there exists a function κ such that the Sonin condition (7) is fulfilled. Thus, the functions κ and k are the Sonin kernels. In what follows, the set of such kernels is denoted by K and referred to as the Kochubei set of the Sonin kernels.
For ( κ , k ) K , Kochubei interpreted the operators defined by the right-hand sides of Equations (13) and (14) as the general fractional integral (GFI) and the general fractional derivative (GFD), respectively:
( I ( κ ) f ) ( t ) : = ( κ f ) ( t ) = 0 t κ ( t τ ) f ( τ ) d τ , t > 0 ,
( D ( k ) f ) ( t ) : = d d t ( k f ) ( t ) = d d t ( I ( k ) f ) ( t ) , t > 0 .
Kochubei also introduced a regularized GFD in the form
( * D ( k ) f ) ( t ) : = ( D ( k ) f ) ( t ) f ( 0 ) k ( t ) , t > 0 .
For the functions with the integrable first order derivative, this GFD can be also represented as follows:
( * D ( k ) f ) ( t ) = ( I ( k ) f ) ( t ) , t > 0 .
The GFD (17) and the regularized GFD (18) with the kernels from the Kochubei set K are the left inverse operators to the GFI defined by (16) [5,8,13]. Moreover, the solutions to the Cauchy problems for the time-fractional ordinary and partial differential equations with the regularized GFD (18) possess some typical features of solutions to the evolution equations; see [8,13,14,15] for more details. For a treatment of some inverse problems for the fractional differential equations with the GFDs, we refer to the publications [16,17,18]. The Adams-type predictor–corrector method for the numerical solution of the fractional differential equations with the GFDs was presented in the recent paper [19].
It is worth mentioning that the properties of the GFI defined as in (16) and of the GFDs defined as in (17) and (18) essentially depend on the additional conditions posed of their kernels. In this sense, there exists not just one, but several theories of the GFIs and the GFDs, each one constructed for the operators with the kernels from its own set.
Recently, in a series of papers [20,21,22,23,24,25,26] by Luchko, another important set of the Sonin kernels was introduced. The kernels from this set are continuous on the real positive semi-axes and possess an integrable singularity of the power function type at the point zero. In Ref. [27] (see also [28,29,30,31,32]), this set of Sonin kernels was called the Luchko set, and the GFIs and the GFDs with these kernels were referred to as the Luchko GFIs and the Luchko GFDs, respectively. In the papers [20,21,22,23,24,25,26], the GFI (16), the GFD (17), and the regularized GFD (18) with the Sonin kernels from the Luchko set were studied on the space of functions that are continuous on the real positive semi-axis and have an integrable singularity of power function type at point zero and its suitable sub-spaces.
In the meantime, the GFI (16), the GFD (17), and the regularized GFD (18) with the Sonin kernels from the Luchko set were already employed in a number of innovative mathematical models of some important physical theories. In particular, in [28,29,30,31,32,33,34], Tarasov used these operators for the formulation of general fractional dynamics, general non-Markovian quantum dynamics, general fractional vector calculus, general non-local continuum mechanics, a non-local probability theory, non-local statistical mechanics, and a non-local gravity theory, respectively. It is also worth mentioning that the theory of the GFIs and the GFDs was recently employed in [35] for the construction of the so-called Scarpi variable-order fractional integrals and derivatives.
In Ref. [26], a construction of a GFD that comprises both the GFD (17) and the regularized GFD (18) was introduced for the first time in the FC literature. Following the notation suggested in [5] for the fractional derivative (1) with the order less than or equal to one, this GFD was called the 1st-level GFD. To define the 1st-level GFD, the Sonin condition (7) was extended to the case of three kernels. Then, a suitable set of kernels that satisfy this extended condition was introduced. For the kernels from this set, the GFI was defined as in (16), whereas the 1st-level GFD was introduced as a composition of a GFI, the first-order derivative, and another GFI. In the case of the power law kernels, this 1st-level GFD is reduced to the Hilfer fractional derivative (1) with the order less than or equal to one.
It is worth mentioning that the construction of the 1st-level GFD proposed in [26] is restricted to the case of the “generalized fractional order” less than or equal to one. In the case of the Hilfer fractional derivative, this corresponds to the value n = 1 in Equation (1). However, in [21], the GFD of arbitrary order and the regularized GFD of arbitrary order were defined and investigated. Moreover, the case of the general FC operators of arbitrary order is very important for applications of this theory, say, for the construction of the general non-local models of the anomalous diffusion-wave processes.
The main aim of this paper is defining and investigating the 1st-level GFD of arbitrary order that comprises both the GFD and the regularized GFD of arbitrary order that were introduced in [21].
The rest of the paper is organized as follows: Section 2 is devoted to a presentation of the basic results derived in the literature for the GFD of arbitrary order and the regularized GFD of arbitrary order with the kernels from the Luchko set. In Section 3, the 1st-level GFD of arbitrary order is introduced and investigated. For this fractional derivative and the GFI defined as in (16), the 1st and the 2nd fundamental theorems of FC are formulated and proved. As a consequence of the 2nd fundamental theorem of FC, a formula for the projector operator for the 1st-level GFD of arbitrary order is derived. In particular, this formula determines the form of the natural initial conditions for the fractional differential equations with the 1st-level GFDs of arbitrary order. Section 4 contains some discussions and directions for further research.

2. The GFI, the GFD, and the Regularized GFD of Arbitrary Order

In this section, we provide the definitions and the basic properties of the GFI of arbitrary order, the GFD of arbitrary order, and the regularized GFD of arbitrary order that were introduced in [21] for the first time.
The GFI, the GFD, and the regularized GFD are a far reaching generalization of the Riemann–Liouville fractional integral (2), the Riemann–Liouville fractional derivative (5), and the Caputo fractional derivative (6), respectively. However, the GFI (16) and the GFDs (17) and (18) with the Sonin kernels cover only the case of the Riemann–Liouville fractional integral I 0 + α , the Riemann–Liouville fractional derivative D 0 + α , and the Caputo fractional derivative * D 0 + α with the order α from the interval ( 0 , 1 ) :
( I ( h α ) f ) ( t ) = ( h α f ) ( t ) = ( I 0 + α f ) ( t ) , 0 < α < 1 , t > 0 ,
( D ( h 1 α ) f ) ( t ) = d d t ( h 1 α f ) ( t ) = d d t ( I 0 + 1 α f ) ( t ) = ( D 0 + α f ) ( t ) , 0 < α < 1 , t > 0 ,
( * D ( h 1 α ) f ) ( t ) = ( h 1 α f ) ( t ) = ( I 0 + 1 α f ) ( t ) = ( * D 0 + α f ) ( t ) , 0 < α < 1 , t > 0 .
The condition 0 < α < 1 in Equations (20)–(22) is a consequence of the fact that the power law functions κ ( t ) = h α ( t ) and k ( t ) = h 1 α ( t ) are the Sonin kernels if and only if the parameter α satisfies this condition.
It is worth mentioning that the Sonin condition (10) can be extended to the case α = 0 or α = 1 , respectively, in the sense of the generalized functions (the function h 0 plays the role of the Dirac δ -function):
h 0 h 1 = h 1 .
The relation (23) immediately leads to the following interpretations of the GFI (16), the GFD (17), and the regularized GFD (18) (or (19)) with the kernels h 0 and h 1 :
( I ( h 0 ) f ) ( t ) = ( h 0 f ) ( t ) = ( I 0 + 0 f ) ( t ) = f ( t ) , t > 0 ,
( I ( h 1 ) f ) ( t ) = ( h 1 f ) ( t ) = ( I 0 + 1 f ) ( t ) = 0 t f ( t ) d t , t > 0 ,
( D ( h 0 ) f ) ( t ) = d d t ( h 0 f ) ( t ) = f ( t ) , t > 0 ,
( D ( h 1 ) f ) ( t ) = d d t ( h 1 f ) ( t ) = d d t 0 t f ( t ) d t = f ( t ) , t > 0 ,
( * D ( h 0 ) f ) ( t ) = ( I ( h 0 ) f ) ( t ) = f ( t ) , t > 0 ,
( * D ( h 1 ) f ) ( t ) = ( I ( h 1 ) f ) ( t ) = 0 t f ( t ) d t = f ( t ) f ( 0 ) , t > 0 .
Because the right-hand side of the Sonin condition (7) is the function h 1 that corresponds to the definite integral (the Riemann–Liouville fractional integral of the order one), the “generalized fractional orders” of the GFI (16), the GFD (17), and the regularized GFD (18) with the Sonin kernels κ and k are less than or equal to one.
However, it is well known that both the Riemann–Liouville fractional integral and the Riemann–Liouville and Caputo fractional derivatives are well defined for an arbitrary order α 0 . To introduce the GFI, the GFD, and the regularized GFD of arbitrary order, in [21], an extension of the Sonin condition (7) for the kernels κ and k has been suggested, and a set of the kernels that satisfy the extended Sonin condition and belong to the suitable spaces of functions has been specified. In this section, we represent those results from [21] that will be needed for the further discussions.
Definition 1 
([21]). Let the functions κ and k defined on the positive real semi-axes satisfy the following conditions:
(L1)
( κ k ) ( t ) = { 1 } < n > ( t ) , n N , t > 0 ,
where
{ 1 } < n > ( t ) : = ( { 1 } { 1 } n t i m e s ) ( t ) = h n ( t ) = t n 1 ( n 1 ) ! .
(L2)
κ C 1 ( 0 , + ) ,
where
C 1 ( 0 , + ) : = { f : f ( t ) = t p f 1 ( t ) , t > 0 , p > 1 , f 1 C [ 0 , + ) } .
(L3)
k C 1 , 0 ( 0 , + ) ,
where
C 1 , 0 ( 0 , + ) : = { f : f ( t ) = t p f 1 ( t ) , t > 0 , 1 < p < 0 , f 1 C [ 0 , + ) } .
The set of the ordered pairs ( κ , k ) of such kernels is denoted by L n .
Remark 1. 
The condition (30) is a generalization of the Sonin condition (7). The kernels ( κ , k ) that satisfy the condition (30) are not the Sonin kernels unless n = 1 . Moreover, in Definition 1, the spaces of functions for the kernels κ and k are specified in the inclusions (31) and (33). Both the condition (30) and the inclusions (31) and (33) are essential components for the construction of a self-contained theory of the GFIs and the GFDs with the kernels from the set L n .
Remark 2. 
Evidently, the inclusion C 1 , 0 ( 0 , + ) C 1 ( 0 , + ) holds true. Any function from C 1 , 0 ( 0 , + ) is continuous on the real positive semi-axes and has an integrable singularity at the point zero. In contrast, the functions from C 1 ( 0 , + ) can have an integrable singularity at the point zero or not. Thus, in general, the kernels κ and k from Definition 1 cannot be interchanged if n 2 . However, in the case n = 1 , both functions κ and k are the Sonin kernels that have an integrable singularity at the point zero and thus κ , k C 1 , 0 ( 0 , + ) if ( κ , k ) L 1 .
Remark 3. 
The kernels of all reasonable time-fractional integrals and derivatives introduced so far in the FC literature belong to the set L n . In particular, the kernels κ ( t ) = h α ( t ) and k ( t ) = h n α ( t ) of the Riemann–Liouville fractional integral of arbitrary order α > 0 defined by Equation (2) and the Riemann–Liouville and Caputo fractional derivatives of arbitrary order α > 0 defined as in Equations (5) and (6), respectively, are from the set L n if the order α satisfies the conditions n 1 < α < n , n N .
Indeed, the inclusions h α C 1 ( 0 , + ) and h n α C 1 , 0 ( 0 , + ) evidently hold true for any α under the conditions n 1 < α < n , n N . The relation
( h α h n α ) ( t ) = h n ( t ) , n 1 < α < n , n N
is a particular case of the formula
( h α h β ) ( t ) = h α + β ( t ) , α , β > 0 ,
that immediately follows from the well-known relation between the Euler beta- and gamma functions:
B ( α , β ) = 0 1 t α 1 ( 1 t ) β 1 d t = Γ ( α ) Γ ( β ) Γ ( α + β ) .
For n = 1 , the condition (30) is just the Sonin condition (7). Thus, the set L 1 contains all Sonin kernels that belong to the space C 1 , 0 ( 0 , + ) (see Remark 2). Because the space C 1 , 0 ( 0 , + ) is a very general one, almost all of the Sonin kernels introduced so far belong to this space. In particular, in [6], the following family of the Sonin kernels was derived:
κ ( t ) = h α ( t ) · κ 1 ( t ) , κ 1 ( t ) = k = 0 + a k t k , a 0 0 , 0 < α < 1 ,
k ( t ) = h 1 α ( t ) · k 1 ( t ) , k 1 ( t ) = k = 0 + b k t k ,
where κ 1 and k 1 are analytical functions, and the coefficients a k , b k , k N 0 satisfy an infinite triangular system of linear equations:
a 0 b 0 = 1 , k = 0 n Γ ( k + 1 α ) Γ ( α + n k ) a n k b k = 0 , n 1 .
Evidently, the kernels κ and k specified by Equations (37) and (38) are from the set L 1 . In Ref. [6], Sonin derived an important particular case of the kernels in forms (37) and (38):
κ ( t ) = t ν / 2 J ν ( 2 t ) , k ( t ) = t ν / 2 1 / 2 I ν 1 ( 2 t ) , 1 < ν < 0 ,
where
J ν ( t ) = k = 0 + ( 1 ) k ( t / 2 ) 2 k + ν k ! Γ ( k + ν + 1 ) a n d I ν ( t ) = k = 0 + ( t / 2 ) 2 k + ν k ! Γ ( k + ν + 1 )
are the Bessel and the modified Bessel functions, respectively. For other examples of the Sonin kernels that belong to the set L 1 , we refer to [6,8,20,22,36,37,38].
In Ref. [21], the set L n of the Luchko kernels was defined for the first time. Many examples of the kernels from L n as well as several methods for their construction based on the kernels from L 1 were presented in [21,27], and we shortly mention here some of them.
The first method employs the relation (30) in the Laplace domain (provided that the Laplace transforms κ ˜ , k ˜ of the kernels κ and k do exist) in the form
κ ˜ ( p ) · k ˜ ( p ) = 1 p n , ( p ) > p κ , k R , n N ,
along with the tables of the direct and inverse Laplace transforms (see, for example, [39,40]). In particular, in [21], this method was used to derive the following kernels from the set L n :
κ ( t ) = t ν / 2 J ν ( 2 t ) , k ( t ) = t n / 2 ν / 2 1 I n ν 2 ( 2 t ) , n 2 < ν < n 1 , n N ,
where the Bessel function J ν and the modified Bessel function I ν are defined as in (41). Please note that for n = 1 , the kernels from (43) are reduced to the Sonin kernels (40).
Let now κ and k be the kernels from the set L 1 , i.e., the kernels that satisfy the Sonin condition (7). Then the kernels κ ^ = κ < n > and k ^ = k < n > satisfy the relation (30)
( κ ^ k ^ ) ( t ) = ( κ < n > k < n > ) ( t ) = ( κ k ) < n > ( t ) = { 1 } < n > ( t ) .
Because of the inclusion κ C 1 ( 0 , + ) , the kernel κ ^ = κ < n > belongs to the space C 1 ( 0 , + ) [41]. However, the inclusion k ^ = k < n > C 1 , 0 ( 0 , + ) does not always hold true. Say, in the case of the Riemann–Liouville fractional derivative with the kernel k ( t ) = h 1 α , the kernel k < n > ( t ) = h n ( 1 α ) ( t ) belongs to the space C 1 , 0 ( 0 , + ) if and only if 1 1 n < α < 1 . Thus, the pair of the kernels κ ^ = κ < n > and k ^ = k < n > is from the set L n if and only if k < n > C 1 , 0 ( 0 , + ) . On the other hand, if k < n > C 1 , 0 ( 0 , + ) and if one of the derivatives d j k < n > d t j , j = 1 , , n 1 belongs to the space C 1 , 0 ( 0 , + ) , we obtain the relation
( κ ^ d j d t j k < n > ) ( t ) = d j d t j ( κ < n > k < n > ) ( t ) = d j d t j { 1 } < n > ( t ) = { 1 } < n j > ( t ) .
Thus, the pair of the kernels κ ^ = κ < n > and d j d t j k < n > belongs to the set L n j (it is easy to see that there exists at most one derivative of the order j { 1 , , n 1 } that satisfies the inclusion d j d t j k < n > C 1 , 0 ( 0 , + ) ).
To illustrate the procedure described above, let us consider two examples:
(1)
κ ( t ) = h 2 / 3 ( t ) , k ( t ) = h 1 / 3 ( t ) . Evidently, κ and k are the kernels from the set L 1 . Then κ < 2 > ( t ) = ( κ κ ) ( t ) = h 4 / 3 ( t ) and k < 2 > ( t ) = ( k * k ) ( t ) = h 2 / 3 ( t ) (see Equation (36)). Because of the inclusion k < 2 > ( t ) = h 2 / 3 ( t ) C 1 , 0 ( 0 , + ) , the kernels κ < 2 > ( t ) = h 4 / 3 ( t ) and k < 2 > ( t ) = h 2 / 3 ( t ) are from the set L 2 .
(2)
κ ( t ) = h 1 / 3 ( t ) , k ( t ) = h 2 / 3 ( t ) . Once again, κ and k are the kernels from the set L 1 . Then κ < 2 > ( t ) = ( κ κ ) ( t ) = h 2 / 3 ( t ) and k < 2 > ( t ) = ( k k ) ( t ) = h 4 / 3 ( t ) . This time, k < 2 > ( t ) = h 4 / 3 ( t ) C 1 , 0 ( 0 , + ) . However, the function d d t k < 2 > ( t ) = d d t h 4 / 3 ( t ) = h 1 / 3 ( t ) is from the set C 1 , 0 ( 0 , + ) . Thus, the kernels κ < 2 > ( t ) = h 2 / 3 ( t ) and d d t k < 2 > ( t ) = h 1 / 3 ( t ) are from the set L 1 .
In Ref. [27], a generalization of the method presented above is suggested. Let the pairs of the kernels ( κ 1 , k 1 ) ,..., ( κ n , k n ) belong to the set L 1 . Then the kernels κ ( t ) = ( κ 1 κ n ) ( t ) and k ( t ) = ( k 1 k n ) ( t ) satisfy the condition (30). The inclusion κ C 1 ( 0 , + ) is also fulfilled. Thus, the pair of the kernels ( κ , k ) is from the set L n if and only if k C 1 , 0 ( 0 , + ) . If k C 1 , 0 ( 0 , + ) and there exists j { 1 , , n 1 } such that d j d t j k C 1 , 0 ( 0 , + ) , the pair of the kernels ( κ , d j d t j k ) is from the set L n j (see the discussions above).
Another important method for the construction of kernels from the set L n , n > 1 based on the Sonin kernels ( κ , k ) L 1 is as follows: ([21]):
κ ^ ( t ) = ( { 1 } < n 1 > κ ) ( t ) , k ^ ( t ) = k ( t ) .
Evidently, the kernels κ ^ , k ^ defined as in (44) satisfy all three conditions from Definition 1 and thus ( κ ^ , k ^ ) L n .
In the rest of this section, for the kernels ( κ , k ) from the set L n , we define the GFI of arbitrary order with the kernel κ , the GFD of arbitrary order with the kernel k, and the regularized GFD of arbitrary order with the kernel k and discuss some of their basic properties.
Definition 2 
([21]). Let ( κ , k ) be a pair of the kernels from the set L n . The GFI of arbitrary order with the kernel κ, the GFD of arbitrary order with the kernel k, and the regularized GFD of arbitrary order with the kernel k are defined by the following formulas, respectively:
( I ( κ ) f ) ( t ) : = ( κ f ) ( t ) = 0 t κ ( t τ ) f ( τ ) d τ , t > 0 ,
( D ( k ) f ) ( t ) : = d n d t n ( k f ) ( t ) = d n d t n ( I ( k ) f ) ( t ) , t > 0 ,
( * D ( k ) f ) ( t ) : = D ( k ) f ( · ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( · ) ( t ) , t > 0 .
It is worth mentioning that under certain conditions, the regularized GFD (47) can be represented in a different form [21]:
( * D ( k ) f ) ( t ) = ( D ( k ) f ) ( t ) j = 0 n 1 f ( j ) ( 0 ) d n j 1 d t n j 1 k ( t ) , t > 0 .
In particular, the representation (48) is valid for f , k C 1 n 1 ( 0 , + ) , where the space of functions C 1 m ( 0 , + ) , m = 0 , 1 , 2 , is defined as follows:
C 1 m ( 0 , + ) : = { f : f ( m ) C 1 ( 0 , + ) } .
For m = 0 , the space C 1 m ( 0 , + ) is interpreted as the space C 1 ( 0 , + ) . This means, in particular, that for n = 1 , the representation (48) is valid without any additional conditions on the kernel k.
Moreover, for f C 1 n ( 0 , + ) , the regularized GFD (47) can be represented as follows [21]:
( * D ( k ) f ) ( t ) = ( I ( k ) f ( n ) ) ( t ) = 0 t k ( t τ ) f ( n ) ( τ ) d τ , t > 0 .
As mentioned in Remark 3, the pair of the kernels κ ( t ) = h α ( t ) , α > 0 and k ( t ) = h n α ( t ) belongs to the set L n provided that n 1 < α < n , n N . This implicates that the GFI (45) with the kernel κ ( t ) = h α ( t ) , α > 0 is the Riemann–Liouville fractional integral of the order α > 0 , whereas the Riemann=-Liouville and the Caputo fractional derivatives of the order α , n 1 < α < n , n N are particular cases of the GFD (46) and the regularized GFD (47) (or (48) or (50)) with the kernel k ( t ) = h n α ( t ) , respectively.
The constructions of the GFD and the regularized GFD of arbitrary order presented in Definition 2 produce the Riemann–Liouville and the Caputo fractional derivatives of the non-integer order α that satisfy the conditions n 1 < α < n , n N . To define them for α = n N 0 , we proceed as in the case n = 1 (see the discussions at the beginning of this section). Namely, we extend the relation (35) to the case α = 0 or α = n N , respectively, in the sense of the generalized functions (the function h 0 plays the role of the Dirac δ -function):
h 0 h n = h n , n N .
Then the GFI (45), the GFD (46), and the regularized GFD (47) with the kernels h 0 and h n can be interpreted as follows:
( I ( h 0 ) f ) ( t ) = ( h 0 f ) ( t ) = f ( t ) , t > 0 ,
( I ( h n ) f ) ( t ) = ( h n f ) ( t ) = ( I 0 + n f ) ( t ) , t > 0 ,
( D ( h 0 ) f ) ( t ) = d n d t n ( h 0 f ) ( t ) = f ( n ) ( t ) , t > 0 ,
( D ( h n ) f ) ( t ) = d n d t n ( h n f ) ( t ) = d n d t n ( I 0 + n f ) ( t ) = f ( t ) , t > 0 ,
( * D ( h 0 ) f ) ( t ) = D ( h 0 ) f ( · ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( · ) ( t ) = f ( n ) ( t ) , t > 0 ,
( * D ( h n ) f ) ( t ) = D ( h n ) f ( · ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( · ) ( t ) = f ( t ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( t ) , t > 0 .
As an example of a new and non-trivial particular case of the GFI (45) of arbitrary order, the GFD (46) of arbitrary order, and the regularized GFD (47) of arbitrary order, respectively, we mention the following operators constructed for the kernels ( κ , k ) L n given by the formula (43) that is valid under the conditions n 2 < ν < n 1 , n N [21]:
( I ( κ ) f ) ( t ) = 0 t ( t τ ) ν / 2 J ν ( 2 t τ ) f ( τ ) d τ , t > 0 ,
( D ( k ) f ) ( t ) = d n d t n 0 t ( t τ ) n / 2 ν / 2 1 I n ν 2 ( 2 t τ ) f ( τ ) d τ , t > 0 ,
( * D ( k ) f ) ( t ) : = D ( k ) f ( · ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( · ) ( t ) , t > 0 .
Because the kernels ( κ , k ) L n of the GFI (45) of arbitrary order, the GFD (46) of arbitrary order, and the regularized GFD (47) of arbitrary order are from the space C 1 ( 0 , + ) , the natural domains of definition for these operators are the space C 1 ( 0 , + ) and its suitable sub-spaces. The basic properties of the GFI (45) easily follow from the properties of the Laplace convolution on the space C 1 ( 0 , + ) (see [41]):
I ( κ ) : C 1 ( 0 , + ) C 1 ( 0 , + ) , κ L n ( mapping property ) ,
I ( κ 1 ) I ( κ 2 ) = I ( κ 2 ) I ( κ 1 ) , κ 1 , κ 2 L n ( commutativity law ) ,
I ( κ 1 ) I ( κ 2 ) = I ( κ 1 * κ 2 ) , κ 1 , κ 2 L n ( semi - group property ) .
According to the axioms of FC suggested in [42], the GFDs should be left-inverse operators to the corresponding GFIs (1st fundamental theorem of FC). Moreover, as shown in [5,23,24], the compositions of the GFIs and the corresponding GFDs applied to a function f are the so-called convolution polynomials subtracted from the function f (2nd fundamental theorem of FC).
In the rest of this section, we present the fundamental theorems of FC for the GFI (45) of arbitrary order, the GFD (46) of arbitrary order, and for the regularized GFD (47) of arbitrary order. For the proofs of these theorems, we refer to [21,25].
Theorem 1 
([21]). Let ( κ , k ) be a pair of the kernels from the set L n .
Then the GFD (46) is a left-inverse operator to the GFI (45) on the space C 1 ( 0 , + ) :
( D ( k ) I ( κ ) f ) ( t ) = f ( t ) , f C 1 ( 0 , + ) , t > 0 ,
and the regularized GFD (47) is a left-inverse operator to the GFI (45) on the space I ( k ) ( C 1 ( 0 , + ) ) :
( * D ( k ) I ( κ ) f ) ( t ) = f ( t ) , f I ( k ) ( C 1 ( 0 , + ) ) , t > 0 ,
where
I ( k ) ( C 1 ( 0 , + ) ) : = { f : f ( t ) = ( I ( k ) ϕ ) ( t ) , ϕ C 1 ( 0 , + ) } .
Theorem 2 
([21]). Let ( κ , k ) be a pair of the kernels from the set L n .
Then the relation
( I ( κ ) * D ( k ) f ) ( t ) = f ( t ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( t )
holds true on the space C 1 n ( 0 , + ) defined as in (49), and the formula
( I ( κ ) D ( k ) f ) ( t ) = f ( t ) , t > 0
is valid for the functions f I ( κ ) ( C 1 ( 0 , + ) ) .
The Equation (68) means that the GFD of arbitrary order is a right-inverse operator to the corresponding GFI on the space I ( κ ) ( C 1 ( 0 , + ) ) . However, in general, this formula does not hold true if we consider the GFD (46) of arbitrary order on its natural domain of definition, namely, on the following space of functions:
C 1 , ( k ) 1 ( 0 , + ) = { f C 1 ( 0 , + ) : D ( k ) f C 1 ( 0 , + ) } .
For the functions from C 1 , ( k ) 1 ( 0 , + ) , a more general result compared to the one given by Equation (68) is valid. This result is provided in the theorem below.
Theorem 3 
([25]). Let ( κ , k ) be a pair of the kernels from the set L n and κ C 1 n 1 ( 0 , + ) .
Then for any function f C 1 , ( k ) 1 ( 0 , + ) , the formula
( I ( κ ) D ( k ) f ) ( t ) = f ( t ) j = 0 n 1 d j d t j I ( k ) f ( 0 ) d n j 1 d t n j 1 κ ( t ) , t > 0
holds true.
The difference in Equations (68) and (70) is caused by the inclusion I ( κ ) ( C 1 ( 0 , + ) ) C 1 , ( k ) 1 ( 0 , + ) that immediately follows from Equation (64) of Theorem 1. Moreover, for any function from the space I ( κ ) ( C 1 ( 0 , + ) ) , a comparison of Equations (68) and (70) leads to the relations
d j d t j I ( k ) f ( 0 ) = 0 , j = 0 , 1 , , n 1 .
It is worth mentioning that the left-hand sides of the relations (71) can be interpreted as the natural initial conditions while dealing with the fractional differential equations with the GFDs of arbitrary order (see [25]). In particular, in the case of the Riemann–Liouville fractional integral I 0 + α with the kernel κ ( t ) = h α ( t ) , n 1 < α < n , n N and the Riemann–Liouville fractional derivative D 0 + α with the kernel k ( t ) = h n α ( t ) , Equation (70) takes the following well-known form [43]:
( I 0 + α D 0 + α f ) ( t ) = f ( t ) j = 0 n 1 d j d t j I 0 + n α f ( 0 ) h α n + j + 1 ( t ) , t > 0 .
Thus, the natural initial conditions for the fractional differential equations with the Riemann–Liouville fractional derivatives are formulated as follows ([43]):
d j d t j I 0 + n α f ( 0 ) = a j , j = 0 , 1 , , n 1 .

3. The 1st-Level GFD of Arbitrary Order

Similar to the case of the Riemann–Liouville fractional derivative and the Caputo fractional derivative, both the definitions, the particular cases, and the formulas for the GFD (46) and for the regularized GFD (47) look very different. Moreover, as a rule, these GFDs as well as the fractional differential equations with these derivatives are treated in the separate publications and by employing the unequal methods. As mentioned in the Introduction, the Hilfer fractional derivative (1) is a unification and a generalization of the Riemann–Liouville fractional derivative and the Caputo fractional derivative. In this section, we introduce the 1st-level GFD of arbitrary order that generalizes the Hilfer fractional derivative (1) to the case of the arbitrary kernels that satisfy an extended Sonin condition and investigate its basic properties. The main particular cases of this derivative are the GFD (46) of arbitrary order and the regularized GFD (47) of arbitrary order. The case of the 1st-level GFD with the generalized order less than or equal to one was treated in [26].
We start with a suitable generalization of the kernels from Definition 1 that we call the kernels of the 1st-level GFDs of arbitrary order.
Definition 3. 
Let the functions κ , k 1 , k 2 : ( 0 , + ) R satisfy the following conditions:
1st condition:
( κ k 1 k 2 ) ( t ) = { 1 } < n > ( t ) = h n ( t ) = t n 1 ( n 1 ) ! , n N , t > 0 ,
2nd condition:
κ C 1 ( 0 , + ) ,
3rd condition:
k 1 , k 2 C 1 , 0 ( 0 , + ) .
The set of such triples ( κ , k 1 , k 2 ) is denoted by L n 1 .
Remark 4. 
As suggested in [26], the denotation L n m stands for the set of the kernels of the mth-level GFDs with the generalized order from the interval ( n 1 , n ) . The GFD (46) of arbitrary order and the regularized GFD (47) of arbitrary order with the kernels from L n treated in Section 2 can be interpreted as the 0th-level GFDs, and thus we set L n 0 : = L n . The case of the kernels from the set L 1 1 was considered in[26].
In what follows, we associate the first kernel from the triple ( κ , k 1 , k 2 ) L n 1 with the GFI defined as in (16), whereas the kernels k 1 and k 2 are assigned to the 1st-level GFD of arbitrary order that will be defined below. Due to the inclusion k 1 , k 2 C 1 , 0 ( 0 , + ) , the kernels k 1 and k 2 can be interchanged. However, for n 2 , the kernels k 1 and k 2 cannot be interchanged with the kernel κ and thus, the triples ( κ , k 1 , k 2 ) are partially ordered. It is worth mentioning that as soon as any two kernels from a triple ( κ , k 1 , k 2 ) are fixed, the 3rd kernel is uniquely determined by the relation (73). This property is a consequence from Theorem 3.1 in [41] that states that the ring R 1 = ( C 1 ( 0 , + ) , + , * ) does not have any divisors of zero. For the kernels k 1 and k 2 , the kernel κ from the triple ( κ , k 1 , k 2 ) is called the 1st-level GFD kernel associated to the kernel pair ( k 1 , k 2 ) .
As an example of the kernels from the set L n 1 we mention the following power law functions
κ ( t ) = h α ( t ) , k 1 ( t ) = h γ ( t ) , k 2 ( t ) = h n α γ ( t ) , t > 0
that satisfy the conditions from Definition 3 if
n 1 < α < n , 0 < γ < n α , n N .
Indeed, the function κ ( t ) = h α ( t ) is from the space C 1 ( 0 , + ) , whereas the functions k 1 ( t ) = h γ ( t ) and k 2 ( t ) = h n α γ ( t ) belong to the space C 1 , 0 ( 0 , + ) due to the inequalities (76). Equation (36) easily leads to the relation
( h α h γ h n α γ ) ( t ) = h α + γ + n α γ ( t ) = h n ( t ) , t > 0 .
Thus, we have the inclusion
( h α ( t ) , h γ ( t ) , h n α γ ( t ) ) L n 1 .
It is worth mentioning that the relation (77) can be also verified by employing the Laplace transform technique. Provided the Laplace transforms of the kernels κ , k 1 , and k 2 do exist, application of the Laplace transform to the condition (73) from Definition 3 leads to the relation
κ ˜ ( p ) · k 1 ˜ ( p ) · k 2 ˜ ( p ) = 1 p n .
In the case of the power law kernels (78), we immediately get the Laplace transform formulas
κ ˜ ( p ) = 1 p α , k 1 ˜ ( p ) = 1 p γ , k 2 ˜ ( p ) = 1 p n α γ
and thus the relation (79) evidently holds true:
1 p α · 1 p γ · 1 p n α γ = 1 p n .
The relation (79) and the tables of the Laplace transforms and the inverse Laplace transforms can be used to determine other triples of the kernels from L n 1 . For other techniques for construction of the Sonin kernels that can be employed for derivation of the 1st-level GFDs kernels, we refer the readers to [37].
Remark 5. 
In Definition 3, the kernels κ, k 1 , and k 2 are the functions from the spaces C 1 ( 0 , + ) and C 1 , 0 ( 0 , + ) , respectively. However, Equation (23) in the sense of generalized functions implicates the following forms of the condition (73) that involve the generalized function h 0 :
κ h 0 k 2 = κ k 2 = h n ,
κ k 1 h 0 = κ k 1 = h n .
Equations (80) and (81) can be interpreted as follows: if one of the kernels k 1 or k 2 of the 1st-level GFD of arbitrary order is set to be the generalized function h 0 , then the inclusion ( κ , k 2 ) L n or ( κ , k 1 ) L n , respectively, holds true, i.e., the pair of the remainder kernels satisfies the conditions from Definition 1.
Let us now proceed with defining the 1st-level GFDs of arbitrary order.
Definition 4. 
Let ( κ , k 1 , k 2 ) be a triple of the kernels from the set L n 1 .
The 1st-level GFD of arbitrary order is defined by the formula
( 1 L D ( k 1 , k 2 ) f ) ( t ) : = I ( k 1 ) d n d t n I ( k 2 ) f ( t ) ,
whereas the corresponding GFI with the kernel κ is provided by the relation (16).
An important particular case of the 1st-level GFD of arbitrary order is the Hilfer fractional derivative (1). Indeed, let us put the power law kernels k 1 and k 2 from the triple (78) into Equation (82). Because the GFI with a power law kernel is nothing else than the Riemann–Liouville fractional integral (2), Equation (82) takes the form
( 1 L D ( h γ , h n α γ ) f ) ( t ) = I 0 + γ d n d t n I 0 + n α γ f ( t ) .
It is easy to see that the operator at the right-hand side of (83) is the Hilfer fractional derivative (1) in a slightly different parametrization suggested in [5]. Indeed, the Hilfer derivative (1) coincides with the operator at the right-hand side of (83) if we set γ = β ( n α ) . The form (83) of the Hilfer fractional derivative with n = 1 was called in [5] the 1st-level fractional derivative. As we see, the operator (82) is a natural generalization of the Hilfer derivative (83) to the case of arbitrary kernels from the set L n 1 .
As in the case of the Hilfer fractional derivative that unifies the Riemann–Liouville and the Caputo fractional derivatives in one formula, the main idea behind the 1st level GFD (82) of arbitrary order is that this derivative contains both the GFD (46) of arbitrary order and the regularized GFD of arbitrary order in form (50) as its particular cases.
Indeed, the interpretation (52) of the GFI with the kernel h 0 and the discussions presented in Remark 5 immediately lead to the following particular cases of the 1st-level GFD:
  • For k 1 = h 0 , the 1st level GFD (82) of arbitrary order is reduced to the GFD (46) of arbitrary order with the kernel k 2 :
    ( 1 L D ( h 0 , k 2 ) f ) ( t ) = I ( h 0 ) d n d t n I ( k 2 ) f ( t ) = d n d t n I ( k 2 ) f ( t ) = ( D ( k 2 ) f ) ( t ) .
  • For k 2 = h 0 , the 1st level GFD (82) of arbitrary order is reduced to the regularized GFD (50) of arbitrary order with the kernel k 1 :
    ( 1 L D ( k 1 , h 0 ) f ) ( t ) = I ( k 1 ) d n d t n I ( h 0 ) f ( t ) = I ( k 1 ) d n d t n f ( t ) = ( * D ( k 1 ) f ) ( t ) .
In the publications devoted to the GFDs and the regularized GFDs, their properties were investigated separately and often by employing unequal methods for the GFDs and the regularized GFDs, respectively. The same statement is valid for the publications devoted to the fractional differential equations with the GFDs and the regularized GFDs. The concept of the 1st-level GFD introduced in Definition 4 opens a gateway for the derivation of formulas that are valid both for the GFDs and for the regularized GFDs. In the rest of this section, we formulate and prove some results for the 1st-level GFDs of arbitrary order, including the 1st and the 2nd fundamental theorems of FC. In particular, these results are valid both for the GFDs and for the regularized GFDs of arbitrary order.
We start with the 1st fundamental theorem of FC for the 1st-level GFDs of arbitrary order.
Theorem 4. 
Let a triple of the kernels ( κ , k 1 , k 2 ) belong to the set L n 1 .
Then the 1st-level GFD (82) of arbitrary order is a left inverse operator to the GFI (16) on the space I ( k 1 ) ( C 1 ( 0 , + ) ) defined as in (66):
( 1 L D ( k 1 , k 2 ) I ( κ ) f ) ( t ) = f ( t ) , f I ( k 1 ) ( C 1 ( 0 , + ) ) , t > 0 .
Proof. 
By definition, any function f from the space I ( k 1 ) ( C 1 ( 0 , + ) ) can be represented in the form
f ( t ) = ( I ( k 1 ) ϕ ) ( t ) = ( k 1 * ϕ ) ( t ) , ϕ C 1 ( 0 , + ) .
For the kernels ( κ , k 1 , k 2 ) L 1 1 , the condition (73) is fulfilled. Then we obtain the following chain of the equations
( 1 L D ( k 1 , k 2 ) I ( κ ) f ) ( t ) = I ( k 1 ) d n d t n I ( k 2 ) I ( κ ) f ( t ) = I ( k 1 ) d n d t n I ( k 2 ) I ( κ ) I ( k 1 ) ϕ ( t ) =
I ( k 1 ) d n d t n ( k 2 κ k 1 ϕ ) ( t ) ( t ) = I ( k 1 ) d n d t n ( { 1 } < n > ϕ ) ( t ) ( t ) = ( I ( k 1 ) ϕ ) ( t ) = f ( t )
which proves Equation (86). □
According to Theorem 1, the formula of type (86) is valid both for the GFD (46) of arbitrary order and for the regularized GFD (47) of arbitrary order. Both these formulas are particular cases of Equation (86). Indeed, setting k 1 = h 0 in (86) leads to Equation (64) for the GFD (46) of arbitrary order with the kernel k = k 2 on the space C 1 ( 0 , + ) , whereas a substitution k 2 = h 0 in (86) results in Equation (65) for the regularized GFD (47) of arbitrary order with the kernel k = k 1 on the space I ( k ) ( C 1 ( 0 , + ) ) . Thus, Theorem 4 covers the results presented in Theorem 1, including the spaces of functions used in its formulation.
The next important result is a formula for a composition of the GFI of arbitrary order with the kernel κ and the 1st-level GFD of arbitrary order with the pair of the associated kernels ( k 1 , k 2 ) . This result is referred to as the 2nd fundamental theorem of FC for the 1st-level GFD of arbitrary order.
Theorem 5. 
Let a triple of the kernels ( κ , k 1 , k 2 ) belong to the set L n 1 .
Then the formula
( I ( κ ) 1 L D ( k 1 , k 2 ) f ) ( t ) = f ( t ) j = 0 n 1 d j d t j ( I ( k 2 ) f ) ( 0 ) d n 1 j d t n 1 j ( κ k 1 ) ( t ) , t > 0
holds valid for any function f I ( k 2 ) n ( C 1 ( 0 , + ) ) under the condition that
κ k 1 C 1 n 1 ( 0 , + ) ,
where the space C 1 n 1 ( 0 , + ) is defined as in (49) and
I ( k ) m ( C 1 ( 0 , + ) ) : = { f C 1 ( 0 , + ) : d m d t m ( I ( k ) f ) ( t ) C 1 ( 0 , + ) } .
Proof. 
First, we determine the null-space of the 1st-level GFD 1 L D ( k 1 , k 2 ) . Because the GFI (16) is an injection [20], we obtain the following chain of relations:
( 1 L D ( k 1 , k 2 ) f ) ( t ) = 0 d n d t n I ( k 2 ) f = 0 ( I ( k 2 ) f ) ( t ) = ( k 2 * f ) ( t ) = j = 0 n 1 a j h j + 1 ( t ) .
Now we apply the GFI with the kernel κ k 1 to the relation deduced above and obtain
( I κ * k 1 f ) ( t ) = ( κ k 1 k 2 f ) ( t ) = ( j = 0 n 1 a j h j + 1 ( · ) κ k 1 ) ( t )
( { 1 } < n > f ) ( t ) = ( j = 0 n 1 a j h j + 1 ( · ) κ k 1 ) ( t ) .
n-times differentiation of the last formula that is allowed under the condition (88) leads to the representation
f ( t ) = j = 0 n 1 a j d n 1 j d t n 1 j ( κ k 1 ) ( t )
and we arrive at the following description of the null-space of the 1st-level GFD of arbitrary order:
k e r 1 L D ( k 1 , k 2 ) = { f : f ( t ) = j = 0 n 1 a j d n 1 j d t n 1 j ( κ k 1 ) ( t ) , a j R , j = 0 , 1 , , n 1 } .
Now we introduce an auxiliary function ϕ as follows:
ϕ ( t ) : = ( I ( κ ) 1 L D ( k 1 , k 2 ) f ) ( t ) .
Because of the inclusion f I ( k 2 ) n ( C 1 ( 0 , + ) ) , the 1st-level GFD D ( k 1 , k 2 ) f does exist and is a function from the space C 1 ( 0 , + ) . Thus, we obtain the inclusion ϕ I κ ( C 1 ( 0 , + ) ) . The application of Theorem 4 results in the formula
( 1 L D ( k 1 , k 2 ) ϕ ) ( t ) = ( 1 L D ( k 1 , k 2 ) I ( κ ) 1 L D ( k 1 , k 2 ) f ) ( t ) = ( 1 L D ( k 1 , k 2 ) f ) ( t )
which can be rewritten as follows:
( 1 L D ( k 1 , k 2 ) ( ϕ f ) ) ( t ) = 0 , t > 0 .
Due to Equation (90) for the null-space of the 1st level GFD, we obtain the representation
ϕ ( t ) = f ( t ) + j = 0 n 1 a j d n 1 j d t n 1 j ( κ k 1 ) ( t ) , t > 0 .
Applying the GFI I ( k 2 ) to the left-hand side of (94), i.e., to the function ϕ , leads to the formula
( I ( k 2 ) ϕ ) ( t ) = ( I ( k 2 ) I ( κ ) D ( k 1 , k 2 ) f ) ( t ) = ( I ( k 2 ) I ( κ ) I ( k 1 ) d n d t n I ( k 2 ) f ) ( t ) =
( k 2 κ k 1 ( d n d t n I ( k 2 ) f ) ) ( t ) = ( { 1 } < n > * ( d n d n I ( k 2 ) f ) ) ( t ) =
( I 0 + n d n d t n I ( k 2 ) f ) ( t ) = ( I ( k 2 ) f ) ( t ) j = 0 n 1 d j d t j ( I ( k 2 ) f ) ( 0 ) h j + 1 ( t ) .
Thus, we arrive at the representation
( I ( k 2 ) ϕ ) ( t ) = ( I ( k 2 ) f ) ( t ) j = 0 n 1 d j d t j ( I ( k 2 ) f ) ( 0 ) h j + 1 ( t ) .
Now we apply the GFI I ( k 2 ) to the right-hand side of Equation (94):
( I ( k 2 ) ϕ ) ( t ) = ( I ( k 2 ) f ( · ) + j = 0 n 1 a j d n 1 j d t n 1 j ( κ k 1 ) ( · ) ) ( t ) =
( I ( k 2 ) f ) ( t ) + j = 0 n 1 a j d n 1 j d t n 1 j ( k 2 κ k 1 ) ( t ) =
( I ( k 2 ) f ) ( t ) + j = 0 n 1 a j d n 1 j d t n 1 j { 1 } < n > ( t ) = ( I ( k 2 ) f ) ( t ) + j = 0 n 1 a j h j + 1 ( t ) .
The last formula along with Equation (95) leads to the following expressions for the coefficients a j :
a j = d j d t j ( I ( k 2 ) f ) ( 0 ) , j = 0 , 1 , , n 1 .
The 2nd fundamental theorem of FC for the 1st-level GFD of arbitrary order (Equation (87) from Theorem 5) immediately follows from Equations (91), (94) and (96), and that completes the proof. □
It is worth mentioning that Theorem 5 contains a very general result that covers many known formulas for particular cases of the 1st-level GFD of arbitrary order.
If we set k 2 = h 0 in Equation (87), then the condition ( κ k 1 k 2 ) ( t ) = { 1 } < n > ( t ) is reduced to the condition ( κ k 1 ) ( t ) = { 1 } < n > ( t ) . Thus, ( κ , k 1 ) L n and the 1st-level GFD is the regularized GFD with the kernel k 1 . The space of functions I ( k 2 ) n ( C 1 ( 0 , + ) ) becomes the space C 1 n ( 0 , + ) , and the condition (88) from Theorem 5 is automatically satisfied:
( κ k 1 ) ( t ) = { 1 } < n > ( t ) C 1 n 1 ( 0 , + ) .
Equation (87) then takes the form
( I ( κ ) * D ( k 1 ) f ) ( t ) = f ( t ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( t ) ,
that is exactly Equation (67) for the regularized GFD with the kernel k 1 from Theorem 2.
In the case of the power law kernels κ ( t ) = h α ( t ) and k 1 ( t ) = h n α ( t ) , n 1 < α n , the regularized GFD is the Caputo fractional derivative * D 0 + α of the order α , and Equation (97) takes the well-known form
( I 0 + α * D 0 + α f ) ( t ) = f ( t ) j = 0 n 1 f ( j ) ( 0 ) h j + 1 ( t ) .
For k 1 = h 0 , the 1st-level GFD is reduced to the GFD with the kernel k 2 , where ( κ , k 2 ) L n . Taking into account the relation κ k 1 = κ h 0 = κ , Equation (87) can be rewritten as follows:
( I ( κ ) D ( k 2 ) f ) ( t ) = f ( t ) j = 0 n 1 d j d t j I ( k 2 ) f ( 0 ) d n j 1 d t n j 1 κ ( t ) , t > 0 .
According to the results presented in Theorem 3, this formula is valid under the condition κ C 1 n 1 ( 0 , + ) . Because of the inclusion ( κ , k 2 ) L n , the space I ( k 2 ) n ( C 1 ( 0 , + ) ) mentioned in Theorem 5 becomes the space C 1 , ( k 2 ) 1 ( 0 , + ) from Theorem 3.
In the case of the power law kernels κ ( t ) = h α ( t ) and k 2 ( t ) = h n α ( t ) , n 1 < α n , the GFD with the kernel k 2 is the Riemann–Liouville fractional derivative D 0 + α of the order α . Equation (98) in this case is well known:
( I 0 + α D 0 + α f ) ( t ) = f ( t ) j = 0 n 1 d j d t j I 0 + n α f ( 0 ) h α n + j + 1 ( t ) , t > 0 .
Another important particular case of Equation (87) concerns the Hilfer fractional derivative (1). To generate the Hilfer fractional derivative, we put the kernels ( κ , k 1 , k 2 ) that are the power law functions defined as in (78) into the formula for the 1st-level GFD of arbitrary order. For these kernels, the GFIs from the definition of the 1st-level GFD of arbitrary order are reduced to the Riemann–Liouville fractional integrals:
( I ( k 1 ) f ) ( t ) = ( I 0 + γ f ) ( t ) , ( I ( k 2 ) f ) ( t ) = ( I 0 + n α γ f ) ( t ) .
Moreover, for the power law kernels, the relations
( h α h γ ) ( t ) = h α + γ ( t ) , d d t h α ( t ) = h α 1 ( t )
hold true. Substituting the formulas presented above into Equation (87), we arrive at the 2nd fundamental theorem for the Hilfer fractional derivative (1) with the order α ( n 1 < α n , n N ) and type β ( 0 β 1 ):
( I 0 + α D 0 + α , β f ) ( t ) = f ( t ) j = 0 n 1 d j d t j ( I 0 + ( 1 β ) ( n α ) f ) ( 0 ) h j ( 1 β ) ( n α ) ( t ) , t > 0 .
This formula was derived for the first time in [3].
Finally, we mention that Equation (87) can be rewritten in terms of the projector operator P 1 L of the 1st-level GFD of arbitrary order:
( P 1 L f ) ( t ) : = f ( t ) ( I ( κ ) 1 L D ( k 1 , k 2 ) f ) ( t ) = j = 0 n 1 d j d t j ( I ( k 2 ) f ) ( 0 ) d n 1 j d t n 1 j ( κ k 1 ) ( t ) , t > 0 .
The coefficients d j d t j ( I ( k 2 ) f ) ( 0 ) by the functions d n 1 j d t n 1 j ( κ k 1 ) ( t ) , j = 0 , 1 , , n 1 at the right-hand side of Equation (100) determine the form of the natural initial conditions for the fractional differential equations with the 1st-level GFDs of arbitrary order. In particular, it is well known that the initial conditions for the fractional differential equations with the Caputo fractional derivatives are posed in the form y ( j ) ( 0 ) = a j , j = 0 , 1 , , n 1 , whereas the initial conditions for the fractional differential equations with the Riemann–Liouville fractional derivatives are formulated as follows: d j d t j I 0 + n α f ( 0 ) = a j , j = 0 , 1 , , n 1 .

4. Conclusions and Directions for Further Research

In this paper, the 1st-level GFDs of arbitrary order were defined and investigated for the first time in the FC literature. These derivatives can be interpreted as a generalization of the Hilfer derivative to the case of some general kernels of the Sonin type.
To define the 1st-level GFDs of arbitrary order, a suitable generalization of the Sonin condition for their kernels was introduced in the form
( κ k 1 k 2 ) ( t ) = h n ( t ) = t n 1 ( n 1 ) ! , n N .
Then we specified a set of the kernels that satisfy this condition and possess an integrable singularity of the power law type at the origin. This set was denoted by L n 1 (kernels of the 1st level general fractional derivatives of the generalized order from the interval ( n 1 , n ) ).
The 1st-level GFDs of arbitrary order are integro-differential operators of the convolution type with the kernels from the set L n 1 defined by
( 1 L D ( k 1 , k 2 ) f ) ( t ) = I ( k 1 ) d n d t n I ( k 2 ) f ( t ) ,
where I ( k 1 ) and I ( k 2 ) are the integral operators of convolution type with the kernels k 1 and k 2 , respectively.
The main advantage of these derivatives is that both the GFDs of arbitrary order of the Riemann–Liouville type and the regularized GFDs of arbitrary order considered in the literature so far are their particular cases. In this paper, some important properties, including the 1st and the 2nd fundamental theorems of fractional calculus for the 1st level GFDs of arbitrary order were formulated and proved.
As for the directions for further research, the notion of the 1st-level GFDs can be extended to the case of the 2nd- and even the mth-level GFDs following the procedure presented in [5] for the case of the power law kernels. For instance, the kernels of the mth-level GFDs of the generalized order from the interval ( 0 , 1 ) should satisfy the generalized Sonin condition in the form
( κ k 1 k m + 1 ) ( t ) = h 1 ( t ) = { 1 } , t > 0 .
As in the case of the 1st-level GFDs, a set L 1 m of the kernels that satisfy this condition and possess an integrable singularity of power law type at the origin can be defined (kernels of the mth-level GFDs of the generalized order from the interval ( 0 , 1 ) ).
A natural definition of the mth-level GFDs with the kernels from L 1 m is as follows:
( m L D ( k 1 , , k m + 1 ) f ) ( t ) : = I ( k 1 ) d d t I ( k n ) d d t I ( k m + 1 ) f ( t ) .
By analogy with the construction presented above, the mth-level GFDs of the generalized order from the interval ( n 1 , n ) with the kernels from the set L n m can be also defined and studied.
All these and further related topics will be considered elsewhere.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hilfer, R. Fractional calculus and regular variation in thermodynamics. In Applications of Fractional Calculus in Physics; Hilfer, R., Ed.; World Scientific: Singapore, 2000; pp. 429–463. [Google Scholar]
  2. Hilfer, R. Mathematical and physical interpretations of fractional derivatives and integrals. In Handbook of Fractional Calculus with Applications; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; Volume 1: Basic Theory, pp. 47–86. [Google Scholar]
  3. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for the solution of fractional differential equations with generalized Riemann-Liouville fractional derivatives. Fract. Calc. Appl. Anal. 2009, 1, 299–318. [Google Scholar]
  4. Djrbashian, M.M.; Nersessian, A.B. Fractional derivatives and the Cauchy problem for differential equations of fractional order. Izv. Akad. Nauk Armjan. SSR. Ser. Mat. 1968, 3, 1–29. (In Russian) [Google Scholar] [CrossRef]
  5. Luchko, Y. Fractional derivatives and the fundamental theorem of Fractional Calculus. Fract. Calc. Appl. Anal. 2019, 23, 939–966. [Google Scholar] [CrossRef]
  6. Sonine, N. Sur la généralisation d’une formule d’Abel. Acta Math. 1884, 4, 171–176. [Google Scholar] [CrossRef]
  7. Abel, N.H. Auflösung einer mechanischen Aufgabe. J. Die Reine Angew. Math. 1826, 1, 153–157. [Google Scholar]
  8. Kochubei, A.N. General fractional calculus, evolution equations, and renewal processes. Integr. Equa. Operator Theory 2011, 71, 583–600. [Google Scholar] [CrossRef]
  9. Kochubei, A.N. General fractional calculus. In Handbook of Fractional Calculus with Applications; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany, 2019; Volume 1: Basic Theory, pp. 111–126. [Google Scholar]
  10. Kochubei, A.N. Equations with general fractional time derivatives. Cauchy problem. In Handbook of Fractional Calculus with Applications; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany, 2019; Volume 2: Fractional Differential Equations, pp. 223–234. [Google Scholar]
  11. Kochubei, A.N.; Kondratiev, Y. Growth Equation of the General Fractional Calculus. Mathematics 2019, 7, 615. [Google Scholar] [CrossRef]
  12. Schilling, R.L.; Song, R.; Vondracek, Z. Bernstein Functions. Theory and Application; De Gruyter: Berlin, Germany, 2010. [Google Scholar]
  13. Luchko, Y.; Yamamoto, M. The general fractional derivative and related fractional differential equations. Mathematics 2020, 8, 2115. [Google Scholar] [CrossRef]
  14. Sin, C.-S. Well-posedness of general Caputo-type fractional differential equations. Fract. Calc. Appl. Anal. 2018, 21, 819–832. [Google Scholar] [CrossRef]
  15. Sin, C.-S. Cauchy problem for general time fractional diffusion equation. Fract. Calc. Appl. Anal. 2020, 23, 1545–1559. [Google Scholar] [CrossRef]
  16. Janno, J.; Kasemets, K. Identification of a kernel in an evolutionary integral equation occurring in subdiffusion. J. Inverse Ill-Posed Probl. 2017, 25, 777–798. [Google Scholar] [CrossRef]
  17. Kinash, N.; Janno, J. Inverse problems for a generalized subdiffusion equation with final overdetermination. Math. Model. Anal. 2019, 24, 236–262. [Google Scholar]
  18. Kinash, N.; Janno, J. An Inverse Problem for a Generalized Fractional Derivative with an Application in Reconstruction of Time- and Space-Dependent Sources in Fractional Diffusion and Wave Equations. Mathematics 2019, 7, 1138. [Google Scholar] [CrossRef]
  19. Bazhlekov, I.; Bazhlekova, E. A predictor-corrector numerical approach to equations with general fractional derivative. Int. J. Appl. Math. 2022, 35, 693–709. [Google Scholar] [CrossRef]
  20. Luchko, Y. General Fractional Integrals and Derivatives with the Sonine Kernels. Mathematics 2021, 9, 594. [Google Scholar] [CrossRef]
  21. Luchko, Y. General Fractional Integrals and Derivatives of Arbitrary Order. Symmetry 2021, 13, 755. [Google Scholar] [CrossRef]
  22. Luchko, Y. Operational Calculus for the general fractional derivatives with the Sonine kernels. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  23. Luchko, Y. Special Functions of Fractional Calculus in the Form of Convolution Series and Their Applications. Mathematics 2021, 9, 2132. [Google Scholar] [CrossRef]
  24. Luchko, Y. Convolution series and the generalized convolution Taylor formula. Fract. Calc. Appl. Anal. 2022, 25, 207–228. [Google Scholar] [CrossRef]
  25. Luchko, Y. Fractional Differential Equations with the General Fractional Derivatives of Arbitrary Order in the Riemann-Liouville Sense. Mathematics 2022, 10, 849. [Google Scholar] [CrossRef]
  26. Luchko, Y. The 1st level general fractional derivatives and some of their properties. J. Math. Sci. 2022. [Google Scholar] [CrossRef]
  27. Tarasov, V.E. General Fractional Calculus: Multi-Kernel Approach. Mathematics 2021, 9, 1501. [Google Scholar] [CrossRef]
  28. Tarasov, V.E. General Fractional Dynamics. Mathematics 2021, 9, 1464. [Google Scholar] [CrossRef]
  29. Tarasov, V.E. General Fractional Vector Calculus. Mathematics 2021, 9, 2816. [Google Scholar] [CrossRef]
  30. Tarasov, V.E. General Non-Local Continuum Mechanics: Derivation of Balance Equations. Mathematics 2022, 10, 1427. [Google Scholar] [CrossRef]
  31. Tarasov, V.E. Nonlocal Probability Theory: General Fractional Calculus Approach. Mathematics 2022, 10, 3848. [Google Scholar] [CrossRef]
  32. Tarasov, V.E. Nonlocal statistical mechanics: General fractional Liouville equations and their solutions. Phys. Stat. Mech. Its Appl. 2023, 609, 128366. [Google Scholar] [CrossRef]
  33. Tarasov, V.E. General non-Markovian quantum dynamics. Entropy 2021, 23, 1006. [Google Scholar] [CrossRef]
  34. Tarasov, V.E. Nonlocal classical theory of gravity: Massiveness of nonlocality and mass shielding by nonlocality. Eur. Phys. J. Plus 2022, 137, 1336. [Google Scholar] [CrossRef]
  35. Garrappa, R.; Giusti, A.; Mainardi, F. Variable-order fractional calculus: A change of perspective. Commun. Nonlinear Sci. Numer. Simul. 2021, 102, 10590. [Google Scholar] [CrossRef]
  36. Hanyga, A. A comment on a controversial issue: A Generalized Fractional Derivative cannot have a regular kernel. Fract. Calc. Anal. Appl. 2020, 23, 211–223. [Google Scholar] [CrossRef]
  37. Samko, S.G.; Cardoso, R.P. Integral equations of the first kind of Sonine type. Intern. J. Math. Sci. 2003, 57, 3609–3632. [Google Scholar] [CrossRef]
  38. Zacher, R. Boundedness of weak solutions to evolutionary partial integro-differential equations with discontinuous coefficients. J. Math. Anal. Appl. 2008, 348, 137–149. [Google Scholar] [CrossRef]
  39. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series: Direct Laplace Transforms; Gordon & Breach: New York, NY, USA, 1992; Volume 4. [Google Scholar]
  40. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series: Inverse Laplace Transforms; Gordon & Breach: New York, NY, USA, 1992; Volume 5. [Google Scholar]
  41. Luchko, Y.; Gorenflo, R. An operational method for solving fractional differential equations. Acta Math. Vietnam. 1999, 24, 207–234. [Google Scholar]
  42. Hilfer, R.; Luchko, Y. Desiderata for Fractional Derivatives and Integrals. Mathematics 2019, 7, 149. [Google Scholar] [CrossRef]
  43. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives. Theory and Applications; Gordon and Breach: New York, NY, USA, 1993. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luchko, Y. On the 1st-Level General Fractional Derivatives of Arbitrary Order. Fractal Fract. 2023, 7, 183. https://doi.org/10.3390/fractalfract7020183

AMA Style

Luchko Y. On the 1st-Level General Fractional Derivatives of Arbitrary Order. Fractal and Fractional. 2023; 7(2):183. https://doi.org/10.3390/fractalfract7020183

Chicago/Turabian Style

Luchko, Yuri. 2023. "On the 1st-Level General Fractional Derivatives of Arbitrary Order" Fractal and Fractional 7, no. 2: 183. https://doi.org/10.3390/fractalfract7020183

Article Metrics

Back to TopTop