key: cord-0890488-c0zd7pig
authors: Abdo, Mohammed S.; Shah, Kamal; Wahash, Hanan A.; Panchal, Satish K.
title: On a Comprehensive Model of the Novel Coronavirus (COVID-19) Under Mittag-Leffler Derivative
date: 2020-05-08
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.109867
sha: 11e9ece57ff751122ccace09987dcd0d0a06f7a0
doc_id: 890488
cord_uid: c0zd7pig

Abstract The major purpose of the presented study is to analyze and find the solution for the model of nonlinear fractional differential equations (FDEs) describing the deadly and most parlous virus so-called coronavirus (COVID-19). The mathematical model depending of fourteen nonlinear FDEs is presented and the corresponding numerical results are studied by applying the fractional Adams Bashforth (AB) method. Moreover, a recently introduced fractional nonlocal operator known as Atangana-Baleanu (AB) is applied in order to realize more effectively. For the current results, the fixed point theorems of Krasnoselskii and Banach are hired to present the existence, uniqueness as well as stability of the model. For numerical simulations, the behavior of the approximate solution is presented in terms of graphs through various fractional orders. Finally, a brief discussion on conclusion about the simulation is given to describe how the transmission dynamics of infection take place in society.

One of the greatest missions given to humankind is to control the environment within which they live. However, the growth of the human population speedily is increasing in relation to different cultures and ways of life. Although Humanity has evolved many species of devices and equipment to obtain a beautiful life. But this excessive development sometimes leads to disasters in the environment. Especially, the use of nutrients, carriages, mobiles, cosmetics, electrified and petroleum equipment, etc, making the environment highly contaminated and infused with viruses.

Newly, the whole world hardship a new coronavirus comprehensive and it was named (COVID- 19) which was claimed to flowed first in Wuhan city, China. It has been considered that the origin of COVID-19 is the transportation from animal to human as numerous infected status claimed that they had been to a local fish and wild animal market in Wuhan on 28 November [1] . Soon, some investigators assured that transportation also occurs person to other [2] . This virus dates back to 1965 when Tyrrell and Bynoe were identified when they passage a virus so-called B814 [3] . This virus is found in the cultures of the organs of the human embryonic trachea acquired by the respiratory system of an adult [4] .

Because infectious diseases are a major threat to humans as well as the country's economy. A proper understanding of the dynamics of the disease plays a significant role in reducing infection in society.

Application of an appropriate strategy against the disease transportation is another challenge. The mathematical modeling tactic is one of the main tools in order to handle these challenges. Many disease models have been evolved in the existing literature which authorizes us to explore and dominance the prevalence of infectious diseases in a better style. Most of these models rely on classical differential equations. However, in the past few years, it is paying attention that fractional differential equations can be applied to model global phenomena with a greater grade of precision and its applications can be found in various fields for instance dynamic, biology, engineering, control theory, economics, finance and in epidemiology.

The Fractional calculus transacts with differentiation and integration involving fractional order, which is further outstanding and beneficial than the ordinary integer order in the explanation of the real-world problems, also in the modeling of real phenomena due to a characterization of the memory and hereditary properties [5, 6] , the integer-order derivative doesn't familiarize the dynamics among two various points. Various types of fractional order or nonlocal derivatives were proposed in the present literature to transact the reduction of a traditional derivative. For instance, based on power-law, Riemann-Liouville introduced the idea of fractional derivative. After that Caputo-Fabrizio in [7] have proposed a new fractional derivative utilizing the exponential kernel. This derivative has a few troubles related to the locality of the kernel. Newly, to overcome Caputo-Fabrizio's problem, Atangana and Baleanu (AB) in [8] have proposed a new modified version of a fractional derivative with the aid of popularized Mittag-Leffler function (MLF) as nonsingular kernel and nonlocal. Since the generalized MLF is used as the kernel and it's guaranteed no singularity, the AB fractional derivative supply a stellar description of memory [9, 10, 11] .

Because of mathematical models, we can know the rate of change in the COVID-19 how the disease can affect people at risk and quarantined. The area dedicated to the study of the biological model of infectious diseases is a warm area for recent research. Numerous studies on mathematical models were presented to the study of stability theory and the results of existence and improvement of biological models, see [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24] . For example, T. Chen, et al., in [12] presented a big mathematical model for simulating the phase-based transmissibility of a novel coronavirus. A. Hussain et al., in [13] showed the existence and uniqueness of results for the fractional model involving fractional derivative in Atangana-Baleanu sense using Schaefer's and Banach's fixed point theorems. Moreover, they applied the Shehu transform and Picard method to explore the iterative solutions and its stability for the proposed fractional model. A conceptual model for the COVID-19 disease, which effectively catches the time line of this virus outbreak was proposed by Lin et. al. in [14] . In [16] A mathematical modeling and dynamics of a novel (2019-nCoV) are formulated under a fractional model with considering the available infection cases within 7 days. A mathematical model considered by Shaikh, et al. [15] to study estimate of the effectiveness of prophylactic measures, prophesying future outbreaks and possibility control strategies of COVID-19.

Due to the success of this operator in modeling infectious diseases, and motivated by the above useful applications of some fractional operators in epidemic mathematical models, in this paper we are studying the dynamics of the novel coronavirus model suggested by Chen et.al., [12] in the form of the system of the nonlinear differential equations involving the AB Caputo fractional derivative:

where ABC D ∝ 0 + denotes the Atangana-Baleanu-Caputo fractional derivative of order ∝. The model is based on the following facts:

1. The bats were split into four closets: S B is a susceptible class of bats, E B is an exposed class of bats, I B is an infected class of bats, and R B is a removed class of bats. 2. The hosts were split into four classes: S H is a susceptible class of hosts, E H is an exposed class of hosts, I H is an infected class of hosts, and R H is a removed class of hosts. 3. The people were split into six classes: S ρ is a susceptible class of people, E ρ is an exposed class of people, I ρ is an infected class of people, A ρ is asymptomatic infected people, R ρ is a population of removed people due to death, and W is the population of the virus in reservoir various.

, and W 0 are the initial values corresponding to the three categories.

The parameters of the model (1) are described as follows: Λ B denote the birth rate of bats, m B is the death rate of bats, and β B is the transmission rate from I B to S B . The parameters 1 wB , 1 γB , 1 w ρ , and 1 γ ρ are the incubation period of bats, the infectious period of bats, the latent period of people, and the infectious period of asymptomatic infection of people, respectively. Λ H is Recruitment from total number of hosts and m H The death rate of hosts. The symbols β H and β BH represent the transmission rate from I H to S H and from I B to S H , respectively. Λ ρ and m ρ are the birth and death rate parameter of people, respectively. In the proposed model β ρ and κ are define the transmission rate from I p to S p and multiple of the transmissibility of A p to I p , respectively. β ω and β ρ means the transmission rate from W to S p and from I p to S p , respectively. δ ρ is the proportion of asymptomatic infection rate of people, γ ρ is the infectious period of symptomatic infection of people and a is the retail purchases rate of the hosts in the market. The shedding coefficients from I ρ to W and the shedding coefficients from A ρ to W are denoted by µ ρ and µ ρ , respectively. Finally, 1 ε is the lifetime of the virus in W and N H is the total number of hosts.

The major aim of the paper is to demonstrate the existence, uniqueness and Ulam stability of solution for the model (1)-(2) by using some fixed point techniques. Moreover, the numerical simulations via the fractional version of Adams Bashfully technique to approximate the ABC fractional operator are performed, which it is shown that the model displays rich dynamical behaviors through graphical representation of numerical solutions. This paper is coordinated as follows: Section 1 transacts with the introduction which contains a survey of the literature. Section 2 consists of some foundation preliminaries related the fractional calculus and nonlinear analysis. The existence and Ulam stability results on a proposed model are obtained in Sections 3, 4. The numerical solution and numerical simulations of the model at hand are presented in Section 5. For the numerical simulation, we use a powerful two step numerical tool called fractional AB method. The concerned numerical method is more powerful than usual Euler method as well Taylor method. Because the mentioned method is faster convergent and stable as compared to Taylor and Euler method which are slowly convergent. For detail about this method, see [32, 33, 34, 35] .

For short, we will use the following notations

For the next analysis, we define Banach space for onward analysis. Let 0 ≤ t ≤ T < ∞, we define the Banach space by using J = [0, T ] as = C(J, R 14 ) under the supremum norm given by

Definition 2.1. [8] Let ∝∈ (0, 1] and σ ∈ H 1 (0, T ). Then the left-sided ABC fractional derivative with the lower limit zero of order ∝ for a function σ is defined by

where N (∝) is the normalization function which is defined as N (∝) = ∝ 2−∝ , 0 <∝≤ 1 and satisfies the result

here Re (∝) > 0 and Γ (.) is a gamma function.

Then the left-sided ABC fractional integral with the lower limit zero of order ∝ for a function σ is defined by

The Laplace transform of ABC fractional derivative of a function σ(t) is given by

The solution of the proposed problem for ∝∈ (0, 1]

is given by

be a Banach space. The operator Π : → is a Lipschitzian if there exists a constant κ > 0 such that

here κ is the Lipschitz constant for Π. If κ < 1 we say that Π is a contraction. Theorem 2.6. [25] (Let M be a closed, convex, non-empty subset of a Banach space . Suppose that E and F map M into and that

Then, there exists Ψ ∈ M such that EΨ + FΨ = Ψ.

Now, we debate the existence and uniqueness results of the model (1)-(2) by utilizing the fixed-point technique. Let us reformulated model (1) in the subsequent appropriate form

We take our system as by using (1)

where

Here the symbol A T denotes the transpose operation. Utilizing Lemma 2.3.1, the model (6) can be turned to the fractional integral equation in the sense of AB fractional integral as follows

Expressing some growth condition and Lipschitzian assumption for existence uniqueness as:

(Hypothesis 1) There exists two constants µ H , η H such that

(Hypothesis 2) There exists constants L H > 0 such that

for each Ψ ∈ and t ∈ [0, T ].

Theorem 3.1. Let the hypotheses 1,2 hold. The Integral equation (8) which equivalent of the considered

We define the operators Π 1 and Π 2 as

where = Π 1 + Π 2 . Now, we offer the proof in several steps as:

Step1

Indeed, From (Hypothesis 1) and equations (9),(10), we obtain

This confirms that

Step1: we show that Π 1 is contraction.

Let Ψ, Ψ * ∈ B λ . By (Hypothesis 2), we have

Step 3: We show that Π 2 is relatively compact.

To prove that, we show that Π 2 is continuous, uniform bounded, and equicontinuous. Since Ψ(t) is continuous, then Π 2 Ψ(t) is continuous. Next, let Ψ ∈ B λ , we have

Hence Π 2 is uniformly bounded on B λ . Finally, we show that Π 2 equicontinuous. Let Ψ ∈ B λ and t 1 , t 2 ∈ [0, T ] such that t 1 < t 2 . Then

As t 1 → t 2 the right-hand side of the above inequality tends to zero. Consequently, by Arzelá-Ascoli theorem, Π 2 is relatively compact and so completely continuous. Thus by Theorem 2.6, the integral equation (8) has at least one solution and consequently the model under consideration has at least one solution. 

Proof. Let the operator Π : → defined by

Let Ψ and Ψ * in and t ∈ [0, T ]. Then

Due to (11) , Π is contraction. Therefore the integral equation (8) has unique solution. Hence our model (1)-(2) has unique solution.

The notion of Ulam stability was initiated by Ulam [26, 27] . Then aforesaid stability has been scrutinized for classical fractional derivatives in many of the research articles, we refer to some of them like [28, 29, 30, 31] . Additionally, since stability is a prerequisite in respect of approximate solution, so we endeavor on Ulam type stability for the model (1) via using nonlinear functional analysis. 

then there exists Ψ ∈ satisfying the model (1) with the following initial condition

such that

= 1, 2, 3, ..., 14, = max( 1 , 2 , 3 , ...., 14 ) T , λ = max(λ 1 , λ 2 , λ 3 , ...., λ 14 ) T , Remark 1. Consider a small perturbation g ∈ C[0, T ] such that g(0) = 0 comply with the following properties.

1. |g(t)| ≤ , for t ∈ [0, T ] and 1 > 0. 2. For t ∈ [0, T ] we have the following model 

satisfies the given relation .

Proof. Thanks to Remark 1 (2), and Lemma 2.3.1, the solution of (15) is given by

Also, we have

It follows from Remark 1(1) that 

That is

Due to (16) , Ψ 0 = Ψ 0 , the equation (17) becomes

Thus by Hypothesis 1 and Lemma 4.1.1, we obtain

Due to Θ 3 < 1. For λ = 2κ 1−Θ3 , we get Ψ − Ψ ≤ λ . Hence the model (1)-(2) is U-H stable.

In this part, we give approximation solutions of the ABC fractional model (1)- (2) . Then the numerical simulations are acquired via the suggested scheme. To this aim, we employ the modified fractional version for Adams Bashforth technique to approximate the fractional integral in the sense AB.

Using the initial conditions and fractional integral operator, we convert model (1) into the integral equations

To procure an iterative scheme, we go ahead with the first equation of the model (19) as follows:

Set t = t r+1 , for r = 0, 1, 2, ..., it follows that (20) Now, we approximate the function P 1 (ζ, S B ) on the interval [t , t +1 ] through the interpolation polynomial as follows

Now, we compute the following integrals I −1,∝ and I ,∝ as follws

and

Put t = h, we get

and

Substituting (22) and (23) into (21), we get

Similarly

Now to present the numerical simulations of the ABC fractional model (1)-(2), we apply the iterative solution contained in (24)-(37). Taking the time (Days) range up to 10 unit. The numerical amounts of the parameters applied in the simulations are specified in Table 1 The birth rate parameter of bats 0.0300 m B

The death rate of bats 0.011 β B

The transmission rate from I B to S B 0.01166

The incubation period of bats 0.00073 The total number of hosts 10 millions In Figures 1 − 14 , we have given a global dynamics of each compartment in the considered model (1) by using the numerical values in Table 1 against different fractional order. The growth and decay of various compartment is different that some compartments are growing with different rate due to fractional order and similar behavior may be observed for decay of some compartments. Lower the order faster is the growing rate and vice versa and similarly the decay also is different at different fractional order. Therefore fractional calculus can help in understanding the transmission dynamics of novel coronavirus-19 disease. Also the concerned numerical scheme can be used as a powerful tools to perform numerical simulation of such complicated model. In Figure 1 as susceptible papulation is decreasing as they are exposing so the density of expose class is going on growing as in Figure 2 . On the other hand the infected bats density is increasing as presented in Figure 3 , while the number of removed bats are also increasing in Figure 4 until stable in the market. On the other hand as the density of susceptible host is decreasing in Figure 5 because they are exposing to infection. Therefore the density of exposed host is growing as in Figure 6 as a results the host catching more infecting and hence the density of host infected is also increasing. This process may be observed from Figure 7 . As a results of infection the removed (either death or get ride from infection) will increasing as in Figure 8 . In same way the density of susceptible people is decreasing in Figure 9 because they are exposing to infection so the density of exposed people will decrease (see Figure 10 ), because the people rapidly catching infection, hence the density of symptomatic infected people will grow as in Figure 11 . The similar behavior one can observed for density of asymptomatic and removed (either death or recovered from infection) in Figure 12 and 13 respectively. Since as the density of infected bats is increasing, also the infection in host papulation is growing. As a result the papulation of virus in reservoir will grow further as presented in Figure 14 .

This manuscript has been devoted to comprehensively investigate a mathematical model for calculating the transmissibility of novel Coronavirus (COVID-19) disease by using nonsingular fractional order derivative. The existence and uniqueness of the considered model has been guaranteed by applying Krasnoselskii and Banach fixed point theorems. Also some stability results of Ulam type have been constructed. With the help of fractional Adams Bashforth method method we have simulated the results corresponding to various fractional orders. The obtained results play important role in developing the theory of fractional analytical dynamic for the current pandemic due to Coronavirus -19 which has badly affected the whole world. From the simulation we noted that the decrease in susceptibility is faster at lower fractional order of the derivative and in same line the increase in infections is also rapid with smaller order. The presented results may be helpful in understanding the present pandemic more comprehensively and can help in taking precautionary measure to reduce the infection to minimum.

Outbreak of Pneumonia of Unknown Etiology in Wuhan China: the Mystery and the Miracle

Homologous recombination within the spike glycoprotein of the newly identified coronavirus may boost cross-species transmission from snake to human

Cultivation of viruses from a high proportion ofpatients with colds

History and recent advances in coronavirus discovery

Theory and Applications of Fractional Differential Equations

A new definition of fractional derivative without singular kernel

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

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

Fractional derivatives with no-index law property: application to chaos and statistics

Decolonisation of fractional calculus rules: breaking commutativity and associativity to capture more natural phenomena

mathematical model for simulating the phase-based transmissibility of a novel coronavirus

On a nonlinear fractional-order model of novel coronavirus (nCoV-2019) under AB-fractional derivative

A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action

Model of COVID-19 Using Fractional Derivative

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

SIR model for Dengue disease with effect of Dengue vaccination

Semi-analytical study of Pine Wilt Disease model with convex rate under Caputo-Fabrizio fractional order derivative

Theory of Modeling and Simulation

On a two-dimensional magnetohydrodynamic problem: modelling and analysis

Numerical solution of SIR model of Dengue fever

New numerical method for ordinary differential equations: Newton polynomial

New dynamical behaviour of the coronavirus (COVID-19) infection system with nonlocal operator from reservoirs to people

The Existence and Analytical Study of Fractional order Dynamics of COVID-19 Under Attangana-Baleanu-Caputo Derivative

Basic theory of fractional differential equations

Modern Mathematics

A Collection of Mathematical Problems

Investigation of Ulam stability results of a coupled system of nonlinear implicit fractional differential equations

On Ulam's stability for a coupled systems of nonlinear implicit fractional differential equations

Ulam stability to a toppled systems of nonlinear implicit fractional order boundary value problem

Existence and Ulam-Hyers stability for Caputo conformable differential equations with four-point integral conditions

Numerical analysis for the fractional diffusion and fractional Buckmaster equation by the two-step Laplace Adam-Bashforth method

Stability analysis for fractional-order partial differential equations by means of space spectral time Adams-Bashforth Moulton method

A Generalization of the Adam-Bashforth Method

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

Numerical approximation of Riemann-Liouville definition of fractional derivative: from Riemann-Liouville to Atangana-Baleanu

Prepared the draft Kamal Shah: Developed the concept and methodology Hanan A. Wahash: Performed the analysis. Satish K. Panchal: Polished the final draft

We are thankful to the reviewers for their constructive suggestions which improved this paper very well.