Abstract:<p>An optimal control problem of a general reaction-diffusion tumor-immune system with chemotherapy is investigated to minimize the tumor burden and side effects as well as treatment costs. Firstly, the existence, uniqueness and some estimates of strong solution to the state system in spatial dimensions <span class="math"><math>n=1,2,3</math></span>are obtained by making use of the semigroup theory and truncation method. Subsequently, we prove the existence of optimal pair by utilizing the minimizing sequence technique. Furthermore, we show the differentiability of the control-to-state mapping and establish the first-order necessary optimality condition. Finally, several numerical simulations are presented to illustrate the practical application of the theoretical results obtained in this work and to validate some clinical observations.</p>
automation & control systems,engineering, electrical & electronic, multidisciplinary,mathematics, interdisciplinary applications
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the optimal control problem of the tumor - immune system reaction - diffusion model under chemotherapy conditions. Specifically, the author aims to minimize the tumor burden, side effects and treatment costs by optimizing the administration strategy of chemotherapy drugs. This research not only considers the interaction between tumor cells and immune cells, but also includes the effects of chemotherapy drugs on normal cells, tumor cells and immune cells. To achieve these goals, a general reaction - diffusion tumor - immune system model including the influence of chemotherapy was constructed in the paper, and the optimal control problem of this model was explored by mathematical methods.
### Main contributions of the paper
1. **First study on the optimal control problem of the diffusion - type tumor - immune model**: As far as the author knows, this is the first study on the optimal control problem of the diffusion - type tumor - immune model.
2. **Generality of the model**: The model used in the research is a general reaction - diffusion tumor - immune model, not a specific - type model, which makes the research results more universal.
3. **Combination of theoretical and numerical analysis**: The paper not only theoretically proves the existence and necessary conditions of the optimal control problem, but also verifies the validity of the theoretical results through numerical simulation.
### Mathematical model
The mathematical model in the paper is described as follows:
\[
\begin{cases}
\frac{\partial N}{\partial t} = D_N \Delta N + N g_1(N) - c_1 N T - N f_1(U), & \text{in } Q_{tf}, \\
\frac{\partial T}{\partial t} = D_T \Delta T + T g_2(T) - c_2 T N - c_3 T I - T f_2(U), & \text{in } Q_{tf}, \\
\frac{\partial I}{\partial t} = D_I \Delta I + s + \frac{\rho I T}{\alpha + T} - c_4 I T - d_1 I - I f_3(U), & \text{in } Q_{tf}, \\
\frac{\partial U}{\partial t} = D_U \Delta U - d_2 U + u, & \text{in } Q_{tf}, \\
\frac{\partial N}{\partial \nu} = \frac{\partial T}{\partial \nu} = \frac{\partial I}{\partial \nu} = \frac{\partial U}{\partial \nu} = 0, & \text{on } \Sigma_{tf}, \\
N(x,0) = N_0, \quad T(x,0) = T_0, \quad I(x,0) = I_0, \quad U(x,0) = U_0, & \text{in } \Omega,
\end{cases}
\]
where:
- \( N(x,t) \) represents the density of normal cells,
- \( T(x,t) \) represents the density of tumor cells,
- \( I(x,t) \) represents the density of immune cells,
- \( U(x,t) \) represents the concentration of chemotherapy drugs,
- \( D_N, D_T, D_I, D_U \) represent the diffusion coefficients of normal cells, tumor cells, immune cells and chemotherapy drugs respectively,
- \( g_1(N), g_2(T) \) represent the growth rates of normal cells and tumor cells respectively,
- \( f_1(U), f_2(U), f_3(U) \) represent the killing rates of chemotherapy drugs on normal cells, tumor cells and immune cells respectively,
- \( u(x,t) \) represents the amount of externally input chemotherapy drugs.
### Objective function
To balance the efficacy, side effects and cost of chemotherapy, the following objective function is defined in the paper:
\[
J(T(u), U(