pith. sign in

arxiv: 2605.20140 · v1 · pith:FQ2M7M3Rnew · submitted 2026-05-19 · 🧮 math.NA · cs.NA

A Novel Stochastic Particle-Field Algorithm for a Reaction-Diffusion-Advection Cancer Invasion Model

Jingyuan Hu , Zhongjian Wang , Jack Xin , Zhiwen Zhang This is my paper

Pith reviewed 2026-05-20 03:13 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords stochastic particle methodreaction-diffusion-advectioncancer invasion modelparticle-in-cell algorithmpositivity preservationerror analysisthree-dimensional simulationhaptotaxis
0
0 comments X

The pith

A stochastic particle-field method solves 3D cancer invasion models while preserving positivity and bounding mass changes.

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

This paper develops a stochastic particle-field algorithm to solve a reaction-diffusion-advection model of cancer invasion in three dimensions. Cell density is represented by empirical measures of particles with variable mass, and chemical concentrations are constructed on a grid using Particle-in-Cell interactions and spectral diffusion solvers. The authors show that particle mass change rates stay bounded in finite time and that the constructed concentrations preserve positivity unconditionally. Numerical experiments back up the error analysis, allowing the first observation of rapid haptotactic cell spread in 3D. A sympathetic reader would value this as a positivity-preserving and stable discretization for simulating complex biological invasion processes.

Core claim

The authors introduce a novel numerical framework for solving a specific biological reaction-diffusion-advection system of cancer growth in three dimensions using particles of variable mass. They adopt empirical particle measures to represent cell density and dynamically construct the concentration fields of multiple related chemical species throughout the 3D domain. Efficient interaction between the particles and the spatial grid is achieved through a Particle-in-Cell algorithm, while diffusion in space is solved rapidly using a spectral method. For this particular system, the rate of change of particle mass remains bounded over finite time intervals. In addition to the inherent positivity,

What carries the argument

Empirical particle measures combined with Particle-in-Cell coupling to represent cell density and construct concentration fields.

If this is right

  • The algorithm enables efficient simulations of cancer invasion in full three-dimensional domains.
  • Cell density and chemical concentrations remain positive without additional restrictions.
  • Theoretical convergence rates are achieved and confirmed in numerical experiments.
  • Rapid cell spread driven by haptotactic flux appears in 3D, matching two-dimensional observations.

Where Pith is reading between the lines

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

  • The bounded mass change could support stable extensions to longer simulation times for tumor progression studies.
  • Similar particle-grid coupling may improve handling of advection terms in other biological reaction-diffusion systems.
  • Direct comparisons with grid-only discretizations could quantify efficiency gains for high-resolution 3D runs.

Load-bearing premise

The approach assumes that the empirical particle measures combined with the Particle-in-Cell coupling accurately discretize the continuous reaction-diffusion-advection system without introducing uncontrolled instabilities or mass leakage in three dimensions.

What would settle it

A numerical run in which particle mass changes at a rate exceeding the bounded value over a finite time interval, or in which negative concentration values appear on the spatial grid.

Figures

Figures reproduced from arXiv: 2605.20140 by Jack Xin, Jingyuan Hu, Zhiwen Zhang, Zhongjian Wang.

Figure 1
Figure 1. Figure 1: Evolution of the cell density u. domain is L = 6.0: E := ∥u1 − u2∥2 max(∥u1∥2, ∥u2∥2) , ∥u∥2 := L 16   X 0≤i,j≤15 u 2 ij   1 2 . (4.2) The numerical convergence of our method as the resolution increases from H = 2 6 to H = 210 is shown in Figure 2a. The spatial convergence order is validated through the slope of log2 E versus log2 H, with a measured value of −2.109, which is higher than the predicted t… view at source ↗
Figure 2
Figure 2. Figure 2: Numerical convergence of our method, and comparison with a second-order [PITH_FULL_IMAGE:figures/full_fig_p028_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of the cell density u and the concentrations v, m, w on the cross￾section at x3 = L/2. tions H = 24 , ..., 2 7 . As shown in [PITH_FULL_IMAGE:figures/full_fig_p029_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Three-dimensional distribution of the cell density [PITH_FULL_IMAGE:figures/full_fig_p030_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Relative error convergence against grid resolution ( [PITH_FULL_IMAGE:figures/full_fig_p030_5.png] view at source ↗
read the original abstract

In this paper, we present a novel numerical framework for solving a specific biological reaction-diffusion-advection system of cancer growth in three dimensions (3D) using particles of variable mass. We adopt empirical particle measures to represent cell density and dynamically construct the concentration fields of multiple related chemical species throughout the 3D domain. Efficient interaction between the particles and the spatial grid is achieved through a Particle-in-Cell (PIC) algorithm, while diffusion in space is solved rapidly using a spectral method. We demonstrate that for this particular system, the rate of change of particle mass remains bounded over finite time intervals. Furthermore, in addition to the inherent positivity preservation of cell density guaranteed by the empirical particle measures, the concentrations constructed by the algorithm are also unconditionally positivity-preserving on the spatial grid. Moreover, we present a rigorous error analysis for the proposed method, and numerical experiments confirm the theoretical convergence rates. To the best of our knowledge, this is the first numerical work to solve this system in three dimensions, wherein a rapid spread of cells driven by haptotactic flux is observed, similar to the behavior documented in the two-dimensional case.

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 introduces a stochastic particle-field algorithm for a 3D reaction-diffusion-advection cancer invasion model. Cell density is represented by empirical particle measures with variable mass, coupled to grid-based chemical concentration fields via a Particle-in-Cell (PIC) scheme, with diffusion handled by spectral methods. The central claims are that particle mass change remains bounded over finite time intervals, concentrations are unconditionally positivity-preserving, a rigorous error analysis is provided, and numerical experiments confirm the predicted convergence rates. The work is positioned as the first 3D simulation of the system, reproducing haptotactic invasion fronts observed in 2D.

Significance. If the bounded-mass and positivity claims together with the error analysis hold for the fully discrete scheme, the approach would provide a useful addition to the numerical toolkit for high-dimensional biological models. The hybrid particle-grid strategy could scale better than pure grid methods for problems with localized cell populations while preserving key structural properties such as positivity.

major comments (2)
  1. [PIC coupling and mass-evolution analysis] The claim that the rate of change of particle mass remains bounded (stated in the abstract and demonstrated for the system) is load-bearing for the overall stability and error analysis. In the section describing the PIC coupling and the subsequent mass-evolution argument, the proof must explicitly control the cumulative mass defect introduced by the 3D interpolation kernels, whose discrete divergence is not identically zero. Without an additional term bounding this defect over the number of time steps required to reach the invasion front, the continuous boundedness result does not automatically transfer to the discrete algorithm.
  2. [Error analysis] In the rigorous error analysis, the convergence statement should include an explicit estimate for the mass-conservation error arising from the PIC interpolation in three dimensions. The current transfer from the continuous system to the fully discrete scheme appears to omit this term, which is necessary to justify the claimed convergence rates under the stated assumptions.
minor comments (2)
  1. [Abstract] The abstract refers to 'multiple related chemical species' without indicating their number or the precise reaction terms; a brief enumeration would improve readability.
  2. [Numerical experiments] In the numerical experiments section, the captions and axis labels of the 3D invasion-front visualizations could be expanded to include the specific grid size, time step, and kernel width used, facilitating direct comparison with the theoretical rates.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and indicate the changes we will make to strengthen the presentation of the mass-evolution analysis and error estimates.

read point-by-point responses

  1. Referee: [PIC coupling and mass-evolution analysis] The claim that the rate of change of particle mass remains bounded (stated in the abstract and demonstrated for the system) is load-bearing for the overall stability and error analysis. In the section describing the PIC coupling and the subsequent mass-evolution argument, the proof must explicitly control the cumulative mass defect introduced by the 3D interpolation kernels, whose discrete divergence is not identically zero. Without an additional term bounding this defect over the number of time steps required to reach the invasion front, the continuous boundedness result does not automatically transfer to the discrete algorithm.

    Authors: We agree that the non-vanishing discrete divergence of the 3D interpolation kernels introduces a mass defect that must be controlled explicitly for the fully discrete scheme. The original mass-evolution argument was derived at the continuous level; we will add a new lemma in the revised manuscript that bounds the cumulative interpolation defect over a finite number of time steps. The bound will exploit the Lipschitz regularity of the kernels, the boundedness of the advection velocity, and the fact that the invasion front is reached in finite time. This additional estimate will be inserted directly into the stability proof so that the discrete particle-mass bound follows from the continuous one with an explicit remainder term. revision: yes

  2. Referee: [Error analysis] In the rigorous error analysis, the convergence statement should include an explicit estimate for the mass-conservation error arising from the PIC interpolation in three dimensions. The current transfer from the continuous system to the fully discrete scheme appears to omit this term, which is necessary to justify the claimed convergence rates under the stated assumptions.

    Authors: We thank the referee for this observation. The error analysis section currently transfers continuous estimates to the discrete setting but does not isolate the mass-conservation error generated by the 3D PIC interpolation. In the revision we will augment the error theorem with an additional term that quantifies this interpolation error. Under the mesh-size and time-step restrictions already present in the manuscript, we will show that the extra term is of order O(Δx + Δt) and can be absorbed into the existing convergence rate without altering the final statement. The revised proof will therefore contain a complete accounting of all discretization contributions. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation uses standard PIC/spectral techniques on external model

full rationale

The paper introduces an empirical-particle PIC algorithm with spectral diffusion for a given reaction-diffusion-advection cancer model. Boundedness of particle-mass rate, unconditional positivity on the grid, and error estimates are asserted from the continuous-system properties together with standard discretization arguments. No quoted step reduces a claimed prediction or bound to a fitted parameter, self-citation chain, or definitional renaming; the central results remain independent of the algorithm's own outputs. This is the normal non-circular outcome for a methods paper that applies established tools to a new setting.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on standard numerical PDE assumptions and the specific form of the biological model; no new free parameters, invented entities, or ad-hoc axioms beyond domain-standard well-posedness are introduced in the abstract.

axioms (1)
  • domain assumption The underlying reaction-diffusion-advection PDE system is well-posed and admits solutions with sufficient regularity for the error analysis to hold.
    Required to justify bounded mass change, positivity preservation, and convergence rates.

pith-pipeline@v0.9.0 · 5735 in / 1306 out tokens · 42021 ms · 2026-05-20T03:13:05.273482+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

  1. [1]

    Chertock and A

    A. Chertock and A. Kurganov, A second-order positivity preserving central-upwind scheme for chemotaxis and haptotaxis models,Numerische Mathematik.111(2008) 169–205

    work page 2008

  2. [2]

    S. Wise, J. Lowengrub, H. Frieboes and V. Cristini, Three-dimensional multispecies nonlinear tumor growth I,Journal of Theoretical Biology.(2008)

    work page 2008

  3. [3]

    Anderson and M

    A. Anderson and M. Chaplain, Continuous and Discrete Mathematical Models of Tumor-induced Angiogenesis,Bulletin of Mathematical Biology.60(1998) 857-899

    work page 1998

  4. [4]

    M. H. Zangooei and J. Habibi, Hybrid multiscale modeling and prediction of cancer cell behavior,PLOS ONE.12(2017)

    work page 2017

  5. [5]

    Anderson, A hybrid mathematical model of solid tumour invasion: the importance of cell adhesion,Mathematical medicine and biology: a journal of the IMA.22(2005) 163–186

    A. Anderson, A hybrid mathematical model of solid tumour invasion: the importance of cell adhesion,Mathematical medicine and biology: a journal of the IMA.22(2005) 163–186

    work page 2005

  6. [6]

    Ayati, G

    B. Ayati, G. Webb and A. Anderson, Computational methods and results for struc- tured multiscale models of tumor invasion,Multiscale Modeling and Simulation.5 (2006) 1–20

    work page 2006

  7. [7]

    Z. Wang, J. Xin and Z. Zhang, A Novel Stochastic Interacting Particle-Field Algo- rithm for 3D Parabolic-Parabolic Keller-Segel Chemotaxis System,Journal of Scien- tific Computing.102(2025)

    work page 2025

  8. [8]

    J. Hu, Z. Wang, J. Xin and Z. Zhang, A fast stochastic interacting particle-field method for 3D parabolic parabolic Chemotaxis systems: numerical algorithms and error analysis,arXiv preprint arXiv:2512.03452.(2025)

    work page arXiv 2025

  9. [9]

    B. Hu, Z. Wang, J. Xin and Z. Zhang, A Stochastic Interacting Particle-Field Algo- rithm for a Haptotaxis Advection-Diffusion System Modeling Cancer Cell Invasion, arXiv preprint arXiv:2407.05626.(2024)

    work page arXiv 2024

  10. [10]

    Keller and L

    E. Keller and L. Segel, Initiation of slime mold aggregation viewed as an instability, Journal of Theoretical Biology.26(1970)

    work page 1970

  11. [11]

    Chertock, Y

    A. Chertock, Y. Epshteyn, H. Hu and A. Kurganov, High-order positivity-preserving hybrid finite-volume-finite-difference methods for chemotaxis systems,Advances in Computational Mathematics.44(2018) 327–350. May 20, 2026 1:37 WSPC/INSTRUCTION FILE V5 32Jingyuan Hu, Zhongjian Wang, Jack Xin and Zhiwen Zhang

    work page 2018

  12. [12]

    Shen and J

    J. Shen and J. Xu, Unconditionally bound preserving and energy dissipative schemes for a class of Keller-Segel equations,SIAM Journal on Numerical Analysis.58(2020) 1674–1695

    work page 2020

  13. [13]

    W. Chen, Q. Liu and J. Shen, Error estimates and blow-up analysis of a finite- element approximation for the parabolic-elliptic Keller-Segel system,arXiv preprint arXiv:2212.07655.(2022)

    work page arXiv 2022

  14. [14]

    J. Liu, L. Wang and Z. Zhou, Positivity-preserving and asymptotic preserving method for 2D Keller-Segel equations,Mathematics of computation.87(2018) 1165-1189

    work page 2018

  15. [15]

    Hillen and K

    T. Hillen and K. Painter, Global existence for a parabolic chemotaxis model with prevention of overcrowding,Advances in Applied Mathematics.26(2001) 280–301

    work page 2001

  16. [16]

    Gnanasekaran, A

    S. Gnanasekaran, A. Columbu, R. Fuentes and N. Nithyadevi, Global existence and lower bounds in a class of tumor-immune cell interactions chemotaxis systems,Discrete and Continuous Dynamical Systems.18(2025)

    work page 2025

  17. [17]

    Vel´ azquez, Point dynamics in a singular limit of the Keller–Segel model 1: Motion of the concentration regions,SIAM Journal on Applied Mathematics.(2004) 1198–1223

    J. Vel´ azquez, Point dynamics in a singular limit of the Keller–Segel model 1: Motion of the concentration regions,SIAM Journal on Applied Mathematics.(2004) 1198–1223

    work page 2004

  18. [18]

    Haˇ skovec and C

    J. Haˇ skovec and C. Schmeiser, Stochastic particle approximation for measure valued solutions of the 2D Keller-Segel system,Journal of Statistical Physics.135(2009) 133–151

    work page 2009

  19. [19]

    Liu and R

    J. Liu and R. Yang, A random particle blob method for the Keller-Segel equation and convergence analysis,Mathematics of Computation.86(2017) 725–745

    work page 2017

  20. [20]

    Haˇ skovec and C

    J. Haˇ skovec and C. Schmeiser, Convergence of a stochastic particle approximation for measure solutions of the 2D Keller-Segel system,Communications in Partial Differ- ential Equations.36(2011) 940–960

    work page 2011

  21. [21]

    Mischler and C

    S. Mischler and C. Mouhot, Kac’s program in kinetic theory,Inventiones mathemati- cae.193(2013) 1–147

    work page 2013

  22. [22]

    Liu and R

    J. Liu and R. Yang, Propagation of chaos for the Keller-Segel equation with a loga- rithmic cut-off,Methods and Applications of Analysis.26(2019) 319–348

    work page 2019

  23. [23]

    S. Khan, J. Johnson, E. Cartee and Y. Yao, Global regularity of chemotaxis equations with advection,Involve, a Journal of Mathematics.9(2015) 119–131

    work page 2015

  24. [24]

    Jin and L

    S. Jin and L. Li, Random batch methods for classical and quantum interacting particle systems and statistical samplings,Active Particles, Volume 3: Advances in Theory, Models, and Applications.(Springer, 2021) 153–200

    work page 2021

  25. [25]

    L. Chen, S. Wang and R. Yang, Mean-field limit of a particle approximation for the parabolic-parabolic Keller-Segel model,arXiv preprint arXiv:2209.01722.(2022)

    work page arXiv 2022

  26. [26]

    Fournier and M

    N. Fournier and M. Tomaˇ sevi´ c, Particle approximation of the doubly parabolic Keller- Segel equation in the plane,Journal of Functional Analysis.285(2023)

    work page 2023

  27. [27]

    Stevens, The derivation of chemotaxis equations as limit dynamics of moderately interacting stochastic many-particle systems,SIAM Journal on Applied Mathematics

    A. Stevens, The derivation of chemotaxis equations as limit dynamics of moderately interacting stochastic many-particle systems,SIAM Journal on Applied Mathematics. 61(2000) 183–212

    work page 2000

  28. [28]

    B. Hu, Z. Wang, J. Xin and Z. Zhang, A Stochastic Genetic Interact- ing Particle Method for Reaction-Diffusion-Advection Equations,arXiv preprint arXiv:2511.12275.(2025)

    work page arXiv 2025

  29. [29]

    Walker and G

    C. Walker and G. F. Webb, Global existence of classical solutions for a haptotaxis model,SIAM Journal on Mathematical Analysis.38(2007) 1694–1713

    work page 2007

  30. [30]

    Grigorev, V

    I. Grigorev, V. Vshivkov and M. Fedoruk, Numerical ”particle-in-cell” methods: the- ory and applications, Reprint 2012. (VSP, 2002)

    work page 2012

  31. [31]

    G. N. Milstein and M. V. Tretyakov, Stochastic numerics for mathematical physics. ( Springer, 2004)

    work page 2004