








key: cord-308115-bjyr6ehq
authors: Baba, Isa Abdullah; Nasidi, Bashir Ahmad
title: Fractional Order Model for the Role of Mild Cases in the Transmission of COVID-19
date: 2020-10-20
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.110374
sha: 
doc_id: 308115
cord_uid: bjyr6ehq

Most of the nations with deplorable health conditions lack rapid COVID-19diagnostic test due to limited testing kits and laboratories. The un-diagnosticmild cases (who show no critical sign and symptoms) play the role as a route that spread the infection unknowingly to healthy individuals. In this paper, we present a fractional order SIR model incorporating individual with mild cases as a compartment to become SMIR model. The existence of the solutions of the model is investigated by solving the fractional Gronwall's inequality using the Laplace transform approach. The equilibrium solutions (DFE & Endemic) are found to be locally asymptotically stable, and subsequently the basic reproduction number is obtained. Also the global stability analysis is carried out by constructing Lyapunov function. Lastly, numerical simulations that support analytic solution follow. It was also shown that when the rate of infection of the mild cases increases, there is equivalent increase in the overall population of infected individuals. Hence to curtail the spread of the disease there is need to take care of the Mild cases as well.

The outbreak of the novel strain of Corona viruses (COVID- 19) started late December in the Wuhan province in China [1] . It became a global pandemic causing the devastating impact in terms of morbidity, infections and fatality in addition to socio-economic disaster. The virus source which is yet to be identified is said to have genetic linkages with Severe Acute Respiratory Syndrome (SARS-COV) but less severe than Middle East Respiratory Syndrome (MERS-COV) [2] .

The virus is transmitted to healthy persons via eyes, mouth and nose when an infected person produced respiratory droplets of cough and sneeze or as a result of contact with contaminated surfaces. The average incubation period from catching the virus to time onset of major symptoms (like fever, cough and sneeze) is between 2 -14days [3] .

As the vaccine is not yet found, the control measures such as: social distancing, quarantine of suspected case, use of personal protective equipment (like face Mask, hand globes, gown), regular hand sanitation using antibacterial agents (like soaps, sanitizer) and imposing the lockdown curfew(when necessary) are the effective intervention that mitigate the transmission of the infection.

To execute these measures effectively, there is need to have an in depth study about the number of persons that each infected individual can infect, meanwhile a mathematical model describing the transmission dynamics of the disease should be established. In this regard,Zhao and Chen [4] developed a Susceptible, Un-quarantined infected, Quarantined infected, Confirmed infected (SUQC) model to characterize the dynamics of COVID-19 and explicitly parameterized the intervention effects of control measures. Similarly, Song et al. [5] developed a mathematical model based on the epidemiology of COVID-19, incorporating the isolation of healthy people, confirmed cases and close contacts.

Tahir et al. [6] developed a mathematical model (for MERS) inform of nonlinear system of differential equations, in which he considered a camel to be the source of infection that spread the virus to infective human population, then human to human transmission, then to clinic center then to care center. However, they constructed the Lyapunov candidate function to investigate the local and global stability analysis of the equilibriums solution and subsequently obtained the basic reproduction number or roughly, a key parameter describing transmission of the infection. Also,Chen et al. [7] developed a Bats-Hosts-Reservoir-People (BHRP) transmission network model for the potential transmission from the infection source (probably bats) to the human infection, which focuses on calculating .

To suit Korean outbreak, Sunhwa and Moran [8] established deterministic mathematical model (in form of SEIHR), in which they estimated the reproduction number and assessed the effect of preventive measures. Similarly, Lin et al. [9] modeled (based on SEIR) the outbreak in Wuhan with individual reaction and governmental action (holiday extension, city lockdown, hospitalization and quarantine) in which they estimated the preliminary magnitude of different effect of individual reaction and governmental action. Also Yang and Wang [10] proposed a mathematical model to investigate the current outbreak of the coronavirus disease (COVID- 19) in Wuhan, China. The model described the multiple transmission pathways in the infection dynamics, and emphasized the role of environmental reservoir in the transmission and spread of the disease. However, the model employed non-constant transmission rates which change with the epidemiological status and environmental conditions and which reflect the impact of the ongoing disease control measures.

Nonlocality is one of the main drivers of interest in fractional calculus applications. There are interesting phenomena that have what are called memory effects, meaning their state does not depend solely on time and position but also on previous state. Such system can be very difficult to model and analyze with classical differential equations, but nonlocality gives fractional derivative built-in ability to incorporate memory effects [11] .Fractional differential equations appear naturally in numerous fields of study including physics, polymer rheology, regular variation in thermodynamics, biophysics, blood flow phenomena, aerodynamics, electrodynamics of complex medium, viscoelasticity, capacitor theory, electrical circuits, electron-analytical chemistry, biology, control theory, and fitting of experimental data [12 -14, 23 -26] . Recently there are many studies on epidemiological disease modeling using fractional order differential equations [27 -33] .

The Rieman-Liouville fractional derivative is mostly used by mathematicians, but this approach is not suitable for real world physical problems since it requires the definition of fractional order initial conditions, which have no physically meaningful explanation yet. Caputo introduced an alternative definition, which has the advantage of defining integer order initial condition for fractional differential equations [15] .

In mostly poor and underdeveloped territories where there is no capacity of rapid diagnostic COVID-19 test due to insufficiency of testing kits, the Mild cases (who usually show no symptoms of the infection due to their strong and active immune system up to their recovery period) play a major role as a route that spread the disease to healthy individuals. Here we build our model by incorporating the population of mild individuals into the compartmental SIR model to become SMIR model in form of system of fractional order differential equations (FODE) in the Caputo sense. It should be noted that to our knowledge no model in literature considered the contribution of the mild cases of COVID -19 in the proliferation of the pandemic.

The paper is organized as follows: section 1is the introduction, section 2 is the preliminary definitions and theorems, the model formulation followed in section 3, section 4 is the stability analysis, section 5 gives numerical simulations and discussions and section 6 gives the conclusion. The notion of convergence of mittag-Leffler function is fully discussed in [15] .

Theorem 1 [19] :The equilibrium solutions of the system

is locally asymptotically stable if all the eigenvalues of the Jacobian matrix evaluated at the equilibrium points satisfy 

Despite the fact that almost 80% of the COVID -19cases are mild who recover naturally (due to stronger and active immune system that fight against the virus), they still play a role as a route of transmission of the infection [2]. The model was formulated based on the assumption that new born of human are recruited into susceptible class ( ) at the rate . The susceptible individual who had contact with an infected at rate can developed mild symptoms and move to mild class ( ). Based on [2], the mild patients with infectivity rate play the role as a routine that spread the infection, those with strong immunity recovered naturally at the rate , While some with critical illness became infected and moved to infectious class ( ) after the incubation period .

The infectious individual may then recovered ( ) or died at the rates and respectively. Figure 1 gives the schematic diagram describing the transmission dynamics of the disease. 

With the total population

The linearity of the Caputo operator yield

We apply the Laplace transform method to solve the Gronwall's inequality (6) with initial condition ( )

Linearity property of the Laplace transform gives (7) to partial fraction gives

Taking the inverse Laplace transform of (8)

Where ( ), ( ) are the series of Mittag-Leffler function (as in definition 4) which converges for any argument, hence we say that the solution to the model is bounded. Thus,

Consider the system (1) through (4) written as

]

Proof

Reference to Picard-Lindelof theorem [20] we establish the following theorem. 

Since and its closed set, then is complete metric space.

The continuous system (10) can be transformed to equivalent integral equations as; (12) is equivalent to Volterra integral equation that solves (10) .

Define an operator in

Now, we need to verify that (13) satisfies the hypothesis of contraction mapping principle.

First to show

Hence the operator maps onto itself.

Secondly, to show that is a contraction, we have

Since by hypothesis , then is a contraction and has a unique fixed point.. Thus, system (10) has unique solution.

To obtain equilibrium solution, we set the system to zero and solve simultaneously as follows;

Considering ( ) ( ) , in (11)- (14) we find the endemic equilibrium,

Consider system (1) then, we have the following Jacobian matrix Clearly, allthe eigenvalues have zero imaginary part ( ) Therefore,

Hence by theorem 1 above, the disease free equilibrium is locally asymptotically stable.

To obtain basic reproduction number which is a key parameter describing the number of secondary infections generated by a single infectious individual, we consider the eigen values above.

For implies that

This threshold quantity which if less than one disease free equilibrium will be stable and if greater than one it is unstable is what we termed as basic reproduction ratio . Hence we define

The endemic equilibrium points can be rewritten in the form of 

To derive the Lyapunov candidate function for fractional order as in [24] , consider the family of quadratic Lyapunov function

And define the Lyapunov candidate function as

Applying lemma 2 [23] above 

Hence by theorem 3 above, the disease free equilibrium is globally asymptotically stable.

Case 2: At the endemic (positive) equilibrium, (28) becomes

We had earlier established that the positive equilibrium is stable if . Now consider the relations

Therefore back substitution the above relations into (30) yields

where ( ( )) ( ) ( ) Hence by theorem 3 above, the endemic equilibrium is globally asymptotically stable.

In this chapter we carry out numerical examples to support the analytic results using parameter values in Table 1 . For the variables we use ( ) ( ) ( ) ( ) When the rate of infection of the mild cases increases, there is equivalent increase in the overall population of infected individuals. Hence to curtail the spread of the disease there is need to take care of the Mild cases as well.

In mostly poor and underdeveloped territories where there is no capacity of rapid diagnostic COVID-19 test due to insufficiency of testing kits, the Mild cases (who usually show no symptoms of the infection due to their strong and active immune system up to their recovery period) play a major role as a route that spread the disease to healthy individuals. Here we build our model by incorporating the population of mild individuals into the compartmental SIR model to become SMIR model in form of system of fractional order differential equations (FODE) in the Caputo sense.

The existence of the solutions of the model was shown by solving the fractional Gronwall's inequality using the Laplace transform approach. Two equilibrium solutions, disease free and endemic were obtained. Both local and global stability of the equilibria were shown to depend on the magnitude of basic reproduction ratio. Numerical simulations were carried out and dynamics of the populations were shown to vary for different values of It was also shown that when the rate of infection of the mild cases increases, there is equivalent increase in the overall population of infected individuals. Hence to curtail the spread of the disease there is need to take care of the Mild cases as well. 

World Health Organization (WHO). Novel Coronavirus-China

Modeling the epidemic dynamics and control of COVID-19 in China

The impact of isolation on the transmission of COVID-19 and estimation of potential second epidemic in China. Preprints (www.preprints.org)

Stability behavior of mathematical model of MERS Corona virus spread in population

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

Estimating the reproductive number and the outbreak size of novel coronavirus (COVID-19) using mathematical model republic of Korea

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

A mathematical for the novel Coronavirus epidemic in Wuhan

What is fractional Calculus. Cantor's Paradise

A survey on existence result for boundary value problem of nonlinear fractional differential equations

The Analysis of fractional differential equations: an application-oriented exposition using differential operators of Caputo type

stability and dynamics of a fractional order Leslie-Gower model

Fractional Differential Equations

Stability analysis of fractional differential system with Rieman-Liouville derivative

Differential Equations: An Introduction to Differential Equations, creative commons, 94105

Differential Equation of Fractional Order: Methods, Results and Problems

stability result for fractional differential equations with application to control processing

Stability analysis of Caputo fractional-order nonlinear systems revisited

Lyapunov function for fractional; order systems

Volterra type Lyapunov function for fractional order epidemic model

Analysis of Caputo fractional-order model for COVID-19 with lockdown

Analysis of meningitis model: A case study of northern Nigeria

Mathematical modeling for adsorption process of dye removal nonlinear equation using power law and exponentially decaying kernels

Modeling chickenpox disease with fractional derivatives: From Caputo to Atangana-Baleanu

A new study on the mathematical modelling of human liver with Caputo-Fabrizio fractional derivative

A mathematical model for COVID-19 transmission by using the Caputo fractional derivative

Epidemiological and clinical characteristics of the COVID-19 epidemic in Brazil

A fractional differential equation model for the COVID-19 transmission by using the Caputo-Fabrizio derivative

On the mathematical model of Rabies by using the fractional Caputo-Fabrizio derivative

A mathematical theoretical study of a particular system of Caputo-Fabrizio fractional differential equations for the Rubella disease model

Analysis of the model of HIV-1 infection of CD4^{+} T-cell with a new approach of fractional derivative

We write to declare our interest in publishing our work titled "Fractional Order Model for the Role of Mild Cases in the Transmission of COVID-19 " with the journal " Chaos, Solitons and Fractals".

