Results in Physics

In this paper, we develop a new mathematical model for an in-depth understanding of COVID-19 (Omicron variant). The mathematical study of an omicron variant of the corona virus is discussed. In this new Omicron model, we used idea of dividing infected compartment further into more classes i


Introduction
The infection casus by Omicron variant of COVID-19 is called Omicron virus.In the last week of November 2021, this new invariant is identify in South Africa.After identification and repaid growth the virus transmits to other countries.As compared to other variant of corona virus omicron is not severe but transmission rate is very high.According to (CDC) Center for Diseases Control says that any one having infected by Omicron virus can transfer the infection to other vaccinated or individuals which do not have any common symptoms [1].Some known symptoms of Omicron are runny nose, body ache, congestion, fatigue, cough etc.In order to minimize the infection rate other countries of the universe stop their flight to infected country South Africa.
The dynamics of corona virus have been studied by many scientists and researchers to control the spread of disease in the population [2][3][4][5].To more effectively minimize the infection in population in future spreading [6][7][8] the scientist of the filed are tried to produce the immune and vaccinated the majority of population.Although new SARS-CoV-2 varieties have emerged with the passage of time, infections are still a problem in many nations.This paper discusses some integral order mathematical model for investigating SARS-CoV-2 infections.For instance, [9] discusses the early infection of disease in China by using a very comprehensive model.By taking into account the actual SARS-CoV-2 cases, the best elimination and control of the disease in Pakistan has been researched in [10].The best and most efficient technique to minimize the infection is isolation and quarantine, which has explored by the authors in [11] using a mathematical modeling approach.The SARS-CoV-2 infection can spread to other uninfected persons extremely quickly.An analysis of the lockout and its effects on disease prevention using mathematical modeling can be found in [12].The authors provided the disease control scenario for a model SEIR approaches utilizing real numerical data from France and Italy [13].Due to its numerous attributes and practical applications to issues in engineering and physics, fractional calculus is attracting the interest of researchers all over the world.The memory, the crossover behavior, and the hereditary qualities can be shown in a fractional order model.In integral differential equations [14], the construction of the operators and applications to epidemiology [15][16][17], applications to follow dynamical equations [18,19], etc., the fractional calculus with various fractional operators has been found.Researchers [20] explored the SARS-CoV-2 stochastic model along time delay.Furthermore, authors [21] have used a model of fractional order to examine the singular and nonsingular kernels of the SARS-CoV-2 infection.[22] discusses the Caputo fractional model for the investigation of SARS-CoV-2 infection.The authors of [23] developed a mathematical model that takes the media campaign role into account while tackling the SARS-CoV-2 infection.[24] studies the modeling of the SAES-CoV-2 in China.The investigation and analysis of SARS-CoV-2 non reported cases are done in [25].Authors in [26] examined how to use the new generalized idea of Caputo fractional order differentiation to arrive at the numerical analysis of the SARS-Cov-2 infecting model.[27] has examined the SARS-CoV-2 disease model with integral and non-integral orders.The time-delayed COVID-19 infection model has been proposed in [28].Researchers [29] have looked at a fractional model to pinpoint when the illness peaked in Brazil.[30] discusses the numerical analysis of the SARS-CoV-2 model.The authors of [31] took into account the affected instances in India, created a fractional model, and then got their results.Authors [32] have discussed a SARS-CoV-2 infection model using immunization.The infection cases in Argentina have been studied using a fractional model proposed in [31].The authors published their numerical simulations and two vaccine models for the SARS-CoV-2 infection in [32].To analyze the infection cases in Spain, a fractional-order SARS-CoV-2 model with optimal control strategies was applied [33].The authors examined the fractional-order SARS-CoV-2 model infection using the idea of the modified Euler approach in [34].The conditions for the SARS-CoV-2 model's global stability in the situation of no disease were proposed by the authors [35].A nonlinear fractional order SARS-CoV-2 model infection has been proposed in [36].To explore the qualitative analysis of the model, a mathematical model for SARS-CoV-2 has been built in [37].Moreover authors [38] have investigated a model on the SARS-CoV-2 for global dynamics.In additions, the concepts of fractional calculus have been used very well recently to investigate COVID-19.Additionally, various models involving concepts of mathematical analysis have been utilized very well.In this regards we refer some work, where various analysis and numerical scheme have been used to study different kinds of problems as [39][40][41][42][43][44][45][46][47][48][49][50].
The primary objective of this paper is to use mathematical modeling to comprehend the behavior of the new SARS-CoV-2 type known as the Omicron.We develop the model using the omicron feature and parameterize it using real data from South Africa.In order to explore the potential existence of several layers or waves, a NCOVID-19 mathematical model for new variant (omicron) is taken into consideration.Sensitivity analysis is used to identify the variables that have the greatest potential to improve or reduce the fundamental reproduction number  0 .
This study present epidemics model for SARS-CoV-2 regarding focus on strategies, control and specific vaccination of population.Proposed model are solve for Reproduction Number, analysis and model description respectively.Numerical simulation are calculated by NSFD Scheme [51].A general group of techniques in numerical analysis known as non-finite difference schemes create a discrete model to provide numerical solutions to differential equations.In the last conclude this paper by 'Conclusion'.

Model formulation
Due mainly to the new SARS-CoV-2 version, also called Omicron virus, which was first discovered in South Africa, people there once more had to adhere to tight SOPs, keep social distances, wear face masks, and other limitations.Because of the new SARS-CoV-2 variant's rapid proliferation in several countries, these nations prohibited their citizens from traveling to South Africa.We are developed a new model to comprehend the dynamical behavior of the COVID-19 cases from South Africa while keeping in mind the Omicron variation.We take into account the entire population, indicating it by  (), and  is divided in  =  +  + .Further classifying it into six distinct epidemiological groups: Those who have been exposed to (), represent the uncontagious population, and () is stand for exposed population, which coming into touch with asymptomatic population, symptomatic population, or omicron variant-infected individuals; asymptomatic population are represent (); and symptomatic population by () (infected individuals population although they have no clear symptoms sign), symptomatic people (population who exhibit clinical symptoms consistent with SARS-CoV-2 infections), having Omicron variant infection is represent by () (who exhibit clear clinical symptom consistent with omicron infection; and have the ability to spread the disease to others whether or not they have received the vaccine), the recovered population is represent by, () (population recovered from (), () and ()).We obtain our model as The flow of infection and different rate from one compartment to another is presented as in Fig. 1.Symbols involved in the model (1) are described in Table 1.

Feasible region, positivity and boundedness
On addition of all equations of the model (1) yields ≤  −  0  From (2), we have All values of non-negative for .
Therefore, all solutions began of the system (1) will remain nonnegative for all value of t is equal to zero.So, the system (1) is well-posed mathematically.
The feasible region for the dynamical analysis is following

Endemic equilibrium (EE)
Also the EE is computed as

Expression for 𝓡 𝟎
In epidemiology there is a factor, which called the basic reproduction number ℛ 0 which describe the control and transmission of infection.Furthermore from ℛ 0 , we can determine how the infection of COVID-19 is going in the population and which is the greatest choice to control the infection from the population.The next generation matrix method is used to find ℛ 0 is given below let  = (, , , ), then form system (1), we have Where .
Also for the DFE, Jacobian of matrix  is follow Hence,  0 is spectral radius of   −1 computed as .
Hence,  0 has three parts, first part represent contact between infected population with uninfected (healthy) population.While second part represent the contact between uninfected population to asymptomatic and third part show infection of symptomatic population.The infection of SARS-CoV-2 can be control in the population whenever the valve of  0 < 1, on the other hand the infection will be spread in the community if  0 > 1.In next theorem, we discus the stability of the system (1) under condition of  0 > 1 and  0 < 1 at disease-free equilibrium.
Using the numerical value of Table 2, we give a 3D profile of  0 in Fig. 2.

Stability analysis
This section has the study of properties of equilibrium point with local and global stability of the system (1).Also, we will investigate the properties of these equilibrium and global analysis of the system (1).

Local stability
In this section, we present the local stability of the model (1) with condition  0 < 1.

Theorem 2. ''If ℛ 0 < 1, infection model COVID-19 is stable locally asymptotically at DFE  0 .''
Proof.To obtained the stability result, the Jacobian Matrix at  0 is given below The characteristic polynomial of the Jacobian matrix at DFE is given by ( 0 − ) = 0, where  is the eigenvalue and I is 6 by 6 identity matrix.Thus,  has eigenvalues given by 1 ,  2 ,  3 ,  4 ,  5 and  6 are negative strictly.Thus, Disease free stability (DEF) is stable locally asymptotically.□

Global stability
In this part of our manuscript, we present global stability at DFE of the system (1).We constructed a function called Lyapunov function to show the model ( 1) is stable asymptotically globally in the next theorem.
Theorem 3. '' If  0 ≤ 0, then infection of system (1) is stable globally asymptotically at  0 .''Proof.To prove this theorem result, we constructed a function called Lyapunov function. ( Here  1 ,  2 ,  3 and  4 are some constants.Now w.r.t ''t'' take derivative of above equation. Putting the values from ( 1) After some basic calculation We obtained Hence,  0 ≤ 0 then the system (1) is globally asymptotically stable.□

Numerical results and discussion
We take into account the infected cases that were recorded in South Africa.The infection are evaluated on a regular basis using a unit per day.Calculated natural birth rate and death rate 0.0155432, which is the mean rate of life of South Africans in 2020-21, are two of the 12 parameters in the given model (1).The method (0) =   0 is used to calculate the birth rate, here (0) represent total population of South Africa in 2021.About (0) = 60140000 people are projected to live in South Africa in 2021.(0) = 60069540 is the uninfected population individual in the non-appearance of the disease, and (0) = 0.062  million, (0) = 0.008 million, (0) = 100, and (0) = 360 are the reported people on 1st November, 2021, with (0) = 0 being the initial given conditions susceptible to fitting data.Let suppose that recovery from infection is zero yet.In order to suitable the system to the indicated numerical data with the time mentioned early, we used the nonlinear least square method.The experiments were carried out up until the required level of accuracy in the model fitting.Here, we establish NSFD scheme [51] for our considered model.Therefore, consider first equation of (1) as Which is decomposed by NSFD scheme as Like (7), by utilizing NSFD method, we present the model (1) as Developing NSFD scheme [51], we plotted the model (1).To understand the dynamics of the system, we used real data of South Africa to simulate.With using numerical value table, we present graphically the dynamical behaviors of various compartments of our model in Figs.3-5 using classical order model.

NSFD Scheme for Fractional Order Model (1)
In this part of our paper, we present some fractional concept for model (1), we used Grunwald-Letnikov method for Caputo derivative as used in [41,42].Moreover, the numerical method used in this paper has the ability to demonstrate the dynamical properties of model (1).Taking a small interval ℎ, our using numerical method preserve the stability of disease-free and endemic equilibrium points.Mostly researchers used different methods like Euler, RK2, and RK4 method but our numerical method have the advantages over it.For instance, the said numerical scheme has less complexity than RK2 and RK4 methods.This numerical scheme has a significant advantage over other scheme that for a large time raising it yield a good numerical results.Definition 4. ''Integral of non-integer order  > 0 of a function  ∶ [0, ∞) → R is defined as provided the integral exists at the right sides''.Moreover, by definition of Caputo, we have From Riemann Liouville operator integral with some fractional order as Following the procedure of [42], we have where Here it is interesting if  = 0, we have from ( 12) From Eq. ( 12) to discretion (13) as we have ] , where  = 0, 1, 2, … , (15) where . The stability rule is obtained from some functions as sin h, cos h, h, etc.Here, we confer our counsel system (1) to non-integral order 0 <  ≤ 1 as follow Considering of (15) and by utilizing Grunwald-Letnikov discrimination method, for system (16) for non-integral order, we present numerical method here.] . ( Here in Figs.6-8, we present graphically the dynamical behaviors of various classes using different fractional orders.

Conclusion
In this manuscript, we constructed a new mathematical SEIVR model for NCOVID-19, which is base on behavior of virus.Our main goal is to estimate the effect of immunization on population.In the presence of presenting no clear infection or symptoms and loss of immunity are the root measure which differentiate this virus from other modeled infectious diseases.On the bases of analytical expression for  0 which is a main element to determine necessary and sufficient conditions for disease-free and endemic equilibrium.We testify our theoretical results with real data for a case of South Africa.
Our main resulting and finding are follow.Vaccination helps to overcome the disease, we provide a mathematical proof that transmission rates are reduce with vaccination.We mathematical justified that the increasing in number of infection in South Africa was followed by vaccination campaign.On the bases of collected data from the parameters of model, we proved that the rate of vaccination for SARS-CoV-2 control spread of disease efficiently.
In near past, the researchers [52] of the fields observe that this virus change their behavior each time.As we know that this virus change their behavior time to time.We are trying to develop a new model for observed new pattern as like beta and omicron in future.

Fig. 3 .
Fig. 3. Dynamical behaviors of susceptible and exposed classes of the considered model corresponding to classical model.

Fig. 4 .
Fig. 4. Dynamical behaviors of symptomatic and symptomatic with no visible symptoms classes of the considered model corresponding to classical model.

Fig. 5 .
Fig. 5. Dynamical behaviors of omicron infected and recovered classes of the considered model corresponding to classical model.

Table 1
Nomenclature and discerption of the model (1).