Separation and stability of solutions to nonlinear systems involving Caputo–Fabrizio derivatives

This work mainly investigates the separation and stability of solutions to nonlinear systems involving Caputo–Fabrizio fractional derivatives. An inequality ensuring the positivity of the fractional derivative at a given point is derived, by which the sufficient conditions for the separation of solutions are obtained. The comparison principle and the inequality for the fractional derivatives of convex functions are obtained, by which the approach of the convex Lyapunov functions is extended effectively to establish the criteria for the stability of solutions in the context of Caputo–Fabrizio fractional derivatives. Applications of the main results are illustrated by using examples.


Introduction
The fractional derivative, in which the order of the derivative is permitted to be noninteger, originated in 1695. Although it is as old as the classical calculus, the world of analysis seems to be always dominated by the latter. It has been thought out of mainstream science for over 300 years, and it is only in the last 40 years that it has come back to life.
It is worth stressing that the recovery of vitality is not just another way of presenting old stories, since the fractional derivatives and fractional differential equations have been applied to numerous fields such as viscoelasticity, signal processing, biology, engineering, neuroscience and materials science in the last four decades, and have become a wellaccepted instrument in the description of complex systems [1][2][3][4][5].
Many important results have been obtained for the investigation of the theory of differential equations with the Riemann-Liouville or Caputo fractional derivatives. For instance, the existence of solutions to the boundary value problems was discussed in [6][7][8][9][10][11][12] for some specific fractional differential or integro-differential equations or inclusions by the use of certain fixed point theorems; and new results were derived in [13] on exact controllability of a class of fractional neutral integro-differential systems. Moreover, the boundary value problems involving the ψ-Caputo fractional derivative were studied in [14].
It is also worth pointing out that the concepts of fractional derivatives are also modified and further developed to adapt to the new requirements of theory and practice. Up to now, there exist several definitions of the fractional derivative in the literature, each being constructed to satisfy various constraints and to be consistent with physical background and experimental data, among which the most popular ones are the Riemann-Liouville and Caputo fractional derivatives. The integral kernels in the previous two kinds of fractional derivatives are singular, which describe the processes involving memory effects.
The retarded effects in many dynamical processes can also be depicted by using integrals with nonsingular kernels, among which the exponential kernels were used in [15,16] to capture the temporal memory. Based on the previous considerations, a new fractional derivative with nonsingular kernel was proposed recently by Caputo and Fabrizio in [17,18], and then the concept of the Atangana-Baleanu derivative with Mittag-Leffler kernel was introduced in [19] and the existence of solutions to some integral equations was discussed in [20,21].
Moreover, the definition of the conformable fractional derivative was also suggested in [22] such that the chain rule holds for the new derivative, and some results for this type of derivative have been already achieved [23][24][25][26].
The Caputo and Fabrizio fractional derivative (CFFD) and the differential equations involving the CFFD have received growing attention. In recent years, a considerable literature has sprung up in the framework of CFFDs. The definition of the CFFD and its properties were deeply explored and further extended in [27][28][29][30][31][32][33][34], among which the definition of the extended CFFD was first introduced in [34]. The well-posedness for the differential systems with the CFFD was well investigated in the published studies [35][36][37][38][39][40][41][42]; among them, the existence of solutions for a coupled system of differential inclusions was discussed in [37] by using some fixed point theorems for the multivalued maps, the existence of solutions and of approximate solutions for some high-order fractional integrodifferential equations were completely considered in [38] and [39][40][41], respectively, by means of fixed point theorems and by virtue of approaches involving α-contractive maps, and the results for the existence and dimension of the set of solutions were obtained in [42] for the second fractional integro-differential inclusion problem with the extended CFFD. And more importantly, the differential equations involving the CFFD have found a wide range of applications [43][44][45][46][47][48][49].
For the above-mentioned equations involving the CFFD, relatively speaking, there are very few results in the published studies for the separation of solutions and for the stability in the Lyapunov sense, and in the context of Caputo factional derivatives, the former has been considered for the case where the state space is one-dimensional [2,[50][51][52]; while the latter has been completely discussed in [53][54][55][56][57][58][59] by using the Lyapunov second's method. It is noted that the direct calculation of the fractional derivative of the Lyapunov functions is far from easy, and the method of the convex Lyapunov functions was recently introduced in [60] and then developed in [61,62], which effectively settles the difficulty.
It is also worth remarking that the Hyers-Ulam stability was discussed in [63][64][65] for the specific integro-differential equations and advection-reaction diffusion system with the Atangana-Baleanu derivative.
Inspired by results and techniques in the aforementioned papers, in this work we will discuss the qualitative properties and stability of solutions to the initial value problem (IVP) with the extend CFFD. Concretely speaking, the first goal is to investigate the separation of trajectories of solutions to the IVP involving CFFDs as follows: where D (α) a denotes the CFFD of order α (defined below); ω and F are vector-valued functions on [a, ∞) and R d+1 , respectively, and F(a, ω a ) = 0. The approach employed here depends on the inequality to be established in this work, which is different from any other methods used to discuss the same problems with the Caputo fractional derivative in the existing studies [2,[50][51][52]. Here the separation of solutions means that the trajectories of solutions with different initial values do not intersect each other.
Based on the previous discussion, under the local Lipschitz condition, we then explore the stability of the zero solution to the system where f is a vector-valued function on R d , and p(t) is a scalar function on [a, ∞) with p(a) = 0. It is also worth remarking that, in the framework of CFFDs and under the local Lipschitz condition, there are no global existence results for the solutions in the published studies due to the restriction of the construction of the integral equations equivalent to the fractional ones under consideration, which causes the difficulty to the discussion of the stability.
To arrive at our second goal, using the comparison principle and the inequality involving the CFFD of convex functions to be established in this work, we follow the idea from [61,66] by extending the vector fields of fractional nonlinear systems and obtain the criteria for the stability of solutions under the local Lipschitz condition. Moreover, we have also to point out that the discussion of stability of autonomous systems (p(t) ≡ 1 in (2)) is trivial (see Remark 4.1 in Sect. 4).
The paper is divided into six parts. The Caputo-Fabrizio fractional derivative and integral are introduced in Sect. 2. We first establish a comparison principle for the CFFDs in Sect. 3, and then obtain an inequality involving CFFDs by which the result for the separation of trajectories of solutions to the IVP (1) for the case of one dimension is derived. We then establish an inequality involving the CFFDs of convex functions in Sect. 4, by which, combined with the comparison principle, we discuss the stability of solutions to the IVP (2). The validity of the main result is illustrated in Sect. 5 by using two examples, and finally, the conclusions are discussed in Sect. 6.

Preliminaries
We introduce the Caputo-Fabrizio fractional derivative and integra in this section, and show their basic properties to be used in what follows.
Given a positive integer d, the d-dimensional Euclidean space is denoted by R d , equipped with the classical norm | · |. Let J a be a left-closed interval with the left endpoint a and Ω be a connected open subset of R d , and denote by C(J a , Ω) the space of all continuous functions from J a into Ω. The gradient vector of a function V on R d is denoted by ∇V .
Moreover, unless otherwise specified, it is always assumed that the exponent α is in the interval (0, 1) throughout this work.
Given an absolutely continuous function ω in the interval [a, b], the Caputo-Fabrizio fractional derivative of order α is defined [17,18] Here the function e α (·) is defined by Using the results obtained by integration by parts on the right-hand side of (3), we next give the weak version of the CFFD, denoted by D (α) a , without requirement for the smoothness of functions.

Definition 2.2
Letting ω be in C(J a , R), the Caputo-Fabrizio fractional integral of order α is defined by Remark 2.2 Here our definition of the fractional integral is slightly different from those in [27,33], by which a harmonic connection between D (α) a I (α) a and I (α) a D (α) a is established.
Proof Letting β = α 1-α and using Remark 2.1 and Definition 2.2, a direct calculation yields and The assertions are thus verified.

Comparison principles and separation of trajectories
The following assertion can be directly verified by using Definition 2.2.
Using Lemma 3.1, we next derive a comparison principle, which will play a key role in the discussion of the stability.

Lemma 3.2 Let ω and υ belong to
Proof The continuities of ω and υ ensure that both D (α) a ω(t) and a ω(a) = 0 = D (α) a υ(a) and using the inequality for CFFD in the conditions, we infer by Lemma 3.1 that and hence by Lemma 2.1 that which is equivalent to which due to the condition ω(a) ≤ υ(a) yields the desired result.
An important inequality is next established to be used to discuss the separation of trajectories of solutions to IVPs.

Lemma 3.3 Let a function χ belong to
Proof From Definition 2.1 and the condition χ(t 1 ) = 0, it follows that This, together with the assumption χ(t) < 0 on [a, t 1 ), verifies the desired assertion.

Remark 3.1 The preceding inequality is strict, which plays a central role in the investigation on the separation of trajectories.
Before starting with the discussion of the separation of trajectories, we here have to consider the existence of solutions.
Lemma 3.4 is easily verified by using Definitions 2.1 and 2.2 and Lemma 2.1. (1) is replaced by that in (3), then the equivalence in Lemma 3.4 does not hold generally any more.

Remark 3.2 If the fractional derivative in problem
And then, using the Lipschitz condition, we obtain the following estimate: Multiplying both sides of the above inequality by exp(-μt) yields which implies that A is a contraction because of L(2(1α) + αμ -1 ) < 1. It follows from the Banach fixed point theorem that the operator has a unique fixed point in C([a, T], R d ), and Lemma 3.4 ensures that this fixed point is the unique solution to problem (1). Finally, the solution is clearly defined on R + due to arbitrariness of choice of T. The proof is complete.

Convex Lyapunov functions and stability
We explore in this section the stability of the solutions to the IVP (2). (2) is stable if, for any > 0, there is a positive number δ such that if |ω a | < δ, then, for the solution ω(t, ω a ) of (2), we have |ω(t, ω a )| < for all t ≥ a. To see this, let p(t) ≡ 1 in (2), and then the system (2) becomes autonomous. Let ω be in C(J a , R d ), then D (α)

Definition 4.1 The zero solution of the equation in
a ω(t) belongs also to C(J a , R d ) with D (α) t ω(a) = 0, which forces f (ω a ) = 0, and thus ω(t) ≡ ω a is the unique solution to (2). Accordingly, in order to discuss the stability of system (2) with p(t) ≡ 1, it must hold that f (ω) ≡ 0 whenever |ω| < δ for some δ > 0, which obviously is trivial.

Definition 4.2 ([67])
A continuous function κ : R + → R + is said to belong to class K if it is strictly increasing and κ(0) = 0. For such a function, we shall usually write κ ∈ K.
Lemma 4.1 (Kirszbraun's theorem, [68]) If S ⊂ R m and f : S → R n is Lipschitz, then f has a Lipschitz extension f E : Given a positive number r, define the subsets S(r) of R d by We now are in a position to present an important inequality involving the CFFD of composition function in this work.

Lemma 4.2 Let V : S(r) → R be a continuously differentiable and convex function. Then for any function ω in C(J a , S(r)), the inequality
holds on J a . Here and in what follows, the dot · denotes an inner product on R d .
Using Definition 2.1, by the inequality just derived, we deduce the following estimate: which verifies the desired inequality.
In order to discuss the stability of solutions to problem (2), let S(r) be defined as in Sect. 4 and we further impose several assumptions on the functions in (2) We now present a result for the stability of solutions by using convex Lyapunov functions. Proof Given a positive number r, choose an arbitrary positive number ε in the interval (0, r). And then it follows from the continuity of V (ω) and V (0) = 0 that there exists a positive number δ < ε such that for any ω a in S(δ), Using condition (A2), we infer by Lemma 4.1 that f has a Lipschitz extension f E : and f E = f on S(r), and using the function f E just obtained and the initial data ω a previously chosen, construct an IVP as follows: and then it follows from Theorem 3.1 that there is a unique solution to the IVP (7), denoted byω(t, ω a ) and defined on R + , and we further claim that it is also the unique solution to the IVP (2) defined on R + . To this end, we first verify that the inequality holds on R + . Now suppose that, contrary to our assertion, the inequality is false, and then, from the condition |ω(a, ω a )| = |ω a | < δ < ε and the continuity of the solutionω(t, ω a ), it follows that there exists a positive number t 1 > 0 such that for each t in the interval [a, t 1 ), with |ω(t 1 , ω a )| = ε, and observing ε < r, it follows from the inequality in (9) and the con- Also, by the inequality in (9) and the second inequality in the conditions, we deduce that for each t in the interval [a, t 1 ), Due to the differentiability and convexity of the function V , Lemma 4.2 and the preceding inequality imply that for each t in the interval [a, t 1 ), Combined with the inequalities in (8)- (11), the chain of inequalities are derived as follows: and for t = t 1 , in particular, and this obvious contradiction verifies the previous assertion that the inequality in (8) holds on R + . The inequality in (8) obviously implies that for each t in R + ,ω(t, ω a ) lies in S(r) from which we get f E ω(t, ω a ) = f ω(t, ω a ) , and thus Therefore, the functionω(t, ω a ) is also a unique solution to the IVP (1) defined on R + .
Again, using the inequality in (8), the stability is consequently verified.

Conclusions
The CFFD was introduced in 2015. Up to now, many results for the theory of differential equations with the CFFD have been obtained, and these equations have also found a large spectrum of applications. But very little has been published on the subject of the separation and Lyapunov stability of solutions to nonlinear systems involving CFFD.
In this work, using the techniques of inequalities we derived a sufficient condition for the separation of solutions to the specific nonlinear equation, and by establishing the inequality for the convex functions, we successfully generalized the method of convex Lyapunov functions to the context of the CFFD and set up the criteria for the stability in the Lyapunov sense for the nonlinear system with the CFFD.
We believe that the inequalities established in this work are also significant on their own. It is expected that the inequality in Lemma 3.3 might be used to further explore the comparison principles, and that in Lemma 4.2 might be utilized to investigate maximum principles for some partial differential equations involving the CFFD.
It is easy to check the validity of the results obtained in this work, and, moreover, it remains to further explore whether or not the results and approaches here can be extended to the systems involving the Mittag-Leffler kernel.