Modeling the effect of vaccinations, hospital beds, and treatments on the dynamics of infectious disease with delayed optimal control and sensitivity analysis
Arjun Kumar,Uma S. Dubey,Balram Dubey
DOI: https://doi.org/10.1140/epjp/s13360-024-05508-1
2024-08-21
The European Physical Journal Plus
Abstract:Immunization plays a vital role in eradicating infectious diseases, typically requiring multiple doses at specific time intervals. This study focuses on developing and analyzing an infectious disease model governed by a six-dimensional system of ordinary differential equations, considering the impact of first and second vaccination doses along with hospital beds and treatment. The model's qualitative behavior is analyzed, including conditions for positive solutions, the invariant region of the solution, equilibrium points, and their stability. When the basic reproduction number ( ) is less than one ( ), the disease will be eradicated; conversely, an epidemic occurs when . Moreover, the transcritical bifurcation of the system is examined using the center manifold theory. Interestingly, backward bifurcation is discovered, and it indicates that the disease is not entirely eradicated even when . We have investigated different bifurcations like saddle-node, transcritical, and Hopf bifurcations of codimension 1, as well as Generalized-Hopf (GH), Cusp (CP), and Bogdanov–Takens (BT) bifurcations of codimension 2. Additionally, a delayed epidemiological model is explored, assuming a lag in vaccination among the susceptible population. A Hopf-bifurcation is observed near the endemic equilibrium point, linked to critical parameter values during the latent period. Moreover, the model is calibrated using the least-squares technique, incorporating coronavirus-infected case data and vaccination information from India and Italy's mass vaccination program between March 1, 2021, and May 30, 2021. Global sensitivity analysis, utilizing the Partial Rank Correlation Coefficient (PRCC), identifies crucial parameters affecting threshold quantities after fitting the model. The study highlights the significance of critical parameters such as the effective transmission rate, rates of first and second-dose vaccinations, and recovery rate due to double-dose vaccination. Further, delayed optimal control measures are determined using Pontryagin's maximal principle to mitigate infection, prevention, and treatment burdens. Numerical simulations are conducted to understand the effect of these delayed control measures on disease progression and demonstrate the insights obtained through analytical investigations. The study indicates that implementing all control strategies effectively reduces the disease burden among the population. Accurate estimation of vaccine efficacy is crucial for disease prevention, underlining the importance of well-planned vaccination strategies. Moreover, the numerical simulations validate all the theoretical findings, emphasizing the validity of this model in a real-world situation. Relying solely on vaccination might not swiftly or completely control the disease. Complementary pharmaceutical and non-pharmaceutical measures are necessary to combat the infection effectively. Further limitations on medical resources could lead to a backward bifurcation. Simulation results suggest that delaying the implementation of control measures could exacerbate epidemic situations.
physics, multidisciplinary