Existence, Stability and Simulation of a Class of Nonlinear Fractional Langevin Equations Involving Nonsingular Mittag–Lefﬂer Kernel

: The fractional Langevin equation is a very effective mathematical model for depicting the random motion of particles in complex viscous elastic liquids. This manuscript is mainly concerned with a class of nonlinear fractional Langevin equations involving nonsingular Mittag–Lefﬂer (ML) kernel. We ﬁrst investigate the existence and uniqueness of the solution by employing some ﬁxed-point theorems. Then, we apply direct analysis to obtain the Ulam–Hyers (UH) type stability. Finally, the theoretical analysis and numerical simulation of some interesting examples show that there is a great difference between the fractional Langevin equation and integer Langevin equation in describing the random motion of free particles.


Introduction
To expound the random motion of particles in fluid after colliding with each other, Langevin raised the famous Langevin equation in 1908.Afterward, many random phenomena and processes were found to be described by the Langevin Equation [1,2].However, the integer-order Langevin equation is unable to meet the accuracy requirements in describing complex viscoelasticity.Thereby, the classical Langevin equation has been extended and modified.Kubo [3,4] put forward a general Langevin equation to simulate the complex viscoelastic anomalous diffusion process.Eab and Lim [5] applied a fractional Langevin equation to describe single-file diffusion.Sandev and Tomovski [6] established a fractional Langevin equation model to study the motion of free particles driven by power-law noise.Furthermore, the stability of the system represents the most important dynamics characteristic.Ulam and Hyers [7,8] proposed a concept of system stability called UH-stability in the 1940s.Over the past decade, there have been many works published (some of which can be found in [9][10][11][12][13][14][15]) on the UH-stability of a fractional system.
It is worth noting that these works on the fractional Langevin system basically involve Caputo or Riemann-Liouville fractional derivatives.In fact, the Caputo or Riemann-Liouville fractional derivatives can produce singularity under some conditions.This makes them difficult to employ as mathematical models of certain physical phenomena.Consequently, a new nonsingular fractional derivative with exponential kernel was raised by Caputo and Febrizio in [16].Furthermore, another new nonsingular fractional derivative with ML-kernel was put forward by Atangana and Baleanu in [17].Since their introduction, these nonsingular fractional derivatives have attracted much attention and research in theory [18][19][20][21] and application [22][23][24][25][26][27].Some new findings on the fractional Langevin equation have been published in recent papers (see [28][29][30][31][32][33][34][35][36]).However, there are a paucity of papers on Ulam-Hyers stability of fractional Langevin system with ML-kernel.
The remaining structure of the manuscript is as follows.Section 2 introduces some fundamental definitions and lemmas.In Section 3, we obtain some criteria on the existence of solutions to the system (1) by using some fixed-point theorems.The UH-type stabilities of (1) are built in Section 4. As applications, we conduct theoretical analysis and numerical simulation on some examples to verify the correctness and effectiveness of our main results in Section 5. Finally, a brief summary is provided in Section 6.
Remark 1.The Caputo fractional derivative of order 0 < α ≤ 1 of a continuous function u : (0, ∞)→R is defined by provided that the right-hand side is pointwise defined on (0, ∞).From Definition 2 and the Caputo fractional derivative, one finds two differences between them.One is that the coefficients are different.The other is that the kernel function is different.The kernel function (t − s) −α of Caputo fractional derivative is singular at s = t, but the kernel function E α − α 1−α (t − s) α of Mittag-Leffler fractional derivative is nonsingular at s = t.Lemma 1. [37] Assume that h ∈ C[0, T].Then, the unique solution of fractional differential equation Remark 2. It follows from Definition 2 and Lemma 1 that ML D α 0 + u(t) ≡ 0 if and only if u(t) ≡ constant.
N (α) = 0, then the fractional differential Langevin Equation (1) is equivalent to the following integral equation Proof.Assume that the function u(t) ∈ C(0, T) is a solution of (1), Then, for t ∈ D, we derive from Lemma 1 that (3) gives From Lemma 1, u(0) = A and (4), we yield (5) Noting the last integral term of (5), we exchange the order of double integrals to get It follows from ( 5) and (6) that Thus, Equation (2) holds, that is, u(t) ∈ C(0, T) is also a solution of integral Equation (2).Furthermore, vice versa, if u(t) ∈ C(0, T) is a solution of integral Equation (2), then one knows that (4) and (3) hold by finding the fractional derivative Ml D α 0 + at both ends of (2).Next, by finding the fractional derivative Ml D β 0 + at both ends of (3), one easily gets the first fractional equation of (1).According to ( 2) and ( 3), we have u(0) = A and Ml D α 0 + u(0) = B. Thus, we verify that u(t) ∈ C(0, T) also satisfies system (1).The proof is completed.

Existence of Solutions
In this section, by applying the following important fixed-point theorems, we emphasize to investigate the existence of solutions for system (1).Lemma 3. [38] Let E be a non-empty closed subset of a Banach space X.
Lemma 4. (Krasnoselskii's fixed-point theorem [39]) Let Y be a non-empty closed convex subset of a Banach space X.Assume that P and Q are two operators satisfying (ii) P is contraction, and Q is continuous and compact.
Then there is at least a solution u * ∈ Y such that u * = Pu * + Qu * .By (2), we take D = [0, T], X = C(D, R).Then X is a Banach space with the norm u = sup t∈D |u(t)|.We shall study the existence and stability of the solution of (1) in (X, • ).In the whole paper, we need the following essential assumption.
(H 1 ) T, α, β, λ, A and B are some constants and satisfy T, Theorem 1. Assume that (H 1 ) holds, and further assume that (H 2 ) and (H 3 ) are also true.
Then, the system (1) has at least one solution u * (t) ∈ X.
Proof.Based on Lemma 2, for all u ∈ X, we define two operators P, Q : X→X as follows: It is easy to see from ( 8) and ( 9) that Pu + Qv ∈ X, ∀ u, v ∈ X.Thus, the condition (i) in Lemma 4 holds.In addition, for all t ∈ D, u, v ∈ X, we have From (H 3 ) and ( 11), we know that P : X→X is contractive.Next, we shall show that Q : X→X is a completely continuous operator by using the Arzelá-Ascoli theorem.Indeed, for all t ∈ D, u ∈ X, we derive from (H 2 ) that where From ( 13), we conclude that, ∀ > 0, Thus, we verify that the condition (ii) is true.Therefore, according to Lemmas 2 and 4, one knows that there exists at least a fixed point u * (t) ∈ X such that u * (t) = (Pu * )(t) + (Qu * )(t), which is a solution of system (1).The proof is completed.
Proof.According to Lemma 2, we define an operator T : X→X as follows: Then, for all u, v ∈ X, we derive from (H 4 ) and (H 5 ) that (15) leads to By ( 16) and (H 5 ), we know that T : X→X is contractive.Thus, it follows from Lemmas 3 and 2 that the operator has a unique fixed point u * (t), which is a unique solution of system (1).The proof is completed.
Definition 5. Assume that for each ε > 0 and each solution z ∈ X of inequality (18), there is a constant C 2 > 0 and a unique solution u ∈ X of system (1) such that then system (1) is called Ulam-Hyers-Rassias stable (abbreviated as UHR-stable).
Definition 6. Assume that for each ε > 0 and each solution z ∈ X of inequality (18), there is a constant C 3 > 0 and a unique solution u ∈ X of system (1) such that then system (1) is called generalized Ulam-Hyers-Rassias stable (abbreviated as GUHR-stable).
Remark 3. A function z ∈ X is a solution of inequality (17) if and only if there exists a function φ ∈ X such that (1) ( ( Theorem 3. If all the conditions of Theorem 2 hold, then the system (1) is UH-stable and also GUH-stable.
Proof.Based on Lemma 2 and Remark 3, the solution z(t) of inequality ( 17) is expressed as By the Theorem 2 and Lemma 2, the unique solution u * (t) of ( 1) satisfies Similar to (15), it follows from ( 19) and ( 20) that Thus, (21) shows that system (1) is UH-stable and also GUH-stable.The proof is completed.Theorem 4. If all the conditions of Theorem 2 hold, then the system (1) is UHR-stable and also GUHR-stable.
Proof.By applying Lemma 2 and Remark 4, the solution z(t) of inequality ( 18) is formulated by Noting that ϕ ≥ 0 is non-decreasing and together with ( 19) and ( 22), we have Thus, (23) shows that system (1) is UHR-stable and also GUHR-stable.The proof is completed.

Applications
In this section, we will apply our results to deal with the existence and stability of solutions for two specific systems.
Example 2. Consider the following nonlinear fractional order Langevin equation Thus, we verify that the conditions (H 1 ), (H 4 ) and (H 5 ) are true.From Theorem 2, we know that Equation (24) has a unique solution u * (t) ∈ C 1 ([0, 1], R).Meanwhile, according to Theorems 3 and 4, we conclude that Equation (24) is stable in the sense of UH, GUH, UHR and GUHR, respectively.

Numerical Simulation
Let v(t) = ( ML D α 0 + − λ)u(t), then Equation ( 1) is transformed into a system of equations as follows: When α = β = 1, the fractional Langevin Equation ( 1) is a classical integer-order differential Langevin equation formulated by The system of equations equivalent to Equation ( 27) is formed as The numerical simulation algorithm in this manuscript is briefly stated as follows: Step 1, by applying Lemma 1, it is similar to the proof of Lemma 2 whereby the fractional Langevin equations ( 26) are transformed into a system of integral equations; Step 2, by calculating the derivative of integer order and simplifying, this system of integral equations becomes a system of delay differential equations of integer order; Step 3, by using ddesd toolbox in MATLAB, this system of delay differential equations of integer order can be numerically simulated.In addition, the corresponding integer-order Langevin Equations (28) are simulated by ode23 toolbox in MATLAB.Next, based on the above algorithm, we numerically simulate and discuss the solutions of ( 26) and (28) corresponding to Example 1 and Example 2, respectively.
Discussions.(a) Under the condition of the same values of system parameters, the simulations of solutions of Example 1 and its corresponding integer-order differential equation are shown in Figure 1 and Figure 2, respectively.Through the comparison of Figures 1 and 2, u(t) is very different in two aspects.On the one hand, although the solutions u(t) of fractional-order and integer-order equations increase monotonically in t ∈ [0, 1], the curvature of u(t) in Figure 1 is much larger than that in Figure 2. On the other hand, when t ∈ [0, 1], we find that 1 ≤ u(t) < 14 in Figure 1 and 1 ≤ u(t) < 3.5 in Figure 2, which is also very different.(b) Under the condition of the same values of system parameters, the simulations of solutions of Example 2 and the corresponding integer-order differential equation are shown in Figure 3 and Figure 4, respectively.Through the comparison of Figures 3 and 4, u(t) is greatly different as follows.In Figure 4, the solution u(t) of integer-order equation increases monotonically in t ∈ [0, 1].However, the solution u(t) of fractional equation is not monotonous and has maxima and minima in t ∈ [0, 1], shown in Figure 3. (c) Under the condition of the same values of system parameters, the simulations of Ulam-Hyers stability of Example 2 are shown in Figure 5.It follows from the images of ε = 0.1 and ε = 0.05 that, when ε→0 + , the solution curve of the inequality (17) almost coincides with that of Equation (25), which shows that Equation ( 25) is UH-stable.

Conclusions
It is well known that the Langevin equation is a powerful tool in describing the random motion of particles in a fluid.In a particularly complex viscous liquid, the integer-order Langevin equation that describes the motion of particles is no longer accurate.Some scholars have started using the fractional Langevin equation as a model to study this problem, and have achieved good results.However, the research results of these fractional Langevin systems are all centered on Riemann-Liouville or Caputo fractional derivatives.Unfortunately, Riemann-Liouville and Caputo fractional derivatives produce singularities under certain conditions, which renders their application difficult in certain physical fields.Interestingly, the fractional derivative with ML-kernel can eliminate the singularity.In this manuscript, we investigate the existence, uniqueness and UH-stability of solutions for the nonlinear fractional-order Langevin equation with ML-kernel.The theoretical analysis and numerical simulations of two examples verify the correctness and effectiveness of our main conclusions.Furthermore, the mathematical theories and methods employed in this paper can be used as a reference for the study of other fractional differential systems.In addition, considering the fact that the Langevin equation is a classical stochastic differential equation, we can further study the fractional random Langevin equation of nonsingular ML-kernel in future work to reveal the influence of random noise on the motion of free particles.

Funding:
The APC was funded by research start-up funds for high-level talents of Taizhou University.

Data Availability Statement:
No data were used to support this study.

Figure 1 .
Figure 1.Numerical simulation of solutions of Example 1.