Intrinsic Discontinuities in Solutions of Evolution Equations Involving Fractional Caputo–Fabrizio and Atangana–Baleanu Operators

There has been considerable recent interest in certain integral transform operators with non-singular kernels and their ability to be considered as fractional derivatives. Two such operators are the Caputo–Fabrizio operator and the Atangana–Baleanu operator. Here we present solutions to simple initial value problems involving these two operators and show that, apart from some special cases, the solutions have an intrinsic discontinuity at the origin. The intrinsic nature of the discontinuity in the solution raises concerns about using such operators in modelling. Solutions to initial value problems involving the traditional Caputo operator, which has a singularity inits kernel, do not have these intrinsic discontinuities.


Introduction
In mathematical models based on evolution equations it is standard to restrict consideration to models whose solutions exist in a space of continuous functions. This is also true in cases of evolution equations with fractional order differential operators, such as the Riemann-Liouville operator [1] and the Caputo operator [2]. Both of these operators are based on integrals with power law kernels that are singular at the origin. In recent years, there has been a great deal of attention focussed on fractional differential operators based on integrals with non-singular kernels. Included in this are the Caputo-Fabrizio (CF) operator [3] and the Atangana-Baleanu in the sense of a Caputo (ABC) operator [4].
The CF and ABC operators have since been employed in numerous modelling applications, including applications to phase transitions [5], fluid and ground water flow [6], cancer treatment [7], and epidemiology [8], among others [9][10][11][12]. Many of the works in this field introduce the model evolution equations by the ad hoc "fractionalisation" of simply replacing integer order derivatives in traditional models with CF or ABC derivatives, without further phenomenological consideration.
On the theoretical side, there has been considerable effort devoted to understanding the interpretations of these differential operators. In a sequence of studies, Tarasov [13][14][15], Ortigueira and Machado [16], and Giusti [17] have shown that the CF differential operator should not be regarded as a fractional order operator and the ABC operator does not extend beyond the Caputo operator. Nevertheless, modelling applications still persist [8], as do numerical studies [18][19][20] and algebraic methods of solution [21]. We also note that, as presented by Hilfer and Luchko [22], there is no absolute agreement on a well defined set of properties to which a fractional derivative must adhere and it is not the aim of this work to interrogate such properties. A greater concern, expressed in the present work, is that CF and ABC operators are not generally suitable for modelling, whenever a solution is sought in a space of continuous functions. In particular, we show that we can construct well formed solutions to initial values problems (IVPs) with CF operators but the solutions have discontinuities at the origin. Many of the works discussed above make extensive use of integral transforms, and more specifically, the Laplace transform. A key point of the present contribution is highlighting that the Laplace transform is not bijective for the class of functions which solve IVPs under the CF operator and careful treatment of the solution near the initial condition is paramount.
The discontinuities in the solution of the IVPs necessitate the more general consideration of the derivatives in the definitions of both the CF and ABC operators. Such generalisations to distributional derivatives are well defined [23]. There is a large body of work concerned with the analysis of such derivatives, both for the case of the CF operator [24] and for the Riemann-Liouville, Caputo, and other fractional derivatives [25][26][27].
The remainder of this paper is organised as follows. We first construct solutions of CF IVPs and show that such solutions must, in general, feature a discontinuity at the origin. We then briefly discuss the impacts of these results on numerical methods for the solution of IVPs involving CF operators and point out a seemingly overlooked simple approach to the numerical evaluation of such equations. Next we repeat this treatment for the ABC operator and again show that solutions, in general, will have a discontinuity at the origin. Finally we consider a more traditional fractional derivative, the Caputo derivative, and show that solutions to Caputo IVPs can not feature such discontinuities.
It is typical and sufficient to take M(α) = 1. The CF operator has been purported to be a fractional derivative when 0 < α < 1 which limits to an integer order derivative as α → 1 − [3,28]. It should also be noted that the derivative in the integral may be considered in the distributional sense; for more information, see [24].
We will consider CF equations in the form of an IVP with Definition 2 (A Caputo-Fabrizio Initial Value Problem). A CF IVP is given by both a CF equation of the form with F(t) a continuous function for t ≥ 0 and an initial value, with u 0 ∈ R.
For now, it will suffice to say that we will consider solutions, u(t), defined over the interval t ∈ [0, ∞) without being overly concerned with the smoothness of such solutions other than asking for the derivative, at least in the distributional sense, to be well defined.
It should be noted that, in general, continuous solutions of this IVP do not exist [29]. Here we will consider a more general form of a solution by taking an ansatz such that where u c is a continuous function, δ is a Dirac delta, and a is an unknown constant. Note that this form of a solution still permits a purely continuous form with a = 0. The form of this Dirac delta is chosen so that we have With the ansatz in Equation (4) it follows that the solution will be of the form We will further simplify this by taking u c (0) = u(0), to give Care has to be taken here due to the fact that the Laplace transform is not bijective, i.e., As such, we can not rely on Laplace transform techniques to find solutions of this form.
Theorem 1. For a CF IVP (Definition 2) assume that a solution in the form of the ansatz (Equation (7)) exists. Then such a solution is given by Proof. The solution is found by assuming that a solution in the form of the ansatz exists, substituting it into the IVP, and showing that the result is then consistent. To find the value of the unknown constant a from the ansatz we first substitute the Equation (4) into Equation (2), which gives for t > 0. Next we take the limit as t approaches 0 from above to give and rearrange to find as F is continuous. To find the continuous part of this solution, we may differentiate the IVP, Equation (10), with respect to t: From Equation (10) we also have Combining these two expressions gives an integer order differential equation for u c , The integral form of this equation is Hence the general form for the solution of the IVP is From Equation (12) we see that the CF IVP only permits a continuous solution in the case where F(0) = 0. This requirement on the existence of continuous solutions has been noted in [29], where the result was obtained via Laplace transforms, and is the case considered in [30]. In all other cases the solution will involve both a continuous component and a step discontinuity at the origin.
This solution of the CF IVP can easily be alternatively verified by taking the CF operator of the solution to recover the original IVP, Here we have simply applied integration by parts. This general solution may alternatively be written as It should be noted that this solution differs from the solution given in other papers, such as [28,31], although it is in agreement with the solution given in [17] for t > 0.

Weakening Continuity Requirements for F(t)
In the above we assumed that F(t) was a continuous function. This is a little restrictive, as it excludes most of the interesting cases where F(t) depends on the function u(t), as u(t) is not a continuous function. To accommodate a discontinuity at the origin we may remove the requirement that F is continuous and assume that we can write where F c is a continuous function and b ∈ R. By following the same methodology as above, we have for t > 0. By taking the limit from above as t tends to 0, we can find the unknown coefficient a.
The continuous part of the solution will again be reduced to the solution of an ODE.
The general solution with a discontinuity at t = 0 for F(t) is thus This solution is completely equivalent to the general solution given above in Equation (17). From this we see that the weakening of the continuity requirement did not effect the solution and we can attempt to solve IVPs of the form with u(0) = u 0 . In the case of IVPs of this form the given solution may only exist in the case that a, found via Equation (22), is real valued. For example with F(u(t), t) = M(α) 1−α u(t), the relation in Equation (22) does not hold with u o = 0, and hence there is no solution of the ansatz form. Furthermore, for a non-linear equation, the resulting ODE for the continuous part of the solution, Equation (23), may not have solutions. Hence we must deal with each non-linear case individually.

Example Solutions for CF Initial Value Problems
We will construct some solutions of simple IVPs to illustrate the forms given above. In each case the validity of the solution can easily be seen by a direct substitution into the original equation. and The solution of this IVP can be found directly via Equation (17) and is Notice that the definition of H(t) ensures that u(0) = u 0 , but the solution has a step at t = 0. This is an illustrative example as it is simple to check against the definition of the CF operator.

Example CF IVP with F(u(t), t) = −u c (t)
It is instructive to consider the differences induced by a discontinuity in F at t = 0. Here we present an example where F is taken to be the continuous part of the solution whilst in the next example we will show the case for F being the full solution.
Consider the CF IVP with and This solution can be found first by solving the ODE for u c , Equation (15), which can be rearranged to obtain, subject to the initial condition u c (0) = u 0 . The continuous part of the solution is thus The discontinuous part of the solution is readily found from Equation (12) and combining the two will give the above solution. Thus the CF IVP has a solution of, We can see that as t → ∞ this solution changes sign and asymptotes to Using the weakened form of the continuity requirement we can consider the IVP of the form and As the right-hand side of this equation is dependent on the solution, we will first find the unknown coefficient from the ansatz via the relation given in Equation (22), with b = a and F c (0) = u c (0) = u 0 . This gives, Following the same procedure as the previous example, we obtain the following ODE for the continuous part of the solution: This can be solved with the initial condition u c (0) = u 0 to give Again, combining the continuous and discontinuous parts of the solution will give the full solution, In contrast to the second example, this solution will remain positive, for u 0 > 0, and will asymptote to 0 as t → ∞.
Again, using the weakened form of the continuity requirement, we can consider the IVP of the form and In this case we see that From the ansatz we also have hence, From Equation (22) we then have the relation This has solutions From Equation (23) we have an ODE for the continuous part of the solution, This ODE has a solution, where W 0 is a Lambert W function [32]. The full solution to the IVP is thus

Numerical Considerations for Equations Involving the CF Operator
Equation (1) purports a memory effect by way of a convolution through time. Discretising the CF operator in this form leads to the unnecessary computation of memory terms. As is shown above, the solution to the IVP (2) may be obtained through the solution of the auxiliary ordinary differential Equation (15) (or in integral form Equation (17)).
The numerics contained within the recent literature are largely restricted to low order numerical methods. From the formulation presented in this work, we suggest that any numerical method appropriate for ODEs may be used to accurately solve IVPs with CF operators, and as such many efficient, highly accurate methods are available to these equations. To the best of the authors' knowledge no preceding work has proposed a numerical method which recovers discontinuous solutions to CF equations. As shown above these equations do not exhibit nontrivial continuous solutions, and as such require numerical methods tailored to recover the discontinuous dynamics.

The Atangana-Baleanu Operator
In a similar manner to the CF operator we will consider another non-singular kernel operator, the Atangana-Baleanu, in the sense of Caputo, (ABC) operator [4].

Definition 3 (The Atangana-Baleanu, in the sense of Caputo, Operator). The ABC operator for
where E α is the Mittag-Leffler function, defined by: x n Γ(αn + 1) .
It is sufficient to take B(α) = 1. The ABC operator is again often seen to be a fractional derivative as we recover an integer order derivative in the case α → 1 − . The use of the ABC operator as a fractional derivative is less contentious than the CF operator as it is non-local in time. As we are considering discontinuous solutions it is again necessary to interpret the derivative in the definition of the ABC operator in a distributional sense.
Again we will consider simple IVPs arising from this operator.

Definition 4 (An ABC Initial Value Problem). An ABC IVP is given by both an ABC equation,
with F(t) a continuous function in t, and an initial condition u(0) = u 0 for some u 0 ∈ R.
To find a solution we consider the ansatz with the integral form of the solution We will further simplify this by taking u c (0) = u(0), to give Theorem 2. For an ABC IVP (Definition 4) assume that a solution in the form of the ansatz (Equation (56)) exists. Then such a solution is given by Proof. Combining Equation (53) and Equation (54) gives for t > 0. Taking the limit as t approaches 0 from above, we obtain Thus the coefficient a is given by provided F(x) is continuous. As u c (t) is a continuous function and the Laplace transform is bijective over the space of continuous functions, Laplace transform techniques can be applied. Therefore we take the Laplace transform of both sides of Equation (58) from t to s domain to obtain having used the result from [1]. By rearranging the equation above, we see that To deal with the s 1−α L{F(t)} on the RHS of Equation (63) we will utilise some results from fractional calculus. The Riemann-Liouville fractional derivative [1] of order 1 − α with 0 ≤ α ≤ 1 is defined by The Laplace transform of the Riemann-Liouville derivative is given by [33], where RL 0 D −α t f (t) is a Riemann-Liouville fractional integral of order α. Furthermore, provided the limits exist, we also have [33,34], and hence provided that lim t→0 + F(t) exists we have The inverse Laplace transform of Equation (63) then gives the following ordinary integro-differential equation for u c (t): The integral form of the solution is obtained immediately from the equation above, the result is Hence the general form for the solution of the IVP is Again, we can alternatively verify this this solution by substituting the solution back to the ABC operator, Note that both sides of Equation (71) We see that inverse Laplace transform of the equation above recovers the IVP: We can note that the continuity requirement on F(t) can be eased in the same manner as the CF operator by considering a discontinuous F(t) such that F(t) = F c (t) + bH(t), where F c (t) is a continuous function. As such, we can attempt to consider IVPs of the form with u(0) = u 0 . Again a solution of the ansatz form will only exist if a real valued constant a can be found such that the following relation holds, where F c , a continuous function, and b, a real valued constant, are found from F(u(t), t) = F c (u(t), t) + bH(t).

Example ABC Initial Value Problems
and The solution of this IVP follows immediately from the Equation (70) and is 3.1.2. Example ABC IVP with F(t) = u c (t) and The solution can be found by independently calculating the continuous and discontinuous parts of the solution. The continuous part of the solution is found by first considering the equation for u c in Laplace space, Equation (63), which can be rearranged to obtain The inverse Laplace transform from s to t domain then gives subject to the initial condition u c (0) = u 0 . The discontinuous part of the solution is found by calculating the coefficient a from Equation (60), giving The solution of the IVP is then as follows:

Example ABC IVP with F(u(t), t) = u(t)
Using the weakened form of the continuity requirement we can consider the IVP of the form and Again the solution is found by considering the continuous and discontinuous parts separately. From Equation (60), we obtain Subject to the initial condition u c (0) = u 0 . The continuous part of the solution can be found via its Laplace transform. Substituting the value for a into Equation (63) gives Inverting the Laplace transform then gives the continuous part of the solution, Combining the continuous and discontinuous parts will then give the full solution,

Singular Kernel Operator Example: Caputo Derivative
Here we attempt to apply the same technique to a fractional derivative with a singular kernel, the Caputo derivative.
Definition 5 (Caputo Derivative). The Caputo derivative of order α with 0 < α < 1, is defined as [2], In order to explore the possibility of discontinuous solutions we will need to allow for the derivative in the Caputo definition to be interpreted in a distributional sense; see [25] for a more detailed exposition of the use of distributional derivatives in Caputo derivatives. We again consider a simple IVP.

Definition 6 (Caputo IVP). A Caputo IVP is given by a Caputo equation,
with an initial condition, u(0) = u 0 . We will assume that F(t) a continuous function in t, and u 0 ∈ R.
In a similar manner as above we consider an ansatz and look for solutions of the form with u c (t) a continuous function, a a constant, and H(t) as defined in the previous sections.

Theorem 3.
For an Caputo IVP (Definition 6) assume that a solution in the form of the ansatz (Equation (93)) exists. Then such a solution does not possess a step discontinuity at t = 0.
Proof. The proof follows from assuming a solution exists in the form of the ansatz and then showing that the only permitted value of the parameter a is zero. From the ansatz we have with u c (t) being the derivative of a continuous function, δ a Dirac delta, and a some unknown constant.
To obtain the value of the unknown constant we will substitute this into the IVP to give Next we take the limit as t → 0 + , The left-hand side of this expression only exists if a = 0, whilst the right-hand side is well defined. As such solutions with a step discontinuity at the origin do not exist for the Caputo derivative.

Conclusions
We have shown that the solutions of initial value problems using both the CF and ABC operators feature, in almost all cases, a discontinuity at the origin. The occurrence of the discontinuity is problematic for the application of the CF and ABC operators in modelling. Very few physical processes are well described by discontinuous functions, and fewer still with the discontinuity at the origin.
This discontinuity also raises issues with the use of these operators as fractional derivatives. Whilst both operators are generalisations of derivatives, in the sense that as α → 1 − we recover an integer order derivative, the lack of smooth solutions is problematic. Many proponents of the use of these operators claim that the non-singular kernel is desirable, but as we have shown here in order for the solution of the IVP to exist the derivative of the solution must be singular, and thus the integrand is still singular.
In the Appendix A we have considered a more generalized ansatz for the solution of CF operator equations. This generalization does not give solutions to IVPs, but we can find solutions that diverge at the origin.
Traditional numerical methods are based on approximating the solution, and derivatives thereof, by their respective discrete counterparts. In the case where the solution exhibits a discontinuity, these approaches attempt to capture an infinite gradient in the same manner as the gradient of a smooth function. As such, resulting schemes are inadequate for capturing the dynamics of solutions admitted by the present IVPs. By decomposing the solution into discontinuous and continuous parts, traditional methods can be modified to approximate the continuous dynamics only.
In the case of the CF operator, the continuous part of the solution follows from a relatively simple integer order differential equation and the vast literature of methods are available to find efficient high order solutions. The numerical solution of the ABC operator equations is more complicated, as the operator does involve a history dependence.
In this work we have concentrated on the often explored case of 0 < α < 1, although extensions to the case α > 1 are possible. Generalisations of the CF and ABC operators for larger values of α exist. To investigate the occurrence of discontinuities in such systems alternate forms of our ansatz would need to be taken.
We have also shown that this type of discontinuity at the origin can not occur in Caputo derivatives. The ansatz approach that we use is applicable to cases where we have derivatives appearing in integrands, such as the CF, ABC, and Caputo operators. This approach would need further modifications to be applicable to Riemann-Liouville type operators where the derivative occurs outside of the integral.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Higher Order Singularities in the the Solution
We could consider a more generalised ansatz, with a higher order singularities, The solution in this case will then be of the form where d i dt i δ(t) denotes i-th derivative of the Dirac delta δ(t) with corresponding unknown constant a i . Inserting this generalised ansatz into Equation (2), one finds for t > 0. In the limit as t → 0 from above, we have since F c is continous. Again, we make use of Leibniz's rule and differentiate Equation (A3) with respect to t; this gives From Equation (A3), one finds Replacing the integral in Equation (A5) with RHS of Equation (A6) and rearranging yields The integral solution of Equation (A7) is given as Hence the general form of the unbounded solution takes the following form Therefore, we see that the unbounded solution is non-unique with the only condition for the coefficients a i s given by Equation (A4). Note that generally the initial value u(0) = u 0 can not be imposed for solution of this form, as it involves a delta function and its distributional derivatives at the origin, unless we force the unknown coefficients a i = 0 for i ∈ Z + . If we set {a 1 , a 2 , ...} = 0, from Equation (A4), it is then required that This allows us to impose the initial condition u(0) = u 0 , and the solution for the IVP in this case reads From this we see that higher order discontinuities at the origin can still produce solutions to the equation, but do not provide solutions for an IVP.