Dynamic complexities in a predator–prey model with prey refuge influenced by double Allee effects

Within the context of a two-dimensional framework encompassing interacting species, an examination is conducted in this study on the double Allee effect and prey refuge, considering both species in the interaction. The stability of the feasible equilibrium of the system and diverse bifurcation patterns including codimension-one and codimension-two bifurcations are scrutinized through theoretical and numerical investigations, which reveals the complex dynamics induced by saturated functional response and double Allee effects. Additionally, one-parameter bifurcation diagrams and two-parameter bifurcation diagrams are constructed to intricately evaluate the system’s dynamics indicative of the presence of multiple attractors like bi-stability and tri-stability. Lastly, the sensitivity analysis is performed to delve into the effect of system parameters on species density, which indicates that the parameter $\eta$ proportional to the conversion rate is the most sensitive parameter. A brief discussion further reveals that the model without double Allee effect reduces dynamic complexity.