pith. sign in

arxiv: 2604.21674 · v1 · submitted 2026-04-23 · 🧮 math.OC

Optimal control of therapies related to an oxytaxis glioblastoma model

Juan J. Forero-Herna\'ndez , Francisco Guill\'en-Gonz\'alez , \'Elder J. Villamizar-Roa This is my paper

Pith reviewed 2026-05-09 21:04 UTC · model grok-4.3

classification 🧮 math.OC
keywords optimal controlglioblastomaKeller-Segel systemoxytaxischemotherapyantiangiogenic therapyparabolic PDEnumerical optimization
0
0 comments X

The pith

Existence of global optimal solutions is proven for controlling glioblastoma growth using chemotherapy and antiangiogenic therapies in an oxytaxis model.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors formulate an optimal control problem for a Keller-Segel type parabolic system modeling glioblastoma cell movement toward oxygen gradients in a two-dimensional bounded domain. They prove that a global optimal solution exists for a cost functional that minimizes both tumor cell density and oxygen concentration, with the two therapies serving as controls. Necessary first-order optimality conditions are derived from the system, and a numerical scheme using the adjoint method and Adam gradient optimization is developed to approximate the therapies. A sympathetic reader would care because the work supplies a rigorous mathematical route to optimize treatment intensities in a biologically detailed PDE model of cancer dynamics driven by oxygen.

Core claim

The paper establishes that for the optimal control problem associated with the oxytaxis glioblastoma model, a global optimal solution exists and first-order necessary optimality conditions can be obtained. The state system is a parabolic Keller-Segel-type model with random diffusion, oxytaxis, and reaction terms between tumor cells and oxygen. The cost functional penalizes both tumor growth and oxygen levels through the two therapeutic controls. A gradient-based numerical method is then used to approximate the optimal therapies in simulations.

What carries the argument

The optimal control problem for the Keller-Segel-type parabolic system with oxytaxis, where the two therapy intensities enter the reaction terms and optimality is characterized through the adjoint system derived from the state equations.

If this is right

  • Global optimal therapy functions exist that achieve the infimum of the tumor-plus-oxygen cost.
  • The first-order optimality conditions characterize the optimal controls via the adjoint state.
  • The adjoint-based gradient together with the Adam optimizer produces computable approximations to those optimal therapies.
  • Numerical experiments on the two-dimensional domain confirm that the computed therapies reduce the cost functional.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same existence and optimality framework could be applied to three-dimensional domains or models that include additional cell types such as immune cells.
  • If the model parameters are calibrated to patient imaging data, the numerical scheme could generate individualized therapy schedules.
  • Similar control analysis might be carried out for other taxis-driven biological processes where two external interventions are available.
  • Direct comparison of the predicted optimal dosing with clinical trial data would test whether the model's oxygen-taxis assumption holds in practice.

Load-bearing premise

The underlying parabolic system admits unique weak-strong solutions for the chosen reaction terms and boundary conditions.

What would settle it

A concrete parameter regime or boundary condition in which the PDE system loses uniqueness of weak-strong solutions, or a numerical simulation in which the adjoint-derived gradient fails to produce a control that reduces the cost functional.

Figures

Figures reproduced from arXiv: 2604.21674 by \'Elder J. Villamizar-Roa, Francisco Guill\'en-Gonz\'alez, Juan J. Forero-Herna\'ndez.

Figure 1
Figure 1. Figure 1: Initial settings for the simulations. We define the initial condition for tumor density u as a Gaussian distribution localized at the center of the domain (see Figure 1b). Denoting by Lx and Ly the horizontal and vertical lengths of the domain, respectively, and by (xc, yc) the center point of the domain, with a width based on the size of the domain, we consider u0(x, y) = 0.6 exp  − (x − xc) 2 + (y − yc)… view at source ↗
Figure 2
Figure 2. Figure 2: First times of evolution. (a) t = 0.25 (b) t = 0.5 (c) t = 1 (d) t = 0.25 (e) t = 0.5 (f) t = 1 [PITH_FULL_IMAGE:figures/full_fig_p019_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: illustrates the evolution of the system over a longer time horizon up to T = 1, where random diffusion becomes significant in conjunction with oxytaxis effects. (a) Initial (b) t = 0.008 (c) t = 0.016 (d) t = 0.032 (e) Initial (f) t = 0.08 (g) t = 0.016 (h) t = 0.032 [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: shows the formation of a cell aggregation point; that is, cells tend to concentrate in regions where the gradient of σ is favorable. This aggregation phenomenon is directly related to the oxytaxis term and its sensitivity parameter χ. In particular, the tumor reaches its maximum value of aggregation at t = 1.664 (see Figure 4c). (a) t = 1.25 (b) t = 1.5 (c) t = 1.664 (d) t = 2 (e) t = 1.25 (f) t = 1.5 (g) … view at source ↗
Figure 5
Figure 5. Figure 5: Evolution over time [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Optimal control compared with initial control. [PITH_FULL_IMAGE:figures/full_fig_p021_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Evolution over time. 0 50 100 150 200 250 300 4 6 8 10 12 Iteration Functional J Functional J Decrease Increase No change (a) Functional evolution (case with lj = 1) 0 50 100 150 200 250 300 3.5 4 Iteration Functional J Functional J Decrease Increase No change (b) Functional evolution (case with lj = 0) -0.0005 -0.0004 -0.0003 -0.0002 -0.0001 0 0.005 0.01 0.015 0.02 0.025 0.03 4.37 4.37 4.37 4.37 pert Func… view at source ↗
Figure 8
Figure 8. Figure 8: Stability of functional and minimum local. [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Optimal control compared with initial control. [PITH_FULL_IMAGE:figures/full_fig_p023_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: shows that the increase in the amount of control c, allows to obtain a greater regulation of both the volume (see Figure 10d) and the L 2 -norm of u (see Figure 10a); this fact demonstrates a greater efficiency of the optimal control strategy. Furthermore, it facilitates enhanced control of the aggregation point (see Figure 10b). 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 5 10 15 20 Time R u 2 Uncontrolled O… view at source ↗
read the original abstract

We propose and analyze an optimal control problem associated with a Keller-Segel type parabolic system with chemoattraction, modeling the glioblastoma growth in a bi-dimensional bounded domain, influenced by the presence of oxygen where the controls are two different (chemotherapy and antiangiogenic) therapies. The model considers the random diffusion of tumor cells and oxygen, the movement of cells towards the oxygen gradient (oxytaxis), and reaction terms describing the interaction between cells and oxygen. We establish a mathematical framework to analyze the existence and uniqueness of weak-strong solution of the model and subsequently we analyze an optimal control problem considering a cost functional that minimizes both the tumor growth and the oxygen concentration. We prove the existence of a global optimal solution and derive necessary first-order optimality conditions. Finally, we propose a methodology for approximating the optimal therapies. We use the gradient of the reduced cost functional through the adjoint scheme, and minimize the cost functional implementing the Adam gradient optimization method. Some numerical experiments are provided to demonstrate the effectiveness of the proposed scheme.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript proposes an optimal control problem for a Keller-Segel-type parabolic system modeling glioblastoma growth in a bounded 2D domain, incorporating oxytaxis toward oxygen gradients. Controls represent chemotherapy and antiangiogenic therapies. The authors construct a framework proving existence and uniqueness of weak-strong solutions to the state system, establish existence of a global optimal solution pair, derive first-order necessary optimality conditions via an adjoint approach, and present a numerical approximation scheme that computes the gradient of the reduced cost functional (penalizing tumor density and oxygen levels) using the adjoint and minimizes it via the Adam optimizer, illustrated by numerical experiments.

Significance. If the well-posedness and optimality results hold, the work supplies a rigorous analytical and computational framework for therapy optimization in a biologically motivated PDE model with chemoattraction. A clear strength is the explicit construction of the adjoint-based gradient for the reduced functional combined with a modern optimizer (Adam) and supporting numerical experiments that demonstrate practical effectiveness. This combination of existence theory, necessary conditions, and reproducible numerics adds value for applied optimal control in cancer modeling, provided the core estimates are verified.

major comments (2)
  1. [Well-posedness section (likely §2 or §3)] The existence and uniqueness of weak-strong solutions to the state system (claimed in the abstract and developed in the well-posedness analysis) is load-bearing for both the optimal-control existence proof and the differentiability needed for the adjoint system. The manuscript must supply the explicit a-priori estimates (energy bounds or maximum-principle arguments) that control the quadratic oxytaxis term together with the chosen proliferation and consumption reaction terms to preclude finite-time blow-up on the bounded 2-D domain; any gap here renders the control-to-state map undefined.
  2. [Optimal control and adjoint section (likely §4)] In the derivation of the first-order optimality conditions (via the adjoint system), the paper assumes sufficient regularity of the control-to-state operator. The precise function spaces for the controls (e.g., L^∞ or L^2) and the states must be stated explicitly so that the linearized system and passage to the limit are justified; without this, the necessary conditions cannot be rigorously obtained.
minor comments (3)
  1. [Abstract] The abstract asserts proofs of existence, uniqueness, optimality conditions, and a numerical scheme; a short sentence indicating the main technical tools (e.g., Galerkin approximation plus compactness or fixed-point arguments) would orient readers without lengthening the abstract.
  2. [Numerical experiments] Numerical experiments section: the description of the Adam implementation should specify the learning rate schedule, iteration count, tolerance, and how the adjoint gradient is discretized to permit full reproducibility of the reported optimal therapies.
  3. [Figures] Figure captions for the simulation results should list the exact parameter values (diffusion coefficients, reaction rates, initial data) used, consistent with the model equations, to strengthen the reproducibility claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and insightful comments on our manuscript. We address each major comment below and plan to incorporate clarifications and additional details in the revised version to strengthen the rigor of our results.

read point-by-point responses

  1. Referee: [Well-posedness section (likely §2 or §3)] The existence and uniqueness of weak-strong solutions to the state system (claimed in the abstract and developed in the well-posedness analysis) is load-bearing for both the optimal-control existence proof and the differentiability needed for the adjoint system. The manuscript must supply the explicit a-priori estimates (energy bounds or maximum-principle arguments) that control the quadratic oxytaxis term together with the chosen proliferation and consumption reaction terms to preclude finite-time blow-up on the bounded 2-D domain; any gap here renders the control-to-state map undefined.

    Authors: We appreciate the referee's emphasis on the foundational role of the well-posedness results. In the manuscript, the existence and uniqueness of weak-strong solutions are established in Section 3 through a combination of Galerkin approximation, a priori estimates, and compactness arguments. The a priori bounds are obtained by testing the tumor cell equation with the density itself, utilizing the 2D Sobolev embedding to control the quadratic oxytaxis term (via |∇c| term bounded by L^2 norms), and leveraging the dissipative structure from the oxygen consumption and logistic proliferation terms to prevent blow-up. These estimates are global in time on the bounded domain. To make this more transparent, we will expand the well-posedness section with a dedicated lemma detailing these energy estimates and maximum principle arguments in the revision. revision: yes

  2. Referee: [Optimal control and adjoint section (likely §4)] In the derivation of the first-order optimality conditions (via the adjoint system), the paper assumes sufficient regularity of the control-to-state operator. The precise function spaces for the controls (e.g., L^∞ or L^2) and the states must be stated explicitly so that the linearized system and passage to the limit are justified; without this, the necessary conditions cannot be rigorously obtained.

    Authors: We agree that explicit specification of the function spaces is crucial for rigor. The controls are chosen in the admissible set U_ad subset L^∞(0,T; L^∞(Ω)) with pointwise bounds, ensuring boundedness and allowing for the application of the adjoint method. The state variables belong to appropriate spaces such as L^2(0,T; H^1(Ω)) ∩ L^∞(0,T; L^2(Ω)) for the tumor density and similar for oxygen. The differentiability of the control-to-state map is justified by the Lipschitz continuity from the well-posedness estimates. We will revise the optimal control section to explicitly state these spaces and provide a brief outline of the linearization and limit passage arguments. revision: yes

Circularity Check

0 steps flagged

No circularity: standard PDE well-posedness and adjoint derivation are independent of the target optimality result

full rationale

The paper first constructs a framework proving existence and uniqueness of weak-strong solutions to the Keller-Segel oxytaxis system (using a-priori estimates, compactness, and handling of the quadratic term on a bounded domain), then invokes this to establish existence of a global optimal control pair and to derive first-order necessary conditions via the adjoint system. These steps rely on external parabolic theory and standard optimal-control arguments rather than any self-definition, fitted-parameter renaming, or self-citation chain that reduces the claimed result to its own inputs by construction. No prediction or optimality condition is shown to be equivalent to a fitted quantity or prior ansatz from the same authors.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Insufficient detail in the abstract to enumerate specific free parameters, axioms, or invented entities; the model is described only at the level of diffusion, oxytaxis, and reaction terms without explicit coefficient values or additional postulates.

pith-pipeline@v0.9.0 · 5490 in / 1084 out tokens · 26574 ms · 2026-05-09T21:04:56.984950+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    & Thompson A

    Anderson A.; Chaplin M.; Newman E.; Steele R. & Thompson A. (2000). Mathematical modelling of tumour invasion and metastasis.Computational and Mathematical Methods in Medicine.2(2), 129–154

    work page 2000

  2. [2]

    & Mark C

    Badal J.; Xueying W.; Sayanti B.; Haiyan T.; Anastasios M. & Mark C. (2009). On immunothera- pies and cancer vaccination protocols: a mathematical modelling approach,J. Theor. Biol.259(4) 820–827

    work page 2009

  3. [3]

    & Winkler M

    Bellomo N.; Bellouquid A.; Tao Y. & Winkler M. (2015). Toward a mathematical theory of Keller- Segel models of pattern formation in biological tissues.Math. Models Methods Appl. Sci.25(09), 1663–1763

    work page 2015

  4. [4]

    Braz e Silva, P.; Guill´ en-Gonz´ alez, F.; Perusato, C.& Rodr´ ıguez-Bellido, M. A. (2023). Bilinear opti- mal control of the Keller–Segel logistic model in 2D-domains.Applied Mathematics & Optimization, 87, 55. A parabolic system with chemoattraction modeling the growth of glioma cells30

    work page 2023

  5. [5]

    Chaplain M. (2000). Mathematical modelling of angiogenesis.Journal of Neuro-Oncology.50(1-2), 37–51

    work page 2000

  6. [6]

    Lolas G.(2006)

    Chaplain M. Lolas G.(2006). Mathematical modelling of cancer invasion of tissue: dynamic hetero- geneity.Networks &; Heterogeneous Media1(13), 399–439

    work page 2006

  7. [7]

    & Rocca, E

    Colli, P.; Gomez, H.; Lorenzo, G.; Marinoschi, G.; Reali, A. & Rocca, E. (2021) Optimal control of cytotoxic and antiangiogenic therapies on prostate cancer growth.Mathematical Models and Methods in Applied Sciences.31(7), 1419–1468

    work page 2021

  8. [8]

    & Liu, B

    Dai, F. & Liu, B. (2019). Optimal control and pattern formation for a haptotaxis model of solid tumor invasion.Journal of the Franklin Institute,356(16), 9364–9406

    work page 2019

  9. [9]

    & de Magalh˜ aes P

    de Araujo A. & de Magalh˜ aes P. (2015). Existence of solutions and optimal control for a model of tissue invasion by solid tumours.J. Math. Anal. Appl.421(1), 842–877

    work page 2015

  10. [10]

    (2022) Optimal control of a Fisher-KPP type model arising in chemotherapy for brain tumors.Pure and Applied Functional Analysis,7(5) 1637-1655

    Fernandez, L. (2022) Optimal control of a Fisher-KPP type model arising in chemotherapy for brain tumors.Pure and Applied Functional Analysis,7(5) 1637-1655

    work page 2022

  11. [11]

    & Su´ arez, A

    Fern´ andez-Romero, A.; Guill´ en-Gonz´ alez, F. & Su´ arez, A. (2022). A glioblastoma PDE–ODE model including chemotaxis and vasculature.ESAIM: Mathematical Modelling and Numerical Analysis, 56(2), 407–431

    work page 2022

  12. [12]

    & Novotn´ y, A

    Feireisl, E. & Novotn´ y, A. (2009).Singular limits in thermodynamics of viscous fluids. Basel: Birkh¨ auser

    work page 2009

  13. [13]

    G´ omez, H. (2017). Quantitative analysis of the proliferative-to-invasive transition of hypoxic glioma cells.Integrative Biology,9, 257–262

    work page 2017

  14. [14]

    & Painter K

    Hillen T. & Painter K. (2009). A user’s guide to PDE models for chemotaxis.J. Math. Biol.58(1-2), 183–217

    work page 2009

  15. [15]

    & Ba, L.J

    Kingma, D.P. & Ba, L.J. 2015 Adam: A Method for Stochastic Optimization.International Con- ference on Learning Representations (ICLR)

    work page 2015

  16. [16]

    F., & Segel, L

    Keller, E. F., & Segel, L. A. (1971). Model for chemotaxis.Journal of theoretical biology,30(2), 225–234

    work page 1971

  17. [17]

    (1969).Quelques m´ ethodes de r´ esolution des probl` emes aux limites non lin´ eaires

    Lions, L., Jean J. (1969).Quelques m´ ethodes de r´ esolution des probl` emes aux limites non lin´ eaires. Paris: Dunod

    work page 1969

  18. [18]

    L.; Rueda-G´ omez, D

    L´ opez-Agredo, J. L.; Rueda-G´ omez, D. A.; Villamizar-Roa,´E. J. (2023). Theoretical and numerical analysis of a parabolic system with chemoattraction modeling the growth of glioma cells.Applied Numerical Mathematics,186, 143–163

    work page 2023

  19. [19]

    L.; Rueda-G´ omez, D

    L´ opez-Agredo, J. L.; Rueda-G´ omez, D. A.; Villamizar-Roa, ´E. J. (2023). Numerical analysis of a mathematical model describing the evolution of hypoxic glioma cells, Comput. Math. Appl. 131 (2023) 138–157

    work page 2023

  20. [20]

    & Kalpathy J

    Menze B.; Bauer S. & Kalpathy J. (2015). The multimodal brain tumor image segmentation bench- mark (BRATS).IEEE Trans. Med. Imaging.34(10), 1993–2024

    work page 2015

  21. [21]

    Ni˜ no-Celis, V; Rueda-G´ omez, D. A. & Villamizar-Roa,´E. J. (2021). Convergence and positivity of finite element methods for a haptotaxis model of tumoral invasion.Comput. Math. Appl.8920–33

    work page 2021

  22. [22]

    Swanson, K. R. & Alvord, E. C. (2000). A quantitative model for differential motility of gliomas in grey and white matter.Cell Proliferation,33(5), 317–329. A parabolic system with chemoattraction modeling the growth of glioma cells31

    work page 2000

  23. [23]

    Simon, J. (1987). Compact sets in the spaceL p(0, T;B).Annali di Matematica Pura ed Applicata, 146(4), 65–96

    work page 1987

  24. [24]

    & Lankeit J

    Winkler M. & Lankeit J. (2020). Facing low regularity in chemotaxis systems.Kahresber. Dtsh. Math.-Ver.122(1), 35–64

    work page 2020