Although there are many results that can be used to treat and prevent Coronavirus Disease 2019 (COVID-19) and Human Immunodeficiency Virus (HIV), these diseases continue to be public health concerns and cause socioeconomic consequences. Following compromised immunity, COVID-19 is considered to be a challenge for people with HIV. People with advanced HIV are considered a vulnerable population at high risk in several case studies that discuss COVID-19 and HIV co-infection. As there is no cure for HIV and there is a chance of contracting COVID-19 again, co-infection continues to pose a problem. The purpose of this study is to investigate the impact of intervention strategies and identify the role of different parameters in risking people living with HIV to death when they get infected with COVID-19. This is achieved through the development and rigorous analysis of a mathematical model that considers a population at risk of death due to COVID-19 and HIV. The model formulation provides a detailed explanation of the transmission dynamics of COVID-19 and HIV co-infection. The solution’s invariant region, positivity, and boundedness were established. The reproduction numbers of the sub-models and the co-infection model were determined. The existence and stability of equilibria, including backward bifurcation for the COVID-19 sub-model, were examined. The epidemiological significance of backward bifurcation is that the condition $${\mathscr {R}}_0$$
R
0
less than 1 for eliminating COVID-19, though necessary, is no longer sufficient. Parametric estimation and curve fitting were performed based on data from Ethiopia. Numerical simulations were employed to support and clarify the analytical findings and to show some parameter effects on COVID-19 and HIV co-infection. Accordingly, the simulations indicated that parameters $$\gamma _c$$
γ
c
, $$\gamma _h$$
γ
h
, $$\epsilon$$
ϵ
, and $$\kappa$$
κ
, related to HIV patients’ exposure to other diseases and the increase in infectiousness, have a positive role in increasing the number of co-infections. On the other hand, an increase in COVID-19 vaccination ($$\xi$$
ξ
) shows the suppression of co-infection cases. In addition, treating co-infected individuals for COVID-19, increasing treatment rates $$\alpha$$
α
and $$\varphi$$
φ
, reduces the death risk of HIV-infected individuals due to the co-infection burden. It was implied that improving vaccine delivery programs and other medical interventions have important contributions to lowering the risk of COVID-19 infection-related fatalities in HIV patients.