# American Institute of Mathematical Sciences

• Previous Article
Predator – Prey/Host – Parasite: A fragile ecoepidemic system under homogeneous infection incidence
• DCDS-B Home
• This Issue
• Next Article
The motion of weakly interacting localized patterns for reaction-diffusion systems with nonlocal effect
January  2021, 26(1): 191-216. doi: 10.3934/dcdsb.2020263

## Chaotic dynamics in a simple predator-prey model with discrete delay

 1 Department of Mathematics, Columbus State University, Columbus, Georgia 31907, USA 2 Department of Mathematics and Statistics, McMaster University, Hamilton, Ontario, Canada L8S 4K1, Canda

* Corresponding author: Gail S. K. Wolkowicz

Received  March 2020 Revised  July 2020 Published  August 2020

A discrete delay is included to model the time between the capture of the prey and its conversion to viable biomass in the simplest classical Gause type predator-prey model that has equilibrium dynamics without delay. As the delay increases from zero, the coexistence equilibrium undergoes a supercritical Hopf bifurcation, two saddle-node bifurcations of limit cycles, and a cascade of period doublings, eventually leading to chaos. The resulting periodic orbits and the strange attractor resemble their counterparts for the Mackey-Glass equation. Due to the global stability of the system without delay, this complicated dynamics can be solely attributed to the introduction of the delay. Since many models include predator-prey like interactions as submodels, this study emphasizes the importance of understanding the implications of overlooking delay in such models on the reliability of the model-based predictions, especially since temperature is known to have an effect on the length of certain delays.

Citation: Guihong Fan, Gail S. K. Wolkowicz. Chaotic dynamics in a simple predator-prey model with discrete delay. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 191-216. doi: 10.3934/dcdsb.2020263
##### References:

show all references

##### References:
Intersections of $\theta(\tau)+2n\pi$ and $\tau \omega_+(\tau), \ n = 0,1,\dots,$. Values of $\tau$ at which the characteristic equation has pure imaginary eigenvalues, and hence candidates for critical values of $\tau$ at which there could be Hopf bifurcations. In both graphs, at all such intersections, transversality holds, since the slope of these curves at these intersections are different. Parameters: $m = 1,\ r = 1,\ K = 1, \ Y = 0.6 .$ (LEFT) $s = 0.02$. For $n = 0$ there are two intersections (i.e. $j_0 = 2$), at $\tau_0^1$ and $\tau_0^2$, but for $n = 1$, and hence $n\geqslant 1$, there are no intersections. (RIGHT) $s = 0.007$. There are two intersections each (i.e. $j_n = 2, \ n = 0,1,2$), at $\tau_n^1$ and $\tau_n^2$, for $n = 0,1$ and $2$, but for $n = 3$, and hence $n\geqslant 3$, there are no intersections. In both (LEFT) and (RIGHT), $E_+$ is asymptotically stable for $\tau\in[0,\tau_0^1)\cup(\tau_0^2,\tau_c)$ and unstable for $\tau\in(\tau_0^1,\tau_0^2)$.
Orbit diagrams. Initial data was taken to be $x(t) = y(t) = 0.1$ for $t\in[-\tau,0]$. However, we found bistability in the portion of the diagram beween the vertical dots and varied the initial data as explained below. Except for the portion between the vertical dots, the rest of the diagram was the same for all of the initial conditions we tried (not shown). (TOP) All local maxima and minima for the $y(t)$ coordinate of the attractor as $\tau$ varies, including kinks. (BOTTOM) Diagram including local maxes and mins for the $y(t)$ coordinate as $\tau$ varies, but with kinks eliminated. There are two saddle-node of limit cycle bifurcations. They occur for $\tau$ approximately equal to 76 and 82, where the curves in the orbit diagrams stop abruptly and there appear to be vertical dots. For $\tau$ between these values, there is is bistability. Two orbitally asymptotically stable periodic orbits (with their maximum and minimum amplitudes shown) and an unstable periodic orbit with amplitudes between them (not shown). The two stable periodic orbits were found by producing this part of the orbit diagram varying $\tau$ forward and then varying it backwards but startng at the last point of the attractor for the previous value of $\tau$
The time series for $y(t)$ when $\tau = 70$, depicting kinks. There are two local maxima and two local minima over each period as shown in Figure 2 (TOP), but only one local maxima and one local minima in Figure 2 (BOTTOM) in which kinks have been removed
Zoom-in of orbit diagram shown in Figure 2 for $\tau\in[75,100]$ including kinks. The vertical dots indicate the boundary of the region of bistability, where the two saddle-node of limit cycle bifurcations occur
Time delay embedding of two orbitally asymptotically stable periodic orbits demonstrating bistability for $\tau = 81$. The one with larger amplitude (dashed) has initial data $x(t) = y(t) = 0.1$, for $t \in [-\tau,0]$, and period approximately $345$. The one with smaller amplitude (solid) has initial data $x(t) = y(t) = 0.1$, for $t \in [-\tau,0)$ and $x(0) = 0.3, \ y(0) = 0.83$ and has period approximately $273.7$
(LEFT) Time series starting from the initial data $x(t) = y(t) = 0.1, t\in[-\tau,0]$ indicating how quickly the orbit gets close to the periodic attractor and (RIGHT) time delay embeddings of the periodic attractors, demonstrating the sequence of period doubling bifurcations initiating from the left at $\tau\approx 83,\ 86$, and $86.6$. Values of $\tau$ selected between these bifurcations: $\tau = 82,\ 85,\ 86.3$, and $86.8$, with periods of the periodic attractor approximately equal to: $340.2,\ 564.6,\ 1144.5$, and $2298.3$, respectively, are shown
(LEFT) Time series starting from the initial data $x(t) = y(t) = 0.1, t\in[-\tau,0]$ indicating how quickly the orbit gets close to the periodic attractor and (RIGHT) time delay embeddings of the periodic attractors, demonstrating the sequence of period halfing bifurcations initiating from the left for values of $\tau$ between $91.5$ and $92$, and at $\tau\approx 92.2, \ 93.2, \ 98.3$. Graphs shown are for values of $\tau$ between these bifurcations: $\tau = 91.95,\ 92,\ 93,\ 96$, and $100$, with periods approximately equal to: $4794.3, \ 2398.3, \ 1211.2,\ 557.6,$ and $295.3$, respectively
Time delay embedding starting at initial data $x(t) = y(t) = 0.1, t\in[\tau,0]$ for $\tau = 90.7$ showing a periodic attractor with period approximately $1800$, having $6(\neq2^n)$ loops for some integer $n$
(LEFT) Time series for $\tau = 90$ starting from the initial data $x(t) = y(t) = 0.1, t\in[-\tau,0]$. (RIGHT) Time delay embedding of the strange attractor for $\tau = 90$. Only the portion of the orbit from $t = 240,000-260,000$ is shown
The strange attractor, for $\tau = 90$, shown in Figure 9 in $(x,y)$-space. Only the portion of the orbit from $t = 240,000$ to $260,000$ is shown
The return map for $\tau = 90$, computing the minimum value of $y(t)$ as a function of the preceding minimum value of $y(t)$ for $y(t)<0.7$ in both cases, using the data in Figure 9
Time series (LEFT) for the solution that converges to a periodic attractor when $\tau = 92$, and (RIGHT) for the solution that converges to a strange attractor when $\tau = 90$, demonstrating that there is no sensitivity to initial data in the former case, but that there is sensitivity in the latter case. Initial data used for the solid curves: $x(t) = y(t) = 0.1$ for $t\in[-\tau,0]$, and for the dotted curves: $x(t) = 0.11$ and $y(t) = 0.1$ for $t\in[\tau,0]$
 [1] Yunfei Lv, Rong Yuan, Yuan He. Wavefronts of a stage structured model with state--dependent delay. Discrete & Continuous Dynamical Systems - A, 2015, 35 (10) : 4931-4954. doi: 10.3934/dcds.2015.35.4931 [2] Rui Hu, Yuan Yuan. Stability, bifurcation analysis in a neural network model with delay and diffusion. Conference Publications, 2009, 2009 (Special) : 367-376. doi: 10.3934/proc.2009.2009.367 [3] Guirong Jiang, Qishao Lu. The dynamics of a Prey-Predator model with impulsive state feedback control. Discrete & Continuous Dynamical Systems - B, 2006, 6 (6) : 1301-1320. doi: 10.3934/dcdsb.2006.6.1301 [4] Nhu N. Nguyen, George Yin. Stochastic partial differential equation models for spatially dependent predator-prey equations. Discrete & Continuous Dynamical Systems - B, 2020, 25 (1) : 117-139. doi: 10.3934/dcdsb.2019175 [5] Qigang Yuan, Jingli Ren. Periodic forcing on degenerate Hopf bifurcation. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2857-2877. doi: 10.3934/dcdsb.2020208 [6] Seung-Yeal Ha, Jinwook Jung, Jeongho Kim, Jinyeong Park, Xiongtao Zhang. A mean-field limit of the particle swarmalator model. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2021011 [7] Linlin Li, Bedreddine Ainseba. Large-time behavior of matured population in an age-structured model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2561-2580. doi: 10.3934/dcdsb.2020195 [8] Michael Grinfeld, Amy Novick-Cohen. Some remarks on stability for a phase field model with memory. Discrete & Continuous Dynamical Systems - A, 2006, 15 (4) : 1089-1117. doi: 10.3934/dcds.2006.15.1089 [9] Naeem M. H. Alkoumi, Pedro J. Torres. Estimates on the number of limit cycles of a generalized Abel equation. Discrete & Continuous Dynamical Systems - A, 2011, 31 (1) : 25-34. doi: 10.3934/dcds.2011.31.25 [10] Peter Benner, Jens Saak, M. Monir Uddin. Balancing based model reduction for structured index-2 unstable descriptor systems with application to flow control. Numerical Algebra, Control & Optimization, 2016, 6 (1) : 1-20. doi: 10.3934/naco.2016.6.1 [11] Jan Prüss, Laurent Pujo-Menjouet, G.F. Webb, Rico Zacher. Analysis of a model for the dynamics of prions. Discrete & Continuous Dynamical Systems - B, 2006, 6 (1) : 225-235. doi: 10.3934/dcdsb.2006.6.225 [12] Johannes Kellendonk, Lorenzo Sadun. Conjugacies of model sets. Discrete & Continuous Dynamical Systems - A, 2017, 37 (7) : 3805-3830. doi: 10.3934/dcds.2017161 [13] Marat Akhmet, Ejaily Milad Alejaily. Abstract similarity, fractals and chaos. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2479-2497. doi: 10.3934/dcdsb.2020191 [14] Didier Bresch, Thierry Colin, Emmanuel Grenier, Benjamin Ribba, Olivier Saut. A viscoelastic model for avascular tumor growth. Conference Publications, 2009, 2009 (Special) : 101-108. doi: 10.3934/proc.2009.2009.101 [15] Ondrej Budáč, Michael Herrmann, Barbara Niethammer, Andrej Spielmann. On a model for mass aggregation with maximal size. Kinetic & Related Models, 2011, 4 (2) : 427-439. doi: 10.3934/krm.2011.4.427 [16] Martin Bohner, Sabrina Streipert. Optimal harvesting policy for the Beverton--Holt model. Mathematical Biosciences & Engineering, 2016, 13 (4) : 673-695. doi: 10.3934/mbe.2016014 [17] Juan Manuel Pastor, Javier García-Algarra, Javier Galeano, José María Iriondo, José J. Ramasco. A simple and bounded model of population dynamics for mutualistic networks. Networks & Heterogeneous Media, 2015, 10 (1) : 53-70. doi: 10.3934/nhm.2015.10.53 [18] Chin-Chin Wu. Existence of traveling wavefront for discrete bistable competition model. Discrete & Continuous Dynamical Systems - B, 2011, 16 (3) : 973-984. doi: 10.3934/dcdsb.2011.16.973 [19] Alba Málaga Sabogal, Serge Troubetzkoy. Minimality of the Ehrenfest wind-tree model. Journal of Modern Dynamics, 2016, 10: 209-228. doi: 10.3934/jmd.2016.10.209 [20] Raghda A. M. Attia, Dumitru Baleanu, Dianchen Lu, Mostafa M. A. Khater, El-Sayed Ahmed. Computational and numerical simulations for the deoxyribonucleic acid (DNA) model. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021018

2019 Impact Factor: 1.27