Boundedness and power-type decay of solutions for a class of generalized fractional Langevin equations

In this paper we study the long-time behavior of solutions for a general class of Langevin-type fractional integro-differential equations. The involved fractional derivatives are either of Riemann–Liouville or Caputo type. Reasonable sufﬁcient conditions under which the solutions are bounded or decay like power functions are established. For this purpose, we combine and generalize some well-known integral inequalities with some crucial estimates. Our ﬁndings are supported by examples and special cases.


Introduction
We consider the following class of nonlinear fractional integro-differential equations: supplied with appropriate initial data. The derivatives D α 0 + , D Clearly, many treated fractional integro-differential equations in the literature become special cases of this general class. It is well-known that x(t) = a 1 and x(t) = a 2 t α−1 , 0 < α < 1, are solutions for C D α 0 + x (t) = 0 and D α 0 + x (t) = 0, respectively, where a 1 and a 2 are constants depending on the initial values of x. Under some sufficient conditions, we prove that the solutions of (1) are bounded by a constant when the fractional derivatives are of Caputo type, C D α 0 + , and decay like t α−1 when the derivatives are of Riemann-Liouville type, D α 0 + . We use fractional calculus tools to treat the terms involving fractional integrals and derivatives and to find appropriate underlying spaces for the solutions. An argument is borrowed from Medved [26] to deal with the singularities inside the fractional derivatives. Also, to treat the associated Volterra integral equations corresponding to (1), we use and generalize some known nonlinear integral inequalities.
From both the theoretical point of view and the application point of view, it is of great importance to have an idea about the behavior of solutions for large values of the time variable. It is known from the definitions of fractional derivatives that all the history of the state is taken into account through a convolution with a singular kernel. Moreover, in our case, the nonlinear term may involve additional singularities. Because of all these features, it is difficult to apply the existing approaches and methods in the literature for integer order to the noninteger order case.
Langevin equation is a widely used model to describe the evolution of physical phenomena in fluctuating environments [10]. Several fractional generalizations of Langevin equations have been suggested to describe dynamical systems in fractal media, see e.g. [5,9,24,25,36,38,39] and the references therein. There is a large volume of literature on existence of solutions for various classes of fractional differential and integro-differential equations (see for instance [1][2][3][4]8,12,17,18,20,35,37] ). Agarwal et al. surveyed in [1] many of these existence results. They focused on initial and boundary value problems for fractional differential equations with Caputo fractional derivatives of orders between 0 and 1 and between 1 and 2. Furati and Tatar considered in [13], the Cauchy-type fractional differential problem: where D α 0 + is the Riemann-Liouville fractional derivative of order α, 0 < α < 1. They used the Schauder fixed point theorem to prove a local existence result in the space C 1−α [0, ∞) (see (9), for some classes of the nonlinearities f and k involving some power functions of t, s and x. The case k ≡ 0 has been studied by Delbosco and Rodino [11], Kilbas et al. [21], and many others. In 2012, Trujillo et al. [3] established the existence and uniqueness of solutions for the nonlinear fractional integro-differential problem: Their main tool is a fixed point theorem for non-self mappings. First, using a suitable substitution, they constructed an equivalent fractional integral equation. Then, they used some fractional integral inequalities and a nonlinear alternative of Leray-Schauder type to achieve their existence result. The uniqueness of the solution was established with the help of the Banach contraction principle.
The long-time behavior of solutions of differential equations has attracted many researchers, see [6,7,14,16,19,27,30,31] . In many cases, the main idea is to establish sufficient reasonable conditions ensuring comparison or similarity with the long-time behavior of solutions of simpler differential equations.
In 2004, Momani et al. [30] discussed the Lyapunov stability and asymptotic stability for solutions of the fractional integro-differential equation (2) for t ≥ a > 0. The assumptions: where γ (t) and δ(t) are continuous nonnegative functions and were imposed. The authors proved that every solution x(t) of (2) satisfies |x(t)| ≤ C 0 (t − a) α−1 where C 0 is a positive constant, and hence the solution of (2) is asymptotically stable. Furati and Tatar [14] showed that the solutions of (2) decay polynomially for some nonlinear functions f and k. When k ≡ 0, they proved in [15] that solutions of the problem exist globally and decay as a power function in the space C α 1−α [0, ∞) defined in (26). The same authors considered in [16] Eq. (2) and found bounds for solutions on infinite time intervals and also provided sufficient conditions assuring decay of power type for the solutions.
In [31], Mustafa and Bȃleanu proved that the solution of the nonlinear fractional differential problem: is asymptotic to o (t aα ) as t → ∞, 0 < 1 − a < α. They assumed that, Bȃleanu et al. discussed in [7] the long-time behavior of solutions to the linear fractional differential problem: where D α 0 + is the Riemann-Liouville fractional derivative of order α. They proved that it has a solution x ∈ C ((0, ∞), R) such that, lim and it has the asymptotic expansion, In 2011, the authors in [6] established under the condition: can be expressed asymptotically as In 2015, Medved and Pospíšsil considered in the paper [27] the case when the right-hand side depends on a Caputo fractional derivative of the solution. They proved that there exists a constant c ∈ R such that any global solution of the initial value problem, where 0 < β < α < 1, is asymptotic to ct β . Also, they considered the case when the right hand side depends on ordinary derivatives up to order n − 1 and Caputo derivatives of fractional orders n − 1 < α j < α < n, j = 1, 2, ..., m, m ∈ N. It has been shown that the solution in this case is asymptotic to bt r with r = max {n − 1, α m }.
In 2016, Kassim et al. studied in [19] the boundedness and asymptotic behavior of solutions for the fractional differential problem: They showed that there exists a positive constant A such that, λ > 0 and ψ 1 , ψ 2 are nondecreasing. In this paper, we consider a much wider class of nonlinearities than the ones considered in the previous papers [14,16,30]. Also, we improve the results in these papers by weakening the imposed conditions. Namely, we consider the problem: with appropriate initial conditions depending on D being the Caputo or the Riemann-Liouville fractional derivative. Clearly, this nonlinearity is much more general than the one in [14,30] and even more general than the ones in [6,7,15,31]. The rest of this paper is organized as follows: In the next section we briefly present the used notation, underlying function spaces, background material and some useful lemmas and inequalities. Section 3 is devoted to the main results of the study of the long-time behavior of solutions of (1).

Preliminaries
In this section, we briefly introduce some basic definitions, notions and properties from the theory of fractional differential equations.
and ess sup | f (t)| is the essential supremum of the function | f (t)|.

Definition 2.3 [22]
We denote by C γ [a, b] , 0 ≤ γ < 1, the following weighted space of continuous functions where Definition 2.4 [22] For n ∈ N and 0 ≤ γ < 1, we denote by C n γ [a, b], the following weighted space of continuously differentiable functions up to order n In

Definition 2.5 The Riemann-Liouville left-sided fractional integral of order α > 0 is defined by
provided the right-hand side exists. We define I 0 The Riemann-Liouville left-sided fractional derivative of order α ≥ 0, is defined by: is the integral part of α. In particular, when α = m ∈ N 0 , it follows from the definition that, D m a + u = D m u. Definition 2.7 The Caputo left-sided fractional derivative of order α ≥ 0, is defined by: where n = [α] + 1 for α / ∈ N 0 and n = α for α ∈ N 0 . In particular, when α = n ∈ N 0 , it follows from the definition that C D 0 a + u = u, C D n a + u = D n u. The next lemma shows that the Riemann-Liouville fractional integral and derivative of the power functions yield power functions multiplied by certain coefficients and with the order of the fractional derivative added or subtracted from the power. Lemma 2.8 [22] If α ≥ 0, β > 0, then The following result is about the composition I α a+ D α a + of the Riemann-Liouville fractional integration and differentiation operators.
The next lemma is an analog of Lemma 2.11 for the Caputo fractional derivative.
for all t ∈ (a, b].
For more details about fractional integrals and fractional derivatives, the reader is referred to the books [22,29,33,34]. We mention here some useful basic inequalities and some linear and nonlinear integral inequalities to be used in the next sections. Lemma 2.13 [28] If λ, ν, ω > 0, then, for any t > 0, we have where C is a positive constant independent of t. In fact, The inequalities in the next lemma are of Jensen's inequalities type.

The constants c i are given by c 0 = c and c
, i = 1, 2, 3.

Remark 2.17
The monotonicity and the ordering requirements of Lemma 2.16, can be dropped using the functions: Note that

Main results
In this section, the long-time behavior of solutions is investigated for the initial value problem composed of Eq. (1) with the initial conditions x(0) = c 0 and I 1−α 0 + x(0 + ) = c 1 , where c 0 , c 1 ∈ R, when the fractional derivatives are of Caputo and Riemann-Liouville types, respectively. The first section is dedicated to the study of the long-time behavior of the solutions for (1) with Caputo fractional derivatives. The case of the Riemann-Liouville derivatives is discussed in Sect. 3.2.
Before presenting our results we need to define the following classes of functions:

Definition 3.3
We say that a function g is of type G if it is continuous nondecreasing on [0, ∞) and positive on (0, ∞).
The above classes are not empty. Examples showing this fact are given in the next sections.

Equations with Caputo fractional derivatives
Equation (1) with the Caputo fractional derivatives of orders 0 ≤ β, γ < α ≤ 1, is considered in this section, that is, We discuss the boundedness of the continuable solutions of (16) in the space C α [0, ∞) defined by: By "continuable", we mean that the solution as it exists locally we may prolongate it to infinity. It will exist globally in time.
We suppose that the functions f and k satisfy the following hypotheses: The main result of this section is stated in the following theorem.

Theorem 3.4 Suppose that the functions f and k satisfy (H 1 ), (H 2 ) and (H 3 ). Then, there exists a positive
constant c ∈ R such that any solution x ∈ C α [0, ∞) of problem (16) satisfies: Proof Applying I α 0 + to both sides of the equation in (16), we obtain from Lemma 2.12, for all t > 0. Taking the derivatives C D β 0 + and C D γ 0 + of (17), gives in view of Lemmas 2.8 and 2.10, By virtue of (H 3 ), we have: Similarly, we get for all t > 0, Using Hölder's inequality, we find: Since p(α − 1) + 1 > 0 and pr + 1 > 0, we see, from Lemma 2.8, that Recalling that α + r − 1 + 1 p = 0, so (18) becomes, for all t > 0. Similarly, we get the following estimates on C D β 0 + x and C D γ 0 + x for all t > 0, where for all t > 0, where It is obvious from the inequalities (19)-(22) that, As g 1 , g 2 and g 3 are nondecreasing functions, it is clear that, Using Lemma 2.14, we obtain: where B 5 = 3 q−1 |c 0 | q and B 6 = 3 q−1 B q 4 . Let u(t) = (z(t)) q . Then we have for all t > 0, From Lemma 2.16, with we deduce that, where M is a positive constant. The desired result follows in virtue of (23).

Remark 3.5
As special case of Theorem 3.4 with β = γ = 0, there exists a positive constant c ∈ R such that any continuable solution x ∈ C α [0, ∞) of the problem: satisfies |x(t)| ≤ c, for all t > 0.

Remark 3.6
It is not hard to see that the conclusion of Theorem 3.4 is still valid if we replace the condition (H 3 ) by the following condition: There are functions k 1 of type H αq−1,q , k 2 of type H 0,1 , q ≥ 1 and f i of Example 3.7 Consider the fractional integro-differential equation, where 0 < α ≤ 1, μ 1 , μ 2 , μ 3 > 0 and 0 ≤ σ 1 ≤ σ 2 q ≤ σ 3 q ≤ 1 . Let All the hypotheses of Theorem 3.4 are fulfilled,

Conclusion
In this study, we considered some classes of fractional integro-differential equations with the Caputo or Riemann-Liouville derivatives in both sides of the equations depending on the solution, its fractional derivatives as well as an integral of a kernel involving the solution or its fractional derivatives. We assumed the continuity of the nonlinearities and the boundedness of these nonlinearities by sums or products of continuous functions of time, in certain Lebesgue spaces, and nondecreasing functions of the states. Under these nonlinear growth conditions on the nonlinearities, we treated initial value problems for which, in general, solutions cannot be found explicitly. We found that their solutions, under these conditions, behave like the solutions of the associated linear fractional differential equations with zero right-hand sides.