Mathematical model to investigate transmission dynamics of COVID-19 with vaccinated class

: The susceptible, exposed, infected, quarantined and vaccinated (SEIQV) population is accounted for in a mathematical model of COVID-19. This model covers the therapy for diseased people as well as therapeutic measures like immunization for susceptible people to enable understanding of the dynamics of the disease’s propagation. Each of the equilibrium points, i.e., disease-free and endemic, has been proven to be globally asymptotically stable under the assumption that R 0 is smaller or larger than unity, respectively. Although vaccination coverage is high, the basic reproduction number depends on the vaccine’s e ff ectiveness in preventing disease when R 0 > 0. The Jacobian matrix and the Routh-Hurwitz theorem are used to derive the aforementioned analysis techniques. The results are further examined numerically by using the standard second-order Runge-Kutta (RK2) method. In order to visualize the global dynamics of the aforementioned model, the proposed model is expanded to examine some piecewise fractional order derivatives. We may comprehend the crossover behavior in the suggested model’s illness dynamics by using the relevant derivative. To numerical present the results, we use RK2 method


Introduction
COVID-19 is a contagious disease that can be easily transmitted from one person to another.At the end of 2019, Wuhan, China, was the site of this terrible infection for the first time.Later, this disease gave rise to a worldwide outbreak, which the World Health Organization (WHO) notified in the first quarter of 2020 [1].Almost all of the countries of the world have experienced widespread transmission of this virus.The aforementioned sickness has currently claimed the lives of almost 50 billion people.COVID-19 has affected more than 500 million people worldwide.We cite [2,3] for a detailed explanation of how diseases are contracted and spread.Because the pandemic has affected the entire world, every nation has implemented its own policies.Some nations have implemented lockdowns, while others have issued laws requiring strict safety measures like wearing face masks and avoiding large gatherings, etc.The countries most impacted by the disease are Brazil, Italy, India, Iran, Spain, USA and the UK.Earlier research suggested that the virus originated in animals, but later research showed that it may also transfer from person to person [4,5].
Numerous infectious diseases are being investigated for potential treatments, controls, cures, etc.Since epidemiology is an important field of medical science, one of the key subfields from this viewpoint is mathematical biology.In contrast to earlier times, researchers are now becoming increasingly interested in biomathematics and bio-math engineering.In order to give a framework for comprehending the dynamics of diverse infectious diseases, mathematical models have been utilized.The most effective tools for studying different diseases are mathematical models [6][7][8].We can create numerous predictions, regulating strategies, etc [9][10][11], for a disease in a community by using mathematical models.Considering these factors, researchers have shown a strong interest in the aforementioned subject, we refer the reader to [12][13][14][15].Many mathematical models have recently been developed by researchers for COVID-19 as well; we provide e.g., see [16][17][18][19].Major public health difficulties and an economic catastrophe are being faced internationally as a result of the aforementioned disease's recent pandemic status and quick global spread.
In order to theoretically improve existing mathematical models for the aforementioned disease, it is necessary to evaluate the influence of recently developed vaccinations with known high efficacy that have weak transmission-controlling measures.For example, Pearson et al. [20] studied COVID-19 immunization in low as well as middle-income countries.they found it to be very cost-effective, according to recent research [21], if a vaccination with high efficacy and low cost is available.This will have a big impact on disease eradication.Prior to pharmaceutical interventions like treatment and immunization, non-pharmaceutical precautionary measures like the self-quarantining of confirmed cases, isolation, using of face masks, hand washing, social distancing, lockdowns, avoiding gatherings, and closing schools have been employed.This, however, is not a permanent solution.The best strategy will thus be to develop a suitable immunization and make it broadly accessible in every civilization at a fair price.Investigation of COVID-19, various researchers have developed many models.In order to do this, we have aimed to develop a mathematical model with a class of people who have received vaccinations and qualitatively analyze it while considering the therapeutic choices, immunization of those who are susceptible, and healthcare for infected/hospitalized people.Our model takes into account certain important biological and epidemiological elements of the disease under study, such as the disease's inhibitory effect, mortality rates from infections and natural causes, birth rate, and various vaccination rates.
Here, we build the dynamical model for COVID-19.Additionally, we depict the interactions between several model compartments in flow charts in Figure 1.Consequently, we first explain our model as follows:  Vaccination rate of quarantine In Figure 1, the model is presented as a flowchart.Epidemiological models have been examined using difference equations or classical order derivatives.Some references are cited as [22,23].The aforementioned operators cannot generate the phenomenon's global dynamical behavior because of their local nature [24][25][26][27].As a result, academics have now taken the majority of models into account when considering fractional order derivatives [28][29][30][31].The growth of the aforementioned field is attributable to its numerous applications to practical issues, particularly in the fields of mathematical biology and dynamical analysis [32][33][34][35][36].As far as we are aware, fractional calculus is being used to describe an increasing number of real-world scenarios; we refer the reader to [37][38][39][40][41][42][43].Additionally, the idea of fractional calculus has been applied in a number of scientific and technological fields [44][45][46].Because of its global nature and memory-preserving methodology, the theory of fractional calculus has been widely used in the mathematical modeling of numerous diseases for further information, see [47][48][49].The conventional analytical and numerical approaches have been upgraded to handle situations with fractional order derivatives (see references [50,51]).To examine diverse dynamical issues involving fractional order derivatives, the classical decomposition, perturbation, and transformation methods, for instance, have been expanded (some references are cited as [52,53]).Traditional numerical techniques, such as RKM tools, Euler and Taylor methods, the Adams-Bashforth method, etc., have been upgraded more and more to cope with fractional order problems (see [54,55]).Recently, a novel idea has been identified that describes the crossover behavior of the dynamics, because the states of multiple evolution processes frequently change suddenly.Due to these effects, typical derivatives, whether classical or fractional, cannot be used to appropriately address the relevant multi-phase behaviors.Therefore, scholars [56][57][58][59][60][61] have suggested this idea, which has the potential to appropriately address the aforementioned behavior to a certain extent.As a result, we also attempt at numerical analysis for our studied mathematical model by looking at multi-phase behaviors under various fractional orders by utilizing some numerical technique based on Newton interpolation polynomials.Discussion and a graphic presentation versus some actual facts are offered.
The provided model (1.1) can be stated as follows in the sense of the new derivative: where PWC D +0 denotes piecewise classical or Caputo derivative with two sub-intervals in [0, T ].
For further clarification, we respectively express the left hand sides of (1.2) as follows: where D +0 and C D φ +0 are classical and Caputo derivatives, respectively.Our manuscript is structured as follows: In the first section, we give introduction.In the Section 2, we give basic results.Stability results are given in the Section 3. Numerical results are given in the Section 4. In the Section 5, the fractional order form of the proposed model is investigated.In the Section 6, we given conclusion.

Basic results
We recollect some fundamental results of fractional calculus.Definition 2.1.[43] Let ϖ > 0, then, the non-integer order integral of a function S : [0, ∞) → R is defined as provided that right-hand side exists.Also, derivative in the Caputo sense is defined as is given by Here, we recollect the piecewise definitions of the fractional order derivative and integral.
t is the classical derivative for 0 < t ≤ t 1 and fractional derivative for t 1 < t ≤ t 2 .Definition 2.4.[61] Let h be a differentiable function, then, the classical, and fractional order piecewise integration is given as such that PF I +0 φ(t) denotes the classical integration for 0 < t ≤ t 1 , and Reimann-Liouville integration for t 1 < t ≤ t 2 .
Lemma 2.5.[61] Consider the problem with the piecewise fractional order derivative PFC 0 Here, we define the Taylor series from [59]; consider where The generalized Taylor series is thus given by 2) The outcomes of fractional calculus are used here.
. By using (2.2), after rearranging the terms, we obtain the general formula shown below. where Assume that φ = 1; the standard method described by (4.2).

Stability analysis
We initially examine our model's feasibility in this part.In order to achieve this, we add each equation in the model (1.1) under the condition that From above, we conclude that dN dt < 0 if N > N 0 .Additionally, we have Disease free equilibrium points: We analyze whether the system (1.1) has an equilibrium point.E 0 = (S 0 , 0, 0, 0) = B d o , 0, 0, 0 is a representation of the model (1.1) disease-free equilibrium.
Endemic equilibrium: From (1.1), lets assume that dS dt = 0, dE dt = 0, dI dt = 0, dQ dt = 0 and dV dt = 0; then, by simple calculation one can find the endemic equilibrium to be as follows:

Expression for R 0
In epidemiology, the idea of the fundamental reproduction number, or R 0 , outlines how illnesses spread and are treated.Both the frequency of the illness in the populace and the best countermeasures to protect the neighborhood's residents from the deadly virus are disclosed by R 0 .The most recent method for determining R 0 is as follows.If we set χ = (E, I), we have from the system (1.
The Jacobian of H for the disease-free equilibrium is and for the disease-free equilibrium, the Jacobian of F is given by Hence We have Hence, the required R 0 is given by If and only if the reproduction number meets the requirement that R 0 < 1, then we may state that the disease-free equilibrium is stable and the population is disease-free.You can halt the spread of an epidemic by exercising prudence.By sitting the numerical values of the variables above, R 0 = 0.98347.This number shows that the spread of COVID-19 in the locality is well managed.On the basis of R 0 , we established the following theorem.Theorem 3.1.(i) If R 0 ≤ 1, then positive equilibria for system (1.1) are not possible.(ii) If R 0 > 1, then there is an endemic equilibrium, often known as a distinct positive (unique) equilibrium E * = (S * , I * , Q * , V * ).Some further results on equilibrium point with global and local stability are presented here for system (1.1).
Proof.The Jacobian matrix at E 0 is given below: From a simple straightforward calculation M 0 has three eigenvalues which are negative; which means that (3.4) has two roots,i.e., one positive root and one negative root.Therefore, the disease-free equilibrium E 0 has an unstable saddle point.When R 0 = 1, to get the result for the disease-free equilibrium E 0 that is globally asymptotically stable, we consider the function called the Lyapunov function.
Taking the derivative of L(E, I) with respect to time t gives Suppose that I = 0 implies that dL(E,I) dt = 0. Furthermore, when R 0 = 1, which implies that E 0 is globally asymptotically stable in Φ. □ In the next theorem, we study the stability at point E * .
Proof.For the system (1.1), the Jacobian matrix is .
We see that the eigen values of M' i.e., Λ 1 = −d 0 and Λ 2 = −(d 0 + v Q ) are negative.Also, we have where Consider the following equation Using the Routh-Hurwitz theorem, we have that A 1 .A 2 − A 3 > 0 by a simple direct calculation which shows that E * is locally asymptotically stable.This completes the proof.□

Numerical simulation and discussion
Using the parameters from Table 2, we conducted the numerical simulations for model (1.1), we describe them in this section.Here, we describe how we simulate our system (1.1) by using the secondorder Runge-Kutta (RK2) method as implemented in [64].For this, we took into consideration Here § = (S, E, I, Q, V).The suggested approach allows us to express the relevant iterative formula as Here, K 1 = hF (t i , x i (t i )), also, h = t i+1 − t i .We simulated our results by using the numerical approach described in (4.2) and use the numbers from Table 2 in order to evaluate with our built-in model (1.1).We use the results of the numerical simulations to demonstrate our theoretical results.The parameter values and their descriptions for the numerical simulations are listed in Table 2.However, some parameter values that are unknown have been assumed for the sake of illustration.Vaccination rate of quarantine 0.0999 [65,67] We applied 300 days of the actual data from Pakistan.The overall population of the nation is approximately equal to 222.60 million, according to the data in [65,67].Nearly 47% of the nation's population was fully vaccinated during the past nine months, while 57% was partially vaccinated.In order to analyze the dynamics of disease propagation and the impact of the vaccination on its cure, we simulated the outcomes for 300 days as shown in Figures 2-6, respectively by using the numerical approach described by (4.2) of RK2 type.Here, we describe the model (1.1) .
Figures 2 through 6 depict the spread of the disease over a 300-day period at the specified speeds.The exposed class and the susceptible class as shown by the actual data that were employed.The infection spread relatively quickly in less than 50 days.As a result, there were also more persons under quarantine.More people had full or partial vaccinations during this time.We can notice that there were more persons who had received vaccinations.

Some numerical treatment via fractional calculus
Here, we develop a numerical scheme based on RK2 method for the model (1.2) as expressed in terms of (1.3).We have developed a numerical method for the two sub-intervals of [0, T ].We express the system (1.2) in terms of 1.3 by using (x) = (S, E, I, Q, V) as follows: (5.1) We may write the aforesaid system (5.1) by using (4.2) for classical order derivative when 0 < t < t 1 and (2.3) with the (2.4) for the fractional order Caputo derivative with t 1 < t < T .Hence, we can write the first equation of the system (5.1) as follows: where h = t i+1 − t i , and Similarly, the equations for the other compartments can be written as follows: and (5.7) Here, we plot the approximated results for two sets of fractional orders in two cases.Case I: When 0 < φ ≤ 0.6, the corresponding solutions exhibiting crossover behavior were obtained as in        Case II: Now we present the numerical results corresponding to different fractional orders such that 0.6 < φ < 1 in Figures 12-16, respectively which show the crossover behaviors of different compartments.Here, we discuss the results of modeling of the piecewise derivative shown in Figures 7-16.We observed a the sudden change after 150 days.So, for various fractional orders in [0, 50) and [50,300], we have plotted the numerical findings for systems (1.3) and (5.1).It is possible to see the declines and rises in several compartments with obvious crossover behavior.This novel idea makes a significant contribution to the understanding of how a system that experiences sudden changes in its dynamics behaves.Additionally, the fractional order derivative is a global operator that furthers our comprehension of the phenomenon.In Figure 17, we have compared some real data of infected cases with the simulated results at given fractional orders.Data applied here were taken from Pakistan from January 1, 2021 to January 30, 2021 (see [65]).We see that the simulated results were very similar to the real data which confirm the authentication of the fractional order model that we have considered.Also, we have compared the simulated data for the vaccinated class for the period of 11 days i.e., March 1, 2022 to March 11, 2022 (see [67]) in Figure 18.Infected cases per day in millions

Figure 2 .
Figure 2. Susceptible compartment's density over the specified time for the model (1.1).

Figure 3 .
Figure 3. Exposed compartment's density over the specified time for the model (1.1).

Figure 4 .
Figure 4. Infected compartment's density over the specified time for the model (1.1).

Figure 5 .
Figure 5. Quarantine compartment's density over the specified time for the model (1.1).

Figure 6 .
Figure 6.Vaccinated compartment's density over the specified time for the model (1.1).
for different fractional orders.

Figure 7 .
Figure 7. Graphical representations of the density of the susceptible compartment in the model (1.3) during the selected period for corresponding values of different fractional orders.

Figure 8 .
Figure 8. Graphical representations of the density of the exposed compartment in the model (1.3) during the selected period for different fractional orders.

Figure 9 .
Figure 9. Graphical representations of the density of the infected compartment in the model (1.3) during the selected period for different fractional orders.

Figure 10 .
Figure 10.Graphical representations of the density of the quarantining compartment in the model (1.3) during the selected period for different fractional orders.

Figure 11 .
Figure 11.Graphical representations of the density of the vaccinated compartment in the model (1.3) during the selected period for different fractional orders.

Figure 12 .
Figure 12.Density of the susceptible compartment in the model (1.3) during the specified period for various fractional orders.

Figure 13 .
Figure 13.Density of the exposed compartment in the model (1.3) during the specified period for various fractional orders.

Figure 14 .
Figure 14.Density of the infected compartment in the model (1.3) during the specified period for various fractional orders.

Figure 15 .
Figure 15.Density of the quarantined compartment in the model (1.3) during the specified period for various fractional orders.

Figure 16 .
Figure 16.Density of the vaccinated compartment in the model (1.3) during the specified period for various fractional orders.

Figure 17 .
Figure 17.Comparison of real and the simulated data plots at the given fractional orders for infected class.

Table 1 .
The parameters of the model (1.1).

Table 2 .
The model's (1.1) parameters and their corresponding values.