Next Article in Journal
Research on Image Encryption Based on Fractional Seed Chaos Generator and Fractal Theory
Next Article in Special Issue
A New Class of Generalized Fractal and Fractal-Fractional Derivatives with Non-Singular Kernels
Previous Article in Journal
Using Laplace Residual Power Series Method in Solving Coupled Fractional Neutron Diffusion Equations with Delayed Neutrons System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Fractional HIV Model under Proportional Hadamard-Caputo Operators

by
Areej A. Almoneef
1,
Mohamed A. Barakat
2,3,* and
Abd-Allah Hyder
4,5,*
1
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
2
Department of Computer Science, College of Al Wajh, University of Tabuk, Tabuk 71491, Saudi Arabia
3
Department of Mathematics, Faculty of Sciences, Al-Azhar University, Assiut 71524, Egypt
4
Department of Mathematics, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
5
Department of Engineering Mathematics and Physics, Faculty of Engineering, Al-Azhar University, Cairo 71524, Egypt
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2023, 7(3), 220; https://doi.org/10.3390/fractalfract7030220
Submission received: 29 January 2023 / Revised: 22 February 2023 / Accepted: 25 February 2023 / Published: 28 February 2023
(This article belongs to the Special Issue Operators of Fractional Integration and Their Applications)

Abstract

:
Modeling human immunodeficiency virus (HIV) via fractional operators has several benefits over the classical integer-order HIV model. The reason is that the fractional HIV model relies not only on the recent status but also on the former conduct of the model. Thus, we are motivated to introduce and analyze a new fractional HIV model. This article focuses on a novel fractional HIV model under the proportional Hadamard-Caputo fractional operators. The study of this model involves the existence and uniqueness (EU) of its solution and the stability examination. We employ Leray–Schauder nonlinear alternative (L-SNLA) and Banach’s fixed point theorems to analyze the EU results. In addition, for this provided model, we develop several forms of Ulam’s stability findings. As a special case of our results, we give and analyze a new fractional HIV model with Hadamard-Caputo operators. Moreover, by appropriate choice of the fractional parameters, the obtained outcomes are valid for analysis of the fractional HIV models formed by several fractional operators defined in the past literature.

1. Introduction

In recent years, mathematical modeling and simulations have been employed as critical tools for predicting the occurrence and severity of infections, as well as for gaining insight into the infection’s behavior, see [1]. These models are salutary instruments that help us understand the dynamics of the immune response to many pathogens. Human immunodeficiency virus (HIV) infection of helper T-cells (HTCs) with therapy is one of the most important models. A retrovirus, the human immunodeficiency virus (HIV) is the cause of acquired immunodeficiency syndrome (AIDS) [2]. HIV targets HTCs, which are the immune system’s most abundant white blood cells and perform a crucial part in safeguarding the human body from illness. If HIV infects a person, the virus assaults HTCs and attempts to multiply. The HIV life cycle [3] is divided into seven stages:
  • Binding: HTCs have receptors on their surfaces that the virus binds to.
  • Fusion: HIV initiates the fusion of its envelope with the helper T-cell (HTC) membrane once it connects to receptors on HTCs. This stage allows the virus to enter the cell.
  • Reverse transcription: The process of transferring genetic information in the form of RNA into DNA through the use of reverse transcription enzyme is known as reverse transcription. HIV can enter the nucleus of HTCs at this stage.
  • Integration: HIV then releases another enzyme called integrase into the nucleus of the HTC after the reverse transcription process. This enzyme is used by the virus to join its DNA with the HTC’s DNA. The virus is still considered inactive at this time, and even sensitive laboratory tests have difficulty detecting it.
  • Replication: HIV may now use the machinery of HTCs to produce viral proteins since it has been integrated into their DNA. It can also manufacture more of its genetic material (RNA) during this time. These two factors make it possible for it to produce additional virus particles.
  • Assembly: Fresh HIV RNA and proteins are sent to the edge of the HTC during the assembly stage, where they mature into immature HIV. In their current state, these viruses are not infectious.
  • Budding: Immature viruses push out of HTCs during the budding stage. They next release a protease that changes the virus’ proteins and turns them into a mature and infectious form.
Antiretroviral therapy (ART) is defined as the use of HIV medications to remedy HIV infection and safeguard the immune system by preventing the virus from growing at specific phases of HIV’s cycle [4].
On the other hand, researchers have used differential equations to solve real-world problems in a variety of fields. For applied researchers, the differential equations theory is well understood. It is well-known in a variety of sectors and is applied in various areas. Various researchers have produced fractional-order differential equations throughout the previous three centuries. They have been widely used to solve real-world problems in a variety of fields; see [5]. Existence theory, stability analysis, numerical and optimization approaches are some of the well-known features that have lately been established. The fractional derivatives presented in fractional calculus were introduced by mathematician researchers. Utilizing fractional operators and applying them to the HIV paradigm has resulted in a variety of research findings. Using numerical simulations, Ding and Ye [6] investigated the stability of balance for a fractional model for HIV infection of HTCs in 2009. The proposed numerical system is intended to be a rough approximation.
Arafa et al. [7] constructed a fractional HIV model and studied the dynamics of HTCs and HIV during primary infection in 2012. This model was solved numerically using the generalized Euler method (GEM). Arafa et al. [8] used the GEM to study and assess the solution attitude for a fractional model of HIV with HTCs in the Caputo style. Based on a Caputo fractional derivative, Arshad et al. [9] investigated the impact of an HIV model on the HTC population in 2017. The numerical approach has been proposed using a finite difference approximation. Khan et al. [10] established the fractional HIV/AIDS model’s stability results and numerical solutions in terms of an Atangana–Baleanu derivative operator. A Caputo fractional derivative HIV model was introduced by Ferrari et al. [11]; it indicated the presence of a reverse transcriptase inhibitor. They proved the model’s EU, and its stability. Using the Mittag–Leffler kernel, Khan et al. [12] examined and evaluated the existence, numerical solutions, and stability of the HIV-TB fractional model. Nazir et al. [13] investigated a model of HIV using the Caputo–Fabrizio fractional operator. Using fixed-point theory, they were able to determine the existence conditions for the solutions. In Hyers–Ulam style, the stability of the relevant solution is likewise demonstrated. The Adomian decomposition method and the Laplace transform approach are used to confirm the model’s approximate solution.
Kongson et al. [14] recently used prominent fixed-point theorems, including Banach and Larey’s nonlinear Schauder’s alternative to investigate the EU of solutions for the fractional HIV model. Furthermore, stability analysis is explored in the perspective of distinct Ulam’s stability. Eventually, they employed a predictor-corrector approach given in [15] to discover approximate fractional HIV model solutions for various fractional derivative orders. Moreover, Albalawi et al. [16] looked at Banach fixed-point theory to analyze EU for the solution of a time-fractional Emden–Fowler-type model and approximation of symmetric solutions utilizing the rational homotopy perturbation method were examined. On the other hand, a numerical investigation of the path-tracking damped oscillatory behavior of a model for the HIV infection of CD4+ T cells was introduced by Shah et al. [17]. A further numerical solution for the HIV model is provided [18]. Readers can view examples of numerical approaches for more information see [19,20].
The motivation of the existing work is to offer and investigate a novel fractional HIV model under the proportional Hadamard-Caputo fractional operators. The EU characteristics for the solution of this model are shown by Banach and Leray’s nonlinear Schauder’s alternative fixed-point forms. Further, different kinds of Ulam’s stability are investigated for the suggested HIV fractional model. As a particular case of our results, we provide and analyze a new fractional HIV model with Hadamard-Caputo operators. Due to the appropriate choice of the fractional parameters, the obtained results are applicable for analyzing some fractional HIV models created by a number of fractional operators in earlier literature.
Our study is divided into different categories, the first of which deals with the introduction of fundamental ideas, notations, and the fractional HIV model associated with the generalized proportional Hadamard-Caputo fractional operator. The second is concerned with examining the resultant model’s EU. The third section defines a few distinct sorts of stability assumptions and uses these assumptions to examine the stability of our model. The final portion concludes final observations and recommendations.

2. Fundamental Instruments

This section covers the basic concepts, essential facts, and notations for generalized fractional derivative and integral operators that will be useful during the rest of our study.
The space of all absolutely continuous functions ζ is referred to as C s [ 1 , T ] . Such that ( s 1 ) -derivative of ζ is absolutely continuous on [1,T].
Definition 1
([5]). If ζ : [ 1 , ) R is a continuous and integrable function and for each k > 0 , then the one-side fractional integral of Hadamard of k t h order is obtained by
k ζ ( ϱ ) = 1 Γ ( k ) 1 ϱ ( ln ϱ ln v ) k 1 ζ ( v ) v d v , k R , k > 0 .
Definition 2
([5]). The Hadamard fractional derivative of k t h order is defined by
D k ζ ( ϱ ) = ( ϱ d d ϱ ) Γ ( n k ) 1 ϱ ( ln ϱ ln v ) n k 1 ζ ( v ) v d v ,
where n = [ k ] + 1 and k > 0 .
Definition 3
([5]). The fractional HCD of order k is obtained by
D ¯ k ζ ( ϱ ) = 1 Γ ( n k ) 1 ϱ ( ln ϱ ln v ) n k 1 [ v d d v ζ ( v ) ] d v v ,
Recently, Barakat et al. [21] employed some highly fascinating fractional derivative and integral fractional operators known as generalized proportional Hadamard-Caputo derivative and integral, respectively, and described by the following definitions.
Definition 4
([21]). The proportional Hadamard-Caputo fractional integral of k order is defined by
P H C J k , σ ζ ( ϱ ) = 1 σ k Γ ( k ) 1 ϱ e σ 1 σ ln ( ϱ v ) ( ln ϱ ln v ) k 1 ζ ( v ) v d v ,
where k > 0 and σ ( 0 , 1 ] .
Definition 5
([21]). The one-side fractional derivative of order k of the proportional Hadamard-Caputo fractional derivative is given by
P H C D k , σ ζ ( ϱ ) = 1 σ l k Γ ( l k ) 1 ϱ e σ 1 σ ln ( ϱ v ) ( ln ϱ ln v ) l k 1 [ D l , σ ζ ( v ) ] d v v ,
where k > 0 , l = [ k ] + 1 , σ ( 0 , 1 ] .
On the other hand, the following important lemmas for our anticipated consequences were recently demonstrated in [21,22].
Lemma 1
([22]). Let r , k C such that R e ( r ) > 0 and R e ( k ) > 0 . Then, for any σ ( 0 , 1 ] we have
  • P H C J k , σ e σ 1 σ ln u ( ln u ) r 1 ( ϱ ) = Γ ( r ) σ k Γ ( r + k ) e σ 1 σ ln ϱ ( ln ϱ ) r + k 1 .
  • P H C D k , σ e σ 1 σ ln u ( ln u ) r 1 ( ϱ ) = σ k Γ ( r ) Γ ( r k ) e σ 1 σ ln ϱ ( ln ϱ ) r k 1 .
Lemma 2
([21]). Assume that all absolutely continuous functions ζ whose derivative of order ( l 1 ) is absolutely continuous on [1,T] are referred to as the set AC l [ 1 , T ] . Let k be a complex number with R e ( k ) > 0 , σ ( 0 , 1 ] ,   l = [ R e ( k ) ] + 1 , ζ L 1 [ 1 , T ] and ( P H C J k , σ ζ ) ( ϱ ) AC l [ 1 , T ] . Then
P H C J k , σ P H C D k , σ ζ ( ϱ ) = ζ ( ϱ ) j = 1 l D k j , σ ζ ( 1 ) Γ ( k j + 1 ) σ k j e σ 1 σ ln ϱ ( ln ϱ ) k j .
Also, the fractional differential equation
P H C D k , σ ζ ( ϱ ) = 0 ,
has a solution
ζ ( ϱ ) = e σ 1 σ ln ϱ s = 0 l 1 M s ( ln ϱ ) s ,
where M s = D s , σ ζ ( 1 ) Γ ( s + 1 ) σ s .

Model Description

Reverse transcription (RT) inhibition occurs prior to the infected cell producing virus, according to a recent model described by Kongson et al. [14]. They looked at the dynamics of the viral population, infected HTCs, and uninfected HTCs. They also divided the infected HTCs into two classes: pre-RT, the stage of affected cells without RT, and post-RT, the stage of affected cells with RT. The fractional HIV model under generalized Caputo derivative is shown as follows
C D k , σ V ( t ) = Ω θ M ( t ) V ( t ) ν 1 V ( t ) + ( β γ + h ) P ( t ) , C D k , σ P ( t ) = θ M ( t ) V ( t ) ( ν 2 + γ + h ) P ( t ) , C D k , σ C ( t ) = ( 1 β ) γ P ( t ) ν 3 C ( t ) , C D k , σ M ( t ) = T ν 3 C ( t ) ϕ M ( t ) ,
where C D k , σ ( . ) represents the generalized Caputo fractional derivative, V ( t ) is the number of HTCs that are susceptible, D ( t ) denotes the number of viral HTCs before RT, C ( t ) represents the concentration of HTCs that are infected and have finished RT, making them able to generate virus, and M ( t ) is the number of viruses. Moreover, consideration was given to the next positive numbers: Ω is the influx rate of HTCs, θ is the infection rate caused by HTC interactions, ν 1 is the proportion of HTCs that naturally die, β is the RT inhibitor’s effectiveness ( β ( 0 , 1 ) ) , γ is the typical rate of HTC conversion from pre-RT to post-RT, h is the rate of reverse transcription failure-induced uninfected cell reversion in infected cells, ν 2 is the infection-related HTC mortality rate, ν 3 is the percentage of HTCs that die after becoming infected, T is the overall quantity of viral particles generated by infected HTCs, and ϕ represents the virus’s elimination rate.
In this paper, we look at the preceding HIV model (8) in terms of the generalized proportional Hadamard-Caputo fractional derivative.
P H C D k , σ V ( t ) = Λ 1 ( t , V , P , C , M ) , P H C D k , σ P ( t ) = Λ 2 ( t , V , P , C , M ) , P H C D k , σ C ( t ) = Λ 3 ( t , V , P , C , M ) , P H C D k , σ M ( t ) = Λ 4 ( t , V , P , C , M ) .
Such that for each t [ 1 , T ] the functions Λ 1 , Λ 2 , Λ 3 , Λ 4 are nonlinearly defined as follows:
Λ 1 ( t , V , P , C , M ) = Ω θ M ( t ) V ( t ) ν 1 V ( t ) + ( β γ + h ) P ( t ) , Λ 2 ( t , V , P , C , M ) = θ M ( t ) V ( t ) ( ν 2 + γ + h ) P ( t ) , Λ 3 ( t , V , P , C , M ) = ( 1 β ) γ P ( t ) ν 3 C ( t ) , Λ 4 ( t , V , P , C , M ) = T ν 3 C ( t ) ϕ M ( t ) .
Having the initial conditions V ( 1 ) = V 0 , P ( 1 ) = P 0 , C ( 1 ) = C 0 , M ( 1 ) = M 0 .

3. Existence and Uniqueness

It is critical that we analyze the model to ensure the existence of its solution. Various ideas and tools can be used to derive this concept. Fixed-point theory is a strong tool for studying the aforementioned requirement. As a result, we use the fixed point theorem to prove the existence of at least one solution and its uniqueness for model (9) in this section of the manuscript.
Take A = C k ( [ 1 , T ] , R ) be the Banach space of absolutely continuous functions ζ from [ 1 , T ] into R with the norm:
| | ζ | | = sup t [ 1 , T ] | ζ ( t ) | , | ζ ( t ) | = | V ( t ) | + | P ( t ) | + | C ( t ) | + | M ( t ) | .
Now, we modify the HIV model (9) to include the following assumptions:
ζ ( t ) = V ( t ) , P ( t ) , C ( t ) , M ( t ) . ζ 0 = V 0 , P 0 , C 0 , M 0 . Λ ( t , ζ ( t ) ) = Λ 1 ( t , V , P , C , M ) , Λ 2 ( t , V , P , C , M ) , Λ 3 ( t , V , P , C , M ) , Λ 4 ( t , V , P , C , M ) ,
and
Λ ( 1 , ζ ( 1 ) ) = Λ 1 ( 1 , V ( 1 ) , P ( 1 ) , C ( 1 ) , M ( 1 ) ) , Λ 2 ( 1 , V ( 1 ) , P ( 1 ) , C ( 1 ) , M ( 1 ) ) , Λ 3 ( 1 , V ( 1 ) , P ( 1 ) , C ( 1 ) , M ( 1 ) ) , Λ 4 ( 1 , V ( 1 ) , P ( 1 ) , C ( 1 ) , M ( 1 ) ) .
By using Lemma 2, (11) and (12). The HIV model (9) can be expressed by the next system
P H C D k , σ ζ ( t ) = Λ ( t , ζ ( t ) ) , ζ ( 1 ) = ζ 0 .
Which is the same as the integral operator that follows
ζ ( t ) = ζ 0 + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v .
Define the operator T : A A such that
( T ζ ) ( t ) = ζ 0 + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v .
The operator T has fixed points if and only if the problem (13), which is analogous to the model (9), has solutions.
For the purposes of further investigation, we assert that the following assumptions are true:
( G 1 )
There is a constant δ Λ > 0 , such that
| | Λ ( t , ζ 1 ( t ) ) Λ ( t , ζ 2 ( t ) ) | | δ Λ | | ζ 1 ( t ) ζ 2 ( t ) | |
( G 2 )
There is a nondecreasing continuous function X : [ 0 , ) [ 0 , ) which satisfies each μ 1 X ( μ ζ ) μ X ( ζ ) , and a function ψ A such that
| | Λ ( t , ζ ( t ) ) | | ψ ( t ) X ( | ζ ( t ) | ) .
( G 3 )
There are some constants ϵ , B > 0 such that
B | | ζ 0 | | + sup t [ 1 , T ] ψ ( t ) X ( ϵ ) ( ln T ) k σ k Γ ( k + 1 ) > 1 .
The following theorem gives the existence of the solution for the HIV model (9) by using the nonlinear Leray–Schauder alternative fixed-point theory.
Theorem 1.
If assumptions ( G 2 ) and ( G 3 ) are fulfilled. Then the HIV model (9) has at least one solution on [ 1 , T ] .
Proof. 
Chose ϵ > 0 such that [ B ] r = { ζ A ; | | ζ | | ϵ } . Utilising ( G 2 ) for each t [ 1 , T ] we get
| ( T ζ ) ( t ) | | | ζ 0 | | + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | Λ ( t , ζ ( t ) ) | | v d v | | ζ 0 | | + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 ψ ( t ) X ( | ζ ( t ) | ) v d v | | ζ 0 | | + sup t [ 1 , T ] ψ ( t ) X ( | ζ ( ϵ ) | ) σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 d v v .
Making use of the fact that e σ 1 σ ln ( t v ) 1 , since σ ( 0 , 1 ] , and 1 v < t T we have
| | T ζ | | | | ζ 0 | | + sup t [ 1 , T ] ψ ( t ) X ( | ζ ( ϵ ) | ) σ k Γ ( k + 1 ) ( ln T ) k .
Hence, bounded balls are mapped into bounded sets in A via the operator T .
On the other hand, we show that T is completely continuous. To do that, choose ζ A and 1 v 1 < v 2 . Such that
| ( T ζ ) ( v 2 ) ( T ζ ) ( v 1 ) | | 1 σ k Γ ( k ) 1 v 2 e σ 1 σ ln ( v 2 v ) ( ln v 2 ln v ) k 1 Λ ( v , ζ ( v ) ) v d v 1 σ k Γ ( k ) 1 v 1 e σ 1 σ ln ( v 1 v ) ( ln v 1 ln v ) k 1 Λ ( v , ζ ( v ) ) v d v | 1 σ k Γ ( k ) [ 1 v 1 e σ 1 σ ln ( v 2 v ) ( ln v 2 ln v ) k 1 | | Λ ( v , ζ ( v ) ) | | v d v + v 1 v 2 e σ 1 σ ln ( v 2 v ) ( ln v 2 ln v ) k 1 | | Λ ( v , ζ ( v ) ) | | v d v 1 σ k Γ ( k ) 1 v 1 e σ 1 σ ln ( v 1 v ) ( ln v 1 ln v ) k 1 | | Λ ( v , ζ ( v ) ) | | v d v ] sup t [ 1 , T ] ψ ( t ) X ( | ζ ( ϵ ) | ) σ k Γ ( k ) 1 v 1 e σ 1 σ ln ( v 2 v ) ( ln v 2 ln v ) k 1 e σ 1 σ ln ( v 1 v ) ( ln v 1 ln v ) k 1 d v v + sup t [ 1 , T ] ψ ( t ) X ( | ζ ( ϵ ) | ) σ k Γ ( k + 1 ) ( ln v 2 ln v 1 ) k 0 as v 2 v 1 .
Therefore, the operator T is completely continuous, according to the Arzel’a–Ascoli theorem. Finally, we prove that the collection of all solutions that satisfy ζ = ρ T ζ for every 0 < ρ < 1 is bounded. Let ζ A be a solution, then for each 1 < t < T we have
| ζ ( T ) | = ρ ( T ζ ) ( t ) | | ζ 0 | | + sup t [ 1 , T ] ψ ( t ) X ( | ζ ( ϵ ) | ) σ k Γ ( k + 1 ) ( ln T ) k .
Hence,
| | ζ | | | | ζ 0 | | + sup t [ 1 , T ] ψ ( t ) X ( | ζ ( ϵ ) | ) σ k Γ ( k + 1 ) ( ln T ) k .
By using the assumption G 3 , B > 0 with | | ζ | | B . Set Q : = { ζ A ; | | ζ | | < B } . It is worth noting that T : Q ¯ A is completely continuous. By the definition of Q ¯ no ζ Q ¯ satisfies | ζ ( T ) | = ρ T ζ . As a result of the (L-SNLA), we can deduce that (9) has a solution on [ 1 , T ] .
We will now show how the Banach fixed-point principle can be used to prove EU of the solution for the HIV model (9).
Theorem 2.
If assumption ( G 1 ) is fulfilled for each ζ 1 , ζ 2 A , t [ 1 , T ] and if
δ Λ ( ln T ) k σ σ k Γ ( k + 1 ) < 1 .
Then, for the HIV model (9) there is a unique solution on [ 1 , T ] .
Proof. 
Let sup v [ 1 , T ] Λ ( v , 0 ) = P ¯ < . Choose B ρ 1 = ζ A : ζ ρ 1 , where
ζ 0 + P ¯ ( ln T ) k σ σ k Γ ( k + 1 ) ρ 1 1 P ¯ ( ln T ) k σ σ k Γ ( k + 1 ) .
we conclude that B ρ 1 is closed, bounded, and convexly set in A .
First, we will show that T B ρ 1 B ρ 1 . Letting ζ B ρ 1 , we get
| ( T ζ ) ( v ) | ζ 0 + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | Λ ( t , ζ ( t ) ) | | v d v ζ 0 + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | Λ ( v , ζ ( v ) ) Λ ( v , 0 ) | | + | | Λ ( v , 0 ) | | v d v ζ 0 + ρ 1 δ Λ + P ¯ σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 d v v ζ 0 + ρ 1 δ Λ + P ¯ σ k Γ ( k ) 1 t ( ln t ln v ) k 1 d v v ζ 0 + ( ρ 1 δ Λ + P ¯ ) ( ln T ) k σ k Γ ( k + 1 ) ρ 1 .
Hence, T B ρ 1 B ρ 1 .
Secondly, we shall prove that T is a contraction operator. For each ζ 1 , ζ 2 A , t [ 1 , T ] , we get
T ζ 1 ( t ) T ζ 2 ( t ) 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | Λ ( v , ζ 1 ( v ) ) Λ ( v , ζ 2 ( v ) ) | | v d v δ Λ σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 ζ 1 ( t ) ζ 2 ( t ) d u v δ Λ ( ln T ) k σ k Γ ( k + 1 ) ζ 1 ζ 2 .
Since δ Λ ( ln T ) k σ k Γ ( k + 1 ) < 1 , there exist ζ B ρ 1 such that ζ = T ζ . Hence, the problem (13), which is analogous to the model (9), has solutions.
Finally, we will demonstrate the fixed point’s uniqueness. Assume that ζ 1 , ζ 2 are two separate fixed points and from (19) we get
| | ζ 1 ζ 2 | | = T ζ 1 ( t ) T ζ 2 ( t ) < | | ζ 1 ζ 2 | | .
This is contradictory. Hence, a fixed point ζ is a unique solution for problem (15). Accordingly, the model (9) has a unique solution on [ 1 , T ] . □
Remark 1.
Hadmard-Caputo integral and differential operators can be obtained by setting σ = 1 in (4) and (5), respectively. These operators are applied to model (9) to produce the Hadamard-Caputo HIV model.
The EU of the Hadamard-Caputo HIV model are given by the following corollaries, which use the nonlinear Leray–Schauder alternative fixed-point theory and the Banach contraction principle, respectively.
Corollary 1.
If assumptions ( G 2 ) and ( G 3 ) are fulfilled. Then, the Hadamard-Caputo HIV model has at least one solution on [ 1 , T ] , where ( G 3 ) is defined as if there are some constants ϵ , B > 0 , such that
B | | ζ 0 | | + sup t [ 1 , T ] ψ ( t ) X ( ϵ ) ( ln T ) k Γ ( k + 1 ) > 1 .
Corollary 2.
If assumption ( G 1 ) is fulfilled for each ζ 1 , ζ 2 A , t [ 1 , T ] and if
δ Λ ( ln T ) k Γ ( k + 1 ) < 1 .
Then, for the Hadamard-Caputo HIV model there is a unique solution on [ 1 , T ] .

4. Results on Ulam–Hyers Stability (UHS)

In this section, we provide some necessary requirements for model (9) to satisfy the assumptions of different types of stability. UHS, extended UHS, Ulam–Hyers–Rassias stability (UHRS), and extended UHRS are examples of these types of stability. Before we show the stability theorems, we need to define a few terms.
Let us take a real positive constant λ and assume that Ψ Λ : [ 1 , T ] R + is a continuous function. The inequalities that will be used to define stability are as follows
P H C D k , σ ζ ( t ) Λ ( t , ζ ( t ) ) λ , t [ 1 , T ] ,
P H C D k , σ ζ ( t ) Λ ( t , ζ ( t ) ) λ Ψ Λ ( t ) , t [ 1 , T ] ,
P H C D k , σ ζ ( t ) Λ ( t , ζ ( t ) ) Ψ Λ ( t ) , t [ 1 , T ] .
Definition 6
([23]). The Equation (13) is seen as being stable under the UH condition if ∃ a solution ζ A for the Equation (13) moreover, there are C Λ > 0 , such that for each x A fulfilling the inequality (20) and for all λ > 0 we have
x ( t ) ζ ( t ) C Λ λ , t [ 1 , T ] ,
where C Λ = max C Λ 1 , C Λ 2 , C Λ 3 , C Λ 4 .
Definition 7
([23]). If there exists Ψ Λ C [ 1 , T ] , with Ψ Λ ( 1 ) = 0 , and for each ζ A fulfilling the inequality (21), then ∃ a solution x A for (13) with
ζ ( t ) x ( t ) Ψ Λ ( λ ) , t [ 1 , T ] ,
where Ψ Λ = max Ψ Λ 1 , Ψ Λ 2 , Ψ Λ 3 , Ψ Λ 4 . Then, the extended UH condition on the problem (13) is said to be stable.
A fundamental attribute that can be used to achieve UHS and extended UHS is now provided.
Lemma 3.
Suppose that k > 0 and σ ( 0 , 1 ] . The subsequent inequality
ζ ( t ) ζ 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v λ ( ln T ) k σ k Γ ( k + 1 ) .
is satisfied if ζ A is a solution of the inequality (20).
Proof. 
Since ζ satisfies (20), then ∃ a function y A (dependent on ζ ) satisfies
y ( t ) λ , y = max ( y 1 , y 2 , y 3 , y 4 ) , t [ 1 , T ] .
And
P H C D k , σ ζ ( t ) = Λ ( t , ζ ( t ) ) + y ( t ) , t [ 1 , T ] , ζ ( 1 ) = ζ 0 0 .
Lemma 2 which states that the solution to problem (27) can be expressed as
ζ ( t ) = ζ 0 + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 y ( t ) v d v .
By using both Equations (28) and (26), we have
ζ ( t ) ζ 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | y ( t ) | | v d v λ ( ln T ) k σ k Γ ( k + 1 ) .
Thus, inequality (25) is established. □
We are now prepared to demonstrate the UHS plus extended UHS.
Theorem 3.
Consider Λ ( t , ζ ( t ) ) C ( [ 1 , T ] , R + ) for each ζ A . Model (9) or problem (13) is stable under UHS and the extended UHS conditions if both ( G 3 ) and Equation (16) are all met.
Proof. 
Consider that x A is a unique solution of the problem (13). Further, consider that λ > 0 and ζ A satisfies the inequality (20). Utilizing Lemma 3 and Equation (13), we obtain
ζ ( t ) x ( t ) ζ ( t ) x 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , x ( t ) ) v d v ζ ( t ) ζ 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | Λ ( t , ζ ( t ) ) Λ ( t , x ( t ) ) | | v d v ζ ( t ) ζ 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v + δ Λ σ k Γ ( k ) 1 t ( ln t ln v ) k 1 | | ζ ( t ) x ( t ) | | v d v λ ( ln T ) k σ k Γ ( k + 1 ) + δ Λ ( ln T ) k σ k Γ ( k + 1 ) ζ ( t ) x ( t ) .
Hence, ζ ( t ) x ( t ) λ C Λ , where
C Λ = ( ln T ) k σ k Γ ( k + 1 ) δ Λ ( ln T ) k .
Therefore, model (9) is (UH)-stable. Additionally, it is implied that model (9) is extended (UH)-stable by setting λ C Λ = Ψ Λ ( λ ) such that Ψ Λ ( 1 ) = 0 . The proof is now complete. □
The following significant result, which provides UHS for the Hadamard-Caputo HIV model, is obtained by using Corollary 1 and Theorem 3.
Corollary 3.
Consider Λ ( t , ζ ( t ) ) C ( [ 1 , T ] , R + ) for each ζ A . Then, the Hadamard-Caputo HIV model is stable under (UH) and the extended (UH) conditions if both ( G 3 ) and Equation (16) with respect to fractional Hadamard-Caputo operators are all met.
Finally, the (UHRS) and extended UHRS assumptions can be realized by model (9) if certain requirements are met. The following characteristics of these stability categories can be seen as follows.
Definition 8
([23]). The Equation (13) is seen as being stable under UHR condition with respect to Ψ Λ , if there exists a solution ζ A for the Equation (13). Moreover, there are M Ψ Λ > 0 , such that for each x A fulfilling the inequality (20) and for all λ > 0 we have
x ( t ) ζ ( t ) λ M Ψ Λ Ψ Λ ( t ) , t [ 1 , T ] .
Definition 9
([23]). The Equation (13) is seen as being stable under extended UHR condition with respect to Ψ Λ , if there exists a solution ζ A for the Equation (13). Moreover, there are M Ψ Λ > 0 , such that for each x A fulfilling the inequality (20) we have
x ( t ) ζ ( t ) M Ψ Λ . Ψ Λ ( t ) , t [ 1 , T ] .
A fundamental lemma that will be explained will be used to demonstrate the UHRS and generalized UHRS result.
Lemma 4.
Consider the following condition:
( G 4 )
 ∃ a non-decreasing operator Ψ Λ A and η Ψ Λ > 0 , such that the following inequality holds:
P H C J k , σ Ψ Λ ( t ) η Ψ Λ Ψ Λ ( t ) t [ 1 , T ] .
If k > 0 and σ ( 0 , 1 ] , and the subsequent inequality
ζ ( t ) ζ 0 P H C J k , σ Λ ( t , ζ ( t ) ) λ η Ψ Λ Ψ Λ ( t )
is satisfied for each ζ A which is a solution of inequality (22).
Proof. 
Since ζ fulfills (21), then ∃ an operator V A (dependent on ζ ) which satisfies
V ( t ) λ Ψ Λ ( t ) , V = max ( V 1 , V 2 , V 3 , V 4 ) , t [ 1 , T ] .
And
P H C D k , σ ζ ( t ) = Λ ( t , ζ ( t ) ) + V ( t ) , t [ 1 , T ] , ζ ( 1 ) = ζ 0 0 .
Hence, the solution to problem (37) can be expressed as
ζ ( t ) = ζ 0 + P H C J k , σ Λ ( t , ζ ( t ) ) + P H C J k , σ V ( t ) .
By using both Equations (36) and (38), we have
ζ ( t ) ζ 0 P H C J k , σ Λ ( t , ζ ( t ) ) P H C J k , σ | | V ( t ) | | λ η Ψ Λ Ψ Λ ( t ) ,
where η Ψ Λ = ( ln T ) k σ k Γ ( k + 1 ) . Thus, inequality (41) is established. □
Lemma 5.
Let σ = 1 in Equations (4) and (5), and consider the following condition:
( G 4 ) :
a non-decreasing operator Ψ Λ A and ∃ η Ψ Λ > 0 , satisfies:
H C J k Ψ Λ ( t ) η Ψ Λ Ψ Λ ( t ) t [ 1 , T ] .
If k > 0 and σ ( 0 , 1 ] , and the subsequent inequality
ζ ( t ) ζ 0 H C J k Λ ( t , ζ ( t ) ) λ η Ψ Λ Ψ Λ ( t ) .
is satisfied for each ζ A which is a solution of inequality (22) with respect to Hadamard-Caputo operators.
Proof. 
The proof can be concluded by adding σ = 1 to the proof of lemma 4. □
Lastly, UHRS and extended UHRS will then be shown for the HIV model (9).
Theorem 4.
Consider Λ ( t , ζ ( t ) ) C ( [ 1 , T ] , R + ) for each ζ A . The HIV model (9) or problem (13) is stable under UHR and extended UHR conditions if ( G 3 ) , ( G 4 ) and Equation (16) are all met.
Proof. 
Choose x A is a unique solution of the problem (13). Also, consider λ > 0 and ζ A satisfies the inequality (21). Utilizing Lemma 4 and Equation (13), we obtain
ζ ( t ) x ( t ) ζ ( t ) x 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , x ( t ) ) v d v ζ ( t ) ζ 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v + 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 | | Λ ( t , ζ ( t ) ) Λ ( t , x ( t ) ) | | v d v ζ ( t ) ζ 0 1 σ k Γ ( k ) 1 t e σ 1 σ ln ( t v ) ( ln t ln v ) k 1 Λ ( t , ζ ( t ) ) v d v + δ Λ σ k Γ ( k ) 1 t ( ln t ln v ) k 1 | | ζ ( t ) x ( t ) | | v d v λ η Ψ Λ Ψ Λ ( t ) + δ Λ ( ln T ) k σ k Γ ( k + 1 ) ζ ( t ) x ( t ) .
Hence, ζ ( t ) x ( t ) λ M Ψ Λ Ψ Λ ( t ) , where
M Ψ Λ = σ k Γ ( k + 1 ) η Ψ Λ σ k Γ ( k + 1 ) δ Λ ( ln T ) k .
Therefore, model (9) is UHR stable. Additionally, it is implied that model (9) is extended UHR stable by setting λ = 1 such that Ψ Λ ( 1 ) = 0 . The proof is now complete. □
The Hadamard-Caputo HIV model’s UHRS and extended UHRS are examined in the following result using Corollary 1 and Theorem 4.
Corollary 4.
Consider Λ ( t , ζ ( t ) ) C ( [ 1 , T ] , R + ) for each ζ A . The Hadamard-Caputo HIV model is stable under UHR and extended UHR conditions if ( G 3 ) , ( G 4 ) and Equation (16) with respect to Hadamard-Caputo operators are all met.

5. Conclusions

Evidently, fractional models have many benefits that classical mathematical models do not. Modeling real phenomena via fractional operators involves memory influences in contrast to modeling by integer-order derivative. This makes fractional mathematical models more beneficial because, in most events, the models rely not only on the current status but also on the previous attitude of the model. This paper has justified the validity of employing the generalized proportional Hadamard-Caputo fractional operators to introduce a new fractional HIV model. The solution of this model has been proven to exist and be unique via Banach’s and (L-SNLA) fixed-point theorems. Further, for the suggested fractional HIV model, diverse kinds of Ulam’s stability results have been provided and analyzed. The technique proposed in this study, combined with the generalized proportional Hadamard-Caputo fractional operators, can be utilized to model and analyze several infectious diseases, like measles, strep throat, COVID-19, and salmonella. Moreover, the obtained existence, uniqueness, and stability results are valid to analyze the solution of the fractional HIV under Hadamard-Caputo fractional operators [24], when σ = 1 , and, in the traditional HIV model, when σ = k = 1 [25]. Recently, Hyder et al. [26] derived a fractional HIV model via the improved fractional operators. They investigated the existence, uniqueness, and stability of the solution on the time interval [0,T]. Obviously, the generalized proportional Hadamard-Caputo operators and the improved fractional operators are alternative fractional operators for generalizing the Hadamard-Caputo fractional operators. Hence, in comparing the results of the current study with the results of Hyder et al. [26], one can conclude that both results are different and give alternative qualitative fractional characteristics for the solution of the fractional HIV model.

Author Contributions

Methodology and conceptualization, A.A.A., M.A.B. and A.-A.H.; data curation and writing—original draft, A.A.A., M.A.B. and A.-A.H.; investigation and visualization, A.A.A., M.A.B. and A.-A.H.; validation, writing—reviewing, and editing, A.A.A., M.A.B. and A.-A.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Khalid University under grant (RGP.2/102/44) and Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R337).

Data Availability Statement

The corresponding authors will provide the data used in this work upon reasonable request.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through Research Groups Program under grant (RGP.2/102/44). The authors would like to acknowledge the Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R337), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nowak, M.A.; May, R.M. Virus Dynamics; Oxford University Press: Oxford, UK, 2000. [Google Scholar]
  2. Haas, G.; Hosmalin, A.; Hadida, F.; Duntze, J.; Debré, P.; Autran, B. Dynamics of HIV variants and specific cytotoxic T-cell recognition in nonprogressors and progressors. Immunol. Lett. 1997, 57, 63–68. [Google Scholar] [CrossRef]
  3. Kirchhoff, F. HIV Life Cycle: Overview. In Encyclopedia of AIDS; Hope, T.J., Stevenson, M., Richman, D., Eds.; Springer: New York, NY, USA, 2021; pp. 1–9. [Google Scholar]
  4. Smith, R.J.; Wahl, L.M. Distinct effects of protease and reverse transcriptase inhibition in an immunological model of HIV-1 infection with impulsive drug effects. Bull. Math. Biol. 2004, 66, 1259–1283. [Google Scholar] [CrossRef] [PubMed]
  5. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  6. Ding, Y.; Ye, H. A fractional-order differential equation model of HIV infection of CD4+ T-cells. Math. Comput. Model. 2009, 50, 386–392. [Google Scholar] [CrossRef]
  7. Arafa, A.A.M.; Rida, S.Z.; Khalil, M. Fractional modeling dynamics of HIV and CD4+ T-cells during primary infection. Nonlinear Biomed. Phys. 2012, 6, 1. [Google Scholar] [CrossRef] [Green Version]
  8. Arafa, A.A.M.; Rida, S.Z.; Khalil, M. A fractional–order model of HIV infection with drug therapy effect. J. Egypt. Math. Soc. 2014, 22, 538–543. [Google Scholar] [CrossRef] [Green Version]
  9. Arshad, S.; Baleanu, D.; Bu, W.; Tang, Y. Effects of HIV infection on CD4+ T-cell population based on a fractional-order model. Adv. Differ. Equations 2017, 2017, 92. [Google Scholar] [CrossRef] [Green Version]
  10. Khan, A.; Gómez-Aguilar, J.F.; Saeed Khan, T.; Khan, H. Stability analysis and numerical solutions of fractional order HIV/AIDS model. Chaos Solitons Fractals 2019, 122, 119–128. [Google Scholar] [CrossRef]
  11. Ferrari, A.J.; Santillan Marcus, E.A. Study of a fractional-order model for HIV infection of CD4+ T-Cells with treatment. J. Fract. Calc. Appl. 2020, 11, 12–22. [Google Scholar]
  12. Khan, H.; Gómez-Aguilar, J.F.; Alkhazzan, A.; Khan, A. A fractional order HIV-TB coinfection model with nonsingular Mittag-Leffler Law. Math. Methods Appl. Sci. 2020, 43, 3786–3806. [Google Scholar] [CrossRef]
  13. Nazir, G.; Shah, K.; Debbouche, A.; Khan, R.A. Study of HIV mathematical model under nonsingular kernel type derivative of fractional order. Chaos Solitons Fractals 2020, 139, 110095. [Google Scholar] [CrossRef]
  14. Kongson, J.; Thaiprayoon, C.; Sudsutad, W. Analysis of a fractional model for HIV CD4+ T-cells with treatment under generalized Caputo fractional derivative. AIMS Math. 2021, 6, 7285–7304. [Google Scholar] [CrossRef]
  15. Erturk, V.S.; Kumar, P. Solution of a COVID-19 model via new generalized Caputo-type fractional derivatives. Chaos Solitons Fractals 2020, 139, 110280. [Google Scholar] [CrossRef] [PubMed]
  16. Albalawi, K.S.; Alkahtani, B.S.; Kumar, A.; Goswami, P. Numerical solution of time–fractional Emden–Fowler–type equations using the rational homotopy perturbation method. Symmetry 2023, 15, 258. [Google Scholar] [CrossRef]
  17. Shah, Z.; Jan, R.; Kumam, P.; Deebani, W.; Shutaywi, M. Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo-Fabrizio Derivative. Molecules 2021, 26, 1806. [Google Scholar] [CrossRef]
  18. Alzaid, S.S.; Alkahtani, B.S.T.; Sharma, S.; Dubey, R.S. Numerical solution of fractional model of HIV-1 infection in framework of different fractional derivatives. J. Funct. Spaces 2021, 2021, 6642957. [Google Scholar] [CrossRef]
  19. Shrahili, M.; Dubey, R.S.; Shafay, A. Inclusion of fading memory to Banister model of changes in physical condition. Discret. Contin. Dyn. Syst.–S 2020, 13, 881–888. [Google Scholar] [CrossRef] [Green Version]
  20. Mishra, M.N.; Aljohani, A.F. Mathematical modelling of growth of tumour cells with chemotherapeutic cells by using Yang–Abdel–Cattani fractional derivative operator. J. Taibah Univ. Sci. 2022, 16, 1133–1141. [Google Scholar] [CrossRef]
  21. Barakat, M.A.; Soliman, A.H.; Hyder, A.-A. Langevin Equations with Generalized Proportional Hadamard-Caputo Fractional Derivative. Comput. Intell. Neurosci. 2021, 2021, 6316477. [Google Scholar] [CrossRef]
  22. Jarad, F.; Alqudah, M.A.; Abdeljawad, T. On more general forms of proportional fractional operators. Open Math. 2020, 18, 167–176. [Google Scholar] [CrossRef] [Green Version]
  23. Rus, I.A. Ulam stabilities of ordinary differential equations in a Banach space. Carpathian J. Math. 2010, 26, 103–107. [Google Scholar]
  24. Jarad, F.; Abdeljawad, T.; Baleanu, D. Caputo-type modification of the Hadamard fractional derivatives. Adv. Differ. Equations 2012, 2012, 142. [Google Scholar] [CrossRef] [Green Version]
  25. Arenas, A.J.; González-Parra, G.; Naranjo, J.J.; Cogollo, M.; De La Espriella, N. Mathematical analysis and numerical solution of a model of HIV with a discrete time delay. Mathematics 2021, 9, 257. [Google Scholar] [CrossRef]
  26. Hyder, A.; Barakat, M.A.; Rizk, D.; Shah, R.; Nonlaopon, K. Study of HIV model via recent improved fractional differential and integral operators. AIMS Math. 2023, 8, 1656–1671. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Almoneef, A.A.; Barakat, M.A.; Hyder, A.-A. Analysis of the Fractional HIV Model under Proportional Hadamard-Caputo Operators. Fractal Fract. 2023, 7, 220. https://doi.org/10.3390/fractalfract7030220

AMA Style

Almoneef AA, Barakat MA, Hyder A-A. Analysis of the Fractional HIV Model under Proportional Hadamard-Caputo Operators. Fractal and Fractional. 2023; 7(3):220. https://doi.org/10.3390/fractalfract7030220

Chicago/Turabian Style

Almoneef, Areej A., Mohamed A. Barakat, and Abd-Allah Hyder. 2023. "Analysis of the Fractional HIV Model under Proportional Hadamard-Caputo Operators" Fractal and Fractional 7, no. 3: 220. https://doi.org/10.3390/fractalfract7030220

APA Style

Almoneef, A. A., Barakat, M. A., & Hyder, A. -A. (2023). Analysis of the Fractional HIV Model under Proportional Hadamard-Caputo Operators. Fractal and Fractional, 7(3), 220. https://doi.org/10.3390/fractalfract7030220

Article Metrics

Back to TopTop