🧬 biology

On competition through growth reduction

This paper analyzes a hierarchical population model formulated as a scalar renewal equation, demonstrating that under general assumptions, the system admits a unique non-trivial stationary birth rate and establishing its stability conditions while relating the results to an equivalent partial differential equation formulation.

Original authors: Carles Barril, Àngel Calsina, Odo Diekmann, József Z. Farkas

Published 2026-07-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Carles Barril, Àngel Calsina, Odo Diekmann, József Z. Farkas

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Imagine a forest where every tree is trying to grow taller. In this forest, the only thing that matters for how fast a tree grows is how much sunlight it gets. But here's the catch: the taller trees cast shadows on the shorter ones. The more tall trees there are, the less light the smaller trees get, and the slower they grow.

This paper is a mathematical story about exactly that kind of forest. The authors want to understand: Can this forest settle into a stable, healthy state where trees keep being born, growing, and dying at a steady pace? And if so, is that state safe from chaos?

Here is the breakdown of their journey, translated into everyday language:

1. Two Ways to Look at the Forest

The authors start by looking at the forest in two different ways, like looking at a crowd from a drone versus looking at it from the ground.

The "Drone" View (The PDE): This is the classic way scientists model populations. They look at the whole forest at once, tracking how many trees exist at every specific height at every moment. It's like a complex, moving map.
The "Ground" View (The Delay Equation): The authors realized there's a simpler way to look at it. Instead of tracking every single tree, they decided to just track the birth rate (how many new saplings are born per minute).
- The Analogy: Imagine you are counting the number of babies born in a city. You don't need to know the height of every adult to know how many babies are coming; you just need to know how many adults are there and how fertile they are.
- The "Delay": There is a time lag. A tree born today doesn't start having babies until it is old enough. So, the number of babies born now depends on the population that existed in the past. This is why they call it a "delay equation."

2. The Big Question: Will the Forest Survive?

The authors asked: "If we start with a forest, will it eventually die out, or will it find a happy, steady rhythm?"

They found that the answer depends on a single number they call $R_0$ (the Basic Reproduction Number).

Think of $R_0$ as the "Baby Boom Potential": It asks, "If a new tree is born into a brand new, empty forest (where it gets 100% of the sunlight), how many babies will it have in its lifetime?"
If $R_0 < 1$ : The tree can't even replace itself. The forest will slowly fade away until there are no trees left. The "zero population" state is stable.
If $R_0 > 1$ : The tree can replace itself and then some. The forest will grow until the competition for light gets so fierce that growth slows down, eventually settling into a steady state.

3. The "One True Rhythm"

One of the most interesting findings is about uniqueness.
Usually, in complex systems, you might have many different stable states (e.g., a forest could be stable at 100 trees, or 500 trees, or 1,000 trees, depending on how you start).

The authors proved that for this specific type of forest (where big trees block light for small ones), there is only one possible stable size for the population.

If the forest survives, it will always settle into the exact same birth rate and tree density. It doesn't matter if you start with a few saplings or a whole grove; the forest will eventually find that one specific "Goldilocks" rhythm where the number of births equals the number of deaths.

4. Is the Forest Safe? (Stability)

Once the forest finds this steady rhythm, is it safe? What if a storm knocks down a few trees, or a few extra seeds blow in? Will the forest bounce back, or will it spiral into chaos?

The authors used a mathematical tool called Linearised Stability (which is like testing how a tightrope walker reacts to a small gust of wind).

They showed that under very general and realistic conditions (trees grow slower when crowded, and bigger trees have more babies), this steady rhythm is locally asymptotically stable.
The Metaphor: Imagine a ball sitting at the bottom of a bowl. If you nudge the ball (a small change in the population), it wobbles a bit but eventually rolls back to the bottom. The forest will self-correct and return to its steady state after small disturbances.

5. The "Special Case" Calculation

To prove this, they looked at a specific, simplified scenario where:

Trees only start having babies once they reach a certain "adult" size.
The growth rate slows down in a specific mathematical way as the forest gets crowded.

In this specific case, they were able to write down an exact formula for the steady state using a special mathematical function called the Lambert W function (a bit like a "super-logarithm" used to solve tricky equations). This confirmed that their theory holds up even when you do the hard math.

Summary

The paper is a mathematical proof that in a hierarchical forest (where big trees shade small ones):

If the trees are fertile enough to survive in an empty forest, the population will grow until it hits a limit.
There is only one specific, stable population size the forest can settle into.
Once it settles there, it is stable; small bumps in the road won't destroy the forest, it will just wobble and return to its rhythm.

They did this by switching from a complicated "map of the whole forest" view to a simpler "counting the babies" view, proving that both views tell the same story.

Technical Summary: Hierarchical Competition through Reduction of Individual Growth

Problem Statement
The paper investigates the dynamics of a size-structured population where individuals are organized hierarchically based on size (interpreted as height). The core biological mechanism is that the growth rate of an individual depends exclusively on the presence of larger individuals, which exert a negative impact (e.g., shading in a forest context). This creates a feedback loop where the environmental condition (light intensity) is determined by the population's size distribution, which in turn dictates individual growth.

The authors aim to analyze the well-posedness, existence, and stability of stationary states for this system. While the classic formulation of such models employs a first-order quasi-linear partial differential equation (PDE) with non-local boundary conditions, the authors note that rigorous stability results (specifically the Principle of Linearised Stability) have not been fully established for the PDE formulation in this specific context. Consequently, the paper focuses on an alternative formulation: a scalar nonlinear renewal equation (a delay equation) for the population birth rate.

Methodology
The study proceeds through the following methodological steps:

Model Formulation:
- PDE Approach: The authors first present the standard PDE formulation involving the population density $u(x,t)$ , growth rate $g(E(x,t))$ , mortality $\mu$ , and fertility $\beta$ . The interaction variable $E(x,t)$ represents the cumulative density of individuals larger than size $x$ .
- Delay Formulation: By tracking cohorts and utilizing the relationship between age and size in a hierarchical system (where larger individuals are necessarily older), the authors derive a scalar renewal equation for the birth rate $B(t)$ . This equation expresses $B(t)$ as an integral over past birth rates, weighted by survival probabilities and fertility functions dependent on the size attained at a given age.
Dynamical Systems Framework:
- The delay equation is treated as an abstract Cauchy problem in a Banach space $X = L^1_\rho(-\infty, 0)$ of weighted birth rate histories.
- The authors establish the well-posedness of the problem by proving that the nonlinear operator defining the delay equation is continuously differentiable (Theorem 1), allowing the application of the Principle of Linearised Stability (PLS) for delay equations.
Steady State Analysis:
- The existence of non-trivial stationary birth rates is analyzed by reducing the delay equation to a fixed-point problem $B = R(B)$ , where $R$ is the net reproduction function.
- Conditions are derived for the existence of a unique positive equilibrium, primarily relying on the basic reproduction number $R_0$ (the expected offspring in a virgin environment) and the behavior of the growth and fertility functions at extreme population densities.
Stability Analysis:
- Trivial Equilibrium: The stability of the zero state (extinction) is determined by the sign of the real part of the roots of the characteristic equation associated with the linearization at $B=0$ .
- Non-Trivial Equilibrium: The authors analyze the local asymptotic stability of the positive steady state. They employ a two-step strategy:
  - First, they explicitly solve the characteristic equation for a piecewise linear fertility function ( $\beta_0$ ), proving that all roots lie in the left half-plane under specific conditions.
  - Second, they use Rouché's Theorem to show that this stability property is preserved for a class of smooth fertility functions that are sufficiently close to the piecewise linear case (Theorem 8).

Key Contributions and Results

Well-Posedness of the Delay Formulation: The paper rigorously proves that the scalar renewal equation formulation is well-posed and generates a nonlinear semigroup, provided the growth and fertility functions have bounded, globally Lipschitzian derivatives (Theorem 2).
Uniqueness of Non-Trivial Equilibrium: Under general assumptions (strictly increasing fertility with size, strictly decreasing growth with competition), the paper demonstrates that at most one non-trivial stationary birth rate exists (Theorem 4).
Stability Criteria:
- The trivial equilibrium is locally exponentially stable if $R_0 < 1$ and unstable if $R_0 > 1$ (Theorem 5).
- Crucially, the paper establishes that if a non-trivial equilibrium exists, it is locally asymptotically stable for a broad class of fertility functions, including smooth perturbations of piecewise linear functions (Theorem 8). This is achieved by proving that the characteristic equation has no roots with non-negative real parts.
Equivalence of Formulations: The authors show in Appendix B that the conditions guaranteeing the existence of stationary population densities in the PDE formulation coincide exactly with the conditions for non-trivial stationary birth rates in the delay formulation.
Semi-Explicit Solutions: For a specific case where fertility is proportional to size and growth follows a Holling Type II-like function, the authors derive explicit expressions for the equilibrium birth rate (using the Lambert W function) and the steady-state size distribution (Section 5.2).

Significance and Claims
The paper claims to provide a rigorous mathematical foundation for the stability of hierarchical population models, a topic where previous PDE-based analyses were often formal or lacked rigorous justification for the Principle of Linearised Stability.

By shifting the focus to the delay formulation, the authors successfully apply the PLS to prove that, under reasonable and general biological assumptions, the system converges to a unique, stable non-trivial equilibrium whenever it exists. The authors explicitly state that while they do not rigorously translate these stability results back to the PDE formulation within this work, the phenomenological equivalence of the two formulations suggests that the stability results for the delay equation likely hold for the PDE as well. The work is dedicated to Professor Glenn F. Webb, acknowledging his foundational contributions to structured population dynamics.

The paper does not propose new experimental applications or future directions beyond the mathematical analysis of the specific model presented. It remains a theoretical contribution to the field of physiologically structured population dynamics.