Abstract:The aim of this paper is to provide a rigorous mathematical analysis of an optimal control problem of a SIR epidemic on an infinite horizon. A state constraint related to intensive care units (ICU) capacity is imposed and the objective functional linearly depends on the state and the control. After preliminary asymptotic and viability analyses, a $\Gamma$-convergence argument is developed to reduce the problem to a finite horizon allowing to use a state constrained version of Pontryagin's theorem to characterize the structure of the optimal controls. Illustrating examples and numeric simulations are given according to the available data on Covid-19 epidemic in Italy.
What problem does this paper attempt to address?
This paper aims to solve an optimal control problem of the SIR epidemic in an infinite - time horizon, in which a state constraint related to the intensive care unit (ICU) capacity is introduced. Specifically, the objective function of the study depends linearly on the state variables (i.e., the number of infected people) and the control variables (such as non - pharmaceutical interventions). Through preliminary asymptotic analysis and viability analysis, the authors use the Γ - convergence method to simplify the problem to a finite - time - horizon problem, so that the Pontryagin maximum principle with state constraints can be used to describe the structure of the optimal control.
### Problems the paper attempts to solve
1. **Optimal control of the SIR epidemic in an infinite - time horizon**:
- Research how to optimally control the SIR epidemic model in an infinite - time horizon to minimize the total cost related to the number of infected people and control measures.
- The specific goal is to find an optimal control strategy \( u(t) \) such that the linear integral cost function \( J_\infty(u, s, i) \) over the entire time horizon is minimized:
\[
J_\infty(u, s, i) = \int_0^\infty (\lambda_1 u(t) + \lambda_2 i(t)) \, dt
\]
where \( s(t) \) and \( i(t) \) represent the density functions of the susceptible and infected populations respectively, and \( u(t) \) represents control measures (such as isolation, social distancing, etc.).
2. **ICU capacity constraint**:
- Introduce a state constraint related to ICU capacity \( i(t) \leq i_M \), where \( i_M \) is the safe upper limit of the ICU, representing the maximum number of infected people that the medical system can handle.
- This constraint ensures that the medical system will not collapse due to too many infected people.
3. **Description of the optimal control structure**:
- Use the Γ - convergence method to simplify the infinite - time - horizon problem to a finite - time - horizon problem, and then use the Pontryagin maximum principle with state constraints to describe the structure of the optimal control.
- Verify and illustrate the theoretical results through numerical simulations and specific examples (such as COVID - 19 epidemic data in Italy).
### Main contributions
- **Mathematical analysis**: Provide strict mathematical analysis to prove the existence and uniqueness of the optimal control problem.
- **Γ - convergence method**: Use the Γ - convergence method to simplify the infinite - time - horizon problem to a finite - time - horizon problem, so that more mature optimal control theories can be applied.
- **Pontryagin maximum principle**: Use the Pontryagin maximum principle with state constraints to describe the structure of the optimal control, providing a theoretical basis for practical applications.
- **Numerical simulation**: Verify the theoretical results through specific numerical simulations and examples (such as COVID - 19 epidemic data in Italy), showing the practical effects of the optimal control strategy.
In conclusion, through strict mathematical analysis and advanced optimal control theories, this paper solves the problem of how to optimally control the SIR epidemic model in an infinite - time horizon when considering ICU capacity constraints.