key: cord-0955829-1992iers authors: Babaei, A.; Ahmadi, M.; Jafari, H.; Liya, A. title: A mathematical model to examine the effect of quarantine on the spread of coronavirus date: 2020-11-30 journal: Chaos Solitons Fractals DOI: 10.1016/j.chaos.2020.110418 sha: 9f57fd2c2a3091c4b8497500712acba8f08e91c3 doc_id: 955829 cord_uid: 1992iers In this study, we propose a mathematical model about the spread of novel coronavirus. This model is a system of fractional order differential equations in Caputo’s sense. The aim is to explain the virus transmission and to investigate the impact of quarantine on decreasing the prevalence rate of the virus in the environment. The unique solvability of the presented COVID-19 model is proved. Also, the equilibrium points and the reproduction number of the proposed model are discussed in two cases with and without considering the quarantine factor. Using the Adams-Bashforth-Moulton predictor-corrector method, some numerical simulations are implemented to survey the behavior of the considered model. The current crisis that many countries are facing is a new virus called coronavirus. In late December 2019, an unknown virus was reported in the city of Wuhan, China. The main concern of Chinese officials was that the number of infected people was increasing exponentially, and this created an epidemic issue 5 in Wuhan [1] . Chinese immunologists who have dealt with this emerging phenomenon believe that the new virus can easily spread in public. The main way to the transmission of coronavirus is through respiratory droplets when people are near to each other [2, 3] . Therefore, quarantine and social distancing seem to be the only appropriate control mechanisms, until a vaccine or some drugs 10 are found for coronavirus disease 2019 . Mathematical biology is one of the most interesting research areas for applied mathematicians. Many theoretical and computational studies are done by the scientists in this field [4, 5, 6, 7, 8] . In recent years, using fractional order operators have provided new aspects for describing mathematical mod- 15 els in biomathematics. The definition of fractional order operators can preserve hereditary and memory traits of a considered variable in a real problem [9, 10, 11, 12, 13, 14] . In many cases of natural biological processes, the present and next state of a system are dependent on its all previous states. Thus, fractional operators are suitable and valuable mathematical tools to a better un- 20 derstanding the behavior of natural systems. This helps researchers to propose more accurate models of various biological phenomena. Authors in literature employed different types of fractional derivatives for these purposes. Caputo and Riemann-Liouville fractional derivatives are the most used definitions by researchers [15, 16, 17, 18, 19, 20, 21] . In recent years, some new definitions of 25 fractional derivatives have been proposed. Caputo [27, 28] used this fractional operator to study some other natural phenomena. COVID-19 causes significant damage to the economies of many countries. Therefore, it is vital to find a working solution to prevent the spread of this virus and to control this disease. That is why many people around the world are working in different disciplines looking for a useful way to control the virus 35 efficiently. In recent months, the main focus of mathematical biology specialists has been on the problems related to this pandemic. These researchers presented some models to study the important factors of virus transmission. By studying these factors, they are trying to take an essential step in this field. Thus, many mathematical models are presented to survey the dynamics of COVID-19 40 infection [29, 30, 33, 34, 35, 36, 37] . In this study, we introduce two fractional order models to analyze the behavior of COVID-19 in society. The main difference between these models is the quarantine factor. First, the model is surveyed according to this factor. Afterwards, the model is investigated without considering this factor. The rest of this paper has the following organization. In Section 2, we provide an initial model of COVID-19 in the presence of the quarantine process and analysis of the model will be examined. This model consists seven compartments, each of them is somehow involved in the virus. In Section 3 we analyze a sub-model of the previous section without considering the quarantine 50 factor. Some numerical simulations are implemented in Section 4 to investigate the effect of quarantine restrictions on the spread of coronavirus. Finally, the main findings and conclusions are presented in Section 5. In this section, we propose a fractional order model for studying the quar-55 antine factor on coronavirus prevalence. Therefore, at first, we review the definition of Caputo fractional derivative [9] as In the following, we use D α insted of C 0 D α t for simplicity. In the survey conducted, we will have seven compartments. People who are susceptible or exposed to the virus are shown by S. These individuals are 60 divided into the following categories: • infected but do not yet have symptoms (A), • infected with symptoms (I), • infected who are quarantined (E q ), • infected who are hospitalized (H), • people who have recovered from the disease (R), • the quarantined susceptible individuals (S q ), and To present the model, the following parameters are introduced. Λ is the birth rate of the population. c is the contact rate. θ shows the transmission rate and β is the probability of transmission per contact. q indicates quarantined rate 70 of susceptible individuals. δ I and δ q represent transition rates of symptomatic infected individuals and quarantined susceptible individuals to the quarantined infected compartment, respectively. γ I , γ A and γ H are the recovery rates of infected people with symptoms, infected people without symptoms and quarantined infected individuals, respectively. µ shows naturally death rates. ρ is 75 the probability of having symptoms among infected individuals. λ shows the release rate of the quarantined uninfected contacts into the community. Finally, α I , α A , α E and α R represent the rates at which the virus removes from the compartments I, A, E and R. According to the above defined parameters and the relations between the 80 considered human categories, we get the following system: The right-hand side of the system (2) has the physical dimension (time) −1 whereas the other side has the dimension (time) −α . To correct this mismatch, we use the approach proposed in [15] . So, we get the modified system as (3) To examine the unique solvability of system (3), first we recall the following 85 lemma. where a ≤ ζ ≤ t, ∀t ∈ (a, b]. Thus, on each hyperplane bounding the nonnegative orthant, the vector field points to Ω + . So, Corollary (1) results Ω + is a positive invariant set. To verify the stability of E 0 , the basic reproduction number denoted by R 0 should be computed. For this aim, the next generation matrix method [5] will be employed. First, based on this method, we define the matrices F and V associated with the model (3) as The required basic reproduction number of model (3) is the spectral radius of the matrix F V −1 . Hence, we get in which Proof. First, we get the Jacobian matrix related to (3) at E 0 as: The DFE is stable if all eigenvalues of the Jacobian matrix J(E 0 ) be negative. The eigenvalues of this matrix are as follows: Since the first five eigenvalues are negative, so it suffices to prove that λ 6 < 0 and λ 7 < 0. Let Then, according to the definition of R 0 , we can rewritten λ 6 in the form: if R 0 < 1 then K 2 > K 2 + 4ηξ(R 0 − 1) and so, it can be concluded λ 6 < 0 if R 0 < 1. According to (6) and (8), for λ 7 , we and the proof is complete. In this section, we consider a sub-model of (3) without considering the quarantine factor and discuss its properties. In fact, we want to see what happens 120 for the basic reproduction number of COVID-19 model in the absence of quarantine. It helps us to better realize the impact of quarantine on the spread of disease. Thus, we can get an apparent comparison between two cases, with and without considering quarantine facilities. For this end, by deleting the variables and parameters related to the quarantine process from (3), we get The category S q of (3) is a subset S. Also, E q in (3) is a subset from the union of A and I. So, the model (3) and the sub-model (10) have the similar main compartments. As a result, the reproduction numbers of these models are the same and obtained as (5). Let E * = (S * (t), I * (t), A * (t), H * (t), R * (t)) is the endemic equilibrium of this 130 model. This point can be determined by solving the system The third equation leads to Put y = I * + θ α A * . Then, we get The first equation in (11) results Substituting (14) into (13) and some direct calculations give Now, from the second equation in (11), we have This relation leads to Hence, due to (15) and (17), the following result can be obtained With substituting (17) into (14) 140 Now, the fourth equation in (11) results Finally, from the fifth equation in (11) and Eqs. (18) , and (20), we can conclude Therefore, all the components of E * can be displayed in terms of I * . Theorem 3. The model (10) has an unique endemic equilibrium if and only if Proof. According to (17) and (19), from (16) and with doing some simple calculations, we have . Hence, rewriting this equation due to the obtained reproduction number R 0 in (5) results So, substituting (22) into Eqs. (18)- (21) gives The above equations show that E * is the endemic point of (10) if and only if R 0 > 1. Table 1 . Also, we suppose the following In the present article, a model for coronavirus disease was proposed to show 190 the effect of quarantine on the spread of the virus. Existence of a positive solution for this model was proved and the stability of the proposed model was studied. Afterwards, to investigate the impact of quarantine, a sub-model was obtained without considering this factor. It was shown that for the reproduction numbers larger than one an endemic equilibrium point exists for this sub-model. The graph of the reproduction number versus the fractional derivative order was plotted for different values of quarantine rate. It was concluded that for smaller values of quarantine rate, the reproduction number has very significant growth. Also, we can observe that the basic reproduction number tends to the values less than 1 when the fractional order decreases. So, for smaller values of α, we [8] J. Amador, A. Gomez-Corral, A stochastic epidemic model with two quar- Novel CoronavirusChina, Disease Outbreak News: Update Interim infection prevention and control recommendations for patients with suspected or confirmed coronavirus disease 2019 (COVID-19) in healthcare settings Emerging coronaviruses: genome structure, replication, and pathogenesis antine states and limited carrying capacity for quarantine Theory and Applications of Fractional Differential Equations Fractional dynamics Fractional Calculus Models and Numerical Methods On the new fractional hy-245 brid boundary value problems with three-point integral hybrid conditions Mathematical analysis of a fractional differential model of HBV infection with antibody immune response Chaos A mathematical model of tuberculosis (TB) transmission with children and adults groups: A fractional model A new fractional analysis on the interaction of HIV with CD4 + T-cells Numerical study for fractional model of nonlinear predator-prey biological population dynamical system A fractional order HIV/AIDS model based on the effect of screening of unaware infectives An efficient computational approach for a fractional-order biological population model with carrying capacity A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions An efficient numerical method for fractional SIR epidemic model of infectious disease by using Bernstein wavelets A new definition of fractional derivative without singular kernel A new fractional SIRS-SI malaria disease model with application of vaccines, antimalarial drugs, and spraying A new study on the mathematical modelling of human liver with Caputo-Fabrizio fractional derivative Analysis of the model of HIV-1 infection of CD4 + T-cell with a new approach of fractional derivative New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model Decolonisation of fractional calculus 290 rules: breaking commutativity and associativity to capture more natural phenomena A new analysis of fractional fish farm model associated with Mittag-Leffler-type kernel Modelling the effects of media coverage and quarantine on the COVID-19 infections in the UK The role of asymptomatic class, quarantine and isolation in the transmission of COVID-19 Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions A mathematical model of covid-19 us-305 ing fractional derivative: Outbreak in india with dynamics of transmission and control Beware of asymptomatic transmission: Study on 2019-nCov preventtion and control measures based on SEIR model A fractional differential equation model for the COVID-19 transmission by using the Caputo-Fabrizio derivative Modelling the spread of COVID-19 with new fractalfractional operators: Can the lockdown save mankind before vaccination? Dynamic analysis of the mathematical model of COVID-19 with demographic effects HIV and shifting epicenters for COVID-19, an alert for some countries Generalized Taylors formula Global existence theory and chaos control of fractional differential equations A predictor-corrector approach for the numerical solution of fractional differential equations On linear stability of predictor-corrector algorithms for fractional differential equations Conceptualization, Methodology, Software, Writing -review & editing A. Ahmadi : Data curation, Writing -original draft Visualization, Investigation, Writing -review & editing Writing -original draft, Software, Validation The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported 335 in this paper.