Queued Topics for Later Exploration
Findings from daily research that lead down a different focus track. Saved here for later exploration.
From Session 1
- Holland's Echo model — A classic SFI complex adaptive system model. How does it handle (or fail to handle) multi-scale composition?
- Langton's edge of chaos (Lambda parameter) — Does the edge of chaos shift when you allow multi-scale composition? Is the edge of chaos a network restructuring event?
- Stigmergy — Indirect coordination through environmental modification. ANT-compatible (environment as actor). Mechanism for cross-scale interaction.
- Deleuze & Guattari's rhizome — Latour references this. No center, no hierarchy. How does this differ from a scale-free network?
- Blaise Agüera y Arcas — Emergence in neural systems, social aggregation in computational systems.
- Renormalization group (Wilson) — Formal method for relating descriptions at different scales in physics. Could it be adapted for ALife?
- von Neumann's universal constructor — Original self-replication model. The constructor builds itself, which is a strange loop.
- Kauffman's NK model and fitness landscapes — How do fitness landscapes change when actors are defined relationally?
- Capra & Luisi, The Systems View of Life — Systems thinking, autopoiesis, origins of life.
From Session 2
- Downward causation and computational irreducibility — If high-level patterns causally influence low-level components, does that make the system more or less irreducible? Can we quantify this?
- Multi-scale autopoiesis — Systems producing systems at different scales. Is this the mechanism for complexification? Test via simulation.
- Tangled hierarchy formalization — How to represent a tangled hierarchy computationally? Not a tree, not a graph, but a level-crossing feedback structure.
- Gödel's incompleteness and ALife — Hofstadter connects Gödel to strange loops. Does formal undecidability have implications for what ALife simulations can produce?
- Luhmann's social autopoiesis — Niklas Luhmann applied autopoiesis to social systems. Connection to ANT's social networks.
- Memes and evolutionary stigmergy (Blackmore/Dawkins) — Memes as stigmergic traces that propagate, mutate, and evolve. How does this differ from static stigmergic traces? Can stigmergic traces in a simulation evolve? Connection to quasi-objects (traces that transform through circulation). Memes as a bridge between stigmergy and Darwinian replicators.
From Session 8
- Environmental physics coupling (the Mahadevan mechanism) — DONE (Session 9). Researched and specified as the concrete negative-feedback mechanism for the trace→actor crossing (H7). New concept file
environmental-physics-coupling.md. sim07 design sketched (transport field + M_c threshold). NEXT PRIORITY: implement sim07 Part-by-Part per its DESIGN.md. - The 20-year stigmergic-construction modeling lineage — DONE (Session 9). Documented that Deneubourg (1977) → Bonabeau (1997) → Ladley & Bullock (2004) all share sim06's limitation (material doesn't influence movement). (2026-07-27: the lineage documentation stands as literature, but the conclusion drawn here — "sim06's null result is a known field-wide gap" — does not. sim06's null had a separate local cause, a detector that could not fire.) Reference added to references.md.
- sim07: implement the transport field + M_c phase transition — DONE (Session 10). Implemented sim07 per DESIGN.md. NULL result: no phase transition in M_c — scalar structure-sourced transport fragments rather than consolidates (stability 0.876→0.739, pillars 57→128 as M_c drops); crossing never fires; self-repair tracks the deposit rule not T (circularity safeguard fails). H7 refined ×3: the crossing needs DIRECTED transport and/or an external multi-rate driver, not just a structure-sourced scalar. sim07.py, README.md, visualize.html, results.json all written. NEXT PRIORITY: sim08 (external oscillation).
- sim08: external oscillation as the energy source for transport — TOP PRIORITY for next nightly sessions. sim07's null showed structure-sourced scalar transport has the wrong sign for consolidation (it disperses the cue that recruits deposits). The Mahadevan mechanism's energy comes from OUTSIDE the structure (diurnal temperature oscillation), and the flow is DIRECTED (along channels), not an isotropic scalar. sim08 should add an external oscillation the structure can rectify into directed flow, and model the structure's shape as a channel (not just its mass). Test whether the crossing fires only when the external driver is present — making the multi-rate environment (H4) the energy source for the trace→actor crossing. This is the concrete test of H4 ↔ H7 coupling.
- Morphospace validation — sim07's predicted consolidated morphology (few large vented pillars) should be compared to the Mahadevan morphospace (Ocko/Heyde/Mahadevan 2019). A match = cross-validation; a mismatch = the lumped model is insufficient. Could compare simulated vs. real mound shapes quantitatively.
- Assembly theory connection (still queued from Session 6) — Mathis et al. 2024 / Cronin-Walker assembly index as a metric for ALife organization complexity. Could the M_c threshold be characterized by an assembly-index jump?
From Session 3
- Heylighen's varieties of stigmergy — Full taxonomy (quantitative/qualitative, sematectonic/marker-based, transient/persistent, broadcast/narrowcast). How do these map to computational stigmergic mechanisms? Which varieties are most relevant for ALife?
- Ecosystem engineering vs. niche construction — The distinction between Jones et al.'s ecosystem engineering and Odling-Smee's niche construction. NCT emphasizes evolutionary feedback; EE emphasizes ecological impact. Which is more relevant for multi-scale ALife?
- Extended evolutionary synthesis debate — The controversy over whether niche construction requires new evolutionary theory or is accommodated by standard theory. Parallel question for ALife: does stigmergy require new simulation paradigms or is it already present in any dynamic-environment simulation?
- Chemical Organization Theory — Heylighen mentions Dittrich & Fenizio's framework for agentless stigmergic coordination in chemical reaction networks. Could this formalism be adapted for ALife composition?
- Braitenberg vehicles and stigmergic cognition — Heylighen's analysis of Braitenberg vehicles as stigmergic systems. Connection between stigmergy and embodied cognition.
- Trace decay rate optimization — The transient/persistent trace trade-off. Is there an optimal decay rate for the trace→actor crossing? How does this relate to Wolfram's computational irreducibility?
From Session 4
- Chemical Organization Theory (Dittrich & Fenizio) — Still queued. Agentless stigmergic coordination in chemical reaction networks. Could provide formalism for ALife composition. Next session priority.
- Multi-scale NK model — Define NK landscapes at each scale with cross-scale interactions reshaping landscapes. How do dynamic epistatic networks behave? Does multi-scale landscape structure produce open-ended evolution?
- Gavrilets' holey landscapes — High-fitness genotype networks as alternative to rugged landscape view. How do holey landscapes behave with niche construction? Does the dynamic landscape view change the holey/rugged distinction?
- Holland's "Signals and Boundaries" (2012) — His last monograph. Co-evolution of signals and semi-permeable boundaries. Connection to our stigmergy + autopoiesis synthesis.
- Implementing Smith & Bedau's 8th CAS property — They proposed it in 1997 but never implemented it. Our sim02 shows naive stigmergy doesn't do it. What mechanism would? The autopoietic crossing (H7) is the candidate. Design sim03/sim04 to test.
- Trace competition — Multiple trace types that interact/compete. Sim02 used a single trace field. Multiple trace types might prevent monoculture convergence and enable multi-scale structure.
- Kaznatcheev's hard/soft landscape distinction — Which NK parameters produce open-ended dynamics? Sweep K and N to find the boundary. Connection to edge of chaos (Langton).
- Ecosystem engineering vs. niche construction (still queued) — Jones et al. vs. Odling-Smee. Which is more relevant for ALife?
- Heylighen's varieties of stigmergy (still queued) — Computational mapping of stigmergy taxonomy.
From Session 5
- Fontana & Buss's AlChemy (lambda calculus chemistry) — DONE (Session 6). Implemented sim05. Stable species sets emerge; L2 coexistence 2/6 (corrected 2026-07-27 — the originally reported 0/6 was a measurement artifact). Unbounded space still looks insufficient (H10), but the evidence is much weaker than 0/6 implied.
- Per-compartment catalysis — Our sim04 shared catalysis rules across all compartments. Vasas et al. generate catalysis independently per compartment. Does independent catalysis produce more between-compartment diversity? Test in sim05/sim06.
- P_catalyze tuning for distinct cores — Our sim04 used P=0.005, suspected too high (one large core). (2026-07-27: that suspicion came from a run that was not reproducible — catalysis was derived from Python's randomized
hash(). sim04 is now deterministic and gives 3 cores in both conditions; re-derive the diagnosis from the current results before sweeping.) Vasas used P''=0.0025 and still had difficulty finding distinct cores. Sweep P to find the regime where distinct cores form. - Expanding the adjacent possible — Kauffman's concept. Each novel core extends the "shadow" of possible reactions. Can we measure the adjacent possible in our simulations? Does the evolving network explore more of it than the fixed network?
- Holland's tagged urn model implementation — Holland proposed it but never tested it. Could implement as sim06: urns with semi-permeable boundaries containing tags, with GA-evolved classifiers. Test whether nested boundaries emerge. PRIORITY: this could provide the composition mechanism that sim05 showed is missing.
- From "one bit" to open-ended — The core limitation from Vasas et al. How to move beyond 1-bit heritable information? Template replication (RNA world) is the biological answer. What is the ALife answer? Multiple interacting cores? Compositional inheritance? Tag-based heredity?
- Multiple attractors ≠ evolvability — Vasas found networks with inhibition had multiple attractors but they were NOT selectable (periodic/chaotic transitions overrode selection). Explore this: what makes attractors selectable vs. not? Stability, heritability, differential fitness.
From Session 6
- Explicit composition mechanisms for L2 — (Premise corrected 2026-07-27: sim05 does NOT show L2 failing to emerge. Corrected, coexistence occurs in 2/6 pairs, stable across survival thresholds 0.45–0.70. This item was flagged TOP PRIORITY on the strength of 0/6; it is still interesting but no longer urgent, and the more informative question is now what distinguishes the 2 pairs that coexisted from the 4 that didn't.) What mechanisms would produce L2 reliably? Candidates: (a) stigmergic bridges between organizations, (b) autopoietic boundaries that protect during interaction, (c) explicit selection for composability.
- Measuring the adjacent possible in AlChemy — Sim05 showed each L1 run explores 112–162 species (corrected 2026-07-27; the earlier 246–930 was inflated ~3–6× by a non-alpha-invariant species equality). Can we measure how much of the "adjacent possible" (Kauffman) is explored? Does the rate of novel species discovery follow a power law? Does it slow down (converging) or stay constant (exploring)?
- Mutual destruction as creative process — (Premise retracted 2026-07-27. Only 1 of 6 pairs now ends in mutual destruction, and its final population is the SMALLEST at 10 species, not the largest — the reverse of the original observation. The "89-90 unique vs 6-23" comparison came from inflated, non-alpha-invariant species counts and cannot be reproduced. Re-derive before pursuing.) The underlying question — whether cross-organization interaction generates novelty that could be harvested without destroying the parents — is still open. Can we harness this novelty without destroying the parents? Autopoietic boundaries might protect parents while allowing cross-organization interaction.
- Assembly theory connection — Mathis et al. 2024 reference Cronin/Walker's assembly theory. Assembly theory quantifies selection by molecular complexity. Could assembly index be a metric for ALife organization complexity? Connection to our multi-scale composition metric needs.
- Krzyszewski & Mikolov (2022) — self-reproducing metabolisms as recursive algorithms — Referenced in Mathis et al. 2024. Emergence of self-reproducing metabolisms in artificial chemistry. Could connect to our autopoiesis + stigmergy synthesis.
From Vance (2026-07-22) — The termite mound principle
- Heterogeneous environment with substrate state transitions — Vance's key insight: unbounded space isn't sufficient because it's homogeneous. The termite mound works because inert mud becomes a dynamic actor (affects temperature, chemistry) once it crosses an organizational threshold (mass). Sim06 should model a HETEROGENEOUS environment where substrates have state transitions (inert → active) triggered by organization. This is the "dynamic landscape" made concrete — not just changing fitness functions, but the environment itself transforming through organism activity. Connects to: niche construction (Session 3), trace→actor crossing (H7), multi-rate environment (Vance's earlier contribution). TOP PRIORITY for sim06 design.
- Multiple fitness attractors at different rates — The termite mound works because multiple selection pressures (temperature, moisture, chemistry, light) operate simultaneously at different rates. One pressure stabilizes while another shifts, preventing convergence. This is the multi-rate environment idea (Vance's earlier contribution) but now grounded in a concrete physical analogy. Sim06 should have multiple interacting gradients, not a single fitness landscape.
- "No termite has ever felt a temperature" (Moltbook thread) — The post Vance references. A termite doesn't experience temperature as a scalar; it experiences the DOWNSTREAM EFFECTS of temperature gradients on its behavior (pheromone evaporation rates, mud plasticity, metabolic rate). This is downward causation (Hofstadter) + stigmergic mediation: the macro-scale environmental factor (temperature) doesn't directly act on the agent — it acts through the stigmergic medium, which the agent DOES experience. Implication for simulation design: agents should not read global state directly; they should experience only local stigmergic traces that are downstream of larger-scale dynamics.
From the 2026-07-27 code review (post-correction questions)
These arose from the construct-validity audit and the rerun. Items 52 and 53 did not exist as questions before the fixes — at 0/6 coexistence and with a broken detector there was nothing to compare.
What distinguishes the coexisting sim05 pairs from the rest? — DONE (Session 12). Corrected, sim05 gives 2/6 L2 coexistence. Analysis of the committed
results.jsonfound: (1) size symmetry is necessary but not sufficient — pair [0,3] is 10v10 and fails; (2) run 3 is lethally self-referential — its top species are fixed-point-like forms that consume other expressions in collision, and it destroys or is destroyed in all 3 pairs; (3) run 1 (20 species) is resilient — coexists with both run 0 and run 2; (4) shared species: pair [0,1] shares 8/10 species and coexists; the other 5 pairs are disjoint — structural overlap is neither necessary nor sufficient for coexistence; (5) the mechanism is collision dynamics, not structure or glue. Removing run 3, coexistence is the majority outcome (2/3). This sharpens H10: the bottleneck may be collision dynamics (dynamical compatibility), not space. Seeconcepts/sim05-coexistence-analysis.md. (Corrected 2026-07-27: an earlier version claimed "zero shared species in any pair" based on truncated top-species data — pair [0,1] actually shares 8/10.)Does NON-SATURATING negative feedback consolidate? — Two attempts at negative feedback both fragmented the structure (sim06 self-maintenance: 66–109 → 219–297 components; sim07 transport: 57 → 128 pillars). Both acted through the pheromone field, whose deposit response
p = base + gain·φ/(1+φ)is flat above φ≈1 — so both destroyed spatial contrast rather than creating it. The refined H7 prescription is negative feedback through a channel that does not saturate: a density cap, a refractory period after deposition, or directional bias along existing walls, acting on deposit probability directly rather than on the cue field. Substantially cheaper than directed transport and it tests the sharper claim. This is the direct test of [[hypotheses/H11]]; see also the 2026-07-27 refinement in [[hypotheses/H7]].Repeat sim06's parameter sweep against the working detector — DONE (Session 12). The Part-8 sweep was run against a detector that could not fire, so its conclusion — "no regime produces the crossing" — was unsupported. A broad sweep (material_decay x deposit_base x phero_follow x maintain_gain x self_maintenance, 2,100 combos, 1,225 unique runs) against the corrected detector found 1,204 combos where the crossing fires — 57% of parameter space. The crossing is not a rare edge case; it is the majority outcome. H7's crossing is reachable within the existing model. Determinism verified (two identical runs produce identical results). Results in
simulations/sim06_termite_mound/output/sweep_crossing_results.json. The crossing tends to fire with higher phero_follow (0.9-0.95) and moderate material_decay (0.001-0.004). This changes H7 from "near miss" to "crossing confirmed in the existing model."Re-derive the P_catalyze diagnosis for sim04. — The suspicion that P=0.005 is too high (producing one large core rather than distinct cores) came from a run that was not reproducible. sim04 is now deterministic and gives 3 cores in both conditions. Re-derive before sweeping — and note that in a 510-species space that exhausts, the fixed-vs-evolving comparison may be uninformative regardless of P (see the 2026-07-27 refinement in [[hypotheses/H9]]).
Should sim05's "L1 organizations" be tested for closure and self-maintenance? — sim05 reports surviving species sets and calls them L1 organizations, drawing an analogy to COT's closure + self-maintenance. It never tests either property. Either implement the test — sim03 already has a (structural, non-flux) version of it — or restate what sim05 measures. This matters because the L1/L2 framing is what connects sim05 to H10 and to the COT literature.
From Session 13 (2026-07-28)
sim09: the curvature channel — a non-saturating rule that RECRUITS as well as LIMITS — TOP PRIORITY for the next nightly session. sim08 confirmed H11's direction (a non-saturating density cap consolidates morphology — pillars 101→52 — where cue-field feedback fragmented), but the cap alone did not fire the crossing: it limits growth without recruiting maintenance, so stability didn't rise. The curvature channel (Calovi et al. 2019) is the one non-saturating channel that does BOTH: depositing at a concavity fills it (limits) AND extends the concavity nearby (recruits further building at the edge). It is also the minimal lumped form of the "directed transport" H7's Session-10 refinement called for — curvature IS directed geometry. sim09 should add a curvature/deposition-edge rule to sim06's Grassé model: loaded termites preferentially deposit at concavities (high local curvature of the material field), excavate/ avoid convexities. Prediction: this consolidates AND the crossing fires, because the channel recruits as well as limits. This is the cheapest remaining candidate that could actually cross, and it is grounded in what real termites do. See
concepts/non-saturating-channels.md. UPDATE (Session 14, 2026-07-29): Grounding complete. Facchini et al. 2020 (J R Soc Interface) built a curvature-only phase-field growth model (no pheromone field) that reproduces real nest morphology, with a phase parameterd(linear instability → walls branch/merge/invade space). Facchini et al. 2024 (eLife) unified curvature≡evaporation flux and confirmed no cement pheromone (2 independent groups). The convex/concave contradiction is resolved (different action components). The growth equation ∂f/∂t = f(1−f)·[(1/2)·Δf + d·Δ²f] gives sim09 its recruit (mean curvature Δf), limit (smoothing d·Δ²f), and surface-restriction (f(1−f)) terms. Public code: github.com/oiluigioi/JRSI_2020_termite_nest. DONE (Session 15, 2026-07-30): DESIGN.md authored atsimulations/sim09_curvature_channel/DESIGN.md— 9 independently-implementable Parts mirroring sim06's structure. The Facchini 2020 growth equation∂f/∂t ≈ f(1−f)·[(1/2)·Δf + d·Δ²f]is adapted to sim06's 2D grid+agent framework: state-gated deposit at convex tips (loaded) / excavate at concavities (unloaded) — the Facchini/Calovi action-component resolution made operational; a linear (non-saturating) deposit-probability routing on curvature; thef(1−f)surface-restriction prefactor as anon_surfacedilation mask; thed-gated biharmonic smoothing as the phase-transition knob (sim09's analog of sim07'sM_c); roughness (curvature std over surface) as the recruit proxy and the channel-adapted crossing criterion 2; baseline_pheromone condition (sim06's saturating rule) as the control. Part 7'sdsweep is the headline phase-transition plot. NEXT: implement Part 1 (GLM, next nightly). UPDATE (Session 17/18, 2026-08-02): sim09 FULLY IMPLEMENTED — all 9 Parts [x]. Part 9 (visualize.html + README.md) shipped; verification passes (selftests OK, run produces results.json, local http server 200 for page/results/sweep). At default params (d=1.0) neither condition crosses — the curvature channel grid-saturates (10000/10000 cells) because the nucleation base floods the grid before curvature routing creates spatial selectivity, so crossing criterion 2 (roughness sustained while mass saturates) cannot fire. Tuned probes (deposit_prob_base=0.01, material_decay=0.002) confirm the predicted consolidation DIRECTION (pillars 25→2 as d rises 0→4, roughness spike at the biharmonic instability) — opposite of sim06/sim07 fragmentation, H11's direction in a 4th mechanism. Perturbation: curvature 1.13× vs baseline 47.34× (baseline inflated by unbounded accumulation). The crossing is now a parameter-regime question, not a mechanism question. NEXT PRIORITY: a broaddeposit_prob_base × material_decay × dsweep in the mass-saturating regime (low nucleation, higher erosion) to located*, plus a spatially-targeted recovery metric distinguishing scar repair from volume restoration. DONE (Session 19, 2026-08-03): The d* sweep ran (100 combos,dpb × decay × d). 0/100 crossed under the original detector. Per-criterion diagnosis: criterion 2's mass-saturation gate (|growth_rate|<0.01) passed 0/100 — it was an unfalsifiable metric-ceiling bug, its threshold ~100× below the Poisson noise floor of a 150-termite deposit process (the sim06 detector-bug lesson repeating). Corrected to a relative-slope plateau (|slope(M)|/mean(M)<0.001over K=16 samples): the crossing now FIRES in the curvature channel at every d∈[0,4] in the tuned probe (non-saturating grid, cells 3123–5754/6400) and does NOT fire in the baseline-pheromone control (same detector, 0/3 — saturating rule never elevates the pheromone cue enough). crossing_step 1550→900, pillars 12→1, roughness 0.44→0.77 as d rises. Determinism verified (0/80 diffs). Honest limitation: the crossing fires at d=0, so the recruit half drives it; the limit half (d-smoothing) consolidates morphology but is not necessary for the verdict. Seedstar_sweep.pyand H7/H11 Session-19 refinements. NEXT PRIORITY: isolate the recruit and limit halves (recruit-only d=0 vs limit-only no-curvature-routing) and build a spatially-targeted recovery metric.Curvature as the minimal form of directed transport — Session 10 concluded sim07's scalar transport needed to be directed (channel geometry carrying cue to building fronts). The curvature channel may BE that minimal directed geometry: depositing at concavities routes building along edges, not away from them. sim09 would unify the "directed transport" and "non-saturating inhibition" candidates into one mechanism — falsifiable: if curvature routes AND recruits, it should fire the crossing where the scalar (sim07) and the cap (sim08) both failed.
From Session 19 (2026-08-03)
Recruit-vs-limit isolation — which half of the curvature channel drives the crossing? — TOP PRIORITY for the next nightly session. The corrected detector fires the crossing at d=0 (no biharmonic smoothing), which means the recruit half (curvature routing + mass plateau) is sufficient for the verdict and the limit half (d-smoothing) is not necessary — it only consolidates morphology (pillars 12→1, crossing_step 1550→900). H11's "recruit as well as limit" refinement (Session 13) is therefore half-supported. A clean test needs two new conditions in sim09: (a) recruit-only — curvature routing ON, d=0 (smoothing OFF); (b) limit-only — d-smoothing ON, curvature routing OFF (termites follow random walks, no curvature-biased movement, but the biharmonic still smooths the field). If recruit-only crosses and limit-only does not, the recruit half is the load-bearing variable and H11's "limit" half is a morphology optimizer, not a crossing requirement. If both cross, the mass-plateau gate is too permissive (the crossing is detecting any stable plateau, not the curvature mechanism specifically). This directly tests whether H7's Session-13 "recruits as well as limits" prescription is necessary or just sufficient. DONE (Session 20, 2026-08-04): The 2×2 factorial (recruit ON/OFF × limit ON/OFF, 4-seed robustness pass) found the recruit half is necessary and almost-sufficient for a stable crossing: recruit-only (d=0) is stable 3/4 seeds (hold 1.00 in 3, 0.65 in the borderline seed); neither (no recruit, no limit) is 0/4. The limit half alone is never stable (0/4 — criteria flicker, hold 0.40–0.55, because the biharmonic shapes convex geometry no agent is routed to; criterion 3
deposits_on_convex_fractionoscillates around 0.60). But the limit half is a stability amplifier: recruit+limit is stable 4/4 where recruit-only is 3/4 — the borderline seed becomes fully stable (hold 1.0) when d>0 is added. So "recruit as well as limit" = recruit necessary + almost-sufficient; limit = stabilizer + morphology optimizer (causal, not strictly necessary). The decisive contrast is recruit ON vs OFF at d=0 (same detector, same regime, only the recruit flag differs). A newstable_crossedmetric (late_hold_rate≥ 0.90) separates the recruit half's stable crossing from the limit half's transient flicker. Determinism verified. Seerecruit_limit_sweep.pyand H7/H11 Session-20 refinements. NEXT PRIORITY: spatially- targeted recovery metric (#60), then L2 composition (#62).Spatially-targeted recovery metric — scar repair vs volume restoration — DONE (Session 24). The grid-wide
recovery = total_material / pre_perturb_totalcannot distinguish "repair at the scar" from "continued growth elsewhere." The baseline's 47.34× "recovery" was the cleanest demonstration — it was unbounded material accumulation, not targeted repair. A spatially-targeted variant (recovery measured in the damaged patch specifically:material_in_patch / pre_perturb_material_in_patch) would make the perturbation acid test decisive. Implemented aspatch_recoveryin sim09.py, plus amirror_recoverycontrol arm (an undamaged same-size region). Result:targeted_repair = patch_recovery − mirror_recoveryis negative in all four conditions (tuned: curvature −1.95, baseline −1.65; default: curvature −0.60, baseline −51.0). Neither channel preferentially repairs the damage site. The scar grows slower than an undamaged mirror (re-nucleation lag). The crossing fires but the structure does not self-repair in the targeted sense — the crossing is a stability claim, not a scar-repair claim. The Session 17 "self-repair" report was an artifact of the grid-wide metric. Determinism verified. Seepatch_recovery_probe.py, H7 Session-24 refinement.The mass-plateau gate as a reusable methodology pattern — The sim06 and sim09 detector-bug corrections share a pattern: a threshold set below the noise floor of the quantity it gates on, making the detector unfalsifiable. sim06's deposit-rate gate could not fire because Grassé positive feedback makes deposit probability rise; sim09's mass-saturation gate could not fire because Poisson window noise sits ~100× above the threshold. Both were caught by computing the metric's ceiling. This is now earned twice and deserves to be a standing methodology rule for any future detector: before running a parameter sweep, compute the noise floor of every gated quantity and verify the threshold sits above it. Could be added to CLAUDE.md §4 step 6 as a checklist item.
Does the crossing compose? — the L2 question with a non-saturating glue — If the curvature channel crosses (it does, Session 19), do two self-maintaining curvature structures compose into a higher-level entity? This is the sim05 L2 question reopened with a non-saturating stigmergic glue — the direct test of H1/H10. sim05's 2/6 coexistence used collision dynamics as the glue; a curvature-channel glue (two structures whose curvature fields interact) might compose more reliably. Candidate sim10 or a sim09 extension: run two curvature-channel structures in adjacent grids with a shared boundary and test whether a composite organization emerges.
From Session 20 (2026-08-04)
The borderline-seed question — what makes seed 123 unstable for recruit-only? — Recruit-only (d=0) is stable in 3/4 seeds but borderline in seed 123 (hold 0.65, still crosses). Recruit+limit (d=1) is stable 4/4 — the limit half rescues seed 123. What is different about seed 123's nucleation trajectory that makes the recruit-only crossing unstable, and is it a morphological difference (initial deposit scatter) or a dynamical one (criterion 3 flickering near threshold)? If it is nucleation scatter, the limit half (smoothing) regularizes it; if it is dynamical, the limit half stabilizes criterion 3 indirectly. Inspect seed 123's history for recruit-only vs recruit+limit: where does hold drop (which criterion flickers), and does d-smoothing fix that criterion specifically? Cheap analysis of the committed sweep JSON; no new runs needed.
A saturating-action control — disentangling "action-based" from "non-saturating" — H11's evidence (sim08 cap, sim09 curvature) is both action-based AND non-saturating simultaneously, so it cannot fully distinguish "action-based" from "non-saturating" as the causal variable (this was flagged in H11's Criticisms section from Session 13). The recruit-vs-limit isolation sharpens this: the recruit half is action-based (curvature routes deposit/excavate selection) and non-saturating (linear gain). A saturating-action control — a deposit-probability routing that saturates (
p = base + gain·curvature/(1+|curvature|)) rather than the linearp = base + gain·curvature— would isolate the two factors. If a saturating recruit half still crosses stably, "action-based" is the load-bearing property; if it degrades to a transient flicker (like limit-only), "non-saturating" is. This is the clean test of H11's core distinction, currently confounded. DONE (Session 21, 2026-08-05): The 2×2×2 factorial (response {linear, saturating} × recruit {ON, OFF} × d {0, 1}, 4-seed robustness pass) found action-based is the primary load-bearing property; non-saturating is a secondary stability amplifier. The saturating action crosses in 8/8 recruit-ON seeds (stable 6/8); the linear action crosses in 8/8 (stable 7/8). The limit half rescues both to 4/4 at d=1. Saturation costs ~0.05 in mean hold rate at d=0 (0.91→0.86) but does not collapse the crossing — criterion 3 (deposits_on_convex_fraction) holds 1.00 for both forms; only the mass-plateau gate (criterion 2p) flickers more under saturation. H11's strict "non-saturating" claim is partially weakened: a saturating action-based channel still crosses stably, but less robustly. The "self-defeating" language applies to cue-based saturating channels (sim06/sim07), not to action-based saturating channels. The three-level causal decomposition: (1) action-based routing = primary, (2) non-saturating response = secondary stability, (3) biharmonic smoothing = tertiary stability + morphology. Seesaturating_action_sweep.pyand H7/H11 Session-21 refinements. NEXT PRIORITY: spatially-targeted recovery metric (#60), then L2 composition (#62).The stable_crossed metric as a reusable methodology pattern — The cumulative
crossedflag (set once criteria hold forCROSSING_PERSISTconsecutive samples, never unset) hides the difference between a crossing that holds and one that flickers on and off. Thelate_hold_rate(fraction of late-window records where all criteria hold) exposes it. This is now earned once (sim09 Session 20: the limit half's transient crossing was invisible until late_hold_rate was computed) and deserves to be a standing metric for any crossing detector: report both the cumulative verdict AND the late-window hold rate. A crossing that fires then degrades is not the same phenomenon as one that holds. Could be added to CLAUDE.md §4 step 6 alongside the metric-ceiling rule (#61).Continuous Game of Life — self-organizing cells at the edge of growth (Guillet & Jülicher 2026) — A continuous-space, continuous-time Game of Life (cGoL) that produces self-replicating, motile, dying cell-like patterns with just 7 parameters. The key finding: a global resource constraint (conservation law) causes the system to self-organize to a phase transition boundary — the "edge of growth" — where morphologies are richest and most life-like. Reference code cloned to
simulations/cGoL_reference/(Julia, FFT-based convolution, GPLv3). Paper: arXiv:2607.27402, to appear in Artificial Life journal.Relevance to our hypotheses:
- H1/H7 (Composition / Trace→Actor Crossing): The cGoL cell patterns have a nucleus+shell structure that emerges from simple convolution rules — a spatially organized, self-maintaining entity. The field L is the "trace"; the emergent cell with homeostatic morphogen concentrations is the "actor." Self-replication and persistence of these cells is a concrete trace→actor crossing. Can we layer H7's crossing detector onto the cGoL cells? Do they satisfy the three operational criteria (persistence, non-reducible dynamics, constraint on agents)?
- H4 (Dynamic Environment): Resource feedback is exactly H4 — the environment participates in a feedback loop. Growth consumes resource → resource depletion retunes parameters → system self-organizes at the phase boundary. A stigmergic medium with its own dynamics.
- H11 (Saturating Channel): The "edge of growth" is a non-saturating channel — resource scarcity acts as feedback that doesn't saturate the way a pheromone field does. The system self-tunes to the transition boundary rather than collapsing.
- H8 (Computational Irreducibility): The phase structure is mapped empirically through extensive simulation — morphologies at the edge of growth can't be predicted from rules alone.
- Multi-scale composition (H1/H10): The cell-like patterns interact, divide, and collide. Whether two such self-maintaining patterns compose into a higher-order structure is directly testable.
The reaction-diffusion interpretation (§4) maps the cGoL onto morphogen concentrations held at homeostatic levels by the nonlinear survival rule — connecting to sim03 (chemical organizations) and sim09 (curvature channel). The "survival rule" is a non-saturating channel that maintains homeostasis.
Next step: port
cGoL_minimal.jlto Python (numpy FFT convolution, ~100 lines), add resource feedback, and test whether the emergent cells satisfy H7's crossing criteria. The minimal Julia implementation uses: (1) two Gaussian FFT convolutions for M and N fields, (2) a sigmoid-based survival rulerule0(M,N,p), (3) explicit Euler time integration. Parameters: p=(0.50, 0.10, 0.23, 0.015, 0.35, 0.26), λ=3.
From Session 21 (2026-08-05)
A truly cue-based saturating action control — completing the 2×2 — The Session 21 saturating-action control tested within the action-based family (linear vs saturating action routing). The remaining cell of the 2×2 is a cue-based non-saturating channel: deposit probability routed on a non-saturating cue field (e.g.
p = base + gain·φwithout saturation, instead ofp = base + gain·φ/(1+φ)). If a non-saturating cue channel crosses, then the action/cue distinction (H11's original framing) is the real divide, not the saturating/non-saturating one. If it does not, the action-based property is confirmed as primary even when the cue is non-saturating. This completes the 2×2: action×{linear,saturating} × cue×{linear,saturating}, isolating which of the two properties (action-based, non-saturating) is truly load-bearing. Cheap: the cue-based condition is sim06 with the deposit rule changed fromφ/(1+φ)to linearφ. DONE (Session 22, 2026-08-06): The cue-based non-saturating control (sim06 withdeposit_responseparameter,cue_response_sweep.py) found the non-saturating cue crosses LESS, not more — the opposite of the action family and opposite to H11's strict prediction. Without self-maintenance: saturating cue 16/16 stable (hold 1.000); linear cue 0/16 stable (hold 0.053). With SM: both 16/16 stable. Seed robustness (4 seeds) confirms. The non-saturating property reverses sign across families: it amplifies stability in the action family (sim09: 7/8 vs 6/8) but destroys it in the cue family (sim06: 0/16 vs 16/16 w/o SM). Mechanism: the linear cuep = base + gain·φclamps to p=1.0 at φ≈1.15, flattening the gradient (mean pheromone drops to 0.467 < 0.5 threshold); the saturating cue'sφ/(1+φ)compression prevents deposit-probability saturation and preserves spatial contrast. The "self-defeating" channel is the non-saturating cue (deposit-probability clamping), not the saturating cue — H11's original framing was backwards for the cue family. Self-maintenance rescues the linear cue (4/4 stable). Seecue_response_sweep.py, H7/H11 Session-22 refinements. NEXT PRIORITY: spatially-targeted recovery metric (#60), then L2 composition (#62).The three-level causal decomposition as a methodology pattern — Sessions 19–21 decomposed the crossing's causal structure into three levels: (1) action-based routing (primary — the causal variable separating crossing from non-crossing), (2) non-saturating response (secondary — stability amplifier), (3) biharmonic smoothing (tertiary — stability amplifier + morphology optimizer). This is a generalizable pattern: when a hypothesis claims two properties matter (H11: action-based AND non-saturating), a single confounded experiment cannot distinguish them; a factorial isolating each property separately, plus a seed-robustness pass with a stable-vs-transient metric, can. The pattern: (a) identify the confounded properties, (b) build a saturating control that holds one constant, (c) run a 2×2×2 factorial, (d) use late_hold_rate to separate stable from transient effects, (e) decompose the result into primary/secondary/ tertiary causal levels. Could be added to CLAUDE.md §4 step 6 alongside the metric-ceiling and stable_crossed rules.
The borderline-seed flip — seed 123 is borderline for linear but stable for saturating — Session 20 found seed 123 is the borderline seed for linear recruit-only (hold 0.65). Session 21 found seed 123 is stable for saturating recruit-only (hold 0.95) — and seeds 42 and 256 are borderline for saturating (holds 0.85, 0.70) but stable for linear. The borderline seeds flip between response curves. This means the linear and saturating forms are not simply "one more stable than the other" — they are fragile to different nucleation trajectories. What makes a seed borderline for one form but not the other? If the nucleation scatter differs, the saturating form's compressed gain may regularize seeds where linear's high gain overshoots, while linear's full gain may stabilize seeds where saturating's compression is too weak. Inspect the borderline seeds' histories: does the hold drop at the same criterion, and does the response curve change which criterion flickers? Cheap analysis of the committed sweep JSON; no new runs needed.
From Session 22 (2026-08-06)
Deposit-probability clamping vs cue-response compression — the two kinds of "saturation" — Session 22's cue-based control revealed that H11's original framing conflated two distinct saturation phenomena: (a) cue-response compression (the
φ/(1+φ)form flattens at high φ — what the saturating cue has) and (b) deposit-probability clamping (the lineargain·φform hits p=1.0 at φ≈1.15, so every high-pheromone cell deposits at 100% — what the non-saturating cue has). The "self-defeating" saturation is (b), not (a): the non-saturating cue clamps to p=1.0 and flattens the gradient; the saturating cue's compression prevents clamping and preserves spatial contrast. This distinction should be formalized: a channel is self-defeating when its response curve saturates the probability (the output), not when it compresses the cue (the input). H11's "self-defeating saturating channel" should be re-read as "self-defeating probability-saturating channel." This is a refinement of the concept, not a new experiment — but it deserves a formal write-up and possibly a concept file, because it changes how the 2×2 should be interpreted. The action family's response curve (base + gain·cvsbase + gain·c/(1+|c|)) saturates only the gain (the routing decision is preserved); the cue family's response curve saturates the probability (the output clamps). That is why the sign reverses.The self-maintenance rescue — is SM necessary or merely sufficient for the non-saturating cue? — Session 22 found self-maintenance rescues the non-saturating cue completely (0/16 → 16/16 stable). But is SM the only mechanism that can rescue it, or would any pheromone-sustaining mechanism work (e.g. slower pheromone decay, higher deposit pheromone, lower diffusion)? If the linear cue's failure is purely "mean pheromone drops below 0.5," then any mechanism that keeps pheromone elevated should rescue it — and SM is just one way to do that. A sweep of pheromone_decay × deposit_pheromone at the linear-cue condition would map the rescue surface. If the rescue is specific to SM (the structure-reemits-pheromone loop), that connects to H7's self-maintenance crossing mechanism; if it is generic (any pheromone elevation), the non-saturating cue's failure is just a parameter-regime issue, not a mechanistic one. Cheap: a small sweep around the linear-cue condition.
The deposit-probability saturation threshold as a predictor — DONE (Session 23). The φ_sat predictor (the input value at which p_deposit first reaches 1.0) was tested as a unifying diagnostic across all four cells of the 2×2. A direct probe (
phi_sat_probe.py) of sim06 (cue) and sim09 (action) at their crossing-proven regimes found the predictor is 50% accurate — no better than chance. It correctly predicts the cue family (saturated→fails, unsaturated→crosses) but fails for the action family: action/linear is saturated (max curvature 2.55 > c_sat 1.165, clamp fraction 1.0%) but still crosses stably. The clamping fraction is tiny everywhere (0–7%). The difference: in the cue family, the deposit probability IS the spatial signal — clamping it destroys the gradient. In the action family, spatial contrast lives in the routing decision (which direction the agent moves), not the deposit probability — the response curve saturates the gain (how hard to deposit), not the routing (where to go). The unifying diagnostic is whether spatial contrast in the routing input survives the response curve, which depends on channel architecture, not just the saturation threshold. Determinism verified. Seephi_sat_probe.pyand H7/H11 Session-23 refinements.
From Session 23 (2026-08-07)
- The two-wire principle — feedback signal and spatial signal on separate channels — Session 23's φ_sat probe found the predictor fails for the action family because spatial contrast survives via the routing decision (which direction to move), not the deposit probability. The cue family puts the feedback signal and the spatial signal on the same wire (the pheromone field → deposit probability → spatial contrast); saturating one destroys the other. The action family puts them on separate wires (curvature → routing decision for spatial contrast; curvature → deposit gain for feedback); saturating one leaves the other intact. This is a generalizable design principle: a self-defeating channel is one where the feedback signal and the spatial signal travel on the same wire. Does this principle hold beyond stigmergic channels? In ACO, the pheromone trail IS both the feedback signal and the spatial signal — but ACO's response function (τ^a·η^β) is unbounded, so it never saturates. In development, morphogen gradients carry positional information (spatial signal) AND feedback (concentration-dependent gene expression) on the same wire — and morphogen saturation is a known developmental pathology. This deserves a concept file and possibly a cross-domain synthesis. Cheap: no new runs; pure synthesis.
From Vance (2026-08-04)
Singh et al. (2025/2026) — MARL-trained weakly electric fish collectives: emergent social behavior from biophysical sensing + individual fitness reward — arXiv:2511.08436. Found via a Bluesky follower (Naomi Saphra is a co-author). The paper is a complete worked example of several things our project has been circling, and it connects to at least five of our hypotheses:
H1 (Multi-scale composition) ↔ emergent collective behavior from individual incentives alone. The paper's central claim: collective foraging, dominance hierarchies, aggression, and context-dependent EOD communication all emerged from individual fitness rewards with no reward for communication, coordination, chasing, or aggression. This is the same "emergence from individual incentives" pattern our project studies, but at a single scale (fish-to-fish). The open question for us: does their framework compose across scales? Their fish are homogeneous agents with the same policy — can heterogeneous policies at different scales produce multi-scale composition?
H7 (Trace→Actor Crossing) ↔ EOD as a stigmergic medium. The EOD is a stigmergic signal: it modifies the electric field (environment), persists briefly, is sensed by conspecifics, and influences their behavior. The paper's Mormyromast "cons-image" (detecting conspecific EOD distortions) IS stigmergic sensing — agents read the environmental trace of another agent's action. The Knollenorgan (long-range conspecific-only sensor) is a dedicated stigmergic channel. The paper shows that ablating the Knollenorgan doesn't affect foraging but reshapes social organization (more aggression, less spacing) — the stigmergic medium is causally efficacious for social structure, not just foraging. This is direct evidence for H4 (dynamic environment as participant, not backdrop) and the ANT claim that the medium is an actant.
H11 (Saturating Channel) ↔ EOD self-cancellation. The Mormyromast has an internal cancellation signal that suppresses the reafferent (self-generated) EOD component — the self-field is ~729× stronger than the conspecific field at 10 cm, so without active cancellation the self-signal would saturate the sensor and mask the conspecific signal entirely. This is a biological instance of our "two-wire principle" (queued topic 73): the self-image and the cons-image travel on separate wires (separate processing channels within the same receptor), so saturation of the self-signal doesn't destroy the cons-specific spatial signal. The paper's "collective sensing" experiment (gating self- vs cons-EOD inputs independently) is exactly the kind of channel-factor decomposition our Session 21–22 factorial experiments did with sim09.
H4 (Dynamic Environment) ↔ the electric field as a shared stigmergic medium. The electric field is not a static backdrop — it's co-determined by all agents' EODs AND the environment (walls, prey distort it). Agents sense not only each other but "how their own and others' EODs are transformed by the shared environment." This is niche construction in the electric domain: agents modify the field they sense through, and the field's distortions carry information about the environment. The field IS the stigmergic medium.
RNN dynamics ↔ cross-scale neural representation. The effective dimensionality of RNN activity scales with group size only when the Knollenorgan (long-range stigmergic channel) is intact — ablate it and dimensionality stays flat at the solo baseline. The social context expands the neural representation space, and this expansion is driven by the stigmergic channel, not by direct interaction. Proximity-dependent correlated latent dynamics (PLSC) between interacting agents' RNN states collapse to zero beyond communication range. This is a potential model for how multi-scale composition could work in a neural system: the stigmergic medium creates a shared subspace between agents that doesn't exist at the individual level — a new dynamical degree of freedom. Could our sim09 curvature structures show a similar dimensionality expansion when two structures interact through a shared curvature field?
Methodological relevance: their in silico intervention design (ablate sensors, silence EODs, change food distribution) is exactly the kind of causal decomposition our project uses. Their "seed selection criterion" (balance biological desiderata across multiple converged policies) is a pattern we could adopt for our sim runs. Their GRU-based actor-critic with recurrent dynamics analysis (PCA, linear decoding, PLSC, power spectrum) is a toolkit we haven't used but could apply to sim09's agent states.
NEXT: This should be a nightly research session topic. The paper deserves a full concept file and a synthesis entry. Key questions: (1) Does the EOD stigmergic medium satisfy our H7 crossing criteria? (2) Can their MARL framework be extended to multi-scale composition (heterogeneous policies at different scales)? (3) Does the two-wire principle (self/cons-image separation in Mormyromasts) generalize to our action/cue channel distinction? (4) Can RNN dimensionality analysis detect when a stigmergic medium creates a new dynamical degree of freedom?
From Session 24 (2026-08-08)
The control-arm methodology pattern — a metric that responds is a description, not a test — Session 24's spatially-targeted recovery metric revealed that every metric in this project needed a control arm to become a test rather than a description. The crossing detector responded to stability (needed the baseline-pheromone control); the recovery metric responded to growth (needed the mirror patch); the φ_sat predictor responded to saturation (needed the action/linear condition). A metric that responds to a phenomenon but cannot distinguish it from confounds is a description, not a test. This is now earned three times (mass-saturation gate, φ_sat predictor, grid-wide recovery) and deserves to be a standing methodology rule: before claiming a metric tests a phenomenon, identify the confound and add a control arm that holds it constant. Could be added to CLAUDE.md §4 step 6 alongside the metric-ceiling rule (#61) and the stable_crossed rule (#65).
Late perturbation after true mass plateau — The current perturbation hits at 60% of steps, when the crossing has fired but the total material is still rising (not truly plateaued). A later perturbation (80-90% of steps, after the mass has equilibrated) may give a different self-repair result: the structure would be at equilibrium, and scar repair would be purely about restoring the damage, not about continuing growth. If the late-perturbation targeted_repair is still negative, the no-self-repair finding is robust; if it becomes positive, the current result is a timing artifact. Cheap: change perturb_at and re-run the probe.
The L2 composition question with a non-saturating glue — Now the top priority. The curvature channel crosses (Session 19), but it does not self-repair (Session 24). Does it compose? Run two curvature-channel structures in adjacent grids with a shared boundary and test whether a composite organization emerges. This is the sim05 L2 question reopened with a non-saturating stigmergic glue — the direct test of H1/H10. (Previously queued as #62; now the next major test after #60 is done.)