key: cord-0924064-dp39eshj authors: Paul, Subrata; Mahata, Animesh; Mukherjee, Supriya; Roy, Banamali title: Dynamics of SIQR epidemic model with fractional order derivative date: 2021-12-16 journal: Partial Differential Equations in Applied Mathematics DOI: 10.1016/j.padiff.2021.100216 sha: f65689d070dcfefa2446bd40238fbc9ebc47c239 doc_id: 924064 cord_uid: dp39eshj The dynamics of COVID-19 (Coronavirus Disease-2019) transmission are described using a fractional order SIQR model. The stability analysis of the model is performed. To obtain semi-analytic solutions to the model, the Iterative Laplace Transform Method [ILTM] is implemented. Real-time data from COVID-19 cases in India and Brazil is employed to estimate the parameters of the fractional order SIQR model. Numerical solutions obtained using Adam-Bashforth-Moulton predictor–corrector technique is compared with those obtained by ILTM. It is observed that the fractional order of the derivatives is more effective in studying the dynamics of the spread of COVID-19 in comparison to integral order of the SIQR model. The first cases of corona virus infection in human species with symptoms similar to common cold [1] were reported in 1965 by Tyrrell and Bynoe. However, after numerous mutations, the virus has ultimately proved to be catastrophic in Wuhan, Hubei Province, China, in 2019. Henceforth, the virus was named SARS-CoV-2 and the disease was declared to be COVID-19. The devastation spread quickly throughout multiple nations, prompting the WHO to identify a pandemic on March 11, 2020 [2] . No sooner had the futile efforts of the doctors and the administrative heads of states succumbed to the wrath of the virus, researchers from all areas of science sprang to activity. Along with medical researches going on in full swing, the study of the transmission of the disease, effect of preventive measures, prediction of future outbreaks and potential control strategies were also investigated extensively. The simple epidemic model by Kermack and Mckendrick in 1927 [3] is among the initial works on mathematical modelling of infectious diseases which later on proved to be the building block of further research [4] [5] [6] . In order to study such mathematical models, the entire population is primarily divided into subclasses, namely, susceptible individuals S(t), infected individuals I(t), the quarantined individuals Q(t) and the recovered individuals R(t). Integral order ordinary differential equations are common among epidemiological modelling of biological systems [7] [8] . However, there are certain reservations regarding the integral order of the differential equations in such models. Such restrictions are taken into account by a comparatively new and emerging area of mathematical calculus, namely the Fractional Order Differential Equations [9] [10] [11] [12] [13] . In fact, fractional-order differential equations (FODEs) and their applications have been intensively A mathematical modelling of infectious disease [26] is critical for better understanding the transmission patterns of the disease and evaluating control strategies. It serves as motivation for mathematical and biological experts to investigate and evaluate the dynamical systems that regulate such diseases in order to anticipate their spread and control in the long term. The four epidemiological divisions that assess the susceptible part of individuals, the infected, the quarantine, and the recovered, have characterized basic definitions of infectious diseases. One of the several variations of the conventional SIR model is the SIQR model. It has been observed that quarantining the infected individuals is a better measure to control the spread of the disease. Over these years, several types of fractional calculus, for example, Riemann-Liouville , Caputo, Caputo-Fabrizio, Katugampola, Atangana-Baleanu, Hadamard etc. have been introduced to study the dynamics of the epidemic models, each displaying certain advantages and disadvantages. Zhang et al. [27] developed the unique asymptotic stability criteria for fractionalorder gene regulation systems with time delay. Wu et al. [28] investigated global uniform asymptotical stability for fractional-order gene regulatory networks with time-varying delays and structured uncertainties. Based on the vector Lyapunov function, Zhang et al. [29] presented a unique stability condition for fractional-order composite systems with time delay. We have used the Caputo-Fabrizio operator because of its possession of a nonlocal and nonsingular exponential kernel and is found to be best suited to analyze the dynamics of COVID-19 [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] . The dynamics of COVID-19 (Coronavirus Disease-2019) transmission are described using a fractional order SIQR model. We investigate the disease's impact using an appropriate mathematical model [SIQR model] in context of the Caputo-Fabrizio fractional differential equation, motivated by early research. The present research purpose is to:  To study local stability as well as global stability of the model system.  Determination of the Basic Reproduction number and Equilibrium points.  The Iterative Laplace Transform Method has been used to generate semi-analytic solutions to the model system.  To confirm the Iterative Laplace Transform Method's predictions, the Adam-Bashforth scheme is applied to the SIQR model.  To confirm the findings and prevent the growth of COVID-19, a numerical simulation was used. The following is an overview of the article's structure:Section 2 discusses some fundamental concepts and findings for the fractional operator and the Laplace transform. Mathematical model with fractional order derivative is discussed in Section 3. Section 4 is focused to a description of the Model's stability analysis and stability criterion. In Section 5, we perform Adam-Bashforth-Moulton predictor-corrector scheme for the SIQR model. In Section 6, numerical simulation and discussion are presented via MATLAB. Section 7 includes the validation of the model with real time data. Finally, Section 8 constitutes the conclusion of the paper. The definitions of fractional differential and integral operators are introduced in this part. . (2.5) In particular, we have . (2.7) Let us consider the case of a Susceptible -Infected -Quarantine -Removed (SIQR) epidemic. The SQIR disease transmission model is derived assuming several strong assumptions [41] . The population (N) is divided into four classes: the susceptible individuals (S), the infected individuals (I), the quarantine individuals (Q) and the recovered individuals (R) at any time t ≥ 0. (3.1) Individuals who are Susceptible in this model are those who are at risk of becoming infected. Quarantine is defined as an infected person who exhibits signs of the disease and is isolated. Quarantined individuals who recover from the disease are considered as Recovered. Table 1 The parameters of the model and their descriptions Interpretations Natural birth rate Contact rate between S and I Infection rate of Q class Recovery rate Natural death rate The differential equations governing the SIQR model are given as: The model (3.2) having fractional order derivatives with the Caputo-Fabriozio operator of order 0 < ≤ 1have been proposed as follows: Proof From the model (3.4), Therefore, the model (3.4) is bounded by . Thus S, I and Q are positive functions and is positively consistent of the system (3.4). The basic reproduction number, indicated by ℜ 0 , is defined as the predicted number of secondary cases created by a single Susceptible individual's infection. The disease-free equilibrium point is locally asymptotically stable when ℜ 0 < 1 and unstable when ℜ 0 > 1. Since the considered model has disease-free equilibrium at ( , 0, 0, 0). So, an effective strategy should be developed, when a pandemic breaks out so that ℜ 0 reduces to less than 1 as soon as possible. Using next generation matrix method [42] [43] , the leading eigen value of the matrix −1 may be used to calculate the reproduction number (ℜ 0 ) where, Therefore, the reproduction number (ℜ 0 ) = ( + ) . (3.10) Consider the model (3.4) along with initial conditions (3.3) . Applying the Laplace transform to both sides of model (3.4), we derive Rearranging, we get Using inverse Laplace transform of the equations (3.12), we get, Using initial conditions, we get the recursive equations given by, 14) To obtain the equilibrium points of the system For further analysis, the Jacobian matrix of the system (3.4) at any equilibrium point ( , , , ) is given by The disease-free equilibrium of the system (3.4) is locally stable if ℜ 0 < 1 and unstable if ℜ 0 > 1. Proof At the equilibrium point 0 = ( , 0, 0, 0), the Jacobian matrix becomes, Therefore, its characteristic equation is, The characteristic roots are -, -, -( + ) and ( +γ) (R0 -1). Since the first three roots are negative and other will be negative if ℜ 0 < 1 and positive if ℜ 0 > 1. Thus the equilibrium point 0 is locally asymptotically stable if ℜ 0 < 1 and unstable if ℜ 0 > 1. If ℜ 0 > 1 , the epidemic equilibrium 1 = ( * , * , * , * ) is locally asymptotically stable. Proof At the equilibrium point 1 = ( * , * , * , * ), the Jacobian matrix becomes, The roots of the characteristic equations ( 1 ) are 1 = − < 0, 2 = ( + ) < 0 and other two satisfies the following quadratic equation, It is noted that the roots of (4.2) will be negative if < 0 and > 0. The disease-free equilibrium point of the system (3.4) is globally asymptotically stable if ℜ 0 < 1 and unstable if ℜ 0 > 1. Proof Considering the suitable Lyapunov function as Ƒ = . Calculative the time fractional derivative of the above function Ƒ( ) = ( ). Since, = ≤ . It follows that, Hence if ℜ 0 < 1, then Ƒ( ) < 0. By LaSalle''s extension to Lyapunov's principle, the disease-free equilibrium point is globally asymptotically stable and unstable if ℜ 0 > 1. ). By LaSalle's Invariance Principle, the epidemic equilibrium 1 is globally asymptotically stable if (2 − * * − * ) ≤ 0. For fractional order initial value problems of any variety, the Adams-Bashforth-Moulton technique is the most used numerical method. Consider the fractional differential equation as follows: ( ) = ( , ( )) , (0) = 0 , = 0,1,2, … , ⌈ ⌉, ∈ ℕ, (5.1) where 0 is the arbitrary real number, > 0 and in the Caputo interpretation, is the fractional differential operator, which is equal to the well-known Volterra integral equation We will explore into numerical simulations of the Caputo-Fabrizio Coronavirus model (3.4) in this section. From the data given in [44] , we use the total initial population of India (0) = 1382339513. Using In the instance of COVID-19 in India, the estimated parametric values are as follows: The estimated parametric values are as follows in Brazil: Fig. 9 depicts that that the number of quarantined individuals grows with time for changes in the values of . Fig. 10 depicts that the number of recovered individuals will increase with time and a decrease in the fractional order leads to a drop in the number of recovered individuals in the Brazilian population. Table 3 for = 0.9 with respect to time (days). The data fitting and model validation of the system (3.4) for Infected and Recovered cases are described in this section. The parametric values are given in Table 2 . Fig. 12 and Fig. 13 depicts the graphical representation of the infected and the recovered individuals respectivelyof the system (3.4) and the real time data of the same reported in India from 1 st July to 15 th July 2020 [44] . The results of model (3.4) thus obtained are in complete agreement with the real time data. The fractional order derivatives using the Caputo-Fabriozio of order 0 < ≤ 1 of the SIQR model were explored in this article. Based on the COVID-19 cases data in India and Brazil, collected upto 1 st August, 2020, we estimated the basic reproduction number ℜ 0 to be 1.7824 and 2.767 respectively. Modeling using fractional-order derivatives is often more efficient than modelling with integer-order derivatives because the option of derivative order gives one more degree of freedom, resulting in a better fit to real-time data with less inaccuracy than the integerorder model. We used the Iterated Laplace Transform Method to solve our proposed model and compared the results to numerical solutions obtained utilizing the Adams-Bashforth-Moulton predictor corrector technique. The parameter values in (3.4) have been estimated using the real time data given in [44] [45] and is presented in Table 2 and Table 3 . In section 4, 0 and 1 are determined, along with their stability analysis. In section 7 we observe that the proposed SIQR model with fractional order derivatives comply with the real time data in case of infected and recovered individuals. Thus, it may be considered to be an effective model to study contagious diseases. According to our understanding of the challenges, the incidence of disease transmission must be reduced, or a large part of the population would be afflicted in a very short span of time. For successful isolation and control of the disease's transmission, common preventative methods include lockdown, curfews, and the designation of containment zones. Cultivation of viruses from a high proportion of patients with colds Contribution to Mathematical Theory of Epidemics Infectious diseases of Humans: Dynamics and central The mathematical theory of infectious disease An epidemiological model with a delay and a nonlinear incidence rate An introduction to compartmental modeling for the budding infectious disease modeler COVID-19 Infection: origin, transmission, and characteristics of human coronaviruses The fractional calculus Existence of solutions to fractional order ordinary and delay differential equations and applications Quantifying undetected COVID-19 cases and effects of containment measures in Italy Transmission dynamics of fractional order typhoid fever model using Caputo-Fabrizio operator A numeric analytic method for approximating a giving up smoking model containing fractional derivatives Discrete Mittag-Leffler kernel type fractional difference initial value problems and Gronwalls inequality Fractional logistic models in the frame of fractional operators generated by conformable derivatives On fractional order cancer model Analytical approximations for a population growth model with fractional order Modeling infectious diseases On fractional order dengue epidemic model Dynamics of hepatitis c virus infection: mathematical modeling and parameter estimation A fractional order SEIR model with vertical transmission The fractional SIRC model and influenza a Impulsive initial value problems for a class of implicit fractional differential equations Two analytical methods for time-fractional nonlinear coupled Boussinesq-Burger's equations arise in propagation of shallow water waves An iterative method for solving nonlinear functional equations Mathematical Modelling and Analysis of Infectious Diseases Novel asymptotic stability criterion for fractional-order gene regulation system with time delay Global uniform asymptotical stability for fractional-order gene regulatory networks with time-varying delays and structured uncertainties Novel stability condition for delayed fractionalorder composite systems based on vector Lyapunov function A new definition of fractional derivative without singular kernel A wavelet based numerical scheme for fractional order SEIR epidemic of measles by using Genocchi polynomials Homotopy perturbation transform method for nonlinear differential equations involving to fractional operator with exponential kernel A new analytical modelling for fractional telegraph equation via Laplace transform Analysis of differential equations involving Caputo-Fabrizio fractional operator and its applications to reaction-diffusion equations The extended fractional Caputo-Fabrizio derivative of order 0 ≤ σ < 1 on CR[0, 1] and the existence of solutions for two higher-order series-type differential equations A fractional model for population dynamics of two interacting species by using spectral and hermite wavelets methods Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? A model for describing the velocity of a particle in Brownian motion by Robotnov function based fractional operator A study of fractional Lotka-Volterra population model using Haar wavelet and Adams-Bashforth-Moulton methods A study on fractional host-parasitoid population dynamical model to describe insect species The Construction of Next-Generation Matrices for Compartmental Epidemic Models Multi-order fractional differential equations and their numerical solution The authors are grateful to the reviewers for their valuable comments and suggestions. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.