Biodiversity in Spatial Rock-Paper-Scissors Games

The struggle for survival between competing bacteria species in a petri dish can create beautiful spiral patterns—until one type conquers all the available space, destroying the entangled structure. Kerr et al. made these observations when they mixed three E. coli populations exhibiting cyclic dominance (Kerr et al., 2002): A beats B beats C beats A. As shown in (Reichenbach et al., 2007; 2008), the mobility of the bacteria, that is, their spatial relocation rate, is crucial for the stability of coexistence. Using this Demonstration, you can simulate a spatial rock-paper-scissors game and observe the role mobility plays.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Random cells on a square lattice with periodic boundary conditions are populated with three species, A (blue), B (yellow), and C (red) that cyclically dominate each other, as in a game of rock-paper-scissors. The simplex below the grid shows you the relative distribution of the species. When the simulation starts, three different reaction types occur at rates that you can adjust with the corresponding sliders:
1. In a reproduction process, a pair of adjacent cells is randomly chosen. If either cell is currently empty, it is populated with the species that occupies the second cell. If the second cell is empty as well, both are left unoccupied.
2. A selection process only occurs if the two cells of a random pair of nearest neighbors are occupied by different species. Then, the cell that is populated with the weaker species (in the sense of the rock-paper-scissors principle described above) is vacated.
3. Whenever a mobility process occurs, two random nearest-neighbor cells simply exchange their respective contents.
If the frequency of the last reaction is low enough, the cells occupied by each of the species, which continuously chase each other, form a spiral-like pattern slowly crawling across the grid. With increasing mobility rate, the size of the spirals grows and the structure becomes more and more coarse-grained, until it eventually dissolves completely as one species is able to conquer the entire grid. This disappearance of biodiversity is inevitable for any setting of the mobility slider. Beyond a critical mobility rate, however, the average duration of coexistence grows slower than linearly in the size of the grid and you will notice that a one-species state is reached much earlier.
The population frequencies of the three species are visualized as a point on the triangular simplex below the grid: close proximity to one of its vertices signals that a large share of the total population belongs to the species indicated next to the vertex. The fraction of grid cells occupied by any of the three species is shown by the slider next to the simplex.
You can see the finer details of the spirals more easily on larger grids, but the formation of any discernible structure takes longer as well. You can also examine the effects of choosing pairs of Moore neighbors (i.e., two cells sharing an edge or a corner) instead of the default von Neumann neighbors (i.e., two cells with a common edge), or try setting a temperature gradient so that reactions occur more often near the center of the grid.
Finally, you can adjust the probability distribution that determines the initial grid population and even include blocked cells, which do not take part in any interaction, to study how the game evolves in less regular topologies.
B. Kerr, M. A. Riley, M. W. Feldman, and B. J. M. Bohannan, "Local Dispersal Promotes Biodiversity in a Real-Life Game of Rock-Paper-Scissors," Nature, 418, 2002 pp. 171–174.
T. Reichenbach, M. Mobilia, and E. Frey, "Mobility Promotes and Jeopardizes Biodiversity in Rock-Paper-Scissors Games," Nature, 448, 2007 pp. 1046–1049.
T. Reichenbach, M. Mobilia, and E. Frey, "Self-Organization of Mobile Populations in Cyclic Competition," J. Theor. Biol., 254(2), 2008 pp. 368–383.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+