Autocatalytic Closure for the Working Mathematician

In which emergence turns out to be a theorem.

1. You Already Know This

You have encountered this shape before. Consider three structures you command:

Generated groups. Given elements a, b in a group G, the subgroup ⟨a, b⟩ is the smallest subgroup of G containing a and b. You specify generators; the group structure determines what must exist.

Eigenspaces. Given a linear operator T on V, the eigenspace Eλ = ker(T − λI) is not chosen—it falls out of the operator. You specify the relation (T); the domain (Eλ) emerges.

Algebraic closure. The algebraic closure Q is not given in advance. You specify the property (every polynomial splits); the field that satisfies it is the smallest such extension. The domain is the answer to a closure condition.

In each case, you did not begin by specifying a domain. You specified a relation or property, and the domain fell out as a fixed point. This is routine mathematics—you do it without noticing the pattern.

Classical Mode Specify Domain → Discover Properties Domain D analyze Properties Given: V. Find: eigenvalues of T. Given: R. Determine: ideal structure. Given: X. Classify: π₁(X). Closure Mode Specify Property → Domain Emerges Relations / Rule close Domain D* ⟨a,b⟩ — generators → subgroup Eₗ — operator → eigenspace Q — "all polys split" → field

The left column is comfortable. The right column is exactly the same mathematics—closure under a specified operation—but you may not have noticed how often you already work in it.

2. A Puzzle

Here is a small reaction network. Molecules are circles; reactions are diamonds. Solid arrows connect reactants to reactions and reactions to products. Dashed red arrows show catalysis—which molecule enables which reaction. Molecules in blue are food: freely available from outside.

f₁ food f₂ food a b c d r₁ r₂ r₃ r₄ r₅ d is never produced Legend Food (free) Molecule Reaction Reactant/Product Catalysis

Six molecules: f1, f2 (food), a, b, c, d. Five reactions. One reaction (r5) needs molecule d as both reactant and catalyst—but nothing produces d. One pair (r1 and r2) forms a mutual dependency: r1 needs b as catalyst; r2 needs a as catalyst and produces b; but a is produced by r1.

The question: which reactions can sustain themselves from the food set alone?

3. The Mathematician's Instinct (And Why It Stalls Here)

Your training tells you: define the domain, then prove theorems about it. Fix a group G, then study its structure. Fix a topological space X, then classify its coverings. The domain comes first; results follow.

Here, this instinct stalls. What "exists" in the system depends on what can be produced. What can be produced depends on which reactions can fire. Which reactions can fire depends on what exists to catalyze them. You feel a circularity—and your training says: stop. Circularity means ill-definition.

But it isn't circular. It's a fixed-point equation on ℘(R).

You're looking for R′ ⊆ R such that R′ = Φ(R′), where Φ keeps exactly those reactions whose preconditions are met within the closure of R′ itself. The apparent circularity is the self-referential character of any fixed-point condition. You solve these routinely in other contexts—you just haven't seen this one before.

The domain is not given. The domain is the answer.

4. The Algorithm (Greatest Fixed Point)

Rather than searching bottom-up through 2|R| subsets, start from the top. Begin with all reactions; iteratively remove any reaction whose preconditions fail; repeat until stable. This is a greatest fixed-point computation: a monotone map on (℘(R), ⊇).

Step through it on the puzzle from Section 2.

Step 0
Click "Next Step" to begin. We start with R' = {r₁, r₂, r₃, r₄, r₅} and food F = {f₁, f₂}.
r₁ f₁ + f₂ → a cat: b r₂ a → b cat: a r₃ f₁ + a → c cat: b r₄ a + b → c cat: c r₅ d + f₂ → a cat: d Closure: {f₁, f₂} R′: {r₁, r₂, r₃, r₄, r₅}

Observe the structure. The operator Φ(R′) = {r ∈ R′ : ρ(r) ⊆ clR′(F) ∧ ∃x ∈ clR′(F), (x,r) ∈ C} is monotone on (℘(R), ⊇): if R′ ⊇ R″ then clR′(F) ⊇ clR″(F), so Φ(R′) ⊇ Φ(R″). Tarski guarantees a unique greatest fixed point. You just computed it.

The algorithm is O(|R|2 · |X|), not 2|R|. This is not an optimisation—it is a structural consequence of the monotonicity of Φ.

5. Now the Definition

Definition (RAF — Hordijk & Steel, 2004)

A Catalytic Reaction System (CRS) is a triple Q = (X, R, C) where X is a finite set of types, R is a set of reactions (each r with reactant set ρ(r) ⊆ X and product set π(r) ⊆ X), and C ⊆ X × R is the catalysis relation.

Given a food set F ⊆ X, a subset R′ ⊆ R is a Reflexively Autocatalytic, F-generated set (RAF) if:

  1. F-generated: ∀r ∈ R′ : ρ(r) ⊆ clR′(F)
  2. Reflexively autocatalytic: ∀r ∈ R′ : ∃x ∈ clR′(F) s.t. (x,r) ∈ C

where clR′(F) is the smallest superset of F closed under products of reactions in R′ whose reactants are available.

Read the definition again. No intrinsic property of the elements of X appears anywhere. The definition references only the maps ρ, π, and C—which are relational structure on a graph. Identity of a type x ∈ X is exhaustively determined by its relational position: what it catalyzes, what produces it, what consumes it. Nothing else.

This means the definition is substrate-independent. Any system that instantiates the same graph of reactant/product/catalysis relations is, formally, the same CRS. The labels change; the structure is identical.

Abstract Chemistry Ecology Knowledge Governance
x ∈ X Molecule type Species Concept / technique Institution / role
r ∈ R Chemical reaction Ecological interaction Derivation / synthesis Policy process
ρ(r) Reactants consumed Inputs required Prerequisite concepts Required inputs
π(r) Products generated Outputs produced New concepts derived Outputs / decisions
(x,r) ∈ C x catalyzes r x enables interaction x enables derivation x authorizes process
F Food molecules Abiotic resources Axioms / observations Constitutional basis
RAF Self-sustaining set Stable ecosystem Self-sustaining field Viable governance
Because the definition is purely relational, it applies to any substrate that instantiates these relations. The proof machinery—Tarski, monotonicity, lattice structure—transfers wholesale. No new theorems needed; only new labels.

6. Emergence, Proved

You do not need to accept "emergence" as a philosophical position. You can prove it as a theorem.

Definition (Irreducible RAF)

A RAF R′ is irreducible (an irrRAF) if it contains no proper non-empty sub-RAF. Equivalently: removing any single reaction from R′ destroys the RAF property.

Theorem (Emergent Property)

Let R′ be an irrRAF. The property "self-sustaining and autocatalytic" holds for R′ but fails for every proper subset R″ ⊊ R′.

Proof. By the definition of irreducibility, no proper non-empty subset of R′ is a RAF. The empty set is trivially not a RAF (it sustains nothing). Therefore self-sustenance is a property of R′ that fails for all R″ ⊊ R′. □

This is the precise algebraic definition of an emergent property: a property P of a set S that holds for S but fails for all proper subsets of S. No philosophy. No hand-waving. A property of the lattice.

The collection of all RAFs of a CRS forms a partially ordered set under ⊆, with the maxRAF at the top. The irrRAFs are the minimal non-empty elements. Each irrRAF is a node where emergence first appears—the smallest subsets at which a collective property crystallizes.

The Phase Transition

Kauffman (1986) conjectured, and Mossel & Steel (2005) proved for certain model classes, that in a random CRS with catalysis probability p, the probability that a RAF exists undergoes a sharp transition near p ∼ 1/n:

catalysis probability p P(RAF exists) 0 1 0 1 threshold p ~ 1/n no RAF no collective property RAF exists irreducible self-sustenance

Below the threshold: no subset of R is self-sustaining. Above: irreducible self-sustenance appears. The transition sharpens with n. This is the moment emergence crystallizes—not as metaphor, but as a phase transition in the lattice of subsets of R.

7. Coda

You didn't need to define the domain first. You needed to specify the relations—the catalysis graph, the reactant/product maps, the food set. The domain fell out as a fixed point of a monotone operator on a finite powerset lattice. Tarski did the rest.

This is how emergence works—not as philosophy, but as algebra. An irreducible RAF is a set whose defining property provably fails for every proper subset. That is a theorem, not an argument. And because the definition is purely relational—no intrinsic properties of X appear anywhere—the theorem holds for any substrate that instantiates the graph. Chemistry, ecology, knowledge, governance: same lattice, same fixed point, same emergence.

References. Hordijk, W. & Steel, M. (2004). Detecting autocatalytic, self-sustaining sets in chemical reaction systems. J. Theor. Biol. 227:451–461. Kauffman, S. A. (1986). Autocatalytic sets of proteins. J. Theor. Biol. 119:1–24. Mossel, E. & Steel, M. (2005). Random biochemical networks. J. Theor. Biol. 233:327–336.

Autocatalytic Closure for the Working Mathematician — B. Stephenson, 2026. No external dependencies.

← Web Annex