Hopf bifurcation analysis in a delayed di ﬀ usive predator-prey system with nonlocal competition and generalist predator

: A delayed di ﬀ usive predator-prey system with nonlocal competition and generalist predators is considered. The local stability of the positive equilibrium and Hopf bifurcation at positive equilibrium is studied by using time delay as a parameter. In addition, the property of Hopf bifurcation is analyzed using the center manifold theorem and normal form method. It is determined that time delays can a ﬀ ect the stability of the positive equilibrium and induce spatial inhomogeneous periodic oscillation of prey and predator population densities.


Introduction
Dynamic systems are widely used in nature, such as in infectious disease models [1,2], vegetation-water models [3], and population models [4,5]. For example, Khan et al. [6] studied the dynamics of pine wilt disease with variable population sizes and showed control strategies for the possible elimination of the infection in the pine tree population. Khan et al. [7] obtained an efficient iterated homotopy perturbation transform method (IHPTM) to solve a mathematical model of HIV infection of CD4+ T cells. Khan et al. [8] proposed a new method to solve partial differential equations arising in the fields of science and engineering. Among the applications of dynamical systems, the predator-prey model is an important research topic [9][10][11][12]. Generally, predator-prey models assume that the prey is the only food source of the predator [13][14][15].
In ecosystems, predators are usually generalist predators. They feed on many types of species and can change their diet to another species when a focal prey population begins to run short [16][17][18][19].
In [17], the authors studied a predator-prey model with a generalist predator. They aimed to achieve biological control by generalists [17]. In [18], the authors studied the spatial pattern formation of a predator-prey model with generalist predator and the harvesting, refuge effect. To better understand the relationship between prey and generalist predators, Upadhyay and Agrawal [19] modified the Leslie-Gower model as dv to describe the growth law of generalist predators. The term dv 2 represents the growth of predators, where d is the mating rate of the predator. e is another source of food for predators. They considered the following model: (1.1) u(t) and v(t) represent the prey and predator densities, respectively. The term ωu A+Bu+v represents is the Beddington-DeAngelis functional response. τ is the gestation delay of the predator. Upadhyay and Agrawal studied the invariance, boundedness, and local and global stability and Hopf bifurcation.
In the real world, the spatial distribution of populations is often inhomogeneous, so it is more practical to use a reaction-diffusion model to describe the relationship between predators and prey. In addition, the reaction-diffusion predator-prey model shows more abundant dynamic properties such as spatial patterns and spatial nonhomogeneous periodic solutions. Considering this factor, Liu et al. [20] proposed the following model with the Crowley-Martin type functional response based on the model (1.1).
In the model (1.2), Liu et al. [20] used the Crowley-Martin-type functional response to reflect the impact of predators on prey. This considered the effect of interference among predators. a, b, and c represent the capture rate, handling time, and magnitude of interference among predators, respectively. They mainly analyzed the instability and Hopf bifurcation induced by time delay [20]. Although they pointed out that time delay may cause spatially inhomogeneous periodic solutions, the numerical simulations did not show stably spatially inhomogeneous periodic solutions. This is because in the delayed reaction-diffusion predator-prey model, the spatial nonhomogeneous Hopf bifurcation curve is usually above the spatial homogeneous Hopf bifurcation curve. In addition, Turing instability cannot occur for the model (1.2). For model (1.2), it is unfortunate that the inhomogeneous distribution of prey and predator in space is not shown. This may be due to the lack of stable spatially inhomogeneous periodic solutions in the delayed reaction-diffusion predator-prey model. In nature, competition within populations exists widely, and this competition is often nonlocal since resources are limited. In [21,22], the authors suggested that the internal competition of the population caused by the natural environment is related not only to the population density at the current location but also to the population density nearby. They measured this effect by weighting and integrating, and modified the u K as 1 K Ω G(x, y)u(y, t)dy. G(x, y) is a kernel function. Wu and Song studied Hopf and state-Hopf bifurcations in a diffusive predator-prey model with a nonlocal effect and delay [23]. Geng et al. studied a diffusive predator-prey model with nonlocal competition, including Hopf, Turing, double-Hopf, and Turing-Hopf bifurcations [24]. The works in [25,26] show that stable spatially inhomogeneous periodic solutions often exist in predator-prey models with nonlocal competition.
Inspired by the above work, we assume there is nonlocal competition in prey and modify the model (1.3) as follows.
where d 1 and d 2 represent the diffusion coefficients of prey and predator, respectively. The term au (1+bu)(1+cv) is the Crowley-Martin-type functional response. Ω is the prey and predator's living region.
Just for the convenience of calculation, we choose Ω = (0, lπ), where l > 0. The boundary condition is the Newman boundary condition, which means that the living region is closed and no prey and predator enter or leave the region. Ω G(x, y)u(y, t)dy represents the nonlocal competition effect, and G(x, y) is the kernel function. We assume that the competition strength among prey individuals in the habitat is the same; we choose G(x, y) = 1 lπ as in previous works [23][24][25][26]. To our knowledge, there is no work about the predator-prey model (1.3) with nonlocal competition in prey and generalist predators from the point of Hopf bifurcation. The aim of this paper is to study the combined effect of time delay and nonlocal competition on model (1.3). Compared with the model (1.2), do new dynamics appear, such as stably spatially inhomogeneous periodic solutions and Turing instability?
This paper is organized as follows: the stability and existence of Hopf bifurcation are studied in Section 2, the property of Hopf bifurcation is analyzed in Section 3, and numerical simulations are given in Section 4. Finally, a short conclusion is presented.

Nondelay model
When τ = 0, the characteristic equation (2.3) is (2.6) Make the following hypothesis: Proof. If (H 1 ) holds, then we can determine that the characteristic roots of (2.5) all have negative real parts. Then, E * (u * , v * ) is locally asymptotically stable. If (H 1 ) holds, then the characteristic roots of (2.5) with k ∈ N have at least one positive real part, but with n = 0 all having a negative real part. This implies that E * (u * , v * ) is Turing unstable. Lemma 2.3. When τ = 0, the Turing instability of E * (u * , v * ) cannot occur for the system (1.2), which lacks the nonlocal competition term.
Proof. When τ = 0, for the system (1.3) without nonlocal competition, the characteristic equation is locally asymptotically stable in the absence of spatial diffusion, which means that a 1 +â < 0. However, T n > 0 and D n > 0 when a 1 +â < 0 holds. Then, the Turing instability of E * (u * , v * ) cannot occur.

Delay model
We can obtain It leads to Let z = ω 2 ; then, (2.8) becomes and the roots of (2.9) are where P n = A 2 n − 2B n − b 2 1 , Q n = B n + C n , and R n = B n − C n . If (H 0 ) and (H 1 ) hold, Q n > 0 (n ∈ N 0 ). By direct calculation, we have (2.10) We have the following lemma.

Numerical simulations
To verify our theoretical results, we give the following numerical simulations. Fix parameters The bifurcation diagram of system (1.3) with the parameter of interference magnitude among predators c is given in Figure 1, where (H 0 ) and (H 1 ) hold. We can see that with an increase in parameter c, the Hopf bifurcation curves decrease, which implies that the stable region of positive equilibrium (u * , v * ) will decrease. This means that increasing the interference magnitude among predators is not conducive to the homogeneous distribution of prey and predators in space. This causes spatial oscillations in prey and predator densities. Since the inhomogeneous Hopf bifurcation curve τ 0 1 is always under the homogeneous Hopf bifurcation curve τ 0 0 . This spatial oscillation of prey and predator density is inhomogeneous. Fixing c = 0.1, we can determine that (u * , v * ) ≈ (9.0000, 1.9048) is the unique positive equilibrium, and (H 1 ) holds. By direct calculation, we have τ * = τ 0 1 ≈ 0.0758 < τ 0 0 ≈ 0.1002, µ 2 ≈ 0.0105, β 2 ≈ −0.0584 and T 2 ≈ 0.8602. If we choose τ = 0.05 < τ * , then (u * , v * ) is locally asymptotically stable for model (1.3) (shown in Figure 2) and model (1.2) (shown in Figure 3). This means that the nonlocal competition term has no effect on the stability of the model (1.3). Prey and predator will coexist in a spatially homogeneous form, and their densities will converge to the positive equilibrium (u * , v * ).  If we choose τ * < τ = 0.09 < τ 0 0 , then (u * , v * ) is unstable, and inhomogeneous periodic solutions exist for the model (1.3) (shown in Figure 4). This means that prey and predator will coexist in the form of spatially inhomogeneous oscillations, and their densities will distribute inhomogeneously in space. To compare our result with the work in [20], we give the numerical simulations of model (1.2) under the same parameter τ = 0.09 in Figure 5. We can see that (u * , v * ) is locally asymptotically stable (1.2). This means that prey and predators will still coexist in a spatially homogeneous form and that their densities will converge to positive equilibrium (u * , v * ).

Conclusions and discussion
Assuming there is nonlocal competition caused by limited resources in prey, we modified the model (1.2) to the model (1.3) with the nonlocal competition in prey, gestation delay in predator, and generalist predator. We study the local stability of the positive equilibrium and Hopf bifurcation at the positive equilibrium by using the time delay as a parameter. Through the center manifold and normal form method, we give some parameters that can determine the direction of Hopf bifurcation, the stability of bifurcating periodic solutions, and the period of periodic solutions.
For the time delay when τ = 0, we determine that nonlocal competition can induce Turing instability of the coexisting equilibrium point (u * , v * ) , but Turing instability cannot occur in the model (1.2) without nonlocal competition. Similar to the results in [20], time delays can affect the stability of the positive equilibrium. The positive equilibrium is stable when the delay is shorter than the critical value but unstable when the delay is longer than the critical value, and prey and predator densities will produce periodic oscillations.
Through the bifurcation diagram, we find that increasing the interference magnitude among predators is not conducive to the homogeneous distribution of prey and predators in space. This causes spatial oscillations in prey and predator densities. Compared with the work in [20], the numerical simulations show that the inhomogeneous Hopf bifurcation curve τ 0 1 is always under the homogeneous Hopf bifurcation curve τ 0 0 for the model (1.3) with nonlocal competition. When the time delay is longer than the critical value τ 0 1 and shorter than τ 0 0 , the prey and predator will coexist in space, but their densities will oscillate periodically in the form of spatially inhomogeneous oscillations for the model (1.3). However, the prey and predator will still coexist in a spatially homogeneous form, and their densities will converge to the positive equilibrium for the model (1.2). The model (1.3) is more realistic than the model (1.2) in the real world since it is difficult for the population density to exist in a completely uniform way. Through our theoretical analysis and numerical simulation, we show that stable inhomogeneous periodic solutions exist when the delay crosses the critical value, which is different from the work in [20].
In the inhomogeneous periodic oscillations of prey and predator densities (Figure 4), we find that the amplitude near the boundary of the region is larger than that at the center of the region. This is because the prey needs to avoid the hunting of predators. They will escape from the interior of the area and move out to the boundary. However, because the living area is closed, they cannot cross the boundary. Thus, they can only move inward again. After reaching the interior of the area, they have to escape the hunting of predators again. This forms the periodic oscillation of this mode.
In applying the predator-prey model, we suggest that scholars consider the nonlocal competition of prey since the dynamics of the models are different in the absence and presence of nonlocal competition. We hope that our model can be applied to ecological environment protection and population control. In future work, we will study the predator-prey model with a nonlocal competition term, nonconstant kernel function, and the Allee effect.

Y. Shao,
Dynamics of an impulsive stochastic predator-prey system with the Beddington-DeAngelis functional response, Axioms,