Dynamics in a predator-prey model with memory e ﬀ ect in predator and fear e ﬀ ect in prey

: The spatial memory e ﬀ ect in predator and fear e ﬀ ect in prey are incorporated in a di ﬀ usive predator-prey model. We are interested in studying the dynamics generated by the memory e ﬀ ect and fear e ﬀ ect, and mainly study the local stability of coexisting equilibrium, the existence of Hopf bifurcation and the property of Hopf bifurcation. Through the numerical simulations, we show that increasing memory-based di ﬀ usion coe ﬃ cient is not conducive to the stability of the coexisting equilibrium, and the fear e ﬀ ect has both stabilizing and destabilizing e ﬀ ect on the coexisting equilibrium under di ﬀ erent parameters.


Introduction
The relationship between predator and prey is an important research content in ecosystem, and many scholars have studied this interaction by differential equation models [1][2][3][4][5][6]. The direct relationship between predator and prey is the consumption of prey. One of the indirect relationship between predator and prey is the fear of predator. When the prey gets the predator's signal (chemical/vocal), they will increase their vigilance time and reduce their foraging [7], such as mule deer v.s. mountain lions [8], elk v.s. wolves [9]. Consider the fear effect, Panday et al. [10] proposed the following model (1.1) All parameters are positive. The biological interpretation of parameters is given in Table 1. They incorporated fear effect by modifying the prey intrinsic growth rate R 0 as R 0 1+Kv . By the scaling Panday et al. [10] introduced time delay of perceiving predator signals to model (1.3), and mainly studied the boundedness, persistence, local and global behavior of the delayed model. In the real world, in addition to the fear effect of the prey, the clever predator also has spatial-memory and cognition [11], which is often ignored in modeling the predator-prey interaction. For example, blue whales rely on memory for migration, which is presented by B. Abrahms et al. [12] and W. F. Fagan [13]. As another example, animals in polar regions usually determine their spatial movement by judging footprints, which record the history of species distribution and movement, involving time delay [14]. Obviously, highly developed animals can even remember the historical distribution or cluster of species in space. Much progress has been made in implicitly integrating spatial cognition or memory [14][15][16][17][18]. To incorporate the memory effect, Shi et al. proposed a single specie model with spatial memory by introducing a additional delayed diffusion term [14]. They supposed that in addition to the negative gradient of the density distribution function at the present time, there is a directed movement toward the negative gradient of the density distribution function at past time [14]. After this pioneering work, some recent works [19][20][21][22][23] about the population model with memory effect have emerged. In [23], Song et al. obtained a computing method for the normal forms of the Hopf bifurcations in the diffusive predator-prey model with memory effect, which is friendly to use. Inspired by the above work, we suppose the predator has spatial-memory diffusion and the prey has fear effect, then modified the model (1.3) as follow All parameters are positive. The biological description of parameters is given in Table 1. The term −d∇(v∇u(t − τ)) represents the memory-based diffusion effect of the predator. The Neumann boundary conditions is used. The aim of this paper is to study the effect of predator's memory-based diffusion and prey's fear on the model (1.4).
The rest of this paper is organized as follows. In Section 2, the stability of coexisting equilibrium and the existence of Hopf bifurcation are considered. In Section 3, the property of Hopf bifurcation is studied. In Section 3, some numerical simulations are given to analyze the effect of spatial-memory and fear effect. In Section 4, a short conclusion is given.

Stability analysis
For simplicity, we choose Ω = (0, lπ). Denote N as positive integer set, and N 0 as nonnegative integer set. It is easy to obtain (0, 0) and (1, 0) are two boundary equilibria of model (1.4). Next, we will give the existence of coexisting equilibrium.
Proof. The coexisting equilibrium of (1.4) is a positive root of the following equations From the second equation, we have u = µ a−bµ . Substitute it into the first equation, we have Then v = 1 2c −1 ± 1 + 4c (a−(1+b)µ) (a−bµ) 2 . Obviously, the conclusion holds.
In this paper, we mainly study the stability of coexisting equilibrium E * (u * , v * ). Linearize model (1.4) where where When τ = 0, the characteristic Eq (2.3) become l 4 . Make the following hypothesis is locally asymptotically stable under the hypothesis (H 1 ).
Proof. If (H 1 ) holds, we can easily obtain that A n > 0 and B n + C n > 0. Then the characteristic roots of (2.4) all have negative real parts. Then E * (u * , v * ) is locally asymptotically stable.

τ > 0
In the following, we assume (H 1 ) holds. Let iω (ω > 0) be a solution of Eq (2.3), then we have We can obtain cosωτ = ω 2 −B n C n , sinωτ = A n ω C n > 0 under hypothesis (H 1 ). It leads to Let p = ω 2 , then (2.5) becomes and the roots of (2.6) are (2.7) The existence of purely imaginary roots of Eq (2.3) can be divided into the following two cases.
3) has a pair of purely imaginary roots ±iω + n at τ j,+ n for j ∈ N 0 and n ∈ M. Otherwise, Eq (2.3) does not have characteristic roots with zero real parts.
Case 2 : a 2 1 + 2a 2 b 1 < 0. This case can be divided into the following two subcases.
3) has two pairs of purely imaginary roots ±iω ± n at τ j,± n for j ∈ N 0 and n ∈ M 1 . Otherwise, Eq (2.3) does not have characteristic roots with zero real parts.
3) has two pairs of purely imaginary roots ±iω ± n at τ j,± n for j ∈ N 0 and n ∈ M 1 . For d > d * and n ∈ M, Eq (2.3) has a pair of purely imaginary roots ±iω + n at τ j,+ n for j ∈ N 0 and n ∈ M. Otherwise, Eq (2.3) does not have characteristic roots with zero real parts.
The ω ± n and τ j,± n are defined as follow We have the following lemma.
Theorem 2.2. Assume (H 1 ) holds, then the following statements are true for model (1.4).
Remark 2.1. In the Theorem 2.2, if E * (u * , v * ) is locally asymptotically stable, the densities of prey and predator will tend to the equilibrium state in the whole region when the initial densities of prey and predator is near E * (u * , v * ). When Hopf bifurcation occurs at (u * , v * ), then the densities of prey and predator will produce periodic oscillation. Especially, spatially homogeneous periodic oscillations may occur when τ near the critical value τ = τ j,+ 0 or τ = τ j,− 0 , and spatially inhomogeneous periodic oscillations may occur when τ near the critical value τ = τ j,+ n or τ = τ j,− n (n > 0).

Property of Hopf bifurcation
In this section, we use the algorithm in [23] to compute the normal form of Hopf bifurcation. We denote the critical value of Hopf bifurcation asτ and the purely imaginary roots as ±iω n of Eq (2.3).

Numerical simulations
In this section, we give some numerical simulations to analyze the effect of spatial memory in predator and fear in prey on the model (1.4). Fix the following parameters  Figure 1. By the Theorem 2.2, we know that (u * , v * ) is locally stable for τ ≥ 0 when d < d * . But when d > d * , the inhomogeneous Hopf bifurcation curves exist. This means that increasing parameter d is not conducive to the stability of the equilibrium (u * , v * ), and the densities of prey and predator will produce spatially inhomogeneous periodic oscillation. 1380. When τ ∈ [0, τ * ), (u * , v * ) is locally stable (Figure 2). By direct calculation, we can obtain K 1 ≈ 0.0166, and K 2 ≈ −0.0699. Then, the Hopf bifurcation is supercritical and the bifurcating periodic solution is stable (Figure 3). At this time, the bifurcating periodic solution is spatially inhomogeneous and with mode-2.

The effect of c
Next, we will study the effect of fear effect c on the model (1.4). Fix the parameters as (4.1), then model (1.4) has a unique coexisting equilibrium (u * , v * ). And a 2 1 + 2a 2 b 1 > 0 when 0 < c < 2.6194. We give the figure of d * with parameter c as in Figure 7. Set parameter d = 0.7 and d = 0.8, we give the bifurcation diagrams of model (1.4) with parameter c as in Figure 8.
When d = 0.7 and τ = 20, increasing parameter c can destroy the stability of the coexisting equilibrium (u * , v * ), and induce spatially inhomogeneous periodic solution ( Figure 9). This means that increasing parameter c is not conducive to the stability of the coexisting equilibrium (u * , v * ).
When d = 0.8 and τ = 9, increasing parameter c can destroy the stability of the coexisting equilibrium (u * , v * ), and induce spatially inhomogeneous periodic solution initially ( Figure 10). But when c is larger enough, increasing parameter c can rule out the spatially inhomogeneous periodic oscillation and stabilize the coexisting equilibrium (u * , v * ) ( Figure 10). This means that increasing parameter c is not conducive to the stability of the equilibrium (u * , v * ), initially. But when c is large, increasing parameter c is conducive to the stability of the coexisting equilibrium (u * , v * ).

Conclusions
In this paper, we incorporate the memory effect in predator and fear effect in prey into a predatorprey model. By using time delay in the memory of predator as bifurcating parameter, we analyze the local stability of coexisting equilibrium, the existence of Hopf bifurcation, and the property of Hopf bifurcation by the method in [23]. Through the numerical simulations, we analyzed the effect of memory effect in predator and fear in prey on the model.
The spatial memory effect plays an important role in the dynamics of the predator-prey model. Through the numerical simulations, we observed that the memory-based diffusion coefficient d has destabilizing effect on the predator-prey model when it is larger than some critical value. In addition. when d crosses the critical value, time delay τ in the memory of predator can affect the stability of the equilibrium (u * , v * ). In the numerical simulations, we observe that the first Hopf bifurcation curve is inhomogeneous bifurcation curve, and homogeneous Hopf bifurcation curve does not exist. This is different from the predator-prey model without the spatial memory effect. When τ crosses the critical value τ * , the densities of prey and predator will produce spatially inhomogeneous periodic oscillation. When τ crosses the second critical value, the spatially inhomogeneous periodic oscillations with different modes exist, but the densities of prey and predator will converge to the spatially inhomogeneous periodic solution corresponding to the first bifurcation curve. This shows that the spatially memory effect in predator can destroy the stability of the coexisting equilibrium, and induce spatially inhomogeneous periodic oscillations.
In addition, the fear effect parameter c in prey can also affect the stability of the coexisting equilibrium (u * , v * ). A small fear effect parameter c means a large birth rate 1 1+cv , then the large birth rate can support fluctuations. Increasing parameter c can destroy the stability of the coexisting equilibrium (u * , v * ), and induce spatially inhomogeneous periodic solution. Hence, we observed the destabilizing effect on the the coexisting equilibrium (u * , v * ). A large fear effect parameter c means a low birth rate, then the low birth rate can not support fluctuations. Increasing parameter c can rule out the spatially inhomogeneous periodic oscillation and stabilize the coexisting equilibrium (u * , v * ). Hence, we observed the stabilizing effect on the the coexisting equilibrium (u * , v * ).