A kinetic theory for age-structured stochastic birth-death processes

Chris D. Greenman,Tom Chou

DOI: https://doi.org/10.1103/PhysRevE.93.012112

2015-06-06

Abstract:Classical age-structured mass-action models such as the McKendrick-von Foerster equation have been extensively studied but they are structurally unable to describe stochastic fluctuations or population-size-dependent birth and death rates. Stochastic theories that treat semi-Markov age-dependent processes using e.g., the Bellman-Harris equation, do not resolve a population's age-structure and are unable to quantify population-size dependencies. Conversely, current theories that include size-dependent population dynamics (e.g., mathematical models that include carrying capacity such as the Logistic equation) cannot be easily extended to take into account age-dependent birth and death rates. In this paper, we present a systematic derivation of a new fully stochastic kinetic theory for interacting age-structured populations. By defining multiparticle probability density functions, we derive a hierarchy of kinetic equations for the stochastic evolution of an ageing population undergoing birth and death. We show that the fully stochastic age-dependent birth-death process precludes factorization of the corresponding probability densities, which then must be solved by using a BBGKY-like hierarchy. However, explicit solutions are derived in two simple limits and compared with their corresponding mean-field results. Our results generalize both deterministic models and existing master equation approaches by providing an intuitive and efficient way to simultaneously model age- and population-dependent stochastic dynamics applicable to the study of demography, stem cell dynamics, and disease evolution.

Statistical Mechanics,Populations and Evolution,Quantitative Methods

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a new fully stochastic kinetic theory to describe the interactions of age - structured populations in the processes of birth and death. Specifically, although the existing classical age - structured mass - action models (such as the McKendrick - von Foerster equation) have been widely studied, they are structurally unable to describe stochastic fluctuations or birth and death rates related to population size. Meanwhile, the stochastic theories dealing with semi - Markov age - dependent processes (for example, using the Bellman - Harris equation) cannot analyze the age structure of the population, nor can they quantify the dynamic changes related to population size. On the contrary, the current theories of population dynamics that include size - dependence (for example, mathematical models that include carrying capacity, such as the Logistic equation) are difficult to extend to consider age - dependent birth and death rates. Therefore, the goal of this paper is to derive a hierarchical kinetic equation for the stochastic evolution of an aging population by defining a multi - particle probability density function, thereby providing an intuitive and effective method to simultaneously simulate age - and population - dependent stochastic dynamics, which is applicable to the research of demography, stem cell dynamics and disease evolution.

The key equations mentioned in the paper are as follows:

1. **Cumulative Probability Distribution**:
\[
Q_n(a_n;t)=\int_0^{a_1}dx_1\int_{x_1}^{a_2}dx_2\cdots\int_{x_{n - 1}}^{a_n}dx_nf_n(x_n;t)
\]

2. **Probability Current**:
\[
J_\pm(a_n;t)=J_\pm^\beta(a_n;t)+J_\pm^\mu(a_n;t)
\]

3. **Main Kinetic Equation**:
\[
\frac{\partial f_n(a_n;t)}{\partial t}+\sum_{j = 1}^n\frac{\partial f_n(a_n;t)}{\partial a_j}=-f_n(a_n;t)\sum_{i = 1}^n\gamma_n(a_i)+\sum_{i = 0}^n\int_{a_i}^{a_{i+1}}\mu_{n+1}(y)f_{n+1}(a_i,y,a_{i+1},\ldots,a_n;t)dy
\]

4. **Boundary Condition**:
\[
f_n(a_1 = 0,a_2,\ldots,a_n;t)=f_{n-1}(a_2,\ldots,a_n;t)\sum_{i = 2}^n\beta_{n-1}(a_i)
\]

These equations form a system similar to the BBGKY hierarchy, allowing for a fully stochastic treatment of age - dependent birth and death processes. The paper also shows explicit solutions in two simple limit cases and compares them with the corresponding mean - field results. Finally, the results of the paper generalize deterministic models and existing master equation methods, providing new tools for the study of demography, stem cell dynamics and disease evolution.