Stability of cycling behaviour near a heteroclinic network model of Rock-Paper-Scissors-Lizard-Spock

Claire M Postlethwaite,Alastair M Rucklidge
DOI: https://doi.org/10.1088/1361-6544/ac3560
2021-10-15
Abstract:The well-known game of Rock--Paper--Scissors can be used as a simple model of competition between three species. When modelled in continuous time using differential equations, the resulting system contains a heteroclinic cycle between the three equilibrium solutions representing the existence of only a single species. The game can be extended in a symmetric fashion by the addition of two further strategies (`Lizard' and `Spock'): now each strategy is dominant over two of the remaining four strategies, and is dominated by the remaining two. The differential equation model contains a set of coupled heteroclinic cycles forming a heteroclinic network. In this paper we carefully consider the dynamics near this heteroclinic network. We are able to identify regions of parameter space in which arbitrarily long periodic sequences of visits are made to the neighbourhoods of the equilibria, which form a complicated pattern in parameter space.
Dynamical Systems
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the stability analysis of dynamical behavior in the extended "Rock - Paper - Scissors - Lizard - Spock" game model. Specifically, the author focuses on the dynamical behavior near the heteroclinic network formed in this model, especially the stability regions of these behaviors in the parameter space. ### Background and Problem Description "Rock - Paper - Scissors" is a classic competition model that can be used to describe the competitive relationship among three species. When modeled using differential equations, there is a heteroclinic cycle in the system, connecting three equilibrium points, and each equilibrium point represents the state where only one species exists. The paper further extends this model by introducing two new strategies: "Lizard" and "Spock", so that each strategy is superior to two other strategies and is also克制 by two other strategies. This extension forms a more complex heteroclinic network. ### Research Objectives 1. **Stability Analysis of Dynamical Behavior**: The author has developed a technique to calculate in numerical continuation software using the previously defined stability concept (called fragmentary asymptotic stability). This enables them to identify that in certain parameter regions, the trajectory will visit the neighborhoods near the equilibrium points in an arbitrarily long periodic sequence, forming complex patterns. 2. **Stability of Complex Sequences**: In addition to analyzing the stability of sub - cycles, the author also studies the behavior of the trajectory approaching the heteroclinic network under specific sequences, which may include visiting the same equilibrium point multiple times and then visiting different equilibrium points. ### Main Contributions 1. **Computational Method**: Proposed an explicit method to calculate the stability of any complex (repeated) sequence in the heteroclinic network, that is, whether the trajectory will be attracted to or repelled from the network. 2. **Application of Numerical Continuation Software**: Adapted the stability conditions of the heteroclinic cycle into a scalar function that can be used in numerical continuation software, thereby calculating the boundaries of stability. 3. **Stability Regions in Parameter Space**: Discovered that there are complex stability regions in the parameter space, which are similar to the sausage - like resonant tongue - shaped regions in piecewise - continuous systems. These regions indicate that even when certain sufficient conditions are not met, the network still has strong attractive properties. ### Methods and Results - **Poincaré Map**: Constructed a Poincaré map to capture the behavior of the trajectory near the heteroclinic network, regardless of which path the trajectory takes in the network. - **Stability Analysis**: Used the Poincaré map and Podvigina's stability conditions to determine the stability regions of the trajectory approaching the heteroclinic network in different patterns. - **Numerical Simulation**: Verified the theoretical results through numerical simulation, showing the dynamical behaviors in different parameter regions, including periodic sequences, quasi - periodic solutions, and irregular behaviors. ### Conclusion The paper reveals the complex stability structure of the heteroclinic network in the parameter space through a detailed analysis of the dynamical behavior of the extended "Rock - Paper - Scissors - Lizard - Spock" model. These results not only enrich the understanding of heteroclinic network dynamics but also provide new tools and methods for research in related fields.