Dynamical analysis of coronavirus disease with crowding effect, and vaccination: a study of third strain

Countries affected by the coronavirus epidemic have reported many infected cases and deaths based on world health statistics. The crowding factor, which we named "crowding effects," plays a significant role in spreading the diseases. However, the introduction of vaccines marks a turning point in the rate of spread of coronavirus infections. Modeling both effects is vastly essential as it directly impacts the overall population of the studied region. To determine the peak of the infection curve by considering the third strain, we develop a mathematical model (susceptible–infected–vaccinated–recovered) with reported cases from August 01, 2021, till August 29, 2021. The nonlinear incidence rate with the inclusion of both effects is the best approach to analyze the dynamics. The model's positivity, boundedness, existence, uniqueness, and stability (local and global) are addressed with the help of a reproduction number. In addition, the strength number and second derivative Lyapunov analysis are examined, and the model was found to be asymptotically stable. The suggested parameters efficiently control the active cases of the third strain in Pakistan. It was shown that a systematic vaccination program regulates the infection rate. However, the crowding effect reduces the impact of vaccination. The present results show that the model can be applied to other countries' data to predict the infection rate.


Introduction
The acronym SARS stands for severe acute respiratory syndrome. Initially, it was discovered in Asia in 2003. From Asia, it spreads to North and South America, as well as Europe. It was, however, contained in 2004. High fever, headache, body pains, diarrhea, dry cough, and pneumonia are symptoms of SARS caused by the coronavirus. The word corona means ''crown.'' There is a small crown of spike proteins visible. The spike proteins bind to the ACE2 receptors on our respiratory system's cells. After the cells combine, the coronavirus injects its viral RNA into our cells, which is coded to produce more proteins and more virus cells. SARS cov2, also called COVID -19, is the cause of the current pandemic. It began in December 2019 in Wuhan, China. The virus was previously found in a wild animal, most likely a bat, and it was discovered in humans via tests conducted during a pneumonia outbreak in Wuhan. It then traveled all over the world, eventually reaching the USA on January 20, 2020. They can also be transmitted via air droplets and are capable of causing difficulties in breathing. Other symptoms can be mild or severe. Everyone is vulnerable, but it affects the elderly the most. People with underlying medical issues, such as cancer or obesity, are housed in senior living homes or rehab facilities. The immunocompromised individuals and those of similar nature are susceptible, and it is possible for those who recovered from the first infection to be infected again. The global death rate is low, yet there are still a lot of people dying. The USA has one of the highest infection rates in the world. On March 11, 2020, it was declared a pandemic with an R0 (reproduction number) of approximately 2.5. The number of fatalities per 100,000 persons in most countries affected by the pandemic is on the increase. Brazil has the highest mortality rate of 264.6 per 100,000, while Vietnam has the lowest death rate of 2.15 per 100,000. With 187.15 deaths per 100,000, the USA is in the lead. COVID -19 is a unique type of virus. The strain is SARS COV2, and symptoms occur two to fourteen days after exposure. As a result, it is easier to spread. COVID -19 has a higher death rate than the flu. It is more infectious. If it persists from weeks to months, one may lose the sense of taste or smell, resulting in long-term consequences such as lung problems. A polymerase chain reaction (PCR) test or an antigen test is available to figure out if one is infected or not. Antigen tests and PCR tests are two of the most common infection tests now available. The antigen test is quick and inexpensive to carry out; however, it is up to 20% inaccurate. The PCR test, on the other hand, is highly potent in terms of detection. It simply takes a few cells to detect the presence of coronavirus. The results can be obtained in a matter of hours. Depending on the type of PCR test one gets, a good PCR test can be 100 percent accurate. A nasopharyngeal swab is inserted into the nose and used to clean the area leading to the throat. It does not affect one's brain regardless of how much it appears to be. The swab is then rotated in each nostril for 15 s to ensure enough germs or bacteria are picked up. An antibody test can be used to determine if one has ever been infected with the coronavirus. Antibodies indicate that the immune system has successfully attacked the virus. It can tell whether or not a vaccine is effective. If, however, one is newly infected, a quarantine or simple isolation (up to ten days) with medical care is adequate for mild cases. For severe cases leading to difficulties in breathing or organ failure, treatment in the (intensive care unit of a) hospital is recommended. Finally, we discuss the prevention for which social distancing, uses of masks, and sanitizers are vital components. We should also wash our hands frequently and get vaccinated. A distance of six feet from sick persons or other individuals is recommended inside and out outside buildings. The importance of masks cannot be overstated. Even in vaccinated areas, cases of old and new strains (variants) have been recorded. But masks are not required outside unless it is windy or crowded places where the spread is likely to be faster. Facemasks prevent the disease from spreading through the air to other humans and objects we may contact. After contacting any surface, we must wash our hands for at least 20 s before touching our face, going to the restroom, or preparing meals. Hand sanitizers containing at least 60% alcohol can be used if we cannot wash our hands. On COVID -19 vaccination, the spike proteins on the virus hold on to ACE2 and force the RNA into our cells, causing more COVID -19 virus to be produced. This process causes other cells to be infected. Our body's genetic material (mRNA) directs the secretion of the same proteins, the T-cells, and the lymphocytes, to target viral RNA. As a result, our bodies combat the virus. Pfizer vaccines are for people aged 12 and up, while Johnson & Johnson, and Moderna are for people aged 18 and over. The vaccines were released in phases, with phase 1A targeting healthcare workers and elders living in communities, followed by others with a lower priority classification. We must be aware that some persons may experience allergic responses to these immunizations. It is not common, but it can happen. According to some conspiracy theorists, COVID does not exist. They believe it is a government myth designed to keep us under control. One theory has it that masks give us COVID because of our inner COVID. Furthermore, that the use of facemasks will kill us due to our CO2. If this is true, the doctors who have been standing in a surgery room for hours would have been long dead. Some people speak of herd immunity. Vaccines, on the other hand, provide herd immunity. Only the immuned people can prevent the virus from spreading presented an S-I-R-M model to determine the COVIDrelated losses [6].
In 2021, Daniel investigated the SEIQCRW model on COVID diffusion with irregular factors of illness and the need for safety practices in Nigeria [7]. For estimation dynamics of the coronavirus disorder 2019 (COVID-19) outbreak, isolation security controls, Prathumwan et al. proposed a SLIQHR model in 2020 [8]. In 2021, Balike investigated the SEIHQR model for reducing the economic impacts of the COVID in the Congo [9]. In 2021, Raslan examined a SEHQIR model for COVID-19 infection projection in Egypt [10]. In 2020, Chen et al. proposed a BHRP model used to simulate the step in the process propagation of a new coronavirus [11].
In 2021, Sinaga et al. suggested an SEIR model in Indonesia and examined the coronavirus's structural analysis [12]. In 2020, Biswas et al. proposed using the SEAIQHR model to explore the COVID-19 contagion in India [13]. The use and misuse of mathematical modeling for communicable diseases management should be investigated, according to James et al. in 2021 [14]. In 2021, Ameen et al. proposed the use of SSLLIPD model to examine the kinetics of COVID-19 by applying a partial computational formula [15]. In estimating emerging coronavirus epidemiology, Uddin et al. proposed the SUQC model in 2020 [16]. In another development, Jiang et al. in 2020 studied a SEIIRMH model for determining the best SARS-CoV-2 removal plan in China, South Korea, and Italy [17]. In 2021, Kahn et al. used a model to guide COVID vaccination programs and testing procedures in nursing homes [18]. For the COVID -stability study employing fast-slow breakdown, Chen et al. developed a subtype of the SIR epidemic model in 2020 [19]. Staying at home, keeping distances, and early identification are the core points in controlling the spread of COVID -19, according to Kim et al. in 2020 [20]. Machado et al. studied the pandemic of coronavirus via complex systems that have characteristics that give rise to the emergence of rare and extreme events [28]. Rajagopal et al. investigated the fractional-order model to predict the dynamics of coronavirus outbreaks [29]. Quaranta et al. studied a multiscale territorial analysis of the pandemic using various models and data-driven approaches in Italy [30]. Different types of dynamical analysis are studied to model coronavirus-like diseases as presented in [31][32][33][34][35][36]. Many mathematical models have been analyzed with the help of different strategies, as illustrated by authors in [37][38][39][40][41][42]. The stability analysis of COVID-19 model is studied under the properties of fractional calculus in [44].
The design of our paper is as follows. In Sect. 2, we discussed the mathematical model and performed its analysis. Then, in the subsections, positivity, boundedness, existence, and uniqueness are examined. In Sect. 3, parameter estimation is presented. Finally, in Sect. 4, the numerical results are outlined to analyze the dynamics of the virus graphically. We then discussed the local and global stability of the model and gave concluding remarks in the last section.

Model formulation
By the theory of the population dynamics, the human population N t ð Þ is categorized into the four subpopulations like SðtÞ: (Susceptible population), IðtÞ: (Infected population), VðtÞ: (Vaccinated population), and RðtÞ: (Recovered or immune population). A dynamic of infection into the population is based on the law of mass action (the rate of change of reacting substance is directly proportional to the product of interacting substances) and nonlinear complex ordinary differential equations (ODEs). The transmission map of the model is shown in Fig. 1.
The physical relevance of parameters of the model is as follows: KN: (the recruitment rate of the population), bI: (the force of infection of virus), 1 1þaI : (the crowding effect of population on the virus), l: (the rate of mortality due to virus or natural of each subpopulation), d 1 : (the rate at which infected population got vaccination during the period of quarantine, or isolation, etc.), d 2 : (the rate at which susceptible population got vaccination under the program launched by World Health Organization (WHO)), c : (the rate at which infected population may recover due to its internal immunity and natural circumstances), and r : (the rate of doses in the population who recovered or got immune after vaccination). Deterministic modeling based on the following assumptions are as follows: Susceptible population directly vaccinated, considering the two types of immunity for an infected population like vaccination, and isolation or quarantine or natural immunity, and recovered population may not be infected again. Other types of interaction are ignored without loss of generality. The system of equations obtained from the transmission map of the virus is as follows: The total dynamics of the system (1) is obtained by adding the four equations as follows: and hence, The feasible region of the system (1) is defined in the following way 2.1 Properties Proof Let us start from the class I t ð Þ, For the function V t ð Þ; the following inequalities hold: We shall define the norm where D k is the domain of k. Using the above norm, the inequalities for the function S t ð Þ are defined as desired.
Theorem 2 (Boundedness) For any time t, the system (1) is bounded and lies in the feasible region Proof By letting the population function as desired.
We shall define the norm k 1 ¼ sup t2D k k t ð Þ j j, and we consider the Banach space [27]. We present here the existence and uniqueness of the solution piece wisely. To obtain such results, we need to verify growth and Lipschitz condition properties. Let us consider the four positive constants M 1 ; M 2 ; M 3 , and 8i ¼ 1; 2; 3; 4. We first verify that For proof, we consider the function f 1 S; I; V; R; t ð Þ , and the following estimations hold By using the same methodology, we get For the function f 3 ; . We have The condition d 2 For the function f 4 ; we have Therefore, the condition of linear growth is verified if \1; as desired.

Equilibria
We determine the equilibria of the system (1) by assuming that the state variables are constant and by putting the right side equal to zero. Equation (1) admits two types of equilibria as follows: (i) coronavirus-free equilibrium = C 1 ¼ S 1 ; I 1 ; ð

Reproduction Number
We determine the reproduction number of the system (1) by using the well-known results like the nextgeneration matrix method after substituting the value of coronavirus-free equilibrium. We get where: The spectral radius of q AB À1 ð Þ is denoted by

Strength number
The extension of the reproduction number is called the strength number. No doubt, the reproduction number has great significance in the field of epidemiology regarding the spread and extinction of disease. The reason behind different techniques used by the epidemiologist is to obtain such numbers (i.e., violation of uniqueness of reproduction number and many more). The primary outcome of this method relies on the prediction of the waves of the spread of disease. The critical thing is that the next-generation matrix method is used to evaluate this number by assuming the coronavirus-free equilibrium into the system (1) by taking the second derivative of infectious classes. Consequently, the transmission and transition matrices are represented by F and G [25], where: Here, A 0 ¼ À bK d 2 þl ð Þ \0 is called the strength number of the system (1). The remarkable conclusion includes the analysis of local maximum, local minimum, and inflection points. Thus, having a negative strength number is an indication that the system (1) will have a single magnitude, either a maximum with two infection points indicating a single wave or a rapid decrease from the coronavirus-free equilibrium. Thus, the infection will rise after a minimum point with the renewal process and then be stabilized or stopped later on, as desired.

Stability analysis
In this section, we test the local and global stability of the system (1), considering the two defined equilibria.
Theorem 3 (Local stability at C 1 Þ The system (1) at Proof The Jacobian matrix obtained from the system (1) is as follows. : The Jacobian matrix at C 1 is as follows Consider J À kI j j¼ 0, and hence Þcþd 1 þl ð Þ \1; R 0 \1: It is seen that C 1 is locally asymptotically stable (LAS), R 0 \1: Theorem 4 (Local stability at C 2 Þ The system (1) at Proof The Jacobian matrix at C 2 of the system (1) is as follows. .
Consider J À kI j j¼ 0, which implies Since A 0 ; A 1 both are positive when R 0 [ 1; by the Routh-Hurwitz criterion for the second order, the C 2 is locally asymptotically stable.
Theorem 5 (Global stability at C 1 Þ The system (1) at Proof By letting the Lyapunov function U : X ! R defined as.

U I ð Þ ¼ ln
Hence, the system is globally asymptotically stable at C 1 .
Theorem 6 (Global stability at C 2 Þ The system (1) at Proof By letting the Lyapunov function Z : X ! R defined as.
Þare positive constants to be chosen later. We have Hence, by Lasalle's invariance principle, C 2 is globally asymptotically stable (GAS) in X.

Second derivative theory of Lyapunov stability
We present an analysis of the second derivative of the associated Lyapunov function of the system (1) to understand the variability of the process [26].
Theorem 7 (Global stability at C 1 Þ The system (1) at It can be seen that For the realistic analysis and to predict the peak of the third strain of coronavirus in Pakistan, we need to use the real cases reported from August 01, 2021, till August 29, 2021, which is taken from the health ministry of Pakistan (www. Covid19. gov. pk) and presented in Table 1.
Before we parametrize the system (1), we need to calculate the fundamental values of the parameters, such as the population's birth rate and death rate. According to the world meter info, the total population of Pakistan is N 0 ð Þ ¼ 220,000,000 in 2021, while the birth rate of susceptible humans represented by K, and the death rate defined by l are shown to be K ' 8903 per day and l ¼ 1=ð67:7 Â 365Þ, with the number l av ¼ 1 67:7 as the average life span in Pakistan [24]. Furthermore, the values of the remaining parameters of the system (1) are presented in Table 2. Using the least-square curve fitting technique, we show the desired fitting in Fig. 2. The estimated value of the reproduction number is R 0 ¼ 1:3333. Hence, the desired values of the transmission rates are instrumental in studying the system (1) graphically.

Numerical results
We consider the system (1) with newly reported coronavirus cases in Pakistan to obtain the numerical results. Time is defined in days, and parameters values presented in Table 2 are fitted using the nonlinear least-square curve technique. The initial conditions      Figure 3 depicts the prediction to determine the peak of active cases. Figure 4 illustrates the regular residue of active cases. Figures 5, 6, 7, 8, and 9 outline the frequency of confirmed cases, deaths, recovered active cases, and percentage of the confirmed case in the pie chart. Figure 10 shows that systematic uses of vaccination may control the increase in the strain of coronavirus. On the other hand, Fig. 11 indicates that the crowding effect decreases the impact of vaccination on the population.

Nonstandard finite difference method
This section aims to provide a nonstandard finite difference discretization of the mathematical model (1). To this effect, we divide the interval [0, L] into M [ N subintervals, respectively, and step sizes h = L/M. The approximate solutions S; I; V; andR of (1) will be denoted as S n ; I n ; V n ; andR n , respectively, for each n ¼ 0; 1; . . .; N: Under the rules, the discretization of a system (1) is presented in [43]. We have The simulation of the system (17)- (20) is done by using objective data estimation. It reveals the beauty of this method which is routed for long-term behavioral analysis of the model. Also, the nonstandard finite difference method commits the model's dynamical properties like positivity, boundedness, consistency, and stability, as shown in Fig. 12.

Concluding remarks
A mathematical model has been developed to study the third strain of coronavirus in Pakistan. The infected cases were taken from August 01, 2021, till August 29, 2021. The model is asymptotically Fig. 11 Impact of crowding effect on active cases stable in the sense of local (if R 0 \1Þ and global (if R 0 [ 1Þ at both equilibria. Moreover, we determine the basic reproduction number, R 0 using the reported active cases during the stipulated time. The leastsquare curve fitting method is used to obtain the realistic parameters and the numerical results graphically. The effect of crowding on the dynamics of coronavirus is also shown. In the same way, due to the development in vaccination programs worldwide, it is one of the significant measures taken to overcome the menace of the new strain, and its efficiency is shown graphically using the developed model. Furthermore, we analyzed the graphical results and found that the increase in coronavirus cases of the third strain can effectively be controlled if the SOPs designed by the World Health Organization (WHO) are followed. We also observed that the active cases could be reduced in Pakistan as shown by the results from our model. Knowing the peak of the third strain is vital for any country and determines the maximum active cases on a particular day. These results can be helpful for planning and valuable to the Ministry of Health (MOH) and the decision-making authority: National Command and Operation Center (NCOC) in Pakistan. This type of analysis can be applied to data from neighboring countries such as India, Iran, and Bangladesh.
Funding We have no funding for this research article.

Declarations
Conflict of interest The authors declare that they have no conflicts of interest to report regarding the present study.
Data availability Data will be available on reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.