Asymptotic stability and diffusion-driven pattern formation in a predator-prey system with two chemicals

1. Problem Statement

The paper investigates a four-component reaction-diffusion system modeling a predator-prey interaction mediated by two chemical substances. Unlike traditional models that rely solely on self-diffusion, this system incorporates cross-diffusion driven by chemical signaling:

• Prey ($u_2$) release a chemical attractant ($v_1$) that induces directed movement (advection) in the predator ($u_1$).
• Predators ($u_1$) release a chemical repellent ($v_2$) that induces directed movement away from the source in the prey.

The system is defined on a bounded domain $\Omega \subset \mathbb{R}^n$ ($n \ge 2$) with homogeneous Neumann boundary conditions. The primary mathematical challenges addressed are:

1. Proving the global existence and boundedness of classical solutions for this coupled parabolic system.
2. Establishing the global asymptotic stability of spatially homogeneous steady states.
3. Analyzing Turing instability (diffusion-driven pattern formation) induced specifically by the cross-diffusion terms, which are known to be mathematically complex in four-equation systems.

2. Mathematical Model

The system is given by:
$$\begin{cases} u_{1t} = \nabla \cdot (d_{11}\nabla u_1 + d_{12}\nabla v_1) + u_1(\mu_1 - \lambda_1 u_1 + \eta_1 u_2) \\ u_{2t} = \nabla \cdot (d_{21}\nabla u_2 - d_{22}\nabla v_2) + u_2(\mu_2 - \lambda_2 u_2 - \eta_2 u_1) \\ v_{1t} = d_3 \Delta v_1 + a_1 u_2 - b_1 v_1 \\ v_{2t} = d_4 \Delta v_2 + a_2 u_1 - b_2 v_2 \end{cases}$$
Where:

• $u_1, u_2$: Predator and prey densities.
• $v_1, v_2$: Concentrations of chemicals produced by prey and predator, respectively.
• $d_{ij}$: Self-diffusion ($d_{11}, d_{21}, d_3, d_4$) and cross-diffusion ($d_{12}, d_{22}$) coefficients.
• $\eta_2$: Predation rate (interaction strength).
• Parabolicity Condition: The authors assume $4d_{11}d_3 > d_{12}^2$ and $4d_{21}d_4 > d_{22}^2$ to ensure the system remains parabolic and well-posed.

3. Methodology

The authors employ a rigorous combination of analytical techniques and numerical simulations:

• Local Existence: Utilizing quasi-linear parabolic theory (Amann's theory) to establish the existence of a unique local classical solution.
• Global Existence & Boundedness:
  • Non-negativity: Proven using the negative part method and Young's inequality.
  • $L^1$ and $L^\infty$ Bounds: Derived via integration of the equations, Cauchy-Schwarz inequalities, and the Moser-Alikakos iteration technique. This step is crucial to prevent blow-up in finite time.
• Asymptotic Stability:
  • Constructed Lyapunov functionals ($E_1(t)$ and $E_2(t)$) tailored to the coexistence and prey-vanishing states.
  • Used Taylor expansions and Cauchy-Schwarz inequalities to show that the time derivative of the Lyapunov functional is negative definite ($\frac{d}{dt}E(t) \le -\epsilon f(t)$), proving exponential convergence to equilibrium.
• Pattern Formation (Turing Instability):
  • Performed linear stability analysis on the spatially homogeneous steady state.
  • Derived the characteristic polynomial of the linearized system. Due to the high complexity of the four-equation system, explicit analytical conditions for Turing bifurcation were deemed intractable.
  • Instead, the authors used numerical simulations (Finite Difference Method) to identify parameter regimes where the real part of the eigenvalues becomes positive, indicating instability.

4. Key Contributions and Results

A. Global Existence and Stability

• Theorem 2.1: Proved that for any non-negative initial data, the system admits a unique, globally bounded classical solution.
• Theorem 2.2 (Coexistence): If the predation rate $\eta_2$ is low ($\eta_2 < \frac{\mu_2 \lambda_1}{\mu_1}$) and cross-diffusion coefficients are sufficiently small, the system converges exponentially to a unique positive coexistence equilibrium.
• Theorem 2.3 (Prey Extinction): If the predation rate is high ($\eta_2 \ge \frac{\mu_2 \lambda_1}{\mu_1}$), the system converges to a semi-trivial equilibrium where the prey goes extinct, and the predator survives on alternative resources (or the model's specific growth terms).

B. Pattern Formation Dynamics

The study reveals that cross-diffusion is the primary driver of spatial patterns in this system; self-diffusion alone maintains stability.

• Numerical Bifurcation: The authors generated Turing bifurcation diagrams by varying the predation rate ($\eta_2$) and interaction parameters ($\eta_1$).
• Pattern Types Observed:
  • Labyrinthine Stripes: Observed at moderate interaction strengths.
  • Hexagonal Spots: Emerged as predation intensity increased, indicating localized aggregation.
  • Mixed Stripe-Spot: Transitions between different morphologies.
  • Irregular/Disordered Patterns: Observed near instability boundaries or with mixed diffusion coefficients.
• Synchronization: A unique finding is that predator and prey densities often exhibit synchronized spatial patterns (peaks and troughs align). This contrasts with classical predator-prey models where they are often out of phase. This synchronization is attributed to the shared chemical cues ($v_1$ and $v_2$) acting as environmental drivers for both species.

C. Parameter Sensitivity

• Predation Rate ($\eta_2$): Increasing $\eta_2$ drives a transition from extended stripes to localized hexagonal spots, suggesting stronger prey avoidance leads to sharper spatial segregation.
• Chemical Production ($\eta_1$): Higher values promote the formation of isolated coexistence spots (territorial behavior).
• Growth/Decay Rates: Low growth rates lead to fine-grained, fragmented microstructures.

5. Significance

• Theoretical Advancement: The paper provides a rigorous proof of global existence for a complex 4-component cross-diffusion system, a class of problems often difficult to analyze due to the lack of maximum principles in higher dimensions.
• Ecological Insight: It demonstrates that chemical mediation fundamentally alters spatial dynamics. The "chemical cross-diffusion" mechanism allows for pattern formation even when classical diffusion would predict a homogeneous state.
• Biological Relevance: The results offer a mathematical explanation for observed biological phenomena, such as bacterial-bacteriophage interactions or algae-zooplankton dynamics, where chemical signaling dictates spatial distribution.
• Control Mechanisms: The identification of stability boundaries suggests that manipulating predation rates or chemical decay rates could theoretically control ecosystem stability, preventing chaotic oscillations or extinction events.

In conclusion, this work bridges the gap between rigorous PDE analysis and ecological modeling, proving that cross-diffusion mediated by chemical signals is a potent mechanism for generating complex, stable, and synchronized spatial patterns in predator-prey ecosystems.