key: cord-0851164-xq6g1xb3
authors: Dimarco, G.; Toscani, G.; Zanella, M.
title: Optimal control of epidemic spreading in presence of social heterogeneity
date: 2021-07-26
journal: Philos Trans A Math Phys Eng Sci
DOI: 10.1098/rsta.2021.0160
sha: 37539d5a703332f71d8bd738089e7a37d7c16986
doc_id: 851164
cord_uid: xq6g1xb3

The spread of COVID-19 has been thwarted in most countries through non-pharmaceutical interventions. In particular, the most effective measures in this direction have been the stay-at-home and closure strategies of businesses and schools. However, population-wide lockdowns are far from being optimal carrying heavy economic consequences. Therefore, there is nowadays a strong interest in designing more efficient restrictions. In this work, starting from a recent kinetic-type model which takes into account the heterogeneity described by the social contact of individuals, we analyze the effects of introducing an optimal control strategy into the system, to limit selectively the mean number of contacts and reduce consequently the number of infected cases. Thanks to a data-driven approach, we show that this new mathematical model permits to assess the effects of the social limitations. Finally, using the model introduced here and starting from the available data, we show the effectivity of the proposed selective measures to dampen the epidemic trends.

Non pharmaceutical interventions (NPIs) had a strong impact in limiting the spreading of the SARS-CoV-2 [4] . The introduced NPIs span from basic personal hygiene precautions to more severe actions such as the mandatory stay-at-home and business and schools closures.

These measures were intended to diminish transmission rates and, to this end, aimed at reducing person-to-person contacts via social distancing. In fact, even if the precise role played by those interventions has still to be carefully understood in view of a deeper comprehension of the transmission mechanisms of the disease [5, 9] , it has however been recognized how the social attitudes of a population is a crucial feature for the spreading of airborne diseases [1, 6, 11, 19, 20] . Thus, it appears that understanding the interplay between social heterogeneity and clinical characteristics of a disease is an essential feature to design efficient and robust countermeasures in the near future.

The classical epidemiological models [26] are usually based on a compartmentalization of the population, denoted in the following by C, in which the agents composing the system are split into epidemiologically relevant states, like susceptible, who can contract the disease, and several other classes of agents who may transmit it. Therefore, to model the transmission of a pathogen through a population, if N i (t) represents the number of individuals in the compartment i ∈ C with |C| = n > 1, one generally considers a system of differential equations that can be written in compact form as

where N(t) is the n-dimensional vector with components N i (t), and P is the n-dimensional vector whose components P i (N), i ∈ C, are the transition rates between compartments, namely the variation of the number of individuals in the i-th compartment due to interactions with individuals in other compartments. Several variations of system (1) have been introduced to forecast the penetration of a virus and to prevent its spread, see e.g. [3, 10, 16, 26] . Generalizations of the mentioned compartmental models include additional factors characterizing the transmission dynamics, such as the introduction of age structure, contacts heterogeneity and mobility, see e.g. [1, 2, 14, 15] .

In this work, we propose a kinetic approach to derive consistent transition rates in a population characterized by a heterogeneous number of contacts together with selective NPIs. Indeed, the number of social contacts has been recently identified by the scientific community as one of the main relevant causes of the potential pathogen transmission [6, 19, 20, 30] . In this respect, the NPIs imposed by the different governments were mainly based on restrictions whose scope was to reduce as much as possible the number of contacts among individuals.

In more details, we aim at giving a deeper understanding of the mitigation effects due to the reduction of the social interactions among individuals in the dynamics of SARS-CoV-2. This is done through a combination of a multiscale mathematical model which takes into account both the statistical information on the distribution of social contacts in a society as well as a control strategy whose target is to point the population towards a given benchmark number of contacts. The underlying theoretical framework that we consider is the one coming from the kinetic theory of collective social phenomena [32] . The diffusion of epidemics is here described as a result of the interactions between a large number of individuals each one having a different social propensity. In this respect, as shown in [17] and [38] , the fundamental tools of statistical physics are capable to provide a useful insight on epidemiological dynamics linking classical compartmentalization approaches with a statistical understanding of the social aspects. In fact, the essential multiscale nature of kinetic theory permits to determine the resulting macroscopic (or aggregate) and measurable characteristics of the evolution of the disease [38] . One of the key directions we pursue in this research is given by the combination of consolidated theoretical and modeling tools [32] with data-driven techniques. More precisely, starting from the data at disposal, we calibrate our mathematical model passing through the determination of the relevant epidemiological parameters of the epidemic from the data. Subsequently, we concentrate ourselves on modeling the lockdown policies through an optimal control approach. In particular, we show how different kinds of control may result in very different mitigation effects, that deeply depend on the heterogeneity in the contact distribution of the population. In the last part, starting from the calibrated model at our disposal, we dedicate ourselves to the study of the effects induced by alternative control strategies.

We stress that, even if in the rest of the paper we mainly refer our results to a compartmental model with exposed and asymptomatic, i.e. the so-called SEIAR compartmentalization, the present approach is quite general, and applies both to simple [3, 10, 16, 26] as well as to more complex epidemiological models [22, 23, 29] . For this reason we will keep the presentation as general as possible.

We discuss in this part a possible improvement of the classical compartmental description of epidemic spreading which takes into account the statistical description of the social behavior of individuals. This permits to determine the distribution of the average number of connections in a system of N interacting individuals, see [17, 38] , and to relate this quantity with the infection dynamics. In this new approach, agents which are supposed to belong to a given compartment in term of response to a given pathogen, are additionally characterized by their daily number of contacts x ≥ 0, whose statistical description is obtained through the distribution functions f i (x, t), t ≥ 0, where i ∈ C denotes the i-th compartment (Susceptible, Infected, Asymptomatic, etc..), such that

In this setting, the time evolution of the functions f i (x, t), i ∈ C is obtained by integrating the compartmentalization epidemiological description with the simultaneous characterization of the formation of social contacts in a society [17] . This merging is expressed by the system

As in (1), P is the n-dimensional vector whose components P i (x, f ), i ∈ C, are the transition rates between compartments (now depending both on the contact variable x and on the vector f (x, t)), while Q is the n-dimensional vector whose components Q i (f i ), i ∈ C are suitable relaxation operators which describe the formation of equilibrium distribution of contacts of the underlying compartments. Details about this emerging equilibrium distribution, together with its main properties, will be provided in the following. Last, τ > 0 measures the time-scale at which these equilibrium distributions of social contacts are reached by relaxation. In general, the description of the time evolution of epidemic spreading in terms of the distributions f i (x, t), as given by (2) , is richer than the classic one obeying to (1). This is due to the fact that, in (2), the system depends on an additional variable allowing a more realistic characterization of the transition rates between compartments. In particular, by assuming that the pathogen transmission does not depend upon the contact among individuals, the proposed description (2) collapses to the standard approach (1) as shown in [17] .

In our model, the transition rates can now be assumed dependent on the intensity of social contacts in the compartments. As detailed later, in fact, they are obtained taking into account, in a very natural way, the heterogeneity of the population in the form of the variance of the distributions f i (x, t) [17, 31] . This point of view provides a better understanding of the role of the different individual behaviors in the spreading of a disease [11] .

In general, the knowledge of the densities f i (x, t) allows to evaluate various macroscopic quantities, independent from the number of contacts, the so-called moments of the distributions, obtained by integration with respect to the variable x ∈ R + . In this respect, the zero-order moment gives the number of individuals N i (t) in each compartment at time t ≥ 0

In the same way, one can quantify the mean numbersx i (t) of contacts of individuals in each compartment at time t ≥ 0 evaluatinḡ

Let now f S , S ∈ C, indicate the contact distribution of susceptible people (people at risk of becoming infected), and let the subset I ⊂ C, |I| < |C| include the compartments of people that may transmit the infection. Since the transition rate P S (x, f (x, t)) contains the details about the spreading of the infection among susceptible individuals, the new formulation (2) allows in particular to characterize this transmission in terms of the intensity of social contacts between the compartments of susceptible and infectious individuals. This intensity is measured at time t ≥ 0 by relating the densities f S (x, t) and f j (x, t), j ∈ I. Following [17] , the transition rate P S is then expressed by introducing the function K(f S , {f j } j∈I ), the so-called local incidence rate, given by

In (5) the contact functions κ j (x, y) are nonnegative functions growing with respect to the number of contacts x and y of the populations of susceptible and infectious agents in the compartment j ∈ I, and such that κ j (x, 0) = 0. In particular, the choice

with constants α, β j > 0, corresponds to consider an incidence rate dependent on the product of the number of contacts of susceptible and infectious people. When α = 1, for example, the incidence rate takes the simpler form

We discuss now the details of the dynamics of formation of the social contacts which are modeled by the operators Q i (f i ) in (2) . These operators are of Fokker-Planck-type with variable coefficient of diffusion and linear drift

coupled with no-flux boundary conditions such that the Fokker-Planck-type equation

is mass and momentum preserving. In (7) σ and µ are positive constants which characterize the intensities of the diffusion and drift terms. The above description through Fokker-Planck operators is obtained thanks to an upscaling of an underlying microscopic dynamics which takes into account the typical aspects of human behaviour for what concerns the social attitudes [17] . In particular, the microscopic description takes into account the interactions due to the basic daily needs such as work duties, schools and the search, in absence of epidemics, of opportunities for socialization. These aspects are embedded in the microscopic model by relying on the so-called prospect theory of Kahneman and Tversky [27] as detailed later.

The kernels of the operators Q i furnish the class of equilibrium densities, that are obtained by solving the first-order differential equations

It is immediate to verify that, in this setting, the equilibrium distributions f ∞ i (x) of given mass N i and local meanx i are the Gamma distributions

characterized by the shape parameter λ = µ/σ. It is important to stress the role of the positive constant λ which is a measure of the heterogeneity of the population, statistically distributed according to (8) . It is also remarkable to notice that the above equilibrium dynamics, satisfying the condition λ > 1, is in agreement with the experimental results in [6] . This study represents one of the first attempts in which a large scale experimental measure of the relation between contact and spread of infectious diseases has been done. To better characterize the equilibrium distribution, we also notice that the condition on λ identifies the Gamma densities for which f ∞ i (x = 0) = 0. Moreover, higher order moments of this class of distributions depend only on N i ,x i and λ. In particular, the local variance is expressed by

Hence, expression (9) permits to relate the large values of λ with small heterogeneity of the population as claimed before.

Before introducing in the model the aspects related to the control policies, we briefly discuss how the Fokker-Planck-type operators (7) can be derived from collisional kinetic equations describing the evolution of social connections [24, 36] . These kinetic equations are obtained by looking to a possible characterization of the personal behavior of a single agent in terms of daily contacts starting from some well-established assumptions. Indeed, the total amount of contacts can be simply view as the result of a repeated upgrading: the daily life of each person is based on a certain number of activities each carrying a certain number of contacts. Moreover, independently of the personal choices or needs, this number contains a certain amount of randomness, that has to be taken into account. These assumptions give, letting x be the microscopic number of social contacts resulting from the interaction,

with η a centered random variable such that E[η 2 ] = σ 2 with bounded moments up to order three, and Φ (·) a transition function of the form

The transition function has been constructed by taking into account that while is very easy to reach a high number of social contacts for need or will, going below a given threshold is very difficult, since contacts are forced by basic activities which cannot be avoided in general. This asymmetry between growth and decrease, as exhaustively discussed in [17, 18, 24, 36, 37] , can be suitably modeled by (11) which plays the role of the so-called value function in the prospect theory of Kahneman and Tversky [27] and implements the advocated asymmetry in the contact distribution formation. Hence, with the above discussed choices, at the aggregate scale the evolution of a contact distribution f i, is given by the Boltzmann-type equation in weak form

where B(x) = 1/x is the interaction kernel assigning a low probability to interactions to individuals exhibiting a large number of contacts and assigning a high probability to interactions when x is small. Letting finally → 0 + (the limit of grazing interactions) one shows that (12) converges, up to extraction of a subsequence, to the Fokker-Planck equation (7), see e.g. [17, 21] . This result is particularly useful since it permits to characterize the equilibrium state of the system, i.e. it permits to find the equilibrium Gamma distribution (8) as discussed in [17] and experimentally found in [6] . Moreover, it allows to analytically consider and study control problems at the kinetic level, i.e. to define effective actions at the microscopic level which scope is to modify the behavior of the single individuals. Subsequently, thanks to well established kinetic theory tools [32] , it is possible to upscale the microscopic actions in order to compute its reflection at the aggregate, measurable, scale of description and of interest.

In this part we model the lockdown measures intended to reduce the social interactions by introducing an additive term u at the level of the microscopic formation of the contacts (10) to limit selectively the social activities. This gives the controlled elementary interaction

The function W (x) is such that if W ≡ 1 the control variable u is independent from the number x ≥ 0 of the social contacts, while if W (x) ∝ x the control is considered heavier for people with a higher number of contacts. Note that in (13) , the size of the controlled variable is tuned by the small parameters and τ , which ensure that the control acts at the right scale with respect to both the grazing variable characterizing the passage from a Boltzmann to a Fokker-Planck dynamics and the relaxation variable τ characterizing the speed at which an equilibrium state for what concerns the distribution of contact is reached, compared to the speed at which the epidemic dynamics travels. Then, the so-called optimal control u * is determined as the minimizer of a given cost functional

subject to the constraint (13) . In (14) the minimum is taken on the set U of all admissible controls. A classical choice for J is a quadratic cost functional of the form

being ν > 0 a penalization coefficient measuring the economic and social cost of the control and x T,i > 0 the desired target number of social contacts for the compartment i ∈ C. According to (15) , the cost of the control increases quadratically with the distance to the desired target of average number of contacts. The presence in (15) of the mean operator is related to the fact that, by hypothesis, a control cannot act on the random fluctuations introduced to mimic the unanticipated daily interactions among agents. Let us observe that other convex cost functionals may be considered as well. However, they often lead to problems whose analytical solution cannot be obtained explicitly. Proceeding now as in [34] the control introduced in (13) can be solved explicitly giving

which generates the optimal constrained law of growth of the social contacts

Now, in presence of the above described microscopic controlled interactions, one can consider as before the construction of a kinetic collisional equation for what concerns the description of the formation of the social structure

where we recall that the frequency of interaction has been chosen such that B(x) = 1/x. Then, thanks to the passage to the grazing limit one can show that the controlled kinetic model converges to a Fokker-Planck type equation which contains an additional drift term, that quantifies the role of the control mimicking lockdown measures. In fact, since by hypothesis 1, we can write

withx ∈ (min(x, x * ), max(x, x * )). Then using (17) one gets

where the remainder R Φ (f i, )(x, t) can be shown to go to zero in the limit → 0 + since by assumption ϕ and its derivatives are bounded in R + and the stochastic variable η is assumed to have bounded moments at least of order 3 as previously stated. We point the reader to [17, 32] for a detailed discussion on this aspect. In the limit → 0 + one gets

since for 1 we have

Integrating back by part (20) , and considering the no flux boundary condition at x = 0 + to preserve mass and momentum, the relaxation equations take the form

where the operator C i (f i ) in (21) is given by

and Q i (f i (x, t)) remains the same as in (7). Note in particular that in (22) the control term does not depend on the time scale τ of the Fokker-Planck operator. Therefore, in presence of social restrictions modeled by the control term introduced above, system (2) is modified as

where the n-dimensional vector C quantifies the effect of the control on the epidemic spreading. According to the previous analysis, the components of the vector control C are given by the operators of type (22) , where the values of the penalization coefficients and of the target values may depend on the specific compartment i ∈ C.

In the next part, we focus on a specific compartmental model for which we give the details of the controlled strategy for the reduction of the social contacts. Remark 2.1. We observe that the cost functional (15) introduced at the level of microscopic interactions can be seen as an instantaneous approximation of a macroscopic functional obtained by considering (16) in the limit , τ → 0 + whose form is

In the above functional we consider the contributions at the final time T and the steady distribution f ∞ i of contacts defined in (8) instead of the instantaneous single interaction among agents.

We consider the case C = {S, E, I, A, R}, i.e. we concentrate on a model with exposed and asymptomatic compartments and we discuss the observable effects of the introduced control at the level of the mean number of connections in two leading examples.

The choice C = {S, E, I, A, R} gives, accordingly to (23), the following system of kinetic equations relating the epidemiological dynamics with the number of daily contacts expressed in terms of a probability distributions f i (x, t)

where, to have a completely definite system, it remains to clarify the role of the parameter ξ which gives the fraction of exposed individuals becoming successively infected and its counterpart (1 − ξ) which gives the fraction of asymptomatic individuals after the exposure to the virus. In addition, the functions γ I and γ A represent the recovery rates of infected and asymptomatic individuals respectively, which give the frequency at which individuals move to the compartment of recovered individuals. In the present setting, the choice of the constants β I , β A > 0 takes into account that the population of the infected compartment is considered recognized by the health system through a swab test and therefore subject to isolation. We collect in the asymptomatic compartment unrecognized infectious that are not subject to isolation. We highlight that we preferred to neglect the compartment of deceased people since they do not influence the spreading of the disease. Anyway, the inclusion of this new compartment to assess the risks related to the infection can easily be considered. Since the Q i and C i components in (24) , coupled with no-flux boundary conditions, are mass preserving, we can integrate the equations of the system with respect to the contact variable x to obtain a macroscopic system of equations describing the evolution of the fraction of individuals in each compartment independently on the number of daily social contacts. This gives the system

which is not closed since the evolution in time of the vector solution N(t) depends on the mean number of contactsx i . Since the Q i components, under the same hypothesis on the boundary fluxes, are also momentum preserving, by multiplying the distributions f i (x, t) by x and integrating with respect to the contact variable, it is possible to recover a further macroscopic set of equations describing the constrained dynamics of the mean number of social contacts in each compartment. However, one immediately realizes that, this second set of macroscopic observables depends upon higher order moments of the distributions f i (x, t).

In other words, in general, the system of macroscopic equations which can be obtained by integration of (24) over the number of social contacts is not closed. In fact, system (25) depends on the time evolution ofx i (t), i = S, E, I, A, R while the system giving the time evolution ofx i (t), i = S, E, I, A, R depend upon second order moments of the f i (x, t), i = S, E, I, A, R. This represents nothing else than the classical so-called closure problem often encountered in kinetic theory [32] . However, letting τ → 0 forces the distributions f i (x, t) to converge towards the Gamma vector distribution of components (8) with mass fractions N i (t), and local mean valuesx i (t), i ∈ C = {S, E, I, A, R}. This amounts to assuming that social dynamics are much faster than epidemic dynamics. In this situation, it is possible to obtain a closed set of macroscopic equations. In particular, the second order moments of the distributions f i (x, t) can be explicitly computed in terms of the mass fractions and mean values resorting to (9) . This gives the system

where

Equations (26) coupled with the system (25) describe in closed form the time evolution of an epidemic where the transition functions from one compartment to another depend upon the mean level of social contacts in the population. In addition, these latter depend on the heterogeneity of the population through the parameter λ = µ/σ characterizing the distributions f i (x, t), while the term (27) models the effects of the lockdown measures.

In the rest of the paper, we focus on two specific benchmark cases for which it is possible to compute explicitly the effect of the social restrictions on the system. In particular, we first concentrate on the case of uniform restrictions, corresponding to W (x) ≡ 1, and second on the case of selective restrictions. This last situation describes an agent system in which restrictions are stronger for individuals exhibiting a larger number of contacts, and it is based on the choice W (x) = x/x i . By direct computation we have

Hence, for small penalizations ν > 0 of the control, the mean number of connections stabilizes towards the values

It is worth to notice that, in both situations, the role of the penalization parameter ν > 0 is limited to a scaling of time, without changing the final goal of the control. Thus, a smaller penalization parameter is reflected in a quicker convergence of the solution towards the target value as expected. We also notice that the selective control, i.e. the one especially operating on individuals with larger social nets, allows to reach for the global population a value below the desired target x T . In particular, the mean contact rate is lower than x T of an amount directly proportional to the degree of heterogeneity of the population. This can be explained by the fact that individuals with larger social networks are the main responsible for the creation of the contact dynamics. This situation suggests that it is important to act on this class of people in order to limit the spread of a disease instead of taking measures with uniform impact over the whole population.

From an operative point of view, it is clear that the application of a uniform control, which is based on restrictive measures uniformly valid for the global population, are simpler to apply with respect to measures which provide for heavier restrictions for the categories that show a greater number of social contacts. In this direction, the previous analysis seems to suggest that a mixture of the two types of controls will present the best compromise between containment results and effective applicability of containment measures.

In this section using (25)-(26), we compare possible interventions for the containment of an epidemic based on the limitation of the mean number of contacts. These non-pharmaceutical measures are designed to drive the population towards a prescribed target dependent on the relevant epidemiological state of agents, namely x T,i , i ∈ C in the considered framework.

To this end, we first estimate the relevant epidemiological parameters of the second order macroscopic SEIAR model taking into account the available dataset provided by John Hopkins University. 1 After the calibration of the model with the data at disposal, we assess the impact of different control strategies leading to the observed epidemic dynamics in different countries where lockdown policies have been implemented.

We will mainly focus on the so-called first wave of infection lasting the first half of 2020 since in this period similar contact restrictions have been proposed in all Western countries. In the following we concentrate on the cases of Italy and United Kingdom.

We consider two countries that adopted the stay-at-home policies during early stages of the infection. In Table 1 we report the days of implementation of such measures together with a sketch of the timeline of the epidemic. In the second column, we highlight the begin of restrictions at the global level for each different country. Anyway, it is worth to mention that often local restriction measures, i.e. isolation of a small portion of a territory, have been implemented as well.

Before proceeding with the estimation of the model parameters some remarks are deserved. First, we stress that the calibration of epidemiological models is generally difficult and require special cares at the numerical level. Furthermore, available data typically represent a lower bound with respect to the true numbers of an infection, especially in its initial phases, due to a limited testing capacity. Here, we do not consider this aspect and we refer to recent attempts to overcome this issue based on tools of uncertainty quantification [1, 38] that can be successfully employed in this context. For related approaches we mention also [7, 8, 13, 35] .

For the calibration of the model (25)- (26) we have adopted a bilevel approach . In particular, in the phase preceding the lockdown, we estimated the epidemic parameters in the unconstrained regime assuming that no restriction of social contacts was activated, i.e. C ≡ 0. We also fixed the known clinical parameters according with values reported in the literature: see [12, 22] . We set then γ E = 1/3.32, γ I = 1/10 and γ A = 2γ I , β A = 2β I and we also assumed that ξ = 0.2, such that the 80% of infections come from asymptomatic or paucisymptomatic individuals, in agreement with the recent serological campaigns promoted in Italy 2 . The proposed choice for β A , β I > 0 is coherent with existing literature, see e.g. [33] .

We then solved a least square problem based on the minimization of the relative L 2 norm of the difference between reported number of infectedN I and recoveredN R , and the theoretical evolution of the model N I (t) and N R (t) for t ∈ [t 0 , t ] with N E (t 0 ) = N I (t 0 ) = N A (t 0 ) = N R (t 0 ) = 1, i.e. we assume that in t 0 the pandemic was nearly started. We also considered x S (t 0 ) =x E (t 0 ) =x A (t 0 ) =x R (t 0 ) = 10, in agreement with the experimentally observed mean number of contacts in a Western country before the pandemic [6] , andx I (t) ≡ 3 for recognized infected, corresponding to the average number of family contacts. Thus once the known parameters have been fixed, we considered the following minimization problem

where θ ∈ [0, 1], · L 2 ([t 0 ,t ]) the relative norm over the time horizon [t 0 , t ] and under the constraints β I ∈ [0, 1/100], λ ∈ [0, 50]. In third and fourth columns of Table 1 we report the estimated epidemiological parameters obtained for θ = 10 −3 to get a better resolution of the infected dynamics. It is worth to notice that, according to these values, Italy manifests higher heterogeneity linked to the contact distribution compared to UK.

Once the relevant epidemiological parameters are estimated, we focus on the subsequent lockdown phase where we look for the optimal value x T in the control term which permits to fit the data. We assume this value x T equal for each compartment and we study the two cases of uniform and selective restrictions, respectively obtained by fixing W (x) = 1 and W (x) = x/x i and C(f ∞ ) of the form (28) in the dynamics. Hence, we solve an optimization problem in the lockdown time interval [t + 1, t f ], for a sequence of time steps t n over a time window of one week

In Figure 1 we show the obtained trends for the registered cases for Italy and United Kingdom. We may observe how for both controls the macroscopic trends are coherent with the observed ones. Therefore, the adopted fitting procedure allows us to compare the estimated target x T,i for mean contacts in the considered control strategies. In particular, we observe how the uniform control implies a stronger contact reduction with respect to the selective one, while a similar cost, computed as the sum J S + J E + J A + J R , is produced by the two containment strategies. In other words, the adoption of selective-type containment strategies based on the average number of contacts allows to mitigate the severity of lockdown policies.

In the previous part, we focused on getting the target mean number of contacts for both the control strategies. We remark that the NPIs implemented by both countries have been mainly uniform in the first wave of the epidemic, meaning that, in the present model framework, the average number of contacts where approximately the ones defined by the red line with circle markers in the second row of Figure 1 . Hence, it is interesting to consider a retrospective analysis where the estimated x T,i , i ∈ C, is instead implemented in the dynamics with selective control. In Figure 2 , we present the results of such study, where the evolution of the disease in presence of selective control with target given by the uniform control is reported. It is evident how the amplitude of the epidemic peak is strongly reduced and the fraction of asymptomatic cases is under control in both cases using this alternative strategy. Furthermore, it is worth to observe that a selective-type control is more effective in reducing the peak of infection in presence of high heterogeneity like in the case of Italy.

The recent spreading of COVID-19 epidemic has been thwarted in most countries through non-pharmaceutical interventions, intended to diminish transmission rates and, to this end, reducing person-to-person contacts via social distancing. Even if the confinement measures appear to be a correct tool to limit the infection spreading, a precise quantification of the effects played by them has still to be carefully understood.

Aiming to investigate the impact of the number of daily contacts in the spread of infectious diseases, in [17] we introduced a mathematical description of the epidemic spreading by integrating the epidemiological dynamics of the classical SIR-type compartmental model with a kinetic term quantifying the population-based contacts.

The kinetic description of the formation of social contacts considered in [17] , based on elementary interactions, turns out to be very flexible with respect to external inputs, like a control forcing the number of daily social contacts towards a fixed target. From the modeling point of view, this strategy leads to integrate the compartmental models with social contacts with a new term quantifying in a precise way the effect of the control strategy in dependence of the population heterogeneity.

The main conclusion of our theoretical analysis, confirmed by numerical experiments, is that a selective control, when compared to a uniform control, and at the same social cost, leads to a marked reduction of the epidemic spreading. This reduction is further linked to the parameter of social heterogeneity, so that the selective control gives a better performance in a population with higher heterogeneity.

The main conclusion of our theoretical analysis, confirmed by numerical experiments, is that a selective control, when compared to a uniform control, and at the same social cost, leads to a marked reduction of the epidemic spreading. This reduction is further linked to the parameter of social heterogeneity, so that the selective control gives a better performance in a population with higher heterogeneity.

Control with uncertain data of socially structured compartmental models

Modeling the spatiotemporal epidemic spreading of COVID-19 and the impact of mobility and social distancing interventions

Infectious Diseases of Humans: Dynamics and Control

Estimating the effects of non-pharmaceutical interventions on the number of new infections with COVID-19 during the first epidemic wave

Assessing mandatory stay-athome and business closure effects on the spread of COVI-19

The French Connection: The First Large Population-Based Contact Survey in France Relevant for the Spread of Infectious Diseases

Hyperbolic compartmental models for epidemic spread on networks with uncertain data: application to the emergence of Covid-19 in Italy

Spatial spread of COVID-19 epidemic outbreak in Italy using multiscale kinetic transport equations with uncertainty

Economic and social consequences of human mobility restrictions under COVID-19

A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2

Effects of information-induced behavioural changes during the COVD-19 lockdowns: the case of Italy

Parameter estimation and uncertainty quantification for an epidemic model

Epidemic modeling in metapopulation systems with heterogeneous coupling pattern: Theory and simulations

Infectious disease dynamics in metapopulations with heterogeneous transmission and recurrent mobility

Mathematical epidemiology of infectious diseases: model building, analysis and interpretation

Kinetic models for epidemic dynamics with social heterogeneity

Kinetic modeling of alcohol consumption

Heterogeneous social interactions and the COVID-19 lockdown outcome in a multi-group SEIR model

Inferring the structure of social contacts from demographic data in the analysis of infectious diseases spread

Fokker-Planck equations in the modelling of socio-economic phenomena

Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures

Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy

Human behavior and lognormal distribution. A kinetic description

A new view of multiscale stochastic impulsive systems for modeling and control of epidemics

The Mathematics of Infectious Diseases

Choices, values, and frames

Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vó

A viral load-based model for epidemic spread on spatial networks

Social contacts and mixing patterns relevant to the spread of infectious diseases

On the spread of epidemics in a closed heterogeneous population

Interacting multiagent systems: kinetic equations and Monte Carlo methods

SUIHTER: a new mathematical model for COVID-19. Application to the analysis of the second epidemic outbreak in Italy

Control of tumour growth distributions through kinetic methods

Epidemic models with uncertainty in the reproduction number

Statistical description of human addiction phenomena

Entropy-type inequalities for generalized Gamma densities

A data-driven epidemic model with social structure for understanding the COVID-19 infection on a heavily affected Italian Province

This work has been written within the activities of GNFM and GNCS groups of INdAM (National Institute of High Mathematics). The research was partially supported by the Italian Ministry of Education, University and Research (MIUR: Project "Optimal mass transportation, geometrical and functional inequalities with applications", and Dipartimenti di Eccellenza Program (2018-2022) -Dept. of Mathematics "F. Casorati", University of Pavia. G.D. would like to thank the Italian Ministry of Instruction, University and Research (MIUR) to support this research with funds coming from PRIN Project 2017 (No. 2017KKJP4X entitled "Innovative numerical methods for evolutionary partial differential equations and applications").