Global Stability of Fractional Order HIV/AIDS Epidemic Model under Caputo Operator and Its Computational Modeling

Abstract

The human immunodeficiency virus (HIV) causes acquired immunodeficiency syndrome (AIDS), which is a chronic and sometimes fatal illness. HIV reduces an individual’s capability against infection and illness by demolishing his or her immunity. This paper presents a new model that governs the dynamical behavior of HIV/AIDS by integrating new compartments, i.e., the treatment class T. The steady-state solutions of the model are investigated, and accordingly, the threshold quantity R0 is calculated, which describes the global dynamics of the proposed model. It is proved that for R0 less than one, the infection-free state of the model is globally asymptotically stable. However, as the threshold number increases by one, the endemic equilibrium becomes globally asymptotically stable, and in such case, the disease-free state is unstable. At the end of the paper, the analytic conclusions obtained from the analysis of the ordinary differential equation (ODE) model are supported through numerical simulations. The paper also addresses a comprehensive analysis of a fractional-order HIV model utilizing the Caputo fractional differential operator. The model’s qualitative analysis is investigated, and computational modeling is used to examine the system’s long-term behavior. The existence/uniqueness of the solution to the model is determined by applying some results from the fixed points of the theory. The stability results for the system are established by incorporating the Ulam–Hyers method. For numerical treatment and simulations, we apply Newton’s polynomial and the Toufik–Atangana numerical method. Results demonstrate the effectiveness of the fractional-order approach in capturing the dynamics of the HIV/AIDS epidemic and provide valuable insights for designing effective control strategies.

Keywords

• fractional modeling
• HIV/AIDS epidemic
• Newton’s polynomial
• reproduction number
• Toufik–Atangana
• Ulam–Hyers stability