🔢 mathematics

Uniform boundedness for the two-dimensional Keller-Segel system with Gompertz growth

This paper investigates a minimal two-dimensional Keller-Segel system with Gompertz-type growth, demonstrating that this weaker damping mechanism, compared to classical logistic terms, is sufficient to prevent cell aggregation and ensure the global existence and uniform boundedness of solutions under suitable initial conditions.

Original authors: Nohayla Alaoui, Mohamed Halloumi, Giuseppe Viglialoro

Published 2026-03-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Nohayla Alaoui, Mohamed Halloumi, Giuseppe Viglialoro

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a crowded dance floor where people (cells) are moving around. In this specific scenario, the people have a strange superpower: they can smell a scent that they themselves are producing. If they smell a lot of it, they rush toward the source. This is called chemotaxis.

In the world of mathematics, this is modeled by the Keller–Segel system. The problem is that this "smell and rush" behavior can get out of hand. If too many people rush to the same spot, they might pile up so high that the model predicts an infinite crowd in a finite amount of time. In math terms, this is called a "blow-up." It's like a traffic jam so severe that the cars literally crush into a single point.

The Problem: How do we stop the pile-up?

Scientists have known for a long time that if you add a "logistic source" (a rule that says "people die or stop reproducing when it gets too crowded"), the pile-up stops. Think of this like a bouncer at a club who kicks people out once the room gets too full. The classic bouncer rule is simple: "If the crowd is $X$ , the death rate is $X^2$ ." This is a very strong bouncer; it stops the crash easily.

But, nature isn't always that simple. Sometimes, the "bouncer" is weaker. Maybe the death rate doesn't spike as fast as the crowd grows. The big question for mathematicians is: How weak can the bouncer be before the crowd still crashes?

The New Discovery: The Gompertz "Bouncer"

This paper introduces a new type of bouncer called the Gompertz growth function.

The Old Bouncer (Logistic): Imagine a bouncer who gets angry and kicks people out very aggressively as the room fills up.
The New Bouncer (Gompertz): This bouncer is more subtle. They don't just kick people out based on the raw number; they kick them out based on how close the room is to its maximum capacity, but in a way that slows down the growth rate exponentially. It's like a bouncer who says, "We are getting full, so let's slow down the entry and gently encourage people to leave," rather than "Get out now!"

The authors found that this Gompertz bouncer is actually weaker than a previous "sub-logistic" bouncer that was studied in 2018. It's a much gentler hand trying to control the crowd.

The Big Question: Can a weak bouncer save the day?

The researchers asked: "If our bouncer is this gentle (Gompertz), can we still prevent the infinite pile-up (blow-up) in a 2D room?"

The Answer is Yes, but with conditions.

They proved that the crowd will stay safe and never crash, provided two things happen:

The room isn't too small: The "carrying capacity" (how many people the room can theoretically hold) must be large enough.
The smell isn't too strong: The chemical signal that pulls people together (the "attraction") must be weak enough compared to the size of the room.

If these conditions are met, even this gentle Gompertz bouncer is strong enough to keep the dance floor from collapsing into a singularity.

The "Secret Sauce" of the Proof

How did they prove this? They used a clever mathematical trick involving a "safety net."

The Energy Function: They created a mathematical "score" for the system. This score combines how crowded the room is and how much the people are moving around.
The Trap: They showed that if the crowd starts to get too wild, this "score" would naturally try to decrease (like a ball rolling down a hill).
The Catch: They proved that as long as the "smell" (attraction) isn't too strong, the gentle Gompertz bouncer is strong enough to keep that "score" from ever getting too high. If the score stays low, the crowd stays safe.

Why Does This Matter?

This isn't just about math puzzles. This model is used to understand cancer.

Tumor cells move toward chemical signals they produce.
Tumor growth often follows the Gompertz pattern (they grow fast at first, then slow down as they run out of space/nutrients).

By proving that this specific, biologically realistic growth pattern prevents the mathematical "crash," the authors give us more confidence that we can model tumor behavior accurately without the math breaking down. It tells us that even with a "gentle" natural limit on growth, the system remains stable, which is a huge relief for scientists trying to predict how tumors spread.

In a Nutshell

Think of the Keller–Segel system as a game of musical chairs where the chairs are moving toward the music.

Without a bouncer: Everyone piles into one chair, and the game breaks.
With a strong bouncer (Logistic): The bouncer kicks people out fast enough to save the game.
With the Gompertz bouncer (This paper): The bouncer is much gentler, almost whispering "slow down." The authors proved that even this whisper is enough to keep the game going forever, as long as the music isn't too loud and the room is big enough.

This is a victory for "gentle" mathematics, showing that you don't always need a heavy hand to keep a complex system from falling apart.

1. Problem Statement

The paper investigates the global existence and uniform boundedness of solutions to the two-dimensional Keller–Segel chemotaxis system coupled with a Gompertz growth source. The classical Keller–Segel model describes the movement of cells ( $u$ ) toward a chemical signal ( $v$ ) they produce. In two dimensions, this system is prone to "finite-time blow-up" (cell aggregation leading to infinite density) if the total initial mass exceeds a critical threshold.

To prevent this, researchers often introduce a logistic source term ( $au - bu^2$ ) representing population growth and death. It is well-established that even a small quadratic damping term ( $-bu^2$ ) is sufficient to ensure global boundedness in 2D.

The specific problem addressed here:
The authors replace the standard logistic source with a Gompertz source term:
$f(u) = \alpha u \ln\left(\frac{K}{u}\right)$
where $\alpha > 0$ is the intrinsic growth rate and $K > 0$ is the carrying capacity.

Biological Context: The Gompertz model is widely used in cancer biology to describe tumor growth, as it better captures the asymmetric deceleration of growth rates compared to the symmetric logistic model.
Mathematical Challenge: The Gompertz term provides a weaker damping effect than the quadratic logistic term and even weaker than the "sub-logistic" sources ( $as - \frac{bs^2}{\ln(\ln(s+e))}$ ) analyzed in previous works (e.g., Xiang, 2018). Specifically, for large $u$ , the Gompertz term decays slower than the quadratic term, making the prevention of blow-up more difficult.

The system under study is:
$\begin{cases} u_t = \Delta u - \chi \nabla \cdot (u \nabla v) + \alpha u \ln\left(\frac{K}{u}\right), \\ \tau v_t = \Delta v - v + u, \end{cases}$
on a bounded smooth domain $\Omega \subset \mathbb{R}^2$ , with homogeneous Neumann boundary conditions. Here, $\tau \in \{0, 1\}$ distinguishes between the parabolic-elliptic ( $\tau=0$ ) and fully parabolic ( $\tau=1$ ) cases.

2. Methodology

The authors employ a combination of energy methods, comparison principles, and Gagliardo–Nirenberg inequalities to establish a priori estimates.

A. Local Existence and Extension Criteria

Local Well-posedness: Using standard fixed-point arguments, they establish the existence of a unique local classical solution $(u, v)$ on a maximal time interval $[0, T_{max})$ .
Blow-up Criterion: They refine the standard blow-up criterion. They prove that if the solution blows up, it must be due to the $L^\infty$ norm of $u$ becoming unbounded. Crucially, they use a barrier method to show that $u$ remains strictly positive (bounded away from zero) as long as it is bounded from above, ruling out singularities caused by $u \to 0$ .

B. The Boundedness Criterion

The core strategy relies on Lemma 2.2, which states that if the following two quantities are uniformly bounded in time, then the solution is globally bounded:

The entropy-like functional: $\int_\Omega u \ln(u) \, dx$
The gradient energy of the signal: $\int_\Omega |\nabla v|^2 \, dx$

C. A Priori Estimates

To prove the boundedness of the quantities above, the authors define a Lyapunov-type functional:
$F(t) = \int_\Omega u \ln(u) \, dx + \tau b \int_\Omega |\nabla v|^2 \, dx$
where $b = \chi/2$ .

Mass Control: They first establish that the total mass $\int_\Omega u$ is bounded by $M = \max\{\int u_0, K|\Omega|\}$ using Jensen's inequality on the concave Gompertz function.
Differential Inequality: By differentiating $F(t)$ with respect to time and applying integration by parts, Young's inequality, and the Gagliardo–Nirenberg inequality, they derive a differential inequality of the form:
$F'(t) \leq C_1 - C_2 F(t)$
This inequality holds provided specific conditions on the parameters are met.
Handling the Parabolic-Elliptic Case ( $\tau=0$ ): When $\tau=0$ , the term involving $\nabla v$ vanishes from $F(t)$ . The authors must separately estimate $\int |\nabla v|^2$ . They achieve this by multiplying the $u$ -equation by $u$ , deriving an $L^2$ estimate for $u$ , and using elliptic regularity to bound $\nabla v$ .

3. Key Contributions

Extension of Damping Theory: The paper demonstrates that the Gompertz growth term, despite being a weaker damping mechanism than the quadratic logistic term or the sub-logistic terms previously studied, is still sufficient to prevent blow-up in 2D Keller–Segel systems.
Novel Conditions for Global Existence: The authors derive explicit conditions on the parameters (specifically the carrying capacity $K$ and the chemotactic sensitivity $\chi$ ) that guarantee global boundedness.
Unified Treatment: The analysis covers both the fully parabolic ( $\tau=1$ ) and parabolic-elliptic ( $\tau=0$ ) cases, providing a comprehensive result for the minimal chemotaxis model.

4. Main Results

Theorem 1.1: Under the assumption that the domain $\Omega$ is $C^{2+\delta}$ and the initial data is sufficiently regular, the system admits a unique global classical solution $(u, v)$ that is uniformly bounded in $L^\infty(\Omega \times (0, \infty))$ , provided that:

Carrying Capacity Condition: $K > e^{-2/\alpha}$ (This ensures the growth term is sufficiently strong/negative for large $u$ ).
Chemotaxis Sensitivity Condition: $\chi M \leq \frac{1}{2 C_{GN}^4}$ , where $M$ is the bound on the total mass and $C_{GN}$ is the Gagliardo–Nirenberg constant.

Implication: If the chemotactic sensitivity $\chi$ is small enough relative to the carrying capacity and initial mass, the Gompertz damping prevents the aggregation instability typical of the classical Keller–Segel model.

5. Significance

Biological Relevance: By validating the Gompertz model within the chemotaxis framework, the paper bridges a gap between mathematical theory and biological reality, particularly in modeling tumor growth where Gompertz kinetics are more accurate than logistic ones.
Mathematical Rigor: The work pushes the boundaries of known damping mechanisms. It shows that even "weak" damping (where the mortality rate grows slower than $u^2$ ) can stabilize the system in 2D, provided the parameters satisfy specific thresholds.
Future Directions: This analysis opens the door to studying other non-logistic, biologically motivated source terms in chemotaxis models, potentially leading to more realistic simulations of cell migration and tumor invasion.

In summary, the paper proves that the Gompertz growth term acts as a stabilizer for the 2D Keller–Segel system, ensuring that cell densities remain bounded and do not collapse into singularities, provided the chemotactic attraction is not too strong relative to the growth parameters.