In Vivo HIV Dynamics, Modeling the Interaction of HIV and Immune System via Non-Integer Derivatives

Abstract

The economic burden of HIV extends beyond the individual level and affects communities and countries. HIV can lead to decreased economic growth due to lost productivity and increased healthcare costs. In some countries, the HIV epidemic has led to a reduction in life expectancy, which can impact the overall quality of life and economic prosperity. Therefore, it is significant to investigate the intricate dynamics of this viral infection to know how the virus interacts with the immune system. In the current research, we will formulate the dynamics of HIV infection in the host body to conceptualize the interaction of T-cells and the immune system. The recommended model of HIV infection is presented with the help of fractional calculus for more precious outcomes. We introduce numerical methods to demonstrate how the input parameters affect the output of the system. The dynamical behavior and chaotic nature of the system are visualized with the variation of different input factors. The system’s tracking path has been numerically depicted and the impact of the viruses on T-cells has been demonstrated. In addition to this, the key factors of the system has been predicted through numerical findings. Our results predict that the strong non-linearity of the system is responsible for the chaos and oscillation, which are so closely related. The chaotic parameters of the system are highlighted and are recommended for the control of the chaos of the system.

1. Introduction

With HIV infection, the virus attacks and weaken the immune system, which is the body’s defense against infections and diseases. The virus targets CD4${}^{+}$ T-cells, which are a type of white blood cell that plays a central role in the immune response. As HIV continues to replicate and destroy CD4${}^{+}$ T-cells, the immune system becomes weaker and less able to fight off infections and diseases. This can lead to a wide range of opportunistic infections and cancers that would not normally affect healthy individuals. Without treatment, HIV can progress to AIDS, which is a condition characterized by severe immune system damage and a high risk of developing life-threatening infections and cancers. Effective treatment for HIV involves antiretroviral therapy (ART), which uses a combination of medications to suppress the virus and prevent disease progression. With effective treatment, people living with HIV can achieve an undetectable viral load, which means that the virus is no longer detectable in their blood and they are much less likely to transmit HIV to others. However, even with effective treatment, HIV can still have long-term effects on the immune system, and people living with HIV may be at increased risk for certain infections and cancers.

The infection of HIV infects the host immune system and damage vital organs like the heart, kidneys, and brain, which finally results in death. As there is no effective treatment for this infection, but some medications can significantly improve the patient’s health, while excessive use of these medications has side effects on the host. ART is being used by more people than ever before. People have been seen to occasionally experience symptoms that resemble sore throats, rashes, headaches, fevers, and influenza. Sometimes, the symptoms are dangerous, such as fever, weight gain, diarrhea, coughing and swollen lymph nodes. If left untreated, these can progress to dangerous conditions. Since the first incidence of HIV was discovered in the early 1980s, millions of people globally have become infected, with HIV becoming the most deadly infection in the world [1,2,3]. HIV continues to be a major cause of death globally, with sub-Saharan Africa suffering the worst effect prevalence range: 12 to 42% [2]. The poorest communities bear the heaviest burden of HIV/AIDS and have the least access to medications and other preventive measures. Therefore, a more effective strategy is needed for the prevention, control, and management of this viral infection, which includes raising awareness and providing the optimum ART drug combinations.

Over the years, in vivo studies for the dynamics of HIV infection have been conducted with the goal of better visualization of how the HIV virus and the cells of the body interact. Such knowledge has shown to be extremely helpful, particularly in the creation of ARTs and in the treatment of HIV. In contrast to the past, when research relied on clinical trials, currently, mathematical models of varying degrees of complexity are utilized to simulate and analyze such interactions [4]. As a result, numerous experts in the field of epidemiological modeling [5,6,7,8] have started to create novel models that could be applied to examine the infection. Numerous models that have been created describe and shine a light on different elements that arise from the interaction between HIV and T-cells. Mathematical models for HIV and the immune system are used to understand the dynamics of the virus and the immune response, and to explore strategies for managing the disease. The models typically involve differential equations that describe the interactions between the virus and the immune cells, as well as the effects of treatments and other interventions [9]. Most of these fundamental models comprise at least three state variables, including HIV virion and uninfected and infected CD4${}^{+}$ T-cells for the purpose of providing the most basic knowledge on regulating the viral propagation. For example, Ref. [10] created and examined the interaction of healthy $CD{4}^{+}$ T-cells and infected $CD{4}^{+}$ T-cells and the concentration of HIV viruses.

It is well known that the models have undergone significant improvements, and several more advanced models have been developed in the literature for biological processes [11,12,13]. A mathematical model of the entire dynamics of HIV infection was developed in [14]. In [4], the model was made to examine how HIV dynamics changed early on in infection, where the findings demonstrated that the initial stage of HIV infection was marked by extremely high viral persistence. A mathematical modeling approach was used by the authors in [15] to demonstrate how migration affects the spread of HIV and AIDS cases. In [16], the researchers investigated an optimal treatment strategy, known as target-oriented treatment (TOT), for HIV infections using a mathematical model. Specifically, they focused on the scenario where a specified number of CD4${}^{+}$ T-cells must be achieved during the treatment period. Although antiretroviral therapy can effectively suppress HIV, it often does not fully restore the immune system to its pre-infection state [17]. Thus, further research is required to focus on ways to improve immune reconstitution in individuals with HIV. Other studies have attempted to understand the function of killer T-cells in blocking viral replication in the human host. In [18], the researchers considered $CD{8}^{+}$ T-cells with HIV viruses and showed how the immune system is influenced by the viruses. Some other dynamics considered three-dimensional dynamics by taking viruses and $CD{4}^{+}$ T-cells with the help of mathematical models. The model’s failure to take the HIV virion into account was a significant flaw. The authors in [7] employed a five-dimensional model of ODE’s to illustrate how $CD{4}^{+}$ T-cells, virion, defense cells, and ARTs are related to one another. The findings highlighted how crucial $CD{8}^{+}$ T-cells are in battling HIV virion during acute infection. It is known that the inability of ARTs to reach all cells and exponential viral multiplication caused the disease to spread more widely. In this study, we opt to conceptualize the intricate phenomena of $CD{4}^{+}$ T-cells, $CD{8}^{+}$ T-cells, and virus cells during HIV infection to show how it affects the immune system. Furthermore, we aim to highlight the main characteristics of the phenomenon and provide effective collaboration with the corresponding departments for the prevention of the infection.

This study is organized as follows. In Section 2, we introduce the HIV dynamics to illustrate how the HIV virus interacts with the immune cells. In Section 3, the basic results of fractional calculus are listed for the analysis of the recommended model of HIV infection. For more precise results, the recommended system of HIV infection is presented through fractional derivative in Section 4. In Section 5, the dynamical behavior of the system is investigated through numerical skills. In Section 6, we provide a conclusion.

2. Formulation of HIV Dynamics

In this section, we will introduce a mathematical model for the infection of HIV to represent the interaction of viruses and immune cells. It is well known that numerous mathematical models have been established in the literature for this viral infection [19,20,21]. The authors presented the dynamics of HIV with different assumptions for the input of the system. According to [22], the HIV transmission phenomenon is as follows:

$$\left\{\begin{array}{ccc}\hfill \frac{dT}{dt}& =& s-{\mu }_T\mathit{T}+r\mathit{T}(1-\frac{\mathit{T}+\mathit{I}}{T_{max}})-k\mathit{VT},\hfill \\ \hfill \frac{dI}{dt}& =& k\mathit{VT}-{\mu }_I\mathit{I},\hfill \\ \hfill \frac{dV}{dt}& =& \mathit{I}N{\mu }_I-\mathit{V}{\mu }_V,\hfill \end{array}\right.$$

(1)

where s represents the rate of recruitment of T-cells into the host body and r is the growth rate of the population of healthy T-cells, the terms ${\mu }_T$, ${\mu }_I$, and ${\mu }_V$ are used to indicate the death rate of T-cells, infected T-cells, and the mortality of HIV viruses, respectively, k indicates the infection rate of healthy T-cells, and N is the number of cells that infected T-cells reproduce. The HIV model from Perelson and Nelson [23] is given by the following system:

$$\left\{\begin{array}{ccc}\hfill \frac{d\mathit{T}}{dt}& =& r\mathit{T}(1-\frac{\mathit{T}}{T_{max}})-\frac{\mathit{VT}\eta }{1+\alpha \mathit{V}},\hfill \\ \hfill \frac{d\mathit{I}}{dt}& =& \frac{\mathit{VT}\eta }{1+\alpha \mathit{V}}-{\mu }_I\mathit{I},\hfill \\ \hfill \frac{d\mathit{V}}{dt}& =& N\mathit{I}{\mu }_I-{\mu }_V\mathit{V}.\hfill \end{array}\right.$$

(2)

In [22], a variable source term $sexp(-\kappa \mathit{V})$ was introduced instead of the constant term for T-cells. This new source term reflects the number of new healthy T-cells generated by the thymus, which is dependent on the concentration of the viral load. As the size of the viral load increases, the production of healthy T-cells is observed to decrease. Thus, the source term is now considered to be a variable instead of a constant. Then, according to [24], we have:

$$\left\{\begin{array}{ccc}\hfill \frac{d\mathit{T}}{dt}& =& sexp(-\kappa \mathit{V})-{\mu }_T\mathit{T}+r\mathit{T}(1-\frac{\mathit{T}}{T_{max}})-\frac{k\mathit{VT}}{1+{\alpha }_1\mathit{V}}-\frac{\alpha \mathit{IT}}{1+{\alpha }_2\mathit{I}},\hfill \\ \hfill \frac{d\mathit{I}}{dt}& =& \frac{k\mathit{VT}}{1+{\alpha }_1\mathit{V}}+\frac{\alpha \mathit{IT}}{1+{\alpha }_2\mathit{I}}-{\mu }_I\mathit{I},\hfill \\ \hfill \frac{d\mathit{V}}{dt}& =& N{\mu }_I\mathit{I}-{\mu }_V\mathit{V}-\frac{k\mathit{VT}}{1+{\alpha }_1\mathit{V}},\hfill \end{array}\right.$$

(3)

where the rate of cellular infection and $\alpha$ is the effectiveness of a protease inhibitor, respectively. In our model, we assume the concentration of $CD{8}^{+}$ and activated $CD{8}^{+}$ T-cells by $\mathit{C}$ and $\mathit{A}$, respectively. The thymus continuously produces $CD{8}^{+}$ T-cells at a rate of ${\aleph }_C$ and naturally eliminates them at a rate of ${\mu }_C$. In the presence of HIV virion, $CD{8}^{+}$ T-cell activation proceeds at a different rate. Furthermore, $CD{8}^{+}$ T-cells are transformed into activated $CD{8}^{+}$ T-cells in the presence of HIV virion, which are produced at a rate of $\eta$ and naturally depleted at a rate of ${\mu }_A$. Thus, the dynamics of HIV infection is given by:

$$\left\{\begin{array}{ccc}\hfill \frac{d\mathit{T}}{dt}& =& {\aleph }_T-{\mu }_T\mathit{T}+r\mathit{T}(1-\frac{\mathit{T}}{T_{max}})-k\mathit{VT}-\alpha \mathit{IT},\hfill \\ \hfill \frac{d\mathit{I}}{dt}& =& k\mathit{VT}+\alpha \mathit{IT}-{\mu }_I\mathit{I}-\gamma \mathit{I}\mathit{A},\hfill \\ \hfill \frac{d\mathit{V}}{dt}& =& N{\mu }_I\mathit{I}-\mathit{V}{\mu }_V-\mathit{VT}k,\hfill \\ \hfill \frac{d\mathit{C}}{dt}& =& {\aleph }_C-\eta \mathit{IC}-{\mu }_C\mathit{C},\hfill \\ \hfill \frac{d\mathit{A}}{dt}& =& \mathit{I}\eta \mathit{C}-{\mu }_A\mathit{A},\hfill \end{array}\right.$$

(4)

where ${\aleph }_T=sexp(-\kappa \mathit{V})$ and the initial condition is given as $\mathit{T}\left(0\right)={\mathit{T}}_0,\mathit{I}\left(0\right)={\mathit{I}}_0$, $\mathit{V}\left(0\right)={\mathit{V}}_0$, $\mathit{C}\left(0\right)={\mathit{C}}_0$, and $\mathit{A}\left(0\right)={\mathit{A}}_0$. The above model (4) will be investigated to visualize the interaction of immune system and HIV viruses.

3. Fundamentals of Fractional Calculus

In this section, we will presented the basic results of fractional theory for the analysis of our system. We will use the well-known Liouville–Caputo and Atangana–Baleanu operators in this study (see [25,26,27]). The basic results related to these operators are given below.

Definition 1. 

If h is taken such that h in a way that $h:{\mathbb{R}}^{+}\to \mathbb{R}$, the integral in fractional Liouville–Caputo structure is as follows:

$$\begin{array}{c}\hfill \left(I_t^{\beta }h\right)\left(t\right)=\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^t{(t-x)}^{\beta -1}h\left(x\right)dx,\end{array}$$

where Υ is used for the gamma function and β is the order of the operator.

Definition 2. 

If h is taken such that $h\in B^m$, the Liouville–Caputo fractional derivative is introduced as given below:

$$\begin{array}{c}\hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }h\right)\left(t\right)=\left(I_t^{m-\alpha }{D}^{m}h\right)\left(t\right)=\frac{1}{{\rm Y}(m-\beta )}{\int }_0^t\frac{h^m\left(x\right)}{{(t-x)}^{\beta +m-1}}dx,\end{array}$$

where $m-1<\beta <m\in M$, $D^m=\frac{d^m}{d{t}^m}$, and ${}^C{D}_t^{\beta }$ is closer to ${\mathrm{h}}^{\prime }\left(t\right)$ as beta gets closer to one.

Secondly, we illustrated the operator of Atangana–Baleanu. It captures the occurrence of biological processes with a non-local and non-singular kernel. The results obtained with this novel derivative have shown to be more reliable than those obtained with other alternatives. The fundamental concepts of this operator are given below.

Definition 3. 

Take h such that $h\in F^1(c,d)$, where β is in the range [0, 1] and $c<d$, then the Atangana–Baleanu operator is given by:

$$\begin{array}{c}\hfill \left({{}_c{}^{\mathrm{ABC}}D}_t^{\beta }h\right)\left(t\right)=\frac{C\left(\beta \right)}{1-\beta }{\int }_c^t{h}^{\prime }(\Phi )E_{\beta }\left[-\beta \frac{{(t-x)}^{\beta }}{1-\beta }\right]dx.\end{array}$$

Definition 4. 

The fractional integral of the Atangana–Baleanu operator for h is as follows:

$$\begin{array}{c}\hfill \left({{}_c{}^{\mathrm{ABC}}I}_t^{\beta }h\right)\left(t\right)=\frac{1-\beta }{C\left(\beta \right)}h\left(t\right)+\frac{\beta }{C\left(\beta \right){\rm Y}\left(\beta \right)}{\int }_c^{t}h\left(x\right){(t-x)}^{\beta -1}dx,\end{array}$$

as the fractional order approaches 0, the original function is attained.

Theorem 1 

(see [28]). If you assume that $\mathrm{h}$ such that $\mathrm{h}\in B[c,d]$, then the following is true:

$$\begin{array}{c}\hfill ∥\left({{}_c{}^{\mathrm{ABC}}D}_t^{\beta }h\right)\left(t\right)∥<\frac{C\left(\beta \right)}{1-\beta }\parallel h\left(t\right)\parallel ,\mathrm{where}\parallel h\left(t\right)\parallel =ma{x}_{c\le t\le d}\left|h\left(t\right)\right|.\end{array}$$

Additionally, the Lipschitz requirement for this ABC derivative is met:

$$\begin{array}{c}\hfill ∥\left({{}_c{}^{\mathrm{ABC}}D}_t^{\beta }{h}_1\right)\left(t\right)-\left({{}_c{}^{\mathrm{ABC}}D}_t^{\beta }{h}_2\right)\left(t\right)∥<{\varphi }_1\parallel h_1\left(t\right)-h_2\left(t\right)\parallel .\end{array}$$

4. Fractional Dynamics of HIV Infection

These days, the theory of fractional calculus is frequently used by researchers to formulate epidemic models for infectious diseases due to its accurate findings. These models have the ability to represent the crossover and non-local nature of real-world problems. More detailed results have been obtained recently for infections diseases through fractional-order modeling. Thus, the dynamics of HIV infection will be formulated in the structure of fractional calculus to obtain more precious results in this section. First, we illustrate the dynamics of HIV infection with the help of the Liouville–Caputo derivative and then with the operator of Atangana–Baleanu in the Liouville–Caputo structure.

4.1. HIV Dynamics via the Liouville–Caputo Operator

In this subsection, we will illustrate the intricate phenomena of HIV infection through the Liouville–Caputo operator. The main aim is to conceptualize the dynamics of HIV and the immune system. We will also present a numerical scheme to highlight the tracking path of the recommended model with the variability of different factors. Our model (4) according to the Liouville–Caputo framework is as follows:

$$\left\{\begin{array}{ccc}\hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }\mathit{T}\right)\left(t\right)& =& {\aleph }_T-{\mu }_T\mathit{T}+r\mathit{T}(1-\frac{\mathit{T}}{T_{max}})-k\mathit{VT}-\alpha \mathit{IT},\hfill \\ \hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }\mathit{I}\right)\left(t\right)& =& k\mathit{VT}+\alpha \mathit{IT}-{\mu }_I\mathit{I}-\gamma \mathit{I}\mathit{A},\hfill \\ \hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }\mathit{V}\right)\left(t\right)& =& N{\mu }_I\mathit{I}-\mathit{V}{\mu }_V-\mathit{VT}k,\hfill \\ \hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }\mathit{C}\right)\left(t\right)& =& {\aleph }_C-\eta \mathit{IC}-{\mu }_C\mathit{C},\hfill \\ \hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }\mathit{A}\right)\left(t\right)& =& \eta \mathit{I}\mathit{C}-\mathit{A}{\mu }_A,\hfill \end{array}\right.$$

(5)

where ${{}_0{}^{\mathrm{LC}}D}_t^{\beta }$ stands for the Liouville–Caputo fractional operator. To show the dynamical behavior of the recommended model (5), we represent a numerical scheme for the fractional Liouville–Caputo non-linear ordinary differential system [28]. For the scheme, we take the following system:

$$\begin{array}{c}\hfill \left({{}_0{}^{\mathrm{LC}}D}_t^{\beta }v\right)\left(t\right)=h(t,v\left(t\right)),\end{array}$$

(6)

following the application of the fundamental theorem, we have:

$$\begin{array}{c}\hfill v\left(t\right)-v\left(0\right)=\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t}h(\xi ,v\left(\xi \right)){(t-\xi )}^{\beta -1}d\xi ,\end{array}$$

(7)

here, we take the time $t=t_{k+1}$, $k=0,1,\dots$, and obtain:

$$\begin{array}{c}\hfill v\left(t_{k+1}\right)-v\left(0\right)=\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t_{k+1}}{(t_{k+1}-t)}^{\beta -1}h(t,v\left(t\right))dt,\end{array}$$

(8)

and

$$\begin{array}{c}\hfill v\left(t_k\right)-v\left(0\right)=\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t_k}{(t_k-t)}^{\beta -1}h(t,v\left(t\right))dt.\end{array}$$

(9)

Equations (8) and (9) imply that:

$$\begin{array}{ccc}\hfill v\left(t_{k+1}\right)& =& v\left(t_k\right)+\underset{{\mathcal{T}}_{\beta ,1}}{\underbrace{\frac{1}{\Gamma \left(\beta \right)}{\int }_0^{t_{k+1}}{(t_{k+t}-t)}^{\beta -1}f(t,v\left(t\right))dr}}\hfill \\ & & -\underset{{\mathcal{T}}_{\beta ,2}}{\underbrace{\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t_k}{(t_k-t)}^{\beta -1}h(t,v\left(t\right))dt}},\hfill \end{array}$$

(10)

where

$$\begin{array}{c}\hfill {\mathcal{T}}_{\beta ,1}=\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t_{k+1}}{(t_{k+1}-t)}^{\beta -1}h(t,v\left(t\right))dt,\end{array}$$

(11)

and

$$\begin{array}{c}\hfill {\mathcal{T}}_{\beta ,2}=\frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t_k}{(t_k-t)}^{\beta -1}h(t,z\left(t\right))dt.\end{array}$$

(12)

Then, $\mathrm{h}(t,v(t\left)\right)$ is approximated in the following manner:

$$\begin{array}{ccc}\hfill P\left(t\right)& \simeq & \frac{t-t_{k-1}}{t_k-t_{k-1}}h(t_k,v_k)+\frac{t-t_k}{t_{k-1}-t_k}h(t_{k-1},v_{k-1})\hfill \\ & =& \frac{h(t_k,v_k)}{f}(t-t_{k-1})-\frac{h(t_{k-1},v_{k-1})}{f}(t-t_k).\hfill \end{array}$$

(13)

Following the simplification, we obtain:

$$\begin{array}{ccc}\hfill {\mathcal{T}}_{\beta ,1}& =& \frac{h(t_k,z_k)}{f{\rm Y}\left(\beta \right)}{\int }_0^{{}_{k+1}}{(t_{k+1}-t)}^{\beta -1}(t-t_{k-1})dt\hfill \\ & & -\frac{h(t_{k-1},v_{k-1})}{f{\rm Y}\left(\beta \right)}{\int }_0^{t_{k+1}}{(t_{k+1}-t)}^{\beta -1}(t-t_k)dt.\hfill \end{array}$$

(14)

Furthermore, we have:

$$\begin{array}{ccc}\hfill {\mathcal{T}}_{\beta ,1}& =& \frac{h(t_k,v_k)}{f{\rm Y}\left(\beta \right)}\left[\frac{2h}{\beta }{t}_{k+1}^{\beta }-\frac{t_{k+1}^{\beta +1}}{\beta +1}\right]\hfill \\ & & -\frac{h(t_{k-1},v_{k-1})}{f{\rm Y}\left(\beta \right)}\left[\frac{f}{\beta }{t}_{k+1}^{\beta }-\frac{1}{\beta +1}{t}_{k+1}^{\beta +1}\right].\hfill \end{array}$$

(15)

Next, we have:

$$\begin{array}{ccc}\hfill {\mathcal{T}}_{\beta ,2}& =& \frac{1}{{\rm Y}\left(\beta \right)}{\int }_0^{t_k}{(t_k-t)}^{\beta -1}[\frac{h(t_k,v_k)}{f}(t-t_{k-1})\hfill \\ & & -\frac{h(t_{k-1},v_{k-1})}{f}(t-t_k)]dt.\hfill \end{array}$$

(16)

Simplifying even more, we obtain:

$$\begin{array}{c}\hfill {\mathcal{T}}_{\beta ,2}=\frac{h(t_k,v_k)}{f{\rm Y}\left(\beta \right)}\left[\frac{f}{\beta }{t}_k^{\beta }-\frac{t_k^{\beta +1}}{\beta +1}\right]+\frac{h(t_{k-1},v_{k-1})}{f{\rm Y}\left(\beta \right)}\left[\frac{1}{\beta +1}{t}_k^{\beta +1}\right].\end{array}$$

(17)

Equations (16) and (17) were combined into Equation (10), and the resulting approximation is shown below:

$$\begin{array}{ccc}\hfill v\left(t_{k+1}\right)& =& v\left(t_k\right)+\frac{h(t_k,v_k)}{f{\rm Y}\left(\beta \right)}\left[\frac{2f{t}_{k+1}^{\beta }}{\beta }-\frac{t_{k+1}^{\beta +1}}{\beta +1}+\frac{h}{\beta }{t}_k^{\beta }-\frac{t_{k+1}^{\beta +1}}{\beta +1}\right]\hfill \\ & & +\frac{h(t_{k-1},v_{k-1})}{f{\rm Y}\left(\beta \right)}\left[-\frac{f}{\beta }{t}_{k+1}^{\beta }+\frac{t_{k+1}^{\beta +1}}{\beta +1}+\frac{t_k^{\beta +1}}{\beta +1}\right].\hfill \end{array}$$

(18)

Through this numerical method, we represent the system (4) for HIV infection. More specifically, the tracking path of the system will be represented in the Liouville–Caputo system through the above numerical scheme.

4.2. HIV Dynamics via the ABC Operator

Here, we will illustrate the recommended model (4) using the Atangana–Baleanu fractional derivative. The impact of various factors on the HIV system will then be visualized using a numerical method for the ABC derivative. Our model (4) is first represented by the Atangana–Baleanu derivative as:

$$\left\{\begin{array}{ccc}\hfill \left({{}_0{}^{\mathrm{ABC}}D}_t^{\beta }\mathit{T}\right)\left(t\right)& =& {\aleph }_T-{\mu }_T\mathit{T}+r\mathit{T}(1-\frac{\mathit{T}}{T_{max}})-k\mathit{VT}-\alpha \mathit{IT},\hfill \\ \hfill \left({{}_0{}^{\mathrm{ABC}}D}_t^{\beta }\mathit{I}\right)\left(t\right)& =& k\mathit{VT}+\alpha \mathit{IT}-{\mu }_I\mathit{I}-\gamma \mathit{I}\mathit{A},\hfill \\ \hfill \left({{}_0{}^{\mathrm{ABC}}D}_t^{\beta }\mathit{V}\right)\left(t\right)& =& N{\mu }_I\mathit{I}-\mathit{V}{\mu }_V-\mathit{VT}k,\hfill \\ \hfill \left({{}_0{}^{\mathrm{ABC}}D}_t^{\beta }\mathit{C}\right)\left(t\right)& =& {\aleph }_C-\eta \mathit{IC}-{\mu }_C\mathit{C},\hfill \\ \hfill \left({{}_0{}^{\mathrm{ABC}}D}_t^{\beta }\mathit{A}\right)\left(t\right)& =& \eta \mathit{I}\mathit{C}-\mathit{A}{\mu }_A,\hfill \end{array}\right.$$

(19)

where $\mathit{T}\left(0\right)={\mathit{T}}_0,\mathit{I}\left(0\right)={\mathit{I}}_0$, and $\mathit{V}\left(0\right)={\mathit{V}}_0$ denote the initial conditions and ${{}_0{}^{\mathrm{ABC}}D}_t^{\beta }$ indicates the ABC fractional derivative. In the next step, we will represent a numerical scheme to illustrate the above dynamics of HIV infection. To do so, we take the following general fractional system:

$$\begin{array}{c}\hfill \left({{}_0{}^{\mathrm{ABC}}D}_t^{\beta }v\right)\left(t\right)=\mathrm{C}(t,v\left(t\right)),\end{array}$$

through mathematical skills, we obtain the following:

$$\begin{array}{cc}v\left(t\right)-v\left(0\right)=\frac{(1-\beta )}{ABC\left(\beta \right)}\mathcal{C}(t,v\left(t\right))\hfill & \\ & \hfill +\frac{\beta }{ABC\left(\beta \right)\times {\rm Y}\left(\beta \right)}{\int }_0^t\mathcal{C}(\Phi ,v(\Phi )){(t-\Phi )}^{\Phi -1}d\Phi .\end{array}$$

(20)

At this moment, $t=t_{b+1}$, $l=0,1,2,\dots$, the result is as follows:

$$\begin{array}{ccc}\hfill v\left(t_{b+1}\right)-v\left(0\right)& =& \frac{1-\beta }{ABC\left(\beta \right)}\mathcal{C}(t_b,v\left(t_b\right))+\frac{\beta }{ABC\left(\beta \right)\times {\rm Y}\left(\beta \right)}\hfill \\ & & \times {\int }_0^{t_{l+1}}\mathcal{C}(\Phi ,v(\Phi )){(t_{b+1}-\Phi )}^{\Phi -1}d\Phi ,\hfill \\ & =& \frac{1-\Phi }{ABC(\Phi )}\mathcal{C}(t_b,v\left(t_b\right))+\frac{\beta }{ABC\left(\beta \right)\times {\rm Y}\left(\beta \right)}\hfill \\ & & \times \sum _{j=0}^l{\int }_{t_j}^{t_{j+1}}\mathcal{C}(\Phi ,v(\Phi )){(t_{b+1}-\Phi )}^{\beta -1}d\beta .\hfill \end{array}$$

(21)

To obtain the estimated value in this case, we utilize the interpolation polynomial on the function $\mathcal{C}(\Phi ,v(\Phi \left)\right)$ on the interval $[t_j,t_{j+1}]$ as follows:

$$\begin{array}{ccc}\hfill \mathcal{C}(\Phi ,z(\Phi \left)\right)& \cong & \mathcal{C}(t_j,z\left(t_j\right))\frac{(t-t_{j-1})}{\mathrm{f}}-\mathcal{C}(t_{j-1},v\left(t_{j-1}\right))\frac{(t-t_j)}{\mathrm{f}},\hfill \end{array}$$

(22)

incorporating this into (21), we have:

$$\begin{array}{c}v\left(t_{b+1}\right)=v\left(0\right)+\frac{1-\beta }{ABC\left(\beta \right)}\mathcal{C}(t_l,v\left(t_l\right))+\frac{\beta }{ABC\left(\beta \right)\times {\rm Y}\left(\beta \right)}\hfill \\ \times \sum _{j=0}^l(\frac{\mathcal{C}(t_j,v\left(t_j\right))}{\mathrm{f}}{\int }_{t_j}^{t_{j+1}}(t-t_{j-1}){(t_{b+1}-t)}^{\beta -1}dt\\ \hfill -\frac{\mathcal{C}(t_{j-1},v\left(t_{j-1}\right))}{\mathrm{f}}{\int }_{t_j}^{t_{j+1}}(t-t_j){(t_{l+1}-t)}^{\beta -1}dt).\end{array}$$

(23)

The results of the aforementioned calculation are as follows:

$$\begin{array}{c}v\left(t_{l+1}\right)=v\left(t_0\right)+\frac{1-\beta }{ABC\left(\beta \right)}\mathcal{C}(t_b,e\left(t_b\right))+\frac{\beta }{ABC\left(\beta \right)}\sum _{j=0}^l[\frac{{\mathrm{f}}^{\beta }\mathcal{C}(t_j,v\left(t_j\right))}{{\rm Y}(\beta +2)}({(b+1-j)}^{\beta }\hfill \\ \times (b-j+2+\beta )-{(b-j)}^{\beta }(b-j+2+2\beta ))\\ \hfill -\frac{{\mathrm{f}}^{\beta }\mathcal{C}(t_{j-1},v\left(t_{j-1}\right))}{{\rm Y}(\beta +2)}({(b+1-j)}^{\beta +1}-{(\beta -j)}^{\beta }(b-j+1+\beta ))].\end{array}$$

(24)

We arrive at the following approximation to the proposed model of HIV infection:

$$\begin{array}{cc}\hfill {\mathit{U}}_v\left(t_{b+1}\right)& ={\mathit{U}}_v\left(t_0\right)+\frac{1-\beta }{ABC\left(\beta \right)}{\mathcal{C}}_1(t_n,v\left(t_b\right))+\frac{\beta }{ABC\left(\beta \right)}\hfill \\ \hfill & \times \sum _{j=0}^b[\frac{{\mathrm{f}}^{\beta }{\mathcal{C}}_1(t_j,v\left(t_j\right))}{{\rm Y}(\beta +2)}\left({(b+1-j)}^{\beta }(b-j+2+\beta )-{(b-j)}^{\beta }(b-j+2+2\beta )\right)\hfill \\ \hfill & -\frac{{\mathrm{f}}^{\beta }{\mathcal{C}}_1(t_{j-1},v\left(t_{j-1}\right))}{{\rm Y}(\beta +2)}({(b+1-j)}^{\beta +1}-{(b-j)}^{\beta }(b-j+1+\beta ))],\hfill \\ \hfill {\mathit{I}}_v\left(t_{b+1}\right)& ={\mathit{I}}_v\left(t_0\right)+\frac{1-\beta }{ABC\left(\beta \right)}{\mathcal{C}}_1(b_b,v\left(t_b\right))+\frac{\beta }{ABC\left(\beta \right)}\sum _{j=0}^b\left[\frac{{\mathrm{f}}^{\beta }{\mathcal{C}}_1(t_j,v\left(t_j\right))}{{\rm Y}(\beta +2)}\right({(b+1-j)}^{\beta }\hfill \\ \hfill & (b-j+2+\beta )-{(b-j)}^{\beta }(b-j+2+2\beta ))-\frac{{\mathrm{h}}^{\beta }{\mathcal{C}}_1(t_{j-1},v\left(t_{j-1}\right))}{{\rm Y}(\beta +2)}{((b+1-j)}^{\beta +1}\hfill \\ \hfill & -{(b-j)}^{\beta }(b-j+1+\beta )\left)\right],\hfill \\ \hfill {\mathit{W}}_v\left(t_{b+1}\right)& ={\mathit{W}}_v\left(t_0\right)+\frac{1-\beta }{ABC\left(\beta \right)}{\mathcal{C}}_1(t_b,v\left(t_b\right))+\frac{\beta }{ABC\left(\beta \right)}\sum _{j=0}^t\left[\frac{{\mathrm{f}}^{\beta }{\mathcal{C}}_1(t_j,v\left(t_j\right))}{{\rm Y}(\beta +2)}\right({(b+1-j)}^{\beta }\hfill \\ \hfill & (b-j+2+\beta )-{(b-j)}^{\beta }(b-j+2+2\beta ))-\frac{{\mathrm{f}}^{\beta }{\mathcal{C}}_1(t_{j-1},v\left(t_{j-1}\right))}{{\rm Y}(\beta +2)}{((b+1-j)}^{\beta +1}\hfill \\ \hfill & -{(b-j)}^{\beta }(b-j+1+\beta )\left)\right].\hfill \end{array}$$

(25)

This numerical method is employed to determine the ABC fractional system solution. This technique will be used to research the fractional system (19) of HIV infection’s transmission phenomena.

5. Numerical Discussions

The interaction between the HIV virus and the immune system is a dynamic process that involves complex interactions between the virus and various components of the immune system. Researchers are continually working to better understand these interactions and to develop new treatments and therapies to manage HIV infection and improve the health outcomes of those living with the virus. Here, our main concern is to conceptualize the dynamical behavior of the system and to examine the solution pathways with the variation of different input factors of the system numerically. Other analyses of the system will be considered in our future research. In this study, we showed the phase portrait, chaos, and the tracking paths in order to comprehend how various elements influence the system. To do this, we perform numerous simulations to understand the role of input factors in the system.

In Figure 1 and Figure 2, the impact of the fractional parameter on the dynamics of HIV infection has been visualized. The fractional values are taken to be 0.95 and 0.75 in these Figures, respectively. We noticed that this parameter has an attractive impact on the dynamics of the system of the HIV infection. The chaotic nature of the system has been highlighted in Figure 3, Figure 4, Figure 5 and Figure 6. In Figure 3, we conceptualize the chaotic nature of the system by changing the value of the input factor ${\aleph }_T$ from 1.0 to $0.01$ and $\beta =0.65$. The chaos of the system is also highlighted in Figure 4 by assuming the input parameter $r=0.3$, ${\aleph }_T=0.1$, while in Figure 5, we considered $\beta =0.9$, ${\aleph }_T=5.0$, and $r=3.0$. These results are obtained by changing the initial values of the state-variable. The chaos in the system is characterized by sensitive dependence on initial conditions, which means that small differences in initial conditions can lead to large differences in the behavior of the system over time. These graphs show that the system exhibits a strong chaotic phenomena, which can be managed using these input settings. The oscillatory and chaotic characteristics of the suggested system are intimately linked by the non-linearity of the system. Policymakers are urged to utilize the fractional parameter, since it has been shown to have a significant impact on the system. The importance of input factors on the output of the system has been shown in Figure 7, Figure 8 and Figure 9. In these simulations, we conceptualize the tracking path behavior with the variation of input factors. We noticed how these factors affect the solution pathways of the system.

Chaos theory focuses on non-linear processes that are challenging to predict or control. The chaotic nature lead to the system instability, and therefore, it is valuable to find out the critical factors responsible for the chaos of the system. It has been observed that the initial conditions of state variables and the non-linearity of the system are the main reasons of this irregular nature of the system. This is vital because they offer important knowledge on the HIV infection mechanism. The system’s chaotic behavior serves as an example of how even minor disturbances may cause substantial systemic alterations. These disorganized charts illustrate the system’s instability and vulnerability to initial circumstances, which renders it unpredictable. According to our research, the non-linearity of the system considerably affects oscillation and chaos, which are both significantly amplified.

6. Concluding Remarks

Mathematical models for HIV and the immune system can help researchers and clinicians to better understand the dynamics of the disease, predict the impact of interventions, and develop new strategies for the management of the infection. In this study, we constructed a novel model for the in vivo dynamics of HIV infection. The interaction of HIV viruses with both HIV-negative and HIV-positive T-cells was considered in our system. The recommended system is discussed in terms of Liouville–Caputo and Atangana–Baleanu fractional operators for more accurate findings. The solution routes and chaotic behavior of HIV infection are displayed with the help of numerical methods. We performed numerous simulations to mimic the system’s most crucial parameters. The solution pathways of the system has been shown numerically. The fractional parameter is found to be sensitive to the system outputs. We have shown a strong relationship between chaos and oscillatory behavior. The impact of input factors has been conceptualized with the help of numerical results and findings. In a future study, we will enhance the model and demonstrate how delay affects the dynamics of HIV infection. In order to demonstrate the effects of immunization and treatment on the system, the model will also be upgraded to include the effect of the therapeutic approach that can restore immune system health and reduce the risk of opportunistic infections and other complications.