key: cord-0818168-ull9uv2s authors: Raza, Ali; Ahmadian, Ali; Rafiq, Muhammad; Salahshour, Soheil; Ferrara, Massimiliano title: An Analysis of a Nonlinear Susceptible-Exposed-Infected-Quarantine-Recovered Pandemic Model of a Novel Coronavirus with Delay Effect date: 2020-12-28 journal: Results Phys DOI: 10.1016/j.rinp.2020.103771 sha: 4336cf6004caf87a146c47b1ad4f14091f275637 doc_id: 818168 cord_uid: ull9uv2s In the present study, a nonlinear delayed coronavirus pandemic model is investigated in the human population. For study, we find the equilibria of susceptible-exposed-infected-quarantine-recovered model with delay term. The stability of the model is investigated using well-posedness, Routh Hurwitz criterion, Volterra Lyapunov function, and Lasalle invariance principle. The effect of the reproduction number on dynamics of disease is analyzed. If the reproduction number is less than one then the disease has been controlled. On the other hand, if the reproduction number is greater than one then the disease has become endemic in the population. The effect of the quarantine component on the reproduction number is also investigated. In the delayed analysis of the model, we investigated that transmission dynamics of the disease is dependent on delay terms which is also reflected in basic reproduction number. At the end, to depict the strength of the theoretical analysis of the model, computer simulations are presented. Coronavirus is a group of viruses that causes illness like the common cold, respiratory syndrome, and severe acute respiratory syndrome. A novel coronavirus (COVID-19) firstly saw in humans. This is the new one that has not been seen in humans previously [1] . The antibiotic has been discovered as a novel covered RNA coronavirus-2 that has been currently named as SARS coronavirus-2 (SARS-COV-2) which has emergent evolution similarly to SARS-COV-3. Infected people have been recorded both in the infirmary and in-home. The world health organization (WHO) has recently declared coronavirus disease 2019 (COVID-19) as a pandemic. Due to the rapid spread of COVID-19, we established that an updated analysis of cases throughout the world might help to identify the dynamics of coronavirus. Many peoples have died around the world due to this novel coronavirus [2] . Arif et al. [3] [4] investigated the stochastic analysis of typhoid fever and hepatitis B diseases. Baleanu et al. [5] presented the dynamics of influenza which is basically the symptoms of coronavirus, in the human population with constant vaccination strategy. Li et al. [6] studied the classification of brain waves in the human population. Lin et al. [7] presented the coronavirus model in the population of China under the assumption of an individual's reactions and government action. Shereen et al [8] investigated the characteristics of coronavirus in the human population. Shim et al. [9] studied the dynamics of coronavirus in the population of South Korea. Tahir et al. [10] presented the dynamics of middle east respiratory syndrome in the human population. Zhao et al. [11] studied control techniques regarding the pandemic of coronavirus. World Health Organization (WHO) declared the outbreak of coronavirus is a global issue. Due to this pandemic, this project is a platform for all researchers to introduce new delayed mathematical modeling for a better understanding of the extinction and persistence of coronavirus. Thus, the biologically consistent mathematical modeling for the transmission and dynamics of COVID-19 is an extensive need around the world. Ivorra et al. [12] investigated the mathematical modeling of coronavirus disease by taking undetected infections. Yang et al. [13] presented a mathematical model for novel coronavirus the case study Wuhan, China. Sameni [14] presented the mathematical modeling of different epidemic diseases and discussed a special case about Covid-19. Rajagopal et al. [15] investigated the fractional order model for the outbreak of coronavirus and predict different results in the field of a fractional derivative. Baleanu et al. [16] presented the coronavirus model by using the Caputo-Fabrizio derivative in the field of a fractional derivative. Goufo et al. [17] investigated the joint dynamical analysis of the equilibria of the mathematical model of HIV and COVID-19. Kouidere et al. [18] studied the dynamics of the coronavirus pandemic model with an optimal control strategy and its transmission. Mandal et al. [19] presented the prediction and control techniques for the coronavirus model with the concept of mathematical modeling. Naveed et al. [20] designed the compartmental modeling of the coronavirus pandemic with delay strategies in the subpopulation of humans as susceptible-exposed-systematical infected-a systematical infected and recovered. Shatanawi et al. investigated the nonstandard computational method on the dynamics of stochastic coronavirus model in which main focus is positivity, boundedness, consistency and stability of aforementioned purpose [24] [25] . Atangana et al. gives an authentic idea for the modeling of coronavirus pandemic in which sea food market, asymptomatic and symptomatic humans are the main sources for the prevalence of virus. How the effective use of facemask may eradicate the dreadful virus [26] . On the other way, great idea was delivered by scientist in which the use of fractal fractional derivative technique and how the lack down strategy may help to cure the virus and better than any type of vaccination [27] . Atangana et al. studied the dynamics of Ebola virus homologic fever in the west countries of Africa. Also, give the suitable precautions for the treatment of such dreadful infection [28] . Raza et al. studied the effect of delay tactics on the transmission dynamics of HIV/AIDS [29] . Due to the unavailability of vaccines for the coronavirus worldwide, delay factors such as social distance, quarantine, travel restrictions, extended holidays, hospitalization, and isolation have contributed to controlling the coronavirus epidemic. That is the reason, we motivated to do mathematical modeling of coronavirus pandemic model with delay strategy. So, we have investigated the delay differential equations in this study due to the delay term. The flow of paper has based on the following sections: In section 2, described the invention of the nonlinear delayed model and its equilibrium points. In section 3, the local stability of the proposed model is discussed. In section 4, the global stability of the proposed model is presented. Before closing, in section 5, computer simulations for the solution of susceptible-exposed-infectedquarantine-recovered pandemic models with a different strategy. In the end, the conclusion of the study is presented. In this section, we considered the total population of humans, presented as . Define a function of population dynamics , for all . The susceptible component of humans is represented as , the exposed component of humans is represented as , the infected component of humans is represented as , the quarantine component of humans is represented as and the recovered component of humans is represented as , as desired. the infection rate of exposed humans who are interaction with susceptible humans), (denotes the rate 1 of exposed humans who are directly moving to quarantine component), (denotes the rate of natural immunity), (denotes the rate of exposed humans who are infected), (denotes the rate of infected humans who are recovered after vaccination or quarantine), (denotes the rate of death of infected 1 humans due to corona virus), (denotes the rate of quarantine humans who are recovered) (denotes 2 the rate of quarantine humans who dies after quarantine because they have a weak immune system or having diseases) and (denotes the natural death rate of humans components). After that, the dynamics with nonlinear delayed equations are as follows: , . (1) . Also, the feasible region 0 = (0), 0 = (0), 0 = (0), 0 = (0), 0 = (0) of the equations (1) to (5) is . Note that, the solutions of a given system are positive and bounded, and lies in the feasible region. At the analysis that the system (1-5) has three types of equilibria as follows: Trivial equilibrium In this section, we employ the next generation matrix (NGM) method to the system (1) to (5), for obtaining reproduction number with calculating the transmission and transition matrices as follows: After that, transmission matrix and transition matrix , along with virus absenteeism equilibrium (VAE) is , 1 where, and V= . Notice that, the spectral radius of , is called reproduction number and . Notice that, two eigenvalues are repeated as, and third eigenvalue is 1 = -< 0, 2 = -< 0 . where, , , . -= 1 ( + + + 1 ) = 2 ( + + 1 ) = 3 By using the second order Routh-Hurwitz Criterion as, , if , 3 + 2 -2 1 > 0 ( + + + 1 ) + ( + + 1 ) - Proof: Considering the Jacobean matrix (JM) for the system (1-5) at is estimated as follows: [ Notice that, the eigen values are, 1 = -< 0, 2 = -( + + 2 ) < 0 ( 1 3 5 -2 5 + 3 5 -4 ) 0 > 1 Thus, we have concluded that all eigenvalues are negative and by Routh Hurwitz (RH) criteria, the given equilibria, is locally asymptotical stable. 2 In this section, we satisfied the well-known theorems at the both equilibria as follows: ( + + 1 ) + 1 -( + + 2 ) + + + - only if At the end, that trajectory of the ⇒ ≤ 0 0 < 1 = 0 = 0 , = = = = 0. system (1) to (5) is unique solution on which is . Consequently, with the assumption of Lasalle = 0 1 invariance principle (LIP), is globally asymptotically stable (GAS). 1 Theorem: The virus incidence equilibrium (VIE), is globally asymptotical stable 2 = ( 1 ,E 1 ,I 1 ,Q 1 ,R 1 ) Proof: Considering the Volterra Lyapunov function , as well-defined: :Ω→ = ( -1 -1 log 1) + ( -1 -1 log 1) + ( -1 -1 log 1) + ( -1 -1 log 1) + , . , and only if At the end, that ⇒ ≤ 0 0 > 1 = 0 if = 1 , = 1 , = 1 , = 1 , = 1 . trajectory of the system (1) to (5) is unique solution on which is . Consequently, with the = 0 2 assumption of Lasalle invariance principle (LIP), is globally asymptotically stable (GAS). 1 In this section, we discussed the simulations of the Eq. (1) to Eq. (5), by using the physical values of parameters presented in Tab. 2. In this part, we aim to present the graph of virus absenteeism equilibrium (VAE) at the given real data as the value of the threshold/ reproduction number is = 0.1702 < 1. Then, the system converges 0 to the . Figure 1 exhibits the solution of the system (1) to (5) with the absence of 1 = ( Λ , 0, 0, 0, 0) delay by using the values of the parameters presented in Tab.2 and touches to true equilibria of the model, as desired. In this example, we depict the graph of virus incidence equilibrium (VIE) at the given real data as the value of the reproduction number is = 1.0222 > 1. Therefore, the system converges to the 0 2 . Fig. 2 exhibits the solution of the system (1) to (5) = 0.4891, 0.4353,0.0621,0.00069721,0.0128) with the absence of delay by using the values of the parameters presented in Tab.2 and converges to true equilibria of the system, as desired. Let, . We can observe that, the effective use of delay tactics may decrease the value of = 1.022 threshold number and even though the given system switch from virus incidence equilibrium to virus absenteeism equilibrium. Thus, the absence of coronavirus is stable. But, Fig. 3 shows that, infection could be controlled due to the effective uses of different kind of delay tactics, as desired. Let now assume for different values of (delay term). We can see that, the infected step by step moves towards zero means infection is controlled. Therefore, Fig. 4 shows that the delay strategy or delay tactics such as vaccination, quarantine, travel restriction, holidays and distancing measure, etc. are a vital role to control the current strain of coronavirus. In the present study, we have analyzed the dynamics of the coronavirus model in the human population with the technique of mathematical modeling of the delayed term. We have categorized the whole population into five components of subpopulations. We have studied the system of nonlinear delayed differential equations for the dynamics of coronavirus in humans. We have proved its stability locally and globally by using well-known theorems. We have examined the validity of delay term on the reproduction number and the infected component of humans. The dynamics of the investigated system depends on the use of delay terms. Furthermore, delay strategies such as social distancing, placement of isolation, extension in holidays, restriction on travels are helpful to overcome the dreadful situation of coronavirus. In the future, our study can be extended in stochastic spatiotemporal and fraction order models, in the sense of delayed models [22] [23] . This type of modeling could also be studied with different fractal fractional derivatives as presented in [30] [31] . Also, numerical algorithms may be designed in the modeling of epidemic diseases using fuzzy differential equations and stochastic differential equations of fractional order [32] [33] [34] [35] . We have concluded that the analysis of delayed mathematical modeling (DMM) plays a significant role in the dynamics of pandemic models. Towards an artificial intelligence framework for data-driven prediction of coronavirus clinical severity Early dynamics of transmission and control of COVID -19: a mathematical modelling study A reliable numerical analysis for stochastic hepatitis B virus epidemic model with the migration effect A reliable stochastic numerical analysis for typhoid fever incorporating with protection against infection Competitive analysis for stochastic influenza model with constant vaccination strategy Emotion classification based on brain wave: a survey. Human-centric Computing and Information Sciences A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and government action COVID-19 infection: origin, transmission and characteristics of human coronaviruses Transmission potential and severity of COVID-19 in South Korea Stability behavior of mathematical model MERS corona virus spread in population. Filomat Modeling the epidemic dynamics and control of COVID-19 outbreak in China Mathematical modeling of the spread of the coronavirus disease 2019 (COVID-19) taking into account the undetected infections. The case of China A mathematical model for the novel coronavirus epidemic in Wuhan Mathematical modeling of epidemic diseases; a case study of the COVID-19 coronavirus A fractional-order model for the novel coronavirus (COVID-19) outbreak A fractional differential equation model for the COVID-19 transmission by using the Caputo-Fabrizio derivative Advances in Difference Equations Hiv and shifting epicenters for COVID-19, an alert for some countries A mathematical modeling with optimal control strategy of transmission of COVID-19 pandemic virus A model-based study on the dynamics of COVID-19: Prediction and control Mathematical analysis of novel coronavirus (2019-nCov) delay pandemic model Modeling the dynamics of novel coronavirus (219-nCov) with fractional derivative A reliable numerical analysis for stochastic gonorrhea epidemic model with treatment effect A reliable numerical analysis for stochastic dengue epidemic model with incubation period of virus An effective numerical method for the solution of a stochastic coronavirus (2019-ncovid) pandemic model. Computers, Materials and Continua Design of nonstandard computational method for stochastic susceptible-infected-treated-recovered dynamics of coronavirus model A: Facemasks simple but powerful weapons to protect against COVID-19 spread: can they have sides effects? Results in Physics Modelling the spread of covid-19 with new fractal-fractional operators: can the lockdown save mankind before vaccination? On the mathematical analysis of Ebola hemorrhagic fever: deathly infection disease in West African countries Modeling the effect of delay strategy on the transmission dynamics of HIV/AIDS disease A: Can transfer function and Bode diagram be obtained from Sumudu transform New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model S: A numerical algorithm for solving non-homogeneous fuzzy differential equations of fractional order D: Numerical solutions of fuzzy differential equations by an efficient Runge-Kutta method with generalized differentiability, Fuzzy Sets and Systems S: An operational matrix based on Legendre polynomials for solving fuzzy fractional-order differential equations Funding Statement: The authors declare that they have no funding for the present study.