Global stability of COVID-19 model involving the quarantine strategy and media coverage effects

In this paper, we build and analyze a mathematical model of COVID-19 transmission considering media coverage effects. Due to transmission characteristics of COVID-19, we can divided the population into five classes. The first class describes the susceptible individuals, the second class is exposed individuals, the third class is infected individuals, the fourth class is quarantine class and the last class is recovered individuals. The existence, uniqueness and boundedness of the solutions of the model are discussed. The basic reproduction number ℛ0 is obtained. All possible equilibrium points of the model are investigated and their local stability is discussed under some conditions. The disease-free equilibrium is local asymptotically stable when ℛ0<1 and unstable when ℛ0>1. The globally asymptotical stability of all point is verified by Lyapunov function. Finally, numerical simulations are carried out to confirm the analytical results and understand the effect of varying the parameters on spread of COVID-19. These findings suggested that media coverage can be considered as an effective way to mitigate the COVID-19 spreading.


Introduction
In December 2019, the world is facing the emergence of a new pandemic, which is called coronavirus disease 2019 . Then, COVID-19 spreads to world widely over the first two months in 2020. There were 492,510 confirmed cases of COVID-19 infection and 22,185 dead cases in world [1,2]. Therefore, it poses a continuing threat to human health because of its high transmission efficiency and serious infection consequences as well, it transmits by direct contact. Many researchers have tried to study and understand the dynamical behavior of COVID-19 through the transmission dynamics and calculate the basic reproduction number of COVID-19. It has become a key quantity to determine the spread of epidemics and control it. For example, in [3], Li et al. conducted a study of the first 425 confirmed cases in Wuhan, China, showing that the reproduction number of COVID-19 was 2.2, and revealed that person to person transmission occurred between close contacts. Other research [4] shows that the reproduction number of COVID-19 becomes 2.90, which is being increasing. In [5], Riou et al. studied pattern of early human to human transmission of COVID-19 in Wuhan, China. In [6], Hellewell et al. investigated the feasibility of controlling 2019-nCoV outbreaks by isolation of cases and contacts. Chen et al. [7], suggested mathematical model for simulation the phase-based transmissibility of novel coronavirus. Bentout et al. [8] developed an susceptible exposed infectious recovered model to estimation and prediction for COVID-19 in Algeria. Belgaid et al. [9] suggested and analysis of a model for Coronavirus spread. Owolabi et al. [10] proposed and analyzed a nonlinear epidemiological model for SARS CoV-2 virus with quarantine class. Flaxman et al. [11] suggested and estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Kennedy et al. [12] suggested a mathematical model involving the effects of intervention strategies on COVID-19 transmission dynamics. Feng et al. [13] studied a COVID-19 model with the effects of media and quarantine in UK. In this present study, we will show effects of the quarantine strategy and media reports on the spread of COVID- 19. We propose a mathematical model for COVID-19 transmission dynamics with the quarantine strategy and media effects. We start the model formulation by denoting the total size of the population by N which is classified further into five classes, the susceptible S (t), the exposed E (t) , the infected I (t) , the hospital quarantined Q (t) and the recovery R (t) at any time t, So, N = S + E + I + Q + R. The exposed class means low-level virus carrier, which is considered to be non infectious. The quarantined class in which the individual who is in the process in hospital, we suppose that only those who treat it will be in contact with the infected population. Accordingly, the flow of corona virus pandemic along with the above assumptions can be representing in the following block diagram: And the corresponding dynamical model has formulated through the nonlinear differential equations as follows, with initial conditions In model (1), the birth rate A is taken into susceptible class and natural death rate of population is given by the parameter d. The susceptible will be infected through sufficient direct contacts with infected people in the absence of media alerts by β 1 , with fraction parameter c, where c ∈ [0, 1] . The term β 2 I m+I reduce the transmission as media continuously alert the susceptible and exposed regarding infected cases and possible preventive measures. Usually, we assume that β 1 ≥ β 2 . As well, we consider the media awareness cannot stop the outbreak of COVID-19 but can aware the population to minimize the transmission risk through half saturation of media constant m. The death due to disease rate µ affecting from infected class only. k represent to a fraction denoting the level of exogenous re-infection. The quarantined rate is given by . And the mean recovery rates of class I, Q are γ i , i = 1, 2, respectively.
It is easy see that the 4 th and 5 th equations are a linear differential equation with respect to variables I(t) and R(t), which are not appear in the other equations of model (1). Hence model (1) can be reduced to the following model: In this paper, we will discuss the dynamics of model (3) with initial conditions This paper is organized as follows. In section 2, we will build the basic properties of model such as (positivity, boundedness of solutions and basic reproduction number). Existence of equilibrium points is presented in section 3. In section, the phenomenon of backward bifurcation is considered. The local and global stability of equilibrium points are studied in sections 4. In section 5, numerical simulation results are given. We conclude this paper with a brief conclusion.

Positivity of solutions
On the positivity of solutions for model (3), we have the following result. Theorem 2.1 Every solution of (3) with initial values (4) is positive as t > 0.
which can be written as thus, so that Similarly, it can be shown that E(t) > 0 and I(t) > 0 for all time t > 0. Hence all solutions of the model (3) remain positive for all non-negative initial conditions, as required.

Boundedness
Theorem 2.2 All solutions of model (1) which initiate in R 5 + are uniformly bounded.
Proof. Define the function N (t) = S (t) + E (t) + I (t) + Q (t) + R (t) and then take the time derivative of N (t) along the solution of model (1) gives dN dt ≤ A − LN. Then, dN dt + LN ≤ A, where L = min{d, d + µ} . Now, it is easy to verify that the solution of the above linear differential inequalities can be written as and N(t) ≤ A L for ∀t > 0. Thus all solutions are uniformly bounded and the proof is complete.

Basic reproduction number
It is easy to see that model (3) always has a disease-free equilibrium P 0 (S 0 , 0, 0), where S 0 = A d . We can calculate the reproduction number R 0 of model (3) by using the next-generation matrix method illustrated by van den Driessche and Watmough in [14].
Consequently, from Theorem 2 of [14], we have the following result. Theorem 2.3 The disease-free equilibrium P 0 of the model (3) is locally asymptotically when R 0 < 1 and P 0 is unstable when R 0 > 1.
The basic reproduction number for COVID-19 infection R 0 measures the average number of new COVID-19 infections generated by a single infected individual in a completely susceptible population [14,15]. Theorem 2.3 implies that COVID-19 can be eliminated from the community (when R 0 < 1) if the initial sizes of the sub-populations of the model (3) are in the attraction basin of the disease-free equilibrium P 0 . To ensure that COVID-19 elimination is independent of the initial sizes of the subpopulations, it is necessary to show that the disease-free equilibrium P 0 is globally asymptotically stable when R 0 < 1.

Existence the COVID-19 equilibria point and backward bifurcation
In this section, we consider the number of equilibrium solutions the model (3). To do so, let P * (S * , E * , I * ) be any arbitrary equilibrium of the model (3). Setting the right sides of the model (3) to zero gives here Since we assume β 1 > β 2 , S * and E * are positive. now, substituting (12) in 3 rd equation of the model (3) and simplifying it, we get where From (15), we can find that D 1 < 0. And D 5 > 0 when R 0 > 1, D 5 < 0 when R 0 < 1. Thus, the number of possible positive real roots the polynomial (12) can have depends on the signs of D 2 , D 3 and The various possibilities for the roots of f (x) can be analyzed using the Descartes Rule of Signs. The various possibilities for the roots of f (x) are tabulated in Table 1. Table 1. Number of possible positive real roots of equation (14).
Number of possible positive real roots has a unique endemic equilibrium if R 0 > 1 and whenever Cases 1, 5, 7 and 8 are satisfied; (ii) could have more than one endemic equilibrium if R 0 > 1 and Cases 2, 3, 4 and 6 are satisfied; (iii) could have 2 or more endemic equilibria if R 0 < 1 and Cases 1-7 are satisfied.
From the 4 th and 5 th equations of model (1) we can determent the values of Q * and R * through The existence of multiple endemic equilibria when R 0 < 1 suggests the possibility of backward bifurcation [16], where the stable disease-free equilibrium co-exists with a stable endemic equilibrium when R 0 < 1. This is can be obtained using Centre Manifold Theory.
Now let β 1 = β * 1 be the bifurcation parameter. Observe that at R 0 = 1, With β 1 = β * 1 the transformed model equation (17) has a simple eigenvalue with zero real part and all other eigenvalues are negative (that is has a hyperbolic equilibrium point). Thus, Center Manifold Theory can be applied to investigate the local dynamics of the transformed system (17) near β 1 = β * 1 . Now the Jacobian matrix of the transformed system evaluated at COVID-19 free equilibrium P 0 is obtained as It is easy to obtain the right eigenvectors of this Jacobian matrix as Similarly, it is possible to obtain the left eigenvectors which are denoted by W = (w 1 , w 2 , w 3 ) = (0, 0, 1). Now proceeding to obtain the bifurcation coefficients a and b as defined in Theorem 4.1 in [17].
First the non-vanishing partial derivatives of the transformed model (17) evaluated at COVID-19 free equilibrium are obtained as so that The sign of the bifurcation parameter b is associated with the following non-vanishing partial derivatives of F(X), also evaluated at the disease free equilibrium P 0 : The bifurcation coefficient b is obtained as Obviously, b is always positive. From Theorem 3.2 the system (17) will exhibit backward bifurcation phenomena if the bifurcation coefficient a is positive. The positivity of a in (22) gives the condition for backward bifurcation, which leads to

Stability analysis
In this section, the stability analysis of the all equilibrium points of model (3) studied as shown in the following theorems by used some criterion.
Theorem 4.1 The COVID-19 equilibrium point P * of the model (3) is locally asymptotically if the following conditions are hold Proof. The Jacobian matrix of model (3) at P * can be written as a 11 0 a 13 a 21 a 22 a 23 a 31 a 32 0 here clearly, the characteristics equation of J (P * ) is given by where furthermore, we have that Now, according to Routh-huewitz criterion P * will be locally asymptotically stable provided that B 1 > 0, B 3 > 0 and ∆ > 0. It is clear that if above conditions (26)-(27) hold.
The purpose of this section is to investigate the global stability by using Lyapunov function for COVID-19 free equilibrium point and COVID-19 equilibrium point respectively. We obtain the result in the following theorems Theorem 4.2 The disease-free equilibrium P 0 is globally asymptotically stable provided that the following condition holds: Proof. Consider the following function clearly, V 0 : R 3 + → R is a continuously differentiable function such that V 0 (S 0 , 0, 0) = 0 and V 0 (S , E, I) > 0, ∀ (S , E, I) (S 0 , 0, 0). Further, we have now, by doing some algebraic manipulation and using the condition (33), we get Obviously, dV 0 /dt = 0 at P 0 = (S , 0, 0), moreover dV 0 /dt < 0 otherwise. Hence dV 0 /dt is negative definite and then the solution starting from any initial point satisfy the condition (33), will approaches asymptotically to COVID-19 free equilibrium point. Hence the proof is complete.
Proof. At the COVID-19 equilibrium point P * = (S * , E * , I * ), S * , E * and I * satisfies the following equations By above equations (4.4) and assumptions now, define the Lyapunov function clearly, by derivative of V 1 we get furthermore, by simplifying the resulting terms, we get that Since the arithmetical mean is greater than, or equal to the geometrical mean, then 2 − x − 1 = 0 if and only if x=y=u=1; y − u + (u 2 − uy) m+I * m+I ≤ 0 for y, u > 0 and y − u + (u 2 − uy) m+I * m+I = 0 if and only if y = u = 1. Therefore,V 1 ≤ 0 for x, y, u > 0 andV 1 = 0 if and only if x = y = u = 1, the maximum invariant set of model (3) on the set (x, y, u) :V 1 = 0 is the singleton (1,1,1). Thus, the COVID-19 equilibrium point P * is globally asymptotically stable if R 0 > 1,by LaSalle Invariance Principle [18]. Hence, the proof is complete.

Numerical simulation
For the parameters values of model (1.1), we can chosen the parameters values from real data available sense Feb. 24 2020 to Apr. 5 2020. The total population of the Iraq for the year 2020 is approximately 40 × 10 6 [19]. The life expectancy in Iraq is approximatily 71.08 [19]. Clearly, we can obtain that the natural death rate d = 3.8545 × 10 −5 per day. The birth rate is estimated from A/d = N, and assumed that this is to be the bound population in the disease absence. So, A = 1541.8 per day and the other parameters of our model shows that in Table 2. We plot the solution trajectories of model (1) with initial point (15,20, 500, 1000, 150) which converges to COVID-19 equilibrium point P * = (1,27,2773,5428,19371), shown that in Figure 2.  Table 2. (1) close all borders with neighboring countries (2) to declare a state of emergency and impose a curfew (3) medical support from the government (4) methodological improvement on the diagnosis and treatment strategy (5) spontaneous household quarantine by citizens (6) more newly-hospitals put into use (7) massive online teaching in postponed semester (8) addition of new diagnosis method clinically diagnosis in Baghdad and some provinces 0.08 In the face of the COVID-19 outbreak, many stringent measures were taken by Iraqi government present in Table 3, to simulate the impact of different government control measures on the number of all S (t), E(t), I(t), Q(t) and R(t). We assumed that some values to contact rates with awareness Table 3, with the other parameters in Table 2 staying still on the all stages.
The following Figure 3 shows the values of S (t), E(t), I(t), Q(t) and R(t) under government measures that above to control of COVID-19 outbreak.  Table 2.
Clearly, from above figure for effect of contact rate Table 3, We obtain that in case decrease the contact rate (social isolation) the reproduction number less than one and the dynamical behavior of model (1.1) still approaches to COVID-19 equilibrium point. Hence, the backward bifurcation is occur. Now, to investigate the effect of the quarantined strategy it is given by on the dynamical behavior of model (1.1) and to control to COVID-19 outbreak in Iraq. We study the impact of this parameter on values of S (t), E(t), I(t), Q(t) and R(t) in follows Table 4 and shows the results in Figure 4. (1) to declare a state of emergency and impose a curfew to reduce the contact between people (2) medical support from the government (3) methodological improvement on the diagnosis and treatment strategy (4) spontaneous household quarantine by citizens (5) addition of new diagnosis method clinically diagnosis in Baghdad and some provinces 4.5 The following Figure 4 shows the values of S (t), E(t), I(t), Q(t) and R(t) under government measures that above to control of COVID-19 outbreak.  Table 2.
Clearly, from above investigate to impact of the quarantined strategy Table 4, when the quarantine strategy increasing we get the number of infected is decrease and other classes are increase. Here, we ask whether the quarantine strategy is the best solution? The answer is possible, but for specific numbers. Whereas, if the quarantine is more than the capacity of the health institutions. We get the dynamical behavior of model (1.1) lose the stability as shown in Figure 5.  Table 2.

Discussion and results
In this research, a mathematical model of COVID-19 transmission has been proposed by compartment the total population into five epidemiological status, namely, susceptible S (t), exposed E(t), infected I(t), quarantine Q(t) and recovered R(t). The model incorporates the impact of social awareness programs conducted by public health officials with quarantine strategy in hospital. It has been noticed that these awareness programs and quarantine strategy result in human behavioral changes in order to avoid risk of disease transmission. The model mainly accounts for the reduction in disease class due to awareness. While we can say the disease goes away due to applied the quarantine it well. The proposed model has two biological equilibrium points are COVID-19 free and COVID-19. The COVID-19 free has been local stability when R 0 < 1. Otherwise when R 0 > 1, the COVID-19 free point becomes unstable and the dynamical behavior of the model converges to COVID-19 equilibruim point. The backward bifurcation occur if R 0 = 1 at the parameter bifurcation β 1 = β * 1 = d ( + γ 1 + d + µ) /cA. As well as the different government control measures have been also discussed. Furthermore, to shown and understand the effect of quarantine rate of disease we have choosed many different value of it say parameter then we have obtained some different results see Table 4 and Figure 4.