Analysis of the Spread of Covid-19 via Atangana-Baleanu Fractional Derivatives

We study the spread of the epidemic via Atangana-Baleanu Fractional Derivatives. We present the mathematical analysis and formulation of a fractional model for the epidemic. The existence and uniqueness of the solution for the proposed model are proved. The study also investigates the existence of disease-free equilibrium and analyzes its stability properties. To validate the theoretical results, we provide a numerical scheme for the fractional model and present various simulation results. These results can serve as a valuable resource in developing strategies to mitigate the spread of the epidemic.


Introduction
The coronavirus epidemic emerged in the Chinese city of Wuhan in December 2019, and it rapidly disseminated to other parts of the world.This virus is primarily transmitted through respiratory infection in humans [2].The disease is characterized by symptoms such as difficulty breathing and fever.The incubation period of the virus ranges from 2 to 14 days before symptoms become evident [2,3,13].In February 2020, Algeria and other countries reported their first confirmed case of the virus in an Italian citizen who had recently arrived.Authorities promptly conducted scans on individuals showing suspicious symptoms or signs of infection.However, despite the daily checks carried out by health authorities, there remains a significant gap between the day of infection and the day of diagnosis.This gap has severe implications for the spread of the disease, as it hinders the accurate assessment of the number of affected individuals and the severity of the disease.Consequently, obtaining precise figures regarding the number of infected and recovering individuals becomes challenging.
In order to study the clinical development and spread of the epidemic, researchers formulate mathematical models that help assess the severity of the disease and determine necessary interventions to mitigate its impact, such as social isolation, quarantine measures, and travel restrictions [17].One widely used model for studying disease dynamics is the classical susceptible-infectedrecovered (SIR) model [14].These models provide insights into the dynamics of disease spread, and the timing of outbreaks and declines, and aid in determining when it is appropriate for individuals to return to work.Typically, these models are formulated as differential or stochastic equation sets.The initial results introducing these models date back to the following references [14].Our model encompasses seven distinct compartments, each representing a different stage or condition related to the epidemic.These compartments include: susceptible individuals, exposed individuals, recovered individuals, quarantined individuals, recovered-exposed individuals and death individuals.By considering these seven compartments, our model aims to provide a comprehensive representation of the different stages and conditions related to the epidemic.In our research, we propose a reformulation of the model presented in reference [13] by incorporating a fractional derivative of the Atangana-Baleanu type.The motivation behind utilizing the Atangana-Baleanu derivative (ABD) in our model is due to its distinct properties, such as a nonlocal and nonsingular kernel.These properties are crucial for accurately capturing the crossover behavior observed in the dynamics of the epidemic.Other fractional derivative operators, such as Riemann-Liouville, Caputo, and Caputo-Fabrizio, do not possess these specific properties, which may or may not adequately describe the future dynamics of the coronavirus.This motivates our choice of the Atangana-Baleanu derivative in our reformulated model [3].
Fractional differential equations have proven to be highly valuable in modeling various biological systems.This can be observed in several research papers that have explored their applications.Some notable papers related to the use of fractional calculus in modeling biological systems include [3,7,11,12,15].These studies have highlighted the effectiveness and relevance of fractional differential equations in capturing the complex dynamics and behaviors exhibited by biological phenomena.Additionally, there are specific papers that focus on approximating the Atangana-Baleanu derivative [4,9].These works further demonstrate the significance of fractional calculus in accurately describing and understanding the dynamics of biological systems.
This paper presents a comprehensive analysis of the epidemic by introducing a mathematical model that captures the dynamics of disease spread, as described in reference [13,18].Subsequently, we propose a fractional model of the epidemic and investigate the existence and uniqueness of solutions for our system.Additionally, we examine the existence of a disease-free equilibrium and analyze its stability properties.
The structure of the article is as follows: In Section 2, we provide definitions and explanations pertaining to fractional calculus.We then proceed to present the mathematical model of the coronavirus, followed by its formulation using fractional derivatives.In Section 3, we conduct a detailed analysis of the fractional model for the coronavirus, focusing specifically on the existence and uniqueness of solutions.We delve into the mathematical properties and characteristics of the model in order to gain deeper insights into its behavior.Moving on to Section 4, we explore the equilibrium points of the model and examine their stability properties.This analysis allows us to understand the long-term behavior of the epidemic and assess the potential for disease control and mitigation strategies.In Section 5, we present a numerical scheme for the fractional model, providing a practical approach to simulate and study its dynamics.We include illustrative results from numerical simulations conducted with various values of the fractional order paramete.These simulations help us visualize and analyze the impact of different parameters on the spread and control of the epidemic.
Finally, in the conclusion, we summarize our findings and offer recommendations based on the insights gained from the study.These recommendations aim to guide future research and inform potential strategies for managing and combating an epidemic.

Preliminaries
In this section, we will cover some of the most important basic concepts and definitions of fractional calculus and epidemic models.

Model description of corona-virus with Atangana-Baleanu fractional derivative
We analyze the corona-virus model discussed by Ben Fradj and Cherif [13] .Let us consider the total population at time t as N (t).We divide the population N (t) into seven distinct subgroups: susceptible individuals S(t), exposed individuals E 1 (t), infected individuals I(t), recovered-exposed individuals E 2 (t), recovered individuals R(t), deceased individuals D(t), and quarantined individuals [13].
The relationships and interactions between these subpopulations can be visually represented through the following compartmental diagram: The model is developed based on a system of nonlinear differential equations, which can be expressed as follows: with the initial conditions where λ 1 is the contact rate between S and E 1 , λ 2 is infection rate, λ 3 is the rate at which susceptible people enter in quarantine, α 1 is the inverse of the average latent time, α 2 is the recovered rate of E 1 , β 1 is the recovered rate, β 2 is the death rate and the total constant population N is give by To incorporate the aforementioned effect into the mathematical representation, we modify the system by replacing the traditional time derivative with the Atangana-Baleanu ordered derivative [6].The resulting formulation is as follows: where ABD 0 D γ t the Antangana-Baleaneau fractional derivative of order γ ∈ (0, 1], the model variables in (1) are non negative, and has appropriate initial conditions.

Existence and uniqueness of solutions via Atangana-Baleanu fractional derivative
The initial value problem (2) can be written as the following matrix from : In order to establish the nonnegative of solutions with the initial conditions, we need also the following lemmas.

Bounded and non-negative solutions
For the proof of the theorem about non-negative solutions and bounded we shall need the following lemma where x i (t) ∈ C[0, T ], i = 1, . . ., 7. The norm of X(t) ∈ C[0, T ] is given by : where T > 0 and ϵ > 0.
Let f : E → R be a function.Then, the fractional differential equation ABD a . ., m − 1, has a solution, and this solution is unique if the following two conditions hold simultaneously: (1) f is bounded and continuous on E.
(2) f satisfies the Lipschitz condition with respect to the second variable, i.e., where L > 0 is a constant independent of t, X and Y .
Theorem 1.The initial value problem (2) possesses a unique solution.Moreover, this solution remains nonnegative and bounded.
Proof.We will divide the proof into three steps: Step 1 Firstly, we will prove that the solution is nonnegative.From equation (2), we have: Based on Lemma (1) and Remark (2), we can conclude that the solution of equation ( 2) is nonnegative.
Step 2 Secondly we prove that solution is bouned.From (2) by using the first equation and adding the first two equations, we get : , by integration of (ABD) we get : By adding the first three equations of (2), and the same procedure we get : with M (t) = S(t) + E 1 (t) + I(t) and ξ 2 = min(λ 3 , α 2 , β 1 + β 2 ), by following the same procedure, we obtain: Hence, we have: Therefore, We can follow the same approach for the remaining variables E 2 , R, D, and Q, and obtain: Step 3 We define the matrix norm as |||.|||,where |||M ||| = ρ(M ) and ρ(M ) represents the largest eigenvalue of the matrix M .Let us consider the bounded and continuous function Assume that X(t) and Y (t) are two distinct solutions of the initial value problem (2), such that: By applying Lemma (2), and since F (X(t)) is continuous and satisfying Lipschitz condition, the initial value problem (2) has a unique solution.

The stability result of the equilibrium point
To determine the equilibrium points of the system (2), we set the Atangana-Baleanu arbitrarily ordered derivatives to zero, i.e., From this, we can deduce that the equilibrium point is given by E * = (0, 0, 0, 0, 0, 0, 0) T .Now, we proceed to evaluate the equilibrium point for the system (2).The Jacobian matrix of the system (2) is determined as follows: Substituting by the equilibrium point E * in the above Jacobian matrix, we get : The Jacobian matrix obtained from the system (2) has eigenvalues equal to zero.Therefore, the equilibrium point E * is determined to be unstable.

Numerical results and discussions
Upon applying the fractional integral operator ABI I to equation ( 2) and utilizing the property (1), we obtain the following numerical results for our fractional model with the ABD derivative: The kernel functions for the fractional epidemic model ( 2) are defined as follows: To solve the fractional epidemic model (2), we employ the numerical procedure outlined in reference [3].The model can be represented as follows: Discretizing the above equation yields the following system: By utilizing the general formula for Adams' method [16,2], we derive the subsequent numerical approach for the Corona virus model 2:

Data fitting and numerical simulations
We present a set of numerical findings.Firstly, we draw upon data from the National Office of Statistics (NOS) in Algeria to provide key demographic insights for the year 2019.The total population during that period stood at 43,400,000, and the average life expectancy was 76.79 years.Additionally, based on the available information, the estimated mortality rate is µ = 0.0048.Subsequently, we analyze data collected in Algeria from March 24, 2020, to June 30, 2020, and deduce the following parameter values: λ 1 = 0.8, λ 2 = 0.02, λ 3 = 0.155, α 1 = 0.01, α 2 = 0.2, β 1 = 0.023, and β 2 = 0.0041.
For the simulation, we have chosen the following initial conditions: The total population at time t = 0 is N (0) = 43, 400, 000, which can be expressed as To determine the number of susceptible individuals at t = 0, we use the formula Moreover, we consider the initial conditions for the infected and exposed populations as follows: I(0) = 70 and E 1 (0) = 200.This observation highlights the crucial role of the exposed population in driving the spread of the infection.As more individuals are exposed to the virus, there is a higher likelihood of transmission to susceptible individuals, leading to a rise in the number of infections.Understanding and monitoring the dynamics of the exposed population are crucial for assessing the potential risk and impact of an infectious disease outbreak.In Figure 4, the plotted curves illustrate the evolution of the number of infected individuals over time.It is evident from the graph that the number of infected people  The dynamics of infection can be influenced by altering the values of two key parameters: the fractional order γ and the quarantine rate λ 3 .By varying these parameters, different patterns in the infection dynamics can be observed.Adjusting γ affects the rate at which infected individuals recover, while modifying λ 3 influences the effectiveness of the quarantine measures in limiting the spread of the infection.
Understanding the impact of these parameters is essential for devising effective strategies to control and manage the spread of the infection.The simulation results provide valuable insights into the behavior of the infected population and the potential outcomes under different parameter settings.
In Figure 7, the graph showcases the evolution of the number of deceased individuals over time.The number of deaths is significantly influenced by changes in the number of infected individuals, as expected.As the infection spreads and the number of infected individuals rises, the number of deaths also increases due to the severity of the disease.
Moreover, the values assigned to the two parameters, the fractional order γ and the quarantine rate λ 3 , have a notable impact on the number of deceased individuals.Altering these parameters can lead to different outcomes in terms of the mortality rate and the overall impact of the infection on the population.
Understanding the relationship between the number of infected individuals and the resulting number of deaths, as well as the role of parameter values, is crucial for predicting and mitigating the consequences of an infectious disease outbreak.This information can aid policymakers and public health officials in making informed decisions to implement effective measures for reducing the impact of the disease on the population.
then the function f is non-decreasing, and if ABD D γ f ≤ 0 for all t ∈ (0, b), then the function f is non-increasing.Definition 4.[11] Let C[0, T ] be the class of continuous functions x(t) defined on the interval [0, T ], satisfying the condition sup|x(t)| < ∞ for t ∈ [0, T ].Additionally, let C([0, T ]) be the class of continuous column vectors X(t) defined as: