A Stochastic Kinetic Type Reactions Model for COVID-19

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 SIS-model (Susceptible ⇒ Infectious ⇒ Susceptible) and the stochastic SIS-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.


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 • 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-andtrue 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 SIS 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 SIS-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 abovementioned question for the SIS-model subjected to lockdown measures. In Section 3, we introduce, and study, the Stochastic SIS-model subjected to lockdown measures (the (SIS) 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 (SISI h r h d h ) Lmodel. 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 (SISI 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 A and B.

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: (1) 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 (2) 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: with N Tot. denoting the total population and σ ≡ µN Tot. . By assuming that the dynamics of R is much faster that those of S and I, we may set dR/dt 0 and system (3) reduces to which corresponds to the model In the relevant literature, model (5) is referred to as the SIS-model (see, for example, [15]). Notice that condition dR/dt 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 Hence, the dynamics of infectious people is governed by the logistic model 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 SIS-model the lockdowninduced decrease rate c(t) 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 = const.). 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: Equation (9) may be referred to as the (SIS) L -model where L stands for Lockdown. 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 witht denoting the normalised timet ≡ t/t L and β ij are real numbers. Note that, the above conditions (2) and (3) impose strict restrictions on the coefficients β ij . 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 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 2 L /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 with β ≡ β 01 = const. 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 with α ≡αt L . As shown in the Appendix A, the exact solution of Equation (15) reads where σ ≡ 2/β. Equation (16) can be rewritten into the form 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 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.

The Second Wave of COVID-19 Pandemic in the (SIS) 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., dI/dt > 0 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 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) 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 timet 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Ĩ, determine the minimum value that I 1 should have so that the number of infectious people is less than I, during an entire period fromt 1 tot f ". We note that within the intervalt 1 ≤t ≤t f the number of infectious people grows quite slowly. Let us decide to interrupt the lockdown measures at timet 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 allt ∈ ∆t, with ∆t <t f −t 1 . It is easily checked that I 1 should satisfy the condition 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Ĩ, withĨ < K/2. The result is easily obtained by writing the initial condition as and by determining the range of values of x such that the inequality is satisfied. After a little algebra, we get To summarise, if the lockdown measures are interrupted at timet 1 , during the successive period of time ∆t =t 2 −t 1 the number of infectious remains inferior to a preset numberĨ 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Ĩ and ∆t such that 50,000 ≤Ĩ ≤ K IT . As shown in Figure 4, if we set the valueĨ = 50,000, the maximum value of I 1 is I 1 9, 500 for t t L , and I 1 700 for t 2t L , respectively. . As shown in the figure, setting the value ofĨ = 50,000, the maximum value of I 1 is I 1 = 9, 500 for ∆t = t L , and I 1 = 700 for ∆t = 2t L , respectively.

The Stochastic (SIS) L -Model
The Stochastic version of Equation (15) reads with ξ(t) denoting a white noise and √ η is the intensity of the noise. δ denotes the Dirac delta function (distribution). Let us now consider the system at the reference state I RS , which is the solution of the deterministic equation subjected to a perturbation of small amplitude δI, i.e., 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 ): the following equations are valid: 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 Solution I RS (t) is given by Equation (17) and the exact solution for δI(t) is easily obtained: Solution (33) allows to check the validity of Equation (31). Indeed, the correlation function δI(t)-noise at the same time reads 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 The solution of Equation (35) provides the expression for the second moment 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 > (37) 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)     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 timet reads [25] ∂ ∂t p( with µ denoting the drift coefficient and D the diffusion coefficient, respectively. Equation (38) must satisfies the initial condition 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).

The (SISI h r h d h ) L -Model
As can be seen, according to the SIS-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: 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: 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 Of course, we have where C h = Total hospital s capacity (43) The dynamical equations for the entire process are then: 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 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 Tot. = const. 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 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 In absence of the lockdown measures (β = 0), we have the following scenarios: In words: • For case (i), there will be a proper epidemic outbreak with an increase of the number of the infectious people; • For case (ii), 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.  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 (SISI 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.

The Stochastic (SISI h ) L Model
If the dynamic is subjected to white noise and related stochastic equations are written as with δ ij 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 Let us choose, for example, I hRS = C h , then I RS (t) = I 0RS exp((1−α 2 /β) 2 /σ)t β exp(−(t−α 2 /β) 2 /σ) 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 with 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 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 RS , I hRS ). We get The correlation functions < δI(t 0 )δI(t) > and < δI h (t 0 )δI h (t) > are easily obtained with 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 where Equation (31) have been taken into account. The solutions of system (61) read 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 hRS = 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.t Correlation Functions These three correlation functions are obtained by solving numerically system (61).

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 whereα ≡ α/t L and the white noiseξ(t) is defined by the relations The order of magnitude for the intensity of the noiseη may be derived by using Equation (A11) reported in Appendix B:η 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 Figures 20 and 21, respectively. Figures 22 and 23 show a comparison between the theoretical predictions (blue line) and real data (black dots) for USA. Figure 22  Now, let us deal with the stochastic processes. Figures 26 and 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.
Figures 28 and 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. Figures 30-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.

Conclusions
Waiting for the mass production of vaccines, we studied two scenarios by means of two mathematical models: the (SIS) L -model and the (SISI 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 (SIS) L -model, refers to a situation in which the role of hospitals is not taken into account. The second model (the (SISI 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 SIS 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 (SISI 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 = deceasedpeople, 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.
where c 0 is a constant independent of time. Finally, we get I(t) = exp(αt − C(t)) c 0 + να 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 (µ ch. ), 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 Country = T Env. ") stands for environment equilibrium; (ii) Identical chemical potentials (µ ch. Country = µ ch. Env. ) 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 −1/2 Country [37], with K Country 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 −1/2 Country (A11) 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.