Pith. sign in

REVIEW 2 major objections 5 minor 2 cited by

Chemotaxis can turn a stable uniform tumour–immune balance into finite-wavelength spatial patterns once a critical sensitivity is crossed.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 23:46 UTC pith:NYXSWHHP

load-bearing objection Clean, carefully scoped analysis of a minimal tumour–immune–chemokine RDC system: the chemotaxis-as-wavenumber-amplified-coupling claim holds, numerics verify it experiment-by-experiment, and the only real soft spot is the already-flagged lack of global L∞ theory. the 2 major comments →

arxiv 2607.03813 v1 pith:NYXSWHHP submitted 2026-07-04 math.AP cs.NAmath.NA

Computational Oncology of Chemotaxis-Driven Tumour--Immune Spatial Patterning and Stability

Zonghao Liu , Jiguang Yu , Lei Su , Louis Shuo Wang , Yang Du , Jingfeng Liu This is my paper
classification math.AP cs.NAmath.NA MSC 35K5792C1765M0835B3692C50
keywords tumour–immune dynamicschemotaxisreaction–diffusionlinear stabilitydispersion relationpattern formationfinite-volume methodpositivity preservation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper builds a minimal three-field continuum model of solid-tumour tissue: tumour density, immune effector density, and a chemokine signal that both recruits immune cells and steers their directed migration. After nondimensionalization the authors prove that densities stay nonnegative, tumour density stays bounded by carrying capacity, and immune and chemokine masses stay controlled. They show that a tumour-free state is stable precisely when baseline immune supply exceeds a simple threshold, and that any coexistence state reduces to a single scalar algebraic equation. Linearizing about coexistence produces a mode-by-mode stability matrix in which chemotaxis enters only as a term that grows with wavenumber; above a critical chemotactic sensitivity this term drives finite-wavelength instability while leaving the uniform mode untouched. A structure-preserving finite-volume scheme then verifies the thresholds, the dominant unstable modes, positivity, and residual consistency. The result supplies a concrete, checkable mechanism for how chemokine-guided immune migration can convert a homogeneous “hot” coexistence into spatially heterogeneous infiltration or exclusion patterns.

Core claim

About a homogeneous tumour–immune coexistence equilibrium, the chemotactic sensitivity appears in the mode-wise stability matrix solely as the wavenumber-amplified coupling +ξ v* μ_k. This coupling leaves the spatially uniform mode unaffected but produces a finite-wavelength instability once ξ exceeds a critical value ξ_c defined by the first zero of the dispersion relation ω_k(ξ) = max Re λ_j(M_k).

What carries the argument

The mode-wise stability matrix M_k obtained by linearizing about coexistence and projecting onto Neumann Laplacian eigenmodes; chemotaxis enters only through the (2,3) entry as +ξ v* μ_k, which is the sole route to finite-wavelength destabilization.

Load-bearing premise

The analysis and long-time numerics assume classical solutions remain bounded on the time interval of interest, yet the paper only proves local existence, positivity, a tumour density bound, and mass control—not global pointwise bounds for the full chemotaxis system.

What would settle it

Compute the eigenvalues of M_k for the paper’s baseline parameters and show that the first mode that crosses Re λ = 0 is not a finite-k mode (or that the measured growth rates of the corresponding discrete-cosine modes in the nonlinear simulation disagree with the predicted dispersion relation).

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper formulates a minimal three-component reaction–diffusion–chemotaxis model for tumour–immune–chemokine interactions on a bounded domain with no-flux boundaries. After nondimensionalization it proves local classical solvability, positivity, a uniform tumour L^∞ bound, and L^1 mass control for immune and chemokine densities; identifies the tumour-free equilibrium and the immune-control threshold σ₀>δ; reduces homogeneous coexistence to a scalar equation F(u)=0; and derives a mode-wise dispersion relation in which chemotaxis enters the stability matrix only as the wavenumber-amplified coupling +ξ v* μ_k, producing finite-wavelength instability above a critical sensitivity ξ_c. A conservative finite-volume scheme with upwind chemotactic flux is used to verify thresholds, dominant modes, sensitivity maps, positivity, grid convergence, and residual consistency in six targeted experiments.

Significance. If the analysis and numerics hold as stated, the paper supplies a clean, mechanistically interpretable link between chemokine-mediated chemotaxis and the transition from homogeneous tumour–immune coexistence to finite-wavelength spatial heterogeneity—a mathematical caricature of immune infiltration versus exclusion. Strengths that should be credited include: (i) honest scoping of what is proved (local classical solutions, tumour bound, mass estimates) versus what is not (global L^∞ for the full chemotaxis system); (ii) an explicit mode-wise matrix in which the chemotactic contribution is isolated and vanishes for the homogeneous mode; (iii) a structure-preserving FV scheme with positivity, mass-balance, and post-hoc residual diagnostics; and (iv) experiments organized as falsifiable checks of each theoretical claim rather than as illustrations. The contribution is incremental relative to the Keller–Segel and tumour–immune literature, but the combination of threshold analysis, dispersion relation, and certified numerics is useful for computational oncology.

major comments (2)
  1. [§6, Appendix B.6, Experiment 4] Appendix B.6 states that for the Section 7 parameter set one has A_k t − p s < 0 in the relevant band, so the constant-term (stationary/Turing) mechanism is inactive and onset is a Hopf-type crossing of a1 a2 = a3. Experiment 4, however, reports saturated spatial fields that appear stationary (immune aggregation peaks at chemokine maxima, amplitude plateaus). If the linear onset is genuinely Hopf, the weakly nonlinear regime should exhibit temporal oscillation at least near threshold. Please either (i) document the imaginary part of the critical eigenvalues and show time series near ξ_c that confirm or rule out oscillation, or (ii) reclassify the crossing for the reported parameters and reconcile the wording in §6–§7 with the observed patterns. This is load-bearing for the claim of chemotaxis-driven pattern formation.
  2. [Abstract; Remark 4.4; §7 Experiments 3–5; §8] Remark 4.4 and Appendix A.6 correctly state that global L^∞ bounds are not proved and that long-time conclusions rest on classical solutions on their interval of existence (or on numerically verified boundedness). The abstract, §7 introduction, and conclusion still present the homogeneous-to-heterogeneous transition and ξ_c as established model properties without always restating this condition. Please make the conditional status of all long-time pattern and stability conclusions explicit in the abstract and in the statements surrounding Experiments 3–5 (e.g., “for classical solutions remaining bounded on [0,T]”), so that the central claim is not read as unconditional global dynamics.
minor comments (5)
  1. [Figure 2] Figure 2 labels use “s0>delta” / “s0<delta” while the text uses σ₀ and δ. Align notation with the nondimensional symbols of Table 2.
  2. [§6.1, Eq. (6.1)] In (6.1) and (B.1) the (2,3) entry is written as σ₁/(1+w*)² + ξ v* μ_k; a short inline reminder that this comes from −ξ v* Δw̃ = ξ v* μ_k W_k would help readers who skip Appendix B.
  3. [Table 3; Eq. (7.1); Experiment 2] Table 3 and the coexistence values (7.1) are given to four decimals; state the root-finding tolerance used for F(u)=0 so that Experiment 2’s 10⁻¹¹ agreement is reproducible.
  4. [§3.5; §7] The 2D extension of the scheme is described in §3.5 but all patterning experiments are 1D. A single 2D snapshot (or an explicit statement that 2D is left for future work) would match the domain dimension advertised in §2.1.
  5. [Keywords; Figure 1; Table 2] Minor typesetting: “finite-volumemethod” and similar missing spaces in the keywords/abstract; “T umour” in Figure 1 caption; consistent use of ell vs script-l for the degradation parameter.

Circularity Check

0 steps flagged

No significant circularity: thresholds, coexistence reduction, and chemotaxis dispersion relation are derived from the PDE by standard linearization and Routh–Hurwitz, with numerics as independent verification rather than fitted prediction.

full rationale

The paper is a self-contained reaction–diffusion–chemotaxis analysis plus structure-preserving numerics. The tumour-free threshold σ₀>δ follows from the linearized growth rate 1−σ₀/δ of the homogeneous tumour mode (Theorem 5.1); coexistence reduces to the scalar F(u)=0 obtained by substituting v*=1−u* and w*=u*[α+γ(1−u*)]/ℓ into the immune equilibrium (Proposition 5.3); and the mode-wise matrix M_k (Eq. 6.1 / B.1) is the standard projection of the linearized system onto Neumann eigenmodes, with chemotaxis entering only as the wavenumber-amplified coupling +ξ v* μ_k. The critical sensitivity ξ_c is defined from the dispersion relation ω_k(ξ)=max_j Re λ_j(M_k) and is not fitted to data. Baseline parameters (Table 3) are chosen for illustration; Experiments 1–6 verify the analytical thresholds, dominant modes, and scheme properties against those same equations, with residual/positivity/convergence diagnostics. Self-citations to Wang/Yu/Liu appear only in related-work and hybrid PDE–ABM context and are not load-bearing for §§4–6. The acknowledged gap (local well-posedness without global L^∞ bounds) is a correctness/scope limitation, not circularity. No step reduces a claimed prediction to its own definition or to a self-citation chain.

Axiom & Free-Parameter Ledger

2 free parameters · 6 axioms · 1 invented entities

The paper is a mechanistic continuum model plus classical PDE analysis and structure-preserving numerics. Load-bearing inputs are standard parabolic theory, comparison principles, Routh–Hurwitz, and a suite of phenomenological reaction/chemotaxis terms. Free parameters are the nondimensional coefficients chosen for numerics; none are fitted to clinical data. The three fields and their couplings are modeling constructs, not new physical entities with independent empirical handles beyond the general chemokine biology cited.

free parameters (2)
  • Baseline nondimensional set (d1,d2,d3,σ0,σ1,δ,β,α,γ,ℓ,L) in Table 3 = d1=0.01,d2=0.10,d3=1.00,σ0=0.10,σ1=0.50,δ=0.30,β=0.40,α=0.50,γ=0.50,ℓ=1.00,L=10
    Hand-chosen for numerical experiments so that a unique interior coexistence equilibrium exists and is homogeneously stable; used to compute ξ_c≈12.47 and all pattern figures. Not fitted to data.
  • Chemotactic sensitivity ξ (and multiples of ξ_c) = ξ_c≈12.47 for baseline set
    Primary bifurcation parameter scanned in Experiments 3–5; critical value is computed from the model, not measured.
axioms (6)
  • standard math Local classical solvability of normally parabolic quasilinear systems (Amann-type theory) for triangular principal part with positive diffusivities
    Invoked for Theorem 4.1 and Appendix A.1; standard background.
  • standard math Quasi-positivity of reaction terms plus no-flux boundaries imply invariance of the nonnegative cone (parabolic comparison / negative-part argument)
    Appendix A.2; used for biological meaningfulness of densities.
  • standard math Routh–Hurwitz criterion for monic cubics governs modal stability of M_k
    Section 6 and Appendix B.4; classical linear algebra.
  • domain assumption Logistic tumour growth, mass-action immune killing, Michaelis–Menten chemokine-enhanced recruitment, linear immune loss, and linear chemokine production/degradation adequately represent the dominant local interactions
    Model formulation (2.1)/(2.3); phenomenological choices standard in mathematical oncology but not derived from first principles.
  • domain assumption No-flux tissue boundaries and, for linear stability, componentwise Neumann conditions so the Neumann Laplacian eigenbasis applies
    Sections 2.2 and 6; stronger Neumann condition is analytical convenience that implies the physical no-flux condition when ∇A·n=0.
  • ad hoc to paper Classical solutions remain in a regime where linearization about coexistence and long-time numerics are meaningful (global L^∞ not proved)
    Remark 4.4 and Appendix A.6 explicitly decline full global boundedness; subsequent stability and pattern claims assume solutions stay classical/bounded.
invented entities (1)
  • Minimal three-component tumour–immune–chemokine RDC system (2.1)/(2.3) no independent evidence
    purpose: Couples logistic growth, immune killing, chemokine production by T and ET, recruitment, and immune chemotaxis in one analytically tractable continuum model
    A modeling construct assembled from known mechanisms; not a new physical particle or force. Independent biological motivation is cited (chemokine guidance, hot/cold tumours) but the specific equations are postulated.

pith-pipeline@v1.1.0-grok45 · 31718 in / 3979 out tokens · 36014 ms · 2026-07-11T23:46:27.195792+00:00 · methodology

0 comments
read the original abstract

Spatial tumour--immune heterogeneity is a key feature of solid-tumour progression, immune infiltration, and immune exclusion. We develop a computational oncology model in which tumour cells, immune effector cells, and a chemokine signal interact through a reaction--diffusion--chemotaxis system on a bounded tissue domain with no-flux boundaries. Chemokine is produced by tumour cells and tumour--immune contact, recruits immune cells, and guides chemotactic migration. After nondimensionalization, we establish positivity, a tumour-density bound, and immune/chemokine mass estimates. We identify the tumour-free equilibrium, derive the immune-control threshold $\sigma_0>\delta$, and reduce coexistence to a scalar equation. Linear stability analysis about coexistence yields a mode-wise dispersion relation in which chemotaxis appears as a wavenumber amplified coupling, producing finite-wavelength instability above a critical sensitivity. A conservative finite-volume scheme with upwind chemotactic flux verifies the thresholds, dominant unstable modes, sensitivity maps, positivity, convergence, and residual consistency.

Figures

Figures reproduced from arXiv: 2607.03813 by Jiguang Yu, Jingfeng Liu, Lei Su, Louis Shuo Wang, Yang Du, Zonghao Liu.

Figure 1
Figure 1. Figure 1: Tumour–immune–chemokine feedback loop used in the computational oncology model. Tumour cells 𝑇 produce chemokine 𝐴, while tumour–immune contacts 𝐸𝑇 provide an additional chemokine source. The chemokine field recruits immune effector cells 𝐸, and its spatial gradient ∇𝐴 drives chemotactic immune migration. Immune effectors suppress tumour cells through local killing. This feedback structure is represented m… view at source ↗
Figure 2
Figure 2. Figure 2: Experiment 1. (a) Tumour amplitude ‖𝑢(𝑡)‖∞ for the decay regime 𝜎0 > 𝛿 and the invasion regime 𝜎0 < 𝛿; (b) tumour mass ∫Ω 𝑢 𝑑𝑥; (c) measured versus predicted growth rate 1 − 𝜎0∕𝛿 (dashed line is identity). The measured rates agree with Theorem 5.1. 17 [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Experiment 2. (a) Convergence of the solution to the analytical coexistence equilibrium in 𝐿2 ; (b) final numerical fields (solid) versus the analytical values (𝑢 ∗ , 𝑣∗ , 𝑤∗ ) of (7.1) (dashed). 19 [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Experiment 3. (a) Dispersion relation max𝑗 Re 𝜆𝑗 (𝑀𝑘 ) versus wavenumber index 𝑘 for 𝜉 ∈ {0, 0.6 𝜉𝑐 , 1.6 𝜉𝑐 } (curves: continuous 𝜇; markers: admissible Neumann modes). (b) Linear-theory modal growth rates (line) versus rates measured from the nonlinear simulation via the discrete cosine transform (squares: active band; crosses: damped modes at the noise floor) for 𝜉 = 1.6 𝜉𝑐 . 21 [PITH_FULL_IMAGE:figure… view at source ↗
Figure 5
Figure 5. Figure 5: Experiment 4. Top row: final tumour (𝑢) and immune (𝑣) fields for subcritical, near-critical, and supercritical chemotactic sensitivity. Bottom row: corresponding spatial amplitude max𝑥 𝑢 − min𝑥 𝑢 versus time. The transition homogeneous → onset → heterogeneity occurs across 𝜉𝑐 ≈ 12.5. 23 [PITH_FULL_IMAGE:figures/full_fig_p023_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Experiment 5. Stability maps of max𝑘≥1 max𝑗 Re 𝜆𝑗 (𝑀𝑘 ) in the planes (𝜉, 𝜎1 ), (𝜉, 𝛿), and (𝜉, 𝑑2 ). Blue: homogeneous coexistence stable; red: chemotaxis-driven patterning; black curve: onset (Re 𝜆 = 0). The equilibrium is recomputed at every grid point. 25 [PITH_FULL_IMAGE:figures/full_fig_p025_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Experiment 6. (a) Grid convergence ‖𝑈ℎ − 𝑈ℎ∕2‖𝐿2 versus 𝑁 (observed order 1.92; dashed (ℎ 2 ) reference). (b) Post-hoc partial differential equation (PDE) residual norms ‖𝑅𝑢‖, ‖𝑅𝑣‖, ‖𝑅𝑤‖ versus time on the finest grid. 27 [PITH_FULL_IMAGE:figures/full_fig_p027_7.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Reflected Optimal Stopping with a Max-Type Payoff: Measure-Valued Stopping Gains and Killed Resolvent Representation

    math.AP 2026-07 conditional novelty 6.0

    For reflected diffusion optimal stopping with G=x1∨αx2, the gain Γ is a signed measure with singular diagonal part −(n⊤an)/(2√(1+α²))σΔ, and V=G−R_C^r Γ (killed at first entry to the stopping set).

  2. Multiplier Sensitivity in Isoperimetric Optimal Control

    math.AP 2026-07 accept novelty 6.0

    Under stated convexity and attainability conditions, the isoperimetric Pontryagin multiplier is the negative supergradient of the constrained value function, and this sensitivity can survive when the modified control-...

Reference graph

Works this paper leans on

75 extracted references · cited by 2 Pith papers

  1. [1]

    Measuring and modelling tumour heterogeneity across scales

    Garrett F Beeghly, Adrian A Shimpi, Robert N Riter, and Claudia Fischbach. Measuring and modelling tumour heterogeneity across scales. Nature Reviews Bioengineering, 1(10):712–730, 2023

    2023

  2. [2]

    Computational frameworks for modelling cancer across scales.Clinical and Translational Discovery, 5(6):e70093, 2025

    Celine Desoyer, Daniel Loibner, Dagmar Brislinger, and Christian Baumgartner. Computational frameworks for modelling cancer across scales.Clinical and Translational Discovery, 5(6):e70093, 2025

    2025

  3. [3]

    Molecular basis of cancer.Pathological Basis of Oral and Maxillofacial Diseases, pages 415–428, 2025

    Lalima Tiwari and Omar Kujan. Molecular basis of cancer.Pathological Basis of Oral and Maxillofacial Diseases, pages 415–428, 2025

    2025

  4. [4]

    The role of the cell–cell interactions in cancer progression.Journal of cellular and molecular medicine, 19(2):283–296, 2015

    Katarzyna Kamińska, Cezary Szczylik, Zofia F Bielecka, Ewa Bartnik, Camillo Porta, Fei Lian, and Anna M Czarnecka. The role of the cell–cell interactions in cancer progression.Journal of cellular and molecular medicine, 19(2):283–296, 2015

    2015

  5. [5]

    The dynamic tumor ecm: Biophysical cues, cellular crosstalk, and disease progression.Current Opinion in Biomedical Engineering, page 100652, 2026

    Omkar Joshi, Hellyeh Hamidi, Mathilde Mathieu, and Johanna Ivaska. The dynamic tumor ecm: Biophysical cues, cellular crosstalk, and disease progression.Current Opinion in Biomedical Engineering, page 100652, 2026

    2026

  6. [6]

    Extracellular matrix in cancer progression and therapy.Medical Review, 2(2):125–139, 2022

    Xiuxiu He, Byoungkoo Lee, and Yi Jiang. Extracellular matrix in cancer progression and therapy.Medical Review, 2(2):125–139, 2022

    2022

  7. [7]

    Giuseppe Giuliani, Dhruv Chavan, Rushil Gaddam, William C Stewart, Zihai Li, Ciriyam Jayaprakash, and Jayajit Das. Immune and non- immune cell fencing of tumor cells is a widespread and functionally relevant spatial pattern in solid cancers.Computational and Structural Biotechnology Journal, 2025

    2025

  8. [8]

    Inferenceonspatialheterogeneityintumormicroenvironmentusing spatial transcriptomics data.Computational and systems oncology, 2(3):e21043, 2022

    AntaraBiswas,BasselGhaddar,GregoryRiedlinger,andSubhajyotiDe. Inferenceonspatialheterogeneityintumormicroenvironmentusing spatial transcriptomics data.Computational and systems oncology, 2(3):e21043, 2022

    2022

  9. [9]

    Reaction-diffusion model as a framework for understanding biological pattern formation.science, 329(5999):1616–1620, 2010

    Shigeru Kondo and Takashi Miura. Reaction-diffusion model as a framework for understanding biological pattern formation.science, 329(5999):1616–1620, 2010

    2010

  10. [10]

    How chemokines organize the tumour microenvironment.Nature Reviews Cancer, 24(1):28–50, 2024

    Thorsten R Mempel, Julia K Lill, and Lukas M Altenburger. How chemokines organize the tumour microenvironment.Nature Reviews Cancer, 24(1):28–50, 2024

    2024

  11. [11]

    Mathematicalmodellingofthespatio-temporalresponseofcytotoxic t-lymphocytes to a solid tumour.Mathematical Medicine and Biology, 21(1):1–34, 2004

    AnastasiosMatzavinos,MarkAJChaplain,andVladimirAKuznetsov. Mathematicalmodellingofthespatio-temporalresponseofcytotoxic t-lymphocytes to a solid tumour.Mathematical Medicine and Biology, 21(1):1–34, 2004

    2004

  12. [12]

    Globalsmoothsolutionsinachemotaxissystemmodelingimmuneresponsetoasolidtumor.Proceedings of the American Mathematical Society, 152(10):4325–4341, 2024

    YoushanTaoandMichaelWinkler. Globalsmoothsolutionsinachemotaxissystemmodelingimmuneresponsetoasolidtumor.Proceedings of the American Mathematical Society, 152(10):4325–4341, 2024

    2024

  13. [13]

    LouisShuoWang,JiguangYu,YeLiang,andJilinZhang.Ellipticcriticalityversusvolterramemoryinindirectchemotaxiscascades.Transport Phenomena, 1(2):20260061, 2026. 34

    2026

  14. [14]

    Globalboundednessofanndchemotactictumorimmuneevasionsystem.DiscreteandContinuousDynamical Systems-B, pages 0–0, 2026

    KaiqiangLiandYingyingLi. Globalboundednessofanndchemotactictumorimmuneevasionsystem.DiscreteandContinuousDynamical Systems-B, pages 0–0, 2026

    2026

  15. [15]

    Reaction, diffusion and chemotaxis in wave propagation.Discrete Contin

    Shangbing Ai, Wenzhang Huang, and Zhi-An Wang. Reaction, diffusion and chemotaxis in wave propagation.Discrete Contin. Dyn. Syst. Ser. B, 20(1):1–21, 2015

    2015

  16. [16]

    Springer, 2022

    Yuanyuan Ke, Jing Li, and Yifu Wang.Analysis of reaction-diffusion models with the taxis mechanism. Springer, 2022

    2022

  17. [17]

    Chemotaxisandreactionsinbiology.JournaloftheEuropeanMathematical Society, 25(7):2641–2696, 2022

    AlexanderKiselev,FedorNazarov,LenyaRyzhik,andYaoYao. Chemotaxisandreactionsinbiology.JournaloftheEuropeanMathematical Society, 25(7):2641–2696, 2022

    2022

  18. [18]

    Integrative mathematical oncology.Nature Reviews Cancer, 8(3):227–234, 2008

    Alexander RA Anderson and Vito Quaranta. Integrative mathematical oncology.Nature Reviews Cancer, 8(3):227–234, 2008

    2008

  19. [19]

    The mathematics of cancer: integrating quantitative models.Nature Reviews Cancer, 15(12):730–745, 2015

    Philipp M Altrock, Lin L Liu, and Franziska Michor. The mathematics of cancer: integrating quantitative models.Nature Reviews Cancer, 15(12):730–745, 2015

    2015

  20. [20]

    From microscopic damage to macroscopic games: a dimensionality reduction of stem cell homeostasis.Transport Phenomena, 1(2):20260037, 2026

    Jiguang Yu, Louis Shuo Wang, Shihan Ban, and Ye Liang. From microscopic damage to macroscopic games: a dimensionality reduction of stem cell homeostasis.Transport Phenomena, 1(2):20260037, 2026

    2026

  21. [21]

    Computationalapproaches to modelling and optimizing cancer treatment.Nature Reviews Bioengineering, 1(10):695–711, 2023

    ThomasOMcDonald,Yu-ChenCheng,ChristopherGraser,PhillipBNicol,DanielTemko,andFranziskaMichor. Computationalapproaches to modelling and optimizing cancer treatment.Nature Reviews Bioengineering, 1(10):695–711, 2023

    2023

  22. [22]

    ST Gopukumar, Dyumn Dwivedi, Mohamed Rahamathulla, Mohammed Muqtader Ahmed, Tanveen Kaur Soni, Praveen Ganesh Natarajan, S Ramkanth, Arpita Mitra, and Uddalak Das. Multiscale predictive cellular modeling: integrating hypothesis grammars, digital twins, and multi-omics for in silico oncology and precision theranostics.Functional & Integrative Genomics, 26(...

    2026

  23. [23]

    The 2019 mathematical oncology roadmap.Physical biology, 16(4):041005, 2019

    RussellCRockne,AndreaHawkins-Daarud,KristinRSwanson,JamesPSluka,JamesAGlazier,PaulMacklin,DavidAHormuth,AngelaM Jarrett, Ernesto ABF Lima, J Tinsley Oden, et al. The 2019 mathematical oncology roadmap.Physical biology, 16(4):041005, 2019

    2019

  24. [24]

    A review of mechanistic learning in mathematical oncology.Frontiers in Immunology, 15:1363144, 2024

    John Metzcar, Catherine R Jutzeler, Paul Macklin, Alvaro Köhn-Luque, and Sarah C Brüningk. A review of mechanistic learning in mathematical oncology.Frontiers in Immunology, 15:1363144, 2024

    2024

  25. [25]

    Optimal harvesting for nonlinear size-structured populations with nonlocal environmental feedback.Mathematics, 14(11):2025, 2026

    Cai Jie, Xiaoyang Chen, Gu Longfei, Chen Jiayao, Chu Nuo, Louis Shuo Wang, Liang Ye, and Yu Jiguang. Optimal harvesting for nonlinear size-structured populations with nonlocal environmental feedback.Mathematics, 14(11):2025, 2026

    2025

  26. [26]

    YunduoLan,Sung-YoungShin,andLanKNguyen.Fromshallowtodeep:theevolutionofmachinelearningandmechanisticmodelintegration in cancer research.Current Opinion in Systems Biology, 40:100541, 2025

    2025

  27. [27]

    Radiation-induced lymphopenia: From mathematical modeling towards mechanistic learning.International Journal of Radiation Oncology* Biology* Physics, 2025

    FrançoisdeKermenguy,DaphnéMorel,MohammedEl-Aichi,DominiqueBarbolosi,EricDeutsch,andCharlotteRobert. Radiation-induced lymphopenia: From mathematical modeling towards mechanistic learning.International Journal of Radiation Oncology* Biology* Physics, 2025

    2025

  28. [28]

    Mechanistic learning with guided diffusion models to predict spatio-temporal brain tumor growth

    Daria Laslo, Efthymios Georgiou, Marius George Linguraru, Andreas M Rauschecker, Sabine Müller, Catherine R Jutzeler, and Sarah Brüningk. Mechanistic learning with guided diffusion models to predict spatio-temporal brain tumor growth. InInternational Workshop on Learning with Longitudinal Medical Images and Data, pages 68–79. Springer, 2025

    2025

  29. [29]

    Multi-strategy hybrid improved intelligent algorithm for solving uav-mtsp.Information Technology and Control, 54(2):413–438, 2025

    Zixin Wang, Danqing Wang, and Jiguang Yu. Multi-strategy hybrid improved intelligent algorithm for solving uav-mtsp.Information Technology and Control, 54(2):413–438, 2025

    2025

  30. [30]

    Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis.Bulletin of mathematical biology, 56(2):295–321, 1994

    Vladimir A Kuznetsov, Iliya A Makalkin, Mark A Taylor, and Alan S Perelson. Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis.Bulletin of mathematical biology, 56(2):295–321, 1994

    1994

  31. [31]

    Modeling immunotherapy of the tumor–immune interaction.Journal of mathematical biology, 37(3):235–252, 1998

    Denise Kirschner and John Carl Panetta. Modeling immunotherapy of the tumor–immune interaction.Journal of mathematical biology, 37(3):235–252, 1998

    1998

  32. [32]

    Amathematicaltumormodelwithimmuneresistanceanddrugtherapy:anoptimalcontrolapproach

    LisetteGDePillisandAmiRadunskaya. Amathematicaltumormodelwithimmuneresistanceanddrugtherapy:anoptimalcontrolapproach. Computational and Mathematical Methods in Medicine, 3(2):79–100, 2001

    2001

  33. [33]

    Reza Mehdizadeh, Seyed Peyman Shariatpanahi, Bahram Goliaei, and Curzio Rüegg. Targeting myeloid-derived suppressor cells in combinationwithtumorcellvaccinationpredictsanti-tumorimmunityandbreastcancerdormancy:Aninsilicoexperiment.Scientificreports, 13(1):5875, 2023

    2023

  34. [34]

    Beyond diagonal noise: A better predator-prey modeling framework with cross-covariance.PLoS One, 21(5):e0350127, 2026

    Jiguang Yu and Louis Shuo Wang. Beyond diagonal noise: A better predator-prey modeling framework with cross-covariance.PLoS One, 21(5):e0350127, 2026

    2026

  35. [35]

    Impactofmildhyperthermiaontumor-immune dynamics explored through mathematical modeling.Scientific Reports, 2025

    AndrésSebastiánBarrera-León,LeoncioRodríguez-Quiñones,andAdrianaMDeMendoza. Impactofmildhyperthermiaontumor-immune dynamics explored through mathematical modeling.Scientific Reports, 2025

    2025

  36. [36]

    Rollingpredictionmodelofclosingpricebasedoneemddatanoisereductionandhgs-delm

    YuanshengGao,LeiLi,andJiguangYu. Rollingpredictionmodelofclosingpricebasedoneemddatanoisereductionandhgs-delm. In2022 International Conference on Data Analytics, Computing and Artificial Intelligence (ICDACAI), pages 255–260. IEEE, 2022

    2022

  37. [37]

    Interactions between the immune system and cancer: a brief review of non-spatial mathematical models.Bulletin of mathematical biology, 73(1):2–32, 2011

    Raluca Eftimie, Jonathan L Bramson, and David JD Earn. Interactions between the immune system and cancer: a brief review of non-spatial mathematical models.Bulletin of mathematical biology, 73(1):2–32, 2011

    2011

  38. [38]

    Rigorous analysis of a nonlocal transport–renewal system for physiologically structured populations.Mathematical Methods in the Applied Sciences, 2026

    Jiguang Yu, Louis Shuo Wang, and Ye Liang. Rigorous analysis of a nonlocal transport–renewal system for physiologically structured populations.Mathematical Methods in the Applied Sciences, 2026

    2026

  39. [39]

    Hybrid model of tumor growth, angiogenesis and immune response yields strategies to improve antiangiogenic therapy.npj Biological Physics and Mechanics, 1(1):4, 2024

    Andreas G Hadjigeorgiou and Triantafyllos Stylianopoulos. Hybrid model of tumor growth, angiogenesis and immune response yields strategies to improve antiangiogenic therapy.npj Biological Physics and Mechanics, 1(1):4, 2024

    2024

  40. [40]

    Continuous and discrete mathematical models of tumor-induced angiogenesis.Bulletin of mathematical biology, 60(5):857–899, 1998

    Alexander RA Anderson and Mark AJ Chaplain. Continuous and discrete mathematical models of tumor-induced angiogenesis.Bulletin of mathematical biology, 60(5):857–899, 1998

    1998

  41. [41]

    Animage-informedcahn–hilliardkeller–segelmultiphasefieldmodelfortumorgrowth with angiogenesis.Applied Mathematics and Computation, 445:127834, 2023

    AbramoAgosti,AGiottaLucifero,andSabinoLuzzi. Animage-informedcahn–hilliardkeller–segelmultiphasefieldmodelfortumorgrowth with angiogenesis.Applied Mathematics and Computation, 445:127834, 2023

    2023

  42. [42]

    Thebreakdownoflinearquasi-cycles:Demographicnoiseandabsorbingboundaries in finite predator–prey systems.Electronic Research Archive, 34(6):4248–4289, 2026

    LouisShuoWang,JiguangYu,YeLiang,andJilinZhang. Thebreakdownoflinearquasi-cycles:Demographicnoiseandabsorbingboundaries in finite predator–prey systems.Electronic Research Archive, 34(6):4248–4289, 2026

    2026

  43. [43]

    Mathematicalmodelsofavasculartumorgrowth.SIAMreview,49(2):179–208,2007

    TiinaRoose,SJonathanChapman,andPhilipKMaini. Mathematicalmodelsofavasculartumorgrowth.SIAMreview,49(2):179–208,2007

    2007

  44. [44]

    Reaction-diffusion model for the growth of avascular tumor.Physical Review E, 65(2):021907, 2002

    SC Ferreira Jr, Marcelo Lobato Martins, and MJ Vilela. Reaction-diffusion model for the growth of avascular tumor.Physical Review E, 65(2):021907, 2002. 35

    2002

  45. [45]

    Avasculargrowth,angiogenesisandvasculargrowthinsolidtumours:Themathematicalmodellingofthestagesoftumour development.Mathematical and computer modelling, 23(6):47–87, 1996

    MarkAJChaplain. Avasculargrowth,angiogenesisandvasculargrowthinsolidtumours:Themathematicalmodellingofthestagesoftumour development.Mathematical and computer modelling, 23(6):47–87, 1996

    1996

  46. [46]

    Physicell: An open source physics-based cell simulator for 3-d multicellular systems.PLoS computational biology, 14(2):e1005991, 2018

    Ahmadreza Ghaffarizadeh, Randy Heiland, Samuel H Friedman, Shannon M Mumenthaler, and Paul Macklin. Physicell: An open source physics-based cell simulator for 3-d multicellular systems.PLoS computational biology, 14(2):e1005991, 2018

    2018

  47. [47]

    Heterogeneous tumour modeling using physicell and its implications in precision medicine

    Miloš Savić, Vladimir Kurbalija, Igor Balaz, and Mirjana Ivanović. Heterogeneous tumour modeling using physicell and its implications in precision medicine. InCancer, Complexity, Computation, pages 157–189. Springer, 2022

    2022

  48. [48]

    LouisShuoWang,JiguangYu,ShijiaLi,andZonghaoLiu.Analysisandmean-fieldlimitofahybridpde-abmmodelingangiogenesis-regulated resistance evolution.Mathematics, 13(17):2898, 2025

    2025

  49. [49]

    Bidirectionalendothelialfeedbackdrivesturing-vascular patterning and drug-resistance niches: a hybrid pde-agent-based study.Bioengineering, 12(10):1097, 2025

    ZonghaoLiu,LouisShuoWang,JiguangYu,JilinZhang,EricaMartel,andShijiaLi. Bidirectionalendothelialfeedbackdrivesturing-vascular patterning and drug-resistance niches: a hybrid pde-agent-based study.Bioengineering, 12(10):1097, 2025

    2025

  50. [50]

    Nature Reviews Immunology, 17(9):559–572, 2017

    NishaNagarsheth,MaxSWicha,andWeipingZou.Chemokinesinthecancermicroenvironmentandtheirrelevanceincancerimmunotherapy. Nature Reviews Immunology, 17(9):559–572, 2017

    2017

  51. [51]

    Role of chemokines in aging and age-related diseases.Mechanisms of Ageing and Development, 223:112009, 2025

    Jitendra Kumar Chaudhary, Ajay Kumar Danga, Anita Kumari, Akshay Bhardwaj, and Pramod C Rath. Role of chemokines in aging and age-related diseases.Mechanisms of Ageing and Development, 223:112009, 2025

    2025

  52. [52]

    Globalwell-posednessandstabilityofnonlocaldamage-structuredlineagemodel with feedback and dedifferentiation.Mathematics, 13(22):3583, 2025

    YeLiang,LouisShuoWang,JiguangYu,andZonghaoLiu. Globalwell-posednessandstabilityofnonlocaldamage-structuredlineagemodel with feedback and dedifferentiation.Mathematics, 13(22):3583, 2025

    2025

  53. [53]

    Harnessing the chemokine system to home car-t cells into solid tumors.Cell Reports Medicine, 3(3), 2022

    Jade Foeng, Iain Comerford, and Shaun R McColl. Harnessing the chemokine system to home car-t cells into solid tumors.Cell Reports Medicine, 3(3), 2022

    2022

  54. [54]

    Type, density, and location of immune cells within human colorectal tumors predict clinical outcome.Science, 313(5795):1960–1964, 2006

    JérômeGalon,AnneCostes,FatimaSanchez-Cabo,AmosKirilovsky,BernhardMlecnik,ChristineLagorce-Pagès,MarieTosolini,Matthieu Camus, Anne Berger, Philippe Wind, et al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome.Science, 313(5795):1960–1964, 2006

    1960

  55. [55]

    Miguel Lopez de Rodas, Maria Villalba-Esparza, Miguel F Sanmamed, Lieping Chen, David L Rimm, and Kurt A Schalper. Biological and clinical significance of tumour-infiltrating lymphocytes in the era of immunotherapy: a multidimensional approach.Nature Reviews Clinical Oncology, 22(3):163–181, 2025

    2025

  56. [56]

    JiguangYu,LouisShuoWang,ZonghaoLiu,andJingfengLiu. Patternsuppressionandrecoveryunderone-wayversustwo-waychemotactic coupling in hybrid partial differential equation–ordinary differential equation models.Transport Phenomena, (0), 2026

    2026

  57. [57]

    Tumour-infiltrating lymphocytes: from prognosis to treatment selection.British journal of cancer, 128(3):451–458, 2023

    Koen Brummel, Anneke L Eerkens, Marco de Bruyn, and Hans W Nijman. Tumour-infiltrating lymphocytes: from prognosis to treatment selection.British journal of cancer, 128(3):451–458, 2023

    2023

  58. [58]

    Approaches to treat immune hot, altered and cold tumours with combination immunotherapies.Nature reviews Drug discovery, 18(3):197–218, 2019

    Jérôme Galon and Daniela Bruni. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies.Nature reviews Drug discovery, 18(3):197–218, 2019

    2019

  59. [59]

    T cell exclusion, immune privilege, and the tumor microenvironment.Science, 348(6230):74–80, 2015

    Johanna A Joyce and Douglas T Fearon. T cell exclusion, immune privilege, and the tumor microenvironment.Science, 348(6230):74–80, 2015

    2015

  60. [60]

    Extracellularmatrixstiffnessandtumor-associatedmacrophagepolarization: new fields affecting immune exclusion.Cancer Immunology, Immunotherapy, 73(6):115, 2024

    Ke-XunYu,Wei-JieYuan,Hui-ZhenWang,andYong-XiangLi. Extracellularmatrixstiffnessandtumor-associatedmacrophagepolarization: new fields affecting immune exclusion.Cancer Immunology, Immunotherapy, 73(6):115, 2024

    2024

  61. [61]

    Algebraic–spectralthresholdsanddiscrete–continuousstabilitytransferinleslie–gowersystems.Electronic Research Archive, 34(1):251–290, 2026

    LouisShuoWangandJiguangYu. Algebraic–spectralthresholdsanddiscrete–continuousstabilitytransferinleslie–gowersystems.Electronic Research Archive, 34(1):251–290, 2026

    2026

  62. [62]

    Cancer immune exclusion: breaking the barricade for a successful immunotherapy.Frontiers in Oncology, 13:1135456, 2023

    Sofia Bruni, Maria Florencia Mercogliano, Florencia Luciana Mauro, Rosalia Inés Cordo Russo, and Roxana Schillaci. Cancer immune exclusion: breaking the barricade for a successful immunotherapy.Frontiers in Oncology, 13:1135456, 2023

    2023

  63. [63]

    Initiationofslimemoldaggregationviewedasaninstability.Journaloftheoreticalbiology,26(3):399–415, 1970

    EvelynFKellerandLeeASegel. Initiationofslimemoldaggregationviewedasaninstability.Journaloftheoreticalbiology,26(3):399–415, 1970

    1970

  64. [64]

    A user’s guide to pde models for chemotaxis.Journal of mathematical biology, 58(1):183–217, 2009

    Thomas Hillen and Kevin J Painter. A user’s guide to pde models for chemotaxis.Journal of mathematical biology, 58(1):183–217, 2009

    2009

  65. [65]

    Analysis framework for stochastic predator–prey model with demographic noise.Results in Applied Mathematics, 27:100621, 2025

    Louis Shuo Wang and Jiguang Yu. Analysis framework for stochastic predator–prey model with demographic noise.Results in Applied Mathematics, 27:100621, 2025

    2025

  66. [66]

    From 1970 until present: The keller-segel model in chemotaxis and its consequences

    D Horstmann. From 1970 until present: The keller-segel model in chemotaxis and its consequences. i.Jahresber. Deutsch. Math.-Verein., 105:103, 2003

    1970

  67. [67]

    Aggregation vs

    Michael Winkler. Aggregation vs. global diffusive behavior in the higher-dimensional keller–segel model.Journal of Differential Equations, 248(12):2889–2905, 2010

    2010

  68. [68]

    Boundednessinaquasilinearparabolic–parabolickeller–segelsystemwithsubcriticalsensitivity.Journal of Differential Equations, 252(1):692–715, 2012

    YoushanTaoandMichaelWinkler. Boundednessinaquasilinearparabolic–parabolickeller–segelsystemwithsubcriticalsensitivity.Journal of Differential Equations, 252(1):692–715, 2012

    2012

  69. [69]

    A damage-structured pde model of stem cell hierarchies: The dual role of dedifferentiation in tissue homeostasis and aging.Plos one, 21(2):e0335163, 2026

    Louis Shuo Wang, Jiguang Yu, and Zonghao Liu. A damage-structured pde model of stem cell hierarchies: The dual role of dedifferentiation in tissue homeostasis and aging.Plos one, 21(2):e0335163, 2026

    2026

  70. [70]

    Keller-segel chemotaxis models: a review.Acta Applicandae Mathematicae, 171(1):6, 2021

    Gurusamy Arumugam and Jagmohan Tyagi. Keller-segel chemotaxis models: a review.Acta Applicandae Mathematicae, 171(1):6, 2021

    2021

  71. [71]

    A finite volume scheme for the patlak–keller–segel chemotaxis model.Numerische Mathematik, 104(4):457–488, 2006

    Francis Filbet. A finite volume scheme for the patlak–keller–segel chemotaxis model.Numerische Mathematik, 104(4):457–488, 2006

    2006

  72. [72]

    Asecond-orderpositivitypreservingcentral-upwindschemeforchemotaxisandhaptotaxismodels

    AlinaChertockandAlexanderKurganov. Asecond-orderpositivitypreservingcentral-upwindschemeforchemotaxisandhaptotaxismodels. Numerische Mathematik, 111(2):169–205, 2008

    2008

  73. [73]

    Conservative upwind finite-element method for a simplified keller–segel system modelling chemotaxis.IMA journal of numerical analysis, 27(2):332–365, 2007

    Norikazu Saito. Conservative upwind finite-element method for a simplified keller–segel system modelling chemotaxis.IMA journal of numerical analysis, 27(2):332–365, 2007

    2007

  74. [74]

    Nonhomogeneous linear and quasilinear elliptic and parabolic boundary value problems

    Herbert Amann. Nonhomogeneous linear and quasilinear elliptic and parabolic boundary value problems. InFunction spaces, differential operators and nonlinear analysis, pages 9–126. Springer, 1993

    1993

  75. [75]

    Springer, 2003

    James Dickson Murray and James Dickson Murray.Mathematical biology: II: spatial models and biomedical applications, volume 18. Springer, 2003. 36

    2003