Exact Diffusion Design for Maximally Collective Stable Turing Patterns in Binary-Complex Mass-Action Networks
Read the paper (PDF, 19 pages) · Source and verification (ZIP) · Mathematical supplement (PDF, 19 pages) · Data and certificate tables · Complete reproducibility package · Zenodo record
Plain-English explanation · Read on this page · Collaboration
Abstract
Multicomponent Turing theories often seek a lower-dimensional unstable subsystem responsible for patterning. For every n ≥ 4, we construct an n-species binary-complex classical mass-action topology for which, throughout every positive-equilibrium realization, every principal subsystem below order n−1 is Hurwitz while one order-(n−1) subsystem is unstable. We prove a general principal-minor diffusion-ray theorem and obtain an exact network law: for any homogeneously stable realization, a stationary crossing occurs along a positive diagonal diffusion ray exactly when d_Z > 8h_Z ∑_{j=2}^{n−2} d_j/h_j. This gives the nonattained fixed-equilibrium contrast infimum, the unit-equilibrium value 8(n−3), and the topology-wide stationary trade-off χ_Dχ_H > 8(n−3). At unit equilibrium, an exact rational design undergoes a primary supercritical bifurcation with locally exponentially stable positive patterns in the fixed-integrated-mass class. A second exact family redistributes heterogeneity between diffusion and equilibrium concentrations so that both contrasts are Θ(√n), matching the topology-specific stationary minimax lower bound in exponent. The proof uses explicit flux-cone design, exhaustive strongly connected block classification, conservation-compatible center-manifold reduction, and inspectable symbolic certificates. The networks are synthetic; the bifurcation, stability, and robustness conclusions are local for each fixed dimension, and the stationary diffusion law does not classify arbitrary wave instability.
The idea in plain English
Diffusion usually smooths out differences, but in a reacting system it can destabilize a uniform state and create a spatial pattern. Many explanations trace that instability to a small group of interacting components. This paper constructs systems in which the smallest unstable group contains all but one of the species.
For every number of species n ≥ 4, the construction uses ordinary mass-action reactions with at most two molecules on either side. Its collective character holds at every positive equilibrium realization of that fixed reaction topology. The networks are mathematical designs, not identified biochemical systems.
The paper derives an exact threshold for stationary instability in terms of diffusion rates and equilibrium concentrations. It then separately proves that particular designs produce small, positive, locally stable patterns near their first stationary bifurcation.
A second design shares the required variation between diffusion rates and concentrations. Both contrasts grow like the square root of the number of species, reaching the best possible growth exponent for stationary instability within this topology. The optimal constant and the best possible stable-pattern threshold remain undetermined.
The precise result
For every n ≥ 4, an explicit n-species binary-complex classical mass-action topology has every principal subsystem of order below n−1 Hurwitz at every positive-equilibrium realization, while one subsystem of order n−1 has a positive real eigenvalue.
Write h_i = 1/x_i* for inverse equilibrium concentrations and D = diag(d_i) for strictly positive diffusion coefficients. At a homogeneously stable realization, a stationary crossing along the diffusion ray occurs exactly when d_Z > 8h_Z ∑_{j=2}^{n−2} d_j/h_j.
The fixed-equilibrium stationary diffusion-contrast infimum is 8h_Z ∑_{j=2}^{n−2} 1/h_j and is not attained. At unit equilibrium it equals 8(n−3). Across this topology, diffusion contrast χ_D and equilibrium-concentration contrast χ_H satisfy χ_Dχ_H > 8(n−3) for stationary instability.
An exact rational unit-equilibrium design on a one-dimensional Neumann interval undergoes a primary supercritical stationary bifurcation. For each fixed n, sufficiently small positive distance past onset gives positive patterns that are locally exponentially stable in H¹ within the fixed-integrated-mass class defined by the semipositive conservation law.
A second exact family adjusts equilibrium concentrations and physical diffusion together. Both contrasts are Θ(√n), matching the topology-specific stationary minimax lower bound in exponent. This does not establish the optimal constant or the complete trade-off frontier.
Why the result holds
An explicit positive flux cone describes every equilibrium realization of the chosen topology. The resulting Jacobian factorization lets the proof classify every possible strongly connected principal block and show that all blocks smaller than n−1 are stable.
The signed principal minors make the determinant along a positive diagonal diffusion ray especially simple. After removing the conservation factor, every higher coefficient is positive, leaving a single first coefficient to decide the unique positive stationary crossing.
For the selected rational design, separate half-plane certificates control every noncritical Fourier mode. A conservation-compatible center-manifold calculation then establishes a positive crossing derivative and negative cubic coefficient, yielding the locally stable supercritical branch.
The second family changes concentrations, diffusion and the conservation gauge consistently. Exact spectral and cubic-sign certificates persist over a specified parameter interval, allowing the two contrasts to be balanced at square-root growth.
Scope and limitations
The topology is synthetic and the reaction law is classical mass action. The result does not classify arbitrary networks, establish biological realizability or provide a universal diffusion threshold.
The exact ray law concerns stationary zero-eigenvalue crossings. It does not rule out every oscillatory or wave instability, and it does not alone prove the existence or stability of nonlinear patterns.
Homogeneous stability is imposed on the stoichiometric subspace because the full Jacobian has a conservation zero. The conserved quantity has zero weight on X₁; fixed integrated mass is not total mass of all species or a global concentration bound.
Stable-pattern and robustness conclusions are local for each fixed dimension and fixed conserved integral. They do not give global attraction, a dimension-uniform basin, or stability for arbitrary perturbations of the conserved mass.
The interval has Neumann boundary conditions and discrete spatial modes. The construction tunes a diffusion multiplier to onset and does not claim intrinsic finite wavelength selection as the domain grows.
The stationary contrast infimum is strict and unattained. The threshold for a locally stable nonlinear branch can be higher, and its exact infimum is not determined. The square-root construction is exponent-optimal only within this topology and stationary setting.
Subsystem theories, topology-based diffusion constraints, monomial equilibrium parametrizations and local bifurcation theory have established predecessors. Exact symbolic checks supplement the proof; they are not human peer review or a formal proof-assistant certificate.
Status and AI use
AI systems assisted with exploratory theorem discovery, symbolic derivation, code generation, adversarial testing, proof organization and prose editing. The author takes responsibility for the content. Exact computations and internal AI checks supplement the mathematical arguments and do not constitute conventional human peer review.
Zenodo version 1.0.11 was published on 12 September 2026. The main paper and mathematical supplement retain August 2026 on their title pages. The companion software has its own separate Zenodo DOI.
Read the paper
Open PDF in a full window · Download the original from Zenodo
Supporting materials
Collaborate on this result
I welcome substantive collaboration on independent verification of the all-dimension spectral and bifurcation arguments, comparison with earlier reaction-network results, and preparation for journal review. Expertise in reaction–diffusion equations, chemical reaction network theory or exact symbolic certification would be especially useful.
Authorship and responsibilities would be discussed in relation to the substantive contributions made.
Cite this paper
Alec Kriebel. (2026). Exact Diffusion Design for Maximally Collective Stable Turing Patterns in Binary-Complex Mass-Action Networks (Version 1.0.11; manuscript dated August 2026) [Preprint]. Zenodo. https://doi.org/10.5281/zenodo.22729355
Background and prior work
- Satnoianu, Menzinger and Maini, Turing instabilities in general systems: General multicomponent theory already allows an unstable principal subsystem as large as n−1. The present construction realizes that endpoint throughout a constrained mass-action topology.
- Anma, Sakamoto and Yoneda, Unstable Subsystems Cause Turing Instability: Established subsystem criteria distinguish stationary and wave mechanisms. The present exact diffusion-ray law addresses the stationary case.
- Diego, Marcon, Müller and Sharpe, Key Features of Turing Systems Are Determined Purely by Network Topology: Earlier graph methods relate reaction-network topology to diffusion constraints and pattern properties. Topology-based Turing design itself is not new here.
- Haas and Goldstein, Turing’s Diffusive Threshold in Random Reaction–Diffusion Systems: Studies diffusion thresholds and collective behavior in random systems, providing context distinct from an exact family of mass-action realizations.
- Conradi, Mincheva and Uecker, Conditions for Spatial Instabilities and Pattern Formation from Monomial Steady State Parameterizations: Develops related steady-state and mass-action methods for spatial instability, with a nine-species application and numerical nonlinear analysis.
- Crandall and Rabinowitz, Bifurcation from Simple Eigenvalues: The local branch construction relies on classical simple-eigenvalue bifurcation hypotheses, together with conservation-compatible parabolic stability theory.