Next Article in Journal
Studying the Effect of Noise on Pricing and Marketing Decisions of New Products under Co-op Advertising Strategy in Supply Chains: Game Theoretical Approaches
Next Article in Special Issue
A Novel Technique to Control the Accuracy of a Nonlinear Fractional Order Model of COVID-19: Application of the CESTAC Method and the CADNA Library
Previous Article in Journal
Overview and Choice of Artificial Intelligence Approaches for Night-Time Adaptive Optics Reconstruction
Previous Article in Special Issue
Age-Specific Mathematical Model for Tuberculosis Transmission Dynamics in South Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Stochastic Kinetic Type Reactions Model for COVID-19

1
Faculté de Sciences, Campus de la Plaine CP 231, Université Libre de Bruxelles (ULB), Bvd du Triomphe, 1050 Brussels, Belgium
2
Institut de Physique de Nice, CNRS, Parc Valrose, Université Côte d’Azur, 06108 Nice, France
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(11), 1221; https://doi.org/10.3390/math9111221
Submission received: 10 April 2021 / Revised: 9 May 2021 / Accepted: 14 May 2021 / Published: 27 May 2021
(This article belongs to the Special Issue Mathematical Modeling and Analysis in Biology and Medicine)

Abstract

:
We propose two stochastic models for the Coronavirus pandemic. The statistical properties of the models, in particular the correlation functions and the probability density functions, were duly computed. Our models take into account the adoption of lockdown measures as well as the crucial role of hospitals and health care institutes. To accomplish this work we adopt a kinetic-type reaction approach where the modelling of the lockdown measures is obtained by introducing a new mathematical basis and the intensity of the stochastic noise is derived by statistical mechanics. We analysed two scenarios: the stochastic S I S -model (Susceptible ⇒ Infectious ⇒ Susceptible) and the stochastic S I S -model integrated with the action of the hospitals; both models take into account the lockdown measures. We show that, for the case of the stochastic SIS-model, once the lockdown measures are removed, the Coronavirus infection will start growing again. However, the combined contributions of lockdown measures with the action of hospitals and health institutes is able to contain and even to dampen the spread of the SARS-CoV-2 epidemic. This result may be used during a period of time when the massive distribution of vaccines in a given population is not yet feasible. We analysed data for USA and France. In the case of USA, we analysed the following situations: USA is subjected to the first wave of infection by Coronavirus and USA is in the second wave of SARS-CoV-2 infection. The agreement between theoretical predictions and real data confirms the validity of our approach.

1. Introduction

The current number of infectious people with SARS-CoV-2 brought us back to the lockdown phase, but tracking and swabs give a clearer vision of the situation. By reading the figures on new cases of positivity, it seems that we are returning back to mid-April 2020. Hospitalisations are also growing, including new hospitalisations in intensive care units. In this context, the situation appears very different compared to the first months of the epidemic, when the prolonged lockdown was necessary, which allowed reduction in the spread of the Coronavirus in a significant way. The percentage ratio between new cases and swabs now stands at around 5% [1]. It is therefore from 30 percent of the “darkest” periods of the pandemic in certain European Countries. However, these data tend to increase further preventive strategies, focusing on tracing and beyond, taking into account that in some countries the “weight” of the pandemic is starting to be felt again on hospitals. So, we are in the middle of the second wave, but the virus circulates among humans and therefore the lockdown period has limited its spread, so much so that we reached a sort of “stagnation” in the months of June and July with a subsequent resumption in August. However, the virus did not change and the clinical manifestations that it induces have not changed [2]. In this second period of COVID-19, we may trust statistical data much more thanks to the increased number of swabs done and the tracing and monitoring of people who have contracted the SARS-CoV-2 infection. Of course, there are still some cases, especially for asymptomatic subjects who unfortunately can transmit the infection: about 20% of cases of contagion can in fact be attributed to subjects who do not show any health problems. This situation also explains why the percentage of positive swabs in the first weeks of the lockdown was around 30 percent, while today it is much lower. Anyway, we are still far from the figures of several thousand of last March-April, but the trend is still on the rise. For this reason, it is essential also to hold the infections through “ad hoc” measures (lockdown, quarantine, isolation). The carelessness during the summer holidays has created a clear epidemiological link, as the number of cases started to rise again in August and now we have to face the situation, also because there are signs from nearby France warning of this problem.
• USA situation. After sustained declines in the number of COVID-19 cases over recent months, restrictions are starting to ease across the United States. The numbers of new cases are either falling or stable at low numbers in some states, but they are surging in many others. Overall, the USA is experiencing a sharp increase in the number of new cases per day and, by late June, had surpassed the peak rate of the spread in early April. In the USA, cases are rising quickly right now. With more than 11 million confirmed cases, the USA has the highest number of infections in the world and the spread of the Coronavirus shows no sign of slowing down. During the spring wave, swabs were mostly limited to confirming cases in people who were already in hospitals, meaning the true scale of that outbreak was not fully captured. However the latest data compiled by the COVID Tracking Project shows the current surge is not just down to increased testing; the number of tests carried out in USA was up by 12.5% week on week, while the number of cases increased by more than 40% [3]. One likely cause is the change of season and colder weather driving people indoors to socialise, where the risk of spread is heightened due to less social distancing and poor ventilation. The strain on hospitals is growing. Because of the change in the level of testing, a better way to compare waves is to look at the number of people being admitted to hospitals because of COVID-19. These data show that roughly the same number of people across the USA were in hospital during the first and second waves of the outbreak. However, there are already more people in hospitals during the current wave—more than 70,000 at the moment. New deaths per day have not begun to climb, but some hospitals’ intensive care units have recently reached full capacity [4]. This wave is hitting every USA region. We can say that the USA as a whole is not in a second wave because the first wave never really stopped. The virus is simply into new populations or is resurgent in places that let down their guard too soon. This time around, every region is seeing a spike in new cases [3].
• France situation. Coronavirus cases in France rose dramatically, but the number of daily deaths declined. The total number of infections stands at around of 1,829,659, with the death toll standing at 42,207, according to the Ministry of Health. The overall downward trend is a positive sign after Coronavirus numbers increased drastically by reaching a peak of 86,852 cases. Some 99 departments throughout France remain in a vulnerable position with 3042 clusters of infection considered serious [5]. France now has the fourth highest number of infections in the world behind the USA, India and Brazil. Despite the French government’s efforts to slow down the spread of the virus with targeted measures such as lockdown and isolation, the rate of infections has continued to rise across the country. Talks regarding the occurrence of a second wave, and how to avoid it, has been around since the end of the spring lockdown, but the current situation seems to suggest it was unavoidable. One month after the lockdown measures were applied in major cities starting in Paris, and three weeks after extending it to almost the whole country, France can see the light at the end of the dark tunnel of the second wave of Coronavirus. The French Health Minister announced that he believes that the peak of the pandemic has passed, and control has been regained. Indeed, looking at the data, after a peak of 86,000 cases as recorded on November 7, the country has moved on to figures of around 30,000 daily infections with a minimum level of 23,000 as recorded on November 13. Figures are returning in line with those of a month ago when the lockdown began [6]. However, the country will remain in a second nationwide lockdown whicch will continue until Dec. 1 unless the number of infections and deaths decrease in a significant way.
• Modelling the spread of SARS-CoV-2 infection. Until the vaccine arrives, the tried-and-true public health measures taken on the last months (social distancing, universal mask wearing, frequent hand-washing and avoiding crowded indoor spaces) are the ways to stop the first wave and thwart a second one. In addition, when there are surges such as those happening now in the world, further reopening plans need to be put on hold. Meanwhile, in order to understand the situation, theoretical studies, mainly based on mathematical models with the use of up-to-date data, must be carried on in order to understand the of this pandemic in the future as well as to provide predictions. In this study, we propose a simple stochastic model with the objective to show how public health interventions can influence the outcome of the epidemic, by limiting and even dampening the of the SARS-CoV-2 epidemic, while waiting for the production and delivery of an effective vaccine. The analysis is performed by assigning the population to different compartments with labels, such as S, I, R and D (Susceptible, Infectious, Recovered, and Deaths). People can move between the compartments. The order of labels shows the flow patterns between compartments; for example S I S means susceptible, infectious and, therefore, Susceptible again. The origin of these models is at the beginning of the twentieth century [7]. The goal of our model is to predict how a disease spreads, or the dynamics of the number of infectious people, as well as the duration of the COVID-19 pandemic. The dynamics of this kind of models is governed by ordinary differential equations (O.D.E.s). Generally, these O.D.E.s are deterministic (see, for example, [8,9,10,11,12]).
In this work the dynamics is governed by stochastic differential equations, which describe reality in a more realistic way. We would like to stress that testing, e.g., swabs, may provoke fluctuations in experimental data, but these kinds of fluctuations are not quoted in the current work. The stochasticity dealt with in this manuscript concerns the dynamics of intrinsic fluctuations which arise spontaneously in a macroscopic system, and which satisfy the laws of statistical mechanics [13]. The main objective of our series of studies is to obtain the correct space-time stochastic differential equations which could describe realistic situations of spread of SARS-CoV2 pandemic in large countries. This goal can only be achieved if we are able to
(1) Model the distribution of hospitals and health institutes in a country;
(2) Model the distribution of the poles of attraction of susceptible people (e.g., shopping centers, workplaces, etc.);
(3) Identify a mechanism able to establish when a pole of attraction is “saturated” with infected people, by proposing alternative poles of attraction;
(4) Modelling the Lockdown and Quarantine measures adopted by the Government of a particular country;
(5) Describe the dynamics of intrinsic (i.e., spontaneous) fluctuations to which a macroscopic system is subjected by methods of statistical mechanics.
At first glance, such a working program would appear to be too ambitious and, to our knowledge, the state-of-the-art of current alternative techniques are unable to resolve the issues listed above. The Kinetic-Type Reaction (KTR) approach proposed by us is very promising and allows to achieve this goal in a relatively simple way. Indeed the kinetic-type reaction approach
(a) Models each actor by a dedicated chemical species that can only be created or destroyed as the result of one, or several, elementary steps;
(b) Allows to determine the dynamics of the system starting from elementary steps;
(c) Thanks to its flexibility, allows to analyse complex situations when several variables are involved, such as R, Q, R h , I h etc.
We anticipate that, by adding extra space-time dependent contributions, inspired by a biological analogy, to the deterministic dynamic equations obtained by applying the KTR model, we have been able to determine the correct space-time stochastic partial differential equations which are suitable to describe the spread of SARS-CoV2 infection not only in single regions, but also in space extended countries. A work and related interactive simulations, where predetermined distributions of poles of attraction are set up, is in progress. In our view, this is the novelty of our approach with respect to models currently proposed in the literature. To the best of the authors’ knowledge, our approach has never been proposed before.
The paper is organised as follows. In Section 2, we study the deterministic S I S -model in presence of the lockdown measures. We shall see that as soon as the lockdown measures are stopped, the of Coronavirus begins to grow back vigorously once a certain threshold is exceeded. In Section 2.1, we ask ourselves the following question: Let us suppose a situation in which a government needs a number x of weeks (e.g., 6 weeks) before being able to deliver an effective vaccine. What is the minimum expected value of infectious people that should be seen, after the adoption of severe lockdown measures imposed on the population, in order to be sure that after x-weeks (e.g., after 6 weeks) the total number of infectious individuals by Coronavirus remains below a threshold pre-set in advance, and acceptable from a medical and political point of view? Section 2.1 answers the above-mentioned question for the S I S -model subjected to lockdown measures. In Section 3, we introduce, and study, the Stochastic S I S -model subjected to lockdown measures (the ( S I S ) L -model). We compute the relevant correlation functions as well as the probability density function. Hospitals and health institutes play a crucial role in hindering the of Coronavirus. In Section 4, we propose a model that accounts for people who are only traced back to hospitalised infectious individuals. This model is referred to as the ( S I S I h r h d h ) L -model. We shall see that beyond a certain threshold of hospital capacities, the coronavirus is not only kept under control, but its tends to diminish in time. Therefore, the combined effect of lockdown measures with ef = ficient hospitals and health measures are able to contain and even dampen the spread of the SARS-CoV-2 epidemic, while waiting for the delivery of a limited number of vaccines in a given population. In Section 5, we study the simplified ( S I S I h ) L stochastic model by computing the relevant correlation functions.
This approach is in line with the procedure normally adopted by the scientific community. When a process is complex, it is customary to start by analysing the simplest model, the SI(R)S model in this case, subjected to lockdown measures and perturbed by white noise. Of course, the role of the lockdown measures is not to shorten the infectious period of an infected individual, but rather to intervene on susceptible people. The noise takes into account the intrinsic fluctuations of the infected individuals, which are governed by the laws of statistical mechanics. These tasks have been accomplished by adopting a kinetic-type reactions approach, where the compartments S, I, and R are seen as “chemically active molecules” (see Section 2), and by showing two theoretical results which we report for the first time in the literature, namely:
(1) The lockdown measures (also applicable to the quarantine and isolation measures) are modelled by introducing a new (complete) mathematical base of functions (see Equation (12));
(2) The intensity of the white noise is derived by the laws of statistical mechanics (see Equation (A11)).
Comparison between theoretical predictions and real data for USA and France can be found in Section 6. Concluding remarks are reported in Section 7. The derivation of the exact solution of Richards equation in presence of the restriction measures as well as an estimation of the intensity of the noise by Statistical Mechanics methods can be found in the Appendices Appendix A and Appendix B.

2. The Deterministic Model in Presence of the Lockdown Measures-the (SIS)L-Model

The current work starts from the following hypothesis commonly supported by the most accredited virologists: the SARS-CoV-2 behaves like other viruses which cause respiratory diseases (see, for instance, ref. [14]). The common cold and influenza, do not give any long-lasting immunity. Such infections do not give immunity upon recovery from infection, and individuals become susceptible again. Hence, according to the above-cited hypothesis we propose the following simplest compartmental model:
S + I μ 2 I
In our model the SARS-CoV-2 infection does not leave any immunity, thus individuals, after recovery, return back into the S compartment. This added detail can be seen by including an R class in the middle of the model
S + I μ 2 I I γ R γ 1 S
In this work, Equation (2) has to be interpreted as follows. The entire process is described by adopting a kinetic-type reactions approach where the lockdown measures are modelled by some kind of inhibitor reactions where susceptible individuals can be trapped into inactive states. In addition, the substrate is associated to the infected people and the product to the recovered people, respectively. From Equation (2), we get the following O.D.E.s for S, R, and I:
d S d t = σ S N T o t . I + γ 1 R d I d t = σ S N T o t . I γ I d R d t = γ I γ 1 R
with N T o t . denoting the total population and σ μ N T o t . . By assuming that the dynamics of R is much faster that those of S and I, we may set d R / d t 0 and system (3) reduces to
d S d t σ S N T o t . I + γ I d I d t = σ S N T o t . I γ I
which corresponds to the model
S + I μ 2 I I γ S
In the relevant literature, model (5) is referred to as the S I S -model (see, for example, [15]). Notice that condition d R / d t 0 was adopted only for simplicity. This allows to obtain relevant results by performing calculations analytically, without having to solve numerically system (3). From Equation (4) we get the conservation law
d S d t + d I d t = 0 or S + I = N T o t .
Hence, the dynamics of infectious people is governed by the logistic model
d I d t = ( σ γ ) I 1 σ N T o t . ( σ γ ) I
or
d I d t = α ˜ I 1 I K with α ˜ σ 1 γ σ ; K N T o t . 1 γ σ
where α ˜ and K denote the linear growing rate of the COVID-19 and the carrying capacity, respectively.
Lockdown measures are mainly based on an isolation of susceptible individuals, eventually with removal of infected people by hospitalisation. In our model, the effect of lockdown measures are taken into account by introducing in the S I S -model the lockdown-induced decrease rate c ( t )
S + I μ 2 I I γ + c ( t ) S
where
c ( t L ) = 0 for t = t L c ( t ) > 0 for t > t L
in which t L denotes the time when the lockdown measures are applied. The meaning of kinetic reactions (9) is the following. Clearly, lockdown measures act on the susceptible individuals but they are unable to affect the chemical capacity (i.e., to decrease the infection capacity) of the Coronavirus. Hence, the value of the kinetic reaction μ must remain constant and the final effect of these lockdown measures is to increase the number of susceptible people “at the expense” of infected people (since S + I = c o n s t . ). According to the model, this can be done only by increasing the value of the kinetic constant γ . Note that if the effect of the block were to decrease the value of the kinetic constant μ , we would get the non-sense result that the carrying capacity of infected persons would increase over time. Finally, the corresponding deterministic differential equations for the COVID-19 model in presence of the lockdown measures reads then:
d S d t = σ S N T o t . I + γ I + c ( t ) I d I d t = σ S N T o t . I γ I c ( t ) I
Equation (9) may be referred to as the ( S I S ) L -model where L stands for L o c k d o w n . The general expression for the lockdown-contribution, c ( t ) , may be obtained as follow.
(1) At t = t L , c ( t L ) = 0 . This corresponds to the requirement that the lockdown measures start at time t = t L ;
(2) c ( t ) must always be a positive function for all t > t L and for all kind of scenarios;
(3) for all t > t L and for all kind of scenarios, c ( t ) must be a function able to bend, downwardly, the trend of the curve of infected people.
After that, we may proceed as usual. We look for a class of functions satisfying the above three requirements under the form of polynomials (or fractions of polynomials) by bearing in mind that the degrees of the polynomials must be chosen in such a way that the curve of infected people can be bent downwards. This requirement must be satisfied by each term of the series.
(4) Ultimately, we have to prove that the obtained series constitutes a complete base of functions. Here, we shall not discuss on the completeness of the basis functions { t ^ i + 1 t ^ j } i , j = 0 , 1 , . This will be subject of a future work.
Finally, the general expression for the lockdown measures, c ( t ) , may be cast into the form
c ( t ^ ) = i , j = 0 β i j t ^ i + 1 t ^ j for t ^ 1
with t ^ denoting the normalised time t ^ t / t L and β i j are real numbers. Note that, the above conditions (2) and (3) impose strict restrictions on the coefficients β i j . Now, we look for a lockdown expression which takes into account only the most relevant terms. The first two terms of the series (12) are
c ( t ^ ) β 00 ( t ^ 1 ) + β 01 t ^ 2 1 t ^ with t ^ 1
In References [16,17,18] we can find a huge amount of fittings of experimental data carried out for different countries and several SARS-diseases. The results of the fittings have shown that the most relevant contribution is expressed by the term ( t 2 t L 2 / t ) while the term ( t t L ) provides an irrelevant contribution in comparison with the previous one. Of course, there are situations where the term ( t t L ) turns out to be significant, e.g., in chemotherapy that starts with a log-kill effect. In this case, c(t) is the therapy-induced death rate, modelled by the term ( t t C ) since it is imposed at the beginning, and finely calibrated, from the outside. Finally, in our case we get
c ( t ^ ) β t ^ 2 1 t ^ with β > 0 and t ^ 1
with β β 01 = c o n s t . denoting the intensity of the lockdown measures. Equation (14) is the expression for the lockdown measures which will be considered in this work. Notice that Equation (14) generalises the lockdown term introduced in Reference [16]. By taking into account Equation (6), the deterministic differential equation for the infectious people, in presence of lockdown measures, reads
d I d t ^ = α I 1 I K β t ^ 2 1 t ^ I
with α α ˜ t L . As shown in the Appendix A, the exact solution of Equation (15) reads
I = A t ^ β exp ( ( t ^ α / β ) 2 / σ ) 1 + A α / K 1 t ^ x β exp ( ( x α / β ) 2 / σ ) d x
where σ 2 / β . Equation (16) can be rewritten into the form
I = I 0 t ^ β exp ( ( 1 α / β ) 2 / σ ) exp ( ( t ^ α / β ) 2 / σ ) 1 + ( I 0 α / K ) exp ( ( 1 α / β ) 2 / σ ) 1 t ^ x β exp ( ( x α / β ) 2 / σ ) d x
with I 0 denoting the value of the total cases at the time when the lockdown measures are applied, i.e., I t ^ = 1 = I 0 . Notice that for large values of the carrying capacity, Equation (17) has the limiting behaviour
I I 0 t ^ β exp ( ( 1 α / β ) 2 / σ ) exp ( ( t ^ α / β ) 2 / σ )
Figure 1 shows two solutions of model (15) for Italy, during the first wave of infection by SARS-CoV-2. The red dotted line refers to the solution without the lockdown measures and the dark dotted line to the solution when the lockdown measures were applied.

2.1. The Second Wave of COVID-19 Pandemic in the ( S I S ) L -Model

The basic reproduction number, R 0 , is defined as the expected number of secondary cases produced by a single (typical) infection in a completely susceptible population. An epidemic occurs if the number of infected individuals increases, i.e., d I / d t > 0
σ S N T o t . I γ I > 0 or σ γ S N T o t . > 1
In this case, R 0 = σ / γ . When lockdown measures are interrupted, the dynamics of the process is governed again by Equation (8). For Italy, in the first wave we have R 0 > 1 , so K > 0 . From Equation (8) we get
σ γ > 1 hence α > 0 and lim t ^ I = K
The above relations mean, from a biological perspective, that the disease will remain permanently endemic into the population. Notice that the inflection point of the logistic curve (8) occurs at time t = t f with
t ^ f = 1 α log K I L 1
with I L denoting the number of infectious people when the lockdown measures start to be applied. It is easily checked that the inflection point, P f , has coordinates P f = ( t ^ f , I f = K / 2 ) . In our forthcoming reasoning, we shall see that the inflection point plays an important role in our analysis [19]. To summarise, after removing the lockdown measurements, say at time t ^ 1 , the number of infectious people will start to rise again from the value I 1 , with I 1 I t = t ^ 1 , following again a growth rate given by the logistic curve (8). This new growth of outbreak is referred to as the second wave of COVID-19. Now, our goal is to answer the following question: “Given a fixed value of infectious people, say I ˜ , determine the minimum value that I 1 should have so that the number of infectious people is less than I ^ , during an entire period from t ^ 1 to t ^ f ”. We note that within the interval t ^ 1 t ^ t ^ f the number of infectious people grows quite slowly. Let us decide to interrupt the lockdown measures at time t ^ 2 . Hence, during the time interval Δ t ^ = t ^ 2 t ^ 1 , the population will no longer be subjected to lockdown measures. To answer the above-mentioned question, firstly we have to determine the range-values of I 1 such that I 1 < I f = K / 2 for all t ^ Δ t ^ , with Δ t ^ < t ^ f t ^ 1 . It is easily checked that I 1 should satisfy the condition
0 I 1 < K 1 + exp ( α Δ t ^ )
Now, we have to require that after a lapse of time Δ t ^ , the number of infectious individuals must be lower than a fixed a priori value I ˜ , with I ˜ < K / 2 . The result is easily obtained by writing the initial condition as
I 1 = K 1 + exp ( α Δ t ^ ) K x with 0 < x 1 1 + exp ( α Δ t ^ )
and by determining the range of values of x such that the inequality
K 1 + ( K / I 1 1 ) exp ( α Δ t ^ ) I ˜ < K 2
is satisfied. After a little algebra, we get
0 I 1 I ˜ 2 + exp α Δ t ^ + exp ( α Δ t ^ ) ( 1 + exp ( α Δ t ^ ) ) ( I ˜ / K ( exp ( α Δ t ^ ) ) + exp ( α Δ t ^ ) + 1 + exp ( α Δ t ^ ) )
To summarise, if the lockdown measures are interrupted at time t ^ 1 , during the successive period of time Δ t ^ = t ^ 2 t ^ 1 the number of infectious remains inferior to a preset number I ˜ if the number of the infectious at the end of the lockdown measures is inferior to I 1 , with I 1 given by the r.h.s. of Equation (25). A schematic representation of the second wave of SARS-CoV-2 infection is shown in Figure 2.
As an example, let us consider the case of the Italian situation. Figure 3 shows I 1 versus I ˜ and Δ t ^ such that 50,000 I ˜ K I T .
As shown in Figure 4, if we set the value I ˜ = 50,000, the maximum value of I 1 is I 1 9500 for t t L , and I 1 700 for t 2 t L , respectively.

3. The Stochastic (SIS)L-Model

The Stochastic version of Equation (15) reads
d I d t ^ = α I 1 I K β t ^ 2 1 t ^ I + ξ ( t ^ )
with ξ ( t ^ ) denoting a white noise
< ξ ( t ^ ) > = 0 < ξ ( t ^ ) ξ ( t ^ ) > = η δ ( t ^ t ^ )
and η is the intensity of the noise. δ denotes the Dirac delta function (distribution). Let us now consider the system at the reference state I R S , which is the solution of the deterministic equation
d I R S d t ^ = α I R S 1 I R S K β t ^ 2 1 t ^ I R S
subjected to a perturbation of small amplitude δ I , i.e.,
I ( t ^ ) = I R S ( t ^ ) + δ I ( t ^ )
Our goal is to compute the relevant statistical correlation functions of this processes, i.e., < δ I ( t ^ ) ξ ( t ^ ) > and < δ I ( t ^ ) δ I ( t ^ ) > . Our task can easily be accomplished by recalling the following theorem (see, for example, [20,21,22,23]):
For a system of Langevin equations describing the temporal evolution of the processes ( x 1 , , x n ) :
d x i d t = f i ( x 1 , , x n ) + ξ i ( t ) i = 1 , , n < ξ i ( t ) > = 0 ; < ξ i ( t ) ξ j ( t ) > = η i j δ ( t t )
the following equations are valid:
< x i ( t ) ξ j ( t ) > = 0 if t > t < x i ( t ) ξ j ( t ) > = 1 2 η i j
Now, we assume that the perturbation δ I is of order η (with η considered as a small parameter) and we expand Equation (26) up to the first order in η . We get
d I R S ( t ^ ) d t ^ = α I R S ( t ^ ) 1 I R S ( t ^ ) K β t ^ 2 1 t ^ I R S ( t ^ ) ( at the zero order ) d d t ^ δ I ( t ^ ) = α 2 α K I R S ( t ^ ) β t ^ 2 1 t ^ δ I ( t ^ ) + ξ ( t ^ ) ( at the first order )
Solution I R S ( t ^ ) is given by Equation (17) and the exact solution for δ I ( t ^ ) is easily obtained:
I R S ( t ^ ) = I 0 R S t ^ β exp ( ( 1 α / β ) 2 / σ ) exp ( ( t ^ α / β ) 2 / σ ) 1 + ( I 0 R S α / K ) exp ( ( 1 α / β ) 2 / σ ) 1 t ^ x β exp ( ( x α / β ) 2 / σ ) d x δ I ( t ^ ) = t β exp α ( t ^ 1 ) 1 / 2 β ( t ^ 2 1 ) 2 α G ( t ^ ) / K ( δ I ( 1 ) + 1 t ^ d x ξ ( x ) x β exp ( α ( x 1 ) + 1 / 2 β ( x 2 1 ) + 2 α G ( x ) / K ) ) with G ( t ^ ) 1 t ^ I R S ( x ) d x
Solution (33) allows to check the validity of Equation (31). Indeed, the correlation function δ I ( t ) -noise at the same time reads
< δ I ( t ^ ) ξ ( t ^ ) > = η 2
where Equation (30) and the properties of Dirac’s delta function have been taken into account. The equation for the second moment < δ I 2 ( t ^ ) > is easily obtained by multiplying both sides of the second equation of system (32) by 2 δ I ( t ) and by taking into account Equation (31). We get
d d t ^ < ( δ I ( t ^ ) ) 2 > = 2 α 2 α K I R S ( t ^ ) β t ^ 2 1 t ^ < ( δ I ( t ^ ) ) 2 > + η
The solution of Equation (35) provides the expression for the second moment
< ( δ I ( t ^ ) ) 2 > = t ^ 2 β exp ( 2 α ( t ^ 1 ) β ( t ^ 2 1 ) 4 α G ( t ^ ) / K ) ( < δ I ( 1 ) 2 > + η 1 t ^ x 2 β exp ( 2 α ( x 1 ) + β ( x 2 1 ) + 4 α G ( x ) / K ) d x )
From solution (33) we also get the correlation function for δ I
< δ I ( 1 ) δ I ( t ^ ) > = t ^ β exp ( α ( t ^ 1 ) β / 2 ( t ^ 2 1 ) 2 α G ( t ^ ) / K ) < ( δ I ( 1 ) ) 2 >
Figure 5 shows the simulation of 200 trajectories of Equation (26) for Italy. The employed algorithm for solving the Langevin equations is an Order-2 Stochastic Runge-Kutta like scheme (referred to as the RK-2 scheme) [24]. Regarding the discretisation in the numerical solution of the Langevin Equation (26) we use 500 time steps to build each stochastic path. The thick black curve is the numerical solution of the deterministic Equation (15). The noise intensity is η = 0.05 and the values of other parameters are shown in the caption of Figure 5.
Figure 6 and Figure 7 show the correlation functions < ( δ I ( t ^ ) ) 2 > and < δ I ( 1 ) δ ( t ^ ) > for Italy, first wave of infection by SARS-CoV-2. The values of the parameters are reported in the figure captions.
The Fokker–Planck equation, associated to Equation (26), describing the time evolution of the probability density function (PDF), p ( I , t ) , for the number of infectious people at time t ^ reads [25]
t ^ p ( I , t ^ ) = I [ μ ( t ^ ) p ( I , t ^ ) ] + D 2 I 2 p ( I , t ^ ) μ ( t ^ ) = α I 1 I K β t ^ 2 1 t ^ ; D = η 2 2
with μ denoting the drift coefficient and D the diffusion coefficient, respectively. Equation (38) must satisfies the initial condition
p ( I , 0 ) = δ ( I I 0 )
To get the PDF, we proceed as follows. First, we simulate paths or realisations of Equation (26) with the indicated parameters. The PDF is then computed numerically by constructing an histogram at each instant of time with a sufficiently number of realisations. The standard number used here is 22,000 paths. The values of the parameters are α = 1.5 , β = 0.1 , respectively, and for the PDF algorithm we use 800 bins in the construction of the corresponding histograms at the indicated times. The intensity of the noise, η = 0.05 , has been estimated according to Equation (A11) (see Appendix B). Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15 represent a set of screenshots at the indicated times t ^ = t / t L .

4. The (SISIhrhdh)L-Model

As can be seen, according to the S I S -model, after the lockdown measures the number of infectious people starts to increase again. Hospitals and health institutions play a crucial role in hindering the spread of the Coronavirus. In this Section, we propose a model which accounts for people who are only traced back to hospitalised infectious individuals. According to our chemical-type reactions approach, 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 [26,27,28,29]) 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 catalyser in the hospital recovery process [8,9]. We propose the following model:
S + I σ 2 I I γ + c ( t ) S I + b k 1 I h k r r h + b I h k d d h + b r h γ 1 S
where the hypothesis that an individual acquires immunity, after having contracted the Coronavirus and being recovered, is not adopted. In Figure 16, b denotes the number of available hospital beds, 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, respectively. According to Figure 16, people, once recovered, are subjected to the same existing lockdown measures as any other people.
Concerning the recovered people, we would like to make clear the following. R stands for the total number of the recovered people (i.e., the number of recovered people previously hospitalised, plus the number of the asymptomatic people, plus the infected people who have been recovered without being previously hospitalised). However, the natural question is: how can we count R and compare this variable with real data? The current statistics, produced by the Ministries of Health of various Countries, concern the number of people released from hospitals. Apart from Luxembourg (where almost the entire population has been subjected to the COVID-19 test), no other Countries are in a condition as 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 the theoretical predictions for R with real data. We then proceed by adopting approximations and by establishing the differential equations where the solution may realistically be subjected to experimental verification. More specifically, first, we assume that R is given by three contributions:
R = r h + r A + r I
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 . Hence, 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. Finally, due to lack of reliable statistics, we are forced to limit ourselves to consider the (very) simplified case
R r h
Of course, we have
I h + b = C h = c o n s t . where C h = Total hospital s capacity
The dynamical equations for the entire process are then:
d d t ^ S = σ S N T o t . I + γ I + γ 1 r h + β t ^ 2 1 t ^ ( I + r h ) d d t ^ I = σ S N T o t . I γ I k 1 I ( C h I h ) β t ^ 2 1 t ^ I d d t ^ I h = k 1 I ( C h I h ) k r I h k d I h d d t ^ r h = k r I h γ 1 r h d d t ^ d h = k d I h
where, for simplicity, the average recovery time delay and the average death time delay have been neglected [8]. In this case I stands for the infectious individuals not hospitalised. From system (44) we get the following conservation law
S + I + I h + r h + d h = N T o t . = c o n s t .

The Deterministic ( S I S I h ) L -Model

To simplify as much as possible the set of O.D.E.s (44), we adopt several hypotheses that will not compromise the validity of our model. First, we assume that that S + I + I h N T o t . = c o n s t . Secondly, let γ γ 1 . Finally, we take into account the current Belgian hospital-protocol: “Only the seriously sick people are hospitalised. The remaining infectious individuals have to be sent home and they must be subjected to quarantine measures”. Hence, I I h and the total number of recovered people, R, is much larger than the total number of recovered people, previously hospitalised (i.e., R r h ). Under these assumptions, the model simplifies to
S + I σ 2 I I γ + c ( t ) S
I + b k 1 I h
with I h + b = C h . Hence, under these assumptions, after hospitalisation, individuals will be removed from the disease, either due to immunisation (e.g., due to vaccination or special health care received) or due to death. The governing O.D.E.s, associated to the model (46), read
d d t ^ S σ S N T o t . I + γ I + β t ^ 2 1 t ^ I d d t ^ I = σ S N T o t . I γ I k 1 I ( C h I h ) β t ^ 2 1 t ^ I d d t ^ I h k 1 I ( C h I h )
The recovered people r h and the deceased people d h may be obtained by solving, respectively, the following O.D.E.s
d d t ^ r h k r I h ; d d t ^ d h = k d I h
From system (48), we obtain
d d t ^ I = α 1 I 1 I K 1 ρ I I h β t ^ 2 1 t ^ I d d t ^ I h = k 1 C h I k 1 I I h with α 1 = σ 1 γ σ k 1 σ C h ; K 1 = N T o t . 1 γ σ k 1 σ C h ; ρ = σ N T o t . k 1
In absence of the lockdown measures ( β = 0 ), we have the following scenarios: I 0 > 0
( i ) if C h < C h C r i t . σ k 1 1 γ σ , the equilibrium state I = K 1 is stable ( i i ) if C h > C h C r i t . σ k 1 1 γ σ , the equilibrium state I = 0 is stable
In words:
• For case (i), there will be a proper epidemic outbreak with an increase of the number of the infectious people;
• For case ( i i ), independently of the initial size of the susceptible population, the disease can never cause a proper epidemic outbreak.
This result highlights the crucial role of the hospitals and the health care institutes:
If the threshold of the hospital capacities exceeds a lower limit, the spread of the Coronavirus tends to decrease over time, and the stable solution corresponds to zero infectious individuals.
Figures 40 and Figure 17 illustrate the situation. Notice that, with the values of parameters reported in the corresponding figure captions, C h C r i t . = 18,434 .
A similar analysis leading to Equation (51), allowing the calculation of the critical threshold for France and US hospital capacities, may also be performed. By summarising, the ( S I S I h ) L -model shows first the crucial role of hospital capacity and second highlights its limits. For instance, for the Italian situation, to obtain a relevant dampening effect of the COVID-19 infection, the capacity of hospitals in Italy would have to increase by about 4 times its current value (it is huge). Hence the need to combine and coordinate the two actions at the same time: to increase the hospitals’ capacity as much as possible and to distribute effective vaccines. The current analysis is mainly addressed to Countries that do not have the possibility to buy and massively distribute vaccines (such as some African countries). In this case, the role of the hospitals becomes crucial and basically it represents the only real remedy to stop the of the pandemic.

5. The Stochastic ( S I S I h ) L Model

If the dynamic is subjected to white noise and related stochastic equations are written as
d d t ^ I = α 1 I 1 I K 1 ρ I I h β t ^ 2 1 t ^ I + ξ 1 ( t ^ ) d d t ^ I h = k 1 C h I k 1 I I h + ξ 2 ( t ^ )
where
< ξ i ( t ) > = 0 ( i = 1 , 2 ) < ξ i ( t ) ξ j ( t ) > = η i j δ i j δ ( t t ) with η 12 = η 21
with δ i j denoting the Kronecker delta. The statistical properties of this processes, i.e., < δ I ( t ^ ) δ I ( t ^ ) > , < δ I ( t ^ ) δ I h ( t ^ ) > and < δ I h ( t ^ ) δ I h ( t ^ ) > , may be obtained, firstly, by determining the reference state. This state satisfies the following O.D.E.s
d d t ^ I R S = α 1 I R S 1 I R S K 1 ρ I R S I h R S β t ^ 2 1 t ^ I R S d d t ^ I h R S = k 1 C h I R S k 1 I R S I h R S
Let us choose, for example, I h R S = C h , then
I R S ( t ^ ) = I 0 R S exp ( ( 1 α 2 / β ) 2 / σ ) t ^ β exp ( ( t ^ α 2 / β ) 2 / σ ) 1 + ( I 0 R S α 2 / K 2 ) exp ( ( 1 α 2 / β ) 2 / σ ) 1 t ^ x β exp ( ( x α 2 / β ) 2 / σ ) d x α 2 σ 1 γ σ C h N T o t . ; K 2 = N T o t . 1 γ σ C h N T o t .
Successively, we have to find the solution of the Fokker–Planck equation, p ( I , I h , t ^ ) , associated to the system (52), which in this case reads
d d t ^ p = I ( μ 1 p ) I h ( μ 2 p ) + D 11 I 2 p + D 12 I I h p + D 22 I h 2 p
with
μ 1 = α 1 I 1 I K 1 ρ I I h β t ^ 2 1 t ^ I ; μ 2 = k 1 C h I k 1 I I h D 11 = 1 2 ( η 11 + η 12 ) ; D 12 = η 12 ( η 11 + η 22 ) ; D 22 = 1 2 ( η 22 + η 12 )
Finally, with the expression for p ( I , I h , t ) we may compute the correlation functions < δ I ( t ^ ) δ I ( t ^ ) > , < δ I ( t ^ ) δ I h ( t ^ ) > and < δ I h ( t ^ ) δ I h ( t ^ ) > with
δ I = I I R S ; δ I h = I h I h R S
However, this is a quite long and tedious procedure. We may circumvent the obstacle by recalling that, as we did in Section 3, it is possible to obtain such correlation functions in a direct way, by using theorem (30) and (31). For this, we have to determine the dynamical O.D.E.s for the fluctuations δ I ( t ^ ) and δ I h ( t ^ ) around the reference state ( I R S , I h R S ). We get
d d t ^ δ I = α 1 ρ C h 2 α 1 K 1 I R S β t ^ 2 1 t ^ δ I ρ I R S δ I h + ζ 1 ( t ^ ) d d t ^ δ I h = k 1 I R S ( t ) δ I h + ζ 2 ( t ^ )
The correlation functions < δ I ( t ^ 0 ) δ I ( t ^ ) > and < δ I h ( t ^ 0 ) δ I h ( t ^ ) > are easily obtained
< δ I ( t ^ 0 ) δ I ( t ^ ) > = t ^ t ^ 0 β exp ( H ( t ^ ) ) exp 2 α 1 K 1 G ( t ^ ) ( < δ I ( t ^ 0 ) 2 > ρ < δ I ( t ^ 0 ) δ I h ( t ^ 0 ) > t ^ 0 t ^ t ^ 0 s β I R S ( s ) exp ( H ( s ) ) exp ( 2 α 1 K 1 k 1 ) G ( s ) d s ) < δ I h ( t ^ 0 ) δ I h ( t ^ ) > = < δ I h ( t ^ 0 ) 2 > exp ( k 1 G ( t ^ ) ) with G ( t ^ ) t ^ 0 t ^ I R S ( s ) d s ; H ( t ^ ) ( α 1 ρ C h ) ( t ^ t ^ 0 ) 1 / 2 β ( t ^ 2 t ^ 0 2 )
In addition, by using again Equation (59), we obtain the O.D.E.s of the correlation functions < ( δ I ( t ^ ) ) 2 > , < ( δ I h ( t ^ ) ) 2 > , and < δ I ( t ^ ) δ I h ( t ^ ) >
d d t ^ < ( δ I ( t ^ ) ) 2 > = 2 α 1 ρ C h 2 α 1 K 1 I R S β t ^ 2 1 t ^ < ( δ I ( t ^ ) ) 2 > 2 ρ I R S < δ I ( t ^ ) δ I h ( t ^ ) > + η 11 ( t ^ ) d d t ^ < ( δ I h ( t ^ ) ) 2 > = 2 k 1 I R S < ( δ I h ( t ^ ) ) 2 > + η 22 d d t ^ < δ I ( t ^ ) δ I h ( t ^ ) > = α 1 ρ C h 2 α 1 K 1 + k 1 I R S β t ^ 2 1 t ^ < δ I ( t ^ ) δ I h ( t ^ ) > ρ I R S < ( δ I h ( t ^ ) ) 2 > + η 12 ( t ^ )
where Equation (31) have been taken into account. The solutions of system (61) read
< ( δ I h ( t ^ ) ) 2 > = exp ( 2 k 1 G ( t ^ ) ) < ( δ I h ( t ^ 0 ) ) 2 > + η 22 t ^ 0 t ^ exp ( 2 k 1 G ( s ) ) d s < δ I ( t ^ ) δ I h ( t ^ ) > = t ^ t ^ 0 β exp ( H ( t ^ ) ) exp ( ( 2 α 1 K 1 + k 1 ) G ( t ^ ) ) ( < δ I ( t ^ 0 ) δ I h ( t ^ 0 ) > t ^ 0 t ^ t ^ 0 s β exp ( H ( s ) ) exp ( ( 2 α 1 K 1 + k 1 ) G ( s ) ) η 12 + ρ I R S ( s ) < ( δ I h ( s ) ) 2 > d s ) < ( δ I ( t ^ ) ) 2 > = t ^ t ^ 0 2 β exp ( 2 H ( t ^ ) ) exp ( 4 α 1 K 1 G ( t ^ ) ) ( < ( δ I ( t ^ 0 ) ) 2 > t ^ 0 t ^ t ^ 0 s 2 β exp ( 2 H ( s ) ) exp ( 4 α 1 K 1 G ( s ) ) η 11 2 ρ I R S ( s ) < δ I ( s ) δ I h ( s ) > d s )
The value of C h (the total Italian hospitals’ capacity) may be obtained by making reference to the data published in [30]. More specifically, in 2017, when there were 518 public hospitals and 482 accredited private ones, in Italy there were 151,646 beds for ordinary hospitalisation in public hospitals ( 2.5 per 1000 inhabitants) and 40,458 in private ones ( 0.7 per 1000 inhabitants), for a total of over 192 thousand beds ( 3.2 per 1000 inhabitants).
The number of public and private beds intended for intensive care was 5090 (a number very close to the 5100 cited by the newspapers these days), about 8.42 per 100,000 inhabitants [30].
Figure 18 shows the correlation function < δ I h ( t ) δ I h ( t ) > , solution of the second O.D.E. in system (61). As can be seen, this is a typical correlation function at equilibrium for a system which is subjected to random fluctuations. This is not surprising as our reference state correspond to the maximum capacity of the hospitals (i.e., I h R S = C h ) so, fluctuations at equilibrium are the only possible ones. The three correlation functions < δ I ( t ) δ I ( t ) > , < δ I h ( t ) δ I h ( t ) > , and < δ I ( t ) δ I h ( t ) > , solutions of system (61), are shown in Figure 19.

6. Comparison between Theoretical Predictions and Real Data for USA and France

To perform a comparison with real data, it is convenient to rewrite Equation (26) in dimensional form
d I d t = α ˜ I 1 I K β t 2 t L 2 t L 2 t I + η ˜ ξ ˜ ( t )
where α ˜ α / t L and the white noise ξ ˜ ( t ) is defined by the relations
< ξ ˜ ( t ) > = 0 < ξ ˜ ( t ) ξ ˜ ( t ) > = δ ( t t )
The order of magnitude for the intensity of the noise η ˜ may be derived by using Equation (A11) reported in Appendix B:
η ˜ K 1 / 2
References [31,32,33,34,35], report links to the official sites for the number of infected individuals in US and France, respectively. As the data show, France is currently subjected to a second wave of Coronavirus while the USA data may induce to think that they are in a full second (or third) wave. However, looking at the behaviour of the infectious curve we may also argue that, in agreement with ref. [4], USA as a whole is not in a second (or third) wave because the first wave never really stopped. The virus is simply into new populations or resurgent in places that let down their guard too soon. So, we are interested in analysing both of these scenarios.
The number of infectious individuals by SARS-CoV-2 in USA and France is given in Figure 20 and Figure 21, respectively. Figure 22 and Figure 23 show a comparison between the theoretical predictions (blue line) and real data (black dots) for USA. Figure 22 refers to the hypothesis that US is still in the first wave of Coronavirus, whereas Figure 23 has been obtained by assuming that US population is currently subjected to the third wave of Coronavirus, respectively. Figure 24 and Figure 25 show a comparison between the theoretical predictions (blue line) and real data (black dots) for France. The values of parameters for both cases, USA and France, are reported in the corresponding figure captions.
Now, let us deal with the stochastic processes. Figure 26 and Figure 27 illustrate the comparison between the theoretical predictions (blue lines) and real data (black dots) for France concerning the first and the second waves of SARS-CoV-2, respectively. The intensity of the noise has been estimated by using Equation (65). As we can see, the predictions of our model are in a fairly good agreement with real data.
In the case of USA, we have analysed two possible scenarios:
(1) USA is still subjected to the first wave of Coronavirus infection;
(2) USA is in the second or third wave of SARS-CoV-2 infection.
Figure 28 and Figure 29 illustrate these two scenarios. More specifically, Figure 28 refers to scenario (1). In order to obtain the best possible agreement with real data, we have increased, arbitrarily, the intensity of the noise up to a value that is 45 times higher than that estimated by Equation (65) (we set η ˜ = 0.009 ). However, we did not reach this objective in a satisfactory way.
Figure 29 refers to scenario (2). In this case the intensity of the noise corresponds to the value estimated by Equation (65) (we have η ˜ = 0.0002 ). Here, we get a fairly good agreement with real data. Figure 30, Figure 31, Figure 32 and Figure 33 show the probability distribution functions for USA-first wave at the beginning of the lockdown measures (55 days) and after 90 days, 200 days, and 250 days, respectively.
Finally, our results lead to two possible interpretations:
(a) USA is experiencing a second (or even third) wave of Coronavirus infection;
(b) USA is still subjected to the first wave of infection. However, given the heterogeneity and vastness of its territory, the variables of the system (i.e., S, I, R, and D) must depend on time as well as on space. In this case, the kind of equations governing the evolution of the compartments must be Stochastic Partial Differential Equations (S.P.D.E.s).
In our opinion, scenario (b) is the correct interpretation.

7. Conclusions

Waiting for the mass production of vaccines, we studied two scenarios by means of two mathematical models: the ( S I S ) L -model and the ( S I S I h r h d h ) L -model. In this work we adopted a kinetic-type reactions approach, where the compartments S, I, and R are seen as chemically active molecules. The analysis of both models have been carried out by introducing (at the best of our knowledge for the first time in the literature) a new complete mathematical basis and by determining the intensity of the stochastic noise by methods of statistical mechanics. The first scenario, described by the ( S I S ) L -model, refers to a situation in which the role of hospitals is not taken into account. The second model (the ( S I S I h r h d h ) L -model) integrates the contribution provided by the hospitals and health care institutes. The dynamics are governed by stochastic differential equations, which are suitable for describing realistic situations. We compute the relevant correlation functions as well as the probability density functions. In the S I S L model, we determined the minimum value of infectious individuals that must be reached after the lockdown measures so that, once these measures are removed, during a time interval Δ t ^ , the number of infectious still stay below a set threshold.
The ( S I S I h r h d h ) L -model shows the crucial role played by the hospitals. More specifically, we showed that the health care institutions can influence the outcome of the outbreak, limiting, and even dampening, the of the Coronavirus. By way of example, we dealt with cases for the United States and France. The comparison between the theoretical predictions with real data confirmed the validity of our models. In particular, for US, we have examined two possible scenarios: USA still subjected to the first wave of infection by Coronavirus and USA in the second (or third) wave of the SARS-CoV-2 infection. We concluded that, for a so vast and heterogeneous country, spatial dependance of the variables can no longer be ignored and the suitable equations governing the evolution of the different compartments, S, I, R, and D = d e c e a s e d p e o p l e , should be Stochastic Partial Differential Equations-type (S.P.D.E.s).
We would like to stress that the aim of this work is to grasp relevant information by adopting a very simple model, the simplest possible one. Once the key effects are understood, it is not difficult to add supplementary details by including extra compartments, e.g., for studying the effect of temporary immunity generated by COVID-19 infection, the effect of the quarantine measures, implications on age-related results, etc. Of course, in these cases, the price to pay is that, contrary to the current work, the analytical analysis can only be performed in a very limited way and the study must be carried out mainly by numerical schemes.

Author Contributions

G.S.: Conceptualisation, formal analysis, calculations, writing & editing, F.M.: numerics PDF, writing & editing, and P.N.: validation, writing & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting reported results: See the list of References below.

Acknowledgments

G.S. is grateful to M. Malek Mansour, of the Université Libre de Bruxelles (ULB), for his useful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Exact Solution of Richards’ Differential Equation with the Lockdown Term

Richards differential equation is an extension of the logistic differential equation, allowing for more flexible S-shaped curves. This curve satisfies the O.D.E. [36]
d I d t = α I 1 I K ( t ) ν
with ν > 0 (with ν > 0 ) denoting a parameter affecting the point where the asymptote maximum growth occurs. In presence of therapy-induced death rate, Richards’ equation takes the form
d I d t = α I 1 I K ( t ) ν c ( t ) I
The goal is to determine the analytic solution of Equation (A2). We describe the procedure to find the exact solution of Equation (A2) since it is not easy to find it in the literature. Let us first rewrite Equation (A2) as
1 I d I d t = α 1 I K ( t ) ν c ( t )
Successively, we introduce new functions n ( t ) and C ( t ) defined as
n ( t ) ln ( I ( t ) ) + C ( t ) where C ( t ) t 0 t c ( s ) d s
with t 0 denoting the initial time. We get
d n d t = 1 I d I d t + d C d t = 1 I d I d t + c ( t ) and d d t ( n ( t ) α t ) = α exp ( ν ( n ( t ) α t ) ) exp ( ν ( C ( t ) α t ) ) K ( t ) ν
By setting
m ( t ) n ( t ) α t
we obtain
m 0 m d m exp ( ν m ) = α t 0 t exp ( ν ( C ( s ) α s ) K ( s ) ν d s
which, after integration, reads
exp ( m ( t ) ) = c 0 + ν α t 0 t exp ( ν ( C ( s ) α s ) K ( s ) ν d s 1 / ν
where c 0 is a constant independent of time. Finally, we get
I ( t ) = exp ( α t C ( t ) ) c 0 + ν α t 0 t exp ( ν ( C ( s ) α s ) K ( s ) ν d s 1 / ν
Solution (A9) may be cast in the following, more convenient, form
I = I 0 exp ( α ( t t 0 ) C ( t ) ) 1 + ν α I 0 ν exp ( α ν t 0 ) t 0 t ( exp ( ν ( C ( s ) α s ) ) / K ( s ) ν d s 1 / ν
with I 0 denoting the number of population I at time t 0 ( I 0 I ( t 0 ) ).

Appendix B. A Rough Estimate of the Intensity of Noise

The intensity of noise for the entire process may be estimated by adopting the concept of Grand Ensemble in statistical mechanics. More specifically, we may establish the following correspondences: Quantity of SARS-CoV-2 corresponds to the Energy of the system, the Quantity of SARS-CoV-2 that can be absorbed or released due to a change of the infectious individuals models the chemical potential ( μ c h . ), and the surface of a country corresponds to the volume of the system, respectively. The temperature of a system (T) is mimicked by the surrounding environment, which is more or less favourable to the spread of Coronavirus depending on whether particular prevention measures are adopted (e.g., use of masks, maintaining a safe distance, tracking of the infection, adoption of lockdown, quarantine and isolation measures, etc.). Hence, in this scenario:
(i) Temperature equilibrium (“ T C o u n t r y = T E n v . ”) stands for environment equilibrium;
(ii) Identical chemical potentials ( μ C o u n t r y c h . = μ E n v . c h . ) stands for equal variation of the quantity of Coronavirus with respect to a change of the infectious people;
(iii) The surface of a country (corresponding to the volume of the system) is kept constant;
(iv) The total number of infectious individuals (corresponding to the “number of particles”) and the quantity of the SARS-CoV-2 in a country (corresponding to the “energy of the system”) are not fixed.
Under the above conditions, we may consider that the system is governed by the Grand Ensemble in statistical mechanics. So, according to this statistics, the relative fluctuations of the number of the infectious people behaves as K C o u n t r y 1 / 2 [37], with K C o u n t r y denoting the capacity of the country’s population. It is easily checked that this implies that the intensity of the noise in Equation (26) is of the order of
η t L K C o u n t r y 1 / 2
with t L denoting the time when the lockdown measures have been applied. We may object that this reasoning is based on equilibrium conditions. However, we would like to point out that the aim of this section is only to provide a (rough) estimate of the order of magnitude of the noise.

References

  1. Available online: https://www.euronews.com/2020/11/18/is-europe-having-a-COVID-19-second-wave-country-by-country-breakdown (accessed on 18 November 2020).
  2. Available online: https://www.ilsole24ore.com/art/covid-analogie-e-differenze-prima-e-seconda-ondata-AD2Igxu?refresh_ce=1 (accessed on 11 October 2020).
  3. Available online: https://www.bbc.com/news/world-us-canada-54966531 (accessed on 12 December 2020).
  4. Available online: https://theconversation.com/the-us-isnt-in-a-second-wave-of-coronavirus-the-first-wave-never-ended-141032 (accessed on 30 June 2020).
  5. Available online: https://www.thelocal.fr/20201028/why-the-second-wave-of-covid-infections-in-france-could-be-worse-than-the-first (accessed on 31 October 2020).
  6. Available online: https://www.iene.mediaset.it/2020/news/coronavirus-italia-coprifuoco-seconda-ondata_929030.shtml (accessed on 16 November 2020).
  7. May, R.M.; Anderson Roy, M. Infectious Diseases of Humans: Dynamics and Control; Oxford University Press: Oxford, UK, 1991; ISBN 0-19-854040-X. [Google Scholar]
  8. Sonnino, G.; Peeters, P.; Nardone, P. Modelling the Coronavirus Second Wave in Presence of the Lockdown and Quarantine Measures. In Proceedings of the CCS2020-Conference on Complex System, Singapore, 27–31 July 2020; Book of Abstracts. p. 278. Available online: https://zenodo.org/record/4419178#.YKki1JMzbZ0 (accessed on 8 May 2021). [CrossRef]
  9. Sonnino, G.; Peeters, P.; Nardone, P. Modelling the of the SARS-CoV-2 in Presence of the Lockdown and Quarantine Measures by a Kinetic-Type Reactions Approach. medRxiv 2021. Available online: https://www.medrxiv.org/content/10.1101/2021.01.14.21249797v1 (accessed on 15 January 2021). [CrossRef]
  10. Carli, R.; Cavone, G.; Epicoco, N.; Scarabaggio, P.; Dotoli, M. Model predictive control to mitigate the COVID-19 outbreak in a multi-region scenario. Annu. Rev. Control 2020, 50, 373–393. [Google Scholar] [CrossRef]
  11. Calafiore, G.C.; Novara, C.; Possieri, C. A Time-Varying SIRD Model for the COVID-19 Contagion in Italy. Annu. Rev. Control 2020, 50, 361–372. [Google Scholar] [CrossRef] [PubMed]
  12. Estrada, E. COVID-19 and SARS-CoV-2. Modeling the present, looking at the future. Phys. Rep. 2020, 869, 1–51. [Google Scholar] [CrossRef] [PubMed]
  13. Sonnino, G.; Mora, F.; Nardone, P. A Stochastic Compartimental Model for COVID-19. medRxiv 2020. Available online: https://www.medrxiv.org/content/10.1101/2020.12.03.20243584v1 (accessed on 4 December 2020). [CrossRef]
  14. Piovani, P. Coronavirus, il pediatra: Questa Epidemia Durerà tre Mesi. Il Messaggero. 24 March 2020. Available online: https://www.ilmessaggero.it/italia/coronavirus_italia_quanto_durera_quando_finira_picco_news-5129105.html (accessed on 24 March 2020).
  15. Moehlis, J.M. An SIS Model-APC/EEB/MOL 514 Tutorial 4: Seasonal Epidemic Models. Available online: https://sites.me.ucsb.edu/~moehlis/APC514/tutorials/tutorial_seasonal/node2.html (accessed on 14 October 2002).
  16. Sonnino, G.; Nardone, P. Dynamics of the COVID-19-Comparison between the Theoretical Predictions and the Real Data, and Predictions about Returning to Normal Life. Ann. Clin. Med Case Rep. 2020, 4, 1–21. [Google Scholar] [CrossRef]
  17. Sonnino, G. The COVID-19-The Infectious Disease Caused by the Latest Discovered Coronavirus (SARS-CoV-2). European Commission. ARES(2020)1530456. Available online: https://europa.eu/european-union/contact/write-to-us_en (accessed on 16 March 2020).
  18. Mathur, N.; Shaw, G. An Empirical Model on the Dynamics of COVID-19 Spread in Human Population; Department of Theoretical Physics, Tata Institute of Fundamental Research: Mumbai, India, 2020; Available online: https://arxiv.org/pdf/2008.06346.pdf (accessed on 13 August 2020).
  19. Coullet, P. The COVID-19 Epidemic as a Simple Dynamical System. Comptes Rendus Mécanique, 0000, 1, No 0, p. 000-000 DOI Unassigned Yet, Under Review. Available online: https://comptes-rendus.academie-sciences.fr/mecanique/ (accessed on 27 June 2020).
  20. Baras, F.; Malek Mansour, M. Problème d’Actualité de la Physique II-Dynamiques des Fluides: Aspects Microscopiques et Simulations; Lecture Note at the Université Libre de Bruxelles (U.L.B.); Université Libre de Bruxelles: Bruxelles, Belgium, 1991. [Google Scholar]
  21. Malek Mansour, M. Processus Stochastiques en Physique; Lecture Note at the Universite’ Libre de Bruxelles (U.L.B.); Université Libre de Bruxelles: Bruxelles, Belgium, 2014. [Google Scholar]
  22. Gardiner, C. Stochastic Methods. A Handbook for the Natural and Social Sciences, 4th ed.; Series in Synergetics (13); Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  23. Van Kampen, N.G. Stochastic Processes in Physics and Chemistry, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  24. Roberts, A.J. Model Emergent Dynamics in Complex Systems; Section 20.4; SIAM: Philadelphia, PA, USA, 2015; ISBN 978161173556. [Google Scholar]
  25. Risken, H. The Fokker-Planck-Equation. In Methods of Solution and Applications, 2nd ed.; Springer Series in Synergetics 18; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
  26. Michaelis, L.; Menten, M.L. Die Kinetik der Invertinwirkung. Biochem. Z 1913, 49, 333–369. [Google Scholar]
  27. Victor, H. Lois Générales de l’Action des Diastases. In Victor Henri. Whonamedit? Hermann: France, Paris, 1903. [Google Scholar]
  28. Bharath, S. Explicit treatment of non Michaelis-Menten and atypical kinetics in early drug discovery. ChemMedChem 2020. [Google Scholar] [CrossRef]
  29. Bharath, S. Words of advice: Teaching enzyme kinetics. FEBS J. 2020. [Google Scholar] [CrossRef]
  30. Available online: https://www.agi.it/fact-checking/news/2020-03-06/coronavirus-posti-letto-ospedali-7343251/ (accessed on 6 March 2020).
  31. Center for Disease Control and Prevention (CDC). CDC COVID Data Tracker, Maps, Charts, and Data Provided by the CDC-Coronavirus Disease 2019 (COVID-19). Available online: https://covid.cdc.gov/covid-data-tracker/#cases_casesper100klast7days (accessed on 23 November 2020).
  32. Available online: https://www.nytimes.com/interactive/2020/us/states-reopen-map-coronavirus.html (accessed on 23 November 2020).
  33. Available online: https://www.data.gouv.fr/en/datasets/chiffres-cles-concernant-lepidemie-de-covid19-en-france/ (accessed on 23 November 2020).
  34. Available online: https://www.rainews.it/ran24/speciali/2020/covid19/europe.php (accessed on 23 November 2020).
  35. Available online: https://www.gouvernement.fr/info-coronavirus/les-actions-du-gouvernement (accessed on 23 November 2020).
  36. Richards, F.J. A Flexible Growth Function for Empirical Use. J. Exp. Bot. 1959, 10, 290–300. [Google Scholar] [CrossRef]
  37. Reif, F. Fundamentals of Statistical and Thermal Physics; Waveland Press, Inc.: Long Grove, IL, USA, 2009; ISBN 978-1-57766-612-7. [Google Scholar]
Figure 1. Solution of Equation (15) for Italy, first wave of infection by SARS-CoV-2. The black and red dotted lines refer to solutions of Equation (15) with and without application of lockdown measures, respectively. The values of parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Figure 1. Solution of Equation (15) for Italy, first wave of infection by SARS-CoV-2. The black and red dotted lines refer to solutions of Equation (15) with and without application of lockdown measures, respectively. The values of parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Mathematics 09 01221 g001
Figure 2. Schematic representation of the second wave of SARS-CoV-2 infection. The Lockdown measures are applied during the period t ^ 0 t ^ 1 . At the end of the confinement measures, the number of infectious dropped to I 1 . During the time interval Δ t ^ = t ^ 2 t ^ 1 , the lockdown measures are not applied and the infection by SARS-CoV-2 (starts growing again ( second wave of Coronavirus infection) until the preset value I ˜ is reached, which corresponds to I ( t ^ 2 ) . During the whole time interval Δ t ^ , the number of infectious individuals still stay below a set threshold I ˜ if I 1 satisfies condition (25).
Figure 2. Schematic representation of the second wave of SARS-CoV-2 infection. The Lockdown measures are applied during the period t ^ 0 t ^ 1 . At the end of the confinement measures, the number of infectious dropped to I 1 . During the time interval Δ t ^ = t ^ 2 t ^ 1 , the lockdown measures are not applied and the infection by SARS-CoV-2 (starts growing again ( second wave of Coronavirus infection) until the preset value I ˜ is reached, which corresponds to I ( t ^ 2 ) . During the whole time interval Δ t ^ , the number of infectious individuals still stay below a set threshold I ˜ if I 1 satisfies condition (25).
Mathematics 09 01221 g002
Figure 3. I 1 versus I ˜ and Δ t ^ in the range 50,000 I ˜ K I T .
Figure 3. I 1 versus I ˜ and Δ t ^ in the range 50,000 I ˜ K I T .
Mathematics 09 01221 g003
Figure 4. As shown in the figure, setting the value of I ˜ = 50,000, the maximum value of I 1 is I 1 = 9500 for Δ t = t L , and I 1 = 700 for Δ t = 2 t L , respectively.
Figure 4. As shown in the figure, setting the value of I ˜ = 50,000, the maximum value of I 1 is I 1 = 9500 for Δ t = t L , and I 1 = 700 for Δ t = 2 t L , respectively.
Mathematics 09 01221 g004
Figure 5. Stochastic equation for Italy. Sample paths of Equations (26) and (27) for Italy-first wave of infection by SARS-CoV-2 with η = 0.05 and for 200 realisations. The values of the other parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively. The black thick curve is the numerical solution of the deterministic Equation (15).
Figure 5. Stochastic equation for Italy. Sample paths of Equations (26) and (27) for Italy-first wave of infection by SARS-CoV-2 with η = 0.05 and for 200 realisations. The values of the other parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively. The black thick curve is the numerical solution of the deterministic Equation (15).
Mathematics 09 01221 g005
Figure 6. Correlation Function < ( δ I ( t ^ ) ) 2 > for Italy, first wave of infection by SARS-CoV-2. This correlation function corresponds to the solution of Equation (35). The values of the parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Figure 6. Correlation Function < ( δ I ( t ^ ) ) 2 > for Italy, first wave of infection by SARS-CoV-2. This correlation function corresponds to the solution of Equation (35). The values of the parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Mathematics 09 01221 g006
Figure 7. Correlation Function for Infectious people < δ I ( 1 ) δ ( t ^ ) > . Plot of Equation (37) with I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Figure 7. Correlation Function for Infectious people < δ I ( 1 ) δ ( t ^ ) > . Plot of Equation (37) with I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Mathematics 09 01221 g007
Figure 8. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 1.129 . The values of the parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Figure 8. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 1.129 . The values of the parameters are I 0 = 5000 , K = 150,000, α = 1.5 , and β = 0.1 , respectively.
Mathematics 09 01221 g008
Figure 9. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 2.0304 .
Figure 9. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 2.0304 .
Mathematics 09 01221 g009
Figure 10. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 3.0647 .
Figure 10. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 3.0647 .
Mathematics 09 01221 g010
Figure 11. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 4.0397 .
Figure 11. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 4.0397 .
Mathematics 09 01221 g011
Figure 12. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 5.0131 .
Figure 12. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 5.0131 .
Mathematics 09 01221 g012
Figure 13. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 6.0375 .
Figure 13. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 6.0375 .
Mathematics 09 01221 g013
Figure 14. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 7.0656 .
Figure 14. Probability density function (PDF), solution of Equation (36) for Italy, first wave of infection by SARS-CoV-2. Screenshoot at time t ^ = 7.0656 .
Mathematics 09 01221 g014
Figure 15. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 7.9942 .
Figure 15. Probability density function (PDF), solution of Equation (36). Screenshoot at time t ^ = 7.9942 .
Mathematics 09 01221 g015
Figure 16. Infectious people and infectious hospitalised people vs. time— C h > C H C r i t . . In this case, there is an epidemic outbreak with an increase of the number of the infectious people. The values of the parameters are: β = 0 , k 1 = 0.00001 , ρ = 0.0001 , K 1 = 150,000, γ = 0.0001 , σ = 2.7651 , and C h = 5100 .
Figure 16. Infectious people and infectious hospitalised people vs. time— C h > C H C r i t . . In this case, there is an epidemic outbreak with an increase of the number of the infectious people. The values of the parameters are: β = 0 , k 1 = 0.00001 , ρ = 0.0001 , K 1 = 150,000, γ = 0.0001 , σ = 2.7651 , and C h = 5100 .
Mathematics 09 01221 g016
Figure 17. Infectious people and infectious hospitalised people vs. time— C h > C H C r i t . . Independently of the initial size of the susceptible population, the disease can never cause a proper epidemic outbreak. The values of the parameters are: β = 0 , k 1 = 0.00001 , ρ = 0.0001 , K 1 = 150,000, γ = 0.0001 , σ = 2.7651 , and C h = 18,436.
Figure 17. Infectious people and infectious hospitalised people vs. time— C h > C H C r i t . . Independently of the initial size of the susceptible population, the disease can never cause a proper epidemic outbreak. The values of the parameters are: β = 0 , k 1 = 0.00001 , ρ = 0.0001 , K 1 = 150,000, γ = 0.0001 , σ = 2.7651 , and C h = 18,436.
Mathematics 09 01221 g017
Figure 18. Correlation function < ( δ I h ( t ) ) 2 > for Italy, first wave of infection by SARS-CoV-2. This correlation function is the solution of the second equation of system (61). Having chosen our reference state I h R S = C h , we get a typical correlation function for random fluctuations at equilibrium.
Figure 18. Correlation function < ( δ I h ( t ) ) 2 > for Italy, first wave of infection by SARS-CoV-2. This correlation function is the solution of the second equation of system (61). Having chosen our reference state I h R S = C h , we get a typical correlation function for random fluctuations at equilibrium.
Mathematics 09 01221 g018
Figure 19. Correlation functions < δ I ( t ) δ I ( t ) > , < δ I h ( t ) δ I h ( t ) > , and < δ I ( t ) δ I h ( t ) > . These three correlation functions are obtained by solving numerically system (61).
Figure 19. Correlation functions < δ I ( t ) δ I ( t ) > , < δ I h ( t ) δ I h ( t ) > , and < δ I ( t ) δ I h ( t ) > . These three correlation functions are obtained by solving numerically system (61).
Mathematics 09 01221 g019
Figure 20. Number of infectious individuals by SARS-CoV-2 in USA.
Figure 20. Number of infectious individuals by SARS-CoV-2 in USA.
Mathematics 09 01221 g020
Figure 21. Number of infectious individuals by SARS-CoV-2 in France.
Figure 21. Number of infectious individuals by SARS-CoV-2 in France.
Mathematics 09 01221 g021
Figure 22. Comparison between theoretical predictions (blue line) and real data (black dots) for USA by assuming that US is still in the first wave of Coronavirus. The values of the parameters are; α ˜ = 0.021 day 1 , β = 0.095 day 2 , K = 25,000,000, and t L = 54 days , respectively.
Figure 22. Comparison between theoretical predictions (blue line) and real data (black dots) for USA by assuming that US is still in the first wave of Coronavirus. The values of the parameters are; α ˜ = 0.021 day 1 , β = 0.095 day 2 , K = 25,000,000, and t L = 54 days , respectively.
Mathematics 09 01221 g022
Figure 23. Comparison between theoretical predictions (blue lines) and real data (black dots) for USA by assuming that US is in the third wave of Coronavirus. The three series of values of the parameters, corresponding to the three waves, are, respectively: ( α ˜ = 0.043 day 1 , β = 0.45 day 2 ); ( α ˜ = 0.021 day 1 , β = 0.008 day 2 ), and ( α ˜ = 0.01 day 1 , β = 0.02 day 2 ). K = 25,000,000 and t L = 54 days .
Figure 23. Comparison between theoretical predictions (blue lines) and real data (black dots) for USA by assuming that US is in the third wave of Coronavirus. The three series of values of the parameters, corresponding to the three waves, are, respectively: ( α ˜ = 0.043 day 1 , β = 0.45 day 2 ); ( α ˜ = 0.021 day 1 , β = 0.008 day 2 ), and ( α ˜ = 0.01 day 1 , β = 0.02 day 2 ). K = 25,000,000 and t L = 54 days .
Mathematics 09 01221 g023
Figure 24. Comparison between theoretical predictions (blue line) and real data (black dots) for France during the first wave of Coronavirus. The values of the parameters are; α ˜ = 0.145 day 1 , β = 0.5 day 2 , K = 110,000, and t L = 53 days , respectively.
Figure 24. Comparison between theoretical predictions (blue line) and real data (black dots) for France during the first wave of Coronavirus. The values of the parameters are; α ˜ = 0.145 day 1 , β = 0.5 day 2 , K = 110,000, and t L = 53 days , respectively.
Mathematics 09 01221 g024
Figure 25. Comparison between the theoretical predictions (blue line) and real data (black dots) for France during the second wave of Coronavirus. The values of the parameters are; α ˜ = 0.035 day 1 , β = 0.5 day 2 , K = 5,000,000, and t L = 278 days , respectively.
Figure 25. Comparison between the theoretical predictions (blue line) and real data (black dots) for France during the second wave of Coronavirus. The values of the parameters are; α ˜ = 0.035 day 1 , β = 0.5 day 2 , K = 5,000,000, and t L = 278 days , respectively.
Mathematics 09 01221 g025
Figure 26. 200 realisations of the stochastic differential equation for France’s first wave. The values of the parameters are those reported in Figure 24. η ˜ = 0.002 , which has been estimated according to Equation (65). The black curve is the solution of the deterministic equation.
Figure 26. 200 realisations of the stochastic differential equation for France’s first wave. The values of the parameters are those reported in Figure 24. η ˜ = 0.002 , which has been estimated according to Equation (65). The black curve is the solution of the deterministic equation.
Mathematics 09 01221 g026
Figure 27. 200 realisations of the stochastic differential equation for France’s second wave. The values of the parameters are those reported in Figure 25. η ˜ = 0.001 , which has been estimated according to Equation (65).
Figure 27. 200 realisations of the stochastic differential equation for France’s second wave. The values of the parameters are those reported in Figure 25. η ˜ = 0.001 , which has been estimated according to Equation (65).
Mathematics 09 01221 g027
Figure 28. Set of paths (200 realisations) of the stochastic differential equation for the USA’s first wave. The values of the parameters are those reported in Figure 22. The intensity of the noise, has been set arbitrarily at the value η ˜ = 0.009 , i.e., 45 times greater than the value estimated by Equation (65).
Figure 28. Set of paths (200 realisations) of the stochastic differential equation for the USA’s first wave. The values of the parameters are those reported in Figure 22. The intensity of the noise, has been set arbitrarily at the value η ˜ = 0.009 , i.e., 45 times greater than the value estimated by Equation (65).
Mathematics 09 01221 g028
Figure 29. Set of paths (200 realisations) of the stochastic differential equation for USA’s three waves model. The values of the parameters are those reported in Figure 23. The intensity of the noise is η ˜ = 0.0002 , which has been obtained by using Equation (65). The black curve is the solution of the deterministic equation.
Figure 29. Set of paths (200 realisations) of the stochastic differential equation for USA’s three waves model. The values of the parameters are those reported in Figure 23. The intensity of the noise is η ˜ = 0.0002 , which has been obtained by using Equation (65). The black curve is the solution of the deterministic equation.
Mathematics 09 01221 g029
Figure 30. USA-Three-Waves Model. US-PDF at the beginning of the lockdown measures (55 days).
Figure 30. USA-Three-Waves Model. US-PDF at the beginning of the lockdown measures (55 days).
Mathematics 09 01221 g030
Figure 31. USA-Three-Waves Model. US-PDF after 90 days.
Figure 31. USA-Three-Waves Model. US-PDF after 90 days.
Mathematics 09 01221 g031
Figure 32. USA-Three-Waves Model. US-PDF after 200 days.
Figure 32. USA-Three-Waves Model. US-PDF after 200 days.
Mathematics 09 01221 g032
Figure 33. USA-Three-Waves Model. US-PDF after 250 days.
Figure 33. USA-Three-Waves Model. US-PDF after 250 days.
Mathematics 09 01221 g033
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sonnino, G.; Mora, F.; Nardone, P. A Stochastic Kinetic Type Reactions Model for COVID-19. Mathematics 2021, 9, 1221. https://doi.org/10.3390/math9111221

AMA Style

Sonnino G, Mora F, Nardone P. A Stochastic Kinetic Type Reactions Model for COVID-19. Mathematics. 2021; 9(11):1221. https://doi.org/10.3390/math9111221

Chicago/Turabian Style

Sonnino, Giorgio, Fernando Mora, and Pasquale Nardone. 2021. "A Stochastic Kinetic Type Reactions Model for COVID-19" Mathematics 9, no. 11: 1221. https://doi.org/10.3390/math9111221

APA Style

Sonnino, G., Mora, F., & Nardone, P. (2021). A Stochastic Kinetic Type Reactions Model for COVID-19. Mathematics, 9(11), 1221. https://doi.org/10.3390/math9111221

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop