pith. sign in

arxiv: 2101.03470 · v1 · submitted 2021-01-10 · 🧬 q-bio.PE · cond-mat.stat-mech

Kinetic theory for structured populations: application to stochastic sizer-timer models of cell proliferation

Mingtao Xia , Tom Chou This is my paper

Pith reviewed 2026-05-24 14:22 UTC · model grok-4.3

classification 🧬 q-bio.PE cond-mat.stat-mech
keywords kinetic theorystructured populationssizer-timer modelcell proliferationgrowth rate stochasticitydemographic stochasticitybirth-death processespopulation dynamics
0
0 comments X

The pith

Kinetic equations unify PDE and birth-death descriptions of age-size structured cell populations by encoding joint growth and demographic stochasticities.

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

The paper derives the full kinetic equations governing the probability density of cells distributed by age and size in stochastic sizer-timer models. These equations incorporate stochasticity both in individual cell growth rates and in the demographic events of birth and death. Averaging the densities produces a second-order PDE that includes growth-rate fluctuations, while taking marginals over the same densities produces a modified birth-death master equation in which age and size modulate the rates of demographic noise. The resulting kinetic framework therefore contains both the deterministic PDE limit and the stochastic master-equation limit as exact special cases obtained by averaging or marginalizing.

Core claim

The central claim is that a single kinetic density whose evolution is governed by the derived integro-differential equations subsumes both the deterministic PDE representation (recovered by averaging) and the birth-death master-equation representation (recovered by marginalization) for structured populations, without requiring additional closure approximations.

What carries the argument

The kinetic density distribution for cell age and size, whose evolution equations encode the joint action of growth-rate stochasticity and demographic stochasticity.

If this is right

  • Averages over the kinetic density recover a second-order PDE that incorporates growth-rate stochasticity.
  • Marginalization over the kinetic density yields a birth-death process whose rates explicitly depend on age and size.
  • The kinetic model contains the deterministic PDE and the stochastic master equation as exact limiting cases.
  • Population statistics for cell proliferation can be computed from a single density that carries both types of noise.

Where Pith is reading between the lines

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

  • Higher-order correlation functions between age, size, and division events could be extracted directly from the kinetic density without separate closures.
  • The same construction might be applied to other structured populations whose individuals carry both continuous traits and discrete events.
  • Numerical solution of the kinetic equations would provide a benchmark for testing whether existing PDE or master-equation approximations remain accurate when both noise sources are simultaneously present.

Load-bearing premise

The joint stochasticities in growth rate and demographics can be encoded in a single kinetic density whose averages and marginals recover the standard PDE and master equation exactly, without hidden approximations or additional closure assumptions.

What would settle it

Individual-based Monte Carlo trajectories of the stochastic sizer-timer process should produce exact numerical agreement with both the averaged PDE moments and the marginal birth-death rates computed from the kinetic density; any systematic discrepancy would falsify the unification.

Figures

Figures reproduced from arXiv: 2101.03470 by Mingtao Xia, Tom Chou.

Figure 1
Figure 1. Figure 1: FIG. 1: A map of boundary condition interdependences for single-density kinetic theory. In (a) we indicate the dependence [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
read the original abstract

We derive the full kinetic equations describing the evolution of the probability density distribution for a structured population such as cells distributed according to their ages and sizes. The kinetic equations for such a "sizer-timer" model incorporates both demographic and individual cell growth rate stochasticities. Averages taken over the densities obeying the kinetic equations can be used to generate a second order PDE that incorporates the growth rate stochasticity. On the other hand, marginalizing over the densities yields a modified birth-death process that shows how age and size influence demographic stochasticity. Our kinetic framework is thus a more complete model that subsumes both the deterministic PDE and birth-death master equation representations for structured populations.

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

0 major / 2 minor

Summary. The manuscript derives kinetic equations governing the joint probability density for structured cell populations in stochastic sizer-timer models, incorporating both demographic stochasticity and stochastic growth rates. Averages over the kinetic density are shown to recover a second-order PDE, while marginalization recovers a modified birth-death master equation; the kinetic description is presented as a unifying first-principles framework that subsumes the deterministic PDE and master-equation representations.

Significance. If the derivations hold without hidden closures, the work supplies a single consistent object (the kinetic density) from which both the deterministic continuum limit and the stochastic demographic process emerge exactly by averaging and marginalization. This supplies a parameter-free bridge between modeling scales and makes falsifiable predictions about how growth-rate fluctuations propagate into population-level statistics.

minor comments (2)
  1. [Abstract] Abstract, final paragraph: the phrasing 'more complete model' is slightly stronger than the technical result (exact subsumption by construction); a neutral phrasing such as 'unifying framework' would align better with the derivation.
  2. [Introduction] The manuscript would benefit from an explicit statement, early in the introduction, of the precise state variables carried by the kinetic density (age, size, growth-rate) and the form of the stochastic increments.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript and for recommending acceptance. The referee's summary accurately captures the central contribution: a kinetic density that yields both the deterministic PDE (via averaging) and the modified birth-death process (via marginalization) without additional closures.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper starts from an explicit joint probability density over cell age, size, and growth-rate state and writes the kinetic evolution equations for that density from individual-level stochastic rules. Averaging the density then produces the deterministic PDE by direct integration, while marginalization produces the birth-death master equation; both recoveries are identities that hold by the definitions of averaging and marginalization. No parameters are fitted to data, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled in. The subsumption claim is therefore the intended mathematical consequence of the joint-density construction rather than a circular reduction of the result to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only abstract available; ledger entries are inferred from the described approach. No free parameters or invented entities are mentioned. The derivation implicitly rests on standard conservation of probability for structured densities.

axioms (1)
  • standard math Probability densities for cell age and size obey a continuity equation that can be closed under the stated stochastic growth and division rules.
    Required for any kinetic derivation of the form described in the abstract.

pith-pipeline@v0.9.0 · 5639 in / 1087 out tokens · 25115 ms · 2026-05-24T14:22:46.372268+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

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

  1. Solving linear-rate ODE hierarchies (like master equations) using closures and operator splitting

    math.NA 2026-05 unverdicted novelty 7.0

    For linear-rate master equations the generating function admits an exact composition-multiplier representation whose Taylor coefficients on any finite window are obtained from a closed lower-triangular ODE of size 2(N...

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages · cited by 1 Pith paper

  1. [1]

    von Foerster H 1959 The Kinetics of Cellular Proliferation, Grune and Stratton 382–407

    work page 1959

  2. [2]

    Taheri-Araghi S, Bradde S, Sauls J T, Hill N S, Levin P A, Paulsson J, Vergassola M and Jun S 2015 Current Biology 25 385–391

    work page 2015

  3. [3]

    Burov S and Kessler D 2018 Bulletin of the American Physical Society 63

    work page 2018

  4. [4]

    Robert L, Hoffmann M, Krell N, Aymerich S, Robert J and Doumic M 2014 BMC Biology 12 17

    work page 2014

  5. [5]

    Perthame B 2008 Introduction to Structured Equations in Biology

    work page 2008

  6. [6]

    Metz J A J and Diekmann O 1986 The Dynamics of Physiologically Structured Populations (Springer)

    work page 1986

  7. [7]

    Sompayrac L and Maaloe O 1973 Nature: New Biology 241 133–135

    work page 1973

  8. [8]

    Huisman O and D’Ari R 1981 Nature 290 797–799

    work page 1981

  9. [9]

    Chandler-Brown D, Schmoller K M, Winetraub Y and Skotheim J M 2017 Current Biology 27 2774–2783

    work page 2017

  10. [10]

    Delarue M, Weissman D and Hallatschek O 2017 PLoS ONE 12 e0182633 13

    work page 2017

  11. [11]

    Wessels J G H 1994 Annual Review of Phytopathology 32 413–437

    work page 1994

  12. [12]

    Modi S, Vargas-Garcia C A, Ghusinga K R and Singh A 2017 Biophysical Journal 112 2408–2418

    work page 2017

  13. [13]

    Xia M, Greenman C D and Chou T 2020 SIAM Journal on Applied Mathematics 80 1307–1335

    work page 2020

  14. [14]

    Bernard E, Doumic M and Gabriel P 2016 Kinetic and Related Models 12 551–571

    work page 2016

  15. [15]

    Greenman C D and Chou T 2016 Physical Review E 93 012112

    work page 2016

  16. [16]

    Chou T and Greenman C D 2016 Journal of Statistical Physics 164 49–76

    work page 2016

  17. [17]

    Greenman C D 2017 Journal of Statistical Mechanics 2017 033101

    work page 2017

  18. [18]

    Vargas-Garcia C A, Soltani M and Singh A 2016 IEEE Life Sciences Letters 2 47–50

    work page 2016

  19. [19]

    Ho P Y, Lin J and Amir A 2018 Annual Review of Biophysics 47 251–271

    work page 2018

  20. [20]

    Kessler D A and Burov S 2017 Physical Review E 96(4) 042139

    work page 2017

  21. [21]

    Nieto C, Vargas-Garcia C and Pedraza J M 2020 bioRxiv:2020.09.29.319251

    work page 2020

  22. [22]

    Durrett R 2005 Cambridge U Press 39 320–353

    work page 2005

  23. [23]

    Popescu D M and Sun S X 2018 Journal of The Royal Society Interface 15 20180086

    work page 2018

  24. [24]

    Auger P, Magal P and Ruan S 2008 Structured Population Models in Biology and Epidemiology vol 1936 (Springer)

    work page 2008

  25. [25]

    (45) Jm,n;m′,n′(t)dt is the probability flux within time [ t,t + dt] from state (m′,n′) to state (m,n ) arising from from cell division

    Xia M, Shao S and Chou T 2020 arXiv:2009.13170 Appendix: conservation of probability We now define probability fluxes Jm,n;m+1,n−1(t) = (m + 1) ∫ dXmdY2n−2dAmdBn−1 ∫ Λ3 dy1dy2ds ˜βm+1,n−1(y1 +y2,y 1,s,t )× ρm+1,n−1(Xm+1[Xm+1 =y1 +y2], Y2n−2, Am+1[Am+1 =s], Bn−1,t ), Jm,n;m−1,n(t) = 2n m ∫ dXmdY2n−2dAmdBn−1 ∫ Λ2 dy1dy2 m∑ i=1 ˜βm−1,n(y1 +y2,y 1,Ai,t )× ρm−1,...

    work page arXiv 2020