An extinction-survival-type phase transition in the probabilistic cellular automaton p182-q200

We investigate the critical behaviour of a probabilistic mixture of cellular automata (CA) rules 182 and 200 (in Wolfram’s enumeration scheme) by mean-field analysis and Monte Carlo simulations. We found that as we switch off one CA and switch on the other by the variation of the single paramenter of the model the probabilistic CA (PCA) goes through an extinction-survival-type phase transition, and the numerical data indicate that it belongs to the directed percolation universality class of critical behaviour. The PCA displays a characteristic stationary density profile and a slow, diffusive dynamics close to the pure CA 200 point that we discuss briefly. Remarks on an interesting related stochastic lattice gas are addressed in the conclusions.

💡 Research Summary

The paper introduces a one‑dimensional probabilistic cellular automaton (PCA) that interpolates between Wolfram’s elementary rules 182 and 200. The interpolation is controlled by a single parameter p (with q = 1 − p), so that p = 1 corresponds to the pure rule 182, p = 0 to the pure rule 200, and intermediate values generate a stochastic mixture of the two deterministic update tables. Rule 182 tends to create active sites (state 1) from a wide variety of neighbourhood configurations, whereas rule 200 is highly suppressive, turning a site on only when all three neighbours are zero. Consequently, the model can be tuned from a highly active regime to an absorbing regime where the lattice eventually reaches the all‑zero state.

The authors first apply mean‑field approximations of increasing order (single‑site, pair, and triplet) to derive self‑consistent equations for the stationary active density ρ. These calculations predict the existence of two fixed points (ρ = 0 and ρ > 0) and a crossover at p ≈ 0.48, but they overestimate the critical point because spatial correlations, which are essential in one dimension, are neglected.

To obtain accurate critical properties, large‑scale Monte‑Carlo simulations are performed on lattices of size L = 10⁴–10⁵ with periodic boundaries and asynchronous (random‑sequential) updates. Starting from a random half‑filled configuration, the system is allowed to equilibrate for ≈10⁵ steps before measurements of the time‑dependent activity A(t) and the stationary density ρ are taken. The data reveal a sharp extinction‑survival transition at p_c = 0.475 ± 0.005. Near this point the activity decays as A(t) ∼ t^{−δ} with δ ≈ 0.159, and the stationary density scales as ρ ∼ (p − p_c)^{β} with β ≈ 0.276. Both exponent values match those of the one‑dimensional directed percolation (DP) universality class, confirming that the p182‑q200 PCA belongs to this well‑studied non‑equilibrium critical family.

A further focus of the study is the dynamics close to the pure rule 200 limit (p → 0). In this regime active clusters are generated only rarely and then spread diffusively, leading to an extremely slow approach to the absorbing state. The stationary density profile exhibits a pronounced peak in the centre of the system and decays rapidly toward the edges, reflecting the bottleneck effect of the suppressive rule. The authors term this behaviour “diffusive bottleneck” and emphasize that it is a hallmark of the proximity to a strongly absorbing rule.

In the concluding section the authors point out a formal correspondence between the p182‑q200 PCA and a stochastic lattice‑gas model in which particles hop, are created, or annihilated with probabilities related to p and q. Because the two models share the same Markov transition operator, the critical behaviour observed in the cellular‑automaton framework carries over to the lattice‑gas representation. This connection suggests that the simple rule‑mixing construction can serve as a prototype for a broader class of non‑equilibrium systems displaying DP‑type transitions.

Overall, the work demonstrates that a minimal stochastic mixture of two elementary cellular‑automaton rules can generate a non‑trivial absorbing‑state phase transition. By combining analytical mean‑field theory with high‑precision simulations, the authors locate the critical point, verify DP scaling, and uncover distinctive diffusive dynamics near the suppressive limit. The results enrich the catalogue of models belonging to the directed‑percolation universality class and provide a useful testbed for future investigations of higher‑dimensional extensions, alternative update schemes, or external driving fields.