Abstract

A classical epidemiological framework is used to qualitatively assess the impact of early detection and treatment on the dynamics of HIV/AIDS. Within this theoretical framework, two classes of infected populations: those infected but unaware of their serological status and those who are aware of their disease status, are considered. In this context, we formulate and analyze a deterministic model for the transmission dynamics of HIV/AIDS and assess the potential population-level impact of early detection in curtailing the epidemic. A critical threshold parameter for which case detection will have a positive impact is derived. Model parameters sensitivity analysis indicates that the number of partners is the most sensitive (in increasing the average number of secondary transmission) parameter. However, the case detection coverage is the main drivers in reducing the initial disease transmission. Numerical simulations of the model are provided to support the analytical results. Early detection and treatment alone are insufficient to eliminate the disease, and other control strategies are to be explored.

1. Introduction

HIV/AIDS has killed more than 25 million people globally since its emergence in 1981, making it one of the most destructive epidemics in recorded history. The disease continues to inflict a significant morbidity, mortality, and social-economic and public health burden. For the estimated 33.3 million people living with HIV after nearly 30 years into a very complex epidemic, the gains are real but still fragile, even as the number of annual AIDS-related deaths worldwide has steadily decreased from the peak of 2.1 million in 2004 to an estimated 1.8 million in 2009 [1].

Various preventative and therapeutic measures have been embarked upon, aiming at combating one of the greatest pandemics in modern times [2, 3]. In sub-Saharan Africa, many infected individuals are unaware of their disease status. Recent randomized control trials have found that treating HIV-positive individuals with antiretroviral drugs reduces the risk of them transmitting the disease to their heterosexual partners by more than 90% [4]. As the current treatment therapy has been proven beyond reasonable doubt to reduce transmission, it is imperative to identify those who are infected and put them on treatment when eligible. It is therefore desirable to encourage voluntary testing that will increase case detection, thereby reducing the number of secondary infections of individuals receiving treatment. There is no explicit mathematical account of the potential population level impact of case detection when treatment is available known to us. Thus, a dynamical system model is formulated in order to assess the trade-off/population-level impact between treatment and early detection of HIV positive. The results are sensitive to parameter values, and for this reason, a deterministic sensitivity analysis is carried out.

The rest of this work is organized as follows. The basic model formulation and its analysis are provided in Section 2. The extended model incorporating case detection and treatment is described and analyzed in Section 3. The model simulation using heuristic parameter values for the purpose of illustration follows in Section 4.

2. Model Formulation and Analysis

We begin by formulating a deterministic sex-structured basic HIV/AIDS model (i.e., without interventions). Individuals are identified as male and female only in connection with features peculiar to their sex. The male to female infectivity rate is greater than the female to male [5]. We also assume that the mixing between individuals is homogeneous; individuals may become HIV-infected only through sexual contacts with HIV infected individuals. Those in the final disease stage are considered too ill to remain sexually active. We ignore important HIV transmission path such as intravenous drug injections, vertical transmission, breast feeding, blood transfusion, and needle sharing. It is also assumed that there is no recruitment of HIV positive. The total heterosexual population is divided into male and female subpopulations with the following epidemiological subgroups (the classification is based on individuals disease status): susceptible male and female (π‘†π‘š,𝑆𝑓), infected male and female (πΌπ‘š,𝐼𝑓), and symptomatic individuals in the final disease (AIDS) stage (π΄π‘š,𝐴𝑓).

New recruits enter the heterosexually active population at constant rates Ξ›π‘š and Λ𝑓 for male and female, respectively (all recruits into the population are assumed susceptible). Male and female susceptible acquire infection at time-dependent rates πœ†π‘š and πœ†π‘“ and become infectious. Infectious individuals exhibit AIDS clinical defining symptoms at rates β„Ž and 𝑧 for male and female, respectively. In the absence of the disease, individuals in the population die of natural death at the rate πœ‡. The disease-induced mortality rate is 𝑑 for both individuals in the infectious and AIDS classes. A full description of the model variables and parameters used in the model is described in Tables 2 and 3, respectively.

Based on our model description and assumptions, we establish the following equations. We note that the red and dash arrows, respectively, in Figures 1 and 2 are the feedback branches which indicate how the male and female subpopulations are coupled (via the force of infection)Maleπ‘‘π‘†π‘šπ‘‘π‘‘=Ξ›π‘šβˆ’ξ€·πœ‡+πœ†π‘šξ€Έπ‘†π‘š,π‘‘πΌπ‘šπ‘‘π‘‘=πœ†π‘šπ‘†π‘šβˆ’(πœ‡+𝑑+β„Ž)πΌπ‘š,π‘‘π΄π‘šπ‘‘π‘‘=β„ŽπΌπ‘šβˆ’(πœ‡+𝑑)π΄π‘š,Female𝑑𝑆𝑓𝑑𝑑=Ξ›π‘“βˆ’ξ€·πœ‡+πœ†π‘“ξ€Έπ‘†π‘“,𝑑𝐼𝑓𝑑𝑑=πœ†π‘“π‘†π‘“βˆ’(πœ‡+𝑑+𝑧)𝐼𝑓,𝑑𝐴𝑓𝑑𝑑=π‘§πΌπ‘“βˆ’(πœ‡+𝑑)𝐴𝑓.(2.1) The forces of infection for male and female are, respectively, given byπœ†π‘š=πœ‚π‘“π›½π‘šπΌπ‘“π‘π‘“,πœ†π‘“=πœ‚π‘šπ›½π‘“πΌπ‘šπ‘π‘š,(2.2) where π›½π‘š, 𝛽𝑓 and πœ‚π‘š, πœ‚π‘“ are, respectively, the probabilities of acquiring HIV and the average number of male and female sexual partners, respectively. The basic model (2.1) with nonnegative initial conditions is epidemiologically meaningful and mathematically well posed. Thus, system (2.1) is dissipative (i.e., all feasible solutions are uniformly bounded [6, 7]). System (2.1) has a disease-free equilibrium (DFE) given by𝐸0=𝑆0π‘š,𝐼0π‘š,𝐴0π‘š,𝑆0𝑓,𝐼0𝑓,𝐴0𝑓=ξ‚΅Ξ›π‘šπœ‡Ξ›,0,0,π‘“πœ‡ξ‚Ά,0,0.(2.3) Using the next-generation operator method [8], the basic model reproduction number 𝑅0 of model (2.1), defined as the number of secondary infections caused by a typical infected individual introduced into the entire susceptible population during his entire period of infectiousness [9], is given by𝑅0=ξƒŽπ›½π‘“πœ‚π‘š(π›½πœ‡+𝑑+β„Ž)π‘šπœ‚π‘“(πœ‡+𝑑+𝑧).(2.4) The expression of 𝑅0 is a geometric mean of the average number of secondary male infections produced by one female, and the average number of secondary female infections produced by one male. From Theorem 2 of van den Driessche and Watmough [8], the following result holds.

Lemma 2.1. The DFE 𝐸0 of system (2.1) is locally asymptotically stable if 𝑅0<1, and unstable if 𝑅0>1.

If 𝑅0<1, then on average an infected individual produces less than one new infection over its infectious period, and the epidemic cannot grow. That is, a small influx of infected individuals would not generate large outbreaks if 𝑅0<1. Conversely, if 𝑅0>1, then each infected individual produces on average more than one new infection, and the disease can invade the population. However, in order for disease elimination to be independent of the initial sizes of the subpopulations of the model when 𝑅0<1, global stability of 𝐸0 is required.

Lemma 2.2. The 𝐸0 is globally asymptotically stable if 𝑅0<1, and unstable otherwise.

Proof. The proof is based on a comparison theorem [10]. The rate of change of the variables representing the infected components of the system (2.1) can be written as βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽπ‘‘πΌπ‘šπ‘‘π‘‘π‘‘π΄π‘šπ‘‘π‘‘π‘‘πΌπ‘“π‘‘π‘‘π‘‘π΄π‘“βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽπΌπ‘‘π‘‘=(πΉβˆ’π‘‰)π‘šπ΄π‘šπΌπ‘“π΄π‘“βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ βˆ’βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽπœ†π‘šπ‘†π‘š0πœ†π‘“π‘†π‘“0⎞⎟⎟⎟⎟⎟⎟⎠,(2.5) where the matrices 𝐹 and 𝑉 are given, respectively, by βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽπ›½πΉ=00π‘šπœ‚π‘“Ξ›π‘šπœ‡π‘π‘“0𝛽0000π‘“πœ‚π‘šΞ›π‘“πœ‡π‘π‘šβŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽβŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ 0000000,𝑉=πœ‡+𝑑+β„Ž000βˆ’β„Žπœ‡+𝑑0000πœ‡+𝑑+𝑧000βˆ’π‘§πœ‡+𝑑.(2.6) Thus, βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽπ‘‘πΌπ‘šπ‘‘π‘‘π‘‘π΄π‘šπ‘‘π‘‘π‘‘πΌπ‘“π‘‘π‘‘π‘‘π΄π‘“βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽπΌπ‘‘π‘‘β‰€(πΉβˆ’π‘‰)π‘šπ΄π‘šπΌπ‘“π΄π‘“βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ .(2.7) Using the fact that the eigenvalues of the matrix πΉβˆ’π‘‰ all have negative real parts, it follows that the linearized differential inequality above is stable whenever 𝑅0<1. Consequently, (πΌπ‘š,π΄π‘š,𝐼𝑓,𝐴𝑓)β†’(0,0,0,0) as π‘‘β†’βˆž. By a comparison Theorem [10], (πΌπ‘š,π΄π‘š,𝐼𝑓,𝐴𝑓)β†’(0,0,0,0) as π‘‘β†’βˆž. Thus, (π‘†π‘š,πΌπ‘š,π΄π‘š,𝑆𝑓,𝐼𝑓,𝐴𝑓)β†’(Ξ›π‘š/πœ‡,0,0,Λ𝑓/πœ‡,0,0) as π‘‘β†’βˆž for 𝑅0<1, and hence, the DFE is globally asymptotically stable if 𝑅0<1.

The above result indicates that HIV could be eliminated from the community if the threshold quantity 𝑅0 can be brought to (and maintained at) a value less than unity. The endemic equilibrium (EE) of (2.1) is given byπΈβˆ—=ξ‚€π‘†βˆ—π‘š,πΌβˆ—π‘š,π΄βˆ—π‘š,π‘†βˆ—π‘“,πΌβˆ—π‘“,π΄βˆ—π‘“ξ‚,(2.8) whereπ‘†βˆ—π‘š=Ξ›π‘šξ€·πœ†βˆ—π‘šξ€Έ+πœ‡,πΌβˆ—π‘š=Ξ›π‘šπœ†βˆ—π‘šξ€·πœ†βˆ—π‘šξ€Έ+πœ‡(πœ‡+𝑑+β„Ž),π΄βˆ—π‘š=β„ŽΞ›π‘šπœ†βˆ—π‘šξ€·πœ†βˆ—π‘šξ€Έ,𝑆+πœ‡(πœ‡+𝑑+β„Ž)(πœ‡+𝑑)βˆ—π‘“=Ξ›π‘“ξ‚€πœ†βˆ—π‘“ξ‚+πœ‡,πΌβˆ—π‘“=Ξ›π‘“πœ†βˆ—π‘“ξ‚€πœ†βˆ—π‘“ξ‚+πœ‡(πœ‡+𝑑+𝑧),π΄βˆ—π‘“=π‘§Ξ›π‘“πœ†βˆ—π‘“ξ‚€πœ†βˆ—π‘“ξ‚,+πœ‡(πœ‡+𝑑+𝑧)(πœ‡+𝑑)(2.9) where πœ†βˆ—π‘š=πœ‚π‘“π›½π‘š(πΌβˆ—π‘“/π‘βˆ—π‘“) and πœ†βˆ—π‘“=πœ‚π‘šπ›½π‘“(πΌβˆ—π‘š/π‘βˆ—π‘š). Solving for πœ†βˆ—π‘š using the values of πΌβˆ—π‘š, πΌβˆ—π‘“ in (2.9) and the value of πœ†βˆ—π‘“, after some lengthy algebraic manipulations, the endemic equilibrium of HIV/AIDS basic model satisfies the following linear equation:π΄πœ†βˆ—π‘šβˆ’π΅=0,(2.10) where𝐴=πœ‡π‘π‘šπ‘π‘“(πœ‡+𝑑+β„Ž)(πœ‡+𝑑+𝑧)+π‘π‘“πœ‚π‘šπ›½π‘“Ξ›π‘š(πœ‡+𝑑+𝑧),𝐡=πœ‡2π‘π‘šπ‘π‘“ξ‚Έπœ‚(πœ‡+𝑑+β„Ž)(πœ‡+𝑑+𝑧)1βˆ’π‘“π›½π‘šΞ›π‘“πœ‚π‘šπ›½π‘“Ξ›π‘šπœ‡2π‘π‘šπ‘π‘“ξ‚Ή(πœ‡+𝑑+β„Ž)(πœ‡+𝑑+𝑧)=πœ‡2π‘π‘šπ‘π‘“ξ€Ί(πœ‡+𝑑+β„Ž)(πœ‡+𝑑+𝑧)1βˆ’π‘…20ξ€»,(2.11)𝐴>0 while 𝐡<0 provided 𝑅0>1, and consequently, the linear system π΄πœ†βˆ—π‘šβˆ’π΅=0 has a unique positive solution πœ†βˆ—π‘š=𝐡/𝐴, whenever 𝑅0>1. The components of the endemic equilibrium πΈβˆ— are then determined by substituting πœ†βˆ—π‘š=𝐡/𝐴 into (2.9). Noting that 𝑅0<1 implies 𝐡>0; thus, for 𝑅0<1, the force of infection πœ†βˆ—π‘š at steady state is negative (and biologically meaningless). Hence, the model has no endemic equilibrium in this case. Thus, we have established the following result.

Lemma 2.3. The HIV/AIDS model (2.1) has a unique positive EE πΈβˆ— whenever 𝑅0>1 and none otherwise.

The uniqueness of πΈβˆ— and the global stability of the DFE imply that the model does not exhibit the phenomenon of backward or subcritical bifurcation where a locally stable EE coexists with a stable DFE when the reproduction number is less than unity. Because the model parameters are taken from different sources, a deterministic sensitivity analysis is carried out using the approach in [11]. Sensitivity indices (of the reproduction number) which measure initial disease transmission allow us to estimate the relative change in a state variable when a parameter changes. The sensitivity indices of 𝑅0 to the parameters for the HIV/AIDS model are given in Table 1. The negative parameters simply means that an increase in that parameter leads to a decrease in the reproductive number. For instance, a 10% increase of the number of sexual partners will lead to a 5% increase of the value of 𝑅0 (initial disease transmission threshold).

3. Analysis of the Model with Interventions

The basic model is extended to include screening (detected classes, π·π‘š, 𝐷𝑓) and treatment classes (π‘‡π‘š, 𝑇𝑓). It is assumed that the number of contacts made by susceptible individuals under treatment is less than or equal to the number of contacts made with an untreated infective due to behavioral change. This is captured via the parameters π‘Ÿπ‘š and π‘Ÿπ‘“ which are both less than or equal to unity (0β‰€π‘Ÿπ‘š, π‘Ÿπ‘“β‰€1). Also, treatment reduces infectivity, accounted herein by the parameters π‘π‘š, 𝑝𝑓<1. The model compartments and flow are depicted in Figure 2, while the additional variables and parameters for the extended model are described in Tables 2 and 3, respectively.

With the above assumptions and terminology, the model is given by the following system of nonlinear equations: Maleπ‘‘π‘†π‘šπ‘‘π‘‘=Ξ›π‘šβˆ’ξ€·πœ†π‘šξ€Έπ‘†+πœ‡π‘š,π‘‘πΌπ‘šπ‘‘π‘‘=πœ†π‘šπ‘†π‘šβˆ’ξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€ΈπΌπ‘š,π‘‘π·π‘šπ‘‘π‘‘=πœŽπ‘šπΌπ‘šβˆ’ξ€·πœ‡+𝑑+πœπ‘šξ€Έπ·π‘š,π‘‘π‘‡π‘šπ‘‘π‘‘=πœπ‘šπ·π‘šβˆ’ξ€·πœ‡+𝑑+π›Όπ‘šξ€Έπ‘‡π‘š,π‘‘π΄π‘šπ‘‘π‘‘=π›Όπ‘šπ‘‡π‘š+β„ŽπΌπ‘šβˆ’(πœ‡+𝑑)π΄π‘š,Female𝑑𝑆𝑓𝑑𝑑=Ξ›π‘“βˆ’ξ€·πœ†π‘“ξ€Έπ‘†+πœ‡π‘“,𝑑𝐼𝑓𝑑𝑑=πœ†π‘“π‘†π‘“βˆ’ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€ΈπΌπ‘“,𝑑𝐷𝑓𝑑𝑑=πœŽπ‘“πΌπ‘“βˆ’ξ€·πœ‡+𝑑+πœπ‘“ξ€Έπ·π‘“,𝑑𝑇𝑓𝑑𝑑=πœπ‘“π·π‘“βˆ’ξ€·πœ‡+𝑑+𝛼𝑓𝑇𝑓,𝑑𝐴𝑓𝑑𝑑=𝛼𝑓𝑇𝑓+π‘§πΌπ‘“βˆ’(πœ‡+𝑑)𝐴𝑓.(3.1) The force of infections of male and female is, respectively, given by πœ†π‘š=πœ‚π‘“π›½π‘šπΌπ‘“+π‘Ÿπ‘“π·π‘“+𝑝𝑓𝑇𝑓𝑁𝑓,πœ†π‘“=πœ‚π‘šπ›½π‘“πΌπ‘š+π‘Ÿπ‘šπ·π‘š+π‘π‘šπ‘‡π‘šπ‘π‘š.(3.2) The DFE of model system (3.1) denoted by 𝐸00 is given by 𝐸00=ξ‚΅Ξ›π‘šπœ‡Ξ›,0,0,0,0,0,π‘“πœ‡ξ‚Ά,0,0,0,0,0.(3.3) Using the next-generation matrix operator [8], the reproduction number of (3.1) is 𝑅𝑇=ξ”π‘…π‘šπ‘‡π‘…π‘“π‘‡,(3.4) where π‘…π‘šπ‘‡=π›½π‘“πœ‚π‘šπœ‡+𝑑+β„Ž+πœŽπ‘šξƒ©π‘Ÿ1+π‘šπœŽπ‘šπœ‡+𝑑+πœπ‘š+π‘π‘šπœπ‘šπœŽπ‘šξ€·πœ‡+𝑑+πœπ‘šξ€Έξ€·πœ‡+𝑑+π›Όπ‘šξ€Έξƒͺ,𝑅𝑓𝑇=π›½π‘šπœ‚π‘“πœ‡+𝑑+𝑧+πœŽπ‘“ξƒ©π‘Ÿ1+π‘“πœŽπ‘“πœ‡+𝑑+πœπ‘“+π‘π‘“πœπ‘“πœŽπ‘“ξ€·πœ‡+𝑑+πœπ‘“ξ€Έξ€·πœ‡+𝑑+𝛼𝑓ξƒͺ.(3.5) Thus, using Theorem 2 of van den Driessche and Watmough [8], the following result holds.

Lemma 3.1. The DFE 𝐸00 of model system (3.1) is locally asymptotically stable if 𝑅𝑇<1, and unstable if 𝑅𝑇>1.

π‘…π‘šπ‘‡ and 𝑅𝑓𝑇 are the reproduction numbers for males and females, respectively (i.e., π‘…π‘šπ‘‡ represents the number of females infected by a single male during his entire period of infectiousness in a population where treatment is available). If the interventions are dropped, that is, πœŽπ‘š=0, πœŽπ‘“=0, and πœπ‘š=0, πœπ‘“=0, then the effective reproduction number 𝑅𝑇 reduces to the basic reproduction number 𝑅0. The endemic equilibrium of (3.1) denoted by πΈβˆ—βˆ— is πΈβˆ—βˆ—=ξ‚€π‘†π‘šβˆ—βˆ—,πΌπ‘šβˆ—βˆ—,π·π‘šβˆ—βˆ—,π‘‡π‘šβˆ—βˆ—,π΄π‘šβˆ—βˆ—,π‘†π‘“βˆ—βˆ—,πΌπ‘“βˆ—βˆ—,π·π‘“βˆ—βˆ—,π‘‡π‘“βˆ—βˆ—,π΄π‘“βˆ—βˆ—ξ‚,(3.6)   π‘†π‘šβˆ—βˆ—=Ξ›π‘šπœ‡+π΄πΌπ‘“βˆ—βˆ—,π‘†π‘“βˆ—βˆ—=Ξ›π‘“πœ‡+π΅πΌπ‘šβˆ—βˆ—,π·π‘šβˆ—βˆ—=πœŽπ‘šπœ‡+𝑑+πœπ‘šπΌπ‘šβˆ—βˆ—,π‘‡π‘šβˆ—βˆ—=πœŽπ‘šπœπ‘šξ€·πœ‡+𝑑+πœπ‘šξ€Έξ€·πœ‡+𝑑+π›Όπ‘šξ€ΈπΌπ‘šβˆ—βˆ—,π‘‡π‘“βˆ—βˆ—=πœŽπ‘“πœπ‘“ξ€·πœ‡+𝑑+πœπ‘“ξ€Έξ€·πœ‡+𝑑+π›Όπ‘“ξ€ΈπΌπ‘“βˆ—βˆ—,πΌπ‘šβˆ—βˆ—=π΄Ξ›π‘“ξ€Ίπ΄π΅Ξ›π‘šΞ›π‘“βˆ’πœ‡2ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έξ€»ξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έξ€·πœ‡π΄π΅Ξ›π‘š+𝐴2π΅Ξ›π‘šΞ›π‘“ξ€Έ,πΌπ‘“βˆ—βˆ—=π΄π΅Ξ›π‘šΞ›π‘“βˆ’πœ‡2ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έπœ‡ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έπ΄+π΄π΅Ξ›π‘š,π·π‘“βˆ—βˆ—=πœŽπ‘“πœ‡+𝑑+πœπ‘“πΌπ‘“βˆ—βˆ—,π΄π‘šβˆ—βˆ—=π›Όπ‘šπœŽπ‘šπœπ‘šξ€·+β„Žπœ‡+𝑑+πœπ‘šξ€Έξ€·πœ‡+𝑑+π›Όπ‘šξ€Έξ€·(πœ‡+𝑑)πœ‡+𝑑+πœπ‘šξ€Έξ€·πœ‡+𝑑+π›Όπ‘šξ€ΈπΌπ‘šβˆ—βˆ—,π΄π‘“βˆ—βˆ—=π›Όπ‘“πœŽπ‘“πœπ‘“ξ€·+π‘§πœ‡+𝑑+πœπ‘“ξ€Έξ€·πœ‡+𝑑+𝛼𝑓(πœ‡+𝑑)πœ‡+𝑑+πœπ‘“ξ€Έξ€·πœ‡+𝑑+π›Όπ‘“ξ€ΈπΌπ‘“βˆ—βˆ—,(3.7) where 𝛽𝐴=π‘šπœ‚π‘“π‘π‘“ξƒ©π‘Ÿ1+π‘“πœŽπ‘“πœ‡+𝑑+πœπ‘“+π‘π‘“πœπ‘“πœŽπ‘“ξ€·πœ‡+𝑑+πœπ‘“ξ€Έξ€·πœ‡+𝑑+𝛼𝑓ξƒͺ,𝛽𝐡=π‘“πœ‚π‘šπ‘π‘šξƒ©π‘Ÿ1+π‘šπœŽπ‘šπœ‡+𝑑+πœπ‘š+π‘π‘šπœπ‘šπœŽπ‘šξ€·πœ‡+𝑑+πœπ‘šξ€Έξ€·πœ‡+𝑑+π›Όπ‘šξ€Έξƒͺ.(3.8) The endemic equilibrium exists provided πΌπ‘šβˆ—βˆ—>0 and πΌπ‘“βˆ—βˆ—>0. Consider πΌπ‘“βˆ—βˆ— given in (3.7), then πΌπ‘“βˆ—βˆ—=π΄π΅Ξ›π‘šΞ›π‘“βˆ’πœ‡2ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έπœ‡ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έπ΄+π΄π΅Ξ›π‘š=π΄π΅Ξ›π‘šΞ›π‘“βˆ’πœ‡2ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έ=ξƒ¬π΄π΅Ξ›π‘šΞ›π‘“πœ‡2ξ€·πœ‡+𝑑+𝑧+πœŽπ‘“ξ€Έξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έξƒ­.βˆ’1(3.9) Substituting 𝐴 and 𝐡 in (3.9), after some rearrangement, we obtain πΌπ‘“βˆ—βˆ—=𝑅2π‘‡βˆ’1.(3.10) It is therefore evident that πΌπ‘“βˆ—βˆ—>0 provided 𝑅𝑇>1. A similar expression can be derived for πΌπ‘šβˆ—βˆ—. Thus, we have established the following result.

Lemma 3.2. The model system (3.1) has a unique positive EE whenever 𝑅𝑇>1 and none otherwise.

We analytically investigate the impact of case detection on HIV/AIDS dynamics. By partially differentiating π‘…π‘šπ‘‡ with respect to the case detection rate πœŽπ‘š, we obtain πœ•π‘…π‘šπ‘‡πœ•πœŽπ‘š=π›½π‘“πœ‚π‘šξ€·πœ‡+𝑑+β„Ž+πœŽπ‘šξ€Έ2(Ξ”βˆ’1),(3.11) where Ξ”=((πœ‡+𝑑+β„Ž)/(πœ‡+𝑑+β„Ž+πœπ‘š))(π‘Ÿπ‘š+π‘π‘šπœπ‘š/(πœ‡+𝑑+π›Όπ‘š)). For Ξ”<1, (3.11) is negative and consequently, early detection will always have a positive impact on the dynamics of HIV/AIDS. If Ξ”=1, then case detection has no impact (this case will only arise if treatment of those eligible does not follow). From epidemiological and demographical standpoint, the threshold parameter Δ≀1 is to be expected.

The sensitivity indices of 𝑅𝑇 are given in Table 4. Next, we numerically investigate the impact of the number of partners, case detection, and treatment on the disease dynamics. Tanzania started HIV/AIDS care and treatment in October 2004, and the target for the first year was to cover 44,000 patients. About 96 care and treatment providing facilities were selected to initiate the services, which included four referral hospitals; Muhimbili, Kilimanjaro Christian Medical Centre (KCMC), Bugando Medical Centre (BMC), and Mbeya Medical Centre. Some of the parameter values (see Table 3) are provided courtesy of the regional medical officer estimated based on data from the Bugando Medical Centre (BMC) and the Sekou Toure Hospital (STH) both located in Mwanza City in northern Tanzania. Others are taken from the literature, and the remaining ones are assumed within realistic range for the purpose of illustration.

Multiple partnerships increase the risk factor of acquiring HIV. When the number of partners is small over a long time period, the rate of infection is minimal. In this case, the disease tends to die down (Figure 3(a)). The disease will persist when multiple and concurrent partnerships are frequent in the community, in which case infections are on the increase over time (Figure 3(b)).

Since it is assumed that all detected HIV-positive individuals are treated if eligible, the shapes of the time trends of detected (Figure 4(a)) and treated individuals (Figure 4(b)) are similar. The slight difference is due to the rate of developing full blown AIDS from the treated class.

Figure 5 depicts the effect of increasing early detection and treatment on the reproduction numbers 𝑅0 and 𝑅𝑇.

4. Conclusion

The dynamic and determinants of the HIV epidemics are multiple and are shaped by the sexual patterns which are related to social, cultural, and economic factors: for example, promiscuity, low and inconsistent condom use, intergenerational sex, concurrent sexual partners, and various opportunistic infections. Screening is a barometer for achieving success in the fight against the epidemic. It is against this background that a simple deterministic HIV/AIDS model which accounts for case detection and antiretroviral therapy was formulated and analyzed. Conditions for the global stability which rules out any possibility for the model to exhibit the phenomenon of backward bifurcation were provided. Therefore, the classical requirement of the reproduction number to be less than unity might be sufficient for disease elimination. However, HIV/AIDS is inherently a multifaceted disease with poverty-drug use-behavioral/attitudinal change and gender inequality are some of the social factors that need to be addressed. Sensitivity results point to the case detection rate as a driving factor in stemming the tide of the epidemic. Thus, increasing voluntary and or mass screening will always have a positive impact on HIV/AIDS control in reducing the disease burden.

Numerical simulations clearly show that early detection and treatment alone are insufficient to eliminate the disease (Figures 4 and 5). Other control strategies such as condom and microbicides used are to be explored. The study is not exhaustive and can be extended in various ways by incorporating a potential (imperfect) vaccine, withdrawal from sexual activity of a fraction of individuals in the AIDS defining stage, and development of drug resistance and superinfection (with different virus strains).

Acknowledgments

A. S. Mbitila acknowledges with thanks the support in part of the Norad’s programme for Master Studies (in mathematical modelling) at the University of Dar es Salaam and the Arusha Technical College for financial support and study leave.