Pattern formation, synchronization, and outbreak of biodiversity in cyclically competing games

Optimal dispersal and diffusion-enhanced robustness in two-patch metapopulations: origin's saddle-source nature matters

A two-patch logistic metapopulation model is investigated both analytically and numerically focusing on the impact of dispersal on population dynamics. First, the dependence of the global dynamics on the stability type of the full extinction equilibrium point is tackled. Then, the behaviour of the total population with respect to the dispersal is studied analytically. Our findings demonstrate that diffusion plays a crucial role in the preservation of both subpopulations and the full metapopulation under the presence of stochastic perturbations. At low diffusion, the origin is a repulsor, causing the orbits to flow nearly parallel to the axes, risking stochastic extinctions. Higher diffusion turns the repeller into a saddle point. Orbits then quickly converge to the saddle's unstable manifold, reducing extinction chances. This change in the vector field enhances metapopulation robustness. On the other hand, the well-known fact that asymmetric conditions on the patches is beneficial for the total population is further investigated. This phenomenon has been studied in previous works for large enough or small enough values of the dispersal. In this work, we complete the theory for all values of the dispersal. In particular, we derive analytically a formula for the optimal value of the dispersal that maximizes the total population.

Response of a three-species cyclic ecosystem to a short-lived elevation of death rate

A balanced ecosystem with coexisting constituent species is often perturbed by different natural events that persist only for a finite duration of time. What becomes important is whether, in the aftermath, the ecosystem recovers its balance or not. Here we study the fate of an ecosystem by monitoring the dynamics of a particular species that encounters a sudden increase in death rate. For exploration of the fate of the species, we use Monte-Carlo simulation on a three-species cyclic rock-paper-scissor model. The density of the affected (by perturbation) species is found to drop exponentially immediately after the pulse is applied. In spite of showing this exponential decay as a short-time behavior, there exists a region in parameter space where this species surprisingly remains as a single survivor, wiping out the other two which had not been directly affected by the perturbation. Numerical simulations using stochastic differential equations of the species give consistency to our results.

The complementary contribution of each order topology into the synchronization of multi-order networks

Higher-order interactions improve our capability to model real-world complex systems ranging from physics and neuroscience to economics and social sciences. There is great interest nowadays in understanding the contribution of higher-order terms to the collective behavior of the network. In this work, we investigate the stability of complete synchronization of complex networks with higher-order structures. We demonstrate that the synchronization level of a network composed of nodes interacting simultaneously via multiple orders is maintained regardless of the intensity of coupling strength across different orders. We articulate that lower-order and higher-order topologies work together complementarily to provide the optimal stable configuration, challenging previous conclusions that higher-order interactions promote the stability of synchronization. Furthermore, we find that simply adding higher-order interactions based on existing connections, as in simple complexes, does not have a significant impact on synchronization. The universal applicability of our work lies in the comprehensive analysis of different network topologies, including hypergraphs and simplicial complexes, and the utilization of appropriate rescaling to assess the impact of higher-order interactions on synchronization stability.

Enhancing coexistence of mobile species in the cyclic competition system by wildlife refuge

We investigate evolving dynamics of cyclically competing species on spatially extended systems with considering a specific region, which is called the "wildlife refuge," one of the institutional ways to preserve species biodiversity. Through Monte-Carlo simulations, we found that the refuge can play not groundbreaking but an important role in species survival. Species coexistence is maintained at a moderate mobility regime, which traditionally leads to the collapse of coexistence, and eventually, the extinction is postponed depending on the competition rate rather than the portion of the refuge. Incorporating the extinction probability and Fourier transform supported our results in both stochastic and analogous ways. Our findings may provide valuable evidence to assist fields of ecological/biological sciences in understanding the presence and construction of refuges for wildlife with associated effects on species biodiversity.

Mortality makes coexistence vulnerable in evolutionary game of rock-paper-scissors

Multiple species in the ecosystem are believed to compete cyclically for maintaining balance in nature. The evolutionary dynamics of cyclic interaction crucially depends on different interactions representing different natural habits. Based on a rock-paper-scissors model of cyclic competition, we explore the role of mortality of individual organisms in the collective survival of a species. For this purpose a parameter called "natural death" is introduced. It is meant for bringing about the decease of an individual irrespective of any intra- and interspecific interaction. We perform a Monte Carlo simulation followed by a stability analysis of different fixed points of defined rate equations and observe that the natural death rate is surprisingly one of the most significant factors in deciding whether an ecosystem would come up with a coexistence or a single-species survival.

Dynamically generated hierarchies in games of competition

Spatial many-species predator-prey systems have been shown to yield very rich space-time patterns. This observation begs the question whether there exist universal mechanisms for generating this type of emerging complex patterns in nonequilibrium systems. In this work we investigate the possibility of dynamically generated hierarchies in predator-prey systems. We analyze a nine-species model with competing interactions and show that the studied situation results in the spontaneous formation of spirals within spirals. The parameter dependence of these intriguing nested spirals is elucidated. This is achieved through the numerical investigation of various quantities (correlation lengths, densities of empty sites, Fourier analysis of species densities, interface fluctuations) that allows us to gain a rather complete understanding of the spatial arrangements and the temporal evolution of the system. A possible generalization of the interaction scheme yielding dynamically generated hierarchies is discussed. As cyclic interactions occur spontaneously in systems with competing strategies, the mechanism discussed in this work should contribute to our understanding of various social and biological systems.

Emergence of oscillatory coexistence with exponentially decayed waiting times in a coupled cyclic competition system

Interpatch migration between two environments is generally considered as a spatial concept and can affect species biodiversity in each patch by inducing flux of population such as inflow and outflow quantities of species. In this paper, we explore the effect of interpatch migration, which can be generally considered as a spatial concept and may affect species biodiversity between two different patches in the perspective of the macroscopic level by exploiting the coupling of two systems, where each patch is occupied by cyclically competing three species who can stably coexist by exhibiting periodic orbits. For two simple scenarios of interpatch migration either single or all species migration, we found that two systems with independently stable coexisting species in each patch are eventually synchronized, and oscillatory behaviors of species densities in two patches become identical, i.e., the synchronized coexistence emerges. In addition, we find that, whether single or all species interpatch migration occurs, the waiting time for the synchronization is exponentially decreasing as the coupling strength is intensified. Our findings suggest that the synchronized behavior of species as a result of migration between different patches can be easily predicted by the coupling of systems and additional information such as waiting times and sensitivity of initial densities.

Robust coexistence with alternative competition strategy in the spatial cyclic game of five species

Alternative strategy is common in animal populations to promote reproductive fitness by obtaining resources. In spatial dynamics of cyclic competition, reproduction can occur when individuals obtain vacant rooms and, in this regard, empty sites should be resources for reproduction which can be induced by interspecific competition. In this paper, we study the role of alternative competition in the spatial system of cyclically competing five species by utilizing rock-paper-scissors-lizard-spock game. From Monte-Carlo simulations, we found that strong alternative competition can lead to the reemergence of coexistence of five species regardless of mobility, which is never reported in previous works under the symmetric competition structure. By investigating the coexistence probability, we also found that coexistence alternates by passing certain degrees of alternative competition in combination with mobility. In addition, we provided evidences in the opposite scenario by strengthening spontaneous competition, which exhibits the reemergence of coexistence similarly. Our findings may suggest more comprehensive perspectives to interpret mechanisms for biodiversity by alternative strategies in spatially extended systems than previously reported.

Invasion-controlled pattern formation in a generalized multispecies predator-prey system

Rock-scissors-paper game, as the simplest model of intransitive relation between competing agents, is a frequently quoted model to explain the stable diversity of competitors in the race of surviving. When increasing the number of competitors we may face a novel situation because beside the mentioned unidirectional predator-prey-like dominance a balanced or peer relation can emerge between some competitors. By utilizing this possibility in the present work we generalize a four-state predator-prey-type model where we establish two groups of species labeled by even and odd numbers. In particular, we introduce different invasion probabilities between and within these groups, which results in a tunable intensity of bidirectional invasion among peer species. Our study reveals an exceptional richness of pattern formations where five quantitatively different phases are observed by varying solely the strength of the mentioned inner invasion. The related transition points can be identified with the help of appropriate order parameters based on the spatial autocorrelation decay, on the fraction of empty sites, and on the variance of the species density. Furthermore, the application of diverse, alliance-specific inner invasion rates for different groups may result in the extinction of the pair of species where this inner invasion is moderate. These observations highlight that beyond the well-known and intensively studied cyclic dominance there is an additional source of complexity of pattern formation that has not been explored earlier.

Nonlinear dynamics with Hopf bifurcations by targeted mutation in the system of rock-paper-scissors metaphor

The role of mutation, which is an error process in gene evolution, in systems of cyclically competing species has been studied from various perspectives, and it is regarded as one of the key factors for promoting coexistence of all species. In addition to naturally occurring mutations, many experiments in genetic engineering have involved targeted mutation techniques such as recombination between DNA and somatic cell sequences and have studied genetic modifications through loss or augmentation of cell functions. In this paper, we investigate nonlinear dynamics with targeted mutation in cyclically competing species. In different ways to classic approaches of mutation in cyclic games, we assume that mutation may occur in targeted individuals who have been removed from intraspecific competition. By investigating each scenario depending on the number of objects for targeted mutation analytically and numerically, we found that targeted mutation can lead to persistent coexistence of all species. In addition, under the specific condition of targeted mutation, we found that targeted mutation can lead to emergences of bistable states for species survival. Through the linear stability analysis of rate equations, we found that those phenomena are accompanied by Hopf bifurcation which is supercritical. Our findings may provide more global perspectives on understanding underlying mechanisms to control biodiversity in ecological/biological sciences, and evidences with mathematical foundations to resolve social dilemmas such as a turnover of group members by resigning with intragroup conflicts in social sciences.

Biodiversity in the cyclic competition system of three species according to the emergence of mutant species

Understanding mechanisms which promote or hinder existing ecosystems are important issues in ecological sciences. In addition to fundamental interactions such as competition and migration among native species, existing ecosystems can be easily disturbed by external factors, and the emergence of new species may be an example in such cases. The new species which does not exist in a current ecosystem can be regarded as either alien species entered from outside or mutant species born by mutation in existing normal species. Recently, as existing ecosystems are getting influenced by various physical/chemical external factors, mutation due to anthropogenic and environmental factors can occur more frequently and is thus attracting much attention for the maintenance of ecosystems. In this paper, we consider emergences of mutant species among self-competing three species in the cyclic dominance. By defining mutation as the birth of mutant species, we investigate how mutant species can affect biodiversity in the existing ecosystem. Through microscopic and macroscopic approaches, we have found that the society of existing normal species can be disturbed by mutant species either the society is maintained accompanying with the coexistence of all species or jeopardized by occupying of mutant species. Due to the birth of mutant species, the existing society may be more complex by constituting two different groups of normal and mutant species, and our results can be contributed to analyze complex ecosystems of many species. We hope our findings may propose a new insight on mutation in cyclic competition systems of many species.

Cooperation enhanced by the coevolution of teaching activity in evolutionary prisoner's dilemma games with voluntary participation

Voluntary participation, as an additional strategy involved in repeated games, has been proved to be an efficient way to promote the evolution of cooperation theoretically and empirically. Besides, current studies show that the coevolution of teaching activity can promote cooperation. Thus, inspired by aforementioned above, we investigate the effect of coevolution of teaching activity on the evolution of cooperation for prisoner's dilemma game with voluntary participation: when the focal player successfully enforces its strategy on the opponent, his teaching ability will get an increase. Through numerical simulation, we have shown that voluntary participation could effectively promote the fraction of cooperation, which is also affected by the value of increment. Furthermore, we investigate the influence of the increment value on the density of different strategies and find that there exists an optimal increment value that plays an utmost role on the evolutionary dynamics. With regard to this observation, we unveil that an optimal value of increment can lead to strongest heterogeneity in agents' teaching ability, further promoting the evolution of cooperation.

Basins of distinct asymptotic states in the cyclically competing mobile five species game

We study the dynamics of cyclic competing mobile five species on spatially extended systems originated from asymmetric initial populations and investigate the basins for the three possible asymptotic states, coexistence of all species, existences of only two independent species, and the extinction. Through extensive numerical simulations, we find a prosperous dependence on initial conditions for species biodiversity. In particular, for fixed given equal densities of two relevant species, we find that only five basins for the existence of two independent species exist and they are spirally entangled for high mobility. A basin of coexistence is outbreaking when the mobility parameter is decreased through a critical value and surrounded by the other five basins. For fixed given equal densities of two independent species, however, we find that basin structures are not spirally entangled. Further, final states of two independent species are totally different. For all possible considerations, the extinction state is not witnessed which is verified by the survival probability. To provide the validity of basin structures from lattice simulations, we analyze the system in mean-field manners. Consequently, results on macroscopic levels are matched to direct lattice simulations for high mobility regimes. These findings provide a good insight into the fundamental issue of the biodiversity among many species than previous cases.

Emergence of unusual coexistence states in cyclic game systems

Evolutionary games of cyclic competitions have been extensively studied to gain insights into one of the most fundamental phenomena in nature: biodiversity that seems to be excluded by the principle of natural selection. The Rock-Paper-Scissors (RPS) game of three species and its extensions [e.g., the Rock-Paper-Scissors-Lizard-Spock (RPSLS) game] are paradigmatic models in this field. In all previous studies, the intrinsic symmetry associated with cyclic competitions imposes a limitation on the resulting coexistence states, leading to only selective types of such states. We investigate the effect of nonuniform intraspecific competitions on coexistence and find that a wider spectrum of coexistence states can emerge and persist. This surprising finding is substantiated using three classes of cyclic game models through stability analysis, Monte Carlo simulations and continuous spatiotemporal dynamical evolution from partial differential equations. Our finding indicates that intraspecific competitions or alternative symmetry-breaking mechanisms can promote biodiversity to a broader extent than previously thought.

Biodiversity in models of cyclic dominance is preserved by heterogeneity in site-specific invasion rates

Global, population-wide oscillations in models of cyclic dominance may result in the collapse of biodiversity due to the accidental extinction of one species in the loop. Previous research has shown that such oscillations can emerge if the interaction network has small-world properties, and more generally, because of long-range interactions among individuals or because of mobility. But although these features are all common in nature, global oscillations are rarely observed in actual biological systems. This begets the question what is the missing ingredient that would prevent local oscillations to synchronize across the population to form global oscillations. Here we show that, although heterogeneous species-specific invasion rates fail to have a noticeable impact on species coexistence, randomness in site-specific invasion rates successfully hinders the emergence of global oscillations and thus preserves biodiversity. Our model takes into account that the environment is often not uniform but rather spatially heterogeneous, which may influence the success of microscopic dynamics locally. This prevents the synchronization of locally emerging oscillations, and ultimately results in a phenomenon where one type of randomness is used to mitigate the adverse effects of other types of randomness in the system.

Extension of a spatial evolutionary coordination game with neutral options

The multiagent evolutionary games on a lattice are equivalent to a kinetic Ising model if the uniform pair interactions are defined by a two-strategy coordination game and the logit rule controls the strategy updates. Now we extend this model by allowing the players to use additional neutral strategies that provide zero payoffs for both players if one of them selects one of the neutral strategies. In the resulting n-strategy evolutionary games the analytical methods and numerical simulations indicate continuous order-disorder phase transitions when increasing the noise level if n does not exceed a threshold value. For larger n the system exhibits a first order phase transition at a critical noise level decreasing asymptotically as 2/ln(n).

Benefits of tolerance in public goods games

Leaving the joint enterprise when defection is unveiled is always a viable option to avoid being exploited. Although loner strategy helps the population not to be trapped into the tragedy of the commons state, it could offer only a modest income for nonparticipants. In this paper we demonstrate that showing some tolerance toward defectors could not only save cooperation in harsh environments but in fact results in a surprisingly high average payoff for group members in public goods games. Phase diagrams and the underlying spatial patterns reveal the high complexity of evolving states where cyclic dominant strategies or two-strategy alliances can characterize the final state of evolution. We identify microscopic mechanisms which are responsible for the superiority of global solutions containing tolerant players. This phenomenon is robust and can be observed both in well-mixed and in structured populations highlighting the importance of tolerance in our everyday life.

A double-edged sword: Benefits and pitfalls of heterogeneous punishment in evolutionary inspection games

As a simple model for criminal behavior, the traditional two-strategy inspection game yields counterintuitive results that fail to describe empirical data. The latter shows that crime is often recurrent, and that crime rates do not respond linearly to mitigation attempts. A more apt model entails ordinary people who neither commit nor sanction crime as the third strategy besides the criminals and punishers. Since ordinary people free-ride on the sanctioning efforts of punishers, they may introduce cyclic dominance that enables the coexistence of all three competing strategies. In this setup ordinary individuals become the biggest impediment to crime abatement. We therefore also consider heterogeneous punisher strategies, which seek to reduce their investment into fighting crime in order to attain a more competitive payoff. We show that this diversity of punishment leads to an explosion of complexity in the system, where the benefits and pitfalls of criminal behavior are revealed in the most unexpected ways. Due to the raise and fall of different alliances no less than six consecutive phase transitions occur in dependence on solely the temptation to succumb to criminal behavior, leading the population from ordinary people-dominated across punisher-dominated to crime-dominated phases, yet always failing to abolish crime completely.

Synchronization and extinction in cyclic games with mixed strategies

We consider cyclic Lotka-Volterra models with three and four strategies where at every interaction agents play a strategy using a time-dependent probability distribution. Agents learn from a loss by reducing the probability to play a losing strategy at the next interaction. For that, an agent is described as an urn containing β balls of three and four types, respectively, where after a loss one of the balls corresponding to the losing strategy is replaced by a ball representing the winning strategy. Using both mean-field rate equations and numerical simulations, we investigate a range of quantities that allows us to characterize the properties of these cyclic models with time-dependent probability distributions. For the three-strategy case in a spatial setting we observe a transition from neutrally stable to stable when changing the level of discretization of the probability distribution. For large values of β, yielding a good approximation to a continuous distribution, spatially synchronized temporal oscillations dominate the system. For the four-strategy game the system is always neutrally stable, but different regimes emerge, depending on the size of the system and the level of discretization.

How turbulence regulates biodiversity in systems with cyclic competition

Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.

Mesoscopic interactions and species coexistence in evolutionary game dynamics of cyclic competitions

Evolutionary dynamical models for cyclic competitions of three species (e.g., rock, paper, and scissors, or RPS) provide a paradigm, at the microscopic level of individual interactions, to address many issues in coexistence and biodiversity. Real ecosystems often involve competitions among more than three species. By extending the RPS game model to five (rock-paper-scissors-lizard-Spock, or RPSLS) mobile species, we uncover a fundamental type of mesoscopic interactions among subgroups of species. In particular, competitions at the microscopic level lead to the emergence of various local groups in different regions of the space, each involving three species. It is the interactions among the groups that fundamentally determine how many species can coexist. In fact, as the mobility is increased from zero, two transitions can occur: one from a five- to a three-species coexistence state and another from the latter to a uniform, single-species state. We develop a mean-field theory to show that, in order to understand the first transition, group interactions at the mesoscopic scale must be taken into account. Our findings suggest, more broadly, the importance of mesoscopic interactions in coexistence of great many species.

Cyclic dominance in evolutionary games: a review

Rock is wrapped by paper, paper is cut by scissors and scissors are crushed by rock. This simple game is popular among children and adults to decide on trivial disputes that have no obvious winner, but cyclic dominance is also at the heart of predator-prey interactions, the mating strategy of side-blotched lizards, the overgrowth of marine sessile organisms and competition in microbial populations. Cyclical interactions also emerge spontaneously in evolutionary games entailing volunteering, reward, punishment, and in fact are common when the competing strategies are three or more, regardless of the particularities of the game. Here, we review recent advances on the rock-paper-scissors (RPS) and related evolutionary games, focusing, in particular, on pattern formation, the impact of mobility and the spontaneous emergence of cyclic dominance. We also review mean-field and zero-dimensional RPS models and the application of the complex Ginzburg-Landau equation, and we highlight the importance and usefulness of statistical physics for the successful study of large-scale ecological systems. Directions for future research, related, for example, to dynamical effects of coevolutionary rules and invasion reversals owing to multi-point interactions, are also outlined.

From pairwise to group interactions in games of cyclic dominance

We study the rock-paper-scissors game in structured populations, where the invasion rates determine individual payoffs that govern the process of strategy change. The traditional version of the game is recovered if the payoffs for each potential invasion stem from a single pairwise interaction. However, the transformation of invasion rates to payoffs also allows the usage of larger interaction ranges. In addition to the traditional pairwise interaction, we therefore consider simultaneous interactions with all nearest neighbors, as well as with all nearest and next-nearest neighbors, thus effectively going from single pair to group interactions in games of cyclic dominance. We show that differences in the interaction range affect not only the stationary fractions of strategies but also their relations of dominance. The transition from pairwise to group interactions can thus decelerate and even revert the direction of the invasion between the competing strategies. Like in evolutionary social dilemmas, in games of cyclic dominance, too, the indirect multipoint interactions that are due to group interactions hence play a pivotal role. Our results indicate that, in addition to the invasion rates, the interaction range is at least as important for the maintenance of biodiversity among cyclically competing strategies.

Interfaces with internal structures in generalized rock-paper-scissors models

In this work we investigate the development of stable dynamical structures along interfaces separating domains belonging to enemy partnerships in the context of cyclic predator-prey models with an even number of species N≥8. We use both stochastic and field theory simulations in one and two spatial dimensions, as well as analytical arguments, to describe the association at the interfaces of mutually neutral individuals belonging to enemy partnerships and to probe their role in the development of the dynamical structures at the interfaces. We identify an interesting behavior associated with the symmetric or asymmetric evolution of the interface profiles depending on whether N/2 is odd or even, respectively. We also show that the macroscopic evolution of the interface network is not very sensitive to the internal structure of the interfaces. Although this work focuses on cyclic predator-prey models with an even number of species, we argue that the results are expected to be quite generic in the context of spatial stochastic May-Leonard models.

Hopf bifurcation and multistability in a system of phase oscillators

We study the phase reduction of two coupled van der Pol oscillators with asymmetric repulsive coupling under an external harmonic force. We show that the system of two phase oscillators undergoes a Hopf bifurcation and possesses multistability on a 2π-periodic phase plane. We describe the bifurcation mechanisms of formation of multistability in the phase-reduced system and show that the Andronov-Hopf bifurcation in the phase-reduced system is not an artifact of the reduction approach but, indeed, has its prototype in the nonreduced system. The bifurcational mechanisms presented in the paper enable one to describe synchronization effects in a wide class of interacting systems with repulsive coupling e.g., genetic oscillators.

Diverging fluctuations in a spatial five-species cyclic dominance game

A five-species predator-prey model is studied on a square lattice where each species has two prey and two predators on the analogy to the rock-paper-scissors-lizard-Spock game. The evolution of the spatial distribution of species is governed by site exchange and invasion between the neighboring predator-prey pairs, where the cyclic symmetry can be characterized by two different invasion rates. The mean-field analysis has indicated periodic oscillations in the species densities with a frequency becoming zero for a specific ratio of invasion rates. When varying the ratio of invasion rates, the appearance of this zero-eigenvalue mode is accompanied by neutrality between the species associations. Monte Carlo simulations of the spatial system reveal diverging fluctuations at a specific invasion rate, which can be related to the vanishing dominance between all pairs of species associations.

Persistent coexistence of cyclically competing species in spatially extended ecosystems

A fundamental result in the evolutionary-game paradigm of cyclic competition in spatially extended ecological systems, as represented by the classic Reichenbach-Mobilia-Frey (RMF) model, is that high mobility tends to hamper or even exclude species coexistence. This result was obtained under the hypothesis that individuals move randomly without taking into account the suitability of their local environment. We incorporate local habitat suitability into the RMF model and investigate its effect on coexistence. In particular, we hypothesize the use of "basic instinct" of an individual to determine its movement at any time step. That is, an individual is more likely to move when the local habitat becomes hostile and is no longer favorable for survival and growth. We show that, when such local habitat suitability is taken into account, robust coexistence can emerge even in the high-mobility regime where extinction is certain in the RMF model. A surprising finding is that coexistence is accompanied by the occurrence of substantial empty space in the system. Reexamination of the RMF model confirms the necessity and the important role of empty space in coexistence. Our study implies that adaptation/movements according to local habitat suitability are a fundamental factor to promote species coexistence and, consequently, biodiversity.

von Neummann's and related scaling laws in rock-paper-scissors-type games

We introduce a family of rock-paper-scissors-type models with Z(N) symmetry (N is the number of species), and we show that it has a very rich structure with many completely different phases. We study realizations that lead to the formation of domains, where individuals of one or more species coexist, separated by interfaces whose (average) dynamics is curvature driven. This type of behavior, which might be relevant for the development of biological complexity, leads to an interface network evolution and pattern formation similar to the ones of several other nonlinear systems in condensed matter and cosmology.

Junctions and spiral patterns in generalized rock-paper-scissors models

We investigate the population dynamics in generalized rock-paper-scissors models with an arbitrary number of species N. We show that spiral patterns with N arms may develop both for odd and even N, in particular in models where a bidirectional predation interaction of equal strength between all species is modified to include one N-cyclic predator-prey rule. While the former case gives rise to an interface network with Y-type junctions obeying the scaling law L∝t1/2, where L is the characteristic length of the network and t is the time, the latter can lead to a population network with N-armed spiral patterns, having a roughly constant characteristic length scale. We explicitly demonstrate the connection between interface junctions and spiral patterns in these models and compute the corresponding scaling laws. This work significantly extends the results of previous studies of population dynamics and could have profound implications for the understanding of biological complexity in systems with a large number of species.

Synchronization of multi-frequency noise-induced oscillations

Using a model system of FitzHugh-Nagumo type in the excitable regime, the similarity between synchronization of self-sustained and noise-induced oscillations is studied for the case of more than one main frequency in the spectrum. It is shown that this excitable system undergoes the same frequency lockings as a self-sustained quasiperiodic oscillator. The presence of noise-induced both stable and unstable limit cycles and tori, as well as their tangential bifurcations, are discussed. As the FitzHugh-Nagumo oscillator represents one of the basic neural models, the obtained results are of high importance for neuroscience.

Effects of competition on pattern formation in the rock-paper-scissors game

We investigate the impact of cyclic competition on pattern formation in the rock-paper-scissors game. By separately considering random and prepared initial conditions, we observe a critical influence of the competition rate p on the stability of spiral waves and on the emergence of biodiversity. In particular, while increasing values of p promote biodiversity, they may act detrimentally on spatial pattern formation. For random initial conditions, we observe a phase transition from biodiversity to an absorbing phase, whereby the critical value of mobility grows linearly with increasing values of p on a log-log scale but then saturates as p becomes large. For prepared initial conditions, we observe the formation of single-armed spirals, but only for values of p that are below a critical value. Once above that value, the spirals break up and form disordered spatial structures, mainly because of the percolation of vacant sites. Thus there exists a critical value of the competition rates p(c) for stable single-armed spirals in finite populations. Importantly though, p(c) increases with increasing system size because noise reinforces the disintegration of ordered patterns. In addition, we also find that p(c) increases with the mobility. These phenomena are reproduced by a deterministic model that is based on nonlinear partial differential equations. Our findings indicate that competition is vital for the sustenance of biodiversity and the emergence of pattern formation in ecosystems governed by cyclical interactions.