Size-structured populations with growth fluctuations: Feynman--Kac formula and decoupling
Pith reviewed 2026-05-18 22:08 UTC · model grok-4.3
The pith
When internal variables decouple from size in both ensembles, size dynamics transform via random time change into a growth-homogeneous process evaluable by Feynman-Kac tilting.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When decoupling occurs in both the lineage and population ensembles, the size dynamics can be transformed, via a random time change, into a growth-homogeneous process, and expectations can be evaluated through an exponential tilting procedure that follows from the Feynman-Kac formula. The paper further characterizes weaker, ensemble-specific decoupling that holds in only one ensemble and supplies a more general interpretation of tilted expectations in terms of the mass-weighted phenotype distribution.
What carries the argument
The Feynman-Kac formula applied to branching processes, which supplies the exponential tilting that converts expectations between lineage and population measures once decoupling conditions are met.
If this is right
- Size dynamics reduce to a growth-homogeneous process after the random time change.
- Expectations of size-dependent observables become computable by exponential tilting of the lineage measure.
- Weaker decoupling still relates lineage and population quantities when it holds in only one ensemble.
- Tilted expectations receive a direct interpretation as averages under the mass-weighted phenotype distribution.
Where Pith is reading between the lines
- The same tilting procedure could simplify numerical sampling of large heterogeneous cell populations in synthetic-biology models.
- The random-time-change reduction may extend to other stochastic growth processes that admit a branching-process description.
- Experimental cell-size and gene-expression time series could test whether observed distributions match the predicted homogeneous process after time re-scaling.
Load-bearing premise
The fluctuating internal variable directly sets each cell's instantaneous growth rate inside a branching process whose generator permits the Feynman-Kac representation.
What would settle it
Simulate the full model with an internal variable that sets growth rate; apply the random time change and tilting only under the stated decoupling conditions and check whether the resulting size distribution matches the direct simulation.
Figures
read the original abstract
We study a size-structured population model in which individual cells grow at a rate determined by a fluctuating internal variable (e.g., gene expression levels). Many previous models of phenotypically heterogeneous populations can be viewed as special cases of this model, and it has previously been observed that the internal variable decouples from cell size under certain conditions. In this work, we generalize these results and connect them to the Feynman-Kac formula, which yields relationships between the lineage dynamics and population distribution in branching processes. To this end, we derive conditions for decoupling, both in the lineage and population ensemble. When decoupling occurs in both ensembles, the size dynamics can be transformed, via a random time change, into a growth-homogeneous process, and expectations can be evaluated through an exponential tilting procedure that follows from the Feynman-Kac formula. We further characterize weaker, ensemble-specific forms of decoupling that hold in either the lineage or the population ensemble, but not both. We provide a more general interpretation of tilted expectations in terms of the mass-weighted phenotype distribution
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies size-structured branching-process models in which individual growth rates are set by a fluctuating internal variable (e.g., gene-expression level). It derives explicit decoupling conditions between the internal variable and cell size that hold separately in the lineage ensemble and in the population ensemble. When both ensembles decouple, the size process is shown to be equivalent, after a random time change, to a growth-homogeneous Markov process whose expectations are obtained by an exponential tilt furnished by the Feynman–Kac representation of the branching-process generator. Weaker, ensemble-specific decoupling statements are also characterized, and the tilted expectations are re-interpreted as expectations under the mass-weighted phenotype distribution.
Significance. If the integrability conditions required by the Feynman–Kac representation are verified and the decoupling derivations are free of hidden restrictions, the work supplies a unified, mathematically transparent framework that recovers many existing models of phenotypically heterogeneous populations as special cases and furnishes practical computational tools (random time change plus exponential tilting) for lineage-to-population relations. The explicit link to the Feynman–Kac formula and the mass-weighted interpretation are genuine strengths.
major comments (2)
- [Model statement and Feynman–Kac application] The central claim that decoupling in both ensembles permits a random-time-change reduction to a growth-homogeneous process whose expectations follow from Feynman–Kac tilting (abstract and model statement) rests on the representation being well-defined. The manuscript does not state or verify the necessary integrability condition (finite exponential moments of the integrated growth functional) on the growth-rate function of the internal variable. When fluctuations are unbounded this condition can fail even if the formal decoupling equations hold, undermining the claimed equivalence between lineage and population expectations.
- [Derivation of decoupling conditions] The derivation of the decoupling conditions (both ensembles) appears to rely on the generator of the joint size–internal-variable process admitting a multiplicative functional given by the integrated growth rate. It is not shown whether these conditions remain non-tautological when the growth-rate function is allowed to depend arbitrarily on the internal state; a concrete counter-example or explicit verification that the conditions are independent of the choice of generator would strengthen the result.
minor comments (2)
- [Notation] Notation for the lineage versus population measures should be introduced once and used consistently; the current abstract and model section alternate between several symbols for the same objects.
- [Introduction] The statement that “many previous models can be viewed as special cases” would benefit from an explicit table or subsection listing which earlier works are recovered and under what parameter restrictions.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments raise important points regarding the mathematical foundations of our results, particularly the conditions for the Feynman-Kac representation and the nature of the decoupling conditions. We address each major comment below and indicate the revisions we will make to the manuscript.
read point-by-point responses
-
Referee: The central claim that decoupling in both ensembles permits a random-time-change reduction to a growth-homogeneous process whose expectations follow from Feynman–Kac tilting (abstract and model statement) rests on the representation being well-defined. The manuscript does not state or verify the necessary integrability condition (finite exponential moments of the integrated growth functional) on the growth-rate function of the internal variable. When fluctuations are unbounded this condition can fail even if the formal decoupling equations hold, undermining the claimed equivalence between lineage and population expectations.
Authors: We agree that the integrability condition is essential for the Feynman-Kac formula to be applicable and for the equivalence to hold rigorously. The original manuscript assumes this implicitly through the existence of the expectations but does not state it explicitly. In the revised version, we will add a clear statement of the assumption: that the growth-rate function r(x) satisfies E[exp(∫ r(X_s) ds)] < ∞ for the relevant time intervals, where X is the internal variable process. This is a standard condition in the theory of Feynman-Kac semigroups for branching processes. We will also note that for bounded fluctuations this holds automatically, and discuss potential failures for unbounded cases, thereby strengthening the claim. revision: yes
-
Referee: The derivation of the decoupling conditions (both ensembles) appears to rely on the generator of the joint size–internal-variable process admitting a multiplicative functional given by the integrated growth rate. It is not shown whether these conditions remain non-tautological when the growth-rate function is allowed to depend arbitrarily on the internal state; a concrete counter-example or explicit verification that the conditions are independent of the choice of generator would strengthen the result.
Authors: The decoupling conditions are obtained by requiring that the action of the joint generator on test functions of the form f(size) * g(internal) separates appropriately when multiplied by the exponential of the integrated growth. This leads to specific PDEs or ODEs that the growth rate must satisfy relative to the internal dynamics, which are not automatically true for arbitrary dependence. To demonstrate this, we will include in the revision a simple example with a two-state internal variable (e.g., on/off gene expression) where the growth rate depends on the state. In this case, the decoupling holds only when the transition rates and growth values satisfy a particular algebraic relation, which is independent of the specific form of the size growth operator. This shows the conditions are non-tautological and depend on the interplay between growth and internal state transitions. We believe this addresses the concern without altering the core derivations. revision: partial
Circularity Check
No significant circularity; derivation applies standard Feynman-Kac to independently derived decoupling conditions
full rationale
The paper states a size-structured branching-process model with growth rate set by an internal fluctuating variable, then derives explicit conditions for decoupling in lineage and population ensembles directly from the generator. When both ensembles decouple, a random time change reduces size dynamics to a growth-homogeneous process whose expectations are obtained by the exponential tilting that the Feynman-Kac formula supplies for multiplicative functionals. This is a standard application of an external theorem (Feynman-Kac) to a model whose generator and decoupling equations are written out explicitly; no parameter is fitted to data and then relabeled a prediction, no self-citation supplies a uniqueness theorem that forces the result, and the integrability conditions required for the martingale are left as an implicit modeling assumption rather than being smuggled in by definition. The derivation chain therefore remains self-contained and does not reduce to its inputs by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The population obeys a branching process whose generator admits a Feynman-Kac representation for expectations along lineages.
- domain assumption The internal variable directly modulates the instantaneous growth rate of cell size.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
When decoupling occurs in both ensembles, the size dynamics can be transformed, via a random time change, into a growth-homogeneous process, and expectations can be evaluated through an exponential tilting procedure that follows from the Feynman-Kac formula.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Under this decoupling, population-level expectations are obtained from lineage-level expectations by an exponential tilting given by the Feynman–Kac formula.
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.
Forward citations
Cited by 1 Pith paper
-
Probing the Critical Point (CritPt) of AI Reasoning: a Frontier Physics Research Benchmark
CritPt benchmark shows state-of-the-art LLMs reach only 5.7% average accuracy on full-scale unpublished physics research tasks, rising to about 10% with coding tools.
Reference graph
Works this paper leans on
-
[1]
Cell size regulation in bacteria.Physical review letters, 112(20):208102, 2014
Ariel Amir. Cell size regulation in bacteria.Physical review letters, 112(20):208102, 2014
work page 2014
-
[2]
Pierre Auger, Pierre Magal, and Shigui Ruan.Structured population models in biology and epidemiology, volume 1936. Springer, 2008
work page 1936
-
[3]
Nicolas Baca¨ er.A short history of mathematical population dynamics, volume 618. Springer, 2011
work page 2011
-
[4]
Gerald R Benitz and James A Bucklew. Large deviation rate calculations for nonlinear detectors in gaussian noise.IEEE Transactions on Information Theory, 36(2):358–371, 1990
work page 1990
-
[5]
Bernard Bercu, Fabrice Gamboa, and Alain Rouault. Large deviations for quadratic forms of stationary gaussian processes.Stochastic Processes and their Applications, 71(1):75–90, 1997
work page 1997
-
[6]
Jean Bertoin. On a feynman–kac approach to growth-fragmentation semigroups and their asymptotic behaviors.Journal of Functional Analysis, 277(11):108270, 2019
work page 2019
-
[7]
Feynman-kac formula for stochastic hybrid systems.Physical Review E, 95(1):012138, 2017
Paul C Bressloff. Feynman-kac formula for stochastic hybrid systems.Physical Review E, 95(1):012138, 2017
work page 2017
-
[8]
W lodzimierz Bryc and Amir Dembo. Large deviations for quadratic functionals of gaussian processes.Journal of Theoretical Probability, 10(2):307–332, 1997. 27
work page 1997
-
[9]
Dynamical large deviations of diffusions.arXiv preprint arXiv:2210.09040, 2022
Johan du Buisson. Dynamical large deviations of diffusions.arXiv preprint arXiv:2210.09040, 2022
-
[10]
Clotilde Cadart, Sylvain Monnier, Jacopo Grilli, Pablo J S´ aez, Nishit Srivastava, Rafaele Attia, Emmanuel Terriac, Buzz Baum, Marco Cosentino-Lagomarsino, and Matthieu Piel. Size control in mammalian cells involves modulation of both growth rate and cell cycle duration.Nature communications, 9(1):3275, 2018
work page 2018
-
[11]
Dynamical large deviations of linear diffusions.Physical Review E, 107(5):054111, 2023
Johan du Buisson and Hugo Touchette. Dynamical large deviations of linear diffusions.Physical Review E, 107(5):054111, 2023
work page 2023
-
[12]
Option pricing and esscher transform under regime switching.Annals of Finance, 1(4):423–432, 2005
Robert J Elliott, Leunglung Chan, and Tak Kuen Siu. Option pricing and esscher transform under regime switching.Annals of Finance, 1(4):423–432, 2005
work page 2005
-
[13]
Klaus-Jochen Engel, Rainer Nagel, and Simon Brendle.One-parameter semigroups for linear evolution equations, volume 194. Springer, 2000
work page 2000
-
[14]
Steady distribution of the incremental model for bacteria proliferation
Pierre Gabriel and Hugo Martin. Steady distribution of the incremental model for bacteria proliferation.arXiv preprint arXiv:1803.04950, 2018
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[15]
F Gamboa, A Rouault, and M Zani. A functional large deviations principle for quadratic forms of gaussian stationary processes.Statistics & probability letters, 43(3):299–308, 1999
work page 1999
-
[16]
Reinaldo Garc´ ıa-Garc´ ıa, Arthur Genthon, and David Lacoste. Linking lineage and population observables in biological branching processes.Physical Review E, 99(4):042413, 2019
work page 2019
-
[17]
Crispin W Gardiner et al.Handbook of stochastic methods, volume 3. springer Berlin, 1985
work page 1985
-
[18]
Arthur Genthon. From noisy cell size control to population growth: When variability can be beneficial.Physical Review E, 111(3):034407, 2025
work page 2025
-
[19]
Extremal events dictate population growth rate inference.arXiv preprint arXiv:2501.08404, 2025
Trevor GrandPre, Ethan Levien, and Ariel Amir. Extremal events dictate population growth rate inference.arXiv preprint arXiv:2501.08404, 2025
-
[20]
Ya¨ ır Hein and Farshid Jafarpour. Asymptotic decoupling of population growth rate and cell size distribution.Physical Review Research, 6(4):043006, 2024
work page 2024
-
[21]
Ya¨ ır Hein and Farshid Jafarpour. Competition between transient oscillations and early stochasticity in exponentially growing populations.Physical Review Research, 6(3):033320, 2024
work page 2024
-
[22]
Po-Yi Ho, Jie Lin, and Ariel Amir. Modeling cell size regulation: From single-cell-level statistics to molecular mechanisms and population-level effects.Annual review of biophysics, 47(1):251–271, 2018
work page 2018
-
[23]
Farshid Jafarpour. Cell size regulation induces sustained oscillations in the population growth rate.Physical Review Letters, 122(11):118101, 2019
work page 2019
-
[24]
Bridging the timescales of single-cell and population dynamics.Physical Review X, 8(2):021007, 2018
Farshid Jafarpour, Charles S Wright, Herman Gudjonson, Jedidiah Riebling, Emma Dawson, Klevin Lo, Aretha Fiebig, Sean Crosson, Aaron R Dinner, and Srividya Iyer-Biswas. Bridging the timescales of single-cell and population dynamics.Physical Review X, 8(2):021007, 2018
work page 2018
-
[25]
Bacterial persistence: a model of survival in changing environments.Genetics, 169(4):1807–1814, 2005
Edo Kussell, Roy Kishony, Nathalie Q Balaban, and Stanislas Leibler. Bacterial persistence: a model of survival in changing environments.Genetics, 169(4):1807–1814, 2005
work page 2005
-
[26]
Ethan Levien, Trevor GrandPre, and Ariel Amir. Large deviation principle linking lineage statistics to fitness in microbial populations.Physical review letters, 125(4):048102, 2020
work page 2020
-
[27]
Ethan Levien, Joon Ho Kang, Kuheli Biswas, Scott R Manalis, Ariel Amir, and Teemu P Miettinen. Stochasticity in mammalian cell growth rates drives cell-to-cell variability independently of cell size and divisions.bioRxiv, pages 2025–06, 2025
work page 2025
-
[28]
Ethan Levien, Jane Kondev, and Ariel Amir. The interplay of phenotypic variability and fitness in finite microbial populations.Journal of the Royal Society interface, 17(166):20190827, 2020
work page 2020
-
[29]
Ethan Levien, Jiseon Min, Jane Kondev, and Ariel Amir. Non-genetic variability in microbial populations: survival strategy or nuisance?Reports on progress in physics, 84(11):116601, 2021
work page 2021
-
[30]
On population growth in a randomly varying environment
Richard C Lewontin and Daniel Cohen. On population growth in a randomly varying environment. Proceedings of the National Academy of sciences, 62(4):1056–1060, 1969
work page 1969
-
[31]
Jie Lin and Ariel Amir. The effects of stochasticity at the single-cell level and cell size control on the population growth.Cell systems, 5(4):358–367, 2017. 28
work page 2017
-
[32]
From single-cell variability to population growth.Physical Review E, 101(1):012401, 2020
Jie Lin and Ariel Amir. From single-cell variability to population growth.Physical Review E, 101(1):012401, 2020
work page 2020
-
[33]
Shev MacNamara and Gilbert Strang. Operator splitting. InSplitting methods in communication, imaging, science, and engineering, pages 95–114. Springer, 2017
work page 2017
-
[34]
Aline Marguet. A law of large numbers for branching markov processes by the ergodicity of ancestral lineages.ESAIM: Probability and Statistics, 23:638–661, 2019
work page 2019
-
[35]
Uniform sampling in a structured branching population.Bernoulli, 25(4A):2649– 2695, 2019
Aline Marguet. Uniform sampling in a structured branching population.Bernoulli, 25(4A):2649– 2695, 2019
work page 2019
-
[36]
Pierre Moral.Feynman-Kac formulae: genealogical and interacting particle systems with applications. Springer, 2004
work page 2004
-
[37]
Takashi Nozoe, Edo Kussell, and Yuichi Wakamoto. Inferring fitness landscapes and selection on phenotypic states from single-cell genealogical data.PLoS genetics, 13(3):e1006653, 2017
work page 2017
-
[38]
Springer Science & Business Media, 2013
Bernt Oksendal.Stochastic differential equations: an introduction with applications. Springer Science & Business Media, 2013
work page 2013
-
[39]
Phenomenology of stochastic exponential growth.Physical review e, 95(6):062406, 2017
Dan Pirjol, Farshid Jafarpour, and Srividya Iyer-Biswas. Phenomenology of stochastic exponential growth.Physical review e, 95(6):062406, 2017
work page 2017
-
[40]
EO Powell. Growth rate and generation time of bacteria, with special reference to continuous culture.Microbiology, 15(3):492–511, 1956
work page 1956
-
[41]
Cell-size control.Current Biology, 31(21):R1414–R1420, 2021
Nicholas Rhind. Cell-size control.Current Biology, 31(21):R1414–R1420, 2021
work page 2021
-
[42]
Convergence of large-deviation estimators.Physical Review E, 92(5):052104, 2015
Christian M Rohwer, Florian Angeletti, and Hugo Touchette. Convergence of large-deviation estimators.Physical Review E, 92(5):052104, 2015
work page 2015
-
[43]
On the krein-rutman theorem and beyond.arXiv preprint arXiv:2305.06652, 2023
Claudia Fonte Sanchez, Pierre Gabriel, and St´ ephane Mischler. On the krein-rutman theorem and beyond.arXiv preprint arXiv:2305.06652, 2023
-
[44]
John T Sauls, Dongyang Li, and Suckjoon Jun. Adder and a coarse-grained approach to cell size homeostasis in bacteria.Current opinion in cell biology, 38:38–44, 2016
work page 2016
-
[45]
Francis R Sharpe and Alfred J Lotka. A problem in age-distribution.The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 21(124):435–438, 1911
work page 1911
-
[46]
David Siegmund. Importance sampling in the monte carlo study of sequential tests.The Annals of Statistics, pages 673–684, 1976
work page 1976
-
[47]
Cell-size control and homeostasis in bacteria
Sattar Taheri-Araghi, Serena Bradde, John T Sauls, Norbert S Hill, Petra Anne Levin, Johan Paulsson, Massimo Vergassola, and Suckjoon Jun. Cell-size control and homeostasis in bacteria. Current biology, 25(3):385–391, 2015
work page 2015
-
[48]
Sorin T˘ anase-Nicola and Pieter Rein Ten Wolde. Regulatory control and the costs and benefits of biochemical noise.PLoS computational biology, 4(8):e1000125, 2008
work page 2008
-
[49]
Hugo Touchette. Introduction to dynamical large deviations of markov processes.Physica A: Statistical Mechanics and its Applications, 504:5–19, 2018
work page 2018
-
[50]
Robust growth of escherichia coli.Current biology, 20(12):1099–1103, 2010
Ping Wang, Lydia Robert, James Pelletier, Wei Lien Dang, Francois Taddei, Andrew Wright, and Suckjoon Jun. Robust growth of escherichia coli.Current biology, 20(12):1099–1103, 2010
work page 2010
-
[51]
G. F. Webb. The semigroup associated with nonlinear age dependent population dynamics. International Journal of Computational Mathematics and Applications, 9:487–498, 1983
work page 1983
-
[52]
G. F. Webb.Theory of Nonlinear Age-Dependent Population Dynamics, volume 89 ofMonographs and Textbooks in Pure and Applied Mathematics. Marcel Dekker, Inc., New York and Basel, 1985
work page 1985
-
[53]
Mingtao Xia, Chris D Greenman, and Tom Chou. Pde models of adder mechanisms in cellular proliferation.SIAM journal on applied mathematics, 80(3):1307–1335, 2020
work page 2020
-
[54]
A unified framework for measuring selection on cellular lineages and traits.Elife, 11:e72299, 2022
Shunpei Yamauchi, Takashi Nozoe, Reiko Okura, Edo Kussell, and Yuichi Wakamoto. A unified framework for measuring selection on cellular lineages and traits.Elife, 11:e72299, 2022. 29
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.