Greenberg–Hastings cellular automaton

The Greenberg–Hastings Cellular Automaton (abbrev. GH model) is a three state two dimensional cellular automaton (abbrev CA) named after James M. Greenberg and Stuart Hastings,[1] designed to model excitable media,[2] One advantage of a CA model is ease of computation. The model can be understood quite well using simple "hand" calculations, not involving a computer.[2] Another advantage is that, at least in this case, one can prove a theorem characterizing those initial conditions which lead to repetitive behavior.[3]

Informal description

As in a typical two dimensional cellular automaton,[1] consider a rectangular grid, or checkerboard pattern, of "cells". It can be finite or infinite in extent. Each cell has a set of "neighbors". In the simplest case, each cell has four neighbors, those being the cells directly above or below or to the left or right of the given cell.[2]

Like this: the b's are all of the neighbors of the a. The a is one of the neighbors of each of the b's.

                                                  b
                                                b a b
                                                  b

At each "time" t=0,1,2,3,...., each cell is assigned one of three "states", typically called "resting" (or "quiescent" ; see excitable medium), "excited", or "refractory".[2] The assignment of states for all cells is arbitrary for t = 0, and then at subsequent times the state of each cell is determined by the following rules.[2]

1. If a cell is in the excited state at time t then it is in the refractory state at time t+1.

2. If a cell is in the refractory state at time t then it is in the resting state at time t+1.

3. If a given cell is in the resting state at time t and at least one of its neighbors is in the excited state at time t, then the given cell is in the excited state at time t+1, otherwise (no neighbor is excited at time t) it remains in the resting state at time t+1.

In this way the whole grid of cells advances from their initial states at t = 0 to their states at t = 1, then to their states at t = 2,3,4, etc., producing a pattern of cells in the various states for each time.

See the first animation in Belousov–Zhabotinsky reaction for a striking example of behavior that can be exhibited by this model. The three states are indicated by different colors.

Mathematical description

To describe the GH model more mathematically, consider the simplest case of a grid of square cells . At each time the cell has a state The type of neighborhood is not important, so long as each cell has some neighbors. In a square grid (as opposed to hexagonal), either four or eight cell neighborhoods work fine. The evolution rule is as follows: If then . If then If and for some neighboring cell then . Otherwise, [1][2]

Relation to a model of Wiener and Rosenblueth

While the GH model has often been compared with the ground breaking model of Wiener and Rosenblueth [4] developed earlier for the same purpose, the analogy is incorrect because the latter is not a CA. See for example,[5] in which it is stated that "The organization of muscle cells, the contraction of muscle, the dependence of the activity of the medium on the activity of its component elements, problems of memory, reliability, and mobility were formulated by Wiener in the form of definitions and theorems for a three-phase threshold-invariant continuous excitable medium." Even this statement is misleading, however. By reading the original paper [4] carefully, it can be seen that neither the time, the medium, nor the state are discrete. This is immediately obvious as regards the time and the medium. Wiener and Rosenblueth do say, though, that "There are three conditions in which any given region of the fiber can exist." They name these as the "active", the "refractory", and the "rest" states. However, in the next paragraph they refine this by identifying the rest state with the number 1, the active state with the number 0, and the refractory state with the open interval (0,1) on the real line. The number assigned to a give point at a given time is called its "epoch number" at that time. Thus, the "state space" is the interval [0,1]. And in the next paragraph we find "..behind every wave front moving freely there will be a band of fixed width within which the recovery process is taking place." A wave therefore consists of a "front", a smooth curve of points with epoch number 0 which moves in the plane with constant speed, followed by a refractory region of points with epoch number in (0,1), of constant width (depending on the velocity), and leaving behind a rest region of points with epoch number 1. This is far from a cellular automata, and is more correctly called a "geometric" model.

Further on in [5] it is stated that

"Automaton models of excitable media were investigated in [9] and [10]. These models are a discrete analog of Wiener's medium." (The models in their citations "[9]" and "[10]" are different from GH.)

Many writers have called the Weiner–Rosenblueth model a cellular automaton. The earliest paper found by Google Scholar with this designation clearly stated is.[6] However, as mentioned above, the continuity of the Wiener-Rosenblueth medium has not so far allowed as precise a theorem about persistence of patterns as the one for GH which is described below. On the other hand, several theorems are stated in [4] which are similar to those proved in,[3] though the proofs given for the theorems in [4] are less clear than those in [3] because of the natures of the respective models.

See also [7] for an often cited computer study based on a model which is similar to that of Wiener and Rosenblueth.

Generating a spiral

One interesting behavior is seen, for a four cell neighborhood and a square grid, when the initial condition consists of a half line of excited cells (1's) going off to infinity and below this half line a half line of refractory cells (2's). The rest of the cells are resting when t=0.[2]

Like this:

...................................................................
....00000000000000000000000000000000000000000000000000000000000....
....00000000000000000000000000000000000000000000000000000000000....
....00000000000000000001111111111111111111111111111111111111111....
....00000000000000000002222222222222222222222222222222222222222....
....00000000000000000000000000000000000000000000000000000000000....
....00000000000000000000000000000000000000000000000000000000000....
................................................................... 

The spiral which is produced can be seen in.[2] This spiral can easily be seen by following the pattern forward for a few iterations "by hand". No computer is necessary.

A theorem

As stated in,[2] and proved in [3] (which also considered n-state models), if the set of initial 1's is finite, then every individual cell oscillates forever through the cycle 0,1,2,0 if and only if at the start there is at least one square of four neighboring cells with one of the following patterns:

                    1 2       2 0       0 1
                    0 0       1 1       2 2

or some reflection or rotation of one of these. These patterns cannot be killed off by their surroundings. It is the "only if" part of this result which is most interesting. If none of these patterns is present at t = 0, then in any bounded region the pattern eventually settles down to all zeros.[2] The key tool in the proof in [3] is a "winding number" which is shown to be invariant for this model.

One easy consequence of the theorem stated above is that if there are no "refractory" cells initially then the pattern will die out in any bounded region (whether the total grid is finite or infinite in extent). This property of an excitable medium was found earlier, in the paper of Wiener and Rosenblueth,[4] and does not hold if there are "holes" in the region.

Notes

  1. ^ a b c G. de Vries; T. Hillen; M. Lewis; J. Miller; B. Schonfisch (2006). "6". A Course in Mathematical Biology: Quantitative Modeling with Mathematical & Computational Methods. SIAM.
  2. ^ a b c d e f g h i j J. M. Greenberg; S. P. Hastings (1978). "Spatial Patterns for Discrete Models of Diffusion in Excitable Media". SIAM Journal on Applied Mathematics. 54 (3): 515–523. doi:10.1137/0134040.
  3. ^ a b c d e C. Greene, James; J. M. Greenberg, Curtis; S. P. Hastings, Stuart (1980). "A combinatorial problem arising in the study of reaction-diffusion equations". SIAM J. Algebr. Discrete Methods. 1 (1): 34–42. doi:10.1137/0601006.
  4. ^ a b c d e N. Wiener; A. Rosenblueth (1946). "The mathematical formulation of the problem of conduction of connected excitable elements, specifically in cardiac muscle". Arch. Inst. Cardiol. Mex. 16 (3): 205–265. PMID 20245817.
  5. ^ a b Letichevskii, A. A.; Reshodko, L. V. (1974). "N. Wiener's theory of the activity of excitable media". Cybernetics. 8: 856–864. doi:10.1007/bf01068458.
  6. ^ M. Gerhardt; H. Schuster; J. Tyson (1990). "A cellular automaton model of excitable media: II. Curvature, dispersion, rotating waves and meandering waves". Physica D. 46: 392–415. doi:10.1016/0167-2789(90)90101-T.
  7. ^ G. K. Moe; W. C. Rheinboldt; J. A. Abildskov (1964). "A computer model of atrial fibrillation". Am. Heart J. 67: 200–220. doi:10.1016/0002-8703(64)90371-0.

References

  • R. Fisch, J. Gravner, D. Griffeath, Metastability in the Greenberg–Hastings model, The Annals of Applied Probability, vol. 3 (1993), 935–967.
  • R. Durrett and J. Steif, Some rigorous results for the Greenberg–Hastings model, Journal of Theoretical Probability vol 4 (1991), 669–690.
  • S. Wolfram, A New Kind of Science, 2003, pg. 1013.