Threshold effect with stochastic fluctuation in bacteria-colony-like proliferation dynamics as analyzed through a comparative study of reaction-diffusion equations and cellular automata.

We report a comparative study on pattern formation between the methods of cellular automata (CA) and reaction-diffusion equations (RD) applying to a morphology of bacterial colony formation. To do so, we began the study with setting an extremely simple model, which was designed to realize autocatalytic proliferation of bacteria (denoted as X ) fed with nutrition (N) and their inactive state (prespore state) P1 due to starvation: X+N-->2X and X-->P1, respectively. It was found numerically that while the CA could successfully generate rich patterns ranging from the circular fat structure to the viscous-finger-like complicated one, the naive RD reproduced only the circular pattern but failed to give a finger structure. Augmenting the RD equations by adding two physical factors, (i) a threshold effect in the dynamics of X+N-->2X (breaking the continuity limit of RD) and (ii) internal noise with onset threshold (breaking the inherent symmetry of RD), we have found that the viscous-finger-like realistic patterns are indeed recovered by thus modified RD. This highlights the important difference between CA and RD, and at the same time, clarifies the necessary factors for the complicated patterns to emerge in such a surprisingly simple model system.