The interaction with predator-prey populations and their possible outcomes are probably the most studied topics in ecology, because of their existence and universal relevance, will continue to be one of the dominant topics in both ecology and mathematical ecology . In spite of the predator-prey theory has undergone significant developments in the past several years, many long-standing mathematical and ecological problems still deserve researchers’ attention    . The functional response of predators to their prey density refers to the change in the density of attached prey per unit time per predator with the change in prey density. It is the average number of prey killed per individual predator per unit of time. For the functional response functions, there are many types. The most commonly used functional response function, which was suggested by Holling (1959) and called Michaelis-Menten type or Holling type II functional re-
sponse . And it is takes the form: . where u represents the density of the prey population, and the positive constants m (units: 1/times) and a (units: 1/prey) describe the effects of capture rate and handling time, respectively, on the feeding rate.
However, Prey-dependent functional responses cannot describe the interaction between predators, and have been challenged by the biology and physiology communities  . Some biologists have argued that in many situations, especially when predators have to search for food (and therefore, have to share or compete for food), the functional response in a prey-predator model should be predator-dependent. There is significant evidence that predator dependence in functional response occurs frequently in laboratories and natural systems  . Given that many experiments and observations show predators do indeed interfere with each other’s activities to trigger competition effects and that prey changes its behavior under increased predator threat, models with a functional predator-dependent response are reasonable alternatives to models with prey-dependent functional response.
Starting with this argument and the traditional prey-only model, to describe the mutual interference of a predator, Beddington (1975) and DeAngelis (1975)     proposed that an individual from a population of over two predators not only allocate time in searching for and processing their prey but also takes time in encountering with other predator. This result in the so-called
Beddington-DeAngelis functional response . The Beddington-DeAngelis is like the well-known Holling type II functional response, but
it has an additional term cv in the denominator modelling mutual interference among predators and has some of the same qualitative features as the ratio-dependent form but avoids some of the singular behaviors of ratio-dependent models at low densities which have been the source of controversy.
The Beddington-DeAngelis functional response and the analogous response with diffusion in a constant environment have received much attention in the literatures  . Such as Crowley and Martin  assumed that interference among predators occurs regardless of whether a particular predator is looking for prey or handling prey and proposed a functional response, which is called the Crowley-Martin-type functional response, and its takes the form:
It is assumed that the predator feeding rate decreases by higher predator density even when prey density is high, and therefore the effects of predator interference in the feeding rate remain important all the time whether an individual predator is handling or searching for a prey at a time. And in . The authors consider a Leslie-Gower predator-prey model with a functional Crowley-Martin response describing predator mutual interference, which its takes the form:
Initially, their analysis focuses on local asymptotic stability and Hopf bifurcation of the spatially homogeneous model based on ODE system for the reaction-diffusion model.
Based on the above discussion, and therefore. In this paper, we consider the following predator-prey model with mutual interference with Beddington-De Angelis functional response:
where and represent population densities of prey and predator at time t, respectively. r, K, s, m, a, b, and h are positive constants, K is the carrying capacity of the prey, and r and s denote to their intrinsic growth rate, respectively,
is a function on the prey population size (h is a measure of the food quality
of the prey for conversion into predator growth depending on the density of the population).
For the model (1.2), we mainly discuss about Hopf bifurcation on the positive equilibrium for the ODE system.
The aim of this article is to show that the diffusive predator-prey model (1.2) shows different spatial, temporal and spatiotemporal patterns across the mechanisms described above.
For simplicity, applying the following scaling:
System (1.2) can be simplified by:
where is a bounded smooth domain and is the outward unit normal vector on . The constants and are the diffusion coefficients of and , which represent the natural dispersive force of movement of the prey and predator density, respectively. The condition of the homogeneous Neumann boundary means that the two species have zero-flux across the boundary . The initial conditions , are smooth functions satisfying .
In section two of this paper, we investigate the asymptotical behaviour of the interior equilibrium and occurrence of Hopf bifurcation of the local ODE system of (1.3).
2. Hopf Bifurcation and Stability Analysis for ODE System
Firstly, we consider the stability of the following ODE system:
When , we can see that the system has a boundary equilibrium , two positive interior equilibrium
if and only if
where and .
By the simple calculations, we note that the boundary equilibrium is unstable, we will study the properties of the positive interior equilibrium. Thus, consider the stability of . Note that the Jacobian of (1.3) at is
For convenience, we denote that and , so if and only if
Thus, we can obtain the following theorem.
Theorem 2.1. Suppose that (H1).
1-1) if , then is locally asymptotically stable.
1-2) if (H2) and hold, then is unstable and there exists at least one periodic orbit for system (1.3).
2) is unstable.
1) The characteristic equation corresponding to is
Therefore, yield that and combined with , we completed the (1-1). Otherwise, if , then is unstable.
2) Apply the same argument to the , we have
thus, is unstable.
In the following, we analyse the existence of Hopf bifurcation at the interior equilibrium , choose s as the bifurcation parameter.
Denote by the two roots of the characteristic Equation (2.1) with and . Obviously, (2.1) has a pair of imaginary roots at . Moreover, . This shows that the transversality condition holds. Thus undergoes Hopf bifurcation of as s passes through the .
Now, to make a better use of Hopf bifurcation theory, We translate the interior equilibrium to the origin by the transformation , . For the sake of convenience, we still denote and by u and v, respectively. Thus, the system (1.3) is transformed into
Then, rewrite the system (2.2) as follows:
By setting the matrix
where and . It is easy to obtain that
When , we have
By the transformation , system (2.3) becomes
Rewrite the system (2.5) in the following polar coordinate form:
And then, from the Taylor expansion of (2.6) at we have:
Now, to determine the stability of the Hopf bifurcation periodic solution, So, we need to calculate the sign of the coefficient , which is given by the following formula:
where all partial derivatives are evaluated at the bifurcation point , and
Noting that , thus we can calculate that
Now by calculation, we can determine the value and sign of in (2.8).
Recall that and , by Poincaré-Andronov-Hopf Bifurcation Theorem, we have the following result.
Assume that (H1) and (H2) hold. The system (1.3) undergoes a Hopf bifurcation at the coexistence equilibrium when . Furthermore:
1) determines the stability of the bifurcated periodic solutions: if then the bifurcating periodic solutions are stable (unstable).
2) determines the directions of the Hopf bifurcation: if then the Hopf bifurcation is supercritical (subcritical).
In this article, we showed that the predator-prey model with mutual interference exhibits a rich and interesting dynamic behavior. We first studied the local stability and Hopf bifurcation in the corresponding ODE system. Then we studied the existence and direction of Hopf bifurcation and the stability of the bifurcating periodic solution in the reaction-diffusion system. The system (1.3) undergoes a Hopf bifurcation at the interior equilibrium when .
The authors are greatly indebted to the anonymous referees for their careful reading and helpful comments, especially the referee who gave constructive suggestions for revision of the manuscript.
 Kuang, Y. (1995) Global Stability in Delay Differential Systems without Dominating Instantaneous Negative Feedbacks. Journal of Differential Equations, 119, 503-532.
 Skalski, G.T. and Gilliam, J.F. (2001) Functional Responses with Predator Interference: Viable Alternatives to the Holling Type II Model. Ecology, 82, 3083-3092.
 Dolman, P.M. (1995) The Intensity of Interference Varies with Resource Density: Evidence from a Field Study with Snow Buntings, Plectrophenax nivalis. Oecologia, 102, 511-514.
 Jost, C. and Ellner, S.P. (2000) Testing for Predator Dependence in Predator-Prey Dynamics: A Non-Parametric Approach. Proceedings of the Royal Society of London. Series B: Biological Sciences, 267, 1611-1620.
 Antwi-Fordjour, K., Parshad, R.D. and Beauregard, M.A. (2020) Dynamics of a Predator-Prey Model with Generalized Holling Type Functional Response and Mutual Interference. Mathematical Biosciences, 326, Article ID: 108407.
 Hwang, T.-W. (2004) Uniqueness of Limit Cycles of the Predator-Prey System with Beddington-DeAngelis Functional Response. Journal of Mathematical Analysis and Applications, 290, 113-122.
 Cantrell, R.S. and Cosner, C. (2001) Effects of Domain on the Persistence of Populations in a Diffusive Food-Chain Model with Beddington-Deangelis Functional Response. Natural Resource Modeling, 14, 335-367.
 Bai, D.Y., Yu, J.S., Fan, M. and Kang, Y. (2020) Dynamics for a Non-Autonomous Predator-Prey System with Generalist Predator. Journal of Mathematical Analysis and Applications, 485, Article ID: 123820.
 Crowley, P.H. and Martin, E.K. (1989) Functional Responses and Interference within and between Year Classes of a Dragonfly Population. Journal of the North American Benthological Society, 8, 211-221.
 Rahmi, E., Darti, I. and Suryanto, A. (2021) A Modified Leslie-Gower Model Incorporating Beddington-DeAngelis Functional Response, Double Allee Effect and Memory Effect. Fractal and Fractional, 5, 84.