key: cord-0822637-9wpxo3jw
authors: Alzahrani, Ebraheem; El-Dessoky, M.M.; Baleanu, Dumitru
title: Mathematical modeling and analysis of the novel Coronavirus using Atangana–Baleanu derivative
date: 2021-04-26
journal: Results Phys
DOI: 10.1016/j.rinp.2021.104240
sha: 5194bcbb3e174b9d4c8f5a9e50ce55d030474766
doc_id: 822637
cord_uid: 9wpxo3jw

The novel Coronavirus infection disease is becoming more complex for the humans society by giving death and infected cases throughout the world. Due to this infection, many countries of the world suffers from great economic loss. The researchers around the world are very active to make a plan and policy for its early eradication. The government officials have taken full action for the eradication of this virus using different possible control strategies. It is the first priority of the researchers to develop safe vaccine against this deadly disease to minimize the infection. Different approaches have been made in this regards for its elimination. In this study, we formulate a mathematical epidemic model to analyze the dynamical behavior and transmission patterns of this new pandemic. We consider the environmental viral concentration in the model to better study the disease incidence in a community. Initially, the model is constructed with the derivative of integer-order. The classical epidemic model is then reconstructed with the fractional order operator in the form of Atangana–Baleanu derivative with the nonsingular and nonlocal kernel in order to analyze the dynamics of Coronavirus infection in a better way. A well-known estimation approach is used to estimate model parameters from the COVID-19 cases reported in Saudi Arabia from March 1 till August 20, 2020. After the procedure of parameters estimation, we explore some basic mathematical analysis of the fractional model. The stability results are provided for the disease free case using fractional stability concepts. Further, the uniqueness and existence results will be shown using the Picard–Lendelof approach. Moreover, an efficient numerical scheme has been proposed to obtain the solution of the model numerically. Finally, using the real fitted parameters, we depict many simulation results in order to demonstrate the importance of various model parameters and the memory index on disease dynamics and possible eradication.

Infectious diseases are disorders caused by different organisms such as viruses, bacteria, fungi or parasites. A number of outbreaks of different infectious diseases have been reported across the globe taking the lives of millions of humans each year. These outbreaks are not only a serious issue for public health, but also having a significant impact on the whole societies, economic interest and political systems and is felt in both developing and developed countries across the world. The world is facing different types of deadly infectious diseases. Although some of these infections have been eradicated from the world or almost wiped out still many diseases such as HIV/AIDS and TB persist with little or even no hope of getting them under control. In addition, the emergence of new infectious diseases is another serious threat to the world. The novel Coronavirus infectious disease or COVID-19 is one of the newly emerged deadly infectious disease.

It was initiated in the mainland China during the end of 2019 and spread to the rest of the word within no time. Many facts are still to explore about the novel COVID-19 pandemic. Initially, the transmission of the COVID-19 infection was observed from animals to humans. The first virus transmission among humans was observed in Guangdong a province of China on January 20 [1, 2] . New research confirmed that this virus can also be transferred to humans through environment. As on 28 August 2020, the Who reports show the cumulative cases confirmed with this infection are higher than 24,021,218 with the total death cases were 821,462 and the recovered cases are more than 15,907,858 [3] . The KIngdom of Saudi Arabia (KSA) is among those Arabian countries seriously affected from the COVID-19 pandemic and has the highest graph of the infected cases. The first case confirmed with COVID-19 is reported in Saudi Arabia on 2nd of March 2020. Later on after few days on 14th of March, the second case was detected when a companion came from Iran. Currently, Saudi Arabia is placed at fourteen position in the list of high confirmed infected cases. As of 28 August, the Kingdom has reported 311,855 confirmed COVID-19 cases, with 286,255 recoveries and 3,785 deaths.

Recently, many approaches are used for the analysis and investigations of the Coronavirus infection and its early eliminations. In this regards, the modeling tools are gaining much attention from researchers point of view. The mathematical modeling is effectively and confidently being used to provide to help the public health decision-making in order to set effective control measure for diseases. These models are helpful to study the dynamical behavior and spatiotemporal patterns of infectious disease and consequently, can be used to set effective control measures to eradicate the disease. Mathematical models are mostly based on differential and integral operators with integer or non-integer orders. A rich literature regarding deterministic compartmental epidemic models have been formulated on the dynamics of effective controlling strategies of the ongoing novel COVID-19 pandemic in different countries and territories in the world. The importance of some control measures on the dynamics and future prediction of COVID-19 outbreak in Pakistan has been presented in [4] . The authors in [4] also have developed the optimal control 2 J o u r n a l P r e -p r o o f problem using optimal control concepts. A deterministic mathematical model incorporating the incidence of infection through environmental route in Ghana is suggested in [5] . The authors carried out the cost effective analysis and used the confirmed COVID-19 cases in Ghana to estimate the model parameters. A new transmission model on COVID-19 describing the dynamical behavior of COVID-19 has been studied in [6] . Additionally, the authors [6] simulated the proposed model for real parameters estimated from infected cases in Nigeria. The effective mitigating measure for the COVID-19 in the population of Canada have been presented in [7] .

Mathematical modeling of real-world physical problems with fractional derivatives gains much attention from the researchers and scientists in recent days. It is due to some serious limitations of classical integer order derivative, which can be overcome through fractional order operators. The crossover behavior and fading memory effects found in many physical and biological processes can be captured only using fractional order derivatives. The model with fractional operators provides a better understanding of disease dynamics. In the recent era, a class of fractional operators based on either singular or nonsingular kernels were studied in literature. But, among these, the operators which gain more interest are known as Caputo derivative developed in [8] , the Caputo-Fabrizio (CF) recently described in [9] and the Atangana-Baleanu (AB) operator developed in the sense of Caputo shown in [10] . These operators were effectively utilized to express many problems in science and other field as well [11, 12, 13] . The Atangana-Baleanu derivative has a nonsingular as well as nonlocal kernel and is the most recent and generalized fractional operator.

The application of ABC derivative can be found in almost every field of science such as [14, 15, 16, 17] and references therein. Recently, a number of the epidemic transmission models describing the COVID-19 dynamics with fractional order have been developed and studied in the literature. The transmission of COVID-19 pandemic using a mathematical model in ABC sense are explored in [18] . Moreover, in the study [18] , the authors have taken the Wuhan infected cases for parameter estimation in order to make a realistic study. Baleanu et al. [19] analyzed transmission procedures of COVID-19 with the help of a compartmental epidemic model in CF sense. The application of the Caputo operator to describe the dynamical behavior of COVID-19 through a new model is recently studied in [20] . A realistic and novel approach based on fractional-fractal operator is used to construct a new compartmental model based on the dynamics of the current COVID-19 pandemic in [21, 22] . The application of well-known classical Caputo-type operator coupled with real data to studying the transmission dynamics of novel COVID-19

can be found in [23, 24] .

In the current study, a fractional epidemic model with environmental transmission of the virus that contributes the further spread of infection among humans population is considered. The model in integer case was studied recently in [25] . The model is developed by taking into consideration the environmental contributions of the latent, infected and asymptomatic infected population. The model under consideration is taken in the from of fractional order differential operators instead of integer order. The ABC derivative is used to derive the proposed epidemic compartmental model and obtain their results. In the first stage, the model construction is briefly described using the integer order differential equations. The classical COVID-19 model is then simulated to provide a better fit to COVID-19 cases in KSA and ultimately to estimate the updated parameters. The present paper has been arranged in various sections as follows: The related concepts of fractional operator and its integral has been suggested in Section 2. The novel COVID-19 model for the real data analysis of Saudi Arabia in integer case is presented in Section 3. The model in fractional case is developed in the 4th Section, while its necessary mathematical aspects are studied in Section 5. The Section 6 presents an iterative scheme of the considered model and the simulation results along with discussion, which are performed in Section 7. Finally, the conclusion is given in Section 8.

The related concepts regarding the fractional operators used in the given paper are shown in the following:

Definition 1. Let y(t) ∈ C n , then the fractional derivative of Caputo type with order ρ in (n − 1, n] such that n ∈ N defined in [8] is given by:

(1)

Clearly, C D ρ t (y(t)) tends to y ′ (t) when ρ → 1.

Definition 2. We define the Atangana-Baleanu derivative defined in [10] and shown by the following:

where, 0 ≤ ρ ≤ 1 and y(t) ∈ C[a, b]. The expression AB(ρ) in the above integral equation defines the normalized function.

Definition 3. The Atangana-Baleanu derivative defined above has the following integral [10] :

Let us express the model proposed in the fractional derivative of ABC case by the following way [16] :

where Υ : D ⊂ R n × R + −→ R n and ρ ∈ (0, 1]. So, to obtain the global asymptotical stability of the model in the Atangana-Baleanu sense, the results for (4) can be used. 3. Integer order model

The present section explores the modeling of the novel Coronavirus. We begin the modeling process be denoting the host population by N (t) dividing into five mutually exclusive epidemiological classes based on dynamics of COVID-19 infection. These classes consist of susceptible S, exposed E, infectious showing symptoms of disease I, infected with no disease symptoms I a and the individuals recovered are denoted R respectively. So, that (iii) Individuals in the exposed class within the incubations which may or may not yet show any disease signs are capable to transmit the disease further.

(iv) The COVID-19 related death is only considered in symptomatically-infected class.

(v) The exposed, symptomatic, asymptomatic COVID-19 population has contributed the virus to the environment and is one of the transmission routes of the infection.

Taking the above discussion and assumptions in account, we formulate the following dynamical model for the COVID-19:

where the additional non-negative initial conditions

In the model described in (5), the individuals in the susceptible class is generated by the birth rate Π, while the death naturally in each class can be shown by µ. The death rate due to COVID-19 infection is denoted by ζ, which only appears in symptomatically-infected class. The parameters η 1 , η 2 and η 3 are the effective transmission rates of infection due to exposed, symptomatically-infected and asymptomatically-infected individuals, respectively. The parameter η 4 denotes the generation of infection due to environment. We denote the incubation period of the individuals by κ. The exposed individuals at the end of incubation period, a proportion τ , where (0 < τ < 1), remain in the I which develops symptoms compartment while the rest join I a class which have no or mild disease symptoms. The parameters δ 1 and δ 2 account for the recovery rates of symptomatically and asymptomatically infected people. The contribution of the virus to the environment due to exposed, and the population of both infected comparers (i.e. I and I a ) are shown respectively by ̟ 1 , ̟ 2 and ̟ 3 . The virus can be removed from the environment using the parameter φ.

We define the following expressions to make the model easier for onward analysis:

Thus, the model (5) can be reformulated as follows:

The estimation of parameters is one of the essential parts for the mathematical models in epidemiology.

Various approaches are used for the estimation of parameters of a biological model having real data. Among these techniques, the nonlinear least square curve fitting is one of them. We used the nonlinear least square curve fitting method in order to determine the best values of the parameters involved in the model under consideration. We use the method of nonlinear least square curve fitting method by considering the infected cases of KSA for a selected period of time as shown in Figure 1 . In parameters estimation, we get the value of the parameter Π and µ using the literature. The remaining parameters have been gotten through the fitting. The procedure used in this paper for estimation is obtained from the work in [4] , while the readers can see the details procedure of this technique in [4] . After using the technique, the desired result for the data fitting is shown in Figure 1 while, Table 1 

The classical integer order epidemic models developed for an infectious having its own importance in epidemiology and its importance cannot be ignored. However, there are certain number of limitations of the models developed via classical differential equations. Such as the non existing of the memory or nonlocal effects, not able to capture the crossover behavior of a physical or a biological process. In the result, the mathematical models developed via the integer order derivatives are not suitable in some cases due to the properties mentioned above. To deal with limitations, the fractional operator specifically the ABC operator comprises the memory effects and the crossover behavior of the model. Therefore, to explore Cumulative COVID-19 infected cases the COVID-19 dynamics more realistically, the proposed model (7) described in Section 2 is reformulated with the replacement of classical derivative by the one having fractional order in ABC sense. Thus, the fractional epidemic model for COVID-19 with the nonlocal kernel is formulated through the following system:

The notation ABC D ρ t represents the time fractional derivative in ABC sense and the parameter ρ ∈ (0, 1] denotes the arbitrary order of the ABC operator. 

The current section presents the basic mathematical features of the ABC Coronavirus transmission model as formulated in (8) . Initially, we investigate the biologically-feasible region for the ABC model.

After this, we evaluate the model equilibria, the basic reproductive quantity and then explore the asymptotic stability of the model at the disease free equilibrium. Finally, the existence and uniqueness are provided. We proceed as follows:

Now, summing the first five equations of the system (8), we are lead to the following equation

Solving equation (9) via Laplace transform, we deduce that lim t→∞ N (t) ≤ Π/µ.

As a result, we derived the biologically feasible region as shown in the following set:

J o u r n a l P r e -p r o o f Journal Pre-proof

This section explores the equilibrium points of the fractional ABC case described in (8) . Also, we evaluate the basic reproductive number R 0 and then prove the asymptotical stability of the disease free equilibrium. Moreover, we will explore impartment features of fractional models known as the existence and uniqueness (EU). The fractional model describing the COVID-19 dynamics (8) posses usually two types of steady states. The one is known as the disease free equilibrium (DFE) and the second one is named as the endemic equilibrium (EE). The DFE is obtained by solving the system by equating the right sides of (8) to zero at the infection free state. Thus, the following expression is derived for DFE of the model (8):

In order to derive R 0 of the model, we follow the next generation matrix approach [27] . Using the approached mentioned earlier, we have the desired results for our system (8) given by:

Thus, using the well-known criteria developed in [27] , we finally, evaluate the desired expression as follows:

We focus on the existence of EE for the epidemic model (8) in ABC sense in this section. The EE is denoted by K * * 1 = (S * * , E * * , I * * , I * * a , R * * , B * * ),

where, the expressions involved are obtain by solving the equations in (8) at steady state.

where,

After substituting the values involved (10) into (12), we obtain the following equation in λ * * :

where;

Thus, the following result regarding the existence of EE for the COVID-19 (8) is parented. 

For the proof of the local asymptotical stability (LAS) of the model around the DFE, the necessary condition to be shown is that for the entire eigenvalues of the matrix (14) lie outside the closed angular sector |arg(Π)| ≤ 1 [28] . Although the Matignon's result was applied to the Caputo case, but it can be extended to systems having ABC fractional order (FO) derivative. To proceed, let us evaluate the corresponding Jacobian matrix of (8) as:

Next, we introduce the following result for the desired stability of the model.

satisfy the necessary condition:

Proof. It is clearly observed from J K0 that the two eigenvalues i.e., −µ and −µ are negative and clearly satisfy the necessary condition (15) for all ρ ∈ (0, 1]. Moreover, after some manipulations, it is easy to show that the remaining four eigenvalues will have negative real parts whenever, R 0 < 1 and thus, following [28] , the arguments of the remaining eigenvalues fulfil the desired condition as described (15) . Hence, K 0 is LAS whenever, R 0 < 1.

Journal Pre-proof

The GAS of the COVID-19 model (8) at the DFE is discussed in the following result. The result for this purpose is stated as follows:

Proof. We introduce the following suitable Lyapunov function in order to proceed

where, C j , for j= 1, · · · , 4, describe unknown positive constants to be chosen later. The time fractional ABC derivative of F (t) leads to the following

Using the system defined in (8), we get

Let us choose the constants values as

and then after some simplifications, we have

It is clear that when R 0 < 1 then ABC D ρ t F (t) is negative, due to the non-negativity of all parameters. Thus,with the help of Theorem 2.1, that the DFE K 0 is GAS in the region Φ. This subsection is focused on proving the EU of fractional compartmental epidemic model (8) . The well-known fixed point theory is utilized for the said purpose. Firstly, the COVID-19 model (8) is shown through a general initial value problem as follows:

In the problem described in (16) 

Additionally, the function F satisfies the Lipschitz condition expressed as below:

Next, we introduce the following theorem regarding the EU of the model solution.

Theorem 5.4. There exists a unique solution of the epidemic model (16) , if the condition shown in (18) holds:

Proof. Taking the ABC integral operator over the problem (16) , which gives a nonlinear Volterra integral equation in the following form:

Let J = (0, T ), and the operator F : C(J , R 6 ) → C(J , R 6 ) be defined as

Thus, the expression in (19) takes the form:

Further, let . J be the supremum norm over J and expressed as:

It is clear that C(J , R 6 ) with the respective norm . J constructs a Banach space. Also, the inequality given in (23) can be easily demonstrated as:

with v ∈ C(J , R 6 ), F(t, x) ∈ C(J 2 , R) such that

Using equation (21), we obtain

Further, utilizing equations (17), (23) along with triangular inequality, the equation (25) leads to the following form:

Finally, we get the following equation:

where,

If the condition given in (18) holds then the operator F will be a contraction. Thus, the desired model shown by (16) will possess unique solution.

It is obvious that the biological models are often complicated due to their nonlinearity. The nonlinearity together with fractional derivative the model becomes more complex and difficult to solve analytically. In this regard, the numerical techniques are often considered for its solution graphically. Different numerical methods are considered for the solution of fractional order numerically. Among these methods, the method developed in [29] based on Lagrange interpolation polynomial of two-step is appropriate for the given problem. We will present the details of the scheme in the present section and then write the scheme for the proposed system. To get the approximate solution of the problem, we use the method in [29] . Considering the Atangana-Baleanu integral described in (3) and applying upon the general problem described in (16), then we have the following:

J o u r n a l P r e -p r o o f

Now, utilizing the polynomial interpolation concept, we approximate the function F(ς, v(ς)) by two-step

Lagrange polynomial over the closed interval [t j , t j+1 ] as follows:

Thus, equation (29) becomes:

By simplifying the above equation (31), the final approximate form can be shown by the below expression:

Using the above scheme (32), we can write the system (8) below that presents the numerical scheme: 

We derive the numerical scheme using the two-step Lagrange interpolation polynomial for the approximate solution of the model involves fractional derivative. Then, we particularize the scheme for our model by giving the scheme shown in (33). As we taken the daily reported cases, so the time unit considered in the numerical solution in days. The respective values of the biological parameters involved in the model obtained through the real cases of KSA data shown in Table 1 are used for the simulation purpose. Firstly, we depict the influence of variation in the memory index (or fractional order ρ) upon the dynamics of all population classes of the proposed model. After that, we analyzed the influence of various model parameters such as ̟ 3 and φ on the COVID-19 incidence in the population of Saudi Arabia.

The dynamics of the susceptible population for various values of fractional order ρ of the ABC derivative are shown in Figure 2(a) . It is observed that susceptible individuals are increasing with the decrease in ρ. The impact of memory index ρ over the exposed, symptomatically and asymptomatically COVID-19

infective individuals is shown in sub-plots (b), (c) and (d) Figure 2 , respectively. It is observed that peaks of population in the exposed, symptomatically and asymptomatically infected classes are decreasing significantly for the smaller values of ρ. The dynamical behavior of the recovered people is depicted in Figure 2 (e). The recovered people shows a decreasing effect with a decrease in the fractional order ρ. The impact of memory index ρ on the dynamics of virus concentration in the environment B is analyzed in Figure 2 (f). Further, Figure 3 (a-d) depicts the influence of the variation in the viral contribution due asymptomatically-infected individuals ̟ 3 on the COVID-19 incidence. We have analyzed this impact for various values of fractional order ρ. In all cases, it is observed that the peaks cumulative symptomatically and asymptomatically infected curves decreased well with the decrease in the parameter ̟ 3 . Moreover, it is observed that the decrease in pandemic peaks is comparatively faster for smaller values of ρ as shown in Figure 3(b-d) . Finally, the impact of parameter φ for various values of ρ is shown in Figure 4 (a-d).

It is observed that infected curves become flatten with the increase in parameter φ by different rates to their baseline values as shown in Figure 4(a) . Furthermore, the decrease in the peaks of symptomatically and asymptomatically infected curves is slightly faster for smaller values of fractional order ρ showing the significance of the memory index on disease dynamics.

In the last few decades, mathematical models have been used as a useful tool to understand the complex dynamics and to determine the future trend of an infectious disease. Although, the epidemic models developed via classical integer order differential systems have its own significance to explore a disease dynamics. But the mathematical model designed through the fractional operators are more useful than the ordinary case, due to data fitting, the memory effects and the crossover behavior etc. We studied a new fractional model for understanding the complex dynamics of Coronavirus in the Kingdom of Saudi Arabia reported cases using non-singular kernel operator. We studied the model and provided its mathematical results in details. We found that the given model is locally asymptotically stable at the infection free state when its threshold below 1. We used the nonlinear least square method for the obtaining of parameters estimations and then provide the good fit to the proposed data. We also proved the existence and uniqueness criteria of the fractional COVID-19 model. The proposed parameters obtained through estimation are considered further to obtain its numerical results. We studied graphically the model equations solution with memory index and found that the solution converges for the arbitrary order. We then 

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

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

Analysis of a mathematical model for covid-19 population dynamics in lagos, nigeria

Quantifying the role of social distancing, personal protection and case detection in mitigating covid-19 outbreak in ontario, canada

Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications

A new definition of fractional derivative without singular kernel

New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model

A fractional model for the dynamics of tuberculosis infection using caputo-fabrizio derivative

New aspects of poor nutrition in the life cycle within the fractional calculus

Chaos in a simple nonlinear system with atangana-baleanu derivatives with fractional order

Non validity of index law in fractional calculus: A fractional differential operator with markovian and non-markovian properties

Modeling and analysis of the fractional hbv model with atanganabaleanu derivative

Stability and lyapunov functions for systems with atangana-baleanu caputo derivative: An hiv/aids epidemic model

A robust study on 2019-ncov outbreaks through non-singular derivative

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

A fractional differential equation model for the covid-19 transmission by using the caputo-fabrizio derivative

A mathematical model for covid-19 transmission by using the caputo fractional derivative

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

The dynamics of covid-19 with quarantined and isolation

A fractional order mathematical model for covid-19 dynamics with quarantine, isolation, and environmental viral load

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

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

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

Stability results for fractional differential equations with applications to control processing

New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models