Optimal control of HIV transmission model with pre-ART counselling and treatment
Main Article Content
Abstract
HIV attacks CD4 cells of the immune system, leading to progressive immune deficiency. Antiretroviral therapy (ART) involves the use of HIV drugs to treat HIV infection and is administered daily to slow disease progression. This paper aims to develop and analyze a mathematical model of HIV transmission that incorporates pre-antiretroviral therapy counselling and HIV treatment to reduce the number of HIV-infected individuals with high-risk behaviours for HIV transmission. A nonlinear dynamical system is constructed, and model parameters are estimated from Indonesia’s annual HIV case data using a genetic algorithm method. The model exhibits two equilibrium points: the disease-free equilibrium and the endemic equilibrium. Stability analysis shows that disease-free equilibrium is globally asymptotically stable when the basic reproduction number is less than one. Optimal control theory is applied to a system that consists of two time-dependent controls, pre-antiretroviral therapy counselling and HIV treatment. Healthcare professionals provide pre-antiretroviral therapy counselling to help people with HIV understand the disease and the benefits of antiretroviral therapy. Pontryagin's maximum principle is employed to derive optimal control conditions. The optimal control problem is numerically solved using the forward–backward sweep method with a fourth-order Runge–Kutta scheme. Three potential strategies were developed and investigated in our simulation. Implementing the two combined controls could significantly reduce the number of HIV-infected individuals and improve overall disease control in the population.
Article Details

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
References
[1] World Health Organization (WHO). HIV and AIDS; 1992. https://www.who.int/news-room/fact-sheets/detail/hiv- aids.
[2] Kemenkes RI. Infodatin HIV. Jakarta: Ministry of Health Republic of Indonesia; 2020.
[3] Division of HIV Prevention. HIV Transmission Topics; 2020. https://www.cdc.gov/hiv/basics/transmission.html.
[4] Burns DN, DeGruttola V, Pilcher CD, Kretzschmar M, Gordon CM, Flanagan EH, et al. Toward an endgame: Finding and engaging people unaware of their HIV-1 infection in treatment and prevention. AIDS Research and Human Retroviruses. 2014 mar;30(3):217-24. doi:10.1089/aid.2013.0274.
[5] Low A, Gavriilidis G, Larke N, B-Lajoie MR, Drouin O, Stover J, et al. Incidence of Opportunistic Infections and the Impact of Antiretroviral Therapy among HIV-Infected Adults in Low- and Middle-Income Countries: A Systematic Review and Meta-analysis. Clinical Infectious Diseases. 2016 jun;62(12):1595-603. doi:10.1093/cid/ciw125.
[6] Nkambule BS, Sambo G, Aydin HZ, Yildiz NG, Aydin K, Yildiz H, et al. Factors associated with HIV-positive status awareness among adults with long term HIV infection in four countries in the East and Southern Africa region: A multilevel approach. PLOS Global Public Health. 2023 dec;3(12):e0002692. doi:10.1371/journal.pgph.0002692.
[7] Brauer F, Castillo-Chavez C, Feng Z. Mathematical Models in Epidemiology. vol. 69 of Texts in Applied Mathematics. New York, NY: Springer New York; 2019. doi:10.1007/978-1-4939-9828-9.
[8] Abidemi A, Fatmawati, Alfiniyah C, Windarto, Nyabadza F, Aziz MHN. Insights into HIV/AIDS transmission dynamics and control in Indonesia — A mathematical modelling study. Partial Differential Equations in Applied Mathematics. 2025 jun;14:101185. doi:10.1016/j.padiff.2025.101185.
[9] Olaosebikan ML, Kolawole MK, Bashiru KA. Transmission Dynamics of Tuberculosis Model with Control Strategies. Jambura Journal of Biomathematics (JJBM). 2023 dec;4(2):110-8. doi:10.37905/jjbm.v4i2.21043.
[10] Sangotola AO, Adigun AJ, Nuga OA, Adeyemo S, Kataboh PK, Akinde OT, et al. An Isolation Model for Tuberculo- sis Dynamics with Optimal Control Application. Communication in Biomathematical Sciences. 2025 jul;8(1):55-65. doi:10.5614/cbms.2025.8.1.4.
[11] Rois MA, Fatmawati, Alfiniyah C, Martini S, Aldila D, Nyabadza F. Modeling and optimal control of COVID-19 with comorbidity and three-dose vaccination in Indonesia. Journal of Biosafety and Biosecurity. 2024 sep;6(3):181-95. doi:10.1016/j.jobb.2024.06.004.
[12] Akanni JO, Abidemi A, Fatmawati F, Chukwu CW. A Non-linear Fractional Model for Analyzing the Impact of Vaccina- tion on the Dynamics of COVID-19 in Indonesia. Jambura Journal of Biomathematics (JJBM). 2025 jun;6(2):109-28. doi:10.37905/jjbm.v6i2.30383.
[13] Huo HF, Feng LX. Global stability for an HIV/AIDS epidemic model with different latent stages and treatment. Applied Mathematical Modelling. 2013 feb;37(3):1480-9. doi:10.1016/j.apm.2012.04.013.
[14] Maimunah, Aldila D. Mathematical model for HIV spreads control program with ART treatment. Journal of Physics: Conference Series. 2018 mar;974(1):012035. doi:10.1088/1742-6596/974/1/012035.
[15] Tabassum MF, Saeed M, Akgül A, Farman M, Chaudhry NA. Treatment of HIV/AIDS epidemic model with verti- cal transmission by using evolutionary Padé-approximation. Chaos, Solitons and Fractals. 2020 may;134:109686. doi:10.1016/j.chaos.2020.109686.
[16] Ibrahim IA, Daniel EE, Danhausa AA, Adamu MU, Shawalu CJ, Yusuf A. Mathematical Modelling of Dynamics of HIV Transmission Depicting the Importance of Counseling and Treatment. Journal of Applied Sciences and Environmental Management. 2021;25(6):893-903. doi:10.4314/jasem.v25i6.1.
[17] Chandra TD, Permata GI. Modeling Hiv/Aids Using Shat Model. Barekeng. 2023 jun;17(2):745-56. doi:10.30598/barekengvol17iss2pp0745-0756.
[18] Fatmawati, Tasman H. An Optimal Treatment Control of TB-HIV Coinfection. International Journal of Mathematics and Mathematical Sciences. 2016;2016:1-11. doi:10.1155/2016/8261208.
[19] Marsudi, Trisilowati, Suryanto A, Darti I. Optimal Control of an HIV Model with Changing Behavior through an Education Campaign, Screening and Treatment. IOP Conference Series: Materials Science and Engineering. 2019 jun;546(5):052043. doi:10.1088/1757-899X/546/5/052043.
[20] Norasia Y, Zulaikha Z, Tafrikan M, Ghani M, Mukama DS. Optimal Control of HIV-1 Spread in Combination with Nutritional Status and ARV-Treatment. International Journal of Computing Science and Applied Mathematics. 2022 sep;8(2):66. doi:10.12962/j24775401.v8i2.13764.
[21] Gurmu ED, Bole BK, Koya PR. Mathematical modelling of HIV/AIDS transmission dynamics with optimal control strat- egy. International Journal of Mathematcs and Computer Research. 2021;9(04):2237-54. doi:10.47191/ijmcr/v9i4.04.
[22] Ministry of Health Republic of Indonesia. Laporan Eksekutif Perkembangan HIV AIDS dan PIMS Tahun 2022. Jakarta; 2023. https://hivaids-pimsindonesia.or.id/download.
[23] Sa’adah A, Sasmito A, Pasaribu AA. Comparison of Genetic Algorithm (GA) and Particle Swarm Optimization (PSO) for Estimating the Susceptible-Exposed-Infected-Recovered (SEIR) Model Parameter Values. Journal of Information Systems Engineering and Business Intelligence. 2024 jun;10(2):290-301. doi:10.20473/jisebi.10.2.290-301.
[24] Central bureau of statistics Indonesia. Population by Province in Indonesia; 2021. https://bekasikab.bps.go.id/id/statistics-table/1/MjQ4MyMx/jumlah-penduduk-menurut-provinsi-di-indonesia-ribu-20122016.html.
[25] Central bureau of statistics Indonesia. The Lifespan of Indonesia; 2021. https://www.bps.go.id/id/statistics- table/2/NTAxIzI=/angka-harapan-hidup-ahh-menurut-provinsi-dan-jenis-kelamin.html.
[26] Brauer F, Castillo-Chavez C. Mathematical Models in Population Biology and Epidemiology. vol. 40 of Texts in Applied Mathematics. New York, NY: Springer New York; 2012. doi:10.1007/978-1-4614-1686-9.
[27] La Salle JP. The Stability of Dynamical Systems. Society for Industrial and Applied Mathematics; 1976. doi:10.1137/1.9781611970432.
[28] Subbaram Naidu D. Optimal control systems. New York: CRC Press; 2002. doi:10.1515/9783110789737-005.