1. Introduction
The essential idea of fractional calculus is to consider functional operators that can be seen as extensions to non-integer orders of the basic repeated differentiation and integration operators [
1]. The first question arising from this idea is: how can such extensions be defined? Interestingly, this question has no unique answer, as there is a rich tapestry of different operators which can be called ”fractional derivatives” or “fractional integrals” from one viewpoint or another.
Among these operators, many are expressed as convolutions with some particular kernel function that determines the nature of the operator. The classical choice of kernel function is a fractional power, giving rise to the original Riemann–Liouville fractional calculus as well as other related definitions such as Caputo and Hilfer [
2,
3]. Other popular choices of kernel function have included Mittag-Leffler functions, products of power and exponential functions, hypergeometric functions, and many others [
4,
5].
In response to the flood of new definitions using assorted kernel functions, mathematicians have realised that many of these are not so different from each other, and that it makes more sense to study general convolution operators in order to discover mathematical properties [
5]. The general operators thus proposed have included convolutions with analytic kernels [
6,
7], Sonine kernels [
8,
9,
10], and other types of general kernels [
11,
12,
13]. Among the studies so far, a particularly interesting direction has been the study of Sonine kernels: these are rigorously posed with suitable function spaces clearly defined, but they do not depend on any explicit real or complex number playing the role of a “fractional order”, instead having two kernel functions which satisfy a sort of inversion-type relation making them suitable for defining fractional integral and derivative operators with fundamental theorem of calculus type relations.
A point of debate within the fractional calculus community has concerned the operators with non-singular kernel functions: that is, functions that do not have a singularity (blowup) at the endpoint of the integral. For example, the function
has a singularity at
for positive values of the parameter
, but the function
is non-singular. Some scientists have claimed that operators with non-singular kernels are not suitable to be called “fractional” [
14,
15], while others have claimed that this approach is faulty [
16]. We prefer to sidestep this terminology debate, pointing out instead that the operators with non-singular kernels have found real-world applications [
17,
18] and are therefore also worth studying from the mathematical point of view [
19], regardless of whether we call them “fractional” or not. Many studies of integral equations, not necessarily in a fractional setting, have not put any particular emphasis on whether a kernel is singular or non-singular [
20,
21].
So far, most publications related to such operators with non-singular kernels have been written by engineers and applied scientists who are concerned mainly with their applications. We, as mathematicians, have noticed some interesting mathematical properties of these operators, such as their function spaces and inversion properties, which we believe deserve to be highlighted and discussed. As a detailed mathematical investigation of these operators is currently missing in the literature, we are taking the opportunity to conduct one and present our findings to be used by mathematicians and scientists alike.
Our work is arranged as follows.
Section 2 serves to correct a (minor but important) misconception regarding function spaces appearing in the previous literature;
Section 3 focuses on how derivative-type operators with non-singular kernels can be rewritten so as not to need differentiability;
Section 4 uses an abstract algebra viewpoint to discover a quasi-inverse relation which can serve as an alternative to the Sonine relation; and finally,
Section 5 provides a summary and concluding remarks. We note once again that we do not intend to enter the debate on whether or not it is appropriate to include operators with non-singular kernels as “fractional derivatives”. This is not our business here: we are simply presenting some objective mathematical facts concerning certain operators and the associated function spaces, and we will allow the results to speak for themselves.
2. A Misconception on Function Spaces
The classical
space of Dimovski [
22,
23] is defined as follows:
where
is arbitrary.
In his construction of the fractional calculus with Sonine kernels [
9] and its generalisation [
10], Luchko made use of a function space
inspired by this
. The general space
was defined in ([
10] Equation (23)), as a natural generalisation of the space
first defined in ([
9] Equation (33)), as follows:
where
in
. The purpose of using
for the fractional operators with Sonine kernels, instead of simply
as was done for the classical fractional operators [
24], was in order to obtain “an integrable singularity at the point zero” [
9]: Reference [
10] explicitly states that functions in
“possess the integrable singularities of the power function type at the origin”. We show here that this is not necessarily the case, and that a different function space must be used if it is desired to exclude non-singular functions in
from the space of kernels.
Theorem 1. For any in , the space is exactly the space .
Proof. It is clear that
. For the converse, let
be arbitrary, say
with
and
. If
, then we are done. Otherwise, let
and
, so that
and
where
so that
. Therefore, the function
f is in the space
, which completes the proof. □
Example 1. The non-singular functions and are in the space because they can be expressed as and with defined by As a result, we see that the function space used by Luchko [
9,
10] for the operators with Sonine kernels does not explicitly exclude non-singular kernels as it was intended to. However, the Sonine condition itself excludes them because if
k and
are both in
, then
with
, so
cannot be either
or
(
) for the Sonine [
9] or generalised Sonine [
10] kernels. So, if
and
k are a pair of Sonine kernels in the space
defined in ([
9] Definition 2) or the space
defined in ([
10] Definition 2), then at least one of
and
k must have a singularity at zero. For this reason, it seems that the results of Luchko [
9,
10] on Sonine kernels are not seriously affected by the realisation of our Theorem 1 above.
Is there any way to replace Luchko’s definition of the
space with another space that successfully excludes all non-singular functions? An obvious choice would be the following set:
but this set is not a vector space because we can always add two functions in
to obtain a function in
. It is also not closed under convolution (rng multiplication), and indeed we can state a stronger result, as in the following theorem.
Theorem 2. Let , and let A be a subset of that is non-trivial (i.e., contains at least one non-zero function). If A contains only singular functions (in other words, if ), or more generally if for any fixed , then A cannot be closed under convolution, and A cannot be an ideal or subrng of .
Proof. It is well known [
25,
26] that the convolution of a function in
and a function in
is always a function in
when
. Therefore, repeatedly taking the convolution of any function in
with itself will eventually yield a function in
for arbitrarily large
. Therefore, any subrng or ideal, or any set closed under multiplication in
, must contain elements of
for arbitrarily large values of
. Since the
spaces are nested, with
if
, this completes the proof. □
3. Equivalent Formulations of Differential Operators with Non-Singular Kernels
Let us consider a non-singular kernel
and define corresponding derivative-type operators as follows:
in both cases for suitable functions
f. (We will discuss suitable function spaces for these operators shortly.)
One of the distinguishing features of non-singular kernels, as discussed in ([
27] Section 4.1) is the fact that operators such as (
1) and (
2) with non-singular kernels can be rewritten without any explicit differentiation and without needing any differentiability assumptions on the function
f. This is done using, respectively, the fundamental theorem of calculus and the method of integration by parts:
This manipulation works whenever
k and
f are both in the space
but the final expressions (
3) and (
4) are valid for any
or any
, respectively, as long as the fixed kernel function
k is in a suitable space such as
. We state this result formally as follows.
Theorem 3. If the kernel function k is in the space , then the formulae on the right-hand sides of Equations (1) and (2) give well-defined mappings and from into . If the kernel function k is in the space , then the formulae on the right-hand sides of Equations (3) and (4) give well-defined mappingsandand these operators are identical to those defined by (1) and (2) when f is in the subspace such that all are defined. Therefore, assuming that , the suitable (largest) function spaces on which to define the operators and are the spaces and , respectively. In particular, no differentiability conditions are required on f in order to define and as functions in .
Proof. Firstly, when
and
, Formula (
2) defines a function in
since
k and
are in the space
(which is closed under convolution), and Formula (
1) defines a function in
by ([
25] Theorem 2.5).
Now let us fix
and consider expressions (
3) and (
4). Since
, the integral in (
3) defines a function in
for any
, while
exists by ([
25] Fact (3)), so we are done for (
3). The difference between (
3) and (
4) is only in the term
, so the latter defines a function in
whenever
and
exists, which is the same as requiring
.
Finally, the equivalence of (
1) and (
2) with (
3) and (
4) is clear from the above arguments, using the fundamental theorem of calculus and integration by parts, which are valid when
k and
f are both differentiable on
with
existing and
. □
Remark 1. Note that the rewriting of Equations (1) and (2) to Equations (3) and (4) is only possible when the kernel k is non-singular, otherwise the terms in (3) and (4) cannot be meaningful. Thus, our reformulation of Riemann–Liouville and Caputo derivative-type operators, to a form that does not require differentiability conditions, is only possible for operators defined using non-singular kernels. Example 2. Taking the kernel function k to be a power function with exponent greater than one, in the formwe find that sinceIn this case, the operators and , although denoted as derivatives, are actually the Riemann–Liouville fractional integrals of order . This reminds us that the Riemann–Liouville fractional integrals of order greater than one are operators with non-singular kernels. Example 3. Taking the kernel function k to be a constant times a Mittag-Leffler function, in the formwith and B being a suitable function satisfying , we certainly have sincewhich is times a continuous function in . In this case, the operators and are precisely the Atangana–Baleanu derivatives, respectively of Riemann–Liouville type and Caputo type. Here, Theorem 3 tells us, as already known [28], that these Atangana–Baleanu derivatives can be defined on larger function spaces than most fractional derivatives, without requiring any differentiability assumptions. The reformulation (4) of the Caputo-type operator gives, in this case, precisely the formula that is called the “modified ABC derivative” in [29]. Example 4. More generally, taking the kernel function k to be a three-parameter Mittag-Leffler function [30] with the second parameter equal to one, in the formwith constant and , we again have sincewhich is times a continuous function in as in the previous example. In this case, the operators and are precisely the Prabhakar fractional derivatives with second parameter equal to one, respectively of Riemann–Liouville type and Caputo type. These particular types of Prabhakar operators were already discussed in ([27] Section 4.1) as the non-singular case of Prabhakar fractional calculus, and the only case where Prabhakar derivatives can be defined on large enough function spaces to not need any differentiability assumptions. The above results can be extended to operators defined using multiple (repeated) derivatives together with convolution integrals. Let us consider the following operators, for
and for suitable functions
k and
f:
and try to rewrite Formulas (
5) and (
6) using the fundamental theorem of calculus and integration by parts as before.
Repeatedly differentiating the integral in (
5) and using the fundamental theorem of calculus, we find:
and finally
where we note that if
, then no differentiation of
f needs to take place in the terms outside the integral.
Repeatedly using integration by parts in (
6), we find:
where in this case some derivatives of
f are involved in the boundary terms regardless of the nature of
k.
Theorem 4. If the kernel function k is in the space , then the formulae on the right-hand sides of Equations (5) and (6) are well-defined functions in provided that f is such that (for ) or (for ), respectively. If the kernel function k is in the space and , then the formulae on the right-hand sides of Equations (7) and (8) give well-defined mappingsandand these operators are identical to those defined by (5) and (6) when f is such that all are defined. Therefore, assuming that and , the suitable (largest) function spaces on which to define the operators and are respectively the space and the subspace of consisting of functions times differentiable at 0. In particular, no differentiability conditions are required on f in order to define as a function in .
Proof. Firstly, it is clear that (
5) defines a function in
when
, and that (
6) defines a function in
when
and
.
If
and
, then
and the finite sum in (
7) reduces to a single term
, so (
7) defines a function in
for all
. The difference between (
7) and (
8) is only in the terms
, so the latter defines a function in
whenever
and
exist.
Finally, the equivalence of (
5) and (
6) with (
7) and (
8) is clear from the above arguments, using the fundamental theorem of calculus and integration by parts, which are valid when
k and
f are both
on
with derivatives up to
th order existing at 0. □
Example 5. Taking the kernel function k to be a power function with exponent greater than n, in the formwe find that with , sinceIn this case, the operators and , although notated as derivatives, are actually the Riemann–Liouville fractional integrals of order . Once again, as in Example 2, this reminds us that the Riemann–Liouville fractional integrals of order greater than one are operators with non-singular kernels. Remark 2. Note that the rewriting of Equations (5) and (6) to Equations (7) and (8) is only possible when the kernel k is non-singular, otherwise the terms in (7) and (8) cannot be meaningful. It is known [25] that the space is contained in . Thus, once again, our reformulation of Riemann–Liouville and Caputo derivative-type operators, to a form that requires less strong differentiability conditions, is only possible for operators defined using non-singular kernels. 4. Mikusiński’s Operational Calculus for Non-Singular Kernels
In the above, we have considered what happens when an operator of fractional derivative type (a combination of a derivative or repeated derivative operator with a convolution integral operator) is permitted to have a non-singular kernel. A rigorous examination of the required assumptions and function spaces for such operators revealed exactly which differentiability assumptions could be relaxed. However, so far we have not considered any sort of inversion relation, as is commonly seen between derivatives and integrals.
In Luchko’s work on Sonine kernels [
9,
10], he imposed Sonine-type conditions in order to ensure inversion relations between his operators. If two kernels
and
k satisfy
then the integral operator
defined by
obeys a fundamental theorem of calculus together with the differential operators
and
defined by (
1) and (
2). Similarly, if
and
k instead satisfy
then the integral operator
defined by
obeys a fundamental theorem of calculus together with the differential operators
and
defined by (
5) and (
6).
But some types of operators with non-singular kernels have meaningful inversion relations between integrals and derivatives even without satisfying a Sonine condition. How can we extend the invertibility condition beyond the Sonine one in order to capture these different behaviours?
To answer this, we need some concepts and results from abstract algebra.
Proposition 1 ([
31])
. The space under the operations of addition and convolution forms a commutative rng (ring without identity) as well as a vector space over . Its field of fractions is both a field and a vector space over . Proof. This is a well-known result in the theory of Mikusiński’s operational calculus. □
Theorem 5. The algebra generated by and the identity element from , i.e., the set of elements with and (where I is the multiplicative identity from ) is a commutative ring with identity as well as a vector space over .
Proof. This is a standard variant of the well-known Dorroh extension [
32] that allows every rng to be extended to a ring. Under the operations
of addition and
of scalar multiplication and
of multiplication, it is immediate that the new space satisfies all the axioms to be both a vector space and a commutative ring with identity. □
Theorem 6. All elements of have multiplicative inverse elements in .
Proof. One of the main results of Heatherly and Huffmann ([
33] Theorem 2.1) was that every element
(considered as a commutative rng under addition and convolution) has a quasi-inverse
, i.e., such that
Adding
to both sides of this equation gives
Therefore, for any
and any
, the element
satisfies
where
is the quasi-inverse of
. □
Let us now discuss why the above results are important.
In Mikusiński’s operational calculus, a function space is interpreted as an rng, in which functions can be added, subtracted, and multiplied, but not divided, and there is a zero (additive identity) but not a one (multiplicative identity). The construction of the Mikusiński field (the field of fractions of this rng, using a theorem of Titchmarsh to determine that there are no zero divisors) enables division of functions, but at the expense of moving into a more abstract space: dividing two functions does not necessarily yield another function, but an abstract field element, called an “operator” by Mikusiński [
31]. In this field, there is a multiplicative identity element, whose convolution with any function yields the same function back again; as a generalised function concept, this identity element is analogous to the Dirac delta from distribution theory.
The full structure of the Mikusiński field is necessary in order to allow inverses of all elements of the starting rng, but it is not necessary if all we need is a multiplicative identity element. Theorem 5 shows that adjoining just one element (and generating a new vector space) is enough to get a commutative ring with a multiplicative identity, and Theorem 6 shows that almost all elements of this ring have inverses: it is not a field, but it is almost as good as being a field.
Furthermore, elements of the ring can be understood as operators without needing to resort to abstract “generalised functions”. Any element can be interpreted as an operator from to in a very natural way: mapping any to . This means that is a very natural setting for defining operators if we do not strictly require that everything must have a multiplicative inverse.
Theorem 7 (Quasi-inverse functions give inverse functional operators).
For any non-singular kernel and any nonzero constant , the operator defined byhas an inverse operator defined by another function , namelywhere for all and the relation between the kernel functions can be expressed as Proof. This is essentially a restatement of Theorem 6. The operator A is precisely multiplication by the element , which has an inverse by Theorem 6. Defining , we have the operator as stated being multiplication by .
Since embeds as a subalgebra in , their fields of fractions similarly embed, and so embeds as a subalgebra of the corresponding algebra generated by I and . Therefore, if the product of and is the identity in , then it is also the identity in the field of fractions of . So, the product of these two elements with any is again f, since multiplication in is commutative and associative.
The inversion relation in reduces to in after subtracting I from both sides. □
Definition 1. We can say two pairs and , each consisting of a constant in and a function in , are quasi-inverse to each other if and Equation (9) is satisfied, i.e., if and . This is a sufficient condition to get well-defined integral operators that are inverse to each other as in Theorem 7. This is essentially an analogue, for non-singular kernels, of the Sonine relation. We have a direct relation between functions κ and k, which can be formulated in terms of their convolutions in a suitable function space, which gives rise to an inversion relation between corresponding operators.
Remark 3. Note that the inversion relation given by Theorem 7 is symmetric: the operators A and are two-sided inverses of each other. This is unlike the usual derivative and integral operators (either classical, or fractional with singular kernels), which have asymmetric inversion relations given by the fundamental theorem of calculus. Thus, we observe an interestingly different structure for operators with non-singular kernels, as previously noted in some other works [27,28,34]. Example 6. Taking k to be a constant function in Theorem 7, we find the operator A to be a linear combination of a function with its integral, the same formulation used for the Caputo–Fabrizio integral:How do we find a suitable function κ satisfying the identity (9) when k is a constant function? The answer can be easily found using Mikusiński’s methodology in the field of fractions :so the kernel function is given byThus, we have proved that the quasi-inverse in of a constant function is a constant times an exponential function, and the inverse of the operator (10) is given by:In the usual Caputo–Fabrizio notation [34,35], we would write and , where and is a function satisfying , to find:The operator A is the Caputo–Fabrizio integral, while the operator is the corresponding derivative of Riemann–Liouville type [36], as we can see from a quick computation:Note that this operator is a two-sided inverse to the Caputo–Fabrizio integral A, on the space and therefore on any subspace such as . The corresponding Caputo-type derivative is not an exact algebraic inverse; instead, it has a Newton–Leibniz type formula together with the integral operator, as discussed in detail in [34]. Remark 4. The Atangana–Baleanu integral is a linear combination of a function with its Riemann–Liouville fractional integral. As the Riemann–Liouville integral is singular, this is not interpretable as an operator in our algebra generated from the rng of continuous functions , and so the Atangana–Baleanu operators are not covered by our Theorem 7. Further work will be necessary to ascertain whether operators like those of Atangana–Baleanu can be covered by an extension of the result of Theorem 7.