![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
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.
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.
1. 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 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 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 (Citation2017). In Wang and Li (Citation2006), 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 (Citation2013), 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 (Citation2009) and given by(1)
(1)
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 ,
is natural death rate, c represents the average number of contacts of an individual per unit time,
and
are probability of disease transmission per contact by an infective in the first stage and in the second stage, respectively.
represent transfer rate of individuals from the asymptomatic phase to the symptomatic phase,
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 Caputo and Fabrizio (Citation2015, Citation2016) and El-Saka (Citation2014). 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 Diwekar (Citation2014) and Wang, Yang, Ma, and Sun (Citation2014). 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. (Citation2009), 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
as given below
(2)
(2)
Taking , 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 (Citation2015) and Baleanu, Mousalou, and Rezapour (Citation2017).
2. 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 (Citation2015) and Losada and Nieto (Citation2015).
Definition 2.1
Let ,
,
then the new fractional order in Caputo derivative sense is recalled as
where the normalization function is denoted by with
. But, if the function does not belong to
, then the derivative can be reformulated as
Remark 2.2
If, we take ,
, then we will get the new Caputo derivative having fractional order as
In addition,
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 (Citation2015).
Definition 2.3
Let . 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 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 as
Definition 2.5
For any function over a set, the Sumudu transform will be given as
is defined by
Definition 2.6
The Sumudu transform of ordinary Caputo fractional-order derivative of a function is given by
In view of Definition 2.6, we recall the Sumudu transform of a function with Caputo-Fabrizio fractional derivative as below:
Definition 2.7
Atangana and Alkahtani (Citation2015) Let be a function for which the Caputo-Fabrizio exists, then the Sumudu transform of the Caputo-Fabrizio fractional derivative of f(t) is given as
3. 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(3)
(3)
Now applying the inverse Sumudu transform on both sides of Equation (3), we obtain
We next obtain the following recursive formula(4)
(4)
And the solution of (4) is provided by
4. Application of fixed-point theorem for stability analysis of iteration method
Let us suppose as a Banach space and P as a self-map of
. Let
be particular recursive procedure. Suppose that, F(P) the fixed-point set of P has at least one element and that
converges to a point
. Let
and define
. If
implies that
, then the iteration method
is said to be P-stable. Analogously, we, therefore, consider that, our sequence
has an upper bound, otherwise there is no possibility of convergence. The iteration will be P-stable, if all these conditions are satisfied for
which is also known as Picard’s iteration.
Theorem 4.1
Let be a Banach space and P be a self-map of
satisfying
for all where
. Suppose that P is Picard
stable . Let us take into account the following recursive formula from (4) connected to (2).
where is the fractional Lagrange multiplier.
Theorem 4.2
Let us defined a self-map P as
is stable in
if
(5)
(5)
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 (6)
(6)
Let consider the first equation of (6) and taking norm of both hand sides, then without loss of generality, we have(7)
(7)
Thanks to the triangular inequality, the right-hand side of the Equation (7) becomes(8)
(8)
Upon further simplification, (8) yields that(9)
(9)
Since both the solutions play the same role, we shall assume in this case that
Replacing this in Equation (9), we obtain the following relation(10)
(10)
Since ,
,
and
are bounded as they are convergent sequence, therefore, we can find four different positive constants, M, L,
and
for all t such that
(11)
(11)
Now considering Equation (10) with (11), we obtain the following(12)
(12)
where f, g and h are functions from
In the same way, we will get(13)
(13)
(14)
(14)
(15)
(15)
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.
5. Uniqueness of the special solution
In this section, we show that the special solution of Equation (2) using the iteration method is unique. We shall first assume that, Equation (2) has an exact solution via which, the special solution converges for a large number m. We consider the following Hilbert space which can be defined as
We now, consider the following operator
The aim of this part is to prove that the inner product of
where ,
,
and
, are special solution of system. However,
(16)
(16)
We shall evaluate the first equation in the system without loss of generality(17)
(17)
Since both solutions play almost the same role, we can assume that,
Then the Equation (17) becomes
Based on the relationship between norm and the inner product, we obtain the following(18)
(18)
Repeating the same fashion, from the second, third and fourth equations of the system (16), we can obtain as follow(19)
(19)
Putting Equations (18) and (19) in Equation (16), we get(20)
(20)
But, for sufficiently large values of with
both the solutions converge to the exact solution, using the topological concept, there exist three very small positive parameters
and
such that
Thus, plugging the exact solution in the right-hand side of Equation (20) and applying the triangular inequality by taking We obtain
As l is a small positive parameter, therefore, on the basis of topological idea, we have
But, it is obvious that
Therefore, we have
Which yields that
In this way, we have completed the proof.
6. 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.
7. Authors contribution
All authors equally contributed this paper and approved the final version.
Funding
The authors received no direct funding for this research.
Acknowledgements
We all authors are really thankful to the reviewer for their useful comments and suggestions.
Additional information
Funding
Notes on contributors
Samia Bushnaq
Samia Bushnaq is an assistant professor in the Department of Basic Sciences, Princess Sumaya University for Technology, Amman 11941, Jordan. Her research interest is devoted to Numerical solutions of FDEs and mathematical analysis and simulation of Biological Models.
Sajjad Ali Khan
Sajjad Ali Khan is a PhD scholar in the department of Mathematics University of Malakand Dir(L), KPK, Pakistan under the supervision of Gul Zaman and Kamal Shah. His area of research is the study of mathematical models of infectious disease using Caputo Febrizo fractional derivative.
Kamal Shah
Kamal Shah is a asistant professor in department of Mathematics University of Malakand Dir(L), KPK, Pakistan since 2007 till date. His research interest is to study FDEs and their numerical solutions, qualitative theory, Fixed point theory and applications, mathematical modeling and simulations.
Gul Zaman
Gul Zaman is an associate professor in the department of Mathematics University of Malakand Dir(L), KPK, Pakistan. His research is devoted to Fluid Mechanics, mathematical modeling of infectious disease and their simulation, existence theory.
References
- Ali, N., Zaman, G., & Alshomrani, A. S. (2017). Optimal control strategy of HIV-1 epidemic model for recombinant virus. Cogent Mathematics, 12(4), 1–10.
- Arqub, O. A., & El-Ajou, A. (2013). Solution of the fractional epidemic model by homotopy analysis method. Journal of King Saud University, 25, 73–81.
- Atangana1, A., & Talkahtani, B. S. (2015). Extension of the resistance, inductance, capacitance electrical circuit to fractional derivative without singular kernel. Advances in Mechanical Engineering, 7(6), 1–6.
- Atangana, A., & Alkahtani, B. S. T. (2015). Analysis of the Keller-Segel model with a fractional derivative without singular kernel. Entropy, 17, 4439–-4453.
- Baleanu, D., Mousalou, A., & Rezapour, S. (2017). A new method for investigating approximate solutions of some fractional integro-differential equations involving the Caputo-Fabrizio derivative. Advances in Difference Equations, 2017, 12.
- Cai, L., Li, X., Ghosh, M., & Guo, B. (2009). Stability analysis of an HIV/AIDS epidemic model with treatment. Journal of Computational and Applied Mathematics, 229, 313–323.
- Caputo, M., & Fabrizio, M. (2015). A new definition of fractional derivative with out singular kernel. Progress in Fractional Differentiation and Applications, 1, 73–85.
- Caputo, M., & Fabrizio, M. (2016). Applications of new time and spatial fractional derivatives with exponential kernels. Progress in Fractional Differentiation and Applications., 2, 1–11.
- El-Saka, H. A. A. (2014). The fractional-order SIS epidemic model with variable population size. Journal of the Egyptian Mathematical Society, 22, 50–54.
- Losada, J., & Nieto, J. J. (2015). Properties of a new fractional derivative without singular kernel. Progress in Fractional Differentiation and Applications, 1, 87–92.
- Toledo-Hernandez, R., Rico-Ramirez, V., Iglesias-Silva, G. A., & Diwekar, U. M. (2014). A fractional calculus approach to the dynamic optimization of biological reactive systems. Part I: Fractional models for biological reactions. Chemecal Engineering Science, 117, 217–228.
- Wang, L., & Li, M. Y. (2006). Mathematical analysis of the global dynamics of a model for HIV infection of CD4+ T-cells. Mathematical Biosciences, 200, 44–57.
- Wang, Z., Yang, D., Ma, T., & Sun, N. (2014). Stability analysis for nonlinear fractional-order systems based on comparison principle. Nonlinear Dynamics, 75(1--2), 387–402.