key: cord-0727374-a7muub05 authors: Gondim, João A.M.; Machado, Larissa title: Optimal quarantine strategies for the COVID-19 pandemic in a population with a discrete age structure date: 2020-08-13 journal: Chaos, Solitons & Fractals DOI: 10.1016/j.chaos.2020.110166 sha: 73178d4c12c75a840d0d0375a5a034296e078f28 doc_id: 727374 cord_uid: a7muub05 Abstract The goal of this work is to study the optimal controls for the COVID-19 epidemic in Brazil. We consider an age-structured SEIRQ model with quarantine compartment, where the controls are the quarantine entrance parameters. We then compare the optimal controls for different quarantine lengths and distributions of the total control cost by assessing their respective reductions in deaths in comparison to the same period without quarantine. The best strategy provides a calendar of when to relax the isolation measures for each age group. Finally, we analyse how a delay in the beginning of the quarantine affects this calendar by changing the initial conditions. At the end of 2019 a novel coronavirus emerged in the city of Wuhan, China. In January 2020 the disease was given the name COVID-19 and, by mid February, China already faced over 60 thousand cases [1] . Many scientists then began to model disease to forecast its worldwide impact ( [2] , [3] ), even influencing policies of many governments. As the disease spread across Europe and the United States, some countries were forced to implement quarantines or even lockdowns to mitigate the harms and were able to do the so called "flattening of the curve", i.e., to postpone and dampen the maximum number of active cases. By mid May, Brazil is entering this stage, with some states declaring stricter quarantine policies. Optimal control theory has been applied to general epidemic models ( [4] , [5] ) as well as to specific disease models such as HIV ( [6] , [7] , [8] ), tuberculosis ([9] , [10] ) and influenza ( [11] ). Recently, a few works applying optimal control theory to the COVID-19 pandemic, such as [12] and [13] , have appeared. This paper focuses on a SEIR model with quarantine as was proposed in [14] , but dividing the population in age groups as in [15] . This is particularly important since COVID-19 has worse consequences on the elderly than it does on younger people. Our goal is to calculate the optimal quarantine strategies numerically for different choices of parameters in the model, which reflect the decisions governments must make when implementing these policies, such as evaluating the economical costs of the quarantine for each of the age groups and when to start implementing the measures. Then, we compare the controls by looking at how they reduce deaths in comparison to the same period without quarantine. The best strategy gave us a calendar of when to relax the measures in each of the age groups. Our model consists of a classical SEIR model with a quarantined class. Besides, we assume that the population has an age structure (see [16] , [17] for models with a continuous age structure and [18], [19] for models with a discrete one). There are three age groups, described in Table 1 . and Q i (t) be the number of susceptible, exposed, infected, recovered and quarantined individuals in each age group at time t ≥ 0, respectively. We assume that the total population All parameters are nonnegative. β ij is the transmission coefficient from age group i to age group j. Typically, it will be assumed that β ij = β ji for all i, j. σ i and γ i are the latency and recovery periods, respectively, for age group i. λ is the exit rate from the quarantine. Our controls are the u i (t), which denote the fraction of susceptible individuals in each age group that are put into quarantine per unit time at t. As such, they satisfy, a priori, However, it is unrealistic to expect an entire population to stay under quarantine for a long time. There are essential workers such as healthcare professionals and police officers that cannot stay at home during these times. As most of these workers are in age group 2, we suppose that all of age groups 1 and 3 can be quarantined (for age group 3, indeed, this is especially important since they are a risk group for the COVID-19 pandemic). Thus, we shall loosen (2) by considering, for example, Evaluating u max is one of the tasks of each government's authorities. In this paper, we fix this parameter at u max = 0.9. This means that Let be the total population of age group i. Adding the equations in the system above, we see that N i (t) is also constant for i ∈ {1, 2, 3}. Since R i (t) only appears in the other equations as a part of N i (t), we substitute the equations for R i i(t) by N i (t). Hence, we may also consider the system of Parameter values are taken from [15] and are listed in Table 2 . The data fitting was performed using an adaptation of a least-squares algorithm from [20] . To see how the numbers of infections and recoveries are distributed in the three age groups, we refer to the data available in [21] , shown in Table 3 . For simplification, we suppose that the difference between the number of cases and the number of deaths represents the number of recoveries. This is not necessarily correct, because some of the patients that we considered as recovered might still carry the disease, but we use this approach due to the scarcity of information regarding recoveries we currently have. The respective distributions are shown in Table 4 . According to [1] , Brazil had as of May 13, 2020 a total of 97,575 active COVID-19 cases. Even though there seems to be a large underreporting in the country [22] , this number will be considered as the total number of infected individuals nonetheless. To estimate the number of exposed cases, we look at data from May 8, 2020, since the mean incubation period of the disease is thought to be around 5 days [23] . Once again according to [1] , at this time Brazil had 76,603 active cases, so this gives us an estimation of 20,972 exposed cases. We also suppose that these cases follow the age distributions of cases from Table 4 . Finally, as of May 8, there were 65,124 recovered cases in Brazil Table 2 : Parameter values (data from [15] ). Value Parameter Value Therefore, our initial time will consist of data from Brazil as of May 8, 2020. The total population is assumed to be 200 million, divided into 40% young people, 50% adults and 10% elderly. We also assume that there are no quarantined individuals when the simulation starts. Since the numbers of exposed, infected and recovered are very small in comparison to the total population, we assume that the initial number of susceptible individuals is equal to the total population of the respective age group. The initial conditions of all variables, rounded to the nearest integers, are listed in Table 5. 6 3. The optimization problem Using system (4), we consider the functional to be minimized as Quarantined 0 0 0 In the formula above, T is the quarantine duration and the parameters B i are the costs of the control. We assume that B i > 0 for i ∈ {1, 2, 3} and that where B ∈ R is the total control cost. Sufficient conditions for the existence of the optimal controls follow from standard results from optimal control theory. For instance, we can use Theorem 2.1 in [24] to show that the optimal control exists. Pontryagin's maximum principle ( [25] , [26] ) establish that optimal controls are solutions of the Hamiltonian system with Hamiltonian function where λ S i , λ E i , λ I i , λ Q i and λ N i are the adjoint variables. These variables must 8 satisfy the adjoint equations where i ∈ {1, 2, 3} and C ∈ {S, E, I, Q, N }. The adjoint system is detailed in (9) below. The adjoint variables also must satisfy the transversality conditions for i ∈ {1, 2, 3}. Finally, the optimality conditions come from solving This results in Since we are considering bounded controls (because of (3)), the u * i are calculated using where u 1 max = u 3 max = 1 and u 2 max = 0.9. Uniqueness of the optimal controls (at least for small enough T ) also follow from standard results, such as Theorem 2.3 in [24] . Numerical solutions of systems (4) and (9) can be found by a forward-backward sweep method [26] . The algorithm starts with an initial guess of the controls u 1 , u 2 and u 3 and then solves (4) forward in time. After this first part, it uses the results and the initial guesses to solve (9) backward in time and new controls are defined following (13) . This process continues until it converges. Quarantines are not just a matter of public health, for they also present economic questions, for example. This means that we must pay close attention to the control costs B i . These numbers reflect how the population is capable of dealing with the quarantine of the respective age group. Smaller values of B i mean that the population can withstand a stricter quarantine without many economical side effects. This is not possible, on the other hand, for bigger values of B i . Since the bigger economic toll of the quarantine lies on the adults (because they form almost all of the economically active population), we assume that B 2 is the greatest of the three values. As the isolation of the young implies closing schools, this educational impact makes B 1 the second highest cost, albeit much smaller than B 2 . How the total cost B is distributed among the age groups determines the shape of the optimal controls. To study this relation, we let B = 5000, with B 2 ∈ [3600, 4600] and B 3 ∈ [10, 110]. For our simulations, we consider B 2 ∈ {3600 + 50k : 0 ≤ k ≤ 20} and B 3 ∈ {10 + 5k : 0 ≤ k ≤ 20}. This means that 441 cost distributions will be analysed for each quarantine length. Our goal now is to compare these distributions by investigating the number of deaths caused by the pandemic at the end of the quarantine for each one of them. As in [15] , the deaths will be calculated as a fraction of the recovered, since there is no disease induced mortality in our model. From Table 3 , we can derive the case fatality rates µ 1 , µ 2 and µ 3 for age groups 1, 2 and 3, respectively. The results are Let D(b 2 , b 3 , t) and R i (b 2 , b 3 , t), i ∈ {1, 2, 3}, denote the cumulative number of deaths due to the disease at time t and the number of recovered individuals in the optimal control problem for age group i, respectively, for the cost distribution with B 2 = b 2 and B 3 = b 3 . By our discussion above, we can write According to [22] , ratio estimates of reported to unreported cases vary from Figure 3 plots the graphs of the optimal controls u 1 (t), u 2 (t) and u 3 (t) for this cost distribution. An interesting feature of these plots is that they provide an "optimal calendar" of when the quarantine should start to be relaxed. This calendar is shown in Table 6 . We now investigate how the length of the quarantine influences its effectiveness by plotting the total number of infections for the cost distribution B 1 = 1390, B 2 = 3600 and B 3 = 10 in Figure 4 . The three curves reach their minimum around the same time that the controls reach zero. For quarantines of 45 and 60 days, the number of cases at the end is, indeed, much smaller than it is in the beginning. However, for the shorter quarantine, there are more cases in the end than the initial total. This fact shows that quarantines cannot be too short, or else the overall situation in the end could be worse than in the beginning. Moreover, the number of cases starts to go up again towards the end in all three cases. To finish this Section, we analyse how the initial conditions affect the optimal control. We can interpret this as a way to see what happens if it takes too long for these measures to be implemented. We do this by considering initial conditions of exposed, infected and recovered twice and four times as much as their original values. As of May 13, 2020, the number of active cases in Brazil doubles every 10 days [1] , so this means waiting 10 or 20 days, respectively, to start the quarantine. In the plots of Figures 5 and 6 , we consider distribution B 1 = 1390, B 2 = 3600 and B 3 = 10 with initial conditions twice and four times their original values, respectively. They provide different relaxation calendars, which are described in Tables 7 and 8. A comparison of Tables 6, 7 and 8 shows that the quarantines have to be much stricter if there is a delay in their implementation. Moreover, as we did in Figure 2 , we look at the reduction in deaths at the end of the quarantine in Figure 6 : Plots of the optimal controls for different quarantine lengths. Initial conditions of exposed, infected and recovered are multiplied by four. Hence, the longer it takes for the quarantine to start, the less effective it is. In this paper we considered an age-structured SEIRQ model, where the quarantine entrance parameters are thought as controls of the system, and we looked for the optimal controls via Pontryagin's maximum principle. After writing down the optimality system, we calculated the optimal controls numerically and analysed how some of the parameters influence the results. These parameters represent the difficult choices authorities must make, such as deciding how many essential workers are allowed to remain circulating, estimating the economic impact of the quarantine and even when to start it. As such choices are made, the optimal controls give guidelines of how to proceed. In Section 4, we considered a constant total control cost and distributed it among the age groups in 441 ways. The distribution with the best results with regard to deaths during the quarantine gave us a calendar of when to relax the isolation measures in the three age groups (for quarantine lengths of 30, 45 and 60 days, respectively): • For the youngs, the date of relaxation was the 16th, the 18th or the 14th day. • For the adults, the relaxation started at the 12th, the 13th or the 11th day. • For the elderly, it started at the 25th, the 39th or the 52nd day. The optimal controls that induce this calendar produced a reduction in the number of deaths of 286, 439 and 462 times, respectively, in comparison to the same period of time without quarantine. However, in the three cases the number of infected cases reached a minimum just before the end of the simulation, so by the time the quarantine ended, the cases were going up again, even becoming bigger than the original values for the shorter length we considered. This shows that the quarantines are not effective if they are not long enough. We also showed that waiting too long to start the quarantine makes the period before the relaxation become longer. This also produced a loss in efficacy, since the reduction of deaths due to the quarantine (in comparison to the "doing nothing" scenario) decreased as the number of initial cases increased. In our model we used data from Brazil as initial conditions and in the parameter fitting. Brazil is a very large country, with many cities at different stages of the pandemic. This means that studies such as this one should be made locally to best suit the characteristics of each city. As the plots of Figures 5 and 6 suggest, the sooner the quarantine is implemented, the shortest the time the controls need to stay at their maximum is. Finally, we acknowledge that not all individuals have the same possibility of fulfilling the quarantine. In countries with severe socioeconomic inequalities such as Brazil, the wealthier citizens have much more resources to go through the isolation period than the poorer [27] . We did not include these factors in our model. However, many countries around the world have provided financial support for those in need [28] in an attempt to mitigate this problem during the pandemic. Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in wuhan, china: a modelling study Optimal control of deterministic epidemics Optimal control of non-autonomous seirs models with vaccination and treatment Optimal control of an hiv immunology model Optimizing chemotherapy in an hiv model Optimal control of the chemotherapy of hiv Optimal control for a tuberculosis model with reinfection and post-exposure interventions Optimal control of treatments in a two-strain tuberculosis model Optimal control for pandemic influenza: the role of limited antiviral treatment and isolation Optimal quarantine strategies for covid-19 control models Modeling the control of covid-19: Impact of policy interventions and meteorological factors Assessing the efficiency of different control strategies for the covid-19 epidemic Disease extinction and disease persistence in age structured epidemic models Mathematical analysis of an age-structured sir epidemic model with vertical transmission Dynamics of a discrete age-structured sis models Global dynamics of a discrete age-structured sir epidemic model with applications to measles vaccination strategies Centro de coordinación de alertas y emergencias sanitarias goberno españa -enfermedad por el coronavirus documentos/Actualizacion_104_COVID-19.pdf Using a delay adjusted case fatality ratio to estimate under reporting The incubation period of coronavirus disease 2019 (covid-19) from publicly reported confirmed cases: estimation and application Optimal control methods applied to disease models Mathematical theory of optimal processes, Routledge Optimal control applied to biological models Lockdowns are fine for the rich, but millions are too poor to shelter from coronavirus Coronavirus bailouts: Which country has the most generous deal ☒ The authors declare that they have no known competing financial interests or personal relationships Conceptualization, Methodology, Software, Formal Analysis, Investigation, Data Curation, Writing (Original draft), Writing Software, Validation, Resources, Writing (Review and editing) The authors would like to thank César Castilho (UFPE) for all the valuable discussions and suggestions during the preparation of this manuscript.