Mathematical analysis of HIV/AIDS infection model with Caputo-Fabrizio fractional derivative

This manuscript is concerned to the existence and stability of HIV/AIDS infection model with fractional order derivative. The corresponding derivative is taken in Caputo-Fabrizio sense, which is a new approach for such type of biological models. With the help of Sumudu transform, some new results are handled. Further for the corresponding results, existence theory and uniqueness for the equilibrium solution are provided via using nonlinear functional analysis and fixed point theory due to Banach. Subjects: Advanced Mathematics; Applied Mathematics; Physics


Introduction
Acquired immune deficiency syndrome (AIDS) has developed into a global pandemic since the first patients were identified in 1981. Human immune deficiency virus HIV is a lenti virus that causes

PUBLIC INTEREST STATEMENT
In most of the situations in mathematical modeling of physical, biological process and phenomenon, the hereditary materials and memory terms cannot be explained accurately by using classical derivatives. As mathematical models are used in physics, biology, Earth science, meteorology and engineering disciplines (such as computer science, artificial intelligence, as well as in the social sciences. Therefore by using arbitrary order derivative one can easily explain the aforesaid terms with excellent way. Further by adopting Caputo Febrizo fractional derivative, the aforementioned terms can be further excellently described. In recent time, the investigation of dynamical systems through fractional derivative of Caputo Febrizo type is a warm area of research. AIDS. This serious disease destroys the immune system of human being which produce life-threatening opportunistic infections in the body. It is reported that 38.6 million people currently live with HIV/AIDS infection, 4.1 million people have been newly infected and 2.8 million AIDS deaths occurred in 2005. HIV/AIDS is an epidemic disease which spreading continuously all over the world and there have been few generator which continue it. In fact numbers of virus in the blood is a great indicator for the stages of the disease. In some situation, these stages are meant to correspond to CD4+ T-cell count ranges. In a normal healthy individual's peripheral blood, the level of CD4+ T-cells is between 800 and 1,200∕mm 3 and once this number fall to 200 or below in an HIV infected patient, the person is classified as having AIDS. The aforesaid viruses transfer from a calamitous illness into a chronic conditions. The said infection led to dramatic change in mobility and mortality from illness. Furthermore, despite these improvement on the biomedical front, the spreading of this epidemic continue and treatment remains unavailable to the overwhelming majority of those who require it. It causes expenditure of a very large amount of money in health care and research and destroy millions of peoples. In the medical field there have been many achievement for health care, but there is still no vaccine available for HIV/AIDS. For the determination of the transmission dynamics of HIV/AIDS disease, a mathematical model is a useful tool which also provide techniques to control the spread of these disease Ali, Zaman, and Alshomrani (2017). In Wang and Li (2006), a simple mathematical model introduced by Perelson for the primary infection with HIV/AIDS. In the field of mathematical modeling for HIV/AIDS infection this model has great importance and also many other models with HIV/AIDS infection have been proposed, which consider this model as their inspiration. Arqub and El-Ajou (2013), made modification in this model to study the evaluation of drug resistance. Sadegh Zebai et al. also present a mathematical model of HIV/AIDS infected and evaluate the stability of the model. According to clinical symptoms or viral load and CD4 C T cell count, 2-6 stages of infection before AIDS can be classified.
Mathematical models are powerful tools in this approach which help us to optimize the use of finite sources or simply to goal (the incidence of infection) control measures more impressively. The classical order model of HIV/AIDS is provided in Cai, Li, Ghosh, and Guo (2009) and given by where S represent number of individuals who are uninfected/suseptible, I represents number of individuals who are in the asymptomatic phase, J represent number of individuals who are in the symptomatic phase and A represent a full-blown AIDS group. The recruitment rate of the population is denoted by k, is natural death rate, c represents the average number of contacts of an individual per unit time, and b are probability of disease transmission per contact by an infective in the first stage and in the second stage, respectively. k 1 represent transfer rate of individuals from the asymptomatic phase to the symptomatic phase, k 2 represent transfer rate of individuals from the symptomatic phase to the AIDS cases, is treatment rate from the symptomatic phase to the asymptomatic phase and d denoted the disease related death rate. The parameter unite can be taken in numbers of cell per cubic millimeters. Recently, to study mathematical models by using fractional order derivative is an important area of research. Because, it has been found that the area involving fractional order differential equations have significant applications in various disciplines of science and technology, we refer few of them in Fabrizio (2015, 2016) and El-Saka (2014). In recent years, the fractional-order models were given much attentions, because the biological models that involved fractional-order derivative are more realistic and accurate as compared to the classical order models, for detail see Toledo-Hernandez, Rico-Ramirez, Iglesias-Silva, and (1) Diwekar (2014) and Wang, Yang, Ma, and Sun (2014). By adapting fractional order derivative in mathematical modeling is a global operator as compared to classical derivative which is local. Motivated by the above work, in this manuscript, we considered the model discussed in Cai et al. (2009), by taking the Caputo-Fabrizio fractional derivative of the system.
The model that we study in this paper is a Caputo-Fabrizio derivative and fractional order such that ∈ (0, 1] as given below Taking = 1, one can get the classical model as given in (8). With the help of Sumudu transform some new results are handled which demonstrate existence and uniqueness of equilibrium solutions to the proposed model. By existence of equilibrium solutions, we concluding the wellposedness of the proposed model. Since the Caputo fractional derivative contains singular kernel which often cannot explain many phenomenons excellently. Therefore, to describe various process and phenomenons of biology, and physical science, the Caputo-Fabrizio derivative and fractional can excellent describes the aforesaid. For further characteristics and features of the aforesaid derivative, one can read the articles refer as Atangana1 and Talkahtani (2015) and Baleanu, Mousalou, and Rezapour (2017).

Preliminaries
We first give the definitions of Caputo-Fabrizio derivative of fractional order. Caputo-Fabrizio derivative with fractional order has been considered with no singular kernel Caputo and Fabrizio (2015) and Losada and Nieto (2015). The connected anti-derivative turns out to be imperative at this instant following to the preface of the novel derivative, which was proposed by Losada and Nieto (2015).
(2) Definition 2.3 Let ∈ (0, 1). Then the fractional integral of order of a function f is defined by Remark 2.4 According to above definition, it has to be noted that the fractional integral of Caputo type of a function of order 0 < ≤ 1 is an average between function and its integral of order one. This therefore imposes The above equation generate an explicit formula for Keep in view the above derivative, reformulated the new Caputo derivative of order 0 < < 1 as

Derivation of the special solution
The aim of this section is to provide a special solution of the system (2) by applying the Sumudu transform on both sides of all equations of (2) together with an iterative method. To get this, we proceed as: Applying Sumudu transform on both sides of proposed model (2), we obtain

Rearranging, we obtain
Now applying the inverse Sumudu transform on both sides of Equation (3), we obtain We next obtain the following recursive formula And the solution of (4) is provided by S(t) = lim n→∞ S (n) (t),

Application of fixed-point theorem for stability analysis of iteration method
Let us suppose (X 1 , || ⋅ ||) as a Banach space and P as a self-map of X 1 . Let y n+1 = g(P, y n ) be particular recursive procedure. Suppose that, F(P) the fixed-point set of P has at least one element and that y n converges to a point p ∈ F(P). Let {x n ⊆X 1 } and define e n = ||x n+1 − g(P, x n )||. If lim n→∞ e n = 0 implies that lim n→∞ x n = p, then the iteration method y n+1 = g(P, y n ) is said to be Pstable. Analogously, we, therefore, consider that, our sequence {x n } has an upper bound, otherwise there is no possibility of convergence. The iteration will be P-stable, if all these conditions are satisfied for y n+1 = Py n which is also known as Picard's iteration.
Theorem 4.1 Let (X 1 , || ⋅ ||) be a Banach space and P be a self-map of X 1 satisfying for all x, y ∈ X 1 where 0 ≤ C, 0 ≤ c < 1. Suppose that P is Picard P−stable . Let us take into account the following recursive formula from (4) connected to (2).
is the fractional Lagrange multiplier.

Theorem 4.2 Let us defined a self-map P as
Proof In the first step of the proof we will show that P has a fixed point. For this, we evaluate the followings for all Let consider the first equation of (6) and taking norm of both hand sides, then without loss of generality, we have Thanks to the triangular inequality, the right-hand side of the Equation (7) becomes Upon further simplification, (8) yields that Since both the solutions play the same role, we shall assume in this case that (6) P(S (n) (t)) − P(S (m) (t)) = S (n) (t) − S (m) (t) ||P(S (n) (t)) − P(S (m) (t))|| ≤ ||S (n) (t) − S (m) (t)|| Replacing this in Equation (9), we obtain the following relation Since I n , S m ,I m and J n are bounded as they are convergent sequence, therefore, we can find four different positive constants, M, L, M 1 and L 1 for all t such that Now considering Equation (10) with (11), we obtain the following where f, g and h are functions from  −1 ( 1− + u M( ) ) .
In the same way, we will get where Thus the nonlinear P-self mapping has a fixed point. We next show that, P satisfies the conditions in Theorem 4.1. Let (12)-(15) hold and therefore using Then the above shows that condition of Theorem 4.1 exist for the nonlinear mapping P. Thus all the conditions in Theorem 4.2 are satisfied for the defined non-linear mapping P. Hence, P is Picard P-stable. ✷

Which yields that
In this way, we have completed the proof.

Conclusion
With the help of Picard successive approximation technique and Banach contraction theorem, we have investigated HIV/AIDS model with fractional order derivative. The arbitrary derivative of fractional order has been taken in the Caputo-Febrizo sense which does not contain singular kernel. The concerned results have been handled by coupling Sumudu transform with mentioned iterative techniques. Further, the existence and uniqueness results for equilibrium solutions have been proved by applying Banach theorem.