[Paper Review] Second order additive invariants in elementary cellular automata

Henryk Fukś|ArXiv.org|Feb 17, 2005

Cellular Automata and Applications18 references3 citations

TL;DR

This paper investigates second-order additive invariants in elementary cellular automata (ECA), extending the study of conserved quantities beyond the standard number-conserving rules. Using Hattori and Takesue's general framework for additive invariants, it identifies ECA rules with second-order invariants, analyzes their fundamental diagrams (which exhibit singularities similar to first-order invariants), and finds that the current decays toward equilibrium with a power-law exponent of approximately -1/2 at critical densities—indicating universal critical dynamics across invariant orders.

ABSTRACT

We investigate second order additive invariants in elementary cellular automata rules. Fundamental diagrams of rules which possess additive invariants are either linear or exhibit singularities similar to singularities of rules with first-order invariant. Only rules which have exactly one invariants exhibit singularities. At the singularity, the current decays to its equilibrium value as a power law $t^α$, and the value of the exponent $α$ obtained from numerical simulations is very close to -1/2. This is in agreements with values previously reported for number-conserving rules, and leads to a conjecture that regardless of the order of the invariant, exponent $α$ seems to have a universal value of 1/2.

Motivation & Objective

• To extend the theory of additive invariants in cellular automata beyond first-order (number-conserving) rules to second-order invariants.
• To identify and classify elementary cellular automata rules that possess second-order additive invariants.
• To analyze the structure of fundamental diagrams in such rules, particularly the presence and nature of singularities.
• To investigate the dynamics at critical densities where the invariant's value equals its critical density, focusing on convergence behavior of the current.

Proposed method

• Adapts Hattori and Takesue's general existence condition for additive invariants to second-order invariants in binary, one-dimensional cellular automata.
• Applies the conservation law condition (equation 6) to derive necessary and sufficient conditions for a function ξ to be a second-order additive invariant.
• Uses the current function J defined via equation (7) to express the flux and verify conservation of the invariant across time steps.
• Employs numerical simulations to compute fundamental diagrams and analyze the time evolution of the current at critical densities.
• Compares dynamics of rules with second-order invariants to rule 184, particularly through spatiotemporal patterns and defect propagation.
• Identifies local transformations (e.g., via rule 60) that relate rules with second-order invariants to known rules like 184 and 226, enabling potential analytical extensions.

Experimental results

Research questions

• RQ1 Which elementary cellular automata rules possess second-order additive invariants, and how can they be systematically identified?
• RQ2 How do the fundamental diagrams of rules with second-order invariants differ from those with first-order invariants, particularly in terms of singularities?
• RQ3 What is the universal behavior of the current's relaxation toward equilibrium at critical density in rules with second-order invariants?
• RQ4 Can the dynamics of rules with second-order invariants be related to known models like rule 184 through local transformations?
• RQ5 Is the power-law decay exponent of the current at critical density universally -1/2 across different orders of additive invariants?

Key findings

• Five elementary cellular automata rules possess second-order additive invariants, with four of them (14, 35, 43, 142) exhibiting fundamental diagrams that are piecewise-linear with singularities.
• The singularities in the fundamental diagrams of these rules are structurally similar to those in first-order invariant rules, appearing only in rules with exactly one invariant.
• At critical density, the current decays toward equilibrium with a power-law behavior j(ρc, ∞) - j(ρc, t) ∼ t^(-1/2), indicating a universal exponent of -1/2 across different invariant orders.
• Rules 43 and 142 are locally transformable into rules 184 and 226 via rule 60, suggesting that their dynamics are closely related and may allow for rigorous derivation of the current's time evolution.
• Rules 14 and 35, though not locally transformable to 184, display spatiotemporal patterns with propagating and annihilating defects similar to rule 184, indicating analogous critical dynamics.