2025-11-18T00:37:13.255978

The simultaneous effect of chemotaxis and alarm-taxis on the global existence and stability of a predator-prey system

Shanmugasundaram, Saha, Fuentes
This study examines a fully parabolic predator-prey chemo-alarm-taxis system under homogeneous Neumann boundary conditions in a bounded domain $Ω\subset \mathbb{R}^n$ with a smooth boundary $\partialΩ$. Under specific parameter conditions, it is shown that the system admits a unique, globally bounded classical solution. The convergence of the solution is established through the construction of an appropriate Lyapunov functional. In addition, numerical simulations are presented to validate the asymptotic behaviour of the solution. The results highlight the significant role of chemotaxis and alarm-taxis coefficients in determining the existence and stability of predator-prey models, as discussed in the literature.
academic

The simultaneous effect of chemotaxis and alarm-taxis on the global existence and stability of a predator-prey system

Basic Information

  • Paper ID: 2510.14728
  • Title: The simultaneous effect of chemotaxis and alarm-taxis on the global existence and stability of a predator-prey system
  • Authors: Gnanasekaran Shanmugasundaram (NIT Tiruchirappalli), Jitraj Saha (NIT Tiruchirappalli), Rafael Díaz Fuentes (University of Cagliari)
  • Classification: math.AP (Analysis of PDEs)
  • Submission Date: October 16, 2025
  • Paper Link: https://arxiv.org/abs/2510.14728

Abstract

This study examines a fully parabolic predator-prey chemo-alarm-taxis system with homogeneous Neumann boundary conditions on a bounded domain Ω⊂Rn\Omega \subset \mathbb{R}^n with smooth boundary ∂Ω\partial\Omega. Under specific parameter conditions, the paper proves the existence of a unique global bounded classical solution to the system. The convergence of solutions is established through the construction of appropriate Lyapunov functionals. Furthermore, the asymptotic behavior of solutions is verified through numerical simulations. The results emphasize the crucial role of chemotaxis and alarm-taxis coefficients in determining the existence and stability of predator-prey models.

Research Background and Motivation

Biological Background

  1. Alarm Call Mechanisms: The research is based on the "thief alarm" hypothesis in marine ecosystems, wherein prey emit signals to attract secondary predators, thereby deterring primary predators. A typical example is the bioluminescent signals emitted by dinoflagellates in response to copepod disturbance, attracting fish as secondary predators.
  2. Chemotaxis Phenomena: The movement behavior of microorganisms in response to chemical gradients, which has widespread applications in ecology, immunology, biology, and medicine.

Mathematical Modeling Requirements

The core problems this research addresses are:

  • How to mathematically model a three-species interaction system simultaneously incorporating chemotaxis and alarm-taxis
  • Under what conditions the solutions of such a system exist globally and remain bounded
  • What the long-term asymptotic behavior of the system is

Limitations of Existing Research

Existing literature primarily focuses on:

  • Studies of single taxis mechanisms (pure chemotaxis or pure alarm-taxis)
  • Local results in two-dimensional spaces
  • Lack of rigorous mathematical analysis simultaneously considering both taxis mechanisms

Core Contributions

  1. First establishment of rigorous mathematical theory for global existence, boundedness, and convergence of a predator-prey system simultaneously incorporating chemotaxis and alarm-taxis
  2. Construction of innovative Lyapunov functionals to analyze the global asymptotic stability of the system
  3. Proof of global classical solutions in n≥1n \geq 1 dimensional spaces with uniform boundedness
  4. Establishment of a complete stability analysis framework including asymptotic behavior of coexistence, extinction, and semi-coexistence states
  5. Provision of numerical verification validating theoretical results in two and three-dimensional spaces

Detailed Methodology

Mathematical Model

The study considers the following coupled nonlinear parabolic PDE system:

u_t = d_1\Delta u - \chi_1\nabla \cdot (u\nabla(vw)) + \mu_1u(1-u+a_1v+a_2w) \\ v_t = d_2\Delta v - \chi_2\nabla \cdot (v\nabla z) + \mu_2v(1-v-a_3u+a_4w) \\ w_t = d_3\Delta w + \xi\nabla \cdot (w\nabla z) + \mu_3w(1-w-a_5u-a_6v) \\ z_t = d_4\Delta z + \alpha v + \beta w - \gamma z \end{cases}$$ where: - $u(x,t)$: secondary predator density - $v(x,t)$: primary predator density - $w(x,t)$: prey density - $z(x,t)$: chemical attractant concentration ### Key Parameter Meanings - $\chi_1$: alarm-taxis coefficient (secondary predator movement toward primary predator-prey density gradient) - $\chi_2$: chemotaxis coefficient (primary predator movement toward chemical signal) - $\xi$: chemical repulsion coefficient (prey movement away from chemical signal) - $d_i$: diffusion coefficients - $\mu_i$: logistic growth coefficients - $a_i$: interspecific interaction coefficients ### Technical Innovations #### 1. Global Boundedness Proof Strategy A multi-step approach is employed: - **Step One**: Establish $L^p(\Omega)$ norm estimates (Lemma 3.1) - **Step Two**: Utilize Neumann heat semigroup techniques to obtain $L^{\infty}(\Omega)$ boundedness (Lemma 3.2) - **Step Three**: Establish gradient estimates $W^{1,\infty}(\Omega)$ (Lemmas 3.4-3.6) #### 2. Lyapunov Functional Construction Corresponding energy functionals are constructed for different equilibrium states: **Coexistence State**: $$E_1(t) = \int_\Omega \left(u-u^*-u^*\ln\frac{u}{u^*}\right) + \Gamma_1\int_\Omega \left(v-v^*-v^*\ln\frac{v}{v^*}\right) + \Gamma_2\int_\Omega \left(w-w^*-w^*\ln\frac{w}{w^*}\right) + \frac{1}{2}\int_\Omega (z-z^*)^2$$ where $\Gamma_1 = \frac{\mu_1a_1}{\mu_2a_3}$, $\Gamma_2 = \frac{\mu_1a_2}{\mu_3a_5}$ ## Main Theoretical Results ### Theorem 1.1 (Global Existence) There exists $\mu > 0$ such that when $\min\{\mu_2,\mu_3\} > \mu$, for non-negative initial data satisfying the conditions, system (1.1) admits a unique classical solution $(u,v,w,z)$ that is uniformly bounded: $$\|u(\cdot,t)\|_{L^{\infty}(\Omega)} + \|v(\cdot,t)\|_{W^{1,\infty}(\Omega)} + \|w(\cdot,t)\|_{W^{1,\infty}(\Omega)} + \|z(\cdot,t)\|_{W^{1,q}(\Omega)} \leq C$$ ### Stability Results The study establishes four classes of asymptotic behaviors: 1. **Coexistence State** (Theorem 1.2): When parameters satisfy specific conditions, solutions converge exponentially to the coexistence equilibrium $(u^*,v^*,w^*,z^*)$ 2. **Secondary Predator Alone State** (Theorem 1.3): Solutions converge to $(1,0,0,0)$ 3. **Semi-Coexistence States** (Theorem 1.4): - Prey extinction state: $(ū,v̄,0,z̄)$ - Primary predator extinction state: $(û,0,ŵ,ẑ)$ ## Numerical Experiments ### Experimental Setup - **Spatial Domain**: $[-0.5,0.5]^n$ ($n=2,3$) - **Grid**: $101×101$ for 2D, $50×50×50$ for 3D - **Initial Conditions**: Radially symmetric bell-shaped functions - **Parameters**: $d_i=\chi_1=\chi_2=\xi=\mu_i=\alpha=\beta=1$, $\gamma=2$ ### Numerical Results Verification Four typical examples verify the theoretical predictions: 1. **Example 5.1**: With $a_i=0.5$, solutions converge to coexistence state $(1.286877, 0.428973, 0.144368, 0.286793)$ 2. **Example 5.2**: Under specific parameters, solutions converge to $(1,0,0,0)$ 3. **Examples 5.3-5.4**: Verify convergence of two semi-coexistence states ## Technical Challenges and Solutions ### Main Technical Difficulties 1. **Nonlinear Term Treatment**: Complexity of the alarm-taxis term $\chi_1\nabla \cdot (u\nabla(vw))$ 2. **Coupled System Analysis**: Strong coupling nature of the four equations 3. **Gradient Estimates**: Technical complexity in high-dimensional spaces ### Key Technical Tools 1. **Gagliardo-Nirenberg Inequalities** 2. **Neumann Heat Semigroup Theory** 3. **Maximal Sobolev Regularity** 4. **Young and Hölder Inequalities** ## Comparison with Related Work ### Historical Development - **Keller-Segel Model** (1970): Foundational work on chemotaxis - **Haskell-Bell Model** (2021): First mathematical modeling of alarm-taxis - **Existing Extensions**: Primarily focused on single mechanisms or low-dimensional cases ### Innovations of This Paper Compared to existing work, this paper is the first to: - Rigorously handle the simultaneous action of both taxis mechanisms - Establish global theory in high-dimensional spaces - Provide a complete stability analysis framework ## Conclusions and Discussion ### Main Conclusions 1. Under appropriate parameter conditions, the chemo-alarm-taxis system admits global bounded classical solutions 2. The long-term behavior of the system is determined by parameter relationships, potentially converging to four different equilibrium states 3. Chemotaxis and alarm-taxis coefficients play crucial roles in determining system stability ### Biological Significance - Validates the mathematical feasibility of the "thief alarm" hypothesis - Explains complex three-species interaction dynamics in nature - Provides theoretical foundation for ecological conservation strategies ### Limitations 1. **Parameter Conditions**: Stability results require relatively strong parameter constraints 2. **Linearization Analysis**: Local stability analysis is based on linearization 3. **Numerical Verification**: Verification only conducted within specific parameter ranges ### Future Directions 1. Relax restrictions on parameter conditions 2. Consider more complex reaction terms 3. Investigate spatial pattern formation 4. Extend to stochastic environments ## In-Depth Evaluation ### Strengths 1. **Theoretical Rigor**: Provides a complete mathematical proof framework 2. **Methodological Innovation**: Cleverly combines multiple analytical techniques 3. **Result Completeness**: Encompasses existence, boundedness, and stability 4. **Numerical Verification**: 2D and 3D numerical experiments support theoretical results ### Weaknesses 1. **Complex Parameter Conditions**: Parameter verification is difficult in practical applications 2. **Biological Realism**: Some mathematical assumptions may be overly idealized 3. **Computational Complexity**: Numerical computation challenges in high-dimensional cases ### Impact - **Theoretical Contribution**: Provides new analytical framework for taxis system research - **Application Value**: Provides important tools for ecological mathematical modeling - **Methodological Value**: Lyapunov functional construction method has generalization potential ### Applicable Scenarios - Marine ecosystem modeling - Signal transduction in cellular biology - Immune system dynamics analysis - Ecological conservation strategy formulation ## References The paper cites 46 related references, covering important works in multiple fields including chemotaxis theory, alarm-taxis modeling, and partial differential equation theory, providing a solid theoretical foundation for the research. --- **Overall Assessment**: This is a high-quality paper of significant theoretical importance in mathematical biology, establishing for the first time a complete mathematical theory for chemo-alarm-taxis systems. The methodology is rigorous, the results are profound, and it has important implications for the development of related fields.