Global stability of three trophic levels predator–prey model with alarm-taxis

This paper is concerned with the three trophic levels predator–prey system with alarm-taxis

$$\begin{aligned} \left\{ \begin{array}{lll} u_{t}=d_{1} \Delta u+u\left( 1-u-\frac{a v}{v+\rho }\right) , &{} x \in \Omega , &{} t>0, \\ v_{t}=d_{2} \Delta v+v\left( \frac{b u}{v+\rho }-\alpha -\frac{c w}{w+\sigma }\right) , &{} x \in \Omega , &{} t>0, \\ w_{t}=d_{3} \Delta w-\chi \nabla \cdot \left( w\nabla (uv)\right) +w\left( \frac{m v}{w+\sigma }-\beta \right) , &{} x \in \Omega , &{} t>0 \end{array}\right. \end{aligned}$$

under homogeneous Neumann boundary condition in smooth bounded domains \(\Omega \subset {\mathbb {R}}^n (n\ge 1)\). We prove that the system possesses a unique global bounded classical solution for all sufficiently smooth initial data. Moreover, we show the large time behavior of the solution with convergence rates and perform some numerical simulations to verify the analytic results.