key: cord-0861196-ai8o8raf
authors: Verma, Harendra; Mishra, Vishnu Narayan; Mathur, Pankaj
title: Effectiveness of lock down to curtail the spread of corona virus: A mathematical model()
date: 2021-01-20
journal: ISA Trans
DOI: 10.1016/j.isatra.2021.01.033
sha: 262a9d20f8c56346b8b97850bcbad4dc658ee7c6
doc_id: 861196
cord_uid: ai8o8raf

In this paper, we have considered a mathematical model that deals with the effectiveness of the measures that may be helpful for reducing the spread of the COVID-19 virus in the society. Here we have illustrated the importance of lock down in controlling and maintaining the spread of the COVID-19 virus. The impact of the virus on the susceptible population has been considered in the model. Also, we have taken into account the susceptible population, which by taking preventive measures viz., by having strong immunity, maintaining social distancing, wearing PPE kits and masks etc., is able to reduce the possibility of getting infected from the virus. Local as well as global stability of the equilibrium points of the model have been studied using Lyapunov function and the geometrical approach techniques. Basic reproduction number has also been obtained by using the next generation matrix. To show the effectiveness of the model, different cases obtained by varying the parameters involved in the model have been considered. A comparison between the actual number of infected cases in India and that obtained by the proposed model, showing the effectiveness of the proposed model, has also been carried out.

Corona virus disease [1] is an infectious disease caused by a newly discovered corona virus. The virus that causes COVID-19 is mainly transmitted through droplets generated when an infected person coughs, sneezes, or exhales. These droplets are too heavy to hang in the air and quickly fall on floors or surfaces. One can get infected by breathing in the virus if he is within close proximity of someone who has COVID-19 or by touching a contaminated surface. People who are infected with COVID-19 disease experience mild, moderate and sometimes no symptoms.

Even though there is no specific treatment for the disease caused by this virus, the doctors, nurses and paramedical staff are putting their best effort to cure and save life of the patient. In doing so they are in a very close proximity with the infected patient and hence are most vulnerable of getting infected. But by wearing PPE kit, masks, gloves etc. and taking proper care they try to protect themselves from the infection. But still, many medical staff are getting infected either because the viral load in the patient is very high or due to their carelessness or not using safety kits.

being has to go out for his earnings the social distancing becomes difficult to follow, whereas the immunity varies from person to person.

As such the only way to minimize the spread of the disease is not only avoiding the interaction of the uninfected person with the infected persons and the persons to whom they have met but also stopping the susceptible persons to visit the places where the infected persons have visited. This is possible only when the movement of the uninfected as well as the population which has somehow come in contact (directly or indirectly) with the infected person is restricted.

To reduce the spread of disease government of many countries took a hard decision to lockdown their countries. Even though lock down is not an easy step for any government as it serves a major disruption of the economy, still no nation was in the position to take the risk of allowing the epidemic to spread and then treating a huge number of patients. The information regarding the possible ways of protection against the disease is propagated through TV, radio, social media etc. As a result, some of the aware people use precautionary measures for their protection.

A lot of mathematical models dealing with different aspects of COVID-19 have come up recently viz., [2] , [3] , [4] , [7] - [10] and references therein. Most of them have considered the interaction of infected population with the susceptible population as a major cause of spread of disease. But, this disease does not only spreads through direct contact between the susceptible to the infected individuals, whereas, these infected individuals act as carriers of the virus and may spread the disease. As such, we have assumed that the disease spreads among the susceptible individuals through transmission of virus that is present not only in the infected individuals but also in the environment.

In this paper, we have considered a mathematical model that deals with the importance of lock down [4] in controlling and maintaining the spread of the COVID-19 virus. We have considered the impact of the virus on the susceptible population instead of the interaction of infected population with the susceptible population. This results in Herds transmission instead of one to one transmission, which is the case in the spread of COVID-19 virus. Also, we have taken into account the fact that a portion of the susceptible population have a lesser possibility of getting infected from the virus due to various reasons like strong immunity, taking proper measures of protection like maintaining social distancing, wearing PPE kits, masks etc..

Taking the preventive measure against the pandemic situation created worldwide due to COVID-19, government of India announced a complete lock down from 25th March, 2020 to protect its citizens from getting infected from the deadly virus. It not only affected the life of common people but also the economic activities of the country very badly. So as to overcome the economic and social problems, the lock down was relaxed in different phases starting from 14th April, 2020 and has not been removed completely yet. With the help of this mathematical model, we are interested in studying the impact of virus during different phases of the lock down in India. We find the error between the number of predicted infected cases with the help of the proposed model and the actual infected cases as per the data provided by Indian Council of Medical Research (ICMR) during the lock down period in India. Further, as this error is under tolerable limits we are able to predict the growth of the virus and when it would be reaching its peak.

Local as well as global stability of the equilibrium points of the model have been studied using Lyapunov function and the geometrical approach techniques. Basic reproduction number, indicating the spread of the disease, has also been obtained by using the next generation matrix. The parameters of the model are firstly approximated by using the best fit least squares approximation technique. Secondly, we have considered different cases by varying the parameters and in each case obtained the reproduction number (R 0 ) showing whether the spread of virus is controlled or not. It turns out that a few parameters play a vital role to maintain the effectiveness of the lock down. Also, a comparative study, between the actual number of infected cases in India as reported by ICMR [26] and that obtained by the proposed model, has been carried out. Both the values are comparable. Lastly, with the help of the proposed model, we have obtained a figure predicting the total number of infected cases, number of days for the virus to spread, reach to its peak and its eradication, in India.

J o u r n a l P r e -p r o o f

We consider that, in the region under consideration the total population is N(t) at time t. The total population is divided into three classes; the susceptible population S (t), the infective population I(t) and recovered population R(t). Also v(t) be free virus present in infected person at any time t. Here it is assumed that the disease COVID-19 spreads through transmission of virus [19] to susceptible individuals by infected person or through virus spread by the infected person by coughing, sneezing, exhaling etc. at the places where he has visited. Growth rate of virus is assumed to be proportional to the number of infected individuals. The immigrants coming from different nations or places where the virus has already spread are more eligible virus carriers. They, even though are infected, stay with susceptible class of population knowingly or unknowingly. Also, a lot of susceptible individuals, to whom these infected individuals contact to, become infected and the chain continues due to which the number of infected individual grow exponentially.

The above discussed model can be represented by the following system of non-linear ordinary differential equations:-

In the above model (2.1), let A be the number of immigrated people in the susceptible class and µ be the rate at which the immigrated population in susceptible class becomes infected, β be the rate at which a person in susceptible class moves to infected class, n p (0 < n p < 1) is the rate at which susceptible class protect themselves from moving to the infected class by taking proper precautions like wearing protective PPE kit or mask, maintaining social distancing, improving their immunity etc., a is the rate at which the infected persons move together with the susceptible persons without taking any preventive measures like putting a mask or maintaining the social distancing. In other words, we can say that a is the measure of how much the lock down has been followed by the public. Lower the value of a higher is the effectiveness of lock down and a = 0 imply complete lock-down strictly followed by the public. δ is the natural rate of death, α is the disease induced death rate, γ is the recovery rate, ac 1 is the self growth rate of the virus due to increase in the infected individuals, c 2 is the natural death rate of the virus and c 3 is its growth rate, for the feasibility [5] of the model it is assumed that c 2 − c 3 > 0. Using the fact that N = S + I + R, the above system reduces to the following system,

(2.2)

It is important that the model system (2.1) is epidemiologically feasible [15] , [16] . We will show that all the variables involved in model system (2.1) are non-negative with positive initial conditions for all time t > 0 as proved in the following theorem: Proof 1. We will show that all the variables of model (2.1) are non-negative for all time t. It follows from the second equation of (2.1) that dI dt ≥ µA − (a + δ + α + γ)I, the above equation can be written as

Hence

µA exp(a + δ + α + γ)s ds,

This shows that I(t) > 0 for all t > 0. Further, from the first equation of (2.1), we have

hence we obtain that Proof 2. The system (2.2) is split into two parts, the human population (i.e., S (t), I(t), R(t)) and free virus (i.e., v(t)). It follows from the first three equations of the system (2.2) that

Multiply the integrating factor e δt of the above equation and integrating from 0 to t, we obtain

From the fourth equation of the model (2.2), and using the fact that

Proceeding in a similar way as done for (3.1), we obtain lim

for all t > 0. From the above discussion, we can see the feasible region of the human population for system (2.2) is

and feasible region of virus population for system (2.2) is

which proves the theorem. 2

Define Ω = Ω H × Ω v . Now, Ω is a positive invariant for the system (2.2). Let Int(Ω) denote the interior of Ω.

The above system (2.2) have three non-negative equilibrium points as follows:

This equilibrium point is obtained by equating to zero the right hand side of the system of equation (2.2) together with I 0 = 0, R 0 = 0 and v 0 = 0, which gives N 0 = A δ .

Here, we obtain the basic reproduction number of model system (2.2) , at E 0 , using next generation matrix approach [12] , [20] . For our model (2.2) F and V are given as:

Let F and G be the Jacobian's of F and V respectively at disease free equilibrium point which are given as:

Now the basic reproduction number R 0 [10] of the model is the spectral radius of the matrix FG −1 , which is

Special Case: Since an increment in the recovered population R(t) depends only on recovery rate γ, thus, in particular, if γ = 0, then

It is easy to observe from the above cases and expressions that R 01 is bigger than R 0 .

Equating to zero the right hand side of the system of equation (2.2) together with R # = 0, which gives

From equation (4.2) and (4.3),

Putting the above value N # and v # in (4.1), we get

By Descartes' rule of signs at least one positive root must exist.

Equating to zero the right hand side of the system of equation (2.2), which gives

and

From equation (4.5), (4.6) and (4.7), we get

Putting the above value R * , N * and v * in equation (4.4), we get In this section we discuss the stability of the model at the three equilibrium points obtained in the preceding section.

In this section, we obtain the local and global stability conditions for disease free equilibrium point E 0 given in section (4.1).

Theorem 3. The disease free equilibrium point E 0 , given in section (4.1), is locally asymptotically stable for R 0 ≤ 1.

Clearly there are repeated roots at λ = −δ < 0. Now by Routh-Hurwtz criterion, the other two roots of the characteristic equation (5.1) will be negative if both P 0 and P 1 are positive. Since c 2 > c 3 thus P 0 > 0 and P 1 > 0 when R 0 < 1. Thus we conclude that the disease free equilibrium point E 0 of the model system (2.2) is locally asymptotically

Theorem 4. The disease-free equilibrium E 0 is globally asymptotically stable for R 0 ≤ 1.

Proof 4. Consider the following Lyapunov function:

On differentiating the above equation w.r.t. t and using (2.2), we get, [14] , we observe that the system (2.2) has the maximum invariant set for L = 0 if and only if R 0 ≤ 1 and I = v = 0. The epidemiological implication of Theorem 4 is that the infected fraction (i.e., I and v) goes to zero in time when R 0 ≤ 1, that is, the COVID-19 eventually disappears from the population. The disease is endemic if the infected fraction remains above a certain positive level for sufficiently large time [23] . System (2.2) can be defined to be uniformly persistent if

for some > 0 for all initial points in int(Ω).

Theorem 5. The model system (2.2) is uniformly persistent in int(Ω) if R 0 > 1.

Proof 5. Ω is positively invariant with respect to the model. The disease free equilibrium is unique and located on the boundary of Ω. On using the result stated in theorem (4.3) in [13] , the unstability of the disease free equilibrium implies the uniform persistence of the model system (2.2). 2

In this section, first, we obtain the condition for the local stability of the endemic equilibrium point E 1 using the Routh-Hurwitz criterian followed by global stability of the point using the Lyapunov's method [6] , [11] .

Local stability of non-removal class equilibrium point Theorem 6. If R 01 > 1 then the endemic equilibrium point is locally asymptotically stable provided ρ i > 0, where i = 1, 2, 3 and ρ 1 ρ 2 − ρ 3 > 0.

Therefore the corresponding characteristic equation is

where λ is the eigen value with the first eigen value λ = −δ and the other eigen values are given by the equation

Since ρ 1 , ρ 2 and ρ 3 are positive and by algebraic calculations we have that ρ 1 ρ 2 − ρ 3 > 0, thus by Routh-Hurwitz criterion all roots of (5.2) are either negative or have negative real parts. Therefore equilibrium point is locally asymptotically stable. 

Proof 7. Consider a positive definite function U such that

where p 1 and p 2 are positive constants and i, n and v are small perturbation from the equilibrium point E 2 , that is,

dU dt will be negative definite provided, for k > 1, the inequalities

hold. On choosing 9 J o u r n a l P r e -p r o o f

the above inequalities may be combined from which the theorem follows. 2

In this section, firstly we obtain the condition for the local stability of the endemic equilibrium point E 2 using the Routh-Hurwitz criterion followed by global stability of the point using geometric approach [18] .

Local stability of endemic equilibrium Theorem 8. If R 0 > 1 then the endemic equilibrium point is locally asymptotically stable provided σ i > 0, where i = 1, 2, 3 and σ 1 σ 2 − σ 3 > 0. Therefore corresponding characteristic equation is

where λ is the eigen value and the one of the eigen value is λ = −δ, and other eigen values are given by

Since σ 1 , σ 2 and σ 3 are positive and algebraic manipulation convey that σ 1 σ 2 −σ 3 > 0, thus by Routh-Hurwitz criterion all roots of (5.3) are either negative or have negative real parts. Therefore equilibrium point is locally asymptotically stable. Global Stability of endemic equilibrium E 2 (I * , R * , N * , v * ) In the following, will apply the second compound matrix techniques and autonomous convergence theorems developed by Smith [22] , LiandMuldowney [17] to demonstrate the global stability of E 2 .

Theorem 9. The unique endemic equilibrium is globally stable if R 0 > 1 provided that the following inequality:

Proof 9. The Jacobian matrix J of model system (2.2) at an arbitrary point (I, R, N, v) is given as

The second compound matrix J [2] is

and

J o u r n a l P r e -p r o o f Next, using the second compound matrix J [2] and Q f is the directional derivative of Q in the direction of the vector field f , we calculate the following [24] , [21] ,

which is obtained as

Since, we have

where,

Now, we define the following norm on R 6

where z ∈ R 6 , with components z i (i = 1, 2, 3, 4, 5, 6) and Case 1. U 1 (z) > U 2 (z) and z 1 , z 2 , z 3 > 0 In this case, z = max{|z 1 |, |z 2 | + |z 3 |}.

Subcase 1.1: |z 1 | > |z 2 | + |z 3 |. Then z = |z 1 | = z 1 and U 2 (z) < |z 1 |. Taking the right hand derivative of z , we obtain

Since |z 3 | < |z 1 | and |z 5 | ≤ U 2 (z) < |z 1 |, we have

Then z = |z 2 | + |z 3 | = z 2 + z 3 and U 2 (z) < |z 2 | + |z 3 |. On taking the right hand derivative of z , we obtain

J o u r n a l P r e -p r o o f

Case 2. U 1 (z) > U 2 (z) and z 1 < 0 < z 2 , z 3 In this case z = max{|z 1 | + |z 3 |, |z 2 | + |z 3 |}.

Subcase 2.1: |z 1 | > |z 2 |.

Then z = |z 1 | + |z 3 | = −z 1 + z 3 and U 2 (z) < |z 1 | + |z 3 |. On taking the right hand derivative of z , we obtain

J o u r n a l P r e -p r o o f

Subcase 2.2: |z 1 | < |z 2 |.

Then z = |z 2 | + |z 3 | = z 2 + z 3 and U 2 (z) < |z 2 | + |z 3 |. On taking the right hand derivative of z , we obtain

The remaining other 12 cases may also be dealt with in a similar way, we omit details. Combining the results of the four case presented here in equations (5.4),(5.5), (5.6) and (5.7) as well as the remaining 12 cases, we get

Therefore,μ(B) < 0 on int(Ω), hence we can obtain that the endemic equilibrium E 2 is globally asymptotically stable. 2

To check the feasibility and effectiveness of the model proposed in (2.1), we perform numerical simulation. Different values of the parameters involved have been estimated with the help of best fit least squares approximation and the reproduction number has been calculated in each case. The estimated values of different parameters of the proposed model are: In the model (2.1), the parameters a, β and n p play a vital role in determining the effectiveness of the lockdown in the spread of the virus. As such, we consider a few cases, where by varying only these parameters and keeping the other parameters fixed, the value of the reproduction number R 0 becomes less than 1, that is, the spread of the disease decreases. Precisely, we fix parameters as A, δ, α, γ, c 1 , c 2 , c 3 , µ and vary a, β and n p as per the following: (1), the reproduction number is 3.34989, which means that the spread of disease is endemic and that is happening in the present situation. 2. In Case (2), n p is increased, that is, when the susceptible class becomes aware and takes some precautions by wearing protective kits, washing hands or increasing immunity and the values of other parameters are kept fixed as in Case (1), then the reproduction number decreases, but spread does not break although it slows down. 3. In Case (3), a is decreased, that is, the rate at which the infected individuals move together with susceptible individuals without maintaining social distancing or taking any preventive measures, is decreased and values of other parameters are same as in Case (2), then reproduction number decreases and becomes < 1, stating that the spread of virus is decreased. 4. Further in Case (4), when β, that is, the transmission rate of the virus is increased and values of other parameters are kept same as in Case (2), then the reproduction number decreases considerably and spread of virus deceases. 5. In the Case (5), if the value of a is further decreased and values of other parameters remain same as in Case

(1), then the reproduction number decreases and becomes less than 1, which implies that the spread of the virus gets downward trend and the infection gets reduced.

In this section, a comparative study of values of the variables obtained from the model (2.1) by suitable choice of parameters with the real data has been achieved with the help of figures. The data of the daily new COVID-19 cases in India has been compared with that obtained by the model (2.1) with proper choice of the parameters involved. The real data regarding the actual number of COVID-19 cases in India have been collected from the official website of ICMR ( [26] ) and World Health Organization (WHO)( [25] ). Here we have considered the data of 68 days with effect from (w.e.f.) 25th March, 2020 to 31st May, 2020 and divide it into 4 phases. Phase-1 w.e.f. 25th March, 2020 to 14th April, 2020 (21 days), Phase 2 w.e.f. 15th April, 2020 to 3rd May, 2020 (19 days), Phase 3 w.e.f. 4th May, 2020 to 17th May, 2020 (14 days) and Phase 4 w.e.f. 18th May, 2020 to 31st May, 2020 (14 days).

During Phase 1, it was a strict lock down and only persons in essential services were allowed to move, due to which, the number of daily reported infected cases were quite low. But to prevent the downfall of the economy it was decided by the government of India to relax the lock down in phases. In Fig.1 cases is quite close to that of actual ones.

The lock down was further relaxed by allowing more and more activities like opening up of offices, shops etc. This allowed the more public to move freely which increased the chances of spread of the virus. Figures 2, 3 and 4 reflect situation of the actual and predicted infected cases during the Phase 2 (15th April, 2020 to 3rd May, 2020), Phase 3 (4th May, 2020 to 17th May, 2020) and Phase 4(18th May, 2020 to 31st May, 2020) respectively. As is seen from the figures the number of infected cases get increased in much higher proportion than that in earlier Phases and an approximate number of infected cases as calculated by the model (2.1) is quite near to the actual number of cases. Daily Reported New Infective Cases

A cumulative affect of the nation wide lock down (complete as well as partial) w.e.f. 25th March, 2020 to 31st May, 2020 (68 days) has been presented in Fig.5 , showing that the error in predictions of the number of infected cases by model (2.1) to that of actual cases during the lock down period is comparable. 

To predict the effectiveness of the lock down, we have plotted above the reproduction number against the variation of the parameters a and n p . Fig.6 shows that as the value of a is increased, that is, as the lock down is relaxed the reproduction number increases which implies that the spread of the virus is pandemic. Fig.7 shows that with the increase in n p , that is, if the self immunity is high and protective measures like social distancing, wearing masks, washing hands etc. are followed, then the reproduction number decreases, that is, the spread of virus decreases. Also, on varying β, that is, if the rate of transmission is increased the reproduction number increases, indicating the increase in the spread of virus as shown in Fig.8 . 

Assuming the situation that the lock down be relaxed further, that is, more and more susceptible population 19 J o u r n a l P r e -p r o o f move out of their residence freely. The possibility of these people coming in contact with the symptomatic or asymptomatic infected persons increases and hence getting infected with COVID-19 virus also increases. With the suitable choice of parameters, using the model (2.1), it is shown in Fig.9 that the spread of COVID-19 virus will come to an end in India by itself but may take time approximately 2 years. As per model, the peak value of the total number of infected population in India is approximately 1 × 10 7 and is expected to reach by the end of January, 2021. If some medication of a vaccine is developed in the mean time, then the spreading of COVID-19 may end up early.

Epidemic analysis of COVID-19 in China by dynamical modeling

Assessment of 21 days lock down effect in some states and overall India: A predictive mathematical study on

Effective lock down and role of hospital-based COVID-19 transmission in some Indian states: An outbreak risk analysis

Modeling the impact of sanitation and awareness on the spread of infectious diseases

Effects of awareness programs by media on the epidemic outbreaks; A mathematical model

Why are hospitals hotbeds of COVID-19 transmission

Epidemic doubling time of the 2019 novel corona virus outbreak by province in mainland China

Epidemiological and clinical features of the 2019 novel corona virus outbreak in China

Preliminary estimation of the basic reproduction number of novel corona virus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak

Modeling and analysis of media's role in popularizing the blue whale game

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Uniform persistence and flows near a closed positively invariant set

Ordinary differential equations

Mathematical epidemiology of infectious diseases, model building, analysis and interpretation

Three basic epidemiological models

Smith's autonomous convergence theorem

A geometric approach to global-stability problem

Virus dynamics

An introduction to the basic reproduction number in mathematical epidemiology

Dynamics of a waterborne pathogen model under the influence of environmental pollution

Some applications of Hausdorff dimension inequalities for ordinary differential equations

Persistence under relaxed point-dissipativity (with application to an endemic model)

Global results for a cholera model with imperfect vaccination

Corona virus disease (COVID-19) outbreak

IJAA

Authors are thankful to the anonymous Referees for their valuable comments which helped a lot in improving the quality of the paper.

The authors declare that they have no conflict of interest in work reported in this paper.