Keller–Segel PDE with Logistic Damping

• The Keller–Segel PDE with logistic damping is characterized by a balance between diffusion, chemotactic drift, and nonlinear logistic decay, leading to precise blow-up and boundedness criteria.
• It illustrates that in 2D any positive logistic damping ensures global boundedness, while in higher dimensions, stricter damping parameters are required to prevent finite-time blow-up.
• Methodologies such as mass-accumulation transformation and energy inequalities are key tools used to derive sharp thresholds and to understand the dynamics of chemotactic aggregation.

The Keller–Segel PDE with logistic damping is a canonical model for chemotactic aggregation incorporating nonlinear density-dependent decay. It exhibits a rich interplay between diffusion, chemotactic drift, and nonlinear population regulation. The classic Keller–Segel dynamics can feature finite-time blow-up—unbounded density growth—depending on system parameters, dimension, and initial data. Introduction of a logistic (density-limiting) sink fundamentally alters this balance, modifying global existence thresholds and, in some scenarios, entirely suppressing blow-up. Recent analyses have established sharp criteria for global-in-time regularity, precise blow-up thresholds, and the critical role of logistic nonlinearity in both parabolic–elliptic and fully parabolic variants across dimensions.

1. Mathematical Formulation and Variants

The standard Keller–Segel system with logistic damping, posed in a domain $\Omega \subset \mathbb R^n$, $n \geq 1$, takes the form

$$\begin{cases} u_t = \Delta u - \nabla\cdot(u \nabla v) + \lambda u - \mu u^\kappa, & \text{in } \Omega,\, t>0 \ 0 = \Delta v - \overline{m}(t) + u, & \text{in } \Omega,\, t>0 \ \text{Neumann }\, \partial_\nu u = \partial_\nu v = 0, & \text{on } \partial\Omega \ u(x,0) = u_0(x)\geq 0 \end{cases}$$

where $u$ is the cell density, $v$ the chemoattractant, $\lambda$ a linear growth rate, $\mu > 0$ the logistic damping coefficient, and $\kappa > 1$ the exponent of the nonlinear sink (Fuest, 2020). The mass-control condition $$\overline{m}(t) = \frac{1}{|\Omega|}\int_\Omega u(x,t)\,dx$$ enforces zero spatial mean for $v$.

Extensions include higher-dimensional, quasi- or fully parabolic settings (with $v_t$ instead of $v$ replacing the elliptic equation), generalized nonlinear diffusions, and the inclusion of source/sink terms or reaction kinetics beyond quadratic decay. The logistic term is biologically motivated as a combination of resource limitation and crowding-induced mortality.

2. Global Existence vs Finite-Time Blow-Up: Sharp Criteria

Dimensional Thresholds and Damping Exponents

The well-posedness and regularity theory is governed by competing mechanisms:

• Diffusive regularization ($\Delta u$)
• Aggregative chemotaxis ($-\nabla\cdot(u\nabla v)$)
• Logistic sink ($-\mu u^\kappa$)

Key results (Fuest, 2020, Liu et al., 16 Apr 2025, Cavallazzi et al., 23 Dec 2025):

• Global boundedness: For all $n\ge2$, if $\kappa>2$, all classical solutions are global and uniformly bounded.
• Critical quadratic damping: For $n\ge4$, $\kappa=2$ is the sharp threshold. If $\kappa<2$ or, for $\kappa=2$, $\mu<(n-4)/n$ ($n\ge5$), radial initial data can lead to finite-time blow-up.
• 3D case: For $n=3$, blow-up exists if $\kappa \in (1, \frac{3}{2})$ ($\lambda\ge0$, any $\mu>0$) (Fuest, 2020). For classical quadratic damping term $-\mu u^2$, the sharp threshold is $\mu=1/3$; for $\mu<1/3$ finite-time blow-up is provable, whereas for $\mu\ge 1/3$ all solutions are global and bounded (Liu et al., 16 Apr 2025).

Main Theorem (Sharpness)

For $n\ge3$, pick $\kappa\in(1,2)$ ($n\ge4$), or $\kappa=2$ with $\mu \in (0, (n-4)/n)$ ($n\ge5$), or $\kappa \in (1,3/2)$ ($n=3$). For any $\lambda\ge0$, there exists radial, monotonic $u_0$ with $\int_\Omega u_0>0$ such that (KS$_\mathrm{log}$) exhibits finite-time blow-up (Fuest, 2020, Liu et al., 16 Apr 2025).

These thresholds reflect the failure of the damping term to offset the non-linearity of the chemotactic drift—leading to superlinear ODE growth in appropriate functionals.

3. Analytical Framework and Key Methodologies

The analysis exploits reduction to a mass-accumulation function for radial solutions, $w(s,t)$, representing accumulated mass up to radius $s$, and transforms the PDE system into an equivalent parabolic equation for $w$ (with superlinear source), together with crucial monotonicity and upper-lower bounds: $$w_s(s,t) = n\,u(s^{1/n}, t), \qquad w(s,t)\leq s\,w_s(s,t)$$ Functional inequalities for $w$ facilitate the construction of test functionals $\phi(t)$ that, via ODE comparison arguments, demonstrate finite-time blow-up when the superlinear source dominates dissipative contributions (see Lemmas 3.1, 3.2, 4.1 in (Fuest, 2020)).

In higher dimensions, blow-up is controlled by the mass-accumulation ODE when $\kappa$ falls below critical. For quadratic damping, exact thresholds relate to the parameter $\mu$ and dimension (Liu et al., 16 Apr 2025, Cavallazzi et al., 23 Dec 2025).

In boundedness regimes, the approach relies on a suite of a priori estimates (energy identities, entropy functions, interpolation inequalities such as Gagliardo–Nirenberg, de la Vallée–Poussin conditions), Moser iteration, and localized or mass-dissipation functionals (Xiang, 2017, Le, 2023, Silva et al., 2022).

4. Effects of Logistic Damping and Extensions

2D Case: Strength of Damping

In two dimensions, any positive logistic damping ($\mu > 0$) ensures global boundedness, regardless of the chemotactic strength or initial data—demonstrating the absolute regularizing effect of quadratic logistic decay (Jin et al., 2018, Xiang, 2017). Moreover, sub-logistic damping of the form $-\mu u^2/\log^p(u+e)$ ($p\le 1$) remains sufficient to prevent blow-up in 2D, even though this is strictly weaker than quadratic decay (Le, 2023, Le, 2023). The criticality for blow-up in higher dimensions is not observed in $n=2$.

Higher Dimensions and Quasilinear Extensions

In $n\ge3$, suppression of blow-up by logistic damping requires a "sufficiently strong" coefficient, with explicit relationships dictated by system parameters (diffusion scaling, chemotactic sensitivity). For quasilinear models with nonlinear diffusion $\phi(u)\sim u^p$ or generalized logistic forms $-\mu u^r$, sharp parameter-dependent thresholds for boundedness vs blow-up can be established. For example, strong enough damping ($\mu>\frac{d-2}{d}\chi$ for the parabolic–elliptic system in $\mathbb R^d$) ensures global regularity (Cavallazzi et al., 23 Dec 2025), and more general results are available for cross-diffusion and attraction-repulsion systems (Chiyo et al., 2022). In the subquadratic case ($1 < r < 2$), the chemotactic sensitivity must be restricted in terms of initial data and the nonlinear exponent (Kang et al., 2024).

5. Weak Solutions, Control, and Numerical Aspects

The existence of global bounded weak solutions is established in a variety of frameworks, including unbalanced optimal transport splitting schemes, which yield explicit parameter-dependent thresholds in the subquadratic regime (Kang et al., 2024). Control-theoretic extensions consider optimal control problems for the chemically-mediated dynamics, with the logistic damping ensuring well-posedness even in weak or very-weak settings and under singular controls (Silva et al., 2022).

Numerical investigation of the traveling-wave regime and pattern formation under strong aggregation reveals that, in 1D, invasion speed remains at the classical Fisher–KPP value irrespective of chemotactic strength, as long as logistic damping is present (Henderson et al., 2023).

6. Blow-Up Mechanisms and New Phenomena

Finite-time blow-up in supercritical or marginally critical regimes is constructed via explicit self-similar solutions, abnormal blow-up rates, and stability analysis of perturbations in weighted $L^2$ spaces (Liu et al., 16 Apr 2025). For the 3D system with quadratic damping, type I blow-up profiles can exist for all $\mu < 1/3$, and a countable family of subcritical law solutions is constructed, clarifying the precise failure of logistic control in certain parameter domains. The analysis exploits stability against perturbations in singular-weighted spaces and the modulation method for unstable directions of the linearized operator.

7. Summary Table: Critical Regimes for Keller–Segel with Logistic Damping

Dimension $n$	Damping Exponent $\kappa$	Sharp Bound on $\mu$	Blow-up/Global Solution
$n=2$	any $\kappa>1$	any $\mu>0$	All solutions global and bounded
$n=3$	$\kappa=2$	$\mu = 1/3$	$\mu<1/3$: blow-up; $\mu\ge 1/3$: global for all initial data
$n\ge4$	$\kappa=2$	$\mu = (n-4)/n$	$\mu < (n-4)/n$: blow-up; $\mu\ge (n-4)/n$: global
$n\ge3$	$\kappa>2$	any $\mu>0$	All solutions global and bounded
$n\ge4$	$1<\kappa<2$	any $\mu>0$	Existence of blow-up solutions

Interpretation: Subquadratic (or critical quadratic) damping is insufficient to prevent blow-up in dimensions $n\ge3$, unless the damping parameter exceeds a sharp threshold. In $n=2$, any logistic damping suffices, and even strictly subquadratic or logarithmically weakened terms ensure boundedness (Fuest, 2020, Liu et al., 16 Apr 2025, Xiang, 2017, Le, 2023).

---

References:

Fuest, J. "Approaching optimality in blow-up results for Keller–Segel systems with logistic-type dampening" (Fuest, 2020). Jin, H.-Y., Xiang, T. "Chemotaxis effect vs logistic damping on boundedness in the 2-D minimal Keller-Segel model" (Jin et al., 2018). Le, M. "Blow-up prevention by sub-logistic sources in Keller-Segel cross diffusion type system" (Le, 2023). Wu, K. et al. "Finite time blowup for Keller-Segel equation with logistic damping in three dimensions" (Liu et al., 16 Apr 2025). Fournier, N., Mishura, Y. "Quantitative approximation of a Keller–Segel PDE by a branching moderately interacting particle system and suppression of blow-up" (Cavallazzi et al., 23 Dec 2025). Xiang, T. "Sub-logistic source can prevent blow-up in the 2D minimal Keller-Segel chemotaxis system" (Xiang, 2017). Braz e Silva et al. "Bilinear optimal control for weak solutions of the Keller-Segel logistic model in $2D$ domains" (Silva et al., 2022). Yang, X.F. et al. "Boundedness in a quasilinear fully parabolic Keller-Segel system of higher dimension with logistic source" (Yang et al., 2015). Cañizo, J.A. et al. "Bounded weak solutions for Keller-Segel equations with generalized diffusion and logistic source via an unbalanced Optimal Transport splitting scheme" (Kang et al., 2024). Chiyo, T., Frassu, S., Viglialoro, G. "A nonlinear attraction–repulsion Keller–Segel model with double sublinear absorptions: criteria toward boundedness" (Chiyo et al., 2022).