Strong Allee effect and basins of attraction in a discrete‐time zoonotic infectious disease model

Motivated by the Feline immunodeficiency virus, the virus that causes AIDS in cat populations, we use discrete‐time infectious disease models with demographic strong Allee effect to examine the impact of the fatal susceptible‐infected (SI) infections on two different types of growth functions: Holling type III or modified Beverton–Holt per‐capita growth function (compensatory density dependence), and Ricker per‐capita growth function with mating (overcompensatory density dependence). The occurrence of the strong Allee effect in the disease‐free equation renders the SI population model bistable, where the two coexisting locally asymptotically stable equilibrium points are either the origin (catastrophic extinction state) and the second fixed point (compensatory dynamics) or the origin and an intrinsically generated demographic period k > 1 population cycle (overcompensatory dynamics). We use the basic reproduction number, ℛ 0 , and the spectral radius, λ k , to examine the structures of the coexisting attractors. In particular, we use MATLAB simulations to show that the fatal disease is not only capable of enlarging or shrinking the basin of attraction of the catastrophic extinction state, but it is also capable of fracturing the basins of attraction into several disjoint sets. Thus, making it difficult to specify the asymptotic zoonotic SI disease outcome in terms of all initial infections. The complexity of the basins of attractions appears to increase with an increase in the period of the intrinsically generated demographic population cycles.