Cellular automaton models for time-correlated random walks: derivation and analysis.

Many diffusion processes in nature and society were found to be anomalous, in the sense of being fundamentally different from conventional Brownian motion. An important example is the migration of biological cells, which exhibits non-trivial temporal decay of velocity autocorrelation functions. This means that the corresponding dynamics is characterized by memory effects that slowly decay in time. Motivated by this we construct non-Markovian lattice-gas cellular automata models for moving agents with memory. For this purpose the reorientation probabilities are derived from velocity autocorrelation functions that are given a priori; in that respect our approach is "data-driven". Particular examples we consider are velocity correlations that decay exponentially or as power laws, where the latter functions generate anomalous diffusion. The computational efficiency of cellular automata combined with our analytical results paves the way to explore the relevance of memory and anomalous diffusion for the dynamics of interacting cell populations, like confluent cell monolayers and cell clustering.

Extracting cellular automaton rules from physical Langevin equation models for single and collective cell migration.

Cellular automata (CA) are discrete time, space, and state models which are extensively used for modeling biological phenomena. CA are "on-lattice" models with low computational demands. In particular, lattice-gas cellular automata (LGCA) have been introduced as models of single and collective cell migration. The interaction rule dictates the behavior of a cellular automaton model and is critical to the model's biological relevance. The LGCA model's interaction rule has been typically chosen phenomenologically. In this paper, we introduce a method to obtain lattice-gas cellular automaton interaction rules from physically-motivated "off-lattice" Langevin equation models for migrating cells. In particular, we consider Langevin equations related to single cell movement (movement of cells independent of each other) and collective cell migration (movement influenced by cell-cell interactions). As examples of collective cell migration, two different alignment mechanisms are studied: polar and nematic alignment. Both kinds of alignment have been observed in biological systems such as swarms of amoebae and myxobacteria. Polar alignment causes cells to align their velocities parallel to each other, whereas nematic alignment drives cells to align either parallel or antiparallel to each other. Under appropriate assumptions, we have derived the LGCA transition probability rule from the steady-state distribution of the off-lattice Fokker-Planck equation. Comparing alignment order parameters between the original Langevin model and the derived LGCA for both mechanisms, we found different areas of agreement in the parameter space. Finally, we discuss potential reasons for model disagreement and propose extensions to the CA rule derivation methodology.