Stability Concepts of Riemann-Liouville Fractional-Order Delay Nonlinear Systems

: First, we set up in an appropriate way the initial value problem for nonlinear delay differential equations with a Riemann-Liouville (RL) fractional derivative. We deﬁne stability in time and generalize Mittag-Lefﬂer stability for RL fractional differential equations and we study stability properties by an appropriate modiﬁcation of the Razumikhin method. Two different types of derivatives of Lyapunov functions are studied: the RL fractional derivative when the argument of the Lyapunov function is any solution of the studied problem and a special type of Dini fractional derivative among the studied problem.


Introduction
Various processes with anomalous dynamics in science and engineering can be formulated mathematically using fractional differential operators ( [1][2][3]). When the Riemann-Liouville (RL) fractional derivative is applied in differential equations, the statement of initial conditions is important. It is worth mentioning that the physical and geometric interpretations of operations of fractional integration and differentiation were suggested by Podlubny [4]. Fractional differential equations in terms of the RL derivative require initial conditions expressed in terms of initial values of fractional derivatives of the unknown function ( [5,6]). In [7], it was shown that the initial conditions for fractional differential equations with RL derivatives expressed in terms of fractional derivatives has physical meaning. In fact, it was shown that for any physically realistic model, zero initial conditions will be found for a continuous loading program or even in the case of a step discontinuity. Nonzero conditions will only be found in the case of an impulse, and this type of process can be found in physics, chemistry, engineering, biology and economics. In the case of zero initial conditions, the RL, Grünwald-Letnikov (GL) and Caputo fractional derivatives coincide ( [4]). For this reason, some authors either study Caputo derivatives or use RL derivatives but avoid the problem of initial values of fractional derivatives by treating only the case of zero initial conditions. This leads to the consideration of mathematical correct problems, but without taking the physical nature of the described process into account. Sometimes, such as in the case of impulse response, nonzero initial conditions appear (see, for example, [7]).
In connection with the main idea of stability properties, we will consider in this paper nonzero initial conditions for RL fractional equations, and we will define in an appropriate way stability properties that are slightly different than those for Caputo fractional differential equations.
Note that stability properties of delay differential equations can be considered by an application of the Lyapunov-Krasovskii method by functionals or by the Razumikhin method by Lyapunov functions. It is worth mentioning that both mentioned methods are applied for stability study of Caputo fractional delay differential equations (see, for example, [8][9][10][11]).
In the case of delay fractional differential equations with the RL fractional derivative, following the idea of initial conditions in ordinary delay differential equations and the above-mentioned idea concerning the initial condition for RL fractional differential equations without any delay, we will set up initial conditions in an appropriate way. Note that any solution of the defined initial conditions with RL fractional derivatives is not continuous at zero (the initial point), which is the same as in the case without any delay. Delay RL fractional differential equations are set up and studied in [12], but the initial condition does not correspond to the idea of the case of delay differential equations with ordinary derivatives (the lower bound of the RL fractional derivative coincides with the left end side of the initial interval).
Asymptotic stability for RL fractional differential equations with delays is studied in [13][14][15], but only the autonomous case is considered. A Lyapunov functional and its integer order derivative is applied. This functional is similar to the one used in the theory of differential equations with ordinary derivative and delay. On one hand, the application of ordinary derivative of the Lyapunov functional is not similar to the used fractional derivatives in the equation; on the other hand, it leads to some restrictions on both the delay and the right side parts of the equation ( [16]). Additionally, in [15], the initial condition is not adequately associated with the RL fractional derivative. RL fractional equations with delays were studied recently in [17,18], but there are unclear parts in the statement of the problem (the lower limit of the RL derivative is different than the initial time point) as well as in the initial condition (the RL fractional integral has no meaning, compared with [5] at the initial time). The Razumikhin method is applied to RL fractional differential equations in [11], but the initial condition is not connected with the RL fractional derivative.
In this paper, the initial value problem for nonlinear delay differential equations with the RL fractional derivative is studied. Based on the arguments in the books [5,6], we set up initial conditions expressed in terms of initial values of fractional derivatives of the unknown function. Any solution of the defined initial conditions with RL fractional derivatives is not continuous at zero (the initial point). We require a new definition for stability excluding a small interval around zero. We define stability in time and generalize Mittag-Leffler stability in time for RL fractional differential equations. The stability properties of the zero solution are studied by Lyapunov functions. An appropriate modification of the Razumikhin method is suggested. Two types of derivatives of Lyapunov functions are applied: the RL fractional derivative when the argument of the Lyapunov function is a solution of the studied problem and the Dini fractional derivative among the studied problem.
The main contribution in the paper could be summarized as follows: the initial conditions connected with the RL fractional derivative are set up in an appropriate way; -new types of stability connected with the type of initial conditions are defined; -the RL fractional modification of the Razumikhin method is presented; -new sufficient conditions for the defined stability are obtained; -two types of fractional derivatives of the Lyapunov functions are used.

Preliminary Notes
In this paper, we will use the following definitions that are well known in the literature ( [5,19]): . (1) The fractional derivatives for scalar functions could be easily generalized to the vector case, by taking fractional derivatives with the same fractional order for all components.
Define the class We will provide some well-known results from the literature (see, for example, [5,19]):

Remark 2.
As it is mentioned in [20] (see Example 1 [20]), a function might not be differentiable at one point in the classical sense, but it is RL differentiable. The positive RL fractional derivative RL 0 D q t m(t) > 0 of order q ∈ (0, 1) only means that the RL fractional integral RL 0 I is monotonously increasing with respect to t and it does not imply that the function m(t) is monotonously increasing. So, we cannot regard RL and Caputo derivatives as the generalization of the ordinary derivative in a rigorous mathematical way. Proposition 4 (Lemma 2.3 [21]). Let m ∈ C 1−q ([0, T], R). Suppose that for any t 1 ∈ (0, T], we have m(t 1 ) = 0 and m(t) < 0 for 0 ≤ t < t 1 . Then, it follows that RL 0 D q t m(t)| t=t 1 ≥ 0.
has a unique solution is the two-parameter Mittag-Leffler function.
Now, based on Proposition 5, we will illustrate the importance of the initial condition when the RL fractional derivative is used in the equation. For simplicity, we will consider the case of equations without any delays.

Remark 3. Consider the scalar linear RL fractional equation
It is well known that RL 0 D q t t q−1 E q,q (t q ) = t q−1 E q,q (t q ), i.e., the solution of the above RL fractional differential equation, is x(t) = ct q−1 E q,q (t q ) where c is a real constant. Now, consider the initial condition x(0) = k where k is a real constant. However, t q−1 E q,q (t q )| t=0+ = ∞ = k. This illustrates that the initial condition of the type x(0) = k is not applicable for RL fractional equations (see, for example, [11]). Now, consider the initial condition t 1−q x(t)| t=0+ = k where k is a real constant. Then, i.e., the initial condition t 1−q x(t)| t=0+ = k has a meaning for the RL fractional derivative with k = c Γ(q) .
The practical definition of the initial condition of fractional differential equations with RL derivatives is based on the following result: (b) If there exists a.e. a limit lim t→0+ 0 I 1−q t m(t) = c ∈ R, and if there exists the limit

Statement of the Problem
Consider the following nonlinear Riemann-Liouville fractional delay differential equation (RLFrDDE) of fractional order q ∈ (0, 1): with initial conditions where We denote the solution of the initial value problem (IVP) for RLFrDDE (2) and (3) by In this paper, we will assume that the function f is such that for any continuous initial function φ the IVP for RLFrDDE (2) and (3) has a solution. Note that some existence and uniqueness results to RL fractional differential equations with delay were obtained in [22].
We will give the basic definitions for stability: (2) and (3) (with the zero initial function) is said to be stable in time if for any number ε > 0 there exist numbers δ > 0 and T ε > 0 depending on ε such that for any initial functions φ ∈ C 0 : ||φ|| 0 < δ, the corresponding solution x(t; φ) of IVP (2) and (3)

Definition 1. The zero solution of RLFrDDE
As an example to discuss stability in time, we consider a scalar RL fractional equation without any delay, whose exact solution is known. Example 1. According to Proposition 5, the scalar initial value problem Therefore, the zero solution is generalized Mittag-Leffler stable with γ = 1 − q, m(u) = u and λ = a.
The zero solution is stable in time because for any ε > 0 there exist δ, T ε > 0 such that At the same time, the zero solution is not stable (in the regular sense) because lim t→0 t q−1 E q,q (−at q ) = ∞ and |t q−1 E q,q (−at q )| < ε cannot be satisfied for values of t close to 0.
For example, if a = 1, Throughout the paper, we shall use the class

Lyapunov Functions and Their Derivatives
One approach to study stability properties of nonlinear RL fractional differential equations is based on the application of Lyapunov functions and an appropriate modification of the Razumikhin method. The first step is to define a Lyapunov function. The second step is to define its derivative among the studied equation.
We will use the following class of functions called Lyapunov functions: , T ≤ ∞, be a given interval, and ∆ ⊂ R n , 0 ∈ ∆ be a given set. We will say that the function V(t, x) : J × ∆ → R + belongs to the class Λ(J, ∆) if V(t, x) is continuous on J/{0} × ∆, and it is locally Lipschitzian with respect to its second argument.
In our study, we will use the Razumikhin condition for the Lyapunov function V ∈ Λ(J, ∆) and any ψ ∈ C([−τ, 0], R n ): We will give a brief overview of the derivatives of Lyapunov functions among solutions of fractional differential equations in the literature. There are three main types of derivatives of Lyapunov functions from the class Λ(J, ∆) used in the literature to study stability properties of solutions of fractional differential Equation (2): (2) and (3). Then, we consider Then, consider (see [8]) where q C r is defined by (1) and The Dini fractional derivative is applicable for continuous Lyapunov functions.
We will study stability properties of the zero solution of RLFrDDE (2) by an application of both defined types of fractional derivatives of Lyapunov functions.

Stability by the RL Fractional Derivative of Lyapunov Functions
We will obtain some sufficient conditions for stability with applications of the RL fractional derivative of Lyapunov functions.
Then, the zero solution of (2) and (3) with the zero initial function is stable in time.
Proof. Let > 0 be an arbitrary number. According to condition (i), there existsT ε > 0 such that inequality (9) holds for t >T ε , x ∈ R n .
From inequality (13) and condition (i), it follows that where T ε = max{T ε ,T ε }. This proves the stability in time of the zero solution.

Remark 7.
The main condition in Theorem 1 is condition (iii), which is connected with any solution of the IVP (2) and (3).

Remark 8.
According to Remark 6, the Lyapunov function V(t, y Remark 9. Condition (i) is different than the corresponding condition for Lyapunov functions for differential equations with ordinary derivatives as well as with Caputo fractional differential equations. This is because of the type of initial condition (see Remark 3).
holds where c > 0, x(t) is a solution of the IVP for RLFrDDE (2) and (3). Then, the zero solution of (2) and (3) with the zero initial function is generalized Mittag-Leffler stable.

Corollary 1.
If all the conditions of Theorem 2 are satisfied, then the zero solution of (2) and (3) is asymptotically stable.

Stability by the Dini Fractional Derivative of Lyapunov Functions
We will study stability by the application of the defined above Dini fractional derivative of Lyapunov functions among the studied delay fractional differential equations.
Initially, we will prove a comparison result for Lyapunov functions.
Lemma 2. Assume: is a solution of the IVP for RLFrDDE (2) and (3) with 2. The function V ∈ Λ([−τ, T], ∆), ∆ ⊂ R n , is such that : (i) There exists an increasing function g ∈ C(R + , R + ) such that the inequality is the Dini fractional derivative defined by (5) and Remark 10. Let us, for simplicity, again consider the RL fractional differential equation without any delay RL 0 D q t x(t) = −x(t) for t > 0, RL 0 I 1−q t x(t)| t=0+ = c with a solution x(t) = ct q−1 E q,q (−t q ) (see Example 1). If we consider the quadratic Lyapunov function V(x) = x 2 , then D + (2) V(t, ψ(0), ψ) ≤ 0 and (ct q−1 E q,q (−t q )) 2 ≤ c 2 is not satisfied. However, if V(x) = t 1−q x 2 then D + (2) V(t, ψ(0), ψ) ≤ 0 and t 1−q (ct q−1 E q,q (−t q )) 2 ≤ c 2 t q−1 is satisfied. This example again illustrates the changes in the applied Lyapunov functions and their conditions in the application of RL fractional derivatives comparatively with the application in Caputo fractional derivatives.
Then, the zero solution of (2) with the zero initial function is stable.
This proves the stability in time of the zero solution.

Conclusions
The nonlinear RL fractional differential equation is studied. The initial value problem is a subject that remains quite up-to-date (see, for example, the books [5,6]). Note the initial condition imposed to study fractional kinetic equations with RL fractional derivative. This point is critical in many physical situations, especially in astrophysical problems and the problem of anomalous subdiffusion ( [23]). A good overview of the physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives was most clearly formulated by Diethelm [19] and it is detailed discussed in [7], where it is shown that initial conditions for RL fractional differential equations have physical meaning, and that the corresponding quantities can be obtained from measurements.
In this paper, some new definition for stability excluding a small interval around zero is defined and studied. These types of stability are called stability in time and generalize Mittag-Leffler stability in time for RL fractional differential equations. The definitions are deeply connected with the singularity at the initial time point. The stability properties of the zero solution are studied by Lyapunov functions. Two types of derivatives of Lyapunov functions: the RL fractional derivative when the argument of the Lyapunov function is a solution of the studied problem and the Dini fractional derivative among the studied problem.