---
title: McKendrick–von Foerster Equations
url: https://www.emergentmind.com/topics/mckendrick-von-foerster-equations
type: topic
---

# McKendrick–von Foerster Equations

The McKendrick–von Foerster Equations constitute a canonical mathematical framework for the dynamics of structured populations, most typically age or size-structured. Formulated as first-order partial differential equations (PDEs) with nonlocal boundary conditions, these models rigorously encapsulate vital demographic processes—aging, death, renewal (birth), and in extensions, fragmentation or fission. Their well-posedness, spectral properties, measure-theoretic generalizations, numerical analysis, and applications in population biology, epidemiology, and ecology have been the subject of a large and deep literature.

## 1. Mathematical Formulation and Abstract Cauchy Problem

The classical McKendrick–von Foerster equation describes the density \(u(t,a)\) of individuals at age \(a \ge 0\) and time \(t \ge 0\), evolving according to
\[
\begin{cases}
\partial_t u(t,a) + \partial_a u(t,a) + B(a) u(t,a) = 0, \quad t>0,\, a>0, \\
u(t,0) = \int_0^{\infty} \beta(a)\,u(t,a)\,da, \quad t>0, \\
u(0,a) = u_0(a), \quad a > 0.
\end{cases}
\]
Here, \(B(a)\) is the death (or exiting) rate at age \(a\), and \(\beta(a)\) is the fertility (renewal) kernel. This system is interpreted on the Banach lattice \(X = L^1(\mathbb{R}_+)\), or more generally in a weighted \(L^1(w\,da)\) if \(B(a)\) is unbounded at infinity.

In operator-theoretic terms, the dynamics are triggered by the abstract Cauchy problem:
\[
\frac{d}{dt} u(t) = (L + A) u(t), \quad u(0) = u_0,
\]
with
\begin{align*}
L u &= -\partial_a u - B(a) u, \quad D(L) = \{u \in L^1 : \partial_a u \in L^1,\, u(0) = 0\}, \\
A u(a) &= \delta_{a=0} \int_0^{\infty} \beta(s) u(s) ds.
\end{align*}
The operator \(L\) is hypodissipative and generates a positive semigroup, while \(A\) is a bounded, rank-one (hence compact) operator [1310.7773].

## 2. Spectral Theory, Long-Time Asymptotics, and the Spectral Gap

The spectral properties of the semigroup generated by \(\Lambda = L + A\) are foundational:

- **Main Spectral Theorem:** Under natural positivity, irreducibility, and regularity conditions on \(B\) and \(\beta\), the spectral bound \(s(\Lambda)\) is an algebraically simple real eigenvalue \(\lambda^* > 0\). There exist strictly positive eigenvectors \(\phi \in D(\Lambda)\), \(\psi \in D(\Lambda^*)\), normalized such that \(\int_0^\infty \psi(a)\phi(a) da = 1\), \(\|\phi\|_{L^1} = 1\).

- **Asymptotic Behavior:** The solution decomposes as
\[
u(t, a) = e^{\lambda^* t} \langle \psi, u_0 \rangle \phi(a) + R(t, a), \quad \|R(t, \cdot)\|_{L^1} \leq C' e^{-(\lambda^* - a) t} \|u_0\|_{L^1},
\]
where \(a < \lambda^*\) tracks a uniform spectral gap. Thus, any nontrivial initial datum converges exponentially in \(L^1\) to the unique asymptotic profile \(\phi\) at the population Malthusian parameter \(\lambda^*\) [1310.7773].

The proof strategy involves: (i) Duhamel/ Dyson–Phillips expansion expressing the semigroup as a sum of the semigroup generated by \(L\) and correction terms, (ii) detailed spectral mapping and Weyl-type theorems, (iii) constructive Krein–Rutman approach for positivity and simple leading eigenvalue, and (iv) explicit spectral gap analysis via compactness of the regularizing component \(A\).

## 3. Measure-Valued Extensions and Generalized Entropy Methods

Modern treatments rigorously extend the McKendrick–von Foerster framework to the space of bounded nonnegative Radon measures \(\mathcal{M}^+([0,\infty))\), which encompasses arbitrary initial data—including singular (e.g., sum of Dirac masses) and weakly convergent limits—in place of function spaces. In this setting, solutions admit distributional sense:
\[
\partial_t \mu_t + \partial_a \mu_t + B(a) \mu_t = 0; \quad \mu_t(\{0\}) = \int_0^\infty b(a)\,\mu_t(da)
\]
[1604.07657].

Key results:
- **Existence and Uniqueness:** For any measure initial datum, the solution exists and is unique.
- **Exponential Convergence:** There exist explicit weights \(\varphi\in C_b([0,\infty))\) and constants \(\alpha>0\) such that
\[
\int_0^\infty \varphi(a)\left|\mu_t(da) - m_0 N(a) da\right| \leq e^{-\alpha t}\int_0^\infty \varphi(a)\left|\mu_0(da) - m_0 N(a) da\right|,
\]
where \(N(a)\) is the stationary age density and \(m_0 = \int P(a)\mu_0(da)\) the preserved mass.
- **Technical Approach:** The generalized entropy method, extended to measures using Reshetnyak’s theorem and convex-analytic recession functions, ensures the contraction property even in the absence of strong continuity in the total variation norm [1604.07657].

## 4. Population-Structured Extensions and Related Frameworks

The McKendrick–von Foerster model generalizes to cover several advanced structured population scenarios:

- **Size- or Trait-Structured Models:** Replacing the structure variable by physiological size or generic trait, and including non-trivial growth \(g(x)\). Growth–fragmentation equations with McKendrick–von Foerster boundary conditions fit in this class. The nonlocal renewal term is then
\[
g(0) u(t,0) = \int_0^\infty B(x) u(t,x) dx
\]
[2210.07950].

- **Two-Sex and Multi-Compartment Systems:** Systems of coupled McKendrick–von Foerster equations model interacting populations, e.g., two-sex (Fredrickson–Hoppensteadt) or age-and-sex-structured human populations. Each subpopulation is governed by a coupled PDE with interlinked nonlocal boundary or source terms [1410.2660, 1806.01770].

- **Stochastic Generalizations:** The deterministic PDE emerges as the first-moment closure of a stochastic, age-dependent birth-death process. Fully stochastic kinetic theories yield BBGKY-type hierarchies governing evolution of multiparticle age distributions, and recover the McKendrick–von Foerster equation in the mean-field, population-independent-rate limit. Doi–Peliti field-theoretic and path-integral approaches formalize these connections [1506.02111, 1512.05432, 1512.05431].

- **Mutation and Fragmentation:** Systems with multiple types or genetic classes, coupled at division by a mutation or fragmentation kernel, admit matrix-valued generalizations. Under singular scaling (rapid cycling), the system reduces in the hydrodynamic limit to a network of coupled ODEs for class sizes [1605.07788].

## 5. Numerical Schemes and Computational Analysis

Discretization of the McKendrick–von Foerster equation and its extensions is an active area of research:

- **Finite Difference Methods:** Both explicit and implicit finite difference schemes (including Crank–Nicolson θ-schemes) are implemented for the age and time variables, with consistency and convergence rates \(O(h+\tau)\) for smooth solutions; proof relies on semigroup theory and stability via resolvent estimates [1410.2660].

- **Escalator Boxcar Train (EBT) Method:** This cohort-based method approximates the structured density by a finite sum of moving Dirac delta functions (“cohorts”) with ODE-prescribed mass and location evolution. Convergence is established in the space of nonnegative Radon measures equipped with the flat (bounded–Lipschitz) metric. This approach generalizes to coupled systems (e.g., two-sex models) and nonlinear source/boundary conditions [1806.01770, 1506.00016].

- **Diffusive and Nonlinear Variants:** Finite-difference and hybrid characteristic–finite-difference schemes for equations with age diffusion and nonlinear, nonlocal (Robin-type) boundary conditions achieve stability and convergence under suitable CFL-type conditions. Nonlocality and nonlinearity in boundary conditions are addressed via stability-with-thresholds and fixed-point theory [2201.10440, 2201.08561].

- **Representation in ODE/DAE Systems:** For structured transit- or compartmental models, connections allow translation of PDE models into ODE systems (e.g., via reduction with gamma or uniform maturation time kernels), as well as state-dependent distributed delay DDEs [1811.05930].

## 6. Applications: Demography, Ecology, Epidemiology, and Beyond

The McKendrick–von Foerster equations underpin rigorous quantitative analysis in structured demography (age, size, or type):

- **Classical Demography:** They yield precise connections to discrete-age Leslie matrix models, Euler–Lotka equilibrium relations, and can incorporate frequency-dependent selection via replicator–age–structured frameworks; this enables game-theoretic analyses reflecting sex ratio, strategy competition, and reproductive value dynamics [1303.1432].

- **Epidemic Modeling:** Representing the age of infection as a structure variable provides a universal reduction for complex compartmental epidemic models, including SEIR-type and multi-type risk classes. The McKendrick–von Foerster representation captures arbitrary sojourn distributions and yields explicit renewal equations for inference and parameter estimation [2007.09622].

- **Cell Proliferation and Cycle Kinetics:** Multi-compartmental Markov models of cell division converge to the McKendrick–von Foerster equation in the scaling limit. Explicit connections to spatially-extended traveling waves, reaction-diffusion models, and structured branching processes are established [2108.02249].

- **Growth-Fragmentation and Marine Ecosystems:** In marine food-web models, the McKendrick–von Foerster equation (and its diffusive extension) encapsulates body-mass–structured size spectra. Rigorous spectral analysis reveals instability of power-law steady states in the absence of diffusion, and identifies mechanisms for restoration of stability [1001.3826].

- **Measure-theoretic Population Dynamics:** Entropy-type methods and measure-valued solutions allow for analysis of singular or heterogeneous initial conditions and track convergence toward stable stationary distributions in a broad array of scalar and coupled models [1604.07657].

## 7. Generalizations, Assumptions, and Limitations

- **Functional Requirements:** Key assumptions on the rates (local boundedness, positivity, irreducibility/mixing conditions) determine well-posedness in \(L^1\) or weighted spaces. Relaxing these conditions, e.g. to unbounded death rates or more singular fertility kernels, typically requires additional weighted norms or robust measure-space techniques [1310.7773, 2210.07950].

- **Boundary and Nonlocal Nonlinearity:** Nonlinear and nonlocal birth boundary conditions (integral, history-dependent, environmental feedback) introduce analytical complications requiring refined fixed-point, Banach contraction, or threshold-based discretization analysis [1406.6392, 2201.10440].

- **Extensions:** Spatial extension leads to age-structured reaction-diffusion equations. Generalization to multi-type, branching, or interacting populations leads to kinetic (BBGKY-type) or master equation hierarchies, with the mean-field McKendrick–von Foerster equation recovered only for population-independent rates [1512.05432, 1506.02111].

- **Spectral and Eigenstructure Limitations:** In the presence of population- or environment-dependent rates, compactness and spectral gap properties may fail. Nonlinear models typically lack the explicit spectral representation and require alternative Lyapunov or entropy-based global analysis [1310.7773, 1604.07657].

---

The McKendrick–von Foerster equations thus provide a rigorous, universal backbone for structured transport, renewal, and fragmentation processes in mathematical biology, with extensions spanning both continuum and measure-theoretic settings, and with deep connections to stochastic process hierarchies and modern numerical analysis [1310.7773, 1604.07657, 1806.01770, 2210.07950, 1506.02111].

Source: https://www.emergentmind.com/topics/mckendrick-von-foerster-equations