Existence, Stability and Slow Dynamics of Spikes in a 1D Minimal Keller--Segel Model with Logistic Growth

TL;DR

This paper investigates the existence, stability, and slow dynamics of localized spike patterns in a 1D Keller--Segel chemotaxis model with logistic growth, using a novel singular limit approach based on small chemoattractant diffusivity.

Contribution

It introduces a new analytical framework for spike pattern analysis in Keller--Segel models with small chemoattractant diffusivity, including stability criteria and slow dynamics derivation.

Findings

N-spike equilibria can be destabilized by zero-eigenvalue crossing or Hopf bifurcation.

Explicit stability range for cellular diffusion rate d_1 is identified.

Derived a differential algebraic system governing slow spike dynamics.

Abstract

We analyze the existence, linear stability, and slow dynamics of localized 1D spike patterns for a Keller--Segel model of chemotaxis that includes the effect of logistic growth of the cellular population. Our analysis of localized patterns for this two-component reaction-diffusion (RD) model is based, not on the usual limit of a large chemotactic drift coefficient, but instead on the singular limit of an asymptotically small diffusivity $d_2=\epsilon^2\ll 1$ of the chemoattractant concentration field. In the limit $d_2\ll 1$, steady-state and quasi-equilibrium 1D multi-spike patterns are constructed asymptotically. To determine the linear stability of steady-state $N$-spike patterns we analyze the spectral properties associated with both the ''large'' ${\mathcal O}(1)$ and the ''small'' $o(1)$ eigenvalues associated with the linearization of the Keller--Segel model. By analyzing a nonlocal eigenvalue problem characterizing the large eigenvalues, it is shown that $N$-spike equilibria can be destabilized by a zero-eigenvalue crossing leading to a competition instability if the cellular diffusion rate $d_1$ exceeds a threshold, or from a Hopf bifurcation if a relaxation time constant $\tau$ is too large. In addition, a matrix eigenvalue problem that governs the stability properties of an $N$-spike steady-state with respect to the small eigenvalues is derived. From an analysis of this matrix problem, an explicit range of $d_1$ where the $N$-spike steady-state is stable to the small eigenvalues is identified. Finally, for quasi-equilibrium spike patterns that are stable on an ${\mathcal O}(1)$ time-scale, we derive a differential algebraic system (DAE) governing the slow dynamics of a collection of localized spikes. Unexpectedly, our analysis is rather closely related to the analysis of spike patterns for the Gierer--Meinhardt RD system.