Dynamics of the optimality control of transmission of infectious disease: a sensitivity analysis

Over the course of history global population has witnessed deterioration of unprecedented scale caused by infectious transmission. The necessity to mitigate the infectious flow requires the launch of a well-directed and inclusive set of efforts. Motivated by the urge for continuous improvement in existing schemes, this article aims at the encapsulation of the dynamics of the spread of infectious diseases. The objectives are served by the launch of the infectious disease model. Moreover, an optimal control strategy is introduced to ensure the incorporation of the most feasible health interventions to reduce the number of infected individuals. The outcomes of the research are facilitated by stratifying the population into five compartments that are susceptible class, acute infected class, chronic infected class, recovered class, and vaccinated class. The optimal control strategy is formulated by incorporating specific control variables namely, awareness about medication, isolation, ventilation, vaccination rates, and quarantine level. The developed model is validated by proving the pivotal delicacies such as positivity, invariant region, reproduction number, stability, and sensitivity analysis. The legitimacy of the proposed model is delineated through the detailed sensitivity analysis along with the documentation of local and global features in a comprehensive manner. The maximum sensitivity index parameters are disease transmission and people moved from acute stages into chronic stages whose value is (0.439, 1) increase in parameter by 10 percent would increase the threshold quantity by (4.39, 1). Under the condition of a stable system, we witnessed an inverse relationship between susceptible class and time. Moreover, to assist the gain of the fundamental aim of this research, we take the control variables as time-dependent and obtain the optimal control strategy to minimize infected populations and to maximize the recovered population, simultaneously. The objectives are attained by the employment of the Pontryagin maximum principle. Furthermore, the efficacy of the usual health interventions such as quarantine, face mask usage, and hand sanitation are also noticed. The effectiveness of the suggested control plan is explained by using numerical evaluation. The advantages of the new strategy are highlighted in the article.


Model structure
Mathematical models of infectious diseases are crucial for comprehending the patterns of disease spread and devising effective disease control strategies.Consequently, when constructing such models, emphasis should be placed on delineating the epidemiological characteristics of the disease and identifying pivotal, modifiable parameters for disease control.Numerous epidemiological models, rooted in diverse disease transmission mechanisms, have been established and documented in the literature, see 34 .These models have played a pivotal role in shaping and crafting control strategies for various diseases.
Let us consider that, for any time t > 0 , the vulnerable populations under the study, say N(t), is divided into five compartments such as Susceptible class S(t) , Acute infected class I 1 (t) , Chronic infected class I 2 (t) , Recov- ered class R(t) , and Vaccinated class V(t) .The humans specified by S(t) are those who are at risk of contracting a specific infectious disease.Thus, the susceptible who are currently infected are moved into two classes namely Acute infected class I 1 (t) and Chronic infected class I 2 (t) means that if people have the time period of disease are less than six month they moved into acute infected class, if people have time period of disease beyond the six month they moved into chronic infected class, and those who are vaccinated at the rate σ move to V(t) .The vaccinated individuals may either move to Recovered class R(t) at the rate φ after receiving an infection as a result of having contact with an infectious person.Moreover, some parameters are also used in the model formulation, such as the u rate of vaccinated people, σ rate of vaccinated individuals become susceptible again, β disease transmission rate,γ reduced the rate of the incubation period, ρ 1 and ρ 2 recovery rate from acute and chronic stages naturally, γ 2 and γ 3 recovery rate of the acute and chronic stage, d untreated death rate, d 1 disease-related death rate, φ a treated person could regain into susceptible since vaccination rate is not perfect and γ 1 people in the acute stage are moved into chronic stage.
The proposed model leverages some key assumptions grounded in biological evidence to illuminate infectious disease dynamics under optimal control strategies.One critical assumption is that the infectious disease system exhibits backward bifurcation, implying complex and non-linear interactions that contribute to sustained transmission.This assumption aligns with empirical observations in infectious disease ecology, where factors such as host heterogeneity and co-infection dynamics often lead to intricate and non-trivial system behavior.Moreover, the incorporation of optimal control strategies into the model reflects the imperative of intervention measures in disease management.Biological evidence supporting this assumption can be drawn from successful instances of interventions like vaccination campaigns and quarantine measures that have demonstrated efficacy in mitigating disease spread.While the model's novelty or modification status is not explicitly stated, its integration of biological evidence and optimal control strategies contributes to advancing our understanding of infectious disease dynamics and intervention planning 35 .
Therefore, based on the aforementioned considerations, the interaction flow pattern between infectious disease states is depicted in Fig. 1.
Thus, the whole human population at any time t is given by The infectious disease model is depicted through a non-linear system of ordinary differential equations while connecting considered compartments in deterministic formation such as.www.nature.com/scientificreports/ The initial conditions are defined as follows, Furthermore, the details of each term of the equations influencing the system given in the model ( 1) are documented in Table 1.

Qualitative analysis of the model
Here, in this section, some significant results of the proposed model, namely positivity, invariant region, equilibrium point, and basic reproduction number are given.

Positivity
Positivity is enforced upon the initial conditions or parameters inherent in the model.It ensures that the values of the initial conditions and the parameters in the proposed model are endured greater than zero or non-negative.Therefore, reorganize the model (1) in terms is explained by

Invariant region
We determined the invariant region where the solution of the proposed model is bounded.Now, consider the total population of the given model as Differentiate (4) respecting the solution of the model ( 1) provides ( 1)  www.nature.com/scientificreports/Therefore, by substituting all the state equations from the model (1) in Eq. ( 5), then one gets Taking the integration on both side of (6) with t → ∞ , one obtain Clearly, is positively invariant, inside which the model is considered to be epidemiological meaningful and mathematically well-posed.

Equilibrium points
Equilibrium points indicate a stable condition in the progression of the disease, where the occurrence of new infections is offset by the combined occurrences of recoveries and deaths.This section outlines the calculation of the disease-free equilibrium point and the endemic equilibrium point of the proposed system.The equilibrium points associated with the model ( 1) by solving dS dt = dI 1 dt = dI 2 dt = dR dt = dV dt = 0 .Accordingly, the disease-free equilibrium point F 0 , is attained by letting for I 1 = 0, I 2 = 0, R = 0, V = 0 such as, Now, the endemic equilibrium point is denoted as

The basic reproduction number
The basic reproduction number R 0 is a prevailing metric for evaluating the contagiousness or transmissibility of a communicable disease.It signifies the average number of individuals that are infected by one infectious individual in a susceptible population.Biologically, it reflects the contagiousness and transmissibility of a pathogen.If R 0 > 1 , it indicates that the disease has the potential to sustain transmission within the population, leading to an outbreak.Conversely, if R 0 < 1 , the infection is unlikely to establish a self-sustaining chain of transmission, and it may eventually die out in the population.Moreover, the reproduction number is estimated by using the method of Darraish and Wathmough 37 , and with the principle of the next-generation matrix, we have Therefore, the reproduction number R 0 for our proposed model ( 1) is the spectral radius of the next generation matrix FV −1 and calculated as, Therefore, from (7), the reproduction number R 0 is calculated as, where The disease will spread in the population if the basic reproduction number R 0 > 1. (

Sensitivity analysis
The degree of influence of input parameters over the dynamics of the infectious disease model is an asset through the launch of sensitivity analysis.The sensitivity indices of the infectious disease model given in (1) are deduced by using the approach of Chitnis et al. 38 .We calculate the normalized forward sensitivity indices of a parameter x of R 0 .Let The estimated sensitivity projected in reproduction number with respect to various parameters is given as Table 2 comprehends the directional projection of the parameters involved in the considered model explaining the viral infectious disease transmission flow.The existence of a positive relationship is voted through the plus sign whereas the negative sign avows the showing of the negative relationship between parameters and the transmission rate.One may be noted that an increase in the reproduction number remains associated with an increase in parameters such as disease transmission rate, reduced incubation period, rate of transfer from the acute stage to the chronic stage, and untreated death rate.Whereas, the reproduction number is estimated to be negatively associated with parameters such as disease-related death rate, and recovery rates belonging to the acute stage and chronic stage.
Figure 2 shows the dynamical behavior of the sensitivity analysis of the basic reproduction number R 0 , in the context of the disease transmission.Moreover, this analysis helps us to understand the significance of various factors that play a vital role in the transmission of infectious diseases.We plotted the reproduction number R 0 as a function of γ 3 , d 1 , γ 1 , γ 2 , β and d. Figure 2a shows the effect of effectual disease-related death rate and recovery rate of chronic infected patients on the reproduction number, we can see that these two parameters may increase to R 0 greater than unity if the value of γ 3 remain up to greater than 1 and the recovery rate also remains greater than 1.The combined effect of transmission and recovery rate on R 0 is shown in Fig. 2b.As one parameter β is directly proportional and γ 1 is inversely related to R 0 , thus we can see that the basic reproduction number value may be reduced below unity even if the transmission rate β approaches 1 and the recovery rate of acute and chronic infected increases.Figure 2c shows the effect of transmission rate β and d.We can see that these two parameters have a considerable impact in reducing or increasing the value of the basic reproduction number R 0 .Furthermore, Fig. 2d shows the effect of recovery rate gamma 2 and disease-related death rate d, and also Fig. 2e shows the inverse relationship between R 0 and parameters.Similarly, Fig. 2f shows the impact of reducing or increasing the value of reproduction number R 0 .

Stability analysis
In this section, we analyze the local and global stabilities of the proposed model at equilibrium points.Local stability is often determined by examining the eigenvalues of the Jacobian matrix evaluated at the equilibrium point.If all eigenvalues have negative real parts, the equilibrium is locally stable.Global stability considers the system's behavior across its entire range, often requiring more advanced mathematical techniques such as Lyapunov analysis.

Local stability analysis
The aim of conducting local stability analysis is to determine whether a minor disturbance in the disease system would lead to the persistence or elimination of the disease.In the realm of epidemiological models, comprehending the concept of local stability at equilibrium points holds significance.It provides insights into the dynamics of infectious diseases and their transmission within a population.The endemic equilibrium point and disease-free equilibrium point, which are locally asymptotically stable, are illustrated in the following theorem.

Theorem 4.1
The disease-free equilibrium point F 0 is locally asymptotically stable when R 0 < 1, however, when Table 2. Sensitivity indices and parameters of the reproduction number R 0 .

Parameters Sensitivity indices Parameters Sensitivity indices
Vol.:(0123456789) www.nature.com/scientificreports/Proof The Jacobin matrix of disease-free equilibrium point F 0 for the proposed model is given as The evaluation of the jacobian at drinking-free equilibrium provides . www.nature.com/scientificreports/With the aid of the Jacobin matrix of the disease-free equilibrium point, the resultant eigenvalues are , and 5 are real and negative.Further, 3 is negative and real if R 0 < 1 .Using the Routh-Hurwitz criterion 39 , it is verified that, each eigenvalue of the polynomial equation contains a non-positive real parts when R 0 < 1 .As a result, the disease-free equilibrium point F 0 is locally asymptotically stable.

Theorem 4.2
The endemic equilibrium point F * will be locally asymptotically stable if the value of R 0 > 1, however, when R 0 < 1 , it is unstable.
Proof The determination of the Linearization for model (1) at the endemic equilibrium point F * is carried out by By using row matrix transformation the following matrix is given by: where, let µ = β(I It is clear that all of eigenvalues values of J(F * ) have negative and real if R 0 > 1 .

Global stability analysis
Global stability in a disease model examines whether, over the entire range of conditions, the system converges to a stable state in terms of disease behavior.It is crucial for understanding the persistent, long-term dynamics and overall prevalence of drinking in a population.The endemic equilibrium point and disease-free equilibrium point, which are locally asymptotically stable, are illustrated in the subsequent theorem.Theorem 4. 3 The disease free equilibrium point F 0 is globally asymptotically stable for R 0 < 1 , otherwise, unstable for R 0 > 1.
Proof The following Lyapunov function has been properly defined and satisfies the condition of being positive definite and its derivative is negative definite.By focusing on infected and susceptible individuals, the Lyapunov function effectively quantifies the balance between disease transmission and population susceptibility, providing insights into the long-term behavior and stability of the epidemic model.To show the result, Firstly, we state the Lyapunov function is specified by, Where f 1 , f 2 and f 3 are positive constant and derivative of L is given by Now, by arranging the above terms yield, (10) Vol.:(0123456789) where F 0 = (1−u)α βγ +d .Further, for S > S 0 and R 0 < 1, then L ′ (t) < 0 .Also if S = S 0 then L ′ (t) = 0 .By Lasla invari- ence principle 40,41 , I 1 = I 2 = 0.The Lyapunov function conditions are all positive definite, the function is positive definite and its derivative is negative definite.Therefore the disease-free equilibrium F 0 is globally asymptotically stable.
Theorem 4. 4 The endemic equilibrium point F * is globally asymptotically stable, if the value of R 0 > 1 , however, when R 0 < 1 , it is unstable.
Proof In order to prove the global stability of the proposed model (1) at the endemic equilibrium point F * , the Castilo Chevez method is used from 42,43 .Further, we consider (1) as: Taking the jacobian and the additive compound matrix of order 2 for the Eq. ( 13), get the following matrix,  .www.nature.com/scientificreports/

Taking integral of l(A), one get
Hence the model ( 1) is globally asymptotically stable.

Numerical results of stability analysis
The Runge-Kutta method is widely employed in epidemiological models for numerical simulation due to its robustness and accuracy in solving differential equations that govern disease dynamics.Epidemiological models often involve a system of coupled differential equations that describe the interactions between various population compartments, such as susceptible, infected, and recovered individuals.The Runge-Kutta fourth-order method is renowned for its high accuracy in approximating solutions of differential equations, providing reliable and precise numerical results.Additionally, RK-4 exhibits strong convergence properties, ensuring stability and efficiency in solving a wide range of mathematical problems.Its versatility makes it suitable for capturing the complex and dynamic nature of infectious disease spread, allowing researchers to explore different scenarios and intervention strategies.The RK method's efficiency and accuracy contribute to the reliability of numerical simulations, aiding in the understanding and prediction of epidemiological dynamics essential for public health decision-making.
In this section, we used the fourth-order Runge-Kutta method to solve the deterministic mathematical model and verify the analysis 45 .Additionally, the use of numerical simulations can provide valuable insights into the dynamic behavior of the model.However, it's important to note that the selection of parameter approximations should be justified and explained in detail.Additionally, it would be useful to discuss any potential limitations or assumptions made in using the Runge-Kutta method and how it affects the accuracy of the results.Variables and parametric values are given in Table 3, which are used for simulations.Furthermore, the time interval is taken from 0 -200 units with initial problems for Susceptible class S(t) , Acute infected class I 1 (t) , Chronic infected class I 2 (t) , Vaccinated class V(t) , and Recovered class R(t) .The applications of the Runge-Kutta method of fourth order on the proposed model lead to the following system:
Step 3: . By running the above findings with the aid of the Matlab software, the graphs presented in the Figs. 3, 4, 5, 6, 7, which shows the dynamic behavior of Susceptible class S(t) , Acute infected class I 1 (t) , Chronic infected class I 2 (t) , Vaccinated class V(t) and Recovered class R(t).
Figures 3, 4, 5, 6, 7 present the dynamics of the stability behavior of all considered classes of the population while taking into account all the employed initial conditions.The magnitude of the susceptible class is varied for values such as 100, 95, 80, and 85 in Fig. 3.The susceptible population increase in starting time after-some time the population decreases rapidly with time.Next, Fig. 4 depicts the stability attainment for varying sizes of the acute infected class such as 20, 18, 16, and 14.Time to gain stability increases with the decrease in the population of the acute infected class.Figure 5. present the relationship between stability and varying extent of ( 18) Vol   the chronic infected class for values such as 15, 12, 9, and 6.In starting the population decreases rapidly after 65 days population is constant with time.In Fig. 6 The recovered class is varied for values as, 4, 5, 3, and 2. The recovered population increase starting time after-some time the population decreases rapidly with time.Lastly, Fig. 7 shows the relationship between stability and vaccinated population for values such as 10, 8, 6, and 4. Vaccinated class increase in starting time after a short time the population decreases rapidly with time.

Optimal control strategy
This section is dedicated to exploring the applicability of the infectious disease model with respect to the optimal control strategy.The employment of optimal control is initiated to intervene in the transmission of infectious disease in a particular community.Control strategies are formulated by employing sensitivity analysis and examining the dynamics of the proposed model.The maximum sensitivity index parameter is (γ 1 , β) whose value is between the interval (0.439, 1) increase in this parameter by 10 percent would increase the threshold quantity by (4.39, 1).Therefore, by using the control variable u 1 (t) , u 2 (t) , and u 3 (t) , we must reduce these values to stop the sickness from spreading.u 1 (t) and u 2 (t) represent (awareness about medication, isolation, and ventilation).We employ the control of another variable u 3 (t) signify (the vaccination rates and quarantine level).Here, our aim is to reduce the infectious population by increasing the number of recovered people and decreasing the number of infected people and vaccinated people by incorporating the time-dependent control parameters, u 1 (t), u 2 (t), and u 3 (t) which are described as • u 1 (t) is the control variable that characterizes the awareness about medication.
• u 2 (t) is the control variable that characterizes the isolation and ventilation.
• u 3 (t) is the control variable that characterizes the Vaccination rates and quarantine level.
The differential equation of the model ( 1), with the inclusion of the control variables, is, now given as along with the initial conditions are The objective of this work is to show that time-dependent parameters minimize the cost of implementing those control measuring techniques 46 .In optimal control of infectious disease models, utilizing three classes of objective functions allows for a comprehensive assessment of intervention strategies.By considering objectives related to minimizing infections, controlling transmission rates, and maximizing the number of recovered people.Therefore, the objective function is specified as where, z 1 , z 2 , z 3 , z 4 , z 5 , and z 6 , represents weight constant associated with each class.The weight constant z 1 , z 2 , z 3 , represents the cost of acute infected people, chronic infected people, and vaccinated people, while z 4 , z 5 , z 6 , represents the control variables.Further, 1 2 z 4 u 2 1 , 1 2 z 5 u 2 2 , 1 2 z 6 u 2 3 describes the awareness about medication, isolation and ventilation, vaccination rates and quarantine level.Our objective is to find the optimal control variables u * 1 , u * 2 , u * 3 such that dependent on control system (19), where we define the control set U as

Existence of the optimal control system
The existence of optimal control is persuaded by considering initial time t = 0 for the fulfillment of each con- dition of the control system.To achieve the goal, the bounded Lebesgue measurable control 47 by considering the initial conditions and ensuring an upper-bound solution for the system.The matter of optimal control is addressed by delving into the Lagrangian and Hamiltonian, whereas the following equation illustrates the optimal control problem within the Lagrangian framework as www.nature.com/scientificreports/ The minimal possible value of optimal control is achieved by providing the Hamiltonian function H in the subsequent form Henceforth in this way the presence of optimal control, we consider the subsequent resultant.First we have to show that (u * 1 , u * 2 , u * 3 ) these control variables are exist.
Proof In order to prove the existence of an optimal control, thru the result in 48 .Subsequently, the state variables and the control variables have all positive values.It is also detected that the control variable set U is convex and closed by the declaration.Furthermore, the control system is bounded which states the compactness of the considered problem.One may notice the convex nature of integrand involves in concerning control set U which guarantees the exist- ence of the optimal control (u * 1 , u * 2 , u * 3 ) .

Optimality condition
To define an optimal approach to the problem ( 19) and (20).First, we investigate the optimal control issues of Eqs. ( 19) and ( 20) using Lagrangian and Hamiltonian functions.The Lagrangian shows the optimal control problem is presented by the subsequent equation The associated Hamiltonian (H) is defined via the following notion F = (F 1 , F 2 , F 3 , F 4 , F 5 ) and = ( 1 , 2 , 3 , 4 , 5 ) .For the minimal value of the Lagrangian, we govern the Hamiltonian (H) for the optimal control problem by way of where, Further, We revive the use of the Pontryagin maximum principle 49,50 for finding the optimal solution to the proposed model (19).Further, the existence of the non-trivial vector functions such as = ( 1 , 2 , 3 , 4 , 5 ) is noticeable, especially when (x * , u * ) is considered as the optimal procedure for the launched control problem in regarding The Maximality condition describe the transversally condition as (t f ) = 0. Theorem 5.2 Given an optimal state (u * 1 , u * 2 , u * 3 ) variables and solution S * , I * 1 , I * 2 , R * , V * of the equivalent state systems (19) and (20), there exist adjoint variables (t) satisfying x = (S, I 1 , I 2 , R, V), = ( 1 , 2 , 3 , 4 , 5 ), F(x, u) = F 1 (x, u), F 2 (x, u), F 3 (x, u), F 4 (x, u), F 5 (x, u).Proof The Pontryagin Maximum Principle is directly applied to produce the adjoint system (25), while (t) = 0 , directly results in the transversal condition.Regarding the class of optimal functions u * 1 , u * 2 and u * 3 , we used ∂H ∂u .In the section that follows, we numerically solve the optimality system.The control system defines the optimality system (19) with the adjoint system (25), boundary circumstances, as well as the best control strategies.Our control strategies, which are illustrated in Figures 8 and 10, to increase the susceptible and recovered population and reduce the infected and vaccinated population were supported by the simulation that was run.

Numerical simulation and discussion
This section is dedicated to delineating the legitimacy of the proposed control scheme through the launch of numerical evaluation.The strategic environment is enriched by involving a diverse and wide range of parametric settings.The foundations of the numerical proceedings are instigated by considering the compartmental initial conditions are taken from 34 for the devised model such as S(0) = 100, I 1 (0) = 20, I 2 (0) = 15, V(0) = 10 and R(0) = 5 .Moreover, the objective of consistency is maintained by taking into account the literature-oriented values of parameters associated with different compartments given in Table 3.We determine the ideal control system (19) to see the impact of medication, ventilation, isolation, and quarantine level.The resolve of the issue of optimization is persuaded through the employment of the Runge-Kutta fourth-order technique along with the application of the transversality criterion within the time interval of [0, 50].Furthermore, the values of weight constant are assumed to be z * 1 = 0.90030, z * 2 = 0.43454, z * 3 = 0.3550, z * 4 = 0.5560, z * 5 = 0.67676 and z * 6 = 0.999 .Graphical representation and table of parameters used in Optimal control analysis are given below: Figures 8, 9, 10 depict the dynamics of comparative effectiveness of the proposed model, with respect to both scenarios that is with control and without control.This description is foreseen for all five compartments with respect to time in days.The gain of the optimal control strategy is noticeable for all compartments of the population.Figures 8 shows the dynamic behavior of Susceptible class S(t) , Acute infected class I 1 (t), Chronic infected class I 2 (t) , and Recovered class R(t) and the Vaccinated class V(t) with and without control interventions.Further, Fig. 8a represents the behavior of number of susceptible population decrease after using the optimal control strategy.Figure 8b Utilizing these control measures will reduce the number of acute infected populations after applying the control.Similarly, Fig. 8c the number of population of chronic infected populations decreased after applying the control.Figure 8d the recovered population rapidly increases and achieves convergence after applying optimal control.Figure 8e the number of the vaccinated population rapidly decreases after applying the control.The recovered population increase, reducing the infected and vaccinated population was evidently supported by the simulation that was run. Figure 10b,c,d are dedicated to highlighting the robustness of control variables with respect to time in days.

Conclusion
This article proposes a novel mathematical model to encapsulate the viral transmission of infectious diseases.The optimal control strategy is launched while considering five related epidemiological compartments.The population under study is divided into five strata such as; Susceptible class S(t) , Acute Infected class I 1 (t) , Chronic Infected class I 2 (t) , Recovered class R(t) and Vaccinated class V(t) .The target of generality is maintained by employing a wide range of parametric settings involving, disease transmission rate from the susceptible class to the acute infected class, the untreated death rate of the susceptible class, the disease-related death rate of the acute infected class, a recovery rate of the chronic infected class.In the procession, rigorous persuasion into the local stability of the system, global stability of the system along with sensitivity analysis, positivity, and invariant region is considered.The implementations of the outcomes of the research are further supported and widened by the inclusion of more popular health interventions such as; quarantine, vaccination, medication, and mandatory use of masks.It is delineated that in the case of stability of the system, the susceptible class increase in starting point after some time population decrease rapidly with respect to time.Moreover, under the norm of stability, the interaction of underlying compartments indicates a higher likelihood of a decrease in the number of Acuteinfected and Chronic-infected individuals.The exploration into the sensitivity parameter indicated an increase of 4.39 units in threshold quantity remained associated with a 10 percent increase in the sensitivity parameter.www.nature.com/scientificreports/approach involved vaccination rates u 3 as time-dependent controls.Graphical results clearly demonstrate the effectiveness of each strategy in reducing both acute and chronic infections, suggesting their potential implementation to mitigate the spread of infectious diseases in the population.The graphical results also revealed that the case where we considered all three control variables, simultaneously, is more effective in reducing the spread of infectious flow.The optimal control strategy advocated the efficacy of the employment of popular health interventions such as; quarantine, face masks, and hand washing, to restrict the transmission chain of the infectious disease.Additionally, our findings concluded that simulations with time-dependent controls are more cost-effective compared to those with time-independent controls.
In the future, it will be interesting to extend the above-developed model to explore the impact of the fractional order approach while maintaining the optimal control in place.It is anticipated that the ability of the fractional order scheme to use available information will be beneficial to further uncover the transmission of diseases in the population.Moreover, rigorous focus will remain on the exploration of the applicability of the proposed scheme on data sets involving varying geographical locations, socio-economic stratifications prevalent in the society along the degree of access to health facilities.Furthermore, it will be engaging to explore the efficacy of the devised model when dealing with multiple infectious flows at the same time.By figuring out how these strategies work together or against each other, we hope to create better plans to control these diseases around the world.

Figure 1 .
Figure 1.Schematic diagram illustrates the transmission dynamics of infectious disease behavior.

Figure 2 .
Figure 2. The graphical result shows the dynamics of various compartmental parameters and their effects on the reproduction number R 0 .

Figure 3 .
Figure 3.The effects of variations in the Susceptible class S(t) are described by simulation results.

Figure 4 .
Figure 4.The effects of variations in the Acute Infected class I 1 (t) are described by simulation results.

Figure 5 .
Figure 5.The effects of variations in the Chronic Infected class I 2 (t) are described by simulation results.

Figure 6 .
Figure 6.The effects of variations in the Recovered class R(t) are described by simulation results.

Figure 7 .
Figure 7.The effects of variations in the Vaccinated class V(t) are described by simulation results.

Table 1 .
Description of each term of the equations of the model (1).

Table 3 .
Parameters and its values.