Research Article | Open Access

Xiao-zhou Feng, Zhi-guo Wang, "Positive Steady States of a Strongly Coupled Predator-Prey System with Holling-() Functional Response", *Journal of Applied Mathematics*, vol. 2013, Article ID 851028, 8 pages, 2013. https://doi.org/10.1155/2013/851028

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

**Academic Editor:**Junjie Wei

#### 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 [6–11], the stability analysis for steady states [12–14], and the steady state problems [15–21]. 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

- 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. View at: Publisher Site | Google Scholar | MathSciNet - 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. View at: Publisher Site | Google Scholar | MathSciNet - 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. View at: Publisher Site | Google Scholar | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - N. Shigesada, K. Kawasaki, and E. Teramoto, “Spatial segregation of interacting species,”
*Journal of Theoretical Biology*, vol. 79, no. 1, pp. 83–99, 1979. View at: Publisher Site | Google Scholar | MathSciNet - H. Amann, “Dynamic theory of quasilinear parabolic systems. III. Global existence,”
*Mathematische Zeitschrift*, vol. 202, no. 2, pp. 219–250, 1989. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Dung, “Cross diffusion systems on
*n*spatial dimensional domains,”*Indiana University Mathematics Journal*, vol. 51, no. 3, pp. 625–643, 2002. View at: Publisher Site | Google Scholar | MathSciNet - 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. View at: Google Scholar - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Y. Lou and W.-M. Ni, “Diffusion, self-diffusion and cross-diffusion,”
*Journal of Differential Equations*, vol. 131, no. 1, pp. 79–131, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - D. Henry,
*Geometric Theory of Semilinear Parabolic Equations*, vol. 840 of*Lecture Notes in Mathematics*, Springer, Berlin, Germany, 3rd edition, 1981. View at: MathSciNet

#### Copyright

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.