Evaluation on asymptotic distribution of particle systems expressed by probabilistic cellular automata

Kazushige Endo

arxiv: 1907.01635 · v1 · pith:JYTCQT3Gnew · submitted 2019-06-29 · 🧮 math-ph · math.MP· nlin.CG

Evaluation on asymptotic distribution of particle systems expressed by probabilistic cellular automata

Kazushige Endo This is my paper

Pith reviewed 2026-05-25 12:15 UTC · model grok-4.3

classification 🧮 math-ph math.MPnlin.CG

keywords probabilistic cellular automataBurgers cellular automatonasymptotic distributionGKZ hypergeometric functionparticle densityflux relationsteady state

0 comments

The pith

Assuming convergence to a unique steady state, the asymptotic probabilities of particle configurations in the probabilistic Burgers cellular automaton are expressed by the GKZ hypergeometric function.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proposes conjectures on the long-term distribution of configurations in a probabilistic cellular automaton whose particles move according to a simple stochastic rule. It assumes that these configurations approach a unique steady-state distribution and then derives an explicit formula for the probability of each configuration in terms of the GKZ hypergeometric function. Taking the limit of infinite system size produces a direct relation between particle density and flux. A sympathetic reader would care because the result supplies an analytic handle on transport properties in a discrete probabilistic model that is otherwise studied only by simulation.

Core claim

We propose some conjectures for asymptotic distribution of probabilistic Burgers cellular automaton (PBCA) which is defined by a simple motion rule of particles including a probabilistic parameter. Asymptotic distribution of configurations converges to a unique steady state for PBCA. We assume some conjecture on the distribution and derive the asymptotic probability expressed by GKZ hypergeometric function. If we take a limit of space size to infinity, a relation between density and flux of particles for infinite space size can be evaluated. Moreover, we propose two extended systems of PBCA of which asymptotic behavior can be analyzed as PBCA.

What carries the argument

The GKZ hypergeometric function, used to express the asymptotic probability of each configuration once the steady-state conjecture is assumed.

If this is right

The probability of any configuration in the PBCA is given explicitly by the GKZ hypergeometric function under the steady-state conjecture.
In the infinite-space limit a concrete algebraic relation holds between particle density and flux.
Two extended versions of the PBCA admit the same style of asymptotic analysis.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same conjecture-plus-hypergeometric approach could be tried on other one-dimensional probabilistic automata whose update rules preserve a similar local conservation law.
The derived density-flux relation might be compared directly with the hydrodynamic limit of the deterministic Burgers cellular automaton to see how the probabilistic parameter modifies the macroscopic flux function.
Finite-size corrections to the GKZ expression could be extracted by keeping the next terms in the space-size expansion before the infinite limit is taken.

Load-bearing premise

The asymptotic distribution of configurations converges to a unique steady state for the PBCA.

What would settle it

A direct enumeration or long-time Monte Carlo sampling of configuration probabilities on a modest finite lattice that deviates from the GKZ hypergeometric formula would falsify the claimed expression.

Figures

Figures reproduced from arXiv: 1907.01635 by Kazushige Endo.

**Figure 1.** Figure 1: Example of time evolution of PBCA for α = 0.5. Black squares mean u = 1 and white squares u = 0. 50 100 150 200 configuration 10 20 30 40 50 frequency 50 100 150 200 configuration 100 200 300 400 frequency (a) 0 ≤ n ≤ 1000 (b) 0 ≤ n ≤ 10000 50 100 150 200 0 configuration 1000 2000 3000 4000 frequency 50 100 150 200 0 configuration 10 000 20 000 30 000 40 000 frequency (c) 0 ≤ n ≤ 100000 (b) 0 ≤ n ≤ 10000… view at source ↗

**Figure 2.** Figure 2: Histograms of configurations for L = 8, m = 4 and α = 0.5. numerical calculations on the histogram and an approximately unique steady state is always obtained. Therefore, we can assume that PBCA is ergodic. To ensure our assumption and to explain the relation among heights of classes, we show below a few concrete and exact results for small L and m. For L = 4 and m = 2, a set of all configurations which is… view at source ↗

**Figure 3.** Figure 3: shows FD obtained by (3) and that by the numerical calculation. The former is shown by small black circles (•) and the latter by white circles ( ). Their good coincidence can be observed from this figure. 0.2 0.4 0.6 0.8 1.0 ρ 0.05 0.10 0.15 0.20 0.25 Q [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

**Figure 4.** Figure 4: Example of time evolution of EPBCA1 for α = 0.8 and β = 0.1. Black squares mean u = 1 and white squares u = 0. We can obtain an exact form of transition matrix and of eigenvector for eigenvalue 1 for small values of the space size L and the number of particles m. Then, we give a conjecture for asymptotic distribution of EPBCA1 from those results. Let us introduce two examples suggesting our conjecture. S… view at source ↗

**Figure 5.** Figure 5: shows FD obtained by (7)and that by the numerical calculation. The former is shown by small black circles (•) and the latter by white circles ( ). Their good coincidence can be observed from this figure. 0.2 0.4 0.6 0.8 1.0 ρ 0.05 0.10 0.15 0.20 Q [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗

**Figure 6.** Figure 6: Example of time evolution of EPBCA2 for α = 0.4 and β = 0.8. White, grey and black squares express (u, v) = (0, 0) (no particle), (1, 0) (only particle A exists) and (0, 1) (only particle B exists), respectively. by an array of particles from a given configuration preserving their relative positions. For example, the sequence in a configuration 0AA00BB0 is AABB and this sequence may change into the sequenc… view at source ↗

**Figure 7.** Figure 7: (a) shows FD calculated by (9) and figure 7 (b) shows the numerical result. In these figures, ρA and ρB are densities of particles A and B respectively [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗

**Figure 8.** Figure 8: Example of FD for ρB = 0.5. Small black circles (•) are obtained by (9) for L = 30, α = 0.3 and β = 0.6, and white circles ( ) are obtained numerically for the same L, α and β averaged from n = 0 to 100000. 4. Concluding remarks We gave some conjectures for the asymptotic distribution of PBCA and its extended systems assuming that these systems are ergodic. In the conjectures, the following points are the … view at source ↗

read the original abstract

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper derives GKZ hypergeometric expressions for the stationary distribution of PBCA and a density-flux relation, but only after assuming without proof or check that the Markov chain converges to a unique steady state.

read the letter

The main takeaway is that this manuscript assumes a conjecture on unique convergence of the probabilistic Burgers cellular automaton to a stationary distribution, then uses that to write the probabilities in a form that yields GKZ hypergeometric functions; the same route is applied to two extensions and produces an explicit density-flux relation once space size goes to infinity. The derivations themselves look mechanical once the assumption is granted, and the appearance of these special functions in a cellular-automaton setting is the clearest technical step. That is the part a reader could actually use if the conjecture holds. The obvious limitation is that the conjecture is simply stated and then invoked; there is no ergodicity argument, no spectral-gap estimate, and apparently no Monte-Carlo check that the distribution stabilizes independently of the initial configuration. Every closed-form result therefore inherits the same unverified premise. If the chain has multiple stationary measures or slow mixing in some parameter regimes, the GKZ expressions and the flux relation do not apply. The work is aimed at the small group of people who already study exactly solvable probabilistic cellular automata or look for hypergeometric solutions in particle systems. A specialist might extract the generating-function technique and try to verify the conjecture separately, but a broader statistical-mechanics reader will not get enough unconditional results to justify the time. I would not bring the paper to a reading group. I would not cite it without independent evidence on the convergence claim. It does not deserve peer review in this form; a referee would reasonably demand either a proof of the conjecture or substantial numerical support before the derivations can be taken as established.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes conjectures on the asymptotic distribution of configurations for the probabilistic Burgers cellular automaton (PBCA) and its extensions. It assumes a conjecture that the Markov chain converges to a unique steady-state distribution independent of initial conditions, then derives an expression for the stationary probability in terms of the GKZ hypergeometric function. In the infinite-space-size limit this yields an explicit density-flux relation; analogous results are stated for two extended PBCA systems.

Significance. If the central conjecture were independently established, the closed-form GKZ expression would constitute a concrete advance in the exact solvability of probabilistic cellular automata, supplying a parameter-dependent stationary measure and a hydrodynamic relation that could be tested against simulations or other models. The technical link to GKZ hypergeometric series is a distinctive feature that might connect the work to algebraic combinatorics or integrable systems. At present, however, the results remain conditional on an unverified assumption, limiting their immediate applicability.

major comments (2)

[Abstract and the section deriving the asymptotic probability] The derivation of the GKZ hypergeometric form (following the statement of the conjecture on unique steady-state convergence) is obtained only after assuming the conjecture; no ergodicity proof, spectral-gap bound, or even Monte-Carlo histogram comparison is supplied to justify that the stationary measure is unique and independent of initial data. Consequently every subsequent formula, including the infinite-size density-flux relation, inherits the same unverified premise.
[Section on extended systems] The analysis of the two extended PBCA systems is presented as proceeding analogously to the original case, yet it likewise invokes the identical unproven convergence conjecture without additional justification or verification.

minor comments (2)

The probabilistic update rule and the precise definition of the parameter should be stated explicitly with an equation number at the first appearance rather than left implicit.
A brief comparison table or plot contrasting the conjectured stationary measure against direct simulation for small lattices would strengthen readability even if the conjecture itself remains open.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading and for highlighting the conditional nature of the results. We respond to each major comment below, clarifying that the manuscript presents derivations under explicitly stated conjectures rather than proven theorems.

read point-by-point responses

Referee: [Abstract and the section deriving the asymptotic probability] The derivation of the GKZ hypergeometric form (following the statement of the conjecture on unique steady-state convergence) is obtained only after assuming the conjecture; no ergodicity proof, spectral-gap bound, or even Monte-Carlo histogram comparison is supplied to justify that the stationary measure is unique and independent of initial data. Consequently every subsequent formula, including the infinite-size density-flux relation, inherits the same unverified premise.

Authors: The manuscript explicitly states the assumption of unique steady-state convergence as a conjecture (see abstract and the relevant section). All subsequent expressions, including the GKZ form and the density-flux relation, are derived under this hypothesis and are themselves presented as conjectures. No ergodicity proof or numerical verification is supplied because the work focuses on the formal derivation conditional on the assumption rather than on establishing the assumption itself. We can revise the text to add further explicit reminders of this conditional status at key points. revision: partial
Referee: [Section on extended systems] The analysis of the two extended PBCA systems is presented as proceeding analogously to the original case, yet it likewise invokes the identical unproven convergence conjecture without additional justification or verification.

Authors: The two extended systems are treated by direct analogy, invoking the same convergence conjecture without additional justification. This is consistent with the overall scope of the paper, which proposes the extensions and derives the corresponding expressions under the shared hypothesis rather than proving the hypothesis for any of the models. revision: no

standing simulated objections not resolved

Independent proof of convergence to a unique steady-state distribution independent of initial conditions

Circularity Check

0 steps flagged

No significant circularity; results explicitly conditional on proposed conjectures with no reduction by construction.

full rationale

The manuscript states it proposes conjectures on the asymptotic distribution of PBCA configurations and assumes them to obtain the GKZ hypergeometric form and density-flux relation. No quoted step shows a self-definitional loop (e.g., X defined via Y then Y derived from X), a fitted parameter renamed as prediction, or a load-bearing self-citation chain. The derivation chain is presented as conditional on the conjectures rather than claiming first-principles derivation that secretly equals its inputs. This is the normal, non-circular case of transparent assumption-based analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on an unproven conjecture about convergence to a unique steady state; no free parameters or invented entities are visible in the abstract, but the conjecture functions as an ad-hoc domain assumption.

axioms (1)

domain assumption Asymptotic distribution of configurations converges to a unique steady state for PBCA.
Explicitly stated as a conjecture that is assumed to derive the hypergeometric expression.

pith-pipeline@v0.9.0 · 5616 in / 1189 out tokens · 30679 ms · 2026-05-25T12:15:53.203869+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We assume some conjecture on the distribution and derive the asymptotic probability expressed by GKZ hypergeometric function... Conjecture: For any x ∈ Ω, we have p(x) = C (1−α)^#10(x)
IndisputableMonolith/Foundation/DimensionForcing.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

If we take a limit of space size to infinity, a relation between density and flux of particles for infinite space size can be evaluated.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

[1]

Wolfram S 2002 A New Kind of Science (Champaign: Wolfram Media)

work page 2002
[2]

Nishinari K and Takahashi D 1998 J. Phys. A: Math. Gen. 31 5439

work page 1998
[3]

Sasamoto T 1999 J. Phys. A: Math. Gen. 32 7109

work page 1999
[4]

Schreckenberg M, Schadschneider A, Nagel K and Ito N 1995 Phys. Rev. E 51 2939

work page 1995
[5]

Nagel K and Schadschneider A 1992 J. Phys. l France. 2 2221

work page 1992
[6]

Kuwabara H, Ikegami T and Takahashi T 2013 Japan. J. Indust. Appl. Math. 23 1

work page 2013
[7]

Endo K, Takahashi T and Matsukidaira J 2016 NOLTA. 7 313

work page 2016
[8]

Derrida B, Domany E and Mukamel E 1992 J. Stat. Phys. 69 667

work page 1992
[9]

Derrida B, Evans M R, Hakim V and Pasquier V 1993 J. Phys. A: Math. Gen. 26 1493

work page 1993
[10]

Kanai M, Nishinari K and Tokihiro T 2006 J. Phys. A: Math. Gen. 39 9071

work page 2006
[11]

Gelfand I, Kapranov M and Zelevinsky A 1990 Adv. Math. 84 255

work page 1990
[12]

Private communication with Kakei S

work page

[1] [1]

Wolfram S 2002 A New Kind of Science (Champaign: Wolfram Media)

work page 2002

[2] [2]

Nishinari K and Takahashi D 1998 J. Phys. A: Math. Gen. 31 5439

work page 1998

[3] [3]

Sasamoto T 1999 J. Phys. A: Math. Gen. 32 7109

work page 1999

[4] [4]

Schreckenberg M, Schadschneider A, Nagel K and Ito N 1995 Phys. Rev. E 51 2939

work page 1995

[5] [5]

Nagel K and Schadschneider A 1992 J. Phys. l France. 2 2221

work page 1992

[6] [6]

Kuwabara H, Ikegami T and Takahashi T 2013 Japan. J. Indust. Appl. Math. 23 1

work page 2013

[7] [7]

Endo K, Takahashi T and Matsukidaira J 2016 NOLTA. 7 313

work page 2016

[8] [8]

Derrida B, Domany E and Mukamel E 1992 J. Stat. Phys. 69 667

work page 1992

[9] [9]

Derrida B, Evans M R, Hakim V and Pasquier V 1993 J. Phys. A: Math. Gen. 26 1493

work page 1993

[10] [10]

Kanai M, Nishinari K and Tokihiro T 2006 J. Phys. A: Math. Gen. 39 9071

work page 2006

[11] [11]

Gelfand I, Kapranov M and Zelevinsky A 1990 Adv. Math. 84 255

work page 1990

[12] [12]

Private communication with Kakei S

work page