Structure and dynamics in the low-density phase of a two-dimensional cellular automaton model of traffic flow

We analyze the structure and dynamics in the low-density phase of the deterministic two-dimensional cellular automaton model of traffic flow introduced in [O. Biham, A.A. Middleton and D. Levine, Phys. Rev. A 46, R6124 (1992)]. The model consists of horizontally-oriented (H) cars that move to the right and vertically-oriented (V) cars that move downward, on a square lattice of size $L$ with periodic boundary conditions. Starting from a random initial state of density $p$, which is equally divided between the H and V-cars, the model exhibits a phase transition at a critical density $p_c$. For $p<p_c$ it evolves toward a free-flowing periodic (FFP) state, while for $p>p_c$ it evolves toward a fully-jammed state or to an intermediate state of congested traffic. In the FFP states, the H and V-cars segregate into homogeneous diagonal bands, in which they move freely without obstruction. To analyze the convergence toward the FFP states we introduce a configuration-space distance measure $D(t)=D_{\parallel}(t)+D_{\perp}(t)$ between the state of the system at time $t$ and the set of FFP states. The $D_{\parallel}(t)$ term accounts for the interactions between homotypic pairs of H (or V) cars, while $D_{\perp}(t)$ accounts for the interactions between heterotypic pairs of H and V-cars. We show that in the FFP states $D(t)=0$, while in all the other states $D(t)>0$. As the system evolves toward the FFP states, there is a separation of time scales, where $D_{\parallel}(t)$ decays very fast while $D_{\perp}(t)$ decays much more slowly. Moreover, the time dependence of $D_{\perp}(t)$ is well fitted by an exponentially truncated power-law decay of the form $D_{\perp}(t)\sim t^{-γ} \exp(-t/τ_{\perp})$, where $τ_{\perp}$ depends on $L$ and $p$. The power-law decay suggests avalanche-like dynamics with no characteristic scale, while the exponential cutoff is imposed by the finite lattice size.

💡 Research Summary

The paper presents a thorough investigation of the low‑density regime of the deterministic two‑dimensional Biham‑Middleton‑Levine (BML) cellular automaton traffic model. In the BML model a square lattice of size L×L hosts two species of cars: horizontally oriented (H) cars that move one cell to the right on odd time steps, and vertically oriented (V) cars that move one cell downward on even time steps. Periodic boundary conditions turn the lattice into a torus, and the total number of cars of each type, as well as the number of H‑cars per row and V‑cars per column, are strictly conserved. The initial configuration is random with overall density p, split equally between H and V (probability p/2 for each species per cell).

When the density lies below a critical value pc, the system invariably evolves toward a free‑flow periodic (FFP) state. In an FFP state the lattice self‑organizes into homogeneous diagonal bands: each counter‑diagonal (i + j mod L) contains either only H‑cars, only V‑cars, or is empty. Consequently there are no adjacent H‑cars in the same row, no adjacent V‑cars in the same column, and no heterotypic pairs on neighboring diagonals that could block motion. The dynamics then becomes perfectly periodic with period T = 2L; after each full cycle every car returns to its original lattice site and moves with unit speed (v = 1). These states are absorbing in the sense that once entered the deterministic dynamics cannot leave them.

To quantify how the system approaches an FFP state the authors introduce a configuration‑space distance \