pith. sign in

arxiv: 1907.10452 · v1 · pith:BPIZCX4Pnew · submitted 2019-07-24 · 🧮 math.OC · q-bio.TO

A distributed control problem for a fractional tumor growth model

Pierluigi Colli , Gianni Gilardi , J\"urgen Sprekels This is my paper

Pith reviewed 2026-05-24 16:49 UTC · model grok-4.3

classification 🧮 math.OC q-bio.TO
keywords distributed optimal controlfractional operatorstumor growth modelCahn-Hilliard systemFréchet differentiabilityadjoint systemoptimality conditions
0
0 comments X

The pith

The control-to-state operator for a fractional tumor growth system is shown to be Fréchet differentiable.

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

The paper considers a system of three evolutionary equations with fractional powers of selfadjoint monotone operators having compact resolvents, generalizing a Cahn-Hilliard phase-field model of tumor growth. It proves Fréchet differentiability of the map taking distributed controls to states, establishes existence for the corresponding adjoint system, and obtains first-order necessary conditions of optimality for a tracking-type cost. A sympathetic reader would care because the result supplies a route to optimizing interventions such as drug administration while the fractional setting allows more flexible diffusion regimes than the Laplacian alone.

Core claim

The central claim is that the control-to-state operator associated with the system of three evolutionary equations involving fractional powers is Fréchet differentiable, that the associated adjoint system admits solutions, and that first-order necessary conditions of optimality hold for the tracking-type cost functional.

What carries the argument

The control-to-state operator, whose Fréchet differentiability is established by means of the fractional powers of the three given operators.

If this is right

  • First-order necessary optimality conditions become available for the distributed control problem.
  • The adjoint system is solvable under the stated operator assumptions.
  • The results cover cases in which the three diffusional operators may differ from one another and from the Laplacian.

Where Pith is reading between the lines

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

  • The same differentiability argument might apply to other phase-field systems once the operator hypotheses are met.
  • Numerical gradient-based optimization schemes could be constructed directly from the derived optimality conditions.
  • The framework could be tested on heterogeneous domains by assigning different fractional orders to the three operators.

Load-bearing premise

The three unbounded linear operators are selfadjoint, monotone, and possess compact resolvents.

What would settle it

An explicit choice of control for which the state map fails to be differentiable or for which the adjoint system has no solution, while the three operators still satisfy selfadjointness, monotonicity, and compact resolvent.

read the original abstract

In this paper, the authors study the distributed optimal control of a system of three evolutionary equations involving fractional powers of three selfadjoint, monotone, unbounded linear operators having compact resolvents. The system is a generalization of a Cahn-Hilliard type phase field system modeling tumor growth that goes back to Hawkins-Daarud et al. (Int. J. Numer. Math. Biomed. Eng. 28 (2012), 3-24). The aim of the control process, which could be realized by either administering a drug or monitoring the nutrition, is to keep the tumor cell fraction under control while avoiding possible harm for the patient. In contrast to previous studies, in which the occurring unbounded operators governing the diffusional regimes were all given by the Laplacian with zero Neumann boundary conditions, the operators may in our case be different; more generally, we consider systems with fractional powers of the type that were studied in the recent work Adv. Math. Sci. Appl. 28 (2019), 343-375 (see arXiv:1906.10874), by the present authors. In the analysis, the Fr\'echet differentiability of the associated control-to-state operator is shown, by also establishing the existence of solutions to the associated adjoint system, and deriving the first-order necessary conditions of optimality for a cost functional of tracking type.

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 / 2 minor

Summary. The manuscript analyzes a distributed optimal control problem for a system of three coupled evolutionary PDEs driven by fractional powers A^α, B^β, C^γ of distinct self-adjoint monotone unbounded linear operators with compact resolvents. This generalizes a Cahn-Hilliard-type tumor growth model. The central results are the Fréchet differentiability of the control-to-state map, existence of solutions to the associated adjoint system, and derivation of first-order necessary optimality conditions for a tracking-type cost functional.

Significance. If the claims hold, the work extends optimal control theory for phase-field tumor models to heterogeneous diffusion regimes via distinct fractional operators, building directly on the authors' prior framework for fractional powers. This could enable more flexible modeling of nutrient/drug distribution in non-uniform tissues. The explicit construction of the adjoint and optimality conditions for the generalized setting is a technical contribution.

major comments (2)
  1. [§3.3, Theorem 3.8] §3.3, Theorem 3.8 (linearized system): The proof that the linearized operator around a solution maps into the required dual spaces for distinct A, B, C relies on the spectral theorem applied separately to each operator, but does not verify that the cross-coupling terms (arising from the nonlinearities) remain bounded when the eigenbases differ; this is load-bearing for the subsequent Fréchet differentiability claim in §4.
  2. [§4.2, Eq. (4.12)] §4.2, Eq. (4.12) (adjoint system): Existence of the adjoint is obtained by transposition from the linearized state equation, yet the argument assumes the same a-priori estimates as the forward problem without additional structural hypotheses (e.g., commutativity or uniform sectoriality) that would guarantee the dual pairing is well-defined for arbitrary α, β, γ; this directly affects the derivation of the optimality condition in Theorem 4.5.
minor comments (2)
  1. [§2] Notation for the fractional powers is introduced in §2 but the precise definition of the domains D(A^α) ∩ D(B^β) etc. is used without explicit reference to the spectral decomposition in the estimates of §3.1.
  2. [§1] The cost functional is stated as tracking type in the abstract and §1, but the precise form (including any regularization terms on the control) appears only in §4.1; a forward reference would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the two major comments point by point below, providing clarifications and indicating planned revisions to strengthen the presentation.

read point-by-point responses

  1. Referee: [§3.3, Theorem 3.8] §3.3, Theorem 3.8 (linearized system): The proof that the linearized operator around a solution maps into the required dual spaces for distinct A, B, C relies on the spectral theorem applied separately to each operator, but does not verify that the cross-coupling terms (arising from the nonlinearities) remain bounded when the eigenbases differ; this is load-bearing for the subsequent Fréchet differentiability claim in §4.

    Authors: We appreciate the referee highlighting this aspect of the proof. The cross-coupling terms arising from the nonlinearities are controlled via the Lipschitz continuity of the nonlinear maps and the regularity of the state solutions in the spaces H^1 and L^2, which yield uniform bounds on the products through standard embeddings; these estimates do not depend on the specific eigenbases of A, B, and C. The spectral theorem is invoked only to define the fractional powers individually, while the boundedness of the linearized operator into the dual spaces follows from the monotonicity and the a-priori estimates already established for the forward system. To address the concern explicitly, we will insert a short paragraph after the application of the spectral theorem in the proof of Theorem 3.8 that verifies the cross terms remain bounded independently of basis alignment. revision: partial

  2. Referee: [§4.2, Eq. (4.12)] §4.2, Eq. (4.12) (adjoint system): Existence of the adjoint is obtained by transposition from the linearized state equation, yet the argument assumes the same a-priori estimates as the forward problem without additional structural hypotheses (e.g., commutativity or uniform sectoriality) that would guarantee the dual pairing is well-defined for arbitrary α, β, γ; this directly affects the derivation of the optimality condition in Theorem 4.5.

    Authors: The transposition method for the adjoint system is justified by the well-posedness of the linearized forward problem (Theorem 3.8), whose estimates rely solely on the self-adjointness, monotonicity, and compact resolvent properties of each operator separately, together with the structure of the coupling. No commutativity between A, B, and C is used; the dual pairing is well-defined because the test functions are dense in the state space and the bilinear form associated with the linearized operator is continuous on the product space for arbitrary positive exponents α, β, γ. We will add a clarifying remark in §4.2 explaining why the same a-priori estimates carry over to the transposed equation without extra hypotheses, thereby supporting the optimality conditions in Theorem 4.5. revision: partial

Circularity Check

0 steps flagged

No circularity: self-citation supports base model only; control-to-state differentiability and adjoint are independent derivations

full rationale

The paper cites its prior work (arXiv:1906.10874) solely for the well-posedness framework of the fractional-power system under the stated operator assumptions (self-adjoint, monotone, compact resolvent). The central results—Fréchet differentiability of the control-to-state map, existence of the adjoint system, and first-order optimality conditions—are established via new analysis in this manuscript and do not reduce by the paper's own equations or definitions to quantities fixed in the cited work. No self-definitional loops, fitted inputs renamed as predictions, or load-bearing uniqueness claims appear. The derivation chain remains self-contained against the external mathematical benchmarks of the prior well-posedness result.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on functional-analytic assumptions about the operators and on the well-posedness of the underlying state system, both inherited from the authors' prior work rather than re-derived here.

axioms (1)
  • domain assumption The three unbounded linear operators are selfadjoint, monotone, and have compact resolvents.
    Invoked in the abstract to justify the fractional powers and the generalization beyond the Laplacian.

pith-pipeline@v0.9.0 · 5781 in / 1210 out tokens · 28649 ms · 2026-05-24T16:49:32.543212+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

48 extracted references · 48 canonical work pages · 2 internal anchors

  1. [1]

    Ainsworth, Z

    M. Ainsworth, Z. Mao, Analysis and approximation of a fractional Cahn–Hilliard equation, SIAM J. Numer. Anal. 55 (2017), 1689-1718

    work page 2017

  2. [2]

    Akagi, G

    G. Akagi, G. Schimperna, A. Segatti, Fractional Cahn–Hilliard, Alle n–Cahnn, and porous medium equations, J. Differential Equations 261 (2016), 2935-2985

    work page 2016

  3. [3]

    Op´ erateurs maximaux monotones et semi-groupes de contractions dans les espaces de Hilbert

    H. Brezis, “Op´ erateurs maximaux monotones et semi-groupes de contractions dans les espaces de Hilbert”, North-Holland Math. Stud. 5, North-Holland, Amsterdam, 1973

    work page 1973

  4. [4]

    Cavaterra, E

    C. Cavaterra, E. Rocca, H. Wu, Long-time dynamics and optimal control of a diffuse interface model for tumor growth, Appl. Math. Optim. (2019), https://doi.org/10.1007/s00245-019-09562-5

    work page doi:10.1007/s00245-019-09562-5 2019

  5. [5]

    Well-posedness, regularity and asymptotic analyses for a fractional phase field system

    P. Colli, G. Gilardi, Well-posedness, regularity and asymptotic analy ses for a fractional phase field system, Asymptot. Anal. , to appear (see also preprint arXiv:1806.04625 [math.AP] (2018), pp. 1-34)

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  6. [6]

    Colli, G

    P. Colli, G. Gilardi, D. Hilhorst, On a Cahn–Hilliard type phase field syst em related to tumor growth, Discrete Contin. Dyn. Syst. 35 (2015), 2423-2442

    work page 2015

  7. [7]

    Colli, G

    P. Colli, G. Gilardi, G. Marinoschi, E. Rocca, Sliding mode control for a phase field system related to tumor growth, Appl. Math. Optim. 79 (2019), 647-670

    work page 2019

  8. [8]

    Colli, G

    P. Colli, G. Gilardi, E. Rocca, J. Sprekels, Vanishing viscosities and e rror estimate for a Cahn–Hilliard type phase field system related to tumor growth, Nonlinear Anal. Real World Appl. 26 (2015), 93-108

    work page 2015

  9. [9]

    Colli, G

    P. Colli, G. Gilardi, E. Rocca, J. Sprekels, Optimal distributed cont rol of a diffuse interface model of tumor growth, Nonlinearity 30 (2017), 2518-2546

    work page 2017

  10. [10]

    Colli, G

    P. Colli, G. Gilardi, E. Rocca, J. Sprekels, Asymptotic analyses an d error estimates for a Cahn–Hilliard type phase field system modelling tumor growth, Discrete Contin. Dyn. Syst. Ser. S 10 (2017), 37-54

    work page 2017

  11. [11]

    Colli, G

    P. Colli, G. Gilardi, J. Sprekels, Optimal velocity control of a visco us Cahn–Hilliard system with convection and dynamic boundary conditions, SIAM J. Control Optim. 56 (2018), 1665-1691. Optimal control of a fractional tumor growth model 33

    work page 2018

  12. [12]

    Colli, G

    P. Colli, G. Gilardi, J. Sprekels, Well-posedness and regularity for a generalized fractional Cahn–Hilliard system, Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl. , to appear (see also preprint arXiv:1804.11290 [math.AP] (2018), pp . 1-36)

    work page arXiv 2018

  13. [14]

    Deep quench approximation and optimal control of general Cahn-Hilliard systems with fractional operators and double obstacle potentials

    P. Colli, G. Gilardi, J. Sprekels, Deep quench approximation and op timal control of general Cahn–Hilliard systems with fractional operators and doub le-obstacle poten- tials, preprint arXiv:1812.01675 [math.AP] (2018), pp. 1-32

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  14. [15]

    Colli, G

    P. Colli, G. Gilardi, J. Sprekels, Well-posedness and regularity for a fractional tumor growth model, Adv. Math. Sci. Appl. 28 (2019), 343-375

    work page 2019

  15. [16]

    Colli, G

    P. Colli, G. Gilardi, J. Sprekels, Longtime behavior for a generalize d Cahn–Hilliard system with fractional operators, preprint arXiv:1904.00931 [mat h.AP] (2019), pp. 1-18

    work page arXiv 1904

  16. [17]

    Colli, G

    P. Colli, G. Gilardi, J. Sprekels, Asymptotic analysis of a tumor gro wth model with fractional operators, in preparation

    work page

  17. [18]

    Conti, A

    M. Conti, A. Giorgini, The three-dimensional Cahn–Hilliard–Brinkm an system with unmatched densities, preprint hal-01559179 (2018), pp. 1-34

    work page 2018

  18. [19]

    Cristini, X

    V. Cristini, X. Li, J. S. Lowengrub, S. M. Wise, Nonlinear simulation s of solid tumor growth using a mixture model: invasion and branching. J. Math. Biol. 58 (2009), 723-763

    work page 2009

  19. [20]

    Multiscale Modeling of Cancer: An I ntegrated Experi- mental and Mathematical Modeling Approach

    V. Cristini, J. S. Lowengrub, “Multiscale Modeling of Cancer: An I ntegrated Experi- mental and Mathematical Modeling Approach”, Cambridge Universit y Press, Leiden (2010)

    work page 2010

  20. [21]

    M. Dai, E. Feireisl, E. Rocca, G. Schimperna, M. Schonbek, Analy sis of a diffuse interface model of multi-species tumor growth, Nonlinearity 30 (2017), 1639-1658

    work page 2017

  21. [22]

    Della Porta, A

    F. Della Porta, A. Giorgini, M. Grasselli, The nonlocal Cahn–Hilliard– Hele–Shaw system with logarithmic potential, Nonlinearity 31 (2018), 4851-4881

    work page 2018

  22. [23]

    Ebenbeck, H

    M. Ebenbeck, H. Garcke, Analysis of a Cahn–Hilliard–Brinkman mo del for tumour growth with chemotaxis, J. Differential Equations, 266 (2019), 5998-6036

    work page 2019

  23. [24]

    Ebenbeck, P

    M. Ebenbeck, P. Knopf, Optimal medication for tumors modeled by a Cahn–Hilliard– Brinkman equation, preprint arXiv:1811.07783 [math.AP] (2018), pp . 1-26

    work page arXiv 2018

  24. [25]

    Ebenbeck, P

    M. Ebenbeck, P. Knopf, Optimal control theory and advance d optimality conditions for a diffuse interface model of tumor growth, preprint arXiv:1903 .00333 [math.OC] (2019), pp. 1-34

    work page 1903

  25. [26]

    Frigeri, M

    S. Frigeri, M. Grasselli, E. Rocca, On a diffuse interface model of tumor growth, European J. Appl. Math. 26 (2015), 215-243. 34 Colli — Gilardi — Sprekels

    work page 2015

  26. [27]

    Frigeri, K

    S. Frigeri, K. F. Lam, E. Rocca, G. Schimperna, On a multi-specie s Cahn–Hilliard– Darcy tumor growth model with singular potentials, Commun. Math Sci. (16) (2018), 821-856

    work page 2018

  27. [28]

    C. G. Gal, On the strong-to-strong interaction case for doub ly nonlinear Cahn– Hilliard equations, Discrete Contin. Dyn. Syst. 37 (2017), 131-167

    work page 2017

  28. [29]

    C. G. Gal, Non-local Cahn–Hilliard equations with fractional dyna mic boundary conditions, European J. Appl. Math. 28 (2017), 736-788

    work page 2017

  29. [30]

    C. G. Gal, Doubly nonlinear Cahn–Hilliard equations, Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire35 (2018), 357-392

    work page 2018

  30. [31]

    Garcke, K

    H. Garcke, K. F. Lam, Global weak solutions and asymptotic limits of a Cahn– Hilliard–Darcy system modelling tumour growth, AIMS Mathematics 1 (2016), 318- 360

    work page 2016

  31. [32]

    Garcke, K

    H. Garcke, K. F. Lam, Well-posedness of a Cahn–Hilliard system m odelling tumour growth with chemotaxis and active transport, European. J. Appl. Math. 28 (2017), 284-316

    work page 2017

  32. [33]

    Garcke, K

    H. Garcke, K. F. Lam, Analysis of a Cahn–Hilliard system with non– zero Dirichlet conditions modeling tumor growth with chemotaxis, Discrete Contin. Dyn. Syst. 37 (2017), 4277-4308

    work page 2017

  33. [34]

    Trends on Applications of Mat hematics to Me- chanics

    H. Garcke, K. F. Lam, On a Cahn–Hilliard–Darcy system for tumo ur growth with solution dependent source terms, in “Trends on Applications of Mat hematics to Me- chanics”, E. Rocca, U. Stefanelli, L. Truskinovski, A. Visintin (ed.), Springer INdAM Series 27, Springer, Cham, 2018, pp. 243-264

    work page 2018

  34. [35]

    Garcke, K

    H. Garcke, K. F. Lam, R. N¨ urnberg, E. Sitka, A multiphase Cah n–Hilliard–Darcy model for tumour growth with necrosis, Math. Models Methods Appl. Sci. 28 (2018), 525-577

    work page 2018

  35. [36]

    Garcke, K

    H. Garcke, K. F. Lam, E. Rocca, Optimal control of treatmen t time in a diffuse interface model of tumor growth, Appl. Math. Optim. 78 (2018), 495-544

    work page 2018

  36. [37]

    Garcke, K

    H. Garcke, K. F. Lam, E. Sitka, V. Styles, A Cahn–Hilliard–Darcy model for tumour growth with chemotaxis and active transport, Math. Models Methods Appl. Sci. 26 (2016), 1095-1148

    work page 2016

  37. [38]

    Hawkins-Daarud, S

    A. Hawkins-Daarud, S. Prudhomme, K. G. van der Zee, J. T. Od en, Bayesian calibra- tion, validation, and uncertainty quantification of diffuse interface models of tumor growth. J. Math. Biol. 67 (2013), 1457-1485

    work page 2013

  38. [39]

    Hawkins-Daruud, K

    A. Hawkins-Daruud, K. G. van der Zee, J. T. Oden, Numerical s imulation of a ther- modynamically consistent four-species tumor growth model, Int. J. Numer. Math. Biomed. Engng. 28 (2012), 3-24

    work page 2012

  39. [40]

    Quelques M´ ethodes de R´ esolution des Probl` em es aux Limites non Lin´ eaires

    J.-L. Lions, “Quelques M´ ethodes de R´ esolution des Probl` em es aux Limites non Lin´ eaires”, Dunod, Gauthier-Villars, Paris, 1969. Optimal control of a fractional tumor growth model 35

    work page 1969

  40. [41]

    Miranville, E

    A. Miranville, E. Rocca, G. Schimperna, On the long time behavior o f a tumor growth model, J. Differential Equations 267 (2019), 2616-2642

    work page 2019

  41. [42]

    J. T. Oden, A. Hawkins, S. Prudhomme, General diffuse-interf ace theories and an approach to predictive tumor growth modeling, Math. Models Methods Appl. Sci. 20 (2010), 477-517

    work page 2010

  42. [43]

    Signori, Optimal distributed control of an extended model o f tumor growth with logarithmic potential, Appl

    A. Signori, Optimal distributed control of an extended model o f tumor growth with logarithmic potential, Appl. Math. Optim. (2018), https://doi.org/10.1007/s00245- 018-9538-1

    work page doi:10.1007/s00245- 2018

  43. [44]

    Signori, Optimality conditions for an extended tumor growth m odel with double obstacle potential via deep quench approach, preprint arXiv:1811 .08626 [math.AP] (2018), pp

    A. Signori, Optimality conditions for an extended tumor growth m odel with double obstacle potential via deep quench approach, preprint arXiv:1811 .08626 [math.AP] (2018), pp. 1-25

    work page 2018

  44. [45]

    Signori, Optimal treatment for a phase field system of Cahn– Hilliard type mod- eling tumor growth by asymptotic scheme, preprint arXiv:1902.0107 9 [math.AP] (2019), pp

    A. Signori, Optimal treatment for a phase field system of Cahn– Hilliard type mod- eling tumor growth by asymptotic scheme, preprint arXiv:1902.0107 9 [math.AP] (2019), pp. 1-28

    work page arXiv 1902

  45. [46]

    Signori, Vanishing parameter for an optimal control problem modeling tumor growth, preprint arXiv:1903.04930 [math.AP] (2019), pp

    A. Signori, Vanishing parameter for an optimal control problem modeling tumor growth, preprint arXiv:1903.04930 [math.AP] (2019), pp. 1-22

    work page arXiv 1903

  46. [47]

    Sprekels, H

    J. Sprekels, H. Wu, Optimal distributed control of a Cahn–Hilliar d–Darcy system with mass sources, Appl. Math. Optim. (2019), https://doi.org/10.1007/s00245-019- 09555-4

    work page doi:10.1007/s00245-019- 2019

  47. [48]

    S. M. Wise, J. S. Lowengrub, H. B. Frieboes, V. Cristini, Three- dimensional mul- tispecies nonlinear tumor growth - I: model and numerical method, J. Theor. Biol. 253 (2008), 524-543

    work page 2008

  48. [49]

    X. Wu, G. J. van Zwieten, K. G. van der Zee, Stabilized second-o rder splitting schemes for Cahn–Hilliard models with applications to diffuse-interfac e tumor-growth models, Int. J. Numer. Meth. Biomed. Engng. 30 (2014), 180-203

    work page 2014