{'fraction_non_alphanumeric': 0.03582659808963997, 'fraction_numerical': 0.019191770756796473, 'mean_word_length': 4.416427889207259, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 8, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 22, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mathematical modelling of COVID-19: A case study of Italy

This manuscript describes a mathematical epidemiological model of COVID-19 to investigate the dynamics of this pandemic disease and we have fitted this model to the current COVID-19 cases in Italy. We have obtained the basic reproduction number which plays a crucial role on the stability of disease free equilibrium point. Backward bifurcation with respect to the cure rate of treatment occurs conditionally. It is clear from the sensitivity analysis that the developments of self immunities with proper maintaining of social distancing of the exposed and asymptomatic individuals play key role for controlling the disease. We have validated the model by considering the COVID-19 cases of Italy and the future situations of epidemicity in Italy have been predicted from the model. We have estimated the basic reproduction number for the COVID-19 outbreak in Italy and effective reproduction number has also been studied. Finally, an optimal control model has been formulated and solved to realize the positive impacts of adapting lock down by many countries for maintaining social distancing.


Introduction
In the history of worldwide disease spreading public faced threat at different times due to appearance of new diseases in the human society, among them some are infectious diseases and some others are vector born diseases. The most notable threatening disease spreading were the Black death in Europe, small pox in Mexico, SARS in China in [2002][2003], MERS in Saudi Arabia in 2012, AIDS, Cancer, Malaria, etc. [7,12,24]. In 2019-2020 the most threatening disease is the coronavirus disease (COVID- 19), which is caused by novel coronavirus (SARS-CoV-2) [30]. As on 29th March, 2020 the COVID-19 spreads throughout the world [26,35,36] and about 199 countries have been affected. As a result WHO declared this disease as the Pandemic. The COVID-19 first was identified in Wuhan, China. Throughout the country up to 29th March, 2020 about 81 439, 3300 people were infected, died, respectively and throughout the whole world the numbers are 679 079 and 31 772, respectively [37]. The coronavirus infections previously occurred in different countries with low degree of harmfulness: those were the SARS-Cov (Severe Acute Respiratory Syndrome) [10] and MERS-CoV (Middle East Respiratory Syndrome) [34]. Still the outbreak of coronavirus is ongoing with higher degree of harmfulness as a result huge number of people are infected throughout the world. It creates huge threat to the global public health as well as the financial system also [6,23,38]. Though the infection was first reported in Wuhan, China in December 2019 but due to human movements it spreads throughout the world [22]. The emergence of COVID-19 occurs in a crucial time because this time is the spring festival time in Asia, so a large number of people move from one place to another places [23] and the disease starts to spread to huge number of people who move from their working place to their home land as a consequence the disease spreads to the whole world.
The infection of this disease spreads among the people through interactions and the intensity of the spreading is very high. So a large number of people is infected in a short span of time. Due to shortfall of medical facility and not proper medicine a large number of people in different countries is not getting proper treatment as a result large number of infected people is dying. To save from this disease the only way is maintaining social distancing from the infected persons and sterilizing whole body specially hand when a person goes outside the home. Hence the minimization of people movements is only way to control to spread the disease. To control the spreading of disease Chinese government adopted lock down policy to minimize the people interactions on 23rd January, 2020 as a result all types of public traffic were suspended in the infected Wuhan and Hubei province [23]. Following the strategies of China the highly affected countries are adopting the same policy.
The phenomenon of backward bifurcation has serious consequences for controlling the disease. For backward bifurcation, the stable disease free equilibrium point co-exists with a stable endemic equilibrium when the basic reproduction number is less than unity. The epidemiological consequence of the phenomenon of backward bifurcation is that having the basic reproduction number less than unity is, although necessary, no longer sufficient for eliminating disease from the community. In recent years, the phenomenon of the backward bifurcation has arisen the interests in disease control (see [1,2,13,15,19,25,39]). In [39], the authors studied the backward bifurcation analysis of an epidemic model with saturated treatment. In [15], Garba et al. presented a deterministic model for the transmission dynamics of a single strain Dengue disease, which exhibits the phenomenon of backward bifurcation. In [1,19], the authors described the phenomenon of backward bifurcation in the deterministic models for the transmission dynamics of Zika virus. Motivated by the above discussion, here we have studied the phenomenon of backward bifurcation in the epidemiological model of COVID-19, which signifies that the making the basic reproduction number less than unity is not enough to eradicate the COVID-19 disease from the community.
To study the COVID-19 disease dynamics and predict the presumed abundance of infected people in Italy, here we have formulated a five dimensional COVID-19 epidemic model. Here, we have divided total population into five classes, namely susceptible, exposed and asymptomatic, quarantine, infectious and recovered class. We have studied the condition of stability of the disease free equilibrium point (DFE) and the condition for backward bifurcation considering the cure rate of treatment as the bifurcation parameter. Also, we have fitted the model to the data on COVID-19 cases in Italy and we shall make some predictions on the epidemic in Italy. We have also reformulated the model as an optimal control problem by introducing the lock down controls adapted by the government of different countries to maintain the social distancing for controlling the epidemic. The effect of lock down will be studied in our present work and we shall investigate the effect of different model parameters, which can minimize the spreading of the disease.
Organization of the manuscript is as follows: In Section 2, we have formulated the model. The positivity and boundedness of solutions are discussed in Section 3. The evaluation of basic reproduction number and the DFE are given in Section 4. Section 5 is devoted for backward bifurcation and in Section 6, we validate the model with real data and estimate the model parameters. The estimation of the basic reproduction number for actual COVID-19 epidemics in Italy and the study of effective reproduction number have been given in Section 7. In Section 8, we have reformulated the model as an optimal control problem and finally in Section 9 some conclusions are summarized.

Model formulation
In formulating the model, we remember that the susceptible people (S) may be infected when they interact with the novel corona virus infected person. The COVID-19 infected person may not being infectious at the time of infection. The symptoms of COVID-19 may appear in as few as 2 days or as long as 14 days, during this incubation period the virus is contagious but the patient does not display any symptoms. So, the patient in that incubation period is able to transmit the virus to susceptible. Therefore, the asymptomatic transmission occurs in COVID-19 infections i.e. SARS-CoV-2 transmissions occur in the community from persons without symptoms [3,11,20,32]. Transmission by persons who are infected but do not have any symptoms can arise from two different infection states: presymptomatic individuals (who are infectious before developing symptoms) and individuals who never experience symptoms (asymptomatic infections, which can be called as never symptomatic) [20]. The presymptomatic persons and the never symptomatic persons both go to infectious class I after identification by the COVID test organized by the suitable authorities. A large number of cases in which the infected persons become infectious but asymptomatic and when susceptible individuals interact with them then infection spreads. The exposed individuals (infected but not infectious) and asymptomatic individuals are jointly included into the class E. The persons having symptoms like fever, cough, and shortness of breath are in quarantine and medical persons investigate them and after proper diagnosis COVID-19 positive patients are isolated. The individuals who stay in quarantine to maintain the social distancing from others have been included in the class Q and those who are diagnosed as COVID-19 positive patients have been included into the class I . Two important points we have to take into consideration. One is the delay in diagnosis of the symptomatic patients for conformation of COVID-19 positive patient and due to this fact some non-diagnosed persons in the Q class transmit the disease among the susceptible persons [17,18,23,26]. Second is the limitation of medical facilities when a large number of people becomes infectious. Since till date there is no proper treatment to recovery from COVID-19, only way is the development of self immunity. For this reason it is essential to provide the support system to the infectious persons at the critical situations. But if a large number of people is being infectious in very short span of time then medical support system suppliers will not be able to provide proper support to all the infectious persons. To include this effect we have adopted the saturated treatment rate to the infectious populations in the form a I (t) 1+bI (t) , where a is the cure rate of treatment and b is the delay in treatment of the infectious persons. Thus, at time t the total human population (N (t))is divided into five class, namely, the susceptible S(t), exposed and asymptomatic E(t), quarantined (Q(t)), infectious I (t) and the recovered R(t) class i.e. N (t) = S(t) + E(t) + Q(t) + I (t) + R(t). We assume that S(t) becomes infected with the interaction of the E(t), Q(t), I (t) classes with the rate of β 1 S(t)I (t)+β 2 S(t)E(t)+β 3 S(t)Q(t) , where the coefficients β 1 , β 2 and β 3 denote the transmission rates of infections from the corresponding classes. Here we also assume that the new recruitment of the populations occurs from the susceptible class only. The flow diagram of the COVID-19 disease dynamics in the population has been displayed in Fig. 1 and the corresponding mathematical formulation of the transmission dynamics of COVID-19 is given below: with initial conditions S(0) > 0, E(0) ≥ 0, Q(0) ≥ 0, I (0) > 0, R(0) ≥ 0. All the model parameters have been interpreted in Table 1. In the Appendix, we have discussed the case for other expressions for treatment rate and analysed the model by introducing an extra compartment for deaths.

Positivity and boundedness of solutions
Here, we shall show the positivity and boundedness of solutions of the model (1), which implies that the proposed model is biologically valid.
is a positively invariant set for the system (1).
Proof. The 1 st component of the solution of the system (1) can be approximated by using the following inequality Integrating the above and using the initial condition, we have Similarly, from the remaining four equations of (1) the following results can be obtained.
By using initial conditions we can conclude that all the solutions are feasible.
Next, adding all the equations of the system (1) we get Integrating, using the initial conditions and taking lim sup as t → ∞ we have lim sup t→∞ N (t) ≤ A µ . Therefore, the set Ω is the positively invariant and attracting for the system (1). Hence the theorem is proved. □ Thus, we have established that all components of the solutions of the system (1) are non-negative and ultimately bounded, which indicates that the model is biologically well behaved.

The basic reproduction number and the DFE
The basic reproduction number plays an important role in epidemiology. It indicates the transmissibility of the COVID-19 in the population. How rapidly the disease spreads is dependent on the value of the basic reproduction number R 0 , which is the number of secondary infections one COVID-19 positive person will produce in the population. Clearly, the system (1) have the disease free equilibrium point ( at which all the disease components are zeros. Now, we shall evaluate the basic reproduction number for the system (1) through next-generation approach [14].
Then Routh-Hurwitz criterion for polynomials implies that DFE E 0 is locally asymptotically stable.
Hence the theorem is proved. □

Backward bifurcation
In this section, we shall establish that the phenomenon of backward bifurcation takes place in the model (1) by considering the cure rate of treatment a as the bifurcation parameter. The method developed by Castillo-Chavez and Song [5,8] will be applied to determine the direction of the bifurcation at the critical value of the cure rate of treatment.
Proof. Let us rewrite the variables in the new form as S = x 1 , E = x 2 , Q = x 3 , I = x 4 and R = x 5 . Then, we rewrite the system (1) as We know that for R 0 = 1 the eigenvalues of the Jacobian matrix J corresponding to the system (1) at DFE E 0 (A/µ, 0, 0, 0, 0) are 0, −µ, −µ, λ ′ , λ ′′ , where λ ′ , λ ′′ are the roots of the following quadratic equation: When R 0 = 1 i.e. when a = a [B B] then β 2 γ 1 +γ 2 +σ 2 +µ < 1, which implies C 1 > 0. Again, by our assumption C 2 must be positive. Hence, J (E 0 ) has a simple zero eigenvalue and all other eigenvalues have negative real parts. Therefore, we can apply Castillo-Chavez and Song bifurcation theorem. Let, for the critical value a [B B] of a, W = (w 1 , w 2 , w 3 , w 4 , w 5 ) t and V = (v 1 , v 2 , v 3 , v 4 , v 5 ), respectively, be the right eigenvector and the left eigenvector corresponding to the zero eigenvalue of the Jacobian matrix J (E 0 ). Then The red and blue lines, respectively, represent the lines of unstable and stable equilibrium points. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) and the coefficient Since the coefficient ψ is obviously a positive number, hence the system (1) undergoes the backward bifurcation if φ > 0. Hence the theorem is proved. □ It is clear from Fig. 2 that there exists a critical value of the basic reproduction number, say R [C] 0 , such that for R [C] 0 < R 0 < 1, the system (1) contains two endemic equilibrium points where lower density is unstable and the upper density is stable. Thus, in this range two stable equilibrium points exist, among them one is disease free and the other is endemic. Hence, the condition R 0 < 1 is not sufficient to eradicate the disease, it depends on initial size of the populations. Again, it is clear from the expression of R 0 that it is a decreasing function of the cure rate of treatment (a), hence as the cure rate of treatment increases the value of R 0 will decrease and consequently the disease will be controlled. One can find three ranges of values of a such as 0 < a < a [C] , a [C] < a < a [B B] and a > a [B B] for which, the DFE is unstable and one endemic is stable (i.e. disease will not be eradicated from the population) for 0 < a < a [C] , two stable equilibrium points exist (i.e. elimination of disease depends on the initial size of the population) for a [C] < a < a [B B] and for a > a [B B] disease will be eradicated from the population.
For the considered values of the parameters if the cure rate of treatment (a) is greater than 32.01 then the disease will be eliminated from the population and for the range 2.501 < a < 32.01 eradication of disease depends on the initial population size and if a < 2.501 then disease will persist in the population. Thus eradication of the disease is highly dependent on the cure rate of treatment and hi-capacity of the cure rate of treatment is the indicative of disease elimination. Studying the bifurcation with respect to other parameters, it can be shown that the transmission rates of infections (β 1 , β 2 and β 3 ) and the auto recovery σ 1,2,3 (i.e. development of self immunity) play crucial role for controlling the disease. Due to the harmfulness of the extremely infectious disease COVID-19 the rate of infection depends on the interactions among the populations. So, lock down in the COVID-19 infected countries must play very important role to control the disease. As lock down is implemented in the population then interactions among the populations will be significantly reduced, which will be reflected in the expression of basic reproduction number (R L 0 ) of the optimal control problem (3).

Model validation, parameter estimation and prediction
In this section, we shall validate the COVID-19 model (1) using real data, estimate the model parameters and finally varying the sensitive model parameters we shall give some predictions. For this purpose here we have considered the coronavirus cases in Italy and we shall estimate the best fitted model parameters.
To estimate the model parameters, we use non-linear least square method. The principle of this method is to minimize the objective function where Σ is the set of all parameters, C [E] J (t) is the cumulative number of real infective cases, C [P] J (t) is the cumulative number of model predicted infective cases at J th observation, and n is the number of sample depending on the model parameters [4]. The cumulative number of infected can be found from the following formula: Using Matlab minimization and the Matlab package fmincon we have estimated the model parameters. To estimate the model parameters we consider the cumulative number of infected cases of Italy from 14th February to 23rd May, 2020. To utilize fmincon software we consider the initial populations size as S(0) = 60461826, E(0) = 5, Q(0) = 5, I (0) = 3, R(0) = 0 among them the E(0) is estimated and Q(0) is assumed and other three have been taken from [37]. The estimated model parameters and their sensitivity indices are given in Table 2. The model is fitted to the cumulative number of infected cases, which has been presented in Fig. 3(a) and the predictive number of cumulative and daily infected cases have been presented in Fig. 3 (c-e). The residuals of the fit is presented in Fig. 3(b), which shows that the residuals are small and random. Thus, we can conclude that the fit is reasonably good [24]. In Fig. 3, we have plotted the time series of infected population for 155 days after that a very few number of daily infected cases will come. It is clear from Fig. 3 that the state of being epidemic in Italy for COVID-19 will continue up to 3rd week of July, 2020 and after that the disease will be controlled. Since for these values of the model parameters R 0 > 1, hence disease will persist in population for long time. Now the goal is to bring it below 1, because an epidemic with a basic reproductive number below 1 will gradually disappear.
The considered model contains large number of parameters and almost all of them are involved in the basic reproduction number R 0 . Now, the problem is to decide which parameters are most influential to control the epidemic. But for experimental design of the epidemic models the determination of the most effective parameters is essential [29]. For this purpose we shall use the normalized forward sensitivity index [9] for R 0 with respect to any characteristic model parameter α, which is denoted by Γ α and is defined by ∂α . The sensitivity indices with respect to the model parameters are given in Table 2. It is clear from the table that the parameters β 1 , β 2 , β 3 , γ 1 enhance the value of R 0 and so they enhance the epidemic whereas other parameters with negative sensitivity indices will reduce the epidemic. The most sensitive parameters are β 2 (the transmission rate of infection from E class) and σ 2 (auto recovery rate of E class).  Table 2, we vary only the transmission rate of infections. Numerically,  we have verified that the transmission rate of infection from E(t) (asymptomatic and exposed) class is most effective but other two transmission rates are not so (see Fig. 4 (a-c)). It is realistic feature because I (t) and Q(t) classes cannot spread much more infections as E(t) can spread, because these two classes can infect only the medical persons and the persons who take care of them. So to minimize the disease spreading our target will be to minimize the contracts among the susceptible persons and the exposed & asymptomatic persons, which is possible if we increase the number of tests randomly in the areas where a single infected person is found, then the carrier of COVID-19 will be identified easily and they will be separated from the susceptible populations.
To control the spreading of disease Italy government adopted the lock down policy on 9th March, 2020 and daily number of cases starts to decrease on 23rd March i.e. after completion of the 14 days incubation period. Since the parameter σ 2 (auto recovery rate of E class) has high negative sensitivity index, hence it can change the disease dynamics rapidly. To justify this we change its value on 23rd March, then the prevalence significantly increases (decreases) as auto recovery rate of E class σ 2 decreases (increases) (see Fig. 5(a)). But the parameter γ 2 is not so effective as its low sensitivity index (see Fig. 5(b)). Thus, the developments of self immunity with the proper maintenance of social distancing of the exposed and asymptomatic individuals play key role for controlling the epidemic disease COVID-19.

Basic reproduction number of the COVID-19 outbreak in Italy
In this section, we shall estimate the basic reproduction number R 0 from the initial growth phase of actual epidemics in Italy. The effective reproduction number R(t) which is time dependent, will also be estimated from daily new COVID-19 cases in Italy.

Estimation of R 0 for actual epidemics
Here by using the initial COVID-19 infected cases in Italy, we shall estimate R 0 by applying the method developed by [27,31]. We suppose that the cumulative number of cases C(t) ∝ ex p(Λt), where Λ is the force of infection. Then the infected compartments can be taken in the following form: where E 0 , Q 0 , I 0 are the constants and we assume that at the earlier phase of the epidemic the abundance of the infected class is very low compared to the total susceptible populations. So, we assume that the number of susceptible populations equals the number of total populations at the time of initial infection and so S(t) = N (t) = A µ . We also assume that at beginning of the epidemics the delay parameter b in treatment is zero.
Putting the above values of S(t), E(t), Q(t) and I (t) in (1) and rearranging we get Determining β 2 from the above three expressions and putting in the expression of R 0 , we obtain . Following [27], we have the relation which is the number of new cases per day ∼ ΛC(t). Now we plot new daily COVID-19 cases against the cumulative number of COVID-19 cases C(t) in Italy from 14th February to 23rd May, 2020 (see Fig. 6). We fit a linear regression model by using the least square method to this exponential growth data. The force of infection Λ is nothing but the slope of the fitted line (see Fig. 6(b)). Thus, we estimate the force of infection as Λ = 0.1232311997657 ± 0.0067125 and from the above relation between the basic reproduction number and the force of infection, we can estimate the basic reproduction number as R 0 = 2.766467198 ± 1.0985.

The effective reproduction number
The basic reproduction number is the average number of secondary infection that is generated by a single infected person. But in reality the reproduction number varies with time for any epidemic system [31]. Depending on the reproduction number the planners update the epidemic control policy time to time. For this purpose here we have estimated the effective reproduction number (R(t)) using current COVID-19 infection data in Italy from 14th February to 23rd May, 2020. The estimation of R(t) can be done from the following renewal equation [27,33] where B(t) is the number of new cases in the tth day and g(τ ) is the generation interval distribution for the COVID-19 disease in Italy. The generation interval distribution g(t) is the probability distribution function of time from infection of a person to the secondary infection case by that person. For simplicity here we assume that b = 0. The generation interval distribution is the combination of the three exponential functions D 1 e −D 1 t , D 2 e −D 2 t and D 3 e −D 3 t in the following form : and τ > 0. The above relation is valid when the force of Using the estimated model parameters, we have calculated the effective reproduction numbers and presented them in Fig. 7. It is clear from the figure that the effective reproduction number diminishes below one after approximately 85 days from 14th February, 2020 and the epidemic with R 0 < 1 will gradually disappear.

The optimal control problem
To minimize the COVID-19 infections every country is adopting the lock down policy. The lock down is a policy that is to minimize the spreading of COVID-19 infections by controlling the movements of populations to maintain social distancing. For applying lock down the interactions among S(t) class with E(t), I (t), Q(t) classes are reduced and consequently the transmission rates of infections must decrease. To include this lock down effect in the model, we have replaced β 1 , β 2 and β 3 by β 1 (1 − u 1 ), β 2 (1 − u 2 ) and β 3 (1 − u 3 ), respectively where 0 ≤ u i ≤ 1(i = 1, 2, 3). Again, for adapting lock down each country is facing huge financial losses. For adapting complete lock down the loss is very huge, so some countries are adapting partial lock down. The u , i s are proportional to how many areas are implemented under lock down or how well lock down is maintaining. If there is no lock down then u , i s are all zero and if total lock down is implemented then u , i s are all 1. The main goal of this study is to minimize the COVID-19 infections as well as the cost or financial loss of implementing the controls. Thus, we reformulate the system (1) as an optimal control problem which is given below: Here, we construct the objective functional J , which has to be minimized, where The constants A 1 , A 2 , A 3 are, respectively, the loss due to the presence of E, Q, and I class and the constants A 4 , A 5 , A 6 , respectively, represent the loss due to implementing these controls. We assume that all the controls will be implemented for the time interval [0, T ] and they will be stopped after the time T . Now, the problem is to find optimal functions u * Proof. The integrand of the objective function J (u 1 , u 2 , u 3 ) is clearly a convex function. Also, the system (3) satisfies the Lipschitz property with respect to the state variables, because the solutions of the system (3) are bounded. It implies the existence of the optimal control functions. Hence the theorem is proved. □ Theorem 6. The optimal control functions u * i (i = 1, 2, 3) which minimize the objective functional J over the region U are u * 1 = max , 1 }} , where, S * , E * , Q * , I * , R * are, respectively, the optimum values of S, E, Q, I, R and (λ * 1 , λ * 2 , λ * 3 , λ * 4 , λ * 5 ) is the solution of the system (4) with the conditions (5) and N * = S * + E * + Q * + I * + R * . Proof. The Lagrangian of the optimal control problem is . Now, we define the Hamiltonian H for the problem which is as follows: In order to determine the adjoint equations with transversality conditions, we apply Pontryagin's maximum principle [16,21,28] with the transversality conditions λ i (T ) = 0, i = 1, 2, 3, 4, 5. Thus, the adjoint variables λ i , i = 1, 2, 3, 4, 5 can be obtained from the following: To determine optimal control functions, we use the optimality conditions ∂ H ∂u i = 0 (i = 1, 2, 3) and we get Here, S * , E * , Q * , I * , R * are, respectively, the optimum values of S, E, Q, I, R and (λ * 1 , λ * 2 , λ * 3 , λ * 4 , λ * 5 ) is the solution of the system (4) satisfying the conditions (5).
Thus, the control functions u * i (i = 1, 2, 3) minimize the objective functional J . Hence the theorem is proved. □ To justify the theoretical findings of the optimal control problem (3), we solve it numerically by using forwardbackward sweep method. This method combines the forward application of a fourth order Runge-Kutta method for the state system (3) with the backward application of a fourth order Runge-Kutta method for the adjoint system (4). Here, we assume T = 1 unit of time. For this simulation we have considered the parametric values which are given in Table 2 10. Effect of the controls u 2 and u 3 on R 0 , when u 1 = 0.01. exposed and asymptomatic, quarantined and infectious populations at any time t ∈ [0, 1] between the system with no control and the system with the controls. The optimal control functions u * i (t) (i = 1, 2, 3) have been represented by Fig. 9(a)-(c).
Again, it can be easily proved that the basic reproduction number R L 0 of the model (3) is It is clear from the expression of the basic reproduction number R L 0 that the quantities u , i s can reduce the value of the basic reproduction number i.e. these quantities must be good controller for the disease spreading. If values of u , i s increase then R L 0 decreases i.e. the lock down must play an important role for controlling the disease. In Fig. 10, we see that the control u 2 is so effective to control the disease, because if u 2 increases, the basic reproduction number significantly decreases and if we make it less than unity for large value of u 2 , then disease will be controlled.

Conclusions
This paper deals with an epidemic compartmental model to investigate the dynamics of pandemic COVID-19 and to predict the future epidemic situations in Italy. Here, we have considered three disease compartments, namely exposed-asymptomatic, quarantine and infectious. To include the effect of limited medical facilities, here we have considered the rate of treatment as saturated type. The basic reproduction number has been obtained for the proposed model, which plays key role for analysing the stability of DFE. The system experiences backward bifurcation with respect to the cure rate of treatment conditionally. The epidemiological implication of this result is that making the basic reproduction number less than unity is not sufficient to eliminate the COVID-19 disease from the community. Thus the basic reproduction number less than unity is only a necessary but not a sufficient condition for eliminating COVID-19, because the stable disease free equilibrium point co-exists with a stable endemic equilibrium. To validate the model and for further predictions we consider the current COVID-19 cases in Italy from 14th February to 23rd May, 2020 and we have estimated the best-fitted model parameters. From the simulations we see that the disease will be controlled in July, 2020 in Italy. The sensitivity analysis shows that the transmission rate of infection from exposed and asymptomatic class plays important role for disease spreading and the development of self immunity of that class is very influential in the model output. We have estimated the basic reproduction number R 0 from the initial growth phase of actual COVID-19 epidemics in Italy and the time dependent effective reproduction number R(t) has also been estimated from daily new COVID-19 cases in Italy.
We have also formulated an optimal control problem by considering the effects in the transmission rates of infections due to implementing lock down policy by many countries. We have solved it both analytically and numerically. The main objective of analysing this optimal control problem is to minimize the COVID-19 positive cases as well as the total costs or financial losses for implementation of controls due to adapting lock down by countries. A comparative study between the system with lock down and the system without lock down has been presented that shows the positive impacts of taking lock down for controlling the epidemic. The effects of lock down on the basic reproduction number have also been shown and it has been shown that the basic reproduction number decreases as the lock down policy is implemented and if we can make it less than unity then the epidemic will disappear. The formulation of optimal control problem is a theoretical modelling and it can be applied by using real data in some countries. The proposed model of COVID-19 has some limitations. Here we consider saturated treatment function to describe the saturation phenomenon of the limited medical resources. If some countries or communities have the potential to adapt very strong medical resources in very short span of time then the saturated treatment function will not be used in the model. Again, we did not consider vaccination term in our model though the COVID-19 vaccines have reached billions of people worldwide.
(ii) If we introduce an extra compartment D for the deaths, then the model equations will be as follows: = d 1 I Thus, if we introduce the death compartment D, then the dynamics of the model will not be changed, because the last equation is a redundant equation as other equations do not contain D.
But at the time of parameter estimation, we will be able to predict the probable number of death cases in current time as well as the future time. This prediction is important biologically, because we can validate the model by comparing the result of model simulated death cases with the exact number of death cases.