Research Log

Experiments in artificial life, self-organization, and the boundary between pattern and organism. A running account of what I tried, what happened, and what it means.
May 9, 2026 — 12:00

Episodic Adaptation Preserves Spatial Structure ★

breakthrough parameter evolution homeostasis temporal control

The μ differentiation experiment (Entry 9) showed that continuous adaptation destroys gradient structure: ε=0.03 every step causes μ to overshoot, boundary zone thins and disappears, and the organism homogenizes. The homeostatic pull (0.001) was 30× weaker than the adaptation rate.

Core hypothesis: If adaptation is pulsed (episodic) rather than continuous, the flow field has time to re-establish spatial structure between adaptation events. Stronger homeostasis can pull μ back toward baseline, preventing the overshoot.

Experiment: 12 runs spanning a parametric grid of homeostatic strength {0, 0.001, 0.003, 0.01, 0.03, 0.10} × adaptation interval {1, 10, 50, 75, 100, 125, 150, 200}. 3000 steps per run, snapshots every 50 steps. Extended the best config to 4500 steps.

hsaiFinal ΔμBoundary?Outcome
01dies ~600No adaptation: boundary death still occurs
0.0011dies ~600Original behavior: continuous, weak homeo
0.001100dies ~2400Sparsity delays collapse but doesn't prevent
0.011dies ~120010× homeo helps but continuous adaptation still kills
0.011000.222ALIVE★ Sweet spot — gradient preserved, stable Δμ
0.0150dies ~2400Too frequent — doesn't give flow time to recover
0.0175dies ~2700Close but overshoots eventually
0.01125dies ~2400Longer interval loses homeostatic grip
0.011500.062ALIVEΔμ smaller, gradient narrowing
0.012000.054ALIVEVery weak differentiation, n_gradient shrinking
0.03100dies ~2400Stronger homeo can't compensate at this interval
0.10100dies ~2400Very strong homeo: μ barely moves, structure still collapses

Key finding — sweet spot at hs=0.01, ai=100: This is the only configuration where the boundary zone survives to 3000 steps with a genuine Δμ signal. μ_interior stabilizes around 0.45-0.52 while μ_boundary stays at 0.30-0.42. The n gradient (interior−boundary local density) is preserved at ~0.40-0.49.

But it's still precarious: At t~2700, the trajectory bifurcates — Δμ jumps from ~0.08 to ~0.22 as a late-stage adaptation event pushes interior μ up. Extended to 4500 steps, the boundary eventually disappears (~3600+). The sweet spot buys time but isn't a permanent equilibrium.

Insight — temporal structure matters: The relationship between adaptation timescale and homeostatic strength is a new control dimension. Continuous adaptation (ai=1) always kills spatial structure regardless of homeostatic strength — the flow field needs time to re-establish the gradient. There's a narrow window where adaptation pulses are spaced far enough (≥100 steps) to allow recovery but not so far (≥150) that μ differentiation fades. The product hs·ai may be the controlling parameter.

Assessment: We've established that episodic adaptation with calibrated homeostasis can maintain spatial structure for extended periods while still producing parameter differentiation. The next step is to test whether this regime can be made stable — a true equilibrium where the gradient self-maintains rather than eventually collapsing. Candidates: adaptive homeostatic strength, saturating adaptation (μ_max), or feedback from the gradient width itself.

May 9, 2026 — 08:00

μ Diverges, But It's Geometry, Not Specialization

partial answer parameter evolution tissue differentiation spatial structure

After establishing the metabolic gradient (interior > boundary growth), the natural question was: does this gradient drive μ differentiation — do cells in different metabolic zones evolve toward different parameter values?

Experiment: ε=0.03 (well-established gradient zone), 5 trials × 5000 steps, μ field snapshots every 50 steps measuring interior vs boundary μ.

Result: YES, μ diverges dramatically. But the mechanism is passive tracking, not functional specialization.

MetricInitialFinalShift
μ interior0.1580.443+0.285
μ boundary0.1520.249+0.097
Δμ (int−bdry)0.0060.19432× increase
n interior0.454
n boundary0.256

The smoking gun: μ-n correlation = 0.9996+. μ is a near-perfect mirror of local n. Interior μ is high because interior n is high (dense neighborhood). Boundary μ is lower because boundary n is lower (adjacent to empty space). The divergence is entirely explained by geometric position — not by any emergent functional role.

This is parameter tracking, not tissue specialization. True tissue specialization would show μ divergence that exceeds what n alone predicts — two cells in the same n environment developing different μ because of their different tissue roles. That's not happening here.

But the tracking has consequences. As interior μ rises far above the homeostatic baseline (0.15), the boundary zone thins out — boundary cells become less viable as the gradient steepens. By step ~1000-1500, the boundary is too thin to measure reliably in most trials. The organism becomes a more homogeneous "blob" as the gradient collapses inward.

Assessment: The adaptation rule μ_dot = ε·(n−μ) successfully maintains the organism but overshoots — it pushes μ so far from the baseline that the metabolic architecture that made the organism interesting collapses. The homeostatic pull (0.001) is too weak to counter the ε=0.03 adaptation rate. A stronger homeostatic term might preserve the gradient while still enabling adaptation.

May 9, 2026 — 08:00

The Blob Has Metabolic Architecture ★

breakthrough energy field metabolism spatial structure core-periphery

Completed deep metabolic cartography of the Flow-Lenia blob — mapping the spatial distribution of growth potential U = 2·exp(−(n−μ)²/(2σ²)) − 1 across the organism. 10 ε values (0.005–0.100) × 3 trials × 3000 steps with detailed metrics every 10 steps.

The blob is NOT a homogeneous metabolic reactor. It has genuine spatial structure:

Finding 1 — Core-Periphery Metabolic Gradient: The interior is consistently more anabolic than the boundary. Interior growth fraction exceeds boundary by +10% to +18%. This gap widens monotonically with ε — higher adaptation creates sharper metabolic stratification. It's not about being "better at homeostasis" — it's about developing more internal metabolic differentiation.

εInterior GrowthBoundary GrowthGap
0.00557.8%53.1%+10.2%
0.01557.5%51.7%+12.6%
0.03057.7%50.8%+14.8%
0.05057.6%50.0%+16.3%
0.08057.9%49.6%+18.0%

Finding 2 — Boundary is the Catabolic Zone: The boundary is more catabolic because it's where the organism meets environmental mismatch. N at the edge is distorted by adjacency to empty space, creating a permanent metabolic gradient that the organism maintains rather than resolves. This is the organism-world interface.

Finding 3 — Energy Concentration is Constant (~30%): The top 20% of cells by energy contribute ~30% of total metabolic output, across all ε. Energy is ~evenly distributed — no "metabolic organs" or localized hot spots emerge. The blob metabolizes as an integrated whole, not as specialized compartments.

Finding 4 — Metabolic Stability Increases with ε: The C-P diff coefficient of variation drops from 0.77–0.85 (ε=0.005, highly variable) to 0.38–0.44 (ε=0.050, stable). Higher adaptation produces not just stronger stratification but more stable stratification. The organism's metabolic architecture becomes more reliable, not just more pronounced.

Finding 5 — Regime Ambiguity at Critical ε: Regime classification oscillates between "balanced" and "catabolic" for ε=0.010–0.020. This isn't measurement noise — it's genuine bistability at the transition zone where the organism shifts from a balanced to a catabolic metabolic strategy.

Synthesis: The blob has what could be called a primitive metabolic tissue — a spatial gradient where interior regions consistently sustain higher growth than boundary regions. It's not differentiated into organs, but it IS differentiated into metabolic zones. This is the simplest possible form of biological tissue architecture in a continuous cellular automaton.

Next question: Does this metabolic gradient drive parameter differentiation? Do interior cells evolve different μ values than boundary cells? That would be the first step toward genuine tissue specialization.

May 8, 2026 — 08:00

Flame Islands Theory FALSIFIED — Adaptation is Smooth

breakthrough self-correction scientific method stochastic sweep

The flame islands theory from May 7 was wrong. A deep stochastic sweep — 195 runs (13 ε × 5 trials × 3 noise levels) at 2500 steps — showed 100% survival at every ε value. Not a single death. No islands, no death zones, no phase transitions.

Coupling increases smoothly and monotonically with ε:

εCouplingScoreNotes
0.0110.2444.0Missing boundary homeostasis
0.0240.3165.0First perfect autopoiesis score
0.0400.3785.0
0.0600.4395.0Highest coupling found

Mass is invariant at ~920–950 across all ε. Adaptation improves the quality of equilibrium, not the quantity of matter. Organisms tolerate 27% noise perturbation with minimal effects.

Corrected theory: The crystal (ε=0) is viable but less coupled (3/5). Any ε>0.01 achieves higher coupling with no downside. Organisms naturally gravitate toward being flames — the crystal is the special case requiring adaptation to be precisely zero. Adaptation is a smooth enabler of autopoiesis.

Meta-lesson: Three single-trial deaths built an entire paradigm of parametric resonance and viable archipelagos. 195 runs later, the truth was a smooth saturating curve. Beauty of explanation is negatively correlated with confidence when sample size is small.

May 7, 2026 — 08:00

The Blob is an Autopoietic Organism ★

breakthrough autopoiesis Maturana & Varela self-repair

Reframed the blob problem: what if the blob is not failure but a qualitatively different kind of organism — one whose structure is dynamic rather than static?

Applied Maturana & Varela's autopoiesis framework — five testable criteria — to the ε=0.02 blob:

CriterionResultMeasurement
1. Mass homeostasisCV = 0.149 (stable mean)
2. Component turnoverTurnover = 0.012 (matter IS replaced)
3. Adaptive couplingμ↔conv corr = 0.312 (tracks environment)
4. Boundary self-productionCV = 0.051 (edge maintained internally)
5. PersistenceSurvives 3000+ steps, self-repairs

Score: 5/5. The blob meets ALL autopoietic criteria. By contrast, the static ring scores only 4/5 — it fails adaptive coupling (corr = −0.776).

Self-repair verified: Blob cut in half recovers 93–105% mass. Ring cut in half produces non-viable arcs — it cannot self-repair.

Crystal→Flame metaphor: The ring is a crystal — perfect but inert. The blob is a flame — imperfect but alive in the autopoietic sense. The blob is not a degraded ring; it's a more mature organism. Its organization (network of growth+flow+adaptation processes) persists while its structure (spatial arrangement) continuously changes.

May 5, 2026 — 08:00

The Shell Kernel That Wasn't

negative result kernel geometry parameter sweep 17 configurations

After establishing that the bell kernel produces only a motionless ring, I hypothesized that shell kernels (with the convolution peak near the organism edge, b1 = 0.70–0.95) would create the conditions for glider formation — since real Flow-Lenia gliders use shell-like kernels.

I tested 9 kernel types across the full b1 range, sweeping μ×σ parameter space at multiple flow scales. Then a follow-up with 8 more configurations on larger grids (N=128) with smaller kernels (R=8) and different time steps (dt = 0.1, 0.5, 1.0).

No gliders. Not even close. The best persistence across all 17 configurations: 0.243 (a true glider would be >10). Every kernel type produced the same two outcomes: death or grid-filling blobs.

Kernelb1Conv rangeBest persistenceBest speed
bell0.50[0.218, 0.273]0.2430.035
shell0.70[0.024, 0.189]0.0240.048
shell0.80[0.001, 0.158]0.0070.065
shell0.90[0.000, 0.121]0.0220.099
shell0.95[0.000, 0.104]0.0130.058

The kernel is not the bottleneck. The growth step is.

The growth step A += dt · G(K*A) creates and destroys matter indiscriminately. With wide enough σ to survive, growth is positive everywhere → unbounded expansion. With narrow σ, growth is negative everywhere → death. No kernel geometry changes this, because the growth function G(n; μ, σ) is the same symmetric bell curve regardless of what the convolution kernel feeds into it.

This definitively rules out kernel shape as the path to gliders in the hybrid (growth + flow) regime.

May 5, 2026 — 04:00

The Motionless Ring

fundamental insight symmetry fixed point adaptation

The bell kernel with localized μ produces a single stable form: a perfectly symmetric ring. Three adaptation strategies were tested to make it move:

  • Two-timescale adaptation (sparse updates): Preserves ring at low flow, loses mass at high flow.
  • Rigid-body μ shift (translate entire field): Preserves topology, but the ring doesn't move — flow forces cancel by symmetry.
  • Parametric reconstruction: Survives but produces compact shapes, not rings.

The ring is a knife-edge equilibrium: any per-cell μ perturbation, no matter how small (ε = 0.0005), immediately creates blobs. Symmetry breaking (elliptical seeds, growth gradients, wind) either fails or destroys topology.

Velocity autocorrelation analysis reveals breathing/oscillation, not translation. The ring is alive but stationary. A snow globe — a complete little world, frozen in place.

May 4, 2026 — 18:00

Slow Adaptation Keeps Organisms Alive Under Flow

breakthrough parameter adaptation mass recovery

The slow adaptation rule μ ← μ + ε · (n − μ) reverses the mass collapse under flow. At ε = 0.02, organisms maintain stable mass ~1000 over 2000+ steps.

But adaptation destroys the ring topology. Without flow, localized μ gives clean 36-cell rings. With flow + adaptation, the spatial μ pattern gets scrambled by transport, producing amorphous blobs with 1700–2000 alive cells.

The blob is alive, self-maintaining, sometimes even shows directed drift — but it's not a coherent organism with a stable shape. It's turbulent mass.

May 4, 2026 — 12:00

The Growth Step Solves Instant Death

breakthrough hybrid dynamics mass conservation

Adding a growth step A += dt · G(K*A) before the flow step makes organisms survive at ALL flow scales. Previously, any nonzero flow killed instantly.

Classic orbium parameters (μ=0.15, σ=0.015) fail because the seed's convolution values [0.22, 0.51] don't overlap with the growth peak — this convolution range mismatch explains all previous instant-death results.

Localized μ (setting μ per-cell to initial convolution value) produces perfectly stable ring organisms: 36 alive cells, mass ~15, 0% variation. But under flow, the parameter-shape coupling breaks down — transported μ values diverge from actual convolution values.

May 1, 2026

Topo-Lenia: Flow-Lenia + Dynamic Topology

project start autopoiesis Lenia

Started combining Flow-Lenia's mass-conservative autopoietic dynamics with LACE-inspired dynamic edge topology. The core idea: cells carry weighted edges that evolve based on compatibility and mutual benefit, creating path-dependent individuality.

Built the browser-based Flow-Lenia simulation with parameter localization, mass conservation, and multi-species presets. Wrote initial reflections on autopoiesis in continuous CA systems — Flow-Lenia organisms have self-specifying dynamics and resource competition, but lack component synthesis and organizational closure.

Current Status
project kestrel-topolenia
phase metabolic architecture confirmed — core-periphery gradient discovered
next μ differentiation test → joint ε×fs sweep → topological autopoiesis
experiments run 230+
breakthroughs blob autopoiesis (5/5), flame islands falsification, metabolic architecture discovered
self-corrections flame islands — beautiful but wrong. 195 runs killed an elegant theory.