## Abstract

In December 2019, a new virus belonging to the coronavirus strain has been discovered in Wuhan, China, this virus has attracted world-wide attention and it spread rapidly in the world, reaching nearly 216 countries in the world in November 2020. In this chapter, we study the fractional incommensurate SIQR (susceptible, infections,quarantined and removed) COVID-19 model with nonlinear saturated incidence rate using Atangana–Baleanu fractional derivatives. The existence and uniqueness of the solutions for the fractional model is proved using fixed point theorem, the model are shown to have two equilibrium point (disease-free and an endemic equilibrium). Some numerical simulations using Euler method are also carried out to support our theoretical results. We estimated the value of the fractional orders and the parameters of the proposed model using the least squares method.. Further, the sensitivity analysis of the parameter is performed as a result, our incommensurate model gives a good approximation to real data of COVID-19.

## 1. Introduction

At the end of December 2019, an unidentified virus was found in Hubei province, China. The responsible virus was later confirmed as a new coronavirus [1]. The World Health Organization (WHO) named the virus as the 2019 novel coronavirus (2019-nCoV) or COVID-19 or 11 February 2020. The first confirmed case of the virus was discovered on 17 November 2019 in Hubei. As of November 3, 2020, more than 47.58 million cases have been reported in 216 countries in the world, which resulted in more than 1,215,180 deaths. Over 34.185 million people recovered [2].

A fractional differential equations are equations where the order of the derivative is a real or complex number. Its first appearance is in a letter written to Guillaume de l’Hôpital by Gottfried Wilhelm Leibniz in 1695. There are several types of fractional derivatives in the literature. We mention them Riemann–Liouville derivative, Caputo derivative, Riesz derivative, Grünwald–Letnikov derivative, Hadamard derivative, Marchaud derivative, Weyl derivative and recently, other types of fractional derivative have appeared, the most important of which are Caputo-Fabrizio derivative and Atangana-Baleanu derivative, the latter has aroused the interest of many researchers in various fields, due to its efficiency in modeling real problems in several areas, namely epidemiological problems. In recent years researchers have been interested in studying some real problems in various fields using fractional calculus, such as physical problems [1, 3, 4], engineering mechanics [5], epidemiological models [6, 7,8,9,10], image processing [11, 12], chaos theory [13, 14] and others.

After the emergence of the Corona pandemic, researchers in various fields tried to study the evolution of this virus, how it is transmitted and provide the appropriate mathematical model for it. Many epidemiological models with fractional (non-integer order) derivative have been developed to understand the transmission of 2019-nCoV and other infectious diseases outbreak from various aspects. We mention them, Singh et al. [15] have been studied the fractional epidemiological model for computer viruses using CF-fractional derivative. The comparison between the classical and fractional SEIR epidemic model of Ebola virus was studied by Area et al. [7]. Tulu et al. [16] developed the Caputo fractional mathematical model of the Ebola virus with a quarantine strategy and the presence of the vaccine by using MCMC method. Dokuyucu and Dutta [8] examined the new Ebola virus model with CF-fractional operators. Kumar et al. [17] developed Caputo–Fabrizio fractional SIRS-SI model for the malaria disease. Singh et al. [6] studied the diabetes epidemic model using the operators proposed by Caputo and Fabrizio. Baleanu et al. [10] studied a time fractional SEIR model for COVID-19 with Caputo-Fabrizio derivative and an auxiliary parameter by evoking asymptomatic infected people and the COVID-19 case in the reservoir. Abdo et al. [18] used Atangana-Baleanu-Caputo operator to study the mathematical model of COVID-19. Thabet et al. [19] suggested a fractional SEIHDR model for COVID-19 using ABC fractional derivative. Zhang et al. [20] applied the fractional SEIRD Model to the real data of the 2019-nCoV for China using Caputo fractional derivative. ud Din et al. [9] studied a general SEIARW fractional model for the transmission of COVID-19 by using the new fractional operators of Caputo and Fabrizio. Higazy [21] suggested a fractional SIDARTHE epidemic model for 2019-nCoV by using CF-fractional operators. Some other outstanding studies of 2019-nCoV by fractional derivative have been made in [22, 23, 24, 25].

Always when a new virus appears, researchers look for effective solutions and methods to control the new virus, including vaccination, isolation, and quarantine. In the absence of the vaccine, the isolation and quarantine strategies remain effective to mitigate and eliminate the impact of the virus (see [26,27, 28,29,30,31,32,33]).

In [34] the authors suggested changing the fractional SIAR model for HIV/AIDS with treatment compartment from classical commensurate model to incommensurate Caputo-Fabrizio fractional model by giving each equation a different order from the other equations. In this chapter we will study the general SQIR epidemiological model for COVID-19 with saturated incidence rate and using ABC-fractional operator. We will compare the results obtained for the fractional model with the classical model by predicting the model using Least Squares method and real data of Mexico from April 12 to October 28, 2020. Then the novelty of this paper is to use a different fractional order for each equation of the SQIR epidemic system with AB-fractional operators, this is fractional incommensurate SQIR model with time fractional derivative and saturated incidence rate.

The paper is organized as follows. In section 2, we present the Background of AtanganaBaleanu derivative and integral operators. The classical and fractional models are formulated in Section 3. The fixed point theorem is applied to prove existence and uniqueness results in section 4. The disease-free equilibrium and their local stability are given in section 5. In section 6 the Euler numerical method is used to determine the numerical solutions of the fractional proposed model. The model parameters are estimated according to real data of Mexico and the influence of fractional order has been shown in section 7. Finally, we gave some remarks and conclusions in section 8.

## 2. Background of Atangana-Baleanu operators

We first give some preliminary on the definitions of fractional derivatives will be used in this paper. [35].

*Let* (), *then the Atangana–Baleanu-Caputo (ABC for short) derivative is given as*
*where B is a function fulfills the condition B*(1) = *B*(0) = 1. *and E*_{α} *is the generalized Mittag–Leffler function given by*
*where* Γ *is Gamma function*.

The corresponding integral is given by the flowing definition.

*The AB-fractional integral associated to the Atangana-Baleanu fractional derivative is defined as*

When *α* → 0 in the above definition we find the initial function. Moreover, we find the Riemann integral if the fractional order *α* turns to 1.

## 3. The SQIR model formulation

We consider in this part the model of COVID-19. The population is divided into four compartments, the first is the susceptible class *S*(*t*), Infected class *I*(*t*), Quarantined class *Q*(*t*) and Recovered class *R*(*t*).

The transfer diagram for this model is described by Figure 1 and the classical version of this model formulated by the flowing system of four nonlinear ODEs :

with the initial conditions *S*(0) = *S*_{0} ⩾ 0, *Q*(0) = *Q*_{0} ⩾ 0, *I*(0) = *I*_{0} ⩾ 0 and *R*(0) = *R*_{0} ⩾ 0.

### 3.1. Fractional model

In this section, we change the classical derivative in (1) by ABC-fractional derivatives mentioned in Section 2.

Now remplacing the classical derivative in (1) by Atangana-Baleanu-Caputo derivative we obtain by moment closure, the following system of ABC-fractional ODEs
with the initial conditions *S*(0) = *S*_{0} ⩾ 0, *Q*(0) = *Q*_{0} ⩾ 0, *I*(0) = *I*_{0} ⩾ 0, *R*(0) = *R*_{0} ⩾ 0 and , *i* = 1, 2, 3, 4, is the ABC-fractional operator of order 0 *< α*_{i} *<* 1.

The total population *N* (*t*) = *S*(*t*) + *I*(*t*) + *R*(*t*) + *Q*(*t*). Then
In the absence of the disease, , this implies
The dynamics of the fractional COVID-19 model will be studied in a biologically feasible set denoted by

## 4. Existence and uniqueness of solutions of the ABC-fractional model

In this section, we prove the existence of the system of solutions by applying the fixedpoint theorem.

Let ℋ = (*C*([0, *T*]))^{4}, and *C*([0, *T*]) be a Banach space of all real-valued, continuous functions on the interval [0, *T*] with the norm
where ‖ · ‖ denote the supremum norm in *C*([0, *T*]).

We use the following notations for the right side of the COVID-19 models
by applying the definition of fractional integral to the Equation (1) and using the properties of fractional calculus we get
To prove the following theorems, we will assume that ‖*S*(*t*) ‖ ≤ *C*_{1}, ‖*Q*(*t*) ‖ ≤ *C*_{2}, *‖I*(*t*) ‖ ≤ *C*_{3} and *‖R*(*t*) *‖* ≤ *C*_{4}, where *C*_{i}, *i* = 1, …, 4, are some positive constants. Denote

*Assume that the following inequality is verified*
*then the kernels G*_{1}, *G*_{2}, *G*_{3} *and G*_{4} *fulfill Lipschitz conditions, furthermore are contraction mappings*.

Proof. Let *S*_{1} and *S*_{2} be two functions of susceptible idividuals, then we have
where *L*_{1} is defined in (4). Similar results for the kernels *G*_{2}, *G*_{3}, *G*_{4}, and *G*_{5} can be obtained using {*Q*_{1}, *Q*_{2}}, {*I*_{1}, *I*_{2}} and {*R*_{1}, *R*_{2}} respectively, as follows
where *L*_{2}, *L*_{3} and *L*_{4} are expressed in (4). Therefore, the Lipschitz conditions are fulfilled for *G*_{2}, *G*_{3}, and *G*_{4}. In addition, since 0 ≤ max{*L*_{1}, *L*_{2}, *L*_{3}, *L*_{4}} *<* 1, the kernels are contractions.

*Assume that the following inequality fulfill:*
*then the COVID-19 model with ABC-fractional operator (2) has a unique solution*.

Proof. The proof of this theorem is based on the fixed point method and using Picard operator. We consider the following notations
Keeping in mind the above notation, system (2) may be expressed as
where
Now using (3) and according to the above notation we get the following fractional integral equation
define
Now using the metric on C (𝕁, *B*_{1}, *B*_{2}, *B*_{3}, *B*_{4}) induced by the norm given as
We define Picard’s operator as
Now let’s try to prove that 𝒫 is a contraction.

Let 𝕫_{1}, *𝕫*_{2} ∈ C (𝕁, *B*_{1}, *B*_{2}, *B*_{3}, *B*_{4}) two functions, then we have
by using the triangular inequality and the Lipschitz condition we will find
If condition (6) fulfill then 𝒫 is a contraction on (𝕁, *B*_{1}, *B*_{2}, *B*_{3}, *B*_{4}). Then by Banach fixed point theorem 𝒫 has a unique fixed point. Hence, the fractional COVID-19 model (2) has a unique continuous solution.

## 5. Basic reproduction number, existence and stability of equilibria

In this section we determine the equilibria of the COVID-19 models, for that we solve :
then the disease-free equilibrium is given by
where and *I*^{0} = *R*^{0} = 0.

Now let’s determine the expression of basic reproduction number denoted by ℛ_{0}, by using the method presented in [36], we get
Furthermore, the system (2) has a unique endemic steady state
such as
where the value of *I** is the positive solution of the following algebraic equation
with
It remains to show that the solution of (17) is real and positive. We note
Let the discriminant of (17) be Δ, so that Δ = *B*^{2} −4*AC*.We conclude the following results for the endemic equilibrium.

*The system (2) has a unique endemic steady state whenever* ℛ_{0} *>* 1 *given by*

*where the expressions of S**, *Q***and R** *are given in (16) and* .

### 5.1. Local stability analysis for the fractional model

In this part we prove the local asymptotic stability of the equilibria of the commensurate model of (2), that is to say in the case where the equations have the same fractional order *α* ∈ (0, 1). Since *R*(*t*) is not present in first three equations, then we can write the system (2) as
To show the local asymptotic stability of Atangana-Baleanu-Caputo model, it suffices to show that all the eigenvalues *λ*_{1}, *λ*_{2} and *λ*_{3} of the jacobian matrix evaluated at the equilibrium point satisfy the condition
this result is mentioned in [37, Lemma 3.1.] and other references [38, 39].

*If* ℛ_{0} *<* 1, *then ℱ*_{0} = (*S*^{0}, *Q*^{0}, *I*^{0}) *is locally asymptotically and if ℛ*_{0} *>* 1, *then* F* = (*S**, *E**, *I**) *is locally asymptotically stable*.

Proof. The Jacobian matrix *J* (*ℱ*_{0}) for the fractional model given by (18), at the equilibrium *ℱ*_{0} is given by:
We find the eigenvalues by solving the following characteristic equation
We obtain the equation
The eigenvalues of matrix *J* (*ℱ*_{0}) are
If *ℛ*_{0} *<* 1 then *λ*_{1}, *λ*_{2} and *λ*_{3} are real negative roots. Then
Hence the disease free equilibrium of the fractional model (18) is locally asymptotically stable.

Now, we prove the local stability of the endemic equilibrium point. We consider the Jacobian matrix *J* (*ℱ**) for the fractional COVID-19 model given by (18), at the equilibrium *ℱ** is given by
The eigenvalues of this matrix can be obtained by solving the flowing characteristic equation
Therefore, we find the following characteristic equation:
where
Let us denote by
its discriminant, the following lemma completes the proof of theorem.

*Asume that* ℛ_{0} *>* 1 *and one of the following conditions is satisfied:*

Δ

*>*0,*c*_{1}*>*0,*c*_{2}*>*0*and c*_{1}*c*_{2}−*c*_{3}*>*0.,

*c*_{1}≥ 0,*c*_{2}≥ 0*and c*_{3}*>*0.

*Then the equilibrium ℱ** *of the fractional COVID-19 model with ABC derivative is locally asymptotically stable*

Proof. The detailed proof of this lemma is similar to [40, Proposition 1.].

## 6. Numerical method

To develop approximation schemes for the fractional model, we will use the Fractional version of Euler Method [41]. To applied this method, using same notation of section 4 and we consider the incommensurate system of ABC-fractional differential equations
We consider *t*_{0} = 0 *< t*_{1} *<* … *t*_{k} … *< t*_{N} = *T* an uniform time discretization of an interval [0, *T*] in *N* + 1 points. Let is the time step size. We denote as the numerical approximation of ℤ ^{i} (*t*_{k}).

Using AB-fractional integral, we have
This implies
At the point *t*_{n+1} = (*n* + 1)*h*, where *h* = *t*_{n+1} − *t*_{n}, we have
after using the fractional Euler method [41], we obtain
where the coefficients , *j* = 0,…,*k*, are given by
According to [41], the numerical scheme (27) for the fractional model (2) is conditionally stable. Now we will give a theorem for estimating the error between exact and numerical solution of the proposed model.

≤ {}

*Let ℤ*(*t*) *and* (ℤ_{k})_{k} *be respectively the exact and approximate solutions of the fractional incommensurate model (2). If* 0 ≤ max {*L*_{1}, *L*_{3}} *<* 1, *where L*_{1} *and L*_{3} *are given in (4), then we have*

*where C is an independent positive constant free of k and h*.

Proof. The proof is a direct conclusion of [41, Theorem 4].

## 7. Numerical simulation and calibration of the fractional model

### 7.1. Estimation of model parameters and best fit of fractional orders

The model calibration problem seeks to estimate the model parameters which, to some extent, to make the error between the observed values and the predicted values as small as possible, we use for that the Least Squares method to produce the values of the parameters and the fractional multi-order. Some model parameters can be known by these real values, on the other hand some parameters are unknown and need to be estimated, the same for the initial conditions of the COVID-19 model.

We consider the model solution *Z*(*t*) = (*S*(*t*), *Q*(*t*), *I*(*t*), *R*(*t*)), depending on the vector of parameters *ξ* = (*β, c, λ, θ, δ, γ, σ*, Λ), the initial conditions *Z*_{0} = (*S*_{0}, *Q*_{0}, *I*_{0}, *R*_{0}) and the fractional orders *α* = (*α*_{1}, *α*_{2}, *α*_{3}, *α*_{4}, *α*_{5}). The initial conditions *S*_{0}, *I*_{0} and *R*_{0} are known from data of COVID-19, but *Q*_{0} is unknown so it is necessary to estimate its value.

The vector of parameters belongs to a set given by
where *l*_{b} and *u*_{b} are of lower and upper vector values for the model parameters respectively.

Let the vector Θ ∈ ℝ^{N} of the observation data at given times *t*_{k}, *k* = 1, … *N* and a prediction vector Ψ(*ξ, Q*_{0}, *α*), the calibration aims at finding a vector such that
The aim is to minimize the objective function
subject to (*ξ, Q*_{0}, *α*) ∈ Θ and Equations (2) to obtain the best estimate of parameters and the vector of fractional orders α.

Now we consider the real data of COVID-19 for Mexico (that are found on the WHO web-page) to estimate the model parameters using the least squares method. In the first time we estimate the parameters for the classical model (1), second for the commensurate model which the four equations of the fractional model have the same order and finally the incommensurate model (2) when each equation has a different order.

The resulting fitted curves for reported infected and recovered cases of COVID-19 in Mexico is shown in Figure 2 for the classical model, in Figure 3 for the commensurate model and in Figure 4 for the incommensurate model, where the best fitted values of the parameters are depicted in Table 3 and the estimated values of the fractional orders of the two fractional models are given in Table 4. Furthermore, the initial conditions are given in the table 2.

From Table 3 we notice that our models give a good approximation of COVID-19 real data in Mexico compared to the classical model and more better than the case in which the orders of the four model equations have the same fractional order. The error of least square for the commensurate and incommensurate models are respectively **3.2931e+2** and **3.2901e+2**, this means that the incommensurate model with ABC-fractional derivative gives a better fit of the real data compared to the commensurate model with ABC-fractional derivative. In general the ABC-fractional derivative gives results that are closer to the real data compared to the classical model.

### 7.2. Sensitivity analysis of model parameters with respect to ℛ_{0}

When a virus spreads quickly it is necessary to look for methods to control the virus, that is, seek to increase or decrease some parameters which have an influence in the transmission of the virus. Sensitivity analysis is one way to determine the importance of each parameter for disease transmission.

Precisely, we use the definition of the sensitivity index which gives the significance of the variable associated with each parameter of the model to reduce the disease, the sensitivity index is defined as in [42].

*The sensitivity index of* ℛ_{0} *with respect to x, defined by*

The sign of each index makes it possible to know whether the parameter increases (positive sign) or decreases (negative sign) the value of ℛ_{0},

Using the expression of ℛ_{0} given in (15) to calculate the sensitivity index for ℛ_{0} with respect to the parameters of the model, we will find
Note that ℛ_{0} does not depend on *σ* so .

From Table 5, we can see that while which means that an increment in Λ, *c, β, λ* and *δ* will cause ℛ_{0} to increase, while a decrement in the recovery rate of infected γ or the quarantine rate of susceptible *θ* will cause ℛ_{0} to decrease.

### 7.3. Impact of fractional order

Based on numerical fractional Euler method derived in section (6), we present the graphical results of the proposed incommensurate fractional model (2) in Figures 5 and 6 to analyze the influence of fractional order. We have considered the initial values given in Table 2 and the estimated parameters given in Table 3.

The obtained figures show the importance of choosing the fractional order in the 2019nCOV model, as well as choosing a different order in each equation. From Figures 5-6, we see that any slight change in the fractional orders gives a change in the dynamics of the infected people and this indicates the effectiveness of this type of derivation. It can also be used to predict the number of people infected with COVID-19 better than the classical model and the fractional commensurate model (see Table 3).

## 8. Conclusion

In this chapter we have studied an incommensurate fractional model of COVID-19 with Atangana-Baleanu fractional derivative, we have proved the existence together with the uniqueness of the solution using the fixed-point theorem. We have determined the equilibrium points of the proposed model and the expression of base reproduction number *R*_{0} as well as the local stability of each equilibrium point. The numerical fractional version of Euler method is used to approximate the solutions of the fractional model. We simulated the models with five different fractional order vector and we compared the results obtained. The model efficiency is shown by comparing the error values between the COVID-19 data in Mexico and the values predicted by the three models studied. The best fit for COVID-19 data in Mexico is established when using incommensurate model with Atangana-Baleanu operator. The Least Squares method is used to determine the best fit to real data and estimate fractional order. Finally, the chapter gives an example of the effectiveness of using the Atangana-Baleanu operators in real problems and the importance of incommensurate propriety in the fractional model.

## Data Availability

None

## Footnotes

*Email addresses:*abdelouahed.allahamou{at}gmail.com (Abdelouahed Alla hamou), elhoussine.azroul{at}gmail.com (Elhoussine Azroul), hammouch.zakia{at}gmail.com (Zakia Hammouch), lamranii{at}gmail.com (Abdelilah Lamrani alaoui)