Research Article | Open Access
Ruofeng Rao, Quanxin Zhu, Jialin Huang, "Existence, Uniqueness, and Input-to-State Stability of Ground State Stationary Strong Solution of a Single-Species Model via Mountain Pass Lemma", Complexity, vol. 2021, Article ID 8855351, 11 pages, 2021. https://doi.org/10.1155/2021/8855351
Existence, Uniqueness, and Input-to-State Stability of Ground State Stationary Strong Solution of a Single-Species Model via Mountain Pass Lemma
In this study, the authors utilize mountain pass lemma, variational methods, regularization technique, and the Lyapunov function method to derive the unique existence of the positive classical stationary solution of a single-species ecosystem. Particularly, the geometric characteristic of saddle point in the mountain pass lemma guarantees that the equilibrium point is the ground state stationary solution of the ecosystem. Based on the obtained uniqueness result, the authors use the Lyapunov function method to derive the globally exponential stability criterion, which illuminates that under some suitable conditions, a certain internal competition is conducive to the global stability of the population, and a certain amount of family planning is conducive to the overall stability of the population. Most notably, the regularity technique of weak stationary solution employed in this study can also be applied to some existing literature related with time-delays reaction-diffusion systems for the purpose of regularization of weak solutions. Finally, an illuminative numerical example shows the effectiveness of the proposed methods.
The logistic system is one of the most classical models in ecology and mathematics, which is very important to the development of ecology . It is usually expressed aswhere represents the density or quantity of population at time , and and represent the intrinsic growth rate of population and environmental capacity, respectively. In 2011, Xiaoling Zou and Ke Wang investigated the long time behavior of the following stochastic ecosystem for single species (, Theorem 2):where and describe the growth rate and the intraspecific competition; measures the intensity of the environmental disturbances. In recent years, model (2) has been widely adopted in many applications (e.g., [3–5]). A large number of facts have shown that the spatial scale and structure of the environment can affect population interaction [6, 7] and community composition . In the landmark document , Kellam gave a large number of observations, which had a profound impact on the study of spatial ecology. First, he linked random walk with diffusion equation. The former is a description of the individual movement of some theoretical biological species, and the latter is a description of the density distribution of biological populations. He uses the data of muskrat transmission in Central Europe to prove that this connection is reasonable for small animals. Second, he combines diffusion with population dynamics and effectively introduces the reaction-diffusion equation into theoretical ecology.
In recent years, many dynamical systems, including reaction systems, have been considered as the theoretical branches of dynamical systems [10, 11, 12]. In addition, the competition within the population is the participation of the adult population, and there is a period from infancy to adulthood. At the same time, this time-delay problem is affected by many stochastic factors such as weather, temperature, humidity, and so on. Besides, in real life, the factors that affect population growth do not change only at a fixed time but also occur randomly. When these factors occur, the system will also change randomly. As is well known, the phenomenon of population clustering is widespread in nature, which is likely to be completely affected by environmental factors and human factors. In this case, the growth curve of mosquitoes or small fish will be different from the previous form. This phenomenon can be expressed as a switch between two environmental states because the switching between different environments is not memory free; therefore, one can use continuous time Markov chain to model the situation of environment switching [13, 14, 15].
Inspired by some ideas and methods of related literature [16–31], particularly [17, 18], we are to investigate the stability of stationary density of a single-species model with diffusion and delayed feedback under natural state. This study has the following highlights:(i)As far as we know, it is the first study to investigate the stability of stationary density of a single-species model with diffusion and delayed feedback under the Dirichlet zero boundary value. And the Dirichlet boundary value can well simulate the fact that the species lives in its biosphere, while the population density tends to zero at the boundary of biosphere due to the harsh condition. Different from exiting literature involved in Neumann zero boundary value which implies that there is no animal migration at the edge of the biosphere, the Dirichlet boundary value of this study admits the fact of animal migration, but no animals under study live on the border for a long time.(ii)It is the first comprehensive application of mountain pass lemma, variational technique, and the Lyapunov function method to derive the unique existence of globally exponentially stable positive stationary solution of a single-species model with diffusion and delayed feedback under the Dirichlet zero boundary value.(iii)The obtained stability criterion illuminates that under some suitable conditions, a certain internal competition is conducive to the overall stability of the population, and a certain amount of family planning is conducive to the overall stability of the population.(iv)Different from many existing literature related with the global stability of discontinuous systems [30–32] and time-delays reaction-diffusion systems [33, 34], the weak stationary solution is regularized in this study. Most notably, the regularity technique of weak stationary solutions can also be applied to such existing literature [24, 25, 33, 34] for the purpose of regularization of weak solutions.
In next sections, the authors have made the following arrangement:
In Section 2, the authors present some descriptions on the ecosystem, and some necessary definitions and lemmas are presented. In Section 3, the authors utilize the existence technique and regularity method employed in  to derive the existence of positive strong stationary solution of the ecosystem. Moreover, utilizing the uniqueness technique used in [18, 35] results in the unique existence of the stationary solution. Finally, the Lyapunov function method is applied to derive the stability criterion. In Section 4, numerical example and comparisons are given. And in final section, conclusions and further considerations are proposed.
For convenience, throughout of this study, we denote by the first positive eigenvalue of Laplace operator in . Denote . Denote by , the norm of , and by , the first positive eigenvalue of Laplace operator in . Besides, we denote for and for matrix .
2. System Descriptions
Denote by the complete probability space with a natural filtration . Let and the random form process be a homogeneous, finite-state Markovian process with right continuous trajectories with generator and transition probability from mode at time to mode at time ,where is the transition probability rate from to and , and .
Consider the following ecosystem with diffusion and delayed feedback,where is the external input, and describe the growth rate and the intraspecific competition, and is the bounded domain with its boundary and is also a domain in (e.g., ). It is also suitable to the case that the species lives in two dimensional plane (, Remark 2.1).
Remark 1. Because only the adult is competitive for survival and there is a mature period from the larva to the adult, we consider the time-delayed system in this study, which is better to simulate this maturity problem. Particularly, the gain coefficient of time-delay feedback can be derived from the statistical data of the observed system.
Throughout this study, we assume (H1) the positive function is only a microperturbation. That is, there exists a positive number small enough such thatwhere with and being a pair of coprime odd numbers. And is continuous and for all . Here, is the Sobolev critical exponent. In addition, for all .
Remark 2. In essence, the restrictive condition with and is set to ensure that the function is odd, so that its primitive function is even. Thereby, we can also assume in H1 that the function is odd.
Let be a stationary solution of system (4) implies that is a solution of the following equation:Of course, each solution of equation (6) must be one of the solutions of system (4).
Lemma 1. Consider the following equation:where is a domain of , and satisfies the following conditions:(a)There exists such that for any given positive number ,(b)If and , or and , then (as ) holds uniformly for (c).
Then, the solution of equation (6) in is the strong solution. In addition, for .
Lemma 2. Let . Then, there is a conclusion that . Besides, if , if , if , and if . In addition, if and if . Besides, if .
Lemma 3 (Mountain pass lemma without the PS condition). Let be a Banach space, , satisfying , and there exists such that . Besides, there is , such that . Let be the set of all paths connecting 0 and , that is,Set
Then, , and possesses a critical sequence on .
Remark 3. Lemma 3 is the mountain pass lemma without the PS condition (e.g., ). If, in addition, satisfies the PS condition, then is a critical value of . Besides, let be the functional corresponding to equation (6), then must be a ground state solution of equation (6) if is a critical point of the functional with defined in (10) of Lemma 3.
3. Main Result
First, we may present the existence of a stationary strong solution of system 2.1. In addition, it is necessary to guarantee that and , which may be proved as follows:
Theorem 1. Suppose the condition H1 holds, and if
Then, there is a ground-state strong stationary solution for system (4).
Proof. Let be a positive stationary solution of system (4), satisfyingwhose functional iswhere is a constant, and . It is obvious that , and a critical point of the functional is corresponding to the solution of equation (12). Next, we claim that satisfies the condition of the mountain road geometry. In fact, obviously .
The microperturbation condition (H1) yields that there are there positive constants with such thatorwhich impliesand then,Moreover,Next, (13) and (18), Poincare inequality and Sobolev embedding theorem yield that there is a positive constant such thatBesides, (11) and small lead towhich together with means that we can setsuch that .
On the other hand, it follows by (15) and thatwhich impliesorWe may select with , and then,Let with be the eigenfunction corresponding to the first positive eigenvalue (e.g., ), and set ; then, if , so that there exists satisfying and . And then, satisfies the condition of the mountain road geometry. According to mountain pass lemma, let be the set of all paths connecting 0 and . That is,SetThen, , and possesses a critical sequence on , say with and in . That is, for any given , there exists big enough such thatandwhere represents such an infinitesimal that when .
(17) yieldsSimilar as the methods of , (24)–(27), employing (29)–(33) results inwhere both are the positive constants, and due to the small , which means the boundedness of . Due to , equation (12) is the subcritical growth. It is a routine proof of the fact that is sequently compact, say in and , which implies . Besides, is the critical point of , so thatIn (35), setting , Lemma 2 leads towhich implies that a.e. . Now, we have proved that and .
Similar as that of , now we claim that the abovementioned is the strong solution.
Indeed, is the nonnegative solution of the following Dirichlet problem:whereIt is easy from the assumptions on to verify that satisfies the conditions (a)–(c); then, Lemma 1 yields that is the strong solution.
Set . Since is a stationary solution of system (4), system (4) is equivalent to the following system:whereObviously, of system (4) is corresponding to the zero solution of system (38).
Equip system (38) with the initial value,Moreover, we give some suitable assumptions as follows:
H2. There are positive numbers , such thatH3. There is a positive number , such that
Remark 4. Everyone knows the fact that the population density of any species must have the bounded below or the species will die out. For example, when the population density of whales is lower than a certain degree, it will be difficult for male and female whales to meet each other in the vast sea, leading to the extinction of the species. Besides, due to the limited resource, the population density of any species must have supper boundedness. So, the condition H2 is a suitable assumption.
There are some techniques on the existence and uniqueness of positive stationary solution in the proofs (, Theorems 1-2), which are also employed in the proofs of (, Theorems 1-2). But the methods used in the proof of Theorem 1 in this study are different from those of both [18, 35]. Besides, we are willing to consider similarly the uniqueness of the positive stationary solution in this study.
Proof. Let both are the stationary solutions of system (4). First, the condition H2 yieldsSince both are the stationary solutions of system (4), then Poincare inequality yieldswhich proves , and the proof is completed.
Theorem 3. Suppose the conditions H1–H3 and (45) hold and if there are positive numbers , such thatthen the null solution of system (38) with the initial value (40) is the globally exponential input-to-state stability; at the same time, is the globally exponential input-to-state stability at the convergence rate , where , , and is the unique positive solution of .
Proof. Consider the following Lyapunov function:where is a positive number for each .
First, H2 yieldsH3 yieldsLet be the weak infinitesimal operator (e.g., ), such thatwhere .
Now, (, Lemma 2) yields thatwhere , This completes the proof.
Remark 5. In this study, we employ mountain pass lemma and variational technique to derive the existence of positive stationary solution, which is different from the methods in . Particularly, ground-state solution is more suitable to practical engineering (e.g., [36–40]). Besides, the equilibrium points of ecosystems with the Neumann boundary value are always constants solutions, while equilibrium points of the system (4) are always the nontrivial solutions of the corresponding elliptic equation, which need the existence criterion of the solutions for the elliptic equation. And it adds the computation complexity of the results obtained in this study. Also, our model and method are different from those in [41–46].
4. Numerical Example
Example 1. In system (4), let , then (, Remark 14). , then , and condition (11) holds. Setwhere , then H1 holds. Let , then direct computation yields , and obviously, (45) holds. Then, Theorems 1-2 yield is the unique stationary solution of system (4), and is the ground state stationary strong solution of system (4).
Moreover, set , and , , and direct computation yields (46) holds for , and . According to Theorem 3, the null solution of system (38) with the initial value (40) is globally exponential input-to-state stability; at the same time, is globally exponential input-to-state stability at the convergence rate .
In Example 1, replacing with and other data unchanged, direct computation yields the convergence rate .
In example 1, replacing with and other data unchanged, direct computation yields the convergence rate .
Remark 7. Table 2 provides that a certain degree of inhibition and competition within the population is beneficial to the overall stability of the population for the natural ecosystem of a single-species model. Figure 2 shows that the bigger the intrapopulation competition intensities, the faster the stability of a single-species system.
In Example 1, replacing with and other data unchanged, direct computation yields the convergence rate .
Remark 9. Different from [2, 3], the population density boundedness of the species was considered in this study due to the important factor (Remark 4 for details). In Example 1, authors assume and , and we can see it from Figures 1–3 that .
In this study, critical point theory and variational methods are utilized to derive the unique existence of stationary solution of the single -species model, which is positive and strong. Moreover, the geometric characteristic of saddle point in the mountain pass lemma guarantees that the positive strong stationary solution is the ground state one. Moreover, the method of Lyapunov function yields the global exponential stability of the ground-state classical positive stationary solution which is the unique stationary solution of the ecosystem.
Besides, impulse control reflects the human intervention in the natural ecosystem (e.g., , Theorem 3); we may consider the next study on impulsive stabilization of the single-species ecosystem of this study. In addition, we propose mathematical conjectures that (, Problem 4) and (, Problem 1) may be correct even in the case of the stochastic differential system. In this study, we propose the mathematical conjecture that under suitable conditions, small diffusions that can make the unique stable equilibrium point of the delayed ordinary differential system become multiple stationary solutions of its corresponding partial differential system, even in the case of the stochastic model. Moreover, how to give a global stability invariance criterion of a stochastic model similar as (, Corollary 3.4)? This is an interesting problem. Finally, we may consider an interesting application to epidemic control (e.g., ). This study is involved in a single-species biological dynamic system. If we investigate the dynamics of a single-species infectious disease model, it is another interesting problem, for the human society may approximately be regarded as a single-species model.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Ruofeng Rao wrote the article. Quanxin Zhu and Huang Jialin guided the writing of the article and gave their contributions on the revision, being the corresponding authors in charge of the communication of this article. All the authors typed, read, and approved the article.
The authors would like to thank the National Natural Science Foundation of China (61773217) and Chinese Science and Technology Department of Sichuan Province for funding this work through the application basic research project (2020YJ0434) and thank Chengdu Normal University for funding this work through major scientific research projects of Chengdu Normal University (CS19ZDZ01).
- L. Chen, X. Meng, and J. Jiao, Biodynamics, Science Press, Beijing, China, 2009.
- X. Zou and K. Wang, “A robustness analysis of biological population models with protection zone,” Applied Mathematical Modelling, vol. 35, no. 12, pp. 5553–5563, 2011.
- W. Ji and G. Hu, “Stability and explicit stationary density of a stochastic single-species model,” Applied Mathematics and Computation, vol. 390, Article ID 125593, 2021.
- X. Yu, S. Yuan, and T. Zhang, “Persistence and ergodicity of a stochastic single species model with allee effect under regime switching,” Communications in Nonlinear Science and Numerical Simulation, vol. 59, pp. 359–374, 2018.
- Y. Jin, “Analysis of a stochastic single species model with allee effect and jump-diffusion,” Advances in Difference Equations, vol. 165, pp. 1–11, 2020.
- G. F. Gause, The Struggle for Existence, Williams & Wilkins, Baltimore, MD, USA, 1935.
- J. Huisman, M. Arrayás, U. Ebert, and B. Sommeijer, “How do sinking phytoplankton species manage to persist?” The American Naturalist, vol. 159, no. 3, pp. 245–254, 2002.
- R. H. MacArthur and E. O. Wilson, The Theory of Island Biogeography, Princeton University Press, Princeton, NJ, USA, 1967.
- J. G. Skellam, “Random dispersal in theoretical populations,” Biometrika, vol. 38, no. 1/2, pp. 196–218, 1951.
- J. Hofbauer and K. Sigmund, Dynamical Systems and the Theory of Evolution, Cambridge University Press, Cambridge, UK, 1988.
- V. Hutson and K. Schmitt, “Permanence and the dynamics of biological systems,” Mathematical Biosciences, vol. 111, no. 1, pp. 1–71, 1992.
- H. L. Smith and H. Thieme, Dynamical Systems and Population Persistence, American Mathematical Society, Providence, RI, USA, 2011.
- R. Rudnicki and K. Pichór, “Influence of stochastic perturbation on prey-predator systems,” Mathematical Biosciences, vol. 206, no. 1, pp. 108–119, 2007.
- R. Rudnicki, “Long-time behaviour of a stochastic prey-predator model,” Stochastic Processes and Their Applications, vol. 108, no. 1, pp. 93–107, 2003.
- B. Gompertz, “On the nature of the function expressive of the law of human mortality and on a new method of determining the value of life contingencies,” Philosophical Transactions of the Royal Society of London, vol. 115, pp. 513–585, 1825.
- Y. Zhao and Q. Zhu, “Stabilization by delay feedback control for highly nonlinear switched stochastic systems with time delays,” International Journal of Robust and Nonlinear Control, vol. 31, no. 8, pp. 3070–3089, 2021.
- R. Rao, “Positive solution for the Dirichlet zero-boundary value problem (In Chinese),” College Mathematics, vol. 26, no. 2, pp. 146–152, 2010.
- R. Rao, J. Huang, and X. Li, “Stability analysis of nontrivial stationary solution and constant equilibrium point of reaction-diffusion neural networks with time delays under Dirichlet zero boundary value,” Neurocomputing, vol. 445, pp. 105–120, 2021.
- R. Rao, “Stability and stabilization of ecosystem for epidemic virus transmission under Neumann boundary value via Impulse control,” Preprints, Article ID 2021020197, 2021.
- W. Lu, Variational Methods in Differential Equations, Science Press, Beijing, China, 2003.
- R. Rao, S. Zhong, and X. Wang, “Stochastic stability criteria with LMI conditions for Markovian jumping impulsive BAM neural networks with mode-dependent time-varying delays and nonlinear reaction-diffusion,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 1, pp. 258–273, 2014.
- R. Rao, “Global stability of a Markovian jumping chaotic financial system with partially unknown transition rates under impulsive control involved in the positive interest rate,” Mathematics, vol. 7, no. 7, 2019.
- M. Zhang and Q. Zhu, “Stability analysis for switched stochastic delayed systems under asynchronous switching: a relaxed switching signal,” International Journal of Robust and Nonlinear Control, vol. 30, no. 18, pp. 8278–8298, 2020.
- U. Humphries, G. Rajchakit, R. Sriraman et al., “An extended analysis on robust dissipativity of uncertain stochastic generalized neural networks with Markovian jumping parameters,” Symmetry, vol. 12, no. 6, p. 1035, 2020.
- P. Chanthorn, G. Rajchakit, U. Humphries, P. Kaewmesri, R. Sriraman, and C. Lim, “A delay-dividing approach to robust stability of uncertain stochastic complex-valued hopfield delayed neural networks,” Symmetry, vol. 12, no. 5, 2020.
- U. Humphries, G. Rajchakit, P. Kaewmesri et al., “stochastic memristive quaternion-valued neural networks with time delays: an analysis on mean square exponential input-to-state stability,” Mathematics, vol. 8, no. 5, 2020.
- P. Chanthorn, G. Rajchakit, J. Thipcha et al., “Robust stability of complex-valued stochastic neural networks with time-varying delays and parameter uncertainties,” Mathematics, vol. 8, no. 5, p. 742, 2020.
- M. S. Ali, M. Usha, Q. Zhu, and S. Shanmugam, “Synchronization analysis for stochastic T-S fuzzy complex networks with Markovian jumping parameters and mixed time-varying delays via impulsive control,” Mathematical Problems in Engineering, vol. 2020, p. 9739876, 2020.
- F. Kong, Q. Zhu, R. Sakthivel, and A. Mohammadzadeh, “Fixed-time synchronization analysis for discontinuous fuzzy inertial neural networks with parameter uncertainties,” Neurocomputing, vol. 422, pp. 295–313, 2021.
- W. Cao and Q. Zhu, “Razumikhin-type theorem for pth exponential stability of impulsive stochastic functional differential equations based on vector Lyapunov function,” Nonlinear Analysis: Hybrid Systems, vol. 39, p. 100983, 2021.
- F. Kong and Q. Zhu, “New fixed-time synchronization control of discontinuous inertial neural networks via indefinite Lyapunov-Krasovskii functional method,” International Journal of Robust and Nonlinear Control, vol. 31, no. 2, pp. 471–495, 2020.
- F. Kong, Q. Zhu, and T. Huang, “New fixed-time stability lemmas and applications to the discontinuous fuzzy inertial neural networks,” IEEE Transactions on Fuzzy Systems, p. 1, 2020.
- R. Song, B. Wang, and Q. Zhu, “Delay‐dependent stability of nonlinear hybrid neutral stochastic differential equations with multiple delays,” International Journal of Robust and Nonlinear Control, vol. 31, no. 1, pp. 250–267, 2021.
- L. Gao, Z. Cao, M. Zhang, and Q. Zhu, “Input-to-state stability for hybrid delayed systems with admissible edge-dependent switching signals,” Journal of the Franklin Institute, vol. 357, no. 13, pp. 8823–8850, 2020.
- R. Rao, Q. Zhu, and K. Shi, “Input-to-State stability for impulsive gilpin-ayala competition model with reaction diffusion and delayed feedback,” IEEE Access, vol. 8, pp. 222625–222634, 2020.
- X. Yang and Q. Zhu, “Stabilization of stochastic functional differential systems by steepest descent feedback controls,” IET Control Theory & Applications, vol. 15, no. 6, pp. 805–813, 2021.
- H.-S. Ding and Q. He, “Existence and asymptotic behavior of positive ground state solutions for nonlinear Kirchhoff problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 90, Article ID 105369, 2020.
- S. Tian, “Non-degeneracy of the ground state solution on nonlinear Schrödinger equation,” Applied Mathematics Letters, vol. 111, Article ID 106634, 2021.
- X. Chen, X. Wu, and C. Tang, “existence of a positive ground state solution for a class of asymptotically 3-linear Kirchhoff-type equations (in Chinese),” Journal of Southwest China Normal University(Natural Science Edition), vol. 44, no. 4, pp. 22–25, 2019.
- Y. Li, G. Li, and C. Tang, “Existence and concentration of ground state solutions for Choquard equations involving critical growth and steep potential well,” Nonlinear Analysis, vol. 200, Article ID 111997, 2020.
- R. Rao and S. Zhong, “Impulsive control on delayed feedback chaotic financial system with Markovian jumping,” Advances in Difference Equations, vol. 50, 2020.
- X.-X. Liao and J. Li, “Stability in Gilpin-Ayala competition models with diffusion,” Nonlinear Analysis: Theory, Methods & Applications, vol. 28, no. 10, pp. 1751–1758, 1997.
- L. Bai and K. Wang, “Gilpin-Ayala model with spatial diffusion and its optimal harvesting policy,” Applied Mathematics and Computation, vol. 171, no. 1, pp. 531–546, 2005.
- Y. Liu and Y. Tao, “Dynamics in a parabolic-elliptic two-species population competition model with cross-diffusion for one species,” Journal of Mathematical Analysis and Applications, vol. 456, no. 1, pp. 1–15, 2017.
- K. Ding and Q. Zhu, “Fuzzy intermittent extended dissipative control for delayed distributed parameter systems with stochastic disturbance: a spatial point sampling approach,” IEEE Transactions on Fuzzy Systems, p. 1, 2021.
- X. Liu and Q. Zhu, “Stochastically globally exponential stability of stochastic impulsive differential systems with discrete and infinite distributed delays based on vector Lyapunov function,” Complexity, vol. 2020, Article ID 7913050, 16 pages, 2020.
Copyright © 2021 Ruofeng Rao et al. 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.