Hopf Bifurcations in a Predator-Prey System of Population Allelopathy with Discrete Delay

A delayed Lotka-Volterra two-species predator-prey system of population allelopathy with discrete delay is considered. By linearizing the system at the positive equilibrium and analyzing the associated characteristic equation, the asymptotic stability of the positive equilibrium is investigated and Hopf bifurcations are demonstrated. Furthermore, the direction of Hopf bifurcation and the stability of the bifurcating periodic solutions are determined by the normal form theory and the center manifold theorem for functional differential equations (FDEs). Finally, some numerical simulations are carried out for illustrating the theoretical results.


t x t x t r a x t a y t y t y t r a G t s x s s a y t
The delay kernel function G(s) may take the so-called "weak" generic kernel function G(s) =   (α > 0), where the "weak" generic kernel implies that the importance of events in the past simply decreases exponentially the further one looks into the past while the "strong" generic kernel implies that a particular time in the past is more iportant than any other [13].When G(s) takes the "weak" generic kernel function and the "strong" generic kernel function G(s) = 2 s e       1 2 a y t t   (α > 0), properties of the stability of the positive equilibrium of system (2) and Hopf bifurcations of nonconstant periodic solutions have been investigated respectively by using the normal form theory and the center manifold reduction for FDEs [14,15].See [5,16] for details.
x t x t r a x t y t y t r a x a y t He [17] and Lu andWang [18] investigated the stability of the positive equilibrium of the system, and they found that the positive equilibrium is globally asymptotically stable for any values of delays τ and η when the coefficients of the system satisfy the condition 11 22 a a  tra-specific competitions, one can refer to [10,11].For Latka-Volterra two species competition model, an important observation made by many works is that increased population of one species might affect the growth of another species by the production of allelopathic toxins or stimulators, thus influencing seasonal succession [19].For instance, Maynard Smith [20] incorporated the effect of toxic substances in a two species Lotka-Volterra competitive system by considering that each species produce a substance toxic to the other but only when the other is present.Then the Lotka-Volterra competitive system was modified into the form: Chattopadyay [21] studied the stability properties of the above system, although the study contains the flaw of ignoring an important delay factor in the system.[22] suggested "In reality, a species needs sometime for maturity to produce a substance which will be toxic (or stimulatory) to the other and hence a delay term in the system arises".Mukhopadhyay, Chattpadhyay and Tapaswi [22] introduced a delay in system (4), which leads the following form When γ 1 > 0, γ 2 > 0, the model system (5) represents an allelopathic inhibitory system, each species producing a substance toxic to the other; when γ 1 < 0, γ 2 < 0, (5) repr esents an allelopathic stimulatory system, each species producing a substance stimulatory to the growth of the other species.
Similar phenomenon also exist in predator-prey model.Rice [19] has suggested that "all meaningful, functional ecological models will eventually have to include a category on allelopathic and other allelochemic effects".To our knowledge, such viewpoint haven't been investigated in predator-prey model so far.On the other hand, some species in the real nature world may produce substances which are toxic or stimulatory to the others while they themselves do not experience any reciprocal effects during the process of predation.For example, some species of poisonous snake release toxic substance to control prey.The production of toxic substance by the predator species will not be instantaneous, but mediated by some time lag, see [7,[19][20][21][22]. From this viewpoint and com-bining the factors appeared above of different type of time delay and allelopathic effect in predator-prey model, we have modified the model of (3).Therefore, by considering that one species produced a substance toxic to the other during the process of predation, but only when the other is present.Then the system (3) can be written as , .
x t x t k x t y t x t y t y t y t k x t y t where , 0   e have investigated the bifurcation behavior on time delay of this modified dynamical system (6).It has also been observed that time delay can drive the competitive system to sustained oscillations, as shown by Hopf bifurcation analysis and limit cycle stability.Hence interaction between the time delay effect produced by delayed toxin can regulate the densities of different competing species in the aquatic ecosystem, thus influencing seasonal successsion, blooms and pulses.To the best of our knowledge no such attempts have been taken to include interaction between the time delay effect produced by delayed toxin in a predator-prey system.Therefore, this research might behelpful to the study of predator-prey model and related problem in biological system.
. W This paper is organized as follows.In Section 2, by linearizing the resulting two-dimensional system at the positive equilibrium and analyzing the associated characteristic equation, it is found that under suitable conditions on the parameters the positive equilibrium is asymptotically stable when the delay is less than a certain critical value and unstable when the delay is greater than this critical value.Meanwhile, according to the Hopf bifurcation theorem for FDEs, we find that the system can also undergo a Hopf bifurcation of nonconstant periodic solution at the positive equilibrium when the delay crosses through a sequence of critical values.In Section 3, to determine the direction of the Hopf bifurcations and the stability of bifurcated periodic solutions occurring through Hopf bifurcations, an explicit algorithm is given by applying the normal form theory and the center manifold reduction for FDEs developed by Hassard, Kazarinoff and Wan [23].To verify our theoretical predictions, some numerical simulations are also included in Section 4.

Stability of Equilibria and Existence of Hopf Bifurcations
The state of equilibria of the system (6) for τ = 0 are as follows: where    .
is a unique positive equilibrium when the condition (H1) 1 2 2 1 0 k k     holds.Throughout this section, we always assume that the condition H(1) holds.
Clearly, the characteristic equation of the linearized system of system (6) at the equilibrium is which has two real roots, 1 , 2 .Therefore, the equilibrium is unstable and is a saddle point of system (6).The linearized system of system (6) at the Therefore, the equilibrium , is an unstable node of system (6).The characteristic equation at the the condition (H1), Equation ( 7) has a negative real root and a positive real root  is unstable and is also a saddle point of system (6) when the condition (H1) is satisfied.
In what follows, we investigate the stability of the positive equilibrium of system (6).

E x y  
Underthe assumption (H1), let x t y t y    .Then system ( 6) is equivalent to the following two dimensional system: and the positive equilibrium of system ( 6) is transformed into the zero equilibrium (0, 0) of system (8).It is easy to see that the characteristic equation of the linearized system of system (8) at the zero equilibrium (0, 0) is where It is well known that the stability of the zero equilibrium (0, 0) of system ( 8) is determined by the real parts of the roots of Equation ( 9).If all roots of Equation ( 9) locate the left-half complex plane, then the zero equilibrium (0, 0) of system ( 8) is asymptotically stable.If Equation ( 9) has a root with positive real part, then the zero solution is unstable.Therefore, to study the stability of the zero equilibrium (0, 0) of system (8), an important problem is to investigate the distribution of roots in the complex plane of the characteristic Equation (9).
For Equation ( 9), according to the Routh-Hurwitz criterion, we have the following result.
Next, we consider the effects of a positive delay τ on the stability of the zero equilibrium (0, 0) of system (8).Since the roots of the characteristic Equation ( 9) depend continuously on τ, a change of τ must lead to a change of the roots of Equation (9).If there is a critical value of τ such that a certain root of (9) has zero real part, then at this critical value the stability of the zero equilibrium (0, 0) of system (8) will switch, and under certain conditions a family of small amplitude periodic solutions can bifurcate from the zero equilibrium (0, 0); that is, a Hopf bifurcation occurs at the zero equilibrium (0, 0).Now, we look for the conditions under which the characteristic Equation ( 9) has a pair of purely imaginary roots, see [24].Clearly, iω(ω > 0) is a root of Equation ( 9) if and only if ω satisfies the following equation: Separating the real and imaginary parts of the above equation yields the following equations: Adding up the squares of the corresponding sides of the above equations yields equations with respect to ω: (11) has no positive real root.Otherwise (11) has an unique positive root, sign it as in the following.From the first equation of ( 10), we know that the value of τ associated with If we define then when Equation ( 9) has a pair of purely imaginary roots ±  0,1, 2, , Let λ(τ) = α(τ) + iω(τ) be a root of Equation ( 9) near τ = j  satisfying α( j  ) = 0 and ω( j  ) = 0  .For this pair of conjugate complex roots,we have the following result.
Proof.Differentiating both sides of Equation ( 9) with respect to τ, and noticing that λ is a function with respect to τ, we have   From the above equation, one can easily obtain cos 2 2 sin 2 , sin 2 2 cos 2 .
j Thus, we have Combining (10) and some simple computations show that This completes the proof.
From the above discussion and the Hopf bifurcation theorem of FDEs [14,23], we can obtain the following results on the stability of the zero equilibrium of system (8); that is, the stability of the positive equilibrium of system (6).
Theorem 2.3.Suppose that the coefficients i , k i  (i = 1, 2) in system (6) satisfy the condition (H1) and 0 , 0 satisfies the condition (H2); then the following results hold. 1) The positive equilibrium is asymptotically stable when and unstable when 0

Properties of Hopf Bifurcations
In the previous section, we studied mainly the stability of the positive equilibrium of system (6) and the existence of Hopf bifurcations at the positive equilibrium In this section, we shall study the properties of the Hopf bifurcations obtained by Theorem 2.3 and the stability of bifurcated periodic solutions occurring through Hopf bifurcations by using the normal form theory and the center manifold reduction for retarded functional differential equations (RFDEs) due to Hassard, Kazarinoff and Wan [23].To guarantee the existence of the above Hopf bifurcations, throughout this section, we always assume that the conditions (H1) and (H2).Under these conditions, for fixed   0,1, 2, j   , let τ = j  + μ; then μ = 0 is the Hopf bifurcation value of system (6) at the positive equilibrium   , E x y   .Since system (6) is equivalent to system (8), in the following discussion we shall consider mainly system (8).
In system (8), let and drop the bars for simplicity of notation.Then system (8) can be rewritten as a system of RFDEs in C([−1, 0], R2) of the form Define the linear operator and the nonlinear operator and respectively, where , and let . By the Riesz representation theorem, there exists a 2 × 2 matrix function η(θ, μ), −1 ≤ θ ≤ 0, whose elements are of bounded variation such that In fact, we can choose where and Then system ( 14) is equivalent to For , define and a bilinear inner product where η(θ) = η(θ, 0).Then A(0) and A * are adjoint operators.In addition, from Section 2 we know that ± 0 j i  are eigenvalues of A(0).Thus, they are also eigenvalues of A * .Let q(θ) is the eigenvector of A(0) corresponding to From the above discussion, it is easy to know that Thus, we can easily obtain in the direction of q  and q  .Note that W is real if t x is real.We consider only real solutions.For solution 0 t x C  of ( 14), since μ = 0, which assures that Using the same notations as in Hassard, Kazarinoff, and Wan [23], we first compute the coordinates to describe the center manifold 0 at μ = 0. Let C t x be the solution of Equation ( 14) when μ = 0. Define that is where   Then it follows from ( 23) that On the center manifold we have z and z are local coordinates for center manifold 0 C It follows together with ( 16) that  ( 1) Comparing the coefficients with (27), we obtain  From ( 19) and ( 23), we have where Substituting the corresponding series into (29) and comparing the coefficients, we obtain From (29), we know that for H z z q f q q f q g z z q g z z q Comparing the coefficients with (30) gives that and From ( 31) and (33), we get Similarly, from (31) and (34), we have In what follows we shall seek appropriate and .
From the definition of A and (31) that and where     Substituting ( 35) and (39) into (37), and noticing that Similarly, substituting (35) and ( 40) into (38), we get It follows from ( 35), ( 36), (42), and (43) that g 21 can be expressed.Thus, we can compute the following values: which determine the quantities of bifurcating periodic solutions at the critical value j ); and 2 determines the period of the bifurcating periodic solutions: the period increase (decrease) if 2 ( 2 T ).Further, it follows from Lemma 2.2 and (44) that the following results about the direction of the Hopf bifurcations hold.

Numerical Simulations
In this section, we give some numerical simulations for a special case of system (6) to support our analytical results in this paper.As an example, we consider system (6) with the coefficients , 1


) can undergo a Hopf bifurcation at the positive equilibrium   , E x y   ; that is, a family of nonconstant periodic solutions can bifurcate from the positive equi-(j = 0, 1, 2, 3, •••).

20 W  and 11 W
   in 21g , we still need to compute them.

 2 
. That is, 2  determines the directions of the Hopf bifurcation: If 2 0 , then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for j determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions in the center manifold are stable (unstable)

Theorem 3 . 1 .
Suppose that (H1), (H2) hold.If ( ), then system (6) can undergo a supercritical (subcritical) Hopf bifurcation at the positive equilibrium when τ crosses through the critical values addition, the bifurcated periodic solutions occurring through Hopf bifurcations are orbitally asymptotically stable on the center manifold if and unstable if