Self-Orthogonal Cellular Automata
Luca Mariot, Federico Mazzone
TL;DR
The paper investigates when the Latin square generated by a bipermutive CA is self-orthogonal, introducing an algebraic criterion based on the invertibility of a stacked transition matrix. By exploiting the circulant structure of this matrix, it reduces self-orthogonality to gcd tests between the local-rule polynomial $p_f(X)$ and $X^{2(d-1)}-1$, with a binary-specific simplification to gcd$(p_f(X),X^{d-1}+1)=1$; irreducibility of $p_f$ suffices in the binary case, and a parity criterion applies when the stacked size is a power of two. Empirical exploration up to diameter $d=6$ shows no nonlinear self-orthogonal CA and highlights the prominence of linear SOCA, while the theoretical results provide a complete linear-characterization and a clear pathway to further questions. The findings link CA-based Latin squares to circulant-matrix invertibility and gcd computations, with potential implications for combinatorial designs, secret sharing, and cryptographic constructions.
Abstract
It is known that no-boundary Cellular Automata (CA) defined by bipermutive local rules give rise to Latin squares. In this paper, we study under which conditions the Latin square generated by a bipermutive CA is self-orthogonal, i.e. orthogonal to its transpose. We first enumerate all bipermutive CA over the binary alphabet up to diameter $d=6$, remarking that only some linear rules give rise to self-orthogonal Latin squares. We then give a full theoretical characterization of self-orthogonal linear CA, by considering the square matrix obtained by stacking the transition matrices of the CA and of its transpose, and determining when it is invertible. Interestingly, the stacked matrix turns out to have a circulant structure, for which there exists an extensive body of results to characterize its invertibility. Further, for the case of the binary alphabet we prove that irreducibility is a sufficient condition for self-orthogonality, and we derive a simpler characterization which boils down to computing the parity of the central coefficients of the local rule.
