Modelling the Spreading of the SARS-CoV-2 in Presence of the Lockdown and Quarantine Measures by a “Kinetic-Type Reactions” Approach ======================================================================================================================================== * Sonnino Giorgio * Peeters Philippe * Nardone Pasquale ## Abstract We propose a realistic model for the evolution of the COVID-19 pandemic subject to the lockdown and quarantine measures, which takes into account the time-delay for recovery or death processes. The dynamic equations for the entire process are derived by adopting a *kinetic-type reactions* approach. More specifically, the lockdown and the quarantine measures are modelled by some kind of inhibitor reactions where susceptible and infected individuals can be *trapped* into inactive states. The dynamics for the recovered people is obtained by accounting people who are only traced back to hospitalised infected people. To get the evolution equation we take inspiration from the Michaelis-Menten’s enzyme-substrate reaction model (the so-called *MM reaction*) where the *enzyme* is associated to the *available hospital beds*, the *substrate* to the *infected people*, and the *product* to the *recovered people*, respectively. In other words, everything happens as if the hospitals beds act as a *catalyzer* in the hospital recovery process. Of course, in our case the reverse *MM reactions* has no sense in our case and, consequently, the kinetic constant is equal to zero. Finally, the O.D.E.s for people tested positive to COVID-19 is simply modelled by the following kinetic scheme *S* + *I* ⇒ 2*I* with *I* ⇒ *R* or *I* ⇒ *D*, with *S, I, R*, and *D* denoting the compartments Susceptible, Infected, Recovered, and Deceased people, respectively. The resulting *kinetic-type equations* provide the O.D.E.s, for elementary *reaction steps*, describing the number of the infected people, the total number of the recovered people previously hospitalised, subject to the lockdown and the quarantine measure, and the total number of deaths. The model foresees also the second wave of Infection by Coronavirus. The tests carried out on real data for Belgium, France and Germany confirmed the correctness of our model. Key words * Mathematical model * COVID-19 * Dynamics of population * Pneumonia ## 1 Introduction Coronavirus disease 2019 (COVID-19) is caused by a new Coronavirus (SARS-CoV-2) that has spread rapidly around the world. Most infected people have no symptoms or suffer from mild, flu-like symptoms, but some become seriously ill and can die. In recent weeks coronavirus has had too many opportunities to spread again. After successfully tamping down the first surge of infection and death, Europe is now in the middle of a second coronavirus wave as it moves into winter [1], [2], [3], [4], [5], [6]. Even though several vaccines for COVID-19 are actually been produced other ways of slowing its spread have to continue to be explored. One way of controlling the disease are the lockdown and the quarantine measures. The lockdown measures are emergency measures or conditions imposed by governmental authorities, as during the outbreak of an epidemic disease, that intervene in situations where the risk of transmitting the virus is greatest. Under these measures, people are required to stay in their homes and to limit travel movements and opportunities for individuals to come into contact with each other such as dining out or attending large gatherings. The lockdown measures are more effective when combined with other measures such as the quarantine. Quarantine means separating healthy people from other healthy people, who may have the virus after being in close contact with an infected person, or because they have returned from an area with high infection rates. Similar recommendations include isolation (like quarantine, but for people who tested positive for COVID-19) and physical distancing (people without symptoms keep a distance from each other). Several governments have then decided that stricter lockdown and quarantine measures are needed to bring down the number of infections. In this work we shall propose interventions which are as targeted as possible. Unfortunately, the greater the number of infections, the more sweeping the measures have to be. Tightening the measures will impact on our society and the economy but this step is needed for getting the coronavirus under control. The aim of this work is to model the dynamics of the infectious, recovered, and deceased people when population is subject to lockdown and quarantine measures imposed by governments. We shall see that the combined effect of the restrictions measures with the action of the Hospitals and Health Institutes is able to contain and even dampen the spread of the SARS-CoV-2 epidemic. The dynamics of the entire process will be obtained by taking into account the theoretical results recently appeared in literature [7] and [8] and by adopting a *kinetic-type reactions* approach. In this framework, the dynamics of the Health Institutes is obtained by taking inspiration from the Michaelis-Menten’s enzyme-substrate reaction model (the so-called *MM reaction* [9], [10], and [11]) where the *enzyme* is associated to the *available hospital beds*, the *substrate* to the *infected people*, and the *product* to the *recovered people*, respectively. In other words, everything happens as if the hospitals beds act as a *catalyzer* in the hospital recovery process [12]. In addition, the time-delay for recovery or death processes are duly taken into account. More specifically, in our model, we have the following 10 compartments: *S* = Number of susceptible people. This number concerns individuals not yet infected with the disease at time *t*, but they are susceptible to the disease of the population; *I* = Number of people who have been infected and are capable of spreading the disease to those in the susceptible category; *I**h* = Number of hospitalised infected people; *I**Q* = Number of people in quarantine. This number concerns individuals who may have the virus after being in close contact with an infected person; *R* = Total number of recovered people, meaning specifically individuals having survived the disease and now immune. Those in this category are not able to be infected again or to transmit the infection to others; *r**h* = Total recovered people previously hospitalised; *D* = Total number of people dead people for COVID-19; *d**h* = Total number of people previously hospitalised dead for COVID-19; *L* = Number of inhibitor sites mimicking lockdown measures: *Q* = Number of inhibitor sites mimicking quarantine measures. In addition, *N*, defined in Eq. (19), denotes the number of total cases. The manuscript is organised as follows. In Section 2 we derive the deterministic Ordinary Differential Equations (ODSs) governing the dynamics of the infectious, recovered, and deceased people. The lockdown and quarantine measures are modelled in Subsection 2.2. The dynamics of the hospitalised individuals (i.e., the infectious, recovered, and deceased people) can be found in Subsection 2.4. As mentioned above, the corresponding ODEs are obtained by considering the *MM reaction model*. The equations governing the dynamics of the full process and the related *basic reproduction number* are reported in Section 3 and Section 4, respectively. It is worth mentioning that our model foresees also the second wave of Infection by Coronavirus. As shown in Section 5, in absence of the restrictive measures and by neglecting the role of the Hospitals and the delay in the reactions steps, our model reduces to the classical *Susceptible-Infectious-Recovered-Deceased-Model* (SIRD-model) [13]. Finally, Section 6 shows the good agreement between the theoretical predictions with real data for Belgium, France and Germany. The last Section 7 presents the conclusions and perspectives of this manuscript. ## 2 Model for COVID-19 in Presence of the Lockdown and Quarantine Measures As mentioned in the Introduction, the population is assigned to compartments with labels *S, I, R D* etc. The dynamics of these compartments is generally governed by deterministic ODEs, even though stochastic differential equations should be used to describe more realistic situations [7]. In this Section, we shall derive the deterministic ordinary differential equations obeyed by compartments. This task will be carried out by taking into account the theoretical results recently appeared in literature [8], [14] and without neglecting the delay in the reactions steps. ### 2.1 Modelling the Susceptible People If a susceptible person encounters an infected person, the susceptible person will be infected as well. So, the scheme simply reads ![Formula][1] ### 2.2 Modelling the Lockdown and Quarantine Measures The lockdown measures are mainly based on the isolation of the susceptible people, (eventually with the removal of infected people by hospitalisation), but above all on the removal of susceptible people. ### Subsection 2.1. Modelling the Lockdown and Quarantine Measures with Chemical Interpretation It is assumed the lockdown and quarantine measures are modelled by some kind of inhibitor reaction where the susceptible people and the infected can be *trapped* into inactive states *S**L* and *I**Q*, respectively. Indicating with *L* and *Q* the Inhibitor sites mimicking the lockdown and the quarantine measures respectively, we get ![Formula][2] In the scheme (2), symbol =⇒ stands for a *delayed reaction* just like *enzyme degradation processes* for instance. Here, *L**max* = *S**L* + *L* hence, if *L ≃ L**Max*, an almost perfect lockdown measures would totally inhibit virus propagation by inhibiting all the susceptible people *S* and the infected people *I*. A not so perfect lockdown measures would leave a fraction of *I* free to spread the virus. The number of inhibitor sites maybe a fraction of the number of the infected people. Fig. 1. shows the behaviour of the lockdown efficiency parameter adopted in our model. For simplicity, we have chosen a parameter which is constant *k**LMax* ≠ 0 inside the time-interval *t*1 ≤ *t* ≤ *t*2 and vanishes outside it. The *inverse Lockdown efficiency parameter* is ![Graphic][3], which is equal to *k**LMax* outside the door and vanishes inside the the interval *t*1 ≤ *t* ≤ *t*2. Finally, from Schemes (1) and (2), we get the O.D.E.s for *S, L, Q*, and *I**Q*: ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F1) Figure 1: Lockdown Efficiency Parameter. For simplicity, in our model the lockdown efficiency parameter kL is a door-step function. This function is constant, KLMax ≠ 0,within the range t1 ≤ t ≤ t2 and zero outside it. ![Formula][4] ![Formula][5] with the *dot* above the variables denoting the *time derivative*. ### 2.3 O.D.E. for the Total Recovered People At the first approximation, the O.D.E. for the *total recovered people R* (i.e. the total individuals having survived the disease) is trivially obtained by considering the following *kinetic scheme*: ![Formula][6] That is, the rate of *R**t* is approximatively proportional to the number of the infected people *I* at time *t* i.e.1. ![Formula][7] where we have introduced the time-delay *t**R* (the number of the recovered people at time time *t* is proportional to the infected people at time *t* − *t**R*). However, it is useful to clarify the following. In Eqs (5), *R* stands for the *total number of the recovered people* (i.e. the number of the recovered people previously hospitalised, plus the number of the asymptomatic people, plus the infected people who have been recovered without being previously hospitalised). The natural question is: *how can we count R and compare this variable with the real data ?*. The current statistics, produced by the Ministries of Health of various Countries, concern the people released from the hospitals. Apart from Luxembourg (where the entire population has been subject to the COVID-19-test), no other Countries are in a condition to provide statistics regarding the total people recovered by COVID-19. Hence, it is our opinion that the equation for *R*, is not useful since it is practically impossible to compare *R* with the experimental data. We then proceed by adopting approximations and to establish the differential equation whose solution can realistically be subject to experimental verification. More specifically: Firstly, we assume that *R* is given by three contributions: ![Formula][8] with *r**h*, *r**A*, and *r**I* denoting the *total number of the recovered people previously hospitalised, the total number of asymptomatic people*, and the *total number of people immune to SARS-CoV-2*, respectively. Secondly, we assume that the two contributions *r**A* and *r**I* are negligible i.e. we set *r**A* ≈ 0 and *r**I* ≈ 0 2. ### 2.4 O.D.E. for the Recovered People in the Hospitals Now, let us determine the dynamics for the recovered people in the hospitals. So, we account people who are only traced back to hospitalised infected people. We propose the following model3: ![Formula][9] with *b**h* denoting the number of available *hospital beds, I* the number of *infected people, I**h* the number of *infected people blocking an hospital bed, r**h* the number of *recovered people previously hospitalised*, and *d**h* the number of *people deceased in the hospital*. Of course, ![Formula][10] The dynamic equations for the processes are then: ![Formula][11] where *t**r* and *t**d* are the *average recovery time delay* and the *average death time delay*, respectively, and we have taken into account Eq. (9) i.e., *b**h* = *C**h* − *I**h*. In general *t**r* ≠ *t**d* ≠ 0. Of course, the variation of *r*(*t*) over a period Δ*t* is: ![Formula][12] ### 2.5 O.D.E. for People Tested Positive to COVID-19 The number of the infected people may be modelled by the following *kinetic scheme* ![Formula][13] The scheme (12) stems from the following considerations 1. If a susceptible person encounters an infected person, the susceptible person will be infected; 2. The infected people can either survive and, therefore, be recovered after an average time-delay *t**R*, or die after an average time-delay *t**D*; 3. The schemes (2) and (8), respectively, have been taken into account. 4. The differential equation for the infected people is reads then ![Formula][14] ### 2.6 O.D.E. for Deaths In this model, we assume that the rate of death is proportional to the infected people, according to the scheme (12). By also taking into account the scheme (2), we get ![Formula][15] and the corresponding O.D.E. for deaths reads ![Formula][16] ## 3 Set of O.D.E.s for the Spread of SARS-CoV-2 when the Lockdown and the Quarantine Measures are Adopted By collecting the above O.D.E.s, we get the full system of differential equations governing the dynamics of the number of the infected people, the total number of the recovered people previously hospitalised and the total number of deceased peopled, when the lockdown and the quarantine measures are adopted ![Formula][17] From Eqs (16) we get ![Formula][18] or, by taking into account that *S* + *S**L* = *S**T ot*., *R* + *r**h* = *R**T ot*., *D* + *d**h* = *D**T ot*., and *I* + *I**Q* + *I**h* = *I**T ot*. we get ![Formula][19] The number of total cases *N* is defined as ![Formula][20] ## 4 The Basic Reproduction Number We note that, in absence of the lockdown and the quarantine measures, the dynamics of the infectious class depends on the following ratio: ![Formula][21] with *N**T ot*. denoting the *Total Population. R* is the *basic reproduction number*. This parameter provides the expected number of new infections from a single infection in a population by assuming that all subjects are susceptible [2], [3]. The epidemic only starts if *R* is greater than 1, otherwise the spread of the disease stops right from the start. ## 5 Comparison with the SIRD model The *Susceptible-Infectious-Recovered-Deceased-Model* (SIRD-model) is one of the simplest compartmental models, and many models may be derived from this basic form. According to the SIRD model, the dynamic equations governing the above compartments read [13] ![Formula][22] It is easily checked that Eqs (16) reduce to Eqs (21) by adopting some assumptions. In particular: 1. The system is not subject to the lockdown and quarantine measures; 2. The average times-delay may be neglected; 3. Hospitals do not enter in the dynamics. Under these assumptions, Eqs (16) reduce to the SIRD equations: ![Formula][23] ## 6 Application of the Model and Appearance of the Second Wave of SARS-CoV-2 Infection Let us now apply our model to the case of a small Country, Belgium, and to other two big Countries, France and Germany. Real data are provided by the various National Health agencies (Belgium - *Sciensano* [15]; France - *Santé Publique France* [16]; Germany -*Robert Koch Institut. Country data from Worldbank*.*org* [17]) and compiled, among others, by European Centre for Disease Prevention and Control (ECDC). It should be noted that this measures does not generally provide the true new cases rate but reflect the overall trend since most of the infected will not be tested [18]. It should also be specified that real data provided by ECDC refer to the *new cases per day*, which we denote by Δ*I**new*(*t*). By definition, Δ*I**new*(*t*) corresponds to the new infected people generated from step ![Graphic][24] solely during 1 day, and *not* to the compartment *I*. Hence, the ECDC data have to be confronted vs the theoretical predictions provided by the solutions for *S*(*t*) and *S**L*(*t*) of our model, according to the relation Δ*I**new*(*t*) = −Δ*S*(*t*) − Δ*S**L*(*t*). The values of the parameters used to perform these comparisons are shown in Table 1. Initial *µ* en *k*1 values have been estimated (fitted) from the measurements using the short period at the start of the pandemic using simple solution valid during that period. *I*(60) (from March 1, 2020). Hospital capacity is evaluated from the different Countries published capacity. However, we are aware that the interpretation may vary from one Country to another. During the first lockdown, Countries have taken various actions to limit Coronavirus spreading (social distancing, wearing masks, reducing high density hotspots etc.). In order to include these measures in a simple way, we assumed that the net effect is to reduce the actual infection kinetic rate *µ* by some constant factor. This is given in the table as *µ* after *L*1. Note that the transition occurs instantaneously in our model hence the sharp drop in the total infected at that time. Other parameters are tuned to account for the actual variability of Δ*I**new* (but not its absolute value) and official number of deaths (*D*(*t*) + *d*(*t*)). The delay for recovery or death processes has been estimated from the measurements of hospitalisation recovery in a Country. For instance, Fig. 2 shows the estimation of the recovery time-delay for Belgium: it corresponds to the *time-interval* between the peak of the new admission and the peak of the recovered people from hospitals. Such a procedure has been adopted for estimating the recovery and death time-delays also for France and Germany. View this table: [Table 1:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/T1) Table 1: List of the Parameters ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F2) Figure 2: Estimation of the time-delay. The time-delays have been estimated by considering the time-interval between the peak of the new admission and the peak of the recovered people from hospitals. This figure corresponds to the Belgian case. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F3) Figure 3: Theoretical solutions for infectious (I), cumulative number of recovered people (R) and deaths (D) for Belgium. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F4) Figure 4: Theoretical solutions for hospitalised infectious (Ih), total recovered (rh) and total deceased (dh) people, previously hospitalised, for Belgium. ### Belgian Case Figs (3) refer to the Belgian case. In particular, Fig (3) shows the solutions of our model for the infectious (*I*), total recovered (*R*) and total deceased (*D*) people. Fig. (4) illustrates the theoretical solutions for hospitalised infectious (*I**h*), the total recovered (*r**h*) and total deceased (*d**h*) people previously hospitalised. Figs (5) and (6) shows the comparison between the theoretical predictions for Δ*I**new*(*t*) and deaths and real data for Belgium (according to the database *Sciensano*). Notice in Fig. 5 the prediction of the *second wave of infection by SARS-CoV-2* ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F5) Figure 5: Comparison between the the-oretical prediction for ΔINew with real data provided by the data base Sciensano, for Belgium. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F6) Figure 6: Comparison between the theoretical solution of our model for Deaths with real data provided by the database Sciensano, for Belgium. ### French Case Figs (7) and (8) shows the comparison between the theoretical predictions for Δ*I**new*(*t*) and deaths and real data for Belgium (according to the database *Santé Publique France*). Notice in Fig. 7 the prediction of the *second wave of infection by SARS-CoV-2* ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F7) Figure 7: Comparison between the theoretical prediction for ΔINew with real data provided by the data base Santé Publique France, for France. ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F8) Figure 8: Comparison between the theoretical solution of our model for Deaths with real data provided by the database Santé Publique France, for France. ### German Case Figs (9) and (10) shows the comparison between the theoretical predictions for Δ*I**new*(*t*) and deaths and real data for Belgium (according to the database *(Robert Koch Institut). Country data from Worldbank*.*org*). Notice in Fig. 9 the prediction of the *second wave of infection by SARS-CoV-2* ![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F9) Figure 9: Comparison between the theoretical prediction for ΔINew with real data provided by the data base (Robert Koch Institut. Country data from Worldbank.org, for Germany. ![Figure 10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249797/F10.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249797/F10) Figure 10: Comparison between the theoretical solution of our model for Deaths with real data provided by the database (Robert Koch Institut. Country data from Worldbank.org, for Germany. ## 7 Conclusions and Perspecives We showed that our model is able to produce predictions not only on the first but also on the second or even the third waves of SARS-CoV2 infections. The theoretical predictions are in line with the official number of cases with minimal parameter fitting. We discussed the strengths and limitations of the proposed model regarding the long-term predictions and, above all, the duration of how long the lockdown and the quarantine measures should be taken in force in order to limit as much as possible the intensities of subsequent SARS-CoV-2 infection waves. This task has been carried out by taking into account the theoretical results recently appeared in literature [7] and without neglecting the delay in the reactions steps. Our model has been applied in two different situations: the spreading of the Coronavirus in a small Country (Belgium) and in big Countries (France and Germany). It is worth noting the *degree of the flexibility* of our model. For example, let us suppose that we need to set up a model able to distinguish old population (over 65 year old) from the young one (with age not exceeding 35 years), by assuming that the older population is twice as likely to get infected by Coronavirus with respect to the younger one. In this case, it is just sufficient to replace the scheme ![Graphic][25] with the scheme ![Formula][26] with *S**Y* and *S**o* denoting the *susceptible young people* and the *susceptible old people*, respectively. Another example could be the following. Let us suppose that we need to distinguish two class of infected individuals: 1. infected people (denoted by *I*1) able to transmit the Coronavirus to susceptible according to the (standard) scheme *I*1 + *S* → 2*I*; 2. Infected people (denoted by *I*2) having the capacity to transmit the virus, say, 7 times higher with respect to the category **1)**. In this case, the corresponding scheme reads: ![Formula][27] It is then easy to write the ordinary differential equations associated to schemes (23) and (24). Let us now consider another aspect of the model. In the Subsection (2.2), we have introduced scheme (2) that models the lockdown measures. As mentioned, such measures are imposed by national governments to all susceptible population. However, we can also take into consideration the hypothesis that these measures are not rigorously respected by the population and this for various reasons: neglect of the problem, depression due to prolonged isolation, lack of confidence in the measures adopted by the Government, desire to attend parties with friends and relatives, etc. Scheme (2) still adapts to describe these kind of situations with the trick of replacing Fig. 1 with a curve that models the *emotional behaviour* of susceptible people. The O.D.E.s read ![Formula][28] where *E* stands for *Emotional*. Finally, we mention that in ref. [14] we have incorporated real data into a stochastic model. The goal is to obtain a comparative analysis against the deterministic one, in order to use the new theoretical results to predict the number of new cases of infected people and to propose possible changes to the measures of isolation. ## Data Availability All data referred to this manuscript are available. ## Footnotes * Emails: gsonnino{at}ulb.ac.be, peeters.philippe{at}gmail.com, pnardon{at}ulb.ac.be * 1 Notice that the first *reaction* in the scheme Eq. (5) is the dynamic equation for the total recovered people adopted in the SIRD-model [13]. * 2 We consider that the SARS-CoV-12 has just appeared for the first time. So, we do not consider the asymptomatic people who are immune to the virus without any medical treatment. * 3 Our model is inspired by Michaelis-Menten’s enzyme-substrate reaction. Of course, the reverse *MM reaction* has no sense in our case and, consequently, the *kinetic constant* is equal to zero. * Received January 14, 2021. * Revision received January 14, 2021. * Accepted January 15, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].Cacciapaglia G., Cot C., and Sannino F., Second wave COVID-19 pandemics in Europe: a temporal playbook, Scientific Reports, 10, Article number: 15514, 2020. 2. [2].Bailey, Norman T. J., (1975), The mathematical theory of infectious diseases and its applications (2nd ed.). London: Griffin. ISBN 0-85264-231-8. 3. [3].Sonia A. and Nunn, C. (2006). Infectious diseases in primates: behavior, ecology and evolution. Oxford Series in Ecology and Evolution. Oxford [Ox- fordshire]: Oxford University Press. ISBN 0-19-856585-2. 4. [4].Vynnycky E. and White R. G., An Introduction to Infectious Disease Modelling. Oxford University Press. ISBN 978-0-19-856576-5 (eds. 2010). 5. [5].Gleick J., Chaos: The Making of a New Science, Viking Press, New York (1987). 6. [6].Coullet P., The Covid-19 epidemic as a simple dynamical system, Comptes Rendus Mécanique, 0000, 1, no , p. 000-000 DOI unassigned yet, Under Review (2020). 7. [7].Sonnino G. and Nardone P., Dynamics of the COVID-19 - Comparison between the Theoretical Predictions and the Real Data, and Predictions about Returning to Normal Life, Annals of Clinical and Medical Case Reports, 4 (Issue 9), pp. 1–21, 2020. ISSN 2639-8109. doi: [http://dx.doi.org/10.47829/ACMCR.2020.4902](http://dx.doi.org/10.47829/ACMCR.2020.4902). 8. [8].Sonnino G, The COVID-19 - The Infectious Disease Caused by the Latest Discovered Coronavirus (SARS-CoV-2), European Commission, ARES(2020)1530456, (March 2020). 9. [9].Michaelis L. and Menten M.L., Die Kinetik der Invertinwirkung. Biochem Z. 49:333–369, 1913 [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000200780700021&link_type=ISI) Henri, Victor (1903). Lois Générales de l’Action des Diastases. Paris: Hermann. Victor Henri. Whonamedit ?. Retrieved 24 May 2011. 10. [10].Srinivasan, Bharath (2020-10-08), Explicit treatment of non Michaelis-Menten and atypical kinetics in early drug discovery. ChemMed-Chem. Retrieved 2020-11-09. doi:10.20944/preprints202010.0179.v1. PMID 33231926. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20944/preprints202010.0179.v1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33231926&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249797.atom) 11. [11].Srinivasan, Bharath (2020-09-27), Words of advice: teaching enzyme kinetics. The FEBS Journal. doi:10.1111/febs.15537. ISSN 1742-464X. PMID 32981225. Substrate Concentration (Introduction to Enzymes). [www.worthington-biochem.com](http://www.worthington-biochem.com). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/febs.15537&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32981225&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249797.atom) 12. [12].Sonnino G., Peeters P., and Mardone P, Modelling the Coronavirus Second Wave in Presence of the Lockdown and Quarantine Measures, International Conference on Complex Systems 2020 (CCS2020)-December 2020 (page 278). DOI reference: [https://doi.org/10.5281/zenodo.4419178](https://doi.org/10.5281/zenodo.4419178) 13. [13].Kermack W. O. and McKendrick, A. G., A Contribution to the Mathematical Theory of Epidemics. Proceedings of the Royal Society A. 115 (772): 700–721 (1927). doi:10.1098/rspa.1927.0118 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1927.0118&link_type=DOI) 14. [14].Sonnino G., Mora F., and Nardone P., A Stochastic Compartmental Model for COVID-19, submitted to publication to the review CHAOS, 2020. 15. [15].[https://www.sciensano.be/en/covid-19-data](https://www.sciensano.be/en/covid-19-data) 16. [16].[https://www.santepubliquefrance.fr/](https://www.santepubliquefrance.fr/) 17. [17].[https://dc-covid.site.ined.fr/en/data/germany/](https://dc-covid.site.ined.fr/en/data/germany/) 18. [18].[https://ourworldindata.org/covid-models](https://ourworldindata.org/covid-models) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-6.gif [7]: /embed/graphic-7.gif [8]: /embed/graphic-8.gif [9]: /embed/graphic-9.gif [10]: /embed/graphic-10.gif [11]: /embed/graphic-11.gif [12]: /embed/graphic-12.gif [13]: /embed/graphic-13.gif [14]: /embed/graphic-14.gif [15]: /embed/graphic-15.gif [16]: /embed/graphic-16.gif [17]: /embed/graphic-17.gif [18]: /embed/graphic-18.gif [19]: /embed/graphic-19.gif [20]: /embed/graphic-20.gif [21]: /embed/graphic-21.gif [22]: /embed/graphic-22.gif [23]: /embed/graphic-23.gif [24]: /embed/inline-graphic-2.gif [25]: /embed/inline-graphic-3.gif [26]: /embed/graphic-34.gif [27]: /embed/graphic-35.gif [28]: /embed/graphic-36.gif