Local and global dynamics of a prey–predator system with fear, Allee effect, and variable attack rate

Fear prompts prey to adopt risk-averse behaviors, such as reduced foraging activity, increased vigilance, and avoidance of areas with high predator presence, which affects its reproduction. In a real scenario, a population requires a minimum density to avoid extinction, known as an Allee threshold. In light of these biological factors, we propose a predator–prey model with (i) a fear effect in a prey population, (ii) an Allee effect in a predator population, and (iii) a non-constant attack rate that modifies the functional response. We ensured the non-negativity and boundedness of the solutions and examined the local and global stability status for each existing steady state solutions. We investigated some deep dynamical properties of the system by varying different parameters, such as cost of fear in prey and strength of the Allee effect in predators and their mortality rate. In codimension one bifurcations, we observed saddle node, Hopf, homoclinic, and coalescence of two limit cycles. Additionally, codimension two bifurcations were observed, including Bautin and Bogdanov Takens bifurcations. To provide a clearer understanding of these bifurcations, we conducted biparametric analysis involving the fear and Allee parameters, as well as the fear parameter and predator mortality rate. Our investigation shows that cost of fear and strength of Allee strongly influences the survival status of the predator. Furthermore, bistability and tristability reveal that the survival and extinction of predator are dependent on the initial population level. Numerical simulations and graphical illustrations are provided to support and validate our theoretical findings.