key: cord-0851716-pjuxiahj authors: Butt, A.I.K.; Ahmad, W.; Rafiq, M.; Baleanu, D. title: Numerical analysis of Atangana-Baleanu fractional model to understand the propagation of a novel corona virus pandemic date: 2021-12-22 journal: Alexandria Engineering Journal DOI: 10.1016/j.aej.2021.12.042 sha: a41822440fd524253714632510caa8d5ca63f3df doc_id: 851716 cord_uid: pjuxiahj In this manuscript, we formulated a new nonlinear SEIQR fractional order pandemic model for the Corona virus disease (COVID-19) with Atangana-Baleanu derivative. Two main equilibrium points F 0 ∗ , F 1 ∗ of the proposed model are stated. Threshold parameter R 0 for the model using next generation technique is computed to investigate the future dynamics of the disease. The existence and uniqueness of solution is proved using a fixed point theorem. For the numerical solution of fractional model, we implemented a newly proposed Toufik-Atangana numerical scheme to validate the importance of arbitrary order derivative ρ and our obtained theoretical results. It is worth mentioning that fractional order derivative provides much deeper information about the complex dynamics of Corona model. Results obtained through the proposed scheme are dynamically consistent and good in agreement with the analytical results. To draw our conclusions, we explore a complete quantitative analysis of the given model for different quarantine levels. It is claimed through numerical simulations that pandemic could be eradicated faster if a human community selfishly adopts mandatory quarantine measures at various coverage levels with proper awareness. Finally, we have executed the joint variability of all classes to understand the effectiveness of quarantine policy on human population. Corona virus disease (COVID-19) is a severe pandemic with an extremely high fatality rate caused by an infectious respiratory virus called SARS-CoV-2. The virus spreads uniformly among individuals and affecting people worldwide through direct and indirect physical contact [1] . The first and overall third outbreak of novel Corona virus disease was experienced at the end of December 2019 in the Wuhan City of China. It spreads rapidly in different parts of China and then worldwide in almost 223 countries of Asia, Australia, America and Europe [4] . World Health Organization (WHO) recorded the global situation and reported that there are more than 232, 636, 622 confirmed cases with 4, 762, 089 deaths by September 29, 2021 . It is remarkable that WHO reported 408, 990 new cases to date while the number of confirmed cases are continuing to increase. However, almost 210, 951, 891 people have been recovered from COVID-19. A total of 6, 136, 962, 861 vaccine doses have been administered by 29 September, 2021 . Currently, the number of highest positive cases encountered in United States of America followed by India, Brazil and Spain respectively. 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 The Caputo fractional derivatives and integrals are the most well known operators, historically been used for modeling many real world problems. Moreover, the Riemann-Lioville derivative operator has a strong relation with the Caputo fractional derivative. When these operators are used to investigate the structure of different models, they may block obtaining better results. The main issue with the Riemann-Lioville and Caputo operators is the singularity property of their kernals [54, 55] . Due to this reason, researchers have felt the need of fractional operators with nonsingular kernals to better understand the dynamics of models. To this extent, many scientists succeeded in presenting fractional operators with nonsingular kernels, for example, Caputo-Fabrizio fractional operator [25, 54, 56] . In the recent past, Atangana and Baleanu [57] proposed a new fractional derivative operator (ABC) with a Mittag-Leffer kernel with one parameter. The main advantage of this operator is that it has a nonlocal and nonsingular kernel. It provides an advantage for those people who work in numerical modeling of real world problems. More recently, new advances and studies in fractional differential equations with ABC derivative operator has been published [28, 29, 30 ]. Among non-pharmaceutical interventions, execution of various quarantine strategies in human communities is a significant public health measure, the most effective and a safe way that historically been utilized to control the spread of all communicable diseases. Physically weak populations from the whole world especially in developing countries are at the highest priority for quarantine process. World Health Organization works closely with researchers around the world to assess the developed global standards, safety norms and the effectiveness of the implemented quarantine policy during the quarantine period. To control the spread of epidemic diseases, many researchers are attracted to develop mathematical models [10] - [16] dealing with quarantine strategies. A several number of articles [17] - [24] have been published recently on the effectiveness of implemented quarantine program. It is worth mentioning that the most effective and a safe way to control the spread of Corona infection is the prospective quarantine of exposed humans. Recently, a research article [11] has been published to understand the effect of quarantine policy on the propagation of Corona virus epidemic. Authors analyzed Corona virus disease through an integer order SEIQR epidemic model. To the best of our knowledge, no one explored a complete mathematical analysis of this model for the fractional order using Atangana-Baleanu derivative operator. We will move in this direction. In this paper, our aim is to promote the applications of Atangana-Baleanu operator to a SEIQR model [11] of a Corona virus disease. We analyzed the developed fractional order initial value problem with Atangana-Baleanu derivative operator to observe class wise transmission dynamics of the infection in a human population. Mathematical modeling of COVID-19 shall work confidentially to predict and comprehend how virus spread. These predictions will also help to understand how the Corona infection might decrease or increase in the future. The proposed model [11] incorporates a class of those humans who have quarantined through appropriate quarantine programme at their level of exposedness. We assumed that the quarantine is perfect; that is, the humans will not get infected during quarantine period and will be healthy at the end of this period. The impact of the quarantine programme on individuals with the help of fractional model will be assessed. Using some known techniques, we shall perform a complete mathematical analysis of this model with an aim of controlling the spread of COVID-19. In Chaos theory, an important area of research is to study various mathematical models with different descriptions and requirements for their numerical stability. Recently, the stability of various epidemic models [12, 62, 63, 64] of Ebola and Corona viruses with respect to the involved parameters was obtained. Some appropriate reliable numerical techniques [68, 67] to solve the governing equations were employed. The dynamics of both the diseases was studied with the help of graphic illustrations. We tend to expect that our projected SEIQR fractional model also contains all vital self-propelling properties such as: well-posedness of the model along with the boundedness and positivity of obtained solutions. It is ascertained that most of the traditional standard methods sometimes become unbounded divergent when implemented to a nonlinear system. Implementation of these numerical methods could cause major problems such as generating oscillations, bifurcations, chaos, negative solutions, or solutions converging to false equilibrium 3 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 states as time grid size increased [62, 65] . Therefore, the best approach to solve our fractional model herein is to develop Toufik-Atangana finite difference scheme [45, 46] . We will prove the dynamic consistency of the developed scheme by performing detailed numerical analysis of the proposed model. Furthermore, we will simulate computationally the dynamics of the spread of Corona virus disease over time t. The numerical analysis highlights the effect of quarantine strategy on the dynamics of COVID-19 at different quarantine levels. The effect of quarantine on susceptible and exposed class is executed reveals that number of susceptible and the exposed humans are decreased continuously by increasing the levels of quarantine. A rapid regress in the population of infected humans is observed when quarantine rates are high. It is noticed that infectivity will approach to zero as time t goes to infinity. However, the number of recovered humans are quickly increased by increasing the coverage levels of proper quarantine. An increment in number of recovered and the decrease in number of infected humans will significantly reduce the risk of Corona virus disease. Consequently, the host population will become healthy if the coverage levels and effectiveness of implemented quarantine strategy are high. Additionally, we will discuss statistically the joint variability among all populations. The paper has been organized into ten sections. The derivation of proposed compartmental model is outlined in section 2. The fractional form of the proposed Corona model is given in section 3. In section 4, we explore in details some mathematical properties of the fractional model. Existence and uniqueness of solutions is proved using some fixed point theorem. The proofs of the positivity and boundedness of solutions are done in section 5. Disease free and endemic equilibrium points of the fractional model are obtained in the section 6. Threshold parameter R 0 is calculated analytically in section 7. We proved the local and global stability of equilibrium points in section 8. Numerical analysis of the fractional model along with graphical illustrations and necessary discussions is performed in section 9. We have developed a Toufik-Atangana numerical scheme for the fractional model to approximate its solution. Moreover, the computational advantages of proposed Toufik-Atangana scheme are discussed in details. We conclude the present work in section 10. Nonlinear models [30]- [37] describing the transmission dynamics of deadly Corona virus play a significant role in the disciplines of epidemiology. These models assist the public health planners and policy makers in many ways. There is a series of various mathematical models in the existing literature [38] - [53] with different assumptions depending on the propagation mechanism of Corona virus disease. In this section, we have considered a new real world integer order SEIQR Corona epidemic model [11] . All of our efforts will be to study the impact of proposed quarantine strategy on the dynamics of COVID-19 through the fractional form of this model. Detail of variables used to model the flow pattern of Corona pandemic within the human population is given below: The total human population N (t) is separated into five classes denoted by S(t), E(t), I(t), Q(t) and R(t), all remain disjoint over the continuous time evolution. The first class denoted by S(t) includes the number of susceptible humans at time instant t. The second class which represents the number of exposed humans at time instant t is denoted by E(t). The third epidemic class characterizes number of infected humans at time instant t who are able to spread the infection and is denoted by I(t). It is not possible for any regime to quarantine the entire population, hence not everyone will be able to be quarantined right away. The fourth class which is denoted by Q(t) gives the number of exposed humans who have been quarantined at time instant t. Finally, the fifth class denoted by R(t) represents the number of recovered humans who have developed complete immunity from the Corona virus disease at time instant t. We assume that the recovered humans will remain in the class R(t) throughout their life. Thus S, E, I, Q, and R are the variables for the Corona model (1) . All of these variables are assumed to be real valued continuously differentiable functions defined on the interval [0, +∞) of t values. Fig. 1 describes a schematic which represents the involvement of Corona infection in the human population. The description of employed parameters and their physical relevance are provided in Table 1 . The SEIQR mathematical model representing the flow pattern of Corona virus in a community is given by the following system of ordinary differential equations: We have assumed that all the initial conditions and involved parameters of the system are nonnegative. We first recall some fundamental notions related to Caputo [54] and Atangana-Baleanu fractional derivatives [57] . The Caputo fractional derivative of order ρ ∈ (n − 1, n), n ∈ N, for a given function v : [a, b] → R, is stated [47, 54] as: [57] is defined as: where E ρ is the well-known one parameter Mittag-Leffler function of form: and B(ρ) > 0 is a normalization function satisfying B(0) = B(1) = 1. The Laplace Transform [57, 58] of (3) is defined as follows: The associated fractional integral [47, 57] with non-local kernal is defined by For more details about the properties of Atangana-Baleanu fractional derivative, we refer the readers to [44] - [52] . has the following unique solution [28, 29, 57 ] Generally, to analyze the dynamical behavior of an infectious disease, the classical integer order models are neither reliable nor more helpful. However, the fractional order models show cooperatively better fit to the real data. Hence, to generate the system (1) for the present work, we replace the classical integer order time derivative D t by fractional Atangana-Baleanu derivative ABC 0 D ρ t given in (3) . Through this reformulation, we will be able to observe memory effects and gain more insights about the pandemic. For t ≥ 0 and ρ ∈ (0, 1], the proposed nonlinear fractional order Corona virus model in the sense of ABC-fractional operator is given as: and it can be written in the following simple form: where v : [0, +∞) → R 5 and F : respectively. Clearly, F i ; i = 1, 2, 3, 4, 5, are functions of t, S, E, I, Q, and R. In the next section we investigate the existence and uniqueness of the solution for system (8) by fixed point theorem. In this section, we prove existence and uniqueness [45, 46] of solution with the help of theorems. Proof. On inspection, the components, is a line segment in R 5 which joins the point v 1 to v 2 as α varies from 0 to 1, and is a compact subset of R 5 . Thus, we may choose a point w ∈ L(v 1 , v 2 , α) using Mean Value Theorem, such that 7 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 where e i is the i th coordinate unit vector in R 5 and JF i (w) is a bounded linear operator. As all the partial derivatives of function F i , i = 1, 2, 3, 4, 5, are bounded, there exits a number M > 0 such that and hence F is Lipschitz continuous. then the problem (8) has a unique solution for the specific initial population if F (0, v(0)) = 0 and Proof. Suppose that F (0, v(0)) = 0 for a specific initial data. We will prove that the function v(t) satisfies Eqs. (8a) and (8b) if and only if it satisfies the relation Let v(t) satisfies Eq. (8a). Applying Atangana-Beleanu fractional integral (5) to both sides of system (8a), that is On making use of the Proposition 3.4 in [57] , we obtain the following non-linear Volterra integral equation Since Implies that the solution representation satisfies the initial data. Secondly, we prove the uniqueness of solution with the condition Let J = (0, T ) and consider the operator ψ : C(J , . 8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 The supremum norm on J , . J is: Clearly, C(J , R 5 ) along with norm . J construct a Banach space. Furthermore, we can show an important inequality proved in the following: Let us define an integral operator ϕ : Here k(t, ξ) : J × J → R is a given function called the kernal of ϕ. It is assumed to be a continuous function on the closed square G = J × J in the tξ-plane, and hence bounded. Then there exist a real number k 0 such that Our target is to show that ϕ is bounded. Consider Since k(t, ξ) and v(ξ) are continuous, therefore with On make use of the definition given in (13), we get 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 Applying the Lipschitz condition (10) along with the result in (14), we deduce that Thus, the operator ψ will be a contraction if the constant is less than one. Hence, by the Banach Contraction Principle there exists a unique solution of the model (8). In this section, we discuss some important properties of the model (7) or equivalently system (8) such as boundedness and positivity of the solutions for t ≥ 0. Thus the obtained solution will be realistic and physically realizable and the problem becomes well-posed. All these computations are done in the following steps. In this subsection, we determine the boundary of solutions (S, E, I, Q, R) of the non-linear system (7) with a non-negative initial data. Our main target is to show that the obtained feasible region in R 5 + which is positively invariant with respect to the fractional model (7). (7) is given by The existence and uniqueness of the solution of model (7) are proved in the previous section, it remains to show that the set Π is positively invariant with initial conditions S(0) ≥ 0, E(0) ≥ 0, The following theorem will be used for the proof of Theorem 4. Theorem 5 The solutions of system (7) are bounded. Proof. Addition of all the five equations in the fractional model (7) gives where Therefore, from equation (16) it follows that Applying the Laplace transform on both sides of above inequality, we obtain 10 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 ](s) and N (0) represents the initial value of the total population. Implies Applying the inverse Laplace, we arrive at where , and E α,β is the Mittag-Leffler function with two parameters α > 0 and β > 0 may be defined by the series provided that s > |Λ| 1/α . For α, β > 0, the Mittag-Leffler function satisfies and for the case α = ρ, β = ρ + 1 and z = −Ωt ρ , we have The Mittag-Leffler function is bounded for all t > 0, possess an asymptotic behavior [57] , introducing the relation (18) in inequality (17), it is obvious that all other variables of the model (7) are bounded in a region Π. The solution remains x = 0 for all t > 0 if x 0 = 0. Moreover, for any nonnegative set of initial conditions in Π, every solution of model (7) in approach asymptotically in finite time t, enters and remains in the region Π. Implies the region Π attracts all solutions in R 5 + . Thus the closed set Π is a positively invariant [12] for our system (7). 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 We have to prove that the involved state variables S, E, I, Q and R are non-negative for any t > 0, that is the solution will remain positive forever corresponding to any non-negative initial data in R 5 + . This property is required to show our model physically realizable. The solution space (S, E, I, Q, R) of the system (7) will remain positive forever with any positive initial costs. Proof. First equation of the model (7) is rearranged to give Since all the solutions are bounded, let E(t) and I(t) are bounded by σ 1 and σ 2 respectively. Then, where c = μ + β 1 σ 2 + β 2 σ 1 is a constant. Applying the Laplace transform in (19), we deduce (1−ρ) . Introducing inverse Laplace transform in the above inequality, we obtain Since both the quantities on right hand side of (20) are positive. Implies the solution S(t) is positive for all t ≥ 0. Similarly, we can easily prove that E(t) ≥ 0, I(t) ≥ 0, Q(t) ≥ 0 and R(t) ≥ 0 for all t ≥ 0 corresponding to any non-negative initial data. Thus the solutions in R 5 + remain positive forever. Proofs of existence and uniqueness, boundedness and the positivity of all solutions show that the problem (8) is mathematically well-posed. Equilibrium points of the proposed fractional model (7) are obtained by solving the system After simple calculations, the SEIQR model (7) has a unique non-negative Corona free equilibrium at the point given by 12 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 and a unique non-negative Corona present equilibrium at the point given by in the epidemiological region Π, where We can determine the dynamics of mathematical models by computing the value of threshold parameter, usually denoted by R 0 . The spread of any infectious disease begins when an infected person enters into a class of completely susceptible humans. Therefore, a threshold parameter gives the total number of new infections produced by an infected human at t ≥ 0. It works confidentially to predict the dynamical behavior of a model. Implies it provides a better information about the future spread or control of any disease. The value of threshold parameter for the model (7) is computed by using the standard next generation matrix approach [12, 62, 63, 64] . If v = (E, I, Q, S, R) T ∈ R 5 + , then the model (7) can be written as The functions, F and V for the rate of new infection terms entering, and the rate of transfer into and out of the exposed, infected, and the quarantined class, respectively are as follows: 13 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 whereas the matrix G representing the transition rate of existing or transported cases is computed as The jacobian of the matrices F and G evaluated at point F * 0 arē The maximum absolute eigenvalue of the positive matrixFḠ −1 is the value of R 0 for the model (7). That is, In this part, the local and global stabilities of a Corona SEIQR fractional model (7) 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 Proof. The Jacobian J for the system (1) at F * 0 can be written as Jacobian matrix J(F * 0 ) has the following eigenvalues: It is to be noted that the parameters in the proposed model are assumed to be positive. Therefore, the eigenvalues λ 1 ,λ 2 and λ 3 are clearly negative. Indeed the quantities μ, and q + μ + d 2 are strictly positive. Now consider Eq. (21d), that is , which is true. Lastly, it has been proved that Therefore, Eq. (21e) implies, λ 5 < 0 if and only if R 0 < 1. Thus all the eigenvalues of J(F * 0 ) when R 0 is less than unity, are negative and hence F * 0 is locally asymptotically stable. Once the dynamical system is stable at disease free point F * 0 , then there will be no pandemic threat. Proof. The Jacobian for the system (1) at F * 1 can be written as 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 where Then the Jacobian matrix J(F * 1 ) has the following eigenvalues: On inspection, it is verified that Implies that λ 4 and λ 5 are negative. Hence F * 1 is locally asymptotically stable. Theorem 9 A Corona free equilibrium F * 0 is globally asymptotically stable (GAS) in the domain Π when R 0 < 1. We construct a candidate [48] Lyapunov function U 0 : Π → R such that, Substituting the values of fractional derivatives from system (1), we have Hence by LaSalle's invariance principle [69] , F * 0 is globally asymptotically stable in Π. As a result, Corona virus disease disappears from the human population. Π if R 0 > 1. To show the global stability of a Corona present steady state F * 1 , we construct a proposed [48] Lyapunov function U 1 : Π → R such that The Atangana-Baleanu fractional derivative of U 1 with respect to time t can be written as Using equations given in the model (1), we obtain that ABC 0 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 and Since all the parameters used in the model (1) are non-negative, we have ABC Hence by LaSalle's invariance principle [69] , it is concluded that F * 1 is globally asymptotically stable in Π. Global stability of F * In this segment, we employ Toufik-Atangana scheme [45, 46] to acquire a numerical solution of fractional model for different values of ρ. The numerical values of involved parameters, given in Table 1 are taken from [11, 12] . Dynamical behavior of Corona virus disease over time t is simulated for various values of fractional order ρ. Moreover, the effectiveness of quarantine programme is analyzed numerically for different coverage levels. The presented numerical results are discussed in detail. In this section, we develop a numerical scheme which is based on a recently developed Toufik-Atangana rule to analyze and predict the numerical stability of a Corona fractional model (8) . To get an iterative scheme, we firstly describe the method briefly and then apply it to the fractional model (8) . By applying the fundamental theorem of fractional calculus, the system (8) can be written as 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 The function F (ξ, v(ξ)) can be approximated over [t j , t j+1 ], using the interpolation polynomial Calculating these integrals we finally get the approximate solution as: Hence, we obtain the following recursive formulae for the model equations: The purpose of this segment is to explain the role of fractional order ρ on the transmission patterns and control of COVID-19 through the model (7) . To this end, we present some graphical results 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 of the proposed fractional model (7) using an iterative finite difference scheme (24) developed by Toufik and Atangana [57] . In Fig. 2 , the effect of arbitrary fractional order ρ on total number of humans in each class is demonstrated. The dynamics of the model (7) is simulated for ρ = 0.2, 0.4, 0.6, 0.8, 1 respectively. It is remarkable that we were at the endemic equilibrium F * 1 for the integer case ρ = 1. A significant decrease in total Corona exposed and the infected cases is observed by decreasing the value of order ρ. However, the total number of humans in susceptible, quarantined and the recovered class increases slowly by reducing the value of ρ from 1. Implies we can reduce the transmission of corona infection in human population by decreasing the value of order ρ of the fractional model. It is noticed that as ρ gets close to zero from the right, the memory effects of the epidemic model decreases and consequently the number of cases in E and I gets close to 0. Moreover, for the case ρ = 1, we obtained the same graphical results as given in [11] where NSFD numerical scheme [68] , was used. In this subsection, we have studied the impact of quarantine on the dynamics of Corona virus pandemic quantitatively with the help of proposed fractional model (7) . The aim was to handle the spread of COVID-19 in a human population through persistent quarantine program. The effect of parameter q 1 (quarantine level) and fractional order ρ versus the total humans in each class is analyzed in Figs. 3-8. The study indicates that there is a significant decrease in the corona exposed and infected humans by increasing the order of fractional derivative and the level of quarantine. It is also noticed that the number of susceptible humans increases as the value of ρ and q 1 increases. Moreover, in Fig. 5 , when ρ = 0.6, the number of recovered humans gradually increases for the quarantine levels q 1 = 0, q 1 = 0.1 and q 1 = 0.3 respectively. The curves of recovered humans eventually goes down for higher quarantine levels over time t, as needed. With the continuous increase in the value of fractional order and quarantine level, the number of exposed, infected, and the recovered humans finally goes to zero whereas the number of susceptible humans goes to the value S * 0 = λ/μ for the case ρ = 1 and q 1 = 0.7. This is a particular state of corona free point F * 0 at which the total population will become healthy. It shows that at all quarantine levels ranging 0% to ≤ 70%, the Corona virus disease will persist in the human population. Indeed, from 70% onward, each population approach the equilibrium state of F * 0 automatically as illustrated in Table 2 . Hence, our study indicates that Corona virus disease could be eradicated if we quarantine at least 70% of the human population, which justifies the title of our paper. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 We preferred Toufik-Atangana scheme over all standard and non-standard methods that solve a given problem only for an integer order. The preference is justified taking into account several numerical aspects of the fractional model, as it works gently for non-integer values of order ρ. Indeed, the fractional order ρ affects the global dynamical behavior of COVID-19 pandemic and provides a detailed information about the model at least in this work. We have studied efficiently the effect of fractional order ρ and the proposed quarantine strategy on disease dynamics over the longer time. Importance of graphic results is attributed to the numerical scheme (24) In the current subsection, the joint variability [12, 62] of all classes is studied numerically. For each pair of populations, we have computed the relationship quantity ρ * and its consequences as illustrated in Table 3 . A numerical study performed in Figs. 3-7 authentically justified the inverse relationship among exposed, infected, recovered and the quarantined humans. However, it is noticed that there is a direct correlation between the susceptible and quarantined class as demonstrated in each of the Figs. 3-7. Subsequently, a continuous rapid increment in the susceptible class will lead the entire human population to Corona free state of F * 0 . Populations Table 3 : Joint variability of populations. A mathematical analysis of the proposed model (7) suggests that we will have to quarantine at least 70% of the human population to control the spread of COVID-19. It looks difficult for all human communities, especially in developing countries due to their limited resources. To make the quarantine strategy effective, firstly, we need to prepare priming people who can help others to think about many ways to reduce issues of loss of freedom, boredom, irritability or enxiety during quarantine period. We should encourage people to make plans such as stucturing of time for virtual socialising and for exercise, or other shared activities. Every government or agency should plan for what will happen as and after quarantine restrictions are lifted. We need to provide a reliable internet access to psychological, social and medical support. Confidential telephone hotlines that provide professional counseling should be provided. The state needs to provide a clean and a safe place to quarantine along with basic supplies such as clothes, water, food, and financial support. The basic necessities such as payment for house rent, electricity, child care and other utilities should be provided. We need to provide people the ability to work remotely. 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 In this study, we assessed the dynamics of a novel Corona virus disease by developing a fractional order SEIQR model. We have used Atangana-Beleanu derivative operator with generalized Mittag-Leffler function to formulate the proposed fractional model. The Atangana-Beleanu derivative operator was preferred due to its non-local and non-singular kernel. The purpose of our study was to focus on a SEIQR model to depict the prevalent characteristics of Corona virus. To confirm the existence and uniqueness of solution of the proposed fractional model, the fixed point theorem is applied. Threshold parameter is computed theoretically and numerically to handle the dynamical behavior of the model. Disease free and endemic equilibrium points were obtained analytically. It is proved that both the equilibria are locally and globally stable in the feasible region. Global stability of equilibrium points is analyzed by the theory of Lyapunov functions. The solution of the model is then investigated by implementing a recently proposed numerical scheme from the existing literature. All the results obtained are simulated by figures and seen to be compatible with the model. We demonstrate the long term dynamical behavior of COVID-19 with the help of numerical simulations. The purpose was to investigate the impact of different values of fractional order ρ. Through the fractional analysis, it is observed that number of susceptible humans increases as the value of fractional order ρ decreases. However, as ρ decreases, the number of humans in exposed, infected, quarantined and the recovered class also decreases. Hence, Corona virus disease decreases slowly and can be controlled by reducing the value of fractional order ρ from 1. These fractional results show that fractional models will be more effective than integer-order ones in truly estimating the optimal quarantine level. As a second approach, it was suggested that the number of Corona infective cases can be reduced by increasing the value of fractional order ρ and levels of quarantine q 1 . Consequently, it is deduce that SEIQR fractional model is more efficient in fitting real data than integer-order SEIQR model. Finally, the joint variability of the quarantined class and the remaining classes is executed statistically. The obtained relationship along with the consequences among all pairs were also observed in all numerical simulations, one can validate our study. Our numerical investigations prove that an earliest proper quarantine measure helps in controlling the Corona pandemic. To make it possible, the government especially in developing countries will have to take serious steps like providing less expensive suitable quarantine strategy to control the Corona virus disease. Mathematical analysis of the fractional model given in this paper makes an attempt to understand the dynamical behavior of a disease. To combat COVID-19 in the human population, we assure that the investigations in the current research work will be beneficial for the decision making and health authorities. In future, we will present a comprehensible depiction of Corona virus disease using different intervention strategies through SEIQHR fractional model with ABC operator. A general fractional order optimal control problem will be analyzed to find the best controls for the employed quarantine and hospitalization measures. Funding is not applicable for this manuscript. The authors declare that they have no known competing interests. 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 It is found that number of susceptible and the recovered humans gradually increases by decreasing the value of order ρ from 1. However, with the decrease in the value of fractional order ρ from 1, the number of exposed and infected humans also decreases. The curves for each of the class S, E, I, Q, R, and N get flattened as ρ reduces from 1 to 0.6. 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 Profiles for each class are obtained using Toufik-Atangana numerical scheme with step size h = 0.01. It is seen that number of susceptible and recovered humans gradually increases by increasing the level of quarantine. However, the number of exposed and infected humans decreases with the increase in quarantine rate q 1 . 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 Profiles for each class are obtained using Toufik-Atangana numerical scheme with step size h = 0.01. It is seen that number of susceptible and recovered humans gradually increases by increasing the level of quarantine. However, the number of exposed and infected humans decreases with the increase in quarantine rate q 1 . 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 Profiles are obtained using Toufik-Atangana numerical scheme with step size h = 0.01. It is seen that number of susceptible humans gradually increases by increasing the level of quarantine whereas the number of exposed and infected humans decreases. One important thing is to be noted that number of recovered humans initially increases (see for q 1 = 0, 0.1, 0.3) and then reduces (see for q 1 = 0.5, 0.7) with the increase in quarantine rate. 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 Profiles are obtained using Toufik-Atangana numerical scheme with step size h = 0.01. It is seen that number of susceptible humans gradually increases by increasing the level of quarantine whereas the number of exposed and infected humans decreases. One important thing is to be noted that number of recovered humans initially increases (see for q 1 = 0, 0.1, 0.3) and then reduces (see for q 1 = 0.5, 0.7) with the increase in quarantine rate. 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 Profiles are obtained using Toufik-Atangana numerical scheme with step size h = 0.01. It is seen that number of susceptible humans gradually increases by increasing the level of quarantine whereas the number of exposed and infected and the recovered humans decreases. The graph shows that we are at disease free state F * 0 when round about 70% of the human population is quarantined. 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 World of Health Organization The continuing 2019-nCoV epidemic threat of novel corona viruses to global health: the latest 2019 novel coronavirus outbreak in Wuhan, China A discrete stochastic model of the COVID-19 outbreak: forecast and control Middle East respiratory syndrome coronavirus: risk factors and determinants of primary, household, and nosocomial transmission Handbook of COVID-19 prevention and treatment, The First Affiliated Hospital COVID-19 and thrombotic or thromboembolic disease: implications for prevention, antithrombotic therapy, and follow-up Care for critically ill patients with COVID-19 Animal communication: when I am calling you, will you answer too The effect of control measures on COVID-19 transmission in Italy: Comparison with Guangdong province in China Isolation, quarantine, social distancing and community containment: pivotal role for old-style public health measures in the novel coronavirus (2019-ncov) outbreak Design of a nonlinear model for the propagation of COVID-19 and its efficient nonstandard computational implementation The impact of quarantine on Covid-19 infections Modelling the effects of media coverage and quarantine on the COVID-19 infections in the UK A model for COVID-19 with isolation, quarantine and testing as control measures Modeling the dynamics of novel coronavirus (COVID-19) via stochastic epidemic model Control strategies for COVID-19 epidemic with vaccination, shield immunity and quarantine: A metric temporal logic approach Simplified model of Covid-19 epidemic prognosis under quarantine and estimation of quarantine effectiveness On a quarantine model of Coronavirus infection and data analysis Global stability of COVID-19 model involving the quarantine strategy and media coverage effects On confinement and quarantine concerns on an SEIAR epidemic model with simulated parameterizations for the COVID-19 pandemic How quarantine and social distancing policy can suppress the outbreak of novel coronavirus in developing or under poverty level countries: a mathematical and statistical analysis 2020 Modeling COVID-19 pandemic with hierarchical quarantine and time delay A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator On q-analogue of meromorphic multivalent functions in lemniscate of Bernoulli domain Reassessments of gross domestic product model for fractional derivatives with non-singular and singular kernels Discrete fractional differences with nonsingular discrete Mittag-Leffler kernels Fractional operators with generalized Mittag-Leffler kernels and their iterated differintegrals Stability analysis and simulation of the novel Corornavirus mathematical model via the Caputo fractional-order derivative: A case study of Algeria On the iterative methods for solving fractional initial value problems: new perspective A fractional-order model of COVID-19 considering the fear effect of the media and social networks on the community Stability analysis of fractional nabla difference COVID-19 model Effect of weather on the spread of COVID-19 using eigenspace decomposition A Caputo power law model predicting the spread of the COVID-19 outbreak in Pakistan A numerical and analytical study of SE(Is)(Ih)AR epidemic fractional order COVID-19 model Advances in Difference Equations (2021) 2021:293 Selection of an effective hand sanitizer to reduce COVID-19 effects and extension of Topsis technique based on correlation coefficient under neutrosophic hypersoft set On the necessity of proper quarantine without lock down for 2019-nCoV in the absence of vaccine Integration by parts and its applications of a new nonlocal fractional derivative with Mittag-Leffler nonsingular kernel On fractional derivatives with exponential kernel and their discrete versions A new fractional model for vector-host disease with saturated treatment function via singular and non-singular operators Mathematical modeling for novel coronavirus (COVID-19) and control Fractional differential equations with Atangana-Baleanu fractional derivative: Analysis and applications New numerical approximation of fractional derivative with nonlocal and non-singular kernel: application to chaotic models Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative On some new properties of fractional derivatives with Mittag-Leffler kernel Analysis of Atangana-Baleanu fractional-order SEAIR epidemic model with optimal control A Fractional SAIDR Model in the Frame of Atangana-Baleanu Derivative, Fractal and Fractional Analysis of some generalized ABCfractional logistic models On the nonlinear dynamical systems within the generalized fractional derivatives with Mittag Leffler kernel Fractional order mathematical modeling of covid-19 transmission Optimal control and sensitivity analysis for transmission dynamics of Coronavirus Fractional differential equations: an introduction to fractional derivatives. In: fractional differential equations, to methods of their solution and some of their applications Fractional integrals and derivatives: theory and applications A new definition of fractional derivative without singular kernel New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model Chaos in a simple nonlinear system with Atangana-Baleanu derivative with fractional order A fractional order mathematical model for COVID-19 dynamics with quarantine, isolation, and environmental viral load The dynamics of COVID-19 with quarantined and isolation Modeling of fractional-order COVID-19 epidemic model with quarantine and social distancing Effect of quarantine on transmission dynamics of Ebola virus epidemic: a mathematical analysis Mathematical analysis to control the spread of Ebola virus epidemic through voluntary vaccination A reliable and competitive mathematical analysis of Ebola epidemic model Ebola virus epidemic disease its modeling and stability analysis required abstain strategies A mathematical analysis of an isothermal tube drawing process Numerical based approach to develop analytical solution of a steady-state melt-spinning model Nonstandard Finite Difference Models of Differential Equations, World Scientific publishing The Stability of Dynamical Systems The authors are grateful to the anonymous reviewers and the handling editor for their constructive comments and suggestions, which lead to improvements of the original manuscript.