returning to normal life - arxiv · giorgio sonnino and pasquale nardone facult e des sciences...

Dynamics of the COVID-19

Comparison between the Theoretical Predictions

and the Real Data, and Predictions about

Returning to Normal Life

Giorgio SONNINO and Pasquale NARDONEFaculté des Sciences

Université Libre de Bruxelles (ULB)Boulevard de Triomphe, Campus Plain CP 231

1050 Brussels, BelgiumEmails: [email protected] and [email protected]

Manuscript published inAnnals of Clinical and Medical Case Reports, 4 (9), 2020

ISSN 2639-8109

Abstract

A new coronavirus disease, called COVID-19, appeared in the Chi-nese region of Wuhan at the end of last year; since then the virus spreadto other countries, including most of Europe. We propose a differentialequation governing the evolution of the COVID-19. This dynamic equa-tion also describes the evolution of the number of infected people for 13common respiratory viruses (including the SARS-CoV-2). We validateour theoretical predictions with experimental data for Italy, Belgium andLuxembourg, and compare them with the predictions of the logistic model.We find that our predictions are in good agreement with the real worldsince the beginning of the appearance of the COVID-19; this is not thecase for the logistic model that only applies to the first days. The secondpart of the work is devoted to modelling the descending phase, i.e. thedecrease of the number of people tested positive for COVID-19. Also inthis case, we propose a new set of dynamic differential equations that wesolved numerically. We use our differential equations parametrised withexperimental data to make several predictions, such as the date whenItaly, Belgium, and Luxembourg will reach a peak number of SARS-CoV-2 infected people. The descending curves provide valuable informationsuch as the duration of the COVID-19 epidemic in a given Country andtherefore when it will be possible to return to normal life. The study ofthe the dynamics of COVID-19 when the population have been subject

1

arX

iv:2

003.

1354

0v7

[q-

bio.

PE]

12

Sep

2020

to less restrictive measures is beyond the scope of this work and it will bematter of future works.

Key words: Mathematical model; COVID-19; Pneumonia

1 Introduction

Viral infections usually affect the upper or lower respiratory tract. Althoughrespiratory infections can be classified according to the causative agent (e.g.the flu), they are mostly clinically classified according to the type of syndrome(e.g., common cold, bronchiolitis, laryngo-tracheo-bronchitis acute, pneumonia).Although pathogens typically cause characteristic clinical manifestations (e.g.,rhinovirus causes the common cold, respiratory syncytial virus [RSV] usuallycauses bronchiolitis), they can all cause many of the most common respiratorysyndromes. The severity of viral respiratory disease is highly variable; seriousillness is more frequent in elderly patients and young children. Morbidity caneither directly result from the infecting agent, or may be indirect. The lattercase can be due to the exacerbation of an underlying cardiopulmonary disease,or a bacterial superinfection of the lung, paranasal sinuses, or middle ear. Themain motivation of this work is to verify, by making theoretical predictions, thatpolitical decisions are truly effective to minimise the number of infected peoplein order to (i) not overload local health services (such as hospitals), and to (ii)gain time to allow research institutes to deliver vaccines or the anti-virals.Tables (1) and (3) respectively provide the experimental data for Italy [1] and forBelgium [2], [3]. They show the number of active people (i.e., people currentlyinfected by SARS-CoV-2), the recovered people, and deaths for COVID-19.We start our theoretical analysis by introducing the definition of the basic re-production number of an infection R0, defined as the number of infected peoplederived from a first case in a population where all the others are susceptible. So,it is not possible to modify R0, in any case, but it is possible to get a differenteffective R1. This parameter is strictly linked to the replication time of a virus,indicated with µi, defined as the time interval after which the number of infectedpeople has increased by R0 times. Fig. 1 schematically represents the diffusiondynamics of the virus. By indicating with N the number of infected people,after n steps we get2:

N = Rn0 (1)

Of course, after n steps, the elapsed time is t =∑ni=1 µi and, if there are M

1More rigourously, in epidemiology, the basic reproduction number of an infection, R0,is the expected number of cases directly generated by one case in a population where allindividuals are susceptible to infection in absence of any deliberate intervention in diseasetransmission (see, for example, [4]).

2In this Section we shall follow the definitions and the expressions reported in standardbooks or thesis dissertation such as, for example, [5] or [6].

2

n1

n2

n3

Figure 1: Schematic dynamics of respiratory virus in the absence of the lock-down measures In this graphics, for illustrative purpose only, we set R0 = 3.However, for SARS-CoV-2, the value of R0 is 2 even at the beginning of theoutbreak in China and Italy. After a period of time µ1, an infected individualcan infect R0 other individuals. In turn, after a period µ2, each of these newlyinfected individuals can infect other R0 people, and so on. After n steps theelapsed time is t =

∑ni=1 µi.

3

outbreaks of infectious viruses, Eq. (1) can be cast into the form3

N = MRt/µ0 (2)

with µ ≡ 1/n∑ni=1 µi. Note that the two parameters R0 and µ are not inde-

pendent (see, for example, [7], [8] or [9])4. It is more convenient to work in theEuler base e rather than in base R0; in the Euler base Eq. (2) provides the lawof growth of a Malthusian population [5].

N = M exp(t/τ) where τ =µ

log(R0)(3)

In literature, τ is referred to as the characteristic time of the exponential trend.So, in the absence of containment measures the number of infected people followsthe exponential law (3). Let us now analyse Eq. (3) in more dept. We havethree possible scenarios:

1. R0 > 1 (as is the current world’s situation). For Italy, for example, beforethe adoption of (severe) containment measures, the value of τ was aboutτ ∼ 3.8 days (and µ ∼ 2.6 days). In this case the number of the infectedpeople increases exponentially.

2. R0 = 1 If the infection-capacity of the virus is of the type one-to-one(i.e., a person infected by SARS-CoV-2 can in turn infects only anotherperson), we get the stationary situation corresponding to N = 1. Thissituation is referred to as the latent situation i.e., the virus is still presentbut does not spread. In this limit case, the SARS-CoV-2 is substantiallyineffective. Scenarios (1) and (2) are illustrated in Fig. 2.

3. 0 < R0 < 1. We may also imagine that the capacity of infection ofSARS-CoV-2 is less than 1. This means that the virus is no longer ableto be spread (e.g., thanks to protective measures, or to the productionof vaccines and anti-virals, or because people who overcame the diseasebecame immune. In this case, the value of τ is negative and the numberof infected people decreases ever time. That is, the infection eventuallydisappears. The rate of decrease of the number of infected people dependson the value of τ . This scenario is depicted in Fig. 3.

1.1 Comparison with the Real Data for COVID-19 beforethe Lockdown Measures

It is understood that the main objective of the lockdown measures established bymost European governments and health organisations is to reduce the ability

3Actually, Eq. (2) applies only if the M outbreaks of the virus are exactly at the same

conditions. In general, the correct expression reads N =∑Mi=1R

t/µ̃i0 , with µ̃i indicating the

replication time of the virus for the i-th outbreak.4In ref. [7], the doubling time is used to calculate R0, by means of the equation R0 =

1 + (γ + ρ)log(2)/µ where γ is the duration of the incubation period, ρ is the duration of thesymptomatic period, and µ is the doubling time (see [7]). In this respect, we would also liketo mention another excellent work recently produced by G. Steinbrecher [9].

4

Days

N

“Latent” situation

1

R0 > 1

0 2 4 6 8 10 12 14

10

20

30

40

R0=1

Figure 2: Situation before the lock-down measures. Number of infectedpeople corresponding to the exponen-tial law. The red line represents thecase R0 > 1, such as the situationbefore the adoption of lockdown mea-sures. The black line corresponds tothe case R0 = 1, the latent situationin which the virus is substantially in-effective.

0 2 4 6 8 10 12 14

0.2

0.4

0.6

0.8

1.0

N

Days

R < 1 0

Figure 3: Number of infected peo-ple corresponding to the exponentiallaw. The red line represents the caseR0 < 1. In this situation the numberof infected people decreases exponen-tially and the virus disappears after afew weeks.

of a virus to spread. From a mathematical point of view, we would like tohave R0 = 1 (or, better, R0 < 1), in Eq. (3) instead of R0 > 1. In practicalterms, this means reducing the frequency of all involuntary contacts with a largenumber of people, reducing unnecessary movements to avoid encounters, and toprolong the closure of schools. Although these measures cannot prevent thespread of the infection in the long term, they can reduce the number of newinfections daily. This has the benefit of leaving room for seriously-ill patients byavoiding to overload the healthcare system. We can easily realise what are theconsequences if the lockdown measures are not set up. To make a comparisonbetween the theoretical predictions and the experimental data in absence oflockdown measures, we have to consider the correct reference period. Morespecifically, we saw that the number of positive cases grows in the course oftime by following the law (3). Hence, at the reference time t0, the number ofpeople infected by the virus is

N0 = M exp(t0/τ) (4)

After a period of time, say t, Eq. (3) reads

N = M exp(t/τ) (5)

Hence,N = N0 exp((t− t0)/τ) (6)

5

Eq. (6) is the equation that we use for comparing the mathematical predictionswith experimental data during the initial phase where the spread of SARS-CoV-2, causing the COVID-19, follows the exponential law, and (t− t0) is ourreference period. For the case of COVID-19 we get (see, for example, [5] or [6])

• All infectious outbreaks are exactly at the same conditions. So, Eq. (2)applies;

• R0 = 2;

• All the µi are equal with each other: µi = const = µ (see also [5] and [7]).

In this case, µ is referred to as the doubling time. So, the doubling time isthe amount of time it takes for a given quantity to double in size or value at aconstant growth rate [8]. If we do not apply the locking measures, the evolutionin the course of time of the number of infected people is best approximated byan exponential curve with R = 2, even though we have to stress that R0 isonly associated with the beginning of the epidemic and, with certain approxi-mations, with the early stages, but not beyond. Fig. 4 and Fig. 5 respectivelyshow the comparison between the theoretical predictions and the experimentaldata for Italy and Belgium before the lockdown measures. We get τ ' 3.8days and µ ' 2.6 days for Italy, and τ ' 5.2 days and µ ' 3.7 days forBelgium. We conclude this Introduction by mentioning that there are severalmethods currently proposed in Literature to derive by mathematical models,the value of R0. For example, in ref. [9], we have a short numerical code, writ-ten in R-programming language for statistical computing and graphics, ableto compute the estimated R0 values for the following 17 infectious diseases:Chickenpox (varicella) (Transmission: Aerosol), Common cold (Transmission:Respiratory Droplets), COVID-19 (Transmission: Respiratory Droplets), Diph-theria (Transmission: Saliva), Ebola - 2014 Ebola outbreak (Transmission::Body fluids, HIV/AIDS (Transmission: Body fluids), Influenza - 1918 pandemicstrain (Transmission:: Respiratory Droplets), Influenza - 2009 pandemic strain(Transmission: Respiratory Droplets, Influenza - seasonal strains (Transmis-sion: Respiratory Droplets), Measles (Transmission: Aerosol), MERS (Trans-mission: Respiratory Droplets), Mumps (Transmission: Respiratory Droplets),Pertussis (Transmission: Respiratory Droplets), Polio (Transmission: Fecaloral route), Rubella (Transmission: Respiratory Droplets), SARS (Transmis-sion: Respiratory Droplets), Smallpox (Transmission: Respiratory Droplets).However, this task is particularly problematic if there are intermediate vectorsbetween hosts, such as malaria.This manuscript is organised as follows. Section (2) determines the dynamicdifferential equation for the COVID-19; Section (2.3) compares the theoreticalpredictions and experimental data for Italy and Belgium. The differential equa-tions providing the evolution of the decrease of the number of people testedpositive for COVID-19 can be found in Section (3); Section (4) concludes. Thecomparison between the theoretical predictions of our model and experimentaldata for Luxembourg, as well as the solution of the differential equations in the

6

From 01/02/2020 to 11/03/2020

Italian situation before the lockdown measures

Days

Lockdown measures: 11/03/2020

Exponential phase

N

26 28 30 32 34 36

2000

4000

6000

8000

Figure 4: Number of infected peo-ple in Italy on the 10th of March2020 (before the adoption of lock-down measures). The blue line cor-responds to the theoretical predic-tions and the black dots correspondto experimental data. The valuesof the parameters τIT and µIT areτIT ' 3.8 days and µIT ' 2.6 days,respectively.

5 10 15 20 25 30

2000

4000

6000

8000

10000

From 29/02/2020 to 29/03/2020

Belgian situation before the

lockdown measures

DaysLockdown measures:

16/03/2020

Exponential phaseN

Figure 5: Exponential phase in Bel-gium. The lockdown measure havebeen adopted on the 16the of March2020 (however, initially not so strictas in Italy).The red line correspondsto the theoretical predictions andthe black dots correspond to exper-imental data. The values of the pa-rameters τBE and µBE are τBE '5.3 days and µBE ' 3.7 days, re-spectively.

descending phase for Luxembourg are reported in Appendix. We may objectthat we are dealing data from countries which have passed the peak of infection,such as South Korea, Iceland or Austria etc. The situation in other Countries,which may have adopted other political decisions about the application of thelockdown measures, may be subject of future works. However, we would like tomention that several authors are currently applying our model to other Coun-tries. In this regard, we have received their pre-prints such as the work cited inRef. [11]. More specifically, we have received a message where our model hasbeen used, with success, to analyse data from UK, USA, NY City, Spain, andMumbai City. We stress the fact that this manuscript deals with the spread ofSARS-CoV-2 until May 16, 2020, as the objective of this work is to study theeffect of the strict lockdown measures. After May 2020, these measures havebeen modified by the various Governments, which have decided to introduceless heavy and much less restrictive lockdown measures5.

5The study of the the dynamics of COVID-19 when the population has been subject to lessrestrictive measures is beyond the scope of this work and it will be matter of future works.

7

2 Modelling the COVID-19 - Virus’ growth

The objective of this section is to determine the coefficients of the evolutionarydifferential equation for the COVID-19 (see the forthcoming Eq. (13)). Wealso determine the generic analytical expression for the time-dependent numberof infected people through fitting techniques validated by the χ2 tests. Thisexpression is proposed after having previously analysed 12 respiratory infectiousdiseases caused by viruses [10], in addition of being solution to the Richard’sdifferential equation.

2.1 General background

Letting N represent population size and t represent time, the Logistic modelmodel is formalised by the differential equation below:

dN

dt= αN

(1− N

K

)(7)

where α > 0 defines the grow rate and K > 0 is the carrying capacity. In thisequation, the early, unimpeded growth rate is modelled by the first term +αN .The value of the rate α represents the proportional increase of the populationN in one unit of time. Later, if the system is closed (i.e. the system is isolatedand, hence, not in contact with a reservoir allowing the system to exchangeindividuals), as the population grows the modulus of the second term, −αN2/K,becomes almost as large as the first, until to saturating the exponential growth.This antagonistic effect is called the bottleneck, and is modelled by the value ofthe parameter K. The competition diminishes the combined growth rate, untilthe value of N ceases to grow (this is called maturity of the population). Thesolution of Eq. (7) is

N(t) =K

1 +B exp (−t/τ)(8)

where B > 0 is a constant related to the value of N(0). It is more convenientto rewrite Eq. (8) in terms of the initial Logistic time t0L defined as

t0L = τ logB (9)

So, Eq. (8) may be cast into the form

N(t) =K

1 + exp (−(t− t0L)/τ)(10)

where τ is linked to the steepness of the curve. Since the environmental condi-tions influence the carrying capacity, as a consequence it can be time-varying,with K(t) > 0, leading to the following mathematical model (see, for example,[12]):

dN

dt= αN

(1− N

K(t)

)(11)

8

More generally, the growth modelling is well described by Richards’ differentialequation (RDE) [13]

dN

dt= αN

(1−

( NK(t)

)ν)(12)

where ν > 0 affects near which asymptote maximum growth occurs. The phe-nomenological logistic function is used to model the evolution of the COVID-19pandemic in different Countries. The logistic model is mainly used in epidemi-ology and provides insights into the transmission dynamics of the virus. Wenote, however, to evaluate the dynamics of transmission of SARS-CoV-2, morerefined models are needed, which take into account specific measures adoptedin each Country [14]. So, let us suppose that the Government decides to adoptthe lockdown measures. After the application of the lockdown measures theequation may be revised to be

dN

dt= αN

(1−

( NK(t)

)ν)− c(t)N (13)

where c(t) takes into account the degree of effectiveness of the lockdown mea-sures.

2.2 Determination of the Carrying Capacity and the Lock-down Coefficient for the COVID-19

According to ref. [15]6 Respiratory viruses remain quiet for months, inactive butviable, within living cells. Then suddenly they activate and become virulent asthey say, the infectious capacity grows to a maximum, after which it decreases.The time duration is about of 2 or 3 months. So we can expect that the epidemicwill soon die out in Italy too. So, there is no valid reason to think that this coro-navirus behaves differently from others [15]. The present work starts from thefollowing hypothesis: the SARS-CoV-2 behaves like the other viruses that causerespiratory diseases. As a consequence, for the COVID-19 case, functions K(t)and c(t) are determined by performing several fittings on the growth rate-trendsof infection capacity of the viruses that mainly affect the respiratory system.More specifically, we considered the following 13 different diseases caused by12 different viruses: Whooping Cough (Pertussis), Swine Flu (H1N1), Bird Flu(Avian Flu H5N1), Enterovirus, Flu in Children, Flu in Adults, Bacterial Pneu-monia, Viral Pneumonia, Bronchitis, Common Cold (Head Cold), Severe acuterespiratory syndrome (SARS), and MERS (Middle East Respiratory Syndrome).In all the examined cases, we took into account the fact that the therapy-induceddeath rate is greater than the baseline proliferation rate, then there is the eradi-cation of the disease. In other words, for the above-mentioned cases the function

6Prof. Roberto Ronchetti is currently working at the Pediatric Clinic of La Sapienza ofthe University of Rome, at the Policlinico Umberto I and at the S. Andrea 24 MarchHospitalwhere he helped to found, dealt with childhood respiratory diseases, and studied bronchiolitisin particular. In these days he has studied (with his collaborators) the data available onSARS-CoV-2 in China, in South Korea and now in Italy.

9

c(t) in Eq. (13) represents the therapy-induced death rate [16]-[24]. Of course,this is an oversimplified model of both the growth and the therapy (e.g., it doesnot take into account the phenomenon of clonal resistance). We empiricallynoticed (according to the χ2 test) that all these viruses have in common thesame growth rate-trend of infected people (of course, each of these behaviourshave their own growth rate parameters). In particular, we get that the trendsof the total number of infected people by respiratory viruses (indicated withN), subject to the therapy-induced death rate, versus time satisfy the followingO.D.E. [10]

dN

dt̂= αN

(1− N

KN

)−(αt̂2 − 1

t̂

)N with t̂ > 1/α1/2 (14)

where we have introduced the dimensionless time t̂ ≡ t/t0. The coefficient

c(t) ≡(αt̂2 − 1

t̂

)(with t̂ > 1/α1/2) (15)

is referred to as the average therapy-induced death rate. In our case the term−c(t)N in the dynamic equation represents the lockdown measures. The lock-down is mainly based on the isolation of the susceptible individuals, eventuallywith the removal of infected people by hospitalisation7. In the idealised case,for αt̂2 > 1, c(t) may be modelled as a linear function of t̂, by getting

c(t) ' αt̂ (16)

As for the epidemiological explanations relating to the various modelling ofc(t) (constant, linear in time etc.), we refer the reader to the well-known andextensive literature on the subject (see, for example, Ref [25] or to the referencescited in [26]). Here, we limit ourselves to give a very intuitive explanation on thephysical meaning of this contribution. Immediately after the lockdown measureshave been adopted, i.e. during the very first initial phase, we expect that c(t̂)is practically constant in time, as these measures have not yet been able to acteffectively. However, after a short period of time, the positive effects of thelockdown measures become increasingly efficient and it is intuitive to expect alinear growth of c(t̂) in time. More specifically, we expect that, after a shorttransition period the coefficient c(t) starts to grow linearly in time. Successively,at the leading order, c(t̂) will be equal in magnitude to the coefficient of thelinear term (in order to balance the growth rate induced by the linear term).This because the lockdown measures will be able to satuarte the exponentialgrowth. Briefly, we expect that the O.D.E. governing the dynamics of the SARS-CoV-2 is of the form (14) where c(t̂) ' αt̂ for t̂ > 1/α1/2. Indeed, for values of

7It is worth mentioning that initially England did not adopt any lockdown measures believ-ing that the British system be a closed system. Basically, it was believed that the system begoverned by a simply logistic equation with a carrying capacity decreasing in time. However,contrary to the expectations, in England the carrying capacity did not decrease in time. Thisinduced the British government to adopt the lockdown measures.

10

time t̂ ' 1, the lockdown term in Eq. (14) is able to balance the exponentialgrow, which is in agreement with our intuitive expectations.From Eq. (14) we get that the time derivative of N vanishes for N = Ns, withNs satisfying the equation

NsKN

=αt̂− αt̂2 + 1

αt̂> 0 (17)

By taking into account the inequality reported in Eq. (15), we get that theO.D.E. (14) is valid in the range

1

α1/2< t̂ <

(1 + α/4)1/2

α1/2+

1

2(18)

Parameters KN and α depend on the virus in question and on the externalconditions (e.g. in our case, the lockdown measures) to which the population issubject. In Eq. (14), the term −N2/KN is the term that tends to saturate theexponential growth. KN is constant (or decreases) in the course of time sincethe non-linear contribution becomes more and more important until saturatingthe exponential growth. In our model, the carrying capacity is kept constant.For large values of the carrying capacity, the solutions of Eq. (14) reach theplateau at the time t̂Max given, approximatively, by the expression

t̂Max ≡tMaxt0' 1α1/2

+1

2(19)

Notice that α is linked to µ. Indeed, as shown in Section 1, during the expo-nential period the COVID-19 grows according to the law (see Eq. (6)):

dN

dt̂' τ̂−1N where τ̂ ≡ τ

t0(20)

Hence, we get

α ' 1τ̂

=log(R0)

2µ̂where µ̂ ≡ µ

t0(21)

We conclude this Section by mentioning that we may easily check that, forsystems having a large carrying capacity, the solution of Eq. (14) is well ap-proximated by the expression

N ' At exp(−(t− t0)2/σ) with σ = 2t20/α (22)

The values of parameters A, t0 and σ depend on the virus in question. It isconvenient to re-write Eq. (22) in dimensionless form

N ' Ât̂ exp(−(t̂− 1)2/σ̂) where Â ≡ At0 ; σ̂ ≡σ

t20(23)

To sum up, according to our model for COVID-19, the ascending behaviour ofthe total cases (i.e., the number of of positive cases plus the cumulative number

11

25 30 35 40 45 50

3

4

5

6

7

8

9n

Days

n0

Italian situation before the lockdown measures

Italian during the transient period after the lockdown measures

From 25/02/2020

to 10/03/2020

From 11/03/2020

to 24/03/2020

LMt

Figure 6: Italian transient period(from the 10th of March 2020 to the24th of March 2020). During thisperiod, the doubling time µ oscil-lates over time. µ0 indicates the(constant) doubling time during theexponential period (for Italy µ0 '2.6 days).

0 5 10 15 20 25

2

3

4

5

6

7

Days

nBelgian situation before the

lockdown measures Belgian during the transient period

after the lockdown measures

n0

From 01/03/2020

to 16/03/2020

From 17/03/2020

to 29/03/2020

LMt

Figure 7: Belgian transient period(from the 17th of March 2020 tothe 29th of March 2020). Duringthis period, the doubling time µoscillates over time. µ0 indicatesthe (constant) doubling time duringthe exponential period (for Belgiumµ0 ' 3.7 days).

of recovered people plus the cumulative number of deaths) is given by the solutionof Eq. (14) for 1/α1/2 ≤ t̂ ≤ t̂Max.Notice that the determination of the O.D.E. (14) is of a fundamental importanceif we wish also to describe the stochastic process (and the associated Fokker-Planck equation) where a white noise is added to this O.D.E. According to theliterature nomenclature, we refer to the differential equation (14) as COVID-19dynamic model8.

2.3 Comparison between the Theoretical Predictions andExperimental Data

For Italy and Belgium one observes two distinct phases related to the dynamicsof the COVID-19, which we classify as before the adoption of the lockdownmeasures and after some days after the adoption of the lockdown measures.The question therefore naturally arises, of whether these two types of regimeare separated by a well-defined transition. We shall see that this is indeed thecase. We may identify three different periods, which may be classified as follows:

1. The exponential period. As seen in Section 1, before the adoption of lock-down measures, the exponential trend is the intrinsic behaviour of the grow

8Viral dynamics is a field of applied mathematics concerned with describing the progressionof viral infections within a host organism (see, for example, [27].)

12

rate of the COVID-19. In this period the doubling time µ is a constantparameter versus time.

2. The transient period. The transient period starts after having applied thesevere lockdown measures. In this period we observe a sort of oscillations(or fluctuations) of µ versus time. In this case the time variation of µ(t)reflects the behaviour of the time effective reproduction number, R(t), de-fined as the number of cases generated in the current state of a population,which does not have to be the uninfected state. Fig. 6 and Fig. 7 show thebehaviour of the parameter µ versus time for Italy and Belgium, respec-tively. The transient period ends when the last step of the exponentialtrend fits real data as good as the linear trend9.

3. The bell-shaped period (or the post-transient period). In the bell-shapedperiod parameter µ is a (typical) function of time obtained by usingEq. (14). Several theoretical models can be used to study the post-transient period (e.g., by using Gompertz’s law [28]). Here, we use twomathematical models: the solution of the differential equation (14) andthe logistic model (see, for example, Ref. [16]), and we compare these twotheoretical models with real data for Italy and Belgium.

Figures (8), (9), and (16) (see Appendix) compare the predictions of our model(blue lines) against the logistic model (red lines) for Italy, Belgium, and Lux-embourg, respectively. Notice that the number of free parameters of these twomodels are exactly the same, since τ and α cannot be chosen arbitrarily. Morespecifically,

a) The logistic model possesses two free parameters : K and t0L. Notice thatparameter τ is not free since it is linked to the doubling time µ;

b) Also our model possesses two free parameter: KN and t0. Notice thatparameter α is linked to the doubling time µ (see Eqs (14) and (20)).

Fig. 8 and Fig. 9 compare the theoretical predictions, with the experimentaldata for Italy and Belgium updated to the 15th of May 2020. The values of theparameters τ , KN , and t0 for Eq. (14) and the parameters τ , t0 and K for thelogistic function are reported in the figure captions. As we can see, for bothCountries Eq. (14) fits well all the real data from the initial days, while thelogistic model applies only to the first data. The curves reach the plateau atthe time tMax given by Eq. (19). By inserting the values of the parameters, weget

tMaxIT ' 80 days and tMaxBE ' 60 days (24)

corresponding to tMaxIT ' 21 April 2020 and tMaxBE ' 2 May 2020 for Italyand Belgium, respectively.

9A numerical condition may be established by using the χ2 test: the fittings of the twotrends are considered both good if, for example, for both trends, the χ2-tests get values ≥ 0.9.

13

0 20 40 60 80 100

50000

100000

150000

200000

From 01/02/2020

to 11/03/2020

From 11/03/2020 to 15/05/2020

Days

Italian situation before the

lockdown measures

Blue full line: COVID-19 dynamic model

Italian situation after 65 days of

lockdown measures

Red full line: Logistic model

N

Figure 8: Situation in Italy on 15 May 2020—before, and 65 days after, theadoption of lockdown measures. The black dots correspond to experimentaldata. The red dotted line corresponds to the situation in Italy before the adop-tion of the lockdown measures. The blue and the red solid lines correspond tothe theoretical predictions for Italy according to the solution of Eq. (14) andthe logistic model, respectively. Solution of Eq. (14) fits well all the experimen-tal data from the initial days (i.e., from the 1st of February 2020), while thelogistic model applies only to the first days. The values of the parameters ofEq. (14) and the logistic function (10) are: τIT ' 3.8 days (µIT = 2.6 days),KITN ' 355250, and t0IT ' 72.5 days for Eq. (14), and τIT ' 3.8 days (µIT = 2.6days), KIT = 225000, t0LIT = 53 days for the Logistic function.

3 Modelling the COVID-19 - The DescendingPhase

Here, for the descending phase is intended the phase where the number of thepositive cases starts to decrease10. So, our model cannot be used for describingalso the descending phase since Nt is the number of the total cases and, duringthe descending phase, Nt tends to reach the plateau. The objective of this sectionis to determine the trend of the curve of positive people during the descendingphase. This task is accomplished by establishing the appropriate equationsfor the recovered people and the deceased people for COVID-19. During thedescent phase the number of active people over time must satisfy a conservationequation. This allows determining the time- evolution for the positive people.

10We define the number of positive people as the number of people tested positive for COVID-19, hence, by excluding the number of the deceased people and the number of people whorecovered.

14

N

From 01/03/2020

to 15/03/2020 From 16/03/2020 to 15/05/2020

Days20 40 60 80

10000

20000

30000

40000

50000

60000



I

Figure 9: Situation in Belgium on 15 May 2020—before, and 60 days after, theadoption of lockdown measures. The black dots correspond to real data. Theblue dotted line corresponds to the situation in Belgium before the adoptionof the lockdown measures. The blue and the red solid lines correspond to thetheoretical predictions for Belgium according to the solution of Eq. (14) and thelogistic model, respectively. Solution of Eq. (14) fits well all the experimentaldata from the initial days (i.e., from the 29th of February 2020), while thelogistic model applies only to the first data. The values of the parametersof Eq. (14) and the logistic function (10) are: τBE ' 5.3 days (µBE = 3.7days), KBEN ' 42626, and t0BE ' 53.4 days for Eq. (14), and τBE ' 5.3 days(µBE = 3.7 days), KBE = 111000, t0LBE = 39.5 days for the Logistic function,respectively. The zone I corresponds to the period before the adoption of thelockdown measures.

In the sequel, we denote with rt, dt, and nt the number of people released fromthe hospital at the time t, the total deaths, and the number of positive individualsat time t, respectively11

3.1 Number of the Recovered People

We start with the recovered people previously hospitalised. Let us suppose thata hospital has 50 patients in intensive therapy, corresponding to its maximumavailability capacity. If the hospital is unable to heal any patient, the growth rateof healed people is necessarily equal to zero. On the other hand, if the hospitalis able to heal a certain number of people, the places previously occupied by

11We draw the attention of the reader that in this manuscript Nt (capital letter) representsthe number of the total cases at time t, whereas nt (small letter) refers to the number of thepositive individuals at time t.

15

the sick people will free and other patients affected by COVID-19 will be ableto be hospitalised. In the latter case, the growth rate of the healed people willrise thanks to the fact that the hospital is able to heal more and more patients.This initial phase may be modelled by introducing into the dynamic equationthe term +γr(t), with r(t) indicating the number of the recovered people at thetime t, previously hospitalised

drtdt

= ζrt (25)

However, it is reasonable to suppose that ζ is constant for low values of rt,whereas, when rt takes more and more large values, ζ is proportional to Ir,with Ir denoting the number of the infected people entering in the hospital(and not the total number of the infected people, which is indicated with nt).Hence,

ζ = αr − βrIr > 0 (26)The sign minus in Eq. (26) is due to the fact that the recovered people willcontinue to grow linearly until when it reaches a maximum limit i.e. until whenthe hospital is no longer able to accept other sick people; this causes a reductionof people who recover. The competition between these two effects diminishesthe combined growth rate. Hence,

Ir <αrβr

(27)

where βr is linked to the hospital’s capacity to accept sick people12. To sum up,

drtdt

= αrr(t)− βrIr(t− θ)r(t) where Ir(t− θ) = r(t) +Dr(t+ θ1) (28)

r(t) is the number of the recovered people, previously hospitalised, at the timet who have been infected, in average, at the time t− θ and Dr(t+ θ1) denotesthe number of deceased people at the hospital at the time t+ θ1 who have beeninfected, in average, at the time t− θ (in general, θ1 6= 0).Clearly, the number of the recovered people, previously hospitalised, at the stepn (i.e. rn), is linked to the total number of the recovered people previouslyhospitalised at the step n (denoted by hn) by the relation

rn = hn − hn−1 or ht =n=t/∆t∑n=1

rn (with ∆t ' 1 day) (29)

where we have set h0 = 0.

3.1.1 Approximated O.D.E. for the Recovered People PreviouslyHospitalised

We assume that all the infected people entering in the hospitals will heal. So

Dr(t+ θ1) ≈ 0 hence Ir(t− θ) ' r(t) (30)12Indeed, βr is proportional to the inverse of the hospital capacity.

16

The final O.D.E. for the recovered people reads then

drtdt' αr

(1− 1

Krrt

)rt (31)

where Kr is the hospital’s capacity, which we assume to be a time-independentparameter. Notoce that, under this approximation, the equation for the numberof recovered people is in itself and independent of the equations for the othervariables (i.e. for nt and dt)

13.

3.1.2 O.D.E. for the Total Recovered People

At the first approximation, the O.D.E. for the total recovered people Rt (i.e. thetotal individuals having survived the disease) is trivially obtained by consideringthat the rate of Rt is approximatively proportional to the number of the infectedpeople nt at time t i.e.

14.dRtdt

= χnt (32)

However, it is useful to clarify the following. In Eqs (29), ht stands for thetotal number of the recovered people previously hospitalised whereas the variableRt in Eq. (32) is the total number of the recovered people (i.e. the number ofthe recovered people previously hospitalised, plus the number of the asymp-tomatic people, plus the infected people who have been recovered without beingpreviously hospitalised). The natural question is: ”how can we count Rt andcompare this variable with the real data ?”. The current statistics, produced bythe Ministries of Health of various Countries, concern the people released fromthe hospitals. Apart from Luxembourg (where the entire population has beensubject to the COVID-19-test), no other Countries are in a condition to providestatistics regarding the total people recovered by COVID-19. Hence, it is ouropinion that the equation for Rt, is not useful since it is practically impossibleto compare Rt with the experimental data.

3.2 Equation for the Deceased People

The rate of deceased people per unit time is modelled by the following dimen-sionless equation

d

dt̂dt = αdn(t−td) − βdn

2(t−td) (33)

The meaning of Eq. (33) is the following. Manifestly, the rate of deaths isproportional to the number of active people. However, individuals infected by

13Eq. (31) models a hospital’s ability to heal people and, by no means, it must be linkedto the number of people tested positive for COVID-19 or to the mortality rate caused by theSARS-CoV-2.

14Notice that Eq. (32) is the dynamic equation for the total recovered people adoptedin the Susceptible-Infectious-Recovered-Deceased-Model (SIRD-model) [29]. The comparisonbetween our model with the SIRD-model will be found soon in Ref. [30].

17

SARS-CoV-2 do not die instantly since the rate of deaths at time t is propor-tional to the people who were infected at an earlier time t− td (td > 0) with tddenoting the time-delay. We indicate with αd the, time-independent, constantproportional to the increase of the deaths dt. The second term, −βdn2(t−td),models the presence of the lockdown measures, having the effect of saturatingthe rate of infected people and, consequently, of deaths. Indeed, in the absenceof lockdown measures, we may approximatively write

d

dt̂dt = αrn(t−td) (34)

with αr denoting a positive constant. The purpose of the lockdown measuresis to decrease the number of infected people, and therefore deaths. We mayassume that the effect of the lockdown measures is proportional to n(t−td) suchas to dampen the linear growth of the mortality rate. In other terms, we get

αr → αr − βrn(t−td) (35)

which combined with Eq. (34) gives Eq. (33).

3.3 Equation for the Positive People

Of course, during the descent phase, the number of active people nt satisfies asimple law of conservation: If we are in the situation where there are no longernew cases of people tested positive for COVID-19 and if we assume that theactive people cannot leave their country of origin (or else, if they do, they willbe rejected by the host Country), then the number of infected people cannot butdecrease either because some people are deceased or because others have beenrecovered. In mathematical terms

nt = n0 − (ht − h0)− (dt − d0) = NMax − ht − dt (36)

with h0, d0 and n0 denoting the values of ht, dt and nt evaluated at the timet = tMax (see Eq. (19) i.e., the time that maximises the number of the totalcases), respectively. It should be noted that the conservation law (36) appliesonly when there are no longer new cases of people tested positive to COVID-1915. Here, by the descending phase we mean the phase where Eq. (36) applies.

15So, Eq. (36) does not apply necessarily as soon as the number nt (the number of peopletested positive for COVID-19) starts to decrease. Indeed, it may happen that nt decreasesbecause, for example, the number of new cases of people tested positive is less than the numberof the people who have recovered in the meantime. Conservation law (36) applies only fromthe moment where the number of new cases of people tested positive is strictly equal to zero.

18

To sum up, the equations describing the descending phase are

d

dt̂rt = αrrt

(1− rt

Kr

)with rt=tMax = r0 (37)

d


2(t−td) with dt=tMax = d0

nt = NMax − ht − dt with n∞ = 0

ht =

n=t/∆t∑n=1

rn where ∆t ' 1 day

with tMax given by Eq. (19). Notice that the first two equations of system (37)are valid also during the ascending-phase. Of course, in this case, the initialconditions are rt=0 = 0, dt=0 = 0 and nt=0 = 0. Hence, during the ascendingphase the evolution equations are

d

dt̂rt = αrrt

(1− rt

Kr

)with rt=0 = 0 (38)

d


2(t−td) with dt=0 = 0

nt = Nt − ht − dt with nt=0 = 0d

dt̂Nt = αNt

(1− Nt

KN

)−(αt̂2 − 1

t̂

)Nt

ht =

n=t/∆t∑n=1

rn where ∆t ' 1 day

According to our expectations, by solving numerically system (38), with goodapproximation, we get

dt ∝ n(t−td) (39)

Next, we find the numerical solution of systems (37)-(38) for Italy and Belgium.A similar analysis for Luxembourg is reported in Appendix.

3.4 Theoretical Predictions for the Descending Phase

In this subsection, we report the numerical solutions of Eqs (37)-(38) for Italyand Belgium. The solution for Luxembourg can be found in the Appendix.Fig. (10) and (11) concern the Italian situation. They show the numericalsolution of Eqs (37)-(38) for the number of recovered people and deaths, respec-tively. These theoretical predictions are plotted against the experimental datareported in the Table (1). According to the theoretical predictions, for Italy weget tITd = 12 days. Fig. (12), illustrates the descendant-phase for Italy.Figs (13) and (14) refer to the Belgian situation. The figures illustrate thenumerical solutions of Eqs (37)-(38) for the number of recovered people anddeaths, respectively. The theoretical predictions are plotted against the experi-mental data reported in the Table (3). According to the theoretical predictions,

19

Days0 20 40 60 80 100 120

20000

40000

60000

80000

100000

120000

140000

th

Figure 10: Italy situation. Theoret-ical predictions (blue line) againstthe experimental data (black cir-cles) for the recovered people.

Days

d t

0 20 40 60 80 100

5000

10000

15000

20000

25000

30000

Figure 11: Italy situation. Theoret-ical predictions (blue line) againstthe experimental data (black cir-cles) for the deceased people.

nIT

Days

Plateau of N The epidemic has decreased significantly

2 months

t

10 20 30 40 50 60

20000

40000

60000

80000

Current situation

Figure 12: The descending phase for Italy. According to the theoretical predic-tions, after two months the lockdown measures may heavily be lightened andwe can return to normal work. The estimated t ime-delay is tITd = 12 days -see Eq. (37)

for Belgium we get tBEd = 8.8 days. Fig. (15) shows the descendant-phase forBelgium.

20

0 20 40 60 80 100

5000

10000

15000

Days

th

Figure 13: Belgian situation.Theoretical predictions (blue line)against the experimental data(black circles) for the recoveredpeople.

0 20 40 60 80

2000

4000

6000

8000

10000

Days

d t

Figure 14: Belgian situation.Theoretical predictions (blue line)against the experimental data(black circles) for the deceasedpeople.

10 20 30 40

10000

20000

30000

nBE

Days


1 month

t

Current situation

Figure 15: The descending phase for Belgium. According to the theoreticalpredictions, after one month the lockdown measures may heavily be lightenedand we can return to normal work. The estimated t ime-delay is tBEd = 8.8days - see Eq. (37).

4 Conclusions

In this work we studied the spread of SARS-CoV-2 until when the strict lock-down measures have been adopted (i.e. until 16th May 200). The dynamics of

21

COVID-19 when the population is under less restrictive lockdown measures willbe subject of future studies. Through fitting techniques previously performed,and validated by the χ2-tests, we determined a general differential equation, hav-ing two-free parameters, which we propose for 13 respiratory infectious diseasescaused by viruses, including SARS-CoV-2. The solution of Eq. (14) providesthe number of the total case in time. Successively, we compared the theoreti-cal predictions, provided by the solution of Eq. (14) and by the logistic model(see Eq. (7)), with the real data for Italy and Belgium (for Luxembourg seeAppendix). We saw that the solution of Eq. (14) is in good agreement with theexperimental data since the beginning of the appearance of the COVID-19; thisis not the case for the logistic model which applies only to the few last days. Wefound the days where the maximum number infected people by COVID-19 willbe reached in Italy and Belgium by parametrising the solution of Eq. (14) withexperimental data: we get, tMaxIT ' 21 April 2020 and tMaxBE ' 2 May 2020for Italy and Belgium, respectively.We also noted, empirically, that the infection process caused by SARS-CoV-2may be divided into three qualitatively different periods; i.e., the exponentialperiod, the transient period and the bell-shaped period (or the post-transientperiod). The solution of Eq. (14) allows defining more precisely these threeperiods. Indeed, we may classify the above periods as follows

The exponential period for 0 ≤ t̂ ≤ t̂LM (40)The transient period for t̂LM < t̂ ≤ t̂flexThe bell-shaped period for t̂ > t̂flex

with t̂LM indicating the dimensionless time when the lockdown measures areapplied and t̂flex the dimensionless inflection point of the solution of Eq. (14),respectively. It is easily checked that, for large values of KN , the value of t̂flexsatisfies, approximatively, the equation

αt̂3flex − 2αt̂2flex + (α− 3)t̂flex + 2 ' 0 with t̂flex ≡tflext0

(41)

Hence, according to Eq. (41), the transient period ended on 31 March 2020 forItaly and on 7 April 2020 for Belgium, respectively. The second part of the workis devoted to modelling the descending phase, i.e. the decrease of the numberof people tested positive for COVID-19. Also in this case, we proposed a newset of dynamic differential equations that we solved numerically. The solutionof Eqs (37) (and Eq. (38)) provided valuable information such as the durationof the COVID-19 epidemic in a given Country and therefore when it will bepossible to return to normal life.

5 Acknowledgments

I am very grateful to Alberto Sonnino from Facebook Calibra and UniversityCollege London for comments on late manuscript, and to Ing. Alessandro Leonefrom the Italian Embassy in Belgium for his suggestions and encouragement.

22

6 Appendix: Comparison between the Theoret-ical Predictions of Eq. (14) and ExperimentalData for Luxembourg

Luxembourg, due to the particularly severe lockdown measures adopted by thegovernment, may be considered, with good approximation, as a closed system(628108 inhabitants, most of them concentrated in only one town). Indeed,right from the start, the city of Luxembourg was literally closed and citizenswere unable to enter and leave the city freely (people who had to enter in thecity for working reasons were obliged to undergo each time the test that, ofcourse, had to result negative).Italy, on the other hand may be considered, with a good approximation, asan open system (60317116 inhabitants dislocated in all the Country). In Italy,especially during the initial phase, the citizens of northern Italy moved freely tothe south of Italy, by train, by plans or by car. Only in a second time the Italiangovernment decided to introduce much more restrictive measures concerning themovement of citizens from one region to another.For the reason mentioned above, it is our opinion that it is very interesting toanalyse these two Countries, Luxembourg and Italy, which are so different witheach other. In this Appendix we report the comparison between the theoreticalpredictions of the COVID-19 model (14) and the real data for Luxembourgupdate to 15 May 2020 (see Fig. (16)). In the columns of table (5) we can findthe number of active people (currently infected by SARS-CoV-2), the numberof recovered people, and the number of deceased people, respectively. Theexperimental data have been found in the databases [31] and [32]. Luxembourgreached its peak on 12 April 2020.

6.1 The Descending Phase for Luxembourg

Figs (17) and (18) refer to the Luxembourg situation. The figures illustratethe numerical solutions of Eqs (37)-(38) for the number of recovered peopleand deaths, respectively. The theoretical predictions are plotted against theexperimental data, which can be found in the Table (5). Fig. (19) shows thedescending phase for Luxembourg.

References

[1] Il Sole 24 Ore, Coronavirus in Italia, i dati e la mappa (2020) - Website:https://lab24.ilsole24ore.com/coronavirus

[2] Sciensano, Epidemologie des maladies infectueuses (2020) - Website:https://epidemio.wiv-isp.be/ID/Pages/2019-nCoV.aspx

[3] Wikipedia, the free encyclopedia, 2020 coronavirus pandemic in Belgium,https://en.wikipedia.org/wiki/2020 coronavirus pandemic in Belgium

23

20 40 60 80

1000

2000

3000

4000



I

From 01/03/2020

to 16/03/2020 From 16/03/2020 to 15/05/2020

Days

Luxembourg situation after the adoption of

the lockdown measures

N

Figure 16: Situation in Luxembourg on 15 May 2020. The black dots correspondto real data. The red dotted line corresponds to the situation in Luxembourgbefore the adoption of the lockdown measures. The blue and the red solidlines correspond to the theoretical predictions for Luxembourg according to thesolution of Eq. (14) and the logistic model, respectively. Solution of Eq. (14)fits well all the experimental data from the initial days (i.e., from the 29th ofFebruary 2020), while the logistic model applies only to the first data. Thevalues of the parameters of Eq. (14) and the logistic function (10) are: τLUX '3.2 days (µLUX = 2.2 days), K

LUXN ' 6880, and t0LUX ' 40 days for Eq. (14),

and τLUX ' 3.2 days (µLUX = 2.2 days), KLUX = 3950, t0LLUX = 30 days forthe Logistic function, respectively. The zone I corresponds to the period beforethe adoption of the lockdown measures.

24

Days20 40 60 80

500

1000

1500

2000

2500

3000

3500

th

Figure 17: Luxembourg situa-tion. Theoretical predictions (blueline) against the experimental data(black circles) for the recovered peo-ple.

20 40 60 80

20

40

60

80

100

d t

Days

Figure 18: Luxembourg situa-tion. Theoretical predictions (blueline) against the experimental data(black circles) for the deceased peo-ple.

5 10 15 20 25

500

1000

1500

2000

2500

3000

nLUX

Days


~ 3 weeks

t

Current situation

Figure 19: The descending phase for Luxembourg. According to the theoreticalpredictions, after one month the lockdown measures may heavily be lightenedand we can return to normal work. The estimated t ime-delay is tLUXd = 15days - see Eq. (37).

[4] History and Epidemiology of Global Smallpox Eradication (Archived 2016-05-10 at the Wayback Machine), a module of the training course Smallpox:Disease, Prevention, and Intervention. The CDC and the World Health

25

Organisation, (2001). Slide 17. This gives sources as Modified from EpidRev, 15, 265-302 (1993), Am J Prev Med,20 (4S): 88-153 (2001), MMWR49 (SS-9); 27-38 (2000).

[5] Murray J.D., Mathematical Biology, Springer-Verlag Berlin, HeidelbergGmbH (1993).

[6] Lusuardi A., Modelli Matematici per lo Studio di Popolazioni Interagenti inun Ecosistema, Thesis dissertation in Mathematics, Università defli Studidi Bologna (2016).

[7] Anderson R.M. and May R.M., Infectious diseases of humans, Oxford Uni-versity Press (1991).

[8] https://populationeducation.org/what-doubling-time-and-how-it-calculated/

[9] Steinbrecher G. (University of Craiova - Romania), A short sample in R-program able to compute the estimated value of R0 for 17 infectious diseases.Private Communication (20th April 2020).

[10] Sonnino G, The COVID-19 - The Infectious Disease Caused by theLatest Discovered Coronavirus (SARS-CoV-2), European Commission,ARES(2020)1530456, (March 2020).

[11] N. Mathus and G. Shaw, Analysis of Indian data on Covid-19 using a dy-namical model, Department of Theoretical Physics, Tata Institute of Fun-damental Research, Mumbai (India) 2020 - Pre-print sent in April 2020.

[12] Yukalov, V. I.; Yukalova, E. P.; Sornette, D., Punctuated evolution dueto delayed carrying capacity. Physica D: Nonlinear Phenomena. 238 (17):1752?1767 (2009).

[13] Richards, F. J., A Flexible Growth Function for Empirical Use, Journal ofExperimental Botany, 10, 290?300 (1959).

[14] Ruffini R., COVID-19 statistics, ICRANet Newsletter. March 2020. Web-site: http://www.icranet.org/index.php?option=com content&task=view&id=1310#1

[15] Piovani P., Coronavirus, il pediatra: Questa epidemia durerà tre mesi, Ilmessaggero, 24 March 2020.

[16] Hilbe J/M., Practical Guide to Logistic Regression, CRC Press, Taylor &Francis Group, LLC (2015).

[17] On Health, Respiratory Illnesses: 13 Types of Lung Infections (Re-viewed by Charles Patrick Davis, MD, PhD on 3/20/2020). Website:https://www.onhealth.com/content/1/respiratory infections causes

26

http://www.icranet.org/index.php?option=com$_$content

[18] Manjarrez-Zavala E. et al., Pathogenesis of Viral Res-piratory Infection, (2013), DOI: 10.5772/54287. Website:https://www.intechopen.com/books/respiratory-disease-and-infection-a-new-insight/pathogenesis-of-viral-respiratory-infection

[19] Huang C. et al., Clinical features of patients infected with 2019 novel coron-avirus in Wuhan, China, The Lancet Journals, 395, Issue10223, P497-506,(Published: 24 January 2020). Website: https://doi.org/10.1016/S0140-6736(20)30183-5

[20] Garrett Nichols W. et al., Respiratory Viruses Other than InfluenzaVirus: Impact and Therapeutic Advances, American Society for Microbiol-ogy - Clinical Microbiology Review, DOI: 10.1128/CMR.00045-07, Website:https://cmr.asm.org/content/21/2/274

[21] Nickbakhsh S. et al., Virus-virus interactions impact the population dy-namics of influenza and the common cold, Proceedings of the NationalAcademy of Science of the United State of America (PNAS), PNAS, 116(52) 27142-27150 (2019); Edited by Burton H. Singer, University of Florida,Gainesville, FL, and approved (2019).

[22] Lauer S.A., The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Ap-plication, Annals of Internal Medicine (ACP), (10 March 2020).Website: https://annals.org/aim/fullarticle/2762808/incubation-period-coronavirus-disease-2019-covid-19-from-publicly-reported

[23] Chen T.-M., et al., A mathematical model for simulating the phase-basedtransmissibility of a novel coronavirus, BMC Part of Springer Nature, In-fectious Diseases of Poverty, 9, Article number: 24 (2020).

[24] Galvani A.P., Xiudong L. and Jewell N.P., Severe Acute Respiratory Sin-drome: Temporal Stability and Geographic Variation in Death Rates andDoubling Time, Emerging Infect ious Diseases, 991-994 (2003) - Website:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3020622/

[25] Sun X., Bao J., and Shao Y., Scientific Report, 6, 22498 (2016). doi:10.1038/srep22498

[26] https://en.wikipedia.org/wiki/Logistic function

[27] Nowak, M. A. and May, R. Virus Dynamics: mathematical principles ofimmunology and virology (2001). ISBN 9780198504177. Archived from theoriginal on 13 February 2015. Retrieved 13 February 2015.

[28] Gompertz, B., On the Nature of the Function Expressive of the Law of Hu-man Mortality, and on a New Mode of Determining the Value of Life Con-tingencies, Philosophical Transactions of the Royal Society, 115, 513?585(1825).

27

[29] Kermack W. O. and McKendrick, A. G., A Contribution to the Mathemat-ical Theory of Epidemics. Proceedings of the Royal Society A. 115 (772),700?721 (1927).

[30] Sonnino G., Nardone P, and Bassios V., A Mathematical Model Describingthe Spreading of the Virus SARS- CoV-2, Causing the COVID-19 in Pres-ence of the Lockdown Mea- sures, to be submitted to the review Chaos,Solitons & Fractals (September 2020).

[31] Site Coronavirus (COVID-19). Le gouvernementluxembourgeois. Ministère de la Sant. Website:https://msan.gouvernement.lu/fr/dossiers/2020/corona-virus.html

[32] Wikipedia: 2020 coronavirus pandemic in Luxembourg. Website:https://en.wikipedia.org/wiki/2020 coronavirus pandemic in Luxembourg

28

Date Active Recovered Deceased Total cases

25-Feb 322 1 10 33326-Feb 400 3 12 41527-Feb 650 45 18 71328-Feb 888 46 21 95529-Feb 1049 50 29 112801-Mar 1577 83 34 169402-Mar 1835 149 52 203603-Mar 2263 160 79 250204-Mar 2706 276 107 308905-Mar 3296 414 148 385806-Mar 3916 523 197 463607-Mar 5061 589 233 588308-Mar 7375 622 366 836309-Mar 9172 724 463 1035910-Mar 10149 1004 631 1178411-Mar 10590 1045 827 1246212-Mar 12839 1258 1016 1511313-Mar 14955 1439 1266 1766014-Mar 17750 1966 1441 2115715-Mar 20603 2335 1809 2474716-Mar 23073 2749 2158 2798017-Mar 26062 2941 2503 3150618-Mar 28710 4025 2978 3571319-Mar 33190 4440 3405 4103520-Mar 37860 5129 4032 4702121-Mar 42681 6072 4825 5357822-Mar 46638 7024 5475 5913723-Mar 50418 7432 6077 6392724-Mar 54030 8326 6820 6917625-Mar 57511 9362 7503 7437626-Mar 62013 10361 8165 8053927-Mar 66414 10950 9134 8649828-Mar 70065 12384 10023 9247229-Mar 73880 13030 10779 9768930-Mar 75528 14620 11591 10173931-Mar 77635 15729 12428 10579201-Apr 80572 16847 13155 11057402-Apr 83049 18278 13915 11524203-Apr 85388 19758 14681 11982704-Apr 88274 20996 15362 12463205-Apr 91246 21815 15887 12894806-Apr 93187 22837 16523 13254707-Apr 94067 24392 17127 13558608-Apr 95362 26491 17669 13942209-Apr 96877 28470 18279 14362610-Apr 98273 30455 18849 14757711-Apr 102253 34211 19899 15227112-Apr 100269 32534 19468 15636313-Apr 103616 35435 20465 15951614-Apr 104291 37130 21067 16248815-Apr 105418 38092 21645 165155

Table 1: Situation in Italy on 15 May 2020. Columns report the number of active people(currently infected by SARS-CoV-2), the number of recovered people, and the number ofdeceased people.

29


16-Apr 106607 40164 22170 16894117-Apr 106962 42727 22745 17243418-Apr 107771 44927 23227 17592519-Apr 108257 47055 23660 17897220-Apr 108237 48877 24114 18122821-Apr 107709 51600 24648 18395722-Apr 107699 54543 25085 18732723-Apr 106848 57576 25549 18997324-Apr 106527 60498 25969 19299425-Apr 105847 63120 26348 19535126-Apr 106103 64928 26644 19767527-Apr 105813 66624 26977 19941428-Apr 105205 68941 27359 20150529-Apr 104657 71252 27682 20359130-Apr 101551 75945 27967 20546301-May 100946 78249 28236 20742802-May 100704 79914 28710 20932803-May 100179 81654 28884 21071704-May 99980 82879 29079 21193805-May 98467 85231 29315 21301306-May 91528 93245 29684 21445707-May 89624 96276 29958 21585808-May 87961 99023 30201 21718509-May 84842 103031 30395 21826810-May 83324 105186 30560 21907011-May 82488 106587 30739 21981412-May 81266 109039 30911 22121613-May 78457 112541 31106 22210414-May 76440 115288 31368 22309615-May 72070 120205 31610 22388516-May 70187 122810 31763 22476017-May 68351 125176 31908 22543518-May 66553 127326 32007 22588619-May 65129 129401 32169 22669920-May 62752 132282 32330 22736421-May 60960 134560 32486 22800622-May 59322 136720 32616 22985823-May 57752 138840 32735 22932724-May 56594 140479 32785 22985825-May 55300 141981 32877 23015826-May 52942 144981 32955 23055527-May 50966 147101 33072 23113928-May 47986 150604 33142 231732

Table 2: Situation in Italy on 15 May 2020. Columns report the number of active people(currently infected by SARS-CoV-2), the number of recovered people, and the number ofdeceased people.

30


29-Feb 1 0 0 101-Mar 1 0 0 102-Mar 6 0 0 603-Mar 13 0 0 1304-Mar 23 0 0 2305-Mar 50 0 0 5006-Mar 109 0 0 10907-Mar 169 0 0 16908-Mar 200 0 0 20009-Mar 239 0 0 23910-Mar 267 0 0 26711-Mar 311 0 3 31412-Mar 396 0 3 39913-Mar 556 0 3 59914-Mar 686 0 4 68915-Mar 881 1 4 88616-Mar 1052 1 5 105817-Mar 1218 20 5 124318-Mar 1441 31 14 148619-Mar 1619 155 21 179520-Mar 2016 204 37 225721-Mar 2485 263 67 281522-Mar 2986 340 75 340123-Mar 3305 350 88 374324-Mar 3737 410 122 426925-Mar 4234 547 178 493726-Mar 5340 675 220 623527-Mar 6398 858 289 728428-Mar 7718 1063 353 913429-Mar 9046 1359 431 1083630-Mar 9859 1527 513 1189931-Mar 10374 1696 705 1277501-Apr 11004 2132 828 1396402-Apr 11842 2495 1011 1534803-Apr 12755 2872 1143 1677004-Apr 13901 3247 1283 1843105-Apr 14493 3751 1447 1969106-Apr 15196 3986 1632 2081407-Apr 16002 4157 2035 2219408-Apr 16482 4681 2240 2340309-Apr 17296 5164 2523 2498310-Apr 18080 5568 3019 2666711-Apr 18686 5986 3346 2801812-Apr 19584 6463 3600 2964713-Apr 19979 6707 3903 3058914-Apr 20094 6868 4157 3111915-Apr 22025 7107 4440 33573

Table 3: Situation in Belgium on 15 May 2020. Columns report the number of activepeople (currently infected by SARS-CoV-2), the number of recovered people, and the numberof deceased people.

31



Table 4: Situation in Belgium on 15 May 2020. Columns report the number of activepeople (currently infected by SARS-CoV-2), the number of recovered people, and the numberof deceased people.

32


29-Feb 1 0 0 101-Mar 1 0 0 102-Mar 1 0 0 103-Mar 1 0 0 104-Mar 1 0 0 105-Mar 2 0 0 206-Mar 3 0 0 407-Mar 4 0 0 408-Mar 5 0 0 509-Mar 5 0 0 510-Mar 7 0 0 711-Mar 7 0 0 712-Mar 26 0 0 2613-Mar 33 0 1 3414-Mar 50 0 1 5115-Mar 76 0 1 7716-Mar 80 0 1 8117-Mar 139 0 1 14018-Mar 201 0 2 20319-Mar 325 6 4 33520-Mar 474 6 4 48421-Mar 656 6 8 67022-Mar 780 10 8 79823-Mar? 857 10 8 87524-Mar 1081 10 8 109925-Mar 1315 10 8 133326-Mar 1434 10 9 145327-Mar?? 1550 40 15 160528-Mar 1773 40 18 183129-Mar 1889 40 21 195030-Mar 1926 40 22 198831-Mar 2115 40 23 217801-Apr?? 2250 40 29 231902-Apr 2417 40 30 248703-Apr 2081 500 31 261204-Apr 2198 500 31 272905-Apr 2268 500 36 280406-Apr 2302 500 41 284307-Apr 2426 500 44 297008-Apr 2488 500 46 303409-Apr 2563 500 52 311510-Apr 2669 500 54 322311-Apr 2708 500 62 327012-Apr 2715 500 66 328113-Apr 2725 500 67 329214-Apr 2740 500 67 330715-Apr 2778 520 69 3373

? Attention on the 23rd the figures instead of being given at 9 am are given at 5 pm.?? Including 1 evacuated from France.

Table 5: Situation in Luxembourg on 15 May 2020. Columns provide the number of activepeople (currently infected by SARS-CoV-2), the number of recovered people, and the numberof deceased people.

33



Table 6: Situation in Luxembourg on 15 May 2020. Columns provide the number of activepeople (currently infected by SARS-CoV-2), the number of recovered people, and the numberof deceased people.

34

1 Introduction1.1 Comparison with the Real Data for COVID-19 before the Lockdown Measures

2 Modelling the COVID-19 - Virus' growth2.1 General background2.2 Determination of the Carrying Capacity and the Lockdown Coefficient for the COVID-192.3 Comparison between the Theoretical Predictions and Experimental Data

3 Modelling the COVID-19 - The Descending Phase3.1 Number of the Recovered People3.1.1 Approximated O.D.E. for the Recovered People Previously Hospitalised3.1.2 O.D.E. for the Total Recovered People

3.2 Equation for the Deceased People3.3 Equation for the Positive People3.4 Theoretical Predictions for the Descending Phase

4 Conclusions5 Acknowledgments6 Appendix: Comparison between the Theoretical Predictions of Eq. (14) and Experimental Data for Luxembourg6.1 The Descending Phase for Luxembourg

returning to normal life - arxiv · giorgio sonnino and pasquale nardone facult e des sciences...

Documents