Issue 
Math. Model. Nat. Phenom.
Volume 15, 2020
Coronavirus: Scientific insights and societal aspects



Article Number  57  
Number of page(s)  42  
DOI  https://doi.org/10.1051/mmnp/2020045  
Published online  30 November 2020 
COVID19 pandemic control: balancing detection policy and lockdown intervention under ICU sustainability^{*}
^{1}
Faculté des Sciences, Université du Québec à Montréal (UQAM),
Montreal,
Québec, Canada.
^{2}
LAMA, Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS,
77454
MarnelaVallée, France.
^{3}
ORFE Department, Princeton University,
Princeton, USA.
^{**} Corresponding author: chi.tran@upem.fr
Received:
21
May
2020
Accepted:
3
November
2020
An extended SIR model, including several features of the recent COVID19 outbreak, is considered: the infected and recovered individuals can either be detected or undetected and we also integrate an intensive care unit (ICU) capacity. We identify the optimal policy for controlling the epidemic dynamics using both lockdown and detection intervention levers, and taking into account the tradeoff between the sanitary and the socioeconomic cost of the pandemic, together with the limited capacity level of ICU. With parametric specification based on the COVID19 literature, we investigate the sensitivities of various quantities on the optimal strategies. The optimal lockdown policy is structured into 4 phases: First a quick and strong lockdown intervention to stop the exponential growth of the contagion; second a short transition to reduce the prevalence of the virus; third a long period with full ICU capacity and stable virus prevalence; finally a return to normal social interactions with disappearance of the virus. The optimal scenario avoids the second wave of infection, provided the lockdown is released sufficiently slowly. Whenever massive resources are introduced to detect infected individuals, the pressure on social distancing can be released, whereas the impact of detection of immune individuals reveals to be more moderate.
Mathematics Subject Classification: 49N90 / 92D30 / 34H05
Key words: Optimal control / lockdown / testing / intensive care units / SARSCoV2
The authors thank Hélène Guérin, Olivier l’Haridon, Olivier Pietquin and Gabriel Turinici for stimulating discussions and feedbacks on early versions of that work, as well as participants of the Machine Learning Paris seminar, of the European Network for Business and Industrial Statistics (ENBIS) seminar, as well as participants of the working group at CMAP, École Polytechnique. R.E. and V.C.T. are partly supported by the Bézout Labex (ANR10LABX58). V.C.T. is also part of ANR Cadence (ANR16CE320007), by the Chair “Modélisation Mathématique et Biodiversité” of Veolia EnvironnementEcole PolytechniqueMuseum National d’Histoire NaturelleFondation X. A.C. gratefully acknowledges funding from NSERC grant 07077.
© The authors. Published by EDP Sciences, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Within a few weeks, the COVID19 took most of the world by surprise, starting with Asia, then Europe and finally America (see [53, 58, 82]). More specifically, the COVID19 turned out to be much more contagious than expected, and official health authorities understood that letting the virus spread freely would probably saturate hospital capacities. In most countries, aggressive quarantine and lockdown measures have been taken, hopefully to stop the spread of the virus and infection, or at least to slow down its diffusion and delay the cases over time to avoid saturation of healthcare systems. The effectiveness of these nonpharmaceutical interventions have been questioned in [41, 52]. It was already suggested in [67], but then emphasized by [32], that these measures should rather aim at ‘flattening the curve’: mitigation seems more realistic than suppression, since eliminating the threat necessarily takes time, and a massive arrival of people in need in hospitals would be dramatic.
The first control used by health authorities was based on quarantines, historically a natural technique to lower social interactions and reduce the spread of the disease, as recalled in [76]. Classically, infected people are isolated, to avoid further transmission, but COVID19 has a potentially high rate of asymptomatic infected individuals, as mentioned in [5, 22, 40, 60] or [84]. This explains why most countries used indiscriminate lockdown. In order to distinguish between the susceptible, but still unaffected, individuals and the either asymptomatic infected or undetected recovered ones, testing is necessary. Piloting these interventions is a key problem, in particular for monitoring exit strategies.
In the recent literature on the COVID19, many papers have focused on the effect of lockdown on the spread of the disease (e.g. [12, 69]) or on the evaluation via simulations of exitstrategies (e.g. [26, 31]). In this paper, we suggest a different approach for evaluating lockdown and exit scenarios based on Optimal Control theory (e.g. [66]). Although the analysis of public health policies using controls is of clear practical value, there has not been a so dense literature on the subject (e.g. [1, 2, 4, 11, 36, 39, 75]). For example, [2] considered the case where it is possible to isolate infectious people (with pure information, in the sense that we can always distinguish infectious, susceptible and recovered), and studies the optimal share of the population of infected people that should be isolated to control the spread of the disease. More recently, [11] considered vaccination, quarantine, screening or health promotion campaigns are possible control in a SIR model, and studied uniqueness and monotonicity properties of mathematical epidemiological models. As said on [75], optimal control theory has proven to be a successful tool in understanding ways to curtail the spread of infectious diseases by devising the optimal diseases intervention strategies. But standard SIR models should be adapted to take into account specific properties of the current pandemic. We consider a modified SIR model, incorporating all the features related to the recent COVID19 outbreak mentioned above, and use optimization algorithms to compute the ‘best’ level of lockdown and detection effort in testing individuals. The criterion to rank the different possible strategies takes into account: (1) the death toll of the epidemic, (2) the economic and social costs of the lockdown, (3) the required effort for detecting infected or immune individuals. Besides, we only consider strategies never exceeding the capacity level U_{max} of the ICU sanitary system, chosen to be dedicated to COVID19 patients.
One of the originality of the model considered here is the particular focus on the outcomes of the intensive care units (ICU) saturation, which in many countries have been put under great pressure and whose capacities are determining factors for the death toll of the disease. We assume a capacity constraint on how many individuals can be treated in ICU at a given time. Once the threshold is reached, the fatality rate for the additional patients sent to ICU rises sharply. For the COVID19 outbreak, we are aware of some works by [7, 25, 65], but none of these works models the ICU capacities. In [3], no ICU capacity is explicitly introduced, but these authors choose a death rate in hospitals that depends linearly on the number of patients. This constraint ishowever at the heart of debates, and in the simulations for evaluating the various possible public health responses (see e.g. [26] or [31]), the ICU capacity is always represented to check whether the different scenarios will be feasible or not. Adding a constraint in the optimization program is thus natural. The nonlinearity induced by the threshold that we introduce provides very interesting phenomena and insights for the control of the disease. Furthermore, we investigate cases where the ICU capacity is increasing and evolving dynamically over time.
A second feature in our work is that we consider two levers for control: the quarantine ratio and the detection effort for infectious and immuneindividuals. One of our motivation is to study qualitatively the balance existing between efficient testing and strong confinement strategies. Actually, two testing procedures can be performed here. A short range one, that can detect, at time t if someone is infected, or not (Polymerase chain reaction or PCR tests, that can detect whether an individual is infected and send her/him in quarantine).And a long range one, that can detect if someone is immune or not (a blood test as in [44]). For this purpose, we introduced a socalled SIDUHR^{+∕−} propagation model for the epidemic, which explicitly allows to take into account the detection ratios of immune and infected individuals.
Our main findings consist in identifying four main phases in the optimally controlled epidemic dynamics. First, a strong and quick lockdown intervention allows to stop the exponential growth of the epidemic spreading and bring the effective reproduction number ℜ_{t} under level 1. Then begins a short transition period during which ℜ_{t} remains below 1 as the prevalence of the virus within the population diminishes progressively until reaching a threshold I^{*}. Since the beginning of the epidemic, the number of individuals requiring ICU is rising and just reaches U_{max} when the virus prevalence attains I^{*}. Then begins a very long period during which the reproduction number ℜ_{t} stabilizes at level 1 while the virus prevalence and the ICU occupation remain stable at levels I^{*} and U_{max}. Finally, when a sufficient proportion of the population has already been infected and the herd immunity level is close enough, a light final release of lockdown measure allows to retrieve a regular level of social interactions and end the epidemic phase. Such a control strategy allows to achieve an epidemic size very close to the herd immunity level, which corresponds to the minimal size of an epidemic, in comparison to a benchmark case without intervention for example. Besides, given the current level of information and confidence on the COVID19 dynamics modelling, the design of robust optimal control strategy should be seen as a minimal requirement. For our model, the enhanced structure of the optimally controlled epidemic dynamics reveals to be very robust to variations of the model specifications within a reasonable range.
The effect of the detection effort within the population reveals to be more subtle. As explained above, two types of detections are considered. Detecting quickly infected individuals and breaking the contamination chains reveals to be key in reducing the economic and sanitary burden of the epidemic. First, it allows to slow down the exponential growth of the epidemic and decrease hereby the required epidemic size for reaching the herd immunity and the induced death toll. It also reduces the lockdown effort on the population, which helps minimize the socioeconomic impact. Optimizing over the detection rate of infected individuals reveals the huge impact on both sanitary and economic burdens of early effort in detection, as well as a stable threshold level in detection. Using the detection of infectious individuals as a control lever allows to keep the reproduction number below 1 while limiting the lockdown effort and the induced socioeconomic cost. Concerning the detection of immune individuals, it allows to reduce significantly the socioeconomic burden, later in time during the epidemic, whereas its impact on the sanitary outcome seems rather limited.
The rest of the paper is organized as follows. In Section 2, we derive our SIDUHR^{+∕−} model from a standard SIR model (see Sect. 2.1), and explain all the components of the dynamics. The infectious individuals may either: (i) be asymptomatic or have mild symptoms that do not require hospitalization, (ii) need hospitalization with or without intensive care. In Section 2.2 we present the optimal control problem and objective function of public health authorities. Possible control levers are discussed in Section 2.2.1, with lockdown, detection efforts, and ICU capacities limitation. In Section 2.2.2 we motivate the choice of the objective function, and discuss optimal control in Sections 2.2.3 and 2.2.4.
In Section 3, we compute the optimal strategy based on a choice of parameters from the COVID19 epidemic in France, but of course the methods and findings can be transposed to any other country or set of parameters. In Section 3.1, the benchmark scenario is described. Then, in Section 3.2 we consider optimal lockdown, and discuss the 4 stage dynamics. In Section 3.3, we provide a sensitivity analysis, and discuss the robustness of our findings in response to changes in the model specifications and parameters. In Sections 3.4 and 3.5, we consider the impact of two other controls: the ICU sustainability level and detection resources. Given the simplicity of the compartmental model considered here, our purpose is not to give detailed prediction of the epidemic but rather to highlight very interesting aspects of the interplay between different levers (lockdown, detection and ICU capacity), provide qualitative description of the outcomes in each scenarios and give rough estimates of the duration and size of the epidemics.
Let us finally emphasize that the level on uncertainty upon the characteristics of the epidemic dynamics is still currently very high. The enhanced conclusions of this paper remain valid within the scope of the chosen model specifications chosen here, which are in line with the current understanding of the epidemic dynamics.
2 Methods
2.1 The SIDUHR^{+∕−} model and its dynamics
Compartmental models, where the population is divided into classes defined by the status of the individuals with respect to the disease or how they are taken in charge by the healthcare system are very popular in mathematical epidemiology (see e.g. [8, 9, 24] for an introduction). One of the simplest and most fundamental of all epidemiological models is the popular SIR model, developed in [47]. The population, of fixed size, is divided into three categories: susceptible (S), infected (I) and recovered (R) individuals. We denote by S_{t}, I_{t} and R_{t} the respective proportions in the populations at time t (measured in days). In order to adapt the SIR model to the COVID19 epidemic and incorporate lockdown and testing strategies with their economic and social impacts, we modify some rates and add compartments to the SIR model.
We want to take into account two important features of the COVID19 pandemic. The first one is the important proportion of asymptomatic patients, and the small use of testing in several countries. We will split the I group in two: I^{−} the nondetected, mostly asymptomatic individuals, but contagious, and I^{+} the detected part of infected individuals. The second feature is the problem highlighted in the ‘flatten the curve’ concept (discussed in [50, 68] or [37]) which means that a goal in the choice of a public health intervention should be to avoid a surge of demand on the health care system and in particular in ICUs. From a modelling perspective, we introduce the proportions H_{t} and U_{t} of individuals using the health care system with or without the need of ICUs. There is an exogenous sustainable limit U_{max} for the ICU system, above which the death rate is significantly higher.
2.1.1 Compartments
More specifically, all possible states can be visualized on Figure 1. The parameters for the dynamics are described in the next section.

S: susceptible, never tested positive, never infected and not dead (where people are by default). They can get infected by infectious and not quarantined individuals, i.e. individuals in I^{−}. When susceptible individuals get infected, they move to I^{−}.

I^{−}: infected nondetected. Those individuals can be asymptomatic, not sick enough to go the hospital, and are not detected. They can either get tested (with type1 tests, and then move to I^{+}), get sicker and then go to the hospital (H) or simply recover and get immune, but still not detected (R^{−})

I^{+}: infected detected (and nonhospitalized). They can either get sicker and then go to the hospital (H) or recover, in which case they end up in the compartment R^{+} since they were tested positive. Here, we assume that the individuals in I^{+} cannot infect anyone since they are strongly isolated.

R^{−}: recovered nondetected. Those can be detected using bloodtype tests (also named type2) and move to R^{+}.

R^{+}: recovered detected. We assume here that all recovered individuals (both R^{−} and R^{+}) are immune long enough with respect to the epidemic duration.

H: hospitalized (and detected). All people entering the hospital get tested. They can either get sicker and then go to intensive care (U), or get immunity after recovering, and then move to R^{+}. Note that H is an intermediary state, before getting really sick, or die. Individuals are not allowed here to go from I (detected or not) to D. Furthermore, we did not distinguish explicitly between mild and severe infected (as in many other models) in order to emphasize the role of detection resources, but, somehow, H could correspond to severe cases.

U: hospitalized in ICU, for people at a more advanced stage (and detected), who need ventilators and other specialized medical care. They arrive from H, and they can either get immunity after recovering, and then move to R^{+}, or die (D).

D: dead, from the disease.
Here, recovered are not contagious anymore and we assume that they get immunity (and therefore, can not get infected again). This is a rather strong assumption, especially on a long term horizon. Reference [16] used a longitudinal study to show that those infected by MERS in South Korea (with a significant degree of immunity postinfection) had immunity that lasted for up to a year, as well as, more recently [6], on historical patterns of three common coronaviruses.
Figure 1
The SIDUHR^{+∕−} model. 
2.1.2 Dynamics
The evolution of the sizes of each compartments is assumed to be modelled by the following system of ordinary differential equations: (2.1)
Let us explain the different parameters. The transmission rate (1 − δ_{t})β is the rate at which an undetected infected individual transmits the disease to a susceptible one. Multiplying this parameter by , which is proportional to the number of pairs that we can form with an S and an I^{−} individual, provides the total infection rate at the level of the population and at time t: . Heuristically, this quantity can also be understood as the fraction of the population infected during a unit time interval at time t. The infection rate β is reduced by a timedependent factor (1 − δ) where δ ∈ [0, 1] comes from social distancing, quarantine, isolation interventions that reduce infectious contact rates. The strength of the lockdown intervention interprets as follows: δ = 0 means no intervention while δ = 1 means no social interaction. We consider two ways of detecting individuals: short range tests allow to discover new infectious individuals in I^{−} and we denote by λ^{1} the rate at which each of these infectious individuals gets detected, long range test allow to find formerly infected and hereby immune individuals that are now in R^{−} and we denote by λ^{2} the corresponding rate. Here we are interested in the case where the parameters δ and λ^{1} are control variables and in the sequel, these parameters are hereby time dependent.
The recovery rate for an infected individual is γ_{IR}. The rate at which an infected individual enters the health system is γ_{IH}. We assume here that both rates γ_{IR} and γ_{IH} are the same for detected or nondetected infected individuals. The individual transition rates from H to U or R^{+}, from U to R^{+} or to D are denoted by γ_{HU}, γ_{HR}, γ_{UR} (U) and γ_{UD} (U) respectively. Notice that the twolater rates γ_{UR}(U) and γ_{UD}(U) are not constant but depend on U_{t}: when the ICU limitation is reached, the extra patients who are not taken in charge by the system can not recover and are exposed to extra death rates, as explained in the Appendix A.
These parameters are chosen following the tracks of [26, 73]: (see Tab. 1). In these papers, some parameters are fixed according to clinical studies and other are estimated from French data. The detail on how these parameters are computed in the present paper is given in Appendix A. However, we emphasize that although we choose French benchmarks, our input is more on methodology and can be of interest for other countries or epidemics as well.
The system described here can be derived from individual based stochastic processes (e.g. [9]) and thus, there is an underlying individual based model where all the rates have a probabilistic and statistical interpretation (see Appendix A), with connection to continuous time Markov chains. As we will discuss later on, those parameters can therefore be interpreted as inverse of transition time lengths.
The initial time t = 0 corresponds to the time at which policies that reduce the infection rate are put in place (δ > 0), or the time at which testing is started (λ^{1} > 0). For the French COVID19 epidemic, it would correspond to March 17th 2020, but again, we emphasize that our methodology remains of interest for other countries or diseases (see e.g. [63] for country specific informations).
2.1.3 Basic and effective reproduction numbers ℜ_{0} and ℜ_{t}
The basic reproduction number ℜ_{0} corresponds to the expected number of individuals directly contaminated by a typical infected individual during the early stage of the outbreak (when the proportion of susceptible individuals is close to 1), see [23]. When ℜ_{0} < 1, the disease doesnot spread quickly enough, resulting in a decay in the number of infected individuals. But when ℜ_{0} > 1, the infected population grows over time.
As recalled in [77], the value of the ℜ_{0} is modeldependent. For the model considered here, only the I^{−} individuals can transmit the disease to infect susceptible individuals. If we see I^{−} as the infected class, and if we merge the classes I^{+}, R^{−}, R^{+}, H, U and D into a large ‘recovered’ class, in this SIR type of model, the dynamics of has infection rate (1 − δ)β and removal rate λ^{1} + γ_{IR} + γ_{IH}, assuming that δ = δ_{0} and are constant over time. In our model, the basic reproduction number is (2.2)
Without intervention (i.e. for (δ, λ^{1}) = (0, 0) over time) and with our chosen parameters provided in Table 1, this corresponds to ℜ_{0} = β∕(γ_{IR} + γ_{IH}) = 3.3 which is in line with [73]. Sensitivity of our results with respect to variations of ℜ_{0} within the range [3.3, 3.6] is discussed in Appendix B.2.
In the case where λ and δ evolve with time, we end up with a dynamic reproduction number that is called effective reproduction number, see [8]. More specifically, at time t ≥ 0, let us define (2.3)
Its interpretation is very close to the static version ℜ_{0}, when assuming S_{0} = 1, which is a reasonable first order assumption since I_{0} was assumed to be extremely small. When ℜ_{t} < 1, the average number of secondary cases started from one primary case with symptom onset on day t, dies out quickly, resulting in a decay in the number of infected individuals in the population. But when ℜ_{t} > 1, the infected population grows over time. So we should expect to have controls on the quarantine level δ_{t} and on the testing rate λ_{t} that will constrain ℜ_{t} to be reduced beyond 1, at least after some starting time.
Another information brought by the ℜ_{t} < 1 is that, when S is smaller than the socalled ‘herd immunity’ threshold, the epidemic enters into a subcritical phase where it goes to extinction. In our case, the herd immunity threshold S^{*} without intervention is around 30% and durable lockdown and detection efforts can allow to increase it significantly until the arrival of a vaccination solution.
2.2 The optimal control problem
In the previous section, we described our epidemiological model in order to characterize the dynamics of the virus spreading processes. We now consider a government acting as a global planer, who wants to control the epidemic dynamics in order to balance the induced sanitary and economic outcomes. In order to mitigate the effects of an epidemic, several parameters in the model are now interpreted as ‘control levers’, and we now focus on the derivation of their optimal dynamic choice. Optimizing intervention strategies is an important policy issue for the management of infectious diseases, such as COVID19. Optimal control in the context of pandemic models has been used since [1, 2], [14] or [74], in the 70’s.
2.2.1 Possible control levers
Quarantines and lockdown
Lockdown measures proved to be effective in controlling the COVID19 outbreak in China, as recalled in [53] or [56]. Quarantine is a rather old technique used to prevent the spread of diseases. It refers to the restriction of movement of anyone (not necessarily sick people), and it usually takes place at home, and may be applied at the individual level. It should be distinguished from isolation, which refers to the restriction of movement of infectious individuals who have a contagious disease. It can be done in hospitals, in dedicated facilities, or at home. Those two are seen as active controls, and are different from contact surveillance, which is passive. Modern quarantine’s goal is to reduce transmission by increasing social distance between people, by reducing the number of people with whom everyone can be in contact with. This includes canceling public gatherings, closing public transportation, etc. A lockdown is intended to stop people from moving between places, and it could involve cancelling flights, closing borders, and shutting down restaurants. The idea is to reduce the flow of people to curb transmission, and we might consider here that lockdown and quarantine measures are almost equivalent.
Reference [74] suggested an optimal quarantine problem, in an SIR model, where a proportion of infected people could be quarantined. He proved that the optimal strategy was either to quarantine all infected individuals (if so, at early stages) or none. In the context of the 2003 SARS pandemic in Singapore, [62] recalls that a Home Quarantine Order was signed, but it was more an isolation of sick people than a lockdown as the ones used to lower possible consequences of COVID19. But options considered are different, from one country to another. Reference [81] mentions that several countries closesd their borders, in most European countries, schools and restaurants were closed for weeks.
Reference [55] tried to quantify the effectiveness of a quarantine strategy, where healthy people are advised to avoid contacts with individuals that might carry the disease, using network based models. Here, the parameter has a direct interpretation, since it is related to the average rate of contact between individuals. In [83] or [51], quarantine and isolation are used as control. In our model, those controls correspond to the δ_{t} component. As in [15], such control can also be related to travel bans, for a weaker form than strict quarantine. Note finally that [7] suggested that a complete lockdown was impossible to reach, and that an upper bound of 70% (of the population) should be considered. In this paper, we do not incorporate any constraint (except that we cannot lockdown more than the entire population), but interestingly, in the numerical simulation, the upper bound we obtained as optimal is very close to this value.
Testing, tracing and isolating
The second most important lever that can be used is detection. As claimed by several governments (see [35]) in the context of the COVID19 pandemic, testing is key to exit lockdown, and mitigate the health and economic harms of the virus. Here, we assume that health authorities can perform two kinds of active detection actions.
Type1 test, that could be an antigen test (Reverse Transcription Quantitative Polymerase Chain Reaction (RTqPCR)  also called molecular or PCR – Polymerase Chain Reaction – test). This test allows to detect whether an individual is currently infected or not. It should be used to determine who selfisolates and for contact tracing. For that test, a sample is collected – usually with a deep nasal swab (and analysed in a laboratory). This is a shortterm test: when performed at time t on an individual, we know whether an individual is infected or not. Together with proper detection means, the use of this test corresponds to the control in our dynamics. This test allows to detect and isolate infectious individuals. In our model, they transfer from compartment I^{−} to compartment I^{+}. It is performed as mandatory when someone arrives at the hospital. Hereby, any individual being at some point in compartment I^{+}, H or ICU ends up in the identified immune compartment R^{+} if he survived the epidemic. We assume for simplicity that those tests do not have false positive or false negative and leave for further research the impact of the test sensitivity on the epidemic dynamics and optimal control. Thus, the tests are all positive on individuals in I^{−} and negative on individuals in S and R^{−}. Other individuals will not be tested.
Type2 test is an antibody test (using serological immunoassays that detects viralspecific antibodies – Immunoglobin M (IgM) and G (IgG) – also called serology or immunity test). It allows to test whether an individual can now be considered immune to the virus. It could potentially be used to issue immunity certificates in order to help restarting the economy quicker. It also reduces the overall uncertainty and potential fear of individuals to be infected, and has a positive impact on the economic consumption. It is also useful for contact tracing purpose as it helps identifying individuals that can not become infectious again. For such a test, a blood or saliva sample is applied to a strip that identifies presence of antibodies. We assume here that immunity lasts long enough with respect to the epidemic duration, and R^{+} regroups the collection of immune detected individuals. It is a longterm one: when performed at time t on an individual, we know if that person had the disease before. The detection using those tests identifies in our model to the control variable λ^{2}. It allows to transfer individuals from compartment R^{−} to compartment R^{+}, but again, people in S will be tested also, as well as people in I^{−}. As for Type1 tests, we will suppose here that these tests do not make any wrong identification.
For type2 tests, a small sample of a patient’s blood – for instance via a pin prick – and the test looks for two specific types of antibody: IgM and IgG. IgM are the first antibodies to be produced by the immune system. They have a halflife of around five days, and they usually appear within five to seven days of infection and peak at around 21 days. Detection of these antibodies suggests theperson has existent or a recent infection. IgG antibodies are more numerous and can be detected around 10 to 14 days after infection. The presence of these antibodies indicates a person has recovered from the virus and is now immune.
‘Raising the line’ and ICU capacities
So far, we have assumed that the sustainable capacity of ICUs was not a control, and should be considered as an exogeneous fixed boundary (denoted U_{max}). As discussed in [34], preventing a health care system from being overwhelmed requires a society to do two things: ‘flatten the curve’ – that is, slow the rate of infection so there aren’t too many cases that need hospitalization at one time – and ‘raise the line’ – that is, boost the hospital system’s capacity to treat large numbers of patients" (see also [10]). Therefore, a natural control variable is the level of the line, or U_{max} as we named it.
Nevertheless, we will not consider U_{max} as a control variable that can be optimized (such as the strength of the lockdown, or the effort in detection), since it is difficult to assess the cost of raising it. Nevertheless, in Section 3.4, we consider the case where health authorities choose to increase that capacity, by making ventilators and other medical material necessary inthe context of the pandemic. But it cannot be increased indefinitely since ICU also require trained personal and even if adjustments can be considered, we assume that it is not realistic to assume that U_{max} can be increased by more than 50%.
2.2.2 Objective function
We now turn to the design of the objective function, trying to take into consideration the sanitary and socioeconomic outcomes of the lockdown and detection policy. In the context of planning vaccination campaigns, [21, 79] suggest a convex quadratic cost function by minimizing both the number of infected individuals in a time horizon and the cost to implement the control policy. In [48], a model for 2009 A/H1N1 influenza in Korea is considered: the goal is to minimize the number of infected individuals and the cost of implementing the control measures, and the cost is taken to be a nonlinear quadratic function. Quarantine and vaccination are considered as control variables in [45], and the optimal control is obtained by minimizing a quadratic cost function. Inspired by those approaches, we will also consider quadratic cost functions. Before going further, let first introduce some relevant quantities characterizing the epidemic phase that will reveal useful in the upcoming analysis.
Quarantined individuals
Let denote by Q the quantity of people concerned by the lockdown policy and defined by: (2.4)
All the individuals in Q_{t} are having the level of contact rate δ_{t} with the population at time t. This proportion of individuals identifies to the people for which we can not say if they already contracted the virus or not. Therefore, they are all concerned with non targeted lockdown strategies as well as possibly randomized testing trials. Recall that, on the opposite, we suppose that infected detected individuals in I^{+} are isolated andhave no social interactions with the population, whereas detected immune individuals suffer no mobility constraints and can reach a level of local interactions similar to the one they had before the beginning of the pandemic.
Global level of social interactions
On the other hand, the global level of social interactions among the population is denoted W_{t} and given by: (2.5)
The quantity W_{t} represents the proportion of social interactions in the population and can be interpreted as a macroscopic labour force level for the economy at time t.
Testing resources
Finally, in order to measure the detection efforts of infected and immune individuals, we introduce the following metrics; which identifies to the number of virologic and antibody tests done in a (somehow quite unrealistic) fully randomized detection trial: (2.6)
Observe that the different formulation between both metrics is due to the fact that individual are automatically tested using virologic detection methods, when admitted at the hospital.
Vaccination arrival date
All health authorities hope for the development of a vaccine for COVID19 in a close future. In order to encompass uncertainty concerns upon the arrival date of the vaccine, we consider the arrival date of the vaccination solution to be a random time with exponential distribution of parameter α, denoted τ. According to recent studies, [19, 46] the creation of such vaccine for a sufficient quantity of individuals can presumably be assumed to require around 500 days. We shall consider a parameter α equal either to 0, 1∕250 or 1∕500 for our numerical experiments. For convenience, we assume here that both a vaccine and a cure simultaneously appear at time τ. Reference [7] also uses such a discount approach, on top of some economic discount rate (assuming that the vaccine will arrive at an expected date of 18 months).
Sanitary cost
For the sanitary cost, it seems rather natural to simply consider a death count metric given by D_{T}. Nevertheless, [61] legitimates the use of an appropriate discounting of health outcomes. About the kind of discount we should use, [78] discusses the use of hyperbolic discount (opposed to classical exponential discount), in the case of social choice regarding health outcome. The hyperbolic discount allows for asymmetry in discounting, and nonstationarity, in the sense that postponing in two years now can be substantially different from postponing in three years (from now) next year. On the other hand, [78] claims that is not necessarily a realistic assumption, so that using a standard timeconsistant exponential discount makes sense in terms of health outcomes. Besides, the longer we can postpone deaths, the more likely a vaccine will be found before. Adding up a exponential discounting factor reveals to be also perfectly consistent with the consideration of random development time of a vaccine solution at large scale, with exponential distribution. Hence, we will consider a sanitary cost of the form: (2.7)
recalling that the random arrival time of the vaccine follows an exponential distribution with parameter α.
Economic and social cost
The quantification of the socioeconomic induced by lockdown restrictions is not easy to quantify properly. More involved modelling may wish to take into account a feedback effect between labor and consumption levels together with the epidemic propagation characterized by the global rate of social interactions within the population. Such approach goes beyond the scope of this paper and we solely intend to quantify the impact of the epidemic prevalence in terms of reduction of social interactions. Here, we would like to capture the loss of productivity and wellbeing among the population due to lockdown situation. In a very general setting, the production function is a function of capital input (K, that usually corresponds to machinery, equipment, buildings, etc.) and labor (W, that represents the total number of people working): the output produced in a given period of time is classically given by a CobbDouglas production function function (as introduced in [18]) K^{1−a}W^{a}, with a ∈ (0, 1). Power coefficients identify to the output elasticities. Here, the function display constant returns to scale, in the sense that doubling the use of capital K and labor W will also double produced output. Assuming that, on the short term (say less than a year) capital remains constant, it means that if at time t, only a proportion W_{t} of individuals can be economically active, the productivity loss is proportional to . Thus, the intertemporal economic cost at time t is proportionalto , or (up to an additive constant) . Nevertheless, several other metrics have been used on the literature in order to account for the economic cost induced by the lockdown situation: [38] or [28] considered a linear function, while [13, 65] or [25] considered a convex function. As for the other cost structures, we choose for simplicity to keep a quadratic structure quantifying the distance between the current global rate of social interactions W given in (2.5) (and identified as a metric for welfare, wellbeing or workingforce) and its optimal value 1. Hereby, the social and economic cost will be given by: (2.8)
where we recall that the discounting rate α encompass the distribution of arrival time of a vaccination solution solution.
Detection cost
The quantification of efforts put into the detection of infectious individuals is a difficult task, as outlined in [57] in the case of children, but several detection means are possible. As mentioned in [72, 42], testing is only part of a strategy. The World Health Organization recommends a combination of measures: rapid diagnosis and immediate isolation of cases, rigorous tracking and precautionary selfisolation of close contacts. Reference [20] compares various techniques used to detect infected people. Recall that we assume that a virologic testing is performed for all individuals entering into the hospital. Besides, proper tracing of contamination chains together with potential random screening of the population should be conducted on the non detected individuals, i.e. on the subpopulation Q_{t}. This group contains the susceptible S_{t}, the undetected infected as well as the undetected immune . As antibodies detection should also be performed on the same group Q_{t}, it seems important to take into account the impact of the size of this group within the population. Besides, as testing resources are sparse, the metric quantifying the effort of detection must take into account the increasing difficulty for scaling massively testing resources. For all these reasons, we pick a cost for detecting infected or immune individuals of quadratic form as follows:
where we recall that and are given in (2.6). Again, we consider here a detection cost, not a testing cost.
2.2.3 Optimal control under ICU capacity constraint
We are now in position to turn to the design of the global objective function. We are looking towards the societal optimum chosen by an omniscient government acting as a global planer and is facing a multiobjective control problem. As in [54], we shall consider a weighted sum of different costs:

the sanitary cost described in (2.9) directly related to the mortality rate of the virus;

An economic and social cost given in (2.8) due to the reduction of social interactions within the population;

A detection cost for identifying infectious and immune individuals; which will be mainly discussed in Section 3.5 below.
Combining these effects, the global objective function is given by: (2.11)
where the w_{i} terms are positive weights associated to each marginal cost. The choice of proper weights w_{i} is a hard ethical task, that we choose to avoid as out main purpose is to highlight the main patterns of the optimally controlled epidemic dynamics. We pick some weight levels for our numerical experiments and study in Appendix B.3 the sensitivity of our findings to strong variations of these weights.
In our approach, the global planer tries to minimize the global cost function given in (2.11) using the lockdown control lever δ together with the detection lever λ = (λ^{1};λ^{2}). Nevertheless, we did not yet take into account the risks and outcomes induced by the saturation of ICU hospital care facilities. The dynamics of the SIDUHR^{+∕−} model already takes into account the strong impact of ICU saturation on the mortality rate on individuals with severe symptoms. Whenever the proportion of individual in U exceeds the upper bound level U_{max}, patients can not be treated correctly so that patients are more likely to transfer to D instead of R^{+}, in comparison to calmer times where U < U_{max}. With such feature, the optimal lockdown and detection strategies should undoubtedly try to limit the duration of periods where the maximal ICU capacity is exceeded. Nevertheless, as we believe that exceeding capacities of ICUs would have a strong negative impact on the hospital staff and even the entire population, we rather choose to impose the nonsaturation of ICU as a required state constraint on the system. Hence, we try are looking for a robust solution to the following control problem: (2.12)
where denotes the set of admissible strategies driving the dynamics of the SIDUHR^{+∕−} model and is given by (2.13)
2.2.4 Numerical approximation of the solution
From a mathematical point of view, the optimal control problem of interest falls into the class of control for deterministic dynamical system over infinite horizon in the presence of additional state constraint. The detailed presentation of the numerical resolution of this problem is presented in Appendix C.2. We briefly recall here the main underlying approximations used in the numerical approach.
In order to optimize the impact on the dynamical structure of the underlying SIDUHR^{+∕−} of the chosen control, we make use of Pontryagin maximum principle. For numerical purposes, we restrict our analysisto the consideration of a control problem with finite horizon. We thus pick a maturity T large enough in order to have a very small remaining level of infected individuals at time T. Besides, whenever the level of infected is small enough, SIR type dynamics are not fully reliable anymore, limiting the impact of such approximation on the derived optimal strategy. In practice, the numerical experiments are indeed not sensitive to the choice of maturity T, as soon as it is picked large enough. For the graphs shown here, we used T = 700 days (almost two years), corresponding to a scenario where a vaccination solution is available within 2 years. We did run computations up to T = 900 looking for a potential impact on the interpretation of the optimal strategy. But since we did not observe relevant ones, we restricted our attention to T = 700 in order to simplify the overall presentation of the results. Besides, in order to encompass the state constraint into the maximum principle based algorithm, we simply represent it under the form of a penalization cost, whenever the ICU constraint is not satisfied. The penalization is strong enough to discourage crossing this threshold for most scenarios. See Appendix C for more details on the implementation.
The algorithm then relies on an iterative procedure. Time is discretized so that the search for an optimal control reduces to the search for one vector per control variable (δ, λ^{1}, λ^{2}), with one value per time step. Starting from an initial guess, at each iteration the approximation for the optimal control is updated based on the expression of the gradient in terms of the solution to a forwardbackward system of ordinary differential equations. In order to ensure that the controls remain between 0 and 1, a projection of the new controls on the interval [0, 1] is performed at each iteration.
3 Discussion
3.1 Benchmark scenario without intervention
Let us first investigate the outcomes of the epidemic dynamics when no lockdown, detection nor quarantine policies are available. As commonly done in the current literature, we do not take into account the reaction of individuals in response to the virus prevalence dynamics, as discussed e.g. in [30]. Anticipating realistic individuals behaviors in such context is still difficult, as no available data allows to tackle this phenomenon properly. We suppose that the virus dynamics has no impact on the level of social interactions of individuals,which remains constant equal to β_{0}. Such analysis provides a very useful nointervention benchmark scenario for the rest of our study. The dynamics of each compartment is provided in Figure 2 while Table 2 provides important epidemic characteristics computed at terminal date T = 700 days.^{1}
As no intervention interferes with the epidemic dynamics, the overall duration of the epidemic is very short as the number of susceptible vanishes exponentially fast. The initial reproduction number ℜ_{0} = 3.3 induces an epidemic lasting around 80 days. The virus prevalence is quite important over that period, leading to the contamination of almost 96% of the population, way above the herd immunity threshold, S^{*}~ 71% in our example. The epidemic peak happens after a one month period, and as much as one third of the population is infected at that date. Around 7% of the population is hospitalized due to the virus, leading to a highly severe overwhelming of the healthcare capacity system for around 60 days. The epidemic induces in total a mortality of almost 1% of the population, while half of deaths occur while the ICU capacity is saturated.
It is worth noticing that our model, with the parameter values of Table 1 that were calibrated on the French epidemic of COVID19 ([26, 73]), gives predictions that are in accordance with the ones made by [70] in absence of interventions and neglecting the ICU capacity constraint (see Appendix A for details).
Figure 2
Evolution of states without control, over 200 days. 
Epidemic characteristics.
3.2 Optimal lockdown policy without detection effort
In order to emphasize the marginal role of each component in the optimal control problem of interest, we first focus on the optimal lockdown strategy and corresponding epidemic trajectory induced without any detection effort. This does not mean that we do not test anyone, we simply assume that only the people entering in the hospital are tested with type1 tests. No testing is done on infected individuals without severe symptoms (λ^{1} = 0) and there no testing is done to detect recovered people after being asymptomatically sick (type2 test, λ^{2} = 0).
Hereby, we are starting with a weaker form of the optimal control problem, where we control only the lockdown intervention policy: (3.1)
with the following choice of parameters and w_{sanitary} = 100000, w_{eco} = 1 and α =0. The sensitivity of the enhanced results to this adhoc parametric choice is discussed in Appendices B.3 and B.4. Besides, once we will have a better understanding about the derived lockdown strategy and its robustness, we will focus on the impact of detection strategies in Section 3.5. The optimal lockdown strategy together with the main characteristics of the pandemic dynamics are presented in Figure 3 and Table 2. The optimal control problem has been solved using the numerical algorithm described in Appendix C.2 and numerical approximation errors undoubtedly remain in the results presented here.
When we follow the lockdown intervention (δ_{t}) proposed bythe optimal response policy, four distinct consecutive phases can be observed:

Quick activation of a strong lockdown: The optimal strategy consists in activating as soon as possible a strong lockdown measure on the entire population as no targeting or detection is considered or possible at that stage. This strategy has been widely used in March and April 2020 as half of the world population was under lockdown. With our parameters, the contact rate over the population is reduced by 80% over a couple of days, in order to stop immediately the exponential growth of infected and infectious individuals. Hereby, the size of the epidemic peak is highly reduced and, after less than one month, only 2% of the population is contaminated at the peak of the epidemic (instead of 34% in the nointervention benchmark case). Of course, at the exact same date, the effective reproduction number ℜ_{t} goes below 1. During this initial short period, the effort in economic activity reduction has been very important: the welfare level W went down by 80%, while at the same time the occupation rate in ICU hospital care system has been constantly rising.

Light lockdown release and prevalence decrease: During the short second phase, the intensity of the lockdown remains intense (above 60%) in order to keep ℜ_{t} below 1. The prevalence of the virus decreases slowly, while the level of ICU admissions is still rising. After the previous phase where urgent decisions had to be taken by the global planner, a short transition period of a couple of days now starts. During this period, the proportion of infectious individuals decreases strongly in order to reach an optimal prevalence level I^{*} in the population, required for the very long next phase.

Long period with stable prevalence and ICU sustainable capacity: The third phase lasts around one year during which the effective reproduction number ℜ_{t} remains stable at level 1. Hereby, the prevalence of the virus remains stable while the ICU capacity remains at the sustainable level U_{max}. During this one year period, the level of contact rate within the population is slowly growing at a regular constant rate, moving gradually from 40% to 80%. The working force is rising back to more normal levels and the economy is slowly restarting: one year after the beginning of the pandemic, around 70% of regular social interactions level is already attained within the population. The mortality rate over this period remains stable, and the death toll grows linearly, almost reaching its terminal value at the end of this period.

Terminal slow progressive release of the lockdown: The last phase is a couple of months long and consists in bringing progressively back to normal the level of social interactions, while the occupation level in ICU is decreasing and the prevalence of the virus is heading towards 0. The effective reproduction number ℜ_{t} is smaller than 1. At some point, the deterministic equation (2.1) are no more valid and the system becomes stochastic, but the extinction of the disease occurs with probability 1.
The optimal structure in four steps for controlling the pandemic dynamics is quite robust and very natural. It is worth observing that although the dynamic reproduction number is not part of the criterion of interest, reaching its optimal trajectory with four successive patterns is key for controlling the dynamic level of admissions in ICU. Using the optimal control policy, the infections by the virus in the population have been spread out over 500 days, providing a sustainable level of admissions in ICU for the sanitary system. In comparison to the nointervention benchmark case, the death toll has been divided by almost 100, thanks to a strong reduction in the level of social interactions. The size of the pandemic only represents 70% of the population in comparison to 96% in the nointervention scenario.
A nice visualisation of the epidemic dynamics is provided by looking at the phase diagram presented in Figure 4. In the nointervention scenario, the epidemic dynamics starts from the point (S, I) = (1, 0) and moves up and left until reaching the herd immunity threshold S^{*} corresponding to the time of epidemic peak. Then, I starts decreasing until the end of the epidemic, point where S is valued around 4% in our example: hence, the final epidemic size is almost 96% of the population. In the optimally controlled scenario, we can clearly identify the 4 phases described above: I is still rising during the first phase until reaching a small epidemic peak level, it starts decreasing during the second transition phase until reaching the optimal level I^{*} corresponding to the level U_{max} of ICU capacity. It remains stable there during the long third phase and finally decreases in the last phase until reaching a susceptible proportion close to the herd immunity level S^{*}.
Of course, both dynamics depend on the value of the parameters, such as the initial virus prevalence , the starting basic reproductive number ℜ_{0}, expectations about possible arrival of vaccine and cure (through parameter α or weights considered in the objective function). But as we can see in Appendix B, and as we will discuss more now, the enhanced results are robust to the parameter specifications.
Figure 3
Evolution of state proportions – S_{t} (susceptible), I_{t} (infected, detected or not, ), R_{t} (recovered and immune, detected or not, ), U_{t} (hospitalized in ICU), D_{t} (dead), ℜ_{t} (dynamic version of the reproduction number), W_{t} (labor force available), Q_{t} (susceptible to be quanrantines, ) and δ_{t} (lockdown intensity) – with optimal control δ (plain line) and without control (dashed line). States vary here from t = 0 (beginning of the pandemic, or at least of possible measures, with δ ≥ 0) to t = 700, and computations are based on a dt = 1∕5 days. 
Figure 4
Phase diagram S_{t} (susceptible)versus I_{t} (infected, detected or not, ), together with the herd immunity threshold S^{*} and the levelI^{*} of infected corresponding to U = U_{max}. 
3.3 Sensitivity of the optimally controlled epidemic dynamics to the model specifications
The calibration of SIR type dynamics to current available data is currently a challenging task, as the quality of available mortality and hospitalization figures is still questionable and incomplete. It is worth noting that one of the main advantages of the SIDUHR^{+∕−} dynamics is the possibility of calibrating its dynamics to mortality, infected detected, hospitalized and ICU data points. Still, such inverse calibration problem remains illposed and the precision range over the calibrated parameters is not satisfying.
In order to demonstrate the robustness of our findings in such context, we provide a sensitivity analysis of our results to modifications in the model specifications. More specifically, we analyse the impact of reasonable variations with respect to the initial prevalence in Appendix B.1, the basic reproduction number ℜ_{0} in Appendix B.2, the balance between sanitary and socioeconomic cost in Appendix B.3, and the arrival date of a vaccination solution in Appendix B.4. Let now present the main outputs of those numerical experiments, while keeping in mind that numerical errors due to the optimization approximation scheme can still alter some numerical results presented in this study.
On Figure B.1, we can visualize the impact of a change in , with a smaller (five times) and larger (also five times) value than the one used in benchmark computations (we used %_{0}). The lower value was suggested from [71], with 1%_{0}, five times lower than our estimate, and an upper value five times upper, of 2.5%, close to the situation two weeks after according to [71]. A smaller value means that measures were taken earlier, but overall, it looks like the overall impact would have been rather small. Note that we observe here that a smaller might cause more deaths (D_{T}), because of some balance with economic gains (economics losses are smaller with a lower ). Observe also that variations in have almostno impact on the ICU saturation dynamics: it might start earlier with a large , but the duration remains more or less the same.
On Figure B.2, we can visualize the impact of a change in ℜ_{0}, with a smaller (− 0.3 i.e. − 10%) and larger (+ 0.3 i.e. + 10%) value than the one used in benchmark computations (ℜ_{0} = 3.3). Such variation in ℜ_{0} is compatible with the literature, see [26, 73]. The overall 4phases global patterns of the optimally controlled epidemics remain identical. The lockdown effort increases slightly but significantly with ℜ_{0}, together with the ICU saturation phase duration, the global death toll and global economic cost. Note that a decrease of ℜ_{0} from 3.6 to 3.0 (attained for example through extra hygiene habits) would decrease ICU saturation time length by almost 25%.
On Figure B.3, we can visualize the impact of a change of the ratio between the sanitary weight w_{sanitary} and the socioeconomic one w_{eco}, with a smaller (half) and larger (twice) value than the one used in benchmark computations. The global patterns in the epidemic dynamics remain globally similar. A larger sanitary weight induces a smaller death toll, a smaller epidemic peak together with a shorter ICU saturation period. The economic gain or loss seems to be mostly sensible during the second half of the epidemic duration.
Figure B.4 provides the shapes of numerically approximated optimal epidemic curves for different beliefs about the expected arrival time of a vaccine (and a cure). In comparison to the benchmark case where no vaccine potential creation is taken into account, the arrival date τ of the vaccine is supposed to follow an exponential distribution with expectation 100, 250 or 500 days, while the optimization is occurring on the time interval [0, T ∧ τ]. The epidemic curves are very similar when the anticipated time for vaccine availability is long enough, confirming the robustness of the approach and main findings. Whenever a vaccination procedure is supposed to arise rather early, the optimal strategy consists in postponing the mortality while keeping a stronger and longer lockdown effort, hoping for a quick vaccination solution. The sanitary burden is hereby smaller during the early stage of the epidemics but rises more strongly at the terminal phase.
Our main overall observation is that the dynamics of the optimally controlled epidemic dynamics always present the same patterns and enlights the robustness of our findings. It divides into four successive phases in a similar way to the one described in the previous main scenario analysis. Besides, the effective reproduction number ℜ_{t} together with the optimal lockdown intervention policy also have the exact similar patterns. Of course small variations in the specifications still have an influence on the exact figures characterizing the balance between sanitary and economic outcomes, the total duration of the epidemic or the pressure on the ICU hospital system, but the overall characteristics remain identical. These observations confirm the robustness of our findings when dealing with the optimal control problem at hand, under specification uncertainty. The more specific dependence on the ICU care system capacity is studied in more details and discussed in the next section.
3.4 On the impact of additional ICU capacities
As discussed earlier, most countries have been urged to take public health measures to ‘flatten the curve’. An important motivation was to avoid that hospital systems get overwhelmed, and more specifically, social distancing was introduced to ensure that the occupation of ICUs with ventilators remain within the capacities of hospitals, dedicated to COVID19 patients. Previously, we integrated an upper limit to ICU – denoted by U_{max} – and in this section, we discuss the impact of a potential increase of that capacity. Of course, increasing capacity does not simply mean getting more masks, beds, ventilators, or medicines for the serious cases. It also means diverting resources, training more staff to work in the ICU, and probably introduce telemedicine solutions for nonCOVID19 patients (who do not require hospitalization) to free up beds.
On Figure 5, we consider a change in U_{max}, where we assume that health authorities are able to increase ICU capacities by 50% over a 2 months period and remain at that level. We can imagine an increase in the number of ventilators or backups with staff coming from other regions that are spared by the epidemics, as was done in China or France for instance. Another possible interpretation is that the medical staff on the front line might be infected in early periods leading to an initial shortage of personal at the beginning followed by a return to normal within a few weeks. Here, we consider a regular (linear) increase of ICU capacities over two months, until reaching a new limit, 50% higher than the original one.
This scenario has a major impact on all quantities. First, we observe that the derived optimal lockdown strategy is able to let the ICU admission level bind the sustainable constraint over time, although it now grows dynamically with time. This corresponds to a longer second phase duration until the stable in prevalence and ICU capacity starts. The final number of dead people is overall smaller than in the original scenario, as the expansion of ICU capacity allows to treat more patients simultaneously. Therefore, the social distancing can be much weaker during the long ICU saturation phase, which also happens to be way shorter. Furthermore, ICU saturation duration is much shorter with a major positive impact on the economy from day 100 until day 400, as the lockdown intervention remains much weaker (after 200 days, the global level of social interactions is around 65% in comparison to less than 50% in our benchmark scenario). This reduces the socioeconomic impact of the epidemic and hereby allows to reach a smaller mortality count at time T. The overall gain at terminal date is massive (even if here, we do not include a cost for increasing capacities). Nevertheless, the augmentation of ICU capacity implies also a higher daily dynamic death toll, so that an early nonanticipated vaccination solution at an intermediate date would induce a higher sanitary cost in comparison to the scenario without ICU capacity rising.
Figure 5
Evolution of states with optimal control δ when ‘raising the line’ (increase of U_{max} by 50%). The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to the scenario with a capacity of U_{max} which increases linearly from time t = 0 until t = 60 days up to 1.5 times the initial capacity and then stays at this level. 
3.5 On the impact of detection resources
We now focus our analysis on the impact and optimization over the detection efforts λ^{1} and λ^{2} of infectious or immune individuals. As stated in [64], mass testing could end the epidemic rapidly, and we now try to emphasize the impact of light or mass testing upon the lockdown optimal intervention policy. The numerical results are presented in Appendix B for ease of readability and we present here the main outputs of our experiments.
3.5.1 Effort in virologic detection
We first quantify on Figure B.5 the impact of introducing additional virologic detection resources with fixed success rate λ^{1}. Namely, we consider a constant daily detection level λ^{1} of 1%_{0}, 1%, 5% or 10%. Of course, the more effort we put into detecting infected people, the better the output. For instance we observe a smaller death toll, a weaker economic impact as well as less ICU saturation. Up to remaining numerical errors, moving the detection capacity λ^{1} from 1% to 10% has a strong impact on the overall outcomes of the epidemic: The epidemic size drops from 70% to 30%, the remaining number of susceptible individuals S_{T} ≃ 70% corresponds to the herd immunity obtained with λ^{1} = 0.1 and δ ≃ 0.25, the duration of ICU capacity saturation is divided by a factor larger than 2, and the overall sanitary and economic costs are divided by 2. The huge detection of infected individuals has a strong impact on the socioeconomic cost for two reasons: First, detection of infected individuals allows to slow down the epidemic propagation. Second, detected infected individuals are later on identified as immune people, i.e. are in compartment R^{+}, so that they are not impacted by social relations limitation, which reduces the overall economic cost. Finally, observe that massive testing has a major impact on most output (and must be massive to really impact, as suggested in [64]). Observe that S_{T} remains rather large, but hopefully any epidemic restart is avoided by the detection effort λ^{1} which must not be released (the patterns of some figures close to terminal date T are due to the numerical approximation of the infinite horizon problem by a finite one).
Another experiment presented in Figure B.6 considers the case where the detection resources are rising linearly from 0% to 20% aver 700 days. In comparison to the corresponding case where detection resources of level λ^{1} = 10% are constantly available, the sanitary death toll induced by a rising λ^{1} is way higher. This confirms the intuition that intense detection resources are required from the beginning of the epidemic. More specifically, we now focus on the optimal dispatch over time of allocation resources: Figure B.7 provides the optimal policy for the detection of infected whenever no lockdown measures are in place, while Figure B.8 provides the optimal combination of detection policy and lockdown intervention. For both cases, the weight w_{prevalence} has been set up equal to 1. As we can observe on both Figure B.7 and Figure B.8, and as stated clearly in [64], massive testing and detection strategies could have been used to end the epidemic rapidly. The optimal strategy consists in having extremely massive and early type1 detection resources, in order to isolate infected individuals and break the contamination chains. Herd immunity is ensured while we keep the detection active so that ℜ_{t} remains below 1, allowing to epidemic size to remain very small, while everyone simply waits for for a vaccination cure.
With that strategy, the peak on ICU is 25% of the upper limit U_{max}, D_{T} is almost null, with no real impact on the economy (there is still an impact since, in order to be safe, we still quarantine a significant proportion of people, above 10%). Note here that the optimal λ_{1}, once we passed the peak, exceeds 25%, which is higher than the scenarios considered earlier. Observe also that massive detection tools induce a constant optimal lockdown level.
3.5.2 Effort in immunity detection
We now turn to the impact of detecting immune individuals, as found by the optimally controlled epidemic dynamics. First, we consider on Figure B.9 the introduction of fixed type2 detection strategy (or testing) and measure the impact on the optimal lockdown policy. Such detection resources allow to transfer individuals from compartment R^{−} to R^{+}. These additional resources have no significant impact on the optimal lockdown strategy and death toll dynamics of the epidemic. The only impact concerns the socioeconomic cost of the pandemic, since identified recovered (and immune) individuals will be allowed to have a regular level of social interactions and won’t be quarantined.
On Figure B.10, we look towards the optimal detection strategy (with w_{immunity} = 1) whenever no lockdown policy is in place, while Figure B.11 focuses on the optimization on both lockdown intervention δ and immunity detection effort λ^{2}. As expected, additional optimization over λ^{2} has few impact on the sanitary burden, but can be used to lower the economic cost. The optimal detection strategy consists in identifying recovered individuals as soon as possible, and requires a large amount of effort in the first part of the epidemic phase. As a consequence, in the presence of immunity detection policy, the level of social interactions during the last phase of the epidemic is significantly higher, thus reducing the induced economic cost.
Finally, on Figure B.12, we optimize over all possible control levers: lockdown measure, virologic detection effort λ^{1} as well as the immunity detection one λ^{2}. We observe an optimal policy and induced epidemic dynamics rather similar to the one obtained in Figure B.8, whenever no immunity detection mean was available. The optimal λ^{2} helps in reducing the overall economic costs by detecting immune people.
4 Conclusion
The present work investigates, using an approach from Optimal Control theory, what could be optimal interventions for limiting the spread of a disease like the COVID19 epidemic. We consider first the case where only the lockdown lever is operated, as was done in France when massive testing was still not an option. Then, in the other scenarios, optimization is done on both lockdown and detection capacity of infected and immune individuals.
The notion of ‘best scenario’ is of course very much impacted by the choice of objective function that is made. The optimization takes into account the tradeoff between lowering the number of deaths and minimizing the economic and social costs. To our knowledge, taking into consideration these diverse aspects into the objective function has not been studied much in the recent literature (see e.g. [3]), although the literature on the topic is growing rapidly.
Another originality of our approach is the consideration of a strict constraint on the ICU capacity level dedicated to COVID19 patients. We observe that in the solutions, the lockdown and detection controls are set so that the occupation of ICU remains a high as possible under U_{max} (due to the economic and social costs). It is really this constraint that shapes the decisions leading to a curve that is ‘flattened’ sufficiently to ensure the sustainability of the health system. This is particularly visible in the scenario where U_{max} is increased of 50% within 2 months at the beginning of the epidemic: we see a relaxation of the lockdown that is permitted by the increased capacity of ICUs.
The best scenarios that we obtain are structured into 4 phases: (1) quick and strong measures to recover the control of the epidemic, (2) relaxation of the epidemic once under control to reach the fluxes imposed by the sustainability of the health care system, (3) the fluxes obtained after the period 2 are kept as long as possible to flatten the curve in order to avoid overwhelming the ICUs, (4) once a herd immunity is reached, the lockdown and testing controls can be lowered.
Without being an explicit target in the control problem, the evolution of the effective reproduction number ℜ_{t} is extremely interesting, and consistent with existing epidemiological literature. Again, without monitoring that indicator, we obtain in all the optimal strategies, ℜ_{t} becomes controlled below 1 after a few days (after the phase 1, and at the price of an important initial lockdown), and remains just below that critical value^{2}. It is allowed to reach and remain stable at level 1 in order to stabilize the ICU occupancy at their sustainable level U_{max}.
Of course, the model considered here is schematic because we wanted: (1) to put the methodology in evidence, (2) to explore easily the different scenarios. Generalizations can be carried in a quite straightforward manner. First, the major drawback of our approach is the necessity to fully observe some underlying latent variables and parameters of the epidemic dynamics, such as the prevalence of the virus. But we observed that the shape of the optimal control is quite robust to a variation of the initial prevalence , and more generally to variations upon lots of parameters in the model. An important issue about that COVID19 pandemic was the poor quality of data in the beginning, especially when trying to compare among countries to get robust estimators of various quantities regarding the disease. In our model, we tuned parameters using values from existing literature, and most of the results we obtained were actually extremely robust to moderate changes of those values.
Second, most individuals will not experience severe symptoms whenever contaminated by the virus. Such analysis calls for the use of differential quarantine strategies based on an estimation of the own risk of each individual. In our approach we use immunity detection procedure for a similar purpose, but less accurate estimates of the risk of each individuals are also available depending on their current and past health characteristics. As discussed in [31], the key ingredient for considering and using a differential lockdown strategy is a proper handling of the risk of bad specification of the group associated to some individuals. A dynamic optimisation over both the dynamic lockdown intervention together with the risk level associated to each misspecification is left for further research. Similarly the risk of wrong identification of infected or immune individuals using parameters λ^{1} and λ^{2} should be also taken into account in the design of a proper and efficient detection strategy. Also, in view of more practical applications for health policy guidance, a vector based version of SIDUHR^{+∕−} should probably be derived, as in [3] for instance. Categories can be splited depending on the age of individuals. Hence, S_{t} will then be , with children, adults and seniors, for instance. Most of the equations of the dynamics become multivariate, and parameters will be vectors, except β that becomes a WAIFW (or Who Acquires Infection From Whom) matrix, as in [80] or [59]. The main interest would be to have target and agespecific controls, with testing and lockdown that can be per age group. The difficult component in this multivariate extension is on the objective function, where the economic component should probably be more related to wellbeing cost of being quarantined, especially for people who do not belong to labor groups.
Finally let us notice that in this paper, as always in epidemic modelling, an important feature is the − ∕+ problem, with undetected individuals. The huge impact of infected detection resources λ^{1} in our numerical experiments emphasize the critical role of tracing and testing resources upon the dynamic of the epidemic. As suggested in [17, 27, 43, 49] (among many others), tracing is clearly a natural tool that should be optimized, to lower the socioeconomic cost of the disease without endangering the entire population (and increase the sanitary cost of the pandemic). This paves the way for further research topics.
Appendix A Parameters of the epidemic dynamics
First, let us notice, as mentioned in the main text, that there is an underlying individual based model for the system (2.1), in which the rates have a probabilistic interpretation, and can be interpreted as inverse of mean times (e.g. [9]), possibly with thinning when the corresponding events only occur for a fraction of the population. Based on these remarks, we can compute the rates in our model with the parameters provided in the French COVID19 literature (see [26, 73]).
A.1 Obtaining the transition rates from the sojourn times and transition probabilities
First, we introduce some variables on which the mechanistic construction of our rates will depend:

p_{a}: conditioned on being infected, the probability of having light symptoms or no symptoms;

p_{U}: conditioned on being mild/severely ill, the probability of needing hospitalization;

p_{H}: conditioned on needing hospitalization, the probability of needing ICU;

p_{d}: conditioned on being admitted to ICU, the probability of dying,
and for the durations:

: if asymptomatic, number of days until recovery;

: if symptomatic, number of days until recovery without hospital;

N_{IH}: (if severe symptomatic) number of days until hospitalization;

N_{HU}: if in H, number of days until ICU;

N_{HR}: if hospitalized but not in ICU, the number of days until recovery;

N_{UD} = 10: if in ICU, number of days until death;

N_{UR} = 20: if in ICU, number of days until recovery.
These variables can be measured quite easily on the data. Then, once these parameters are chosen, we can propose the following transition rates for our model:

γ_{IR}, the daily individual transition rate from I to R, is assumed to be of the form:
where on the right hand side, the first term is associated to mild symptomatic while the second term is associated to severe ones.

For the transition rates from U to R or D, the nonlinearity induced by the constraint on the ICU capacity U_{max} implies that it is easier to give these rates at the population level rather than at the individual level. For the transition from U to R, the individual rate is denoted by γ_{UR}(U) and the population rate is then γ_{UR}(U) ⋅ U. We assume here that:
This means that when the ICU capacity U_{max} is reached, new patients can not be taken in charge any more and since they correspond to individuals with severe forms of the disease, they can not recover until the ICU occupation is lowered.

As for the transition rate from U to R, the transition rate from U to D is nonlinear. We denote by γ_{UD}(U) the daily death rate of an individual needing ICU, and by γ_{UD}(U) × U the transition at the population level, that we assume is of the form: (A.1)
The hypothesis is that when ICU gets saturated, the patients who are not taken in charge experience an extra death rate.
A.2 Parameter values
For the infection parameter β, we fix the value of ℜ_{0} and then compute the corresponding value of β using (2.2). Reference [73] use ℜ_{0} = 3.3 (in absence of lockdown intervention and of testing policy), which we also recover from the early French dynamics ([82]) using the methods described in [77] (not shown).
The other transition rates appearing on Figure 1 can be obtained once the probabilities and durations detailed above are fixed. For this, we choose our parameters on the basis of the works by [26, 73] for the French COVID19 epidemics. We emphasize that despite of this choice, the methodology developed here is general and can apply to any country of disease.
The line(a) in Table A.1 can be obtained as follows: the model of [26] for IledeFrance is structured into 3 ageclasses: children (25%), adults (60%) and seniors (15%). Using the agespecific transition probabilities given in this reference and averaging over the age classes, we obtain the values announced:
However, itappears in our simulations that the probability p_{a} = 0.52 to be asymptomatic or only with light symptoms does not provide realistic orders for the proportion of dead individuals. In [73], they estimate the proportion of being hospitalized when infected at 2.6%. This proportion corresponds to (1 − p_{a}) p_{H} in our model, yielding a probability p_{a} = 0.85 when we use the estimate p_{H} = 0.17 of [26].
In France, at the beginning of the lockdown on March 15th 2020, there was roughly 12,000 beds in ICUs, which provides an estimate for the value of U_{max}. The additional mortality rate for individuals in ICUs when U_{t} > U_{max} is chosen 20 times more than the usual death rate in (A.1). This factor is calibrated from the observation by [73] that the mean time to death in hospitals is bimodal, with a group dying quickly after 0.67 days on average and a second group who die after longer time periods of mean 13.2 days. Extrapolating we assumed that ‘urgent’ cases die 20 times faster than ‘usual cases’.
For the value of , [33, Tab. 1) report a total number of infected of 3%, with a 95% credible interval [1.1%;7.8%], as of 28th March 2020 (which is two week after our t = 0). Because of the exponential increase, a crude estimate would be below 1% two weeks before. The number of reported cases in France on March 15 was around 5,000. According to [71], when Italy had 5,000 reported cases, an estimation of 60,000 cases was suggested, which would imply a close to 1%_{0}. In our benchmark scenario, we kept that 5%_{0} value for .
To check the validity of our parameters, we simulated our model in absence of intervention (δ = 0, λ^{1} ≡ 0 and λ^{2} ≡ 0) and compare the results with the prediction of [70]. Because in the latter reference they do not model the ICU saturation, we also removed the threshold U_{max} and the associated additional mortality (which amounts to considering U_{max} = +∞). Doing this we obtain similar results as those of [70].
Appendix B Sensitivities to the model specifications
In this appendix, we provide sensitivity graphs, where parameters are changed, to assess the robustness of our conclusions. Do not forget that numerical approximations errors can affect a too precise interpretation of the following numerical illustrations.
We consider the following benchmark scenario. The parameters for the dynamics are specified in Table 1. In the control problem, the benchmark situation is the one described in Section 3.2. To be precise, for the numerical implementation, we use the cost function is the one defined in (C.2) with the following choice of parameters: α =0, w_{sanitary} = 10^{5}, w_{econ} = 1, w_{prevalence} = 1, w_{immunity} = 1, w_{ICU} = 5.10^{4} and T = 700. Hence the optimal control problem for which we compute an approximate solution numerically is: (B.1)
with the set of measurable functions from [0, T] to [0, 1].
B.1 Impact of the initial virus prevalence
Figure B.1
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of . The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to a smaller weight for while the dotted line corresponds to a larger . 
B.2 Impact of the basic reproductive number ℜ_{0}
Figure B.2
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of ℜ_{0}. The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to a smaller value of ℜ_{0} (− 0.3) while the dotted line corresponds to a larger value of ℜ_{0} (+ 0.3). 
B.3 Impact of balancing sanitary and socioeconomic costs
Figure B.3
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of w_{sanitary}. The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to a smaller weight for w_{sanitary} (half) while the dotted line corresponds to a larger weight (twice). 
B.4 Impact of the vaccination arrival date anticipation α
Figure B.4
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of α, the discount rate that reflects the expected arrival time of a vaccine (and a cure). The plain line is the benchmark scenario discussed in Section 3.2; the dotted line corresponds to α = 1∕500, dashed line corresponds to α = 1∕250 and the mixed line corresponds to α = 1∕100. 
B.5 Impact of additional constant effort in virologic detection λ^{1}
Here, we add virologic effort for a constant value of λ^{1} and we compare the result of the optimization over δ for several values of λ^{1}. The optimalcontrol problem for which we compute an approximate solution numerically is: (B.2)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.5
Evolution of states with optimal control δ for three different values of λ^{1}, or (δ^{*}, λ^{1}, 0). The plain line is the benchmark scenario discussed in Section 3.2, with λ^{1} = 0; and then three scenarios, with λ^{1} = 1% (dashed line), λ^{1} = 5% (dotted line)and λ^{1} = 10% (mixed line). 
B.6 Impact of additional rising effort in virologic detection λ^{1}
Here, we add virologic effort for a value of λ^{1} which increases linearly in time, . The optimalcontrol problem for which we compute an approximate solution numerically is: (B.3)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.6
Evolution of states with optimal control δ for two different patterns for λ^{1}. The dotted line corresponds to the case λ^{1} = 1% and the mixed line is the increasing scenarios for (from 0% and 20% over two years). 
B.7 Optimizing over effort in virologic detection without lockdown intervention
Here, we compare the situation without intervention at all and the situation with intervention through λ^{1}. The optimal control problem for which we compute an approximate solution numerically is: (B.4)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.7
Evolution of states with optimal control (0, λ^{1}, 0), where the plain line corresponds to the benchmark scenario discussed in Section 3.1 (with no intervention δ = λ^{1} = λ^{2} = 0); and evolution with optimal control λ^{1} (when δ = λ^{2} = 0). 
B.8 Optimizing over both lockdown intervention δ and effort in virologic detection λ^{1}
Here, we compare the situation with optimal lockdown policy versus the situation in which we optimize over both lockdown and virologic effort. The latter corresponds to the following optimal control problem, for which we compute an approximate solution numerically: (B.5)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.8
Evolution of states with optimal controls (δ, λ^{1}, 0), where the plain line corresponds to the benchmark scenario (δ, 0, 0) discussed in Section 3.2 (plain line); and evolution with joint optimal controls (δ, λ^{1}, 0) (dashed line). 
B.9 Impact of additional effort in immunity detection λ^{2}
Here, we addvirologic effort for a constant value of λ^{2} and we compare the result of the optimization over δ for several values of λ^{2}. The optimalcontrol problem for which we compute an approximate solution numerically is: (B.6)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.9
Evolution of states with optimal control δ for three different values of λ^{2} given a fixed level of λ^{1} (here 1%), with λ^{2}= 0 (dotted line), λ^{2} = 1%_{0} (dashed line) and λ^{2} = 1% (mixed line),i.e. (δ^{*}, λ_{1}, λ_{2}). 
B.10 Optimizing over effort in immunity detection λ^{2} without lockdown intervention
Here, we compare the situation without intervention at all and the situation with intervention through λ^{2}. The optimal control problem for which we compute an approximate solution numerically is: (B.7)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.10
Evolution of states with optimal controls (0, 0, λ^{2}), where the plain line corresponds to the benchmark scenario discussed in Section 3.1 (with no intervention δ = λ^{1} = λ^{2} = 0); and evolution with optimal control λ^{2} (when δ = λ^{1} = 0). 
B.11 Optimizing over both lockdown intervention δ and effort in immunity detection λ^{2}
Here, we compare the situation with optimal lockdown policy versus the situation in which we optimize over both lockdown and immunity detection. The latter corresponds to the following optimal control problem, for which we compute an approximate solution numerically: (B.8)
with the set of measurable functions from [0, T] to [0, 1].
Figure B.11
Evolution of states with optimal controls (δ, 0, λ^{2}), where the plain line corresponds to the benchmark scenario (δ, 0, 0) discussed in Section 3.2 (plain line); and evolution with joint optimal controls (δ, 0, λ^{2}) (dashed line). 
B.12 Optimizing over lockdown intervention δ, and efforts in virologic detection λ^{1} and immunity detection λ^{2}
Here, we compare the optimization on δ only versus the optimization (C.1) over (δ, λ^{1}, λ^{2}).
Figure B.12
Evolution of states with optimal controls (δ, λ^{1}, λ^{2}), where the plain line corresponds to the benchmark scenario (δ, 0, 0) discussed in Section 3.2 (plain line); and evolution with joint optimal controls (δ, λ^{1}, λ^{2}) (dashed line). 
Appendix C Numerical resolution of the optimal control problem
We discuss in this section the numerical resolution of the optimal control problem (2.12) of interest.
The original problem has infinite horizon, but for numerical purpose, we focus on the equivalent version over a sufficiently large time horizon T. The proximity between both optimal control problems is well established in the literature, see e.g. Section 4.1 in [29], while we observe empirically very few impact by the choice of maturity T, as soon as T is chosen large enough in order for the epidemic phase to be over (i.e. such that the proportion of infected people is small enough at time T). With respect to the random time of vaccination τ, such approximation corresponds to optimizing over the time interval [0, τ∧ T] instead of [0, τ].
Similarly, the problem (2.12) is an optimal control problem with state constraint on the ICU capacity, and we choose to replace for ease of numerical tractability the state constraint by a penalization cost with a sufficiently large weight w_{ICU} > 0. Hence, we focus numerically on an approximate version of the original optimal control, which is given by: (C.1)
where the approximate objective function is given by (C.2)
C.1 Optimality condition by Pontryagin’s maximum principle
To alleviate the notation and use more standard notations from the literature, the problem can be recast in the following form, where a denotes the control and X the state: Minimize
where we denoted ρ_{UR}(u) = γ_{UR}(u)u and ρ_{UD}(u) = γ_{UD}(u)u in order to alleviate the notation.
Let H be the Hamiltonian defined by
Pontryagin’s maximum principle leads, at least informally, to the following optimality condition: if is optimal, then for every t ∈ [0, T],
where solve the forwardbackward ODE system: (C.3)
The backward equation can be written as follows (dropping the hats to alleviate the notation):
C.2 Numerical algorithm
The numerical method is an iterative procedure to compute the optimal control . Starting from an initial guess a^{(0)}, at iteration k ≥ 0, in order to compute a^{(k+1)}, we first compute the solution (X^{(k)}, Y^{(k)}) to the forwardbackward ODE system corresponding to (C.3) but with control a^{(k)} instead of â, i.e.,
where θ^{(k)} > 0 is a step size (which may depend on the iteration in order to adjust the convergence rate).
Finally, we define, for each component i of the control,
where π_{[0,1]} denotes the projection on the interval [0, 1].
If the optimization is to be performed over only one component of the control (e.g. over δ only), then the updates are done only for this component instead of the whole vector of controls.
References
 A. Abakuks, An optimal isolation policy for an epidemic. J. Appl. Probab. 10 (1973) 247–62. [CrossRef] [Google Scholar]
 A. Abakuks, Optimal immunisation policies for epidemics. Adv. Appl. Probab. 6 (1974) 494–511. [CrossRef] [Google Scholar]
 D. Acemoglu, V. Chernozhukov, I. Werning and M. Whinston, A multirisk SIR model with optimally targeted lockdown. NBER Working Paper 27102 (2020) 1–38. [Google Scholar]
 F. Agusto and A. Adekunle, Optimal control of a twostrain tuberculosisHIV/AIDS coinfection model. Biosystems 119 (2014) 20–44. [CrossRef] [PubMed] [Google Scholar]
 J.A. AlTawfiq, Asymptomatic coronavirus infection: MERSCoV and SARSCoV2 (COVID19). Travel Med. Infect. Dis. 101608 (2020). [CrossRef] [PubMed] [Google Scholar]
 R. Aldridge, D. Lewer, S. Beale, A. Johnson, M. Zambon, A. Hayward, E. Fragaszy and n. null, Seasonality and immunity to laboratoryconfirmed seasonal coronaviruses (HCoVNL63, HCoVOC43, and HCoV229E): results from the flu watch cohort study. Wellcome Open Res. 5 (2020) 52. [CrossRef] [PubMed] [Google Scholar]
 F. Alvarez, D. Argente and F. Lippe, A simple planning problem for COVID19 lockdown. Natl. Bureau Eco. Res. 26981 (2020) 1–35. [Google Scholar]
 R. Anderson and R. May, Infectious Diseases of Humans: Dynamics and Control. Oxford University Press, Oxford (1991). [Google Scholar]
 F. Ball, T. Britton, C. Larédo, E. Pardoux, D. Sirl and V. Tran, Stochastic Epidemic Models with Inference, edited by T. Britton and E. Pardoux, Lecture Notes in Mathematics. Springer (2019). [Google Scholar]
 E. Barclay, The US doesn’t just need to flatten the curve. it needs to “raise the line”. Available from https://www.vox.com/2020/4/7/21201260/coronavirususachartmaskshortageventilatorsflattenthecurve (2020). [Google Scholar]
 H. Behncke, Optimal control of deterministic epidemics. Optim. Control Appl. Methods 21 (2000) 269–285. [Google Scholar]
 D.W. Berger, K.F. Herkenhoff and S. Mongey, An SEIR infectious disease model with testing and conditional quarantine. Workingpaper 26901. National Bureau of Economic Research, Inc. (2020). [Google Scholar]
 J. Bernstein, A.W. Richter and N. Throckmorton, Covid19: A view from the labor market. Federal Reserve Bank of Dallas Working Paper 2010 (2020). [Google Scholar]
 L. Bobisud, Optimal control of a deterministic epidemic. Math. Biosci. 35 (1977) 165–174. [CrossRef] [Google Scholar]
 M. Chinazzi, J.T. Davis, M. Ajelli, C. Gioannini, M. Litvinova, S. Merler, A. Pastore y Piontti, K. Mu, L. Rossi, K. Sun, C. Viboud, X. Xiong, H. Yu, M.E. Halloran, I.M. Longini and A. Vespignani, The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID19) outbreak. Science 368 (2020) 395–400. [CrossRef] [Google Scholar]
 P. Choe, R. Perera, W. Park, K. Song, J. Bang, E. Kim and M. Oh, MERSCoV antibody responses 1 year after symptom onset, South Korea, 2015. Emerg. Infect. Dis. 23 (2017) 1079–1084. [CrossRef] [Google Scholar]
 S. Clémençon, V. Tran and H.D. Arazoza, A stochastic SIR model with contacttracing: large population limits and statistical inference. J. Biol. Dyn. 2 (2008) 391–414. [Google Scholar]
 C.W. Cobb and P.H. Douglas, A theory of production. Am. Eco. Rev. 18 (1928) 139–165. [Google Scholar]
 J. Cohen, Vaccine designers take first shots at COVID19. Science 368 (2020) 14–16. [CrossRef] [Google Scholar]
 J. Cohen and K. Kupferschmidt, Countries test tactics in ‘war’ against COVID19. Science 367 (2020) 1287–1288. [CrossRef] [Google Scholar]
 A.R. da Cruz, R.T.N. Cardoso and R.H.C. Takahashi, Multiobjective dynamic optimization of vaccination campaigns using convex quadratic approximation local search, in Evolutionary MultiCriterion Optimization, edited by R.H.C. Takahashi, K. Deb, E.F. Wanner and S. Greco. Springer, Berlin (2011) 404–417. [Google Scholar]
 M. Day, COVID19: identifying and isolating asymptomatic people helped eliminate virus in Italian village. Br. Med. J. 368 (2020). [Google Scholar]
 O. Diekmann, J. Heesterbeek and J. Metz, On the definition and the computation of the basic reproduction ratio r_{0} in models for infectious diseases in heterogeneous populations. J. Math. Biol. 28 (1990) 365–382. [Google Scholar]
 O. Diekmann, H. Heesterbeek and T. Britton, Mathematical Tools for Understanding Infectious Disease Dynamics. Princeton Series in Theoretical and Computational Biology. Princeton University Press, New Jersey (2012). [Google Scholar]
 R. DjidjouDemasse, Y. Michalakis, M. Choisy, M.T. Sofonea and S. Alizon, Optimal COVID19 epidemic control until vaccine deployment. Preprint medRxiv 20049189v1 (2020). [Google Scholar]
 L.D. Domenico, G. Pullano, C. Sabbatini, P.Y. Boëlle and V. Colizza, Expected impact of lockdown in IledeFrance and possible exit strategies. Preprint medrxiv 20063933v1 (2020). [Google Scholar]
 K. Eames and M. Keeling, Contact tracing and disease control. Proc. R. Soc. London B 270 (2003) 2565–2571. [CrossRef] [Google Scholar]
 M.S. Eichenbaum, S. Rebelo and M. Trabandt, The macroeconomics of epidemics. Working Paper 26882. National Bureau of Economic Research (2020). [Google Scholar]
 R. Elie, Finite time Merton strategy under drawdown constraint: a viscosity solution approach. Appl. Math. Optim. 58 (2008) 411–431. [CrossRef] [Google Scholar]
 R. Elie, E. Hubert and G. Turinici, Contact rate epidemic control of COVID19: an equilibrium view. Preprint arXiv:2004.08221 (2020). [Google Scholar]
 T. Evgeniou, M. Fekom, A. Ovchinnikov, R. Porcher, C. Pouchol and N. Vayatis, Epidemic models for personalised COVID19 isolation and exit policies using clinical risk predictions. Preprint medRxiv 20074054v1 (2020). [Google Scholar]
 N. Ferguson, D. Laydon, G. NedjatiGilani, N. Imai, K. Ainslie, M. Baguelin, S. Bhatia, A. Boonyasiri, Z. Cucunubá, G. CuomoDannenburg, A. Dighe, I. Dorigatti, H. Fu, K. Gaythorpe, W. Green, A. Hamlet, W. Hinsley, L.C. Okell, S. van Elsland and A.C. Ghani, Impact of nonpharmaceutical interventions (NPIs) to reduce COVID19 mortality and healthcare demand. Imperial College COVID19 Response Team 9 (2020). [Google Scholar]
 S. Flaxman, S. Mishra, A. Gandy, J. Unwin, H. Coupland, a.Z. Thomas A Mellan, T. Berah, A. Ghani, C.A. Donnelly, S. Riley, L.C. Okell, M.A.C. Vollmer, N.M. Ferguson and S. Bhatt, Estimating the number of infections and the impact of nonpharmaceutical interventions on COVID19 in 11 European countries. Imperial College COVID19 Response Team 13 (2020). [Google Scholar]
 C. Gelardi, Colonialism made puerto rico vulnerable to coronavirus catastrophe. Available from: https://www.thenation.com/article/politics/puertoricocoronavirus/ (2020). [Google Scholar]
 K. Gostic, A.C. Gomez, R.O. Mummah, A.J. Kucharski and J.O. LloydSmith, Estimated effectiveness of symptom and risk screening to prevent the spread of COVID19. eLife 9 (2020) e55570. [Google Scholar]
 D. Greenhalg, Some results on optimal control applied to epidemics. Math. Biosci. 88 (1988) 125–158. [CrossRef] [Google Scholar]
 S.K. Gudi, K. Undela, R. Venkataraman, U.V. Mateti, M. Chhabra, S. Nyamagoud and K.K. Tiwari, Knowledge and beliefs towards universal safety precautions to flatten the curve during novel coronavirus disease (nCOVID19) pandemic among general public in India: Explorations from a national perspective. Preprint medRxiv 20047126v1 (2020). [Google Scholar]
 V. Guerrieri, G. Lorenzoni, L. Straub and I. Werning, Macroeconomic implications of COVID19: Can negative supply shocks cause demand shortages? Working Paper 26918, National Bureau of Economic Research (2020). [Google Scholar]
 E. Hansen, and T. Day, Optimal control of epidemics with limited resources. J. Math. Biol. 62 (2011) 423–451. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 X. He, E. Lau, P. Wu, X. Deng, J. Wang and X. Hao, Temporal dynamics in viral shedding and transmissibility of COVID19. Nat. Med. 26 (2020) 672–5. [CrossRef] [PubMed] [Google Scholar]
 J. Hellewell, S. Abbott, A. Gimma, N. Bosse, C. Jarvis, T. Russell, J. Munday, A. Kucharski, J. Edmunds, C.C. W. Group, S. Funk and R. Eggo, Feasibility of controlling COVID19 outbreaks by isolation of cases and contacts. Lancet Glob. Health 8 (2020) 488–496. [CrossRef] [Google Scholar]
 J. Hellewell, S. Abbott, A. Gimma, N.I. Bosse, C.I. Jarvis, T.W. Russell, J.D. Munday, A.J. Kucharski, W.J. Edmunds, F. Sun, S. Flasche, B.J. Quilty, N. Davies, Y. Liu, S. Clifford, P. Klepac, M. Jit, C. Diamond, H. Gibbs, K. [van Zandvoort], S. Funk, and R.M. Eggo, Feasibility of controlling COVID19 outbreaks by isolation of cases and contacts. Lancet Glob. Health 8 (2020) e488–e496. [CrossRef] [Google Scholar]
 T. House and M. Keeling, The impact of contact tracing in clustered populations. PLoS Comput. Biol. 6 (2010) e1000721. [CrossRef] [PubMed] [Google Scholar]
 L. Huang, Y. Shi, B. Gong, L. Jiang, X. Liu, J. Yang, J. Tang, C. You, Q. Jiang, B. Long, T. Zeng, M. Luo, F. Zeng, F. Zeng, S. Wang, X. Yang and Z. Yang, Blood single cell immune profiling reveals the interferonMAPK pathway mediated adaptive immune response for COVID19. Preprint medRxiv 20033472v1 (2020). [Google Scholar]
 D. Iacoviello and G. Liuzzi, Optimal control for SIR epidemic model: A two treatments strategy. 2008 Mediterranean Conference on Control and Automation  Conference Proceedings, MED’08 (2008) 842–847. [CrossRef] [Google Scholar]
 F. Jiang, L. Deng and L. Zhang, Review of the Clinical Characteristics of Coronavirus Disease 2019 (COVID19). J. Gen. Inter. Med. 35 (2020) 1545–9. [CrossRef] [Google Scholar]
 W. Kermack and A. McKendrick, A contribution to the mathematical theory of epidemics. Proc. R. Soc. London A 115 (1927) 700–721. [Google Scholar]
 S. Kim, J. Lee and E. Jung, Mathematical model of transmission dynamics and optimal control strategiesfor 2009 A/H1N1 influenza in the Republic of Korea. J. Theor. Biol. 412 (2017) 74–85. [CrossRef] [PubMed] [Google Scholar]
 I. Kiss, D. Green and R. Kao, Infectious disease control using contact tracing in random and scalefree networks. J. R. Soc. Interface 3 (2013) 55–62. [CrossRef] [PubMed] [Google Scholar]
 S.M. Kissler, C. Tedijanto, M. Lipsitch and Y. Grad, Social distancing strategies for curbing the COVID19 epidemic. Preprint medRxiv 20041079v1 (2020). [Google Scholar]
 C.C. Ku, T.C. Ng and H.H. Lin, Epidemiological benchmarks of the COVID19 outbreak control in China after Wuhan’s lockdown: A modelling study with an empirical approach. SSRN Electron. J. (2020) 3543589. [Google Scholar]
 A. Kucharski, P. Klepac, A. Conlan, S. Kissler, M. Tang, H. Fry, J. Gog, J. Edmunds and C.C. working group, Effectiveness of isolation, testing, contact tracing and physical distancing on reducing transmission of SARSCoV2 in different settings. Preprint medRxiv 20077024v1 (2020). [Google Scholar]
 A. Kucharski, T. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk and R. Eggo, Early dynamics of transmission and control of COVID19: a mathematical modelling study. Lancet Infect. Dis. 20 (2020) 553–58. [CrossRef] [PubMed] [Google Scholar]
 A. Kumar and P.K. Srivastava, Vaccination and treatment as control interventions in an infectious disease model with their cost optimization. Commun. Nonlinear Sci. Numer. Simul. 44 (2017) 334–343. [CrossRef] [Google Scholar]
 C. Lagorio, M. Dickison, F. Vazquez, L.A. Braunstein, P.A. Macri, M.V. Migueles, S. Havlin and H.E. Stanley, Quarantinegenerated phase transition in epidemic spreading. Phys. Rev. E 83 (2011) 026102. [CrossRef] [Google Scholar]
 S. Lai, N.W. Ruktanonchai, L. Zhou, O. Prosper, W. Luo, J.R. Floyd, A. Wesolowski, M. Santillana, C. Zhang, X. Du, H. Yu and A.J. Tatem, Effect of nonpharmaceutical interventions for containing the COVID19 outbreak in China. Preprint medRxiv 20029843v1 (2020). [Google Scholar]
 W. Liu, Q. Zhang, J. Chen, R. Xiang, H. Song, S. Shu, L. Chen, L. Liang, J. Zhou, L. You, P. Wu, B. Zhang, Y. Lu, L. Xia, L. Huang, Y. Yang, F. Liu, M.G. Semple, B.J. Cowling, K. Lan, Z. Sun, H. Yu, and Y. Liu, Detection of Covid19 in children in early January 2020 in Wuhan, China. New Engl. J. Med. 382 (2020) 1370–1371. [CrossRef] [Google Scholar]
 P. Magal and G. Webb, Predicting the number of reported and unreported cases for the COVID19 epidemic in South Korea, Italy, France and Germany. Preprint medRxiv 20040154v1 (2020). [Google Scholar]
 J. Mossong, N. Hens, M. Jit, P. Beutels, K. Auranen and R. Mikolajczyk, Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 5 (2008) e74. [Google Scholar]
 H. Nishiura, T. Kobayashi, A. Suzuki, S.M. Jung, K. Hayashi, R. Kinoshita, Y. Yang, B. Yuan, A. Akhmetzhanov, N. Linton and T. Miyama, Estimation of the asymptomatic ratio of novel coronavirus infections (COVID19). Int. J. Infect. Dis. 94 (2020) 154–155. [Google Scholar]
 A. Oliver, A normative perspective on discounting health outcomes. J. Health Serv. Res. Policy 18 (2013) 186–189. [CrossRef] [PubMed] [Google Scholar]
 P.L. Ooi, S. Lim and S.K. Chew, Use of quarantine in the control of SARS in Singapore. Am. J. Infect. Control 33 (2005) 252–257. [CrossRef] [Google Scholar]
 M.G. Pedersen and M. Meneghini, A simple method to quantify countryspecific effects of COVID19 containment measures. Preprint medRxiv 20057075v1 (2020). [Google Scholar]
 J.Peto, COVID19 mass testing facilities could end the epidemic rapidly. BMJ 368 (2020). [Google Scholar]
 F. Piguillem and L. Shi, Optimal COVID19 quarantine and testing policies. EIEF Working Papers Series (2004) 2020. [Google Scholar]
 L. Pontryagin, G. Boltyanskii, R. Gamkrelidze and E. Mishchenko, Mathematical Theory of Optimal Processes, CRC Press, New York (1964). [Google Scholar]
 N. Qualls, A. Levitt and N. Kanade, Community mitigation guidelines to prevent pandemic influenza  United States. MMWR 66 (2017) 1–34. [Google Scholar]
 M.L. Ranney, V. Griffeth and A.K. Jha, Critical supply shortages  the need for ventilators and personal protective equipment during the COVID19 pandemic. New Engl. J. Med. 382 (2020) e41. [CrossRef] [Google Scholar]
 L. Roques, E. Klein, J. Papaix, A. Sar and S. Soubeyrand, Effect of a onemonth lockdown on the epidemic dynamics of COVID19 in France. Preprint medRxiv 20074054v1 (2020). [Google Scholar]
 J. Roux, C. Massonnaud and P. Crépey, COVID19: Onemonth impact of the French lockdown on the epidemic burden. Available from: https://www.ehesp.fr/wpcontent/uploads/2020/04/ImpactConfinementEHESP20200322v11.pdf (2020). [Google Scholar]
 A. Sachdeva and A. Sheth, COVID19, panic now!! a call to action because the numbers are deceptive. SSRN 3563419 (2020). [Google Scholar]
 M. Salathé, C.L. Althaus, R. Neher, S. Stringhini, E. Hodcroft, J. Fellay, M. Zwahlen, G. Senti, M. Battegay, A. WilderSmith, I. Eckerle, M. Egger and N. Low, COVID19 epidemic in Switzerland: on the importance of testing, contact tracing and isolation. Swiss medical weekly 150 (2020) w20225. [PubMed] [Google Scholar]
 H. Salje, C.T. Kiem, N. Lefrancq, N. Courtejoie, P. Bosetti, J. Paireau, A. Andronico, N. Hoze, J. Richet, C.L. Dubost, Y.L. Strat, J. Lessler, D.L. Bruhl, A. Fontanet, L. Opatowski, P.Y. Boelle and S. Cauchemez, Estimating the burden of SARSCoV2 in France. Available from: https://halpasteur.archivesouvertes.fr/pasteur02548181 (2020). [Google Scholar]
 S.P. Sethi, Optimal quarantine programmes for controlling an epidemic spread. J. Operat. Res. Soc. 29 (1978) 265–268. [CrossRef] [Google Scholar]
 O. Sharomi and T. Malik, Optimal control in epidemiology. Ann. Operat. Res. 251 (2017) 55–71. [CrossRef] [Google Scholar]
 E. Tognotti, Lessons from the history of quarantine, from plague to influenza A. Emerg. Infect. Dis. 19 (2013) 254–259. [CrossRef] [PubMed] [Google Scholar]
 P. Trapman, F. Ball, J.S. Dhersin, V. Tran, J. Wallinga and T. Britton, Inferring R0 in emerging epidemicsthe effect of common population structure is small. J. R. Soc. Interf. 13 (2016) 20160288. [CrossRef] [Google Scholar]
 M. van der Pol and J. Cairns, A comparison of the discounted utility model and hyperbolic discounting models in the case of social and private intertemporal preferences for health. J. Eco. Behav. Organ. 49 (2002) 79–96. [CrossRef] [Google Scholar]
 E. Verriest, F. Delmotte and M. Egerstedt, Control of epidemics by vaccination, Vol. 2 of Proceedings of the 2005, American Control Conference (2005) 985–990. [CrossRef] [Google Scholar]
 J. Wallinga, P. Teunis and M. Kretzschmar, Using Data on Social Contacts to Estimate Agespecific Transmission Parameters for Respiratoryspread Infectious Agents. Am. J. Epidemiol. 164 (2006) 936–944. [CrossRef] [Google Scholar]
 S. Wong, A. Vaughan, C. QuiltyHarper and L. Liverpool, Covid19 news: Us not involved in global WHO plan to tackle pandemic. New Scientist April 24, 2020. [Google Scholar]
 World Health Organization, Coronavirus disease 2019 (COVID19). Available from: https://www.who.int/emergencies/diseases/novelcoronavirus2019 (2020). [Google Scholar]
 X. Yan and Y. Zou, Optimal and suboptimal quarantine and isolation control in SARS epidemics. Math. Comput. Model. 47 (2008) 235–45. [CrossRef] [PubMed] [Google Scholar]
 X. Zhou, Y. Li, T. Li and W. Zhang, Followup of the asymptomatic patients with SARSCoV2 infection. Clin. Microbiol. Infect. 26 (2020) P957–959. [CrossRef] [PubMed] [Google Scholar]
In all other figures of this article, but that one, we consider a 700 day horizon. Here, the epidemic phase is over within two months so we shorten the time frame on Figure 2 for ease of readability.
All Tables
All Figures
Figure 1
The SIDUHR^{+∕−} model. 

In the text 
Figure 2
Evolution of states without control, over 200 days. 

In the text 
Figure 3
Evolution of state proportions – S_{t} (susceptible), I_{t} (infected, detected or not, ), R_{t} (recovered and immune, detected or not, ), U_{t} (hospitalized in ICU), D_{t} (dead), ℜ_{t} (dynamic version of the reproduction number), W_{t} (labor force available), Q_{t} (susceptible to be quanrantines, ) and δ_{t} (lockdown intensity) – with optimal control δ (plain line) and without control (dashed line). States vary here from t = 0 (beginning of the pandemic, or at least of possible measures, with δ ≥ 0) to t = 700, and computations are based on a dt = 1∕5 days. 

In the text 
Figure 4
Phase diagram S_{t} (susceptible)versus I_{t} (infected, detected or not, ), together with the herd immunity threshold S^{*} and the levelI^{*} of infected corresponding to U = U_{max}. 

In the text 
Figure 5
Evolution of states with optimal control δ when ‘raising the line’ (increase of U_{max} by 50%). The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to the scenario with a capacity of U_{max} which increases linearly from time t = 0 until t = 60 days up to 1.5 times the initial capacity and then stays at this level. 

In the text 
Figure B.1
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of . The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to a smaller weight for while the dotted line corresponds to a larger . 

In the text 
Figure B.2
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of ℜ_{0}. The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to a smaller value of ℜ_{0} (− 0.3) while the dotted line corresponds to a larger value of ℜ_{0} (+ 0.3). 

In the text 
Figure B.3
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of w_{sanitary}. The plain line is the benchmark scenario discussed in Section 3.2; the dashed line corresponds to a smaller weight for w_{sanitary} (half) while the dotted line corresponds to a larger weight (twice). 

In the text 
Figure B.4
Evolution of states with optimal controls (δ, 0, 0) – without any testing – for three different values of α, the discount rate that reflects the expected arrival time of a vaccine (and a cure). The plain line is the benchmark scenario discussed in Section 3.2; the dotted line corresponds to α = 1∕500, dashed line corresponds to α = 1∕250 and the mixed line corresponds to α = 1∕100. 

In the text 
Figure B.5
Evolution of states with optimal control δ for three different values of λ^{1}, or (δ^{*}, λ^{1}, 0). The plain line is the benchmark scenario discussed in Section 3.2, with λ^{1} = 0; and then three scenarios, with λ^{1} = 1% (dashed line), λ^{1} = 5% (dotted line)and λ^{1} = 10% (mixed line). 

In the text 
Figure B.6
Evolution of states with optimal control δ for two different patterns for λ^{1}. The dotted line corresponds to the case λ^{1} = 1% and the mixed line is the increasing scenarios for (from 0% and 20% over two years). 

In the text 
Figure B.7
Evolution of states with optimal control (0, λ^{1}, 0), where the plain line corresponds to the benchmark scenario discussed in Section 3.1 (with no intervention δ = λ^{1} = λ^{2} = 0); and evolution with optimal control λ^{1} (when δ = λ^{2} = 0). 

In the text 
Figure B.8
Evolution of states with optimal controls (δ, λ^{1}, 0), where the plain line corresponds to the benchmark scenario (δ, 0, 0) discussed in Section 3.2 (plain line); and evolution with joint optimal controls (δ, λ^{1}, 0) (dashed line). 

In the text 
Figure B.9
Evolution of states with optimal control δ for three different values of λ^{2} given a fixed level of λ^{1} (here 1%), with λ^{2}= 0 (dotted line), λ^{2} = 1%_{0} (dashed line) and λ^{2} = 1% (mixed line),i.e. (δ^{*}, λ_{1}, λ_{2}). 

In the text 
Figure B.10
Evolution of states with optimal controls (0, 0, λ^{2}), where the plain line corresponds to the benchmark scenario discussed in Section 3.1 (with no intervention δ = λ^{1} = λ^{2} = 0); and evolution with optimal control λ^{2} (when δ = λ^{1} = 0). 

In the text 
Figure B.11
Evolution of states with optimal controls (δ, 0, λ^{2}), where the plain line corresponds to the benchmark scenario (δ, 0, 0) discussed in Section 3.2 (plain line); and evolution with joint optimal controls (δ, 0, λ^{2}) (dashed line). 

In the text 
Figure B.12
Evolution of states with optimal controls (δ, λ^{1}, λ^{2}), where the plain line corresponds to the benchmark scenario (δ, 0, 0) discussed in Section 3.2 (plain line); and evolution with joint optimal controls (δ, λ^{1}, λ^{2}) (dashed line). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.