key: cord-0940017-s1fp3ddt
authors: DarAssi, Mahmoud H.; Safi, Mohammad A.; Khan, Muhammad Altaf; Beigi, Alireza; Aly, Ayman A.; Alshahrani, Mohammad Y.
title: A mathematical model for SARS-CoV-2 in variable-order fractional derivative
date: 2022-02-03
journal: Eur Phys J Spec Top
DOI: 10.1140/epjs/s11734-022-00458-0
sha: c783da659d28a19a67c8f196edccc2f9bd41488b
doc_id: 940017
cord_uid: s1fp3ddt

A new coronavirus mathematical with hospitalization is considered with the consideration of the real cases from March 06, 2021 till the end of April 30, 2021. The essential mathematical results for the model are presented. We show the model stability when [Formula: see text] in the absence of infection. We show that the system is stable locally asymptotically when [Formula: see text] at infection free state. We also show that the system is globally asymptotically stable in the disease absence when [Formula: see text] . Data have been used to fit accurately to the model and found the estimated basic reproduction number to be [Formula: see text] . Some graphical results for the effective parameters are drawn for the disease elimination. In addition, a variable-order model is introduced, and so as to handle the outbreak effectively and efficiently, a genetic algorithm is used to produce high-quality control. Numerical simulations clearly show that decision-makers may develop helpful and practical strategies to manage future waves by implementing optimum policies.

The corona virus affected and effecting the populations throughput the world by providing severe infection and death. The early and effective strategy to overcome this infection needed to follow the rule and regulations suggested by the World health organization to prevent the people from further reinfection. It is well known that every researchers from their respective field of research provided great efforts for the elimination of this infection. Many countries were the infection provided many deaths with so many waves of infection. Pakistan is one of among those where there are some waves occurred for this COVID-19 infection from which the population suffered. It has been reported the number of reported cases and deaths in each wave and seems the increase in the cases. There has been reported the number of death percentage in wave 1 and 2 to be 1 %.

The formulation of mathematical models in epidemiology and engineering areas has proven his role to be important; see [1] [2] [3] . Many researchers in different pera e-mail: altafdir@gmail.com (corresponding author) spective studied the coronavirus infection. While the mathematicians and biologists who worked on the modeling of the COVID-19 infection studied and formulated the mathematical models with real cases of different countries and presented about the peak of the infection and its eliminations of infection. In this regard, some mathematical models are highlighted here that formulated for this infection are [4] [5] [6] [7] [8] [9] [10] . The initial cases reported in China for the COVID-19 infection have been investigated in [4] through a mathematical model. The real cases from Spain and Italy are considered in the modeling of infection in the form of an SEIR model which has been studied in [5] . The COVID-19 infection dynamics and its prediction through a mathematical model of an SEIR have been investigated in [6] . The data from Italy are considered, and a mathematical model is formulated for the COVID-19 infection and its analysis has been considered in [7] . The impact of social distancing and other features that can be considered essential for the decrease in the COVID-19 infection has been studied by the authors in [8] . They formulated a mathematical model for the South African infected cases. A comparison study for the dynamics of coronavirus infection has been discussed in [9] . The number of real cases from Mexico population through a mathematical model is considered in [10] . The real infected cases from KSA are considered by the authors using mathematical modeling approach and obtained the results regarding the disease eliminations in the country [11] . The modeling and the contrail measure of infection of COVID-19 have been discussed in [12] . The disease eliminations and its control through optimal control theory are suggested by the authors in [13] . Addressing the model with lockdown measure for the COVID-19 infection is discussed in [14] . Some other interesting work related to the modeling of COVID-19 and related results regarding the disease are suggested in [15] [16] [17] . Moreover, the authors in [18] [19] [20] [21] [22] considered the COVID-19 infection models in fractional derivatives and some recommendations regarding the minimization of infection are given in the form of lockdown and control strategies. Moreover, the applications of fractional calculus and the numerical methods for the physical and biological problems have been discussed in [23] [24] [25] [26] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] . For example, the HIV infection model with antiviral drug therapy is considered in [23] . A fractional SEIR model using wavelet method has been proposed in [24] . Some biological models with analytical and numerical techniques are discussed in [25] . A fractional model in ABC operator in astrophysics is discussed in [26] . A recent study is conducted for the understanding of coronavirus infection and its controls [27] . A mathematical model in non-singular order to investigate the coronavirus infection is considered in [28] . Nonetheless, some studies show that in some conditions, time-varying fractional-order models are able to better demonstrate behavior of systems [40] [41] [42] [43] [44] [45] [46] [47] . However, now, there are few studies on time-varying fraction diseases dynamics, and it demands more studies in this field. Motivated by this in the current research, we propose and study a variable-order fractional model of SARS-CoV-2 dynamic.

In this work, a new mathematical model in the presence of hospitalization class is considered to better understand the COVID-19 infection in the presence of hospitalization and then using the concept of variable fractional order to study its dynamics. First, we will study the model in integer order and study the detailed mathematical results, and then, we extend the proposed model to have the results for the variable fractionalorder model. We discuss in details the formulation of the model in Sect. 2. The stability results and the model endemic equilibria have been shown in Sect. 3. Estimations of the model parameters and its numerical results are discussed briefly in Sect. 4 and 5, respectively. In Sect. 6, using the approach of variable order and its related results are discussed. All the results are summarized and concluded in Sect. 7.

Due to the COVID infection and its emerging in the society that producing many cases again and again in the society and creates a lot of difficulties for the humans as well as the government. Therefore, the eradications of the infections and to prevent the populations further, government present policies in light to World health organization. Therefore, here, our aims by formulating a novel mathematical model in the presence of hospitalization infection to understand in an effective way the dynamics of the COVID-19 infection with real cases reported from Pakistan. We denote the population of humans by N (t) and split into it further into six compartments: The susceptible S(t), exposed E(t), infected with visible symptoms I(t), infected with no visible symptoms A(t), the individuals that are hospitalized/ quarantined H(t) and those recovered from infection is R(t). We present the model below in the form of differential equations based on the above discussions

subject to the initial conditions that are non-negative. The parameters given in the model (1) are defined as:

The populations of healthy people that are susceptible to get corona virus infection are generated through the birth rate Π, while the natural death rate for each compartment in the model is shown by ω. The contact rate that generates the infection is given by ρ, while the infection that generated through the individuals with no clinical symptoms is given by the parameter σ. The death due to disease in the compartments of E, I, A, and H is shown by ω 0 , ω 1 , ω 2 , and ω 3 . The parameters ψ and τ distribute the infection in symptomatic and asymptomatic. The infected people are hospitalized at the rate ν. The recovery rate from the infected, infected with no clinical symptoms, and hospitalized people is shown by γ i for i = 1, 2, 3, respectively. We write the model (2) in the absence of the last equation, given by

We have the total dynamics for the system (2) shown by

Eur. Phys. J. Spec. Top.

Thus

and further

We consider the region given below for the system (2)

which is feasible biologically and all the solutions contained in it.

The equilibria for the system (2) in case if there is no infection in the population denoted by E 0 are shown below

Next, we compute the basic reproduction number R 0 , using analysis given in [48] , and present the results given by

where l 1 = (ψ + ω + ω 0 ), l 2 = (ω + ω 1 + γ 1 + ν), l 3 = (ω + ω 2 + γ 2 ), and l 4 = (γ 3 + ω 3 + ω). The following expression is the required basic reproduction number for our considered system and is given by

In the following, the stability analysis for the considered system (2) at the disease-free equilibrium E 0 is carried out.

The system (2) at the disease-free case E 0 is locally asymptotically stable whenever R 0 < 1.

Proof Computing the Jacobian at the disease-free case E 0 is given by

We get easily the two roots −ω and −l 4 . The remaining eigenvalues containing negative real parts can be determined through the equations given by

where

It can be observed that f i for i = 1, 2, 3 are positive whenever R 0 < 1, and furthermore, it can be obtained easily the conditions of Routh-Hurtwiz criteria, f i > 0 and f 1 f 2 − f 3 > 0, for i = 1, 2, 3. Therefore, it can be concluded that the system given in (2) at E 0 when R 0 < 1 is locally asymptotically stable.

We state and prove the following theorem to show that the system (2) at E 0 is globally asymptotically stable.

The system (2) at E 0 is locally asymptotically if R 0 ≤ 1.

Proof Define the Lyapunov function given by

where ξ m > 0 for m = 1, 2, 3, 4, that can be adjusted their values at later stage. Taking d dt of L, we have

Using the equations from (2) and making some setting, we have

Let us consider the values for the constants ξ 1 , ξ 2 , and ξ 3 . We choose as follows: ξ 1 = l 2 , ξ 3 = ρ, and ξ 3 = ρσl 2 /l 3 . Then, we have the following results finally:

Here, R 0 < 1, and so, we have L (t) is negative and so the model (2) is globally asymptotically stable at E 0 .

Here, we present and obtain the endemic equilibria of the system (2) by describing its endemic equilibrium by E 1 = (S * , E * , I * , A * , H * ), and so, we have the below

Consider the above values in the following equation:

getting the following:

We can see that Φ 1 > 0, while Φ 2 can be positive depending on the value of R 0 . When R 0 < 1, so we getting I * = −Φ 2 /Φ 1 , and thus, there is no equilibria. The requirements for the positive endemic equilibria are to have R 0 > 1, and hence, for the system (2), a unique positive endemic equilibrium exists.

To understand effectively the dynamics of an epidemic model, it is necessary or useful to have the real data to demonstrate the model parameters and to obtain reasonable results for the epidemic model. In this regard, we consider the cases from March 06, 2021, till April 30, 2021 for the finding of parameter values; we can refer the readers for the real data [49] . The population of Pakistan considered here is N (0) = 220000000, and we estimate the parameters Π and ω as follows: Π ≈ 8903 per day using the expression N (0) = Π/ω, and ω = 1/(67.7 × 365), where ω = 1/(67.7) is the average life span in Pakistan. The parameters other than these can be obtained through the model data fitting. The data considered in this paper from March 06, 2021to April 30, 2021. We use the curve fitting technique to get the values of the parameters, as shown in Table 1 , and the possible fitting of the model to the data is shown in Fig. 1a with the estimated basic reproduction number R 0 ≈ 1.2036. Figure 1b describes the behavior of the model versus data for long time level. It follows from Fig. 1(a) that the model behaves good to the real cases. Thus, the model parameters considered in Table 1 are realistic and can better study the dynamics of the COVID-19 infection.

Here, we consider the system (2) with the consideration of the cases reported from March to the end of April in Pakistan to understand the parameters impact on the model. The cases consider here are in daily basis, so the time unit is to be per day. In numerical results, we consider the parameters values shown in Table 1 . We assuming that that there is no recovery from the infection yet and so R(0) = 0, and no enough information about the asymptomatic and hospitalized data, so we consider A(0) = H(0) = 0. We consider the disease contact rate ρ and σ with different values on the dynamics of infected people shown in Fig. 2. In Fig. 2a , we see when decreasing the value of ρ, the number of infected people decreases faster. It can concluded that the prevention from the infection using the instructions suggested by World Health organization should be follows to decrease the disease burden from the populations. Washing hands, social distances, using masks, avoid gathering, self-quarantine, and restrictions to home can better reduce the infection. In Fig. 2b , one can see that the infection decreases when decreasing the value of σ. The parameter σ that denote the contacts of healthy people with asymptomatic infected people (individuals do now show clinical symptoms) can decrease the infection risk in the population. It should be noted that asymptomatic infected does not show visible disease symptoms and can increase the infection. It is documented that the number of infected individuals due to COVID-19 has high percentage of infection due to asymptomatic infection. Therefore, it is important for the government to make increase the testing facility to identify the number of asymptomatic infected people. The impact of the parameters ψ , τ and ν has been shown graphically in Fig. 3 . Decreasing the values of these parameters, we have decrease in the number of infected people. It has been shown in the literature that fractional calculus is able to produce more accurate results in modeling natural systems and processes. As a result, in this part, we present a novel variable-order fractional model of SARS-CoV-2 dynamics in Pakistan. Also, we consider an additional component, namely, hospitalized people (H). The fractional model with variable order is given as follows:

γ H is the recovery rate of quarantined infected individuals and δ I is the rate of symptomatic infected individual's transition to the quarantined infected component. Also, q(t) indicates the time-varying fractional-order derivative. In this study, the Caputo fractional derivative based on the predictor-corrector method is applied for numerical calculations of the fractional model. It is one of the most widely used methods for analyzing fractional differential equations in a chaotic manner (see [52, 53] for more details on the predictor-corrector method). To investigate the effects of the time-varying fractional derivative, we present the results of the model with two different fractional derivatives. Figure 4 shows the number of individuals when q 1 (t) = 0.9 and q(t) = 0.9 + 0.1 1+exp(−t) . As it is evident from this figure, when the value of fractional derivative varies by time, the results are different, and it can dramatically affect the behavior of the model.

Herein, we propose a genetic optimization approach to find optimal decision rules. Genetic optimization is an appropriate technique to handle such problems, since there is a trade-off between all possible values. Table 3 lists the parameters and configuration of the applied multi-objective genetic algorithm.The fractional-order derivative is considered to vary over time as q(t) = 0.9+

We consider δ I as the decision variable. The following is the cost functions that have been chosen:

We have chosen this cost function to reduce the accumulated infected people. In what follows the results of the optimization are provided. Figure 5 shows the value of the cost function from the genetic algorithm optimization. As it is shown in this figure, we could choose the optimal controllers, which Fig. 6 . As can be seen, the maximum number of infected people has decreased when compared to the case where no optimal control effort was made. Also, the decreasing rate of infectious people is noticeably quicker.

Here, in the present investigation, a attempt has been made to study a new mathematical model for the dynamics of novel coronavirus. We first formulated a new mathematical model for the better understanding of the virus by considering the hospitalized class. After formulation the model, we presented in details the model analysis for the integer-order case and obtained the stability results in detailed. We found the local and global asymptotical stability of the model when R 0 < 1.

Eur. Phys. J. Spec. Top. Fig. 6 The number of individuals with optimal control Furthermore, we estimated the parameters using the real cases of coronavirus from Pakistan and obtained the results for the data fitting. The realistic parameters for the considered data gave the basic reproduction number R 0 = 1.2036. Using further parameter values, some graphical results are obtained that indicate the infection can be minimized by following the reduction in contact rate, asymptomatic contact, rapid testing of individuals to isolate the individuals with no clinical symptoms. The numerical results for the integer cases can be useful for the minimization of infection by following the suggestions given in figures in numerical results section. Finally, a variable-order fractional model of SARS-CoV-2 in Pakistan was introduced. Then, through an evolutionary algorithm, the control of the variable-order fraction model was studied. Through numerical simulations, it was shown that using the proposed optimal strategy, the number of infected people has decreased significantly, and the third wave is controlled more efficiently.

A new Hepatitis B model in light of asymptomatic carriers and vaccination study through Atangana-Baleanu derivative

Heat and mass transfer in MHD Williamson nanofluid flow over an exponentially porous stretching surface, Case Stud

Comparative analysis of (Zinc ferrite, Nickel Zinc ferrite) hybrid nanofluids slip flow with entropy generation

Modeling the dynamics of novel coronavirus (2019-ncov) with fractional derivative

A modified seir model to predict the covid-19 outbreak in spain and italy: simulating control scenarios and multi-scale epidemics

On an interval prediction of covid-19 development based on a seir epidemic model

Monitoring italian covid-19 spread by a forced seird model

Modelling the potential impact of social distancing on the covid-19 epidemic in south africa. Computational and mathematical methods in medicine

Mathematical models for covid-19 pandemic

Modeling and prediction of covid-19 in mexico applying mathematical and computational models

Dynamics of fractional order covid-19 model with a case study of saudi arabia

Mathematical modeling for novel coronavirus (covid-19) and control

Modeling the impact of non-pharmaceutical interventions on the dynamics of novel coronavirus with optimal control analysis with a case study

Modelling the spread of covid-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Chaos

Global stability and cost-effectiveness analysis of covid-19 considering the impact of the environment: using data from ghana

Novel fractional order sidarthe mathematical model of covid-19 pandemic

Dynamic analysis of a mathematical model with health care capacity for covid-19 pandemic

A fractional-order model describing the dynamics of the novel coronavirus (COVID-19) with nonsingular kernel

Threshold condition and non pharmaceutical interventionss control strategies for elimination of COVID-19

On the necessity of proper quarantine without lock down for 2019-nCoV in the absence of vaccine

A fractional-order SIRD model with timedependent memory indexes for encompassing the multifractional characteristics of the COVID-19

A novel fractional mathematical model of COVID-19 epidemic considering quarantine and latent time

Analytical and numerical study of the HIV-1 infection of CD4+ T-cells conformable fractional mathematical model that causes acquired immunodeficiency syndrome with the effect of antiviral drug therapy

Analytical and numerical solutions of mathematical biology models: the Newell-Whitehead-Segel and Allen-Cahn equations

A numerical algorithm for the solutions of ABC singular LaneEmden type models arising in astrophysics using reproducing kernel discretization method

Mathematical modeling and optimal control of the COVID-19 dynamics

A vigorous study of fractional order COVID-19 model via ABC derivatives

Incorporating fast and intelligent control technique into ecology: a Chebyshev neural network-based terminal sliding mode approach for fractional chaotic ecological systems

Fuzzy adaptive control technique for a new fractional-order supply chain system

A Caputo-Fabrizio fractional-order model of HIV/AIDS with a treatment compartment: sensitivity analysis and optimal control strategies

Deep recurrent neural networks with finite-time terminal sliding mode control for a chaotic fractional-order financial system with market confidence

Artificial macro-economics: a chaotic discrete-time fractionalorder laboratory model

Spectral entropy analysis and synchronization of a multi-stable fractionalorder chaotic system using a novel neural network-based chattering-free sliding mode technique

Coexisting attractors in a fractional order hydro turbine governing system and fuzzy PID based chaos control

Fractal-fractional neuro-adaptive method for system identification

A new RBF neural network-based fault-tolerant active control for fractional time-delayed systems

Tracking control and stabilization of a fractional financial risk system using novel active finite-time fault-tolerant controls

A new fractionalorder hyperchaotic memristor oscillator: dynamic analysis, robust adaptive synchronization, and its application to voice encryption

Comments on time-varying fractional order

Synchronization of different fractional order chaotic systems with timevarying parameter and orders

Electric vehicle power propulsion system control based on time-varying fractional calculus: implementation and experimental results

Controlling and synchronizing a fractional-order chaotic system using stability theory of a time-varying fractional-order system

Intelligent parameter identification and prediction of variable time fractional derivative and application in a symmetric chaotic financial system

On the variable-order fractional memristor oscillator: Data security applications and synchronization using a type-2 fuzzy disturbance observer-based robust control

On the development of variable-order fractional hyperchaotic economic system with a nonlinear model predictive controller

King algorithm: a novel optimization approach based on variable-order fractional calculus with application in chaotic financial systems

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

A dynamical study of SARS-COV-2: a study of third wave

A new predictor-corrector method for fractional differential equations

A predictorcorrector approach for the numerical solution of fractional differential equations

Acknowledgements The authors extend their appreci-