Regenerate the retracted figures, and the v1.1 material that followed

.github/workflows/model-contracts.yml

@@ -187,6 +187,54 @@ jobs:
                  " -- a check was added or silently lost; update EXPECTED_CHECKS")
           cat("R verifier ran complete.\n")
 
+      # PP-SORP-01's producer.  The ledger quotes K = 29.72 as the concentration
+      # ratio at which the shipped Henry form and Langmuir diverge by 10%; a
+      # number nobody can re-run is the defect PP-62-11 records and the P1 Codex
+      # raised on #23.  Run here because this job already has R and deSolve.
+      - name: Reproduce the Henry/Langmuir divergence bound
+        run: Rscript analysis/henry_langmuir_bound.R --report hl-verify.json
+
+      # Same receipt discipline as above, and for the same reason: the producer
+      # reports SKIPPED rather than failing when deSolve is missing, so a job
+      # whose setup step was dropped would exit 0 having measured nothing.
+      - name: Assert the divergence-bound producer ran everything
+        shell: Rscript {0}
+        run: |
+          EXPECTED_CHECKS <- 3L
+          if (!file.exists("hl-verify.json"))
+            stop("no receipt: the divergence-bound producer did not complete")
+          r <- jsonlite::fromJSON("hl-verify.json", simplifyVector = TRUE)
+          for (k in c("checks_run", "failures", "skipped", "complete"))
+            if (is.null(r[[k]])) stop("receipt is missing '", k, "'")
+          if (!isTRUE(r$complete))
+            stop("receipt does not declare itself complete: complete=",
+                 format(r$complete))
+          int1 <- function(k) {
+            v <- r[[k]]
+            if (!is.numeric(v) || length(v) != 1L || is.na(v) ||
+                v < 0 || v != trunc(v))
+              stop("receipt field '", k, "' is not a non-negative integer ",
+                   "scalar: ", paste(format(v), collapse = ", "))
+            as.integer(v)
+          }
+          ran <- int1("checks_run"); fail <- int1("failures"); skip <- int1("skipped")
+          if (skip > 0L)
+            stop("producer skipped ", skip, " check(s): ",
+                 paste(unlist(r$skips), collapse = "; "))
+          if (fail > 0L) stop("producer reported ", fail, " failure(s)")
+          if (ran != EXPECTED_CHECKS)
+            stop("expected ", EXPECTED_CHECKS, " checks, receipt says ", ran,
+                 " -- update EXPECTED_CHECKS deliberately")
+          # AND THE QUOTED NUMBER ITSELF, because "3 checks passed" does not say
+          # WHICH bound they passed about.  The ledger cites 29.72; if the
+          # producer starts emitting a different K the row is stale and this is
+          # the only place that would notice.
+          if (!is.numeric(r$K_transient) || abs(r$K_transient - 29.72) > 0.01)
+            stop("K_transient drifted from the value PP-SORP-01 quotes (29.72): ",
+                 format(r$K_transient))
+          cat(sprintf("divergence-bound producer complete; K_transient = %.2f\n",
+                      r$K_transient))
+
   manuscript-build:
     runs-on: ubuntu-latest
     timeout-minutes: 25

README.md

@@ -82,9 +82,10 @@ than a simulation output. At the shipped constants `I0 = 1.0` and `T_cpm = 5.0`:
 | `ΔH_mel` at the reported M = 1.44 | −0.720 | **1.155** |
 
 `β_ion` — the one parameter Table 2 tabulates per species, and the one the entire sign convention is
-written around — biases acceptance by **one part in 10⁵** for exactly the species whose radial
-stratification is the headline result. The melanin term biases it by **15.5%**: four orders of
-magnitude more, through a coefficient of `0.5` hard-coded at its call site in `compute_delta_H`,
+written around — biases acceptance by **one part in 10⁵** *for one role of the two negatively
+signed species occupying a site*, which is not the term's reach: signed by role it reaches
+7.505e-2. The melanin term biases acceptance by **15.5%**, larger by about an order of
+magnitude (9.6 in ΔH), through a coefficient of `0.5` hard-coded at its call site in `compute_delta_H`,
 appearing in no table and in no configuration file.
 
 Radiation still reaches the dynamics. It reaches them **indirectly**, because `melanin_drive` is
@@ -400,7 +401,8 @@ the **core**, not the wall; `β_ion` is negative for the two radiotropic fungi a
 *B. subtilis*.
 
 Spatial sorting is dominated by the hand-specified adhesion matrix `J`. Within the radiation-derived
-pathway, melanin dominates `β_ion` by four orders of magnitude (see the melanin finding above).
+pathway, melanin exceeds `β_ion` by about an order of magnitude — 9.6 in ΔH — not by four (see
+the melanin finding above).
 `tests/deterministic_radiation.jl` demonstrates the inertness directly: it must *deliberately amplify*
 $\beta$ to make multi-seed drift detectable, because the production values are not. The local ΔH sign
 per species is a **contract**; the collective drift is a **result**; the test keeps them separate.
@@ -710,7 +712,7 @@ centimetres or grams, because no conversion to any of those exists.
 | Result | Value | Status |
 |---|---|---|
 | *C. sphaerospermum* melanin, MCS 100 | 1.44 | **Dimensionless field units.** Mean of a dimensionless field over occupied sites. The ordering among the three producers is set by the input `α_M`; the magnitude and spatial distribution are emergent. |
-| Melanin acceptance bias at that value | 1.155 (15.5%) | **Dimensionless.** Computed from the shipped `0.5` coefficient and `T_cpm = 5.0`; four orders of magnitude above the `β_ion` bias of 1.000010. The dominant radiation-derived term, and undeclared until 2026-08-15. |
+| Melanin acceptance bias at that value | 1.155 (15.5%) | **Dimensionless.** Computed from the shipped `0.5` coefficient and `T_cpm = 5.0`; about an order of magnitude above the `β_ion` bias AT THE TERM'S REACH, 1.0151 (a 10.2x ratio of excesses). Not above 1.000010, which is one ROLE of one pair of species occupying a site and gives 15,500x -- comparing against it is the withdrawn comparison. The dominant radiation-derived term, and undeclared until 2026-08-15. |
 | Membrane integrity after 100 MCS | m = 0.779 | **Dimensionless.** No physical dose exists: `Ḋ = 1.0` is a placeholder, `dt_rd = 0.5` is an uncalibrated MCS→second knob, `seconds_per_mcs` is `NaN`, and `accrue_dose!` raises rather than run without it. There is no value of Gy this corresponds to. |
 | Permeability ratio | P_eff / P₀ = e ≈ 2.72 | **Dimensionless, and exact by construction.** A closed-form function of four hand-set constants ($\alpha_P \dot{D} \Delta t_{rd} n_{MCS} = 1$); the simulation computes no permeability in cm s⁻¹. The former "0.010 → 0.027 cm s⁻¹" was the Nafion-117 literature prior `P₀` rescaled by that factor, and a literature prior is not a calibration. |
 | Interior contaminant, unweighted node mean | c / c_ext = 0.024 | **Dimensionless ratio.** Measures **non-penetration, not depletion**: `c(t=0) ≡ 0` everywhere, so the interior never held contaminant to remove. The mean is unweighted over radial nodes, not volume-weighted. The embedded figure labels this "% depleted"; the figure is wrong and the model supports only non-penetration. |