key: cord-1055840-2nllnex6 authors: Wintachai, Phitchayapak; Prathom, Kiattisak title: Stability analysis of SEIR model related to efficiency of vaccines for COVID-19 situation date: 2021-04-15 journal: Heliyon DOI: 10.1016/j.heliyon.2021.e06812 sha: f32b620f6d4f1cf07244d325e074d1180f08d21d doc_id: 1055840 cord_uid: 2nllnex6 This work is aimed to formulate and analyze a mathematical modeling, [Formula: see text] model, for COVID-19 with the main parameters of vaccination rate, effectiveness of prophylactic and therapeutic vaccines. Global and local stability of the model are investigated and also numerical simulation. Local stability of equilibrium points are classified. A Lyapunov function is constructed to analyze global stability of the disease-free equilibrium. The simulation part is based on two situations, the US and India. In the US circumstance, the result shows that with the rate of vaccination 0.1% per day of the US population and at least 20% effectiveness of both prophylactic and therapeutic vaccines, the reproductive numbers [Formula: see text] are reduced from 2.99 (no vaccine) to less than 1. The same result happens in India case where the maximum reproductive number [Formula: see text] in this case is 3.38. To achieve the same infected level of both countries, the simulation shows that with the same vaccine's efficiency the US needs a higher vaccination rate per day. Without vaccines for this pandemic, the model shows that a few percentages of the populations will suffering from the disease in the long term. Coronavirus disease is a severe acute respiratory disease caused by a coronavirus 2 (SARS-CoV-2) that is a new member of the genus Beta coronavirus and family Coronaviridae [1, 2] . The virus primarily spreads from person to person through droplet, airborne, and contact transmission. The clinical symptoms of SARS-CoV-2 infected patients had mild, moderate, and severe symptoms such as fever, dry cough, difficulty breathing, fatigue, new loss of taste or smell, nausea, diarrhea, pneumonia, and respiratory symptom [3, 4] . The critically severe conditions such as chronic medical illness, organ dysfunctions, and death have been frequently reported in elderly patients and people with immunodeficiencies [5, 6] . However, many SARS-CoV-2 infected patients are minimally symptomatic or asymptomatic [7, 8, 9] . The outbreak of SARS-CoV-2 started in China and then transmitted to humans and animals [10, 11] . Nowadays, the virus has recently caused epidemics around the world in more than 215 countries with 46,403,652 confirmed cases and 1,198,569 mortalities, as of November 2, 2020 [12] . According to a recent report from the World Health Organization on SARS-CoV-2 outbreak, the number of confirmed cases in America is higher than the number of confirmed cases in Europe, * Corresponding author. E-mail address: kp525914@ohio.edu (K. Prathom). South-East-Asia, Eastern Mediterranean, Africa, and Western Pacific, respectively. The top two countries reporting the most confirmed cases are the United States of America (9,032,465 cases) and India (8,229,313 cases). As COVID-19 cases keep increasing, predictions of the number of infected cases and the termination of COVID-19 are worth it to study. Mathematical model of infectious diseases is a crucial tool that has been used to study dynamics of how diseases spread. A mathematical model can predict the future situation of an outbreak and evaluate the best strategy to control spreading diseases. There are many different types of mathematical models for predicting an epidemic infection. One of them is called compartment models. Compartment model is an interesting tool for COVID-19 situation. It is a powerful mathematical model for understanding the complex dynamics of epidemics. In this work we construct a well known model called model with 4 compartments of susceptible population , exposed population , infectious population , and recovered population . The model is suitable for disease transmission which an infected individual needs a short time period to be an infectious. Many researches have been studied by adapting model to forecast dynamics of endemic and epidemic such as Dengue Fever [13, 14, 15] , Ebola [16, 17] , Middle East Respiratory Syndrome (MERS) [18, 19] Severe Acute Respiratory Syndrome (SARS) [20, 21] , to name a few. According to [22, 23, 24] , COVID-19 has an average incubation period of 11.5 days before spreading of the viruses, so model is suitable for predicting COVID-19 situation. model have been adapted by adding strategy parameters such as social distancing and face mask using to control and predict COVID-19 situation in several researches [25, 26, 27, 28, 29, 30] . Vaccine administration is a highly effective method of preventing and reducing viral infections [31] . Even though there is no vaccine or a specific antiviral for the treatment of patients infected with SARS-CoV-2 available, several vaccines against SARS-CoV-2 such as a live attenuated vaccine, inactivated vaccines, subunit vaccines, DNA and RNA vaccines, and vector vaccines have been developed [32, 33] . Vaccination and optimal control are key points to control an epidemic situation as discussed in [34, 35, 36, 37] . In this study, we use SEIR model equipped with effectiveness of vaccination to forecast COVID-19 situation when a vaccine comes out. There are two main types of vaccine in our model prophylactic and therapeutic vaccines. Prophylactic vaccine is a preventing vaccine and therapeutic vaccine is a vaccine that is administrated after infection [38] . In Section 2, we formulate a model and investigate all equilibrium points together with their global and local stability of the model. Section 3 is mathematical simulation part. In particular, we applied recorded parameters of US and India circumstances to our model and predicted the potential of COVID-19 in both countries when vaccines come out. We consider the 4-compartment model called which ( ), ( ), ( ), and ( ) are the fractions of susceptible, exposed, infectious, and recovered populations, respectively, at the time . The trivial solution ≡ 0, ≡ 0, ≡ 0 and ≡ 0 is out of our interest. The system of differential equations related to the schematic diagram in Fig. 1 is as follows: under the conditions that The density ( ) at the time is the faction of susceptible numbers, ( ), ( ) and ( ) are similar. The explanation of variables and parameters in (1) are presented in Table 1 . A motivation for the model is that vaccination rate per day ( ) cannot terminate the flow of the system immediately since the whole population cannot be vaccinated at once. A person can get vaccinated once he or she is susceptible, exposed, or infectious. In the first equation of System (1), the rate of change in Death rate of exposed population plus 0 2 Death rate of infectious population plus 0 susceptible depends on the numbers of vaccinated humans, , and non-vaccinated humans, (1 − ) . Note that by the fundamental existence-uniqueness theorem for nonlinear systems, the nonlinear system (1) has a unique solution set ( ( ), ( ), ( ), ( )). To guarantee that the densities ( ), ( ), ( ), and ( ) in the model (1) are non-negative at any time > 0, we prove the following lemma. Proof. To prove this lemma we use the fact that a function with (0) ≥ 0 is a non-negative function if | = * ≥ 0 when ( * ) = 0; i.e., the function is non-decreasing at * . By the condition (2), there is such that ( ) ≥ 0 on 0 ≤ < and ( ) = 0. Based on the first equation of Model (1), we have It implies that ( ) ≥ 0 for any ≥ 0. Next, let be the time such that ( ) ≥ 0 on 0 ≤ < and ( ) = 0. By the third equation of (1), we have Since and are non-negative on [0, ], it follows by the second equation of (1) that ]. This implies that Equations (3) and (4) imply that | = ≥ 0, so ( ) ≥ 0 for any ≥ 0. It is easy to check that ( ) ≥ 0 when ( ) ≥ 0. Since , , and are non-negative for > 0, it is obvious that ( ) ≥ 0 for ≥ 0. □ From this lemma we can conclude that the set [0, ∞) 4 is positive invariant with respect to the model (1) and it attracts all solutions of the model. Equilibrium points of the system can be found by setting = 0, = 0, = 0, and = 0 of (1); i.e., solving the following system: A disease-free equilibrium is an equilibrium when there is no spread of the disease; i.e., ≡ 0 ≡ . By solving (5), the disease-free equilibrium is unique in the form for fixed parameters 0 , , and 0 . Apart from the disease-free equilibrium, others equilibrium points (endemic equilibrium) of the model can be found by solving (5) be a representation of the equilibrium point in the form of (6) or (7). The following theorem describes stability of the equilibrium point (8). Proof. Consider the Jacobian matrix of the model (1) with respect to the equilibrium point (8) which is as follows: where 1 = 3 * * + ( + 1 + 2 + 2 * ) 2 = + 1 + 2 + 1 2 + 2 * + ( + + 1 + 2 ) 3 = + + 1 + 2 + . It is obvious that 1 > 0, 2 > 0, and 3 > 0. Since 1 , 2 , 3 are positive real numbers, it follows that all solutions of Equation (10) have negative real parts. Therefore, the equilibrium point of the model (1) is locally asymptotic stable. □ Using the matrices generation method [39] , the basic reproductive number,  0 , is the dominant eigenvalue (the spectral radius) of and Hence, the basic reproductive number,  0 , corresponding to the disease-free equilibrium (6) is in the form By (11) and (12), we note here that the dominant eigenvalues of ] where 2 , , , and are parameters defined in Table 1 . Note that is a 1 × 2 matrix of positive real components. It is easy to check that where and are defined in (11) and (12), respectively. Equation (14) implies that Next, let We note here that  is a zero matrix only at the disease-free equilibrium. By using (1), we have Define the Lyapunov function  as follows: [42] 0.32 [43] 1/11.5 per day [22] 1/11.5 per day [22] 0.0686 per day [12] 0.0686 per day [ Since −1 is a 1 × 2 matrix of positive real components and  is a non-negative matrix, it follows that  ≥ 0 and we also have that  = 0 if and only if = 0 and = 0. This implies that  is positive definite. Moreover, by (17) and (15) we obtain Since  < 0 if  0 < 1, it follows that the disease-free equilibrium (6) is globally asymptotic stable [40] . On the other hand, if  0 > 1, then  > 0 which implies that the equilibrium is unstable. Note that in the case of  0 = 1, we can conclude that the equilibrium is locally stable since  = 0. □ The maximum reproductive number  0 with respect to the diseasefree equilibrium (6) occurs when there is zero vaccination ( = 0); that is, Based on Table 2 , the maximum  0 of Case I is 2.99 and the maximum  0 of Case II is 3.38. By increasing the vaccination rate ( ), the values of  0 are decreasing corresponding the effectiveness of prophylactic ( ) and therapeutic ( , ) vaccines, see Fig. 2 . We note here as an example that = 0.4 means 40% effectiveness of prophylactic when applied to susceptible ( ); i.e., if 100 people in are administered a prophylactic vaccine, it will be 40 people recovered. Fig. 2 shows that if the vaccination rate ( ) is under 0.0002 or 0.02% per day of the populations, the basic reproductive numbers are still higher than 1 no matter how much vaccine's effectiveness is. If the rate of vaccination per day is increased (0-5% in the US and 0-1% in India), see Fig. 3 , the infection rate in both countries will be decreased. This implies that the vaccination rate are play an important role to terminate the pandemic. However, the vaccine efficacy is also important, the effectiveness can vary the risk of illness [48] . SARS-CoV-2 is a zoonotic infection that has transmitted from a vertebrate to a human [49] . During outbreak, the virus infections in humans have been reported at higher rates than animals infections. SARS-CoV-2 may use animals and humans as reservoirs for reemerging, similar to SARS coronavirus [50, 51] . Thus, coronavirus disease may be a re-emerging viral diseases which is a diseases that has been observed previously within a population. To completely control SARS-CoV-2 infection, strategies for increasing vaccination rates is interesting to investigate for effective infection prevention and control of the disease. We need a sufficient vaccination rate depending on the power of vaccines and several doses of vaccine might be recommended. Based on Fig. 2 , if we had a vaccine effectiveness higher than 20%, it would suffice to proceed 0.1% of the populations per day to reduce the basic reproductive numbers to be under 1. Moreover, with the same vaccine's efficiency of 70% of prophylactic and 60% of therapeutic, the US need higher rate of vaccination than India to flatten the curve as seen in Fig. 3. Fig. 3 . Fractions of US infectious cases and India infectious cases over time if we had 70% and 60% effectiveness of prophylactic and therapeutic vaccines, respectively, with different rates of vaccination, 0% (no vaccination process), 0.1%, 1%, and 5% (only US) per day of each population. The equilibrium point related to the US and India situations can be computed by using Equation (7) . With the vaccination rate 0.1% per day of the US population ( = 0.001) and 90% efficiency of prophylactic and therapeutic vaccines, the equilibrium point corresponding to the fixed parameters in Table 2 of the US case is ( * , * , * , * ) = (0.0344, 0, 0, 1.3026). If there is no vaccine, the equilibrium point of the US case is (0.1486, 0.0003, 0.0004, 1.1876), that is the disease will not die out eventually. In the long term, there are about 0.04% infectious of the US population. India's case has ( * , * , * , * ) = (0.0532, 0, 0, 2.4032) for = 0.001 and 90% vaccines' efficiency and it has ( * , * , * , * ) = (0.2145, 0.0005, 0.0006, 2.2407) for no vaccines. Similarly to the US, a few percentages (0.06%) of India's population are infectious in the long term if there is no vaccine. The difference between efficiency of prophylactic and therapeutic vaccines in human SARS-CoV-2 infection treatment is depicted as in Fig. 4 . The effectiveness of both vaccines was set to the same values. The results showed that prophylactic vaccine has higher efficiency than therapeutic vaccine in both the US and India. Prophylactic vaccine will stimulate the immune system and then produce long-lived memory lymphocytes [52, 53] . Subsequently, the immune system can rapidly respond to virus infection, leading to a reduction of infected cases. The constructed model in this manuscript gives a future perspective when we have a vaccine for COVID-19. The simulation shows that having an effective vaccine significantly flatten the peak of infectious population. The model illustrates that having a vaccine does not immediately terminate the pandemic. It takes a period of time depending on the effectiveness of the arrival vaccine and the vaccination rate. Under the same vaccine's effectiveness, the simulation shows that the US need vaccination rate per day higher than the rate used in India to achieve the same result. According to the formula of  0 in Section 2, we see that the vaccination rate and the efficiency of vaccines play an important role to reduce the value of  0 . The theoretical results have confirmed that when the reproductive number  0 of the pandemic is less than 1, the COVID-19 situation will be under control; i.e., the model is stable. The equilibrium point of the model for specific parameters gives the stationary flow of the pandemic situation in the long term that a few percentages of the considered populations will be infectious if we have no vaccine for COVID-19. Author contribution statement P. Wintachai: Contributed reagents, materials, analysis tools or data; Wrote the paper. K. Prathom: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Wrote the paper. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. This work was supported by Walailak University grant no. WU64219. The authors declare no conflict of interest. No additional information is available for this paper. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in A novel coronavirus from patients with pneumonia in China Pathological findings of Covid-19 associated with acute respiratory distress syndrome Clinical features of patients infected with 2019 novel coronavirus in Clinical characteristics of elderly patients with Covid-19 in Hunan province, China: a multicenter, retrospective study If you are immunocompromised, protect yourself from Covid-19 Clinical characteristics of 24 asymptomatic infections with Covid-19 screened among close contacts in Nanjing A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster Estimation of the asymptomatic ratio of novel coronavirus infections (Covid-19) SARS-CoV-2 infection in farmed minks, the Netherlands Sars-CoV-2 in fruit bats, ferrets, pigs, and chickens: an experimental transmission study Coronavirus disease (Covid-19) dashboard SEIR model for transmission of Dengue fever in Selangor Malaysia Analysis of a Dengue disease transmission model Mathematical model of Dengue virus with primary and secondary infection A modified SEIR model for the spread of Ebola in western Africa and metrics for resource allocation Modeling the spread of Ebola with SEIR and optimal control Applying discrete SEIR model to characterizing MERS spread in Korea Mathematical model for MERS-COV disease transmission with medical mask usage and vaccination A double epidemic model for the Sars propagation Mathematical modelling of Sars and other infectious diseases in China: a review The incubation period of coronavirus disease 2019 (Covid-19) from publicly reported confirmed cases: estimation and application Incubation period of Covid-19: a rapid systematic review and meta-analysis of observational research Incubation period of the coronavirus disease 2019 (Covid-19 SEIR model for Covid-19 dynamics incorporating the environment and social distancing A simulation of a Covid-19 epidemic based on a deterministic SEIR model SEIR modeling of the Covid-19 and its dynamics Analysis of a mathematical model for Covid-19 population dynamics in Stability analysis and numerical simulation of SEIR model for pandemic Covid-19 spread in Indonesia Can the Covid-19 epidemic be controlled on the basis of daily test reports? Understanding how vaccines work Therapeutic and vaccine strategies against Sars-CoV-2: past, present and future Covid-19: progress in diagnostics, therapy and vaccination Determination of an optimal control strategy for vaccine administration in Covid-19 pandemic treatment On vaccination controls for the SEIR epidemic model Vaccination strategies based on feedback control techniques for a general SEIR-epidemic model Possible treatment of Covid-19 with a therapeutic vaccine Reproduction numbers of infectious disease models Estimation of infection rate and predictions of disease spreading based on initial individuals infected with Covid-19 Covid-19 pandemic scenario in India compared to China and rest of the world: a data driven and model analysis, medRxiv Vaccine effectiveness: How well do the flu vaccines work Coronavirus disease 2019 (Covid-19): situation report A review of studies on animal reservoirs of the Sars coronavirus Covid-19 pandemic: current knowledge about the role of pets and other animals in disease transmission Induction and maintenance of long-term immunological memory following infection or vaccination The authors would like to thank Walailak University for financial support under grant WU64219. We also thank the reviewers for invaluable comments and suggestions that greatly improved the manuscript.