Tutorials › Quantitative Systems Pharmacology (QSP) › Bounded and Unbounded Population Growth

Cell Population Dynamics · Tutorial 177 of 799

Bounded and Unbounded Population Growth

Learn how net growth and density-dependent feedback determine whether a cell population remains bounded or grows without limit.

Intermediate 10 min read

What You'll Learn

  • Derive the exponential trajectory from a constant per-capita net growth rate.
  • Identify when a linear cell-count model is bounded or unbounded.
  • Derive the logistic trajectory and its conditions for boundedness.
  • Translate birth, death, and crowding parameters into a carrying capacity.
  • Diagnose when apparent plateaus do not establish density regulation.

Why Boundedness Depends on the Growth Law

A cell count can rise rapidly without establishing that growth will continue without limit. A model with a constant net growth rate predicts exponential increase whenever proliferation exceeds loss. A density-regulated model instead reduces per-capita growth as the population grows, which can produce a finite limiting population. The distinction matters when extrapolating cell-count trajectories: the same short-term increase may be consistent with very different long-term behavior.

Here, \(N(t)\) denotes the number of cells at time \(t\), measured in cells. We consider deterministic, continuous-time models with \(N(0)=N_0>0\). A trajectory is bounded if there is a finite constant \(M\) such that \(N(t)\leq M\) for every \(t\geq 0\). It is unbounded if no such constant exists. These definitions concern the model's mathematical prediction, not a claim that a real culture can contain infinitely many cells.

Definition: A constant-net-growth model has the form \(\frac{dN}{dt}=gN\), where \(g\) is a constant net per-capita growth rate (time\(^{-1}\)). A logistic density-regulated model has the form \(\frac{dN}{dt}=rN(1-N/K)\), where \(r\) is the low-density per-capita growth rate (time\(^{-1}\)) and \(K\) is a positive carrying capacity (cells).

Constant Net Growth: A Linear ODE With Exponential Trajectories

In the constant-net-growth model, \(g\) is proliferation minus loss per cell per unit time, assumed constant and independent of population size. The model also assumes a well-mixed population and no external cell input or removal beyond the processes represented in \(g\). Separating variables and integrating gives \(\ln(N(t)/N_0)=gt\), so the trajectory is \(N(t)=N_0e^{gt}\). Although the differential equation is linear in \(N\), its solution is exponential.

$$ \frac{dN}{dt}=gN \quad\Longrightarrow\quad \frac{dN}{N}=g\,dt \quad\Longrightarrow\quad N(t)=N_0e^{gt} $$

The sign of \(g\) determines boundedness. If \(g>0\), then \(e^{gt}\) grows without limit as \(t\) increases, so \(N(t)\) is unbounded. If \(g=0\), the population remains at \(N_0\). If \(g<0\), the population decreases toward zero and is bounded above by \(N_0\). This last case describes declining cell abundance, not extinction at a finite time: the continuous exponential trajectory remains positive for every finite \(t\).

Worked Example: Constant Net Growth

Suppose a hypothetical cell population starts at \(N_0=200{,}000\) cells. Its constant proliferation rate is \(0.32\) day\(^{-1}\), and its constant loss rate is \(0.20\) day\(^{-1}\). Under the model assumptions, the net rate is \(g=0.32-0.20=0.12\) day\(^{-1}\).

After 14 days, the predicted count is \(N(14)=200{,}000e^{(0.12)(14)}=200{,}000e^{1.68}\). Using \(e^{1.68}\approx5.366\), this gives \(N(14)\approx1.073\times10^6\) cells.

Because \(g>0\), this model predicts an unbounded trajectory if continued indefinitely. The calculation is a model-based extrapolation; it does not establish that the constant-rate assumption remains biologically plausible as the population expands.

Density Regulation: A Finite Bound From Negative Feedback

The logistic model makes per-capita growth depend on population size: \(N^{-1}dN/dt=r(1-N/K)\). For \(r>0\) and \(K>0\), per-capita growth is positive when \(0<N<K\), zero at \(N=K\), and negative when \(N>K\). These assumptions represent a stable density-dependent limitation: crowding or resource limitation increasingly offsets growth as the population approaches capacity. The model has no immigration, discrete division events, or explicit delay in that feedback.

Separating variables and integrating the logistic equation yields the trajectory below. The constant is set by \(N(0)=N_0\), giving the stated solution.

$$ N(t)=\frac{K}{1+\left(\frac{K}{N_0}-1\right)e^{-rt}} $$
Theorem: For the logistic model with \(r>0\), \(K>0\), and \(N_0>0\), the trajectory remains positive and bounded for all \(t\geq0\). If \(N_0<K\), it increases toward \(K\); if \(N_0>K\), it decreases toward \(K\); and if \(N_0=K\), it remains constant.

The solution makes the bound precise. If \(N_0<K\), then \(K/N_0-1>0\); the denominator is at least 1, so \(0<N(t)\leq K\). If \(N_0>K\), then \(-1<K/N_0-1<0\). Since \(0<e^{-rt}\leq1\), the denominator lies between \(K/N_0\) and 1, so \(K\leq N(t)\leq N_0\). In both cases the denominator remains positive, ruling out a singularity and establishing a finite upper bound. As \(t\) tends to infinity, \(e^{-rt}\) tends to zero and \(N(t)\) tends to \(K\).

Worked Example: Logistic Growth Below Capacity

Consider a hypothetical population with \(N_0=200{,}000\) cells, \(r=0.40\) day\(^{-1}\), and \(K=1{,}000{,}000\) cells. The assumptions \(r>0\), \(K>0\), and \(N_0>0\) hold, so the logistic boundedness result applies.

At day 10, \(N(10)=\frac{1{,}000{,}000}{1+(1{,}000{,}000/200{,}000-1)e^{-(0.40)(10)}}\). The ratio \(K/N_0-1\) is 4, and \(e^{-4}\approx0.01832\). Therefore, \(N(10)\approx1{,}000{,}000/[1+4(0.01832)]\approx931{,}700\) cells.

The predicted count is below capacity and approaches \(1{,}000{,}000\) cells as time increases. Under this model, the population cannot exceed \(K\) when it starts below \(K\). This is a consequence of the assumed density feedback, not a conclusion that can be drawn from the day-10 count alone.

Reading Boundedness From Birth and Death Parameters

A useful density-regulated form separates low-density proliferation and loss from a crowding penalty: \(\frac{dN}{dt}=(b_0-d_0)N-cN^2\). Here \(b_0\) and \(d_0\) are nonnegative low-density rates in time\(^{-1}\), and \(c>0\) is a density-feedback coefficient in cells\(^{-1}\) time\(^{-1}\). The quadratic term means that the per-capita net rate, \(b_0-d_0-cN\), declines as \(N\) rises.

Write \(r=b_0-d_0\). When \(r>0\), the equation is logistic with \(K=r/c\), so a positive initial population remains bounded. When \(r=0\), the equation becomes \(dN/dt=-cN^2\), whose solution is \(N(t)=N_0/(1+cN_0t)\); it declines and is bounded by \(N_0\). When \(r<0\), both \(rN\) and \(-cN^2\) are negative for \(N>0\), so the population decreases and is likewise bounded. Thus \(c>0\) prevents unbounded positive growth in all three cases, although only \(r>0\) gives a positive equilibrium capacity.

Worked Example: Conditions for a Positive Capacity

Suppose a hypothetical population has low-density proliferation rate \(b_0=0.30\) day\(^{-1}\), loss rate \(d_0=0.10\) day\(^{-1}\), and crowding coefficient \(c=2.0\times10^{-7}\) cells\(^{-1}\) day\(^{-1}\). The net low-density rate is \(r=b_0-d_0=0.20\) day\(^{-1}\), which is positive, and \(c>0\).

The corresponding capacity is \(K=r/c=(0.20\ {\rm day}^{-1})/(2.0\times10^{-7}\ {\rm cells}^{-1}{\rm day}^{-1})=1.0\times10^6\) cells. The units reduce to cells. Because the parameters meet the conditions \(r>0\) and \(c>0\), the model predicts bounded trajectories approaching this capacity for any positive initial count.

If instead \(b_0\leq d_0\) while \(c\) remains positive, there is no positive capacity \(K=r/c\). The population does not grow toward a positive equilibrium; it declines or, at \(r=0\), decreases through the crowding term. Capacity should therefore not be reported as a positive biological quantity when the parameter conditions do not support it.

Diagnostics, Alternatives, and Limits

A finite sequence of increasing counts cannot by itself establish that growth is unbounded: a density-regulated population can increase for a long time while still approaching a finite limit. Conversely, a short apparent plateau does not prove that the limit is a stable carrying capacity. Measurement variation, changing proliferation or death rates, depletion of a resource, and changes in the culture environment can all produce patterns that a constant-parameter model misses.

Inspect whether the assumed mechanism is consistent with the setting. For constant net growth, a plot of log count against time should be approximately linear if the model is appropriate, because \(\ln N(t)=\ln N_0+gt\). For logistic growth, estimates of per-capita growth should decline with density and approach zero near the proposed \(K\). Systematic curvature in the log-count plot challenges constant \(g\); continued positive per-capita growth near the fitted capacity challenges the logistic interpretation. These diagnostics can reveal mismatch but do not, by themselves, identify its cause.

ModelKey assumptionLong-term prediction for positive initial count
Constant net growth, \(dN/dt=gN\)Net per-capita rate is constant and independent of densityUnbounded if \(g>0\); bounded if \(g\leq0\)
Logistic, \(dN/dt=rN(1-N/K)\)Per-capita growth declines linearly with density, with \(r>0\) and \(K>0\)Bounded and approaches \(K\)
Quadratic density regulation, \(dN/dt=rN-cN^2\)A positive crowding coefficient reduces net growth as density risesBounded for \(c>0\); approaches a positive capacity only if \(r>0\)

The logistic model is a compact alternative to constant net growth when density dependence is plausible, but its simplicity is also a limitation. It treats feedback as immediate and smooth, and its fitted \(K\) is conditional on the modeled conditions. A population can exceed an estimated capacity temporarily if assumptions change or if the model is not appropriate; the theorem guarantees boundedness only under the stated logistic model and parameter conditions. More detailed mechanisms may be needed when resources, compartments, or time-varying rates matter.

Common Mistake

A common error is to reason that “the ODE is linear, so cell count grows linearly in time.” The correction is to distinguish the equation's dependence on \(N\) from the shape of its solution: \(dN/dt=gN\) is linear in \(N\), but solving it gives \(N(t)=N_0e^{gt}\). If \(g>0\), that solution is unbounded. Boundedness follows only when the rate or another modeled mechanism prevents indefinite positive growth; it is not implied by calling an equation linear.

Key takeaway: For a positive initial population, constant positive net per-capita growth gives unbounded exponential increase. Logistic density regulation with \(r>0\) and \(K>0\) gives bounded trajectories approaching \(K\). The conclusion depends on the growth law and its assumptions, not simply on whether observed counts rise.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. For \(dN/dt=-0.05N\) day\(^{-1}\) and \(N_0=10{,}000\) cells, is the trajectory bounded? Give its long-term limit.
  2. A logistic model has \(r=0.3\) day\(^{-1}\), \(K=500{,}000\) cells, and \(N_0=700{,}000\) cells. Is the trajectory bounded, and does it increase or decrease initially?
  3. In \(dN/dt=rN-cN^2\), let \(r=0.12\) day\(^{-1}\) and \(c=3.0\times10^{-7}\) cells\(^{-1}\) day\(^{-1}\). Calculate the positive capacity and state whether the trajectory is bounded for positive \(N_0\).
  4. Why does observing a plateau over a short interval not prove that the population has a stable carrying capacity?

Solutions. (1) The solution is \(N(t)=10{,}000e^{-0.05t}\) cells. It is bounded above by \(10{,}000\) cells and tends to zero. (2) Since \(r>0\), \(K>0\), and \(N_0>0\), the trajectory is bounded and approaches \(K\). Because \(N_0>K\), \(dN/dt=rN(1-N/K)<0\) initially, so it decreases. (3) \(K=r/c=(0.12)/(3.0\times10^{-7})=400{,}000\) cells. Since \(r>0\) and \(c>0\), positive initial populations are bounded and approach this capacity. (4) A short plateau could reflect measurement variation, changing rates, or a temporary environmental effect. It does not establish the persistent density-dependent feedback assumed by the logistic model.

Clinical Biostats