`Journal of Applied MathematicsVolume 2013 (2013), Article ID 851028, 8 pageshttp://dx.doi.org/10.1155/2013/851028`
Research Article

## Positive Steady States of a Strongly Coupled Predator-Prey System with Holling-() Functional Response

1College of Science, Xi'an Technological University, Xi'an 710032, China
2Institute of Mathematics, Shaanxi Normal University, Xi'an 710062, China

Received 25 December 2012; Accepted 2 March 2013

Copyright © 2013 Xiao-zhou Feng and Zhi-guo Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper discusses a predator-prey system with Holling-() functional response and the fractional type nonlinear diffusion term in a bounded domain under homogeneous Neumann boundary condition. The existence and nonexistence results concerning nonconstant positive steady states of the system were obtained. In particular, we prove that the positive constant solution is asymptotically stable when the parameter k satisfies some conditions.

#### 1. Introduction

In this paper, we are interested in the positive steady states of the strongly coupled predator-prey system with Holling- functional response. The specific system is as follows: where ; is a bounded domain in with smooth boundary ; is the outward directional derivative normal to ; and stand for the densities of the prey and predator; the given coefficients , , , , and are positive constants. The term is named Holling- functional response [1, 2]. In the second equation, the fractional type nonlinear diffusion term models a situation in which the population pressure of the predator species weakens in high-density areas of the prey species. For more precise details, we can refer to [3, 4]. Paper [3] discusses a strongly coupled predator-prey system with nonmonotonic functional response, the existence and nonexistence results concerning nonconstant positive steady states of the system were proved by degree theory. Paper [4] considers the positive steady states for a prey-predator model with some nonlinear diffusion terms, and the sufficient conditions for the existence of positive steady state solutions were obtained by bifurcation theory.

In recent years, there has been considerable interest in the dynamics of strongly coupled reaction-diffusion systems with cross-diffusion. We point out that most efforts have concentrated on the Lotka-Volterra competition system which was proposed first by Shigesada et al. [5]. Since their pioneering work, many authors have studied population models with cross-diffusion terms from various mathematical viewpoints, for example, the global existence of time-depending solutions [611], the stability analysis for steady states [1214], and the steady state problems [1521]. In this paper, we mainly consider the existence of solutions of (1). The research method refers to [3, 22, 23].

For convenience of the research, we write (1) as the following form: where

The main work of this paper is to study the effects of the fractional type nonlinear diffusion pressures on the existence of nonconstant positive steady states of (1). Here, a positive solution means a smooth solution with both and being positive. We will demonstrate that the cross-diffusion pressure may help forming more patterns. Obviously, for system (1), one notes that when , there holds , so that the only nonnegative solutions to (1) are and . Consequently, (1) does not have any positive solution. On the other hand, when , the unique positive constant solution to (1) is ; that is,

The organization of our paper is as follows. In Section 2, we establish a priori upper and lower bounds for positive solutions of (1). In Section 3, we use a degree theory to develop a general result to enable one to conclude the existence or nonexistence of nonconstant steady-state solutions or patterns as the index of positive constant steady states changes. In Section 4, we establish the existence of nonconstant positive solutions to (1) for a large range of diffusion and cross-diffusion coefficients. Meanwhile, we prove that the positive constant solution is asymptotically stable for different ranges of parameters.

#### 2. Upper and Lower Bounds for Positive Solutions

The main purpose of this section is to give a priori upper and lower positive bounds for positive solutions of (1). Firstly, we cite two known results.

Harnack Inequality (see [24]). Let and be a positive classical solution to in subject to the homogeneous Neumann boundary condition. Then, there exists a positive constant such that .

Maximum Principle (see [25]). Suppose that . If satisfies and , then .

Theorem 1. Let and , be fixed constants. Assume that , and there exists a positive constant , such that any positive solution of (1) satisfies

Proof. Let and be defined as in (3), and denote that Then, (1) becomes
For the first equation of (9), by the Maximum Principle, we have
The function satisfies where . It is easy to verify that the norm then, is bounded by a constant depending only on and . By the Harnack inequality, we have where . It follows that
Integrating the equations of (1) over by parts and making use of the boundary conditions, we have For any , we denote that . Then, from (15) and (10), we have Along with (14), we have Similarly, consider the equation of as follows: where . Then, is bounded by a constant depending only on and . By the Harnack Inequality, there exists a positive constant such that . Then, Thus, along with (10), (14), and (17), we can complete the proof.

Theorem 2. Assume that . Let , be fixed constants. There exists a constant , such that any positive solution of (1) satisfies provided that and .

Proof. Since , there exists , such that ; that is, It follows that . Consequently, by in Theorem 1, we have where and .
In the following, we mainly prove that as . Now, we suppose that claim is not true; then, there exists a sequence with and . And the positive solution of (1) corresponding to is such that By (14), we have , and then, ; furthermore,   uniformly holds as . Then, the first equation of (1) becomes the following:
According to (7), there exists such that ,  , where is a positive constant which does not depend on . For each given in problem (24), it follows from estimate that , where . Let ; then, by Sobolev Imbedding Theorems, we get where . We choose such that Imbedding is compact, and along with elliptic equation regularity theory, there exists the subsequence of , which is still denoted by , and exists such that uniformly holds in . On the other hand, , and when , the limit of (24) becomes the following problem: Applying Maximum Principle to problem (26) and noting that , we have .
However, when , the limit of (24) becomes the following problem: which implies that for some nonnegative constant . Since , by letting and noting that ,  , and , we also have .
By a similar argument as that in (24), for the second equation of (1), we can prove that there exists a subsequence in , such that . Since , then . Dividing the second equation of (1) by , and integrating over , we have Let , and note that and ; then, . This contradiction to the assumption completes the proof.

#### 3. A Result on Degree Theory

In this section, we obtain nonexistence of nonconstant positive solutions to (1) as . Meanwhile, by degree theory, a general result to establish the existence of nonconstant positive solutions to (1) in the next section is proved.

Denote and . We will fix and take as bifurcation parameters, the dependence of will often be suppressed. Define Since and is positive for all nonnegative , exists. Hence, is a positive solution to (1) if and only if where is the inverse of in . As is a compact perturbation of the identity operator , the Leray-Schauder is well defined if for all .

By Theorems 1 and 2, there exists a positive constant such that if in , then . Note that can be taken as a continuous function for . By the invariance property of the Leray-Schauder degree, we then conclude that does not depend on if , and it also does not depend on , if changes continuously in without touching the surface .

For the case , by degree invariance, we need only consider a special ; say with large . For this we can use the following nonexistence result. To compare the existence regions of (1) with and without cross-diffusion, we give a nonexistence result stronger than what is needed here.

Theorem 3. Let be given. Then, there exists a positive constant such that when and , (1) does not have any nonconstant positive solution.

Proof. First, by (7), there exists a positive constant such that a classical solution to (1) satisfies , provided that . Now, we write as average of over , where Multiplying the equation for of (1) by and integrating over by parts, we have for some positive constant and , where and lie between and , and , respectively. Similarly, we have for some positive constant and , where and are the same as in (33). Adding (33) and (34), we obtain It follows from the Poincaré Inequality that Since we can choose such that , we may also choose and sufficiently small such that . Consequently, by (36), which implies that constant, and, hence, constant if . Thus, we complete the proof of the theorem.

In the following, we only calculate when all solutions to are positive constant solutions in .

Let be the eigenvalues of the operator in with zero flux boundary condition, be the eigenspace corresponding to in ,  ,  , an orthonormal basis of , and . Then, , where . This decomposed method is similar to that of [22].

Let be a positive root to . We can calculate where is a linear mapping from to itself.

Denote that .

Lemma 4. Let be a positive root to , and assume that for all . Then,

Proof. If is invertible, then the index of at is defined as index, where is the number of eigenvalues of with negative real parts. The is then equal to summation of the indexes over all solutions to in , provided that on . Since for each integer is invariant under , and is an eigenvalue of on if and only if is an eigenvalue of following matrix: Since ,   is invertible. Therefore, the number of eigenvalues with negative real parts of on is odd if and only if , and therefore

Theorem 5. Assume that with . Then, for every , there exists a constant such that for every ,

Proof. Since (1) does not have any positive solution in , when , we then conclude that
Next, to complete the proof of Theorem 5. We need only to calculate the degree for the case . In this case, by Theorem 3, all positive solutions to are the unique positive constant solution to , denoted by , when so that ,  . This implies that when is sufficiently large, for all , and thus . It follows from Lemma 4 that . This completes the proof.

Remark 6. The change of degree when passes the borderline is due to the appearance (disappearance) of a positive constant steady-state bifurcating from .

#### 4. Existence of Nonconstant Positive Solutions

In this section, we establish the existence of positive nonconstant solutions for (1). In particular, we show that for certain ranges of parameters where (1) does not have any positive nonconstant steady state, our model can still produce patterns. The idea is as follows. First we calculate the index of at positive constant steady states. Suppose that the sum of all these indices is not equal to the degree stated in Theorem 5. Then, in for must have a nonconstant positive solution, which also solves (1).

In the following, we always assume that , and (1) has a unique positive constant solution , when is a positive solution to . We can get the following results at by simply calculating

To calculate the roots of , we will restrict our attention to large . Note that where

The sign of the trace is determined by and the sign of the determinant .

Hence, we will discuss separately the following cases:(i) , obviously , then ;(ii) :(iia)if , then ;(iib)if , then .

##### 4.1. The Case

In this subsection, we consider local stability of the constant steady state for evolution dynamics where is a positive constant solution to by (5).

Theorem 7. Let   . Then the positive constant solution is asymptotically stable with respect to the dynamics (47). Consequently, in a small neighborhood of  , (1) does not have any nonconstant positive solution.

Proof. The linearization of (47) at takes the form Denote that the corresponding linear operator . For each eigen-pair    of on with no flux boundary condition, is a solution to (48) if and only if is an eigen-pair of the matrix . is invariant under the operator . From (44), we have Then Hence, we conclude that the two eigenvalues and of have negative real parts. More specifically, we have the following:
(i) For , , if , then ; if , then
(ii) For , as is increasing with respect to and as , it follows that if , then if , since and , then for some positive constant that does not depend on . The previous arguments show that there exists a positive constant , which does not depend on , such that , . Consequently, the spectrum of , which consists of eigenvalues, lies in . It then follows from Theorem of [26, page 98] that the constant steady-state is asymptotically stable to (47).

##### 4.2. The Case
###### 4.2.1.

Theorem 8. Let . Then the positive constant solution is asymptotically stable with respect to the dynamics (47).

Proof. Similar to the proof of Theorem 7.

###### 4.2.2.

In this case, is the only positive constant solution to . By fixing the diffusion coefficients (for prey) and using the diffusion coefficients and (for predator) as bifurcation parameters, we will show that (1) can create nonconstant positive solutions. We want to emphasize that it is caused by the presence of cross-diffusion which has a more complex role than that of the diffusion coefficients and .

Theorem 9 (existence with suitable and ). Assume that and , define and as (46).
(i) Suppose that and are given such that for some positive even integer . There exists a positive constant such that if , then (1) has at least one nonconstant positive solution.
(ii) Suppose that is given such that for some positive even integer . Then, for any given , there exists a positive constant such that if , (1) has at least one nonconstant positive solution.

Proof. Denote by , with , the two roots to , and then From (45), we see that Suppose that for some positive even integer ; then, there exists a positive constant such that , and we have Hence, if , is equivalent to , since is even. It follows from Lemma 4 that Consequently, has at least one nonconstant positive solution that is different from the constant function . Otherwise, the degree of in would be for all large enough , which would contradict Theorem 5. This proves, the first assertion of the theorem, and the second assertion is similarly proved.

Remark 10. For to be positive, it is necessary and sufficient to have When this inequality holds, is also positive provided that is large, and we then can adjust to and make the assumptions in (i) or (ii) of theorem hold.

Remark 11. In fact, if , then . This case is similar to case (i). Furthermore, the positive constant solution is also asymptotically stable with respect to the dynamics (47).

#### Acknowledgments

This project is supported by the NSF of China under Grant 11001160, the Scientific Research Plan Projects of Shaanxi Education Department (no. 12JK0865), and the President Fund of Xi’an Technological University (XAGDXJJ1136).

#### References

1. W. Wang and J. H. Sun, “On the predator-prey system with Holling-($n+1$) functional response,” Acta Mathematica Sinica, vol. 23, no. 1, pp. 1–6, 2007.
2. J. A. Johnson, “Banach spaces of Lipschitz functions and vector-valued Lipschitz functions,” Transactions of the American Mathematical Society, vol. 148, pp. 147–169, 1970.
3. X. F. Chen, Y. W. Qi, and M. X. Wang, “A strongly coupled predator-prey system with non-monotonic functional response,” Nonlinear Analysis. Theory, Methods & Applications, vol. 67, no. 6, pp. 1966–1979, 2007.
4. T. Kadota and K. Kuto, “Positive steady states for a prey-predator model with some nonlinear diffusion terms,” Journal of Mathematical Analysis and Applications, vol. 323, no. 2, pp. 1387–1401, 2006.
5. N. Shigesada, K. Kawasaki, and E. Teramoto, “Spatial segregation of interacting species,” Journal of Theoretical Biology, vol. 79, no. 1, pp. 83–99, 1979.
6. H. Amann, “Dynamic theory of quasilinear parabolic systems. III. Global existence,” Mathematische Zeitschrift, vol. 202, no. 2, pp. 219–250, 1989.
7. Y. S. Choi, R. Lui, and Y. Yamada, “Existence of global solutions for the Shigesada-Kawasaki-Teramoto model with strongly coupled cross-diffusion,” Discrete and Continuous Dynamical Systems A, vol. 10, no. 3, pp. 719–730, 2004.
8. L. Dung, “Cross diffusion systems on n spatial dimensional domains,” Indiana University Mathematics Journal, vol. 51, no. 3, pp. 625–643, 2002.
9. L. Dung, L. V. Nguyen, and T. T. Nguyen, “Shigesada-Kawasaki-Teramoto model on higher dimensional domains,” Electronic Journal of Differential Equations, vol. 2003, pp. 1–12, 2003.
10. Y. Lou, W.-M. Ni, and Y. Wu, “On the global existence of a cross-diffusion system,” Discrete and Continuous Dynamical Systems, vol. 4, no. 2, pp. 193–203, 1998.
11. A. Yagi, “Global solution to some quasilinear parabolic system in population dynamics,” Nonlinear Analysis. Theory, Methods & Applications, vol. 21, no. 8, pp. 603–630, 1993.
12. Y. Kan-on, “Stability of singularly perturbed solutions to nonlinear diffusion systems arising in population dynamics,” Hiroshima Mathematical Journal, vol. 23, no. 3, pp. 509–536, 1993.
13. K. Kuto, “Stability of steady-state solutions to a prey-predator system with cross-diffusion,” Journal of Differential Equations, vol. 197, no. 2, pp. 293–314, 2004.
14. Y. Wu, “The instability of spiky steady states for a competing species model with cross diffusion,” Journal of Differential Equations, vol. 213, no. 2, pp. 289–340, 2005.
15. K. Kuto and Y. Yamada, “Multiple coexistence states for a prey-predator system with cross-diffusion,” Journal of Differential Equations, vol. 197, no. 2, pp. 315–348, 2004.
16. Y. Lou and W.-M. Ni, “Diffusion vs cross-diffusion: an elliptic approach,” Journal of Differential Equations, vol. 154, no. 1, pp. 157–190, 1999.
17. Y. Lou, W.-M. Ni, and S. Yotsutani, “On a limiting system in the Lotka-Volterra competition with cross-diffusion,” Discrete and Continuous Dynamical Systems A, vol. 10, no. 1-2, pp. 435–458, 2004.
18. M. Mimura, Y. Nishiura, A. Tesei, and T. Tsujikawa, “Coexistence problem for two competing species models with density-dependent diffusion,” Hiroshima Mathematical Journal, vol. 14, no. 2, pp. 425–449, 1984.
19. K. Nakashima and Y. Yamada, “Positive steady states for prey-predator models with cross-diffusion,” Advances in Differential Equations, vol. 1, no. 6, pp. 1099–1122, 1996.
20. W. H. Ruan, “Positive steady-state solutions of a competing reaction-diffusion system with large cross-diffusion coefficients,” Journal of Mathematical Analysis and Applications, vol. 197, no. 2, pp. 558–578, 1996.
21. K. Ryu and I. Ahn, “Positive steady-states for two interacting species models with linear self-cross diffusions,” Discrete and Continuous Dynamical Systems A, vol. 9, no. 4, pp. 1049–1061, 2003.
22. P. Y. H. Pang and M. X. Wang, “Strategy and stationary pattern in a three-species predator-prey model,” Journal of Differential Equations, vol. 200, no. 2, pp. 245–273, 2004.
23. P. Y. H. Pang and M. Wang, “Qualitative analysis of a ratio-dependent predator-prey system with diffusion,” Proceedings of the Royal Society of Edinburgh A, vol. 133, no. 4, pp. 919–942, 2003.
24. C.-S. Lin, W.-M. Ni, and I. Takagi, “Large amplitude stationary solutions to a chemotaxis system,” Journal of Differential Equations, vol. 72, no. 1, pp. 1–27, 1988.
25. Y. Lou and W.-M. Ni, “Diffusion, self-diffusion and cross-diffusion,” Journal of Differential Equations, vol. 131, no. 1, pp. 79–131, 1996.
26. D. Henry, Geometric Theory of Semilinear Parabolic Equations, vol. 840 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 3rd edition, 1981.